土壤侵蚀 的估算方法作者：牛健平时间：2011年10月11日北京天合数维科技有限公司(CONTENT)一、 所需数据与参数 31、 所需数据 32、 所需中间参数 32.1、 水土保持因子P 32.2、 地标覆盖因子C 32.3、 地形因子LS 42.4、 土壤可视性因子K 42.5、 降水侵蚀因子R 43、 所需参数 53.1、 潜在土壤侵蚀量Ap 53.2、 现实土壤侵蚀量Ar 53.3、 土壤保持量Ac 54、 指标结果参数 54.1、 保护土壤肥力的经济效益Ef 64.2、 减少土地废弃的经济效益Es 64.3、 减轻泥沙淤积的经济效益En 6二、 处理流程 71、 DEM数据的处理 81.1、 坡长 L 81.2、 百分比坡度a 81.3、 地形因子LS 92、 气象数据 92.1、 月降雨量Pi的计算 92.2、 土壤侵蚀力指标R 103、 土壤类型数据 104、 遥感影像数据 105、 土壤理性化数据 11三、 所需参数的计算 11四、 指标结果参数计算 11、所需数据与参数在计算的过程中，总共涉及到的数据有地形数据、遥感影像数据、气象数据、 土壤类型数据、土壤理性化数据以及统计数据，涉及到的中间参数有水土保持因 子P，地标覆盖因子C，地形因子LS，土壤可视性因子K，降水侵蚀因子R，所 需要的参数有潜在土壤侵蚀量Ap，现实土壤侵蚀量Ar，土壤保持量Ac，指标结 果参数有保护土壤肥力的经济效益Ef，减少土地废弃的经济效益Es，减轻泥沙 淤积的经济效益En。
1、 所需数据在进行土壤侵蚀的估算过程中，需要以下数据：A、 地形数据；B、 遥感影像数据；C、 气象数据，主要是降雨量数据；D、 土壤类型数据；E、 土壤理性化数据；F、 统计数据2、 所需中间参数在数据处理的过程中，所涉及到的中间参数与计算公式如下2.1、 水土保持因子P按照游松财的方法，水田的P值取0.15,其他土地利用方式基本没有采取水 土保持措施，因此取值为1.002.2、 地标覆盖因子C地表覆盖因子是根据地面植被覆盖状况不同而反映植被对土壤侵蚀影响的因素，与土地利用类型、覆盖度密切相关C值的估算采用如下公式:c — 1, fc = 0c = lc = 0.6508-0.3436lg fC,0 < fC <78.3%c = 0, fC > 78.3%式中，fC表示植被覆盖度，由遥感影像通过NDVI计算2.3、 地形因子LS通过数字高程模型(DEM)，计算获得坡长和坡度，然后根据黄炎和等建 立的方程式，获得LS的空间分布特征LS = 0.08 L0.35 a 0.6式中，L，坡长(m)； a，百分比坡度2.4、 土壤可视性因子K采用陈明华等人建立的土壤可蚀性K值的计算公式：K = 10- 3 (160. 80 - 2. 31X 1 + 0. 38X 2 + 2. 26X 3 + 1. 31X 4 + 14. 67X5)式中：K 一可蚀性K值(美国习用单位)；X 1、X 2、X 3、X 4、X 5分别表示 细砾(1 〜3mm )%、细砂(0. 05〜0. 25mm )%、粗粉粒(0. 01 〜0. 05 mm )%、 细粉粒(0. 005〜0. 01mm )%、有机质(10g/kg)。
2.5、 降水侵蚀因子R采用周优建和黄炎和等人(1997,2000,2002)根据实测数据提出的R值计算 式，该值考虑了月降水量R =寸(-1.5527 + 0.1792Pi)i =1式中，R为降雨侵蚀力指标(j.cm/m2.h)，Pi为月降雨量(mm)月降水量P.： 根据气象站点观测的降水点位及数据，进行克里金插值获得3、 所需参数在计算的过程中主要涉及到潜在土壤侵蚀量Ap，现实土壤侵蚀量Ar以及土 壤保持量Ac三个参数，可以通过以上5个中间参数来获得，其计算公式如下3.1、 潜在土壤侵蚀量Ap潜在土壤侵蚀量不考虑地表覆盖类型和土地管理因素，即c=1，p=1，此时， USLE的形式为：A = R x K x LS式中，Ap为潜在土壤侵蚀量（吨/平方千米），R为降雨侵蚀力指标，K为土壤可 侵蚀因子，LS为坡长坡度因子3.2、 现实土壤侵蚀量Ar现实土壤侵蚀量主要考虑了地表覆盖类型和土地管理因素，其计算式为：A = R x K x LS x C x P式中，Ar为现实土壤侵蚀量（吨/平方千米），R为降雨侵蚀力指标，K为土壤可 侵蚀因子，LS为坡长坡度因子，C为地表覆盖因子，P为土壤保持措施因子。
3.3、 土壤保持量Ac由上两式可得到土壤保持量：平式中，Ap为潜在土壤侵蚀量（吨/平方千米），Ar为现实土壤侵蚀量（吨/平方千米）， Ac为土壤保持量（吨/平方千米）4、 指标结果参数指标结果参数总共涉及到的数据有保护土壤肥力的经济效益Ef、减少土地废 弃的经济效益Es和减轻泥沙淤积的经济效益En三个参数，其计算公式与所需参 数值见下4.1、保护土壤肥力的经济效益Ef土壤侵蚀使大量土壤营养物质流失，主要表现在N、P、K等营养物质的 流失Ef = （A xD xC xRxC） + Z A.C.F^/10000（' = N,P,K）i式中，E为保护土壤肥力的经济效益（元/ a）； A，土壤保持量（t/hm2a）； D， f c e土壤密度（t/m3）； C土壤有机质含量（kg/t）； R，薪柴转化为土壤有机质的 比例，一般为2； C，薪柴的平均价格（元/kg）； q，土壤中N、P、K的纯含量； Pi，N、P、K的价格注：N、P、K的价格：参照国家林业局2008年4月28日发布的《中华人民共和 国林业行业标准一森林生态系统服务功能评估规范》，氯化钾价格为2200元/吨， 磷酸二铵价格为2400元/吨。
元素名称价格（元/吨）N509.04P563.52K1151.7表1、NPK价格表4.2、 减少土地废弃的经济效益Es根据土壤保持量和土壤表土平均厚度（0.6m）来推算因土壤侵蚀而造成的废 弃土地面积，再根据机会成本法计算因土地废弃而失去的年经济价值Es=Ac!P!0.6XB!10000式中，Es，减少土地废弃的经济效益（元/a）； Ac，土壤保持量（t/a）； P，土壤 的容重（t/m3）； B，湿地年均收益（元/hm2）注：土壤容量取26.5吨/立方米，湿地年均收益根据崔丽娟（2004年）研究推 算，湿地年均收益为245.5元/平方米4.3、 减轻泥沙淤积的经济效益En按照我国主要流域的泥沙运动规律，全国土壤侵蚀流失的泥沙有24%淤积于 水库、江河、湖泊，这部分泥沙直接造成了水库江河、湖泊蓄水量的下降，在一定程度上增加了干旱、洪涝灾害发生的机会，因此可根据蓄水成本计算损失价值En二Ac!pX24%XC式中，En,减轻泥沙淤积的经济效益（元/a）； Ac, 土壤保持量（t/a）； P, 土壤 的容重（t/m3）； C,水库工程费用（元/m3）注：土壤容量取26.5吨/立方米，核算过程中采用的水库清淤工程费用为6.94 元/立方米，此数据来源于硕士论文。
二、处理流程总体处理流程如下图所示:潜在土壤侵蚀量 （AP） NPK价格Pi NPK纯含量Ci 土壤容量P 土壤密度De 土壤有机质含量Co 湿地年均收入B 水库工程费用Z土壤保持量（）图1、数据处理流程图1、DEM数据的处理运用DEM数据得到的最后结果是求出了地形因子LS,在求LS的过程中，需 要坡长L和百分比坡度a，其求解主要在Arcmap中完成1.1、坡长L坡长的求解主要参考汤国安老师的ArcGIS地理信息系统空间分析实验教程 中的第十一章中的水文分析来完成，Arcmap中建立模型如下图所示：图2、无洼地DEM生成模型图3、生成坡长L模型1.2、百分比坡度a百分比坡度a是用slope来实现的图4、百分比坡度的生成1.3、地形因子LS地形因子LS计算的模型为:图5、地形因子LS的计算模型2、气象数据根据气象站点观测的降水点位及数据，进行克里金插值获得月降雨量数据 Pi，再根据Pi计算土壤侵蚀力指标R2.1、月降雨量Pi的计算2.2、土壤侵蚀力指标R土壤侵蚀力指标 P 是在 Arcmap 下 Spatial Analyst|Raster Caculator…中 计算而来图6、栅格计算3、土壤类型数据土壤可蚀性因子K是通过土壤类型数据直接对其进行赋值。
4、遥感影像数据通过遥感影像数据获得地表覆盖因子C，首先求NDVI，NDVI的求解在ENVI 中实现，主菜单中transform—ndvi，算归一化植被指数接着根据NDVI数据 来求解植被覆盖度FC,其计算公式为：FC=(NDVI-NDVI_MIN)/(NDVI_MAX-NDVI_MIN)式中：NDVI_MIN为NDVI的最小值，NDVI_MAX为NDVI的最大值由此，我们可以根据公式c=lfC二°c 二* 二 °.65°&°.3436lgG° < fC< 78.3%c 二 °,f^C〉78.3%来计算 C 值，在 Arcmap 中 Spatial Analyst|Raster Caculator…实现C = con([fc.img] == 0, 1, _([fc.img] > 0 & [fc.img] <1.?83,0.6508-0.3^36 * LoglOftfc.img])),{0})5、土壤理性化数据按照游松财的方法，水田的P值取0.15，其他土地利用方式基本没有采取水 土保持措施，因此取值为1.00三、所需参数的计算所需要的参数有潜在土壤侵蚀量Ap，现实土壤侵蚀量Ar，土壤保持量Ac， 在进行了上述的计算以后，其计算比较简单，根据公式，在SpatialAnalyst|Raster Caculator…中输入对应公示即可实现。
最后对Ac进行汇总统计，在Spatial Analyst|Zonal statistics中进行Spatial Analyst T Layer：插~DistanceDensity...Interpolate to RasterSurface- AnalysisCell Statistics...Neighborhood Statistics...Zonal Statistics...四、指标结果参数计算这一步只是简单的数值计算，在此不再详述，但是在计算的过程中注意单位的统一注：在进行统计Ac总和的时候，要乘以图像分辨率还有单位一定要统一！。