基于事件的新疆山地-绿洲-荒漠系统多类型干旱传播三维追踪

《Journal of Hydrology: Regional Studies》：Event-based 3D tracking of multi-type drought propagation in Xinjiang’s Mountain-Oasis-Desert system

　　

研究区域：新疆（XJ）位于中国西北部，是典型干旱区，具有山地-绿洲-荒漠（MOD）耦合系统特征。研究重点：常规二维分析无法完整捕捉干旱MOD系统中多类型干旱的三维时空动态或完整传播链。为克服此局限，研究人员使用标准化降水蒸散指数（SPEI）、标准化径流指数（SRI）、标准化土壤水分指数（SSMI）和标准化生态干旱指数（SEDI）分别表征气象干旱（MD）、水文干旱（HD）、农业干旱（AD）和生态干旱（ED），并开发了一个基于事件的3D框架，将时空连通性聚类与有向无环图（DAG）网络分析相结合。该框架追踪干旱事件的生长、迁移和消散，同时重建1960–2024年新疆地区同步和滞后的跨类型传播联系。新的区域水文见解：共识别出231次3D干旱事件，包括55次MD、61次HD、83次AD和32次ED事件；总体上，它们表现出功能区聚集和基于质心的定向迁移，其轨迹大致与地形-河流-绿洲景观格局一致。DAG分析揭示了43个传播事件组，以线性级联模式为主（56%），平行和复杂模式分别占21%和23%。传播表现出路径依赖的滞后结构和不对称的跨系统耦合，而网络拓扑在2000年后变得更加复杂。此外，干旱传播敏感性表现出明显的三分空间格局：新疆北部以大气强迫为主（P_Atmo=24.1%），新疆南部以水文-农业反馈为主（P_Nexus≈20%），新疆东部以最高的生态响应敏感性为主（P_Eco=14%）。这些发现为区域干旱分析建立了基于事件的三维网络视角，并为干旱MOD系统中的非线性传播机制提供了新见解。

干旱作为全球变暖背景下最具破坏性的自然灾害之一，通常并非静态孤立事件，而是通过水文循环动态传播：降水亏缺引发气象干旱（MD），经土壤水分消耗和径流减少演变为农业干旱（AD）与水文干旱（HD），最终导致植被生产力下降和生态系统退化，即生态干旱（ED）。这种跨类型转化被称为干旱传播。新疆（XJ）地处中国西北内陆，是典型干旱区，具有由高山产流、绿洲用水、荒漠耗散构成的“山地-绿洲-荒漠”（MOD）耦合系统。其水源-利用-耗散结构为追踪冰冻圈、灌溉绿洲与荒漠生态系统的干旱信号传播提供了天然试验场。然而，既有研究多基于点尺度或二维时间序列，难以刻画干旱事件在三维时空中的连续演化；且多集中于“MD→HD”或“MD→AD”等二元关系，缺乏MD、HD、AD、ED全链条的系统分析；对非线性传播拓扑（平行、跳跃、反馈）及空间分异特征亦认识不足。针对上述三个关键缺口，研究人员提出了一个集成三维时空聚类与有向无环图（DAG）网络分析的事件级干旱传播识别框架，并应用于1960–2024年XJ的MOD系统，旨在回答三个问题：（1）不同类型干旱的三维时空演化特征是什么？（2）控制干旱传播的非线性拓扑结构和滞后机制是什么？（3）MOD系统如何塑造干旱传播的空间分异？该研究发表于《Journal of Hydrology: Regional Studies》。

在方法上，研究人员基于多源数据集（包括TPDC气象格网数据、ERA5-Land径流与土壤水分再分析数据、MODIS NDVI/GPP遥感产品），分别采用标准化降水蒸散指数（SPEI）、标准化径流指数（SRI）、标准化土壤水分指数（SSMI）和熵权欧氏距离（EWED）构建的标准化生态干旱指数（SEDI）表征MD、HD、AD与ED；利用基于连通性的三维时空聚类算法（含连通分量标记CCL和时空重叠追踪）识别三维干旱事件；通过同步与滞后匹配准则（含时空回溯投影STBP）建立事件间传播关系，构建DAG网络并利用深度优先搜索（DFS）提取传播链，计算传播概率（PP）、传播时间（PT_onset/PT_peak/PT_end）、传播位移矢量（PDV）和格网触发概率（P_trigger），进而识别大气强迫、水文-农业反馈与生态响应三类敏感性热点。

4.1 多类型干旱指数的时空变异性
通过Mann-Kendall检验和Sen斜率分析发现，1960–2024年SPEI呈“北湿南干”格局，新疆北部和东部（如天山北坡诸流域、吐鲁番-哈密盆地）显著变干；SSMI在大部分流域呈持续下降趋势；SRI趋势复杂，伊犁河流域和昆仑山区显著变湿，而塔里木河下游显著变干；2000–2024年SEDI高分辨率显示，显著变干区集中在塔里木河上中游农牧交错带和天山中东段，而天山北坡高寒草甸、伊犁河谷及塔里木河下游河岸林呈生态改善趋势。这表明MOD系统内不同地貌和生态功能单元对干旱的响应存在高度异质性。

4.2 三维干旱事件的时空动态
利用三维时空聚类共识别出231个事件（55个MD、61个HD、83个AD、32个ED）。MD事件平均历时4.0个月，平均面积3.8×10⁵ km²，质心集中于盆地内部，2000年后高强度MD频率增加；HD事件多集中高山区和源头区，但高烈度事件质心向塔里木河中下游聚集；AD事件呈显著绿洲聚集性，平均历时5.8个月，最长30个月，持久性最强；ED事件沿河道分布，平均历时4.4个月，空间扩展有限。迁移分析表明，MD迁移距离最长（最大约3800 km），方向多元；HD和AD以向东迁移为主，AD沿绿洲带形成明显东移走廊；ED主要在河道廊道附近局部振荡。典型案例（如2022年MD事件）验证了三维识别模型能较好重建实际干旱的生命史。

4.3 干旱传播特征与网络结构
通过DAG网络分析，在231个事件中，135个事件参与形成43个传播事件组，96个为孤立事件。传播概率矩阵显示，HD→AD的传播概率最高（39.3%），其次为MD→AD（36.4%）和MD→HD（25.5%）；AD→HD存在8.4%的反馈。时滞分析揭示MD驱动路径存在路径依赖差异：MD→HD和MD→ED的起始滞后较短（约1.8–1.9个月），MD→AD更长（平均2.8个月）；HD→AD起病快（1.0个月）但恢复慢（4.7个月）；AD→HD起始滞后最长且离散（6.5个月）。网络结构演化显示，1960–1999至2000–2024期间，MD→HD→AD级联链显著增强，ED作为末端受体加入，复杂网络模式占比从18.2%升至28.6%，表明拓扑结构趋于复杂。

4.4 干旱传播的空间异质性与敏感性
基于PDV和P_trigger分析，大气强迫路径（MD→HD/AD）表现出从山区向盆地汇聚的格局，新疆北部阿尔泰山区触发概率最高（P_Atmo平均24.1%）；水文-农业反馈路径（HD?AD）在塔里木河流域形成高风险绿洲廊道（P_Nexus约20%），尤其四条源流流域超过22%；生态响应路径（MD/HD/AD→ED）集中于下游河道和荒漠-绿洲交错带，吐鲁番-哈密盆地生态响应敏感性最高（P_Eco平均14.0%）。这一“北-中-南-东”分区模式与MOD系统的产流-用水-耗散功能分区高度对应。

讨论部分，研究人员基于典型事件拓扑重构，将传播模式归纳为级联（56%）、平行（21%）和复杂网络（23%）三类：级联模式反映水文缓冲与滞后；平行模式在极端气象强迫下出现“跳跃”和多受体同步触发；复杂网络模式则体现多源叠加与非线性反馈。结合MOD垂直分带，研究人员提出物理概念模型：上游高山区受冰冻圈调节，北疆降水-融雪混合补给使径流对气象波动更敏感，南疆冰川融水补给使系统更稳健；中游绿洲区由于强水文-土壤耦合和高蒸发需求，形成水文-农业反馈热点；下游荒漠区生态响应存在区域差异，南疆塔里木河下游因生态输水和深层地下水缓冲而敏感性较低，东疆吐鲁番-哈密盆地缺乏外源补给，生态脆弱性高。方法论上，该研究将干旱传播分析从时间序列依赖推进到事件对象与网络拓扑层面，但存在ERA5-Land再分析偏差、匹配阈值敏感性、自然与人为驱动解耦不充分等局限。

研究结论可概括为：（1）多类型三维干旱事件呈现功能区聚集和定向迁移，231个事件中MD跨区流动性最强，AD历时最长，ED沿河迁移受限；（2）传播以级联模式为主（56%），平行（21%）和复杂网络（23%）共存；HD→AD传播概率最高（39.3%），不同驱动路径具有不对称滞后；2000年后网络拓扑复杂化；（3）传播敏感性呈三区格局：北疆大气强迫主导（P_Atmo=24.1%），南疆水文-农业反馈主导（P_Nexus≈20%），东疆生态响应最敏感（P_Eco=14%）。该事件级三维网络框架为干旱MOD系统的非线性传播机制与分区风险管理提供了科学基础。