Spatial Distributions, Pollution Assessment, and Qualiﬁed Source Apportionment of Soil Heavy Metals in a Typical Mineral Mining City in China

: Daye is a city in China known for its rich mineral resources, with a history of metal mining and smelting that dates back more than 3000 years. To analyze the spatial distribution patterns, ecological risk, and sources of heavy metals (Cd, Co, Cr, Cu, Mn, Ni, Pb, and Zn) in soils, 213 topsoil samples were collected in the main urban area of Daye in September 2016. The mean concentrations of Cd, Cu, Pb, and Zn were higher than the corresponding background values, with the mean concentration of Cd being almost seven times its background value. Spatially, the high concentrations of Cd, Mn, Pb, and Zn were mainly concentrated in the southeastern part of the region due to nonferrous metal mining and smelting. However, the high concentrations of Co and Cu were concentrated in the central part of the study area, resulted from copper mining and smelting. The data of the geoaccumulation index showed that the contamination levels ranged from no pollution (Co, Cr, Mn, and Ni) to heavy contamination (Cd, Cu, and Pb). Ecological risk assessment showed that Cd posed a high, serious, and even severe ecological risk in 53.78% of the area of Daye. According to the results of the principal component analysis, mineral exploitation and smelting involving a variety of minerals (ES_M), mining exploitation, and smelting of copper ore (ES_C), and natural sources are the three main sources of heavy metals in these soils. Furthermore, the absolute principal component scores showed that 69.21% and 23.17% of the heavy metal concentrations were ascribed to ES_M and ES_C, respectively.


Introduction
With the increases in intense urbanization, industrialization, and human population, soil heavy metal pollution has been significant over the past several decades because of the toxicity and undegradability of these pollutants [1][2][3][4]. Many studies have shown that the accumulation of heavy metals in the soil degrades the quality of farmland, food crops, and humans' living environment and threatens human health via the food chain or dermal contact [5,6].
In China, there are 171 varieties of mineral resources, and the proven reserves of mineral resources in China constitute 12% of the total mineral resources in the world [7]. Existing studies have shown that heavy metal soil pollution is related to mining activities, including mining, smelting, and transportation. For example, Liang et al. reported that 26.05% of the heavy metals in soil were a result of mining activities in a typical coal mining city, Lianyuan, in the Hunan province of China [8]. Wang et al. reported that 65.7% of soil samples taken near a tungsten mine in the Dayu County of Jiangxi province, China, were polluted with heavy metals [9]. Xu et al. investigated paddy soils from lead-zinc mining

Study Area and Data
The city of Daye (latitude 29 • 40 -30 • 15 N, longitude 114 • 31 -115 • 20 E) is located in the southeastern area of the Hubei province along the south bank of the middle reaches of the Yangtze River. Due to its rich mineral resources, there are still numerous mines (coal, copper, iron, and gold) and smelting factories around Daye. To assess the soil heavy metal pollution, a total of 213 topsoil samples were collected in the main urban area of Daye during September 2016. The sampling points were randomly distributed in the study area based on a regular grid of 0.7 × 0.7 km, and each grid had at least one sampling point. The topsoil (0-20 cm) was collected by mixing 5 sub-samples obtained in different directions using a stainless steel hand spade. The spatial site of each sampling point was recorded using a hand-held positioning system (GPS). All soil samples were air-dried at room temperature, then passed through a 0.15 mm nylon sieve. These fine fractions were digested using a ternary mixture of concentrated acids consisting of 6ml HCO 4 , 5 mL HF, and 3 mL HNO 3 [17]. The concentrations of eight heavy metal elements (Cd, Co, Cr, Cu, Mn, Ni, Pb, and Zn) were measured using inductively coupled plasma mass spectrometry (ICP-MS). Quality assurance and quality control for metals in Sustainability 2018, 10, 3115 3 of 16 soils sample were carried out by the 20% analysis (in duplicate) of randomly chosen samples and was found to be within the accepted standard reference materials (GSS-3), which were obtained from the Center of National Standard Reference Material of China. The results were considered acceptable when the coefficients of variation were within 5%. All soil samples were analyzed at the Key Laboratory of Arable Land Conservation (middle and lower reaches of the Yangtze River), Chinese Ministry of Agriculture, Wuhan city. Land cover data were obtained by visual interpretation of contemporaneous Google remote sensing data. As shown in Figure 1, the mines are distributed throughout the southern and northern portions of the mountain forest area; some mineral smelters are distributed on the edge of the urban area; the factories, including those for mechanical and electrical manufacturing, textiles, clothing, electronics, new materials, and food and beverages, are located in the northern portion of the study area. The spatial distributions of soil sampling points and land cover are shown in Figure 1. The results were considered acceptable when the coefficients of variation were within 5%. All soil samples were analyzed at the Key Laboratory of Arable Land Conservation (middle and lower reaches of the Yangtze River), Chinese Ministry of Agriculture, Wuhan city. Land cover data were obtained by visual interpretation of contemporaneous Google remote sensing data. As shown in Figure 1, the mines are distributed throughout the southern and northern portions of the mountain forest area; some mineral smelters are distributed on the edge of the urban area; the factories, including those for mechanical and electrical manufacturing, textiles, clothing, electronics, new materials, and food and beverages, are located in the northern portion of the study area. The spatial distributions of soil sampling points and land cover are shown in Figure 1.

Data Analysis and Spatial Distribution Maps
Standard statistical analysis was carried out to describe the soil heavy metal concentrations using Excel software (Version 14, Microsoft, Inc., Redmond, Washington, DC, USA). Then, Pearson correlation analysis based on sampling points was used as an indicator of the relationships between different heavy metals in soils using IBM SPSS Statistics software (Version 19.0, Armonk, NY, USA). Ordinary Kriging (OK) interpolation was performed on the concentrations of soil heavy metals using ArcGIS software (Version 10.2, ESRI, Inc., Redlands, CA, USA). The standard sample-based cross-validation technique was implemented to assess the spatial interpolation accuracy. Based on the leave one-out cross validation technique, two accuracy indicators were computed from the pairs of "estimated-observed" soil heavy metal concentrations at the sampling points: The Pearson correlation coefficient (r) and the mean absolute error (MAE). The r and MAE were used to measure the strength of the linear relation and mean absolute deviation between the estimated and the observed soil heavy metal concentrations, respectively. For an accurate spatial interpolation, r should be close to 1, and the MAE should be as small as possible. In addition, to remove any quantitative differences among the various heavy metals in the soils, an error rate (ER) was defined as the ratio of the MAE over the mean value of the heavy metals.

Methods for the Assessment of Pollution and Potential Ecological Risk
The geoaccumulation index ( ) was used to assess metal contamination levels because is the most popular index for pollution evaluation based on single heavy metals in soils [18]; is defined as follows [19]:

Data Analysis and Spatial Distribution Maps
Standard statistical analysis was carried out to describe the soil heavy metal concentrations using Excel software (Version 14, Microsoft, Inc., Redmond, Washington, DC, USA). Then, Pearson correlation analysis based on sampling points was used as an indicator of the relationships between different heavy metals in soils using IBM SPSS Statistics software (Version 19.0, Armonk, NY, USA). Ordinary Kriging (OK) interpolation was performed on the concentrations of soil heavy metals using ArcGIS software (Version 10.2, ESRI, Inc., Redlands, CA, USA). The standard sample-based cross-validation technique was implemented to assess the spatial interpolation accuracy. Based on the leave one-out cross validation technique, two accuracy indicators were computed from the pairs of "estimated-observed" soil heavy metal concentrations at the sampling points: The Pearson correlation coefficient (r) and the mean absolute error (MAE). The r and MAE were used to measure the strength of the linear relation and mean absolute deviation between the estimated and the observed soil heavy metal concentrations, respectively. For an accurate spatial interpolation, r should be close to 1, and the MAE should be as small as possible. In addition, to remove any quantitative differences among the various heavy metals in the soils, an error rate (ER) was defined as the ratio of the MAE over the mean value of the heavy metals.

Methods for the Assessment of Pollution and Potential Ecological Risk
The geoaccumulation index (I geo ) was used to assess metal contamination levels because I geo is the most popular index for pollution evaluation based on single heavy metals in soils [18]; I geo is defined as follows [19]: where C n and B n are the concentration of heavy metal, n, in the soil and the geochemical background concentration of the corresponding heavy metal, n, respectively, in the Hubei province. Unlike calculations in previous studies, I geo was calculated based on the spatial distribution of various heavy metals instead of the concentrations in soil samples. That is, I geo was determined for each heavy metal at each location, and the overall picture of heavy metal pollution in the study area can be obtained in the present study. Based on the I geo value for each heavy metal, the heavy metal pollution was classified into 7 categories: I geo ≤ 0 (practically uncontaminated), 0 < I geo ≤ 1 (uncontaminated to moderately contaminated), 1 < I geo ≤ 2 (moderately contaminated), 2 < I geo ≤ 3 (moderately to heavily contaminated), 3 < I geo ≤ 4 (heavily contaminated), 4 < I geo ≤ 5 (heavily to extremely contaminated), and I geo > 5 (extremely contaminated). The potential ecological risk index (RI) was used to quantitatively express the potential ecological risk posed by a variety of toxic heavy metals. The RI is calculated by the following formulas [20,21]: where E i is the potential ecological risk factor for a given heavy metal, i, T i is the toxic coefficient of a heavy metal, i, C i is the concentration of the heavy metal, and C 0 is the background value of the corresponding heavy metal in the Hubei province. Based on previously published research, the toxic coefficients of Cd, Co, Cr, Cu, Mn, Ni, Pb, and Zn are 30, 5, 2, 5, 1, 5, 5, and 1, respectively [22]. Based on E i or the RI, the potential ecological hazard posed by heavy metals was classified into 5 grades (see Table 1). The calculations introduced in this section were completed using Excel software (Version 14, Microsoft Inc., Redmond, Washington, DC, USA).

Principal Component Analysis (Pca) and Absolute Principal Component Scores (APCS)
PCA/APCS is an effective multivariate factor analysis model for the identification of pollution sources. It was first applied to estimate sources of particulate matter in Boston by Thurston and Spengler [23]. The PCA/APCS receptor model can quantitatively determine not only the load of each pollutant attributed to each source, but also the average contribution of each source to each pollutant and to the pollution in each sample. Currently, the model is mainly used in research on air pollution and water pollution [24][25][26][27]. In this study, PCA was performed on the spatial distributions of heavy metal accumulation in soils for source identification using factor extraction with eigenvalues > 1 after varimax rotation. Then, APCS were used to estimate the source contributions to each heavy metal in the following steps: Step 1: The first step of determining APCS is the normalization of all heavy metal accumulation values as Z ik using the following equation: where A ik is the accumulation of the heavy metal, i, at the location, k, andĀ i and σ i are its mean value and standard deviation, respectively; Sustainability 2018, 10, 3115 5 of 16 Step 2: By introducing an artificial sample, with accumulations equal to zero for all heavy metals, the true zero for each factor score is calculated: Step 3: The factor scores of those normalized heavy metal accumulation values are obtained by the PCA method. The APCS of each heavy metal at each location are obtained by subtracting the factor scores for this artificial sample from the factor scores of each location; and Step 4: The source contributions are derived by regression using the following equation: where A i is the accumulation of the heavy metal, i, (b 0 ) i is the constant term of multiple regression for the heavy metal, i, b pi is the coefficient of multiple regression of the source, p, for the heavy metal, i, and APCS p is the scaled value of the rotated factor, p, for the considered location. APCS p × b pi is the contribution of the source, p, to A i . The mean of the product, APCS p × b pi , at all locations represents the average contribution of the sources. The calculations introduced in this section were completed using SPSS software (Version 19.0 for Windows, IBM Inc., Armonk, NY, USA) and Excel software (Version 14, Microsoft Inc., Redmond, WA, USA).

Descriptive Statistics and Spatial Distributions of Soil Heavy Metals
Based on the available samples, the statistical characteristics of the heavy metal concentrations in the Daye region are presented in Table 2. As shown in Table 2, the percentage of sampling points where the concentration exceeds the background (PbB) was greater than 50% for Cd, Cu, Pb, and Zn, indicating serious Cd, Cu, Pb, and Zn pollution in the study area. On the other hand, although the mean concentrations of Co and Mn were less than their background values, their PbB values exceeded 30%, indicating that part of the study area was polluted by these two heavy metals. The heavy metals in soils were also studied in different mines in China. As shown in Table 3, the concentration distribution patterns of Cd in this study were similar to the values found in the Nandan County, which is located in the north of the Guangxi province as a mining town [29]. However, the mean concentration of Cd was much lower than those reported from some large lead-zinc mining areas, such as the lead-zinc mine in Qixiashang of Nanjing city with a mean of 17.56 mg/kg [30], and the lead-zinc mine with a mean of 15.56 mg/kg in the Lanping county, Yunnan province, China [31].
Compared to soil from other mining areas, the mean value of Cu in Daye soils was markedly three to six times higher than the values reported in other urban soils [8,[31][32][33][34]. However, the mean value of Cu was still lower than the value in soils of the mining area in Nanjing city [30,35]. The measured mean values of Mn, Pb, and Zn were extremely higher in the soils of lead-zinc or Mn mining areas than the values presented in this study. For example, the mean concentrations of Mn, Pb, and Zn were 14015.31 mg/kg, 2232.16 mg/kg, and 5201.86 mg/kg, respectively, in the soils of the lead-zinc mine of Nanjing city [35]. The mean values of Pb and Zn were 1093.03 mg/kg and 867.08 mg/kg, respectively, in the soils of a lead-zinc mine in the Guangdong province [10]. The mean concentrations of Mn, Pb, and Zn were 8853.21 mg/kg, 1211.29 mg/kg, and 685.36 mg/kg, respectively, in a Mn mine in Xiangtan city, Hunan province [36]. According to the comparisons above, the mean concentration of Cu in soils was higher than most of the other mining areas in China. However, the concentrations of Cd, Mn, Pb, and Zn were lower than other mining areas, especially lead-zinc mines, in China. Thus, the concentrations of heavy metals in soils may be closely related to the local mineral type.
The Pearson's correlation coefficients between the concentrations of heavy metals in soils based on sampling points are summarized in Table 4. The concentration of Cd, Mn, Pb, and Zn showed a high significant positive relationship with each other, which may suggest a common origin. The correlation coefficient between Cu and Co is 0.62, indicating a common source. Additionally, the correlation coefficient between Cu and Co is 0.51, also suggesting a common origin. Furthermore, although there are some other significant positive relationships, such as Cr and Cd, Mn and Cr, Zn and Co, and so on, the correlation coefficients between them are not high. Thus, the results of the correlation analysis based on the sampling points cannot clearly indicate the sources of pollution.
The spatial distributions of the Cd, Co, Cr, Cu, Mn, Ni, Pb, and Zn concentrations in the study area are presented in Figure 2. To intuitively understand the severity of the heavy metal pollution levels in the soil of the Daye region, in the maps of Figure 2, the yellow, blue, and red colors represent heavy metal concentrations that are close to, lower than, and higher than their background values, respectively. As shown in Figure 2, the maps of the Cd, Cu, Pb, and Zn distributions are almost entirely covered with red and yellow colors, indicating that the concentrations of these four heavy metals were higher than their background values over almost the entire study region. Regarding the Cr and Ni maps, most of the study region was covered in blue, indicating that only a small part of the region was polluted with Cr and Ni. Based on the Co and Mn maps, it was concluded that nearly half of the study region was polluted with these elements. Furthermore, the spatial distribution patterns of all heavy metals were characterized by decreasing concentrations from the central and southeastern parts toward the southern and northern parts. These concentrations peaked in the central and southeastern part, which may be related to the mining and smelting operations in that region (see Figure 1).  Cu 0.28 ** 0.62 ** 0.2 ** 1 Mn 0.43 ** 0.32 ** 0.2 ** 0.32 ** 1 Ni 0.15 * 0.26 ** 0.51 ** 0.06 0.24 ** 1 Pb 0.86 ** -0.06 0.09 0.14 * 0.52 ** 0.13 1 Zn 0.8 ** 0.19 ** 0.25 ** 0.4 ** 0.65 ** 0.3 ** 0.9 ** 1 ** p < 0.01 (2-tailed), * p < 0.05 (2-tailed).    Table 5 summarizes the cross-validation results for the spatial distributions of heavy metal concentrations obtained by OK interpolation. These results show that the r and ER values of the eight heavy metal concentrations ranged from 0.75 to 0.99 (p < 0.01) and from 4.75% to 18.42%, respectively, indicating high spatial interpolation accuracy (note the high r and low ER values). Hence, the spatial interpolation results were valuable for subsequent analyses.

Pollution and Ecological Risk Assessment
Obviously, the spatial I geo distribution pattern of each heavy metal is similar to the spatial pattern of the corresponding heavy metal distribution in soils, as shown in Figure 2. Therefore, it is not necessary to present the I geo for each heavy metal spatial distribution here. According to the quantitative I geo results for the eight heavy metals (Figure 3), 74.23%, 77.16%, 30.46%, and 94.62% of the study area was contaminated by Cd, Cu, Pb, and Zn, respectively, at various pollution levels. Approximately 0.90%, 3.82%, and 0.44% of the study area was heavily contaminated by Cd, Cu, and Pb, respectively. Additionally, approximately 38.40%, 31.72%, 6.86%, and 29.09% of the study area was moderately contaminated by Cd, Cu, Pb, and Zn, respectively. These findings demonstrate that the soils of the study area were partially seriously polluted by some heavy metals (e.g., Cd, Cu, Pb, and Zn).
As shown in Figure 3b, Cd posed a moderate, high, serious, and severe potential ecological risk in approximately 23.30%, 36.80%, 15.39%, and 1.59% of the study area, respectively, whereas almost all of the study area was free from ecological risk posed by other heavy metals, especially Co, Cr, Mn, Ni, and Zn. Thus, it can be concluded that the potential ecological hazard due to heavy metals in Daye soils was mostly derived from the accumulation of Cd. As presented in Figure 3c, the RI values exceeded 150 throughout the central part of the study area. The spatial distribution pattern of the RI was similar to the pattern of Cd, indicating that the potential ecological hazard due to soil heavy metals was mainly related to Cd in Daye City.  Table 5 summarizes the cross-validation results for the spatial distributions of heavy metal concentrations obtained by OK interpolation. These results show that the r and ER values of the eight heavy metal concentrations ranged from 0.75 to 0.99 (p < 0.01) and from 4.75% to 18.42%, respectively, indicating high spatial interpolation accuracy (note the high r and low ER values). Hence, the spatial interpolation results were valuable for subsequent analyses.

Pollution and Ecological Risk Assessment
Obviously, the spatial distribution pattern of each heavy metal is similar to the spatial pattern of the corresponding heavy metal distribution in soils, as shown in Figure 2. Therefore, it is not necessary to present the for each heavy metal spatial distribution here. According to the quantitative results for the eight heavy metals (Figure 3), 74.23%, 77.16%, 30.46%, and 94.62% of the study area was contaminated by Cd, Cu, Pb, and Zn, respectively, at various pollution levels. Approximately 0.90%, 3.82%, and 0.44% of the study area was heavily contaminated by Cd, Cu, and Pb, respectively. Additionally, approximately 38.40%, 31.72%, 6.86%, and 29.09% of the study area was moderately contaminated by Cd, Cu, Pb, and Zn, respectively. These findings demonstrate that the soils of the study area were partially seriously polluted by some heavy metals (e.g., Cd, Cu, Pb, and Zn).
As shown in Figure 3b, Cd posed a moderate, high, serious, and severe potential ecological risk in approximately 23.30%, 36.80%, 15.39%, and 1.59% of the study area, respectively, whereas almost all of the study area was free from ecological risk posed by other heavy metals, especially Co, Cr, Mn, Ni, and Zn. Thus, it can be concluded that the potential ecological hazard due to heavy metals in Daye soils was mostly derived from the accumulation of Cd. As presented in Figure 3c, the RI values exceeded 150 throughout the central part of the study area. The spatial distribution pattern of the RI was similar to the pattern of Cd, indicating that the potential ecological hazard due to soil heavy metals was mainly related to Cd in Daye City. (a)

Source Identification and Apportionment By PCA/APCS
The PCA results obtained by applying varimax rotation to the concentrations of heavy metals in the soil are shown in Table 6. Three principal components (PCs) with eigenvalues higher than 1 (before and after rotation) were extracted; PCA led to a reduction in the initial dimension of the dataset to three components that explain 87.21% of the data variation. The first PC (PC1) explains 38.98% of the total variance and with high loadings for Cd, Mn, Zn, and Pb. As shown in Figure 2, the hot spots of Cd, Mn, Zn, and Pb were all located in the southeastern part of the study area, where many mines and smelters of gold, silver, copper, iron, molybdenum, and sulfur, among other minerals, are located. Thus, PC1 can be considered a source of mineral exploitation and smelting involving a variety of minerals (ES_M). The second PC (PC2), dominated by Co and Cu and moderately by Cd, accounts for 25.74% of the total variance. As shown in Figure 2, the concentrations of Co and Cu have similar spatial distribution patterns characterized by an increasing accumulation from the north and south to the center of the study area. The only hot spot is located in the central part of the study area, where the oldest in-use copper mine (Tonglu Mountain Copper Mine) in China, with related tailings and smelters, is located. Cobalt (Co) resources are mostly associated with copper deposits. Thus, PC2 is mainly related to the mining

Source Identification and Apportionment By PCA/APCS
The PCA results obtained by applying varimax rotation to the concentrations of heavy metals in the soil are shown in Table 6. Three principal components (PCs) with eigenvalues higher than 1 (before and after rotation) were extracted; PCA led to a reduction in the initial dimension of the dataset to three components that explain 87.21% of the data variation. The first PC (PC1) explains 38.98% of the total variance and with high loadings for Cd, Mn, Zn, and Pb. As shown in Figure 2, the hot spots of Cd, Mn, Zn, and Pb were all located in the southeastern part of the study area, where many mines and smelters of gold, silver, copper, iron, molybdenum, and sulfur, among other minerals, are located. Thus, PC1 can be considered a source of mineral exploitation and smelting involving a variety of minerals (ES_M). The second PC (PC2), dominated by Co and Cu and moderately by Cd, accounts for 25.74% of the total variance. As shown in Figure 2, the concentrations of Co and Cu have similar spatial distribution patterns characterized by an increasing accumulation from the north and south to the center of the study area. The only hot spot is located in the central part of the study area, where the oldest in-use copper mine (Tonglu Mountain Copper Mine) in China, with related tailings and smelters, is located. Cobalt (Co) resources are mostly associated with copper deposits.
Thus, PC2 is mainly related to the mining exploitation and smelting of copper ore (ES_C). The third PC (PC3) is strongly correlated with Cr and Ni. The concentrations of Cr and Ni were lower than or approximately equal to their background values in Hubei soil, indicating that these elements mainly originate from natural sources.
Next, the APCS receptor model was applied to the three extracted components to quantify the contributions of the sources to each spatial grid. SPSS software was used to develop suitable programs and outputs for carrying out the apportionment procedure. The results of the statistical analysis are presented in Table 7. The high values of R 2 for all heavy metals demonstrate that the source apportionment results meet the modeling requirements. Figure 4 shows the spatial distributions of various heavy metal concentrations that resulted from the three major sources. Notably, in the results of the source apportionment based on APCS, the source contribution estimates are not constrained to be non-negative [26,37,38]. Thus, the negative values in the results of the APCS source apportionment, including the negative values in Figure 4, indicate relatively small contributions. Furthermore, the contributions from non-specified sources (the forth to eighth components in Table 6) were negligible because the estimated values of b 0 in Equation (5) were usually very small. Based on the percentage of individual heavy metals apportioned to each source (Table 7), the fraction of the contribution of each source to the total heavy metal concentration in soils was estimated (see the last line in Table 7).     According to the results in Table 7 and Figure 4, in Daye City, ES_M accounted for an average of 51.97%, 82.69%, 70.95%, and 69.83% of Cd, Mn, Pb, and Zn, respectively. The high contributions made by the ES_M source were located in the western and southeastern parts of the study area, where various kinds of mines and smelters are distributed. The ES_M source contributed 69.21% of the total soil heavy metals in the study area. Therefore, ES_M is the largest source of soil heavy metal pollution in Daye. On the other hand, ES_C accounted for an average of 79.65% and 83.4% of the Co and Cu, respectively. The hotspots of the contributions by ES_C are located in the western and central parts of the study area, which have a long history of copper mining and related smelting enterprises. This again indicates that the second greatest source of heavy metal pollution was directly related to the mining and smelting of copper.
Meanwhile, as is shown in Figure 4, the high contributions for all heavy metals by source 1 and source 2 were located in or around the mines and smelters (also refer to the spatial distributions of land cover in Figure 1). It is worth noting that there are the same spatial distribution patterns of contribution maps for all heavy metals from the same source. For example, Figure 4a,d,g,j,m,p,s,v shows spatial distributions of the various heavy metals' concentrations contributed by ES_M. Those maps clearly point out three locations of ES_M sources (the north, west, and southeast parts of the study area) and four kinds of heavy metals mainly contributed by the ES_M source (Cd, Mn, Pb, and Zn). Similarly, Figure 4b,e,h,k,n,q,t,w indicates that ES_C is located in the center of the study area, and it is mainly responsible for the accumulations of Cu and Co. In addition, as is shown in Figure 4j, ES_M also contributed Cu to the soil. In addition, a part of Cd, Mn, and Zn in the soils were contributed by ES_C, as shown in Figure 4b,q,w; thus, the contribution maps of the identified sources of heavy metals might be a powerful tool to analyze the land covers that contaminate surrounding soils.

Conclusions
The present study identified the pollutant characteristics and sources of heavy metals in soils of Daye, China. It was based on a spatial distribution analysis, pollution assessment, and PCA/APCS methods. There are two main conclusions: (1) The soils in Daye are seriously polluted by the heavy metals, Cd, Cu, Pb, and Zn. The mean concentrations of Cd, Cu, Pb, and Zn were higher than the corresponding background values, with the mean concentration of Cd being almost seven times more than the background value. The results of the correlation analysis showed that Cd, Mn, Pb, and Zn have a high significant positive relationship with each other. In addition, the correlation between Cu and Co was highly significant and positive. The I geo results showed that the contamination levels ranged from no pollution (Co, Cr, Mn, and Ni) to moderate contamination (Cd, Cu, Pb, and Zn) and heavy contamination (Cd, Cu, and Pb). More than 90% of the study area was polluted by at least one kind of heavy metal, and the accumulation of Cd posed the most potential ecological risk among all kinds of heavy metals; and (2) There were three main factors influencing the concentrations of soil heavy metals: ES_M, ES_C, and natural sources. Among these sources, ES_M was foremost, explaining 57.21% of the total variance. The APCS showed that 69.21% and 23.17% of the heavy metal concentrations were ascribed to ES_M and ES_C, respectively. Those sources were significantly related to the mines and smelters scattered in the study area.
Based on this study, we suggest a few steps to be taken to solve this serious polluted soil issue: (1) Soil remediation is urgently needed in Dye; (2) mines and smelters in the studied area should be controlled strictly to prevent further pollution of heavy metals in the soils; (3) the concentration of heavy metal in local agricultural products should be strictly detected before entering the market; and (4) local residents should avoid direct contact with soils in heavily polluted area as far as possible.