气候变化下青藏高原植被绿度与恢复力的对比格局及其环境关联

《Science of Remote Sensing》：Contrasting patterns of vegetation greenness and resilience across the Tibetan Plateau under climate change

　　

植被绿度增加，但推断恢复力（inferred resilience）未表现出一致性改善。绿化现象与推断恢复力的上升和下降同时发生。植被绿度与时间持久性（temporal persistence）具有不同的环境关联。活动层厚度（active layer thickness, ALT）与森林绿度相关联，但与时间持久性无关。

气候变化下青藏高原植被绿度与恢复力的对比格局——论文解读

青藏高原是全球海拔最高且对气候变化极为敏感的区域之一，近几十年的升温速率显著高于全球平均，由此引发的冰川退缩、多年冻土退化和活动层增厚等冰冻圈变化，深刻改变了区域土壤热状况和水文过程，并触发了高寒生态系统在植被生产力和分布上的显著响应。然而，以往针对青藏高原的研究多聚焦于植被绿度或生产力，对生态系统稳定性和恢复能力的关注相对不足。卫星反演的绿度指数（如NDVI）主要反映植被活动的平均水平，并不能直接表达生态系统在扰动后的恢复动态。在气候持续变化和极端事件增多的背景下，植被绿度上升可能伴随着对干旱等胁迫的脆弱性增加。因此，亟需在区域尺度上同时刻画植被绿度与恢复力的时空变化，并厘清二者的环境驱动因子是否一致。

为此，研究人员以2005—2020年为研究时段，以青藏高原自然植被为对象，对比分析了植被绿度与推断恢复力的时空格局及其与环境因子的关联。植被绿度采用核归一化差异植被指数（kNDVI）表征，推断恢复力则基于近红外植被反射率（NIRv）的一阶自回归系数[AR(1)]。研究结果表明，区域平均绿度显著上升，而推断恢复力并未表现出同步的、一致性的改善；绿化既可以伴随恢复力上升，也可以伴随恢复力下降；绿度与恢复力的主导环境因子存在明显差异。这一发现说明，广泛绿化不能简单等同于生态系统恢复能力的增强，将绿度与恢复力相结合才能更完整地评估高原生态系统的响应。

在方法上，研究人员利用了2000—2020年的MODIS遥感产品、再分析土壤温度数据、高分辨率气候数据和活动层厚度（ALT）数据，并将分析限制在2000年和2020年土地覆盖类型一致的自然植被像元上，分草地和森林分别建模。恢复力指标AR(1)基于60个月滑动窗口估计，窗口敏感性通过36、60、84和120个月窗口进行检验。趋势显著性采用傅里叶相位随机化替代检验，以考虑滑动窗口带来的时间依赖性。环境归因采用极端梯度提升（XGBoost）结合Shapley加法解释（SHAP）模型，预测因子包括平均降水量（Pmean）、平均气温（Tmean）、潜在蒸散发（PET）、干旱指数（AI）、温度变率（Tvar）、降水变率（Pvar）和平均土壤温度（STmean）。模型使用五折交叉验证（按像元分组）评估预测表现，并用留一变量分析和SHAP排序稳定性检验验证归因结果的稳健性。样本来源上，全植被、草地和森林分别有299、295和293个采样像元，产生约5.6万—5.7万个月度窗口观测。

研究结果如下：

3.1 植被绿度与恢复力的时间动态对比

区域平均kNDVI在2005—2020年间显著增加：全植被肯德尔τ=0.778，草地τ=0.694，森林τ=0.782（均P<0.001）。相比之下，NIRv衍生的AR(1)在三个植被域中均无显著的全时期趋势。分时段看，2005—2010年森林AR(1)显著下降（τ=-0.781，P<0.001），指示恢复力上升；2011—2020年各植被域AR(1)趋势均不显著。这表明绿度上升与恢复力变化在时间上并不同步。

3.2 空间格局与联合趋势

2005—2020年间，43.3%的合格像元kNDVI显著上升，3.6%显著下降；而AR(1)显著下降和显著上升的像元分别仅占13.6%和6.4%，80.0%的像元无显著趋势。联合趋势分析显示，仅9.3%的像元出现绿度与AR(1)同时显著变化，其中5.7%为绿化伴随恢复力上升，2.8%为绿化伴随恢复力下降。分时期看，2005—2010年联合显著变化面积为25.0%，2011—2020年降至14.1%。在后期森林中，绿化伴随AR(1)上升的像元比例（8.9%）略高于伴随下降的比例（6.5%），而草地中两者比例相近。这说明绿化与恢复力变化的组合存在明显的时空异质性。

3.3 环境关联

环境模型对kNDVI的预测能力在全植被和草地中较高（交叉验证R²分别为0.640和0.563），森林中较低（0.361）；对AR(1)的预测能力则相对较低（全植被0.267，草地0.153，森林0.403）。SHAP归因表明，Pmean是全区和草地绿度的首要关联因子，Tmean是森林绿度的首要因子；AI则是三个植被域AR(1)的首要关联因子，Tvar和STmean也贡献显著。空间上，Pmean主导了83.4%面积的绿度解释；AI和Tvar分别主导了40.2%和37.0%面积的AR(1)解释。ALT与森林kNDVI呈显著正相关（R²=0.46，P=0.001），但与AR(1)无显著关系，表明活动层变化可能更多影响绿度而非时间持久性。

讨论部分指出，绿度与恢复力反映的是不同的生态系统属性：kNDVI主要描述冠层活动和绿度水平，AR(1)则刻画植被异常的持续性，用于推断恢复动态。二者对环境因子的响应不同，意味着仅用绿度评估生态系统状态会产生偏差。研究还强调，AR(1)是恢复力的间接代理指标，受外部强迫、随机变异和估计窗口等影响，解释时需保持谨慎；同时，CO₂、物种组成、干扰历史等未纳入模型，可能贡献了未解释的变异。

研究结论可概括为：2005—2020年青藏高原植被绿度显著上升，但推断恢复力未呈现出一致的改善。绿度与恢复力的环境驱动因子不同，活动层厚度与森林绿度相关，但与恢复力指标无关。因此，广泛绿化不应被解释为恢复力同步增强的证据；将绿度与恢复力相关的时间持久性相结合，才能为气候变化下青藏高原生态系统的评估提供更全面的依据。