Effects of afforestation of agricultural land with grey alder(Alnus incana(L.)Moench)on soil chemical properties,comparing two contrasting soil groups

O??erts Nikodemus ,Dārta Kaupe,Imants Kuku?s,Guntis Brūmelis,Raimonds Kasparinskis,Iluta Dau?kane and Agita Treimane

Abstract

Keywords: Afforestation, Grey alder, Soil properties,Umbrisols, Arenosols, Topsoil

Background

Abandonment of agricultural land and gradual overgrowing with trees and shrubs are common processes across Europe, and were particularly rapid in Eastern Europe after the collapse of the Soviet system (Prishchepov et al. 2012; Estel et al. 2015). Abandonment of agricultural practice and establishment of forest cause changes in the environment, including changes in soil.Data on changes of soil properties in some cases are inconsistent or incomplete, and thus create a need for further research on effect of land use change on ecosystem development. In some studies an increase in soil organic carbon (SOC) stock was observed after afforestation of agricultural lands (Del Galdo et al. 2003;Hooker and Compton 2003; Mao et al. 2010; Armolaitis et al. 2013). However, there have also been cases where the SOC in soil decreased during the first years of afforestation (Ross et al. 1999; Paul et al. 2002; Chen et al.2004). In the advanced stages of forest development,when SOC begins to accumulate in the litter horizon,the total SOC content in soil increases (Kuku?s et al.2015). During forest development, the increase in SOC content is accompanied by an increase in soil nitrogen content (Li et al. 2012; Holubík et al. 2014). Other studies have reported a decrease in magnesium (Mg), calcium (Ca) and potassium (K) concentration in soil,which may affect soil fertility and thus productivity of tree stands (Berthrong et al. 2009). There are inconsistent results on changes in phosphorous (P) levels after land transformation (MacDonald et al. 2012; Deng et al.2016). After afforestation, the content and behaviour of P forms in soil differ within former agricultural lands,depending on soil properties. Changes in labile P content after afforestation were reported to depend on soil characteristics like clay content and pH, while changes in total P were not observed (MacDonald et al. 2012).Changes in SOC and nutrient element content in soil during forest development depend on other factors like previous land use (grasslands, herbaceous plants, etc.),soil tillage and cultivation, tree species, tree stand age,relief, climate and soil properties like clay content and pH value (Paul et al. 2002; Jandl et al. 2007; Laganière et al. 2010; Uri et al. 2014; Cukor et al. 2017). Afforestation of agricultural land with deciduous trees like alder(Alnus) species is common in temperate (Anthelem et al. 2001) and hemi-boreal biomes of Europe (Uri et al.2011). It is considered that afforestation of abandoned lands with grey alder has a potential to increase biomass supply and carbon (C) capture in northern Europe. It is well known that grey alder (Alnus incana (L.) Moench)in symbiosis with Frankia group Actinomycetes captures nitrogen from the atmosphere and litter deposition causes accumulation of nitrogen (N) in soil, along with increase of organic carbon content (Bārdule and Lazdi?? 2010; Uri et al. 2011, 2014). Soil conditions can also affect accumulation of C, N and other chemical elements. Therefore, the aim of the study was to determine changes in soil properties during the natural afforestation of abandoned agricultural lands with grey alder,depending on soil group. We paid attention to development of spatial heterogeneity of soil properties during the afforestation from concentric tree patches.

Materials and methods

Characteristics of study site

The study was conducted in the central part of Latvia, in Vidzeme upland (Fig. 1). Mean annual precipitation in Latvia is 750-800 mm, of which about 500 mm falls in warm seasons. Mean temperature during the warm period is 16.9°C and 5.4°C during the cold period.Mean snow depth is 22-26 cm (Briede and Kore?ska 2018).The study territory is located in the hemiboreal forest zone (Ahti et al. 1968). The landscape of the study territory is a mosaic of agricultural land, forests dominated by spruce, pine, birch and aspen, and former agricultural land mostly overgrown by grey alder and spruce (Ruskule et al. 2016).

Arenosols, Podzols, Cambisols, Retisols, Luvisols, Stagnosols, Gleysols, Umbrisols and Histosols are the most dominant soils in the studied region, which formed on glacigenic (loam, sandy loam and sand), glaciofluvial(sand, coarse sand) and peat deposits (Kasparinskis et al.2017). Afforestation with grey alder was studied in two sites. In one, Dystric Arenosols formed on glacigenic deposits of loamy sand in the soil parent material and sand in the topsoil. In the other, afforestation occurred on Endostagnic Umbrisols formed by loamy sand in topsoil,and loam, silt loam and sandy loam in soil parent material. The depth of the Ap horizon varied from 24 to 40 cm in the Arenosols, and from 30 to 50 cm in the Umbrisols. Agricultural practice in both sites was discontinued in the mid-1990s. Historically these lands were managed as meadows and croplands.

Study methods

Three patches of grey alder were selected in each territory delineated by a soil group. We selected tree patches that contained only grey alder. The mean area of grey alder patches was 81 m2. Lists of plants species were made within and just outside (2-3 m) the selected tree patches prior to soil sampling. Age of the grey alder trees was determined from tree cores obtained with a Presler’s auger.

Soil textural classes were determined according to FAO Guidelines for Soil Description (Jahn et al. 2006).Soil profiles were described and soil groups were determined according to the international WRB soil classification (IUSS Working Group WRB 2015).

Soil samples for chemical analysis were collected at the centre of patches, and at distances of 2 and 4 m in north, east, south and west directions from the centre to the edges of the patches (Fig. 2). The patch edges were only recently colonised and this zone contained only a few small alder trees. Control samples were collected in the grassland 1-2 m outside the tree crown and alder rooting zone, at a distance about 3-4 m from the sample located at 4 m from the patch centre. Soil samples were collected from small pits with size 20 cm×20 cm at depths 0-10 cm and 11-20 cm. In each grey alder patch,the soil profile was described in pits. A total of 32 soil samples were collected from each alder patch.

Three replicates were taken from each collected soil sample and analysed in the soil laboratory of the University of Latvia Faculty of Geography and Earth Sciences.Air dried samples were broken up and sieved through a 2-mm mesh. Particle sizes were determined by pipette analysis (van Reeuwijk 1995; Cools and De Vos 2010),pH value was measured with an “Adrona pH meter AM1605” in 1 mol·L-1BaCl2solution (Tan 2005; Cools and De Vos 2010),content of total carbon(Ctot) and total nitrogen (Ntot) was determined by an “EuroVector EuroEA”analyser(Tan 2005;Cools and De Vos 2010)and content of exchangeable cations (Ca2+, Mg2+, K+, Fe3+,Al3+)in 1 mol·L-1BaCl2solution by an atomic absorption spectrometer “Perkin Elmer Analyst 200” (Tan 2005;Cools and De Vos 2010). Available phosphorous (Pl) content was determined by Mehlich 3 extraction method with an “Thermo Scientific iCAP7000 ICP-OES” spectrometer(Gutiérrez Boem et al.2011;AgroEcoLab 2016).

Data analysis

The Shapiro-Wilk test was used to test for normal distributions. For statistical analysis, data on soil physical and chemical properties not conforming to normal distribution (pHBaCl2, Ntot, K+, Fe3+and Al3+) underwent lognormal distribution transformation. Firstly, main and interaction effects of soil, horizon and location under the canopy were tested by 3-way ANOVA.Then, statistically significant (p ＜0.05) differences between soil types for each soil horizon in the grassland (control) sites were determined by Tukey’s one-way ANOVA. A one-way ANOVA with Tukey’s post hoc test was used to test for significant differences in soil properties between locations for each soil type and horizon. SPSS PASW Statistics 18 software was used for the analysis.

Results

Tree patches

Tree age in the centre of the patches was up to 20 to 25 years and the youngest trees (4 to 10 years old) were at the edges, demonstrating concentric vegetative spread around trees that colonised by seed. Average height of trees in the centre of the patches was 13.2±0.5 m on sandy Dystric Arenosols and 12.2±0.6 m on sandy loam Endostagnic Umbrisols. Average diameter of these trees was 16.2±0.7 cm and 14.8±0.7 cm, respectively. Height and diameter of the trees gradually decreased from centre to the edges of patches. Vegetation in the study site on Arenosols was dominated by Elytrigia repens,

Dactylis glomerata and Festuca ovina, and on Umbrisols soil by Elytrigia repens, Aegopodium podagraria, Vicia cracca and other plants (Table 1).

Descriptive statistics

Three-way ANOVA indicated significant main effects of soil, horizon and location on chemical properties, except for soil on Ntotand location on Ca2+and pH (Table 2).Significant interaction effects of soil horizon occurred for Ntot, K+, Ca2+, Al3+and pH, and of soil location for Ntot, Ctot, P, Mg2+, K+and Al3+. For the grassland sampling sites, which can indicate soil properties before treecolonisation, there were no significant differences (1-way ANOVA) in Ntotcontent between Dystric Arenosols and Endostagnic Umbrisols (Table 3). Sandy Arenosols had higher pHBaCl2and Mg2+and Fe3+and Al3+content,while Ctotand P content was higher in sandy loam Endostagnic Umbrisol soil.

Table 1 Characterisation of the dominant vegetation in the study sites,listed by decreasing cover

Effect of location under the tree canopy

One-way ANOVA for each soil type and soil layer showed that soil pHBaCl2, Ca2+and Mg2+content after 25 years of afforestation. Ntotcontent in topsoil (1-10 cm layer) was significantly higher (by 29.4%) in the centre of the grey alder patches on sandy Dystric Arenosol soil than in the periphery of wood patches and in grassland (Fig. 3a). In the deeper 11-20 cm soil layer,the difference in Ntotcontent was greater, being 60%higher in the patch centre than in the grassland. In loamy clay Endostagnic Umbrisol soil, the Ntotcontent did not significantly differ between sampling sites (Fig.3a). In all studied sites Ntotcontent in the top layer of the Ap horizon (0-10 cm) was higher than in the deeper layer (11-20 cm).

Ctotcontent in the 0-10 cm and 11-20 cm layers of the sandy Dystric Arenosols was significantly higher at the centre of tree patches than in sampling sites closer to the edges and grasslands (Fig. 3b). Afforestation with grey alder increased Ctotcontent by 37.9% in the 0-10 cm and by 88.1% in the 11-20 cm layer. However, Ctotcontent did not significantly differ between sampling locations in the sandy loam Endostagnic Umbrisols (Fig.3b). Both of the studied soils had higher Ctotcontent in the 0-10 cm layer compared to the 11-20 cm layer.

Soil in the centre of tree patches in sandy Dystric Arenosols also had significantly higher K+content than in the periphery and grassland (Fig. 3c). K+content of soil in the centre of wood patch was two times higher than in the grassland. This spatial pattern was not observed for the Endostagnic Umbrisols. Mean content of P in the 0-10 cm layer in the centre of patches was higher than in the grassland for both soils (Fig. 3d), but this difference was significant (by 84.9%) only for sandy Dystric Arenosols.

The only significant difference in Fe3+concentration between sampling sites was observed in Endostagnic Umbrisols in the 0-10 cm layer, where Fe3+concentration in the centre of wood patch was higher by about 140% than in the other sites (Fig. 4a). Aluminium was the only studied element that had higher content in the deeper layer (11-20 cm) than in the top layer (0-10 cm)(Fig. 4b). In Endostagnic Umbrisols, the Al3+concentration in the 11-20 cm layer tended to be higher toward the patch centre (Fig. 4b).

Discussion

The observed interaction of soil and horizon effects was not surprising considering impact of soil texture on soil chemistry and processes (Vanmechelen et al. 1997).However, the significant interaction effect of soil andlocation under the tree canopy clearly shows the importance of soil type in understanding of the role of grey alder in afforestation. The space-for time study clearly showed impact of afforestation by grey alder on soil properties, with the effect being stronger in the centre of tree patches where colonisation was initiated, and weaker toward the edges where trees were younger.

Table 2 Three-way ANOVA effects of soil,horizon and location under grey alder canopy on soil chemistry

Table 3 Mean values of chemical properties of sandy Dystric Arenosols and sandy loam Endostagnic Umbrisols in the grassland sites

Sandy Dystric Arenosols, which are relatively poor in nutrients (P and K+) and organic matter, accumulated carbon and nutrient elements (nitrogen, potassium,phosphorous) during the afforestation of grasslands at a higher rate than in the more fertile Endostagnic Umbrisols. Many other studies have reported the role of grey alder afforestation in soil enrichment with nitrogen (Uri et al. 2014; Innangi et al. 2017; Morozov et al. 2018) and carbon (Uri et al. 2017). The changes in topsoil (0-10 cm) are caused by the accumulation of alder leaf litter and formation of easy decomposable and nitrogen-rich biomass (Innangi et al. 2017). The development of rootnodule symbioses with Frankia in roots of alder (Huss-Danell and Lundmark 1988) is affected by soil properties like pH value and nitrogen content. Increased nitrogen content in soil can affect nodulation of Frankia causing decreased nitrogen accumulation (Bond et al. 1954). In our case, there were no differences between the studied soils in N content in grassland, suggesting that this factor did not cause differences between the soil types in N fixation and uptake with plant roots. In contrast to the Dystric Arenosols soils, Endostagnic Umbrisols soil in grassland had significantly higher Ctot, P and K+content.Experimental studies have shown that nitrogen fixation by grey alder can be affected by phosphorous and potassium content in soil (Ekblad and Huss-Danell 1995). Low soil P availability limits N2fixation in alder stands (Uliassi and Ruess 2002), while high P content in soil has an positive effect on nodulation and N2-fixation rate (Ekblad and Huss-Danell 1995; Uliassi and Ruess 2002; Gentili and Huss-Danell 2003). In our study more rapid accumulation of Ntotand available P in soil after afforestation of agricultural lands with grey alder occurred in the soil type with lower phosphorus content(Fig. 3a and d).

Part of the nitrogen released from alder leaf litter may have been taken up by the herb layer plants under the tree canopy. This might have caused part of the differences in N accumulation in soil between the soil groups.The herb layer was visually denser on the richer Endostagnic Umbrisols, but we did not investigate nitrogen uptake by herb layer plants. Previously, it has been reported that nitrogen accumulation by plants can cause a decrease in soil nitrogen content (Mileti? et al. 2012;Podwika et al. 2018)after afforestation.

Our results showed that leaf litter accumulation not only had significant effect on N accumulation in soil, but also on K+and available P content in the upper layer of soil (Fig. 3c and d). The decreased Kexchcontent in the centre of alder patches in the Endostagnic Umbrisols soils might be due to herb layer plant uptake. The observed P accumulation in soil after afforestation with grey alder has also been reported in other studies (Uri et al.2003;Taleshi et al.2009),which can be explained by microbiological and mycorrhizal activity (Ingestad 1987;Giardina et al. 1995; Arveby and Granhall 1998). Alder has also been reported to increase availability of phosphorus via root secretion of phosphate enzymes(Giardina et al. 1995). However, the increase in available P content was statistically significant only in the sandy Dystric Arenosols and not in the Endostagnic Umbrisols.

The increased Al3+concentration in the centre of grey alder patches in the deeper horizon of Endostagnic Umbrisols might be due to root exudates that indirectly affect soil reaction by solubilising unavailable aluminium(Skyllberg 1996). However, in our study a significant change in pH value was not observed in topsoil after afforestation with grey alder. The literature contains contrasting reports on effect of alder afforestation on soil pH value. In some cases afforestation with grey alder has been observed to cause soil acidification (Uri et al. 2014;Podwika et al. 2018), while other studies showed an increase of soil reaction and exchangeable base content(Huss-Danell and Lundmark 1988).

During the afforestation of grasslands with grey alder high spatial heterogeneity in soil properties was found,indicated by significant differences between tree patches and grassland, and also within tree patches. However,with creation of a closed cover of alder over time, the soil properties in grassland can be expected to converge.Formation of soil properties after afforestation is also affected by soil bacterial communities (Preem et al. 2012),and the developmental trajectories of these may differ depending on soil properties. Further studies are needed in soil ecology to evaluate effect of tree roots on soil physical properties, moisture conditions, biochemical and other processes.

Conclusions

The results showed increased content of Ctot, Ntot, available P and exchangeable K+concentration in topsoil after afforestation of agricultural lands with grey alder(Alnus incana (L.) Moench) during the time span of 25 years on sandy soils poor in organic matter soils. The effect in more fertile soils with higher organic matter content was less evident and might have been confounded by herb layer plant uptake. Afforestation of agricultural lands with grey alder had effects on physical-chemical and biochemical processes in soils: increased P availability to plants in soil, and changes in Al3+and Fe3+ion mobility. Our results show that initial soil properties in grasslands significantly interact with effects due to afforestation with grey alder, which adds complexity to producing models of biogeochemical processes and creates difficulty in extrapolation of results to other soil types to guide land management.

Mosaic type afforestation was shown to create higher spatial variability of soil properties, which needs to be considered during soil mapping and in choosing the correct mapping scale.

Abbreviations

Ctot: Total carbon; Ntot:Total nitrogen;SOC: Soil organic carbon; Ca2+, Mg2+,K+,Fe3+,Al3+: Exchangeable cations.

Acknowledgments

Not applicable.

Authors’contributions

Study concept and design:Nikodemus O, Kaupe D, Brūmelis G. Analysis and interpretation of data: Kaupe D, Kuku?s I, Dau?kane I,Treimane A. Drafting of the manuscript: Nikodemus O, Kaupe D. Critical revision of the manuscript for important intellectual content: Kuku?s I, Brūmelis G, Kasparinskis R.Statistical analysis: Kaupe D, Kasparinskis R. Obtained funding: Nikodemus O.Study supervision: Nikodemus O. The author(s)read and approved the final manuscript.

Funding

The study was supported by University of Latvia grant Nr. AAp2016/B041//Zd2016/AZ03 within the project “Climate change and sustainable use of natural resources”.

Availability of data and materials

Research data are not shared.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that there are no competing interests.

Author details

1Faculty of Geography and Earth Sciences, University of Latvia, Jelgavas iela 1,Riga LV-1004, Latvia.2Faculty of Biology, University of Latvia, Jelgavas iela 1,Riga LV-1004, Latvia.

Received: 20 December 2019 Accepted: 9 June 2020