还剩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经验公式计算Wischmeier1969;式中pi和p分别是月均和年均降雨量mm计算得到各站点在2000-2007年平均降雨侵蚀力,然后利用Kriging空间内插方法对34个站点包括@@@@@站点进行插值,得到流域水平降雨侵蚀力图层,最后得到流域30m×30m的R因子栅格图层(图2)K为土壤可蚀性因子(t·hm-2·h·MJ-1·mm-1·hm-2)它是衡量土壤抗蚀性的指标,用于反映土壤对侵蚀的敏感性K表示标准小区单位降雨侵蚀力引起的单位__上的土壤侵蚀量由于缺乏各土壤类型的结构系数和渗透性等级数据,因此选择侵蚀/生产力影响模型EPIC的公式计算流域各类型土壤的K因子值,EPIC的计算公式为式中SAN、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的栅格图层,利用ARCGIS
9.2中的RasterCalculator模块把K值赋给土壤类型,得到K因子图层(图3)LS为坡长坡度因子(无量纲),其中L为坡长因子,被定义为坡长的幂函数S为坡度因子,LS表示在其他条件不变的情况下,某给定坡长和坡度的坡面上土壤流失量与标准径流小区典型坡面上土壤流失量的比值,它对土壤侵蚀起加速作用(Renardetal.1991)在ARCGIS
9.2中利用流域的DEM数据提取出山脊线,然后利用Distan__模块计算每格到山脊线的垂直距离,以此作为每个栅格的近似坡长,然后采用Wischmeier和__ith的坡长因子计算公式(Wischmeierand__ith,1978)以上各式中,λ为水平坡长,α为坡长指数,
22.3为标准小区的坡长(m),θ为利用DEM提取的坡度得到流域30m×30m的L因子图层S因子采用McCool的计算公式以及刘宝元的陡坡计算公式(McCool,1987;LiuBY,1994)利用流域DEM数据,在ARCGIS
9.2中Spatial____ysis模块中的Su_________ysis功能提取坡度专题图,并转化为弧度单位,然后在RasterCalculator模块中利用公式计算得到流域30m×30m的S因子图层C——植被覆盖与作物管理因子C因子由人工与天然植被的覆盖度和植被类型决定由于平顶山地区尚无C值的资料库可以利用,而一般认为归一化植被指数比较真实地反映了植被的分布,因而选择求取的归一化植被分布指数(图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.3436l__式中C——植被盖度C=0,则等于c=
78.3%;C=1,则等于c=
0.1%表明当植被盖度
78.3%时,地表的侵蚀量都是极其微弱,其侵蚀量基本上可以忽略不计但植被盖度
0.1%时,它的减蚀作用基本没有反映因而上式可以表示为。