Page 65 - 《广西植物》2023年第8期
P. 65
8 期 余江洪等: 国家重点保护野生植物的保护现状及潜在分布区预测分析 1 4 0 7
(http: / / www.mee.gov.cn) 和世界自然保护区数据 于 5 的物种进行模型运算(Pearsonꎬ 2006ꎻ Tang et
库(https: / / www.protectedplanet.net/ ) 的文件ꎬ绘制 al.ꎬ 2018)ꎬ772 种 国 家 重 点 保 护 野 生 植 物 用 于
了中国自然保护区图层ꎬ包括 464 个国家级自然 MaxEnt 模型的预测分析ꎮ
保护区和 806 个省级自然保 护 区ꎮ 利 用 ArcGIS 基于 MaxEnt 模型的预测分析ꎬ我们使用了两
v.10.6 软件将所确定的热点网格与自然保护区 种方法来构建模型ꎬ即对有 5 ~ 29 个分布点的物种
(国家级和省级) 图层叠加ꎬ评估当前保护网格对 采用刀切法(Pearsonꎬ 2006)ꎬ对大于等于 30 个分
国家重点保护野生植物的保护成效ꎬ并确定保护 布点的物种采用交叉验证法( Tang et al.ꎬ 2018)ꎮ
空缺ꎮ 当热点网格分布有自然保护区时ꎬ则认为 我们将 75%的分布数据设置为训练集ꎬ其余 25%
该网格 中 的 物 种 受 到 保 护ꎬ 否 则 即 为 保 护 空 缺 设置为测试集ꎬ其他设置设为默认值ꎮ 模型预测
(Hou et al.ꎬ 2010ꎻ Chi et al.ꎬ 2017ꎻ Yang et al.ꎬ 的结 果 采 用 受 试 者 工 作 曲 线 ( receiver operating
2021ꎻ Xue et al.ꎬ 2021)ꎮ 我们最终统计了热点网 characteristic curveꎬ ROC 曲线) 进行检验ꎬ评价标
格、保护成效和保护空缺中国家重点保护野生植 准为 ROC 曲 线 与 横 坐 标 所 围 面 积 ( AUC 值)
物的种类ꎮ (Kumar et al.ꎬ 2014)ꎬ AUC 值 为 0. 50 ~ 1ꎬ 其 中
1.4 气候数据获取及处理 0.70 ~ 1 表示模型性能好或极好 ( Swetsꎬ 1988)ꎮ
研究中所使用的环境变量包括气候因子和地 为了保证预测结果的准确性ꎬ仅选取 AUC 值大于
形因子(海拔)ꎮ 当前(1960—1990 年ꎬ version 1.4) 0.70 的物种进行后续分析ꎮ 利用 ArcGIS v.10.6 软
和未 来 ( 2070ꎬ CMIP5 ) 的 19 个 气 候 因 子 从 件提取分辨率为 10 min 的不同网格单元内的各个
WorldClim 数据库 ( https:/ / www. worldclim. org / ) 中 物种的适生值(0 ~ 1)ꎬ适生值越高代表物种在该
下载ꎮ 其中ꎬ19 个气候数据在 ArcGIS v.10.6 中使用 地区存在的概率越大ꎮ Yang 等(2013) 研究表明ꎬ
掩膜剪裁工具裁剪并提取到中国的气候数据ꎬ然后 适生值大于 0.40 为中度潜在分布区或高度潜在分
将其转换为 ASCⅡ格式的数据ꎮ 未来的气候数据选 布区ꎮ 因 此ꎬ 我 们 将 同 一 网 格 中 不 同 物 种 大 于
择 IPCC 第 5 次评估报告中 BCC ̄CSM 1 ̄1 模型下未 0.40 的适生值进行求和得到每个网格最终的适生
来四种气候变暖情景中代表 CO 最低和最高排放 值ꎬ以呈现最终的结果ꎮ
2
情景的 RCP2.6 和 RCP8.5ꎮ 海拔数据从 DIVA ̄GIS
2 结果与分析
(http:/ / www.diva ̄gis.org)数据库中获取ꎮ
在最大熵模型构建时ꎬ环境变量的多重共线
性可 能 会 导 致 模 型 过 度 拟 合 ( Grahamꎬ 2003ꎻ 2.1 物种丰富度及热点的分布格局
Pearson et al.ꎬ 2006)ꎬ从而影响拟合结果的准确 国家重点保护野生植物物种丰富度空间分布
性ꎮ 因此ꎬ本研究对构建模型的气候数据进行预 格局的研究表明ꎬ68.70%的网格分布有国家重点保
处理ꎬ利用 SPSS 13.0 软件对环境变量进行多重共 护野生植物ꎬ物种丰富度最高的网格主要集中分布
线性分析、检验变量之间的相关性ꎬ当 2 个变量相 在中国西南部和南部ꎬ包括四川中部、云南南部和
关性≥0.85 时ꎬ仅保留对预测概率贡献较大的变 东南部、广西北部、广东北部、海南等(图 1:A)ꎮ 根
量并用于模型运算ꎮ 经过相关性分析后ꎬ最终筛 据前 5%丰富度算法共确定了 200 个多样性热点网
选出 9 个环境变量:bio2(平均昼夜温差)、bio3( 等 格ꎬ这些热点网格主要分布在中国西南部、南部以
温性)、bio6( 最冷月最低温度)、bio7( 年温度变化 及东部地区ꎬ如四川中部、云南西北部、滇黔桂交界
差值)、bio8( 雨季平均温度)、bio12( 年降雨量)、 处(滇东南、桂西南和黔西南)、桂北黔南、鄂渝湘交
bio15(降雨量变异系数)、bio17(最干季降雨量) 和 界处、华东地区、广东北部和中部地区与海南等地
海拔ꎮ (图 1:B)ꎮ 本研究所确定的热点包括 81.98%(846
1.5 潜在分布区预测分析 种)的国家重点保护野生植物(图 1:B)ꎮ
本研究采用最大熵算法(MaxEnt v.3.4.1)对当 2.2 多样性热点的保护成效和空缺分析
前和未来气候条件下国家重点保护野生植物的潜 保护成效分析结果表明ꎬ171( 85.50%) 个热
在分布区进行建模ꎮ 物种分布点小于 5 会导致模 点网格被自然保护区所覆盖ꎬ主要位于四川中部、
型的预测结果不可靠ꎬ我们仅对分布点大于或等 云南西北部、滇黔桂交界处( 滇东南、桂西南和黔