整整100年前，米兰科维奇提出地球轨道参数变化，包括偏心率、倾角以及岁差，影响了地球表面获得的太阳辐射纬度和季节分配，驱动了地球气候准周期变化(Milankovitch, 1920)。过去50年中，米兰科维奇假说不但成了一种古气候研究的理论范式，而且也为古气候研究提供了天文定年方法。
从1980年代中期开始，古气候****利用深海沉积物的碳、氧同位素记录来集成重建新生代气候历史。迄今最有影响的集成曲线是2001年 UCSC古海洋学教授Jim Zachos团队发表在Science上的论文(Zachos et al., 2001a)。这篇文章很大程度上奠定了学界对新生代气候演化的整体认识。但当时，天文定年时间标尺刚刚迈进了古近纪门槛(Zachos et al., 2001b)，新生代早期的时间标尺主要依赖于生物地层和磁性地层，还有不少分歧。
在过去20年里，全球新获得的深海沉积钻孔极大地弥补了这些缺憾（图1），天文时间标尺逐渐跨过新生代，向古生代延伸。新的集成曲线呼之欲出。
图1 CENOGRID钻孔分布（Westerhold et al., 2020附件）
最近，这项工作由德国不来梅大学海洋环境科学中心(MARUM)的Thomas Westerhold领衔完成，论文9月11日在Science发表(Westerhold et al., 2020)。
他们在超过1000个深海沉积钻孔中挑选出14个钻孔，仔细检查并修正了这些岩芯的拼接方式，选择两个长寿的有孔虫属Cibicidoides和Nuttallides的氧、碳同位素记录重建气候历史。根据初步的时间标尺，补充测试了部分晚中新世到早始新世的样品，来保证足够的时间分辨率。他们收集了所有记录已有的天文时间标尺，并且把这些时间标尺统一调整到La2010b的轨道方案，最终获得了一条经过天文调谐定年的，连续覆盖整个新生代全球气候参考曲线CENOGRID (CENOzoic Global Reference benthic foraminifer carbon and oxygen Isotope Dataset)（图2）。
图2 CENOGRID数据集（Westerhold et al., 2020）
数据集由23629个数据点组成，时间分辨率渐新世以来高达2ka，古新世和始新世为4.4ka，估计年代误差古新世-始新世为10万年，渐新世-中中新世为5万年，晚中新世-更新世为1万年。这是全球第一条完整覆盖新生代的高清晰度同位素地层参考曲线。
全球气候是一个复杂的动态系统，在万年到百万年尺度上，对准周期性的天文强迫有复杂的非线性响应。为研究CENOGRID的时域特征，研究团队进行了重现分析(recurrence analysis, 也称递归分析)。重现分析可以揭示系统的非线性动力学过程，以及非线性相互作用信息，重现图是对时间序列的内部结构及可预测性的可视化。
CENOGRID氧同位素的重现图揭示了4种截然不同的方块区域（图3）。每个方块对应着气候在特定状态下反复循环。新生代气候据此可以分成4个状态，称为热室、温室、冷室和冰室状态。热室状态是从56 Ma持续到47 Ma，温度比现在高10度，还出现了氧同位素和碳同位素同步负漂的极热事件。温室状态有两个时期，分别是66 Ma到56 Ma，以及47 Ma到34 Ma。这两个温室状态温度相似，但碳同位素值和二氧化碳浓度截然不同。34 Ma前后始新世-渐新世界线是新生代最显著的转换事件。冷室状态从34 Ma持续到3.3 Ma，以13.9 Ma可以分为两个阶段。3.3 Ma之后为冰室状态，地球的气候受北半球冰盖消长控制。这四种状态的划分，与之前对新生代气候历史的认识大体相符，但是重现图首次提供了统计稳健的客观证据。
图3 CENOGRID氧同位素重现图清晰显示了新生代气候的不同状态（Westerhold et al., 2020）
演化谱显示（图4），在13.9 Ma以前，包括地球热室、温室以及冷室状态的第一阶段，气候周期以偏心率为主，意味着低纬过程驱动着气候演化，推测偏心率调控岁差，影响了季风降水的季节性，进而作用于全球水汽和能量的分布。而倾角周期在暖室热室和第一阶段的冷室状态中表现不显著，可能和当时冰盖的规模小，缺乏高纬的放大效应所致。随着高纬变冷和冰盖增长，13.9 Ma之后地轴倾斜度的信号逐渐增强，到3.3 Ma成为冰室地球气候系统的主导周期。
图4 CENOGRID的演化谱（改自Westerhold et al., 2020）
作者们还计算了CENOGRID曲线重现分析的确定性（图5），定量描述系统的可预测性。当确定性接近0，表示系统是随机的，不可预测，接近1，表明系统确定。结果显示，温室和热室地球比冷室和冰室地球更可预测。34 Ma南极冰盖出现，确定性参数显著降低，地球气候系统非线性程度大大增强。在冷室地球第一阶段中，南极冰量相对较少的25-14 Ma，确定性也相对较高。碳、氧同位素比较而言，6 Ma之前，北极冰量增加，碳同位素的确定性高于氧同位素，主要原因可能是碳同位素主要受制于低纬过程，受极地冰量影响较小。直到6 Ma之后，冰冻圈才深刻地影响了碳循环。到3.3 Ma之后，氧同位素显示出强烈的偏心率周期，确定性有所增强。
图5 CENOGRID重现分析的确定性（改自Westerhold et al., 2020）
相对而言，热室地球比温室地球更加不易预测，主要原因是极热事件显示了强烈的非线性过程，放大了天文强迫。另外，值得注意的是，47 Ma之后，确定性波动幅度越来越大，直到34 Ma到达临界点，变成了不可预测的状态。据此，作者推测，两极的冰量不仅仅定义了地球的基本气候状态，同时影响了气候系统对天文辐射响应的可预测性。
CENOGRID团队还同时公开了对未来气候的预测：目前人为造成的全球变暖的速度远远超过了在新生代任何时候的自然气候波动，并且有可能将地球气候从目前的冰室推向热室状态（图6）。“预计的人为变暖将比这要大得多，IPCC预测，如果‘一切照旧’，2300年全球气温达到5000万年以来的最高水平。”
图6 CENOGRID团队对未来气候的预测（图片来源：Thomas Westerhold）
主要参考文献
Milankovitch M. Théorie Mathématique des Phénomènes Thermiques produits par la Radiation Solaire [M]. Paris: Gauthier-Villrs, 1920.（链接）
Westerhold T, Marwan N, Drury A J, et al. An astronomically dated record of Earth’s climate and its predictability over the last 66 million years[J]. Science, 2020, 369(6509): 1383-1387.（链接）
Zachos J, Pagani M, Sloan L, et al. Trends, rhythms, and aberrations in global climate 65 Ma to present[J]. Science, 2001a, 292(5517): 686-693.（链接）
Zachos J C, Shackleton N J, Revenaugh J S, et al. Climate response to orbital forcing across the Oligocene-Miocene boundary[J]. Science, 2001b, 292(5515): 274-278.（链接）
（撰稿：唐自华/新生代室）
