RESUMO
Revealing the spatial features and source of associated potentially toxic elements (PTEs) is crucial for the safe use of selenium (Se)-rich soils. An integrative risk assessment (GRRRA) approach based on geostatistical analysis (GA), random forest (RF), and receptor models (RMs) was first established to investigate the spatial distribution, sources, and potential ecological risks (PER) of PTEs in 982 soils from Ziyang City, a typical natural Se-rich area in China. RF combined with multiple RMs supported the source apportionment derived from the RMs and provided accurate results for source identification. Then, quantified source contributions were introduced into the risk assessment. Eighty-three percent of the samples contain Cd at a high PER level in local Se-rich soils. GA based on spatial interpolation and spatial autocorrelation showed that soil PTEs have distinct spatial characteristics, and high values are primarily distributed in this research areas. Absolute principal component score/multiple line regression (APCS/MLR) is more suitable than positive matrix factorization (PMF) for source apportionment in this study. RF combined with RMs more accurately and scientifically extracted four sources of soil PTEs: parent material (48.91%), mining (17.93%), agriculture (8.54%), and atmospheric deposition (24.63%). Monte Carlo simulation (MCS) demonstrates a 47.73% probability of a non-negligible risk (RI > 150) caused by parent material and 3.6% from industrial sources, respectively. Parent material (64.20%, RI = 229.56) and mining (16.49%, RI = 58.96) sources contribute to the highest PER of PTEs. In conclusion, the GRRRA method can comprehensively analyze the distribution and sources of soil PTEs and effectively quantify the source contribution to PER, thus providing the theoretical foundation for the secure utilization of Se-rich soils and environmental management and decision making.