观潮 · 模式预报与气象 AI 观察
云南气象干旱向1米深土层垂向传播：最深层要等4.64个月，触发阈值降到-1.70，传播概率只剩40%
> 气象干旱和土壤干旱不同步：降水亏缺要先过入渗、蒸散这些陆面过程的筛子才进得了土壤，而且不同深度的土壤不是一回事。这篇用 1981—2020 年 ChinaMet 气象数据和 ERA5-Land 三层土壤水分（0–7、7–28、28–100 cm），在云南省把这条垂向链条量化成了三个数：SPEI 向三层SSI 的平均传播时间 1.41、2.08、4.64 个月，平均传播概率 51%、48%、40%，有效格点的平均 SPEI 触发阈值 −1.32、−1.46、−1.70——越深的土，等得越久、越难被触发、要求的气象干旱越强。同期还有一个不对称：SPEI-1 长期趋势不显著（通过显著性检验的格点仅 20%），三层 SSI 却都显著下降（显著格点 92%、87%、75%）。方法上用相关分析定传播时间、AIC 优选的 Copula 定传播概率、0.60 条件概率标准定触发阈值，再用地理探测器拆环境背景——浅层传播看潜在蒸散发（q=0.496），深层传播看降水（q=0.391），且双因子组合的解释力全部高于单因子。作者同时交底：土壤湿度来自再分析资料而非实测，传播时间是统计特征时间尺度而非水分真实传递时间，解释力不等于因果。
缺的是哪一块
干旱传播研究已经走完三步：多尺度指数匹配与相关分析回答「何时响应」；Copula 函数、贝叶斯网络等概率方法回答「发生可能性多大、临界条件在哪」；再往后是传播速率、动态传播时间和非平稳关系的刻画。但一个结构性缺口在数据侧：多数研究用的是表层土壤湿度或综合根区土壤湿度，很少在统一框架下比较不同土层的传播特征——气象水分亏缺由浅入深的响应差异因此刻画不出来。内蒙古已有分层研究证实事件特征、传播概率和触发阈值存在明显垂向差异，而云南这种干湿季分明、岩溶地貌发育、水热梯度复杂的区域可能进一步放大土层间与区域间的差异。
原文图 1：(a) 云南在中国的位置，(b) 高程（6674 m 至 −54 m，西北高东南低），(c) 土地利用——林地占主导，滇中盆地交错分布耕地与城乡建设用地。
数据与分区怎么搭
气象侧用 ChinaMet 0.01° 月降水、气温、潜在蒸散发（1981—2020），对气候水分平衡量 D=P−E 做 1–24 个月滑动累积、按历月非参数标准化（Gringorten 经验概率 + 标准正态反函数）得到 SPEI-1～24；土壤侧用 ERA5-Land 约 0.1° 逐日体积含水量，聚成月值后按历月标准化得三层 SSI，云南 19 个地面站的独立评价支持其 7–28 cm 层能较好表征 40 cm 根区土壤水分的时间变化。干旱判定：SPEI<−1 为气象干旱、SSI<−1 为土壤干旱；事件识别用三阈值游程理论（识别阈值 R₁=−0.5、单月事件保留阈值 R₂=−1.0、事件合并阈值 R₀=0），逐格点统计频次、历时、烈度、峰值亏缺。
分区以 129 个县级单元为对象，用高程、坡度、土壤类型、NDVI、多年平均年降水量、年均气温、土地利用七个变量做带空间邻接约束的 Ward 层次聚类，切成四个地理分区：
分区 · 定位 · 环境特征
Ⅰ区 · 岩溶林地区（滇东及滇东北丘陵高原） · 多年平均年降水量全省最低，866 mm
Ⅱ区 · 核心农业区（滇中及中部高原盆地） · 平均高程全省最低 1402 m，降水、气温、NDVI 均为四区最高
Ⅲ区 · 高寒农林区（滇西北高寒高海拔山区） · 平均高程 2868 m、坡度 25.4° 为最高，年均气温最低 10.9 ℃
Ⅳ区 · 南部森林区（南部及西南部季风森林区） · 平均坡度最低 15.3°，地形平缓、森林分布广
传播时间由逐格点 SPEI-k 与 SSIj 的 Pearson 相关系数最大值对应的累积尺度确定；传播概率用高斯、Clayton、Frank、Student-t 四种 Copula 按 AIC 优选（结果统一选了高斯型）计算「气象干旱条件下该土层出现土壤干旱」的条件概率；触发阈值定义为条件概率达到 0.60 时对应的最弱 SPEI 值，从 −0.50 以 0.01 步长扫到 −3.00，全程达不到 0.60 的格点记为无有效触发阈值——0.60 这条线是操作性标准，定得再高，大量深层格点就拿不到有效阈值。
结果数字
长期趋势：气象干旱没显著变干，土壤在显著变干
Theil-Sen 斜率 + Mann-Kendall 检验（p<0.05）的全省结果：
指数 · 年际 Sen 斜率 · 月际 Sen 斜率 · 显著下降格点占比
SPEI-1 · −0.0056 a⁻¹（不显著，p=0.08） · −0.0061 a⁻¹（不显著） · 20%
SSI1 · −0.0183 a⁻¹ · −0.0188 a⁻¹ · 92%
SSI2 · −0.0192 a⁻¹ · −0.0195 a⁻¹ · 87%
SSI3 · −0.0206 a⁻¹ · −0.0198 a⁻¹ · 75%
同一个省、同一段时间，降水蒸散平衡的长期趋势过不了显著性检验，土壤水分却大面积显著变干，而且越深斜率越陡。分区上，Ⅰ、Ⅱ区三层 SSI 都在大范围显著下降；Ⅲ区和Ⅳ区的显著下降范围随土层加深收缩——Ⅲ区从 SSI1 的 57% 缩到 SSI3 的 17%，Ⅳ区从 98% 缩到 73%。2009—2013 年西南持续干旱期在序列里清晰可见：浅层土壤水分 2009—2010 年最低、随后略有回升、与 SPEI-1 同步，深层则持续偏低。
原文图 2：1981—2020 年云南省区域平均的 SPEI-1 与三层 SSI 年际变化，红色虚线为 Theil-Sen 趋势。SPEI-1 斜率 −0.0056 a⁻¹ 置信区间跨零（p=0.08），三层 SSI 斜率 −0.0183 至 −0.0206 a⁻¹、均 p<0.001。
原文图 3：月尺度 SPEI-1 与三层 SSI 序列（−2.5 至 2.5 的锯齿）与趋势线。2009—2013 年持续干旱期内多数月份为负，SSI3 大部分月值为负——深层土壤干旱的持续时间更长。
传播三要素：越深越慢、越难、越苛刻
层次 · 平均传播时间 · 平均传播概率 · 平均触发阈值
SSI1（0–7 cm） · 1.41 个月 · 51% · −1.32
SSI2（7–28 cm） · 2.08 个月 · 48% · −1.46
SSI3（28–100 cm） · 4.64 个月 · 40% · −1.70
空间上三条指标各有走向：SSI1、SSI2 的较长传播时间集中在滇西北高海拔区，SSI3 的较长传播时间更多出现在滇中及其周边；SSI3 的有效触发阈值区域明显收缩，Ⅳ区及Ⅲ区北部大量格点无有效阈值，深层触发阈值的空间连续性弱。干旱事件特征也随深度改形：随土层加深，事件频次减少、平均历时和烈度增加、峰值亏缺幅度降低——单月亏缺减弱而持续性和累积性增强，深层土壤水分有明显的「惰性」，干旱不是被单月亏缺高频打出，而是长期亏缺累积出来的。
原文图 6：SPEI-1 与三层 SSI 干旱事件的频次、平均历时、平均烈度、峰值亏缺空间分布（前四列）加分区箱线图。SPEI-1 各项分布相对均一，SSI 事件的空间差异随土层加深增大，SSI3 最为显著。
原文图 7：SPEI 向三层 SSI 的干旱传播时间空间分布与分区统计。SSI1、SSI2 全省以浅色（短历时）为主，SSI3 的深色长历时区集中出现在滇中及滇西北；分区箱线图上Ⅱ区 SSI3 传播时间中位数最高。
原文图 8：干旱传播概率空间分布与分区统计。SSI1、SSI2 的高概率区（红）集中在滇中，SSI3 全省整体转淡；箱线图显示各区 SSI3 传播概率中位数全面低于浅层。
原文图 9：干旱传播触发阈值空间分布（灰色 N/A 为条件概率达不到 0.60 的格点）与有效阈值覆盖率。SSI3 图上Ⅳ区大片灰色，分区覆盖率 SSI3 依次为 75%（Ⅰ区）、87%（Ⅱ区）、60%（Ⅲ区）、39%（Ⅳ区），而 SSI1 普遍在 78%–100%。
分区对比：两种深层干旱响应
Ⅱ区（核心农业区）与Ⅳ区（南部森林区）的对比最鲜明：Ⅱ区深层土壤干旱呈低频次、长历时、高烈度特征，SSI3 传播时间较长但有效触发阈值覆盖率较高，是典型的长期累积型响应；Ⅳ区深层干旱频次较高而持续性、累积性较弱，SSI3 传播概率和有效阈值覆盖率都偏低。Ⅲ区的南北差异最突出——北部高海拔与南部林地之间传播概率和有效阈值分布明显不同。Ⅰ区（岩溶区）呈局地高低值交错、空间不连续的形态，作者给出的潜在解释是岩溶裂隙、浅薄土层和快速入渗改变了降水补给路径与储水过程，但聚类变量里没有岩性或岩溶发育程度，也未直接观测裂隙流和地下水补排。
谁在解释传播差异：地理探测器
因子探测器的 q 值（解释力）按土层分化得很清楚：
传播指标 · 浅层（SSI1）最强因子 · 深层（SSI3）最强因子
传播概率 · 潜在蒸散发 0.496 · 降水 0.391
触发阈值 · 潜在蒸散发 0.323 · 降水 0.100（跨尺度敏感性存疑）
传播时间 · 气温 0.206 · 降水 0.253
风险探测器给出方向：SSI1 触发阈值随潜在蒸散发等级升高总体趋近于 0——高蒸散背景下较弱的气象干旱就能触发浅层土壤干旱；SSI2、SSI3 传播时间随降水等级升高总体缩短。交互探测的结果是：降水+潜在蒸散发组合对三层传播概率的解释力均最高，降水+气温组合对三层传播时间均最高；触发阈值上 SSI1 以气温+潜在蒸散发组合最高、SSI2 和 SSI3 以降水+潜在蒸散发组合最高。所有传播指标的最强双因子解释力都高于相应单因子，表现为双因子增强或非线性增强。
原文图 10：八个环境因子对九项传播指标的 q 值热力图。潜在蒸散发列在传播概率/触发阈值的浅层行里最深红（0.496、0.323），降水列在 SSI3 传播概率行（0.391）和传播时间行（0.253）最深；NDVI、坡度、土壤类型、土地利用的 q 值普遍在 0.2 以下。
作者承认的限制
- 土壤湿度是再分析不是实测。ERA5-Land 由 ERA5 大气强迫场驱动 HTESSEL 陆面模式模拟生成，受大气强迫、模式参数化和网格尺度限制；对云南复杂地形与岩溶背景，裂隙入渗、浅薄土层与地下水补排等局地过程难以充分刻画——SSI 传播关系主要反映区域尺度特征，绝对体积含水量不作主要结论依据（全国尺度评价也显示其绝对含水量可能存在系统偏差）。
- 传播时间是统计量不是物理量。传播时间由 Pearson 相关最大值对应的 SPEI 累积尺度确定，反映线性统计联系最强的特征累积尺度，不等同于水分在土壤剖面中的真实物理传递时间；对非线性响应的刻画有限。
- 解释力不是因果。地理探测器只反映候选因子与传播指标空间差异的一致程度，不能作确定性归因。
- 因子体系有缺口。缺少灌溉、实际蒸散发、根系吸水以及高海拔地区积雪冻融等过程数据；NDVI 只取 2001—2020 年平均值表征空间背景，不用于刻画植被的年际变化与动态反馈。
- 0.60 是操作性标准。阈值判定所用的条件概率标准并不统一，取 0.60 是为了兼顾传播关系识别与空间有效覆盖——若采用更高标准，大量深层格点将无法获得有效阈值；传播特征数值不视为普适参数。
一句话判断
气象干旱传导到土壤里是三个门槛同时抬高：云南 40 年数据给出的垂向规律是 1.41→2.08→4.64 个月的等待、51%→48%→40% 的概率、−1.32→−1.46→−1.70 的触发阈值——浅层的命门在蒸散发，深层的命门在降水，分区预警阈值若不分土层不分区域，等于把深层干旱当浅层干旱管。