淇澳岛红树林湿地地下水-沉积物系统中微生物群落的溯源解析：多样性缓冲与垂直分异

《Water》：Enhancing Energy Dissipation in Piano Key Weirs Using Curved Inlet and Outlet Keys

　　

红树林湿地是陆地和海洋环境之间的关键生态过渡带，然而其地下水-沉积物系统中微生物群落的空间分布和环境来源仍缺乏充分表征。研究人员利用基于SourceTracker2贝叶斯算法的微生物源追踪技术，对中国珠江口淇澳岛一条红树林湿地样带的22个样品进行了分析——这些样品包括一口淡水井、潮上带和过渡带、向海的红树林湿地测压管，以及来自红树林植被区和相邻无植被区的两个60厘米沉积柱。源谱的香农多样性在整个梯度范围内被严格约束（均值3.15±0.26；变异系数8.4%），且在样品组间无显著差异（Kruskal-Wallis检验H=7.38，p=0.194），表明存在强烈的多样性缓冲效应。相比之下，源组成呈现高度结构化特征（PERMANOVA伪F=4.14，R2=0.56，p=0.0001）。红树林柱状沉积物以沉积物相关源（9.98%）、根际源（8.58%）和土壤源（6.69%）为主，而邻近的无红树林柱则以淡水沉积物源（17.52%）和海相沉积物源（9.29%）为主。将104个参考源聚合成八个生态类别后发现，植被和土壤衍生源占红树林柱的29.5%，但在无红树林柱中仅占20.1%（p=0.035），而淡水来源则呈现相反模式（15.3%对32.2%，p=0.002）。废水和宿主相关源在红树林湿地测压管中合计贡献了39.3%，揭示了明显的人为影响印记。深度剖面分析显示，仅在红树林柱中存在系统性垂向结构，海相和未分化沉积物源随深度递减（β=?0.377和?0.115%·cm?1；R2=0.71和0.76；p=0.034和0.025），而植被和土壤衍生源随深度递增（β=+0.375%·cm?1，R2=0.75，p=0.026）；无红树林柱在垂向上呈均质状态。这些结果量化了红树林地下水系统中陆地、海洋和人为微生物输入的相互作用，并为监测沿海含水层中的污染物迁移和生态系统完整性提供了基线。由于该设计仅包含单次湿季采样的单一断面，且深度区间存在伪重复，因此这些发现应被视为空间明确的基线，而非对红树林微生物组的普适性描述。

一、研究背景与科学问题

红树林湿地作为地球上生产力最高、生态意义最重大的海岸生态系统之一，在碳封存、海岸保护和养分循环等关键服务功能中扮演着不可替代的角色。然而，尽管红树林沉积物的微生物生态学已受到广泛关注，关于这些微生物群落的环境来源仍存在重大认知空白：沉积物微生物组主要通过原位筛选海水传播的繁殖体、从陆地土壤源招募、还是通过红树林根系和凋落物的植物接种进行组装？人为微生物输入在塑造群落结构中发挥了多大作用？微生物源贡献如何沿环境梯度和沉积物剖面发生空间变化？这些问题之所以悬而未决，是因为传统的扩增子测序虽能有力刻画群落组成，却难以区分可能含有系统发育相似类群的不同源环境。为此，研究人员首次将定量微生物源追踪技术应用于珠江口淇澳岛红树林湿地的地下水-沉积物连续体，旨在填补这一知识空白。

二、主要关键技术方法

研究人员采用SourceTracker2贝叶斯算法进行微生物源追踪，该算法基于吉布斯采样构建混合模型，估算候选源环境对各汇样品的比例贡献。源环境来自涵盖104种环境类型的精选参考文库，包括陆地源（土壤、根际、植物相关）、水生源（淡水、海水、沉积物）、湿地特异性源、人为源（废水、活性污泥、人和动物肠道）、地下环境及专门栖息地。此外，研究人员将这104个参考源进一步聚合为八个互斥的生态类别：植被与土壤、淡水、海洋与河口、地下水与地下、废水与工程、宿主相关与粪便、食品与发酵、以及未分化沉积物。所有22个汇样品均来自淇澳岛一条完整水文梯度样带，涵盖淡水井、潮上带、过渡带、红树林湿地四个水文区，以及一个红树林植被区和一个相邻无植被潮滩的配对沉积柱（各分六层）。

三、研究结果

3.1 微生物源多样性沿环境梯度保持缓冲

微生物源谱的香农多样性在整个样带上高度受限，均值为3.15±0.26，变异系数仅为8.4%。Kruskal-Wallis检验未检测到样品组间的显著差异（H=7.383，p=0.194），最大对比（红树林柱与无红树林柱）也不显著（Mann-Whitney U=30，p=0.065）。这表明α多样性本身在该系统中是一个不敏感的监测指标，强烈的组成差异代表的是贡献源的种类更替而非丰富度的增减。

3.2 不同湿地带的主导微生物源存在显著差异

红树林柱以沉积物相关群落（9.98%）、根际微生物（8.58%）和土壤细菌（6.69%）为主；无红树林柱则以淡水沉积物（17.52%）、沉积物（10.28%）和海相沉积物（9.29%）为主；红树林湿地测压管由废水相关微生物主导（13.52%），人类肠道源位列第三（6.70%）。聚合为八个生态类别后，每个湿地带均有独特的主导类别：淡水井为植被与土壤源（45.1%），潮上带为海洋与河口感潮源（50.7%），过渡带为宿主相关与粪便源（45.6%），红树林湿地则为废水加宿主相关源（合计39.3%）。人为信号呈现空间异质性而非均匀分布，废水贡献在相距数米的测压管间从5.59%到19.89%不等，指向离散的局部输入。

3.3 红树林柱存在垂向结构，无红树林柱则无

在红树林柱中，未分化沉积物源从表层14.66%降至底层8.26%（β=?0.115%·cm^?1，R²=0.755，p=0.025）；海相源降幅最大，从12.48%降至1.97%（β=?0.177%·cm^?1，R²=0.697，p=0.039）；相反，根际、植物和土壤源随深度单调递增。在类别层面，植被与土壤源从表层14.8%升至底层38.3%（β=+0.375%·cm^?1，R²=0.753，p=0.026），而海洋与河口感潮源从27.8%降至4.9%（β=?0.377%·cm^?1，R²=0.712，p=0.034）。无红树林柱则无此可比结构，三种主导源均未随深度显著变化，表明植被覆盖而非生物扰动才是垂向结构化的驱动因素。

3.4 多元分析识别出两个主成分变异轴

层次聚类将样品分为三个良好分辨的簇：六个无红树林柱层段、六个红树林柱层段以及地下水测压管。PERMANOVA证实六个预定义样品组解释了大部分组成方差（伪F=4.140，R²=0.564，p=0.0001）。主成分分析显示，PC1和PC2分别解释总变异的19.1%和15.1%，前十个成分累计占90.3%。PC2是主要的生态兴趣轴：它将红树林柱层段从表层到底层有序排列，并将高废水负荷的红树林湿地测压管置于相反极端。

3.5 生态源类别解析样带空间结构

聚合生态类别提供了该系统空间结构的最清晰总结：植被存在与沉积物微生物源池组成相关（红树林柱植被与土壤源为无红树林柱的1.5倍，淡水源仅为其一半）；潮汐梯度位置与海相组分占比相关；靠近居民点和池塘基础设施则与人为组分占比相关。值得注意的是，没有任何类别随向海位置单调变化，表明该系统并非简单的线性海陆梯度，而是由植被覆盖、局部潮汐淹没频率和离散人为输入共同作用的镶嵌格局。

四、讨论与结论总结

讨论部分强调，SourceTracker2检测的是源环境的组成印迹而非活跃繁殖体本身，但水文输送机制、人为信号的空间异质性以及富集类别的合理性共同支持将其解释为活跃繁殖体供应。α多样性与组成变化的解耦警示不应将多样性指数作为海岸生态系统健康指标。植被依赖的垂向结构化表明红树林恢复不仅改变表层沉积物微生物组，还重塑整个沉积柱至至少60厘米深的组织方式，这对蓝碳核算具有直接意义。人为信号在沉积柱中的衰减提供了红树林沉积物过滤服务的首个定量指示。研究局限性包括单次湿季采样、单个断面的伪重复设计、缺乏同步理化测量等，因此结论应视为空间明确的基线和可检验假说，而非普适性描述。

研究结论部分翻译如下：本研究获得了四项数据支撑的结论：（1）微生物源谱的香农多样性在整个梯度上均匀一致（3.15±0.26；CV 8.4%），样品组间无显著差异（H=7.383，p=0.194），而源组成高度结构化（PERMANOVA R²=0.564，p=0.0001）；因此多样性指数在该系统中是不敏感的变化指标。（2）红树林植被区与相邻无植被区沉积物携带不同的源特征，植被与土壤源在红树林柱中富集（28.8%对20.1%，p=0.035），淡水源在无植被柱中富集（30.0%对15.3%，p=0.002）。（3）仅在有植被柱中出现源组成的垂向结构化，海相源表面富集（β=?0.377%·cm^?1，p=0.034），植被源深层富集（β=+0.375%·cm^?1，p=0.026），而无植被柱垂向均匀。（4）废水和宿主相关源在红树林湿地测压管中合计贡献39.3%，具有显著的点间异质性，表明离散的局部输入；该信号在沉积柱中减弱约两到三倍。