Mapping soil erosion risk using multivariate geostatistics

The assessment of erosion risk at regional scale requires modelling spatial variability of soil erodibility. The traditional soil cartographic approach, in many cases has been considered unsatisfactory for this purpose. This paper presents an approach in which two global risk indices are estimated and mapped using multivariate geostatistics and GIS. The proposed method takes into account some indicators, for which critical threshold values must be defined. This paper presents an application to a study area in Sardinia (Italy), representative of the Mediterranean zones at high erosion risk.