《Discover Life》：Climatic niche stasis of Pinus rzedowskii (Pinaceae), a relictual and endangered species
编辑推荐：
气候变化数千年来塑造了物种分布，物种可能向更适宜栖息地迁移，或持续存留于气候稳定的避难所中。本研究旨在识别历史与未来气候变化情景下墨西哥稀有濒危物种松属热泽多斯基松（Pinus rzedowskii）生态位的变化，以甄别其持续存活的适宜位点。研究人员利用三种算
气候变化数千年来塑造了物种分布，物种可能向更适宜栖息地迁移，或持续存留于气候稳定的避难所中。本研究旨在识别历史与未来气候变化情景下墨西哥稀有濒危物种松属热泽多斯基松（Pinus rzedowskii）生态位的变化，以甄别其持续存活的适宜位点。研究人员利用三种算法，基于一个全球环流模型（Global Circulation Model）与两条对比的共享社会经济路径（SSP1-2.6与SSP5-8.5），分别构建现生、末次盛冰期（LGM）及2071–2100年未来时段的生态位共识模型（Ecological Niche Consensus Models）。基于上述模型估算生态位质心与生态位重叠度（Schoener’s D）。研究发现降水季节性（Precipitation Seasonality）对该物种生态位贡献最大，其后依次为干旱指数（Aridity Index）与温度季节性（Temperature Seasonality）。适宜区面积呈从LGM到现生下降、再至未来略增的趋势；LGM生境质量最低，2071–2100时段最高。现生与未来生态位空间构型相似，但与LGM不同；然而质心与重叠度指示生态位随时间高度稳定。研究人员结论为：南马德雷山脉（Sierra Madre del Sur）尤其米却肯亚区（Michoacán subprovince）经历气候停滞，充当松属热泽多斯基松的气候潜在避难所，其分布中心 temporal variation 最小。
研究背景方面，气候作为影响物种分布的关键生态因子，其波动可引发物种纬向与海拔向位移导致种群瓶颈，或形成气候稳定避难所使物种长期留存。松类（Conifers）演化历史长、适应性强，但墨西哥多数高山松种分布狭窄且濒危，气候变化可能进一步压缩其分布并加剧灭绝风险。已有生态位模型（Ecological Niche Models, ENMs）多关注广布种或破碎种群的现生与未来情景，对过去气候重建及窄域濒危松种（如松属热泽多斯基松 Pinus rzedowskii）跨尺度评估不足。该物种为墨西哥米却肯州南马德雷山脉煤尔科曼山（Sierra de Coalcomán）特有残遗濒危松，仅存17个破碎种群、生殖力极低而遗传变异较高，是检验“气候变迁导致栖息地丧失与局地灭绝”假说的理想系统。研究人员由此提出假说：因其当前分布受限破碎，未来气候变化将缩减其潜在分布，并借ENM回溯LGM、对照现生与2071–2100年两SSP路径，甄别决定变量、历史趋势与最适存续区，为优先保护提供基线。本文发表于《Discover Life》。
关键技术方法上，研究人员整合文献与野外调查得331条松属热泽多斯基松出现记录，按1 km2像元去重保留27个独立点；环境层取30弧秒（~1 km2）分辨率CHELSA现生与未来（2071–2100）、PaleoClim的LGM生物气候变量，剔除BIO19并加入德马东干旱指数（De Martonne Aridity Index, AI = P/(t+10)），经Pearson相关（r≤0.8）与方差膨胀因子（VIF>10）筛选得温度季节性（BIO04）、最冷月均最低气温（BIO06）、降水季节性（BIO15）及AI四变量；校准区限定南马德雷米却肯亚区。采用biomod2包GLM、GAM、MAXNET三算法共识集成（EMwmean，仅纳TSS≥0.7单模），5折交叉验证10重复、1:10伪缺失配比，用ROC、TSS、Kappa评估；历史与未来直接转移现生参数并加clamping mask；LGM用NCAR CCSM4、未来用MPI-ESM配SSP1-2.6与SSP5-8.5；连续概率按最小训练存在阈值（Minimum Training Presence, MTP）二值化，依P50/P75分四级适宜度；借适宜度加权质心与Schoener’s D量化重叠，Haversine算质心位移。
研究结果部分，变量贡献与模型验证小节指出，BIO15贡献41.02最高，AI 25.63次之，BIO04为22.84，BIO06为20.33；集成模型EMwmean表现最佳（TSS=0.950，ROC=0.990，Kappa=0.724）。生态位地理投影小节显示，现生适宜区占南马德雷山脉7.2%（81000 ha），高适宜占2.7%集中于煤尔科曼山中部；LGM适宜面积比最大（19.7%）但全为低适宜；未来F-1.26适宜12.2%（高适宜5.4%），F-5.85为9.5%（低适宜为主），均值F-AV 12.2%但低适宜占6.6%；空间上现生与未来几乎重合，LGM偏南。生态位时间动态小节表明，适宜面积LGM→现生降约12.5%，现生→未来增约5%；生境质量（高+中适宜之和）LGM最低、未来最高；质心LGM→现生西北移5.4 km，现生→未来南移0.5 km；Schoener’s D：LGM-现生0.869、现生-未来0.936、LGM-未来0.867，均>0.86，证明生态位高度停滞。
讨论与结论翻译部分，研究人员指出降水季节性（BIO15）、干旱指数与温度季节性长期界定该种生态位，与近缘残遗松种（如Pinus pinceana、Pinus maximartinezii）及多项墨西哥松杉研究一致，湿度/干旱是松属分化驱动。空间上生态位构型跨时段稳定，LGM广而低质、现生缩而高质、未来微扩，质心位移<5.4 km，重叠度>0.86，符合气候停滞（climatic stasis）而非迁移假说；南马德雷米却肯煤尔科曼山因白垩纪古老地貌与间冰期缓冲充当潜在气候避难所（climatic potential refuge），分布中心变率最小。研究人员结论：松属热泽多斯基松生态位在LGM、现生及2071–2100双SSP下均表现停滞，南马德雷尤其米却肯亚区为气候潜在避难所，其中心分布区 temporal variation 最低，保护应优先锁定该高稳定核心区。
