Next Article in Journal
Toward an Improved Air Pollution Warning System in Quebec
Previous Article in Journal
The Potential of Remedial Techniques for Hazard Reduction of Steel Process by Products: Impact on Steel Processing, Waste Management, the Environment and Risk to Human Health
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heavy Metals in Agricultural Soils of the Lihe River Watershed, East China: Spatial Distribution, Ecological Risk, and Pollution Source

1
Guangzhou Marine Geological Survey, 477 Huanshi East Road, Guangzhou 510075, China
2
School of Geography and Ocean Science, Nanjing University, 163 Xianlin Road, Nanjing 210023, China
3
College of Forestry, Nanjing Forestry University, 159 Longpan Road, Nanjing 210037, China
4
Institute of Land and Urban-Rural Development, Zhejiang University of Finance & Economics, Hangzhou 310018, China
5
College of Materials and Environmental Engineering, Hangzhou Dianzi University, Hangzhou 310018, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(12), 2094; https://doi.org/10.3390/ijerph16122094
Submission received: 29 March 2019 / Revised: 31 May 2019 / Accepted: 4 June 2019 / Published: 13 June 2019

Abstract

:
Concentrations of cadmium, chromium, copper, nickel, lead, and zinc in agricultural soils at 32 sites in the Lihe River Watershed of the Taihu region, East China, and their potential ecological risks and possible sources were investigated. Enrichment factor analysis demonstrated enrichment in the order Cd > Pb > Zn > Cu > Ni > Cr. The potential ecological risk index and risk assessment code analyses indicated that, of the metals studied, Cd posed the most significant ecological risk in the study area. Statistical analyses, GIS mapping, and enrichment factor analysis suggested that Cd, Pb, Cu, and Zn were derived mainly from anthropogenic sources, including agricultural, industrial, and vehicular emissions, while Cr and Ni were mainly from natural sources. Positive matrix factorization revealed that Cd, Cr, Cu, Ni, Pb, and Zn were sourced from industrial and vehicular emissions (73.7%, 21.3%, 71.4%, 20.3%, 75.0%, and 62.2%, respectively), the agricultural sector (26.3%, 36.3%, 6.8%, 38.9%, 15.7%, and 6.9%, respectively), and parent materials (0%, 42.4%, 21.8%, 40.8%, 9.2%, and 30.9%, respectively). It was recommended that strategies be implemented to reduce industrial point-source pollution.

Graphical Abstract

1. Introduction

Rapid urbanization and industrialization over recent decades have resulted in heavy metal contamination of agricultural soils, and heavy metal pollution of agricultural soil has become a worldwide environmental concern [1,2,3,4,5]. Heavy metal accumulation in soils can significantly affect human health through their environmental persistence, bioaccumulation, and transfer through food chains; consequently, control of pollution sources is required to minimize their impact [6]. The effective assessment of ecological risks and source apportionment is a crucial step in this regard.
Numerous analytical techniques have been applied in assessing ecological risks posed by heavy metals in agricultural soils, based on their concentration, distribution, and speciation [7,8,9]. Wang et al. [9] used the potential ecological risk index (PERI) to quantify degree of metal pollution and assess the ecological risk of heavy metals in soils. Although the PERI considers the toxicity and total concentrations of heavy metals, a degree of subjective assessment is involved and chemical speciation is ignored [8,10,11], even though toxicity and mobility are strongly dependent on the latter [12,13]. The risk assessment code (RAC) was introduced because of the advantage of comprehensive analysis of heavy metals and other components in soil, providing a better indication of relationships between bioavailability and mobility, as well as environmental risk [8,14]. Until now, there has been no systematic or integrated research focusing on ecological risk assessment of heavy metal contamination in the Lihe River Watershed in East China. In addition, the relevance of acid-soluble fractions of heavy metals, and source apportionment in this region, have not yet been investigated.
It is generally considered that heavy metals originate from two primary sources: natural (e.g., weathering of parental material) and anthropogenic (e.g., metalliferous industries, mining, vehicle exhausts, and agronomic practices) [1,2,15,16]. Existing methods of heavy metal source apportionment include descriptive statistical methods (coefficient of variation (CV) and correlation analysis), geographic information system (GIS) mapping, enrichment factor (EF) analysis, and positive matrix factorization (PMF). The CV, correlation, and EF analysis approaches determine roughly which heavy metals may be homologous, and which may be from anthropogenic or natural sources [17,18,19]. GIS mapping can estimate and quantify the spatial distributions of heavy metals [2,20,21,22,23], and provide additional information to distinguish point from non-point sources and identify potential pollutant sources [2,24,25]. However, these methods provide approximations only, and cannot definitively quantify the relative source contributions. PMF is therefore more widely used, because it estimates source component spectra based on a weighted-least-squares method and decomposes sample data into two matrices, namely the ‘factor-contribution’ and ‘factor-component’ matrices or spectra [26]. In previous studies, CV and correlation analysis, GIS mapping, EF analysis, and PMF have commonly been used alone, whereas the present study represents an attempt to integrate and apply all methods to a single site, thereby enabling cross-verification and the generation of more accurate and reliable results.
Unlike other studies, a research story in our study from pollution distribution to pollution assessment to sources apportionment was clearly and integrally described at sites in the Lihe River Watershed of the Taihu region, East China. Additionally, from the specific implementation details, a novel systematic risk assessment and a combined method of source apportionment for heavy metals were conducted. The aims were to (1) investigate the concentrations and spatial distribution patterns of Cd, Cr, Cu, Ni, Pb, and Zn in agricultural soils of the Lihe River Watershed; (2) assess the ecological risks posed by these metals based on combined PERI and RAC methods; and (3) identify the natural and anthropogenic sources of the metals and also calculate their contribution rate by combining CV and correlations analysis, GIS mapping, EF analysis, and PMF. The results provide insights into the management of heavy metal pollution in the Lihe River Watershed, and this study serves as a reference for other regions both in China and worldwide.

2. Materials and Methods

2.1. Study Area

The study area, in the Lihe River Watershed, is located to the west of Taihu Lake, Yixing City, Jiangsu Province, China. It includes two towns, Hufu and Dingshu, and is bounded by the coordinates 31°09′00″–31°20′31″ N and 119°42′00″–119°56′20″ E. Lihe River watershed has a surface area of ~260 km2, and agricultural land makes up 41% of the area [27]. There is a high density of industrial activity in the area, including a ceramics factory, a refractory materials plant, and a chemical plant, among others.
This study area was chosen because a previous study involving potential ecological risk assessment of heavy metal pollution in estuaries of the 24 main rivers flowing in and out Taihu Lake revealed that the Lihe River estuary was the most seriously contaminated [14]. However, there has been scant research into risk assessment and the source apportionment of heavy metals in this high-risk area.

2.2. Sample Collection and Preparation

Soils in crop fields were sampled at 0–10 cm depth from 32 randomly selected sites throughout the study area in May 2016 (Figure 1). The crop grown in the soil was wheat. Each sample of 0.5–1.0 kg was divided into 5–9 subsamples. The samples were air-dried at room temperature, ground, and passed through a 2 mm nylon sieve to remove stones and plant roots. The fine soil powders were then stored in polythene zip bags [28].

2.3. Chemical Analysis

2.3.1. Soil Properties

Measurements of the particle size distribution were conducted using a MALVERN 2000 Laser particle size analyzer. The samples first received ultrasonic treatment in a 20% (NaPO3)6 solution, and then 3 cycles were measured with the analyzer to obtain 1/4 interval distribution of granularity data [3,4,29,30]; pH values of the soils were determined using a 10 g sample mixed with 25 mL of distilled water, which was stirred and left at room temperature (25 °C) for 0.5 h and then determined by a pH electrode [31]; loss-on-ignition method was used to analyze the soil organic matter (OM). Soil samples with constant weight were burned at 550 °C in a muffle in this process [32]; cation exchange capacity (CEC) was measured by the extraction method, using ammonium acetate (NH4OAc) and potassium chloride (KCl) [32]. A conductivity meter (Mettler Toledo) was used to measure the soil electrical conductivity (EC). Firstly, a suspension of 1:5 soil–water was configured and centrifuged for 30 min, and the clear supernatant extract was then used to test conductivity [33].

2.3.2. Total Heavy Metal Concentrations

Total heavy metal concentrations were determined by HCl/HNO3/HF/HClO4 extraction. Approximately 100 mg of each sample was digested with 3 mL 37% HCl, 1 mL 65% HNO3, 6 mL 65% HF, and 0.5 mL 65% HClO4, using a two-stage digestion program involving heating for 10 min to 200 °C, and 15 min at 200 °C [34]. After cooling, the digestion solutions were evaporated to near dryness and the residue was dissolved in 1.0 mL 65% HNO3, after which deionized water (20 mL) was added. These solutions were stored in 25 mL high-density polyethylene vials at 4 °C until instrumental analysis [34]. Concentrations of Cd, Zn, Pb, Cu, Cr, and Ni in each solution were determined by inductively coupled plasma–mass spectrometry (ICP–MS; Perkin Elmer SCIEX, Elan 9000, Perkin Elmer Inc., America). Instrumental conditions are listed in Table S1 of the Supplementary Materials.

2.3.3. Acid-Soluble Heavy Metal Fractions

Acid-soluble heavy metal fractions were determined by 0.1 mol L–1 HCl extraction [34]. Approximately 5 g of soil was placed in a pre-weighted 15 mL plastic centrifuge tube and the exact mass was recorded. HCl (5 mL, 0.1 mol L–1) was added and the tube was shaken for 4 h on an oscillating table. After shaking, the mixture was allowed to stand for 24 h before concentrations of Cd, Zn, Pb, Cu, Cr, and Ni were determined by ICP–MS.

2.3.4. Quality Assurance and Control

Quality assurance (QA) and quality control (QC) involved analysis of replicate samples, method blanks, and standard reference materials, with duplicate analyses. Three replicate samples were analyzed, with relative standard deviations (RSD) for metal concentrations being no greater than 5% in all cases. All analyses were duplicated, and results are expressed as mean concentrations. Standard reference material GBW07405 (National Research Center for Standards in China) was used to verify analytical accuracy. All reagents were analytical grade or better. The lab glassware (bottles, tubes, etc.) was pre-cleaned by soaking in 10% HNO3 for at least two days, followed by soaking in and rinsing with de-ionized water prior to use.

2.4. Assessment of Soil Contamination

2.4.1. Potential Ecological Risk Index

The PERI, introduced by Hakanson (1980), is used to assess the potential ecological risk of heavy metals in soil. This method not only assesses the pollution status of the soil, but also combines ecological and environment effects with toxicology, providing a better indication of potential risks denoted by index values. The potential ecological risk factor of a given metal (Ei) is defined as:
E i = T i × ( C i / C 0 )
Equation (2) was used to calculate the total risk index (RI) of a sampling site, as follows:
R I = i = 1 n E i
where Ci is the concentration of metal i in soil, C0 is the background concentration, Ti is the biological toxicity factor of an individual element (Ti values are as follows: Cu = Pb = Ni = 5, Zn = 1, Cr = 2, and Cd = 30 [8,35,36], Ei is the potential ecological risk factor of metal i, and RI is the total potential ecological risk index of metals i–n. The PERI of heavy metals was categorized according to five levels, as shown in Table 1.

2.4.2. Risk Assessment Code

Assessment of soil heavy metal pollution was mainly based on the total concentration currently, and the speciation of heavy metals was often neglected. In fact, the acid-soluble fraction was the real reason for the risk. Therefore, methods of ecological risk assessment based on acid-soluble heavy metals, such as risk assessment code, have provided a new idea for ecological risk assessment of heavy metal pollution. Perin et al. put forward this method in 1985 [37]. The RAC was used to assess the risk and mobility of non-stable fractions of the metals [8], expressed as the percentage acid-soluble fraction (Table 2).

2.5. Source Identification and Apportionment

2.5.1. Source Identification

GIS Mapping

GIS is a spatial information analysis technology that has developed rapidly in recent years. It has a powerful spatial analysis function, and is widely used in source identification. The greatest advantage of this technology is that it can convert huge amounts of data information into a more intuitive spatial distribution map in the process of source identification, and utilize the spatial distribution maps to qualitatively identify the pollution sources. In general, GIS mapping can distinguish point source pollution from non-point source pollution [38]. In spatial analysis using ArcMap 10.0 software, the Kriging method was used for spatial interpolation.

Coefficient of Variation (CV)

When comparing the degrees of dispersion of two groups of data, it is not appropriate to use the standard deviation (SD) if the dimensions of the two groups are different. The CV overcomes this problem because it involves the ratio of the standard deviation of raw data to the mean value, as follows:
CV = SD Mean
Coefficient of variation can be divided into three grades: low variation (CV < 15%), medium variation (15% < CV < 36%), and high variation (CV > 36%) [39].
Coefficient of variation of heavy metal concentration can reflect the interference degree of human activities to environmental media. Generally speaking, environmental media that is strongly affected by human activities may produce a higher coefficient of variation of heavy metal concentration [39]. There are similarities between source identification by calculating coefficient of variation and by using GIS mapping. A high coefficient of variation and an uneven spatial distribution of heavy metals concentration can indicate a high proportion of human sources.

Correlation Analysis of Heavy Metals

Correlation analysis is a relatively simple and convenient method for source identification of heavy metal pollution, and it is one of the most widely used methods currently. Through the correlation analysis of heavy metals, some effective information about the pollution sources has been obtained. It can provide some important information for the further analysis of pollution sources, and also can verify the results of source apportionment obtained by other methods. In general, when the concentrations of two heavy metals are significantly positively correlated, the two heavy metals may be homologous, and when they are negatively correlated, they may have different sources [40]. In this study, Pearson correlation analysis was used in SPSS software. Remarkably, most multivariate statistical methods require variables to conform to the normal distribution. Thus, the normality of the distribution of each variable was checked by analyzing kurtosis and skewness statistical tests before Pearson correlation analysis [41,42].

Enrichment Factors (EF)

Enrichment factor (EF) is a useful tool for determining the degree of anthropogenic heavy metal pollution [11,23]. It was used to describe metal enrichment in agricultural soil, because it is necessary to reduce the effect of grain size on contamination through the normalization of concentration data with respect to a conservative element [9,43]. Iron is inert in migration processes and is derived mainly from natural lithogenic sources, so it is generally recognized as being an appropriate normalization element [44]. The EF values for elements in soil samples were thus calculated as follows [44]:
E F = ( C i / C F e ) s a m p l e ( C i / C F e ) c r u s t
where Ci is the element concentration in the sample or continental crust, and CFe is the concentration of the reference element (Fe). Background values for soils in Jiangsu Province were taken as continental crust values: Fe 3.02%, Cd 0.126 mg kg−1, Cr 77.8 mg kg–1, Cu 22.3 mg kg–1, Ni 26.7 mg kg–1, Pb 26.2 mg kg–1, and Zn 62.6 mg kg–1 [45]. An EF value of 0.5–1.5 indicates that a metal is derived mainly from crustal materials or natural weathering processes, while EF > 1.5 indicates that a significant portion originated from non-crustal or anthropogenic processes [46].

2.5.2. Source Apportionment by Positive Matrix Factorization (PMF)

The above methods only qualitatively identify the major sources of heavy metals in soil. However, it is impossible to quantitatively calculate the contribution rate of each type of pollution source using them. PMF analysis can achieve the goal of source apportionment.
PMF analysis employed the US EPA PMF 5.0 model for soil pollutant source apportionment, as developed by Paatero and Tapper [47]. The PMF is a typical receptor model and can be used without source composition as an input. The notable features of this approach are non-negativity constraints in obtaining physically realistic meanings, and the use of uncertainty to weight each data point individually [48]. Furthermore, PMF can supplement values that are missing or below detection limits to ensure the reliability of each data point [26].
The PMF model defines an n × m original data matrix X, where n represents the number of samples and m the number of chemical species, which can be factorized into two sub-matrices, namely G (n × p) and F (p × m), with an unexplained part E (n × m) [26,49], as follows:
X = G + F
x i j = k = 1 p g i k f k j + e i j
where xij is the concentration of the jth chemical species measured in the ith sample; gik is the contribution of source k to the ith sample; fkj is the concentration of the jth chemical species in source k; and eij is the residual for each sample and species. The residual error matrix, eij, is obtained by minimizing the function Q, where
Q = i = 1 n j = 1 m ( e i j u i j ) 2
and uij is the uncertainty of the jth chemical species for sample i. The model can be optimized by including the uncertainty, u, for each sample, calculated using the measurement uncertainty (MU) and detection limit (MDL). When the metal concentration is ≤MDL, u is calculated as
u = 5 6 × MDL
and when the concentration is >MDL, we have
u = ( MU × concentration ) 2 + ( MDL ) 2

3. Results and Discussion

3.1. Concentrations of Heavy Metals in Soils

3.1.1. Descriptive Statistics of Soil Properties

Descriptive statistics of soil properties are shown in Table 3. The pH of the soil ranged from 4.89 to 7.45, with an average of 5.81. The proportion of samples in which pH was less than 7 among the 32 soils was 93.75%. Soil cation exchange capacity (CEC) is an important factor reflecting soil buffer capacity, and also an important basis for evaluating soil nutrient holding capacity, improving soil equality, and rational fertilization. The average soil cation exchange capacity in the study area was 18.41 cmol·kg−1, which belongs to a medium-upper level. Soil electrical conductivity (EC) is an index for determining water-soluble salts in soil, and it is a factor to determine whether salt ions in soil restrict crop growth or not. Different plants have different optimum intervals for soil conductivity. The average conductivity of soil in the study area was 1.59 ms/cm. Soil organic matter content (OM) is an index reflecting soil fertility. According to soil fertility classification recommended by the second national soil census, if the organic matter of a soil is less than 1%, it is organic matter deficient soil. There was no sample in which the organic matter was less than 1%, which indicated that the soil fertility in the study area was relatively good. Soil particle size was related to soil adsorption capacity of heavy metals. Generally, the smaller the particle size of soil is, the stronger the adsorption capacity and the higher the heavy metal concentration would be. The average soil particle size in the study area was 24.37Φ.

3.1.2. Descriptive Statistics of Total Heavy Metal Concentrations

Basic statistics related to the total heavy metal concentrations in soils from the Lihe River Watershed and other selected areas are summarized in Table 4, together with background values of soil in Jiangsu Province and standard values (Soil Environmental Quality Standard Two, SEQ-II, GB 15618-1995, grade II for agricultural land). SEQ-II values are threshold values for the protection of human health and agricultural production in China. A wide range of total heavy metal concentrations was observed in soils across the study area (Table 4), and mean Cd, Cu, Pb, and Zn concentrations were all higher than their respective background values (the Cd concentration was four times the background value). Mean Cr and Ni concentrations were lower than their background values. The Cd, Cr, Cu, Ni, Pb, and Zn concentrations exceeded background values in 100%, 3%, 31%, 25%, 94%, and 94% of the 32 soil samples, respectively. Comparison of mean metal concentrations with SEQ-II values indicates that Cd pollution is the most serious, with a concentration almost twice the SEQ-II value, indicating that anthropogenic sources have a direct effect on soil Cd concentrations [50]. The CV values of metal concentrations followed the order Cd (0.72) > Pb (0.52) > Cu (0.42) > Zn (0.39) > Ni (0.20) > Cr (0.17).
In comparisons with other regions (Table 4) the mean concentration of Cd in the Lihe River Watershed (0.57 mg kg–1) was higher than that in the Yangtze and Pearl river systems (0.19 and 0.32 mg kg–1, respectively [9,51], while concentrations of the other metals were both higher and lower. This indicates that the problem of high Cd levels in soil of the Lihe River Watershed requires attention.
Spatial distributions of the metals are shown in Figure 2, revealing that the distributions of Cd, Pb, Cu, and Zn are largely similar, being higher in central-eastern and southern parts of the study area. Cr and Ni are distributed relatively uniformly over the whole study area, although they occur at slightly higher levels in the west and lower levels in the east.

3.1.3. Acid-Soluble Heavy Metal Concentration

The ratio of acid-soluble concentration to total concentration of six heavy metals can reflect the mobility of heavy metals in soil. Descriptive statistics of acid-soluble heavy metal concentrations are shown in Table 5. The order of acid-soluble concentrations of the six heavy metals was Pb > Cu > Cd > Zn > Ni > Cr, which indicated that the mobility of Pb and Cu was relatively high, while that of Ni and Cr was low. The result was consistent with the study by Dario et al. [51]. This suggested that for Cr and Ni, the input ratio of anthropogenic sources was relatively low, while it was high for Pb and Cu. It is noteworthy that the mobile heavy metals in soil are easily absorbed by plants, and eventually enter the human body through the food chain, which poses a great threat to the human body. Therefore, Pb, Cu, Cd, and Zn in the soil of the study area are the key heavy metals that posed a threat to the human body and need to be paid attention, while the threat of Cr and Ni to the human body is relatively small.

3.1.4. Correlation between Heavy Metal Concentration and Soil Properties

The correlation between heavy mental concentration and soil properties is shown in Table 6. It was suggested that the total concentration of Cr was significantly correlated with CEC and pH. The available concentrations of Cr and Ni were positively correlated with CEC at 0.01 level. However, for other heavy metals, there was no correlation between the concentration and soil properties. Thus, Cr and Ni were more affected by soil properties. Whether this result indicates that Cr and Ni derived mainly from parent material is a question worthy of further verification.

3.2. Risk Assessment of Heavy Metals in Soil

3.2.1. Potential Ecological Risk Indices

Potential ecological risks of the metals were calculated using Equations (1) and (2), with results presented in Table 7. Mean potential ecological risks are ranked as Cd > Pb > Cu > Ni > Zn > Cr. The Ei values (Equation (1)) for Cr, Cu, Ni, Pb, and Zn are all <40, indicating low potential ecological risk, whereas Cd (Ei = 131.9) poses a considerable risk. The mean RI value (Equation (2)) for the study area indicates a moderate risk, to which the Ei value for Cd contributes 86.8% of the potential ecological risk.

3.2.2. Risk Assessment Code

The acid-soluble fraction indicates the degree of mobility of heavy metals in soil, further reflecting ecological risk, and follows the order Pb > Cu > Cd > Zn > Ni > Cr (Table 8). According to the RAC classification (Table 2), Pb, Cu, and Cd exhibit the highest mobility and therefore the greatest potential bio-available risk, with 43–68% solubility. The lower acid solubility of Zn and Ni indicate medium risk, with low risk for Cr. Overall, the environmental risk of the bio-available fraction, based on RAC values, decreases in the order Pb > Cu > Cd > Zn > Ni > Cr, indicating that Pb, Cu, and Cd pose the greatest ecological risk in the Lihe River Watershed.

3.3. Identification and Apportionment of Pollution Sources

3.3.1. Identification of Pollution Sources

GIS Mapping of Heavy Metals

GIS mapping can assist in distinguishing point from non-point sources, and in identifying potential pollutant sources [2,20,21,24]. The higher Cd, Pb, Cu, and Zn concentrations tend to be clustered in circular patterns around specific points (Figure 2), indicating point-source pollution. In contrast, Cr and Ni concentrations are relatively uniform, with no obvious low or high values in the study area, indicating geological background or non-point sources.

Coefficient of Variation of Metal Concentrations

Generally, the greater the CV, the higher the contribution of anthropogenic sources to metals in soils [40]. Therefore, the calculated CV values (Table 4) indicate that Cd, Pb, Cu, and Zn (CV = 0.39–0.72) are derived mainly from anthropogenic sources, while Cr and Ni (CV = 0.17–0.20) are mainly crustal in origin.

Correlations between Metals

Inter-elemental relationships can provide information regarding pollution sources [52], as the heavy metals of most interest commonly have a shared origin [53]. Before correlation analysis, the normality of the distribution of heavy metal concentration was checked [42]. Through analysis in SPSS, our statistical data were conformed to normal distribution. Pearson correlation coefficients for metal concentrations (Table 9) indicate highly positive correlations (p ≤ 0.01) between Cd and Pb (r = 0.84), Cr and Ni (r = 0.68), Cu and Zn (r = 0.76), and Cu and Pb (r = 0.62), indicating that these pairs of metals may have the same origin.

Enrichment Factors

The statistics of the EFs for all the metals analyzed were evaluated to further clarify the potential sources of metals in soils in the study area. The EF values for Cd are highest, with a mean of 15.56, and values for Pb and Zn are moderate, with means of 2.63 and 2.58, respectively. Mean EF values (and CV values) follow the order Cd > Pb > Zn > Cu > Ni > Cr. Cd, Pb, Zn, and Cu are thus considered to be enriched (EF = 1.75–15.56), while Cr and Ni show relatively low enrichment (EF = 1.19–1.53). Considered together, the mean EF and CV values (Table 10) indicate that Cd, Pb, Cu, and Zn have anthropogenic origins, whereas Cr and Ni are derived from the crust.

3.3.2. Source Apportionment by PMF

The PMF model was run using between three and six factors, and each run was initialized with different starting points. A random-seed mode with 20 random starting points was selected and 3–6 factors were examined. PMF analysis identified four appropriate factors among the soils (Table 11). Tentative identification of the factor profiles was based on elemental markers in the different sources.
Source composition profiles based on a four-factor solution are provided in Table 11. Factor 1 is dominated by Cd and accounts for 52.11% of the Cd source. The EF value for Cd is the highest of the six metals, averaging 15.56, which indicates Cd is significantly enriched in the soil. In addition, the results of the PERI and RAC assessments indicate that Cd poses a high ecological risk. The widely distributed ceramics industries in the study area use large amounts of Cadmium Yellow pigment, which may be an important source of the Cd pollution. It is clear, therefore, that the first factor represents a source related to the ceramics industry.
The second factor is dominated by Cu, Pb, and Zn. Concentrations of these metals are higher than background values, with mean EF values of 1.75, 2.63, and 2.58 (Table 8), respectively, which indicate moderate enrichment in soil. The spatial distributions of Cu, Pb, and Zn in the study area are uneven (Figure 2). Previous studies have indicated that industrial and vehicular sources are mainly responsible for Cu, Pb, and Zn inputs to soils and the environment [54], with Cu being mainly from machinery manufacturing plants, Zn (a hardness additive) from tire dust, and Pb from coal combustion and automobile exhaust emissions [53,55]. Since the phasing out of leaded petrol, the input of Pb to the environment has reduced, but vehicle brakes and tires still release Pb, meaning that vehicular emissions remain an important source of pollution [54]. The development of coal-fired power stations and heavy vehicle traffic have accompany industrialization in the study area, so the second factor can be attributed to industrial activity and vehicular emissions.
In the third factor, all of the elements have some loading, except for Ca. The large-scale use of agricultural chemicals (e.g., manure, fertilizers, and pesticides) may cause the enrichment of many elements in soil, including Cd, Cr, Cu, Ni, Pb, Zn, Co, Mn, K, and Mg [26,53,55,56]. Some of these elements, such as Cu, Ni, Pb, Zn, Co, Mn, K, and Mg, are important nutrients in crop production, and are provided primarily by commercial fertilizers to improve the soil fertility. Some, such as Cd, Cr, and Pb, are not essential for plant growth and are not intentionally introduced. PMF analysis thus indicates that the third factor represents sources based on agronomic practices.
For the fourth factor, Ca, Co, Cr, Fe, K, Mg, Mn, and Ni received higher weightings than other metals. These elements, in particular Ca, K, Mg, and Mn, are associated with crustal material. Minerals in the topsoil of the study area comprise mainly quartz, potassium feldspar, and anorthoclase, which contain Al2O3, CaO, MgO, SiO2, and K2O [26], and these oxides are released during weathering. Concentrations of Cr and Ni are slightly lower than background values, and their uniform distribution in the study area and low CVs indicate that they have low spatial variability. The mean EF values of Cr and Ni are ≤1.5, indicating minimal enrichment. In addition, PERI and RAC assessments indicate that Cr and Ni pose low or zero potential ecological risk. This group of elements is thus likely to originate from natural sources, such as mineral weathering and atmospheric precipitation. This is consistent with the findings of previous studies [7,8,26,53,56].
The present results concerning heavy metal source identification are largely consistent with those of previous studies. Xue et al. [26] combined PMF with a geo-statistical approach in Changxing County, Zhejiang Province, and found that Cd, Cu, Pb, and Zn were mainly from anthropogenic sources, while Cr and Ni were mainly from natural sources. Dai et al. [56] used multivariate analysis and geostatistics to identify the sources of heavy metals in Laiwu City, Shandong Province, Eastern China, and found that Cd, Cu, Pb, and Zn, were dominated by industrial, agricultural, and vehicular sources, while Cr and Ni came from natural sources and were controlled by parent materials. Peris et al. [57] performed a multivariate statistical analysis of agricultural soils in the European Mediterranean region, and found that Ni was mainly associated with parental rocks, while Cd, Cu, Pb, and Zn were connected with anthropogenic activity. Using a combination of GIS mapping and multivariate statistical analysis, Cai et al. [53] demonstrated that Ni and Cr originate from natural sources, whereas Cd and Pb are mainly derived from anthropogenic activity. Li et al. [7] used correlation and multivariate statistical analysis to show that Cu, Cd, Pb, and Zn concentrations in soils were mainly anthropogenic in origin, whereas Cr and Ni concentrations were mainly influenced by natural factors, such as parent materials.
Source contribution rates were also determined by PMF (Figure 3). Factor 1 (representing industrial activity related to the ceramics industry) and factor 2 (industrial and vehicular emissions) were combined. The total contributions of industrial and vehicular emissions to Cu, Pb, and Zn levels in the agricultural soils are 71.4%, 75.0%, and 62.2%, respectively, while agricultural sources contribute 6.8%, 15.7%, and 6.9%, respectively, and parent materials contribute 21.8%, 9.2%, and 30.9%, respectively. For Cr and Ni, the contribution rates of parent material are highest, followed by agricultural activity and industrial and vehicular emissions. For Cd, industrial and vehicular emissions, and agricultural sources account for 73.7% and 26.3%, respectively. There is no natural source of Cd in the study area, and because Cd pollution is overwhelmingly identified as the most significant ecological risk factor in the study area, it should have the highest priority in controls imposed on heavy metal pollution by industrial sources.

4. Conclusions

This study analyzed the spatial distribution, ecological risks, and sources of the heavy metals Cd, Cr, Cu, Ni, Pb, and Zn in agricultural soils in the Lihe River Watershed. Average concentrations of Cd, Cu, Pb, and Zn were all higher than background values, with the mean Cd concentration being four times higher and almost twice the SEQ-II value. Average Cr and Ni concentrations were lower than background values. Cd, Cu, Pb, and Zn exhibit similar spatial distributions, with concentrations clearly higher in the central-eastern and southern parts of the study area. In contrast, Cr and Ni have relatively uniform distributions across whole study area, although slightly higher in the west and slightly lower in the east. Based on the results of PERI and RAC assessments, it was concluded that Cd poses a considerable ecological risk, and it was identified as the most significant factor affecting the ecological environment of the study area.
Source identification based on the combined methods of GIS mapping and CV, correlation, and EF analyses indicates that Cd, Pb, Cu, and Zn are derived mainly from anthropogenic sources, including agricultural, industrial, and vehicular emissions, while Cr and Ni originate mainly from natural sources. Source apportionment by PMF analysis indicates that for Cd, Cr, Cu, Ni, Pb, and Zn in the agricultural soils, industrial and vehicular emissions contribute 73.7%, 21.3%, 71.4%, 20.3%, 75.0%, and 62.2%, respectively; agricultural sources contribute 26.3%, 36.3%, 6.8%, 38.9%, 15.7%, and 6.9%, respectively; and parent materials contribute 0%, 42.4%, 21.8%, 40.8%, 9.2%, and 30.9%, respectively. Important strategies should be implemented to reduce the use of chemical fertilizers and pesticides and, even more importantly, control industrial point-source pollution to reduce ecological risks associated with Cd pollution.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/16/12/2094/s1, Table S1: Operating parameters for ICP-MS (ELAN 9000, PerkinElmer SCIEX) for the determination of elemental concentrations, Figure S1: The distribution of heavy metal concentration in the study area.

Author Contributions

S.Z., S.W. and G.W. conceived and designed the experiments and supervised the quality of the article; C.W. and L.C. performed the experiments; Z.X. and Z.C. provided the payment for publishing the article and put forward suggestions for revision of the article analyzed the data; and L.C. wrote the paper.

Funding

This research was funded by the Comprehensive Geological Survey in the Coastal Area of Northeastern Hainan Island (DD20190308), the National Natural Science Foundation of China (41771243), the Comprehensive Geological Survey of Beibu Gulf and Other Key Coastal Zones (DD20160146), the National Key Research and Development Plan (2017YFD0800305), the Special Fund for Research in the Public Interest of the Ministry of Land and Resources (201511001-03) and Guangdong Natural Science Foundation Project (2018A0303130063).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, T.; Liu, X.M.; Li, X.; Zhao, K.; Zhang, J.B.; Xu, J.M.; Shi, J.C.; Dahlgren, R.A. Heavy metal sources identification and sampling uncertainty analysis in a field-scale vegetable soil of Hangzhou, China. Environ. Pollut. 2009, 157, 1003–1010. [Google Scholar] [CrossRef] [PubMed]
  2. Shao, S.; Hu, B.; Fu, Z.; Wang, J.; Lou, G.; Zhou, Y.; Jin, B.; Li, Y.; Shi, Z. Source Identification and Apportionment of Trace Elements in Soils in the Yangtze River Delta, China. Int. J. Environ. Res. Public Health 2018, 15, 1240. [Google Scholar] [CrossRef] [PubMed]
  3. Shamuyarira, K.K.; Gumbo, J.R. Assessment of heavy metals in municipal sewage sludge: A case study of Limpopo province, South Africa. Int. J. Environ. Res. Public Health 2014, 11, 2569–2579. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, L.; Zhou, S.; Shi, Y.; Wang, C.; Li, B.; Li, Y.; Wu, S. Heavy metals in food crops, soil, and water in the Lihe River Watershed of the Taihu Region and their potential health risks when ingested. Sci. Total Environ. 2017, 615, 141. [Google Scholar] [CrossRef] [PubMed]
  5. Laidlaw, M.A.S.; Filippelli, G.M.; Brown, S.; Paz-Ferreiro, J.; Reichman, S.M.; Netherway, P.; Truskewycz, A.; Ball, A.S.; Mielke, H.W. Case studies and evidence-based approaches to addressing urban soil lead contamination. Appl. Geochem. 2017, 83, 14–20. [Google Scholar] [CrossRef]
  6. Lu, A.; Wang, J.; Qin, X.; Wang, K.; Han, P.; Zhang, S. Multivariate and geostatistical analyses of the spatial distribution and origin of heavy metals in the agricultural soils in Shunyi, Beijing, China. Sci. Total Environ. 2012, 425, 66. [Google Scholar] [CrossRef] [PubMed]
  7. Li, C.; Wang, F.; Cao, W.; Pan, J.; Lv, J.; Wu, Q. Source analysis, spatial distribution and pollution assessment of heavy metals in sewage irrigation area farmland soils of Longkou city. Environ. Sci. 2017, 38, 1018–1027. [Google Scholar]
  8. Ke, X.; Gui, S.; Huang, H.; Zhang, H.; Wang, C.; Guo, W. Ecological risk assessment and source identification for heavy metals in surface sediment from the Liaohe River protected area, China. Chemosphere 2017, 175, 473. [Google Scholar] [CrossRef]
  9. Wang, H.; Wang, J.; Liu, R.; Yu, W.; Shen, Z. Spatial variation, environmental risk and biological hazard assessment of heavy metals in surface sediments of the Yangtze River estuary. Mar. Pollut. Bull. 2015, 93, 250–258. [Google Scholar] [CrossRef]
  10. Liu, E.; Wang, X.; Liu, H.; Liang, M.; Zhu, Y.; Li, Z. Chemical speciation, pollution and ecological risk of toxic metals in readily washed off road dust in a megacity (Nanjing) China. Ecotoxicol. Environ. Saf. 2019, 173, 381–392. [Google Scholar] [CrossRef]
  11. Maanan, M.; Saddik, M.; Maanan, M.; Chaibi, M.; Assobhei, O.; Zourarah, B. Environmental and ecological risk assessment of heavy metals in sediments of Nador lagoon, Morocco. Ecol. Indic. 2015, 48, 616–626. [Google Scholar] [CrossRef]
  12. Perin, G.; Fabris, R.; Manente, S.; Wagener, A.R.; Hamacher, C.; Scotto, S. A five-year study on the heavy-metal pollution of Guanabara Bay sediments (Rio de Janeiro, Brazil) and evaluation of the metal bioavailability by means of geochemical speciation. Water Res. 1997, 31, 3017–3028. [Google Scholar] [CrossRef]
  13. Duodu, G.O.; Goonetilleke, A.; Ayoko, G.A. Potential bioavailability assessment, source apportionment and ecological risk of heavy metals in the sediment of Brisbane River estuary, Australia. Mar. Pollut. Bull. 2017, 117, 523–531. [Google Scholar] [CrossRef] [PubMed]
  14. Cao, C.; Zhang, Q.; Ma, Z.B.; Wang, X.M.; Chen, H.; Wang, J.J. Fractionation and mobility risks of heavy metals and metalloids in wastewater-irrigated agricultural soils from greenhouses and fields in Gansu, China. Geoderma 2018, 328, 1–9. [Google Scholar] [CrossRef]
  15. Kachenko, A.G.; Singh, B. Heavy Metals Contamination in Vegetables Grown in Urban and Metal Smelter Contaminated Sites in Australia. Water Air Soil Pollut. 2006, 169, 101–123. [Google Scholar] [CrossRef]
  16. Franco-Uría, A.; López-Mateo, C.; Roca, E.; Fernández-Marcos, M.L. Source identification of heavy metals in pastureland by multivariate analysis in NW Spain. J. Hazard. Mater. 2009, 165, 1008–1015. [Google Scholar] [CrossRef] [PubMed]
  17. Liang, Y.; Jiang, X.; Li, G.; Ren, H.; Yang, F.; Cao, Q. Content analysis and risk evaluation of heavy metal elements in soil and ginseng root from nonforested land base with cultivated ginseng in Jilin Province. J. Plant Res. Environ. 2015, 24, 68–76. [Google Scholar]
  18. Brady, J.P.; Ayoko, G.A.; Martens, W.N.; Goonetilleke, A. Enrichment, distribution and sources of heavy metals in the sediments of Deception Bay, Queensland, Australia. Mar. Pollut. Bull. 2014, 81, 248–255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Chen, J.B.; Gaillardet, J.; Bouchez, J.; Louvat, P.; Wang, Y.N. Anthropophile elements in river sediments: Overview from the Seine River, France. Geochem.Geophys. Geosyst. 2015, 15, 4526–4546. [Google Scholar] [CrossRef]
  20. Li, X.; Lee, S.L.; Wong, S.C.; Shi, W.; Thornton, I. The study of metal contamination in urban soils of Hong Kong using a GIS-based approach. Environ. Pollut. 2004, 129, 113. [Google Scholar] [CrossRef]
  21. Zhang, X.Y.; Lin, F.F.; Wong, M.T.F.; Feng, X.L.; Wang, K. Identification of soil heavy metal sources from anthropogenic activities and pollution assessment of Fuyang County, China. Environ. Monit. Assess. 2009, 154, 439. [Google Scholar] [CrossRef] [PubMed]
  22. Maas, S.; Scheifler, R.; Benslama, M.; Crini, N.; Lucot, E.; Brahmia, Z.; Benyacoub, S.; Giraudoux, P. Spatial distribution of heavy metal concentrations in urban, suburban and agricultural soils in a Mediterranean city of Algeria. Environ. Pollut. 2010, 158, 2294–2301. [Google Scholar] [CrossRef] [PubMed]
  23. Tume, P.; Roca, N.; Rubio, R.; King, R.; Bech, J. An assessment of the potentially hazardous element contamination in urban soils of Arica, Chile. J. Geochem. Explor. 2018, 184, 345–357. [Google Scholar] [CrossRef]
  24. Xie, Y.; Chen, T.B.; Lei, M.; Yang, J.; Guo, Q.J.; Song, B.; Zhou, X.Y. Spatial distribution of soil heavy metal pollution estimated by different interpolation methods: Accuracy and uncertainty analysis. Chemosphere 2011, 82, 468. [Google Scholar] [CrossRef] [PubMed]
  25. Tóth, G.; Hermann, T.; Szatmári, G.; Pásztor, L. Maps of heavy metals in the soils of the European Union and proposed priority areas for detailed assessment. Sci. Total Environ. 2016, 565, 1054–1062. [Google Scholar] [CrossRef] [PubMed]
  26. Xue, J.L.; Zhi, Y.Y.; Yang, L.P.; Shi, J.C.; Zeng, L.Z.; Wu, L.S. Positive matrix factorization as source apportionment of soil lead and cadmium around a battery plant (Changxing County, China). Environ. Sci. Pollut. Res. 2014, 21, 7698–7707. [Google Scholar] [CrossRef] [PubMed]
  27. Li, H.P.; Huang, W.Y.; Yang, G.S.; Liu, X.M. Non-point pollutant concentrations for different land uses in Lihe River watershed of Taihu Region China. Environ. Sci. 2006, 26, 243–247. [Google Scholar]
  28. Li, Q.; Chen, Y.; Fu, H.; Cui, Z.; Shi, L.; Wang, L.; Liu, Z. Health risk of heavy metals in food crops grown on reclaimed tidal flat soil in the Pearl River Estuary, China. J. Hazard. Mater. 2012, 227, 148–154. [Google Scholar] [CrossRef]
  29. Gao, W.; Du, Y.; Wang, D.; Gao, S. Distribution Patterns of Heavy Metal in Sediment and Their Influence on the Environment Quality of the Intertidal Flat of Luoyuan, Bay, Fujian Coast. Environ. Sci. 2012, 33, 3097–3103. [Google Scholar]
  30. Cheng, P.; Gao, S.; Bokuniewicz, H. Net sediment transport patterns over the Bohai Strait based on grain size trend analysis. Estuar. Coast. Shelf Sci. 2004, 60, 203–212. [Google Scholar] [CrossRef]
  31. Huang, M.; Zhou, S.; Sun, B.; Zhao, Q. Heavy metals in wheat grain: Assessment of potential health risk for inhabitants in Kunshan, China. Sci. Total Environ. 2008, 405, 54–61. [Google Scholar] [CrossRef] [PubMed]
  32. Kwak, J.I.; Nam, S.-H.; Kim, S.W.; Bajagain, R.; Jeong, S.-W.; An, Y.-J. Changes in soil properties after remediation influence the performance and survival of soil algae and earthworm. Ecotoxicol. Environ. Saf. 2019, 174, 189–196. [Google Scholar] [CrossRef] [PubMed]
  33. Xie, X.; Pu, L.; Shen, H.; Wang, X.; Zhu, M.; Ge, Y.; Sun, L. Effects of soil reclamation on the oat cultivation in the newly reclaimed coastal land, eastern China. Ecol. Eng. 2019, 129, 115–122. [Google Scholar] [CrossRef]
  34. Hu, X.; Zhang, Y.; Luo, J.; Wang, T.; Lian, H.; Ding, Z. Bioaccessibility and health risk of arsenic, mercury and other metals in urban street dusts from a mega-city, Nanjing, China. Environ. Pollut. 2011, 159, 1215–1221. [Google Scholar] [CrossRef] [PubMed]
  35. Lu, S.G.; Bai, S.Q. Contamination and potential mobility assessment of heavy metals in urban soils of Hangzhou, China: Relationship with different land uses. Environ. Earth Sci. 2010, 60, 1481–1490. [Google Scholar] [CrossRef]
  36. Saravanan, P.; Krishnakumar, S.; Pradhap, D.; Silva, J.D.; Arumugam, K.; Magesh, N.S.; Srinivasalu, S. Elemental concentration based potential ecological risk (PER) status of the surface sediments, Pulicat lagoon, Southeast coast of India. Mar. Pollut. Bull. 2018, 133, 107–116. [Google Scholar] [CrossRef] [PubMed]
  37. Gan, G.; Liu, W.; Qiu, Y.; Li, K.; Hou, H.; Peng, P.; Liao, B. Heavy metal pollution and ecological risk assessment of the paddy soils in a smelting area in Central Hunan. Environ. Chem. 2013, 32, 132–137. [Google Scholar]
  38. Huang, Y.; Li, T.; Wu, C.; He, Z.; Japenga, J.; Deng, M.; Yang, X. An integrated approach to assess heavy metal source apportionment in peri-urban agricultural soils. J. Hazard. Mater. 2015, 299, 540–549. [Google Scholar] [CrossRef]
  39. Lv, J.; Liu, Y.; Zhang, Z.; Dai, J. Factorial kriging and stepwise regression approach to identify environmental factors influencing spatial multi-scale variability of heavy metals in soils. J. Hazard. Mater. 2013, 261, 387–397. [Google Scholar] [CrossRef]
  40. Wang, S.; Cai, L.; Wen, H.; Luo, J.; Wang, Q.; Liu, X. Spatial distribution and source apportionment of heavy metals in soil from a typical county-level city of Guangdong Province, China. Sci. Total Environ. 2019, 655, 92–101. [Google Scholar] [CrossRef]
  41. Naifar, I.; Pereira, F.; Zmemla, R.; Bouaziz, M.; Elleuch, B.; Garcia, D. Spatial distribution and contamination assessment of heavy metals in marine sediments of the southern coast of Sfax, Gabes Gulf, Tunisia. Mar. Pollut. Bull. 2018, 131, 53–62. [Google Scholar] [CrossRef] [PubMed]
  42. Reimann, C.; Filzmoser, P.; Fabian, K.; Hron, K.; Birke, M.; Demetriades, A.; Dinelli, E.; Ladenberger, A. The concept of compositional data analysis in practice:Total major element concentrations in agricultural and grazing land soils of Europe. Sci. Total Environ. 2012, 426, 196–210. [Google Scholar] [CrossRef] [PubMed]
  43. Huang, P.; Li, T.G.; Li, A.C.; Yu, X.K.; Hu, N.J. Distribution, enrichment and sources of heavy metals in surface sediments of the North Yellow Sea. Cont. Shelf Res. 2014, 73, 1–13. [Google Scholar] [CrossRef]
  44. Hu, X. Size Distribution and Source Apportionment of Airborne Metallic Elements in Nanjing, China. Aerosol Air Qual. Res. 2013, 13, 1796–1806. [Google Scholar] [CrossRef] [Green Version]
  45. Hu, W.; Wang, H.; Dong, L.; Huang, B.; Borggaard, O.K.; Bruun Hansen, H.C.; He, Y.; Holm, P.E. Source identification of heavy metals in peri-urban agricultural soils of southeast China: An integrated approach. Environ. Pollut. 2018, 237, 650–661. [Google Scholar] [CrossRef] [PubMed]
  46. Zhang, J.; Liu, C.L. Riverine Composition and Estuarine Geochemistry of Particulate Metals in China—Weathering Features, Anthropogenic Impact and Chemical Fluxes. Estuar. Coast. Shelf Sci. 2002, 54, 1051–1070. [Google Scholar] [CrossRef]
  47. Paatero, P.; Tapper, U. Positive matrix factorization: A non-negative factor model with optimal utilization of error estimates of data values. Environmetrics 1994, 5, 111–126. [Google Scholar] [CrossRef]
  48. Paatero, P. Least squares formulation of robust non-negative factor analysis. Chemom. Intell. Lab. Syst. 1997, 37, 23–35. [Google Scholar] [CrossRef]
  49. Wang, C.; Wu, S.; Zhou, S.L.; Wang, H.; Li, B.; Chen, H.; Yu, Y.; Shi, Y. Polycyclic aromatic hydrocarbons in soils from urban to rural areas in Nanjing: Concentration, source, spatial distribution, and potential human health risk. Sci. Total Environ. 2015, 375, 527–528. [Google Scholar] [CrossRef]
  50. Chai, M.; Shi, F.; Li, R.; Shen, X. Heavy metal contamination and ecological risk in Spartina alterniflora marsh in intertidal sediments of Bohai Bay, China. Mar. Pollut. Bull. 2014, 84, 115. [Google Scholar] [CrossRef]
  51. Di Giuseppe, D.; Vittori Antisari, L.; Ferronato, C.; Bianchini, G. New insights on mobility and bioavailability of heavy metals in soils of the Padanian alluvial plain (Ferrara Province, northern Italy). Geochemistry 2014, 74, 615–623. [Google Scholar] [CrossRef]
  52. Tang, W.; Shan, B.; Zhang, H.; Mao, Z. Heavy metal sources and associated risk in response to agricultural intensification in the estuarine sediments of Chaohu Lake Valley, East China. J. Hazard. Mater. 2010, 176, 945. [Google Scholar] [CrossRef] [PubMed]
  53. Cai, L.; Xu, Z.; Ren, M.; Guo, Q.; Hu, X.; Hu, G.; Wan, H.; Peng, P. Source identification of eight hazardous heavy metals in agricultural soils of Huizhou, Guangdong Province, China. Ecotoxicol. Environ. Saf. 2012, 78, 2–8. [Google Scholar] [CrossRef] [PubMed]
  54. Imperato, M.; Adamo, P.; Naimo, D.; Arienzo, M.; Stanzione, D.; Violante, P. Spatial distribution of heavy metals in urban soils of Naples city (Italy). Environ. Pollut. 2003, 124, 247–256. [Google Scholar] [CrossRef]
  55. Duzgoren-Aydin, N.S.; Weiss, A.L. Use and abuse of Pb-isotope fingerprinting technique and GIS mapping data to assess lead in environmental studies. Environ. Geochem. Health 2008, 30, 577–588. [Google Scholar] [CrossRef] [PubMed]
  56. Dai, B.; Shu, L.J.; Cheng, Z.J.; Lu, Z.Z.; Yang, L.; Jia, Z.R. Assessment of Sources, Spatial Distribution and Ecological Risk of Heavy Metals in Soils in a Typical Industry-based City of Shandong Province, Eastern China. Environ. Sci. 2015, 36, 507. [Google Scholar]
  57. Peris, M.; Recatalá, L.; Micó, C.; Sánchez, R.; Sánchez, J. Increasing the Knowledge of Heavy Metal Contents and Sources in Agricultural Soils of the European Mediterranean Region. Water Air Soil Pollut. 2008, 192, 25–37. [Google Scholar] [CrossRef]
Figure 1. Location of the study area and distribution of sampling points.
Figure 1. Location of the study area and distribution of sampling points.
Ijerph 16 02094 g001
Figure 2. Spatial distributions of Cd, Cr, Cu, Ni, Pb, and Zn concentrations (mg/kg) in the study area. The spatial distributions of elements were determined with the kriging interpolation method.
Figure 2. Spatial distributions of Cd, Cr, Cu, Ni, Pb, and Zn concentrations (mg/kg) in the study area. The spatial distributions of elements were determined with the kriging interpolation method.
Ijerph 16 02094 g002
Figure 3. Relative source contributions to the heavy metals analyzed in the present study. I&TE, AA, and PM represent industrial and vehicular emissions, agricultural activity, and parent materials, respectively.
Figure 3. Relative source contributions to the heavy metals analyzed in the present study. I&TE, AA, and PM represent industrial and vehicular emissions, agricultural activity, and parent materials, respectively.
Ijerph 16 02094 g003
Table 1. Classification of potential ecological risk index (PERI).
Table 1. Classification of potential ecological risk index (PERI).
Assessment CriterionPERI
LowModerateConsiderableHighVery High
Ei<4040–8080–160160–320≥320
RI<150150–300300–600≥600
Table 2. Classification of risk assessment code (RAC).
Table 2. Classification of risk assessment code (RAC).
RiskAcid Soluble Fraction (%)
No risk<1
Low risk1–10
Medium risk11–30
High risk31–50
Very high risk51–75
Table 3. Soil properties including pH, cation exchange capacity (CEC), electrical conductivity (EC), organic matter (OM), and particle size.
Table 3. Soil properties including pH, cation exchange capacity (CEC), electrical conductivity (EC), organic matter (OM), and particle size.
Soil PropertiesRangeMean ValueStandard DeviationCV
pH4.85–7.455.810.690.12
CEC (cmol·kg−1)11.76–22.9918.412.600.14
EC (ms·cm−1)0.67–3.051.590.650.41
OM (%)1.84–5.353.580.940.26
Average particle size (Φ)8.45–52.0624.378.370.34
CV: Coefficient of Variation.
Table 4. Total heavy metal concentrations in soil from the Lihe River Watershed and other selected areas for comparison (unit mg kg–1).
Table 4. Total heavy metal concentrations in soil from the Lihe River Watershed and other selected areas for comparison (unit mg kg–1).
Heavy MetalsThe Lihe River WatershedBackground Values of Soil in Jiangsu Province SEQ-IIThe Yangtze River Estuary, China
RangeMean ± SDMedianFirst QuartileThird QuartileCVMean
Cd0.28–2.620.57 ± 0.410.450.400.590.720.130.30.19
Cr39.71–79.6654.60 ± 9.0453.9547.6459.610.1777.820079.1
Cu11.68–46.8322.71 ± 9.4421.1716.0922.960.4222.310024.7
Ni16.31–34.1124.35 ± 4.8824.0521.1027.170.2026.75031.9
Pb21.26–141.7239.78 ± 20.6136.7529.1043.130.5226.230023.8
Zn57.14–238.0693.29 ± 38.5082.6176.0893.450.3962.625082.9
ReferencesPresent studyChen et al. [4]Chen et al. [4]Wang et al. [9]
Table 5. The available concentration of six heavy metals in the soil.
Table 5. The available concentration of six heavy metals in the soil.
Heavy MetalsRangeMean Value ± Standard DeviationCoefficient of Variation
Cd0.11–1.540.25 ± 0.2598.16
Cr0.60–1.851.27 ± 0.2620.01
Cu7.45–34.0014.85 ± 7.0947.76
Ni1.16–4.182.64 ± 0.7729.29
Pb12.51–91.8426.84 ± 13.6150.73
Zn9.32–59.1018.92 ± 11.4460.46
Table 6. The correlation between heavy mental concentration and soil properties.
Table 6. The correlation between heavy mental concentration and soil properties.
Heavy MetalspHCECECOMAverage Particle Size
Cd-s0.2220.1230.2870.1320.033
Cr-s0.444 *0.449 **0.3160.266−0.223
Cu-s0.2290.0750.1460.0760.103
Ni-s0.0410.0400.024−0.098−0.131
Pb-s0.0980.1150.2550.182−0.038
Zn-s−0.008−0.117−0.086−0.226−0.113
Cd-sa0.1990.0810.2160.1220.072
Cr-sa0.2150.627 **0.3010.285−0.334
Cu-sa0.1230.0490.0750.0670.184
Ni-sa−0.1200.457 **0.1300.032−0.154
Pb-sa0.0440.1240.2470.194−0.062
Zn-sa−0.062−0.242−0.111−0.309−0.104
Notes: “Cd-s” denotes the concentration of Cd in soil; “Cd-sa” denotes the available concentration of Cd in soil. It was same as “Cr, Cu, Ni, Pb, and Zn.” In addition, “**” indicates correlation was significant at the 0.01 level (two-tailed), while “*” indicates correlation was significant at the 0.05 level (two-tailed). CEC: cation exchange capacity; EC: electrical conductivity; OM: organic matter.
Table 7. Descriptive statistics of Ei and RI for heavy metals in the study area.
Table 7. Descriptive statistics of Ei and RI for heavy metals in the study area.
ParameterHeavy MetalMeanMaxMinSDCVRisk Classification
EiCd131.9604.464.1695.260.722Considerable risk
Cr1.4032.0481.0210.2320.166Low risk
Cu5.09310.502.6192.1170.416Low risk
Ni4.5606.3883.0540.9140.200Low risk
Pb7.59227.044.0573.9330.518Low risk
Zn1.4903.8030.9130.6150.413Low risk
RI(Multi-metal)152.0649.5300-60099.790.656Moderate risk
SD: Standard Deviation; CV: Coefficient of Variation.
Table 8. Acid-soluble fraction of heavy metals in soil and their risk classifications.
Table 8. Acid-soluble fraction of heavy metals in soil and their risk classifications.
Heavy MetalsAcid Soluble Fraction (%)Risk Classification
Cd43.02 ± 10.35High risk
Cr2.36 ± 0.46Low risk
Cu65.19 ± 10.29Very high risk
Ni11.25 ± 4.05Medium risk
Pb67.72 ± 6.97Very high risk
Zn19.79 ± 4.72Medium risk
Table 9. Pearson correlation coefficients for heavy metals (n = 32).
Table 9. Pearson correlation coefficients for heavy metals (n = 32).
Heavy MetalsCdCrCuNiPbZn
Cd1−0.0120.450 a0.1210.841 a0.079
Cr 10.0530.677 a−0.194−0.045
Cu 10.0400.617 a0.763 a
Ni 1−0.013−0.103
Pb 10.345
Zn 1
a Correlation is significant at the 0.01 level (two-tailed).
Table 10. Descriptive statistics for the enrichment factors (EF) of heavy metals and their coefficient of variation (CV) values.
Table 10. Descriptive statistics for the enrichment factors (EF) of heavy metals and their coefficient of variation (CV) values.
Heavy MetalMean ± SDRangeCV
Cd15.56 ± 11.687.92–74.030.75
Cr1.19 ± 0.120.95–1.450.10
Cu1.75 ± 0.800.88–3.880.45
Ni1.53 ± 0.181.29–2.160.12
Pb2.63 ± 1.461.22–9.630.56
Zn2.58 ± 1.281.51–7.840.50
Table 11. Source contribution for different elements by positive matrix factorization (PMF).
Table 11. Source contribution for different elements by positive matrix factorization (PMF).
ElementsProfile Contribution (mg/kg) Percentage Contribution (%)
Factor 1Factor 2Factor 3Factor 4Factor 1Factor 2Factor 3Factor 4
Ca196.842.9 × 10−155.9 × 10−62220.708.141.2 × 10−162.4 × 10−791. 86
Cd0.300.120.152.0 × 10−1352.1121.6226.273.5 × 10−11
Co0.924.717.248.174.3822.3734.4338.82
Cr0.2111.2319.5222.820.4020.8836.3042.43
Cu1.5114.031.484.756.9264.446.8121.82
Fe254.283477.76417.37712.21.4219.4735.9343.18
K69.382557.13587.25308.90.6022.1931.1346.07
Mg5.3 × 10−5491.51982.341624.21.7 × 10−615.8731.7152.43
Mn19.4652.6845.14140.267.5520.4617.5354.46
Ni0.504.339.279.702.1018.1938.9440.77
Pb9.7418.906.013.5325.5149.5015.749.25
Zn056.526.2428.10062.216.8730.93

Share and Cite

MDPI and ACS Style

Chen, L.; Wang, G.; Wu, S.; Xia, Z.; Cui, Z.; Wang, C.; Zhou, S. Heavy Metals in Agricultural Soils of the Lihe River Watershed, East China: Spatial Distribution, Ecological Risk, and Pollution Source. Int. J. Environ. Res. Public Health 2019, 16, 2094. https://doi.org/10.3390/ijerph16122094

AMA Style

Chen L, Wang G, Wu S, Xia Z, Cui Z, Wang C, Zhou S. Heavy Metals in Agricultural Soils of the Lihe River Watershed, East China: Spatial Distribution, Ecological Risk, and Pollution Source. International Journal of Environmental Research and Public Health. 2019; 16(12):2094. https://doi.org/10.3390/ijerph16122094

Chicago/Turabian Style

Chen, Lian, Genmei Wang, Shaohua Wu, Zhen Xia, Zhenang Cui, Chunhui Wang, and Shenglu Zhou. 2019. "Heavy Metals in Agricultural Soils of the Lihe River Watershed, East China: Spatial Distribution, Ecological Risk, and Pollution Source" International Journal of Environmental Research and Public Health 16, no. 12: 2094. https://doi.org/10.3390/ijerph16122094

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