Next Article in Journal
Phenolic Compounds Regulating the Susceptibility of Adult Pine Species to Bursaphelenchus xylophilus
Next Article in Special Issue
Estimation of Aboveground Vegetation Water Storage in Natural Forests in Jiuzhaigou National Nature Reserve of China Using Machine Learning and the Combination of Landsat 8 and Sentinel-2 Data
Previous Article in Journal
A Proposal for a Forest Digital Twin Framework and Its Perspectives
Previous Article in Special Issue
Crown Profile Modeling and Prediction Based on Ensemble Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Quantile Mixed-Effects Models for Prediction of the Maximum Crown Width of Fagus sylvatica L., Pinus nigra Arn. and Pinus brutia Ten.

by
Dimitrios I. Raptis
1,*,
Vassiliki Kazana
1,
Stavros Kechagioglou
2,
Angelos Kazaklis
3,
Christos Stamatiou
1,
Dimitra Papadopoulou
1 and
Thekla Tsitsoni
2
1
Department of Forest and Natural Environment Sciences, International Hellenic University, 66100 Drama, Greece
2
Department of Forestry and Natural Environment, Aristotle University of Thessaloniki, P.O. Box 262, 54124 Thessaloniki, Greece
3
OLYMPOS-Centre for Integrated Environmental Management, Kalamaria, 55132 Thessaloniki, Greece
*
Author to whom correspondence should be addressed.
Forests 2022, 13(4), 499; https://doi.org/10.3390/f13040499
Submission received: 14 February 2022 / Revised: 18 March 2022 / Accepted: 21 March 2022 / Published: 23 March 2022
(This article belongs to the Special Issue Modeling of Forest Tree and Stand Parameters)

Abstract

:
In the current study, a novel approach combining quantile regression with nonlinear mixed-effects (QR-NLME) modeling was applied to predict the maximum crown width (cwmax) of three economically important forest species—the European beech (Fagus sylvatica L.), the black pine (Pinus nigra Arn.), and the Calabrian pine (Pinus brutia Ten.) at tree level. A power QR-NLME model was fitted first to a dataset including 1414 European beech trees obtained from 29 randomly distributed sample plots, 770 black pine trees from 25 sample plots, and 1880 Calabrian pine trees from 41 sample plots in Greece, to predict the cwmax at tree level. Additionally, a nonlinear mixed-effects model (NLME) was fitted to the same dataset to predict the average crown width at tree level for all species. In the second stage, the crown competition factor (CCF) was estimated based on the population average response of the cwmax predictions. The proposed approach presented sound results when compared with the outcomes of relevant models from other regions fitted to open-grown tree data, and therefore, it can be well implemented on clustered data structures, in cases of absence of open-grown tree data.

1. Introduction

Estimation of the basic crown dimensions at the tree or stand level is of high importance in forest management science [1]. This is mainly due to the fact that the tree crown is closely related to the total tree growth [2], the competition level [3], and the basic attributes at the stand level, which further affect the stand structure and composition, including the understory vegetation. Assuming that a tree’s growing space (GS) is approximately circular-shaped [4], the prediction of the GS of the tree is based on the crown width dimensions. However, the trees growing under competition status are expected to develop shorter crown widths in comparison to those growing as open trees. In order to capture the tree crown’s biological tendency, open-grown trees free of competition have been typically used in crown width modeling, leading to the establishment of the maximum crown width (cwmax) dimension. In this context, cwmax has been defined as the theoretical maximum crown diameter of an ith tree, with a specific DBHi, growing in stand conditions, which would be observed had the specific tree been grown as an open-grown tree free of competition. According to Krajicek et al. [3], the mathematical expression of the cwmax may lead to the evaluation of the competition at stand level, through the use of the crown competition factor (CCF) concept.
CCF has been advocated as an objective competition measure in terms of evaluating the optimal stocking density and the consequent application of thinning strategies [5]. When trees within a stand can no longer achieve the maximum crown area without being limited by the neighboring trees, CCF values exceed 100%, and the competition level increases [6]. For CCF values higher than 100%, trees have less crown area available for maximum tree growth, but stand growth is still maintained at high levels due to the development of efficient foliage on the top of the remaining trees on-site [7]. In this sense, a CCF value equal to 100% implies that the tree at the stand has the maximum area available as if it would have been grown in an open-spaced area with no competition [8]. As a consequence, silvicultural treatments may affect significantly the CCF value levels within a stand. The removal of a small number of large-diameter trees reduces the competition level to the same degree as would the removal of many small-diameter trees [9].
From a practical perspective, the importance of the CCF determination is not limited only to the applied silvicultural treatments and, consequently, to the forest management planning. Due to its mathematical expression, it has been used as a significant regressor in forest growth modeling, in an effort to capture the competition status within a forest stand (e.g., [10,11,12]). In addition, based on the fact that the CCF is a distance-independent competition measure, there is no need for prior estimations of within-tree distances, and therefore, it can be easily estimated as long as the cwmax value is provided, which is clearly an advantage in modeling processes.
However, several problems have been reported with regard to the cwmax determination. For example, the cwmax value for a well-defined tree is not always available, whereas the cwmax modeling process is often a relatively complex task, particularly due to field difficulties related to the extensive traveling required in order to identify trees that meet the definition of an “open-grown tree” [13]. Furthermore, the definition of the typical “open-grown tree” is not quite clear, and this may introduce certain subjectivity [1]. Alternative modeling approaches for the prediction of the cwmax of certain forest species have been already attempted. Condés and Sterba [14] used a large dataset from the second Spanish National Forest Inventory to propose a modeling method for cwmax by defining “open-grown” trees as those whose residuals fell within both the upper 90th percentile of a general, double-logarithmic model developed to predict the crown width and the tree height based on the DBH, and the lower 10th percentile of the log of total height.
The majority of the applied methods for cwmax model development have used the DBH as the basic independent variable. In cases in which isolated open-grown trees have served as primary sampling units, the modeling process involved mainly the ordinary least squares (OLS) approaches usually presenting relatively reliable results. However, in most cases, the field data for forest models are obtained from sample plot installations of various sizes and shapes, and it is now well established that for multiple measurements from individual sampling units, the OLS method is inappropriate. This is because it leads to underestimates of the covariance matrix of the parameter estimates and the residual variance of the OLS equation [15,16]. While a number of methods have been evaluated as alternatives to solve this problem, the mixed-effects modeling approach is advocated as one of the most advanced tools for the analysis of clustered data [17]. Therefore, the correlation between measurements should be identified even in the case of using alternatives methods for analyses, such as quantile regression (QR).
Based on the importance of the cwmax determination for forest management planning and silvicultural treatment purposes at the stand level, the aims of the current study were (i) the development of cwmax prediction models for three economically important species—the European beech, the black pine, and the Calabrian pine—by applying a novel quantile nonlinear mixed-effects modeling approach; (ii) the estimation of the corresponding CCF values at the species level; (iii) the comparison of the current study predictions with the outcomes of relevant models fitted to open-grown trees of other regions, in order to evaluate the results at species biological tendency level.

2. Materials and Methods

2.1. Data Sampling

The raw field data were obtained from 29 sample plots installed in the pure European beech forests of Rodopi Mountain-Range National Park (RMRNP), located in the Region of Eastern Macedonia and Thrace, in northern Greece, between 24°10′ E and 24°50′ E, 41°20′ N and 41°32′ N (Figure 1). The forests of this region account for 34% of the total growing stock of Greece’s forest area [18]. The sample plots were randomly distributed using the “create random points” command of the ArcGis 10.2 (ESRI, Redlands, CA, USA) environment by using as base a well-defined shapefile of species expansion. Each circular-shaped sample plot covered a total area of 1000 m2 (0.1 ha), and it was installed in the field using a fixed radius of 17.8 m at vertical projection. For all living tress the total height (h), the diameter at breast height (DBH), the crown base height (cbh), and the crown width (cw) were measured with a LaserAce 1000 Rangefinder (Trimble) instrument. The crown width was measured in two vertical directions using the Canopy Spread module of the same instrument. The total sample included 1414 individual European beech trees over a wide range of density and competition status.
The second dataset was composed of 770 trees obtained from 25 sample plots located in the pure black pine natural forest stands of the Volakas Drama area, in northern Greece, between 24°00′ E and 24°03′ E, 41°16′ N and 41°19′ N (Figure 1). For the field measurements, the same methodology as above was followed, but the installed sample plots were half the size, covering 500 m2 (0.05 ha) each and thus corresponding to a 12.6 m fixed radius per plot at vertical projection.
The third dataset was composed of 1880 Calabrian pine trees obtained from a total of 41 randomly distributed sample plots in northern Greece, of which 21 and 10 plots were located, respectively, in the peri-urban forests of Drama (between 24°09′ E and 24°11′ E, 41°09′ N and 41°10′ N) and Xanthi (24°52′ E and 24°54′ E, 41°08′ N and 41°09′ N) cities in the Region of Eastern Macedonia and Thrace and 10 plots in the peri-urban forest of Thessaloniki city in the Region of Central Macedonia, extending between 22°58′ E and 22°59′ E, 40°37′ N and 40°39′ N (Figure 1). The randomly distributed sample plots were circular-shaped of 1000 m2 (0.1 ha), each with a fixed radius of 17.8 m at vertical projection. Within these established sample plots, the same variables as those described previously for the other 2 datasets were measured at tree level, using similar equipment.
The European beech stands in the study area are managed for wood production. Thinnings in the form of selective cuttings are applied periodically between 20 and 30 years according to the approval by the Forest Service management plan. This silvicultural system, as part of the sustainable forest management concept, aims primarily at soil protection by ensuring continuous forest cover over time. The same system is applied in the natural black pine forests growing in the wider area. In the case of the Calabrian pine peri-urban forests, a different management system is followed aiming primarily for fire protection and recreation activity. This system is based on different combinations of undergrowth clearings, pruning of dead crown lower parts, and some thinnings at specific points (on both sides of roads, paths, ridges) in order to prevent fire from crowning.

2.2. Analysis

The analysis was partly based on the approach of Russell and Weiskittel [1], using a QR technique, which allows the estimation of the response variable, that is, the crown width, for any percentile of the data [19]. According to the same authors, large datasets consist of various crown widths corresponding to both, open- and forest-grown trees. In order to estimate the maximum potential crown width for a given species at a specific DBH, the 99th percentile can be fitted to represent the cwmax value that could be safely assumed similar to open-grown trees. However, in the current research, due to the clustered structure of the raw data and the inherent correlation between measurements, it was preferred to combine QR with nonlinear mixed-effects models (NLMEs), which have effectively been used for analyzing similar grouped data structures [20,21]. The proposed QR-NLME model is currently based on the implementation of the stochastic approximation of the expectation-maximization (SAEM) algorithm [22] and the asymmetric Laplace (AL) distribution, as it has been described by Galarza et al. [20] through the qrNLMM package in R (cran) open source statistical software. In addition, based on the sample size, the 90th percentile was finally fitted, which mostly corresponds to suggestions by Condés and Sterba [14].
Despite the fact that the linear model is the most commonly used for cwmax predictions [23], a nonlinear (power) model of the mathematical expression y = axb was selected in the current study, because the specific form has presented increased accuracy for crown width modeling of European beech [24], black pine [25], and Calabrian pine trees [26], and it can be easily linearized through logarithmic transformation. Russell and Weiskittel [1] also used the same model form for cwmax modeling of various species both conifers and hardwoods. In addition, a simple mixed-effects model was separately fitted to the same datasets aiming at an average (mean) crown width response. A detailed description of the mixed-effects fitting approach is provided by Calama and Montero [17]. In both those modeling procedures, it was assumed that both model parameters were composed of a random part, as long as convergence was feasible, following the suggestions of Pinheiro and Bates [27].
The mixed-effects modeling approach provides predictions for the population average (marginal) and for the subject-specific response. For the latter, additional data are required in order to calibrate the mixed-effects model at the local level (sample plot or stand). However, in some cases, the required data to localize the random part are difficult to obtain [28,29], and the population average response can be used instead, by setting the random effects equal to zero. Since data for open-grown trees are rarely available, the fixed part of the mixed-effects model was considered as a reasonable solution due to the low inter-subject variability, after percentile separation. In addition, the utility of the mixed-modeling approach in the specific case aimed at the population average response in an effort to capture the biological cwmax tendency at the species level, and therefore, local calibration was not necessary. The final cwmax models were compared graphically with those that were proposed by Hasenauer [30] for European beech and black pine in Austria, and by Kara and Topaçoğlu [31] for Calabrian pine in Turkey.
For the estimation of the CCF value at the final stage, the fixed part of the cwmax model was used as a potential unbiased solution at the population level. The statistical analysis was performed in the R (cran) open source statistical software environment, along with the qrNLMM, ggplot2, lmfor, and NLME libraries.

Mathematical Expression

The general formulation of the mixed-effects models is [17,20]
y i = f ( Φ i , x i ) + e i
and
Φ i = A i β p + B i b i
where y is the response variable of the ith subject, Φi is a parameter vector of dimension r, and f is the nonlinear function form. The ei term represents the independent and normally distributed vector of random errors. The parameter vector can be defined through Equation (2), where βp is the regression coefficient corresponding to the pth quantile, Ai and Bi are design matrices of dimensions r × d and r × q, for the fixed and random effects specific to each plot, and bi is a q-dimensional random effects vector associated with the ith subject. Hence, the pth quantile regression for the nonlinear mixed-effects model can be expressed as [20]
Q p = ( y i j | x i j ,   β p ) = f ( Φ i ,   x i j ) = f ( A i β p + B i b i ,     x i j ) + e i j
In addition, for the specific case, the mathematical form of the basic mixed-effects model was
c w = ( β 1 + b 1 ) D B H ( β 2 + b 2 ) + e
where β1, β2 is the fixed part and b1, b2 is the random part of the model, cw is the crown width (m), and DBH is the diameter at breast height (cm).
Assuming that the vertical projection of a tree crown is approximately circular shaped, the general maximum crown projection area (cpa) of the ith tree crown can be estimated from the following type:
c p a i = π 4 c w i   m a x 2     o r     c p a i = π 4 ( β 1 D B H i β 2 ) 2
where β1, β2 is the fixed part of the quantile nonlinear mixed-effects model (QR-NLME) when no auxiliary information is available, and cpai is the vertical crown projection area of ith tree (m2).
The fixed part of the model may also be localized per plot through a limited sample of data, as one of the main attributes of mixed-effects models, which is not recommended in this specific case.
The CCF of a forest stand including n trees, growing on a total area A, is expressed as the quotient between the sum of the potential (maximum) crown projection areas cpa1 + cpa2 + … cpan of all living trees and the total area A, multiplied by 100 [6] as follows:
C C F = 1 A i = 1 n c p a i 100
where CCF is the crown competition factor (dimensionless); A is the total stand area (m2), which may represent the plot area [9,32]; n is the total number of living trees.

3. Results

The descriptive statistics of the total sample (n = 4064) per species are presented in Table 1.

3.1. Convergence

The QR-NLME model successfully converged in all cases using the default values for maximum iterations (MaxIter = 500) and for the cut point (c = 0.25), which determines the percentage of initial iterations with no memory (Figure 2). It should be noted that convergence is not always ensured, and different starting values of parameter combinations should be tested instead.
Based on the variance of the parameters at plot level (Figure 3 and Table 2), it was assumed that both parameters of the average (mean) mixed-effects nonlinear model contain random parts, and the model was not simplified further. In addition, for the same model, a power form variance function was used in order to stabilize the error’s variance, since it presented the best fitting ability during the preliminary research stage. The variance function was of the form V a r ( e ) = σ 2 D B H 2 k , where σ2 is the residual variance of the estimated model, and k is the variance function coefficient. The average model converged for the three species and the residual plots revealed no heteroscedastic trends across the total range of the fitted values (Figure 4), after weighting. The value of the k coefficient is also presented in Table 2.

3.2. Fitting Procedure

The results of the fitting procedure of both models for all species are presented in the following Table 2.
All model parameters were significantly different from zero according to the approximate t-criterion (p < 0.001 for all species). The Fitting index from Table 2, which is similar to R2, lead to the conclusion that the mixed-effects model with the DBH as the sole independent variable managed to explain the 74.5%, 68.8%, and 74.8% of the crown width variance in European beech, black pine, and Calabrian pine, respectively, which is a very important finding in terms of practical use in field conditions. Both, the error and the prediction bias, ranged to low levels (Table 2).

3.3. Graphical Representation

Figure 5 depicts the graphical performance of the developed cwmax models separately for each of the three forest species. For comparison reasons, the regression (dashed/dotted) lines of the equivalent suggested cwmax models for the same species were added within the same graphs. It should be noted, however, that the selected models were developed using data from open-grown trees in Austria and Turkey, in contrast with the current research, in which the tree data of various densities and competition status were clustered at the sample plot level. The population responses (fixed effects) of the average (mean) models are graphically presented in dashed-line style. As it can be observed, a nonlinear relationship between the crown width and the DBH is revealed, especially for the European beech tree data, and to a lesser extent for the Calabrian pine tree data. The basic model of power form presented increased flexibility and managed to perform well against almost linear distributed data, as it can be concluded from the right (b) section of Figure 5, which corresponds to the black pine tree data.

3.4. The Crown Competition Factor (CCF) Estimation

Using the population-average cwmax model and Equations (5) and (6), CCF was estimated for each species at the sample plot level. The results are presented in the following Table 3.
The European beech presented the highest levels of competition in the specific area, where it reaches its optimal growth. For the pine species, the competition levels were relatively close, with black pine forest stands presenting slightly higher CCF values than those of Calabrian pine.

4. Discussion

In order to deal with the problem of the correlation of measurements from the same subjects, Sun et al. [21] introduced the quantile regression for the linear mixed-effects model into crown profile modeling, with encouraging results. Furthermore, the quantile regression provides the opportunity of estimating the conditional distribution of dependent variables while assessing the effects of regressors at various quantiles [33]. The same authors compared the quantile regression method and the mixed-effects modeling technique for height prediction of two forest species in Turkey. Cao and Wang [34] also compared the two methods to calibrate a taper equation for loblolly pine (Pinus taeda L.) species in Louisiana. To the best of our knowledge, this is the first time that these novel methods are linked together with a nonlinear model form in order to predict the cwmax dimension at tree level, accounting for clustered data structures from sample plots. This approach is also recommended in the case of hierarchical measurements over time, which is common in the case of permanent sample plots installation. From the modeling perspective, there are no similar approaches on cwmax prediction, and the results must be related to relevant studies of different modeling approaches.
Hasenauer [30] developed cwmax models for European beech and black pine in Austria, while Kara and Topaçoğlu [31] suggested similar models for the Calabrian pine in Turkey. The results of the current study present some notable similarities with the outcomes of the reported models, as can be seen in Figure 5. The differences appearing between the models’ predictions were expected since the data from the previous studies represented samples from completely different regions in Austria and Turkey. Similar findings were reported by Russell and Weiskittel [1], who compared their predictions with relevant models from other regions in the US. In addition, the open conditions for open-grown trees could have affected the crown width expression because those types of trees are usually exposed to climatic factors (wind, snow, etc.), compared with stand-grown trees, which are certainly affected by the competition status for light. The largest differences were observed on the European beech cwmax curves, mainly for small-diameter trees. The comparison between the Calabrian pine models leads to the conclusion that the corresponding curves had almost identical shapes, with the current model to provide slightly increased cwmax predictions, over a wide range of diameter classes. For the black pine, the cwmax curves are relatively similar, and the differences are very small, which might have a genetic explanation [35]. Hence, the proposed method seems to produce reliable cwmax predictions for the two conifer species, while additional research is needed in the case of European beech, for which data from the same region must be compared.
For the development of the average crown width model, the mixed-effects approach presented sound results, explaining the largest part of crown variation according to the corresponding fitting index (Table 2). The increased flexibility of the simple mixed-effects model has also been reported by other researchers [24,36], and it is based on the fact that the random part of the model was able to capture the effect of other factors rather than DBH, which were not included as additional regressors into the fitting model form. From that perspective, the random effects can be correlated with variables at a tree or stand level [37], leading to the development of generalized models that may, in turn, be expanded to mixed-effects and explain the functional mechanisms of crown allometry. The inherent simplicity of the proposed average crown width models is rather an advantage for practical use purposes since the DBH is the most common tree variable measured in forest inventories and management plans. Transferability on potential model applications should be based on the restrictions dictated by the model conditions, particularly those referred to the specific stand densities, as they are described in Table 1.
In order to select the 90% quantile for the QR analyses, different percentage values of quantiles were tested. Since CCF is assumed independent from site index [9], the value of 400 has been assumed as a reasonable upper bound, representing a stable condition between growth and mortality [9,14], although significantly higher values have been reported only as maximum values. In this context, Arney [9] reported a maximum CCF value of 701 for Douglas-fir (Pseudotsuga menziesii) naturally regenerated forests. The 99% QR led to unreasonable mean CCF values, especially for European beech species, which is further endorsed by the fact that these forests are subject to various management regimes. In addition, the 99% quantile could have been affected by outliers, common in large datasets, since no smoothing method has been applied. Similar findings were derived using 95% as quantile input. The 90% quantile presented reasonable results, which were further verified for all species according to Figure 5.
The analysis showed increased competition status for the European beech trees, while the corresponding values for the pine species were found to be at lower levels (Table 3). Since the CCF is supposed to be independent of the stand structure [13], those discrepancies may be attributed to the differences in the estimated cwmax for the studied species, provided that the stem density alone led to opposite conclusions (see the number of trees per ha in Table 1). From Table 2, it is inferred that, for example, a tree with 40 cm DBH is expected to develop a 12.52 m cwmax if it is a beech open-grown tree, 8.71 m if it is a black pine tree, and 9.41 m if it is a Calabrian pine tree. Hence, this may potentially lead to increased CCF values for the European beech plots, as it is depicted in Table 3. On the other hand, the difference in the competition status between black pine and Calabrian pine was expected to be larger, due to the differences between the management regimes, which are typically applied for the two species. Black pine forests are usually managed for wood production [38], while the Calabrian pine peri-urban forests are mainly managed for protection and recreation [39]. However, the adoption of selective cuttings as the primary silvicultural treatment in black pine forests has ensured the continuity of the forest cover, as well as the regeneration of this type of forest in a sustainable way.
Critical issues during the CCF application are the threshold values that have been proposed in order to schedule silvicultural treatments to regulate the stand competition in the frame of sustainable forest management. Arney [9] proposed the threshold value of 200 for the crown recession of coastal Douglas fir and 300 for suppression mortality. Castaldi et al. [40] reported increased values of CCF for the same species in Italy, reaching 619 in some forest stands. Since there are no similar values for the species of the current study, the CCF measure must be correlated with mortality and productivity aspects per species in a featured research study.
The application of the SAEM algorithm on a relatively large dataset was a complex process; thus, some convergence problems emerged. Indeed, for the Calabrian pine dataset, the model converged after 393 iterations within a total processing time of 36.7 min. For the European beech dataset, the model converged after 458 iterations and 30.3 min, and for the black pine after 454 iterations and 19.1 min of processing time. In addition, the model failed to converge in two cases, and the procedure had to start from the beginning.

5. Conclusions

The nonlinear quantile mixed-effects approach developed reliable cwmax models for two conifer species in northern Greece—the black pine (Pinus nigra Arn.) and the Calabrian pine (Pinus brutia Ten.). The proposed method can be safely applied in cases of the absence of open-grown tree data, while at the same time, it rectifies the problem of correlated observations from clustered data structures. In the case of the European beech (Fagus sylvatica L.) species, the proposed method should be applied with caution, and additional research is required.

Author Contributions

Conceptualization, D.I.R.; methodology, D.I.R. and V.K.; software, D.I.R.; formal analysis, D.I.R., V.K. and A.K.; investigation, D.I.R.; resources, D.I.R. and V.K.; data curation, D.I.R., S.K., C.S. and D.P.; writing—original draft preparation, D.I.R. and V.K.; writing—review and editing, A.K. and T.T.; project administration, V.K.; funding acquisition, V.K. and A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partly funded by the European Regional Fund and the Eastern Macedonia and Thrace Region of Greece NSRF (ESPA) 2007-2013, MIS 430008.

Data Availability Statement

The data are available on reasonable request.

Acknowledgments

The authors would like to thank the three anonymous reviewers of the manuscript, who provided insightful comments that significantly improved our paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Russell, M.; Weiskittel, A. Maximum and Largest Crown Width Equations for 15 Tree Species in Maine. North. J. Appl. For. 2011, 28, 84–91. [Google Scholar] [CrossRef] [Green Version]
  2. Monserud, R.A.; Sterba, H. A basal area increment model for individual trees growing in even- and uneven-aged forest stands in Austria. For. Ecol. Manag. 1996, 80, 57–80. [Google Scholar] [CrossRef]
  3. Krajicek, J.E.; Brinkman, K.A.; Gingrich, S.F. Crown competition—A measure of density. For. Sci. 1961, 7, 35–42. [Google Scholar]
  4. Foli, E.G.; Alder, D.; Miller, H.G.; Swaine, M.D. Modelling growing space requirements for some tropical forest tree species. For. Ecol. Manag. 2003, 173, 79–88. [Google Scholar] [CrossRef]
  5. Fernandez-Moya, J.; Urbán-Martínez, I. Estimation of crown competition factor for hybrid walnut (Juglans × intermedia) Mj209xRa planted forests in Spain. Ann. Silvic. Res. 2020, 44, 24–29. [Google Scholar] [CrossRef]
  6. Pretzsch, H. Forest Dynamics, Growth and Yield; Springer: Berlin/Heidelberg, Germany, 2010; p. 664. [Google Scholar]
  7. Hall, R.B. Use of the crown competition factor concept to select clones and spacings for short-rotation woody crops. Tree Physiol. 1994, 14, 899–909. [Google Scholar] [CrossRef]
  8. Strub, M.R.; Vasey, R.B.; Burkhart, H.E. Comparison of diameter growth and crown competition factor in Loblolly Pine plantations. For. Sci. 1975, 21, 427–431. [Google Scholar]
  9. Arney, J.D. A modeling strategy for the growth projection of managed stands. Can. J. For. Res. 1985, 15, 511–518. [Google Scholar] [CrossRef]
  10. Rijal, B.; Weiskittel, A.R.; Kershaw, J.A. Development of height to crown base models for thirteen tree species of the North American Acadian Region. For. Chron. 2012, 88, 60–73. [Google Scholar] [CrossRef] [Green Version]
  11. Temesgen, H.; Lemay, V.; Mitchell, S.J. Tree crown ratio models for multi-species and multi-layered stands of southeastern British Columbia. For. Chron. 2005, 81, 133–141. [Google Scholar] [CrossRef]
  12. Boisvenue, C.; Temesgen, H.; Marshall, P. Selecting a small tree height growth model for mixed-species stands in the southern interior of British Columbia, Canada. For. Ecol. Manag. 2004, 202, 301–312. [Google Scholar] [CrossRef]
  13. Weiskittel, A.R.; Hann, D.W.; Kerhsaw, J.A., Jr.; Vanclay, J.K. Forest Growth and Yield Modeling; Wiley: Oxford, UK, 2011. [Google Scholar]
  14. Condés, S.; Sterba, H. Derivation of compatible crown width equations for some important tree species of Spain. For. Ecol. Manag. 2005, 217, 203–218. [Google Scholar] [CrossRef]
  15. West, P.W.; Ratkowsky, D.A.; Davis, A.W. Problems of hypothesis testing of regressions with multiple measurements from individual sampling units. For. Ecol. Manag. 1984, 7, 207–224. [Google Scholar] [CrossRef]
  16. West, P.W.; Davis, A.W.; Ratkowsky, D.A. Approaches to regression analysis with multiple measurements from individual sampling units. J. Stat. Comput. Simul. 1986, 26, 149–175. [Google Scholar] [CrossRef]
  17. Calama, R.; Montero, G. Interregional nonlinear height–diameter model with random coefficients for stone pine in Spain. Can. J. For. Res. 2004, 34, 150–163. [Google Scholar] [CrossRef] [Green Version]
  18. Kazana, V.; Kazaklis, A.; Raptis, D.; Stamatiou, C. A combined multi-criteria approach to assess forest management sustainability: An application to the forests of Eastern Macedonia & Thrace Region in Greece. Ann. Oper. Res. 2020, 294, 321–343. [Google Scholar] [CrossRef]
  19. Koenker, R.; Hallock, K.F. Quantile regression. J. Econ. Perspect. 2001, 15, 143–156. [Google Scholar] [CrossRef]
  20. Galarza, C.E.; Castro, L.M.; Louzada, F.; Lachos, V.H. Quantile regression for nonlinear mixed effects models: A likelihood based perspective. Stat. Pap. 2018, 61, 1281–1307. [Google Scholar] [CrossRef]
  21. Sun, Y.; Gao, H.; Li, F. Using Linear Mixed-Effects Models with Quantile Regression to Simulate the Crown Profile of Planted Pinus sylvestris var. Mongolica Trees. Forests 2017, 8, 446. [Google Scholar] [CrossRef] [Green Version]
  22. Delyon, B.; Lavielle, M.; Moulines, E. Convergence of a stochastic approximation version of the EM algorithm. Ann. Stat. 1999, 27, 94–128. [Google Scholar] [CrossRef]
  23. Paine, D.P.; Hann, D.W. Maximum Crown-Width Equations for Southwestern Oregon Tree Species; Research Paper 46; Forest Research Laboratory, Oregon State University: Corvallis, OR, USA, 1982; p. 21. [Google Scholar]
  24. Sharma, R.P.; Vacek, Z.; Vacek, S. Individual tree crown width models for Norway spruce and European beech in Czech Republic. For. Ecol. Manag. 2016, 366, 208–220. [Google Scholar] [CrossRef]
  25. Raptis, D.; Kazana, V.; Kazaklis, A.; Stamatiou, C. A Crown Width-Diameter Model for Natural Even-Aged Black Pine Forest Management. Forests 2018, 9, 610. [Google Scholar] [CrossRef] [Green Version]
  26. Avsar, M.D. The relationships between diameter at breast height, tree height and crown diameter in Calabrian pines (Pinus brutia ten) of Baskonus mountain, Kahramanmaras, Turkey. J. Biol. Sci. 2004, 4, 437–440. [Google Scholar]
  27. Pinheiro, J.C.; Bates, D.M. Mixed-Effects Models in S and S-PLUS; Spring: New York, NY, USA, 2000. [Google Scholar]
  28. Bronisz, K.; Zasada, M. Comparison of Fixed- and Mixed-effects Approaches to Taper Modeling for Scots Pine in West Poland. Forests 2019, 10, 975. [Google Scholar] [CrossRef] [Green Version]
  29. de Miguel, S.; Mehtätalo, L.; Shater, Z.; Kraid, B.; Pukkala, T. Evaluating marginal and conditional predictions of taper models in the absence of calibration data. Can. J. For. Res. 2012, 42, 1383–1394. [Google Scholar] [CrossRef]
  30. Hasenauer, H. Dimensional relationships of open-grown trees in Austria. For. Ecol. Manag. 1997, 96, 197–206. [Google Scholar] [CrossRef]
  31. Kara, F.; Topaçoğlu, O. Onset of canopy closure for black pine, Turkish red pine and Scots pine forests. J. For. Sci. 2018, 64, 224–229. [Google Scholar]
  32. Alegria, C. Simulation of silvicultural scenarios and economic efficiency for maritime pine (Pinus pinaster Aiton) wood-oriented management in centre inland of Portugal. For. Syst. 2011, 20, 361–378. [Google Scholar]
  33. Özçelik, R.; Cao, Q.V.; Trincado, G.; Göçer, N. Predicting tree height from tree diameter and dominant height using mixed-effects and quantile regression models for two species in Turkey. For. Ecol. Manag. 2018, 419–420, 240–248. [Google Scholar] [CrossRef]
  34. Cao, Q.V.; Wang, J. Evaluation of Methods for Calibrating a Tree Taper Equation. For. Sci. 2015, 61, 213–219. [Google Scholar] [CrossRef] [Green Version]
  35. Scaltsoyiannes, A.; Rohr, R.; Panetsos, K.P.; Tsaktsira, M. Allozyme frequency distributions in 5 European populations of black pine (Pinus nigra Arnold). 1. Estimation of genetic variation within and among populations. 2. Contribution of isozyme analysis to the taxonomic status of the species. Silvae Genet. 1994, 43, 20–30. [Google Scholar]
  36. Temesgen, H.V.J.; Monleon, V.J.; Hann, D.W. Analysis and comparison of nonlinear tree height prediction strategies for Douglas-fir forests. Can. J. For. Res. 2008, 38, 553–565. [Google Scholar] [CrossRef] [Green Version]
  37. Raptis, D.I.; Kazana, V.; Onisiforou, N.; Stamatiou, C.; Kazaklis, A. Height Allometry of Pinus nigra Arn. in Troodos National Forest Park, Cyprus. Sustainability 2021, 13, 5998. [Google Scholar] [CrossRef]
  38. Raptis, D.I.; Kazana, V.; Kazaklis, A.; Stamatiou, C. Mixed-effects height–diameter models for black pine (Pinus nigra Arn.) forest management. Trees 2021, 35, 1167–1183. [Google Scholar] [CrossRef]
  39. Christopoulou, O.; Polyzos, S.; Minetos, D. Peri-urban and urban forests in Greece: Obstacle or advantage to urban development? Manag. Environ. Qual. 2007, 18, 382–395. [Google Scholar] [CrossRef]
  40. Castaldi, C.; Vacchiano, G.; Marchi, M.; Corona, P. Projecting Nonnative Douglas Fir Plantations in Southern Europe with the Forest Vegetation Simulator. For. Sci. 2017, 63, 101–110. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Distribution of the sample plots in the study areas.
Figure 1. Distribution of the sample plots in the study areas.
Forests 13 00499 g001
Figure 2. Graphical summary of successful convergence for the fixed-effect parameters of the cwmax model, variance components of the random effects, and nuisance parameters, extracted from the qrNLMM package for European beech (a), black pine (b), and Calabrian pine (c) data. The vertical dashed lines delimit the beginning of the almost sure convergence, as defined by the cut-point parameter (𝑐 = 0.25).
Figure 2. Graphical summary of successful convergence for the fixed-effect parameters of the cwmax model, variance components of the random effects, and nuisance parameters, extracted from the qrNLMM package for European beech (a), black pine (b), and Calabrian pine (c) data. The vertical dashed lines delimit the beginning of the almost sure convergence, as defined by the cut-point parameter (𝑐 = 0.25).
Forests 13 00499 g002
Figure 3. The 95% confidence intervals on the average power model parameters per sample plot for European beech (a), black pine (b), and Calabrian pine (c) data, based on the nlsList fit (NLME package). A considerable variability among plot parameter estimates for all species can be observed, revealing the need for introducing random effects to account for the between-plot variation.
Figure 3. The 95% confidence intervals on the average power model parameters per sample plot for European beech (a), black pine (b), and Calabrian pine (c) data, based on the nlsList fit (NLME package). A considerable variability among plot parameter estimates for all species can be observed, revealing the need for introducing random effects to account for the between-plot variation.
Forests 13 00499 g003aForests 13 00499 g003b
Figure 4. Standardized residual plot of the weighted average mixed-effect model, for European beech (a), black pine (b), and Calabrian pine (c) data. The large dots represent the mean values of the residuals according to 10 classes. The thick and thin vertical lines represent the 95% confidence interval of class mean and of one observation (mean ± SD), respectively. The red color indicates the lines that do not cross the baseline at y = 0.
Figure 4. Standardized residual plot of the weighted average mixed-effect model, for European beech (a), black pine (b), and Calabrian pine (c) data. The large dots represent the mean values of the residuals according to 10 classes. The thick and thin vertical lines represent the 95% confidence interval of class mean and of one observation (mean ± SD), respectively. The red color indicates the lines that do not cross the baseline at y = 0.
Forests 13 00499 g004
Figure 5. Comparative graphs of the developed models for cwmax prediction (solid line), along with a general model for the mean crown width (dashed line), for European beech (a), black pine (b), and Calabrian pine (c) tree data. The cwmax models by Hasenauer [30] for European beech and black pine in Austria, and by Kara and Topaçoğlu [31] for Calabrian pine in Turkey, are also graphically presented.
Figure 5. Comparative graphs of the developed models for cwmax prediction (solid line), along with a general model for the mean crown width (dashed line), for European beech (a), black pine (b), and Calabrian pine (c) tree data. The cwmax models by Hasenauer [30] for European beech and black pine in Austria, and by Kara and Topaçoğlu [31] for Calabrian pine in Turkey, are also graphically presented.
Forests 13 00499 g005
Table 1. Descriptive statistics of the total sample.
Table 1. Descriptive statistics of the total sample.
Fagus Sylvatica L. (n = 1414)MinMeanMaxSD
DBH (cm)0.224.4088.514.04
h (m)1.518.2734.46.85
cw (m)0.97.6019.62.91
Number of trees per ha220488880166
Basal area (m2·ha−1)11.5930.3248.528.75
Pinus nigra Arn. (n = 770)
DBH (cm)3.129.1777.212.71
h (m)1.916.4524.84.76
cw (m)0.34.7815.02.27
Number of trees per ha3206161060165
Basal area (m2·ha−1)26.8748.9482.5012.85
Pinus brutia Ten. (n = 1880)
DBH (cm)0.728.1577.99.84
h (m)1.714.4421.74.12
cw (m)0.46.4215.02.31
Number of trees per ha210498750129
Basal area (m2·ha−1)11.3932.8760.8812.71
Table 2. Model parameters and fitting statistics for all species.
Table 2. Model parameters and fitting statistics for all species.
Model Form β1DBHβ2
Fagus sylvaticaPinus nigraPinus brutia
Maximum Crown Width Model
β1 (SE)2.38356 (0.0692)0.53990 (0.0377)0.96984 (0.0423)
β2 (SE)0.44960 (0.0083)0.75397 (0.0231)0.61605 (0.0091)
var(b1)0.299720.046700.15569
var(b2)0.003920.013970.09540
cov(b1,b2)−0.03149−0.02432−0.03528
σ20.070320.051510.02862
AIC6017.143069.266375.88
Log likelihood−3002.57−1528.63−3181.94
Percentile909090
Average (Mean) Crown Width Model
β1 (SE)1.49139 (0.0810)0.31579 (0.0275)0.54825 (0.0356)
β2 (SE)0.52723 (0.0148)0.81158 (0.0245)0.71199 (0.0174)
var(b1)0.118630.004890.02402
var(b2)0.003380.002930.00413
cov(b1,b2)−0.01821−0.00265−0.00799
σ20.590470.067510.04278
k0.215970.474670.47917
AIC5201.232561.585474.78
Log likelihood−2593.61−1273.79−2730.39
Fitting index0.744590.687860.74843
Root mean square error1.47021.26491.0194
Bias0.0000−0.0068−0.0033
Table 3. The estimated CCF values per sample plot, based on the population average approach for cwmax estimation.
Table 3. The estimated CCF values per sample plot, based on the population average approach for cwmax estimation.
SpeciesCrown Competition Factor (CCF)
Min.MeanMax.SD
Fagus sylvatica218.3378.5598.689.9
Pinus nigra138.8245.2388.858.9
Pinus brutia86.5214.2376.465.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Raptis, D.I.; Kazana, V.; Kechagioglou, S.; Kazaklis, A.; Stamatiou, C.; Papadopoulou, D.; Tsitsoni, T. Nonlinear Quantile Mixed-Effects Models for Prediction of the Maximum Crown Width of Fagus sylvatica L., Pinus nigra Arn. and Pinus brutia Ten. Forests 2022, 13, 499. https://doi.org/10.3390/f13040499

AMA Style

Raptis DI, Kazana V, Kechagioglou S, Kazaklis A, Stamatiou C, Papadopoulou D, Tsitsoni T. Nonlinear Quantile Mixed-Effects Models for Prediction of the Maximum Crown Width of Fagus sylvatica L., Pinus nigra Arn. and Pinus brutia Ten. Forests. 2022; 13(4):499. https://doi.org/10.3390/f13040499

Chicago/Turabian Style

Raptis, Dimitrios I., Vassiliki Kazana, Stavros Kechagioglou, Angelos Kazaklis, Christos Stamatiou, Dimitra Papadopoulou, and Thekla Tsitsoni. 2022. "Nonlinear Quantile Mixed-Effects Models for Prediction of the Maximum Crown Width of Fagus sylvatica L., Pinus nigra Arn. and Pinus brutia Ten." Forests 13, no. 4: 499. https://doi.org/10.3390/f13040499

APA Style

Raptis, D. I., Kazana, V., Kechagioglou, S., Kazaklis, A., Stamatiou, C., Papadopoulou, D., & Tsitsoni, T. (2022). Nonlinear Quantile Mixed-Effects Models for Prediction of the Maximum Crown Width of Fagus sylvatica L., Pinus nigra Arn. and Pinus brutia Ten. Forests, 13(4), 499. https://doi.org/10.3390/f13040499

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