基于多源遥感数据解析中国多种干旱类型对总初级生产力的影响

《Agronomy》：Unraveling the Responses of Gross Primary Productivity to Multiple Drought Types Across China Using Multi-Source Remote Sensing Data

　　

干旱是气候变化下影响陆地生态系统功能的主要因素之一。然而，不同干旱指数间的相互作用及其对总初级生产力（GPP）的响应机制仍不清楚。研究人员利用2003—2020年中国多源遥感数据，采用标准化降水蒸散指数（SPEI）、温度条件指数（TCI）、土壤湿度条件指数（SMCI）和植被条件指数（VCI）表征干旱时空动态，并量化干旱指数与GPP的响应机制。偏相关分析与结构方程模型（SEM）用于区分不同干旱类型对GPP的直接和间接效应。结果表明，2003—2020年GPP呈波动上升趋势，介于163.71至458.56 g C·m?2之间。2005年较2004年，GPP在华北下降14.47%、华中13.93%、华东12.05%、华南8.45%、西北4.21%、东北4.15%，而西南微增0.34%。不同类型干旱特征各异：SMCI干旱频率44.35%、平均持续75.64天；SPEI频率最低（12.54%）、强度最高（0.14）；TCI频率最高（60.22%）、持续最长（98.46天）；VCI严重度最低（1.96）。TCI被认为是GPP变化的主要驱动因子，各区均呈显著正相关（r=0.45–0.88，p<0.001），而SPEI影响较弱。“TCI→VCI→GPP”路径为全国尺度主导路径。该研究为评估未来气候变化下GPP响应特征提供方法。

研究背景方面，干旱作为气候变化下陆地生态系统失调的主要驱动因子，具有高频、广覆盖、长历时和形成机制复杂等特点，严重威胁全球水安全、区域社会经济与生态稳定。传统地面观测衍生的干旱指数（如SPI、SPEI）受站点稀疏与空间插值不确定性的限制，难以直接捕捉大尺度植被生理与土壤水分状态；而单一指数往往孤立提取某类水文气候信号，导致对GPP影响的强度、方向与主导驱动认识不一致。现有研究多依赖两两相关或单因子诊断，无法在共线性强的中介变量网络中剥离独立与交互效应。为此，研究人员在《Agronomy》发表该研究，集成2003—2020年中国多源遥感与气象再分析数据，构建气象—农业—土壤水分—热异常四维干旱框架，旨在厘清多维干旱对GPP的联合调控路径。

关键技术方法上，研究人员选取MODIS派生NDVI、地表温度（LST）、土壤湿度（SM）及降水（P）、潜在蒸散（PET）数据，统一重投影至5 km网格并将日值聚合为8天尺度以匹配GPP产品；计算SPEI（基于P与PET水量平衡）、VCI（NDVI极值标准化）、TCI（LST极值标准化）、SMCI（SM极值标准化）；采用多阈值游程理论在各网格独立识别干旱事件并统计频率、历时、强度与严重度；通过偏相关分析（控制其余三指数）提取GPP与各干旱指数的条件关联；构建结构方程模型（SEM，最大似然估计）将四指数作为外生变量、GPP为内生变量，量化直接/间接效应与“TCI→VCI→GPP”等路径系数，以χ2/df<3、RMSEA≤0.08、CFI≥0.9为拟合准则。

研究结果部分，3.1节时序动态显示，全国GPP在163.71–458.56 g C·m^?2间波动上升，2005年多数区显著回落（华北?14.47%等，西南+0.34%），SPEI无一致干湿趋势（西北均值0.11、负值频率6%远低于东北56%），VCI与TCI多数区2013年附近偏高，SMCI在东北/华北上升、华南下降。3.2节空间异质性表明，SPEI干旱低频（12.54%）短历时（32.47天）但高强度（0.14）高严重度（8.0），集中于西北与高纬；VCI频率44.87%、历时30.44天、严重度1.96，随湿润度降低而加重；TCI频率最高（60.22%）、历时最长（98.46天），东南向西北增强；SMCI频率44.35%、历时75.64天，呈西北—东南梯度。3.3节偏相关显示，TCI与GPP在各区均最强正相关（东北r=0.88、华东0.83、西北0.80，p<0.001），SPEI偏弱不显著，VCI在华中（0.40）/华南（0.42）显著，SMCI在西南与TCI/VCI协同显著；TCI与SPEI、SMCI多呈负向，反映水热权衡。3.4节SEM表明，TCI对GPP直接路径系数0.42–0.98（东北/华北0.98最高、华中0.42最低），VCI为关键中介（TCI→VCI系数0.38–0.97，西南0.97），SMCI与SPEI主要通过TCI/VCI间接作用；全国主导路径为“TCI→VCI→GPP”，模型R2介于0.43（华中）至0.86（西南），除西南拟合稍弱（RMSEA=0.40）外其余达标。

讨论部分指出，GPP年际变率由温度主导而非水分直接约束，TCI在干旱半干旱区通过光合效率、物候与呼吸平衡施加一阶控制；水分因子经VCI或TCI间接调制，湿润区VCI为气候信号集成中介，西南多约束耦合但仍以温度路径为主；单指数（尤SPEI）不足以表征GPP动态，需多维复合框架。结论部分翻译如下：基于2003—2020年多源遥感数据，研究人员系统考察中国GPP对多种干旱类型的集成响应。GPP对干旱扰动高度敏感，响应非仅由水分限制支配，而是以温度主导、植被中介、水分调制的耦合机制。各干旱类型发生特征差异重塑区域水热约束，引致GPP对气候异常恢复力空间分异。（1）GPP波动上升（163.71–458.56 g C·m^?2），2005年典型旱年多数区同步下降（华北?14.47%等，西南+0.34%）；（2）SMCI低频长历（44.35%，75.64天），SPEI最低频（12.54%）但最高强度（0.14）严重度（8.0），TCI最高频（60.22%）最长历（98.46天），VCI最低严重度（1.96）；（3）TCI与GPP全域显著正相关（r=0.45–0.88，p<0.001），SPEI较弱；（4）“TCI→VCI→GPP”为全国主导路径，R2 0.43–0.86，证实温度主驱、植被中介、水分次调的响应框架。该研究为未来气候变化下GPP响应评估提供方法学支撑。