《Journal of Hydrology: Regional Studies》：Mechanistic drivers and spatiotemporal thresholds of socioeconomic drought in the East River Basin, China: A causal-informed explainable AI approach
编辑推荐：
研究区域
东江流域（ERB）是华南特大城市的关键水源地。
研究焦点
为解决数据稀缺和极端赤字期间传统模型的预测崩溃问题，研究人员为分布式标准化社会经济干旱指数（DSSEDI）开发了一个综合诊断-预测框架。通过利用可获取的多源数据，克服了部门特定的监
研究区域
东江流域（ERB）是华南特大城市的关键水源地。
研究焦点
为解决数据稀缺和极端赤字期间传统模型的预测崩溃问题，研究人员为分布式标准化社会经济干旱指数（DSSEDI）开发了一个综合诊断-预测框架。通过利用可获取的多源数据，克服了部门特定的监测限制，实现了高保真DSSEDI预测（在极端赤字期间保持纳什-苏特克利夫效率为0.774）。通过将因果卷积Transformer与Shapley加法解释（SHAP）相结合，该框架通过双因果约束确保了物理一致性。至关重要的是，引入分位数回归损失函数（q=0.5）在捕获非高斯极端干旱峰值方面显著优于基于均方误差的模型。
新的水文见解
该分析揭示了东江流域潜在的社会经济干旱机制，该流域表现出以10至11个月水文气象滞后为特征的持续气候记忆。一个关键发现是识别出土壤水分悖论，即前期高土壤水分在城市热岛胁迫下通过增强蒸散发反而加速干旱发展。此外，该建模框架隐式诊断出新丰江水库的调度，并揭示了显著的社会经济拮抗作用，证明刚性工业需求是持续的干旱放大器。识别出不同的土壤水分拐点（城市区域750 mm，农业区域600 mm）可作为适应性管理的科学触发条件。这些发现建立了一种阈值驱动的治理范式，以增强城市化区域流域的干旱韧性。
干旱是气候变化加剧的复杂气候现象，通过持续降水赤字、增强蒸散发和土壤水分耗竭引发级联影响。干旱分为气象、农业、水文和社会经济四类，前三类已有标准化指数（如标准化降水指数SPI、标准化土壤水分指数SSMI、标准化径流指数SRI）监测，而社会经济干旱自2015年才引起关注。社会经济干旱是指区域供水系统无法满足人口和经济活动需求，导致水资源短缺。在东江流域（ERB），2004年和2009年的严重社会经济干旱使东莞工业产出增长下降了5.4%。然而，现有社会经济干旱指数均存在局限：多元标准化可靠性与韧性指数（MSRRI）依赖水库数据，社会经济干旱指数（SEDI）固定需求阈值且忽略基础设施调蓄，标准化供水需水指数（SWSDI）不能解析空间异质性。新提出的分布式标准化社会经济干旱指数（DSSEDI）通过时空降尺度解决了空间约束，但受限于高分辨率观测数据稀缺，传统统计方法难以捕捉非线性驱动机制。此外，传统机器学习模型如XGBoost等缺乏因果约束，且基于均方误差（MSE）优化，容易学习虚假相关并过度平滑极端干旱峰值。因此，亟需开发一个既能高保真预测又具有物理可解释性的综合框架。
针对上述问题，研究人员开发了一个集诊断与预测于一体的因果卷积Transformer（CC-Transformer）-Shapley加法解释（SHAP）框架。该框架采用双因果约束（因果卷积和掩蔽自注意力）和分位数回归损失函数（q=0.5），在数据稀缺条件下实现了DSSEDI的高分辨率合成与预测。在东江流域的验证表明，该模型在极端干旱期间保持了0.774的纳什-苏特克利夫效率（NSE），显著优于传统模型。研究揭示了土壤水分悖论、水库调节信号、社会经济拮抗等机制，并识别出关键土壤水分阈值。这些发现为城市化流域的干旱韧性提供了阈值驱动的治理范式。该研究成果发表于《Journal of Hydrology: Regional Studies》。
研究人员所用的主要关键技术方法包括：基于1992—2013年东江流域多源网格化数据（水文气象变量：降水、径流、土壤水分、温度、潜在蒸散发；社会经济变量：人口、GDP，来自国家青藏高原数据中心、TerraClimate、国家地球系统科学数据中心及广东省统计局），通过Pearson相关分析确定最优累积滞后（1—12个月）输入集；构建因果卷积Transformer模型，以因果卷积和掩蔽自注意力实现双因果约束，并使用分位数回归损失（q=0.5）替代MSE；采用空间分层划分（80%训练、10%验证、10%测试）检验空间外推；利用SHAP进行驱动贡献解耦，并与随机森林（RF）和XGBoost基准比较。DSSEDI目标值依据Meng等(2025)方法生成。
3.1 线性解耦与非线性建模的必要性
通过线性相关分析，研究人员发现径流、降水、土壤水分等水文气象变量与DSSEDI具有较强的线性相关（径流在1个月滞后时相关系数r=0.741），而社会经济变量（人口POP和GDP）与DSSEDI的线性相关很弱（|r|≤0.217）。这一“线性解耦”现象表明，人类活动主要通过复杂的非线性途径影响干旱演变，传统线性统计方法无法揭示其机制，因此必须采用非线性模型。
3.2 CC-Transformer模型的优越预测性能
将CC-Transformer与随机森林（RF）和XGBoost对比，在全样本评估中，CC-Transformer取得相关系数R=0.981、NSE=0.962，优于RF（R=0.816，NSE=0.659）和XGBoost（R=0.948，NSE=0.896）。在极端干旱情景（DSSEDI<-1.5）下，RF和XGBoost的NSE降为负值（-7.947和-4.263），发生性能崩溃；而CC-Transformer仍维持NSE=0.774和R=0.895。同时，CC-Transformer训练与测试性能差异极小（ΔR=0.004），表明其具备稳健的空间外推能力和泛化性。
3.3 掩蔽实验揭示功能非冗余性
通过在独立测试集上对单个和组合变量进行掩蔽（以均值替换），发现掩蔽降水（Pre1）导致最严重的预测崩塌（ΔNSE=14.94），超过径流和土壤水分；掩蔽GDP或POP也造成显著损失（ΔNSE>11.5），尽管它们线性相关微弱，说明社会经济变量通过非线性机制发挥独特作用。联合掩蔽中，降水-径流-土壤水分的联合损失（ΔNSE=22.86）小于三者个体损失之和，表明水文变量间存在部分补偿效应；而GDP-POP联合损失（ΔNSE=12.24）超过个体平均损失，体现社会经济因素的复合胁迫机制。值得注意的是，降水与径流高度相关（r=0.92），但二者掩蔽损失差异（ΔNSE=3.05）证明统计相关不等于功能等价，降水缺失无法由径流完全补偿。
3.4 CC-Transformer-SHAP动态驱动贡献解耦
通过CC-Transformer-SHAP分析，在全部情景和干旱情景下，土壤水分（SM1）均被识别为核心驱动因子，且贡献为负。这一负贡献源于东江流域广泛分布的砖红壤特有的“土壤水分悖论”：前期高土壤水分在城市化热岛效应增强的蒸散发需求下被快速消耗，反而加速干旱发展。GDP在干旱情景中持续显示负贡献，反映高GDP城市区域工业用水刚性需求成为干旱放大器；降水（Pre1）表现出缓解干旱的正贡献且旱季增强；径流（Runoff1）贡献稳定，但水库调度导致局部负效应；温度（Tem10）在干旱情景下出现反转，证实了热岛效应的放大；人口（POP）通过优先供水、承载力约束和节水效率三重协同产生稳定正贡献。进一步分析验证了模型能隐式捕捉新丰江水库调节信号，径流SHAP值的脉冲变化与水库调度阶段（蓄水、中性、放水）同步，证明框架无需直接运行数据即可诊断人为干预。
总结讨论部分，研究人员通过季节分离和空间聚集的诊断，揭示了驱动贡献的显著时空异质性。湿润季节（4—9月）SM1表现出阈值行为：低土壤水分对应正贡献（缓解干旱），高土壤水分对应负贡献（加剧干旱），且城市区域该悖论因热岛效应而放大；干燥季节（10月—次年3月）则出现反转，高土壤水分可缓解干旱。空间上，高度城市化的深圳在SM1约750 mm处出现贡献由正转负的拐点，而农业区河源的拐点约为600 mm，表明城市化将自然调节阈值抬高了约150 mm。基于此，研究提出了“双轨”适应策略：在城市区域利用海绵城市措施调节地表能量平衡和热岛效应；在农业区域依据600 mm阈值优化精准灌溉；同时优化水库在湿润季末（如9月）的蓄水调度以缓冲后续干旱，并对工业用水实施配额管理。讨论部分亦指出三点局限：社会经济数据时间粒度粗（GDP和人口为年度数据经均匀增长率降尺度）、部分生物物理反馈仍为半定量、训练数据截至2013年；未来方向包括将水库调节转化为显式参数、利用WRF耦合模式验证土壤水分-降水反馈、开展跨流域多尺度验证，以及引入随机降雨生成器扩展数据长度以支撑概率性干旱风险评估。
研究结论部分翻译如下：(1) CC-Transformer模型对社会经济干旱表现出优越的预测保真度和稳健的空间外推能力。通过集成因果卷积与掩蔽自注意力，并采用分位数回归损失函数优化，该框架确保时间因果一致性，克服了传统MSE模型固有的峰值过度平滑问题。关键的是，在极端干旱情景下，所提模型维持了0.774的高NSE，有效防止了传统树模型（RF和XGBoost）中出现的灾难性性能崩溃（负NSE）。这一稳健性凸显了该框架在高影响事件期间提供可靠早期预警的实际效用。(2)基于SHAP的可解释性揭示城市化东江流域深刻的社会经济拮抗作用。虽然人口主导的节水通常产生协同缓解效应，但工业用水需求刚性由于工业生产线的连续要求而成为持续的干旱放大器。值得注意的是，在没有直接水库运行输入的情况下，该框架成功隐式诊断出人为干预信号，这体现在径流SHAP值随新丰江水库调度强度的脉冲演变和相位同步，证明了模型解构复杂社会水文过程的能力。(3)社会经济干旱驱动因子的贡献表现出显著的时空异质性和非线性阈值突变。结果为精准定制的双轨适应策略提供了科学依据。在深圳等高城市化区域，识别的SM1拐点约750 mm标志着土壤水分从缓冲器转变为干旱加速器的关键转折，该阈值应作为海绵城市调节的科学触发条件；河源等农业区域表现出更高的缓冲能力，具有约600 mm的明显阈值，在此水平上必须优化精准灌溉以提高用水效率。除土壤水分反馈外，水治理必须利用水库调度的动态缓冲能力，具体而言，湿润季末最大化蓄水可以中和极端干旱信号并稳定由经济和人口增长驱动的基线需求。这种综合方法兼顾季节性反转和空间解耦，为全球城市化流域提高干旱韧性提供了可扩展的范式。
