水土流失计算方法

1、RUSLE 模型是通过对通用土壤流失方程 USLE 模型的改进得到的。RUSLE 与 USLE 具有相同的数学表达式：A=R·K·LS·C·P式中，A 为年均土壤侵蚀量（t·hm-2·a-1）,主要指由降雨和径流引起的坡面细沟或细沟间侵蚀的年均土壤流失量；R 为降雨侵蚀力因子（MJ·mm·hm-2·h-1·a-1）,它反映降雨引起土壤流失的潜在能力。本方案基于月平均降雨量和年平均降雨量的Wischmeier经验公式计算(Wischmeier, 1969)；式中pi和p分别是月均和年均降雨量(mm

2、)。计算得到各站点在2000-2007年平均降雨侵蚀力，然后利用Kriging空间内插方法对34个站点(包括站点)进行插值，得到流域水平降雨侵蚀力图层，最后得到流域30 m×30 m的R因子栅格图层（图2）。K为土壤可蚀性因子（t·hm-2·h·MJ-1·mm-1·hm-2）,它是衡量土壤抗蚀性的指标，用于反映土壤对侵蚀的敏感性。K表示标准小区单位降雨侵蚀力引起的单位面积上的土壤侵蚀量。由于缺乏各土壤类型的结构系数和渗透性等级数据，因此选择侵蚀/生产力影响模型EPIC 的公式计算流域各类型土壤的K因子值，EPIC的计算公式为：式中,S

3、AN 、SIL、CLA和C是砂粒、粉粒、粘粒和有机碳含量(%)，其中SN1=1-SAN/100。由公式II计算得到流域各土壤类型的K值如表2所示。表2 流域各土壤类型K因子值(单位:t·hm-2·h·MJ-1·mm-1·hm-2)土壤类型棕壤褐土石灰性褐土粗骨土红粘土草甸风沙土石灰土潮土红壤K值0.290.280.260.250.300.260.280.230.26将流域土壤类型图数字化，然后生成30m×30m的栅格图层，利用ARCGIS9.2中的Raster Calculator 模块把K值赋给土壤类型，得到K因子图层（图3）。LS

4、为坡长坡度因子（无量纲），其中L为坡长因子，被定义为坡长的幂函数。S为坡度因子，LS表示在其他条件不变的情况下，某给定坡长和坡度的坡面上土壤流失量与标准径流小区典型坡面上土壤流失量的比值，它对土壤侵蚀起加速作用（Renard, et al. 1991）。在 ARCGIS9.2 中利用流域的DEM数据提取出山脊线，然后利用Distance模块计算每格到山脊线的垂直距离，以此作为每个栅格的近似坡长，然后采用Wischmeier和Smith的坡长因子计算公式（Wischmeier and Smith，1978）：2 / 4以上各式中，为水平坡长，为坡长指数，22.3为标准小区的坡长（m），为利用DE

5、M提取的坡度,得到流域30m×30m的L因子图层。S因子采用McCool的计算公式以及刘宝元的陡坡计算公式（McCool，1987；Liu B Y，1994）。 利用流域DEM数据，在ARCGIS9.2中Spatial Analysis 模块中的Surface Analysis功能提取坡度专题图，并转化为弧度单位，然后在Raster Calculator模块中利用公式计算得到流域 30m×30m的S因子图层。 C植被覆盖与作物管理因子。C因子由人工与天然植被的覆盖度和植被类型决定。由于平顶山地区尚无C值的资料库可以利用，而一般认为归一化植被指数比较真实地反映了植被的分布，因

6、而选择求取的归一化植被分布指数（图1）与植被盖度图（图6）进行对比统计分析。二者的相关性达到0.77，二者线性回归方程为：Vc = 108.49Ic + 0.717式中：VC植被覆盖度（%）；IC归一化植被指数值。植被盖度求得之后，下一步将对实验区具体计算RUSLE模型中植被因子。由RUSLE模型对K值的定义，在标准状况下Y标准=100K，100为t/mh2到g/m2单位转换值，植被因子值为：C=Y/Y 标准，K若取平均值，则：C = 0.6508 - 0.3436 lgc式中：C植被盖度。C=0，则等于c=78.3%；C=1，则等于c=0.1%。表明当植被盖度> 78.3%时，地表的侵蚀量都是极其微弱，其侵蚀量基本上可以忽