137



地球に於ける氷期・間氷期の繰り返しが長い時間スケー

地球の軌道要素の変化と自転軸の歳差による極域へのルでの日射量の変化に駆動されていることは広く知ら日射量の変動が要因しているのではないかという説がれている. この日射量変化は惑星間の重力相互作用に形成されて行った[2]. ひどく大雑把な言い方をすれ

よる地球の軌道要素変化と自転軸の運動によって持たば, 極域への日射量が増えれば氷床は縮小し, 極域へらされるが, この原理は火星に於いても同様に機能し,の日射量が減れば氷床は拡大するであろうという予想

火星の大気上端への日射量変動を引き起こしているはである. 地球上では白い氷床が太陽光を反射してしまずである. 地球に於ける日射量変動が氷床の存在を介うことによりますます気温を下げてしまう効果 (｢氷して気候システムに大きな影響を与えて来たように,床・アルベドフィードバック｣ と呼ばれる) が働くこ

火星に於いても日射量変動が氷床 (極冠) を介して気ともあり, 日射量変動と氷床変動の相関可能性はかな

候変動を支配している可能性がある. 本稿では惑星のり高いと思われていた.

公転運動と自転運動が原因する日射量変動研究の現状

こうした一連の定性的議論を詳細な計算により初めを概観し, それらが火星の気候研究にとってどのようて定量化したのがセルビアの科学者 M. Milankovitch

な意味を持ち得るのかを考えてみたい.

である. Milankovitch は地球へ入射する日射量の緯度分布と季節変化を当時としては驚くべき高精度で計算

し, また, 当時得られていた最高精度の公転軌道変化理論を用いて日射量の長周期変化を詳しく計算した[3].

地球の話から始めよう. 現在から遡ること過去約二Milankovitchの研究の大半は1920年代から1930年代に百万年の時代は第四紀と呼ばれている. 第四紀の気候かけて遂行され, 1950年代辺りまでは地質学的なデーは氷期・間氷期サイクルという周期変動により特徴付タを解読するための時間目盛としてその結果を使おう

けられていると言って良い. 氷期とは極域に大きな氷とする地質学者も存在した. 当時は現在のように放射床が存在して全惑星の気候を寒冷たらしめる時期であ性同位元素による精密な年代測定法が確立しておらず,り, 間氷期とは氷床が縮小して気候がやや温暖になる

氷期に関連する地質データを解釈する上での年代尺度時期である. 第四紀に於いては氷期と間氷期が数万年が Milankovitch の日射量変動の理論計算値くらいしかから十万年の時間スケールで交互に訪れ, 地球の寒暖なかったからである. 当時のやり方は地質学的なデーを司って来た[1]. この繰り返しを本稿では氷期・間タから推測される気候条件の変化を Milankovitch の日氷期サイクルと呼ぶことにする. 第四紀の氷期・間氷射量変動理論計算値と定性的に比較して, データの年期サイクルの存在は19世紀に於いて既に知られていた.

代を推算するというものであった. 当然予想されるよその原因については19世紀後半から20世紀前半にかけ

うにこの方法による年代推定の誤差や不確定性は大きて J. Adhemar, J. Croll, L. Pilgrim らにより議論され,

い. やがて1960年代になると放射性同位元素による年   ! "#$%&' ()%$*+%&,"%-.&%+.&,/01*2$3
4*&53+,.67&
38.$%

138

日本惑星科学会誌Vol.13.No.3,2004

代測定法がその精密さを増して来る. そうなると年代れる物理過程であり, 火星に於いても同様に発生し尺度として日射量変動の理論計算値を使う必要はもはている. 惑星の公転・自転運動の長年変化に関する研やなくなってしまう. その上に Milankovitch が計算し究は伝統的な天体力学の範疇にあり, 膨大な知識と精た日射量変動と地質学的なデータの微妙な喰い違いが密な理論体系の蓄積がある. 次節ではその一部の概略問題にされたりして, 氷期・間氷期サイクルの要因ををざっと眺め, 惑星の軌道要素の変動が日射量の変化

地球の軌道要素の変動に基づいて議論する風潮はいつに至るまでの道筋を駆け足でおさらいしてみることにの間にか消えかけてしまった[4].

したい.

けれども1970年代から1980年代にかけて海洋底や氷その前に余談として, ミランコビッチ・サイクルに床を掘削した堆積物のサンプルが採取されるようにな関してよく為される誤解について言及しておこう. そると, その中に残された気候変動指標の変化がれは ｢ミランコビッチ・サイクルとは第四紀の氷床変Milankovitch の理論的な日射量変化の周期と良く似て動そのものである｣ という説である. これはまったくいるという報告が続々となされるようになった[5].の誤りである. ミランコビッチ・サイクルとはあくま

気候変動指標とは, 例えば海洋底から取得したサンプで日射量変動の代名詞であり, 氷床変動に代表されるル中の酸素同位体比の値や浮遊性生物化石の種の構成,ような気候変動を指すものではない. 日射量変動の周あるいは氷床コア中の気泡から得られる二酸化炭素濃期は天体の公転運動と自転運動の力学から理論的に導度などである[6]. CLIMAP, SPECMAP などという大かれるものに過ぎず, それらが惑星の気候変動にどの規模プロジェクトにより得られた海底堆積物の多くかような影響を与えるかについては原理的に関与しない.

ら Milankovitch の理論と調和的な周期の気候変動指標

専門書や学術論文の中でもこの点に関して誤った記述が得られたのである. こうした結果は第四紀の氷期・を載せているものが存在するので, 読者は十分に注意間氷期サイクルが日射量変動に駆動されたものであるする必要がある.という Milankovitch 仮説の正しさを示すものとなった.
地球の気候システムを理論的にモデル化する研究も進

展し, わずか数％の日射量変動でも氷期・間氷期サイクルのような気候の大変動を十分に持たらし得るとい摂動計算などと書くと仰々しく聞こえるが, 要するう主張も強い説得力を持つようになった. こうして氷に私達の目的である惑星への日射量変動計算のために期・間氷期サイクルの主要因を惑星運動起源の日射量は, 万有引力に支配された質点系としての惑星の公転変動に求める考え方は見事に復活して一般に定着し,運動, および有限形状を持つ天体の自転軸の運動が計

このような周期的日射量変化はミランコビッチ・サイ算できれば良い. 広く知られているように天体が二個クル (Milankovitch Cycles) と呼ばれるようになった.

のみの場合, 即ち重力二体問題には厳密な解析解 (い近年では第四紀のものに限らず, 堆積物に見られる周わゆるケプラー運動) が知られているが, 天体が三個期的な環境変動の要因を解釈する場合に ｢これはミラ以上になるとそれらの運動を表す解析的な解は存在せンコビッチ・サイクルの反映ではないのか｣ という可ず, 物事は急に複雑になる. 太陽系惑星の軌道は楕円能性を考えることはごく当然となっているように思えであるとしばしば表現されるが, これも ｢近似的にはる. この辺りの事情に関心のある読者は書籍[7,8]に楕円に近い｣ というだけの話であり, 実際の惑星軌道ある記事を参照されると良いだろう.

は互いの重力の影響によって常に揺れ動いており, 一

ミランコビッチ・サイクルの原因である地球の軌道定の楕円形状を保っているわけではない. しかし太陽要素と自転軸の変化は惑星間に働く重力により支配さ系惑星の運動に於いては太陽重力の影響が他の力に比火星の日射量変動と気候／伊藤139

べて圧倒的に大きい. 太陽重力の大きさを1とすると-3非常に少なくて済む. この方法論もまた古くから知ら惑星間の重力相互作用の大きさは最大でも10 程度でれており, 永年摂動論と呼ばれている. 最先端の永年ある. これの意味するところは, 太陽系惑星の運動は摂動論の精度は非常に高く, 惑星の運動方程式を計算

ほぼ太陽重力に支配されており, それに対して惑星間機で直接解いた結果 (数値積分) とも良い一致を示しの重力相互作用が微小な摂動として働いているというている[9].

ことである. 私達は重力二体問題に関してケプラー運

ちなみにこうした計算結果がどのくらいの期間にわ動という厳密な解を知っている. 従って, 重力二体問たり有効であるのかという問いに答えるのは簡単なこ

題からさほど外れてはいない惑星運動の第一近似解ととではない. 数百万年や数億年という長期間の惑星運

してケプラー運動を採用し, それに対して惑星間の重動に関しては計算結果と比較できる観測事実も無いの力相互作用による補正を少しずつ積み重ねて行けばで, 具体的に計算結果の有効期間を示すのは難しいの(逐次近似), 最終的にはかなり正確な惑星運動の解をである. けれども, 複数の計算方法から得られる結果得ることが出来る. このような方法で ｢解けない問題｣が数億年以上にわたり良い一致を示すという事実(例の近似解を得る手法は摂動論と呼ばれ, 惑星の運動にえば[10])や, そもそも惑星運動の方程式の解はとて限らず今や物理学全般で広く使われている概念である.

も滑らかなので長期にわたり高精度での求解が可能でところで惑星の軌道要素変動は大雑把に言って二種あることなどから, 数百万年間程度での惑星の日射量類に分類される. まずその年その日の楕円軌道の上変動計算に関しては特に問題となるような計算誤差はで惑星が時々刻々と位置を変える運動. いわゆる狭い発生しないものと考えて良いだろう.

意味での公転運動であり, 強大な太陽重力に起因するさて, 惑星に入射する日射量変動を計算するには公ものである. この時間スケールは地球であれば一年,転軌道の変化を知るだけでは十分でない. 例えば地球火星であれば二年弱である. 次は惑星が日々それにの場合, 自転軸は公転軌道面に対して垂直ではなく約沿って公転する軌道自体の形状や方向が変化する運動

23.4度傾いている. また自転軸は慣性系に対して静止で, 惑星間の重力相互作用が原因している. この現象しているわけではなく, 公転軌道面の法線のまわりをの時間スケールは長く, 一般に数万年から数十万年で周期約2万6千年で周回している. これが歳差運動と呼ある. 惑星の気候変動に影響を与えると思われる日射ばれるものである. 惑星に入射する日射量の変動はこ量変動を持たらすものは後者であり, 前者に比べの自転軸の傾きと歳差運動に大きな影響を受けているて非常にゆっくりとした変化を行う. 言葉を変えれば,ので, 公転軌道要素に加えて惑星の自転軸の運動を知っ一公転周期程度では惑星の軌道形状や方向の変化は実て初めて惑星への日射量変動を計算できるようになる.に微々たるものだということである. 理科年表や天文

惑星の公転運動と同様に, 惑星自転軸の歳差運動の周年刊を見るまでもなく, 地球の近日点の位置や離心率期も惑星の自転周期や公転周期に比べてかなり長いのの値は昨年も今年も来年もほとんど同じである. 地球が普通である. 地球の歳差運動の周期2万6千年に対しの近日点が現在の場所から大きく移動するには何千年て公転周期は一年, 自転周期は一日, という具合であもの時間が必要となる. この性質に着目し, 惑星の運る. また, 惑星の形状 (正確に言えば重力ポテンシャ動方程式を公転周期について平均化するというテクニッル) は完全に球対称ではないものの, 赤道部分が微妙クが存在する. 要するにの効果を無視するのだが,に膨らんだ軸対称の回転楕円体として扱えばかなり良これをやると惑星の運動に関与する項の数が激減し,い近似になることがわかっている. 従って, 軸対称な方程式はとても簡単な形になる. こうして得られた方惑星の自転運動に関する運動方程式を立て, それを公程式を用いて摂動論的に近似解を求めれば, 計算量は転周期に関して平均して簡略化してしまえば, 公転軌

140

日本惑星科学会誌Vol.13.No.3,2004

道要素の永年摂動を計算した際と同様にして惑星自転軸の長年変動を計算することが可能になる.
以上のような手続きを踏み, 摂動論や数値積分などの解法を駆使して惑星の公転運動と自転運動の時間変化を追うことで, 私達は過去や未来の特定の時点に於ける惑星軌道の形状と方向, および自転軸の傾きと方向と知ることが可能になる. それさえわかってしまえ
ば, 後は幾何学的な関係によって惑星上の任意の地点に於ける日射量の強さと時間変化を計算することが出来る. 1930年代に Milankovitch が行ったのはこのような定量的計算であった.


惑星軌道面の形状や自転軸の方向とひと口に言っても, それを表現するのに必要な変数は数多い. しかし長い周期での日射量変動に直接影響するものはさほど

   !
"!
#$!
%&%
'()*+,-.
/0123456 7 
89:;
<"!
)=%
!
<>('<'?
?
!
@%'()*+,-2ABCD+EB F 3 ./02GHIJIFKF L 7 MFJN 

多くない. 特に重要なものは惑星軌道の離心率, 近日の夏の日射量は減り, 低緯度では日射量が増える. つ点の方向を表す角度 (近日点経度), 自転軸の傾きをまり赤道傾角項は緯度帯毎の日射量のコントラストを表す角度 (赤道傾角), それに自転軸の方向を表す角変える効果を持つのである. これに対して気候的歳差度 (ここでは歳差角と呼ぶ) である. これらの量がど項は季節のコントラストに影響を与える (図1). 自のようにして長い時間スケールでの日射量変動に関与

転軸の歳差により, 惑星が最も太陽に近くなる地点するのかを詳細に説明するには長い数式が必要となる(近日点) での自転軸の方向は刻々と移り変わる. もが, 本稿ではそれを行わずに簡略な定性的説明を試み

しも或る年に北半球の夏の時点で惑星が近日点にあれてみよう.

ば, 北半球での夏の日射量は非常に大きくなるであろ気候変動を駆動し得る日射量変動は大雑把に言ってう. 一方, その年の冬は惑星は遠日点にあり, 北半球次の二種類に分類できる (図1). 第一は自転軸の傾きが受ける日射は少なくなる. 逆の場合も同様であり,(赤道傾角) の変化による日射量変動で, 赤道傾角項

北半球の冬の時点で地球が近日点にあれば北半球での1と呼ばれる . 二番目は自転軸の方向 (歳差角) と近夏の日射量は比較的小さくなるし, 冬は近日点に当た日点の位置および離心率の変化による日射量変動で,るために日射量は比較的大きくなる. 惑星の軌道が楕

気候的歳差項と呼ばれる.

円であることと自転軸が歳差運動することにより, 季

赤道傾角項の働きは直感的にわかりやすい (図1).

節毎の日射量のコントラストが長い時間スケールで変地球の赤道傾角は主として公転軌道面の運動に依存し動するのである. 赤道傾角項と気候的歳差項の意味のて変化するが, 赤道傾角が大きい場合, 即ち地球の自より詳しい説明については[2,11]などを参照されたい.

転軸が大きく傾いている場合には高緯度の夏の日射量

こうした事柄を頭に入れて, 地球と火星の離心率,が相対的に大きくなり, 低緯度の夏の日射量が相対的

赤道傾角, 気候的歳差 (ここでは離心率×sin [近日に小さくなる. 逆に赤道傾角が小さい場合には高緯度点の方向]の値) の時間変化の計算結果とその周波数 !"#$%"&'()*+' , -./0
12%30245678火星の日射量変動と気候／伊藤141
 !" #$%&'()" *+,-.'/)01*+,-.'/)" 234
5( 67 89:;<:=>'?@A)BCDEFGHIJKLMNOBPQRSTUVWXYO"#$%&"*+,-.":=>" #$%&"
*+,-." :=>"Z[\]^_B` abcdefgB

解析結果を見てみよう (図2). 地球の離心率の振動に扁平な惑星に対して他天体から重力トルクが働くこ(図2) は約10万年と約40万年の特徴的周期を持ち,とにより発生するが, 地球以外の惑星については太陽赤道傾角には約4万年の周期が卓越している (図2とからの重力トルクが圧倒的に強い. 地球だけは例外で,). 気候的歳差は離心率と近日点経度に依存し, 地月からの重力トルクが太陽からの二倍ほどもある[12].

球の場合には2万3千年付近と1万9千年付近に強い周期

すなわち地球の自転軸は月からの強い重力トルクによ性を持つ (図2と). 火星の赤道傾角 (図2と)り激しい歳差運動を強制されていることになる. 火星と気候歳差項 (図2と) は地球のそれに比べてだにはそのような大きな衛星がないので, 歳差運動の原いぶゆっくりと変動する. 自転軸の歳差運動は重力的因となる天体は太陽のみとなる. 火星の赤道傾角と気142

日本惑星科学会誌Vol.13.No.3,2004

候的歳差の変動周期が地球のそれらに比べて随分と長一般に, 惑星間の重力相互作用により惑星の自転軸のいのはこの理由に拠る[13].

みならず惑星の公転軌道面も時々刻々と歳差する. 自以上のような軌道要素の変化から計算される日射量転軸の歳差運動の周期は惑星の力学的偏平率や月や太変動の様子が図2である. ここでは[14]の定式化を陽と言った外部天体からの重力トルクの大きさに依存

用い, 北半球の夏至に於ける北緯65度の地点での日平し, 公転軌道面の歳差運動の周期は惑星間の重力相互均日射量(W/㎡)を地球と火星のそれぞれについて計作用 (各惑星の相互配置や質量) に依存する. この両算した. 長期の気候変動という観点からは北半球で言者, すなわち惑星の自転軸の歳差周期と公転軌道面のえば北緯65度付近での夏半年の日射量変動 (夏至の日歳差周期が何らかの理由で近い値になると共鳴が発生

平均日射量の傾向もこれとほぼ似たようなもの) を議し, 惑星の自転軸は非常に大きな振幅で振動するよう論することが多いので, 本稿でもそれに従った. 地球になる[16]. これが自転軌道共鳴という現象であり,に於けるこの時期のこの緯度帯の日射量には赤道傾角現在の火星はまさにこの条件下にあるのである. 火星項よりも気候的歳差項の影響が強く, そのことは周波の日射量変動に於いては気候的歳差項に比べて赤道傾数領域に於ける図2にも良く示されている. ここに角項の寄与が大きく卓越する (図2)のもこの理図を示すことはしないが低緯度域では赤道傾角項の周由による. 前述したように地球の場合には質量の大き期性は更に弱く, 日射量の時間変化はもっぱら気候的な月に伴われているおかげで自転軸の歳差周期が短か歳差項の変動を反映したものとなる.

く, 幸いにして自転軌道共鳴が発生して赤道傾角が大図2に示された地球と火星の日射量変動の様子を変動するようなことはない. けれどももし月が存在し見るとすぐに気が付くことがあろう. 火星の日射量変なかったり, 存在しても現在よりもずっと小さくて軽動の振幅が地球のそれに比べてかなり大きいというこいものだったりしたら, 地球の赤道傾角も火星と同様とである. 大雑把に言って図2に於ける地球への日な大変動を経験する可能性がある. 現在の地球の穏や射量の平均値は約500W/㎡, 最大変動振幅が±約50かな気候状態のひとつの要因が赤道傾角振動の振幅のW/㎡ (平均値に対して±約10％の振幅) であるのに小ささであるとすれば, 大きな月の存在が地球の気候

対して火星のそれは平均値が約250W/㎡, 最大変動振を穏やかなものに保つために重要な役割を果たして来幅が±約100W/㎡であり, 平均値に対して±約40％もたのだと言えるかもしれない[17].

2の振幅を持っている . 地球と火星の平均日射量の違

火星の自転軸歳差運動や自転軌道共鳴については

いは太陽から各惑星までの平均的距離比の逆二乗の反1970年代から W. Ward らによって精力的に研究がな映結果だと思えば良いが (太陽からの平均距離比は地されて来た[13,18,19,20]. 最近では自転軌道共鳴の理球/火星〜1/1.5), 振幅の比は赤道傾角の振幅の違いが論が拡張され, 太陽系外惑星の自転軸安定性と居住可大きく反映されている. 図2にあるように火星の赤能性 (habitability) に関する議論に適用されたりして道傾角は地球のそれに比べて極めて大きな振動を行う.

いる[21].

図2での火星の赤道傾角の最大値はせいぜい36度から37度程度だが, より長い時間スケールではこれより

も遙かに大きな値を取ることがわかっている. 最新の計算結果によると火星の赤道傾角は過去46億年にわた

前節までに述べたように, 地球や火星に到達する日り40度近い平均値を持ちながら激しく振動し, その最射量の分布や時間変化については天体力学的手法によ

大値は80度を超えるという報告までなされている[15].

り精密に計算が可能である. そして地球に於ける±10この原因はいわゆる自転軌道共鳴というものである.

％程度の日射量変動が第四紀の氷期・間氷期サイクル  !"#$% &'()*+,-./
01234 56(--789
:"%

火星の日射量変動と気候／伊藤143を駆動したのであれば, 火星に於ける±40％もの日射測の域を完全には出ていないと言っても過言ではない.

量変動が大きな気候変動を持たらした可能性はかなり

日射量に関しても, 私達が精密に計算できるのは惑星高いのではないかという期待が発生するのは当然であの大気上端に入射する太陽からのエネルギー放射量にる. 火星の気候変動に関する理論的予測や地形データ過ぎない. 地球に関して言えばそこから下には厚い大から推測した気候進化の研究は古くから行われて来た

気があり, 惑星を広く覆う海洋があり, 陸面があり,が, 近年の火星探査衛星による精密な画像データを基

地下水があり, 高緯度域には氷床がある. 大気上端に礎として, こうした議論は精密さを増している. 最近届いた日射量の変動が氷床の変動として発現するまで

では遂に火星の日射量変動と氷期・間氷期サイクルのにはこうした複雑な気候システムの要素をひとつずつ関係に直接言及した論文が出版されるに至っている.

介して行くことになり, その過程は単純なものでは決文献[22]は Mars Global Surveyer (MGC) により撮影してない. 言葉を変えれば, 日射量変動が氷床変動にされた火星北極域にある氷と塵の堆積構造に着目し,直結して反映されるという一般的保証はどこにも無い堆積物の輝度変動と北極域へ入射する夏の日射量変動のである. 実際のところ第四紀に於ける地球の氷床変の関係を比較した. 堆積速度に関する客観的なデータ

動の最も強い周期性は約10万年であるが, この周期はが得られていないためにこの比較は主観的なパターン

日射量変動に於いては非常に弱いパワーしか持ってい認識となるが, 著者はその結果を元にして過去約百万ないという事実は文献[5]の時代から知られて来た大年間の日射量変動がこの堆積構造の要因であると堂々問題であり, これに対しては現在に至るまで未だ決定と主張している. また文献[23]は火星探査機からの様々的な回答は得られていない. 長い研究の歴史と膨大ななデータを改めて総合的に解釈し, 火星の中緯度域ま量の観測データが蓄積されている地球の気候に関してで広がる塵と氷の堆積物が日射量変動に起因するもの

も, 大気上端に届く日射量変動と気候変動 (氷床変動)であり, 過去210万年前から40万年前までは火星に氷との関係はさほど自明ではないのである. 火星の場合期が訪れていたと言う. この論文の著者の説によると,も同様であり, 火星の気候に関する知見もそれにりに

火星の氷期は極域の日射量が大きくなって大気中の水蓄積されてはいるものの (例えば[24]などを見よ),蒸気と塵が中緯度域へ運ばれることにより発生する.

大気上端に入射する日射量の変動が計算できたからと即ち極域が暖かい時期が火星の氷期であり, 地球のそ言って火星の気候変動を理解できたことにはまったくれとは逆だとする説である.

ならない. しかも, 火星の表層気候を司る物理過程はこうした報告を聞くと火星表層で繰り広げられてい地球のそれとは大きく異なる. 地球上には巨大な熱容る気候変動と日射量変動の相関が具体的に手に取れる量を誇る海洋が存在し, 氷期・間氷期サイクルに大きようでわくわくする面もあるが, 現時点での議論を聞な関与をして来たと予想されている. 少なくとも現在

くにあたり私達は慎重な態度を崩してはいけないであの火星には海洋は存在しないから, 気候変動に関するろう. 火星探査機から送られて来る堆積構造写真の解

この違いは極めて大きいであろう. また地球上では大像度がいくら高いとは言え, 私達はそうした堆積物の気中の塵が雨で洗い流され (rainout), その大部分が具体的な形成速度 (時間的な堆積レート) も知らなけ

海底に堆積する現象が発生し得るが, 火星ではそれもれば, 堆積のメカニズム自体も知らない. そのような期待できず, 多量の塵が大気中に放出されて気候に強堆積構造の時間スケールが日射量変動のそれと全く異い影響を与え続けて来た可能性がある. 更に地球ではなっている可能性すらあるわけである. 上で紹介した氷期・間氷期サイクル中にも大気の主要成分が大きくような報告の結果は ｢火星の気候変動が地球と同様に変わるということは全くと言って良いほど無かったと日射量変動に駆動されていてほしい｣ という希望的観思われるが, 火星に於いては大気の主成分である二酸144

日本惑星科学会誌Vol.13.No.3,2004

化炭素が極冠にあるドライアイスと平衡状態にあるために, 条件によっては大気全体の気圧や二酸化炭素濃[11] 伊藤孝士, 2002: ｢全地球史解読｣ , 137, 東京大学出版会.

度が甚大に変化し得る. こうした違いを始めとし, 火[12] Ito, T. et al., 1995: J. Geophys. Res. 100, 15147.

星の大気, 極冠, 陸面, あるいは地下水 (地下氷) な[13] Ward, W.R., 1974: J. Geophys. Res. 79, 3375.

どの複雑な相互作用体系としての気候システムに関し[14] Berger, A.L., 1978: J. Atmos. Sci. 35, 2362.て私達の造詣が地球並みあるいはそれ以上に深くなっ[15] Laskar, J. et al., 2004: Icarus, 170, 343.て初めて, 日射量変動が火星の気候にどのように影響[16] Laskar, J. and Robutel, P., 1993: Nature 361, 608.

して来たかを正確に知ることが出来るようになろう.

[17] Laskar, J. et al., 1993: Nature 361, 615.

そこへ到達する道程はまだまだ長く険しいものである[18] Ward, W.R., 1979: J. Geophys. Res. 84, 237.

ように思われる.

[19] Ward, W.R. et al., 1979: J. Geophys. Res. 84,243.

[20] Bills, B.G., 1990: J. Geophys. Res. 95, 14137.
[21] Atobe, K. et al., 2004: Icarus 168, 223.

本稿の執筆に際しては東京大学気候システム研究セ[22] Laskar, J. et al., 2002: Nature 419, 375.

ンターの阿部彩子さんと東京大学大学院理学系研究科[23] Head, J.W. et al., 2003: Nature 426, 797.の田近英一さんから多数の助言を頂いた. この場を借[24] Nakamura, T. and Tajika, E., 2003: Geophys. Res.

りてお二人に深く御礼を申し上げたい.


[1] 阿部彩子, 2002:｢全地球史解読｣ , 234, 東京大学出版会.
[2] 中島映至, 1980: 気象研究ノート 140, 81.
[3] Milankovitch, M., 1941:     ,:Koniglich Serbische Academie.
[4] 増田耕一, 1993: 気象研究ノート 177, 223.
[5] Hays, J.D. et al., 1976: Science 194, 1121.
[6] 阿部彩子, 増田耕一, 1996: 岩波口座地球惑星科
学11 ｢気候変動論｣ , 103, 岩波書店.
[7] Berger, A.L. et al., 1984:  
,D. Reidel.
[8] 安成哲三, 柏谷健二, 1992: ｢地球環境変動とミラ
ンコヴィッチ・サイクル｣ , 古今書院.
[9] Laskar, J., 1988: Astron. Astrophys. 198, 341.
[10] Ito, T. and Tanikawa, K., 2002: Mon. Not. R.
Astron. Soc. 336, 483.

Lett. 30, 1685, doi:10.1029/2002GL01672.