Page 52 - 《广西植物》2024年第11期

２ ０ ３ ６ 广 西 植 物 ４４ 卷
１.３ 自然保护区分布和植物园引种栽培记录 量>２５０ 的网格作为样本相对较为完善的取样单元ꎮ
根据 生 态 环 境 部 发 布 的 自 然 保 护 区 名 录 以上分析使用 Ｒ 软件包 ＫｎｏｗＢＲ(Ｌｏｂｏ ｅｔ ａｌ.ꎬ ２０１８)
(ｈｔｔｐ:/ / ｗｗｗ.ｍｅｅ.ｇｏｖ.ｃｎ) 和«江西省生物多样性保 来完 成ꎮ 从 世 界 气 候 数 据 集 ２. １ 版 ( Ｆｉｃｋ ＆
护战略与行动计划(２０１３—２０３０ 年)»列举的江西省 Ｈｉｊｍａｎｓꎬ ２０１７)获取年平均温度、年平均降水量等
自然保护地名录绘制江西省 １５ 个国家级自然保护 １６ 个环境变量ꎮ 通过相关分析ꎬ以皮尔逊(Ｐｅｒｓｏｎ)
区和 ３２ 个省级自然保护区的矢量图层ꎮ 将标本分 相关系数<０.７ 为标准ꎬ最终保留年平均温度、温度
布点与自然保护区图层叠加ꎬ统计自然保护区有分 的季节性变化、最冷月的最低温、温度的年际变化
布的重点保护野生植物种类数量ꎮ 根据«中国迁地 范围、年平均降水量、最干月的降水量、降水的季节
栽培植物志名录» (黄宏文ꎬ２０１４)、«中国迁地栽培 性变化、最干季降水量和海拔 ９ 个相关系数<０.７ 的
植物大全»(黄宏文ꎬ２０１７)、已出版的«中国迁地植 环境变量运行 Ｒ－Ｋ 模型分析ꎮ 以上分析利用 Ｒ

物志»(黄宏文ꎬ２０１４—２０２１)、«中国栽培植物名录» 软件 包 ａｐｅ (Ｐａｒａｄｉｓ ｅｔ ａｌ.ꎬ ２００４)、 ＭｕＭＩｎ ( Ｂａｒｔｏｎꎬ
(林秦文ꎬ２０１７)ꎬ确认和统计植物园引种栽培的重 ２００９)和 ｒａｓｔｅｒ(Ｈｉｊｍａｎｓ ＆ ｖａｎ Ｅｔｔｅｎꎬ ２０１４)来完成ꎮ

点保护植物的种类和数量ꎮ
１.４ 重点保护植物物种丰富度的空间分布估计 ２ 结果与分析
将江西省以 １ ｋｍ × １ ｋｍ 为空间分析单元划
分为 １ ６８０ 个网格ꎮ 通过空间矢量图层和 １５ ０６５ ２.１ 重点保护植物的多样性以及在自然保护区和
条分布点数据(３６１ 种) 叠加分析ꎬ计算每个网格 植物园的分布
中出现的分布点数量和物种数量ꎮ 分布点数量使 整合形成的江西省重点保护野生维管植物总
用 ｌｏｇ 对数转换并栅格化作图ꎬ表达每个网格区 计有 ４０７ 种ꎬ分属 ８５ 科 ２０８ 属ꎮ 其中ꎬ列入« 国家
１０
域的 调 查 强 度ꎮ 以 广 泛 使 用 的 物 种 丰 富 度 重点保护野生植物名录»(２０２１)的有 １４８ 种ꎬ列入
(ｓｐｅｃｉｅｓ ｒｉｃｈｎｅｓｓꎬ ＳＲ) ( Ｐｒｅｎｄｅｒｇａｓｔ ｅｔ ａｌ.ꎬ １９９３ꎻ «江西省重点保护野生植物名录» (２００５) 的有 ３０５
Ｄａｖｉｅｓ ＆ Ｃａｄｏｔｔｅꎬ ２０１１) 为指标ꎬ计算每个网格出 种(图 １:Ａ)ꎮ 在 ４０７ 种植物中ꎬ石松类和蕨类植
现的物种数量作为物种丰富度的观察值ꎮ 利用地 物 ９ 科 １０ 属 １６ 种ꎬ裸子植物 ６ 科 １７ 属 ２６ 种ꎬ被
统计克里金回归(ｒｅｇｒｅｓｓｉｏｎ－ｋｒｉｇｉｎｇꎬ Ｒ－Ｋ) 模型估 子植物 ７０ 科 １８１ 属 ３６５ 种( 图 １:Ｂ－Ｄꎻ附录 １)ꎮ
计每个网格的物种丰富度ꎬ作为每个网格物种丰 在 ４０７ 种植物中ꎬ６０.７％的种类(２４７ 种) 在江西省
富度的估计值ꎮ Ｒ－Ｋ 模型的基本技术是最佳线性 １４ 个国 家 级 自 然 保 护 区 ( ２３０ 种ꎬ 占 总 种 数 的
无偏估计的空间版本ꎮ 它将目标变量的值作为确 ５６.５％)和 １８ 个省级自然保护区(１３３ 种ꎬ占总种
定性(即空间自相关性)和随机成分的函数进行建 数的 ３２.７％) 都有分布ꎻ７０.５％的种类(２８７ 种) 在
模 ( Ｈｅｎｇｌ ｅｔ ａｌ.ꎬ ２００７ꎻ Ｈｅｎｇｌꎬ ２００９ꎻ Ｃｒｅｓｓｉｅ ＆ 国内 ６７ 个植物园有迁地保护记录( 图 １:Ｅ、Ｆꎻ附
Ｗｉｋｌｅꎬ ２０１１)ꎮ 用物种丰富度估计值减去观察值 录 １)ꎮ 重点保护植物在各自然保护区分布的数量
作为网格单元调查完整性和知识空缺的评估指 由多到少依次是庐山(１３０ 种)、井冈山(８７)、伊山

标ꎬ并对相关结果进行空间作图ꎮ (７１)、九连山(５８)、齐云山(４６)、官山(４５)、马头
由于 Ｒ－Ｋ 模型分析首先需要选择调查相对完 山(４５)、岩泉(３７)、五指峰(３５)、阳际峰(３３)、玉
善的网格来计算半方差函数(ｓｅｍｉ￣ｖａｒｉｏｇｒａｍ)ꎬ因此 京山( ２８)、 武 夷 山 ( ２７)、 南 风 面 ( ２３) 和 云 居 山
我们使用拟合的物种累积曲线来描述每个网格的 (２０)ꎬ其他 １８ 个自然保护区分布的重点保护植物

物种丰富度与调查强度间的关联关系(Ｌｏｂｏ ｅｔ ａｌ.ꎬ 种类数量少于 ２０ 种(附表 １)ꎮ 重点保护植物在各
２０１８ꎻ Ａｌｖｅｓ ｅｔ ａｌ.ꎬ ２０１９)ꎮ 综合从这一过程中得出 植物园分布的数量由多到少依次是中国科学院武
的 ３ 个度量指标:(１)完整性ꎬ即通过物种累积曲线 汉植 物 园 ( ＷＢＧꎬ ２２０ 种 )、 华 南 国 家 植 物 园
估算出的物种丰富度百分比ꎻ(２)比率ꎬ即网格内出 (ＳＣＢＧꎬ２０５)、杭州植物园( ＨＢＧꎬ１８５)、南京中山

现记录数与观察到的物种丰富度之间的比率ꎻ(３) 植物园(ＮＢＧꎬ１７４)、上海辰山植物园( ＣＢＧꎬ１６０)、
斜率ꎬ即观察到的累积物种丰富度与分布点数量之 西双版纳热带植 物 园 ( ＸＴＢＧꎬ１６０)、昆明植物园
间的斜率ꎮ 按照(Ｓｏｕｓａ￣Ｂａｅｎａ ｅｔ ａｌ.ꎬ ２０１４)的方法ꎬ (ＫＢＧꎬ１４２)、 桂林植物园(ＧＸＩＢꎬ１３８)、深圳仙湖植
选择完整性>８５％、比率>４、斜率<０.０２、网格记录数 物园( ＦＬＢＧꎬ１３５) 、东莞植物园( ＤＧＢＧꎬ１３０)、 赣