Model averaging for mapping topsoil organic carbon in France

Abstract The soil organic carbon (SOC) pool is the largest terrestrial carbon (C) pool and is two to three times larger than the C stored in vegetation and the atmosphere. SOC is a crucial component within the C cycle, and an accurate baseline of SOC is required, especially for biogeochemical and earth system modelling. This baseline will allow better monitoring of SOC dynamics due to land use change and climate change. However, current estimates of SOC stock and its spatial distribution have large uncertainties. In this study, we test whether we can improve the accuracy of the three existing SOC maps of France obtained at national (IGCS), continental (LUCAS), and global (SoilGrids) scales using statistical model averaging approaches. Soil data from the French Soil Monitoring Network (RMQS) were used to calibrate and evaluate five model averaging approaches, i.e., Granger-Ramanathan, Bias-corrected Variance Weighted (BC-VW), Bayesian Modelling Averaging, Cubist and Residual-based Cubist. Cross-validation showed that with a calibration size larger than 100 observations, the five model averaging approaches performed better than individual SOC maps. The BC-VW approach performed best and is recommended for model averaging. Our results show that 200 calibration observations were an acceptable calibration strategy for model averaging in France, showing that a fairly small number of spatially stratified observations (sampling density of 1 sample per 2500 km2) provides sufficient calibration data. We also tested the use of model averaging in data-poor situations by reproducing national SOC maps using various sized subsets of the IGCS dataset for model calibration. The results show that model averaging always performs better than the national SOC map. However, the Modelling Efficiency dropped substantially when the national SOC map was excluded in model averaging. This indicates the necessity of including a national SOC map for model averaging, even if produced with a small dataset (i.e., 200 samples). This study provides a reference for data-poor countries to improve national SOC maps using existing continental and global SOC maps.