气候变化背景下中国羊肚菌属物种分布区的最大熵模型预测研究学术报告
一、 研究团队与发表信息
本研究由来自中国西北大学（第一单位）和青海大学（第二单位）的联合团队完成。主要作者包括曹玉婷、卢昭平等，通讯作者为商倩楠和李忠虎。该研究以题为“Maximum Entropy Modeling the Distribution Area of Morchella Dill. ex Pers. Species in China under Changing Climate”的论文形式，于2022年7月8日发表于学术期刊《Biology》（2022年，第11卷，第1027页）。该期刊由MDPI出版，本论文遵循知识共享署名（CC BY）许可协议开放获取。
二、 学术背景与研究目的
本研究属于物种分布模型（Species Distribution Model， SDM）与保护生物学交叉领域，具体聚焦于气候变化生态学。研究对象是羊肚菌属（Morchella Dill. ex Pers.）真菌，这是一类具有重要食用和药用价值的珍稀子囊菌，在中国拥有丰富的野生资源。然而，过度采挖、生境破坏以及气候变化对其野生种群的生存构成了严重威胁。与动植物相比，关于气候变化对真菌地理分布影响的研究相对匮乏。因此，科学评估和预测气候变化对羊肚菌属物种分布的影响，对于其资源保护和可持续利用具有紧迫性。
本研究的主要目的是：1）识别影响中国野生羊肚菌属物种地理分布的关键环境因子；2）利用最大熵（MaxEnt）模型，模拟并预测羊肚菌在历史时期（末次盛冰期、中全新世）、当前以及未来（2050年代、2070年代）不同气候变化情景下的潜在地理分布变化；3）为未来气候变化背景下羊肚菌野生资源的保护和合理利用提供科学依据和理论参考。
三、 详细研究流程
本研究遵循了物种分布模型研究的标准流程，主要包括数据收集与处理、模型构建与验证、结果分析与解释三个核心阶段，具体步骤如下：
1. 物种分布数据收集与处理： 研究团队通过野外调查和文献检索，共收集到288个羊肚菌物种的分布点记录。为确保模型质量，避免空间自相关导致的过拟合，作者对数据进行了严格清洗：首先剔除重复记录，随后采用缓冲区方法，以环境变量2.5弧分（约5公里）的空间分辨率为基准，剔除了相互距离在5公里以内的重叠分布点。最终，保留了180个有效分布点用于后续建模。这些点主要集中在中国西南、西北等地区。
2. 环境变量获取与筛选： 研究共考虑了28个环境变量，涵盖气候、地形和土壤三大类。气候数据（19个生物气候变量，Bio1-Bio19）来自WorldClim数据库；地形数据（海拔、坡向）和土壤数据（7个变量，如土壤pH值、有机碳含量等）来自联合国粮农组织等的和谐世界土壤数据库。所有数据均统一空间分辨率为2.5弧分，并转换为ASCII格式以便在MaxEnt模型中使用。研究涵盖了五个时期：末次盛冰期（LGM）、中全新世（MH）、当前（Current）、2050年代和2070年代。未来气候数据采用了IPCC第五次评估报告中的CCSM4.0气候模型，并考虑了三种不同的温室气体排放情景（RCP2.6, RCP4.5, RCP8.5）。
为降低环境变量间的多重共线性对模型结果的干扰，作者进行了严谨的变量筛选。首先，利用MaxEnt模型初步运行，评估各环境变量的贡献率。其次，使用SPSS软件对全部环境变量进行皮尔逊相关性分析。结合贡献率排名和相关性系数（设定阈值为|0.8|），筛选出最终用于建模的18个环境变量。当两个变量高度相关时，仅保留贡献率较高的一个。这一步骤确保了输入模型的环境变量相对独立，提高了模型的稳定性和解释力。
3. MaxEnt模型构建与参数设置： 将筛选后的180个分布点数据和18个环境变量图层导入MaxEnt 3.4.1软件进行建模。模型参数设置如下：随机选取25%的分布点作为测试集用于验证模型精度，剩余75%作为训练集用于构建模型；采用子抽样法，进行10次重复运算以提高结果的稳健性；选择“最大训练灵敏度加特异性”作为阈值规则；输出格式为Cloglog。其他参数保持默认。模型性能通过受试者工作特征曲线下的面积（Area Under the Curve， AUC）值进行评估，AUC值越接近1表明模型预测性能越好。
4. 适生区划分与分布中心迁移分析： 模型运行后，得到羊肚菌在不同时期和情景下的存在概率图（0-1）。根据IPCC第五次评估报告的建议，研究者使用ArcGIS 10.2软件中的重分类工具，采用自然断点法将适生区划分为四个等级：非适生区（0 ≤ 概率 ≤ 0.13）、低适生区（0.13 < 概率 ≤ 0.35）、中适生区（0.35 < 概率 ≤ 0.63）和高适生区（0.63 < 概率 ≤ 1）。同时，利用SDM Tools工具箱，基于二值化的适生区分布图，计算了羊肚菌在不同时期分布区的质心地理坐标，从而直观展示其分布中心随时间和气候变化的迁移轨迹。
四、 主要研究结果
1. 模型精度评估： 所有时期（过去、当前、未来）MaxEnt模型的AUC值均较高。训练数据的平均AUC值在0.899至0.908之间，测试数据的平均AUC值在0.839至0.858之间。根据通用的评价标准（AUC>0.8为良好，>0.9为优秀），本研究所构建的模型预测结果可靠，具有较高的准确性。
2. 主导环境因子识别： 通过刀切法（Jackknife）分析各环境变量的贡献率发现，影响羊肚菌地理分布的最主要环境因子在不同时期具有一致性。贡献率排名前四的因子始终是：最干季度降水量（Bio17）、海拔（Elevation）、最冷季度平均温度（Bio11）和年均温（Bio1）。这四个因子的累计贡献率在不同时期均超过70%（范围在70.6%至80.0%之间），表明它们是决定羊肚菌分布的关键限制因子。其中，Bio17的贡献率最高（约31%-34%），凸显了水分条件，特别是干旱季节的降水，对羊肚菌生存的极端重要性。
单因子响应曲线进一步揭示了羊肚菌生存的适宜环境范围：最干季度降水量（Bio17）在11.32–77.78毫米之间，海拔在1480.78–3827.03米之间，最冷季度平均温度（Bio11）在-5.98–9.32°C之间，年均温（Bio1）在6.16–16.98°C之间。这些数值范围明确了羊肚菌偏好凉爽、湿润的中高海拔环境，是一种典型的喜低温、喜湿型真菌。
3. 潜在地理分布格局及其时空变化： * 当前分布格局： 模型模拟显示，当前中国羊肚菌的总适生区面积约为405.82万平方公里，占国土面积的42.34%。高适生区主要集中在云南、四川、甘肃、陕西、新疆等省份，这与已知的羊肚菌主要产区相符。 * 历史时期变化： 与当前相比，末次盛冰期（LGM）和中全新世（MH）的总适生区面积均有所减少，分别减少了12.43%和1.79%。特别是在LGM时期，适生区向中国西南部收缩，西北和东北地区的适生区显著减少。 * 未来时期变化： 在所有未来气候情景下（2050s和2070s的RCP2.6, RCP4.5, RCP8.5），羊肚菌的总适生区面积均呈现增加趋势。与当前相比，增加幅度在2.10%到6.71%之间。一个显著的趋势是，适生区表现出向高纬度、高海拔地区扩张的态势，具体表现为向中国东北部（黑龙江、吉林、内蒙古）和西北部（新疆、青海、甘肃）扩展。同时，在华中和华东部分地区（如湖南、湖北、江西、安徽等），适生区面积出现不同程度的退化。这种“北扩南缩”的分布格局变化，与许多物种为适应气候变暖而向高纬度/高海拔地区迁移的普遍规律一致。在不同排放情景下，扩张的程度和具体区域略有差异，例如在RCP8.5（高排放）情景下，2070年代的扩张面积最大。
4. 分布中心迁移轨迹： 分布质心计算结果表明，当前羊肚菌的分布中心位于陕西省西南部（106.67°E, 34.48°N）。从末次盛冰期到中全新世，其分布中心从重庆北部向西北方向迁移了约305.70公里至陕西西南部。在未来所有气候情景下，分布中心均表现出向东北方向迁移的趋势，预计将迁移至甘肃省甚至宁夏回族自治区境内。这进一步从几何中心移动的角度印证了羊肚菌适生区整体北移的结论。
五、 研究结论与意义
本研究得出结论：1）降水（特别是最干季度降水）、海拔和温度是影响中国羊肚菌属物种分布的最关键环境因子，定义了其喜冷凉、湿润、中高海拔的生态位。2）MaxEnt模型能很好地模拟羊肚菌在中国的潜在分布。3）在未来气候变暖背景下，羊肚菌的潜在适生区总面积将增加，但其空间格局将发生显著变化，表现为向中国东北和西北地区扩张，同时在部分中东部地区萎缩，分布中心向东北方向迁移。
本研究的科学价值在于：首次系统地利用物种分布模型（SDM）模拟了气候变化对中国重要经济真菌——羊肚菌属物种地理分布格局的历史变迁和未来趋势，填补了该领域的研究空白。其应用价值在于：为制定科学的羊肚菌野生资源保护策略提供了直接依据。研究指出，未来应重点关注东北和西北新扩张的潜在适生区，将其作为资源保护和人工抚育的重点区域；同时，对当前核心产区（如西南地区）需加强保护，以应对局部生境退化风险。这为在气候变化背景下实现生物多样性保护和资源可持续利用提供了重要参考。
六、 研究亮点
七、 其他有价值的内容
研究在讨论部分也坦诚地指出了本工作的局限性，这为后续研究指明了方向：1）模型未考虑物种扩散限制和人为传播的影响。羊肚菌的孢子传播能力及其与共生植物的专性关系可能限制其实际扩散速度，使其无法像模型预测的那样快速占据新的理论适生区。2）模型未纳入人类活动（如土地利用变化、采挖压力）和种间相互作用等难以量化的因子。3）将海拔作为独立变量，未深入探讨其与温度、降水等气候因子之间的时空耦合关系对物种分布的综合影响。作者建议，未来的研究可以尝试整合物种迁移过程到模型中，以更真实地模拟物种随环境变化的动态过程。这些思考体现了研究的严谨性和前瞻性。