Incorporating topographic factors in nonlinear mixed-effects models for aboveground biomass of natural Simao pine in Yunnan, China

Abstract

A total of 128 Simao pine trees (Pinus kesiya var. langbianensis) from three regions of Pu’er City, Yunnan Province, People’s Republic of China, were destructively sampled to obtain tree aboveground biomass (AGB). Tree variables such as diameter at breast height and total height, and topographical factors such as altitude, aspect of slope, and degree of slope were recorded. We considered the region and site quality classes as the random-effects, and the topographic variables as the fixed-effects. We fitted a total of eight models as follows: least-squares nonlinear models (BM), least-squares nonlinear models with the topographic factors (BMT), nonlinear mixed-effects models with region as single random-effects (NLME-RE), nonlinear mixed-effects models with site as single random-effects (NLME-SE), nonlinear mixed-effects models with the two-level nested region and site random-effects (TLNLME), NLME-RE with the fixed-effects of topographic factors (NLMET-RE), NLME-SE with the fixed-effects of topographic factors (NLMET-SE), and TLNLME with the fixed-effects of topographic factors (TLNLMET). The eight models were compared by model fitting and prediction statistics. The results showed: model fitting was improved by considering random-effects of region or site, or both. The models with the fixed-effects of topographic factors had better model fitting. According to AIC and BIC, the model fitting was ranked as TLNLME > NLMET-RE > NLME-RE > NLMET-SE > TLNLMET > NLME-SE > BMT > BM. The differences among these models for model prediction were small. The model prediction was ranked as TLNLME > NLME-RE > NLME-SE > NLMET-RE > NLMET-SE > TLNLMET > BMT > BM. However, all eight models had relatively high prediction precision (>90 %). Thus, the best model should be chosen based on the available data when using the model to predict individual tree AGB.

Keywords

Aboveground biomass / Mixed-effects models / Regional effect / Site quality effect / Topographic factors / Pinus kesiya var. langbianensis

References

[1] Abrams MD, Knapp AK, Hulbert LC. A ten-year record of aboveground biomass in a Kansas tallgrass prairie: effects of fire and topographic position. Am J Bot, 1986, 73(10): 1509-1515.

[2] Alves LF, Vieira SA, Scaranello MA, Camargo PB, Santos FAM, Joly CA, Martinelli LA. Forest structure and live aboveground biomass variation along an elevational gradient of tropical Atlantic moist forest (Brazil). For Ecol Manage, 2010, 260(5): 679-691.

[3] Asner GP, Hughes RF, Varga T, Knapp D, Kennedy-Bowdoin T. Environmental and biotic controls over aboveground biomass throughout a tropical rain forest. Ecosystems, 2009, 12(2): 261-278.

[4] Basuki T, van Laake PE, Skidmore AK, Hussin YA. Allometric equations for estimating the above-ground biomass in tropical lowland Dipterocarp forests. For Ecol Manage, 2009, 257(8): 1684-1694.

[5] Bennie J, Hill MO, Baxter R, Huntley B. Influence of slope and aspect on long-term vegetation change in British chalk grasslands. J Ecol, 2006, 94(2): 355-368.

[6] Blackard JA, Finco MV, Helmer EH, Holden GR, Hoppus ML, Jacobs DM, Lister AJ, Moisen GG, Nelson MD, Riemann R, Ruefenacht B. Mapping U.S. forest biomass using nationwide forest inventory data and moderate resolution information. Remote Sens Environ, 2008, 112: 1658-1677.

[7] Brown S, Lugo AE. Biomass of tropical forests: a new estimate based on forest volumes. Science, 1984, 223: 1290-1293.

[8] Brown IF, Martinelli LA, Thomas WW, Moreira MZ, Ferreira CA, Victoria RA. Uncertainty in the biomass of Amazonian forests: an example from Rondonia, Brazil. For Ecol Manage, 1995, 75: 175-189.

[9] Budhathoki CB, Lyneh TB, Guldin JM (2005) Individual tree growth models for natural even-aged Shortleaf pine. Gen. Tech. Rep. SRS-92. U.S. Department of Agriculture, Forest Service, Southern Research Station. Asheville, pp 359–361

[10] 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.

[11] Calegario N, Maestri R, Leal CL, Daniels RF. Growth estimate of Eucalyptus stands based on nonlinear multilevel mixed-effects model theory. Ciência Florestal, 2005, 15: 285-292.

[12] Case B, Hall RJ. Assessing prediction errors of generalized tree biomass and volume equations for boreal forest region of west-central Canada. Can J For Res, 2008, 38: 878-889.

[13] Chapin FS III, Bloom AJ, Field CB, Waring RH. Plant responses to multiple environmental factors. Bioscience, 1987, 37(1): 49-57.

[14] Chapin FS III, Matson PA, Mooney HA. Principles of Terrestrial Ecosystem Ecology. 2002, New York: Springer, 529.

[15] Chen DS, Li FR, Sun XM, Jia WW. Models to predict knot size for Larch plantation using linear mixed model. Sci Silvae Sin, 2011, 47(11): 121-128.

[16] Chojnacky DC. Allometric scaling theory applied to FIA biomass estimation. Proceeding of the third annual forest inventory and analysis symposium, GTR NC-230. 2002, Minot: North Central Research Station, Forest Service USDA, 96 102

[17] Compilation Committee of Yunnan Forest Yunnan Forest. 1986, Kunming, Beijing: Yunnan Science and Technology Press, China Forestry Publishing House

[18] de Castilho CV, Magnusson WE, de Araujo RNO, Luizão RCC, Luizão FJ, Lima AP, Higuchi N. Variation in aboveground tree live biomass in a central Amazonian Forest: effects of soil and topography. For Ecol Manage, 2006, 234(1/3): 85-96.

[19] DeWalt SJ, Chave J. Structure and biomass of four lowland Neotropical forests. Biotropica, 2004, 36(1): 7-19.

[20] Dorado FC, Diéguez-Aranda U, Anta MB, Rodríguez MS, von Gadow K. A generalized height-diameter model including random components for Radiata pine plantations in Northwestern Spain. For Ecol Manage, 2006, 229: 202-213.

[21] Fahey TJ, Knapp AK. Principles and standards for measuring primary production. 2007, New York: Oxford University Press, 288

[22] Fang Z, Bailey RL. Nonlinear mixed effects modeling for slash pine dominant height growth following intensive silvicultural treatments. For Sci, 2001, 47: 287-300.

[23] Fehrmann L, Lehtonen A, Kleinn C, Tomppo E. Comparison of linear and mixed-effect regression models and a K-nearest neighbour approach for estimation of single-tree biomass. Can J For Res, 2008, 38(1): 1-9.

[24] Fu LY (2012) Nonlinear mixed effects model and its application in forestry. Doctor of Thesis, Chinese Academy of Forestry, Beijing, China

[25] Fu LY, Sun H. Individual crown diameter prediction for Cunninghamia lanceolata forests based mixed effects models. Sci Silvae Sin, 2013, 49(8): 65-74.

[26] Fu LY, Zeng WS, Tang SZ, Sharma RP, Li HK. Using linear mixed model and dummy variable model approaches to construct compatible single-tree biomass equations at different scales –A case study for Masson pine in Southern China. J For Sci, 2012, 58(3): 101-115.

[27] Fu L, Zeng W, Zhang H, Wang G, Lei Y, Tang S. Generic linear mixed-effects individual-tree biomass models for Pinus massoniana in southern China. Southern Forests, 2014, 76(1): 47-56.

[28] Garber SM, Maguire DA. Modeling stem taper of three Central Oregon species using nonlinear mixed effects models and autoregressive error structures. For Ecol Manage, 2003, 179: 507-522.

[29] Gregoire TG, Schabenberger O. A nonlinear mixed-effects model to predict cumulative bole volume of standing trees. J Appl Stat, 1996, 23: 257-271.

[30] Gregoire TG, Schabenberger O, Barrett JP. Linear modelling of irregularly spaced, unbalanced, longitudinal data from permanent- plot measurements. Can J For Res, 1995, 25: 137-156.

[31] Guilley E, Herve J, Nepveu G. The influence of site quality, silviculture and region on wood density mixed model in Quercus petraea Liebl. For Ecol Manage, 2004, 189: 111-121.

[32] Hall DB, Clutter M. Multivariate multilevel nonlinear mixed effects models for timber yield predictions. Biometrics, 2004, 60: 16-24.

[33] Henry M, Picard N, Trotta C, Manlay RJ, Valentini R, Bernoux M, Saint-André L. Estimating tree biomass of sub-Saharan African forests: a review of available allometric equations. Silva Fenn, 2011, 45: 477-569.

[34] Jenkins JC, Chojnacky DC, Heath LS, Bird-se RA. National-scale biomass estimators for United States tree species. For Sci, 2003, 49: 12-35.

[35] Jiang LC, Li FR. Application of mixed effects models on forestry modelling. 2014, Beijing: Science Press, 143.

[36] Kang M, Dai C, Ji W, Jiang Y, Yuan Z, Chen HYH. Biomass and its allocation in relation to temperature, precipitation, and soil nutrients in Inner Mongolia grasslands, China. PLoS One, 2013, 8 7 e69561

[37] Laird NM, Ware JH. Random effects models for longitudinal data. Biometrics, 1982, 38: 963-974.

[38] Lappi J, Bailey RL. A height prediction model with random stand and tree parameters: an alternative to traditional site methods. For Sci, 1988, 34: 907-927.

[39] Laumonier Y, Edin A, Kanninen M, Munandar AW. Landscape-scale variation in the structure and biomass of the hill Dipterocarp forest of Sumatra: implications for carbon stock assessments. For Ecol Manage, 2010, 259: 505-513.

[40] Laurance WF, Fearnside PM, Laurance SG, Delamonica P, Lovejoy TE, de Merona RJM, Chambers JQ, Gaston C. Relationship between soils and Amazon forest biomass: a landscape-scale study. For Ecol Manage, 1999, 118: 127-138.

[41] Lei XD, Li YC, Xiang W. Individual basal area growth model using multi-level linear mixed model with repeated measures. Sci Silvae Sin, 2009, 45(1): 74-80.

[42] Leites LP, Robinson AP. Improving taper equations of loblolly pine with crown dimensions in a mixed-effects modeling framework. For Sci, 2004, 50: 204-212.

[43] Li CM. Application of mixed effects models in forest growth model. Sci Silvae Sin, 2009, 45(4): 131-138.

[44] Li CM (2010) Application of mixed effects models in forest growth models. Doctor of Thesis, Chinese Academy of Forestry, Beijing, China

[45] Li J (2010) Dynamics of biomass and carbon stock for young and middle aged plantation of Simao pine (Pinus kesiya var. langbianensis). Doctoral Dissertation, Beijing Forestry University, Beijing, China

[46] Li CM. Study on dominant height growth of fir plantation based on a nonlinear mixed modeling approach for longitudinal data. For Res, 2011, 24(1): 68-73.

[47] Li CM. The simultaneous equation system of total volume in fir plantation. Sci Silvae Sin, 2012, 48(6): 80-88.

[48] Li YX, Jiang LC. Modeling wood density with two-level linear mixed effects models for Dahurian larch. Sci Silvae Sin, 2013, 49(7): 91-98.

[49] Lieth H, Whittaker RH. Primary productivity of the biosphere. 1975, New York: Springer, 237 263

[50] López IF, Lambert MG, Mackay AD, Valentine I. The influence of topography and pasture management on soil characteristics and herbage accumulation in hill pasture in the North Island of New Zealand. Plant Soil, 2003, 255(2): 421-434.

[51] Mahadev S, John P. Height-diameter equations for boreal tree species in Ontario using a mixed-effects modeling approach. For Ecol Manage, 2007, 249: 187-198.

[52] Malhi Y, Wood D, Bakers TR, Wright J, Phillips O, Cochrane T, Neir P, Chave J, Almeida S, Arroyo L, Higuchi N, Killeen TJ, Laurance SG, Laurance WF, Lewis SL, Monteagudo A, Neill DA, Vargas PN, Pitman NCA, Quesada CA, Salomão R, Silva JNM, Lezama AT, Terborgh J, Matínez RV, Vinceti B. The regional variation of aboveground live biomass in old-growth Amazonian forests. Global change Biol, 2006, 12(7): 1107-1138.

[53] Matthew BR, Harold EB, Ralph LA. Biomass partitioning in a miniature-scale loblolly pine spacing trial. Can J For Res, 2009, 39: 320-329.

[54] McEwan RW, Lin YC, Sun IF, Hsieh CF, Su SH, Chang LW, Song GZM, Wang HH, Hwong JL, Lin KC, Yang KC, Chiang JM. Topographic and biotic regulation of aboveground carbon storage in subtropical broad-leaved forests of Taiwan. For Ecol Manage, 2011, 262: 1817-1825.

[55] Muukkonen P. Generalized allometric volume and biomass equations for some tree species in Europe. Eur J For Res, 2007, 126(2): 157-166.

[56] Nanos N, Calama R, Montero G, Gil L. Geostatistical prediction of height/diameter models. For Ecol Manage, 2004, 195: 221-235.

[57] Návar J. Allometric equations for tree species and carbon stocks for forests of northwestern Mexico. For Ecol Manage, 2009, 257: 427-433.

[58] Paoli GD, Curran LM, Slik JWF. Soil nutrients affect spatial patterns of aboveground biomass and emergent tree density in southwestern Borneo. Oecologia, 2008, 155(2): 287-299.

[59] Pearce HG, Anderson WR, Fogarty LG, Todoroki CL, Anderson SAJ. Linear mixed-effects models for estimation biomass and fuel loads in shrublands. Can J For Res, 2010, 40(10): 2015-2026.

[60] Phillips OL, Malhi Y, Higuchi N, Laurance WF, Núnez PV, Vásquez RM, Laurance SG, Ferreira LV, Stern M, Brown S, Grace J. Changes in the carbon balance of tropical forests: evidence from long-term plots. Science, 1998, 282: 439-442.

[61] Picard N, Saint-André L, Henry M (2012) Manual for building tree volume and biomass allometric equations: from field measurement to prediction. Food and Agricultural Organization of the United Nations, Rome, and Centre de Coopération Internationale en Recherche Agronomique pour le Développement, Montpellier

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

[63] Polis GA. Why are parts of the world green? Multiple factors control productivity and the distribution of biomass. Oikos, 1999, 86(1): 3-15.

[64] Popescu SC. Estimating biomass of individual pine trees using airborne lidar. Biomass Bioenergy, 2007, 31: 646-655.

[65] Popescu SC, Wynne RH, Nelson RF. Measuring individual tree crown diameter with lidar and assessing its influence on estimating forest volume and biomass. Can. J. Remote Sens, 2003, 29: 564-577.

[66] Ramananantoandro T, Rafidimanantsoa HP, Ramanakoto MF. Forest aboveground biomass estimates in a tropical rainforest in Madagascar: new insights from the use of wood specific gravity data. J For Res, 2015, 26(1): 47-55.

[67] Rose R, Rosner LS, Ketchum JS. Twelfth-year response of Douglas-fir to area of weed control and herbaceous versus woody weed control treatments. Can J For Res, 2006, 36(10): 2464-2473.

[68] Shen ZH, Fang JY. Niche comparison of two Fagus species based on the topographic pattern of their populations. Acta Phytoecologica Sinica, 2001, 25: 392-398.

[69] Shen ZH, Zhang XS, Jin YX. Gradient analysis of the influence of mountain topography on vegetation pattern. Acta Phytoecologica Sinica, 2000, 24: 430-435.

[70] Simard M, Zhang K, Rivera-Monroy VH, Ross MS, Ruiz PL, Castañeda-Moya E, Twilley RR, Rodriguez E. Mapping height and biomass of mangrove forests in Everglades national park with SRTM elevation data. Photogramm. Eng Rem Sens, 2006, 72: 299-311.

[71] Southwest Forestry College, Forestry Department of Yunnan Province Iconographia Arbororum Yunnanicorum. 1988, Kunming: Yunnan Science and Technology Press, 572.

[72] Stegen JC, Swenson NG, Enquist BJ, White EP, Philips OL, Jørgensen PT, Weiser MD, Mendoza AM, Vargas PN. Variation in aboveground forest biomass across broad climatic gradients. Global Ecol Biogeogr, 2011, 20(5): 744-754.

[73] Takyu M, Aiba SI, Kitayama K. Changes in biomass, productivity and decomposition along topographical gradients under different geological conditions in tropical lower montane forests on Mount Kinabalu. Borneo. Oecologia, 2003, 134(3): 397-404.

[74] Temesgen H, Monleon VJ, Hann DW. Analysis and comparison of nonlinear tree height prediction strategies for Douglas–fir forests. Can J For Res, 2008, 38: 553-565.

[75] Ter-Mikaelian MT, Korzukhin MD. Biomass equations for sixty-five North American tree species. For Ecol Manage, 1997, 97: 1-24.

[76] Valencia R, Condit R, Muller-Landau HC, Hernandez C, Navarrete H. Dissecting biomass dynamics in a large Amazonian forest plot. J Trop Ecol, 2009, 25: 473-482.

[77] Vuong QH. Likelihood ratio tests for model selection and non-nested hypotheses. Econometrica, 1989, 57(2): 307-333.

[78] Wang HL (2003) Study on stand growth models for Pinus kesiya var. langbianensis natural secondary woodland. Master Thesis, Kunming: Southwest Forestry College

[79] Wang CK. Biomass allometric equations for 10 co-occurring tree species in Chinese temperate forests. For Ecol Manage, 2006, 222: 9-16.

[80] Weiskittel AR, Maguire DA, Monserud RA. Response of branch growth and mortality to silvicultural treatments in coastal Douglas- fir plantations: implications for predicting tree growth. For Ecol Manage, 2007, 251: 182-194.

[81] Weiskittel AR, Hann DW, Kershaw JA Jr, Vanclay JK. Forest growth and yield modeling. 2011, Chichester: Wiley, 415

[82] Wen QZ, Zhao YF, Chen XM, Yang ZX, Ai JL, Yang XS. Dynamic study on the values for ecological service function of Pinus kesiya forest in China. For Res, 2010, 23(5): 671-677.

[83] West PW. Tree and forest measurement, 2009 2 Berlin: Springer, 192

[84] Wu ZL, Dang CL. The biomass of Pinus kesiya var. langbianensis stands in Pu’er district, Yunnan. J Yunnan Univ (Nat Sci Ed), 1992, 14(2): 161-167.

[85] Wu Y, Strahler AH. Remote estimation of Crown size, stand density, and biomass on the Oregon transect. Ecol Appl, 1994, 4(2): 299-312.

[86] Yan J, Luo YJ, Zheng DF, Wang SC. Source appointment of differences in biomass estimates of Eucalypt Plantation. Sci Silvae Sin, 2014, 50(2): 92-98.

[87] Yue F, Yang B. Study on carbon sink of Pinus kesiya forests. Jiangsu Agric Sci, 2011, 39(5): 467-469.

[88] Zeng WS, Tang SZ, Xia ZS, Zhu S, Luo HZ. Using linear mixed model and dummy variable model approaches to construct generalized single- tree biomass equations in Guizhou. Forest Research, 2011, 24(3): 285-291.

[89] Zhang YJ, Borders BE. Using a system mixed effects modeling method to estimate tree compartment biomass for intensively managed Loblolly pines - an allometric approach. For Ecol Manage, 2004, 194: 145-157.

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

[91] Zianis D, Muukkonen P, Mäkipää R, Mencuccini M (2005) Biomass and stem volume equations for tree species in Europe. Silva Fennica Monographs 4, Tammer-Paino Oy: Tampere, Finland

[92] Zou WT, Zeng WS, Zhang LJ, Zeng M. Modeling crown biomass for four pine species in China. Forests, 2015, 6: 433-449.