Sergio de Tomás Marín1*†
Javier Galán Díaz2†
Jesús Rodríguez-Calcerrada3†
Iván Prieto4†
Enrique G. de la Riva1,4†- 1Department of Ecology, Brandenburgische Technische Universität Cottbus-Senftenberg, Cottbus, Germany
- 2Department of Pharmacology, Pharmacognosy and Botany, Faculty of Pharmacy, Universidad Complutense de Madrid, Madrid, Spain
- 3Functioning of Forest Systems in a Changing Environment Research Group, Universidad Politécnica de Madrid, Madrid, Spain
- 4Ecology Department, Faculty of Biology and Environmental Sciences, Universidad de León, León, Spain
Introduction: Functional trait-based approaches are extensively applied to the study of mechanisms governing community assembly along environmental gradients. These approaches have been classically based on studying differences in mean values among species, but there is increasing recognition that alternative metrics of trait distributions should be considered to decipher the mechanisms determining community assembly and species coexistence. Under this framework, the main aim of this study is to unravel the effects of environmental conditions as drivers of plant community assembly in sub-Mediterranean ecotones.
Methods: We set 60 plots in six plant communities of a sub-Mediterranean forest in Central Spain, and measured key above- and belowground functional traits in 411 individuals belonging to 19 species, along with abiotic variables. We calculated community-weighted mean (CWM), skewness (CWS) and kurtosis (CWK) of three plant dimensions, and used maximum likelihood techniques to analyze how variation in these functional community traits was driven by abiotic factors. Additionally, we estimated the relative contribution of intraspecific trait variability and species turnover to variation in CWM.
Results and discussion: The first three axes of variation of the principal component analyses were related to three main plant ecological dimensions: Leaf Economics Spectrum, Root Economics Spectrum and plant hydraulic architecture, respectively. Type of community was the most important factor determining differences in the functional structure among communities, as compared to the role of abiotic variables. We found strong differences among communities in their CWMs in line with their biogeographic origin (Eurosiberian vs Mediterranean), while differences in CWS and CWK indicate different trends in the functional structure among communities and the coexistence of different functional strategies, respectively. Moreover, changes in functional composition were primarily due to intraspecific variability.
Conclusion: We observed a high number of strategies in the forest with the different communities spreading along the acquisitive-conservative axis of resource-use, partly matching their Eurosiberian-Mediterranean nature, respectively. Intraspecific trait variability, rather than species turnover, stood as the most relevant factor when analyzing functional changes and assembly patterns among communities. Altogether, our data support the notion that ecotones are ecosystems where relatively minor environmental shifts may result in changes in plant and functional composition.
1 Introduction
The functional trait-based approach is a broadly-known tool applied to the study of mechanisms governing community assembly along environmental gradients (Garnier and Navas, 2012; de la Riva et al., 2018a; Dong et al., 2020; Kermavnar et al., 2022), being the initial step towards developing a more mechanistic understanding of how the environment influences ecosystem structure and function (Wieczynski et al., 2019). The classical approach for assessing the assembly processes within a given community have mostly aimed at quantifying the trait composition of plant communities, i.e., the presence of species with particular trait values within communities (e.g. Maire et al., 2012; de la Riva et al., 2016a; Hejda et al., 2019) by focusing on the first two moments of a trait distribution weighted by species abundance, i.e. community-weighted mean and variance of their trait distributions (Garnier et al., 2004; Violle et al., 2012). Mean and variance trait values are only meaningful when traits have a normal (gaussian) distribution within the community. However, the distribution of trait values within communities often deviates from the symmetric normal distribution, for example as a result of rapid environmental changes or asymmetric competition favoring a limited portion of the community with specific trait values (Enquist et al., 2015; Le Bagousse-Pinguet et al., 2017; Gross et al., 2021). Thus, there is an increasing recognition that alternative aspects of trait distributions should be considered to decipher the mechanisms determining community assembly and species coexistence, such as the third and fourth community-weighted moments (skewness and kurtosis) (Le Bagousse-Pinguet et al., 2017; Danet et al., 2018; Wieczynski et al., 2019; Gross et al., 2021). The community-weighted skewness and kurtosis bring light on the shape of trait distributions in naturally assembled communities (Figure 1) (Gross et al., 2017; Le Bagousse-Pinguet et al., 2021). The skewness quantifies the asymmetry of the trait distribution, where extreme skewness values are found when rare species within a community have infrequent trait values (Le Bagousse-Pinguet et al., 2021). The kurtosis, in turn, is a measure of the evenness of trait distributions, where extreme negative kurtosis values reflect an even distribution of trait values in the community or even bimodal trait distributions, and extreme positive kurtosis values reflect a reduced range of trait values (Gross et al., 2021; Le Bagousse-Pinguet et al., 2021). Thus, combining the skewness and kurtosis of a trait distribution with the mean value for a given trait, provides a more complete understanding of how environmental changes shape the functional diversity and species assembly of plant communities (Wieczynski et al., 2019; Gross et al., 2021). Few works have combined the study of abundance-weighted trait moments, i.e. community weighted mean, skewness and kurtosis, for a wide range of traits at the whole plant level, together with potential environmental drivers simultaneously (Wieczynski et al., 2019), despite the importance of looking simultaneously at multiple traits, since different plant organs may respond differently to different selection pressures (Carvajal et al., 2019). In fact, to the best of our knowledge, this approach has never been performed in ecotones or with belowground traits, as it is done in the present study.
Figure 1 Schematic representation of the shifts in the shape of a normal trait distribution (A) produced by skewness (B: positive; C: negative) and kurtosis (D: high; E: low).
The variation of the indices describing the position of the dominant trait values in the trait space and shape of the trait distribution do not reflect the source of variation along environmental gradients. That is, changes in the functional structure of plant communities may arise due to the replacement of species (species turnover) and/or intraspecific trait variability (ITV). These two sources of variation may contribute in different degrees to the overall change in the functional structure of plant communities (Albert et al., 2010; de Bello et al., 2011; Pérez-Ramos et al., 2012; de la Riva et al., 2016a). Hence, the estimation of the trait distribution moments at the whole-plant level, i.e. including both above- and belowground traits, and the role of species turnover and ITV in the variation of these distribution moments among communities, may reveal different responses of plants to environmental changes, offering higher insights into the assembly processes taking place in communities along environmental gradients (Lepš et al., 2011; Kichenin et al., 2013; de la Riva et al., 2016a).
Ecotones are transitional areas between adjacent biocenosis with high levels of biodiversity, as they host species from the adjacent biocenosis as well as unique species endemic to the ecotone (Hufkens et al., 2009; Vila-Viçosa et al., 2020b; Shea et al., 2022). Ecotones are of especial interest for conservation (Carlson et al., 2014; Rawal et al., 2018; Lopes et al., 2020) for their ecological uniqueness (Danso Marfo et al., 2019), and for contributing to stabilize adjacent ecosystems (Pogue and Schnell, 2001). Further, community assembly processes in ecotones are very dynamic, as slight modifications in environmental conditions lead to substantial changes in plant composition (Ruiz-Labourdette et al., 2012; Brice et al., 2020), which makes them ideal study systems to test such processes (Rahman et al., 2021). However, understanding which environmental factors govern the process of species assembly in transitional areas remains challenging (Kitagawa et al., 2020; Rahman et al., 2021) and could be key to understanding how these communities will respond to future climatic conditions.
Sub-Mediterranean areas are transition areas between Mediterranean and temperate climates. The Iberian Peninsula, in the south of Europe, represents one of the largest ecotones of this kind (Sánchez de Dios et al., 2009), being the low-latitudinal edge for the distribution of many central European species (Gil et al., 2010; Dorado-Liñán et al., 2017; Hernández-Lambraño et al., 2021). Sub-Mediterranean regions host Mediterranean and Eurosiberian species as well as many endemic species (Vila-Viçosa et al., 2020a; Box, 2021; Cantoral et al., 2023), creating a mosaic of unique plant communities with contrasting assembly processes (de la Riva et al., 2023). Sub-Mediterranean regions in the Iberian Peninsula are very sensitive to slight environmental alterations (Sánchez de Dios et al., 2009; Vila-Viçosa et al., 2020a), which makes them highly vulnerable to ongoing climatic change and thus critical areas for conservation (Benito Garzón et al., 2008; Cantoral et al., 2023). Understanding trait-based assembly rules in sub-Mediterranean regions can be critical for understanding how species coexist in these and other ecotones, which could be important to predict the trajectory of ecotone plant communities facing increasingly warmer and drier conditions (Benito Garzón et al., 2008; Frenette-Dussault et al., 2013; Pereira et al., 2021).
The broad aim of this study is to unravel the effects of environmental conditions as drivers of the assembly of plant communities in a sub-Mediterranean forest in Central Spain. To this end, we measured four leaf, two stem and three root traits in a total of 411 individuals of 19 woody plant species in six plant communities. We then calculated three trait distribution moments weighted by the relative abundance of each species in each community, i.e., the community-weighted mean (CWM), skewness (CWS) and kurtosis (CWK), of major axes of variation at the aboveground (i.e. Leaf Economics Spectrum or LES, Wright et al., 2004) and belowground levels (i.e. Root Economics Spectrum or RES; Prieto et al., 2015) across plant communities, and assessed their relationship with environmental variables to estimate which, and to what extent, environmental drivers rule the assembly of sub-Mediterranean plant communities. Additionally, we estimated the proportion of variation in CWMs attributable to species turnover and to ITV to have a detailed view of assembly processes occurring in these communities. Using this novel approach, we will gain a great understanding on the functional structure of sub-Mediterranean plant communities and, ultimately, will be able to guide conservation efforts to maintain their proper functioning (Walker, 1995).
2 Materials and methods
2.1 Study area
The study was conducted in the ‘Hayedo de Montejo’ (Montejo hereafter), a 125 ha. forest located in the ‘Sistema Central’ range in Spain (41°7’N, 3°30’W), between 1250 and 1550 m a. s. l. Orientation is predominantly northeast, although it ranges from south to north in some locations. Slope varies from nearly flat at higher elevations to ca. 58% in the steepest areas. Soil depth ranges from 30 cm in upper areas to 150 cm down in the valley. The dominant bedrock is micaceous gneiss and soils on top are acidic with a sandy-loam texture (López Santalla et al., 2003; Gil et al., 2010). The climate is sub-Mediterranean, with Mediterranean influence in the east and a continental climate in the west. Mean annual precipitation was 858.8 mm and mean annual temperature was 9.7°C for the period 1994-2021 (Montejo meteorological station). Summers are hot with a marked dry period between July and August, typical of Mediterranean climate.
The Hayedo de Montejo is a sub-Mediterranean, mixed beech-oak forest (Fagus sylvatica-Quercus pyrenaica), with a high diversity of plant species and communities coexisting in a relatively small area. In this forest, Eurosiberian species coexist with typically Mediterranean species in the overstory and the understory (Sánchez de Dios et al., 2009; Gil et al., 2010; Rivas-Martínez et al., 2011), which harbors some herbaceous species of special interest for conservation for their rarity, such as Paris quadrifolia L., Narcissus pseudonarcissus ssp. pseudonarcissus L. or Aconitum napellus L. (Hernández Bermejo et al., 1982; Moreno et al., 2005; Gil et al., 2010). The forest is dominated by the temperate tree species European beech (Fagus sylvatica L.; hereafter beech) and sessile oak (Quercus petraea (Matt.) Liebl.), both at the southernmost limit of their distribution range, and the Mediterranean Pyrenean oak (Quercus pyrenaica Willd.), at the core of its distribution range (Rubio-Cuadrado et al., 2020; Rubio-Cuadrado et al., 2021; de Tomás Marín et al., 2023).
This forest has been traditionally managed as a wooded pasture (Pardo and Gil, 2005). After centuries of forest exploitation, mainly for firewood and cattle grazing, these traditional land uses ceased in the early 1960s (López Santalla et al., 2003; Pardo and Gil, 2005) – although, today, it is still relatively easy to find a few cows grazing in the forest. These land-use changes transformed the former open woodland structure with dispersed trees and scarce tree recruitment into the dense secondary forest that we find today (Pardo and Gil, 2005; Rubio-Cuadrado et al., 2021). The Hayedo de Montejo forest is included in the ‘Ancient and Primeval Beech Forests of the Carpathians and Other Regions of Europe’ list of European beech forests included in the UNESCO World Heritage Site list. Furthermore, it is embedded in the ‘Sierra del Rincón’ (Rincón mountain range), designated as a Reserve of the Biosphere by UNESCO in 2005.
2.2 Sampling design
Three forest inventories (FI) were carried out in Montejo in the years 1994, 2005 and 2015, respectively. In each FI, 125 circular plots of 30-m diameter were systematically established in a 100 x 100 m square grid covering the entire forest. Based on the last FI, we divided the forest into six different areas (communities hereafter) according to the distribution of the three most abundant tree species, i.e., beech, Pyrenean oak and sessile oak (see also de la Riva et al., 2023). An additional community dominated by shrubs where trees are scarce and dispersed was included. These communities were classified based on the abundance of the three main species recorded in the FI of 2015 as follows: (I) Fagus sylvatica, (II) Quercus pyrenaica or (III) Quercus petraea, when >60% of the abundance in terms of tree density corresponded to one of these three species, respectively; (IV) mixed forest 1 (transition areas between F. sylvatica and Q. pyrenaica), (V) mixed forest 2 (transition area with co-dominance of F. sylvatica, Q. pyrenaica and Q. petraea) and (VI) a shrubland community (Supplementary Material, Appendix S1, Figure S1).
Ten random plots from the FIs within each plant community (I to VI) were selected for sampling. In each plot, four parallel 20-m-long transects perpendicular to the slope were randomly set, two on each side of the plot, with a minimum distance of 2 m between transects. Species composition (i.e. identity) and abundance (i.e. % cover) were recorded for each woody species intercepted by the transect line (abundance was measured as the length in meters of canopy of each woody species projecting on the transect line; total abundance may thus exceed 20 m due to overlaying vegetation layers). For leaf, stem and root trait measurements, we selected five individuals per species and community, except for Huber values, where n=3 (see below for description of this trait). Traits were measured in species that together made up at least 90% of the maximum cumulative cover within a community. With this criterion, a total of 19 species were selected across communities: Adenocarpus hispanicus, Crataegus monogyna, Cytisus purgans, Cytisus scoparius, Erica arborea, Erica australis, Fagus sylvatica, Genista florida, Hedera helix, Ilex aquifolium, Juniperus communis, Lavandula stoechas, Prunus avium, Quercus petraea, Quercus pyrenaica, Rosa sp., Rubus sp., Sorbus aria and Sorbus aucuparia. All species occurred at least in two plant communities, if not more, except Cytisus purgans, that was found only in the shrubland community (Supplementary Material, Appendix S2, Table S2).
2.3 Plant sampling and trait measurements
Plant sampling was carried out in June and July 2021, during the peak of the growing season. From this material, five aboveground and three belowground functional traits related to water and nutrient acquisition were measured (see functional roles in Table 1): Leaf Dry Matter Content (LDMC; leaf dry mass per unit of water-saturated fresh mass; mg g-1), Specific Leaf Area (SLA; leaf area per unit of leaf dry mass; m2 kg-1), leaf carbon isotopic composition (δ13C; ‰), leaf carbon to nitrogen ratio (leaf C:N), Stem Dry Matter Content (SDMC; stem dry mass per unit of water-saturated fresh mass; mg g-1), Huber value (Hv; sapwood cross-sectional branch area to supported leaf area ratio; cm2 cm-2), Specific Root Area (SRA; root area per unit of root dry mass; m2 kg-1), average root diameter (Rdi; mm), and root C:N ratio (root C:N). Leaf δ13C is a proxy of intrinsic water use efficiency (iWUE, Prieto et al., 2018; Ma et al., 2023). The Huber value reflects the mass investment in xylem area with respect to leaf area (Mencuccini et al., 2019); it is thus a measure of how species adjust their leaf and sapwood areas to water supply (Carter and White, 2009; Zhang et al., 2019) and is related to groundwater depth extraction (Carter and White, 2009). Thus, we used this trait as a proxy of the hydraulic architecture of the plants.
Table 1 List of the nine functional traits considered in this study, their abbreviations, units and their functional roles.
For leaf and stem measurements, two-year-old branches with fully expanded, sun-exposed mature leaves from five healthy adult individuals per species and community were sampled. A variable number of leaves were measured depending on the average leaf size of each species. For root measurements, fine roots (<2 mm diameter) from five individuals per species and community were collected within the first 20-30 cm of soil by excavating next to the plant stem base to ensure that roots belonged to the selected individual. Roots were stored fresh in a cooler after collection and taken to the lab on the same day, where they were rinsed with distilled water to eliminate adhered soil particles, and then frozen until they were measured. The Hv was measured at regional level, i.e. selecting a total of three individuals per species for all the species, except for F. sylvatica and Q. pyrenaica, for which we measured 5 individuals per community (de Tomás Marín et al., 2023).
Leaf and root material was scanned with an EPSON® V850 PRO scanner at a resolution of 600 dpi. Leaf area was measured from the scanned images using the ImageJ software (Schneider et al., 2012). Root area and diameter were measured from the scanned images using the WinRHIZO 2009 software (Regent Instruments Inc., Quebec, Canada). After scanning, roots were oven-dried at 60°C for 48h prior to measurements. Leaf carbon isotopic composition and C and N concentrations were measured simultaneously with a THERMO/Finnigan MAT V isotope ratio mass spectrometer, coupled to a THERMO Flash EA 1112 elemental analyzer via a THERMO/Finnigan Conflo IV- interface. Leaf carbon isotopic composition is expressed in the conventional delta notation (δ13C) relative to VPDB (Vienna PeeDee Belemnite standard). All plant material was collected, stored and processed following the protocols detailed by Pérez-Harguindeguy et al. (2013). For a detailed protocol of sample harvesting and trait measurements see also de Tomás Marín et al. (2023), de la Riva et al. (2016a) and de la Riva et al. (2016b).
2.4 Abiotic variables
In each selected plot, soil organic matter content (SOM), total soil N, soil nitrate and soil phosphate were determined from the top 15 cm of soil. Soil samples were collected with a hand trowel in three points in the plot, one of the samples corresponding to the center of the plot and the other two approximately 5 m above and below the center, and mixed in a plastic bag to homogenize the soil. Then soils were taken to the lab where they were dried at 50°C to constant weight (~72 h) and sieved at 2 mm before analyses. Additionally, slope, orientation, altitude and soil depth were also measured in the same plots (Supplementary Material, Appendix S3, Table S3). For the description of the methods employed to measure total soil N see Rodríguez-Calcerrada et al. (2019); for the methods employed to measure the rest of abiotic variables see de Tomás Marín et al. (2023).
2.5 Data analyses
As there is a great variety of ways in which plants combine a set of traits to achieve a successful performance (Díaz et al., 2016; Iozia et al., 2023), we focused on trait syndromes rather than on individual functional traits. To this end, we conducted a Principal Component Analysis (PCA) including the nine functional traits considered in this study, measured in 411 individuals belonging to the 19 sampled species, to describe trait trade-offs or trait syndromes (Albert et al., 2011). The first three principal components (PCs) had eigenvalues greater than one and represented three key dimensions of plant ecological strategies (Supplementary Material, Appendix S4, Table S4). That is, PC1, PC2 and PC3 reflected clear spectrums related to the Leaf Economics Spectrum (PC1, LES, Wright et al., 2004), Root Economics Spectrum (PC2, RES, de la Riva et al., 2018b) and hydraulic architecture (PC3, HyArq, Tyree and Ewers, 1991), respectively. We then used the scores of each individual in these three PCA axes to estimate the community-weighted moments and in further analyses.
2.5.1 Community weighted moments along the environmental gradient
We calculated three informative moments related to the dominant trait values (mean) and the shape of the trait distribution (skewness and kurtosis) of the individual PCA scores for each of the three PC axes (LES, RES and HydArq) kept for analyses (see mathematical equations in Supplementary Material, Appendix S5). To scale up from species to community level, these moments were weighted by the relative abundance of each species in each community, resulting in the ‘community-weighted moments’ of the trait distributions. The community-weighted mean (CWM, Garnier et al., 2004) gives a representation of the dominant trait values within a community. This metric is directly related to Grime’s mass-ratio hypothesis (Grime, 1998) and considers that the traits of the most abundant species in the community have the largest influence in ecosystem processes (Bílá et al., 2014). The community-weighted skewness (CWS) quantifies the degree of asymmetry of the trait distribution within the community (Wieczynski et al., 2019; Gross et al., 2021). In our case, the CWS quantifies the degree of asymmetry of the individual values in the three spectrum axes (PC1, PC2 and PC3 scores) and reflects the presence of subordinate species within a community with trait values close to one of the extremes of the PC scores range (Le Bagousse-Pinguet et al., 2017). The community-weighted kurtosis (CWK) refers to the evenness of the distribution of individual values in the three spectrum axes (PC1, PC2 and PC3 scores) (Wieczynski et al., 2019; Gross et al., 2021), where a peaked distribution indicates low levels of diversity along the spectrums and a flat distribution indicates high levels of diversity.
To analyze the effect of the abiotic environment on the variation of the three community-weighted moments (CWM, CWS and CWK) of the different plant ecological dimensions in our communities (LES, RES and HyArq), we built multiple linear regressions models. The continuous abiotic variables (independent variables) and the discrete variable ‘community’ (factor) were introduced in the models as predictors and each of the three distribution moments of the selected PCs, separately, as response variables. The lm function from the base R package stats was used to perform the linear regression models analyses (R Core Team, 2021). To identify the best predictors, we applied the Akaike information criterion corrected for small sample size (AICc) on the full factorial model, i.e., including all the abiotic variables considered in the study (Table 2) and the communities. The selection of the best model (ΔAIC < 2) was performed using the function dredge from the MuMIn R package (Barton and Barton, 2015). The strength of the relationship between the predictors and the response variable for the selected models (lowest AICc) for each community-weighted moment and PC axis, was given by the adjusted R2, obtained with the summary function, also from the stats package (R Core Team, 2021). To assess the % of variance explained by each predictor in each model, we obtained the sum of squares of each predictor with the Anova function from the R package car (Fox and Weisberg, 2019) and calculated the proportion that the sum of squares of each variable was relative to the total sum of squares. To determine differences among the six communities for each community-weighted moment in each of PCs spectrum (LES, RES and HydArq), we ran one-way ANOVAs followed by post-hoc multiple pairwise comparisons (Tukey’s test), to detect differences between pairs of communities.
Table 2 Best-fitted linear regression models. Models are presented for each plant ecological strategy (LES, RES and HydArq) and for each of the Community Weighted moments (CWM, CWS, CWK).
2.5.2 Relative importance of intraspecific trait variability and species turnover
We determined the relative contribution of species turnover, i.e., the sum of the contributions of species occurrence and species abundance, and ITV to the overall change in the functional structure of our communities following the method developed by Lepš et al. (2011), using the function traitflex.anova implemented in the R software. For a detailed explanation of the method and mathematical calculations, please refer to the Appendix S6 in the Supplementary Materials and to Lepš et al. (2011). Variance partitioning was estimated on the CWMs calculated for the first two PC axes (PC1 and PC2), which represent the LES and RES, respectively. We estimated two partitions in parallel, one with “community” as factor and the other with those abiotic variables selected by AICc in the multiple regression linear models (section 2.5.1), to assess the effects of the type of community and of the environment on the variation of CWMs, respectively. We could not carry out the partition for PC3 (HyArq) because Huber values were measured at the regional level, not at the community level.
All analyses were performed in R (version 4.1.2; R Core Team, 2021) interfaced with R Studio (RStudio Team, 2021).
3 Results
3.1 Trade-offs among traits
According to the loadings on the PCA axes (Supplementary Material, Appendix S4, Table S4), the first three axes of variation (PC1, PC2 and PC3) were related to trait trade-offs within plant organs and thus to plant ecological dimensions (Figure 2). The first PCA axis accounted for 28.1% of the overall variation and was representative of the aboveground resource acquisition strategy or Leaf Economics Spectrum (LES, Wright et al., 2004). Positive values of the PC1 represented species with high SLA, indicative of fast resource uptake strategies at aboveground level, while in the opposite end of the spectrum were species with trait values indicative of the resource conservation strategy (high LDMC, leaf C:N, leaf δ13C and SDMC). The second PC axis explained 18.4% of the variability and was related to the Root Economics Spectrum (RES) (Prieto et al., 2015; de la Riva et al., 2018a). Within this axis, positive scores were associated to species with high SRA and negative scores to species with high Rdi. The third PC axis explained 13.7% of the variance and was associated to the Huber value, and thus to the plant’s hydraulic architecture (Mencuccini et al., 2019). Within this axis, positive scores were associated with species with high Huber values, i.e. conservative species.
Figure 2 Principal Component Analysis (PCA). Results show the first two PC axes for the nine functional traits considered in this study measured on 411 individuals from 19 woody species (trait abbreviations are in Table 1; species list is in Supplementary Material, Appendix S2, Table S2). Symbols correspond to the 411 individuals. The type of community which observations belong to is shown with different symbols and colors.
3.2 Relationships among environmental variables and functional structure (community-weighted moments) and differences in the functional structure among communities
Overall, the models showed a high predictive power (Table 2) with all adjusted R2 values over 0.45. Our results showed that the type of community was the most important factor determining differences in the functional structure, with environmental factors having a secondary role (Table 2; Appendix S7). The effect of the different abiotic variables depended on the plant ecological dimension analyzed; altitude and orientation were relevant factors determining variation in the LES (aboveground level), while soil depth and nutrients were more important factors belowground, for RES, and HyArq (only soil depth in this case, Table 2).
Regarding variations in the functional structure among communities, we found clear differences in their community weighted moments related to the main axes of variation (Figure 3). The CWMs showed clear trends in relation to the position of the communities along the acquisitive-conservative axis of resource uptake. Mediterranean communities, i.e., Shrubland and Q. pyrenaica communities, had the lowest CWM scores in the LES (PC1 axis), i.e., these communities were composed by species with high tissue dry matter content and low SLA, and did not differ from each other, indicating a conservative resource-use strategy in both communities (Figure 3). In contrast, the Eurosiberian communities (dominated by either Q. petraea or F. sylvatica) showed the highest CWM scores in the LES axis (Figure 3); the two mixed communities showing intermediate CWM scores in the LES (Figure 3). In relation to the RES (Figure 3), the Shrubland community showed significantly higher CWM scores than the rest of the plant communities, indicating a more acquisitive resource use strategy belowground. F. sylvatica community had the second highest CWM scores, while Q. petraea community showed the lowest scores in the RES. Again, mixed communities (Mixed 1 and 2) showed intermediate CWM scores along the RES axis (Figure 3). With regard to the hydraulic strategy, we found an overall decreasing trend in CWM scores from Mediterranean to Eurosiberian communities (with the exception of Q. petraea community, which had similar CWM scores than Mediterranean communities; Figure 3).
Figure 3 Boxplots of the community-weighted moments (CWM: community weighted mean, CWS: community weighted skewness and CWK: community weighted kurtosis) of each plant ecological strategy (LES, RES and HyArq, related to the principal components 1, 2 and 3, respectively; see materials and methods section for a detailed description of the different plant ecological strategies) for the six communities studied (see section ‘2.2 Sampling design’ and Supplementary Material, Appendix S2 for the composition of each community). LES, Leaf Economics Spectrum; RES, Root Economics Spectrum; HyArq, Hydraulic Architecture. The line inside the box represents the mean value, the box limits represent the SE and the whiskers represent the SD. Different letters indicate significant differences (p < 0.05) among communities. Community-weighted moments are significantly different from zero when its 95%-confidence interval does not overlap with zero (dashed lines).
We observed strong differences among communities on the community-weighted skewness (CWS) scores on the three plant ecological dimensions (LES, RES and HyArq), although we did not detect a clear pattern regarding the biogeographic origin of the communities (Mediterranean vs Eurosiberian). We found a high asymmetry in the scores’ distribution (high positive CWSs) in the two Mixed communities and in the Q. petraea community for the LES, in Shrubland community for the RES and in F. sylvatica community for HyArq. By contrast, F. sylvatica and the Mixed 2 (F.sy-Q.py-Q.pe) communities showed a high negative asymmetry (very low values) for the RES.
Similarly, community-weighted kurtosis (CWK) varied among communities depending on the plant ecological dimension considered. CWK values greater than 0 were found in the Mixed 2 (F.sy-Qpy-Qpe) community for the LES, in the Q. pyrenaica community for the RES, and in the Mixed 1 (F.sy-Qpy) community for the HyArq. By contrast, CWK scores for the Shrubland, Q. petraea and F. sylvatica communities did not differ from 0 for any plant ecological strategy (except for the HyArq axis in F. sylvatica community).
3.3 Contribution of ITV and species turnover to variation in the CWMs for the LES and RES ecological strategies
The percentage of variance explained by both ITV and species turnover was higher for the LES than for the RES (Figure 4; Supplementary Material, Appendix S8, Table S8.1). For both LES and RES, results showed a stronger influence of ITV than of species turnover in the CWMs’ variation among communities. This was especially evident for the RES, for which the effect of species turnover was particularly small compared to that of ITV: close to 61.3% of the variation was explained by ITV against 9.5% explained by turnover (Figure 4; Supplementary Material, Appendix S8, Table S8.1).
Figure 4 Decomposition of total variability in individual trade-offs among traits. LES corresponds to the trade-offs among the study traits related to the Leaf Economics Spectrum. The same applies for the RES in relation to the study traits related to the Root Economics Spectrum.
4 Discussion
4.1 Trait distribution and environmental effects on community assembly
We found large variation in the functional trait composition considered in this study. Our results showed that functional traits at community level (functional parameters) were mainly driven by differences among plant communities, while abiotic factors exerted a minor role on community weighted trait distributions. In this sub-Mediterranean ecosystem, we observed in a recent study (de la Riva et al., 2023) that environmental suitability explained species segregation better than their functional composition, being the six woody plant communities segregated spatially and functionally partly in relation to microenvironmental conditions, mostly soil depth at aboveground level, and aspect and soil nutrients at belowground level.
Overall, the CWMs showed a trend that was in line with the acquisitive-conservative axis of resource uptake (Reich, 2014; Díaz et al., 2016; de la Riva et al., 2016b) at aboveground level. Plant communities dominated by Mediterranean species (Q. pyrenaica, Shrubland and Mixed 1 and 2) are dominated by species with low SLA and SRA, and high LDMC, thick root diameters (Rdi), and high leaf and root C:N. These communities are found in harsh areas which are characterized by low soil depths with stony soils, steep slopes or high insolation (de la Riva et al., 2023; de Tomás Marín et al., 2023), conditions that probably limit water availability. In fact, this water shortage promoted the prevalence of communities with a conservative hydraulic architecture (i.e. water efficient with high leaf δ13C and Huber values), commonly associated with conservative water use strategies at the expense of a slower plant growth (Prieto et al., 2018; Mencuccini et al., 2019). On the other hand, communities dominated by Eurosiberian species (Q. petraea and F. sylvatica communities) are found in areas of Montejo with deeper soils —which likely increase soil water storage— and potentially lower evapotranspiration due to lower insolation. These communities had CWM values associated with faster growth and acquisitive resource-use strategies (i.e., higher SLA and SRA, and lower tissue dry matter content, Rdi, C:N ratios and Huber values).
When belowground strategies were considered, however, Shrubland community was the most acquisitive (high CWM) but also showed the highest positive skewness, indicating the presence of subordinate species with conservative strategies. This community occupies flat areas with shallow soils (and thus lower water availability) and high soil nutrient contents in Montejo (de la Riva et al., 2023). In these areas, greater aridity would promote the convergence of traits towards conservative strategies (Cornwell and Ackerly, 2009; de la Riva et al., 2018b; Carrascosa et al., 2023). Although shallow soils prevent the accumulation of big amounts of water, the flat terrain favors water accumulation in the soil surface during short periods. Under these conditions, with a high-enough nutrient concentration, fast resource-uptake belowground strategies would give plants a competitive advantage to absorb the water and nutrients accumulated in shallow soils (Ryel et al., 2008; Schenk, 2008; Querejeta et al., 2021), and to grow during short wet pulses (Querejeta et al., 2018; Carvajal et al., 2019; de la Riva et al., 2021b). Shrub species from Mediterranean arid environments are able to produce cheaper roots in terms of carbon investment to structural components, resulting in roots with high SRA and high absorptive efficiency (de la Riva et al., 2021b). Further, adaptation to low soil water availability during the summer months also explains the negative kurtosis of Shrubland for HyArq, indicative of an evenness in the distribution of trait values related to the plant water balance, promoting the coexistence of functionally contrasting species (Gross et al., 2021). Overall, these results suggest a niche segregation and/or complementarity in water-use strategies that could minimize competition and favor coexistence (Bello et al., 2019; Illuminati et al., 2022), which could lead to a more efficient use of water in the Shrubland community. These results are in line with a recent study (Illuminati et al., 2022), in which the authors found a clear niche segregation and complementarity in functional traits related to water-use strategy of species coexisting in a Mediterranean shrubland community, and suggested that the presence of species with contrasting water-use strategies may represent a key driver in the assembly of dryland communities.
At the other end of the environmental gradient (see also de la Riva et al., 2023), we found the F. sylvatica community with predominantly acquisitive strategies at the whole plant level. In fact, species with a profligate water-use strategy, such as F. sylvatica (de Tomás Marín et al., 2023), dominated the lowest areas of Montejo where deeper soils are found (de la Riva et al., 2023), which could explain the shift in the hydraulic dimension along the soil depth gradient better than soil nutrients. This profligate use of water and nutrients of F. sylvatica makes it a strong fast-growing species that promotes strong competitive exclusion (Meyer et al., 2003; Rodríguez-Calcerrada et al., 2011; Castaño-Santamaria et al., 2021). In this regard, the highest differences in the skewness values of the RES and HyArq dimensions imply that the asymmetric competition promoted by F. sylvatica selects for species with very different functional profiles to that of the dominant species (Le Bagousse-Pinguet et al., 2017). Thus, sub-dominant species seem to be able to buffer competition with F. sylvatica by displaying contrasting morphological resource uptake strategies, and possessing different hydrological architecture. As a matter of fact, deeper soils may favor the ecohydrological niche segregation among coexisting species (Silvertown et al., 2015), allowing species to buffer competition by the uptake of water from different soil depths.
In relation to the Mixed communities, overall and looking at the CWM, they occupied an intermediate position along the Economics Spectrum (LES, RES and HydArq) and along the environmental gradient in Montejo. These communities, however, exhibited an overall positive asymmetric distribution (CWS = 0.62 for Mixed 1 and CWS = 0.27 for Mixed 2), i.e. a higher frequency of species with traits associated to conservative strategies (CWS > 0), and especially so for the LES dimension. These mixed communities assembled following a competitive hierarchy (de la Riva et al., 2023) where species with acquisitive strategies aboveground are excluded in favor of those with more conservative traits that confer higher competitive advantage in harsher environments (Gross et al., 2013; Carmona et al., 2019). This is further supported by the high kurtosis in both mixed communities for the LES and HyArq, respectively, reflecting a selection of conservative trait values under strong competition (Le Bagousse-Pinguet et al., 2017).
Further, we found that the PC scores for the three plant ecological dimensions were normally distributed across species in the Q. pyrenaica community, with a predominance of conservative traits. On the other hand, Q. petraea community showed the most conservative traits at root level, opposite to the aboveground level, and deviated from normality for the LES (CWS significantly different from 0), which suggests contrasting resource uptake strategies at leaf level in the less dominant species with regard to the dominant aboveground acquisitive profile of the community. Thus, and similarly to the mixed communities, the departure from normality in the LES dimension in the Q. petraea community (indicated by its positive CWS values) suggests a process of hierarchical competition (de la Riva et al., 2023) in which species with aboveground conservative strategies are favored. Interestingly, we observed a mismatch between CWM root and leaf strategies. This mismatch can be explained by the ability of Q. petraea to take up water from deep soil layers to cope with drought stress during dry periods (Zapater et al., 2011; Bello et al., 2019; Staszel et al., 2022). This habit would be consistent with a vertical decoupling between water and nutrient uptake (Querejeta et al., 2021), explaining its conservative fine root traits in the topsoil (<20 cm) associated with a lower capacity for nutrient uptake.
To the best of our knowledge, our study is the first one that incorporates the use of community-weighted moments to disentangle the mechanisms and processes driving community assembly in a sub-Mediterranean ecotone. In this regard, the use of community-weighted moments of a trait distribution allowed us to get a deeper insight into the mechanisms driving community assembly in this sub-Mediterranean forest that would be otherwise overlooked, supporting the need to use metrics-based indices beyond the mean for a complete understanding of the underlying assembly patterns of plant communities (Enquist et al., 2015; Le Bagousse-Pinguet et al., 2017; Wieczynski et al., 2019; Guerin et al., 2022).
4.2 Disentangling the role of species turnover and intraspecific variability in shaping the trait distribution between communities
Our results from the analysis of variance decomposition showed that variations in community functional traits along Montejo were mainly explained by changes in intraspecific trait variability (ITV above 60%) rather than by species turnover, in line with previous findings (Jung et al., 2014; Siefert and Ritchie, 2016; Niu et al., 2020; Guo et al., 2022; or Spitzer et al., 2023) but in sharp contrast with others (de Bello et al., 2011; Pérez-Ramos et al., 2012; Kichenin et al., 2013; de la Riva et al., 2016a; Weemstra et al., 2021). The contrasting results found in the literature suggest that the relative contribution of interspecific vs intraspecific variability to the functional structure among plant communities vary with the spatial scale considered and the range of environmental variation. For instance, in this study, 17 out of the 19 species occurred in at least two communities if not more, which could explain the low contribution of species turnover. These results add further support to the hypothesis that ITV gains importance at smaller spatial scales (Lajoie and Vellend, 2015; Siefert et al., 2015; Pichon et al., 2022), when environmental conditions are less variable and interspecific trait variation is low (Cordlandwehr et al., 2013; Petruzzellis et al., 2017). Our results add up to the increasing evidence of the preponderant role that above- and belowground ITV plays on community assembly, suggesting that, in transitional areas such as ecotones, species are not solely filtered based on their mean trait values but also by their trait variability (Siefert, 2012). In fact, ITV is higher across the whole forest than at the community level in the Hayedo de Montejo (de la Riva et al., 2023), supporting the notion that adaptive shifts in traits within species allow them to establish in different communities along transitional areas.
Despite a majority of studies focusing on aboveground traits (e.g. Lepš et al., 2011; Pérez-Ramos et al., 2012; Siefert, 2012; Kichenin et al., 2013; Jung et al., 2014; Siefert et al., 2015; Siefert and Ritchie, 2016; Henn et al., 2018; Niu et al., 2020; Chelli et al., 2021; Guo et al., 2022; Pichon et al., 2022), evidence showing the contribution of belowground intraspecific trait variability to community assembly is starting to accumulate (de la Riva et al., 2016a; Navarro-Fernández et al., 2016; Weemstra et al., 2021; Hogan et al., 2023; Spitzer et al., 2023). In this study we observed that the extent and sources of variation of root traits followed a similar pattern to those of leaves, even with a relatively higher importance of ITV for roots than for leaves (Figure 4; Supplementary Material, Appendix S8, Table S8.1). A recent study showed that ITV explained the mean trait variation of morphological root traits rather than of chemical root traits (Spitzer et al., 2023), while previous studies carried out in Mediterranean woody communities showed higher ITV for the number of arbuscular mycorrhizas vesicles than for SRL (Navarro-Fernández et al., 2016). This indicates that ITV differs between root traits, as found at the aboveground level (Lepš et al., 2011; Kichenin et al., 2013). Considering root trait syndromes at unified dimensions help to depict the role of ITV vs species turnover in community assembly, because traits co-vary from trade-offs between different plant functions, producing syndromes that influence species’ fitness (Albert et al., 2011).
5 Conclusions
In this study, we discerned the assembly patterns of woody plant communities in a sub-Mediterranean forest, using a whole-plant, trait-based approach. We observed a high number of strategies in the forest, reflecting coordination in the change of above- and belowground traits, with the different communities spreading along the acquisitive-conservative axis of resource-use, partly matching their Eurosiberian vs Mediterranean nature, respectively. Additionally, we showed that ITV is a highly relevant factor that should be incorporated when analyzing functional changes and assembly patterns in plant communities, both above- and belowground, especially at small spatial scales. We also found slight relationships between changes in trait syndromes and topographic and soil factors. This is probably due to the small spatial scale and sampling area of the study and to the potential biotic interactions occurring among the different plant species, which may prevent us from detecting stronger functional responses to environmental conditions (See Figure 5 for a conceptual illustration of the main findings of the study).
Figure 5 Schematic representation of the main findings of this study. Overall, the communities in Montejo spread along the conservative-acquisitive axis of resource use, partly matching their biogeographic origin (Mediterranean vs Eurosiberian). The positive skewness indicates a general shift of the Montejo forest towards more conservative trait values. The low kurtosis, shown by a flattened trait distribution curve, suggests the coexistence of different plant ecological strategies. Furhtermore, the relative contribution of ITV to the changes in the functional structure of Montejo is much higher than that of turnover, both at above- and belowground levels. Additionally, the type of community is the most relevant factor determining differences in the functinoal structure among communities, with environmental variables having a minor role. Shr, Shrubland community; Qpy, Q. pyrenaica-dominated community; Mixed 1, community co-dominated by Fagus sylvatica and Q. pyrenaica; Mixed 2, community co-dominated by F. sylvatica, Q. pyrenaica and Q. petraea; Qpe: Q. petraea-dominated community; Fsy, F. sylvatica-dominated community; ITV, Intraspecific Trait Variability.
Altogether, these findings represent an important step towards a more general understanding of the processes driving the assembly and functional composition of sub-Mediterranean plant communities. Our data support the notion that ecotones are fragile regions where minor environmental shifts may result in abrupt changes in ecosystem structure and composition. Thus, understanding which and how abiotic factors produce changes in the composition and functional structure of these transitional regions (e.g. species turnover vs intraspecific variability) may help to their protection in the face of global change. Our results set the ground for future research of community assembly processes and mechanisms of plant communities in these and other ecotone regions.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
ST: Data curation, Formal Analysis, Writing – original draft, Writing – review & editing. JG: Data curation, Formal Analysis, Writing – review & editing. JR-C: Writing – review & editing. IP: Conceptualization, Methodology, Writing – review & editing. ER: Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by the German Research Foundation (Deutsche Forschungsgemeinschaft), being part of the project “the Functional Frontier among Mediterranean and Eurosiberian Plant Communities” (ECOFUMER, 441909701). ER and JG are respectively supported by a María Zambrano and a Margarita Salas fellowships funded by the Spanish Ministry of Universities and European Union-Next Generation Plan. IP acknowledges funding from a Ramón y Cajal contract (RYC2021-033081-I) funded by the Ministry of Science and Innovation and co-funded by the European Union-Next Generation Plan funded by European Union-NextGenerationEU.
Acknowledgments
The authors of this study would like to thank Mario Vega and all the technical staff of the “Hayedo de Montejo” for their always good willingness to facilitate our works in the forest and the interest shown in this study. We also thank Guillermo G. Gordaliza, Eva Miranda and Carlos Magro for data on FIs. Many thanks to Claudia Buchwald for her help and guidance in the lab analysis carried out for this research. Also thanks to Diego Alejandro Melo Prieto for his great help both in the field and in the lab, and to Pardis Golabvand and Marina Tsioli for their help with the lab analyses. We would also like to thank Professor Dr. Francisco Lloret for his advice on field ecological methodology and to his biology student at the time of the performance of this study Anna Ticó i Pifarré for her great help and support with the field work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1303022/full#supplementary-material
Abbreviations
CWM, community-weighted mean; CWS, community-weighted skewness; CWK, community-weighted kurtosis; LES, Leaf Economics Spectrum; RES, Root Economics Spectrum; HyArq, hydraulic architecture; ITV, intraspecific trait variability; SOM, Soil organic matter.
References
Albert, C. H., Grassein, F., Schurr, F. M., Vieilledent, G., Violle, C. (2011). When and how should intraspecific variability be considered in trait-based plant ecology? Perspect. Plant Ecol. Evol. System. 13 (3), 217–225. doi: 10.1016/j.ppees.2011.04.003
Albert, C. H., Thuiller, W., Yoccoz, N. G., Soudant, A., Boucher, F., Saccone, P., et al. (2010). Intraspecific functional variability: extent, structure and sources of variation. J. Ecol. 98 (3), 604–613. doi: 10.1111/j.1365-2745.2010.01651.x
Barton, K., Barton, M. K. (2015). Package ‘MuMIn’; Version 1.43.17. Available at: https://cran.r-project.org/web/packages/MuMIn/MuMIn.pdf (Accessed January 5th, 2023).
Bello, J., Hasselquist, N. J., Vallet, P., Kahmen, A., Perot, T., Korboulewsky, N. (2019). Complementary water uptake depth of Quercus petraea and Pinus sylvestris in mixed stands during an extreme drought. Plant Soil 437, 93–115. doi: 10.1007/s11104-019-03951-z
Benito Garzón, M., Sánchez de Dios, R., Sainz Ollero, H. (2008). Effects of climate change on the distribution of Iberian tree species. Appl. Vege. Sci. 11 (2), 169–178. doi: 10.3170/2008-7-18348
Bílá, K., Moretti, M., de Bello, F., Dias, A. T., Pezzatti, G. B., Van Oosten, A. R., et al. (2014). Disentangling community functional components in a litter-macrodetritivore model system reveals the predominance of the mass ratio hypothesis. Ecol. Evol. 4 (4), 408–416. doi: 10.1002/ece3.941
Box, E. O. (2021). Why does Mediterranean vegetation seem so diverse? Tools Landscape-Scale Geobot. Conserv., 287–313. doi: 10.1007/978-3-030-74950-7_15
Brice, M. H., Vissault, S., Vieira, W., Gravel, D., Legendre, P., Fortin, M. J. (2020). Moderate disturbances accelerate forest transition dynamics under climate change in the temperate–boreal ecotone of eastern North America. Global Change Biol. 26 (8), 4418–4435. doi: 10.1111/gcb.15143
Cantoral, A. L., Alonso-Redondo, R., García-González, M. E., Penas, Á., del Río, S. (2023). Phytosociological analysis of the endemic Quercus faginea forests of the Iberian Peninsula. Plant Biosystems-An Int. J. Deal. All Aspects Plant Biol. 1-51, 419–436. doi: 10.1080/11263504.2023.2165554
Carlson, B. Z., Renaud, J., Biron, P. E., Choler, P. (2014). Long-term modeling of the forest–grassland ecotone in the French Alps: implications for land management and conservation. Ecol. Appl. 24 (5), 1213–1225. doi: 10.1890/13-0910.1
Carmona, C. P., de Bello, F., Azcárate, F. M., Mason, N. W., Peco, B. (2019). Trait hierarchies and intraspecific variability drive competitive interactions in Mediterranean annual plants. J. Ecol. 107 (5), 2078–2089. doi: 10.1111/1365-2745.13248
Carrascosa, A., Silvestre, M., Morgado, L., Azcárate, F. M. (2023). Diversity and mean specific leaf area of Mediterranean woody vegetation changes in response to summer drought across a double stress gradient: The role of phenotypic plasticity. J. Vege. Sci. 34 (2), e13180. doi: 10.1111/jvs.13180
Carter, J. L., White, D. A. (2009). Plasticity in the Huber value contributes to homeostasis in leaf water relations of a mallee Eucalypt with variation to groundwater depth. Tree Physiol. 29 (11), 1407–1418. doi: 10.1093/treephys/tpp076
Carvajal, D. E., Loayza, A. P., Rios, R. S., Delpiano, C. A., Squeo, F. A. (2019). A hyper-arid environment shapes an inverse pattern of the fast–slow plant economics spectrum for above-, but not below-ground resource acquisition strategies. J. Ecol. 107 (3), 1079–1092. doi: 10.1111/1365-2745.13092
Castaño-Santamaria, J., López-Sánchez, C. A., Obeso, J. R., Barrio-Anta, M. (2021). Structure, environmental patterns and impact of expected climate change in natural beech-dominated forests in the Cantabrian Range (NW Spain). For. Ecol. Manage. 497, 119512. doi: 10.1016/j.foreco.2021.119512
Chelli, S., Ottaviani, G., Simonetti, E., Campetella, G., Wellstein, C., Bartha, S., et al. (2021). Intraspecific variability of specific leaf area fosters the persistence of understorey specialists across a light availability gradient. Plant Biol. 23 (1), 212–216. doi: 10.1111/plb.13199
Cordlandwehr, V., Meredith, R. L., Ozinga, W. A., Bekker, R. M., van Groenendael, J. M., Bakker, J. P. (2013). Do plant traits retrieved from a database accurately predict on-site measurements? J. Ecol. 101 (3), 662–670. doi: 10.1111/1365-2745.12091
Cornwell, W. K., Ackerly, D. D. (2009). Community assembly and shifts in plant trait distributions across an environmental gradient in coastal California. Ecol. Monogr. 79 (1), 109–126. doi: 10.1890/07-1134.1
Danet, A., Anthelme, F., Gross, N., Kéfi, S. (2018). Effects of indirect facilitation on functional diversity, dominance and niche differentiation in tropical alpine communities. J. Vege. Sci. 29 (5), 835–846. doi: 10.1111/jvs.12667
Danso Marfo, T., Datta, R., Vranová, V., Ekielski, A. (2019). Ecotone dynamics and stability from soil perspective: Forest-agriculture land transition. Agriculture 9 (10), 228. doi: 10.3390/agriculture9100228
de Bello, F., Lavorel, S., Albert, C. H., Thuiller, W., Grigulis, K., Dolezal, J., et al. (2011). Quantifying the relevance of intraspecific trait variability for functional diversity. Methods Ecol. Evol. 2 (2), 163–174. doi: 10.1111/j.2041-210X.2010.00071.x
de la Riva, E. G., Marañón, T., Pérez-Ramos, I. M., Navarro-Fernández, C. M., Olmo, M., Villar, R. (2018a). Root traits across environmental gradients in Mediterranean woody communities: are they aligned along the root economics spectrum? Plant Soil 424, 35–48. doi: 10.1007/s11104-017-3433-4
de la Riva, E. G., Marañón, T., Violle, C., Villar, R., Pérez-Ramos, I. M. (2017). Biogeochemical and ecomorphological niche segregation of Mediterranean woody species along a local gradient. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.01242
de la Riva, E. G., Pérez-Ramos, I. M., Tosto, A., Navarro-Fernández, C. M., Olmo, M., Marañón, T., et al. (2016a). Disentangling the relative importance of species occurrence, abundance and intraspecific variability in community assembly: a trait-based approach at the whole-plant level in Mediterranean forests. Oikos 125 (3), 354–363. doi: 10.1111/oik.01875
de la Riva, E. G., Prieto, I., de Tomás Marín, S., Rodríguez-Calcerrada, J., Golabvand, P., Galán Díaz, J. (2023). Living at the edge: the functional niche occupation of woody plant communities in the submediterranean ecotone. Ann. Bot. 132 (3), 471–484. doi: 10.1093/aob/mcad138
de la Riva, E. G., Prieto, I., Marañón, T., Pérez-Ramos, I. M., Olmo, M., Villar, R. (2021b). Root economics spectrum and construction costs in Mediterranean woody plants: The role of symbiotic associations and the environment. J. Ecol. 109 (4), 1873–1885. doi: 10.1111/1365-2745.13612
de la Riva, E. G., Querejeta, J. I., Villar, R., Pérez-Ramos, I. M., Marañón, T., Galán Díaz, J., et al. (2021a). The Economics Spectrum drives root trait strategies in Mediterranean vegetation. Front. Plant Sci. 12, 2662. doi: 10.3389/fpls.2021.773118
de la Riva, E. G., Tosto, A., Pérez-Ramos, I. M., Navarro-Fernández, C. M., Olmo, M., Anten, N. P., et al. (2016b). A plant economics spectrum in Mediterranean forests along environmental gradients: is there coordination among leaf, stem and root traits? J. Vege. Sci. 27 (1), 187–199. doi: 10.1111/jvs.12341
de la Riva, E. G., Violle, C., Pérez-Ramos, I. M., Marañón, T., Navarro-Fernández, C. M., Olmo, M., et al. (2018b). A multidimensional functional trait approach reveals the imprint of environmental stress in Mediterranean woody communities. Ecosystems 21, 248–262. doi: 10.1007/s10021-017-0147-7
de Tomás Marín, S., Rodríguez-Calcerrada, J., Arenas-Castro, S., Prieto, I., González, G., Gil, L., et al. (2023). Fagus sylvatica and Quercus pyrenaica: Two neighbors with few things in common. For. Ecosyst. 10, 100097. doi: 10.1016/j.fecs.2023.100097
Díaz, S., Kattge, J., Cornelissen, J. H., Wright, I. J., Lavorel, S., Dray, S., et al. (2016). The global spectrum of plant form and function. Nature 529 (7585), 167–171. doi: 10.1038/nature16489
Dong, N., Prentice, I. C., Wright, I. J., Evans, B. J., Togashi, H. F., Caddy-Retalic, S., et al. (2020). Components of leaf-trait variation along environmental gradients. New Phytol. 228 (1), 82–94. doi: 10.1111/nph.16558
Dorado-Liñán, I., Zorita, E., Martínez-Sancho, E., Gea-Izquierdo, G., Di Filippo, A., Gutiérrez, E., et al. (2017). Large-scale atmospheric circulation enhances the Mediterranean East-West tree growth contrast at rear-edge deciduous forests. Agric. For. Meteorol. 239, 86–95. doi: 10.1016/j.agrformet.2017.02.029
Enquist, B. J., Norberg, J., Bonser, S. P., Violle, C., Webb, C. T., Henderson, A., et al. (2015). Scaling from traits to ecosystems: developing a general trait driver theory via integrating trait-based and metabolic scaling theories. Adv. Ecol. Res. 52, 249–318. doi: 10.1016/bs.aecr.2015.02.001
Fox, J., Weisberg, S. (2019). An R Companion to Applied Regression. 3rd ed. (Thousand Oaks CA: Sage). Available at: https://socialsciences.mcmaster.ca/jfox/Books/Companion/.
Frenette-Dussault, C., Shipley, B., Meziane, D., Hingrat, Y. (2013). Trait-based climate change predictions of plant community structure in arid steppes. J. Ecol. 101 (2), 484–492. doi: 10.1111/1365-2745.12040
Freschet, G. T., Roumet, C., Comas, L. H., Weemstra, M., Bengough, A. G., Rewald, B., et al. (2021). Root traits as drivers of plant and ecosystem functioning: current understanding, pitfalls and future research needs. New Phytol. 232 (3), 1123–1158. doi: 10.1111/nph.17072
Garnier, E., Cortez, J., Billès, G., Navas, M. L., Roumet, C., Debussche, M., et al. (2004). Plant functional markers capture ecosystem properties during secondary succession. Ecology 85 (9), 2630–2637. doi: 10.1890/03-0799
Garnier, E., Navas, M. L. (2012). A trait-based approach to comparative functional plant ecology: concepts, methods and applications for agroecology. A review. Agron. Sustain. Dev. 32, 365–399. doi: 10.1007/s13593-011-0036-y
Gil, L., Náger, J. A., González Doncel, I., Aranda García, I., Gonzalo Jiménez, J., López de Heredia, U., et al. (2010). El Hayedo de Montejo: Una Gestión Sostenible (Comunidad de Madrid: Universidad Politécnica de Madrid. Ed. Madrid: 1036 Dirección General de Medio Ambiente), ISBN: ISBN: 978-84-451-3218-0.
Grime, J. P. (1998). Benefits of plant diversity to ecosystems: immediate, filter and founder effects. J. Ecol. 86 (6), 902–910. doi: 10.1046/j.1365-2745.1998.00306.x
Gross, N., Börger, L., Soriano-Morales, S. I., Le Bagousse-Pinguet, Y., Quero, J. L., García-Gómez, M., et al. (2013). Uncovering multiscale effects of aridity and biotic interactions on the functional structure of Mediterranean shrublands. J. Ecol. 101 (3), 637–649. doi: 10.1111/1365-2745.12063
Gross, N., Le Bagousse-Pinguet, Y. L., Liancourt, P., Berdugo, M., Gotelli, N. J., Maestre, F. T. (2017). Functional trait diversity maximizes ecosystem multifunctionality. Nat. Ecol. Evol. 1 (5), 132. doi: 10.1038/s41559-017-0132
Gross, N., Le Bagousse-Pinguet, Y., Liancourt, P., Saiz, H., Violle, C., Munoz, F. (2021). Unveiling ecological assembly rules from commonalities in trait distributions. Ecol. Lett. 24 (8), 1668–1680. doi: 10.1111/ele.13789
Guerin, G. R., Gallagher, R. V., Wright, I. J., Andrew, S. C., Falster, D. S., Wenk, E., et al. (2022). Environmental associations of abundance-weighted functional traits in Australian plant communities. Basic Appl. Ecol. 58, 98–109. doi: 10.1016/j.baae.2021.11.008
Guo, A., Zuo, X., Zhang, S., Hu, Y., Yue, P., Lv, P., et al. (2022). Contrasting effects of plant inter-and intraspecific variation on community trait responses to nitrogen addition and drought in typical and meadow steppes. BMC Plant Biol. 22 (1), 90. doi: 10.1186/s12870-022-03486-z
Hejda, M., Štajerová, K., Pergl, J., Pyšek, P. (2019). Impacts of dominant plant species on trait composition of communities: comparison between the native and invaded ranges. Ecosphere 10 (10), e02880. doi: 10.1002/ecs2.2880
Henn, J. J., Buzzard, V., Enquist, B. J., Halbritter, A. H., Klanderud, K., Maitner, B. S., et al. (2018). Intraspecific trait variation and phenotypic plasticity mediate alpine plant species response to climate change. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.01548
Hernández Bermejo, J. E., Costa Tenorio, M., Sáinz Ollero, H., Clemente Muñoz, M. (1982). Catálogo florístico del hayedo de Montejo de la Sierra (provincia de Madrid). Lagascalia 11 (1), 3–65.
Hernández-Lambraño, R. E., de la Cruz, D. R., Agudo, J.Á.S. (2021). Effects of the climate change on peripheral populations of hydrophytes: a sensitivity analysis for European plant species based on climate preferences. Sustainability 13 (6), 3147. doi: 10.3390/su13063147
Hogan, J. A., Xu, H., Baraloto, C. (2023). Intraspecific trait variation and species turnover in successional tropical forests: assessing trait imputation for community-weighted means. Plant Ecol. 224 (5), 463–477. doi: 10.1007/s11258-023-01314-4
Hufkens, K., Scheunders, P., Ceulemans, R. (2009). Ecotones in vegetation ecology: methodologies and definitions revisited. Ecol. Res. 24, 977–986. doi: 10.1007/s11284-009-0584-7
Illuminati, A., Querejeta, J. I., Pías, B., Escudero, A., Matesanz, S. (2022). Coordination between water uptake depth and the leaf economic spectrum in a Mediterranean shrubland. J. Ecol. 110 (8), 1844–1856. doi: 10.1111/1365-2745.13909
Iozia, L. M., Crisafulli, V., Varone, L. (2023). Climatic variations along an aridity gradient drive significant trait intraspecific variability in Mediterranean plant species. J. Arid Environ. 217, 105042. doi: 10.1016/j.jaridenv.2023.105042
Jung, V., Albert, C. H., Violle, C., Kunstler, G., Loucougaray, G., Spiegelberger, T. (2014). Intraspecific trait variability mediates the response of subalpine grassland communities to extreme drought events. J. Ecol. 102 (1), 45–53. doi: 10.1111/1365-2745.12177
Kermavnar, J., Kutnar, L., Marinšek, A. (2022). 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. Plant Ecol. 1-14, 229–242. doi: 10.1007/s11258-021-01203-8
Kichenin, E., Wardle, D. A., Peltzer, D. A., Morse, C. W., Freschet, G. T. (2013). Contrasting effects of plant inter-and intraspecific variation on community-level trait measures along an environmental gradient. Funct. Ecol. 27 (5), 1254–1261. doi: 10.1111/1365-2435.12116
Kitagawa, R., Masumoto, S., Nishizawa, K., Kaneko, R., Osono, T., Hasegawa, M., et al. (2020). Positive interaction facilitates landscape homogenization by shrub expansion in the forest–tundra ecotone. J. Vege. Sci. 31 (2), 234–244. doi: 10.1111/jvs.12818
Lajoie, G., Vellend, M. (2015). Understanding context dependence in the contribution of intraspecific variation to community trait–environment matching. Ecology 96 (11), 2912–2922. doi: 10.1890/15-0156.1
Le Bagousse-Pinguet, Y., Gross, N., Maestre, F. T., Maire, V., de Bello, F., Fonseca, C. R., et al. (2017). Testing the environmental filtering concept in global drylands. J. Ecol. 105 (4), 1058–1069. doi: 10.1111/1365-2745.12735
Le Bagousse-Pinguet, Y., Gross, N., Saiz, H., Maestre, F. T., Ruiz, S., Dacal, M., et al. (2021). Functional rarity and evenness are key facets of biodiversity to boost multifunctionality. Proc. Natl. Acad. Sci. 118 (7), e2019355118. doi: 10.1073/pnas.2019355118
Lepš, J., de Bello, F., Šmilauer, P., Doležal, J. (2011). Community trait response to environment: disentangling species turnover vs intraspecific trait variability effects. Ecography 34 (5), 856–863. doi: 10.1111/j.1600-0587.2010.06904.x
Lopes, T. R., Moura, L. B., Nascimento, J. G., Junior, L. S. F., Zolin, C. A., Duarte, S. N., et al. (2020). Priority areas for forest restoration aiming at the maintenance of water resources in a basin in the Cerrado/Amazon ecotone, Brazil. J. South Am. Earth Sci. 101, 102630. doi: 10.1016/j.jsames.2020.102630
López Santalla, A., Navarro, F. P., Náger, J. A., Sánchez, L. G. (2003). El aprovechamiento tradicional del monte y sus efectos sobre la vegetación en el “Hayedo de Montejo” (Madrid). Cuadernos la Sociedad Española Cienc. Forestales 16), 109–114.
Ma, W. T., Yu, Y. Z., Wang, X., Gong, X. Y. (2023). Estimation of intrinsic water-use efficiency from δ13C signature of C3 leaves: Assumptions and uncertainty. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1037972
Maire, V., Gross, N., Börger, L., Proulx, R., Wirth, C., Pontes, L. D. S., et al. (2012). Habitat filtering and niche differentiation jointly explain species relative abundance within grassland communities along fertility and disturbance gradients. New Phytol. 196 (2), 497–509. doi: 10.1111/j.1469-8137.2012.04287.x
Mencuccini, M., Rosas, T., Rowland, L., Choat, B., Cornelissen, H., Jansen, S., et al. (2019). Leaf economics and plant hydraulics drive leaf:wood area ratios. New Phytol. 224 (4), 1544–1556. doi: 10.1111/nph.15998
Meyer, P., Tabaku, V., v. Lüpke, B. (2003). Die Struktur albanischer Rotbuchen-Urwälder–Ableitungen für eine naturnahe Buchenwirtschaft: structural characteristics of Albanian Beech (Fagus sylvatica L.) Virgin Forests–Deductions for semi-natural forestry. Forstwissenschaftliches Centralblatt 122 (1), 47–58. doi: 10.1046/j.1439-0337.2003.02041.x
Moreno, J.Á.C., de las Heras Puñal, P., Estébanez, N. L., Martín, M. J. R. (2005). Caracterización paisajística y ecológica de la Sierra del Rincón (Madrid) (Centro de Investigaciones Ambientales de la Comunidad de Madrid “Fernando González Bernáldez”. Ed. Comunidad de Madrid). doi: 10.13140/2.1.4707.2322
Navarro-Fernández, C. M., Pérez-Ramos, I. M., de la Riva, E. G., Vera, J. R., Roumet, C., Villar, R., et al. (2016). Functional responses of Mediterranean plant communities to soil resource heterogeneity: a mycorrhizal trait-based approach. J. Vege. Sci. 27 (6), 1243–1253. doi: 10.1111/jvs.12446
Niu, K., Zhang, S., Lechowicz, M. J. (2020). Harsh environmental regimes increase the functional significance of intraspecific variation in plant communities. Funct. Ecol. 34 (8), 1666–1677. doi: 10.1111/1365-2435.13582
Pardo, F., Gil, L. (2005). The impact of traditional land use on woodlands: a case study in the Spanish Central System. J. Historic. Geogr. 31 (3), 390–408. doi: 10.1016/j.jhg.2004.11.002
Pereira, S. C., Carvalho, D., Rocha, A. (2021). Temperature and precipitation extremes over the Iberian Peninsula under climate change scenarios: a review. Climate 9 (9), 139. doi: 10.3390/cli9090139
Pérez-Harguindeguy, N., Diaz, S., Garnier, E., Lavorel, S., Poorter, H., Jaureguiberry, P., et al. (2013). New handbook for standardized measurement of plant functional traits worldwide. Aust. J. Bot. 61, 167–234. doi: 10.1071/BT12225_CO
Pérez-Ramos, I. M., Roumet, C., Cruz, P., Blanchard, A., Autran, P., Garnier, E. (2012). Evidence for a ‘plant community economics spectrum’ driven by nutrient and water limitations in a Mediterranean rangeland of southern France. J. Ecol. 100 (6), 1315–1327. doi: 10.1111/1365-2745.12000
Petruzzellis, F., Palandrani, C., Savi, T., Alberti, R., Nardini, A., Bacaro, G. (2017). Sampling intraspecific variability in leaf functional traits: Practical suggestions to maximize collected information. Ecol. Evol. 7 (24), 11236–11245. doi: 10.1002/ece3.3617
Pichon, N. A., Cappelli, S. L., Allan, E. (2022). Intraspecific trait changes have large impacts on community functional composition but do not affect ecosystem function. J. Ecol. 110 (3), 644–658. doi: 10.1111/1365-2745.13827
Pogue, D. W., Schnell, G. D. (2001). Effects of agriculture on habitat complexity in a prairie-forest ecotone in the Southern Great Plains of North America. Agricult. Ecosyst. Environ. 87 (3), 287–298. doi: 10.1016/S0167-8809(01)00150-5
Prieto, I., Querejeta, J. I., Segrestin, J., Volaire, F., Roumet, C. (2018). Leaf carbon and oxygen isotopes are coordinated with the leaf economics spectrum in Mediterranean rangeland species. Funct. Ecol. 32 (3), 612–625. doi: 10.1111/1365-2435.13025
Prieto, I., Roumet, C., Cardinael, R., Dupraz, C., Jourdan, C., Kim, J. H., et al. (2015). Root functional parameters along a land-use gradient: evidence of a community-level economics spectrum. J. Ecol. 103 (2), 361–373. doi: 10.1111/1365-2745.12351
Querejeta, J. I., Prieto, I., Torres, P., Campoy, M., Alguacil, M. M., Roldán, A. (2018). Water-spender strategy is linked to higher leaf nutrient concentrations across plant species colonizing a dry and nutrient-poor epiphytic habitat. Environ. Exp. Bot. 153, 302–310. doi: 10.1016/j.envexpbot.2018.06.007
Querejeta, J. I., Ren, W., Prieto, I. (2021). Vertical decoupling of soil nutrients and water under climate warming reduces plant cumulative nutrient uptake, water-use efficiency and productivity. New Phytol. 230 (4), 1378–1393. doi: 10.1111/nph.17258
Rahman, I. U., Afzal, A., Iqbal, Z., Hashem, A., Al-Arjani, A. B. F., Alqarawi, A. A., et al. (2021). Species distribution pattern and their contribution in plant community assembly in response to ecological gradients of the ecotonal zone in the Himalayan region. Plants 10 (11), 2372. doi: 10.3390/plants10112372
Rawal, R. S., Rawal, R., Rawat, B., Negi, V. S., Pathak, R. (2018). Plant species diversity and rarity patterns along altitude range covering treeline ecotone in Uttarakhand: conservation implications. Trop. Ecol. 59 (2), 225–239.
R Core Team (2021). R: A Language and Environment for Statistical Computing (Vienna, Austria: R Foundation for Statistical Computing), ISBN: ISBN: 3-900051-07-0.
Reich, P. B. (2014). The world-wide ‘fast–slow’ plant economics spectrum: a traits manifesto. J. Ecol. 102 (2), 275–301. doi: 10.1111/1365-2745.12211
Rivas-Martínez, S., Penas, Á., Díaz González, T. E., Ladero Álvarez, M., Asensi Marfil, A., Díez Garretas, B., et al. (2011). Mapa de series, geoseries y geopermaseries de vegetación de España (Memoria del mapa de vegetación potencial de España). Parte II. Itinera Geobot. 18 (2), 425–800.
Rodríguez-Calcerrada, J., Nanos, N., Del Rey, M. C., de Heredia, U. L., Escribano, R., Gil, L. (2011). Small-scale variation of vegetation in a mixed forest understorey is partly controlled by the effect of overstory composition on litter accumulation. J. For. Res. 16 (6), 473–483. doi: 10.1007/s10310-010-0237-2
Rodríguez-Calcerrada, J., Salomón, R., Barba, J., Gordaliza, G. G., Curiel Yuste, J., Magro, C., et al. (2019). Regeneration in the understory of declining overstory trees contributes to soil respiration homeostasis along succession in a sub-Mediterranean beech forest. Forests 10 (9), 727. doi: 10.3390/f10090727
Rubio-Cuadrado, Á., Camarero, J. J., Gordaliza, G. G., Cerioni, M., Montes, F., Gil, L. (2020). Competition overrides climate as trigger of growth decline in a mixed Fagaceae Mediterranean rear-edge forest. Ann. For. Sci. 77, 1–18. doi: 10.1007/s13595-020-01004-5
Rubio-Cuadrado, Á., Camarero, J. J., Rodríguez-Calcerrada, J., Perea, R., Gómez, C., Montes, F., et al. (2021). Impact of successive spring frosts on leaf phenology and radial growth in three deciduous tree species with contrasting climate requirements in central Spain. Tree Physiol. 41 (12), 2279–2292. doi: 10.1093/treephys/tpab076
Ruiz-Labourdette, D., Nogués-Bravo, D., Ollero, H. S., Schmitz, M. F., Pineda, F. D. (2012). Forest 980 composition in Mediterranean mountains is projected to shift along the entire elevational gradient 981 under climate change. J. Biogeography 39 (1), 162–176. doi: 10.1111/j.1365-9822699.2011.02592
Ryel, R. J., Ivans, C. Y., Peek, M. S., Leffler, A. J. (2008). Functional differences in soil water pools: a new perspective on plant water use in water-limited ecosystems. Prog. Bot. 69, 397–422. doi: 10.1007/978-3-540-72954-9_16
Sánchez de Dios, R., Benito-Garzón, M., Sainz-Ollero, H. (2009). Present and future extension of the Iberian submediterranean territories as determined from the distribution of marcescent oaks. Plant Ecol. 204 (2), 189–205. doi: 10.1007/s11258-009-9584-5
Schenk, H. J. (2008). The shallowest possible water extraction profile: a null model for global root distributions. Vadose Zone J. 7 (3), 1119–1124. doi: 10.2136/vzj2007.0119
Schneider, C. A., Rasband, W. S., Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9 (7), 671–675. doi: 10.1038/nmeth.2089
Shea, M. E., Mladenoff, D. J., Clayton, M. K., Berg, S., Elza, H. (2022). Pattern of tree species co-occurrence in an ecotone responds to spatially variable drivers. Landscape Ecol. 37 (9), 2327–2342. doi: 10.1007/s10980-022-01485-x
Shipley, B., Vu, T. T. (2002). Dry matter content as a measure of dry matter concentration in plants and their parts. New Phytol. 153 (2), 359–364. doi: 10.1046/j.0028-646X.2001.00320.x
Siefert, A. (2012). Incorporating intraspecific variation in tests of trait-based community assembly. Oecologia 170 (3), 767–775. doi: 10.1007/s00442-012-2351-7
Siefert, A., Ritchie, M. E. (2016). Intraspecific trait variation drives functional responses of old-field plant communities to nutrient enrichment. Oecologia 181, 245–255. doi: 10.1007/s00442-016-3563-z
Siefert, A., Violle, C., Chalmandrier, L., Albert, C. H., Taudiere, A., Fajardo, A., et al. (2015). A global meta-analysis of the relative extent of intraspecific trait variation in plant communities. Ecol. Lett. 18 (12), 1406–1419. doi: 10.1111/ele.12508
Silvertown, J., Araya, Y., Gowing, D. (2015). Hydrological niches in terrestrial plant communities: a review. J. Ecol. 103 (1), 93–108. doi: 10.1111/1365-2745.12332
Spitzer, C. M., Sundqvist, M. K., Wardle, D. A., Gundale, M. J., Kardol, P. (2023). Root trait variation along a sub-arctic tundra elevational gradient. Oikos 2023 (1), e08903. doi: 10.1111/oik.08903
Staszel, K., Lasota, J., Błońska, E. (2022). Effect of drought on root exudates from Quercus petraea and enzymatic activity of soil. Sci. Rep. 12 (1), 7635. doi: 10.1038/s41598-022-11754-z
Team, R. (2021). RStudio: integrated development for R 2020 (Boston, MA: RStudio, PBC). Available at: https://posit.co/.
Tyree, M. T., Ewers, F. W. (1991). The hydraulic architecture of trees and other woody plants. New Phytol. 119 (3), 345–360. doi: 10.1111/j.1469-8137.1991.tb00035.x
Vila-Viçosa, C., Arenas-Castro, S., Marcos, B., Honrado, J., García, C., Vázquez, et al. (2020a). Combining satellite remote sensing and climate data in species distribution models to improve the conservation of Iberian White Oaks (Quercus L.). ISPRS Int. J. Geo-Inform. 9 (12), 735. doi: 10.3390/ijgi9120735
Vila-Viçosa, C., Gonçalves, J., Honrado, J., Lomba, Â., Almeida, R. S., Vázquez, et al. (2020b). Late Quaternary range shifts of marcescent oaks unveil the dynamics of a major biogeographic transition in southern Europe. Sci. Rep. 10 (1), 1–12. doi: 10.1038/s41598-020-78576-9
Violle, C., Enquist, B. J., McGill, B. J., Jiang, L. I. N., Albert, C. H., Hulshof, C., et al. (2012). The return of the variance: intraspecific variability in community ecology. Trends Ecol. Evol. 27 (4), 244–252. doi: 10.1016/j.tree.2011.11.014
Walker, B. (1995). Conserving biological diversity through ecosystem resilience. Conserv. Biol. 9 (4), 747–752. doi: 10.1046/j.1523-1739.1995.09040747.x
Weemstra, M., Freschet, G. T., Stokes, A., Roumet, C. (2021). Patterns in intraspecific variation in root traits are species-specific along an elevation gradient. Funct. Ecol. 35 (2), 342–356. doi: 10.1111/1365-2435.13723
Wieczynski, D. J., Boyle, B., Buzzard, V., Duran, S. M., Henderson, A. N., Hulshof, C. M., et al. (2019). Climate shapes and shifts functional biodiversity in forests worldwide. Proc. Natl. Acad. Sci. 116 (2), 587–592. doi: 10.1073/pnas.1813723116
Wright, I. J., Reich, P. B., Westoby, M., Ackerly, D. D., Baruch, Z., Bongers, F., et al. (2004). The worldwide leaf economics spectrum. Nature 428 (6985), 821–827. doi: 10.1038/nature02403
Zapater, M., Hossann, C., Bréda, N., Bréchet, C., Bonal, D., Granier, A. (2011). Evidence of hydraulic lift in a young beech and oak mixed forest using 18O soil water labelling. Trees 25 (5), 885–894. doi: 10.1007/s00468-011-0563-9
Zhang, J., He, N., Liu, C., Xu, L., Chen, Z., Li, Y., et al. (2020). Variation and evolution of C:N ratio among different organs enable plants to adapt to N-limited environments. Global Change Biol. 26 (4), 2534–2543. doi: 10.1111/gcb.14973
Zhang, Z., Huang, M., Zhao, X., Wu, L. (2019). Adjustments of leaf traits and whole plant leaf area for balancing water supply and demand in Robinia pseudoacacia under different precipitation conditions on the Loess Plateau. Agric. For. Meteorol. 279, 107733. doi: 10.1016/j.agrformet.2019.107733
Keywords: community assembly, community-weighted moments, ecotone, functional structure, intraspecific variability, mixed-species forest, sub-Mediterranean community, trait-based
Citation: de Tomás Marín S, Galán Díaz J, Rodríguez-Calcerrada J, Prieto I and de la Riva EG (2023) Linking functional composition moments of the sub-Mediterranean ecotone with environmental drivers. Front. Plant Sci. 14:1303022. doi: 10.3389/fpls.2023.1303022
Received: 10 October 2023; Accepted: 17 November 2023;
Published: 08 December 2023.
Edited by:
Stef Bokhorst, VU Amsterdam, NetherlandsReviewed by:
Jalal Kassout, National Institute for Agricultural Research, MoroccoMaria Amélia Martins-Loução, University of Lisbon, Portugal
Copyright © 2023 de Tomás Marín, Galán Díaz, Rodríguez-Calcerrada, Prieto and de la Riva. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sergio de Tomás Marín, dtm.sergio@gmail.com
†ORCID: Sergio de Tomás Marín, orcid.org/0000-0003-1549-0017
Javier Galán Díaz, orcid.org/0000-0002-3284-3059
Jesús Rodríguez-Calcerrada, orcid.org/0000-0001-7408-7653
Iván Prieto, orcid.org/0000-0001-5549-1132
Enrique G. de la Riva, orcid.org/0000-0002-3393-8375