Next Article in Journal
Seasonality of Suicides among Victims Who Are under the Influence of Alcohol
Next Article in Special Issue
Effects of Soil Aggregate Stability on Soil Organic Carbon and Nitrogen under Land Use Change in an Erodible Region in Southwest China
Previous Article in Journal
Partial Smoking Ban and Secondhand Smoke Exposure in Japan
Previous Article in Special Issue
In Situ Persulfate Oxidation of 1,2,3-Trichloropropane in Groundwater of North China Plain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cd, Cu, and Zn Accumulations Caused by Long-Term Fertilization in Greenhouse Soils and Their Potential Risk Assessment

1
Agro-Environmental Protection Institute/Key Laboratory for Environmental Factors Control of Agro-product Quality Safety, Ministry of Agriculture and Rural Affairs, Tianjin 300191, China
2
College of Natural Resources and Environment, South China Agricultural University, Guangzhou 510642, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(15), 2805; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152805
Submission received: 13 June 2019 / Revised: 30 July 2019 / Accepted: 30 July 2019 / Published: 6 August 2019

Abstract

:
The intense management practices in greenhouse production may lead to heavy metal (HM) accumulations in soils. To determine the accumulation characteristics of HM and to evaluate possible HM sources in greenhouse soils, thirty typical greenhouse soil samples were collected in Shouguang District, Shandong Province, China. The results indicate that the Cd, Cu, and Zn concentrations are, respectively, 164.8%, 78.6%, and 123.9% higher than their background values. In the study area, Cd exhibits certain characteristics, such as wide variations in the proportion of its exchangeable form and the highest mobility factor and geo-accumulation index, which are indicative of its high bioavailability and environmental risk. In addition, there is a significant positive correlation between pairs of Cd, P, soil organic carbon, and cultivation age. Combined with principal component analysis, the results indicate the clear effects that agricultural activities have on Cd, Cu, and Zn accumulation. However, Cr, Ni, and Pb have a significant correlation with soil Fe and Al (hydr)-oxides, which indicates that these metals mainly originate from parent materials. This research indicated that long-term intensive fertilization (especially the application of chemical fertilizers and livestock manure) leads to Cd, Cu, and Zn accumulation in greenhouse soils in Shouguang. And the time required to reach the maximum permeable limit in agricultural soils for Cd, Cu, and Zn is 23, 51, and 42 years, respectively, based on their current increasing rates.

1. Introduction

For the last two decades in China, greenhouse vegetable cultivation has increased from 1.395 million hm2 to over 4.1 million hm2 [1,2]. Greenhouse cultivation has contributed significantly to vegetable production. However, this intensive agricultural practice is the cause of several environmental problems due to its high multiple cropping index and productivity. Agricultural activies (e.g., through excessive use of fertilizers, manures, agrochemicals, and irrigation) are possibly the main contributor to heavy metal (HM) accumulation in greenhouse soils [3,4], which ultimately affects human health via food production [5,6].
The Shouguang District is the largest base for greenhouse vegetable production in China, with an annual vegetable production of 4 × 106 tons [7]. Intensive management practices (e.g., more fertilizers and agrochemicals application, high multiple cropping index) in the greenhouse industry in this region have lasted for more than three decades, which has affected HM spatial distributions [8] and led to excessive accumulation of HMs [9] in greenhouse soils. Higher HM accumulations in greenhouse soils than in arable soils were also reported [4,10]. Among the various risk assessment methods used to evaluate potential environmental risks of HMs accumulated in soils, geo-accumulation index (Igeo) is widely used, which indicated higher environmental risks of Cd, Cu and Zn than other HMs in greenhouse soils [11,12,13].
However, Igeo or other assessment methods, such as the ecological risk index and Nemerow multi-factor pollution index, are insufficient indicators of HM bioavailability, mobility, or toxicity. Meanwhile, the majority of previous studies on greenhouse soils in Shouguang or other regions have mainly focused on the total HM content with little emphasis on their chemical speciation. In general, actual environmental risks of HMs depend strongly on the chemical speciation of the HMs in the ambient environment [14,15,16,17]. According to the sequential chemical extraction method proposed by Tessier, et al. [18], the exchangeable fraction represents weakly bound metals and considered as the fraction that is most mobile and readily absorbed by plants. The carbonate-bound fraction is susceptible to changes in pH, whereas the Fe-Mn oxide-bound and organically bound fractions are unstable when exposed to strong reductive and oxidation conditions, respectively. The residual fraction is insoluble under natural conditions, with extremely low bioavailability. Thus, the sum of the exchangeable, carbonate-bound, Fe-Mn oxide-bound, and organically bound fractions is used to characterize the potential bioavailability of a HM, and the sum of the exchangeable and carbonate-bound fractions is known as the mobility factor (MF) [19,20,21]. Thus, the chemical fractions of HMs, which correspond to their chemical binding forms, can provide more useful information when attempting to assess the environmental risks of HMs.
In general, soil HMs derive from parent materials and can be further enriched by human activities [22]. The main anthropogenic sources of HMs in greenhouse soils include chemical fertilizers, manure, sewage sludge, agricultural chemicals and irrigation water [4]. In addition, how to take appropriate actions to reduce HM inputs to soils requires information on the importance and extent of HM contaminations from above-mentioned sources. The most common method to identify main HM sources is to use multivariate statistical analysis, such as correlation analysis and principal component analysis (PCA), to classify metals and selected edaphic parameters (e.g., Al2O3 and Fe2O3) into groups and perform a simple comparison with background values (BVs) of metals [6,23,24,25,26,27].
Therefore, in this study, in order to determine HM (including Cd, Cr, Cu, Ni, Pb, and Zn) pollution characteristics and accurately assess their potential risks and clarify their possible sources in Shouguang greenhouse soils, a field survey was conducted with following objectives: (1) analyze the HM concentrations and fractions in greenhouse soils, (2) assess HM accumulations and potential environmental risks, and (3) identify HM sources for targeted reduction of HM inputs to greenhouse soils in Shouguang.

2. Materials and Methods

2.1. Study Area Description

Shouguang District (36°41′–37°19′ N, 118°32′–119°10′ E) is located in the northern part of the Shandong Province, which is in the lower reaches of the Xiaoqing River and along the southwest coast of Laizhou Bay, Bohai Sea. Shouguang District comprises a total land area of 2.28 × 103 ha, which accounts for 1.43% of the total area of the Shandong Province. Shouguang belongs to the warm temperate monsoon area and has a continental climate, which is characterized by a mean annual precipitation of 608 mm, a frost-free period of 179 to 227 days, and a mean annual temperature of 12.4 °C. The parent soil materials are mainly alluvial and marine sediments, with lake sediments in certain areas, and the soil type is cinnamon soil. Shouguang has unique conditions for vegetable production due to its flat terrain, fertile soil, abundant sunlight, proper temperature, abundant underground freshwater resources, and suitable climate. In the study area, cultivation ages of greenhouses vary from 1 to 30 years. Meanwhile, N-P-K compound fertilizers and livestock manures were frequently applied (Table S1), leading to excessive fertilizer application and unbalanced nutrient inputs in greenhouse soils [28].

2.2. Sampling and Pretreatment

Thirty sites in the southwestern region of Shouguang were randomly selected for greenhouse surface soil (0–20 cm) sampling, mainly distributed in the greenhouse cultivation area of Shouguang (Figure 1). Based on The Technical Specification for Soil Environmental Monitoring (HJ/T 166-2004) [29], large debris, stones, and plastic films were removed from all samples after they were air-dried at room temperature. Samples were then sieved through a 2-mm polyethylene sieve and portions of the sample were ground and sieved through a 0.075-mm polyethylene sieve for further analysis.

2.3. Soil Analysis

Soil pH was measured in a 1:5 (w/v) ratio of soil to 0.01 mol·L−1 CaCl2. For soil organic carbon (SOC) measurement, soil samples were first treated with 2 mol·L−1 HCl to remove inorganic carbon, and then measured by total organic carbon analyzer (muti N/C 3100, Analytik Jena, Jena, Germany).
The total contents of HMs, that is Cd, Cu, Zn, Pb, Ni, and Cr, as well as major elements like Fe, Al, K, P in the soils were determined by digesting the soil samples according to the method recommended by USEPA [30] in a 1:10 (w/v) ratio of soil to HNO3. A soil reference material (GBW07402) was used for quality control during analysis.
The sequential chemical extraction was performed based on the method proposed by Tessier et al. [18], which defined HM fractions as water-soluble/exchangeable form (F1), carbonate-bound form (F2), Fe-Mn oxide-bound form (F3), organically bound form (F4), and residual form (F5).
Cu, Zn, Pb, Ni, Cr, Fe, Al, K, and P concentrations in all extracts (including digestion and sequential chemical extraction) were determined by an Inductively Coupled Plasma Optical Emission Spectrometer (ICP-OES, ELAN DRC-e, Perkin Elmer, Waltham, MA, USA) due to their high concentrations in solutions. Cd concentration was measured by an Inductively Coupled Plasma Mass Spectrometry (ICP-MS, Optima 5300 DV, Perkin Elmer) due to its low concentration in the extracts.

2.4. Parameter Estimation

To assess the risk of HM pollution in greenhouse soils in the Shouguang District, methods recommended by the Environmental Quality Evaluation Standard for Farmland of Greenhouse Vegetable Production (HJ333-2006) [31], which is implemented by the State Environmental Protection Administration of China, as well as the geo-accumulation index ( I geo ) [32] were used.

2.4.1. Environmental Quality Evaluation Standard for Farmland of Greenhouse Vegetables Production (HJ333-2006)

According to pollutant toxicological characteristics, the HJ333-2006 method divides the evaluation indices into two categories: strict and general control indicators. The strict control indicator is evaluated according to the single quality index ( SI ) while the general control indicator is evaluated by the comprehensive quality index ( CI ). Cadmium, As, Hg, Pb, and Cr belong to the strict control indicator while Cu, Zn, and Ni belong to the general control indicator. Table 1 lists the limits of these HMs defined in the HJ333-2006 method.
The evaluation parameters mentioned above are calculated using the following equations:
SI =   single   measured   value   /   single   standard   value
CI   =   [ ( ( Average   SI )   ^   2 + ( Max   SI )   ^   2   ) ] / 2
The single measured value represents the value of a certain HM in soil and the single standard value represents the corresponding HM limit in soil. The contamination categories, based on the HJ333-2006 method, are classified as follows: safety ( SI or CI < 0.7), slight safety (0.7 < SI or CI   < 1), and pollution ( SI or CI   > 1).

2.4.2. Geo-accumulation Index

The geo-accumulation index ( I geo ), originally proposed by Müller et al. [32], is widely used to assess the potential risks of metals that derive from anthropogenic activities. The I geo is calculated with the following equation:
I geo = log 2 [ C n / ( 1 . 5   ×   B n ) ]
where C n denotes the concentration of metal n in the soil and B n represents the corresponding background value. In this study, the metal concentrations in the cinnamon soils of Shouguang were used as the background values (Table 2). A constant of 1.5 was used to analyze natural fluctuations and detect minuscule anthropogenic influences. A positive I geo indicates that there is a large HM contribution from anthropogenic activities. The contamination categories, based on the I geo , are further classified as follows: uncontaminated ( I geo < 0), uncontaminated to moderately contaminated (0 < I geo < 1), moderately contaminated (1 < I geo < 2), moderately to heavily contaminated (2 < I geo < 3), heavily contaminated (3 < I geo < 4), heavily to extremely contaminated (4 < I geo < 5), and extremely contaminated ( I geo > 5).

2.5. Data Analysis

Data analysis was performed using the SPSS 23.0 statistical software (International Business Machines Corporation, Armonk, NY, USA). The sampling map and I geo of each site were charted by ArcGIS 10.2 software (Environmental Systems Research Institute, Redlands, CA, USA). Pearson correlation analysis was used to determine relationships between different HMs, as well as relationships between HM concentrations and selected basic edaphic parameters, such as SOC, total phosphorus (P), total potassium (K), and aluminum and iron (Al2O3 and Fe2O3, i.e., the main constituents of geogenic and pedogenic materials in soils and measured by digesting soil samples), the data obtained are normally distributed and have 3 replications. To further identify associations among and common sources of metals, PCA was performed with a varimax rotation, which can facilitate interpreting the results by minimizing variable numbers with high loading on each component. Variables were not considered to identify source categories until their factor loadings (absolute magnitude) were >1.0.
The data that support the findings of this study are available from the corresponding author and first author on request.

3. Results and Discussion

3.1. HM Concentrations in Greenhouse Soils

The mean concentrations of Cd, Cr, Cu, Ni, Pb, and Zn in the greenhouse soils were 0.21, 70.7, 42.1, 23.0, 18.9 and 144.4 mg·kg−1, respectively (Table 2), which do not exceed the limits set by HJ333-2006 (Table 1). However, there are 5 samples exceed the Cd limit with an excessive rate of 16.1%, while other HMs in all samples do not exceed the limits. And compared with the corresponding natural background values in Shouguang (in mg·kg−1: Cd = 0.08; Cu = 23.6; Zn = 64.5) [10], the mean concentrations of Cd, Cu, and Zn significantly increased by 164.8%, 78.6%, and 123.9%, respectively. On the other hand, their coefficients of variation (CVs) were also relatively high, with a decreasing order of Cd (52.0%) > Cu (35.5%) > Zn (30.1%). Manta, et al. [33] suggested that anthropogenic activities lead to high HM concentrations and CVs, which suggests that Cd, Cu, and Zn accumulations in Shouguang are possibly due to agricultural activities. Meanwhile, the mean concentrations of Cr, Ni, and Pb in greenhouse soils were lower than or similar to their background values (in mg·kg−1: Cr = 66.0; Ni = 28.6; Pb = 26.0). Their low CVs (< 20%) further imply a limited influence from agricultural activities. Similarly, Cd, Cu, and Zn accumulations in greenhouse soils also exist in the Wuwei District and Yellow River Irrigation Region in northwest China [4,34]. However, compared with other regions in China, the mean Cd concentration in the Shouguang greenhouse soils was relatively low, only higher than those in the Nanjing and Wuqing Districts, while Cu and Zn concentrations were higher, except for the Zn in Nanjing and Yunnan (Table 3). In addition, compared with Spain and Turkey (due to rather few data sets of HM contents in greenhouse soils were obtained in other countries), the average Cd concentration, as well as the concentrations of Ni and Pb, in Shouguang were still relatively low while Cu and Zn exhibited opposite trends (Table 3). Overall, Cd, Cu, and Zn show general accumulations in greenhouse soils, with Shouguang exhibiting lower Cd accumulation and higher Cu and Zn accumulations.

3.2. Potential HM Contamination Risks in Greenhouse Soils

To further evaluate the HM pollution status, both the SI and CI indices (Table 4) were calculated. The SI values for six HMs in greenhouse soils from Shouguang were all less than 0.7, decreasing in the order of Cd (0.61) > Zn (0.53) > Cu (0.42) > Ni (0.42) > Pb (0.38). Based on the SI and CI (0.537) results, the contamination categories in Shouguang can be classified as safe and suitable for pollution-free vegetable production. Nonetheless, Cd exhibited a relatively higher environmental risk in the study area. As mentioned above, Cd had an over-standard rate of 16.1% and its maximum concentration was 43.7% higher than the limit (HJ333-2006), which is consistent with results from Zhang, et al. [44] and Wang, et al. [9] conducted in Baiyin (Gansu Province) and Shandong Province in China, respectively.
The I geo results of each sampling site were showed in Figure 2, to clarify the HM accumulations and distributions in the greenhouse soils in Shouguang. In Figure 2, spatial distributions of I geo for Cd, Cu, and Zn are distinctly different from Cr, Ni, and Pb, as well as the fact that Cd, Cu and Zn have relatively high spatial variations and high I geo . Combined with the HM concentration distributions, Cd showed the highest contamination risk, with 30% of the sampling sites characterized as moderately contaminated ( I geo > 1), while Cu and Zn showed a slight contamination risk with over 50% of the sampling sites characterized as uncontaminated to moderately contaminated (0 < I geo < 1). This is similar to the results reported in Tian et al. [12], who found that Cd accumulation was the most likely environmental problem in greenhouse soils of Shandong and Jiangsu provinces. Meanwhile, Cd, Cu, and Zn contaminations caused by anthropogenic sources has also been reported as a serious problem in greenhouse soils in southwestern China [13]. However, I geo values were negative for Cr (except one site), Ni, and Pb at all sampling sites (Figure 2), indicating no Cr, Ni, or Pb pollutions.

3.3. Chemical Fractions of HMs in Greenhouse Soils

To further clarify environmental risks of HMs, the migration ability and bioavailability of HMs that is closely related to their chemical fractions in soils should be concern. Figure 3 shows the results of sequential chemical extraction, which are expressed as percentages of the sum of each chemical fraction.
In Shouguang, the dominant chemical fraction for Cd was the residual fraction (33.7% on average), followed by the exchangeable fraction (24.7% on average), carbonate-bound fraction (23.7% on average), Fe-Mn oxide-bound fraction (4.4% on average), and organically bound fraction (3.47% on average). Compared with other HMs, the proportion of Cd associated with the exchangeable fraction was the highest, ranging from 16.6–36.1%, which may have resulted from Cd’s strong polarization ability that led to its weak adsorption on the interface of soil organic matter [45]. Previous studies have also indicated that a high proportion of exchangeable Cd exists in most agricultural soils (including greenhouse soils), which would easily lead to higher environmental risks [46,47,48,49].
For Cu, the dominant chemical fraction was also the residual fraction (48.8% on average), followed by the organically bound fraction (16.5% on average), Fe–Mn oxide-bound fraction (14.9% on average), carbonate-bound fraction (10.9% on average), and exchangeable fraction (8.86% on average). Compared with other HMs, Cu is easily bound to oxidizable substances (including organic matter and sulfides) [50,51,52]. However, although Cu is a potentially less mobile metal, it can be released into surrounding environment with exposure to strong oxidative, reductive, or acidic soil conditions.
The majority of Zn was bound with the Fe-Mn oxide fraction (41.6% on average), followed by the residual fraction (33.8% on average), carbonate-bound fraction (17.8% on average), organically bound fraction (5.32% on average), and exchangeable fraction (1.48% on average). Zinc was easily bound to clay minerals, carbonates, or hydrous oxides [53,54,55,56,57], so that the highest proportion of Zn bound to Fe-Mn oxide was also found in Shouguang.
As the low-risk HMs, the majority of the Pb (78.7% on average) and Cr (86.9% on average) were found in the residual fraction, which is indicative of their strong combinations within mineral crystal structures and low migration risks in natural conditions [54]. In addition, Ni was mainly associated with the residual (32.5% on average) and Fe–Mn oxide-bound fractions (28.6% on average). Nickel had the larger proportion of residual fraction due to the ease with which it binds to aluminosilicate minerals [58]. Furthermore, Ni belongs to the same group on the periodic table as Fe and Mn, such that Ni can exist in the Fe-Mn oxide lattice via isomorphic substitution or adsorption and precipitation onto the Fe-Mn oxide surface [21].
The sum of four fractions, except residual fraction, for Cd, Cu, Zn, Pb, Cr, and Ni is 66.3%, 51.2%, 66.2%, 21.3%,13.1%, and 67.5%, respectively. These imply that the potential bioaccessibility increases in the order of Cr < Pb < Cu < Zn < Cd < Ni. Meanwhile, mobility factor (MF), sum of exchangeable and carbonate-bound fractions, show increasing order of Cr (0.00%) < Zn (19.2%) < Cu (19.8%) < Pb (20.4%) < Ni (23.4%) < Cd (48.4%). Although Ni has the highest potential bioavailability, it mainly exists in fractions like Fe-Mn oxide-bound and organically bound fractions. However, for Cd, it exhibits a high potential bioavailability and MF, which indicates a clearly high Cd bioavailability in Shouguang and its potential threat to human health and ecosystem safety.

3.4. Identification of the Main Sources Causing HM Accumulations

Since HMs in greenhouse soils could originate from different sources, such as sewage sludge, livestock manure, inorganic fertilizers, pesticides and other agricultural chemicals, irrigation water [22,59], identifying pollution sources is the key to reducing HM accumulations in greenhouse soils.
As shown in Figure 3, proportions of dominant Cd, Cu, and Zn fractions at different sampling sites varied widely. However, for Cr, Ni and Pb, the proportion for each fraction at the different sampling sites mostly showed little change. As previous studies have indicated, the variations in the HM fraction distributions may have mainly resulted from anthropogenic pollution influences. For example, Zong et al. [57] reported that the proportions of the HM fractions varied in a wide range, except for Cr with pedogenic origins. Guo et al. [46] also found that there was a clear spatial heterogeneity in the HM fraction distributions in the different regions of Changchun, which was caused by differences in the pollution sources and soil pH. Additionally, Acosta et al. [53] proposed that different pollution sources affect the HM fraction distributions in Spain. Furthermore, the differences in the pollution sources may result in the differences in proportions of a certain HM fraction [46,60]. Organically bound Cu and Fe–Mn oxide-bound Zn were the major chemical fractions of Cu and Zn in swine and cow manure, respectively [61,62], and higher amounts of exchangeable Cd were mostly found in soils with frequent applications of chemical fertilizers [63]. In Shouguang, it has problems with excessive application of fertilizers [28], mainly as livestock manures and N-P-K compound fertilizer (Table S1). Thus, the varied proportions of organically bound Cu, Fe-Mn oxide-bound Zn, and exchangeable Cd should be related to agricultural activities, i.e., excessive fertilizer (livestock manures and chemical fertilizer) application. Additionally, for Cr, Ni and Pb, high residual fraction proportions and their little spatial heterogeneity in chemical fraction distribution should indicate their parent rock origins [57].
Furthermore, inter-element relationships can provide interesting information on metal element sources [64]. Pearson correlation coefficients were calculated between six HMs and related soil properties (Fe2O3, Al2O3, total K, total P, SOC, and cultivation age), and the results are shown in Table 5. A significant positive correlation (p < 0.05) was observed between Cd and Zn (r = 0.379) while significant positive correlations (p < 0.01) occur between Cu and Zn (r = 0.788), Cr and Ni (r = 0.686), Cr and Pb (r = 0.759), and Ni and Pb (r = 0.918). Significant positive correlations were also observed between Cd, Cu, and Zn with SOC and total P, with the following correlation coefficients: 0.419 (Cd-SOC, p < 0.05), 0.440 (Cu-SOC, p < 0.05), 0.578 (Zn-SOC, p < 0.01), 0.535 (Cd-P, p < 0.01), and 0.402 (Zn-P, p < 0.05). These findings are similar to the results of previous studies [1,65,66]. Furthermore, significant positive correlations (p < 0.01) between Cr, Ni, and Pb with Al2O3 and Fe2O3 were observed, with correlation coefficients of 0.783 and 0.732 for Cr with Al2O3 and Fe2O3, 0.697 and 0.641 for Ni with Al2O3 and Fe2O3, and 0.790 and 0.729 for Pb with Al2O3 and Fe2O3, respectively. Since positive correlations between different HMs imply similar origins [3], it is speculated that Cd, Cu, and Zn in greenhouse soil in Shouguang derive from an identical source while Cr, Ni, and Pb originate from different sources.
Since a single correlation analysis is not sufficient to determine the HM sources [67], PCA was also used. Table S2 lists PCA results and Figure 4 shows principal component loadings. In this study, due to the limitation of an eigenvalue >1, two principal components (PC1 and PC2) were obtained, which cumulatively explained 69.147% of the total variance in the greenhouse soils. After a varimax rotation, PC1 can explain 43.163% of the total variance, which is strongly and positively related to Cr (0.895), Ni (0.840), Pb (0.911), Al (0.925), Fe (0.884) and K (0.682). The PC2 component was characterized by high positive factor loadings for Cd (0.699), Cu (0.791), Zn (0.865) and P (0.670) and can explain 23.984% of the total variance. The high positive factor loading on an identical component strongly implies that these metals have a common source.
PC1 can be representative of geogenic and pedogenic contributions because mean concentrations of Cr, Ni, and Pb did not exceed or approach background values. Numerous previous studies [23,68,69] have shown that concentrations of Cr and Ni in soils are closely related to their contents in parent materials, which is consistent with PCA results of this study. In Shouguang, soil parent material was mostly alluvium, in which Ni and Cr concentrations were mainly affected by pedogenic processes [70]. The study area is far away from main roads and any factories or mining regions. The soils were also kept in a closed greenhouse environment, which prevented them from being contaminated by dry and wet atmospheric Pb deposition.
Thus, Pb in non-polluted soils should mainly be affected by parent materials [71]. Based on correlation analysis results, significant positive correlations of Cr, Ni, and Pb with Al2O3 and Fe2O3 (Table 5) further indicate natural origins for Cr, Ni, and Pb in greenhouse soils since Al2O3 and Fe2O3 are main constituents of geogenic and pedogenic materials. Thus, both PCA and correlation analysis results demonstrate natural origins for Cr, Ni, and Pb.
PC2 can be considered as anthropogenic origin, since mean concentrations of Cd, Cu, and Zn in greenhouse soils were respectively 165%, 78.6%, and 124% higher than their background values. In agricultural soils, these metals may originate from agricultural activities, such as irrigation, fertilization, and pesticide application [8,72]. According to field investigation results (Table S1) and researches conducted by Zeng et al. [28], main chemical fertilizers applied in Shouguang are ammonium dihydrogen phosphate, urea, calcium superphosphate and N-P-K compound fertilizers, while organic fertilizers were mainly livestock manures, all leading to accumulations of SOC and P in soils. Significant positive correlations between Cd, Cu, and Zn with SOC and P (Table 5) further suggested that livestock manures and phosphorus fertilizers are main sources of Cd, Cu, and Zn in Shouguang. Moreover, significant positive correlations between Cd and P and the cultivation age (Table 5) further indicate that Cd accumulation in soils results from long-term phosphorus fertilizer application due to high Cd concentrations (up to 30 mg·kg−1) in phosphorus fertilizers [73,74,75]. In China, livestock manures is a significant HM contributor, which accounts for approximately 55%, 69%, and 51% of total Cd, Cu, and Zn accumulations in agricultural soils, respectively [59]. Bai et al. [41] also reported that Cd, Cu, and Zn accumulation in greenhouse soils in Siping, China are mainly due to application of organic and chemical fertilizers. An investigation of four cities in Shandong Province by Liu et al. [7] showed that fertilizers and agrochemicals are the main sources of Cd, Cu, and Zn accumulation in soils while the concentrations of Cr and Ni in the soils are mainly affected by the parent rocks. Similar results have also been obtained by Nicholson et al. [73], who reported that livestock manure and chemical fertilizers account for approximately 41%, 43%, and 42% of the total Cd, Cu, and Zn accumulated in agricultural soils in England and Wales. Thus, PC2 represents agricultural origin.
PCA results combined with correlation analysis results show that Cd, Cu, and Zn mainly derive from agricultural inputs, especially as livestock manures and chemical fertilizers which was also the main inputs in greenhouse soil in Shouguang (Table S1). Thus, annual input of heavy metals through livestock manures and compound fertilizers to greenhouse soils in Shouguang was calculated (Table S3). Suppose current Cd concentration and its increasing input rate were respectively 0.21 mg·kg−1 and 0.004 mg·kg−1·yr−1 in greenhouse soils in Shouguang, the time required to reach the maximum permeable Cd limit (i.e., 0.3 mg·kg−1 for soil from HJ333-2006) is only approximately 23 years, whereas, for Cu and Zn, the time required is approximately 51 and 42 years, indicating a relative higher environmental risk of Cd in greenhouse soils in Shouguang

4. Conclusions

In Shouguang District, Cd, Cu, and Zn, exhibited varying degrees of accumulations (i.e., increases by 165%, 78.6%, and 124%, respectively) compared with their background values. In addition, Cd and P concentrations increased with the cultivation age and 16.1% of the sampling sites exceeded the standard Cd limit defined in HJ333-2006. Despite these results, the overall environmental quality of the greenhouse soils in Shouguang is safe for production of non-polluted vegetables, based on the criteria defined in HJ333-2006.
Cadmium had the highest MF among six HMs and exhibited variable proportions of exchangeable fraction in Shouguang, indicating its high bioavailability and environmental risk. In addition, both the correlation analysis and PCA results suggest that agricultural activities, especially application of livestock manures and chemical fertilizers, result in Cd, Cu, and Zn accumulations in greenhouse soils, whereas parent rocks mainly determine Cr, Ni, and Pb concentrations.
In Shouguang, the time required to reach the maximum permeable Cd, Cu and Zn limits (HJ333-2006) in agricultural soils is approximately 23, 51 and 42 years respectively if keeping present annual Cd, Cu and Zn input rates from livestock manures and compound fertilizers. To reduce their accumulations in Shouguang greenhouse soils, significant attention should be paid to perform effective pollution management policies, such as formulation and implementation of stricter measures on fertilizer production and application rate. This study provide evidence that current farming practices have significant contributions to heavy metal accumulations in soils in Shouguang district.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1660-4601/16/15/2805/s1, Table S1: Basic information of sampling sites of this study; Table S2: Total variance explained for heavy metal contents based on PCA analysis; Table S3: Annual heavy metal input to greenhouse soils in Shouguang.

Author Contributions

Conceptualization, Y.C.; Data curation, Z.L.; Formal analysis, Z.L.; Funding acquisition, L.W.; Methodology, J.M.; Project administration, Y.L.; Writing—Original draft, Z.L.; Writing—Review & editing, Z.L. and Md. S.I.

Funding

National Key Research and Development Program of China (2016YFD0800304); Natural Science Foundation of China (41701355); Natural Science Foundation of Tianjin City (17JCQNJC08100).

Acknowledgments

The authors would like to thank Ying-xuan Deng, Yu-ling Ma and Long Wang for soil sampling and pretreatment, and are also grateful to the anonymous referees for their helpful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, H.Z.; Li, H.; Wang, Z.; Zhou, L.D. Accumulation Characteristics of Copper and Cadmium in Greenhouse Vegetable Soils in Tongzhou District of Beijing. Procedia Environ. Sci. 2011, 10, 289–294. [Google Scholar] [CrossRef]
  2. Wu, X.; Ge, Y.; Wang, Y.; Liu, D.; Gu, B.; Ren, Y.; Yang, G.; Peng, C.; Cheng, J.; Chang, J. Agricultural carbon flux changes driven by intensive plastic greenhouse cultivation in five climatic regions of China. J. Clean. Prod. 2015, 95, 265–272. [Google Scholar] [CrossRef]
  3. Gil, C.; Boluda, R.; Ramos, J. Determination and evaluation of cadmium, lead and nickel in greenhouse soils of AlmerÃ-a (Spain). Chemosphere 2004, 55, 1027–1034. [Google Scholar] [CrossRef] [PubMed]
  4. Bai, L.Y.; Zeng, X.B.; Su, S.M.; Duan, R.; Wang, Y.N.; Gao, X. Heavy metal accumulation and source analysis in greenhouse soils of Wuwei District, Gansu Province, China. Environ. Sci. Pollut. Res. Int. 2015, 22, 5359–5369. [Google Scholar] [CrossRef] [PubMed]
  5. Gao, Y.F.; Duan, Z.Q.; Heng-Fu, H. Investigation and Evaluation of Heavy Metal Contamination of Greenhouse Soils in Tai Lake Region. Soils 2007, 39, 910–914. [Google Scholar]
  6. Martín, J.A.R.; Ramos-Miras, J.J.; Boluda, R.; Gil, C. Spatial relations of heavy metals in arable and greenhouse soils of a Mediterranean environment region (Spain). Geoderma 2013, 200, 180–188. [Google Scholar] [CrossRef]
  7. Liu, P.; Zhao, H.J.; Wang, L.L.; Liu, Z.H.; Wei, J.L.; Wang, Y.Q.; Jiang, L.H.; Liang, D.; Zhang, Y.F. Analysis of Heavy Metal Sources for Vegetable Soils from Shandong Province, China. J. Integr. Agric. 2011, 10, 109–119. [Google Scholar] [CrossRef]
  8. Liu, Q.; Du, Z.Y.; Shi, Y.X.; Pang, X.G. Spatial distribution charicteristics of soil heavy metals in vegetable growing area based on GIS in Shouguang city, Shandong province. Trans. CSAE 2009, 25, 258–263. (In Chinese) [Google Scholar]
  9. Wang, J.; Mi, W.; Song, P.; Xie, H.; Zhu, L.; Wang, J. Cultivation Ages Effect on Soil Physicochemical Properties and Heavy Metal Accumulation in Greenhouse Soils. Chin. Geogr. Sci. 2018, 28, 717–726. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, P.; Yang, L.; Yu, S.F.; Liu, Z.H.; Wei, J.L.; Wang, X.J.; Wang, Y.Q. Evaluation on Environmental Quality of Heavy Metal Contents in Soils of Vegetable Greenhouses in Shouguang City. Res. Environ. Sci. 2008, 21, 66–71. (In Chinese) [Google Scholar]
  11. Tian, K.; Huang, B.; Xing, Z.; Hu, W. In situ investigation of heavy metals at trace concentrations in greenhouse soils via portable X-ray fluorescence spectroscopy. Environ. Sci. Pollut. Res. 2018, 25, 1–12. [Google Scholar] [CrossRef] [PubMed]
  12. Tian, K.; Hu, W.; Xing, Z.; Huang, B.; Jia, M.; Wan, M. Determination and evaluation of heavy metals in soils under two different greenhouse vegetable production systems in eastern China. Chemosphere 2016, 165, 555–563. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, H.; Huang, B.; Dong, L.; Hu, W.; Akhtar, M.S.; Qu, M. Accumulation, sources and health risks of trace metals in elevated geochemical background soils used for greenhouse vegetable production in southwestern China. Ecotoxicol. Environ. Saf. 2017, 137, 233–239. [Google Scholar] [CrossRef] [PubMed]
  14. Andrew, B.; Cave, M.R.; Joanna, W.; Fordyce, F.M.; Bewley, R.J.F.; Graham, M.C.; Ngwenya, B.T.; Farmer, J.G. Determination of the bioaccessibility of chromium in Glasgow soil and the implications for human health risk assessment. Sci. Total Environ. 2010, 409, 267–277. [Google Scholar]
  15. Duong, T.T.T.; Lee, B.K. Partitioning and mobility behavior of metals in road dusts from national-scale industrial areas in Korea. Atmos. Environ. 2009, 43, 3502–3509. [Google Scholar] [CrossRef]
  16. Gleyzes, C.; Tellier, S.; Astruc, M. Fractionation studies of trace elements in contaminated soils and sediments: A review of sequential extraction procedures. Trends Anal. Chem. 2002, 21, 451–467. [Google Scholar] [CrossRef]
  17. Pyeong-Koo, L.; Byoung-Young, C.; Min-Ju, K. Assessment of mobility and bio-availability of heavy metals in dry depositions of Asian dust and implications for environmental risk. Chemosphere 2015, 119, 1411–1421. [Google Scholar]
  18. Tessier, A.; Campbell, P.G.C.; Bisson, M. Sequential extraction procedure for the speciation of particulate trace metals. Anal. Chem. 1979, 51, 844–851. [Google Scholar] [CrossRef]
  19. Narwal, R.P.; Singh, B.R. Effect of Organic Materials on Partitioning, Extractability and Plant Uptake of Metals in an Alum Shale Soil. Water Air Soil Pollut. 1998, 103, 405–421. [Google Scholar] [CrossRef]
  20. Salbu, B.; Krekling, T. Characterisation of radioactive particles in the environment†. Analyst 1998, 123, 843–850. [Google Scholar] [CrossRef]
  21. Jing, W.Q.; Nan, Z.R.; Liu, X.W.; Wang, S.L.; Li, Y.; Huang, H.; Zhao, C.C. Speciation and bioavailability assessment of cu, zn and ni in argricultural soils from Jinchang, Gansu, China. Environ. Chem. 2010, 29, 220–225. (In Chinese) [Google Scholar]
  22. Martín, J.A.R.; Arias, M.L.; Corbí, J.M.G. Heavy metals contents in agricultural topsoils in the Ebro basin (Spain). Application of the multivariate geoestatistical methods to study spatial variations. Environ. Pollut. 2006, 144, 1001–1012. [Google Scholar] [CrossRef] [PubMed]
  23. Facchinelli, A.; Sacchi, E.; Mallen, L. Multivariate statistical and GIS-based approach to identify heavy metal sources in soils. Environ. Pollut. 2001, 114, 313–324. [Google Scholar] [CrossRef]
  24. Chai, Y.; Guo, J.; Chai, S.; Cai, J.; Xue, L.; Zhang, Q. Source identification of eight heavy metals in grassland soils by multivariate analysis from the Baicheng-Songyuan area, Jilin Province, Northeast China. Chemosphere 2015, 134, 67–75. [Google Scholar] [CrossRef] [PubMed]
  25. Loska, K.; Wiechuła, D. Application of principal component analysis for the estimation of source of heavy metal contamination in surface sediments from the Rybnik Reservoir. Chemosphere 2003, 51, 723–733. [Google Scholar] [CrossRef]
  26. Martínez, J.; Llamas, J.F.; Miguel, E.D.; Rey, J.; Hidalgo, M.C. Soil contamination from urban and industrial activity: Example of the mining district of Linares (southern Spain). Environ. Geol. 2008, 54, 669–677. [Google Scholar] [CrossRef]
  27. Romic, M.; Romic, D. Heavy metals distribution in agricultural topsoils in urban area. Environ. Geol. 2003, 43, 795–805. [Google Scholar] [CrossRef]
  28. Zeng, X.B.; Bai, L.Y.; Li, L.F.; Su, S.M. The status and changes of organic matter, nitrogen, phosphorus and potassium under different soil using styles of Shouguang of Shangdong Province. Acta Ecol. Sin. 2009, 29, 3737–3746. (In Chinese) [Google Scholar]
  29. State Environmental Protection Administration, People’s Republic of China. The Technical Specification for Soil Environmental Monitoring (HJ/T 166–2004); China Environmental Science Press: Beijing, China, 2004. [Google Scholar]
  30. USEPA. Test Methods for Evaluating Soild Waste, Method 3050B: Acid Digestion of Sediments Sludges and Soils (Revision 2); U.S. Government Printing Office (GPO): Washington, DC, USA, 1996. [Google Scholar]
  31. State Environmental Protection Administration, People’s Republic of China. Environmental Quality Evaluation Standard for Farmland of Greenhouse Vegetables Production (HJ333-2006); China Environmental Science Press: Beijing, China, 2007. [Google Scholar]
  32. Müller, M.; Anke, M. Distribution of cadmium in the food chain (soil-plant-human) of a cadmium exposed area and the health risks of the general population. Sci. Total Environ. 1994, 156, 151–158. [Google Scholar] [CrossRef]
  33. Manta, D.S.; Angelone, M.; Bellanca, A.; Neri, R.; Sprovieri, M. Heavy metals in urban soils: A case study from the city of Palermo (Sicily), Italy. Sci. Total Environ. 2002, 300, 229–243. [Google Scholar] [CrossRef]
  34. Kong, X.; Cao, J.; Tang, R.; Zhang, S.; Dong, F. Pollution of intensively managed greenhouse soils by nutrients and heavy metals in the Yellow River Irrigation Region, Northwest China. Environ. Monit. Assess. 2014, 186, 7719–7731. [Google Scholar] [CrossRef] [PubMed]
  35. Topcuoglu, B. Effects of Repeated Applications of Sewage Sludge and MSW Compost on the Bioavailability of Heavy Metals in Greenhouse Soil. Pol. J. Environ. Stud. 2005, 14, 217–222. [Google Scholar]
  36. Sungur, A.; Soylak, M.; Özcan, H. Chemical fractionation, mobility and environmental impacts of heavy metals in greenhouse soils from Çanakkale, Turkey. Environ. Earth Sci. 2016, 75, 334. [Google Scholar] [CrossRef]
  37. Huang, S.W.; Jin, J.Y. Status of heavy metals in agricultural soils as affected by different patterns of land use. Environ. Monit. Assess. 2008, 139, 317–327. [Google Scholar] [CrossRef] [PubMed]
  38. Hu, W.; Chen, Y.; Huang, B.; Niedermann, S. Health Risk Assessment of Heavy Metals in Soils and Vegetables from a Typical Greenhouse Vegetable Production System in China. Hum. Ecol. Risk Assess. Int. J. 2014, 20, 1264–1280. [Google Scholar] [CrossRef]
  39. Shi, J.; Zhang, N.M. The Distributing Character of Heavy Metals and Its Pollution Estimate in Greenhouse soils of Yunnan Province. J. Yunnan Agric. Univ. 2010, 25, 862–867. (In Chinese) [Google Scholar]
  40. Qin, L.L.; Zhang, B.Y.; Wang, W.X. Specificity and potential ecological risk of heavy metals pollution in greenhouse vegetable soils—A case study of Lanzhou city. Guangdong Agric. Sci. 2013, 40, 175–178. (In Chinese) [Google Scholar]
  41. Bai, L.Y.; Zeng, X.B.; Lian-Fang, L.I.; Chang, P.; Shu-Hui, L.I. Effects of Land Use on Heavy Metal Accumulation in Soils and Sources Analysis. J. Integr. Agric. 2010, 9, 1650–1658. [Google Scholar] [CrossRef]
  42. Song, Z.T.; Zhao, Y.J.; Zhou, Q.W.; Liu, X.W.; Zhang, T.L. Applications of Geostatistical Analyses and Stochastic Models to Identify Sources of Soil Heavy Metals in Wuqing District, Tianjin, China. Environ. Sci. 2016, 37, 2756–2762. (In Chinese) [Google Scholar]
  43. Chu, C.J.; Wang, Z.H.; Zhou, J.F.; Peng, S.L.; Yang, F.L. Contents and accumulation characteristics of heavy metals in open and protected vegetable soils. J. Northwest A F Univ. (Nat. Sci. Ed.) 2013, 41, 125–132. (In Chinese) [Google Scholar]
  44. Zhang, S.Q.; Guo, R.Y.; Cao, J.; Nan, Z.R.; Liu, W.T. Characteristic of Nutrient Accumulation and Heavy Metal Pollution Assessment in Greenhouse Soils in Baiyin City. J. Agro-Environ. Sci. 2010, 29, 711–716. (In Chinese) [Google Scholar]
  45. Kumari, J.; Paul, B. Geochemical and environmental risk assessment of hazardous elements in river sediment. Environ. Earth Sci. 2017, 76, 731. [Google Scholar] [CrossRef]
  46. Guo, P.; Xie, Z.L.; Jun, L.I.; Kang, C.L.; Liu, J.H. Relationships between fractionations of Pb, Cd, Cu, Zn and Ni and soil properties in urban soils of Changchun, China. Chin. Geogr. Sci. 2005, 15, 179–185. [Google Scholar] [CrossRef]
  47. Hao, Z.; Wu, C.; Wang, J.; Zhang, X. The effect of simulated acid rain on the stabilization of cadmium in contaminated agricultural soils treated with stabilizing agents. Environ. Sci. Pollut. Res. 2018, 25, 17499–17508. [Google Scholar]
  48. Kosolsaksakul, P.; Farmer, J.G.; Oliver, I.W.; Graham, M.C. Geochemical associations and availability of cadmium (Cd) in a paddy field system, northwestern Thailand. Environ. Pollut. 2014, 187, 153–161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Memoli, V.; Eymar, E.; Garcíadelgado, C.; Esposito, F.; Panico, S.C.; De, A.M.; Barile, R.; Maisto, G. Soil element fractions affect phytotoxicity, microbial biomass and activity in volcanic areas. Sci. Total Environ. 2018, 636, 1099–1108. [Google Scholar] [CrossRef]
  50. Lu, Y.; Zhu, F.; Chen, J.; Gan, H.; Guo, Y. Chemical fractionation of heavy metals in urban soils of Guangzhou, China. Environ. Monit. Assess. 2007, 134, 429. [Google Scholar] [CrossRef] [PubMed]
  51. Yang, H.F.; Wang, Y.B.; Huang, Y.J. Chemical fractions and phytoavailability of copper to rape grown in the polluted paddy soil. Int. J. Environ. Sci. Technol. 2015, 12, 2929–2938. [Google Scholar] [CrossRef]
  52. Zou, X.; Yuan, T.; Zhu, Y.; Zhang, X.; Feng, S.; Shen, Z.; Wang, W. Heterogeneous distribution of copper in different grain size and density fractions of contaminated surface sediment from Nansi Lake (China). Environ. Geol. 2007, 51, 813–820. [Google Scholar] [CrossRef]
  53. Acosta, J.A.; Faz, A.; Kalbitz, K.; Jansen, B.; Martínez-Martínez, S. Partitioning of heavy metals over different chemical fraction in street dust of Murcia (Spain) as a basis for risk assessment. J. Geochem. Explor. 2014, 144, 298–305. [Google Scholar] [CrossRef]
  54. Ma, X.; Zuo, H.; Tian, M.; Zhang, L.; Meng, J.; Zhou, X.; Chang, X.; Liu, Y. Assessment of heavy metals contamination in sediments from three adjacent regions of the Yellow River using metal chemical fractions and multivariate analysis techniques. Chemosphere 2016, 144, 264–272. [Google Scholar] [CrossRef] [PubMed]
  55. Pueyo, M.; Sastre, J.; Hernández, E.; Vidal, M.; Lópezsánchez, J.F.; Rauret, G. Prediction of trace element mobility in contaminated soils by sequential extraction. J. Environ. Qual. 2003, 32, 2054. [Google Scholar] [CrossRef] [PubMed]
  56. Rodríguez, L.; Ruiz, E.; Alonsoazcárate, J.; Rincón, J. Heavy metal distribution and chemical speciation in tailings and soils around a Pb-Zn mine in Spain. J. Environ. Manag. 2009, 90, 1106–1116. [Google Scholar] [CrossRef] [PubMed]
  57. Zong, Y.; Xiao, Q.; Lu, S. Chemical fraction, leachability, and bioaccessibility of heavy metals in contaminated soils, Northeast China. Environ. Sci. Pollut. Res. 2016, 23, 1–8. [Google Scholar]
  58. Li, H.; Ji, H. Chemical speciation, vertical profile and human health risk assessment of heavy metals in soils from coal-mine brownfield, Beijing, China. J. Geochem. Explor. 2017, 183, 22–23. [Google Scholar] [CrossRef]
  59. Luo, L.; Ma, Y.; Zhang, S.; Wei, D.; Zhu, Y.G. An inventory of trace element inputs to agricultural soils in China. J. Environ. Manag. 2009, 90, 2524–2530. [Google Scholar] [CrossRef] [PubMed]
  60. Liu, G.; Wang, J.; Zhang, E.; Hou, J.; Liu, X. Heavy metal speciation and risk assessment in dry land and paddy soils near mining areas at Southern China. Environ. Sci. Pollut. Res. 2016, 23, 8709–8720. [Google Scholar] [CrossRef]
  61. Lv, B.; Xing, M.; Yang, J. Speciation and transformation of heavy metals during vermicomposting of animal manure. Bioresour. Technol. 2016, 209, 397–401. [Google Scholar] [CrossRef]
  62. Zeng, X.; Xiao, Z.; Zhang, G.; Wang, A.; Li, Z.; Liu, Y.; Hua, W.; Zeng, Q.; Liang, Y.; Zou, D. Speciation and bioavailability of heavy metals in pyrolytic biochar of swine and goat manures. J. Anal. Appl. Pyrolysis 2018, 132, 82–93. [Google Scholar] [CrossRef]
  63. Qiutong, X.; Mingkui, Z. Source identification and exchangeability of heavy metals accumulated in vegetable soils in the coastal plain of eastern Zhejiang province, China. Ecotoxicol. Environ. Saf. 2017, 142, 410. [Google Scholar] [CrossRef]
  64. Rodríguez, J.A.; Nanos, N.; Grau, J.M.; Gil, L.; López-Arias, M. Multiscale analysis of heavy metal contents in Spanish agricultural topsoils. Chemosphere 2008, 70, 1085–1096. [Google Scholar] [CrossRef] [Green Version]
  65. Gjoka, F.; Felixhenningsen, P.; Wegener, H.R.; Salillari, I.; Beqiraj, A. Heavy metals in soils from Tirana (Albania). Environ. Monit. Assess. 2011, 172, 517. [Google Scholar] [CrossRef] [PubMed]
  66. 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–74. [Google Scholar] [CrossRef] [PubMed]
  67. Zhou, L.; Yang, B.; Xue, N.; Li, F.; Seip, H.M.; Cong, X.; Yan, Y.; Liu, B.; Han, B.; Li, H. Ecological risks and potential sources of heavy metals in agricultural soils from Huanghuai Plain, China. Environ. Sci. Pollut. Res. 2014, 21, 1360–1369. [Google Scholar] [CrossRef]
  68. Nanos, N.; Martín, J.A.R. Multiscale analysis of heavy metal contents in soils: Spatial variability in the Duero river basin (Spain). Geoderma 2012, 189, 554–562. [Google Scholar] [CrossRef]
  69. Qiang, L.; Liu, J.; Wang, Q.; Yang, W. Source identification and availability of heavy metals in peri-urban vegetable soils: A case study from China. Hum. Ecol. Risk Assess. Int. J. 2016, 22, 1–14. [Google Scholar] [CrossRef]
  70. Mihailović, A.; Budinski-Petković, L.; Popov, S.; Ninkov, J.; Vasin, J.; Ralević, N.M.; Vasić, M.V. Spatial distribution of metals in urban soil of Novi Sad, Serbia: GIS based approach. J. Geochem. Explor. 2015, 150, 104–114. [Google Scholar] [CrossRef]
  71. Salah, E.A.; Turki, A.M.; Mahal, S.N. Chemometric Evaluation of the Heavy Metals in Urban Soil of Fallujah City, Iraq. J. Environ. Prot. 2015, 06, 1279–1292. [Google Scholar] [CrossRef] [Green Version]
  72. Pan, L.B.; Ma, J.; Wang, X.L.; Hou, H. Heavy metals in soils from a typical county in Shanxi Province, China: Levels, sources and spatial distribution. Chemosphere 2016, 148, 248–254. [Google Scholar] [CrossRef]
  73. Nicholson, F.A.; Smith, S.R.; Alloway, B.J.; Carltonsmith, C.; Chambers, B.J. An inventory of heavy metals inputs to agricultural soils in England and Wales. Sci. Total Environ. 2006, 311, 205–219. [Google Scholar] [CrossRef]
  74. Dharma-Wardana, M.W.C. Fertilizer usage and cadmium in soils, crops and food. Environ. Geochem. Health 2018, 219, 1–21. [Google Scholar] [CrossRef] [PubMed]
  75. Huang, Q.Q.; Liu, X.; Zhang, Q.; Liu, X.H.; Qiao, Y.H.; Su, D.C.; Cheng, J.S.; Li, H.F. Evaluating the environmental risk and the bioavailability of Cd in phosphorus fertilizers. Environ. Sci. Technol. 2016, 39, 156–161. (In Chinese) [Google Scholar]
Figure 1. The distribution of sampling sites in the study region (denoted by green circles).
Figure 1. The distribution of sampling sites in the study region (denoted by green circles).
Ijerph 16 02805 g001
Figure 2. Spatial variations of I geo index for HMs in greenhouse soils in Shouguang.
Figure 2. Spatial variations of I geo index for HMs in greenhouse soils in Shouguang.
Ijerph 16 02805 g002
Figure 3. Chemical fractions of HMs (Cd, Cr, Cu, Ni, Pb, and Zn) in greenhouse soils in Shouguang.
Figure 3. Chemical fractions of HMs (Cd, Cr, Cu, Ni, Pb, and Zn) in greenhouse soils in Shouguang.
Ijerph 16 02805 g003
Figure 4. Factor loadings for two principal components after varimax rotation.
Figure 4. Factor loadings for two principal components after varimax rotation.
Ijerph 16 02805 g004
Table 1. The HM limits based on the Environmental Quality Evaluation Standard for Farmland of Greenhouse Vegetables Production.
Table 1. The HM limits based on the Environmental Quality Evaluation Standard for Farmland of Greenhouse Vegetables Production.
Soil pHCdCrCuNiPbZn
mg·kg−1
Soil pH < 6.50.30150504050200
6.5 < Soil pH < 7.50.302001005050250
Soil pH > 7.50.402501006050300
Table 2. Descriptive statistics for soil properties and HM contents in the study area.
Table 2. Descriptive statistics for soil properties and HM contents in the study area.
ParameterRange (mg·kg−1)Mean (mg·kg−1)Median (mg·kg−1)SDCVSkewnessKurtosisBackground Value a (mg·kg−1)
Cd0.09–0.480.210.170.1152.01.120.390.08
Cr37.0–109.170.768.416.122.70.31−0.1766.0
Cu20.8–76.942.139.814.935.50.730.1123.6
Ni16.2–31.223.021.93.7516.30.46−0.6228.6
Pb14.5–25.718.918.22.8915.30.67−0.1226.0
Zn75.9–246.1144.4139.043.530.10.760.5664.5
pH6.6–7.867.397.520.385.1−0.95−0.33-
SOC (%)1.4–3.02.02.13.7118.20.09−0.23-
SOC Soil organic carbon (%), SD standard deviation, CV coefficient of variation (in %). a Background values for soils in Shouguang from Liu et al. [10].
Table 3. Comparisons of HM concentrations in greenhouse soils from various regions.
Table 3. Comparisons of HM concentrations in greenhouse soils from various regions.
Study AreaCdCrCuNiPb ZnReference
mg·kg−1
SpainWestern Almería1.20- a- a38.669.9- a[3]
Granada and Almería1.1050.330.236.068.9133[6]
TurkeyAntalya Aksu- a- a25.016.645.088.0[35]
Çanakkale1.0785.546.755.519.684.6[36]
ChinaNorth China Plain0.6487.929.9- a29.482.4[37]
Nanjing, Jangsu0.15- a52.2- a53117[38]
Yunnan0.4089.946.618.319.848.6[39]
Lanzhou, Gansu0.2665.5-a12.323.9- a[40]
Siping, Jilin0.4767.53725.218.087.7[41]
Wuwei, Gansu0.4267.933.928.620.885.3[4]
Wuqing, Tianjin0.1462.524.532.124.371.6[42]
Henan- a54.432.037.7- a49.0[43]
Shouguang, Shandong0.2170.742.123.018.9144.4This study
a Not available.
Table 4. The single and comprehensive quality indices for HMs in greenhouse soils in Shouguang.
Table 4. The single and comprehensive quality indices for HMs in greenhouse soils in Shouguang.
ElementCdCrCuNiPbZn
Single Quality IndexMean0.6150.3150.4210.4190.3780.526
SD0.3350.0870.1490.0950.0580.184
Range0.224–1.500.148–0.4750.208–0.7690.269–0.6240.291–0.5130.253–0.984
Over-standardrate (%)16.1300000
Comprehensive Quality Index0.537
Table 5. The Pearson Correlation coefficients (r) between the HM contents and soil property parameters in greenhouse soils in Shouguang.
Table 5. The Pearson Correlation coefficients (r) between the HM contents and soil property parameters in greenhouse soils in Shouguang.
ItemCdCrCuNiPbZnFe2O3Al2O3KPSOCCultivation Age
Cd1
Cr−0.0481
Cu0.3360.1391
Ni−0.0220.686 **0.1341
Pb−0.0970.759 **0.0750.918 **1
Zn0.379 *0.0130.788 **0.2210.0941
Fe2O3−0.0430.732 **0.0380.641 **0.729 **−0.0281
Al2O3−0.1180.783 **0.0300.697 **0.790 **0.0280.923 **1
K−0.0050.673 **0.3230.3300.468 **0.2230.544 **0.575 **1
P0.535 **−0.1290.194−0.080−0.1830.402 *−0.136−0.1410.0651
SOC0.419 *0.2240.440 *0.2310.1420.578 **0.0900.1440.379 *0.566 **1
Cultivation age0.647 **−0.111−0.097−0.200−0.214−0.0250.001−0.091−0.1320.560 **0.0721
* Correlations are significant at the 0.05 level. ** Correlations are significant at the 0.01 level.

Share and Cite

MDPI and ACS Style

Liao, Z.; Chen, Y.; Ma, J.; Islam, M.S.; Weng, L.; Li, Y. Cd, Cu, and Zn Accumulations Caused by Long-Term Fertilization in Greenhouse Soils and Their Potential Risk Assessment. Int. J. Environ. Res. Public Health 2019, 16, 2805. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152805

AMA Style

Liao Z, Chen Y, Ma J, Islam MS, Weng L, Li Y. Cd, Cu, and Zn Accumulations Caused by Long-Term Fertilization in Greenhouse Soils and Their Potential Risk Assessment. International Journal of Environmental Research and Public Health. 2019; 16(15):2805. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152805

Chicago/Turabian Style

Liao, Zhongbin, Yali Chen, Jie Ma, Md. Shafiqul Islam, Liping Weng, and Yongtao Li. 2019. "Cd, Cu, and Zn Accumulations Caused by Long-Term Fertilization in Greenhouse Soils and Their Potential Risk Assessment" International Journal of Environmental Research and Public Health 16, no. 15: 2805. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152805

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