Estimating crown width in degraded forest: A two-level nonlinear mixed-effects crown width model for Dacrydium pierrei and Podocarpus imbricatus in tropical China

Abstract Tropical forest degradation makes a major contribution to greenhouse gas emission. Crown width (CW) is one of the important predictors in forest growth and yield models that provide basic data for assessment of forest degradation. Precise method of estimating tree crown for two tropical tree species (Dacrydium pierrei Hickel and Podocarpus imbricatus Bl), which are the major species in the degraded coniferous mixed forests in the tropical China, is necessary. These forests play a pivotal role in maintaining ecosystem functions, but are under the threat of severe degradation in recent years, and none of the studies has provided focus to these forests. We developed a nonlinear mixed-effects CW model using the permanent sample plot data acquired from D. pierrei and P. imbricatus forests. A number of tree- and stand-level variables were evaluated for their potential contribution to the CW variations, and included only highly significant ones in the model. The random effects at the levels of both sample plots and stands with different site quality class (blocks) were included in the CW model through mixed-effect modeling, and resulting model is therefore a two-level nonlinear mixed-effects model. Leave-one-out cross-validation was applied to evaluate the models. Among several predictor variables, diameter at breast height (DBH), height-to-DBH ratio (HDR), and height to crown base (HCB) contributed relatively highly to the CW variations. Dummy variable was introduced into the model to differentiate CW variations of two tree species. Correlations of CW and predictor variables significantly decreased when random effects at both the block and sample plot levels were included. We calibrated the nonlinear mixed effects CW model following the empirical best linear unbiased prediction theory, using four strategies of selecting CW trees per sample plot (largest, medium-sized, smallest trees and randomly selected trees) and fifteen sample sizes (one to fifteen trees). The prediction accuracy increased with increasing number of trees per sample plot, except the smallest trees, but the largest increase occurred with three largest trees used in calibration. This article emphasized more on modeling methodology, which can be applied to construct CW models for any forest elsewhere including degraded forest in the tropics.

[1]  D. Pothier,et al.  Impact of dominant tree dynamics on site index curves , 2003 .

[2]  R. Monserud,et al.  A basal area increment model for individual trees growing in even- and uneven-aged forest stands in Austria , 1996 .

[3]  Francis E. Putz,et al.  Critical need for new definitions of “forest” and “forest degradation” in global climate change agreements , 2009 .

[4]  Henrik Meilby,et al.  Site-specific height growth models for six common tree species in Denmark , 2009 .

[5]  Derek F. Sattler,et al.  Differences in crown characteristics between black (Picea mariana) and white spruce (Picea glauca) , 2012 .

[6]  Juha Lappi,et al.  A non-linear hierarchical mixed model to describe tree height growth , 2006, European Journal of Forest Research.

[7]  T. Pukkala,et al.  Predicting spatial distribution of direct radiation below forest canopies , 1991 .

[8]  S. Garman,et al.  Comparison of five canopy cover estimation techniques in the western Oregon Cascades , 2006 .

[9]  Wang Xinjie,et al.  Linear Mixed-Effects Models to Describe Individual Tree Crown Width for China-Fir in Fujian Province, Southeast China , 2015, PloS one.

[10]  Guillermo Trincado,et al.  A multilevel individual tree basal area increment model for aspen in boreal mixedwood stands , 2009 .

[11]  D. Hann,et al.  Crown profile equations for stand-grown western hemlock trees in northwestern Oregon , 2003 .

[12]  Margarida Tomé,et al.  A generalized nonlinear mixed-effects height-diameter model for Eucalyptus globulus L. in northwestern Spain. , 2010 .

[13]  Liyong Fu,et al.  A generalized nonlinear mixed-effects height to crown base model for Mongolian oak in northeast China , 2017 .

[14]  V. Grimm,et al.  Animal species diversity driven by habitat heterogeneity/diversity: the importance of keystone structures , 2004 .

[15]  John S. Kush,et al.  Crown and basal area relationships of open-grown southern pines for modeling competition and growth , 1992 .

[16]  Klaus von Gadow,et al.  A generalized height–diameter model including random components for radiata pine plantations in northwestern Spain , 2006 .

[17]  D. Bragg.  A local basal area adjustment for crown width prediction , 2001 .

[18]  R. Sharma,et al.  Individual tree crown width models for Norway spruce and European beech in Czech Republic , 2016 .

[19]  Hubert Hasenauer,et al.  A crown ratio model for Austrian forests , 1996 .

[20]  E. Vonesh,et al.  Linear and Nonlinear Models for the Analysis of Repeated Measurements , 1996 .

[21]  Y. Lei,et al.  Method of Estimating Degraded Forest Area: Cases from Dominant Tree Species from Guangdong and Tibet in China , 2020, Forests.

[22]  Shouzheng Tang,et al.  Comparison of seemingly unrelated regressions with error-in-variable models for developing a system of nonlinear additive biomass equations , 2015, Trees.

[23]  B. Parresol,et al.  Additivity in tree biomass components of Pyrenean oak (Quercus pyrenaica Willd.) , 2003 .

[24]  YangYuqing,et al.  Comparison of different methods for fitting nonlinear mixed forest models and for making predictions , 2011 .

[25]  Ram P. Sharma,et al.  Modelling crown width–diameter relationship for Scots pine in the central Europe , 2017, Trees.

[26]  Timothy G. Gregoire,et al.  Linear modelling of irregularly spaced, unbalanced, longitudinal data from permanent-plot measurements , 1995 .

[27]  Guangxing Wang,et al.  Modelling a system of nonlinear additive crown width models applying seemingly unrelated regression for Prince Rupprecht larch in northern China , 2017 .

[28]  Rafael Calama,et al.  Interregional nonlinear height-diameter model with random coefficients for stone pine in Spain , 2004 .

[29]  L. Fu,et al.  A generalized interregional nonlinear mixed-effects crown width model for Prince Rupprecht larch in northern China , 2017 .

[30]  Nonlinear Models for Repeated Measurement Data , 1996 .

[31]  Turan Sönmez.  Diameter at breast height-crown diameter prediction models for Picea orientalis. , 2009 .

[32]  J. Breidenbach,et al.  Modeling height-diameter relationships for Norway spruce, Scots pine, and downy birch using Norwegian national forest inventory data , 2015 .

[33]  Hubert Hasenauer,et al.  Biased predictions for tree height increment models developed from smoothed ‘data’ , 1997 .

[34]  Hailemariam Temesgen,et al.  Analysis and comparison of nonlinear tree height prediction strategies for Douglas-fir forests , 2008 .

[35]  Hans Pretzsch,et al.  The single tree-based stand simulator SILVA: construction, application and evaluation , 2002 .

[36]  D. Bates,et al.  Mixed-Effects Models in S and S-PLUS , 2001 .

[37]  Robert L. Bailey,et al.  Nonlinear Mixed Effects Modeling for Slash Pine Dominant Height Growth Following Intensive Silvicultural Treatments , 2001 .

[38]  LiYun,et al.  Evaluation of nonlinear equations for predicting diameter from tree height , 2012 .

[39]  Gregory S. Biging,et al.  Evaluation of Competition Indices in Individual Tree Growth Models , 1995, Forest Science.

[40]  William N. Venables,et al.  Modern Applied Statistics with S-Plus. , 1996 .

[41]  Rafael Calama,et al.  Multilevel linear mixed model for tree diameter increment in stone pine (Pinus pinea): a calibrating approach , 2005 .

[42]  Millenium Ecosystem Assessment.  Ecosystems and human well-being: synthesis , 2005 .

[43]  Yingjuan Su,et al.  Population genetic variation, differentiation and bottlenecks of Dacrydium pectinatum (Podocarpaceae) in Hainan Island, China: implications for its conservation , 2010 .

[44]  D. Bates,et al.  Nonlinear mixed effects models for repeated measures data. , 1990, Biometrics.

[45]  T. Pukkala,et al.  Relationship between radiation interception and photosynthesis in forest canopies: effect of stand structure and latitude , 1989 .

[46]  G. M. Bonnor.  A Comparison of Photo and Ground Measurements of Canopy Density , 1968 .

[47]  Christina L. Staudhammer,et al.  Individual Tree-Based Diameter Growth Model of Slash Pine in Florida Using Nonlinear Mixed Modeling , 2013 .

[48]  David L. R. Affleck,et al.  Additivity of nonlinear tree crown width models: Aggregated and disaggregated model structures using nonlinear simultaneous equations , 2018, Forest Ecology and Management.

[49]  E. Paoletti,et al.  Ozone Amplifies Water Loss from Mature Trees in the Short Term But Decreases It in the Long Term , 2019 .

[50]  Sun Yun-xiao.  Study on Biomass and Net Primary Productivity of Podocarpus imbricatus Plantation in Jianfengling,Hainan Island , 2004 .

[51]  Shouzheng Tang,et al.  Nonlinear mixed-effects crown width models for individual trees of Chinese fir (Cunninghamia lanceolata) in south-central China , 2013 .

[52]  Zhang Wei-yin.  The groups and features of tropical forest vegetation of Hainan Island , 2002 .

[53]  S. Meng,et al.  Improved calibration of nonlinear mixed-effects models demonstrated on a height growth function. , 2009 .

[54]  William A. Bechtold,et al.  Using crown condition variables as indicators of forest health , 2004 .