Nonlinear mixed-effects height to crown base and crown length dynamic models using the branch mortality technique for a Korean larch (Larix olgensis) plantations in northeast China

Abstract

Korean larch (Larix olgensis) is one of the main tree species for afforestation and timber production in northeast China. However, its timber quality and growth ability are largely influenced by crown size, structure and shape. The majority of crown models are static models based on tree size and stand characteristics from temporary sample plots, but crown dynamic models has seldom been constructed. Therefore, this study aimed to develop height to crown base (HCB) and crown length (CL) dynamic models using the branch mortality technique for a Korean larch plantation. The nonlinear mixed-effects model with random effects, variance functions and correlation structures, was used to build HCB and CL dynamic models. The data were obtained from 95 sample trees of 19 plots in Meng JiaGang forest farm in Northeast China. The results showed that HCB progressively increases as tree age, tree height growth (HT growth) and diameter at breast height growth (DBH growth). The CL was increased with tree age in 20 years ago, and subsequently stabilized. HT growth, DBH growth stand basal area (BAS) and crown competition factor (CCF) significantly influenced HCB and CL. The HCB was positively correlated with BAS, HT growth and DBH growth, but negatively correlated with CCF. The CL was positively correlated with BAS and CCF, but negatively correlated with DBH growth. Model fitting and validation confirmed that the mixed-effects model considering the stand and tree level random effects was accurate and reliable for predicting the HCB and CL dynamics. However, the models involving adding variance functions and time series correlation structure could not completely remove heterogeneity and autocorrelation, and the fitting precision of the models was reduced. Therefore, from the point of view of application, we should take care to avoid setting up over-complex models. The HCB and CL dynamic models in our study may also be incorporated into stand growth and yield model systems in China.

Keywords

Larix olgensis plantation / Height to crown base / Crown length / Branch mortality technique / Nonlinear mixed-effects models

References

[1] Assmann E. The principles of forest yield study: studies in the organic production, structure, increment and yield of forest stands, 1970, Oxford: Pergamon Press.

[2] Calama R, Montero G. Multilevel linear mixed-model for tree diameter increment in stone pine (Pinus pinea): a calibrating approach. Silva Fenn, 2005, 39(1): 37-54.

[3] Corral-Rivas S, Alvarez-Gonzalez JG, Crecente-Campo F, Corral-Rivas JJ. Local and generalized height-diameter models with random parameters for mixed, uneven-aged forests in Northwestern Durango, Mexico. For Ecosyst, 2014, 1(6): 1-9.

[4] Dong LB, Liu ZG, Bettinger P. nonlinear mixed-effects branch diameter and length models for natural Dahurian larch (Larix gmelini) forest in northeast China. Trees, 2016, 30: 1191-1206.

[5] Fahey TD, Cahill JM, Senllgrove TA, Heath LS. Lumber and veneer recovery from intensively managed young growth Douglas-fir US. For Serv Res, 2001, 437: 1-31.

[6] Fang ZX, Bailey RL. Nonlinear mixed effects modeling for slash pine dominant height growth following intensive silvicultural treatments. For Sci, 2001, 38: 2274-2286.

[7] Fu LY, Sun H, Sharma RP, Lei YC, Zhang HR, Tang SZ. Nonlinear mixed-effects crown width models for individual trees of Chinese fir (Cunninghamia lanceolata) in south-central China. For Ecol Manag, 2013, 302: 210-220.

[8] Fu LY, Zhang HR, Lu J, Zang H, Lou MH, Wang GX. Multilevel nonlinear mixed-effect crown ratio models for individual trees of Mongolian Oak (Quercus mongolica) in northeast China. PLoS ONE, 2014, 8: 1-20.

[9] Fu LY, Zhang HR, Sharma RP, Pang LF, Wang GX. A generalized nonlinear mixed-effects height to crown base model for Mongolian oak (Quercus mongolica) in northeast China. For Ecol Manag, 2017, 384: 34-43.

[10] Fujimori T. Dynamics of crown structure and stem growth based on knot analysis of a hinoki cypress (Chamaecyparis obtusa). For Ecol Manag, 1993, 56: 57-68.

[11] Garber SM, Monserud RA, Maguire DA. Crown recession patterns in three conifer species of the Northern Rocky Mountains. For Sci, 2008, 54: 633-646.

[12] Gonzalez-Benecke CA, Gezan SA, Samuelson LJ, Cropper WP. Estimating Pinus palustris tree diameter and stem volume from tree height, crown area and stand-level parameters. J For Res, 2014, 25: 43-52.

[13] Hann DW. Revised volume and taper equations for six major conifer species in Southwest Oregon, 2014, Corvallis: Oregon State University.

[14] Hann DW, Hanus ML (2002) Enhanced height-growth-rate equations for undamaged and damaged trees in Southwest Oregon. Oregon State University, Forest Research Laboratory, Research Contribution 41

[15] Hann DW, Hanus ML. Evaluation of nonspatial approaches and equation forms used to predict tree crown recession. Can J For Res, 2004, 34: 1993-2003.

[16] Hann DW, Walters DK, Scrivani JA. Incorporating crown ratio into prediction equations for Douglas-fir stem volume. Can J For Res, 1987, 17: 17-22.

[17] Hasenauer H, Monserud RA. A crown ratio model for Austrian forests. For Ecol Manag, 1996, 84(1–3): 49-60.

[18] Houllier F, Leban JM, Colin F. Linking growth modeling to timber quality assessment for Norway spruce. For Ecol Manag, 1995, 74: 91-102.

[19] Hynynen J. Modelling tree growth for managed stands, 1995, For Res Inst Res: Finnish 576p

[20] Jiang LC, Liu RL. Segmented taper equations with crown ratio and stand density for Dahurian Larch (Larix gmelinii) in Northeastern China. J For Res, 2011, 22: 347-352.

[21] Jozsa LA, Middleton GR (1994) A discussion of wood quality attributes and their practical implications. Forintek Canada Corp, Vancouver, B.C. Spec. Publ. SP

[22] Kraft G. Beiträge zur Lehre von den Durchforstungen, Schlagstellungen und Lichtungshieben, 1884, Hannover: Klindworth.

[23] Krajicek JE, Brinkman KA, Gingrich SF. Crown competition–a measure of density. For Sci, 1961, 7: 35-42.

[24] Kuprevicius A, Auty D, Achim A, Caspersen JP. Quantifying the influence of live crown ratio on the mechanical properties of clear wood. Forestry, 2013, 86(3): 361-369.

[25] Ledermann T. A non-linear model to predict crown recession of Norway spruce (Picea abies [L.] Karst.) in Austria. Eur J For Res, 2011, 130(4): 521-531.

[26] Leites LP, Robinson AP, Crookston NL. Accuracy and equivalence testing of crown ratio models and assessment of their impact on diameter growth and basal area increment predictions of two variants of the Forest Vegetation Simulator. Can J For Res, 2009, 39: 655-665.

[27] Li FR. Modeling crown profile of Larix olgensis tree. Sci Silvae Sin, 2004, 40(5): 16-24. (in Chinese with English abstract)

[28] Lindstrom JM, Bates DM. Nonlinear mixed effects models for repeated measures data. Biometrics, 1990, 46(3): 673-687.

[29] Liu ZG, Li FR. Crown structure model and three dimensional computer graphics simulation for Mongolian Scots pine plantation. Sci Silvae Sin, 2009, 6(45): 54-61. (in Chinese with English abstract)

[30] Lu J, Li FR, Zhang HR, Zhang SG. A crown ratio model for dominant species in secondary forests in Mao’er mountain. Sci Silvae Sin, 2011, 6(47): 70-76. (in Chinese with English abstract)

[31] Maguire DA, Hann DW. A stem dissection technique for dating branch mortality and reconstructing past crown recession. For Sci, 1987, 33: 858-871.

[32] Maguire DA, Hann DW. A sampling strategy for estimating past crown recession on temporary growth plots. For Sci, 1990, 36: 549-563.

[33] Maguire DA, Hann DW. Constructing models for direct prediction of 5-year crown recession in south-western Oregon Douglas-fir. Can J For Res, 1990, 20: 1044-1052.

[34] Mäkinen H. Effect of stand density on the branch development of silver birch (Betula pendula Roth) in central Finland. Trees, 2002, 16: 346-353.

[35] Meng SX, Huang S. Improved calibration of nonlinear mixed-effects models demonstrated on a height growth functions. For Sci, 2009, 55(3): 239-248.

[36] Mitchell KJ. Dynamics and simulated yield of Douglas-fir. For Sci Monogr, 1975, 17: 39-45.

[37] Monserud RA, 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(1–3): 57-80.

[38] Nguyen TT, Biber P, Pretzsch H. Analysis and management of stand dynamics of Vietnamese dipterocarp forests by applying a dynamic growth model. Ann For Sci, 2012, 69(5): 581-601.

[39] Pinheiro JC, Bates DM. Approximations to the log-likelihood function in the nonlinear mixed-effects model. J Comput Graph Stat, 1995, 4: 12-35.

[40] Pinheiro JC, Bates DM. Mixed-effects models in S and S-PLUS, 2000, New York: Springer 528p

[41] Popoola FS, Adesoye PO. Crown ratio models for Tectona ģrandis stands in Osho forest reserve, Oyo state, Nigeria. J For Sci, 2012, 28(2): 63-67.

[42] Pretzsch H. Modelling the crown competition of Norway spruce and Beech in pure and mixed stands. Allg Forst und Jagdztg, 1992, 163: 203-213.

[43] Pretzsch H, Biber P, Dursky J. The single tree based stand simulator SILVA. Construction, application and evaluation. For Ecol Manag, 2002, 162: 3-21.

[44] Rijal B, Weiskittel AR, Kershaw JA. Development of height to crown base models for thirteen tree species of the North American Acadian region. For Chron, 2012, 88: 60-73.

[45] Ritchie MW, Hann DW (1985) Equations for predicting basal area increment in Douglas-fir and Drand fir. Oregan State University, Forest Research Laboratory, Research Bulletin

[46] Ritchie MW, Hann DW. Equations for predicting height to crown base for fourteen tree species in southwest Oregon, 1987, Corvallis: Oregon State University, Forestry Research Laboratory.

[47] Ritchie MW, Zhang JW, Hamilton TA. Aboveground tree biomass for Pinus ponderosa in Northeastern California. Forests, 2013, 4: 179-196.

[48] Russell MB, Weiskittel AR. Maximum and largest crown width equations for fifteen tree species in Maine. North J Appl For, 2011, 28: 84-91.

[49] Russell MB, Weiskittel AR, Kershaw JAJ. Comparing strategies for modeling individual tree height and height-to-crown base increment in mixed-species Acadian forests of northeastern North America. Eur J For Res, 2014, 133: 1121-1135.

[50] Sattlerderek F, LeMay V. A system of nonlinear simultaneous equations for crown length and crown radius for the forest dynamics model SORTIE-ND. Can J For Res, 2011, 41(8): 1567-1576.

[51] Saud P, Lynch TB, Anup KC, Guldin JM. Using quadratic mean diameter and relative spacing index to enhance height-diameter and crown ratio models fitted to longitudinal data. Forestry, 2016, 89: 215-229.

[52] Seanm G, Robert M, Maguire DA. Crown recession patterns in three conifer species of the Northern rock mountains. For Sci, 2008 54 6 633

[53] Sharma RP, Bilek L, Vacek Z, Vacek S. Modelling crown width–diameter relationship for Scots pine in the central Europe. Trees, 2017, 31(6): 1-15.

[54] Soares P, Tomé M. A tree crown ratio prediction equation for eucalypt plantations. Ann For Sci, 2001, 58(2): 193-202.

[55] Spathelf P. Reconstruction of crown length of Norway spruce (Picea abies (L.) Karst.) and Silver fir (Abies alba Mill.)—technique, establishment of sample methods and application in forest growth analysis. Ann For Sci, 2003, 60: 833-842.

[56] Spiecker H. Growth variation and environmental stresses: long-term observations on permanent research plots in Southwestern Germany. Water Air Soil Poll, 1990, 54(2): 247-256.

[57] Stage AR (1973) Prognosis model for stand development. USDA Forest Service, Intermountain Forest and Range Experiment Station. Ogden, UT Research Paper INT-137

[58] Temesgen H, Lemay V, Mitchell SJ. Tree crown ratio models for multispecies and multi-layered stands of southeastern British Columbia. For Chron, 2005, 81(1): 133-141.

[59] Toney C, Reeves MC. Equations to convert compacted crown ratio to uncompacted crown ratio for trees in the interior west. West J Appl For, 2009, 24(2): 76-82.

[60] Valentine HT, Ludlow AR, Furnival GM. Modeling crown rise in even-aged stands of Stika spruce or loblolly pine. For Ecol Manag, 1994, 69: 189-197.

[61] Valentine HT, Makela A, Green EJ, Amateis RL, Mäkinen H, Ducey MJ. Models relating stem growth to crown length dynamics: application to loblolly pine and Norway spruce. Trees, 2012, 26: 469-478.

[62] Wang Y, LeMay VM, Baker TG. Modeling and prediction of dominant height and site index of Eucalyptus globulus plantation using a nonlinear mixed-effects model approach. Can J For Res, 2007, 37: 1390-1403.

[63] Weiskittel AR, Crookston NL, Radtke PJ. Linking climate, gross primary productivity, and site index across forests of the western United States. Can J For Res, 2011, 41: 1710-1721.

[64] Wensel LC, Daugherty PJ (1985) CACTOS user’s guide: the California conifer timber output simulator, version 2.0. California Forest Yield Cooperative, University of California, Department of Forest Resource Management, Berkeley, CA Res Nite 10

[65] Wykoff WR. A basal area increment model for individual conifers in the Northern Rocky Mountains. For Sci, 1990, 36(4): 1077-1104.

[66] Wykoff WR, Crookston NL, Stage AR (1982) Users’ guide to the stand prognosis model. US For Serv Gen Tech Rep. INT-GTR-133

[67] Young YQ, Huang SM. Estimating a multilevel dominant height- age model from nested data with generalized errors. For Sci, 2011, 57(2): 102-116.

[68] Zhao DH, Wilson M, Borders BE. Modeling response curves and testing treatment effects in repeated measures experiments: a multilevel nonlinear mixed-effects model approach. Can J For Res, 2005, 35: 122-132.

[69] Zumrawi AA, Hann DW (1989) Equations for predicting the height to crown base of six species in the Central Western Willamette Valley of Oregon. Forest Research Laboratory, Oregon State University, Research Paper 52