This paper studies the assessment of sensitivity to land degradation of Deliblato sands (the northern part of Serbia), as a special nature reserve. Sandy soils of Deliblato sands are highly sensitive to degradation (given their fragility), while the system of land use is regulated according to the law, consisting of three zones under protection. Based on the MEDALUS approach and the characteristics of the study area, four main factors were considered for evaluation: soil, climate, vegetation and management. Several indicators affecting the quality of each factor were identified. Each indicator was quantified according to its quality and given a weighting of between 1.0 and 2.0. ArcGIS 9 was utilized to analyze and prepare the layers of quality maps, using the geometric mean to integrate the individual indicator map. In turn, the geometric mean of all four quality indices was used to generate sensitivity of land degradation status map. Results showed that 56.26% of the area is classified as critical; 43.18% as fragile; 0.55% as potentially affected and 0.01% as not affected by degradation. The values of vegetation quality index, expressed as coverage, diversity of vegetation functions and management policy during the protection regime are clearly represented through correlation coefficient (0.87 and 0.47).

The erosion potential model was applied to estimate the soil erosion status of rural settlements during the years 1971 and 2011. We used univariate and bivariate local Moran’s I indices to detect and visualize the spatial clustering of settlements with respect to changes in erosion intensity and agricultural land use, as well as their mutual spatial correlation. The study area was differentiated into four statistically significant clusters using the calculated bivariate local Moran’s I indices. The statistical analysis examined the two largest clusters, i.e., the high–high and low–low clusters, and the results of the research indicate that the first four principal components explained 70.50% and 73.47% of the total variance, respectively. In the high–high cluster, the low rates of erosion reduction (average Index Z = 98) in the most significant types of rural settlements were determined according to demographic indicators (i.e., the higher population vitality and population density, the smaller share of the old population and the lower average age of the population) and the large proportion of arable land and Neogene sediments. In the low–low cluster, high erosion reduction rates were detected (average index Z = 64). In this cluster, the more statistically significant influence of natural conditions in combination with demographic–agrarian processes (i.e., the larger share of the old population, the higher average age of the population, the lower vitality index and deagrarization) were decisive factors in changing erosion intensity.

Climate change has a direct impact on flash floods, and indirectly on the environment, society, and economy, due to the rapid development and difficulty of predicting this hydrological phenomenon. The main objective of this study is to assess the potential flash flood hazard areas in the Likodra watershed (218.62 km²), one of the most vulnerable parts to flash floods in Serbia, using the flash flood potential index (FFPI) and analytic hierarchy process (AHP) method. Recurring events from 1995 to this day and the devastating impact on settlements of the analyzed area show that this territory is extremely vulnerable. The data used include hydrological statistics (maximum daily rainfall) and spatial data on watershed geographical characteristics (slope, soils, land use, vegetation, drainage density) obtained or derived from various sources (maps, satellite images, digital databases) which were integrated into the GIS environment. The results indicate a severe flash flood hazard level, with high flash flood susceptibility classes occupying 76.20%, 87.78%, and 91.73% of the area, depending on the considered criteria and weights assigned to them.

Soil erosion is a product of natural and anthropogenic factors and, at the same time, an economic and environmental concern. One of the methods applied to calculate the intensity of erosion is the erosion potential method (EPM), with two possible procedures for determining the average erosion coefficient of an area: analytical and graphical. Using GIS and EPM methods, without field observations of erosion, based on cartographic materials and satellite images, erosion maps were created for 1970 and 2018, for part of the peri-urban area of Belgrade. Based on the created erosion maps, the values of the mean coefficients of erosion, as well as the arithmetic means for the study area, were determined for the settlements. The aim of the study is to assess the reliability of the mean coefficient of soil erosion, obtained from the erosion map created from the cartographic materials and satellite images, without field observations of erosion. Thus, the obtained values of the mean erosion coefficient were compared with the values obtained from the erosion map with field observation and the values obtained by the analytical procedure. Statistical analysis (F test) for 1970 and for 2018 determined a high degree of reliability (p < 0.05) of the mean erosion coefficients of the area obtained from erosion maps that were created from cartographic materials and satellite images without field observation. Regardless of the procedure for determining the mean erosion coefficient, a significant decrease in soil losses was observed, from 10.64 to 5.97 t ha⁻¹year⁻¹ (average annual specific production of sediments, year 1970 and 2018, respectively).

This paper presents the results of an erosion processes research on the territory of the morphological unit of Sutilovacki stream, the prevention measures of same processes, as well as the justification of socio-economic investments and a sensitivity analysis of the economic efficiency parameters. The research includes a period from 2017-2032 year, whereby the evaluation methods of professor Gavrilovic are applied in erosion process intensity assessment, on the whole research area, and universal soil loss equation - USLE, on the agricultural land area. Different rates of soil erosion vulnerability are established on the whole research area by applying these methods. Present erosion processes are prevented by the predicted model of future production from the aspect of soil resource protection in the area of the city of Belgrade, by professor Zlatic (1994). Erosion soil losses are significantly diminished and reduced under a tolerable limit, along with the accomplishment of socio-economic justification of investment, examined by the dynamic methods of internal rate of yield, term for the repayment of invested funds, cost-benefit ratio, and net present value. In addition, an analysis of internal rate of yield sensitivity and term for repayment of invested funds are executed. The obtained data indicate justification of the investment.