Assessing above-ground biomass-functional diversity relationships in temperate forests in northern Mexico

Studies on the relationship between biodiversity and ecosystem productivity have suggested that species richness and functional diversity are the main drivers of ecosystem processes. Several patterns on this relationship have been found, including positive, unimodal, negative, and neutral trends, keeping the issue controversial. In this study, taxonomic diversity and functional diversity as drivers of above-ground biomass (AGB) were compared, and the mechanisms that influence biomass production were investigated by testing the complementarity and the mass-ratio hypotheses. Using data from 414 permanent sample plots, covering 23% of temperate forest in the Sierra Madre Oriental (México), we estimated the above-gound biomass (AGB), taxonomic and functional diversity indices, as well as community weighted mean values (CWM) for three functional traits (maximum height, leaf size and wood density) for trees ≥7.5 cm DBH, in managed and unmanaged stands. To compare taxonomic diversity differences between managed and unmanaged stands we carried out a rarefaction analysis. Furthermore, we evaluated the relationship between AGB and taxonomic and functional diversity metrics, as well as CWM traits throught spatial autoregressive models. We found a hump-shaped relationship between AGB and species richness in managed and unmanaged forests. CMW of maximum height was the most important predictor of AGB in both stands, which suggested that the mechanism underlaying the AGB-diversity relationship is the dominance of some highly productive species, supporting the mass-ratio hypothesis. Above-ground biomass was significantly correlated with three of the five functional diversity metrics, CWM maximum height and species richness. Our results show the importance of taking into account spatial autocorrelation in the construction of predictive models to avoid spurious patterns in the AGB-diversity relationship. Species richness, maximum height, functional richness, functional dispersion and RaoQ indices relate with above-ground biomass production in temperate mixed-species and uneven-aged forests of northern Mexico. These forests show a hump-shaped AGB-species richness relationship. Functional diversity explains better AGB production than classical taxonomic diversity. Community weighted mean traits provide key information to explain stand biomass in these forests, where maximum tree height seems to be a more suitable trait for understanding the biomass accumulation process in these ecosystems. Although the impact of forest management on biodiversity is still debated, it has not changed the AGB-diversity relationships in the forests of the Sierra Madre Occidental, Mexico.


Background
One of the most important topics in community ecology has been understanding the mechanisms that underlie ecosystem processes. Particularly, the aboveground biomass (AGB)-biodiversity relationship has been extensively studied in different ecosystems with controversial conclusions. For example, some studies developed across the world with forest inventory data, suggest a positive relationship between AGB and species diversity (Gamfeld et al. 2013;Forrester and Bauhus 2016); however, negative (Reich et al. 2004), unimodal (Wang et al. 2012) and null relationships (Grace et al. 2016) have also been reported.
Two hypotheses have been widely postulated to understand the effects of biodiversity on ecosystem processes, the complementarity hyphotesis proposed by Tilman (1997) and the mass-ratio hypothesis by Grime (1998). The complementarity hyphotesis proposes that species richness and functional diversity, as well as the difference in functional traits between species, maximize the variety of resource use strategies, and that less overlap of niches along the resource axes increases the functioning of the ecosystem through the partition of the niche. In natural ecosystems, it has been suggested that at very low levels of productivity, species richness is limited, and only a few species are adapted to survive under nutrientpoor conditions. Conversely, at high levels of productivity, competition for light and space leads to competitive exclusion of species, with only a few been dominants. On the other hand, Grime (1998) postulated that the productivity of a community is determined by the relative abundance of the functional traits values of the dominant contributors to the plant biomass. According to this hypothesis, ecosystem properties such as biomass should be predictable from the weighted mean community values (CWM) of traits with links to the capture, use, and release of resources at the individual and ecosystem levels (Díaz et al. 2007). Likewise, the assumptions of Grime's hypothesis consider that functional aspects of the dominant species have a greater influence on ecosystem processes because few highly competitive species increase interspecific competition and, eventually, decrease complementary processes and, in consequence, decrease functional diversity. However, this relationship and their underlying mechanisms are still debated (Kirby and Potvin 2007;Schumacher and Roscher 2009;Schuldt et al. 2014;Isbell et al. 2015;Zhang and Chen 2015;Cheng et al. 2018).
One of the main approaches for understanding the relationship between species diversity and ecosystem productivity has been through species richness (Díaz and Cabido 2001). Although species richness might provide a reasonable measure, many of the processes that regulate the functioning of ecosystems are difficult to assign to each species, and it is not possible to determine their relative contributions to ecosystem functioning and productivity. Accordingly, species richness assumes that all species are equivalent and does not consider functional differences among species (i.e. functional diversity), key determinants of ecosystem functioning (Tilman 2001).
The importance of functional diversity (FD) in explaining AGB production has recently been recognized and debated Conti and Díaz 2013;Becknell and Powers 2014). Functional diversity is a component of biodiversity that encompasses the characteristics of species and predicts more accurately the functioning of the ecosystem (Petchey and Gaston 2002;Hooper et al. 2005;Cianciaruso et al. 2009). Moreover, FD links morphological, physiological and phenological variation at the individual level with ecosystem processes and patterns ). Functional diversity metrics are based on traits, which are correlated with environmental filters. In relation to the functioning of organisms, the traits modulate their suitability, allow understanding the interactions between organisms and the components of their environment, and link the different levels of organization of ecological systems (McGill et al. 2006;Violle et al. 2007).
Some authors have proposed that functional diversity and functional traits may be more reliable predictors of ecosystem productivity than taxonomic diversity, which represents species-based diversity indices (including species richness, Shannon and Simpson diversity indices), as plant functional traits are often used as proxies to determine whether species have different ecological strategies for resource acquisition (Vitousek and Hooper 1993;Lawton and Brown 1994;Cornelissen et al. 2003;McGill et al. 2006). However, despite efforts to link functional diversity of plant communities to productivity, the relationships are still not entirely clear.
Most of the studies on this topic have been carried out in relatively simple ecosystems, such as grasslands and monocultures; however, no records have been found in "natural" forest ecosystems with a long history of management as the uneven-aged and multispecies temperate forests of Durango, Mexico. The forests of Durango are the most important forests in Mexico because of their extent and economical value. These forests are of particular interest, not only because they represent a unique ecosystem (with 20 pine and 43 oak species and some species of the genera Arbutus, Fraxinus, Cupressus, Juniperus and occasionally Abies and Pseudotsuga) (González-Elizondo et al. 2012), but also because they are owned and managed by local communities in the form of collective land ownerships known as "Ejidos", or indigenous land holdings known as "Comunidades" (Madrid et al. 2009). Most of the forests in Durango are uneven-aged, and conifer species occur as mixtures with broadleaf species. This irregularity refers to the spatial arrangement of trees (vertical and horizontal irregularity) and the variation in the age structure of trees and stands. The management of these mixed uneven-aged forests is more complex than the management of even-aged forests, and one of the most pressing research problems is the development of models to define sustainable harvests. In that sense, no study has attempted to evaluate the effect of tree species diversity on forest productivity.
In this study, we address two topics widely debated in the current literature. On the one hand, we test the predictive power of taxonomic and functional diversity on above-ground biomass prediction in a managed forest landscape. On the other hand, based on the calculation of the taxonomic and functional diversity metrics in managed and unmanaged forests, we evaluate the effect of the complementarity and the mass-ratio hypotheses. In particular, we tested if functional diversity metrics will be better predictors of forest biomass, and whether, since forests are managed under a wood-oriented forestry model, the dominant species will play a predominant role in the AGB-diversity relationship of the ecosystem (mass-ratio hypothesis).

Study area and inventory data
The study was conducted in the temperate forests of Durango, Mexico (22°20′49″ to 26°46′33″ N; 103°46′ 38″ to 107°11′36″ W), which occupies about 23% of the Sierra Madre Occidental mountain range (Fig. 1). The elevation above sea level varies between 1200 and 3200 m (average 2264 m). Precipitation ranges from 443 to 1452 mm, with an anual average of 917 mm, whereas the mean annual temperature varies from 8.2 to 26.2°C, with an anual average of 13.3°C (Silva-Flores et al. 2014). Dominat vegetation types are pine-oak forests. Some major tree groups occur in these forests, with important tree families including Pinaceae, Fagaceae, and Betulaceae; dominant shrub families include Rosaceae and Ericaceae.
Durango is the most important forestry state in Mexico, with 4.9 million ha of natural forests (mostly managed) and an average timber production of 1.8 million m 3 in the last decade (29.7% of the national forest production) (SEMARNAT 2016). These forests are shaped by more than one century of timber-oriented forestry, mainly by selective removals for sustainable timber production, but also for the maintenance of biological diversity and uneven-aged stand structures (Padilla-Martínez et al. 2020).
The data were obtained from 414 permanent sample plots (PSP) used to monitor the growth and yield of Durango's forests (Crecente-Campo et al. 2014), of which 365 are placed in managed stands and 49 in unmanaged stands. The plots were mostly measured between 2011 and 2015, and cover the main forest types and the current diameter distributions of commercial forests in Durango. The plots are 50 m × 50 m in size and were distributed by systematic sampling, with a variable grid ranging from 3 to 5 km. The following variables by species were recorded: diameter at breast height (DBH, measured in cm at 1.3 m above ground level), total tree height (m), height to the live crown (m), azimuth (°) and radius (m) from the centre of the plot to all trees equal or larger than 7.5 cm on diameter at breast height. The following stand variables were calculated from the tree data recorded in each plot: stems per hectare (N), stand basal area (G, m 2 •ha − 1 ), quadratic mean diameter (dg, cm), and dominant height (H0, m). Dominant height was calculated from the proportion of the 100 pine trees with the largest diameter per hectare. The PSP represent two broad conditions based on management history and structural attributes: (1) unmanaged stands (more than 50 years of non-management history); and (2) managed stands, with components of structural heterogeneity in form of multi-layered canopy.

Biomass estimates
Above-ground biomass (AGB) was considered as measure of stand productivity. AGB of live trees in each plot was estimated using species-specific allometric equations based on DBH and total height (h) developed specifically for all major Durango coniferous and broadleaf tree species (Vargas-Larreta et al. 2017). We used the biomass equations reported by Rojas- García et al. (2015) for less frequently occurring species, or for those tree species for which no biomass equations were developed in the study carried out by Vargas-Larreta et al. (2017).

Functional traits
We measured three functional traits for the species present in each individual plot: maximum tree height (Hmax, m), wood density (WD, g•cm − 3 ), and leaf size (LS, cm). Functional traits Hmax and WD are the most important and commonly studied traits, and are accurate predictors of stand biomass (Ruiz-Benito et al. 2014;Wu et al. 2015). Moreover, we selected these traits because they are important drivers of plant growth and survival (Poorter et al. 2008;Wright et al. 2010;Rüger et al. 2012), and hence for standing above-ground biomass, biomass productivity and carbon stocks (Pérez-Harguindeguy et al. 2013). Tree height was measured in the field, whereas the data for wood density were extracted from the literature for the major tree species in Mexico (Ordoñez . LS values were compiled from Farjon and Filer (2013) and García and González (2003).
The CWM represents the most probable attribute that an individual would have if drawn at random from the Rao's index (RaoQ) is one of the most common indices of funtional diversity, and estimates the abundanceweighted variance of the dissimilarities between all species pairs (Leps et al. 2006). As proposed by Mason et al. (2005) and by Villéger et al. (2008), FRic represents the amount of functional space filled by the community, and it was measured using the convex hull volume method (Cornwell et al. 2006). FEve was measured using a minimum spanning tree which links all the species in the multidimensional functional space (Villéger et al. 2008). Functional dispersion (FDis) was calculated as the weighted distancies from a weighted centroid in multitrait space using plotspecific trait values (Laliberté and Legendre 2010); and functional divergence index (FDiv), was defined as the variance in the attribute values of the species present at a site, and represents the variation in trait values, weighted by the abundance of each species in the community (Mason et al. 2003). FDiv has no units, and its value ranges between 0 and 1, were 0 = low or no trait diversity, and 1 = high trait diversity. The formula for CWM traits and for all the functional diversity indices described are presented in Table 1.
Functional diversity metrics (FDis, FRic, FEve, FDiv, and RaoQ) and CWM traits (CWM_Hmax, CWM_LS, and CWM_WD) were calculated by using the 'dbFD' function of Laliberté and Legendre (2010) of the FD package in R (R Core Team 2020). The taxonomic diversity indices (Shannon-Wiener index; H′, and Simpson's index; λ) were calculated for each plot with the 'vegan' package in R (R Core Team 2020).

Statistical analyses
To compare observed species richness of woody plants between managed and unmanaged stands, we carried out a rarefaction approach using iNEXT (Colwell et al. 2012;Chao et al. 2014;Hsieh et al. 2016) on individualbased rarefaction for a comparable number of individuals (6926). To carry out rarefaction analysis we use the minimum sample size approach.
As the visual inspection of AGB plotted against taxonomic diversity and functional diversity measures clearly indicated unimodal relationships, we fitted quadratic models to assess how stand biomass related to diversity (S, H′, λ), mean trait values (CWM) and functional diversity of the stands. Due to the presence of spatial autocorrelation among biomass data (Moran's I = 0.151, p < 0.0001), we considered both regression and generalized least squares (GLS) models (Dormann Table 1 CWM and functional diversity indices used to describe the functional structure of the forests evaluated

Index Expression Description
CWM CWM x = ∑ p i x i CWMx is the CWM for the trait x, p i is the relative abundance of species i in the community, and x i is the trait value for the species i.

RaoQ
RaoQ ¼ P S i¼1 P S j¼1 d ij p i p j p i and p j are the relative abundances of species i and j, and the d ij values are the dissimilarities between species i and j in the community.

FEve
FEev ¼ is the Euclidean distance between species i and j, the species involved is branch l, and w i is the relative abundance of species i. PEW is the partial weighted evenness; S is number of species.
a j is the abundance of species j and z j is the distance of species j to the weighted centroid c, where c is the weighted centroid in the idimensional space, and x ij the attribute of species j for trait i.
FDiv FDiv ¼ 2 π arctanð5VÞ with V = ∑ p i (lnx i − ln x) 2 and p i ¼ ai P ai 5 is a scaling factor used to define the index over a range of 0-1, and V is the weighted variance of trait x. Inx i is the trait value for the species i; a i is the relative abundance of species i in the community.

et al. 2007
). For the latter, we included a Gaussian correlation structure based on the distance between plots. Confidence intervals were calculated using bootstrap resampling and goodness-of-fit was assessed through Pearson's R 2 for the regression model, and the squared correlation measure r 2 in the GLS case (Smith 2020). All analyses were performed in R (R Core Team 2020) using package nlme (Pinheiro et al. 2020).

Results
Minimum, maximum and mean values of diversity metrics and functional diversity indices, as well as of AGB for the analysed plots are summarized in Table 2. Species richness was greater in managed stands than in unmanaged ones (average 7.5 and 6.78, respectively); as well as the mean values of Shannon-Wiener (H′) and Simpson (λ) indices (Table 2). In the same sense, the rarefaction analysis showed significant differences in species richness estimations between managed (53.22, CI: 52.7-54.47) and unmanaged stands (46.81, CI: 44.94-48.69). Meanwhile, the CWM mean values of Hmax, LS and WD were higher in managed stands than in unmanaged stands; the same case was for all functional diversity measures (Table 2). Unmanaged stands accumulated 16.9% less AGB than managed stands (Table 2); the minimun AGB value (12.24 Mg•ha − 1 ) was found in the unmanaged stands, whereas the maximum AGB value (378.25 Mg•ha − 1 ) was recorded in managed stands; on average, managed stands had more biomass than unmanaged stands (average 137.8 and 114.42 Mg•ha − 1 , respectively).

AGBtaxonomic diversity relationship
Above-ground biomass values varied considerably with species richness (S), and a significant relationship between AGB and number of species was found for all plots together (Table 3). It was observed a slight trend for the entire data set, through a smooth unimodal curve of AGB respect to species richness (Fig. 2). The relationship was quadratic (2nd order polynomial regression: R 2 = 0.04, p < 0.0001; Table 3). Furthermore, the 2nd order GLS model allowed confirming that AGB was linked to a nonlinear relationship with species richness, showing in both cases the well-known hump-shaped trend between biomass and diversity; i.e., AGB was highest at medium species richness values and then decreased with the stand number of species (Fig. 2). Average AGB among species richness values ranged from 69.6 to 160.6 Mg•ha − 1 , with the highest AGB mean values associated with values from 7 to 9 species per plot.
Above-ground biomass was unrelated to Shannon (H′) and Simpson's (λ) indices (figures not shown); in fact, stand biomass was not affected by any of these two taxonomic diversity indices when all plots were pooled together. AGB variability increased with the Shannon's index value and were higher at medium H′ values (1.3-1.7). Heighest AGB value (378.25 Mg•ha − 1 ) was associated with a Simpson's value of 0.75.

AGBfunctional diversity relationships
The second order GLS model showed that functional diversity metrics outperform traditional taxonomic diversity indices (S, H′, λ), but not CWM traits, to explain the variation of AGB. AGB was significantly related to CWM_Hmax (r 2 = 0.14; p < 0.0001), and CWM_WD (r 2 = 0.02; p < 0.03); no significant relationship between AGB and CWM_LS was observed (Table 3). Above-ground biomass was significantly correlated with three of the five functional diversity metrics when all plots were pooled together ( Table 3). The overall relationship between AGB and functional diversity was quadratic; i.e., biomass increased with functional diversity values and then decreased afterwards. In all cases AGB was highest at medium FD values.
AGB in relation to diversity metrics in managed and unmanaged forest stands AGB in relation to taxonomic diversity Table 4 presents the summary of the 2nd order GLS model fitted to AGB-taxonomic diversity metrics, in both managed and unmanaged stands. Among the taxonomic diversity metrics, only species richness (S) in managed stands was a significant predictor of AGB (Table 4). Figures 3, 4 and 5 show the relationship between AGB and taxonomic diversity obtained from the two fitted models. The 2nd order polynomial regression model showed a unimodal relationship for AGB against species richness, Shannon-Wiener and Simpson indices in unmanaged stands; in managed stands this unimodal trend was observed only for species richness (Fig. 3). However, taking into account spatial autocorrelation the unimodal relationships disappear in all cases, except for species richness.
The highest AGB values were related to species richness values between 7 and 9; 1.5 and 0.75 to Shannon-Wiener and Simpson indices, respectively, in managed and unmanaged stands. The above suggests that, without considering the spatial autocorrelation, the forest management does not affect the AGB-species richness relationship in these forests.

AGB in relation to CWM traits
The AGB-CWM_Hmax relationship was positive and highly significant (Table 4), showing the same trend in managed and unmanaged stands (Fig. 4). With the 2nd order polynomial model, a slight negative relationship between AGB and wood density (CWM_WD) in both type of stands was observed (Fig. 4), but not with the GLS model. The above indicates that stands with hard wood species have lower productivity than stands with soft wood species. In managed and unmanaged stands the AGB-CWM_LS relationship did not show a clear trend.

AGB in relation to functional diversity
Relationships between AGB and functional diversity were significant only for three indices (FDis, Fric and RaQ) in managed stands, and for all five indices in unmanaged stands (Table 4).
All functional diversity metrics, except FDiv, showed a unimodal relationship with AGB in both stand types, when such relationship was described with the 2nd order polynomial model. The metrics FDis, FEve, FRic and RaoQ showed the same behavior with the GLS model. In unmanaged stands, the relationships did not show any clear trends with the GLS model (Fig. 5).

Discussion
Numerous studies have documented the relationships between biodiversity and ecosystems biomass production, with quite different conclussions including positive, negative, hump-shaped and null relationships (Adler et al. 2011;Flynn et al. 2011;Barrufol et al. 2013;Zhang and Chen 2015). Lehman and Tilman (2000) and Barrufol et al. (2013), for instance, point out that plant diversity increases when the productivity of the community also increases, whereas Kenkel et al. (2001) found that increasing plant diversity does not influence forest biomass production. Zhang et al. (2012) performed a meta-analysis of more than 50 studies and found a monotonic increase of forest biomass with plant species diversity. In addition, Vilà et al. (2007) report a positive relationship between species richness and woody productivity across different environmental conditions in Spain. Despite this large number of studies, it  has not been possible to establish a definitive pattern that is valid for the pine-oak forests of the Sierra Madre Occidental (Corral-Rivas et al. 2019). In this study we analyze the effect of taxonomic diversity and functional diversity on above-ground biomass by employing 414 permanent plots within uneven-aged temperate forests. The results exhibited a highly variable AGB-species richness relationship, but we identified the hump-shaped (unimodal) relationship between species richness and the accumulation of total biomass (AGB), which has been frequently observed in mature ecosystems (Mittelbach et al. 2001). It is in line with the most recent global meta-analyses and empirical findings across different biomes, which reveal a consistent positive concave-down effect of biodiversity on forest productivity across the world (Liang et al. 2016). Our findings reveal that, when species richness is relatively low, AGB is also low, but increases to a certain level of species richness where the biomass is maximum; after this optimal point, AGB trends to be lower. AGB values were the highest at medium species richness values (number of species in plots between 7 and 9) and then decreased afterwards. These results are similar to those obtained by Corral-Rivas et al.
(2019), who reported an optimum level of species richness between 5 and 10 species associated with the maximum volume production in temperate forests in northern Mexico. Our findings are also in line with the prevailing conclusions that the AGB-diversity is unimodal, which means that productivity peaks at intermediate levels of diversity (Adler et al. 2011;Barrufol et al. 2013;Zhang and Chen 2015;Ding et al. 2019), but are contrary to Wu et al. (2015), who found that forest biomass in temperate forest in China was positively related (R 2 = 0.08) with especies richness; Ammer (2019), who pointed out that, although it is difficult to generalize a definitive pattern, the shape of the curve relating forest productivity to tree species diversity represents an asymptotic curve, and that it seems to be similar across biomes (Liang et al. 2016). Adler et al. (2011) mentioned that the theoretical basis for a hump-shaped stand productivity-species richness relationship has been challenged, and that methodological differences confound cross-study comparisons of empirical evidences (Gotelli and Colwell 2001). We present empirical evidence that the classical hump-shaped productivity (measured as AGB)species richness may be found in the temperate mixed-species and uneven-aged forests of the Sierra Madre Occidental in Mexico.
The highest biomass level was found close to 380 Mg•ha − 1 for the three taxonomic diversity metrics used, and it is related to 7, 0.9 and 2.1 values of species richness (S), Shannon-Wiener (H′) and Simpson (λ) indices, respectively (Fig. 2). The mechanisms behind these relationships have been discussed in terms of species facilitation and competition. When biomass is relatively low, species richness increases due to interspecific facilitation (Guo 2007); whereas when biomass accumulates to a certain point, competition leads to lower diversity (Weiner and Thomas 2001). This pattern could be consequence of the interspecific interactions at community level, where niche differentiation and facilitative interactions between species improve the utilization of available limiting resources (Wang et al. 2016). In this study, when species richness increased, AGB accumulation had a peak, and then reached a decline. Interestingly, AGB peak corresponded with the average values of species richness, Shanon-Wiener and Simpson indices ( Table 2, Fig. 2). These results are consistent with Wang et al. (2016), who state that, when species diversity reaches a certain degree, the productivity may also peak at intermediate diversity leves due to saturation in resource utilization.
Nevertheless, it is important to point out that this relationship could be influenced by others factors besides diversity, such as environmental variables, whose effects on productivity and species richness change along altitudinal gradients (Wu et al. 2015). In this sense, in additional analysis we observed a unimodal relationship between AGB and altitude (r 2 = 0.02; F 2,411 = 3.31; p < 0.03), a negative trend between AGB and temperature (r 2 = 0.04; F 2,411 = 8.81; p ≤ 0.0001), while the relationship between AGB and precipitation was positive (r 2 = 0.21; F 2,411 = 47.61; p < 0.0001). Because these relationships were not objective of this study, future research should include the impact assess of climatic variables on AGB to confirm these observations. We compare our results with some studies that reported positive or hump-shaped AGB-diversity relationships in terms of correlation degree. We found that the effects of species richness on AGB production yield up 10% across the temperate pine-oak northern Mexico (Table 4), which is greater than 6% reported by Wu et al. (2015) for boreal and temperate forest across China. For similar forest ecosystems to the one considered in this study, low R 2 values have been reported; for example, Cavanaugh et al. (2014), Potter and Woodall (2014), and Wardle et al. (2006), who reported 0.16, 0.00-0.25, and 0.21, respectively. Wu et al. (2015) states that the R 2 values between diversity and biomass reported are generally less than 0.3, and that more than 70% of variation in forest biomass can not be explained by biodiversity, which coincides with the results obtained in this study.
The above discussed confirms that the AGB-species richness relationship is highly complex (Li et al. 2018), and it is because Adler et al. (2011) highlights that productivity is a poor predictor of plant species richness. Such a conclusive statement has led to the use of species-level functional traits or functional-trait based approaches to better link biodiversity with forest productivity.
We evaluated the effect of forest management on diversity and AGB production, contrasting managed and unmanaged stands. We were unable to find a clear influence of management on species richness, CWM traits, or functional diversity values, because minimum and maximum values of species richness, as well as all other taxonomic diversity and functional diversity measures used were similar in both managed and unmanaged stands. These findings extend those of Duguid and Ashton (2013), confirming that there is no clear influence of forest management on overstory plant biodiversity when contrasting managed and unmanaged stands. Our results are opposite to the findings by Paillet et al. (2010), who carried out a study to evaluate biodiversity differences between managed and unmanaged forests in Europe, that species richness was only marginally higher in unmanaged than in managed forests, which means that the effect of management on biodiversity is not clearly identificable.
We attribute the heigher AGB and species richness in managed stands, to harvesting regimes. Selective felling is the dominant management type in the Sierra Madre Occidental; it is a management method that prescribes harvest of selected trees within a planning horizon defined by the cutting cycle (ranging up 70 years). A dominant characteristic of this method is the regulation of harvesting regimes to guarantee the natural state and ecological diversity of the forests (Torres-Rojo et al. 2016). This type of management (known as close-tonature-silviculture) has maintained a diverse structure with small, medium and large trees in a multi-layered canopy, and it has originated uneven-aged multispecies stands. That is compelling evidence that current forest management does not negatively affects the AGBdiversity relationships in these forests.
In the last two decades, there has been a growing consensus regarding the importance of using functional traits as predictors of forest productivity measured through biomass (Díaz and Cabido 2001). Furthermore, it has been widely observed that patterns in functional traits can be used to predict relationships in species communities, because they represent the evolutionary processes of species; and the mean values of the traits at the community level are interpreted as results of the niche complementarity or the mass ratio (Cavender-Bares et al. 2009;Götzenberger et al. 2012;HilleRisLambers et al. 2012). In this sense, this study used functional diversity metrics, as well as the community weighted means (CWM) of three functional traits as drivers of the AGB in managed and unmanaged forest stands. CWM traits better correlated with AGB (Table 4) than did taxonomic diversity or FD metrics, in both managed and unmanaged stands, indicating, as Roscher et al. (2012) highlight, that the temperate forests evaluated containing dominant species with different trait combinations reached high productivity. Our analyses showed that AGB was best explained by CWM_Hmax, which indicates that dominant height is the most important trait for high productivity. This has sense, because the dominant height attained by a forest stand at a specified age (known as site index), is the most widely used indicator of potential productivity in forest ecosystems (Hägglund 1981;Vanclay 1994).
Our results suggest that the mass-ratio and environmental filtering hypotheses are not mutually exclusive (Table 3). However, the mass-ratio hypothesis explains to a greater extent the relationship between biomass and functional diversity, since it postulates that the production rate of an ecosystem process, e.g. AGB, is mainly determined by a set of dominant species (Grime 1988). Therefore, under the assumptions of the hypothesis, that is, considering that the environment is uniform, the woody species community has the highest biomass production by a set of species with dominant traits (Pinus cooperi, P. durangensis, P. arizonica, Quercus sideroxyla, and Cupressus lusitanica in unmanaged stands; Pinus durangensis, Quercus sideroxyla, P. arizonica, P. herreae, and P. cooperi in manage stands) which, in our study, represent 63% and 65% of the total AGB, respectively. Furthermore, the null relationship between AGB and functional evenness (FEve) confirms the low contribution of environmental filtering and the facilitating processes in biomass production. Studies carried out by Vile et al. (2006), and Sonkoly et al. (2019) in grass ecosystems in France and Hungary, respectively, as well as studies carried out in tropical ecosystems in Brazil (Duarte et al. 2010) andCosta Rica (Finegan et al. 2015) and temperate forests in China (Yuan et al. 2018) found similar patterns regarding the mechanisms that influence biomass.
On the other hand, our results suggest that forest management is favoring the permanence and abundance of species with commercial importance. The relationship between forest management and the increase in stand biomass as a consequence of the reduction in intra and interspecific competition has been widely documented (Binkley 2004;Juodvalkis et al. 2005). For example, Pretzsch (2005) found that reducing the basal area of beech stands may, to some extent, increase productivity. However, favoring the establishment of a set of high commercial value species may be reducing functional features, functional diversity, and consequently the mechanisms of environmental filtering. Consequently, the null relationship between biomass and functional diversity metrics could be being masked by forest management.
Finally, we must highlight the importance of taking into account the spatial autocorrelation of biomass between plots (Dormann et al. 2007;Fortin and Dale 2009). As shown in Figs. 3, 4 and 5, the inclusion of spatial autocorrelation reduced the effect of the explanatory variables on AGB, even leading to the elimination of the effect, highlighting the existence of spurious relationships. Additionally, it is observed that the incorporation of the spatial autocorrelation reduces the uncertainty on the relationships, i.e., smaller confidence intervals.

Conclusions
Species richness does not influence biomass production in the temperate mixed-species and uneven-aged forests of northern Mexico. These forests show the classic hump-shaped AGB-species richness relationship, with biomass accumulation increasing at low to intermediate levels of species plant diversity and decreasing at high species richness. Functional diversity explains better AGB production than classical taxonomic diversity metrics. Multi-trait functional diversity, but not classical taxonomic diversity, provides key information to explain stand biomass in these forests, where the single trait maximum tree height seems to be the the better option for understanding the biomass accumulation process in these ecosystems. Altohugh the impact of forest management on biodiversity is still debated, it has not changed the AGB-diversity relationships in the forests of the Sierra Madre Occidental, Mexico.