Next Article in Journal
Defining Post-COVID Symptoms (Post-Acute COVID, Long COVID, Persistent Post-COVID): An Integrative Classification
Next Article in Special Issue
Impact of Environmental Airborne Manganese Exposure on Cognitive and Motor Functions in Adults: A Systematic Review and Meta-Analysis
Previous Article in Journal
Psychosocial Effects and Use of Communication Technologies during Home Confinement in the First Wave of the COVID-19 Pandemic in Italy and The Netherlands
Previous Article in Special Issue
A Positive Relationship between Exposure to Heavy Metals and Development of Chronic Diseases: A Case Study from Chile
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hazardous Heavy Metals Accumulation and Health Risk Assessment of Different Vegetable Species in Contaminated Soils from a Typical Mining City, Central China

1
School of Environmental Studies, China University of Geosciences, Wuhan 430074, China
2
Hubei Key Laboratory of Mine Environment Pollution Control and Remediation, School of Environmental Science and Engineering, Hubei Polytechnic University, Huangshi 435003, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(5), 2617; https://doi.org/10.3390/ijerph18052617
Submission received: 19 January 2021 / Revised: 22 February 2021 / Accepted: 24 February 2021 / Published: 5 March 2021
(This article belongs to the Special Issue Metal Exposure and Health Risk Assessment 2020)

Abstract

:
Heavy metal poisoning has caused serious and widespread human tragedies via the food chain. To alleviate heavy metal pollution, particular attention should be paid to low accumulating vegetables and crops. In this study, the concentrations of five hazardous heavy metals (HMs), including copper (Cu), chromium (Cr), lead (Pb), cadmium (Cd), and arsenic (As) were determined from soils, vegetables, and crops near four typical mining and smelting zones. Nemerow’s synthetical pollution index (Pn), Potential ecological risk index (RI), and Geo-accumulation index (Igeo) were used to characterize the pollution degrees. The results showed that soils near mining and metal smelting zones were heavily polluted by Cu, Cd, As, and Pb. The total excessive rate followed a decreasing order of Cd (80.00%) > Cu (61.11%) > As (45.56%) > Pb (32.22%) > Cr (0.00%). Moreover, sources identification indicated that Cu, Pb, Cd, and As may originate from anthropogenic activities, while Cr may originate from parent materials. The exceeding rates of Cu, Cr, Pb, Cd, and As were 6.7%, 6.7%, 66.7%, 80.0%, and 26.7% among the vegetable and crop species, respectively. Particularly, vegetables like tomatoes, bell peppers, white radishes, and asparagus, revealed low accumulation characteristics. In addition, the hazard index (HI) for vegetables and crops of four zones was greater than 1, revealing a higher risk to the health of local children near the mine and smelter. However, the solanaceous fruit has a low-risk index (HI), indicating that it is a potentially safe vegetable type.

1. Introduction

Heavy metal (HM) contamination, including copper (Cu), chromium (Cr), lead (Pb), cadmium (Cd), and arsenic (As), is primarily attributed to anthropogenic activities, such as mining and metal smelting, traffic and factory emissions, ore exploitation dust fall, and sewage irrigation [1,2,3,4,5]. The soil quality of arable and cultivated lands near industrial and abandoned multiple metal mine is concerning. As reported by the Ministry of Environment Protection (MEP, China, 2014) [6], about 16.1% of total soils, including 19.4% of the arable lands and 36.3% of enterprise surrounding soils was over the Grade II environment quality standard (GB15618-1995) [7]. According to previous studies, agricultural soils were mainly polluted by anthropogenic activities such as mining exploitation, metal smelting, industrial processes, traffic emissions, and agricultural pesticides and fertilizers [8,9,10,11]. The soil-plant-human route is a dominant pathway for HMs in vegetables and crops. So far, many studies have shown that plant roots can uptake and accumulate these contaminants and gradually transfer them to the edible parts, ultimately posing a potential risk to human health through the food chain [12,13,14,15]. Besides, direct foliar uptake has been also considered as the direct pathway of atmospheric HMs such as Pb, Hg, and partly Cd which can be absorbed by stomata pores of leaves [16]. Therefore, the direct consumption of leaf vegetables can cause a severe threat to the health of inhabitants [17,18]. Some essential nutrients stored in the body are gradually depleted due to the consumption of foods containing HMs. Chronic exposure to Pb, Cd, and As can lead to various disorders, such as nervous system symptoms, kidney dysfunction and liver damage, lung and skin cancer, and bone fractures [3,19].
It is noteworthy that, currently, the Chinese government is paying considerable attention to food security, especially the concentration of HMs in vegetables and crops [12,20]. Many studies have investigated the heavy metal pollution of vegetables near tungsten mining in Gannan of Jiangxi [21]; Zhuzhou smelter of Hunan [22]; Jijie Pb/Zn smelter of Yunnan [23]; Tonglushan copper ming of Hubei in China [3], and other countries [24,25,26,27]. These findings revealed that most of the vegetables in the exposed areas exceed the maximum permissible concentration, and even caused a higher health risk than those in non-contaminated areas. It was reported that not only could the vegetables uptake amounts of HMs, but also elevated contents of HMs were found in food crops such as rice [12,28], wheat [29], and maize [30]. Particularly, Pb and Cd exhibited greater transfer capacity and ultimately accumulated in the grain [16,18]. However, most of the studies only considered a limited number of species in vegetables or certain hazardous metals. Little attention has been paid to the health risk of dietary intake to multiple hazardous metals. In addition, researchers have gradually drawn attention to how to limit HMs content in edible parts [31,32,33]. To move this issue, a new collection of low accumulating HMs species has been proposed. This ensures a reasonably low concentration of HMs in edible parts of vegetables or crop plants, even though they are cultivated on contaminated soil. Thus, this is one available pathway to reduce health risk, which may be due to the accumulation variation of different species and genotypes [33].
Daye City has a mining history of over 3000 years and has played a significant role in the national economic development of China as one of the largest copper producing areas. According to the National Ministry of Agriculture’s Environmental Protection Research and Monitoring Institute, heavy metal contamination of arable soils was extremely serious and mostly spread in villages and towns including Luojiaqiao, Jinhu, Dajipu, Chengui, and Jinshandian [34]. Previous studies have shown that Daye City’s soils, waters, and agricultural products have been polluted to higher degrees due to a vast number of industrial activities [3,34,35]. However, heavy metal exposure in agricultural soils typically occurs from a number of causes, not just mining activities. To date, studies to identify the sources of HMs mainly focused on a specific contaminated region or a small range of mining and smelting areas [16,21]. There are few studies on possible sources of HMs in Daye City [35]. Therefore, it is important to investigate some representative sites near the sensitive zones.
The main objectives of this study are to: (1) estimate heavy metals pollution levels and relevant impacts of soils, vegetables, and crops; (2) identify the sources of HMs in different study zones; (3) compare the multiple hazardous metals accumulation in vegetables and crops for the selection of low accumulating species; and (4) assess the dietary intakes and health risks of HMs in different food types to local children. The evaluated results will provide a scientific basis and theoretical support for environmental management in a typical mining city. Furthermore, the results also ensure human health and food safety and will provide valuable technical guidance to agricultural planting pattern.

2. Materials and Methods

2.1. Study Areas

The study areas (E 114°49′13″–114°57′9″, N 30°4′13″–30°13′29″) are located in Daye City which is situated in the middle reaches of the Yangtze River and in the southeast of Hubei Province, central China. The climate is the subtropical monsoon, and the average annual temperature is 17 °C. In addition, the average annual precipitation is 1382.6 mm with abundant rainfall. The four distinct seasons and superior water-heat conditions are conducive to the growth of crops.
Considering the relationship between heavy metals distributions and anthropogenic emissions, the study areas can be classified into four categories of vegetable fields near the non-ferrous metals smelter (A zone), abandoned copper mine (B zone), limestone quarry (C zone), and iron mine (D zone) (Figure 1). Moreover, they are located in four townships: Xialu District, Jinhu town, Huandiqiao town, and Jinshandian town. A zone is located in Xialu District, whose north part is adjacent to the open pit and concentrator. B zone is located in Jinhu town, which has a mining area of 7.8 km2 and a long history of open-cast mining. C zone is located at the eastern foot of Taiping Mountain of Huandiqiao town, which has been affected by quarrying and transportation for years. D zone is located in Jinshandian town, which is one of Wuhan Iron and Steel (Group) Company’s large underground iron mines.

2.2. Vegetable Samples and Rhizosphere Soils Collection

Vegetable and crop samples and rhizosphere soil samples were collected from April to October 2019. Ninety sampling sites were randomly selected near the mining areas. Each sampling site was recorded by using a global position system (GPS) (Figure 1). Ninety samples represented fifteen common species of vegetables, including four species of solanaceous fruits, five species of leafy vegetables, two species of rhizome geophytes, Legumes, and one species of nuts and cereals (Table 1), respectively. At each site, five sub-samples were collected and mixed to form an about 1.0 kg of the cross-sectional sample [27].
For each vegetable or crop sample at its corresponding site, we collected rhizosphere soil sample. In order to obtain a representative sample, five sub-samples were randomly collected at each site and bulked together to form a composite sample [3]. Approximately 500 g of soil samples were stored and transported to the laboratory immediately.

2.3. Sample Preparation and Analysis

The soil samples were air-dried for two weeks at room temperature, then grounded and sieved through a 100-mesh polyethylene sieve. All ground soil samples were stored in air-tight plastic bags before analysis. Each soil sample (~0.25 g) was prepared in a digestion tube with 10 mL of optimal levels of acid (HCl/HNO3/HF = 6/2/2). Then, they were heated at 190 °C for 30 min by using a microwave digestion oven (Mars5, CEM Corporation, Matthews, NC, USA) [36]. After cooling, the digested solution was filtered and diluted to 50 mL with 0.2% HNO3.
The vegetable samples were thoroughly washed by using tap water and ultrapure water and then weighed. After blanching at 105 °C for 30 min, samples were oven-dried at 60 °C to a constant weight for obtaining water percent, then ground to <1.0 mm by using a grinder. Each ash sample (0.5 g) was obtained using muffle furnace at 600 °C for 5 h, digested with 10 mL of ultrapure HNO3, cooled and diluted with deionized water to 50 mL [3].
Soil pH levels were determined in a water/soil suspension with a ratio of 2.5:1 (v/w) [36]. The concentrations of heavy metal in soil and vegetable samples were determined by using inductively coupled plasma optical emission spectrometry (ICP-OES 7300, Thermo Fisher Scientific, Waltham, MA, USA). In order to ensure analytical precision and accuracy, sample replicates, reagent blanks, and standard reference materials (GBW 08301 and GBW 10015) from the Chinese Academy of Measurement Sciences were carried out simultaneously in each batch. The element recovery rate of 89–108% showed better accuracy. The relative standard deviations (RSD) were less than 5%. The limit of detection (LOD) for Cu, Cr, Pb, Cd, and As are 0.05, 0.05, 0.02, 0.002, and 0.002 mg/kg, respectively.

2.4. Data Analysis

2.4.1. Nemerow’s Synthetical Pollution Index (Pn)

The level of soils contamination was assessed by the single-factor (Pi) and Nemerow’s synthetical pollution index (Pn) [37]. These indices were computed by Equations (1) and (2):
P i = C i / S i
where Pi is the single-factor index; Ci is the measured concentration of each heavy metal (mg/kg), and Si is the standard value of each element (mg/kg). In this study, the standard of GB15618-2018 was applied as agricultural soil assessment criteria.
P n = ( P imax 2 + P ¯ 2 ) / 2
where Pn is Nemerow’s synthetical pollution index, P ¯ is the average of single-factor index, P imax is the maximum of single-factor index, and Pn is classified as safety (Pn ≤ 0.7, Class I), precaution (0.7 < Pn ≤ 1.0, Class II), slightly polluted (1.0 < Pn ≤ 2.0, Class III), moderately polluted (2.0 < Pn ≤ 3.0, Class IV), and seriously polluted (Pn > 3.0, Class V).

2.4.2. Potential Ecological Risk Index (RI)

The potential ecological risk index (RI) was conducted to evaluate multiple HMs risk from soils, the index is a comprehensive consideration of HMs contents, environmental influence, and biotoxicity [38]. The RI is computed by Equations (3) and (4):
E r i = T r i × C r i = T r i × C s i / C n i
RI = i = 1 n E r i
where E r i represents the potential ecological risk index of a certain element, C s i is the actual concentration of each heavy metal, T r i is the toxic response coefficients and the values of Cu, Cr, Pb, Cd, and As are 5, 2, 5, 30 and 10, respectively. The RI classifications are low risk (RI < 150), moderate risk (150 ≤ RI < 300), considerable risk (300 ≤ RI < 600), and high risk (RI ≥ 600). The single element risk degrees ( E r i ) are low ecological risk ( E r i < 40), moderate ecological risk (40 ≤ E r i < 80), considerable ecological risk (80 ≤ E r i < 160), high ecological risk (160 ≤ E r i < 320), and very high ecological risk ( E r i ≥ 320).

2.4.3. Geo-Accumulation Index (Igeo)

The geo-accumulation index (Igeo) was applied to evaluate spatial distribution pollution levels of different heavy metals in the four research areas, the index was calculated using the following Equation (5) [39]:
I geo = log 2 ( C n / K · B n )
where Cn is the concentrations of heavy metals (mg/kg), Bn is the soil background values in Hubei province (Table 1), K is introduced as the constant 1.5 to minimize some possible influence of background values due to geogenic variation. The contamination degrees of Igeo are classified as unpolluted (Igeo ≤ 0), slightly to moderately polluted (0 < Igeo ≤ 1), moderately polluted (1 < Igeo ≤ 2), moderately to heavily polluted (2 < Igeo ≤ 3), heavily polluted (3 < Igeo ≤ 4), heavily to extremely polluted (4 < Igeo ≤ 5), and extremely polluted (Igeo > 5).

2.4.4. Health Risk Assessment

The health risk of heavy metals via vegetable consumption can be evaluated by target hazard quotients (THQ) and hazard index (HI) which was proposed by the US Environmental Protection Agency (USEPA 2000). The THQ was calculated using Equation (6) to assess the non-carcinogenic risks of individual heavy metal for local children. Moreover, The HI was computed to estimate the potential risk of multiple heavy metals by using the following Equation (7):
THQ = EF × ED × C i × F IR RfD × BW × AT
HI = n = 1 i THQ n ; i = 1 , 2 , 3 n
where Ci is the concentration of heavy metals in the edible part of vegetables (mg/kg, Fresh weight); FIR is the ratio of food-intake for which the value of local children was 94 g/d [40]; EF is exposure frequency (365 d/a); ED is exposure duration (76.4 years, life expectancy in China according to World Health Statistics 2018); BW represents the average of body weight (18.7 kg for an adolescent); AT represents the average exposure time for non-carcinogenic effects (assuming 76.4 years, 365 days a year); RfD is the daily reference dose of heavy metals (μg/(kg·d)), the values of RfD for Cu, Cr, Cd, As were 40, 3.0, 1.0, 0.3, respectively [41]. Although Pb element had not been required in USEPA, Pb could exert a bad influence on the central nervous system of adolescents. Thus, the RfD for Pb was 3.5 μg/(kg·d) in present study [22].
If the value of THQ or HI is ≤1.0, there is no obvious health risk to residents; on the contrary, the non-carcinogenic risks are likely to occur when THQ or HI is >1.0, and the human health risk increases with the increase of THQ or HI.

2.4.5. Statistical Analysis

Excel 2007 (Microsoft Office, Las Vegas, NV, USA), SPSS 23.0 (IBM Corporation, New York, NY, USA), and Origin 8.5, 2020b (OriginLab Corporation, Northampton, MA, USA) were used to conduct experimental data processing and statistical analysis. The Kolmogorov–Smirnov (K-S) test was used to determine whether the HMs’ contents were normally distributed in all areas [42]. A significant variance analysis was employed by using a one-way ANOVA test (LSD at p < 0.05 and p < 0.01). Pearson correlation analysis was used to clarify the relationships between total heavy metal concentrations and vegetable species in the different sites. Principal component analysis (PCA) and cluster analysis (CA) were conducted to identify possible sources in the study areas. ArcGIS 10.6 software (Environmental Systems Research Institute, Redlands, CA, USA) was used to evaluate spatial distribution levels with Igeo values.

3. Results

3.1. Heavy Metal Concentrations in Soils

The concentrations of HMs in soils are shown in Table 2. The pH of the soil in the four study areas showed a distinct pattern, with the mean values increasing in order as follow: D zone (pH = 5.05) < A zone (pH = 6.71) < B zone (pH = 7.03) < C zone (pH = 7.72). Most of the results of K-S test confirmed that the concentrations of HMs in the four zones passed the K-S test for normality (p > 0.05). Except Cr, all the soil HMs contents were higher than the background values of soil in Hubei Province. Compared with the soil environmental quality standard of agricultural land (GB 15618-2018), agricultural soils in the exposed zones were identified as contaminated areas by Cu, Cd, As and Pb, so the total excessive rate of them followed a decreasing order of Cd (80.00%) > Cu (61.11%) > As (45.56%) > Pb (32.22%) > Cr (0.00%).
In the A zone, the mean contents of Cu and Cd were significantly higher than the national standard values in all sampling sites, and the excessive rate followed the order of Cu (100%) = Cd (100%) > As (71.4%) > Pb (61.9%) > Cr (0.0%). In the B zone, the excessive rate values of HMs decreased in the following order of Cu (100%) > Cd (88.46%) > As (53.85%) > Pb (42.31%) > Cr (0.0%). In the C zone, the excessive rate of HMs increased in the following order of Cu (30.0%) < As (60.0%) < Cd (85.0%). However, in the D zone, only the mean contents of Cu, Pb, and Cd exceeded the soil standard with the excessive rates of 8.7%, 21.74%, and 47.83%, respectively. Compared to the four zones, almost all the HMs concentrations of soil near the non-ferrous metals smelter, abandoned copper mine and limestone quarry were significantly higher (p < 0.05) than those of iron mine, especially Cu, Cd, and As showed more extreme levels of pollution. Furthermore, the coefficients of variation (CV%) for Cu, Cd, Pb, and As were over 50% in the four study areas, indicating that a larger influence of anthropogenic activities contributed to the heavy metals in the soils.

3.2. Contaminated Evaluation of HMs in Soils

3.2.1. Nemerow’s Synthetical Pollution Index (Pn) and Potential Ecological Risk (RI)

In order to obtain an accurate contaminated status of different zones, the values of the single-factor index and Nemerow’s synthetical pollution index were calculated (Table S1). The most soil samples were polluted to varying levels ( p i   > 1) except for Cr. The results indicate that the A and B zones were highly polluted by Cd ( P i   = 12.98) and Cu ( P i   = 12.25). As illustrated in Figure 2a, the P n values ranged from 0.56 to 18.44, exposing multiple ranges from safety to serious pollution. Thus, the pollution index ( P n ) of four regions decreased in the following order: A zone > B zone > C zone > D zone.
As shown in Figure 2b, the higher potential ecological risk was found in A and B zones in comparison with C and D zones. More specifically, in the A zone, the single ecological risk index ( E r i ) of Cd was higher than 320 ( E r i = 389.4), indicating a higher contribution of Cd to potential ecological risk (Figure S1a). In the B zone, except Cd, the E r i value of Cu was larger than 40 ( E r i   = 63.0), posing a moderate potential ecological risk though RI being lower risk level than non-ferrous metals smelter (Figure S1b). In the C zone, only the single ecological risk index of Cd was more than 40 ( E r i   = 55.2), suggesting a moderate potential ecological risk (Figure S1c), while D zone corresponded little sensitivity of biological community to HMs in soils because all the E r i  values were less than 40 except Cd (Figure S1d). Based on the results of E r i and RI in the study areas, Cd and Cu in four categories of vegetable garden soil exhibited higher ecological risk.

3.2.2. Geo-Accumulation Index (Igeo)

The spatial distributions of Igeo values of HMs in soils are shown in Figures S2–S5. In this study, we found that the Igeo values were greater than zero to some degree in the three different categories of vegetable fields excluding the D zone. Among the five hazardous metals, the pollution levels of Igeo-Cd showed larger spatial variation in Figure 3, suggesting that different anthropogenic activities can contribute to its distributions. In contrast, the negative Igeo value of Cr is observed in Figure S3, the concentration of Cr was obviously below the background values, which indicated the study area free from Cr contamination.

3.3. Heavy Metals in Vegetables and Crops

The concentrations of HMs in edible parts of vegetables and crops in the study areas are presented in Table S2. The mean contents of Cu, Cr, Pb, Cd, and As from vegetables and crops were 4.05, 0.22, 0.41, 0.15, and 0.26 mg/kg (FW) in the exposed areas (A–C zones), respectively, which was 2–8 folds higher than those in D zones (p < 0.05). Additionally, several larger coefficients of variation (CV% ≥ 100%) of HMs were displayed in the study areas, suggesting that a wide variation of values could be related to the difference of vegetable species.
In present study, the maximum allowable concentration (MAC) of food standard (GB2762-2017) is presented in Table S3. Compared to 15 species of vegetables and crops, the highest Cu contents were found in brassica campestris (Figure 4a), which almost surpassed the MAC value, while the higher Cr content was determined in peanut (Figure 4b). Moreover, the species exceeding rates of Cu and Cr were 6.7%. According to the MAC, Pb and Cd contents in most leaf vegetables were found above the criterion line as shown in Figure 4c,d, whereas relatively low concentration was found in tomato, bell pepper, white radish, and asparagus bean. In contrast, several crop types including peanut, ipomoea, maize and soya bean showed the higher accumulation ability for Pb and Cd. So, the percentage of vegetables and crops that exceed the MAC value is 66.7% for Pb and 80.0% for Cd. For As, only 26.7% species exceed the MAC, and the contents of As in the vegetables are listed in Figure 4e, the higher value can be observed in leaf vegetables, especially in edible amaranth.
Overall, in tomato, bell pepper, white radish, and asparagus bean, five hazardous metals contents were rarely observed, indicating low accumulation characteristics to multiple hazardous metals. Furthermore, the observed HMs concentrations of all collected samples in the contaminated areas (A–C zones) were significantly higher than those in D zone.

3.4. Source Analysis of HMs in Different Contaminated Soils

Pearson correlation analysis was applied to clarify the relationship of soil HMs in the contaminated zones. As shown in Figure 5, the total amount of several HMs in soils showed a significantly positive correlation between pairs of elements in the four study areas. The correlation coefficients of Cu-Pb, Cu-Cd, Cu-As, Pb-Cd, Pb-As, and Cd-As in the A zone were 0.91, 0.76, 0.87, 0.75, 0.89, and 0.77 (p < 0.01), respectively. Furthermore, the correlation coefficients of Cu-Cr and Pb-Cd in the B zone were 0.64 and 0.68 (p < 0.01), Cr-Pb and Pb-As were 0.46 and 0.47 (p < 0.05), respectively. The positive correlation coefficients of Cu-Pb (0.72) and Cd-As (0.45) were also observed in the C zone, and Cu-Cr (0.91) was significantly correlated in the D zone, which indicated that these elements were derived from the same sources. Nevertheless, the negative correlation coefficients of Cr-Cd, Cr-As, and Pb-Cd in C and D zones were calculated to be −0.56, −0.49, and −0.48, which suggested that these elements could originate from different sources. As shown in Figure S6, the correlation coefficients between soils and foodstuffs were not found at significant range (p < 0.05), indicating that they were barely or negatively correlated with each other. The possible explanation is that the enrichment of HMs in foods is affected not only by soil HMs, but also by vegetable and crop species or other complex factors [22,29].
Principal component analysis (PCA) and cluster analysis (CA) were used to further identify the sources among the five heavy metals in the study areas. The results of the PCA are shown in Figure S7 and Table S4. Two principal components were extracted in the four study zones because the eigenvalues were greater than 1. In the A zone, PC1 explained 69.70% of the variance and positively related to Cu, Pb, Cd, and As, while PC2 accounted for 22.13% of the variance, which was dominated by Cr. In the B zone, PC1 explained 49.13% of the variance and depended partially on Pb, Cd, and As; meanwhile, PC2, which was dominated by Cu and Cr, accounted for 23.14% of the variance. In the C zone, PC1 occupied 45.50% of the variance and loaded positively on Cd and As, additionally, PC2 accounted for 29.73% of the variance and was dominated by Cu and Pb. In the D zone, PC1 occupied 43.67% of the variance and depended positively on Cu and Cr, while PC2 accounted for 27.52% of the variance and was dominated by Pb. Based on the results of the previous spatial distribution, the concentrations of several HMs (Cu, Pb, Cd, and As) were higher than the soil background values of Hubei province. Thus, the PCA results revealed that these HMs were derived from anthropogenic activities, Cr could originate from natural process.
The dendrogram results obtained from CA for HMs are shown in Figure S8. Two main clusters were distinguished as following: (1) Cu-Pb-Cd-As, (2) Cr in the A and C zone; (1) Pb-Cd-As, (2) Cu-Cr in the B and D zones, respectively. It revealed similar results as Pearson correlation and principal component analysis.

3.5. Health Risk Assessment of Dietary Intake

The HI and THQ were applied to assess the non-carcinogenic risk of food intake which was considered as the main exposed risk pathway. In Table S5, the HI values for all zones were over 1, indicating a higher non-carcinogenic risk of dietary intake for local children. These findings are consistent with the results of earlier research [31,40]. In addition, the HI values of four investigated zones followed a decreasing order: B zone > A zone > C zone > D zone. As shown in Figure 6, the THQ values of solanaceous fruits were lower than those of the other five types. In contrast, the highest HI values were observed in leaf vegetables posing the highest potential health risk for local children.
The health risk contribution rates of different elements in the four zones are shown in Figure S9. Among the five hazardous metals, As was the primary contributor for health risk in all the four zones, and the result of contribution rate were 67.71% to A zone, 72.31% to B zone, 67.87% to C zone, and 43.48% to D zone.

4. Discussion

4.1. Pollution and Source Identification in Soils

As numerous studies reported, heavy metal pollution of soil is a severe environmental concern, especially in areas where agricultural soils are close to anthropogenic pollution, such as mining and smelting zones [1,8], industrial suburban areas [43], e-waste dismantling areas and traffic roadsides [16,44,45]. Our results indicated that the contents of Cu and Cd in all sites near the Tonglushan mine and copper smelter were lower than those in previous studies [40]. The level of soil pollution is affected by source distance, Li et al. [46] and Sun et al. [47] also presented that sampling sites closer to the mine were more seriously polluted by Cu, Cd, Pb, and As, and as expected, the vegetables and crops were also contaminated by these toxic metals to a greater extent, which threatened the health of residents. In addition to the variation of sampling sites, the possible reason is that these toxic metals, especially Cu and Cd, are easily absorbed by crop uptake or phytoextraction from soils. Furthermore, the contents of Cd and Cu in current study were significantly higher than those of other areas in Daye City (Cd: 1.41 mg/kg and Cu: 105 mg/kg) [35]. The pollution of HMs was serious in the agricultural soils, so the study zones chosen could be more vulnerable areas impacted by mining and smelting. In comparison with other countries, the concentrations of Cu, Pb, and Cd were much higher than those in Bangladesh [4]. However, in Germany and Northern Tunisia, the soil Pb concentration was 10-fold higher than that in this study [24,26].
The four study zones are representative areas polluted by HMs from different anthropogenic activities. In the A zone, an elevated concentration of Cd was observed in soil-vegetable system. Many studies showed that non-ferrous smelter could lead to a large amount of Cd emitting into the soils through rainfall and deposition [48]. According to this study, Cu, Pb, Cd, and As may originate from smelting activities. Bi et al. [16] also reported that the Cd, Pb, Cu, and As contents in soil and road dust exceeded the referenced values near industrial facilities in Shanghai, suggesting that one of the major sources of HMs pollution were industrial sources, including metal and metal compound application. It is evident in the B zone that mining operations have polluted vast fields of agricultural land. Pb, Cd, and As are known to be the associated elements of Cu, and primarily derived through waste water discharge, spoil disposal, and dust fall from copper mining processes. In the C zone, enrichment of As can be related with agricultural activities, as the superfluous use of pesticides, herbicide, phosphate fertilizers, and livestock manure are also important sources for As in soil [49,50]. In addition, vehicle exhaust emissions and road dust can be considered as other major sources for Pb and Cd in soils [51,52]. Several cement plants and lime quarries were installed along the roadside so the contribution of transportation to Pb cannot be ignored. Similar findings were found in the D zone, the major sources of Pb could be considered as atmospheric deposition and automobile exhaust. According to PCA results, Cr may originate from natural sources related to parent materials. HMs, one of the massive rock elements, can be released into the environment by rock weathering, rain erosion and volcanic eruption. Cr, an important source of mineral element, can reach into soil-foodstuffs-human chain by radicular system absorption.
Based on the above mentioned argument, the distribution patterns of HMs in mining and smelting areas were similar to those in industrial areas, which were also regarded as the most sensitive regions for HMs contamination. Due to multiple sources of heavy metals emission in different regions, it is difficult to identify the source of one heavy metal pollution. Yet we might at least find a great link between emissions of heavy metals and anthropogenic sources.

4.2. Comparisons of HMs Accumulation in Different Vegetable and Crop Species

Previous studies have reported that elevated concentration of HMs was found near the multiple-metal mining and smelter, e-waste dismantling area, and industrial sites [16,23,45,53]. In addition, Rahmdel et al. [17] and Mandal et al. [20] reported that higher concentrations of Cu and Pb in leaf vegetables were observed in Iran and India than those in the present study, in which samples were collected near the rural districts of industrial town. Comparatively, in China, the average concentrations of HMs in vegetables were mostly similar to those in the peri-urban areas, but most of them were significantly higher than those from Pb/Zn smelters in Gejiu of Yunnan province [23], and seven areas in Zhuzhou of Hunan Province [22]. Therefore, we inferred that this was possibly linked to the industrial activities of local copper mining and smelter operations in Daye City. For over 3000 years, the Tonglushan Mine has become an inactive copper mine. Moreover, the Daye Nonferrous Metal Smelting Plant has also been releasing numerous toxic metals as by-products into the atmosphere and soil for almost 50 years. Along these lines, Cai et al. [40] measured the amount of contamination of food crops grown near the Tonglushan Mine of Daye City, which found that the contents of Cd and As in rice and most mine-affected vegetables were above the maximum permissible limits. Similar results also were obtained in Daye copper smelter. Du et al. [35] also indicated that in the vicinity of the active Pb-Zn mine and Cu smelter, a higher environmental risk can occur to rice and vegetable due to Cd exposure.
Referring to the sources of heavy metals in vegetables and crops, previous studies have revealed relationship between physiochemical properties of soil and the contents of HMs in foods. They have established certain migration and transport models for predicting bioaccumulation factors of HMs in soil-vegetable or soil-crop system [14,54]. In this study, five heavy metals, especially Pb, Cd, and As, were highly enriched in leaf vegetables, while barely accumulated in fruit vegetables collected from the heavy contaminated sites. Similar to our findings, Zhou et al. [55] and Hu et al. [56] also reported that leaf vegetables can uptake more heavy metals than rootstalk and fruit vegetables. Generally, the contents of Cd in leaf vegetables were nearly 3 times than those in non-leaf vegetables. By Comparing the contamination levels of vegetables in the open field with those in greenhouse, Li et al. [57] found that the accumulation of Cd in most vegetables was mainly from soil rather than atmospheric deposition. This result further confirmed that Cd is typically a labile trace element in the soil, and many vegetables, especially leaf vegetables, are more capable of absorbing it. The higher translocation rates of Cd were also found even when vegetables were grown in low-Cd soils [13,58]. Moreover, the accumulation of Cd in vegetables is closely related to the amount of soil extractable Cd depending primarily on the properties of the soil, such as pH, organic matter, cation exchange capacity and clay content [13,59,60]. However, this is difficult for Pb to migrate through the soil-to-crops due to its low solubility and bioavailability [30,61]. In order to explore the influence of airborne heavy metals on vegetable species, Harrison et al. [62] confirmed that Pb in the atmosphere contributed as much as 85% to spinach pollution, while Cd only contributed 23%. In addition, Bi et al. [16] further demonstrated that atmospheric Pb via foliar uptake is the primary source for the accumulation of Pb in leaf vegetables by using 208Pb/206Pb and 207Pb/206Pb ratios of isotopic composition. In terms of As accumulation, the higher concentration of As was found in leaf vegetables than fruits vegetables [34]. In general, Cr and Cu were readily absorbed by roots and their distribution was ordered as following root > stem > leaf [63]. Similarly, the higher concentrations of Cr were observed in rice and wheat [12]. However, our results indicated that the elevated Cr contents in ipomoea and peanut were accumulated mainly from soils, implying high potential health risk via the consumption of rhizome geophytes.
Vegetables and crops grown in contaminated farmlands could have higher accumulation capabilities for HMs than those cultivated in the safe areas [47]. However, we found that tomato, bell pepper, white radish, and asparagus bean have lower heavy metal contents than leaf vegetables, suggesting that these species showed low accumulation characteristics for hazardous heavy metals. In general, the accumulation of HMs in vegetables could be affected by types and species as well as spatial distribution of HMs in soils [64]. Low accumulation abilities may be governed by the expression of specific genes [33]. Previous studies have provided genetic basis for screening and breeding low accumulating species such as rice, wheat, Chinese cabbage, and tomato [31,33,65]. On the basis of the results of this research, it is believed that foodstuffs with low accumulation are found in contaminated soils, suggesting that these genotypes can be treated as reliable native cultivars.

4.3. Health Risk Assessment for Local Children

Contaminated soil particles and foods are the most important exposures for humans, which can lead to potential health risks for residents. Based on non-carcinogenic risk assessment proposed by USEPA [66,67,68], the average daily dose (ADD) was dependent on 3 factors including ingestion, inhalation, and dermal contact. Generally, children are more sensitive to hazardous elements via non-dietary ingestion than adults because of physiological activities, such as rubbing eyes and sucking fingers [37,69,70]. Furthermore, food ingestion was the leading contributor to children’s ADD of HMs in comparison with soil ingestion and dermal [71]. Similar findings have been demonstrated by Huang et al. [2] and Cai et al. [40]. In this study, local children living around mining and smelter areas could suffer from adverse health impacts because the HIs of leaf vegetables were higher than 1. Interestingly, the HI values of solanaceous fruits and legumes were less than 1, suggesting that the presence of low accumulation foods was healthy for local children.
Previous studies indicated the higher accumulation ability for Cu, Pb, and Cd in leaf, rhizome vegetables, and rice grain [2,22,64]. In this study, As had a higher THQ value than other hazardous heavy metals. A similar report also confirmed that As exposure through vegetable ingestion has significant adverse effects on residents [23]. Thus, much attention should be paid to As although there is a lack of detailed criteria of chemical form. It is a well-established fact that the availability and toxicity of HMs depend on speciation and chemical states of HMs. Such as, As (III), binding to certain enzymes and proteins from phytochelatins, is recognized as a more poisonous form than As (V), and it is readily absorbed via plant roots from the soil [72]. Dhanker et al. [72] documented that As (V) can be reduced into As (III) as a result of arsenate reductase inside most plant tissues. In general, detoxification of As can be achieved via methylation reaction with low toxic intermediate organic arsenic chemicals including monomethylarsonic acid (MMA), dimethylarsinic acid (DMA), and trimethylarsine oxide (TMAO) [73]. However, Zhou et al. [74] found that monomethylarsonic acid (MMAIII) and dimethylarsinic acid (DMAIII) in the human body may be more or less cytotoxic to As (III), which were recognized as human hepatocytes impairer.
Cr also exists in two common types of oxidation states including Cr (III) and Cr (VI). In comparison, Cr (VI) originates from some anthropogenic compounds in the processing of industry and smelter, which is the most toxic form due to higher solubility, mobility and oxidizability, thus Cr (VI) can enter into human body via food chain and eventually damage tissues. According to previous studies, the high Cr levels had been found in several foods, posing potential human health risk [75,76,77]. Due to the elevated toxicity in different chemical states, a strict permissible risk index for foods should be developed for metal states like Cr (VI) and As (III). To sum up, it is necessary to clear up the relationship between element chemical form and human health risk.

5. Conclusions

In this study, four typical vegetable fields around the non-ferrous metals smelter, abandoned copper mine, limestone quarry, and iron mine in Daye City were selected as study zones. The results showed that the pollution of non-ferrous metals smelter and abandoned copper mine was serious, and heavy metals pollution was mainly attributed to mining and metal smelting activities. The main polluted metals were followed a decreasing order of Cd (80.00%) > Cu (61.11%) > As (45.56%) > Pb (32.22%) > Cr (0.00%). According to the results of spatial distribution and PCA, these five elements could be classified into two groups: (1) PC1 (Cu, Pb, Cd and As) which may originate from anthropogenic activities; (2) PC2 (Cr) which was mostly associated with the parent material of soil. Compared with vegetable and crop species, the contents of HMs in tomato, bell pepper, white radish, and asparagus bean were relatively low or closed zero, possibly due to their low accumulation characteristics. However, the elevated contents of Cd, Pb, and As were observed in peanut, ipomoea, maize, soya bean, and most leaf vegetables.
In addition, the results of health risk assessment revealed that the HI values of four zones through consumption of vegetables and crops were >1, revealing that the local children near mining and smelter regions suffer a higher non-carcinogenic risk. In contrast, the consumption of solanaceous fruits could cause a low potential health risk, which suggested that low accumulation species should be selected and applied on contaminated agricultural soils. Thus, adjustment of planting structure or conversion of land use should be adopted to protect human health and ensure food safety.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/18/5/2617/s1, Table S1: Results of single-factor index and Nemerow’s synthetical pollution index. Table S2: Concentration (mg/kg, Fresh weight) of heavy metals in the foodstuffs of the four categories of vegetable gardens. Table S3: The assessment criteria in edible parts of different vegetables (mg/kg). Table S4: Principal component analysis for heavy metals of soils. Table S5: Indices of non-carcinogenic risk for children. Figure S1: Potential ecological risk index (RI) of single element in the four zones: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine. Figure S2: Spatial distribution pollution levels of Igeo-Cu in the four zones: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine. Figure S3: Spatial distribution pollution levels of Igeo-Cr in the four zones: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine. Figure S4: Spatial distribution pollution levels of Igeo-Pb in the four zones: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine. Figure S5: Spatial distribution pollution levels of Igeo-As in the four zones: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine. Figure S6: Correlation coefficients of heavy metals between soils and foodstuffs. Figure S7: Principal component analysis for heavy metals of soils in the four zones. Figure S8: Cluster analysis (CA) results for five heavy metals of soils (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine. Figure S9: The health risk contribution rates of five different elements in the four zones.

Author Contributions

Z.W.: Writing-original draft, Methodology, Software; J.B.: Conceptualization, Methodology; T.W.: Investigation; H.T.M.: Formal analysis; W.K.: Writing-review & editing; J.Z.: Funding acquisition; C.Z.: Data curation; W.X.: Supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (NO. 42077312), the Open Foundation of Hubei Key Laboratory of Mine Environment Pollution Control & Remediation (2019105), the Foundation of State Key Laboratory of Loess and Quaternary Geology (SKLLQG1326), Key Scientific Research Project of Hubei Polytechnic University (19XJK01Z), and Coopertion Project of Hubei Environmental Remediation & Goverance Techological Research Co., Ltd. And Hubei Polytechnic University (KY2018-180).

Institutional Review Board Statement

Ethical review and approval were waived for this study, due to this study did not involve biological humans or animals experiment and data.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

We gratefully thank to the editors and the reviewers for valuable comments that have improved the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhuang, P.; McBride, M.B.; Xia, H. Health Risk from Heavy Metals Via Consumption of Food Crops in the Vicinity of Dabaoshan Mine, South China. Sci. Total Environ. 2009, 407, 1551–1561. [Google Scholar] [CrossRef] [PubMed]
  2. Huang, Y.; Chen, Q.Q.; Deng, M.H.; Japenga, J.; Li, T.Q.; Yang, X.E.; He, Z. Heavy Metal Pollution and Health Risk Assessment of Agricultural Soils in a Typical Peri-Urban Area in Southeast China. J. Environ. Manag. 2018, 207, 159–168. [Google Scholar] [CrossRef] [PubMed]
  3. Cai, L.M.; Xu, Z.C.; Qi, J.Y.; Feng, Z.Z.; Xiang, T.S. Assessment of Exposure to Heavy Metals and Health Risks among Residents near Tonglushan Mine in Hubei, China. Chemosphere 2015, 127, 127–135. [Google Scholar] [CrossRef]
  4. Ahmed, M.; Matsumoto, M.; Kurosawa, K. Heavy Metal Contamination of Irrigation Water, Soil, and Vegetables in a Multi-Industry District of Bangladesh. Int. J. Environ. Res. 2018, 12, 531–542. [Google Scholar] [CrossRef]
  5. Vatanpour, N.; Malvandi, A.M.; Talouki, H.H.; Gattinoni, P.; Scesi, L. Impact of Rapid Urbanization on the Surface Water’s Quality: A Long-term Environmental and Physicochemical Investigation of Tajan River, Iran (2007–2017). Environ. Sci. Pollut. Res. Int. 2020, 27, 8439–8450. [Google Scholar] [CrossRef] [PubMed]
  6. MEP of China. National Soil Pollution Survey Bulletin. 2014. Available online: http://www.zhb.gov.cn/gkml/hbb/qt/201404/t20140417-270670.htm (accessed on 10 February 2020).
  7. Zhao, F.J.; Ma, Y.; Zhu, Y.G. Soil Contamination in China: Current Status and Mitigation Strategies. Environ. Sci. Technol. 2015, 49, 750–759. [Google Scholar] [CrossRef] [PubMed]
  8. Li, Z.Y.; Ma, Z.W.; van der Kuijp, T.J.; Yuan, Z.W.; Huang, L. A Review of Soil Heavy Metal Pollution from Mines in China: Pollution and Health Risk Assessment. Sci. Total Environ. 2014, 468, 843–853. [Google Scholar] [CrossRef]
  9. Anju, M.; Banerjee, D.K. Multivariate Statistical Analysis of Heavy Metals in Soils of a Pbezn Mining Area, India. Environ. Monit. Assess. 2012, 184, 4191–4206. [Google Scholar] [CrossRef]
  10. Gu, Y.G.; Li, Q.S.; Fang, J.H.; He, B.Y.; Fu, H.B.; Tong, Z.J. Identification of Heavy Metal Sources in the Reclaimed Farmland Soils of the Pearl River Estuary in China Using a Multivariate Geostatistical Approach. Ecotox. Environ. Saf. 2014, 105, 7–12. [Google Scholar] [CrossRef]
  11. Facchinelli, A.; Sacchi, E.; Mallen, L. Multivariate Statistical and Gis-Based Approach to Identify Heavy Metal Sources in Soils. Environ. Pollut. 2001, 114, 313–324. [Google Scholar] [CrossRef]
  12. Chen, L.; Zhou, S.L.; Shi, Y.X.; Wang, C.H.; Li, B.J.; Li, Y.; Wu, S.H. Heavy metals in Food Crops, Soil, Andwater in the Lihe River Watershed of the Taihu Region and Their Potential Health Risks When Ingested. Sci. Total Environ. 2018, 615, 141–149. [Google Scholar] [CrossRef]
  13. Xiao, W.; Ye, X.Z.; Zhang, Q.; Chen, D.; Hu, J.; Gao, N. Evaluation of Cadmium Transfer from Soil to Leafy Vegetables: Influencing Factors, Transfer Models, and Indication of Soil Threshold Contents. Ecotox. Environ. Saf. 2018, 164, 355–362. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, J.R.; Li, H.Z.; Zhou, Y.Z.; Dou, L.; Cai, L.M.; Mo, L.P.; You, J. Bioavailability and Soil-to-Crop Transfer of Heavy Metals in Farmland Soils: A Case Study in the Pearl River Delta, South China. Environ. Pollut. 2018, 235, 710–719. [Google Scholar] [CrossRef]
  15. Kumar, S.; Prasad, S.; Yadav, K.K.; Shrivastava, M.; Gupta, N.; Nagar, S.; Bach, Q.V.; Kamyab, H.; Khan, S.A.; Yadav, S.; et al. Hazardous Heavy Metals Contamination of Vegetables and Food Chain: Role of Sustainable Remediation Approaches—A Review. Environ. Res. 2019, 179, 108792. [Google Scholar] [CrossRef] [PubMed]
  16. Bi, C.J.; Zhou, Y.; Chen, Z.L.; Jia, J.P.; Bao, X.Y. 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, 1349–1357. [Google Scholar] [CrossRef]
  17. Rahmdel, S.; Rezaei, M.; Ekhlasi, J.; Zarei, S.H.; Akhlaghi, M.; Abdollahzadeh, S.M.; Sefidkar, R.; Mazloomi, S.M. Heavy Metals (Pb, Cd, Cu, Zn, Ni, Co) in Leafy Vegetables Collected from Production Sites: Their Potential Health Risk to the General Population in Shiraz, Iran. Environ. Monit. Assess. 2018, 190, 650. [Google Scholar] [CrossRef]
  18. Huang, L.K.; Wang, Q.; Zhou, Q.Y.; Ma, L.Y.; Wu, Y.J.; Liu, Q.Z.; Wang, S.; Feng, Y. Cadmium Uptake from Soil and Transport by Leafy Vegetables: A Meta-Analysis. Environ. Pollut. 2020, 264, 114677. [Google Scholar] [CrossRef] [PubMed]
  19. Jia, Z.M.; Li, S.Y.; Wang, L. Assessment of Soil Heavy Metals for Eco-environment and Human Health in a Rapidly Urbanization Area of the Upper Yangtze Basin. Sci. Rep. 2018, 8, 3256. [Google Scholar] [CrossRef] [Green Version]
  20. Mandal, R.; Kaur, S. Impact of Environmental Pollution on Trace Elements in Vegetables and Associated Potential Risk to Human Health in Industrial Town Mandi-Gobindgarh (India). Chemosphere 2019, 219, 574–587. [Google Scholar] [CrossRef]
  21. Chen, M.; Li, F.G.; Hu, L.W.; Yang, T.; Yang, Q.; Tao, M.X.; Deng, Y.W. Heavy Metal Pollution in Topsoil and Vegetables in the Typical Mining Area near Gannan, Jiangxi Province, China. Environ. Eng. Sci. 2019, 10, 1307–1314. [Google Scholar] [CrossRef]
  22. Li, X.Y.; Li, Z.G.; Lin, C.J.; Bi, X.Y.; Liu, J.L.; Feng, X.; Zhang, H.; Chen, J.; Wu, T.T. Health Risks of Heavy Metal Exposure through Vegetable Consumption near a Large-Scale Pb/Zn Smelter in Central China. Ecotox. Environ. Saf. 2018, 161, 99–110. [Google Scholar] [CrossRef] [PubMed]
  23. Guo, G.H.; Zhang, D.G.; Wang, Y.T. Probabilistic Human Health Risk Assessment of Heavy Metal Intake Via Vegetable Consumption around Pb/Zn Smelters in Southwest China. Int. J. Environ. Res. Public Health 2019, 16, 3267. [Google Scholar] [CrossRef] [Green Version]
  24. Antoniadis, V.; Shaheen, S.M.; Boersch, J.; Frohne, T.; Laing, G.D.; Rinklebe, J. Bioavailability and Risk Assessment of Potentially Toxic Elements in Garden Edible Vegetables and Soils around a Highly Contaminated Former Mining Area in Germany. J. Environ. Manag. 2017, 186, 192–200. [Google Scholar] [CrossRef] [PubMed]
  25. Douay, F.; Pelfrene, A.; Planque, J.; Fourrier, H.; Richard, A.; Roussel, H.; Girondelot, B. Assessment of Potential Health Risk for Inhabitants Living near a Former Lead Smelter. Part 1: Metal Concentrations in Soils, Agricultural Crops, and Homegrown Vegetables. Environ. Monit. Assess. 2013, 185, 3665–3680. [Google Scholar] [CrossRef]
  26. Boussen, S.; Soubrand, M.; Bril, H. Transfer of Lead, Zinc and Cadmium from Mine Tailings to Wheat (Triticum Aestivum) in Carbonated Mediterranean (Northern Tunisia) Soils. Geoderma 2013, 192, 227–236. [Google Scholar] [CrossRef]
  27. Ji, K.; Kim, J.; Lee, M.; Park, S.; Kwon, H.J.; Cheong, H.K.; Jang, J.Y.; Kim, D.S.; Yu, S.; Kim, Y.W.; et al. Assessment of Exposure to Heavy Metals and Health Risks among Residents near Abandoned Metal Mines in Goseong, Korea. Environ. Pollut. 2013, 178, 322–328. [Google Scholar] [CrossRef] [PubMed]
  28. Vatanpour, N.; Feizy, J.; Talouki, H.H.; Es’haghi, Z.; Scesi, L.; Malvandi, A.M. The High Levels of Heavy Metal Accumulation in Cultivated Rice from the Tajan River Basin: Health and Ecological Risk Assessment. Chemosphere 2020, 245, 125639. [Google Scholar] [CrossRef]
  29. Chen, H.Y.; Yuan, X.Y.; Li, T.Y.; Hu, S.; Ji, J.F.; Wang, C. Characteristics of Heavy Metal Transfer and Their Influencing Factors in Different Soil-Crop Systems of the Industrialization Region, China. Ecotox. Environ. Saf. 2016, 126, 193–201. [Google Scholar] [CrossRef]
  30. Bi, X.Y.; Feng, X.B.; Yang, Y.G.; Li, X.D.; Shin, G.P.Y.; Li, F.L.; Qiu, G.L.; Li, G.H.; Liu, T.Z.; Fu, Z.Y. Allocation and Source Attribution of Lead and Cadmium in Maize (Zea Mays L.) Impacted by Smelting Emissions. Environ. Pollut. 2009, 157, 834–839. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, W.T.; Zhou, Q.X.; Sun, Y.B.; Liu, R. Identification of Chinese Cabbage Genotypes with Low Cadmium Accumulation for Food Safety. Environ. Pollut. 2009, 157, 1961–1967. [Google Scholar] [CrossRef]
  32. Liu, W.T.; Zhou, Q.X.; Zhang, Y.L.; Wei, S.H. Lead Accumulation in Different Chinese Cabbage Cultivars and Screening for Pollution-Safe Cultivars. J. Environ. Manag. 2010, 91, 781–788. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, Z.M.; Tan, X.Q.; Mei, X.Q.; Li, Q.S.; Zhou, C.; Wang, L.L.; Ye, H.J.; Yang, P. Low-Cd Tomato Cultivars (Solanum Lycopersicum L.) Screened in Non-Saline Soils Also Accumulated Low Cd, Zn, and Cu in Heavy Metal-Polluted Saline Soils. Environ. Sci. Pollut. Res. Int. 2018, 25, 27439–27450. [Google Scholar] [CrossRef] [PubMed]
  34. Yang, J.; Lv, F.H.; Zhou, J.C.; Song, Y.W.; Li, F. Health Risk Assessment of Vegetables Grown on the Contaminated Soils in Daye City of Hubei Province, China. Sustainability 2017, 9, 2141. [Google Scholar] [CrossRef] [Green Version]
  35. Du, P.; Xie, Y.F.; Wang, S.J.; Zhao, H.H.; Zhang, Z.; Wu, B.; Li, F.S. Potential Sources of and Ecological Risks from Heavy Metals in Agricultural Soils, Daye City, China. Environ. Sci. Pollut. Res. Int. 2015, 22, 3498–3507. [Google Scholar] [CrossRef]
  36. Wang, M.L.; Liu, R.H.; Lu, X.Y.; Zhu, Z.Y.; Wang, H.L.; Jiang, L.; Wu, Z.H. Heavy Metal Contamination and Ecological Risk Assessment of Swine Manure Irrigated Vegetable Soils in Jiangxi Province, China. Bull. Environ. Contam. Toxicol. 2018, 100, 634–640. [Google Scholar] [CrossRef]
  37. Sawut, R.; Kasim, N.; Maihemuti, B.; Hu, L.; Abliz, A.; Abdujappar, A.; Kurban, M. Pollution Characteristics and Health Risk Assessment of Heavy Metals in the Vegetable Bases of Northwest China. Sci. Total Environ. 2018, 642, 864–878. [Google Scholar] [CrossRef]
  38. Hakanson, L. An Ecological Risk Index for Aquatic Pollution Control. A Sedimentological Approach. Water Res. 1980, 14, 975–1001. [Google Scholar] [CrossRef]
  39. Müller, G. Heavy Metals in Sediment of the Rhine-Changes since 1971. Umsch. Wiss. Tech. 1979, 79, 778–783. [Google Scholar]
  40. Cai, L.M.; Wang, Q.S.; Luo, J.; Chen, L.G.; Zhu, R.L.; Wang, S.; Tang, C.H. Heavy Metal Contamination and Health Risk Assessment for Children near a Large Cu-Smelter in Central China. Sci. Total Environ. 2019, 650, 725–733. [Google Scholar] [CrossRef]
  41. USEPA. Human Health Risk Assessment. In Regional Screening Level-Summary Table; United States Environmental Protection Agency: Washington, DC, USA, 2018. [Google Scholar]
  42. Tian, K.; Biao, H.; Xing, Z.; Hu, W.Y. Geochemical Baseline Establishment and Ecological Risk Evaluation of Heavy Metals in Greenhouse Soils from Dongtai, China. Ecol. Indic. 2017, 72, 510–520. [Google Scholar] [CrossRef]
  43. Yang, Q.Q.; Li, Z.Y.; Lu, X.N.; Duan, Q.N.; Huang, L.; Bi, J. 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. [Google Scholar] [CrossRef]
  44. Luo, C.; Liu, C.; Wang, Y.; Liu, X.; Li, F. Heavy Metal Contamination in Soils and Vegetables near an E-Waste Processing Site, South China. J. Hazard. Mater. 2011, 186, 481–490. [Google Scholar] [CrossRef]
  45. Zhang, T.; Ruan, J.J.; Zhang, B.; Lu, S.Y.; Gao, C.Z.; Huang, L.F.; Bai, X.Y.; Xie, L.; Gui, M.W.; Qiu, R.L. Heavy Metals in Human Urine, Foods and Drinking Water from an E-Waste Dismantling Area: Identification of Exposure Sources and Metal-Induced Health Risk. Ecotox. Environ. Saf. 2019, 169, 707–713. [Google Scholar] [CrossRef]
  46. Li, N.; Kang, Y.; Pan, W.J.; Zeng, L.X.; Zhang, Q.Y.; Luo, J.W. Concentration and Transportation of Heavy Metals in Vegetables and Risk Assessment of Human Exposure to Bioaccessible Heavy Metals in Soil near a Waste-incinerator Site, South China. Sci. Total Environ. 2015, 521–522, 144–151. [Google Scholar] [CrossRef]
  47. Sun, Z.H.; Hu, Y.N.; Cheng, H.F. Public Health Risk of Toxic Metal (Loid) Pollution to the Population Living near an Abandoned Small-Scale Polymetallic Mine. Sci. Total Environ. 2020, 718, 137434. [Google Scholar] [CrossRef]
  48. Ma, C.Y.; Cai, D.J.; Yan, H. Soil Cd Pollution and Research Progress of Treatment Techniques. Henan Chem. Ind. 2013, 30, 17–22, (In Chinese with English Abstract). [Google Scholar]
  49. Bhattacharya, P.; Welch, A.H.; Stollenwerk, K.G.; McLaughlin, M.J.; Bundschuh, J.; Panaullah, G. Arsenic in the Environment: Biology and Chemistry. Sci. Total Environ. 2007, 379, 109–120. [Google Scholar] [CrossRef] [PubMed]
  50. Huang, S.S.; Liao, Q.L.; Hua, M.; Wu, X.M.; Bi, K.S.; Yan, C.Y.; Chen, B.; Zhang, X.Y. Survey of Heavy Metal Pollution and Assessment of Agricultural Soil in Yangzhong District, Jiangsu Province, China. Chemosphere 2007, 67, 2148–2155. [Google Scholar] [CrossRef] [PubMed]
  51. Imperato, M. Spatial Distribution of Heavy Metals in Urban Soils of Naples City (Italy). Environ. Pollut. 2003, 124, 247–256. [Google Scholar] [CrossRef]
  52. Moryani, H.T.; Kong, S.Q.; Du, J.K.; Bao, J.G. Health Risk Assessment of Heavy Metals Accumulated on PM2.5 Fractioned Road Dust from Two Cities of Pakistan. Int. J. Environ. Res. Public Health 2020, 17, 7124. [Google Scholar] [CrossRef] [PubMed]
  53. Li, Y.Y.; Wang, H.B.; Wang, H.J.; Yin, F.; Yang, X.Y.; Hu, Y.J. Heavy Metal Pollution in Vegetables Grown in the Vicinity of a Multi-Metal Mining Area in Gejiu, China: Total Concentrations, Speciation Analysis, and Health Risk. Environ. Sci. Pollut. Res. Int. 2014, 21, 12569–12582. [Google Scholar] [CrossRef]
  54. Hu, B.F.; Xue, J.; Zhou, Y.; Shao, S.; Fu, Z.Y.; Li, Y.; Chen, S.C.; Qi, L.; Shi, Z. Modelling Bioaccumulation of Heavy Metals in Soil-Crop Ecosystems and Identifying Its Controlling Factors Using Machine Learning. Environ. Pollut. 2020, 262, 114308. [Google Scholar] [CrossRef] [PubMed]
  55. Zhou, H.; Yang, W.T.; Zhou, X.; Liu, L.; Gu, J.F.; Wang, W.L.; Zou, J.L.; Tian, T.; Peng, P.Q.; Liao, B.H. Accumulation of Heavy Metals in Vegetable Species Planted in Contaminated Soils and the Health Risk Assessment. Int. J. Environ. Res. Public Health 2016, 13, 289. [Google Scholar] [CrossRef] [Green Version]
  56. Hu, W.Y.; Huang, B.; Tian, K.; Holm, P.E.; Zhang, Y.X. Heavy Metals in Intensive Greenhouse Vegetable Production Systems Along Yellow Sea of China: Levels, Transfer and Health Risk. Chemosphere 2017, 167, 82–90. [Google Scholar] [CrossRef] [PubMed]
  57. Li, B.; Wang, Y.H.; Jiang, Y.; Li, G.C.; Cui, J.H.; Wang, Y.; Zhang, H.; Wang, S.C.; Xu, S.; Wang, R.Z. The Accumulation and Health Risk of Heavy Metals in Vegetables around a Zinc Smelter in Northeastern China. Environ. Sci. Pollut. Res. Int. 2016, 23, 25114–25126. [Google Scholar] [CrossRef] [PubMed]
  58. Chen, H.; Huang, Q.; Liu, L.; Cai, P.; Liang, W.; Li, M. Poultry Manure Compost Alleviates the Phytotoxicity of Soil Cadmium: Influence on Growth of Pakchoi (Brassica Chinensis L.). Pedosphere 2010, 20, 63–70. [Google Scholar] [CrossRef]
  59. Vandecasteele, B.; Willekens, K.; Laing, G.D.; Waes, J.V.; Tack, F.M.G. Effect of Feedstock, Organic Matter Content and pH on Cd, Zn and Mn Availability in Farm Compost Based on Bark and Wood Chips. Acta Hortic. 2014, 1018, 661–668. [Google Scholar] [CrossRef]
  60. Yang, Y.; Chen, W.; Wang, M.; Peng, C. Regional Accumulation Characteristics of Cadmium in Vegetables: Influencing Factors, Transfer Model and Indication of Soil Threshold Content. Environ. Pollut. 2016, 219, 1036–1043. [Google Scholar] [CrossRef]
  61. Clemens, S. Toxic Metal Accumulation, Responses to Exposure and Mechanisms to Tolerance in Plants. Biochimie 2006, 88, 1707–1719. [Google Scholar] [CrossRef]
  62. Harrison, R.M.; Chirgawi, M.B. The Assessment of Air and Soil as Contributors of Some Trace Metals to Vegetable Plants. I. Use of a Filtered Air Growth Cabinet. Sci. Total Environ. 1989, 83, 13–34. [Google Scholar] [CrossRef]
  63. Fan, Y.; Li, Y.L.; Li, H.; Cheng, F.Q. Evaluating Heavy Metal Accumulation and Potential Risks in Soil-Plant Systems Applied with Magnesium Slag-Based Fertilizer. Chemosphere 2018, 197, 382–388. [Google Scholar] [CrossRef]
  64. Ji, Y.; Wu, P.; Zhang, J.; Zhang, J.; Zhou, Y.; Peng, Y.; Zhang, S.; Cai, G.; Gao, G. Heavy Metal Accumulation, Risk Assessment and Integrated Biomarker Responses of Local Vegetables: A Case Study Along the Le’an River. Chemosphere 2018, 199, 361–371. [Google Scholar] [CrossRef]
  65. Zhou, Q.; Yang, Y.C.; Yang, Z.Y. Molecular Dissection of Cadmium-Responsive Transcriptome Profile in a Low Cadmium-Accumulating Cultivar of Brassica Parachinensis. Ecotox. Environ. Saf. 2019, 176, 85–94. [Google Scholar] [CrossRef]
  66. USEPA. Risk Based Concentration Table; United States Environmental Protection Agency: Washington, DC, USA, 2000.
  67. USEPA. Risk Assessment Guidance for Superfund. Human Health Evaluation Manual, (Part a); Office of Emergency and Remedial Response [Epa/540/1-89/002]; United States Environmental Protection Agency: Washington, DC, USA, 1989; Volume I.
  68. USEPA. Exposure Factors Handbook; Office of Research and Development, National Center for Environmental Assessment [Epa/600/P-95/002f a-C]; United States Environmental Protection Agency: Washington, DC, USA, 1997.
  69. Tao, X.Q.; Shen, D.S.; Shentu, J.L. Bioaccessibility and Health Risk of Heavy Metals in Ash Fromthe Incineration of Different E-Waste Residues. Environ. Sci. Pollut. Res. Int. 2015, 22, 3558–3569. [Google Scholar] [CrossRef]
  70. White, P.D.; Marcus, A.H. The Conceptual Structure of the Integrated Exposure Uptake Biokinetic Model for Lead in Children. Environ. Health Persp. 1998, 106, 1513. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Cao, S.Z.; Duan, X.L.; Zhao, X.G.; Wang, B.B.; Ma, J.; Fan, D.L.; Sun, C.Y.; He, B.; Wei, F.S.; Jiang, G.B. Health Risk Assessment of Various Metal(Loid)S Via Multiple Exposure Pathways on Children Living near a Typical Lead-Acid Battery Plant, China. Environ. Pollut. 2015, 200, 16–23. [Google Scholar] [CrossRef] [PubMed]
  72. Dhankher, O.P.; Pilon-Smits, E.A.; Meagher, R.B.; Doty, S. Biotechnological Approaches for Phytoremediation. Plant Biotechnol. Agric. 2012, 2012, 309–328. [Google Scholar]
  73. Vahter, M. Species Differences in the Metabolism of Arsenic Compounds. Appl. Organomet. Chem. 1994, 8, 175–182. [Google Scholar] [CrossRef]
  74. Zhou, X.; Yoshida, K.; Kuroda, K.; Endo, Y.; Endo, G. Effects of Cysteine on the Cytotoxicity of Arsenic Compounds. Arch. Environ. Con. Tox. 2003, 45, 324–330. [Google Scholar] [CrossRef]
  75. Bakkali, K.; Martos, N.R.; Souhail, B.; Ballesteros, E. Determination of Heavy Metal Content in Vegetables and Oils from Spain and Morocco by Inductively Coupled Plasmass Spectrometry. Anal. Lett. 2012, 45, 907–919. [Google Scholar] [CrossRef]
  76. Pajevic, S.; Arsenov, D.; Nikolic, N.; Borisev, M.; Orcic, D.; Zupunski, M.; Mimica-Dukic, N. Heavy Metal Accumulation in Vegetable Species and Health Risk Assessment in Serbia. Environ. Monit. Assess. 2018, 190, 459. [Google Scholar] [CrossRef] [PubMed]
  77. Hussain, R.; Luo, K.L.; Liang, H.D.; Hong, X.P. Impact of the Coal Mining-Contaminated Soil on the Food Safety in Shaanxi, China. Environ. Geochem. Health 2019, 41, 1521–1544. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Locations of the investigated areas and sampling sites in Daye City.
Figure 1. Locations of the investigated areas and sampling sites in Daye City.
Ijerph 18 02617 g001
Figure 2. Boxplot of pollution index (a) and risk index (b) for heavy metals in the four zones.
Figure 2. Boxplot of pollution index (a) and risk index (b) for heavy metals in the four zones.
Ijerph 18 02617 g002
Figure 3. Spatial distribution pollution levels of Igeo-Cd in the four zones: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine.
Figure 3. Spatial distribution pollution levels of Igeo-Cd in the four zones: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine.
Ijerph 18 02617 g003
Figure 4. Heavy metals in edible parts of vegetable and crop species in the four zones: (a) Cu; (b) Cr; (c) Pb; (d) Cd; (e) As.
Figure 4. Heavy metals in edible parts of vegetable and crop species in the four zones: (a) Cu; (b) Cr; (c) Pb; (d) Cd; (e) As.
Ijerph 18 02617 g004
Figure 5. Correlation plots of heavy metals in the four zones: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine; * p < 0.05, ** p < 0.01.
Figure 5. Correlation plots of heavy metals in the four zones: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine; * p < 0.05, ** p < 0.01.
Ijerph 18 02617 g005
Figure 6. Target hazard quotients (THQ) and hazard index (HI) of different vegetable and crop types: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine.
Figure 6. Target hazard quotients (THQ) and hazard index (HI) of different vegetable and crop types: (a) Non-ferrous metals smelter; (b) Abandoned copper mine; (c) Limestone quarry; (d) Iron mine.
Ijerph 18 02617 g006
Table 1. Vegetable species used in the experiments.
Table 1. Vegetable species used in the experiments.
NO.TypesVegetable SpeciesHarvest Season
1Solanaceous fruitsTomato (Lycopersicum esculentum Mill.)Spring and summer
2 Bell pepper (Capsicum annuum L.)Spring and summer
3 Cucumber (Cucumis sativus L.)Spring and summer
4 Eggplant (Solanum melongena L.)Spring and summer
5Leafy vegetablesEdible amaranth (Amaranthus tricolor L.)Spring and summer
6 Lettuce (Lactuca sativa L.)Spring and summer
7 Pakchoi (Brassica Chinensis L.)Spring and summer
8 Indian Lettuce (Lactuca sativa L. var longifoliaf. Lam)Spring and summer
9 Brassica campestris (Brassica campestris L. var. purpurea Baileysh)Spring and summer
10NutsPeanut (Arachis hypogaea Linn.)Autumn and winter
11Rhizome geophytesIpomoea (Ipomoea batatas L.)Autumn and winter
12 White radish (Raphanus sativus L.)Autumn and winter
13CerealsMaize (Zea mays L.)Autumn and winter
14LegumesAsparagus bean (Vigna unguiculata Linn.)Spring and summer
15Soya bean (Glycine max Linn.)Autumn and winter
Table 2. Characteristics of heavy metals in soils of the four typical zones.
Table 2. Characteristics of heavy metals in soils of the four typical zones.
Sampling Areas pHCu (mg/kg)Cr (mg/kg)Pb (mg/kg)Cd (mg/kg)As (mg/kg)
Non-ferrous metals smelter (n = 21)Range6.51–6.87116.17–858.9761.33–94.8227.04–219.510.44–7.514.49–75.54
Mean ± SD6.71 ± 0.13475.39 ± 226.1072.27 ± 7.55136.89 ± 59.993.80 ± 2.1042.48 ± 21.26
CV/%1.9447.5610.4543.8255.2650.05
K-Sp0.200 *0.200 *0.0060.190 *0.144 *
a Excessive rate %100.000.0061.90100.0071.4
Abandoned copper mine (n = 26)Range6.58–7.35213.36–1446.1066.21–131.7033.25–223.480.15–3.284.09–86.04
Mean ± SD7.03 ± 0.30729.39 ± 383.8395.69 ± 19.37107.76 ± 56.641.34 ± 0.7841.96 ± 21.68
CV/%4.2552.6220.2452.5658.2151.67
K-Sp0.062 *0.200 *0.000 0.200 *0.200 *
a Excessive rate %100.000.0042.3188.4653.85
Limestone quarry (n = 20)Range7.56–7.8374.53–123.755.10–99.7828.01–103.250.40–2.062.93–45.98
Mean ± SD7.72 ± 0.1194.45 ± 14.9355.00 ± 31.2260.57 ± 19.771.10 ± 0.4526.45 ± 11.20
CV/%1.4115.8156.7632.6440.9142.34
K-Sp0.134 *0.200 *0.200 *0.200 *0.200 *
a Excessive rate %30.000.000.0085.0060.00
Iron mine (n = 23)Range4.76–5.4220.17–64.4849.12–80.7736.09–76.250.05–1.750.79–39.24
Mean ± SD5.05 ± 0.2236.31 ± 11.0258.23 ± 7.7252.34 ± 13.530.45 ± 0.418.33 ± 9.49
CV/%4.2930.3513.2625.8591.11113.93
K-Sp0.125 *0.0050.011 0.0210.000
a Excessive rate %8.700.0021.7447.830.00
b Background values 30.78626.70.1712.3
c Soil pollutant standards 50, 100150, 200, 25070, 90, 120, 1700.3, 0.640, 30, 25
Total excessive rate % 61.110.0032.2280.00 45.56
K-Sp is the significance level of Kolmogorov-Smirnov test for normality, * p > 0.05. a Refer to soil pollutant standards. b Soil background values in Hubei province (Soil background values of Chinese elements 1990). c Soil environment quality: Risk control standard for soil contamination of agricultural land (GB15618-2018 Ministry of Environment Protection).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Bao, J.; Wang, T.; Moryani, H.T.; Kang, W.; Zheng, J.; Zhan, C.; Xiao, W. 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. https://doi.org/10.3390/ijerph18052617

AMA Style

Wang Z, Bao J, Wang T, Moryani HT, Kang W, Zheng J, Zhan C, Xiao W. Hazardous Heavy Metals Accumulation and Health Risk Assessment of Different Vegetable Species in Contaminated Soils from a Typical Mining City, Central China. International Journal of Environmental Research and Public Health. 2021; 18(5):2617. https://doi.org/10.3390/ijerph18052617

Chicago/Turabian Style

Wang, Zhen, Jianguo Bao, Tong Wang, Haseeb Tufail Moryani, Wei Kang, Jin Zheng, Changlin Zhan, and Wensheng Xiao. 2021. "Hazardous Heavy Metals Accumulation and Health Risk Assessment of Different Vegetable Species in Contaminated Soils from a Typical Mining City, Central China" International Journal of Environmental Research and Public Health 18, no. 5: 2617. https://doi.org/10.3390/ijerph18052617

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop