Abstract

The spatial patterns, ecological risks, and sources of heavy metals (HMs), including Pb, Zn, Mn, Cu, Cd, Hg, and As in river sediments, were identified around a lead-zinc mine of Danzhai, Guizhou, China. The concentrations of selected HMs and their coefficient variations indicated that the river sediments around this typical lead-zinc mine were obviously contaminated with HMs. Anthropogenic activities had further enhanced the accumulation of HMs. The higher contents of the most common selected HMs were mainly distributed in the area close to the lead-zinc mine. Based on the combined evaluations of the single factor pollution index, geo-accumulation index, and potential ecological risk index, it indicated that the ecological risks of Hg, Cd, Zn, and Pb were high or extremely high, and of Mn, Cu, and As were slight or none in the sediments around this lead-zinc mine. It was found that lead-zinc mining and smelting activities, coal mining activities, and agricultural activities (livestock and poultry breeding) are the primary sources of selected HMs, based on the results of correlation analysis together with principal component analysis (PCA) and positive matrix factorization (PMF) model. The pollution of HMs in the river sediments around a lead-zinc mine was predominantly caused by lead-zinc mining and smelting activities. Therefore, for environmental persistence, lead-zinc mining and smelting activities should be given careful consideration and under close surveillance.

1. Introduction

Heavy metals (HMs) are regarded as serious contaminations in water systems due to their environmental persistence, toxicity, nonbiodegradable properties, and accumulative behaviors [1, 2]. HMs pollution in water systems has caused wide concerns throughout the world, especially in developing countries where the economic activities are extensive [3]. Sediments are an important sink in aquatic ecosystems for HMs emitted from anthropogenic sources and will directly affect the overlying water and aquatic organisms [4, 5]. The concentrations of HMs in sediments can reflect the quality of the water column to some extent [6]. HMs can release into the overlying water and result in some secondary contamination problems when the environmental conditions change [7]. The spatial pattern of HMs can vividly present the concentrations of HMs (high or low) in the selected study area. The risk assessment of HMs would provide theoretical support for risk management. Numerous evaluation methods, such as geo-accumulation index (Igeo), potential ecological risk index (PERI), single factor pollution index (Pi), fuzzy coefficient method, gray clustering method, etc., have been applied to evaluate the contamination of HMs in the environment [812]. The selection of evaluation methods should consider their emphasis and the research purpose [13]. In general, multiple evaluation methods were conducted to obtain comprehensive results [5, 14].

The contents of HMs in the sediments around the areas with extensive economic activities mainly depend on anthropogenic activities such as rapid industrialization, mining, and smelting mineral resources, nonjudicial use of agrochemical, urbanization, agricultural production, etc. [15, 16]. These activities had resulted in that HMs were accumulated in soil and water systems and thus declined the environment quality over the past few decades in China [13, 1721]. In addition to anthropogenic sources, natural processes such as atmospheric sedimentation, rainfall, and rock weathering are important sources for HMs in the environment [22]. Therefore, it is necessary to identify the different sources of HMs in sediments [23], which could provide certain guidelines for the government to protect the aquatic ecosystems that are polluted with HMs [24].

There are multiple methods used to investigate the potential sources of HMs in soil and water environment, which include correlation analysis, principal component analysis (PCA), positive matrix factorization, geographic information systems, enrichment factor, etc. [23, 2529]. Although these methods can identify the common characteristics of selected HMs quickly via classification or dimension reduction, they cannot calculate the contribution rates of each heavy metal [3033]. In addition to these methods, positive matrix factorization (PMF) provides an ideal receptor that has been used to characterize the sources of contaminations in soil, sediments, and particulate matter [3436].

Danzhai, located in the southeast of Guizhou province, China, abounds with minerals such as lead, zinc, coal, Mercury, gold, etc. A long-term of more than 600 years of Hg mining activities have experienced in Danzhai, causing serious Hg pollution in soil [37, 38]. Moreover, lead, zinc, coal, and gold mining and smelting activities had developed rapidly in recent years. Large amounts of gangues and mine tailings are controlled stockpiled stacked optionally near mining sites. The acid mine drainage (AMD) generated from these solid wastes contains a high level of HMs and enters into the soil and aquatic environment under the impact of rainfall. Once the AMD enters the aquatic environment, the HMs will be accumulated in sediments, resulting in serious contamination of HMs in large mining areas [13, 39]. Therefore, we hypothesized that the sediments around mining and smelting sites are polluted with HMs seriously, with a certain spatial distribution pattern. Given the above consideration, the main objectives of this research were to (1) characterize the spatial patterns of HMs (Pb, Zn, Mn, Cu, Cd, Hg, and As) in the sediments around a typical lead-zinc mine, (2) evaluate the ecological risks of selected HMs to the aquatic ecosystem, and (3) accurately apportion the sources of HMs and identify their contributions in the study area for guiding restoration strategies using PMF model.

2. Materials and Methods

2.1. Site Description and Sample Collection

The lead-zinc mine is located southeast of Guizhou, China (Figure 1). Mining and smelting activities started in 2005. This region has complicated topography and a variety of landforms. The general altitudes range from 700 to 1100 above sea level. This area has a subtropical monsoon humid climate, with annual mean temperatures of 12.6–17.2°C with annual precipitations of 1259.7–1508.4 mm.

Surface sediments with a depth of 5 cm from 8 sites were collected in April 2019. The samples were transferred to the laboratory within 6 hours with preserving them in self-styled bags under 4°C. Site S1 is a livestock rearing area and is located in the upstream region of the lead-zinc mine. Sites S2 and S4 are just under the accumulation area of slag. Site S3 is the paddy field just under the accumulation area of slag. Sites S5–S8 are located in the river around this lead-zinc mine.

2.2. Chemical Analyses and Data Analysis

The preparation, digestion, and determination of sediment samples have been described in our previous work [5]. Simply, the sediment samples were dried freeze, dislodged the stones, debris, and organisms, and ground to 200-mesh. The ground samples were dispelled with HNO3, HF, and HClO4 under 25°C for 12 hours. Then these acids were volatilized at 160°C until solidified material became light yellow or white. The solidified material was dissolved with HNO3 and diluted with ultrapure water for determination [5]. Flame atomic absorption spectrometry (Perkin Elmer AAnalyst 800, USA) was employed for determining Pb, Cu, Zn, and Mn, graphite furnace atomic absorption spectrometry was used for measuring Cd, and cold atomic fluorescence spectrometry (Model III, Brooks Rand, USA) was employed for determining As and Hg, in prepared solution [5].

The recovery rates of Chinese national standard material (GSS-8) were determined to be 94 ± 2%, 90 ± 4%, 97 ± 3%, 103 ± 5%, 98 ± 4%, 112 ± 7%, and 106 ± 2% for Pb, Zn, Mn, Cu, Cd, Hg, and As, respectively. The relative standard deviation (%RSD) of all samples was less than 6%. The procedure of method blanks (stochastic 15% of samples) was conducted as detecting real samples, and contaminations were not found during analysis.

The single factor pollution index (Pi) [12, 40], geo-accumulation index (Igeo) [9], and potential ecological risk index (PERI) [8, 41, 42] were employed to evaluate the pollution status of HMs in the sediments around this lead-zinc mine. A detailed description of these three evaluation methods is provided in supplementary materials. The source of HMs in sediments around a lead-zinc Mine was apportioned by positive matrix factorization (PMF 5.0) [43], which was presented in supplementary materials.

3. Results and Discussions

3.1. Overview and Spatial Patterns of HMs

Table 1 displays the descriptive data of HMs in the sediments around this lead-zinc mine. The average concentrations of Pb, Zn, Cd, Mn, and Hg were 2.91, 4.32, 1.05, 3.79, and 18.6 times greater than their own corresponding background values, which indicated that the sediments were obviously polluted with these metals. Based on the values of coefficient variations for HMs, the degree of variation can be classified as high variation (>36%), moderate variation (16%–36%), and low variation (<16%) [44]. Although both the concentrations of Cu and As were lower than their own corresponding background levels, the coefficient variations of these two metals were high. The values of coefficient variations for Cd showed a moderate variation. The other metals showed a high variation and indicated that anthropogenic effects could result in the accumulating of HMs in sediments [45].

The concentrations of all HMs in sediments showed significant variations among the sampling sites (Figure 2). The comparatively high levels of Pb, Hg, Cu, and Zn in the sediments were generally exhibited in Sites S3, S4, and S5. The contents of Mn, S3, S4, and S5 had low values, indicating an opposite tendency of the distribution characteristics of Pb, Zn, Cu, and Hg. The peak concentrations of Cd can be found in sites S1 and S3–S7. The higher concentrations of As were mainly distributed in sites S2, S3, and S5. The contents of Zn, Cd, and Hg in sediments in all sampling sites were higher than their own corresponding background levels. The contents of Pb in sties S7 and S8, Mn in sites S3–S5, S5, S6, and S8, Cu in sites S1, S2, S4, and S6–S8, As in all sites, did not exceed their own corresponding background values. It is indicated that the higher levels of some HMs (Pb, Hg, Cu, and Zn) were mainly distributed in the area near this lead-zinc mine. These results also demonstrated that the accumulation degrees of Cd and Mn were complex and that there was nonpollution with As in our study area.

3.2. Ecological Risk

As shown in Figure 3(a), the Pi of Hg in sediments around this lead-zinc mine peaked in site S3 (57.62), followed by S5 (23.88), S4 (21.57), S8 (20.81), S1 (7.23), S2 (6.05), S7 (5.64), and S6 (5.63). All the values of Pi for Hg in every site were higher than 3, which indicated that the sediments around this lead-zinc mine were heavily polluted with Hg. The Pi of Zn in the sediments around the lead-zinc mine were higher than 3 in sites S3, S4, and S5, showing heavy pollution with Zn in these sites. The Pi of Zn in the sediments of sites S1, S7, and S8 were between 2 and 3, which indicated that these sediments were moderately polluted with Zn. Sites S2 and S6 were slightly polluted with Zn because the Pi of Zn in the two sites were between 1 and 2. The Pi of Cd in the sediments of all sites except site S2 (1.35) were calculated to be more than 3, indicating that the sediments around this lead-zinc mine were heavily polluted with Cd. The Pi of Pb in sediments of sites S3, S4, and S5 were higher than 3, sites S2, S6, and S7 were between 1 and 2, and sites S1 and S8 were lower than 1. It was found that sites S3, S4, and S5 showed heavy pollution, sites S2, S6, and S7 showed slight pollution, and sites S1 and S8 showed nonpollution with Pb. The sediments in most sites around lead-zinc mine were nonpollution or slight pollution caused by Cu and Mn because their Pi were lower than 1 or 2. The Pi of As in all sites were lower than 1, indicating that there were nonpollution caused by As. The mean Pi of Hg, Zn, Cd, Pb, Mn, Cu, and As in all sites were calculated to be 18.55, 4.32, 3.80, 2.91, 1.05, 0.99, and 0.36, respectively. This indicates that the sediments around the lead-zinc mine are heavily polluted with Hg, Zn, and Cd, moderately polluted with Pb, slightly polluted with Mn, and nonpollution with Cu and As as a whole.

By comparing the current concentrations of HMs with their corresponding background values, the Igeo can figure out the contamination status of HMs in environment [2]. Therefore, the Igeo were calculated for understanding the pollution level of HMs with respect to natural environment in the sediments around this lead-zinc mine. Although there are no background concentrations of HMs recorded in the sediments in Danzhai, the Igeo can be calculated as a function of background contents of soil elements in Guizhou Province, China. The overall value of Igeo ranged from −0.78 to 2.39 for Pb, −0.26 to 2.50 for Zn, −2.63 to 0.55 for Mn, −1.42 to 0.80 for Cu, −0.15 to 1.67 for Cd, 1.91 to 5.26 for Hg, and −3.04 to −1.18 for As. The rank order of mean values of Igeo for these HMs were Hg (mean 3.12) > Cd (mean 1.26) > Zn (mean 1.12) > Pb (mean 0.54) > Mn (mean −0.78) = Cu (mean −0.78) > As (mean −2.21) (Figure 3(b)). This suggested that the sediments around this lead-zinc mine were heavily polluted with Hg, moderately polluted with Cd and Zn, slightly polluted with Pb, and nonpollution with Mn, Cu, and As. The status of Hg contamination in the sediment of each site was ranked in the order of S3 (extreme pollution) > S4, S5, and S8 (heavy pollution) > S6 and S7 (moderate pollution). The surface sediments were moderate pollution of Cd throughout the study area, expect for site S2 (nonpollution). The degrees of Zn contamination were determined to be heavy pollution in sites S3, S4, and S5, slight pollution in sites S1, S6, S7, and S8, and nonpollution in site S2. The degree of Pb contamination was determined to be moderate to heavy pollution in sites S5, moderate pollution in S3 and S4, slight pollution in sites S6 and S7, and nonpollution in sites S1, S2, and S8. All sampling sites were determined to be nonpollution with Mn, Cu, and As except for S1 (slight pollution with Mn) and S5 (slight pollution with Cu). It was found that all sites were contaminated with Hg and Cd. In addition, the sediments near the lead-zinc mine were also contaminated with Pb and Zn.

According to the mean values of

of each element, the highest single risk index was Hg (742.15), followed by Cd (114.06), Pb (14.55), Cu (4.96), Zn (4.32), As (3.63), and Mn (1.05) (Figure 3(c)). This demonstrated that the ecological risks induced by Pb, As, Zn, Mn, and Cu were low for the sediments around lead-zinc mine. The values of

for Cd in all sites were between 80 and 160 except sites S2 (40.61 ± 4.0), indicating Cd caused considerable potential ecological risk for the sediments around lead-zinc mine. The values of

for Hg were calculated to be between 160 and 320 in sediments of sites S1, S2, S6, and S7, and more than 320 in the other sites. This demonstrated that Hg caused a high risk for half of the study area and extremely high risk for the other half. Besides, both Hg and Cd are significant polluted elements in the sediments around this lead-zinc mine.

The PERI is the sum of potential risk for selected HMs, which can figure out their combined toxicology, ecological, and environmental effects based on HMs contents [8]. PERI values in site S3 (2463 ± 48), S4 (1031 ± 41), S5 (1152 ± 88), and S8 (955 ± 68) were higher than 600, demonstrating that HMs in these areas had extremely high risk. For the rest of the sampling sites, PERI values were equal to or more than 300, which suggested that these sites posed considerable risk with HMs. Hg and Cd had the highest contributions to PERI (77.02% and 19.08%, respectively), which induced high PERI values in the sediments of this study area, especially in the sites near mining sites, e.g., sites S3 and 4. The contribution rate of the other individual heavy metal was calculated to be lower than 1, and their sum was 2.23%. It was found that the sediments in the study area were polluted at different degrees with various HMs based on the evaluation with Pi, Igeo, and PERI. Hg, Cd, Zn, and Pb should be given careful consideration among these HMs. There were slightly polluted or nonpolluted with As, Mn, and Cu in the sediments around this lead-zinc mine.

3.3. Correlation Analysis

If the HMs came from similar single or multiple sources, their contents in sediments would present positive relationships and vice versa. Relationships among selected HMs were thus analyzed by using Spearman’s correlation. It showed that there were positive relations among almost all HMs except Mn (Table S3). What's noted is that there were statistically significant correlations between As, Zn, Cu, and Pb. This showed that these HMs might come from the same sources related to anthropogenic activities, including mining and smelting, industrial production, and coal-fired power, etc. [31, 46, 47]. In our study area, the lead-zinc mine is the obvious pollution source for HMs. Therefore, we concluded that As, Zn, Pb, and Cu in the sediments were mainly from lead-zinc mining and smelting factory. Negative correlations were found between Mn and Pb, Zn, and Cd, indicating that lead-zinc mining and smelting factory is not the source of Mn for the sediments in the study area around the lead-zinc mine. The spatial distribution characteristics of HMs in sediments (Figure 2) confirmed this point. The contents of Mn were higher in the sites where the contents of As, Zn, Pb, and Cu were lower. Positive relationships were found between Hg and the other metals, especially Zn, with a significant coefficient of 0.717

. This indicated that lead-zinc mining and smelting factory was not the only source of Hg; some other sources such as large-scale Hg mining activities and atmospheric deposition of Hg might also affect the concentrations of Hg for the sediments in the study area [31, 48].

3.4. PCA for HMs in Sediments

PCA can describe the common features in contents and spatial distribution of selected HMs and based on these features to recognize their origins. The first three components contributed 92% to the variance and all their initial eigenvalues are higher than 0.5. The initial eigenvalues of the other four components are lower than 0.5 (Table S4), indicating that these components need not be discussed because they can explain little for the variability of HMs. As shown in Figure 4 and Table S5, principal components 1 (PC1) was predominantly controlled by Pb (0.87), Zn (0.63), Cu (0.96), and As (0.78) and these HMs accounted for 58.46% of the total variance. The study area is abundant in lead-zinc mineral resources and extensive mining such as activated slag stacking optionally for the long term, resulting in the sediments contaminated with Pb and Zn. The association phenomena of As, Zn, Pb, and Cu deposits were observed in the study area, where belongs to the low-temperature mineralization area. As they originated from the low-temperature, hydrothermal fluids can facilitate the epithermal deposits of these metals [49]. Therefore, positive relationships were found among these HMs, and they would exhibit common characters such as originating from the same source, enrichment, and transport. We thus concluded that PC1 was mainly from mining and smelting activities of lead-zinc mine.

Principle components 2 (PC2) was dominated by Hg (0.89) and Zn (0.69) and explained 18.61% of the total variance (Table S5). Guizhou province is located in the center of the circum-Pacific, where the content of Hg of natural soils is generally high because of high background values of Hg [48]. More importantly, Guizhou is a significant site for producing Hg in the world. Although all the large-scale mining activities for producing Hg had been ceased in the study area, the environment had been seriously contaminated with Hg already [48]. In addition, Guizhou province is rich in zinc ores, which contain a high level of Hg. The smelting activities of Zn can discharge a large amount of Hg into the environment [48, 50]. Therefore, it was found that PC2 mainly came from soil surface runoffs, Hg mining activities, and smelting of Zn.

Principle components 3 (PC3) can explain 14.99% of the total variance and was dominated by Cd (0.95). Generally, Cd is the indicator element for agricultural activities, especially for livestock and poultry breeding [51]. Besides, activities of mining and smelting in the lead-zinc mine can discharge a large amount of Cd into the environment [52]. The spatial distribution characteristics of Cd showed that higher values of Cd concentrations were observed mainly in the region near the livestock rearing area and close to lead-zinc mine, which supported these viewpoints. Therefore, the origins of Cd is mostly related to livestock and poultry breeding and lead-zinc mine.

3.5. Source Apportionment by PMF

Except for correlation analyses and PCA, the PMF model was run to identify where the HMs come from to the sediment more accurately. The minimum Q true value was found for regulating the residual matrix E and determining the reasonable factor numbers to obtain a rational result of PCA. S/N ratios categorized all the HMs to be “strong.” The factor number was set to 3, 4, and 5, and the PMF model run 20 times. The Q true value was smallest (3.62) if the factor number was set to be 3. The fitting coefficients between the determined contents and the model predicted concentration were shown in Table S6. Although the fitting coefficient of Hg was relatively low, and of the other HMs was above 0.70, which suggested that the fitting effects of PMF could elucidate the potential information of original data and then met the research needs.

Factor 1 was characterized by Mn and As (Figure 5(a)), and their corresponding contributions were calculated to be 64.2% and 45.7% (Table 2). Compared with the background level of Mn, the concentrations of Mn in the sediments around this lead-zinc mine were higher, and their coefficient of variations was 59% (Table 1). This demonstrated that anthropogenic activities had significant effects on the accumulation of Mn in the sediments around this lead-zinc mine. Mn contamination was correlated with the coal mining activities, which is consistent with the fact that Guizhou province is rich in coal sources [53]. All the contents of As in the study area were lower than its own background value, and the coefficient of variation of As was lower than that of the other HMs (Table 1), indicating that As originated from the parental materials of soil. Factor 1 was thus identified as an industrial source associated with coal mining.

Factor 2 was predominantly characterized by Pb, Hg, Zn, and As (Figure 5(b)), and their contributions were calculated to be 75.2%, 73.8%, 73.3%, and 51.4%, respectively (Table 2). We could ignore the contribution of As because As originated from the parental materials of soil. The average contents of Pb, Hg, and Zn were greater than their own corresponding geological baseline, and the coefficient of variation of these HMs was higher than 70% (Table 1), indicating external sources were responsible for these HMs. In the study area rich in mineral resources with Pb, Zn, and Hg [48], extensive mining activities could result in the accumulation of these HMs in the sediments around this lead-zinc mine. Therefore, factor 2 was related to mining and smelting activities.

Factor 3 was mainly characterized by Cd (Figure 5(c)), its contribution was 62.7% (Table 2). The contents of Cd in the sediments around the lead-zinc mine were greater than its own background value, and their coefficient of variations was calculated to be 28% (Table 1). This suggested that the contamination of Cd in the sediments around this lead-zinc mine was caused by human activities. It has been identified that mining and smelting activities of lead-zinc mine could strongly raise the contents of Cd in the environment [5457]. In addition, Cd could be recognized as indicator elements for livestock and poultry breeding [51]. So it is was found that factor 3 was mainly induced from industrial and agricultural sources.

The results obtained from the PMF model showed that factor 1 (coal mining activities), factor 2 (lead-zinc mining and smelting activities), and factor 3 (lead-zinc mine and agricultural activities), respectively, contributed 24.1%, 48.8%, and 27.1% to the total sources of HMs in the sediments around this lead-zinc mine (Table 2). Compared with agricultural sources, the contamination of HMs in the sediments of the survey region was predominantly caused by lead-zinc mining and smelting activities.

4. Conclusions

Except for Cu and As, the concentrations of selected HMs were greater than their own corresponding background levels. The concentrations of HMs together with their coefficient variations indicated that the river sediments around a typical lead-zinc mine were obviously contaminated with HMs and anthropogenic activities had significant effects on the accumulating of HMs. The concentrations of all HMs in the sediments showed significant variations among the sampling sites, and the higher contents of almost all HMs were mainly found in the sites near the lead-zinc mine.

The sediments around this lead-zinc mine were polluted to different degrees with various HMs based on different evaluation methods. Based on Pi, the sediments around this lead-zinc mine were heavily polluted with Hg, Zn, and Cd, moderately polluted with Pb, slightly polluted with Mn, and nonpolluted with As and Cu. The Igeo indicated that the sediments were determined as heavily polluted with Hg, moderately polluted with Zn and Cd, slightly polluted with Pb, and nonpolluted with Mn, Cu, and As. The PERI found that the ecological risks of HMs were predominantly caused by Hg and Cd. In sum, Hg, Cd, Zn, and Pb should be given careful consideration, and there was slight pollution or nonpollution with Mn, As, and Cu in the sediments around this lead-zinc mine.

Correlation analyses suggested that lead-zinc mining and smelting factory was the major source of Pb, As, Zn, Cu, and Hg in the sediments around this lead-zinc mine, and some other sources would also affect the contents of Hg. PCA suggested that Pb, Cu, As, and Zn predominantly came from mining and smelting activities of lead-zinc mine. Hg and Zn mainly came from soil surface runoff, Hg mining activities, and refining of Zn. Cd was mostly associated with livestock and poultry breeding and lead-zinc mine. The PMF model found that Mn originated from coal mining activities, Pb, Zn, and Hg were from mining and smelting activities of lead-zinc, and Cd originated from lead-zinc mine and agricultural activities. In summary, the HMs contamination in the sediments around this lead-zinc mine was predominantly caused by lead-zinc mining and smelting activities.

Data Availability

The data on heavy metal contents of sampling sites, results of the evaluation, correlation analysis, PCA, and PMF are included within this article and the supplementary material.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was financially supported by Science and Technology Department of Guizhou Province (No. Qian ke he ji chu [2020]1Z038) and the Doctoral Scientific Research Foundation of Guizhou Normal University for 2014 (No. 11904/0514154).

Supplementary Materials

Supplementary Materials: (a) Evaluation methods of contamination status of HMs. (b) Description of source apportionment by PMF. (c) Tables S1 and S2: the classification and indications of Igeo and PERI. (d) Table S3: the results of correlation analysis. (e) Tables S4 and S5: the results of PCA. (f) Table S6: the results of PMF. (Supplementary Materials)