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 会导致模                           点网格被自然保护区所覆盖ꎬ主要位于四川中部、
            型的预测结果不可靠ꎬ我们仅对分布点大于或等                              云南西北部、滇黔桂交界处( 滇东南、桂西南和黔
   60   61   62   63   64   65   66   67   68   69   70