利用自适应三维聚类方法刻画中国小时极端降水时空动态特征

《Journal of Hydrology: Regional Studies》：Characterizing spatiotemporal dynamics of hourly extreme precipitation in China with an adaptive 3D clustering approach

　　

研究区域：中国，其范围涵盖东亚季风区及其西北内陆过渡带，在这些区域中，极端降水事件（Extreme Precipitation Events, EPEs）表现出强烈的空间异质性，并构成重大的水文灾害风险。
研究重点：本研究调查了1985–2024年间中国EPEs的时空演变特征与传播模式。核心目标在于超越基于站点、二维或日尺度的描述，在连续三维空间中重建完整的降水系统轨迹。为此，研究开发了一种基于三维密度聚类（3D-DBSCAN）的自适应三维时空聚类框架，并将其应用于欧洲中期天气预报中心第五代再分析数据集（ERA5）的小时降水数据。将极端降水格点映射到标准化的经度—纬度—时间空间中，同时利用k距离参数优化和物理约束过滤来减少主观参数选择，并避免对快速移动或间歇性事件的人为分割。随后，使用降水加权质心量化事件的移动距离和方向，并对典型传播模式进行分类。
区域新水文见解：该框架识别出独立的EPEs，并揭示了五种主导传播模式。中国的EPEs表现出明显的西北向移动偏好，而向南事件则稀少且空间范围有限。西北向事件，尤其是超长距离事件，表现出最大水文极端性，具有更长的持续时间、更大的影响面积和更高的累积量级。这些轨迹表明极端降水活动存在区域性的“西移北跳”趋势，该趋势与东亚夏季风、西太平洋副热带高压和中纬度西风带相关联。

论文解读：中国小时极端降水时空动态的自适应三维聚类研究

一、研究背景与问题

在全球自然灾害中，极端气候事件是造成气象灾害人员伤亡和经济损失的首要因素。政府间气候变化专门委员会（IPCC）第六次评估报告（AR6）明确指出，全球大部分地区的极端降水事件（Extreme Precipitation Events, EPEs）正变得更加频繁和强烈，对人类社会经济与生态系统构成严重威胁。中国是受EPEs影响最严重的国家之一，其降水事件空间分布极不均匀，夏季EPEs主要集中在长江中下游、华南和华北东部等地。气候变化进一步放大了短时EPEs的极端性与突发性，例如2021年7月河南郑州极端降水事件造成398人死亡或失踪，直接经济损失达1200.6亿元人民币。当前关于EPEs的研究主要侧重于频率、影响面积、持续时间和强度等单个统计特征，且长期高时间分辨率观测数据匮乏，以往分析多依赖日尺度或更粗分辨率数据。虽然日尺度数据能刻画EPEs的总体特征，但会平滑掉降水日内过程的动态细节，难以表征短时、快速移动或间歇性演变降水系统的连续三维时空演化，从而限制了对EPEs整体动力学和物理机制的深入理解。近年来，三维识别方法逐渐应用于气候极端事件研究，但现有方法仍存在关键参数缺乏客观标准、主观性强、对快速移动事件易产生错误分割等问题。因此，亟需一种灵活且自适应的聚类框架，以在不依赖严格空间重叠或人为连接规则的前提下，刻画连续演变的极端降水系统。

二、研究内容与主要方法

针对上述问题，研究人员利用欧洲中期天气预报中心（ECMWF）提供的第五代全球气候再分析数据集（ERA5）小时降水数据，时间范围为1985年1月1日至2024年12月31日，空间范围70°E–138.5°E、0°N–58°N，分辨率为0.25°×0.25°。研究提出了自适应三维时空聚类框架，该框架基于三维密度聚类算法（3D-DBSCAN），包含以下主要关键技术方法：（1）极端降水阈值确定：将小时降水≥0.1 mm定义为有效降水小时，采用基于百分位数的方法，选取99.9th百分位数作为各格点的极端降水阈值，以严格捕捉可能引发严重灾害的高影响EPEs。（2）时空映射与标准化：将超过阈值的格点映射到标准化经度—纬度—时间三维坐标系，引入时空缩放因子τ，实现空间与时间坐标的无量纲统一。（3）自适应参数优化：采用k距离自适应方法，通过计算每个降水点到其第k近邻的欧氏距离，动态确定候选邻域半径（Eps）和最小样本数（Min_samples），并通过聚类稳定性判断最优参数组合。（4）三维聚类与物理约束过滤：在标准化三维空间中进行密度聚类，核心点邻域内满足最小样本数则形成时空连通簇，代表独立EPEs；随后应用持续时间（T_min=6 h）、空间范围（N_min=16个格点，约10,000 km2）和地理有效性约束，排除短时孤立或局地小尺度噪声事件。（5）移动特征刻画：采用降水加权质心量化事件移动距离和方向，利用大圆距离公式和方位角计算；结合K-means聚类和内部效度指标（Davies-Bouldin指数、轮廓系数、惯性下降率）确定最优移动模式数量。（6）轨迹密度估计：通过对轨迹线性插值生成密集点云，使用自适应核密度估计（KDE）和Scott准则确定带宽，并将连续密度场栅格化为0.1°×0.1°分辨率，采用二维高斯卷积滤波平滑伪影。研究人员依据上述框架开展系统敏感性实验，确定最优参数配置：稳定性阈值S=5、时空缩放因子τ=4。

三、研究结果

（1）1985–2024年中国EPEs总体特征
通过自适应三维聚类框架，共识别出799次独立EPEs。统计显示，中国EPEs以中等持续时间和平均强度为主，大规模、高量级事件发生频率较低。持续时间高度集中于24–120 h区间，占85.9%，其中48–72 h最频繁；超长事件（>240 h）仅占0.6%。平均强度呈单峰分布，87.7%的事件集中在13–16 mm/h，峰值位于14–15 mm/h。影响面积分布跨度大（2.34×10⁴至1.01×10⁶ km²），分布相对均匀。累积量级呈明显右偏分布，63.6%的事件低于1×10⁷ km²·mm/h，极端事件占比仅5.3%。

（2）EPEs空间分布特征
EPEs空间分布呈显著异质性，从东南沿海向西北内陆递减。高频热点区主要位于华南沿海及邻近岛屿，累积频次达60–80次；次级活跃区覆盖长江中下游和淮河流域（30–50次），并向北延伸至山东和辽东半岛；而西北内陆地区为极端少发区（<5次），并存在从四川盆地至华北平原的显著过渡带。

（3）EPEs时空移动特征
移动方向表现出显著的不对称性，西—西北（WNW）和西北（NW）方向为主导移动方向，共占所有事件的26.0%，而正南向事件仅9次（1.1%）。移动距离在250–500 km区间出现频次峰值（221次），超1500 km的事件有93次。方向与距离存在系统耦合：西至西北方向在长距离事件中占绝对优势，而东向和南向事件移动距离普遍低于750 km。基于移动方向和距离的联合聚类，划分为五种移动模式：类型1（西北向超长距离）、类型2（西北向中长距离）、类型3（西南向短距离）、类型4（东南向短距离）和类型5（东北向中距离）。类型1平均移动距离1933.0 km，平均持续时间132.2 h，累积量级高达3.35×10⁷ km²·mm/h，破坏性最强；类型3和类型4向南移动，持续时间短、影响范围有限；类型5东北向事件内部变率大，平均强度较高。

（4）EPEs移动模式空间分布与典型案例
西北向事件（类型1和2）主要集中于华南沿海省份，合计占所有EPEs的40.4%，且具有最高累积量级。类型1的典型事件（2018年9月11–17日）受台风“山竹”影响，从开阔洋面生成后向西北移动3143.6 km，持续145 h，影响广东、广西等五个省份，累积量级5.51×10⁷ km²·mm·h^?1。类型5集中于华北内陆，频率最高（27.4%），其代表事件（2021年7月19–22日）质心从河南向河北东北移动308.4 km，平均强度达到14.89 mm/h。南向事件类型3主要局限于雷州半岛和海南岛，类型4集中于江淮流域，后者影响面积大但平均强度低。此外，框架成功捕捉到快速移动事件（如2024年10月31日–11月2日事件，39小时移动1674.4 km）以及2012年7月21–22日北京复杂阶段性演变事件，证明该方法能保持快速、短时及复杂演化降水系统的完整性。

四、讨论与结论

敏感性分析表明，时空缩放因子τ=4在平衡时间与空间权重上最优，τ过小会过度分割快速移动系统，τ过大会错误合并独立事件；稳定性阈值S=5提供最保守稳定的聚类结果；物理约束阈值T_min=6 h和N_min=16格点与中尺度对流系统（MCSs）和Mesoscale-β尺度（20–200 km）相匹配，参数邻域内结果稳健。移动模式的合理性归因于多尺度环流系统协同调控：类型1和2的长距离西北向移动主要受东亚夏季风（EASM）和西太平洋副热带高压（WPSH）西伸北抬的共同驱动，超长距离事件常与台风或热带低压相伴；类型5受中纬度西风系统和太行山、燕山地形强迫影响；类型3和4则反映局地热力过程和海陆风、地形抬升等中小尺度过程。中国EPEs雨带“西移北跳”现象与“西部增湿、北方降水恢复”的年代际转型趋势一致，其成因包括WPSH增强、面积扩大和西伸，以及北大西洋多年代际振荡（AMO）正位相和太平洋年代际振荡（PDO）负位相协同作用，加上温室气体增暖背景下热带对流对海温异常的敏感性增加。中高纬度环流调整，如极地放大效应导致西风急流减弱弯曲，有利于准静止大振幅罗斯贝波列形成，与异常西伸的WPSH在西北—华北交汇，配合青藏高原东北侧地形强迫，产生强辐合上升区，引导降水系统西北向移动并增加西北内陆EPEs。数据与方法局限性方面，与CMORPH和MSWEP产品的交叉比较显示，区域尺度EPEs主要结论较为稳健，但绝对事件数、局部影响面积和轨迹细节依赖于数据产品和分辨率。与26连通性方法的对比表明，自适应3D-DBSCAN框架在事件连续性与分离度之间实现更好平衡，尤其适合快速移动和间歇性系统。研究结论可概括为：（1）提出了结合k距离自适应优化与时空标准化映射的自适应三维时空聚类方法，最优参数为S=5、τ=4，显著降低主观参数依赖。（2）识别出中国799次EPEs，持续时间和平均强度呈集中分布，影响面积分布较宽，累积量级右偏。（3）揭示五种移动类型，发现显著“西北向”移动偏好和南向规避特征，其中西北向超长距离类型灾害风险最高，其动力学归因于WPSH异常西伸与EASM共同作用，以及中纬度西风带动力引导，建立了通往内陆的稳定水汽输送通道。该研究弥补了日尺度分析在描述降水系统动态过程方面的不足，为中国极端降水移动规律提供了新视角，对揭示水文响应机制和支撑精细防灾减灾具有科学意义与应用价值。