式中：
?——多年平均年降雨侵蚀力（MJ mm hm-2 h-1 a-1）??；
-2-1-1???——第k个半月的降雨侵蚀力（MJ mm hm h a）； 半月??
k——一年的24个半月，即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）反映了土壤颗粒被水力分离和搬运的难易程度，是评价土壤对侵蚀敏感程度的重要指标，主要与土壤质地、有机质含量、土体结构、渗透性等土壤理化性质有关。采用如下公式进行计算。
KEPIC??0.2?0.3exp?-0.0256ms(1-msilt/100)???[msilt/(mc?msilt)]0.3 ??1?0.25orgC/[orgC?exp(3.72?2.95orgC)]? ??1?0.7(1?ms/100)/{(1?ms/100)?exp[?5.51?22.9(1?ms/100)]}?
K??-0.01383?0.51575KEPIC??0.1317式中：
K——土壤可蚀性因子（t hm2 h hm-2 MJ-1 mm-1）； mc——粘粒（<0.002 mm）百分含量（%）； msilt——粉粒（0.002~0.05 mm）百分含量（%）； ms——砂粒（0.05~2 mm）百分含量（%）； orgC——有机碳的百分含量（%）。
坡长和坡度因子（L、S）反映了坡长、坡度等对土壤侵蚀的影响，按照如下公式计算。 L=(
) 22.13
??
??
m=???(1+??) β=(sin???0.089)/[3.0×(sin??)0.8+0.56]
?10.8sin??0.03 ??5.14??16.8sin??0.5 5.14????10.20??S???21.91sin??0.96 10.20????28.81???9.5988 ??28.81?
式中：
L——坡长因子； S——坡度因子； m——坡长指数； θ——坡度（°）；
5
?——坡长（m）。
植被覆盖因子（C）反映了生态系统对土壤侵蚀的影响，是控制土壤侵蚀的积极因素。水田、湿地、城镇和荒漠分别赋值为0，0，0.01和0.7，其余各生态系统类型按不同植被覆盖度进行赋值（表A.1）。旱地的植被覆盖因子按照以下公式计算：
C=0.221?0.595logc
式中，C为旱地的植被覆盖因子，c为小数形式的植被覆盖度。
表A.1 不同植被覆盖的C值
植被覆盖度（%） 生态系统类型 <10 森林 灌丛 草地 乔木园地 灌木园地 0.10 0.40 0.45 0.42 0.40 10~30 0.08 0.22 0.24 0.23 0.22 30~50 0.06 0.14 0.15 0.14 0.14 50~70 0.02 0.085 0.09 0.089 0.087 70~90 0.004 0.040 0.043 0.042 0.042 >90 0.001 0.011 0.011 0.011 0.011 A.3 防风固沙量
采用修正风蚀方程RWEQ进行评价。 防风固沙量：
SR=SL潜-SL
潜在风蚀量：
22?z??z/s潜?SL潜=2QMAX潜?eS潜QMAX潜=109.8?WF?EF?SCF?K'? S潜=150.71?WF?EF?SCF?K'?实际风蚀量：
?0.37112?z?(z/s)2Q?eMAXS2S?150.71?(WF?EF?SCF?K'?C)?0.3711SL?Qmax?109.8?WF?EF?SCF?K'?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为土壤可蚀因子；SCF为土壤结皮因子；C为植被覆盖因子。
气象因子WF计算方法如下：
WF?Wf?式中：
6
?g?SW?SD
WF——气候因子，单位为kg/m，12个月WF总和得到多年年均WF； Wf——各月多年平均风力因子，ρ为空气密度，g为重力加速度； SW——各月多年平均土壤湿度因子，无量纲； SD——雪盖因子，无量纲。
土壤可蚀性因子EF的计算方法如下：
EF?29.09?0.31sa?0.17si?0.33?sa/cl??2.59OM?0.95Caco3 100
式中：sa为土壤粗砂含量（0.2 mm～2 mm）（%）； si——土壤粉砂含量（%）； cl——土壤粘粒含量（%）； OM——土壤有机质含量（%）；
Caco3——碳酸钙含量（%），可不予考虑。 土壤结皮因子SCF的计算方法如下：
11?0.0066?cl??0.021?OM2SCF?式中：
cl——土壤粘粒含量； OM——有机质含量。
?2
植被覆盖因子C的计算方法如下：
C?eai(SC)
式中：
SC——植被覆盖度； ai——不同植被类型的系数。 地表粗糙度因子K’的计算方法如下：
0 . 934 ? ? ( 1 . 86 K 2 . 41 K 0 . 127 Crr ) r r
? e K '
??H?Kr?0.2?
L式中：
Kr——土垄糙度，以Smith-Carson方程加以计算，单位cm； Crr——随机糙度因子，取0，单位cm； L——地势起伏参数；
ΔH——距离L范围内的海拔高程差。 A.4 生境不可替代性指数
7
2
首先选择指示物种，根据历史数据，以县为单元确定每个指示物种的分布区。使用Marxan选址运算模型，按约束条件进行迭代计算得到全国生物多样性保护优先区域。
评估单元为评估区内的各县级行政单位，物种选择全国境内有记录分布的国家一级物种、二级物种和其他有重要保护价值的物种，参照中国动物志和植物志统计这些物种在每个评估单元中的出现数量。利用不可替代性指数为评估单元赋值，该数值在0-100分布，数值越高表示该单元对保护生物多样性的价值（不可替代性）越大。
Marxan模型的迭代运算目标函数为：
∑??????????????+??????∑??????????????????????+∑??????????????????????+??????????????????????????????????????(??) ① ② 式中：
①——规划单元总成本；
②——保护体系边界总长度修正值； ③——未达到保护目标的补偿值； ④——超出成本阈值的补偿值。
③ 8
④
生态系统服务功能评估 - 图文
