一方、更新世初期には気候の歳差運動の周期には弱い変動しか見られず（図1a）、氷河期の周期は赤道傾角の周期が支配的であり、これがいわゆる41 kyr問題を構成している。 41 kyr問題を解決する機構として、夏の日射の強さと夏の長さの反相関、北半球（NH）と南半球（SH）の氷床の変動の間にある20 kyr成分のオフセットなどが提案されている。 歳差運動の信号はその短周期性から赤道傾角に比べ検出が困難であるが、最近の研究では、軌道に依存しない年代モデルを持つ古気候代替指標を用いて、前期更新世の歳差運動の変動が確認されている。 また、MPT以前の時代には、サイクルの非対称性（すなわち、速い氷床後退期）の存在が推察されている。 MPT期以前の氷床後退期の決定における歳差の役割は不明であるが、最近、北大西洋の氷筏の記録を用いて、MPT期以前でも気候の歳差の極小時に退氷が発生することが示された。 しかし、41 kyr氷期サイクルの終了時期や間氷期の期間を決定する上で、赤道傾角と歳差の相対的な寄与は定量的には解明されてはいない。
私達は今回、動的な三次元NH氷床モデルIcIESと大循環モデルMIROCの結果に基づく気候パラメタリゼーションを組み合わせたモデル IcIES-MIROC を採用した（補足図1）。 このモデルでは、天文学的強制力の変化によって引き起こされるNH陸地の気候変動を離散的なGCMスナップショットを用いた一連の実験に基づく気候パラメタリゼーションを使って推定する。 氷床モデルには気候強制力として地表温度が与えられる。 この設定により、速いフィードバック（水蒸気、雲、海氷、植生のフィードバック）と遅いフィードバック（アルベド/温度/氷床、経過率/温度/氷床のフィードバック）の両方を表現できる。 これまで、このモデルは過去400 kyrにおける100 kyr氷期サイクルを再現することに成功している。 本研究ではIcIES-MIROCで用いられている気候パラメタリゼーションを改良し、NH高緯度の植生フィードバックを考慮した（補足図1、2、補足表1）。 このモデルを用いて、私達はδ18Oシグネチャー内に41 kyrの氷期サイクルが明確に観測される約1.6-1.2Ma（海洋同位体段階、MIS、53から36まで；表1）の氷期サイクルをシミュレーションした（図1dの青枠の部分）。 この時期には離心率（図1a）、赤道傾角（図1c）、δ18Oシグネチャーの41 kyr成分の振幅が最大となり、大気中の二酸化炭素濃度は後期更新世のそれに匹敵する（図1c）。 このことは、41 kyrの氷期サイクルを駆動する天文的な強制力の役割を調べるためにこの時代が理想的であることを示している。 我々は氷床後退期のタイミング、氷河期の時系列の形状、氷床の形状、41 kyr氷河期の周期性を決定する上での赤道傾角と歳差の役割について議論する。