冰川期火山活动如何触发大西洋经向翻转环流(AMOC)的崩溃与重建

本研究通过结合冰芯硫酸盐记录与气候模型模拟，揭示了大型赤道火山喷发如何通过大气-海洋环流变化和海气浮力通量，引发大西洋经向翻转环流(AMOC)的千年尺度突变，为理解冰期气候突变机制及现代气候 tipping points 风险提供了新视角。

在地球漫长的气候历史中，末次冰期（约11.5万至1.17万年前）是一个充满剧烈波动的时代。格陵兰冰芯记录揭示了一系列惊人的快速变暖事件，即Dansgaard-Oeschger（D-O）事件，其间格陵兰地区气温可在短短几十年内飙升10-15°C。这些反复出现的、持续数百至数千年的气候突变，被认为是晚第四纪冰期-间冰期旋回中最显著的特征之一。科学家们普遍认为，D-O事件与北大西洋经向翻转环流（AMOC）在强弱状态之间的切换密切相关。AMOC作为全球海洋输送带的关键组成部分，其变化能重新分配巨大的热量，从而引发全球范围的气候重组。然而，一个核心谜团始终悬而未决：是什么触发了AMOC这些突如其来的崩溃与恢复？尽管轨道强迫、温室气体浓度变化和冰盖动力学等因素都被认为可能参与其中，但火山活动作为潜在触发器的角色长期以来却备受争议，且缺乏基于物理过程的模型证据。

传统观点认为，大型火山喷发主要影响年际到十年尺度的气候，例如导致全球短暂降温。但对于火山活动能否驱动千年尺度的气候变迁，科学界一直存在疑问。主要的挑战在于三个方面：过去气候和火山记录定年的准确性、可用数据集中对历史喷发的低估、以及连接火山强迫与长期气候响应的过程性理解的缺乏。尽管改进的测年技术和更高分辨率的冰芯记录为解决前两个问题带来了希望，但第三个问题——火山活动如何通过复杂的海洋-大气耦合过程最终导致AMOC的突变——仍需借助复杂的气候模型来深入探索。正是在这一背景下，发表在《科学进展》（Science Advances）上的这项研究应运而生，旨在填补这一关键空白。

为了回答火山活动是否以及如何触发冰期气候突变这一问题，研究人员开展了一项雄心勃勃的模型研究。他们首先利用来自格陵兰和南极的冰芯记录，建立了全新世和海洋同位素3阶段（MIS 3，约5.9万至2.7万年前，D-O事件最活跃的时期）的火山硫酸盐沉降序列。通过广义极值分布（GEV）分析，他们估算了不同规模火山事件（如1991年皮纳图博、1815年坦博拉和1257年萨马拉斯火山喷发）的发生频率。在此基础上，研究团队将基于冰芯记录重建的“真实”火山强迫引入到美国国家大气研究中心（NCAR）的社区气候系统模型第4版（CCSM4）中，进行了大量冰期条件下的海气耦合模型模拟。他们设计了一系列精巧的数值实验，包括在数千年长度的控制模拟中，定期施加不同强度的赤道火山气溶胶强迫（如230 Tg、115 Tg和58 Tg，分别对应萨马拉斯、坦博拉和皮纳图博级别的喷发），以观察D-O循环的响应。更重要的是，为了评估气候系统内部变率（即“噪声”）与火山扰动之间的相互作用，研究者在D-O振荡的不同阶段（如间冰阶的早、中、晚期和冰阶阶段）启动了大型集合模拟。每个阶段都包含两组各50个成员的集合：一组施加单次大型火山爆发强迫，另一组则不施加，仅通过微小的初始条件扰动来表征内部变率。通过比较这两组集合中AMOC状态演变的差异，研究人员得以量化火山活动在特定气候背景下触发状态转换的概率。此外，他们还借助一个简化的随机动力系统模型，对复杂模型模拟出的AMOC行为进行了理论阐释，深化了对潜在动力机制的理解。

主要技术方法

本研究主要依赖于气候模型模拟和古气候数据同化分析。关键技术包括：1) 利用社区气候系统模型第4版（CCSM4）进行长时间尺度的冰期气候模拟；2) 基于格陵兰和南极冰芯（如NGRIP冰芯）中的硫酸盐沉积记录，重建过去火山活动的时序和强度；3) 应用广义极值分布（GEV）统计模型分析大型火山喷发的重现期；4) 采用大型集合模拟方法，通过扰动初始场来分离火山强迫和内部气候变率的影响；5) 结合简化的非线性动力系统模型，对复杂模型结果进行理论框架下的诠释。

RESULTS

Volcanic impacts on interstadial and stadial climate

研究结果显示，火山喷发对AMOC的影响具有强烈的状态依赖性。在典型的D-O循环中，暖的间冰阶（Interstadial）可持续数千年，期间AMOC强盛；而冷的冰阶（Stadial）则相对短暂。模型模拟表明，非常大型的赤道火山喷发（如230 Tg气溶胶）能够显著改变AMOC的演化轨迹。通过分析在不同D-O阶段启动的大型集合，研究人员发现，当气候系统处于间冰阶晚期、即将自然转向冰阶的“临界点”附近时，其对火山扰动的敏感性最高。此时，一次大型火山喷发可以显著提高AMOC从强状态崩溃至弱状态的概率。集合模拟显示，火山强迫后，AMOC通常会先出现一个短暂的增强（在几年内），随后在几十年内减弱至低于未受扰动集合的平均水平。这一“先增强后减弱”的响应模式是触发潜在状态转换的关键。有趣的是，即使在敏感的晚期间冰阶，也并非所有遭受火山扰动的模拟都会发生AMOC崩溃，部分模拟在经历初始扰动后仍会恢复到原有的强AMOC状态。这表明，气候系统的内部变率在火山触发突变的过程中扮演着双重角色：它既可能放大扰动效应，促使系统跨越临界点，也可能抵消扰动，增强系统的稳定性，尤其是在系统接近 tipping point 时。

Postvolcanic teleconnections and AMOC strengthening

为了理解火山喷发如何引发AMOC的初始增强和后续减弱，研究人员详细分析了北大西洋地区大气和海洋变量的演变。火山喷发后，全球尺度的辐射强迫赤字导致地表显著降温。这种冷却效应改变了大气环流，特别是在北大西洋地区，诱发了一个负位相的北大西洋涛动（NAO）型态，表现为冰岛低压和亚速尔高压的同时减弱。负NAO相减少了向高纬度地区的湿气输送，导致在关键的深水形成区（如格陵兰-冰岛-挪威海，GIN Seas）净淡水通量（P - E + R，即降水减蒸发加径流）减少，从而增加了海表盐度（SSS）。更咸、更冷的海水密度增大，有利于海洋深层对流，这是AMOC初始增强的直接原因。同时，次极地环流（Subpolar Gyre）的增强也促进了来自副热带盐度较高海水的向北输送，进一步强化了盐度正反馈。然而，这种增强是短暂的。几年后，持续的高纬度冷却促进了海冰的扩张。海冰的生长和夏季融化会在海洋上层形成稳定的卤跃层（Halocline），抑制深层对流。此外，后期北大西洋地区海表热浮力通量输入增加（尤其在冬季）和卤浮力通量输入增加（在夏季），共同导致海洋表层浮力增大，进一步阻碍对流。这些过程共同作用，最终引发AMOC的减弱，并可能启动海冰-深水形成的负反馈循环，推动系统向冰阶状态转变。

Dynamical systems perspectives on postvolcanic decadal climate variability

为了从更宏观的视角理解AMOC在火山强迫后几十年到上百年的长期演化行为，研究者借助了一个简化的随机动力系统模型。该模型将AMOC的动态抽象为一个具有双稳态（强AMOC和弱AMOC）的非线性系统，并引入了随机噪声来代表气候内部变率。通过调整模型参数（如控制参数-CO₂浓度、噪声强度和火山扰动幅度），这个简单模型成功地复现了复杂气候模型中观察到的许多关键特征，例如火山扰动后AMOC的类正弦响应、状态转换的概率性、以及高噪声水平下系统对火山触发“失灵”（即不易发生转换）的现象。分析表明，在本文模拟的冰期CO₂水平（230 ppm）下，气候系统可能处于一个稳定的强AMOC状态附近，但其稳定性受到内部变率的持续挑战。火山扰动可以看作一个巨大的“推力”，有可能将系统推过势垒（Basin Boundary），落入弱AMOC的吸引域。然而，系统内部固有的噪声同样可以产生类似的效果（噪声诱导跃迁），或者相反，在火山扰动后帮助系统“爬回”原来的稳定状态。这种理论框架将复杂的物理过程提炼为更本质的动力机制，强调了在评估气候 tipping points 时，考虑外部强迫与内部变率相互作用的重要性。

研究结论与意义

本研究通过创新的“古气候数据-复杂模型-理论模型”相结合的研究策略，有力地论证了大型赤道火山喷发作为冰期千年尺度气候突变（特别是D-O事件）潜在触发器的物理可能性。研究揭示，火山活动并非直接导致AMOC崩溃，而是通过一系列级联的大气-海洋遥相关过程（如引发负NAO相、改变北大西洋盐度和海冰分布），显著改变了AMOC的稳定性，特别是在气候系统处于临界状态（如间冰阶晚期）时，能够大大提高状态转换的概率。更重要的是，研究强调了气候系统内部变率的关键作用，它既可以促进也可以抑制火山触发的突变，这使得气候响应呈现出概率性而非确定性的特点。

这项研究的发现具有深远的意义。首先，它为解读冰芯等古气候档案中记录的气候突变事件提供了一个新的解释框架，将火山活动这一外部强迫与气候系统的内部动力学联系起来。其次，研究中所揭示的AMOC对巨大扰动的响应机制，以及系统接近 tipping point 时表现出的复杂行为，为评估现代和未来气候变化背景下AMOC的稳定性提供了宝贵的历史相似型。在当前 anthropogenic 变暖可能导致AMOC减弱甚至崩溃的担忧日益加剧的背景下，理解其在过去极端气候条件下的行为模式至关重要。本研究指出，未来的研究应更关注火山活动集群（而非孤立事件）的可能影响，以及气候背景（如冰盖规模、温室气体浓度）如何调制系统对火山强迫的敏感性。总之，这项工作深化了我们对地球气候系统脆弱性和恢复力的认识，凸显了在非线性系统中，微小扰动可能引发巨大后果的潜在风险。