高耦合但低协调：快速绿化下中国黄土高原气候-植被相互作用

《Ecological Indicators》：High coupling but low coordination: Climate–vegetation interactions on the Chinese Loess Plateau under rapid greening

　　

在干旱和半干旱地区，大规模植被恢复被广泛用于防治土地退化。然而，人为导致的绿化是否与气候动态同步以实现长期可持续性仍不清楚。中国黄土高原作为一个重要的生态恢复区域，提供了理想的案例研究。研究人员采用耦合协调度（CCD）模型和1990—2020年的多源数据集，研究了气候与植被之间的时空相互作用。结果显示，研究区呈现显著的暖湿化趋势，并伴随大量绿化和碳储量增加。一个明显的错配出现：系统表现出高耦合但低协调（耦合度>0.80，协调指数<0.15），表明存在持续的轻度失调（耦合协调度<0.35）。这表明旱地恢复中可能存在系统性风险，即人为驱动的植被扩张可能超过气候反馈。主导驱动因子从自然生产力转变为土地利用变化，并最近转变为降水约束。这些结果表明，水资源可用性是绿化的关键约束，植被覆盖增加并不必然保证可持续性。因此，生态恢复政策可能需要从以绿化为主导的方法转向以综合评估为导向的框架，考虑热点地区的恢复质量、不显著地区的生态系统稳定性和冷点地区的水文气候适宜性，以支持区域生态恢复的长期可持续性。

**研究背景与意义**

气候变化加剧全球水文循环，改变降水格局，尤其对干旱和半干旱地区造成更大生态压力。在水资源有限的生态系统中，气候变异与植被动态的相互作用是决定生态系统恢复力与可持续性的核心因素。大规模生态恢复工程旨在逆转土地退化，但人类诱发的绿化是否与气候动态同步演化，仍是全球性挑战。中国黄土高原（CLP）拥有高度可蚀的黄土和长期人为干扰历史，是全球研究生态脆弱性和恢复力的理想区域。自1999年实施退耕还林（Grain for Green）工程以来，该区域经历了显著绿化。然而，快速植被扩张也引发对生态水文后果和长期可持续性的担忧。以往研究多聚焦气候对植被的单向影响或造林对水文的影响，缺乏对气候-植被子系统间协同与协调发展的定量评估。因此，研究人员引入耦合协调度（CCD）模型，旨在量化CLP气候-植被系统的时空演变，识别驱动机制，为大规模恢复的可持续性提供适应性管理建议。

**研究方法概览**

研究人员收集1990—2020年多源数据集，包括气候因子（降水、温度、日照、蒸发、土壤水分）、归一化植被指数（NDVI）、土地利用变化和碳储量。空间数据统一投影到Albers等积圆锥投影，重采样至1 km分辨率。采用InVEST模型“碳储存与固持”模块估算碳储量；运用空间自相关和热点分析评估NDVI空间集聚；构建包含气候和植被子系统的综合评价指标体系，结合熵权法确定权重，计算耦合度（C）、协调指数（T）和耦合协调度（CCD）。使用ArcGIS 10.8进行空间分析，Python 3.8执行CCD建模，Origin 2024b可视化，SPSS 25进行显著性检验。地理探测器模型用于识别驱动因子及其交互作用。

**研究结果**

*3.1 气候因子变化*：1990—2020年，CLP呈现显著暖湿化趋势。年降水量（PRE）从438.7±130.4 mm增至507.5±183.3 mm；年均气温（TEM）从8.0±2.9°C升至8.5±2.9°C；表层土壤水分（SWC）从20.4%±3.8%恢复至21.6%±3.9%；蒸发量（EVP）持续下降；日照时数（SSD）先降后升。空间上，PRE、TEM和SWC自西北向东南增加，SSD和EVP相反。

*3.2 NDVI与土地利用变化*：NDVI在30年间显著上升（R²=0.8），但空间异质性持续存在，约60%区域仍为低中覆盖（NDVI<0.4）。NDVI呈显著集聚格局（Moran's I>0.87），热点和冷点面积分别由30.3%和38.1%降至28.1%和30.1%，不显著面积由31.3%增至41.9%。土地利用转换以耕地转为林地和草地为主，以及未利用地生态恢复。净变化为耕地减少21,500 km²、未利用地减少15,800 km²，林地增加17,300 km²、草地增加7,700 km²、水域增加600 km²、建设用地增加11,600 km²。

*3.3 碳储量变化*：区域碳储量（CS）从1990年的3.95×10⁹ t增至2020年的4.14×10⁹ t。草地和林地CS持续增加，耕地CS下降，与退耕还林直接相关。空间上呈现中高CS区扩张、低CS区收缩的趋势；CS热点区先减后增，冷点区持续收缩，不显著区逐步增加。

*3.4 耦合协调分析*：耦合度（C）保持在0.80—0.86高位，表明气候与植被子系统相互作用强烈；但协调指数（T）仅0.11—0.15，导致CCD处于0.28—0.35，属轻度失调。空间上，Iveg、T和CCD均呈东南高、西北低，C呈相反格局。所有指标排序为热点区>不显著区>冷点区。敏感性分析（±10%权重扰动和等权重方案）确认结果稳健。

*3.5 驱动机制分析*：地理探测器显示，PRE、NDVI、土地利用类型（LUT）、NPP和CS是解释CCD空间变异的持续重要因子，LUT相对重要性在2000年后显著增加。分区域看，热点区主要受LUT（0.44）、CS（0.43）、NDVI（0.41）驱动；不显著区以NPP（0.31）、PRE（0.25）为主；冷点区则以EVP（0.24）、PRE（0.12）为主。双因子交互作用普遍强于单因子，且与PRE、LUT、NDVI和CS的组合增强效应明显。

**讨论总结**

讨论部分指出，CLP植被恢复受气候变异和人类活动共同调节，但两者作用具有阶段性和空间异质性。尽管暖湿化有利于植被生长，降水是NDVI增加的最直接气候驱动因子，但在低覆盖区，降水增加可能加剧土壤侵蚀，不必然促进植被。NDVI与PRE、TEM、SWC正相关，与SSD在年际尺度负相关。人类活动驱动的土地利用变化加速绿化，但也加剧景观破碎化，降低植被对气候压力的恢复力。高耦合低协调状态反映气候反馈滞后和人为干预主导，植被恢复可能已超出水资源承载力。时间上主导因子呈“自然生产力→土地利用变化→降水约束”的三阶段演变，表明水限制日益突出。空间上热点区相对稳定，冷点区受蒸发和土壤水分制约，河谷地区受城市扩张影响，不显著区呈现复杂风险。政策上建议从绿化导向转向协调发展导向，构建包含气候适宜性、水资源、碳汇、土地利用结构和植被生长的综合评价体系；热点区维持森林质量，不显著区优化土地利用和加强水土保持，冷点区选用耐旱节水乡土植被。研究也存在局限：InVEST模型依赖土地利用类别，未考虑物种和年龄差异；年尺度分析无法捕捉季节动态与极端事件；遥感数据存在不确定性。

**研究结论**

中国黄土高原经历了明显的暖湿化趋势，表现为年降水量（PRE）和年均气温（TEM）增加，年蒸发量（EVP）下降。大规模生态恢复和土地利用转变增加了CLP绿化过程中的空间异质性。气候-植被系统在整个研究期表现出典型的高耦合但低协调模式。这一模式表明，植被恢复依赖人为干预而非自然气候驱动的演替，使生态系统极易受到极端气候事件或人为维护减少的影响，尤其是在缺水的西北冷点区。系统协调的主导驱动因子呈现明显的三阶段转变：从1990年的自然生产力指标，到2000—2010年的土地利用变化指标，最后到2020年的水分约束指标。空间上，协调水平呈现东南-西北梯度，热点区保持相对稳定，不显著区面临复杂生态风险，冷点区受刚性水分限制约束。总体而言，CLP植被恢复的长期可持续性最终取决于气候条件，特别是水分可用性。因此，建议生态恢复工作重点关注恢复林质量、平衡人地关系、优化植被和土地利用结构、增强水土保持功能。将耦合协调度纳入生态恢复评估可进一步支持差异化管理，促进区域生态恢复的长期可持续性。然而，仍存在一些不确定性：遥感产品和模型模拟可能引入数据和参数误差，而年尺度的气候-植被耦合分析可能限制对季节动态和极端气候影响的表征。未来研究应结合实地观测以完善气候-植被耦合协调度评估。