亏缺灌溉耦合生物炭增强盐渍土机械团聚体稳定性：细菌群落确定性装配与网络稳健性的关联

《Environmental Technology & Innovation》：Biochar and deficit irrigation are associated with enhanced saline soil mechanical stability and shifts in bacterial assembly

　　

生物炭（biochar）与亏缺灌溉（deficit irrigation, DI）联合处理与微生物群落装配（microbial assembly）中更强的确定性影响（deterministic influence）相关。生物炭（30 t ha^-1）与增强的微生物网络稳健性（network robustness）并降低脆弱性（vulnerability）相关。装配中更强的确定性影响和网络稳健性与团聚体稳定性（aggregate stability）相关。

研究背景方面，土壤盐渍化威胁干旱和半干旱区农业可持续性，钠诱导的团聚体破碎会恶化土壤结构，并损害微生物群落和养分循环。内蒙古河套灌区蒸发/降水比极高，表土盐分积累和结构脆弱问题突出。亏缺灌溉（deficit irrigation, DI）可节水，但减少淋洗，可能加剧盐分积累和土壤微生物渗透胁迫；生物炭（biochar）可通过吸附 Na⁺ 和改善土壤理化性质缓解盐害。以往研究多关注生物炭单独的直接理化效应，较少考虑水分管理与微生物途径；而微生物是团聚体稳定性的关键驱动者。DI 与生物炭交互对微生物群落装配和共现网络架构的级联效应尚不清楚。研究人员因此假设，DI 与生物炭整合可形成协同策略：生物炭缓解 DI 可能加重的盐胁迫，DI 的受控水分条件优化生物炭功能和微生物活性；该协同不仅通过加和理化效应，还主要通过重构土壤微生物组增强团聚体稳定性。研究以归一化随机性比率（normalized stochasticity ratio, NST）<0.5 表示确定性装配增强等框架检验假设，并提出：水-炭协同通过缓解盐胁迫，使土壤微生物组朝确定性装配和稳健网络转变，从而增强团聚体稳定性。研究问题包括：DI 与生物炭交互是否使微生物群落装配由随机向确定性转变；若如此，是否与共现网络稳定性增强相关；这些微生物性质是否比灌溉或生物炭输入本身与团聚体稳定性关联更强。研究人员进一步预测，若水-炭协同通过确定性装配和网络稳健性发挥作用，则微生物网络性质应与平均重量直径（mean weight diameter, MWD）相关且独立于处理因素，并且改善应限于机械稳定性而非水力稳定性。

研究人员开展了一项为期2年的田间试验，地点位于内蒙古巴彦淖尔五原科技小院试验站，时间为2022和2023年向日葵生长季，土壤为盐碱砂壤土。试验采用两因素全因子设计：灌溉水平为充分灌溉（full irrigation, FI，100% ETc）和亏缺灌溉（DI，60% ETc），生物炭用量为0、15和30 t ha^-1，共18个小区、3次重复。采样深度0–20 cm，并在幼苗、开花和成熟期采样。研究表明，DI+B30处理与盐渍土机械团聚体稳定性增强、细菌装配确定性影响增强和细菌网络稳健性增加相关；水稳性团聚体未表现平行改善；真菌群落未显示显著响应。该研究的意义在于将水分管理、生物炭改良、细菌群落装配、网络稳健性和机械团聚体稳定性联系起来，同时强调观察性田间设计限制因果推断，相关路径仍需操控实验验证。

开展研究的关键技术方法可概括为：2年田间试验，样本队列来源为2022—2023年内蒙古巴彦淖尔五原科技小院试验站盐碱砂壤土向日葵根区0–20 cm土壤；采用FI/DI与生物炭0、15、30 t ha^-1两因素全因子设计，18小区3重复；用干筛和湿筛分离团聚体并计算MWD、几何平均直径（geometric mean diameter, GMD）和团聚体破坏率（percentage aggregate destruction, PAD）；进行16S rRNA V3–V4和ITS2测序，DADA2获得扩增子序列变异（amplicon sequence variants, ASVs）；用NST、中性群落模型（neutral community model, NCM）和iCAMP框架分析装配；用MENA构建共现网络；用随机森林（random forest, RF）和结构方程模型（structural equation modeling, SEM）分析关联。

结果部分：3.1 Soil aggregate stability（土壤团聚体稳定性）表明，DI+B30增加了机械稳定大团聚体比例，使MWD较DI+B0增加14.4%；>2 mm和0.5–2 mm粒级质量比例提高，FI下B30未显著改变机械团聚体分布或MWD；DI+B15未显著提高MWD，且显著低于DI+B30。水稳性方面，DI下B30使WR_0.25由35.2%增至38.6%，但不显著；FI下B30未改变WR_0.25；B30在DI和FI下分别使水稳性团聚体GMD降低12.0%和6.9%；PAD无显著差异。结论是生物炭在亏缺灌溉下增强机械团聚体稳定性，而水稳性响应不一致。

3.2 Restructuring of microbial community and diversity（微生物群落与多样性重构）中，3.2.1 Composition of soil bacterial and fungal communities（土壤细菌和真菌群落组成）显示，细菌群落存在较高比例未分类类群，优势属为 Sphingomonas、Nocardioides 和 Lysobacter；真菌以 Ascomycota 和 Basidiomycota 为主，属水平包括 Aspergillus、Mucor、Acaromyces、Cercospora 和 Pseudogymnoascus。生物炭显著改变细菌组成，DI+B30使 Actinobacteria 相对丰度较DI+B0增加13.3%，使 Chloroflexi 减少36.6%，对真菌影响较小。Sankey分析显示，DI下生物炭使核心OTU较B0增加51.9%（B15）和60.2%（B30）；FI下B15使核心OTU减少55.3%，B30增加63.7%，表明存在性状过滤过程。

3.2.2 Shift toward deterministic assembly（向确定性装配转变）中，NCM、NST和iCAMP分析表明，DI+B30细菌NCM解释方差R²=0.584，高于DI+B0；DI+B30的Nm大于DI+B0，DI+B0扩散限制更强。NST由DI+B0的0.58±0.02降至DI+B30的0.43±0.02，向更强确定性影响转变。iCAMP显示，DI+B0中同质选择为28.4%、异质选择14.0%、扩散限制15.6%、同质化扩散9.5%、漂移32.5%；DI+B30中分别为13.6%、36.9%、3.9%、13.6%和18.0%。漂移仍占18.0%，说明随机和确定性过程同时存在。冗余分析中7个土壤因子共同解释42.7%的细菌群落变异，EC_1:5最强，R²=0.21，其次为土壤含水量（soil water content, SWC）R²=0.17和土壤有机碳（soil organic carbon, SOC）R²=0.14；DI+B30样本向低EC_1:5和高SOC分离。NST与EC_1:5正相关r=0.82，与SWC负相关r=-0.74，与SOC负相关r=-0.69。结论是低盐、适度水分和有机碳增加与确定性装配相关。

3.2.3 Complexity and stability of bacterial and fungal networks（细菌和真菌网络复杂性与稳定性）中，MENA网络分析显示，生物炭增加细菌网络节点和边，DI+B30平均度和聚类系数显著高于DI+B0；真菌网络参数无显著变化，P=0.06–0.12。DI+B30使细菌网络脆弱性较DI+B0降低44.2%，真菌网络脆弱性降低32.9%；50%随机物种移除下细菌网络稳健性增加13.1%，正连接增加49.3%。真菌响应有限，表明网络稳健性及其与团聚体稳定的联系主要由细菌驱动。

3.3 Microbial properties as correlates of mechanical aggregate stability（微生物性质与机械团聚体稳定性的关联）中，RF模型解释MWD变异的72.5%，NST的%IncMSE最高，为21.5，其次为异质选择18.7、网络稳健性16.3、生物炭用量15.0和灌溉制度9.8；但变量重要性差异不大，交叉验证R²为0.67–0.76，说明模型存在不稳定性，结果应视为提示性。SEM拟合较好，CMIN/DF=1.81，CFI=0.97，RMSEA=0.04；水-炭交互与确定性装配β=0.75，与网络稳健性β=0.22；确定性装配与网络稳健性β=0.68；网络稳健性与MWDβ=0.40；确定性装配到MWD路径不显著β=0.07。原文认为结果与网络稳健性中介装配状态和团聚体稳定性的统计关系一致，但观察性设计限制因果推断，SEM应视为探索性。

讨论部分：4.1 Water-Biochar interactions are associated with shifts in microbial assembly（水-生物炭交互与微生物装配转变相关）指出，DI+B30与装配平衡向确定性影响转变相关，NST下降、异质选择增加，但漂移仍占18.0%，NCM拟合中等，说明并非完全体制转变。RDA显示EC_1:5、SWC和SOC是主要环境过滤器；生物炭在DI和FI下均降低EC_1:5和交换性Na⁺，但仅DI+B30增加MWD，说明单独缓解盐分不足。PICRUSt2功能预测仅为探索性，不能证明EPS合成或基因表达。高剂量生物炭可稳定能量供应，减少随机漂移；异质选择通过微生境异质性竞争排除增强。FI+B30增加细菌OTU但提高孔隙水EC，再动员Na⁺，减少>2 mm水稳团聚体18.4%。B15在DI下NST=0.51，未显著改善MWD，提示剂量依赖而非离散阈值。

4.2 Microbial network robustness links biochar–water synergy to aggregate stability（微生物网络稳健性连接生物炭-水协同与团聚体稳定性）表明，生物炭重塑微生物共现网络并增加复杂性。DI+B30虽使MWD增加14.4%，但水稳性无对应改善，水稳性GMD降低12.0%，反映干筛和湿筛机制不同，即改善的是机械抗性而非水力稳定性。网络模块度和路径效率变化可能反映跨模块资源交换。Nocardioides 与 Proteobacteria 的稳定正关联可能支持营养循环相关双功能子网络。脆弱性降低44.2%说明生物炭缓冲盐诱导的连接破坏，保留促进团聚体形成的微生物功能。FI+B30增加聚类但降低模块度，未改善水稳大团聚体；B15稳健性增益小。真菌无统计显著响应，主要因菌丝对盐渍土渗透胁迫敏感、生物炭可能偏好细菌定殖、ITS2分辨率有限，因此网络稳健性和团聚体稳定联系主要归因于细菌。

4.3 Microbial-engineering synergy: from functional taxa to smart remediation（微生物工程协同：从功能类群到智能修复）指出，Actinobacteria 对DI+B30强正响应，相对丰度增加13.3%，并呈DI+B30 > FI+B30 > DI+B0梯度；其作为耐胁迫生态工程师可能贡献团聚体稳定性。该嗜碱类群具有活跃的H⁺-ATPase质子运输系统。与 Bacteroidota 共现可能与生物炭芳香化合物降解和necromass碳池形成有关，疏水胞外多糖可能增强0.5–2 mm团聚体抗侵蚀。真菌响应有限，进一步支持主要细菌驱动。

研究结论翻译：为期两年的田间试验表明，亏缺灌溉与高剂量生物炭联合施用与盐渍土机械稳定性增强及细菌微生物组重构相关。具体而言，该组合与细菌装配中更强的确定性影响、异质选择贡献增加和更高的细菌网络稳健性共变。在随机森林和SEM分析中，这些微生物性质与平均重量直径统计相关，但预测因子重要性差异不大，应谨慎解释。水稳性团聚体的解耦响应表明，该稳定化限于机械抗性，并未改善水力稳定性。真菌群落也未显示平行响应。研究人员得出结论，水-生物炭协同在统计上与细菌装配和网络性质相关，这些性质又与机械团聚相关。SEM分析应视为产生假设而非检验假设。该田间试验的观察性设计排除了因果归因，所提出的路径需要通过针对性操控实验进行验证。