干旱西北中国卫星观测植被绿度变化的季内不对称性及其气候关联

《Frontiers in Plant Science》：Intra-seasonal asymmetry in satellite-observed vegetation greenness changes and their climatic associations in arid Northwest China

　　

干旱生态系统对气候变化高度敏感，然而区域绿化是否伴随整个生长季植被绿度的一致变化仍不清楚。研究人员利用1999—2018年SPOT/VEGETATION归一化差分植被指数（NDVI）和气象数据，采用月际NDVI差值指标（VNDVI）刻画干旱西北中国卫星观测植被绿度的季内变化，并检验其与同期及前期气候条件的关系。研究将5—10月定义为区域比较的通用分析窗口，其中5—7月代表早至峰值生长季，8—10月代表晚生长季。尽管多数月份NDVI呈正长期趋势，VNDVI趋势却表现出显著季内不对称性：6月和7月VNDVI显著增加，而10月显著减少，表明区域绿化伴随的是月际绿度变化的重新分配，而非生长季内的均匀增强。这些趋势亦具明显空间异质性：5月负VNDVI趋势集中于阿尔泰山和西天山，6—7月正趋势在中天山及塔里木盆地北缘最显著。VNDVI与气候变量的关系随月份和区域变化：降水通常与月际绿度变化正相关，而温度和辐射在若干月份和区域呈负相关。层次分解（hierarchical partitioning）进一步揭示，气候变量对NDVI变异的独立贡献比例因植被类型而异，其中农田在早生长季的温度贡献占比尤其大。

研究背景与问题提出：干旱与半干旱生态系统因水分限制而对气候变率高度敏感，卫星遥感归一化差分植被指数（NDVI）广泛用于记录区域绿化，但既有研究多关注年或生长季均值、积分或物候转换日期，可能掩盖生长季内部的实质性变化。区域绿化未必意味着各月绿度均匀增加，而可能表现为特定时段增幅更强、峰值提前或月际差值改变。在干旱西北中国（ARNC），强气候梯度、复杂地形与植被覆盖空间异质性并存，既往绿化研究多聚焦年或季尺度，季内月际绿度变化的不对称性、空间分异、植被类型差异及其与同期/前期气候的统计关联尚未明晰。该研究基于此，结合1999—2018年SPOT/VEGETATION NDVI与网格气象数据，引入月际NDVI差值（V_NDVI=NDVI(t)?NDVI(t?1)）作为描述性指标，量化NDVI与V_NDVI的时空趋势，并评估其与气温、降水、地表太阳辐射的同期及前期（滞后0—3月）偏相关关系，采用层次分解比较不同植被类型下气候变量对年际NDVI变异的独立解释贡献。论文发表于《Frontiers in Plant Science》。

主要技术方法概括：研究人员使用1 km月度SPOT/VEGETATION NDVI（最大合成法降云），以1999—2018年5—10月为通用窗口；气候数据取自CHELSA V2.1约1 km网格的月气温、降水、辐射；土地覆盖采用ESA CCI重分类稳定像元（农田、草地、灌木、森林、稀疏植被）。趋势分析用Theil–Sen斜率+Mann–Kendall检验；气候关联先按像元-目标月-变量筛选滞后0—3月中绝对Pearson相关系数最大者，再做控制其余两变量的偏相关；相对重要性用R语言hier.part包层次分解，得温度、降水、辐射独立贡献占比及模型R2。所有分析在1 km植被掩膜与稳定土地覆盖掩膜下像元级进行。

研究结果：

3.1 干旱西北中国V_NDVI的时空格局与趋势

3.1.1 季节格局与长期变化：区域多月NDVI均正趋势（5月0.013、6月0.029、7月0.042、8月0.029、9月0.024 decade^-1，10月正但不显著），但V_NDVI趋势不对称——6月（+0.016）、7月（+0.013 decade^-1）显著正，10月（?0.019 decade^-1）显著负，5月弱负不显著。说明绿化主要来自早—峰值段月际增量放大，晚季9—10月降幅扩大。

3.1.2 空间格局：5月负V_NDVI趋势在阿尔泰、西天山；6—7月正趋势集中中天山、塔里木北缘；8月弱异质，9—10月天山/阿尔泰负趋势扩展，10月最清楚。

3.2 V_NDVI与同期/前期气候关联

3.2.1 年际偏相关：早季降水对V_NDVI多正关联，辐射多负；5月V_NDVI与前期温、降水正，6月与前期温负、前期降水正，7月与前期降水正、前期辐射负。晚季8月V_NDVI与前期辐射负；10月与前期降水负、前期辐射负。NDVI与V_NDVI气候关联在晚季常不一致。

3.2.2 空间异质：早季中天山/塔北降水-V_NDVI正关联重叠正趋势区；阿尔泰5月V_NDVI与前期温正、前期降水负共存；晚季空间符号更乱。

3.3 植被类型间气候对NDVI变异的独立贡献：区域模型R2=0.78，5月与10月温度占独立贡献约60%。草地早季温度占比>60%（R2=0.83）；灌木季节波动大（R2=0.64）；农田5月温度占约72%但R2仅0.57（灌溉管理等未解释）；森林温、降水贡献交替（R2=0.69）；稀疏植被贡献分散且R2最低（0.54）。草地、森林气候耦合强，农田、稀疏植被受非气候因素干扰大。

讨论浓缩：季节均值绿化掩盖了月际轨迹重组——ARNC绿化实质是早—峰值段（尤6—7月）月际增量放大、晚秋（10月）月际降幅加深的不对称再分配，非均匀生理加速。V_NDVI趋势由相邻两月NDVI趋势差决定，故6月正趋势反映6月比5月增得多、10月负趋势反映10月比9月增得少。气候关联上，早期降水正向、辐射负向符合水分限制逻辑但仅为统计偏相关；温度在山地可能兼表雪融/土温正向与高蒸发需求负向，不能径判生理机制。空间上，阿尔泰早季温正雨负或涉雪相/云辐射协变；中天山/塔北正V_NDVI趋势与前期降水正关联空间重叠，但不可等同因果。植被类型间，草地温控突出、农田早期温控高占比但低R2提示管理因子重要、稀疏植被低R2反映背景反射与局地水情干扰。

结论翻译：1999—2018年干旱西北中国区域绿化并未在生长季内均匀表达，而是表现为月际植被绿度变化的显著再分配。V_NDVI在6月和7月显著增加、10月下降，说明长期NDVI增益更多集中于早—峰值生长季而非晚秋。这些时间偏移在空间上异质，阿尔泰、天山与塔里木北缘格局对比明显。气候关联随月份、区域、植被类型变化：降水多与早季绿度变化正相关，温度与辐射呈季节拮抗关系。层次分解显示气候变量相对重要性与联合解释力因植被类型而异，草地和森林的气候—NDVI耦合强于农田与稀疏植被。综上，季节均值绿化可掩盖季内绿度轨迹的实质重组，解析月际变化为水分限制生态系统植被—气候关系提供了更高时间分辨率的框架。