Concentrations, Distribution, Sources and Ecological Risk Assessment of Trace Elements in Soils from Wuhan, Central China

This study aimed to determine the concentration levels, potential sources and ecological risks of eleven trace elements, namely Cr, Fe, Co, Ni, Cu, As, Sb, Cd, Zn, Hg and Pb, in the soil from Huangpi district, Wuhan, Central China. Soil samples were collected from eighteen sites at soil depths of 1–10 and 10–20 cm and analyzed using Inductively Coupled Plasma-Mass Spectrometer ICP-MS (Thermo X SERIES 2, Scientific and Innovative Technology Co. Ltd., Beijing, China). The recorded mean concentration of the elements were in a decreasing order of Fe > Co > Cr > Ni > Pb > Cu > As > Cd > Sb > Zn > Hg. The mean concentration of trace elements, soil pH and total organic carbon (TOC) were higher at a soil depth of 1–10 cm. The obtained mean concentration of Cr, Co, As, Cd, Ni, Cu, Hg and Pb were above the soil background values of Wuhan and Hubei Province. The mean concentration values of Co, Ni and Cd, exceeded the recommended FAO (Food and Agriculture Organization)/ISRIC (International Soil Reference and Information Centre) (2004) and WHO/FAO (2001) values. Pearson’s correlation analysis illustrated that there was a strong and significant correlation between trace elements, whereas, a weak positive and negative correlation between elements and soil properties (pH and TOC). The principal component analysis (PCA) and cluster analysis (CA) result indicated that the concentration of trace elements in Huangpi soil were originated from anthropogenic sources. Potential ecological risk index (RI) of this study revealed that there is a high ecological risk of trace elements in the soil. Enrichment factor (EF) and geo-accumulation index (Igeo) of trace elements for this study indicated that the study area is strongly contaminated with Cd and Co. Generally, the finding of this research showed that Huangpi soil is contaminated.


Sample Collection and Pretreatment
Soil samples were collected from eighteen sampling sites of different land use types (barren land, farmland, paddy field, plastic greenhouse) at soil depths of 1-10 and 10-20 cm. The sample sites were located using a global positioning system (GPS). For each site, three replicates were taken to make up a composite soil. The collected samples were packed in polyethylene containers, labeled and transported to the laboratory for analysis. All samples were dried using a benchtop lab vacuum freeze dryer (Xinzhi Biotechnology Co., Ltd., Ningbo, China) at −40 • C for 24 h, the samples were ground, sieved through 0.6 mm mesh nylon sieve and stored in a refrigerator at −20 • C for further analysis [12,36].

Sample Extraction and Analysis
About 0.1-0.15 g of soil samples were measured using a digital analytical balance (METTLER TOLEDO, Columbus, OH, USA), placed into Teflon vessels, digested with 4 mL of nitric acid (HNO 3 , 63%), 2 mL of hydrogen fluoride (HF ≥ 40%) and perchloric acid (HClO 4 70-72%) using a microwave digester (ETHOS ONE, Milestone, Leutkirch im Allgau, Germany). The digester was operated for 2 h at a controlled pressure, temperature and output power according to [37]. The digested samples were heated on the heating plate at a temperature of 135 • C for 2 h and cooled to room temperature. Then the final volume was topped up to 50 mL using double distilled water and filtered using 0.22 um membrane filter paper for analysis [22,36] and the filtered samples were stored in plastic bottles at −4 • C to minimize volatilization and biodegradation [38]. Finally, the samples were analyzed using ICP-MS (Thermo XSERIES2, Beijing, Scientific and Innovative Technology Co. Ltd., Beijing, China). In addition, selected soil properties (soil pH and TOC) were measured. Soil pH was measured according to [22], using a pH meter (METTLER TOLEDO) and the total organic carbon (TOC) of the soil was measured by [39] method using a TOC analyzer (Elementar GmbH, Langenselbold, Germany).

Standard Preparations and Calibration Curves
All chemicals and reagents used were analytical grades with purity of 99%. Five working solutions with concentrations of 0.1, 0.5, 1, 5, 10 and 20 µg/L were prepared from the stock solutions of 1000 µg/mL using double distilled water with 10% nitric acid (HNO 3 , 63%). The working solutions were analyzed before running samples to check the accuracy and reliability of the instrument. Calibration curves produced for all elements indicated that the obtained curves had r 2 values greater than 0.994.

Statistical Analysis
Statistical Package for Social Sciences (SPSS) version 20 (IBM, Armonk, NY, USA) and Microsoft Offices Excel 2013 (Microsoft Corporation, Albuquerque, NM, United States) were used to analyze the data. Pearson's correlation coefficient was used to analyze the relationships between elements, soil pH and TOC. The potential sources of trace elements were identified using factor analysis (PCA and cluster analyses).

Quality Analysis and Quality Control
To ensure the quality of the experiment all reagents and chemicals used were analytical-reagent grades. All plastic containers were soaked with 10% HNO 3 (63%) for 12 h and washed three times with double distilled water to remove other contaminants, dried in oven at 60-65 • C for 24 h. The relative standard deviations (RSD) of three consecutive measurements of the standard solutions were used to determine the precision of ICP-MS. The obtained percentage of relative standard deviation (% RSD) was less than 10% indicating a good precision of the instrument. The recovery value of elements ranged between (90-104%).

Potential Ecological Risk Index (RI)
Potential ecological risk index helps to evaluate the pollution of trace elements in the soil [40]. It is the sum of the ecological risk factor of a single element in the sample [41]. The following formulas given by [41] were used to calculate potential ecological risk indices of the elements.
where C i f is the contamination factor for trace elements; C i S is the measured concentrations value of elements in the soil; C i R is the background reference values of trace elements in the soil; C d is a degree of contamination; E i f is the ecological risk index of a single element; T i f is the toxicity coefficient of an element and RI represents the total potential risk index of elements [40] (Table 1).

Enrichment Factor (EF)
Enrichment factor (EF) was evaluated to determine the degree of anthropogenic factors on trace elements accumulation in the soil. It reflects the disturbance degree of human activities on the natural environment [44]. EF is calculated as follows: where C i is the concentration of trace elements in the sample; C ref is the concentration of reference element in the sample; B i is the background value of interest element and B ref is the background value of reference elements in the study area. Fe was selected as the reference element for this study because Fe is a major sorbent and it is a quasi-conservative tracer of natural elements in fluvial and coastal sediments [12]. Wuhan soil background was used as the reference. The following EF class were given for elements; EF < 2, 2 ≤ EF < 5, 5 ≤ EF < 20, 20 ≤ EF < 40 and EF ≥ 40 indicated low, moderate, significant, very high and extremely high enrichment factor, respectively.

Geo-accumulation Index (I geo )
The geo-accumulation indexes (I geo ) of elements were calculated to identify the degree of contamination and to compare it with the pre-industrial level [45]. This method classifies the pollution level of elements in terms of seven enrichment classes. It is calculated using the following formula: where I geo is the geo-accumulation index, C n is the measured concentration of an element (mg/kg) in the sample, B n is the geochemical background value of an element (mg/kg) and 1.5 is the factor used for lithological variations of elements. Seven classes and contamination intensity are given for I geo : I geo ≤ 0 (uncontaminated), 0 < I geo ≤ 1 (uncontaminated to moderately contaminated), 1 < I geo ≤ 2 (moderately contaminated), 2 < I geo ≤ 3 (moderately to strongly contaminated), 3 < I geo ≤ 4 (strongly contaminated), 4 < I geo ≤ 5 (strongly to extremely contaminated) and I geo > 5 (extremely contaminated) [45].

Spatial Distribution of Selected Soil Properties (pH and TOC)
It has been reported that soil properties such as soil pH and total organic carbon (TOC) are the most important factors that influence cation mobility and regulate the solubility of trace elements in the soil. The obtained mean result of soil pH and TOC for this study were presented in Table 2. As illustrated, the mean soil pH results ranged from 4.20-6.87 with a mean value of 5.71. The comparison between individual samples (Table 3) indicated the soil pH ranged from 4.12-8.07 with an average of 5.87 at the soil depth of 1-10 cm and 4.28-6.52 with an average of 5.55 at the soil depth of 10-20 cm. The highest soil pH value (8.07) was recorded in the greenhouse soil from Changdi site at a soil depth of 1-10 cm. The main reason for this might be connected with relatively low precipitation amount and less leaching of base-forming cation in the green house. The mean comparison between the two depths indicated that the soil pH at 1-10 cm (5.87) is higher than the soil depth of 10-20 cm (5.55) (Figure 2). The result is in line with the findings of [46,47], whereas the result is different from the findings of [48]. According to the classification of soil pH grade, the pH value of (<5) indicated slightly acidic, (5-6.5) mildly acidic, (6.5-7.5) neutral, (7.5-8.5) mildly alkaline and (>8.5) indicated strongly alkaline [49]. According to these classifications, Huangpi soil is classified as slightly acidic, mildly acidic, neutral and mildly alkaline.
The average total organic carbon of soil ranged from 0.65-2.41 with the mean value of 1.71. The highest mean value of TOC (2.41) was recorded in Dujiatian site (farmland), whereas the lowest TOC (0.65) was recorded from barren land at the Zhoujiawan site. TOC ranged from 0.63-2.60 with a mean value of 1.74 at soil depth of 1-10 cm and 0.66-2.71 with a mean value of 1.71 at the soil depth of 10-20 cm. The highest mean was recorded at the soil depth of 1-10 cm ( Figure 2). A similar result was reported by [47,50]. Numerous studies indicated that TOC of soil decreases with soil depth [50,51]. The main reason for this might be biological activity in the top layer of soil [51].

Concentrations and Distribution of Trace Elements
Eleven trace elements in Huangpi district soil were studied. The results of descriptive statistics; mean, maximum, minimum, standard deviation and Skewness were presented in Table 2 and individual results of all samples and mean for both depths were presented in Table 3 and Figure 2, respectively. The results indicated that the concentrations of the elements showed variation between samples and within samples at different soil depths. All elements except As, Sb, Hg and Cd were detected in all samples at soil depths of 1-10 and 10-20 cm with detection frequency of 100% for Cr, Fe, Co, Ni, Cu, Zn and Pb; and 72%, 89%, 78% and 42% for As, Sb, Cd, Hg respectively. The mean concentration of trace elements at soil depth of 1-10 cm were in decreasing order of Fe > Co > Cr > Ni > Pb > Cu > Cd > As > Zn > Sb > Hg, whereas, Fe > Co > Cr > Ni> Pb > Cu > As > Cd > Zn > Sb > Hg at soil depth 10-20 cm ( Figure 2).The highest mean concentration values for all elements were recorded at a soil depth of 1-10 cm. The results of this study was in line with those of [52][53][54]. According to Camobreco [55] the highest accumulation of trace elements in the surface layer of the soil might be due to a high sorption capacity of trace elements which results from soil chemical reactions between soil solid phases, including silicate clays, hydroxides and oxides of elements. Another study by Rahaman [53] also indicated that trace elements were found abundantly in the surface layer and the value decreased with an increase in soil depth with few exceptions. Converse to the results obtained in this study, the report by [56] from Kenya, indicated that subsurface soil accumulates high concentration of trace elements than surface soil due to soil leaching.    The highest mean of Cr was recorded in the sample from Dujiatian (farmland), whereas the lowest was recorded at Tangjiawan (barren land). The obtained mean result of Cr (140.10 mg/kg) was higher than the soil background value of Wuhan (90 mg/kg) and Hubei Province (86 mg/kg), but less than soil background of China (200 mg/kg). In comparison to the permissible limit of FAO (Food and Agriculture Organization)/ISRIC (International Soil Reference and Information Centre) (2004) (250 mg/kg) and USEPA (1983) (1000 mg/kg) the mean concentration value of Cr (140.10 mg/kg) was lower. However, the individual mean result of single site indicated that the sample from Dujiatian tea farmland (321.73 mg/kg) was above the permissible limit of FAO/ISRIC (2004). As compared to the other findings in Wuhan, the obtained mean value of Cr for this study was higher than the finding of [32] (85 mg/kg), but less than the finding by [33]  The highest mean of Cr was recorded in the sample from Dujiatian (farmland), whereas the lowest was recorded at Tangjiawan (barren land). The obtained mean result of Cr (140.10 mg/kg) was higher than the soil background value of Wuhan (90 mg/kg) and Hubei Province (86 mg/kg), but less than soil background of China (200 mg/kg). In comparison to the permissible limit of FAO (Food and Agriculture Organization)/ISRIC (International Soil Reference and Information Centre) (2004) (250 mg/kg) and USEPA (1983) (1000 mg/kg) the mean concentration value of Cr (140.10 mg/kg) was lower. However, the individual mean result of single site indicated that the sample from Dujiatian tea farmland (321.73 mg/kg) was above the permissible limit of FAO/ISRIC (2004). As compared to the other findings in Wuhan, the obtained mean value of Cr for this study was higher than the finding of [32] (85 mg/kg), but less than the finding by [33] (152.78 mg/kg). In comparison to the other studies from other places, the mean concentration of Cr for this study was higher than the findings from Cuba (85.9 mg/kg) [59], Bangladesh (53.7, 34.2 mg/kg) [60], India (8.01 mg/kg) [61], Tanzania (7.68 mg/kg) [10], Brazil (20.61 mg/kg) [62], Pakistan (5.86 mg/kg) [63], Iran (48.08, 53.21 mg/kg) [64], along Chao River in China (118 mg/kg) [65] and Northern Pakistan (29.94 mg/kg) [66].
The maximum mean value of Fe (55,398.01 mg/kg) was recorded at Dujiatian farmland site, whereas the lowest value was recorded at Tangjiawan site in the soil from a barren land. The recorded mean value of Fe was below the soil background value of Hubei Province and China (29,400 mg/kg). As compared to the other studies, the mean concentration of Fe (27,304.9 mg/kg) from this study was higher than a study from Bangladesh [11] (1800 mg/kg).The highest concentration of Co (55,580.36 mg/kg was obtained at the Dujiatian site (farmland), whereas, the lowest value (4355 mg/kg) was recorded at Tangjiawan site (barren land). The obtained mean Co was higher than the soil background value of Wuhan, Hubei Province and China ( Table 2). The mean concentration value of Co (22,656.94 mg/kg) was above the permissible limit of FAO/ISRIC (2004) and WHO/FAO (2001) ( Table 2). Compared to other findings the mean concentration of Co was higher than the study in Cuba (9.16 mg/kg) [59], Brazil (7.44 mg/kg) [62], Pakistan (7.56 mg/kg) [63], Wuhan (China) (16 mg/kg) [32] and (16.37 mg/kg) [33], Iran (38.5, 16.51 mg/kg) [64], Chao River China (17.5 mg/kg) [65] and Northern Pakistan (36.76 mg/kg) [66].
The highest value of As was recorded at Zhulinyuan site from a barren land, whereas the lowest was recorded at Hanjiafan from a paddy field. The mean concentration result of As was higher than soil background values of Wuhan and Hubei province ( Table 2). The average concentration of As recorded was below the permissible limit of USEPA (1983). However, the individual result of each site indicated that the concentrations of As in samples from Fengdouhu, Changdi, Zhulinyuan, Zhoujiawan and Zhujiashan sites were above the soil background value of China (30 mg/kg). The mean comparison of As concentration of this study with the other studies are presented in Table 4. The obtained mean concentrations of Sb (0.58 mg/kg) was lower than the soil background value of Hubei (1.65 mg/kg) and China (10 mg/kg). The mean result of Sb in this study was lower than the finding from Brazil (13.81 mg/kg) [62], whereas higher than the finding from Iran (0.22, 0.21 mg/kg) [64]. The maximum value of Cd was recorded in the sample from Zhujiashan from a greenhouse, whereas the lowest was recorded at Tangjiawan from barren land. The mean concentration of Cd (15.44 mg/kg) was above the soil background values of Wuhan, Hubei Province and China. In addition the concentration of Cd in this study passes the permissible limit of FAO/ISRIC (2004), WHO/FAO (2001) USEPA (1983). In comparison to the other studies the concentration of Cd for this study was lower than the finding of [62] (13.81 mg/kg) and [67] (49 mg/kg) ( Table 4).
The mean concentration of Zn for this was lower than the soil background value of Wuhan, Hubei and China. As compared to the permissible limit of FAO/ISRIC (2004) (500 mg/kg), WHO/FAO (2001) (300 mg/kg), USEPA (1983) (500 mg/kg) and China 250 mg/kg the mean concentration of Zn in this study was lower.
The mean result of Hg was higher than the soil background value of Wuhan, Hubei Province and China ( Table 2), but below, however, the permissible limits of FAO/ISRIC (2004), WHO/FAO (2001) and USEPA (1983). However, the result of the individual site indicated that sample from Dujiatian farm (1.13 mg/kg) exceeded the permissible limit of WHO/FAO (2001) (1 mg/kg). In comparison to the other studies, the obtained mean concentration of Hg was lower than the finding of in soil from Cuba [59], Bangladesh [60], Brazil [62], Iran [64] and soil along the Chao River in China [65].

Relationships of Trace Elements and Selected Soil Properties (pH and TOC)
Pearson's correlation was applied to analyse the relationships between trace elements, soil pH and TOC. The obtained result indicated that there was a significant positive correlation between trace elements at a significant level of 0.01 and 0.05 (2-tailed) whereas, there was a negative and weak correlation between trace elements and selected soil propertis (pH and TOC). Strong positive correlations of Cr with Fe, Co, Ni, Cu, and Hg, Fe with Co, Ni, Cu and Hg, Co with Ni, Cu, and Hg, Cu with Hg, As with Cd, Sb with Cd, Zn and Pb and Zn with Pb were observed at a significance level of 0.01 and 0.05 (Table 5). Moderate and weak correlations between trace elements were also obtained ( Table 5). The relationships between trace elements and soil properties indicated that soil pH has a weak positive correlation with Fe, Ni, Cu, As, Sb, Cd and Pb, whereas a negative correlation with Cr, Co, Zn and Hg. The same finding was reported by [47], for Cr, As, Hg and Pb from Guangzhou (China). TOC had a negative correlation with Co, As, Sb, Cd and Hg, and a weak correlation with Cr, Fe, Ni, Cu, Zn and Hg. The same result was reported on Cd by [69]. Weak and negative correlations were observed between soil pH and TOC. A significant positive correlation among elements indicates a common origin. Moderate correlations among elements show those elements share a common source but they are not totally from the same sources, and a weak correlation among those elements indicates they have different origins. A negative correlation among elements and soil properties indicates no influence of soil properties (pH and TOC) on the distribution and concentration of trace elements in the soil.

Factor Analysis
Principal Component Analysis (PCA) is a dimensional reduction tool that is used to reduce large and complex data to a small set of variables which makes it easy for interpretation [70]. It is used to identify correlated variables having common sources [71]. PCA result for trace elements is presented in Table 6. The significant principal component is selected based on the basis of Varimax rotation of Kaiser Criteria with an eigenvalue of greater than 1 [9]. The result indicated that eigenvalues greater than 1 gave a total cumulative value of 86.029%. The variables were correlated with two principal components. Two-component factors were extracted and the result indicated that the first component with an eigenvalue of 7.167 and with a variance of 65.125% was highly correlated with the high loadings of Cr, Fe, Co, Ni, Cu and Hg (Table 6). According to PCA l values, high and strong positive loading of the elements Cr, Fe, Co, Ni and Cu was connected to anthropogenic factors. The samples were collected from different agricultural fields (paddy field, farm land, vegetable field, plastic greenhouse and barren land). Comparison between fields indicated that the highest value of Cr, Fe, Co, Ni, Cu and Hg was recorded in farmland field. This indicated that the potential sources for these elements might be associated with agricultural inputs (fertilizers, pesticides, herbicides and fungicides) and use of wastewater for irrigation. Different studies have indicated that phosphate fertilizers are highly used in China [6]. A previous study [6] reported that phosphate fertilizers are the main sources of trace elements in the environment. PCA 2 gave eigenvalue of 2.296 and % variance of 20.877 with high loading of As, Sb, Cd Zn and Pb. The PCA 2 might be connected with both anthropogenic and natural factors (parent minerals, weathering processes) and different point and nonpoint sources (application fertilizer, mining, industrial discharge, using wastewater for irrigation). A study by [19] indicated that As might have originated from parent minerals.

Cluster Analysis
Cluster Analysis (CA) is a method used to group data according to their similarity [19]. Hierarchical Cluster Analysis (HCA) was applied based on the Wards method and Euclidean distance methods [9]. The HCA result (Figure 3) indicated that elements were clustered into two major clusters. The first cluster includes Cr, Fe, Cu, Ni, Co and Hg. The second cluster includes As, Cd, Zn, Pb and Sb. HCA indicated that there is a close cluster between Cr and Fe, Ni and Cu, Co and Hg, As and Cd and Zn and Pb. This implied that those elements were from the same sources. Pearson's correlation analysis also pointed out that those elements had a strong positive correlation which suggested that the elements have the same common origin. The moderate cluster between elements ( Figure 3) indicated that these elements share similar sources.

Contamination Level and Ecological Risk of Trace Elements
Contamination factor, degree of contamination, ecological risk factor, potential ecological risk index and enrichment factor and the geo-accumulation index of trace elements were evaluated and the results presented in Table 7. The obtained result of (C i f) indicated that Co > Cd > Ni > Pb > Cu > Cr > As > Hg > Fe > Sb > Zn. According to the obtained value of (C i f), the contamination level of trace elements are grouped as low contamination from Zn (0.03), Fe (0.93) and Sb (0.35), moderate contamination from Cr, Ni, Cu, As, Hg and Pb and a high contamination of soil from Co and Cd. The degree of contamination (Cd) for elements for this the study was 1560.2, which indicates a high degree of soil contamination. Ecological risk factors (E i f) for single elements were in the decreasing order of Co > Cd > Hg > Ni > Pb > As > Cu > Sb > Cr > Zn > Fe. According to the classification given for ecological risk factor for single elements in Table 1, Cr, Ni, Cu, As, Zn, Pb, Fe and Sb had low-risk factors, Hg had a moderate risk factor and Cd and Co had an extremely high risk factors. Thus, Cd, Hg and Co can cause high risk to human and environment in Huangpi district. The same result was reported on Case label Num

Contamination Level and Ecological Risk of Trace Elements
Contamination factor, degree of contamination, ecological risk factor, potential ecological risk index and enrichment factor and the geo-accumulation index of trace elements were evaluated and the results presented in Table 7. The obtained result of (C i f ) indicated that Co > Cd > Ni > Pb > Cu > Cr > As > Hg > Fe > Sb > Zn. According to the obtained value of (C i f ), the contamination level of trace elements are grouped as low contamination from Zn (0.03), Fe (0.93) and Sb (0.35), moderate contamination from Cr, Ni, Cu, As, Hg and Pb and a high contamination of soil from Co and Cd. The degree of contamination (C d ) for elements for this the study was 1560.2, which indicates a high degree of soil contamination. Ecological risk factors (E i f ) for single elements were in the decreasing order of Co > Cd > Hg > Ni > Pb > As > Cu > Sb > Cr > Zn > Fe. According to the classification given for ecological risk factor for single elements in Table 1, Cr, Ni, Cu, As, Zn, Pb, Fe and Sb had low-risk factors, Hg had a moderate risk factor and Cd and Co had an extremely high risk factors. Thus, Cd, Hg and Co can cause high risk to human and environment in Huangpi district. The same result was reported on Cd and Hg from East Dongting and Honghu Lake in Hunan Province, China [37]. The potential risk index (RI) for this study was 5352.03, which indicated a high ecological risk due to these trace elements. The obtained enrichment factor of elements revealed that soil is enriched with Co and Cd, moderately enriched with Ni, Cu, Hg and Pb; and less enriched with Cr, Fe As and Sb (Table 7). This indicated that there is a high rate of anthropogenic disturbance in Huangpi soil. Extremely high enrichment for Cd was also reported by [72] and [73]. The obtained I geo values of the elements in soil were in decreasing order of Co > Cd > Ni >Cu > Pb > Cr > As > Hg > Fe > Sb > Zn. According to the seven classes proposed by [45], the obtained I geo result revealed that the soils under study were extremely contaminated with Cd and Co, moderately contaminated with Cr, Ni, Cu and Pb and less contaminated with As, Zn, Hg, Fe and Sb. The obtained I geo result for Cd was in line with those of [72,74,75].

Conclusions
Eleven trace elements in soil from Huangpi district were studied. All trace elements except As, Sb, Hg and Cd were detected in all samples at both soil depths. The obtained mean concentration of Cr, Ni, Cu, As, Cd, Pb were above the soil background value of Wuhan and Hubei Province. Pearson's correlation result indicated that there was a significant positive correlation among trace elements, whereas, weak and negative correlations between trace elements and soil properties (pH and TOC). PCA, HCA and EF of the soil indicated that anthropogenic factors are the major sources of trace elements in Huangpi soil. The result of contamination factor (C i f ) for trace elements were in decreasing order of Co > Cd > Ni > Pb > Cu > Cr > As > Hg > Fe > Sb > Zn. The obtained C d values indicated that there is high degree of soil contamination. The result of ecological risk factor (E i f ) of elements were in the decreasing order of Co > Cd > Hg > Ni > Pb > As > Cu > Sb > Cr > Zn > Fe. Potential ecological risk index (RI) result for the studied soil (RI > 600) indicated that there is a high ecological risk of elements in the study area. Enrichment factor (EF) showed that there is extremely high enrichment from Cd and Co, moderate from Ni, Cu, Hg and Pb, and low from Cr, Fe, As, Sb and Zn. The geo accumulation index (I geo ) results point out that the study area is highly contaminated with Cd and Co. From the obtained result it's concluded that Huangpi soil is contaminated with trace elements. Therefore, more attention should be given and remediation action should be set to minimize the concentration and ecological impacts of trace elements in the study area. Further research should be conducted on human health impact of trace element in Huangpi soil.