1980-2023年岗日嘎布山冰川物质平衡变化及其驱动因素

《Advances in Climate Change Research》：Changes in glacier mass balance and their drivers in the Kangri Karpo Mountains from 1980 to 2023

　　

位于青藏高原东南部的岗日嘎布山冰川正经历加速物质损失，然而对其物质平衡演变和气候驱动因素的长期过程基评估仍然有限。在此，研究人员应用分布式COSIPY（COupled Snowpack and Ice surface energy and mass balance model in PYthon）能量质量平衡模型，以高分辨率气象数据为驱动，重建了1980–2023年岗日嘎布山561条冰川的物质平衡。模型性能通过实地花杆观测和大地测量物质平衡估计进行了评估。与全部561条冰川的大地测量估计的区域对比显示，2000–2010年和2010–2020年相关系数分别为0.76和0.81。模拟结果显示区域持续物质损失，至2023年累积物质平衡达到?11.9 ± 3.5 m w.e.（米水当量）。年均物质平衡从1980–2000年的?0.1 ± 0.1 m w.e. 每年下降至2001–2023年的?0.4 ± 0.1 m w.e. 每年，表明2000年后冰川损耗明显加剧。能量平衡诊断表明，这种加剧主要由增强的短波能量吸收驱动：两个时期之间的平均净短波辐射从32.3 W/m²增至37.6 W/m²，熔化能量从17.7 W/m²增至20.6 W/m²。敏感性实验进一步识别出气温是主导气候控制因素。在雅弄冰川上，+3 K升温使物质平衡减少了2300 mm，而?3 K降温使其增加了1400 mm。相比之下，降水减少50%使物质平衡减少了1200 mm，而增加50%仅产生了900 mm的增益。这些结果表明，区域冰川物质损失主要由大气变暖驱动，并通过雪雨相变、地表反照率降低和低海拔消融加剧而放大，为理解青藏高原东南部最近的冰川变化以及预测未来冰川演变和水文响应提供了过程基基础。

论文解读：岗日嘎布山冰川物质平衡变化及其驱动因素（1980–2023年）

一、研究背景与问题

冰川作为气候变化的重要指示器，在全球变暖背景下正加速退缩，导致海平面上升、自然灾害频发并影响区域水资源。青藏高原东南部（SETP）的岗日嘎布山位于印度季风强烈影响的区域，是海洋型冰川发育中心之一，其对温度、降水和能量通量变化极为敏感。然而，现有研究多依赖遥感手段获取冰川面积、厚度和物质变化的净结果，难以解析季节性物质积累与消融过程、物质平衡组分分配以及驱动年际变率的气候控制因素。此外，大地测量方法只能提供离散时段内的阶段性估计，缺乏连续年序列来揭示长期趋势和加速转折点。尽管已有研究应用温度指数模型（如OGGM）进行区域重建，但这类方法简化了地表能量交换，无法归因于具体物理过程。因此，针对该区域，特别是岗日嘎布山的长期连续冰川物质平衡演变及关键气候驱动因素的相对贡献，尚缺乏系统的过程基评估。

二、研究内容与意义

研究人员应用分布式COSIPY能量质量平衡模型，结合高分辨率气象再分析数据和冰川编目数据，重建了1980–2023年岗日嘎布山561条冰川的连续物质平衡序列。通过花杆观测和大地测量估计验证模型性能，并量化了时空变率、高程依赖性及气候敏感性。研究发现，区域冰川物质损失持续且加速，2000年后年均损失从?0.1±0.1 m w.e. 每年加剧至?0.4±0.1 m w.e. 每年，累积物质平衡至2023年达?11.9±3.5 m w.e.。能量平衡诊断表明，净短波辐射增强是消融加剧的主要近因，而气温通过融-反照率反馈和雪雨相变等机制成为主导气候控制因素。该研究首次在季风主导区域实现了区域尺度的过程基物质平衡重建，为理解近期冰川变化和预测未来演变及水文响应提供了基础。论文发表在《Advances in Climate Change Research》。

三、主要关键技术方法

研究采用的关键技术方法包括：（1）高分辨率气象强迫数据：以第三极气象强迫数据集（TPMFD，1/30°逐时）为主，辅以高亚洲精细再分析（HAR v2）云量数据，通过垂直梯度分析和双线性插值降尺度至冰川尺度，并结合COSIPY辐射校正模块处理短波辐射地形效应。（2）分布式COSIPY能量质量平衡模型：耦合地表能量平衡与多层雪/冰方案，显式解析短波穿透、液态水渗透与滞留、重冻结释热等过程，模拟表面物质平衡与内部物质平衡。（3）模型校准与验证：基于帕隆94号冰川花杆消融观测（2010–2017）和雅弄冰川大地测量物质平衡数据（2000–2020）进行迭代参数优化，利用其余560条冰川独立大地测量估计进行区域验证。样本来源于青藏高原东南部岗日嘎布山区域，无额外队列来源。

四、研究结果

**4.1 模型校准与验证** 通过帕隆94号冰川花杆观测校准，模拟与实测物质平衡的皮尔逊相关系数为0.858，均方根误差（RMSE）为0.591 m w.e.；雅弄冰川模拟与大地测量估计在2000–2020年吻合良好。区域尺度上，全部561条冰川的模拟与大地测量物质平衡相关系数在2000–2010年为0.761，2010–2020年为0.812，表明校准后的COSIPY模型能可靠重现从单点到区域的物质平衡变率。

**4.2 年物质平衡** 模拟显示，1980–2000年年均物质平衡为?0.1±0.1 m w.e.每年，2001–2023年下降至?0.4±0.1 m w.e.每年，表明2000年后冰川损失显著加剧。累积物质平衡持续下降，至2023年达到?11.9±3.5 m w.e.。

**4.3 物质平衡组分与季节循环** 表面融化是主导消融组分，2000年后显著增强；降雪是主要物质来源但不足以抵消融化损失。季节循环显示，降雪主要集中在3–5月，消融峰值出现在6–8月（夏季），两时期对比表明夏季融化加剧是近期加速损失的主要驱动力，冬季积累变化较小。

**4.4 能量通量** 净短波辐射（Net SW）是最大正能量输入，净长波辐射（Net LW）为持续辐射汇。消融季节（6–9月）净短波辐射和净辐射分别上升至67.1和64.5 W/m²，驱动强烈消融（66.9 W/m²）。后期（2001–2023）较前期（1980–2000）净短波辐射从32.3增至37.6 W/m²，净辐射从17.2增至22.2 W/m²，熔化能量从17.7增至20.6 W/m²，表明增强的短波吸收是后期区域融化加剧的主要驱动。

**4.5 雅弄冰川高程依赖性物质平衡变化与平衡线高度（ELA）位移** 雅弄冰川上，低海拔区域消融加剧最为显著，平衡线高度（ELA）从4805 m上升至4894 m。由于冰川面积-高程分布（hypsometry）集中在现代ELA附近，这一微弱上升导致积累区面积比（accumulation-area ratio）缩小，消融区扩大。组分分析显示，低-中海拔融化增加最强，多数高程带降雪减少，低海拔融化增强是物质平衡负值化的主要成因。

五、讨论总结

不确定性分析表明，模型参数中雨雪阈值和新鲜雪反照率是控制物质平衡模拟的主导因素，冰反照率影响次之，而粗糙度参数影响极小。云量数据来源（HAR v2）可能引起长波辐射估计的不确定性，但不影响区域物质损失加速的结论。碎屑覆盖仅占冰川面积的4.7%，且厚度较薄，其忽略可能导致局部偏差，但不会改变气候强迫为主导控制因素的基本结论。

通过相关分析和单因子敏感性实验，研究人员发现：年尺度上，冰川物质平衡与降水（r=0.81）和反照率（r=0.74）显著正相关，与气温（r=?0.53）和短波辐射（r=?0.57）显著负相关；月尺度上，冬季和春季以降雪驱动积累为主，暖季以气温和辐射驱动消融为主。敏感性实验显示，雅弄冰川物质平衡对温度响应呈强非线性，+3 K升温导致物质损失?2300 mm，而?3 K降温仅增加+1400 mm；降水变化的影响不对称且有限，50%降水减少造成?1200 mm损失，50%增加仅带来+900 mm增益；短波辐射增加20%引起?600 mm损失。这些结果表明，气温是主导控制因素，降水增加难以抵消变暖引起的消融增强。

区域机制分析将高程依赖性变化与大气变暖直接联系：变暖抬升冻结高度，使更多季风降水由雪转为雨，减少积累并降低反照率，导致低海拔地表吸收更多短波辐射，形成正融-反照率反馈，并导致ELA上升和积累区收缩。这些机制共同解释了2000年后的加速损失，且由于雅弄冰川具有代表性，可推广至整个岗日嘎布区域。

方法论进展方面，本研究首次在季风主导区域应用分布式能量质量平衡模型对所有561条冰川进行连续重建，提供了三个传统方法无法实现的优势：连续年序列（1980–2023）精确定位加速转折点；显式能量通量归因识别净短波辐射为近因；过程级组分分配和单因子敏感性分离出气温通过反馈机制成为主导控制。

六、研究结论翻译

本研究应用分布式COSIPY模型，以高分辨率气象强迫为驱动，重建了1980–2023年岗日嘎布山561条冰川的物质平衡，这是该季风主导区域首次区域性、基于物理过程的物质平衡重建。通过花杆观测和全部561条冰川独立大地测量估计的验证表明，校准后的模型能够重现单冰川和区域尺度的物质平衡变率。主要发现总结如下：

（1）1980–2023年期间，岗日嘎布山冰川经历了持续且加速的物质损失。加速在2000年后明显：年均物质平衡从1980–2000年的?0.1 ± 0.1 m w.e. 每年下降至2001–2023年的?0.4 ± 0.1 m w.e. 每年，累积物质平衡至2023年达到?11.9 ± 0.35 m w.e.。

（2）增强的表面消融是近期退化的主要原因，在能量平衡层面由增加的净短波辐射驱动。表面消融主导冰川损耗，净短波辐射上升是后期更强消融的主要驱动力——并由雪盖收缩和降水由雪转为雨造成的地表反照率下降所放大。

（3）气温是主导气候控制因素。敏感性实验表明，变暖产生的响应强度远大于降水或短波辐射的等效扰动，且单独增加降水无法抵消变暖引起的物质损失。变暖还通过加剧低海拔消融、平衡线高度上升和积累区收缩进一步放大了物质损失。通过整合区域模型重建、独立验证和过程级归因，本研究加深了对青藏高原东南部海洋型冰川如何响应持续变暖的理解，并为其未来演变和水文影响的预测提供了过程基基础。