青藏高原中部高程依赖性融雪贡献的突变

《Journal of Hydrology: Regional Studies》：Abrupt shift in elevation-dependent snowmelt contribution in the Central Tibetan Plateau

　　

研究人员的研究区域是青藏高原长江源区。研究重点在于探究高山流域中融雪驱动径流的高程依赖性是否会在近期气候变暖下发生突变。研究人员应用了一个分布式水文模型，该模型受径流观测和卫星衍生的积雪数据的约束，旨在研究过去四十年间高程-径流-融雪关系的变化。针对该区域的新水文见解是：研究结果揭示了大约在2000年发生了一次突然的水文重组，当时春季冻结层高度（H₀）达到约4750米海拔，接近统计得出的阈值（H_crit = 4750±150米）。超过阈值后，高程对总径流的方差解释减少，但对融雪贡献的方差解释增加。全流域的融雪贡献也从长期增加转变为持续下降。这一转变与高海拔变暖加速、雪线上升、降雪比例下降、蒸散发增强以及模拟活动层增厚和基流增加相吻合。气候逆转实验表明，在较冷条件下，转变前的状态可以恢复，但存在多年的滞后。这些发现表明，变暖可引发高山流域中高程依赖的水文阈值转变，并强调冻结层高度是检测冰冻圈影响源区新兴变化的有效诊断指标。

**论文解读文章**

**研究背景与问题**

融雪是高寒山区水文循环的核心组成部分，维持着全球众多雪域流域的径流、生态系统及下游供水。山区作为关键的水塔，其水文过程受温度、降雪比例及高程依赖的积雪存储与融化强烈控制，因而对气候变化高度敏感。在气候变暖背景下，冷季降水形态、积雪积累与融化时机的改变已开始影响山区源区河流的径流季节性及量级，对水安全和生态系统稳定构成重要影响。青藏高原（TP）作为全球最大最高的山地水塔之一，其变暖速率（1979-2020年约0.34°C/十年）超过同期全球均值，并伴随广泛的冰冻圈与水文变化，包括积雪覆盖改变、降水相态转变及径流 regime 变化。现有研究表明，这些水文气候与生态水文变化具有强烈的高程结构：极端降水增幅在高海拔显著放大（Song等, 2026），变暖信号在4300米以下呈现高程依赖性（Xu等, 2020），植被绿度与物候响应亦随高程变化。然而，多数研究聚焦于长期趋势，对高程-径流关系是否会出现类似阈值的突变性重组，以及融雪贡献对径流的高程依赖性如何因持续变暖而骤变，尚缺乏直接定量证据。一个可能的机制是春季冻结层高度（H₀，即0°C等温线平均海拔）的上升：当H₀超过流域特定临界高程时，变暖驱动的雪线迁移与降水相态转变可能大幅重组径流生成的空间结构，削弱高程与总径流的历史耦合，同时增强对融雪贡献的高程控制，并级联影响积雪储量、融水释放、蒸散发及冻土相关水文过程。因此，该研究旨在探究高山流域中高程依赖性融雪贡献是否存在突变、春季冻结层高度是否为气候触发因子、以及多过程耦合如何驱动并维持这一转变。

**研究内容与结论**

研究人员以青藏高原长江源区高海拔流域为研究对象，结合物理机制驱动的分布式可变下渗能力（VIC）陆面水文模型、高分辨率气象强迫数据集及多源遥感与地面观测，开展了系统研究。主要结论包括：① 在高程-径流与高程-融雪贡献关系中检测到约2000年发生的显著突变（Altitudinal Snowmelt Contribution Shift, ASCS）；② 春季冻结层高度（H₀）在2000年越过统计临界阈值（H_crit≈4750±150米），且ΔR²（融雪贡献高程解释力与径流高程解释力之差）与H₀间呈现明确的两阶段阈值行为；③ 高海拔变暖加速、雪线上升、降雪比例下降、活动层增厚及基流增加、蒸散发增强等多过程协同作用，共同驱动并维持了转变；④ 气候逆转实验表明该转变在冷胁迫下可逆，但恢复存在6-7年的滞后，表明流域具有水文记忆。

**技术方法**

本研究主要采用以下关键技术方法：(1) 分布式VIC陆面水文模型（可变下渗能力模型），空间分辨率0.1°×0.1°，模拟地表水文过程及融雪积累与融化（VIC标准雪模块），并启用冻土选项；采用SCE-UA全局优化算法进行参数率定。(2) 高分辨率气象强迫数据：第三极气象强迫数据集（TPMFD），空间分辨率1/30°（≈0.033°），时间分辨率1小时，涵盖1979-2022年，已在高山复杂地形验证。(3) 多源观测验证数据：长江源直门达（1979-1997及2014-2022）与奔子栏（2008-2020）日径流观测用于模型率定与独立验证；MODIS逐日无云500米积雪覆盖产品（2000-2020）用于独立验证模拟积雪面积。(4) 突变检测与阈值分析：采用Pettitt检验与标准正态同质性检验（SNHT）检测高程-径流与高程-融雪贡献关系的断点；通过bootstrap重采样估计H_crit的置信区间；构建气候逆转实验（以1980-1990年冷时期序列替换末期强迫）测试可逆性及恢复时间尺度。

**研究结果**

**3.1 模型性能**

研究人员通过率定期（1980-1997，NSE=0.73, R2=0.75, PBIAS=-2.4）与验证期（直门达2014-2022：NSE=0.74, R2=0.76, PBIAS=-3.6；奔子栏2008-2020：NSE=0.78, R2=0.79, PBIAS=-2.1）验证了VIC模型对径流与积雪动态的合理再现能力。与MODIS积雪覆盖对比（2000-2020）显示整体精度0.767，加权F1值0.755。针对1998-2013年观测空白的伪缺验证表明，重建径流的中位NSE为0.774，中位绝对年径流误差13.8%，重建误差可量化约束。模型可靠再现了各高程带积雪面积的高程梯度与季节性演变（MODIS比对相关系数0.981-0.996），为后续分析奠定基础。

**3.2 高程融雪贡献转变（ASCS）**

研究人员发现，约1999-2000年起高程与径流的关系急剧减弱，同时高程与融雪贡献的关系显著增强。以ΔR2（R2(SCR, mm) – R2(Runoff)）为主断点指标，Pettitt检验与SNHT均识别出2000年为显著断点，蒙特卡洛扰动（5000次）中位断点仍为2000年，95%区间1998-2004年，83.8%-84.3%识别在1999-2001年，证实断点稳健。1980-1999年间，年径流与高程中等相关（R2=0.45），融雪贡献（SCR）与高程相关为0.27（mm）和0.51（%）；2000年后，径流R2降至0.12，而SCR的R2升至0.58（mm）和0.71（%），表明高程-径流-融雪关系发生突变成组（ASCS）。全流域SCR也从前一时期的上升趋势（+0.78×10⁹ m³/十年）转变为下降趋势（-0.31×10⁹ m³/十年）。

**3.3 春季冻结层高度阈值（H_crit）**

ASCS的出现与春季冻结层高度（H₀，4-6月0°C等温线平均海拔）越过约4750米临界值密切相关。H₀呈显著上升趋势（+9.3 m/yr，p<0.001），从1980年代初约4600米升至2020年代初4900米以上。将年ΔR2对H₀进行LOESS拟合与分段回归，显示明确的两阶段模式：H?<4750米时ΔR2近零，超过后快速增大（分段回归R2=0.513，p=0.00069，斜率变化显著p=0.0158）。Bootstrap分析得H_crit=4750±150米（95%CI:4600-4900米）。这表明一旦H?超过流域特定阈值，高程对年径流的解释力快速减弱，而对融雪贡献的高程控制则加强，支持水文阈值响应解释。

**3.4 转变的驱动因子**

研究人员通过多指标协同分析揭示ASCS受一套相互耦合的冰冻圈与水文过程驱动。

**（1）高海拔变暖加速与雪线上升**
高程分层分析表明，变暖速率在3000-4500米仅从+0.20增至+0.24°C/十年，但在4500-5000米从+0.22增至+0.50°C/十年，5000-5800米从+0.25增至+0.55°C/十年，最强加速在4700-4800米（接近H_crit）。相应，年最大雪线从1980-1999年约5080米升至2000-2022年约5180米，上升速率从+2.9增至+6.1 m/yr（加速112%）；积雪持续时间从46天减至35天，缩短速率从-0.31增至-1.02天/年。

**（2）降水相态转变与积雪减少**
年固态降水比例（SnowFrac）从1980-1999年约25%降至2000-2022年约20%，下降速率从-0.3×10^-3/年加速5倍至-1.5×10^-3/年。季节分解显示春、秋季变化最显著：3月降雪比例从60%降至47%，4月从45%降至23%，10月从30%降至18%，11月从30%降至17%；1-2月仍以雪为主（>85%），反映雪季缩短、春季消融提前、秋季积雪推迟。

**（3）活动层增厚与基流增加**
模拟活动层厚度（ALT）在2000年前缓慢增加（+0.29 cm/十年），平均约2.01米；2000年后加速至+2.52 cm/十年，2022年约2.26米。ALT增厚主要集中在4500米以上，最大增幅（+14.8 cm/十年）出现在5500-6000米带。基流从1980-1999年平均约45毫米升至2000-2022年约91毫米，且ALT与基流呈显著正相关（Kendall τ=0.48，p=0.001），表明正反馈路径。

**（4）蒸散发增强与反馈机制**
实际蒸散发（ET）在2000年前略降（近480毫米），2000年后快速增加（+21.4 mm/十年，p<0.001），至2022年增加约50毫米（相对增量约10%）。Mann-Kendall检验支持此趋势。ET增加与H?>H_crit及水文重组同步，尤其在低、中海拔地区增强的大气损失削弱了可用于径流的水量，高海拔地区ET增加与积雪减少共同强化了融雪贡献的高程对比。

**3.5 气候逆转实验**

研究人员通过以1980-1990年冷时期序列替换末期11年强迫，考察系统可逆性。结果显示，所有指标（R2(Runoff)、R2(SCR, mm)、ΔR2）在强迫切换后7年内重返历史目标带（R2(Runoff)=0.349，R2(SCR, mm)=0.265，ΔR2=-0.050）。指数拟合给出e-folding时间τ约2.5-2.8年，90%恢复时间t₉₀约6-7年，表明ASCS状态在冷条件下降可逆，但恢复存在滞后。

**讨论与结论**

讨论部分强调，2000年长江源区的突变反映了气候变暖下融雪供给的阈值性重组，而非简单线性响应。春季冻结层高度（H₀）上升是主要气候触发因子，其跨越统计临界高程带（H_crit≈4750±150米）后，ΔR²-H₀关系呈现明确两阶段行为，表明转变非渐进过程。多过程协同驱动机制包括：高海拔变暖加速（4500米以上变暖速率增至+0.50-0.55°C/十年，最强加速在H_crit附近）、雪线上升（年最大雪线从5080米升至5180米，上升速率加快）、降雪比例加速下降（SnowFrac从25%降至20%，春季月下降幅度达13-22个百分点）、活动层增厚（ALT从2.01米增至2.26米，与基流正相关）及蒸散发增强（ET从略降转为+21.4 mm/十年）。这些变化共同反映了雪储量压缩、降水相态转变、水道径流路径重组与大气水损失增加，而非单一过程作用。气候逆转实验表明转变可逆但存在6-7年滞后，可能反映积雪重建、季节缓冲恢复及地下储水路径调整所需时间。研究还讨论了不确定性：TPMFD强迫误差、1998-2013年径流观测空缺、VIC模型未显式模拟永久冻土退化与冰川融化过程等，但多源验证的一致性支持转变本身的稳健性。H_crit应视为统计阈值带，非通用值。该研究提供了可转移的诊断框架，指出春季H₀、雪线位置、降雪比例等可作为预测水塔流域突变逼近的早期指示指标。

研究结论翻译如下：本研究显示，持续变暖已重组青藏高原长江源区融雪-径流过程的高程依赖性。通过将分布式水文模型与长期气象强迫、遥感产品及径流观测相结合，研究人员识别出约2000年高程-径流-融雪关系的突然转变。该转变与春季冻结层高度（H₀）越过统计临界高程带（H_crit≈4750±150米）一致，表明流域经历了阈值性水文重组而非线性变暖响应。检测到的转变与一系列冰冻圈与水文变化相关联：高海拔变暖抬升了春季H₀与雪线位置，降雪比例下降，季节性积雪储量被压缩，总融雪贡献减少；同时，蒸散发增加及地下径流分配变化（反映为活动层增厚与基流增强）可能强化了历史高程-径流关系的减弱。气候逆转实验进一步表明，此 regime 转变在冷胁迫下物理可逆，但恢复延迟（约6-7年），意味着流域具有可测量的水文记忆。这些发现直接相关于长江源区，但亦暗示类似的高程触发水文重组可能发生在其他具有可比冰冻圈环境的高山流域。具体阈值与水文响应将取决于盆地特定控制因素（地形、气候 regime、冰川影响与冻土条件）。在此背景下，监测春季H₀、雪线位置、降雪比例及相关冰冻圈边界可为识别持续变暖下逼近可用水突变的流域提供实用基础。