关于在热带河流水质指数中考虑土壤侵蚀速率的必要性：来自哥伦比亚塞萨尔河的见解

《Ecohydrology & Hydrobiology》：On the need to account for soil erosion rates in water quality indices for tropical Rivers: Insights from Cesar River, Colombia

　　

大量研究表明，发展中国家的土壤侵蚀正达到临界速率，其场外效应（例如，水体中沉积物浓度增加）可能因全球变暖而加剧。在大多数情况下，水质指数（WQI）无法区分自然负荷和人为负荷，可能导致评估和管理决策的偏差。然而，这种区分在热带河流中尤为重要，因为强降雨、高侵蚀速率和沉积物连通性强烈影响水质。本研究介绍了一种整合沉积物负荷的水质指数（WQIssy），这是一种对沉积物敏感的度量，它将对自然侵蚀的定量输入纳入水质指标。在操作上，WQIssy将RUSLE-GGS（Rosas & Gutierrez, 2020）土壤侵蚀估算与基于主成分分析（PCA）的子指数相结合。WQIssy被应用于塞萨尔河（约20,000 km2），这是下马格达莱纳河（哥伦比亚）的主要支流之一。为此，研究人员分析了来自15个监测站（2023年）的11个水质参数，以及来自DynQual模型的9个虚拟站。研究结果表明，自然侵蚀增加了塞萨尔河下游段的浊度（>100–500 NTU），改变了河流中下游的水质指数分类，其中多个站在WQIssy尺度上从“可接受”变为“一般”或“差”类别。这些发现强调了侵蚀驱动的沉积物动态可以显著改变热带河流水质解释的事实。通过区分自然和人为输入，WQIssy减少了过度诊断，并支持在易发生侵蚀的热带流域中开展成本效益高的河流修复和流域管理。

**论文解读文章**

**研究背景与问题**
土壤侵蚀是自然环境过程，但人类活动（如城市化、农业、采矿）可加剧其强度。在发展中国家，土壤侵蚀已达到临界水平，并可能因全球变暖而进一步恶化，尤其是在热带地区。传统的水质指数（WQI）大多基于专家判断或固定权重，未能区分自然侵蚀与人为污染对水质的贡献，这可能导致评估偏差和管理决策失误。例如，高浊度可能源于自然侵蚀而非人为排放，但常规指数常将其归因于污染，从而引发过度诊断和非必要的修复措施。这一问题在热带河流中尤为突出，因其强降雨、高侵蚀速率和沉积物连通性显著影响水质。现有研究虽已提出多种沉积物质量指数（如污染负荷指数PLI、地质累积指数Igeo等），但这些指标通常与水质指数脱节，且未系统整合自然沉积物负荷。因此，亟需开发一种能区分自然与人为输入的综合性水质指标，以支持热带河流的修复规划与适应性管理。

**研究方法**
本研究基于哥伦比亚塞萨尔河（约22,000 km2）的监测数据，开发了整合沉积物负荷的水质指数（WQIssy）。关键方法包括：
1. **RUSLE-GGS模型**：采用Rosas和Gutierrez（2020）提出的RUSLE-GGS框架，估算流域尺度的土壤侵蚀速率（以t/h/y为单位），并通过沉积物输送比（SDR）函数转换为站点特异性沉积物产量（SSY）。模型校准基于2011–2020年数据，并与实测站点的悬浮沉积物浓度（SSY*）进行交叉验证。
2. **主成分分析（PCA）**：基于2023年15个地面监测站（GS）和9个虚拟站（VS）的11个水质参数（pH、溶解氧DO、总悬浮固体TSS、总氮TN、总磷TP、大肠杆菌FC、总溶解固体TDS、生化需氧量BOD、水温T、化学需氧量COD、电导率EC），以及SSY**（转换为浊度NTU），通过PCA提取主成分，并基于第一主成分（PC1）的载荷计算各参数的经验权重。
3. **WQIssy构建**：将PCA权重与标准化子指数（基于IDEAM方法和NSF曲线）结合，最终通过加权聚合得到WQIssy值。SSY**通过线性回归转换为浊度，并纳入子指数计算。

**研究结果**
**5.1 参数选择与侵蚀模型结果**
RUSLE-GGS模型得出塞萨尔河流域平均侵蚀速率为104.69 t/h/y（2011–2020年），其中57%的流域为自然侵蚀（<1 t/h/y），22%为中度侵蚀（1–10 t/h/y），21%为严重侵蚀（>10 t/h/y）。基于±30%偏差阈值，144个侵蚀速率样本中，样本135-SDR4（SDR4为SDR函数之一）与Canoas大桥站实测值最吻合，被选为代表样本。SSY**沿流域变化范围为29,506至190,000吨，下游段浊度显著增加（>100–500 NTU），尤其在塞萨尔河与Zapatosa沼泽交汇处（GS-15和VS-7）。

**5.2 参数子指数估计**
所有子指数（0–1尺度）显示空间异质性：pH子指数接近1；DO子指数中到高，但部分站（GS-5、GS-8、VS-5–6）下降；BOD和COD子指数分散，GS-3、GS-4、GS-11低值；TSS和TDS子指数变异性大，反映侵蚀和扩散输入差异；FC子指数在多数站低值，指示微生物污染；EC子指数中低值，GS-3、VS-7、VS-8临界；T子指数变异性较低。关键站（GS-4、GS-5、GS-11、VS-4、VS-5）多处参数同时低值。

**5.3 PCA权重推导**
PCA提取四个特征值>1的主成分，PC1和PC2解释54.99%的变异性。PC1中TP、TN、COD、EC载荷较高（0.12–0.14），SSY**载荷为0.09，而FC、TDS、pH权重较低（0.01–0.04）。相比WQI_IDEAM的均匀权重（0.15–0.17），PCA权重更反映数据实际结构，突出了营养盐和有机负荷的主导作用，并纳入了T、BOD、SSY**等额外变量。

**5.4 水质评估**
WQI_PORH（2013）将中游段（GS-4至GS-6、VS-4）评为“很差”，而WQI_IDEAM（2023）评为“可接受”。WQIssy则显示“差”到“一般”，捕捉到自然沉积物负荷的影响。下游段（GS-12至GS-15、VS-6、VS-7），WQIssy识别出“差”状态，而WQI_IDEAM仍为“可接受”，体现了WQIssy对累积污染压力的敏感性。在Zapatosa沼泽，三种指数均显示“一般”状态。

**讨论与结论**
讨论部分指出，WQIssy通过整合自然沉积物负荷，提高了对热带河流水质评估的代表性，避免将自然过程误判为人为污染。研究强调，塞萨尔河的水质不仅受废水排放影响，还受沙矿开采、土地利用和自然侵蚀的驱动。WQIssy的应用有助于区分自然与人为主导的河段，为制定差异化修复目标（如侵蚀控制措施 vs. 污染减排）提供依据。此外，WQIssy具有可复制性和可扩展性，可适应其他热带流域。但研究也指出数据限制（如缺乏详细污染源清单和长期监测数据）以及TSS-浊度关系受颗粒粒径影响的不确定性，建议未来工作纳入沉积物学调查。

**结论部分**：本研究提出了一种整合沉积物负荷的水质指数（WQIssy），旨在为数据稀缺的热带发展中国家提供河流修复评估的技术框架。该指标结合了PCA经验权重和RUSLE-GGS自然沉积物负荷（转换为浊度）。结果表明，WQIssy能够扩展传统指数的解释能力，更真实地反映自然过程对水质的贡献，并强调沉积物相关浊度可显著改变水质分类，因此需区分自然变异与人为压力。WQIssy支持基于地质和水文过程的适应性水资源管理与河流修复规划，有助于平衡生态系统功能、土壤侵蚀输入与人类用水需求。未来研究可整合生态流量等参数，并向区域和国家尺度推广。