全球变化背景下南方红豆杉地域分布变化.docxVIP

全球变化背景下南方红豆杉地域分布变化 摘要：气候是影响植物栖息地的重要因素之一，预测气候变化对植物潜在分布范围变动的影响，对促进植物资源的可持续利用具有重要意义。基于最大熵Maxent模型结合11个环境变量，预测2050s三种气候情景下（RCP2.6、RCP4.5和RCP8.5）南方红豆杉（Taxus chinensis var.mairei）在中国的潜在地理分布状况，分析影响其分布的主要因素，探讨其分布格局的改变对我国亚热带北界的指示意义。结果表明：（1）南方红豆杉的适宜栖息地（生境指数P0.2）主要分布在我国亚热带暖温带季风区，绝大部分核心栖息地（生境指数P0.6）分布在秦岭大巴山以南地区；（2）Jackknife测试结果显示，最冷季降水量（bio19）、气温平均日较差（bio2）、气温年变化范围（bio7）、最暖季平均温度（bio10）和海拔（Elev）对南方红豆杉空间分布适宜性影响最大；（3）随气候变化，2050s南方红豆杉有沿纬度向北和海拔向上迁移的趋势，并且我国亚热带北界受气候变化的影响将逐渐向北移动。 据政府间气候变化专门委员会（IPCC）第五次评估报告（AR5）预测，全球气候将持续变暖，到21世纪末全球平均表面温度可能比1986—2005年间升高0.3～4.8℃物种分布模型（Species Distribution Models,SDMs）已成为研究全球变化背景下物种与环境之间相关生态问题的一个极其重要的工具南方红豆杉是我国亚热带常绿阔叶林、常绿与落叶阔叶混交林中分布的特有种和重要的濒危孑遗树种之一，是国家Ⅰ级重点保护植物1 研究方法与数据来源1.1 数据来源和变量筛选南方红豆杉地理分布信息数据通过三种渠道收集：（1）从全球生物多样性数据库（GBIF,）下载；（2）从中国数字植物标本馆（CVH,）查找；（3）查阅其他相关研究文献选取三种环境变量类型，包括19个生物气候变量、1个地形变量和2个土壤变量。为确定不同气候条件下物种的分布，选用由全球气候模型（GCM）生成的当前（1960—1990年）和未来CCSM1.2 研究方法1.2.1 模型运行和精度评估在Arc GIS 10.2中将上述11个环境变量图层以中国行政矢量边界数据为基础，采用掩膜提取方法进行裁剪。将南方红豆杉的地理分布数据（csv格式）和11个环境变量数据（bil格式）分别导入Maxent 3.4.1版本中进行模型运算。参数设置：创建受试者工作特征曲线（ROC）、显示预测结果并采用Jackknife预测各环境变量的贡献率；设置全部样本点的25%用于评估模型，其余75%用于训练模型；最大背景值设置为10000，迭代运行500次；模型其余参数选择默认。利用ROC分析方法进行Maxent预测精度的检测，ROC取值范围为[0,1],ROC越接近于1表示模型预测结果越好。一般认为，ROC值0.5表示预测失败，0.5～0.7表示预测较差，0.7～0.9表示一般或中等，0.9～1.0表示非常好1.2.2 预测结果处理在Arc GIS 10.2采用重分类工具将Maxent模拟的南方红豆杉当前（1960—1990年）和未来（2041—2060年）栖息地可能分布的概率即生境指数（P）进行划分为：不适宜栖息地（P≤0.2）和适宜栖息地（P0.2）。适宜栖息地又分为三类，包括边缘栖息地（0.2～0.4）、中等栖息地（0.4～0.6）和核心栖息地（0.6）。为估计南方红豆杉的分布对气候变化的响应，对栅格数据重投影并通过栅格计算器来计算未来三种气候情景下南方红豆杉适宜栖息地分布范围的扩散或缩减情况，使用SDMtoolbox 2.4计算适宜栖息地分布质心随气候变化迁移的情况。为进一步分析南方红豆杉栖息地分布与海拔之间的关系，将其适宜栖息地（P0.2）按照中国DEM数据（来源于地理空间数据云，）划分为8个等级：200 m、200～500 m、500～800 m、800～1200 m、1200～1600 m、1600～2000 m、2000～2500 m和2500 m，计算每一等级的面积并进行制图分析，统一采用WGS＿1984＿Albers投影坐标系。2 结果分析2.1 模型预测精度评价Maxent预测的结果显示ROC训练集AUC=0.950，测试集AUC=0.901 (SD=0.017），表明Maxent对南方红豆杉地理分布范围预测效果良好，具有很高的精确度和可信度。Jackknife测试了11个环境变量对模型预测结果的贡献率。贡献较大的5个环境变量（贡献率90%）包括最冷季降水量（bio19）、气温平均日较差（bio2）、气温年变化范围（bio7）、最暖季平均温度（bio10）和海拔（Elev）（表1）。温度（贡献率=