Determination of high-risk factors and related spatially influencing variables of heavy metals in groundwater

Identifying high-risk factors (heavy metals (HMs) and pollution sources) by coupling receptor models and health risk assessment model (HRA) is a novel approach within the field of risk assessment. However, this coupled model ignores the contribution of spatial differentiation to high-risk factors, resulting in the assessment being subjective. Taking Dongting Plain (DTP) as an example, a coupling framework by jointly using the positive matrix factorization model (PMF), HRA, Monte Carlo simulation, and geo-detector was developed, aiming to identify high-risk factors in groundwater, and further explore key environmental variables influencing the spatial heterogeneity of high-risk factors. The results showed that at least 82.86 % of non-carcinogenic risks and 97.41 % of carcinogenic risks were unacceptable for people of all ages, especially infants and children. According to the relationships among HMs, pollution sources, and health risks, As and natural sources were defined as high-risk HMs and sources, respectively. The interactions among Holocene thickness, oxidation-reduction potential, and dissolved organic carbon emerged as the primary drivers of spatial variability in high-risk factors, with their combined explanatory power reaching up to 74%. This proposed framework provides a scientific reference for future studies and a practical reference for environmental authorities in developing effective pollution management measures.