Next Article in Journal
Resistance Mechanisms to 2,4-D in Six Different Dicotyledonous Weeds Around the World
Previous Article in Journal
Marginal Grapevine Germplasm from Apulia (Southern Italy) Represents an Unexplored Source of Genetic Diversity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Potentially Toxic Element Content in Arid Agricultural Soils in South Iran

by
Tahereh Moghtaderi
1,
Ata Shakeri
2 and
Andrés Rodríguez-Seijo
3,4,*
1
Landscape and Green Spaces Organization of Shiraz Municipality Shiraz, Shiraz 71366-18947, Iran
2
Department of Applied Geology, Faculty of Earth Science, Kharazmi University, 49 Mofatteh Avenue, P.O.Box 15614, Tehran, Iran
3
Interdisciplinary Centre of Marine and Environmental Research (CIIMAR), University of Porto, Terminal de Cruzeiros do Porto de Leixões, Av. General Norton de Matos s/n, 4450-208 Matosinhos, Portugal
4
Department of Biology, Faculty of Sciences, University of Porto, Rua do Campo Alegre, 4169-007 Porto, Portugal
*
Author to whom correspondence should be addressed.
Submission received: 16 March 2020 / Revised: 5 April 2020 / Accepted: 10 April 2020 / Published: 14 April 2020
(This article belongs to the Section Soil and Plant Nutrition)

Abstract

:
Potentially toxic elements (PTE) are considered to be dangerous threats, both for human health and the environment. Here, the contamination level, sources, and ecological risks posed by PTE were investigated in 19 topsoils from agricultural lands in the Bandar Abbas County on the southern coast of Iran. The soil fraction <63 μm was used for the analysis of the pseudototal contents of PTE (arsenic, cadmium, cobalt, chromium, copper, lead, manganese, molybdenum, nickel, vanadium, and zinc). The results were analyzed using principal component analysis (PCA) and Pearson’s correlations, different pollution indices [enrichment factor (EF), contamination factor (Cf), and modified degree of contamination (mCd), and the potential ecological risk index (PERI). The results of the pollution indices showed that the contamination levels of the studied agricultural soils are low to moderate, while the average PERI values indicate that the ecological risks of PTE range from low to high in the soils of the studied area. Correlations and PCA analyses suggest that studied elements have three possible origins: geogenic (chromium, cobalt, copper, manganese, vanadium, and zinc), anthropogenic (cadmium, nickel, phosphorus and lead) and atmospheric (arsenic) sources. These findings are helpful for the long-term assessment of PTE in the Bandar Abbas County.

1. Introduction

Soil plays an essential role in human development as an environment that supports life, a natural buffer of contaminants, a substrate for food production, and a source of biodiversity. However, soils are also highly impacted by anthropogenic activities as a result of industrialization and human development and are sometimes neglected in the environmental guidelines of national governments [1,2,3,4,5,6]. Agricultural soils, in particular, are essential for humans due to their role in food production, but are sometimes highly degraded at the physical, chemical, and biological levels by inorganic (potentially toxic elements, PTE) and organic contaminants [5,6]. These contaminants reach the soils due to agricultural activities (fertilization, pesticides, etc.) and/or the agricultural area’s proximity to sources of soil contamination, such as industrial or urbanized areas, energy production infrastructure or landfills. Moreover, agricultural soil contamination is not only an environmental issue, since the potential transfer of contaminants from soils to crops can also pose a safety concern based on their risk to humans [7,8].
Iran has experienced widespread development in the last four decades, including fast urbanization of large areas, industrialization, and intensive agriculture in several regions. These changes have sometimes been accompanied by overlooked environmental depreciation [9,10,11,12,13]. In this sense, Bandar Abbas is one of Iran’s most industrialized regions, hosting a major shipping hub, zinc-smelting industry, and an important petrochemical industry, including several petroleum industries and oil fields [9]. At the same time, it is an important agricultural area with the cultivation of crops such as wheat, barley, and sugar beet.
Soils in agricultural areas can be polluted by contaminants from industrial as well as urban developments, creating a risk for humans through the consumption of food cultivated in these soils, as indicated by Doabi et al. [12,13] for other regions of Iran. Although some works have highlighted the influences of anthropogenic activities on street dust [14] and industrial soils [9] on the study region, available information about the contamination levels and potential sources of PTE in the agricultural soils is still scarce.
This work is the first detailed research about the PTE accumulation in the farmlands of the Bandar Abbas region, an important region in Iran with more than 0.6 million inhabitants. The objectives of this study were: (1) to determine the main physicochemical soil properties (soil texture, pH, organic carbon, and cation exchange capacity), and the levels of several PTE [arsenic (As), cadmium (Cd), cobalt (Co), chromium (Cr), copper (Cu), manganese (Mn), molybdenum (Mo), nickel (Ni), lead (Pb), vanadium (V), and zinc (Zn)] in the agricultural soils from the study area; (2) to estimate the degree of contamination through several environmental indicators (contamination factor, enrichment factor, potential ecological risk index, and the modified degree of contamination), and (3) to identify the possible origin (natural or anthropogenic) of the studied elements.

2. Materials and Methods

2.1. Study Area

This study was performed in Bandar Abbas County South District (BAC); which covers 4063 km2 and is located at 27° 00′–27° 30′ N, 55° 54′–56° 29′ E [9]. The city of Bandar Abbas is the capital of the Hormozgan province and is the seventh most populated city in Iran. Figure 1 provides the location map and the selected samples in the study area. The climate of the study area is semi-arid to arid with an average precipitation of 250 mm and an annual temperature of 27.2 °C [9,15]. It is recognized as one of the agricultural and industrial pole in Iran, with significant economical activities, such as fishing and commercial ports, petrochemical and steel, zinc and aluminum complexes, cement factories, etc. [10,16]. The major fertilizers used in this area are superphosphate, ammonium phosphate, and potassium chlorate fertilizers. The geology of the study area is mainly composed by limestone, shale, marl, sandstone, and conglomerate, while the soils are mainly entisoils and arid soils [9].

2.2. Soil Collection and Sample Characterisation

A total of 19 topsoil samples (0–10 cm sampling depth) were gathered from agricultural lands from BAC. Sampling sites were selected in such a way as to cover an impacted area based on recognized pollution sources such as industrial complex zones, landfills, etc. following the methodology developed for previous work with industrial areas [9]. Details about sampling methodology were previously indicated for industrial areas [9].
After the collection and once in the laboratory, samples were air-dried, sieved (2 mm nylon sieve) to eliminate impurities (e.g., stones, agricultural plastics, roots, or cigarette butts) following the methodology indicated by Keshavarzi et al. [14], ground with an agate mortar and pestle, and finally passed through a 63-μm nylon sieve. The soil fraction <63 μm was selected for further metal analysis [17]. According to Zheng et al. [18] and Doabi et al. [12], the health risks associated with fine particles are more significant than those for coarser fractions. The physical and chemical properties of the studied samples were measured using standard methods for the <2 mm fraction: soil pH (1:5, soil:water ratio), particle size distribution determined by Pippete method [19], organic carbon content was measured by titration method [20], and the cation exchange capacity (CEC) [21].
The contents of the studied PTE [aluminum (Al), As, Cd, Co, Cr, Cr, iron (Fe), Mn, Mo, Ni, phosphorous (P), Pb, scandium (Sc), V, and Zn] were measured using a mixture of different acid combinations (HF, HClO4, HNO3, and H2O2) following the methodology indicated in Moghtaderi et al. [9]. The concentrations of the selected elements were evaluated using ICP–MS (Agilent, 7700x, Santa Clara, CA, USA). Details about internal quality standards were previously indicated [9].

2.3. Determination of Pollution Indices

In accordance with several previous studies [17,22,23,24], different environmental indices were used to assess soil contamination on the basis of PTE concentrations, including the contamination factor (Cf), enrichment factor (EF), the modified degree of contamination (mCd), and the potential ecological risk index (PERI). These indexes are widely used in different soils since can be used to assess the presence/intensity of a contaminant deposition on sediments or surface soils, allowing a quick and robust assessment of soil contamination e.g., [8,9,13,17,22,23,24]. As indicated by [25], these calculation methods allow converting numerical results into broad descriptive bands of pollution ranging from low to high intensity.
Cf was employed to evaluate the potential contamination of studied soils by individual elements [26]. It was calculated according to Hakanson (1980) using the following equation (Equation (1)):
Cf = Mx/Mb
where Mx and Mb are the contents of each PTE in the soil and reference environment, respectively is divided into four categories (Table S1).
Abrahim [27] presents a modified and comprehensive version of Hakanson’s [26] equation for the computation of the modified degree of contamination (mCd), following Shakeri et al. [28] and Krzysztof et al. [22] (Equation (2)):
mCd = ∑Cf/n
where n is the number of 11 PTE. The different degrees of mCd are shown in Table S1.
PERI was used to determine the toxicity of PTE in the soils [26,29]). In accordance with this methodology, the potential ecological risk factor (Eir) of each PTE and the PERI of multiple elements can be calculated by Equation (3):
PERI = ∑ni=1 Eir = ∑ni=1 Cf × Tir
where Cf and Tr are the contamination factor and toxic-response factor, respectively, for each PTE (i.e., each element of i). The toxic-response factors for As, Cd, Cr, Cu, Ni, Pb and Zn are 10, 2, 5, 5, 30, 5, and 1. PERI indicates the sensitivity at the biological level to the PTE and explains the potential ecological risk caused by the overall contamination [30,31]. The degree of ecological risk of each element (Eir) and the Potential Ecological risk Index (PERI) classes are shown in Table S1 [30,31,32,33].
The enrichment factor (EF) relative to background was computed to measure the fractions of PTE in soil that originated from human sources/anthropogenic vs. lithogenic. It is supposed that the contents of Al, Fe, Sc, and Ti in the soil can be described solely by geogenic source [23,34,35,36]. In this study, the EF was measured by the following equation (Equation (4)):
EF = [M]/[Sc]soil/[M]/[Sc]background
where M is the element under consideration and Sc is a reference element. The range of EF from 0.5 to 2 regarded as natural variation, while ratios higher than 2 represent enrichment due primarily to human inputs [31,37]. Finally, the EF values for each toxic element are classified in five categories that are shown in Table S1.

2.4. Data Analysis

Pearson’s correlation analysis was carried out using R (v. 3.1.3, R Core Team, Vienna, Austria) to determine the associations between the different PTE. A Principal component analysis (PCA) was also carried out using SPSS 19.0 (Armonk, NY: IBM Corp.) to explore potential sources of toxic elements in agricultural soil samples from the study area. The PCA is an efficient statistical technique used to identify probable sources of PTE [38]. Varimax rotation with Kaiser Normalization was used to simplify the data and their interpretation [10,11].

3. Results and Discussion

3.1. Soil Physicochemical Properties and PTE Levels in Studied Soils

Physicochemical properties of studied soils are shown in Table 1. Briefly, soils are loam and silt loam according to the USDA classification, with slightly to moderately alkaline soils. Organic matter content was lower (up to 1.13%), while CEC values were lower. The values are similar to the samples from industrial soils previously characterized [9], although, in the agricultural soils, more alkaline values were observed but with less organic carbon content.
Regarding the studied PTE, only Cr and Ni showed higher levels that exceeded the Iranian Environmental Quality Standards for agricultural lands (200 and 110 mg·kg−1, respectively) [13], and also the Canadian soil quality guidelines (CSQG) for agricultural soils (64 and 45 mg·kg−1, respectively) [39]. Although these elevated levels can be related to the geogenic origin, according to the higher background levels of these elements (74.5 and 77 mg·kg−1, respectively), it’s also widely known that Cr and Ni can be related to fertilizers, limestone and manures [40,41,42,43]. It was interesting also to advise that our soils have higher levels of P, probably by the widespread use of sulfate and phosphate fertilizers as indicated by Keshavarzi al. for Iranian soils [44] that can reinforce the origin of Cr and Ni in our soils.
The PTE levels obtained in this studied were compared with data reported by other studies carried out for agricultural soils in Iran (Table 2). When compared with the results of previous studies in Iran e.g., [13,40,41,43,44,45,46,47], average values are similar or lower to indicate by them, probably due to some of these studies were carried out near to industrial areas e.g., [40,41]. Results showed significant enrichment of Cr and Ni in studied agricultural lands, although with similar values as indicated by Doabi et al. [13] for agricultural areas in Kermanshah province, Western Iran. When compared with other agricultural areas from arid/semiarid environments e.g., [3,48,49,50,51,52], our studied soils showed similar values to studies from other regions. In general, the contamination level of our samples was low, and higher contents seem to be related to inputs of fertilizers and pesticides [42,53], with a low influence from industrial or urban sources.

3.2. Assessment of Potentially Toxic Elements Pollution

The EF, Cf, mCd, and PERI indexes were used to determine the possibility of PTE pollution in the studied soils. The ranges of average EF values for the agricultural areas for Ni, Cr, and Cd show that the soils are moderately enriched by these elements (Figure 2). In contrast, the average EF values for the other elements show that the agricultural soil samples were minimally enriched. PTE with low EF values (<2) in the agricultural soil samples indicate that those elements originated from geogenic causes. Geological and lithological maps support this since, in the northeast (NE) of the study area, tuffs and shale were observed [9]. Higher concentrations of PTE have also been observed for similar environments like our study [31], where Cr and Cd showed moderate to severe enrichment.
Contamination factor (Cf) was determined based on the background concentrations of PTE in BAC soils (Figure 3). The highest Cf values were observed in agricultural soil samples collected from (AG5), (AG8), and (AG16). The mCd values in agricultural soil samples varied from 0.83 to 1.63, with an average of 1.22, illustrating that the degree of contamination of these soil samples is very low.
PERI indicates the sensitivity of diverse functional groups to toxic substances and also quantifies the PERI posed by the PTE [32,54,55]. The results for the measured selected elements in BAC soil samples are given in Table 3. The mean values of the PTE in the samples were, in decreasing order, Cd > As > Ni > Cu > Pb > Co > Cr > Zn. Only the values for Cd were in the range of moderate to considerable ecological risk. Based on these data, five of the sampling locations in the agricultural area pose a high potential ecological risk (potential ecological risk index >130).

3.3. Multivariate Statistical Analysis

Potential correlations between the studied PTE and the soil physicochemical parameters were identified through Pearson’s correlation analysis (Table 4). A significate and positive correlation were found between Sc and Al, Fe, Cd, Cr, Cu, Mn, Mo, Ni, V, and Zn. Except for As and Pb, all of the PTE concentrations are highly associated with each other, indicating probable similar origins. As indicated by Peris et al. for agricultural soils in Spain [53], Cr and Mn are also highly correlated with Fe, Co, and Ni suggesting an affinity between them (r > 0.55; p < 0.01). Element concentrations were not well correlated with physicochemical properties, suggesting that the effect of the physicochemical properties on the total contents of studied elements is not very important in the study area. Similar results were also indicated by Manta et al. [56] or Acosta et al. [57] for neutral–alkaline soils, like our samples (soil pH = 7.3–8.45).
The inter-element relationships give information useful for identifying the sources of PTE in environmental matrices [58]. In the present study, the PCA was used to identify different element sources and anthropogenic activities (Table 5).
According to the PCA analyses, three components were described to classify the possible sources of elements. These components explain 77% of the variance in the data. Component 1, which explains 42.9% of the variance, was significantly loaded by Fe, Co, Al, Cu, Cr, Sc, V, Mn, and Zn. Component 2 explained 19.45% of the variance and was loaded by Cd, Ni, P, and Pb. Finally, component 3 explained 14% of the variance and was composed of clay, Pb and As. Also, Cr received approximately similar loading factors in PC1 and PC2. Hence, it exhibited a quasi-independent behaviour in the agricultural soil samples, indicating that Cr might have mixed sources.
In general, it seems that a “geogenic factor” is the most significant source of PTE in the studied soils, as indicated by component 1 and previously indicated by Moghtaderi et al. (2018) for industrial areas [9]. Component 2 could be called as “anthropogenic factor” and was constituted by Cd, Ni, P, and Pb with a source from anthropogenic origins, such as industrial/transport emissions or agricultural practices (e.g., agricultural machinery, pesticides, and fertilizers) [13,36,42,48,59]. The presence of P with these PTEs in this component reinforces the possibility that agricultural malpractice could be the source of these elements, such as phosphate fertilizers and uncontrolled uses of pesticides [8,48,51,59,60]. Although Cd levels were lower, uncontrolled use of phosphate fertilizers could be an explanation for this common origin on Component 2 [59,61], as also observed by Cheraghi et al. [60] for agricultural farmlands in the Hamadan province (Western Iran) or by Micó et al. [59] for agricultural soils in Spain. Interestingly, component 3 was formed by clay and As. This combination could have an “atmospheric origin”, as indicated by Li et al. [48] for arid agricultural soils in China, where arsenic is related to fine soil fractions (e.g., clay fraction) with an origin of wind erosion.

4. Conclusions

The concentrations of 11 PTEs and their pollution level in agricultural soils from BAC were quantified in this study. In general, the contamination level of our samples was low, and seems to be related to the agricultural activities, with a low influence from industrial or urban sources. According to the different analysis, there are mixed natural/geogenic and anthropogenic sources in the studied soils, mainly to the malpractices of application of phosphate fertilizers. In general, chromium and nickel are priority PTEs for control in the studied soils. Furthermore, the potential impact of these PTEs on plant uptake or water resources should be considered for further works. Finally, this work provides a basis for the implementation of public policies to reduce the input of PTE in agricultural soils, and therefore, to reduce their long term accumulation, both in soils and in their possible uptake by biota.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4395/10/4/564/s1, Table S1: Indices of pollution assessment for studied potentially toxic elements.

Author Contributions

Conceptualization, T.M. and A.S.; Methodology, T.M. and A.S.; Validation and Formal Analysis, T.M., A.S., and A.R.-S.; Investigation, T.M. and A.S.; Resources, T.M. and A.S.; Writing—Original Draft Preparation, T.M. and A.S.; Writing—Review & Editing, T.M., A.S., and A.R.-S.; Supervision, A.S.; Funding Acquisition, A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The APC was funded by ARS.

Acknowledgments

We thank the research committee of Science and Research Branch, Islamic Azad University, Tehran, Iran and the Environment Protection organization of Hormozgan province for financial support and providing the required equipment. This research was also supported by the Strategic Funding UID/Multi/04423/2019 (CIIMAR) through national funds provided by FCT (Foundation for Science and Technology) and the European Regional Development Fund (ERDF), in the framework of the PT2020 Partnership Agreement. A.R.-S would like to express their thanks to the FCT (Foundation for Science and Technology) for their individual research contract (CEECIND/03794/2017).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Huang, S.S.; Liao, Q.L.; Hua, M.; Wu, X.M.; Bi, K.S.; Yan, C.Y.; Zhang, X.Y. Survey of heavy metal pollution and assessment of agricultural soil in Yangzhong district, Jiangsu Province, China. Chemosphere 2017, 67, 2148–2155. [Google Scholar] [CrossRef]
  2. Huo, X.N.; Li, H.; Sun, D.F.; Zhou, L.D.; Li, B.G. Multi-scale spatial structure of heavy metals in agricultural soils in Beijing. Environ. Monit. Assess. 2010, 164, 605–616. [Google Scholar]
  3. Rodríguez-Martín, J.A.; Ramos-Miras, J.J.; Boluda, R.; Gil, C. Spatial relations of heavy metals in arable and green house soils of a Mediterranean environment region (Spain). Geoderma 2013, 200–201, 180–188. [Google Scholar] [CrossRef]
  4. Benhaddya, M.L.; Hadjel, M. Spatial distribution and contamination assessment of heavy metals in surface soils of HassiMessaoud, Algeria. Environ. Earth Sci. 2014, 71, 1473–1486. [Google Scholar] [CrossRef]
  5. Ullah, Z.; Naz, A.; Saddique, U.; Khan, A.; Shah, W.; Muhammad, S. Potentially toxic elements concentrations and human health risk assessment of food crops in Bajaur Agency, Pakistan. Environ. Earth Sci. 2017, 76, 482. [Google Scholar] [CrossRef]
  6. Muhammad, S.; Ullah, R.; Jadoon, I.A.K. Heavy metals contamination in soil and food and their evaluation for risk assessment in the Zhob and Loralai valleys, Baluchistan province, Pakistan. Microchem. J. 2019, 149, 103971. [Google Scholar] [CrossRef]
  7. Tóth, G.; Hermann, T.; Da Silva, M.R.; Montanarella, L. Heavy metals in agricultural soils of the European Union with implications for food safety. Environ. Int. 2016, 88, 299–309. [Google Scholar] [CrossRef]
  8. Khalifa, M.; Gad, A. Assessment of Heavy Metals Contamination in Agricultural Soil of Southwestern Nile Delta, Egypt. Soil Sediment Contam. 2018, 27, 619–642. [Google Scholar] [CrossRef]
  9. Moghtaderi, T.; Mahmoudi, S.; Shakeri, A.; Masihabadi, M.H. Heavy metals contamination and human health risk assessment in soils of an industrial area, Bandar Abbas–South Central Iran. Hum. Ecol. Risk Assess. 2018, 24, 1058–1073. [Google Scholar] [CrossRef]
  10. Moghtaderi, T.; Mahmodi, S.; Shakeri, A.; Masihabadi, M.H. Contamination evaluation, health and ecological risk index assessment of potential toxic elements in the surface soils (case study: Central Part of Bandar Abbas County). J. Soil Water Conserv. 2019, 8, 51–65. [Google Scholar]
  11. Moghtaderi, T.; Aminiyan, M.M.; Alamdar, R.; Moghtaderi, M. Index-based evaluation of pollution characteristics and health risk of potentially toxic metals in schools dust of Shiraz megacity, SW Iran. Hum. Ecol. Risk Assess. 2019, 25, 410–437. [Google Scholar] [CrossRef]
  12. Doabi, S.A.; Karami, M.; Afyuni, M.; Yeganeh, M. Pollution and health risk assessment of heavy metals in agricultural soil, atmospheric dust and major food crops in Kermanshah province, Iran. Ecotoxicol. Environ. Saf. 2018, 163, 153–164. [Google Scholar] [CrossRef] [PubMed]
  13. Doabi, S.; Karami, M.; Afyuni, M. Heavy metal pollution assessment in agricultural soils of Kermanshah province, Iran. Environ. Earth Sci. 2019, 78, 70. [Google Scholar] [CrossRef]
  14. Keshavarzi, B.; Abbasi, S.; Moore, F.; Mehravar, S.; Sorooshian, A.; Soltani, N.; Najmeddin, A. Contamination Level, Source Identification and Risk Assessment of Potentially Toxic Elements (PTEs) and Polycyclic Aromatic Hydrocarbons (PAHs) in Street Dust of an Important Commercial Center in Iran. Environ. Manag. 2018, 62, 803–818. [Google Scholar] [CrossRef] [PubMed]
  15. Bakhtiari, B.; Kermanib, M.N.; Bordbar, M.H. Rain Gauge Station Network Design for Hormozgan Province in Iran. Desert 2013, 18, 45–52. [Google Scholar]
  16. Dadras, M.; Shafri, H.Z.; Ahmad, N.; Pradhan, B.; Safarpour, S. Spatio-temporal analysis of urban growth from remote sensing data in Bandar Abbas city, Iran. Egypt. J. Remote Sens. Space Sci. 2015, 18, 35–52. [Google Scholar] [CrossRef] [Green Version]
  17. Keshavarzi, B.; Tazarvi, Z.; Rajabzadeh, M.A.; Najmeddin, A. Chemical speciation, human health risk assessment and pollution level of selected heavy metals in urban street dust of Shiraz, Iran. Atmos. Environ. 2015, 119, 1–10. [Google Scholar] [CrossRef]
  18. Zheng, N.; Liu, J.; Wang, Q.; Liang, Z. Health risk assessment of heavy metal exposure to street dust in the zinc smelting district, Northeast of China. Sci. Total Environ. 2010, 408, 726–733. [Google Scholar] [CrossRef]
  19. Gee, G.W.; Bauder, J.W. Particle-size Analysis. In Methods of Soil Analysis; Klute, A., Ed.; SSSA: Madison, WI, USA, 1986; pp. 383–411. [Google Scholar]
  20. Gaudette, H.E.; Flight, W.R.; Toner, L.; Folger, D.W. An inexpensive titration method for the determination of organic carbon in recent sediments. J. Sediment. Petrol. 1974, 44, 249–253. [Google Scholar]
  21. Hendershot, W.H.; Duquette, M. A simple barium chloride method for determining cation exchange capacity and exchangeable cations. Soil Sci. Soc. Am. J. 1986, 50, 605–608. [Google Scholar] [CrossRef]
  22. Krzysztof, L.; Wiechula, D.; Korns, I. Metal contamination of farming soils affected by industry. Environ. Int. 2004, 30, 159–165. [Google Scholar]
  23. Li, P.; Lin, C.; Cheng, H.; Duan, X.; Lei, K. Contamination and health risks of soil heavy metals around a lead/zinc smelter in southwestern China. Ecotoxicol. Environ. Saf. 2015, 113, 391–399. [Google Scholar] [CrossRef] [PubMed]
  24. Xiang, J.; Wang, J.; Chen, L.; Ling, Z.; Han, J.; Li, Q.; Wang, Q. Distribution, Source Identification, and Assessment of Potentially Toxic Elements in the Sediment Core from the Estuarine Region of the Golmud River to the Qarhan Salt Lake, Qinghai, China. Minerals 2019, 9, 506. [Google Scholar] [CrossRef] [Green Version]
  25. Abrahim, G.M.S.; Parker, R.J. Assessment of heavy metal enrichment factors and the degree of contamination in marine sediments from Tamaki Estuary, Auckland, New Zealand. Environ. Monit. Assess. 2008, 136, 227–238. [Google Scholar] [CrossRef]
  26. Hakanson, L. An ecological risk index for aquatic pollution control. A sedimentological approach. Water Res. 1980, 14, 975–1001. [Google Scholar] [CrossRef]
  27. Abrahim, G.M.S. Holocene Sediments of Tamaki Estuary. Characterization and Impact of Recent Human Activity on an Urban Estuary in Auckland, New Zealand. Ph.D. Thesis, University of Auckland, Auckland, New Zealand, 2005; p. 361. Available online: http://hdl.handle.net/2292/1419 (accessed on 2 April 2020).
  28. Shakeri, A.; Moore, F.; Modabberi, S. Heavy metal contamination and distribution in the Shiraz Industrial Complex Zone soil, South Shiraz, Iran. World Appl. Sci. J. 2009, 6, 413–425. [Google Scholar]
  29. Tang, R.; Ma, K.; Zhang, Y.; Mao, Q. The spatial characteristics and pollution levels of metals in urban street dust of Beijing, China. Appl. Geochem. 2013, 35, 88–98. [Google Scholar] [CrossRef]
  30. Paramasivam, K.; Ramasamy, V.; Suresh, G. Impact of sediment characteristics on the heavy metal concentration and their ecological risk level of surface sediments of Vaigai river, Tamilnadu, India. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2015, 137, 397–407. [Google Scholar] [CrossRef]
  31. Shakeri, A.; Shakeri, R.; Mehrabi, B. Contamination, toxicity and risk assessment of heavy metals and metalloids in sediments of Shahid Rajaie Dam, Sefidrood and Shirinrood Rivers, Iran. Environ. Earth Sci. 2016, 75, 679. [Google Scholar] [CrossRef]
  32. Luo, W.; Lu, Y.; Giesy, J.P.; Wang, T.; Shi, Y.; Wang, G.; Xing, Y. Effects of land use on concentrations of metals in surface soils and ecological risk around Guanting Reservoir, China. Environ. Geochem. Health 2007, 29, 459–471. [Google Scholar] [CrossRef]
  33. Islam, S.; Ahmed, K.; Mamun, H.; Masunaga, S. Potential ecological risk of hazardous elements in different land-use urban soils of Bangladesh. Sci. Total Environ. 2015, 512–513, 94–102. [Google Scholar] [CrossRef] [PubMed]
  34. Zhao, N.; Lu, X.W.; Chao, S.G. Level and contamination assessment of environmentally sensitive elements in smaller than 100μm street dust particles from Xining, China. Int. J. Environ. Res. Public Health 2014, 11, 2536–2549. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Bourennane, H.; Douay, F.; Sterckeman, T.; Villanneau, E.; Ciesielski, H.; King, D.; Baize, D. Mapping of anthropogenic trace elements inputs in agricultural topsoil from Northern France using enrichment factors. Geoderma 2010, 157, 165–174. [Google Scholar] [CrossRef]
  36. Rodríguez-Seijo, A.; Andrade, M.L.; Vega, F.A. Origin and spatial distribution of metals in urban soils. J. Soils Sedim. 2017, 17, 1514. [Google Scholar] [CrossRef]
  37. Hernandez, L.; Probst, A.; Probst, J.L.; Ulrich, E. Heavy metal distribution in some French forest soil, evidence for atmospheric contamination. Sci. Total Environ. 2003, 312, 195–219. [Google Scholar] [CrossRef] [Green Version]
  38. Mirzaei Aminiyan, M.; Baalousha, M.; Mousavi, R.; Mirzaei Aminiyan, F.; Hosseini, H.; Heydariyan, A. The ecological risk, source identification, and pollution assessment of heavy metals in road dust: A case study in Rafsanjan, SE Iran. Environ. Sci. Pollut. Res. 2018, 25, 13382–13395. [Google Scholar] [CrossRef]
  39. CCME. Canadian Soil Quality Guidelines for the Protection of Environment and Human Health; Canadian Council of Ministers of the Environment: Winnipeg, MB, Canada, 2018. [Google Scholar]
  40. Esmaeli, A.; Moore, F.; Keshavarzi, B.; Jaafarzadeh, N.; Kermani, M. A geochemical survey of heavy metals in agricultural and background soils of the Isfahan industrial zone, Iran. Catena 2014, 121, 88–89. [Google Scholar] [CrossRef]
  41. Pourang, N.; Noori, A.S. Heavy Metals Contamination in Soil, Surface Water and Groundwater of an Agricultural Area Adjacent to Tehran oil Refinery, Iran. Int. J. Environ. Res. 2014, 8, 871–886. [Google Scholar]
  42. Kabata-Pendias, A. Trace Elements in Soils and Plants, 4th ed.; CRC Press: Boca Raton, FL, USA, 2010; pp. 41–42. [Google Scholar]
  43. Mohammadpour, G.A.; Karbassi, A.R.; Baghvand, A. Origin and spatial distribution of metals in agricultural soils. Global J. Environ. Sci. Manag. 2016, 2, 145–156. [Google Scholar]
  44. Keshavarzi, B.; Moore, F.; Ansari, M.; Rastegari Mehr, M.; Kaabi, H.; Kermani, M. Macronutrients and trace metals in soil and food crops of Isfahan Province, Iran. Environ. Monit. Assess. 2015, 187, 4113. [Google Scholar] [CrossRef]
  45. Rahmanipour, F.; Marziaoli, R.; Bahrami, H.A.; Fereidouni, Z.; Rahimi Bandarabadi, S. Assessment of soil quality indices in agricultural lands of Qazvin Province, Iran. Ecol. Indic. 2014, 40, 19–26. [Google Scholar] [CrossRef]
  46. Mirzaei, R.; Teymourzade, S.; Sakizadeh, M.; Ghorbani, H. Comparative study of heavy metals concentration in topsoil of urban green space and agricultural land uses. Environ. Monit. Assess. 2015, 187, 741. [Google Scholar] [CrossRef] [PubMed]
  47. Ahmadi, M.; Jorfi, S.; Azarmansuri, A.; Jaafarzadeha, N.; Hosein Mahvi, A.; Darvishi Cheshmeh Soltani, R.; Akbari, H.; Akhbarizadeh, R. Zoning of heavy metal concentrations including Cd, Pb and As in agricultural soils of Aghili plain, Khuzestan province, Iran. Data Brief 2017, 14, 20–27. [Google Scholar] [CrossRef] [PubMed]
  48. Li, Y.; Gou, X.; Wang, G.; Zhang, Q.; Su, Q.; Xiao, G. Heavy metal contamination and source in arid agricultural soil in central Gansu Province, China. J. Environ. Sci. 2008, 20, 607–612. [Google Scholar] [CrossRef]
  49. Kouchou, A.; El Ghachtouli, N.; Duplay, J.; Ghazi, M.; Elsass, F.; Thoisy, J.C.; Bellarbi, M.; Ijjaali, M.; Rais, N. Evaluation of the environmental and human health risk related to metallic contamination in agricultural soils in the Mediterranean semi-arid area (Saiss plain, Morocco). Environ. Earth Sci. 2020, 79, 131. [Google Scholar] [CrossRef]
  50. Carkovic, A.B.; Calcagni, M.S.; Vega, A.S.; Coquery, M.; Moya, P.M.; Bonilla, C.A.; Pastén, P.A. Active and legacy mining in an arid urban environment: Challenges and perspectives for Copiapó, Northern Chile. Environ. Geochem. Health 2016, 38, 1001–1014. [Google Scholar] [CrossRef]
  51. Yang, P.; Mao, R.; Shao, H.; Gao, Y. An investigation on the distribution of eight hazardous heavy metals in the suburban farmland of China. J. Hazard Mater. 2009, 167, 1246–1251. [Google Scholar] [CrossRef] [PubMed]
  52. Micó, C.; Peris, M.; Sánchez, J.; Recatalá, L. Heavy metal content of agricultural soils in a Mediterranean semiarid area: The Segura River Valley (Alicante, Spain). Span. J. Agric. Res. 2006, 4, 363–372. [Google Scholar] [CrossRef] [Green Version]
  53. 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]
  54. 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. 2014, 48, 616–626. [Google Scholar] [CrossRef]
  55. Qing, X.; Yutong, Z.; Shenggao, L. Assessment of heavy metal pollution and human health risk in urban soils of steel industrial city (Anshan), Liaoning, Northeast Chin. Ecotoxicol. Environ. Saf. 2015, 120, 377–385. [Google Scholar] [CrossRef] [PubMed]
  56. Manta, D.D.; 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]
  57. Acosta, J.A.; Faz, A.; Martinez-Martinez, S. Identification of heavy metal sources by multivariable analysis in a typical Mediterranean city (SE Spain). Environ. Monit. Assess. 2010, 169, 519–530. [Google Scholar] [CrossRef] [PubMed]
  58. Lu, X.; Wang, L.; Li, L.-.Y.; Lei, K.; Huang, L.; Kang, D. Multivariate statistical analysis of heavy metals in street dust of Baoji, NW China. J. Harzard. Mat. 2010, 173, 744–749. [Google Scholar] [CrossRef] [PubMed]
  59. Micó, C.; Recatalá, L.; Peris, M.; Sánchez, J. Assessing heavy metal sources in agricultural soils of an European Mediterranean area by multivariate analysis. Chemosphere 2006, 65, 863–872. [Google Scholar] [CrossRef]
  60. Cheraghi, M.; Lorestani, B.; Merrikhpour, H. Investigation of the Effects of Phosphate Fertilizer Application on the Heavy Metal Content in Agricultural Soils with Different Cultivation Patterns. Biol. Trace Elem. Res. 2012, 145, 87. [Google Scholar] [CrossRef]
  61. Lugon-Moulin, N.; Ryan, L.; Donini, P.; Rossi, L. Cadmium content of phosphate fertilizers used for tobacco. Agron Sustain. Dev. 2006, 26, 151–155. [Google Scholar] [CrossRef]
Figure 1. Study area and location of the sampling stations (orange dots).
Figure 1. Study area and location of the sampling stations (orange dots).
Agronomy 10 00564 g001
Figure 2. Box plot of the enrichment factor for studied PTE in agricultural soils.
Figure 2. Box plot of the enrichment factor for studied PTE in agricultural soils.
Agronomy 10 00564 g002
Figure 3. Box plot of the contamination factor for studied PTE in agricultural soils.
Figure 3. Box plot of the contamination factor for studied PTE in agricultural soils.
Agronomy 10 00564 g003
Table 1. Summary statistics of physical and chemical properties and potentially toxic elements (PTE) levels (mg·kg−1) for the agricultural soil samples of Bandar Abbas County South District (BAC).
Table 1. Summary statistics of physical and chemical properties and potentially toxic elements (PTE) levels (mg·kg−1) for the agricultural soil samples of Bandar Abbas County South District (BAC).
ParameterUnitsMinimumMaximumAverageStandard DeviationSkewnessKurtosisCoefficient of VariationIEQSCSQGWMSBackground
pH-7.38.457.730.330.61−0.344.25 --
Organic Carbon%Bdl1.130.420.350.85−0.0382.51 --
CECmeq 100g−15.415.638.732.710.850.670.29
Sand%148043.5318.790.27−0.4243.16 --
Silt 126842.0517.89−0.10−1.2142.54 --
Clay 82014.423.86−0.18−0.8626.79 --
Almg·kg−118,15636,70226,51654240.37−0.720.46 80,0002800
As 4.48.86.531.28−0.06−0.9419.58-124.77.65
Cd 0.070.340.150.071.673.0843-1.41.10.17
Co 9.225.5514.644.191.010.9928.63-406.910.64
Cr 61237126.6850.860.61−0.3740.15110644274.5
Cu 14.839.0524.397.020.55−0.6828.77200631430.89
Fe 21,79040,18030,74654630.06−0.6617.77--47,20040,857
Mn 509886694.71115.82−0.04−1.2516.67--418879.70
Mo 0.71.350.960.250.52−1.525.64-51.81.06
Ni 65299.5125.3257.071.653.7445.54110451877
P 41288629.43141.880.23−0.9122.54- --
Pb 6.213.58.821.291.3222.75-702510.8
Sc 6149.582.460.37−0.8925.65--9.510.08
V 58.511887.0318.370.13−0.9821.1-1306098.08
Zn 35.260.350.027.46−0.68−0.2614.925002506265.09
IEQS—Iranian Environmental Quality Standards for agricultural soils [13], CSQG Canadian Soil Quality Guidelines [39], WMS—World average values in soils [42], CEC Cation exchange capacity.
Table 2. Comparison of studied PTE in agricultural soils between the present work and other studies from arid and semiarid environments (mean values, mg·kg−1).
Table 2. Comparison of studied PTE in agricultural soils between the present work and other studies from arid and semiarid environments (mean values, mg·kg−1).
LocationSoil TypenAlAsCdCoCrCuFeMnMoNiPPbSbVZnReference
Bandar Abbas, Southern IranAGS1926,5166.530.1514.64126.6824.3930,746694.710.96125.32629.438.8-87.0350.02This study
Granada and Almeria provinces, South-East of SpainAGS142--0.4-29.625.7---26.913.625.6--65.7[3]
Kermanshah province, Western IranAGS53----79.2141.2125,936559.06-131.46----74.62[13]
Hamedan Province, Western IranAGS58------37,690403.37 -21.98-2.43--[43]
Teheran area, Northwestern IranAGS141----87.15--665.06-43.60-35.08-67.14130.36[41]
Industrial area Isfahan, Center IranAGS10553,000-0.4314.785.935.728,000649.9-66.2-34.6--111.5[40]
Isfahan Province, Center IranAGS114-----35.2835,600---375032.14--100.96[44]
Qazvin Province, North-West IranAGS71--0.1418.934.5------17.36---[45]
Semnan Province, Northern IranAGS47--0.3-83.724.3---33.3-18.1--80.4[46]
Aghili plain, Khuzestan province, Western IranAGS55-2.810.29--------6.12---[47]
Gansu Province, ChinaAGS5-8.8--40.117.1-----23.30---[48]
Gansu Province, ChinaAGS-i13-11.17--38.8227.2-----21.44---[48]
Saiss plain, Northern MoroccoAGS-i ----7855---31----119[49]
Copiapó, Northern ChileAGS7-32.1-13.3-20921,1801049---50.4--176[50]
Hebei Province, Eastern ChinaAGS100-6.160.15-57.7721.22---25.04-18.80--69.96[51]
Alicante province, Eastern SpainAGS29--0.387.928.321.615,274320-23.7-19.6--57.8[52]
Castellón province, Eastern SpainAGS77--0.327.733.336.616,915385 19.5 55.8 78.5[53]
n—number of samples, AGS—Agricultural Soil, AGS-I—Agricultural irrigation soil.
Table 3. Potential ecological risk and potential ecological risk index (PERI) of PTE in the agricultural soils.
Table 3. Potential ecological risk and potential ecological risk index (PERI) of PTE in the agricultural soils.
Degree of Ecological Risk of Each ElementPotential Ecological Risk Index
AsCdCoCrCuNiPbZn
Maximum13.47101.738.726.909.5820.459.601.21156.22
Mean10.0041.715.003.715.989.065.601.0082.05
Minimum6.7313.963.141.643.634.223.520.70156.22
Table 4. Pearson’s correlation analysis of PTE in agricultural soils.
Table 4. Pearson’s correlation analysis of PTE in agricultural soils.
ScCdCrCuFeNiZnAlCoMnMoVAs_PbCECpHOCSandSiltClay
Sc1
Cd0.651
Cr0.66−0.651
Cu0.88−0.560.511
Fe0.96−0.590.690.871
Ni0.86−0.580.790.760.831
Zn0.60NCNC0.640.660.461
Al0.87NCNC0.840.860.750.781
Co0.95−0.60.740.870.930.940.590.861
Mn0.67−0.460.610.760.760.55NC0.510.621
Mo−0.840.72−0.75−0.638−0.81−0.72NC−0.630.76−0.521
V0.95−0.680.730.820.960.770.5190.760.880.67−0.881
AsNCNCNCNCNCNCNCNCNCNCNCNC1
PbNC0.63−0.55NCNCNCNCNCNCNCNC−0.46NC1
CEC0.51−0.610.460.49NC0.48NCNC0.52NCNC0.50NCNC1
pHNCNC0.52NCNCNCNCNCNC0.62NC0.51NCNCNC1
OC−0.49NCNC-−0.58−0.48NCNC−0.53−0.47−0.55NCNCNCNCNCNC1
SandNCNCNCNCNCNCNC−0.58NCNCNCNCNCNCNCNCNC1
Silt0.48NCNCNCNCNCNC0.680.50NCNCNCNCNC0.48NCNC−0.981
Clay−0.61NC−0.58NC−0.68NCNCNC−0.48−0.740.60−0.67NC0.517NC−0.610.46NCNC1
NC No correlation; bold letter, stands for highly significant correlations (p < 0.01); Italic letter stands for significant correlations (p < 0.05).
Table 5. Principal component analysis for physicochemical properties and PTE in BAC soils.
Table 5. Principal component analysis for physicochemical properties and PTE in BAC soils.
Agricultural Soils
PC1PC2PC3
Fe0.96−0.04−0.22
Co0.94−0.200.01
Al0.940.100.24
Cu0.91−0.170.02
V0.90−0.17−0.29
Sc0.89−0.02−0.19
Mo−0.780.110.31
Zn0.720.430.14
Mn0.70−0.15−0.45
Cr0.640.47−0.33
OC−0.530.130.12
P−0.010.890.26
Cd−0.290.870.25
Ni0.020.65−0.02
Pb−0.080.710.49
CEC0.49−0.510.16
Clay−0.53−0.050.78
As0.190.230.76
pH0.31−0.34−0.66
% of Variance42.9219.4514.00
Cumulative %42.9262.3776.37
Bartlett’s sphericity tests (p = 0); Kaiser–Meyer–Olkin (KMO = 0.715).

Share and Cite

MDPI and ACS Style

Moghtaderi, T.; Shakeri, A.; Rodríguez-Seijo, A. Potentially Toxic Element Content in Arid Agricultural Soils in South Iran. Agronomy 2020, 10, 564. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10040564

AMA Style

Moghtaderi T, Shakeri A, Rodríguez-Seijo A. Potentially Toxic Element Content in Arid Agricultural Soils in South Iran. Agronomy. 2020; 10(4):564. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10040564

Chicago/Turabian Style

Moghtaderi, Tahereh, Ata Shakeri, and Andrés Rodríguez-Seijo. 2020. "Potentially Toxic Element Content in Arid Agricultural Soils in South Iran" Agronomy 10, no. 4: 564. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10040564

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