土壤流失模型(RUSLE)土壤流失评价.docx

《土壤流失模型(RUSLE)土壤流失评价.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、9.4土壤流失模型（RUSLE）：土壤流失评价水土流失(waterLossandSoilErosion)是当今世界面临的一个重大环境问题，是指在水力、重力、风力等外力的作用下，水土资源和土地生产力的破坏和损失，包括了土壤侵蚀(SoilErosion)及水的流失两个方面。水土流失是一个受多种因素驱动作用的自然过程，由降雨作为主要侵蚀动力的即为水力侵蚀，其作用多体现于具有一定坡度的山地、丘陵地区，以降水直接冲走表层土壤为主要表现形式。当表层土壤受侵蚀的程度超过一定的安全阈值时，就会形成或引发滑坡、泥石流等剧烈的土壤流失过程，降低

2、区域土壤生产力，严重影响生态系统的平衡。由于这一过程表现出明显的斑块性和异质性，难以用定点实验进行追踪监测，故多以预测估算的方式对之进行研究。9.4.1RUSLE模型RUSLE模型(TheRevisedUniversalSoilLossEquation)是美国农业部于1997年在通用土壤流失模型USLE(TheuniversalSoilLossEquation)的基础上修订建立并正式实施的一种适用范围更广的修正模型。自颁布之后即在美国得到了广泛的应用，目前美国各州都己有成功应用的案例。该模型也被世界各国的研究者借鉴，于20世

3、纪90年代被引入中国。RUSLE模型目前己在国内外的土壤侵蚀预测研究中得到了非常广泛的应用。其计算表达式为：(1)图1RUSLE模型的技术流程图A为预测土壤侵蚀量(estimatedaveragesoillossintonsperacreperyear)，主要指由降雨和径流引起的坡面细沟或细沟间侵蚀的年均土壤流失量。R为降雨侵蚀力因子(rainfall-runofferosivityfactor)(MJ·mm·hm-2·h-1·a-1）,它反映降雨引起土壤流失的潜在能力。在USLE中，它被定义为降雨动能和最大30分钟降雨强度

4、的乘积。K为土壤可蚀性因子(soilerodibilityfactor)(t·hm2·h·MJ-1·mm-1·hm-2)，它是衡量土壤抗蚀性的指标，用于反映土壤对侵蚀的敏感性。K表示标准小区单位降雨侵蚀力引起的单位面积上的土壤侵蚀量。LS为坡长坡度因子(L=slopelengthfactor,s=slopesteepnessfactor)(无量纲)，其中L为坡长因子，被定义为坡长的幂函数。S为坡度因子，LS表示在其他条件不变的情况下，某给定坡长和坡度的坡面上土壤流失量与标准径流小区典型坡面上土壤流失量的比值，它对土壤侵蚀起加

5、速作用。C为覆盖与管理因子(cover-managementfactor)(无量纲)，它指在其他因子相同的条件下，在某一特定作物或植被覆盖下的土壤流失量与耕种后的连续休闲地的流失量的比值。该因子衡量植被覆盖和经营管理对土壤侵蚀的抑制作用。P为水土保持措施因子(supportpracticefactor)(无量纲)，它指采取水土保持措施后的土壤流失量与顺坡种植的土壤流失量的比值。9.4.2RUSLE模型所需数据搜集进行RUSLE模型的计算需要不同的数据，对所需数据的描述如下：1)研究区大比例尺地形图地形图；2)研究区降雨量监测

6、数据；3)研究区土壤类型图；4)研究区土壤普查资料（包括土壤机械组成、粒级含量、有机质含量等数据）；5)研究区数字高程数据(DEM)；6)研究区不同时段归一化植被指数NDVI数据；7)研究区不同时段土地利用栅格数据9.4.3RUSLE模型中因子的计算9.4.3.1R因子的计算R因子的计算一般通过基于月平均降雨量和年平均降雨量的Wischmeier公式计算[1]：(2)其中pi和p分别是月均和年均降雨量(mm)。得到R因子的栅格图如图2所示。图2R因子的栅格图9.4.3.2K因子的计算K的估算公式根据不同的土壤组成与搜集数据的

7、多少而变化。如根据美国粒径分级制进行计算的Whischmeier公式[2]：(3)其中，(据美国粒径分级制)；OM=土壤有机质含量；S=结构系数；P=渗透性等级。得到研究区K因子的栅格图如图3所示。图3K因子的栅格图9.4.3.3LS因子的计算在ARCGIS中利用流域的DEM数据提取出山脊线，然后利用Distance模块计算每个栅格到山脊线的垂直距离，以此作为每个栅格的近似坡长，然后采用Wischmeier和Smith提出的坡长因子计算公式[3]：(4)(5)(6)以上各式中，λ为水平坡长，α为坡长指数，22.13为标准小区

8、的坡长(m)，θ为利用DEM提取的坡度。得到研究区的L因子的栅格图如图4所示。S因子采用McCool的计算公式以及刘宝元的陡坡计算公式[4][5]：(7)利用研究区DEM数据，在ARCGIS中SpatialAnalyst模块中的SurfaceAnalysis功能提取坡度专题图，并转化为弧度