- 1School of Public Health and Management, Chongqing Three Gorges Medical College, Chongqing, China
- 2Department of Environment and Food Hygiene, Chongqing Wanzhou District Center for Disease Control and Prevention, Chongqing, China
Heavy metal contamination in agricultural ecosystems poses significant risks to human health. Metals accumulating in soil-crop systems can be absorbed and concentrated by crops. Chronic exposure through the consumption of contaminated crops may lead to severe health consequences. This study investigated eight common heavy metals (Cd, Ni, As, Cu, Hg, Pb, Cr, Zn) in the soil-crop system of Wanzhou District, Chongqing City, and performed a probabilistic risk analysis using Monte Carlo simulation. The study results show that for non-carcinogenic risks, the total hazard quotient for adults (2.365) and children (1.176) both exceed the safety threshold of 1.0, with As contributing significantly to population health risk. For carcinogenic risks, the total carcinogenic risk for adults (2.28 × 10−3) and children (1.11 × 10−3) significantly surpass the unacceptable risk level of 1.0 × 10−4, with As, Cr, Cd, and Ni presenting particularly high carcinogenic risks to the population. Additionally, the sensitivity analysis revealed that the concentration (C) of heavy metals in crops is a key exposure parameter influencing the health risks for both adults and children. Given the significant contributions of As, Cr, Cd, and Ni to health risks, these metals should be prioritized for monitoring and control. Long-term intake of crops with excessive heavy metal content increases health risks, highlighting the urgent need to address heavy metal pollution in urban fragmented farmland soil-crop systems.
1 Introduction
Rapid urbanization and industrialization have made soil heavy metal (HM) pollution an important environmental issue. Soil pollution from HM is characterized by uneven distribution, easy accumulation, non-biodegradability, and low volatility, which lead to the long-term persistence of HM pollutants in the external environment (1, 2). Data from the national communique of soil pollution survey by the Ministry of Environmental Protection of China and the Ministry of Land and Resources in 2014 revealed that the exceedance rate of HMs in arable land soil in China was 19.4%, and HMs were the main pollutants. The exceedance rates for cadmium (Cd), nickel (Ni), arsenic (As), copper (Cu), mercury (Hg), lead (Pb), chromium (Cr), and zinc (Zn) were 7.0, 4.8, 2.7, 2.1, 1.6, 1.5, 1.1, and 0.9%, respectively (3). HM pollution not only deteriorates soil quality and reduces crop yields but also poses a serious threat to human health through bioaccumulation and the food chain (4–7).
HM exposure in the human body mainly occurs through three routes: the respiratory tract, digestive tract, and skin, with dietary intake through the digestive tract accounting for over 90% of the total exposure (8, 9). Long-term consumption of crops contaminated with HMs can adversely affect various human systems, including the nervous, skeletal, circulatory, and immune systems (10, 11). It can also alter enzymatic activity, disrupt hormone secretion, and cause cancer (12–14). Various crops provide the energy, protein, and nutrients required for human physiological activities, and their quality and safety are directly related to human health. In some Asian and African regions, there are similar severe HM pollution problems like those in our study area. For example, in Pakistan, studies on drinking water, soil, and vegetables have found that in some local areas, the Cd, Cr, and Ni content in drinking water is above safe limits. Also, vegetables are contaminated with Cr, Ni, Cd, and Pb. Long-term consumption of these vegetables can increase health risks, especially for children (15–17). A study in East Cameroon shows that the average content of Fe, Al, Cr, Mn, Ba, Zn, Pb, Ni, Cu, Co, As, and Se (mg/kg) is higher than the FAO–set safety limits. The carcinogenic risk index indicates that Cr has the highest cancer risk, with children being more vulnerable than adults (18). Studies in India and Bangladesh also show similar HM pollution problems (19–22). This demonstrates the global impact of HM pollution and its serious burden on public health. In China, relevant studies also indicate the severity of heavy metal pollution. For instance, a study conducted in Wenling City, Zhejiang, showed severe Cd contamination in rice, with Cd, Cr, and As being the chief contributors to its overall carcinogenic risk (23). Similarly, a study in Shanghai found contamination with Pb, Cd, and As in vegetables (24). Furthermore, an analysis of HM pollution in vegetables from urban and suburban areas in China revealed numerous samples with excessive levels of Cd, Hg, and Pb in urban and suburban vegetables, indicating the seriousness of HM pollution (25). Economic development and the improvement of living standards have increased the variety and quantity of various crops consumed and the risk of HM exposure from crops (26).
The assessment of soil environmental quality typically focuses on HMs in the soil, and methods used for evaluation include the single pollution index, Nemerow's comprehensive pollution index, geo-accumulation index, potential ecological risk index, and target hazard quotient index (27–30). However, these methods mainly rely on background values of HMs and screening values, focusing more on the pollution status of individual or multiple HMs in the soil (31, 32), and are unable to conduct a comprehensive pollution assessment within the soil-crop system. In practice, HM pollution often involves both soil and crops, necessitating paired sampling of soil and crops and the use of integrated methods to assess the environmental risks of HMs. This involves considering the mutual influence between HMs in soil and agricultural products on the environmental quality of farmland soil. The Impact Index of Comprehensive Quality (IICQ), comprising the Impact Index of Comprehensive Quality in soil (IICQs) and the Impact Index Comprehensive Quality in agricultural products (IICQap), enables a comprehensive evaluation of soil environmental quality and agricultural product safety related to HMs (33). This approach has been successfully used to assess the HM pollution situation in the soil-crop system (34–36).
Wanzhou District in Chongqing City is an ecologically sensitive region and lies at the heart of the Three Gorges Reservoir (37). The rapid development of urbanization and industrialization around the Three Gorges Reservoir has significantly increased the discharge of industrial and urban wastewater, along with the emission of other pollutants. Elevations in the amount of pollutant emissions facilitate the entry and accumulation of HMs in agricultural soils, posing significant risks to local agriculture and public health (38, 39). Wanzhou District is a typical karst region characterized by a high geological background. Carbonate rocks serve as an important reservoir for HM elements, which can be enriched in the soil through weathering and enter the human body through the food chain, posing a threat to human health (40). Thus, we speculate that this region may experience heavy metal pollution, potentially impacting local population health via the soil–crop system. This study investigates HM pollution in soils and crops in the fragmented farmlands of Wanzhou District, Chongqing, and their potential risks to local food safety. The objectives of this study are: (1) to determine the pollution levels and spatial distribution of HMs in Wanzhou District; (2) to analyze the levels of HM pollution in the soil-crop system; (3) to explore the possible sources of HMs using correlation analysis and Mantel test analysis; (4) to assess the potential health risks of HMs in crops to adults and children. This study enhances our understanding of heavy metal pollution in local soils and crops, as well as the associated health risks. The research results not only help evaluate the heavy metal environmental exposure level, but also provide a scientific basis for subsequent policymaking. Through this study, we aim to reduce the risk of heavy metal environmental exposure and minimize the health hazards caused by the pollution.
2 Materials and methods
2.1 Study area and sample collection
The study area is located in Wanzhou District (30°3′50″-310′18″N, 107°2′22″-108°3′52″E), in the northeastern part of Chongqing, at the heart of the Three Gorges Reservoir region. The study area has an approximate area of 3,456.41 km2, with an administrative division that includes 14 sub-districts, 27 towns, and 11 townships, and a total population of roughly 1.55 million people. The district experiences a subtropical monsoon climate, with distinct seasons, an average annual temperature of 15–24°C, and an average annual precipitation of 1,000–1,200 mm, with the flood season concentrated from May to September. By 2022, the region had a total agricultural acreage of 83,908 hectares, with paddy fields accounting for 30,622 hectares (36.49%) and dry land making up 53,286 hectares (63.51%). The main soil types include purplish soil, paddy soil, and yellow soil. The research focused on fragmented cropland in the district. Based on preliminary soil monitoring data, the sampling sites covered five representative towns: Lihe, Longju, Longsha, Yanshan, and Luotian town. In 2023, 50 crop samples, including four crop types (26 grains, nine vegetables, three leaf vegetables, 12 tubers) and their rhizosphere soil samples, were collected from the study area. Figure 1 shows the distribution of sampling spots. Soil sample collection was conducted in accordance with the “Soil Environmental Monitoring Technical Specifications” (HJ/T166-2004) (41). Soil samples were collected from a depth range of 5–20 cm using the plum blossom grid method to ensure uniform sampling over an area of 1 m2. Subsampling was performed using the quartering method, retaining 1 kg of soil per sample. Additionally, we also collected samples of the corresponding crops grown on the soil. After mixing, we used the quartering method for subsampling, retaining 500 g as one sample, and properly recorded the sampling process.
2.2 Sample testing analysis, methods, and quality control
In the laboratory, soil samples were dried, ground, and sieved through a 100-mesh sieve. Subsequently, an HCl-HNO3-HF-HClO4 solution was added, and the samples were heated at 150°C for 48 h to perform digestion (42–45). Grain samples were rinsed with purified water, dried, threshed, dehulled, and ground to a 40-mesh size and then microwave-digested to determine the concentration of HMs. Vegetable, leaf vegetable, and tuber samples were cleaned with purified water, homogenized, and microwave-digested to determine the concentration of HMs. A programmed temperature ramp was applied to gradually increase the temperature to 185°C, followed by 40-min microwave digestion. The soil pH value was determined using the potentiometric method as specified in the NY/T 1377-2007 (46). The Hg content in soil was measured using the catalytic thermal decomposition-cold atomic absorption spectrophotometry method (HJ 923-2017) (47). The concentration of HMs such as Pb, Cd, Cr, As, Cu, Zn, and Ni in soil were measured using the aqua regia extraction-inductively coupled plasma mass spectrometry (ICP-MS) method (HJ 803-2016) (48). The presence of Hg in crops was assessed using direct mercury measurement (Method II of GB 5009.17-2021) (49), while Pb, Cd, Cr, As, Cu, Zn, and Ni were measured using the inductively coupled plasma mass spectrometry (ICP-MS) method (Method I of GB 5009.268-2016) (50).
Soil and crop detection indicators both adopt standard methods. Instruments used are the Inductively Coupled Plasma Mass Spectrometer NexION 2000 (PerkinElmer, USA), the Fully Automatic Direct Mercury Analyzer MAX-S (Leiberte), and the Multi–parameter Tester pH Meter (Mettler–Toledo). The main reagents used are nitric acid [Thermo Fisher Scientific (China) Co., Ltd.] and hydrochloric acid (Chengdu Kelong Chemicals Co., Ltd., chemically pure). All instruments are verified/calibrated, and all standards are within the validity period.
For quality assurance (QA) and quality control (QC), all sampling and testing personnel were uniformly trained, all equipment was calibrated and verified, all standard materials were within their valid period, and quality control measures such as laboratory blanks, parallel samples, and recovery tests were employed during testing. The results of all quality control measures were within the acceptable range.
2.3 Evaluation methods
2.3.1 Index of geoaccunulation
The geoaccumulation index (Igeo) was introduced in the early 1960s to assess the elemental contamination level in sediment. Igeo values were calculated using the following equation (51, 52):
where Csoil is the measured HM concentration in the sample, Cbi is the average HM concentration in the study area, k is the empirical coefficient, which, in general, is 1.5.
The Igeo can be classified into different grades: Class 0, Igeo ≤ 0, uncontaminated; Class 1, 0 < Igeo ≤ 1, mildly contaminated; Class 2, 1 < Igeo ≤ 2, moderately contaminated; Class 3, 2 < Igeo ≤ 3, moderately to heavily contaminated; Class 4, 3 < Igeo ≤ 4, heavily contaminated; Class 5, 4 < Igeo ≤ 5, heavily contaminated to extremely contaminated; Class 6, Igeo ≥ 5, extremely contaminated.
2.3.2 Bioaccumulation factor
The bioaccumulation factor (BCF) reflects the migration characteristics of HMs with soil-crop systems and can be calculated using the following equation (53, 54):
where Capi is the HM contents in the edible parts of crops, and Csoil denotes HM content in the rhizosphere soil samples.
2.3.3 IICQ index
The Impact Index of Comprehensive Quality (IICQ) was used to assess the relative effects of HM contamination on soil quality at particular sites under individual or composite pollution conditions. Compared to previous assessment methods, IICQ is more scientific and rational as it incorporates HM background values, HM screening values in soil, and permissible limits of HM contents in crops. This approach allows for the simultaneous consideration of soil environmental quality and crop quality, facilitating a quantitative and comprehensive estimation of HMs in soil-crop systems. Furthermore, IICQ can help compare regional soil environmental quality as it can be easily visualized on maps at different scales. In this study, the IICQ of soil-crop systems included the Impact Index of Comprehensive Quality in soil (IICQs) and the Impact Index Comprehensive Quality in agricultural products (IICQap). These indices were calculated as follows (33, 53, 55):
(1) Relative impact equivalent (RIE) in soil:
Where Pssi is the single pollution index of element I; N is the number of elements; Csoil and Csi are the content and screening value of element I, respectively; n is the oxidation number of element i (Pb,2; Hg,2; As,5; Cr,3; Cd,2; Cu,2; Ni,2; Zn,2).
(2) The degree of deviation between detected and background value (DDDB) in soil:
where Psbi is the ratio of HM measured value to its background value in soil for element i, and Cbi is the background value of element i.
(3) The degree of deviation of the soil environmental quality standard from the background value (DDSB):
Where the symbols correspond with those in the equations mentioned above.
(4) The quality index of agricultural products (QIAP):
where Papi is the ratio of HM measured value to its threshold limit value of crops for element i and Capi is the measured content of element i in crops, and Cisi is the threshold limit value of element i in crops.
(5) The impact index of comprehensive quality (IICQ):
where IICQ comprises IICQs and IICQap, which are the respective impact indices for soil and agricultural products. X and Y denote the number of soil measurements exceeding the soil standard values and soil background values, respectively. Meanwhile, Z denotes the number of agricultural products exceeding the pollutant limit standards; k is the factor of background correction, commonly set as five.
Based on IICQs, IICQap, and IICQ, the HM pollution risk was classified into five grades as follows: clean (IICQ ≤ 1), slight pollution (1 < IICQ ≤ 2), light pollution (2 < IICQ ≤ 3), moderate pollution (3 < IICQ ≤ 5), and heavy pollution (IICQ > 5).
2.4 Health risk assessment
Health risk assessment method proposed by the United States Environmental Protection Agency can be used to assess the health risks posed by HMs. Since the carcinogenic and non-carcinogenic effects caused by environmental factors are different in risk characterization, we divided them into carcinogenic risk and non-carcinogenic risk when conducting quantitative risk estimation. The hazard quotient (HQ) and cancer-risk index (CR) were adopted to evaluate the potential non-carcinogenic risks and carcinogenic risks from HMs in crops, respectively. Multiple HMs may have a synergistic effect, so the total hazard quotient (THQ) and total carcinogenic risk (TCR) of crops were determined based on the sum of the HQ and CR, respectively. These were calculated through the following equation (56–58):
where Capi was the HM content in crops (mg·kg−1), IR was the ingestion rate of crops (kg·day−1), EF was the exposure frequency (day·year−1), ED was the exposure duration (year), BW was the body weight (kg), AT was the average exposure time for carcinogenic (non-carcinogenic) effects (day), RfD (oral reference dose) was the corresponding reference toxicity threshold dose (mg·kg−1·day−1), and SF was slope factors (kg·day·mg−1).
HQ and THQ values < 1.0 were considered safe, but values exceeding 1.0 suggested potential adverse effects on human health. CR and TCR values above 1.0 × 10−4 indicated unacceptable risks, values below 1.0 × 10−6 implied no obvious hazard, and values between 1.0 × 10−6 and 1.0 × 10−4 were acceptable. In order to avoid the uncertainties from using single-point deterministic parameters, which may result in either overestimation or underestimation of risks, Monte Carlo simulation with 10,000 iterations was employed to evaluate the population health risk. The probability distributions and values for the mentioned parameters are detailed in Supplementary Table S1.
2.5 Statistical analysis
All spatial distribution maps were generated using ArcGIS 10.2. Statistical analyses were performed using Excel 2016 and IBM SPSS 27. RStudio and Origin 2024 were used for data visualization. Additionally, Monte Carlo simulation was performed using Crystal Ball software.
3 Results
3.1 Characteristics of HMs in soil
The physicochemical properties and amounts of potentially toxic elements (PTEs) in soil samples are shown in Table 1. Among the eight HMs detected, only Hg (20%), As (48%), Cd (28%), Cu (8%), Ni (4%), and Zn (70%) exceeded the background values in Chongqing, with Zn showing the most severe exceedance (59). According to the Soil Environmental Quality—Risk Control Standard for Soil Contamination of Agricultural Land in China (GB 15618-2018) (60), only 14% of the soil samples for Cd exceeded the risk screening values in the study area, indicating an enrichment of soil Cd in the region (Supplementary Table S2). This enrichment is likely associated with rapid urbanization and industrialization, which causes pollutants to enter the soil through atmospheric deposition, leading to elevated levels (61). The coefficient of variation (CV) is used to assess whether HMs are evenly distributed and can reflect the influence of anthropogenic activities (62). The CVs for Pb, Hg, As, Cr, Cd, Cu, Ni, and Zn were 16.8, 103, 29.8, 13.9, 43.6, 36.0, 16.9, and 34.3%, respectively. The CV for Hg was >1, indicating significant spatial variability in its distribution, suggesting anthropogenic activities are the primary source of this variation.

Table 1. Descriptive statistics for soil potentially toxic elements and soil physicochemical properties (n = 50; units in mg·kg−1).
3.2 Spatial distribution characteristics of HMs contamination in soil
The spatial distributions of the contents of HMs (Pb, Hg, As, Cr, Cd, Cu, Ni, Zn) in the study area were obtained using the Kriging interpolation method (Figure 2). The spatial distribution of Pb and As showed a clear similarity, mainly concentrated in the northern and southwestern parts of the study area. Sampling points with high Cr and Ni concentrations were mainly distributed in parts of the northwest, southwest, and southeast, with their spatial distribution showing similar patterns. The areas with high Cu and Zn values were both concentrated in the eastern part of the study area, with their spatial distribution showing similar patterns. Areas with high Hg and Cd concentrations were relatively few, occurring only in sporadic spots. The area with a high value for Hg was located in the northwest, while the area with a high value for Cd was in the southwest.
3.3 HMs contamination in soil
Figure 3 and Supplementary Table S3 show the Igeo of eight HMs in 50 soil samples. The mean and standard deviation of the Igeo of Pb, Hg, As, Cr, Cd, Cu, Ni and Zn were −1.08 ± 0.25, −1.15 ± 0.91, −0.61 ± 0.45, −0.91 ± 0.24, −0.80 ± 0.59, −1.22 ± 0.45, −0.91 ± 0.28, and −0.42 ± 0.50, respectively. The max Igeo values of Pb, Cr, and Ni were below 0, suggesting that these metals did not contaminate the soil in the study area. As, Cd, Cu, and Zn showed mild contamination in 12, 6, 2, and 12% of the samples, respectively. Mild and moderate Hg contamination was found in 4 and 6%, respectively. These findings indicate that the soil in the study area was mainly contaminated by Hg, As, Cd, Cu, and Zn, with Hg having a higher contamination risk.
3.4 Crops HM contents and bioaccumulation factor evaluation
The characteristics of HMs contents in crops are shown in Table 2. The content of Pb in crops was within the safety limits specified by the Chinese National Food Safety Standards for Contaminants in Food (63), and the contents of Cu and Zn in crops were below the thresholds set by the Chinese Health Standards for Cu and Zn limits in food (64, 65). In crops, a small number of samples for Hg, As, and Cd showed standard exceedances, with exceedance rates of 6, 6, and 8%, respectively. However, Cr and Ni exhibited significant exceedances, with more than half of the samples exceeding the standards for these elements, which requires attention. Compared with other crops, the concentrations of As, Cu, Zn, and Hg in grains were significantly higher (P < 0.05), which may be related to differences in HM accumulation patterns among crops.
The BCF serves as a valuable indicator for assessing the ability of crops to uptake HMs from the soil, providing a quantitative evaluation of their cumulative characteristics and potential harm to crops (66). The study showed significant differences in the enrichment capacity for eight HMs among different crops (Figure 4 and Supplementary Table S4). These differences were largely attributed to the unique growth habits of each crop and their varying tolerance levels to different HMs. Additionally, the absorption of different HMs by crops may involve synergistic or antagonistic effects, further influencing the enrichment capacity of different crops for HMs. Among all crops, grains demonstrated a higher enrichment capacity for eight HMs, which may be related to their longer growth cycles. The enrichment capacity for As, Cu, Zn, and Hg in grains was significantly higher than in other crops (P < 0.05). Furthermore, this study found that the average enrichment coefficients for Cd across four types of crops were relatively high, likely due to the higher content of Cd in the local soil (61).
3.5 IICQ index evaluation
The average IICQ value for the soil-crop system was 2.97, with pollution levels classified as clean, slight, light, moderate, and heavy in 6, 18, 40, 26, and 10% of the samples, respectively. The mean IICQs value for soil HMs was 1.15, with no heavy pollution in all samples. The soil was clean, slightly polluted, moderately polluted, and heavily polluted in 54, 30, 6, and 10% of the samples, respectively. For crops, the mean IICQap value was 1.82, with no samples showing heavy pollution. The majority of the samples had light pollution (52%), followed by clean (26%), moderate pollution (12%), and slight pollution (10%; Supplementary Table S5). Crops HMs contributed to approximately 61.28% of the overall IICQ for the soil-crop system, explaining crops as the main source of HM exposure. However, soil still contributed 38.72% of HM exposure. The study found differences in the average IICQap values among various types of crops. The average IICQap value for leaf vegetables was 0.83, indicating no pollution. The average IICQap values for vegetables and tubers were 1.38 and 1.20, respectively, indicating slight pollution levels. The average IICQap value for grains was 2.36, indicating light pollution levels (Supplementary Table S6). The spatial distribution of the IICQ values (Figure 5) revealed higher IICQs in the northwest and southwest regions, while higher crop IICQap values were mainly observed in the southeast. The spatial correlation between IICQs and IICQap was low, which may be influenced by the types of crops cultivated and their growth habits (67), causing the spatial distribution of IICQs and IICQap to be inconsistent.
3.6 Correlation analysis of HMs in soil-crop systems
Pearson's correlation analysis results of eight HMs in the soil-crop system are shown in Figure 6 and Supplementary Table S7. HMs in soil and crops were significantly correlated at the 0.01 or 0.05 levels. The correlation analysis results showed that the pH value had a weak positive correlation with Cu in soil (P < 0.01), a moderate positive correlation with Zn in soil (P < 0.01), and a weak negative correlation with Hg in crops (P < 0.05). Pb in the soil had a moderate positive correlation with Cr and Ni in soil (P < 0.01), a strong positive correlation with As in soil (P < 0.01), and a weak positive correlation with Cd (P < 0.01), Cu (P < 0.05), and Hg (P < 0.05) in soil, and a weak negative correlation with Ni in crops (P < 0.05). Cd in the soil had a weak positive correlation with Cr and Ni in soil (P < 0.05) and a weak negative correlation with Ni in crops (P < 0.05). As in soil had a weak positive correlation with Cr (P < 0.05), Cu (P < 0.05), and Ni (P < 0.01) in soil, and a weak negative correlation with Cr (P < 0.01), Ni (P < 0.01), and Zn (P < 0.05) in crops. Cr in soil had a very strong positive correlation with Ni in soil (P < 0.01). Cu in the soil had a strong positive correlation with Zn in soil (P < 0.01). Hg in the soil had a weak negative correlation with Ni in crops (P < 0.05). Pb in crops had a weak positive correlation with Cu in crops (P < 0.05). Cd in crops had a moderate positive correlation with Cr, Cu, Zn, and Ni in crops (P < 0.01). As in crops had a moderate positive correlation with Cr and Ni in crops (P < 0.01), a weak positive correlation with Cu in crops (P < 0.01), and a strong positive correlation with Zn in crops (P < 0.01). Cr in crops had a moderate, strong, and extremely strong positive correlation with Cu (P < 0.01), Zn (P < 0.01), and Ni (P < 0.01) in crops, respectively. Cu in crops had a strong and moderate positive correlation with Zn (P < 0.01) and Ni (P < 0.01) in crops, respectively. Zn in crops had a strong positive correlation with Ni in crops (P < 0.01).

Figure 6. Correlation matrix of eight HMs of soil and crop in soil-crop systems. S represents soil, and C represents crop.
The results of the Mantel test for the eight HMs in the soil-crop system are shown in Figure 7. The right half of the graph shows the results of the Mantel test analysis, displaying the correlation between the concentrations of HMs in soil and crops. The left half of the graph shows the correlation of the concentrations of HMs in soil. The figure displays the correlations among different factors within the same matrix and between different matrices. There was a significant correlation (P < 0.01) between As in soil and eight HMs in crops, but the strength of the correlation was weak (r < 0.2). Other HMs in the soil did not show significant correlations with HMs in crops (P > 0.05, r < 0.2).

Figure 7. Mantel test of eight HMs of soil and crop in soil-crop systems. S represents soil, and C represents crop.
3.7 Health risk assessment of HMs in crops based on Monte Carlo model
3.7.1 Uncertainty analysis of non-carcinogenic health risk
For non-carcinogenic risks, the rank of the average HQ values was As > Zn > Cd > Pb > Hg > Cr > Cu > Ni for both children and adults (Figures 8a–h). Among these HMs, the average HQ values of As were 1.899 and 0.923 for adults and children, with about 59.9 and 29.88% of crops samples exceeding 1.0, respectively (Figure 8c). The average HQ values for other HMs were lower than 1.0, with only a few HQ values for Cd and Cr exceeding 1.0. Besides, the average THQ value for adults was 2.365 [90% confidence interval (CI): 0.656–6.104], and the average THQ value for children was 1.176 (90% CI: 0.324–3.028; Supplementary Table S8), which were far higher than the guideline threshold (1.0). The excess rates were 80.80% for adults and 40.88% for children (Figure 8i).

Figure 8. The probability distribution for health risk index of non-carcinogenic HMs in crops. Hazard quotient (HQ) for (a) Pb, (b) Hg, (c) As, (d) Cr, (e) Cd, (f) Cu, (g) Ni, (h) Zn, and index for (i) total hazard quotient (THQ) [the red and blue dashed lines indicated the mean values, while the black indicated the guideline value for non-carcinogenic risk (1.0)].
3.7.2 Uncertainty analysis of carcinogenic health risk
For carcinogenic risks, the average CR values followed the order of Ni > Cr > Cd > As > Pb for both adults and children (Figures 9a–f). The average CR values of the above HMs (except Pb) all exceeded the unacceptable risk of 1.0 × 10−4. However, only about 14.39% (adults) and 0.24% (children) of CR values for Pb were acceptable (Figure 9a). Additionally, the average TCR value for adults was 2.28 × 10−3 (90% CI: 8.33 × 10−4 to 4.65 × 10−3), and the average TCR value for children was 1.11 × 10−3 (90% CI: 4.06 × 10−4 to 2.27 × 10−3). Nearly 100% of the TCR values surpassed 1.0 × 10−4 (Supplementary Table S8).

Figure 9. The probability distribution for health risk index of carcinogenic HMs in crops. Carcinogenic risk (CR) for (a–e), and index for (f) total carcinogenic risk (TCR) [the red and blue dashed lines indicated the mean values, while the black indicated the guideline value for carcinogenic risk (10−4). The green dashed lines presented the acceptable carcinogenic risks value that of 10−6].
3.7.3 Sensitivity analysis of Monte Carlo model
Sensitivity analysis was conducted to estimate the weight of each parameter in health risk assessment. The larger the sensitivity value, the greater the impact on risk. The results indicated that the concentration (C) of HMs in crops was a crucial exposure parameter influencing the health risks for adults and children. For the HQ value, the exposure frequency also had an impact, contributing 18.0% to the variance of HQ in adults and 17.7% in children. Parameters such as body weight and oral reference dose (RfD) were negatively correlated with health risk. The carcinogenic risk was also influenced by the cancer slope factor (SF) and exposure frequency (EF). The slope factor contributed to 48.8 and 49.6% of the variance in CR for adults and children, respectively, while exposure frequency contributed to 17.2% for adults and 17.1% for children. Body weight contributed < 10% to health risk (Supplementary Tables S9, S10). Among the analyzed HMs, As contributed the most to the total non-carcinogenic risk, with 90.29% for adults and 90.40% for children. In terms of TCR, Ni, Cr, and Cd contributed more than As, with Ni at 44.40% for adults and 44.53% for children, Cr at 44.43% for adults and 43.26% for children, and Cd at 32.79% for adults and 33.88% for children (Figure 10). Overall, As, Ni, Cr, and Cd contributed significantly to the total health risk.

Figure 10. Sensitivity analyses of HMs for (a) total hazard quotient and (b) total carcinogenic risk.
4 Discussion
4.1 Factors influencing the distribution of heavy metal concentrations in crops and soil
In our study, the concentration of Zn in soil samples was higher than the background values in Chongqing in many cases and the average concentration of Zn in a study on HM pollution in Chinese agricultural and urban soils over the past 20 years; however, it was lower than the concentration of Zn in the HM risk assessment of the soil-crop system along the Yangtze River in Nanjing (53, 68). Greater attention is required for the enrichment of Cd in soil. This enrichment is related to the widespread occurrence of carbonate rocks rich in Cd in the southwestern provinces of China and the impact of urbanization, which leads to the irrigation of sewage and excessive use of chemical fertilizers and pesticides. Unlike the concentrations of HM in the soil, the exceedance of Cr and Ni in crops is more severe. In this study, the content of Cr in crops was higher than that in crops in Gansu, China, and lower than that in crops near non-ferrous metal mining areas in Yunnan (69, 70). The content of Ni in crops was higher than that in grains in Anhui Province, China, and lower than that in vegetables planted in roadside soils in Pakistan (71, 72). This finding indicates that the concentration of HMs in crops is related to the type and growth characteristics of the crops and correlates with the pollution levels of HMs in soil (73–76). Further analysis revealed a significant link between soil arsenic (As) and its accumulation in graminous crops, particularly rice. This may relate to the structure of the iron plaque on rice roots. In As-contaminated soils, this plaque might enhance arsenic uptake (77). Under flooded conditions, soil As (V) is easily reduced to As (III), which has weaker adsorption to soil minerals than As (V). This leads to higher soil As concentrations and easier uptake by rice (78).
4.2 Health risks for different populations through crop ingestion
HMs in crops can be transferred to the human body through the food chain, leading to several health risks. The accumulation of various HMs in crops is a relatively complex process influenced by external environmental and anthropogenic factors. HMs enter crops through contaminated agricultural soil, sewage, atmospheric deposition, and fertilizers containing HMs (79, 80). Through the absorption mechanism of crop root systems, HMs are transferred to various parts of the plant, such as tubers, leaves, fruits, and seeds, by specific transport proteins (81). Long-term consumption of crops contaminated with HMs poses health risks to human health, including kidney damage, osteoporosis, neurological damage, cardiovascular diseases, and increased cancer risk (82–86). The unique physiological characteristics of children predispose them to the toxic effects of HMs. Long-term intake of crops contaminated with HMs can have a negative impact on their cognitive development and physical growth (87, 88). In this study, the average THQ for adults and children were 2.365 and 1.176, respectively, which were greater than the standard value of 1.0, suggesting that consuming crops is associated with a higher non-carcinogenic risk. The average THQ in this study was much higher than the values of 0.0441 for adults and 0.324 for children near the coal mines in Jinzhong City, Shanxi Province (89). Compared to the THQ values in this study, a study on HM contamination in crops in Iran reported a lower value for adults (1.67) but a slightly higher value for children (1.28) (90). Only As had an average HQ value exceeding the standard value of 1.0 in 59.90% of samples for adults, while As had an average HQ value close to the standard value of 1.0 in 29.88% of samples for children. This finding indicates that As poses a significant health risk to adults, but the health risk to children should not be ignored. Cd (adults: 0.5%, children: 0.05%), Cr (adults: 0.34%, children: 0.06%), Zn (children: 0.05%) had average HQ values below the standard value of 1.0, suggesting that Cd, Cr, and Zn may not pose a health risk to the population at most exposure levels, but the cumulative effect still warrants attention (91, 92).
The average TCR in this study was 2.28 × 10−3 for adults and 1.11 × 10−3 for children, both significantly higher than the unacceptable level of 1.0 × 10−4, indicating that HM pollution in crops is prone to carcinogenic risk. The TCR values for both adults and children were lower than the values of 2.92 × 10−3 for adults and 4.96 × 10−3 for children near the Leng-shuikeng mine in Jiangxi (93). However, the TCR for children in this study was higher than that reported in agricultural soil in Lanzhou (94). For individual HMs, As (adults: 72.61%, children: 45.47%), Cr (adults: 95.72%, children: 85.22%), Cd (adults: 83.15%, children: 61.83%), and Ni (adults: 96.79%, children: 87.20%) had average CR values greater than the standard value of 1.0 × 10−4, indicating that these metals pose a higher carcinogenic risk to the population and require attention (25, 95–97).
Several studies show that adults are at greater health risk than children due to higher crop consumption and longer exposure time (98–100). However, other studies demonstrate that children are at greater health risk than adults due to their lower body weight and poorer tolerance to toxic HMs (101–103). HMs affect public health through various routes, such as crops, drinking water, soil, and air pollution (104). Therefore, efforts to identify and control the key exposure factors in crops will reduce the population health risks.
4.3 The impact sensitivity analysis using Monte Carlo simulation
Health risk assessments for HM hazards based on Monte Carlo simulation involved many uncertainties. Differences in regions, ethnicities, and ages often lead to variations in model parameters (e.g., exposure frequency and body weight), which can cause deviations in actual health risk levels (105–107). Identifying exposure parameters that significantly impact health risk assessments is essential for reducing uncertainties in the risk assessment process. Reporting the influence and sensitivity of the assessment results can lead to a more accurate evaluation of the impact of heavy metals on population health risks. This study focused on the differences in health risks between adults and children and selected parameters related to human exposure for Monte Carlo simulation and sensitivity analysis. The findings showed no significant differences between the parameters. For both non-carcinogenic and carcinogenic risk assessments, crop intake was found to be the most sensitive parameter. This finding is consistent with previous studies on rice HM pollution in the Minxi region and vegetable HM pollution in Chongqing (56, 61). Sensitivity analysis can help identify key exposure parameters, and when combined with Monte Carlo simulation, it can effectively assess the health risks of HM pollution in different population groups.
5 Conclusions
This study provides a systematic evaluation of HM pollution in the fragmented farmland soil-crop system in Wanzhou District, Chongqing City. The results show that grains have a greater capacity to accumulate HMs compared to other crops. In the soil-crop system, over 70% of the samples had IICQ values indicating light pollution levels or higher. Both children and adults are exposed to high levels of non-carcinogenic and carcinogenic risks. As was the chief contributor to non-carcinogenic risk, while Ni, Cr, Cd, and As significantly contributed to the carcinogenic risk. Considering the concentration of HMs, pollution risks, and health hazards, As, Ni, Cr, and Cd should be listed for enhanced monitoring. Finally, we recommend strengthening the monitoring of HMs in crops, adopting proactive and effective regulatory strategies, and reasonably adjusting the dietary intake of residents. These measures are crucial for reducing the hazard of HMs to public health.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
LT: Writing – original draft, Formal analysis, Methodology, Funding acquisition, Visualization. AY: Software, Writing – original draft, Data curation, Formal analysis. MF: Formal analysis, Writing – original draft, Data curation. DZ: Supervision, Writing – original draft, Methodology. JL: Data curation, Investigation, Writing – original draft. HL: Supervision, Investigation, Writing – original draft. CL: Resources, Formal analysis, Investigation, Writing – original draft. YT: Formal analysis, Investigation, Writing – original draft, Resources. LC: Writing – review & editing, Investigation, Funding acquisition.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Science and Technology Research Program of Chongqing Municipal Education Commission (Grant No. KJQN202202709), Wanzhou District, Chongqing City Science and Health Joint Medical Research Project (Grant No. wzstc-kw2022018) and Chongqing public health key specialty (discipline) project Support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Gen AI was used in the creation of this manuscript.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1637595/full#supplementary-material
References
1. Du Y, Chen L, Ding P, Liu L, He Q, Chen B, et al. Different exposure profile of heavy metal and health risk between residents near a Pb-Zn mine and a Mn mine in Huayuan county, South China. Chemosphere. (2019) 216:352–64. doi: 10.1016/j.chemosphere.2018.10.142
2. Fei X, Lou Z, Lv X, Ren Z, Xiao R. Pollution threshold assessment and risk area delineation of heavy metals in soils through the finite mixture distribution model and Bayesian maximum entropy theory. J Hazard Mater. (2023) 452:131231. doi: 10.1016/j.jhazmat.2023.131231
3. Yang Q, Li Z, Lu X, Duan Q, Huang L, Bi J, et al. A review of soil heavy metal pollution from industrial and agricultural regions in China: pollution and risk assessment. Sci Total Environ. (2018) 642:690–700. doi: 10.1016/j.scitotenv.2018.06.068
4. Kumar S, Prasad S, Yadav KK, Shrivastava M, Gupta N, Nagar S, et al. Hazardous heavy metals contamination of vegetables and food chain: role of sustainable remediation approaches - a review. Environ Res. (2019) 179:108792. doi: 10.1016/j.envres.2019.108792
5. Tóth Z, Dombos M, Hornung E. Urban soil quality deteriorates even with low heavy metal levels: an arthropod-based multi-indices approach. Ecol Appl. (2023) 33:e2848. doi: 10.1002/eap.2848
6. Asaiduli H, Abliz A, Abulizi A, Sun X, Ye P. Assessment of soil heavy metal pollution and health risks in different functional areas on the northern slope of the eastern Tianshan mountains in Xinjiang, NW China. Int J Environ Res Public Health. (2023) 20:4843. doi: 10.3390/ijerph20064843
7. Chen L, Chang N, Liu H, Li T, Huang F, Zeng Y, et al. Multiomics provides insights into the impacts of microplastics on heavy metal(loid) accumulation in lettuce under simulated acid precipitation. ACS Sustain Chem Eng. (2025) 13:7667–81. doi: 10.1021/acssuschemeng.5c03114
8. Zheng F, Guo X, Tang M, Zhu D, Wang H, Yang X, et al. Variation in pollution status, sources, and risks of soil heavy metals in regions with different levels of urbanization. Sci Total Environ. (2023) 866:161355. doi: 10.1016/j.scitotenv.2022.161355
9. Chen X-X, Liu Y-M, Zhao Q-Y, Cao W-Q, Chen X-P, Zou C-Q. Health risk assessment associated with heavy metal accumulation in wheat after long-term phosphorus fertilizer application. Environ Pollut. (2020) 262:114348. doi: 10.1016/j.envpol.2020.114348
10. Zhang X, Wang X, Cao X, Xiao G, Miao H. Heavy element contents of vegetables and health-risk assessment in China. Sci Total Environ. (2022) 828:154552. doi: 10.1016/j.scitotenv.2022.154552
11. Wang B, Gao F, Qin N, Duan X, Li Y, Cao S, et al. A comprehensive analysis on source-distribution-bioaccumulation-exposure risk of metal(loid)s in various vegetables in peri-urban areas of Shenzhen, China. Environ Pollut. (2022) 293:118613. doi: 10.1016/j.envpol.2021.118613
12. Lee N-W, Wang H-Y, Du C-L, Yuan T-H, Chen C-Y, Yu C-J, et al. Air-polluted environmental heavy metal exposure increase lung cancer incidence and mortality: a population-based longitudinal cohort study. Sci Total Environ. (2022) 810:152186. doi: 10.1016/j.scitotenv.2021.152186
13. Rehman K, Fatima F, Waheed I, Akash MSH. Prevalence of exposure of heavy metals and their impact on health consequences. J Cell Biochem. (2018) 119:157–84. doi: 10.1002/jcb.26234
14. Wang M, Yu Q. Association between blood heavy metal concentrations and skin cancer in the national health and nutrition examination survey, 2011–2018. Environ Sci Pollut Res. (2023) 30:108681–93. doi: 10.1007/s11356-023-29674-4
15. Khatoon N, Ali S, Hussain A, Huang J, Li Z, Yu Z, et al. Human health risks of heavy metals in water and vegetables: a deterministic and probabilistic study in Gilgit-Baltistan, Pakistan. Environ Geochem Health. (2025) 47:88. doi: 10.1007/s10653-025-02389-w
16. Khatoon N, Ali S, Hussain A, Huang J, Yu Z, Liu H. Evaluating the carcinogenic and non-carcinogenic health risks of heavy metals contamination in drinking water, vegetables, and soil from Gilgit-Baltistan, Pakistan. Toxics. (2024) 13:5. doi: 10.3390/toxics13010005
17. Khatoon N, Ali S, Hussain A, Huang J, Yu Z, Liu H. Human health risks assessment of toxic metals via water, food, and soil: a case study of northern areas (Ghizer and Gilgit) of Pakistan. Results Eng. (2024) 24:103333. doi: 10.1016/j.rineng.2024.103333
18. Ngwese SN, Mouri H, Akoachere RAI, McKinley JM, Candeias C. Sources, pathways, and health risk assessment of harmful elements in soil and crops: a case study from east Cameroon. Expo Health. (2025). doi: 10.1007/s12403-025-00714-3
19. Thakkar D, Valand M, Vachhrajani K. Assessment of seasonal variations in soil heavy metal concentrations and potential health risks in Gujarat, India. Environ Geochem Health. (2024) 46:391. doi: 10.1007/s10653-024-02170-5
20. Singh G, Patel N, Jindal T, Ranjan MR. Heavy metal contamination in soils and crops irrigated by Kali river in Uttar Pradesh, India. Bull Environ Contam Toxicol. (2021) 107:931–7. doi: 10.1007/s00128-021-03349-7
21. Islam MM, Karim MdR, Zheng X, Li X. Heavy metal and metalloid pollution of soil, water and foods in Bangladesh: a critical review. Int J Environ Res Public Health. (2018) 15:2825. doi: 10.3390/ijerph15122825
22. Hasan GMMA, Das AK, Satter MA. Accumulation of heavy metals in rice (Oryza sativa. L) grains cultivated in three major industrial areas of Bangladesh. J Environ Public Health. (2022) 2022:1836597. doi: 10.1155/2022/1836597
23. He M, Shen H, Li Z, Wang L, Wang F, Zhao K, et al. Ten-year regional monitoring of soil-rice grain contamination by heavy metals with implications for target remediation and food safety. Environ Pollut. (2019) 244:431–9. doi: 10.1016/j.envpol.2018.10.070
24. Bi C, Zhou Y, Chen Z, Jia J, Bao X. Heavy metals and lead isotopes in soils, road dust and leafy vegetables and health risks via vegetable consumption in the industrial areas of Shanghai, China. Sci Total Environ. (2018) 619–620:1349–57. doi: 10.1016/j.scitotenv.2017.11.177
25. Hu N-W, Yu H-W, Deng B-L, Hu B, Zhu G-P, Yang X-T, et al. Levels of heavy metal in soil and vegetable and associated health risk in peri-urban areas across China. Ecotoxicol Environ Saf. (2023) 259:115037. doi: 10.1016/j.ecoenv.2023.115037
26. Liu X, Song Q, Tang Y, Li W, Xu J, Wu J, et al. Human health risk assessment of heavy metals in soil–vegetable system: a multi-medium analysis. Sci Total Environ. (2013) 463–464:530–40. doi: 10.1016/j.scitotenv.2013.06.064
27. Zhang J, Yang R, Li YC, Peng Y, Wen X, Ni X. Distribution, accumulation, and potential risks of heavy metals in soil and tea leaves from geologically different plantations. Ecotoxicol Environ Saf. (2020) 195:110475. doi: 10.1016/j.ecoenv.2020.110475
28. Wang Y, Xu W, Li J, Song Y, Hua M, Li W, et al. Assessing the fractionation and bioavailability of heavy metals in soil–rice system and the associated health risk. Environ Geochem Health. (2022) 44:301–18. doi: 10.1007/s10653-021-00876-4
29. Wang Z, Chen Y, Wang S, Yu Y, Huang W, Xu Q, et al. Pollution risk assessment and sources analysis of heavy metal in soil from bamboo shoots. Int J Environ Res Public Health. (2022) 19:14806. doi: 10.3390/ijerph192214806
30. Chen R, Han L, Liu Z, Zhao Y, Li R, Xia L, et al. Assessment of soil-heavy metal pollution and the health risks in a mining area from southern Shaanxi province, China. Toxics. (2022) 10:385. doi: 10.3390/toxics10070385
31. Cai L-M, Wang Q-S, Wen H-H, Luo J, Wang S. Heavy metals in agricultural soils from a typical township in Guangdong Province, China: occurrences and spatial distribution. Ecotoxicol Environ Saf. (2019) 168:184–91. doi: 10.1016/j.ecoenv.2018.10.092
32. Dong B, Zhang R, Gan Y, Cai L, Freidenreich A, Wang K, et al. Multiple methods for the identification of heavy metal sources in cropland soils from a resource-based region. Sci Total Environ. (2019) 651:3127–38. doi: 10.1016/j.scitotenv.2018.10.130
33. Wang Y, Liu C, Zhou D, Cheng H. A new approach for evaluating soil heavy metal impact: a comprehensive index combined soil environmental quality and agricultural products quality. J Agro-Environ Sci. (2016) 35:1225–32. doi: 10.11654/jaes.2016.07.001
34. Chen R, Cai X, Ding G, Ren F, Wang Q, Cheng N, et al. Ecological risk assessment of heavy metals in farmland soils in Beijing by three improved risk assessment methods. Environ Sci Pollut Res. (2021) 28:57970–82. doi: 10.1007/s11356-021-14695-8
35. Liu S, Zhao H, Wu K, Zhang Z, Hou Y, Chen T, et al. Evaluation of heavy metal distribution characteristics of agricultural soil-rice system in a high geological background area according to the influence index of comprehensive quality (IICQ). Environ Sci Pollut Res. (2020) 27:20920–33. doi: 10.1007/s11356-020-08453-5
36. Yu Y, Yang L, Yang Y, Chen H, Bo L, Huang B, et al. A collaborative approach for metal pollution assessment in production system of plastic-shed vegetables near industrial areas. Bull Environ Contam Toxicol. (2024) 113:7. doi: 10.1007/s00128-024-03916-8
37. Zhang M, Cheng L, Yue Z, Peng L, Xiao L. Assessment of heavy metal(oid) pollution and related health risks in agricultural soils surrounding a coal gangue dump from an abandoned coal mine in Chongqing, southwest China. Sci Rep. (2024) 14:18667. doi: 10.1038/s41598-024-69072-5
38. Liang A, Wang Y, Guo H, Bo L, Zhang S, Bai Y. Assessment of pollution and identification of sources of heavy metals in the sediments of Changshou lake in a branch of the Three Gorges reservoir. Environ Sci Pollut Res. (2015) 22:16067–76. doi: 10.1007/s11356-015-4825-8
39. Wang T, Pan J, Liu X. Characterization of heavy metal contamination in the soil and sediment of the Three Gorges reservoir, China. J Environ Sci Health Part A. (2017) 52:201–9. doi: 10.1080/10934529.2016.1246931
40. Liu J, Li X, Zhang P, Zhu Q, Lu W, Yang Y, et al. Contamination levels of and potential risks from metal(loid)s in soil-crop systems in high geological background areas. Sci Total Environ. (2023) 881:163405. doi: 10.1016/j.scitotenv.2023.163405
41. State Environmental Protection Administration of China. The Technical Specification for Soil Environmental Monitoring, HJ/T 166-2004 (2004) (in Chinese).
42. Cheng JW, Cai SW, Huang MQ, Yan X. The migration, transformation, and risk assessment of heavy metal forms in the rhizosphere soil of plants in the Zunyi manganese mine area, Guizhou. Environ Chem. (2022) 41:2833–41. doi: 10.7524/j.issn.0254-6108.2021050606
43. Marković J, Jović M, Smičiklas I, Pezo L, Šljivić-Ivanović M, Onjia A, et al. Chemical speciation of metals in unpolluted soils of different types: correlation with soil characteristics and an ANN modelling approach. J Geochem Explor. (2016) 165:71–80. doi: 10.1016/j.gexplo.2016.03.004
44. Miletić A, Vesković J, Lučić M, Onjia A. Monte Carlo simulation of source-specific risks of soil at an abandoned lead-acid battery recycling site. Stoch Environ Res Risk Assess. (2024) 38:3313–29. doi: 10.1007/s00477-024-02747-7
45. Onjia A, Huang X, Trujillo González JM, Egbueri JC. Editorial: chemometric approach to distribution, source apportionment, ecological and health risk of trace pollutants. Front Environ Sci. (2022) 10:1107465. doi: 10.3389/fenvs.2022.1107465
46. Ministry of Agriculture of the People's Republic of China. Determination of pH in Soil, NY/T 1377–2007. (2007) (in Chinese).
47. Ministry of Environmental Protection of the People's Republic of China. Oil and Sediment-Determination of Total Mercury-Catalytic Pyrolysis-Cold Atomic Absorption Spectrophotometry, HJ 923–2017 (2017) (in Chinese).
48. Ministry of Environmental Protection of the People's Republic of China. Soil and Sediment-Determination of Aqua Regia Extracts of 12 Metal Elements-Inductively Coupled Plasma Mass Spectrometry, HJ 803–2016. (2016) (in Chinese).
49. National Health Commission of the People's Republic of China and State Administration for Market Regulation. Determination of Total Mercury and Organic Mercury in Foods, GB 5009.17–2021. (2021) (in Chinese).
50. National Health and Family Planning Commission of the People's Republic of China and the China Food and Drug Administration. Determination of Multiple Elements in Food, GB 5009.268–2016. (2016) (in Chinese).
51. Muller G. Index of geoaccumulation in sediments of the Rhine river. GeoJournal. (1969) 2:108–18.
52. Muller G. Schwermetallbelstung Der Sedimente Des Neckars Und Seiner Nebenflusse: Eine Estandsaufnahme. Germany: Chemiker-Zeitung (1981).
53. Wan M, Hu W, Wang H, Tian K, Huang B. Comprehensive assessment of heavy metal risk in soil-crop systems along the Yangtze river in Nanjing, southeast China. Sci Total Environ. (2021) 780:146567. doi: 10.1016/j.scitotenv.2021.146567
54. Drozdova I, Alekseeva-Popova N, Dorofeyev V, Bech J, Belyaeva A, Roca N, et al. A comparative study of the accumulation of trace elements in Brassicaceae plant species with phytoremediation potential. Appl Geochem. (2019) 108:104377. doi: 10.1016/j.apgeochem.2019.104377
55. Gan Y, Huang X, Li S, Liu N, Li YC, Freidenreich A, et al. Source quantification and potential risk of mercury, cadmium, arsenic, lead, and chromium in farmland soils of Yellow river delta. J Clean Prod. (2019) 221:98–107. doi: 10.1016/j.jclepro.2019.02.157
56. Lü Q, Xiao Q, Guo Y, Wang Y, Cai L, You W, et al. Pollution monitoring, risk assessment and target remediation of heavy metals in rice from a five-year investigation in western Fujian region, China. J Hazard Mater. (2022) 424:127551. doi: 10.1016/j.jhazmat.2021.127551
57. Hu B, Shao S, Fu T, Fu Z, Zhou Y, Li Y, et al. Composite assessment of human health risk from potentially toxic elements through multiple exposure routes: a case study in farmland in an important industrial city in East China. J Geochem Explor. (2020) 210:106443. doi: 10.1016/j.gexplo.2019.106443
58. Miletić A, Lučić M, Onjia A. Exposure factors in health risk assessment of heavy metal(loid)s in soil and sediment. Metals. (2023) 13:1266. doi: 10.3390/met13071266
59. Bao L, Deng H, Jia Z, Li Y, Dong J, Yan M, et al. Ecological and health risk assessment of heavy metals in farmland soil of northwest Xiushan, Chongqing. Geol China. (2020) 47:1625–36. doi: 10.12029/gc20200602
60. Ministry of Ecology and Environment of the People's Republic of China and State Administration for Market Regulation. Soil Environmental Quality Risk Control Standard for Soil Contamination of Agricultural Land, GB 15618–2018. (2018) (in Chinese).
61. Shi J, Yang Y, Shen Z, Lin Y, Mei N, Luo C, et al. Identifying heavy metal sources and health risks in soil-vegetable systems of fragmented vegetable fields based on machine learning, positive matrix factorization model and Monte Carlo simulation. J Hazard Mater. (2024) 478:135481. doi: 10.1016/j.jhazmat.2024.135481
62. Fei X, Lou Z, Xiao R, Ren Z, Lv X. Source analysis and source-oriented risk assessment of heavy metal pollution in agricultural soils of different cultivated land qualities. J Clean Prod. (2022) 341:130942. doi: 10.1016/j.jclepro.2022.130942
63. National Health Commission of the People's Republic of China and State Administration for Market Regulation. Contaminant Limits in Food, GB 2762–2022. (2022) (in Chinese).
64. Ministry of Health of the People's Republic of China. Tolerance Limit of Copper in Foods, GB 15199–1994. (1994) (in Chinese).
65. Ministry of Health of the People's Republic of China. Tolerance Limit of Zinc in Foods, GB 13106–1991. (1991) (in Chinese).
66. Zhou Y, Wang Q, Wang C, Liu F, Song Y, Guo Z, et al. Heavy metal pollution and cumulative effect of soil-crop systems around typical enterprises in Xiong'an new district. Environ Sci. (2021) 42:5977–87. doi: 10.13227/j.hjkx.202103255
67. Liu H, Wang H, Zhou J, Zhang Y, Wang H, Li M, et al. Environmental cadmium pollution and health risk assessment in rice–wheat rotation area around a smelter. Environ Sci Pollut Res. (2024) 31:433–44. doi: 10.1007/s11356-023-31215-y
68. Yuan X, Xue N, Han Z. A meta-analysis of heavy metals pollution in farmland and urban soils in China over the past 20 years. J Environ Sci. (2021) 101:217–26. doi: 10.1016/j.jes.2020.08.013
69. Luo Y, Wang Z, Zhang Z-L, Zhang J-Q, Zeng Q-P, Tian D, et al. Contamination characteristics and source analysis of potentially toxic elements in dustfall-soil-crop systems near non-ferrous mining areas of Yunnan, southwestern China. Sci Total Environ. (2023) 882:163575. doi: 10.1016/j.scitotenv.2023.163575
70. Xie N, Kang C, Ren D, Zhang L. Assessment of the variation of heavy metal pollutants in soil and crop plants through field and laboratory tests. Sci Total Environ. (2022) 811:152343. doi: 10.1016/j.scitotenv.2021.152343
71. Wu D, Liu H, Wu J, Gao X, Nyasha NK, Cai G, et al. Bi-directional pollution characteristics and ecological health risk assessment of heavy metals in soil and crops in Wanjiang economic zone, Anhui province, China. Int J Environ Res Public Health. (2022) 19:9669. doi: 10.3390/ijerph19159669
72. Mabood F, Hadi F, Jan AU, Ditta A, Islam Z, Siddiqui MH, et al. Assessment of Pb and Ni and potential health risks associated with the consumption of vegetables grown on the roadside soils in district Swat, Khyber Pakhtunkhwa, Pakistan. Environ Monit Assess. (2022) 194:906. doi: 10.1007/s10661-022-10627-2
73. Wang Z, Bao J, Wang T, Moryani HT, Kang W, Zheng J, et al. Hazardous heavy metals accumulation and health risk assessment of different vegetable species in contaminated soils from a typical mining city, central China. Int J Environ Res Public Health. (2021) 18:2617. doi: 10.3390/ijerph18052617
74. Liu Y, Zhang R, Pan B, Qiu H, Wang J, Zhang J, et al. Uptake of heavy metals by crops near a mining field: pathways from roots and leaves. Chemosphere. (2023) 322:138215. doi: 10.1016/j.chemosphere.2023.138215
75. Islam EU, Yang X, He Z, Mahmood Q. Assessing potential dietary toxicity of heavy metals in selected vegetables and food crops. J Zhejiang Univ Sci B. (2007) 8:1–13. doi: 10.1631/jzus.2007.B0001
76. Zhu S, Li S, Shao M, Qiao J, Zhu Y. Source identification and migration fate of heavy metals of soil-groundwater system in a thousand-year cultivation region. Environ Geochem Health. (2024) 46:345. doi: 10.1007/s10653-024-02130-z
77. Hu Y, Huang Y, Liu Y. Genotypic differences in iron plaque formation on root surface and their effect on arsenic uptake and translocation in rice (Oryza sativa L) grown in arsenic contaminated soil. Asian J Ecotoxicol. (2013) 8:923–30. doi: 10.7524/AJE.1673-5897.20121102001
78. Jia Y, Huang H, Chen Z, Zhu Y-G. Arsenic uptake by rice is influenced by microbe-mediated arsenic redox changes in the rhizosphere. Environ Sci Technol. (2014) 48:1001–7. doi: 10.1021/es403877s
79. Wang C-C, Zhang Q-C, Yan C-A, Tang G-Y, Zhang M-Y, Ma LQ, et al. Heavy metal(loid)s in agriculture soils, rice, and wheat across China: status assessment and spatiotemporal analysis. Sci Total Environ. (2023) 882:163361. doi: 10.1016/j.scitotenv.2023.163361
80. Lin H, Wang Z, Liu C, Dong Y. Technologies for removing heavy metal from contaminated soils on farmland: a review. Chemosphere. (2022) 305:135457. doi: 10.1016/j.chemosphere.2022.135457
81. Liu X, Gu S, Yang S, Deng J, Xu J. Heavy metals in soil-vegetable system around E-waste site and the health risk assessment. Sci Total Environ. (2021) 779:146438. doi: 10.1016/j.scitotenv.2021.146438
82. Khoshakhlagh AH, Mohammadzadeh M, Bamel U, Gruszecka-Kosowska A. Human exposure to heavy metals and related cancer development: a bibliometric analysis. Environ Sci Pollut Res. (2023) 30:109867–88. doi: 10.1007/s11356-023-29939-y
83. Rzymski P, Tomczyk K, Rzymski P, Poniedziałek B, Opala T, Wilczak M. Impact of heavy metals on the female reproductive system. Ann Agric Environ Med. (2015) 22:259–64. doi: 10.5604/12321966.1152077
84. Pan Z, Gong T, Liang P. Heavy metal exposure and cardiovascular disease. Circ Res. (2024) 134:1160–78. doi: 10.1161/CIRCRESAHA.123.323617
85. Mwelwa S, Chungu D, Tailoka F, Beesigamukama D, Tanga C. Biotransfer of heavy metals along the soil-plant-edible insect-human food chain in Africa. Sci Total Environ. (2023) 881:163150. doi: 10.1016/j.scitotenv.2023.163150
86. Yang S, Zhou Q, Sun L, Sun Y, Qin Q, Song K, et al. A prospective health risks analysis of regulatory limits for heavy metals in rice from representative organizations and countries worldwide: are they protective? Sci Total Environ. (2023) 904:167130. doi: 10.1016/j.scitotenv.2023.167130
87. Zeng X, Xu X, Qin Q, Ye K, Wu W, Huo X. Heavy metal exposure has adverse effects on the growth and development of preschool children. Environ Geochem Health. (2019) 41:309–21. doi: 10.1007/s10653-018-0114-z
88. Peng J, Zhang S, Han Y, Bate B, Ke H, Chen Y. Soil heavy metal pollution of industrial legacies in China and health risk assessment. Sci Total Environ. (2022) 816:151632. doi: 10.1016/j.scitotenv.2021.151632
89. Yang X, Cheng B, Gao Y, Zhang H, Liu L. Heavy metal contamination assessment and probabilistic health risks in soil and maize near coal mines. Front Public Health. (2022) 10:1004579. doi: 10.3389/fpubh.2022.1004579
90. Sanaei F, Amin MM, Alavijeh ZP, Esfahani RA, Sadeghi M, Bandarrig NS, et al. Health risk assessment of potentially toxic elements intake via food crops consumption: Monte Carlo simulation-based probabilistic and heavy metal pollution index. Environ Sci Pollut Res. (2021) 28:1479–90. doi: 10.1007/s11356-020-10450-7
91. Wang J, Su J, Li Z, Liu B, Cheng G, Jiang Y, et al. Source apportionment of heavy metal and their health risks in soil-dustfall-plant system nearby a typical non-ferrous metal mining area of Tongling, eastern China. Environ Pollut. (2019) 254:113089. doi: 10.1016/j.envpol.2019.113089
92. Yang S, He M, Zhi Y, Chang SX, Gu B, Liu X, et al. An integrated analysis on source-exposure risk of heavy metals in agricultural soils near intense electronic waste recycling activities. Environ Int. (2019) 133:105239. doi: 10.1016/j.envint.2019.105239
93. Huang Y-N, Dang F, Li M, Zhou D-M, Song Y, Wang J-B. Environmental and human health risks from metal exposures nearby a Pb-Zn-Ag mine, China. Sci Total Environ. (2020) 698:134326. doi: 10.1016/j.scitotenv.2019.134326
94. Chen S, Gao Y, Wang C, Gu H, Sun M, Dang Y, et al. Heavy metal pollution status, children health risk assessment and source apportionment in farmland soils in a typical polluted area, northwest China. Stoch Environ Res Risk Assess. (2024) 38:2383–95. doi: 10.1007/s00477-024-02685-4
95. Yu P, Han Y, Wang M, Zhu Z, Tong Z, Shao X, et al. Heavy metal content and health risk assessment of atmospheric particles in China: a meta-analysis. Sci Total Environ. (2023) 867:161556. doi: 10.1016/j.scitotenv.2023.161556
96. Cui W, Mei Y, Liu S, Zhang X. Health risk assessment of heavy metal pollution and its sources in agricultural soils near Hongfeng lake in the mining area of Guizhou province, China. Front Public Health. (2023) 11:1276925. doi: 10.3389/fpubh.2023.1276925
97. Kan X, Dong Y, Feng L, Zhou M, Hou H. Contamination and health risk assessment of heavy metals in China's lead–zinc mine tailings: a meta–analysis. Chemosphere. (2021) 267:128909. doi: 10.1016/j.chemosphere.2020.128909
98. Singh PK, Shikha D, Saw S. Evaluation of potential toxic heavy metal contamination in soil, fly ash, vegetables and grain crops along with associated ecological and health risk assessment of nearby inhabitants of a thermal power station in Jharkhand (India). Environ Sci Pollut Res. (2023) 30:7752–69. doi: 10.1007/s11356-022-22638-0
99. Li X, Li Z, Lin C-J, Bi X, Liu J, Feng X, et al. Health risks of heavy metal exposure through vegetable consumption near a large-scale Pb/Zn smelter in central China. Ecotoxicol Environ Saf. (2018) 161:99–110. doi: 10.1016/j.ecoenv.2018.05.080
100. Tang L, Deng S, Tan D, Long J, Lei M. Heavy metal distribution, translocation, and human health risk assessment in the soil-rice system around Dongting lake area, China. Environ Sci Pollut Res. (2019) 26:17655–65. doi: 10.1007/s11356-019-05134-w
101. Wan F, Jiang N, Yu L, Zang K, Liu S, He W, et al. Heavy metal ecological-health risk assessment under wheat–maize rotation system in a high geological background area in eastern China. Sci Rep. (2022) 12:17912. doi: 10.1038/s41598-022-22608-z
102. Cai L-M, Wang Q-S, Luo J, Chen L-G, Zhu R-L, Wang S, et al. Heavy metal contamination and health risk assessment for children near a large Cu-smelter in central China. Sci Total Environ. (2019) 650:725–33. doi: 10.1016/j.scitotenv.2018.09.081
103. Wang Y, Cheng H. Soil heavy metal(loid) pollution and health risk assessment of farmlands developed on two different terrains on the Tibetan plateau, China. Chemosphere. (2023) 335:139148. doi: 10.1016/j.chemosphere.2023.139148
104. Tóth G, Hermann T, Da Silva MR, Montanarella L. Heavy metals in agricultural soils of the European union with implications for food safety. Environ Int. (2016) 88:299–309. doi: 10.1016/j.envint.2015.12.017
105. Orosun MM. Assessment of arsenic and its associated health risks due to mining activities in parts of north-central Nigeria: probabilistic approach using Monte Carlo. J Hazard Mater. (2021) 412:125262. doi: 10.1016/j.jhazmat.2021.125262
106. Jalali M, Karimi Mojahed J. Assessment of the health risks of heavy metals in soils and vegetables from greenhouse production systems in Iran. Int J Phytoremediation. (2020) 22:834–48. doi: 10.1080/15226514.2020.1715917
Keywords: heavy metals, crops, IICQ index, Monte Carlo simulation, health risks
Citation: Tong L, Yang A, Fan M, Zheng D, Li J, Liu H, Li C, Tang Y and Cheng L (2025) Assessment of heavy metal sources and health risks in soil-crop systems of fragmented farmland. Front. Public Health 13:1637595. doi: 10.3389/fpubh.2025.1637595
Received: 29 May 2025; Accepted: 08 July 2025;
Published: 31 July 2025.
Edited by:
Mohsen Saeedi, University Canada West, CanadaReviewed by:
Antonije Onjia, University of Belgrade, SerbiaLi Chen, Lanzhou University, China
Nafeesa Khatoon, Zhengzhou University, China
Copyright © 2025 Tong, Yang, Fan, Zheng, Li, Liu, Li, Tang and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Longfeng Cheng, aGFpbmFzdGlva291QGhvdG1haWwuY29t
†These authors have contributed equally to this work and share first authorship