A Dynamical Model Based on the Chapman–Richards Growth Equation for Fitting Growth Curves for Four Pine Species in Northern Mexico

Abstract

Tree growth models describe the growth and development of forest ecosystems by considering how the dimensions of each simulated tree change within a certain time. These models have commonly used three growth parameters that describe various biological processes and behaviours, considering a sigmoid growth function: (i) the upper asymptote (θ1), which is the maximal yield indicated by a final dimension (such as the maximal stem diameter); (ii) the maximum specific growth rate (θ2), defined as the slope of the tangent at the inflexion point; and (iii) the time elapsed (θ3), defined by the intercept of this tangent with the abscissas. To the best of our knowledge, however, associations between the three parameters have not been documented for tree species. Using diameter growth data from pine trees located in typical mixed and uneven-aged pine-oak forests in the Sierra Madre Occidental, Mexico, our study aims were: (i) to quantify the putative associations between the three growth parameters and (ii) to test the accuracy of a proposed Hybrid Chapman-Richards growth model based on associations between the three growth parameters, but including only one single parameter, relative to the widely used Generalized Algebraic Difference Approach (GADA) based on the Chapman-Richards, Lundqvist and Hossfeld models and the Hybrid Weibull Model. For statistical comparison of the quality of the models, we used the mean relative percentage error, root mean square error, coefficient of determination and Akaike information criterion to assess the quality of the fit. Although the quality of the five growth models studied was similar, from a practical point of view, the proposed Hybrid Chapman-Richards Model (CR-H) is easier to apply than the other models and has a lower data collection and computational cost. The parameter of CR-H can be easily obtained, by measuring just the dominant trees, especially in coniferous forests with irregular ages. Moreover, in contrast to the Chapman-Richards-GADA factor χ0, when θ2 is assumed to be site-specific, the CR-H has always a closed-form solution.