塔里木河下游胡杨林对生态输水的地表淹没梯度与滞后植被覆盖度响应

《Forests》：Surface Inundation Gradients and Lagged Fractional Vegetation Cover Responses of Populus euphratica Oliv. Forests to Ecological Water Conveyance in the Lower Tarim River, China

　　

摘要翻译：为量化生态输水与胡杨（*Populus euphratica* Oliv.）林植被覆盖度（FVC）之间的空间及滞后关联，研究人员整合了1995–2025年的Landsat影像、输水记录、地下水观测数据以及多期随机森林分类结果，研究区位于塔里木河下游。在动态年度响应窗口内识别了新形成或扩展的地表水体，并以年度地表淹没扩展面积、多年淹没频率以及森林像元至反复淹没水面的距离为特征参量。上述指标采用Theil–Sen趋势及Lag 0–Lag 3相关分析与FVC建立关联。2010年之后，水库泄水量、年度地表淹没扩展面积及地下水状况总体趋于改善。平均FVC以每年0.182个百分点的速率增加，且改善与退化在空间上交错分布。相对低至中等频率的淹没与显著的FVC正向趋势相关联，FVC改善集中在距反复淹没水面0–0.5 km范围内。FVC与水库泄水量、地下水埋深及年度地表淹没扩展面积的关联在滞后2–3年时最强。综上，这些地表淹没指标补充了水库泄水量与地下水记录，并支持对胡杨林FVC变化开展空间差异化监测。

论文解读：

研究背景与科学问题：塔里木河下游地处塔克拉玛干沙漠与库鲁克沙漠之间，是维护绿洲稳定与区域生态安全的重要生态屏障。胡杨（*Populus euphratica* Oliv.）作为该区域优势河岸树种，其生长主要依赖河川径流、地下水以及间歇性地表水补给。自2000年起实施的生态输水工程促进了河道水文恢复与胡杨生境改善，但输水量本身无法完全解释森林恢复的空间异质性。已有评价多结合输水量、水文响应与植被恢复指标，但水库泄水量和输水历时仅反映工程输水强度，无法表征水体是否抵达侧向渠道、漫滩或林分斑块；地下水埋深仅刻画监测点的地下水条件；植被覆盖度（FVC）和归一化差异植被指数（NDVI）虽能反映植被结果，却未定位地表水的空间可达性。因此，输水地表水的空间扩展范围、重现频率及其与胡杨林的距离如何影响森林FVC的空间分异，成为亟待解决的关键问题。

研究内容与结论：研究人员基于1995–2025年Landsat长时序影像，结合水库泄水量、生态输水日期、地下水埋深观测和多期胡杨林分类结果，构建了多时相共识胡杨林掩膜，提取了输水响应期内新形成或扩展的地表水体，计算了年度地表淹没扩展面积、多年淹没频率以及森林像元至反复淹没水面的距离，并与胡杨林FVC进行Theil–Sen趋势分析和Lag 0–Lag 3滞后相关分析。结果表明：2010年后水库泄水量、年度地表淹没扩展面积和地下水条件总体改善，胡杨林年平均FVC以每年0.182个百分点的速率显著上升，但改善与退化在空间上交错分布；相对低至中等频率的淹没区域表现出更明显的FVC改善，且改善集中在距反复淹没水面0–0.5 km范围内；FVC与水库泄水量、地下水埋深及年度地表淹没扩展面积的关联在滞后2–3年时最强。这些发现表明，地表淹没扩展面积、淹没频率和距离指标分别从扩展范围、重现频率和空间邻近度三个维度表征了地表水可达性，补充了泄水量与地下水记录，为胡杨林恢复的长期监测与输水管理提供了空间差异化依据。

关键技术方法：研究使用了Landsat 5/7/8/9地表反射率影像，在Google Earth Engine（GEE）平台上进行云与云影去除、反射率定标、波段协调及跨传感器一致性校正。采用随机森林（Random Forest, RF）分类器提取1995、2000、2005、2010、2015、2020和2025年共七期土地覆盖类型，通过≥3期重复识别原则构建多时相共识胡杨林掩膜；利用修正土壤调节植被指数（MSAVI）和像元二分模型反演逐年FVC；基于改进归一化差异水体指数（MNDWI）、NDVI、近红外、短波红外1和亮度指数综合阈值法提取新形成或扩展的地表水体；使用Theil–Sen中值斜率与Mann–Kendall检验分析FVC趋势；通过Kruskal–Wallis检验、Mann–Whitney U检验及Benjamini–Hochberg错误发现率（FDR）校正比较不同淹没频率类别间的差异，并采用空间分块聚类自助法和全局Moran’s I评估空间稳健性。样本来源包括野外样点、Google Earth历史影像和Landsat目视解译；地下水数据来自英苏、喀尔达伊和依干不及玛三个典型监测断面。

研究结果：

1. 生态输水驱动的水文过程与地表淹没格局：2000–2009年水库泄水量总体偏低且年际变率大；2010–2020年平均泄水量从2.27×10⁸ m³增至5.62×10⁸ m³，三个监测断面平均地下水埋深从5.71 m降至4.55 m；2021–2025年泄水量保持中高水平。年度地表淹没扩展面积在2000–2009年平均约193.3 km²，2010–2020年增至451.0 km²，2017年达峰值820.3 km²，2021–2025年降至342.2 km²。淹没频率呈明显河道约束格局，低频率和相对低频率淹没面积占54.1%和18.7%，相对高频率和高频率仅占6.7%和2.0%；检测到的地表淹没与487.7 km²的胡杨林分析区重叠，占18.2%，其中88.5%集中于低至中等频率类别。

2. 胡杨林FVC变化特征：七期随机森林分类的平均总体精度（OA）为88.0%，平均Kappa系数为0.840；胡杨林生产者精度（PA）和用户精度（UA）均值分别为78.4%和76.4%。1995–2025年胡杨林年均FVC多年均值为37.29%，Theil–Sen斜率为0.182个百分点/年，Mann–Kendall检验显示显著上升趋势（p = 0.041）。空间上，FVC呈沿河条带状和斑块状分布，以30%–60%的中低值为主。FVC改善区占41.8%，退化区占38.3%，无明显变化区占19.9%，改善与退化交错分布。

3. 地表淹没模式与胡杨林生长响应：不同淹没频率类别间FVC增长差异显著。无淹没区净改善指数为?1.7个百分点；相对低频率淹没区改善比例最高（63.0%），净改善指数为42.7个百分点，平均Sen斜率最大（0.77个百分点/年）；高频率淹没区改善比例降至47.2%，平均Sen斜率仅为0.07个百分点/年。沿距离梯度，FVC改善主要集中在0–0.5 km范围内，该距离带内改善像元占53.7%，退化像元占29.4%，净改善指数为+24.3个百分点；超过0.5 km后改善比例下降、退化比例上升。

4. 滞后相关分析：FVC与水库泄水量在滞后2–3年显著正相关，与地下水埋深在滞后0–3年显著负相关（2年滞后时最强），与年度地表淹没扩展面积在滞后2–3年显著正相关。在最优滞后下，年度地表淹没扩展面积（Lag 3）的R²为0.433，水库泄水量（Lag 3）为0.281，地下水埋深（Lag 2）为0.449。Spearman和去趋势分析支持水库泄水量Lag 2、地下水埋深Lag 1–2、地表淹没扩展面积Lag 2–3的结果。

讨论总结：研究指出，年度地表淹没扩展面积与水库泄水量并非完全同步，表明输水量变化并未直接转化为比例性的侧向水体扩展。相对低至中等频率的淹没更有利于FVC改善，而最高频率类别的改善较弱，可能与高频率淹没区域的土壤盐分和立地条件差异有关。距离梯度显示，距反复淹没水面越近，地表水可达性越高，浅层入渗和局部地下水抬升可能增加根区水分有效性。滞后2–3年的相关模式说明水文变化与FVC动态之间存在紧密的时间关联，年际气候变率可能构成额外背景影响。研究建议，对远离反复淹没水面且FVC下降的区域优先开展野外验证；对确认水文失联的区域可评估侧向渠道或漫滩重连、清除局部阻水障碍、浅沟连通及调整输水历时等措施。相对低至中等淹没频率可作为间歇性输水的参考条件；0–0.5 km响应带、低至中等淹没频率及2–3年滞后为未来监测提供了经验参考范围。

研究结论：研究人员基于1995–2025年Landsat长时序影像，结合年度地表淹没扩展面积、多年淹没频率和森林像元至反复淹没水面的距离，评价了塔里木河下游胡杨林对生态输水的响应。主要结论如下：（1）2000–2025年生态输水响应期内，新形成或扩展的地表水体呈现明显的时空变化，2010年后水库泄水量和年度地表淹没扩展面积总体增加，但年际变化并非完全同步；多年高频率淹没区主要沿主河道、连通侧向渠道和低洼尾闾区分布。（2）胡杨林FVC在1995–2025年呈总体波动上升趋势，多年均值为37.29%，Theil–Sen斜率为0.182个百分点/年，上升趋势显著；像元尺度上改善、稳定和退化并存，空间异质性明显。（3）地表淹没频率和距离与胡杨林FVC响应格局相关联，相对低至中等频率淹没区的FVC改善强于高频率类别，FVC改善集中在距反复淹没水面0–0.5 km范围内，超过0.5 km后明显减弱。（4）胡杨林FVC与水库泄水量、地下水埋深和年度地表淹没扩展面积的关联在滞后2–3年时最强。年度地表淹没扩展面积、多年淹没频率和距离指标分别从范围、重现与邻近度表征地表水可达性，补充了水库泄水量与地下水记录，为生态输水条件下胡杨林恢复的长期监测与评价提供了空间信息。