Influence of the soil sealing on the geoaccumulation index of heavy metals and various pollution factors

Soil sealing belongs to the most destructive and damaging processes to the soil environment. Soil sealing interrupts or greatly restricts the exchange of matter and energy between the biosphere, hydrosphere, and atmosphere and the soil environment. The aim of this study was to compare the content of heavy metals (Cd, Cr, Cu, Hg, Fe, Ni, Pb, Zn) of Ekranic Technosols by applying indicators such as geoaccumulation index (I geo), enrichment factor (EF), and pollution load index (PLI), which allowed to determine quantitatively the impact of the soil sealing degree on the content of heavy metals and to distinguish natural from anthropogenic sources of origin of heavy metals. In general, 42 soils from different parts of the city of Toruń (NW Poland) were sampled and divided into three groups according to the degree of soil sealing: completely sealed with asphalt or concrete (A), semi-permeable (partially sealed with cobblestones and concrete paving slabs (B)), and reference (non-sealed) (C). The results indicate that the artificial sealing in urban areas slightly affects the content of heavy metals in soils. However, based on PLI, I geo, and EF, it was found that the sealing has influence on soil properties and unsealed soil is the most exposed to the accumulation of pollutants.


Introduction
Understanding the complexity of the functioning of soil systems and the interaction with human activity is particularly important in urban areas. Due to a continuous human existence and activities, these soils are very often truncated or buried under transported material, as well as highly compacted and contaminated. As a consequence, their hydrological characteristics, degree of compaction and disturbances, and relocations of material from original horizons are extremely different from those of natural soil types.
The concentration of heavy metals in soils, associated with litogenesis and pedogenesis, depends on the mineralogical composition of parent material and the direction and pace of the process of soil formation, which determines the distribution of trace elements in the soil profile (Luo et al. 2012a;Wong et al. 2006). Moreover, technogenic activity applies all of the abovementioned processes. The study of heavy metal deposition and accumulation is of increasing interest because of the awareness that heavy metals present in soil may have negative consequences on the human health and whole environment. Heavy metals may enter into aquatic ecosystems from anthropogenic sources, such as industrial waste water discharges, sewage waste water, fossil fuel combustion, and atmospheric deposition. Urban soils are changed as a result of activity of various processes, i.e., sealing, compaction, storage, and mixing (Plak et al. 2015;Wei et al. 2013). Sealing the soil belongs to the most Responsible editor: Zhihong Xu destructive processes to the soil environment. In principle, it is an irreversible process, which is defined as the destruction of soil cover by partial or total application of the soil impermeable layer. Total sealing of the soil in urban areas is caused by concrete or asphalt (Siebielec et al. 2015). Soils can also be sealed with a semi-permeable surface, e.g., concrete paving, which allows partial penetration of air and water (Nestroy 2006) and allows saving the certain features of soil (Piotrowska-Długosz and Charzyński 2015). In the European Union, on average 51% of artificial surfaces are sealed, but the composition varies considerably among Member States, depending on the dominant settlement structures and the intensity of the interpretation of artificial surfaces (Siebielec et al. 2015). In urban areas, as a direct result of human activity, the impact of industrial and municipal building, industry and communication, and artifacts deposited in soil, heavy metals are present in high concentrations Madrid et al. 2004;Puskas and Farsang 2009;Xia et al. 2011).
In order to determine the degree of contamination of surface soil levels with heavy metals, their concentration is compared with the content of the local geochemical background, or by calculation of geoaccumulation index (I geo ), enrichment factor (EF), or pollution load index (PLI) (Horckmans et al. 2005;Muller 1969;Tomlinson et al. 1980). The main factors influencing the distribution of heavy metals in the soil profile are the content of humus and iron and manganese oxides, as well as texture, soil pH, and the rinsing processes, accompanied by the movement of water within the soil profile (Luo et al. 2012b;Madrid et al. 2008). Metals of anthropogenic origin are more mobile than those of pedogenic origin. In urban areas, the nature of the metals accumulation and their genesis in the soil is often ambiguous, and all the abovementioned processes, which are widely described in the literature, are distorted when the soil is sealed (Kabata-Pendias and Pendias 1999;Reimann et al. 2005;Sauerwein 2011;Siebielec et al. 2015).
The aim of this study was to compare the content of heavy metals (Cd, Cr, Cu, Hg, Ni, Pb, Zn) of Ekranic Technosols from Toruń, Poland, and adjacent, non-sealed soils. We assumed that the soils sealed for 30 to 40 years affect the diversity of heavy metal content in two separate groups of the completely and partially sealed soils while compared with the non-sealed soil. Moreover, we assumed that the use of indicators such as I geo , EF, and PLI allows to determine quantitatively the impact of the soil sealing on the content of heavy metals.
The problem of uncontrolled propagation of the sealed surface in Poland, especially in the cities, is remotely explored, and the processes occurring in the sealed soils are not comprehensively studied. Moreover, few researchers in the world do their research on geochemical environment of the sealed soil, focusing primarily on the changes in the physical properties or biological activity (lithium) (Greinert 2015;Puskás and Farsang 2009;Zhao et al. 2012). The presented studies show the unique issues related to the impact of (completely or partially) sealed soil on disturbance of heavy metal circulation, which is of great importance in terms of quality of life and health of urban residents, the condition of the environment in cities, and the penetration of pollutants into groundwater. The research material gives extensive knowledge about the impact of the sealing process on soil geochemistry.

Materials and methods
Description of the study area and soil sampling The study sites were located within the area of Toruń (18.609°E, 53.020°N), the city in the North West Poland (Fig. 1). Majority of natural and technogenically transformed soils in the studied area have texture of medium sand according to USDA Soil Taxonomy (Soil Survey Staff 2014), but some have texture of coarse and fine sand. The moderate climate of the region is between the marine type of Western Europe and the continental type of Eastern Europe. The average annual temperature for the whole year was 8.4°C, while the annual rainfall was 535 mm (Source: data of the Institute of Meteorology and Water Management: stat.gov.pl).
In general, 42 soils from various parts of the city (Table 1) were sampled between April of 2012 and September 2013 ( Fig. 1). In each place, soils sealed with asphalt or concrete, fulfilling criteria of technic hard material (IUSS Working Group WRB 2015) (completely sealed, group A; 15 samples), semi-permeable materials (cobblestones and concrete paving slabs, partially sealed, group B; 11 samples), and reference (non-sealed, group C; 16 samples) soils were selected (Fig. 2). They were situated in direct proximity to ensure similar history of technogenic disturbances, the same environmental conditions, and soil texture. Samples from both categories of sealed soils were collected from the topmost original horizon that was left after the process of construction of sealing surface. Soil samples were collected from a depth of 15-25 cm or 10-20 cm. Depth of sampling depended on the thickness of the technic hard material combined with the ballast layer used for stabilization of the pavement construction. Reference sites without sealing were located on lawns of green-belts next to pavements or roads (about 1 m) to ensure maximum soil comparability. Samples were collected from the same depth as in the sealed soils. The sealing was performed between 30 and 40 years ago, in 1970 or 1980.
Chemical and physical analysis of the soil After air-drying at room temperature and sieving (<2 mm), the soil samples were analyzed for physical and chemical properties. Physicochemical properties were determined using standard methods (Van Reeuwijk 2002). Each soil sample was analyzed in triplicate. Texture was assessed by the areometric method combined with the sieve method. The pH in 1 mol l −1 KCl was measured using the potentiometric method in 1:2.5 soil:solution suspensions; total organic carbon (OC TOT ) and total nitrogen (N TOT ) content were determined using a dry combustion CN analyzer (Vario Max CN). Available phosphorus (P) was assayed using the vanadium-molybdenum method.

Heavy metal analysis
The total content of heavy metals (denoted by TOT ) in the soil samples was determined using the ICP MS technique before dissolving with aqua regia.

Data analysis
I geo (Muller 1969) is computed using the following equation: where C n is the measured concentration of the element in environment and B n is the geochemical background value in soil. According to Muller (1969), the I geo for each metal is calculated and classified as 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).
The EF calculation is expressed below as where C x is the concentration of the element of interest and C ref is the concentration of reference element for normalization. Enrichment factors of heavy metals were calculated for each soil sample relative to the background values of abundance of chemical elements in the local parent rock, choosing Fe as the reference element (according to Kabata-Pendias and Pendias (1999) −12.9 g kg −1 Fe).
Geochemical background was determined on the basis of Czarnowska (1996), Kabata-Pendias and Pendias (1999), Lis and Pasieczna (1995) and Pasieczna (2003). There have been analyzed uncultivated areas, taking into account the different geological provinces according to classification proposals used for geochemical map of the world (Damley 1995). The area of Toruń is situated in the northern province of the Polish Lowland, where sand is mostly present. Five contamination categories are recognized on the basis of the enrichment factor: EF < 2 states deficiency to minimal enrichment, EF = 2-5-moderate enrichment, EF = 5-20significant enrichment, EF = 20-40-very high enrichment, and EF > 40 extremely high enrichment (Duzgoren-Aydin et al. 2006;Sezgin et al. 2003).
PLI (Tomlinson et al. 1980) is expressed as follows: where CF is the contamination factor obtained by calculating between each metal's concentration and its background value. When PLI is greater than 1, it means that contamination exists; however, if PLI is less than 1, there is no metal contamination.

Statistical data evaluation
Data were analyzed using non-parametric tests (Statistica 6.1, Stat Soft. Inc.). Moreover, principal component analysis (PCA) of two factors and cluster analysis (CA) using the Euclidean distance were applied.

Results and discussion
The basic properties of soils The physical and chemical properties of the studied soils are presented in Table 2. The studied soils represented poorly diversified texture distribution, and according to the USDA, they represented mainly sand and loamy sand. In the studied profiles, there was a significant presence of the skeleton, which was mostly the building rubble. There is no difference between the soil groups in the average pH values (Table 2). Soil sealing significantly reduced the content of OC TOT and N TOT compared to non-sealed soil. The C:N ratio ranged widely from 9.3 to 28.9 and was not significantly different for the reference and sealed soils, or for both categories of soil sealing ( Table 2). Distribution of organic carbon content showed high variability, which is quite typical for the anthropogenic soil (e.g., replenishing a new layer or mixing of soil material) (Greinert 2015). Higher levels of OC TOT and N TOT in non-sealed soils in relation to other comparable groups can be primarily explained by the lack of barriers in the supply of organic matter but also higher biological activity in the non-sealed soil. The content of available phosphorus, P AV , also showed large range reaching from 47 to 456 mg kg −1 , but no clear trend appears between the compared groups of soil.

Relation of soil sealing and heavy metal content
Heavy metal concentrations in urban soils in Toruń are given in Table 3. The total content of heavy metals (Cd, Cr, Cu, Ni, Pb, Zn) in soils of all groups does not exceed the size of the content permitted in Poland (Regulation of the Minister of Environment on soil quality standards and earth quality standards, Dz.U.02.165.1359 dated 04.10.2002). The results were also compared with the permissible content of heavy metals in the soil in the Netherlands and Germany, and no exceeded values were detected (NGG 1995, BBodSchV, 1999. Analysis of heavy metal content in the three soil groups studied in Toruń arranged soil in the following order: group C (non-sealed) > group A (completely sealed) > group B (partially sealed) ( Table 3).
The level of geochemical background of heavy metals in relation to the completely sealed soil exhibits higher values for Cu, Pb, Zn, and Hg (Table 3). The highest concentration of heavy metals in all examined  Comparison of heavy metal content in sealed soils with other cities in Poland, Germany, Hungary, and Russia shows that the concentration of the analyzed elements in the urban soils of Toruń is on a medium level (Table 4).
The I geo index is used to assess heavy metal contamination in urban soils by comparing current and preindustrial concentrations, although it is not always easy to reach pre-industrial sediment layers. It is also employed in pollution assessment of heavy metals in urban road dust.
In this study, the background geochemical composition of the city soil types (Czarnowska 1996;Kabata-Pendias and Pendias 1999;Lis and Pasieczna 1995;Pasieczna 2003) is chosen as the local background value for calculating the I geo values, which are presented in Table 5. The constant equal 1.5 allows to analyze natural fluctuations in the content of a given substance in the environment and to detect very small anthropogenic influences.
The average value of I geo decreases for soils from group A in the following order: Pb > Hg > Cu > Cd = Zn > Ni > Cr, for soils from group B: Pb = Cd > Hg = Zn > Cu > Ni > Cr, a n d f o r s o i l s f r o m g r o u p C a s f o l l o w s : Pb > Hg > Zn > Cu > Cd > Ni > Cr. Despite negative values of I geo for Cr, Ni, Zn, and Cd in all types of studied soils, they should not be simply classified as Buncontaminated^because such values may be connected with low background concentration of Vistula valley, but not with lack of technogenic accumulation of pollutants. For Pb and Hg I geo values in reference soil and completely sealed soil belonged to the category of Buncontaminated to moderately contaminated,^while the average value of I geo for partially sealed soil was indicated as Buncontaminated.^This means that the results obtained for each of the elements are considerably different. In the case of  (Czarnowska 1996;Kabata-Pendias and Pendias 1999;Lis and Pasieczna 1995;Pasieczna 2003)  Cu, I geo in the group of sealed soils had a value higher than zero (Buncontaminated to moderately contaminated^). Generally, low average values of I geo in all examined groups of soil in Toruń have indicated a small level of soil contamination with heavy metals. Simultaneously, it is worth noting that on the basis of the maximum value of I geo for Cu, Cd, Pb, Hg, and Zn, sealed soil is classified as Bmoderately to heavily contaminated.M oreover, the contamination of Pb and Zn are associated (Fig. 4). It probably means that they have the same source, i.e., cars, but accumulation occurred at different chronological periods, which could be associated with the motor fuel quality (unleaded petrol). Comparison of average I geo value among the three groups of soil indicates a following order: partially sealed soil (group B) < non-sealed soil (group C) < completely sealed soil (group A).

Environmental risk assessment of studied soils
A great number of environmental risk assessment is known in the literature (Hilton et al. 1985;Tomlinson et al. 1980;Verca and Dolence 2005); therefore, the character of the analyzed area was taken into account, and two factors, EF and PLI, were chosen for the evaluation of the soil sealing.
The EF factor evaluated heavy metal pollution according to the content of heavy metals, but it hardly distinguished their source, chemical activity, or biological availability. The main advantage of EF application is the possibility to compare the results with the data presented by other researchers.
Calculation of EF was based on the standardization of a measured element against a reference element. A reference element is often the one characterized by low occurrence variability, such as Al, Fe, Ti, Si, Sr, and K (Duzgoren-Aydin 2007; Sezgin et al. 2003). EF values less than 5.0 are not considered significant, because such small enrichments may arise from differences in the composition of local soil material and reference soil used in EF calculations (Sezgin et al. 2003). However, there is no accepted pollution ranking system or categorization of degree of pollution on the enrichment ratio and/or factor methodology.
The calculated values of EF are presented in Table 6. In the completely sealed soil, average EF values decreased in the following order: Pb > Hg > Cu > Cd > Ni > Zn > Cr, while in the partially sealed soil, as follows: Pb > Cd > Hg > Zn > Cu > Ni > Cr. The average value of the EF for nonsealed soil showed the following diminishing order: Pb > Hg > Zn > Cu > Cd > Ni > Cr, similar to the order described for I geo . The analysis of EF values have indicated that the content of chromium in all examined groups of soil states Bdeficiency to minimal enrichment,^whereas the average values of EF for Zn, Ni and Cd have been in the range from 2 to 5, what is called the level of Bmoderate enrichment.^In the completely sealed soil, the average value of the EF for Hg, Cu, and Pb has been greater than 5, which clearly shows that these elements are derived from anthropogenic sources. In partially sealed soil, only Pb has  Gr. A-completely sealed soil, Gr. B-partially sealed soil, Gr. C-non-sealed soil  shown anthropogenic origin, and in the reference soil, average EF index above 5 has been observed for Pb and Cu. It seems, therefore, that EF can also be an effective indicator used to distinguish the natural origin of anthropogenic sources of heavy metals.
The PLI index was evaluated to assess the mutual contamination effects of the seven metals measured in this study.
Mean of PLI (Table 6) was the highest in the non-sealed soils and arranged in the decreasing order as follows: group C > group A > group B, which showed that the non-sealed soils occurred the most contaminated ones. Furthermore, the largest concentrations of pollutants such as polycyclic aromatic hydrocarbons (PAHs) and the greatest biological activity were detected in non-sealed soils in Toruń (Mendyk and Charzyński 2016;Piotrowska-Długosz and Charzyński 2015). The sealing of the soil accompanied limitation of the migration of heavy metals in the soil profile due to the reduction of the supply of the surface levels, limitation of water movement, increase in the specific density of the soil, and the reduction of the degree of aeration. Hence, we can conclude that higher concentrations of Hg, Cu, and Pb have been derived from anthropogenic sources before sealing the soil. It was observed that the content of heavy metals in all soil groups has often been significantly different. This could depend mainly on the type of substrate (e.g., debris, garbage, waste) and the origin of the pollutants such as sewage, sludge, and road traffic, which are characteristics for urban soils (Luo et al. 2012b). Similar conclusions were presented by Bezuglova et al. (2016) who studied sealed and non-sealed soil in Rostov-on-Don. Total contamination factor Z c (Revich et al. 1985) was used to show that the sealing of urban soils reduces the vertical migration of heavy metals in the soil profile. The determination of the level of risk of soil contamination with heavy metals depends on their content in the parent rock and on the technogenic deposition in the surface layers. In soils of Rostov-on-Don developed on loess, the content of heavy metals in parent rock was higher than the geochemical background (Bezuglova et al. 2016).
Data were evaluated using classical statistical methods. They did not show a normal distribution according to the Shapiro-Wilk test; thus, they were analyzed using non-parametric Kruskal-Wallis test, but any statistically significant differences were detected. Moreover, PCA of two factors, CA analysis using the Euclidean distance (Lee et al. 2006;Zheng et al. 2008), and Spearman's correlation coefficients were applied.  In PCA analysis, all heavy metals were represented by the first two principal components, which accounted for 79.9% of the total variance. Nevertheless, factors 1 and 2, which explained 65.4 and 14.5% of total variance, respectively, did not differentiate types of soil on the basis of the content of heavy metals (Fig. 3). The result of CA analysis is illustrated in Fig. 4, on which two distinct clusters can be identified. Cluster I contained Cd, Cr, Cu Ni, and Hg, while the long distance between Cu and the other four heavy metals suggested that this cluster could be further divided into two sub-clusters. Cluster II contained Pb and Zn. Concentrations of Cu, Pb, and Hg were significantly higher than the background values of Poland, and concentrations of Cr, Cd, Ni, and Zn were comparable to the background values. Therefore, previously described EF and I geo indexes show that the distribution of Cu, Pb, and Hg in urban soils of Toruń was affected mostly by anthropogenic sources, while Cr, Cd, Ni, and Zn were mainly from natural sources.
Moreover, significant Spearman's correlations between heavy metal content and physical properties of the urban soil in Toruń (Table 7) were determined mainly in nonsealed soil and then for completely sealed soil. In partially sealed soil, the amount of correlations was the least. However, partial or complete sealing did not differentiate the soil in terms of heavy metal content, which confirmed the idea that the soil is heterogeneous under technopressure and remains in contrast to the soil formed under natural conditions. Similarly, the results analyzed in the population of non-sealed soil were the evidence of their high heterogeneity. Soil sealing interrupts the exchange between the soil system and other ecological compartments, including the biosphere, hydrosphere, and atmosphere, which affects processes in the water cycle, biogeochemical cycles, and energy transfers (Siebielec et al. 2015). However, in case of the impact of sealing of urban soils on the content of heavy metals, we proved only slight effect. The sources of Cd, Cu, Cr, Hg, Ni, Pb, and Zn in urban soils of Toruń are extremely heterogenous, and thus, situation is rather complicated. Anthropogenic sources such as vehicle exhaust, household waste, and construction activities have made the heavy metal concentrations higher than their background values. Besides, as Toruń is a city with a long history, and the historical buildings such as palaces and temples are all well-preserved, the historical use of heavy metals in pigments, wood preservation, and brassware would also play an important role for their accumulation in urban soils around.

Conclusion
To assess the status of environmental quality in the completely and partially sealed and non-sealed soils, three indicators, I geo , EF, and PLI, and multifactorial statistical analysis, were used. The results indicate that the artificial sealing in urban areas slightly affects the content of heavy metals in soils. However, based on PLI, it was found that the sealing has influence on soil properties and non-sealed soil is the most exposed to the accumulation of pollutants. Determination of EF index has also become an effective indicator used to distinguish the natural origin from anthropogenic sources of heavy metals. In our studies, it was found that in completely sealed soils, heavy metals (Hg, Cu, Pb) were of anthropogenic origin. Similar pattern was detected for I geo ; thus, for completely sealed soils, it decreased in the following order: Pb > Cu > Hg > Cd > Zn > Ni > Cr, for the partially sealed soil: Pb > Cd > Hg > Zn > Cu > Ni > Cr, while for the reference soil: Pb > Hg > Zn > Cu > Cd > Ni > Cr.
Open Access This article is distributed under the terms of the Creative Comm ons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.