气候变化与生态输水驱动下的塔克拉玛干沙漠绿化：水文?生态机制与生态风险

《Ecological Indicators》：Greening of the Taklamakan Desert: Effects of climate change and ecological water conveyance

　　

干旱区植被对气候变暖的响应存在高度异质性，塔克拉玛干沙漠绿洲?荒漠过渡带在极端干旱环境下出现的“绿化”现象，其水文?生态驱动机制尚缺乏多因子定量归因分析。研究人员以塔克拉玛干沙漠及周边绿洲系统为对象，整合多源遥感数据、野外实测资料、随机森林（RF）模型和地理探测器（GeoDetector），量化植被动态并分离气候变化、生态输水、地下水动态与土地利用变化的贡献。结果表明，在气候变暖与生态输水共同作用下，研究区绿洲呈现显著绿化趋势：自2000年以来，农田面积以约678 km²·yr^?1速率扩张，胡杨（Populus euphratica）林地和柽柳（Tamarix chinensis）灌丛分别以约389 km²·yr^?1和1579 km²·yr^?1速率增长，地表水域面积年均增加约68 km²，生态系统质量指数（EQI）和核归一化植被指数（kNDVI）均显著上升。驱动因子分析进一步揭示，冬季增温与地下水变化是植被绿化的首要驱动因子，而气候?水文因子交互作用进一步强化了植被响应；然而，局地复合干旱（compound drought）呈增强趋势，可能削弱增湿微弱的收益并威胁生态安全。同时，农田持续侵占天然河岸生境，乔木?灌木群落相互转化加速，暗示即使在水分供给改善条件下，乔木生境仍可能退化且群落长期稳定性下降。该研究定量揭示了气候变暖与生态输水背景下绿洲绿化的时空格局与驱动机制，可为生态输水调控、农业用水优化及荒漠河岸植被分区保育提供科学依据。

**论文解读**

**一、研究背景与问题提出**

全球变暖背景下，干旱半干旱区陆地生态系统正承受日益加剧的热胁迫、水文波动和持续性干旱的多重压力。塔克拉玛干沙漠位于中国西北塔里木盆地，四周高山环绕，年降水量不足50 mm，而年潜在蒸发量超过2500 mm，是世界上生态系统最脆弱的区域之一。该区域水资源几乎完全依赖周边高山冰雪融水补给的内陆河流，形成“山地?绿洲?荒漠”的独特景观格局。荒漠河岸林（desert riparian forests）对水分波动极度敏感，抗干扰能力弱，自我恢复能力有限；20世纪下半叶，上游过度引水叠加持续干旱曾导致塔里木河下游断流、地下水大幅下降、胡杨林大面积衰亡，绿洲生态系统濒临崩溃。为此，中国政府于2000年启动生态输水工程（Ecological Water Conveyance Project），向下游干涸河道周期性输水以恢复河道及沿岸绿洲。

然而，现有研究存在两方面明显不足：一是缺乏多因子定量归因分析，气候变化与生态输水等关键驱动因子的各自贡献尚未被清晰量化；二是受限于长时序、区域尺度高分辨率动态监测数据的匮乏，对绿洲演替及其驱动机制的认识仍不完整。基于此，该项研究聚焦核心科学问题：在气候变化和生态输水双重作用下，塔克拉玛干沙漠绿洲“绿化”现象的水文?生态机制是什么？绿化趋势是否伴随着潜在的生态风险？

**二、研究内容与重要意义**

研究人员整合多源遥感数据、野外实测资料、随机森林模型与地理探测器分析，系统量化了2000年以来研究区绿洲植被、农田扩张、地表水动态和地下水条件的长时序变化，识别了农田与天然植被绿化的主导气候与水文驱动因子，重点考察了地下水变化、冷季增温和干旱过程的作用，并评估了绿洲绿化是否伴随农田侵占、河岸植被结构转变和生态系统稳定性下降等潜在生态风险。该研究为理解极端干旱区绿洲生态系统对气候变化与人为水资源调控的综合响应机制提供了重要科学依据，对生态输水调控、农业用水优化和荒漠河岸植被分区保育具有直接的指导意义。

**三、关键技术与方法概述**

在技术方法上，研究人员主要依托以下手段：（1）基于Landsat 5/7/8无缝影像合成数据，利用随机森林算法进行土地覆盖分类（land?cover mapping），以量化农田、乔木和灌丛的多年来际变化；（2）基于地质云平台国家地下水监测数据与MODIS影像，构建随机森林回归模型反演地下水埋深，并通过2018?2022年训练、2023?2025年独立检验评估模型的时间迁移能力；（3）利用30 m分辨率的河流与湖泊地表面积数据集（SARL）分析地表水域动态；（4）基于ERA5?Land数据构建复合干旱指标——将30天累积降水转化为标准化降水干旱指数（SPI），将5天滑动平均温度转化为标准化热力指数（SHI），两者结合为复合日热力干旱指数（CDHI，compound daily heat index），以阈值法识别复合干旱事件；（5）利用1 km逐日全天候土壤水分数据集，以五日均值百分位法识别骤发干旱（flash drought），并仅保留持续6?18个候（约30?90天）的事件；（6）基于MODIS产品构建生态系统质量指数（EQI），综合叶面积指数（LAI）、植被覆盖度（FVC）和总初级生产力（GPP）三者归一化等权合成；同时计算核归一化植被指数（kNDVI）作为补充指标；（7）采用Mann?Kendall趋势检验与Sen斜率结合，分析各要素时空变化趋势及其显著性；（8）运用地理探测器（GeoDetector）量化单因子及双因子交互对植被变化空间分异的解释力。所有变量分析前统一重采样至500 m空间分辨率并对齐投影与网格系统，数据时段为2003?2022年。

**四、研究结果**

（一）土地覆盖分类与变化

基于Landsat影像和随机森林分类器，研究人员获得了2000?2024年研究区农田、灌丛和乔木覆盖（林地）的分类图，模型决定系数R²为0.8349，Kappa系数为0.8095。年际变化分析显示，农田以678.07 km²·yr^?1的速率扩张，但增速在2015年前后开始放缓。天然植被中，胡杨代表的乔木群落面积以389.04 km²·yr^?1的速率增加，柽柳为主的灌丛群落面积年增1578.94 km²。地表水域面积以68.3826 km²·yr^?1的速率扩大。与此同时，自2000年以来，农田扩张分别挤占了约8256.503 km²的乔木林地和3996.711 km²的灌丛。群落间转变方面，乔木转为灌丛的面积以年均70.3972 km²增加，灌丛转为乔木的面积以年均75.65 km²增加。

（二）地下水动态

结合MODIS影像建立的随机森林回归反演模型，在独立时间验证期（2023?2025年）获得R² = 0.73、RMSE（均方根误差）= 4.82 m的性能。2001?2024年地下水反演结果显示，研究区地下水位整体呈上升趋势，平均每年上升约0.07 m。

（三）气候变化

骤发干旱方面，2003?2024年塔克拉玛干沙漠骤发干旱发生次数总体呈减少趋势。气温方面，2000?2024年年均温升幅为0.0289 °C·yr^?1，夏季升幅0.019 °C·yr^?1，冬季升幅高达0.0699 °C·yr^?1，表明增温主要由冬季驱动；空间上，总体增温在东北部较温和，其余区域更显著，夏季增温西强东弱，冬季增温自西向东减弱。降水方面，年均降水增加0.02 mm·yr^?1，夏季增加0.065 mm·yr^?1，冬季微减0.0012 mm·yr^?1；夏季降水在西南部增加较明显，冬季降水在中部区域上升而西、东部下降。复合干旱方面，综合多指标年际趋势显示，复合干旱整体呈加剧态势，加剧区域主要位于研究区西部和西南部。

（四）生态系统变化

2001?2024年EQI和kNDVI的MK趋势分析一致表明，塔克拉玛干沙漠生态系统状况总体改善。进一步分析复合干旱趋势与kNDVI变化的空间关系发现，在复合干旱加剧的区域植被绿化并未被全面压制，部分绿洲与河岸带的kNDVI仍呈正向趋势，提示地下水恢复、生态输水和土地利用相关的绿洲扩张可能在局地对植被起到缓冲复合干旱胁迫的作用。

（五）驱动因子分析

基于更高分辨率（30 m）的Landsat?kNDVI产品，随机森林模型结果显示：对农田kNDVI而言，冬季温度变化和年均温变化影响力最大，其次为地下水变化；对天然植被kNDVI而言，年均温变化、冬季温度变化和地下水变化亦为最主要预测因子，表明温度相关因子和地下水动态是制约长时序绿度趋势的关键约束条件。地理探测器分析进一步显示，地下水变化与冬季温度的单因子解释力（q值）相对最高，而双因子交互分析表明，冬季温度与地下水变化、地下水变化与复合干旱总天数、年均温变化与地下水变化等因子组合的q值约为0.05?0.07，呈现非线性增强或双变量增强效应，说明气候?水文多因子耦合对植被健康趋势的空间分异具有更强的联合解释力。

**五、讨论与结论**

讨论部分指出，研究区2000?2024年整体绿化与气候增暖同步，但冬季增暖尤为突出，达到0.699 °C·decade^?1，超过近期全球平均增温率和西北干旱区已有报道的增温率，凸显了内陆干旱绿洲生态系统的高气候敏感性。冬季增暖可能通过减少土壤冻结深度、促进春季提前解冻和融雪、改变物候期等机制，改善旱春水热条件，延长有效生长季，这可以解释冬季温度在驱动分析中的重要性；然而，若早季蒸发需求上升而地下水或地表水补给不足，也可能加剧后期水分胁迫。复合干旱加剧意味着高温与降水亏缺的同步冲击可能通过提升潜在蒸散发、增大植被蒸腾需求和加剧本本已脆弱的光合?水力系统热力风险，削弱微弱增湿的正面效应。

在生态输水与土地利用影响方面，地表水域扩大和地下水位上升需要放在区域水平衡背景下理解：塔里木盆地极度干旱，局地降水远不足以支撑大规模农田扩张，绿洲农业主要依赖上游来水（冰雪融水）、渠系引水、水库调蓄、生态输水和灌溉回归水等调控输入。因此，地下水位上升不能简单视为降水驱动的自然“变湿”信号，而更可能是生态输水、渠系调控和灌溉渗漏等综合作用的结果。农田持续侵占河岸生境、乔木?灌丛相互转化加速的现象进一步表明，绿洲面积扩张并不自动等价于生态质量的同比例改善；水分供给增加初期可促进植被生长，但长期积水可能诱发缺氧胁迫、盐分再分布和土壤次生盐渍化，危及成熟乔木生存，而灌木凭借更快的资源利用策略在竞争中可能占据优势，进一步压缩乔木更新空间，导致群落功能退化和稳定性下降。

研究结论明确指出：在气候变暖与生态输水共同作用下，塔里木河下游河岸绿洲呈现显著绿化趋势。自2000年以来，农田以约678 km²·yr^?1速率扩张，胡杨林和柽柳灌丛面积分别以约389 km²·yr^?1和1579 km²·yr^?1速率增加，地表水域面积年均扩大约68 km²，EQI和kNDVI均显著上升，反映区域生态系统状况整体改善。冬季增温与地下水变化是植被绿化的首要驱动因子，气候与水文因子的交互效应进一步强化了植被响应。然而，局地干旱风险依然存在，绿洲仍面临潜在水分胁迫；农田对天然河岸生境的持续侵占以及种间竞争的加剧，已导致植被结构退化，暗示即使在增水条件下乔木生境仍可能恶化，群落长期稳定性趋于下降。该研究从时空格局与驱动机制层面为生态输水调控、农业水资源优化和荒漠河岸植被分区保护提供了定量科学支撑。