Next Article in Journal
Ethnobotany of the Himalayas: Safeguarding Medical Practices and Traditional Uses of Kashmir Regions
Next Article in Special Issue
Enhancement of the Aroma Compound 2-Acetyl-1-pyrroline in Thai Jasmine Rice (Oryza sativa) by Rhizobacteria under Salt Stress
Previous Article in Journal
Vincristine in Combination Therapy of Cancer: Emerging Trends in Clinics
Previous Article in Special Issue
The Role of Plant Growth-Promoting Bacteria in Alleviating the Adverse Effects of Drought on Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contrasting Effects of Forest Type and Stand Age on Soil Microbial Activities: An Analysis of Local Scale Variability

by
Anna Walkiewicz
1,*,
Andrzej Bieganowski
1,
Adrianna Rafalska
1,
Mohammad I. Khalil
2,3 and
Bruce Osborne
2
1
Institute of Agrophysics, Polish Academy of Sciences, Doświadczalna 4, 20-290 Lublin, Poland
2
UCD School of Biology and Environmental Science and UCD Earth Institute, University College Dublin, Belfield, 4 Dublin, Ireland
3
School of Applied Sciences and Technology, Prudence College Dublin, 22 Dublin, Ireland
*
Author to whom correspondence should be addressed.
Submission received: 13 July 2021 / Revised: 25 August 2021 / Accepted: 26 August 2021 / Published: 31 August 2021

Abstract

:

Simple Summary

The potentially important role of forests in climate change mitigation suggests a strong need for a more detailed understanding of these ecosystems. Besides climatic conditions, diverse forest vegetation creates varied conditions for the activity of soil microorganisms, and particular attention should be focused on a comprehensive study on the influence of different forest types on microbial activities. We conducted an experiment on six different forest soils (two coniferous, two deciduous, and two mixed sites comprising trees of different ages) collected from the same region (Lublin Upland, Poland) to assess the relationship between forest type and seasonal changes in microbial parameters. The annual mean values of the soil microbial indicators suggest that the mature deciduous stand was the most sustainable in microbial activities among the forest soils investigated. The diversity of the forest environment and the multifactorial dependence of the microbiological activity of forest soils necessitates further research in this field, especially using the same soil types. An understanding of forest ecosystem functioning can also be useful for forest management.

Abstract

Understanding the functioning of different forest ecosystems is important due to their key role in strategies for climate change mitigation, especially through soil C sequestration. In controlled laboratory conditions, we conducted a preliminary study on six different forest soils (two coniferous, two deciduous, and two mixed sites comprising trees of different ages) collected from the same region. The aim was to explore any differences and assess seasonal changes in soil microbial parameters (basal respiration BR, microbial biomass Cmic, metabolic quotient qCO2, dehydrogenase activity DHA, and Cmic:Corg ratio). Indicator- and forest-specific seasonality was assessed. In addition to litter input, soil parameters (pH, nutrient content, texture and moisture) strongly regulated the analyzed microbial indicators. PCA analysis indicated similarity between mature mixed and deciduous forests. Among annual mean values, high Cmic and DHA with simultaneously low qCO2 suggest that the mature deciduous stand was the most sustainable in microbial activities among the investigated forest soils. Research on the interrelationship between soil parameters and forest types with different tree ages needs to be continued and extended to analyze a greater number of forest and soil types.

1. Introduction

An understanding of the microbial activity in forest soils is important given its role in the maintenance of the biodiversity of these ecosystems whilst also contributing to carbon (C) storage and greenhouse gas (CO2, CH4, N2O) fluxes [1,2,3]. Consequently, forests with appropriate management practices are often an important part of global biodiversity maintenance and underpin many strategies for climate change mitigation. Globally, forests cover 31% of all land surface [2], including 35% of the area of Europe, with a coverage of about 227 million ha in Europe [4]. Considering their types, coniferous stands dominate in Europe (46%), followed by broadleaved (37%) and mixed (17%) stands. In terms of the age of European forests, sites between 20 and 80 years dominate (accounting for 40%), while 12% and 18% represent sites below 20 years and over 80 years, respectively [5]. Forests grow in a range of climatic conditions and soil types [6], which is associated with various levels of biodiversity [2].
The soil ecosystem determines the growth environment for plants and forms a reservoir of available water and nutrients for trees, whose roots connect with soil microbial communities [7,8]. The activity of soil microbes reflects the functioning and condition of the forest ecosystem, since soil microorganisms are highly sensitive to changes and react faster to disturbances than many multicellular organisms [9,10]. Soil quality can be assessed using various microbial indicators [11]. Basal respiration (BR) is one of the oldest parameters, reflecting the biological activity of soil microorganisms and C availability for microbial growth. This includes the catabolic processes of aerobic microorganisms associated with organic matter decomposition and nutrient release or immobilization in the soil [1,10,12,13]. Soil microbial biomass (Cmic) is an important ecological indicator of soil fertility and quality [14] and refers to the soil fraction that is responsible for the regulation of the transformation of organic matter and the cycling of nutrients and energy [10]. The metabolic quotient (qCO2) is a sensitive indicator of the likely constraints on the soil microbial community and can be used to assess ecosystem development and soil disturbance [15,16]. A low qCO2 value indicates that a greater part of C is built into the biomass and a smaller part is emitted as CO2. Among the soil microbial enzymes, the activity of intracellular dehydrogenase (DHA) is a bio-indicator of soil quality, health, and degradation and provides information about overall microbial activity in soils, because the enzyme is linked to viable cells [10,17,18,19]. Soil dehydrogenases act as carriers of electrons and hydrogen from the substrate to the acceptors during the early stages of soil organic matter oxidation; thus, the activity of these enzymes is also considered an indicator of oxidative metabolism in soils [18,19]. The Cmic:Corg ratio can be used as a measure of the deterioration of soil quality [20], which determines the availability of soil C. In turn, C is necessary for growth and microbial activity [16] and is closely associated with the capacity of the soil to support populations of microorganisms [21]. Investigations that take into account this set of contrasting parameters allow a comprehensive assessment of soil (in our case—forest soil) quality.
Soil microbial biomass and activity in soils of different forest types show seasonal dynamics and are regulated by changes in climate, soil chemical and physical properties—mainly temperature, moisture, texture, pH, and nutrient availability [19,20,22,23,24,25,26,27,28,29,30]. The seasonality of BR may be highly dependent on the dynamics of soil temperature variability [31,32,33]. High microbial biomass was reported in pine forests in spring [34] and in a fir stand in summer [35]. A high Cmic in autumn may be an effect of a large amount of litter, which provides a source of nutrients and energy for soil microorganisms [36]. The metabolic quotient (qCO2) was higher in a cedar forest during the winter season and lower in a beech forest in the spring season [37], which demonstrates that changes in C use efficiency can occur during the year. Soil DHA is stimulated by a rise in temperature [29], with the decomposition of litter [38] reaching its lowest values in winter. Among the different physical parameters, soil texture influences microorganisms, since it regulates hydrology, aeration, and soil organic matter. Clay soils have a higher Cmic and DHA than sandy soils [7,23,26,27,39].
A specific aspect of the forest ecosystem is the regulation of soil microbial activity by trees of different ages representing different species, which produce diverse litters and root exudates and thus change soil hydrological and physical properties [40,41,42,43,44,45,46,47]. The effect of trees on soil microbial communities was observed in unfavorable site conditions with a low nutrient content [8]. It was reported that site conditions and tree species are very important factors defining the structure of the soil microbial community [48]. As well as the underlying soil type, differences in soil respiration may be an effect of the influence of different tree stands on the soil specific and on soil microorganisms [49]. Through the different ways that trees of different ages regulate the soil environment (pH, litter quality), this also results in changes in the enzymatic activity of soil [17,23,43]. Soil microbial parameters may also be significantly affected by nutrient (C, N, and P) content regulated by trees of different age; for instance, young trees need more N but less P than old trees [44,46]. Soil biota participation in C cycle results in their strong relationship with C sequestration and soil-atmosphere GHG exchange. The microbial contribution to C sequestration is governed by the interactions between microbial (the amount of biomass, the community structure, and microbial by-products) and soil properties (texture, clay mineralogy, and pore-size distribution) [50]. In addition to CO2 emission through the respiration process, the activity of different soil microorganisms, such as methanotrophs or nitrifiers, causes CH4 and N2O uptake.
To conduct preliminary research and better characterize the microbial populations associated with soils of diverse forests and vegetation, a comprehensive approach was taken in this study encompassing soils under different species and age of trees located in the same region (under similar weather conditions, Lublin Upland in Poland). We focused on a set of biochemical parameters related to microbial biomass/activity and C-use efficiency, that, in combination, allowed for a better interpretation of their role(s) in contrasting forest soils at the local scale. We hypothesized that interactions between climatic and soil parameters, forest type, and tree age significantly differentiate soil microbial indicators. We assumed higher values of soil microbiological biomass and activity and lower C use efficiency may be determined by:
(1)
Higher temperatures in summer and spring;
(2)
Higher water availability in spring than in autumn and winter (while drought conditions in summer may reduce microbial activity);
(3)
Higher water storage and nutrient (N, P, C) availability in heavy vs. light soils;
(4)
Lower water and higher nutrient demand of younger trees;
(5)
The presence of leaf litter, which is most abundant in mature deciduous forests, and which could maintain a higher soil temperature in winter than in sites without litter and protect against excessive water loss in summer.
The aim of the study was to determine and utilize the microbial parameters: basal respiration, soil microbial biomass, metabolic quotient, dehydrogenase activity, and the Cmic:Corg ratio for characterizing microbial activity and C use efficiency in six different forest soils in the same region (two coniferous, two deciduous, and two mixed sites) under trees of different age.

2. Materials and Methods

2.1. Site and Soil Characteristics

Soil samples were collected from six different forest sites located in the Lublin Upland, Poland. The basic characteristic of the forests examined are presented in Table 1.
During the sampling period (from July 2018 to April 2019), the annual mean temperature in the study area was 8.5 °C and the total annual rainfall was 330 mm. In the studied region, the highest mean temperature reached 20.7 °C in summer, while the lowest value (−0.4 °C) was recorded in winter. Based on the regional long-term data (2013–2018), the mean air temperature was 8.5 °C and mean total precipitation was 556 mm. Data for individual months are shown in Figure 1. The average temperature of the soils in situ were: 19.3 °C in summer, 9.9 °C in autumn, 2.5 °C in winter and 9.5 °C in spring.
Basic soil properties (Table 2) were determined for samples collected from a depth of 0–15 cm (n = 3) in July 2018 before the main experiment.
Soil texture was measured by the laser diffractometer Mastersizer 2000 with the Hydro G dispersion unit (Malvern Ltd., Malvern, UK). The following settings were used: two light sources of laser (633 nm) and diode (466 nm), stirrer speed 700 rpm, and pump speed 1750 rpm, soil refractive index 1.52 and soil absorption index 0.1 [51]. Soil pH was determined potentiometrically at room temperature in a soil and water slurry, after the soil had settled, with a ratio of soil:water of 1:2.5 w/w. The total C and N contents were determined in soils ground in a mortar using a dry combustion method (an oven temperature of 1020 °C; a Thermo Scientific Flash 2000 Organic Elemental Analyzer). Soil inorganic (SIC) and organic carbon (SOC) contents were determined with a TOC-VCPH analyzer (Shimadzu, Kyoto, Japan). Phosphorus (total P) concentration was determined by ICP-OES (Inductively Coupled Plasma Optical Emission Spectrometry) using a Thermo Scientific iCAP Series 6500 equipped with a charge injection device (CID) detector [52]. Soil temperature was measured during sampling using a TDR soil multimeter FOM/mts (n = 8). Undisturbed soil cores (100 cm3 brass cylinders) were tested after equilibration on a Richard chamber (Soilmoisture Equipment Corp., USA) in order to determine soil bulk density as the mass of a dry sample divided by the volume of the sample.

2.2. Experimental Design

The study was conducted on soil samples collected at monthly intervals in summer 2018 (July, August), autumn 2018 (September, October, November), winter 2018/2019 (January, February), and spring 2019 (April). The samples (surface layer, i.e., 0–15 cm depth after removal of litter) were collected at five randomly selected sites in each forest site, mixed, and thoroughly homogenized, and combined into one representative sample for the season and for each forest type. Representative samples were sieved through a 2 mm mesh, frozen (−20 °C) and next stored in the dark at 4 °C for four weeks before the microbial analyses to minimize microbial activity. To represent field conditions, the analyses were conducted at soil moisture levels corresponding to the natural conditions at the time of sampling (Table A1). The measurement procedures were performed after a two-day preincubation to stimulate microorganism activity (5 g of air-dry soils was placed in 60 cm3 glass vessels and incubated in the dark at 25 °C).

2.3. Analysis of Microbial Parameters

Soil basal respiration (BR, expressed in µg CO2-C/g/h) was measured after 2 h incubation at 25 °C. Soil microbial biomass (Cmic) was determined using the substrate-induced respiration (SIR) method. Emission of CO2 was measured after enrichment of the soil with a glucose solution (10 mg per gram of soil) that provided an easily available source of carbon and energy [53]. After 2 h incubation with shaking at 25 °C in a water bath, the CO2 produced was collected by injection and measured chromatographically using a Shimadzu GC-14A (Shimadzu Corp., Kyoto, Japan). Microbial biomass content was recalculated according to Šimek and Kalčík, 1998 [54]. The activity of soil dehydrogenases (DHA) was determined with the Casida et al., 1964 [55] method using triphenyl tetrazolium chloride (TTC); this analysis was based on the amount of triphenyl formazan (TPF) produced after a 20 h incubation of the soil samples at 30 °C. After filtration, DHA (mg TPF/g/20 h) was measured spectrophotometrically at an absorbance of 485 nm (UV-1601PC, Shimadzu Corp., Kyoto, Japan). In total, 216 soil samples were used in this study to assess the microbial indicators. Based on the measured parameters, the Cmic:Corg ratio and BR:Cmic ratio (the metabolic quotient qCO2 expressed as µg CO2-C/mg Cmic/h) were calculated [56].

2.4. Calculation and Statistical Analysis

The values of the microbiological parameters were calculated for each season and each forest type. Based on the seasonal indicators, the average annual values of BR, Cmic, qCO2, DHA, and Cmic:Corg were calculated for each study site. The value of BR and qCO2 was corrected for CO2 dissolution in soil solution according to the Henderson–Hasselbach equation.
The results were statistically analyzed with Statistica 13 software (StatSoft Inc., Tulsa, OK, USA). A non-parametric Kruskall–Wallis test was used to evaluate the significance (at the 5% level) of the differences in soil parameters between the seasons in each forest and between the annual values for each forest separately. Principal component analyses (PCA) were carried out to examine the relationships among soil microbial indicators (BR, Cmic, qCO2, DHA, Cmic:Corg), tree age, and soil parameters (sand, silt, and clay contents, moisture, pH, N, C, and P concentration, C:N, N:P, and C:P).

3. Results

The studied microbiological parameters of the forest soils showed different annual mean values (Figure 2). The preliminary results of seasonal variability of the tested microbial parameters are shown in Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7.

3.1. Annual Mean of Microbiological Parameters of Forest Soils

The annual mean values of each analyzed microbiological parameter are shown in Figure 2.
The highest BR was recorded in the soil from the young coniferous forests. It was about 20% lower for the rest of the studied forest soils. The mean annual Cmic had significantly higher values for the silty soils collected from the young coniferous forest and from the mature deciduous forest. It was about 40% lower for the sandy soil and the mature coniferous forest. The mean qCO2 value in silty soil from deciduous mature forest was significantly lower than in other stands, and the highest values were found in the sandy soils from the mature coniferous forest. The annual soil DHA values decreased in the following order: young deciduous > mature deciduous > mature and young coniferous and mixed. The highest annual Cmic:Corg ratio was found for the silty soil of the young coniferous and middle-aged mixed forest, and values for these mature forests were approximately 50% lower.

3.2. Seasonal Changes in Soil Basal Respiration (BR)

Seasonal changes in BR were observed in all studied forest soils (Figure 3), although the pattern varied and differences between seasons were not significant.
In the mature coniferous forest on the sandy soil, BR was about 40% higher in winter compared to the other seasons. A higher BR in winter was also observed in the young coniferous forest on the silty soil; a similar level of BR was observed in autumn, but lower values were recorded in summer and spring. For the sandy soil from the middle-aged mixed forest, a higher BR was recorded in samples collected in summer; BR was about 24% lower in the soil collected in autumn and spring. In samples from the mature mixed forest on silty soil, the highest BR was found in winter and the lowest in summer. Silty soils of both deciduous forests had the highest BR value in autumn and the lowest value in spring, which did not differ significantly from the level observed in summer and winter.

3.3. Seasonal Changes in Microbial Biomass (Cmic)

The seasonal changes in microbial biomass (Cmic) are shown in Figure 4.
In all soils collected from the coniferous and mixed forests, the highest soil microbial biomass (Cmic) was observed in winter. The young coniferous forest soil was the only case where the lowest Cmic was found in autumn and was almost two-fold lower than in samples collected in summer, autumn, and winter. The sandy (from mature coniferous and middle-aged mixed stands) and silty soils (from deciduous forests) showed the lowest Cmic in summer. The lowest Cmic in the silty soil from the mature mixed forest was recorded in spring, although this value did not differ significantly from the results obtained in samples from summer and autumn. The highest Cmic was recorded in samples from autumn only in the soils from the deciduous forests. Inconsistency due to the large standard deviation should be considered in the case of young coniferous and mature deciduous stands.

3.4. Seasonal Changes in Soil Metabolic Quotient (qCO2)

The results of the metabolic quotient (qCO2) are presented in Figure 5.
In summer, soils from most forests (mature coniferous, mature and young deciduous, middle-aged mixed) showed the highest metabolic quotient (qCO2). For the sandy soil collected from the mature coniferous forest, qCO2 was also high in autumn, while similar and significantly lower values were found in winter. For the young coniferous forest, the highest qCO2 was recorded in the autumn, and the values for the other seasons were approximately 3 times lower. The soils from the deciduous forests and the mature mixed forest had the lowest qCO2 values throughout all seasons. The lowest qCO2 values were noted in the silty soils collected from these forests in autumn (both deciduous forests), winter (mature deciduous and mature mixed site), and spring (young deciduous site). The greatest (approximately 4-fold) difference in qCO2 values between the highest (in summer) and lowest (in the other seasons) was recorded in the sandy soil from the middle-aged mixed forest.

3.5. Seasonal Changes in Soil Dehydrogenase Activity (DHA)

Significant seasonal changes in DHA occurred in all soils, apart from those collected from the mature coniferous forest (Figure 6).
Values for DHA in the sandy soil of the mature coniferous forest did not differ significantly across all seasons, although the highest numbers were observed in soil samples collected in spring. The silty soil from the young deciduous forest had particularly low values that were approximately 10 times higher in spring compared to the other seasons. In the soils from both mixed forests, the highest DHA was observed in samples from summer and autumn. The silty soil of the mature deciduous forest was the only soil that reached the highest DHA in summer. The value was three times higher than in the other seasons. The lowest DHA was recorded in the soils from winter in all the studied forests.

3.6. Seasonal Changes in Soil Cmic:Corg Ratio

The seasonal changes in the Cmic:Corg ratio in the forest soils are shown in Figure 7.
In soils collected from both coniferous forests, the highest Cmic:Corg ratios were found for samples collected in winter, reaching values that were approximately two times greater than those from the younger forest. The lowest Cmic:Corg values were found in the sandy soils of the mature coniferous forest collected in summer and in the silty soil of the younger coniferous forest collected in autumn. For the mature mixed forest, the highest Cmic:Corg values were found in winter. For the sandy soil of the middle-aged mixed forest, the highest Cmic:Corg ratios were recorded in autumn and winter and were about 3 times lower than those collected in the summer. The Cmic:Corg ratio in the silty soil of the mature deciduous forest did not differ significantly across all seasons, but the highest value was noted in autumn and the lowest in summer and spring. Similarly, in the silty soil collected from the young deciduous forest, the highest Cmic:Corg ratio was recorded in the autumn samples. It was about half that of the samples collected in summer.

3.7. Principal Component Analysis

Principal component analysis (Figure 8), including a total of five microbial indicators (BR, Cmic, DHA, qCO2, Cmic:Corg) and twelve parameters (tree age, sand, silt, and clay contents, pH, C, N, P content, moisture, C:P, N:P and C:N) in six forests, confirmed the high sensitivity of different microbial parameters to soil characteristics. The analysis generated two principal components. The first and second principal components (PC1 and PC2) explained 64.67% of the total variability of the data set. PC1 was negatively correlated with DHA, Cmic, pH, P, moisture, silt and clay content, N, Corg, and N/P and positively correlated with BR, qCO2, Cmic:Corg, sand content, tree age, C:N, and C:P. PC2 was negatively correlated with Cmic:Corg, DHA, qCO2, BR, sand content, pH, P, moisture and positively correlated with Cmic, tree age, clay and silt content, N, Corg, N:P, C:P, and C:N.
The PCA analysis of forest soils clearly distinguished different forest types and tree age, although it failed to separate the mixed mature forest from the mature deciduous forest (Figure 9).
The mature coniferous and middle-aged mixed forests had positive PC1 scores, whereas the young coniferous and deciduous forests and the mature mixed and deciduous forests had negative PC1 scores. Clusters of samples from the mature forests had mostly positive PC2 scores, while those from the younger forests were found in the negative PC2 fields.

4. Discussion

The potentially important role that forests may have in climate change mitigation, especially through C sequestration, suggests the need to identify the variability and drivers of soil microbial activity in different forest types and stand ages. It is well-known that interactions among vegetation, climate and soil conditions have a combined effect on soil microbes. It was previously reported that forest type may affect soil microbial structure and activity through the presence of different tree species of different ages, and that this controls the soil biota [7,8,22,23]. Diverse stands also produce diverse litters, produce different amounts and type of root exudates, and influence soil hydrology and nutrient content [40,41,42,43,44,45,46,47]. As revealed by the PCA (Figure 9), different stands generally have contrasting impacts on soil microbiology. Interestingly, this analysis indicates that stand age is also an important factor regulating microbial activity. Although research on the same soils and different trees would show more precise differences, this preliminary study on forest soils from the same region (with similar weather conditions) presents results that are valuable and are important for understanding local scale variation.

4.1. Soil Basal Respiration (BR)

Seasonal differences in BR were most clearly visible in the soils collected from the coniferous forests but were not evident in the other forested sites. Unexpectedly, the highest BR was recorded in the winter samples (and in the autumn soil from the young forest), which may be the result of the thawing process on the growth and activity of soil bacteria [57,58]. However, these results should be interpreted taking into account the litter layer and the snow cover that protect the soils from complete freezing. The samples taken were stored at 4 °C before tests and analyzed in 25 °C according to standard procedures. Increasing temperatures may result in the enhanced growth of some bacteria, while during thawing, C and N can be released [59]. Low BR during spring (Figure 3) may be a result of quite low soil temperatures (about 10 °C) [31,32,33], while a high soil BR suggests the rapid decomposition of organic residues and the stimulation of heterotrophic microorganisms in coniferous forests.
A study on different soils in Germany reported that forest type had a significant effect on BR only in nutrient-poor sites, and BR in soils collected from a pure beech forest was 65% higher than in soils under Douglas fir [8]. In our study, soil collected from the young coniferous forest had a high mean annual BR (Figure 2). Malchair and Carnol (2009) [60] also showed higher BR values in soils under two spruce forests, which was explained by low substrate quality resulting in lower C use efficiency, as well as a low quality of litter or high organic matter content.

4.2. Soil Microbial Biomass (Cmic)

The seasonal dynamics of Cmic relates to changes in soil microbial biomass. The lowest Cmic values for most of the analyzed soils were found in the summer samples and were especially low in the soils of the mature coniferous and middle-aged mixed forests (Figure 4). This may be a result of low moisture combined with higher temperatures, which may limit the growth and activity of soil microorganisms [61,62]. A positive correlation between soil moisture and Cmic is shown in Figure A1. Although other studies have shown that the soil microbial biomass may be resistant to drought [62], the high metabolic quotient (qCO2) in the soil collected in the summer from the young mixed forest (Figure 5) may suggest a greater degree of stress than in other seasons in this stand. The highest Cmic values in the deciduous forests were recorded in the soils sampled in autumn (Figure 4); this may be related to the large amount of litter, which provides nutrients and energy for topsoil microorganisms during this season [36]. It is notable that the autumn Cmic in the soil under the mature deciduous forest was about three times higher than in the soil from both coniferous forests, which we assume was a consequence of the thicker litter layer in the deciduous site [34,63,64]. There was also an approximately three-fold larger amount of tree biomass in the mature deciduous forest compared to the coniferous forests (Table 1). Soil microbial biomass may also increase with the diversity of trees [65] and was high in the mature mixed forest analyzed in the present study (Figure 2).
In our study, the soils under the coniferous and mixed forests showed the highest Cmic values in winter, as they may have been protected and more active under litter and snow cover [66]. During freezing-thawing events that could occur locally, microbial cells are destroyed, and damaged cells release nutrients that can be used by the surviving microbes, which can be highly active during soil thawing [67]. A high microbial biomass in the winter season was also observed in pine and oak stands, which was explained by the larger proportion of fungi in the biomass [68,69,70]. Similar to our study, a lower Cmic was observed in summer than in winter, which may suggest that the water availability may have limited Cmic more than temperature [68,69].
The mean annual Cmic decreased in the following order: young coniferous forest > mature deciduous forest > mature mixed forest > young deciduous forest > middle-aged mixed forest > mature coniferous forest (Figure 2). Some authors suggest that the Cmic content is often higher in soils of older forests compared to younger-aged forests [71,72,73], which in our experiment was confirmed only in the soil of the deciduous and mixed stands. In contrast, the Cmic value in the coniferous forests was significantly higher in the soil of the young forest compared to the older one. Texture also regulates Cmic, and in our study, silty soils showed a higher Cmic than the sandy soils (Figure 2). Texture influences microbial biomass since it affects the distribution of micro-, mezo- and macropores that create different microenvironments and determines water holding capacity, nutrient availability and oxygen concentration [7,26,27,39]. Soil with a higher clay fraction can maintain larger microbial communities due to the greater absorption and accumulation of organic C, and the maintenance of a higher water storage capacity [23,34,74,75].

4.3. Microbial Metabolic Quotient (qCO2)

The seasonality of the qCO2 value reflects variations in the bioenergetic status of microbial biomass [76] and the proportion of Corg that can be easily metabolized [77]. In soils from the mature coniferous forest and from the middle-aged mixed forest, the qCO2 value was almost four times higher in summer than in the other seasons (Figure 5). It is known that a higher qCO2 is a result of lower metabolically efficient microbial activity due to the effects of stress or ecosystem disturbances [15,78,79]. One of the unfavorable factors limiting the soil microbial activity may be summer droughts, which have been reported in European forests with increasing frequency [80,81,82,83]. In our study, the lower the moisture, the higher the qCO2 value that was observed (Figure A1). An explanation for the summer sensitivity of the qCO2 of the soils from the mature coniferous forest may be their sandy texture (Table 2), which store less water than silty soils. We also observed a higher qCO2 in the soils collected in autumn from both coniferous forests (Figure 5). Malchair & Carnol, (2009) [60] showed that, in September, the highest qCO2 was found in soils from a spruce forest compared with other forest types (beech, oak, Douglas fir) because of the poor litter quality of coniferous stands. Moreover, a higher qCO2 may also be linked with a lower C use efficiency in young stands [84,85] and a higher BR in a spruce forest [60].
The annual summary shows a clear relationship between the forest stand and the qCO2 value (Figure 2); qCO2 decreased in the following order: mature coniferous > middle aged mixed > young coniferous > young deciduous > mature mixed > mature deciduous. Almost a four-fold higher qCO2 was recorded in topsoil under coniferous compared to mixed forest in a study conducted in Italy by Agnelli et al. (2001) [86]. Higher seasonal changes (Figure 5) and higher annual qCO2 values (Figure 2) in soils from coniferous compared to deciduous stands may be a result of the presence of decomposition-resistant litter [87] associated with waxes, phenolic acids, and lignin contained in coniferous litter [88,89,90].

4.4. Soil Dehydrogenase Activity (DHA)

All preliminarily analyzed soils, except the sample collected from the mature coniferous forest, showed seasonal variations in DHA (but not always significant, Figure 6), reflecting changes in the overall activity of soil microorganisms. These differences may be associated with temperature changes. As reported by Wolińska & Stępniewska, (2012) [29], the highest DHA was observed at 20–30 °C, which may correspond to temperatures characteristic of summer and early autumn in the studied period (Figure 1). Moreover, this temperature range is comparable to the optimal value for the growth and activity of soil microorganisms [91]. DHA can also be stimulated by compounds released during litter decomposition [38], which is more intensive under wet conditions. Based on the meteorological data, we assume that the summer may have created near optimal conditions for DHA activity, since July 2018 was also distinguished by the highest rainfall (Figure 1). As a consequence, soil from the mature deciduous forest had the highest DHA value in this season (Figure 6), in contrast to the soil collected in the summer from the young deciduous stand where there was less aboveground biomass (Table 1) and a low litter level.
Tree species composition has been suggested as one of the main factors influencing soil DHA variations in forests [19], since the presence of vegetation may indirectly affect DHA by increasing fine root biomass in the soil [92] and influence litter amount and quality [93]. Mean annual DHA was significantly higher in the soils under deciduous forests and significantly lower in soils under the coniferous and mixed forests (Figure 2). With the increasing dominance of conifers, the conditions for organic matter decomposition were less favorable and associated with decreasing soil pH [17,23]. Soils under coniferous forests have an acidic pH, which may result in lower DHA compared to soils under deciduous forests [23]. The highest DHA was observed in the soil under the young deciduous forest (Figure 2), which was also characterized by the highest pH of all studied sites (Table 1). Similar to our results, a higher DHA in soils with a higher pH was confirmed in other studies on forest soils [22,23,28]. However, the impact of soil pH may be determined by the range of values to which the microbes are exposed. In our study, the soil pH varied from 4.25–5.37, while the optimum pH for DHA was reported to be about 7 [94,95,96].

4.5. Soil Cmic:Corg

The studied soils showed seasonal differences in Cmic:Corg, although the differences were not always significant between all seasons (Figure 7). High Cmic:Corg indicates a favorable environment for microbial growth, whereas a low ratio may be related to poor quality of the soil organic matter [24]. In our study, higher qCO2 values were observed under higher soil moisture (Figure A1). In the summer, the Cmic:Corg values were at similar levels in most soils, although higher values were found in the soil collected from the young coniferous forest. This may suggest that different forests may create similar soil conditions during the summer season. In our study, winter was the season with the highest Cmic:Corg ratio in the coniferous and mixed forests, which could be a result of a higher soil moisture despite the lower temperatures [97].
The annual summary shows that the highest Cmic:Corg in our study was found in the younger stands of the coniferous and mixed forests (Figure 2), which confirms that these stands created favorable conditions for microbial growth. Other studies have shown contrasting values for the soil Cmic:Corg ratio in different forest types, with the lowest value in coniferous forests and the highest value in a deciduous forest [98] that reflects the contribution of microbial biomass to soil organic C [99].
Differences in the allocation of C, N, and P (Table 2, Figure 10) in the different forest types may also affect the microbial parameters. The uptake of C, N, and P normally constrains the growth of a range of organisms. Uptake from the environment is a key function, as these are the most important elements building living organisms [100], and the ratios of C, N, and P are key indicators of the nutritional status of soils [25].
Tree age may be also connected with soil nutrient ratios [40,41,42,44,45,46,47]. The contribution of N to nutrient content was lower in the younger forests in comparison to the mature stands (Table 2, Figure 10), which may be a consequence of the higher N demand of young compared to older trees [44]. An opposite trend is observed in the case of P, and soils associated with young forests may contain more P than old sites (Table 2, Figure 10), which may be related to the increase in P uptake with stand age [46]. European forest soils were reported to have C:N ratios below 20, with a mean C/N ratio of 16.5 and 15.3 for Cambisols and Luvisols, respectively [40]. These values were close to the mean C:N ratios in our soils (Table 2).

5. Conclusions

We concluded that the presented preliminary studies on forest soils from the same region showed seasonal and annual diversification of microbial parameters. Among the annual mean values, a high Cmic and DHA with a simultaneously low qCO2 suggests that the deciduous mature stand is the most sustainable in terms of microbial activities of the investigated forest soils. Based on the microbial parameters, we may predict that this ecosystem may be the most efficient in C use, considering the contribution of soil microbes in mitigating GHG emissions and enhancement of C sequestration, since it creates favorable conditions for microbiological activity. We revealed similarities between mature mixed and deciduous forests, which may indicate similar conditions for microbial activity. The soils may be more active in summer, autumn and spring than in winter, as confirmed by enzymatic (DHA) activity, as an effect of higher temperatures and litter, especially in soils from the mature deciduous forest and from both mixed forests. Since litter provides nutrients for soil microorganisms in the surface layer, we suggest that this may explain the highest Cmic in the soils collected in autumn from the deciduous forests. A low Cmic and high qCO2 were usually recorded in summer, which suggests that the conditions during this period limited the growth of soil microorganisms, probably due to low moisture. The results are valuable and are important for understanding local scale variation. The diversity of the forest environment and the multifactorial dependence of the microbiological activity of forest soils warrant further research in this area, especially in soil-tree interactions and their relationship with C sequestration and GHGs exchange.

Author Contributions

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

Funding

The research was partly financed by the Polish National Centre for Research and Development within the ERA-NET CO-FUND ERA-GAS Program (contract number ERA-GAS/I/GHG-MANAGE/01/2018) “GHG-Manage”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest

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

Appendix A

Table A1. Soil moisture content at the time of sampling tested soils [% v/v].
Table A1. Soil moisture content at the time of sampling tested soils [% v/v].
Forest TypeSummerAutumnWinterSpring
Coniferous Mature2.4%1.9%11.4%10.6%
Coniferous Young10.6%4.9%26.7%4.9%
Deciduous Mature4.3%6.7%22.0%20.5%
Deciduous Young14.8%25.6%24.3%32.2%
Mixed Mature11.0%5.6%32.8%11.5%
Mixed Middle-aged2.7%12.9%14.4%8.7%
Figure A1. Correlations between soil moisture content and microbial parameters based on mean values for each season for all studied forests (n = 24) (ac).
Figure A1. Correlations between soil moisture content and microbial parameters based on mean values for each season for all studied forests (n = 24) (ac).
Biology 10 00850 g0a1

References

  1. Babur, E.; Dindaroglu, T. Seasonal changes of soil organic carbon and microbial biomass carbon in different forest ecosystems. In Environmental Factors Affecting Human Health; Uher, I., Ed.; IntechOpen Limited: London, UK, 2020; pp. 1–21. [Google Scholar]
  2. FAO; UNEP. The State of the World’s Forests 2020. In Forests, Biodiversity and People; FAO: Rome, Italy, 2020. [Google Scholar] [CrossRef]
  3. Feng, H.; Guo, J.; Han, M.; Wang, W.; Peng, C.; Jin, J.; Song, X.; Yu, S. A review of the mechanisms and controlling factors of methane dynamics in forest ecosystems. For. Ecol. Manag. 2020, 455, 117702. [Google Scholar] [CrossRef]
  4. Forest Europe. 2020: State of Europe’s Forests; Forest Europe: Bratislava, Slovakia, 2020. [Google Scholar]
  5. Forest Europe. 2015: State of Europe’s Forests; Forest Europe: Madrid, Spain, 2015. [Google Scholar]
  6. EEA. Technical Report No 9/2006 European Forest Types Categories and Types for Sustainable Forest Management Reporting and Policy, 2nd ed.; European Environment Agency: Copenhagen, Denmark, 2007.
  7. Binkley, D.; Fisher, R.F. Ecology and Management of Forest Soils, 5th ed.; Wiley-Blackwell: Hoboken, NJ, USA, 2020. [Google Scholar]
  8. Lu, J.Z.; Scheu, S. Response of soil microbial communities to mixed beech-conifer forests varies with site conditions. Soil Biol. Biochem. 2021, 155, 108155. [Google Scholar] [CrossRef]
  9. Muscolo, A.; Panuccio, M.R.; Mallamaci, C.; Sidari, M. Biological indicators to assess short-term soil quality changes in forest ecosystems. Ecol. Indic. 2014, 45, 416–426. [Google Scholar] [CrossRef]
  10. Nielsen, M.N.; Winding, A. Microorganisms as indicators of soil health. In NERI Technical Report No. 388; National Environmental Research Institute: Silkeborg, Denmark, 2002; p. 84. [Google Scholar]
  11. Goberna, M.; Sánchez, J.; Pascual, J.A.; García, C. Surface and subsurface organic carbon, microbial biomass and activity in a forest soil sequence. Soil Biol. Biochem. 2006, 38, 2233–2243. [Google Scholar] [CrossRef]
  12. Morales, D.; Vargas, M.M.; de Oliveira, M.P.; Taffe, B.L.; Comin, J.; Soares, C.R.; Lovato, P. Response of soil microbiota to nine-year application of swine manure and urea. Cienc. Rural 2016, 46, 260–266. [Google Scholar] [CrossRef] [Green Version]
  13. Salazar, S.; Sánchez, L.E.; Alvarez, J.; Valverde, A.; Galindo, P.; Igual, J.M.; Peix, A.; Santa-Regina, I. Correlation among soil enzyme activities under different forest system management practices. Ecol. Eng. 2011, 37, 1123–1131. [Google Scholar] [CrossRef]
  14. Luo, P.; Han, X.; Wang, Y.; Han, M.; Shi, H.; Liu, N.; Bai, H. Influence of long-term fertilization on soil microbial biomass, dehydrogenase activity, and bacterial and fungal community structure in a brown soil of northeast China. Ann. Microbiol. 2015, 65, 533–542. [Google Scholar] [CrossRef] [Green Version]
  15. Santos, C.A.; Krawulski, C.C.; Bini, D.; Filho, T.G.; Knob, A.; Medina, C.C.; Filho, G.A.; Nogueira, M.A. Reclamation status of a degraded pasture based on soil health indicators. Sci. Agric. 2015, 72, 195–202. [Google Scholar] [CrossRef] [Green Version]
  16. Susyan, E.A.; Wirth, S.; Ananyeva, N.D.; Stolnikova, E.V. Forest succession on abandoned arable soils in European Russia—Impacts on microbial biomass, fungal-bacterial ratio, and basal CO2 respiration activity. Eur. J. Soil Biol. 2011, 47, 169–174. [Google Scholar] [CrossRef]
  17. Błońska, E.; Lasota, J.; Zwydak, M. The relationship between soil properties, enzyme activity and land use. For. Res. Pap. 2017, 78, 39–44. [Google Scholar] [CrossRef] [Green Version]
  18. Januszek, K.; Błońska, E.; Długa, J.; Socha, J. Dehydrogenase activity of forest soils depends on the assay used. Int. Agrophys. 2015, 29, 47–59. [Google Scholar] [CrossRef] [Green Version]
  19. Quilchano, C.; Marañon, T. Dehydrogenase activity in Mediterranean forest soils. Biol. Fertil. Soils 2002, 35, 102–107. [Google Scholar] [CrossRef]
  20. Fterich, A.; Mahdhi, M.; Mars, M. Seasonal changes of microbiological properties in steppe soils from degraded arid area in Tunisia. Arid Land Res. Manag. 2014, 28, 49–58. [Google Scholar] [CrossRef]
  21. Deng, S.P.; Parham, J.A.; Hattey, J.A.; Babu, D. Animal manure and anhydrous ammonia amendment alter microbial carbon use efficiency, microbial biomass, and activities of dehydrogenase and amidohydrolases in semiarid agroecosystems. Appl. Soil Ecol. 2006, 33, 258–268. [Google Scholar] [CrossRef]
  22. Błońska, E. 2010 Enzyme Activity in Forest Peat Soils. Folia For. Pol. Ser. A For. 2010, 52, 20–25. [Google Scholar]
  23. Błońska, E.; Lasota, J.; Gruba, P. Effect of temperate forest tree species on soil dehydrogenase and urease activities in relation to other properties of soil derived from loess and glaciofluvial sand. Ecol. Res. 2016, 31, 655–664. [Google Scholar] [CrossRef] [Green Version]
  24. Cheng, F.; Peng, X.; Zhao, P.; Yuan, J.; Zhong, C.; Cheng, Y.; Cui, C.; Zhang, S. Soil microbial biomass, basal respiration and enzyme activity of main forest types in the Qinling Mountains. PLoS ONE 2013, 8, e67353. [Google Scholar] [CrossRef]
  25. Condron, L.M.; Turner, B.L.; Cade-Menun, B.J. Chemistry and dynamics of soil organic phosphorus. In Phosphorus: Agriculture and the Environment; Sims, J.T., Sharpley, A.N., Eds.; ASA, CSSA, and SSSA: Madison, WI, USA, 2005; pp. 87–121. [Google Scholar]
  26. Hamarashid, N.H.; Othman, M.A.; Hussain, M.A.H. Effect of soil texture on chemical compositions, microbial populations and carbon mineralization in soil. Egypt. J. Exp. Biol. 2010, 6, 59–64. [Google Scholar]
  27. Hemkemeyer, M.; Dohrmann, A.B.; Christensen, B.T.; Tebbe, C.C. Bacterial preferences for specific soil particle size fractions revealed by community analyses. Front. Microbiol. 2018, 9, 149. [Google Scholar] [CrossRef] [Green Version]
  28. Moeskops, B.; Buchan, D.; Sleutel, S.; Herawaty, L.; Husen, E.; Saraswati, R.; Setyorini, D.; De Neve, S. Soil microbial communities and activities under intensive organic and conventional vegetable farming in West Java, Indonesia. Appl. Soil Ecol. 2010, 45, 112–120. [Google Scholar] [CrossRef]
  29. Wolińska, A.; Stępniewska, Z. Dehydrogenase Activity in the Soil Environment. In Dehydrogenases; Canuto, R.A., Ed.; IntechOpen Limited: London, UK, 2012; pp. 183–210. [Google Scholar]
  30. Zavišić, A.; Yang, N.; Marhan, S.; Kandeler, E.; Polle, A. Forest soil phosphorus resources and fertilization affect ectomycorrhizal community composition, beech P uptake efficiency, and photosynthesis. Front. Plant Sci. 2018, 9, 463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Niinistö, S.M.; Silvola, J.; Kellomäki, S. Soil CO2 efflux in a boreal pine forest under atmospheric CO2 enrichment and air warming. Glob. Chang. Biol. 2004, 10, 1363–1376. [Google Scholar] [CrossRef]
  32. Pregitzer, K.S.; King, J.S.; Burton, A.J.; Brown, S.E. Responses of tree fine roots to temperature. New Phytol. 2000, 147, 105–115. [Google Scholar] [CrossRef] [Green Version]
  33. Zhou, Y.; Li, M.H.; Cheng, X.B.; Wang, C.G.; Fan, A.N.; Shi, L.X.; Wang, X.X.; Han, S. Soil respiration in relation to photosynthesis of Quercus mongolica trees at elevated CO2. PLoS ONE 2010, 5, e15134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Hackl, E.; Buchmann, G.; Zechmeister-Boltenstern, S. Soil microbial biomass and rhizosphere effects in natural forest stands. Phyton 2000, 40, 83–90. [Google Scholar]
  35. Bolat, İ.; Kara, Ö.; Tunay, M. Effects of seasonal changes on microbial biomass and respiration of forest floor and topsoil under Bornmullerian firstand. Eurasian J. For. Sci. 2015, 3, 1–13. [Google Scholar] [CrossRef] [Green Version]
  36. Thoms, C.; Gleixner, G. Seasonal differences in tree species’ influence on soil microbial communities. Soil Biol. Biochem. 2013, 66, 239–248. [Google Scholar] [CrossRef]
  37. Babur, E.; Dindaroğlu, T.; Solaiman, Z.M.; Battaglia, M.L. Microbial respiration, microbial biomass and activity are highly sensitive to forest tree species and seasonal patterns in the Eastern Mediterranean Karst Ecosystems. Sci. Tot. Environ. 2021, 775, 145868. [Google Scholar] [CrossRef]
  38. Błońska, E.; Piaszczyk, W.; Staszel, K.; Lasota, J. Enzymatic activity of soils and soil organic matter stabilization as an effect of components released from the decomposition of litter. Appl. Soil Ecol. 2021, 157, 103723. [Google Scholar] [CrossRef]
  39. Sessitsch, A.; Weilharter, A.; Gerzabek, M.H.; Kirchmann, H.; Kandeler, E. Microbial population structures in soil particle size fractions of a long-term fertilizer field experiment. Appl. Environ. Microbiol. 2001, 67, 4215–4224. [Google Scholar] [CrossRef] [Green Version]
  40. Cools, N.; Vesterdal, L.; De Vos, B.; Vanguelova, E.; Hansen, K. Tree species is the major factor explaining C:N ratios in European forest soils. For. Ecol. Manag. 2014, 311, 3–16. [Google Scholar] [CrossRef]
  41. Lovett, G.M.; Weathers, K.C.; Arthur, M.A. Control of nitrogen loss from forested watersheds by soil carbon:nitrogen ratio and tree species composition. Ecosystems 2002, 5, 712–718. [Google Scholar] [CrossRef]
  42. Lovett, G.M.; Weathers, K.C.; Arthur, M.A.; Schultz, J.C. Nitrogen cycling in a northern hardwood forest: Do species matter? Biogeochemistry 2004, 67, 289–308. [Google Scholar] [CrossRef]
  43. Lucas-Borja, M.E.; Hedo, J.; Cerdá, A.; Candel-Pérez, D.; Viñegla, B. Unravelling the importance of forest age stand and forest structure driving microbiological soil properties, enzymatic activities and soil nutrients content in Mediterranean Spanish black pine (Pinus nigra Ar. ssp. salzmannii) forest. Sci. Total Environ. 2016, 562, 145–154. [Google Scholar] [CrossRef]
  44. Schultz, H.; Sigrid, H.; Stange, C.F. Species-specific differences in nitrogen uptake and utilization by six European tree species. J. Plant Nutr. Soil Sci. 2011, 174, 28–37. [Google Scholar] [CrossRef]
  45. Wang, F.; Li, Z.; Xia, H.; Zou, B.; Li, N.; Liu, J.; Zhu, W. Effects of nitrogen-fixing and non-nitrogen-fixing tree species on soil properties and nitrogen transformation during forest restoration in southern China. Soil Sci. Plant Nutr. 2010, 56, 297–306. [Google Scholar] [CrossRef] [Green Version]
  46. Wu, H.; Xiang, W.; Chen, L.; Ouyang, S.; Xiao, W.; Li, S.; Forrester, D.I.; Lei, P.; Zeng, Y.; Deng, X.; et al. Soil phosphorus bioavailability and recycling increased with stand age in Chinese fir plantations. Ecosystems 2020, 23, 973–988. [Google Scholar] [CrossRef]
  47. Zhou, Y.; Boutton, T.W.; Wu, X.B. Soil C:N:P stoichiometry responds to vegetation change from grassland to woodland. Biogeochemistry 2018, 140, 341–357. [Google Scholar] [CrossRef]
  48. Scheibe, A.; Steffens, C.; Seven, J.; Jacob, A.; Hertel, D.; Leuschner, C.; Gleixner, G. Effects of tree identity dominate over tree diversity on the soil microbial community structure. Soil Biol. Biochem. 2015, 8, 219–227. [Google Scholar] [CrossRef]
  49. Khlifa, R.; Paquette, A.; Messier, C.; Reich, P.B.; Munson, A.D. Do temperate tree species diversity and identity influence soil microbial community function and composition? Ecol. Evol. 2017, 7, 7965–7974. [Google Scholar] [CrossRef] [Green Version]
  50. Buragohain, P.; Nath, D.J.; Phonglosa, A. Role of microbes on carbon sequestration. Intern. J. Microb. Res. 2019, 1, 1464–1468. [Google Scholar]
  51. Bieganowski, A.; Ryżak, M.; Sochan, A.; Makó, A.; Barna, G.; Hernádi, H.; Beczek, M.; Polakowski, C. Laser diffractometry in the measurements of soil and sediment particle size distribution. Adv. Agron. 2018, 151, 215–279. [Google Scholar]
  52. Kitowski, I.; Sujak, A.; Wiacek, D.; Strobel, W.; Rymarz, M. Trace element residues in eggshells of Grey Heron (Ardea cinerea) from colonies of East Poland. North-West. J. Zool. 2014, 10, 346–354. [Google Scholar]
  53. Anderson, J.; Domsch, K. A physiological method for the quantitative measurement of microbial biomass in soils. Soil Biol. Biochem. 1978, 10, 215–221. [Google Scholar] [CrossRef]
  54. Šimek, M.; Kalčík, J. Carbon and nitrate utilization in soils: The effect of long-term fertilization on potential denitrification. Geoderma 1998, 83, 269–280. [Google Scholar] [CrossRef]
  55. Casida, L.E.; Klein, D.A.; Santoro, T. Soil Dehydrogenase Activity. Soil Sci. 1964, 98, 371–376. [Google Scholar] [CrossRef]
  56. Agnelli, A.; Ascher, J.; Corti, G.; Ceccherini, M.T.; Nannipieri, P.; Pietramellara, G. Distribution of microbial communities in a forest soil profile investigated by microbial biomass, soil respiration and DGGE of total and extracellular DNA. Soil Biol. Biochem. 2004, 36, 859–868. [Google Scholar] [CrossRef]
  57. Meisner, A.; Rousk, J.; Bååth, E. Prolonged drought changes the bacterial growth response to rewetting. Soil Biol. Biochem. 2015, 88, 314–322. [Google Scholar] [CrossRef] [Green Version]
  58. Koponen, H.T.; Bååth, E. Soil bacterial growth after a freezing/thawing event. Soil Biol. Biochem. 2016, 100, 229–232. [Google Scholar] [CrossRef]
  59. Schmiel, J.; Balser, T.C.; Wallenstein, M. Microbial stress-response physiology and its implications for ecosystem function. Ecology 2007, 88, 1386–1394. [Google Scholar] [CrossRef]
  60. Malchair, S.; Carnol, M. Microbial biomass and C and N transformations in forest floors under European beech, sessile oak, Norway spruce and Douglas-fir at four temperate forest sites. Soil Biol. Biochem. 2009, 41, 831–839. [Google Scholar] [CrossRef]
  61. Bérard, A.; Sassi, M.B.; Kaisermann, A.; Renault, P. Soil microbial community responses to heat wave components: Drought and high temperature. Clim. Res. 2015, 66, 243–264. [Google Scholar] [CrossRef]
  62. Gillespie, L.M.; Fromin, N.; Milcu, A.; Buatois, B.; Pontoizeau, C.; Hättenschwiler, S. Higher tree diversity increases soil microbial resistance to drought. Commun. Biol. 2020, 3, 377. [Google Scholar] [CrossRef]
  63. Prevost-Boure, N.C.; Maron, P.A.; Ranjard, L.; Nowak, V.; Dufrene, E.; Damesin, C.; Soudani, K.; Lata, J.C. Seasonal dynamics of the bacterial community in forest soils under different quantities of leaf litter. Appl. Soil Ecol. 2011, 47, 14–23. [Google Scholar] [CrossRef]
  64. Ananyeva, N.D.; Susyan, E.A.; Gavrilenko, E.G. Determination of the soil microbial biomass carbon using the method of substrate-induced respiration. Eurasian Soil Sci. 2011, 44, 1215–1221. [Google Scholar] [CrossRef]
  65. Thoms, C.; Gattinger, A.; Jacob, M.; Thomas, F.M.; Gleixner, G. Direct and indirect effects of tree diversity drive soil microbial diversity in temperate deciduous forest. Soil Biol. Biochem. 2010, 42, 1558–1565. [Google Scholar] [CrossRef]
  66. Henry, H.A.L. Soil freeze-thaw cycle experiments: Trends, methodological weaknesses and suggested improvements. Soil Biol. Biochem. 2007, 39, 977–986. [Google Scholar] [CrossRef]
  67. Koponen, H.T.; Jaakkola, T.; Keinänen-Toivola, M.M.; Kaipainen, S.; Tuomainen, J.; Servomaa, K.; Martikainen, P.J. Microbial communities, biomass, and activities in soil as affected by freeze thaw cycles. Soil Biol. Biochem. 2006, 38, 1861–1871. [Google Scholar] [CrossRef]
  68. Alvarez, E.; Torrado, V.M.; Marcos, M.L.F.; Diaz-Ravia, M. Microbial biomass and activity in a forest soil under different tree species. Electron. J. Environ. Agric. Food Chem. 2009, 8, 878–887. [Google Scholar]
  69. Díaz-Raviña, M.; Acea, M.J.; Carballas, T. Seasonal changes in microbial biomass and nutrient flush in forest soils. Biol. Fertil. Soils 1995, 19, 220–226. [Google Scholar] [CrossRef]
  70. Voříšková, J.; Brabcová, V.; Cajthaml, T.; Baldrian, P. Seasonal dynamics of fungal communities in a temperate oak forest soil. New Phytol. 2014, 201, 269–278. [Google Scholar] [CrossRef] [PubMed]
  71. Bauhus, J.; Paré, D.; Côté, L. Effects of tree species, stand age and soil type on soil microbial biomass and its activity in a southern boreal forest. Soil Biol. Biochem. 1998, 30, 1077–1089. [Google Scholar] [CrossRef]
  72. Graham, M.H.; Haynes, R.J. Organic matter status and the size, activity and metabolic diversity of the soil microflora as indicators of the success of rehabilitation of mined sand dunes. Biol. Fertil. Soils 2004, 39, 429–437. [Google Scholar] [CrossRef]
  73. Yan, J.; Zhang, D.; Zhou, G.; Liu, J. Soil respiration associated with forest succession in subtropical forests in Dinghushan Biosphere Reserve. Soil Biol. Biochem. 2009, 41, 991–999. [Google Scholar] [CrossRef]
  74. Chodak, M.; Niklińska, M. Effect of texture and tree species on microbial properties of mine soils. Appl. Soil Ecol. 2010, 46, 268–275. [Google Scholar] [CrossRef]
  75. Neto, M.S.; Scopel, E.; Corbeels, M.; Cardoso, A.N.; Douzet, J.M.; Feller, C.; Piccolo, M.C.; Cerri, C.C.; Bernoux, M. Soil carbon stocks under no-tillage mulch-based cropping systems in the Brazilian Cerrado: An on-farm synchronic assessment. Soil Tillage Res. 2010, 110, 187–195. [Google Scholar] [CrossRef]
  76. Lu, H.; Li, Z.; Fu, S.; Méndez, A.; Gascó, G.; Paz-Ferreiro, J. Combining phytoextraction and biochar addition improves soil biochemical properties in a soil contaminated with Cd. Chemosphere 2015, 119, 209–216. [Google Scholar] [CrossRef] [PubMed]
  77. Jia, G.M.; Cao, J.; Wang, C.; Wang, G. Microbial biomass and nutrients in soil at the different stages of secondary forest succession in Ziwulin, northwest China. For. Ecol. Manag. 2005, 217, 117–125. [Google Scholar] [CrossRef]
  78. Anderson, T.H.; Domsch, K.H. The metabolic quotient for CO2 (qCO2) as a specific activity parameter to assess the effects of environmental conditions, such as pH, on the microbial biomass of forest soils. Soil Biol. Biochem. 1993, 25, 393–395. [Google Scholar] [CrossRef]
  79. Novak, E.; Carvalho, L.A.; Santiago, E.F.; Portilho, I.I.R. Chemical and microbiological attributes under different soil cover. Cerne 2017, 23, 19–30. [Google Scholar] [CrossRef]
  80. Zang, C.; Hartl-Meier, C.; Dittmar, C.; Rothe, A.; Menzel, A. Patterns of drought tolerance in major European temperate firest trees: Climatic drivers and levels of variability. Glob. Chang. Biol. 2014, 20, 3767–3779. [Google Scholar] [CrossRef] [PubMed]
  81. Xu, X.; Schimel, J.P.; Janssens, I.A.; Song, X.; Song, C.; Yu, G.; Sinsabaugh, R.L.; Tang, D.; Zhang, X.; Thornton, P.E. Global pattern and controls of soil microbial metabolic quotient. Ecol. Monogr. 2017, 87, 429–441. [Google Scholar] [CrossRef] [Green Version]
  82. Böhme, L.; Langer, U.; Böhme, F. Microbial biomass, enzyme activities and microbial community structure in two European long-term field experiments. Agric. Ecosyst. Environ. 2005, 109, 141–152. [Google Scholar] [CrossRef]
  83. Bérard, A.; Sassi, M.B.; Renault, P.; Gros, R. Severe drought-induced community tolerance to heat wave. An experimental study on soil microbial processes. J. Soils Sediments 2012, 12, 513–518. [Google Scholar] [CrossRef]
  84. Anderson, T.H.; Domsch, K.H. Application of eco-physiological quotients (qCO2 and qD) on microbial biomasses from soils of different cropping histories. Soil Biol. Biochem. 1990, 22, 251–255. [Google Scholar] [CrossRef]
  85. Bauhus, J.; Khanna, P.K. The significance of microbial biomass in forest soils. In Going Underground—Ecological Studies in Forest Soils; Rastin, N., Bauhus, J., Eds.; Research Signpost: Trivandrum, India, 1999; pp. 77–110. [Google Scholar]
  86. Agnelli, A.; Ugolini, F.C.; Corti, P.; Pietramellara, G. Microbial biomass-C and basal respiration of fine earth and highly altered rock fragments of two forest soils. Soil Biol. Biochem. 2001, 33, 613–620. [Google Scholar] [CrossRef]
  87. Wardle, D.A.; Ghani, A. A critique of the microbial metabolic quotient (qCO2) as a bioindicator of disturbance and ecosystem development. Soil Biol. Biochem. 1995, 27, 1601–1610. [Google Scholar] [CrossRef]
  88. Priha, O.; Smolander, A. Microbial biomass and activity in soil and litter under Pinus sylvestris, Picea abies and Betula pendula at originally similar field afforestation sites. Biol. Fertil. Soils 1997, 24, 45–51. [Google Scholar] [CrossRef]
  89. Rahman, M.M.; Tsukamoto, J.; Rahman, M.M.; Yoneyama, A.; Mostafa, K.M. Lignin and its effects on litter decomposition in forest ecosystems. Chem. Ecol. 2013, 29, 540–553. [Google Scholar] [CrossRef]
  90. Smolander, A.; Loponen, J.; Suominen, K.; Kitunen, V. Organic matter characteristics and C and N transformations in the humus layer under two tree species, Betula pendula and Picea abies. Soil Biol. Biochem. 2005, 37, 1309–1318. [Google Scholar] [CrossRef]
  91. Wolińska, A.; Stępniewska, Z. Microorganisms abundance and dehydrogenase activity as a consequence of soil reoxidation process. In Soil Tillage & Microbial Activities; Miransari, M., Ed.; Research Singpost: Kerala, India, 2011; pp. 111–143. [Google Scholar]
  92. Józefowska, A.; Woś, B.; Pietrzykowski, M. Tree species and soil substrate effects on soil biota during early soil forming stages at afforested mine sites. App. Soil Ecol. 2016, 102, 70–79. [Google Scholar] [CrossRef]
  93. Hansson, K.; Olsson, B.A.; Olsson, M.; Johansson, U.; Kleja, D.B. Differences in soil properties in adjacent stand of Scots pine, Norway spruce and silver birch in SW Sweden. For. Ecol. Manag. 2011, 262, 522–530. [Google Scholar] [CrossRef]
  94. Brzezińska, M.; Stępniewski, W.; Stępniewska, Z.; Przywara, G. Effect of oxygen deficiency on soil dehydrogenase activity in a pot experiment with Triticale cv. Jago vegetation. Int. Agrophys. 2001, 15, 145–149. [Google Scholar]
  95. Ros, M.; Hernandez, M.; Garcia, C. Soil microbial activity after restoration of a semiarid soil by organic amendments. Soil Biol. Biochem. 2003, 35, 463–469. [Google Scholar] [CrossRef]
  96. Włodarczyk, T.; Stępniewski, W.; Brzezińska, M. Dehydrogenase activity, redox potential, and emissions of carbon dioxide and nitrous oxide from Cambisols under flooding conditions. Biol. Fertil. Soils 2002, 36, 200–206. [Google Scholar] [CrossRef]
  97. Almeida, D.; Filho, O.K.; Almeida, H.C.; Gebler, L.; Felipe, A.F. Soil microbial biomass under mulch types in an integrated apple orchard from Southern Brazil. Sci. Agric. 2011, 68, 217–222. [Google Scholar] [CrossRef]
  98. Stolnikova, E.V.; Ananyeva, N.D.; Chernova, O.V. The microbial biomass and its activity and structure in the soils of old forests in the European Russia. Eurasian Soil Sci. 2011, 44, 437–452. [Google Scholar] [CrossRef]
  99. Anderson, T.H.; Domsch, K.H. Ratios of microbial biomass carbon to total organic carbon in arable soils. Soil Biol. Biochem. 1989, 21, 471–479. [Google Scholar] [CrossRef]
  100. Sardans, J.; Rivas-Ubach, A.; Peñuelas, J. The C:N:P stoichiometry of organisms and ecosystems in a changing world: A review and perspectives. Perspect. Plant Ecol. Evol. Syst. 2012, 14, 33–47. [Google Scholar] [CrossRef]
Figure 1. Meteorological parameters in the study region in 2018 and 2019, based on information from the weather station of the Institute of Agrophysics PAS located in Lublin.
Figure 1. Meteorological parameters in the study region in 2018 and 2019, based on information from the weather station of the Institute of Agrophysics PAS located in Lublin.
Biology 10 00850 g001
Figure 2. Annual mean microbial parameters of soils: (a) BR—basal respiration; (b) Cmic—microbial biomass; (c) qCO2—metabolic quotient; (d) DHA—dehydrogenase activity, (e) and Cmic:Corg ratio) collected from different forests (average values ± standard deviation). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each parameter and for each forest).
Figure 2. Annual mean microbial parameters of soils: (a) BR—basal respiration; (b) Cmic—microbial biomass; (c) qCO2—metabolic quotient; (d) DHA—dehydrogenase activity, (e) and Cmic:Corg ratio) collected from different forests (average values ± standard deviation). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each parameter and for each forest).
Biology 10 00850 g002
Figure 3. Basal respiration (BR) in soil samples collected in summer, autumn, winter, and spring from different forest sites (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Figure 3. Basal respiration (BR) in soil samples collected in summer, autumn, winter, and spring from different forest sites (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Biology 10 00850 g003
Figure 4. Soil microbial biomass (Cmic) in samples collected in summer, autumn, winter, and spring from different forests (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Figure 4. Soil microbial biomass (Cmic) in samples collected in summer, autumn, winter, and spring from different forests (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Biology 10 00850 g004
Figure 5. The metabolic quotient (qCO2) in soil samples collected in summer, autumn, winter, and spring from different forests (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Figure 5. The metabolic quotient (qCO2) in soil samples collected in summer, autumn, winter, and spring from different forests (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Biology 10 00850 g005
Figure 6. Dehydrogenase activity (DHA) in soil samples collected in summer, autumn, winter, and spring from different forests (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Figure 6. Dehydrogenase activity (DHA) in soil samples collected in summer, autumn, winter, and spring from different forests (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Biology 10 00850 g006
Figure 7. Cmic:Corg ratio in soil samples collected in summer, autumn, winter, and spring from different forests (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Figure 7. Cmic:Corg ratio in soil samples collected in summer, autumn, winter, and spring from different forests (average values ± standard deviation; n = 3). (a) Coniferous Young (silt loam); (b) Coniferous Mature (loam sand); (c) Mixed Middle-aged (sandy loam); (d) Mixed Mature (silt loam); (e) Deciduous Young (silt loam); (f) Deciduous Mature (silt loam). Different letters indicate significant differences (a non-parametric Kruskal–Wallis test, p < 0.05, was carried out separately for each season and for each forest).
Biology 10 00850 g007
Figure 8. Principal component analysis (PCA)—biplot presentation between microbial and physicochemical parameters of six forest soils.
Figure 8. Principal component analysis (PCA)—biplot presentation between microbial and physicochemical parameters of six forest soils.
Biology 10 00850 g008
Figure 9. PCA plot of the properties of soils from mature and young coniferous, deciduous, and mixed forests.
Figure 9. PCA plot of the properties of soils from mature and young coniferous, deciduous, and mixed forests.
Biology 10 00850 g009
Figure 10. Contribution of Corg, Ntot, and P content in the analyzed forest soils.
Figure 10. Contribution of Corg, Ntot, and P content in the analyzed forest soils.
Biology 10 00850 g010
Table 1. Basic characteristic of the forests examined.
Table 1. Basic characteristic of the forests examined.
Forest TypeCoordinatesTree SpeciesContribution of Dominant Tree Species *Tree Age *Trees Biomass *
[%][year][m3 ha−1]
Coniferous Mature51.109675, 22.455327Scot pine Pinus sylvestris10010083
Coniferous Young51.036295, 22.993784Norway spruce Picea abies1003668
Deciduous Mature51.034831, 22.993789Oak Quercus robur
Aspen Populus tremula
70
30
71
63
193
87
Deciduous Young51.217270, 22.670177Sweet cherry Prunus avium1001226
Mixed Mature51.220832, 22.676460Scot pine Pinus sylvestris
Oak Quercus robur
80
20
85
85
306
65
Mixed Middle Aged51.111152, 22.453757Oak Quercus robur
Scot pine Pinus sylvestris
Black alder Alnus glutinosa
80
10
10
59
59
59
199
35
21
* Source: based on the Forest Data Bank, www.bdl.lasy.gov.pl (accessed on 30 August 2021).
Table 2. Properties of the different forest soils (0–15 cm depth; n = 3 avg. ± SD).
Table 2. Properties of the different forest soils (0–15 cm depth; n = 3 avg. ± SD).
Forest TypeConiferous MatureConiferous YoungDeciduous MatureDeciduous YoungMixed MatureMixed Middle Aged
Soil TypeDystric CambisolAlbic Dystric LuvisolLeptic CambisolDystric CambisolDystric CambisolAlbic Luvisol
Textureloamy sandsilt loamsilt loamsilt loamsilt loamsandy loam
Clay Content [%]1.6 ± 0.115.6 ± 0.115.4 ± 0.134.8 ± 0.105.1 ± 0.152.0 ± 0.08
Silt Content [%]22.9 ± 1.2972.4 ± 0.3569.7 ± 0.5364.3 ± 1.4069.4 ± 0.7627.7 ± 0.82
Sand Content [%]75.5 ± 1.4022.1 ± 0.3824.9 ± 0.6230.8 ± 1.4325.4 ± 0.6670.3 ± 0.89
Bulk Density [g cm−3]1.200.941.141.231.131.23
pHH2O4.37 ± 0.0064.82 ± 0.0154.80 ± 0.0065.34 ± 0.0064.22 ± 0.0214.79 ± 0.010
N [%]0.074 ± 0.0160.083 ± 0.0040.108 ± 0.0170.107 ± 0.0020.163 ± 0.0040.063 ± 0.001
Corg [%]1.20 ± 0.3630.91 ± 0.0681.73 ± 0.6701.14 ± 0.0652.25 ± 0.1030.82 ± 0.074
C/N16.12 ± 1.3911.04 ± 0.8515.71 ± 3.6810.66 ± 0.4313.81 ± 0.3013.17 ± 0.95
Ptot [mg/kg]277.0 ± 6.74273.6 ± 9.41316.7 ± 23.81456.8 ± 159.5458.2 ± 20.36264.2 ± 14.55
C/P42.4 ± 12.9432.3 ± 2.5151.2 ± 18.0318.6 ± 0.7847.2 ± 2.2929.7 ± 2.63
N/P2.60 ± 0.562.93 ± 0.163.21 ± 0.461.77 ± 0.0163.41 ± 0.162.25 ± 0.04
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Walkiewicz, A.; Bieganowski, A.; Rafalska, A.; Khalil, M.I.; Osborne, B. Contrasting Effects of Forest Type and Stand Age on Soil Microbial Activities: An Analysis of Local Scale Variability. Biology 2021, 10, 850. https://0-doi-org.brum.beds.ac.uk/10.3390/biology10090850

AMA Style

Walkiewicz A, Bieganowski A, Rafalska A, Khalil MI, Osborne B. Contrasting Effects of Forest Type and Stand Age on Soil Microbial Activities: An Analysis of Local Scale Variability. Biology. 2021; 10(9):850. https://0-doi-org.brum.beds.ac.uk/10.3390/biology10090850

Chicago/Turabian Style

Walkiewicz, Anna, Andrzej Bieganowski, Adrianna Rafalska, Mohammad I. Khalil, and Bruce Osborne. 2021. "Contrasting Effects of Forest Type and Stand Age on Soil Microbial Activities: An Analysis of Local Scale Variability" Biology 10, no. 9: 850. https://0-doi-org.brum.beds.ac.uk/10.3390/biology10090850

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