2000–2023年羊卓雍错流域水位变化驱动因素分析

《Hydrology》：Analysis of the Drivers of Water-Level Changes in the Yamdrok Lake Basin During 2000–2023

　　

摘要：羊卓雍错（Yamdrok Lake）是位于青藏高原南部的一个典型封闭内陆湖泊，对气候变化和人类活动高度敏感。尽管已有研究报告了湖水位下降趋势，但自2005年以来急剧下降的机制仍存在争议。特别是，诸如抽水蓄能电站运行等人类活动的干扰效应尚未得到严格表征，且在多种因素综合影响下的定量归因仍然有限。本研究旨在系统识别2000–2023年间羊卓雍错水位变化的驱动因素，并利用Shapley R2分解量化气候变化和人类活动对水位波动的相对统计解释力。研究人员整合了24年的水文、气象、遥感和水电运行数据，采用累积距平分析、多元线性回归和Shapley R2分解，建立了一个包含降雨量、蒸发量、温度、冰川融水量以及由年发电量代表的水电运行强度的归因模型。通过Shapley R2分解，针对不同时期（包括水电运行期和非运行期）量化了这些因素对年际水位变化（ΔH）的相对统计解释力，该解释力代表各因子在回归框架内对ΔH方差的解释比例，而非绝对的物理体积贡献。结果表明，2000–2023年间羊卓雍错年均水位显著下降，累计下降5.1米，且2005年后下降加速。主导控制因素从早期的降雨-蒸发模式转变为以气温升高为主导的系统性水分亏缺模式。气温以0.03 °C yr-1（p < 0.05）的速率显著上升，构成了湖水位长期下降的根本驱动力。在水电运行期（2000–2014），五个因子共同解释了年际水位变化方差的84.51%，其中水电运行强度（32.32%）、气温（27.98%）和冰川融水量（26.54%）是统计解释力最大的三个贡献因子。在非运行期（2015–2023），无论是对2022–2023年缺失的冰川融水数据进行外推还是将其排除在分析之外，气温始终是最重要的单个贡献因子，占解释方差的34.3%至49.5%。变暖通过两条途径影响湖水位：它既直接增强湖面蒸发，又同时促进流域内冰川持续退缩。2000年至2021年间，冰川面积减少了约30.5 km2，从而削弱了冰川融水补给的长期韧性。因此，湖泊系统已从动态平衡转向一种水分损失超过水分输入的持续状态。总体而言，羊卓雍错的水位变化是渐进式变暖引起的水分亏缺与叠加的水电运行干扰共同作用的结果。通过将观测延长至2023年并纳入动态冰川面积和水电运行指标，本研究阐明了主导驱动因素的时序转变，修正了先前将水位下降主要归因于降雨减少的解释，并为青藏高原湖泊水资源安全评估和气候变化适应提供了科学依据。

羊卓雍错水位变化驱动机制研究解读

一、研究背景与科学问题

羊卓雍错位于青藏高原南部，是典型的封闭内陆湖泊，对气候变化与人类活动高度敏感。青藏高原多数内陆湖泊近几十年来呈现扩张趋势，但羊卓雍错却表现出与众不同的水位下降态势。此前的研究虽已报告其水位下降，但关于2005年后急剧下降的机制存在显著分歧：部分研究认为降雨减少是主因，另一些则强调人类活动的影响。然而，现有工作普遍缺乏对抽水蓄能电站运行等人类活动干扰效应的严格刻画，也未能在多因素综合作用下进行定量归因。此外，冰川融水作为重要补给源，其动态变化及其与气候变暖的耦合关系在以往研究中关注不足。因此，系统识别2000–2023年间羊卓雍错水位变化的驱动因素，并量化气候变化与人类活动的相对贡献，对于理解青藏高原湖泊水文响应机制及保障区域水资源安全具有重要意义。该论文发表于《Hydrology》期刊。

二、研究方法概述

研究人员集成24年（2000–2023年）的水文、气象、遥感及水电运行数据，主要采用了三项关键技术方法：首先，应用累积距平分析法初步揭示各单因子与水位变化的阶段性对应关系；其次，构建多元线性回归模型，以年际水位变化（ΔH）为因变量，将降雨量、蒸发量、气温、冰川融水量以及年发电量（代表水电运行强度）作为自变量；最后，运用Shapley R²分解方法，定量评估各因子对ΔH方差的相对统计解释力。气象数据来源于浪卡子国家气象站，水位与蒸发数据来自白地水文站，冰川面积基于遥感影像提取，水电运行数据来自羊卓雍错抽水蓄能电站运行统计。

三、主要研究结果

3.1 气候变化特征

研究发现，2000–2023年间羊卓雍错流域年均气温呈显著上升趋势（p < 0.05），升温速率为0.03 °C/年，2013–2023年间增温尤为明显。年降雨量呈不显著下降趋势（p > 0.05），平均每年减少约1.63 mm。值得注意的是，年蒸发量呈现显著下降趋势（p < 0.01），降幅约为7.31 mm/年，这与同期气温显著上升形成了“蒸发悖论”。研究人员指出，这一现象在青藏高原已被多项研究证实，主要归因于风速降低、水汽压增加和太阳辐射减弱等因素抵消了升温带来的蒸发需求增加。

3.2 水文变化特征

羊卓雍错年均水位在2000–2023年间呈现显著的长期下降趋势。2000–2004年水位保持较高水平，2004年达到峰值后持续大幅下降，2008年起低于多年平均值。至2023年，水位降至研究期内最低点，累计下降达5.1米。尽管个别年份（如2016–2018年）出现小幅回升，但整体下降趋势未变。

3.3 水文气象相关性分析

累积距平分析显示，降雨与水位在2000–2004年和2019–2022年期间呈现较一致的同相位变化，但在2005–2015年间出现明显背离，表明降雨单独无法完全解释水位波动。蒸发与水位在2012年前大致呈反向趋势，但此后水位在蒸发减少的情况下仍持续下降。气温累积距平与水位变化呈现一定程度的反向对应关系。Pearson相关分析进一步表明，年际水位变化（ΔH）与气温呈显著中等负相关（r = -0.494, p = 0.017），而与降雨和蒸发的相关性均不显著。月尺度分析揭示了水位的季节性循环：1–6月持续下降，7月起迅速回升，9–10月达峰值后回落，这一模式与季风降雨的季节分布高度吻合，且水位峰值滞后降雨峰值约两个月。

3.4 人为驱动因素分析

羊卓雍错抽水蓄能电站总装机容量112.5 MW。2000–2014年间，年发电量先升后降，2006年达历史峰值6.55×10⁸ kW·h，此后急剧下滑，2014年趋近于零。研究人员强调，由于缺乏年抽水量和净补水量的数据，年发电量仅用作表征电站运行扰动强度的间接指标，不能直接等同于净耗水量。

3.5 冰川融水分析

基于度日因子法（DDF）估算的冰川融水量显示，2000–2021年间流域年冰川融水量在1.75×10⁹至2.12×10⁹ m³之间波动，呈微弱下降趋势（约-4.2×10⁶ m³/年）。冰川面积从2000年到2021年减少了约30.5 km²。这表明冰川系统已越过临界点——融化面积的缩减效应超过了单位面积融化速率的增加效应，导致后期（2018年后）总融水量下降。独立验证显示，2000–2016年间流域冰体积损失约0.72 km³。对于2022–2023年缺失的冰川融水数据，研究人员采用CMFD 2.0气象强迫数据集经偏差校正后估算，结果显示这两年的融水量（分别为3.88×10⁹和3.70×10⁹ m³）远高于历史平均水平。

3.6 影响因素贡献分析

在水电运行期（2000–2014年），多元回归模型的决定系数R²为0.8451，即五个因子共同解释了年际水位变化方差的84.51%。Shapley R²分解结果显示，各因子按贡献排序为：水电运行强度（32.32%）> 气温（27.98%）> 冰川融水量（26.54%）> 蒸发量（7.03%）> 降雨量（6.13%）。在非运行期（2015–2023年），无论是否纳入2022–2023年估算的冰川融水数据，气温始终是贡献最大的因子（34.3%–49.5%），但由于样本量极小（n=7或9），模型未通过显著性检验，结果应视为探索性发现。

四、讨论与结论

讨论部分指出，本研究存在若干局限性：未能分离季节性积雪融水和冻土退化对径流过程的影响；缺乏抽水蓄能电站的年抽水量和净输水量数据，无法进行基于水量平衡的物理归因；非运行期样本量过小限制了统计推断的可靠性。研究人员强调，“蒸发悖论”并不削弱气温作为水位下降根本驱动力的结论——变暖通过直接增强蒸发需求和间接加速冰川退缩两条途径影响湖水位。Shapley R²分解得到的贡献率反映的是统计解释力而非物理体积贡献，这一点在解读水电运行强度变量的贡献时需要特别留意。

研究结论翻译如下：

1. 1.
   气候变暖引起的水热失衡是羊卓雍错水位长期下降的根本驱动力。 2000–2023年间，湖水位累计下降约5.1米，2005年后衰退加速。气温以0.03 °C yr^-1（p < 0.05）的速率显著上升，并与年际水位变化呈显著负相关（r = -0.494, p = 0.017）。变暖一方面直接加剧了湖面蒸发和流域蒸散的水分损失，另一方面驱动了流域内冰川持续萎缩（2000–2021年间冰川面积减少约30.5 km²）。冰川融水补给从2006–2012年的暂时增强阶段转变为2018年后的减弱趋势，随后因2022–2023年的异常温暖条件而急剧增加。降雨与水位的关系在2005年后显著脱钩，表明降雨不再是控制长期水位变化的主导因子。从根本上说，羊卓雍错水位下降反映了流域水热系统在气候变暖下从动态平衡向净水分亏缺状态的转变，而非单纯由降雨减少所致的过程。
2. 2.
   抽水蓄能电站运行是影响水位波动的额外扰动，而非简单的线性耗水项。 在2000–2014年的发电期，由年发电量代表的水电站运行强度对年际水位变化的统计解释贡献占32.32%，是该时期模型中解释力最大的单一变量。该百分比代表在多元回归框架内，发电量年际变异对ΔH方差统计解释的比例，而非水电运行导致的净湖水损失比例。虽然2015年后大规模发电停止，湖水位并未停止下降，但水位波动的变异性显著减小，进一步表明电站运行主要是放大了年际水位变异性，而非长期下降趋势的独立驱动力。
3. 3.
   本研究与先前将降雨减少认定为主因的研究之间的差异，主要源于研究时段和驱动因素识别框架的不同。 先前研究通常基于1970–2000年左右的数据，当时气温上升相对缓慢，降雨年际变异的短期调节作用仍可掩盖变暖的渐进累积效应，因此降雨被识别为最明显的驱动因子。相比之下，本研究覆盖2000–2023年，捕捉到了2005年后特别是2013年后变暖加剧、冰川退缩加速、降雨调节作用减弱的临界转型期，使得气温的长期、先前不那么显著的作用变得更加突出。本研究中获得的贡献排序——水电运行强度 > 气温 > 冰川融水 > 蒸发 > 降雨——应被理解为这些因子在2000–2014年特定时期内对年际水位变异性的统计解释力排序，而非它们对总水位下降的绝对物理贡献排序。气候变暖通过气温升高和冰川退缩为水位持续下降提供了背景条件，而电站运行则在这一背景趋势上施加了显著的年内和年际扰动，两者共同塑造了2000–2014年羊卓雍错水位的演变过程。