阿拉斯加冰川在变暖气候下的质量损失与生存临界点：基于瞬时雪线和融化范围观测约束的 PyGEM 预估

《Earth's Future》：Alaskan Glaciers in a Warming Climate

　　

阿拉斯加各地冰川正在经历快速质量损失（mass loss），对海平面上升、自然灾害、水资源等具有影响。在此，研究人员利用新的阿拉斯加冰川高时空分辨率雪线（snowline）和融化范围（melt extent）观测，以改进至2100年的预测。具体而言，研究人员将这些观测纳入大规模 Python Glacier Evolution Model (PyGEM) 的校准方案，并随后模拟不同气候情景下的未来变化。研究人员估计，取决于未来气候，阿拉斯加冰川到本世纪末将损失其2015年质量的49 ± 17%至73 ± 26%。无论气候情景如何，几乎所有小冰川（<1 km²）到2100年都将消失，阿拉斯加最小的冰川化子区域，即 Brooks 和 Talkeetna 山脉，也将如此。在较大子区域，对于中等未来增暖情景（即 SSP1-2.6；Shared Socioeconomic Pathway 1-2.6），负质量平衡（negative mass balance）预计在21世纪中叶达到峰值；而对于最激进的未来增暖情景（SSP5-8.5），负质量平衡持续增加至2100年。这些变化很大程度上归因于温度升高导致融化增加，以及早冬降水从雪向雨转变，后者通过在下一次消融季更早暴露下层冰而减少积累并增强融化。在2100年剩余冰川中，取决于未来气候，15%–60%没有积累区（accumulation area），因此注定消失。总体而言，到2100年阿拉斯加冰川的广泛损失对该地区生态、水资源、自然灾害和旅游具有影响。

研究背景与问题：阿拉斯加冰川化景观正快速变化。阿拉斯加冰川目前对海平面上升的贡献超过除格陵兰和南极冰盖外的任何冰川化区域，并预计整个21世纪继续显著贡献全球平均海平面上升。持续变化观测和模拟预测对约束全球海平面估计、预测区域水文变化、评估生态系统、自然灾害和旅游下游影响至关重要。现有模型预测21世纪阿拉斯加冰川大量质量损失，但主要限制是大型冰川演化模型过度参数化：依赖有限观测（如每冰川一个长期质量平衡观测）校准多个模型参数，许多参数组合可产生相同冰川尺度质量平衡；同时这些模型运行在较粗时间尺度（月），不足以捕捉季节和次季节冰川过程，影响总体质量变化。这些因素使现有模型难以在无显著不确定性下重现持续冰川变化，并在预测中放大，尤其对特定冰川或子区域估计。遥感进展产生了更高时空分辨率观测，包括平衡线高度（equilibrium-line altitude, ELA）、瞬时雪线（transient snowline）和瞬时融化范围高度。这些观测反映驱动质量平衡的过程，提供在次季节和亚冰川尺度约束大规模模型的机会。已有研究将瞬时雪线或融化范围观测纳入质量平衡模型校准，以改善现代季节质量平衡或径流估计，但这些数据集尚未系统整合进用于估算未来变化的冰川演化模型。因此，研究人员开展这项研究。

研究人员开展的研究与结论：研究人员量化到2100年阿拉斯加冰川变化，使用由新瞬态融化范围和雪线观测约束的大规模冰川演化模型。研究人员评估模型重建次季节融化范围高度和雪线观测能力，随后解释四个未来气候情景下 ELA 预测，评估阿拉斯加各子区域冰川脆弱性，量化驱动变化的关键过程。结果显示，2015—2100年阿拉斯加冰川将损失49 ± 17%至73 ± 26%的2015年质量，取决于未来气候；几乎所有 <1 km² 小冰川及 Brooks、Talkeetna 山脉到2100年将消失；中等增暖下负质量平衡在21世纪中叶达峰，而强增暖下持续增至2100年。研究人员将变化归因于温度升高导致融化增加，以及早冬降水由雪转雨减少积累并使冰面更早暴露。到2100年，剩余冰川中15%–60%无积累区，注定消失。该研究意义在于强调可靠高分辨率观测约束对冰川演化建模的价值，并提供可推广至全球冰川化区域的框架。

关键技术方法（不超过250字）：研究人员采用 Hugonnet et al. 的2000—2019年冰川尺度大地测量质量平衡（geodetic mass balance）观测，以及 Wells, Rounce, & Fahnestock 的2016年中—2024年次季节融化范围高度和瞬时雪线观测（覆盖99%面积≥2 km²冰川，n=3,018，约85%冰川化面积）。气候数据为 ERA5 日气温和降水（2000—2024）及 CMIP6 12个全球环流模型（general circulation model, GCM）与四个共享社会经济路径（Shared Socioeconomic Pathway, SSP）共48个实现。模型为 Python Glacier Evolution Model (PyGEM)，动力学用 Open Global Glacier Model (OGGM)。校准用马尔可夫链蒙特卡洛（Markov chain Monte Carlo, MCMC）分阶段校准温度偏差、降水因子和雪的度日因子；<2 km² 小冰川用 eXtreme Gradient Boosting (XGBoost) 预测参数。以积累区比率（accumulation area ratio, AAR）11年移动平均为0定义临界点。

研究结果：
4.1 现代模拟（Present-Day Simulations）：通过2000—2024日模拟评估模型对雪线和融化范围高度重建，并检验日月一致性。与约407,000雪线观测强一致（r²=0.82，偏差?80 m）；与约444,000融化范围高度观测一致性稍低（r²=0.63，偏差?25 m），早春4月误差较大，10月数据因 SAR 局限排除。日月模拟经日到月温度调整后累计质量平衡差<1%，说明可用日校准、月运行。使用参数均值、排除冰川动力学对现代质量平衡影响小，但大幅降低计算成本。
4.2 与质量平衡观测比较（Comparison With Mass Balance Observations）：通过 Gulkana、Wolverine、Lemon Creek 三个冰川的连续冰川学测量，模拟2000—2024年年和季节质量平衡。年均差分别为0.16、0.06、?0.13 m w.e.，年 MAE 为0.33—0.57 m w.e.；Lemon Creek 季节质量平衡偏差较大，可能因输入雪线观测噪声较大（中位不确定性17%海拔范围）。
4.3 未来模拟（Future Simulations）：通过月模拟到2100，阿拉斯加冰川损失49 ± 17%至73 ± 26%的2015年质量，取决于情景。Brooks 和 Talkeetna 即使 SSP1-2.6 也损失至少96%；St. Elias、Wrangell、Western Chugach、Western Alaska 等最抗性子区域 SSP1-2.6 损失36%—57%；SSP5-8.5 下12个子区域中9个损失>75%。SSP1-2.6 有67%（n=18,172）至 SSP5-8.5 有91%（n=24,720）冰川消失；<1 km² 消失79%—97%，1—10 km² 消失25%—76%，10—100 km² 消失0.8%—22%，>100 km² 几乎不全部消失。ELA 上升295—433 m，AAR 从40%降至24%（SSP1-2.6）或7%（SSP5-8.5），积累区减少40%—83%。
4.4 超过临界点的冰川（Glaciers Beyond Tipping Points）：通过11年移动平均 AAR 达到0定义临界点，2100年剩余冰川15%—60%无积累区，注定消失。多数冰川本世纪上半叶达临界点，因 <2 km² 小冰川占数量约89%但面积约12%。到2100年，72%—96%冰川按数量达临界点，但仅代表现代面积10%—34%。至少1,041条冰川可能最终达到新平衡，但2100年仍不平衡，平均损失2015年质量66%、面积47%。
4.5 模型对校准数据敏感性（Model Sensitivity to Calibration Data）：通过仅用2017—2021年融化和雪线子集校准，2000—2024年累计质量变化与全数据差3.2%；雪线 r²=0.79、偏差?94 m，融化范围高度 r²=0.60、偏差?49 m；校准期与未校准期准确度类似，说明结果对校准数据稳健，额外数据不预期进一步改善。
4.6 变化气候下冰川质量平衡组分（Components of Glacier Mass Balance in a Changing Climate）：通过评估融化、积累、降水、再冻结，夏季比融化几乎全部子区域增加；再冻结到2100年显著减少，尤其 Brooks Range；积累减少、降水增加，尤其晚夏，表明消融季延长，过渡季雪变雨。温度轨迹约2050年开始分歧，到2100年组分差异50%—100%。径流可下降至2015年值的100%，影响依赖冰川补给的水资源与生态。

讨论部分：5.1 与预测比较（Comparisons With Projections）：本预测与 GloGEM、OGGM、PyGEM 先前研究相似；Eastern Alaska range 少6%—7%质量损失；与185条最大阿拉斯加冰川研究一致。5.2 观测质量变化将继续加速（Observed Mass Changes Will Continue to Accelerate）：遥感观测显示质量损失率增加；预测到至少2050年继续增加，2045—2055年 SSP2-4.5 为?131 ± 58 Gt yr^?1，SSP5-8.5 为?153 ± 44 Gt yr^?1，比2015—2024年增加18%—51%。SSP2-4.5 世纪末率下降，SSP5-8.5 保持高。ELA 上升率3.7—5.4 m yr^?1。AAR 下降大于先前估计。5.3 更广泛区域影响（Implications on the Broader Region）：广泛冰川损失影响生态系统、水资源、自然灾害、旅游；冰川径流影响海岸水文、海洋、生物地球化学和生态；鲑鱼栖息地影响复杂；水电存在“峰值水”，主要水电区域约2050—2075年接近峰值径流；冰川湖和溃决洪水、旅游也受影响。

研究结论翻译：研究人员预测阿拉斯加冰川到2100年变化，评估区域和子区域尺度。模型框架包括瞬态融化范围和雪线观测及长期质量平衡数据。模型与观测雪线强一致（r²=0.82），与冰川学测量相比季节质量平衡准确（Gulkana 和 Wolverine 平均差在0.23 m w.e.内）。未来模拟显示广泛质量损失，49 ± 17%至73 ± 26% 到2100年相对于2015年。伴随 ELA 上升和积累区消失。预测区域 ELA 可上升达433 m，导致 SSP5-8.5 下仅7%当前冰面积在2100年 ELA 以上。2100年质量平衡评估显示增强融化与雪向雨降水转变，尤其晚夏，是质量损失主要贡献。2100年以后，估计2100年剩余冰川15%—60%过了临界点，注定消失，取决于未来气候。最终，发现强调高分辨率观测数据对大规模模型校准的重要性，特别是准确捕捉和代表季节过程。