京津冀地区过去四十年（1984–2024）水体透明度演变：水文扩张与人为干预驱动

《Ecological Indicators》：Evolution of water clarity across the Beijing-Tianjin-Hebei region over the past four decades (1984–2024): Driven by hydrological expansion and anthropogenic intervention

　　

摘要翻译：鉴于水资源短缺加剧和快速城市化，刻画京津冀（BTH）地区地表水环境的长期演变对区域协调发展至关重要。本研究利用1984–2024年Landsat影像，量化了京津冀地区1105个水体的Secchi盘深度（Z_sd）动态，并通过相关性分析和广义线性模型（GLMs）识别了主要驱动因素。结果表明，京津冀地区水体以小型为主（1.33 ± 10.32 km2），基线水透明度普遍较低（平均Z_sd = 0.62 ± 0.22 m）。然而，检测到普遍长期变清趋势：过去四十年中60.0%的水体（N = 663）透明度显著增加，而仅4.2%（N = 46）显著下降（P < 0.05）。GLM分析揭示，水域面积扩张是Z_sd变异的主要驱动因素，持续解释超过45%的长期变化。这种扩张与河北省降水变率密切相关，而在北京和天津则主要由南水北调工程驱动。土地覆盖转变作为次要驱动因素，解释了27.2%的观测方差，其中裸地减少和不透水面增加，加之强化环境治理，有效缓解了沉积物负荷。相反，某些地区（尤其是河北）的集约化农田可能仍产生局部影响，农业径流可能部分限制透明度改善。总体而言，这些发现表明战略性调水和生态恢复与区域地表水质量改善相关，而局部农业压力需要更有针对性的管理干预。本研究为理解自然和人为驱动因素如何交互塑造这一高度缺水地区长期水透明度动态提供了新见解，并为知情的水资源和土地覆盖管理及区域可持续发展提供了坚实的科学基础。

论文解读文章

**研究背景与意义**

京津冀（BTH）地区是中国典型的水资源短缺区域，地表水有限，以小型浅水水体为主。这些水体对维持区域生态平衡和社会经济可持续发展至关重要，尤其考虑到北京和天津作为国家政治经济中心，河北作为主要农业生产基地。然而，由于水体规模和深度有限，京津冀地区的水生生态系统本质上较为脆弱，易受人类活动和气候变率影响。近年来，快速城市化、工业化、人口增长和集约农业加剧了营养盐和沉积物输入，导致水体透明度下降和富营养化问题。大型水文干预措施，特别是2014年通水的南水北调工程（SNWDP），显著改变了区域水资源可得性和水文连通性。尽管已有研究调查了少数水体的水质变化，但缺乏对整个京津冀水体网络长期趋势及其驱动因素的系统性区域评估。因此，开展本研究对于提升科学认知、支持可持续水资源管理和区域协调发展至关重要。

**研究内容与结论**

研究人员利用1984–2024年Landsat卫星影像（TM、ETM+、OLI、OLI2），基于QAARGB半解析算法反演了1105个水体的Secchi盘深度（Z_sd），并结合Spearman相关分析和广义线性模型（GLMs）量化了自然和人为驱动因素的相对贡献。研究得出以下主要结论：（1）京津冀地区水体以小型为主（平均面积1.33 km2），基线透明度普遍较低（平均Z_sd = 0.62 ± 0.22 m）；（2）60.0%的水体（N = 663）呈现显著的长期变清趋势，仅4.2%（N = 46）显著变浊（P < 0.05）；（3）水域面积扩张是Z_sd变异的主要驱动因素，解释超过45%的长期变化，其中河北主要受降水变率影响，北京和天津主要受南水北调工程驱动；（4）土地覆盖转变（裸地减少、不透水面增加）作为次要驱动因素，解释约27.2%的方差；（5）集约农田（尤其在河北）产生局部负面影响，农业径流部分限制透明度改善。该研究发表在《Ecological Indicators》，为区域水资源和土地覆盖管理提供了坚实的科学依据。

**主要技术方法**

本研究采用以下关键技术方法：（1）基于Landsat系列传感器（TM、ETM+、OLI、OLI2）的QAARGB半解析算法反演Z_sd，该算法利用红、绿、蓝波段估算固有光学性质（IOPs）和漫衰减系数（K_d）；（2）跨传感器一致性校正，通过同一天过境的影像对（ETM+-OLI、TM-ETM+、OLI-OLI2）建立比值因子，将TM、ETM+和OLI2的Z_sd调整至OLI参考基准；（3）使用Google Earth Engine（GEE）平台处理12515幅影像，采用月最大合成和3×3像元窗口避免邻近效应；（4）统计驱动分析，包括Spearman相关分析、广义线性模型（GLMs）和方差膨胀因子（VIF）共线性诊断。样本来源：2017–2025年期间在6个水体（白洋淀、官厅、金海、潘家口、沙河、王快）采集145个原位Z_sd测量值，其中50对用于验证（±7天时间窗口）。

**研究结果**

**模型验证与跨传感器一致性**：QAARGB算法在OLI传感器上表现良好（R2 = 0.85，RMSE = 0.46 m，MAPE = 30.1%）。跨传感器校正显示TM和ETM+的Z_sd系统偏低，通过比值因子（Z_sd^TM = 0.629 Z_sd^OLI，Z_sd^ETM+ = 0.510 Z_sd^OLI，Z_sd^OLI2 = 0.980 Z_sd^OLI）实现一致化。

**时空格局**：1984–2024年间，平均Z_sd范围为0.17–2.20 m，总体均值0.62 ± 0.22 m。北京（0.74 ± 0.22 m）> 河北（0.62 ± 0.23 m）> 天津（0.54 ± 0.13 m）。26.2%水体Z_sd < 0.50 m，73.8%高于0.5 m。城市水体（0.64 ± 0.23 m）比农村水体（0.59 ± 0.20 m）稍清澈。趋势分析显示60.0%水体显著变清，4.2%显著变浊，变清趋势广泛分布，变浊集中西北部。

**驱动因素分析**：Spearman相关分析显示水域面积（WA）、不透水面比例（IF）和GDP普遍与Z_sd正相关，而温度、降水和农田比例（CF）多负相关。GLM分析表明环境因素整体解释超过60%的年际变异，其中水域面积贡献最大（北京43.9%、天津53.6%、河北44.2%），裸地比例（BF）次之（北京28.6%、天津16.9%、河北21.7%），不透水面和农田贡献较小。GDP在GLM中贡献微弱（0.3–1.1%），表明其影响主要通过土地覆盖变化间接作用。

**气候与调水工程效应**：河北水域面积与降水显著相关（r = 0.55，p < 0.01），而北京和天津无显著相关，表明非气候驱动因素主导。南水北调工程自2014年后对区域水资源贡献增加（2023年约占14%），与水域面积扩张时间吻合。调水后降水与水域面积的相关性增强，显示调水通过补给地下水增强了水文响应。

**土地覆盖转变效应**：裸地减少在51.4%水体中促进透明度改善，不透水面增加与污水处理设施改进和环境管理强化相关。农田比例在河北贡献约2.7%的Z_sd变异，与农业径流导致的营养盐输入和藻类增殖有关。

**讨论与结论**

讨论部分强调了水域面积作为主要驱动因素揭示了维持生态水位的重要性，例如白洋淀的生态水位范围约7.8–8.6 m。南水北调工程已超越单纯供水功能，成为区域水质改善的重要驱动机制。土地覆盖变化方面，城市化伴随的环境治理（如污水处理）促进了透明度提升，但农田压力在河北仍需针对性管理。研究局限性包括早期Landsat数据缺失、验证样本有限、未显式量化调水工程对单个水体的影响、以及以面积代替体积的近似。未来研究应整合多参数水质监测和机器学习方法。

结论部分重申：京津冀地区水体透明度总体改善，水域面积扩张是主要驱动因素（解释>45%），河北以降水驱动为主，北京和天津以南水北调等人类干预为主；土地覆盖转变（裸地减少、不透水面增加）解释约27.2%方差；农田影响局部显著（约2.7%）；GDP作为宏观指标间接通过土地覆盖变化起作用。这些发现为区域水资源管理和可持续发展提供了科学依据。