《Rangeland Ecology & Management》：Drivers of Grassland Biodiversity: Disentangling the Soil-Biome Interaction Across Physiographic Regions in Nepal’s Protected Areas
编辑推荐：
保护区内（protected areas, PAs）的草地生态系统是至关重要的生物多样性储存库，然而其植被结构与土壤性质在环境梯度上的相互作用仍未得到充分探索。本研究考察了尼泊尔沿地貌梯度的植物多样性、生物量和土壤特征，涵盖了喜马拉雅山脉、中山丘陵和低海拔特莱
保护区内（protected areas, PAs）的草地生态系统是至关重要的生物多样性储存库，然而其植被结构与土壤性质在环境梯度上的相互作用仍未得到充分探索。本研究考察了尼泊尔沿地貌梯度的植物多样性、生物量和土壤特征，涵盖了喜马拉雅山脉、中山丘陵和低海拔特莱（Terai）地区的保护区。特莱是沿尼泊尔-印度边境的低地带，由平坦冲积平原和沼泽草地组成。研究人员采用分层随机抽样方法，于2024年4月至2025年2月共采集288份土壤样品。植物区系组成出现显著分化：喜马拉雅地区表现出最高的物种丰富度和多样性。多变量分析确定了不同的植物群落，特莱以低地禾草为主，而喜马拉雅和中山丘陵以高海拔物种为主。Mantel检验显示，喜马拉雅地区土壤-植物多样性关系强烈，中山丘陵中度相关，而特莱地区相关微弱且不显著，表明后者生物地球化学与多样性动态存在解耦。土壤有机碳（soil organic carbon, SOC）储量在中山丘陵达到峰值，广义加性模型（generalized additive models, GAMs）确定土壤容重（bulk density, BD）为碳变化的主要预测因子。这些结果表明，草地结构和功能的生态驱动因素并非一致，而是随海拔显著变化，强调需要针对保护区内部本地化驱动因素制定特定区域的保护策略。
草地覆盖了地球上约40%的陆地表面，在全球碳循环、生物多样性维持和生态系统服务提供中发挥关键作用。草地土壤储存了超过三分之一的陆地碳库，但其完整功能受到土地利用变化、气候变迁和人为干扰的威胁。保护区（protected areas, PAs）作为生物多样性的庇护所，尤其对维持草地生态系统的完整性和功能至关重要。尼泊尔拥有从低地特莱（Terai）平原到高海拔喜马拉雅山脉的保护区网络，涵盖多样的草地类型。然而，以往研究多聚焦于单一生态系统或单点，缺乏对海拔梯度上土壤-植物关系系统比较；同时，草地生态系统功能生态学，尤其是生物量分配、碳动态及与土壤性质的机制联系，仍存在大量知识空白。为此，研究人员选取了尼泊尔三个不同地貌区的保护区——喜马拉雅山脉的Dhorpatan狩猎保护区（Dhorpatan Hunting Reserve, DHR）、中山地区的Shivapuri Nagarjun国家公园（Shivapuri Nagarjun National Park, SNNP）和低地特莱的奇旺国家公园（Chitwan National Park, CNP），旨在揭示草地群落结构、生物量分配、土壤性质与碳储量的变化规律，并识别驱动土壤有机碳（soil organic carbon, SOC）储量的关键因子。研究结果表明，植物多样性与群落组成呈现显著的空间异质性；土壤-植物耦合强度随海拔而变化，由喜马拉雅强耦合过渡至特莱的解耦状态；SOC储量在中山地区最高，且主要由土壤容重（bulk density, BD）预测。这些发现意味着草地结构和功能的生态驱动因素并非统一，而是随海拔急剧变化，从而强调需要针对保护区内局部驱动因素制定分区特异性的保护策略。该论文发表于《Rangeland Ecology》。
研究人员采用分层随机抽样方法，于2024年4月至2025年2月从三个保护区采集288个土壤样品（144个用于土壤性质分析，144个用于测定容重），并在每个1 m×1 m样方中收获植物以测定地上和地下生物量。实验室分析包括SOC（Walkley-Black湿氧化法）、全氮（total nitrogen, TN）、有效磷、有效钾、pH、电导率、机械组成和水分等。基于群落数据计算了Shannon多样性指数（H）、Simpson优势度指数（D）、Margalef丰富度指数（R）和Pielou均匀度指数（J）。统计方法包括层次聚类、方差分析（ANOVA）、主成分分析（PCA）、相关网络分析、Mantel检验和广义加性模型（generalized additive models, GAMs）。样本来自三个保护区的固定样地：每个保护区16个1-ha样方，每样方4个1-m2子样方。
**Plant diversity patterns and structure of community（植物多样性格局与群落结构）**
研究发现Shannon多样性指数（H）以DHR最高（1.90±0.10），其次为SNNP（1.86±0.11），CNP最低（1.78±0.11）。Simpson优势度（D）和Pielou均匀度（J）趋势类似，Margalef丰富度（R）同样以DHR最高（1.64±0.20）。重要值指数显示不同科的相对重要性存在显著差异。物种热图显示，CNP以低地禾草（如Phragmites karka、Saccharum spontaneum）为主，而DHR和SNNP以高海拔物种（如Meconopsis paniculata、Gentiana depressa）为主；7个共有种出现于所有保护区，18个物种为特定区域分布。
**Vertical stratification of biomass and soil physicochemical properties（生物量与土壤理化性质的垂直分层）**
生物量分配差异显著：DHR的地上生物量最高，而CNP的地下生物量和凋落物生物量最高。多数土壤性质随深度增加而降低，但BD随深度增加；各保护区之间所有土壤指标均存在极显著差异（P<0.001）。SOC储量以SNNP最高，CNP最低。
**Multivariate relationships among soil variables（土壤变量的多元关系）**
主成分分析揭示不同保护区土壤性质变异模式不同。DHR的PC1与SOC和BD密切相关；SNNP的PC1与SOC和BD呈负相关；CNP的PC1与砂粒含量、SOC、水分和黏粒负相关，PC2与pH、TN、有效磷和有效钾正相关，并与BD负相关。
**Plant-soil relationships: direct and multivariate linkages（植物-土壤关系：直接与多元联系）**
相关网络分析显示，DHR中SOC、TN、有效磷、有效钾和土壤质地构成高度正相关簇，并与BD负相关；SNNP中SOC-TN-海拔正相关簇与BD负相关；CNP中SOC-海拔-TN簇与黏粒相关养分正相关。Mantel检验表明，DHR中植物多样性与土壤因子呈强正相关（r≥0.4，P<0.05），SNNP呈中等显著相关（r=0.2–0.4，P<0.05），而CNP呈弱且不显著相关（r<0.2，P>0.05），说明低地保护区土壤与植物多样性之间存在生物地球化学解耦。
**Modeling the drivers of SOC stocks（SOC储量的驱动因子建模）**
广义加性模型显示，在SNNP中BD是SOC储量的首要预测因子（P<0.001，R2adj=0.528），解释超过50%的碳变异；在DHR和CNP中，TN是显著预测因子（P<0.05），而生物量和有效磷在所有保护区均无显著作用。
讨论部分总结了土壤-植物耦合强度沿海拔的变化及其机制。DHR受低温和养分限制，植物群落结构由土壤因子（SOC、TN、粉粒）强烈驱动；SNNP处于中等耦合状态，黏粒含量对植物-土壤关系具有调节作用；CNP因放牧和踩踏等自上而下控制，导致植被与土壤性质解耦。SOC储量的驱动因子随区域而异，BD在SNNP中的强预测力反映了土壤团聚体对碳的物理保护，而历史土地利用遗留效应可能削弱了生物量与当前碳库的关联。研究还指出，单一时间快照和缺乏高分辨率气候数据是局限性，但对完整生态系统的原位测量为保护提供了关键基线。
研究结论（翻译自原文Implications部分）：这项研究强调了保护区作为生物多样性热点和生态避难所的关键作用，捕捉了尼泊尔主要地貌区植被多样性、生物量与土壤性质之间的复杂相互作用。结果表明，每个保护区都拥有尼泊尔草地生物多样性的独特组成部分。科学上，这挑战了沿海拔梯度具有统一草地群落组装规则的假设，证明地貌异质性在保护区内产生了不同的生态策略。植物多样性与土壤性质（尤其是SOC和BD）之间强烈且位点特定的相关性突出表明，土壤完整性是保护区内生态系统结构的根本驱动因子。喜马拉雅和中山保护区紧密的土壤-植物关系表明，其保护至关重要，既能保护巨大生物多样性，也能保护关键的地下生物量。相比之下，奇旺国家公园（CNP）的解耦关系意味着其管理必须优先考虑水文和干扰等外部因素以维护生态价值。关键启示是，保护政策必须采用地貌特异性的策略，摒弃一刀切方法。例如，DHR和SNNP需要侧重于保护土壤完整性，而CNP则需要管理人为压力和水文条件。此外，确定BD为SOC的关键预测因子为快速碳储量评估提供了实用工具，增强保护区在气候减缓策略中的作用。未来研究应优先在保护区内开展长期监测和操控实验，以追踪生态对气候变化的响应，并扩展包括动物多样性和传粉网络的研究，最终将成果转化为与当地社区合作的管理策略，以确保尼泊尔保护区草地的长期韧性。
