Assessing the 3D distribution of soil organic carbon by integrating predictions of water and tillage erosion into a digital soil mapping-approach: a case study for silt loam cropland (Belgium)

Although agricultural intensification has generally increased crop yields, it also resulted in a range of environmental issues. These include increased erosion rates and declined soil organic carbon (SOC) stocks. In order to improve our understanding on how erosion impacts the overall SOC storage capacity of croplands, this study analyses the 3D distribution of SOC as a function of water and tillage erosion in a conventionally ploughed field in the Belgian silt loam region. We present a novel methodological framework to integrate the output of an advanced erosion model as a co-variate within a Digital Soil Mapping (DSM) approach with the objective to create detailed SOC maps. More precisely, we combined (i) the Water and Tillage Erosion Model and Sediment Delivery Model (WaTEM/SEDEM), simulating spatial patterns of soil erosion and sediment deposition due to water and tillage erosions, with (ii) a SOC sampling campaign, resulting in a SOC spatial distribution model that considers both types of erosion. The results show that, as compared to plateaus, SOC stocks are nearly half as large along eroding slopes (i.e. convex slope positions affected by tillage erosion and steep slopes affected by water erosion). Yet, they are up to twice as large in areas characterized by sediment deposition (i.e. concave positions due to tillage erosion and foot slope positions due to water erosion). Our results further show that tillage erosion has a significant influence on the SOC stocks in the top 0.7 m and in particular on the top 0.4 m. The influence of water erosion is less strong but mostly significant along the entire depth profile. Overall, this work demonstrates the relevance of considering different erosion processes when aiming to predict spatial patterns of SOC. Considering the top 1 m and a WaTEM/SEDEM application at a resolution of 10 m our 3D SOC modelling approach obtained a coefficient of determination (R2) of 0.62, a relative root mean square error (RRMSE) of 30.5 % and a relative mean absolute error (RMAE) of 26.0 %. While future work may likely lead to further improvements e.g. a more detailed SOC sampling network along the foot slopes and thalwegs in order to obtain a more spatially detailed prediction of the area characterized by depositions due to water erosion, we demonstrate the great potential of existing erosion and deposition models in doing so.