Ground Vegetation in Pinus sylvestris Forests at Different Successional Stages following Clear Cuttings: A Case Study

The impact of intensive forestry on various components of ecosystems has become the main subject of public and scientific debate in many regions in recent years. Forest ground vegetation is considered one of the most consistent and biodiversity-rich indicators of a certain stage of successional forest development. Therefore, changes in this forest component can potentially show the risks of forest damage due to clear-cutting and recovery trends. This study was carried out to identify the ground vegetation species diversity, including species composition and cover, also ground vegetation species relations with organic layer (forest floor) and upper mineral soil parameters at the different successional stages of the Pinus sylvestris L. stand development, including 1–2-year-old clear-cuts, and 6–130 years old stands. This study identified that the herb and dwarf shrub species were more light-demanding in the 2-year-old clear-cuts, as well as in the 6-year and 10-year old P. sylvestris stands compared to the middle-aged and mature forest stands. The dominant ground vegetation species, characteristic for the Pinetum vaccinio-myrtillosum forest type, were negatively dependent on the forest floor mass; they also had negative correlations with the concentrations of total P, K, Ca, and Mg in the forest floor and upper mineral soil but had positive correlations with the soil pH values and total N. The developed regression models of the percentage cover of mosses, herbs and dwarf shrubs according to the P. sylvestris stand age highlight the stabilization of the increase in the moss cover about 30 years after clear-cutting, with no clear trend for vascular species. The herbs and dwarf shrub species were highly variable during the stand rotation due to the species-specific characteristics and random factors rather than due to the influence of stand age. In this study, relatively short-term changes in ground vegetation species composition and percentage cover were determined after clear-cutting, but an important aspect is that new ground vegetation species appeared in the open areas, creating the potential for increasing species diversity. The clear-cutting system supports different species and numbers of herbs and mosses at different stages of stand development, which potentially increases the overall vegetation species diversity of the ecosystem.


Introduction
Forest management of high intensity, such as clear-cutting, directly impacts vegetation species diversity at the ecosystem level [1]. Although clear-cutting has recently received a lot of interest and opposition from public organizations, this type of silvicultural system remains an important and effective practice in Lithuania [2,3]. In recent years, the areas of mature forest stands have varied in a range of 439-466 thou. ha, and about 19 thou. ha of mature stands was harvested annually by clear-cutting in Lithuania [3]. Moreover, clearcutting accounted for 64% of the total area of felling in the State forests in 2021, and these fellings also prevailed in the Pinus sylvestris L. forests, especially in the private sector [3].
In the boreal and hemiboreal forest zone, the territories of managed mature forests are rapidly decreasing, and they are replaced by young forests after logging [4]. Overall, forest 6-year stand, the number of atypical species diminished, and the total species richness decreased to 25 species. The peak of P. sylvestris forest-related but light-demanding species continued until a 10-years old forest stand, and 31 species were observed at this age of the stand. The decreasing trend of species richness starts in the 15-year stand, and the number of species was in the following order: 28 species (15-year stand), 25 species (30-year stand), 17 species (70-year stand), and 12 species (mature forest) ( Table 1).

Ground Vegetation Cover at Different Forest Succession Stages
The assessment of percentage coverage of all ground vegetation species showed that the mosses were the dominant vegetation, covering more than 89% in the 30-and 70-year stands, and mature forests ( Table 2). The most frequent moss species in the mature stand were Hylocium splendens (Hedw.) Schimp.  Table 1). In the 30-and 70-year P. sylvestris stands, Pleurozium schreberi prevailed with the coverage of 50.7% and 67.8%, respectively, followed by Hylocium splendens (30.3% and 13.7%), and Ptilium cristacastrensis (5.6% and 10.2%). The clear-cutting tended to decrease the moss cover at the early successional stage from 56% in the 1-year-old clear-cut to 18% in the 2-year-old clear-cut, and 25% in the 6-year stand ( Table 2). In the 2-year-old clear-cut, new species such as Pohlia nutans (Hedw.) Lindb., Polytrichum commune L., and Polytrichum juniperinum Hedw. were detected. Furthermore, Pohlia nutans with a cover of 13.9% was fixed as dominated in the 6-year stand ( Table 1). The coverage of the mosses recovered up to 53.0% and 59.8% at the 10-year and 15-year stand, respectively (Table 2). Pleurozium schreberi (33.6%) prevailed in the 10-year stand, and Cirriphyllum piliferum (Hed.) Grout. (29.0%) in the 15-year stand ( Table 1).  Among the species of herbs and dwarf shrubs, Vaccinium myrtillus and Vaccinium vitisidaea were the dominants with the highest percentage cover of 10.6% and 5.8%, respectively, in the mature P. sylvestris stand, representing Pinetum vaccinio-myrtillosum forest type (Table 1). In comparison to the mature forest stand, the lower total coverage (14,5%) of dwarf shrubs was found in the 1-year-old clear-cut, and the next year, in the 2-year-old clear-cut, the coverage of herbal plants increased (23.9%) ( Table 2). Vaccinium myrtillus and Vaccinium vitis-idaea remained the most common species with the highest prominence value in the 1-year-old clear-cut (Table 1). The coverage of Vaccinium myrtillus and Vaccinium vitis-idaea was 6.6-6.8% and 1.9-2.2%, respectively, in the 1-2-year-old clear-cuts. The coverage of Pteridium aquilinum was 3.5-10.7 times higher in the 1-2-year-old clear-cuts than in the mature P. sylvestris forest. Some specific species, such as Luzula pilosa, Senecio vulgaris, and Trientalis europaea, were observed in the clear-cuts. The highest prominence value for Vaccinium myrtillus and Senecio vulgaris was fixed in the 2-year-old clear-cut. Significant differences in the species composition were observed in the 6-year stand, where Calluna vulgaris was the most frequent and dominant species with a coverage of 25.9% (Table 1). Other common species obtained were Vaccinium vitis-idaea, Rubus idaeus, Pteridium aquilinum, Festuca ovino, and Vaccinium myrtillus. The highest coverage of herbs and dwarf shrubs (59.0%) was found in the 10-year stand (Table 2). In this stand, the most frequent species were Calluna vulgaris (23.4%), Vaccinium vitis-idaea (13.4%), and Vaccinium myrtillus (7.8%) ( Table 1). The species Calluna vulgaris had the highest prominence value in the 6-year and 10-year stands. In the 10-year stand, the coverage of both Vaccinium myrtillus and Vaccinium vitis-idaea was higher in the 6-year stand. Vaccinium vitis-idaea was identified as the most frequent in the 15-year and 30-year stands, covering 24.8% and 6.3%, respectively, and it was found in all research sites. Calluna vulgaris, covering 2.3%, was identified as the less frequent species, found in 25% of the plots of the 15-year stand (Table 1). In the 15-year stand, the frequent species were Festuca ovino, Melampyrum pratense L., and Vaccinium myrtillus. In the 30-year stand, the coverage of Vaccinium vitis-idaea was followed by Vaccinium myrtillus, Festuca psammophila (Hack. exČelak.) Fritsch, and Luzula Pilosa. In the 70-year stand, both Vaccinium myrtillus and Vaccinium vitis-idaea were found in all research plots, but Vaccinium myrtillus coverage was higher than Vaccinium vitis-idaea.
The regression models of the percentage cover of mosses, in addition to herbs and dwarf shrubs, according to the P. sylvestris stand age, were developed ( Figure 2). Regardless of the successional status of individual species, the increase in the moss cover stabilized approximately 30 years after clear-cutting. Despite the relatively high variability of the data, the faster growth of herbs and dwarf shrubs was fixed at the early stage of the stand; this trend changed insignificantly in the 30-50-year-old stand and showed a slight decrease trend in the mature forest. The plant community obtained in 2-year-old clear-cut has a higher Shannon species diversity index value and was assigned as the most diverse in comparison with the 1-year old clear-cut and the P. sylvestris stand of 6-130 years old ( Table 3). The lowest diversity of the plants was identified in the 70-year stand. The obtained results were comparable with the total species richness values obtained at all studied sites (see Table 1). However, the highest equitability was found for the 10-year stand followed by the 2-year-old clear-cut and 6-year stand.
For herbs and dwarf shrubs, the highest Ellenberg indicator values for light were obtained in the 2-year-old clear-cut, also in the 6-year and 10-year P. sylvestris stands, followed by a consistent downward trend towards the oldest studied stands ( Table 4). The Ellenberg indicator values for temperature inconsistently increased, while the values for the soil moisture showed a decreasing trend with increasing stand age. The Ellenberg values for nitrogen were higher in the 1-2-year-old clear-cuts compared to the older stands.
For mosses, there were no clear trends of Ellenberg indicator values for light, temperature, moisture, and pH over the entire stand rotation (Table 4). However, the highest Ellenberg indicator values for light and soil moisture were fixed in the 10-and 15-year stands. The values for pH were higher in the stands older than 10 years old compared to the stands of an early stand development stage. Table 3. Shannon species diversity index (H ), and equitability (E) of ground vegetation species at different Pinus sylvestris sites, including 1-2-year-old clear-cuts and 6-130 years old stands. Different letters a, b, and c indicate significant differences among sites (p < 0.05).

Clear-Cut (Years)
Age of Pinus sylvestris Stand (Years)

The Ground Vegetation Relations with the Forest Floor and Mineral Soil Parameters
The species richness of the mosses had relatively strong relation with the forest floor mass (R 2 = 0.733), i.e., the lower number of the mosses was identified under the higher forest floor mass ( Figure 3). The species richness of herbs and dwarf shrubs also showed a downward trend with the increase in forest floor mass. The higher forest floor mass tended to decrease the cover of herbs and dwarf shrubs, while no relation with the moss cover was obtained. The cover of Vaccinium myrtillus moderately (0.50 ≤ r ≤ 0.75) positively correlated with forest floor pH and total N ( Table 5). The moderately negative correlations were found between the cover of Vaccinium myrtillus and organic C, total P, K, Ca, and Mg in forest floor, and total P, K, Ca, and Mg, in addition to P 2 O 5 , K 2 O, and mobile Ca in the upper (0-10 cm) mineral soil layer. The cover of Vaccinium vitis-idaea moderately positively correlated with the pH value obtained in the forest floor and upper mineral soil layer, as well as with the organic C and total N of the mineral soil (Table 5). Moderately negative correlations were found between the Vaccinium vitis-idaea and total P, K, Ca, and Mg obtained in both forest floor and mineral soil, as well as with P 2 O 5 , K 2 O, and mobile Ca of mineral soil.
Other ground vegetation species were less related to the chemical parameters of the forest floor and upper mineral soil (Table 5). However, for example, Luzula pilosa showed a moderately (0.50 ≤ r ≤ 0.75) positive relationship with forest floor pH and a weak (0.30 ≤ r ≤ 0.50) relationship with organic C and total N of the forest floor. Moderately negative correlations were found between the cover of Luzula pilosa and total P, K, and Ca obtained in both forest floor and upper mineral soil layer, as well as with total Mg, P 2 O 5 , and mobile Ca obtained in mineral soil. The relationship between the cover of Hylocium splendens with total P, K, and Ca in forest floor and mineral soil, and also with total Mg, K 2 O, and mobile Ca and Mg in mineral soil mostly showed a correlation coefficient about r = 0.5. Similar relationships were identified between the cover of the Ptilium crista-castrensis and Dicranum polysetum with the chemical parameters of the forest floor and upper mineral soil layer.

The Ground Vegetation Biomass
The highest total biomass of ground vegetation was found in the 10-and 15-year P. sylvestris stand (Table 6). In the 30-, 70-and 130-year stands, the total biomass of ground vegetation was 1.6-3.6 times lower than in the 10-and 15-year stands. The total biomass of total ground vegetation in the 6-year stand was different from the remaining values. The highest mass of mosses was obtained in the 1-year-old clear-cut, which was comparable with the values obtained in the mature P. sylvestris stand. The lowest mass of mosses was found in the 2-year-old clear-cut, and then increased consistently up to the 30 years age. The biomass of herbs and dwarf shrubs showed large variability, but the highest values were obtained in the 10-and 15-year stand. The obtained results indicated that the biomass of herbs and dwarf shrubs tended to be more species-specific than depended on stand age.
The herbs and dwarf shrubs in different sites consisted of similar species groups with the highest biomass proportion of the dominant species Vaccinium myrtillus and Vaccinium vitis-idaea with slight exceptions at the 6-year stand ( Figure 4A). In the period from 6to 15-year stand, a larger proportion of Calluna vulgaris biomass was observed. Mosses consisted mainly of the same species groups with the highest biomass proportion of the Pleurozium schreberi, Hylocium splendens, and Dicranum sp. (Figure 4B). Table 5. Relations between the cover of ground vegetation species and forest floor (FF) and mineral soil (MS) chemical parameters (Pearson's r and p-value).

Discussion
Ground vegetation species richness and diversity varied between different study sites, including 1-2-year-old clear-cuts, and the P. sylvestris stands of different ages recovered after clear-cutting. The results of this study indicated that the most abundant moss species were Hylocium splendens, Pleurozium schreberi, and Ptilium crista-castrensis (see Table 1), and well represented the moss species typical for the Pinetum vaccinio-myrtilliosa forest type [11,28]. Previous studies also stated that Ptilium crista-castrensis and Hylocomium splendens were dominated species in the comparable P. sylvestris forest stands [29]. The dominant species of vascular plants were Pteridium aquilinum at an early stage after clear-cutting and both, Vaccinium myrtillus and Vaccinium vitis-idaea, which were obtained at all studied sites of different stand ages (see Table 1). Total species richness tended to increase at the clear-cuts, and this lasted up to 10 years after clear-cutting. However, the species of vascular plants showed a decreasing trend, and gradually the prevailing moss species was recorded. The predominance of light-tolerant vascular plants in the boreal forest zone was identified by several studies [30,31].
Our findings were consistent with previous studies that showed differences in ground vegetation cover and species diversity along a stand age gradient [10,11,18]. The species composition and coverage vary since the vegetation colonizes the available resources after clear-cutting-induced disturbance and consequently reaches the pre-cut forest level. Vegetation succession begins with the removal of overstory vegetation, which initiates the growth of stress-tolerant vegetation species, and vegetation species composition tend to become associated with the forest at later stages of stand development [28,[32][33][34].
In our study, the regression models were developed based on the moss and vascular plant cover at different P. sylvestris stands aged up to 130 years, including clear-cutting (see Figure 2). The obtained trends were comparable with those developed for the Pinus contorta Douglas. ex Loudon stands in Ref. [19]. The obtained trends indicated compositional stability of the species in the middle-aged stand, mainly between 40 and 60 years of age. Similarly, along with the decreased cover of vascular species, increased cover of moss species was reported. According to Ref. [11], the moss cover in the P. sylvestris stands during the stand rotation is not uniform. Although a decrease in moss cover in the first years after felling could be distinguished, the moss cover already recovers in young stands. Furthermore, the moss cover in the middle-aged stands tends to be higher than in the mature stands.
The current study found that the mosses were the dominant vegetation, covering more than 89% with the highest cover of Hylocium splendens, Pleurozium schreberi and Ptilium crista-castrensis in the P. sylvestris stands older than 30 years old (see Tables 1 and 2). However, in the clear-cuts and young forest stands, representing an early successional stage of forest development, the mosses covered approximately 18-56%. The clear-cuts were characterized by the occurrence of new species, such as Pohlia nutans, Polytrichum commune, and Polytrichum juniperinum. Overall, a period of 10-15 years after clear-cutting can be considered as the period of successful recovery of the moss cover in the Pinetum vaccinio-myrtilliosa forest type. In accordance with the present results, previous studies have demonstrated that the cover of vascular plants was 14-40%, and the moss cover was 80-85% in the mature P. sylvestris stands growing at the Vaccinio-myrtillosum forest type [27,29]. At the early stage of stand development after clear-cutting, the increase or even doubling in vascular plant diversity is detected [34,35]. At later stages, the diversity of mosses and lichen species increases, but this process is longer and requires 5-6 years to find the significant increase [36]. Furthermore, interrelationships between different vegetation groups, such as mosses and vascular plants, are often studied as an indicator, having a great influence on the dominant plants in the ecosystem [37,38]. The mentioned study revealed that all studied mosses suppressed the regeneration of vascular plants, but differences between the moss species were identified [38]. Therefore, the reduction in vascular species richness often leads to a higher richness and diversity of non-vascular species [39].
Among all abiotic factors, intensive forest management, including clear-cutting, directly impacts light and soil nutrient availability, and these are the factors of great significance for recovering plant communities. The results of our study indicated that light was an important factor for herbs and dwarf shrubs. Therefore, the highest Ellenberg indicator values for light were obtained at the early stages of stand development, and these values decreased with the increase in stand age (see Table 4). However, the highest Ellenberg indicator values for light for the mosses were found in the relatively young to middle-aged stands compared to the clear-cuts and mature forests. Our data revealed that herbs and dwarf shrubs were more related to changes in soil temperature than to changes in moisture compared to mosses. These findings may be explained by the fact that the mosses are characterized as tolerant to shade stress and less dependent on soil moisture and fertility than vascular plants [40,41]. Overall, understory diversity increases with the increase in resource availability [42][43][44]. For example, the higher abundance of Deschampsia flexuosa on clear-cuts and the sites at early stages of forest development after clear-cutting was detected due to the increased nitrogen availability [45][46][47][48][49].
It is interesting to note that higher forest floor mass responded to the lower number of the moss species and tended to decrease the number of herbs and dwarf shrubs species (see Figure 3). However, the higher forest floor mass caused the lower percentage cover of herbs and dwarf shrubs but did not change the moss cover. It was already known from previous studies that moss colonization is associated with structural elements of the forest floor [50]. Overall, as noted by Ref. [13], most species typical for older forests could successfully survive in the clear-cuts. By increasing spatial diversity, clear-cutting creates the opportunity for populations of vegetation species other than bryophytes, which are more typical for the old-growth forest floor [10]. Our study also indicated large variability among the groups of the dominant herbs and dwarf shrubs, also moss species at different P. sylvestris sites (see Figure 4). For example, clearcutting greatly reduced Vaccinium myrtillus cover, which recovered slightly after 10 years and continued to increase in older stands. These results seem to be consistent with other studies which found a reduction in Vaccinium myrtillus abundance soon after clear-cutting and full recovery to mature stand level 60-80 years after clear-cutting [51]. Despite the ability to survive after forest harvesting, Vaccinium myrtillus was negatively affected by clear-cutting and this was well documented [47,49,52]. We also found the decrease in the cover of Pleurozium schreberi and Hylocomium splendens as a response to clear-cuttings, and these moss species are often defined as sensitive to final forest harvest [53].
When comparing the relations between the abundance of selected ground vegetation species and the chemical parameters of forest floor and upper mineral soil, the dominant species such as Vaccinium myrtillus, Vaccinium vitis-idaea, and Hylocium splendens showed moderately strong positive correlations with the soil pH values and total N concentrations (see Table 5). However, the soil chemical parameters such as total P, K, Ca, and Mg often responded in negative correlation to the mentioned species. No significant relations with soil nutrients were found for the Pleurosium shreberi, and low correlation was found for Dicranum polysetum and Melampyrum pratense. Accordingly, we identified species relation to chemical parameters as species-specific response.
To conclude, the potential changes and variations in ground vegetation species indices in the P. sylvestris sites of different successional stages following clear-cuttings were investigated in this study. We have focused on the main maturity groups of P. sylvestris forest stands, typical for the hemiboreal forests: the clear-cuts of the first two years after felling, young stands of the 6, 10, and 15 years, middle-aged stands of 30 years, premature stands of 70 years, and mature stands 130 years old. This study has raised an important question about the relations between the ground vegetation species and chemical parameters of forest floor and upper mineral soil. The study findings indicated that the ground vegetation of the clear-cuts is more light-demanding, while the most shaded vegetation was obtained in the middle-aged stands. The ground vegetation species, obtained in the P. sylvestris stands assigned to the Pinetum vaccinio-myrtillosum forest type, showed the dependence on the parameters of the forest floor and upper mineral soil. These relations with the specific soil factors, along with the light factor, provide additional information in the evaluation of the ground vegetation diversity during the stand rotation following clear-cutting. It should be recognized that our data and insights do not give a full-scale view of ground vegetation development over the P. sylvestris stand rotation, but the results highlight the need to explore the links among a variety of forest ecosystem components. Additionally, the obtained specific findings could help to harmonize the aspects of intensive and close-to-nature forest management, taking into account the diversity of ground vegetation, and could help in decision making. Previous studies have shown that, for example, the close-to-nature system is not fully compatible with sustainable forest management in terms of biodiversity [5].
The trends identified in this study summarize several insights: firstly, changes in ground vegetation species composition and percentage cover following clear-cutting were identified as a relatively short-term process; secondly, forest ecosystems affected by clearcutting, as an external disturbance, are forced to start a succession process, but vegetation changes tended to be related to the forest floor and upper mineral soil properties, the composition of the dominant forest-related species, and dependencies among the species of different groups; thirdly, new ground vegetation species occurring in the open-areas result in increased species richness, which has a positive effect on the overall ecosystem biodiversity. Therefore, clear-cutting system supports different species and a number of herbs and mosses at different stages of stand development. Forest ecosystems may lose certain vegetation species because some of them are limited by the frequency of disturbance rather than by the severity of disturbance [5].

Study Sites
Almost 34% of Lithuania is covered by forests [3]. Formally, Lithuanian forests belonged to the transitional hemiboreal forest zone of Europe with a prevalence of mixed deciduous and coniferous stands [54]. This study analyzes Scots pine (Pinus sylvestris L.) stands, which are the predominant coniferous forests (34.6%) in Lithuania [3].
The study on the understory vegetation in pure Pinus sylvestris forests at different successional stages following clear-cuttings was performed at the Trakai Regional Division of State Forest Enterprise in south-eastern Lithuania (54 • 38'16.19" N; 24 • 56'3.59" E) in June-September 2020 and 2021. The 1-year-old and 2-year-old clear-cuts, and 6-, 10-, 15-, 30-, 70-, and 130-year-old managed P. sylvestris stands were selected for this study. The 1-year-old clear-cut sites were established after the clear-cutting of the mature P. sylvestris forest, cultivated as even-aged stands, in 2019. The stems and logging residues, except stumps and roots, were removed from the cutting site. The 1-year-old clear-cut site was reestimated the next year, obtaining the 2-year-old clear-cut site. More specifically, 2-year-old clear-cut site was reforested by 2-year-old P. sylvestris seedlings. Other study sites in the stands, representing different stand ages, were selected in the adjacent forest sites with a known forest management history. A pure P. sylvestris stand, which has recovered after clear-cutting managed mature stands, was chosen for the study. All selected sites were similar in climatic and edaphic conditions. The average annual temperature in Lithuania is 7.4 • C (1991-2020 climatic normal). The hottest month in Lithuania is July; the coldest is January. The average annual precipitation is 695 mm (1991-2020 climatic normal) [55]. All selected study sites were characterized by the Pinetum vaccinio-myrtilliosa forest type [28]; the oligotrophic mineral soil of a normal moisture regime (Nbl) [56]; and nutrient-poor sandy Arenosols with the light soil texture [57]. Typically, the underbrush species in the Pinetum vaccinio-myrtilliosa forest type are Sorbus aucuparia L., Frangula alnus Mill., and Juniperus communis L. [28]. The mean characteristics are given in Table 7.

Ground Vegetation Sampling and Analyses
The vegetation was assessed in July-August 2020 and 2021 when the vegetation cover was the highest and most diverse. For the assessment of ground vegetation, forest floor, and mineral soil parameters, a 2000 m 2 circular plot was established at each site. In 2020, the permanent plot was established one year after clear-cutting (1-year-old clear-cut), and the investigations were repeated in the same plot 2 years after clear-cutting (2-year-old clear-cut). The temporary plots were established in the 6-, 10-, 15-, 30-, 70-, and 130-year-old P. sylvestris stands. The vegetation assessment was carried out on five sample plots of 10 × 10 m at each site, located in the directions of 0°, 90°, 180°, and 270°azimuths from the centre and one selected in the centre of the plot. At each site, the vegetation assessment was performed on 25 vegetation quadrates (in 1 × 1 m subplots). The frame of 1 dm 2 square mesh was used to assess the cover (%) of each species in the vegetation quadrate visually. Then, the mean value of each species cover was calculated per site. The frequency of each species was determined as follows: the number of subplots in which the species was detected divided by a total number of subplots and multiplied by 100 [59]. The detailed scheme of the sample site, identifying the places for the ground vegetation assessment and soil sampling, was given in Ref. [60]. The ground vegetation species were recorded separately for moss layer and herb together with dwarf shrub layers < 0.5 m and > 0.5 m [61]. The tree layer was recorded for characteristics of the selected forest sites. The total richness of ground vegetation species was estimated as the total number of species obtained at each site, i.e., in the 1-year-old and 2-year-old clear-cuts, and 6-, 10-, 15-, 30-, 70-, and 130-year-old P. sylvestris stands. The mean richness of ground vegetation species was estimated as the mean of the species number obtained in every five sample plots within each study site (n = 5).
Aboveground biomass of ground vegetation was harvested within 25 × 25 cm squares from four subplots, located systematically in each of five sample plots within each site (1-year old and 2-year-old clear-cuts, and 6-, 10-, 15-, 30-, 70-, and 130-year-old P. sylvestris stands). All the living vegetation within a square was cut using scissors close to the ground surface. The harvested vegetation was put into plastic bags and transported to the laboratory. In the laboratory, litter and logging residues were removed from the samples, leaving only the living parts of the ground vegetation species. The moss species were separated into five groups: Pleurozium schreberi, Hylocium splendens, Ptilium crista-castrensis, Dicranum sp., and other mosses. The herbs and shrubs were separated into four groups: Vaccinium myrtillus, Vaccinium vitis-idaea, Calluna vulgaris, and other plants. The samples were oven-dried at 105 • C to a constant mass and weighed.

Soil Sampling and Analyses
The organic layer (forest floor) was sampled within a 25 × 25 cm square using a metallic frame. The composite samples of the forest floor were obtained from four subsamples (n = 4), oven-dried at 105 • C to a constant mass, and weighted. Mineral soil was sampled from 0-10 cm by a metallic soil auger in August-September 2020. Four composite samples were combined from 10 subsamples collected systematically in each sample plot. The forest floor and mineral soil samples were analysed for pH, which was determined in a 0.01 M CaCl 2 suspension (ISO 10390:2005); organic carbon (Corg) was determined using a dry combustion method with a total carbon analyzer Analytic Jena multi EA 4000, Jena, Germany, and total nitrogen (N) by the Kjeldahl method (ISO 11261). The concentration of mineral N was determined by the spectrometric method (ISO 14256-2) in 1 M KCl extraction. Soil available K 2 O and P 2 O 5 were determined in ammonium lactate-acetic acid extractant [62] by the Shimadzu UV 1800 spectrophotometer (Shimadzu, Kyoto, Japan) and JENWAY PFP7 flame photometer, respectively. Soil available Ca and Mg were determined by atomic absorption spectrometer AAnalyst 200 (PerkinElmer, Waltham, MA, USA) (laboratory method LVP D-13:2016).

Calculations and Statistical Analysis
To assess the species diversity at different sites, the Shannon diversity index (H ) and equitability (E) were calculated [63]. The Shannon diversity index (H ) was calculated by the Equation (1): here, n i is the coverage of the ith species; N is the total coverage of vegetation in the site. Equitability (E) was calculated by the Equation (2): here H'max = ln(N); N is the total coverage of vegetation in the site. Ellenberg's original values, quantified originally for Central Europe, were used to identify the differences in the ecological performance of species under specific conditions in the clear-cuts and in the stands of different ages [45]. The average cover of vascular plants and mosses was used to calculate the mean weighted ecological scale values of Ellenberg by the Equation (3): here Ex-weighted average of environmental factor, x i -indicator value of ith species ecological factor, n i -the average coverage of ith species, and N-total coverage of species that are not indifferent to that ecological factor. To calculate and draw a Venn diagram, we used an online program (http://bioinformatics. psb.ugent.be/webtools/Venn/, accessed on 30 August 2022). For the normality of the variables, the Shapiro-Wilk test was used, and the normality hypothesis was rejected with a 0.05 significance level. Therefore, the Kruskal-Wallis analysis of variance (ANOVA) test was used to ascertain the significant differences in variables between the sites. Pearson correlation was used to measure the strength of a linear correlation between two variables. Data are presented as means ± standard error (SE). Different letters next to the mean values are statistically significantly different at p < 0.05 between the sites. Statistical analyses were conducted using STATISTICA 12.0 (StatSoft. Inc., TULSA, OK, USA, 2007) software.