Individual Tree Diameter and Height Growth Models for 30 Tree Species in Mixed-Species and Uneven-Aged Forests of Mexico

: Lack of knowledge of individual tree growth in species-rich, mixed forest ecosystems impedes their sustainable management. In this study, species-speciﬁc models for predicting individual diameter at breast height (dbh) and total tree height (h) growth were developed for 30 tree species growing in mixed and uneven-aged forest stands in Durango, Mexico. Growth models were also developed for all pine, all oaks, and all other species of the genus Arbutus (strawberry trees). A database of 55,158 trees with remeasurements of dbh and h of a 5-year growth period was used to develop the models. The data were collected from 217 stem-mapped plots located in the Sierra Madre Occidental (Mexico). Weighted regression was used to remove heteroscedasticity from the species-speciﬁc dbh and h growth models using a power function of the tree size independent variables. The ﬁnal models developed in the present study to predict dbh and total tree height growth included size variables, site factors, and competition variables in their formulation. The developed models ﬁtted the data well and explained between 98 and 99% and of the observed variation of dbh, and between 77 and 98% of the observed variation of total tree height for the studied species and groups of species. The developed models can be used for estimating the individual dbh and h growth for the analyzed species and can be integrated in decision support tools for management planning in these mixed forest ecosystems.


Introduction
Accurate estimates of both current resource levels and the expected resource changes from implementing selected management alternatives are needed to develop tools to support forest management informed decision-making (e.g., [1]). Forest growth models describe the change in size of a tree or another tree attribute or stand-level variable with age. They are used to analyze and estimate different relationships in forest stand development such as species composition, site characteristics, tree competitive status, and silvicultural management (e.g., [2][3][4]). Such models, together with stand regeneration, harvesting, and mortality models, are important tools in long-term forest management systems because they can simulate stand development and production under various management alternatives (e.g., [3,[5][6][7]).
Empirical forest tree growth models predict growth as a function of tree attributes (such as tree size, crown ratio, vigor and tree competition) and plot or stand level variables (such as age site index and stand density) [8][9][10][11][12]. In particular, individual-tree diameter and height growth models are a fundamental component of forest growth and yield frameworks and constitute very important tools for sustainable forest management planning [5,7]. Although many factors such as tree size, neighborhood competition and environmental variables, have been proposed as being important in explaining patterns of tree growth, their relative contribution for different species and ecosystems are still subject to debate (e.g., [5,[13][14][15]). Improving our knowledge of the controls influencing tree growth is central to better understand forest community organization and dynamics [15], and is also important to the development of sustainable management of forest ecosystems [16].
Understanding the factors that should be included into growth models to explain patterns in mixed-species forest with high species richness is critical for modeling ecological dynamics in such stands [16]. In uneven-aged, mixed-species natural forests, tree-to-tree variation in growth can be considerable within and among species, as well as spatially within a stand (e.g., [17]). This variability results from the structural and compositional complexity that often characterize these forests (e.g., [13,[17][18][19]). In this context, an improved understanding of the factors regulating variability in growth rates can lead to better anticipation of individual tree success and community dynamics [17]. Such understanding can help in developing biometric tools to support forest management decision-making in mixed-species stands, including the development of management strategies based on natural stand dynamics (e.g., [18,20,21]). Nevertheless, factors influencing tree growth in structurally complex forests remain poorly understood [18]. Most previous studies of competition and tree growth have focused on plantations, managed forests, and even-aged, relatively young natural forests (e.g., [22,23]). In contrast, studies analyzing individual tree growth in mixed forest ecosystems are relatively more scarce (e.g., [13,17,18,24]), but have received increased attention in the last decades due to a worldwide trend in managing such forests to improve the biodiversity and ecosystems services and to assure long-term sustainable forest resources [3,25].
The Sierra Madre Occidental in the Mexican State of Durango constitutes a hotspot of forest biodiversity, being home to about five million hectares of species-rich and uneven-aged forest stands [26]. It is the first forest reserve in Mexico with an estimate of 422.1 million m 3 of timber volume, and represents 22.4% of the total estimated production at the national level [27,28]. In 2015, Durango produced 28.5% of the total timber volume harvested in Mexico [3]. These forests are dominated by irregular and multispecies stands [29], and have been subject to selective harvesting for more than a century to provide a mix of services to local communities [30]. Although the main objective of forestry in Durango has been timber production, increased emphasis is being given to the concepts of ecological sustainability and biodiversity. Consequently, the development of forest growth models for mixed and uneven-aged stands is imperative to be used as important tools to support the informed decision-making in this type of forest ecosystems. To date in Durango, forest growth studies have been conducted on site quality [31][32][33], stand density management [28], dbh growth models [3], tree competition [34,35], tree volume and biomass equations [27,36,37] and tree diameter-height relationships [38]. However, most of the previous studies dealt with only one or a few tree species and there have been very limited efforts to develop individual-tree dbh or height growth models for most of the commercial tree species in mixed-species forests. Definitely, we are not aware of studies that have developed individual-tree height growth models for the main tree species in Durango using data of stem-mapped and re-measured plots. In this situation, more research and studies become necessary to better understand stand growth and predict the effect of silvicultural practices in these species-rich and uneven-aged forests that are among the richest in the country. Therefore, the objective of this study was to develop new species-specific models for predicting individual-tree dbh and total tree height growth for 30 commercial species growing in the mixed and uneven-aged stands in Durango, Mexico. Some generic growth equations were also developed for all pines, all oaks, and all strawberry species. Such generic group equations can provide useful estimates for a few remaining tree species with still not enough remeasurements of dbh and total tree height in the area of study.

Study Area and Data
The data used in this study were collected in the Sierra Madre Occidental in the State of Durango, Northwestern Mexico, located between the geographic coordinates 22 • 20 and 26 • 46 N latitude and from 103 • 46 to 107 • 11 W longitude. The data are part of a network of stem-mapped re-measured plots, which were established to monitor the stand dynamics of Durango's forests [38]. Plot establishment began in 2007 and re-measurements started in 2012, after a five-year interval. Each plot covers an area of 50 × 50 m. The plots are generally distributed systematically, with a variable grid ranging from 3 to 5 km (Figure 1). The tree number, genus, species, dbh outside bark (dbh; cm), total tree height (h; m), and other variables of all trees with dbh ≥ 7.5 cm were recorded. The installation and re-measurement of the plots followed the methodology developed by Corral-Rivas et al. [39]. Only undisturbed plots without intermediate harvesting between measurements were selected and considered for this study. A total of 217 re-measured plots were available for the development of the individual-tree dbh and total tree height growth models. The research plots included a variety of stand site conditions, tree sizes and competition levels. A description of the main plot-level variables for the re-measured plots is summarized in Table 1. of studies that have developed individual-tree height growth models for the main tree species in Durango using data of stem-mapped and re-measured plots. In this situation, more research and studies become necessary to better understand stand growth and predict the effect of silvicultural practices in these species-rich and uneven-aged forests that are among the richest in the country. Therefore, the objective of this study was to develop new species-specific models for predicting individual-tree dbh and total tree height growth for 30 commercial species growing in the mixed and uneven-aged stands in Durango, Mexico. Some generic growth equations were also developed for all pines, all oaks, and all strawberry species. Such generic group equations can provide useful estimates for a few remaining tree species with still not enough remeasurements of dbh and total tree height in the area of study.

Study Area and Data
The data used in this study were collected in the Sierra Madre Occidental in the State of Durango, Northwestern Mexico, located between the geographic coordinates 22°20′ and 26°46′ N latitude and from 103°46′ to 107°11′ W longitude. The data are part of a network of stem-mapped re-measured plots, which were established to monitor the stand dynamics of Durango's forests [38]. Plot establishment began in 2007 and re-measurements started in 2012, after a five-year interval. Each plot covers an area of 50 × 50 m. The plots are generally distributed systematically, with a variable grid ranging from 3 to 5 km ( Figure 1). The tree number, genus, species, dbh outside bark (dbh; cm), total tree height (h; m), and other variables of all trees with dbh ≥ 7.5 cm were recorded. The installation and re-measurement of the plots followed the methodology developed by Corral-Rivas et al. [39]. Only undisturbed plots without intermediate harvesting between measurements were selected and considered for this study. A total of 217 re-measured plots were available for the development of the individual-tree dbh and total tree height growth models. The research plots included a variety of stand site conditions, tree sizes and competition levels. A description of the main plot-level variables for the re-measured plots is summarized in Table 1.    The re-measured dbh (dbh 2 ) and total tree height (h 2 ) after a 5-year period were obtained from 55,158 sampled trees of 30 species (each species included more than 100 trees for both variables

Procedures for Developing the Species-Specific Individual-Tree Growth Equations
We used the following empirical model form as starting point of the model selection (1). The equation can be also considered a form of the growth equation proposed by Wykoff [40] for mixed conifer stands in the Northern Rocky Mountains. It predicts the observed dbh and total tree height for a growth period (5-year in our study) as a function of tree initial size (dbh or height), stand level variables and tree competition as predictor variables.
where y is the dependent variable (dbh -in cm-or total tree height -in m-) re-measured in 2012; T are initial tree size variables (dbh or total tree height) measured in 2007; S represents the effect of stand variables that commonly include stand density and dominant height; C represents the effect of tree competition; α, β, γ and δ are the regression coefficients of the equation; and ε is the model error term. Initially, a regression procedure was used to select the definitive tree and stand variables for the individual tree dbh and height growth models. The significance level of 5% was used for the selection of candidate predictive variables in the model. These predictive variables were the measurements at the beginning of the selected 5-year growth period. They included initial tree dbh (dbh 1 , cm); initial tree total height (h 1 , m); the sum of the basal area of the trees larger than the subject tree (BAL, proposed by Wykoff et al. [41]), a modified BAL (BALMOD, proposed by Schröder and Gadow [42]) which includes the relative spacing index (i.e., BALMOD = (1 − (1 − (BAL/G))/IH, where IH = 100/H 0 N 1/2 ); the stand quadratic mean diameter (QMD, cm); stand density (N, number of trees per hectare); the total basal area of the plot (G, m 2 ha −1 ); stand dominant height (H 0 , m); and crown competition factor (CCF) [43]. H 0 was calculated from the proportion of the 100 thickest trees of pines per hectare. Data from trees with broken or dead tops were not used to calculate H 0. These predictive variables, along with their various combinations were tested separately by each species or species group using the lm and nls functions of R software environment [44]. In order to avoid problems related to multicollinearity, the value of the variance inflation factor (VIF) was computed for the groups of variables that best predicted the dependent variables. Candidate variables with a VIF larger than 10 were excluded, as suggested by Myers [45]. Non-significant parameters at 5% were also removed in the modeling process (i.e., they were fixed at 1 in the final model fitting). Therefore, a little bit different models were used for such species, resulting in the following two dbh and height models as the best candidate model formulations for describing both dbh and height growth in general among the models tested: Forests 2020, 11, 429 6 of 16 where dbh 2 is the diameter at breast height after the selected 5-year growth period (cm), dbh 1 is the initial tree diameter (cm), H 0 is the stand dominant height (m), BAL is the sum of the basal area of the trees larger than the subject tree proposed by Wykoff et al. [41], h 2 is the total tree height after a 5-year growth period (m), h 1 is the initial total tree height (m), a 0 -b 3 are the regression coefficients, and ε is the model error term. The models were first fitted under the assumption that the errors within equation were independent and identically distributed [46]. As only two measurements were available for each tree in the dataset, autocorrelation among measurements within a tree was not considered in the present study. However, residual analysis revealed the presence of heteroscedasticity in most species-specific dbh and total tree height growth equations. To deal with this problem weighted regression was used, modelling the error variances as a power function of the independent size variables (e.g., [46,47]). For the diameter growth model, the error variance was modelled with a power function of dbh 1 as independent variable (i.e., σ 2 i = dbh k 1 ). Using the method suggested by Park [48], squared residuals (ε 2 i ) were regressed against dbh 1 to obtain an estimate of k as follows:ε 2 i =γxdbh k 1 , this expression was linearized by the use of the natural logarithm: lnε 2 i = lnγ + kln dbh 1 and parameters were then estimated by linear least squares regression; the species-specific k value was subsequently included in the weighting factor (w i = 1/dbh k i1 ). The same was done for the species-specific height growth models (3), but using h 1 as the independent variable.

Model Comparison
Evaluation of the performance of the dbh and height growth equations was based on numerical and graphical analysis of the residuals. Two goodness-of-fit statistics were calculated: the root mean square error (RMSE) and the coefficient of determination (R 2 ). These statistics were computed on the basis of the unweighted residuals because these are the values of most interest by users of the equations. The expressions of these statistics are: where y i ,ŷ i and y are the observed, predicted, and mean values of the dependent variable, respectively; and n and p are the total number of observations and the number of parameters used to fit the models, respectively. Plots of studentized residuals against predicted values of the dependent variables were used for visual evaluation of the models [49,50]. These graphs are very important for indicating areas for which the equations provide especially poor or good predictions.

Diameter Growth Equations
The parameter estimates and the goodness-of-fit statistics for predicting dbh after a 5-year growth period with model 2 for each analyzed tree species are presented in Table 3. Most of the parameter estimates were significant at α = 0.05. All non-significant parameters were dropped from model 2. The initial tree diameter was the most significant and important predictor of dbh 2 . In all cases the coefficients a 0 and a 1 were positive, and a 1 was less than 1.0 in all cases, indicating that dbh 1 limits diameter growth by itself. The term H 0 was positive in all species except for Pinus engelmannii and Quercus arizonica (see a 2 in Table 3), resulting in increasing diameter with increasing stand dominant height. On the other hand, coefficient a 3 multiplying the BAL competition index was negative in all species except for Pinus engelmannii and P. herrerae, indicating a decrease in diameter growth as individual tree competition increases. The graphs of studentized residuals against predicted dbh values showed no evidence of heterogeneous variance over the full range of predicted values for all species after fitting the models with weighted regression. In accordance with Huang et al. [50], plots of RMSE were obtained across tree sizes to evaluate the accuracy and precision of diameter predictions of model 2 for the developed equations. Figure 2 shows the plots for the three generic equations developed in this paper (species groups: all pines, all oaks, and all strawberries). It can be inferred that the final equations performed quite well for all sampled diameter classes. All strawberries model produced larger prediction errors for large-sized trees in comparison with all oaks and all pines models. n = number of trees used in model fitting, a1-a3 = parameter estimates, ns = these parameters were not significant at p < 0.05, and were therefore fixed at 1 in the final model fitting, resulting in a little different model for these species, RMSE = root mean square error, R 2 = coefficient of determination.
The graphs of studentized residuals against predicted dbh values showed no evidence of heterogeneous variance over the full range of predicted values for all species after fitting the models with weighted regression.
In accordance with Huang et al. [50], plots of RMSE were obtained across tree sizes to evaluate the accuracy and precision of diameter predictions of model 2 for the developed equations. Figure 2 shows the plots for the three generic equations developed in this paper (species groups: all pines, all oaks, and all strawberries). It can be inferred that the final equations performed quite well for all sampled diameter classes. All strawberries model produced larger prediction errors for large-sized trees in comparison with all oaks and all pines models.

Height Growth Equations
The parameter estimates and the goodness-of-fit statistics for predicting total tree height after a 5-year growth period with model 3 for each species are presented in Table 4. Most of the parameter estimates were significant at α = 0.05. All non-significant parameters were removed from model 3. The initial total tree height was the most significant and important predictor of height after a 5-year growth period. In all cases, the coefficients b0 and b1 were positive, b1 was less than 1.0 in all cases (except for Quercus urbanii), indicating than h1 limits height growth by itself. The term H0 was positive in all species except for Quercus urbanii and Cupressus lusitanica (see b2 in Table 3), resulting in

Height Growth Equations
The parameter estimates and the goodness-of-fit statistics for predicting total tree height after a 5-year growth period with model 3 for each species are presented in Table 4. Most of the parameter estimates were significant at α = 0.05. All non-significant parameters were removed from model 3. The initial total tree height was the most significant and important predictor of height after a 5-year growth period. In all cases, the coefficients b 0 and b 1 were positive, b 1 was less than 1.0 in all cases (except for Quercus urbanii), indicating than h 1 limits height growth by itself. The term H 0 was positive in all species except for Quercus urbanii and Cupressus lusitanica (see b 2 in Table 3), resulting in increasing height growth with increasing stand dominant height. On the other hand, coefficient a 3 multiplying the BAL competition index was negative in all species except for Cupressus lusitanica, indicating a decrease in height growth as individual tree competition increases.
The developed equations fitted the data well and generally explained between 77 and 98% of the observed variation of the total tree height for all studied species. Average RMSE values of total tree height of pines, other conifers, oaks, and other broadleaf species were 1.  n = number of trees used in model fitting, b 0 -b 3 = parameter estimates, ns = these parameters were not significant at p < 0.05, and were therefore fixed at 1 in the final model fitting, resulting in a little different model for these species, RMSE = root mean square error, R 2 = coefficient of determination.
The graphs of studentized residuals against predicted total height values showed no evidence of heterogeneous variance over the full range of predicted values for all species after fitting the models with weighted regression. Figure 3 shows the plots of RMSE obtained across tree height classes for the three generic equations developed in this paper to predict total tree height after a 5-year growth period with model 3. It can be seen that the final models performed quite well for all sampled height classes. The prediction errors for the all strawberries equation can be only evaluated for trees less than 25 m in height due to these tree species occupy always the subdominant vertical layers when mature and are small-sized trees [51]. prediction errors for the all strawberries equation can be only evaluated for trees less than 25 m in height due to these tree species occupy always the subdominant vertical layers when mature and are small-sized trees [51].

Discussion
The species-specific growth models developed in the present study to predict dbh and total tree height of individual trees after a 5-year growth period included size variables (dbh1, h1), site factors (i.e., H0), and competition variables (BAL) in their formulation. In our study, the most important contribution of the variance explained by equations 2 and 3 came from size variables, followed by competition and site factors, respectively. These variables have also been selected by other authors for individual dbh, basal area and height increment models (e.g., [11,15,40,47,52,53]), although the most important variable has largely varied between ecosystems and species (e.g., [5,15,54,55]).
Our observed significant positive role of tree size on tree growth for all species agrees with studies that have observed that larger trees tend to grow faster (e.g., [15,17,56]). For example Stephenson et al. [57] presented a global analysis showing that lager trees had higher growth rates than smaller trees. Higher growth rates of larger trees are likely associated with increased access to resources; for example, larger trees generally have prominent position in the canopy, greater leaf area, and therefore have more access to light and other resources (e.g., [17]). In contrast, other studies have reported that growth rates increase to a maximum with tree size and then either levels off or decreases for the largest trees, as a result of reduced efficiency, accumulated stress, or senescence [16,17,55,[58][59][60]. We did not observe a decrease in tree growth for the largest trees for our observed diameter range, which included diameters of up to 50-70 cm for most of the species. Some authors have found declines above this diameter range, sometimes at ages above 50 years (e.g [55,61]). Such very mature and larger trees are relatively scarce in most forests under management schemes and our sample was not big enough to find them.
Significant competition and site effects associated with BAL and H0 were observed for most of the studied tree species both for individual dbh and height growth models. For dbh growth models, only five tree species showed no competition effect, while a non-significant site effect was found for

Discussion
The species-specific growth models developed in the present study to predict dbh and total tree height of individual trees after a 5-year growth period included size variables (dbh 1 , h 1 ), site factors (i.e., H 0 ), and competition variables (BAL) in their formulation. In our study, the most important contribution of the variance explained by equations 2 and 3 came from size variables, followed by competition and site factors, respectively. These variables have also been selected by other authors for individual dbh, basal area and height increment models (e.g., [11,15,40,47,52,53]), although the most important variable has largely varied between ecosystems and species (e.g., [5,15,54,55]).
Our observed significant positive role of tree size on tree growth for all species agrees with studies that have observed that larger trees tend to grow faster (e.g., [15,17,56]). For example Stephenson et al. [57] presented a global analysis showing that lager trees had higher growth rates than smaller trees. Higher growth rates of larger trees are likely associated with increased access to resources; for example, larger trees generally have prominent position in the canopy, greater leaf area, and therefore have more access to light and other resources (e.g., [17]). In contrast, other studies have reported that growth rates increase to a maximum with tree size and then either levels off or decreases for the largest trees, as a result of reduced efficiency, accumulated stress, or senescence [16,17,55,[58][59][60]. We did not observe a decrease in tree growth for the largest trees for our observed diameter range, which included diameters of up to 50-70 cm for most of the species. Some authors have found declines above this diameter range, sometimes at ages above 50 years (e.g [55,61]). Such very mature and larger trees are relatively scarce in most forests under management schemes and our sample was not big enough to find them.
Significant competition and site effects associated with BAL and H 0 were observed for most of the studied tree species both for individual dbh and height growth models. For dbh growth models, only five tree species showed no competition effect, while a non-significant site effect was found for thirteen tree species. In contrast, for height growth models a non-significant competition effect associated with BAL was found for eleven tree species, and only four of the studied species showed non-significant site factor effects. This appears to indicate that competition is a more important key driving factor for individual dbh growth in the study area than site factor (i.e., H 0 ), while the latter appears to be more appropriate for predicting height growth, as it significantly contributes to describing this variable. This higher sensitivity of diameter growth to competitive stress is probably explained because of a higher dependence of this process on light interception, whereas height growth is less affected (e.g., [62][63][64]).
Our results agree with previous studies that have shown improvements in diameter growth prediction for different species when a competition index was included (e.g., [36,55,[65][66][67][68]). Competition is an important driver of tree growth since the degree of shadowing emanating from competing neighbors substantially determines the amount and quality of solar radiation that is available for the target tree [64,69], together with the role of neighboring plants in competing for growing space and resources including water, and soil nutrients [17,70]. The consideration of competition after tree size for explaining diameter growth in our study agrees with previous studies that have also found that competition indices were the second most important variables, after tree size, in explaining individual tree diameter growth, for a variety of tree species. For example, in the studies of Wykoff [40], Monserud and Sterba [11] or Crecente-Campo et al. [46], the main variables in explaining growth were size variables (dbh or crown ratio), followed by the BAL competition index. In contrast, other studies have found competitive position, measured by BAL or modifications of BAL, to be the strongest individual predictor of diameter growth in both even-aged and uneven-aged stands [5,58,71,72].
In our study, BAL index explained a higher variance of tree growth than BALMOD or CCF indices. These results are similar to the findings of Crecente-Campo et al. [46], who found that BAL showed better results than BALMOD in terms of predicting basal area growth of Scots pine in NW Spain, unlike results from other studies where the latter has outperformed BAL (e.g., [43,73]), including a previous study in even-aged Pinus cooperi stands in the area of study by Corral-Rivas et al. [36] where BALMOD was selected.
In this study, significant site effects associated with H 0 were found for seven out of the eight developed individual dbh growth models for the pine species. Also, a significant H 0 effect was observed for the height growth models in all pine species, and indeed in almost all the studied tree species. Thus, H 0 appears to be an appropriate variable for predicting individual dbh and height growth for the study area, as it significantly contributes to describe individual tree growth within different stands. This finding is very interesting because more than 80% of Durango's forests are covered by mixed and uneven-aged stands, for which the most traditional growth modeling variables such as site index and age which appear to be important in describing the growth of individual trees within different stands, at least in even-aged forests [47]) are generally not available [74]. H 0 has also previously been included in generalized h-d models developed for some of the tree species studied in this paper in the same study area (e.g., [38,75]). Therefore, the inclusion of H 0 in the models might be used to indicate the stage of development of a stand, thus compensating better sites with lower age [11,76]. Nevertheless, the effect of dominant height on growth varied between species, similar to previous studies that have documented that the effect of site quality on tree growth was not consistent across species (e.g., [5]).
Individual tree species in forests have adapted to specific site conditions and are capable of responding to a range of complex environmental challenges to improve their chances of survival, reproduction, and growth. Some species are shade tolerant, others are light demanding. Some require high humidity and others are adapted to dry conditions [26]. In our study some trees species did not show significant competition or site factor effects, and therefore it is interesting to examine the possible reasons for the observed lack of sensitivity to those factors.
The species Quercus urbanii, Arbutus arizonica, A. tessellata, A. xalapensi and A. madrensis showed a non-significant effect of tree competition on dbh growth. According to Lujan-Soto et al. [51], these tree species occupy always the subdominant vertical layer when are mature, and consequently are shade-tolerant and possibly less sensitive to competition for light. In addition, Quercus crassifolia, Q. fulva, Q. laeta, Q. mcvaughii, Q. urbanii, Abies durangensis, Arbutus arizonica and A. xalapensis showed non-significant competition effects (i.e., BAL) on height growth, also being species adapted to the subdominant strata [51]. Similarly, in the area of study, the previous study by Quiñonez et al. [3] found that oak species seemed to have an advantageous competition in the mixed-species stands, since their growth was not significantly influenced by competition from other species. These results agree with observations from other studies in mixed-species conifer forests that have noted a lower sensitivity to competition for the more shade-tolerant species (e.g., [13,15,55,64]).
Non-significant site factor effects (i.e., H 0 ) on diameter growth were observed for Pinus herrerae, Quercus durifolia, Q. jonesii, Q. laeta, Arbutus bicolor and the group of Arbutus spp., which are tree species that are found growing in a vast geographic area of contrasting soils and site conditions [26]. The large plasticity of these tree species may suggest that such variation in dbh growth was not possible to be explained with dominant height [77,78], given their successful capacity of adaptation to very varied site quality conditions. On the other hand, the species Abies durangensis, Alnus firmifolia, and Cupressus lusitanica, which are species found in a very limited range of site conditions, also showed no significant competition effects on diameter growth. These are the most shade-tolerant species considered in our study which are growing only in a few plots with high humidity [26]. Similarly, Quercus radiata and Alnus firmifolia, which occurred on a very limited range of site conditions [26], showed no significant effects of site quality on tree height growth. The fact that these tree species occur on a very limited range of site conditions, possibly limited the inclusion of dominant height as a diameter growth predictor variable, given the low variability in this variable covered by the natural distribution of these species in our study area.
Future studies expanding the range of observed environmental conditions should also analyze the role of climatic and edaphic variables on tree growth, allowing to integrate the role of climatic or soil limitations on tree growth predictions (e.g., [16,17,58,79,80]) for improving our decision making capacity under climatic or management scenarios.
Although individual-tree dbh and height growth models as the ones developed in this paper are key components in many forest growth and yield prediction systems, their use would have restricted applicability without other model components such as ingrowth and mortality [11]. Hopefully, these model components will be developed for the studied tree species in the future. Thus, a reliable prediction system will be available for managing these natural species-rich and irregular forests in Durango, Mexico. However, the growth models developed in our study represent useful operational tools to support decision making for a comprehensive variety of tree species in Durango's forests. They may provide good simulation of growth (diameter, height, volume) for short term projections, provide detailed information about stand structure development (diameter and height distribution), calculate biomass and carbon stocks projection, and allow considerations of a wide variety of silvicultural treatments/prescriptions, among others.

Conclusions
The predictive variables used in the models to predict dbh and total tree height after a 5-year growth period were tree initial dbh and total tree height, respectively, the dominant height as a site factor (H 0 ), and the sum of the basal area in trees larger than the subject tree as a competition index (BAL). Some of the studied trees species did not show significant competition or site factor effects in the developed models, being the species-specific shade-tolerance and plasticity, respectively, the main reasons for the observed lack of sensitivity to those factors. Because these new species-specific growth models were developed with a large and unique dataset provided by and extensive network of stem-mapped and re-measured plots that covers a large variation of individual tree growth, they can be recommended not only for the forests of Durango but also for the surrounding regions with similar growth conditions. For other tree species that exist in the study area, but that were not considered in the study due to the lack of enough data, the generic models fitted on the combined data sets could be used for dbh and total tree height growth prediction. It is recommended that further studies on individual-tree growth for the studied species and species groups are continued with the aim to develop a reliable prediction system for managing these natural species-rich forests. Future models should analyze the role of species-specific ecological aspects on individual tree growth, once this knowledge has been generated.