GPC: 地图约束模拟揭示1980-2018年青藏高原多年冻土面积消失近18万平方公里
在全球变暖的背景下,青藏高原这个“亚洲水塔”和全球最大低纬度多年冻土区,正经历着剧烈而复杂的地下变化。冻土退化不仅影响高原水文循环、生态系统,还可能通过碳释放影响全球气候。然而,传统陆面过程模型(LSM)因参数校准困难,模拟结果往往差异巨大。
以往冻土模拟研究多采用“站点校准—区域外推”策略:用单个钻孔或活动层站点的土壤温度数据优化参数,再推至整个高原。但这种方法极易受站点代表性限制。在干旱站点优化出的参数,放到湿润站点模拟就变差。
Permalab最新的工作提出地图约束校准(map-based calibration)方法。以2010年高质量的青藏高原冻土分布图(Cao et al. 2023)为“全域标尺”,在超过13,000个敏感过渡网格(transitional cells) 上,利用贝叶斯优化框架同时调整26个关键参数(包括植被反照率、土壤孔隙度、底边界温度等),避开了传统二元地图带来的严重异差同效问题,同时保证参数在空间上的普适性。
独立验证显示,地图约束模拟(CaliMAP)在83个活动层站点和291个钻孔上全面领先。而传统站点校准(CaliTGL、CaliTSH)和默认参数(CaliFree)在多数区域表现明显逊色,部分情况下站点校准表现甚至不如默认参数。
利用优化参数集合(13成员)驱动1979-2018年连续模拟,研究给出高可信度的QTP冻土演变全景:
- 面积大幅收缩:1980年1.262 ± 0.048 × 10⁶ km² → 2018年1.085 ± 0.049 × 10⁶ km²,总减少约17.7万 km²,平均速率44.4 × 10³ km²/10a,处于以往同类模拟研究(36.2–92.0 × 10³ km²/10a)的较低端,更为保守可靠。
- 热力学指标普遍退化:活动层加深0.04 m/10a,冻土顶升温0.10 °C/10a,零年振幅深度上升0.52 m/10a。
- 空间异质性显著:连续冻土区边缘、冈底斯山脉、三江源东部退化最剧烈;而西北部(喀喇昆仑区域)1990年代末出现明显大气冷却(与Karakoram异常一致),导致局部新生冻土、冻土面积短暂回升,显著减缓了整个高原的退化趋势。
这是方法上首次将高质量分布图成功应用于高维参数化过程模型(30个可调参数),为冻土建模社区提供了可复制的地图约束范式。提供了更准确的历史冻土退化情况,对青藏高原基础设施规划、碳源汇评估、水资源安全具有直接参考价值。(图 文: 嵇海龙, 南卓铜)

引用:
Ji H, Wu X, Zhao S, Nan Z. Constrained simulation of permafrost thermal changes from 1980 to 2018 on the Qinghai-Tibet Plateau. Global and Planetary Change. 2026. https://doi.org/10.1016/j.gloplacha.2026.105542.