1.本发明涉及生态保护与恢复技术领域，具体而言，涉及一种生态系统恢复力的评价方法。
背景技术：
2.生态系统恢复力是生态系统受到扰动后恢复到稳定状态的能力，包括维持其重要特征，如生物组成、结构、生态系统功能和过程速率的能力。保持和提高生态系统恢复力是区域可持续发展的基础。受气候变化和人类活动等影响，生物多样性减少、生态系统生产力下降等导致了区域乃至全球生态系统恢复力的不断下降，成为制约生态系统可持续发展能力因素之一。
3.生态系统的多稳态机制导致系统状态受到外部干扰时会发生突变而进入不理想的状态，因而通过定量评估提取出恢复力的主要影响因子从而确定生态系统恢复力是恢复力定量评价的重要途径，而目前存在评价方法不够全面，导致恢复力定量评价结果不够准确，因此，提出一种生态系统恢复力的评价方法。
技术实现要素：
4.本发明的目的在于提供一种生态系统恢复力的评价方法，以解决上述背景技术中存在的问题。
5.本发明的实施例是这样实现的：本技术实施例提供了一种生态系统恢复力的评价方法，其包括以下步骤：根据不同类型的生态系统选取若干一级指标，分别将若干一级指标划分为若干的二级指标，并归为恢复力评价因子；将所述二级指标进行归一化处理，计算公式如下：；式中：e 为恢复力水平；n 为评价因子的个数；k 为恢复力评价因子（k=1,2,3
…
）； fk 为第k 个恢复力评价因子的水平值；wk 为第k 个恢复力评价因子的权重，并采用层次分析法对k 进行赋值。
6.在本发明的一些实施例中，上述所述一级指标包括生态存储指标，且所述生态存储指标包括：生态系统净初级生产力指标，并基于casa 光能利用率模型进行计算，其计算公式如下：；式中，apar(x,t)表示像元x 在t 月吸收的光合有效辐射，表示像元x 在t 月的实际光能利用率；上述apar(x,t)的计算公式如下：
；式中：sol(x,t)表示t月在像元x处的太阳总辐射量，fpar(x,t)植被层对入射光合有效辐射的吸收比例，常数0.5表示植被所能利用的太阳有效辐射占太阳总辐射的比例；上述fpar(x,t)的计算公式如下：或；式中：ndvi
i,max
和ndvi
i,min
分别对应第i种植被类型的ndvi最大和最小值；fpar
min
和fpar
max
的取值与植被类型无关，分别为0.001和0.95； sr
i,max
和sr
i,min
分别对应第i种植被类型ndvi的95%和5%下侧百分位数；上述sr(x,t)由以下公式表示：；在上述fpar(x,t)的两计算公式，取其加权平均或平均值作为fpar(x,t)的值，即：；在公式：中，；式中：和表示低温和高温对光能利用率的胁迫作用；为水分胁迫影响系数，反映水分条件的影响；是理想条件下的最大光能利用率；上述由以下公式表示：；式中：topt(x) 为植物生长的最适温度，定义为某一区域一年内ndvi 值达到最高时的当月平均气温℃，并当某一月平均温度小于或等于-10℃时，其值取0；上述由以下公式表示：；式中，当某一月平均温度t(x, t)比最适温度t
opt
(x)高10℃或低13℃时，月平均温度t(x, t)为最适温度t
opt
(x)时，该月的值等于最适温度t
opt
(x)值的一半；上述由以下公式表示：；式中：et为区域实际蒸散量（mm）； pet为区域潜在蒸散量（mm）；
植被指数ndvi指标，其计算公式如下：ndvi=（nir-r）/（nir+r）以及；式中：nir 为近红外波段，r 为可见光红波段；ndviveg 为完全植被覆盖地表所贡献的信息，ndvisoil 为无植被覆盖地表所贡献的信息。
7.在本发明的一些实施例中，上述所述一级指标还包括生态系统服务指标，且所述生态系统服务指标包括：水土保持指标，并使用修正通用水土流失方程的功能算法进行计算，其模型结构为：a
c = a
p
ꢀ‑ꢀ
a r
＝ r
×ꢀk×ꢀ
l
×s×
（1
‑ꢀ
c）；式中： ac为土壤保持量，a
p
为潜在土壤侵蚀量，ar为实际土壤侵蚀量，r为降水因子，k为土壤侵蚀因子，l、s为地形因子，c为植被覆盖因子；其中，r－降水因子表示为：；式中：pi为月均降水量，p为年均降水量,α=0.3589，β= 1.9462；k－土壤侵蚀因子表示为：k = f
csand
×fcl-si
×forgc
×fhisand
；而上式中：f
csand = 0.2 + 0.3 exp [-0.0256 ms (1-msilt/100)]；f
cl-si = [msilt/(mc + msilt)]
0.3
；f
orgc = 1
‑ꢀ
0.25orgc / [orgc + exp(3.72
ꢀ‑ꢀ
2.95 orgc)]；f
hisand = 1-0.7(1-ms/100)/{(1-ms/100)+exp[-5.51+22.9(1-ms/100)]}；式中：ms为土壤粗砂含量，msilt为土壤粉砂含量，mc为土壤粘粒含量，orgc为有机碳含量；l、s—地形因子表示为：l = (λ/22.13)m；s = 10.8sinθ + 0.03
ꢀꢀꢀꢀꢀꢀ
θ＜5
°
；16.8sinθ— 0.5
ꢀꢀꢀꢀꢀꢀ5°
≤θ＜10
°
；21.91sinθ— 0.96
ꢀꢀꢀꢀ
θ≥ 10
°
；式中：λ为坡长（m），m为坡长指数，θ为坡度（
°
）；水源涵养指标，并采用水量平衡方程来计算，其计算公式为：；式中：tq为总水源涵养量（m3），p 为降雨量（mm），r 为地表径流量（mm），et 为蒸散发（mm），a 为 i类生态系统面积（km2），i为研究区第 i类生态系统类型，j为研究区生态系统类型数。
[0008]
在本发明的一些实施例中，上述所述一级指标还包括生态敏感性指标，且所述生态敏感性指标包括：水土流失敏感性指标，计算公式如下：；
式中：ssi 为i 空间单元水土流失敏感性指数，评价因子包括降雨侵蚀力、土壤可蚀性、坡长坡度、地表植被覆盖；土地沙化敏感性指标，计算公式如下：；式中：di为i评价区域土地沙化敏感性指数；ii、wi、ki、ci分别为评价区域干燥指数、起沙风天数、土壤质地和植被覆盖的敏感性等级值；石漠化敏感性指标，计算公式如下：；式中：si为i评价区域石漠化敏感性指数；di、pi、ci为i评价区域碳酸盐出露面积百分比、地形坡度和植被覆盖度；di是区域单元范围内碳酸盐出露面积占单元总面积的百分比；pi根据评价区数字高程（dem）在地理信息系统下进行处理和分级；ci的数据来源和处理方法参照土地沙化敏感性。
[0009]
相对于现有技术，本发明的实施例至少具有如下优点或有益效果：通过本技术中的方式来构建生态系统恢复力测度指标体系，并展开实证分析，旨在为区域可持续发展研究提供理论与实践指导，为生态修复提供方向；相比目前存在的生态系统恢复力的评价方式，本技术中根据不同的生态系统而选取若干指标，其囊括的范围更大、更加全面，能够较精准的评定出不同生态系统的恢复力。
具体实施方式实施例
[0010]
本技术实施例提供了一种生态系统恢复力的评价方法，具体方式如下：目前大量学者探究了生态系统恢复力的影响因素，并积累了较多的区域案例。但整体而言，目前对生态系统恢复力影响因素的理解依然有限，缺少一致的观点和表征生态系统恢复力的指标体系，而已有文献中的恢复力指标没得到充分利用，根据大量文献对生态恢复的影响因素分析和结合不同地方的实际情况，例如四川省，本实施例中通过林地生态系统和草地生态系统进行进一步的说明，而林地选取8个一级指标，33个二级指标对四川省林地生态系统恢复力进行评价，以及4个修正指标进行评价；将二级指标进行归一化处理，并采用限制系数法计算恢复力水平。计算公式如下：；式中：e 为恢复力水平；n 为评价因子的个数；k 为恢复力评价因子（k=1,2,3
…
）； fk 为第k 个恢复力评价因子的水平值；wk 为第k 个恢复力评价因子的权重，并采用层次分析法对k 进行赋值。
[0011]
具体评价因子如下：林地恢复力水平评价因子构成
而草地选取8个一级指标，32个二级指标对四川省草地生态系统恢复力进行评价，以及4个修正指标进行评价；具体评价因子如下：草地恢复力水平评价因子构成而上述各指标的数据可通过以下的渠道获得：指标数据来源
而上述的一级指标中，地形地貌：地形地貌及其特征是植被恢复能力的关键，一般地，高山、极高山、高原地区，气候恶劣、空气稀薄、极端气候频发，海拔越高植被恢复能力越差，同时，太阳辐射对植被生长过程进行光合作用较大，阴坡半阴坡地区的光照时间短于阳坡，植被恢复能力相对差，地形还会对径流产生影响，从而影响地形湿度。
[0012]
气候条件：气候主要是通过温度、光照和降水对植物的光合作用及呼吸作用的速率及其它生理过程产生影响。在温度适宜和阳光照射充沛条件下，植物光合作用的速率与其水分的吸收成正比。在湿度适宜条件下，植物呼吸作用的速率与温度成正比。
[0013]
蒸散发：蒸散发量为植物蒸腾量与株间土壤蒸发量之和。蒸散发是一个复杂的物理和生理过程，与作物种类、生育期、生长状况、土壤含水量及气象条件有密切关系。蒸散发采用彭曼公式进行计算。
[0014]
干旱指数，为年蒸发能力和年降水量的比值。
[0015]
公式为：r=e0/p；是反映气候干旱程度的指标。对林草、农田植被生长产生影响。
[0016]
干燥度指数表征一个地区干湿程度的指标，也称干燥度指数，一般以某个地区水分收支与热量平衡的比值来表示, 其倒数称为湿润指数(humidity index, hi)。干燥度表达公式如下：k=w0/r；式中：k为干燥度；w0水面可能蒸发量；r为同期降水量。以k值为1.0的等值线来区分湿润地区和半湿润地区的指标。k《1为湿润地区，1≤k≤1.25为半湿润地区。这条等值线大致相当于秦岭——淮河线。但此线以北的大、小兴安岭，长白山地，胶东半岛的k值亦小于1。k《1的地区表示降水量大于可能蒸发量，植被为森林；1≤k≤1.25为半湿润区，植被为森林草原，草甸草原，土壤中有一些石灰质积累，有些地区有盐渍现象，旱患频率较大，1.25＜k≤4为半干旱地区，包括广大的西北地区和青藏高
原，半干旱地区与干草原区相当。 k》4为干旱地区，植被为干荒漠，发展农业需进行灌溉;土壤中含有较多的可溶性盐类。
[0017]
土壤条件：土壤条件包括土壤理化性质。土壤有机质含量和土壤全氮以及土壤腐殖质是土壤肥力的重要指标，土壤肥力越好，越利于林草植被的生长。
[0018]
土壤湿度表征的是土壤含水量，一般土壤含水量越高，越利于植被的生长。
[0019]
土壤结皮指数：生物结皮又称生物土壤结皮、土壤微生物结皮等，其是由微细菌、真菌、藻类、地衣、苔藓等隐花植物及其菌丝、分泌物等与土壤砂砾粘结形成的复合物，是干旱半干旱区重要的地表覆盖类型（40%以上）。生物结皮也是干旱半干旱沙漠最具有特色的生物景观之一，生物结皮的存在对沙漠的固定、土壤表面的物理化学生物学特性、土壤抗风蚀水蚀等方面具有重要意义。生物结皮也是沙漠植被演替的先锋种，对促进沙漠植被演化具有重要作用。土壤结皮有利于抵抗风蚀，防止土地沙化，是干旱半干旱区域植物生态恢复力的重要指标。
[0020]
人工干扰：人类的生产和活动一般对植被生长会产生一定的影响，依据土地利用类型，将建设用地、耕地、园地和其他土地利用类型作为负向因子进行影响力分析，对建设用地和耕地做欧氏距离分析，距离越远，影响越小；园地对周边植被的影响相对较小，因此，将园地和非园地分为两类区；最后将所有图层做空间叠加，得到人工干扰度。
[0021]
生态存储：生态存储是影响生态系统恢复力的重要因素，恢复力是由其生态系统中生物和环境决定的，因而保证系统的生态存储具有重要意义。生态存储通常分为内部存储和外部存储。内部存储是指受到干扰后，生态系统所存在的原有生物及其环境的重新组合的能力，外部存储是指受到干扰后，借助系统外围或没有受到干扰的区域内生物进行重新组合的能力。对生态系统恢复力进行评价，要确定生态系统内外部存储的数量、种类、分布以及其他因素的影响。
[0022]
生态存储表征生态系统的能量和物质基础，一般的，生态存储越丰富，生态系统的抵抗力越高，受到破坏时越容易恢复。针对森林生态系统，主要考虑林业数据中的优势树种、郁闭度、龄级、平均胸径、单位面积株数等因子；针对草地生态系统，主要考虑植被净初级生产力、植被覆盖度因子。
[0023]
生态系统净初级生产力（npp）可基于casa光能利用率模型计算。casa模型认为npp是由植物光合作用与其对光能利用率的大小共同决定的。所以，casa模型中npp的估算可以由植物的光合有效辐射（apar）和实际光能利用率（ε）两个因子来表示，其估算公式如下：；式中，apar(x,t)表示像元x在t月吸收的光合有效辐射（g c
·
m-2
·
month-1
），ε(x,t)表示像元x在t月的实际光能利用率（g c
·
mj-1
）。
[0024]
apar的估算：apar的值由植被所能吸收的太阳有效辐射和植被对入射光合有效辐射的吸收比例来确定。
[0025]
；
式中：sol(x,t)表示t月在像元x处的太阳总辐射量，fpar(x,t)植被层对入射光合有效辐射的吸收比例，常数0.5表示植被所能利用的太阳有效辐射占太阳总辐射的比例。
[0026]
fpar的估算：由于在一定范围内，fpar与ndvi之间存在着线性关系，这一关系可以根据某一植被类型ndvi的最大值和最小值以及所对应的fpar最大值和最小值来确定。
[0027]
；式中：ndvi
i,max
和ndvi
i,min
分别对应第i种植被类型的ndvi最大和最小值。fpar与比值植被指数（sr）也存在着较好的线性关系，可由以下公式表示：；式中，fpar
min
和fpar
max
的取值与植被类型无关，分别为0.001和0.95； sr
i,max
和sr
i,min
分别对应第i种植被类型ndvi的95%和5%下侧百分位数， sr(x,t)由以下公式表示：；通过对fpar-ndvi和fpar-sr所估算结果的比较发现，由ndvi所估算的fpar比实测值高，而由sr所估算的fpar则低于实测值，但其误差小于直接由ndvi所估算的结果，因此我们可以将二者结合起来，取其加权平均或平均值作为估算fpar的估算值：；光能利用率的估算：光能利用率是在一定时期单位面积上生产的干物质中所包含的化学潜能与同一时间投射到该面积上的光合有效辐射能之比。环境因子如气温、土壤水分状况以及大气水汽压差等会通过影响植物的光合能力从而调节植被的npp。
[0028]
；式中：t
ε1
(x,t)和t
ε2
(x,t)表示低温和高温对光能利用率的胁迫作用；w
ε
(x,t)为水分胁迫影响系数，反映水分条件的影响；ε
max
是理想条件下的最大光能利用率（g c/mj）。
[0029]
温度胁迫因子的估算：t
ε1
(x, t)的估算：其反映在低温和高温时植物内在的生化作用对光合的限制而降低第一性生产力。
[0030]
；式中：t
opt
(x)为植物生长的最适温度，定义为某一区域一年内ndvi值达到最高时的当月平均气温（℃）；当某一月平均温度小于或等于-10℃时，其值取0。 t
ε2
(x, t)的估算：表示环境温度从最适温度t
opt
(x)向高温或低温变化时植物光能利用率逐渐变小的趋势，这是因为低温和高温时高的呼吸消耗必将会降低光能利用率，生长在偏离最适温度的条件下，其光能利用率也一定会降低。
[0031]
；当某一月平均温度t(x, t)比最适温度t
opt
(x)高10℃或低13℃时，月平均温度t(x, t)为最适温度t
opt
(x)时，该月的t
ε2
(x, t)值等于最适温度t
opt
(x)值的一半。
[0032]
水分胁迫因子的估算：水分胁迫影响系数w
ε
(x,t)反映了植物所能利用的有效水
分条件对光能利用率的影响，随着环境中有效水分的增加，w
ε
(x,t)逐渐增大,它的取值范围为0.5(在极端干旱条件下)到1(非常湿润条件下)。
[0033]
；式中：et为区域实际蒸散量（mm）； pet为区域潜在蒸散量（mm）；最大光能利用率的确定：月最大光能利用率ε
max
的取值因不同的植被类型而有所不同，在casa模型中全球植被的最大光能利用率为0.389 g c
·
mj-1
。
[0034]
全国尺度上利用casa模型计算npp时主要采用modis250米每16天合成的ndvi数据产品和反射率数据产品，以及modis1千米每8天合成的陆地表面温度(lst)数据产品，并利用气象站点太阳辐射量数据作为遥感估测模型的补充，通过克里金插值得到250米分辨率日照时数数据，以实现npp估算由点向面的推演和扩展。在数据条件允许的情况下，可使用更高时空分辨率数据进行计算。
[0035]
植被覆盖度地表植被覆盖是影植被恢复能力的一个重要因素，在植被覆盖高的地区，一般生物多样性丰富，生态系统受损后恢复能力强；相反，植被稀少区域生物多样性低下、土壤较贫瘠，生态系统恢复能力弱。
[0036]
归一化植被指数ndvi，又称标准化植被指数，定义为近红外波段nir（0.7-1.1微米）与可见光红波段r （ 0.4-0.7微米）数值之差和这两个波段数值之和的比值，即公式：0.7微米）数值之差和这两个波段数值之和的比值，即公式：式中：ndvi
veg
为完全植被覆盖地表所贡献的信息，ndvi
soil
为无植被覆盖地表所贡献的信息。
[0037]
运用地理信息系统软件进行图像处理，获取植被ndvi影像图，进而计算植被覆盖度。由于大部分植被覆盖类型是不同植被类型的混合体，所以不能采用固定的ndvi
soil
和ndvi
veg
值，通常根据ndvi的频率统计表，计算ndvi的频率累积值，累积频率为5%的ndvi值为ndvi
soil
，累积频率为58%的ndvi值为ndvi
veg
。
[0038]
生态系统服务（ecosystem services）是指人类从生态系统获得的所有惠益，一般的，生态系统服务功能越重要，其生态系统恢复力就越强。主要考虑生物多样性、水土保持、水源涵养生态系统服务功能。
[0039]
生物多样性维护功能重要性生物多样性维护功能重要性在生态系统、物种和遗传资源三个层次进行评价。
[0040]
在生态系统层次，将原真性和完整性高，需优先保护的森林、灌丛、草地、湿地、荒漠、内陆湿地、海洋生态系统评定为生物多样性维护极重要区，其他需保护的生态系统评定为生物多样性维护重要区。
[0041]
在物种层次，以具有重要保护价值的物种为保护目标，将国家重点保护野生动植物和列入iucn红色名录中的极危、濒危物种的集中分布区、极小种群野生动植物的主要分布区评定为生物多样性维护极重要区，将省级重点保护物种等其他具有重要保护价值物种的集中分布区评定为生物多样性维护重要区。
[0042]
在遗传资源层次，将重要的野生农作物、水产、畜牧重要种质资源的主要天然分布
区评定为生物多样性维护极重要区。
[0043]
修正通用水土流失方程（usle）的水土保持服务功能算法。
[0044]
模型结构：a
c = a
p
ꢀ‑ꢀ
a r
＝ r
×ꢀk×ꢀ
l
×s×
（1
‑ꢀ
c）；式中： ac为土壤保持量，a
p
为潜在土壤侵蚀量，ar为实际土壤侵蚀量，r为降水因子，k为土壤侵蚀因子，l、s为地形因子，c为植被覆盖因子。
[0045]
r－降水因子；式中：pi为月均降水量，p为年均降水量, α= 0.3589，β= 1.9462。
[0046]
k－土壤侵蚀因子k = f
csand
×fcl-si
×forgc
×fhisandfcsand = 0.2 + 0.3 exp [-0.0256 ms (1-msilt/100)]f
cl-si = [msilt/(mc + msilt)]
0.3forgc = 1
‑ꢀ
0.25orgc / [orgc + exp(3.72
ꢀ‑ꢀ
2.95 orgc)]f
hisand = 1-0.7(1-ms/100)/{(1-ms/100)+exp[-5.51+22.9(1-ms/100)]}式中：ms为土壤粗砂含量，msilt为土壤粉砂含量，mc为土壤粘粒含量，orgc为有机碳含量。
[0047]
l、s—地形因子l = (λ/22.13)ms = 10.8sinθ + 0.03
ꢀꢀꢀꢀꢀꢀ
θ＜5
°
16.8sinθ— 0.5
ꢀꢀꢀꢀꢀꢀ5°
≤θ＜10
°
21.91sinθ— 0.96
ꢀꢀꢀꢀ
θ≥ 10
°
式中：λ为坡长（m），m为坡长指数，θ为坡度（
°
）。
[0048]
水源涵养功能重要性评估通过降雨量减去蒸散量和地表径流量得到的水源涵养量，评价生态系统水源涵养功能的相对重要程度。降雨量较大、蒸散量及地表径流量较小的区域，水源涵养功能重要性较高。森林、灌丛、草地和湿地生态系统质量较高的区域由于地表径流量小，水源涵养功能相对较高。一般地，将累积水源涵养量最高的前50%区域划分为极重要区。在此基础上，结合大江大河源头区、饮用水水源地等边界进行适当修正。
[0049]
水源涵养是生态系统（如森林、草地等）通过其特有的结构与水相互作用，对降水进行截留、渗透、蓄积，并通过蒸散发实现对水流、水循环的调控，主要表现在缓和地表径流、补充地下水、减缓河流流量的季节波动、滞洪补枯、保证水质等方面。以水源涵养量作为生态系统水源涵养功能的评估指标。
[0050]
评估模型采用水量平衡方程来计算水源涵养量，计算公式为：；式中：tq为总水源涵养量（m3），p 为降雨量（mm），r 为地表径流量（mm），et 为蒸散
发（mm），a 为 i类生态系统面积（km2），i为研究区第 i类生态系统类型，j为研究区生态系统类型数。
[0051]
数据准备（1）数据来源与获取根据上述模型，水源涵养功能重要性评估需收集生态系统类型数据集、气象数据集和蒸散发数据集等，具体信息如下：水源涵养功能重要性评估数据表数据预处理降雨量因子：根据气象数据集处理得到。在 excel 中计算出区域所有气象站点的多年平均降水量，将这些值根据相同的站点名与 arcgis 中的站点（点图层）数据相连接（join）。在 spatial analyst 工具中选择 interpolate to raster 选项，选择相应的插值方法得到降水量因子栅格图。
[0052]
地表径流因子：降雨量乘以地表径流系数获得，计算公式如下：；式中：r为地表径流量（mm），p为多年平均降雨量（mm），α为平均地表径流系数，如表如下。
[0053]
各类型生态系统地表径流系数均值表
蒸散发因子：根据国家生态系统观测研究网络科技资源服务系统网站提供的产品数据。原始数据空间分辨率为 1km，通过 arcgis 软件重采样为 250m 空间分辨率，得到蒸散发因子栅格图。
[0054]
生态系统面积因子：根据全国生态状况遥感调查与评估成果中的生态系统类型数据集得到。原始数据为矢量数据，通过arcgis软件转为 250m 空间分辨率的栅格图。
[0055]
模型运算将各因子统一成 250m 分辨率的栅格数据，在arcgis栅格计算器（spatial analyst
→
raster calculator）中，根据公式计算得到生态系统水源涵养量。
[0056]
（7）生态敏感性生态敏感性包括了水土流失敏感性、石漠化敏感性和土地沙化敏感性。一般的，生态越敏感其生态系统恢复力越低，反之则越高。
[0057]
①
水土流失敏感性评价根据土壤侵蚀发生的动力条件，水土流失类型主要有水力侵蚀和风力侵蚀。以风力侵蚀为主带来的水土流失敏感性将在土地沙化敏感性中进行评价；本节主要对水动力为主的水土流失敏感性进行评价，根据原国家环保总局生态功能区划技术规范的要求，并结合研究区的实际情况，选取降水侵蚀力、土壤可蚀性、坡度坡长和地表植被覆盖等评价指标，并根据研究区的实际对分级评价标准作相应的调整。将反映各因素对水土流失敏感性的单因子分布图，用地理信息系统技术进行乘积运算，公式如下：式中：ssi为i空间单元水土流失敏感性指数，评价因子包括降雨侵蚀力（ri）、土壤
可蚀性（ki）、坡长坡度（lsi）、地表植被覆盖（ci）。不同评价因子对应的敏感性等级值见水土流失敏感性的评价指标及分级赋值表。
[0058]ri
—降水侵蚀力值可利用降水资料计算的中国100多个城市的r值，采用内插法，用地理信息系统绘制r值分布图。
[0059]
lsi—坡度坡长因子：对于大尺度的分析，坡度坡长因子ls是很难计算的。这里采用地形的起伏大小与土壤侵蚀敏感性的关系来估计。在评价中，可以应用地形起伏度，即地面一定距离范围内最大高差，作为区域土壤侵蚀评价的地形指标。然后用地理信息系统绘制区域土壤侵蚀对地形的敏感性分布图。
[0060]ki
—土壤质地因子：可用雷诺图表示。通过比较土壤质地雷诺图和k因子雷诺图，将土壤质地对土壤侵蚀敏感性的影响分为5级。根据土壤质地图，绘制土壤侵蚀对土壤的敏感性分布图。
[0061]ci
—覆盖因子：地表覆盖因子与潜在植被的分布关系密切。根据植被分布图的较高级的分类系统，将覆盖因子对土壤侵蚀敏感性的影响分为5级，并利用植被图绘制土壤侵蚀对植被的敏感性分布图。
[0062]
水土流失敏感性的评价指标及分级赋值在数据条件具备的条件下也可采用通用水土流失方程（usle）计算评价区土壤侵蚀量的空间分布值，根据土壤侵蚀量大小进行水土流失敏感性分级。
[0063]
土地沙化敏感性评价根据原国家环保总局生态功能区划技术规范的要求，并结合研究区的实际情况，选取干燥指数、起沙风天数、土壤质地、植被覆盖度等评价指标，并根据研究区的实际对分级评价标准作相应的调整。
[0064]
根据各指标敏感性分级标准及赋值（土地沙化敏感性评价指标及分级表），利用地理信息系统的空间分析功能，将各单因子敏感性影响分布图进行乘积运算，得到评价区的土地沙化敏感性等级分布图，公式如下：；式中：di为i评价区域土地沙化敏感性指数；ii、wi、ki、ci分别为评价区域干燥指数、起沙风天数、土壤质地和植被覆盖的敏感性等级值。
[0065]
土地沙化敏感性评价指标及分级ii
—干燥指数：表征一个地区干湿程度，反映了某地、某时水分的收入和支出状况。
[0066]
式中：∑10℃指日温≥10℃持续期间活动积温总和；r为同期降水量（mm）。
[0067]
wi—起沙风天数：风力强度是影响风对土壤颗粒搬运的重要因素。已有研究资料表明，砂质壤土、壤质砂土和固定风砂土的起动风速分别为6.0、6.6和5.1m/s，选用冬春季节大于6m/s起沙风天数这个指标来评价土地沙化敏感性。根据研究区各气象站点的气象数据，在地理信息系统中利用插值生成土地沙化对起沙风天数敏感性的单因素评价图。
[0068]ki
—土壤质地：不同粒度的土壤颗粒具有不同的抗蚀力，粘质土壤易形成团粒结构，抗蚀力增强；在粒径相同的条件下，沙质土壤的起沙速率大于壤质土壤的起沙速率；砾质结构的土壤和戈壁土壤的风蚀速率小于沙地土壤；基岩质土壤供沙率极低，受风蚀的影响不大。以土壤质地图为底图，在地理信息系统中得出土壤质地对土地沙化敏感性的单因素评价图。
[0069]ci
—植被覆盖度地表植被覆盖是影响沙化敏感性的一个重要因素，在水域、冰雪和植被覆盖高的地区，不会发生土壤的沙化；相反，地表裸露、植被稀少都会使土壤沙化的机会增加。因此，植被覆盖是评价土地沙化敏感性的又一重要指标。
[0070]
式中：ndvi
veg
为完全植被覆盖地表所贡献的信息，ndvi
soil
为无植被覆盖地表所贡献的信息。
[0071]
运用地理信息系统软件进行图像处理，获取植被ndvi影像图，进而计算植被覆盖度。由于大部分植被覆盖类型是不同植被类型的混合体，所以不能采用固定的ndvi
soil
和ndvi
veg
值，通常根据ndvi的频率统计表，计算ndvi的频率累积值，累积频率为2%的ndvi值为ndvi
soil
，累积频率为98%的ndvi值为ndvi
veg
。
[0072]
③
石漠化敏感性评价石漠化敏感性主要取决于是否为喀斯特地形、地形坡度、植被覆盖度等因子。根据各单因子的分级及赋值（见下表格），利用地理信息系统的空间叠加功能，将各单因子敏感性影响分布图进行乘积计算，得到石漠化敏感性等级分布图，公式如下：；式中：si为i评价区域石漠化敏感性指数；di、pi、ci为i评价区域碳酸盐出露面积百分比、地形坡度和植被覆盖度；di是区域单元范围内碳酸盐出露面积占单元总面积的百分
比；pi根据评价区数字高程（dem）在地理信息系统下进行处理和分级；ci的数据来源和处理方法参照土地沙化敏感性。
[0073]
土地石漠化敏感性评价指标及分级生态系统连通性图斑大小即景观面积(ta)，单位: ha。
[0074]
景观面积决定了景观的范围以及研究和分析的最大尺度，也是计算其它指标的基础。图斑面积越大对于维护高数量的物种，维持稀有种、濒危种以及生态系统的稳定具有重要作用。
[0075]
连片度代表了各类生态景观的连接程度，一般的，越连片的生态系统其抵抗力越强，生态系统越稳定，生态系统生物多样性维护能力越强。
[0076]
地块形状指数地块形状指数代表生态景观类型的聚集程度，一般生态景观的地块形状指数越大、地块形状越规则（越接近于圆形），其聚集程度越高，生态系统越稳定。
[0077]
结果：依据集成评价结果，四川省林地生态系统恢复力综合评分74.51分，整体表现为生态系统恢复力较高。其中，高值区分布于龙门山区，秦巴山区及大小凉山，低值区主要分布于城镇及其周边区域，川西干热河谷地带，及高山低植被覆盖区。
[0078]
根据分级标准，四川省全域处于中等和较高水平的恢复力所占比重最大，其中，恢复力达到强或较强水平以上的区域面积为141362.88 km2，占区域总面积的29.03%，主要分布在门山区，秦巴山区及大小凉山，该区域降雨丰富，植被覆盖度高，土壤肥力好，光热条件好，利于植被再生长；恢复力为一般水平的面积为42046.57km2，占比8.64%，主要分布在四川盆地内植被覆盖度较高的丘陵、山区，以及盆周山地人为活动较为频繁区域；恢复力为弱和较弱的区域面积为303418.89 km2，占比62.33%，主要分布在四川盆地内植被覆盖度低的丘陵区，人类活动频繁的城镇及其周边区域，川西高山、高寒区。具体评价结果见图（见林地生态系统恢复力分级表）。
[0079]
林地生态系统恢复力分级表对于本领域技术人员而言，显然本技术不限于上述示范性实施例的细节，而且在
不背离本技术的精神或基本特征的情况下，能够以其它的具体形式实现本技术。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本技术的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本技术内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。