观潮 · 模式预报与气象 AI 观察
石羊河流域干热复合事件 26 年增多：每多一次气象干旱 NDVI 减少 0.0064，民勤增速最显著
> 高温与干旱撞在同一段时间里，对植被的打击大于两个事件单独影响之和。这篇用石羊河流域 5 个国家气象站 2000—2025 年 5—9 月气温降水资料加 5 个县区 NDVI 数据，按「90 分位高温阈值连续 3 天 + Ped 干旱指数判定干旱月」识别复合高温干旱事件，再把 NDVI 拆成趋势项与气象波动项做归因。结果：复合事件年平均 1.96 次且波动增多，民勤线性倾向率 0.738 次/10 a（p<0.01）；复合事件每多一次，全流域气象 NDVI 减少 0.0064；2023 年频次达 4 次峰值，对应气象 NDVI 跌到 26 年最低。作者自己点明两条硬限制：全流域只有 5 个国家站，回归只放了气温降水两个因子。
缺的是哪一块
复合高温干旱事件的环流特征、成因机制和演变趋势已有不少研究，北半球尺度上干热复合事件对高纬度植被多呈促进作用、对中低纬度多呈抑制作用，且滞后时间较短——2010 年俄罗斯高温干旱导致作物歉收、2020 年美国西部热浪干旱使植被生产力严重下降都有案可查。石羊河流域的问题在于：地处祁连山北麓、青藏高原东北侧，地形阻挡水汽输送，常年降水稀少、干旱和高温事件频繁，是气候变化敏感区和生态脆弱带，但针对该流域复合高温干旱事件及其对植被影响的研究此前缺乏。2023 年流域罕见强高温热浪叠加春旱、春末夏初旱及持续伏旱，经济损失 83 815 万元。
这篇补的是事件识别 → 时空特征 → 植被响应的完整链条，关键方法动作是把 NDVI 里人类生态工程的贡献（趋势项）剥离出去，只拿气象驱动的波动项（气象 NDVI）去对复合事件频次做相关与回归。
事件怎么定义、数据怎么选
复合高温干旱事件的判定分两步：
- 高温事件：把单站逐年生长季日最高气温升序排列，取第 90 百分位为该年阈值，多年平均得该站高温总阈值；日最高气温连续 3 天超过总阈值记为一次高温事件。
- 干旱月：Ped 气象干旱指数 X>0 定义为干旱月（该指数此前已在石羊河流域验证，识别出的干旱事件与真实干旱过程基本吻合）。一个干旱月内出现一次及以上高温事件，记为一次复合高温干旱事件；年内生长季所有复合事件之和为当年频次。
数据：流域内仅有天祝（乌鞘岭站）、古浪、凉州、民勤、永昌 5 个国家气象站——2004 年后陆续建立的区域站资料年限差异大，未纳入分析。NDVI 为 500 m 分辨率，经云层、太阳高度角、气溶胶等校正，全流域值由 5 个县区平均。
气象 NDVI 的分离：NDVI = 趋势项（5 年滑动平均提取，代表生态保护与修复措施推动的稳定增长）+ 气象项（气象要素驱动的波动）+ 随机噪声（微小可忽略）。趋势检验按蒙特卡罗模拟，趋势系数绝对值大于 0.306、0.365、0.443 分别对应 α=0.1、0.05、0.01 的显著性水平。
原文图 1：石羊河流域地理及气象站点分布。流域地处青藏高原、黄土高原与蒙新高原交汇过渡带，海拔自西南祁连山边坡的 2500–4850 m 降至东北部 1300–2000 m，五个国家气象站（天祝、古浪、凉州、民勤、永昌）标注于图上。
结果数字
复合事件：空间差异小、时间增多明显
年平均频次全流域 1.96 次，空间格局古浪（2.08）> 天祝（2.00）> 凉州（1.96）> 民勤（1.92）> 永昌（1.85）——空间差异小，与流域范围不大、常受同一天气系统控制有关。时间上不均匀：2023 年多达 4 次（各地均达 4 次），2009、2010、2016、2018、2020、2024、2025 年均为 3 次；2020 年后呈显著增多趋势。
站点/区域 · 年均频次（次） · 线性倾向率（次/10 a） · 显著性
民勤 · 1.92 · 0.738 · p<0.01
凉州 · 1.96 · 区间内 · p<0.05
全流域 · 1.96 · 区间内 · p<0.05
天祝 · 2.00 · 区间内 · p<0.1
古浪 · 2.08 · 区间内 · 不显著
永昌 · 1.85 · 0.280（最小） · 不显著
全流域及各地线性倾向率落在 0.280–0.738 次/10 a 之间。背景是流域气温以 0.421 ℃/10 a 显著增温、降水仅以 6.874 mm/10 a 微弱增加——增温快、补水慢，复合事件频发的结构没变。
原文图 2：2000—2025 年生长季石羊河流域复合高温干旱事件年平均频次的空间分布。古浪最高 2.08 次，永昌最低 1.85 次，全流域均值 1.96 次，柱间差异相对较小。
原文图 3：复合高温干旱事件频次的年际变化。2023 年全流域及各地均达 4 次，为 26 年峰值；2020 年后各站频次整体抬升，民勤增多趋势通过 p<0.01 检验。
NDVI 与气象 NDVI：一个显著增、一个纯波动
生长季 NDVI 全流域年均 0.2922，空间差异大：天祝 0.4881 > 凉州 0.3067 > 永昌 0.2779 > 古浪 0.2615 > 民勤 0.1268——南部山区森林草原广布对高值，北部民勤荒漠戈壁对低值。全流域及各地 NDVI 均呈波动增大趋势（线性倾向率 0.007–0.054/10 a，全部 p<0.01），作者把这归因于 2010 年启动的石羊河流域治理工程（持续调水、退耕还林、禁牧封山）。
剥离趋势项后的气象 NDVI 则完全不同：年际波动明显、无清晰规律，线性倾向率 −0.003~0.003/10 a 且全部不显著。全流域气象 NDVI 2012 年最大 0.0185，2002 年次大 0.0180，2023 年最小 −0.0344——恰是复合事件频次峰值年。
原文图 4：生长季年平均 NDVI 的空间分布。天祝 0.4881 最高，民勤 0.1268 最低，与海拔、地形决定下的植被类型差异对应。
原文图 5：生长季 NDVI 年际变化。柱状为全流域均值，折线为各站；2010 年后明显增大，2023 年出现偏低点，与当年复合事件峰值对应。
原文图 6：气象 NDVI（剔除趋势项后的气象驱动波动分量）年际变化。全流域及各站围绕零轴波动，2012 年最大 0.0185、2023 年最小 −0.0344，线性趋势均不显著。
频次对气象 NDVI 的定量影响
复合事件频次与气象 NDVI 全流域及各地均为负相关：全流域（R=−0.476）和古浪（R=−0.478）在 p<0.05 水平显著，天祝（R=−0.370）较显著，永昌（R=−0.084）、凉州（R=−0.287）、民勤（R=−0.304）不显著。
区域 · 每多 1 次复合事件的气象 NDVI 减量
全流域 · 0.0064
古浪 · 0.0116
天祝 · 0.0070
凉州 · 0.0040
永昌 · 0.0014
民勤 · 0.0013
原文图 7：复合高温干旱事件年频次与年平均气象 NDVI 的关系，六个面板分别为全流域、永昌、民勤、凉州、古浪、天祝，图中标注线性回归方程与相关系数 R。全流域拟合 y=−0.0064x+0.0071（R=−0.476），古浪斜率最大（−0.0116，R=−0.478）。
结论指向明确：复合高温干旱事件是制约该流域 NDVI 动态变化的关键气象因素，对作物生长发育有抑制作用；准确预报复合事件后，可用线性趋势方程对 NDVI 做线性范围内的预测。
作者承认的限制
- 站点太少、代表性有限。研究仅依赖流域内 5 个国家气象站，站点数量少且空间分布不均，「其结果只反映了流域整体平均情况，难以精准刻画复合高温干旱事件的局地差异与空间细节」。作者写明未来非常有必要引入更高密度的区域气象观测数据。
- 回归因子不全。线性回归模型只考虑高温与干旱（气温、降水）两个核心气象因子，未纳入蒸发、日照等其他影响 NDVI 的因子，可能导致回归系数估计偏差——作者认为这也可能是气象 NDVI 呈微弱增大趋势的原因。
- 人类活动未单独定量。未具体分析 2010 年以来持续调水、退牧还草、退耕还林等生态工程对 NDVI 的人为驱动作用，作者把这条与上一条并列为「本研究的主要局限」。
一句话判断
26 年 5 站资料把「干热同季更频繁」落成了具体数字——民勤每 10 年多 0.738 次、2023 年 4 次峰值对应气象 NDVI 26 年最低、每多一次复合事件全流域气象 NDVI 平均减 0.0064；但 5 个国家站撑不起局地细节、两因子回归拿不下全部波动，作者把这两条限制写在了结论旁边，用这篇方法套用时得先掂量自己手里的站网密度。