Next Article in Journal
Genome-Wide Identification and Analysis of Cell Cycle Genes in Birch
Next Article in Special Issue
Climate Effects on Black Spruce and Trembling Aspen Productivity in Natural Origin Mixed Stands
Previous Article in Journal
Influence of Topographic Conditions on Teak Growth Performance in Mountainous Landscapes of Lao PDR
Previous Article in Special Issue
Biscogniauxia Charcoal Canker—A New Potential Threat for Mid-European Forests as an Effect of Climate Change
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Species Mixing Proportion and Aridity Influence in the Height–Diameter Relationship for Different Species Mixtures in Mediterranean Forests

1
Sustainable Forest Management Research Institute, Universidad de Valladolid-INIA, Av. da Madrid 44, 34004 Palencia, Spain
2
ECM Ingeniería Ambiental, S.L., C/Curtidores 17, 34003 Palencia, Spain
3
Integrated Remote Sensing Studio, Department of Forest Resources Management, University of British Columbia, 2424 Main Mall, Vancouver, BC V6T1Z4, Canada
4
Southern Swedish Forest Research Centre, Swedish University of Agricultural Sciences, P.O. Box 190, 23422 Lomma, Sweden
5
School of Computer Science, College of Science and Engineering, National University of Ireland, H91TK33 Galway, Ireland
*
Author to whom correspondence should be addressed.
Adapted from my Ph.D. Thesis, E.T.S de Ingenierías Agrarias, Universidad de Valladolid, 34004 Palencia, Spain.
Submission received: 28 October 2021 / Revised: 15 December 2021 / Accepted: 10 January 2022 / Published: 14 January 2022
(This article belongs to the Special Issue Climate Change Effect on Mixed-Species Forest Management)

Abstract

:
Estimating tree height is essential for modelling and managing both pure and mixed forest stands. Although height–diameter (H–D) relationships have been traditionally fitted for pure stands, attention must be paid when analyzing this relationship behavior in stands composed of more than one species. The present context of global change makes also necessary to analyze how this relationship is influenced by climate conditions. This study tends to cope these gaps, by fitting new H–D models for 13 different Mediterranean species in mixed forest stands under different mixing proportions along an aridity gradient in Spain. Using Spanish National Forest Inventory data, a total of 14 height–diameter equations were initially fitted in order to select the best base models for each pair species-mixture. Then, the best models were expanded including species proportion by area (mi) and the De Martonne Aridity Index (M). A general trend was found for coniferous species, with taller trees for the same diameter size in pure than in mixed stands, being this trend inverse for broadleaved species. Regarding aridity influence on H–D relationships, humid conditions seem to beneficiate tree height for almost all the analyzed species and species mixtures. These results may have a relevant importance for Mediterranean coppice stands, suggesting that introducing conifers in broadleaves forests could enhance height for coppice species. However, this practice only should be carried out in places with a low probability of drought. Models presented in our study can be used to predict height both in different pure and mixed forests at different spatio-temporal scales to take better sustainable management decisions under future climate change scenarios.

1. Introduction

Mixed stands have been evidenced as a smart strategy to adapt forest to climate change since they have been shown to provide positive productive and ecosystem service outcomes compared with monocultures [1,2,3,4,5]. Here, diversity plays a significant role in tree growth, environmental adaptation and improvement of meeting demands for goods and services [6,7,8]. Interactions between species may vary with climatic conditions [9] and, hence, climate–growth sensitivity should be explicitly taken into account in the modelling process to consider the impact of climate change [10]. Hence, development of models that consider known species-mixing effects on tree and stand productivity and the potential to improve the resiliency of forestry under expected climate change is essential to properly design the initial composition and subsequent spatial arrangement and management for mixed-species stands [11].
In this context, estimating diameter at breast height and total tree height is fundamental to both developing and applying many growth and yield models [12,13]. The height–diameter relationship (H–D) allows us to describe stand characteristics and development over time, mean height estimation, stand stability, site index, and growth [14]. Many height–diameter equations have been developed using only diameter at breast height as the independent variable for estimating total height [15,16,17,18,19,20]. However, the height–diameter relationship varies from stand to stand [14], being significantly influenced by other tree and stand characteristics, such as tree health and vigor, site quality and stand development stage, stand density or competition, and species mixture in a stand [12,21,22]. We refer to this type of model as “generalized H–D model” [23]. Thus, the inclusion of additional explanatory variables in the H–D models would make the H–D relationships generalizable over large areas [12,14,24,25,26]. In Europe, generalized height–diameter models have been fitted and used since the 1930s [27,28,29,30,31,32,33,34], though there are still a limited number of studies focused on tree H–D relationship in mixed forests [22,25,35,36,37,38]. In this line, recent studies evidence the strong influence of stand structure and species composition on allometric relationships [39,40,41]. However, the magnitude of variations in allometric relationships due to species-mixing is a subject of intense debate, with evidence of adverse and positive effects that depend on multiple factors operating simultaneously [2,39]. Increasing the accuracy of tree height estimates in mixed stands has important implications, because differences in tree morphology directly affect crown competition, stem volume, biomass production, mechanical stability and wood quality predictions [37,42].
On the other hand, climatic conditions are considered one of the main drivers in the biodiversity–productivity relationships showing a wide range of effects across forests globally [5]. However, for instance resistance and relative resilience to drought in mixed stands might depend on depending on species identity, functional diversity in the mixture, site conditions and the developmental stage of the stands [43]. At the same time, in some recent studies, climate variables such as mean annual temperature and precipitation have been found to have a significant effect on tree allometry [44,45]. Hence, H–D models usually designed to be climate-independent will fail to evaluate future scenarios linked to climate change. Previous studies proposed to include of the concept of limiting resources in the tree allometric relationship. Accordingly, [46] added in the modelling approach terms representing abiotic stress, i.e., water availability and [47] included the plant competition symmetry. Tree interactions might be modulated along stress gradients with facilitation being more common in conditions of high abiotic stress (arid sites) relative to more benign abiotic conditions [46]. Therefore, a trend towards higher H–D relationship in mixed stands compared with pure could be expected in the drier and hotter sites, also observed in drought resilience [43].
In the present context of global warming, it is essential to analyze and understand how the H–D relationship is influenced by mixed-species dynamics and environment conditions in order to provide a future management guideline to adapt a Mediterranean forest to climate change. Thus, the inclusion of additional factors on the H–D relationship in mixed forests, such as species proportions and climatic variables, needs to be further considered and analyzed. In this study, we hypothesized that species identity and species composition would define the H–D relationship between in mixed stands compared with pure. We also hypothesized that tree height in mixtures could be benefited under drought conditions compared with pure stands. In this context, this paper focused on analyzing how the identities of competitors and aridity influence the height–diameter relationship in Mediterranean mixed forests. To achieve that, new generalized H–D models that allow us to estimate tree total height in mixed forests along a gradient of different mixtures and species proportions were fit along an aridity gradient in Spain.

2. Materials and Methods

2.1. Data

In this study, data from the second edition of the Spanish National Forest Inventory (SNFI) were used. The SNFI is a systematic sample where plots are located at the nodes of a kilometer square grid and remeasured every ten years. SNFI plots are composed of four concentric sample circles of 5 m, 10 m, 15 m and 25 m radius, where the diameter and height of all trees over 7.5 cm, 12.5 cm, 22.5 cm and 42.5 cm breast height diameter are recorded, respectively [48]. Stand-level variables were calculated using tree measures and expansion factors. Annual mean temperature (T) and annual total precipitation (P) were derived from Worldclim 2 [49] for the selected SNFI plots (Supplementary Table S1). With the aim of analyzing how species composition influences H–D relationship in mixed stands, a broad range of mixing conditions were considered in our study, selecting both pure and mixed plots (Supplementary Figure S1). We focused on 13 main tree species: Aleppo pine (Pinus halepensis Mill.), Black pine (Pinus nigra Arnold), Maritime pine (Pinus pinaster Ait.), Stone pine (Pinus pinea L.), Scots pine (Pinus sylvestris L.), Pinus mugo (Pinus uncinata Turra), Beech (Fagus sylvatica L.), Portuguese oak (Quercus faginea Lam.), Holm oak (Quercus ilex L.), Sessile oak (Quercus petraea Matt. Liebl.), Pyrenean oak (Quercus pyrenaica Willd.), Pedunculate oak (Quercus robur L.) and Cork oak (Quercus suber L.). First, SNFI mixed plots (composed by two main species) were selected based on the following criteria: (1) the main species are included in the pool of studied species; (2) the proportion of both species (based on basal area) being higher than 90% and (3) the proportion of each main species being higher than 15%. SNFI plots were considered as pure if the proportion of basal area of the main species was 90% or higher. Only pure stands close to mixed stands were taking into account for creating the study dataset in order to guarantee similar site conditions. To achieve that, a buffer of 3 km of diameter was created around the dissolved polygon of mixed plots and pure plots, and final plot selection was made inside the intersection surface. All spatial processes were carried out in QGis software [50]. We discarded all species compositions with less than 50 mixed plots in the dataset to keep similar site conditions, which results in 29 representative species compositions being analyzed.

2.2. Model Fitting

2.2.1. Base H–D Model Selection

Model fitting was carried out in two steps. In a first step, basic height–diameter (H–D) models were fitted for each species and species compositions studied. A total of 14 expressions of H–D equations most used from the forestry literature (Table 1) were considered. These base candidate models were fitted for each species and mixture composition using the optimize.curve_fit function from the scikit-learn [51] library in Python. To validate models and check overfitting, a cross validation bootstrap process was carried out. Data were split into train and test sets in an 80/20 proportion for each pair of mixture composition. The train set was used to obtain parameter estimates, since the test set was used for model evaluation performance. Confidence intervals for each parameter were obtained by using a nonparametric bootstrap procedure, as described by Robinson and Froese [52]. The number of bootstrap replicates was set to 1000. Following this bootstrap procedure, root mean squared error (RMSE) and Akaike Information criterion (AIC) were obtained for each model. The best H–D base model was selected for each species in a mixture composition based on parameter significance and model performance (lower AIC and RMSE). Then, selected H–D base models were expanded to allow the inclusion of competition and climate variables in the H–D relationships.

2.2.2. Extension of the Base H–D Models

In a second step, the previously selected H–D base models were expanded including several independent variables representing stand development, competition, species proportion and climatic conditions. Since we considered these variables to mostly affect asymptotic parameter (βo) of the base H–D models, this parameter was further expanded both as a linear (Equation (1)) and multiplicative (Equation (2)) form as follows:
βo= α0 + α1·BAL + α2·Dq + α3·mi + α4·M + α5·Ho
βo= α0 · BALα1 · Dqα2 · miα3 · Mα4 · Hoα5
where βo is asymptotic parameter of base H–D models presented in Table 1; BAL is the basal area of trees larger than a subject tree, in m2·ha−1; Dq is the quadratic mean diameter, in cm; mi is the species proportion by area of species i, described below; M is the De Martonne aridity index, in mm·°C−1; and Ho is the dominant height, in m.
Note that Dq and Ho were tested both at plot and species (Dqi, Hoi) level since preliminary analysis showed a high correlation between these variables at species level with height estimation, especially for Quercus species.
Then, the best H–D model was selected for each species and mixture, considering variable significance and model parsimony using AIC corrected for small sample size (AICc). After that, we fitted non-linear mixed effects models to consider hierarchical and spatial data structure of forest inventory data. Numerous studies have applied the mixed-effects models to describe H–D relationships and have alleviated this lack of independence of error terms, improving the model fitting and prediction accuracy [53,54,55,56]. To consider random effects in the expanded H–D models, the following procedure was carried out: Firstly, we tested plot as random variable based on a similar study [37]. This variable was included iteratively into the different parameters of the expanded H–D models to determine the best model for each pair of species–species composition to allow their comparison, models were fitted using maximum likelihood (ML). Secondly, among the pool of fitted models, we identified the best expanded mixed-effects H–D model for each pair species-mixture composition considering model parsimony using an information-theoretic approach. The model with the lowest AICc and greater Akaike weight (Wi) was considered the best and most parsimonious model for the observed data relative to the set of alternative models. Finally, the best-fitting expanded H–D functions were re-fitted using the unbiased restricted maximum likelihood method (REML) with R-package ‘nlme’ [53]. To ensure that our modelling strategy has accounted for heteroscedasticity, models were assessed using residuals plots. We assessed the contribution of forest stand and climate variables by looking at the significance of their respective parameters.

2.2.3. Mixing Proportion Influence on the H–D Relationship in Mixed Forest Stands

We used species proportion by area (mi) to characterize stand composition and to identify the significant effects of species interactions on H–D relationships in mixed stands. The species-specific growing space occupied is relevant for calculation of the mixing proportions, stand density and quantification of mixing effects on growth [57,58]. Thus, species mixing proportions might be calculated to avoid bias in the quantification of the net total mixing effect, as well as in the relative importance of under- or overyielding by species, due to differences in the potential densities of the species [59]. Competition equivalence coefficient (CEC) compares species-specific growing space requirements of a species with their value in mixed stands, and it is calculated as the ratio of potential carrying capacity of both species in pure stands (by means of maximum stand density index [60]–SDImax). Since recent studies have shown evidence that SDImax varies with climate [61,62], in this study we obtained climate-dependent CECs (Supplementary Table S2) to calculate species mixing proportions for each species in mixed stands using climate-dependent Maximum Size–Density Relationship (MSDR) models presented by Rodriguez-de-Prado et al. [63].
Therefore, species proportion by area was calculated for each species based on SDI and equivalence coefficients following similar studies [57,77,78,79]:
m i = SDI i SDI i +   SDI j · CEC j , i
CEC j , i = SDI max , j SDI max , i
where mi is the species proportion by area of species i; SDIi is the stand density index of species i; SDIj is the stand density index of species j; CECj,i is the competition equivalence coefficient for species j, taking as reference species i; SDImax,i is the climate-dependent maximum stand density index of species i; and SDImax,j is the climate-dependent maximum stand density index of species j.
To compare the magnitude of the mixture effect on tree height, the M(ratio) was calculated, following a similar approach to that of [41]. For each of the studied species and mixtures, tree height for a given diameter variation according to the species competition environment was estimated for pure stands and mixed stands using the fitted models for a reference tree of 30 cm diameter, and the other independent variables were fixed to the mean (Supplementary Table S1). M(ratio) was calculated as follows:
M ( ratio ) i , j , m i   = H ( mix ) i , j , m i H ( pure ) i , j , m i  
where H(mix) is the estimated mean total height of species i in mixture j, at a species mixing proportion mi, and H(pure) is the estimated mean total height of species i in mixture j at mi = 1. Thus, we quantified the relative change on tree height for a tree of 30 cm diameter in response to inter-specific competition along a gradient of different species mixing proportions (mi from 0.2 to 0.8). Therefore, when M(ratio) > 1, the estimated total tree height would be higher in mixed than in pure stands, and lower tree height is expected in mixed stands when M(ratio) < 1.

2.2.4. Climatic Influence on H–D Relationships in Mixed Forests Stands

The De Martonne Aridity Index (M) was included in the models as a measure of aridity (M = P/(T + 10); mm °C−1) [80]. Including this index into H–D models allowed us to estimate total tree height along an aridity gradient. Accordantly with M(ratio), C(ratio) was here calculated to analyze potential differences of total tree height for a given tree diameter (30 cm) between arid and humid environments. All independent variables (excepting M, and mi = 0.5) were fixed to the mean values for each species and mixture (Supplementary Table S1). C(ratio) was proposed and calculated in this study as follows:
C ( ratio ) i , j   = H ( arid ) i , j H ( humid ) i , j
where H(arid) is the total tree height estimated using fitted models for species i in mixture j along the aridity gradient, from smallest to the highest values of the De Martonne Index. H(humid) is the estimated total tree height predicted at the highest De Martonne Index found. C(ratio) depicts the magnitude of relative tree height for a given diameter higher (C(ratio) > 1) or lower (C(ratio) < 1) in more arid than humid places, respectively.

3. Results

3.1. New H–D Models for Mixed Forests in Spain

Coefficients for the final best generalized mixed-effects H–D models for each species and mixture are presented in Table 2. Results indicated that M1 [64] was the most frequently chosen base H–D model among conifers, whereas M12 [74] was the most chosen among broadleaves. An exception was found for Pinus sylvestris and Pinus uncinata, where M3 [66] and M5 [68] showed better fits for H–D relationships, respectively. All the parameters of the expanded H–D models obtained by bootstrap procedure were significant (p-value < 0.05), finding some differences regarding the independent variables entering the models for coniferous and broadleaves. For the coniferous species, competition (BAL) and species proportion by area (mi) had a significant effect explaining height in the models. For these species, the competition at plot level (by means of Dq) was included in the models more times when mixed with conifers, but Dqi (at species level) was included more often when mixed with broadleaves. Indeed, for all the Quercus species analyzed, Dqi and Hoi were included into the H–D expanded models as independent variables. Note that, for almost all the species and species compositions of broadleaved species, the De Martonne Index (M) was significant with a positive sign, meaning a height increment by humid conditions. For all pairs of species mixtures, expanded H–D models significantly improved the goodness of fit, in terms of AIC and RMSE, compared with the base H–D models. Among the analyzed species, Pinus pinea and Quercus ilex selected models showed the smallest RMSE with values ranging 0.80–1.06 m. On the contrary, the highest RMSE values were found for Pinus sylvestris, Pinus uncinata, Quercus petraea, Quercus pyrenaica and Quercus robur with values close to 1.5 m.

3.2. Species Mixing Influence on H–D Relationships

Figure 1 shows the height variation at the same diameter size under different species mixing proportions in mixed forest stands for the studied species. A general trend was found for Pinus species (Figure 1a–f), with taller trees in pure than in mixed stands (M(ratio) < 1). However, an exception was found in for Pinus sylvestris in Pinus sylvestris-Pinus pinaster mixture (Figure 1e), with slightly higher trees at the same diameter size in the mixture than pure stands of this species. Pinus uncinata height was nearly insensitive to proportional area change and mixed and pure forest allometry (Figure 1f). On the other hand, an inverse trend was found for broadleaved species, where trees were comparatively taller (~2–8%) in mixed stands than in pure ones (Figure 1g–m) at the same diameter size. Among them, Fagus sylvatica, Quercus faginea and Quercus robur height was comparatively stable under different species mixing, finding similar estimated tree heights both in pure and mixed forests (Figure 1g,h,l). Among the conifers, the highest differences in height between mixed and pure stands were observed for Pinus nigra (Figure 1b) when mixed with Quercus faginea and Quercus ilex (~10% reduction in height in mixtures). Regarding broadleaved species, Quercus ilex showed a high difference between mixed and pure stands when mixing with Pinus halepensis (Figure 1i), with taller trees at the same diameter in mixed stands. This trend was also found for the Quercus suber and Pinus pinaster mixture (Figure 1m); Quercus petraea also presented significant differences in height with respect to pure stands when mixed with Fagus sylvatica (Figure 1j).

3.3. Aridity Influence on H–D Relationships

The variation of height at the same diameter size under different aridity conditions, by terms of the De Martonne Aridity Index, is presented in Figure 2. According to our results, taller trees could be found in more humid conditions for almost all the analyzed species and species mixtures. Based on estimated height along an aridity gradient, results also indicated that conifers studied here are less sensitive to changes in aridity than broadleaved species (Figure 2). Among the studied species and species compositions, aridity influence was nearly insignificant for Pinus nigra, Pinus pinaster and Quercus suber (Figure 2b,c,m). Two different trends were found regarding climatic influence on height for conifers. On one hand, Mediterranean conifers such as Pinus halepensis (Figure 2a) and Pinus pinea (Figure 2d) are approximately 5% higher in more arid than humid conditions (C(ratio) < 1) when mixed with Quercus ilex, mainly. Oppositely, estimated heights for species living at higher altitudes such as Pinus sylvestris (Figure 2e) and Pinus uncinata (Figure 2f) presented higher values at more humid environments (C(ratio) > 1), although these differences were nearly insignificant (~1–2%). Broadleaved species (Figure 2g–m) seemed to be highly influenced by aridity according to their height estimations, excepting Quercus suber. Among them, the highest differences in height between arid and humid conditions were found for Quercus faginea in the mixture with Pinus sylvestris (Figure 2h) and Quercus robur in mixture with Fagus sylvatica (Figure 2l), with values close to 15% (height difference in humid compared drought sites). Figure 3 shows total tree height estimations (y-axis) along the diameter distribution of each species in each mixture (x-axis) under different values of dominant height. Here, we can see that trends in aridity and height found for each species and mixture are constant under different tree sizes and stand developmental stages. In the case of Pinus halepensis, taller trees were found in more arid than in less arid places for all sizes and developmental stages (Figure 3a) when mixed with Quercus ilex. On the other hand, height estimations for Quercus robur when mixed with Fagus sylvatica indicated that, under the three simulated stand developmental stages, taller trees were found in more humid places than in more arid places (Figure 3b).
Aridity–total height relationship for different sizes and developmental stages are available for the remaining species of study in Supplementary Figure S2.

4. Discussion

This study presents novel non-linear height–diameter models for predicting tree height in Mediterranean mixed forests under different species mixing proportions and aridity conditions. Here, we evidenced that stand composition, target species, identity of competitors and aridity conditions modulate H–D relationship for the species and mixtures analyzed, resulting in positive, neutral or negative effects on H–D relationships. Compared with existing models to estimate tree height in Spain [81], our work contributes towards H–D models accounting for aridity conditions and species composition in the estimates in addition to tree competition and plot-level factors. Tree height estimates accurately represent species mixing proportions and aridity conditions will improve derived upscaling volume estimations from tree to stand level. Hence, our results could be important to support management and policy decision for Mediterranean forests under the context of climate change, especially for mixed stands.

4.1. Total Tree Height Response to Species Mixing Proportions

The information on tree heights obtained from H–D models is essential in forest inventories for computing tree volumes and evaluate tree slenderness which could be related to stand stability [82]. Our results showed that tree volume seemed to be higher in pure stands for conifer, but with lower stability (high slenderness), increasing the risk of damage in windstorms. On the contrary, slenderness and tree volume increased in mixtures for broadleaves species, which could be key to stimulating tree development for Mediterranean coppice stands. Taking into account species mixing proportions in mixed forest studies is key to understanding differences in tree and stand variables of a given species [79]. In this line, our results suggested that conifers studied here could be taller at the same diameter size in pure than in mixed forest stands (Figure 1a–f). We hypothesized that it may be caused by species-specific traits, such as pines are pioneer and shade intolerant species; thus, they are very sensitive to growing under broadleaves cover. This depends primary on species identity, i.e., shade tolerance is caused for plant ontogeny and influenced by numerous biotic and abiotic factors [83]. Shade tolerance represents how species intercept light and how efficiently it is converted into photosynthates; hence, we would expect that mixtures of different shade tolerance could affect tree light competition, growth and H–D relationship. Accordingly, [84] showed that species-specific shade tolerance may affect tree allometry and be age-mediated.
Among the conifers studied here, Pinus nigra showed the biggest differences in height between mixed and pure stands (Figure 1b), especially when mixed with oaks. It is known that although this is a high-shade-tolerant species compared with other Mediterranean pines, water limitation could drastically limit its distribution in Mediterranean plant communities [85,86]. An exception among conifers studied here was found in Pinus sylvestris mixed with Pinus pinaster, with higher trees at the same diameter size found in mixed than in pure stands, although the M(ratio) keeps close to 1 along all the mixing proportions gradient (Figure 1e). A similar trend was found by Riofrio et al. [37] for this species and mixture composition from SNFI data. Another related explanation to these patterns could be due to differences in species growth rates. When the mixture occurs between two species with different growth rates, the estimated height for the pure stand may be higher at the same stage of development [87]. This seems logical because all trees in a regular pure stand grow very similarly and in a regular way, maintaining the upper average and competing for light availability. This is much more evident for pine in conifer–broadleaved mixtures, where the latter do not grow as much in height and the crown structure is not pyramidal [36], increasing inter-specific competition for pine species. Competition may also explain our results found for Quercus petraea when mixed with Fagus sylvatica (Figure 1j), with great differences in height at the same diameter size between mixed and pure stands for oak species. In mixture, we hypothesized that the growth rate of the light-demanding oak is supposed to be higher than the shade tolerant [3]. Thus, an increase of beech proportion in mixture could reduce intra-specific competition for oak and enhance height for this species. However, this contradicts Ligot et al. [88] observations in forests of Central Europe, where oak trees were systematically outcompeted by beech in the mixture, preventing also oak regeneration. Focusing on broadleaved species, and contrary to conifers, a common trend was found where higher trees at the same diameter size were found in mixed stands (Figure 1g–m). Among other reasons, this inverse trend could be due to stand stratification promoted by silvicultural history. Silviculture can accelerate the growth of different species in mixtures, affecting stand composition, structure and dynamics [89]. In some stands, one of the species in the mixture may historically experience repeated fires, grazing or silvicultural treatments (cutting) to favor the other species. This can be the case of the Quercus-Pinus mixtures analyzed, where H(mix)/H(pure) ratios higher than 1 for broadleaved in mixed stands may suggest a positive effect of the mixture. Note that the proposed M(ratio) index relied on species plasticity in tree allometry in response to competition, which, according to Pretzsch (2014) [2], can be defined as the capacity of the species to change tree allometry under different competitive environments. As broadleaf coppice forests are wide extended in Mediterranean areas, we think that the mixture with conifers may reduce the intra-specific competition and thus may facilitate the apparition of taller trees. In this sense, our results may have a relevant importance for Mediterranean coppice stands, suggesting that introducing or combing conifers in broadleaved forests, could enhance the height–diameter relationship, obtaining taller trees at the same diameter for coppice species. This practice would also increase the ecosystem services and productivity; thus, it has been commonly used in reforestation [90]. Differences in trends between conifers and broadleaved species studied here could be also explained based on initial SNFI selected plots, which could have different stand structures and development stage of the species, as shown in Supplementary Figures S3 and S4. It seems that in conifer mixtures, the size of the trees seems to be different in strict Mediterranean species (Pinus halepensis, Pinus pinaster, Pinus pinea) than in montane species (Pinus sylvestris, Pinus nigra, Pinus uncinata), probably as a consequence of the character of greater shade intolerance. In addition, as the stand distribution for a species moves away from its ideal ecological niche, competition with other more adapted and opportunistic species appears. If the second or third species is equally or better adapted, it leads to the formation of mixed stands. This can be the reason why Quercus species may perform better in mixtures (mainly broadleaved-conifer as Quercus ilex-Pinus halepensis or Quercus suber-Pinus pinaster) than pines, since they naturally form this kind of mixed forests. Differences in water uptake depth appears to be stimulated in mixed forests because Quercus ilex surrounded by Pinus halepensis explores deeper water sources than in monospecific formations [91]. We observed that shade-tolerant species as Fagus sylvatica, Quercus robur and Pinus uncinata, were no affected by mixture, regardless species composition. This may be due to shade tolerant species having high plasticity for certain traits, particularly for morphological features which optimize the light capture [83]; therefore, a high canopy cover promoted by mixed stand seems to be practically irrelevant.

4.2. Total Tree Height Response to Aridity

We observed a significant climatic influence on the H–D relationship for most of the species under study, in line with similar research in the topic [44]. Including climatic related variables in the H–D models could be the key to design forest prescriptions under different present and future environmental conditions [92]. A basic assumption in tree allometry studies is that, as long as site conditions are homogenous, trees of a given species, at a given location, with the same diameter at breast height, would have the same height [78]. However, climatic conditions change over locations and time, and they could contribute towards species-specific dynamics [93,94]. In this context, it is well known that height–diameter relationships depend heavily on local environmental conditions [95]. In the present scenario of climate change, the fluctuating climatic conditions may have a positive or negative influence on the H–D relationship [44]. We observed that for most of the mixtures analyzed here, total tree height along a diameter range was reduced under drought conditions (Figure 2). However, two exceptions opposed to this trend were observed: Pinus halepensis-Quercus ilex (Figure 2a) and Pinus pinea-Pinus pinaster (Figure 2d). Although in the first case it may be due to a reduction in intra-specific competition for Pinus halepensis during drought conditions, in the second case, it could be due to a facilitation process, i.e., a complementarity in the use of resources [96]. Due to its stronger plastic character and its potential adaptability, Pinus halepensis has been demonstrated to be the most suitable species in terms of tree growth in arid sites of Spain [97]. A more efficient use of below-ground water resources due to different root stratification could be also speculated for these mixtures, which allows for a better exploration of the soil profile [91,98]. Therefore, these species could use more resources to strengthen the root system than the aerial system as aridity increases.
Potential changes in height due to changes in aridity may be higher in Mediterranean broadleaved than conifers studied here, suggesting that the allometry of oak species in Spain may be highly affected by future climate change. We hypothesized that the differences between species on water use efficiency may underlay this process. Pinus species can maintain a relatively stable leaf water potential by strict stomatal control during drought events (isohydric strategy). On the other hand, broadleaves do not have a discernible threshold of minimum water potential response, i.e., a light stomatal control under drought conditions (anisohydric behavior) [96,97,98,99]. Consequently, higher competition for water use resources in arid sites is expected for broadleaves when mixed with conifers [100]. Our results are in line with this hypothesis, as shown especially for Quercus faginea when mixed with Pinus sylvestris (Figure 2h). As previously observed in other oak–pine mixture compositions [101], the impact of soil water deficit on species may be worsened in mixtures, but a mixture of species with different growth sensitivities to the seasonality of the drought periods might help to buffer these effects.

5. Conclusions

Data from the Spanish National Forest Inventory (SNFI) were used to estimate the H–D relationship, covering thirteen species in twenty-nine different species compositions in Spain. The modified generalized mixed model was more adequate for modeling the height–diameter relationship of the tree species groups over the base H–D models. We found a significant contribution of climate and stand species composition variables in modifying the H–D relationship. Regarding species mixing proportions, taller trees at the same diameter size were found more commonly in monospecific stands for conifers and mixed stands for the broadleaved species studied here. Based on aridity conditions, in general, taller trees at the same diameter are supposed to be found in more humid conditions, with exceptions found for Pinus pinea and Pinus halepensis species in species-specific mixtures. The broadleaved and conifer mixed stands showed different patterns. Broadleaves species are more sensitive to prone drought sites. The use of forest inventory data allowed us to take into account a wide range of different species and conditions (climatic, soils, competition, species composition or mixtures). H–D models presented in this study permit us to save time and effort by predicting tree heights instead of performing direct measurements. The use of new programming techniques such as machine learning and deep learning applied to more precise data (LiDAR and satellite imaginery) could help us obtaining more precise height estimations in the future within the scope of further research. H–D models presented in our study can be an important tool to estimate height at different spatio-temporal scales in order to take better sustainable management decisions under future climate change scenarios.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/f13010119/s1: Table S1: Stats (mean and range) of the main tree and plot variables characterizing selected SNFI plots used in our study; Table S2: Competition Equivalence Coefficients for the most representative species compositions in Spain; Figure S1: Location of the Spanish National Forest Inventory plots selected to fit the H–D models for the 29 different mixtures analyzed in this study; Figure S2: Aridity influence in terms of the De Martonne Index (M), on H–D relationship for different stand dominant height (Ho); Figure S3: Mean values of mean quadratic diameter by species (Dgi) for the selected SNFI plots for the different conifer–conifer (a), conifer–broadleaved (b) and broadleaved–broadleaved (c) mixtures analyzed; Figure S4: Mean values of dominant height by species (Hoi) for the selected SNFI plots for the different conifer–conifer (a), conifer–broadleaved (b) and broadleaved–broadleaved (c) mixtures analyzed.

Author Contributions

Conceptualization, D.R.d.P., F.B. and C.H.d.A.; methodology, D.R.d.P., J.R. and J.A.; data curation and formal analysis, D.R.d.P. and J.M.; writing—original draft preparation, D.R.d.P.; writing—review and editing, D.R.d.P.; supervision, J.R., J.A., C.H.d.A. and F.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the SPANISH MINISTRY OF ECONOMY AND COMPETITIVENESS (MINECO), through an Industrial PhD Grant (grant number DI-15-07722) and the Project “CLU-2019-01-iuFOR Institute Unit of Excellence” of the University of Valladolid, funded by the Junta de Castilla and co-financed by the European Union (ERDF “Europe drives our growth”). In addition, this project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 734907.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Susana Miguel Barreiro for revising the complete text. In addition, authors would also thank to MSCA-RISE “3D-NEONET: Drug Discovery and Delivery NEtwork for ONcology and Eye Therapeutics” project for the opportunity of working close to the National University of Galway, School of Computer Science in the development and use of new programming techniques applied to sustainable forest management. Finally, the authors would like to thank to the Ministry for Ecological Transition and the WorldClim team for sharing and providing the data used in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Forrester, D.I. The spatial and temporal dynamics of species interactions in mixed-species forests: From pattern to process. For. Ecol. Manag. 2014, 312, 282–292. [Google Scholar] [CrossRef]
  2. Pretzsch, H. Canopy space filling and tree crown morphology in mixed-species stands compared with monocultures. For. Ecol. Manag. 2014, 327, 251–264. [Google Scholar] [CrossRef] [Green Version]
  3. Pretzsch, H.; del Río, M.; Ammer, C.; Avdagic, A.; Barbeito, I.; Bielak, K.; Brazaitis, G.; Coll, L.; Dirnberger, G.; Drössler, L.; et al. Growth and yield of mixed versus pure stands of Scots pine (Pinus sylvestris L.) and European beech (Fagus sylvatica L.) analysed along a productivity gradient through Europe. Eur. J. For. Res. 2015, 134, 927–947. [Google Scholar] [CrossRef] [Green Version]
  4. Pretzsch, H.; Schutze, G. Effect of tree species mixing on the size structure, density, and yield of forest stands. Eur. J. For. Res. 2016, 135, 1–22. [Google Scholar] [CrossRef]
  5. Jactel, H.; Emmanuel, G.; Lars, D.; Forrester, D.; Mason, W.; Morin, X.; Pretzsch, H.; Castagneyrol, B. Positive biodiversity–productivity relationships in forests: Climate matters. Biol. Lett. 2018, 14, 20170747. [Google Scholar] [CrossRef]
  6. Larjavaara, M. A review on benefits and disadvantages of tree diversity. Open. For. Sci. J. 2008, 1, 24–26. [Google Scholar] [CrossRef] [Green Version]
  7. Gamfeldt, L.; Snall, T.; Bagchi, R.; Jonsson, M.; Gustafsson, L.; Kjellander, P.; Ruiz-Jaen, M.C.; Froberg, M.; Stendahl, J.; Philipson, C.D.; et al. Higher levels of multiple ecosystem services are found in forests with more tree species. Nat. Commun. 2013, 4, 1340. [Google Scholar] [CrossRef] [PubMed]
  8. Ivetić, V.; Devetaković, J.; Nonić, M.; Stanković, D.; Šijačić-Nikolić, M. Genetic diversity and forest reproductive material-from seed source selection to planting. iForest 2016, 9, 801. [Google Scholar] [CrossRef] [Green Version]
  9. Condés, S.; del Río, M. Climate modifies tree interactions in terms of basal area growth and mortality in monospecific and mixed Fagus sylvatica and Pinus sylvestris forests. Eur. J. For. Res. 2015, 134, 1095–1108. [Google Scholar] [CrossRef]
  10. Bosela, M.; Merganičová, K.; Torresan, C.; Cherubini, P.; Fabrika, M.; Heinze, B.; Höhn, M.; Kašanin-Grubin, M.; Klopčič, M.; Mészáros, I.; et al. Modelling Future Growth of Mountain Forests Under Changing Environments. In Climate-Smart Forestry in Mountain Regions; Tognetti, R., Smith, M., Panzacchi, P., Eds.; Springer: Cham, Switzerland, 2022; Volume 40. [Google Scholar] [CrossRef]
  11. Pretzsch, H.; Zenner, E.K. Toward managing mixed-species stands: From parametrization to prescription. For. Ecosyst. 2017, 4, 19. [Google Scholar] [CrossRef] [Green Version]
  12. Sharma, M.; Yin Zhang, S. Height–diameter models using stand characteristics for Pinus banksiana and Picea mariana. Scand. J. Forest. Res. 2004, 19, 442–451. [Google Scholar] [CrossRef]
  13. Blozan, W. Tree measuring guidelines of the eastern native tree society. Bull. East. Nat. Tree. Soc. 2006, 1, 3–10. [Google Scholar]
  14. Curtis, R.O. Height-diameter and height-diameter-age equations for second-growth Douglas-fir. For. Sci. 1967, 13, 365–375. [Google Scholar]
  15. Wykoff, W.R.; Crookston, N.L.; Stage, A.R. User’s guide to Stand Prognosis Model, USDA Forest Service. Int. For. Rang. Exp. Stat. Gen. Tech. Rep. 1982, 1982, 112. [Google Scholar]
  16. Huang, S.; Titus, S.J.; Wiens, D.P. Comparison of nonlinear height–diameter functions for major Alberta tree species. Can. J. For. Res. 1992, 22, 1297–1304. [Google Scholar] [CrossRef]
  17. Zhang, S.; Burkhart, H.E.; Amateis, R.L. The influence of thinning on tree height and diameter relationships in loblolly pine plantations. South. J. Appl. For. 1997, 21, 199–205. [Google Scholar] [CrossRef] [Green Version]
  18. Fang, Z.; Bailey, R. Height–diameter models for tropical forests on Hainan Island in southern China. For. Ecol. Manag. 1998, 110, 315–327. [Google Scholar] [CrossRef]
  19. Peng, C. Nonlinear Height-Diameter Models for Nine Boreal Forest Tree Species in Ontario; Ontario Ministry of Natural Resources: Peterborough, ON, Canada, 1999; OFRI-Rep. 155. [Google Scholar]
  20. Jayaraman, K.; Zakrzewski, W. Practical approaches to calibrating height–diameter relationships for natural sugar maple stands in Ontario. For. Ecol. Manag. 2001, 148, 169–177. [Google Scholar] [CrossRef]
  21. Calama, R.; Cañadas, N.; Montero, G. Inter-regional variability in site index models for even-aged stands of stone pine (Pinus pinea L.) in Spain. Annal. For. Sci. 2003, 60, 259–269. [Google Scholar] [CrossRef] [Green Version]
  22. Sharma, R.P.; Breidenbach, J. Modeling height-diameter relationships for Norway spruce, Scots pine, and downy birch using Norwegian national forest inventory data. For. Sci. Technol. 2015, 11, 44–53. [Google Scholar] [CrossRef]
  23. Ciceu, A.; Garcia-Duro, J.; Seceleanu, I.; Badea, O. A generalized nonlinear mixed-effects height–diameter model for Norway spruce in mixed-uneven aged stands. For. Ecol. Manag. 2020, 477, 118507. [Google Scholar] [CrossRef]
  24. Sánchez, C.A.L.; Varela, J.G.; Dorado, F.C.; Alboreca, A.R.; Soalleiro, R.R.; González, J.G.Á.; Rodríguez, F.S. A height-diameter model for Pinus radiata D. Don in Galicia (Northwest Spain). Ann. For. Sci. 2003, 60, 237–245. [Google Scholar] [CrossRef]
  25. Temesgen, H.; van Gadow, K. Generalized height–diameter models—An application for major tree species in complex stands of interior British Columbia. Eur. J. For. Res. 2004, 123, 45–51. [Google Scholar] [CrossRef]
  26. Mugasha, W.; Mauya, E.; Njana, A.; Karlsson, K.; Malimbwi, R.; Ernest, S. Height-diameter allometry for tree species in tanzania mainland. Int. J. For. Res. 2019, 2019, 4832849. [Google Scholar] [CrossRef] [Green Version]
  27. Lang, A. Bestandeseinheitshö henkurven der Württembergischen Forsteinrichtungsanstalt. Allg. Forst Jagdztg. 1938, 114, 168–176. [Google Scholar]
  28. Kramer, H. Die genauigkeit der massenermittlung nach dem Reihenverfahren-zu dem gleichlautenden beitrag von Oberforstmeister von Laer. Forst Holzwirt. 1964, 19, 140–141. [Google Scholar]
  29. Von Laer, W. Die Genauigkeit der Massenermittlung nach dem ‘‘Reihenverfahren’’ [Accuracy of volume assessment using a sequential approach]. Forst. Holzwirt. 1964, 139–140. [Google Scholar]
  30. Kennel, E. Waldinventur und Holzaufkommensprognose für Bayern [Forest resource assessment and projections for Bavaria]. Allge. Forstzeit. 1972, 27, 615–616. [Google Scholar]
  31. Nagel, J. Einheitshö henkurvenmodell für Roteiche. Allg Forst Jagdztg. 1991, 1, 16–18. [Google Scholar]
  32. Hui, G.; von Gadow, K. Zur Entwicklung von Einheitshöhenkurven am Beispiel der Baumart Cunninghamia lanceolata. Allg. Forst. Jagdztg. 1993, 164, 218–220. [Google Scholar]
  33. Diéguez-Aranda, U.; Barrio-Anta, M.; Castedo-Dorado, F.; Álvarez-González, J. Relación altura-diámetro generalizada para masas de Pinus sylvestris L. procedentes de repoblación en el noroeste de España. Investig. Agraria. Sist. Recur. For. 2005, 14, 229–241. [Google Scholar] [CrossRef] [Green Version]
  34. Sánchez-González, M.; Cañellas, I.; Montero, G. Generalized height-diameter and crown diameter prediction models for cork oak forests in Spain. Investig. Agraria. Sist. Recur. For. 2007, 16, 76–88. [Google Scholar] [CrossRef] [Green Version]
  35. Temesgen, H.; Zhang, C.; Zhao, X. Modelling tree height–diameter relationships in multi-species and multi-layered forests: A large observational study from Northeast China. For. Ecol. Manag. 2014, 316, 78–89. [Google Scholar] [CrossRef]
  36. Liu, C.L.C.; Kuchma, O.; Krutovsky, K.V. Mixed-species versus monocultures in plantation forestry: Development, benefits, ecosystem services and perspectives for the future. Glob. Ecol. Conserv. 2018, 15, e00419. [Google Scholar] [CrossRef]
  37. Riofrío, J.; del Río, M.; Maguire, D.A.; Bravo, F. Species mixing effects on height–diameter and basal area increment models for scots pine and maritime pine. Forests 2019, 10, 249. [Google Scholar] [CrossRef] [Green Version]
  38. Xie, L.; Widagdo, F.R.A.; Dong, L.; Li, F. Modeling Height–Diameter Relationships for Mixed-Species Plantations of Fraxinus mandshurica Rupr. and Larix olgensis Henry in Northeastern China. Forests 2020, 11, 610. [Google Scholar] [CrossRef]
  39. Forrester, D.I.; Benneter, A.; Bouriaud, O. Diversity and competition influence tree allometric relationships—Developing functions for mixed-species forests. J. Ecol. 2017, 105, 761–774. [Google Scholar] [CrossRef] [Green Version]
  40. Pretzsch, H. The Effect of Tree Crown Allometry on Community Dynamics in Mixed-Species Stands versus Monocultures. A Review and Perspectives for Modeling and Silvicultural Regulation. Forests 2019, 10, 810. [Google Scholar] [CrossRef] [Green Version]
  41. Del Río, M.; Bravo-Oviedo, A.; Ruiz-Peinado, R.; Condés, S. Tree allometry variation in response to intra- and inter-specific competitions. Trees 2018, 33, 121–138. [Google Scholar] [CrossRef]
  42. Uzquiano, S.; Barbeito, I.; San Martín, R.; Ehbrecht, M.; Seidel, D.; Bravo, F. Quantifying Crown Morphology of Mixed Pine-Oak Forests Using Terrestrial Laser Scanning. Remote Sens. 2021, 13, 4955. [Google Scholar] [CrossRef]
  43. Pardos, M.; del Río, M.; Pretzsch, H.; Jactel, H.; Bielak, K.; Bravo, F.; Brazaitis, G.; Defossez, E.; Engel, M.; Godvod, K.; et al. The greater resilience of mixed forests to drought mainly depends on their composition: Analysis along a climate gradient across Europe. For. Ecol. Manag. 2021, 481, 118687. [Google Scholar] [CrossRef]
  44. Fortin, M.; Van Couwenberghe, R.; Perez, V.; Piedallu, C. Evidence of climate effects on the height-diameter relationships of tree species. Ann. For. Sci. 2019, 76, 1–20. [Google Scholar] [CrossRef] [Green Version]
  45. Condés, S.; Aguirre, A.; del Río, M. Crown plasticity of five pine species in response to competition along an aridity gradient. For. Ecol. Manag. 2020, 473, 118302. [Google Scholar] [CrossRef]
  46. Maestre, F.T.; Callaway, R.M.; Valladares, F.; Lortie, C.J. Refining the stress-gradient hypothesis for competition and facilitation in plant communities. J. Ecol. 2009, 97, 199–205. [Google Scholar] [CrossRef]
  47. Schwinning, S.; Weiner, J. Mechanisms determining the degree of size asymmetry in competition among plants. Oecologia 1998, 113, 447–455. [Google Scholar] [CrossRef] [PubMed]
  48. Herrero, C.; Bravo, F. Can we get an operational indicator of forest carbon sequestration?: A case study from two forest regions in Spain. Ecol. Indic. 2012, 17, 120–126. [Google Scholar] [CrossRef]
  49. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  50. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project. 2021. Available online: http://qgis.osgeo.org (accessed on 27 October 2021).
  51. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn Res. 2011, 12, 2825–2830. [Google Scholar]
  52. Robinson, A.P.; Froese, R.E. Model validation using equivalence tests. Ecol. Model. 2004, 176, 349–358. [Google Scholar] [CrossRef]
  53. Pinheiro, J.; Bates, D.; DebRoy, S.; Sarkar, D.; R Core Team. Nlme: Linear and Nonlinear Mixed Effects Models. R Package Version 3.1-152. Available online: https://CRAN.R-project.org/package=nlme (accessed on 27 October 2021).
  54. Vargas-Larreta, B.; Castedo-Dorado, F.; Álvarez-González, J.G.; Barrio-Anta, M.; Cruz-Cobos, F. A generalized height–diameter model with random coefficients for uneven-aged stands in El Salto, Durango (Mexico). Forestry 2009, 82, 445–462. [Google Scholar] [CrossRef] [Green Version]
  55. Crecente-Campo, F.; Corral-Rivas, J.J.; Vargas-Larreta, B.; Wehenkel, C. Can random components explain differences in the height–diameter relationship in mixed uneven-aged stands? Ann. For. Sci. 2014, 71, 51–70. [Google Scholar] [CrossRef] [Green Version]
  56. Harrison, X.A.; Donaldson, L.; Correa-Cano, M.E.; Evans, J.; Fisher, D.N.; Goodwin, C.E.; Robinson, B.S.; Hodgson, D.J.; Inger, R. A brief introduction to mixed effects modelling and multi-model inference in ecology. PeerJ 2018, 6, e4794. [Google Scholar] [CrossRef] [Green Version]
  57. Sterba, H.; del Rio, M.; Brunner, A.; Condes, S. Effect of species proportion definition on the evaluation of growth in pure vs. mixed stands. For. Syst. 2014, 4, 547–559. [Google Scholar] [CrossRef] [Green Version]
  58. Riofrío, J.; del Río, M.; Bravo, F. Mixing effects on growth efficiency in mixed pine forests. For. Int. J. For. Res. 2017, 90, 381–392. [Google Scholar] [CrossRef] [Green Version]
  59. Kearsley, E.; Moonen, P.C.; Hufkens, K.; Doetterl, S.; Lisingo, J.; Bosela, F.B.; Boeckx, P.; Beeckman, H.; Verbeeck, H. Model performance of tree height-diameter relationships in the central Congo Basin. Ann. For. Sci. 2017, 74, 1–13. [Google Scholar] [CrossRef] [Green Version]
  60. Reineke, L.H. Perfection a stand-density index for even-aged forest. J. Agric. Res. 1933, 46, 627–638. [Google Scholar]
  61. Bravo-Oviedo, A.; Condés, S.; del Río, M.; Pretzsch, H.; Ducey, M.J. Maximum stand density strongly depends on species-specific wood stability, shade and drought tolerance. For. Int. J. For. Res. 2018, 91, 459–469. [Google Scholar] [CrossRef]
  62. Aguirre, A.; Del Rio, M.; Condés, S. Intra- and inter-specific variation of the maximum size-density relationship along an aridity gradient in Iberian pinewoods. For. Ecol. Manag. 2018, 411, 90–100. [Google Scholar] [CrossRef]
  63. Rodriguez-de-Prado, D.; San Martín, R.; Bravo, F.; de Aza, C.H. Potential climatic influence on maximum stand carrying capacity for 15 Mediterranean coniferous and broadleaf species. For. Ecol. Manag. 2020, 460, 117824. [Google Scholar] [CrossRef]
  64. Cañadas, N.; García, C.; Montero, G. Relación Altura-Diámetro Para Pinus Pinea L. en el Sistema Central. Congreso de Ordenación y Gestión Sostenible de Montes, Santiago de Compostela; Secretaría de Educación Pública: Centro Histórico, Mexico, 1999; pp. 139–153. [Google Scholar]
  65. Del Rio, M. Thinning Regime and Production Model for Pinus sylvestris L. in the Iberian and Central Ranges [Spain]. Ph.D. Thesis, Centre for Forest Research, INIA, Madrid, Spain, 1999. [Google Scholar]
  66. Gaffrey, D. Forstamts-und Bestandesindividuelles Sortimentierungsprogramm als Mittel zur Plannung, Aushaltung und Simulation; Diplomarbeit Forscliche Facultät, Universität Göttingen: Göttingen, Germany, 1988. [Google Scholar]
  67. Michailoff, I. Zahlenmässiges Verfahren für die Ausführung der Bestandeshöhenkurven. Forstwiss. Cent. Tharandter Forstl. Jahrb. 1943, 6, 273–279. [Google Scholar]
  68. Monness, E.N. Diameter distributions and height curves in even-aged stands of Pinus sylvestris L. Medd. Nor. Inst. Skogforsk. 1982, 36, 1–43. [Google Scholar]
  69. Nilson, A. Pidev metsakorraldus-mis see on. Pidev metsakorraldus. EPMÜ Metsandusteaduskonna Toim. 1999, 32, 4–13. [Google Scholar]
  70. Meyer, H.A. A mathematical expression for height curves. J. Forest. 1940, 38, 415–420. [Google Scholar]
  71. Pearl, R.; Reed, L.J. On the rate of growth of the population of the United States since 1790 and its mathematical representation. Proc. Natl. Acad. Sci. USA 1920, 6, 275–288. [Google Scholar] [CrossRef] [Green Version]
  72. Ratkowsky, D.A.; Reedy, T.J. Choosing near-linear parameters in the fourparameter logistic model for radioligand and related assays. Biometrics 1986, 42, 575–582. [Google Scholar] [CrossRef]
  73. Richards, F.J. A flexible growth function for empirical use. J. Exp. Biol. 1959, 10, 290–300. [Google Scholar] [CrossRef]
  74. Schumacher, F.X. A new growth curve and its application to timber yield studies. J. For. 1939, 37, 819–820. [Google Scholar]
  75. Seber, G.A.F.; Wild, C.J. Nonlinear Regression; Wiley: New York, NY, USA, 1989; 768p. [Google Scholar]
  76. Zeide, B. Analysis of growth equations. For. Sci. 1992, 39, 594–616. [Google Scholar] [CrossRef]
  77. Del Río, M.; Pretzsch, H.; Alberdi, I.; Bielak, K.; Bravo, F.; Brunner, A.; Condés, S.; Ducey, M.J.; Fonseca, T.; von Lüpke, N. Characterization of the structure, dynamics, and productivity of mixed-species stands: Review and perspectives. Eur. J. For. Res. 2016, 135, 23–49. [Google Scholar] [CrossRef]
  78. Pretzsch, H.; Río, M.; Schütze, G.; Ammer, C.; Annighöfer, P.; Avdagić, A.; Barbeito, I.; Bielak, K.; Brazaitis, G.; Coll, L.; et al. Mixing of Scots pine (Pinus sylvestris L.) and European beech (Fagus sylvatica L.) enhances structural heterogeneity, and the effect increases with water availability. For. Ecol. Manag. 2016, 373, 149–166. [Google Scholar] [CrossRef] [Green Version]
  79. Riofrío, J.; del Río, M.; Pretzsch, H.; Bravo, F. Changes in structural heterogeneity and stand productivity by mixing Scots pine and Maritime pine. For. Ecol. Manag. 2017, 405, 219–228. [Google Scholar] [CrossRef]
  80. De Martonne, E. L’indice d’aridité. Bull. Assoc. Géographes Français 1926, 3, 3–5. [Google Scholar] [CrossRef]
  81. Bravo, F.; Álvarez-González, J.G.; del Río, M.; Barrio-Anta, M.; Bravo-Oviedo, A.; Calama, R.; Castedo-Dorado, F.; Crecente-Campo, F.; Condés, S. Growth and Yield Models in Spain: Historical Overview, Contemporary Examples and Perspectives; Instituto Universitario de Investigación en Gestión Forestal Sostenible: Madrid, Spain, 2012. [Google Scholar]
  82. Mehtätalo, L.; de-Miguel, S.; Timothy, G. Gregoire. Modeling height-diameter curves for prediction. Can. J. For. Res. 2015, 45, 826–837. [Google Scholar] [CrossRef] [Green Version]
  83. Valladares, F.; Niinemets, Ü. Shade Tolerance, a Key Plant Feature of Complex Nature and Consequences. Annu. Rev. Ecol. Evol. 2008, 39, 237–257. [Google Scholar] [CrossRef] [Green Version]
  84. Franceschini, T.; Schneider, R. Influence of shade tolerance and development stage on the allometry of ten temperate tree species. Oecologia 2014, 176, 739–749. [Google Scholar] [CrossRef]
  85. Martín-Benito, D.; Cherubini, P.; del Río, M.; Cañellas, I. Growth response to climate and drought in Pinus nigra Arn. trees of different crown classes. Trees 2008, 22, 363–373. [Google Scholar] [CrossRef]
  86. Savi, T.; Casolo, V.; Dal-Borgo, A.; Rosner, S.; Torboli, V.; Stenni, B.; Bertoncin, P.; Martellos, S.; Pallavicini, A.; Nardini, A. Drought-induced dieback of Pinus nigra: A tale of hydraulic failure and carbon starvation. Conserv. Physiol. 2020, 7, coz012. [Google Scholar] [CrossRef]
  87. Pretzsch, H.; Bielak, K.; Block, J.; Bruchwald, A.; Dieler, J.; Ehrhart, H.P.; Kohnle, U.; Nagel, J.; Spellmann, H.; Zasada, M.; et al. Productivity of mixed versus pure stands of oak (Quercus petraea (M att.) Liebl. and Quercus robur L.) and European beech (Fagus sylvatica L.) along an ecological gradient. Eur. J. For. Res. 2013, 132, 263–280. [Google Scholar] [CrossRef]
  88. Ligot, G.; Balandier, P.; Fayolle, A.; Lejeune, P.; Claessens, H. Height competition between Quercus petraea and Fagus sylvatica natural regeneration in mixed and uneven-aged stands. For. Ecol. Manag. 2013, 304, 391–398. [Google Scholar] [CrossRef]
  89. Pretzsch, H.; Block, J.; Dieler, J.; Dong, P.H.; Kohnle, U.; Nagel, J.; Spellmann, H.; Zingg, A. Comparison between the productivity of pure and mixed stands of Norway spruce and European beech along an ecological gradient. Annal. For. Sci. 2010, 67, 712. [Google Scholar] [CrossRef] [Green Version]
  90. Pausas, J.G.; Bladé, C.; Valdecantos, A.; Seva, J.P.; Fuentes, D.; Alloza, J.A.; Vilagrosa, A.; Bautista, S.; Cortina, J.; Vallejo, R. Pines and oaks in the restoration of Mediterranean landscapes of Spain: New perspectives for an old practice—A review. Plant Ecol. 2004, 171, 209–220. [Google Scholar] [CrossRef]
  91. Vicente, E.; Vilagrosa, A.; Ruiz-Yanetti, S.; Manrique-Alba, À.; González-Sanchís, M.; Moutahir, H.; Chirino, E.; Del Campo, A.; Bellot, J. Water Balance of Mediterranean Quercus ilex L. and Pinus halepensis Mill. Forests in Semiarid Climates: A Review in A Climate Change Context. Forests 2018, 9, 426. [Google Scholar] [CrossRef] [Green Version]
  92. Pan, Y.; Birdsey, R.A.; Fang, J.; Houghton, R.; Kauppi, P.E.; Kurz, W.A.; Phillips, O.L.; Shvidenko, A.; Lewis, S.L.; Canadell, J.G. A large and persistent carbon sink in the world’s forests. Science 2011, 333, 988–993. [Google Scholar] [CrossRef] [Green Version]
  93. Ng’andwe, P.; Chungu, D.; Yambayamba, A.M.; Chilambwe, A. Modeling the height-diameter relationship of planted Pinus kesiya in Zambia. For. Ecol. Manag. 2019, 447, 1–11. [Google Scholar] [CrossRef]
  94. Ruiz-Peinado, R.; Pretzsch, H.; Löf, M.; Heym, M.; Bielak, K.; Aldea, J.; Barbeito, I.; Brazaitis, G.; Drössler, L.; Godvod, K. Mixing effects on Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies (L.) Karst.) productivity along a climatic gradient across Europe. For. Ecol. Manag. 2021, 482, 118834. [Google Scholar] [CrossRef]
  95. Özçelik, R.; Yavuz, H.; Karatepe, Y.; Gürlevik, N.; Kiriş; R. Development of ecoregion-based height-diameter models for 3 economically important tree species of southern Turkey. Turk. J. Agric. For. 2014, 38, 399–412. [Google Scholar] [CrossRef]
  96. Grossiord, C. Having the right neighbors: How tree species diversity modulates drought impacts on forests. New Phytol. 2020, 228, 42–49. [Google Scholar] [CrossRef] [Green Version]
  97. Martínez del Castillo, E.; Tejedor, E.; Serrano-Notivoli, R.; Novak, K.; Saz, M.Á.; Longares, L.A.; De Luis, M. Contrasting Patterns of Tree Growth of Mediterranean Pine Species in the Iberian Peninsula. Forests 2018, 9, 416. [Google Scholar] [CrossRef] [Green Version]
  98. Vergarechea, M.; Calama, R.; Pretzsch, H.; Alday, J.; Río, M. Short-and long-term growth response to climate in mixed and monospecific forests of Pinus pinea and Pinus pinaster. Eur. J. For. Res. 2021, 140, 387–402. [Google Scholar] [CrossRef]
  99. Grossiord, C.; Granier, A.; Ratcliffe, S.; Bouriaud, O.; Bruelheide, H.; Chećko, E.; Forrester, D.I.; Dawud, S.H.; Finér, L.; Pollastrini, M. Tree diversity does not always improve resistance of forest ecosystems to drought. PNAS 2014, 111, 14812–14815. [Google Scholar] [CrossRef] [Green Version]
  100. Fernández-de-Uña, L.; Rossi, S.; Aranda, I.; Fonti, P.; González-González, B.D.; Cañellas, I.; Gea-Izquierdo, G. Xylem and Leaf Functional Adjustments to Drought in Pinus sylvestris and Quercus pyrenaica at Their Elevational Boundary. Front. Plant Sci. 2017, 8, 1200. [Google Scholar] [CrossRef] [PubMed]
  101. Toïgo, M.; Vallet, P.; Tuilleras, V.; Lebourgeois, F.; Rozenberg, P.; Perret, S.; Courbaud, B.; Perot, M. Species mixture increases the effect of drought on tree ring density, but not on ring width, in Quercus petraeaPinus sylvestris stands. For. Ecol. Manag. 2015, 345, 73–82. [Google Scholar] [CrossRef]
Figure 1. Variation of height under different species’ mixing proportions (mix-axis) based on the M(ratio) (y-axis). Letters a-m indicate the species name in order to make reference to them along the text. Note: M(ratio) represents whether height in mixed forests is higher (M(ratio) > 1) or lower (M(ratio) < 1) than pure stands. All independent variables (excepting mi) have been fixed to the mean of each species.
Figure 1. Variation of height under different species’ mixing proportions (mix-axis) based on the M(ratio) (y-axis). Letters a-m indicate the species name in order to make reference to them along the text. Note: M(ratio) represents whether height in mixed forests is higher (M(ratio) > 1) or lower (M(ratio) < 1) than pure stands. All independent variables (excepting mi) have been fixed to the mean of each species.
Forests 13 00119 g001
Figure 2. Variation of height along the aridity gradient (x-axis), with Arid being the lowest M value and Humid the highest M found in the distribution of each species in each mixture composition. Letters a-m indicate the species name in order to make reference to them along the text. Note: C(ratio) ratio (y-axis) represents whether height in mixed forests is higher (C(ratio) > 1) or lower (C(ratio) < 1) in more arid than humid places. All independent variables (excepting M and mi = 0. 5) have been fixed to the mean values for each species.
Figure 2. Variation of height along the aridity gradient (x-axis), with Arid being the lowest M value and Humid the highest M found in the distribution of each species in each mixture composition. Letters a-m indicate the species name in order to make reference to them along the text. Note: C(ratio) ratio (y-axis) represents whether height in mixed forests is higher (C(ratio) > 1) or lower (C(ratio) < 1) in more arid than humid places. All independent variables (excepting M and mi = 0. 5) have been fixed to the mean values for each species.
Forests 13 00119 g002
Figure 3. Aridity influence in terms of the De Martonne Index (M), on H–D relationship for different stand dominant height (Ho) in (a) Pinus halepensis mixed with Quercus ilex and (b) Quercus robur mixed with Fagus sylvatica, depending on parameter significance from selected H–D models. Diameter at breast height (d, in centimeters) is represented in x-axis, since total tree height (h, in meters) is represented in y-axis.
Figure 3. Aridity influence in terms of the De Martonne Index (M), on H–D relationship for different stand dominant height (Ho) in (a) Pinus halepensis mixed with Quercus ilex and (b) Quercus robur mixed with Fagus sylvatica, depending on parameter significance from selected H–D models. Diameter at breast height (d, in centimeters) is represented in x-axis, since total tree height (h, in meters) is represented in y-axis.
Forests 13 00119 g003
Table 1. Base tree Height–Diameter equations used in this study. h = total tree height; d: diameter at breast height; Ho = dominant height of the plot; do = dominant diameter of the plot, β0, β1, and β2: model parameters.
Table 1. Base tree Height–Diameter equations used in this study. h = total tree height; d: diameter at breast height; Ho = dominant height of the plot; do = dominant diameter of the plot, β0, β1, and β2: model parameters.
ModelReferenceFormula
M1[64] h = 1.3 + ( β 0 ( 1 d 1 do ) + ( 1 Ho 1.3 ) 1 / 2 ) 2
M2[65] h = β 0 e β 1 ( 1 d 1 do )
M3[66] h = 1.3 + ( Ho 1.3 ) e β 0 ( 1 do d ) + β 1 ( 1 do 1 d )
M4[67] h = Ho e β 0 ( 1 d 1 do )
M5[68] h = 1.3 + ( β 0 ( 1 d 1 do ) + ( 1 Ho 1.3 ) 1 / 3 ) 3
M6[69] h = Ho 1 β 0 ( 1 ( do d ) β 1 )
M7[16] h = β 0 e β 1 e ( β 2 d )
M8[70] h = β 0 ( 1 e ( β 1 d ) )
M9[71] h = β 0 ( 1 + β 1 e ( β 2 d ) )
M10[72] h = β 0 ( 1 + β 1 1 d β 2 )
M11[73] h = β 0 ( 1 e β 1 d ) β 2
M12[74] h = 1.3 + β 0 e β 1 d
M13[75] h = β 0 e e ( β 1 ( d β 2 ) )
M14[76] h = β 0 e ( β 1 e ( β 2 d β 3 ) )
Table 2. Estimated model coefficients for the selected H–D models for the different species and species composition.
Table 2. Estimated model coefficients for the selected H–D models for the different species and species composition.
SpeciesMixtureModelForm *α0 (Int)α1 (BAL)α2 (Dq)α2 (Dqi)α3 (mi)α4 (M)α5 (Ho)α5 (Hoi)β0β1RMSEAIC
Pinus halepensisPinus nigraM3A0.8740−0.0097 −0.00920.1385 −5.62541.07625950.45
Pinus pinasterM1A1.2130−0.0402 0.0431 1.20086675.29
Pinus pineaM1A −0.0125 0.03901.2598 1.06235022.50
Quercus fagineaM1A 0.02941.0935 1.04672319.73
Quercus ilexM1A1.3583−0.0151 −0.02001.4876−0.0307 1.119512,939.80
Pinus nigraPinus halepensisM5A0.9181−0.01680.0797 1.12186870.48
Pinus pinasterM1A1.0318−0.01120.0440 0.6216 1.24791,265,284
Pinus sylvestrisM1A1.4454−0.01370.0474 1.251035,767.45
Quercus fagineaM1A0.7447 −0.03192.1911 1.25868690.03
Quercus ilexM1A −0.0068 −0.02092.8952 1.374519,658.38
Pinus pinasterPinus halepensisM1A −0.01730.0734 1.5235 1.22395962.40
Pinus nigraM1A1.8800−0.01830.0499 1.076112,580.34
Pinus pineaM1M 0.23140.3722 1.046410,159.86
Pinus sylvestrisM1A2.7801−0.0132−0.0203 1.211921,720.17
Quercus ilexM1M 2.82660.3382 1.18135850.48
Quercus pyrenaicaM1A −0.0104 0.03191.6618 1.07018513.35
Quercus suberM3A−0.6969 0.01370.1964 13.94931.07514130.34
Pinus pineaPinus halepensisM5M 0.0997 0.81001898.24
Pinus pinasterM1A4.0383−0.0221 1.5266−0.1198 0.98006111.83
Quercus ilexM3A1.3421 −0.0162 −0.0079 −8.42470.88003512.50
Quercus suberM1A2.1130−0.0258 −0.04601.2023 0.93001979.70
Pinus sylvestrisFagus sylvaticaM3A−0.1063 0.0037 7.19871.610033,678.46
Pinus nigraM5A0.8792−0.00680.0465 1.190042,508.63
Pinus pinasterM3A0.1674 −0.0846 5.78481.340029,369.08
Pinus uncinataM5A 0.0271 1.620022,112.05
Quercus fagineaM5M 0.12860.6005 1.250011,709.53
Quercus ilexM1A−1.1988−0.0104 0.02411.63560.0178 1.330014,518.65
Quercus petraeaM5A 1.3877 1.53005481.21
Quercus pyrenaicaM3A 0.0038 4.36521.400033,274.37
Pinus uncinataPinus sylvestrisM5A 0.0244 1.520016,426.61
Fagus sylvaticaPinus sylvestrisM5A −0.0049 0.02430.6819 2.080011,466.10
Quercus petraeaM5A0.9408−0.00660.0184 1.93009802.89
Quercus pyrenaicaM1M 0.32930.5965−0.2352 2.10009292.34
Quercus roburM12M −0.1038 −0.05090.12441.0300 −9.48272.46005726.19
Quercus fagineaPinus halepensisM12M −0.2481 1.3855 −4.27071.0900355.92
Pinus nigraM12M −0.3677 0.2853 1.1564 −7.01281.05001709.34
Pinus sylvestrisM12M −0.3601 0.3284 1.0817 −8.42751.34003439.68
Quercus ilexM1A0.7578 0.8303 1.24004311.78
Quercus ilexPinus halepensisM12M −0.2772−0.04250.1651 1.1925 −6.07630.96005638.96
Pinus nigraM12M −0.3333 0.1989 1.2466 −6.51640.85004446.14
Pinus pinasterM12M −0.1802 0.1547 1.0835 −5.97810.89002131.00
Pinus pineaM12M −0.1910 0.1413 1.1488 −6.69680.95002836.90
Pinus sylvestrisM12M −0.3877−0.01430.2363 1.2026 −6.39410.98005086.48
Quercus fagineaM1A −0.0282 1.16360.0407 0.87004927.66
Quercus ilexQuercus pyrenaicaM12M −0.2510 1.4862 −6.86550.9300838.59
Quercus suberM12A −0.2571−0.03240.2214 1.0963 −7.25681.07005404.67
Quercus petraeaPinus sylvestrisM12M −0.3962 0.3045 1.1671 −9.13251.4800986.33
Fagus sylvaticaM12M −0.2111−0.03640.1916 1.0490 −7.33591.69003202.77
Quercus pyrenaicaFagus sylvaticaM1M 0.1170 1.51009281.96
Pinus pinasterM1A −0.03680.1155 1.26001811.55
Pinus sylvestrisM1A −0.0091 0.01631.2693 1.39008372.09
Quercus ilexM5A 1.4330 1.32001653.03
Quercus roburM12M −0.2799 0.2404 1.0418 −7.70231.8000883.50
Quercus roburFagus sylvaticaM12M −0.2329 0.3245 0.9650 −11.73812.58001000.90
Quercus pyrenaicaM12M −0.2167 0.1838 1.0630 −8.18761.98002032.19
Quercus suberPinus pinasterM1A −0.0332 2.65350.0252 1.09002779.01
Pinus pineaM1M 0.36200.8264 0.99003862.59
Quercus ilexM1A1.4835−0.03860.0578 0.8232 1.12008671.20
* Form refers to the way whether the asymptotic parameter (β0) is expanded in an Additive (A) or Multiplicative (M) way. See Section 2.2.2 for further explanation. Note: Only coefficients of significative (p-value < 0.05) variables are presented for each model.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rodríguez de Prado, D.; Riofrío, J.; Aldea, J.; McDermott, J.; Bravo, F.; Herrero de Aza, C. Species Mixing Proportion and Aridity Influence in the Height–Diameter Relationship for Different Species Mixtures in Mediterranean Forests. Forests 2022, 13, 119. https://0-doi-org.brum.beds.ac.uk/10.3390/f13010119

AMA Style

Rodríguez de Prado D, Riofrío J, Aldea J, McDermott J, Bravo F, Herrero de Aza C. Species Mixing Proportion and Aridity Influence in the Height–Diameter Relationship for Different Species Mixtures in Mediterranean Forests. Forests. 2022; 13(1):119. https://0-doi-org.brum.beds.ac.uk/10.3390/f13010119

Chicago/Turabian Style

Rodríguez de Prado, Diego, Jose Riofrío, Jorge Aldea, James McDermott, Felipe Bravo, and Celia Herrero de Aza. 2022. "Species Mixing Proportion and Aridity Influence in the Height–Diameter Relationship for Different Species Mixtures in Mediterranean Forests" Forests 13, no. 1: 119. https://0-doi-org.brum.beds.ac.uk/10.3390/f13010119

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop