The multi-assumption architecture and testbed (MAAT v1.0): R code for generating ensembles with dynamic model structure and analysis of epistemic uncertainty from multiple sources

Abstract. Computer models are ubiquitous tools used to represent systems across many scientific and engineering domains. For any given system, many computer models exist, each built on different assumptions and demonstrating variability in the ways in which these systems can be represented. This variability is known as epistemic uncertainty, i.e. uncertainty in our knowledge of how these systems operate. Two primary sources of epistemic uncertainty are (1) uncertain parameter values and (2) uncertain mathematical representations of the processes that comprise the system. Many formal methods exist to analyse parameter-based epistemic uncertainty, while process-representation-based epistemic uncertainty is often analysed post hoc, incompletely, informally, or is ignored. In this model description paper we present the multi-assumption architecture and testbed (MAAT v1.0) designed to formally and completely analyse process-representation-based epistemic uncertainty. MAAT is a modular modelling code that can simply and efficiently vary model structure (process representation), allowing for the generation and running of large model ensembles that vary in process representation, parameters, parameter values, and environmental conditions during a single execution of the code. MAAT v1.0 approaches epistemic uncertainty through sensitivity analysis, assigning variability in model output to processes (process representation and parameters) or to individual parameters. In this model description paper we describe MAAT and, by using a simple groundwater model example, verify that the sensitivity analysis algorithms have been correctly implemented. The main system model currently coded in MAAT is a unified, leaf-scale enzyme kinetic model of C3 photosynthesis. In the Appendix we describe the photosynthesis model and the unification of multiple representations of photosynthetic processes. The numerical solution to leaf-scale photosynthesis is verified and examples of process variability in temperature response functions are provided. For rapid application to new systems, the MAAT algorithms for efficient variation of model structure and sensitivity analysis are agnostic of the specific system model employed. Therefore MAAT provides a tool for the development of novel or toy models in many domains, i.e. not only photosynthesis, facilitating rapid informal and formal comparison of alternative modelling approaches.

[1]  Graham D. Farquhar,et al.  An Empirical Model of Stomatal Conductance , 1984 .

[2]  I. E. Woodrow,et al.  A Model Predicting Stomatal Conductance and its Contribution to the Control of Photosynthesis under Different Environmental Conditions , 1987 .

[3]  G. Farquhar,et al.  Effect of temperature on the CO2/O2 specificity of ribulose-1,5-bisphosphate carboxylase/oxygenase and the rate of respiration in the light , 1985, Planta.

[4]  Jens Kattge,et al.  Temperature acclimation in a biochemical model of photosynthesis: a reanalysis of data from 36 species. , 2007, Plant, cell & environment.

[5]  P. Jones,et al.  The Twentieth Century Reanalysis Project , 2009 .

[6]  Zong-Liang Yang,et al.  Technical description of version 4.5 of the Community Land Model (CLM) , 2013 .

[7]  Cajo J. F. ter Braak,et al.  Equifinality of formal (DREAM) and informal (GLUE) Bayesian approaches in hydrologic modeling? , 2009 .

[8]  Stephen D. Sebestyen,et al.  Biophysical drivers of seasonal variability in Sphagnum gross primary production in a northern temperate bog , 2017 .

[9]  S. V. Caemmerer,et al.  Biochemical models of leaf photosynthesis. , 2000 .

[10]  Ming Ye,et al.  Variance-based global sensitivity analysis for multiple scenarios and models with implementation using sparse grid collocation , 2015 .

[11]  Stephen W Pacala,et al.  Optimal stomatal behavior with competition for water and risk of hydraulic impairment , 2016, Proceedings of the National Academy of Sciences.

[12]  F. Molteni,et al.  The ECMWF Ensemble Prediction System: Methodology and validation , 1996 .

[13]  D. Baldocchi An analytical solution for coupled leaf photosynthesis and stomatal conductance models. , 1994, Tree physiology.

[14]  I. Colinprentice,et al.  A simulation model for the transient effects of climate change on forest landscapes , 1993 .

[15]  Atul K. Jain,et al.  Forest water use and water use efficiency at elevated CO2: a model‐data intercomparison at two contrasting temperate forest FACE sites , 2013, Global change biology.

[16]  Reto Knutti,et al.  The use of the multi-model ensemble in probabilistic climate projections , 2007, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[17]  R. Knutti,et al.  Climate model genealogy , 2011 .

[18]  J. Whitaker,et al.  The Twentieth Century Reanalysis Project 3 , 2011 .

[19]  Dan Bruhn,et al.  Tracking the origins of the Kok effect, 70 years after its discovery. , 2017, The New phytologist.

[20]  Keith Beven,et al.  Equifinality, data assimilation, and uncertainty estimation in mechanistic modelling of complex environmental systems using the GLUE methodology , 2001 .

[21]  Suan Chin Wong,et al.  A simple calibrated model of Amazon rainforest productivity based on leaf biochemical properties , 1995 .

[22]  Peter A. Vanrolleghem,et al.  Uncertainty in the environmental modelling process - A framework and guidance , 2007, Environ. Model. Softw..

[23]  Ray Leuning,et al.  Modelling Stomatal Behaviour and and Photosynthesis of Eucalyptus grandis , 1990 .

[24]  Atul K. Jain,et al.  Using ecosystem experiments to improve vegetation models , 2015 .

[25]  I. C. Prentice,et al.  Global climatic drivers of leaf size , 2017, Science.

[26]  P. Ciais,et al.  Spatiotemporal patterns of terrestrial gross primary production: A review , 2015 .

[27]  Pierre Friedlingstein,et al.  Uncertainties in CMIP5 Climate Projections due to Carbon Cycle Feedbacks , 2014 .

[28]  Stan D. Wullschleger,et al.  Biochemical Limitations to Carbon Assimilation in C3 Plants—A Retrospective Analysis of the A/Ci Curves from 109 Species , 1993 .

[29]  P. Struik,et al.  C3 and C4 photosynthesis models: An overview from the perspective of crop modelling , 2009 .

[30]  Stephen Sitch,et al.  Effects of parameter uncertainties on the modeling of terrestrial biosphere dynamics , 2005 .

[31]  Ian J. Wright,et al.  Impacts of trait variation through observed trait–climate relationships on performance of an Earth system model: a conceptual analysis , 2012 .

[32]  Ying Ping Wang,et al.  A comparison of three different canopy radiation models commonly used in plant modelling. , 2003, Functional plant biology : FPB.

[33]  Mark G. Tjoelker,et al.  Modelling respiration of vegetation: evidence for a general temperature‐dependent Q10 , 2001 .

[34]  Susanne von Caemmerer,et al.  Coupled response of stomatal and mesophyll conductance to light enhances photosynthesis of shade leaves under sunflecks. , 2016, Plant, cell & environment.

[35]  James F. Reynolds,et al.  Coordination theory of leaf nitrogen distribution in a canopy , 1993, Oecologia.

[36]  Michael Bahn,et al.  Inter‐specific variation of the biochemical limitation to photosynthesis and related leaf traits of 30 species from mountain grassland ecosystems under different land use , 1999 .

[37]  Emil L. Smith THE INFLUENCE OF LIGHT AND CARBON DIOXIDE ON PHOTOSYNTHESIS , 1937, The Journal of general physiology.

[38]  P. Green Reversible jump Markov chain Monte Carlo computation and Bayesian model determination , 1995 .

[39]  Paola Annoni,et al.  Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index , 2010, Comput. Phys. Commun..

[40]  Florian A Busch,et al.  Changes in the chloroplastic CO2 concentration explain much of the observed Kok effect: a model. , 2017, The New phytologist.

[41]  Richard Harding,et al.  A canopy conductance and photosynthesis model for use in a GCM land surface scheme , 1998 .

[42]  R. McMurtrie,et al.  Long-Term Response of Nutrient-Limited Forests to CO"2 Enrichment; Equilibrium Behavior of Plant-Soil Models. , 1993, Ecological applications : a publication of the Ecological Society of America.

[43]  Allison R Kolbe,et al.  Mesophyll conductance in Zea mays responds transiently to CO2 availability: implications for transpiration efficiency in C4 crops. , 2018, The New phytologist.

[44]  Reto Knutti,et al.  Challenges in Combining Projections from Multiple Climate Models , 2010 .

[45]  Pierre Martre,et al.  The Coordination of Leaf Photosynthesis Links C and N Fluxes in C3 Plant Species , 2012, PloS one.

[46]  Lianhong Gu,et al.  Reliable estimation of biochemical parameters from C₃ leaf photosynthesis-intercellular carbon dioxide response curves. , 2010, Plant, cell & environment.

[47]  J. R. Evans Photosynthesis and nitrogen relationships in leaves of C3 plants , 2004, Oecologia.

[48]  C.J.T. Spitters,et al.  Separating the diffuse and direct component of global radiation and its implications for modeling canopy photosynthesis Part II. Calculation of canopy photosynthesis , 1986 .

[49]  James F. Reynolds,et al.  Modelling photosynthesis of cotton grown in elevated CO2 , 1992 .

[50]  Keith Beven,et al.  A manifesto for the equifinality thesis , 2006 .

[51]  G. Collatz,et al.  Physiological and environmental regulation of stomatal conductance, photosynthesis and transpiration: a model that includes a laminar boundary layer , 1991 .

[52]  Stephen Sitch,et al.  A roadmap for improving the representation of photosynthesis in Earth system models. , 2017, The New phytologist.

[53]  Scott L. Painter,et al.  Managing complexity in simulations of land surface and near-surface processes , 2016, Environ. Model. Softw..

[54]  Keith Beven,et al.  Facets of uncertainty: epistemic uncertainty, non-stationarity, likelihood, hypothesis testing, and communication , 2016 .

[55]  Karsten Lehmann,et al.  Selecting a climate model subset to optimise key ensemble properties , 2017 .

[56]  Alain Dassargues,et al.  Conceptual model uncertainty in groundwater modeling: Combining generalized likelihood uncertainty estimation and Bayesian model averaging , 2008 .

[57]  Ashehad A. Ali,et al.  A global scale mechanistic model of photosynthetic capacity (LUNA V1.0) , 2016 .

[58]  Markus Reichstein,et al.  Improving canopy processes in the Community Land Model version 4 (CLM4) using global flux fields empirically inferred from FLUXNET data , 2011 .

[59]  Paolo De Angelis,et al.  Reconciling the optimal and empirical approaches to modelling stomatal conductance , 2011 .

[60]  Fred L. Ogden,et al.  GSSHA: Model To Simulate Diverse Stream Flow Producing Processes , 2004 .

[61]  Susan E. Trumbore,et al.  Models of soil organic matter decomposition: the SoilR package, version 1.0 , 2012 .

[62]  Saman Razavi,et al.  What do we mean by sensitivity analysis? The need for comprehensive characterization of “global” sensitivity in Earth and Environmental systems models , 2015 .

[63]  F. Woodward,et al.  The relationship of leaf photosynthetic traits – Vcmax and Jmax – to leaf nitrogen, leaf phosphorus, and specific leaf area: a meta-analysis and modeling study , 2014, Ecology and evolution.

[64]  J. Randerson,et al.  Plant responses to increasing CO2 reduce estimates of climate impacts on drought severity , 2016, Proceedings of the National Academy of Sciences.

[65]  Ming Ye,et al.  Global sensitivity analysis in hydrological modeling: Review of concepts, methods, theoretical framework, and applications , 2015 .

[66]  Edward B. Rastetter,et al.  Modeling for Understanding v. Modeling for Numbers , 2017, Ecosystems.

[67]  M. Jansen Analysis of variance designs for model output , 1999 .

[68]  Dmitri Kavetski,et al.  A unified approach for process‐based hydrologic modeling: 1. Modeling concept , 2015 .

[69]  Eric Klopfer,et al.  Modeling for Understanding , 2000 .

[70]  Ming Ye,et al.  A new process sensitivity index to identify important system processes under process model and parametric uncertainty , 2017 .

[71]  Lawren Sack,et al.  Optimal plant water economy. , 2017, Plant, cell & environment.

[72]  C. Jacobs,et al.  Direct impact of atmospheric CO2 enrichment on regional transpiration , 1994 .

[73]  Jeroen P. van der Sluijs,et al.  A framework for dealing with uncertainty due to model structure error , 2004 .

[74]  Denis Loustau,et al.  Temperature response of parameters of a biochemically based model of photosynthesis. II. A review of experimental data , 2002 .

[75]  F. Woodward,et al.  The impact of alternative trait-scaling hypotheses for the maximum photosynthetic carboxylation rate (Vcmax ) on global gross primary production. , 2017, The New phytologist.

[76]  Ming Ye,et al.  Towards a comprehensive assessment of model structural adequacy , 2012 .

[77]  R. Knutti,et al.  Robustness and uncertainties in the new CMIP5 climate model projections , 2013 .

[78]  Carl J. Bernacchi,et al.  Improved temperature response functions for models of Rubisco‐limited photosynthesis , 2001 .

[79]  A. Jameson,et al.  Optimum Aerodynamic Design Using the Navier–Stokes Equations , 1997 .

[80]  X. Li,et al.  Application of Bayesian Model Averaging in the Reconstruction of Past Climate Change Using PMIP3/CMIP5 Multimodel Ensemble Simulations , 2016 .

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

[82]  Susanne von Caemmerer,et al.  Temperature responses of mesophyll conductance differ greatly between species. , 2015, Plant, cell & environment.

[83]  Ichiro Terashima,et al.  Resistances along the CO2 diffusion pathway inside leaves. , 2009, Journal of experimental botany.

[84]  J. Randerson,et al.  Technical Description of version 4.0 of the Community Land Model (CLM) , 2010 .