Nemerow Pollution Index and Inter-Variable Relationships in Esuk Oro River, Southern Nigeria

技术摘要：尼日利亚南部埃苏克·奥罗河（Esuk Oro River）的内梅罗夫污染指数与变量间关系

问题陈述
本研究针对先前对埃苏克·奥罗河评估中存在的一个关键方法论缺陷进行了探讨，即此前报告的“内梅罗夫污染指数”并非无量纲的浓度与限值的比值。原有的表格数据错误地将几乎所有水质变量都归类为受污染状态。此外，先前的分析未能区分具有上限筛选值（如营养盐）的变量与具有下限目标值（如溶解氧）的变量，从而导致了错误的污染驱动因素结论。作者指出，水质指数中的归一化错误会导致误导性的分类，因此有必要使用站点平均值和适当的基准值对单因子指数及复合内梅罗夫指数进行重新计算。

研究方法
本研究利用了 2022 年 10 月至 2023 年 2 月期间采集自埃苏克·奥罗河三个站点（上游、中游和下游）的水质数据。研究采用了以下分析方法：

• 数据处理： 计算了九项指标的站点平均值：电导率 (EC)、总溶解固体 (TDS)、总悬浮固体 (TSS)、溶解氧 (DO)、生化需氧量 (BOD)、碱度、氯化物、硝酸盐和磷酸盐。
• 指数计算：
  • 单因子指数 ($P_i$)： 对于具有上限值的指标，$P_i$ 计算为测量站点的平均值 ($C_i$) 与筛选值 ($S_i$) 的比值。对于具有下限目标的溶解氧，则使用反比关系 ($5/C_i$)，以确保当浓度低于 5 mg L⁻¹ 时，$P_i > 1$。
  • 复合内梅罗夫指数 ($P_N$)： 使用公式 $P_N = \sqrt{(P_{max}^2 + P_{mean}^2)/2}$ 进行计算，其中 $P_{max}$ 是最大单因子指数，$P_{mean}$ 是九个指数的算术平均值。
• 相关性分析： 直接从五个月度平均值中重新计算了所有变量对的皮尔逊积矩相关系数 ($r$)。显著性水平在 $\alpha = 0.05$（双侧检验）下进行评估。
• 基准值： 应用了操作性筛选值和 WHO 指南（例如：磷酸盐为 5 mg L⁻¹，硝酸盐为 50 mg L⁻¹）。温度、pH 值和酸度未纳入指数计算，因为简单的比值法不适用于基于范围的变量。

主要贡献
本工作的核心贡献在于纠正了污染指数的计算方法。作者成功实现了以下目标：

1. 重新计算指数： 提供了修正后的单因子指数和复合内梅罗夫指数，准确反映了测量值与操作基准值的比率。
2. 确定主导驱动因素： 确定了磷酸盐是复合污染指数的主要贡献者，这一发现此前在原始分析中被掩盖。
3. 澄清溶解氧逻辑： 纠正了对溶解氧的解释，证明高浓度（观测值为 20.50–29.91 mg L⁻¹）在反比法下并不构成污染超标，这与原研究的暗示相反。
4. 重新评估相关性： 提供了一套直接源自月度平均值的新探索性相关性结果，而非直接复制可能存在缺陷的原始数据。

结果

• 污染状态： 复合内梅罗夫指数 ($P_N$) 分别为 3.015（上游）、3.118（中游）和 3.200（下游）。所有数值均超过 1.0，表明存在污染，但其量级是由特定变量驱动的。
• 主导变量： 磷酸盐在所有站点均为主要影响因素，其单因子指数范围在 4.172 至 4.438 之间。下游站点的磷酸盐比率最高（4.438）。
• 次要超标情况： 碱度在所有站点均略微超过其操作基准（1.110–1.162）。上游总悬浮固体 (TSS) 略微超过其基准（1.054）。
• 未超标情况： EC、TDS、BOD、氯化物和硝酸盐的比值均低于 1.0。
• 相关性： 碱度–硝酸盐 ($r=0.987$)、碱度–氯化物 ($r=0.967$) 以及氯化物–硝酸盐 ($r=0.939$) 之间存在显著的正相关关系。相反，TDS 与硝酸盐 ($r=-0.936$)、TDS 与氯化物 ($r=-0.923$) 以及 TDS 与碱度 ($r=-0.919$) 之间存在强烈的负相关关系。

意义与声明
论文谦逊地将研究结果定位为筛选性评估，而非对饮用水安全性或临床健康风险的最终判定。作者强调，该指数旨在识别营养盐富集是埃素克·奥罗河污染信号的主要驱动力。

研究明确指出，较小的样本量（$n=5$ 个月度观测值）限制了建立因果关系或来源归属的能力。相关性分析被视为探索性假设，用于未来的监测，而非证明存在共同来源。作者总结道，修正后的分析强调了以下必要性：

1. 在数据准备过程中保留单位。
2. 使用无量纲的归一化比率。
3. 区分上限值与下限目标值。

最后，论文呼吁对异常高的溶解氧和磷酸盐水平进行原始野外数据的核实，绘制局部污水和农业输入的分布图，并建立完整的 12 个月重复监测计划，以验证这些筛选结果。