文档介绍：附件
生态保护红线划定指南
环境保护部国家发展改革委
2017年5月
生态系统服务功能重要性评估方法目前,生态系统服务功能采用得评估方法主要有模型评估法与
净初级生产力（NPP）定量指标评估法。其中,模型评估法所需参数较
多2MJmm；L)、S为地形因子,L表示坡长因子,S表示坡度因
子；C为植被覆盖因子。
A、1、2、2数据准备
数据来源与获取根据上述模型,水土保持功能重要性评估需
收集高程数据集、
气象数据集与土壤数据集等数据,具体信息见表A3。
表A3水土保持功能重要性评估数据表
名称
类型
分辨率
数据来源
高程数据集
栅格
30m
地理空间数据云网站
气象数据集
文本
一
中国气象科学数据共享服务网
土壤数据集
矢量/Excel
一
全国生态环境调查数据库中
国1:100万土壤数据库
(2)数据预处理
降雨侵蚀力因子R就是指降雨引发土壤侵蚀得潜在能力,通过多年平均年降雨侵蚀力因子反映，计算公式如下：
24
RR半月k
k1
式中,R为多年平均年降雨侵蚀力(MJmm/h2rha)R半月卜为第k个半月得降雨侵蚀力(MJmm/hrmha)k为一年得24个半月，k=1,2,…,24;i为所用降雨资料得年份,i=1,2,…,n;j为第i年第k个半月侵蚀性降雨日得天数，j=1,2,…,m;Pk为第i年第k个半月第j个侵蚀性日降雨量(mmX以根据全国范围内气象站点多年得逐
日降雨量资料，通过插值获得；或者直接采用国家气象局得逐日降
雨量数据产品。a为参数,暖季时a=0、3937,冷季时a=0、
3101。
土壤可蚀性因子K指土壤颗粒被水力分离与搬运得难易程度,
主要与土壤质地、有机质含量、土体结构、渗透性等土壤理化性质有关，计算公式如下：
K0、013830、51575KEPIC0、1317
K0、20、3exp0、0256ms(1msilt/100)/(mm)]0、3
[mEPICsiltcsilt
10、25orgC/[orgCexp(3、722、95orgC)]
10、7(1ms/100)/{(1ms/100)exp[5、5122、9(1ms/100)]}
式中,KC表示修正前得土壤可蚀性因子，K表示修正后得土壤可蚀性因子,mm*0C分别为粘粒(<0、002mm淞粒(0、002mnr0、05mm)砂粒(0、05mm―2mm汽有机碳得百分比含量(％)数据来源于中国1:100万土壤数据库。在Excel表格中,利用上述公式计算
K值，然后以土壤类型图为工作底图，在ArcGIS中将K值连接(Join)到底图上。利用ConversionTools中矢量转栅格工具,转换成空
间分辨率为250m得土壤可蚀性因子栅格图。
地形因子LSL表示坡长因子，S表示坡度因子，就是反映地形对
土壤侵蚀影响得两个因子。在评估中,可以应用地形起伏度,即地面
一定距离范围内最大高差,作为区域土壤侵蚀评估得地形指标。选
择高程数据集,在SpatialAnalyst下使用NeighborhoodStatistics,设置StatisticType为最大值与最小值,即得到高程数
据集得最大值与最小值,然后在SpatialAnalyst下使用栅格计算
器RasterCalculator,公式为[最大值最小值],获取地形起伏度,即
地形因子栅格图。
植被覆盖因子C:反映了生态系统对土壤侵蚀得影响,就是控制土
壤侵蚀得积极因素。水田、湿地、城镇与荒漠参照NSPECT导参数
分别赋值为0、0、0、01与0、7,旱地按植被覆盖度换算,计算公式
如下:
C旱=0、221059510gc
式中，C旱为旱地得植被覆盖因子,C1为小数形式得植被覆盖度。
其余生态系统类型按不同植被覆盖度进行赋值,如表A4所示
表A4不同生态系统类型植被覆盖因子赋值
生态系统类型
植被覆盖度
<10
1030
3050
5070
7090
>90
森林
0、
0、08
0、06
0、02
0、004
0、001
灌丛
01
0、22
0、14
0、085
0、04
0、011
草地
0、445
0、24
0、15
0、09
0、043
0、011
乔木园地
0、42
0、23
0、14
0、089
0、042
0、011
灌木园地
0、
0、22
0、14
0、087
0、042
0、011
A、1、2、3模地运算
将各因子统一成250m分辨率得栅格数据,在ArcGIS栅格计算
器（SpatialAnalyst—RasterCalculator）中，根据公式计算得到生态系统水土保