Abstract

Various studies have shown that soils surrounding mining areas are seriously polluted by heavy metals. In this study, 58 topsoil samples were systematically collected throughout the coal mining city Wuhai, located within the Inner Mongolia Autonomous Region of China. The concentrations of As, Hg, Cr, Ni, Cu, Zn, Cd, and Pb in these samples were measured and statistically analyzed. The mean concentrations of all heavy metals were lower than their Grade I values defined by the Chinese Soil Quality Standard. However, the mean concentrations of individual heavy metals in many samples exceeded their background values. The spatial distribution of heavy metals was analyzed by the ordinary kriging interpolation method. The positive matrix factorization model was used to ascertain contamination sources of the eight heavy metals and to apportion the contribution of each source. The most severely polluted area was the Wuhushan mine site in the Wuda district of Wuhai. Our results showed that coal mining strongly affected heavy metal contamination of the local soils. Results of source apportionment indicated that contributions from industrial activities, atmospheric deposition, agricultural activities, and natural sources were 31.3%, 26.3%, 21.9%, and 20.5%, respectively. This clearly demonstrates that anthropogenic activities have markedly higher contribution rates than natural sources to heavy metal pollution in soils in this area.

1. Introduction

Heavy metal contamination has attracted increasing concern over the past decades because it represents a potential hazard to human health and to the environment [13]. In particular, heavy metal pollution in soils has become a serious environmental problem around the world, linked to rapid urbanization and industrialization [4, 5]. Heavy metals in soils are completely undegradable and are easily accumulated in organisms through the food chain [6, 7]. In addition, they can be transported to distant downwind areas through atmospheric diffusion and to rivers through surface runoff under rainfall condition [8]. Therefore, soil pollution by heavy metals has become a serious threat to human health and regional ecological systems.

It is commonly recognized that natural and anthropogenic activities are the two major origins of heavy metals. Natural sources are mainly related to trace elements in the crust [9], while anthropogenic sources are mainly related to metal smelting, fuel and coal combustion, coal mining, and agricultural fertilization [10, 11]. To identify sources and understand the spatial variation of heavy metals in soils, multivariate statistics combined with geostatistical methods were generally used [1214]. Typically, geostatistical methods are applied to study the spatial distribution characteristics of heavy metals in soils. The ordinary kriging method is the main interpolation method frequently used to study the spatial distribution characteristics of heavy metals in soils [12, 15]; it is based on measured data and a semivariogram model to predict unknown sites [12, 15, 16]. The US Environmental Protection Agency’s positive matrix factorization (EPA’s PMF) model is widely applied to evaluate the sources of heavy metals in soils. The PMF model analyzes and quantifies the contributions of different pollution sources to samples. There have been numerous studies that systematically analyzed the pollution status of heavy metals in soils by combining geostatistical methods with the EPA’s PMF model [14, 17, 18].

Coal mining, grinding, coal transportation, and coal combustion processes are known to transfer several heavy metals to the topsoil, leading to contamination of heavy metals in mining regions [1921]. Tang et al. [18] reported concentration levels of heavy metals in soils in Fuxin, China. The mean concentrations of As, Cd, Cr, Ni, and Zn all exceeded the soil background values of these heavy metals. In fact, the mean concentration of Cd was nine times higher than its background value. In contrast, the mean concentrations of Cu, Hg, and Pb were lower than their respective background values. Sun et al. [22] studied the concentrations of heavy metals in soils of Tangshan, China. The mean concentrations of Zn, Pb, Cd, Hg, and Cu were higher than their background values, while those of Cr, Ni, and As were lower than their background values. None of them exceeded their corresponding Grade II values defined by the Chinese Soil Environmental Quality Standard. Hua et al. [23] measured the concentrations of heavy metals in soils of Daye, China. The mean concentrations of Cd, Cu, Pb, and Zn were higher than their background values, while Cr, Ni, and Mn concentrations were not. The mean concentration of Cd was seven times higher than its background value. Liang et al. [14] investigated average contents of twelve heavy metals (Zn, Cd, Cu, Hg, Pb, Sb, As, Mo, V, Fe, and Cr); except for Mn, all exceeded their background levels in soils in Hunan Province. In fact, in Lianyuan, the Hg level was higher than its corresponding Grade II value defined by the Chinese Soil Environmental Quality Standard. Reza et al. [16] researched the average contents of heavy metals in a mining area of northwest India. The mean concentrations of Cr, Cd, Ni, and Pb were 112.3 mg/kg, 2.60 mg/kg, 87.5 mg/kg, and 183.1 mg/kg, respectively. Overall, the degree of soil polluted by heavy metals varies significantly in different coal mining cities. Furthermore, the sources contribution was dominated by anthropogenic activities. In contrast to other land-use areas [17, 24, 25], studies of heavy metal pollution in soils of coal mining areas are relatively few at present. Therefore, the existing information on heavy metals in soils of coal mining areas is clearly inadequate for researchers and managers to comprehensively evaluate their contamination levels and potential risks.

Wuhai is a typical coal mining city within the Inner Mongolia Autonomous Region, China. There are 3.147 billion tons of coal reserves in Wuhai, of which more than 0.5 billion tons of coal have been excavated. The development of chemical, coal, iron, and steel production and other industries has likely led to the deterioration of soil quality in Wuhai [26]. However, the pollution status of heavy metals in soils in Wuhai has not been investigated to date.

The main objectives of this study were (1) to report contamination levels of eight heavy metals, namely, As, Hg, Cr, Ni, Cu, Zn, Cd, and Pb; (2) to ascertain spatial distribution characteristics using geostatistics and the ordinary kriging interpolation method; and (3) to quantitatively identify sources of heavy metals in soils by using the EPA’s PMF model.

2. Materials and Methods

2.1. Study Area

Wuhai (39°02′30″–39°54′5″N; 106°36′25″–107°08′05″E) is located in the southwest of the Inner Mongolia Autonomous Region, adjacent to the Ningxia Hui Autonomous Region, and at the junction of three deserts—the Wulanbuhe, the Kubuqi, and the Maowusu. The Yellow River crosses the region from south to north along its eastern edge. The city covers an area of around 1754 km2, encompassing three administrative districts—Haibowan, Wuda, and Hainan. The city terrain is the highest on its east and west sides and the lowest at its center. Average annual precipitation and evaporation are 160 and 3500 mm, respectively. The diurnal temperature (−36.6°C–40.2°C) varies greatly, which is characteristic of its typical continental climate. Wuhai has abundant coal resources, comprising mainly bituminous and coking coals. Its industries are centered on coal mining, chemistry, electric power generation, and building materials. As an emerging energy city, its heavy metal soil contamination is closely related to its rapid industrial development.

2.2. Sample Collection and Concentration Determination

A total of 58 surface soil samples were collected from Wuhai during 25 to 28 July, 2019 (Figure 1). Sampling sites were selected using a grid, with a density of about one sample every 3–5 km [27, 28]. To ensure that the sampling sites represented the soil quality of the entire monitoring area, the land-use status and geological information were documented to establish context. The sampling sites covered the coal mining area, both sides of the Yellow River, near the main highway, beside a factory, as well as farmland, grassland, and desert areas. All of the soil samples were collected from the 0–20 cm topsoil horizon, using a stainless spade [29]. They were immediately sealed in polyethylene plastic bags to avoid contamination and labeled. The latitude and longitude of each sampling point were recorded using a GPS, and any relevant environmental information pertaining to the sample site was recorded.

Sand (0.05–2 mm) was the dominant soil particles in the 0–20 cm soil profiles in the study area. Thus, all samples were dried in open air at room temperature in a laboratory and sieved through a 2 mm sieve, after removing gravel and plant materials [30]. The processed samples were sealed in polyethylene plastic bags for concentration determination. The determination of Cr, Ni, Cu, Zn, Cd, and Pb was carried out according to the method outlined in HJ803-2016; Hg was determined according to the method outlined in DZ/T0279.17-2016; and Arsenic (As) was determined according to the method outlined in DZ/T0279.13-2016. The concentrations of Cr, Cd, Pb, Cu, Zn and Ni were analyzed by inductively coupled plasma mass spectrometry; while Hg and As were measured using atomic fluorescence spectrometry. The accuracy of the elemental analyzes was verified using standard reference materials: GSD-31, GSS-32, GSS-33 and GSS-34. The results showed that the recoveries of these eight heavy metals were within the range 94.6%–106.8%. Each set of ten samples included one duplicate sample, and relative standard deviations were lower than 15% for all batch treatments. Method detection limits (MDLs) were 1, 0.0005, 5, 1, 1, 4, 0.02, and 0.5 mg/kg for As, Hg, Cr, Ni, Cu, Zn, Cd, and Pb, respectively.

2.3. Data Analysis
2.3.1. Descriptive Statistical Analysis

SPSS 24.0 software (IBM Corp., Armonk, NY, USA) was used for basic descriptive statistical analysis of the sample data. This provided the maximum, minimum, mean, standard deviation, coefficient of variation, skewness, and kurtosis and allowed us to assess the normality for the sampled heavy metal concentrations. The correlations between heavy metals were assessed using Pearson’s correlation method.

2.3.2. Geostatistical Analysis

Geostatistical analysis is based on regional variable theory and is used to study distribution characteristics and variation with geographic properties. Meanwhile, a semivariogram model was used to interpolate unknown points, as an important complementary tool [31]. In this study, ArcGIS 10.2 (Esri, Inc., Redlands, CA, USA) was used for calculating the spatial distribution characteristics of soil heavy metals using the ordinary kriging interpolation method [32]. The original concentration data, which did not meet the conditions of a normal distribution, were normalized using a logarithm transformation, and a normality test was carried out using a Q-Q plot. The semivariogram model and kriging calculation method are defined as follows [33]:where is semivariogram, is the step length or distance, is the number of pairs of sample points separated by . is the measured value of heavy metals in the soil at the regionalized variable position , is the concentration of the heavy metal to be tested, is the number of sampling points, and is the set of weight coefficients. The semivariogram model includes linear, spherical, exponential, and Gaussian components. In ordinary kriging space interpolation, the optimal fitting model improves the accuracy of interpolation.

2.4. Positive Matrix Factorization Model

The EPA’s PMF 5.0 is a modified factor analysis receptor model, proposed by Paatero [34]. It has been successfully applied to determine the source apportionment of various types of environmental pollution. In our study, EPA’s PMF 5.0 was used to determine the source apportionment of heavy metals in local soils. According to the EPA’s PMF 5.0 user guide,where is the measurement matrix of the heavy metal element in number of samples; is a contribution matrix of the source factor for number of samples; is a source profile of heavy metal element for the source factor; and refers to the residual value for the metal element in number of samples. The minimum value of the objective function Q can be computed by the following formula:where refers to the uncertainty in the heavy metal element for sample . The remarkable feature of PMF is that it uses uncertainty to analyze the quality of all individual concentration data points. If the concentration of heavy metal does not exceed the MDL value, then the uncertainty is calculated using the following formula:while if the concentration of heavy metal exceeds its corresponding MDL value, then the uncertainty is calculated using this formula:

3. Results and Discussion

3.1. Levels of Heavy Metal Contamination

The statistical characteristics of concentrations of As, Hg, Cr, Ni, Cu, Zn, Cd, and Pb in all soil samples are presented in Table 1. In this study, the soil background values for the Inner Mongolia Autonomous Region were used as evaluation criteria [35]. In addition, the concentrations of heavy metals in soils were compared with Grade I and Grade II values of the Chinese Soil Environmental Quality Standard (GB15618-1995). The mean concentrations of As, Hg, Cr, Ni, Cu, Zn, Cd, and Pb in soils were 9.31 mg/kg, 0.12 mg/kg, 61.7 mg/kg, 24.7 mg/kg, 19.6 mg/kg, 55.2 mg/kg, 0.16 mg/kg, and 28.0 mg/kg, respectively. Generally, the mean concentrations of heavy metals in the soils were all lower than their Grade I values defined in the Chinese Soil Environmental Quality Standard, although they exceeded their background values. The exceedance ratios of As, Hg, Cr, Ni, Cu, Zn, Cd, and Pb were 100%, 29.3%, 100%, 98.3%, 96.6%, 39.6%, 100%, and 100%, respectively. The most serious heavy metal contaminations were associated with Hg and Cd; these species had mean concentrations that exceeded their background values by 4.3 and 4.2 times, respectively. This is consistent with previously published statements that Hg and Cd produce the most serious heavy metal contaminations in China [14, 22, 36].

The coefficients of variation of As, Hg, Cr, Ni, Cu, Zn, Cd, and Pb were 24.7%, 540%, 15.9%, 22.6%, 71.5%, 81.2%, 48.3%, and 139%, respectively (Table 1). The variability of heavy metal concentrations was ranked: Hg > Pb > Zn > Cu > Cd > As > Ni > Cr. Therefore, Hg, Pb, Zn, Cu, and Cd had higher coefficients of variation and wider concentration ranges than other heavy metals. Soil samples revealed very heterogeneous spatial distributions of these metals, indicating that their concentrations were strongly influenced by anthropogenic sources [37]. In contrast, As, Cr, and Ni had relatively low coefficients of variation, suggesting that these heavy metals were less affected by anthropogenic sources.

Table S1 shows a comparison of recently published heavy metal concentrations in soils from typical coal mining cities and nonmining cities throughout China. When compared to data for other typical coal mining cities, our results indicate that soil concentrations of As, Cr, Ni, Cu, and Pb represent moderate contamination levels; the soil concentration of Hg represents a high contamination level, while soil concentrations of Zn and Cd represent low contamination levels. When compared to other studies conducted in nonmining cities, the concentrations of As, Hg, Cr, Ni, and Pb represent moderate contamination levels, while Cu, Zn, and Cd represent low contamination levels. Additionally, a comparison of heavy metal concentrations in soils was also made between studies conducted in typical coal mining cities and nonmining cities. The results show that concentrations of Hg in mining cities were slightly lower than those in nonmining cities, but concentrations of As, Cr, Ni, Cu, Zn, Cd, and Pb were not significantly different among cities.

The above discussion indicates that heavy metals in soils of Wuhai have been contaminated to different degrees, with Hg and Cd pollution being more severe than for other heavy metals. Hg, Pb, Zn, Cu, and Cd showed more spatial variation than other heavy metals. According to our comparison of concentrations of heavy metals in Chinese coal mining cities and nonmining cities, heavy metal concentrations in soils of Wuhai represent moderate contamination. Therefore, it is necessary to analyze the spatial distribution characteristics to determine the heavy metal sources within this region.

3.2. Spatial Distributions of Heavy Metal Concentrations

The spatial distributions of heavy metal concentrations in the soils of Wuhai are shown in Figure 2. Generally, the spatial trends of heavy metals within the study area are as follows: the degree of heavy metal pollution is the highest in the middle to western parts of Wuhai, with lower degrees of pollution to the south and north. In the midwestern part of Wuhai, there are mainly coalfields and a major highway. In the south of Wuhai, there are various infrastructure and land-use types, including a gas company, a cement plant, a petrol station, and various expressways, as well as large areas of farmland and grassland. In the north of Wuhai, there are expressways and plastic manufacturing and coal chemical engineering plants. Heavy pollution of Hg, Ni, Cu, Zn, Cd, and Pb was found in samples that were located in the midwestern part, proximal to the Wuhushan coalfield in the Wuda district of Wuhai. In this region, there are coal mining and coal transportation activities. Clearly, coal mining activities made a significant contribution to heavy metal contamination of soils in the midwestern part of Wuhai.

Spatial distributions of Cu and Zn in soils were highly similar, indicating that Cu and Zn may originate from the same source. Large areas having soils contaminated with Cu and Zn were situated in the midwestern, northern, and southwestern parts of Wuhai. In the midwestern area, high concentrations of Cu and Zn were proximal to the coalfield and the highway, close to coal chemical industries and plastic manufacturing enterprises in the north, and close to a gas company, a cement plant, a gas station and expressways in the south. We speculate that Cu and Zn contamination of the soil were caused by industrial emissions, vehicle exhaust, and mechanical deterioration particles, which were later dispersed and deposited by atmospheric processes [13, 38, 39].

The Pb hotspots were located in midwestern and southwestern parts of Wuhai, as well as within the central Haibowan district of Wuhai. These areas were proximal to the major highway and surrounding coalfields. The highest concentration area was associated with the highway. Vehicle exhaust and coal dust were likely important sources of Pb [40, 41]. This indicates that Cu, Zn, and Pb inputs to soils mainly occurred via atmospheric deposition.

The hotspots of Cd in soils were mainly distributed in the midwestern and southwestern parts of Wuhai. In the midwestern area, Cd was associated with the highway and surrounding coalfield. Therefore, coal mining and its associated activities are likely the primary contributors to Cd contamination [42]. Elsewhere, the Cd distributions coincided with Pb distributions, especially in the southwest. This suggests a similar dispersion of Cd and Pb, related to a similar source. Vehicle exhaust could contribute both Cd and Pb inputs to soils [13, 43]. However, other large-scale areas with moderate concentrations of Cd were present in southern to central areas of Wuhai. These areas host numerous industrial activities, linked to a gas company, a gas station, and a cement plant. These industries have been previously linked to Cd contamination in soils [44]. Cd and Zn had zonal distributions toward the southwest, along the Yellow River, possibly related to fluvial processes. Pavolvic et al. reported that both Cd and Zn have high potential for adhering to clay minerals and are readily transported by rivers. This may lead to their accumulation in riparian soils [45].

The hotspots of Hg were in the midwestern part of Wuhai. In this region, Hg had an extreme value of more than 60 times the background value, far exceeding its Grade II value. Given the long history of coal mining in this area, it is likely that large spoil quantities have been generated, which release Hg-rich effluents into the surrounding environment [22]. Moreover, other moderately polluted areas were scattered in central and southern parts of Wuhai, associated with the gas company, the gas station, and the expressways. Previous studies have indicated that petrochemical plants are a major source of Hg [46]. Clearly, the Hg distribution was predominantly influenced by coal mining and the local petrochemical industry [42, 47].

Arsenic (As) concentrations increased continually from northeast to southwest, with hotspots of Arsenic (As) situated on the southwest limit of Wuhai within farmland and a village. A number of reports have identified Arsenic (As) as originating from agricultural pesticides and fertilizers [48, 49]. In this region, Arsenic (As) contamination in soils was also likely related to agricultural activities.

High concentrations of Cr were documented in the central Hainan district of Wuhai, related to both the coalfield and the local cement plant. A Ni hotspot, located in the midwestern part of Wuhai, was also associated with coal mining. In addition to the hotspot in the midwestern area, another Ni hotspot occurred within the southern to central desert area, having no obvious point sources of Ni. The coefficients of variation and spatial variances of Cr and Ni were both relatively low. This indicates that Cr and Ni levels in soils were primarily controlled by parent materials [50].

Overall, the spatial distributions of Cu, Zn, and Pb were mainly affected by atmospheric deposition of airborne particles; Hg was mainly related to coal mining; Cd was mainly related to industrial activities; Arsenic (As) was mainly related to agricultural activities; and Cr and Ni were mainly related to parent materials of the soils. Furthermore, spatial distribution of soil heavy metals is rarely influenced by transport in the atmosphere since sand has the highest fractions in the research area. The sources of these heavy metals are discussed in detail below.

3.3. Source Apportionment of Heavy Metals

Pearson’s correlation coefficients of the eight heavy metals are shown in Table S2. Strong correlations were found between Pb and Cu, Zn, Cd (), while pairs Ni and Cr, Cu and As were slightly correlated (). Hg did not display a correlation with any other heavy metals. Significant correlations between metal elements in soils suggest that they had similar contamination sources. For example, Cu, Zn, and Pb may be derived from the same source (r > 0.9), as reported by Ma et al. [13]. Given the complexity of heavy metal pollution in soils related to human activities, the conclusions based on our Pearson’s correlation coefficient analysis are not comprehensive. Therefore, the EPA’s PMF model was used to identify and quantify the contribution rates of different pollution sources in the study area. Results are shown in Figure 3.

Factor 1 was dominated by Hg (74.0%) and Cd (43.5%). The high coefficient of variation for Hg (540%) indicates that Hg pollution is significantly affected by human activities. It has been reported that Hg in soils mainly originates from industrial sources, especially petrochemical production and coal combustion [46, 51, 52]. Similarly, Cd also enters the environment through human activities, via exhaust gas and coal combustion emissions, owing to the stable bond between Cd and organic matter and sulfide in coal [53, 54]. These sources are consistent with the distributions observed in our study. A large number of coal mines are exploited around Wuhai to provide energy for thermal power, iron and steel production, and cement production plants. Coal is also extensively used as a raw material in the coal chemical industry and in metal smelting [55]. In general, the concentrations of Hg and Cd in soils are mainly related to industrial activities, especially coal mining, coal chemical engineering, and metal smelting. Therefore, Factor 1 reflects the input from industrial activities.

Factor 2 was dominated by Arsenic (As) (69.7%). Arsenic (As) is commonly used in fertilizers as an additive [56, 57]; therefore, it is routinely introduced into soils during fertilizer applications. Arsenic (As) also is an effective ingredient in various herbicides and pesticides used regularly in agricultural production [58, 59]. Animal excrement may be another source of Arsenic (As), because antibiotics used for poultry contain it [60]. Clearly, Arsenic (As) is mainly linked to agricultural activities, occurring within the large area of farmland and grassland in the southwestern part of the study area. These long-term activities may lead to ongoing enrichment of Arsenic (As) in local soils, which will aggravate local Arsenic (As) pollution. Therefore, Factor 2 reflects the input of agricultural activities.

Factor 3 was dominated by Cr (51.2%) and Ni (38.6%). Our descriptive statistical analysis revealed small coefficients of variation for both Cr and Ni, with similar concentrations measured between background and sample values. Differences in concentration between Cr and Ni may be related to variability in the parent materials of the soils. Previous studies have shown that variations in Cr and Ni in different soils are inherited from the parent material [15, 22, 61]. Hence, Factor 3 reflects inputs from mainly natural sources.

Factor 4 was weighted on Pb (78.4%), Zn (44.7%), and Cu (32.0%), which are all affected by human activities, according to our descriptive statistical analysis. Many previous studies reported that the main sources of Pb in soils are vehicle exhaust emissions and coal combustion [6264]. Gao et al. found that Pb and Zn were derived from atmospheric deposition of dust from open pit coal mines [65]. Li et al. found that Cu, Zn, and Pb were all released in the smelting process of sulfide minerals (sphalerite and galena) [15]. Given the abundant coal resources and the long-term mining and processing operations in the study area, several heavy metals contained in the coal enter the environment as dust. In addition, heavy metals contained in flue gas are discharged into the atmosphere by many coal-burning industries; these enter the soil via wet and dry deposition processes. Furthermore, growth in the transport sector within the study area has resulted in a large amount of vehicle emissions containing Pb; this Pb constantly enters the soil via atmospheric deposition. Thus, Factor 4 reflects the contribution of atmospheric deposition.

To summarize, Hg and Cd point to sources of industrial activities, while Cu, Zn, and Pb are indicators of sources of atmospheric deposition. Moreover, arsenic (As) is indicative of the sources of agricultural activities, while Cr and Ni point to natural sources. Based on the factor fingerprint of each heavy metal, the overall percent contribution of each source was computed. Industrial activities (31.3%) contributed most of the heavy metal content to Wuhai soils, followed by atmospheric deposition (26.3%), agricultural activities (21.9%), and natural sources (20.5%). Overall, anthropogenic sources were the predominant factors, reflecting inputs from coal mining, coal chemical engineering, metal smelting, coal burning, vehicle emissions, fertilization, pesticide spraying, and other activities, yielding a total contribution of approximately 80%.

4. Conclusions

This study investigated contamination levels, spatial distribution characteristics, and percentage contributions of various sources of heavy metals in soils in Wuhai, China, by using descriptive statistical analysis, ordinary kriging interpolation, and EPA’s PMF model. Generally, mean concentrations of heavy metals in soils in Wuhai exceeded their background values, with exceedance ratios of As, Hg, Cr, Ni, Cu, Zn, Cd, and Pb being 100%, 29.3%, 100%, 98.3%, 96.6%, 39.6%, 100%, and 100%, respectively. Their coefficients of variation ranged from 15.9% to 540%. The most serious heavy metal contaminations in this area were associated with Hg and Cd. According to comparisons between soils from typical coal mining and nonmining cities throughout China, Wuhai soils were moderately contaminated by heavy metals. Spatial distributions and source apportionment suggest that heavy pollution levels of Hg, Ni, Cu, Zn, Cd, and Pb proximal to the Wuhushan coalfield of the Wuda district of Wuhai are linked to significant inputs from coal mining activities. Meanwhile, Cu, Zn, and Pb contaminations were mainly derived from fuel and coal combustion followed by atmospheric deposition. Hg and Cd contaminations were mainly derived from industrial activities, including coal mining, coal chemical engineering, and metal smelting. The Arsenic (As) contamination was mainly related to agricultural activities, with inputs from fertilizers, herbicides, pesticides, and animal excrement. In contrast, parent materials played a vital role in producing Cr and Ni enrichments in these soils. Generally, anthropogenic sources were the predominant factors influencing heavy metal pollution in this area.

Data Availability

The numerical data used to support the findings of this study are included within the article and the associated supporting information.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was cosupported by the National Key R&D Program of China (2017YFC0504401), the National Natural Science Foundation of China (21777188), and the Fundamental Research Funds for the Central Universities in China (2015QH02). Q. B. was financially supported by Yue Qi Young Scholar Project, China University of Mining & Technology, Beijing (2017QN15).

Supplementary Materials

Table S1: summary of concentrations of soil heavy metals collected from typical coal mining cities and nonmining cities in recent years (mg/kg). Table S2: correlation coefficients between different heavy metals. (Supplementary Materials)