观潮 · 模式预报与气象 AI 观察

江苏气象干旱传成农业干旱平均要 3 个月:25 年 62 站记录给出四级触发概率 0.51、0.65、0.75、0.83

2026-10-14

> 气象干旱是降水亏缺,农业干旱是土壤水分亏缺,两者之间隔着一个传递过程——隔多久、多大降水亏缺才会把土壤拖进干旱,决定农业干旱预警能提前多久发布。这篇论文用 1990—2025 年江苏 62 个气象站资料加 ERA5-Land 土壤湿度,以 SPI 和 SSI 分别表征两类干旱,联游程理论识别干旱事件、最大皮尔逊相关系数定传播时间、Copula 函数算触发概率与阈值、地理探测器归因空间分异。核心数字:传播时间集中在 3—5 个月(占全年月份 83.3%),站点尺度以 3 个月为主(61.3%),南短北长;轻、中、重、极端气象干旱条件下农业干旱发生概率的全省平均分别为 0.51、0.65、0.75、0.83;降水是传播时间空间分异的主导因子(q=0.56),降水与土壤湿度的交互解释力最高(q=0.74)。作者同时交代了两条可能动摇结论适用范围的限定:结果对雨养旱作农田解释力更强,水稻田没算;灌溉的调节作用因资料获取困难没有纳入。

缺的是哪一块

干旱传播研究在长江流域(传播时间约 2 个月)、全国尺度(气象—农业、农业—地下水、气象—地下水平均传播时间分别为 4.1、15.1、24.1 个月)、海河流域都做过,蒙古高原的多类型干旱级联阈值、沁河流域气象—水文干旱阈值也有量化。江苏的位置特殊:地处南北气候过渡带、长江—淮河流域下游,水热条件和农业环境空间差异显著,是全国重要粮食产区,但既有研究多停留在单一类型干旱的时空演变或具体作物的水分响应,气象干旱到农业干旱的传播时间、触发概率和传播阈值在省内缺乏系统性量化。2022 年全省春夏连旱叠加长江流域夏秋冬连旱,是这个问题现实性的直接注脚。这篇要补的是一套多维分析体系:多类型干旱指数 + 传播时间 + 传播阈值 + 空间分异归因。

四个工具各干一件事

数据是 62 站逐日气温、降水、相对湿度、实际水汽压及海拔(1990—2025 年),配同期 ERA5-Land 逐日土壤含水量(0—7、7—28、28—100 cm 三层,0.1° 分辨率,双线性插值到站点)。SSI 用 0—100 cm 深度土壤湿度,三层按 0.07、0.21、0.72 加权合成——权重对应作物根系密度随深度的递减分布,水稻根系集中 0—20 cm、冬小麦返青前 0—40 cm、玉米可扎到 1 m 以下,0—100 cm 已覆盖大部分根系吸水。SPI 用伽马分布拟合,SSI 用广义极值(GEV)分布拟合。

方法链是四段:

方法 · 干什么 · 关键设定

游程理论 · 从 SPI/SSI 序列里识别干旱事件,提取频次、历时、强度 · 三级截断 R0=0、R1=−0.5、R2=−1;1 个月事件且指数大于 R2 的剔除;间隔 1 个月且间隔月指数小于 R0 的合并

最大皮尔逊相关系数 · 定传播时间:SSI-1 与 SPI1—SPI12 逐站算相关,最大相关对应的 SPI 尺度即传播时间 · —

Copula 函数 · 建 SPI 与 SSI 联合分布,算条件概率与触发阈值 · 3 种边缘分布(正态、极值、GEV)过 KS 检验优选;5 种 Copula(Gaussian、t、Clayton、Frank、Gumbel)按平方欧几里得距离选优;阈值以 −0.1 步长从 −0.50 逆向迭代,概率达 0.50 时的 SPI 值为阈值

地理探测器 · 归因传播时间空间分异:因子 q 值 + 交互探测 · 因变量为站点传播时间,解释变量为降水、土壤湿度、平均气温、饱和水汽压差(VPD)、海拔

干旱等级统一分四档:轻度(−1.0~−0.5)、中度(−1.5~−1.0)、重度(−2.0~−1.5)、极端(≤−2.0)。

先看底账:气象干旱在减弱,农业干旱在增强

1990—2025 年,江苏气象干旱事件数每站 28—43 次,平均历时 2.70—4.46 个月,平均强度 3.30—4.99;事件数、历时、强度都呈波动减弱,速率分别 −0.06/a、−0.019/a、−0.023/a。农业干旱事件数 16—33 次,明显更少,但平均历时 3.41—8.24 个月、平均强度 3.98—9.19,两项都高于气象干旱,且呈波动上升(+0.005/a、+0.03/a、+0.03/a)。2004 年是农业干旱最突出的一年,平均历时 10.13 个月、强度 12.38。

图1 江苏省气象干旱时空分布与年际变化
图1 江苏省气象干旱时空分布与年际变化

图1 气象干旱平均历时(a,d)、事件频数(b,e)、平均强度(c,f)的空间分布与年际变化:北部徐州、宿迁、连云港一带事件数最高(43 次),南京、镇江一带历时与强度更高,年际趋势线向下。

图2 江苏省农业干旱时空分布与年际变化
图2 江苏省农业干旱时空分布与年际变化

图2 农业干旱对应三要素的空间分布与年际变化:苏南苏州、无锡、常州事件数较高,扬州、徐州、盐城一带历时与强度更高,三条趋势线均向上——气象干旱减弱的同时农业干旱在增强。

「气象干旱减弱、农业干旱增强」这组反向趋势是全文最有业务含义的发现。作者的解释链:土壤水分受降水之外的多重调节——深层补给缓慢、植被蒸腾耗水放大脆弱性,降水要先满足植被吸收,剩余的才能补深层土壤,土壤水分恢复因此滞后;同时江苏近几十年明显增暖,作物蒸散量对气温和日照敏感,增温背景下作物需水量可能增加。一降一升,农业干旱风险的结构在变。

传播时间:3—5 个月占 83.3%,南短北长

逐月看,SSI-1 与多时间尺度 SPI 的相关系数最大值在 0.657—0.854,高值集中在 SPI2—SPI7;按最大相关对应的 SPI 尺度,各月传播时间 2—7 个月,其中 3—5 个月的月份占全年 83.3%。季节上春季和夏季平均传播时间均为 4.67 个月,冬季 3.67 个月,秋季最短 3.0 个月。

图3 SSI-1与多时间尺度SPI的逐月相关系数
图3 SSI-1与多时间尺度SPI的逐月相关系数

图3 各月份 SSI-1 与 SPI1—SPI12 的相关系数热力图:星号标出该月最大相关所在的 SPI 时间尺度(即传播时间),高相关区集中在 SPI2—SPI7。

空间上,各站最大相关系数 0.475—0.684、平均 0.612;传播时间以 3 个月为主(占站点 61.3%),不超过 3 个月的站点占 72.6%,总体南短北长——苏南水分补给足、水循环快,土壤水分对近期降水异常响应快,苏北降水相对少,土壤水分受长期累积和其他陆面因素共同调节。

图4 最大相关系数与传播时间的空间分布
图4 最大相关系数与传播时间的空间分布

图4 (a) 1990—2025 年各站 SPI 与 SSI 最大相关系数(0.475—0.684);(b) 传播时间空间分布与站点占比饼图——3 个月占 61.3%,4 个月占 12.9%,5 个月占 9.7%。

触发概率与阈值:等级每升一档,概率抬一格、阈值深一截

Copula 条件概率的结果是一条单调上升的阶梯:

气象干旱等级 · 农业干旱发生概率(全省平均) · 站点变化范围

轻度 · 0.51 · 0.43—0.62

中度 · 0.65 · 0.55—0.77

重度 · 0.75 · 0.63—0.86

极端 · 0.83 · 0.70—0.95

空间上概率呈中部较高、南北两侧较低:淮安、扬州及南京周边整体偏高,徐州、连云港及苏州部分地区偏低——苏南低可能与水网密布有关,苏北低可能与作物类型和土壤蓄水能力有关。

图5 四级气象干旱触发农业干旱的概率分布
图5 四级气象干旱触发农业干旱的概率分布

图5 轻度、中度、重度、极端气象干旱条件下农业干旱发生概率的空间分布:从 (a) 到 (d) 颜色整体由绿转红,淮安、扬州、南京周边概率最高。

阈值一端同样单调。触发轻度农业干旱需要 SPI 达 −1.00~−0.50,中度 −2.80~−1.40,重度 −3.60~−2.00,极端 −3.90~−2.60——农业干旱等级越高,需要的气象水分亏缺越深。能凑够 0.50 触发概率的站点也随等级锐减:轻度 60 站、中度 58 站、重度 34 站、极端 31 站。阈值高的地区(数值更接近 0)意味着较弱的气象干旱就能触发,苏中及沿江地区因此对传播相对敏感;苏北西北部、东部沿海及苏南部分地区阈值低,传播抵抗力强。

图6 四级农业干旱的传播阈值空间分布
图6 四级农业干旱的传播阈值空间分布

图6 轻度、中度、重度和极端农业干旱的 SPI 触发阈值:灰色点为未拟合或未达 0.50 概率的站点,等级越高灰色点越多。

归因:降水主导,降水×土壤湿度交互最强

地理探测器的因子解释力排序:降水 q=0.56 最高,土壤湿度 0.53、平均气温 0.39 次之,海拔 0.18 和 VPD 0.08 较低。交互探测更关键——所有双因子组合的解释力都高于相应单因子,其中降水与土壤湿度交互 q=0.74 最高,降水×温度 0.69、土壤湿度×温度 0.67。传播时间的空间差异不由单一因子决定,而是水分输入、土壤储水和蒸散需求协同调节的结果;气温的作用还受水分可利用程度制约——土壤水分严重受限时,实际蒸散未必随气温同步上升。

图7 驱动因子及其交互作用对传播时间的解释力
图7 驱动因子及其交互作用对传播时间的解释力

图7 (a) 五个驱动因子的 q 值:降水 0.56、土壤湿度 0.53、气温 0.39、海拔 0.18、VPD 0.08;(b) 交互探测矩阵:降水×土壤湿度 0.74 为全场最高。

作者承认的限制

论文在讨论部分写了整整两大段边界,逐条如下:

限制 · 原文口径

作物过程没算 · 「尚未结合具体作物的生育期、需水量和产量资料」,水稻田的农业干旱还受田面水层保持能力、灌溉水源条件和田间水分管理影响

适用范围 · 结果「对雨养或旱作农田及以土壤水分亏缺为主要限制因子的农业干旱过程解释力更强,对于灌溉农田,尤其是水稻」,需后续结合水层保持、灌溉、作物类型、生育期和产量再研究

插值偏差 · 站点降水与 ERA5-Land 土壤湿度的双线性插值配准「仍存在一定偏差」,未来需实测墒情、多源遥感土壤湿度与不同再分析数据交叉验证

灌溉缺失 · 「由于区域尺度灌溉资料获取困难且空间分布存在较大差异,本文未考虑其对干旱传播过程的调节作用,这可能带来一定的不确定性」

土地利用变化 · 退耕还林、城市扩张等通过改变地表反照率和蒸散发影响传播,「作用方向和强度具有明显的区域差异」,未纳入

一句话判断

对江苏这样的水网区,气象干旱预警到农业干旱预警之间隔着大约 3 个月,且极端气象干旱有 0.83 的概率把农业干旱拖下水——但这套阈值出自雨养口径的土壤水分平衡,搬到灌溉稻区之前要先打折。

在观潮中打开(含配图与交互)→