《Ecology》：Species-specific responses to environmental heterogeneity underpin spatial variation in annual plant communities
编辑推荐：
空间异质的环境可以通过物种对其直接环境的响应差异来促进物种共存。然而，物种对异质环境的特异性响应如何在群落组成促进共存的机制中贡献变异，目前仍不清楚。在本研究中，研究人员采用互置田间实验，考察了粗木质残体（CWM）分别提供的物理遮蔽和土壤调理两种细尺度空间变异
空间异质的环境可以通过物种对其直接环境的响应差异来促进物种共存。然而，物种对异质环境的特异性响应如何在群落组成促进共存的机制中贡献变异，目前仍不清楚。在本研究中，研究人员采用互置田间实验，考察了粗木质残体（CWM）分别提供的物理遮蔽和土壤调理两种细尺度空间变异来源，如何影响冬季一年生植物群落的组成以及群落中三种优势物种的发芽和生长表现。研究发现，CWM的土壤调理效应和物理遮蔽效应各自独立地促进了植物组成的空间变异。植物个体附近CWM的存在与否解释了群落组成约10%的变异，并对所有三种目标物种的发芽和生物量产生了显著的物种特异性效应。植物群落中更大比例的变异（约50%）由实验区块这一较粗尺度空间单元解释，表明较粗尺度的环境因子对物种组成的影响大于CWM带来的局部细尺度效应（约10%）。本研究提供了实证证据，表明物种对空间异质性的特异性响应与群落组成的空间变异相一致，符合空间共存机制运行的必要前提。细尺度空间异质性对整体物种组成的空间变异有贡献，但其效应相对小于较粗尺度的空间异质性。综合来看，研究结果表明植物群落组成是通过细尺度和较粗尺度空间过程的共同作用而构建的。
论文解读
一、研究背景与科学问题
理解环境空间变异与自然群落多样性维持之间的联系，是植物群落生态学的核心议题之一。从理论层面而言，空间异质的环境比均质环境能支持更多植物物种，因为不同环境需求的物种可以在不同地点成功或失败。现代共存理论（MCT）将空间异质性在稳定共存中的作用形式化，提出了空间相对非线性、空间储存效应和空间生长-密度协方差等空间机制，这些机制可以同时作用，防止低密度物种被竞争排斥，从而维持物种共存。尽管每种空间机制通过不同途径运作，但所有机制都要求物种对空间变异的环境产生不同响应。这些差异可导致物种在异质景观中的空间分离，进而引起群落组成的空间变异。
然而，目前物种性能对空间环境异质性的响应与植物群落组成的研究通常是独立进行的，物种层面的性能表现与群落层面的组成变异之间的联系在很大程度上仍停留在理论层面。缺乏将物种性能与群落组成联系起来的实证证据，意味着必须推断二者之间的关系，这可能导致高估或低估某一产生空间变异的特定过程对群落组成的影响及其在介导物种共存中的作用。因此，检验物种对空间异质性的特异性性能响应是否能转化为可检测的群落组成变异，至关重要。
西澳大利亚York gum-jam林地中物种丰富的一年生植物群落非常适合解决空间异质性、物种个体性能及其对群落组成影响之间缺失的联系。先前研究表明，该系统中的物种对不同环境梯度有不同的响应方式，观测研究也识别出群落组成对局部环境特征（主要是水分、遮荫、粗木质残体（CWM）和土壤磷）的响应。然而，物种层面和群落层面响应对局部空间变异的关系仍不明确。尽管有明确证据表明物种性能沿环境梯度存在特异性变异，但最近对空间储存效应的检验未能检测到该机制在自然群落中的强证据，这一令人意外的结果促使研究人员更全面地审视空间共存机制及其对群落结构构建的集体贡献。
本研究聚焦于CWM，因为它可以容易且一致地被操纵，以评估群落组成和个体物种的性能。CWM创造两种非互斥的环境变异类型：一是物理遮蔽（截留水分、提供遮荫、缓冲热和风），二是通过分解进行土壤调理作为养分来源。由于CWM在景观中斑块状分布，其物理和土壤相关效应可以独立地产生环境中持续的细尺度异质性。研究人员期望这些效应会不同地影响物种性能，特别是那些对热和水胁迫耐受性较差的物种，从而创造空间共存机制运行的必要条件。
二、研究设计与技术方法
研究人员在澳大利亚西南部小麦带Caron Dam自然保护区（29°34′14.25″ S，116°18′38.02″ E）的York gum-jam林地中开展了研究。研究系统气候为半干旱到地中海式，冬季（7-9月）生长季降水约110毫米/月，平均温度约20°C。2021年生长季总降水量约120毫米，2022年约190毫米。植被以稀疏的York gum（Eucalyptus loxophleba）和jam（Acacia acuminata）树冠层为主，林下为物种丰富的冬季一年生草本植物与散生的多年生草本和禾草的混合群落。
研究人员于2020年建立了七个实验区块（约25×25 m2），每个区块内设置32个1米线性样方，分配到六种处理之一：CWM（自然位置的粗木质残体，代表历史土壤预调理）、对照（无CWM区域，无土壤调理）、APS+ASC（将CWM移入无CWM区域，新增物理遮蔽+新增土壤调理）、APS（放置PVC管，仅物理遮蔽）、PSC（移除自然CWM，仅保留历史土壤预调理）、APS+PSC（移除CWM并替换为PVC管，物理遮蔽+历史土壤预调理）。该嵌套实验设计使研究人员能够区分CWM相关的细尺度空间变异与区块间的较粗尺度空间变异。
研究聚焦三种常见共存的冬季一年生植物：Trachymene ornata（Araliaceae，与CWM相关）、Trachymene cyanopetala（Araliaceae，生境灵活性较强）、Goodenia rosea（Goodeniaceae，开放全日照环境）。物种层面评估发芽率和地上生物量，发芽率以每个样方中出土幼苗数占播种种子数的比例衡量，地上生物量作为生长和繁殖潜力的指标。群落层面在盛花期（2020-2022年8-9月）沿每个样方的20厘米样线进行垂直标尺法调查，记录接触标尺的植物个体数量和身份。2022年9月采集土壤样品分析植物可利用磷、钾、有机碳、硝酸盐、铵和pH值。
统计分析使用R 4.3.2版本。物种层面使用广义线性混合模型（GLMM；glmmTMB包）分析发芽数据（采用hurdle框架处理零膨胀），使用线性混合效应模型（LMER；lme4包）分析生物量数据。群落层面使用冗余分析（RDA）和偏RDA检验CWM和区块对植物群落组成的影响，通过变异 partitioning 量化各因素解释的变异比例，使用置换检验（999次）评估显著性。
三、研究结果
物种特异性性能
发芽概率对CWM的响应在目标物种间存在差异，且跨年份定性一致。CWM对G. rosea的发芽概率无显著影响。T. ornata在土壤调理处理下的发芽概率显著高于G. rosea（2021年：emmean ± SE = 1.20 ± 0.57；事后两两比较：z = ?1.95，p = 0.03）。T. cyanopetala在土壤调理效应（2022年：emmean ± SE = 1.24 ± 0.61；z = ?2.04，p = 0.04）和物理遮蔽效应（2022年：emmean ± SE = 1.25 ± 0.61；z = 2.05，p = 0.04）下的发芽概率均显著高于G. rosea。发芽概率由固定处理效应（约9%方差）和随机区块效应（4%-6%）共同解释。条件于发芽的幼苗数量在物种间和处理间无显著差异。
单位生物量响应由“PS+SC”模型最佳解释。T. cyanopetala对土壤调理效应的生物量响应跨年份不一致，2021年边际正响应（emmean ± SE = 0.35 ± 0.27；t = 1.31，p = 0.19），2022年边际显著负响应（emmean ± SE = ?0.49 ± 0.28；t = ?1.74，p = 0.08）。T. ornata在2021年对土壤调理呈边际正响应（emmean ± SE = 0.39 ± 0.31；t = 1.25，p = 0.21），2022年达到统计显著（emmean ± SE = 0.69 ± 0.34；t = 2.06，p = 0.04）。因此，2022年T. ornata对土壤调理的生物量响应显著比T. cyanopetala更正向（emmean ± SE = 1.18 ± 0.44；t = ?2.69，p = 0.01）。2021年单位生物量由固定处理效应（24.6%）和区块随机效应（约5%）解释。
群落组成
CWM和区块在三年中均显著解释植物群落组成（2020年：调整R2 = 0.51，伪F = 2.93，p = 0.001；2021年：调整R2 = 0.26，伪F = 1.65，p = 0.001；2022年：调整R2 = 0.29，伪F = 1.74，p = 0.002）。与CWM相关的细尺度空间变异在三年中一致解释了一小部分组成变异（2020年：8.3%，伪F = 2.21，p = 0.01；2021年：12.1%，伪F = 2.12，p = 0.007；2022年：11.4%，伪F = 2.07，p = 0.01）。与实验区块相关的较粗尺度空间变异解释了更多物种组成变异，跨年份占53.7%-69%。
群落组成变异的驱动机制
利用2021和2022年的实验结果，偏RDA显示土壤调理效应在2021年解释了3.8%（伪F = 2.05，p = 0.003）和2022年解释了3.4%（伪F = 1.81，p = 0.13）的物种组成变异，而物理遮蔽效应分别解释了3.1%（伪F = 1.66，p = 0.03）和3.3%（伪F = 1.73，p = 0.02）。各年份物理遮蔽和土壤调理的加性效应小于Q2中CWM对群落组成的总影响（2021年12.1%，2022年11.4%）。
群落组成的大部分变异在较粗尺度的实验区块水平上观察到，与土壤养分和pH的空间更替相关。在每个区块内，CWM和开放群落组之间的组成差异通常与CWM样方中较高的土壤C和/或N相关。区块水平较高的pH和土壤P与群落组间更大的组成差异相关。
四、讨论与结论
本研究发现物种在发芽和生物量方面对细尺度和较粗尺度环境变异存在种间差异。细尺度上，这些差异与CWM产生的土壤调理和物理遮蔽效应相关的土壤和微气候变异有关。CWM解释了约10%的局部植物组成变异，而约50%的组成变异与跨研究地点较粗尺度的土壤空间变异（特别是养分和pH）相关。这些发现表明局部和较粗尺度过程同时介导局部植物物种多样性，多尺度空间共存机制可能在同一植物群落中运行。
物种特异性响应方面，目标物种在CWM功能效应下的发芽概率存在差异，Trachymene同属物种在CWM存在下相对于G. rosea具有发芽优势，这与先前观测中Trachymene物种在CWM附近更丰富的现象一致。Trachymene物种具有未成熟胚，由光和热提示完成发育，这种双重休眠发芽策略在不利年份提供保护，CWM的存在可能缓冲这种保守的发芽策略。与发芽的物种特异性差异不同，生物量响应未观察到跨物种的一致模式，表明CWM产生的空间变异主要影响发芽阶段而非成体生长阶段。年际气候变异可能影响CWM调理土壤中的物种性能，2022年较湿润年份观察到T. cyanopetala和T. ornata之间的生物量差异。
群落响应方面，CWM相关的细尺度空间异质性与植物群落组成的空间变化相对应，三年中一致解释了一定比例的组成变异。然而，较粗尺度的环境变异解释了更大比例的组成变异，表明局部植物群落由多尺度同时运作的空间变异塑造。这些群落层面模式与物种层面分析中观察到的物种对CWM不同功能效应的特异性响应一致，表明物种对细尺度环境异质性的差异响应有助于群落组成的空间构建。CWM的土壤调理和物理结构效应各自独立影响种间性能和局部群落的物种组成变异，每种效应都满足空间过程介导共存所需的生态条件。
较粗尺度空间变异的影响方面，群落组成受区块间较粗尺度空间变异的影响强于CWM相关的细尺度变异，可能反映多种非互斥机制：有限扩散导致区块间组成差异更大；较粗尺度的环境异质性，特别是土壤养分和pH的区块水平空间更替，也贡献于群落组成变异。较粗尺度的土壤肥力变异可能与CWM介导的局部共存机制相互作用，在背景土壤P和pH较高的区块中，CWM相关的局部土壤C和N增加倾向于在同一区块内产生更大的群落组成差异。
研究结论表明，CWM的物理遮蔽和土壤调理效应独立地促进共存物种间的生态位分化，塑造植物群落并可能促进空间共存。不同空间尺度的空间异质性共同影响群落组成，未来研究需要建立物种特异性过程与跨尺度空间结构之间的明确联系，方差分解方法为量化多种共存机制同时运作提供了有前景的途径。
