本次采用中国土壤流失方程CSLE计算土壤侵蚀模数,计算公式为:
A=RKLSBET
式中,式中:A—土壤侵蚀模数。t·hm-2a-1;R—降雨侵蚀力因子,MJ·mm·hm-2·h-1·a-1;K—土壤可蚀性因子,t·hm2·h·hm-2·MJ·mm-1;L—坡长因子,无量纲;S—坡度因子,无量纲;B—植被覆盖度与生物措施因子,无量纲;E—工程措施因子,无量纲;T—耕作措施因子,无量纲。
本次分享R因子计算过程,采用ArcGIS计算R因子首先需要降雨量数据,研究区降雨量插值计算详见第二篇文章。
01降雨侵蚀力因子R计算公式
R因子计算根据所获取的降雨资料分别采用不同计算公式:
(1)当获得逐日降雨资料时,按如下公式计算:
式中:
Rr—第r日的降雨侵蚀力因子,MJ·mm/(hm2·h);
Pr—第r日的日降雨量,mm,需单日降雨量≥12mm,否则按0计;
α,β—计算降雨侵蚀因子统计系数,无量纲;
Pd12—日降雨量≥12mm的日降雨量,mm;
Py12—日降雨量≥12mm的年降雨量,mm;
获得逐日降雨侵蚀力因子后,可累加得多日、月、年降雨侵蚀力因子。
(2)当获得月降雨资料时,按如下公式计算:
式中:
Rm—第m月的降雨侵蚀力因子,MJ·mm/(hm2·h);
Pm—第m月的降雨量,mm。
获得月降雨侵蚀力因子后,可累加得多月、年降雨侵蚀力因子。
(3)可获得年降雨资料时,按如下公式计算:
式中:
Rn—年降雨侵蚀力因子,MJ·mm/(hm2·h);
Pn—年降雨量,mm。
(4)可获得多年平均降雨资料时,按如下公式计算:
式中:
Rd—多年平均降雨侵蚀力因子,MJ·mm/(hm2·h);
Pd—多年平均降雨量,mm。
由中国气象数据网下载数据为累年月值降雨量,可分别计算12个月的R因子和多年平均R因子,本次只计算多年平均R因子做示范,按照公式(4)计算
02降雨侵蚀力因子计算
打开插值后的降雨量数据,打开Spatial Analyst Tools—地图代数—栅格计算器工具
输入计算公式
点击确定,即可计算得出内黄县R因子图层。
好了,降雨侵蚀力因子计算完成,R因子计算很简单,重点是降雨量数据的插值计算,需要注意的是不同的降雨资料选用不同的计算公式,如果学习过程中缺少行政区边界数据,可留言需要的行政区,后续我也会分享如何根据影像数据自己制作区域边界。