土壤侵蚀不仅会引起耕地生产力下降、河床抬升、泥沙淤积阻塞河道等生态环境问题，也会对人们正常的生产生活产生威胁。生态系统的土壤保持量（吨/公顷/年），是土壤潜在侵蚀量与实际侵蚀量的差值。
地理遥感生态网平台提供的生态系统服务——土壤保持量分布数据，采用通用土壤流失方程模拟。在具体计算时，需利用已有实测土壤侵蚀数据对模型模拟结果进行验证，并修正参数（式2–4）。
土壤保持量：
Asoilerosion=Apotentialerosion−AacturalerosionAsoilerosion=Apotentialerosion−Aacturalerosion （2）
实际土壤侵蚀量：
Aacturalerosion=R×K×L×S×CAacturalerosion=R×K×L×S×C （3）
潜在土壤侵蚀量：
Apotentialerosion=R×K×L×SApotentialerosion=R×K×L×S （4）
式中Aactual erosion 为单位面积实际土壤侵蚀量，单位为t/(hm2·a)；Apotential erosion 为单位面积潜在土壤侵蚀量，单位为t/(hm2·a)；Asoil conservation 为单位面积土壤保持量，单位为t/(hm2·a)；R为降雨侵蚀力因子，用多年平均年降雨侵蚀力指数表示；K为土壤可蚀性因子，表示为标准小区下单位降雨侵蚀力形成的单位面积上的土壤流失量。L为坡长因子（无量纲）；S为坡度因子（无量纲）；C为植被覆盖因子（无量纲）。其中，降雨侵蚀力因子R来源于Fouriner指数（式5）：
（5）
式中R为降雨侵蚀力因子，单位为（MJ mm hm-2 h-1 a-1）；j为月降水量（mm）；J为年降水量（mm）；i为月份。
土壤可蚀性因子K的计算公式为：
（6）
式中K为土壤可蚀性因子，表示为标准小区下单位降雨侵蚀力形成的单位面积上的土壤流失量。mc 、msilt 、ms 和orgC分别为粘粒（＜0.002mm）、粉粒（0.002–0.05mm）、砂粒（0.05–2mm）和有机碳的百分含量（％）。
坡长–坡度因子通过下述模型从DEM中提取。坡长因子L、坡度因子S计算方法如下（式7–10）：
（7）
（8）
（9）
（10）
式中L为坡长因子（无量纲），m为坡长指数，λ为坡长（m）；S为坡度因子（无量纲），θ为坡度（°）；植被覆盖管理因子C根据植被覆盖类型和植被覆盖度查询文献资料或咨询专家；水土保持措施因子P来源于文献或专家咨询。
原文链接：https://bbs.csdn.net/forums/gisrs?spm=1001.2014.3001.6682