CN108388566A - 一种植物物种丰富度分布的空间模拟方法 - Google Patents

Abstract

本发明公开了一种植物物种丰富度分布的空间模拟方法。该方法的实现主要包括三个步骤：首先，在收集和预处理研究区的自然保护区植物物种数据、土地利用/覆盖数据、年平均生物温度和降水数据、DEM数据和土壤数据的基础上，实现各参数因子的空间化处理和统计分析；其次，分别建立植物物种丰富度与土地利用/覆盖面积、环境因子和土地利用/覆盖类型结构指数三大类生境因子之间相关性的最佳方程；最后，实现三大类因子最佳相关方程的综合集成，构建植物物种丰富度的空间模拟方法，实现植物物种丰富度的空间模拟分析。该方法不仅能够实现中宏观尺度上的植物物种丰富度的空间分布模拟，而且具有很好模拟精度和很强的可操作性。

Description

一种植物物种丰富度分布的空间模拟方法

技术领域

本方法涉及一种植物物种丰富度分布的空间模拟方法，适用于全球变化及生态效应研究领域。

背景技术

气候条件、地形地貌、土壤、土地利用/覆盖变化等生境因子之间的互相作用决定着植物物种的空间分布格局。植物通过截留水分与养分循环来影响水土保持以及地表元素循环过程，其空间分布不仅直接决定其他物种生存所必需的物质和能量的空间分布状况，而且植物物种丰富度在一定程度上决定了其他物种的多样性。因此，如何在对气候条件、地形地貌、土壤、土地利用/覆盖变化等生境因子进行综合考虑的情况下，构建植物物种丰富度的空间模拟方法，不仅是在宏观尺度上定量分析和揭示植物物种丰富度的空间格局变化的关键方法和技术，也是区域、全国乃至全球尺度上开展生物多样性和生态安全空间分析研究的焦点内容之一。

植物物种与生境因子的相关关系研究，一直是国内外生态学领域科研究的热点问题。不管是从单一植物物种到种群及群落的定点监测、网络观测、野外调研，还是从单一生境因子到多个生境因子与植物物种分布之间的相关关系分析，都取得了长足进展。譬如，最具代表性的种-面积曲线方法，定量揭示了物种数量随取样面积增加而变化的规律（Arrhenius O, Species and area. Journal of Ecology,1921, 9: 95-99）。Stephenson研究表明气温和降水是影响植物物种类型及其分布的主要气候因子（Stephenson N L.Climatic control of vegetation distribution: the role of the water balance.American Naturalist, 1990, 135(5): 649-670）。Chapin II等人研究表明土地利用/覆盖作为全球变化的主要载体，直接影响地球系统的生物多样性，对陆地生态系统来说，土地利用/覆盖已成为致使物种多样性丧失、改变植物物种空间分布格局的主要驱动因子之一（Chapin II IFS, Zavaleta ES, Eviner VT, et al. Consequences of changingbiodiversity. Nature, 2000, 405: 234–242）。尤其是近些年来，国内外在植物物种空间分布与生境因子相关关系方面开展了大量的研究工作。但主要集中在利用样方采样和历史观测等数据，进行植物物种与生境因子之间关系的统计分析，进而对局地格局上的植物物种丰富度与生境因子之间的关系进行解释，而如何运用GIS空间分析方法，在对离散的植物物种观测数据与相应生境因子之间的相关关系分析的基础上，进行区域尺度乃至更大尺度上的植物物种分布的空间模拟分析研究则相对较少。另外，在考虑生境因子与植物物种空间分布的关系时，虽然有的研究选择了多个生境因子，但在分析过程中，往往分别统计单一生境因子与植物物种之间的相关关系，而缺乏系统地同时考虑多个生境因子与植物物种分布关系的综合分析。

另外，在对国家自然保护区的植物物种数量、土地利用/覆盖及其他环境因子分析的过程中发现，自然保护区土地利用/覆盖动态变化强度低于非保护区，而且自然保护区的植物物种丰富度的空间分布与自然保护区内特定环境因子、土地利用/覆盖类型及其结构指数具有很好的相关性特征。因此，本发明针对上述现有方法的局限性，综合考虑气候条件、地形地貌、土壤、土地利用/覆盖变化等生境因子，利用目前植物物种资料最全的自然保护区调研、观测和统计数据，发明了一种植物物种丰富度分布的空间模拟方法，并实现植物物种丰富度的空间模拟分析，对我国生物多样性保护战略实施、国际生物多样性公约谈判、以及生物多样性和生态系统服务政府间科学政策平台（IPBES）合作研究具有重要意义。

发明内容

本发明旨在针对自然保护区的植物物种丰富度的空间分布与自然保护区内特定环境因子、土地利用/覆盖类型面积及其结构指数具有很强的相关关系这一规律，以自然保护区的植物物种数据作为基准模拟情况和样本数据，充分利用现有自然保护区的植物物种调研、观测和统计数据，在对其植物物种分布丰富度与土地利用/覆盖类型（耕地、林地、草地、湿地、荒漠、建设用地）面积、环境因子（年平均降水、年平均生物温度、地形起伏度、土壤厚度、土壤有机质）、土地利用/覆盖类型结构指数（斑块连通性、破碎度、生态多样性）等三大类生境因子之间的相互作用和影响机理进行定量统计分析的基础上，筛选和确定最佳相关分析方程，构建植物物种丰富度的空间模拟分析模型，进而实现植物物种丰富度的空间模拟分析。为了实现上述目的，本发明实现的关键技术方案包括以下几个步骤：

第一步、数据收集与预处理

1-1）在收集研究区国家自然保护空间范围图及其植物物种数量的基础上，运用GIS空间制图功能，实现各个自然保护区边界矢量数字化，并结合土地利用/覆盖空间数据、DEM数据、气候因子数据和土壤数据，构建植物物种丰富度空间模拟分析模型所需基础参数的空间数据库。

1-2）运用高精度曲面建模（HASM）方法，分别实现各自然保护区年平均生物温度、年平均降水、土壤有机质、土壤厚度的空间模拟。

1-3）运用GIS空间分析方法，分别求算各个自然保护区类内每一种土地利用/覆盖类型的面积，及其结构的斑块连通性、破碎度和生态多样性指数，以及计算地形起伏度指数。

第二步、构建植物物种丰富度与生境因子相关性的最佳方程

2-1）分别采用线性相关和多种非线性相关方程对自然保护区的植物物种丰富度与土地利用/覆盖类型（耕地、林地、草地、湿地、荒漠、建设用地）面积之间的相关性进行定量求解，并以相关系数最高者对应的拟合方程作为二者之间的最佳相关方程。

2-2）分别采用线性相关和多种非线性相关方程对自然保护区的植物物种丰富度与环境因子（年平均生物温度、年平均降水、土壤有机质、土壤厚度及地形起伏度）之间的相关关系进行拟合求算，并以相关系数最高者对应的拟合方程作为二者之间的最佳相关方程。

2-3）分别采用线性相关和多种非线性相关方程对自然保护区的植物物种丰富度与土地利用/覆盖结构指数（斑块连通性、破碎度及生态多样性）的相关关系进行拟合求算，并以相关系数最高者对应的拟合方程作为二者之间的最佳相关方程。

第三步、构建植物物种丰富度空间模拟分析模型

3-1）运用第二步获得的植物物种丰富度与土地利用/覆盖类型面积、环境因子和土地利用/覆盖结构指数三大类型因子的最佳相关拟合方程，利用空间相关分析方法，建立植物物种丰富度与生境因子之间的最佳空间相关方程，进而构建研究区的植物物种丰富度空间模拟模型。

3-2）基于GIS空间模拟分析平台和运用各生境因子的空间栅格数据，调试和运行植物物种丰富度空间模拟模型，并在利用部分国家自然保护区现有植物物种丰富度数据进行精度验证和误差分析的基础上，实现研究区植物物种丰富度的空间模拟分析。

具体实施方式

为了使本发明的目的、技术方案更加清晰，以下对本发明进一步详细说明。应当理解，此处所描述的具体实时方式仅仅用以解释本发明，并不用于限定本发明。

本发明公开了一种植物物种丰富度分布的空间模拟方法，包括以下步骤：

第一步、数据收集与预处理；

第二步、构建植物物种丰富度与生境因子相关性的最佳方程；

第三步、构建植物物种丰富度空间模拟分析模型。

以下对具体步骤进行详述：

研究区数据收集与预处理。模型构建所需的基础数据包括研究区植物物种数据、自然保护区边界、土地利用/覆盖数据、气候数据、土壤数据和DEM数据。首先，运用GIS空间制图功能，实现各自然保护区边界的空间矢量化，并结合文献资料、网络搜索、野外调研等获取的植物物种数据，构建研究区自然保护区及其植物物种数据的空间数据库；其次，运用高精度曲面建模（HASM）方法，对气候数据、土壤有机质和土壤厚度数据进行高精度空间模拟，获得年平均生物温度、年平均降水、土壤有机质和土壤厚度的空间栅格数据。HASM方法是为了解决现有曲面建模方法的误差问题，运用曲面论基本定理进行曲面建模的方法。针对空间网格,,，基于高斯方程HASM方法可表达为等式约束的最小二乘问题（Yue TX. Surface Modelling: High Accuracy and High Speed Methods. NewYork: CRC Press, 2011）：

(1)

式中，和分别为采样矩阵和采样向量，为采样点个数。

再次，利用GIS空间统计分析方法，实现研究区内各自然保护区内各土地覆盖类型（耕地、林地、草地、湿地、荒漠、建设用地）的面积、年平均生物温度、年平均降水、地形起伏度、土壤有机质和土壤厚度进行空间统计求算；最后，选取相关的指数模型，分别对研究区内各自然保护区的土地利用/覆盖结构指数（斑块连通性、破碎度和多样性）进行计算，理论公式分别如下：

(2)

(3)

（4）

式中，为土地覆盖斑块连通性指数，为第i种土地覆盖类型的第j个斑块的面积占研究区域总面积的比率，为第i种土地利用/覆盖类型第j个斑块的面积，为第i种土地利用/覆盖类型第j个斑块的周长，为正六边形周长的平方与其面积的比值，为土地利用/覆盖的景观破碎度，为土地利用/覆盖的总类型数，为某种土地利用/覆盖类型的斑块数，为土地利用/覆盖景观多样性指数，为第i种土地利用/覆盖类型面积的概率，即其占总研究区土地利用/覆盖面积的比例，为土地利用/覆盖类型总数目，为观测尺度，为研究区总面积，=2.71828，为时间变量。

植物物种丰富度与生境因子相关性的最佳方程构建。在对植物物种丰富度分布与生境因子之间的相关性进行定量分析的过程中，将植物物种丰富度分别与各土地利用/覆盖类型面积、环境因子和土地利用/覆盖结构指数三大类因子进行空间相关分析，进而确定与各大类因子之间的相关性拟合的最佳方程。

首先，对自然保护区的植物物种丰富度与土地利用/覆盖类型（耕地、林地、草地、湿地、荒漠、建设用地）面积之间的相关性分别进行多种回归拟合，在各个拟合方程进行对比分析的基础上，确定二者之间的最佳相关回归方程为：

(5)

式中，代表基于土地利用/覆盖类型面积求解的植物物种丰富度，为植物物种丰富度与各种土地利用/覆盖类型之间的方程系数，为常量，其余参数含义同上。

其次，对自然保护区的植物物种丰富度与环境因子（年平均降水、年平均生物温度、地形起伏度、土壤厚度、土壤有机质）之间的相关性分别进行多种回归拟合，在各个拟合方程进行对比分析的基础上，确定二者之间的最佳相关回归方程为：

(6)

式中，代表基于环境因子求解的植物物种丰富度，为植物物种丰富度与各种环境因子之间的方程系数，为常量，其余参数含义同上。

再次，对自然保护区的植物物种丰富度与土地利用/覆盖结构指数（斑块连通性、破碎度及生态多样性）之间的相关性分别进行多种回归拟合，在各个拟合方程进行对比分析的基础上，确定二者之间的最佳相关回归方程为：

(6)

式中，代表基于土地利用/覆盖类型结构指数求解的植物物种丰富度，为植物物种丰富度与各种土地利用/覆盖结构指数之间的方程系数，和为常量，其余参数含义同上。

植物物种丰富度空间模拟分析模型的构建。首先，在分别获得植物物种丰富度与土地利用/覆盖类型面积、环境因子和土地利用/覆盖结构指数三大类型因子的最佳相关拟合方程的基础上，对三大类因子的最佳拟合方程进行综合集成，建立植物物种丰富度与生境因子之间的最佳空间相关方程，实现研究区植物物种丰富度空间分布的模拟分析。其理论计算公式可表达为：

(7)

式中，为t时段空间模拟单元（模拟栅格）上的植物物种丰富度；为t时段基于土地利用/覆盖类型面积模拟的空间模拟单元（模拟栅格）上植物物种丰富度的模拟值；为t时段基于环境因子变量模拟的空间模拟单元（模拟栅格）上植物物种丰富度的模拟值；为t时段基于土地利用/覆盖类型结构指数模拟的空间模拟单元（模拟栅格）上植物物种丰富度的模拟值；分别为植物物种丰富度与土地利用/覆盖类型面积、环境因子和土地利用/覆盖结构指数三大类型因子之间的最佳空间相关方程的权重系数。

其次，基于GIS空间模拟分析平台和各生境因子的空间栅格数据，调试和运行方程（7），利用国家自然保护区现有植物物种丰富度数据进行模拟，并对精度进行验证和误差分析，最终实现研究区植物物种丰富度的空间模拟分析。

利用本发明构建了植物物种丰富度空间模拟方法，对青藏高原维管植物物种丰富度的空间分布进行了模拟测试分析。测试结果显示，利用方程（7）对青藏高原维管植物物种与土地利用/覆盖类型面积、环境因子和土地利用/覆盖结构指数三大类型因子之间的复相关系数高达0.93，而且满足0.01的显著性相关检验。误差分析表明，该方法对青藏高原35个国家自然保护区维管植物物种丰富度的模拟误差仅约为2种/km²。模拟结果显示，对整个青藏高原范围内的维管植物物种丰富度分布的空间模拟值与该区域的维管植物物种丰富度的实际分布情况具有很好的一致性。以上案例分析表明，该方法不仅能够实现中宏观尺度上的植物物种丰富度的空间分布模拟，而且具有很好模拟精度和很强的可操作性。

Claims (4)

1.一种植物物种丰富度分布的空间模拟方法，主要包括以下三个步骤。

2.第一步、数据收集与预处理

1）在收集研究区国家自然保护空间范围图及其植物物种数量的基础上，运用GIS空间制图功能，实现各个自然保护区边界矢量数字化，并结合土地利用/覆盖空间数据、DEM数据、气候因子数据和土壤数据，构建植物物种丰富度分布空间模拟分析模型所需基础参数的空间数据库；

2）运用高精度曲面建模（HASM）方法，分别实现各自然保护区年平均生物温度、年平均降水、土壤有机质、土壤厚度的空间模拟；

3）运用GIS空间分析方法，分别求算各个自然保护区类内每一种土地利用/覆盖类型的面积，及其结构的斑块连通性、破碎度和生态多样性指数，以及计算地形起伏度指数。

3.第二步、构建植物物种丰富度与生境因子相关性的最佳方程

1）分别采用线性相关和多种非线性相关方程对自然保护区的植物物种丰富度与土地利用/覆盖类型（耕地、林地、草地、湿地、荒漠、建设用地）面积之间的相关性进行定量求解，并以相关系数最高者对应的拟合方程作为二者之间的最佳相关方程；

2）分别采用线性相关和多种非线性相关方程对自然保护区的植物物种丰富度与环境因子（年平均生物温度、年平均降水、土壤有机质、土壤厚度及地形起伏度）之间的相关关系进行拟合求算，并以相关系数最高者对应的拟合方程作为二者之间的最佳相关方程；

3）分别采用线性相关和多种非线性相关方程对自然保护区的植物物种丰富度与土地利用/覆盖结构指数（斑块连通性、破碎度及生态多样性）的相关关系进行拟合求算，并以相关系数最高者对应的拟合方程作为二者之间的最佳相关方程。

4.第三步、构建植物物种丰富度空间模拟分析模型

1）运用第二步获得的植物物种丰富度与土地利用/覆盖类型面积、环境因子和土地利用/覆盖结构指数三大类型因子的最佳相关拟合方程，利用空间相关分析方法，建立植物物种丰富度与生境因子之间的最佳空间相关方程，进而构建研究区的植物物种丰富度空间模拟模型；

2）基于GIS空间模拟分析平台和运用各生境因子的空间栅格数据，调试和运行植物物种丰富度空间模拟模型，并在利用部分国家自然保护区现有植物物种丰富度数据进行精度验证和误差分析的基础上，实现研究区植物物种丰富度的空间模拟分析。