生态系统服务功能指标计算方法

1、附录 A（规范性附录）生态系统服务功能指标计算方法A.1 水源涵养量计算生态系统通过拦截滞蓄降水，增强土壤下渗、蓄积， 涵养土壤水分、调节地表径流和补充地下水所增加的水资源总量。通过水量平衡方程计算：?=?()-3?=1?×?- ?- ? ×10式中：? 水源涵养量（m3/a）；? 产流降雨量（ mm/a）；?Ri 地表径流量（mm/a）；ETi 蒸散发量（ mm/a）；Ai i 类生态系统的面积（m2）；i 第 i 类生态系统类型；n 生态系统类型总数。A.2土壤保持量基于通用土壤流失方程（USLE ）计算。?= ?×?×?×?×(

2、 1 - ?)?式中：? 土壤保持量，单位为t/a；?R 降雨侵蚀力因子；K 土壤可蚀性因子；L 坡长因子（无量纲） ；S 坡度因子（无量纲） ；C 植被覆盖因子（无量纲）；降雨侵蚀力因子 R、土壤可蚀性因子 K 、坡长坡度因子 L 、S 及植被覆盖因子 C 的计算方法如下。降雨侵蚀力因子（R）是降雨引发土壤侵蚀的潜在能力，计算公式如下。?=24?=1半月 ?1?1.7265?半月 ?=?=1 ?=0(?,?,?)4式中：?MJ mm hm-2h-1 -1）；a 多年平均年降雨侵蚀力（? 第 k 个半月的降雨侵蚀力（MJ mm hm -2h-1 a-1）；半月 ?k 一年的 24 个半月，即

3、k 1,2,24；i 所用降雨资料的年份，即i 1,2,n；j 第 i 年第 k 个半月侵蚀性降雨日的天数，即j 1,2,m；Pi,j,k 第 i 年第 k 个半月第 j 个侵蚀性日降雨量（mm ）；为参数，暖季 0.3937，冷季 0.3101。可以根据全国范围内气象站点多年的逐日降雨量资料，通过插值获得。土壤可蚀性因子（K ）反映了土壤颗粒被水力分离和搬运的难易程度，是评价土壤对侵蚀敏感程度的重要指标，主要与土壤质地、有机质含量、土体结构、渗透性等土壤理化性质有关。采用如下公式进行计算。K EPIC 0.2ssilt0.30.3exp - 0.0256m (1- m /100) msilt

4、 / (mc msilt )10.25orgC / orgC exp(3.722.95orgC )10.7(1 ms /100) /(1ms /100)exp 5.51 22.9(1 ms /100)K-0.01383 0.51575 K EPIC0.1317式中：2-2-1-1K 土壤可蚀性因子（t hm h hmMJmm ）；msilt 粉粒（ 0.0020.05 mm ）百分含量（%）；ms 砂粒（ 0.052 mm ）百分含量（%）；orgC 有机碳的百分含量（%）。坡长和坡度因子 （ L 、S）反映了坡长、 坡度等对土壤侵蚀的影响，按照如下公式计算。?L = (22.13)?m =

5、? (1 + ?)0.8+ 0.56= (sin ?0.089 )/ 3.0 ×(sin ?)10.8 sin0.035.1416.8 sin0.55.1410.20S0.9610.2028.8121.91sin9.598828.81式中：L 坡长因子；S 坡度因子；m 坡长指数；坡度（ °）；5 坡长（ m）。植被覆盖因子（C）反映了生态系统对土壤侵蚀的影响，是控制土壤侵蚀的积极因素。水田、湿地、城镇和荒漠分别赋值为 0，0， 0.01 和 0.7，其余各生态系统类型按不同植被覆盖度进行赋值（表 A.1 ）。旱地的植被覆盖因子按照以下公式计算：C = 0.221 -0.5

6、95logc式中， C为旱地的植被覆盖因子，c为小数形式的植被覆盖度。表 A.1不同植被覆盖的 C值植被覆盖度（ % ）生态系统类型<101030305050707090>90森林0.100.080.060.020.0040.001灌丛0.400.220.140.0850.0400.011草地0.450.240.150.090.0430.011乔木园地0.420.230.140.0890.0420.011灌木园地0.400.220.140.0870.0420.011A.3防风固沙量采用修正风蚀方程RWEQ 进行评价。防风固沙量：SR=SL潜-SL潜在风蚀量：S潜 = 2 z QMA

7、X潜 e z/ s潜 2LS潜2QMAX 潜 =109.8 WFEFSCFK '0.3711S潜 =150.71 WFEFSCFK '实际风蚀量：SL2 2z QMAX e (z/ s)2SS 150.71 (WFEFSCFK' C)0.3711Qmax109.8 WFEFSCFK ' C式中， SR为固沙量（ t km -2 a-1）；SL 潜为潜在风力侵蚀量（ t km-2 a-1）；SL 为实际风力侵蚀量（t km-2 a-1）；QMAX 为最大转移量（ kg/m)；Z 为最大风蚀出现距离 （ m)；WF 为气候因子（kg/m）；K；为地表糙度因子； EF

8、 为土壤可蚀因子； SCF 为土壤结皮因子； C 为植被覆盖因子。气象因子 WF 计算方法如下：WFWfSWSDg式中：6WF 气候因子，单位为 kg/m， 12 个月 WF 总和得到多年年均WF ；Wf 各月多年平均风力因子，为空气密度， g 为重力加速度；SW 各月多年平均土壤湿度因子，无量纲；SD 雪盖因子，无量纲。土壤可蚀性因子 EF 的计算方法如下：29.09 0.31sa0.17 si 0.33 sa / cl 2.59OM 0.95Caco3EF100式中： sa 为土壤粗砂含量（0.2 mm 2 mm）（ %）；si 土壤粉砂含量（%）；cl 土壤粘粒含量（%）；OM 土壤有机

9、质含量（%）；Caco3 碳酸钙含量（%），可不予考虑。土壤结皮因子SCF 的计算方法如下：SCF10 . 0066 cl20. 021OM 21式中：cl 土壤粘粒含量；OM 有机质含量。植被覆盖因子C 的计算方法如下：Cea i ( SC)式中：SC 植被覆盖度；i 不同植被类型的系数。a地表粗糙度因子K的计算方法如下：K '(1 .86 K r2 .41 K r0 .9340 .127 Crr)e2K r0.2HL式中：Kr 土垄糙度，以 Smith-Carson 方程加以计算，单位cm；Crr 随机糙度因子，取 0，单位 cm；L 地势起伏参数；H距离 L 范围内的海拔高程差。A.4生境不可替代性指数7首先选择指示物种，根据历史数据，以县为单元确定每个指示物种的分布区。使用 Marxan 选址运算模型，按约束条件进行迭代计算得到全国生物多样性保护优先区域。评估单元为评估区内的各县级行政单位，物种选择全国境内有记录分布的国家一级物种、二级物种和其他有重要保护价值的物种，参照中国动物志和植物志统计这些物种在每个评估单