NCC | 全球变暖背景下南极冰盖流域临界点风险制图

NCC|全球变暖背景下南极冰盖流域临界点风险制图

Fig. 1 | Risk map of Antarctic ice catchment basins.

Fig. 1 | Risk map of Antarctic ice catchment basins.

Fig. 4 | Retreat of ice fronts with increasing warming.

Fig. 4 | Retreat of ice fronts with increasing warming.

Fig. 5 | Risk assessment of critical threshold temperatures and respective sealevel commitments.

Fig. 5 | Risk assessment of critical threshold temperatures and respective sealevel commitments.

一、研究背景与科学意义

1.1 南极冰盖稳定性与气候临界点

南极冰盖是地球系统中质量最大的冰体，其总质量相当于全球海平面上升潜力约60米。在全球变暖背景下，南极冰盖正通过一系列放大反馈机制加速质量损失，这使其成为气候系统中最关键的"临界点元素"（tipping element）之一。所谓临界点，是指当气候强迫超过特定阈值后，系统内部自维持的反馈过程将导致不可逆的、 often abrupt 的状态转变。

古气候记录表明，在过去仅比现代高1-3°C的间冰期温暖期，南极冰盖曾贡献数米量级的全球海平面上升，暗示西南极（WAIS）和东南极（EAIS）的海洋性冰盖区域曾发生大规模退缩。特别是末次冰消期，南极冰盖退缩导致的融水脉冲事件曾造成海平面以每世纪约0.7米的速率上升。当前全球平均气温已较工业化前水平上升约1.3°C，而观测数据显示西南极阿蒙森海湾区（Amundsen Sea Embayment）已出现由海洋驱动冰架减薄导致的接地线（grounding line）后退和冰流加速，显示该区域可能已处于不稳定状态。

1.2 冰盖动力学的复杂性与尺度问题

南极冰盖的动力响应并非均一。传统上，学界常将南极冰盖视为单一或两个（东西南极）临界点元素，但这种粗粒度视角掩盖了不同流域（drainage basins）固有的地形和气候异质性。实际上，南极冰盖由18个主要排水流域组成，每个流域具有独特的床基地形（bed topography）、冰厚分布、海洋强迫条件和冰流动力学特征。

其中最关键的动力机制包括：

• • 海洋性冰盖不稳定性（MISI, Marine Ice-Sheet Instability）：在冰床向海倾斜（retrograde sloping bed）的区域，接地线后退导致冰通量增加，形成正反馈
• • 冰崖不稳定性（MICI, Marine Ice Cliff Instability）：高大垂直冰崖的机械崩塌可能加速冰量损失，但观测约束仍有限
• • 表面消融-高程反馈：表面融化降低冰面高程，使其暴露于更暖的大气温度，加剧融化
• • 冰川均衡调整（GIA, Glacial Isostatic Adjustment）：冰量减少导致地壳抬升，可能稳定接地线位置

本研究的核心科学问题在于：不同流域在全球变暖梯度下的稳定性行为是否存在本质差异？哪些流域呈现渐进式衰退（gradual decline），哪些存在明确的温度阈值（tipping points）？

二、研究方法与数据基础

2.1 数值模型架构

研究团队采用**并行冰盖模型（PISM, Parallel Ice Sheet Model）**v1.0版本进行长期平衡态模拟。该模型为热-力学耦合、多热力学状态的冰流模型，采用浅冰近似（SIA）与浅架近似（SSA）的混合应力平衡方案，确保从冰盖内部（垂直剪切流主导）到冰流（底滑主导）再到冰架（快速塞流）的平滑过渡。

关键模型配置参数：

• • 水平分辨率：16 km（采用亚网格尺度线性插值处理接地线迁移）
• • 垂向网格：非均匀分层，底部20 m至顶部100 m
• • 滑动定律：广义伪塑性幂律（指数q=0.75），基于Mohr-Coulomb准则计算屈服应力
• • 崩解参数化：基于伸展率的"特征崩解"（eigencalving）方案，最小冰架厚度阈值50 m

地球系统反馈：

• • 冰川均衡调整：采用粘弹性地壳形变模型，上地幔粘度设为10²¹ Pa·s，岩石圈密度3300 kg/m³
• • 气候边界条件：使用度日模型（PDD）计算表面融化（雪/冰融化系数分别为3/8 mm PDD⁻¹），积累率基于RACMO v2.3p2再分析数据，并引入温度与海拔的指数型气候校正（Clausius-Clapeyron关系，降水敏感度5%/K）

海洋强迫：

• • 冰架底部融化采用PICO（Potsdam Ice-shelf Cavity Model）模型
• • 驱动数据：1975-2012年观测海洋温盐数据（Schmidtko数据集）
• • 南极区域海洋温度变化与全球平均温度变化的缩放系数：0.7（基于ECHAM5/MPIOM等气候模式的长期平衡实验）

2.2 准静态升温实验设计

本研究采用准静态（quasi-static）方法论，这是识别真实临界点的前提条件：

1. 1. 初始状态：从近似工业化前（pre-industrial）的平衡态启动，该状态经过250 kyr的spin-up（包括Bedmap2地形、RACMO气候强迫等），与现代观测冰几何形态具有合理一致性（冰面高程RMSE=137.2 m，表面速度RMSE=171.9 m/yr）
2. 2. 升温协议：以0.0001°C/yr的极慢速率线性增加全球平均温度（慢于冰盖典型响应时间尺度），确保系统始终接近平衡态
3. 3. 分阶段平衡：在每个整度升温水平（1°C, 2°C...直至完全去冰）暂停升温，运行至少20 kyr（高变暖情景>6°C时运行50 kyr）直至体积变化可忽略，获得真实稳态
4. 4. 流域尺度分析：基于MEaSUREs定义的18个排水流域边界（IMBIE-3标准），分别追踪每个流域的海平面相关冰体积（volume above flotation, VAF）变化

三、关键研究发现：渐进衰退与临界点动力学

3.1 动态响应的二分性

研究揭示南极冰盖对变暖的响应存在两种截然不同的动力学模式：

渐进式衰退（Gradual Decline）：如Abbot/Venable、西南极半岛（George VI）和Ronne流域，冰体积随升温几乎线性减少，或在多个阶段逐步丧失，无明显突变阈值。

临界点动力学（Tipping Dynamics）：如Dronning Maud Land、Enderby Land、Amery、West/Denman、Totten/Moscow、Filchner（东南极）以及Thwaites/Pine Island（西南极）。这些流域在较宽的升温范围内变化微小，一旦超过特定阈值，发生突发的大规模体积损失，在极小的额外升温下即接近完全丧失。

混合模式：部分流域如Ross West (Siple Coast)、Cook/Ninnis/Mertz、Victoria Land在1-3°C出现突变后转为渐进衰退；而Brunt/Riiser-Larsen、Ross East (Byrd)、Getz等在低升温下渐进衰退，但在>6°C时出现阈值行为。

3.2 临界温度的识别标准

研究采用峰度 prominence 方法识别临界温度：计算每度升温导致的相对冰损失率（%/°C），识别超过初始体积15%显著性阈值的峰值。对于双峰结构（如Amery、West/Denman、Getz），将两度区间均视为临界。

西南极高风险区：

• • Thwaites/Pine Island：临界阈值<1°C，当前全球已升温1.3°C，暗示该区域可能已跨越临界点。一旦触发，长期承诺损失约70%初始体积（~0.9 m海平面当量，SLE），主要由MISI驱动，接地线在逆向倾斜床基上不稳定后退，直至正向坡度区域才减缓。
• • Ross West (Siple Coast)：阈值<1°C，承诺损失约1.0 m SLE
• • Ronne：第一阈值<1°C，第二阈值6-7°C

东南极关键区：

• • Cook/Ninnis/Mertz（Wilkes次冰盆）：第一阈值2-3°C，承诺损失40%初始体积（1.2 m SLE）；第二阈值6-7°C，额外损失30%（0.9 m）。机制上，沿海"冰塞"（ice plug）的微小扰动可触发Wilkes盆地的长期自放大退缩。
• • Victoria Land：与Wilkes盆地动力连通，在4-5°C时部分冰量通过该盆地排入南大洋，额外贡献~0.5 m SLE
• • Totten/Moscow、Filchner、Ross East (Byrd)：阈值6-8°C，但海平面承诺量巨大（Filchner达11.8 m SLE）

四、综合风险评估与气候政策含义

4.1 燃烧余烬风险图（Burning Embers）

研究构建了南极版的"燃烧余烬"图（Fig.3），以温度-海平面承诺矩阵展示各流域风险：

极高风险区（深红）：

• • Thwaites/Pine Island、Ross West：当前气候条件下已处于或跨越临界阈值
• • 这些流域占西南极海平面相关冰体积的40%（2.1 m SLE）已在1-2°C升温下被"承诺"（committed）

高风险区（橙红）：

• • Cook/Ninnis/Mertz：在SSP2-4.5情景下（本世纪末升温~2-3°C）面临临界风险
• • Abbot/Venable、西南极半岛：阈值3-5°C，在SSP2-4.5及更高排放情景下可能发生

中高风险区（黄色）：

• • 东南极大部分流域（Filchner、Amery、Totten等）：阈值>6°C，但累积海平面承诺>26 m

4.2 临界点级联与相互作用

研究强调南极冰盖应被视为相互作用的临界点网络而非单一整体。关键发现包括：

1. 1. 东西南极不对称性：西南极整体比东南极更脆弱，但东南极的Wilkes盆地和Aurora次冰盆在适度升温（2-5°C）下即存在风险，贡献~5 m潜在海平面上升。
2. 2. 临界点的不可逆性：由于滞后效应（hysteresis），一旦跨域阈值，即使温度回落至工业化前水平，冰盖也无法完全恢复。例如，西南极冰盖除非温度较工业化前降低至少1°C，否则无法重新生长至现代范围。
3. 3. 时间尺度的分离：部分变化可能在未来数十年触发，但完全展开需数世纪至数千年。这种"承诺效应"意味着当前决策将锁定未来多代人的海平面上升命运。

五、研究局限性与不确定性分析

5.1 模型结构不确定性

未包含的物理过程：

• • MICI（海洋性冰崖不稳定性）：本研究未纳入此机制，因其过程理解不足且观测约束有限。若MICI被激活，可能显著降低多个流域的临界阈值并加速冰损失。近期研究对21世纪MICI的重要性存在争议。
• • 冰-反照率反馈：表面融化导致的反照率变化在低升温情景中可能较次要，但在>6°C情景下重要性增加。
• • 完全耦合的海洋-大气相互作用：使用离线气候边界条件（固定气候模式缩放因子）无法捕捉冰-气候系统的双向反馈，如融水输入对南极底层水形成的影响。

5.2 参数敏感性

研究通过敏感性实验验证了主要结果的稳健性：

• • GIA参数：上地幔粘度降低一个量级（10²⁰ Pa·s）或考虑弹性岩石圈响应，可使阈值向高温方向偏移<1°C，但不改变定性行为
• • PDD参数：融化系数或温度标准差变化在合理范围内，导致的温度偏移<1°C
• • 降水敏感度：4-6%/K的变化范围使阈值在相同温度区间内保持稳定
• • 动力参数：SSA速度增强因子、滑动律指数、冰下融水衰减率等变化主要影响达到平衡的速率，而非最终稳态（在2°C和7°C实验中，所有参数组合最终收敛于两个邻近平衡态之一）

5.3 空间分辨率限制

16 km分辨率足以捕捉MISI等关键反馈的大尺度动力学，但可能无法精确解析接地线的小尺度几何形态和局部海洋环流对冰架底部的强迫。不过，对于识别平衡态之间的临界阈值（而非瞬态路径），分辨率敏感性相对较低。

六、讨论与科学启示

6.1 对IPCC评估的补充

本研究为IPCC第六次评估报告（AR6）提供了重要补充。AR6指出到2100年南极可能贡献3-34 cm海平面上升（可能范围），到2300年不排除7-14 m。本研究通过平衡态承诺（equilibrium commitment）视角，量化了不同升温水平下各流域的长期"锁定"损失，强调了即使升温停滞，海平面仍将持续上升数个世纪。

6.2 早期预警指标的启示

由于南极冰盖由多个具有不同临界温度的子系统组成，寻找"正确尺度"的早期预警指标至关重要。对于Thwaites/Pine Island等已接近或跨越临界点的流域，应重点关注接地线位置、冰架底部融化速率、冰流速度等观测指标；而对于东南极流域，需监测沿海冰塞的稳定性和次冰盆尺度的海洋热含量变化。

6.3 气候治理含义

研究结果对《巴黎协定》温控目标具有直接政策含义：

• • 1.5°C目标：即使实现此目标，西南极部分流域（Thwaites/Pine Island、Ross West）的长期崩溃可能已无法避免，需为至少0.9-1.0 m的长期海平面上升做准备
• • 2°C门槛：Cook/Ninnis/Mertz等东南极关键流域将面临临界风险，额外承诺~1.2 m海平面上升
• • >3°C情景：将导致多流域级联临界点跨越，累积承诺可能超过数米

七、未来研究方向

7.1 过程模型改进

• • 开发显式解析MICI过程的模型，并约束其触发条件（冰崖高度阈值、冰架崩解前缘力学）
• • 引入冰-海洋-大气的全耦合相互作用，特别是融水羽流对冰架空腔环流的影响
• • 改进冰下沉积物（till）力学和冰下湖泊排水过程的参数化

7.2 观测与数据需求

• • 加强东南极Wilkes盆地、Aurora次冰盆的床基地形测绘（当前数据仍存在较大不确定性）
• • 建立多流域尺度的长期海洋-冰架相互作用观测网络
• • 发展冰流速度的实时监测技术以检测临界减速（critical slowing down）等预警信号

7.3 情景与风险评估扩展

• • 开展全球变暖"过冲"（overshoot）情景下的冰盖恢复实验，量化临时超调是否导致永久性损失
• • 开发基于机器学习的快速替代模型（emulators），以覆盖更广泛的参数不确定性和气候情景
• • 将南极临界点与全球沿海影响模型耦合，评估不同流域贡献的时间-空间差异对特定沿海社区的影响

八、结论

本研究通过高分辨率冰盖模型和创新的准静态实验方法，首次系统性地绘制了南极18个主要冰盖流域在全球变暖梯度下的临界点风险图谱。核心结论如下：

1. 1. 南极冰盖是多临界点系统：不同流域呈现从渐进衰退到突发临界点的多样化动力学行为，不能简化为单一临界元素。
2. 2. 西南极已处临界边缘：Thwaites/Pine Island、Ross West和Ronne流域的临界阈值低于1°C，当前1.3°C的全球升温可能已触发其长期崩溃，承诺损失2.1 m海平面当量。
3. 3. 东南极存在中等升温风险：Cook/Ninnis/Mertz等Wilkes盆地相关流域在2-3°C下存在明确临界点，贡献~1.2 m潜在海平面上升；而Totten、Filchner等流域在>6°C下面临风险，但累积承诺>26 m。
4. 4. 不可逆性与承诺效应：由于滞后行为和长响应时间尺度，即使未来降温，已触发的冰损失也将在数百年至数千年尺度上持续，对沿海 civilization 构成前所未有的长期挑战。

作者与单位信息

主要作者：

• • Ricarda Winkelmann（通讯作者）：德国波茨坦气候影响研究所（PIK）地球韧性科学部、马克斯·普朗克地球人类学研究所整合地球系统科学部、波茨坦大学物理与天文研究所
• • Julius Garbe：德国波茨坦气候影响研究所地球韧性科学部、波茨坦大学物理与天文研究所
• • Jonathan F. Donges：德国波茨坦气候影响研究所地球韧性科学部、马克斯·普朗克地球人类学研究所整合地球系统科学部、斯德哥尔摩大学斯德哥尔摩韧性研究中心
• • Torsten Albrecht：德国波茨坦气候影响研究所地球韧性科学部、马克斯·普朗克地球人类学研究所整合地球系统科学部

资助机构： 欧盟地平线2020计划（TiPACCs、PROTECT、OCEAN ICE、ClimTip项目）、德国研究基金会（DFG）、德国联邦教育与研究部（BMBF）、欧洲研究委员会（ERC）、美国国家航空航天局（NASA）、美国国家科学基金会（NSF）等。

数据与代码可用性： 模拟数据与代码均遵循开放科学原则，通过Zenodo等平台公开获取。

以上是对文章的详细解读，如有不当之处欢迎批评指出！

END

声明：欢迎转载、转发。气象学家公众号转载信息旨在传播交流，其内容由作者负责，不代表本号观点。文中部分图片来源于网络，如涉及内容、版权和其他问题，请联系小编处理。