Individual-tree diameter growth model for Korean pine plantations based on optimized interpolation of meteorological variables

Abstract

To explore the influence of meteorological variables on the growth of Korean pine (Pinus koraiensis Sieb. et Zucc.) plantations and provide a scientific reference for the production and management of Korean pine, three approaches to interpolate meteorological variables during the growing season (i.e., May–September) were compared in Heilongjiang Province, China. Optimized meteorological variable interpolation results were then combined with stand and individual tree variables, based on data from 56 sample plots and 2886 sample trees from Korean pine plantations in two regions of the province to develop an individual-tree diameter growth model (Model I) and an individual-tree diameter growth model with meteorological variables (Model II) using a stepwise regression method. Moreover, an individual-tree diameter growth model with regional effects (Model III) was developed using dummy variables in the regression, and the significance of introducing these dummy variables was verified with an F-test statistical analysis. The models were validated using an independent data set, and the predictive performance of the three models was assessed via the adjusted coefficient of determination (

R_a^2

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${R}_{a}^{2}$$\end{document}

) and root mean square error (RMSE). The results suggest that the growth increment in tree diameter of Korean pine plantations was significantly correlated with the natural logarithm of initial diameter (ln D), stand basal area (BAS), logarithmic deformation of the stand density index (ln SDI), ratio of basal area of trees larger than the subject tree to their initial diameter at breast height (DBH) (BAL/D), and the maximum growing-season precipitation (Pgmax). The individual-tree diameter growth models of Korean pine plantations developed in this study will provide a good basis for estimating and predicting growth increments of Korean pine forests over larger areas.

Keywords

Mixed interpolation / Korean pine plantations / Individual-tree diameter growth model / Regional effects

References

[1] Aber J, Neilson RP, McNulty S, Lenihan JM, Bachelet D, Drapek RJ. Forest processes and global environmental change: predicting the effects of individual and multiple stressors. Bioscience, 2001, 51(9): 735-751.

[2] Adame P, Hynynen J, Canellas I, del Rio M. Individual-tree diameter growth model for rebollo oak (Quercus pyrenaica Willd.) coppices. For Ecol Manag, 2008, 255(3–4): 1-1022.

[3] Brooks JR, Flanagan LB, Ehleringer JR. Responses of boreal conifers to climate fluctuations: indications from tree-ring widths and carbon isotope analyses. Can J For Res, 1998, 28: 524-533.

[4] Burkhart HE. Reed D, Soares P. Suggestion for choosing an appropriate level for modelling forest stands. Modelling forest systems, 2003, Wallingford: CAB International 3 10

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

[6] Campos JCC, Leite HG. Mensura florestal: perguntas e respostas, 2006, Editora: UFV 636

[7] Cawley GC, Talbot NL. Fast exact leave-one-out cross-validation of sparse least-squares support vector machines. Neural Netw, 2004, 17(10): 1467-1475.

[8] Chen DS, Huang XZ, Zhang SG, Sun XM. Biomass modeling of larch (Larix spp.) plantations in China based on the mixed model, dummy variable model, and Bayesian hierarchical model. Forests, 2017 8 8 268

[9] Crookston NL, Rehfeldt GE, Dixon GE, Weiskittel AR. Addressing climate change in the forest vegetation simulator to assess impacts on landscape forest dynamics. For Ecol Manag, 2010, 260(7): 1198-1211.

[10] Daly C. Guidelines for assessing the suitability of spatial climate data sets. Int J Climatol, 2006, 6(26): 707-721.

[11] Daniels LD, Veblen TT. Spatiotemporal influences of climate on altitudinal treeline in northern Patagonia. Ecology, 2004, 85(5): 1284-1296.

[12] Delrieu G, Wijbrans A, Boudevillain B, Faure D, Bonnifait L, Kirstetter PE. Geostatistical radar-raingauge merging: a novel method for the quantification of rain estimation accuracy. Adv Water Resour, 2014, 71: 110-124.

[13] Deslauriers A, Morin H, Begin Y. Cellular phenology of annual ring formation of Abies balsamea in the Quebec boreal forest (Canada). Can J For Res, 2003, 33(2): 190-200.

[14] Fabbio G, Cantiani P, Ferretti F, di Salvatore U, Bertini G, Becagli C, Chiavetta U, Marchi M, Salvati L. Sustainable land management, adaptive silviculture, and new forest challenges: evidence from a latitudinal gradient in Italy. Sustainability, 2018 10 7 2520

[15] Fox JC, Ades PK, Bi HQ. Stochastic structure and individual-tree growth models. For Ecol Manag, 2001, 154(1–2): 261-276.

[16] Genton MG, Furrer R. Analysis of rainfall data by robust spatial statistic using S+SPATIALSTATS. J Geogr Inf Decis Anal, 1998, 2(2): 116-126.

[17] Getis A, Ord JK. The analysis of spatial association by use of distance statistics. Geogr Anal, 1992, 24(3): 189-206.

[18] Gordon N, Will S. Effect of climate on lodgepole pine stem taper in British Columbia, Canada. Forestry, 2012, 85(5): 579-587.

[19] Hamann A, Wang TL. Potential effects of climate change on ecosystem and tree species distribution in British Columbia. Ecology, 2006, 87(11): 2773-2786.

[20] Hay L, Viger R, McCabe G. Precipitation interpolation in mountainous regions using multiple linear regression. Hydrol Water Resourc Ecol Headwaters, 1998, 248: 33.

[21] Jin XJ, Pukkala T, Li FR, Dong LH. Optimal management of Korean pine plantations in multifunctional forestry. J For Res, 2017, 28(5): 1027-1037.

[22] Joly D, Brossard T, Cardot H, Cavailhes J, Hilal M, Wavresky P. Temperature interpolation by local information; the example of France. Int J Climatol, 2011, 31(14): 2141-2153.

[23] Lawrence S, Davis J. Forest management: to sustain ecological, economic, and social values, 2005, Long Grove: Waveland Press 816

[24] Lhotka JM, Loewenstein EF. An individual-tree diameter growth model for managed uneven-aged oak-shortleaf pine stands in the Ozark Highlands of Missouri, USA. For Ecol Manag, 2011, 261(3): 770-778.

[25] Li CM. Individual tree diameter increment model for Chinese fir plantation based on two-level linear mixed effects models. Scientia Silvae Sinicae, 2012, 48(3): 66-73. (in Chinese)

[26] Li XF, Hong XL. Research on the use of dummy variables method to calculate the family of site index curves. For Res, 2012, 10(2): 215-219. (in Chinese)

[27] Li X, Liu XT, Wei JT, Li Y, Tigabu M, Zhao XY. Genetic improvement of Pinus koraiensis in China: current situation and future prospects. Forests, 2020 11 2 148

[28] Littell JS, McKenzie D, Kerns BK, Cushman S, Shaw CG. Managing uncertainty in climate-driven ecological models to inform adaptation to climate change. Ecosphere, 2011 2 9 102

[29] Liu Q, Li FR, Dong LH. Individual tree growth model of larch plantation based on crown competition factors. Bull Bot Res, 2014, 34(4): 547-553.

[30] Ma W, Lei XD, Xu G, Yang YJ, Wang QJ. Growth models for natural Quercus mongolica forest-I. Diameter growth model. J Northwest A F Univ (Nat Sci Educ), 2015, 2: 16. (in Chinese)

[31] Mabvurira D, Miina J. Individual-tree growth and mortality models for Eucalyptus grandis (Hill) Maiden plantations in Zimbabwe. For Ecol Manag, 2002, 161(1): 231-245.

[32] McDowell NG, Beerling DJ, Breshears DD, Fisher RA, Raffa KF, Stitt M. The interdependence of mechanisms underlying climate-driven vegetation mortality. Trends Ecol Evol, 2011, 26(10): 523-532.

[33] Meng SX, Huang S, Lieffers VJ, Nunifu T, Yang Y. Wind speed and crown class influence the height-diameter relationship of lodgepole pine: nonlinear mixed effects modeling. For Ecol Manag, 2008, 256(4): 570-577.

[34] Miles J. Everitt BS, Howell DC. Tolerance and Variance Inflation Factor. Encyclopedia of statistics in behavioral science, 2005, Chichester: Wiley 2055 2056

[35] Mitchell KJ. SYLVER: modelling the impact of silviculture on yield, lumber value, and economic return. For Chronicle, 1988, 64(2): 127-131.

[36] Nalder IA, Wein RW. Spatial interpolation of climatic normals: test of a new method in the Canadian boreal forest. Agric For Meteorol, 1998, 92(4): 211-225.

[37] Parresol BR. Assessing tree and stand biomass: a review with examples and critical comparisons. For Sci, 1999, 45(4): 573-593.

[38] State Forestry Administration (2014) The eighth forest resource survey report. https://211.167.243.162:8085/8/index.html

[39] Subedi N, Sharma M. Individual-tree diameter growth models for black spruce and jack pine plantations in northern Ontario. For Ecol Manag, 2011, 261(11): 2140-2148.

[40] Uzoh FC, Oliver WW. Individual tree diameter increment model for managed even-aged stands of ponderosa pine throughout the western United States using a multilevel linear mixed effects model. For Ecol Manag, 2008, 256(3): 438-445.

[41] Vanclay JK (1988) PLATIPUS physiology: Design of a plantation growth and conversion simulator incorporating silviculture and wood quality. In: Ek AR, Shifley SR, Burk TE (eds) Forest growth modelling and prediction. Minneapolis, pp 998–1005

[42] Vanclay JK. Modelling forest growth and yield: applications to mixed tropical forests, 1994, Wallingford: CAB International 16 17

[43] Vicente-Serrano SM, Saz-Sánchez MA, Cuadrat JM. Comparative analysis of interpolation methods in the middle Ebro Valley (Spain): application to annual precipitation and temperature. Clim Res, 2003, 24(2): 161-180.

[44] Worbes M. Annual growth rings, rainfall-dependent growth and long-term growth patterns of tropical trees from the Caparo Forest Reserve in Venezuela. J Ecol, 1999, 87(3): 391-403.

[45] Wu XW, Yan LM. Setting parameters and choosing optimum semivariogram models of ordinary Kriging interpolation—a case study of spatial interpolation to January average temperature of Fujian province. Geo-Inf Sci, 2007, 9(3): 104-108. (in Chinese)

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

[47] Yang YQ, Titus SJ, Huang SM. Modeling individual tree mortality for white spruce in Alberta. Ecol Model, 2003, 163(3): 209-222.

[48] Yue WZ, Xu JH, Xu LH. A study on spatial interpolation methods for climate variables based on geostatistics. Plateau Meteorol, 2005, 24(6): 974-980. (in Chinese)

[49] Yun HS, Um MJ, Cho WC, Heo JH. Orographic precipitation analysis with regional frequency analysis and multiple linear regression. J Korea Water Resour Assoc, 2009, 42(6): 465-480.

[50] Zeng WS, Duo HR, Lei XD, Chen XY, Wang XJ, Pu Y, Zou WT. Individual tree biomass equations and growth models sensitive to climate variables for Larix spp. in China. Eur J For Res, 2017, 136(2): 233-249.

[51] Zhang HP, Li FR, Dong LH, Liu Q. Individual tree diameter increment model for natural Betula platyphylla forests based on meteorological factors. Chin J Appl Ecol, 2017, 6: 1851-1859. (in Chinese)

[52] Zhao DH, Borders B, Wilson M. Individual-tree diameter growth and mortality models for bottomland mixed-species hardwood stands in the lower Mississippi alluvial valley. For Ecol Manag, 2004, 199(2–3): 307-322.

[53] Zhao LF, Li CM, Tang SZ. Individual-tree diameter growth model for fir plantations based on multi-level linear mixed effects models across southeast China. J For Res, 2013, 18(4): 305-315.

[54] Zhao YH, Lyu HC, Zhen Z, Li FR, Wei QB. Interpolation optimization of meteorological factors and its correlation analysis with the larch NPP in Heilongjiang Province, China. J Nanjing For Univ (Nat Sci Ed), 2018, 42(03): 1-9. (in Chinese)

[55] Zhen Z, Li FR, Liu ZG, Liu C, Zhao YL, Ma ZH, Zhang LJ. Geographically local modeling of occurrence, count, and volume of downwood in Northeast China. Appl Geogr, 2013, 37(4): 114-126.