气候变化下中国濒危兰科植物独花兰（Changnienia amoena）与其专性传粉者分布区变化的对比研究

《Biology》：Contrasting Range Shifts of an Endangered Orchid Changnienia amoena and Its Obligate Pollinator Under Climate Change in China Yue Wang, Songwen Guo, Jingxin Zhou and Guangfu Zhang

　　

特殊的植物-传粉者相互作用极易受到气候变化的影响，然而预测其耦合响应仍具挑战性。濒危兰科植物独花兰（Changnienia amoena）及其唯一有效的传粉者——三带熊蜂（Bombus trifasciatus），是中国研究此类动态变化的典型系统。然而，受方法论差异的影响，现有关于两者潜在分布的研究存在显著分歧。研究人员应用Biomod2中的集合建模（Ensemble modeling）方法，整合多种算法，预测了两种物种当前及未来（2041–2100年，基于SSP1-2.6、SSP2-4.5和SSP5-8.5情景）的分布。模型基于全面的occurrence记录（C. amoena 123条，B. trifasciatus 43条）及关键环境变量构建。年降水量（Annual precipitation, Bio12）主要影响C. amoena的分布，而最干季度降水量（Precipitation of the driest quarter, Bio17）是B. trifasciatus的关键因子。气候变化下，C. amoena的适宜栖息地预计将略有收缩并向东南方向移动，而B. trifasciatus的气候适宜性可能增加。两物种的生态位重叠（Niche overlap）预计将减少，表明未来可能存在空间错配。本研究首次通过集合建模预测了该濒危兰科植物及其传粉者的当前适宜栖息地和未来变化。结果表明，这种特化的互利共生关系表现出相反的分布区转移，强调了气候驱动的相互作用中断风险。这些发现为制定针对C. amoena的保护策略提供了关键的科学依据，强调需考虑其专性传粉者的空间生态学。

研究背景与意义

在全球变暖及极端天气事件频发的背景下，物种分布、物候及适应性正经历显著改变，这对特有种及狭域分布的物种构成了严峻挑战。兰科（Orchidaceae）植物因其菌根特异性、特化的传粉系统及低种子萌发率等特性，对气候变化尤为脆弱。与此同时，传粉昆虫，特别是专性传粉者，也因气候波动面临生理耐受阈值被突破及生活史紊乱的风险。植物与传粉者之间若形成时空上的物候失配（Phenological mismatch），将严重降低授粉效率，进而威胁种群存续。尽管已有研究尝试将传粉者数据纳入物种分布模型（Species Distribution Models, SDMs），但多数依赖单一算法（如最大熵模型MaxEnt），且常忽略地形与人为因子，导致预测结果存在不确定性。独花兰（Changnienia amoena）作为中国特有的濒危兰科植物，被列为国家二级重点保护野生植物，其仅依赖三带熊蜂（Bombus trifasciatus）进行有效传粉。由于两者高度的专一性，该系统成为研究气候变化对特化种间互作影响的理想对象。然而，关于两者的潜在分布预测在中国尚存争议，且缺乏基于稳健算法的系统性评估。因此，本研究旨在通过集合建模（Ensemble modeling）解析气候变化下C. amoena与B. trifasciatus的潜在分布动态及生态位匹配度，为濒危兰科植物的协同保护提供理论支撑。该研究成果已发表于《Biology》杂志。

关键技术方法

研究人员首先通过野外调查、标本数据库（NSII、GBIF）及文献检索，收集了C. amoena和B. trifasciatus的分布记录，并利用SDMtoolbox进行空间稀疏化处理以消除采样偏差。环境变量涵盖气候（WorldClim）、地形（海拔、坡度、坡向）及人为（人类影响指数Human Influence Index）三大类共23个初始变量，经相关性分析筛选出9个关键变量用于建模。研究采用Biomod2平台构建集合模型，整合人工神经网络（ANN）、最大熵（MAXENT）、随机森林（RF）等10种算法，基于75%的训练集与25%的验证集进行模型训练与评估。模型性能通过受试者工作特征曲线下面积（AUC）和真实技巧统计值（TSS）进行评价，筛选出表现优异的模型进行集合预测。未来气候情景选用CMIP6框架下BCC-CSM2-MR模型输出的三种共享社会经济路径（SSP1-2.6、SSP2-4.5、SSP5-8.5）。利用ArcMap进行地理空间分析，并根据最大敏感度与特异度之和（MaxSSS）阈值划分适宜等级。此外，利用ENMTools计算Schoener's D和Hellinger's I指数以量化两物种的生态位重叠度。

研究结果

3.1 模型性能

利用Biomod2构建的集合模型显示，C. amoena的最终数据集包含123条记录和9个环境变量，B. trifasciatus包含43条记录和9个变量。经过严格的筛选标准（C. amoena: AUC > 0.9且TSS > 0.7；B. trifasciatus: AUC > 0.8且TSS > 0.7），排除了过拟合模型后，最终分别保留了7个和5个单模型用于集合建模。结果显示，C. amoena的集合模型AUC值为0.978，TSS值为0.885；B. trifasciatus的集合模型AUC值为0.958，TSS值为0.807。这表明集合模型优于单一模型，具有更高的预测准确性和可靠性。

3.2 环境变量的贡献

环境变量对C. amoena分布的影响中，气候变量占比最高（84.96%），其次为地形变量（11.92%）和人为变量（3.12%）。对于B. trifasciatus，气候变量同样占主导地位（71.76%），其次是地形（23.97%）和人为因素（4.27%）。具体而言，年降水量（Bio12，贡献率40.92%）是驱动C. amoena分布的关键因子，而最干季度降水量（Bio17，贡献率40.23%）则是B. trifasciatus的主要限制因子。这揭示了降水相关指标是调控这一专性互作系统空间格局的首要气候因子。

3.3 当前潜在适宜分布

在当前气候条件下，C. amoena的适宜栖息地面积（中、高适宜区）为108.06 × 10⁴km²，占中国国土面积的11.26%，其中高适宜区面积为39.25 × 10⁴km²。其分布主要集中在江苏南部、安徽南部、江西北部、浙江北部等华东地区，以及河南南部、湖南中部、湖北中部等华中地区和重庆、四川等西南地区。相比之下，B. trifasciatus的适宜栖息地面积（低、中、高适宜区）为232.02 × 10⁴km²，占中国国土面积的24.17%，高适宜区面积为15.95 × 10⁴km²。这表明在当前气候下，传粉者的分布范围远广于其寄主植物。

3.4 未来潜在适宜分布的变化

在未来气候情景下，两物种的分布呈现出截然不同的变化趋势。C. amoena的适宜栖息地面积在所有时期和所有情景下均呈现缩减趋势，且分布重心向东南方向迁移。相反，B. trifasciatus的适宜栖息地在所有情景下均表现出扩张趋势，其分布重心则向西部移动。这种反向的移动轨迹预示着两者之间将产生空间错位。

3.5 生态位重叠分析

通过计算Schoener's D和Hellinger's I指数发现，在当前气候条件下，两物种的生态位重叠处于中等水平。然而，在所有未来的气候情景（SSPs）中，这两个指数均呈现下降趋势。这意味着随着气候变化的加剧，C. amoena与其专性传粉者B. trifasciatus在空间上的重合度将越来越低，潜在的生态位分离将导致两者难以维持现有的共生关系。

讨论与结论

本研究通过高精度的集合建模揭示了气候变化对中国特有濒危兰科植物独花兰及其专性传粉者三带熊蜂的差异化影响。研究发现，尽管两者目前存在显著的生态关联，但在未来气候变暖的背景下，独花兰倾向于向东南退缩，而三带熊蜂则向西部扩张，这种反向的分布区转移将导致两者生态位重叠度显著降低。这种空间上的“解耦”（Decoupling）现象将直接导致授粉效率下降，从而进一步限制独花兰的种子结实率和种群恢复能力。这一发现挑战了传统保护生物学中仅关注单一物种保护的范式，强调了在气候变化背景下，对于依赖特化传粉系统的植物，必须采取“物种-传粉者”联合保护策略。研究结果不仅为C. amoena的未来就地保护（In-situ conservation）选址和迁地保护（Ex-situ conservation）规划提供了科学依据，也为全球其他面临类似风险的兰科植物及专性互利共生系统的保护管理提供了重要的参考范式。