促进生态服务功能提升的植被覆盖度恢复上限计算方法技术

技术编号:38989160 阅读:13 留言:0更新日期:2023-10-07 10:19
本发明专利技术涉及生态恢复工程规划技术领域,具体涉及一种促进生态服务功能提升的植被覆盖度恢复上限计算方法,包括:获取生态工程区的多组植被及服务数据;计算各组的生态系统综合服务指数IES;以FVC为横坐标,以IES为纵坐标,绘制FVC与IES的关系散点图,并利用滑动窗口沿横坐标正方向依次滑动,得到n个子集;对每个子集中的散点进行局部回归,得到n个局部回归函数;计算n个局部回归函数的中心点值,并用曲线拟合连接,得到FVC与IES的关系曲线;根据FVC与IES的关系曲线估计FVC的上限点;利用重采样方法估计FVC的上限区间。本发明专利技术可弥补生态系统服务对植被覆盖度的响应上限的缺陷,提高结果准确度。准确度。准确度。

【技术实现步骤摘要】
促进生态服务功能提升的植被覆盖度恢复上限计算方法


[0001]本专利技术涉及生态恢复工程规划
,更具体的说是涉及一种促进生态服务功能提升的植被覆盖度恢复上限计算方法。

技术介绍

[0002]植树造林是世界各国普遍采用的一种强有力的生态恢复手段。植被是影响生态系统服务的重要因素,在生态管理实践中,管理者倾向于通过改善植被达到提升生态系统服务的目标。然而,生态系统服务对植被覆盖度的响应是一个复杂的过程,尽管各种研究发现植被覆盖度与生态系统服务之间存在相关关系,但其线性或非线性特征仍不清楚,广泛植树造林而不根据当地降水条件和其他生态限制来选择树种与最佳覆盖度是不可持续的,反而会使生态系统服务功能下降,并造成社会经济后果。因此植被覆盖度应以维持一定水平为宜,生态恢复的人为干预存在一个合理上限,超过该上限,生态系统功能不会显着增加,反而可能会下降。寻找维持最大生态系统服务功能的植被覆盖度的确定方法显得尤为重要,是一个重要的科学问题。
[0003]目前常用的上限确定方法多依托统计分析方法,但这些传统方法存在一些不足。比如:探索性统计分析类的方法只能洞察上限的数量和一般形式的关系,也即只能用来选择模型,随后还需另选合适模型进行拟合。推断性统计分析类的方法有的可能无法检测到低水平的突变点,而有的虽然数据要求相对较低,但总是只找到单个上限,存在忽略考虑真正的函数关系是否包含上限以及包含多少个上限的问题。
[0004]因此,如何实现对突变点数量的判断以及拟合模型的选择,是本领域技术人员亟需解决的问题。
专利
技术实现思路

[0005]有鉴于此,本专利技术提供了一种促进生态服务功能提升的植被覆盖度恢复上限计算方法,可弥补生态系统服务对植被覆盖度的响应上限的缺陷,提高结果准确度。
[0006]为了实现上述目的,本专利技术采用如下技术方案:
[0007]一种促进生态服务功能提升的植被覆盖度恢复上限计算方法,包括以下步骤:
[0008]获取生态工程区的多组植被及服务数据,每组植被及服务数据分别包含植被覆盖度FVC及对应的多个生态系统服务指标;
[0009]分别对每组中的多个所述生态系统服务指标进行加权求和,得到各组的生态系统综合服务指数IES;
[0010]以FVC为横坐标,以IES为纵坐标,绘制FVC与IES的关系散点图,并利用滑动窗口沿横坐标正方向依次滑动,得到n个子集;
[0011]对每个子集中的散点进行局部回归,得到n个局部回归函数;
[0012]计算n个局部回归函数的中心点值,并用曲线拟合连接,得到FVC与IES的关系曲线;
[0013]根据FVC与IES的关系曲线估计FVC的上限点;
[0014]利用重采样方法估计FVC的上限区间。
[0015]进一步的,每组数据中的生态系统服务指标包括:碳固存服务、土壤保持服务和水源涵养服务;
[0016]其中,所述碳固存服务以生态系统碳储量为指标,计算公式为:
[0017];
[0018]式中,C为总碳存储量;C
a
为地上碳存储量;C
b
为地下碳存储量;C
s
为土壤碳存储量;C
d
为死亡有机质碳存储量;
[0019]土壤保持服务指标以土壤保持量作为指标,计算公式为:
[0020];
[0021]式中,S、S
p
和S
v
分别对应表示土壤保持量、潜在土壤侵蚀和实际土壤侵蚀;R表示降水侵蚀因子;K表示土壤侵蚀因子;L表示坡长因子;F表示坡度因子;C表示植被覆盖因子;P表示管理因子;
[0022]水源涵养服务以产水量为指标,计算公式为:
[0023];
[0024]式中,W为产水量;P为降水量;ET为蒸散量。
[0025]进一步的,获取生态工程区的多组植被及服务数据之后,还包括:
[0026]检验植被覆盖度FVC与各生态系统服务指标之间是否存在相关关系。
[0027]进一步的,所述生态系统综合服务指数IES的计算过程包括:
[0028]将各组中的多个所述生态系统服务指标进行处理,消除量纲差异,处理公式为:
[0029];
[0030]式中,Z
r
为第r个消除量纲差异后的生态系统服务指标,r为生态系统服务指标;ES
r
表示第r个生态系统服务数据,ES
r,mean
表示各组数据中第r个生态系统服务指标的平均值;
[0031]利用变异系数法对各生态系统服务指标进行赋权,得到生态系统综合服务指数IES,计算公式为:
[0032];
[0033];
[0034];
[0035]式中,σ
r
为第r个生态系统服务指标的变异系数,r=1、2、3;R为生态系统服务指标的个数;为第r个生态系统服务指标的平均值,D
r
为第r个生态系统服务指标的权重;IES为生态系统综合服务指数。
[0036]进一步的,n个子集的获得过程包括:
[0037]以FVC为横坐标、IES为纵坐标,绘制FVC与IES的关系散点图,得到(FVC1,IES1)~
(FVC
i
,IES
i
)~(FVC
n
,IES
n
)的二维散点图,其中FVC1≤FVC
i
≤FVC
n
,i为散点图中第i个点,i=1、2、3、

、n,n为总散点个数;
[0038]利用一个宽度为s的滑动窗口,将以FVC
i
为中心s为宽度的附近的点组成为一个子集,i=1、2、3、

、n;
[0039]将该滑动窗口沿着横坐标的正方向依次滑动,得到n个子集;其中,滑动窗口的宽度计算方法如下:
[0040];
[0041]式中,s为滑动窗口的宽度,g为子集中散点的个数占总散点个数的比例,n为总散点个数,g取值在0.2~0.7之间。
[0042]进一步的,局部回归函数的中心点值的计算过程,包括:
[0043]利用权重来调整子集中散点对局部回归函数的贡献,局部回归的计算公式如下:
[0044];
[0045];
[0046];
[0047];
[0048]其中,为第i个拟合出的生态系统综合服务指数,即局部回归函数的中心点值;a和b为局部回归函数的参数;FVC
ij
为第i个子集中的第j个点的FVC,j=1、2、3、

、s;为第i个子集中FVC平均值;IES
ij
为第i个子集中第j个点的生态系统综合服务指数;为第i个子集中生态系统综合服务指数的平均值;为第i个子集中第j个点的权重。
[0049]进一步的,所述根据FVC与IES的关系曲线估计FVC的上限点,本文档来自技高网
...

【技术保护点】

【技术特征摘要】
1.一种促进生态服务功能提升的植被覆盖度恢复上限计算方法,其特征在于,包括以下步骤:获取生态工程区的多组植被及服务数据,每组植被及服务数据分别包含植被覆盖度FVC及对应的多个生态系统服务指标;分别对每组中的多个所述生态系统服务指标进行加权求和,得到各组的生态系统综合服务指数IES;以FVC为横坐标,以IES为纵坐标,绘制FVC与IES的关系散点图,并利用滑动窗口沿横坐标正方向依次滑动,得到n个子集;对每个子集中的散点进行局部回归,得到n个局部回归函数;计算n个局部回归函数的中心点值,并用曲线拟合连接,得到FVC与IES的关系曲线;根据FVC与IES的关系曲线估计FVC的上限点;利用重采样方法估计FVC的上限区间。2.根据权利要求1所述的促进生态服务功能提升的植被覆盖度恢复上限计算方法,其特征在于,每组数据中的生态系统服务指标包括:碳固存服务、土壤保持服务和水源涵养服务;其中,所述碳固存服务以生态系统碳储量为指标,计算公式为:;式中,C为总碳存储量;C
a
为地上碳存储量;C
b
为地下碳存储量;C
s
为土壤碳存储量;C
d
为死亡有机质碳存储量;土壤保持服务指标以土壤保持量作为指标,计算公式为:;式中,S、S
p
和S
v
分别对应表示土壤保持量、潜在土壤侵蚀和实际土壤侵蚀;R表示降水侵蚀因子;K表示土壤侵蚀因子;L表示坡长因子;F表示坡度因子;C表示植被覆盖因子;P表示管理因子;水源涵养服务以产水量为指标,计算公式为:;式中,W为产水量;P为降水量;ET为蒸散量。3.根据权利要求1所述的促进生态服务功能提升的植被覆盖度恢复上限计算方法,其特征在于,获取生态工程区的多组植被及服务数据之后,还包括:检验植被覆盖度FVC与各生态系统服务指标之间是否存在相关关系。4.根据权利要求1所述的促进生态服务功能提升的植被覆盖度恢复上限计算方法,其特征在于,所述生态系统综合服务指数IES的计算过程包括:将各组中的多个所述生态系统服务指标进行处理,消除量纲差异,处理公式为:;式中,Z
r
为第r个消除量纲差异后的生态系统服务指标,r为生态系统服务指标;ES
r
表示第r个生态系统服务数据,ES
r,mean
表示各组数据中第r个生态系统服务指标的平均值;利用变异系数法对各生态系统服务指标进行赋权,得到生态系统综合服务指数IES,计算公式为:
;;;式中,σ
r
为第r个生态系统服务指标的变异系数,r=1、2、3;R为生态系统服务指标的个数;为第r个生态系统服务指标的平均值,D
r
为第r个生态系统服务指标的权重;IES为生态系统综合服务指数。5.根据权利要求1所述的促进生态服务功能提升的植被覆盖度恢复上限计算方法,其特征在于,n个子集的获得过程包括:以FVC为横坐标、IES为纵坐标,绘制FVC与IES的关系散点图,得到(FVC1,IES1)~(FVC
i<...

【专利技术属性】
技术研发人员:刘晓曼袁静芳杨子颢章世喜侯鹏刘成周倩高兵兵王超王健付卓黄艳吕娜彭阳尤春赫
申请(专利权)人:生态环境部卫星环境应用中心
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1