Variance-based sensitivity analysis of BIOME-BGC for gross and net primary production

Parameterization and calibration of a process-based simulator (PBS) is a major challenge when simulating gross and net primary production (GPP and NPP). The large number of parameters makes the calibration computationally expensive and is complicated by the dependence of several parameters on other parameters. Calibration can be simplified by first identifying those parameters for which GPP and NPP are most sensitive. For an appropriate application of a PBS, a sensitivity analysis is an essential step. Sensitivity analysis based on local derivatives (i.e., one-at-a-time analysis) does not examine the PBS behaviour over the whole parameter space. This study therefore implements a variance-based sensitivity analysis (VBSA) addressing the full range of PBS input. A VBSA is also independent of non-linearity in a PBS. This paper performs a VBSA of the process-based simulator BIOME-BGC for GPP and NPP output in a Douglas-fir stand at the Speulderbos forest site, The Netherlands. The results show that GPP and NPP are highly sensitive to the following parameters: fraction of leaf nitrogen in Rubisco, the ratio of fine root carbon to leaf carbon, the ratio of carbon to nitrogen in leaf and fine root, the leaf and fine root turnover, the water interception coefficient and soil depth. GPP and NPP are particularly sensitive to the ratio of fine root carbon to leaf carbon that is responsible for leaf area index development. The study concludes that a VBSA analysis provides a reliable and useful approach for a sensitivity analysis of process-based simulators with a complicated structure in the parameters.

[1]  S. Running,et al.  Contrasting Climatic Controls on the Estimated Productivity of Global Terrestrial Biomes , 1998, Ecosystems.

[2]  N. A. S. Hamm,et al.  Variance-based sensitivity analysis of the probability of hydrologically induced slope instability , 2006, Comput. Geosci..

[3]  Pebesma Edzer,et al.  Managing uncertainty in integrated environmental modelling frameworks: The UncertWeb framework , 2012 .

[4]  J. Oakley.  Eliciting Gaussian process priors for complex computer codes , 2002 .

[5]  S. Dekker.  Modelling and monitoring forest evapotranspiration. Behaviour, concepts and parameters , 2000 .

[6]  Ben Bond-Lamberty,et al.  Reimplementation of the Biome-BGC model to simulate successional change. , 2005, Tree physiology.

[7]  Godefridus M. J. Mohren,et al.  Inter-specific competition in mixed forests of Douglas-fir (Pseudotsuga menziesii) and common beech (Fagus sylvatica) under climate change — a model-based analysis , 2010, Annals of Forest Science.

[8]  Yaseen T. Mustafa.  Improving forest growth estimation: Bayesian networks for integrating satellite images and process-based forest growth models , 2012 .

[9]  Peter E. Thornton,et al.  Parameterization and Sensitivity Analysis of the BIOME–BGC Terrestrial Ecosystem Model: Net Primary Production Controls , 2000 .

[10]  Maosheng Zhao,et al.  A Continuous Satellite-Derived Measure of Global Terrestrial Primary Production , 2004 .

[11]  J. Berry,et al.  A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species , 1980, Planta.

[12]  Peter E. Thornton,et al.  Modeling and measuring the effects of disturbance history and climate on carbon and water budgets in evergreen needleleaf forests , 2002 .

[13]  K. Pilegaard,et al.  Oxidized nitrogen and ozone interaction with forests. II: Multi‐layer process‐oriented modelling results and a sensitivity study for Douglas fir , 2004 .

[14]  R. Keane,et al.  Ecophysiological parameters for Pacific Northwest trees. , 2004 .

[15]  Peter A. Coppin,et al.  Parameter estimation in surface exchange models using nonlinear inversion: how many parameters can we estimate and which measurements are most useful? , 2001 .

[16]  Mirco Migliavacca,et al.  Constraining model simulations of GPP using satellite retrieved leaf chlorophyll , 2012, 2012 IEEE International Geoscience and Remote Sensing Symposium.

[17]  G. Churkina,et al.  Parameter estimation and validation of the terrestrial ecosystem model Biome-BGC using eddy-covariance flux measurements , 2009 .

[18]  The nitrogen cycle and its influence on the European greenhouse gas balance , 2011 .

[19]  A. Saltelli,et al.  Importance measures in global sensitivity analysis of nonlinear models , 1996 .

[20]  Benedict M. Mutua,et al.  Sensitivity analysis and calibration of the Modified Universal Soil Loss Equation (MUSLE) for the upper Malewa Catchment, Kenya. , 2013 .

[21]  Markus Reichstein,et al.  CO2 balance of boreal, temperate, and tropical forests derived from a global database , 2007 .

[22]  W. Bouten,et al.  Soil water dynamics and long-term water balances of a Douglas Fir stand in the Netherlands , 1994 .

[23]  I. Hajnsek,et al.  EAGLE 2006 - multi-purpose, multi-angle and multi-sensor in-situ and airborne campaigns over grassland and forest. , 2009 .

[24]  Inter- and intra-annual variations in canopy fine litterfall and carbon and nitrogen inputs to the forest floor in two European coniferous forests , 2013, Annals of Forest Science.

[25]  D. Simpson,et al.  Deposition and emissions of reactive nitrogen over European forests: A modelling study , 2006 .

[26]  N. Gobron,et al.  Diagnostic assessment of European gross primary production , 2008 .

[27]  W. Kustas,et al.  Satellite retrievals of leaf chlorophyll and photosynthetic capacity for improved modeling of GPP , 2013 .

[28]  Willem Bouten,et al.  Simulation of carbon and water budgets of a Douglas fir forest , 2001 .

[29]  Park S. Nobel,et al.  Physicochemical and Environmental Plant Physiology , 1991 .

[30]  S. Running,et al.  8 – Generalization of a Forest Ecosystem Process Model for Other Biomes, BIOME-BGC, and an Application for Global-Scale Models , 1993 .

[31]  F. Bosveld,et al.  Water storage and evaporation as constituents of rainfall interception , 1998 .

[32]  Piero Toscano,et al.  Remote sensing of annual terrestrial gross primary productivity from MODIS: an assessment using the FLUXNET La Thuile data set , 2013 .

[33]  P. Reich,et al.  Generality of leaf trait relationships: a test across six biomes: Ecology , 1999 .

[34]  Gérard Dedieu,et al.  Methodology for the estimation of terrestrial net primary production from remotely sensed data , 1994 .

[35]  M. Harmon,et al.  Leaf-Litter Decomposition in the Picea/Tsuga Forests of Olympic National Park, Washington, U.S.A. , 1990 .

[36]  J. Domec,et al.  Effects of age-related increases in sapwood area, leaf area, and xylem conductivity on height-related hydraulic costs in two contrasting coniferous species , 2011, Annals of Forest Science.

[37]  Emil Cienciala,et al.  Application of BIOME-BGC model to managed forests: 1. Sensitivity analysis , 2006 .

[38]  R. Edmonds.  Decomposition rates and nutrient dynamics in small-diameter woody litter in four forest ecosystems in Washington, U.S.A. , 1987 .

[39]  P. Thornton,et al.  Ecosystem model spin-up: Estimating steady state conditions in a coupled terrestrial carbon and nitrogen cycle model , 2005 .

[40]  Dan Cornford,et al.  Managing uncertainty in integrated environmental modelling: The UncertWeb framework , 2013, Environ. Model. Softw..

[41]  Saltelli Andrea,et al.  Global Sensitivity Analysis: The Primer , 2008 .

[42]  G. Mohren,et al.  Modeling daily gas exchange of a Douglas-fir forest: comparison of three stomatal conductance models with and without a soil water stress function. , 2000, Tree physiology.

[43]  C. C. Grier,et al.  Above- and below-ground net production in 40-year-old Douglas-fir stands on low and high productivity sites , 1981 .

[44]  P. Reich,et al.  Fine root decomposition rates do not mirror those of leaf litter among temperate tree species , 2009, Oecologia.