水库调节改变了澜沧江-湄公河流域气象–水文干旱的传播

《Agricultural Water Management》：Reservoir regulation alters meteorological–hydrological drought propagation in the Lancang–Mekong River Basin

　　

理解梯级水库如何改变气象干旱向水文干旱的传播对于跨境水管理至关重要。以澜沧江-湄公河流域（1960–2016年）为研究对象，研究人员整合了收敛交叉映射（Convergent Cross Mapping, CCM）、Copula-贝叶斯框架（Copula-Bayesian framework）和主导性分析（dominance analysis），以量化传播动态、动态阈值及潜在驱动机制。结果表明，水库调节在空间上重新配置了水文干旱风险。在旱季，中游和下游子流域经历显著的干旱衰减（旱季标准化径流指数（SRI）增加高达0.85），而上游干旱却矛盾地加剧（高达63%）。此外，气象干旱传播为水文干旱所需时间在全流域范围内延长，表明水文响应对同期气象强迫的依赖性减弱，同时下游存在明显的季节性不对称。关键的是，水库运行对应于触发下游水文干旱所需的本地气象干旱阈值降低（中游子流域最高达0.6个标准化降水指数（SPI）单位），并与中等条件下这种跨类型传播的概率降低相关。根据主导性分析，水库的存在改变了既定的产流机制，显著增强了本地降水对下游径流驱动的相对贡献（从13.5%增加到32.4%）。这些发现凸显了水利基础设施如何改变干旱风险的空间格局，为适应性水库运行提供了坚实的科学基础。

**论文解读：水库调节对澜沧江-湄公河流域气象–水文干旱传播的影响**

**一、研究背景与问题**

干旱是一种由严重缺水引起的自然灾害，对人类社会、财产和生态系统构成重大威胁。不同于突发性灾害，干旱发展缓慢、持续时间长且影响范围广，不仅导致农业减产、水资源短缺、粮食不安全和经济损失，还对生态系统稳定性和社会经济的可持续发展产生长期影响。在全球气候变化加速和社会经济快速发展的背景下，干旱已成为严峻挑战。根据美国气象学会（AMS）的分类，干旱通常分为气象干旱、水文干旱、农业干旱和社会经济干旱四类。气象干旱（以降水亏缺为标志）是最基本的类型，其他类型干旱通常源于气象干旱但存在时间延迟，这一水分亏缺在不同干旱类型间传递的过程被称为干旱传播。研究干旱传播对于理解干旱动态演化过程和影响机制具有重要的理论和实践价值。

早期研究主要依赖相关性分析，近年来发展了更复杂的统计方法，如Copula函数、条件分布模型和因果分析方法（如收敛交叉映射，CCM）。然而，确定干旱传播阈值时，传统研究常采用固定阈值，缺乏坚实物理基础且对极端干旱表现不佳。此外，现有研究多聚焦于自然驱动的干旱传播，对人类活动（尤其是水库调节）的关注不足。水库运行可显著改变气象-水文干旱间的关联，但对其如何系统性改变上下游传播关系、季节性特征和空间响应模式仍缺乏系统理解。澜沧江-湄公河流域（LMRB）是亚洲重要的跨境河流，水电资源丰富，截至2022年已建有多座水库，且数量预计持续增加。这种密集的水库开发将显著改变自然干旱传播过程，但针对该流域的定量研究仍然有限。

**二、研究目的与意义**

为填补上述研究空白，本研究选取LMRB为研究区，整合标准化降水指数（SPI）和标准化径流指数（SRI），结合游程理论（Run theory）、CCM、Copula-贝叶斯框架（Copula-Bayesian framework）、受试者工作特征曲线（ROC曲线）以及基于Shapley值分解的主导性分析（dominance analysis），系统量化梯级水库调节对气象-水文干旱传播的影响机制。研究成果将为跨境水资源的适应性配置、针对性干旱风险评估和预警系统开发提供坚实的科学依据。该论文发表在《Agricultural Water Management》。

**三、关键技术方法**

本研究采用的主要关键技术方法包括：（1）基于游程理论识别干旱事件并提取持续时间和强度特征；（2）应用收敛交叉映射（CCM）量化气象干旱向水文干旱的传播时间，通过分析不同时间尺度SPI与SRI之间的交叉映射技能确定传播滞后；（3）结合Copula函数与贝叶斯公式，构建联合概率分布计算干旱传播条件概率，并引入ROC曲线和混淆矩阵客观确定最优传播阈值；（4）采用基于Shapley值分解的主导性分析，定量解构上游入流（上游SRI）与本地降水（本地SPI）对下游径流变异性的相对贡献。数据来源包括：英国气候研究单元（CRU）的全球降水量时间序列产品（1960–2016年），流域6个主要水文站的日径流观测数据（1960–2016年），以及全球水库与坝数据库（GRanD）的水库属性数据。研究区基于HydroSHEDS Level 6子流域划分和6个水文站位置，将流域分为7个子流域（剔除河口区子Ⅶ），最终选取子流域Ⅰ–Ⅵ进行分析，其中Ⅰ–Ⅱ为上游、Ⅲ–Ⅳ为中游、Ⅴ–Ⅵ为下游。水库期划分以各子流域内最具调控作用的水库投运年份为界。

**四、研究结果**

**3.1 水文干旱风险的时空重构**

基于1960–2016年各子流域SRI计算，揭示了水文干旱的时空演变。旱季中下游子流域（尤其Ⅳ、Ⅴ、Ⅵ）在库后旱季SRI显著上升（Ⅳ上升0.85），而上游Ⅰ–Ⅲ变化小于0.15。游程理论分析表明，库后中下游子流域干旱持续时间和强度显著降低（Ⅳ持续时间减少29.3%，Ⅳ和Ⅵ强度分别从3.54降至2.71、3.89降至2.42），而上游Ⅰ–Ⅲ干旱强度增加（Ⅰ从3.21升至4.82）。结论：水库建设缓解了中下游水文干旱，但加剧了上游干旱，呈现空间异质性。

**3.2 传播时间延长与季节性不对称**

CCM分析表明，所有子流域的干旱传播时间均显著延长：上游Ⅰ、Ⅱ、Ⅳ、Ⅴ从3、4、7、3个月增至8个月；Ⅲ和Ⅵ从2、3个月增至6个月。季节性分析显示，库后上游和中游湿季传播时间从平均4.1个月增至6.3个月（Ⅳ延长103%）；下游Ⅴ和Ⅵ湿季传播时间缩短（分别减少35%和34%），而旱季显著延长（分别增加207%和93%），产生明显下游季节性不对称。此外，SPI–SRI因果耦合强度发生空间重组：上游Ⅰ、Ⅲ显著减弱（Δρ?=-0.25、-0.09），中下游Ⅳ、Ⅴ、Ⅵ显著增强（Δρ?=0.13、0.17、0.27）。结论：水库运行改变了干旱传播的时间尺度及其空间和季节分布。

**3.3 干旱传播概率与阈值的变化**

基于Copula-贝叶斯框架计算旱季气象干旱触发水文干旱的条件概率P(HD|MD)。上游Ⅱ概率显著增加（ΔP?=0.20，p<0.001），中下游Ⅳ、Ⅴ、Ⅵ概率显著降低（ΔP?=-0.17、-0.14、-0.17，均p≤0.001）。ROC曲线最优SPI阈值分析表明：旱季中下游中等和严重干旱阈值加深或消失（Ⅳ和Ⅵ的中等干旱阈值从-0.5、-0.7变为无可靠阈值NT；Ⅴ严重阈值从-1.4变为NT），表明缓冲效应强；上游Ⅰ严重和极端阈值升高（ΔSPI=0.8和0.9），Ⅱ出现新阈值，表明传播更易发生。湿季中下游极端干旱阈值均消失（NT）。结论：水库运行并非全流域抑制干旱传播，而是下游降低触发概率和阈值（缓冲），上游增加。

**五、讨论与结论**

讨论部分指出，水库运行典型地削峰补枯，增加流域蓄水量和水文记忆，使得气象亏缺必须先消耗额外蓄水才能表现为径流亏缺，从而延长传播时间并削弱降水-径流耦合。基于Shapley值分解的主导性分析显示：上游子流域Ⅱ在自然期上游入流贡献76.6%，库后增至86.5%，本地降水贡献降至13.5%，表明人工控制主导；下游子流域Ⅵ自然期模型解释力低（R2=0.20），库后R2升至0.54，上游入流相对贡献从86.5%降至67.6%，本地降水（SPI）从13.5%增至32.4%，说明水库拦截了上游缺水信号并增强了本地降水驱动作用，这解释了旱季传播时间延长和阈值降低的物理机制。研究存在数据偏差、模型结构、方法特定不确定性等限制，未来需耦合分布式水文过程模拟和动态运行策略。

研究结论翻译如下：
本研究评估了梯级水库调节如何改变澜沧江-湄公河流域（LMRB）气象-水文干旱的传播。通过结合收敛交叉映射（CCM）、Copula-贝叶斯建模和主导性分析，研究人员量化了水库运行对干旱传播的影响及其内在机制。主要结论如下：
1. 水库建设显著改变了水文干旱的时空格局。旱季，中下游子流域经历明显的干旱缓解，表现为干旱持续时间和强度降低。相反，特定上游子流域的水文干旱特征加剧。
2. 水库运行普遍延长了全流域的干旱传播时间。上游和中游地区在湿季和旱季均出现传播延迟。然而，下游地区表现出明显的季节性差异：旱季传播时间显著延长，而湿季传播时间缩短。
3. 水库调节了干旱触发条件。在受调控的下游地区，触发水文干旱所需的标准化降水指数（SPI）阈值变得更负，表明水库有效缓冲了中等气象水分亏缺。相反，上游子流域阈值升高，反映了其对严重和极端干旱的更高脆弱性。
4. 梯级水库重塑了径流驱动因素的相对权重。上游水文日益受水库泄流主导（上游入流方差贡献达86.5%）。下游，由于干流水库的缓冲效应，虽然上游入流仍是主要驱动因素（从86.5%降至67.6%），但本地降水的相对贡献显著增加（从13.5%增至32.4%）。这一机制转变解释了两个观测到的下游现象：本地气象-水文干旱传播延迟，以及对上游缺水的纵向传递敏感性降低。