两万年前的地球是什么样子?加拿大和美国北部被厚达数公里的冰盖压在下,北欧埋在另一块冰盖之下,海平面比今天低约 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 万年的长周期:
其中 $a$、$b$ 为轨道半长轴与半短轴。偏心率越大,地球在近日点与远日点接收的日照差异越悬殊。
黄赤交角(obliquity)——地轴相对于轨道面的倾角,在约 22.1°—24.5° 之间摆动,周期约 4.1 万年。倾角决定季节对比的强弱:倾角越大,两极夏季获得的日照越多、冬季越暗。
岁差(precession)——地轴像陀螺般缓慢摇摆,周期约 1.9 万和 2.3 万年。它决定地球抵达近日点时北半球处于什么季节:如今近日点在 1 月初(北半球冬季),约 1.1 万年后将挪到北半球夏季。岁差与偏心率的乘积(气候岁差参数,常写作 )调制季节日照的振幅。
三者叠加,改变的是瞬时日照
其中 为太阳常数(约 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 重建的迁徙路线沿白令陆桥等间冰期前的陆地通道展开,冰期地理学的预测正被分子证据逐段证实与修正。
参考文献
- Milanković, M. (1941). Kanon der Erdbestrahlung und seine Anwendung auf das Eiszeitenproblem. Royal Serbian Academy, Belgrade. (米兰科维奇周期理论的集大成之作)
- Hays, J. D., Imbrie, J. & Shackleton, N. J. (1976). Variations in the Earth's orbit: pacemaker of the ice ages. Science, 194, 1121–1132. (谱分析证实冰期记录中的三个轨道周期)
- 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 深海氧同位素标准曲线与轨道调谐年表)
- 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 万年温度与温室气体记录)
- EPICA Community Members (2004). Eight glacial cycles from an Antarctic ice core. Nature, 429, 623–628. (EPICA 冰芯把记录扩展到 80 万年、八个冰期旋回)
- 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₂—风化壳机制假说)