In order to assess the rainfall erosivity in the Pannonian basin, several parameters which describe distribution, concentration and variability of precipitation were used, as well as 9 extreme precipitation indices. The precipitation data is obtained from the European Climate Assessment and Dataset project for the period 1961-2014, for 8 meteorological stations in northern Serbia, 5 in Hungary and 1 in eastern Croatia. The extreme values of precipitation were calculated following the indices developed by the ETCCDI. RclimDex software package was used for indices calculation. Based on statistical analysis and the calculated values, the results have been presented with Geographic Information System (GIS) to point out the most vulnerable parts of the Pannonian basin, with regard to pluvial erosion. This study presents the first result of combined rainfall erosivity and extreme precipitation indices for the investigated area. Results of PCI indicate presence of moderate precipitation concentration (mean value 11.6). Trend analysis of FI (mean value 22.7) and MFI (mean value 70.2) implies a shift from being largely in the low erosivity class, to being completely in the moderate erosivity class in the future, thus indicating an increase in rainfall erosivity for most of the investigated area (except in the northwestern parts). Furthermore, the observed precipitation extremes suggest that both the amount and the intensity of precipitation are increasing. The knowledge about the areas affected by strong soil erosion could lead to introducing effective measures in order to reduce it. Long term analysis of rainfall erosivity is a significant step concerning flood prevention, hazard mitigation, ecosystem services, land use change and agricultural production.

The paper aims to provide an overview of the most important parameters (the occurrence, frequency and magnitude) in Vojvodina Region (North Serbia). Monthly and annual mean precipitation values in the period 1946–2014, for the 12 selected meteorological stations were used. Relevant parameters (precipitation amounts, Angot precipitation index) were used as indicators of rainfall erosivity. Rainfall erosivity index was calculated and classified throughout precipitation susceptibility classes liable of triggering soil erosion. Precipitation trends were obtained and analysed by three different statistical approaches. Results indicate that various susceptibility classes are identified within the observed period, with a higher presence of very severe rainfall erosion in June and July. This study could have implications for mitigation strategies oriented towards reduction of soil erosion by water.

Estimation of rainfall erosivity (RE) and erosivity density (ED) is essential for understanding the complex relationships between hydrological and soil erosion processes. The main objective of this study is to assess the spatial–temporal trends and variability of the RE and ED in the central and southern Pannonian Basin by using station observations and gridded datasets. To assess RE and ED, precipitation data for 14 meteorological stations, 225 grid points. and an erosion model consisting of daily, monthly, seasonal, and annual rainfall for the period of 1961–2014 were used. Annual RE and ED based on station data match spatially variable patterns of precipitation, with higher values in the southwest (2100 MJ·mm·ha−1·h−1) and southeast (1650 MJ·mm·ha−1·h−1) of the study area, but minimal values in the northern part (700 MJ·mm·ha−1·h−1). On the other hand, gridded datasets display more detailed RE and ED spatial–temporal variability, with the values ranging from 250 to 2800 MJ·mm·ha−1·h−1. The identified trends are showing increasing values of RE (ranging between 0.20 and 21.17 MJ·mm·ha−1·h−1) and ED (ranging between 0.01 and 0.03 MJ·ha−1·h−1) at the annual level. This tendency is also observed for autumn RE (from 5.55 to 0.37 MJ·mm·ha−1·h−1) and ED (from 0.05 to 0.01 MJ·ha−1·h−1), as for spring RE (from 1.00 to 0.01 MJ·mm·ha−1·h−1) and ED (from 0.04 to 0.01 MJ·ha−1·h−1), due to the influence of the large-scale processes of climate variability, with North Atlantic Oscillation (NAO) being the most prominent. These increases may cause a transition to a higher erosive class in the future, thus raising concerns about this type of hydro-meteorological hazard in this part of the Pannonian Basin. The present analysis identifies seasons and places of greatest erosion risk, which is the starting point for implementing suitable mitigation measures at local to regional scales.

The Fruška Gora Mt., as a dominant orographic complex in the Pannonian plain, was selected for a pioneer geodiversity quantification study area due to its unique geology and soil properties. The methodology is based on the geodiversity quantification assessment approach of Serrano and Ruiz-Flaño (Geogr Helv 62:140-147, 2007). It employed a 500 × 500 m grid approach on several maps (lithological, geomorphological, topographical, and pedological) at scales of 1:50.000 to 1:300.000, together with a 30-m resolution digital elevation model for deriving sub-indices and a topographic roughness. The geodiversity index values (Gd) indicate that the highest geodiversity sites are found on the north, north-east and south-western part of the investigated mountain: in steep-sided valleys, along the horst and loess cliffs facing the Danube River. The obtained results are compared with the already recognized in situ geosite location network. This approach can be applied in the given area for geoheritage protection, conservation, and promotion at different levels (from local to national level). Following the results of this study, the criteria for the definition of conservation areas with abiotic significance should be considered, as there is no legal protection of any kind for the areas with the highest geodiversity index values outside the National Park area. Also, it is a potentially effective tool for supporting decision-making processes regarding the management and conservation of natural areas or regions at different scales with further possible applications in Serbia and elsewhere in Europe.

<p>A detailed analysis of extreme heatwave events in Serbia from the biometeorological point of view is presented in this study.  For this purpose, the newly developed Heat Wave Magnitude Index daily (HWMId), was used on Physiologically equivalent temperature (PET) for Serbia. A series of daily maximum air temperature, relative humidity, the wind was used to calculate PET for the investigated period 1979–2019. HWMId is defined as the maximum magnitude of the heatwaves in a year. Here, the heatwave is characterized as 3 consecutive days with maximum PET above the daily threshold for the reference period 1981–2010. The analysis revealed that during the investigated period the most intensive heat waves occurred in 2007, 2012 and 2015. HWMId values for 2007 were in the range of 8 to 23 indicating extreme heat stress, while for the other two events the values were not as high. Hourly temperatures revealed that the PET values during the day were as high as 55°C. Thus, the mitigation and adaptation to extreme temperature events are of vital importance for humans and their everyday activities. Future investigation should be oriented towards a way to deal with the oppressive heat. Additionally, more research is needed in order to explain and predict these catastrophic events. The main focus of future activities will be on determining the physical causes which lead to the occurrence of extreme heatwaves.</p><p>Keywords: Heat Wave Magnitude Index daily, Physiologically equivalent temperature, Serbia, heat waves</p><p>Acknowledgment: This research is supported by <strong>EXtremeClimTwin</strong> project funded from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 952384</p>

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−1 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−1. 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−1·h−1·yr−1, which is 58% higher than the mean R-factor obtained from RUSLE (330 MJ mm ha−1·h−1·yr−1). 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.

n this study we present an in-depth description of the colorimetric values for the lowest section of the Dukatar Loess Palaeosol Sequence (LPS) pedocomplex S5. Formed during the Marine Isotope Stage (MIS) 13-15, it represents the oldest pedocomplex exposed at the base of the Titel loess plateau (TLP), near the confluence of the Tisa and Danube rivers in Vojvodina (northern Serbia). The results of low-field magnetic susceptibility measurements (χlf) were compared to colour properties (obtained by conventional methods as well as instrumental measuring) and quantified Soil Development Indices (SDI). Of these measurements we found that the Redness Index (RI1) yielded the most useful results, as this index appears most sensitive to lithological changes and soil development intensity. It was also observed that a high level of correlation existed between χlf, and a* chromaticity. The initial results of this study highlight the utility of colorimetric methods as an interdisciplinary tool when evaluating the presence of ferromagnetics, and the application of rock magnetism to the Middle and Upper Pleistocene LPS of the Middle Danube Basin. The presented approach can be used to observe the evolution of climatic and ecological conditions in the given study area, and for establishing correlations between sites extending over the Eurasian LPS provinces.