Parameterization of the 3-PG model for Quercus mongolica by using tree-ring data and Bayesian calibration

Abstract

Although Quercus mongolica is a widely distributed, economically and ecologically important deciduous tree in northern China, models to accurately predict stand growth at a regional scale are limited. The physiological process model (3-PG) has the potential to predict stand growth dynamics under varying site conditions and climate change scenarios. Here, we used field inventory, tree ring sampling, and Bayesian calibration to parameterize a model for Q. mongolica. Stand volume and productivity were then predicted under present conditions and three future climate scenarios (RCP26, RCP45 and RCP85). Our results demonstrated that after Bayesian calibration, the posterior ranges of the sensitivity parameters aphaCx, wSx1000 and pRn accounted for 34%, 45% and 65%, respectively, of their prior range. Calibration and validation results revealed a strong correlation between predicted and measured values (R² > 0.87, P < 0.01), with < 20% bias for all growth indicators. Stand volume was projected to increase by 145% and productivity by 80% by the year 2100 under the RCP85 scenario, although these projections may vary across regions. The present study developed a tailored set of 3-PG model parameters for Q. mongolica, based on a comprehensive range of climate conditions, stand structure, and age classes. These parameters offer a scientific basis to accurately predict growth of other monospecific oak or mixed-species stands.

Keywords

Quercus mongolica / 3-PG model / Bayesian calibration / Productivity / Growth forecast

References

[1] Almeida AC, Landsberg JJ, Sands PJ. Parameterisation of 3-PG model for fast-growing Eucalyptus grandis plantations. For Ecol Manag, 2004, 193(1–2): 179-195

[2] Boisvenue C, Running SW (2006) Impacts of climate change on natural forest productivity–evidence since the middle of the 20th century. Glob Change Biol 12(5):862–882. https://doi.org/10.1111/j.1365-2486.2006.01134.x

[3] Bose AK, Scherrer D, Julio Camarero J, Ziche D, Babst F, Bigler C, Bolte A, Dorado-Liñán I, Etzold S, Fonti P, Forrester DI, Gavinet J, Gazol A, de Andrés EG, Karger DN, Lebourgeois F, Lévesque M, Martínez-Sancho E, Menzel A, Neuwirth B, Nicolas M, Sanders TGM, Scharnweber T, Schröder J, Zweifel R, Gessler A, Rigling A. Climate sensitivity and drought seasonality determine post-drought growth recovery of Quercus petraea and Quercus robur in Europe. Sci Total Environ, 2021, 784 147222

[4] Bouwman M, Forrester DI, den Ouden J, Nabuurs GJ, Mohren GMJ. Species interactions under climate change in mixed stands of Scots pine and pedunculate oak. For Ecol Manag, 2021, 481 118615

[5] Bryars C, Maier C, Zhao DH, Kane M, Borders B, Will R, Teskey R. Fixed physiological parameters in the 3-PG model produced accurate estimates of loblolly pine growth on sites in different geographic regions. For Ecol Manag, 2013, 289: 501-514

[6] Caldeira DRM, Alvares CA, Campoe OC, Hakamada RE, Guerrini IA, Cegatta ÍR, Stape JL. Multisite evaluation of the 3-PG model for the highest phenotypic plasticity Eucalyptus clone in Brazil. For Ecol Manag, 2020, 462 117989

[7] Carrasco G, Almeida AC, Falvey M, Olmedo GF, Taylor P, Santibañez F, Coops NC. Effects of climate change on forest plantation productivity in Chile. Glob Chang Biol, 2022, 28(24): 7391-7409

[8] Cheng W, Dan L, Deng XZ, Feng JM, Wang YL, Peng J, Tian J, Qi W, Liu Z, Zheng XQ, Zhou DM, Jiang SJ, Zhao HP, Wang XY. Global monthly gridded atmospheric carbon dioxide concentrations under the historical and future scenarios. Sci Data, 2022, 9(1): 83

[9] Coops NC, Hember RA, Waring RH. Assessing the impact of current and projected climates on Douglas-Fir productivity in British Columbia, Canada, using a process-based model (3-PG). Can J for Res, 2010, 40(3511-524

[10] Couture JJ, Meehan TD, Kruger EL, Lindroth RL. Insect herbivory alters impact of atmospheric change on northern temperate forests. Nat Plants, 2015, 1: 15016

[11] Dye PJ, Jacobs S, Drew D. Verification of 3-PG growth and water-use predictions in twelve Eucalyptus plantation stands in Zululand. South Africa for Ecol Manag, 2004, 193(1–2): 197-218

[12] Esprey LJ, Sands PJ, Smith CW. Understanding 3-PG using a sensitivity analysis. For Ecol Manag, 2004, 193(1–2235-250

[13] Feikema P, Morris J, Beverly CR, Lane PNJ, Baker TG. Using 3PG to simulate longterm growth and transpiration in Eucalyptus regnans forests. Int Congr Environ Model Softw, 2010, 326: 20

[14] Feng F, Wang KC. Merging high-resolution satellite surface radiation data with meteorological sunshine duration observations over China from 1983 to 2017. Remote Sens, 2021, 13(4): 602

[15] Fontes L, Bontemps JD, Bugmann H, Van Oijen M, Gracia C, Kramer K, Lindner M, Rötzer T, Skovsgaard JP. Models for supporting forest management in a changing environment. For Syst, 2010, 19: 8-29

[16] Forrester DI, Tang XL. Analysing the spatial and temporal dynamics of species interactions in mixed-species forests and the effects of stand density using the 3-PG model. Ecol Model, 2016, 319: 233-254

[17] Forrester DI, Hobi ML, Mathys AS, Stadelmann G, Trotsiuk V. Calibration of the process-based model 3-PG for major central European tree species. Eur J for Res, 2021, 140(4847-868

[18] Forrester DI, Mathys AS, Stadelmann G, Trotsiuk V. Effects of climate on the growth of Swiss uneven-aged forests: combining >100 years of observations with the 3-PG model. For Ecol Manag, 2021, 494 119271

[19] Fu LY, Zhang HR, Lu J, Zang H, Lou MH, Wang GX (2015) Multilevel nonlinear mixed-effect crown ratio models for individual trees of Mongolian oak (Quercus mongolica) in NorthEast China. PLoS One 10(8):e0133294. https://doi.org/10.1371/journal.pone.0133294

[20] Gang CC, Zhang YZ, Wang ZQ, Chen YZ, Yang Y, Li JL, Cheng JM, Qi JG, Odeh I. Modeling the dynamics of distribution, extent, and NPP of global terrestrial ecosystems in response to future climate change. Glob Planet Change, 2017, 148: 153-165

[21] Gao WQ, Lei XD, Fu LY, Duan GS, Zhou ML, Cao J. Radial growth response of two oaks to climate at their disparate distribution limits in semiarid areas, Beijing China. Ecosphere, 2020, 11(2 e03062

[22] Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Statist Sci, 1992, 7(4): 457-472

[23] Gertner GZ, Fang SF, Skovsgaard JP. A Bayesian approach for estimating the parameters of a forest process model based on long-term growth data. Ecol Model, 1999, 119(2–3): 249-265

[24] Gong XF, Yuan DY, Zhu LJ, Li ZS, Wang XC. Long-term changes in radial growth of seven tree species in the mixed broadleaf-Korean pine forest in NorthEast China: are deciduous trees favored by climate change?. J Forestry Res, 2024, 35(1): 70

[25] Gonzalez-Benecke CA, Teskey RO, Martin TA, Jokela EJ, Fox TR, Kane MB, Noormets A. Regional validation and improved parameterization of the 3-PG model for Pinus taeda stands. For Ecol Manag, 2016, 361: 237-256

[26] González-García M, Almeida AC, Hevia A, Majada J, Beadle C. Application of a process-based model for predicting the productivity of Eucalyptus nitens bioenergy plantations in Spain. GCB Bioenergy, 2016, 8(1194-210

[27] Hung TT, Almeida AC, Eyles A, Mohammed C. Predicting productivity of Acacia hybrid plantations for a range of climates and soils in Vietnam. For Ecol Manag, 2016, 367: 97-111

[28] Johnsen K, Samuelson L, Teskey R, McNulty S, Fox T. Process models as tools in forestry research and management. For Sci, 2001, 47(1): 2-8

[29] Klesse S, Babst F, Lienert S, Spahni R, Joos F, Bouriaud O, Carrer M, Di Filippo A, Poulter B, Trotsiuk V, Wilson R, Frank DC. A combined tree ring and vegetation model assessment of European forest growth sensitivity to interannual climate variability. Glob Biogeochem Cycles, 2018, 32(8): 1226-1240

[30] Landsberg J, Sands P. Physiological ecology of forest production: principles, processes and models, 2011, USA, Academic Press

[31] Landsberg JJ, Waring RH. A generalised model of forest productivity using simplified concepts of radiation-use efficiency, carbon balance and partitioning. For Ecol Manag, 1997, 95(3): 209-228

[32] Li J. Flora of China. Harv Pap Bot, 2007, 13: 301-302

[33] Li WY, Gu WC. Study on phenotypic diversity of natural population in Quercus mongolica. Sci Silvae Sin, 2005, 41(1): 49-56in Chinese

[34] Liu CJ, Zheng XM, Ren Y. Parameter optimization of the 3PG model based on sensitivity analysis and a Bayesian method. Forests, 2020, 11(12): 1369

[35] Luo YQ, Weng ES, Wu XW, Gao C, Zhou XH, Zhang L. Parameter identifiability, constraint, and equifinality in data assimilation with ecosystem models. Ecol Appl, 2009, 19(3571-574

[36] Morris MD. Factorial sampling plans for preliminary computational experiments. Technometrics, 1991, 33(2): 161-174

[37] National Forestry and Grassland Administration. China forest resources report (2014–2018), 2019, Beijing, China, China Forestry Publishing Housein Chinese

[38] Nölte A, Yousefpour R, Hanewinkel M. Changes in sessile oak (Quercus petraea) productivity under climate change by improved leaf phenology in the 3-PG model. Ecol Model, 2020, 438 109285

[39] Peng CH. Growth and yield models for uneven-aged stands: past, present and future. For Ecol Manag, 2000, 132(2–3): 259-279

[40] Pérez-Cruzado C, Muñoz-Sáez F, Basurco F, Riesco G, Rodríguez-Soalleiro R. Combining empirical models and the process-based model 3-PG to predict Eucalyptus nitens plantations growth in Spain. For Ecol Manag, 2011, 262(6): 1067-1077

[41] Perkins D, Uhl E, Biber P, Du Toit B, Carraro V, Rötzer T, Pretzsch H (2018) Impact of climate trends and drought events on the growth of oaks (Quercus robur L. and Quercus petraea (Matt.) liebl.) within and beyond their natural range. Forests 9(3):108. https://doi.org/10.3390/f9030108

[42] R Core Team (2022) R: A language and environment for statistical computing.

[43] Rollinson CR, Liu Y, Raiho A, Moore DJP, McLachlan J, Bishop DA, Dye A, Matthes JH, Hessl A, Hickler T, Pederson N, Poulter B, Quaife T, Schaefer K, Steinkamp J, Dietze MC. Emergent climate and CO₂ sensitivities of net primary productivity in ecosystem models do not agree with empirical data in temperate forests of eastern North America. Glob Change Biol, 2017, 23(7): 2755-2767

[44] Sands PJ, Landsberg JJ. Parameterisation of 3-PG for plantation grown Eucalyptus globulus. For Ecol Manag, 2002, 163(1–3): 273-292

[45] Schiermeier Q. Droughts, heatwaves and floods: how to tell when climate change is to blame. Nature, 2018, 560(7716): 20-22

[46] Subedi S, Fox T, Wynne R. Determination of fertility rating (FR) in the 3-PG model for loblolly pine plantations in the southeastern United States based on site index. Forests, 2015, 6(9): 3002-3027

[47] ter Braak CJF, Vrugt JA. Differential evolution Markov Chain with snooker updater andfewer chains. Stat Comput, 2008, 18(4): 435-446

[48] Terrer C, Jackson RB, Prentice IC, Keenan TF, Kaiser C, Vicca S, Fisher JB, Reich PB, Stocker BD, Hungate BA, Peñuelas J, McCallum I, Soudzilovskaia NA, Cernusak LA, Talhelm AF, Van Sundert K, Piao SL, Newton PCD, Hovenden MJ, Blumenthal DM, Liu YY, Müller C, Winter K, Field CB, Viechtbauer W, Van Lissa CJ, Hoosbeek MR, Watanabe M, Koike T, Leshyk VO, Polley HW, Franklin O. Nitrogen and phosphorus constrain the CO₂ fertilization of global plant biomass. Nat Clim Change, 2019, 9(9): 684-689

[49] Touhami I, Chirino E, Aouinti H, El Khorchani A, Elaieb MT, Khaldi A, Nasr Z. Decline and dieback of cork oak (Quercus suber L.) forests in the Mediterranean basin: a case study of Kroumirie Northwest Tunisia. J For Res, 2020, 31(5): 1461-1477

[50] Trotsiuk V, Hartig F, Forrester DI. r3PG–An r package for simulating forest growth using the 3-PG process-based model. Meth Ecol Evol, 2020, 11(11): 1470-1475

[51] Van Oijen M, Rougier J, Smith R. Bayesian calibration of process-based forest models: bridging the gap between models and data. Tree Physiol, 2005, 25(7): 915-927

[52] van Oijen M. Bayesian methods for quantifying and reducing uncertainty and error in forest models. Curr for Rep, 2017, 3(4): 269-280

[53] Wang LM, Ren XW, Liu YQ. Geographic distribution of deciduous oaks in China. J Beijing for Univ, 1985, 7(2): 57-69in Chinese

[54] Wang TL, Wang GY, Innes JL, Seely B, Chen BZ. ClimateAP: an application for dynamic local downscaling of historical and future climate data in Asia Pacific. Front Agr Sci Eng, 2017, 4(4): 448

[55] Wang B, Waters C, Anwar MR, Cowie A, Liu DL, Summers D, Paul K, Feng PY. Future climate impacts on forest growth and implications for carbon sequestration through reforestation in southeast Australia. J Environ Manage, 2022, 302(Pt A 113964

[56] Wang JC, Li SH, Li Q, Guo YL, Han YJ, Liu Q, Zhao YQ. Interpreting Larix principis-rupprechtii radial growth in northern China based on the Vaganov-Shashkin model. J for Res, 2024, 35(147

[57] Wei L, Marshall JD, Link TE, Kavanagh KL, Du EH, Pangle RE, Gag PJ, Ubierna N. Constraining 3-PG with a new δ¹³C submodel: a test using the δ¹³C of tree rings. Plant Cell Environ, 2014, 37(182-100

[58] Wu CY, Chen DS, Shen JP, Sun XM, Zhang SG. Estimating the distribution and productivity characters of Larix kaempferi in response to climate change. J Environ Manag, 2021, 280 111633

[59] Xenakis G, Ray D, Mencuccini M. Sensitivity and uncertainty analysis from a coupled 3-PG and soil organic matter decomposition model. Ecol Model, 2008, 219(1–2): 1-16

[60] Xie YL, Wang HY, Lei XD. Application of the 3-PG model to predict growth of Larix olgensis plantations in northeastern China. For Ecol Manag, 2017, 406: 208-218

[61] Xie YL, Wang HY, Lei XD. Simulation of climate change and thinning effects on productivity of Larix olgensis plantations in NorthEast China using 3-PG_mix model. J Environ Manage, 2020, 261 110249

[62] Xu K, Wang XP, Liang PH, An HL, Sun H, Han W, Li QY. Tree-ring widths are good proxies of annual variation in forest productivity in temperate forests. Sci Rep, 2017, 7(1): 1945

[63] Xu K, Wang XP, Liang PH, Wu YL, An HL, Sun H, Wu P, Wu X, Li QY, Guo X, Wen XS, Han W, Liu C, Fan DY. A new tree-ring sampling method to estimate forest productivity and its temporal variation accurately in natural forests. For Ecol Manag, 2019, 433: 217-227

[64] Yu SL, Ma KP, Chen LZ. Preliminary discussion on the origin of Quercus mongolica forest in North China. Guihaia, 2000, 20(2131-137in Chinese

[65] Yu XZ, Zhong LH, Zhou H, Gong L, Zhao Y, Wei L. Tree-ring based forest model calibrations with a deep learning algorithm. For Ecol Manag, 2024, 569 122154

[66] Zhang J, Ge Y, Chang J, Jiang B, Jiang H, Peng CH, Zhu JR, Yuan WG, Qi LZ, Yu SQ. Carbon storage by ecological service forests in Zhejiang Province, subtropical China. For Ecol Manag, 2007, 245(1–3): 64-75

[67] Zobitz JM, Desai AR, Moore DJP, Chadwick MA. A primer for data assimilation with ecological models using Markov Chain Monte Carlo (MCMC). Oecologia, 2011, 167(3): 599-611

RIGHTS & PERMISSIONS

Northeast Forestry University