从群落演替速率与全球变化强度关系推断群落组成层面的适应性

《Environmental Microbiology Reports》：Inferring Adaptation at Community Composition Level From the Relationship Between Succession Rate and Global Change Intensity

　　

为响应人为环境变化，生物群落常发生组成变化，传统上认为这种变化代表适应性。然而，该假说及其机制仍鲜有深入研究。本研究在欧亚草原开展了一项为期10年的氮沉降实验，在围封与割草两种管理策略下设置了9个氮添加强度。随着强度从0–20 g N m?2 yr?1增加至50 g N m?2 yr?1，植物与土壤微生物群落的演替速率均呈先升后降的趋势，表明高强度的氮沉降超出了群落适应的能力范围。此外，群落组成变异的确定性组分演替速率与土壤理化性质的演替速率也随氮强度增加呈先升后降趋势，且二者呈正相关，表明适应性主要由确定性生态过程驱动。同时，割草促进了植物群落的适应，而围封促进了微生物群落的适应，二者均通过抑制随机过程实现。总体而言，为使生物群落在组成层面适应环境变化，应将变化强度控制在阈值以下，并实施合理策略，无论是确定性策略（如恢复土壤理化条件）还是随机性策略（如添加适应物种的种子）。

一、研究背景与科学问题

群落生态学的核心目标之一是探究人为环境变化与生态系统内生物之间的关系。这种关系存在两个对立面：一方面，环境变化可能过滤掉适应度较低的物种；另一方面，生物可通过三个层面的过程实现适应——生理适应（表型可塑性，时间尺度短于一代）、群落组成层面的适应（物种组成更替，时间尺度为数代至数十代）以及遗传适应（进化时间尺度）。以往研究多关注生理与遗传层面的适应，而对群落组成层面是否真正实现完全适应缺乏实验验证。现代生态学将慢性人为压力（如增温、氮富集）下的定向组成变化纳入生态演替框架，演替速率成为量化适应性组成更替的关键指标。理论上，物种库中不同物种占据不同生态位，环境变化后群落组成会相应调整，但物种库物种数量有限，群落组成可通过组成更替适应的环境变化范围也有限。因此，随着全球变化强度增加，群落时间变化速率可能呈现先升后降的非线性模式：强度极低时，轻微组成变化即可实现适应；强度增加时，组成变化加剧，演替速率显著上升；但若强度超过所有物种的适应范围，群落短期内达到非适应稳定状态，后续变化极小，整体时间周转速率下降。换言之，群落组成层面的适应受限于环境变化强度的阈值。基于群落时间变化速率（即群落演替速率）与全球变化强度的关系，可推断特定强度下群落组成的时间变化是否完全代表适应。此外，确定性过程（生态过滤、种间竞争）与随机过程（随机出生/死亡、扩散/定殖）均可能驱动群落组成动态，但哪种过程是组成适应的主要机制尚不清楚，特定管理策略能否通过促进确定性过程或抑制随机过程来加速群落适应也有待探究。

二、研究内容与结论

研究人员在内蒙古典型草原（代表欧亚草原大部分区域的植物区系与生态特征）开展了为期10年的氮沉降实验，设置9个氮添加强度（0–50 g N m^?2 yr^?1）与两种管理策略（围封与割草），共18个处理，每处理6个重复样方。研究同时调查植物群落组成与土壤微生物生物标志物（PLFA）图谱，回答三个关键问题：其一，随氮沉降强度从低到高增加，植物与土壤微生物群落演替速率是否先升后降，即群落层面适应是否在较低强度下发生而在高强度下受限；其二，群落组成层面的适应主要通过促进确定性过程还是抑制随机过程实现；其三，管理策略是否影响演替速率及群落适应，其关键机制是促进确定性过程还是抑制随机过程。

研究得出以下主要结论：植物与土壤微生物群落的演替速率与氮添加强度呈二次曲线关系，即随强度增加先升后降；确定性组分的演替速率同样与氮强度呈显著二次关系，而随机组分的演替速率大多无此规律；土壤理化性质（整合pH、速效氮含量与含水量）的演替速率也与氮强度呈显著二次关系，且与群落整体及确定性组分的演替速率显著正相关，但与随机组分无显著相关；割草促进了植物群落的适应（植物群落演替速率在割草下更高），围封促进了微生物群落的适应（微生物群落演替速率在围封下更高），二者均通过抑制随机过程实现。该研究的重要意义在于：首次从群落演替速率与全球变化强度关系的角度实验验证了群落组成层面适应的存在与阈值；揭示了确定性过程是组成适应的主导机制；阐明了管理策略通过调节随机过程影响适应速率的途径；为维持生物多样性与生态系统功能提供了将全球变化强度控制在适应阈值以下、并合理实施确定性或随机性策略的管理建议。

三、主要关键技术方法

研究采用时间衰减关系（TDR）量化群落演替速率，即对群落组成对数β相似性与对数时间距离进行线性回归，回归斜率即为群落组成演替速率。确定性组分与随机组分的量化采用改进的计算框架：以处理内重复样方在起始年份的平均Bray–Curtis相异度为参考基线，计算跨年相异度与参考基线的差值分别作为确定性组分与随机组分，再对各组分与时间距离回归求斜率作为相应演替速率。土壤理化性质演替速率以标准化后的pH、速效氮、含水量计算Bray–Curtis相异度后回归求得。此外，采用Chase（2010）的零模型方法计算标准效应值（SES），以验证实验方法的可靠性。统计分析包括限制性置换多元方差分析（PERMANOVA）、非度量多维尺度分析（NMDS）、双因素重复测量方差分析、单样本t检验及线性/二次回归，所有分析在R软件v4.3.2的vegan包中完成。

四、研究结果

（一）群落组成演替速率与氮强度的二次关系

限制性PERMANOVA显示管理策略、氮添加强度与处理时间对植物和土壤微生物群落组成均有显著影响。NMDS排序图显示不同氮强度与管理策略下群落演替轨迹存在差异。时间衰减关系显示大多数处理下群落组成呈现显著时间周转模式。最关键的发现是，两种策略下植物与土壤微生物群落的演替速率均与氮添加强度呈二次曲线关系（割草下植物群落除外），即随强度增加先升后降，峰值出现在20 g N m^?2 yr^?1左右，50 g N m^?2 yr^?1时急剧下降。

（二）确定性组分与随机组分的演替速率

群落组成变异的确定性组分在多数处理下呈显著时间周转模式，其演替速率与氮强度呈显著二次关系，与群落整体演替速率的模式一致。随机组分的演替速率则大多无线性或二次关系：割草下植物群落随机组分与氮强度呈显著线性关系，围封下微生物群落随机组分与氮强度呈显著线性关系，而围封下植物群落与割草下微生物群落的随机组分演替速率与氮强度无任何关系。零模型SES演替速率与氮强度呈边缘显著二次关系，进一步支持了实验方法的可靠性。

（三）土壤理化性质演替速率及其与群落演替速率的关系

双因素重复测量方差分析显示氮添加显著提高土壤速效氮含量、降低土壤pH，割草进一步提高速效氮含量并降低土壤含水量。土壤理化性质在所有处理下均呈显著时间周转模式，其演替速率与氮强度呈显著二次关系。土壤理化性质演替速率与植物（围封下）和微生物（围封与割草下）群落整体及确定性组分的演替速率显著正相关，但与随机组分的演替速率无显著相关。

（四）管理策略对演替速率的影响

在0–20 g N m^?2 yr^?1强度范围内，植物群落组成的演替速率在割草下显著高于围封，而植物群落确定性组分与随机组分的演替速率均在围封下显著高于割草。土壤微生物群落组成的演替速率在围封下显著高于割草，确定性组分无显著差异，随机组分在围封下显著低于割草。

五、讨论与结论

研究结果表明，植物与微生物群落的演替速率随氮强度增加先升后降，表明0–20 g N m^?2 yr^?1范围内群落可通过组成调整实现有效适应，而50 g N m^?2 yr^?1时适应能力崩溃。演替速率下降可归因于本地物种库的耗竭：欧亚草原土壤长期缺氮且呈中性或弱碱性，物种库中极少物种经历过并适应了最高强度氮添加导致的过量氮与土壤酸化，群落迅速稳定为物种贫乏状态，时间周转速率降低。

确定性过程而非随机过程是群落组成适应全球变化的主导机制。依据有三：其一，确定性组分演替速率与氮强度呈二次关系，与群落整体演替模式一致，而随机组分无此规律；其二，土壤理化性质演替速率与氮强度呈二次关系，且与群落整体及确定性组分演替速率显著正相关，与随机组分无相关，表明氮添加通过改变土壤环境对物种存活与表现施加强定向过滤；其三，零模型SES演替速率与氮强度的二次关系进一步支持了实验结果。氮添加显然通过改变土壤理化指标确定性筛选物种，而非促进或抑制物种的随机出生/死亡与迁移/定殖。

管理策略通过调节随机生态过程影响群落适应速率。割草下植物群落演替速率更高，其随机组分演替速率低于围封且实际为负值，因为割草每年移除地上植物并阻碍植物物种的随机定殖。围封下微生物群落演替速率更高，其随机组分演替速率低于割草且实际为负值，可能因为围封排除放牧干扰、促进地上生物量与凋落物积累，改变了微生物的扩散与定殖途径，从而降低土壤微生物组成变化的随机组分。随机出生/死亡与迁移/定殖过程会减缓定向适应过程，因此能够抑制这些随机过程的管理策略可促进生物群落对全球变化的确定性适应。

综上所述，从群落演替速率与全球变化强度关系的角度，研究人员成功推断出特定氮添加强度下生物群落是否通过物种组成变化实现适应，识别出确定性适应机制，并阐明了管理策略对该适应过程的影响。该分析框架有助于未来研究其他类型全球变化下不同管理策略中多种生物群落的适应响应与机制。研究建议，为通过物种组成变化促进群落组成层面的适应，应将全球变化强度维持在适应阈值以下，并合理实施确定性策略（如恢复土壤环境）与随机性策略（如添加适应物种种子）。