生态系统服务—土壤侵蚀强度空间分布/降雨侵蚀力

      土壤侵蚀不仅会引起耕地生产力下降、河床抬升、泥沙淤积阻塞河道等生态环境问题，也会对人们正常的生产生活产生威胁。生态系统的土壤保持量（吨/公顷/年），是土壤潜在侵蚀量与实际侵蚀量的差值。

       地理遥感生态网平台提供的生态系统服务——土壤保持量分布数据，采用通用土壤流失方程模拟。在具体计算时，需利用已有实测土壤侵蚀数据对模型模拟结果进行验证，并修正参数（式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来源于文献或专家咨询。