Spatial soil information in forests is crucial to assess ecosystem services such as carbon storage, water purification or biodiversity. However, spatially continuous information on soil properties at adequate resolution is rare in forested areas, especially in mountain regions. Therefore, we aimed to build high-resolution soil property maps for pH, soil organic carbon, clay, sand, gravel and soil density for six depth intervals as well as for soil thickness for the entire forested area of Switzerland. We used legacy data from 2071 soil profiles and evaluated six different modelling approaches of digital soil mapping, namely lasso, robust external-drift kriging, geoadditive modelling, quantile regression forest (QRF), cubist, and support vector machines. Moreover, we combined the predictions of the individual models by applying a weighted model averaging approach. All models were built from a large set of potential covariates which included e.g. multi-scale terrain attributes and remote sensing data characterizing vegetation cover.
Model performances, evaluated against an independent dataset were similar for all methods. However, QRF achieved the best prediction performance in most cases (18 out of 37 models), while model averaging outperformed the individual models in five cases. For the final soil property maps we therefore used the QRF predictions. Prediction performance showed large differences for the individual soil properties. While for fine earth density the R² of QRF varied between 0.51 and 0.64 across all depth intervals, soil organic carbon content was more difficult to predict (R² = 0.19-0.32). Since QRF was used for map prediction, we assessed the 90% prediction intervals from which we derived uncertainty maps. The latter are valuable to better interpret the predictions and provide guidance for future mapping campaigns to improve the soil maps.