Variation in floristic and trait composition along environmental gradients in the herb layer of temperate forests in the transition zone between Central and SE Europe

Species- and trait-environment linkages in forest plant communities continue to be a frequent topic in ecological research. We studied the dependence of floristic and functional trait composition on environmental factors, namely local soil properties, overstory characteristics, climatic parameters and other abiotic and biotic variables. The study area comprised 50 monitoring plots across Slovenia, belonging to the EU ICP Forests monitoring network. Vegetation was surveyed in accordance with harmonized protocols, and environmental variables were either measured or estimated during vegetation sampling. Significant predictors of species composition were identified by canonical correspondence analysis. Correlations between plant traits, i.e. plant growth habit, life form, flowering features and CSR signature, were examined with fourth-corner analysis and linear regressions. Our results show that variation in floristic composition was mainly explained by climatic parameters (mean annual temperature, mean annual precipitation), soil properties (pH) and tree layer-dependent light conditions. Trait composition was most closely related with tree layer characteristics, such as shade-casting ability (SCA, a proxy for light availability in the understory layer), tree species richness and tree species composition. Amongst soil properties, total nitrogen content and soil texture (proportion of clay) were most frequently correlated with different species traits or trait states. The CSR signature of herb communities was associated with tree layer SCA, soil pH and mean annual temperature. The floristic composition of the studied herb-layer vegetation depended on temperature and precipitation, which are likely to be influenced by ongoing climate change (warming and drying). Trait composition exhibited significant links to tree layer characteristics and soil conditions, which are in turn directly modified by forest management interventions.


Introduction
Herb-layer vegetation harbours the highest diversity of vascular plant species in temperate forests (Gilliam 2007;Thrippleton et al. 2016). Although not directly targeted during forest management interventions, this forest stratum plays a key role in the functioning of forest ecosystems . It may serve as a good ecological indicator and its species composition results from a complex interplay between numerous environmental drivers and spatiotemporal resource gradients (McEwan and Muller 2011;Sipoš et al. 2020). Identifying the ecological determinants of community composition has become a crucial step towards gaining a better understanding of how forest vegetation will respond to the challenges imposed by global environmental change, e.g. increasing CO 2 levels, nutrient deposition leading to either eutrophication or acidification and altered disturbance regimes (Suding et al. 2008;Amatangelo et al. 2014;Dalle Fratte et al. 2019).
Temperate forest ecosystems currently cover around 16% of the global forest area (Hansen et al. 2010). They account for a wide variety of forest types, including thermophilous deciduous forests in Mediterranean subzones (Franklin 1988). Temperate forest understories are inhabited by plant species with distinctive life-history strategies that manifest as functional traits adapted to local environmental conditions (Durak and Durak 2021). Forest understory plant species vary in growth and life form, leaf lifespan, phenology, the prevalence of clonal reproduction, successional status and the mechanisms of seed dispersal (Burton et al. 2020). In mature forest stands the closed tree layer canopy induces a strong filtering effect by limiting light availability at the forest floor. Species composition and diversity and the cover of herb-layer communities might be additionally filtered by other environmental stressors, such as low nutrient availability on silicate bedrock or a water deficit on drier sites (lower soil moisture and air humidity) (Catorci et al. 2011). Similarly, coldinduced stress at higher altitudes favours species with a notable degree of conservatism, e.g. smaller leaves and a shorter, rosette-forming stature. Survival, growth and reproduction in resource-limited conditions can be achieved by a broad spectrum of ecological strategies, spanning from long-living organisms (trees) to annual herbaceous species (Silvertown et al. 1993).
Plant functional traits are tightly linked to environmental factors (Terwei et al. 2016;Burton et al. 2017). In fact, species traits can inform us about the leading dimensions of morphological, phenological and physiological adaptations that enable species to persist in the forest herb layer (Candeias and Fraterrigo 2020).
For example, many typical forest species have evolutionarily developed the specific strategy of finishing their life cycle before canopy closure in the deciduous overstory (e.g. early spring ephemerals) (Muller 1978;Sercu et al. 2021). Other forest specialists (mainly late-successional, shade-tolerant herbaceous perennials confined to closed forest interiors and rarely found in other habitats, such as forest edges or canopy gaps; Brunet et al. 2011;Schmidt et al. 2014) start to grow and flower later in the season. Such temporal niche differentiation importantly contributes to species coexistence in forest understories. Forest plants can tolerate a low light environment and exhibit a distinct functional profile, e.g. higher specific leaf area, larger seeds and a variety of storage or clonal organs (Bierzychudek 1982;Durak and Durak 2021). Furthermore, according to the study of Marinšek et al. (2015), habitat specialists in mixed broadleaved deciduous forests have proven to thrive primarily at high altitudes and on colder, drier sites, and the age of first flowering is later than that of generalists. They also have a smaller seed diameter, lower leaf dry matter content, lower mean canopy height and higher specific leaf area than generalists (Marinšek et al. 2015). In terms of Grime's CSR theory (Grime 1977), late-successional forest understories are characterized by stress-tolerant or competitive species, whilst disturbances increase the ruderal component of the community-level CSR signature (Eler et al. 2018).
In this study, we set out to determine (1) which environmental variables most significantly impact the floristic composition of the herb layer and (2) how species' life-history traits and community CSR signature are related to environmental variables. Using systematically sampled vegetation data distributed throughout Slovenia, we hypothesized that the floristic and functional composition of herb-layer vegetation are linked to different environmental factors.

Study area and vegetation sampling
The study was conducted in 50 different forest sites across Slovenia. Although Slovenia is a small country in terms of geographic extent, it exhibits very diverse natural conditions. The selected forest sites differed significantly in terms of geological bedrock, climate, soil type, tree species composition, stand structure and other abiotic and biotic factors. The study area extends from lowland forests below 200 m a.s.l. to highaltitude forests reaching almost 1,500 m a.s.l. The study area is characterized by a broad climatic gradient. According to the WorldClim database (Fick and Hijmans 2017), for the period 1970-2000 mean annual temperature ranged from 3.2 to 11.7°C and mean annual precipitation ranged from 791 to 2499 mm.
We present a brief overview of the main forest vegetation types in Slovenia (Fig. 1), which is geographically and ecologically in the transition zone between Central and SE Europe. Located at the intersection between the Mediterranean, Alpine and continental climatic regimes, Slovenia exhibits great climate variability and can be divided into six phytogeographic regions (Wraber 1969): Alpine, Pre-Alpine, Dinaric, Pre-Dinaric, Sub-Mediterranean and Sub-Pannonian. The riparian forests, primarily comprising lowland woodlands located within the scope of larger streams and bodies of standing water, are mainly dominated by willows (Salix sp.), alders (Alnus glutinosa, A. incana), ashes (Fraxinus excelsior, F. angustifolia) and pedunculate oak (Quercus robur). In the hilly/submontane areas above floodplain and riverine forests, mixed forests of sessile oak (Quercus petraea) and European hornbeam (Carpinus betulus) are the predominant forest type. Most of the mid-altitude mountain areas are covered by forests formed by European beech (Fagus sylvatica), with admixtures of noble broadleaves (e.g. Acer pseudoplatanus, Fraxinus excelsior, Ulmus glabra, Prunus avium) and conifers (Abies alba, Picea abies). In the Alpine region (which encompasses the Julian Alps, Karavanke Mountains, Kamnik-Savinja Alps and Pohorje Mountains), various European beech forests mixed with Norway spruce (Picea abies), silver fir (Abies alba) and European larch (Larix decidua) reach the timberline up to the belt of the dwarf mountain pine (Pinus mugo). Scots pine (Pinus sylvestris) forests can be found throughout the country on acidic and nutrient-poor soils, but also on mixed and carbonate (dolomite) bedrock with shallow soils. Minor areas of Austrian pine (Pinus nigra) forests grow on extreme sites with a warmer microclimate. The Sub-Mediterranean region (positioned in the southwestern and western part of Slovenia with a Sub-Mediterranean climate) is covered with forest and shrubby vegetation of thermophilous broadleaf species (e.g. Ostrya carpinifolia, Fraxinus ornus, Sorbus aria, Quercus cerris, Q. pubescens). A similar type of forest is also present all over the country on sunexposed, south-facing slopes with predominant limestone and dolomite bedrock. The described natural forest types are intersected by numerous forms of secondary and other forest communities (often characterized by even-aged stands with a homogeneous structure), which are a consequence of managementmoderated tree species composition and related abiotic (edaphic) conditions (Kutnar et al. 2012;Bončina et al. 2021).
For the study we chose forest stands where a monitoring plan was established as a part of the ICP network (De Vries et al. 2000). Thirty-nine plots belong to the Level I scheme and 11 plots are designed for the Intensive Monitoring Programme (Level II). Level I plots are systematically distributed in a 16 km 9 16 km grid across Slovenia (Kutnar 2011). Level II plots were established in 11 different locations, representative of the heterogeneity of Slovenian forest types (Urbančič et al. 2016;Kermavnar and Kutnar 2020). We collected plant community and environmental data in 50 monitoring plots of 400 m 2 in size. All vascular plants in the tree (woody species with height [ 5 m), shrub (woody species taller than 0.5 m and shorter than 5 m) and herb layer (all herbaceous species and woody species with height \ 0.5 m) were recorded, following harmonized ICP protocols for ground vegetation sampling (Canullo et al. 2016). The abundance of each plant species was estimated based on a modification of Barkman's method (Barkman et al. 1964). Detailed information on the study area and the sampling design can be found in Kermavnar et al. (2021).

Environmental data
For each study plot, we collected data on soil properties, forest stand characteristics, climatic parameters and other environmental variables. Soils were sampled in each monitoring plot, according to the ICP methodology (Cools and De Vos 2016). Soil samples were taken for different layers distinguishing between organic soil layers (O f , O l , O h ) and mineral soil layers (0-5 cm, 5-10 cm, 10-20 cm, 20-40 cm and 40-80 cm). The following ecologically relevant variables were measured for soil samples from different layers (Table 1): pH (measured in CaCl 2 ), soil organic carbon (SOC, g/kg), total nitrogen content (g/kg), proportion of clay (derived from soil texture analysis) and soil moisture content (%) (Kermavnar et al. 2021). The values were then averaged across layers to obtain single plot-level values. Tree layer characteristics were also considered: tree species richness, the proportion of broadleaves and tree layer shade-casting ability (SCA), which was calculated as weighted averages, using tree species relative abundances and shade production values (evaluated based on the leaf area index) provided by Leuschner and Meier (2018). Low SCA values indicate relatively sparse canopy (e.g. Pinus sylvestris-dominated stands), whilst high SCA values denote dense canopy casting deep shade (e.g. Fagus sylvatica or Carpinus betulus-dominated stands). Two climatic parameters, mean annual temperature (MAT) and mean annual precipitation (MAP) were sourced from the World-Clim database (Fick and Hijmans 2017). The environmental matrix also contained data regarding plot altitude, heat load index (Hi) expressing the relative amount of heat received through solar radiation considering local slope aspect and slope steepness (Geiger 1966), the cover of woody debris (%) and the cover of the moss layer (%) on developed soil -both were visually estimated during the vegetation surveys. Fourteen environmental variables in total were obtained for each plot either from data collected in the field or from external sources (Table 1). Some additional information on the acquisition and/or calculation of environmental data is given in Kermavnar et al. (2021). The correlation matrix for explanatory variables can be found in the Supplementary materials (Appendix 1).

Functional traits data
Life-history traits for the recorded plant species were taken from the national flora (Martinčič et al. 2007). Data on plant growth habit (PGH; distinguishing between six growth habits using our own expert-based botanical knowledge derived from Martinčič et al. 2007), plant life form (PLF; Raunkiaer 1934) and the start, end and duration of flowering phenology (Martinčič et al. 2007) were collected. In addition, each species was characterized in terms of CSR ecological strategies, using the approach outlined by Pierce et al. (2017). Data on the C-, S-and R-score were taken from the PLADIAS open-source database (https://pladias.cz/en/; Chytrý et al. 2021). Information compiled on selected plant traits is presented in Table 2.

Data analysis
Variation in floristic composition between monitoring plots was examined with canonical correspondence analysis (CCA; ter Braak 1987). CCA is a unimodal ordination technique with the ordination axes constrained as linear combinations of environmental predictor variables (ter Braak 1987). We applied a manual forward selection procedure to select the main environmental variables, beginning with an interceptonly model. Model choice was made based on Akaike's information criterion as a measure of model quality. The significance of the final suite of constraints was determined using an ANOVA-like permutation test with 999 permutations. Interdependence between constraining variables was checked with a variance inflation factor. Species abundances were ln(x ? 1) transformed, and rare species (present only in one plot) with minor abundance were excluded beforehand. To increase the interpretability of CCA outputs, all plots/sites were classified based on the tree layer species composition and parent material. First, we analysed the percentage cover of each tree species occurring in the tree layer. Species were then ranked with respect to their abundances, and the ratio between broadleaves and conifers was calculated. If the plot was dominated by a single tree species, it was classified as a monodominant stand (i.e. Fagus sylvatica and Picea abies in our dataset). If the tree layer was in favour of broadleaves other than Fagus sylvatica, it was classified into the group of other broadleaves. If the tree layer was in favour of conifers other than Picea abies, it was classified into the group of other conifers. Within this scheme, we also included geological criteria for plots dominated by Fagus sylvatica, as these communities depicted a wide forest soil and geological spectra. Furthermore, each plot was also classified based on the phytosociological typology of forest site types (Kutnar et al. 2012;Bončina et al. 2021). As a result, we distinguished between five different site groups: Fagus sylvaticadominated sites on carbonate bedrock (12 plots), Fagus sylvatica-dominated sites on non-carbonate bedrock (8 plots), sites dominated by other broadleaf tree species (15 plots), Picea abies-dominated sites (6 plots) and sites dominated by other coniferous tree species (9 plots). For each group we also identified herb-layer indicator species by implementing a statistical method proposed by Dufrêne and Legendre (1997).
Fourth-corner analysis (Dray and Legendre 2008) was employed to evaluate correlations between environmental variables and life-history traits. This analysis combines three matrices: table R (plots 9 environmental variables), table L (species 9 plots) and table Q (species 9 traits). In the case of a community CSR signature, the associations between the community-weighted mean (CWM; Lavorel et al. 2008) for C-score, S-score and R-score and environmental variables were examined with univariate linear regression models. The coefficient of determination (R 2 ) was used as a measure of goodness of fit.
Prior to CCA and fourth-corner analysis, environmental variables were checked for their distribution and in the event of a significant departure from normality, appropriate transformations were used. Each predictor was relativized to vary between zero and one, meaning that all explanatory variables were on the same scale. Significance was determined at P = \ 0.05 in all analyses, which were conducted in R statistical software version 3.5.   Raunkiaer (1934): chamaephyte, geophyte, hemicryptophyte, phanerophyte and therophyte (Martinčič et al. 2007) Onset

Floristic composition
Collectively, the environmental variables included in the final CCA model explained 15.6% of the variation in plant community composition. The constraining variables in the final model were mean annual temperature (MAT), mean annual precipitation (MAP), tree layer shade-casting ability (SCA) and soil pH (Fig. 2) in descending order of importance. The first and the second CCA axis explained 4.6% and 4.5% of the variation, respectively. MAT was most strongly correlated with the first axis (r = -0.92), whereas MAP (r = 0.76), tree layer SCA (r = 0.49) and soil pH (r = 0.58) exhibited a stronger correlation with the second CCA axis (Fig. 2).

Trait composition
A total of 28 associations (out of 196 possible) between life-history traits (trait states) and environmental variables were detected as significant (p \ 0.05). Thirteen of them were positive and 15 were negative correlations (Fig. 3). Three environmental variables (mean annual precipitation -MAP, heat load index -Hi, woody debris cover) did not affect any of investigated life-history traits. Forbs as a plant growth habit were positively influenced by soil pH, total soil nitrogen and tree layer shade-casting ability (SCA). Graminoids showed a negative association with the proportion of broadleaves in the tree layer and tree layer SCA. Shrubs were positively correlated with the proportion of clay in the soil and moss cover but were negatively impacted by tree layer SCA. Tree layer richness increased the frequency/abundance of tree species in the herb layer (Fig. 3).
Chamaephytes as a plant life form exhibited a negative correlation with soil moisture. Hemicryptophytes were positively affected by soil organic carbon, soil nitrogen and altitude. In contrast, this life form was negatively affected by the proportion of clay in the soil, tree layer species richness and mean annual temperature. Phanerophytes showed the reverse pattern. These were positively associated with the proportion of clay, tree layer species richness and mean annual temperature but exhibited a negative correlation with soil nitrogen and altitude (Fig. 3).
In our dataset, phenological traits were most tightly linked to tree layer characteristics. The onset of flowering was negatively correlated with the proportion of broadleaves in the tree layer and tree layer SCA. The same was true for the end of flowering, which was also negatively correlated with tree species richness. This explanatory variable showed a negative correlation with the duration of flowering (Fig. 3).
The analysis of linear regressions between community-weighted means (CWM) of C-, S-and R-scores and individual environmental variables showed four significant relationships (Fig. 4). The C-score was positively related to tree layer SCA (p \ 0.05, R 2 = 0.08). The community-level S-score decreased with increasing soil pH (p \ 0.05, R 2 = 0.11) and increasing tree layer SCA (p \ 0.01, R 2 = 0.14). The CWM of the R-score in herb-layer communities was negatively linked to mean annual precipitation (p \ 0.01, R 2 = 0.20) (Fig. 4).

Discussion
This study highlighted the environmental determinants of floristic and trait composition in assemblages of vascular plant species in the herb layer of temperate forests across Slovenia. The dominant factors which explained the floristic composition are similar to those reported by Leuschner and Lendzion (2009). Our results showed that floristic composition is driven by climatic parameters (mean annual temperature, mean annual precipitation), soil properties (soil pH) and forest stand characteristics (tree layer SCA) regulating light availability at the forest floor. Light conditions were also an important environmental factor for herblayer composition in Hungarian temperate forests (Márialigeti et al. 2016) and light availability drives biodiversity changes in herb-layer communities across Europe (Depauw et al. 2019). Meaningful overstory influences on the herb and shrub communities of mature forests were confirmed in the studies of  Table 1 McKenzie et al. (2000) and Dölle et al. (2017). In our previous study (Kermavnar et al. 2021) we illustrated the complete absence of a climatic signal for herblayer species diversity (species richness, herb-layer cover and evenness). In contrast, the present results indicate the profound effect of climatic parameters on herb-layer species composition. Different aspects of plant community structure (composition vs. diversity) are thus affected by different environmental factors. Besides a strong influence of abiotic factors, biotic factors like dispersal ability, competition or herbivory can explain forest species distribution patterns on a finer spatial scale (Leuschner and Lendzion 2009).
The novelty of our approach lies in the parallel assessment of floristic and trait-based perspectives. Although taxonomic and functional components are often closely linked to each other, broadening the focus from taxonomy to functional traits is widely accepted (Terwei et al. 2016;Tardella et al. 2019). It is noteworthy, however, that some simple traits (yet easily available for the majority of species, e.g. plant life and growth form) allow fewer general conclusions, as plant ecologists seldomly integrate wholeplant traits in comparison to numerical traits describing the global spectrum of plant form and function (Durak and Durak 2021).
One issue worth discussing is covariation between selected environmental variables, which can greatly confound the interpretation of their real effects on species and trait composition. For example, tree layer composition and diversity are often linked to soil properties and climatic parameters. Forest soil properties, however, also depend on the dominant tree species, macroclimate and local topographic factors (Weigel et al. 2019). A reliable distinction between individual influences of different explanatory factors is thus difficult and (in)direct effects may be masked by myriad possible interactions amongst them. Interdependence between predictors could lead to some bias and accurate predictions might be even more complicated if forest management or legacy of past land use (trajectory of successional development) Fig. 4 Linear regression analysis between communityweighted mean of C-(C-strategists), S-(S-strategists) and R-(R-strategists) scores (%) and significant environmental variables: Tree_SCA = shade-casting ability of the tree layer, soil pH and MAT = mean annual temperature. The grey ribbon around the curve indicates 95% confidence interval plays an important role in structuring herb-layer assemblages (Brunet et al. 2011;Depauw et al. 2019). Forest community composition is driven by environmental variables that are under the direct control of the actual silvicultural management (Tinya et al. 2021). Notwithstanding, based on the correlation matrix (Appendix 1), we assume that the degree of such confounding effects is rather low in our dataset. The highest correlations are found for variables within the same group of predictors than across predictor groups (e.g. soil parameters generally show a higher association with each other compared to the associations with the tree layer or climatic factors).
We selected a set of functional traits capturing the main axes of trait variation in forest herb-layer vegetation. Categorical traits such as plant life and growth forms have been shown to be a good surrogate for explaining the variability in key traits describing fundamental trade-offs between leaf economics, plant size and reproduction (Durak and Durak 2021). In our case, forbs, graminoids and shrubs were correlated with different soil properties and tree layer characteristics. These features of forest stands are impacted by management disturbances, as shown by Ramovs and Roberts (2005). The same environmental factors influenced plant growth habits in different ways. For example, most forest forbs benefited from the higher shade induced by dense tree layer closure, whilst graminoids exhibited the opposite pattern. Dölle et al. (2017) showed higher proportion of graminoids and shrubs in the sessile oak (Quercus petraea) stands compared to beech-dominated (Fagus sylvatica) plots. Improved light conditions support the expansion of plant species adapted to a high-light environment (Kermavnar et al. 2021). Many herb-layer communities in our coniferous plots experienced prominent domination of grasses with pervasive clonal reproduction abilities (some examples: Calamagrostis villosa and Deschampsia flexuosa in Picea abies forest, Sesleria autumnalis and Brachypodium rupestre in Pinus nigra forest, Molinia caerulea subsp. arundinacea in Pinus sylvestris forest; see Table 1 in Kermavnar and Kutnar 2020). The differential influence is even more evident for plant life form. The majority of environmental factors (except for SOC) that were positively correlated with hemicryptophytes induced negative control on phanerophytes and vice versa. Such patterns suggest excluding interactions between these two most abundant life forms in our dataset. Tardella et al. (2019) reported that topographic factors (altitude and slope aspect) induced strong environmental filtering which changed the relative abundance of herb-layer traits and functional richness of life forms in Sub-Mediterranean coppiced woods. In contrast, our results suggest that heat load index (combining both slope aspect and slope steepness) was not significant in predicting any of the dependent variables. This is likely because we analysed a broad spectrum of forest types across Slovenia rather than comparing communities on smaller (regional, local) spatial scale.
The flowering attributes (start, end and duration) of herb-layer species were exclusively associated with tree layer characteristics. The results indicate that herb-layer communities in deciduous forests are characterized by the early onset and short duration of flowering. Vernal geophytes are known for their adaptation of producing flowers before the deciduous overstory canopy develops and consequent light availability at the forest floor is reduced (Whigham 2004), reflecting phenological asynchrony between spring and summer plant species. This group of species also have a very short flowering period compared to other species and can be additionally differentiated in terms of some key functional traits related to leaf economics and plant size spectra (e.g. specific leaf area, leaf dry matter content, plant height, seed mass) (Durak and Durak 2021;Cubino et al. 2021). In addition, phenological traits of understory communities perhaps depend more on species (i.e. are species-specific) and not only on canopy closure or tree species composition.
CSR signature was associated with tree layer SCA, soil pH and mean annual temperature. Tree layer SCA (a proxy for light conditions) was the most important factor. This is in line with the study of Catorci et al. (2011), who found that the variation in trait data in the herb layer of Sub-Mediterranean forests was best explained by photosynthetically active radiation. Looking at the patterns for our plots, the C-score increased along the tree layer SCA gradient, suggesting that competitors are more abundant in shaded understory environments. Given that 40% of our plots were dominated by Fagus sylvatica shade-tolerant herbaceous species typically occurring in beech understory may have contributed to the observed trend. Beech forest species show a rather narrow spectrum of ecological niches (higher requirements for soil nutrients and moisture) and share similar functional traits associated with a competitive and/or stress-tolerant strategy but limited dispersal ability (Willner et al. 2009). Understories of broadleaved deciduous European forests (such as beech forests on non-acid soils or ravine forests) more commonly host species with acquisitive leaf economics (Cubino et al. 2021). Such a pattern can also be interpreted with respect to shade-tolerant competitive tree species (e.g. Fagus sylvatica) in the herb layer, facilitated by overstory seed sources. Herb layer in pure beech stands is often dominated by juvenile trees, mainly beech and competitive herbaceous species (Dölle et al. 2017).
Light availability was also decisive for the abundance of S-strategists. However, in contrast to the common belief that the main abiotic stress in forests is manifested in low light availability, we found that the S-score decreased along the tree layer SCA gradient. One possible explanation for the deviation from the prediction might be because our plots with high tree layer SCA values were characterized by higher soil pH (the correlation can be recognized on Fig. 2), e.g. different beech forests on carbonate bedrock. In addition, environmental constraints often do not exert selection on single traits but instead exert selection on multiple traits simultaneously (Candeias and Fraterrigo 2020). We can conclude that the main environmental stress within the studied dataset was induced by low soil pH, whilst lack of light seems to be less expressed as a constraining factor.
The relative importance of a ruderal strategy decreased along the temperature gradient, suggesting the larger impact of natural and anthropogenic disturbances (e.g. windthrows, pest infestations, drought stress and associated tree mortality, sleet (ice break), salvage logging) in high-altitude compared to lowland forests. Additionally, Picea abies sites were mostly characterized by competitive graminoids (e.g. grasses, such as Calamagrostis villosa or Deschampsia flexuosa), which also exhibit substantial ruderal behaviour as they are able to spread after local canopy disturbance through means of extensive vegetative propagation.
Overall, tree layer SCA was recognized as an environmental factor with the highest number of significant associations with either floristic or trait composition (eight correlations in total). Another important characteristic related to forest stands was tree layer richness, which proved significant in five occasions, followed by mean annual temperature (four correlations). Amongst the other potential predictors, tree species composition, soil pH, soil nitrogen content and soil clay content (all exhibiting three correlations in total) were also strong predictors for herb-layer vegetation composition. Given that the main environmental determinants of herb-layer composition were linked to stand characteristics (we did not include other potential factors, such as vertical structure, which has been shown to be deterministic for plant composition in mature forests; Dölle et al. 2017) and forest soil properties (but not accounting for litterlayer attributes; e.g. Mölder et al. 2008), appropriate forest management interventions should be applied to support the conservation of local site-adapted tree layer species composition, richness and structure without significantly impacting soil layers and consequently the preservation of disturbance-sensitive herbaceous species. Understanding the determinants of herb-layer composition and diversity provides a theoretical basis for sustainable management practices (Wei et al. 2020). Clearly, efforts to preserve biodiversity and promote long-term sustainability in temperate forests cannot ignore understory species (e.g. Š ebesta et al. 2017).
In addition, the climatic signal evidenced in the significant effect of mean annual temperature and precipitation on floristic composition suggests that plant communities in temperate forests are likely sensitive to ongoing climate change, i.e. warming and drying (Amatangelo et al. 2014). Since functional composition largely reflects floristic composition, similar effects of future climate shifts can be expected for community trait patterns (Felsmann et al. 2018). As highlighted by Cubino et al. (2021), trait variation of understory plant communities at the continental scale is driven by macroclimate and soil factors and modulated by overstory canopy cover. All these factors (or their surrogates) were addressed also in our study. Such analyses further warrant the importance of assessing both biodiversity and functional responses to increase our understanding of the possible impact of global change on the herb layer (Depauw et al. 2019). Climate changes act as triggers for natural and anthropogenic forest disturbances. Based on a subset of the plots included in the present study, Kutnar et al. (2019) demonstrated that more disturbed sites showed larger temporal changes in species composition compared to less-disturbed sites, suggesting that forest disturbances caused significant compositional shifts over a 10-year period across Slovenian forest ecosystems. reviewers for their constructive comments on the previous version of this manuscript and Klemen Eler for his valuable suggestions regarding the statistical analyses. We are grateful to Philip Jan Nagel for the linguistic revision of the manuscript.
Funding The research was funded by the Slovenian Ministry of Agriculture, Forestry and Food (JGS tasks 1/2 and 1/3) and by the Slovenian Research Agency (research core funding No. P4-0107). J.K. was funded by the Slovenian Research Agency for his Young Researcher Programme (Contract no. 1000-18-0404).
Availability of data and material Not applicable.
Code availability Not applicable.

Declarations
Conflict of interest The authors of this manuscript have no conflict of interest to declare.
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://creativecommons.org/licenses/by/4.0/.