长期生态输水驱动下塔里木河下游生态水文变化：地下水、相对地表盐度与植被生产力的关联

《Water》：Ecohydrological Changes Associated with Long-Term Ecological Water Conveyance in the Lower Tarim River: Linking Groundwater, Relative Surface Salinity, and Vegetation Productivity

　　

摘要：为重新连接极度干旱的内陆河廊道，实施了长期生态输水，但其生态效益不能仅凭输水量推断。本研究利用2000—2023年的年输水记录、地下水观测、生长季基于Carnegie–Ames–Stanford Approach（CASA）模型的净初级生产力（NPP）、Landsat衍生的归一化差异盐度指数（NDSI）、气候数据以及四种操作性NPP-NDSI状态，评估了塔里木河下游的生态水文变化。年均输水量为4.08×108 m3，2017年达到峰值12.09×108 m3。平均地下水埋深从2010年的7.49 m降至2023年的4.25 m。包含当前年份的多时间窗口输水量与地下水埋深的负相关关系强于同年输水量（W3: rho = ?0.743; W4: rho = ?0.782）；Newey–West检验在滞后1–3期仍显著，移动块自助法区间对W3和W4排除了零值，而严格的前三年指标则较弱。高NDSI面积从1795.04 km2减少至42.19 km2，而较高生产力面积从91.81 km2增加至1127.81 km2。配对增强型专题制图仪加/陆地成像仪（ETM+/OLI）分析显示高NDSI对应性（r = 0.980），但OLI减ETM+的平均偏差为?0.0122；相应的2010年跨传感器调整量占绝对原始2010–2015年平均NDSI变化的63.9%。因此，跨时期NDSI类别面积和采样状态轨迹被视为对传感器敏感的定量描述。空间支撑敏感性也很重要：功能受限的低盐度状态（FLLS）是精细对齐网格上2023年最大的状态（53.93%），而较高生产力低/中NDSI状态（HPLN）在聚合至500 m后成为主导（66.57%）；因此，不做出关于2023年主导状态的尺度不变性论断。固定2000年土地覆盖的CASA运行未产生实质性的NPP变化，因为共同掩膜中表示的土地覆盖转变发生在共享相同NDVI端元和εmax的类别之间。总体而言，研究结果支持异步指标变化和基于状态的监测框架，但时间依赖性、不规则的地下水采样、残余空间自相关、跨传感器不连续性和空间支撑限制了因果和精细尺度的推断。

一、研究背景与科学问题

在全球干旱区，由于水资源过度开发和气候变化，许多内陆河流廊道发生断流，导致地下水位下降、土壤盐渍化加剧和植被退化。生态输水工程被广泛用于恢复此类受损生态系统，但其长期生态效应难以单纯依据输水量来判断。塔里木河下游作为典型的极度干旱内陆河，自2000年起实施长期生态输水，然而，地下水、地表盐度和植被生产力之间的协同演变关系尚不明确，尤其是缺乏将多指标耦合进行系统性评估的研究框架。现有研究多关注单一指标的响应，未能揭示各生态水文组分间的异步变化特征及其对监测尺度的敏感性。因此，亟需建立一套综合评估方法，以诊断生态输水的实际效果，并为适应性管理提供科学依据。

二、研究内容与主要结论

研究人员基于2000年至2023年长达24年的观测数据，系统分析了塔里木河下游在长期生态输水驱动下的生态水文演变过程。主要结论包括：年均输水量为4.08×10⁸ m³，地下水埋深从2010年的7.49 m显著上升至2023年的4.25 m；包含当前年份的多时间窗口累计输水量（如W3和W4）与地下水埋深的负相关性远强于同年输水量，表明地下水对输水具有累积性响应而非即时响应。地表相对盐度指标（NDSI）高值区面积从1795.04 km²锐减至42.19 km²，而植被生产力（NPP）高值区面积从91.81 km²增至1127.81 km²。然而，这两类指标的变化并非完全同步，表明盐度降低并不必然伴随生产力提升。基于NPP和NDSI交叉分类建立的四种生态水文状态（HPLN、PSE、WSS、FLLS）显示，2023年主导状态取决于分析尺度：在精细对齐网格上，功能受限的低盐度状态（FLLS）占比最大（53.93%）；而在500米本征分辨率下，较高生产力低/中NDSI状态（HPLN）成为主导（66.57%）。这一发现强调了空间支撑对状态判定的关键影响。此外，固定2000年土地覆盖参数的敏感性实验表明，土地覆盖变化本身并未通过改变模型参数路径来驱动NPP趋势，因为涉及的类别间共享相同的NDVI端元和最大光能利用率（ε_max）。

三、关键技术方法

研究人员采用了多种遥感与统计分析技术。首先，利用Carnegie–Ames–Stanford Approach（CASA）光能利用率模型，基于MODIS NDVI和ERA5-Land气象再分析数据估算生长季净初级生产力（NPP）。其次，基于Landsat Collection 2地表反射率计算归一化差异盐度指数（NDSI），作为相对地表盐度的光谱代理。为处理跨传感器不一致性，研究人员对Landsat 7 ETM+与Landsat 8 OLI进行了近日期配对的交叉定标分析，量化了传感器偏移对NDSI时序的影响。在统计推断方面，采用Spearman秩相关、Newey–West异方差自相关一致（HAC）估计和移动块自助法（Moving-block bootstrap）来处理时间序列中的依赖性问题，并对地下水监测网络的不平衡采样进行了详细的月度覆盖清单分析和严格公共井面板敏感性检验。最后，通过将分析结果聚合至500米空间分辨率，评估了空间支撑对生态水文状态判定的影响。

四、研究结果

4.1 长期生态输水格局与水文气候背景

2000至2023年间，年输水量波动剧烈，平均为4.08×10⁸ m³，2017年达峰值12.09×10⁸ m³。六个遥感分析年份的平均降水量极低（15.79–32.37 mm），维持了极度干旱的气候背景。

4.2 地下水变化及其与输水的关联

地下水埋深总体变浅，均值从2010年的7.49 m降至2023年的4.25 m。同年输水量与地下水埋深的弱相关（rho = ?0.2786）表明，地下水对输水的响应具有滞后性和累积性。包含当前年份的三年滚动平均（W3）和四年滚动平均（W4）均显示出强负相关（rho分别为?0.7429和?0.7821），且通过了Newey–West HAC检验和移动块自助法检验。然而，严格的前三年平均输水量（AnteW3）的关联性较弱，其自助法置信区间跨越零值，说明精确的滞后结构尚未被唯一识别。

4.3 生产力与相对地表盐度指标的异步变化

高NDSI面积从2000年的1795.04 km²骤降至2023年的42.19 km²，而较高生产力面积从91.81 km²增至1127.81 km²。尽管两者均呈现改善趋势，但盐度降低的面积远大于生产力提升的面积，表明指标间存在显著的异步性。跨传感器分析显示，OLI相比ETM+的NDSI平均偏差为?0.0122，该偏差占2010–2015年原始NDSI变化的63.9%，因此跨时期NDSI面积被定性为对传感器敏感的估计。

4.4 生态水文状态重组

在精细对齐网格上，水盐胁迫状态（WSS）面积占比从2000年的70.69%降至2023年的1.65%，而HPLN从3.61%升至44.40%。FLLS在2023年占53.93%，成为最大状态。但在500米分辨率下，HPLN占比跃升至66.57%，FLLS降至约33%。因此，研究人员不做关于2023年主导状态的尺度不变性论断。

4.5 多时期采样状态序列

六期采样状态序列分析揭示了多种演变路径。最大的一类是WSS→FLLS且未观察到HPLN的序列，占46.66%，表明大量像元经历了从高盐低产到低盐低产的转变，但未跨越生产力阈值。FLLS→HPLN序列占21.71%，而持续HPLN仅占2.91%。这些序列对起始年份和分析窗口高度敏感。

4.6 2023年低盐度状态间的地下水差异

HPLN区域内的监测井平均地下水埋深为3.94 m，显著浅于FLLS区域的5.25 m（Mann–Whitney U检验，p = 0.002，效应量 = 0.702）。协变量调整后的普通最小二乘（OLS）回归显示，HPLN相对于FLLS的地下水埋深偏浅1.32 m。然而，残差Moran‘s I显著，表明空间自相关未被完全消除，因此该对比被保留为基于站点而非空间独立因果推断的证据。

4.7 生态水文状态与水管理启示

状态框架提供了差异化监测的诊断基础。HPLN区域应优先维护水文连通性并持续监测；FLLS区域需在追加输水前进行根区水分和植被状况的实地核查；WSS区域则应验证土壤电导率和地下水条件。由于状态判定对尺度敏感，建议以500米本征分辨率作为区域HPLN/FLLS验证的主要基础。

五、讨论与结论

讨论部分指出，地下水与输水的关联支持累积性时间对应关系，但无法唯一确定延迟补给的物理机制。NDSI与NPP的空间不匹配应被解释为指标行为的异步性，而非特定水盐机制的分离证据。状态序列对时间窗口和空间支撑的敏感性表明，离散的快照观测不能等同于连续的生态历史。研究强调，环境流量决策需考虑流量大小、时机、持续时间和频率等多维属性，而本研究的贡献在于提供一个诊断性框架，用以优先安排实地验证需求，而非制定最优输水量方案。

研究结论：2000至2023年间，塔里木河下游地下水变浅、高NDSI面积减少、较高生产力面积增加，但这些指标变化不同步。依赖感知的分析支持地下水埋深与包含当前年份的多时间窗口输水量（W3/W4）之间的负关联，而严格的前期输水证据较弱。固定2000年土地覆盖的CASA运行表明，土地覆盖变化驱动的参数路径并未产生NPP趋势。跨传感器比较显示，传感器调整量占原始2010–2015年NDSI变化的63.9%；500米重复分析显示，2023年主导状态从精细网格上的FLLS转变为本征NPP支撑下的HPLN。这些敏感性分析量化了主要的时间、传感器和尺度不确定性。高NDSI面积从1795.04 km²降至42.19 km²的描述性观察被保留为对传感器敏感的结果，而非环境变化的传感器不变性量级。