Soil erosion represents a complex ecological issue that is present on a global level, with negative consequences for environmental quality, the conservation and availability of natural resources, population safety, and material security, both in rural and urban areas. To mitigate the harmful effects of soil erosion, a soil erosion map can be created. Broadly applied in the Balkan Peninsula region (Serbia, Bosnia and Herzegovina, Croatia, Slovenia, Montenegro, North Macedonia, Romania, Bulgaria, and Greece), the Erosion Potential Method (EPM) is an empirical erosion model that is widely applied in the process of creating soil erosion maps. In this study, an innovation in the process of the identification and mapping of erosion processes was made, creating a coefficient of the types and extent of erosion and slumps (φ), representing one of the most sensitive parameters in the EPM. The process of creating the coefficient (φ) consisted of applying remote sensing methods and satellite images from a Landsat mission. The research area for which the satellite images were obtained and thematic maps of erosion processes (coefficient φ) were created is the area of the Federation of Bosnia and Herzegovina and the Brčko District (situated in Bosnia and Herzegovina). The Google Earth Engine (GEE) platform was employed to process and retrieve Landsat 7 Enhanced Thematic Mapper plus (ETM+) and Landsat 8 Operational Land Imager and Thermal Infrared Sensor (OLI/TIRS) satellite imagery over a period of ten years (from 1 January 2010 to 31 December 2020). The mapping and identification of erosion processes were performed based on the Bare Soil Index (BSI) and by applying the equation for fractional bare soil cover. The spatial–temporal distribution of fractional bare soil cover enabled the definition of coefficient (φ) values in the field. An accuracy assessment was conducted based on 190 reference samples from the field using a confusion matrix, overall accuracy (OA), user accuracy (UA), producer accuracy (PA), and the Kappa statistic. Using the confusion matrix, an OA of 85.79% was obtained, while UA ranged from 33% to 100%, and PA ranged from 50% to 100%. Applying the Kappa statistic, an accuracy of 0.82 was obtained, indicating a high level of accuracy. The availability of a time series of multispectral satellite images for each month is a crucial element in monitoring the occurrence of erosion processes of various types (surface, mixed, and deep) in the field. Additionally, it contributes significantly to decision-making, strategies, and plans in the domain of erosion control work, the development of plans for identifying erosion-prone areas, plans for defense against torrential floods, and the creation of soil erosion maps at local, regional, and national levels.

Soil erosion is a global problem that negatively affects the quality of the environment, the availability of natural resources, as well as the safety of inhabitants. Soil erosion threatens the functioning of urban areas, which was the reason for choosing the territory of the Master Plan of Belgrade (Serbia) as the research area. The calculation of soil erosion loss was analyzed using the G2 erosion model. The model belongs to a group of empirical models and is based on the synthesis of the equation from the Revised Universal Soil Loss Equation (RUSLE) and the Erosion Potential Method (EPM). The estimation of soil degradation was analyzed in two time periods (2001 and 2019), which represent the time boundaries of the management of the Master Plan of Belgrade. The novel approach used in this research is based on using the land cover inventory as a dynamic indicator of the urbanization process. Land cover was identified using remote sensing, machine learning techniques, and the random forest algorithm applied to multispectral satellite images of the Landsat mission in combination with spectral indices. Climatic parameters were analyzed on the basis of data from meteorological stations (first scenario, i.e., 2001), as well as on simulations of changes based on climate scenario RCP8.5 (representative concentration pathways) concerning the current condition of the land cover (second scenario). A comparative analysis of the two time periods identified a slight reduction in total soil loss. For the first period, the average soil loss value is 4.11 t·ha⁻¹·y⁻¹. The analysis of the second period revealed an average value of 3.63 t·ha⁻¹·y⁻¹. However, the increase in non-porous surfaces has led to a change in the focus of soil degradation. Increased average soil loss as one of the catalysts of torrential flood frequencies registered on natural and semi-natural areas were 43.29% and 16.14%, respectively. These results are a significant contribution to the study of soil erosion in urban conditions under the impact of climate change.

Erosion is a significant environmental challenge in Serbia, shaped by natural and human factors. Pronounced relief, fragile geological substrate, a developed hydrographic network, and a climate characterized by an uneven distribution of precipitation throughout the year make this area prone to activating erosion processes and flash floods whenever there is a significant disruption in ecological balance, whether due to the removal of vegetation cover or inadequate land use. Researchers have recorded approximately 11,500 torrents in Serbia, most of which were activated during the 19th century, a period of significant social and political change, as well as intensive deforestation and the irrational exploitation of natural resources. By the mid-19th century, the effects of land degradation were impossible to ignore. As the adequate assessment of soil erosion intensity is the initial step in developing a prevention and protection strategy and the type and scope of anti-erosion works and measures, this article presents the path that the anti-erosion field in Serbia has taken from the initial observations of erosion processes through the first attempts to create the Barren Land Cadastre and Torrent Cadastre to the creation of the Erosion Potential Method (EPM) and its modification by Dr. Lazarević that resulted in the creation of the first Erosion Map of SR Serbia in 1971 (published in 1983). In 2020, a new Erosion Map of Serbia was created with the application of Geographic Information System (GIS) technologies and based on the original method by Professor Slobodan Gavrilović—the EPM—without the modifications introduced by Lazarević. We compared the 1983 and 2020 erosion maps in a GIS environment, where the change in soil erosion categories was analyzed using a confusion matrix. The updated erosion maps mirror the shift in methodology from a traditional approach (Lazarević’s modification) to the modern GIS-based method (Gavrilović’s original EPM) and reflect technological improvements and changes in land use, conservation practices, and environmental awareness.

Soil erosion assessment in Southeast Europe’s mountainous regions often lacks systematic field validation, limiting confidence in model-based predictions. This study integrates the Erosion Potential Model (EPM) with remote sensing and field verification across 26,570 km² in the Federation of Bosnia and Herzegovina (FBiH) and Brčko District (BD). We developed a two-stage framework: initial GIS-based assessment using digital elevation models, soil maps, climate data, CORINE Land Cover, and Landsat imagery, followed by field calibration at 190 representative sites. Spectral indices (NDVI, BSI) provided dynamic corrections for vegetation cover and visible erosion features. Field validation significantly improved model performance; the erosion coefficient increased from Z = 0.21 to Z = 0.24, while discriminatory power improved AUC from 0.82 to 0.85, with corresponding gains in overall accuracy from 0.78 to 0.84 and F1-score from 0.78 to 0.85. The field-validated model estimated mean annual sediment production of 546.60 m³·km⁻²·year⁻¹, with total erosion material production of 14,074,940.2 m³·year⁻¹. Field calibration revealed substantial spatial redistribution, with medium-to-excessive erosion categories expanding by 30.37%, affecting 1319.12 km² requiring priority intervention. The Kappa coefficient (0.81) confirms high classification reliability. This field-validated framework enables evidence-based identification of degradation hotspots and provides actionable guidance for soil conservation planning in geomorphologically heterogeneous, data-limited regions.

Soil erosion represents one of the key forms of environmental degradation, with pronounced impacts on agriculture, water resources, and ecosystem stability. This study analyzes the spatial distribution of erosion intensity and the production of eroded material in the Sarajevo Canton area using the Erosion Potential Method. The method was enhanced through the application of Geographic Information Systems (GIS), remote sensing, and field surveys, enabling the quantification of erosion processes and the creation of a highly reliable Soil Erosion Map. The analysis integrated a Digital Elevation Model (DEM), pedological characteristics, climatic parameters, land use data (CORINE Land Cover 2018), and vegetation index (NDVI) and bare soil index (BSI). Spatial calculations were performed in raster format with a spatial resolution of 100 m. Model validation was carried out by comparing preliminary and final erosion maps after field verification, with accuracy metrics calculated based on the confusion matrix. The results indicate a high overall classification accuracy (Overall Accuracy = 0.943; Kappa = 0.88), confirming the model's reliability and the stability of erosion process categorization. The study area is predominantly characterized by very weak erosion (about 70% of the total area), while strong and excessive erosion forms are spatially limited but clearly expressed. The average specific production of eroded material is 339.31 m³·km⁻²·yr⁻¹, indicating moderate land degradation across most of the territory. The observed increase in areas affected by excessive erosion in the final mapping phase is a result of methodological improvements, higher spatial resolution, and more precise field verification rather than a real expansion of the processes. The results confirm that the integration of remote sensing, GIS analysis, and field research provides a reliable approach for assessing erosion processes and forms a sound basis for planning soil and water protection measures. This approach contributes significantly to improving erosion monitoring systems and supports sustainable management of natural resources.

Acid deposition cause changes of soil properties, primarily decrease of soil pH. The influence of S i N deposition on terrestrial ecosystems was studied in the catchment of river Rasina, where forest ecosystems occupy 51% of the area. Lower soil pH initiates an increase of trace elements mobility in soil solution. According to this, monitoring of soil acidification status is very important in terms of conservation of surface and ground water quality. Acid deposition in Rasina catcment decrease since 1980. and since 1990. does not exceed a critical value. VSD model simulation indicates that forest ecosystems are still under their negative impact, and that forest ecosystems are in so called recovery delay time (RDT). According to the adopted scenario until 2100 the VSD model predicts a further reduction of soil pH by 0.6 units and in smaller areas soil pH increase by 0.2 units. It can be concluded that the current reduction of acid deposition, for the period 1980-2009, are not sufficient for reaching a steady state of forest ecosystems in Rasina catchment.