拉萨河流域溶解态金属的空间分布格局与污染指数：与溶解性有机质的关系

《Toxics》：Spatial Patterns and Pollution Indices of Dissolved Metals in the Lhasa River: Associations with Dissolved Organic Matter

　　

摘要：研究人员于2025年5月对拉萨河32个采样点进行了十种溶解态金属的空间变化及其与溶解性有机质（DOM）关联的评估，采用污染指数、相关性分析、主成分分析（PCA）和溶解性有机碳（DOC）测定等方法。Cd、Co、Cu、Mn、Ni和Zn的局部最大值出现在L8点，其中Mn、Cd和Zn分别达到59.51、0.230和6.58 μg/L，而经多重检验校正后，As、Pb和Sb在河段尺度上的差异得到支持。在所采用的筛选参考值下，所有单因子指数均低于1，而内梅罗指数（Nemerow index）范围为0.0689至0.6980。DOC平均值为1.067 mg/L，与九种金属呈正相关（R2 = 0.132–0.330；Benjamini–Hochberg（BH）校正后q < 0.05），但Cr未显示有支持的相关关系。在敏感性检验中剔除L8后，仅As和Sb保留显著的DOC关联，表明表观的多金属关联部分依赖于热点位点。这些观测结果支持对局部最大值和下游As–Sb富集进行针对性监测，但并未确立DOM介导的迁移机制或量化来源贡献。

拉萨河流域溶解态金属空间分布与污染评估及其与溶解性有机质关联的系统研究解读

一、研究背景与科学问题

河流是陆源物质输送的关键通道，同时维系水生生态系统并支撑人类活动，其水质对自然过程和人为输入的叠加效应高度敏感。重金属因具有环境持久性、不可降解性和生物累积潜力而备受关注，某些元素还可通过食物链传递，对生态系统完整性和人类健康构成长期威胁。然而，金属的环境行为和毒性不仅取决于总浓度，还受化学形态和生物可利用性影响，因此仅基于总浓度的评估可能无法充分反映实际环境风险。溶解性有机质（DOM）作为天然水体中复杂溶解有机化合物的混合物，可通过络合、吸附及介导金属氧化还原转化等方式与金属离子相互作用，从而改变金属的分配和迁移行为，但这种相互作用随DOM组成和分子大小而变化。

青藏高原被誉为"亚洲水塔"和"第三极"，是长江、黄河、澜沧江和雅鲁藏布江等多条亚洲重要河流的发源地，对下游供水和水安全具有关键意义。冰川退缩、降水变化和空气污染物长距离传输可能增加高原河流对环境变化的敏感性。然而，现有高原河流研究多集中于水化学或金属污染的单一方面，DOM空间变化与不同金属分布及污染的系统性关联尚不清楚。拉萨河作为雅鲁藏布江的重要支流，是拉萨地区生活、农业和生态用水的重要水源，其流域兼具青藏高原地质背景与日益增强的人类活动。此前研究主要关注水化学或单一污染物，DOM组成与多种金属之间的空间关联仍不明确。因此，研究人员在2025年5月开展了单次采样调查，旨在描述金属分布与污染指数、考察潜在来源关联，并检验DOC、光学指数和荧光组分与金属的协变关系。

二、研究设计与主要技术方法

研究人员于2025年5月在拉萨河流域32个采样点采集表层水样，其中上、中、下河段分别设10、13和9个站点。溶解态金属浓度采用电感耦合等离子体质谱法（ICP-MS）测定，DOM表征采用溶解性有机碳（DOC）分析、紫外-可见吸收光谱和三维激发-发射矩阵荧光光谱（EEM）结合平行因子分析（PARAFAC）。统计分析包括Kruskal–Wallis检验、Mann–Whitney U检验、Pearson相关分析、主成分分析（PCA）、单因子污染指数和内梅罗综合污染指数评估，并进行了Benjamini–Hochberg（BH）和Holm多重检验校正及敏感性分析。

三、主要研究结果

3.1 拉萨河金属空间分布

金属浓度呈现显著但不均衡的空间变异性。Mn变异系数最大（534%），其次为Cd（193%）和Zn（183%），而Cr变异较小（30%）。Cd、Co、Cu、Mn、Ni和Zn在上游L8点同时达到站点最大值，分别为0.230、0.122、9.58、59.51、1.133和6.58 μg/L，Zn在邻近的L9和L10点也较高，表明存在局部浓度最大值而非连续河段富集。经BH校正后，As、Pb和Sb的河段差异得到支持：下游As中位数为9.933 μg/L，显著高于上中游；Sb下游高于中游；Pb在上游和下游均高于中游。其余七种金属未显示有支持的河段差异。这一结果区分了局部多元素最大值与下游As和Sb的广泛富集。

3.2 拉萨河金属潜在来源

全数据集分析显示Cd、Co、Cu、Mn、Ni和Zn之间存在强正相关，但剔除L8后该模式未一致保留，Cd–Mn仍为正相关（r = 0.746），而Mn–Ni和Cu–Zn分别减弱至r = ?0.023和?0.020。PCA保留三个主成分，解释91.48%方差，PC1对Mn、Ni、Cd、Co、Cu和Zn具有高载荷，但L8占PC1站点得分平方和的91.61%，表明六金属模式强烈受该站点浓度极值影响。PC2与As和Sb关联最强，PC3以Cr为主导。敏感性分析显示标准化无法消除热点影响，因此无法进行定量来源归因。

3.3 重金属污染评估

在所采用的元素特异性参考值下，所有单因子污染指数均低于1，Cu、Sb和As的平均P_i值最高。内梅罗综合污染指数范围为0.0689至0.6980，最大值出现在L8点。尽管所有站点均低于0.7的指数边界，但该筛选分类并非毒性评估或整体I类水质达标的证据，尤其因为Cr(VI)未单独测定。

3.4 DOM特征与空间变化

DOC浓度范围为0.50至1.85 mg/L，平均值1.07 mg/L，表明DOM丰度相对较低。SUVA₂₅₄平均值为3.77 L/mg/m，表明中等至较高芳香性，且从上游至下游总体增加。FI和BIX平均值分别为1.71和1.07，表明陆源和自生源均有贡献。PARAFAC解析出六个荧光组分，其中C1为短波类腐殖质或类富里酸物质，C2为陆源或土壤类腐殖质组分，C3为微生物影响或M峰类腐殖质组分，C4为类色氨酸蛋白组分，C5为可能关联微生物处理或人为输入的类腐殖质组分，C6为类酪氨酸蛋白组分。C1在所有站点占主导，表明拉萨河DOM以低DOC浓度、混合陆源和自生源输入以及类腐殖质荧光组分占优势为特征。

3.5 DOM特征与金属的关系

DOC经BH校正后与As、Cd、Co、Cu、Mn、Ni、Pb、Sb和Zn九种金属正相关（q = 0.0058–0.0457），Sb的R²最大（0.330），Cr无支持关系。显著模型仅解释13.2–33.0%的变异。剔除L8后仅As和Sb保留显著关联，表明多金属正相关部分依赖于热点位点。SUVA₂₅₄与As正相关、与Cr负相关，FI与Cr正相关，C1和C2与Pb负相关，C5与Pb正相关，但这些描述仅反映协变关系，不能证明结合或迁移机制。BIX与Cu、Co、Ni无显著相关。

3.6 研究局限性

由于仅在2025年5月采样一次，研究结果主要反映当时条件，不能推广至其他季节或年份。部分分析受L8点高浓度影响，未来需开展跨季节和不同水文条件的重复采样，并结合水质和金属形态测定。

四、讨论与结论

本研究识别出L8点Cd、Co、Cu、Mn、Ni和Zn的局部最大值，而非沿河的连续富集。所有P_i值均低于1，P_n范围为0.0689至0.6980，但这些指数不能确立生态安全或整体I类达标。六金属协变模式部分依赖L8，其剔除改变了PCA保留主成分数量。DOC平均值为1.067 mg/L，与九种金属正相关，但仅解释13.2–33.0%的变异，剔除L8后仅As和Sb保留显著关联。研究证据支持对局部金属最大值和下游As–Sb富集进行针对性跟踪监测，并需通过重复采样、规范的QA/QC记录以及水化学和形态测定来检验所推断的关联。