Next Article in Journal
The Impact of Metal-Based Nanoparticles Produced by Different Types of Underwater Welding on Marine Microalgae
Previous Article in Journal
Self-Reported Assessment of the Socio-Economic Impact of Anticancer Chemotherapy-Related Neurotoxicity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Risk Assessment of Heavy Metal in Farmlands and Crops Near Pb–Zn Mine Tailing Ponds in Niujiaotang, China

1
College of Resources and Environmental Engineering, Guizhou University, Guiyang 550025, China
2
Key Laboratory of Karst Georesources and Environment, Ministry of Education, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Toxics 2023, 11(2), 106; https://doi.org/10.3390/toxics11020106
Submission received: 17 December 2022 / Revised: 16 January 2023 / Accepted: 19 January 2023 / Published: 22 January 2023
(This article belongs to the Section Metals and Radioactive Substances)

Abstract

:
To accurately evaluate the pollution and risk of heavy metals in crops and farmlands near mines, we determined the contents of Cr, Ni, Cu, As, Cd, Pb, and Zn in 10 farmland soil sampling sites and six crops (pak choi, rice, spring onion, radish, Chinese cabbage, Chrysanthemum coronarium) in an area near the Niujiaotang Pb–Zn mine in Duyun City, China. Four evaluation methods were compared, including the potential ecological hazard index, Nemeiro comprehensive pollution assessment, risk assessment code, and the ratio of secondary phase to primary phase methods. The average concentration of As, Cd, Pb, and Zn exceeded the soil environmental background levels in Niujiaotang and Guizhou Province. Cd exceeded the standard substantially, and Zn pollution accumulation was the most evident. Heavy metal contamination of crops was in the order pak choi > Chinese cabbage > spring onion > paddy > radish > Chrysanthemum coronarium, whereas heavy metal concentration in crops were in the order Zn > As > Cr > Cd > Ni > Pb. The levels of all heavy metals except Cu exceeded Chinese food hygiene standards. Carcinogenic and non-carcinogenic chemicals in crops present significant risks to adults and children. Risk evaluation considering the morphological contents of heavy metals rather than their total concentration was more accurate for environmental quality assessment of agricultural soils. Samples should be collected at different times to study the spatial and temporal distribution, and further studies on the migration transformation of heavy metals between the tailings pond-soil-crop should be conducted.

1. Introduction

China is one of the world’s largest producers and consumers of Pb and Zn. The cumulative production of Pb and Zn in China until 2007 was approximately 6.69 and 12.59 Mt, respectively, and approximately 1.62 Mt of Pb and 3.32 Mt of Zn were released into the environment due to mining, processing, and smelting activities. During mining and smelting, heavy metals can be transported to the surrounding environment through atmospheric deposition, surface water contamination (including acid mine drainage), and groundwater seepage, thereby posing serious ecological risks [1,2]. Guizhou Province is famous for its rich mineral resources and long history of nonferrous metal mining and smelting [3]. Despite the economic benefits of mining, this activity causes heavy metal pollution of nearby soils and crops. Heavy metal pollution is insidious, long-term, and irreversible, and heavy metals that enter the soil cannot be degraded and thus, they persist in the soil for a long time and accumulate, eventually causing harm to the ecological environment and human health via food chain contamination [4,5].
The Niujiaotang deposit is located in the north-south tectonic deformation zone of Guiding at the southwest edge of the Yangzi Plate. It is a Cd-rich Zn deposit and has the highest Cd content in Guizhou Province [6]. The Pb–Zn mine tailings in Niujiaotang are large in volume and mostly stored in constructed tailing ponds. They are mainly composed of gray, off-white, and gray-brown Pb–Zn fine-grained tailings [7]. During the natural accumulation of tailings, influenced by atmospheric rainfall, the tailings sand will undergo epigenetic geochemical reactions with carbon dioxide and oxygen in the air in the water–rock system, which activates the heavy metals in the tailings sand [8,9]. The Cd content in the Niujiaotang Pb–Zn ore in Duyun, Guizhou is unusually high, generally ranging from 2284 to 9850 μg/g, with a maximum of 13,400 μg/g, and its reserves have reached the scale of large deposits (medium-large Zn and Cd reserves greater than 5000 t [10]. Sun et al. analyzed the contamination characteristics of different heavy metals (Pb, Zn, Cr, Cu, Cd, and Hg) and the Pb isotopic composition in natural soil (layers A and C) in the Shuikou Mountain Pb–Zn mining area (and nearby areas) of Hunan, a typical Pb–Zn mining area in China. Their results showed that the soil in layer A was clearly contaminated with heavy metals owing to the influence of Pb–Zn beneficiation and smelting activities. Pan et al. reported that the Niujiaotang Cd-rich Zn deposit is a polymetallic sulfide deposit with high Zn and Cd contents in its ore.
In recent years, heavy metal contamination of agricultural soils and crops near mining areas and the associated risks have attracted widespread attention from various sectors and researchers worldwide. It is important to analyze and understand the morphological characteristics of heavy metals in the soil to identify their sources and determine their biological effectiveness [11]. The migration capacity of heavy metals is closely related to their mobile forms, which affect their cycling and ecotoxicity [12]. For instance, heavy metals in the exchangeable state have high mobility and are easily absorbed by organisms, and this is the most unstable and toxic form, which causes direct toxicity. The reducible and oxidizable states can be transformed into biologically effective forms in reducing and oxidizing environments, respectively, and their potential harm to the environment is known as potential toxicity. The residual state of heavy metals is usually inert and does not participate in chemical reactions, does not change with modifications in the external environment, and has low bioavailability [13]. The results of this study indicated that the total content of heavy metals in the soil does not completely determine their environmental behaviors and biomorphic effects, and plants do not fully absorb and utilize all heavy metals from the soil, only forms in the biologically effective state (mainly the exchangeable and carbonate-bound states) are absorbed. In contaminated areas, the harm to humans upon consumption of heavy metal-contaminated rice is even higher than that caused by ingestion of heavy metal-contaminated water, and long-term consumption of the former can easily lead to chronic poisoning [14].
In this study, we aimed to elucidate the distribution of heavy metals in farmlands and crops near the Duyun Niujiaotang Pb–Zn mine through sampling and experimental analysis. The pollution level, pollution risk, and crop health risk were evaluated based on the total concentration and distribution of heavy metals. We expect our results to contribute to risk assessments of soils and crops near mines and the development of prevention and mitigation measures to ensure food security, human health, and the sustainability of mining areas.

2. Materials and Methods

2.1. Study Area

The Niujiaotang Cd-rich Pb–Zn deposit is located in Duyun City, Guizhou Province (latitude: 26.226643; longitude: 107.666317), which has a subtropical monsoon humid climate, moderate cold in winter and a coldest average daily temperature of 5.6 °C in January. There is no severe summer heat, and the average daily temperature in the hottest month (July) is 24.8 °C. Rainfall is heavy, with an annual average of 1431.1 mm. Rain and heat occur in the same season. The annual average temperature is 16.1 °C and the frost-free period is approximately 300 days. The soil in the study area is generally severely polluted with heavy metals such as Zn and Pb and the heavy metals in the main agricultural products substantially exceed the standards, posing a serious health threat to nearby residents.
According to a preliminary investigation, Bagu town has rich mineral resources, and the main known deposits are of Pb–Zn, limonite, and barite. The Fanjia River Pb–Zn ore in the Mapo Pb–Zn section with shallow metal reserves of 500,000 t belongs to the national medium-sized deposits, with grades of 0.84–11.3%, whereas the geological reserves of the Pb–Zn ore in the Qingshui pond section and Duniu Pb–Zn section amount to approximately 200,000 t, with grades of 0.57–10.8%. There are nine processing plants and 34 mining holes in the town, with a daily output of 20,000 t of ore, producing 0.08 million t of concentrated powder. The currently known reserves of this deposit are more than 350,000 t of Zn (Zn grade of 5.85–24.48%) and 5299.14 t of Cd (large scale or above) [6].

2.2. Sample Collection and Preprocessing

The sampling points were located along the Caiyuan River in farmland areas near the Pb–Zn mine, and the samples were collected during an abundant water period (August 2021) and a flat-water period (November 2021). In total, seven surface water samples and 10 crop samples (including pak choi, rice, spring onion, radish, Chinese cabbage, and Chrysanthemum coronarium) were collected, along with 10 inter-root soil samples. When collecting soil samples, a 3 × 3 m² area was selected and five soil samples were taken by the plum method and mixed into one point sample. Ultrapure water was used in the experiments, and all reagents were of analytical grade. The crop and soil samples were placed into self-sealing bags and water samples were collected into polyethylene bottles, which were all transported to the laboratory. To prepare the samples, the crop samples were washed with tap water after manual shaking to remove large dirt particles on the surface. The crop samples were divided into four parts: leaves, stems, roots, and fruits. All samples were rinsed with deionized water and placed into an oven at 105 °C for 30 min, dried at 75 °C to constant weight, crushed, ground through a 100-mesh sieve, and then set aside. Subsequently, we removed stones, biological debris, and plant fragments from the soil samples, air-dried them in a cool and dry environment, ground them, and passed them through 60-mesh and 200-mesh sieves. Finally, the samples were then bagged and sealed. The surface water samples were filtered through a 0.45-μm membrane, acidified with ultrapure nitric acid to pH < 2 and stored in sealed polyethylene plastic bottles at 4 °C until analysis. The pH, Eh, and EC were measured by Multi3630, WTW, Germany. Reagent blanks, three parallel soil samples and the national standard soil sample GSS-37 were used for quality control during the analysis. The correlation coefficients of the calibration curves of the heavy metal standard solutions were all greater than 0.999, and the experiments were conducted in three replicates. The model of the oven was DHG-9140A produced by Shanghai Huitai Instrument Manufacturing Co. The limits of detection (LODs) were 0.01 mg L−1 for Cd, 0.001 mg L−1 for As, 0.005 mg L−1 for Pb, 0.005 mg L−1 for Ni, 0.02 mg L−1 for Zn, 0.01 mg L−1 for Cu, and 0.005 mg L−1 for Cr. The sampling locations are shown in Figure 1.

2.3. Sample Analysis

Nitric acid-hydrofluoric acid was added to the soil samples, ratio 3:1 (v/v) for high-temperature digestion. The concentrations and chemical forms of heavy metals (Cr, Ni, Cu, As, Cd, Pb, and Zn) in soil and crop solutions were determined by inductively coupled plasma-mass spectrometry (ICP-MS). A procedural blank, parallel experiment, and national standard soil sample (GSS-37) were used for quality control during the analytical process. Ultrapure water was used in the experiments, and all reagents were of analytical grade. The vessels were soaked in 20% HNO3 solution for more than 24 h, washed with ultrapure water, and then dried.
Determination of the chemical forms of heavy metal elements is an important basis for the study of cyclic patterns such as elemental migration and transformation. The extraction steps were conducted according to the sequential extraction scheme of the European Community Standards Agency (BCR sequential extraction procedure), as follows. The samples were sequentially leached under acetic acid solution (CH3COOH) to isolate the weak acid-extracted state (F1), hydroxylamine hydrochloride solution (NH2OH·HCl) to isolate the reducible state (F2), and ammonium acetate (NH4Ac) to isolate the oxidizable state (F3). The total concentration of heavy metals and the contents of the first three forms were determined by ICP-MS (Thermo Fisher). The sum of the first three forms was subtracted from the total concentration to obtain the content of the residual state, which was denoted as F4.

2.4. Pollution and Ecological Risk Assessment

There are numerous methods to evaluate soil heavy metal risk, the most robust and widely used being the index evaluation method [15]. This method is based on the total concentrations of heavy metals, and it can visually reflect the correlation between measured and background concentrations of heavy metals to evaluate their risk in the soil. According to their evaluation criteria, index evaluation methods can be divided into single pollution [16], Nemeiro [17], geological accumulation [18], and potential ecological risk index methods [19]. Ecological risk evaluation, which is based on the total amount of soil heavy metals, often does not effectively represent the chemical activity and bioavailability of heavy metals. In contrast, evaluation based on the morphology of heavy metals can more realistically predict the ecological risk of soil heavy metals and provide a more scientific basis for pollution prevention and control [20]. Therefore, morphology-based evaluation methods, including the commonly used risk assessment code (RAC) method [21] and ratio of secondary phase to primary phase (RSP) method [22], have become more popular. To select the method that can most objectively and comprehensively reflect the pollution level and ecological risk of soil heavy metals, different evaluation methods should be compared, and the environmental effects and behavioral characteristics of heavy metals should be comprehensively evaluated.

2.4.1. Potential Ecological Risk Index (RI) Method

The potential ecological RI method is an internationally recognized method widely used for the evaluation of soil heavy metal pollution [23]. It combines environmental chemistry, biotoxicology, and ecology to quantitatively classify the degree of potential heavy metal hazard [24]. The RI can be calculated as follows:
  C f i = C s i C n i    
E r i = T r   i × C f i
R I = i = 1 n E r i = i = 1 n T r   i × C f i = i = 1 n T r   i × C s i C n i   ,
where Cif is the contamination coefficient of heavy metal i relative to the reference value; Cis is the measured concentration of heavy metal i; Cin is the reference value for heavy metal i, which is the background value in Niujiaotang; Eir is the environmental risk index of heavy metal i; Tir is the toxicity response coefficient of heavy metal i [25]; and RI is the multi-element environmental risk composite index. Reference and toxicity coefficient for various heavy metals in Table 1.

2.4.2. Nemeiro Comprehensive Pollution Assessment

The Nemeiro comprehensive pollution index (PN) highlights the measured degree of soil contamination by the maximum contamination to reflect the overall contamination status of the soil [26]. It is calculated as follows:
P N = P i m a x 2 + P i a v g 2 2
where Pimax and Piavg are the maximum and average pollution index values of heavy metal i in the soil, respectively.

2.4.3. RAC Method (Risk Assessment Code)

The RAC is mainly related to F1, and it is calculated as follows:
R A C = C F 1 C T × 100 % ,
where RAC is the mass fraction of F1 (weak acid extracted state) in the total mass, %; CF1 is the concentration of F1 (weak acid extracted state), mg·kg−1; and CT is the sum of the four BCR forms, mg·kg−1 [22].

2.4.4. RSP Method (Ratio of Secondary Phase to Primary Phase)

The RSP method divides the soil into primary and secondary phases, and it evaluates the degree of soil heavy metal contamination by calculating the ratio of the secondary phase (weathering products of primary minerals and exotic secondary materials) to the primary phase [27]. The RSP is calculated as follows:
R S P = M   s e c   M   p r i m   ,
where Msec is the content of heavy metals in the secondary phase, mg·kg−1; and Mprim is the residual state (primary phase) content, mg·kg−1.
The above four soil heavy metal risk evaluation methods are widely used in China and other countries, and the industry-accepted evaluation criteria for each method are detailed in Table 2.

2.5. Human Health Risk Assessment

The human health risks posed by heavy metals in agricultural products were evaluated using internationally recognized chemical carcinogen and chemical non-carcinogen risk evaluation models. According to the International Agency for Research on Cancer (IARC), the EPA Integrated Risk Information System (IRIS) database, and the classification system compiled by the World Health Organization (WHO), Cd and As are chemical carcinogens, whereas Pb and Zn are non-carcinogens.
The human health risk for a carcinogenic chemical can be calculated as follows:
R i = A D D i × S F i
where Ri is the individual carcinogenic annual risk attributed to carcinogenic chemical i, a−1; ADDi is the average daily exposure dose to the carcinogenic chemical i, mg·(kg·d)−1; and SFi is the carcinogenic intensity factor of i, (mg/kg/day)−1. The grading criteria are as follows: Ri < 1 × 10−6 y−1 indicates negligible risk to human health; 1 × 10−6 < Ri < 1×10−4 a−1 indicates acceptable risk; and Ri > 1 × 10−4 a−1 indicates significant risk to human health.
The human health risk for non-carcinogenic chemicals can be calculated as follows:
H = A D D i R f D ,
where H is the individual non-carcinogenic annual risk attributed to non-carcinogenic chemicals, y−1; and RfD is the reference dose of non-carcinogenic chemical exposure, mg·(kg·d)−1. The grading criteria are as follows: H ≤ 1 indicates that the ingestion of crop products does not pose a health hazard; 1 < H ≤ 10 indicates a high health hazard; and H > 10 indicates a risk of chronic poisoning [28].
ADDi was calculated as follows:
A D D i = C i × I R × E F × E D B W × A T   ,
where Ci is the content of heavy metal i in a crop sample, mg/kg; IR is the daily crop consumption per capita (kg/d); EF is the annual exposure days to heavy metals, d/y; ED is the average exposure period, y; BW is the average human body weight, kg; AT is the life expectancy, d1. The parameters for the health risk evaluation model in Equations (7)–(9) are shown in Table 3.

3. Results and Discussion

The physicochemical properties of soils affect the biological effectiveness of heavy metals mainly by influencing their morphology [35]. In this context, pH is one of the most important soil physicochemical properties because it affects the activity of heavy metals. Under acidic conditions, heavy metals are more mobile and more easily absorbed by plants [36]. Figure 2b shows that the average soil pH in the study area was 7.298, and that the soil was alkaline in all sample collection points, except point 3. The exposed substratum in the study area is Cambrian Balang Formation shale with a small amount of silt-fine sandstone and a thin layer of tuff in the upper part [37]. Different parent rocks lead to different soil physicochemical properties. Therefore, the difference in parent rocks was considered the main reason for the pH difference between the sampling points. Soil organic matter ranged from 2.822% to 7.572% in Figure 2a. Table 4 shows that the concentrations of heavy metals in the water bodies did not exceed the national surface water quality standard for Class V. From Figure 3, it can be seen that the average value of Eh in the water samples was −32.9 mV (the range of Eh was −8.6 to −69.4); the average value of EC was 1096.571 S·cm–1 (the range of EC was 275 to 1487); the average value of pH of the water samples was 7.59 (the range of pH was 7.12 to 8.28), which was basically weakly alkaline. Pan showed that the reason for the alkaline nature of the water in the Oxnard tailing pond mine is mainly related to the neutralizing and buffering effect of carbonate rock on the acidic wastewater produced by sulfide ores [38].

3.1. Heavy Metal Contents and Geochemical Characteristics of Agricultural Soils

The contents of seven heavy metals in agricultural soils in the study area are presented in Table 5 and Figure 4. The heavy metal contents in soil samples from different sampling points varied widely. The average concentrations (and ranges) of As, Cd, Pb, and Zn were 36.096 (19.581–75.149) mg/kg, 20.980 (3.145–58.748) mg/kg, 210.503 (73.679~796.233) mg/kg, 1730.283 (788.423~2382.207) mg/kg, which were 2.927, 91.216, 9.144, and 15.346 times the background values of soils in Niujiaotang, and 1.805, 31.787, 5.980, and 17.390 times the background values in Guizhou Province, respectively. Among the heavy metals, Cd exceeded the standard substantially, and Zn, which was the most abundant heavy metal, presented the greatest accumulation.
The content and distribution of heavy metals in soils are related to their sources. Soils near tailings ponds are often contaminated via atmospheric deposition, ground diffusion, and infiltration. Pb–Zn smelting and the associated wastewater sludge are the main sources of As, Cd, Pb, and Zn contamination to nearby areas [39]. The sampling sites in this study were distributed near Pb–Zn mines, and the ancient Zn refining method used in the early years of the Guizhou Zn mines produced pollutants such as Cd, Pb, and Zn, which were released into the atmosphere. In addition, exhaust from motor vehicles on the roads near tailings ponds increases the atmospheric content of heavy metals, such as Cd, Pb, and Zn, which can enter the soil environment through atmospheric deposition [40]. Guizhou is located in the southwest karst region, with a warm and rainy climate and high leaching potential due to weathering. The weathering of rocks releases heavy metals and as such can be a major heavy metal source [41]. Moreover, heavy metals in cultivated soils exceeded the standard values, whereas the concentrations in natural soils were below the background values, which indicated that agricultural cultivation also is a major source of soil heavy metals. Therefore, the soil pollution with As, Cd, Pb, and Zn was attributed to Pb–Zn mining activities, transportation, parent rock weathering, and agricultural farming.

3.2. Soil–Plant Heavy Metal Migration and Enrichment Characteristics

3.2.1. Morphological Distribution of Heavy Metals in Agricultural Soils

Figure 5 shows the average morphological distribution of each element based on morphological analysis of 10 BCR sequentially extracted rhizospheric soil samples. The extraction steps of the chemical forms of heavy metal elements were conducted according to the sequential extraction scheme of the European Community Standards Agency (BCR sequential extraction procedure). The residual state of Zn was the lowest (68.248%), followed by Cd (88.717%). In contrast, the residual states of Pb, As, Cu, Ni, and Cr were all above 90%, which indicated poor migration ability, that is, they are not easily absorbed by organisms and, therefore, are less harmful to the environment. The weak acid extractable states of Zn and Cd accounted for 10.358% and 5.442%, respectively, whereas those of the other elements were less than 1%. The oxidizable state of Zn was 11.737%, whereas those of the other elements were less than 4%. Zn showed various morphological distribution characteristics. The weak acid extractable, reducible, and oxidizable states accounted for similar proportions, with averages of 10.358%, 9.657%, and 11.737%, respectively. The smaller proportion of the reducible state might be related to the redox environment of the soil. [42] have shown that changes in soil redox conditions may cause oxidizable and reducible Zn to be released from the soil. The lowest Cd content was observed for the oxidizable state, which might be related to the low uptake of Cd by organic matter [43].
The biologically effective state of heavy metals is the sum of their first three forms, which are prone to migratory transformation and re-release into the environment and, therefore, pose serious ecological hazards upon plant or animal intake [44]. The order of biologically effective state contents of the heavy metals in this study was Zn (31.75%) > Cd (11.28%) > Pb (7.99%) > Cu (3.02%) > Ni (2.91%) > Cr (1.55%) > As (0.55%). Among them, Zn and Cd presented relatively high proportions of the biologically effective state (>10%), with a high migration capacity and potential hazard. Therefore, in agricultural production, Zn and Cd are more easily absorbed and used by plants.
The content of weak acid-extractable Cd increased with the increase in overall Cd concentration. Nevertheless, pH and other soil physicochemical features affected the Cd content to a higher extent, which is consistent with results reported by [45]).

3.2.2. Heavy Metal Pollution of Crops

Since the farmlands in the study area were contaminated with heavy metals, the crops growing on them also had a potential risk of heavy metal contamination. Table 6 shows that the exceedance rates of heavy metals (Cr, Ni, As, Cd, Pb, Zn) in pak choi near the mine area ranged from 11.015% to 49.985% (refer to GB2762-2017, GB15199-1994, and GB13106-1991). The exceedance rates of the heavy metals in rice ranged within 0.152–48.042%, whereas those in Chinese cabbage ranged from 0.473% to 112.433%. Among the heavy metals, the rate of Cu was exceeded only in Chrysanthemum coronarium and pak choi in points 5 and 6, indicating that these crops have a strong Cu enrichment ability. Pb was exceeded in all crops in all points, except for radish in point 10. The concentration of Zn, Cd, and Pb in pak choi were significantly higher than those in the other five crops, and the exceedance rate of Zn in pak choi in point 3 was higher than that in all other crops.
The results showed that the overall exceedance rate in the crops was in the order rice > pak choi > Chinese cabbage > spring onion > Chrysanthemum coronarium > radish. Zn had the highest concentration in all six crops, and its average content in descending order was: pak choi > Chinese cabbage > spring onion > rice > radish > Chrysanthemum coronarium, which may be related to the growth habits of the crops. Zn is an essential plant element, which can explain the highest effective content of Zn and its easy absorption by plants [46]. The difference between the mean contents of Cr and As was not large, and As > Cr for all crops except pak choi in sites 5 and 6 and rice in site 2. The effective state of As accounted for the lowest proportion of total As, but the content of As in the crops was higher than that of Cr, indicating that the crops in this study presented a higher absorption rate for As than for Cr. Moreover, Cd exceeded the standard in all crops because it is easily absorbed by crop plants [47].
There was no clear pattern in the concentration of Cu, Cr, and As in the different crops. The generally low content of Ni was attributed to its low content in the soils. The mean concentration of Pb was high in soils, but generally low in crops because Pb has relatively weak mobility and does not easily accumulate in plants [48,49]. The Cr concentration was generally low in soils, but often exceeded the standard in crops, especially in pak choi, and this was attributed to its strong soil mobility [50,51]. Contamination with heavy metals such as Cd, Zn, and As can substantially inhibit the growth and development of crop plants and animals, and long-term ingestion by local residents through dietary exposure can pose a considerable health risk [52]. Therefore, the cultivation of pak choi and rice, which presented high heavy metal enrichment capacity, should be avoided in the area, and radish and spring onion should be given preference so as to reduce the human health risks.

3.2.3. Enrichment and Transport of Heavy Metals

The enrichment factor is an indicator of heavy metal enrichment in crops, and the ratio of the heavy metal content in the edible parts of a crop to that in the soil is used to analyze the soil heavy metal uptake and cumulative effect of the crop [53]. Larger enrichment factor values indicate that crops are more likely to uptake the heavy metal from the farmland, that is, they present greater biological effectiveness [54]. The metal enrichment coefficients of the edible part of rice, a typical crop near the tailings area, followed the order Cr (1.643) > As (0.172) > Cu (0.138) > Ni (0.137) > Zn (0.020) > Cd (0.010) > Pb (0.005), which indicated that the enrichment capacity of the edible part of rice for Cr in the environment was significantly higher than that for other heavy metals. The enrichment factor values of Zn, Cd, and Pb were all less than 0.1, indicating weak enrichment ability.
The transport coefficient is the ratio of the concentration of a heavy metal in the aboveground part of a crop to that in the roots. It can reflect the migration and transfer capacity of heavy metals from the roots to the aboveground part of plants [55]. The transport coefficients for heavy metals from root to edible part of rice were in the order As (0.316) > Cr (0.280) > Cu (0.159) > Ni (0.088) > Zn (0.071) > Pb (0.014) > Cd (0.013), which indicated that the transport capacity for As from root to edible part in rice was higher than that for other heavy metals.

3.3. Ecological Environment and Human Health Risk Assessment

3.3.1. Ecological Risk Assessment

The results of evaluation of the soil environment near the Pb–Zn mining area using the potential ecological risk and Nemeiro indices are shown in Figure 6. The single factor contamination index values of soil heavy metals followed the order Cd > Zn > Pb > As > Cu. Thus, the agricultural soil near the Pb–Zn mining area was severely contaminated with Cd and Zn, which presented mean single factor contamination index values of 91.216 and 15.346, respectively. Pb, As, and Cu presented single pollution index values > 1, with different degrees of pollution.
According to Nemeiro comprehensive pollution assessment, As, Cd, Pb, and Zn pollution was high in 10 sampling sites, and the Nemeiro index values followed the order Cd > Pb > Zn > As > Cu > Ni.> Cr. While Ni and Cu presented slight exceedances, other heavy metals greatly exceeded the standard values. The highest Nemeiro index was recorded in sampling point 1, which is located near the leachate flow path of the Pb–Zn mine, and the highest concentrations of Cd, Pb, and Zn were also noted in point 1, which may be the main reason for the high Nemeiro index at this site.
The potential ecological index of heavy metals followed the order Cd > Pb > As > Zn > Cu > Ni > Cr, with As and Zn posing medium risk, Pb posing high risk, and Cd posing very high risk. Therefore, the potential environmental impacts of As, Cd, Pb, and Zn should be carefully considered as they might be considerably high. Therefore, mitigation and management strategies should be developed and adopted to mitigate and prevent potential risks. Points 3, 4, and 10 were at moderate ecological risk, and the remaining seven points were at very high ecological risk, among which points 1 and 10 were at the highest and lowest risk, respectively.
The results of soil heavy metal pollution risk evaluation based on RAC and RSP are shown in Figure 7. The results of the RSP method (refer to Table 2) indicated that only Cd in soils near the tailings pond (point 3) was associated with heavy pollution, and the proportion of its biologically effective state was 83.19%, which represented a substantial threat to the surrounding environment. The other heavy metals were not associated with pollution levels, and Cr, Ni, Cu, and As were observed mainly in the stable residue state, which represents a lower potential threat to the environment. According to the RAC results, Zn posed a medium risk in 60% of the sampling points, namely in those near tailings, and low risk in the remaining points. Cd posed medium risk at point 3, and low risk in 80% of the points. Pb exerted no risk in 80% of the points. Ni posed low risk in 50% of the sampling points. Cr was only at low risk at point 10, the rest of the points were at no risk, whereas Cu and As presented no risk. The risk levels of the seven heavy metals in the farmlands were in the order Zn > Cd > Ni > Pb > Cu > Cr > As, which is consistent with the conclusion that Pb, As, Cu, Ni, and Cr have the highest residual state and poor migration ability. Together, Zn and Cd posed the greatest environmental threats.
The binding strength of different forms of heavy metals to the fixed medium varies, and different forms also present variability in chemical behavior and toxicity. Lower binding strength and weaker mobility lead to higher environmental risks [56]. The distribution characteristics of different forms of heavy metals in the environment can help estimate their toxicity, mobility, and potential ecological risks [57]. When total and morphological contents of heavy metals are used to assess environmental pollution and ecological risks, different specific values in the total evaluation method can lead to high variability of the results.
It is reported that an increase in the total heavy metal concentration was associated with an increase in the stable metal content (i.e., the residual state) [58]. However, the effect on active states, such as the exchangeable state, was minimal. Therefore, the exchangeable and carbonate-binding states should be carefully examined when evaluating ecological risk, as they have strong migration and toxicity. Evaluation of the states of heavy metals provides a more objective and accurate analysis of the potential ecological risk of pollutants to the environment.

3.3.2. Health Risk Assessment

Heavy metals in the soil can be absorbed by crop plants through the roots and accumulate in different organs and tissues. If the edible parts accumulate large amounts of heavy metals, they can cause serious human health effects via the food chain [59]. The risks of carcinogenic chemicals (Cd, As) and non-carcinogenic chemicals (Pb, Zn) in six crops to the health of adults and children are shown in Figure 8. The results show that the risks of As and Cd in the six crops were higher than the maximum acceptable risk level of 1 × 10−4 for both adults and children, which indicated that the accumulated As and Cd posed serious health risks. The intake of Pb through pak choi by children exceeded the maximum acceptable risk level and posed a risk of chronic toxicity. The risks of Pb to both adults and children were lower than 10 (except for the intake of radish by adults, which presented a low risk), but the possibility of negative effects was high. The health risks of Zn in pak choi and Chinese cabbage for children were above the maximum acceptable risk levels, and Zn in the other crops (below 10 for adults and children) presented a high possibility of health risks for children and adults.

4. Conclusions

The average pH of the investigated agricultural soils was 7.298 (4.58–8.06), which indicated weak alkalinity. The organic matter ranged within 2.822–7.572%, which indicated a low soil fertility. The concentrations of heavy metals in the Caiyuan River and leachate met the standard of surface water class V. The average pH of the water body was 7.587, and Eh and EC ranged from −69.4 to −8.6 mV and from 275 to 1487 μS·cm–1, respectively.
In farmlands, Cd exceeded the standard more substantially than the other metals, and Zn presented the highest abundance and accumulation and thus, the most evident pollution (maximum concentration of 2382.207 mg/kg). The average concentrations of As, Cd, Pb, and Zn were 2.927, 91.216, 9.144, and 15.346 times the soil environmental background values in Niujiaotang, and 1.805, 31.787, 5.980, and 17.390 times the soil environmental background values in Guizhou Province, respectively. The extraction steps of the chemical forms of heavy metal elements were conducted according to the sequential extraction scheme of the European Community Standards Agency (BCR sequential extraction procedure). The proportions of Zn states in agricultural fields were in the order residue > oxidizable > reducible > weak acid extractable state, whereas those of Cd states were in the order residue > weak acid extractable > reducible > oxidizable state.
The degree of heavy metal contamination of crops followed the order pak choi > Chinese cabbage > spring onion > rice > radish > Chrysanthemum coronarium, and the degree of individual heavy metal contamination of crops was in the order Zn > As > Cr > Cd > Ni > Pb. The highest content of Zn in the crops might be related to its higher soil concentration relative to other metals. Except for Cu, all heavy metals exceeded the limits of China’s food hygiene standards. Carcinogenic chemicals (Cd, As) and non-carcinogenic chemicals (Pb, Zn) in six crops posed significant health risks to adults and children, and Zn and Pb in pak choi and Zn in Chinese cabbage posed a risk of chronic poisoning to children. The enrichment capacity of the edible part of rice for Cr was significantly higher than that for the other heavy metals, and the transport capacity from the root system to the edible part of rice was higher for As than that for the other heavy metals.
Comparative analysis of the results of the four selected methods revealed that the farmlands were severely polluted. The Nemeiro index values were in the order Cd > Pb > Zn > As > Cu > Ni > Cr, and the potential ecological index values followed Cd > Pb > As > Zn > Cu > Ni > Cr. The RAC and RSP evaluation results showed that, together, farmland and crops had a high degree of pollution and posed a high ecological risk. The discrepant results obtained from evaluations based on total or morphological heavy metal contents indicated that the method based on morphological contents can more objectively and accurately evaluate the potential ecological risks of pollutants to the environment. Samples should be collected at different times to study the spatial and temporal distribution, and further studies on the migration transformation of heavy metals between the tailings pond-soil-crop should be conducted.

Author Contributions

All authors contributed to the study conception and design. The first draft of the manuscript, material analysis and the drawing of the diagram were completed by Q.L. All authors commented on previous versions of the manuscript. The first draft revision was completed by the corresponding author Z.H. Experimental operation was done by Q.L., Y.T. and M.Y. Data collation and sample collection were performed by Q.L. and H.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (2018YFC1801705, 2018YFC1801701), the National Natural Science Foundation of China (U1612442), and the Guizhou Talent Base Project (RCJD2018-21).

Data Availability Statement

The authors declare that data supporting the findings of this study are available within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shen, W.; Xiong, K.; Gao, Y.; Quan, M.; Peng, H.; Yang, T.; He, L.; Bao, K. Distribution of Potential Harmful Trace Elements and Potential Ecological Risk in the Jiulongchi Wetland of Fanjing Mountain, Southwest China. Int. J. Environ. Res. Public Heal. 2020, 17, 1731. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wright, C.A.; Moore, B.; Vervoort, J.; Beutel, M. Tracking long-distance atmospheric deposition of trace metal emissions from smelters in the upper Columbia River valley using Pb isotope analysis of lake sediments. Environ. Sci. Pollut. Res. Int. 2018, 25, 5501–5513. [Google Scholar]
  3. Zheng, N.; Wang, Q.; Zheng, D. Health risk of Hg, Pb, Cd, Zn, and Cu to the inhabitants around Huludao Zinc Plant in China via consumption of vegetables. Sci. Total. Environ. 2007, 383, 81–89. [Google Scholar] [CrossRef]
  4. Kwon, J.C.; Nejad, Z.D.; Jung, M.C. Arsenic and heavy metals in paddy soil and polished rice contaminated by mining activities in Korea. Catena 2017, 148, 92–100. [Google Scholar] [CrossRef]
  5. Wang, Y.; Wang, R.; Fan, L.; Chen, T.; Bai, Y.; Yu, Q.; Liu, Y. Assessment of multiple exposure to chemical elements and health risks among residents near Huodehong lead-zinc mining area in Yunnan, Southwest China. Chemosphere 2017, 174, 613–627. [Google Scholar] [CrossRef] [PubMed]
  6. Lin, Y.; Cook, N.J.; Liu, T.; Ciobanu, C.L.; Gao, W.; Yang, Y. The Niujiaotang Cd-rich zinc deposit, Duyun, Guizhou province, southwest China: Ore genesis and mechanisms of cadmium concentration. Miner. Depos. 2012, 47, 683–700. [Google Scholar]
  7. Zhang, J.; Wei, H.; Yang, R.; Gao, J.; Qu, Y. Study on Distribution Characteristics of Heavy Metals in Tailings from Niujiaotang Lead-zinc Mine Area, Duyun City, Guizhou Province. Nonferrous Met. Eng. 2018, 8, 122–127. [Google Scholar]
  8. Nassiri, O.; Rhoujjati, A.; EL Hachimi, M.L. Contamination, sources and environmental risk assessment of heavy metals in water, sediment and soil around an abandoned Pb mine site in North East Morocco. Environ. Earth Sci. 2021, 80, 1–20. [Google Scholar] [CrossRef]
  9. Candeias, C.; Ávila, P.F.; Da Silva, E.F.; Ferreira, A.; Durães, N.; Teixeira, J.P. Water–Rock Interaction and Geochemical Processes in Surface Waters Influenced by Tailings Impoundments: Impact and Threats to the Ecosystems and Human Health in Rural Communities (Panasqueira Mine, Central Portugal). Water Air Soil Pollut. 2015, 226, 23. [Google Scholar] [CrossRef]
  10. Jin, S.R.; Ye, L.; Yang, D.; Hu, Y.; Wei, C.; Dai, D.; Li, Z. A Discussion on the Pb-Zn Mineralization Regularity and Ore Prospecting Targets in the Kaili-Duyun Area, Guizhou Province, China. Acta Mineral. Sin. 2018, 38, 675–683. [Google Scholar]
  11. Tang, X.; Lu, H.; Cao, Z.; Xie, J. Morphological characteristics of homozygous wild rice phytoliths and their significance in the study of rice origins. Sci. China Earth Sci. 2021, 65, 107–117. [Google Scholar] [CrossRef]
  12. Peng, J.-F.; Song, Y.-H.; Yuan, P.; Cui, X.-Y.; Qiu, G.-L. The remediation of heavy metals contaminated sediment. J. Hazard. Mater. 2009, 161, 633–640. [Google Scholar] [CrossRef] [PubMed]
  13. Li, Y.; Arocena, J.M.; Zhang, Q.; Thring, R.W.; Li, J. Heavy metals and nutrients (carbon, nitrogen, and phosphorus) in sediments: Relationships to land uses, environmental risks, and management. Environ. Sci. Pollut. Res. 2017, 24, 7403–7412. [Google Scholar] [CrossRef] [PubMed]
  14. Tian, H.-J.; Feng, J.; Zhang, L.-M.; He, J.-Z.; Liu, Y.-R. Ecological drivers of methanotrophic communities in paddy soils around mercury mining areas. Sci. Total. Environ. 2020, 721, 137760. [Google Scholar] [CrossRef]
  15. Bin, H.E.; Yun, Z.J.; Shi, J.B.; Jiang, G.B. Research progress of heavy metal pollution in China: Sources, analytical methods, status, and toxicity. Sci. Bull. 2013, 58, 134–140. [Google Scholar]
  16. Weike, Z.; Liao, Y.; Zhou, S.; Zhou, B. Ecological remediation strategy for urban brownfield renewal in Sichuan Province, China: A health risk evaluation perspective. Sci. Rep. 2022, 12, 4300. [Google Scholar]
  17. Xiang, M.; Li, Y.; Yang, J.; Lei, K.; Li, Y.; Li, F.; Zheng, D.; Fang, X.; Cao, Y. Heavy metal contamination risk assessment and correlation analysis of heavy metal contents in soil and crops. Environ. Pollut. 2021, 278, 116911. [Google Scholar] [CrossRef]
  18. Bordean, D.-M.; Pirvulescu, L.; Poiana, M.-A.; Alexa, E.; Cozma, A.; Raba, D.; Borozan, A.; Misca, C.; Morar, A.; Obistioiu, D.; et al. An Innovative Approach to Assess the Ecotoxicological Risks of Soil Exposed to Solid Waste. Sustainability 2021, 13, 6141. [Google Scholar] [CrossRef]
  19. Salar, R.; Moghaddam, S.S.; Nouri, A.; Aqdam, K.K. Urbanization influences the distribution, enrichment, and ecological health risk of heavy metals in croplands. Sci. Rep. 2022, 12, 3868. [Google Scholar]
  20. Sajid, R.; AliShah, I.; Tulcan, R.X.S.; Rashid, W.; Sillanpaa, M. Contamination, exposure, and health risk assessment of Hg in Pakistan: A review. Environ. Pollut. 2022, 301, 118995. [Google Scholar]
  21. Mazumder, P.; Das, A.; Khwairakpam, M.; Kalamdhad, A.S. A comprehensive insight into ecological risk assessment and remediation of metal contaminated coal mine soil: Towards a cleaner and sustainable environment. J. Clean. Prod. 2021, 324, 129185. [Google Scholar] [CrossRef]
  22. Asare, E.A.; Assim, Z.; Wahi, R.; Fianko, J.R. Eco-toxic risk assessment and source distribution of trace metals in surface sediments of the coastal and in four rivers estuary of Sarawak. Beni-Suef Univ. J. Basic Appl. Sci. 2022, 11, 1–17. [Google Scholar] [CrossRef]
  23. Panghal, V.; Singh, A.; Kumar, R.; Kumari, G.; Kumar, P.; Kumar, S. Soil heavy metals contamination and ecological risk assessment in Rohtak urban area, Haryana (India). Environ. Earth Sci. 2021, 80, 1–20. [Google Scholar] [CrossRef]
  24. Xianfei, H.; Wu, X.; Tang, X.; Zhang, Z.; Ma, J.; Zhang, J.; Liu, H. Distribution characteristics and risk of heavy metals and microbial community composition around the Wanshan mercury mine in Southwest China. Ecotoxicol. Environ. Saf. 2021, 227, 112897. [Google Scholar]
  25. Sayma, S.; Khan, A.A.; Khan, M.R.; Ansari, S.A.; Khan, A. Assessment of Phytoremediation Potential of Seven Weed Plants Growing in Chromium- and Nickel-Contaminated Soil. Water Air Soil Pollut. 2021, 232, 209. [Google Scholar]
  26. Guo, B.; Hong, C.; Tong, W.; Xu, M.; Huang, C.; Yin, H.; Lin, Y.; Fu, Q. Health risk assessment of heavy metal pollution in a soil-rice system: A case study in the Jin-Qu Basin of China. Sci. Rep. 2020, 10, 11490. [Google Scholar] [CrossRef]
  27. Gao, J.; Wang, L. Ecological and human health risk assessments in the context of soil heavy metal pollution in a typical industrial area of Shanghai, China. Environ. Sci. Pollut. Res. 2018, 25, 27090–27105. [Google Scholar] [CrossRef]
  28. Deep, R.; Kumar, A.; Tripti; Risk, S.K.M.H. Assessment of Children Exposed to the Soil Containing Potentially Toxic Elements: A Case Study from Coal Mining Areas. Metals 2022, 12, 1795. [Google Scholar]
  29. Hou, S.; Tang, L.; Zheng, N.; Yang, Y.Y. Bioaccessibility and health risk assessment of heavy metals in vegetables of typical mining area. Acta Sci. Circumstantiae 2018, 38, 343–349. [Google Scholar]
  30. Yang, Q.-W.; Xu, Y.; Liu, S.-J.; He, J.-F.; Long, F.-Y. Concentration and potential health risk of heavy metals in market vegetables in Chongqing, China. Ecotoxicol. Environ. Saf. 2011, 74, 1664–1669. [Google Scholar] [CrossRef]
  31. Response, E.P.A.O. Superfund Public Health Evaluation Manual; EPA: Washington, DC, USA.
  32. Khezerlou, A.; Dehghan, P.; Moosavy, M.-H.; Kochakkhani, H. Assessment of Heavy Metal Contamination and the Probabilistic Risk via Salad Vegetable Consumption in Tabriz, Iran. Biol. Trace Element Res. 2020, 199, 2779–2787. [Google Scholar] [CrossRef] [PubMed]
  33. IRIS, U.E.O.N. Zinc and Compounds (CASRN 7440-66-6); | IRIS | US EPA: Washington, DC, USA, 2005. [Google Scholar]
  34. Environmental Protection Agency. Supplemental Guidance for Developing Soil Screening Levels for Superfund Sites; US EPA: Washington, DC, USA, 2002.
  35. Zhang, J.; Li, C.; Li, G.; He, Y.; Yang, J.; Zhang, J. Effects of biochar on heavy metal bioavailability and uptake by tobacco (Nicotiana tabacum) in two soils. Agric. Ecosyst. Environ. 2021, 317, 107453. [Google Scholar] [CrossRef]
  36. Ma, Y.; Rajkumar, M.; Zhang, C.; Freitas, H. Inoculation of Brassica oxyrrhina with plant growth promoting bacteria for the improvement of heavy metal phytoremediation under drought conditions. J. Hazard. Mater. 2016, 320, 36–44. [Google Scholar] [CrossRef] [PubMed]
  37. Lin, Y.E.; Li, C.Y.; Liu, T.G.; Li, X.B.; Pi, D.H. The environmental impact of Cadmium-rich Pb-Zn deposit for the example of Niujiaotang Cd rich Znic deposit, Duyun, Guizhou, China. Adv. Earth Sci. 2004, 67, 456–460. [Google Scholar]
  38. Pan, Z.P. The pollution characteristic of water during mining in Cadmium-rich Pb-Zn ore area as exemplifed by the NiujiaotangG Cd-rich znic deposit, Duyun, Guizhou Province. Acta Mineral. Sin. 2008, 28, 401–406. [Google Scholar]
  39. Hua, L.; Yang, X.; Liu, Y.; Tan, X.; Yang, Y. Spatial Distributions, Pollution Assessment, and Qualified Source Apportionment of Soil Heavy Metals in a Typical Mineral Mining City in China. Sustainability 2018, 10, 3115. [Google Scholar] [CrossRef] [Green Version]
  40. Huang, Z.; Liu, C.; Zhao, X.; Dong, J.; Zheng, B. Risk assessment of heavy metals in the surface sediment at the drinking water source of the Xiangjiang River in South China. Environ. Sci. Eur. 2020, 32, 1–9. [Google Scholar] [CrossRef] [Green Version]
  41. Qin, G.; Han, Z.; Tan, Y.; Xiao, H.; Yang, M. Chemical characteristics and quality evaluation of shallow groundwater in a typical karst plateau basin of Guizhou province: A case study of Houzhai river in Puding county. Environ. Chem. 2017, 36, 858–866. [Google Scholar]
  42. Wang, J.; Zhang, W.; Ren, E. Chemical state analysis of lead and zinc in soil of xitieshan mining area in the northern margin of qaidam basin. IOP Conf. Series: Earth Environ. Sci. 2018, 186, 012066. [Google Scholar] [CrossRef]
  43. Chlopecka, A.; Bacon, J.R.; Wilson, M.J.; Kay, J. Forms of Cadmium, Lead, and Zinc in Contaminated Soils from Southwest Poland. J. Environ. Qual. 1996, 25, 69–79. [Google Scholar] [CrossRef]
  44. Luo, G.; Han, Z.; Xiong, J.; He, Y.; Liao, J.; Wu, P. Heavy metal pollution and ecological risk assessment of tailings in the Qinglong Dachang antimony mine, China. Environ. Sci. Pollut. Res. 2021, 28, 33491–33504. [Google Scholar] [CrossRef]
  45. Haili, M.; Long, C.; Yang, Y.; Zhou, H.T. Speciation and Distribution Characteristics of Heavy Metals in Soil of Guizhou Niujiaotang Lead-Zinc Mining Area. Environ. Sci. Technol. 2014, 37 (Suppl. S2), 15–19. [Google Scholar]
  46. Kza, B.; Tappero, R.; Ruytinx, J.; Branco, S.; Liao, H.-L. Disentangling the role of ectomycorrhizal fungi in plant nutrient acquisition along a Zn gradient using in-situ X-ray imaging. Sci. Total Environ. 2021, 801, 149481. [Google Scholar]
  47. Wahsha, M.; Fontana, S.; Nadimi-Goki, M.; Bini, C. Potentially toxic elements in foodcrops (Triticum aestivum L., Zea mays L.) grown on contaminated soils. J. Geochem. Explor. 2014, 147, 189–199. [Google Scholar] [CrossRef]
  48. Hidayat, U.; Khaliq, M.; Ullah, N.; Iqbal, A.; Fozia; Ullah, I. Health Risk Assessment and Multivariate Statistical Analysis of Heavy Metals in Vegetables of Khyber Pakhtunkhwa Region, Pakistan. Biol. Trace Elem. Res. 2022, 200, 3023–3038. [Google Scholar]
  49. Mohsen, J.; Amir, M. Heavy metal contents, soil-to-plant transfer factors, and associated health risks in vegetables grown in western Iran. J. Food Compos. Anal. 2022, 106, 104316. [Google Scholar]
  50. Mengxue, W.; 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. [Google Scholar]
  51. Min, M.; Yang, L.; Wei, B.; Cao, Z.; Yu, J.; Liao, X. Plastic shed production systems: The migration of heavy metals from soil to vegetables and human health risk assessment. Ecotoxicol. Environ. Saf. 2021, 215, 112106. [Google Scholar]
  52. Deep, R.; Kumar, M.S. Risk assessment of potentially toxic elements in soils and vegetables around coal-fired thermal power plant: A case study of Dhanbad, India. Environ. Monit. Assess. 2020, 192, 699. [Google Scholar]
  53. Lu, M.; Liang, Y.; Lakshmanan, P.; Guan, X.; Liu, D.; Chen, X. Magnesium application reduced heavy metal-associated health risks and improved nutritional quality of field-grown Chinese cabbage. Environ. Pollut. 2021, 289, 117881. [Google Scholar] [CrossRef]
  54. Fang, F.M.; Wang, L.; Xie, H.; Wang, H. Enrichment characteristics and health risk assessment of heavy metals in vegetables in Sanshan District of Wuhu City of Anhui Province. J. Agro-Environ. Sci. 2010, 29, 1471–1476. [Google Scholar]
  55. Muhammad, M.S. Human health risk from Heavy metal via food crops consumption with wastewater irrigation practices in Pakistan. Chemosphere 2013, 93, 2230–2238. [Google Scholar]
  56. Pukalchik, M.A.; Terekhova, V.A.; Karpukhin, M.M.; Vavilova, V.M. Comparison of Eluate and Direct Soil Bioassay Methods of Soil Assessment in the Case of Contamination with Heavy Metals. Eurasian Soil Sci. 2019, 52, 464–470. [Google Scholar] [CrossRef]
  57. Sun, Y.-Q.; Xiao, K.; Wang, X.-D.; Lv, Z.-H.; Mao, M. Evaluating the distribution and potential ecological risks of heavy metal in coal gangue. Environ. Sci. Pollut. Res. 2020, 28, 18604–18615. [Google Scholar] [CrossRef]
  58. Ruiz-Huerta, E.A.; Armienta-Hernández, M.A.; Dubrovsky, J.G.; Gómez-Bernal, J.M. Bioaccumulation of heavy metals and As in maize (Zea mays L.) grown close to mine tailings strongly impacts plant development. Ecotoxicology 2022, 31, 447–467. [Google Scholar] [CrossRef]
  59. Yazdi, M.; Kolahi, M.; Kazemi, E.M.; Barnaby, A.G. Study of the contamination rate and change in growth features of lettuce (Lactuca sativa Linn.) in response to cadmium and a survey of its phytochelatin synthase gene. Ecotoxicol. Environ. Saf. 2019, 180, 295–308. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Map of the study area and sampling locations.
Figure 1. Map of the study area and sampling locations.
Toxics 11 00106 g001
Figure 2. Physical and chemical properties of soil samples. (a) Organic matter content of sampling points. (b) pH of sampling points. (sampling points: (1, 3, 5, 6) pak choi, (2) rice; (4) spring onion; (7, 10) radish; (8) Chinese cabbage; (9) Chrysanthemum coronarium).
Figure 2. Physical and chemical properties of soil samples. (a) Organic matter content of sampling points. (b) pH of sampling points. (sampling points: (1, 3, 5, 6) pak choi, (2) rice; (4) spring onion; (7, 10) radish; (8) Chinese cabbage; (9) Chrysanthemum coronarium).
Toxics 11 00106 g002
Figure 3. Physical and chemical properties of water samples: (a) pH; (b) redox potential Eh; (c) electrical conductivity EC.
Figure 3. Physical and chemical properties of water samples: (a) pH; (b) redox potential Eh; (c) electrical conductivity EC.
Toxics 11 00106 g003
Figure 4. Heavy metal contents in farmlands. Sampling points are the same as in Figure 2.
Figure 4. Heavy metal contents in farmlands. Sampling points are the same as in Figure 2.
Toxics 11 00106 g004
Figure 5. Chemical proportions of four characteristic heavy metals. Chemical species of: (a) Zn; (b) Cd; (c) As; (d) Pb. Sampling points are the same as in Figure 2.
Figure 5. Chemical proportions of four characteristic heavy metals. Chemical species of: (a) Zn; (b) Cd; (c) As; (d) Pb. Sampling points are the same as in Figure 2.
Toxics 11 00106 g005
Figure 6. Nemeiro comprehensive pollution assessment and potential ecological risk assessment of soils. (a) PN of heavy metals; and (b) RI of heavy metals.
Figure 6. Nemeiro comprehensive pollution assessment and potential ecological risk assessment of soils. (a) PN of heavy metals; and (b) RI of heavy metals.
Toxics 11 00106 g006
Figure 7. (a) Ratio of secondary phase to primary phase (RSP); and (b) risk assessment coding (RAC) evaluation results of soil samples.
Figure 7. (a) Ratio of secondary phase to primary phase (RSP); and (b) risk assessment coding (RAC) evaluation results of soil samples.
Toxics 11 00106 g007
Figure 8. Human health risk values of As, Cd, Pb and Zn in six crops for: (a) adults; and (b) children.
Figure 8. Human health risk values of As, Cd, Pb and Zn in six crops for: (a) adults; and (b) children.
Toxics 11 00106 g008
Table 1. Reference (Cin) and toxicity coefficient (Tir) for various heavy metals.
Table 1. Reference (Cin) and toxicity coefficient (Tir) for various heavy metals.
ElementCin (mg/kg)Tir
As12.3310
Cd0.2330
Cr66.62
Cu33.185
Ni29.395
Pb23.025
Zn112.751
Table 2. Evaluation criteria of the four methods for the evaluation of heavy metals in the soil.
Table 2. Evaluation criteria of the four methods for the evaluation of heavy metals in the soil.
MethodRisk LevelsLiterature
Potential ecological risk index (RI) methodRI < 150150 ≤ RI < 300300 ≤ RI < 600600 ≤ RI/[23]
LowModerateHighVery high/
Nemeiro index methodPN ≤ 0.70.7 < PN ≤ 11 < PN ≤ 22 < PN ≤ 3PN > 3[26]
CleanWarning LineMildModerateSevere
Risk assessment code (RAC)RAC < 1%RAC < 10%10% ≤ RAC < 30%30% ≤ RAC < 50%≥50%[22]
NoLowModerateHighVery high
Ratio of secondary phase to primary phase (RSP)RSP ≤ 11 < RSP ≤ 22 < RSP ≤ 3RSP > 3/[27]
NoMildModerateHigh/
Table 3. Parameters in the human health risk assessment model.
Table 3. Parameters in the human health risk assessment model.
ParameterAdultsChildrenSource
IR0.3550.233[29]
EF350350[30]
ED246
BW61.819.2
ATED × 365ED × 365[31]
RfDPb = 0.0035[32]
Zn = 0.3[33]
SFAs = 1.5; Cd = 6.1[34]
Table 4. Concentrations of heavy metals in water samples and respective national standard values for Class V surface water (mg/kg).
Table 4. Concentrations of heavy metals in water samples and respective national standard values for Class V surface water (mg/kg).
Heavy MetalsCrNiCuZnAsCdPb
River 14.5790.2330.33117.450.445N.D0.272
River 23.6590.0261.1986.8340.372N.D0.298
Leachate3.1464.732.21549.820.051N.D0.241
Paddy water4.0513.2921.996100.10.0530.2780.278
Ditch 13.9143.641.888253.30.0480.2080.290
Ditch 23.62.8331.839134.4N.D0.310.248
Ditch 33.83.3391.83442.080.95N.D0.328
Surface water standard (Class V)10020100020010010100
Table 5. Contents of heavy metals in the studied soils (mg/kg).
Table 5. Contents of heavy metals in the studied soils (mg/kg).
Type of ValueCrNiCuAsCdPbZn
Maximum69.9933.4051.3475.1558.75796.222382.21
Minimum32.8311.9322.5219.583.1573.68788.42
Average50.0024.9734.4036.1020.98210.501730.28
Median49.4725.7134.0635.9521.24144.912078.39
Standard deviation values9.106.857.7215.2015.63201.09630.02
Content in tailings1183.374328.82346.0330.6469.25200.935470.01
Background in Niujiaotang66.6029.3933.1812.330.2323.02112.75
Background in Guizhou Province95.9039.1032.0020.000.6635.2099.50
Table 6. Heavy metal concentrations in crops (mg/kg).
Table 6. Heavy metal concentrations in crops (mg/kg).
CropSampling PointCrNiCuAsCdPbZn
Pak choi12.9742.6144.06813.6099.2770.663323.529
32.0521.9114.86413.0025.1841.714779.324
517.9627.63912.31617.77815.06511.441618.290
662.52516.89913.10611.95811.2620.599305.666
Rice249.04221.7133.44314.2305.5611.982206.058
Spring onion42.6250.8705.16013.4352.7972.723216.351
Chinese cabbage82.9972.6634.69115.7365.6721.700317.910
Chrysanthemum coronarium92.8461.53112.99413.7933.8912.31162.066
Radish72.1421.5173.07614.7961.6310.594181.095
102.1371.6780.70612.0160.7910.25524.330
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Q.; Han, Z.; Tian, Y.; Xiao, H.; Yang, M. Risk Assessment of Heavy Metal in Farmlands and Crops Near Pb–Zn Mine Tailing Ponds in Niujiaotang, China. Toxics 2023, 11, 106. https://doi.org/10.3390/toxics11020106

AMA Style

Li Q, Han Z, Tian Y, Xiao H, Yang M. Risk Assessment of Heavy Metal in Farmlands and Crops Near Pb–Zn Mine Tailing Ponds in Niujiaotang, China. Toxics. 2023; 11(2):106. https://doi.org/10.3390/toxics11020106

Chicago/Turabian Style

Li, Qinyuan, Zhiwei Han, Yutong Tian, Han Xiao, and Miao Yang. 2023. "Risk Assessment of Heavy Metal in Farmlands and Crops Near Pb–Zn Mine Tailing Ponds in Niujiaotang, China" Toxics 11, no. 2: 106. https://doi.org/10.3390/toxics11020106

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