Snow avalanches are one of the most devastating natural hazards in the highlands that often cause human casualties and economic losses. The complex process of modeling terrain susceptibility requires the application of modern methods and software. The prediction of avalanches in this study is based on the use of geographic information systems (GIS), remote sensing, and multicriteria analysis—analytic hierarchy process (AHP) on the territory of the Šar Mountains (Serbia). Five indicators (lithological, geomorphological, hydrological, vegetation, and climatic) were processed, where 14 criteria were analyzed. The results showed that approximately 20% of the investigated area is highly susceptible to avalanches and that 24% of the area has a medium susceptibility. Based on the results, settlements where avalanche protection measures should be applied have been singled out. The obtained data can will help local self-governments, emergency management services, and mountaineering services to mitigate human and material losses from the snow avalanches. This is the first research in the Republic of Serbia that deals with GIS-AHP spatial modeling of snow avalanches, and methodology and criteria used in this study can be tested in other high mountainous regions.

The municipality of Štrpce (Southern Serbia) is an area located within Šar Mountain National Park, which is of great ecological importance. Due to the vicinity of settlements, it is necessary to analyze the terrain's susceptibility to natural hazards. The main goal of this research was to determine locations that are highly vulnerable at times of natural hazards (such as earthquakes, erosion, torrential flooding, snow avalanches, and forest fires). The first step in this research was to analyze seismic hazards for a 475 years return period (VII–VIII MCS for the observed area), which was possible by means of Geographic Information Systems. The second step was to determine the intensity of erosion and total sediment production using the Erosion Potential Model. The third step was related to the analysis of the potential of torrential floods using the Flash Flood Potential Index. The Avalanches Potential Index method was used as the fourth step. The fifth step included the analysis of a terrain susceptibility to the occurrence of forest fires. Following the five criteria analysis, weight coefficients were assigned to each of the analyzed parameters by using the Analytical Hierarchy Process (AHP), which provided results of the total susceptibility to natural hazards of the territory of Štrpce. Results indicated that over 45% of the municipality is highly or very highly susceptible to various natural hazards. This article represents a significant step toward a better understanding of natural hazards and it provides a unique knowledge basis for establishing the management and mitigation guidelines and measures, not only within the researched area but at regional and national levels as well.

The use of weather satellite recordings has been growing rapidly over the last three decades. Determining the patterns between meteorological and topographical features is an important scientific job. Cloud cover analysis and properties can be of the utmost significance for potential cloud seeding. Here, the analysis of the cloud properties was conducted by means of Moderate Resolution Imaging Spectroradiometer (MODIS) satellite recordings. The resolution of used data was 1 km² within the period of 30 years (1989–2019). This research showed moderate changing of cloudiness in the territory of Serbia with a high cloudiness in February, followed by cloudiness in January and November. For the past three decades, May has been the month with the highest cloudiness. The regions in the east and south-west, and particularly in the west, have a high absolute cloudiness, which is connected with the high elevation of the country. By means of long term monitoring, the whole territory of Serbia was analyzed for the first time, in terms of cloudiness. Apart from the statistical and numerical results obtained, this research showed a connection between relief and clouds, especially in the winter season. Linear regression MK (Mann-Kendall test) has proven this theory right, connecting high elevation sides with high absolute cloudiness through the year.

The Western Balkans (WB) region is highly prone to water erosion processes, and therefore, the estimation of rainfall erosivity (R-factor) is essential for understanding the complex relationships between hydro-meteorological factors and soil erosion processes. The main objectives of this study are to (1) estimate the spatial-temporal distribution R-factor across the WB region by applying the RUSLE and RUSLE2 methodology with data for the period between 1991 and 2020 and (2) apply cluster analysis to identify places of high erosion risk, and thereby offer a means of targeting suitable mitigation measures. To assess R-factor variability, the ERA5 reanalysis hourly data (0.25° × 0.25° spatial resolution) comprised 390 grid points were used. The calculations were made on a decadal resolution (i.e., for the 1990s, the 2000s, and the 2010s), as well as for the whole study period (1991–2020). In order to reveal spatial patterns of rainfall erosivity, a k-means clustering algorithm was applied. Visualization and mapping were performed in python using the Matplotlib, Seaborn, and Cartopy libraries. Hourly precipitation intensity and monthly precipitation totals exhibited pronounced variability over the study area. High precipitation values were observed in the SW with a >0.3 mm h⁻¹ average, while the least precipitation was seen in the Pannonian Basin and far south (Albanian coast), where the mean intensity was less than an average of 0.1 mm h⁻¹. R-factor variability was very high for both the RUSLE and RUSLE2 methods. The mean R-factor calculated by RUSLE2 was 790 MJ mm ha⁻¹·h⁻¹·yr⁻¹, which is 58% higher than the mean R-factor obtained from RUSLE (330 MJ mm ha⁻1·h⁻¹·yr⁻¹). The analysis of the R-factor at decadal timescales suggested a rise of 14% in the 2010s. The k-means algorithm for both the RUSLE and RUSLE2 methods implies better spatial distribution in the case of five clusters (K = 5) regarding the R-factor values. The rainfall erosivity maps presented in this research can be seen as useful tools for the assessment of soil erosion intensity and erosion control works, especially for agriculture and land use planning. Since the R-factor is an important part of soil erosion models (RUSLE and RUSLE2), the results of this study can be used as a guide for soil control works, landscape modeling, and suitable mitigation measures on a regional scale.

Soil erosion and wildfires are frequent natural disasters that threaten the environment. Identifying and zoning susceptible areas are crucial for the implementation of preventive measures. The Šar Mountains are a national park with rich biodiversity and various climate zones. Therefore, in addition to protecting the local population from natural disasters, special attention must be given to preserving plant and animal species and their habitats. The first step in this study involved collecting and organizing the data. The second step applied geographic information systems (GIS) and remote sensing (RS) to evaluate the intensity of erosion using the erosion potential model (EPM) and the wildfire susceptibility index (WSI). The EPM involved the analysis of four thematic maps, and a new index for wildfires was developed, incorporating nine natural and anthropogenic factors. This study introduces a novel approach by integrating the newly developed WSI with the EPM, offering a comprehensive framework for assessing dual natural hazards in a single region using advanced geospatial tools. The third step involved obtaining synthetic maps and comparing the final results with satellite images and field research. For the Šar Mountains (Serbia), high and very high susceptibility to wildfires was identified in 21.3% of the total area. Regarding soil erosion intensity, about 8.2% of the area is affected by intensive erosion, while excessive erosion is present in 2.2% of the study area. The synthetic hazard maps provide valuable insights into the dynamics of the erosive process and areas susceptible to wildfires. The final results can be useful for decision-makers, spatial planners, and emergency management services in implementing anti-erosion measures and improving forest management in the study area.