高寒采矿基质首个生长季早期植被恢复启动土壤物理生境重建与微生物群落再组装

《Plants》：Early Revegetation Initiates Soil Physical Habitat Reconstruction and Microbial Community Reassembly Within the First Growing Season in an Alpine Mining Substrate

　　

高寒矿区生态恢复受低温、贫瘠基质和有限水分供应的制约，然而植物、土壤与微生物之间的协调早期响应仍所知甚少。研究人员在西藏一处高寒砂石矿区开展了一项为期90天的原位植被恢复实验，选用三种原生多年生禾草——垂穗披碱草(Elymus nutans)、冰草(Agropyron cristatum)和中华羊茅(Festuca sinensis)，设置单播与混播处理。研究人员对植物生长与生理、土壤理化性质、养分、酶活性及微生物群落进行了测定评估。所有三种先锋物种在第一个生长季内均建立了可测量的生物量，尽管不同种植配置间的生长和生理响应存在差异。土壤毛管孔隙度(capillary porosity, CP)、含水量及持水能力表现出最一致的改善，而酶活性和养分响应则更具处理依赖性。细菌和真菌群落经历了显著重组，推断的共现网络(co-occurrence network)表现出时间拓扑重构以及对模拟靶向节点移除的敏感性。Mantel分析显示，细菌群落变异与CP、非毛管孔隙度(non-capillary porosity, NCP)和总磷(total phosphorus, TP)显著相关，而真菌群落变异与所测土壤变量无显著关联。总体上，90天的先锋植被恢复在高寒采矿基质中启动了可测量的植物-土壤-微生物恢复。物理和水文性质较大多数养分指标表现出更清晰的早期响应，而土壤孔隙特征与细菌群落组成之间的观测关联提示，物理生境变化与微生物群落再组装之间存在潜在联系。

**研究背景与存在问题**

青藏高原高寒生态系统在全球具有重要生态地位，但海拔超过3500 m、年均气温低、降水有限等气候约束会减缓养分循环并限制退化高寒土壤中植物的建植。不断扩张的矿业活动造成大面积地表扰动，包括植被覆盖丧失、土壤结构破坏和养分库耗竭，仅依赖自然演替难以在管理相关时间尺度上实现快速恢复。已有恢复成效评估多集中于地上盖度或生物量等指标，对地下过程关注不足，且很少在同一研究中同时考察植物、土壤与微生物的协调早期响应。尤其是，在高寒矿区早期恢复阶段，植被建植引起的土壤物理性质变化是否与细菌和真菌群落重组相关联，仍不清楚。为此，研究人员提出三个假设：H1，原生先锋禾草可在第一个生长季内建立可测量生物量并维持生理活性；H2，早期植被建植快速改善采矿基质物理和水文性质，而总体养分状况变化相对有限；H3，先锋植被建植诱导土壤细菌和真菌群落发生可测量的重组，且与基质微环境特别是土壤物理结构变化相关联。

**研究内容与总体结论**

研究人员在西藏拉萨市达孜区唐嘎乡卡尔青建筑用砂石矿（29°50′48″ N，91°33′12″ E）开展了一项为期90天的原位恢复实验，选用垂穗披碱草、冰草和中华羊茅三种原生多年生禾草，设置3个单播、4个混播处理及1个裸地对照(CK)，每处理3个重复样方。结果表明，90天先锋植被恢复足以在高寒采矿基质中启动可测量的植物-土壤-微生物恢复；土壤物理和水文性质呈现最一致的早期响应，而养分和生化指标响应更具处理依赖性；细菌和真菌群落发生显著重组，细菌群落变异与毛管孔隙度(CP)、非毛管孔隙度(NCP)和总磷(TP)显著相关。这一发现提示，物理生境变化与微生物群落再组装之间存在潜在联系，对高寒矿区恢复实践具有参考意义。

**关键技术方法**

研究采用的主要方法包括：原位植被恢复实验设计（三种原生禾草、8个处理、每处理3个重复，样方3 m×3 m）；完整土柱法测定土壤水分、孔隙度、持水能力等物理性质；16S rRNA基因V4区和真菌ITS1区扩增子测序（MGISEQ-2000RS测序仪）分析微生物群落组成；DESeq2差异丰度分析；基于Spearman相关构建共现网络并进行鲁棒性分析；Mantel检验与Pearson相关分析环境因子和微生物群落结构的关系。土壤样品采集自0–20 cm表层，每样方按五点采样法混合为一个复合样品；研究对象为高寒砂石矿采矿基质。

**研究结果**

**2.1 植物生长表型与生理指标**：通过单因素方差分析和Tukey HSD检验比较各处理，发现光合色素含量呈物种特异性响应；C1处理垂穗披碱草地上生物量最高（76.3 g·m^?2），C4处理地下生物量最高（169.7 g·m^?2）；Pearson相关分析显示叶绿素a(Chl a)与丙二醛(MDA)含量呈显著负相关，表明较高Chl a与较低膜脂过氧化水平相关。

**2.2 土壤理化性质**：通过完整土柱法和理化测定发现，种植处理土壤含水量(SWC)和CP显著高于CK，CP为CK的1.42–1.54倍，而NCP较CK低27.43–54.86%；饱和含水量、毛管持水量和田间持水量均显著高于CK；土壤容重在各种植处理中均低于CK，其中C4和C7差异显著。物理和水文性质表现出相对一致的响应，而化学性质更具处理特异性。

**2.3 土壤养分与酶活性**：通过比色法和分光光度法测定发现，蔗糖酶活性在C1最高（38.14 mg·g^?1·24 h^?1），约为CK的40.6倍；脲酶活性在C1、C3、C6显著高于CK；碱性磷酸酶仅在C1显著高于CK；部分处理铵态氮(NH₄⁺-N)、硝态氮(NO₃^?-N)和速效钾(AK)显著高于CK；各种植处理TP均高于CK（1.27–1.94 g·kg^?1 vs 0.98 g·kg^?1），但土壤有机碳(SOC)、速效磷(AP)和全氮(TN)与CK无显著差异。

**2.4 土壤微生物多样性与分类组成**：通过α多样性分析和PCoA/PERMANOVA发现，多数种植处理的细菌和真菌Chao1丰富度显著高于CK，而Shannon多样性无显著差异；PERMANOVA显示细菌和真菌群落组成在处理间差异显著，但PERMDISP提示真菌的组内离散度差异可能影响结果解读；优势细菌门为放线菌门(Actinomycetota，13.0–30.3%)、变形菌门(Pseudomonadota，16.8–28.6%)和酸杆菌门(Acidobacteriota，3.4–10.3%)，真菌以子囊菌门(Ascomycota，77.0–90.5%)为主。

**2.5 差异丰度微生物类群与时间共现网络动态**：通过DESeq2分析鉴定出101个差异丰度属，其中85个在种植组富集、16个耗尽；Fusarium在C1、C2、C4中较CK显著降低83.68–88.18%；含潜在植物促生菌株的Sphingomonas、Bradyrhizobium和Microvirga等属在种植处理中显著富集；基于Spearman相关构建的阶段特异性共现网络显示，节点数从30天的62个增至90天的103个，正相关比例从88.89%升至95.79%，但网络密度和模块性下降，靶向节点移除下相对自然连通度快速下降。

**2.6 与早期微生物群落重组相关的土壤环境因子**：通过Mantel检验和Pearson相关分析发现，细菌群落结构与CP、NCP和TP显著相关（p < 0.05），而所有测定土壤变量均与真菌群落结构无显著关联，表明细菌群落的早期重组与土壤孔隙特征和TP变化联系更为紧密。

**讨论总结与研究结论**

讨论指出，90天恢复期内土壤物理和水文性质响应迅速而养分指标变化有限，表明基质功能恢复存在异步性：根系通过生物孔隙形成和颗粒位移重构孔隙结构，但团聚体稳定性与孔径分布未直接测定，结构重建程度尚难全面判定。Fusarium属含多种土传病原菌，其丰度下降提示植被建植可能改变潜在病原菌组成，但致病性未被直接评估；Sphingomonas、Bradyrhizobium和Microvirga等属中部分菌株具有促生潜力，其功能角色仍需实验验证。共现网络的拓扑重组和节点移除敏感性表明早期网络结构脆弱，且网络边代表统计关联而非经验证实的生态互作。细菌和真菌群落对环境因子响应的差异提示二者遵循不同的组装轨迹。研究结论强调：一个生长季（90天）的先锋植被恢复足以在高寒采矿基质中启动可测量的植物-土壤-微生物恢复；土壤物理和水文性质呈现最清晰的早期响应，而总体养分状况变化有限；细菌和真菌群落发生显著重组，但推断网络仍对模拟靶向节点移除敏感；细菌群落变异与CP、NCP和TP显著相关，其方向性和因果性尚待确定；90天恢复仅启动了生境重建而非完成生态系统恢复，单播与混播配置的长期优劣无法在如此短的时间尺度上判定。