《Diversity》:Mountaintop-Extinction Risk for Polylepis tarapacana: An Elevational-Ceiling and No-Escape Framework in the South-Central Andes
编辑推荐:
持续的变暖正在使山地物种向高海拔迁移,然而已经处于山体上限的树木面临着有限的垂直避难所。本研究评估了世界最高海拔树种Polylepis tarapacana在其横跨中南安第斯山脉(秘鲁、玻利维亚、智利和阿根廷)的分布区内,随着气候变暖是否能维持在气候适宜的地形
持续的变暖正在使山地物种向高海拔迁移,然而已经处于山体上限的树木面临着有限的垂直避难所。本研究评估了世界最高海拔树种Polylepis tarapacana在其横跨中南安第斯山脉(秘鲁、玻利维亚、智利和阿根廷)的分布区内,随着气候变暖是否能维持在气候适宜的地形中,并在秘鲁南部的塔克纳(Tacna)进行了野外验证。目的是在延伸至2100年的CMIP6(Coupled Model Intercomparison Project Phase 6)情景下,量化该物种保留多少适宜栖息地,以及任何损失中有多少比例离开了可达的更高海拔避难所。研究人员采用基于GBIF(Global Biodiversity Information Facility)出现记录和野外样地、经共线性过滤的生物气候预测因子及采样偏差校正背景校准的MaxEnt(Maximum Entropy)模型(以随机森林(Random Forest, RF)为基准),在CMIP6二十一世纪情景下作为多全球气候模式(Global Climate Model, GCM)集合进行投影,并针对148个野外制图林分多边形(12,874 ha)进行外部验证。一种新颖的海拔上限和“无逃逸(no-escape)”框架量化了每种情景下无任何可达更高海拔避难所的适宜面积损失。模型对适宜栖息地的判别良好(外部AUC,Boyce)。在最高排放路径下,塔克纳的适宜栖息地从(2030s)收缩至(2090s),高适宜等级崩溃;在各部分损失情景下,适宜性中心向高海拔迁移达 m。到2081–2100年,在所有测试的传播缓冲区内无逃逸比例达到,这是迄今为止关于世界最高海拔树木山顶灭绝风险的首次定量证据。这些结果主张在该垂直逃逸路线关闭之前,采取预期的迁地(ex situ)保护以及对最高海拔阿尔蒂普拉诺(Altiplano)避难所的秘鲁–玻利维亚跨境保护。
研究背景与问题提出
持续变暖导致山地物种向高纬度和高海拔迁移,这是气候变暖最具一致性的生物信号之一。对于已受限于宿主山体范围最高海拔界限的生物而言,这种向上追踪受到上方有限陆地的约束:随着等温线上升,气候适宜地形面积向顶峰收缩,对于最高栖息的分类群可能完全消失,即“灭绝电梯(escalator to extinction)”效应。高海拔山地物种因此成为气候变化的标志性哨兵,而威胁的严重性取决于物种为保持在适宜气候内需做出的向上偏移与其地形物理能提供向上偏移之间的平衡,这一点很少被明确量化。木本物种在树线上限面临全球性的海拔天花板约束。Polylepis属(Rosaceae)形成地球上最高海拔的森林,介于约3000至5000米之间,维持水文调节、土壤保持、碳储存及特有的特有生物区系。其中Polylepis tarapacana尤为特殊,在秘鲁和玻利维亚的阿尔蒂普拉诺高原可达约5200米,是已知最高海拔的被子植物树种。其局限于孤立有限山脉的狭窄海拔带,使其成为几乎无垂直空间逃离变暖气候的典型案例,同时持续的破碎化和土地利用压力已在侵蚀其种群,且高海拔Polylepis森林对气候极为敏感。相关性物种分布模型(Species Distribution Models, SDMs)已成为预测气候变化重分布适宜栖息地的标准工具,最大熵模型(MaxEnt)和随机森林(RF)分类器是应用最广泛的算法。然而对于Polylepis tarapacana,分布模型仍然稀缺且通常基于范围的国别子集校准;截断范围内的校准会偏倚估计的生态位及衍生投影;且二元验证无法揭示向高海拔迁移的种群是否能在上方找到适宜地面或其能否到达。当前变暖远快于地质过去的气候偏移,而长寿树木迁移缓慢,等温线上升速度约比树线推进快一个数量级,这种错配导致从个体到种群乃至物种水平的损失。山顶灭绝的补充概念——高海拔物种因其气候空间升至可用地形顶部之上而丧失——主要针对热带山地鸟类等动物阐述,很少在树木中操作化,几乎从未用独立野外数据验证。因此,世界最高树木是否需要向上偏移及其地形能否供应,这一核心问题仍未得到解答。该研究旨在填补此空白,以Polylepis tarapacana为对象,操作化海拔上限与无逃逸框架,并结合野外实证评估其山顶灭绝风险。该研究论文发表在《Diversity》。
主要关键技术方法
研究人员采用可重复的单配置工作流程耦合Polylepis tarapacana的全范围物种分布建模与新海拔上限分析。发生数据与环境数据来源于GBIF(Global Biodiversity Information Facility)的四国范围记录,经CoordinateCleaner过滤及物种特定海拔过滤(3500–5300米),并加入秘鲁塔克纳地区的16个野外样地,经spThin空间稀疏至1公里最小距离得60个空间独立存在位点用于校准;背景点采用目标群组背景校正采样偏差。环境预测因子包括WorldClim 2.1的19个生物气候变量、SoilGrids 2.0的5个土壤层及SRTM衍生的地形变量,经Pearson滤波(|r|)与方差膨胀因子(Variance Inflation Factor, VIF)逐步选择(<5)保留14个预测因子。未来气候数据取自WorldClim CMIP6档案,涵盖5个全球气候模式(GCM:ACCESS-CM2, EC-Earth3-Veg, GFDL-ESM4, MPI-ESM1-2-HR, MRI-ESM2-0)、4个共享社会经济路径(Shared Socioeconomic Pathways, SSP:126, 245, 370, 585)及4个时期(2021–2040/2041–2060/2061–2080/2081–2100)共80个情景。建模采用MaxEnt(maxnet通过ENMeval调参,特征类Hinge,正则乘数3)与RF(1000棵树,多数类下采样)算法,并进行空间分块五折交叉验证(blockCV,块大小约140公里)。通过多GCM集合均值与标准差构建不确定性,采用多元环境相似性表面(Multivariate Environmental Similararity Surfaces, MESS)识别新气候并钳制投影。新颖的后处理模块量化无逃逸(no-escape)面积:基于当前与未来二元适宜性栅格及重采样至适宜性栅格的SRTM海拔,在平方移动窗口内计算无更高(或等)海拔未来适宜细胞存在时的无逃逸细胞,评估1、10、25、50公里梯度传播缓冲区。野外验证采用塔克纳148个野外制图林分多边形(12,874公顷)作为存在及Polylepis探查的不存在掩膜,并在16个清查样地记录生长形态(单干与多干)作为气候胁迫的形态学代理。
研究结果
3.1. Field-Validated Current Distribution(野外验证的当前分布)
在空间分块交叉验证下,MaxEnt模型能较好区分适宜与不适宜条件(AUC,TSS,Boyce),RF表现相似(AUC,TSS,Boyce)。针对野外制图林分多边形的外部验证显示MaxEnt具有高敏感性()、中等特异性()、TSS、AUC及非常高的外部Boyce指数(),证实模型可靠排序栖息地质量。当前适宜性表面显示范围级及塔克纳裁剪后的分布。塔克纳16个清查样地记录457个个体,81%为多干,这种以矮曲林(krummholz)为主的形态是典型的生理耐受上限边缘树木特征,独立定性证实了上限处的气候胁迫。
3.2. Projected Habitat Shifts(预测的栖息地变化)
塔克纳Polylepis tarapacana的气候适宜栖息地在四个CMIP6时期随排放严重性上升而收缩。相对于当前适宜面积3306平方公里(适宜性中心4385米),SSP126下轨迹温和稳定(至),SSP245下中等恶化(2030s为,2081–2100为),SSP370下严重(2030s为,2070s为,2081–2100为),SSP585下灾难性:2030s丧失,中叶,2061–2080骤降至,2081–2100降至——在塔克纳maxSSS阈值内无适宜栖息地。伴随面积下降的是适宜性中心一致向高海拔位移——SSP126为米,SSP245为米,SSP370为米(2081–2100),SSP585在2070s已达米;SSP585世纪末海拔挤压完成且中心不再定义。高度适宜等级(≥0.5)从3551骤降至20平方公里(SSP585 2081–2100)。面积损失与向上中心位移证实假设(i):最高排放路径下地形天花板驱动适宜面积的净损失而非单纯重分布。
3.3. Mountaintop-Extinction: No-Escape Fraction(山顶灭绝:无逃逸比例)
在当前适宜面积中,SSP585下2081–2100所有四个敏感缓冲区(1、10、25、50公里)对应为无逃逸细胞——因该情景下塔克纳适宜栖息地完全丧失的算术结果。无逃逸信号在前一十年已显著(SSP585 2061–2080:1公里缓冲,10公里,25公里,50公里),标志2061–2080/2081–2100过渡为灾难性步骤。中间SSP245世纪末变暖比例适中(1公里,10公里,25公里,50公里);SSP370世纪末高(1公里缓冲,10公里,25公里,50公里);SSP126各缓冲与各时期基本可忽略(≤6.2%)。因此Polylepis tarapacana的山顶灭绝风险在高排放下世纪末为总量且缓冲不变,该情景下2070年代为十年驱动,低至中间强迫下适中——证实假设(ii)无逃逸比例随排放严重性急剧上升。无逃逸细胞空间分类显示中间SSP245世纪末变暖高海拔核心存续具可达上坡避难所,而SSP585下每当前适宜细胞变为无逃逸。
讨论部分总结
讨论部分指出研究回答了开篇三个问题:预测栖息地损失严重并按排放严重性单调有序,高排放下塔克纳适宜面积完全丧失;无逃逸比例从SSP126可忽略升至SSP585整个当前范围,一旦无区域内避难所则与缓冲无关;81%多干矮曲林形态证实现实种群已在气候边缘。讨论将“灭绝电梯”从主要记录移动山地动物类群操作化为首个针对树种的细胞水平山顶灭绝风险度量(无逃逸比例),使K?rner海拔天花板概念在细胞水平定量化。对于塔克纳Polylepis tarapacana,SSP245下2081–2100适度比例无上坡避难所(1公里,10公里),SSP585升至整个适宜面积且无区域内上坡避难所;伴随适宜中心上移及树线推进慢于等温线上升一个量级,支持近零传播缓冲为更现实参考。结果与广泛分布建模文献一致而非异常,近完全高排放损失反映物种极端海拔局限而非模型悲观。形态学多干主导提供独立定性支持种群处于耐受上限,但非气候唯一诊断(可能含干扰)。方法论上10公里缓冲为乐观上限,1公里为审慎参考;网格平滑可能低估局部微气候避难所,无逃逸比例最好读作气候暴露度量而非人口灭绝预测;残留不确定性来自土壤地形恒定及GCM间变异性;存在点样本适中(60)限制四国生态位解析度但外部Boyce高确认排序可靠。保护含义:高排放下2081–2100无区域内避难所,优先转向迁地(ex situ)保护(种子库、种源园试验)及与玻利维亚跨境协调(其核心范围留存更久适宜条件);中间SSP245最高海拔部分存续为保护及辅助迁移试验(如Yucamane火山北坡及>4500米边境山脊)及缓解土地利用压力。Polylepis tarapacana成为中南安第斯保护及国际(秘鲁–玻利维亚–智利–阿根廷)高安第斯树线生态系统协调的旗舰种。
结论部分翻译
山顶灭绝主要用于研究足以跟随气候上山的动物。本工作表明相同威胁可被测量,并基于野外数据立足于树木。地球最高海拔树种Polylepis tarapacana的栖息适宜性在中南安第斯全范围建模,那些全范围投影随后用于评估其在秘鲁南部塔克纳地区(该物种秘鲁据点)气候驱动的山顶灭绝风险。虽许多树种有气候驱动的范围投影,研究人员未见针对Polylepis tarapacana可比的细胞水平、野外验证的无逃逸评估。研究贡献三要素:一无逃逸度量区分未来栖息地损失中上坡仍有适宜地面可达与无剩余的情况,使海拔天花板可逐细胞测量;二在全物种分布拟合的MaxEnt和RF模型并以塔克纳独立野外制图林分检验;三形态证据,大多数清查树木为低矮多干矮曲林,即耐受上限边缘种群的习性。三个引导问题指向同一方向。各情景下适宜栖息地均预期收缩并上移;最强排放下世纪结束前在塔克纳几近消失,该损失是真实而非移至更高地面,因无更高地面可占。无可达上方避难所的栖息地比例随排放上升快速增加,直至最强变暖覆盖整个当前范围且不再依赖物种能传播多远;最温和路径下保持微小。矮曲林形态增加独立迹象表明种群已处气候极限。随之两慎示:此处识别的约束是几何的,由地形及气候走向设定而非单一模型运行,难以否定;即便如此,它衡量对气候的暴露而非灭绝本身;该暴露是否变为衰退取决于传播、招募与存活,本研究未解决。在此限度内,Polylepis tarapacana作为其他范围受限山地物种的哨兵。当 Summit-bound 植物的气候空间从其生长土地顶部抬起,无传播能追随。绘制何处逃逸仍可能与何处不可,指向中南安第斯值得保护的避难所,并提供可应用于其他海拔受限物种(含树线树木)的方法。实践中指向双轨策略:低排放路径下保护存续的高海拔核心并尝试辅助迁移至最高可用地面;高排放路径下无区域内避难所剩余时,优先跨种源迁地种子库及秘鲁–玻利维亚协调以保障留存最久适宜气候的阿尔蒂普拉诺种群。
需要我帮你把这篇论文解读中涉及的专业术语与英文缩写整理成一份对照表方便你查阅吗?