Evaluating Metal(loid)s Contamination in Soil of a Typical In-Dustry Smelting Site in South Central China: Levels, Possible Sources and Human Health Risk Analysis

: Soil pollution by metal(loid)s caused by smelting activities is a severe problem posing a great threat to environmental and human health. In this study, the concentrations, sources and human health risks posed by six potentially toxic elements (Cr, Mn, Zn, Pb, Cd, and As) were determined in the soil of a typical alloy smelting site in South Central China. The results showed that the concentrations of metal(loid)s were in the descending order of Mn > Cr > Zn > Pb > As > Cd. Addi-tionally, the selected elements were classiﬁed into different pollution degrees using geo-accumulation index and pollution load index. The entire study area had a high pollution level with relatively severe and extensive contamination by Cr, Cd, and Mn. The combined application of principal component analysis and positive matrix decomposition model revealed that the major sources of these elements include smelting activities (48.68%), waste residue stacking (22.95%), and natural sources (28.37%). According to the results of the human health risk assessment, the non-carcinogenic risk was insigniﬁcant as a whole; however, the carcinogenic risk had an unacceptable level. Among them, Cr was the main driver of carcinogenic risk, which needs special attention.


Introduction
The metal(loid)s contamination of soils has become an increasingly severe concern in many parts of the world, it has a profound effect on the ecological environment and human health [1][2][3][4]. The metal(loid)s in the environment usually come from industrial activities and agricultural production [1,[5][6][7]. As one of the largest global metal producers and consumers [8,9], China has had rapid industry development, creating much wealth from non-ferrous metal smelting in the past few decades [8]. However, smelting activities have also brought a heavy burden to the surrounding environment, especially to soil [10]. Fly ash, wastewater, and slag, which are side products of smelting processes, such as material preparation, smelting, and casting processes, contain toxic elements that enter the soil environment through rainfall, sedimentation, leakage, and other means, thus causing severe pollution [11,12]. Due to the obvious toxicity, difficult degradation, and outstanding persistence of metal(loid)s [13,14], it deteriorates soil quality and enters the human body through food intake, skin contact, and other ways, inducing diseases such as cancer and harming human health [15][16][17][18].
Mounting evidence indicates that people's exposure to these potentially toxic metal(loid) elements will adversely affect many human organs [19]. For instance, Cr(VI) is a harmful carcinogenic pollutant with high toxicity, which causes diseases of the respiratory tract, Accordingly, the present study focuses on: (1) determining the contamination levels and distribution of metal(loid)s in the selected smelter area; (2) identifying the possible sources of metal(loid)s in the soil; (3) assessing the human health risks through different exposure pathways. The outcomes of this study could provide basic data for metal(loid)s pollution control and for policy decisions based on the accurate confirmation of regional pollutants. In addition, this study provides a reference case for the systematic evaluation of soil contamination, which can be used for other related studies in similar smelting sites around the world.

Study Area and Sampling
The study area, which is an abandoned smelting factory located in Hunan Province, South Central China, belongs to the subtropical monsoon humid climate, has four distinct seasons. Founded in 1958, the plant was once one of the key ferroalloy enterprises in China, with an annual production capacity of 148,000 tons. Production stopped recently and the factory closed down. The industrial site primarily consists of the NO. 3 smelting workshop, NO. 4 smelting workshop and blast furnace smelting workshop, as shown in Figure 1, which mainly produced Fe-Mn products and chromium-related products. Due to the simple traditional smelting process, remedial nature of environmental treatment measures and the random discharge of industrial pollutants, the soil in the study area has been heavily polluted.

Sample Pretreatment and Analysis
The soil samples were aired-dried at room temperature in an indoor ventilated place, ground thoroughly with an agate mortar, and passed through a 200-mesh (0.75 μm) nylon sieve. To measure the concentration of metal(loid)s, 0.1 g of each soil sample was placed in a crucible and digested in 2 mL HNO3 and 2 mL HCl on a hot plate for 1 h at 110 ℃. Afterwards, 1 mL HF and 1 mL HClO4 were added, the temperature was increased to 250 ℃, and digestion was continued until the white smoke disappeared. Once the sample was allowed to cool down, it was diluted with 2 mL of HCl (1 + 1) and transferred into a 10 mL polyethylene container. The resulting solution was diluted with deionized water to Based on the production history and prevalent processes at the time of operation, 40 sampling points were set up on the site, as shown in Figure 1. The surface soil (20-50 cm) was collected within a single sampling site using a shovel, and then mix them thoroughly. Subsequently, 1000 g of the mixed sample was put it into a sealed bag. Before loading soil samples into the sealed bags, stones, plant debris, and other sundries on the soil surface were discarded. After all the sampling work was completed, it was sent to the laboratory for analysis.

Sample Pretreatment and Analysis
The soil samples were aired-dried at room temperature in an indoor ventilated place, ground thoroughly with an agate mortar, and passed through a 200-mesh (0.75 µm) nylon sieve. To measure the concentration of metal(loid)s, 0.1 g of each soil sample was placed in a crucible and digested in 2 mL HNO 3 and 2 mL HCl on a hot plate for 1 h at 110°C. Afterwards, 1 mL HF and 1 mL HClO 4 were added, the temperature was increased to 250°C, and digestion was continued until the white smoke disappeared. Once the sample was allowed to cool down, it was diluted with 2 mL of HCl (1 + 1) and transferred into a 10 mL polyethylene container. The resulting solution was diluted with deionized water to the scale mark, and the concentrations of Pb, and Cd were determined by ICP-MS (Thermo fisher X2, Bremen, Germany), while Cr, Mn, Zn and As were analyzed by ICP-AES (ICAP6300, USA). All readings were taken in three replicates.

Quality Assurance and Quality Control
The accuracies of the soil sample analysis were controlled using the standard reference materials (GBW07429 and GBW07404) of the National Certified Reference Materials Research Center of China. All analyses were carried out in triplicate, and the recoveries of each element were were within the range of 89-106%. The relative standard deviations of all elements ranged from 1% to 7%, and the analytical precision was within ±10%.

Pollution Indices
In order to evaluate the pollution level of selected elements in the soil of the study area, I geo and PLI was introduced.
I geo is a classical method widely used for the soil pollution level of selected elements [40,41]. The relevant equation is as follows: where C n represents the measured concentration value; B n represents the background value of elements, and in this study, it denotes the soil background value of Hunan Province [42]; K is the change of background value caused by diagenesis (generally k = 1.5). I geo can be classified into seven categories, as shown in Table S1. PLI is a screening tool to evaluate the extent of pollutants of a particular location [43], which has been widely used to measure the cumulative pollution degree [44,45]. The following formula was used to calculate the PLI: where CF represents the pollution factor of single element, C n is the measured concentration value, and B n represents the background value. The values of PLI are generally classified into six categories, as shown in Table S1.

Multivariate Statistical Analysis
PCA is a dimension reduction method widely used to distinguish the different origin of metal(loid)s groups [36]. PMF combined with PCA could quantitatively identify pollution sources and their contributions [44]. PMF is an effective source analysis model, which mainly uses the minimum iterative double multiplication to determine the main pollution sources and contribution rate. PMF mainly used uncertainty to weight all data [34,36]. The relevant equation is as follows: where x ij represents the measured concentration of the selected element; g ik represents the contribution matrix of the k th source factor; f kj represents the source profile of element for the k th source factor; e ij represents the residual for each sample; i and j are the number of samples and chemical elements. The minimum value of the objective function Q was calculated by the following equation: The uncertainty was calculated using the following formula: where Unc refers to the uncertainty, and MDL is the detection limit. If the concentration value ≤ MDL, the uncertainty is calculated using Equation (6); if the concentration value > MDL, Equation (7) is used.

Human Health Risk Assessment
The human health risk assessment was conducted according to the Chinese guideline (HJ 25.3-2019) [39] and the "Exposure Factors Handbook of US Environmental Protection Agency" [37], including the carcinogenic risk assessment and the non-carcinogenic risk assessment. The health risk assessments of the present study mainly considered three different pathways: ingestion, dermal contact and inhalation. The cancer risk (CR) was used to evaluate the carcinogenic risk, while the hazard quotient (HQ) index was employed to evaluate the non-carcinogenic risk. Equations (8)- (11) were used to calculate the carcinogenic risk, while Equations (12)- (15) were applied to calculate the non-carcinogenic risk.
CR n = CR ois + CR dcs + CR pis (11) where C sur denotes the measured value of metal(loid)s; CR ois , CR dcs , and CR pis represent carcinogenic risk via ingestion, dermal contact, and inhalation, respectively; OISER ca , DCSER ca , PISER ca represent the soil exposure dose (carcinogenic) of three exposure pathways, and their equations are given in Table S2. The total carcinogenic risk (CR n ) is the sum of risks posed by the three pathways. The other corresponding parameter values were summarized in Tables S2 and S3. According to previous studies, when CR n < 10 −6 , the carcinogenic risk can be ignored; when 10 −6 < CR n < 10 −4 , it is generally considered acceptable; and when CR n > 10 −4 , it is considered unacceptable [6,14].
H I = HQ ois + HQ dcs + HQ pis (15) where HQ ois , HQ dcs , and HQ pis are the hazard quotients of exposure to metal(loid)s via ingestion, dermal contact, and inhalation, respectively; OISER nc , DCSER nc , and PISER nc represent soil exposure dose (non-carcinogenic) via three pathways, as shown in Table S2. The hazard index (HI) represents the sum of HQ of different pathways for individual elements, indicating the total risk, whereas other parameters are presented in Tables S2 and S3. When HI < 1, it is generally considered that there is no significant non-carcinogenic risk, and when HI > 1, it indicates that there is a non-carcinogenic human health risk.

Statistics Used
The descriptive statistical analysis was conducted using SPSS 25.0 (IBM, Armonk, NY, USA) and Microsoft Excel 2010 (Microsoft, Redmond, WA, USA). Descriptive data were calculated, such as the minimum, maximum, mean and coefficients of variation. The Pearson's correlation and PCA analysis were also carried out by SPSS 25.0. Origin 2018 (OriginLab, Northampton, MA, USA) was used to convert the data into graphics to make the data more visual. The PMF model analysis was calculated by the EPA PMF 5.0 (United States Environmental Protection Agency, Washington, DC, USA) to explore the sources of metal(loid)s. All spatial distribution maps were made through kriging interpolation using Surfer 16.0 (Golden Software LLC, Golden, CO, USA).

Concentration of Metal(loid)s in the Soil
The average concentrations of Cr, Mn, Zn, Pb, Cd, and As in the study area soil were 1930, 13 (Table 1). Furthermore, the excess level in descending order was Cd > Cr > Mn > Pb > Zn > As. The coefficient of variation (CV) can reflect the effect of soil properties by the surrounding environment [46]. Generally, CV is classified into four types: CV < 20%, low variability; CV from 20% to 50%, moderate variability; CV from 50% to 100%, high variability; CV > 100%, exceptionally high variability [47,48]. The CV values of Cr, Mn, Zn, Pb and Cd were all greater than 100%, and this exceptionally high variability indicated that they were strongly impacted by external factors and human activities ( Table 1). The CV of As in the soil was 47.01%, exhibiting moderate variability, indicating that the pollution caused by As might be influenced by human activities, but was not excessively strong.
In order to figure out the pollution level, the metal(loid)s were evaluated by the I geo and PLI indices. The ranges of I geo for Cr, Mn, Zn, Pb, Cd, and As were from −0.61 to 7.40, −3.01 to 6.73, −0.66 to 3.98, −0.93 to 5.63, −2.57 to 6.80, and −1.31 to 1.54, respectively ( Figure 2). The average I geo value of the six elements in the soil decreased in the following order: Cr (3.37) > Cd (3.24) > Mn (2.94) > Pb (1.51) > Zn (0.81) > As (0.03). Based on the average values, it could be found that, in terms of Cr and Cd, the soil was classified as heavily polluted, and 52.5% of Cr and 61.1% of Cd exhibited I geo values greater than 3 ( Figure 2). Based on the Mn, Pb, and Zn average values, the soil was evaluated as moderately polluted to heavily polluted, moderately polluted, and unpolluted to moderately polluted, respectively. The average I geo value of As showed that industrial soils were unpolluted to moderately polluted, the pollution degree was the lightest, and only one sample had an I geo value above 1. The results showed that six elements caused pollution at different levels, among which Cr, Mn and Cd pollutions were comparatively serious, while As pollution was only slight ( Figure 2). The PLI values of Cr, Mn, Zn, Pb, Cd, and As were 15.49, 11.54, 2.64, 4.29, 14.15, and 1.53, respectively, and the PLI values were also in the descending order of Cr > Cd > Mn > Pb > Zn > As. Overall, the PLI values in all of the soil samples varied from 0.75 to 39.53 with an average of 7.83, showing that the industrial soil was strongly polluted by metal(loid)s.

Spatial Distribution of Metal(loid)s
The spatial distribution of metal(loid)s in this study area are presented in Figure 3. The spatial distribution of Zn, Pb, and Cd exhibited similar patterns, and particularly high concentrations of them accumulated near the blast furnace smelting workshop in the southwest corner and the NO.3 smelting workshop (Section 304), which indicated that the source of Zn, Pb, and Cd might be related to the former industrial activities in these two areas. The spatial distribution of Mn was concentrated in its slag landfill in the southeast and the NO.3 smelting workshop (Section 304). The hot spots of Cr were located in the chromium slag yard in the central area. However, the spatial pattern of As was distinct from the other metal(loid)s; its spatial variation was not significant, which was consistent with the analysis of CV. The concentration of As in the site was relatively uniform and was slightly higher in Section 304 in the north. Overall, most of the hot spots of these

Spatial Distribution of Metal(loid)s
The spatial distribution of metal(loid)s in this study area are presented in Figure 3. The spatial distribution of Zn, Pb, and Cd exhibited similar patterns, and particularly high concentrations of them accumulated near the blast furnace smelting workshop in the southwest corner and the NO. 3 smelting workshop (Section 304), which indicated that the source of Zn, Pb, and Cd might be related to the former industrial activities in these two areas. The spatial distribution of Mn was concentrated in its slag landfill in the southeast and the NO. 3 smelting workshop (Section 304). The hot spots of Cr were located in the chromium slag yard in the central area. However, the spatial pattern of As was distinct from the other metal(loid)s; its spatial variation was not significant, which was consistent with the analysis of CV. The concentration of As in the site was relatively uniform and was slightly higher in Section 304 in the north. Overall, most of the hot spots of these elements were clearly associated with the NO. 3 smelting workshop (Section 304), which had mainly produced chromium-related products, including roasting, reduction, calcination, smelting and other processes, with Section 304 mainly responsible for the smelting process. Therefore, the smelting process of Section 304 could provide a reference for determining the source of these elements.

Statistical Analysis
Correlation analysis (CA) and PCA can reflect the correlation among metal(loid)s levels to some extent and distinguish the different groups of metal(loid)s [44], which was used to reveal the association of the source of metal(loid)s in the study area. A strong correlation can reveal the close relationship between elements, indicating that these elements may have common sources and similar migration pathways [35,49], which can provide a reference for the traceability of metal(loid)s. The relationships among Cr, Mn, Zn, Pb, Cd, As, and pH were examined by Pearson's correlation analysis, as shown in Table  2. The correlations among Pb, Zn and Cd were highly significant at the p < 0.01 level (r > 0.77). Similarly, Mn with Cd and Zn were strongly correlated at the p < 0.01 level (r > 0.53), indicating that these elements might share a common pollution source. And the correlation between As and Cd was significant at the p < 0.05 level. In particular, Cr was correlated with none of the other metal(loid)s, but had a significant correlation with pH, indicating that it has a unique source and is related to pH.

Statistical Analysis
Correlation analysis (CA) and PCA can reflect the correlation among metal(loid)s levels to some extent and distinguish the different groups of metal(loid)s [44], which was used to reveal the association of the source of metal(loid)s in the study area. A strong correlation can reveal the close relationship between elements, indicating that these elements may have common sources and similar migration pathways [35,49], which can provide a reference for the traceability of metal(loid)s. The relationships among Cr, Mn, Zn, Pb, Cd, As, and pH were examined by Pearson's correlation analysis, as shown in Table 2. The correlations among Pb, Zn and Cd were highly significant at the p < 0.01 level (r > 0.77). Similarly, Mn with Cd and Zn were strongly correlated at the p < 0.01 level (r > 0.53), indicating that these elements might share a common pollution source. And the correlation between As and Cd was significant at the p < 0.05 level. In particular, Cr was correlated with none of the other metal(loid)s, but had a significant correlation with pH, indicating that it has a unique source and is related to pH. Prior to PCA, KMO and Bartlett's spherical tests (KMO = 0.679, Sig = 0.000 < 0.001) were performed on the data to determine their suitability for PCA. Three principal components (PC1, PC2, and PC3) (eigenvalues > 1) were obtained from the data for the metal(loid)s in the soil, representing almost 81.04% of the total variance of explanatory variables, with 47.01 % of the variance explained by PC1, 19.10 % by PC2, and 14.93 % by PC3 (Table 3; Figure 4). Loading of Cd, Zn, Pb, and Mn were shown by PC1, that of Cr by PC2, and that of As by PC3. According to their strong correlations in PC loadings, six elements could be classified into Group 1 (Cd, Zn, Pb, and Mn), Group 2 (Cr), and Group 3 (As) ( Table 3; Figure 4). Prior to PCA, KMO and Bartlett's spherical tests (KMO = 0.679, Sig = 0.000 < 0.001) were performed on the data to determine their suitability for PCA. Three principal components (PC1, PC2, and PC3) (eigenvalues > 1) were obtained from the data for the metal(loid)s in the soil, representing almost 81.04% of the total variance of explanatory variables, with 47.01 % of the variance explained by PC1, 19.10 % by PC2, and 14.93 % by PC3 (Table 3; Figure 4). Loading of Cd, Zn, Pb, and Mn were shown by PC1, that of Cr by PC2, and that of As by PC3. According to their strong correlations in PC loadings, six elements could be classified into Group 1 (Cd, Zn, Pb, and Mn), Group 2 (Cr), and Group 3 (As) (Table 3; Figure 4).

PMF Modeling
The PMF model were utilized to better identify the sources of metal(loid)s in this study area and quantify their corresponding contribution of sources to each element. To ensure high accuracy of program results, the Q value was minimized and different factors (2, 3, and 4) were selected to control the residual matrix E, and the base runs were performed 20 times to find the most ideal solution. When the number of factors was 3, the Q value was the smallest (Q = 0.5), and the residuals of each element were between −3 and 3. The R 2 values for each element were as follows: Cr (0.97), Mn (0.61), Zn (0.62), Pb (0.74), Cd (0.63), and As (0.99). The R 2 of all six elements were greater than 0.6, suggesting high

PMF Modeling
The PMF model were utilized to better identify the sources of metal(loid)s in this study area and quantify their corresponding contribution of sources to each element. To ensure high accuracy of program results, the Q value was minimized and different factors (2, 3, and 4) were selected to control the residual matrix E, and the base runs were performed 20 times to find the most ideal solution. When the number of factors was 3, the Q value was the smallest (Q = 0.5), and the residuals of each element were between −3 and 3. The R 2 values for each element were as follows: Cr (0.97), Mn (0.61), Zn (0.62), Pb (0.74), Cd (0.63), and As (0.99). The R 2 of all six elements were greater than 0.6, suggesting high correlations, therefore, the positive matrix factorizing model (PMF) was applicable, and the results could be considered reliable.
The operation results of PMF 5.0 were shown in Figure 5, which was consistent with the grouping result of PCA. The first factor (Factor 1) was mainly dominated by Cr, accounted for 22.95% of the total variance and Cr had the highest concentrations (its contribution was 87.0%) in Factor 1. Factor 2 was characterized by Cd, Zn, Pb, and Mn (their contributions were 82.6%, 43.7%, 68.3%, and 82.2%, respectively), accounting for 48.68% of the total variance, which was the largest among the three factors. Factor 3 was mainly dominated by As (the contribution was 83.8%), which accounted for 28.37% of the total variance. The operation results of PMF 5.0 were shown in Figure 5, which was consistent with the grouping result of PCA. The first factor (Factor 1) was mainly dominated by Cr, accounted for 22.95% of the total variance and Cr had the highest concentrations (its contribution was 87.0%) in Factor 1. Factor 2 was characterized by Cd, Zn, Pb, and Mn (their contributions were 82.6%, 43.7%, 68.3%, and 82.2%, respectively), accounting for 48.68% of the total variance, which was the largest among the three factors. Factor 3 was mainly dominated by As (the contribution was 83.8%), which accounted for 28.37% of the total variance.

Cancer Risk and Hazard Quotient
Based on the guideline (HJ 25.3-2019) [39], the study area belonged to category II industrial land, and its human health risk assessment only considered adults. Three elements, Cr, Cd, and As, were associated with carcinogenic risk, and all six elements were associated with non-carcinogenic risk. However, the method we selected (Chinese guideline and USEPA model) did not estimate the relevant parameters of Pb and was not fit for the risk assessment of Pb. Therefore, only the other five elements were evaluated. Table 4 and Figure 6 show

Cancer Risk and Hazard Quotient
Based on the guideline (HJ 25.3-2019) [39], the study area belonged to category II industrial land, and its human health risk assessment only considered adults. Three elements, Cr, Cd, and As, were associated with carcinogenic risk, and all six elements were associated with non-carcinogenic risk. However, the method we selected (Chinese guideline and USEPA model) did not estimate the relevant parameters of Pb and was not fit for the risk assessment of Pb. Therefore, only the other five elements were evaluated. Table 4 and Figure 6 show the adult CR and HQ results under different exposure routes. CR of Cr ranged from 1.15 × 10 −5 to 2.97 × 10 −3 (average = 5.40 × 10 −4 ), indicating an unacceptable carcinogenic risk; the CR of Cd ranged from 5.25 × 10 −10 to 3.46 × 10 −7 (average = 6.75 × 10 −8 ), indicating no obvious carcinogenic risk of Cd; the CR of As ranged from 6.00 × 10 −6 to 4.33 × 10 −5 (average = 1.66 × 10 −5 ), indicating an existing but acceptable carcinogenic risk. Overall, the average total carcinogenic risks of Cr, Cd, and As (CR n = 1.85 × 10 −4 ) suggested that the smelting area posed an unacceptable carcinogenic risk of the surrounding populations. Cr consists of the highest carcinogenic risk contribution (97.30%), most of the computed carcinogenic values for Cr were greatly more than the threshold limit of 1 × 10 −4 , whereas As and Cd rarely exceeded the limit. Based on the exposure pathways, the relationship among the carcinogenic risk of the three exposure pathways was CR pis > CR dcs > CR ois (Figure 6 and Table 4), and the inhalation route was the main source of carcinogenic risk.

Pollution Sources
The CA, PCA, and PMF results were combined with the concentration and spatial distribution of metal(loid)s to identify the source of each element. Factor 1 was mainly defined by Cr. Generally, Cr in the soil was considered to be from a natural source and significantly influenced by the parent material. Many previous studies on industrial areas have also obtained this point of view [14,50,51]. However, the concentration of Cr in the The ranges of HI value for Cr, Mn, Zn, Cd, and As were from 6.28 × 10 −2 to 1.62 × 10 1 , 1.47 × 10 −3 to 1.25 × 10, 6.69 × 10 −4 to 1.67 × 10 −2 , 2.32 × 10 −4 to 1.53 × 10 −1 , and 1.24 × 10 −1 to 8.95 × 10 −1 , respectively. Adults under different exposure pathways showed that the order of average non-carcinogenic risk (HI) of metal(loid)s was Cr (2.93 × 10) > As (3.43 × 10 −1 ) > Mn (2.16 × 10 −1 ) > Cd (2.98 × 10 −2 ) > Zn (2.92 × 10 −3 ). The average values of HI for four metal(loid)s were less than 1 × 10, indicating that these elements did not pose a non-carcinogenic health risk, except Cr. Additionally, the contribution of the three exposure routes to total HI was HQ ois > HQ pis > HQ dcs , and the ingestion route was the main source of non-carcinogenic risk.

Pollution Sources
The CA, PCA, and PMF results were combined with the concentration and spatial distribution of metal(loid)s to identify the source of each element. Factor 1 was mainly defined by Cr. Generally, Cr in the soil was considered to be from a natural source and significantly influenced by the parent material. Many previous studies on industrial areas have also obtained this point of view [14,50,51]. However, the concentration of Cr in the smelter we studied was far beyond the supply of natural sources. According to the previous descriptive analysis, the concentration of Cr was much higher than the background value (more than 28.81 times) and showed exceptionally high variability. Thus, Factor 1 could be excluded from the soil-forming parent material source. Based on the spatial distribution results, Cr was enriched near the chromium slag yards. Due to a historical reason, these two large chromium slag yards were formed near the NO. 3 smelting workshop to stack the waste slag produced using the process of chromium-related products, with 200,000 tons of chromium slag accumulated. This chromium slag infiltrated into deep soil and groundwater through leaching and was gradually enriched in the soil. Additionally, according to the CA results, the source of Cr had a certain correlation with pH. Some studies have also shown that pH is an essential factor in the migration of Cr in the soil [52]. The pH level can affect the chemical transformation of Cr in the soil and change its occurrence form, subsequently impacting its migration capacity in the soil. For example, carbonate bound chromium in the soil can easily be transformed into watersoluble chromium at low pH, increasing its migration ability in the soil environment to increase the migration range and depth of Cr [52]. Therefore, it could be inferred that Cr pollution was mainly caused by the centralized stacking of chromium slag; thus, Factor 1 could be considered due to the stacking of waste slag.
Factor 2 was characterized by Cd, Zn, Pb, and Mn. According to the spatial distribution, the contents of Cd, Pb, and Zn near the blast furnace smelting workshop and the NO. 3 smelting workshop (Section 304) were much higher than their background values. The blast furnace smelting workshop mainly smelted pig iron, manganese residual, and other Fe-Mn products, while Section 304 was also mainly responsible for the smelting process. Some previous studies have several Pb, Zn, and Cd in the ash produced in the material preparation and purification of the smelting process; these ashes can increase the content of Pb, Zn, and Cd in the topsoil [53]. It is generally considered that the main source of Pb, Zn, and Cd is industrial smelting [8]. Although many previous studies reported that Mn in the soil might be controlled by the soil parent material [50,54,55], in the smelting site of our study, manganese ore, as an important raw material, easily accumulated in the soil during the smelting process. In the process of iron tapping, a large amount of dust removal ash and Mn-rich slag were produced. This fly ash and waste residue could enter the soil through atmospheric sedimentation, direct discharge and infiltration [56], resulting in the enrichment of Mn. Obviously, smelting activities were a key pollution source of these four elements in the soil across the study area. Therefore, Factor 2 could be ascribed to the source smelting activities.
Factor 3 was mainly dominated by As. Many previous studies attributed the source of As to agricultural input [36,57], and considered that the use of pesticides and chemical fertilizers leads to the enrichment of As [36]. However, the study area has been used as industrial land for more than 60 years, and most of the surrounding land was industrial land. Therefore, we believe that As in the study area was less likely to come from agricultural input. According to the previous analysis, the As content generally fluctuated around the soil background value of Hunan Province, and it was slightly higher than this background value. In addition, CV was 47.01%, and it differed from other exceptionally high elements, indicating that it was less impacted by human activities in the study area. Thus, we tend to think that As comes from natural input. Some previous studies have also confirmed that As in the soil may come from soil parent material [46,58]. Therefore, Factor 3 could be considered to be affected by a natural source.

Human Health Risk
According to the International Agency for Research on Cancer guidelines [59], Cr poses a significant carcinogenic risk [32,59]. The risk assessment results in this study also indicate the leading contribution of Cr in carcinogenic risk. Therefore, in this smelting site, significant attention should be paid to the risk and migration changes of Cr. Among the three exposure routes, inhalation by respiration is the most important exposure route, indicating that the downwind direction of the smelter might has higher carcinogenic risks than the other areas. Therefore, when entering or passing near the site, people should wear dust masks to prevent inhalation of soil dust particles. Additionally, during subsequent urban construction, residential areas, parks, and other densely populated areas should not be built downwind of the smelting site during subsequent urban construction.
Among the non-carcinogenic risks, only the average hazard quotient of Cr exceeds the threshold, showing an unacceptable non-carcinogenic risk as a whole. However, we found that although other elements have no significant non-carcinogenic risk as a whole, but the HI values of Mn at some sites exceeded the threshold. It showed that, while overall, the region did not present a significant non-carcinogenic risk, some sites might show considerable non-carcinogenic risks for humans, such as slag yard stacking areas with highly concentrated pollutants and waste residues. Therefore, more attention should be paid to sampling sites that exceeded the acceptable hazard quotient in the study area. For instance, for Cr, in 45% of the sites, the HI value exceeded 1, and the maximum HI value reached 16.17; thus, requiring significant attention. Additionally, the risk assessment results showed that oral intake was the main exposure route of non-carcinogenic risk. However, our study area was an industrial land, and the main people concerned were adults. Previous studies also showed that geophagism (soil ingestion) was rarely observed in few adults [32]; thus, avoiding agricultural development in the study area can inhibit the risks brought by the ingestion pathway. Therefore, we believe that the inhalation route needs more attention.

Conclusions
The present study evaluated an industrial smelting site in South Central China for contamination by metal(loid)s, their source, and associated human health risks. Although the smelting site had been dormant for several years, soil samples with high concentrations of metal(loid)s were still detected. The average concentrations of Cr, Mn, Zn, Pb, and Cd were much higher than the background value of the Hunan Province. As evidenced by the I geo and PLI analysis, the pollution order was Cr > Cd > Mn > Pb > Zn > As. Among the six elements, the most heavily polluted by Cr and Cd, whereas it was only slightly polluted by As. The study area was highly contaminated by metal(loid)s as a whole, especially near the smelting workshop and slag yards. The statistical analysis revealed three primary sources of elements, with the main source of Mn, Zn, Pb, and Cd being smelting activities. In contrast, Cr was attributed to the stacking of chromium slag, and As was related to a natural source. Among them, smelting activities constituted the main pollution source with the highest contribution to the human health risk in the study area. Moreover, human health risk assessment indicated that there was no severe non-carcinogenic risk at the time of evaluation, but the carcinogenic risk was unacceptable across the entire study area. Cr was the major element affecting CR and HQ, which significantly contributed to the noncarcinogenic and carcinogenic risks for local residents. Therefore, priority should be given to the remediation of Cr contamination. The existing carcinogenic and non-carcinogenic risks of Cr should be promptly addressed and require appropriate control and protection strategies. Additionally, environmental management and control must find measures to cut off the main receptor exposure pathways-namely, inhalation, and targeted pollution source control is urgently necessary.