Contamination Levels and the Ecological and Human Health Risks of Potentially Toxic Elements (PTEs) in Soil of Baoshan Area, Southwest China

: The primary goals of this study were to reveal the environmental status of potentially toxic elements (PTEs) and their ecological risks, as well as their associated health risks in the Baoshan area, southwest China, which has been surveyed with the scale of 1:250,000 geochemical mapping. Based on a comparison of the PTE concentrations with the soil environmental quality of China and the enrichment factor (EF), geo-accumulation index (I geo ), contamination factor (C f ), and potential ecological risk indexes ( E ir and PERI), as well as the potential non-carcinogenic hazard indices (HI and CHI) and carcinogenic risks indices (TCR and CTCR), the following conclusions were drawn: The PTE concentrations in the surface soil samples that were collected from the investigated area (1.65% sites) exceeded the risk intervention values (RIV) for soil contamination of agricultural land of China. Cadmium (Cd) and mercury (Hg) posed higher ecological risks than other PTEs (arsenic (As), chromium (Cr), lead (Pb), copper (Cu), nickel (Ni), and zinc (Zn)), which was highlighted by their toxic response factor. Arsenic was the main PTE with a non-carcinogenic risk (19.57% sites for children and 0.25% sites for adults) and the only PTE that carries a carcinogenic risk (2.67% sites for Children and 0.76% sites for adults) to humans in the research area. Children are more vulnerable to health risks when compared to adults because of their behavioral and physiological traits. Geological genesis was responsible for the high concentrations, ecological risk, and health risk distribution patterns of the examined PTEs. Even though the present research highlights several important aspects related to PTE pollution in the research area, further investigations are needed, especially in mining areas.


Introduction
Soil is the foundation of food production systems and is a critical component for national food security. Soil quality has become a key concern internationally due to competing demands between people, resources, and the environment [1][2][3][4][5][6][7]. Specifically, the potentially toxic elements (PTEs) such as arsenic (As), cadmium (Cd), chromium (Cr), mercury (Hg), lead (Pb), copper (Cu), nickel (Ni), and zinc (Zn) in soil have received significant attention in the past few decades due to their widespread distribution and severe toxicity to plants, animals, and humans. Oral ingestion (such as unintentional/intentional hand-oral ingestion), respiratory inhalation (as dust via respiratory inhalation), and dermal contact are the main routes through which PTEs from soil can be transferred into the human body [8][9][10][11][12][13]. Excessive intake or long-term exposure to PTEs may cause adverse effects to

Study Area
The Baoshan area is located in Yunnan province, Southwest China (Figure 1a), has a low latitude and moderate elevation, and is tempered by a mild subtropical climate with short, mild, and dry winters and warm, wet summers. Variations in the average temperature are closely connected to altitude, while the western and southern parts of the study area experience the highest amount of rainfall [54,55,64]. The main land-use types in the study area are forested land, dry cultivated land, paddy land, grassland, and garden plots, with few water areas or other land types (residential quarters, industrial and mining land, and unused land) (Figure 1c) [54,55]. Geologically, Baoshan is located in the Sanjiang area of Yunnan province, which consists of the Tengchong, Baoshan, and Lanping-Simao blocks (Figure 1b) [65][66][67]. The rock types of the Baoshan Block are mainly Paleozoic to Mesozoic sedimentary rocks and Early Paleozoic and Late Mesozoic-Early Cenozoic magmatic rocks [65,68]. Precambrian to Quaternary strata and sporadic granite and diabase intrusions are found in the study area ( Figure 1d). Basalt occurs in the Carboniferous and Jurassic rock sequences that are mainly distributed in the western and northeastern parts of the study area (Figure 1d), whereas mineral deposits, which often lead to PTE pollution in soil, are widely distributed across the study area, which include large Cu-Pb-Zn and Hg workings (Figure 1d) [54,55].

Samples Collection, Preparation, Chemical Analysis and Quality Control
Sample collection, preparation, chemical analysis, and quality control in the study area have been described in detail in a prior study [54].
In this study, a total of 1574 surface soil (0-20 cm) and 425 deep soil (150-180 cm or soil from the C horizon if the soil depth is less than 150 cm) composed datasets were used [54], which were collected in 2016 as part of the NGSLQ project in the Baoshan area (6220 km 2 ) using a systematic regular grid-sampling method [69]. In the research area, surface soils were collected at a density of 1 sample/km 2 and deep soils were collected at a density of 1 sample/4 km 2 . To assess sampling variation, about 1% of the sample sites were re-sampled. The samples were disaggregated with a wooden hammer during air drying and were sieved using a 20 mesh screen (0.84 mm) after drying. They were further processed with a −200 mesh screen (0.074 mm) for analysis. Surface soil samples from 4 km 2 and deep soil samples from 16 km 2 were composited for analysis to lower the analytical cost [4].
The total content of 54 elements or indictors, including PETs (such as As, Cd, Cr, Hg, Pb, Cu, Ni, and Zn), and the pH of the composited soil samples were analyzed at the Central Analytical Laboratory of the Institute of Geophysical and Geochemical Exploration (IGGE). Samples were analyzed as dry weight. As was analyzed by means of Hydride Generation-Atomic Fluorescence Spectrometry (HG-AFS) after digestion by aqua regia; Cd, Cu, Ni, Pb, and Zn were analyzed by Inductively Coupled Plasma-Mass Spectrometry (ICP-MS) after digestion by hydrofluoric acid (HF), nitric acid (HNO 3 ), perchloric acid (HClO 4 ), and aqua regia; Cr was analyzed by X-Ray Fluorescence Spectrometry (XRF) after powder pellet formation; Hg was analyzed by Cold Vapor-Atomic Fluorescence Spectrometry (CV-AFS) after digestion by aqua regia; pH was analyzed by Ion-Selective Electrodes (ISE) dissolved in CO 2 -distilled water. During the routine analysis, internal and external controls were introduced to ensure accuracy and precision. Briefly, Certified Reference Materials (CRMs) and Blind Reference Materials (BRMs) were analyzed with the samples simultaneously. Approximately 5% of the samples were chosen at random and were re-analyzed. Within each batch of 50 samples, 4 BRMs were placed in a cipher and examined at the same time. The analytical methods, detection limits and detailed analytical procedures met the analytical requirements set in the NGSLQ project for all elements in all samples [69].

Pollution Estimation and Risk Assessment
Certain indicators have been used in environmental studies on PTEs in soil to assess the severity of pollution, potential ecological risk, and potential health risks to human beings that they cause [70][71][72].

Pollution Evaluation Indicators
To assess soil contamination levels, important indices such as the enrichment factor (EF), geo-accumulation index (I geo ), and contamination factor (C f ) are widely used [5,52,[73][74][75]. The indicators are calculated using the Equations below for completeness: The EF is a valuable geochemical instrument for assessing the degree of PTE contamination in soil, where "soil" denotes surface soil, "background" denotes background, "C i " is the concentration of the selected PTE, and "C r " denotes the concentration of the reference element.
where "C i " denotes PTE concentration in surface soil, "C b " denotes the background PTE concentration. This indictor represents the impact or contribution of an individual PTE on soil contamination. The C f can be divided into four categories: C f < 1: low contamination; 1 < C f < 3: reasonable contamination; 3 < C f > 6: considerable contamination; C f > 6: high or very high contamination [72]. The background level was measured using the deep soil from the 1:250,000 scale land quality geochemistry data from the Baoshan area, since it is considerably less likely to be impacted by human activity [54,79]. The median values of the PTEs in the deep soil of the Baoshan area were used as the background values to calculate Equations (1)-(3).

Ecological Risk Assessment
The potential ecological risk (E i r ) of a given PTE and the potential ecological index (PERI), which indicates the sum of potential ecological risks, are commonly employed at the regional and/or global levels to assess the potential impact of the pollutants in soil on ecosystems [72].

Human Health Risk Assessment
A variety of preventative methods have been developed to evaluate the health risks caused by contaminants in order to maintain public health. The potential human health risk models raised by the US EPA, which include carcinogenic and non-carcinogenic risks [80], have proven successful and have been used globally [5,13,52,74,81].
The direct exposure of children and adults to PETs through soils was taken into account when calculating the human health risk.
The Equations. below were used to calculate chronic daily intake (CDI) of PTEs exposure ingestion (CDI ingestion ), inhalation (CDI inhalation ), and dermal contact (CDI dermal ) [74]. The exposure factors for health risk assessment mentioned in the following Equations. are listed in Table 1. Note: In the case of Cr, 1/6 of the total concentration was used in the risk calculations for both non-carcinogenic and carcinogenic risk [83,84].
Non-carcinogenic hazards are typically characterized by the hazard quotient (HQ). The HQ is defined as the quotient of the chronic daily intake divided by the toxicity threshold value, which is referred to as the reference dose (R f D) of a specific PTE. The HQ of a single PTE is determined by Equation (9). Equation (10) was used to analyze the health risk related to PTE exposure for carcinogenic risk [80]. The R f D and slope factor (SF) through each PTE exposure route are listed in Table 2.
where the subscript "i" in Equations (9) and (10) indicate the three main PTE intake routes: ingestion, inhalation, and dermal contact. Only As, Cd, Cr, Ni, and Pb were calculated for the carcinogenic risk of the SF is given. Note: "-" not available.
For an individual PTE in soil, the possible non-carcinogenic and carcinogenic risks are expressed as the hazard index (HI) and total cancer risk (TCR). HI is the sum of the HQ of an individual PTE through the three intake routes. TCR is the sum of the CR through the three intake routes. The Equations. for the calculations are listed below: HI = HQ ingestion + HQ inhalation + HQ dermal (11) TCR = CR ingestion + CR inhalation + CR dermal (12) However, because HI and TCR refer to the risks of a single PTE, their cumulative values were calculated to provide a more complete picture; this shows the entire impact of all PTEs across all exposure pathways. The cumulative hazard index (CHI) and cumulative total cancer risk are calculated using the following Equations [86]: If HI ≤ 1, there is no significant risk of non-carcinogenic effects, and the magnitude of risk increases as HI increases; HI > 1 is thought to represent a significant risk of noncarcinogenic effects [82]. CR < 10 −6 is considered to be negligible; 10 −6 < CR < 10 −4 represents an acceptable risk that is under control; CR >10 −4 is considered unacceptable [80].

Data Statistical Analysis and Mapping
The descriptive statistics of PTEs such as minimum value, maximum value, mean value, the median, and standard deviation (SD) were conducted by SPSS 21.0 (IBM, Armonk, NY, USA). ArcMap TM 10.2 (Esri Inc., Redlands, CA, USA) was used to examine the spatial distribution pattern and variability of the PTE contamination statues, ecological risk, and health risk. Table 3 shows the descriptive statistics of the analyzed PTE concentrations and pH values of the collected surface soils from the Baoshan area, as well as local, regional, national, and relevant reference guideline values. The pH values of the surface soils from the research area ranged from 4.13 to 8.33 (median: 6.46; mean: 6.45). The registered median values for the concentrations of As, Cd, Cr, Cu, Hg, Ni, Pb, and Zn were 16, 0.19, 113, 41, 0.10, 51, 34, and 100 mg/kg, respectively. All PTEs were above their national background values. The registered median values for the concentrations of As, Cd, Cr, Cu, Hg, Ni, Pb, and Zn were 1.78, 1.73, 2.26, 2.16, 5.56, 2.32, 1.62, and 1.67 times higher than their corresponding national background values, respectively [87]. With the exception of Pb, all of the studied PTEs were above their regional background values (background values of Yunnan province) [88]. In particular, the median values of the Cd, Hg, and Cr concentrations were 2.29, 2.00, and 1.66 times higher than their corresponding regional background values, respectively. However, no noticeable enrichments (>1.5 times) or deficits (<0.8 times) were identified when compared to the local background values of the studied PTEs [54]. When the PTE contents in the surface soil from the Baoshan area were compared to those found in the Australian surface soil [89], the American soil A layer [90], and the surface soil of the European continent [91] (Table 3), the contents of the PTEs in the surface soil of the study area were about 1.5 times higher or more than those in the surface soil of the American, Australian, and European continents. This fact indicates that the study area is a typical geological background area with high PTE levels, and the elevated PTE concentrations in the surface soil are mainly due to geological genesis [54]. The PTE contents in the surface soil of the study area are similar to those in the Sacramento Valley of Northern California, USA [60]. This area is similar to the Baoshan area, which is widely characterized with ultra-basic rocks and mineral deposits (such as As, Hg, Zn, Cu, Au, and Ag) [60]. The concentrations of As, Cd, Cu, Hg, Pb, and Zn in the surface soil collected from the study area were significantly lower than those found in the soil around the mining areas in other research studies [5,92] and in the surface soils deposited in areas with carbonate rocks as the main parent material [93]. Meanwhile, the Cr, Ni, and Cu concentrations in the present study area were determined to be significantly lower when compared to those in the surface soil in areas where ultra-basic rocks are the main parent materials [94]. Table 3. Descriptive statistics of PTEs in the surface soil of the study area and the local, regional, and national background values and related guidelines. The number of sites and the percentage from which the collected surface soils registered values of PTEs (As, Cd, Cr, Cu, Hg, Ni, Pb, and Zn) that exceeded the risk screening values (RSV) and risk intervention values (RIV) for soil environmental quality-risk control standards for soil contamination of agricultural land in China-are listed in Table 4. If the concentration of a PTE in agriculture land soil is lower than or equal to the concentration of RSV, then the ecological risk of this PTE in the soil is low, whereas if the concentration of a PTE in agriculture land soil is higher than the RSV, then the ecological risk for this PTE in agricultural soil is high. In the latter case, the soil should be environmentally monitored. If the concentration of PTE is higher than the RIV, strict intervention and controls should be taken [95]. The results showed that the PTEs in certain surface soil sites in the study area exceed the RSV. The site percentage of surface soils with values that exceeded the RSV were as follows: As (26.37%), Cd (17.41%), Cu (16.71%), Cr (14.04%), Ni (11.8%), Pb (3.30%), Zn (3.3%), and Hg (0.25%) ( Table 4). In case of the RIV values, only a few sites with As, Cd, and Pb were above the limit. The site percentage of surface soils with values exceeding the RSV were as follows: As (0.76%), Cd (0.57%), and Pb (0.38%). To demonstrate the distribution of PTEs at an environmental quality level, figures are plotted by assigning the sample site values to 2 km × 2 km grids (this is the sample analysis density). The spatial distribution of the sites with Cr, Ni, and Cu values that were higher than the RSV and RIV corresponds to the basalt and carbonatite distribution in the studied area. The spatial distribution of the samples with Cu, Pb, and Cd concentrations exceeding the RSV and RIV is compatible with the distribution of Cu-Pb-Zn deposits and a portion of the carbonatite rock distribution. The distribution of the sites with As and Hg concentrations exceeding the RSV and RIV is compatible with the distribution of Hg and Cu-Pb-Zn deposits, as well as a portion of the carbonatite rock distribution (Figure 2). This finding shows that soils formed from basalt and carbonatite rocks, as well as ore parent materials or those that are formed near mining areas, have higher PTE concentrations [2,[33][34][35][36][37]54,55,57,60,62].

Evaluation by Pollution Indices
The values and descriptive statistics of the selected pollution indices are shown in Table 5. For this study, we used Sc as the reference element because of its abundance and stability in the upper continental crust, its uniformity across the entirety of the study area, and its effective use in prior studies [94,96,97]. Appl. Sci. 2022, 12, x FOR PEER REVIEW 11 of 26

Evaluation by Pollution Indices
The values and descriptive statistics of the selected pollution indices are shown in Table 5. For this study, we used Sc as the reference element because of its abundance and  There was no significant enrichment or contamination with the studied PTEs in the surface soils in the study area, based on the average (mean) and median values of EF, I geo , and C f that were obtained (Table 5). With the exception of Cr, the rest of the PTEs studied in this area showed significant to extremely high enrichment by the indicated EF in some locations. Of the elements with EF values that were >5, 2.16% (34 sites) of As, 4.51% (71 sites) of Cd, 0.44% (7 sites) of Cu, 4.00% (63 sites) of Hg, 0.06% (only 1 site) of Ni, 1.97% (31 sites) of Pb, and 0.64% (10 sites) of Zn, these were classified as having a significant contamination or higher level. Specifically, the maximum EF values of As (49.51), Cd (32.14), Cu (16. (19.38) were also determined. Contour maps using inverse distance-weighted (IDW) interpolation were created to analyze the locations of the contaminated PTEs and their affected region, since IDW is a simple, intuitive, and efficient interpolation method and demonstrates good results when the distribution of the sample sites is uniform (Figures 3-5). The spatial distributions of Cu, Pb, and Zn with EF values > 5 (significant contamination or higher levels) coincided with the Cu-Pb-Zn deposits. Cd showed a similar distribution pattern, but it also coincided with the carbonatite rock distribution. The distributions of As and Hg with EF values > 5 are consistent with the distribution of the Hg and Cu-Pb-Zn deposits and a portion of carbonatite rocks (Figure 3). Similar distributions were found using I geo and C f values (Figures 4 and 5).

Assessing the Environmental Ecological Risks of the Surface Soil
To assess the potential impact of a given PTE in the soil on ecosystems, (E i r ) was used in this study. The PTEs in the surface soil of the research area did not demonstrate a risk degree higher than "moderate risk" according to the average (mean) and median values E i r obtained at all sampling sites. While Cd and Hg posed higher ecological risks when compared to the other PTEs by the upper quarter values, Cd and Hg showed examined values of 64.36 and 61.17, respectively, due to the higher toxic response factor. Furthermore, the highest E i r values of Hg (2426.50), Cd (1229.28), As (485. 25), and Pb (343.60) were examined, due to their greater toxic response factors of 40, 30, 10, and 5 when compared to Cu (five), Ni (five), Cr (two), and Zn (one), respectively. The environmental ecological risks of all of the studied PTEs can be ranked as: Cd > Hg > As ≈ Pb > Cr = Ni = Zn. Overall, 0.64% (10 sites) of As, 16.07% (253 sites) of Cd, 16.39% (258 sites) of Hg, and 0.51% (eight sites) of Pb were shown to have E i r values> 80, and 7.05% (11 sites) of all of the studied PTEs with PERI values > 300 were at or above the level of "considerable protentional ecological risk".
Contour maps determined by IDW interpolation were produced to discover the spatial distributions of the possible ecological risk of these studied PTEs ( Figure 6). The distributions of As, Cd, Hg, and Pb, which had E i r values > 80 and PERI values > 300, were compatible with the distribution of Hg and Cu-Pb-Zn deposits in the study area. The high PTE concentrations and their associated ecological risks in the soils surrounding Hg and Cu-Pb-Zn deposits have also been found by other researchers [5,52].

Non-Carcinogenic Risks to Humans
The median values of non-carcinogenic risk assessment indictors, such as CDI, HQ, HI, and CHI are listed in Table 6. The median values of the non-carcinogenic CDI and HQ of all of the studied PTEs by direct oral ingestion were substantially higher when compared to those by dermal contact and respiratory inhalation. The median values of daily non-carcinogenic exposure and the health risks associated with the three different pathways were ranked as follows: CDI ingestion > CDI dermal > CDI inhalation and HQ ingestion > HQ dermal > HQ inhalation . For children, the order of the median values of the non-carcinogenic risks of PTEs from the three pathways was as follows: As > Pb > Cr >Ni > Cu > Hg > Zn > Cr, while for adults it was As > Cr > Pb > Ni > Cu > Cd > Zn > Hg. In the case of children, the non-carcinogenic health risks were higher than those in adults for all of the PTEs, with the exception of dermal contact, which has a lower non-carcinogenic daily exposure.
There were no identified non-carcinogenic risks for humans in case of the individual PTEs by the three pathways or by the sum of them, based on the median values, in either children or adults. However, the conclusions are not the same when we consider all sampling points in the research area. Through direct oral ingestion, the majority of the sites represent no significant risk of non-carcinogenic consequences for children and adults. While 19.38% (305 sites) of As and 0.51% (eight sites) of Pb pose a non-carcinogenic risk to children and 0.13% (two sites) of As pose a non-carcinogenic risk to adults through the direct oral ingestion pathway. For all of the pathways, 19.57% (308 sites) of As and 0.51% (eight sites) of Pb pose non-carcinogenic risks for children and 0.25% (four sites) of As pose a non-carcinogenic risk for adults, respectively. All of the total studied PTEs through the three pathways pose non-carcinogenic risks for children and adults were 32.08% (505 sites) and 0.57% (nine sites), respectively. No sites, on the other hand, represented a non-carcinogenic risk to children and adults through inhalation or dermal contact (Table 7). Figure 7 shows the distribution of the non-carcinogenic risk of PTEs to children and adults determined by IDW interpolation. The distribution of sites with significant noncarcinogenic risk to humans is in accordance with the distribution of the Hg and Cu-Pb-Zn deposits. Children are more vulnerable to environmental contaminants when compared to adults due to behavioral and physiological differences [98]. As poses the highest noncarcinogenic risk when compared to the other studied PTEs because it has the lowest non-carcinogenic daily intake reference dose permitted by oral ingestion.  Figure 7. The distribution of counter maps of the non-carcinogenic risks (HI) and cumulative noncarcinogenic risk (CHI) determined by inverse distance-weighted (IDW) interpolation.

Carcinogenic Risks to Humans
Only the carcinogenic risks of As, Cd, Cr, Ni, and Pb to humans is examined in this chapter because only the carcinogenic slope factors of these researched PTEs have been previously determined in the literature. Table 8 shows the median values of carcinogenic risk assessment indictors, such as CDI and CR, for individually documented PTEs by individual routes, TCR for individual PTEs by total pathways, and CTCR for the cumulative values of these recorded PTEs by the total three pathways. Similar to non-carcinogenic risk, the median values of the carcinogenic CDI for all of the studied PTEs by direct oral ingestion was much higher than that by dermal contact and respiratory inhalation. The median values of the daily carcinogenic exposure of the three different pathways were ranked as follows: CDI ingestion > CDI dermal > CDI inhalation . Only the carcinogenic slope factors of As and Cr were provided by direct oral ingestion. Children have a substantially higher carcinogenic risk from As and Cr exposure than adults. Adults have higher carcinogenic hazards than children when it comes to respiratory inhalation and dermal contact. For children, the order of the median values of carcinogenic risk for the sum of the documented PTEs from the three overall pathways was as follows: As > Cr > Pb >Ni > Cd, and for adults, it was As > Cr > Ni >Pb > Cd. The carcinogenic risk of children to all of the documented PTEs by the total pathways were higher than that of adults.
When considering all of the sample points of the study area, the carcinogenic risk by oral ingestion were considered to be acceptable for As (1532 sites, counting for 97.33%) and Ni (1571 sites, counting for 99.81%) for children and As (1561 sites, counting for 99.17%) and Ni (1569 sites, counting for 99.68%) for adults, respectively. Additionally, the carcinogenic risks of one site for As and three sites for Cr to children, and two sites for As and five sites for Cr to adults were considered to be negligible. For As, only 41 sites (2.60%) and 11 sites (0.70%) represented an unacceptable carcinogenic risk to children and adults via the oral ingestion pathway, respectively. For inhalation, all of the sites represented a negligible carcinogenic risk to both children and adults. With the exception of 2.67% (42 sites) and 0.76% (12 sites) of As, which pose a carcinogenic risk to children and adults, respectively, all of the other known PTEs pose a negligible or acceptable carcinogenic risk through all three pathways. In terms of the cumulative carcinogenic risk of all of the recorded PTEs via the three total pathways, 3.24% (51 sites) and 0.83% (13 sites) were considered to be unacceptable for children and adults, respectively. Arsenic contributed the most to the carcinogenic risk of all of the PTEs (Table 9).
Contour maps determined by IDW interpolation are presented to find the carcinogenic risk distributions of As and all of the PTEs via the three pathways ( Figure 8). The distribution of places with unacceptable carcinogenic risk coincides with the distribution of Hg and Cu-Pb-Zn deposits, just as it does with non-carcinogenic risk. In addition, As had the highest slope factor of all of the examined PTEs, meaning that it poses the highest carcinogenic risk to humans. Although the overall carcinogenic risk in the study area is acceptable, residents living near the mine should be aware of the potential carcinogenic risk. The PTEs in the soil around mining areas demonstrate considerable non-carcinogenic and carcinogenic risks for local inhabitants, and this has been verified by a previous study by analysis of a large number of soil data around the mining area [5].
acceptable, residents living near the mine should be aware of the potential carcinogenic risk. The PTEs in the soil around mining areas demonstrate considerable non-carcinogenic and carcinogenic risks for local inhabitants, and this has been verified by a previous study by analysis of a large number of soil data around the mining area [5].  Figure 8. The distribution of counter maps of the carcinogenic risks (TCR) and cumulative carcinogenic risk (CTCR) determined by inverse distance-weighted (IDW) interpolation. Table 9. The sites number of sites (percentage) of carcinogenic risk levels to humans from soil PTEs.

Conclusions
The following conclusions can be drawn from the above research: (1) The PTE values determined from partial surface soil samples from sites in the study area exceeded the RSV and RIV of soil environmental quality-the risk control standard for soil contamination of agricultural land of China, showing ecological risks; (2) Cd and Hg posed higher ecological risks than the other PTEs (As, Cr, Pb, Cu, Ni, and Zn) highlighted by their high toxic response factor; (3) Arsenic is the main PTE that posed a non-carcinogenic risk, and the only PTE that poses a carcinogenic risk to humans. Children are more vulnerable to health risks than adults because of their behavioral and physiological traits; (4) Geological genesis was responsible for the high concentrations, ecological risk, and health risk distribution patterns of the examined PTEs. The ecological and health risks of the PTEs in the surface soil that had developed from basalt and carbonatite rocks, particularly near mineral resources, were discovered in the study area. Future administration and research are needed, especially in mining areas.  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. The data are not publicly available due to the data are unified management by the project administrator.