青海-西藏高原湖泊沉积物中金属(类金属)的十年尺度对比

《Environmental Science and Ecotechnology》：Decade-scale contrasts in sediment metal(loid)s across Qinghai–Xizang Plateau lakes

　　

青海-西藏高原的湖泊沉积物捕捉了高山生态系统中变化的金属(类金属)动态，该生态系统独特地易受气候变化和人为压力的影响。然而，包含定量源解析的全高原生态风险评估仍然稀缺，留下了阻碍目标生态系统保护和精确政策干预的关键空白。研究人员在此提出一项跨110个湖泊在两个十年期（2014年和2024年）的沉积物地球化学高原尺度分析，利用混合正定矩阵因子分解-随机森林（PMF-RF）框架。砷（As）作为主要诊断哨兵，在两个十年中持续超过基线阈值。尽管平均铬（Cr）浓度从85.7 mg kg^-1降至68.0 mg kg^-1，但更广泛的多金属库表现出显著的区域稳定性，低基线生态风险以镉（Cd）毒性为主，同时在托素湖、尕海湖和西里沟湖存在局部高风险异常。混合源解析表明，交通排放、跨境大气输入和岩性风化的混合足迹日益驱动地球化学变异，在十年内从44%激增至66%，而独立的砷源保持约15%不变。这些模式揭示，虽然自然风化约束区域背景，但局部人为信号在新兴基础设施走廊周围迅速增强。该预测框架建立了一个可转移的模板，用于追踪脆弱高山源头中的污染物动态，为全球针对性的、源特异性缓解提供了重要基线。

### 研究背景、问题与目标

青藏高原作为“亚洲水塔”，其湖泊沉积物是记录气候变化和人类活动影响的敏感指标。然而，现有研究存在三个关键空白：第一，多数研究局限于局部案例，无法揭示全高原尺度的空间异质性；第二，数据来源和方法不一致，阻碍了跨时期的可比性；第三，依赖单一源解析方法，难以区分人为输入、自然风化和气候驱动的再释放。为填补这些空白，研究人员开展了全高原范围的沉积物重金属（类金属）生态风险评估与定量源解析研究。

### 研究内容与结论

研究人员汇编了青藏高原110个湖泊在2014年和2024年两个十年期的沉积物重金属（类金属）浓度数据（Cr、Ni、Cu、Zn、As、Cd、Pb），共152个观测记录。采用混合正定矩阵因子分解-随机森林（PMF-RF）框架进行源解析，并结合地积累指数（Igeo）和潜在生态风险指数（PERI）评估污染与生态风险。结果表明：砷（As）作为主要诊断哨兵，两个时期均持续超过基线阈值；多金属整体呈区域稳定性，但镉（Cd）毒性主导低基线生态风险，在托素湖、尕海湖和西里沟湖存在局部高风险异常。混合源（交通排放、跨境大气输入和岩性风化）的贡献从44%增至66%，而独立砷源稳定在约15%。模式揭示自然风化约束区域背景，但人为信号在基础设施走廊附近迅速增强。该研究建立的预测框架为追踪脆弱高山源头污染物动态提供了可转移模板，为全球针对性、源特异性缓解策略奠定基线。论文发表在《Environmental Science and Ecotechnology》。

### 主要关键技术与方法

数据来源于《青藏高原湖泊科学考察报告》及文献检索，涵盖110个湖泊，按2014年和2024年两个时期进行时间统一化处理（表面沉积物0–20 cm层）。关键技术方法包括：（1）地积累指数（Igeo）和潜在生态风险指数（PERI）结合蒙特卡洛模拟评估污染与风险；（2）Mantel检验和相关性分析筛选环境因子与重金属（类金属）的关联；（3）正定矩阵因子分解（PMF）解析四个潜在源因子；（4）随机森林（RF）回归模型以PMF因子得分为预测变量，对源-元素关系进行非线性细化与交叉验证评估。

### 研究结果

#### 3.1 重金属（类金属）的空间分布

通过对比两个时期的浓度统计和空间分布图，研究人员发现砷（As）呈异质性增强，上尾扩张和变异系数增大表明热点区域驱动区域变异；铬（Cr）浓度均值从85.7 mg kg^-1降至68.0 mg kg^-1，趋向均质化；镉（Cd）、铜（Cu）、镍（Ni）、铅（Pb）和锌（Zn）整体稳定，但局部有变化。空间上，砷在西藏北部和柴达木盆地富集扩张；铬和镍在岩性控制下形成北部富集带；铅和锌在那曲等地增加，可能与青藏铁路交通排放有关。

#### 3.2 生态风险及其空间分布

地积累指数（Igeo）显示大多数金属（类金属）处于无污染至轻度污染级别，无显著全样本时期变化。潜在生态风险指数（PERI）表明整体风险低，镉（Cd）是主要贡献者（因其毒性响应系数高）。高风险点主要位于北部湖泊（托素湖、尕海湖、西里沟湖），2014年西部湖泊（如Beroze、Jiduocuo）风险较高，但2024年西部风险降低，可能与区域湿润化和湖泊扩张导致的稀释效应有关。北部高风险的持续存在与当地未受管控的矿区和交通排放相关。

#### 3.3 源解析

##### 3.3.1 相关性分析

2014年，铬（Cr）、铜（Cu）、镍（Ni）、铅（Pb）强正相关，砷（As）与其他元素无显著相关，社会经济特征（SECs）与砷显著相关，生物群落特征（BCs）与所有元素正相关。2024年，镉（Cd）与多元素相关增强，沉积物特征（SCs）与多数元素显著相关，表明人为输入和沉积物性质调节作用增强。

##### 3.3.2 PMF和RF模型

PMF解析出四个源因子：（1）混合源（2014-FAC1/2024-FAC2）：包含交通排放、跨境大气输入和岩性风化，贡献从44%增至66%；（2）砷主导源（2014-FAC2/2024-FAC1）：煤炭燃烧和冶炼，贡献率稳定在~15%；（3）铅主导源（2014-FAC3/2024-FAC4）：长距离大气传输；（4）过渡源（2014-FAC4/2024-FAC3）：从地质源铬转变为交通/采矿相关铜。PMF-RF混合框架通过交叉验证（R² 0.927–0.985）证实了源因子的稳健性和非线性细化能力。

##### 3.3.3 模型评估

四因子解在统计拟合、稳定性和环境可解释性上最优，残差分布对称，旋转参数（Fpeak）扫描显示解稳定。RF保留了PMF解析的源特征，并加强了主导元素贡献（如铬在2014年第四因子从76.8%增至80.7%）。PMF和RF互补，PMF更适于初级源分配，RF细化元素特异性贡献。

### 总结讨论与结论

讨论部分指出研究的局限性：依赖于两个统一化时期而非连续监测，数据异质性带来不确定性，生态评估仅基于总浓度而缺乏形态数据。源解析结果为模型支持的推断，而非直接观测证据，需结合同位素、大气沉降和过程模型进一步验证。

研究结论翻译如下：结果表明，两个统一化时期之间的差异主要表现为砷（As）的热点行为和空间再分布，而非多金属均匀富集。地积累指数（Igeo）指示总体污染水平低且时间稳定，但局部生态风险热点持续存在；镉（Cd）主导潜在生态风险指数（PERI），尽管砷在基准超标方面更突出。源解析解决了四个源类型，时期间类型恒定但相对强度变化，混合人为-自然影响在关键枢纽周围增强。PMF-RF框架在全高原尺度上提供互补且内部一致的源诊断。总体上，该高山湖泊系统中可检测的人类世信号主要表现为在广泛自然源主导背景下的人为局部叠加。确定的优先元素、热点区域和源类型为风险导向监测和源敏感性管理提供了实用基础。更广泛而言，这些发现对亚洲水塔的环境安全维度具有重要意义，因为高原湖泊作为气候脆弱源头环境中的污染物积累与再释放敏感哨兵。尽管本沉积物结果不应视为下游水质广泛恶化的直接证据，但它们表明高原水资源保护应更重视湖泊中持续或扩展风险信号的早期预警、热点导向监测和源针对性干预。特别是，热点及其明显的北扩需要优先监测，镉应作为核心风险指标（因其主导PERI），北部湖泊和德令哈等枢纽区域需持续管理关注。保护高原湖泊质量也关乎维护亚洲水塔的长期环境韧性。