Shrub richness is primarily driven by climate conditions in Southwestern European woodlands

The shrub species richness in Spanish forests is mainly linked to climatic variables and the importance of the groups of variables scarcely differs among forest types. Forest surrounding the Mediterranean Basin exhibit the highest levels of shrub richness. Shrub species account for a high proportion of the plant diversity in Spanish forests and are a determinant factor in forest dynamics and ecosystem functionality. To investigate the relative importance of climatic, forest stand features, soil and topographic variables in explaining shrub richness in Spanish forests and if the relative importance of these four groups of variables reflects variations among forest types. We used the Spanish National Forest Inventory and a boosted regression trees approach to identify which climatic, soil, stand and topographic variables (N = 19 variables) are related to the richness of shrub species in Spanish woodlands. The shrub species richness is mainly related to climatic variables followed by soil variables whereas stand and topographic variables play a minor role. The importance of the groups of variables scarcely differs among forest types although forests located around the Mediterranean Sea display the highest levels of shrub richness. Shrub richness in Spain is primarily driven by climatic and soil variables, both at country and forest-type scales. Forests surrounding the Mediterranean Basin account for the highest richness of shrub species but are also those most threatened by global change. Therefore, special attention must be paid to the monitoring and assessment of these forest ecosystems.


Introduction
The Mediterranean basin is considered a hotspot of diversity as it harbors many endemic species (Myers et al. 2000;Hewitt 2011;Molina-Venegas et al. 2017;Di Pasquale et al. 2020). As in other hotspots, the native species of the Mediterranean basin are threatened by global change, although these countries are wealthier and can contribute with more resources behind environmental conservation alleviating the impact of global change on biodiversity (Habel et al. 2019). Therefore, it would seem appropriate to implement and develop programs for monitoring and conservation of biodiversity as well as to identify those factors driving species occurrence and potential threats (Molina-Venegas et al. 2017;Chandler et al. 2017).
Shrub layer is an important component in the provision of different services and functions in forest ecosystems, such as food and refugia for many organisms (Focardi and Tinelli 2005;Gastón et al. 2019), pollination (Rodríguez and Kouki 2015), fruits and berries for human and animal consumption (Miina et al. 2016), restoration of damaged environments (Torroba-Balmori et al. 2015;Alday et al. 2016;López-Marcos et al. 2020), the reduction of run-off (Mohammad and Adam 2010), reversal of desertification (Maestre et al. 2009), carbon storage (Pasalodos-Tato et al. 2015) or provision of bioenergy (González-González et al. 2017). Additionally, shrubs have an important role in the natural forest dynamics, especially during the first stages of tree development. With this regard, several studies have highlighted the importance of shrub presence on the establishment and development of new seedlings in areas with intense, severe drought (Smit et al. 2008;Moreno-Fernández et al. 2019), but also under milder climatic conditions where water stress intensifies from the lack of soil structure (Torroba-Balmori et al. 2015;Alday et al. 2016), although the relationship between shrubs and trees seems to be speciesspecific (Gómez-Aparicio et al. 2004;Madrigal-González et al. 2014;Cruz-Alonso et al. 2020). Under such conditions, shrubs provide shade and protection for seedlings against browsing damage caused by domestic and wild ungulates. The protection of seedlings by shrubs against browsing has also been reported in temperate forests (Jensen et al. 2012). Shrub clearance, however, has long been common practice in Mediterranean areas to reduce competition with trees, to facilitate forest operations, increase the agricultural area, improve grassland for cattle, and mitigate fire risk by reducing fuel load and breaking the vertical continuity of the fuel (Pérez-Ramos et al. 2008;Lasanta et al. 2018). Moreover, under the current context of global warming, the shrub stratum gains importance as shrubs depict more tolerance to drought and aridity than tree species (Liu et al. 2011;Bateman et al. 2018) and can be one of the last barriers against desertification and soil erosion in water-limited environments (Maestre et al. 2009).
Given the importance of shrub species in forest processes and functions (López-Marcos et al. 2019), it is of great interest to identify and describe the composition and dynamics of shrub communities along environmental gradients, as well as to identify the response of shrub species to environmental conditions (Olthoff et al. 2016(Olthoff et al. , 2021Bhatta et al. 2018;Jordan et al. 2020;Das Gupta and Pinno 2020). Due to the importance of the shrub stratum, the assessment of habitats included in the Natura 2000 Network considers the shrub composition as well as its maintenance across different ecosystems (Cantarello and Newton 2008;Silva et al. 2019Silva et al. , 2020. In fact, Member States are obliged by Article 11 of the Habitats Directive to undertake surveys and inventories which should provide the basis for assessment, the shrub strata being a key component. Spain is an important country with regard to the study of plant diversity because apart from being located within a biodiversity hotspot, it contains four biogeographical regions -Atlantic, Alpine, Mediterranean and Macaronesian (see sections below for further description). These biogeographical regions are, in turn, divided into several subclimates, mainly because of the altitudinal ranges that give rise to different forest types. Most of the studies addressing shrub species richness in forests around the Mediterranean basin hotspot have focused on a single biogeographical region (Gracia et al. 2007;Torras et al. 2008;Martín-Queller et al. 2011;Zilliox and Gosselin 2014) or the transition between two regions (Olthoff et al. 2016(Olthoff et al. , 2018(Olthoff et al. , 2021. However, the relationships among forest-related variables as well as forest processes may vary when modifying the study scale and data source (Olthoff et al. 2018;Vicente-Serrano et al. 2020). In this line, Moreno-Fernández et al. (2020) linked shrub richness to stand variables in monospecific forests of the main tree species but did not consider environmental variables such as climatic, soil or topography. In this regard, shrub species and the relationships between environmental/stand predictors and understory species richness can reveal variations among forest types (Härdtle et al. 2003;Duan et al. 2019;López-Marcos et al. 2019). However, as far as we know, no study to date has addressed shrub richness for the whole Spain and the variables driving this ecological pattern using an explicit spatial analysis.
The main aim of this study is to analyse the spatial patterns of shrub richness and to identify the climatic, forest stand features, soil and topographic variables related to shrub species richness in Spanish woodlands. We aim to answer the following research questions: What is the relative importance of climatic, forest demographic, soil and topographic variables in explaining shrub richness in Spanish forests? Does the relative importance of these four groups of variables reflect variations among forests types?

Dataset
In this study we use the Third Spanish National Forest Inventory (NFI), which was carried out between 1997 and 2007. The NFI consists of permanent plots established in forest areas (coverage of tree species greater than or equal to 10 %) on a 1 × 1 km grid. The plots follow a nested design with four concentric subplots of 5, 10, 15 and 25 m radii in which adult trees are measured according to a minimum diameter at breast height (dbh) threshold (7.5, 12.5, 22.5 and 42.5 cm, respectively) (Alberdi et al. 2016;Olthoff et al. 2016Olthoff et al. , 2018. As in most NFIs (Alberdi et al. 2010), the Spanish NFI identifies the most common and relevant shrub species (or taxa) within the 10-m-radius subplot (219 of which there are 169 species, 48 genus, and 2 subfamily subgroups). The shrub species list was elaborated taking into consideration both floristic aspects (such as dominance, distribution area), ecological aspects (such as species bioindicators or naturalness) as well as other factors of interest (such as biomass, wildfires or livestock browsing). Technical limitations, ubiquity, shrub conservation status or invasiveness have also been considered. It is worth noting that, in general, the individual shrubs are identified at species level but when species identification is not possible, only the genus is recorded. The Spanish NFI considers shrubs such as woody perennial plants, generally less than 5 meters in height at maturity and without a definite crown. However, for practical reasons, vigorous herbaceous climbing plants are also included in this list, so they can be recorded. For each NFI plot, we calculated the shrub richness as the sum of all shrub records within the 10-m-radius subplot. We refer to each shrub record gathered in the NFI plots as shrub species regardless of the taxa (genus or species). Appendix 2 Table 5 contains the list of shrub species considered in the Third NFI. Climbing species included in the shrub species list of the NFI, such as Hereda helix L., were not considered in this study. For the purposes of this study, we considered plots with the presence of at least one adult tree (dbh ≥ 7.5 cm) to eliminate the gap between 0 and 7.4 cm (plots without adult trees but with presence of seedlings and/saplings) in the mean plot dbh, obtaining a continuous variable.
Spain contains four biogeographical regions (Roekaerts 2002). The Atlantic region is located in Northern Spain, in areas subjected to the influence of the Atlantic Ocean. The Alpine region is that within and surrounded by the Pyrenees while the Macaronesian region corresponds to the Canary Islands. Finally, the Mediterranean region accounts for most of the area of the country. Due to the large extent of this region there is intra-variability of climate conditions, ranging from arid conditions in southern locations to subalpine and oro-mediterranean conditions in some of the Mediterranean montane systems.
Appendix 2 Table 5 presents the location of the four biogeographical regions as well as maps for mean annual temperature and mean cumulative rainfall derived from World-Clim2 (Fick and Hijmans 2017). The dataset comprised 81,017 plots of which 4.1 % were located in the Alpine, 17.9 % in the Atlantic, 3.0 % in the Macaronesian and 75.0 % in the Mediterranean biogeographical region (Table 1).
We assigned forest type to each plot according to the Spanish National Forest Map (MITECO 2013). However, this map considers a very large number of types (ca. 70), so to facilitate the interpretation and identification of gradients by forest type we merged them into the following 10 broader spectrum forest types: Subalpine and oro-mediterranean conifers, Mediterranean conifers, Macaronesian conifers, Deciduous broadleaves, Evergreen broadleaves, Mixed stands of conifers and broadleaves, Open woodlands, Productive plantations and Unclassified plots.
For each biogeographical region and forest type, we calculated the mean number of shrub species per plot.

Candidate variables to explain the shrub richness in Spain
We considered the following climatic, stand, soil and topographic variables as drivers of shrub richness in Spain (Table 2): * Climatic variables derived from WorldClim2 (Fick and Hijmans 2017), which is available at https:// www. world clim. org: i) mean annual temperature (Temp, °C), ii) mean cumulative rainfall (Rainfall, mm), iii) and iv) temperature and rainfall variability as the ratio of the standard deviation of monthly values and the mean, this is the coefficient of variation (Var_Temp and Var_Rainfall, respectively, adimensional).
* Stand variables calculated from the NFI: v) plot basal area (BA, m 2 ha -1 ), vi) mean plot dbh (DBH, cm), vii) the number of adult trees per ha (N).

Boosted regression trees
We used boosted regression trees (hereafter, BRT) to relate the number of shrub species (Nshrub) in the i-th plot (i = 1, … , 81,017) to the 19 predictors described above (see Table 3 for the mean values and dispersion statistics). BRT is a machine learning method that combines many trees and can be thought of as a regression model in which each term is a tree (Elith et al. 2008). In contrast to other machine learning techniques, such as random forests or bagging, BRT is based on boosting (De'ath 2007). Thus, for each new tree, the focus is on reducing the variation in the response that cannot be explained by previous trees (refer to Elith et al. (2008) for more details). BRT has several advantages over traditional approaches: no data transformation or elimination of outliers is required, it is relatively insensitive to the effects of multicollinearity, and complex nonlinear relationships and interactions among variables can be described (Schmidt and Drake 2011;Dedman et al. 2017). Parameterizing a BRT model requires the tuning of learning rate (the contribution of each tree added in the model), tree complexity (number of nodes in each tree), and bag fraction (percentage of the data used to fit the model at each step) (Elith et al. 2008;Soykan et al. 2014). The first two parameters determine the number of trees in the model. In this work, we used a learning rate of 0.01, a tree complexity of 4 and a bag fraction of 0.6. The number of trees was set when the predictive deviance reached the minimum value for a given learning rate to reduce the overfitting. These values of learning rate and bag fraction ensured that the BRT was fitted with more than 1000 trees (Elith et al. 2008). We selected a Poisson structure of errors because the response variable (Nshrub) is count data. The relative importance of predictors was calculated to estimate the influence of the four major groups of variables on shrub richness (Friedman 2001). Finally, we estimated the percentage of the deviance explained by the BRTs as follows: where MTD is the mean total deviance and RD is the residual deviance. We fitted one BRT for the entire data set and one BRT for each of the above defined forest types. We carried out the BRT analyses using the "dismo" package in R 4.0.5 environment (R Core Team 2021).

Spatial patterns of shrub species in south-western European forests
Around 88 % of the plots contained at least one shrub species. Similar values were found when analyzing the data by biogeographical region, except for the Macaronesian region, where the percentage of plots with shrub species dropped to 75 %. The ten species or genera with the greatest presence We found that the plots with the largest number of shrub species were mainly located in the Mediterranean bioregion (mean = 4.29 species per plot) followed by the Atlantic (3.00), Alpine (2.44) and Macaronesian (1.40) bioregions (Table 1). In Fig. 1 it can be seen that several shrub richness hotspots exist in forests close to the Mediterranean Sea and in southern Spain.
Mediterranean conifers followed by Mixed stands of conifers and broadleaves and by Evergreen broadleaves emerged as the forest types with the highest richness of shrub species (richness above 4 species per plot), whereas the two Macaronesian types (broadleaves and conifers) and the Open woodlands display the lowest number of shrub species per plot,  (Table 3).

Variables driving shrub species occurrence in Spanish forests
The final BRT was fitted with 8100 trees and explained 38 % of the deviance (Table 4). Overall, climatic variables, followed by soil variables, were found to have the greatest influence on shrub richness. The variability of the rainfall (Var_Rainfall, relative importance = 15%), Temperature (12 %) and Altitude (12 %) emerged as the variables with the greatest influence on shrub richness while the soil bulk density, coarse fragments, aspect, tree density and clay content only contributed marginally to the BRT model (less than 2 %) (Fig. 2). The partial plot of the centered fitted function for rainfall variability presented two maxima (Fig. 3), the first one ranged from 0.2 to 0.4 and the second peaked around 0.7. These two peaks coincide spatially with areas of high shrub richness. The highest values for rainfall variability (around 0.7) occur in the South of Spain whereas eastern Spain presented values of between 0.2 and 0.4 (Fig. 5 in Appendix 1).

Fig. 1 Shrub richness in the NFI plots and the biogeographical regions
The fitted function for Temperature reached a maximum at 13-17 °C, i.e., shrub richness increased at warmer locations. The warmest Spanish areas tend to correspond to the South-West and those near the East coast of the Iberian Peninsula as well as in the Balearic Islands (Mediterranean region) and the Canary Islands (Macaronesian region) (Fig. 6 in Appendix 1). Similarly, altitude was negatively related to shrub richness (see Fig. 7 in Appendix 1 for a digital elevation model of the study area). As regards the relationship between rainfall and shrub richness, the fitted function for this variable suggests that shrub richness tended to reach a maximum in forest growing in sites with moderate-low rainfall, around 600-700 mm per year. The highest levels of rainfall mainly occur in northern Spain (Atlantic and Alpine biogeographical regions) and, to a lesser extent, in some mountain ranges of the Mediterranean region ( Fig. 8 in Appendix 1). These forests present low shrub richness values (Fig. 1). The number of shrub species also reached a maximum at moderate-low values of nitrogen (300 cg kg -1 ) and sand (200-300 g kg -1 ). Nitrogen follows a N-S pattern (high values in the North and low in the South) ( Fig. 9 in Appendix 1), which partially overlaps the pattern for shrub richness (Fig. 1). The regions with the highest of shrub species values are located in areas with low or moderate sand content ( Fig. 1 and Fig. 10 in Appendix 1). However, it should be noted that the relative importance of sand is lower than that of the previously-mentioned variables. Finally, shrub richness decreased at higher plot DBH, i.e., shrub richness maximized in young stands.
The interactions between Rainfall variability and Temperature, and Rainfall variability with Rainfall, emerged as the two main interactions (Appendix 3 Table 6). The analysis of these two interactions reveals that shrub richness reached a maximum in forests located in warm and dry sites with low rainfall variability (Fig. 4), which is supported by findings from visual exploration (Fig. 1) and analyses of individual effects, i.e., partial plots (Fig. 3).

Importance of groups of variables by forest type
As regards the BRT by forest type, we found that overall, the climatic variables group explained the most variance for all the forest types followed by the soil variables, except for Mixed stands of conifers and broadleaves and Open woodlands. In these latter two forest types this order was reversed, although the variance explained by soil and

Discussion
This is the first study addressing shrub richness in woodlands at the Spanish national level. Among all the variables studied, our results reveal that climatic conditions (particularly Fig. 3 Partial plots of the centered fitted functions for the predictors with the greatest relative importance of the boosted regression trees for the data as a whole. In brackets, the relative importance of each variable. Var _Rainfall = rainfall variability; Temperature = mean annual temperature (°C); Alt = altitude (m asl); Rainfall = mean cumulative rainfall (mm); Nitrogen (cg/kg); Sand (g/kg); DBH = mean plot dbh (cm) Fig. 4 Three-dimensional partial dependence plots for the two strongest interactions in the boosted regression trees model rainfall variability) followed by soil features correlate best with shrub richness, explaining around 38 % of the variance. Similarly, Corney et al. (2006) found that the field vegetation in temperate forests is primarily structured by geo-climatic factors. In contrast to our findings, Wang et al. (2017) reported that shrub richness in a dryland ecosystem of Northwest China is mainly limited by soil attributes rather than by climatic conditions. Jordan et al. (2020) found that shrub richness is associated with ecohydrological predictors such as slope, elevation, the number of days with dry lower soil layers or the number of days with wet upper soil layers, although these relationships were weak. This contradiction can be attributable to the influence of the spatial scale covered and the climatic heterogeneity existing in the area studied. At a large scale, the effect of climatic variables may predominate whereas the importance of other variables such as edaphic and topographical may become more significant at smaller spatial scales (Olthoff et al. 2016(Olthoff et al. , 2018. Several other studies have also highlighted the role of variation in rainfall as a determining factor in plant richness and floristic composition in rainforests (Oliveira-Filho and Fontes (2000) using the rainfall distribution ratio as the porportion between the mean precipitation of the dry and rainy seasons) and semi-deciduous forests (Kristiansen et al. (2011) using the coefficient of variation of monthly precipitation). The importance of including the variation in rainfall, however, becomes more apparent on analysing the interactions with Temperature and Rainfall, which reveal that shrub species mainly grow in coastal Mediterranean forests with low rainfall and high temperature. The high species richness occurring in sites with these harsh conditions may be attributable to the fact that hydric stress is more limiting to the trees than to the shrubs (Martín-Queller et al. 2011;Bateman et al. 2018). Adaptation to environmental heterogeneity and wildfires in the Mediterranean basin contributes to explaining why this region displays the highest species richness in comparison with other parts of the country (Cowling et al. 1996). Additionally, the relative mildness of past climates may have provided refugia for some species and new ones may have emerged (Postigo Mijarra et al. 2009;Verdú and Pausas 2013;Molina-Venegas et al. 2017).
Nitrogen was the soil variable most strongly correlated to shrub richness. In concordance with this finding, López-Marcos et al. (2019) found a significant and strong association between the hemicryptophytes growing in the understory and soil fertility in mixed pine forests in Spain. Das Gupta and Pinno (2020) found that soil nutrients are linked to understory richness in boreal forests, but that the direction of the effect depends on the type of forest. Similarly, Wang et al. (2017) found that shrub species were negatively related to soil nutrients whereas Härdtle et al. (2003) argued that nutrient supply only has a weak influence on ground species richness in moist German forests. In contrast, our results indicate a non-linear (bell-shape relationship) between soil nitrogen and shrub richness, which suggests that shrub richness decreases beyond moderate nitrogen content values. The difference between our results for nitrogen content and those of previous studies may be explained by i) the wider ecological gradient at a large spatial considered in our study compared to the finer scales with a narrower ecological gradient of previous studies (see Olthoff et al. 2016) and ii) the statistical approach used in this study (boosted regression trees), which allowed the identification of non-linear relationships Our findings suggest that the number of shrub species reaches a maximum at lower DBH values and decreases progressively as the DBH increases. At younger stages of tree development, trees and shrub display specific interactions of facilitation and competition (Gómez-Aparicio et al. 2004;Madrigal-González et al. 2014;Cruz-Alonso et al. 2020). However, when the tree stratum increases in size the shrubs do not interfere with mature trees (Vilà et al. 2007) but the amount of light reaching the lower strata drops and may hamper the development of the shrub layer (Barbier et al. 2008). In this regard, several studies have found that shrub richness is negatively related to tree cover (basal area and/or tree density) (Martín-Queller et al. 2011;Kepfer-Rojas et al. 2017;Moreno-Fernández et al. 2020). However, the contribution of basal area to explaining shrub richness at country scale was relatively minor and the role of tree density was residual. This result agrees with those of Gracia et al. (2007) and Cruz-Alonso et al. (2021) for temperate forests in the central Pyrenees and low land forests in Mediterranean areas, respectively. We hypothesized that the residual effect of basal area on shrub richness can be due to the structural heterogeneity of Mediterranean forests which ranges from openwoodlands or with low densities to temperate and high montane.
The climatic, soil, stand and topographic variables explained around 38 % of the variance, while the remaining unexplained variance may be attributable to unmeasured variables. In this regard, management practices, forest conservation status, microbial diversity as well as disturbances such as wildfires are expected to affect the composition of the shrub stratum (Dias et al. 2016;Silva et al. 2019Silva et al. , 2020Das Gupta and Pinno 2020).
In contrast to the findings reported by Härdtle et al. (2003), our results indicate that the importance of the groups of variables driving shrub richness does not differ markedly among forest types. The number of shrub species, however, differed among the different forest types. The mean shrub richness by forest type supports the notion that the Mediterranean region contains the greatest shrub richness. Single species forests of Evergreen broadleaves (e.g. Quercus ilex L., Quercus suber L., Ceratonia siliqua L., Olea europaea) and Mediterranean conifers (e.g. Pinus pinaster mesogeensis Ait., Pinus pinea L., Pinus halepensis Mill., Abies pinsapo Boiss., Juniperus spp) tend to have a higher number of shrub species than forests of Subalpine and oro-mediterranean conifers (e.g., Pinus sylvestris L., P. nigra or Abies alba) or Deciduous broadleaves (e.g., Castanea sativa Mill., Fagus sylvatica L., Betula spp., Quercus petraea (Matt.) Liebl., Quercus robur L.) (Moreno-Fernández et al. 2020). Despite the low tree density in Open woodlands (agrosilvicultural systems mainly dominated by Q. ilex and Q. suber), this forest type presented a low number of shrub species. However, land-use practices such as shrub clearing to promote pasture reduce the shrub cover and therefore species richness (Plieninger et al. 2004;López-Sánchez et al. 2016). As regards the low number of shrubs in Deciduous broadleaf forests, it is worth noting that these stands, such as those of F. sylvatica, are very closed, which hampers the establishment of the shrub cover (Barbier et al. 2008).
Macaronesia is considered to be one of the richest areas in Europe, containing many endemic endangered plant species (Romeiras et al. 2016;Patiño et al. 2016). However, our findings reveal that this region and the forest types found there (Macaronesian broadleaves and conifers) have the lowest shrub richness values of all the regions of the country. This may be because a lot of species belong to the tree stratum according to the NFI classification (e.g., Ilex canariensis Poir., Myrica faya Ait. or Persea indica (L.) Spreng) or to the herb stratum, or do not grow in forests (Arozena and Panareda 2013).
In this study, we used a single inventory, i.e., a static approach, which does not allow demographic changes to be detected. The results identify gaps in our knowledge that must be addressed in future studies, such as identifying changes in shrub species assemblages or the replacement of the tree layer by the shrub stratum as a consequence of more arid conditions. Although NFIs provide the largest source of forest data at country level, they present certain limitations. As in most NFIs, the Spanish NFI only identifies a selection of shrub species (or taxa) according to a national list (Torras et al. 2009;Alberdi et al. 2010Alberdi et al. , 2018Olthoff et al. 2021). Hence, it does not provide a complete inventory and several rare and not frequent species are not recorded, which influences total richness although it should not modify general tendencies. Finally, the fact that some specific information is gathered at genus level restricts more in-depth analyses to establish possible associations between overstory and understory species or the estimation of alpha, beta and gamma diversity. In this regard, future research should aim to infer the difference in shrub richness and diversity derived from the NFI and other data sources.

Conclusions
Our results highlight the importance of considering the spatial patterns of biotic and abiotic variables as well as the interaction between these variables along with non-linear relationships for the assessment of plant-species richness. Shrub richness in Spain is primarily driven by climatic and soil variables, both at country and forest-type scales. Of the four biogeographical regions, the Mediterranean region is that which contains the greatest diversity of shrub species growing in Spanish forests. Additionally, this region, as well as the transition zone between the Mediterranean and Euro-Siberian regions, is prone to suffer the effects of climate change, reduction in rainfall and more frequent occurrence of extreme high-temperature events (Giorgi and Lionello 2008). Therefore, we acknowledge the importance of monitoring programs and specific studies undertaken in Mediterranean forest ecosystems. Moreover, future studies should address the impacts of global change on the shrub stratum in the Mediterranean region.
Appendix 1 Maps of the mean predictors of shrub richness.