中国草地生态补偿政策下中亚山地草原的绿化与稳定性响应分化

《Land Degradation & Development》：Greening and Divergent Stability Responses in Central Asian Mountain Grasslands Under the Grassland Ecological Compensation Policy of China

　　

摘要：中亚山地草原被公认为全球重要的生物多样性热点和关键牧业资源，目前正经历由气候变化和放牧实践演变驱动的深刻转型。本研究探讨了在天山中部巴音布鲁克草原，两个对比时间段——中国草地生态补偿政策实施前（2000–2010年）和实施后（2011–2024年）——生态系统稳定性的多维动态。研究采用了准时间观测设计，整合了25年（25?years）卫星衍生植被指数与野外监测数据。结果如下：（1）研究区域72.7%的面积观测到绿化趋势，主要在2011年后；此外，尽管整体恒定性（constancy）增加，但大多数牧场类型的抵抗力（resistance）下降，春季和冬季牧场的恒定性增强更显著，夏季和秋季牧场则出现更多褐化信号，表明绿化和时间恒定性的改善并不必然赋予抵御极端气候事件的能力增强；（2）气候变量仅解释年际生产力变异的3.1%，且2011–2024年期间的气候调整残差（climate-adjusted residuals）持续高于2000–2010年，表明非气候因素（可能包括与经济社会变化相关的放牧压力重新分配以及政策引导的季节性转场重组）的作用；（3）生产力和恒定性沿海拔和水热梯度呈单峰模式，而抵抗力呈现相反响应；此外，恒定性-抵抗力相关性在中度环境条件下达到峰值。这些结果强调，多维稳定性评估（而非仅依赖绿化趋势）为评价草地管理成效提供了更全面的框架；研究结果还凸显了在持续气候变化背景下，需要采用季节特异性监测和放牧策略来保护高寒草原的功能性。

**论文解读：中国草地生态补偿政策下中亚山地草原的绿化与稳定性响应分化**

**研究背景与科学问题**

中亚山地草原是全球重要的生物多样性热点和牧业资源，但近年来气候变化与放牧管理实践变革正深刻影响其生态系统功能。传统上，温度曾是限制高山草原生产力的主要因素，但随着气温升高，水分可利用性因蒸发需求增加和降水不确定性而日益受限。在此背景下，植被绿化趋势并不必然代表生态系统稳定性的实质性改善，因为气候变化不仅影响平均生产力，还改变其年际变异。生态系统稳定性是一个多维概念，包括生产力（productivity，反映植被功能平均水平）、恒定性（constancy，描述生产力围绕长期均值的时序稳定性）和抵抗力（resistance，生态系统在极端逆境条件下维持功能的能力）。这些维度可能呈现非平行响应，其两两相关性沿环境梯度系统变化，对理解生态系统对干扰的响应至关重要。中国于2011年实施的草地生态补偿政策（Grassland Ecological Compensation Policy, GECP）在八个省份推行放牧改革，减少载畜量并限制某些放牧方式，同时社会经济变化促使牧民缩短迁徙路线，重塑了季节性牧场的利用时空格局。尽管GECP实施后植被绿化已有报道，但绿化是否对应生态系统稳定性多维度的协同变化，以及这些变化是否因季节牧场类型和环境梯度而异，仍不清楚。天山是中亚最大山脉，目前关于该山脉草地稳定性的研究相对较少，且多聚焦生产力对气候的响应，忽视了季节性牧场系统的管理现实。巴音布鲁克草原位于天山中部，具有显著的海拔梯度和完善的季节性牧场系统，是检验上述问题的理想场所。本研究利用25年（2000–2024年）卫星核归一化植被指数（kernel Normalized Difference Vegetation Index, kNDVI）数据与野外监测数据，比较GECP前（2000–2010年）和后（2011–2024年）两个时期的生产力、恒定性和抵抗力，旨在回答三个关键问题：（1）生产力、恒定性和抵抗力在两个时期是否一致变化，且变化是否超出气候趋势所能解释的范围？（2）这些响应用于不同季节牧场类型是否一致？（3）稳定性成分间的两两相关性如何沿海拔和水热梯度变化？

**主要技术方法**

研究人员采用准时间观测设计，整合多源数据：MODIS MOD13Q1 NDVI产品（250?m分辨率，16天合成）经kNDVI非线性泛化计算生产力，并与巴音布鲁克草地生态系统研究站4个固定样地（样地来源：自1984年建立，2004–2024年每年5月中旬至8月下旬采集地上生物量，共881个有效样本）的实测生物量验证（R2=0.61, p<0.001）。气候数据来自国家地球系统科学数据中心（月降水P和月均温T，2000–2024年），SRTM高程数据（90?m分辨率），季节性牧场边界由和静县林业和草原局提供（2024年矢量数据）。数据分析包括：Theil-Sen斜率估计与Mann-Kendall检验检测kNDVI趋势；单因素和双因素ANOVA比较不同时期和牧场类型间的稳定性指标；偏相关分析、气候趋势分析和气候调整残差分析评估气候与非气候因素贡献；皮尔逊相关与Loss曲线模型分析稳定性成分间关系；分段线性回归揭示沿环境梯度（海拔100?m、温度1°C、降水100?mm间隔）的变化。

**研究结果**

**3.1 kNDVI的时空趋势**：2000–2024年研究区年均kNDVI总体呈增加趋势，21.4%的草地显著改善，仅1.4%显著下降。2011年后绿化信号增强，改善面积从2000–2010年的31.3%升至2011–2024年的82.6%，褐化面积相应从68.7%降至17.4%。该结果表明GECP实施后绿化趋势显著，但需结合其他稳定性维度解读。

**3.2 生产力、恒定性和抵抗力在季节性牧场类型间的变异**：空间上，生产力和恒定性呈北高南低梯度，抵抗力呈斑块分布，高值集中于北部和西部。不同牧场类型存在基线差异：秋季和冬季牧场生产力和恒定性较高，春季和春夏秋牧场较低，生长期牧场抵抗力最高。2011–2024年期间，恒定性增加而抵抗力下降，恒定性的增加仅出现在春季和冬季牧场，秋季和夏季牧场出现更多褐化信号。双因素ANOVA显示牧场类型是生产力、恒定性和抵抗力的主要变异来源（所有p<0.001），而时期效应仅对恒定性显著（F=5.041, p<0.05），时期×牧场类型交互作用不显著（抵抗力边缘显著p=0.085），表明时间变化在牧场类型间大致一致。

**3.3 生产力、恒定性和抵抗力的气候驱动因素**：2000–2024年温度显著上升（R2=0.25, p<0.05），降水趋势不显著（R2=0.01, p>0.05）。偏相关分析显示生产力与温度、降水均无显著关系。包含温度和降水的多元线性回归仅解释生产力年际变异的3.1%。气候调整残差在整个研究区和各季节牧场类型中，2011–2024年期间均高于2000–2010年，其中秋季和冬季牧场增加显著。这表明植被变化不能完全归因于气候，非气候因素（如放牧压力重新分布）可能起重要作用。

**3.4 生产力、恒定性和抵抗力沿环境梯度间的关系**：区域尺度上，生产力与恒定性正相关，与抵抗力负相关，恒定性-抵抗力也为负相关，且2011–2024年期间相关性更强。沿环境梯度，生产力和恒定性随海拔、温度、降水呈单峰（驼峰）模式，峰值分别出现在2921?m和2998?m、7.63°C和7.44°C、264?mm和273?mm；抵抗力则先降后升，谷值在3426?m、5.41°C、304?mm。两时期峰值位置略有偏移。恒定性-生产力相关性在海拔3575?m、温度7.27°C处最强，恒定性-抵抗力相关性在海拔3131?m、温度7.19°C处最强，且2011–2024年期间相关性普遍弱于2000–2010年。

**讨论与结论**

讨论部分指出，生产力和恒定性沿南北梯度与区域水热约束一致，而抵抗力受微地形控制更明显。尽管绿化趋势明显，但恒定性和抵抗力呈现分离响应：恒定性提高反映年际变异减小，但抵抗力下降表明生态系统抵御极端事件能力减弱，这与物种多样性降低和群落组成变化有关。统一的强迫气候不能解释牧场类型间的季节不对称性，夏季和秋季牧场褐化更明显、春季和冬季牧场恒定性增强的可能原因包括：社会经济变化导致牧民缩短迁徙路线，集中放牧于低海拔暖季牧场，而政策引导的暖季放牧结合冷季舍饲进一步延长暖季停留时间，增加夏季和秋季牧场局部放牧压力。但研究人员未直接测量载畜量或牲畜移动，因此放牧再分布仅为合理推测。

研究结论翻译：本研究考察了巴音布鲁克草原在GECP前（2000–2010年）和GECP后（2011–2024年）期间生产力、恒定性和抵抗力的动态。显著绿化主要发生在2011年之后；然而，三个指标呈现分离响应：大多数牧场类型生产力和恒定性增加，而抵抗力下降；这一发现表明绿化并不必然反映整体稳定性的改善。气候变量仅解释生产力年际变异的3.1%，且2011–2024年期间气候调整残差持续较高，表明非气候因素可能对观测变化有贡献。季节不对称性明显，春季和冬季牧场恒定性增强更显著，夏季和秋季牧场褐化信号更多，这与暖季放牧压力可能重新分布的观察一致。生产力和恒定性沿环境梯度在中度条件下达到峰值，而抵抗力呈现相反模式。这些发现强调了在持续气候和土地利用变化下，采用季节特异性和多维方法监测和管理山地草原的重要性。