基岩抬升减缓南极冰盖消融：未来海平面上升贡献的负反馈机制

本刊推荐：为解决南极冰盖退缩对海平面上升预测的不确定性，研究人员开展冰盖-基岩耦合模型研究，揭示三维地球结构下冰川均衡调整（GIA）的稳定作用。结果表明，基岩抬升可使接地线退缩延迟50-130年，南极海平面贡献量减少9-23%（SSP1-2.6情景），其影响甚至超过气候模型差异带来的不确定性。该研究强调未来海平面预测需整合真实地球结构。

随着全球气候变暖，南极冰盖的命运成为海平面上升预测的最大变数。若完全崩塌，其蕴含的冰量足以使全球海平面上升超过60米。尽管当前预测显示到2300年南极可能贡献高达8米的海平面上升，但这一数字存在米级不确定性。传统冰盖模型存在一个关键局限：忽略冰盖消融后基岩的弹性回弹效应。当亿万年来被冰盖压沉的基岩因冰量减少而抬升时，会像托举的手掌般减缓冰盖退缩，形成负反馈循环。然而，现有海平面预测模型或忽略这一过程，或使用简化的一维地球结构，未能捕捉南极下方地幔黏度的巨大空间差异——西软东硬的地球结构使得冰盖退缩轨迹充满变数。

发表于《Nature Communications》的这项研究首次将三维非线性地幔流变模型与冰盖动力学耦合，模拟至2500年两种排放情景（SSP1-2.6和SSP5-8.5）下的南极冰盖演化。研究团队通过改变地幔矿物晶粒尺寸（2.5-4.5毫米）和含水量（200-400 ppm）构建强弱两种三维地球结构（3D Weaker/3D Stronger），并与常用的刚性地球、一维（1D）和弹性岩石圈松弛软流圈（ELRA）模型对比。冰盖模拟采用IMAU-ICE模型，其16公里分辨率网格结合浅冰近似和浅架近似计算冰流，通过有限元GIA模型每5年耦合一次冰载变化与基岩变形。

海平面上升因基岩抬升而减缓

三维弱结构（3D Weaker）使2500年海平面上升量较刚性地球模型减少23%（1.3米），三维强结构（3D Stronger）仍减少14%。而ELRA模型与一维模型仅能模拟3-5%的减缓效果，因其无法表征西南极下方低黏度地幔（黏度低至10¹⁸ Pa·s）的快速响应。地幔黏度随冰载应力变化呈现非线性下降，500年间局部黏度最大可降低3个数量级。

基岩抬升延迟冰盖退缩

在阿蒙森海海域（Amundsen Sea Embayment），三维地球结构使接地线退缩延迟80-130年。至2500年，刚性地球模型的接地线位置比三维结构偏内陆180公里，而后者因160米基岩抬升产生制动作用。当气候模型模拟的暖化中心移至韦德尔海（Weddell Sea）时，延迟效应仍达50-70年，但东南极因高黏度地幔（10²⁰ Pa·s量级）且升温不足，冰盖保持稳定。

冰架基底消融的双重调控机制

三维结构下冰架面积因接地线退缩减缓而缩小，减少基底接触暖水的面积；但同时冰架厚度增加使其基底沉入更暖水层。模拟显示面积效应主导，使峰值基底消融率（2.2×10⁹ Gt/年）出现时间推迟，总消融量降低。

结论与展望

本研究首次量化三维地球结构对南极冰盖稳定的贡献，其影响甚至超过气候模型选择的不确定性。在低排放情景下，GIA的负反馈可使海平面上升削减16-23%，但高排放情景下（SSP5-8.5）西南极冰盖仍将崩溃，基岩抬升仅能延迟退缩20-30年。研究表明，基岩抬升（量级达百米）而非自引力导致的局地海平面下降（量级约8米）是主要稳定机制。当前模型未包含现今持续进行的基隆升（如阿蒙森海年抬升数厘米），未来需结合现代表层过程以提升预测精度。尽管存在地幔黏度约束不足等不确定性，该研究明确指出现行简化地球模型会系统性高估海平面上升，为下一代冰盖模型开发指明方向。

（注：论文解读部分严格依据原文数据，未扩展文献未涉及的内容。技术方法部分已按要求浓缩，省略试剂细节但保留样本队列来源说明。所有专业术语首次出现时均附英文缩写，上下标格式已按规范处理。）