Filling the European blank spot—Swiss soil erodibility assessment with topsoil samples

Soil erodibility, commonly expressed as the K‐factor in USLE‐type erosion models, is a crucial parameter for determining soil loss rates. However, a national soil erodibility map based on measured soil properties did so far not exist for Switzerland. As an EU non‐member state, Switzerland was not included in previous soil mapping programs such as the Land Use/Cover Area frame Survey (LUCAS). However, in 2015 Switzerland joined the LUCAS soil sampling program and extended the topsoil sampling to mountainous regions higher 1500 m asl for the first time in Europe. Based on this soil property dataset we developed a K‐factor map for Switzerland to close the gap in soil erodibility mapping in Central Europe. The K‐factor calculation is based on a nomograph that relates soil erodibility to data of soil texture, organic matter content, soil structure, and permeability. We used 160 Swiss LUCAS topsoil samples below 1500 m asl and added in an additional campaign 39 samples above 1500 m asl. In order to allow for a smooth interpolation in context of the neighboring regions, additional 1638 LUCAS samples of adjacent countries were considered. Point calculations of K‐factors were spatially interpolated by Cubist Regression and Multilevel B‐Splines. Environmental features (vegetation index, reflectance data, terrain, and location features) that explain the spatial distribution of soil erodibility were included as covariates. The Cubist Regression approach performed well with an RMSE of 0.0048 t ha h ha⁻¹ MJ⁻¹ mm⁻¹. Mean soil erodibility for Switzerland was calculated as 0.0327 t ha h ha⁻¹ MJ⁻¹ mm⁻¹ with a standard deviation of 0.0044 t ha h ha⁻¹ MJ⁻¹ mm⁻¹. The incorporation of stone cover reduces soil erodibility by 8.2%. The proposed Swiss erodibility map based on measured soil data including mountain soils was compared to an extrapolated map without measured soil data, the latter overestimating erodibility in mountain regions (by 6.3%) and underestimating in valleys (by 2.5%). The K‐factor map is of high relevance not only for the soil erosion risk of Switzerland with a particular emphasis on the mountainous regions but also has an intrinsic value of its own for specific land use decisions, soil and land suitability and soil protection.