气候驱动下稀土元素超富集蕨类 Blechnopsis orientalis 在中国的潜在分布动态及其对植物采矿与保护策略的启示——基于优化MaxEnt与CMIP6多情景预测

《Environmental and Sustainability Indicators》：Climate-driven shifts in the distribution of the rare earth element hyperaccumulator fern Blechnopsis orientalis: Implications for phytomining and conservation strategies in China

　　

全球气温上升正在改变物种分布，但蕨类植物（pteridophytes）的长期动态仍缺乏记录。在此，研究人员预测 Blechnopsis orientalis 的潜在栖息地动态；该蕨类具有卓越的稀土元素（rare earth element, REE）富集能力，并在绿色植物采矿（phytomining）中具有应用前景。研究覆盖中国当代与未来气候情景。基于646条有效分布记录和10个综合筛选的环境预测因子，研究人员应用参数优化最大熵模型（Maximum Entropy, MaxEnt）（正则化乘数(regularization multiplier, RM)=1，特征组合(feature combination, FC)=QHPT，受试者工作特征曲线下面积(area under the receiver operating characteristic curve, AUC)=0.947）评价栖息地适宜性。未来趋势模拟选取四个代表性CMIP6（Coupled Model Intercomparison Project Phase 6，耦合模式比较计划第六阶段）情景（SSP126、SSP245、SSP370和SSP585，SSP即Shared Socioeconomic Pathways，共享社会经济路径），并分析三个时间片段（2050s、2070s和2090s）。主要发现：(1) Bio9（Mean Temperature of Driest Quarter，最干季均温）和Bio12（Annual Precipitation，年降水量）是该物种分布的主要驱动因子，累计贡献86.7%，证实温度和水分为核心限制因子。(2) 与多数稀有植物普遍范围收缩相反，B. orientalis 在所有气候情景下呈显著净扩张，最大预测扩张面积达59.64 × 10⁴ km²，而收缩面积低于0.83 × 10⁴ km²。(3) 适生栖息地质心总体向西北迁移，位移距离随情景严重性增加，在SSP585-2090s达108.75 km。这些发现表明B. orientalis对气候变暖具有强适应能力，其稀土超富集功能可能通过蒸腾依赖机制对其气候生态位施加额外约束。基于这些综合结果，研究人员提出三个系统保护与开发方案：(i) 优先将中国华南西南部相互连通的省际边界区作为野生种群就地保护核心功能区；(ii) 将中国华南东部内陆省际交界区发展为稀土植物采矿技术应用和推广的官方试点；(iii) 沿种群质心迁移轴建立专门种质保存基地。

研究背景与问题：全球变暖驱动植物向高纬高海拔迁移，但蕨类长期分布动态记录不足。蕨类孢子扩散强，可能比种子植物更易与气候平衡，但配子体无角质层保护，对水热胁迫敏感，响应机制可能不同。相比被子植物，蕨类分布预测不足，东亚尤甚。Blechnopsis orientalis 广布热带亚热带亚洲，中国南方常绿阔叶林常见。它是稀土元素（REE）超富集植物，组织REE浓度达1735–3850 μg g^?1，可形成纳米独居石，是绿色植物采矿（phytomining）候选。但研究多集中于生理生化药理，缺乏量化：气候如何控制其分布，未来气候变化如何重塑适生区。现有文献局限：无中国潜在分布建模；蕨类建模常缺参数优化，默认MaxEnt易过拟合；缺少CMIP6多情景多时段预测。研究人员提出三假设并开展研究，旨在揭示REE超富集植物气候变化响应，为保护、采矿选址和迁地保护提供依据。

研究方法：数据来自CVH、GBIF、CNKI和Google Scholar，原始2525点，经去重和2.5′×2.5′网格稀释后得646条有效记录（GBIF 535、CVH 89、CNKI 16、Google 6）。环境变量来自WorldClim 2.1、ArcGIS地形、MODIS NDVI和人类足迹指数（HFI），经Spearman相关和jackknife筛选10个。用ENMeval优化MaxEnt，RM 0.5–4.0，FC含L/Q/P/T/H，以delta AICc和OR10选择RM=1、FC=QHPT；75%训练、25%验证，10000背景点，10次重复，AUC/TSS/Kappa评价。未来用BCC-CSM2-MR和SSP126/245/370/585，2050s/2070s/2090s。阈值0.0994二值化，Natural Breaks分四类，分析稳定/收缩/扩张及质心。

研究结果：
3.1 模型优化与精度评估：默认RM=1、FC=LQHP的delta AICc=68.32，欠拟合；优化RM=1、FC=QHPT的delta AICc=0，OR10为0.2950，较默认0.2966降0.53%，AUC均为0.947。优化模型更简约，降低过拟合，增强泛化和生态可解释性。
3.2 分布决定因素：Bio9贡献52.6%，Bio12贡献34.1%，累计86.7%；气候因子总贡献97.0%。响应曲线显示Bio9 < ?6.10°C时概率极低，>约5.80°C时高；Bio12 <906 mm低，>约1027.23 mm高；最适条件为Bio9 ≥5.80°C且Bio12 ≥1027.23 mm。
3.3 当前适生分布：当前总适生134.24×10⁴ km²，占13.98%；高适生9.00×10⁴ km²（0.94%），主要在台湾、广东中南部、福建、海南、云南、西藏局部；中适生58.79×10⁴ km²；低适生66.45×10⁴ km²。
3.4 未来情景时空变化：四个SSP下总适生面积普遍增加。SSP126总适生166.62–174.22×10⁴ km²，高适生18.70–29.97×10⁴ km²；SSP245总适生179.64–193.19×10⁴ km²，高适生24.02–25.80×10⁴ km²；SSP370总适生170.81–193.35×10⁴ km²，高适生27.54–35.39×10⁴ km²；SSP585总适生182.07–194.20×10⁴ km²，高适生26.00–37.41×10⁴ km²。高适生向江西、湖南等扩展，SSP585-2090s集中于台湾、福建、浙江南部、江西、广东南部、广西东北、湖南南部、云南南部、西藏东南。
3.5 适生面积时间动态：所有情景净扩张。最大收缩0.83×10⁴ km²，最大扩张59.64×10⁴ km²；稳定率69%–82%，扩张率在SSP585-2090s达30.64%，收缩率≤0.49%。
3.6 质心迁移：当前质心位于广西桂林市全州县绍水镇（110°51′E, 25°56′N）。总体向西北迁移，SSP126有短暂西南波动。SSP585情景下最大单次位移108.75 km，后续继续向西北，2090s到达湖南洞口县魏家桥镇附近。

讨论总结：
4.1 模型优化：优化参数显著降低欠拟合和假阳性率，建议蕨类分布建模标准化参数优化。
4.2 温度降水主导：Bio9和Bio12合计贡献86.7%，最干季水热是关键。除蕨类基础生态需求外，REE超富集可能额外约束：独居石形成依赖蒸腾驱动REE吸收和磷酸盐沉淀，最干季低温干旱抑制气孔导度和蒸腾。阈值与华南离子吸附型REE矿区气候重叠，但权重不能量化，需温度控制与REE添加实验及非超富集蕨比较。
4.3 当前分布与地质REE背景：高适生仅0.94%且离散，与华南离子吸附型REE矿带空间重叠。HFI贡献1.3%，宏观非主导，但局地高HFI需关注破碎化。
4.4 未来净扩张与西北质心：多数稀有植物收缩，B. orientalis净扩张，稳定69%–82%。孢子扩散和生理可塑性可能贡献。适当浓度REE增强抗氧化酶，可能间接提高耐热，但未实验验证。净扩张不等于无风险，扩散限制可能使实际殖民低于潜在。质心迁移是气候适生几何位移，非实际种群迁移；成功殖民需微生境和土壤REE。
4.5 保护与植物采矿：就地保护优先稳定区，广西—广东—福建交界为核心。植物采矿试点可考虑扩张区新增高适生，如江西东部、湖南南部；先小规模试验在福建西部、江西南部、广东北部。气候适生不等于采矿适生，需土壤REE背景筛选；B. orientalis生物量低于Dicranopteris，可能适合独居石型REE回收。迁地种质沿质心迁移路径（广西全州县→湖南城步县→绥宁县→洞口县）布局。
4.6 局限与未来：全国高分辨率土壤REE数据缺乏，未叠加REE富集土壤；未纳入扩散限制；未纳入土壤理化性质；仅用单一GCM。未来应编土壤地球化学、便携XRF、数字土壤制图、扩散模型和移植实验、多GCM集合。

研究结论翻译：使用优化MaxEnt模型，研究人员首次量化了REE超富集蕨类B. orientalis在中国分布的气候影响。温度（Bio9）和水分（Bio12）是主导因素，合计贡献86.7%。该物种在所有CMIP6情景下净扩张（最大扩张面积59.64×10⁴ km²，最大收缩仅0.83×10⁴ km²），适生栖息地质心向西北迁移约108.75 km。Bio9主导可能反映蕨类最干季水热需求，REE吸收蒸腾依赖可能额外约束其气候生态位。方法上，优化模型（RM=1，FC=QHPT）较默认参数降低假阳性率0.53%，建立蕨类分布建模可重复优化范式。建议将广西—广东—福建交界作为优先就地保护区；在广东北部、福建西部、江西南部等地开展REE植物采矿预试验；沿泉州（福建省）—城步苗族自治县（湖南省）—绥宁县（湖南省）—洞口县（湖南省）轴线布置种质资源库。该研究未考虑繁殖体扩散限制和土壤有效REE空间异质性，扩张机制与REE超富集的因果关系需通过温度控制和REE添加交叉实验验证。