Generalized Nonlinear Mixed-Effects Individual Tree Diameter Increment Models for Beech Forests in Slovakia

Forests - Tập 10 Số 5 - Trang 451
Ram P. Sharma1, Igor Štefančík2, Zdeněk Vacek1, Stanislav Vacek1
1Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, Kamýcká 129, 16500 Prague 6, Czech Republic
2National Forest Centre, Forest Research Institute, T.G. Masaryka 22, SK-960 01 Zvolen, Slovakia

Tóm tắt

Individual tree growth and yield models precisely describe tree growth irrespective of stand complexity and are capable of simulating various silvicultural alternatives in the stands with diverse structure, species composition, and management history. We developed both age dependent and age independent diameter increment models using long-term research sample plot data collected from both monospecific and mixed stands of European beech (Fagus sylvatica L.) in the Slovak Republic. We used diameter at breast height (DBH) as a main predictor and other characteristics describing site quality (site index), stand development stage (dominant height and stand age), stand density or competition (ratio of individual tree DBH to quadratic mean diameter), species mixture (basal area proportion of a species of interest), and dummy variable describing stand management regimes as covariate predictors to develop the models. We evaluated eight versatile growth functions in the first stage using DBH as a single predictor and selected the most suitable one, i.e., Chapman-Richards function for further analysis through the inclusion of covariate predictors. We introduced the random components describing sample plot-level random effects and stochastic variations on the diameter increment, into the models through the mixed-effects modelling. The autocorrelation caused by hierarchical data-structure, which is assumed to be partially reduced by mixed-effects modelling, was removed through the inclusion of the parameter accounting for the autoregressive error-structures. The models described about two-third parts of a total variation in the diameter increment without significant trends in the residuals. Compared to the age independent mixed-effects model (conditional coefficient of determination, R c 2 = 0.6566; root mean square error, RMSE = 0.1196), the age dependent model described a significantly larger proportion of the variations in diameter increment ( R c 2 = 0.6796, RMSE = 0.1141). Diameter increment was significantly influenced differently by covariate predictors included into the models. Diameter increment decreased with the advancement of stand development stage (increased dominant height and stand age), increasing intraspecific competition (increased basal area proportion of European beech per sample plot), and diameter increment increased with increasing site quality (increased site index) and decreased competition (increased ratio of DBH to quadratic mean diameter). Our mixed-effects models, which can be easily localized with the random effects estimated from prior measurement of diameter increments of four randomly selected trees per sample plot, will provide high prediction accuracies. Our models may be used for simulating growth of European beech irrespective of its stand structural complexity, as these models have included various covariate variables describing both tree-and stand-level characteristics, thinning regimes, except the climate characteristics. Together with other forest models, our models will be used as inputs to the growth simulator to be developed in the future, which is important for decision-making in forestry.

Từ khóa


Tài liệu tham khảo

Garcia, 1994, The state-space approach in growth modeling, Can. J. For. Res., 24, 1894, 10.1139/x94-244

Amaro, A., Reed, D., and Soares, P. (2003). Modelling forest systems, CABI Publishing.

Vanclay, J.K. (1994). Modelling forest growth and yield: Applications to mixed tropical forests, CAB International.

Pretzsch, H. (2009). Forest dynamics, growth and yield: from measurement to model, Springer Verlag.

Weiskittel, A.R., Hann, D.W., Kershaw, J.A., and Vanclay, J.K. (2011). Forest growth and yield modeling, Wiley.

Hasenauer, H.E. (2006). Sustainable forest management: growth models for Europe, Springer-Verlag.

Adame, 2008, Individual-tree diameter growth model for rebollo oak (Quercus pyrenaica Willd.) coppices, For. Ecol. Manag., 255, 1011, 10.1016/j.foreco.2007.10.019

Soares, 2010, Modelling annual individual-tree growth and mortality of Scots pine with data obtained at irregular measurement intervals and containing missing observations, For. Ecol. Manag., 260, 1965, 10.1016/j.foreco.2010.08.044

Subedi, 2011, Individual-tree diameter growth models for black spruce and jack pine plantations in northern Ontario, For. Ecol. Manag., 261, 2140, 10.1016/j.foreco.2011.03.010

Vospernik, 2017, Possibilities and limitations of individual-tree growth models—A review on model evaluations, Die Bodenkultur: J. Land Manag., Food Environ., 68, 103, 10.1515/boku-2017-0010

Canellas, 2006, Distance independent tree diameter growth model for cork oak stands, For. Ecol. Manag., 225, 262, 10.1016/j.foreco.2006.01.002

Hasenauer, H. (2006). Concepts within tree growth modeling. Sustainable forest management: Growth models for Europe, Springer Verlag.

Sterba, 2008, Comparing an individual tree growth model for Pinus halepensis Mill. in the Spanish region of Murcia with yield tables gained from the same area, Eur. J. For. Res., 127, 253, 10.1007/s10342-007-0201-7

Sterba, 1997, Applicability of the forest stand growth simulator PROGNAUS for the Austrian part of the Bohemian Massif, Ecol. Model., 98, 23, 10.1016/S0304-3800(96)01934-5

Hasenauer, H. (2006). The tree growth model MOSES 3.0. Sustaianble forest management, Growth models for Europe, Springer Verlag.

Pinheiro, J.C., and Bates, D.M. (2000). Mixed-effects models in S and S-PLUS, Springer.

Saud, 2016, Using quadratic mean diameter and relative spacing index to enhance height-diameter and crown ratio models fitted to longitudinal data, Forestry, 89, 215, 10.1093/forestry/cpw004

Fu, 2017, A generalized interregional nonlinear mixed-effects crown width model for Prince Rupprecht larch in northern China, For. Ecol. Manag., 389, 364, 10.1016/j.foreco.2016.12.034

Hall, 2001, Modeling and prediction of forest growth variables based on multilevel nonlinear mixed models, For. Sci., 47, 311

Fu, 2013, Nonlinear mixed-effects crown width models for individual trees of Chinese fir (Cunninghamia lanceolata) in south-central China, For. Ecol. Manag., 302, 210, 10.1016/j.foreco.2013.03.036

Sharma, 2015, Modeling height-diameter relationships for Norway spruce, Scots pine, and downy birch using Norwegian national forest inventory data, For. Sci. Technol., 11, 44

Sharma, R.P., Vacek, Z., Vacek, S., Podrázský, V., and Jansa, V. (2017). Modelling individual tree height to crown base of Norway spruce (Picea abies (L.) Karst.) and European beech (Fagus sylvatica L.). PLoS ONE, 12.

Sharma, 2016, Individual tree crown width models for Norway spruce and European beech in Czech Republic, For. Ecol. Manag., 366, 208, 10.1016/j.foreco.2016.01.040

Bosela, 2016, The effects of climate warming on the growth of European beech forests depend critically on thinning strategy and site productivity, Agric. For. Meteorol., 222, 21, 10.1016/j.agrformet.2016.03.005

Keitel, 2007, Potential risks for European beech (Fagus sylvatica L.) in a changing climate, Trees, 21, 1

Knoke, 2008, Admixing broadleaved to coniferous tree species: a review on yield, ecological stability and economics, Eur. J. For. Res., 127, 89, 10.1007/s10342-007-0186-2

Barna, M., Ján, K., and Bublinec, E. (2011). Beech and Beech Ekosystems of Slovakia, VEDA.

(2019, May 23). Green report. Správa o lesnom hospodárstve v Slovenskej republike za rok 2016. Available online: https://www.enviroportal.sk/environmentalne-temy/vplyvy-na-zp/lesnictvo/dokumenty/spravy-o-lesnom-hospodarstve-v-slovenskej-republike.

Vacek, 2015, Changes in horizontal structure of natural beech forests on an altitudinal gradient in the Sudetes, Dendrobiology, 73, 33, 10.12657/denbio.073.004

Petritan, 2007, Effects of shade on growth and mortality of maple (Acer pseudoplatanus), ash (Fraxinus excelsior) and beech (Fagus sylvatica) saplings, Forestry, 80, 397, 10.1093/forestry/cpm030

Bolte, 2010, Climate change impacts on stand structure and competitive interactions in a southern Swedish spruce–beech forest, Eur. J. For. Res., 12, 261, 10.1007/s10342-009-0323-1

Bosela, 2015, Different mixtures of Norway spruce, silver fir, and European beech modify competitive interactions in central European mature mixed forests, Can. J. For. Res., 45, 1577, 10.1139/cjfr-2015-0219

Pretzsch, 2016, Mixing of Scots pine (Pinus sylvestris L.) and European beech (Fagus sylvatica L.) enhances structural heterogeneity, and the effect increases with water availability, For. Ecol. Manag., 373, 149, 10.1016/j.foreco.2016.04.043

Hanewinkel, 2013, Climate change may cause severe loss in the economic value of European forest land, Nat. Clim. Chang., 3, 203, 10.1038/nclimate1687

Boncina, 2007, Effects of selective thinning on growth and development of beech (Fagus sylvatica L.) forest stands in south-eastern Slovenia, Ann. For. Sci., 64, 47, 10.1051/forest:2006087

Vacek, 2018, Effect of thinning regimes on development and growth of crop trees in Fagus sylvatica stands of Central Europe over 50 years, Dendrobiology, 79, 141, 10.12657/denbio.079.013

Štefančík, I. (2019, May 23). Available online: http://sclib.svkk.sk/sck01/Record/000499076.

Sharma, 2011, Modelling dominant height growth from national forest inventory individual tree data with short time series and large age errors, For. Ecol. Manag., 262, 2162, 10.1016/j.foreco.2011.07.037

Hossfeld, J.W. (1822). Mathematik für Forstmänner, Őkonomen und Cameralisten, Nabu Press.

Cieszewski, 2002, Comparing fixed- and variable-base-age site equations having single versus multiple asymptotes, For. Sci., 48, 7

Zhao, 2012, Crown ratio and relative spacing relationships for loblolly pine plantations, Open J. For., 2, 101

Sharma, 2017, Modelling crown width-diameter relationship for Scots pine in the central Europe, Trees, 31, 1875, 10.1007/s00468-017-1593-8

Fonseca, 2017, A silvicultural stand density model to control understory in maritime pine stands, iForest, 10, 829, 10.3832/ifor2173-010

Schelhaas, 2018, Species-specific, pan-European diameter increment models based on data of 2.3 million trees, For. Ecosyst., 5, 21, 10.1186/s40663-018-0133-3

West, P.W. (2009). Tree and forest measurement, Springer.

Zeide, 1993, Analysis of growth equations, For. Sci., 39, 594

Zeide, 1989, Accuracy of equations describing diameter growth, Can. J. For. Res., 19, 1283, 10.1139/x89-195

Gyawali, 2015, Individual tree basal area growth models for Chir pine (Pinus roxberghii Sarg.) in western Nepal, J. For. Sci., 61, 535, 10.17221/51/2015-JFS

Barreiro, 2006, Age-independent difference equations for modelling tree and stand growth, Can. J. For. Res., 36, 1621, 10.1139/x06-065

Richards, 1959, A flexible growth function for empirical use, J. Exp. Bot., 10, 290, 10.1093/jxb/10.2.290

Neyman, J. (1961). Statistical Problems in Dynamics of Exploited Fisheries Populations. Proceedings of the Fourth Berkeley Symposium on Mathematical Statistics and Probability, University of California Press.

Bertalanffy, 1957, Quantitative laws in metabolism and growth, Quart. Rev. Biol., 32, 217, 10.1086/401873

Weibull, 1951, A statistical distribution function of wide applicability, J. Appl. Mech., 18, 293, 10.1115/1.4010337

Gompertz, 1825, On the nature of the function expressive of the law of human mortality and on a new model of determining life contingencies, Phil. Trans. R. Soc., 115, 513, 10.1098/rstl.1825.0026

Korf, 1939, A mathematical definition of stand volume growth law (In Czech), Lesnicka Prace, 18, 337

1936, Skogsforsö ksastaltens gallringsforsök i tallskog (Thinning experiments in pine forest conducted by the forest experiment station), Medd. fran Statens Skogsforsöksanstalt, 29, 1

Levakovic, 1935, An analytical form of growth law, Glasnik za Sumske Pokuse (In Serbo-Croat.), 4, 283

Fu, 2017, Modelling a system of nonlinear additive crown width models applying seemingly unrelated regression for Prince Rupprecht larch in northern China, For. Ecol. Manag., 386, 71, 10.1016/j.foreco.2016.11.038

Sharma, R.P., Vacek, Z., and Vacek, S. (2018). Generalized nonlinear mixed-effects individual tree crown ratio models for Norway spruce and European beech. Forests, 9.

Vonesh, E.F., and Chinchilli, V.M. (1997). Linear and nonlinear models for the analysis of repeated measurements, Marcel Dekker.

SAS Institute Inc. (2012). SAS/ETS1 9.1.3 User’s Guide, SAS Institute Inc.

Littell, R.C., Milliken, G.A., Stroup, W.W., Wolfinger, R.D., and Schabenberger, O. (2006). SAS for mixed models, SAS Institute. [2nd ed.].

Nakagawa, 2013, A general and simple method for obtaining R2 from generalized linear mixed-effects models, Methods Ecol. Evol., 4, 133, 10.1111/j.2041-210x.2012.00261.x

Akaike, 1974, A new look at statistical model identification, IEEE Trans. Autom. Control, 19, 716, 10.1109/TAC.1974.1100705

Calama, 2004, Interregional nonlinear height-diameter model with random coefficients for stone pine in Spain, Can. J. For. Res., 34, 150, 10.1139/x03-199

Soares, 2010, A generalized nonlinear mixed-effects height-diameter model for Eucalyptus globulus L. in northwestern Spain, For. Ecol. Manag., 259, 943, 10.1016/j.foreco.2009.11.036

Sharma, 2018, Modelling individual tree height-diameter relationships for multi-layered and multi-species forests in central Europe, Trees, 33, 103, 10.1007/s00468-018-1762-4

Carmean, 1989, Height growth and site index curves ofr jack pine in north central Ontario, Can. J. For. Res., 19, 215, 10.1139/x89-030

Goelz, 1992, Development of a well-behaved site index equation-Jack pine in North central Ontario, Can. J. For. Res., 22, 776, 10.1139/x92-106

Huang, 1994, An age-independent individual tree height prediction model for boreal spruce-aspen stands in Alberta, Can. J. For. Res., 24, 1295, 10.1139/x94-169

Burkhart, H.E., and Tomé, M. (2012). Modeling forest trees and stands, Springer.

Fang, 1998, Height-diameter models for tropical forests on Hainan Island in southern China, For. Ecol. Manag., 110, 315, 10.1016/S0378-1127(98)00297-7

Sharma, 2007, Height-diameter equations for boreal tree species in Ontario using a mixed-effects modeling approach, For. Ecol. Manag., 249, 187, 10.1016/j.foreco.2007.05.006

Monserud, 1996, A basal area increment model for individual trees growing in even- and uneven-aged forest stands in Austria, For. Ecol. Manag., 80, 57, 10.1016/0378-1127(95)03638-5

Cienciala, 2016, Discerning environmental factors affecting current tree growth in Central Europe, Sci. Total Environ., 573, 541, 10.1016/j.scitotenv.2016.08.115

Zhao, 2013, Individual-tree diameter growth model for fir plantations based on multi-level linear mixed effects models across southeast China, J. For. Res., 18, 305, 10.1007/s10310-012-0352-3

Dormann, 2013, Collinearity: A review of methods to deal with it and a simulation study evaluating their performance, Ecography, 36, 27, 10.1111/j.1600-0587.2012.07348.x

2002, Radial, height and volume increment variation in Picea abies (L.) Karst. Stands with varying thinning intensities, Scand. J. For. Res., 17, 304, 10.1080/02827580260138062

Sharma, 2016, Modeling individual tree height to diameter ratio for Norway spruce and European beech in Czech Republic, Trees, 30, 1969, 10.1007/s00468-016-1425-2

Wonn, 2001, Height: Diameter ratios and stability relationships for four northern rocky mountain tree species, West. J. Appl. For., 16, 87, 10.1093/wjaf/16.2.87

Kim, 2016, Impact of thinning intensity on the diameter and height growth of Larix kaempferi stands in central Korea, For. Sci. Technol., 12, 77

Binkley, 2004, Thinking about efficiency of resource use in forests, For. Ecol. Manag., 193, 5, 10.1016/j.foreco.2004.01.019

Bayer, 2013, Structural crown properties of Norway spruce (Picea abies [L.] Karst.) and European beech (Fagus sylvatica [L.]) in mixed versus pure stands revealed by terrestrial laser scanning, Trees, 27, 1035, 10.1007/s00468-013-0854-4

Pretzsch, 2014, Canopy space filling and tree crown morphology in mixed-species stands compared with monocultures, For. Ecol. Manag., 327, 251, 10.1016/j.foreco.2014.04.027

Pretzsch, 2010, Comparison between the productivity of pure and mixed stands of Norway spruce and European beech along an ecological gradient, Ann. For. Sci., 67, 712, 10.1051/forest/2010037

Sterba, 2014, Effect of species proportion definition on the evaluation of growth in pure vs. mixed stands, For. Syst., 23, 547, 10.5424/fs/2014233-06051

Pretzsch, 2015, Representation of species mixing in forest growth models. A review and perspective, Ecol. Model., 313, 276, 10.1016/j.ecolmodel.2015.06.044

Sharma, 2017, Modelling tree crown-to-bole diameter ratio for Norway spruce and European beech, Silva Fenn., 51, 1740, 10.14214/sf.1740

Kozak, 2003, Does cross validation provide additional information in the evaluation of regression models?, Can. J. For. Res., 33, 976, 10.1139/x03-022

Yang, 2004, An evaluation of diagnostic tests and their roles in validating forest biometric models, Can. J. For. Res., 34, 619, 10.1139/x03-230