基于Biomod2集成模型预测兰州百合生态分布格局及其关键环境驱动因子研究

本研究针对气候变化背景下特色药用植物兰州百合（Lilium davidii var. willmottiae）适宜生境不明的问题，通过整合物种分布模型（SDM）与多源环境数据，利用Biomod2集成建模框架精准预测其潜在分布区。研究发现温度季节性（bio4）、海拔（Elev）和最冷月最低温（bio6）是核心限制因子，首次揭示新疆等非传统分布区具有开发潜力，为珍稀药用植物资源保护与可持续利用提供了科学依据。

在全球气候变化加剧的背景下，植物物种的地理分布和生态适应性正面临严峻挑战。极端天气频发、水循环紊乱等现象不仅直接影响温度、降水和光照等关键气候因子，更通过改变物种的生长周期、繁殖策略和栖息地可用性，对生态系统稳定性产生深远影响。作为生态系统的初级生产者，植物的生存状况与碳循环、水循环等自然过程的稳定性密切相关。尤其对于生态位狭窄、扩散能力有限的物种，环境变化可能导致其适宜生境急剧缩减甚至局部灭绝。在这一背景下，如何准确预测环境变化对特定植物物种地理分布的影响，成为生物多样性保护领域的迫切课题。

兰州百合（Lilium davidii var. willmottiae）作为中国特有的"药食同源"珍稀植物，具有重要的文化和经济价值，被认定为国家地理标志产品。其鳞茎富含黄酮、生物碱、多糖等多种活性成分，具有免疫调节、降血糖血脂、抗氧化等功效。然而，高强度的人类活动和持续的环境变化使得兰州百合面临栖息地退化和遗传资源流失的风险。当前研究多集中于栽培技术、化学成分分析等领域，对其潜在适宜生境的预测研究尚属空白。为此，研究团队通过整合多源地理分布数据与环境变量，采用Biomod2集成建模框架，对兰州百合在中国的潜在适宜分布区进行精准模拟，并深入解析关键环境驱动因子。

本研究主要运用了以下关键技术方法：首先通过野外调查（2010-2019年）、GBIF数据库和文献检索获取物种分布记录，经空间细化处理得到120个有效分布点；其次从WorldClim数据库获取19个生物气候变量（bio1-bio19），结合地形因子（海拔、坡度、坡向）共22个初始环境变量，通过相关性分析和贡献率评估筛选出8个关键变量；最后采用Biomod2平台整合12种算法（包括RF、GAM、MaxEnt等），通过75%:25%的数据分割和伪缺失点生成策略构建集成模型（EM），采用AUC（0.855）和TSS（0.957）等指标评估模型性能。

环境变量筛选

通过Pearson相关性分析和预建模贡献率评估，从22个初始变量中筛选出8个关键环境变量：海拔（Elev）、坡度（Slope）、坡向（Aspect）、等温性（bio3）、温度季节性（bio4）、最冷月最低温（bio6）、最干月降水量（bio14）和降水季节性（bio15）。这些变量间相关性较低（|r|<0.8），有效避免了多重共线性问题。

模型性能与可靠性

随机森林（RF）模型表现最佳，广义可加模型（GAM）和广义提升模型（GBM）次之，而表面包络（SRE）模型精度最低（TSS=0.561）。集成模型（EM）综合9个高性能模型的结果，AUC达0.855，TSS达0.957，表明模型具有优异的预测能力。

环境变量贡献率与响应曲线

温度季节性（bio4）、海拔（Elev）和最冷月最低温（bio6）被识别为核心环境驱动因子。响应曲线显示：温度季节性的最适范围为519.65-937.41（标准差×100），峰值在725.89；海拔适宜区间为1586-2967米，最优值2340米；最冷月最低温耐受范围为-17.78°C至1.90°C，最适值约-4.49°C。

当前潜在适宜分布区

模型预测结果显示，高度适宜区（≥0.75）面积约35.79×10⁴km²，主要集中在甘肃兰州及周边地区，并向陕西、河南等省延伸；中度适宜区（0.5-0.75）面积23.18×10⁴km²；低度适宜区（0.25-0.5）面积39.47×10⁴km²。值得注意的是，新疆地区也被预测为潜在适宜分布区。

本研究通过Biomod2集成建模明确了兰州百合的生态适宜性主要受温度季节性、海拔和冬季极端低温的共同调控。其中，温度季节性（bio4）反映了物种对年温度波动幅度的适应性，适宜的中等波动范围（519.65-937.41）可能与其物候发育（如休眠打破、开花调控）密切相关；海拔偏好（1586-2967米）体现了该物种对山地温凉气候的依赖；而最冷月最低温（bio6）的耐受范围（-17.78°C至1.90°C）则揭示了其较强的抗寒能力与冬季低温需求。这些发现不仅与兰州百合的传统分布区高度吻合，还创新性地指出新疆等非传统产区可能具备开发潜力，为资源扩繁提供了新思路。

该研究的科学意义在于：首次系统量化了兰州百合的生态位特征，为珍稀药用植物的保护性开发提供了范式；建立的高精度分布模型可作为气候变化背景下物种适应性管理的动态工具；提出的新疆引种可能性需结合土壤特性、水资源约束等因子开展实地验证。未来研究应加强土壤因子（pH、养分等）的整合，开展气候变化情景下的动态预测，并通过田间试验验证新产区的生态适应性。本研究发表于《Scientific Reports》，为特色农业发展与生物资源保护提供了重要科学依据。