Next Article in Journal
Unhealthy Lifestyle and Nutritional Habits Are Risk Factors for Cardiovascular Diseases Regardless of Professed Religion in University Students
Next Article in Special Issue
Compliance with Tuberculosis Screening in Irregular Immigrants
Previous Article in Journal
Fair but Risky? Recycle Pricing Strategies in Closed-Loop Supply Chains
Previous Article in Special Issue
Relative Density of Away from Home Food Establishments and Food Spend for 24,047 Households in England: A Cross-Sectional Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan 430074, China
2
Sino-Africa Joint Research Center, Chinese Academy of Sciences, Wuhan 430074, China
3
Wuhan Botanical Garden, University of Chinese Academy of Sciences, Beijing 100049, China
4
Department of management engineering, Shangrao Vocational and Technical College, Shangrao 334109, China
5
Wuhan Britain-China School, Wuhan Foreign Languages School, Wuhan, 430022 China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(12), 2873; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15122873
Submission received: 20 November 2018 / Revised: 7 December 2018 / Accepted: 11 December 2018 / Published: 14 December 2018
(This article belongs to the Special Issue IJERPH: 15th Anniversary)

Abstract

:
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.

1. Introduction

Trace elements are identified as a dangerous group of environmental pollutants, due to their persistence, non-degradability and toxicity to living organisms [1,2]. Unlike many organic pollutants, which eventually degrade to carbon dioxide and water, trace elements tend to accumulate in the environment, especially in soil and sediments [3]. Trace elements have specific gravity greater than 5 g/cm3 and they include elements such as Cd, Zn, Pb, Hg, Sb, Cr, Co, As, Ni Cu, Mo, and Mg [4,5]. Trace elements like Cd, Cr, Pb, Hg, Ni and As have been listed as the most dangerous elements and priority control pollutants by the United States Environmental Protection Agency (USEPA) [6,7,8]. Trace elements can originate from natural sources (e.g., parent materials, weathering of rocks, volcanic eruptions and soil erosion) and anthropogenic sources (e.g., industrialization, urbanization, vehicular emissions, mining activities, smelting, burning of fossil fuels and agricultural inputs such as fertilizers, pesticides, herbicides and fungicides [9,10]. However, currently anthropogenic activities significantly accelerate the accumulation of trace elements in the environment [11,12].
Soil is one of the environmental compartments which plays an important role in plant growth, development and other ecosystem services [8,13]. On the other hand, it also serves as a sink for different environmental pollutants like toxic trace elements [14,15,16]. Several studies from different countries point out the toxic effects of trace elements on soil, plants and animals [1,3,12]. Excessive accumulation of trace elements in soil can reduce soil microorganism levels which results in soil quality degradation [16,17], and reduction of agricultural productivity [3,15,18]. It has been reported that foodstuffs and domestic feeds in Asian countries are severely contaminated with trace metals [19]. According to China national census of pollution report, more than 1.5 million sites in China have been exposed to toxic trace elements, 20 million hectares of agricultural lands have been polluted [20] and over 12 million tons of grains are contaminated by toxic elements per year [7]. About 1.5 × 105 km2 of cultivated lands is polluted by Cd [21]. In addition to the above reports, a survey conducted in 2014 indicated that 16.1% of the sampled lands were contaminated by Hg, As and Pb [6].
Plants not only absorb essential nutrients, they also absorb toxic elements from the environmental compartments [22], which results in a bio-accumulation effect [14,23]. Consequently, consumption of contaminated plants and animals can cause health problems [24]. Particularly, Cd is amongst trace metals that increased international concern due its carcinogenic effect [21]. The adverse effects of toxic trace elements on human beings have been reported from different countries of the world [24,25]. For example: As causes dermal lesions, skin cancer, peripheral vascular diseases and peripheral neuropathy [26,27], Cd can result in kidney dysfunction, hypertension, lung cancer, bone fractures, prostatic hyperplasia and adenocarcinomas, Cu can cause Alzheimer’s, prion disease [28,29], and Pb can affect the endocrine system, immune system, skeletal, circulatory system and nervous system [30].
Due to the rapid economic growth in China, there are heterogeneous anthropogenic activities which increase the accumulation of trace elements in the environment [31]. Agricultural inputs (fertilizers, herbicides, pesticides and fungicides), the use of wastewater for irrigation, urbanization, industrialization and construction of road networks are the major anthropogenic activities in Wuhan, Hubei Province [32,33,34]. Huangpi is one of the districts in Wuhan, in which the above anthropogenic activities are commonly observed [20]. However, there is lack of data and research works on the status of trace elements in Huangpi district. Previously conducted studies in Wuhan were focused only on few trace elements such as Cd, Cr, Pb Ni, As, Cu and Hg. Therefore, it’s important to investigate the status of the above listed and other additional trace elements like Co, Sb, Fe, Zn in the soil of Huangpi district. Thus, the objectives of this study were: (1) to determine the concentration level of trace elements, (2) to identify the potential sources of trace elements, (3) to evaluate pollution status and ecological risks of trace elements in the soil and (4) to determine the effect of soil pH and total organic matter on the concentration of trace elements.

2. Materials and Methods

2.1. Study Area

Soil samples were collected from Huangpi district, Wuhan, Central China. Huangpi is one of thirteen districts of Wuhan, located on northern outskirts of Wuhan between 30°52’30’’ N and 114°22’30’’ E (Figure 1). The area of the district is about 2261 km2, with a population of 1,107,565 [31]. From the total area of the land 56.12%, 18.29%, 0.22%, 7.09%, 0.21%, 23% and 4.08% are covered by cultivated land, forest land, grassland, settlements and industrial sites, land for transport, water area and unused land, respectively [35].

2.2. 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].

2.3. 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 (HNO3, 63%), 2 mL of hydrogen fluoride (HF ≥ 40%) and perchloric acid (HClO4 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).

2.4. 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 (HNO3, 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 r2 values greater than 0.994.

2.5. 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).

2.6. 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% HNO3 (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%).

2.7. Methods of Evaluating Contamination Level and Ecological Risk of Trace Elements

2.7.1. 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.
C f i = C S i / C R i
C d = Σ f i C f i
E f i = C f i T f i
R I = Σ f i / E f i
where Cif is the contamination factor for trace elements; CiS is the measured concentrations value of elements in the soil; CiR is the background reference values of trace elements in the soil; Cd is a degree of contamination; Eif is the ecological risk index of a single element; Tif is the toxicity coefficient of an element and RI represents the total potential risk index of elements [40] (Table 1).
Due to the absence of soil background reference values of trace elements in Huangpi district, Wuhan soil background values were used as a reference: Cr (90), Ni (40), Cu (35), Zn (100), As (15), Cd (0.20), Hg (0.15), Pb (35) mg/kg [42] and Hubei Province soil background values were used for Co (15.4) [32], Sb (1.65 mg/kg) [43] and Fe (29,400 mg/kg) [40]. Tif values given by [41] were Cr = 2, As = 10, Cd = 30, Ni = Cu = Pb =5, Zn =1, Hg = 40 and Co = 2, Fe = 0 and Sb = 15.

2.7.2. 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:
E F = ( C i C r e f ) s a m p l e ( B i B r e f ) b a c k g r o u n d
where Ci is the concentration of trace elements in the sample; Cref is the concentration of reference element in the sample; Bi is the background value of interest element and Bref 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.

2.7.3. Geo-accumulation Index (Igeo)

The geo-accumulation indexes (Igeo) 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:
I geo = Log 2   C n 1.5   B n
where Igeo is the geo-accumulation index, Cn is the measured concentration of an element (mg/kg) in the sample, Bn 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 Igeo: Igeo ≤ 0 (uncontaminated), 0 < Igeo ≤ 1 (uncontaminated to moderately contaminated), 1 < Igeo ≤ 2 (moderately contaminated), 2 < Igeo ≤ 3 (moderately to strongly contaminated), 3 < Igeo ≤ 4 (strongly contaminated), 4 < Igeo ≤ 5 (strongly to extremely contaminated) and Igeo > 5 (extremely contaminated) [45].

3. Result and Discussion

3.1. 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].

3.2. 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] (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 recorded mean concentration of Ni (117.80 mg/kg) was higher than the soil background value of Wuhan, Hubei province and China. In addition, the mean value of Ni (117.80 mg/kg) was above the permissible limit of FAO/ISRIC (2004) (100 mg/kg), WHO /FAO (2001), whereas below the permissible limit of USEPA (1983) (500 mg/kg). The highest value of Cu was recorded at Dujiatian (farmland), whereas the lowest was obtained from the soil of Tangjiawan site (grassland soil). According to the obtained result, the mean concentration of Cu recorded was below the permissible limit of FAO/ISRIC (2004), WHO/FAO (2001), USEPA (1983) and China (Table 2). However, the individual result of each site indicated that the concentration of Cu at Dujiatian site exceeded the permissible limit of FAO/ISRIC (2004), WHO/FAO (2001) and USEPA (1983). As compared to the other studies the mean concentration value of Cu (60.73 mg/kg)was less than the finding in Brazil (111.54 mg/kg) [62], Wuhan, China (60.85 mg/kg) [33], Hubei Province, China (386mg/kg), Democratic Republic of Congo (10,320 mg/kg) [67], Iran (100.84 mg/kg) [64], but the result was higher than the findings in soil from Cuba (43.10 mg/kg) [59], Bangladesh (20.6 mg/kg) [11], India (52.72 mg/kg) [61], Tanzania (5.62 mg/kg) [10], Pakistan (18.12 mg/kg) [63], Northern Pakistan (35.28 mg/kg) [66] and in soil along the Chao River in China (46.5 mg/kg) [65].
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].
The mean concentration of Pb was below the recommended value of FAO/ISRIC (2004) and USEPA (1983), however, above the permissible limit of WHO/FAO (2001) (Table 2). The concentration of Pb in samples from Changdi, Zhujiashan, Zhulinyuan, Lishuwan, Xinyang, Bomogang and Tianjiaxiaowan site exceeded the soil background value of China (80 mg/kg) (Table 2).

3.3. 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.

3.4. 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.

3.5. 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.

3.6. 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 (Cif) indicated that Co > Cd > Ni > Pb > Cu > Cr > As > Hg > Fe > Sb > Zn. According to the obtained value of (Cif), 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 (Eif) 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 Igeo 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 Igeo 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 Igeo result for Cd was in line with those of [72,74,75].

4. 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. The mean concentration values of Cd and Co exceed FAO/ISRIC (2004), WHO/FAO (2001) and USEPA (1983) recommended values of trace elements in soil. The mean concentration of Ni was above the permissible limit of FAO/ISRIC (2004), WHO/FAO (2001), but less than the permissible limit of USEPA (1983). 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 (Cif) for trace elements were in decreasing order of Co > Cd > Ni > Pb > Cu > Cr > As > Hg > Fe > Sb > Zn. The obtained Cd values indicated that there is high degree of soil contamination. The result of ecological risk factor (Eif) 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 (Igeo) 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.

Author Contributions

Author Contributions: Sampling, Q.X.; and X.J.; Statistical Analysis, T.G.; Writing—Original Draft Preparation, A.W.T.; Writing-Review and Editing, J.W.

Funding

This research was funded by the Sino-Africa Joint Research Center, Chinese Academy of Sciences, grant number [Y623321K01].

Acknowledgments

The authors would like to thank all members of Pollution Ecology Laboratory for their support during laboratory analysis. In addition, the authors are grateful to all friends that helped to edit final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hashmi, M.Z.; Yu, C.; Shen, H.; Duan, D.; Shen, C.; Lou, L. Risk Assessment of Heavy Metals Pollution in Agricultural Soils of Siling Reservoir Watershed in Zhejiang Province, China. Bio. Med. Res. Int. 2013. [Google Scholar] [CrossRef]
  2. Lu, C.; Wu, Y.; Hu, S.; Zhang, X.; Fu, Y. Distribution and Transport of Residual Lead and Copper Along Soil Profiles in a Mining Region of North China. Pedosphere 2016, 26, 848–860. [Google Scholar] [CrossRef]
  3. Asati, A.; Pichhode, M.; Nikhil, K. Effect of Heavy Metals on Plants: An Overview. Int. J. App. Innov. Eng. Manag. 2016, 5, 56–66. [Google Scholar]
  4. Arao, T.; Ishikawa, S.; Murakami, M.; Abe, K.; Maejima, Y.; Makino, T. Heavy metal contamination of agricultural soil and countermeasures in Japan. Paddy Water Environ. 2010, 8, 247–257. [Google Scholar] [CrossRef]
  5. Masson, P.; Dalix, T.; Bussière, S. Determination of major and trace elements in plant samples by inductively coupled plasma-mass spectrometry. Commun. Soil Sci. Plant Anal. 2010, 41, 231–243. [Google Scholar] [CrossRef]
  6. Delang, C.O. Environmental & Socio-economic Studies Causes and distribution of soil pollution in China. Environ. Socio-Econ. Studies. 2017, 5, 1–17. [Google Scholar]
  7. Xu, Y.; Liang, X.; Xu, Y.; Qin, X.; Huang, Q.; Wang, L. Remediation of Heavy Metal-Polluted Agricultural Soils Using Clay Minerals: A. Review. Pedosphere 2017, 27, 193–204. [Google Scholar] [CrossRef]
  8. Wang, M.; Zhang, H. Accumulation of Heavy Metals in Roadside Soil in Urban Area and the Related Impacting Factors. Int. J. Environ. Res. Public Health. 2018, 15, 1064. [Google Scholar] [CrossRef]
  9. Ma, L.; Yang, Z.; Li, L.; Wang, L. Source identification and risk assessment of heavy metal contaminations in urban soils of Changsha, a mine-impacted city in Southern China. Environ. Sci. Pollut. Res. 2016, 23, 17058–17066. [Google Scholar] [CrossRef] [PubMed]
  10. Kibassa, D.; Kimaro, A.A.; Shemdoe, R.S. Heavy metals concentrations in selected areas used for urban agriculture in Dar es Salaam, Tanzania. Sci. Res. Essays Acad. J. 2013, 8, 1296–1303. [Google Scholar]
  11. Rahman, M.S.; Biswas, P.K.; Al Hasan, S.M.; Rahman, M.M.; Lee, S.H.; Kim, K.H.; Rahman, S.M.; Islam, M.R. The occurrences of heavy metals in farmland soils and their propagation into paddy plants. Environ. Monit. Assess. 2018. [Google Scholar] [CrossRef] [PubMed]
  12. Fosu-Mensah, B.Y.; Addae, E.; Yirenya-Tawiah, D.; Nyame, F. Heavy metals concentration and distribution in soils and vegetation at Korle Lagoon area in Accra, Ghana. Cogent. Environ. Sci. 2017, 3, 1–14. [Google Scholar] [CrossRef]
  13. Marrugo-negrete, J.; Pinedo-hernández, J.; Díez, S. Assessment of heavy metal pollution, spatial distribution and origin in agricultural soils along the Sinú River Basin, Colombia. Environ. Res. 2017, 154, 380–388. [Google Scholar] [CrossRef] [PubMed]
  14. Fernandes Azevedo, B.; Barros Furieri, L.; Peçanha, F.M.; Wiggers, G.A.; Frizera Vassallo, P.; Ronacher Simões, M. Toxic Effects of Mercury on the Cardiovascular and Central Nervous Systems. J. Biomed. Biotechnol. 2012, 2012, 1–11. [Google Scholar] [CrossRef] [PubMed]
  15. Chauhan, G.; Chauhan, P.U.K. Risk Assessment of Heavy Metal Toxicity Through Contaminated Vegetables From Waste. Int. J. Adv. Technol. Eng. Sci. 2014, 51, 444–460. [Google Scholar]
  16. Olatunji, O.S.; Opeolu, B.O.; Fatoki, O.S.; Ximba, B.J. Heavy metals concentration levels in selected arable agricultural soils in South Western Nigeria. Int. J. Phys. Sci. 2013, 8, 421–427. [Google Scholar]
  17. Chang, C.Y.; Yu, H.Y.; Chen, J.J.; Li, F.B.; Zhang, H.H.; Liu, C.P. Accumulation of heavy metals in leaf vegetables from agricultural soils and associated potential health risks in the Pearl River Delta, South China. Environ. Monit. Assess. 2014, 186, 1547–1560. [Google Scholar] [CrossRef]
  18. Tangahu, B.V.; Sheikh Abdullah, S.R.; Basri, H.; Idris, M.; Anuar, N.; Mukhlisin, M. A review on heavy metals (As, Pb, and Hg) uptake by plants through phytoremediation. Int. J. Chem. Eng. 2011. [Google Scholar] [CrossRef]
  19. Yuan, G.L.; Sun, T.H.; Han, P.; Li, J.; Lang, X.X. Source identification and ecological risk assessment of heavy metals in topsoil using environmental geochemical mapping: Typical urban renewal area in Beijing, China. J. Geochem. Explor. 2014, 136, 40–47. [Google Scholar] [CrossRef]
  20. Hu, H.; Jin, Q.; Kavan, P. A Study of Heavy Metal Pollution in China: Current Status, Pollution-control Policies and Countermeasures. Sustainability 2014, 6, 5820–5838. [Google Scholar] [CrossRef]
  21. Yu, H.; Wang, J.; Fang, W.; Yuan, J.; Yang, Z. Cadmium accumulation in different rice cultivars and screening for pollution-safe cultivars of rice. Sci. Total Environ. 2006, 370, 302–309. [Google Scholar] [CrossRef] [PubMed]
  22. Boke, A.; Megersa, N.; Teju, E. Quantitative Determination of the Heavy Metal Levels in the Wild Edible Plant Parts and their Corresponding Soils of the Central and Western Regions of the Oromia State, Ethiopia. J. Environ. Anal. Toxicol. 2015. [Google Scholar] [CrossRef]
  23. Kim, R.Y.; Yoon, J.K.; Kim, T.S.; Yang, J.E.; Owens, G.; Kim, K.R. Bioavailability of heavy metals in soils: Definitions and practical implementation—A critical review. Environ. Geochem. Health 2015, 37, 1041–1061. [Google Scholar] [CrossRef] [PubMed]
  24. Ogunlaja, O.O.O.; Moodley, R.; Baijnath, H.; Jonnalagadda, S.B. Nutritional evaluation, bioaccumulation and toxicological assessment of heavy metals in edible fruits of FicussurForssk (Moraceae). J. Environ. Sci. Health 2017, 52, 84–91. [Google Scholar] [CrossRef] [PubMed]
  25. Hussain, A.; Abbas, N.; Arshad, F.; Akram, M.; Khan, Z.I.; Ahmad, K. Effects of diverse doses of Lead (Pb) on different growth attributes of Zea-Mays L. Agric. Sci. 2013, 4, 262–265. [Google Scholar]
  26. Barrachina, A.C.; Carbonell, F.B.; Beneyto, J.M. Arsenic uptake, distribution, and accumulation in tomato plants: Effect of arsenite on plant growth and yield. J. Plant Nutr. 1995, 18, 1237–1250. [Google Scholar] [CrossRef]
  27. Kamunda, C.; Mathuthu, M.; Madhuku, M. Health risk assessment of heavy metals in soils from witwatersrand gold mining basin, South Africa. Int. J. Environ. Res. Public Health 2016. [Google Scholar] [CrossRef]
  28. Wang, M.Y.; Chen, A.K.; Wong, M.H.; Qiu, R.L.; Cheng, H.; Ye, Z.H. Cadmium accumulation in and tolerance of rice (Oryza sativa L.) varieties with different rates of radial oxygen loss. Environ. Pollut. 2011, 159, 1730–1736. [Google Scholar] [CrossRef]
  29. Lin, Y.S.; Caffrey, J.L.; Lin, J.W.; Bayliss, D.; Faramawi, M.F.; Bateson, T.F. Increased risk of cancer mortality associated with cadmium exposures in Older Americans with low Zinc Intake. J. Toxicol. Environ. Health 2013, 76, 1–15. [Google Scholar] [CrossRef]
  30. Nieboer, E.; Tsuji, L.J.S.; Martin, I.D.; Liberda, E.N. Human biomonitoring issues related to lead exposure. Environ. Sci. Process. Impacts 2013, 15, 1824. [Google Scholar] [CrossRef]
  31. Hagberg, A. Industrial Wastewater Treatment and other Environmental Problems in Wuhan—Is Swedish Technology a Solution? Uppsala University: Uppsala, Sweden, 2007. [Google Scholar]
  32. Gong, M.; Wu, L.; Bi, X.Y.; Ren, L.M.; Wang, L.; Ma, Z.D.; Bao, Z.Y.; Li, Z.G. Assessing heavy-metal contamination and sources by GIS-based approach and multivariate analysis of urban-rural topsoils in Wuhan, central China. Environ. Geochem. Health 2010, 32, 59–72. [Google Scholar] [CrossRef] [PubMed]
  33. Qu, M.K.; Li, W.D.; Zhang, C.R.; Wang, S.Q.; Yang, Y.; He, L.Y. Source Apportionment of Heavy Metals in Soils Using Multivariate Statistics and Geostatistics. Pedosphere 2013, 23, 437–444. [Google Scholar] [CrossRef]
  34. Li, Z.; Zhang, G.; Liu, Y.; Wan, K.Y.; Zhang, R.; Chen, F. Soil Nutrient Assessment for Urban Ecosystems in Hubei, China. PLoS ONE 2013, 8, e75856. [Google Scholar] [CrossRef] [PubMed]
  35. Lv, J.; Ma, T.; Dong, Z.; Yao, Y. Temporal and Spatial Analyses of the Landscape Pattern of Wuhan City Based on Remote Sensing Images. Int. J. Geo-Inf. 2018, 7, 340. [Google Scholar] [CrossRef]
  36. Sekabira, K.; Oryem Origa, H.; Basamba, T.A.; Mutumba, G.; Kakudidi, E. Assessment of heavy metal pollution in the urban stream sediments and its tributaries. Int. J. Environ. Sci. Tech. 2010, 7, 435–446. [Google Scholar] [CrossRef] [Green Version]
  37. Makokha, V.A.; Qi, Y.; Shen, Y.; Wang, J. Concentrations, Distribution, and Ecological Risk Assessment of Heavy Metals in the East Dongting and Honghu Lake, China. Expo. Health 2016, 8, 31–41. [Google Scholar] [CrossRef]
  38. Weldegebriel, Y.; Chandravanshi, B.S.; Wondimu, T. Concentration levels of metals in vegetables grown in soils irrigated with river water in Addis Ababa, Ethiopia. Ecotoxicol. Environ. Saf. 2012, 77, 57–63. [Google Scholar] [CrossRef] [PubMed]
  39. Wlakey, A. A critical examination of a rapid method for determination of organic carbon in soils—Effect of variations in digestion conditions and of inorganic soil constituents. Soil Sci. 1947, 63, 251–257. [Google Scholar]
  40. Liu, H.; Liu, G.; Zhou, C.; Yuan, Z.; Da, C. Geochemical speciation and ecological risk assessment of heavy metals in surface soils collected from the Yellow River Delta National Nature Reserve, China. Hum. Ecol. Risk Assess. 2017, 23, 1585–1600. [Google Scholar] [CrossRef]
  41. Hakanson, L. An ecological risk index for aquatic pollution control.a sedimentological approach. Water Res. 1980, 14, 975–1001. [Google Scholar] [CrossRef]
  42. Yang, Z.; Wang, Y.; Shen, Z.; Niu, J.; Tang, Z. Distribution and speciation of heavy metals in sediments from the mainstream, tributaries, and lakes of the Yangtze River catchment of Wuhan, China. J. Hazard Mater. 2009, 166, 1186–1194. [Google Scholar] [CrossRef]
  43. He, M.; Wang, X.; Wu, F.; Fu, Z. Antimony pollution in China. Sci. Total Environ. 2012, 421–422, 41–50. [Google Scholar] [CrossRef] [PubMed]
  44. Jiao, X.; Teng, Y.; Zhan, Y.; Wu, J.; Lin, X. Soil heavy metal pollution and risk assessment in Shenyang industrial district, Northeast China. PLoS ONE 2015, 10, e0127736. [Google Scholar] [CrossRef] [PubMed]
  45. Muller, H.G. Mechanical Properties, Rheology, and Haptaesthesis of Food. J. Texture Stud. 1969, 1, 38–42. [Google Scholar] [CrossRef] [PubMed]
  46. Ogundiran, M.B.; Osibanjo, O. Mobility and speciation of heavy metals in soils impacted by hazardous waste. Chem. Speciat. Bioavailab. 2009, 21, 59–69. [Google Scholar] [CrossRef] [Green Version]
  47. Liu, L.; Lu, J.; Zhang, Z.; Zheng, H.; Gao, X.; Zhang, W. Heavy Metals Contamination in Greenhouse Soils and Vegetables in Guanzhong, China. J. Encapsulation Adsorption Sci. 2014, 4, 80–88. [Google Scholar] [CrossRef]
  48. Bai, Y.; Wang, M.; Peng, C.; Alatalo, J.M. Impacts of urbanization on the distribution of heavy metals in soils along the Huangpu River, the drinking water source for Shanghai. Environ. Sci. Pollut. Res. 2016, 23, 5222–5231. [Google Scholar] [CrossRef] [PubMed]
  49. Jiang, X.; Lu, W.X.; Zhao, H.Q.; Yang, Q.C.; Yang, Z.P. Potential ecological risk assessment and prediction of soil heavy-metal pollution around coal gangue dump. Nat. Hazards Earth Syst. Sci. 2014, 14, 1599–1610. [Google Scholar] [CrossRef]
  50. Fang, C.; Moncrieff, J.B. The variation of soil microbial respiration with depth in relation to soil carbon composition. Plant Soil 2005, 268, 243–253. [Google Scholar] [CrossRef]
  51. Celik, I. Land-use effects on organic matter and physical properties of soil in a southern Mediterranean highland of Turkey. Soil Tillage Res. 2005, 83, 270–277. [Google Scholar] [CrossRef]
  52. Behbahaninia, A.; Mirbagheri, S.A.; Khorasani, N.; Nouri, J.; Javid, A.H. Heavy metal contamination of municipal effluent in soil and plants. J. Food Agric. Environ. 2009, 7, 851–856. [Google Scholar]
  53. Rahaman, A.; Sadia Afroze, J.; Bashar, K.; Farhad Ali, M.; Razib Hosen, M. A Comparative Study of Heavy Metal Concentration in Different Layers of Tannery Vicinity Soil and Near Agricultural Soil. Am. J. Anal. Chem. 2016, 7, 880–889. [Google Scholar] [CrossRef]
  54. Li, J.; Lu, Y.; Yin, W.; Gan, H.; Zhang, C.; Deng, X.; Lian, J. Distribution of heavy metals in agricultural soils near a petrochemical complex in Guangzhou, China. Environ. Monit. Assess. 2009, 153, 365–375. [Google Scholar] [CrossRef] [PubMed]
  55. Camobreco, V.J.; Richards, B.K.; Steenhuis, T.S.; Peverly, J.H.; McBride, M.B. Movement of heavy metals through undisturbed and homogenized soil columns. Soil Sci. 1996, 161, 740–750. [Google Scholar] [CrossRef]
  56. Akenga, T.; Sudoi, V.; Machuka, W.; Kerich, E. Heavy Metal Concentrations in Agricultural Farms in Homa Hills Homa Bay County, Kenya. Int. J. Sci. Res. 2016, 5, 1664–1669. [Google Scholar]
  57. Yabe, J.; Ishizuka, M.; Umemura, T. Current Levels of Heavy Metal Pollution in Africa. J. Vet. Med. Sci. 2010, 72, 1257–1263. [Google Scholar] [CrossRef] [Green Version]
  58. Mungai, T.M.; Owino, A.A.; Makokha, V.A.; Gao, Y.; Yan, X.; Wang, J. Occurrences and toxicological risk assessment of eight heavy metals in agricultural soils from Kenya, Eastern Africa. Environ. Sci. Pollut. Res. 2016, 23, 18533–18541. [Google Scholar] [CrossRef]
  59. Accioly, A.; Montero, A.; Ugarte, O.M.; do Nascimento, C.W.; de Aguiar Accioly, A.M.; Biondi, C.M.; da Silva, Y.J. Background concentrations and reference values for heavy metals in soils of Cuba. Environ. Monit. Assess. 2015, 187, 4198. [Google Scholar]
  60. Rahman, S.H.; Khanam, D.; Adyel, T.M.; Islam, M.S.; Ahsan, M.A.; Akbor, M.A. Assessment of Heavy Metal Contamination of Agricultural Soil around Dhaka Export Processing Zone (DEPZ), Bangladesh: Implication of Seasonal Variation and Indices. Appl. Sci. 2012, 2, 584–601. [Google Scholar] [CrossRef] [Green Version]
  61. Rakesh Sharma, M.S.; Raju, N.S. Correlation of Heavy Metal contamination with Soil properties of Industrial areas of Mysore, Karnataka, India by Cluster analysis. Int. Res. J. Environ. Sci. 2013, 2, 22–27. [Google Scholar]
  62. De Almeida Júnior, A.B.; do Nascimento, C.W.A.; Biondi, C.M.; de Souza, A.P.; Barros, F.M.; do Rêgo Barros, F.M. Background and reference values of metals in soils from Paraíba state, Brazil. Rev. Bras. Cienc. Sol. 2016, 40, 1–13. [Google Scholar] [CrossRef]
  63. Tariq, S.R.; Shafiq, M.; Chotana, G.A. Distribution of Heavy Metals in the Soils Associated with the Commonly Used Pesticides in Cotton Fields. Scientifica 2016. [Google Scholar] [CrossRef] [PubMed]
  64. Soltani, N.; Keshavarzi, B.; Moore, F.; Sorooshian, A.; Ahmadi, M.R. Distribution of potentially toxic elements (PTEs) in tailings, soils, and plants around Gol-E.-Gohar iron mine, a case study in Iran. Environ. Sci. Pollut. Res. 2017, 24, 18798–18816. [Google Scholar] [CrossRef] [PubMed]
  65. Qin, F.; Ji, H.; Li, Q.; Guo, X.; Tang, L.; Feng, J. Evaluation of trace elements and identification of pollution sources in particle size fractions of soil from iron ore areas along the Chao River. J. Geochem. Explor. 2014, 138, 33–49. [Google Scholar] [CrossRef]
  66. Rehman, I.U.; Ishaq, M.; Ali, L.; Khan, S.; Ahmad, I.; Din, I.U. Enrichment, spatial distribution of potential ecological and human health risk assessment via toxic metals in soil and surface water ingestion in the vicinity of Sewakht mines, district Chitral, Northern Pakistan. Ecotoxicol. Environ. Saf. 2018, 154, 127–136. [Google Scholar] [CrossRef] [PubMed]
  67. Pourret, O.; Lange, B.; Bonhoure, J.; Colinet, G.; Decrée, S.; Mahy, G.; Séleck, M.; Shutcha, M.; Faucon, M.-P. Assessment of soil metal distribution and environmental impact of mining in Katanga (Democratic Republic of Congo). Appl. Geochem. 2016, 64, 43–55. [Google Scholar] [CrossRef]
  68. Gao, L.; Gao, B.; Zhou, Y.; Xu, D.; Sun, K. Predicting remobilization characteristics of cobalt in riparian soils in the Miyun Reservoir prior to water retention. Ecol. Indic. 2017, 80, 196–203. [Google Scholar] [CrossRef]
  69. Khaledian, Y.; Pereira, P.; Brevik, E.C.; Pundyte, N.; Paliulis, D. The Influence of Organic Carbon and pH on Heavy Metals, Potassium, and Magnesium Levels in Lithuanian Podzols. Land Degrad. Dev. 2017, 28, 345–354. [Google Scholar] [CrossRef]
  70. Constantin, C. Principal Component Analysis—A Powerful Tool in Computing Marketing Information. Bull. Transilv. Univ. Brasov Ser. V Econ. Sci. 2014, 7, 25–30. [Google Scholar]
  71. Yang, J.; Teng, Y.; Song, L.; Zuo, R. Tracing sources and contamination assessments of heavy metals in road and foliar dusts in a typical mining city, China. PLoS ONE 2016, 11, e0168528. [Google Scholar] [CrossRef]
  72. Abdulqaderismaeel, W.; Kusag, A. Enrichment Factor and Geo-accumulation Index for Heavy Metals at Industrial Zone in Iraq. IOSR J. Appl. Geol. Geophys. 2015, 3, 2321–2990. [Google Scholar]
  73. Zhang, Z.; Juying, L.; Mamat, Z. Sources identification and pollution evaluation of heavy metals in the surface sediments of Bortala River, Northwest China. Ecotoxicol. Environ. Saf. 2016, 126, 94–101. [Google Scholar] [CrossRef] [PubMed]
  74. Gupta, S.; Jena, V.; Matic, N.; Kapralova, V.; Solanki, J.S. Assessment of Geo-Accumulation Index of Heavy Metal and Source of Contamination By Multivariate Factor Analysis. Int. J. Hazard Mater. 2014, 2, 18–22. [Google Scholar]
  75. Alshahri, F.; El-Taher, A. Assessment of Heavy and Trace Metals in Surface Soil Nearby an Oil Refinery, Assessment of Heavy and Trace Metals in Surface Soil Nearby an Oil Refinery, Saudi Arabia, Using Geoaccumulation and Pollution Indices. Arch. Environ. Contam. Toxicol. 2018, 75, 390–401. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Study area. (The number (1–18) represents sampling sites: S1 (Tangjiawan), S2 (Fengdouhu) S3 (Erpaiqu), S4 (Changdi), S5 (Zhujiashan), S6 (Tujiadun), S7 (Zhulinyuan). S8 (Zhoujiawan), S9 (Lishuwan), S10 (Xinyang), S11 (Leqianwan), S12 (Bomogang), S13 (Hanjiafan), S14 (Wanjiatian), S15 (Hongguanshanxiawan), S16 (Dujiatian), S17 (Tianjiaxiaowan), S18 (Tianjiaxiaowan).
Figure 1. Study area. (The number (1–18) represents sampling sites: S1 (Tangjiawan), S2 (Fengdouhu) S3 (Erpaiqu), S4 (Changdi), S5 (Zhujiashan), S6 (Tujiadun), S7 (Zhulinyuan). S8 (Zhoujiawan), S9 (Lishuwan), S10 (Xinyang), S11 (Leqianwan), S12 (Bomogang), S13 (Hanjiafan), S14 (Wanjiatian), S15 (Hongguanshanxiawan), S16 (Dujiatian), S17 (Tianjiaxiaowan), S18 (Tianjiaxiaowan).
Ijerph 15 02873 g001
Figure 2. The mean concentrations of trace elements and selected soil properties (pH and TOC) at soil depth of (1–10 and 10–20 cm).
Figure 2. The mean concentrations of trace elements and selected soil properties (pH and TOC) at soil depth of (1–10 and 10–20 cm).
Ijerph 15 02873 g002
Figure 3. Dendrogram of hierarchal cluster analysis for trace elements in soil using Ward methods.
Figure 3. Dendrogram of hierarchal cluster analysis for trace elements in soil using Ward methods.
Ijerph 15 02873 g003
Table 1. The standard given for the contamination factor, degree of contamination, ecological risk factor and risk index of elements.
Table 1. The standard given for the contamination factor, degree of contamination, ecological risk factor and risk index of elements.
CifContamination LevelCd ClassDegree
Contamination Level
EifPollution DegreeRIRisk Degree
Cif < 1Low contamination factorCd < 8Low Eif < 40Low riskRI < 150Low ecological risk
1 ≤ Cif < 3Moderate contamination8 ≤ Cd < 16Moderate 40 ≤ Eif ≤ 80Moderate risk150 ≤ RI < 300Moderate ecological risk
3 ≤ Cif < 6Considerable contamination factor 16 ≤ Cd < 32Considerable80 ≤ Eif <160Considerable risk300 ≤ RI < 600Considerable ecological risk
Cif ≥ 6Very high contaminationCd ≥ 32Very high160 ≤ Eif < 320high riskRI > 600Very strong
Eif ≥ 320Extremely high
Source: [41].
Table 2. Descriptive statistics of trace elements (mg/kg) and selected soil properties (pH and TOC) of soil from Huangpi district and background value of Wuhan, Hubei and China.
Table 2. Descriptive statistics of trace elements (mg/kg) and selected soil properties (pH and TOC) of soil from Huangpi district and background value of Wuhan, Hubei and China.
Elements MinMaxMeanSD Skewness(a)(b)(c)(d)(e)(f)
Cr66.56321.73140.158.841.7959086200250-1000
Fe13,583.0455,398.0127,304.910,705.11.148-29,40029,400---
Co7,244.465,4621.9122,656.9410,317.81.578-15.4405050-
Ni51.18210.63117.850.930.7014037.350100100500
Cu26.09139.9860.7330.061.1973530.7100100100100
AsND47.5815.5817.680.8261512.330--75
SbND1.540.580.420.662-1.6510---
Cd0.0777.6215.4423.841.690.20.1720.5530.7
Zn1.534.813.320.94−0.54410083.6250500300300
HgND1.130.150.263.4580.150.080.7-21
Pb38.81117.974.1620.220.3953526.78015050200
pH4.26.875.710.73−0.543-6.5----
TOC0.652.411.720.47−0.815------
ND: Not detected; Max = Maximum; Min=minimum; SD = Standard deviation. (a), Wuhan soil background value [42], (b), Hubei province soil background [37], (c), china soil background value [27], (d), Food and Agriculture Organization (FAO)/International Soil Reference and Information Centre (ISRIC) (2004) [57] (e), World Health Organization (WHO)/Food and Agricultural Organization (FAO) (2001) [12], (f), United States and Environmental Protection Agency (USEPA) (1983) [58].
Table 3. The concentration of trace elements (mg/kg), soil pH and TOC for each sample at soil depths of (1–10 and 10–20 cm).
Table 3. The concentration of trace elements (mg/kg), soil pH and TOC for each sample at soil depths of (1–10 and 10–20 cm).
SCrFeCoNiCuAsSbCdZnHgPbPHTOC
Depth of 1–10 cm
194.2014,491.110,133.973.6439ND0.58ND2.850.1553.425.211.59
2169.2028,732.122,910.7163.310931.781.0934.614.350.3670.715.911.75
3163.1328,633.923,232.1160.86122.110.6313.824.380.5177.957.861.71
4251.5440,557.733,375293.7512342.931.7575.135.650.19106.448.072.6
5266.2546,817.336,740.4253.0813863.831.95124.596.570.58151.065.392.14
6112.6920,019.216,653.997.662ND0.2ND1.970.1961.885.731.51
7224.5039,65031,433.318116560.821.78108.465.40.31141.336.291.03
8137.7723,946.435,053.6120.273343.640.6528.454.28ND75.776.570.63
9124.0421,75022,846.296.92931.010.70.714.13ND94.176.142.58
10130.1720,17520,925120.08533.240.521.673.670.08106.756.582.31
11138.5628,903.924,884.6100.194310.020.525.22.58ND74.235.861.14
12149.0229,151.829,098.2122.866118.280.7112.894.030.08101.074.792.23
1388.5115,87510,846.257.952ND0.04ND2.75ND65.776.482.25
14172.5029,033.328,258.3125.4260ND0.71ND3.47ND83.055.761.71
15110.3623,071.415,116.170.31361.3NDND2.65ND65.434.61.46
16328.6651,901.855,580.4203.7513123.030.133.083.731.4157.65.051.4
17138.1720,908.325,366.7109.33475.520.874.784.930.3387.174.121.55
18107.0516,642.919,991.18135ND0.56ND3.48ND74.85.331.68
Depth of 10–20 cm
138.9312,675435528.71130.690.220.151.13ND24.26.012
2165.1746,416.723,508.3207.676929.741.0330.733.84ND77.345.891.85
39526,758.313,491.787.085115.430.365.532.6ND40.595.461.81
4125.4233,083.317,733.3127.54230.590.3911.982.27ND55.285.661.5
5145.0936,517.921,232.1138.935727.241.1330.653.06ND65.515.472
661.4716,57510,041.748.78185.360.080.41.09ND30.766.071.42
7130.6334,928.620,723.2107.233034.340.5318.092.820.0594.465.90.93
8118.7531,00028,875100.84244.490.5825.823.20.1879.415.660.66
999.1320,884.617,682.782.6946ND0.45ND2.54ND67.796.32.24
10113.320,232.117,598.2101.3442ND0.29ND3.16ND78.186.422.17
1170.0817,044.611,857.165.8712NDND0.741.01ND39.46.521.08
12106.1625,50021,00090.363413.60.395.342.47ND62.594.41.94
1365.6213,182.78,005.7749.5457NDND2.522.38ND49.655.362.04
14135.4523,392.922,187.5106.52502.730.140.392.64ND58.495.351.99
1593.7119,682.712,990.459.9123NDND0.291.93ND52.134.651.71
16314.8158,894.253,663.5207.0214916.050.142.314.080.8670.45.152.71
17112.2320,446.421,982.193.394312.540.567.063.71ND78.784.281.46
18146.5225,50026,276.8106.52660.541.080.594.790.2696.255.381.28
S: Sample, ND: Not detected.
Table 4. Comparison of the median concentrations (mg/kg) result of trace elements for this study with the other studies.
Table 4. Comparison of the median concentrations (mg/kg) result of trace elements for this study with the other studies.
Trace Elements
CountriesCrFeCoNiCuAsSbCdZnHgPbReference
China (Wuhan)140.1027,304.922,656.94119.1260.7315.580.5815.443.320.1574.16This study
Bangladesh-18,000-24.4820.068.34----0.85[11]
Cuba85.9-9.1669.5743.120.32-0.52100.295.414.22[59]
Bangladesh53.730,404-48.1604073.1-0.0072209486.649.66[60]
China (Xihu district)53.3-7.3222.938.7--0.387139-70[8]
China (Hubei)----38635.4-2.59--120[34]
India8.0132.12-10.8652.72---44.72--[61]
Tanzania 7.68---5.62--0.2233.18-14.32[10]
Brazil 20.6120,273.757.4412.86111.54-13.8138.31224.29954.88-[62]
Pakistan 5.86-7.5622.1618.12--0.59--16.18[64]
Wuhan (China)85-163434--0.2850.1133[32]
Wuhan (China) 152.78-16.3752.8760.85--3.9886.4-30.17[33]
Democratic republic of Congo-4.59902010,32029-49726-135[67]
Iran48.08101,588.8938.574.69100.84160.220.1672.9622.3010.80[64]
Iran53.2133,428.1316.5171.5638.9511.770.210.2271.9131.4515.07[64]
Soil along Chao River (China)118-17.52046.56.07-0.161130.36020.3[65]
Northern Pakistan29.94-36.7626.6135.28--2.04101.76 4.69[66]
Miyun Reservoir (China)--10.7–38.74--------[68]
Table 5. Pearson’s correlation coefficient between individual trace elements and soil properties.
Table 5. Pearson’s correlation coefficient between individual trace elements and soil properties.
ElementsCrFeCoNiCuAsSbCdZnHgPbpHTOC
Cr1
Fe0.961 **1
Co0.920 **0.861 **1
Ni0.874 **0.917 **0.730 **1
Cu0.904 **0.895 **0.750 **0.838 **1
As0.535 *0.692 **0.489 *0.725 **0.544 *1
Sb0.340.4450.2350.634 **0.4260.733 **1
Cd0.4220.576 *0.2730.662 **0.514 *0.902 **0.854 **1
Zn0.618 **0.609 **0.599 **0.719 **0.631 **0.635 **0.770 **0.609 **1
Hg0.849 **0.772 **0.809 **0.588 *0.763 **0.269-0.0030.10.3631
Pb0.3880.4120.3750.495 *0.4480.592 **0.759**0.694 **0.786 **0.0221
pH−0.070.008−0.1650.1730.070.2220.1390.202−0.072−0.1820.0171
TOC0.1480.059−0.0240.1750.292−0.318−0.048−0.1620.0870.118−0.0260.0071
** Correlation is significant at the 0.01 level (2-tailed); * Correlation is significant at the 0.05 level (2-tailed).
Table 6. Rotated component matrix of trace elements in soil.
Table 6. Rotated component matrix of trace elements in soil.
ElementsComponent
12
Cr0.9490.282
Fe0.8860.408
Co0.9040.198
Ni0.7370.579
Cu0.8460.377
As0.3570.805
Sb0.0770.949
Cd0.1780.911
Zn0.4380.748
Hg0.948−0.103
Pb0.1300.848
Eigenvalues % of Variance
Cumulative %
7.167
65.152
65.152
2.296
20.877
86.029
Extraction Method: Principal Component Analysis. Rotation Method: Varimax with Kaiser Normalization.
Table 7. Contamination Factor (Cif), Degree of Contamination (Cd), Ecological Risk Factor (Eif), Risk Index (RI), Enrichment Factor (EF) and Geo-accumulation Index (Igeo) of trace elements.
Table 7. Contamination Factor (Cif), Degree of Contamination (Cd), Ecological Risk Factor (Eif), Risk Index (RI), Enrichment Factor (EF) and Geo-accumulation Index (Igeo) of trace elements.
ElementsCifEifEFDegree of EFIgeoContamination Level
Cr1.563.111.75Low enrichment0.05Uncontaminated to moderately contaminated
Fe0.9301Low enrichment−0.69Uncontaminated
Co1471.232942.461584.12Extremely high enrichment9.94Extremely contaminated
Ni2.9814.893.44Moderate enrichment0.99Uncontaminated to moderately contaminated
Cu1.748.682.13Moderate enrichment0.21Uncontaminated to moderately contaminated
As1.0410.391.36Low enrichment−0.53Uncontaminated
Sb0.355.270.38Low enrichment−2.09Uncontaminated
Cd77.222316.6190.4Extremely high enrichment5.69Extremely contaminated
Zn0.030.030.04Low enrichment−5.5Uncontaminated
Hg1402.02Moderate enrichment−0.58Uncontaminated
Pb2.1210.592.99Moderate enrichment0.5Uncontaminated to moderately contaminated
     Cd = ΣCif =1560.2, RI = ΣEif = 5352.03

Share and Cite

MDPI and ACS Style

Tadesse, A.W.; Gereslassie, T.; Xu, Q.; Tang, X.; Wang, J. Concentrations, Distribution, Sources and Ecological Risk Assessment of Trace Elements in Soils from Wuhan, Central China. Int. J. Environ. Res. Public Health 2018, 15, 2873. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15122873

AMA Style

Tadesse AW, Gereslassie T, Xu Q, Tang X, Wang J. Concentrations, Distribution, Sources and Ecological Risk Assessment of Trace Elements in Soils from Wuhan, Central China. International Journal of Environmental Research and Public Health. 2018; 15(12):2873. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15122873

Chicago/Turabian Style

Tadesse, Ababo Workineh, Tekleweini Gereslassie, Qiang Xu, Xiaojun Tang, and Jun Wang. 2018. "Concentrations, Distribution, Sources and Ecological Risk Assessment of Trace Elements in Soils from Wuhan, Central China" International Journal of Environmental Research and Public Health 15, no. 12: 2873. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15122873

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop