中国河流水质指标的多尺度敏感性及其空间异质性驱动机制解析

《Ecological Indicators》：Scale-sensitive river water quality indicators reveal spatially heterogeneous drivers across China

　　

富营养化仍然是中国河流水质的首要制约因素。北方河流的水质明显差于南方河流。景观组成主导了不同地区的水质变化。北方管理应优先考虑农田约束控制。南方应优先考虑城市扩张和社会经济调控。本研究开发了一个多尺度分析框架，整合了机器学习、路径分析和探索性阈值分析，以调查中国河流水质的驱动因素。利用覆盖2021–2024年中国主要流域的全国监测数据集，研究人员考察了异质性区域间河流水质的空间格局。具体而言，本研究旨在：（1）刻画河流水质的区域格局并识别差异化管理的优先污染物；（2）确定异质性区域间的主导驱动因素及其敏感空间尺度；以及（3）阐明多尺度响应路径和候选管理参考范围，以支持区域特异性和尺度敏感性的水质治理。通过明确区域特异性驱动因素、尺度依赖性响应和候选响应转换范围，本研究为中国差异化的河流水质治理提供了基础。研究结果表明，全国河流水质总体保持相对稳定，而南北之间持续存在显著的区域差异。南方水质普遍优于北方，而富营养化仍是全国水质改善的主要制约因素。多模型评估揭示了驱动因素与水质指标之间的非线性关系。驱动因素在河岸缓冲区尺度上对污染相关指标（总氮TN、总磷TP和高锰酸盐指数CODMn）影响最强，而综合水质指数（富营养化指数EI和水质指数WQI）在子流域尺度上得到更好的解释。景观因素是两区域水质的主导驱动因素。社会经济影响在南方更强，而自然影响在北方更为显著。通过整合多尺度驱动分析、路径分析和候选响应转换点，本研究能够更精确地识别水质管理的优先驱动因素和有针对性的干预路径。

论文解读文章

一、研究背景与科学问题

河流作为支撑生态系统稳定和社会经济发展的关键淡水资源，其水质状况直接关系到清洁水的可获得性。然而，气候变化、快速景观演变以及点源与面源污染的叠加效应持续威胁着河流水质安全。全球范围内，超过2000个子流域在水环境方面日益脆弱，其中亚洲面临的风险最为严峻。过去二十年间，中国通过协调污染控制、流域管理和生态修复，显著改善了河流水质，但这些改善呈现出明显的空间异质性：华北部分区域仍面临持续的富营养化风险，西部河流的浊度则较东部有所上升。

既有研究存在若干关键局限。小流域研究通常依赖密集监测网络，在微空间尺度如河岸缓冲带上刻画水质动态，其结论受局地条件强烈影响，难以推广至更广区域；大尺度研究虽能揭示子流域层面的宏观格局，却可能掩盖重要的局地异质性。此外，传统线性方法难以捕捉多驱动因素间的非线性交互作用，且许多水质管理项目仍以单一指标为中心，忽视了污染物耦合及多指标综合控制的必要性。这些局限阻碍了对河流水质尺度依赖性驱动因素的系统理解，也制约了有效的跨区域管理策略的制定。因此，亟需一个整合多源数据、阐明多尺度机制并能够转化为可操作政策的分析框架。

二、研究设计与技术方法

研究人员构建了一套多尺度分析框架，主要采用以下关键技术方法：其一，随机森林（Random Forest, RF）、支持向量机（Support Vector Machine, SVM）和逐步多元线性回归（Stepwise Multiple Linear Regression, SMLR）三种建模方法，用于识别各水质指标最具解释力的空间尺度并排序预测因子重要性；其二，Shapley加性解释（SHAP）分析，用于量化候选驱动因子的相对贡献及效应方向；其三，共现网络分析，用于识别驱动因子与水质指标间的显著关联和潜在耦合关系；其四，偏最小二乘结构方程模型（PLS-SEM），用于检验连接多尺度驱动因子与污染物负荷及综合水质响应的假设路径；其五，非参数变点分析（nCPA），用于识别水质响应沿环境梯度的突变位置，作为管理参考范围。

研究数据来源于中国生态环境部2021年1月至2024年12月的月度监测数据，覆盖中国九大流域的725个监测站点。选取的水质指标包括溶解氧（DO）、电导率（EC）、高锰酸盐指数（COD_Mn）、总磷（TP）、总氮（TN）、氨氮（NH₃?N）以及铜（Cu）、锌（Zn）、汞（Hg）等重金属。以秦岭—淮河线为界划分南北区域，北方包括内陆河流域、松辽流域、黄河流域、海河流域和淮河流域，南方包括长江流域、东南诸河流域、珠江流域和西南诸河流域。空间尺度涵盖95个子流域以及1、3、5、7、10 km五种河岸缓冲区距离。

三、主要研究结果

3.1 中国南北河流水质的时空格局

研究期间河流水质年际变化有限，但夏季出现明显的季节性恶化。空间上，富营养化成为全国范围的主导水质问题，内陆河流域水质相对较好。北方河流水质普遍劣于南方：聚类分析显示，北方流域河流的富营养化指数（EI）显著高出34.68%，水质指数（WQI）低8.59%。贝叶斯层次模型（BHM）分析得出，北方/南方COD_Mn几何浓度比为1.602（95% CrI: 1.486–1.722），TN为1.581（1.455–1.720），TP为1.194（1.061–1.340），三项指标P（北方>南方）均≥0.998。重金属污染指数（NPI）总体低于1.0，表明研究期间重金属污染风险较低。南北水质差异部分归因于水资源禀赋和人口压力的区域对比：北方地表水资源仅为南方的12.04%，却集中了全国40%以上的人口。

3.2 水质对驱动因子的尺度依赖性响应

模型评估显示，随机森林在南北两区域均提供最强的解释力。北方区域TN、TP和COD_Mn的最优解释尺度为7-km缓冲区，而南方区域峰值敏感性出现在5-km缓冲区。这一对比表明，富营养化相关污染物的最敏感空间尺度受区域背景条件调控：北方河流水量较低、稀释输送能力减弱，污染物需经历更长的空间路径才能有效汇入河道；南方密集河网和充沛径流压缩了这一空间足迹。COD_Mn表现出最强的尺度敏感性，北方R2达0.65，南方达0.62。EI和WQI在子流域尺度上解释力最高，表明综合水质与流域尺度环境背景的关联强于局地干扰。

3.3 景观格局作为河流水质的主导驱动因素

随机森林模型识别出景观格局是南北区域河流水质的首要决定因素，其中景观组成贡献最为突出（41.40%–43.40%）。建设用地（CON）和农田（FAR）的扩张与多项污染指标正相关：北方FAR扩张是WQI下降的最强预测因子，南方CON增加与TP和COD_Mn浓度升高相关。非生态用地与生态用地比（NEL/EL）在南北区域均持续重要，与TN、TP、COD_Mn和EI呈显著正相关，反映“强源—弱汇”结构增强营养盐和有机质的径流输出。林地（FOR）和草地（GRA）覆盖增加可通过改善土壤稳定性和增强入渗减少养分流失。水域面积（WAT）在南北区域呈现相反效应：北方较大的水域面积与较低TN相关，可能因厌氧微环境促进反硝化；南方水域扩张与COD_Mn正相关，可能与复合干旱—高温事件下初级生产和藻类生长增强有关。

3.4 自然与社会经济因素对区域水质变化的差异性影响

北方水质对自然背景条件更为敏感，自然因素平均贡献14.05%的解释力。坡度（Slo）与地形湿度指数（TWI）对COD_Mn呈现相反效应，反映地形—水文对有机污染的调控。二氧化氮（NO₂）与TN正相关，表明大气沉降可能是地表水氮输入的额外来源，大气氮沉降可解释地表水TN的24.5%。南方社会经济因素作用更为突出，平均贡献17.80%。人口密度（Pop）是水质变化的关键预测因子之一，但其效应并非严格线性：在经济发展和治理能力较高的区域，人口增长可能对水质产生积极影响，呈现环境库兹涅茨曲线模式，表明Pop在较发达区域可能代表“管理能力”而非单纯的污染压力指标。

3.5 从河岸缓冲区到子流域：尺度特异性管理策略

PLS-SEM结果支持“源—输移—响应”框架：缓冲区尺度驱动因子对TN、TP和COD_Mn等污染物负荷呈现显著直接路径，而EI和WQI所代表的综合水质状态主要与子流域尺度驱动因子关联。南方区域子流域驱动因子对污染物负荷也存在较弱但显著的直接路径。非参数变点分析识别出候选管理参考范围：北方污染相关指标在7-km缓冲区尺度上围绕NO₂=18 μg/m3、NEL/EL=3和Slo=11°出现响应转换；子流域尺度综合水质状态的候选转换点出现在NEL/EL=3和FAR=70%。南方污染相关指标在5-km缓冲区尺度上围绕Tem=22 °C、CON=10%和WAT=9%出现候选响应转换；子流域尺度Pop=207人/km2和SHDI=0.8被识别为综合水质状态的候选参考点。

3.6 综合水质治理的启示

北方水质与农业景观结构和自然背景条件关联更强，南方则更受建设用地扩张和社会经济活动塑造。一刀切的管理方式可能无法应对区域特定挑战，中国河流水质治理应转向区域特异性、尺度敏感性和关键驱动因子导向的方法。

四、讨论与结论

讨论部分指出，缓冲区尺度上气候和空气质量变量的表观重要性并不意味着其效应空间局限于该区域，而是反映了其对河流系统影响最可检测的敏感空间范围。这些驱动因子虽无法直接控制，但可通过间接管理干预加以缓解：北方氮氧化物削减需整合农业管理改善、燃煤和工业排放管控以及排放源空间优化；南方升温效应可通过增加河岸植被覆盖和遮荫来缓解。未来管理应优先关注低坡度农业区域，加强河岸缓冲带和源头控制。

研究结论表明，2021至2024年间中国河流水质总体保持相对稳定，但南北区域差异显著。南方水质普遍优于北方，富营养化仍是全国水质改善的主要制约因素。驱动因子在河岸缓冲区尺度上对污染相关指标影响最强，综合水质指数则在子流域尺度上得到更好解释。景观因素是两区域水质的主导驱动因素，南方社会经济影响更强，北方自然影响更为显著。研究整合多尺度驱动分析、路径分析和候选响应转换点，实现了水质管理优先驱动因素和靶向干预路径的更精确识别，加深了对中国河流水质动态的理解，也为其他空间异质性显著区域的水质治理提供了科学依据。