地形排水阻抗与山地-高山生态系统中有机质组成及汞分配的关系

《Environmental Pollution》：Topographic drainage impedance is associated with organic matter compositions and mercury partitioning in a montane-alpine ecosystem

　　

本研究考察了地形位置、有机质组成与山地-高山土壤中甲基汞（MeHg）库及其分配之间的关联。研究人员在贡嘎山东坡一条海拔1666–4503米的样带上，于12个采样点各采集了三根60厘米长的土壤柱芯，这些点位覆盖冲积平原、山间盆地和山脊-肩部地形。一个由平面曲率和坡度计算得出的排水阻抗指数（DI）被用于表征相对于下坡排水的潜在侧向水流汇聚情况。研究发现，DI与水提取有机质的加权平均碳名义氧化态（NOSC）呈负相关（校正R2 = 0.31, p = 0.034）。NOSC被用作所检测水提取有机质的基于分子式的分子描述符。在进行探索性的样本层级比较时，21个点位-深度分析单元被划分为低NOSC组（NOSC < -0.2；n = 4）和高NOSC组（NOSC ≥ -0.2；n = 17）。低NOSC组含有较大比例的类脂质分子式贡献和较小比例的类植物化学物质分子式贡献。盆地表层位点的总甲基汞（MeHg）库更大，而山脊位点则具有更大的可操作定义的水提取态甲基汞份额。这些观测结果将DI确定为用于筛选甲基汞库及其分配空间变异的一个候选协变量。未来研究可通过包含直接水文和氧化还原测量的独立样带来检验所提出的过程机制。

论文解读：《地形排水阻抗与山地-高山生态系统中有机质组成及汞分配的关系》

一、研究背景与科学问题

汞（Hg）污染对自然生态系统和人类健康构成严重威胁，其中甲基汞（MeHg）因其在食物网中的强生物富集和生物放大效应而备受关注。尽管水生生态系统的汞循环已被广泛研究，但陆地土壤作为全球最大的汞储库，其如何影响无机汞和甲基汞向下游环境的输送仍不明确。特别是山区森林土壤，既是大气汞沉降的汇，也可能是通过水文输出向下游释放汞的源。然而，在地形陡峭、水文条件、植被和土壤生物地球化学性质在短距离内剧烈变化的复杂景观中，控制汞形态转化、滞留和迁移的机制尚不清楚。现有研究表明，甲基汞与土壤有机质（SOM）的结合较弱，这可能是其比无机汞更具流动性的原因之一，但水提取态甲基汞的比例及其与土壤化学、微生物群落和地形的协变关系在野外研究中鲜有报道。因此，理解地形如何通过影响有机质组成来调控甲基汞的库容和分配，对于评估山地生态系统中汞的环境风险至关重要。

二、研究方法与技术路线

为探究上述问题，研究人员在位于青藏高原东缘的贡嘎山东坡，沿一条约3000米海拔梯度（1666至4503米）设置了12个采样点，涵盖冲积平原（Zone A）、山间盆地（Zone B）和山脊-肩部（Zone C）三种地貌单元。2022年7月，在每个点位采集三根60厘米深度的土壤柱芯，并按0–10、10–20、20–40和40–60厘米分层，共获得144个土壤样品。研究采用的主要关键技术方法包括：1）利用30米分辨率的哥白尼数字高程模型（DEM）提取坡度（β）、平面曲率（C_plan）等地形属性，并构建了一个新的排水阻抗指数（DI = ?C_plan/tan β）来综合表征地形对侧向水流汇聚和下坡排水的阻碍程度；2）测定全土总汞（THg）和甲基汞（MeHg）浓度，并通过温和的水提取操作定义水提取态总汞（dTHg）和水提取态甲基汞（dMeHg）；3）使用傅里叶变换离子回旋共振质谱法（FT-ICR-MS）对水提取有机质（WEOM）进行分子层面的表征，计算碳名义氧化态（NOSC）、芳香性指数（AI_mod）及各类分子式的相对丰度；4）通过氨基酸生物标志物估算微生物坏死物碳（MNC）；5）结合16S rRNA基因扩增子测序进行FAPROTAX功能预测，并对部分样品进行靶向宏基因组学分析以鉴定携带hgcA基因的汞甲基化微生物群落。

三、主要研究结果

3.1 地形特征定义了三个地貌分区

研究人员通过平面曲率和坡度的二维空间分析，将12个采样点清晰划分为三个地貌组：冲积平原（Zone A）位于较平坦或缓坡地带，平面曲率多为正值；山间盆地（Zone B）主要分布在平面曲率为负值的侧向汇聚地形区；山脊-肩部（Zone C）则占据发散且坡度较陡的位置。排水阻抗指数（DI）将这些地形属性整合为一个单一轴，所有四个盆地站点均出现在DI值最高的五个点位中。

3.2 总汞和甲基汞在不同地貌区的分布与水提取分配特征

研究发现，土壤THg和MeHg浓度在中海拔地区最高，并随深度增加而降低。地表层THg和MeHg的平均浓度在盆地（Zone B）最高，分别为109.29 ng g^?1和0.633 ng g^?1，显著高于冲积平原和山脊区域。这表明山间盆地是最大的表层甲基汞库。然而，水提取态甲基汞分数（dMeHg/MeHg）在山脊（Zone C）却显著高于盆地（21.38% vs. 7.32%，p = 0.0238），揭示了拥有最大总甲基汞库的区域并非同时拥有最大的可移动甲基汞比例。

3.3 排水阻抗与WEOM氧化态及有机质组成的关联

研究结果表明，DI与表层土壤水提取有机质的加权平均NOSC呈显著负相关（校正R2 = 0.313, p = 0.034），即排水越受阻（DI越高），WEOM的碳氧化态越低。此外，NOSC与微生物坏死物碳占SOC的比例（MNC/SOC）呈正相关（校正R2 = 0.500, p = 0.006），而与含硫分子式（CHOS）的比例呈负相关（校正R2 = 0.321, p = 0.032）。这些发现将地形指标与有机质的分子组成和热力学特性联系了起来。

3.4 排水阻抗与有机质组成与汞库及分配的关联

DI与log₁₀转换的表层MeHg浓度呈显著正相关（校正R2 = 0.272, p = 0.047），表明高排水阻抗的地形更有利于甲基汞的积累。虽然水提取态MeHg分数与MNC/SOC的正向关联未达到统计显著性，但整个观察序列（DI→NOSC→MNC/CHOS→MeHg）构成了一个连贯的协变模式。

3.5 不同NOSC操作分组的有机质特征差异

根据NOSC阈值（?0.2）划分的低NOSC组（n=4，均为盆地B1-B4表层样）相比高NOSC组，具有显著更低的AI_mod（0.313 vs. 0.359）、更多的类脂质分子式（32.06% vs. 14.95%）和更少的类植物化学物质分子式（62.67% vs. 76.95%）。MNC/SOC和基于强度的CHOS相对丰度虽名义上有差异，但未通过多重检验校正。

3.6 对未来筛选与实地测试的启示

DI作为一个易于从数字高程模型中推导的参数，在本样带中与WEOM组成和MeHg库表现出关联。它可作为独立样带研究中的一个候选协变量进行评价。同时，盆地与降水峰值区域的重合暗示了气候与地形的协同作用，需要通过更窄气候带的重复样带来分离这些影响因素。

四、讨论与结论

本研究在贡嘎山样带上的观测结果形成了一条清晰的因果链假设：地形排水阻抗（DI）通过影响土壤水分状况和氧化还原条件，驱动了水提取有机质（WEOM）的分子组成趋于更低的碳氧化态（NOSC）和更高的脂质含量，这种有机质组成的变化进而影响了甲基汞（MeHg）的生成、固定与分配。具体而言，排水不畅的盆地环境有利于形成较大的总MeHg库，但其水提取态比例较低，意味着MeHg被更牢固地结合在有机质上；相反，排水良好的山脊环境中总MeHg库较小，但水提取态比例更高，显示出更强的潜在流动性。

研究结论：在贡嘎山这条样带上，更大的排水阻抗与更低NOSC的WEOM、不同的有机质组成、更大的总MeHg库以及盆地站点较小的水提取态MeHg分数相关联。因此，排水阻抗指数（DI）是一个可用于筛选甲基汞库及其分配空间变异的候选协变量。其普适性应在未来通过包含直接水文和氧化还原测量的独立样带研究中进行验证。该论文发表于《Environmental Pollution》。