Next Article in Journal
Ecology and Esthetics, Esthetic Ecology and the Ecological Esthetic in the Landscape: Contributions to the Apparent TongueTwister
Next Article in Special Issue
Identifying Influencing Factors of Agricultural Soil Heavy Metals Using a Geographical Detector: A Case Study in Shunyi District, China
Previous Article in Journal
Divining the Future: Making Sense of Ecological Uncertainty in Turkana, Northern Kenya
Previous Article in Special Issue
Sources of and Control Measures for PTE Pollution in Soil at the Urban Fringe in Weinan, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differences in the Content of Zn Fractions in the Profiles of Soils from Allotment and Domestic Gardens in South-Eastern Poland

by
Iwona Makuch-Pietraś
1,* and
Anna Wójcikowska-Kapusta
2
1
Department of Nature Conservation and Landscape Ecology, Institute of Agricultural Science, Land Management and Environmental Protection, University of Rzeszów in Rzeszów, 35-959 Rzeszów, Poland
2
Institute of Soil Science, Engineering and Environmental Science, University of Life Science in Lublin, 20-033 Lublin, Poland
*
Author to whom correspondence should be addressed.
Land 2021, 10(9), 886; https://doi.org/10.3390/land10090886
Submission received: 6 July 2021 / Revised: 17 August 2021 / Accepted: 20 August 2021 / Published: 24 August 2021

Abstract

:
The aim of the research was to show the distribution of fractions as well as bioavailability and the total forms of Zn in the profiles of soils from domestic gardens and family allotment gardens in six cities in south-eastern Poland. Results found that the level of Zn total form varied in the ranges from A horizon: 12.75–154.75 mg·kg−1 in sandy soils and 18.20–104.00 mg·kg−1 in silty soils. Accumulation of metals was assessed using concentration indices, Igeo, and the Cav/Ct and BF indices of bioavailable forms. The analysis took into account the role of organic matter as an important component in binding the analyzed metals in soil horizons subjected to long-term horticultural cultivation. In the two groups of sandy and silty soils distinguished according to their particle size distribution, horticultural treatments were found to exert a greater impact on sandy soils. Additionally, higher contents of the examined element were stated in the humus horizons, as indicated by the high values of concentration and Igeo indices showing high Zn pollution in the soils. The content of bioavailable zinc forms was significantly high, especially in soils with a higher metal content.

1. Introduction

Allotment gardens are one of the most interesting types of urban greenery, along with domestic gardens, which serve similar functions in detached housing estates. In addition to their recreational function, they give owners the opportunity to grow their own vegetables and fruit [1,2]. Large-scale multi-year investigations have proved heavy metal pollution as a problem that may affect these areas due to their location within city limits. However, most studies mainly focus on large cities or centers located in highly industrialized areas [3,4,5,6,7,8,9] where soils are actually contaminated with heavy metals. In the case of smaller cities or located in less industrialized areas, it is difficult to find items describing a case study. This is an indication to analyze the soils of allotment gardens in terms of heavy metal content in areas with less urbanization pressure than agglomerations. The problem was also found in south-eastern Poland, which represented mostly the agricultural area.
Additionally, attention should be paid to the specific properties, structure, and long-term use of gardens, which contribute to the transition of natural soils into anthropogenically transformed ones characteristic of a long-cultivated hortisol horizon [10]. The effects of such changes include increased organic matter content, changes in the sorption complex, and an improved soil pH value [11]. Due to the high impact of urban and industrial pollution, the topsoil may contain higher metal concentrations [2,11,12,13,14]. The distribution of metals in the profile may vary depending on its structure. An important role in the metal distribution is played by the well-developed deeper humus horizon with a large amount of organic matter and a more extensive sorption complex. Another differentiating property is the particle size distribution of the soil. Based on the above information, as well as the huge popularity of both home and allotment gardens, the experiment for the article was developed. An important aspect of this experiment is the separation of two groups based on soil composition.
To show differences between soil types and divergences in their organic matter content, not only the total content but also metal fractions are worth determining, using an appropriate sequential extraction such as BCR. Furthermore, simple methods based on widely used indices of contaminant concentrations and soil pollution can help to demonstrate the differences between the anthropogenic and natural sources of elements [15], which is extremely important for garden owners producing vegetables intended for consumption [16].
The aspects discussed above became a reason to try to consider the problem of heavy metals content in one work, which is a new contribution to the literature on the subject.
For proper assessment of the contamination risks associated with the horticultural use of urban soils in cities of south-eastern Poland, the following research objectives were formulated: (i) estimation of Zn content in soils from urban gardens, in comparison with soils from arable fields located outside city limits, using various contamination indices, (ii) estimation of bioavailable forms through fractionation of metals and determination of bioavailable forms, by means of the Rinkis method, (iii) determination of the distribution of trace elements in the soil profile, with particular emphasis on organic matter-bound forms, by means of the sequential extraction BCR procedure [17,18].

2. Materials and Methods

2.1. Study Site

The research was conducted in the south-east of Poland. The general climatic conditions are characterized by an average temperature of around 9 ℃ and precipitation of 600 mm [19]. There are some differences between cities in terms of physiography and geology. Biała Podlaska is situated on the Łuków Plain, which is a flat, sandy area in the outflow zone of the glacial-river waters of the Warcia glaciation. Lubartów is located in the region of Lubartów Heights. This area is a denuded moraine surface with gravel inselbergs. Tarnobrzeg is a city located in the Vistula Lowland, which is part of the Sandomierz Basin. The valley is filled with quaternary river sediments with a thickness of several meters. Next to the flood terrace, there is a higher sandy terrace and a terrace covered with loess. Lublin is located on the Świdnik Plateau, which is part of the Lublin Upland. It is built of a thick series of loess resting on the upper cretaceous layers, degraded moraine, glaciofluvial sands, and gravels. Przemyśl is part of the Przemyśl Foothills, which is the outer part of the Carpathian Mountains, built of East-Carpathian flysch. Zamość lies in the Zamość Basin area, which is an extensive denudation depression, formed in less resistant upper cretaceous marls and chalk [20].

2.2. Sampling Site

Six cities (Figure 1) located in south-eastern Poland were selected for the study. A survey was conducted to find appropriate study sites in each town. The survey enabled the selection of one allotment garden in a complex of Family Allotment Gardens, one domestic garden at a detached house, and a control site (arable field) outside the city limits. Both the allotments and domestic gardens were selected based on an interview with owners with at least 20-years of horticultural activity in the selected criterion. A similar criterion was used for the crop field. An important element in the choice of the gardens was the mode of horticultural cultivation, that is, the use of only organic fertilizers, often produced by the owner, representing high crop culture. The crop field was selected after an interview with the owner, whose arable land was subjected to extensive cultivation by crop rotation. Moreover, the arable field was located outside the city, with no high-traffic transportation routes nearby. Potential heavy metal contamination was a criterion for the selection of the cities. The south-east of Poland is mainly an agricultural region, with a poorly developed industrial infrastructure.
Two groups of soils were distinguished according to the particle size distribution from six cities. One group comprised sandy soils group of light soils. These sand and loamy sand soils were sampled in Biała Podlaska, Lubartów, and Tarnobrzeg (Table S1). The second group was represented by soils from the group of medium soils, made from silt, mostly with the particle size distribution of silt loam, sampled in Lublin, Przemyśl, and Zamość (Table S2). The light sandy soils sampled from the analyzed cities represent podzolic, brunic, and gleying types, whereas the medium silty soils are luvisols [21]. Soil classification and properties are presented in Supporting Materials.

2.3. Soil Sampling and Basic Soil Properties

A soil profile was made in each domestic, allotment garden, and in the arable field. Samples were taken for analysis from each genetic horizon. Samples were composited of several samples from each horizon. The material was air-dried and then passed through a 2-mm mesh sieve.
Basic physicochemical properties were determined in the soil material. The particle size distribution was determined by means of the Bouyoucos and Casagrande areometric method, modified by Prószyński. Soil fractions and formations were identified according to PTGleb. [22]. Soil pH was determined potentiometrically in water and in a 1 M KCl solution. The content of organic matter was measured using the Tiurin method and the total nitrogen content was determined by the Kjeldahl method. The sorption capacity of the soils was calculated as CEC extracted with ammonium acetate at pH = 7.0, and the hydrolytic acidity was determined by the Kappen method [23]. All samples were analyzed in triplicate, and the final content was always considered as the average value of three measurements made in repeatability conditions.

2.4. Heavy Metal Determination and Sequential Extraction Procedure

The total copper and zinc content in the soil samples was determined by extraction using a mixture of concentrated acids (nitric and perchloric acid 1:1). To demonstrate the relationship between zinc contents and organic matter in the soils subjected to intensive horticultural use, the sequential extraction procedure proposed by BCR [17,18] was applied and the content of four fractions was determined. The scheme of the procedure is presented in Table 1. Additionally, to assess the potential bioavailability of zinc, bioavailable forms of this element were determined by the Rinkis method in 0.1 M HCl [23]. All forms of the element were determined by the ICP-AES method, using Leeman PS 950 apparatus. Digestion was performed in the DKL 20 Fully Automatic Digestion Unit (VELP SCIENTIFICA). The reagents used in the analyses were analytically pure, and the water was deionized. The calibration curves and the concentration measurements were performed for each metal separately using the certified single element standards for ICP, manufactured by Inorganic Ventures (Lakewood, NJ, USA). All samples were digested and analyzed in triplicate, and the final content was always considered as the average value of three measurements made in repeatability conditions.

2.5. Indices of Pollution

Based on the results, the indices of the distribution and concentration of zinc in the analyzed soils were calculated [15,24]:
C f = Cs Cn
Cs—concentration of the heavy metal;
Cn—background level of heavy metal.
The geoacumulation index (Igeo) was calculated from the total content of heavy metal determined in the profile horizons in relation to the parent material, using the following formula [25]:
I geo = log 2   [ C 1.5   B ]
C—concentration of the heavy metal in the horizon
B—background level of heavy metal.
The following classes were distinguished based on the value of the index:
  • ≤0: unpolluted;
  • 0–1: unpolluted to moderately polluted;
  • 1–2: moderately polluted;
  • 2–3: moderately to highly polluted;
  • 3–4: highly polluted;
  • 4–5: highly to extremely highly polluted;
  • ≥5: extremely highly polluted;
The ratio of bioavailable Zn forms to the total content of these elements was calculated [26].
( C av / C t ) · 100
Cav—content of the bioavailable form of the element;
Ct—total content of the element.
Based on the results of the sequential extraction procedure, the BF value was calculated using the following formula [27]:
BF = ( C bio · 100 ) / C total
where:
Cbio—concentrations of bioavailable forms of the metal in the studied soils, i.e., forms contained in fractions I–III of the sequential extraction of soil
Ctotal—the total concentration of the metal in the studied soil.

2.6. Statistical Analysis

For the topsoil levels, statistical analyses were carried out to check whether the grouping variables, that is, particle size distribution and land use, discriminate between the properties and parameters of the soils. To assess the effect of soil properties on the content and forms of heavy metal, Pearson’s correlation coefficients between selected data were calculated. The results were analyzed in the Statistica 13.3 program at a significance level of α ≤ 0.05. Additionally, principal component analysis was performed using Canoco 5 software.

3. Results

3.1. Basic Soil Properties

The mode of land use in relation to the light soils was a crucial determinant of such properties as pH, sorption capacity with its components, or organic carbon content. The physicochemical properties of the humus horizons were more favorable in the horticultural samples, compared to the arable soil samples; however, no such relationships were found in the medium soils (Tables S1 and S2).

3.2. Zn Content

3.2.1. Total Zn Content

All topsoil layers of the sandy soils exhibited a higher total zinc level than their parent material, indicating a potential impact of the location and land use. The highest content was determined in the soil from the domestic garden in Biała Podlaska and from both types of gardens in Tarnobrzeg (Table S3).
Regardless of the mode of land use, all medium soils were characterized by higher total zinc contents in the humus horizons, compared to that in the parent material. In some study sites, substantially higher Zn accumulation was observed in the surface layer of the allotment soil in Przemyśl and in the soil from the domestic garden and arable field in Zamość (Table S4). With regard to the study sites with the sandy soils, the distribution of the zinc content in their profiles was more variable and the profiles differed from each other, but no specific pattern was observed (Tables S3 and S4).

3.2.2. Bioavailable Zn Forms

The humus horizons in the sandy soils sampled from the domestic gardens were characterized by high amounts of zinc soluble in 1 mol HCl·dm−3, with the highest content, that is, 114.5 and 115.2 mg·kg−1, determined in the soil from the domestic garden in Biała Podlaska. In terms of the threshold values, the humus horizons of the sandy soils were characterized by high zinc content, with the exception of the arable field soil in Biała Podlaska. In the group of medium soils, high zinc contents were detected in the domestic garden soil in Przemyśl as well as the allotment garden and arable field soil in Zamość. The other soils exhibited a moderate abundance of this component [28]. The content of bioavailable zinc declined with the profile depth in all soils. In terms of the mode of land use, the profiles of the horticultural soils were characterized by a higher level of this element (Table S3).
In the case of the silty soils, the highest content of zinc soluble in 1 mol HCl·dm−3 was exhibited by the humus horizons in the garden soils, and the values decreased proportionally to the profile depth. In turn, the humus horizons of the arable field soils contained substantially lower contents of bioavailable Zn, likewise the other genetic horizons in the profiles (Table S4).

3.2.3. Indices of Pollution

Regardless of the mode of land use, the values of the tested zinc Cf index ranged from 0.63 to 17.34 in the sandy soils. A substantial decrease in the index value was noticed in the case of the garden soil profiles, whereas smaller differences between these values were found for the agricultural soils (Table S3).
The silty soil profiles differed substantially in the values of the zinc concentration index, depending on the mode of land use. The highest values of this indicator were found in the soil profiles from the domestic gardens (0.62–11.30), whereas the allotment soil samples were characterized by its lowest values (0.48–2.45). The zinc concentration index in the arable soils ranged from 0.80 to 4.00, and its value decreased with depth (Table S4).
High Igeo values for zinc in the light soils were calculated only in the case of the topsoil layers, especially in the domestic garden in Biała Podlaska (the highest value: 3.53), which indicates high pollution of the soils. The values in the other soils ranged between 0.00 and 0.80 and classified them as unpolluted (Table S3).
In the group of the silty soils, high values of the Igeo index were obtained in the case of the domestic garden soil from Zamość (1.39–2.91) classified as moderately polluted and in the Ap horizon of the arable soil (1.42) representing the same pollution class. The other profiles exhibited index values in the range from 0.00 to 0.90 and represented the unpolluted soil class (Table S4).
The humus horizons of the silty garden soils exhibited the highest values of the ratio of the bioavailable zinc form to its total content (over 60), whereas lower values were reported for the Ap horizons of the arable soil (up to 40). Regardless of the mode of land use, the solubility of zinc declined with depth (Table S3).
The Ap horizons of the silty soils were characterized by high values of the ratio of the bioavailable zinc form to the total content (even above 100 in the domestic garden in Lublin). Zinc solubility decreased with depth in a majority of the soil profiles (Table S4).
The BF index in the group of light soils was lower in the arable soils. In most profiles, it decreased with depth, which proves the accumulation of bioavailable forms in the topsoil horizons and indicates an important role that may be played by humus in garden soils (Table S3). No such tendencies were found for the silty soils. The BF index values were similar in the entire profile and slightly higher in the case of the garden soils (Table S4).

3.2.4. Sequential Extraction of Zn

The content of water-soluble and exchangeable zinc fraction I in the sandy soils was in the range of 0.96–32.11 mg·kg−1 (Table S3). The fraction is important due to its bioavailability; hence, humus horizons should contain its optimal amounts. The topsoil horizons of the garden soils were characterized by greater amounts of this fraction in comparison with the top layers of the arable soils. The soils from the domestic garden in Biała Podlaska exhibited the highest content of this fraction. The content of exchangeable and water-soluble zinc declined in soil samples from sites subjected to both modes of land use (Figure 2, Table S3).
The percentage proportion of the zinc fraction bound to iron and manganese oxides in the light soils ranged from 7.85 to 41.38% of the total content (Figure 2), which was in the range of 1.09–53.24 mg·kg−1 (Table S3). In the group of garden soils, the humus horizons from the domestic garden in Biała Podlaska had the highest levels of this form of zinc, and this trend was also observed in the other profiles of the garden soils. The topsoil layers in the arable soils exhibited lower contents of zinc fraction II (Figure 2, Table S3).
The content of Zn fraction III in the sandy soils ranged from 1.44 to 24.06 mg·kg−1 (Table S3), which accounted for 5.69–27.67% of the total content (Figure 2). The humus-associated fraction is the best indicator in sequential extraction, showing the relationship between heavy metal contamination and high organic matter content. The humus horizons in the garden soil profiles exhibited a higher level of humus-associated zinc fraction III than the deeper horizons. As in the case of the fractions discussed above, the highest content was determined in the humus horizons of the domestic garden soil from Biała Podlaska. The arable soils contained substantially smaller amounts of this Zn fraction; however, there was a similar trend as in the case of the horticultural soil (Figure 2, Table S3).
The content of residual fraction IV in the light soils was in the range of 1.17–47.04 mg·kg−1 (Table S3), with the percentage proportion in the total content varying in a wide range of 6.21–82.82% (Figure 2). The residual fraction in the humus horizons of the garden soils represented a lower percentage than the other fractions, especially in two study sites, that is, Biała Podlaska and Tarnobrzeg, probably due to the greater proportion of the humus-associated fraction. In comparison with the humus horizon, the parent material exhibited a greater proportion of the residual Zn fraction, except for the arable soils sampled in Biała Podlaska and Tarnobrzeg (Figure 2). The content of fraction IV in most profiles declined with depth (Table S3).
The percentage proportion of zinc fraction I in the total content in the silty soils ranged from 0.62 to 27.22% (Figure 3) and its content was in the range of 0.33–5.95 mg·kg−1 (Table S4). The content of this zinc fraction decreased with depth in the soils from both types of land use. The mode of land use did not have a clear effect on the content of this fraction (Figure 3, Table S4).
The percentage of zinc fraction II bound to iron and manganese oxides ranged from 1.57 to 25.94% in the garden soils from 1.03 to 19.64% in the agricultural soils (Figure 3). Its content was in the range of 0.29–11.02 mg·kg−1 (Table S4). In both types of land use, the content of the zinc fraction bound to iron and manganese oxides in the silty soil profiles was found to vary and was mainly determined by the soil formation process (Table S4).
Regardless of the land use, the content of zinc fraction III was in the range of 1.54–7.05 mg·kg−1 (Table S4), which represented from 4.41% to 24.04% (Figure 3). In contrast to the light soils, the content of this fraction in the humus horizons was not substantially higher than that in the other genetic horizons. There was also no effect of the mode of land use on the content of the humus-bound fraction in the analyzed soils (Table S4).
In the medium soils, the content of the residual zinc fraction in the humus horizons was similar and higher than in the deeper horizons, regardless of the mode of land use. As in the case of the total zinc content, the humus horizon in the domestic garden soil in Zamość contained the highest amounts of the residual form of this element (Table S4). Regardless of the type of land use, the level of this fraction was in the range of 1.97–87.68 mg·kg−1 (Table S4). The percentage proportion of residual zinc in the total content was slightly higher in the arable soils (47.31–91.09%) than in the garden soils (30.14–88.84%) (Figure 3).

3.3. Statistical Analysis

In the sandy soils, the organic carbon content was significantly correlated with the zinc content indices: Cf (r = 0.860 at p < 0.01), Ca/Ct (r = 0.689 at p < 0.05), and Igeo (r = 0.890 at p < 0.01). It was also correlated with the total Zn content, bioavailable forms, and fraction II of this element. Significant correlations were also found for the total content and bioavailable form with Cf (r = 0.934 and 0.965 at p < 0.001) and Igeo (r = 0.915 and 0.954 at p < 0.001). There was also a statistically significant correlation between the total content and the bioavailable zinc form (r = 0.982 at p < 0.001) and between fraction II and pH measured in water (r = 0.669 at p < 0.05) (Table 2).
In comparison to the sandy soils, the organic matter content in the medium soils had no significant effect on the content of zinc, although the total and bioavailable forms exhibited a positive correlation (r = 0.825 at p < 0.01 and r = 0.741 at p < 0.05, respectively). Additionally, these forms were characterized by a significantly high coefficient of correlations with each other (r = 0.942 at p < 0.001) and with Cf (r = 0.909 and 0.962, respectively at p < 0.001) and Igeo (r = 0.886 at p < 0 01 and 0.915 at p < 0.001). The correlation analysis revealed only one statistically significant correlation between the silt fraction content and Ca/Ct (r = 0.855 at p < 0.01) (Table 3).
The PCA results for the parameters of soil Zn are presented in Table 4 and Figure 4. The first component, which explains 34.90% of the total variance, clearly discriminated most of the parameters of soil physicochemical properties and zinc content. Moreover, the content of total zinc and its bioavailable form, the Cf and Ca/Ct indices, and the Corg content were discriminated by the second component, although to a lesser extent. These variables are located close to each other and have the same direction, which indicates their positive and strong correlation. Among the analyzed forms of zinc, residual fraction IV exhibited an opposite direction and sense of the vector, which indicated a negative correlation with the other parameters and a low effect of the physicochemical properties on this form of the element. The BF index expressing the mobility of Zn forms was positively correlated with fractions I and II, which were characterized by the highest bioavailability in soil. In turn, fraction III was negatively correlated with these parameters and positively correlated with the reaction and content of the silt fraction, which allows the conclusion that it is dependent on these parameters. The analyzed soil types form two groups: the medium soils are located closer to the first principal component, while the light soils are more closely related to the second principal component. In the sandy soils, there is a lesser effect of the components on the agricultural land use, which forms a group located closer to the axis of the second variable axis. No similar relationships were found in the case of the silty soils.

4. Discussion

As indicated by literature reports on the problem of soil contamination with heavy metals in gardens, it is generally believed that such sites are contaminated by substances from neighboring houses, roads, and urban or industrial areas [11]. As shown by the author, zinc is one of the most common contaminants, and their form present in soils is influenced by, for example, soil physicochemical properties such as pH, organic matter content, calcium carbonate, and silt fraction. Moreover, it is worth emphasizing that horticultural soils may exhibit higher contents of organic matter than arable soils, which promotes the migration of heavy metals [11].
Pollution assessment has been carried out in the humus horizons of the sandy and silty soils which were characterized with natural zinc contents (except for the domestic garden in Biała Podlaska and the domestic garden in Zamość) [29,30]. This differentiation indicates the possibility of local changes in the accumulation of heavy metals.
The total zinc content decreased with the profile depth in all light and silty soils. This differentiation indicates the possibility of local changes in the accumulation of heavy metals by Gorlach et al. [31] Kabata-Pendias [32], and Kawałko and Chodak [33]. This tendency was reflected in the high zinc concentration indices in the topsoil layers. As suggested by Świercz and Smorzewska [34], high values of this index may be related to many anthropogenic factors and geological backgrounds. In the case of soils located in the area of the influence of the city, the anthropogenic factor becomes very significant. The second used index was Igeo which was characterized by quite high values in the humus horizons and similar results were reported by Ciupa et al. [35], Świercz and Smorzewska [34], and Świercz and Zajęcka [36]. This has been a confirmation of the concentration ratio used earlier. In the case of soils where horticulture is used, land use may be an additional factor influencing the Zn content. In terms of the mode of land use, the zinc content was higher in the horticultural than agricultural soils only in the group of the sandy soils. These results were in agreement with those reported by Kelly and Thornton [37], who found accumulation of this element in intensively cultivated soils. Similar conclusions were formed by Kabała et al. [38], who compared soils from horticultural and agricultural land use. Additionally, the zinc content was positively correlated with organic matter, which was also demonstrated by many authors [1,39,40,41,42,43,44,45,46]. From the literature on the problem of heavy metals in garden soils, it is difficult to find items describing the differences between types, or mechanical groups. Whereas in the arable soils, were stated substantial differences in the total zinc content between the soil types. This was also confirmed by many authors in their analyses of the content of this element in various cultivated soils [32,39,42,46,47,48,49,50,51].
An interesting aspect of the work is the sequential analysis, showing the diversity of forms associated with the organic substance, as well as mobile fractions, which result from a greater accumulation of Zn. The surface horizon of the soils made of sands was characterized by a higher content of mobile fractions, although the fraction associated with the organic substance was lower. The residual fraction dominated in the surface horizons of soils made of silt, regardless of the way of use. Attention was drawn to the distribution of fractions in the deeper levels of the profile, which reflected the natural soil formation process. These results are supported by many literature positions. Similarly, Moćko and Wacławek [52] reported the dominance of the residual zinc fraction and the organic matter-bound fraction in garden soils from Opole, whereas mobile fractions II and I represented the lowest proportion. In China, soils where vegetables were grown exhibited the dominance of the organic matter-bound fraction, whereas the proportion of the exchangeable and residual fractions was highly variable. The content of mobile fractions (I and II) may vary depending on the conditions prevailing in soils [53]. In arable soils, Herencia et al. [54] showed constant amounts of the residual fraction in different fertilization systems. After the application of compost, the sum of fractions I, II, and III and the amount of the organic matter-bound fraction were found to increase. These results prove the high affinity of zinc for the organic carbon content. Industries may be another determinant of the content of mobile zinc fractions, as confirmed by Kaasalainen and Yli-Halla [55] in their investigations of arable soils of Finland. In the case of soils of south-eastern Poland, it could be stated that the lowest proportion of the Zn fractions, especially fractions I and II, in the analyzed soils may be associated with the limited impact of the industry-related factor. In turn, the aspect of the urbanization impact should be emphasized. As reported by Imperato et al. [56], urban soils were characterized by a smaller share of individual fractions in the total zinc content. The organic fraction was the smallest of all mobile fractions, compared to the residual fraction. As reported by other researchers, urban soils were characterized by the highest amounts of mobile fractions in contrast to arable and garden soils [57,58].
Similar conclusions were formulated by Dąbkowska-Naskręt et al. [59] in their study of urban soils from city parks, which were characterized by the highly variable mobility of Zn. The BF index is a very interesting example of reflecting the mobility of heavy metals. The BF index reported by Bielicka-Giełdoń et al. [27] had a value of 50–70% and was more comparable with the light soils, where it reached as high as 90% in the topsoil layer. The bioavailability of zinc in the analyzed medium soils significantly decreased and was below 40% also in the humus horizons. This indicates strong dependence of zinc mobility on the mineral part of soils.

5. Conclusions

Based on research, it was found that the vast majority of the analyzed soils had natural zinc levels. Sites selected for the study were cities with a potentially lower degree of urban pressure, however, at several research points, there was an increased content of the investigated metal.
Moreover, the accumulation of zinc in profiles was reflected by the concentration index, which was higher in the humus horizons, which indicates the influence of the land use combined with the location of the researched sites, which was confirmed by the higher total content in garden soils than in cultivated soils. Moreover, the concentration indices were higher in the sandy than silty soils. The Igeo index was similar, which exhibited quite high values in the humus horizons, and in some points classified the soil as even highly polluted in some study sites (Biała Podlaska—domestic garden; Zamość—domestic garden and arable field), which was confirmed by the comparison of the total content with permissible soil contamination with heavy metals. Based on the above information, it can be concluded that the location of the soil had an important role, and the humus level developed during horticulture may contribute to the accumulation of heavy metals.
To show the connection of soil organic matter formed by the horticulture with Zn content, the sequential extraction was used to distinguished fractions with the best bioavailability, through fractions accumulated in humus, to the residual fraction. The residual fraction of Zn dominated in both soil groups. The horticultural use effect of the soil was associated with greater accumulation of zinc in the humus horizons, in particular in the light soils, which was confirmed by the statistically significant correlations with the organic carbon content.
Generally, the availability of zinc determined on the basis of the BF index was highly influenced by the particle size distribution, especially in the case of light soils, where the BF values were higher in the humus horizons. This allows a conclusion that horticultural production in urban areas requires more attention in the case of sandy soils, as they may contribute to further migration of heavy metals.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/land10090886/s1, Table S1: Basic physicochemical properties of sandy soils, Table S2: Basic physicochemical properties of silty soils, Table S3: Content of the forms and values of Zn pollution indices in sandy soils, Table S4: Content of the forms and values of Zn pollution indices in silty soils.

Author Contributions

Conceptualization, I.M.-P. and A.W.-K.; methodology, I.M.-P. and A.W.-K.; software, I.M.-P.; validation, I.M.-P. and A.W.-K.; formal analysis, I.M.-P.; investigation, I.M.-P. and A.W.-K.; resources, I.M.-P. and A.W.-K.; data curation, I.M.-P. and A.W.-K.; writing—original draft preparation, I.M.-P. and A.W.-K.; writing—review and editing, I.M.-P. and A.W.-K.; visualization, I.M.-P.; supervision, I.M.-P.; project administration, I.M.-P.; funding acquisition, A.W.-K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Informed consent was obtained from all subjects in.

Data Availability Statement

Not applicable.

Acknowledgments

The authors gratefully acknowledge the support of Alicja Strzałka, Jerzy Michalczuk and Bernadetta Ortyl for help and willingness.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burghardt, W.; Schneider, T. Bulk density and content, density and stock of carbon, nitrogen and heavy metals in vegetable patches and lawns of allotments gardens in the northwestern Ruhr area, Germany. J. Soils Sediments 2018, 18, 407–417. [Google Scholar] [CrossRef]
  2. Leitão, E.; Cameira, M.R.; Costa, H.D.; Pacheco, J.M.; Henriques, M.J.; Martins, L.L.; Mourato, M.P. Environmental quality in urban allotment gardens: Atmospheric deposition, soil, water and vegetable assessment at Lisbon city. Water Air Soil Pollut. 2018, 229, 1–22. [Google Scholar] [CrossRef]
  3. Calmuc, V.A.; Calmuc, M.; Arseni, M.; Topa, C.M.; Timofti, M.; Burada, A.; Iticescu, C.; Georgescu, L.P. Assessment of heavy metal pollution levels in sediments and of ecological risk by quality indices, applying a case study: The Lower Danube River, Romania. Water 2021, 13, 1801. [Google Scholar] [CrossRef]
  4. Golia, E.E.; Papadimou, S.G.; Cavalaris, C.; Tsiropoulos, N.G. Level of contamination assessment of potentially toxic elements in the urban soils of Volos City (Central Greece). Sustainability 2021, 13, 2029. [Google Scholar] [CrossRef]
  5. Herbón, C.; Barral, M.T.; Paradelo, R. Potentially toxic trace elements in the urban soils of Santiago de Compostela (Northwestern Spain). Appl. Sci. 2021, 11, 4211. [Google Scholar] [CrossRef]
  6. Horváth, A.; Csáki, P.; Szita, R.; Kalicz, P.; Gribovszki, Z.; Bidló, A.; Bolodár-Varga, B.; Balázs, P.; Winkler, D. A Complex soil ecological approach in a sustainable urban environment: Soil properties and soil biological quality. Minerals 2021, 11, 704. [Google Scholar] [CrossRef]
  7. Pecina, V.; Brtnicky, M.; Balkova, M.; Hegrova, J.; Buckova, M.; Baltazar, T.; Licbinsky, R.; Radziemska, M. Assessment of soil contamination with potentially toxic elements and soil ecotoxicity of botanical garden in Brno, Czech Republic: Are urban botanical gardens more polluted than urban parks? Int. J. Environ. Res. Public Health 2021, 18, 7622. [Google Scholar] [CrossRef]
  8. Petukhov, A.; Kremleva, T.; Petukhova, G.; Khritokhin, N. Biochemical responses of medicinal plant Tussilago farfara L. to elevated heavy metal concentrations in soils of urban areas. Toxics 2021, 9, 171. [Google Scholar] [CrossRef] [PubMed]
  9. Morillo, E.; Romero, A.S.; Madrid, L.; Villaverde, J.; Maqueda, C. Characterization and sources of PAHs and potentially toxic metals in urban environments of Sevilla (Southern Spain). Water Air Soil Pollut. 2008, 187, 41–51. [Google Scholar] [CrossRef]
  10. Charzyński, P.; Bednarek, R.; Hudańska, P.; Świtoniak, M. Issues related to classification of garden soils from the urban area of Toruń, Poland. Soil Sci. Plant Nutr. 2018, 64, 132–137. [Google Scholar] [CrossRef]
  11. Alloway, B.J. Contamination of soils in domestic gardens and allotments: A brief overview. Land Contam. Reclam. 2004, 12, 179–187. [Google Scholar] [CrossRef]
  12. Bechet, B.; Joimel, S.; Jean-Soro, L.; Hursthouse, A.; Agboola, A.; Leitão, T.E.; Costa, H.; Cameira, M.R.; Le Guern, C.; Schwartz, C.; et al. Spatial variability of trace elements in allotment gardens of four European cities: Assessments at city, garden, and plot scale. J. Soils Sediments 2018, 18, 391–406. [Google Scholar] [CrossRef] [Green Version]
  13. Bretzel, F.; Calderisi, M.; Scatena, M.; Pini, R. Soil quality is key for planning and managing urban allotments intended for the sustainable production of home-consumption vegetables. Environ. Sci. Pollut. Res. 2016, 23, 17753–17760. [Google Scholar] [CrossRef] [PubMed]
  14. Giusti, L. Heavy metals in urban soils of Bristol (UK). Initial screening for contaminated land. J. Soils Sediments 2011, 11, 1385–1398. [Google Scholar] [CrossRef]
  15. Weissmannová, H.D.; Pavlovský, J. Indices of soil contamination by heavy metals—methodology of calculation for pollution assessment (minireview). Environ. Monit. Assess. 2017, 189, 616. [Google Scholar] [CrossRef] [PubMed]
  16. Barbieri, M. The importance of enrichment factor (EF) and geoaccumulation index (Igeo) to evaluate the soil contamination. J. Geol. Geophys. 2016, 5, 237. [Google Scholar] [CrossRef]
  17. Thomas, R.P.; Ure, A.M.; Davidson, C.M.; Littlejohn, D. Three stage sequential extraction Procedure for the determination of metals in river sediments. Anal. Chem. Acta 1994, 286, 423–429. [Google Scholar] [CrossRef]
  18. Ure, A.M.; Quevauviller, P.H.; Muntau, H.; Griepink, B. Speciation of heavy metals in soils and sediments. An account of the improvement and harmonization of extraction techniques undertaken under the auspices of the BCR of the Commission of the European Communitie. Int. J. Environ. Anal. Chem. 1993, 51, 135–151. [Google Scholar] [CrossRef]
  19. Klimat Polski 2020. IMGW. Available online: http://klimada.mos.gov.pl/zmiany-klimatu-w-polsce/tendencje-zmian-klimatu/ (accessed on 17 August 2021).
  20. Solon, J.; Borzyszkowski, J.; Bidłasik, M.; Richling, A.; Badora, K.; Balon, J.; Brzezińska-Wójcik, T.; Chabudziński, Ł.; Dobrowolski, R.; Grzegorczyk, I.; et al. Physico-geographical mesoregions of Poland: Verification and adjustment of boundaries on the basis of contemporary spatial data. Geogr. Pol. 2018, 91, 143–170. [Google Scholar] [CrossRef]
  21. Soil Science Society of Poland, Commission on Soil Genesis, Classification and Cartography. Polish Soil Classification; Wydawnictwo Uniwersytetu Przyrodniczego We Wrocławiu, Polskie Towarzystwo Gleboznawcze: Wrocław-Warszawa, Poland, 2019. [Google Scholar]
  22. Polish Standard. Soil and Mineral Soil Materials—Sampling and Determination of Particle Size Distribution; PN-R-04032; Polish Committee for Standardization: Warszawa, Poland, 1998.
  23. Ostrowska, A.; Gawlinski, S.; Szczubiałka, Z. Methods of Analysis and Evaluation of Properties of Soils and Plants; Institute of Environmental Protection: Warsaw, Poland, 1991. [Google Scholar]
  24. Weissmannová, H.D.; Pavlovský, J.; Chovanec, P. Heavy metal contaminations of urban soils in Ostrava, Czech Republic: Assessment of metal pollution and using principal component analysis. Int. J. Environ. Res. 2015, 9, 683–696. [Google Scholar]
  25. Kowalska, J.; Mazurek, R.; Gąsiorek, M.; Setlak, M.; Zaleski, T.; Waroszewski, J. Soil pollution indices conditioned by medieval metallurgical activity—A case study from Krakow (Poland). Environ. Pollut. 2016, 218, 1023–1036. [Google Scholar] [CrossRef] [PubMed]
  26. Dąbkowska–Naskręt, H.; Kędzia, W. Mobilność miedzi w uprawnych czarnych ziemiach kujawskich. Zesz. Naukowe. Pol. Akad. Nauk. Kom. Nauk. Przy Prezydium PAN Człowiek Sr. 1996, 14, 51–56. [Google Scholar]
  27. Bielicka-Giełdoń, A.; Ryłko, E.; Żamojć, K. Distribution, bioavailability and fractionation of metallic elements in allotment garden soils using the BCR sequential extraction procedure. Pol. J. Environ. Stud. 2013, 22, 1013–1021. [Google Scholar]
  28. Zalecenia Nawozowe, C.I. Liczby Graniczne do Wyceny Zawartości w Glebach Makro- i Mikroelementów; IUNG Puławy. Seria P(44): Puławy, Poland, 1990. [Google Scholar]
  29. Rozporządzenie Ministra Środowiska. w Sprawie Sposobu Prowadzenia Oceny Zanieczyszczenia Powierzchni Ziemi; Dz.U. z dn. 5 września 2016 Poz. 1395; Środowiskaed, M., Ed.; Ministry of the Environment of Poland: Warsaw, Poland, 2016.
  30. Kabata-Pendias, A.; Pendias, H. Biogeochemistry of Trace Elements; PWN: Warsaw, Poland, 1999. [Google Scholar]
  31. Gorlach, E.; Brydak, K.; Gambuś, F. Distribution of heavy metals in soil profiles of the Cracow region. Pol. J. Soil Sci. 1993, 34, 35–42. [Google Scholar]
  32. Kabata–Pendias, A. Zawartość metali ciężkich w glebach uprawnych Polski. Pamięt Puł. 1981, 74, 101–111. [Google Scholar]
  33. Kawałko, D.; Chodak, T. Zawartość niektórych metali ciężkich w glebach i warzywach ogródków działkowych Wrocławia. Zesz. Probl. Postępów Nauk. Rol. 1996, 434, 949–954. [Google Scholar]
  34. Świercz, A.; Smorzewska, E. Variations in the zinc and lead content in surface layers of urban soils in Kielce (Poland) with regard to land use. J. Elem. 2015, 20, 449–461. [Google Scholar]
  35. Ciupa, T.; Suligowski, R.; Kozłowski, R. Trace metals in surface soils under different land uses in Kielce city, south-central Poland. Environ. Earth Sci. 2020, 79, 14. [Google Scholar] [CrossRef]
  36. Świercz, A.; Zajęcka, E. Accumulation of heavy metals in the urban soils of the city of Skarżysko-Kamienna (Poland) with regard to land use. Carpathian J. Earth Environ. 2018, 13, 249–266. [Google Scholar] [CrossRef]
  37. Kelly, J.; Thornton, I. Urban geochemistry: A study of the influence of anthropogenic activity on the heavy metal content of soils in traditionally industrial and non-industrial areas of Britain. Appl. Geochem. 1996, 11, 363–370. [Google Scholar] [CrossRef]
  38. Kabała, C.; Chodak, T.; Szerszeń, L.; Karczewska, A.; Szopka, K.; Fratczak, U. Factors influencing the concentration of heavy metals in soils of allotment gardens in the city of Wrocław. Fresenius Environ. Bull 2009, 18, 622–630. [Google Scholar]
  39. Chojnicki, J.; Czarnowska, K. The changes of the contents of total and readily soluble phosphorus and Zn, Cu, Pb, Cd in agricultural soils under intensive cultivation. Soil Sci. Ann. 1993, 44, 99–111. [Google Scholar]
  40. Dąbkowska–Naskręt, H.; Kobierski, M. Metale ciężkie i ich mobilność w czarnych ziemiach aglomeracji Inowrocławia. Zesz. Probl. Postępów Nauk. Rol. 1998, 460, 259–267. [Google Scholar]
  41. Grzebisz, W.; Cieśla, L.; Komisarek, J.; Potarzycki, J. Geochemical assessment of heavy metals pollution of urban soils. Pol. J. Environ. Stud. 2002, 11, 493. [Google Scholar]
  42. Klimowicz, Z.; Melke, J. The content of heavy metals in soils in the vicinity of traffic roads using chosen stretches of road as examples. Soil Sci. Ann. 2000, 51, 36–45. [Google Scholar]
  43. Madrid, L.; Diaz-Barrientos, E.; Reinoso, R.; Madrid, F. Metals in urban soils of Sevilla: Seasonal changes and relations with other soil components and plant contents. Eur. J. Soil Sci. 2004, 55, 209–217. [Google Scholar] [CrossRef]
  44. Römkens, P.; Salomons, W. Cd, Cu and Zn solubility in arable and forest soils: Consequences of land use changes for metal mobility and risk assessment. Soil Sci. 1998, 163, 859–871. [Google Scholar] [CrossRef]
  45. Ruiz-Cortés, E.; Reinoso, R.; Díaz-Barrientos, E.; Madrid, L. Concentrations of potentially toxic metals in urban soils of Seville: Relationship with different land uses. Environ. Geochem. Health 2005, 27, 465–474. [Google Scholar] [CrossRef] [Green Version]
  46. Terelak, H.; Tujaka, A. Występowanie pierwiastków śladowych w glebach użytków rolnych województwa podkarpackiego. Zesz. Probl. Postępów Nauk. Rol. 2003, 493, 245–252. [Google Scholar]
  47. Czarnowska, K. Total content of heavy metals in parent rocks as reference background levels of soils. Soil Sci. Ann. 1996, 47, 43–50. [Google Scholar]
  48. Terelak, H.; Piotrowska, M.; Motowicka-Terelak, T.; Stuczyński, T.; Budzyńska, K. Zawartość metali ciężkich i siarki w glebach użytków rolnych Polski oraz ich zanieczyszczenie tymi składnikami. Zesz. Probl. Postępów Nauk. Rol. 1995, 418, 45–59. [Google Scholar]
  49. Terelak, H.; Tujaka, A.; Motowicka-Terelak, T. Trace element content (Cd, Cu, Ni, Pb, Zn) in farm—Land soils in Poland. Arch. Ochr. Środ. 2001, 27, 159–174. [Google Scholar]
  50. Ullrich, S.M.; Ramsey, M.H.; Helios–Rybicka, E. Total and exchangeable concentrations of heavy metals in soils near Bytom, an area of 180 Pb/Zn mining and smelting in Upper Silesia, Poland. Appl. Geochem. 1999, 14, 187–196. [Google Scholar] [CrossRef]
  51. Właśniewski, S. Pierwiastki śladowe w glebach ogrodów działkowych Rzeszowa. Zesz. Probl. Postępów Nauk. Rol. 2003, 493, 279–287. [Google Scholar]
  52. Moćko, A.; Wacławek, W. Three—Step extraction procedure for determination of heavy metals availability to vegetables. Anal. Bioanal. Chem. 2004, 380, 813–817. [Google Scholar] [CrossRef]
  53. Chao, W.; Liu, X.; Li, Z.; Pei, W.; Guo, Z. Pb, Cu, Zn and Ni concentration in vegetables in relation to their extractable fractions in soils in suburban areas of Nanjing, China. Pol. J. Environ. Stud. 2007, 16, 199–207. [Google Scholar]
  54. Herencia, J.F.; Ruiz, J.C.; Morillo, E.; Melero, S.; Villaverde, J.; Maqueda, C. The effect of organic and mineral fertilization on micronutrient availability in soil. Soil Sci. 2008, 173, 69–80. [Google Scholar] [CrossRef]
  55. Kaasalainen, M.; Yli-Halla, M. Use of sequential extraction to assess metal partitioning in soils. Environ. Pollut. 2003, 126, 225–233. [Google Scholar] [CrossRef]
  56. Imperato, M.; Adamo, P.; Naimo, D.; Arienzo, M.; Stanzione, D.; Violante, P. Spatial distribution of heavy metals in urban soils of Naples city (Italy). Environ. Pollut. 2003, 124, 247–256. [Google Scholar] [CrossRef]
  57. Davidson, C.M.; Urquhart, G.J.; Ajmone-Marsan, F.; Biasioli, M.; Costa Duarte, A.; Diaz-Barrientos, E.; Grčman, H.; Hossack, I.; Hursthouse, A.; Madrid, L.; et al. Fractionation of potentially toxic elements in urban soils from five European cities by means of a harmonised sequential procedure. Anal. Chim. Acta 2006, 565, 63–72. [Google Scholar] [CrossRef] [Green Version]
  58. Mossop, K.F.; Davidson, C.M. Comparison of original and modified BCR sequential extraction procedures for the fractionation of cooper, iron, lead, manganese and zinc in soils and sediments. Anal. Chim. Acta 2003, 478, 111–118. [Google Scholar] [CrossRef] [Green Version]
  59. Dąbkowska-Naskręt, H.; Różański, S.; Bartkowiak, A. Forms and mobility of trace elements in soils of park areas from the city of Bydgoszcz, north Poland. Soil Sci. Ann. 2016, 67, 73–78. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map of sampling in the study area.
Figure 1. Map of sampling in the study area.
Land 10 00886 g001
Figure 2. Percentage content of Zn fractions in sandy soils.
Figure 2. Percentage content of Zn fractions in sandy soils.
Land 10 00886 g002
Figure 3. Percentage content of Zn fractions in silty soils.
Figure 3. Percentage content of Zn fractions in silty soils.
Land 10 00886 g003
Figure 4. Principal component analysis (PCA) of two indices for the zinc content and parameters.
Figure 4. Principal component analysis (PCA) of two indices for the zinc content and parameters.
Land 10 00886 g004
Table 1. Scheme of the sequential analysis of metals according to BCR.
Table 1. Scheme of the sequential analysis of metals according to BCR.
FractionExtracion Solvent
IExchangeable and acid soluble0.11 M CH3COOH, pH = 2
IIReducible0.1 mol·dm−3 NH2OH · HCl, pH = 2
IIIOrganic30% H2O2 + 1 mol·dm−3 CH3COONH4, pH = 2
IVResidualHNO3 + HClO4
Table 2. Values of Pearson’s correlation coefficients of zinc forms and indices with the physicochemical properties of sandy soils.
Table 2. Values of Pearson’s correlation coefficients of zinc forms and indices with the physicochemical properties of sandy soils.
Properties of Studied SoilZn TotalZn Avail.pH H2OpH KClCrop<0.002
Zn total1.0000.982 *0.6160.6270.879 **0.009
Cf0.934 *0.965 *0.3810.3700.860 **−0.259
Zn available0.982 *1.0000.5990.5950.929 *−0.111
Ca/Ct0.3850.5090.6440.5950.689 ***−0.189
Igeo0.915 *0.954 *0.3890.3920.890 **−0.292
I fraction−0.226−0.081−0.121−0.206−0.029−0.551
II fraction0.4700.5820.669 ***0.6470.784 ***−0.196
III fraction−0.0350.111−0.070−0.1140.241−0.197
IV fraction−0.158−0.323−0.304−0.243−0.4940.395
BF0.1580.3230.3040.2430.494−0.395
Significance level: *—0.001; **—0.01; ***—0.05.
Table 3. Values of Pearson’s correlation coefficients of zinc forms and indices with the physicochemical properties of silty soils.
Table 3. Values of Pearson’s correlation coefficients of zinc forms and indices with the physicochemical properties of silty soils.
Properties of Studied SoilZn TotalZn Avail.pH H2OpH KClCrop<0.002
Zn total1.0000.942 *0.2060.1970.825 **−0.135
Cf0.909 *0.962 *0.2870.2720.638−0.356
Zn available0.942 *1.0000.1900.2330.741 ***−0.422
Ca/Ct0.3770.648−0.1020.0940.285−0.855 **
Igeo0.886 **0.915 *0.3050.3340.636−0.328
I fraction−0.611−0.352−0.218−0.152−0.539−0.577
II fraction−0.220−0.132−0.334−0.0980.042−0.480
III fraction−0.616−0.521−0.713 *−0.613−0.386−0.302
IV fraction0.5250.3300.3910.2220.3200.592
BF−0.525−0.330−0.391−0.222−0.320−0.592
Significance level: *—0.001; **—0.01; ***—0.5.
Table 4. Summary of the Zn results in PCA.
Table 4. Summary of the Zn results in PCA.
StatisticAxis 1Axis 2Axis 3Axis 4
Eigenvalues0.34900.31960.09770.0942
Explained variation (cumulative)34.9066.8676.6386.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Makuch-Pietraś, I.; Wójcikowska-Kapusta, A. Differences in the Content of Zn Fractions in the Profiles of Soils from Allotment and Domestic Gardens in South-Eastern Poland. Land 2021, 10, 886. https://doi.org/10.3390/land10090886

AMA Style

Makuch-Pietraś I, Wójcikowska-Kapusta A. Differences in the Content of Zn Fractions in the Profiles of Soils from Allotment and Domestic Gardens in South-Eastern Poland. Land. 2021; 10(9):886. https://doi.org/10.3390/land10090886

Chicago/Turabian Style

Makuch-Pietraś, Iwona, and Anna Wójcikowska-Kapusta. 2021. "Differences in the Content of Zn Fractions in the Profiles of Soils from Allotment and Domestic Gardens in South-Eastern Poland" Land 10, no. 9: 886. https://doi.org/10.3390/land10090886

APA Style

Makuch-Pietraś, I., & Wójcikowska-Kapusta, A. (2021). Differences in the Content of Zn Fractions in the Profiles of Soils from Allotment and Domestic Gardens in South-Eastern Poland. Land, 10(9), 886. https://doi.org/10.3390/land10090886

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