古气候学方法论：重建过去万年气候的系统方法

第 1 章：古气候学方法论：重建过去万年气候的系统方法

引言：从自然档案中解读气候密码

古气候学（Paleoclimatology）作为一门重建地质历史时期气候状态的学科，其核心挑战在于如何从已然改变的自然记录中提取可靠的气候信息。地球的气候系统在漫长的地质演化过程中留下了形形色色的印记——从极地冰盖中封存的古老空气泡，到深海沉积物中累积的微体化石壳体；从石笋中层层叠加的矿物晶体，到树木年轮中记录的光谱信息。这些自然档案（natural archives）如同一部部无言的编年史，等待着科学方法论的钥匙去开启它们所承载的气候密码。

过去万年的气候记录具有特殊的科学意义。这一时间尺度涵盖了从末次冰期（Last Glacial Maximum）向全新世（Holocene）的转变过程，记录了全新世暖期（Holocene Climatic Optimum）、中世纪暖期（Medieval Climate Anomaly）和小冰期（Little Ice Age）等重大气候波动，同时也承载着人类世（Anthropocene）到来前夕的背景状态。理解这一时期的气候变化机制，不仅对于认识自然气候变率（natural climate variability）至关重要，更为评估当代人类活动对气候系统的影响提供了不可或缺的历史参照。

然而，从自然档案到定量化的气候参数，这一转化过程涉及复杂的物理、化学和生物学原理，以及一系列精密的实验技术与分析方法。本章将系统阐述古气候学的方法论框架，涵盖气候代用指标（climate proxies）的类型与原理、绝对与相对年代学方法、高分辨率年代标尺的建立、数据同化与气候再分析技术，以及不确定性的量化与传播分析。通过这一方法论铺垫，读者将获得理解后续各章节证据解读所需的基础知识与技术储备。

1.1 气候代用指标：自然档案的物理化学基础

气候代用指标（climate proxies）是指那些能够间接反映过去气候状况的物理、化学或生物指标。由于我们无法直接观测地质历史时期的气候要素，必须借助这些“代理”来重建过去的气候状态。每一类代用指标都基于特定的科学原理，将气候信息编码于可测量的物质属性之中。理解这些原理是正确解读古气候记录的前提。

1.1.1 冰芯记录：极地与高山气候档案库

冰芯（ice cores）堪称地球上最为连续、最有时间分辨率的自然档案之一。极地冰盖和高山冰川在积累过程中，将大气中的气体、尘埃、化学物质以及周围环境的温度信息一同封存于冰层之中。位于格陵兰和南极洲的冰芯记录可追溯至距今80万年以上，而中国青藏高原和秘鲁安第斯山等高山冰芯则记录了区域气候的详细信息。

冰芯分析所依据的基本原理在于：冰层本身记录了积累时的温度信息，而这一信息的载体是水的稳定同位素比值。当水分子在蒸发、输送和凝结过程中经历相变时，轻水（H₂¹⁶O）和重水（H₂¹⁸O、H₂D）的行为存在显著差异。重水分子具有较高的分子量，蒸气压较低，因此在蒸发过程中倾向于留在液相，而在凝结过程中则更容易进入固相。这意味着高纬度或高海拔地区降水中的 δD 和 δ¹⁸O 值相对于标准平均海水（Standard Mean Ocean Water, SMOW）偏负。大量观测数据表明，现代大气降水中的稳定同位素比值与当地气温存在显著的线性相关关系，这一关系构成了利用冰芯 δ 值重建古温度的理论基础。

然而，冰芯同位素记录并非温度的直接量度，而是受多种因素综合影响的复杂信号。北大西洋地区冰芯记录的 δ¹⁸O 值与温度的相关系数可达0.7至0.9，而在南极地区这一关系则相对较弱。此外，冰芯中保存的大气气泡提供了直接分析古大气成分的可能，包括二氧化碳（CO₂）、甲烷（CH₄）和一氧化二氮（N₂O）等温室气体浓度。2005年南极EPICA冰芯项目获得的80万年记录显示，大气CO₂浓度与全球温度在冰期-间冰期旋回中呈现高度一致的协同变化，CO₂浓度在冰期约为180 ppmv，而在间冰期可达280 ppmv。

冰芯中的化学记录同样承载着丰富的气候信息。钠（Na⁺）和氯（Cl⁻）离子主要来源于海盐气溶胶，其浓度变化反映了源区海冰范围和大气环流模式的变化；钙（Ca²⁺）离子则与亚洲沙尘输送密切相关，青藏高原冰芯中的钙记录揭示了东亚冬季风在冰期的显著增强。火山喷发产生的硫酸盐（SO₄²⁻）信号可以在冰芯中清晰识别，形成精确的年代标记层（volcanic ash layers），如1991年皮纳图博火山喷发和1993年胡安德夫卡火山喷发都在最近的冰芯记录中留下了可辨认的信号。

1.1.2 树轮年代学：陆地气候的年分辨率记录

树木年轮（tree rings）是另一种高精度、高分辨率的自然档案。树木在生长过程中形成层状结构，每年一轮，这一特性使得精确的年表建立成为可能。树轮年代学（dendrochronology）的科学基础在于：树木生长受环境因子（主要是温度和降水）控制，不同年份形成的年轮在宽度、密度和细胞结构上存在差异，这些差异被忠实地记录于年轮之中。

树轮宽度（tree-ring width）是最直观的生长指标。在树木的生理极限区域，如高山林线或干旱区，年轮宽度主要受夏季温度控制；而在半干旱区，年轮宽度则更多地反映水分供应状况。现代树轮气候学通过建立树木生长与气候要素之间的统计模型，可以定量重建过去的气温和降水变化。中国祁连山地区的祁连圆柏（Juniperus przewalskii）树轮序列已延伸至距今5000年以上，其中距今3585年至公元前150年的序列与《竹书纪年》等历史文献记录高度吻合。

树轮密度（tree-ring density）是比宽度更为敏感的指标。早材（early wood）形成于生长季早期，细胞壁薄、密度低；晚材（late wood）形成于生长季后期，细胞壁厚、密度高。最大晚材密度（maximum latewood density, MXD）是欧洲树轮气候学中应用最广泛的指标之一，芬兰、瑞士等国的MXD序列已成功重建了北半球过去千年的夏季温度变化。欧洲7国联合开展的“千年树轮气候学计划”（PAGES 2k Network）综合了196个树轮序列，重建的北半球温度曲线显示出中世纪暖期（950-1250 CE）和小冰期（1400-1850 CE）的清晰特征。

除宽度和密度外，树轮中还存在更为精细的化学信号。树轮纤维素（cellulose）中的 δ¹³C 值反映了大气 CO₂ 浓度和植物水分利用效率的变化，δ¹⁸O 值则与大气湿度和降水氧同位素组成相关联。近年来，高分辨率的激光剥蚀-质谱技术（laser ablation MC-ICP-MS）使得在单个年轮内部获取多个同位素数据点成为可能，将时间分辨率提升至季节甚至月度尺度。

树轮年代学的优势在于其精确的定年能力和季节级的时间分辨率，但其局限同样明显：单棵树木的记录长度有限，必须通过交叉定年（cross-dating）技术将不同树木的序列拼接为长序列；且树轮记录仅保存于有树木生长的区域，覆盖范围受到生态条件的严格限制。

1.1.3 海洋与湖泊沉积物：长时间尺度的连续记录

海洋沉积物（marine sediments）和湖泊沉积物（lacustrine sediments）构成了百万年级别时间尺度气候记录的主体。这些沉积物连续堆积，形成近乎不间断的档案序列。通过对沉积物岩芯（sediment cores）的分析，可以获得关于古气候、古海洋和古生态的综合性信息。

深海沉积物中的有孔虫（foraminifera）壳体是最重要的古气候代用指标之一。浮游有孔虫在生长过程中，其方解石壳体中的 Mg/Ca 比值与海水温度呈正相关关系，这是因为 Mg²⁺ 取代 Ca²⁺ 进入方解石晶格是一个吸热过程，温度升高有利于 Mg 的掺入。现代实验室的 Mg/Ca 测定精度可达 ±1%，对应于约 ±0.3°C 的温度不确定性。南海 ODP 1143 站位的 Mg/Ca 记录成功重建了过去80万年的表层海水温度变化，揭示了冰期-间冰期旋回中高达5-6°C的温度振幅。

有孔虫壳体的 δ¹⁸O 值同样蕴含丰富信息。壳体 δ¹⁸O 受两方面因素控制：一是生长水体的温度，二是水体本身的 δ¹⁸O 值（后者又与全球冰量变化相关）。深海有孔虫壳体的 δ¹⁸O 记录包含了这两个信号的叠加，需要通过数学分解才能分离出温度信号和冰量信号。古气候学家发展了多种分离方法，其中基于 β 分异因子的热力学计算和基于全球冰量模型的反演是两种常用策略。

沉积物中的陆源组分同样记录了气候信息。粒度分析可以指示河流输入和风力搬运的变化，粗粒层（coarse layers）通常对应于洪水事件或强风天气。沉积物颜色（sediment color）的变化反映了铁的价态和含量，可以指示沉积环境的氧化还原条件。中国黄土高原的洛川剖面和西峰剖面记录了250万年的风成沉积历史，黄土-古土壤序列与深海氧同位素记录高度可比，构成了陆相气候记录的经典范例。

湖泊沉积物具有区域代表性强的优势，特别适合重建局地或区域气候。湖泊沉积速率通常高于海洋环境，可以保存更精细的气候信号。纹层（laminae）发育良好的湖泊沉积物可以提供年分辨率甚至季节分辨率的记录。云南洱海、青海湖等中国湖泊的沉积记录揭示了东亚季风在末次冰期和全新世的演化历史。湖泊沉积物中的色素（pigments）、硅藻（diatoms）和介形虫（ostracods）等生物指标，以及总有机碳（TOC）、总氮（TN）等地球化学指标，共同构成了多指标相互验证的证据体系。

1.1.4 洞穴沉积物：石笋中的稳定同位素档案

洞穴沉积物（speleothems），尤其是石笋（stalagmites），是大陆地区高分辨率古气候记录的重要来源。洞穴碳酸钙的生长源于地表降水通过土壤和岩溶系统向洞穴的渗透。当富含 Ca²⁺ 和 CO₃²⁻ 的水溶液在洞穴中释放 CO₂ 并蒸发时，CaCO₃ 沉淀析出，形成层层叠加的碳酸钙沉积。每一层石笋记录了降水渗入时的同位素组成和地球化学信息。

石笋的 δ¹⁸O 值是应用最广泛的指标。洞穴滴水的 δ¹⁸O 值主要继承自大气降水的 δ¹⁸O，而后者受降水时的温度、源区水汽同位素组成和降水量效应（amount effect）控制。在东亚季风区，夏季风从海洋带来大量水汽，降水 δ¹⁸O 值偏负；冬季风控制时期降水稀少且 δ¹⁸O 值偏正。因此，石笋 δ¹⁸O 记录的高频变化可以反映季风强度的波动。南京葫芦洞石笋记录的 δ¹⁸O 序列与格陵兰冰芯记录的 δ¹⁸O 序列在千年尺度上高度一致，为季风-高纬度气候遥相关提供了有力证据。

石笋中的微层（microlayers）或微层带（microbanding）可以提供季节到年分辨率的时间标尺。在某些石笋中，荧光微层（fluorescent microlayers）的出现具有年周期性，可能与植被生长季的降水事件相关。通过纹层计数建立的石笋年代标尺，其精度可达 ±0.5% 以内，优于多数放射性同位素定年方法。

石笋的 U/Th 定年是所有古气候载体中定年精度最高的，将在下一节详述。U/Th 定年可以给出石笋生长段的绝对年龄，结合精确的纹层计数，便可建立绝对年代控制的石笋年表。中国石笋氧同位素研究在重建亚洲季风演化历史方面做出了突出贡献，贵州董哥洞和湖北三宝洞的石笋记录成为东亚古气候研究的基准序列。

1.1.5 花粉分析：古植被与古气候的桥梁

孢粉学（palynology）通过分析沉积物中保存的植物孢子（spores）和花粉（pollen）来重建古植被和古气候。花粉具有坚硬的外壁（exine），由孢粉素（spoporopollenin）构成，能够在沉积物中保存数百万年。不同植物类群的花粉形态具有种的特异性，通过形态鉴定可以确定其植物来源，进而推断当时的植被组成。

孢粉-气候定量重建的方法论经历了从定性到定量的发展过程。早期研究通过孢粉组合的变化定性描述植被演替和气候波动。现代孢粉学则发展了多种定量重建方法，其中转移函数（transfer function）是最成熟的技术路线。转移函数基于现代孢粉-气候数据集，建立两者之间的数学关系，然后应用于化石孢粉数据以获得古气候参数估计值。反距离加权（Weighted Averaging Partial Least Squares, WA-PLS）回归是目前应用最广的转移函数方法，其原理是通过加权平均和偏最小二乘回归，将孢粉组合转换为气候参数。

中国学者基于全国范围的现代孢粉数据集，发展了多种气候要素的转移函数，包括年平均温度、最热月平均温度、最冷月平均温度、年降水量等。泥河湾盆地和长白山的孢粉记录显示，全新世大暖期（Megathermal）期间研究区域的年均温度比现代高约2-3°C，降水也有显著增加。然而，孢粉-气候关系并非一对一对应，同一气候条件下可以存在不同的植被类型，因此孢粉定量重建具有不可避免的不确定性。

1.2 年代学方法：建立时间框架的基石

年代学（chronology）是古气候学的方法论基石。没有可靠的年代控制，任何气候重建都将失去时间参照，不同地点的记录也无法进行对比和综合分析。古气候年代学采用绝对定年（absolute dating）和相对定年（relative dating）两大类方法，前者给出以年为单位的绝对年龄，后者则确定地层序列的相对新老关系。

1.2.1 放射性碳定年：全新世定年的金标准

放射性碳定年（radiocarbon dating）由美国化学家威拉德·利比（Willard Libby）于1949年发明，开创了第四纪年代学的新纪元。利比因此荣获1960年诺贝尔化学奖。放射性碳定年的原理在于：大气上层通过宇宙射线与氮原子的核反应产生放射性碳同位素 ¹⁴C，后者氧化为 CO₂ 并参与全球碳循环，被所有活着的有机体吸收。当生物死亡后，其体内的 ¹⁴C 不再得到补充，开始以指数衰减：$N(t) = N_0 e^{-\lambda t}$，其中 λ = ln(2)/5730 年 = 0.000120968 年⁻¹。通过测量样品中残留的 ¹⁴C 含量，可以计算出生物死亡至今的时间。

然而，¹⁴C 年龄并非日历年龄，两者之间存在系统偏差。这一偏差的来源有二：一是大气 ¹⁴C 产生率并非恒定，而是受太阳活动和地磁场强度变化影响；二是碳循环过程中存在同位素分馏效应，需要对分馏进行校正。自1986年国际校正曲线（IntCal）首次发布以来，¹⁴C 校正已成为放射性碳定年的标准程序。当前使用的 IntCal20 校正曲线基于树轮、珊瑚和沉积物记录的综合分析，覆盖了全部 14,400 年的可校正范围。

实际应用中，放射性碳定年面临多重挑战。首先是样品污染问题：古老样品可能吸附年轻碳或损失老碳，导致定年结果偏离真值。其次是“半衰期”问题：利比最初使用的 5568 年半衰期与现在公认的 5730 年存在差异，历史文献中的 ¹⁴C 年龄可能基于旧半衰期，需要转换。第三是海洋储库效应（marine reservoir effect）：海洋与大气之间的碳交换存在时间滞后，导致海洋生物的 ¹⁴C 年龄可能比同步的大气年龄老数百年至数千年，这一效应在沿海地区尤其显著。

为获得准确的日历年龄，古气候学家通常采用以下策略：选择保存良好、不易受污染的样品（如植物种子、炭屑、泥炭等陆生材料）；对同一地层的多样品进行定年以评估一致性；利用贝叶斯统计方法整合多个 ¹⁴C 年龄和地层约束以优化年代模型。近年来，加速器质谱（AMS）技术的应用使得所需样品量大幅减少，定年精度显著提高，一块数毫克的植物种子即可获得可靠的 ¹⁴C 年龄。

1.2.2 铀系定年：超越万年的精确计时器

铀系定年（uranium-series dating）基于铀同位素的放射性衰变链。自然界存在两个独立的铀衰变系：铀-238（²³⁸U）衰变系和铀-235（²³⁵U）衰变系。²³⁸U 的半衰期为 4.468 × 10⁹ 年，最终衰变为稳定的铅-206（²⁰⁶Pb）；²³⁵U 的半衰期为 7.04 × 10⁸ 年，最终衰变为 ²⁰⁷Pb。铀系定年利用这两个衰变系中间产物的活度比值计算样品年龄。

对于石笋、珊瑚等碳酸盐样品，²³⁰Th/²³⁴U 法是最常用的铀系定年方法。样品初始只含铀不含钍，因此可以用 ²³⁰Th/²³⁴U 比值计算年龄：

$$t = \frac{1}{\lambda_{230}} \ln\left(1 + \frac{^{230}\mathrm{Th}}{^{234}\mathrm{U}} \times \frac{\lambda_{230}}{\lambda_{234}}\right)$$

其中 λ₂₃₀ 和 λ₂₃₄ 分别是 ²³⁰Th 和 ²³⁴U 的衰变常数。

铀系定年的优势在于其高精度和高确定性。230Th 的半衰期约为 75,700 年，适用于从数百年到约 50 万年年龄范围的样品测定。对于年轻的石笋样品（<5 万年），定年精度可达 ±1% 以内，即 1000 年的样品误差不超过 10 年。这一精度远超放射性碳定年，特别适合建立高精度年代标尺。

然而，铀系定年面临两个主要挑战。其一是初始钍污染问题：大多数铀系定年样品为碳酸盐，其初始钍含量极低，但如果样品吸附了初始钍或碎屑钍（detrital thorium），就会导致年龄高估。现代质谱技术（MC-ICP-MS）通过同时测量 ²³⁰Th、²³²Th 和 ⁸⁷Sr/⁸⁶Sr 比值，可以有效区分继承钍和碎屑钍，显著提高了定年准确性。其二是铀系不平衡问题：铀系定年假设样品形成时处于铀系平衡状态，即 ²³⁴U/²³⁸U 比值等于平衡值。如果样品形成时铀系处于不平衡状态（如某些珊瑚样品），则需要测量 ²³⁴U/²³⁸U 比值并引入额外的校正。

中国石笋古气候学的快速发展得益于铀系定年技术的进步。南京师范大学和美国明尼苏达大学等实验室建立了多个高精度石笋年代标尺，为东亚季风演化研究奠定了坚实的时间框架。2012年发表的“石笋记录的更新世气候变化”综合研究，整合了全球 63 个石笋记录，首次提供了热带和亚热带地区高分辨率的冰期-间冰期气候记录。

1.2.3 光释光定年：石英与长石中的地质时钟

光释光定年（optically stimulated luminescence dating, OSL）是第四纪沉积物定年的核心技术之一。其原理基于矿物的光释光特性：石英和长石晶体在地质历史中累积了辐射损伤，当这些矿物被埋藏后，其晶格中的电子陷阱被持续填充。用光或热刺激样品时，累积的电子从陷阱中释放，以光子形式释放储存的能量，产生释光信号。释光强度与矿物最后一次曝光（ exposure to sunlight）以来累积的辐射剂量成正比。

对于风成沉积（如沙漠砂、黄土）和水成沉积（如河流、湖泊沉积），光释光定年可以给出沉积物的埋藏年龄。这一技术特别适合缺乏有机碳的沉积物，因为这些样品无法用放射性碳定年。青藏高原东北缘黄土沉积物的光释光定年显示，黄土的堆积从约 130 万年前开始，与深海氧同位素阶段（MIS）16 以后的冰期旋回同步。

光