Using topographical attributes to evaluate gully erosion proneness (susceptibility) in two mediterranean basins: advantages and limitations

Empirical multivariate predictive models represent an important tool to estimate gully erosion susceptibility. Topography, lithology, climate, land use and vegetation cover are commonly used as input for these approaches. In this paper, two multivariate predictive models were generated for two gully erosion processes in San Giorgio basin (Italy) and Mula River basin (Spain) using only topographical attributes as independent variables. Initially, nine models (five for San Giorgio and four for Mula) with pixel sizes ranging from 2 to 50 m were generated, and validation statistics were calculated to estimate the optimal pixel size. The best models were selected based on model performance using the area under the receiver operating characteristic (AUC) curve and the generalized cross-validation. The best pixel size was 4 m in the San Giorgio basin and 20 m in the Mula basin. The finest resolution was not necessarily the best; rather, the relationship between digital elevation model resolution and size of the landform was important. The two selected models showed an excellent performance with AUC values of 0.859 and 0.826 for San Giorgio and Mula, respectively. The Topographic Wetness Index and the general curvature were identified as key topographical attributes in San Giorgio and Mula basins, respectively. Both attributes were related to the processes observed in the field and described in the literature. Finally, maps of gully erosion susceptibility were produced for each basin. These maps showed that 22 and 20 % of San Giorgio and Mula basins, respectively, present favourable conditions for the development of gullies.     

A GIS-based approach for gully erosion susceptibility modelling: a test in Sicily, Italy

The aim of this study is to analyze the susceptibility conditions to gully erosion phenomena in the Magazzolo River basin and to test a method that allows for driving the factors selection. The study area is one of the largest (225 km²) watershed of southern Sicily and it is mostly characterized by gentle slopes carved into clayey and evaporitic sediments, except for the northern sector where carbonatic rocks give rise to steep slopes. In order to obtain a quantitative evaluation of gully erosion susceptibility, statistical relationships between the spatial distributions of gullies affecting the area and a set of twelve environmental variables were analyzed. Stereoscopic analysis of aerial photographs dated 2000, and field surveys carried out in 2006, allowed us to map about a thousand landforms produced by linear water erosion processes, classifiable as ephemeral and permanent gullies. The linear density of the gullies, computed on each of the factors classes, was assumed as the function expressing the susceptibility level of the latter. A 40-m digital elevation model (DEM) prepared from 1:10,000-scale topographic maps was used to compute the values of nine topographic attributes (primary: slope, aspect, plan curvature, profile curvature, general curvature, tangential curvature; secondary: stream power index; topographic wetness index; LS-USLE factor); from available thematic maps and field checks three other physical attributes (lithology, soil texture, land use) were derived. For each of these variables, a 40-m grid layer was generated, reclassifying the topographic variables according to their standard deviation values. In order to evaluate the controlling role of the selected predictive variables, one-variable susceptibility models, based on the spatial relationships between each single factor and gullies, were produced and submitted to a validation procedure. The latter was carried out by evaluating the predictive performance of models trained on one half of the landform archive and tested on the other. Large differences of accuracy were verified by computing geometric indexes of the validation curves (prediction and success rate curves; ROC curves) drawn for each one-variable model; in particular, soil texture, general curvature and aspect demonstrated a weak or a null influence on the spatial distribution of gullies within the studied area, while, on the contrary, tangential curvature, stream power index and plan curvature showed high predictive skills. Hence, predictive models were produced on a multi-variable basis, by variously combining the one-variable models. The validation of the multi-variables models, which generally indicated quite satisfactory results, were used as a sensitivity analysis tool to evaluate differences in the prediction results produced by changing the set of combined physical attributes. The sensitivity analysis pointed out that by increasing the number of combined environmental variables, an improvement of the susceptibility assessment is produced; this is true with the exception of adding to the multi-variables models a variable, as slope aspect, not correlated to the target variable. The addition of this attribute produces effects on the validation curves that are not distinguishable from noise and, as a consequence, the slope aspect was excluded from the final multi-variables model used to draw the gully erosion susceptibility map of the Magazzolo River basin. In conclusion, the research showed that the validation of one-variable models can be used as a tool for selecting factors to be combined to prepare the best performing multi-variables gully erosion susceptibility model.