基于气象前兆的骤发性干旱渐进式预警系统：以中国东江流域为例

《Journal of Hydrology: Regional Studies》：A progressive early warning system for flash drought based on meteorological precursors: A case study of the Dongjiang River Basin, China

　　

研究区域： 中国东南部的东江流域，是珠江三角洲的关键水源区，属亚热带季风气候，干湿转换频繁。 研究焦点： 骤发性干旱（flash drought）因其发生迅速而威胁水安全和粮食安全，然而业务化预警仍主要依赖单一指标。本研究基于旬尺度（pentad-scale）ERA5-Land根区土壤湿度数据（1950–2025年），开发了一套渐进式预警系统。通过共线性诊断、时滞相关分析和随机森林重要性评估筛选前兆变量，继而采用受试者工作特征曲线（Receiver Operating Characteristic, ROC）分析，以ROC曲线下面积（Area Under the Curve, AUC）量化预测技能，并在两阶段回归—分类框架内确定预警阈值。 区域水文新认识： 研究识别出骤发性干旱发生前三个旬的五个前兆变量。提前两个旬的降水距平（Plag2，AUC = 0.74）与提前一个旬的降水距平（Plag1，AUC = 0.92）构成主要预警链，分别提前2个和1个旬发布L1和L2级预警。三个较弱的前兆变量（Plag3、ET/PETlag3、Tlag2；AUC均≈0.5）被整合为三分类支持条件评分（LOW/MEDIUM/HIGH）。该系统的误报率在训练期为6.3%、验证期为8.5%，显著低于单一指标基线（16.5%–18.6%）和逻辑回归基线（16.3%–22.3%），且具有更高的Heidke技能评分（Heidke Skill Score, HSS）（0.370 vs. 0.262–0.269）。结果在10组骤发性干旱定义参数组合下保持稳健，对2019年事件的验证表明预警及时、连贯且无误报。

论文解读

一、研究背景与问题

干旱是世界上最具破坏性的自然灾害之一，每年造成大量农作物减产和经济损失。传统上，干旱被视为一种缓慢演变的现象，其发展过程通常持续数月甚至数年。然而，近年来一种被称为骤发性干旱（flash drought）的特殊干旱类型日益受到国际科学界的关注。与传统缓慢发展的干旱不同，骤发性干旱以异常快速的加剧过程为特征，土壤湿度可在短短数周内从接近正常状态骤降至严重亏缺水平。这种突发性和快速性使得骤发性干旱对农业生产、水资源供给和生态系统的影响难以预测和缓解，留给灾害管理机构的响应时间极为有限。在全球变暖背景下，蒸发需求的增加与骤发性干旱频率和空间范围的显著上升相关联，这一趋势在中纬度农业区和季风影响区尤为突出，对粮食安全和区域水资源管理构成日益严重的威胁。

尽管骤发性干旱风险已得到越来越多的认识，但业务化早期预警系统仍不完善。现有指标如蒸散发胁迫指数（Evaporative Stress Index, ESI）、标准化蒸散发比率（Standardized Evaporative Stress Ratio, SESR）、短时间尺度标准化降水蒸散发指数（Standardized Precipitation Evapotranspiration Index, SPEI）以及土壤湿度百分位阈值等，主要用于骤发性干旱的回顾性识别而非早期预警。现有方法存在三个重要局限：第一，多数方法依赖单一指标阈值判断，难以捕捉复合前兆信号；第二，预警等级往往凭经验设定，阈值很少针对特定气候区域进行优化和验证；第三，骤发性干旱与非干旱样本之间的严重不平衡很少被明确处理，可能使预警性能评估产生偏差。

东江流域位于中国东南部珠江东部水系，是粤港澳大湾区的主要水源地，对香港、深圳、广州等大城市的供水具有战略意义。该流域属亚热带季风气候，降水季节变率大，需水量持续增长。尽管气候背景湿润，但该流域仍易受骤发性干旱影响。历史记录显示2004年和2019年曾发生严重骤发性干旱事件，其中2004年夏秋骤发性干旱导致干流流量骤降和水库蓄水持续消耗，威胁广东和香港的供水安全。然而，该流域迄今缺乏针对性的骤发性干旱早期预警框架，限制了水资源管理者及时采取有针对性的干旱应对措施的能力。

二、研究内容与主要结论

针对上述问题，该研究旨在基于气象前兆为东江流域构建一套渐进式骤发性干旱预警系统。具体目标包括：（1）利用共线性诊断、时滞相关分析和随机森林变量重要性评估识别关键前兆变量及其最优提前时间；（2）采用类别加权策略、受试者工作特征曲线（ROC）分析和约登指数（Youden's Index, YI）确定统计上稳健的预警阈值；（3）利用历史骤发性干旱事件构建并评估具有统计可解释性的渐进式预警系统。

研究得出以下主要结论：识别出跨越骤发性干旱发生前三个旬的五个前兆变量。提前两个旬的降水距平（P_lag2，AUC = 0.74）和提前一个旬的降水距平（P_lag1，AUC = 0.92）构成主要预警链，分别提前2个和1个旬发布L1和L2级预警。三个较弱的前兆变量（P_lag3、ET/PET_lag3、T_lag2；AUC均≈0.5）被整合为三分类支持条件评分（LOW/MEDIUM/HIGH）。该系统在训练期和验证期的误报率分别为6.3%和8.5%，显著低于单一指标基线（16.5%–18.6%）和逻辑回归基线（16.3%–22.3%），Heidke技能评分（HSS）达0.370，优于对比基线（0.262–0.269）。结果在10组骤发性干旱定义参数组合下保持稳健，对2019年事件的验证表明预警及时、连贯且无误报。该论文发表于《Journal of Hydrology: Regional Studies》。

三、主要技术方法

研究人员采用了以下关键技术方法：基于ERA5-Land再分析数据集（欧洲中期天气预报中心ECMWF发布，空间分辨率0.1°×0.1°，时间跨度1950–2025年），提取气象和土壤湿度变量并聚合至旬尺度（5天），通过厚度加权平均计算0–100cm根区土壤湿度；采用标准化距平处理消除单位、空间异质性和季节背景差异；通过Spearman秩相关矩阵进行共线性诊断（阈值|ρ|>0.70）；开展时滞相关分析（滞后0–5个旬）筛选与土壤湿度损失率显著相关的前兆变量；构建回归型随机森林模型（500棵树）进行变量重要性评估，综合置换重要性、消融ΔR2和Spearman|ρ|三种指标确定最终前兆变量；采用类别加权策略处理样本不平衡；通过ROC分析和约登指数确定预警阈值；在训练期（1950–2010）和独立验证期（2011–2025）严格分离的框架下完成模型构建与性能评估。研究区域为东江流域，经流域边界掩膜后保留318个有效网格单元。

四、研究结果

骤发性干旱事件识别与时空分布特征。 研究人员基于土壤湿度百分位定义骤发性干旱：土壤湿度百分位从40百分位以上在不超过4个旬内降至20百分位以下，且随后至少连续3个旬维持在20百分位以下。利用训练期数据建立了土壤湿度百分位分布，并应用于验证期以避免信息泄漏。

关键气象前兆变量识别。 通过共线性诊断、时滞相关分析和随机森林重要性评估的逐步筛选程序，研究人员从9个候选气象变量中识别出5个关键前兆变量。其中，提前两个旬的降水距平（P_lag2）和提前一个旬的降水距平（P_lag1）具有最强的预测能力（AUC分别为0.74和0.92），构成主要预警链；提前三个旬的降水距平（P_lag3）、蒸散发比（ET/PET_lag3）和提前两个旬的气温距平（T_lag2）预测能力较弱（AUC均≈0.5），作为支持条件使用。

预警阈值确定与渐进式预警系统构建。 研究人员采用类别加权策略处理骤发性干旱与非干旱样本的不平衡问题，通过ROC分析和约登指数确定统计上最优的预警阈值。构建了两阶段回归—分类框架：首先基于主要前兆变量（P_lag2和P_lag1）确定L1和L2两级预警，分别提前2个和1个旬发布；再将三个较弱前兆变量整合为三分类支持条件评分（LOW/MEDIUM/HIGH），用于增强或抑制预警信号。

预警系统性能评估。 该渐进式预警系统在训练期误报率为6.3%，验证期为8.5%，显著低于单一指标基线（16.5%–18.6%）和逻辑回归基线（16.3%–22.3%）。Heidke技能评分（HSS）为0.370，优于对比基线的0.262–0.269。系统在10组骤发性干旱定义参数组合下均保持稳健性能。

典型事件验证。 研究人员将构建的预警系统应用于2019年骤发性干旱事件，结果表明系统能够及时、连贯地发布预警信号，且未出现误报，验证了系统的实际应用能力。

五、讨论与结论

讨论部分。 研究人员讨论了预警结构的物理意义：降水亏缺是骤发性干旱发生的主要驱动因素，提前1–2个旬的降水距平能够有效指示土壤湿度的快速下降；蒸散发比和气温距平作为辅助条件，反映了大气蒸发需求对土壤水分消耗的贡献。研究还检验了系统对骤发性干旱定义参数组合的敏感性，结果表明预警性能在10组参数组合下保持稳健，说明系统对定义参数的选择不敏感。研究的不确定性主要来源于ERA5-Land再分析数据本身的误差、土壤湿度百分位分布基于训练期估计可能引入的偏差，以及验证期样本量有限等问题。未来工作可考虑引入动态预报产品以延长有效预警提前时间，并将系统推广至其他气候区域进行适用性检验。

结论部分翻译。 该研究以东江流域为研究区，基于1950–2025年ERA5-Land根区土壤湿度数据，开发了一套基于气象前兆的渐进式骤发性干旱预警系统。主要结论如下：（1）通过共线性诊断、时滞相关分析和随机森林重要性评估，识别出五个跨越骤发性干旱发生前三个旬的关键气象前兆变量；（2）提前两个旬的降水距平（P_lag2，AUC = 0.74）和提前一个旬的降水距平（P_lag1，AUC = 0.92）构成主要预警链，可分别提前2个和1个旬发布L1和L2级预警；（3）三个较弱前兆变量（P_lag3、ET/PET_lag3、T_lag2）被整合为三分类支持条件评分，用于增强预警的可靠性；（4）该渐进式预警系统的误报率（训练期6.3%，验证期8.5%）显著低于单一指标和逻辑回归基线，Heidke技能评分（HSS）达0.370，优于对比方法；（5）系统在10组骤发性干旱定义参数组合下保持稳健，对2019年事件的验证表明预警及时、连贯且无误报。该研究提出的渐进式预警框架为东江流域骤发性干旱业务化预警提供了技术支撑，并为其他气候区域构建类似预警系统提供了可推广的方法论参考。