土壤侵蚀不仅会引起耕地生产力下降、河床抬升、泥沙淤积阻塞河道等生态环境问题,也会对人们正常的生产生活产生威胁。生态系统的土壤保持量(吨/公顷/年),是土壤潜在侵蚀量与实际侵蚀量的差值。
地理遥感生态网平台提供的生态系统服务——土壤保持量分布数据,采用通用土壤流失方程模拟。在具体计算时,需利用已有实测土壤侵蚀数据对模型模拟结果进行验证,并修正参数(式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来源于文献或专家咨询。