《Atmosphere》：Long-Term Volcanic Signal in 21st-Century Climate Projections with a 25-Member Stochastic Ensemble Using SOCOL-MPIOM
编辑推荐：
未来火山喷发在CMIP6模拟中基本被忽略，从而增加了21世纪气候投影的不确定性。研究人员利用SOCOL-MPIOM气候模型开展了为期80年（2020–2100）的25成员随机集合模拟，驱动场为SSP3-7.0强迫情景，并在每个集合成员中引入五次随机分布的热带喷
未来火山喷发在CMIP6模拟中基本被忽略，从而增加了21世纪气候投影的不确定性。研究人员利用SOCOL-MPIOM气候模型开展了为期80年（2020–2100）的25成员随机集合模拟，驱动场为SSP3-7.0强迫情景，并在每个集合成员中引入五次随机分布的热带喷发——三次强喷发、一次中等喷发和一次弱喷发（以下称SV试验）。研究人员通过将SV试验结果与相同人为强迫情景下的无火山基线模拟进行对比，分析火山影响。80年模拟期内的五次喷发使主要气候指标的全球和年平均尺度趋势与无火山基线相比在统计上无法区分。然而，在区域和季节尺度上，火山活动可显著改变无火山模拟的结果。例如，在北欧地区，火山喷发可产生高达1.0 K的冬季温度增暖（约为参考试验增暖的30%）以及36 mm yr?1的年降水量亏缺。这强调了将火山喷发纳入未来气候投影以实现更精准预测的必要性。SV集合的概率分析显示，年最大日温度（TXx）以"较可能"概率在全球16%的陆地上超过+0.5 K的扰动，这一扰动在标准CMIP6结果中并不存在。由于本研究的情景组合针对21世纪合理火山活动范围的上限，这些超阈值区域代表的是近最大值估计而非最可能估计。
研究背景与问题：火山喷发是地球气候系统的重要自然强迫因子，但当前气候模型对其在未来气候投影中的作用重视不足。爆破式火山喷发将二氧化硫（SO
2)注入平流层，经氧化形成硫酸盐气溶胶层，通过增强行星反照率、抑制地表短波辐射以及吸收太阳入射辐射和地球长波辐射，造成下平流层增温而对流层和地表冷却。1991年皮纳图博（Pinatubo）火山喷发作为仪器观测时代最大的平流层硫注入事件，产生了约0.5 K的地表和对流层降温，伴随全球降水减少和臭氧柱扰动。此外，气溶胶引发的热带下平流层加热会加剧赤道至极地温度梯度，强化极地涡旋，有利于喷发后首个冬季出现正位相北极涛动（Arctic Oscillation, AO），使北半球大陆增暖并在欧亚大陆北部形成异常西风平流。然而，现有研究多聚焦于灾难性单一喷发事件，对21世纪多次中等强度喷发叠加的累积气候效应量化不足。大多数CMIP6积分要么忽略未来火山强迫，要么仅指定恒定背景气溶胶气候态，低估了可能气候轨迹的不确定性范围。同时，尚未有世纪尺度的随机火山集合研究在全大气层范围内耦合了交互式臭氧化学与规定的火山气溶胶强迫。火山强迫叠加在变暖背景态上可能在北极海冰边缘和季风带等阈值敏感系统中产生非线性响应，而气候变化本身也会改变火山气溶胶的平流层寿命和辐射效率，这是单一喷发实验无法捕捉的又一非线性来源。
研究开展与核心结论：为填补上述空白，研究人员利用SOCOL-MPIOM化学-气候模型，在SSP3-7.0排放情景下开展了2020–2100年80年期的25成员随机集合模拟。该研究首次在全耦合化学-气候模型中表征了世纪尺度的火山气候信号，同时解析了交互式平流层臭氧化学、规定气溶胶辐射和异相化学强迫、海洋-海冰耦合以及区域气候极端事件。每个集合成员包含五次独立抽样幅度和时间的喷发。通过减去无火山基线来分离火山信号，并在背景臭氧和温度趋势、垂直和纬向结构、区域地表气候、冰冻圈反馈、由气候变化检测和指标专家组（Expert Team on Climate Change Detection and Indices, ETCCDI）定义的气候极端事件以及暖阈值概率超阈值等方面进行了检验。研究发表在《Atmosphere》期刊，其重要意义在于建立了21世纪随机火山活动对气候影响的概率表征框架，揭示了火山强迫在区域和季节尺度上引入的、当前标准气候投影未能捕捉的实质性扰动。
关键技术方法：研究人员采用的SOCOL-MPIOM模型通过OASIS3界面将大气化学-气候模型SOCOL与海洋-海冰模型MPIOM耦合，每日交换场数据。大气组分基于ECHAM5.4环流模型构建，扩展了MEZON交互式平流层化学方案，处理41种化学物种在140种气相、光解和异相反应中的相互作用，谱分辨率为T31（约3.75°×3.75°），垂直39层从地表延伸至0.01 hPa（约80 km）。海洋组分MPIOM采用约3°水平分辨率的双极网格，40个垂直层。平流层硫酸盐气溶胶尺度分布（SAD）和光学特性被规定以考虑辐射强迫和异相臭氧化学。
实验设计以单个无火山强迫的基线模拟为对照，在SSP3-7.0情景下运行2020–2100年。火山集合（SV）包含25个独立成员，每个成员在2020–2100年间实现五次热带喷发。情景构建依据Toohey等人编制的10至20世纪喷发记录，将火山喷发按最大平流层气溶胶光学深度（AOD）分为强（AOD>0.10）、中（0.05
火山信号定义为各成员与基线的差值，空间和时间分析采用集合平均信号。统计显著性通过25个成员异常值的一样本t检验评估，同时使用双侧Student's t检验（5%显著性水平）评估各格点信号显著性，采用普通最小二乘法估计长期线性趋势并用Mann-Kendall检验评估显著性。概率超阈值图基于25成员分布构建，计算超过给定阈值τ的成员比例，并以50%（"较可能"）和75%（"高信度"，遵循IPCC AR6稳健性标准）等值线表示。
气候变化极端指数采用ETCCDI的子集：年最大日最高温度（TXx）、年最小日最低温度（TNn）、日较差（DTR)、暖日比例（TX90p）、霜冻日数（FD）、夏日数（SU）和年总降水量（PRCPTOT）。区域平均覆盖八个关键域：北极（>66°N）、北欧（50–70°N, 10–40°E）、中亚（35–55°N, 60–90°E）、东亚（25–45°N, 100–140°E）、热带（30°S–30°N）、北美（25–60°N, 130–60°W）、南半球（<30°S）和南极（<60°S）。此外，研究人员对20°N以北的月近地表温度（T2m）和大气表面气压（APS）分别进行冬季（DJF）和夏季（JJA）的经验正交函数（EOF）分析，将火山信号投影到主导模态以量化其激发或抑制程度。
研究结果：
背景趋势：火山与非火山模拟的对比。2020–2100年间，SV集合平均与基线在所有四个变量上紧密跟踪；反复随机喷发未在SSP3-7.0下对气候系统的长期轨迹留下可检测印记。平流层臭氧在两种试验中均增加（SV为+0.009 ppmv dec
?1，基线为+0.006 ppmv dec
?1），与《蒙特利尔议定书》下的臭氧层恢复一致。对流层臭氧分别以?1.440和?1.459 ppbv dec
?1下降，平流层温度分别以?0.420和?0.411 K dec
?1降低，地表温度分别以+0.338和+0.341 K dec
?1上升。五次喷发——构成21世纪合理火山活动上限——在SSP3-7.0下未在任何全球平均变量的百年轨迹上留下可检测印记，证实全球尺度火山信号未表现为长期趋势偏移。
火山信号的时间演变与垂直结构。2020–2100年间火山信号显示显著的年际变率，叠加于弱的正平流层偏移和弱的地表负偏移之上。平流层温度异常平均为+0.181 K，2060年达峰值+0.626 K，与气溶胶辐射加热一致。全球平均信号在约18 km处达峰值+0.4 K，热带纬度经历最大增暖（+0.7–1.0 K）。对流层10–6 km处存在?0.3至?0.5 K的显著冷却。臭氧垂直结构呈双层：热带下平流层（16–20 km）为负异常（峰值?0.02 ppmv），极地中平流层（30–45 km）为正异常（2100年NH极地达+0.05 regional）。纬向平均截面显示，热带损耗（达?40 ppbv）局限于70–200 hPa层，30°N以北和50 hPa以上信号反转，NH极地中平流层达+40 ppbv。
空间分布与季节依赖性。总臭氧柱异常（ΔTOZ）呈现热带减少2–5 DU、高纬度增加3–8 DU的带状结构。地表温度异常在北半球高纬度以正异常为主，北极和北欧大陆达+0.5至+1.5 K。季节分解揭示两种物理机制：冬季（DJF）北欧和欧亚大陆增暖+1.0至+2.0 K（72%域点可检测），与平流层遥相关路径一致——气溶胶加热热带下平流层，强化极涡并驱动异常西风波涛输送；冬季北欧降水弱正异常+2.0 mm season
?1。夏季（JJA）平流层路径 inactive，热带和副热带陆地冷却?0.5至?1.0 K，北欧JJA降水减少?10.6 mm season
?1，年度净亏缺?36 mm yr
?1。
冰冻圈反馈机制。高纬度地表响应受冰冻圈反馈链控制：火山减少地表短波辐射（与北极海冰面积相关系数r=?0.50，p<0.001），云量异常与地表短波相关（r=?0.63，p<0.001）。北极海冰范围平均增加+0.44%（p<0.001），最大异常在巴伦支海和卡拉海。表面反照率与北极最低温度相关达r=?0.89（p<0.001）: 火山地表冷却驱动海冰和积雪扩张，通过提高反射率放大初始冷却。
气候极端事件中的火山信号。全球平均ETCCDI值很小，但区域结构高度异质。北欧在所有10个指数上显示显著异常，包括TXx +0.74 K、TNn +0.37指数为+0.37 K、DTR +0.10 K、SU +12%、FD ?7.3%、ID ?13%、TR +26%、TX90p +3.15%、TN10p ?1.79%、PRCPTOT ?36 mm yr
?1——分析中最大的区域降水异常。北极显示最大的夜间增暖（TNn +0.66 K），通过冰冻圈反馈放大。中亚以夜间冷却为主（TNn ?0.73 K），与北欧形成动力学对比。东亚显示日间增暖（TXx +0.34 K）。热带为唯一日间一致冷却区域（TXx ?0.18 K）。北美未出现任何统计显著响应。
概率评估。25成员集合允许将火山信号表达为超概率。对于TXx超过+0.5 K，50%阈值被全球约16%的陆地超过，集中在北北欧欧亚和北极地区；75%轮廓覆盖北欧和的西伯利亚。提升至+1.0 K，50%概率面积降至4%；+1.5 K时无格点达标。北极海冰+1%增加概率在巴伦支海、卡拉海和拉普捷夫海部分区域超过50%，75%轮廓覆盖约1%的北极域。海冰响应空间相干性高于TXx，反映气溶胶短波减少与海冰面积之间紧密物理耦合（r=?0.89）。
EOF分析：北半球温度和气压模态。DJF温度第一EOF解释19.5%年际方差，呈现北大西洋-欧洲扇区负载荷和俄罗斯北部及北极正载荷的区域偶极结构。火山异常投影（ΔPC1_T2m）在DJF平均为负（?3.69 K），表明集合系统性地使近地表温度场向暖北极、相对冷却北大西洋-欧洲模式偏移。APS第一EOF解释32.4%方差，呈现标准的AO气压空间环形结构。火山异常投影（ΔPC1_APS）在DJF和JJA均在零附近振荡，无任何季节统计显著平均偏移。T2m的显著非零ΔPC1_T2m与APS的非显著ΔPC1_APS之对比表明，北半球冬季温度场的主导火山影响并非源于一致的环形环流偏移，而是由海冰和反照率反馈驱动的热力学信号。
讨论与结论总结：该研究确立了火山活动对21世纪气候嵌入的长期信号的特征化。尽管全球和年平均尺度上火山信号与无火山基线无法区分，但区域和季节尺度效应显著且物理一致。北欧响应最为一致，整合冬季平流层遥相关增暖和夏季气溶胶抑制对流的干旱，形成最大的区域降水亏缺。北极以热力学反馈为主：气溶胶短波减少扩张海冰，提高地表反照率放大初始冷却，反馈链中反照率与最低温度异常相关系数达r=?0.89。EOF分析首次证明该暖北极信号不投影到AO环流模态，将冰冻圈放大路径与平流层遥相关分离。北美无可检测信号反映竞争动力学路径的真实抵消，而非集合规模不足，构成了中等火山强迫地理触及范围的物理约束。全球尺度上，对流层臭氧趋势在SV集合与基线间无法区分，平流层极区富集不传播至最低平流层，排除了增强平流层-对流层交换作为对流层臭氧预算修改机制的假说。平流层臭氧呈现双峰异常——热带损耗与极区富集配对，随卤素负荷下降向2100年增强，使火山臭氧响应从异相化学破坏向动力学极地传输转变。概率框架量化了这一影响的实际边界：TXx超过+0.5 K在16%全球陆地较可能（>50%概率）发生。由于情景针对合理火山活动上限，这些为近最大值估计，但确立了随机火山活动引入区域实质性、物理一致且可概率量化的扰动，这是当前CMIP6投影未能捕捉的。
