跳转到内容
← 返回核心概念
地球系统地球科学 · 古气候学15 分钟阅读

冰期与米兰科维奇周期

Ice Ages and Milankovitch Cycles

两万年前的地球是什么样子?加拿大和美国北部被厚达数公里的冰盖压在下,北欧埋在另一块冰盖之下,海平面比今天低约 120—130 米,白令海峡的位置是一座连接西伯利亚与阿拉斯加的干冷陆桥。而这一切的启动按钮,不在地球上,而在地球绕太阳运行的轨道里。塞尔维亚学者米兰科维奇(Milutin Milanković)用一支笔和几十…

米兰科维奇周期冰期轨道强迫氧同位素末次盛冰期

两万年前的地球是什么样子?加拿大和美国北部被厚达数公里的冰盖压在下,北欧埋在另一块冰盖之下,海平面比今天低约 120—130 米,白令海峡的位置是一座连接西伯利亚与阿拉斯加的干冷陆桥。而这一切的启动按钮,不在地球上,而在地球绕太阳运行的轨道里。塞尔维亚学者米兰科维奇(Milutin Milanković)用一支笔和几十年的计算指出:冰期与间冰期的轮回,听从三个天文参数的指挥。这套理论如今被称为米兰科维奇周期(Milankovitch cycles)。关于冰川如何雕刻地貌、冰期概念的发现史,见 glaciation-ice-ages;本文聚焦轨道驱动本身的机制与证据。

破除误解:冰期的钥匙不在"寒冬",也不在太阳变暗

两个直觉都错了。第一,冰期不是因为"冬天特别冷、雪下得特别大"。极地冬天再冷,降雪量也有限——冷空气能容纳的水汽很少。真正的开关在夏天:只要夏季凉到化不完冬天的积雪,雪就会年年有剩余,积少成多,最终堆成冰盖。冰期的触发词是"凉夏",不是"寒冬"。

第二,驱动冰期的不是太阳本身变亮变暗。太阳输出在万年级尺度上相当稳定;变化的是地球接收阳光的时空分配——阳光在一年里何时到达、落在哪个纬度。米兰科维奇的全部工作,就是精确计算这种分配的节律。

历史脉络:一个在战俘营里开始的计算

用天文解释冰期的想法比米兰科维奇更早:19 世纪就有学者提出轨道变化可能是冰期的幕后推手,苏格兰自学出身的科学家克罗尔(James Croll)在 1875 年把它系统化,但苦于没有可靠的年代证据,始终停留在猜想。

真正把猜想变成定量理论的是米兰科维奇(1879—1958)。他出生于多瑙河畔的达利(今克罗地亚境内),在维也纳取得土木工程博士学位,后任贝尔格莱德大学应用数学教授。第一次世界大战爆发时,他作为塞尔维亚预备役军官在奥匈帝国被拘,拘禁于布达佩斯的几年里,他开始了太阳辐射与地球气候的数学理论。这项研究在 1941 年结出果实——600 余页的巨著《地球日照典则及其在冰期问题上的应用》(Kanon der Erdbestrahlung,通称《日照典则》)在贝尔格莱德出版。命运几乎把它埋没:书稿付印仅数日,1941 年 4 月德军轰炸贝尔格莱德,印刷厂中弹,所幸书版幸存,战时得以重印,但流布极少。

更冷的冷遇来自学界。没有精确的古气候年表,就无法检验理论预测的冰期时间。直到 1976 年,海斯(James Hays)、英布里(John Imbrie)与沙克尔顿(Nicholas Shackleton)在《科学》杂志发表著名论文:他们对深海岩芯的氧同位素记录做谱分析,找到了 10 万年、4.1 万年、2.3 万年三组周期,与轨道参数精确对应。被搁置半个多世纪的《日照典则》,终于被证明是对的。

核心机制:三个天文旋钮

米兰科维奇理论包含三个周期性变化的参数:

偏心率(eccentricity)——地球公转轨道椭圆的扁圆程度,在约 0.005 到 0.06 之间变化(现今约 0.017),主周期约 10 万年,另有约 40 万年的长周期:

e  =  1b2a2e \;=\; \sqrt{1-\frac{b^{2}}{a^{2}}}

其中 $a$$b$ 为轨道半长轴与半短轴。偏心率越大,地球在近日点与远日点接收的日照差异越悬殊。

黄赤交角(obliquity)——地轴相对于轨道面的倾角,在约 22.1°—24.5° 之间摆动,周期约 4.1 万年。倾角决定季节对比的强弱:倾角越大,两极夏季获得的日照越多、冬季越暗。

岁差(precession)——地轴像陀螺般缓慢摇摆,周期约 1.9 万和 2.3 万年。它决定地球抵达近日点时北半球处于什么季节:如今近日点在 1 月初(北半球冬季),约 1.1 万年后将挪到北半球夏季。岁差与偏心率的乘积(气候岁差参数,常写作 esinϖe\sin\varpi)调制季节日照的振幅。

三者叠加,改变的是瞬时日照

F  =  S0(rˉr)2coszF \;=\; S_{0}\left(\frac{\bar{r}}{r}\right)^{2}\cos z

其中 S0S_{0} 为太阳常数(约 1360 瓦/平方米),$r$ 为日地距离,$z$ 为太阳天顶距。关键结论有两条:其一,轨道变化对全球年均总日照的影响微乎其微,它做的是"重新分配"而非"增减总量";其二,真正控制冰盖存亡的是北半球高纬度夏季日照——北半球陆地集中,冰盖只能在那里生长,而 65°N 附近夏季日照的摆动幅度可达约两成。三个旋钮何时同向叠加出"凉夏",冰期便在何时叩门。

证据链:从深海岩芯到冰芯

理论再好,也要档案作证。第一条证据链在海底。浮游与底栖有孔虫的碳酸钙壳体记录着海水的氧同位素(δ¹⁸O):冰期时,大量偏轻的海水被锁进冰盖,海洋整体变"重",有孔虫壳的 δ¹⁸O 随之升高。δ¹⁸O 曲线因此成了全球冰量的代理。1976 年的"节拍器"论文正是靠它验证了轨道周期。1984 年,英布里领导的 SPECMAP 计划把多条深海岩芯的 δ¹⁸O 叠合成标准曲线,并用轨道周期反向校正年代标尺——这套"轨道调谐"(orbital tuning)的年表,至今仍是第四纪研究的时间基准(后经 Martinson 等 1987 年细化,并由 Lisiecki 与 Raymo 2005 年的 LR04 叠加曲线延展到 500 多万年)。深海记录还把气候史切成一个个海洋氧同位素阶段(MIS):奇数为暖、偶数为冷,我们身处 MIS 1,末次盛冰期在 MIS 2。

第二条证据链在冰盖里。冰芯中的气泡直接封存古代大气,冰本身的同位素记录当地温度。1999 年发表的南极东方站(Vostok)冰芯把气候史回溯到 42 万年前;2004 年发表的 EPICA 冰芯更达到 80 万年、覆盖八个冰期旋回。冰芯给出的最震撼图景是:温度与大气 CO₂ 浓度在冰期—间冰期之间严丝合缝地同步起落——冰期 CO₂ 低至约 180—190 ppm,间冰期高至约 280—300 ppm。轨道给节拍,碳循环放大音量,两者合奏出冰期旋回。两类档案还互相校验:深海 δ¹⁸O 提供全球冰量与长时间连续性,冰芯提供大气成分的直接采样,石笋(见本域《岩溶与洞穴》)则提供陆地上的绝对定年——三套独立证据在轨道周期上汇合,是古气候学最坚固的部分。

争议与边界:10 万年问题与中更新世转型

米兰科维奇理论并不完美,它的伤口很清楚。约 80 万年以来,冰期节律由 10 万年的偏心率周期主导,可偏心率对日照的调制是三者中最弱的——这就是著名的"10 万年问题"。更早一段历史同样棘手:大约 120 万—70 万年前,冰期旋回从 4.1 万年节奏切换为 10 万年节奏(中更新世转型),而轨道参数本身并未发生相应变化。这说明地球系统内部——冰盖的厚度与惰性、冰下风化壳的磨蚀、CO₂ 的长期下降——会"整流"和放大微弱的轨道信号。2019 年 Willeit 等人的模拟提出,CO₂ 下降叠加冰盖下松散沉积物被磨光,共同促成了这次转型,但细节仍在争论。

另一条边界是"轨道驱动必要而不充分"。2016 年 Ganopolski 等人提出,冰期启动需要夏季低日照与低 CO₂ 同时满足——日照降到阈值以下而 CO₂ 偏高,冰盖就拒绝生长。轨道给了气候系统一张时刻表,但车开不开,还取决于碳循环与冰盖自身的物理。

末次盛冰期的世界

把时钟拨回末次盛冰期(LGM,约 2.65 万—1.9 万年前):全球平均气温比今天低约 4—6°C,大气 CO₂ 只有约 185 ppm。劳伦泰德冰盖覆盖加拿大与美国北部,斯堪的纳维亚冰盖覆盖北欧,冰量把海平面压低约 120—130 米。退缩的海洋把大陆架晾成陆地:白令陆桥连接起西伯利亚与阿拉斯加,早期人类得以走进美洲;东南亚的巽他大陆架连成广袤的"巽他古陆";澳大利亚与新几内亚则合并为"萨胡尔"大陆。从欧洲草原到西伯利亚,寒冷干燥的"猛犸象草原"供养着成群的猛犸象、野牛与驯鹿,也供养着追踪它们的狩猎者。今天的世界地图,是约 1.17 万年前冰期结束、海平面回涨之后才画成的——我们熟悉的海岸线、岛屿与海峡,全都是间冰期的产物。

为什么今天仍然重要

第一,米兰科维奇周期仍是气候模型的"试金石"。模型若不能复现轨道驱动下的冰期—间冰期温度与冰量变化,就没有资格预测未来。

第二,它回答了"下一次冰期何时来"。2002 年 Berger 与 Loutre 计算发现,由于当前偏心率正处于低位,即使没有人类干扰,本次间冰期也可能异常漫长,自然状态下或可持续约 5 万年。2016 年 Ganopolski 等人进一步指出,按自然条件下一次冰期启动约在 5 万年后,而人为排放的 CO₂ 很可能把它再推迟数万年。与此同时,末次间冰期(MIS 5e,约 12.9 万—11.6 万年前)比工业革命前暖约 1—2°C、海平面高约 6—9 米,它提醒我们:即便只把今天的变暖维持住,冰盖的长期响应也足以重塑海岸线。

第三,也是最尖锐的对照:过去 80 万年里,大气 CO₂ 始终在约 180—300 ppm 的自然区间内随冰期节拍起伏;而今天已超过 420 ppm,完全越出了米兰科维奇节律所熟悉的范围。轨道节拍器仍在按原速运转,但气候系统已被注入了一种它从未在这样的速度下经历过的强迫。理解冰期如何开始与结束,正是理解我们正在进行的这场无先例实验的基线。

跨域连接

  • 傅里叶分析:米兰科维奇理论的平反是一场谱分析的胜利——把深海 δ¹⁸O 序列分解到频率域,10 万、4.1 万、2.3 万年三组谱峰与轨道周期一一对应且相位吻合,这种多重匹配极难由巧合伪造,是天文驱动最硬的证据形态。
  • 统计学:"轨道调谐"本质上是带物理先验的统计定年——在沉积速率的约束下把同位素曲线与轨道目标曲线做最优对齐;它的力量与风险同源于此:年表精度依赖先验,独立测年(如石笋铀钍年龄)因此成为检验调谐偏差的关键外部裁判。
  • 米卢廷·米兰科维奇:理论的命运与提出者的命运交织——一战拘禁中起步的计算、1941 年战火中几乎被炸毁的巨著、身后近二十年才到来的证实,是"超前于证据的理论如何存活"的经典案例,也说明科学的接受从来不只取决于对错,还取决于证据链何时补齐。
  • 农业革命:末次冰期终结、全新世开启了一个万余年罕见的暖稳窗口,农业、城市与文字才相继出现——文明的全部时间深度都嵌在一个间冰期里,气候稳定期是文明的隐藏地基。
  • 古 DNA 革命:冰期的海平面与冰盖进退写进了人类基因组——古 DNA 重建的迁徙路线沿白令陆桥等间冰期前的陆地通道展开,冰期地理学的预测正被分子证据逐段证实与修正。

参考文献

  1. Milanković, M. (1941). Kanon der Erdbestrahlung und seine Anwendung auf das Eiszeitenproblem. Royal Serbian Academy, Belgrade. (米兰科维奇周期理论的集大成之作)
  2. Hays, J. D., Imbrie, J. & Shackleton, N. J. (1976). Variations in the Earth's orbit: pacemaker of the ice ages. Science, 194, 1121–1132. (谱分析证实冰期记录中的三个轨道周期)
  3. Imbrie, J., Hays, J. D., Martinson, D. G., et al. (1984). The orbital theory of Pleistocene climate: support from a revised chronology of the marine δ¹⁸O record. In Milankovitch and Climate (Berger et al., eds.), Reidel, 269–305. (SPECMAP 深海氧同位素标准曲线与轨道调谐年表)
  4. Petit, J. R., Jouzel, J., Raynaud, D., et al. (1999). Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica. Nature, 399, 429–436. (东方站冰芯 42 万年温度与温室气体记录)
  5. EPICA Community Members (2004). Eight glacial cycles from an Antarctic ice core. Nature, 429, 623–628. (EPICA 冰芯把记录扩展到 80 万年、八个冰期旋回)
  6. Clark, P. U., Dyke, A. S., Shakun, J. D., et al. (2009). The Last Glacial Maximum. Science, 325, 710–714. (末次盛冰期的时间界定、冰盖范围与海平面)

延伸阅读

  • Imbrie, J. & Imbrie, K. P. (1979). Ice Ages: Solving the Mystery. Harvard University Press. (冰期理论百年争论的经典科学史)
  • Berger, A. & Loutre, M.-F. (2002). An exceptionally long interglacial ahead? Science, 297, 1287–1288. (自然情景下本次间冰期的超长寿命)
  • Ganopolski, A., Winkelmann, R. & Schellnhuber, H. J. (2016). Critical insolation–CO₂ relation for diagnosing past and future glacial inception. Nature, 529, 200–203. (冰期启动的日照—CO₂ 双阈值与下一次冰期的推迟)
  • Lisiecki, L. E. & Raymo, M. E. (2005). A Pliocene-Pleistocene stack of 57 globally distributed benthic δ¹⁸O records. Paleoceanography, 20, PA1003. (LR04 叠加曲线:500 多万年全球冰量基准)
  • Willeit, M., Ganopolski, A., Calov, R. & Brovkin, V. (2019). Mid-Pleistocene transition in glacial cycles explained by declining CO₂ and regolith removal. Science Advances, 5, eaav7337. (中更新世转型的 CO₂—风化壳机制假说)