本次采用中国土壤流失方程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—耕作措施因子，无量纲。
地形（坡度坡长）因子是土壤侵蚀模型中重要的参数，地形因子在DEM基础上生成，坡度因子计算现在已经比较成熟，计算方法采用如下公式：
操作步骤如下：
在ArcMAP中加载DEM数据（数据获取见上一篇文章：免费下载DEM数据并镶嵌拼接），选中系统工具箱—Spatial Analyst Tools—表面分析—坡度工具。
输入DEM数据，Z因子取值参考下图，点击确定。
结果如下图所示
选中系统工具箱—Spatial Analyst Tools—地图代数—栅格计算器，输入公式：
Con("Slope_NH" < 5,10.8 * Sin("Slope_NH") 0.03,Con("Slope_NH" < 10,16.8 * Sin("Slope_NH")-0.5,21.9 * Sin("Slope_NH")-0.96))
结果如下
随后将计算数据导出即可。
,
