Original Articles
Geostatistical mapping and quantitative source apportionment of potentially toxic elements in top- and sub-soils: A case of suburban area in Beijing, China

https://doi.org/10.1016/j.ecolind.2020.106085Get rights and content

Highlights

  • Most V, Cr, Ni and As were endogenous, while Cd, Zn, Pb and Hg were anthropogenic and exogenous.

  • Urban emission was identified as a major source for 8 PTEs in suburban soils.

  • Relative contributions of three sources to 8 PTEs were estimated by PMF.

Abstract

The risk assessment and source identification for potentially toxic elements (PTEs) in soils, particularly agricultural soils from megacities, are significant for environmental protection and pollution control. In this study, an intensive sampling (4127 topsoil samples and 994 subsoil samples) was conducted in the Shunyi District, Beijing, which is a suburban area with extensive cropland cover and has been impacted by the megacity over several decades. Concentrations and distributions of 8 PTEs, including V, Cr, Ni, As, Cd, Zn, Pb and Hg, were determined, and their possible sources were quantitatively assessed by principal component analysis (PCA), redundancy analysis (RDA), positive matrix factorization (PMF) analysis, and anthropogenic contribution ratio method. Among 8 PTEs, Zn, V and Cr exhibited significantly high concentrations in soils, with means of 68.29, 68.19 and 52.13 mg/kg, respectively, followed by Pb (23.84 mg/kg), Ni (22.91 mg/kg), As (8.30 mg/kg), Cd (0.15 mg/kg) and Hg (0.05 mg/kg). RDA and PCA demonstrated that the rock weathering was a significant source of V, Cr, Ni and As, and the local emissions and atmospheric deposition respectively contributed most of Cd, Zn and Pb, and of Hg in soils. This source category was confirmed the spatial variations of anthropogenic contribution ratios to individual PTEs. PMF results showed that the local emissions contributed 96.3% of Cd, 44.4% of Zn and 32.0% of Pb in soils, and the atmospheric source carrying urban pollutants amounted to 78.7–80.2% of Hg. In this case, several effective analysis methods have been successfully applied to quantify the impact of a megacity to PTEs in suburban soils. These results improve understanding of the contamination status of PTEs in suburban soils from Beijing megacity, and provide basis for policymaker regarding environmental protection and pollution control.

Introduction

Soil, as an important component of terrestrial ecosystems, plays a crucial role in biochemical transformation, the cycling of elements, and the filtration of water (Luo et al., 2012). It also serves as a reservoir for potentially toxic elements (PTEs), which could be transported from atmosphere, biomass and hydrosphere (Lv et al., 2014). Because of the environmental persistence, high toxicity and great bioaccumulation, the presence of PTEs in soils may threaten ecological safety and endanger the health of organisms through food chain (Akinwunmi et al., 2017, Khan et al., 2008, Lin et al., 2017, Tóth et al., 2016). In addition, PTE exposure could directly reduce immune function through suppressing immune system, as well as through indirect effects to the microbiota, both of which would increase risks from potential diseases (Keesing et al., 2010). Therefore, global concern has been raised over PTE contamination in soils (Chen et al., 2016a, Cortada et al., 2018, Gbadamosi et al., 2018, Ravankhah et al., 2017). Studies focusing on the occurrence and pattern of PTEs in soils would improve understandings for their environmental loads and prediction of human exposure risks.

Both natural and anthropogenic sources could contribute PTEs to soil environment (Yuan et al., 2013). Natural source of PTEs is generally a result of parent rock weathering or pedogenesis (Lin et al., 2017, Wang et al., 2019a), and anthropogenic pollutants are directly emitted from human activities (Wang et al., 2019b, Xiao et al., 2019, Zhao et al., 2014). In developing countries, the rapid economic development sharply increased anthropogenic emissions of PTEs, which would be transported atmospherically and deposited into nearby soils (Jiang et al., 2017). This may result in soil contamination at regional and even global scales. In China, a series of reports found that soils were obviously contaminated by PTEs, which increased from 185 mg/kg in 1990 to 400 mg/kg in 2015 (Chen et al., 2015, Pan et al., 2018, Zhang et al., 2018b). Over the past decades, the fast urban expansion in China has dramatically increased the industrial and municipal wastewater discharges, from 17.6 to 32.1 billion tons and 35.1 to 65.9 billion tons, respectively (Oyang and Wang, 2000, Pan et al., 2018). These sources, together with the overwhelming traffic emissions, resulted in a substantial input of PTEs into urban soils, and therefore, much attention has been paid to the urban soil contamination in some Chinese megacities (Jin et al., 2019, Liu et al., 2019, Lv and Liu, 2019). Beijing, the capital of China with a population over 20 million, is experiencing high rate of urbanization and associated soil contamination in urban areas (Chen et al., 2010, Liu et al., 2016, Wang et al., 2012), in which the PTE concentrations were detected as high at ~300 mg/kg (Chen et al., 2015). A recent study in Beijing by Yu et al. (2019) noted that a number of PTEs, including Hg, Cd, Pb, As, Zn, Cr, Ni, Co, and V would be dispersed outside with soil particles, and even be transported from urban areas to remote areas along a certain pathway. This finding may indicate that the particle-borne PTEs would be deposited into suburban soils with various distances from cities. However, few attention has been paid to PTEs in suburban soils, which may serve as a sink area for the urban emissions (Lin et al., 2017). This presents a knowledge gap regarding the impact of urban emissions on PTEs in suburban soils.

Because the human activities in suburban areas would also release PTEs, the levels and patterns of PTEs in suburban soils may not correspond to the contribution from urban transport. In Beijing, the Shunyi District is a typical suburban area. Although soils from the Shunyi District were reportedly contaminated by Beijing urban emissions (Wang et al., 2018), applications of chemical fertilizer, sewage sludge, and livestock manure for agricultural purpose would increase the soil PTE levels in this area (Luo et al., 2009, Xu and Zhang, 2017, Yang et al., 2013). In addition to the anthropogenic activities, the soil parent materials would significantly affect PTE concentrations in suburban soils as well (Guagliardi et al., 2018). However, the relative importance of rock weathering or pedogenesis to PTEs was unclear or previously neglected (Lin et al., 2017, Wang et al., 2019a).

To evaluate the various contributions to PTEs in soils, the geostatistical mapping has been used to identify the spatial features and possible sources of PTEs in soils (Li et al., 2014, Li and Feng, 2012, Men et al., 2018, Xie et al., 2011). Recently, a newly-developed method combining geostatistical analyses, multivariate analyses, redundancy analysis and robust geostatistics, has been successfully applied in discriminating the soil PTEs inherited from rock weathering from those derived by anthropogenic sources (Lin et al., 2017, Wang et al., 2019a). This method may provide knowledge of source apportionment of PTEs in soils from the Shunyi District, and further qualify the relative contributions from different sources. It was also interesting to compare this method with a traditional source appointment method, such as positive matrix factorization (PMF) model (Lv, 2019, Lv and Liu, 2019, Zhang et al., 2018c). In this study, the goals are: (1) to investigate concentrations and spatial distributions of PTEs in topsoils and subsoils from the Shunyi district, Beijing; (2) to discriminate anthropogenic contributions from natural source of rock weathering; and (3) to apportion quantitatively the urban emission for PTEs in suburban soils.

Section snippets

Study area and sampling

Our study area covers the whole area of the Shunyi District (40°00′~40°18′N, 116°28′~116°58′E, 1020 km2) located in the northeast of Beijing. The main topography of this district is plain, with mild terrain slopes (0–3°) from north to south. The soils in the study area are comprised mainly by moisture soil (64.1%) and gray cinnamon soil (23.7%), both of which varies widely in soil textures (Fig. 1), ranging from sand soil (7.1%) and silty sand soil (50.3%) to clay soil (1.8%). The bedrock are

Concentrations of PTEs in soils

The descriptive statistics of 8 PTE concentrations in top- and subsoil samples are summarized in Table S2. In topsoils, the mean concentrations of V, Cr, Ni, As, Cd, Zn, Pb and Hg were 68.19, 52.13, 22.91, 8.30, 0.15, 68.29, 23.84 and 0.05 mg/kg, respectively. These values were observed at similar levels compared with the background concentrations of PTEs in soils in Beijing (Table S2). In subsoils, the mean concentrations of V, Cr, Ni, As, Cd, Zn, Pb and Hg were 67.89, 52.09, 22.85, 8.25,

Source interpretation

V, Cr, Ni and As were suggested to be derived from the same source and tightly associated with Fe2O3, Al2O3 through PCA/RDA analyses. As Al2O3, SiO2 and Fe2O3 are the primary products of rock weathering, the levels of Al2O3/SiO2 ratios and Fe2O3 concentrations are expected to reflect the degree of rock weathering or pedogenesis (Lin et al., 2017, Wang et al., 2019a). As shown in the distribution maps (Fig. 3), the high contents of soil V, Cr, Ni and As were overlapped with those of Al2O3/SiO2

Conclusion

A total of 4127 top- and 994 sub-soil samples were collected from the Shunyi District, Beijing. The mean concentrations of V, Cr, Ni, As, Cd, Zn, Pb and Hg in topsoils were 68.19, 52.13, 22.91, 8.30, 0.15, 68.29, 23.84 and 0.05 mg/kg, respectively. PCA and RDA analyses demonstrated that the rock weathering was a main source of V-Cr-Ni-As in topsoils, and local emissions and atmospheric deposition were the dominant contributors for Cd-Zn-Pb and Hg, respectively. The source-PTE classification in

CRediT authorship contribution statement

Xu-Chuan Duan: Software, Data curation, Writing - original draft. Hong-Hui Yu: Formal analysis. Tian-Rui Ye: Investigation, Resources. Yong Huang: Investigation. Jun Li: Investigation, Visualization. Guo-Li Yuan: Supervision, Conceptualization, Methodology, Writing - review & editing. Stefano Albanese: Software, Validation.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This study was financially supported by the National Natural Science Foundation of China (41872100, 41372249) and the Fundamental Research Funds for the Central Universities (2652018160, 2652018158). We sincerely thank the project members from Beijing Institute of Geo-exploration Technology for their help in collecting the samples.

References (68)

  • X. Fei et al.

    Comprehensive assessment and source apportionment of heavy metals in Shanghai agricultural soils with different fertility levels

    Ecol. Indic.

    (2019)
  • L. Gao et al.

    Aquatic environmental changes and anthropogenic activities reflected by the sedimentary records of the Shima River, Southern China

    Environ. Pollut.

    (2017)
  • M. Gabarrón et al.

    Use of multivariable and redundancy analysis to assess the behavior of metals and arsenic in urban soil and road dust affected by metallic mining as a base for risk assessment

    J. Environ. Manage.

    (2018)
  • M.R. Gbadamosi et al.

    Distribution of radionuclides and heavy metals in the bituminous sand deposit in Ogun State, Nigeria–a multi-dimensional pollution, health and radiological risk assessment

    J. Geochem. Explor.

    (2018)
  • I. Guagliardi et al.

    Geochemical sources of vanadium in soils: evidences in a southern Italy area

    J. Geochem. Explor.

    (2018)
  • Q. Guan et al.

    Source apportionment of heavy metals in agricultural soil based on PMF: a case study in Hexi Corridor, northwest China

    Chemosphere

    (2018)
  • L. Hakanson

    An ecological risk index for aquatic pollution control. A sedimentological approach

    Water Res.

    (1980)
  • Y. Jiang et al.

    Source apportionment and health risk assessment of heavy metals in soil for a township in Jiangsu Province, China

    Chemosphere

    (2017)
  • Y. Jin et al.

    Assessment of sources of heavy metals in soil and dust at children's playgrounds in Beijing using GIS and multivariate statistical analysis

    Environ. Int.

    (2019)
  • S. Khan et al.

    Health risks of heavy metals in contaminated soils and food crops irrigated with wastewater in Beijing, China

    Environ. Pollut.

    (2008)
  • P. Lazo et al.

    Spatial distribution and temporal trend of airborne trace metal deposition in Albania studied by moss biomonitoring

    Ecol. Indic.

    (2019)
  • W. Li et al.

    The identification of ‘hotspots’ of heavy metal pollution in soil-rice systems at a regional scale in eastern China

    Sci. Total Environ.

    (2014)
  • X. Li et al.

    Multivariate and geostatistical analyzes of metals in urban soil of Weinan industrial areas, Northwest of China

    Atmos. Environ.

    (2012)
  • X. Li et al.

    Heavy metal contamination of urban soil in an old industrial city (Shenyang) in Northeast China

    Geoderma

    (2013)
  • J. Liang et al.

    Spatial distribution and source identification of heavy metals in surface soils in a typical coal mine city, Lianyuan, China

    Environ. Pollut.

    (2017)
  • Y. Lin et al.

    Source identification of potentially hazardous elements and their relationships with soil properties in agricultural soil of the Pinggu district of Beijing, China: multivariate statistical analysis and redundancy analysis

    J. Geochem. Explor.

    (2017)
  • E. Liu et al.

    Chemical speciation, pollution and ecological risk of toxic metals in readily washed off road dust in a megacity (Nanjing), China

    Ecotoxicol. Environ. Saf.

    (2019)
  • J. Liu et al.

    Emission control priority of PM2.5-bound heavy metals in different seasons: a comprehensive analysis from health risk perspective

    Sci. Total Environ.

    (2018)
  • R. Liu et al.

    Spatial pattern of heavy metals accumulation risk in urban soils of Beijing and its influencing factors

    Environ. Pollut.

    (2016)
  • L. Luo et al.

    An inventory of trace element inputs to agricultural soils in China

    J. Environ. Manage.

    (2009)
  • X.S. Luo et al.

    Trace metal contamination in urban soils of China

    Sci. Total Environ.

    (2012)
  • J. Lv

    Multivariate receptor models and robust geostatistics to estimate source apportionment of heavy metals in soils

    Environ. Pollut.

    (2019)
  • J. Lv et al.

    An integrated approach to identify quantitative sources and hazardous areas of heavy metals in soils

    Sci. Total Environ.

    (2019)
  • J. Lv et al.

    Multivariate geostatistical analyses of heavy metals in soils: Spatial multi-scale variations in Wulian, Eastern China

    Ecotoxicol. Environ. Saf.

    (2014)
  • Cited by (39)

    View all citing articles on Scopus
    View full text