Implications of Soil Potentially Toxic Elements Contamination, Distribution and Health Risk at Hunan’s Xikuangshan Mine

: A ﬁeld survey was conducted to determine the pollution grade, sources, potential ecological risk, and health risk of soil potentially toxic elements (PTEs) in Xikuangshan Mine (XKS), the largest antimony (Sb) deposit in the world. A total of 106 topsoil samples were collected from 6 sites in XKS to measure the concentrations of PTEs Cr, Zn, Cd, Pb, As, Hg, and Sb. The results show that the average concentrations of these elements at all six sites were generally greater than their corresponding background values in Hunan province, especially Sb, Hg, and As. Correlation and principal component analyses suggested that Cd, Zn, Pb, Hg, and Sb were primarily released from mining and other industrial and human activities, while Cr and As were mainly impacted by the parent material from pedogenesis. A risk index analysis showed that, overall, sites were at very high ecological risk, and Sb is the highest ecological risk factor, followed by Cd and Hg. According to health risk assessment, oral ingestion is the main non-carcinogenic and carcinogenic risk exposure route. The higher potentially non-carcinogenic and carcinogenic risks happen to the local children who live in the vicinity of mining area. It revealed that the mining and smelting processes of XKS have negatively inﬂuenced the local people, therefore, we should pay increasing attention to this practical issue and take effective measures to protect the ecology of XKS.


Introduction
Potentially Toxic Elements (PTEs) are considered as toxic and persistent pollutants in soil which remain in local ecosystems for a long time and pose a negative threat to human health [1,2]. When PTEs in soil are transformed from solid to ionic form and methylated, they will accumulate in organisms, threatening the safety of animals and humans through the food chain [3].
Xikuangshan Mine (XKS) is well known as the World Antimony Capital and has been developed into the largest industrial and population centre in Hunan province over the past hundreds of years, especially the 21st century [4]. Approximately 3.20 billion cubic meters of Sb-containing exhaust gas has been dispersed into the atmosphere in recent years [5]. The soil Sb contents near XKS ranged from 100 to 5045 mg kg −1 [6]. As a consequence, the soil in the XKS has been extensively contaminated with high concentrations of Sb and other associated metals due to the rapid expansion of the mining industry and urbanization [7]. Sb and its compounds are listed as priority pollutants by the U.S. Environmental Protection Agency (USEPA). High levels of toxic metals may pose health risks to humans.
Due to historical reasons, the century-old massive mining operations in the XKS area were directly discharged without treatment and waste residues were deposited at will,

Sample Analysis
The soil samples were naturally air-dried, after removing the shell fragments and grounded with an agate mortar, the solid was passed through the 20-mesh and 100-mesh nylon mesh sieve. Soil pH was measured using a soil/water ratio of 1:2.5. Organic carbon content was determined by the dichromate oxidation method. Soil samples were digested using a concentrated acid mixture (HNO3-H2O2-HF) for PTE analysis. Briefly, approximately 0.5 g of soil were digested with 8 mL 65% HNO3, 5 mL 30% H2O2, and 3 mL 30% HF in Teflon vessels. The vessels were sealed and incubated in microwave system at 150 °C for 10min, at 180 °C for 10 min, and then at 200 °C for 15 min until a clear liquid was obtained. After digestion, samples were evaporated near to dryness in a sand bath at 120 °C for 3 h. After cooling below 100 °C, digested samples were made up to 10 mL volume with Milli-Q water [20,21]. The extract was then stored at 4 °C until it was analyzed. The

Sample Analysis
The soil samples were naturally air-dried, after removing the shell fragments and grounded with an agate mortar, the solid was passed through the 20-mesh and 100-mesh nylon mesh sieve. Soil pH was measured using a soil/water ratio of 1:2.5. Organic carbon content was determined by the dichromate oxidation method. Soil samples were digested using a concentrated acid mixture (HNO 3 -H 2 O 2 -HF) for PTE analysis. Briefly, approximately 0.5 g of soil were digested with 8 mL 65% HNO 3 , 5 mL 30% H 2 O 2 , and 3 mL 30% HF in Teflon vessels. The vessels were sealed and incubated in microwave system at 150 • C for 10min, at 180 • C for 10 min, and then at 200 • C for 15 min until a clear liquid was obtained. After digestion, samples were evaporated near to dryness in a sand bath at 120 • C for 3 h. After cooling below 100 • C, digested samples were made up to 10 mL volume with Milli-Q water [20,21]. The extract was then stored at 4 • C until it was analyzed. The concentrations of Cr, Zn, Cd, and Pb were measured with inductively coupled plasma atomic emission spectrometry (Thermo Scientific, iCAP Qc). The concentrations of As, Hg, and Sb were determined using a dual-channel atomic fluorescence photometer (Beijing Haiguang instrument Co., AFS-2100). For quality control (QC), standard reference soils (GBW07407; Center for Certified Reference Materials, China) and blank controls were The single-factor pollution index (P i ) and Nemerow synthetic pollution index (P N ) were calculated to assess the degree of metal pollution in the soils [22]. The single-factor pollution index was calculated as follows: where P i is the single-factor pollution index of element i in soil; C i is the on-the-spot survey concentration of element i; and C o is the element i background value of Hunan soil (mg kg −1 ) (Table S1). P i > 1 represents element pollution. The Nemerow synthetic pollution index was calculated using the following formula: where P N is the Nemerow synthetic pollution index for all assessed samples; P MAX is the maximum of the single-factor pollution indices for all samples; and P AVE represents the arithmetic average of the single-factor pollution indices for all samples. PTEs pollution was classified into five classes or grades based on the Nemerow index (Table S2).

Potential Ecological Risk Assessment
Hakanson potential ecological hazard index method was used to assess the ecological risk of soil in the six mining areas [23]. The method is the current soil ecological risk assessment of PTEs, not only considering the soil element content, also being linked to the ecological environment effect and toxicology, and the difference of regional background values of PTEs. The calculation formula is as follows: where E i is the potential risk index of element i; Ti is the toxicity response parameter of element i (the toxicity response parameters of Cr, Zn, Cd, Pb, As, Hg, and Sb are 2, 1, 30, 5, 10, 40, and 7, respectively) [23,24]; P i is the single-factor pollution index of i. RI is the comprehensive potential ecological risk index of various PTEs. Potential ecological risk was classified into five grades based on the E i and RI index (Table S3).

Calculation of PTEs Intake
Cr, Zn, Cd, Pb, As, Hg, and Sb all have chronic non-carcinogenic health risks, while Cr, Cd, Pb, and As also exhibit carcinogenic risks. PTEs in the soil can enter the human body through three exposure routes: oral ingestion, dermal absorption, and inhalation. The daily average intake from the three exposure routes were calculated using ADI oral ingestion, ADI dermal absorption, and ADI inhalation (mg/kg/d) models, shown as follows: • Oral ingestion: Processes 2021, 9, 1532 5 of 18 where ADI o is the average daily intake of PTEs from oral ingestion (mg/kg/day), C i is the element concentration found in the soil (mg/kg), IR is the ingestion rate (kg/day), CF is the conversion factor (kg/mg), EF is the exposure frequency (day/year), ED is the exposure duration (year), BW is the body weight of the exposed individual (kg), AT is the time period over which the dose is averaged (day), and SIR is the ingestion rate of soil (mg/day).
• Dermal absorption: where ADI D is the average daily intake of PTEs from dermal absorption (mg/kg/day), SA is the exposed skin surface area (cm 2 ), AF is the adherence factor (mg/cm 2 /day), and ABS is the dermal absorption factor (dimensionless).
• Inhalation: where ADI R is the average daily intake of PTEs from inhalation (mg/kg/day), PM 10 is the content of inhalable suspended particulate matter in the air (mg/cm 3 ), DAIR is the human daily air respiration (m 3 /day), and PIAF is the retention ratio of inhaled particles in human body (dimensionless). FSPO is the proportion of particles from soil in the air (dimensionless). The definitions and values of other parameters proposed and revised by the USEPA are shown in Table S4 [25,26].

Human Health Risk Assessment
The human health risk assessment results of PTEs were divided into two types: one is the non-carcinogenic effect with index HQ i from a single element (i) and HI from multiple PTEs; the other is carcinogenic effect, expressed in terms of carcinogenic risk (CR). The specific calculation model formula is as follows: where i is a non-carcinogenic or carcinogenic element; j is a certain exposure route; ADI ij is the average daily exposure of PTEs i under the exposure route in (mg/kg/day). RfD is the reference dose (mg/kg/day); SF is the carcinogenic risk slope factor under different exposure routes (per mg/kg/day); and the reference dose values and slope factor reference values for various exposure routes are shown in Table S5.
For multiple hazardous substances, if the value of HI < 1, it is unlikely to create adverse health effects for exposed populations. If the value of HI > 1, it is meaning a potential risk associated with the contaminant [27]. CRI < 1 × 10 −6 specified in the "Technical Guidelines for Risk Assessment of Contaminated Sites" issued by the Ministry of Environmental Protection of China [28], indicating that there is no negative impact on the health of the population. CRI lying between 1 × 10 −4 and 1 × 10 −6 are generally considered an acceptable range, depending on the situation and circumstances of exposure. CRI > 1 × 10 −4 are viewed as unacceptable.

Statistical Analysis
Several software were employed in this study, including Arc GIS for preparing the sample distribution map and spatial interpolation mapping, SPSS (IBM SPSS Statistics 23) for correlation analysis among multiple variables, Origin (Origin 2019) for boxplot analysis of PTEs content distribution, and principal component analysis among multiple variables for each sample point.

Soil Properties and PTEs Concentrations
The characteristics of the soils sampled at the six sites are presented in Table 1. The pH values of the soils ranged from 4.29 to 5.33. The soils at YJT, SZS, and SYA were acidic. The low soil pH at these three sites can be attributed to decades of irrigation with waste water derived from mining. The average OC content ranged from 26.41 to 38.67 g kg -1 , with the highest and lowest levels found in the SZS and SXS soils, respectively. As shown in Figure 2, the average concentrations of the PTEs in the topsoil samples (0-20 cm) from the six sites were generally greater than their corresponding background levels in Hunan (indicated by the black lines in Figure 2), especially Sb, Hg, and Cd, which exceeded the standard rate by 100%, 100%, and 97.17%, respectively. In particular, the Sb concentrations at YJT, SZS, SYA, SXN, SXS, and LWC were 75. 13, 825.33, 199.57, 901.57, 567.45, and 154.73 times higher than the background value (2.98 mg kg −1 ), respectively ( Figure 2f). Greatest concentration of Sb was observed in tailings from SXS (8591.16 mg kg −1 ), which was close to the previously reported concentration [11]. The relatively high Sb concentrations were found in the sampling sites rather than the paddy soil in the vicinity of the northern mine (5.91−322.35 mg kg −1 ) [12]. Compared to the soil which was indirectly exposed to pollution sources, such as the areas near the mines, smelters and tailings, the mean concentration of Sb in Macleay catchment coastal floodplain (located in northeastern New South Wales) was 11.2 mg/kg [29]. Mining and smelting activities produce residues that frequently contain high levels of Sb in the form of sulfides, such as stibnite (Sb 2 S 3 ), which are readily oxidized when exposed to the environment [30]. This is consistent with the remarkably high level of Sb pollution observed in the smelter sites of SZS, SXN, and SXS. Moreover, it is well-known that the common sources of Hg are fossil fuel combustion (especially coal burning), industrial fumes, vehicle exhausts, and sewage sludge [31]. In our study, the three sites (SZS, SXN, and SXS) showing high Hg levels are located downwind of XKS, suggesting that coal-burning exhausts, industrial fumes, and vehicle exhausts may contribute to the increasing Hg contamination.
In contrast, the average concentrations of Cr at the six sites were much lower. In the SZS and LWC sites, the concentrations of Cr in soil were only 98.31% and 98.36% of the background value (67 mg kg −1 ), respectively. Moreover, Cr displayed low CV values (Table S1) and quite homogeneous distributions across the entire study area, indicating that the Cr is mainly affected by pedogenesis while Sb and Hg are heavily influenced by anthropogenic activities. Overall, the results indicated that the concentration of PTEs, especially Sb, Hg, and As in the soil samples of XKS is relatively serious and requires adequate attention. The environmental behavior of PTE is critically dependent on their chemical form, which influences their mobility, bioavailability, and toxicity to biota. Determination of the total concentrations of metals in soils is not enough to predict the hazard of these elements [32]. It reported that despite the higher Sb concentration in XKS soil, rice accumulated more As than Sb, indicating the higher phytoavailability of As [12]. We caution that the exchangeable fraction of Sb, Hg, and As should be verified, however, we have obtained the element with high risk at present, and the results lay the foundation for the next experiments. the Cr is mainly affected by pedogenesis while Sb and Hg are heavily influenced by anthropogenic activities. Overall, the results indicated that the concentration of PTEs, especially Sb, Hg, and As in the soil samples of XKS is relatively serious and requires adequate attention. The environmental behavior of PTE is critically dependent on their chemical form, which influences their mobility, bioavailability, and toxicity to biota. Determination of the total concentrations of metals in soils is not enough to predict the hazard of these elements [32]. It reported that despite the higher Sb concentration in XKS soil, rice accumulated more As than Sb, indicating the higher phytoavailability of As [12]. We caution that the exchangeable fraction of Sb, Hg, and As should be verified, however, we have obtained the element with high risk at present, and the results lay the foundation for the next experiments.  (Table S1).

Pollution Index Evaluation of PTEs
The pollution indices of the metals in the soils ( Figure 2) and the pollution grades (Table S2) were calculated. Pi varied greatly from PTEs and sites. Only four Pi (14.69% of the entire sample set) were <1, while the majority were >1. The Pi of Sb had a maximum value of each site. Sb pollution was the most serious element in soil, consistent with the finding of the wastewater, dust, and tailings produced from discharge points at the antimony mineral processing plant of XKS [14]. For all the metals at SXN and SXS, Pi > 1, and it was also found SXN had the highest PSb (902.57 ± 910.98, average ± standard deviation, hereafter the same), indicating an extensive Cr/Zn/Cd/Pb/As/Hg/Sb multi-pollution at  (Table S1).

Pollution Index Evaluation of PTEs
The pollution indices of the metals in the soils ( Figure 2) and the pollution grades (Table S2) were calculated. P i varied greatly from PTEs and sites. Only four P i (14.69% of the entire sample set) were <1, while the majority were >1. The P i of Sb had a maximum value of each site. Sb pollution was the most serious element in soil, consistent with the finding of the wastewater, dust, and tailings produced from discharge points at the antimony mineral processing plant of XKS [14]. For all the metals at SXN and SXS, P i > 1, and it was also found SXN had the highest P Sb (902.57 ± 910.98, average ± standard deviation, hereafter the same), indicating an extensive Cr/Zn/Cd/Pb/As/Hg/Sb multi-pollution at SXN  six sampling sites and a high potential risk from the ecological environment, especially at SNS, SXN, and SZS.

Correlation Analysis of PTEs in Soil
Inter-element relationships can provide helpful information about metal sources and pathways. Pearson correlation coefficients were calculated between seven PTEs at the six sampling sites, as shown in Table 2. The high correlations suggested that the indicated contaminants have originated from a common source [33]. A significantly positive correlation (p < 0.01) was found with the following elemental pairs: Zn-Cd (0.405), Zn-Pb (0.985), and Cd-Pb (0.384). This indicates the possible same sources of these elements. It is in agreement with other authors who also found that Zn exhibited positive and significant correlations with Pb and Cd, indicating their common origins and similar pathways in XKS [34]. On the other hand, a highly positive correlation was found between Cd and Hg (0.445), Cd and Sb (0.344), and Hg and Sb (0.545) at a 0.01 level, which implies these PTEs were supplied by similar pollution sources. In addition, it was also found that Cr is only positively correlated with As at p < 0.01 and negatively correlated with Hg and Sb at p < 0.05.

Principal Component Analysis of PTEs in Soil
Principal component analysis (PCA) was used to identify the original sources of the seven PTEs by evaluating their correlation [35]. Finally, the initial dimension dataset was narrowed down into two major components PC1 and PC2 (Figure 3), which can explain 67.1% of the data variation. PC1 explained 47.0% of the total variance and loaded heavily on Sb, Cd, Pb, Hg, and Zn; meanwhile, PC2 accounted for 20.1% and was dominated by Cr and As. PC1 included Cd, Zn, Pb, Hg, and Sb, which accounted for most of the variability. The load of Cd was up to 0.909, and its correlation with Zn, Pb, Hg, and Sb was extremely significant (p < 0.01, Table 2), which suggests a shared pollution source for the five PTEs. Compared with the background values of Hunan soil, the remarkably elevated concentrations of Cd, Zn, Pb, Hg, and Sb in soils of XKS indicated that anthropogenic sources are likely to be the pollution origins of these metals (Figure 1). Researchers found that smelting slag is an important source of Sb pollution in nearby farmland [36]. PC1 could be defined as an anthropogenic component related to mining and industry practices. The antimony ore mining and smelting industry activities in the study area producing wastewater, waste gas, waste residue are the main sources of Sb, Cd, and Hg [37]. Mining and industrial activities and other human factors are affecting the concentrations of Sb, Cd, Pb, Hg, Zn in the soil.
PC2 included Cr and As. PC2 correlated strongly with Cr, which has the highest loading value (0.909). This correlation was confirmed by the contamination analysis ( Table 2). The coefficient of variation of soil Cr showed non-point source contamination, suggesting that there is no significant anthropogenic input of Cr. Cr concentrations in soils did not decrease with depth, which was derived from parent materials [38]. Moreover, Cr is strongly associated with the aluminosilicate phases and Fe-oxide species of the soils, which are important products of parent rock weathering [39]. A significant correlation was found between Cr and As in YJT (Table S6), while no relationship was found in SZS, SXN, SXS, and LWC. PC2 could be a lithogenic component, as the variabilities of soil Cr and As were controlled by parent material and pedogenic processes. The metal concentrations and their correlations reflect the strong impact on parent materials in the pollution of the surrounding soils. Since As exhibited a positive correlation with Cr in XKS (Table 2), and widespread elevated As values were found in all sites except YJT (Figure 1), its sources may be complex, such as from various emissions, e.g., agricultural practices or industrial contamination [40]. This could be attributed to anthropogenic inputs of Cd, Zn, Pb, Hg, and Sb which disturbed the original correlation with Cr and As.
Processes 2021, 9, x FOR PEER REVIEW 9 of 19 PC2 included Cr and As. PC2 correlated strongly with Cr, which has the highest loading value (0.909). This correlation was confirmed by the contamination analysis (Table 2). The coefficient of variation of soil Cr showed non-point source contamination, suggesting that there is no significant anthropogenic input of Cr. Cr concentrations in soils did not decrease with depth, which was derived from parent materials [38]. Moreover, Cr is strongly associated with the aluminosilicate phases and Fe-oxide species of the soils, which are important products of parent rock weathering [39]. A significant correlation was found between Cr and As in YJT (Table S6), while no relationship was found in SZS, SXN, SXS, and LWC. PC2 could be a lithogenic component, as the variabilities of soil Cr and As were controlled by parent material and pedogenic processes. The metal concentrations and their correlations reflect the strong impact on parent materials in the pollution of the surrounding soils. Since As exhibited a positive correlation with Cr in XKS (Table  2), and widespread elevated As values were found in all sites except YJT (Figure 1), its sources may be complex, such as from various emissions, e.g., agricultural practices or industrial contamination [40]. This could be attributed to anthropogenic inputs of Cd, Zn, Pb, Hg, and Sb which disturbed the original correlation with Cr and As.

Evaluation of Potential Ecological Risks of PTEs
To estimate the potential ecological risk posed by PTEs in the soils from sampling sites, we calculated the potential ecological risk factor (Ei) of each individual element and the comprehensive index (RI) of the seven PTEs, as shown in Table 3 and Figure 4.
According to Hakanson's potential ecological risk classification standard (Table S3), the result indicated that the pollution of the evaluated metals was in the following order: Sb > Cd > Hg > As > Pb > Cr > Zn. This result was in agreement with the study of Tang [15]. For each individual element, Sb, Cd, and Hg posed very high ecological risk, and there was low risk for As, Pb, Cr, and Zn. The Ei of Cr, Zn, Pb, and As in six sites were all

Evaluation of Potential Ecological Risks of PTEs
To estimate the potential ecological risk posed by PTEs in the soils from sampling sites, we calculated the potential ecological risk factor (E i ) of each individual element and the comprehensive index (RI) of the seven PTEs, as shown in Table 3 and Figure 4.
According to Hakanson's potential ecological risk classification standard (Table S3), the result indicated that the pollution of the evaluated metals was in the following order: Sb > Cd > Hg > As > Pb > Cr > Zn. This result was in agreement with the study of Tang [15]. For each individual element, Sb, Cd, and Hg posed very high ecological risk, and there was low risk for As, Pb, Cr, and Zn. The Ei of Cr, Zn, Pb, and As in six sites were all below 40, placing these metals at low ecological risk levels. For Hg, the average E i in YJT and LWS was 172.12 and 257.96, respectively, which places it at a moderate pollution level. Meanwhile, the average E i in SXS and LWC for Cd reached 2505.19 and 1615.96, respectively, which is considered as a very high ecological risk. As the metal with the most potential risk, Sb posed the highest ecological risk in SXN, SZS, and SXS, reaching 6317.99, 5784.30, and 3979.18, respectively, followed by SYA (1403.32). The E Sb of soils surrounding the tailings reservoir in the north of XKS reached 10548.42, indicating very high ecological risk [14]. Based on the average RI values (Figure 4) of all seven PTEs, the analysis showed that, overall, sites were at very high ecological risk. Further calculations showed that 72.64% of all soil samples were at very high ecological risk. Spatially, the pollution level at the six sampling sites were ordered as SXN > SXS > SZS > LWC > SYA > YJT. The results show that due to the long-term mining, beneficiation, and production activities of SXS and SXN enterprises, the ecological risk of PTEs in the surrounding soil is high. To conclude, all sites are experiencing high ecological risk, with Sb presenting the highest risk and Zn the least. The values shown are average ± standard deviation. Different letters in the same column denote significant differences (p < 0.05).

Health Risk Assessment of PTEs for Adults and Children
 Non-carcinogenic risk for adults and children These results for non-carcinogenic risk for adults and children via oral ingestion, dermal absorption, and inhalation pathways are all presented in terms of HQi, as shown in Tables 4-6 and the total non-carcinogenic hazard index (HI) in Figure 5. RfD values are presented in Table S5 and ADI values in Tables S7 to S9. When HQi and HI values exceed 1, this indicates a concern for the potential non-

Health Risk Assessment of PTEs for Adults and Children
• Non-carcinogenic risk for adults and children These results for non-carcinogenic risk for adults and children via oral ingestion, dermal absorption, and inhalation pathways are all presented in terms of HQ i , as shown in Tables 4-6 and the total non-carcinogenic hazard index (HI) in Figure 5. RfD values are presented in Table S5 and ADI values in Tables S7-S9. When HQ i and HI values exceed 1, this indicates a concern for the potential noncarcinogenic effects. The numerical results were higher than the research of heavy metal in agricultural surface soil for the Guangzhou-Foshan urban zone [41]. The HQ of oral ingestion was higher than dermal absorption and inhalation, suggesting that oral ingestion presented a high risk to the health for adults and children. Therefore, for the non-carcinogenic risks, we would first reduce the hazard from the diet. The HQ Sb via oral ingestion was significantly greater than other metals, primarily because of the high concentrations in soils. The average HQ Sb via oral ingestion for children decreased in the following order: Sb > As > Cr > Pb > Hg > Cd > Zn, the same as for adults. Individually, the HQ Sb and HQ As via oral ingestion for children at SZS, SYA, SXN, SXS, and LWC were greater than 1 because of their high concentration or low RfD values. Moreover, the HQ Sb via dermal absorption for adults and children at SZS, SXN, and SXS was greater than 1 and significantly higher than YJT, SYA, and LWC. The ratio of HQ Sb and HQ As (>1) accounted for 8.69% of all the HQ values, indicating that the accumulation of Sb and As would be the main cause of chronic diseases based on their high HQ values. Antimony and its compounds are poisonous materials with acute and chronic toxicity, leading to damage to the liver, skin, lungs, and even the cardiovascular system [42,43]. Thus, because of the hazards posed by these elements, more attention should be paid to this health topic.
The spatial distributions of HI of all interpolated PTEs were similar for both adults and children ( Figure 5). Whereas HI for different populations varies partially, generally in the order of children > adults, the results concur with a previous study [44]. Non-carcinogenic risk values of mining areas in the southern provinces (e.g., Hunan and Guizhou) are higher than other provinces of China, whereas most of the northern and western mining areas' values are less than 1 [45]. The average HI value for the seven metals at all six sites through oral ingestion exposure pathways were greater than 1, which is significantly higher than the safe level. This suggested that oral ingestion was the most significant non-carcinogenic health risk contributor and should be highly stressed.

•
Carcinogenic risk for adults and children Due to the lack of carcinogenic slope factors for Zn, Hg, and Sb, only Cr, Cd, Pb, and As were estimated for their carcinogenic risks. The results of excess lifetime cancer risks were presented in Tables 4-6, and the total cancer risk index (CRI) was shown in Figure 6. More detailed ADI values were listed in Tables S7-S9. The cancer risk was ranked as As > Cr > Cd > Pb, implying that As was the main pollutant source to stimulate cancer. The CR As via oral ingestion pathways were higher than other metals. The CR Cr via dermal absorption and inhalation pathways were higher than other metals. It is worth noting that the CR As and CR Cr via oral ingestion, dermal absorption, and inhalation pathways for adults and children at six sites were greater than 10 −6 , meaning As and Cr are regarded as ultimately causing cancer risks to the population.
The CRI was calculated by summing the individual CR i across all exposure pathways. The results exhibited that the CRI at all six sites through oral ingestion, dermal absorption, and inhalation pathways were totally higher than the acceptable 10 −6 . The acceptable range of CRI for adults and children through dermal absorption and inhalation pathways is between 1 × 10 −4 and 1 × 10 −6 . The ratio of CRI for adults and children >1 × 10 −4 were accounted for 37.26% via oral ingestion, particularly in SXS, indicating that oral ingestion was the severely dangerous carcinogenic health risk factor. A relative research reported that there is a high chance that the non-carcinogenic and carcinogenic effects may occur in smelting site of south mine (SXS site in our research), which needs more attention [46]. Overall, the total cancer risk value of XKS has exceeded the acceptable value, implying great carcinogenic risk. Therefore, the antimony mining areas, smelting plant, and tailings are selected as the priority control mining areas. Management systems and the implementation of reduction projects need to be concentrated in XKS.
Comprehensive analysis of the non-carcinogenic and carcinogenic risks of PTEs in the soil showed that children are more susceptible to PTEs than adults, and the health risks caused by oral ingestion are much higher than dermal absorption and inhalation. This is consistent with a former study that indicates that children have the highest noncarcinogenic and carcinogenic risk in mining areas [47]. It is closely related to the living habits of children, so that the probability of active contact with the soil is much higher than that of adults. Moreover, children have a weak awareness of hygiene, and the rate of ingestion through the hand-mouth way is higher than that of adults [48]. In addition, human sustainably exposure to PTEs in the atmosphere and water have adverse effects including cancer, skin lesions, and peripheral vascular disease, so local governments and the public need to pay more attention and vigilance to ameliorate this urgent situation [49].
The classification standards of PTEs in soil have been clearly defined, but there was no clear specification for other PTE classification standards that prove to have toxic effects on the soil-plants-human body system [50]. It is not accurate to evaluate the environment as a whole only from the level of PTEs pollution in the soil. Therefore, it is recommended to regularly monitor and evaluate the contents and health risks of PTEs in soil and crops in XKS, and the local government should strengthen the remediation of PTEs to help avoid, mitigate and prevent the health hazards of local residents. It is worthy of attention to further improve the soil-plant-human body standard system when conditions permit. Moreover, special attention should be paid to the remediation of polluted soils. It is crucial to take the effective and sustainable method to remediate the PTEs contaminated soils.

Conclusions
(1) The XKS mining areas were highly contaminated. A descriptive statistical analysis determined that the average concentrations of Cr, Zn, Cd, Pb, As, Hg, and Sb at all six sites were generally higher than their corresponding background levels of Hunan. Sb, Hg, and As are relatively more severe and need adequate attention, especially at SZS, SXN, and SXS. (2) A multivariate statistical analysis showed that two main contamination sources were influencing the PTE concentrations. Cd, Zn, Pb, Hg, and Sb were primarily influenced by mining and industrial activities and other human factors, while Cr and As were mainly influenced by the parent material and pedogenic processes. (3) An RI analysis showed that XKS was at a very high ecological risk. Sb was the highest risk factor, followed by Cd and Hg. According to the health risk assessment, the non-

Conclusions
(1) The XKS mining areas were highly contaminated. A descriptive statistical analysis determined that the average concentrations of Cr, Zn, Cd, Pb, As, Hg, and Sb at all six sites were generally higher than their corresponding background levels of Hunan. Sb, Hg, and As are relatively more severe and need adequate attention, especially at SZS, SXN, and SXS. (2) A multivariate statistical analysis showed that two main contamination sources were influencing the PTE concentrations. Cd, Zn, Pb, Hg, and Sb were primarily influenced by mining and industrial activities and other human factors, while Cr and As were mainly influenced by the parent material and pedogenic processes. (3) An RI analysis showed that XKS was at a very high ecological risk. Sb was the highest risk factor, followed by Cd and Hg. According to the health risk assessment, the non-carcinogenic and carcinogenic risks of PTEs in the soil show that children are more susceptible to PTEs than adults. Oral ingestion is the main non-carcinogenic and carcinogenic risk exposure route. The findings of this study found that Sb, As, and Cr pose a certain health risk to local residents, and effective measurements are needed to remediate the toxic metals.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/pr9091532/s1, Table S1: PTEs concentration in the soil of the study area (n = 106), Table S2: Five grade comprising the P N index; Table S3: Grading standards of Hakanson potential ecological risk, Table S4: Values of variables for human health risk assessment models, Table S5: Summary of reference dose (RfD) and cancer slope factor (SF) values for each element from different exposure pathways, Table S6: Correlation coefficients between the PTEs concentrations in the study area, Table S7: Average daily intake (ADI) values for adults and children in soil from XKS mining area for non-carcinogenic risk calculations, Table S8: ADI values for adults and children in soil from XKS mining area for carcinogenic risk calculations.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.