[发明专利]一种土壤有机碳侵蚀动态变化诱因分析方法

权利要求书

1.一种土壤有机碳侵蚀动态变化诱因分析方法，其特征在于，包括如下步骤：

S1，收集目标区域的地面降雨数据和土壤剖面数据，同时获取降雨栅格数据、植被栅格数据、土地利用类型栅格数据以及地形栅格数据；

S2，待收集完成后，结合多源时空数据挖掘及时空建模分析方法进行土壤侵蚀因子建模，得到各土壤侵蚀因子时空分布特征图；

S3，基于得到的所述土壤侵蚀因子时空分布特征图，计算得到土壤侵蚀空间分布图；

S4，基于所述的土壤侵蚀空间分布图，结合土壤有机碳栅格数据，得到土壤有机碳侵蚀空间分布特征图，并利用时间序列回归模型，分析土壤有机碳侵蚀动态变化速率；

S5，根据计算的土壤有机碳侵蚀动态变化速率结果，结合因子贡献分析，得到环境因子对土壤有机碳侵蚀动态变化速率的贡献程度；

S6，根据土壤有机碳侵蚀动态变化速率和环境因子之间的相互联系，得到环境因子对土壤有机碳侵蚀动态变化速率影响过程中的内部互作关系。

2.如权利要求1所述的一种土壤有机碳侵蚀动态变化诱因分析方法，其特征在于：所述地面降雨数据来自于地面气象站，土壤剖面数据为人工采样所得，具体包括土壤类型数据、土壤粒径数据以及土壤有机碳数据；所述降雨栅格数据为TRMM提供的空间分辨率为0.25°，时间分辨率为1天的降雨数据；所述植被栅格数据包括MODIS卫星提供的空间分辨率为250米，时间分辨率为年的MOD44B产品，以及MODIS卫星提供的空间分辨率为500米，时间分辨率为年的MCD12Q1产品；所述土地利用栅格数据来自空间分辨率为1公里的中国多时期土地利用土地覆被遥感监测数据集；所述地形栅格数据来自SRTM提供的DEM数据，空间分辨率为90米。

3.如权利要求2所述的一种土壤有机碳侵蚀动态变化诱因分析方法，其特征在于，在所述S2步骤中，多源时空数据挖掘及时空建模分析的具体方法为：先通过星地多源数据融合算法，结合非参数的二维傅里叶变换和同位协同克里格算法，获取高精度的降雨融合数据，并基于降雨侵蚀力模型计算降雨侵蚀力因子空间分布图；利用Cubist结合随机森林模型的数字挖掘技术，结合土壤可蚀性因子计算模型，模拟土壤可蚀性因子空间分布图；利用土地利用类型栅格数据和植被栅格数据，结合植被盖度，估算植被覆盖因子空间分布图；利用地形栅格数据，结合坡长坡度因子计算模型，得到坡长坡度因子空间分布图；利用土地利用类型栅格数据，结合经验值，得到水土保持措施因子空间分布图。

4.如权利要求1所述的一种土壤有机碳侵蚀动态变化诱因分析方法，其特征在于：在所述S3步骤中，利用所述S2步骤中获取的降雨侵蚀力因子空间分布图、土壤可蚀性因子空间分布图、植被覆盖因子空间分布图、坡长坡度因子空间分布图、水土保持措施因子空间分布图，结合修正土壤流失方程，得到目标区域的土壤侵蚀空间分布图。

5.如权利要求1所述的一种土壤有机碳侵蚀动态变化诱因分析方法，其特征在于：在所述S4步骤中，土壤有机碳侵蚀空间分布特征图由以下公式得到：

土壤有机碳侵蚀＝土壤侵蚀×土壤有机碳

其中，土壤侵蚀由S3步骤得到，土壤有机碳是由SoilGrids250m提供的空间分辨率为250米的土壤有机碳数据。

6.如权利要求1所述的一种土壤有机碳侵蚀动态变化诱因分析方法，其特征在于：在所述S5步骤中，利用增强回归树分析法，分析气候因素、地形因素、土壤因素、植被因素对土壤有机碳侵蚀动态变化速率的贡献程度。