This study assessed soil erosion intensity and soil properties across the Crveni Potok catchment in Serbia, a region of diverse morphology, geology, pedology, and vegetation. Soil samples were collected using a regular grid approach to identify the underlying factors contributing to erosion and the most vulnerable areas. Based on ¹³⁷Cs activities and the profile distribution (PD) model, severe erosion (>10 t ha⁻¹ y⁻¹) was predicted at nearly 60% of the studied locations. The highest mean erosion rates were detected for the lowest altitude range (300–450 m), Rendzic Leptosol soil, and grass-covered areas. A significant negative correlation was found between the erosion rates, soil organic matter, and indicators of soil structural stability (OC/clay ratio and S_t), indicating that the PD model successfully identifies vulnerable sites. The PD and RUSLE (revised universal soil loss equation) models provide relatively similar mean erosion rates (14.7 t ha⁻¹ y⁻¹ vs. 12.7 t ha⁻¹ y⁻¹) but significantly different median values (13.1 t ha⁻¹ y⁻¹ vs. 5.5 t ha⁻¹ y⁻¹). The model comparison revealed a positive trend. The observed inconsistencies were interpreted by the models’ spatiotemporal frameworks and RUSLE’s sensitivity to input data quality. Land use stands out as a significant factor modifying the variance of erosion rate, highlighting the importance of land management practices in mitigating erosion.

Despite presenting a practical approach for the characterization of the environmental risk of potentially toxic elements (PTEs) derived from steel production facilities, the analysis of the spatial distribution of bioavailable PTEs concentrations in the soil is frequently overlooked in the management of polluted sites. In this study, the diethylenetriaminepentaacetic acid (DTPA)-extractable forms of PTEs were investigated in soils surrounding the largest Serbian steel production plant. The correlation and geostatistical analysis indicated their pronounced variability suggesting the anthropogenic origin of most investigated elements, apparently from the steel production facility. The detailed visualization of variables and observations derived by self-organizing maps (SOMs) revealed the homologies in PTEs’ distribution patterns, implying the common origin of some elements. These observations were confirmed by principal component analysis (PCA) and positive matrix factorization (PMF). The аpplied approach supports a comprehensive assessment of contaminated sites' ecological and health risks and provides a basis for soil remediation.

Contemporary research has given much attention to explaining and managing gully erosion processes, which significantly compromise environmental quality. The need to reliably quantify soil redistribution has led to the development of the (nuclear) 137Cs technique, herein applied to obtain such estimates within the gully in the Crveni potok catchment, an erosion-prone area of Eastern Serbia. The study also aimed to find relations among the erosion rates and soil physicochemical properties. The erosion rate was significantly correlated (p<0.05) only with the particle size distribution, specifically being negatively correlated with silt and positively correlated with the fine sand fraction. The loss of silt fraction could lead to erosion expansion and sedimentation outside the gully, implying that remediation measures (e.g., reestablishing a vegetation cover) must be taken to prevent the advancement of soil redistribution issues.