原文来源：Rongstock R, Li H, Lehmann A, et al., 2026. Global change factors differ in effect when acting alone and in a multi-factor background. Nature Communications, 17, 425-436.
https://doi.org/10.1038/s41467-025-68155-9.
科学背景
人类活动正在深刻重塑生态系统过程。在土壤生态系统中，多种全球变化因子往往同时存在并相互作用。与土壤相关的人为影响主要源于城市化、农业和工业活动，由此产生了生物、化学和物理的多种因子，包括气候变暖、干旱、盐渍化、富营养化以及污染物累积等。这些因子并非独立发挥作用，其组合效应可能放大、削弱甚至逆转单一因子的影响，从而限制了基于单因子实验的结论向复杂真实情境外推。
已有研究显示，全球变化因子的叠加往往带来更强的负向生态效应，然而我们对多因子胁迫背景下因子主导性的认识仍然有限。现有实验多局限于少数因子或采用因子数量梯度设计，难以回答在多因子胁迫背景下究竟哪些特定因子在主导制约土壤功能，从而影响我们对多重全球变化因子效应的机制理解。
城市生态系统往往同时暴露于多种全球变化因子，故而具有高度复杂性，正逐渐成为生态保护与恢复研究中的焦点。其全球变化因子组合方式复杂且空间异质性强，使城市土壤修复面临独特挑战。因此，优先识别并消除在多因子情境中产生最强负效应的因子，将为制定针对性环境管理策略和生态修复方案提供重要依据。
科学内涵
研究思路：
该研究采用减法设计，以弥补传统单因子添加以及因子数量梯度设计的不足。在设置控制组和单因子添加处理（+）的基础上，同时比较了六因子全添加（All）以及在六因子背景下去掉某一因子（All−i）的效应，从而将同一因子置于独立作用与多因子胁迫背景两种情境中进行比较。微宇宙实验土壤为白浆土（Albic Luvisol），采集自柏林自由大学草地实验站，处理方式包括：对照、六个单因子（升温、干旱、盐分、氮沉降、轮胎磨损颗粒、表面活性剂）、六因子叠加，以及六个“All−i”组合，总计14个处理方式，每种处理设10个重复。响应变量包含三种胞外酶活性（β-D-葡萄糖苷酶、酸性磷酸酶、N-乙酰-β-葡萄糖胺苷酶）、置入木牙签的分解率、土壤团聚体平均重量直径、水稳性团聚体比例、滴水穿透时间。统计方法上，该研究采用单因素方差分析与Tukey检验判断处理对各响应变量的统计显著性，采用标准差归一化（z-score）后的效应量配合自举法（n=100）呈现添加效应与去除效应的对应关系；利用冗余分析（RDA），将各土壤响应变量的总方差划分为三个组成部分，即独立方差（单因子效应）—共享方差（多个因子共同解释的交互影响）—残差（不可解释部分，通常为生物噪声或测量误差）。实验设计见图3（文中所有图表编号为原文编号）。
核心发现：
（1）不同全球变化因子单独作用时会存在不同的效应，可呈正向、负向和中性的单一因子效应。盐分对β-D-葡萄糖苷酶活性、酸性磷酸酶活性和分解过程均产生负面影响（见图1a、e、b、f、d、h）。盐分处理中，水稳性团聚体的百分比（见图2b，e）有所下降，但团聚体平均重量直径未受影响（见图2a，d）。变暖增强了N-乙酰-β-氨基葡萄糖苷酶活性（见图1c，g）和水稳性团聚体。干旱的加入对土壤团聚体平均重量直径产生了负面影响，此外，干旱还增加了滴水穿透时间（见图2c，f），然而干旱却增加了分解速率（见图1d，h）。微塑料和表面活性剂有增强β-D-葡萄糖苷酶活性的趋势（见图1a、e）但未呈显著影响（见表1）。氮沉积没有产生可测量的影响。
研究认为盐分可能通过多种途径降低酶活性，包括微生物细胞的渗透胁迫导致细胞脱水并释放胞内酶、改变酶蛋白活性位点的离子构象或者特定离子毒性导致微生物生长所需养分失衡，从而影响后续酶合成。盐分影响水稳性团聚体的可能机制则包含盐分上升致团聚体崩解或是抑制了参与团聚体形成与稳定的真菌活动。没有水分限制的条件下，变暖的单因子正效应被归因于能够提高微生物活性，如使真菌菌丝生长并增强黏结物质的分泌，从而提高水稳性团聚体。干旱对土壤团聚体平均重量直径的负效应体现为宏团聚体损失，干旱还延长了滴水穿透时间，但并未影响水稳性团聚体。单独施加微塑料、表面活性剂和氮沉降时并未产生显著影响。检测到微塑料和表面活性剂对β-葡萄糖苷酶活性具有正效应趋势，可能是由于这些添加物被某些特定微生物作为碳源利用，从而导致微生物群落结构发生变化。化学表面活性剂十二烷基硫酸钠也可被微生物快速矿化，可能在6周培养期后已被降解。氮沉降的效应不明显则可能由于研究中采用的添加剂量较低所致。
（2）与对照组相比，六因子联合处理的综合作用对大多数响应变量而言始终表现出显著的负面影响；单一因子的影响强度低于六因子联合作用的影响。单一干旱因子的影响是个例外，平均重量直径和滴水穿透时间在六因子联合处理下，与干旱单因子处理相比，并未表现出显著不同的效应。
该研究认为环境变化导致的新的变化因子组合（如土地利用变化）往往会使多种响应变量的效应方向趋于一致，即使某些单一因子本身表现的效应方向不同。然而，将两个以上全球变化因子结合的实验仍然较为少见，而且在预测多重全球变化因子效应时，不仅因子数量重要，因子的类型及其差异性同样发挥作用。
（3）当从全体因子背景中扣除某些因子时，其影响可能受到该因子本身、其与其他因子的相互作用或无法解释的变异性（生物噪声或测量误差）的影响。该研究将这些变异来源表示为共享方差和残差方差，其中残差方差包含了测量误差或生物噪声。对于酶活性，因子之间的共享方差（相互作用）普遍较低（0～12%），而未解释的方差比例较高（35.1～61.5%），指示生物噪声对该模式有较强影响。在分解、土壤团聚体平均重量直径和滴水穿透时间方面，研究显示因子间存在相互作用（共享方差，26.1～39.8%），且分解和土壤团聚体平均重量直径受生物噪声或测量误差的影响较小（分别为10.4%和27.0%）。然而，生物噪声或测量误差对滴水穿透时间的影响贡献更大（46.7%）。
图1 酶活性（β-D-葡萄糖苷酶、N-乙酰-β-氨基葡萄糖苷酶、酸性磷酸酶）和分解受多种全球变化因子影响。
图2 平均重量直径、水稳性团聚体和滴水穿透时间受多种全球变化因子影响。
图1a–d和图2a–c的左侧部分（白色背景）显示了对照组和单因子处理组的结果。虚线表示对照组的平均值。圆点代表测量值（每个处理组有10个重复），黑色实心圆点表示平均值，穿过该点的线条表示95%置信区间。图1和图2的右侧部分比较了六因子联合处理组合（虚线表示平均值）与从组合中减去该单因子后的处理组合。采用冗余分析（RDA）和方差分解量化每个因子的相对贡献及其相互作用，以及残差（生物噪声或测量误差），结果显示在图1e–h和图2d–f的右上角。为了将因子添加和去除处理的影响方向可视化，将因子添加的效应量（通过自举法计算，n=100）绘制在x轴上，将因素去除的效应量绘制在y轴上。不同的全球变化因子用不同的颜色表示。位于红色对角线上的因子效应表明该因子在多因子背景下的效应方向和强度与单因子效应预期的一致。位于灰色阴影区域的因子效应表明，该因子与其他因子结合后产生的负向影响比单因子效应的预期更大。深灰色区域（右上角）表示单因子效应为正，但与其他因子结合后却产生负向影响（单因子效应量为正，去除该因子后的效应量也为正）。原始数据以Source Data文件的形式提供。
从全体变化因子背景中去除盐分因子后，除土壤团聚体平均重量直径外所有响应变量均有显著改善。在酶活性方面，盐分独立解释了N-乙酰-β-氨基葡萄糖苷酶总方差的18.4%、β-D-葡萄糖苷酶活性的13%、酸性磷酸酶的27.7%；盐分独立解释了分解效应总方差的47.0%，水稳性团聚体的31.4%、滴水穿透时间的6.7%。与全体因子都存在的处理组相比，去除变暖的多因子处理组显示出的更高的酶活性。变暖因子独立解释了土壤酶活性总方差的12.7～15.9%。然而，在所有其他响应变量中，变暖因子解释的总方差不到1%。消除干旱改善了除水稳性团聚体外所有响应变量的表现，进而推论多因子背景下的干旱效应为负面。干旱不是解释分解的主要因子（1.5%），而盐分解释了大部分方差。干旱解释了土壤团聚体平均重量直径方差的44.3%，表明其对该响应变量有强烈且稳健的影响。干旱对滴水穿透时间的影响最大，解释了其14.2%的方差。氮沉降、表面活性剂和微塑料的去除均未对响应变量产生影响。
（4）一些全球变化因子在单因子添加与其从多因子背景中移除时表现出不同的效应大小和方向。
盐分作为单一因子时多表现为负效应，而从多因子背景中去除盐分则会产生相应的正向影响。对N-乙酰-β-氨基葡萄糖苷酶活性而言，单因子盐分无显著影响，但从多重全球变化因子背景去除盐分后，该酶活性显著升高（效应量=0.30，95%置信区间：0.21～0.39，p<0.001）。与对照组相比，盐分显著降低酸性磷酸酶活性（效应量=-2.24，95%置信区间：-3.45～-1.47，p<0.001），而在多因子处理中去除盐分后，该酶活性明显升高（效应量=2.83，95%置信区间：2.07～3.64，p<0.001）。盐分也降低了β-D-葡萄糖苷酶活性（效应量=−0.76，95%置信区间：−1.12～−0.38，p<0.001），而在多因子背景中去除盐分后，该酶活性显著升高（效应量=0.87，95%置信区间：0.48～1.31，p<0.001）。盐分对木材分解有负向影响（效应量=−1.16，95%置信区间：−0.19～−0.13，p<0.001），而去除盐分后，分解率虽有提升，但幅度较小（效应量=0.47，95%置信区间：0.32～0.53，p<0.001）。在土壤结构和性质方面，盐分的影响呈中性或负向，显著降低了水稳性团聚体的比例（效应量=−11.91，95%置信区间：−16.71～−6.08，p=<0.001）。在去除盐分的多因子处理组中，水稳性团聚体含量与对照组相当。由于对照组滴水穿透时间较低，无法明显观察到盐分对该指标的单因子影响，但在多因子背景下去除盐分后，水滴吸收时间有所延长（效应量=17.56，95%置信区间：5.26～29.63，p=<0.050）。
单独考虑变暖因子时，其影响为正向或中性，但从多因子背景下去除变暖因子后，与微生物活性相关的响应变量均有所增加。单因子变暖并不影响β-D-葡萄糖苷酶活性，但从多因子背景下去除变暖因子后，该酶活性显著升高（效应量=0.87，95%置信区间：0.46～1.22，p<0.001）；酸性磷酸酶呈现相同趋势，单因子变暖无显著效应，在去除变暖因子后该酶活性显著增加（效应量=2.14，95%置信区间：1.18～2.88，p<0.001）。变暖处理显著提高了N-乙酰-β-氨基葡萄糖苷酶的活性（效应量=0.23，95%置信区间：0.08～0.38，p<0.0245），从多因子背景中移除变暖因子后，该酶活性也有显著增加（效应量=0.28，95%置信区间：0.17～0.36，p<0.001）。变暖处理组的木材分解速率与对照组相比有显著提高（效应量=0.17，95%置信区间：0.13～0.20，p<0.001），但从多因子背景中去除变暖因子后，分解速率未发生显著变化。在土壤结构指标中，仅水稳性团聚体对变暖有响应，变暖处理显著提高其含量（效应量=7.75，95%置信区间：3.39～13.35，p=0.0205），而多因子处理组中去除变暖因子后，并未观察到显著影响。
干旱的单因子效应总体偏中性，但从多因子背景去除后影响呈正向。单因子干旱不影响酶活性，但与全体因子处理组相比，去除干旱后β-D-葡萄糖苷酶（效应量=0.68，95%置信区间：0.31～1.08，p=<0.00671）、酸性磷酸酶（效应量=2.14，95%置信区间：1.18～2.88，p=<0.001）及N-乙酰-β-氨基葡萄糖苷酶（效应量=0.43，95%置信区间：0.30～0.55，p=<0.001）的活性均显著升高。与对照组相比，单因子干旱处理组的木材分解率显著提高（效应量=0.35，95%置信区间：0.25～0.40，p<0.001），去除干旱的多因子背景下，木材分解率有轻微正向变化（效应量=0.09，95%置信区间：−0.07～0.15，p<0.001）。在土壤结构方面，单因子干旱显著降低了土壤团聚体平均重量直径（效应量=−0.50，95%置信区间：−0.61～−0.33，p<0.001），而在多因子背景中去除干旱后，土壤团聚体平均重量直径明显增加（效应量=0.64，95%置信区间：0.47～0.84，p=<0.001）。干旱显著延长了滴水穿透时间（效应值=5.45，95%置信区间：3.25～7.80，p=<0.001）。当剔除干旱因子后，其他全球变化因子组合的效应与对照组无显著差异。
图3 研究设计与可能结果
a 处理因子包括变暖、干旱、盐分、氮沉降、微塑料（轮胎磨损颗粒）和表面活性剂。设置10个无处理重复作为对照。每个单因子添加（+）、六因子全添加（All）以及六因子组合中逐一去除某一因子（All−i）的处理均设置10个重复。实验共包含140个实验单元。b 实验单元为装有30g土壤的Falcon离心管，每根管中放入一根预先称重的木牙签，用于测定分解速率。培养6周后，测定各响应变量。c 因子添加与因子去除效应方向的可视化：在x轴上绘制因子添加的自举法效应量（n=100），在y轴上绘制因子去除的效应量。
启示1：全球变化因子在单独添加时的效应大小和方向，无法预测其在多因子相互作用背景中的真实生态作用，甚至在效应方向上也可能发生改变。将某一因子从多因子组合中移除时，引发的响应不仅取决于该因子的独立效应，还受到与其他因子之间交互作用的调控。因此，基于单因子实验推断多重全球变化情境下的生态响应明显具有局限性，甚至可能在产生误导。该研究强调，只有将全球变化因子置于多因子背景中加以解析，才能准确评估其在真实生态系统中的功能性作用及其对生态过程的调控意义。
启示2：在多重全球变化因子难以同时缓解的现实条件下，环境管理和生态修复的关键在于识别并优先干预那些在多因子情境中对生态功能产生主导性约束的因子。正如该研究结果显示，不同因子在多因子组合中对土壤功能的贡献并不均衡，一些因子的移除可显著改善多项生态响应，而其他因子的去除则几乎不改变系统状态。因此，城市土壤修复和减压管理不应仅依据单因子效应强度或因子数量，而应考虑目标因子在多因子背景下移除后的效应变化，确定优先干预对象，以在资源受限条件下实现最大的修复效益。
——解读人：张原
——责任编辑：童曼
