天文観測的手法における小惑星4ベスタの研究／長谷川361

天文観測的手法における小惑星4ベスタの研究

長谷川 直

1

2014年7月24日受領，2014年8月22日受理．

（要旨）小惑星 4 ベスタは様々な天文観測的手法に用いて観測し尽くされている天体であり，Dawn 探査機がランデブーした現在，最も物理情報が取得されている小惑星であると言っても過言ではない．筆者はこれまでベスタに関わる天文観測的手法を用いた研究を行い，1）3 ミクロン帯の分光観測から，ベスタ表層に含水鉱物が存在していることを発見，2）内側小惑星帯の V 型小惑星のライトカーブ観測から，それらが破片としてベスタから飛び出した年代を示し，3）位相関数の観測から，ベスタ表層の密度を示し，更に衝効果の原因を明らかにした．

1. 過去におけるベスタの観測的研究の歴史

1945 年以前に，ベスタは精度はともかくとして小惑星の基本物理情報である直径・アルベド・自転周期・スペクトルが測定されていたことになる．

1807 年（江戸時代文化 4 年）3 月 29 日に H. W. Olbers

1951 年には，検出器に光電管を用いた非常に精度に発見された小惑星 4 ベスタは 1801 〜 1807 年の 5 年間の高いライトカーブ観測が行われ，自転周期が高精度の間に発見された 4 大小惑星（1 ケレス・2 パラス・3に計測された [5]．1967 年には，高精度のライトカー

ジュノー・4 ベスタ）の最後に発見された小惑星である．

ブデータを用いて，ベスタの自転軸のあり得る２つのベスタは地球からみて，一番明るくなる小惑星であり，解が示され（小惑星帯の小惑星のライトカーブデータ

条件が良ければ，肉眼で見ることも可能である（筆者のみから自転軸を求める場合は原理的に１つの解に絞は肉眼で見たことがある）．一番明るく見えるが故に，ることが難しい），また，位相角の変化に対する明るベスタはこれまで考えられる様々な観測手法で観測がさの変化を示す位相関数も示された [6]．Gehrels の示行われている．

した自転軸の解の１つは最終的に Dawn 探査機によっ望遠鏡にマイクロメーターを取り付けて，その視直て得られた自転軸の結果とほぼ一致している．1968

径を測るという方法を用いて，4 大小惑星も発見から年には，ベスタ近傍に近接した他の小惑星の軌道のふ直ぐに天体の大きさを測る物理観測が試みられている．

らつきから，その質量が求められた [7]．1970 年には，発見から 18 年後の 1825 年（江戸時代文政 8 年）には，

中間赤外線で天体の熱輻射量から直径を推定する方法

ベスタの直径が測定されている [1]．1883 年（明治 16から，その直径が高精度に求められた [8]．また，同年）には，ベスタ表層のアルベドと明るさの変化が観

年にベスタの可視光域の反射スペクトルが高精度に観測されている [2]．その後，1929 年（昭和 4 年）にベスタ

測され，その表面が分化隕石と似ていることが示されの可視光域のスペクトル（太陽光スペクトルで割られた [9]．1973 年には，自身で観測したデータに加えて，た反射スペクトルの形で掲載されていないが）と自転それまで得られた高精度のライトカーブデータを用い周期が測定されている [3]．1934 年には，位相角の変て，ベスタの南半球に巨大なクレーターがあることが化によるベスタの偏向度の変化も測定されている [4]．

示唆された [10]．Taylor のこの研究は後に Hubble 宇

1. 宇宙航空研究開発機構宇宙科学研究所
hasehase@isas.jaxa.jp■2014遊星人Vol23-4_3校.indd

361

宙望遠鏡の観測や Dawn 探査機の探査で，その正しさを証明されることになる．1975 年には，近赤外波長

2014/12/13

14:04:53

362

日本惑星科学会誌 Vol. 23, No. 4, 2014

域の分光観測が行われ，可視光と合わせた研究で，そな物質に覆われていて，レゴリス層は地質学的に未成の表面が分化隕石のグループの１つである HED 隕石

熟であることが示された [27]．また同年に，補償光学（ホワルダイト・ユークライト・ダイオジェナイト）とを用いて，ベスタの撮像が行われている [28]．1997 年

似ていることが示された [11, 12]．Larson 等と Johnsonには，Hubble 宇宙望遠鏡によりベスタの観測が行われ，等のこの研究によって，HED 隕石の隕石学的な研究

大きさ・形状・自転軸や北半球のアルベド分布が直接とベスタの天文観測的な研究が結びついたことになる．

的に示された [29]．南半球にある巨大クレーターが直1977 年には，スペックル干渉法により，ベスタを点接的に確認されたことより，ベスタ上の巨大クレータ

源でなく，面光源として，捉えることに成功した [13]． ー→クレーターから排出された V 型小惑星群→地球
また，同年にマイクロ波でベスタの観測に成功した

近傍まで来ている V 型近地球型小惑星→地球に落下[14]．1979 年には，自転位相によるカラー差が存在ししてきている HED 隕石と、ベスタと HED 隕石を結ぶていることが示され [15]，自転位相による偏向度の変線が繋がったことになる．

化も示された [16]．Blanco 等と Degewij 等の研究から1998 年には，ISO 赤外線宇宙天文台で取得された中ベスタ表層の異方性が示され，自転周期が決定された

間・遠赤外線域測光データと過去得られた測光データ（通常の小惑星は表層のカラー差の影響より形状の効を組み合わせて，ベスタ表層の熱慣性と中間・遠赤外果の方が大きい為に、ライトカーブ形状は１周期ダブ

線域の輻射率が示された [30]．2010 年には，1997 年のルピークであるが，ベスタは表層のカラー差の影響のHubble 宇宙望遠鏡の観測で得られていなかったベス方が形状の効果より大きかった為にシングルピークでタの南半球のアルベドマップが得られた [31]．2011 年

あった）．また，同年にレーダーによる観測も成功しには，Swift γ線バースト観測衛星や IUE 紫外線天文ている（但し，この時点ではレーダーによる３次元形

衛星，Hubble 宇宙望遠鏡の紫外光域のスペクトルから，状は得られていない）[17]．また，同年に小惑星の軌道ベスタの表層が宇宙風化作用を受けていないことを示

データを用いて，個数は少ないが，ベスタ族の存在がした [32]．2013 年には，過去の熱輻射のデータから，初めて示された [18]．以上，1950 年以降の電子技術のベスタの表層レゴリスの平均粒径が求められた [33]．

発展により，ベスタの様々な物理的情報が高精度で得そして，2011 年には，Dawn 探査機がベスタに到着・られた．

ランデブーし，その探査を基にしてベスタの形状・質

1980 年には，IUE 紫外線天文衛星でベスタの紫外量・密度・表面形態・内部構造・表層物質組成・表層

光域のスペクトルが取得された [19]．1986 年には，年代等々を調べられ [34]，現在も詳細な解析が進んでIRAS 赤外線天文衛星の全天サーベイにベスタが検出いる．

された [20]．1988 年には，スペックル干渉法によって，ベスタの観測は，地球からみて一番明るく，小惑星

ベスタの 3 軸回転楕円体の大きさが示された [21]．

帯で 3 番目に大きい小惑星（準惑星セレスを小惑星と1989 年には，マイクロ波観測からベスタ表層がレゴしてみた時に）故に，その時代での最新鋭の観測技術・

リスに覆われていることが示された [22]．1990 年には， 手法を用いて，観測が行われており，その結果，上記
小惑星のその当時最新の軌道データを用いて，ベスタの様に様々な研究成果がでている．ベスタは小惑星帯族のメンバーが本格的に示された [23]．1991 年には，の小惑星で一番物理的性質が判明していると言っても近地球型小惑星でベスタと似たスペクトルを持つ V

過言ではない．

型小惑星が見つかった [24]．なお，Cruikshank 等は観

本論文では，このようなベスタを巡る天文観測的研

測した V 型の近地球型小惑星をベスタ起源で無いと究の激しい競争の中で，筆者が進めてきたベスタ及び論文では述べているが，後述の様々な発見により，現それに関わる小惑星の研究についてまとめる．

在ではベスタ起源であると考えられている．また同年，
ベスタの掩蔽現象が観測されている [25]．1993 年には，
ベスタ族メンバー，及び，ベスタ近傍の内側小惑星帯

2. 最もよく調べられた小惑星である
ことを活かした研究に V 型小惑星があることが発見された [26]．1996 年に天文観測では天体の明るさの絶対値をどのように決は，レーダー観測によって，ベスタ表層が玄武岩質的定しているかと言うと，物差しで長さを測ることと同■2014遊星人Vol23-4_3校.indd362

2014/12/13

14:04:53

天文観測的手法における小惑星4ベスタの研究／長谷川

様に，標準星と呼ばれる明るさが分かっている星の明

363


るさと比較することによって，未知の天体の明るさを決定している．この標準星は可視光・近赤外域では数

ると標準星の数は少なくなる．理由の１つとして，可視光から波長の長いところまで，明るさを正確に予測する為には，恒星のモデルが存在しなければならない 

多く存在している [e.g., 35] が，中間・遠赤外線域になが，特定のスペクトル型の恒星しかモデルが存在していないことが挙げられる．また，これらの星は中間・

遠赤外線域になると明るさが暗くなってくるというのも理由の１つである．更に，もう１つの理由としては，「ベガ型星」のようにデブリディスクが付随されている天体は中間・遠赤外線域に明るさの超過が見られることがある．この超過は過去の観測からどの天体に超過の有無があるか天体毎にわかっているが，その超過分をモデルで予測することが事実上不可能の為に，明るさを正確に予測できない．よって，中間・遠赤外線












   

図1： 遠赤外線標準星天体のフラックス．フラックスの大きい順に，青線：天王星・海王星．赤線：1 ケレス・4 ベスタ・7
イリス・511 ダビダ・47 アグライヤ．緑線：α Booアーク
トゥルス・β Andアンドロメダ座ベータ星・α Cmaシリウス・やぎ座オメガ星・Θ Umiこぐま座シータ星

域の明るい標準星（例えば 100 ミクロンで 1Jy を超える

表1： ベスタの中間赤外線域での観測値．モデル値はStandard
Thermal Modelによるもの．

天体）は 10 天体程しか存在していない．

band
8.7 um18.8 um

一方で，電波観測では，火星・天王星・海王星といった惑星が，標準星として，使用されている．これら

observed flux
121 Jy468 JyModel
146 Jy420 Jy惑星の標準星は標準星としてはとても明るく，例えば，
一番暗い海王星でも，100 ミクロンで 300 Jy 程の明る前での小惑星標準星整備の為に，あかり赤外線天文衛

さがある．即ち，遠赤外線域では，標準星として使用

星打ち上げ以前の 2003 年に国立天文台ハワイ観測所

できる惑星と恒星の明るさに 1.5 桁ほどのギャップがの 8.2m すばる望遠鏡の COMICS 装置で中間赤外線領出来ていることになる（図 1）．

域のベスタの測光観測を行った
（表 1）
[37]．

この明るさのギャップを埋める為に中間・遠赤外線

8.7 ミクロンと 11.8 ミクロンの観測の結果，それぞ域では小惑星を標準星として使用することが提案され，れの観測と簡易熱モデル予測値の差は 10-20 ％精度あ実際に ISO 赤外線宇宙天文台では遠赤外線域の標準星

り，過去行われた研究と比較すると多少大きめな数値として使用された [e.g., Müller & Lagerros 1998]．中になった．しかしながら，地上の中間赤外線の絶対値

間・遠赤外線域の小惑星は太陽光の反射でなく，太陽の測光精度は良くても 10 ％ほどであり，簡易熱モデ

光を吸収してその熱輻射で輝いている．熱輻射によるルをしたことによることから考えると観測と予測値は

明るさは見かけの断面積に大きく依存する．しかしな一致していると考えて良い．

がら，小惑星は一般的に形状が歪であり，自転軸や絶そこで，ベスタも含めた物理的素性のよく分かって対的な大きさを正確に測ることが困難な為に，精度高いる 55 個の小惑星を標準星として用いて，筆者を含く明るさを予測することは難しい様に思われる．但し， めた研究グループでは，あかり赤外線天文衛星の遠赤ベスタのように形状や絶対的な大きさ・自転ベクトル外サーベイの恒星カタログ [38] の較正，及び，小惑星が分かっている小惑星は地球から見た断面積が精度良カタログ [39] の較正を行った．異なる手法で較正されく予測できる為に，高い精度で中間・遠赤外線域の明た IRAS 赤外線天文衛星や WISE 赤外線天文衛星で得るさを予測が可能である．

られたそれぞれの直径の値とあかり赤外線天文衛星でベスタのフラックス予測の確からしさ（大凡 5 % 程得られた小惑星のカタログの直径の値が 10% 以内の度の確からしさ）は確かめられている [e.g., 36] が，自精度で一致していることが確認されている [40] ．この■2014遊星人Vol23-4_3校.indd363

2014/12/13

14:04:54

364

日本惑星科学会誌 Vol. 23, No. 4, 2014





    





    






 


 


 












    

    





 



 

 

 

 

 

 
 

 

 

 
















    



図2： 過去のベスタの3ミクロン帯の観測． △は[42]，○は[43]のデータ．これらのデータだけみると，3ミクロン帯の観測が無い様に見えるが，◆の[44]と＋の[45]のデータと組み合わせると，3ミクロン帯の吸収があるように見える．

図3： ベスタの3ミクロン帯の観測．黒線が2003年3月1日，赤線
が2003年3月2日のスペクトルデータ． 3月2日のデータに3ミクロン帯の吸収があることがわかる．

ことはあかり赤外線天文衛星で用いた較正法の正しさ介した様に分化した小惑星として知られていたこともを示している．

あり，3 ミクロン帯の吸収があることは『この当時』はこのように，小惑星ベスタは明るさの絶対値を決定あり得ないことと暗黙の内に考えられていた．それ故，する為に，非常に有用な小惑星の１つであるというこ見いだした結果はそれと矛盾する結果であった．但し，とが出来る．それを裏付けように，ベスタはあかり赤

組み合わせたデータには観測誤差があり，また，ベス外線天文衛星以降の天文衛星である Herschel 宇宙望

タの過去の 3 ミクロン観測が対象波長帯を同時刻に取遠鏡の 3 つの機器の PACS・SPIRE・HIFI の絶対値較

得できない多色測光のデータの為に，ベスタ自身の変正にも使用されている [41]．

光度（可視光域で最大差で約 0.2 等くらい）の影響も捨てきれなかった．この問題を解決する為には 3 ミクロ

3. 3ミクロン帯観測に関わる研究

ン帯の分光観測を行う以外は方法が無かった．
そこで，筆者は 1.9～3.5 ミクロンの波長帯での分光

筆者は小惑星の含水鉱物の将来の研究（あかり赤外観測を 2003 年にハワイ・マウナケアにある英国合同

線天文衛星による近赤外分光サーベイ観測計画立案の天文センターの 3.8mUKIRT 望遠鏡の CGS4 装置を使準備）の為に，その時点で公開されていた 3 ミクロン

用して行った．観測の結果，大凡数 % 程度の吸収が帯の多色測光・分光データを中心に検証を行っていた．

あることが判明した（図 3）．この結果はすぐさま雑誌

検証のポイントとしては，単純に１つの論文のデータに投稿したが，吸収レベルが % レベルと言うこととを眺めるだけでなく，異なる波長域の多色測光・分光

上述の様な理由から，なかなかすんなりとは受理はさ

データを組み合わせて，3 ミクロン帯の吸収の有無をれなかったが，観測から半年程度で世の中に公開され調べたことである．ベスタの 3 ミクロン帯の観測は 2た [46]．

つの論文 [42, 43] を単独に見てみると（図 2 の△印と○3 ミクロン帯の吸収としては成因がいくつか考えら

印をそれぞれ単独にみると），それぞれ 3 ミクロン帯

れた．１つめとしては太陽風による陽子の打ち込みにの吸収は無い様に見える．しかしながら，2.5 ミクロ

よるベスタ表層上で含水鉱物の生成説 [46]，２つめはンまでの連続的なスペクトル [44, 45] を組み合わせて

炭素質コンドライトのような天体衝突によってもたらみると，明らかに少なくとも 5 % 程の 3 ミクロン帯のされた外因説 [46]，３つ目は炭素質コンドライト天体

吸収があることが判明した（図 2）．ベスタは序章で紹族起源の惑星間塵がベスタに降り積もったという外因■2014遊星人Vol23-4_3校.indd364

2014/12/13

14:04:54

天文観測的手法における小惑星4ベスタの研究／長谷川365

説が考えられた [47]．論文を書いた当初では３つの説い場所が所々に存在し，それらは炭素質コンドライトの内，筆者は 2・3 つめの説が有力と考えていた．当の衝突の結果と考えられている [55, 56]．そのような時は，他の天体上で陽子の打ち込みによる含水鉱物の衝突の繰り返しで，宇宙風化の蓄積が無くなっている

生 成 が 確 認 さ れ て い な い こ と， ホ ワ ル ダ イ ト 内 にとも考えられている．

CM2 コンドライトの欠片が見つかっていること，ベ

但し、Dawn 探査機は磁力計を搭載していなかった

スタでの衝突速度を考えると衝突圧力的に含水鉱物が為に，ベスタに磁場が存在しているかどうかは現時点

消えないことからである．しかしながら，10 年経っ

では実際の所はっきりと確定していない．よって，少た現在では，研究の進歩があり，必ずしもそうとは言なくとも含水鉱物の太陽風成因説をはっきりと否定す

えない状況にある．

る事ができない．3 ミクロン帯のスペクトルの形から１つ目の重要な研究は 2009 年に，月での含水鉱物

その起源を考えると，ベスタの磁場が太陽風に影響をの発見である [48, 49]．元々月での含水鉱物の存在は与えるほどは存在していないという説に一票を投じる

月起源隕石での炭素質コンドライトの発見から推定しことになるかと思われる．

ていた [46] が，成因としてはそれではなく，磁場のな
い月表面に太陽風によって陽子が打ち込まれたことが成因として考えられた．その後の研究で，3 ミクロン

4. ライトカーブ観測によるベスタの巨大衝突に関わる研究

帯の含水鉱物の吸収の形が，蛇紋石やスメクタイトと序章でも述べたが，ベスタは族を形成し，その近傍言った層状珪酸塩鉱物の吸収の位置（ピークが 2.7 ミクには V 型小惑星が多く存在しているが，それら V 型小ロン付近）が異なり，太陽風による陽子の打ち込みに惑星はベスタで起きた巨大衝突によってベスタから飛よって人工的に作った鉱物と吸収の位置（ピークが 2.8

び出してきた破片と考えられている [e.g., 29]．よって，ミクロン付近）が一致していること [50] から，太陽風破片である V 型小惑星の研究を行うことによって，ベ起源説が決定的であると考えられている．２つ目の重スタの衝突に関する情報が得られる可能性がある．そ要な研究は，Dawn 探査機によって，3 ミクロン帯のこで筆者は，ベスタが存在している内側小惑星帯に数

スペクトルが取得されたことである [51]．ちなみにだ

多く存在している V 型小惑星のライトカーブ観測をが，実は 2003 年の筆者の研究結果はそれまでは色眼行った [57]．

鏡で見られていた感はあったが，Dawn 探査機によって，その正しさが証明されたことになった [52]．

2003～2005 年の間に東京大学木曽観測所の 1.05 m
シュミット望遠鏡と 0.3 m 望遠鏡，宮坂天文台の 0.25

地上望遠鏡で大気吸収がある故に連続的にスペクト

m/0.36 m 望遠鏡，紀美野町みさと天文台の 1.05 m 望ルが取得出来ないないが，Dawn 探査機での観測は真遠鏡，国立天文台岡山天体物理観測所の 0.5 m 望遠鏡，

空中下で行われたので，3 ミクロン帯のスペクトルがハワイ・マウナケアにあるハワイ大の 2.24 m 望遠鏡

連測的に取得された．その結果を見ると，実は月の 3の 7 台の望遠鏡を用いて，22 個の V 型小惑星のライト

ミクロン帯のスペクトルと同様に太陽風による陽子のカーブを取得し，その内 19 個の天体の周期が判明し打ち込みによって人工的に作った鉱物と吸収に位置がた（表 2）．なお，副産物として，V 型小惑星を観測し一致していることが分かった．このことはベスタ表層た視野に写った V 型以外の 15 個の小惑星が検出され，の 3 ミクロン帯の吸収の成因として，太陽風による陽そのうち，13 個の自転周期が判明した．自ら取得し子の打ち込みもその候補として挙げることができる．

た 19 個の V 型小惑星の自転周期に加え，ライトカー

しかし一方で，Dawn 探査機に以外の 2 つの研究：

ブデータベース [58] に公開されている V 型小惑星のデベスタ表層の可視・近赤外のスペクトルが赤化していータを集めて，総計 59 天体の V 型小惑星の角速度分

な い こ と か ら の 磁 場 の 存 在 を 示 唆 す る 研 究 [53] と

布を調べた
（図 4）
．

HED の磁場測定から現在も磁場が存在する可能性を

小惑星の角速度分布は特に直径 40 km 以上の天体示唆する研究 [54] からは，ベスタ表層の太陽風の影響

ではマクスウェル分布と一致することが知られているは否定されている．また，Dawn 探査機の結果から，[59]．衝突履歴があるとマクスウェル分布になるとい

ベスタ表面には可視光で見ると局所的にアルベドの暗う理論的な研究 [60] から，直径 40 km 以上の天体では■2014遊星人Vol23-4_3校.indd

365

2014/12/13

14:04:54

366

日本惑星科学会誌 Vol. 23, No. 4, 2014

体は YORP 効果が効くとされる 40 km 以下で構成されており，また，衝突破壊年代もそれぞれ 2.5 Gyr, 3  



Gyr, 1 Gy とされており，YORP 効果が十分働くだけのタイムスケールである．



他の族と要因から考えると，V 型小惑星の角速度分布

V 型小惑星の殆どが直径 40 km という事実と前述ののマクスウェル分布と一致しない結果の要因は Gy く

らいのタイムスケールで YORP 効果が効いた結果と推定することができる．



ベスタ上で大規模な衝突破壊が何時行われたかとい










     

図4： 内側小惑星帯のV型小惑星の角速度分布．実線が観測デー
タ．破線が観測データをマクスウェル分布でフィットした
線．一致していないことがわかる．

う研究は幾つかある．1996 年にベスタ族の数値シミュレーション研究が行われ，族形成年代が大凡 1 Gy前とされた [64]．2003 年に小惑星帯にあるダストバン
ドは 5-250 Myr 前に衝突破壊された族がその源とされたが，ベスタ族にはそれに対応するダストバンドがな
いので，それより古いと考えられている [65]．2005 年

表2： 観測したV型小惑星．副次的に観測された小惑星はスラッ
シュの後に示す．

望遠鏡
Kiso 1.05m

V型小惑星/副次的に観測された小惑星
1933 2011 2508 2511 2640 2795 3307
3657 3900 4005 4147 4977 8645 10285
10320 / 477 1455 3192 6664 10389
10443 11321 18950 22034 29976 41051
46121

Kiso 0.3m

2511 2653 4796 6331

Miyasaka 0.25m

2653

Miyasaka 0.36m

1933 2011 3900 4005 4434 4796 6331
8645 / 477 11321 89481

Misato 1.05m

4383 4796

Okayama 0.5m

6331

UH 2.28m

4383

天体衝突の履歴が残っているということが知られていに小惑星帯のサイズ分布の再現研究から約 20 % の確率で 3.5 Gyr 前にベスタ族が形成されたとされた [66]．
2007 年にベスタに近接する小惑星の軌道シミュレー
ションを行い族形成年代を大凡 1.2 Gyr 前とした [67]．
2008 年に軌道進化シミュレーションから少なくとも 1
Gy 前にベスタ族が形成されているとした [68]．Dawn
探査機で観測された南半球の巨大クレーターのクレーター個数の観測から，２つの巨大クレーターが 1 Gyr前に形成されたと見積もった [69]．これらの結果は内側小惑星帯の V 型小惑星が Gy くらいの年代で現在のサイズの小惑星となり，その結果 YORP 効果が働いたという結論と矛盾しない．

5. 衝効果に関わる研究る．しかしながら，V 型小惑星の角速度分布はマクス
ウェル分布と Kolmogorov-Smirnov 検定で凡そ 90% の位相角が 0 度に近づいた時に急激に明るくなる現象

確率で一致しない結果になった．このことから，これ「衝効果」は 1895 年に土星のリングで確認されており

ら V 型小惑星達が衝突時に獲得した衝突履歴が消え[70]，1922 年に月 [71] で，1967 年に火星で確認されててしまったことが考えられる．

いる [72]．小惑星での衝効果の現象は 1956 年に確認さ他の族の角速度分布と比較してみると，コロニス族

れ [73]，ベスタ自身も 1967 年に確認されている [74]．

[61]，フローラ族 [62]，マリナ族 [63] も同様にマクスウ

衝効果の原因として影効果と干渉性後方散乱の２つ

ェル分布と一致しない結果になっている．それぞれのが考えられている [e.g., 74]．影効果は表面にある凸凹論文で，マクスウェル分布と一致しない理由として，が明るさに寄与するという効果であり，小惑星では位

YORP（Yarkovsky–O'Keefe–Radzievskii–Paddack）効

相角が 7,8 度から現れる現象である．干渉性後方散乱

果によって，角速度分布がマクスウェル分布から外れは表層が明るく，かつ，位相角が 1～2 度以内の場合たと述べられている．実際に，これらの族の多くの天に現れる現象である [75]．衝効果による干渉性後方散■2014遊星人Vol23-4_3校.indd366

2014/12/13

14:04:54

天文観測的手法における小惑星4ベスタの研究／長谷川

367


表3： ベスタの測光観測と分光観測．



 





望遠鏡
位相角 [度]
Sagamihara 0.064m 0.2, 0.3, 0.8, 0.9, 1.1, 1.2, 1.3, 3.1, 3.2,
4.1, 4.2, 6.0, 6.1, 6.8, 6.9, 7.0, 9.4, 9.5,
10.2, 10.3, 10.4, 12.4, 12.5, 12.8, 12.9,
15.3, 22.6, 22.7, 22.8, 23.2, 23.7, 23.8
Miyasaka 0.36m
0.1, 0.2, 0.4, 0.8, 0.9, 1.1, 1.2, 1.3, 1.4,
1.6, 1.7, 2.6, 2.7, 10.3, 10.6, 10.7, 10.8,
13.5, 13.6, 16.0, 17.8, 17.9, 21.0, 21.1,
21.2
Nishiharima 0.076m 0.2, 0.4, 18.1, 19.4, 23.2, 23.7, 23.8
Maidanak 0.6m
0.1
Nishiharima 2.0m
20.3, 20.9
Okayama 1.88m
24.9














波長帯
B（CB有）
B（CB無）
Rc（CB有）Rc（CB無）
z’
（CB有）
z’
（CB無）

アルベド
0.353
0.321
0.407
0.335
0.305
0.298

空隙率
0.575
0.7
0.4
0.5
0.65
0.7

密度 [kg/m3]
1300
900
1900
1600
1100
900



 

表4： ベスタの反射率とHapkeモデルによる空隙率および密度．

 

図5： ベスタのＲバンドの位相関数．点は観測データ，点線はフィットした影効果のみのHapkeモデル，実線はフィッ
トした影効果と干渉性後方散乱のHapkeモデル，点線はフィットした影効果と干渉性後方散乱（Bc0=0）のHapkeモ
デル．影効果と干渉性後方散乱のHapkeモデルが良く一致しているのがわかる．
ɇ\wȧąʱȧ̑ĊƝɗɮŦŭʫ~ŧǖķ}?{qǮˤ˭yŏɮ
ɇ\wȧąʱȧ̑ĊƝɗɮŦŭʫ~ŧǖķ}?{qǮˤ˭yŏɮ ̒のŧǖķƈƇŧóȶɃʱȧƽ}
3 バンドで行った．位相角 ?Ǯˤ˭ötwĔąʱȧʣtq01̎ʧ̏
0.1～23.8 度の間で，測光ŧǖķƈƇŧóȶɃʱȧƽ}
データを取得したが，Rosetta?Ǯˤ˭ötwĔąʱȧʣtq01̎ʧ̏
探査機と Dawn 探査機で，ȧąʱȧ~'5K}µÐ°xʣtq̒ðɔʲ̚  ƛ}˰x̑ȧą®Ò̒

乱の有無を調べる為には位相角が 1 度以内の観測が必ɇ\wȧąʱȧ̑ĊƝɗɮŦŭʫ~ŧǖķ}?{qǮˤ˭yŏɮȧąʱȧ~'5K}µÐ°xʣtq̒ðɔʲ̚  ƛ}˰x̑ȧą®Ò須である．但し，ベスタの低位相角の観測は位相角がĲƫiw\}x̑o}®Ò«ðɔ˱ǔżē|ɇ\q01̒ȧąʱȧ~'5K}µÐ°xʣtq̒ðɔʲ̚  ƛ}˰x̑ȧą®Ò
それらのデータも位相関数導出に用いた[77, 78]．

1.2 度までの観測しか無く，そもそも 1.2 度未満での位相関数の変化が分からなかった．序章でも述べたが，ベスタはあらゆる観測手法で観測が行われ，Dawn 探査機のランデブーも行われ，それらによってベスタの«Ĳƫjqȃ'AD7EE3ǉǻȉy3H@ǉǻȉx̑ðɔʲƛˍ^®Ò«
位相角30 度を超えるデータを取得されているので，
?Ǯˤ˭ötwĔąʱȧʣtq01̎ʧ̏̒

ŧǖķƈƇŧóȶɃʱȧƽ}«Ĳƫjqȃ'AD7EE3ǉǻȉy3H@ǉǻȉx̑ðɔʲ ƛˍ^®Ò«Ĳƫiw\}x̑o}®Ò«ðɔ˱ǔżē|ɇ\q01̒«Ĳƫjqȃ'AD7EE3ǉǻȉy3H@ǉǻȉx̑ðɔʲƛˍ^®Ò« ʱȧ|twƫqʄǸ~Ŏ̗|ɣl̒ȶɃƵŝ}Ĳƫ}ȱ|ƔȦƳƪǚ

観測によって得られた結果は図 5 に示す．物理情報Ĳƫiw\}x̑o}®Ò«ðɔ˱ǔżē|ɇ\q01̒ ʱȧ|twƫqʄǸ~Ŏ̗|ɣl̒ȶɃƵŝ}Ĳƫ}ȱ|ƔȦƳƪǚ

ǓÝaʑƺiq3B=7Ä®Ìx¹­¯jw\07 9 1̒ǓÝaʑƺiq3B=7Ä®Ìx¹­¯jw\07
9 1̒


Hapke モ
の取得の為に干渉性後方散乱が考慮された
ʱȧ|twƫqʄǸ~Ŏ̗|ɣl̒ȶɃƵŝ}Ĳƫ}ȱ|ƔȦƳƪǚ
デルでフィットしている[e.g., 79]．
ǓÝaʑƺiq3B=7Ä®Ìx¹­¯jw\07 9 1̒


知りうる物理情報は殆ど分かっているが，実は 1 度以下の低位相角での位相関数の振る舞いが実は判明して

ここで，wλは単粒子散乱能，Pλ (α) は Henyey-

いなかった．

 

2006 年初頭にベスタの位相角が 0.1 度になる事に気がついた筆者はベスタの低位相角から高位相角まで広hhx̑
wλ~ĨɶŪǓÝʔ̑Pλ(α)はHenyey-Greenstein位相関数̑Kλ(α,θλ)はhhx̑wλ~ĨɶŪǓÝʔ̑Pλ(α)はHenyey-Greenstein位相関数̑Kλ(α,θλ)はGreenstein 位相関数，K
λ (α, θλ) は表面の凸凹に関する
表面の凸凹に関する補正項̑
r0λ~×ʶ}Ơxɣi̒hhx̑
wλ~ĨɶŪǓÝʔ̑Pr0λ(α)はHenyey-Greenstein位相関数̑Kλ(α,θλ)は表面の凸凹に関する補正項̑
λ~×ʶ}Ơxɣi̒補正項，r0λは下記の式で示される．
表面の凸凹に関する補正項̑
r0λ~×ʶ}Ơxɣi̒

い範囲の位相関数を取得する為の観測を行った．なお，

 BSHλ(α)yBCBλ(α)は×ʶ}Ơxɣi̒

ベスタは軌道傾斜角が地球の軌道傾斜角と異なっている為に，位相角が 0.1 度近くになるのは百年に 1 回程 ×ʶ}Ơxɣi̒ BSHλ(α)yBCBλ(α)はSHλ(α)(α)とyB
BCBλ(α)は×ʶ}Ơxɣi̒CBλ (α) は下記の式で示される．
BBSHλ


しかない．



JAXA 宇宙科学研究所屋上に仮設置した 0.064 m 望

遠鏡と宮坂天文台の 0.36 m 望遠鏡，兵庫県立大学西 BS0λyhSλ~ƧĠǸ}́}ǅƓyĦûƓx[̑BC0λyhCλ~ƔȦƳƪǚǓÝ}́

はりま天文台の 0.6 m 望遠鏡にガイド望遠鏡として同}ǅƓyĦûƓx[̒

架されている 0.076 m 望遠鏡，ウズベキスタン・マイ¹­¯jwƫq¶ÊÃÒ«`m̑3B=7Ä®Ìɬ˷ȾȘh

ダナク天文台にある 0.6 m 望遠鏡を用いて測光観測を，

yaxb̒

兵庫県立大学西はりま天文台の 2.0 m なゆた望遠鏡とBS0λyhSλ~ƧĠǸ}́}ǅƓyĦûƓx[̑BC0λyhCλ~ƔȦƳƪǚǓÝ}́
B}ǅƓyĦûƓx[̒C 0λ と hCλ は干渉性後方散乱の項の振幅と半値幅であ

国立天文台岡山天体物理観測所の 1.88 m 望遠鏡を使って分光観測を行った [76]（表 3）．測光観測は B, Rc, z'BS 0λ と hSλ は影効果の項の振幅と半値幅であり，
  BS0λyhSλ~ƧĠǸ}́}ǅƓyĦûƓx[̑BC0λyhCλ~ƔȦƳƪǚǓÝ}́
}ǅƓyĦûƓx[̒¹­¯jwƫq¶ÊÃÒ«`m̑3B=7Ä®Ìɬ˷ȾȘhる．
¹­¯jwƫq¶ÊÃÒ«`m̑3B=7Ä®Ìɬ˷ȾȘhyaxb̒
フィットして得られたパラメータからまず，Hapke
yaxb̒


■2014遊星人Vol23-4_3校.indd367

2014/12/13

14:04:55 BS0λyhSλ~ƧĠǸ}́}ǅƓyĦûƓx[̑BC0λyhCλ~ƔȦƳƪǚǓÝ}́
}ǅƓyĦûƓx[̒368

日本惑星科学会誌 Vol. 23, No. 4, 2014¹­¯jwƫq¶ÊÃÒ«`m̑3B=7Ä®Ìɬ˷ȾȘhyaxb̒

を行ったアイデア勝負の王道から外れた研究であったモデルより空隙率を求めることができる．
と言えよう．しかしながら，これら研究がベスタの物
理的性質解明に一役買ったということは言うことがで

ρλは空隙率，dl/ds は構成粒子の最大・最小直径比できるかと思う．

ある．空隙率と構成粒子の最大・最小直径比を用いると表層レゴリス層の空隙率を求めることが出来る．

謝 辞

各バンドで異なるが 0.4～0.7 の間での空隙率であった（表 4）．ベスタのレゴリス層の真密度をユークライ
トと同程度と仮定するとその密度は 900 － 1900[kg m-

3] となった．

本研究は様々な方にサポートされて実行された．
標準星の研究を進めるにあたっては，藤原英明氏，
片 坐 宏 一 氏， 藤 吉 拓 哉 氏，T. G. Mueller 氏，M.

各バンドのアルベドが下がると，各バンドのバルク
密度が小さくなることもわかった（表 4）．レゴリスのバルク密度は表面に近いほど，小さくなることが分か

Cohen 氏，R. Moreno 氏 , 大坪貴文氏，関口朋彦氏，山村一誠氏にお世話になった．
3 ミクロンの研究を進めるにあたっては村川幸史氏，っている [80]．よって，アルベドが低いバンドでのバ

石 黒 正 晃 氏， 野 中 秀 紀 氏， 高 遠 徳 尚 氏，Chris J.

ルク密度の情報は，より表面付近のバルク密度を示し

Davis 氏，上野宗孝氏，廣井孝弘氏にお世話になった．

ている可能性がある．

V 型小惑星の研究を進めるにあたっては，宮坂正大氏，

次に衝効果の成因について探った．衝効果で，影効三戸洋之氏，猿楽祐樹氏，小澤友彦氏，黒田大介氏，

果のみ・影効果と干渉性後方散乱の 2 パターンでフィ吉田道利氏，柳澤顕史氏，清水康広氏，長山省吾氏，ットをおこなった．その結果，影効果と干渉性後方散戸田博之氏，沖田喜一氏，河合誠之氏，北里宏平氏，

乱でフィットした方が良く一致していることが分かり， 吉住千亜紀氏，西原説子氏，縫田明理氏，森真知子氏，
影効果のみだと，特に位相角１度未満のデータを説明

関口朋彦氏，石黒正晃氏，阿部琢美氏，安部正真氏にできなかった．このことからベスタは位相角 1 度未満

お世話になった．

では干渉性後方散乱の効果が色濃く出ることが分かった．

衝効果の研究を進めるにあたっては，宮坂正大氏，
時政典孝氏，十亀昭人氏，M. A. Ibrahimov 氏，吉田干渉性後方散乱の強度はベスタの場合はアルベドと相関があることが分かった．このことは干渉性後方散乱が，アルベドが高い粒子，即ち，透明度の高く，反射率の高い粒子の多重散乱で引き起こされている事を示している．

二美氏，尾崎忍夫氏，石黒正晃氏，安部正真氏，黒田
大介氏にお世話になった．
過去の入手困難な文献検索と本原稿執筆と校正において，臼井文彦氏には大変お世話になった．
査読者として，原稿を丁寧に読んでいただき，大変有益なコメントをくださった廣井孝弘氏に心より感謝いたします．

6. まとめ

本研究は JAXA あかり赤外線天文衛星と国立天文
筆者は天文観測的手法観測的手法で，
1） 3 ミクロン帯の分光観測から，ベスタ表層に含水鉱物の吸収があることを発見した．
2） V 型小惑星の観測で，それらが破片として，ベス
タから飛び出した年代を示した．
3） 位相関数の観測から，ベスタ表層の密度を示し，

台すばる望遠鏡によって実施された観測結果を使用した．英国合同天文センターと国立天文台との国際協力の枠組みによって英国 UKIRT 望遠鏡を使用させていただいた．また，宮坂天文台，東京大学木曽観測所，
国立天文台岡山天体物理観測所，兵庫県立大学西はり
ま天文台，紀美野町みさと天文台，ハワイ大学天文学

衝効果の原因を明らかにした．

研究所，ウルグベク天文研究所マイダナク天文台の望筆者が行ったこれらの研究結果は，最新の観測技術

遠鏡を使用して，観測が行われた．

を用いての王道の研究というよりかは，既知のデータから未知であることを探しだし，それに対して，観測■2014遊星人Vol23-4_3校.indd368

本研究は JAXA 宇宙科学研究所スペースプラズマ
共同利用のサポートを受けている．

2014/12/13

14:04:55

天文観測的手法における小惑星4ベスタの研究／長谷川369(Houston: Lunar and Planetary Institute), 7.

参考文献[29] Thomas, P.C. et al., 1997, Science 277, 1492.
[30] Müller, T.G. & Lagerros, J.S.V., 1998, A&A 338, 340.

[1] Aspin, J. 1825, A Familiar Treatise on Astronomy:[31] Li, J.-Y. et al., 2010, Icarus 208, 238.

Explaining the General Phenomena of Celestial Bodies[32] Li, J.-Y. et al., 2011, Icarus 216, 640.

(London: Samuel Leigh).

[33] Gundlach, B. & Blum, J., 2013, Icarus 223, 479.

[2] Harrington, M.W., 1883, Amer. J. Sci. 26, 464.

[34] Russell, C.T. et al., 2012, Science 336, 684.

[3] Bobrovnikoff, N.T., 1929, Lick Obs. bull. 407, 18.

[35] Landolt, A.U., 1992, AJ 104, 340.

[4] Lyot, B., 1934, Bull. Astron. 67, 3.

[36] Müller, T.G. & Lagerros, J.S.V., 2002, A&A 381, 324.

[5] Stephenson, C.B., 1951, AJ 114, 500.

[37] Fujiwara, H. et al., Private Comm..

[6] Gehrels, T., 1967, AJ 72, 929.

[38] Yamamura, I. et al., 2010, http://www.ir.isas.jaxa.jp/[7] Hertz, H.G., 1968, Science 160, 299.
[8] Allen, D.A., 1970, Nature 227, 158.

ASTRO-F/Observation/PSC/Public/RN/AKARI-FIS_
BSC_V1_RN.pdf.

[9] McCord, T. et al., 1970, Science 168, 1445.

[39] Usui, F. et al., 2011, PASJ 63, 1117.

[10] Taylor, R.C., 1973, AJ 78, 1131.

[40] Usui, F. et al., 2014, PASJ 66, 56.

[11] Larson, H.P. & Fink, U., 1975, Icarus 26, 420.

[41] Müller, T.G. et al., 2013, Exp. Astron. 37, 253.

[12] Johnson, T.V. et al., 1975, ApJ 197, 527.

[42] Lebofsky, L.A., 1980, AJ, 85, 573.

[13] Worden, S.P. et al., 1977, Icarus 32, 450.

[43] Feierberg, M.A. et al., 1985. Icarus 63, 183[14] Conklin, E.K. et al., 1777, in: Comets, asteroids,[44] Bell, J.F. et al., 1988, Proc. Lunar Planet. Sci. Conf.

meteorites (Toledo: University of Toledo), 261.
[15] Blanco, C. & Catalano, S., 1979, Icarus 40, 359.
[16] Degewij, J. et al., 1979, Icarus 40, 364.
[17] Dickel, J.R., 1979, in Asteroids, ed. T. Gehrels(Tucson: University of Arizona Press), 212.
[18] Williams, J.G., 1979, in Asteroids, ed. T. Gehrels(Tucson: University of Arizona Press), 1040.

XIX, 57.
[45] Le MouÈlic, S. et al., 2001, Proc. Lunar Planet. Sci.
Conf. XXXII, 1094.
[46] Hasegawa, S. et al., 2003, Geophys. Res. Lett. 30,
2123[47] Hasegawa, S. et al., 2004, Proc. Lunar Planet. Sci.
Conf. XXXV, 1458.

[19] Butterworth, P.S. et al., 1980, Nature 287, 701.

[48] Clark, R.N., 2009, Science 326, 562.

[20] Matson, D.L., 1986, Infrared Astronomical Satellite[49] Pieters, C.M. et al., 2009, Science 326, 568.

Asteroid and Comet Survey (Pasadena: Jet Propulsion[50] Ichimura, A.S. et al., 2012, EPSL 345-348, 90.

Laboratory, California Institute of Technology).

[51] De Sanctis M.C. et al., 2013, Metrotit. Planet. Sci. 48,[21] Drummond, J.D. et al., 1988, Icarus 73, 1.

2166.

[22] Johnston, K.J. et al., 1989, AJ 98, 335.

[52] Clark, B.E., 2012, Nature 491, 45.

[23] Zappalà, V. et al., 1990b, AJ 100, 2030.

[53] Marchi, S. et al., 2010, ApJ 721, L172.

[24] Cruikshank, D.P. et al., 1991, Icarus 89, 1.

[54] Fu, R.R. et al., 2012, Science 338, 238.

[25] Dunham, D.W. et al., 1991, in: Comets, asteroids,[55] McCord, T.B., 2012, Nature 491, 83.

meteorites 1991 (Houston, Lunar and Planetary[56] Reddy, V. et al., 2012, Icarus 221, 544.

Institute), 54.

[57] Hasegawa, S. et al., 2014a, PASJ 66, 53.

[26] Binzel, R.P. & Xu S., 1993, Science 260, 186.

[58] Warner, B.D. et al., 2009, Icarus 202, 134.

[27] Mitchell, D.L. et al., 1996, Icarus 124, 113.

[59] Pravec, P. et al., in Asteroids III, ed. W.F. Bottke,[28] Dumas, C. & Hainaut, O.R., 2006, in Workshop on
Evolution of Igneous Asteroids: Focus on Vesta and the
HED Meterorites, ed. D. W. Mittlefehldt & J. J. Papike■2014遊星人Vol23-4_3校.indd

369

A. Cellino, P. Paolicchi, & R.P. Binzel (Tucson:
University of Arizona Press), 113.
[60] Salo, H., 1987, Icarus 70, 37.

2014/12/13

14:04:55

370

日本惑星科学会誌 Vol. 23, No. 4, 2014[61] Slivan, S.M. et al., 2008, Icarus 195, 226.
[62] Kryszczynska,́
A. et al., 2012, A&A 546, A72.
[63] Kim, M.-J. et al., 2014, AJ 147, 56.
[64] Marzari, F. et al., 1996, A&A 316, 248.
[65] Nesvorný, D. et al., 2003 ApJ 591, 486.
[66] Bottke, W.F. et al., 2005, Icarus 175, 111.
[67] Carruba, V. et al., 2007, A&A 465, 315.
[68] Nesvorný, D. et al., 2008, Icarus 193, 85.
[69] Marchi, S. et al., 2012, Science 336, 690.
[70] Seeliger, H., 1895, Abh. Bayer. Acad. Wiss. Math.
Naturwiss. Kl. 18, 1.
[71] Barabashev, N.P., 1922, Astron. Nachr. 217, 445.
[72] O'Leary, B.T., 1967, AJ 149, 147.
[73] Gehrels, T., 1956, ApJ 123, 331.
[74] Muinonen, K. et al., 2002, in Asteroids III, ed. W.F.
Bottke, A. Cellino, P. Paolicchi, & R.P. Binzel (Tucson:
University of Arizona Press), 123.
[75] Mishchenko, M.I. et al., 2006, Appl. Opt. 45, 4459.
[76] Hasegawa, S. et al., 2014b, PASJ 66, 89.
[77] Fornasier, S. et al., 2011, A&A 53, L9.
[78] Li, J.-Y. et al., 2013, Icarus 226, 1252.
[79] Hapke, B., 2002, Icarus 157, 523.
[80] Houston, W.N. et al., 1974, Proc. ﬁfth Lunar Planet.
Sci. conf. Suppl.5, Geochim. Cosmochim. Acta 3,
2361.

■2014遊星人Vol23-4_3校.indd

370

2014/12/13

14:04:55