Effect of tree species identity and related habitat parameters on understorey bryophytes – interrelationships between bryophyte, soil and tree factors in a 50-year-old experimental forest

Overstorey tree species influence both soil properties and microclimate conditions in the forest floor, which in turn can induce changes in ground bryophyte communities. The aim of the study was to investigate the effect of tree species identity and the most important habitat factors influencing understorey bryophytes. We assessed the effect of 14 tree species and related habitat parameters, including soil parameters, vascular plant presence and light intensity on bryophytes in monospecific plots covered by nearly fifty-year-old trees in the Siemianice Experimental Forest (Poland). The canopy tree species determined bryophyte species richness and cover. The strongest differences were observed between plots with deciduous and coniferous trees. Soils with a more acidic pH and lower content of macronutrients supported larger bryophyte coverage. We also found a positive correlations between vascular plants and availability of light as well as bryophyte species richness. Tree species identity and differences in habitat conditions in the forest floor lead to changes of ground bryophyte richness, cover and species composition. Consequently, the changes in the dominant tree species in the stand may result in significant repercussions on ground bryophyte communities. We indicated that the introduction of alien tree species, i.e. Quercus rubra, has an adverse effect on bryophyte communities and suggested that the selection of tree species that contribute to the community consistent with the potential natural vegetation is highly beneficial for maintaining ground bryophyte biodiversity.

It would be hard not to mention the main plant components of the forest, which are trees. Trees constitute the majority of biomass in forest complexes and are responsible for the functioning of the forest interior (Ellison et al. 2005). They are particularly important for soil-forming processes since leaf litter and root decomposition are the most important source of soil organic matter (Augusto et al. 2003). Tree species composition has been proven to shape the chemical properties of litter and thereby could strongly influence bryophyte communities indirectly via alternation of microclimatic conditions in a forest stand (e.g. Arno et al. 2012).
Several studies have been conducted to elucidate bryophyte distribution and richness in forests to give insight into the different factors that influence their growth (e.g. Gosselin et al. 2017;Jansová and Soldan 2006;Király and Ódor 2010;Vellak and Ingerpuu 2005). Most of the studies concern epiphytic bryophytes, and thus the effect of tree species on bryophyte diversity and composition is relatively well recognised. The main factors are related to tree age, inclination of the trunk or branches, bark structure, bark pH and water capacity (e.g. Fritz et al. 2009;Király and Ódor 2010). On the other hand, little is known about the tree effect on ground bryophytes. Moreover, most of the previous studies were conducted in natural old-growth or mature managed forests that comprised a mixed composition of tree species, which certainly makes it difficult to estimate the direct effect of tree species on bryophytes. As regards ground bryophytes, Jagodziński et al. (2018) studied the effect of tree species on bryophytes in the disturbed ecosystem of a reclaimed lignite mine spoil heap and concluded that the main factors affecting bryophyte diversity are light availability, soil pH and C/N ratio. Another factor was reported by Weibull and Rydin (2005) who found that the amount of litter influenced species richness of bryophytes growing on granite boulders in mixed forests. Márialigeti et al. (2009) recognised substrate availability, litter cover, stand age, and stand structural diversity as the factors influencing the development of a bryophyte layer in temperate-mixed forests; however, soil properties were not investigated in this study. Stand structure attributes were also recognised as a crucial determinant of ground bryophyte communities in a naturally regenerated and planted forest in China (Yan et al. 2013). Finally, soil salinity, volumetric water content and vascular plant species cover proved to significantly affect the composition of bryophytes in Scots pine forests (Tilk et al. 2018). As a general rule, coniferous and deciduous trees provide different habitat conditions with respect to light intensity, moisture availability, litter deposition, soil pH and fertility (e.g. Augusto et al. 2002;Márialigeti et al. 2009). For these reasons, we suppose that coniferous and deciduous trees will harbour different bryophyte assemblages.
In the present study, we analysed the effect of different tree species growing in monospecific stands in an experimental forest on understorey bryophyte communities. Siemianice Experimental Forest (western Poland) is a research area unique on a global scale that comprises monoculture plots with 14 different tree species planted in 1970 and 1971 on the same type of soil substrate. The evaluation of the effect of tree species and other habitat factors on bryophytes in such an experimental forest seems to be highly beneficial compared to natural forest communities since the experimental setup allows for a convenient reduction of confounding factors such as differences in tree stand age, microclimate, topographic factors, admixture of other tree species and other external factors. To the best of our knowledge, there have been no studies concerning the direct effect of tree species on understorey bryophytes under such conditions. We assume that the 50-year period from planting the trees is sufficient to establish differentiated microhabitat conditions and create different ecological niches for bryophytes. Consequently, we hypothesise that tree species, albeit indirectly, has a major influence on bryophyte species composition and diversity in the forest floor.
The overall aims of the study were (1) to determine if bryophyte species richness and cover varies between plots with different tree species, between plots with deciduous and coniferous trees as well as between plots with native and alien to the flora of the region tree species; (2) to recognise the most important habitat factors influencing species richness and composition of understorey bryophyte communities; (3) to determine to what extent bryophyte species composition captured in ecological indicator values corresponds to the measured habitat variables; (4) to recognise tree species which are of particular conservation value in terms of preservation of the biodiversity of bryophytes in the forest interior.

Study area and sampling
The study was conducted in the Siemianice Experimental Forest, W Poland (51º14.87ʹ N, 18º06.35ʹ E, 150 m a. s. l.). According to the updated Köppen-Geiger climate classification (Kottek et al. 2006), the study area is classified under a temperate-oceanic climate (Cfb) with a mean temperature of 8.2 °C and mean annual precipitation of 591 mm.
Sampling plots of size 1 × 1 m were located in the central part of each forest plot to avoid edge effects. Bryophyte and vascular plant species richness and cover were evaluated (see section hereafter). Then, 3 soil subsamples were collected to a depth of ca. 20 cm using a shovel; organic matter (O horizon) was removed before sampling. Collected subsamples were placed in foil bags and homogenised to create a single composite sample representative for each plot. Soils were transferred into new foil bags and stored at 4 ºC. Each sample was collected and processed using a sterile kit to avoid contamination. After sampling at a given plot, the shovel was cleaned, sterilized with 70% ethanol, and treated with fire. Altogether, 53 soil samples were collected on 3 and 4 June 2018. Soils were used for analyses of chemical properties (see section hereafter; see also Stefanowicz et al. 2021).
Bryophyte and vascular plant species richness, cover, and light intensity Bryophyte and vascular plant species composition in each sampling plot were determined as follows. Each bryophyte sample was taken apart and studied in detail under a microscope to avoid overlooking any mixed species or hidden fragments of mosses and liverworts. The bryophyte material was deposited in the OSTR herbarium (University of Ostrava), the nomenclature of bryophytes follows Hill et al. (2006), while vascular plant species were identified according to Rutkowski (2014). Cover of bryophytes and vascular plants was estimated on a percentage scale within each sampling plot using digital photos of vegetation. The borders of each sampling plot were clearly marked with a coloured cord; then a Nikon D5300 Digital Camera attached to portable camera tripod was used to photograph vegetation cover. The photos were taken from 1.5 m above the ground at a downward angle of 90º with the same field of view, resolution and other settings. A bubble level was used to ensure that the tripod, camera, and resulting images were exactly vertical. Subsequently, the coverage was estimated manually by using Motic Images Plus 2.0 software (Hong Kong, Asia) and converted into a percentage of the plot surface (Rożek et al. 2020).

Light intensity
Light intensity (µmol m -2 s -1 ) was measured using a PAR Quantum Sensor (Kipp & Zonen, The Netherlands). The measurements were taken in the middle of a cloudy day for 5 randomly selected locations within each plot at a height of 150 cm above the ground. Average values were treated as a single observation in subsequent analyses. The measurements were taken simultaneously during a very short time period (from 10:00 to 12:00) at all plots (Rożek et al. 2020).

Chemical analyses of soils
Soil samples were passed through a sieve (2 mm) and subjected to analyses for chemical properties (Stefanowicz et al. 2021). The pH in H 2 O (1:5; w:v) was measured with a Hach Lange HQ40D multi meter (ISO 10390:1994). The total and organic C content was analysed with a analyser RC-612 (Leco;ISO 10694:1995) and total N with a Foss Tecator 2300 Kjeltec Analyzer Unit following soil digestion in H 2 SO 4 with Kjeltabs (K 2 SO 4 + CuSO 4 · 5H 2 O) in accordance with application AN 300. Samples for analyses of other elements, i.e. total Ca, K, Mg, and P, were digested with a Foss Tecator Digestor 40 unit in hot concentrated HClO 4 . Samples intended for analyses of exchangeable Ca 2+ , K + , and Mg 2+ were shaken (Laboratory Shaker type 358S, elpan) three times for 1 h in 0.1 M BaCl 2 (PN-EN ISO 11260:2011). Concentrations of elements in the extracts were analysed using flame atomic absorption spectrometry (Varian AA280 FS), with the exception of total P, which was measured with a Hach-Lange DR 3800 spectrophotometer, using the molybdenum-vanadate method. Content of available P (Olsen P) was extracted with 0.5 M NaHCO 3 (1:10; w:v) and measured with a Dionex ICS-1100 ion chromatograph (Olsen and Sommers 1982). N-NH 4 and N-NO 3 were extracted by shaking in H 2 O (1:10; w:v) and measured using Dionex DX-100 and Dionex ICS-1100 ion chromatographs. Certified reference materials were used for quality assurance of available/exchangeable element concentration analysis: ISE-859, ISE-912, and ISE-995 (WEPAL). For calibration purposes, Six Cation Standard II and Seven Anion Standard II (Dionex) were used (Chmolowska et al. 2019;Stefanowicz et al. 2017Stefanowicz et al. , 2018Stefanowicz et al. , 2021.

Calculations and statistical analyses
First we assessed homoscedasticity through a Brown-Forsythe test. Second we carried out a one-way ANOVA with bryophyte species richness and cover being the dependent variables and the identity of tree species the categorical predictor. Given a significant ANOVA statistic, we finally carried out Tukey's HSD tests for unequal sample sizes to determine whether the differences between particular tree species were significant. Student's t-tests were performed in order to verify significant differences in the above-mentioned parameters between plots representing deciduous and coniferous trees. Prior to the analyses, the distribution normality of variables was verified by means of the Kolmogorov-Smirnov test. The nonparametric Mann-Whitney U test was applied for comparison between plots representing trees of native and alien origin.
Permutational multivariate analysis of variance (PERMANOVA) was performed to test for differences in bryophyte species composition between plots representing different tree species, deciduous and coniferous trees as well as native and alien trees (Anderson 2001). Due to unbalanced design type III Sum of Squares (SS) were used for partitioning of multivariate variation. Pair-wise comparisons among all pairs of tree species were calculated as multivariate pseudo-t statistics and p values obtained using permutation procedure. The analyses were based on the matrix of the presence/absence of bryophyte species using the Jaccard coefficient, with 999 permutations for each test. Subsequently, we evaluated which species were most responsible for differentiating bryophyte communities between plots representing deciduous and coniferous trees as well as native and alien trees, using similarity percentage routine (SIM-PER; Clarke 1993). The analysis was based on the square-rooted matrix of bryophytes frequency in plots representing individual tree species; Bray-Curtis similarity index was used.
The matrix with mean values of habitat parameters and frequency of bryophyte species in a given plot type was created (n = 14). The plot type represents the plot with a given tree species. The frequency of bryophyte species was considered percentage of plots of a given type in which a given species occurs. Firstly, non-metric multidimensional scaling (NMDS) was used to visualise the pattern of similarities between different plot types in terms of species composition. All habitat parameters were passively fitted to the ordination space. Subsequently, factor analysis (FA) based on principal component analysis (PCA) was applied to obtain uncorrelated factors representing habitat parameters. The analysis was based on the matrix with mean values of habitat parameters in a given plot type. The factors with eigenvalues > 1 were extracted according to Kaiser criterion and varimaxrotated to facilitate their interpretation. Then canonical correspondence analysis (CCA) was used to relate the frequency of bryophyte species to factors representing habitat conditions. A Monte Carlo permutation test based on 9999 random permutations was performed in order to assess the statistical significance of canonical axes (ter Braak and Smilauer 2002).
Next, factor analysis (FA) was applied on the data matrix including habitat parameters in individual sampling plots (n = 53). Then, we conducted forward stepwise multiple linear regression analysis (with a threshold of P < 0.05 to entry and remove) to investigate the effect of Factors derived from factor analysis on bryophyte species richness and cover. The analysis was based on the data matrix for individual sampling plots (n = 53). The procedure constitutes a combination of the forward selection and backward elimination. The initial models included only a regression constant and a predictor with the lowest input statistic (p-to-enter) was firstly entered into the model. Prior to the analysis, the potential collinearity of the predictors was checked by calculating the variance inflation factors (VIFs). A detailed residual analysis was performed to validate the regression model and to detect outliers. If studentised residuals greater than 3 in absolute value (corresponding to points more than 3 standard deviations from the fitted model) were detected, the outliers were excluded from the analysis. The Durbin-Watson statistic was calculated to evaluate the potential presence of a serial correlation of residuals.
The ecological indicator values were assigned to particular bryophyte species according to Ellenberg et al. (1991) and modified by Hill et al. (2007) classification. Based on bryophyte species composition, the mean indicator values were calculated for each plot. Next, associations between these values and selected habitat parameters were tested with Pearson correlation coefficients. The analysis was done to determine to what extent bryophyte species composition in the studied experimental forest captured in the mean ecological indicator values corresponds to the measured habitat variables. Following the Brown-Forsythe test to assess the equality of variances, Student's t-tests (P < 0.05) were performed in order to verify significant differences in ecological indicator values (L -light, F -moisture, R -reaction, N -nitrogen, human impact (hemeroby)) between plots representing deciduous and coniferous trees. Prior to the analyses, the distribution normality of variables was verified by means of the Kolmogorov-Smirnov test. Additionally, the ecological indicator values for forest bryophytes according to Bernhardt-Römermann et al.

Results
Effect of tree species identity on bryophyte species richness and cover Bryophyte species richness differed significantly between plots representing particular tree species (Fig. 1). The mean number of species per plot was significantly higher in Quercus robur, Picea abies, and Pinus sylvestris plots than in Q. rubra plots (ANOVA; F = 2.89; P = 0.005). The cover of bryophytes also differed between plots, being the highest in P. abies and Larix decidua plots (ANOVA; F = 7.52; P < 0.001). The plots with coniferous trees were characterised by significantly higher bryophyte species richness and cover than those with deciduous trees (Student-t tests; P < 0.05). Concerning the comparison in terms of tree species of different origin, the plots with native trees were characterised by significantly higher bryophyte cover (U-Mann Whitney test; P < 0.05). Bryophyte species richness was also 1 3 higher in these plots; however, the differences were not significant.

Bryophyte species composition
Altogether, 27 bryophyte species from 18 genera were recorded (Table S1). More specifically, one species of liverworts, 12 species representing 5 families of acrocarpous mosses and 14 species representing 6 families of pleurocarpous mosses were identified. The total number of species found in particular plots ranged from 1 (Quercus rubra) to 19 (Picea abies). As regards grouping of plots into deciduous and coniferous, from a total 27 recorded species, 19 were non-specific and occurred in both plot types; plots with deciduous trees harboured 4 exclusive species (i.e. Amblystegium serpens, Aulacomnium androgynum, Plagiothecium succulentum, and Rhytidiadelphus squarrosus), whereas plots with coniferous trees also had 4 exclusive species (i.e. Ceratodon purpureus, Dicranoweisia cirrata, Dicranum montanum, and Dicranum tauricum).
PERMANOVA results showed that species composition differed significantly between plots representing different tree species (F = 1.63, P < 0.001), between plots with deciduous and coniferous trees (F = 4.02, P < 0.001) and between plots with native and alien trees (F = 2.44, P < 0.001). Pairwise comparisons between plots representing different tree species showed that most of significant differences were recorded between deciduous and coniferous tree species (Table 1). NMDS ordination Fig. 1 Bryophyte species richness and cover (mean ± SE) in the studied plots in relation to: particular tree species (A), deciduous and coniferous trees (B), native and alien tree species (C). The results of one-way ANOVA (A), Student-t test (B), and U-Mann Whitney test (C) are provided. The various letters above the whiskers indicate statistically significant differences according to Tukey's HSD test (P < 0.05) diagram also clearly separated the plots with coniferous trees from plots representing deciduous trees (Fig. S1). The differences were the most pronounced for the coniferous trees P. menziesii and L. decidua that differed significantly from 7 and 4 different deciduous tree species, respectively ( Table 1). The plots representing different coniferous tree species did not differ significantly between each other and the average similarity of bryophyte species composition was always greater than 20%. Regarding differences within deciduous tree species, Q. robur differed significantly from A. platanoides and Q. rubra plots.
Based on the SIMPER analysis we identified which bryophyte species contributed disportionally to the differentiation of deciduous and coniferous tree plots as well as native and alien tree plots (Table S2). Sixteen bryophyte species proved to more abundant in plots with coniferous trees, whereas only eleven proved to be associated to a greater extent with plots with deciduous trees. Only nine species turned out to be more abundant in plots with alien tree species.

Relationships between habitat properties and bryophytes
The relationships between bryophyte species composition and factors representing habitat conditions were determined by means of CCA (Fig. 2). The first axis and all canonical axes taken together were statistically significant (F = 2.0, P = 0.020 and F = 1.5, P = 0.017, respectively). The ordination diagram distinguished plots with coniferous trees (the lower left side of the graph), with the exception of Pinus nigra plots, from plots representing deciduous trees. Moving from the left to the right side of the diagram, soil pH and macronutrient concentrations increase. Plots with Acer spp. and Tilia cordata were characterised by the highest K, Mg, and Ca concentrations in soil and the highest soil pH; the occurrence of four bryophyte species, i.e. Rhytidiadelphus squarrosus, Plagiomnium affine, Atrichum undulatum, and Amblystegium serpens, was most strongly associated with these plots. The second clear association occurs in the lower left part of the graph, where species preferring high C, P and N-NO 3 contents in soil (high values of These species were present in most plot types (at least in 11 different plot types) and could be considered as indifferent with respect to measured habitat parameters and tree species. Based on the whole data matrix, factor analysis reduced 16 variables to four factors with eigenvalues > 1 that jointly explained 76.66% of the total variation (Table 2). Factor 1 was associated with soil pH and essential macronutrient levels (total and exchangeable K, Mg, Ca). Factor 2 was related to soil fertility and referred to contents of total and organic carbon and total and NO 3 -form nitrogen. Factor 3 was negatively associated with vascular plant species  richness and cover. Light was highly negatively correlated with Factor 4, whereas Olsen P positively.
The results of the multiple stepwise regression analysis are presented in Table 3. A forward stepwise procedure with four Factors derived from factor analysis as predictor factors and bryophyte species richness as the dependent variable revealed that only Factors 3 and 4 were included in the model (F = 10.38; P < 0.001). Both Factors were negatively related to bryophyte species richness. This means that vascular plants and light positively influenced bryophyte species richness whereas high concentrations of P negatively affected the number of bryophyte species. As regards bryophyte cover, Factors 1 and 3 were included in the model (F = 11.63; P < 0.001). As for species richness, vascular plant species richness and cover were positively related to bryophyte cover; whereas Factor 1, related to soil pH and essential macronutrient concentrations, showed a negative effect.

Bryophyte indicator values
The results showed a significant positive correlation (P < 0.05) between mean indicator value R (reaction) and soil pH (Table 4; Fig. S2). Moreover, indicator value R was also significantly related to both total and exchangeable Ca concentrations in soil as well as macronutrient concentrations (K and Mg). Both nitrate form of nitrogen (N-NO 3 ) and ammonium (N-NH 4 ) were significantly and positively correlated with mean indicator value N. Soil pH and concentrations of Ca, K and Mg and carbon content were also positively related to mean indicator value N. A significant relation was not recorded for mean indicator value L (light) and neither light nor other soil parameters. light (-0.72), Olsen P (0.81) 6.61 Table 3 Result of forward stepwise multiple regression analysis for the effect of Factors related to habitat parameters derived from factor analysis on bryophyte species richness and cover (for Factor characteristics see Table 2). Standardized β coefficients follow from standardization of variables to a mean of 0 and a standard deviation of 1. Concerning the comparison of ecological indicator values between plots with deciduous and coniferous trees, the latter ones were characterised by significantly lower F (moisture) and R (reaction) values (Student-t tests; P < 0.05). The N (nitrogen) values were also lower in plots with coniferous trees; however, the differences were not significant. For both deciduous and coniferous plots, L (light) values were highly variable and the mean L was nearly identical in plots with deciduous and coniferous trees. The indicator values of human impact (hemeroby) also did not differ significantly between these two plot types (Fig. S3).
Based on the overall occurrence frequency of bryophyte species, the comparison of the proportions of bryophytes representing ecological indicator values according to Bernhardt-Römermann et al. (2018) between plots with deciduous and coniferous trees gave consistent results. The greater frequencies of bryophytes preferring lower pH, drier places and less nutrients in soil were found in plots with coniferous trees. Moreover, in these plots a greater proportion of species largely restricted to closed forests were noted compared to plots with deciduous trees (Fig. S4).

Discussion
In the present study we assessed for the first time the impact of 14 tree species growing in monoculture plots and that of related habitat parameters on bryophyte species composition, richness and abundance. Our study showed that both overstorey tree species and various habitat parameters including certain soil factors, light conditions and vascular plant layer had a great impact on ground bryophyte communities that emerged in the experimental forest 50 years after its establishment.

Effect of light conditions
Light intensity had a positive impact on bryophyte species richness. Light availability affects bryophyte species composition (e.g. Humphrey et al. 2002;Tinya et al. 2009) due to the different light requirements of particular species. Forest floor bryophytes have relatively low light demands. They have an evolutionarily set slower rate of photosynthesis and also due to the fact that they lack stomata and cuticle, they cannot directly control water loss and thus they are dependent on the water available in the environment (Proctor 2008). The effect of canopy tightness and the resulting shading may contribute to retaining moisture, but on the other hand, also decrease the amount of water reaching the forest floor. Moreover, light reaching the ground can directly affect soil temperature and thereby increase evaporation and decrease humidity on the microhabitat level (von Arx et al. 2012). Probably due to these complex relationships, there is no consensus about the definite effect of light availability on the ground bryophyte layer in the literature. Some studies suggested no significant relationship between light conditions and bryophyte species richness (e.g. Humphrey et al. 2002;Mills and MacDonald 2004). Tinya et al. (2009) found a significant relation between total cover of ground floor bryophytes and light, while species richness was not dependent on light. On the other hand, Moora et al. (2007) found that bryophyte species richness was positively related to local light availability. We found that light intensity enhanced bryophyte richness and significantly affects  Ellenberg et al. (1991) and modified by Hill et al. (2007)  1 3 species composition. We believe that it is directly related to the effect of tree species. The highest bryophyte species richness were observed in plots with conifers, such as Pinus spp., Larix decidua and Picea abies, as well as from Quercus robur plots. The highest light intensities were also found under these species of coniferous trees (see also Rożek et al. 2020).
In contrast, plots with Abies alba, Tilia cordata, Q. rubra, Fagus sylvatica, and Acer spp. were characterised by high shading conditions which could limit the growth and occurrence of certain light-requiring species, such as Dicranum scoparium, Polytrichum juniperinum, Scleropodium purum, and Ceratodon purpureus (see Table S1). This suggests that light availability might decrease below a critical threshold for some bryophytes, which results in a lower species richness in the most shaded plots. Interestingly, the mean L (light) indicator values based on species composition did not differ significantly between deciduous and coniferous plots. It is worth noting, however, that in coniferous plots where the light intensity was the highest (see also Rożek et al. 2020), the values of the F (moisture) indicator were much lower than in deciduous plots. Such a result suggests that habitat moisture, which indirectly results from light conditions in the forest floor, is a more important factor determining the occurrence of particular species. A second possible explanation is related to the previous observations made by Ewald (2009), who found that Ellenberg indicator values for ground bryophytes did not significantly predict measured habitat parameters in mountain forests and bryophyte indication proved to be particularly poor for light. Such a result could be also caused by a high proportion of indifferent species with a wide range of tolerance to light conditions. Finally, one should keep in mind that the current bryophyte species composition may result more from past light regimes than from the present light conditions, a phenomenon associated with the slow change of species pool similarly as it has been suggested for understory vascular plant communities (Thomas et al. 1999).

Importance of soil chemical parameters
Soil pH and content of macronutrients in soil proved to affect bryophyte cover. This effect seems to be directly related to the influence of tree species since the plots of particular species differed significantly in terms of soil pH and base cation concentrations (Stefanowicz et al. 2021). More acidic soil pH and lower concentrations of earth-alkaline cation concentrations are often found under coniferous trees compared to deciduous ones (Augusto et al. 2002(Augusto et al. , 2003Nordén 1994). The same concerns the studied experimental forest, in which plots with conifers had significantly lower soil pH (Rożek et al. 2020;Stefanowicz et al. 2021). This is reflected in the increased cover of bryophytes in plots with coniferous trees. However, this also depends on particular coniferous tree species. The largest coverage of bryophytes was recorded in the plots with Pinus spp., Picea abies, and Larix decidua; whereas plots with Pseudotsuga menziesii and Abies alba were characterised by relatively low bryophyte coverage, which was comparable to those recorded in plots with deciduous trees. This is in line with the acidifying capabilities of coniferous species reported by Augusto et al. (2003) who found that P. sylvestris and P. abies have the greatest soil acidifying capabilities; whereas P. menziesii and A. alba had moderate effect on soil acidity. The shady conditions under the latter two species may additionally limit the coverage of bryophytes. Plots of Acer spp. and Tilia cordata trees were characterised by the highest soil pH and the concentrations of base cations (Stefanowicz et al. 2021). These tree species have been recognised to have low soil acidifying ability (Augusto et al. 2003), and bryophyte cover was relatively low under these deciduous trees. The cover of bryophytes may be limited in these plots, firstly, by a lower abundance or lack of species confined to extremely acidic and acidic substrates (e.g. Pleurozium schreberi, Polytrichum juniperinum, Dicranum scoparium), and, secondly, by a large amount of litter constituting a mechanical barrier that limits the amount of soil substrate available for bryophyte development. Although we recorded species that previously were found to be associated with leaf litter, e.g. Brachythecium rutabulum (see Bates 2008), this species was frequently noted in plots with both deciduous and coniferous trees. Apart from the distinct effect of soil pH on the coverage of bryophytes, the species composition also clearly responded to this factor. This is demonstrated by significant differences in species composition between plots with deciduous and coniferous trees, significantly lower R (reaction) indicator values for plots with coniferous trees and most importantly by a positive correlation of mean R (reaction) indicator values with soil pH. The higher pH of the forest bottom under coniferous trees is certainly also related to the rate of nutrient decomposition in the soil and the consequent amount of humus. A lower proportion of nutrients reduces the diversity of vascular plants in the herb layer, which do not compete significantly with bryophytes and results in a higher proportion of bryophyte species in these areas. We did not find a significant impact of organic carbon and nitrogen contents in soil on bryophyte diversity and abundance. This may be due to relatively slight differences in the content of these elements between plots with different tree species (Rożek et al. 2020;Stefanowicz et al. 2021) that would potentially result in changes of total species number. On the other hand, we found a significant correlation between organic C, NH 4 and NO 3 forms of nitrogen concentrations in soil and mean N (nitrogen) indicator value based on bryophyte species composition in plots. This shows that species richness and abundance of bryophyte cover may not be a sufficient determinant of the effect of soil fertility conditions, while changes in species composition are of key importance. This is a known phenomenon previously reported in the case of vascular plants in forests, in which gradual replacement of oligotrophic species by eutrophic species as a response to eutrophication was observed with a simultaneous absence of significant changes in species richness (Dirnböck et al. 2014).

Interaction between ground bryophytes and understory vascular plant layer
An interesting result concerns the positive relationship between the diversity/cover of bryophytes and vascular plants inhabiting forest floor. Many studies revealed that vascular plants negatively affect the diversity and abundance of bryophytes (e.g. Bergamini et al. 2001;Löbel et al. 2006;Virtanen et al. 2000;Yan et al. 2013). Such an effect was most frequently explained by competition between these two groups. However, most studies concerned grassland communities in which the herbaceous plant layer is, as a rule, significantly more developed than in the temperate forests. As regards forests, Fojcik et al. (2019) also revealed the negative impact of herbaceous plants on the total coverage of bryophytes in a managed temperate mixed forest in Poland. Similarly, Turkington et al. (1998) showed that the addition of NPK fertilizer to a nutrient-poor boreal forest ecosystem resulted in an intensive growth of vascular plants, mainly grasses, which at the same time caused a decrease in bryophyte cover. Contrastingly, in our study, we observed a positive relationship between vascular plant cover and bryophyte cover. This fact can be explained in two different ways. Firstly, vascular plant cover was not dominated by expansive grasses and sedges that could easily eliminate bryophytes from the forest floor (Chmura and Sierka 2007), and vascular plant cover was relatively low in the plots (see Rożek et al. 2020). This, in turn, is in line with the observations of Ingerpuu et al. (2005), who showed that even in grassland communities, positive interactions between bryophytes and vascular plants prevail when density of plant cover is low. Secondly, the plots with coniferous trees were richer in both bryophytes and vascular plants, while deciduous plots were characterised by a thick leaf litter cover acting as a mechanical and/or chemical barrier that limits the development of both forest-floor bryophytes and vascular plant species (Evans et al. 2012;Startsev et al. 2008). Soils in the plots with coniferous trees were characterised by lower macronutrient concentrations compared to deciduous ones, and this could indicate that asymmetry of competition between bryophytes and vascular plants is rather low as it is under nutrient limitation conditions (Rydin 2008). The higher cover of bryophytes in the plots with coniferous trees could also be associated with generally lower soil fertility since the total biomass of mosses, as a rule, decreases with increasing fertility (Mäkipää et al. 2000). Further to the above, vascular plants can have a positive effect on bryophytes by providing better microclimatic conditions for their growth. It has been suggested that this mechanism occurs in temperate forests in Hungary (Márialigeti et al. , 2016.

Comparison of effects between deciduous/coniferous and native/alien trees
The higher species richness and cover of bryophytes under coniferous trees compared to deciduous trees supports that in forests dominated by conifers, bryophytes are an important component in terms of the biomass they provide (Jiang et al. 2018;Mäkipää et al. 2000), in contrast to deciduous forests, where forest floor bryophytes are often of little significance (Startsev et al. 2008). The higher bryophyte coverage in coniferous plots can certainly be related to the more favourable understorey microhabitat conditions 1 3 under conifers (Bartemucci et al. 2006), while litter from deciduous trees proved to have a physical and/ or chemical inhibitory effect Startsev et al. 2008). Differences in bryophyte communities between conifer-dominated and deciduousdominated stands clearly confirm the occurrence of this phenomenon (e.g. Bartels et al. 2018). Nevertheless, a given tree species is of great importance in shaping the specific microhabitat. For example, in the case of bryophyte coverage, conifers differ from each other in their effect since we recorded significantly lower cover of bryophytes in Abies alba and Pseudotsuga menziesii plots compared to Picea abies and Larix decidua plots. Species richness within plots of deciduous trees also depends on tree species. Quercus robur was clearly the leader, both in terms of the mean species richness per plot and the total number of recorded species, which were comparable to the most species-rich plots of conifers.
Two main tree species that are involved in forming a stand consistent with the natural mixed oak-pine forest (Querco roboris-Pinetum), i.e. Quercus robur and Pinus sylvestris, provided the greatest species richness of ground bryophytes. Furthermore, bryophyte composition in the plots representing these two tree species is relatively diverse. As many as 10 species were exclusive for Q. robur plots, 3 species for P. sylvestris plots, and 8 were common to both plot types. Altogether this accounts for 77.8% of all species recorded in the experimental forest. This shows that the promotion in forestry of tree species forming plant communities consistent with the natural habitat is highly beneficial for the biodiversity of ground bryophytes.
Our results showed that, in general, alien to Poland tree species (i.e. Quercus rubra, Pseudotsuga menziesii, and Pinus nigra) negatively affected bryophyte communities. The strongest adverse effect on diversity and abundance of bryophytes was recorded in Q. rubra plots. It is known that this invasive tree has a negative effect on native plant species richness and abundance both in old-growth and secondary forests (Chmura 2013;Woziwoda et al. 2014). In addition, Q. rubra has a strong influence on soil properties (Ferré and Comolli 2020). With respect to the impact of Q. rubra on bryophytes, the literature reports are rather poor and divergent. Woziwoda et al. (2017) found a high number of native epiphytic bryophytes on the bark of Q. rubra and concluded that this tree species could successfully provide a functional alternative to oaks native to Poland (Q. robur and Q. petraea). On the other hand, Jagodziński et al. (2018) reported fewer epiphytes in Q. rubra stands than in Q. robur in a reclaimed lignite mine spoil heap. As regards ground bryophytes, a report concerning a disturbed habitat showed that Q. rubra did not have a negative impact on ground bryophytes and the number of species was similar to that under the Q. robur canopy . Contrary to this, we found a significantly lower number of bryophyte species in Q. rubra compared to Q. robur plots. We can explain the negative influence of Q. rubra by the limiting effect of thick leaf litter accumulation on the forest floor as well as the high shading generated by this species. The leaf litter of Q. rubra decomposes very slowly (Dobrylovska 2001), and involves high production of barely decomposable material which accumulates on the forest floor and which is probably responsible for the sparse bryophyte cover. A similar observation was made by Woziwoda et al. (2017) in Pinus sylvestris-Quercus rubra secondary forests. In our opinion, the adverse effect of Q. rubra should be considered strong because after 50 years in the whole experimental forest, we recorded as many as 18 different bryophyte species in Q. robur plots and only one species, i.e. Hypnum cupressiforme, in Q. rubra plots. As regards the second alien species, i.e. P. menziesii, the negative effect concerns only the cover of bryophytes compared to native coniferous trees such as Picea abies and Larix decidua; mean species richness of bryophytes did not differ significantly from the remaining coniferous and deciduous trees. Finally, we did not observe significant differences in either species richness or cover of bryophytes between P. nigra (alien) and P. sylvestris (native) plots, although the species composition was slightly different between plots of these two tree species.

Effect of tree species identity on bryophytes
Despite the fact that different tree species grow in the immediate vicinity in one experimental forest, we have observed a great effect of the tree species on bryophyte species composition, richness and abundance of ground bryophytes. Therefore we can assume that the changes in dominant tree species in a forest stand may alter many dependent bryophyte species. This can be extremely important in the context of climate 1 3 changes that may affect bryophytes directly, but also indirectly through changes in the stand species composition. Thus the range contraction of certain forest tree species would have also serious consequences for ground bryophytes. With regard to temperate forests coniferous species, such as L. decidua, P. abies, P. sylvestris, and pioneer deciduous trees, e.g. B. pendula, are predicted to have a vast range contraction in central and eastern Europe (Schueler et al. 2014;Dyderski et al. 2018). Since we recorded both the greatest bryophyte species richness and abundance under these tree species, parallel reduction of ground bryophyte diversity can be expected. A great species richness of bryophytes was also associated with Q. robur trees, which also proved to be a refuge for rare epiphytic bryophytes and lichens (Király and Ódor 2010;Kubiak and Osyczka 2017); nevertheless although Q. robur tends to gain in range as the climate warms, it will lose much of its currently occupied range . Consequently, in addition to bioclimatic factors, which were recognised as principal drivers of bryophyte diversity, causing changes in their distribution ranges (He et al. 2016), it is also worth considering the concurrent effect of tree species identity in future research.
The identity of the tree species also had a pronounced influence on bryophyte species composition. We identified as many as 15 pairs of tree species that differed significantly in bryophyte species composition. This indicate that the process occurring in the structure of bryophyte communities, including the replacement of certain species by others, may be of great importance especially when species richness itself would be a poor indicator of the ground bryophyte conversion induced by changes in stand composition. The abundance of bryophytes associated mostly with coniferous trees (cf Fig. 2) may also experience a decline as a consequence of regression of certain coniferous species from forests in the temperate zone. Consequently, even if a local diversity will not decrease significantly, strong changes in ground bryophyte composition will probably occur.
Two species of acrocarpous mosses considered as expansive were also recorded, i.e. Dicranum tauricum and Dicranoweisia cirrata. They were recorded only in a single plots of Pinus sylvestris and Pinus nigra. The fact that the species were not found in any other surrounding study plots indicates that they are at the beginning of their colonisation of this area and it can be assumed that they have settled the study area only recently. The spread of these species and the range of occupied habitats could increase with time after initial colonisation, depending on an invader's ability to establish in new communities (Söderström 1992). It is also surprising that both species, which almost always colonise the bark of trees, were not observed growing epiphytically there nor on the rest of the fallen branches or bark, but directly on the forest litter. The increase of abundance in Central-East Europe, such as in Poland, has only been reported in last two decades (Stebel and Plášek 2001;Stebel et al. 2012). Such expansive species are increasingly viewed as a significant component of global change and in many cases one of the major drivers of current biodiversity loss (Didham et al. 2007). Nevertheless, in the context of global climate changes, even highly dispersive organisms like expansive bryophytes could be not equipped to fully track the rates of ongoing changes in the course of the next decades (Zanatta et al. 2020).

Implications for management
Forest integrity is considered to be important conservation point that guarantees desirable characteristics such as natural biodiversity, stand structure and continuity (Frego 2007). Although its defining criteria are still under discussion, bryophytes are proposed as important components of forest integrity, because they play a vital role in soil development, nutrient biogeochemical cycling, and ecological succession and indicate changes associated with specific forest processes. Our results provide data that can be useful in reasonably promoted forest regeneration treatments aimed at biodiversity conservation both in degraded forests and after clear-cuts. The inclusion of certain species into the stand may have a beneficial effect on bryophyte communities. Firstly, by means of changing the microhabitat conditions, which could be more favourable for bryophyte species occurrence and/or provide a heterogeneity of habitat factors that would promote the appearance of more diverse assemblages. Secondly, we indicated that bryophyte communities regenerate much faster under certain tree species after clear-cutting. Some trees ensure higher bryophyte species richness and faster entry of species characteristic for a given plant community. Providing suitable conditions for bryophytes to colonise and establish in young forest stands after logging is a key factor 1 3 for bryophyte diversity and species composition in mature forests. We showed that the preference for tree species forming plant communities inherent to the natural habitat in forestry is highly beneficial for the biodiversity of ground bryophytes. Consequently, we support the idea of maintaining the selection of Quercus robur and Pinus sylvestris for plantings in managed forests planted on habitats typical for mixed coniferous forests (Querco roboris-Pinetum), which constitute a potential natural vegetation in the area of 13.64% of Poland (Matuszkiewicz 2008). This is also in agreement with the currently promoted forest management strategies aimed at diversification of tree species in order to increase and/or maintain a high level of understorey biodiversity.

Conclusions
The experimental forest founded after clearing, stump removal, and ploughing gave us an opportunity to trace the colonisation processes of ground bryophytes from nearby mixed coniferous forest communities over a period of 50 years. Due to the persistent close connection of the experimental plots with the neighbouring mature forests, it is now possible to compare formed bryophyte assemblages and to estimate the impact of 14 tree species on them. Our study showed that overstorey tree species had a great impact on bryophyte communities. Therefore, we can assume that changes in the dominant tree species in the stand may result in significant changes in ground bryophyte communities. The effect of tree species identity on bryophyte species composition seems to be even greater than on species richness. This indicate that the process occurring in the structure of bryophyte communities may be of great importance especially when species richness itself would be a poor indicator of the ground bryophyte conversion induced by changes in stand composition. Moreover, the related soil chemistry and other habitat parameters, including light conditions and the vascular plant layer, also proved to be of great importance to the richness and cover of bryophytes in the forest floor. Nevertheless, the effect of tree species and habitat parameters is a network of mutual connections which, only when interpreted jointly, can explain the changes taking place in bryophyte communities. Generally, we found significantly higher bryophyte species richness and cover in the plots representing coniferous trees. Soil pH, concentration of macronutrients (Ca, K, and Mg) in soils, vascular plant layer and light conditions, proved to be the key habitat factors influencing the occurrence of ground bryophytes, and the effects of these parameters are related to the tree species effect. We also indicated that the introduction of alien tree species, in particular Q. rubra, may have an adverse effect on bryophyte communities. Finally, we conclude that the selection of tree species (i.e. Q. robur and P. sylvestris) that form plant communities consistent with the natural habitat (Querco roboris-Pinetum) in forestry is highly beneficial for maintaining ground bryophyte biodiversity.

Data availability
The data generated or analysed during this study are included in this published article (and its supplementary information files). Additional data are available from the corresponding author on request.

Declarations
Conflicts of interest/Competing interests The authors have no conflicts of interest to declare that are relevant to the content of this article.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.