干旱区生态系统中被忽视的灌丛结构恢复驱动非对称性恢复

《Earth's Future》：Structural Recovery of Overlooked Shrublands Drives Asymmetric Restoration in Dryland Ecosystems

　　

当前对干旱区生态系统的遥感监测 fundamentally 受限于对植被指数（"绿度"）的依赖，这类指数难以解混混合像元信号，且无法捕捉对恢复力至关重要的非光合结构组分。这种"绿度偏差"掩盖了稀疏植被的动态变化，导致在区分真正恢复与短暂绿化方面存在显著不确定性。研究人员开发了一种新型时空框架，将连续变化检测与分类（Continuous Change Detection and Classification, CCDC）与光谱混合分析（Spectral Mixture Analysis, SMA）相结合，以中国温带干旱区为例，解混亚像元尺度的光合植被（Photosynthetic Vegetation, PV）、非光合植被（Non-Photosynthetic Vegetation, NPV）和土壤组分比例。通过重建21年连续的端元丰度时间序列，研究人员绘制了此前未被检测到的稀疏灌丛分布，揭示其构成覆盖近半数植被区的优势"景观骨架"。研究结果揭示了生态系统转变的深刻非对称性：退化通常表现为植被-土壤耦合的突发性崩溃，而恢复则遵循独特的"结构优先"滞后路径。该恢复路径由NPV的逐渐累积所驱动，NPV稳定基质并促进后续冠层扩展。此外，轨迹分析表明这些恢复力机制严格依赖于基质条件，恢复路径在沙质、砾质和盐碱基质间存在显著差异。通过超越像元级绿度评估，实现亚像元结构表征，本研究为理解塑造干旱区的非对称过程提供了新见解，并为监测气候与人为压力下的适应性恢复提供了稳健的、空间显式的工具。

干旱区覆盖地球近40%的陆地表面，在全球生物地球化学循环和人类福祉方面具有根本性重要作用，但同时也面临着降水变率加剧、温度上升和干旱强化等前所未有的压力，加上不可持续的地下水开采、过度放牧和耕地扩张等人为干扰的叠加，使其面临潜在崩溃风险。尽管大规模生态恢复工程已广泛实施，但恢复结果高度异质，存在超出水文与土壤限制导致灌木死亡和次生退化的情况。这些对比鲜明的结果揭示了干旱区生态系统在退化与恢复之间转换的关键知识空白，以及绿度能否可靠反映生态系统恢复力的问题。

现有研究的核心局限在于，大多数大规模评估继续将绿度作为生态系统状况的代理指标，这种概念性混淆掩盖了观测到的绿化究竟反映真正的结构恢复还是仅反映短暂的功能响应，从而限制了对恢复阶段、恢复力阈值和状态转换的诊断能力。研究人员假设干旱区稀疏灌丛周围的恢复普遍遵循"结构优先"的演替过程，即非光合结构组分在强冠层绿化显现之前就已积累或稳定。为此，研究人员开发了基于轨迹的框架，首先通过光合植被、非光合植被和土壤基质的长期动态映射年际土地覆被，进而识别生态系统状态，最终重建代表性温带干旱区的状态转换长期轨迹。

研究区域选择中国西北部的河西走廊，该地区涵盖沙质沙丘、砾质戈壁、盐碱洼地以及广泛的稀疏灌丛和绿洲-荒漠生态交错带，为研究水分胁迫和缓慢植被恢复如何塑造生态系统功能提供了天然实验室。研究人员利用美国地质调查局（USGS）的Landsat Collection 2 Level-2 Tier 1地表反射率数据，时间跨度为2000至2020年，共计26,157景影像，通过FMASK算法质量控制波段去除云、云影和积雪影响，并剔除物理无效反射率值（<0或>1），预处理后在大部分区域保留了超过400个有效观测值，实现约每16天一个有效观测的高时间密度。

在技术方法层面，研究人员将传统的基底-植被-暗像元（Substrate–Vegetation–Dark, SVD）框架扩展为适用于异质干旱区的五端元方案：PV代表"快速"功能响应，NPV捕捉"缓慢"的结构生物量累积（如立枯凋落物、木质残体），沙质（Sand, SL）和盐碱（Saline, SA）组分表征驱动物理和化学反馈的不同土壤限制，暗色物质（Dark, DA）则解释阴影和土壤湿度变率。研究人员采用基于影像内在维度的标准化提取程序，利用主成分分析和凸几何理论从第一至第六主成分构成的特征空间顶点中识别潜在纯净像元，并结合ECOsis软件包的Fractional Cover Simulated VSWIR Data set Version 2作为NPV光谱参考源，以解决中等分辨率下纯净NPV端元难以一致识别的问题。为保持多年代际记录的时间一致性，分别针对Landsat TM/ETM+和OLI传感器建立了独立光谱库。完全约束线性光谱混合分析（Fully Constrained Linear Spectral Mixture Analysis, FCLS）在Google Earth Engine（GEE）平台上实现，模型平均均方根误差仅为0.0094，有效区分了绿洲（PV主导）、裸地（SL或SA主导）和稀疏灌丛（高NPV和DA分数特征）。

研究人员进一步应用CCDC-SMA算法对五端元分数的21年时间序列进行分割，使用谐波函数建模原始时间序列，当五个端元的观测值与模型值偏差连续超过预定阈值达8次观测时识别断点，确保检测到的变化代表持续性转变而非短期波动。模型生成描述长期趋势和季节模式的谐波参数，包括截距、斜率和前三项正弦/余弦系数。基于功能和基质框架，研究人员定义了三个功能层级：基底控制的荒漠、作为"生态桥梁"的稀疏灌丛（进一步分为沙质/砾质环境的旱生稀疏灌丛和盐碱/混合土壤的盐生稀疏灌丛）、以及高功能的高生产力密集植被。通过21年数据中的重复植被-基质组合推断六种基质类型（沙质、盐碱、砾质、沙-盐复合、盐-砾复合和混合基质），最终重建连续生态系统轨迹并分类为恢复轨迹（如SED）、退化轨迹（如SE）和复杂三状态轨迹（PR）。

研究结果首先验证了分类方法的可靠性，CCDC-SMA土地覆被分类平均总体精度达0.838，Kappa系数0.823，后处理提升至0.849和0.835。稀疏灌丛（SXS和SHS）占植被区面积的53.7%，而非此前认为的以密集冠层为主；SXS的NPV与PV振幅比值高达3.0，明确其为NPV主导型植被类型。2000至2020年间，SXS面积增加最大（+5,060.67 km²），反映出从荒漠向低覆盖度灌丛系统的结构性转变。研究区域内恢复轨迹（16,029.8 km²，64.3%）远多于退化轨迹（5,253.0 km²，21.1%），其中灌丛定植（Shrub Establishment in Desert, SED）是最普遍的轨迹类型（10,463.0 km²）。SED通常在初始恢复干预后滞后3-7年出现，并随距离现有斑块增加呈幂律衰减，表明存在强烈的局部促进效应和自组织恢复模式。

结果部分还揭示了显著的状态转换非对称性和基质依赖性。退化轨迹（如渐进退化PD、突发荒漠化AD）表现为植被-土壤耦合的急剧崩溃，PV振幅分别下降417和122，DA下降321，SA下降509；而恢复轨迹则是渐进的累积过程，渐进恢复（PR）的PV和NPV振幅增益最大（分别为+632和+533）。退化轨迹的PV振幅变化绝对值平均为恢复轨迹的1.37倍，NPV振幅差异达3.47倍，体现了"快速崩溃-缓慢恢复"的内在非对称性。基质因素强烈调制植被-土壤反馈：砾质基质上植被振幅增益最小且NPV甚至下降，而盐碱和沙质基质上PV和NPV均显著增加；沙质基质上PV和NPV下降最为陡峭，表明其对持续胁迫的高脆弱性。

在讨论部分，研究人员首先阐述了稀疏灌丛的功能作用与自我强化恢复机制。研究识别出覆盖研究区近三分之一植被面积的稀疏灌丛是当前全球评估中的"缺失状态"，约39,978 km²的稀疏灌丛被现有数据产品误分类为裸地。这些NPV主导的结构通过增加地表粗糙度减少风蚀、通过深根系机械锚固移动基质、通过"肥岛"效应维持生产力和促进逐步恢复，构成了干旱区恢复力的结构基础。其次，研究人员系统阐述了退化与恢复的非对称性和滞后动态：退化是阈值介导的不可逆过程，而恢复必须通过NPV积累的滞后路径逐步重建。突发崩溃（AD）表现为PV和NPV的同步急剧下降；渐进退化（SE）反映结构缓冲能力的逐渐耗竭；灌丛定植（SED）体现NPV持续上升而PV维持低位的关键滞后阶段；而促进型冠层扩展（SG）则是结构支撑下功能突破的实现。最后，研究人员强调了基质依赖的恢复策略：强制绿化的突发绿化（AG）轨迹在转型后呈现PV和NPV的负斜率，表明其不稳定性和地下水资源耗竭；而SG轨迹维持功能与结构组分的显著正斜率，证明稀疏灌丛是比密集造林更稳定的参照状态。沙质基质上的恢复遵循SED的机械稳定路径，砾质基质依赖生物土壤结皮的保护，盐碱基质则通过耐盐灌木的"肥岛"实现化学缓冲，这要求建立基质匹配的差异化恢复框架。

研究结论部分指出，有效的干旱区管理需要超越表层绿度，理解真正主导景观恢复力的结构机制。通过CCDC-SMA解混非光合与光合组分，本研究揭示干旱区恢复并非由快速绿化驱动，而是由缓慢的结构生物量累积所驱动，这种累积稳定了基质并促进了功能再生。研究结果凸显了覆盖近半数植被面积的稀疏灌丛的关键而常被忽视的作用——这些系统作为恢复轨迹中的优势"景观骨架"，在裸荒漠与稳定绿洲之间发挥着关键的生态桥梁功能。研究进一步证明，生态系统转变本质上是"结构优先"的滞后的，而非对称的：退化表现为功能的快速崩溃，恢复则需要NPV累积提供支持后续功能（PV）恢复的必要支架。此外，这些路径严格依赖于基质条件，意味着统一的造林措施在不兼容地形中往往失败。这些见解挑战了基于绿度的干旱区生产力传统评估范式，强调了稀疏灌丛等结构组分在增强生态系统应对气候变化和人为压力恢复力方面的重要性，为水受限环境下的有效恢复策略设计提供了宝贵指导。