Equifinality, sloppiness, and emergent structures of mechanistic soil biogeochemical models

Abstract Biogeochemical models increasingly consider the microbial control of carbon cycling in soil. The major current challenge is to validate mechanistic descriptions of microbial processes and predicted system responses against experimental observations. We analyzed soil biochemical models of different complexity regarding parameter identifiability using information geometry, i.e. a model is geometrically interpreted as a manifold embedded in data space. The most complex model (PECCAD) was used as a test case to reveal parsimonious process formulations. All models showed sloppiness, i.e. most individual parameter values cannot be inferred from the observed data. We derived a less complex model formulation of PECCAD with effective inferable parameter combinations and identified structural model limitations. The complexity of identified effective models was systematically reduced with decreasing information content of data. Our results suggest that information geometry provides a powerful approach to derive effective descriptions of relevant biogeochemical processes and reduce structural model uncertainty.

[1]  A. Porporato,et al.  A theoretical analysis of microbial eco-physiological and diffusion limitations to carbon cycling in drying soils , 2014 .

[2]  W. Nowak,et al.  Model selection on solid ground: Rigorous comparison of nine ways to evaluate Bayesian model evidence , 2014, Water resources research.

[3]  Wayne M Getz,et al.  Making ecological models adequate. , 2018, Ecology letters.

[4]  G. Ågren,et al.  Quality or decomposer efficiency - which is most important in the temperature response of litter decomposition? A modelling study using the GLUE methodology , 2010 .

[5]  Peng Qiu,et al.  Experimental design and model reduction in systems biology , 2018, Quantitative Biology.

[6]  Jens Timmer,et al.  Cause and cure of sloppiness in ordinary differential equation models. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[7]  Mark K. Transtrum,et al.  Information Geometry Approach to Verification of Dynamic Models in Power Systems , 2018, IEEE Transactions on Power Systems.

[8]  Kris Villez,et al.  Global parameter optimization for biokinetic modeling of simple batch experiments , 2016, Environ. Model. Softw..

[9]  Amilcare Porporato,et al.  Soil carbon and nitrogen mineralization: Theory and models across scales , 2009 .

[10]  James V. Stone Bayes' Rule: A Tutorial Introduction to Bayesian Analysis , 2013 .

[11]  Jim W. Hall,et al.  Sensitivity analysis of environmental models: A systematic review with practical workflow , 2014, Environ. Model. Softw..

[12]  Christopher R. Myers,et al.  Python Unleashed on Systems Biology , 2007, Computing in Science & Engineering.

[13]  C. Sierra,et al.  A general mathematical framework for representing soil organic matter dynamics , 2015 .

[14]  T. Egli,et al.  Kinetic models for the growth of Escherichia coli with mixtures of sugars under carbon-limited conditions. , 1998, Biotechnology and bioengineering.

[15]  Jessica A. M. Moore,et al.  Multiple models and experiments underscore large uncertainty in soil carbon dynamics , 2018, Biogeochemistry.

[16]  Philipp Kügler,et al.  Moment Fitting for Parameter Inference in Repeatedly and Partially Observed Stochastic Biological Models , 2012, PloS one.

[17]  Michael Hucka,et al.  MOCCASIN: converting MATLAB ODE models to SBML , 2016, Bioinform..

[18]  Yujie He,et al.  Explicitly representing soil microbial processes in Earth system models , 2015 .

[19]  P. Kügler,et al.  Parameter identification for chemical reaction systems using sparsity enforcing regularization: a case study for the chlorite-iodide reaction. , 2009, The journal of physical chemistry. A.

[20]  Saman Razavi,et al.  Global sensitivity analysis for high-dimensional problems: How to objectively group factors and measure robustness and convergence while reducing computational cost , 2019, Environ. Model. Softw..

[21]  Steven D. Allison,et al.  Building Predictive Models for Diverse Microbial Communities in Soil , 2017 .

[22]  M. B. Beck,et al.  Assessing local structural identifiability for environmental models , 2017, Environ. Model. Softw..

[23]  Lianhong Gu,et al.  Microbial dormancy improves development and experimental validation of ecosystem model , 2014, The ISME Journal.

[24]  Max D. Morris,et al.  Factorial sampling plans for preliminary computational experiments , 1991 .

[25]  Mark K. Transtrum,et al.  Manifold boundaries give "gray-box" approximations of complex models , 2016 .

[26]  Philipp Kraft,et al.  Rejecting hydro-biogeochemical model structures by multi-criteria evaluation , 2017, Environ. Model. Softw..

[27]  Nuno Carvalhais,et al.  Implications of the carbon cycle steady state assumption for biogeochemical modeling performance and inverse parameter retrieval , 2008 .

[28]  William H. Press,et al.  Numerical Recipes 3rd Edition: The Art of Scientific Computing , 2007 .

[29]  J. Lynch,et al.  The turnover of organic carbon and nitrogen in soil. , 1990 .

[30]  Frank Nielsen,et al.  An Elementary Introduction to Information Geometry , 2018, Entropy.

[31]  H. Piaggio Differential Geometry of Curves and Surfaces , 1952, Nature.

[32]  S. Chib,et al.  Understanding the Metropolis-Hastings Algorithm , 1995 .

[33]  P. Reichert,et al.  Biogeochemical model of Lake Zürich : sensitivity, identifiability and uncertainty analysis , 2001 .

[34]  Dane Jordan Bjork,et al.  An Algorithm for Symbolic Computing of Singular Limits of Dynamical Systems , 2018 .

[35]  Otto Richter,et al.  Microbial growth in soil and nitrogen turnover : A theoretical model considering the activity state of microorganisms , 1998 .

[36]  C. Poll,et al.  Regulation of bacterial and fungal MCPA degradation at the soil–litter interface , 2010 .

[37]  Shu-lin Chen,et al.  Evaluation of a soil greenhouse gas emission model based on Bayesian inference and MCMC: Model uncertainty , 2013 .

[38]  B. Iooss,et al.  A Review on Global Sensitivity Analysis Methods , 2014, 1404.2405.

[39]  C. Poll,et al.  Micro-scale modeling of pesticide degradation coupled to carbon turnover in the detritusphere: model description and sensitivity analysis , 2013, Biogeochemistry.

[40]  H. Engl,et al.  Inverse problems in systems biology , 2009 .

[41]  O. Monga,et al.  Emergent Properties of Microbial Activity in Heterogeneous Soil Microenvironments: Different Research Approaches Are Slowly Converging, Yet Major Challenges Remain , 2018, Front. Microbiol..

[42]  O. Monga,et al.  Simulating microbial degradation of organic matter in a simple porous system using the 3-D diffusion-based model MOSAIC , 2013 .

[43]  Thomas J. Snowden,et al.  Methods of Model Reduction for Large-Scale Biological Systems: A Survey of Current Methods and Trends , 2017, Bulletin of Mathematical Biology.

[44]  A. V. Bruggen,et al.  Modeling wave-like dynamics of oligotrophic and copiotrophic bacteria along wheat roots in response to nutrient input from a growing root tip , 2005 .

[45]  Matthew J. Smith,et al.  Microbial models with data‐driven parameters predict stronger soil carbon responses to climate change , 2015, Global change biology.

[46]  Thomas Kätterer,et al.  Modeling coupled enzymatic and solute transport controls on decomposition in drying soils , 2016 .

[47]  Sabine Attinger,et al.  Multiresponse, multiobjective calibration as a diagnostic tool to compare accuracy and structural limitations of five coupled soil‐plant models and CLM3.5 , 2013 .

[48]  Steven D. Allison,et al.  The Michaelis–Menten kinetics of soil extracellular enzymes in response to temperature: a cross‐latitudinal study , 2012 .

[49]  Z Zi,et al.  Sensitivity analysis approaches applied to systems biology models. , 2011, IET systems biology.

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

[51]  Ricky Chachra,et al.  Structural susceptibility and separation of time scales in the van der Pol oscillator. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[53]  Ryan N. Gutenkunst,et al.  Sloppiness, Modeling, and Evolution in Biochemical Networks , 2007 .

[54]  C. Poll,et al.  Micro-scale modelling of carbon turnover driven by microbial succession at a biogeochemical interface , 2008 .

[55]  C. Poll,et al.  Modeling coupled pesticide degradation and organic matter turnover: From gene abundance to process rates , 2016 .

[56]  Mark K Transtrum,et al.  Geometry of nonlinear least squares with applications to sloppy models and optimization. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[57]  Mark K Transtrum,et al.  Why are nonlinear fits to data so challenging? , 2009, Physical review letters.

[58]  W. Nowak,et al.  A Primer for Model Selection: The Decisive Role of Model Complexity , 2018 .

[59]  Michael C. Abbott,et al.  Maximizing the information learned from finite data selects a simple model , 2017, Proceedings of the National Academy of Sciences.

[60]  H. Künsch,et al.  Practical identifiability analysis of large environmental simulation models , 2001 .

[61]  Enli Wang,et al.  Confidence in soil carbon predictions undermined by the uncertainties in observations and model parameterisation , 2016, Environ. Model. Softw..

[62]  Trevor Hastie,et al.  Statistical Learning with Sparsity: The Lasso and Generalizations , 2015 .

[63]  W. Riley,et al.  A total quasi-steady-state formulation of substrate uptake kinetics in complex networks and an example application to microbial litter decomposition , 2013 .

[64]  D. Makowski,et al.  To what extent do uncertainty and sensitivity analyses help unravel the influence of microscale physical and biological drivers in soil carbon dynamics models? , 2018, Ecological Modelling.

[65]  E. Davidson,et al.  Sensitivity of decomposition rates of soil organic matter with respect to simultaneous changes in temperature and moisture , 2015 .

[66]  George C. Verghese Getting to the gray box: Some challenges for model reduction , 2009, ACC.

[67]  J. Sethna,et al.  Parameter Space Compression Underlies Emergent Theories and Predictive Models , 2013, Science.

[68]  Gaurav Venkataraman,et al.  Identifiability, reducibility, and adaptability in allosteric macromolecules , 2017, The Journal of general physiology.

[69]  Xiaobo Zhou,et al.  Global Sensitivity Analysis , 2017, Encyclopedia of GIS.

[70]  C. Poll,et al.  Evidence for the importance of litter as a co-substrate for MCPA dissipation in an agricultural soil , 2016, Environmental Science and Pollution Research.

[71]  W. Parton,et al.  Analysis of factors controlling soil organic matter levels in Great Plains grasslands , 1987 .

[72]  Karline Soetaert,et al.  Inverse Modelling, Sensitivity and Monte Carlo Analysis in R Using Package FME , 2010 .

[73]  A. Franzluebbers,et al.  Soil Carbon and Nitrogen Mineralization , 2004 .

[74]  Francesca Pianosi,et al.  Global Sensitivity Analysis of environmental models: Convergence and validation , 2016, Environ. Model. Softw..

[75]  Zhongkui Luo,et al.  Uncertain future soil carbon dynamics under global change predicted by models constrained by total carbon measurements. , 2017, Ecological applications : a publication of the Ecological Society of America.

[76]  S. Allison,et al.  Evolutionary-Economic Principles as Regulators of Soil Enzyme Production and Ecosystem Function , 2010 .

[77]  Mehdi Gharasoo,et al.  Introduction of a new platform for parameter estimation of kinetically complex environmental systems , 2017, Environ. Model. Softw..

[78]  Matthew J. Smith,et al.  Oscillatory behavior of two nonlinear microbial models of soil carbon decomposition , 2013 .

[79]  Clemens Zarzer On Tikhonov regularization with non-convex sparsity constraints , 2009 .

[80]  Sean C. Warnick,et al.  A unified view of Balanced Truncation and Singular Perturbation Approximations , 2015, 2015 American Control Conference (ACC).

[81]  Michael T. Stobb,et al.  A local and global sensitivity analysis of a mathematical model of coagulation and platelet deposition under flow , 2018, PloS one.

[82]  S. Trumbore,et al.  An Uncertain Future for Soil Carbon , 2008, Science.

[83]  Li Zhang,et al.  Parameter identifiability, constraint, and equifinality in data assimilation with ecosystem models. , 2009, Ecological applications : a publication of the Ecological Society of America.

[84]  R. B. Jackson,et al.  Toward an ecological classification of soil bacteria. , 2007, Ecology.

[85]  P. Reed,et al.  Sensitivity-guided reduction of parametric dimensionality for multi-objective calibration of watershed models , 2009 .

[86]  André G. van Turnhout,et al.  A toolbox to find the best mechanistic model to predict the behavior of environmental systems , 2016, Environ. Model. Softw..

[87]  Mark K Transtrum,et al.  Model reduction by manifold boundaries. , 2014, Physical review letters.

[88]  Dario Vretenar,et al.  "Sloppy" nuclear energy density functionals: Effective model reduction , 2016, 1606.08617.

[89]  Eric A Davidson,et al.  Rate my data: quantifying the value of ecological data for the development of models of the terrestrial carbon cycle. , 2013, Ecological applications : a publication of the Ecological Society of America.

[90]  Andrew White,et al.  The Limitations of Model-Based Experimental Design and Parameter Estimation in Sloppy Systems , 2016, PLoS Comput. Biol..

[91]  Christopher R. Myers,et al.  Universally Sloppy Parameter Sensitivities in Systems Biology Models , 2007, PLoS Comput. Biol..

[92]  J. Harte,et al.  Microbial community-level regulation explains soil carbon responses to long-term litter manipulations , 2017, Nature Communications.

[93]  W. Rappel,et al.  Systematic reduction of a detailed atrial myocyte model. , 2017, Chaos.

[94]  Peng Qiu,et al.  Bridging Mechanistic and Phenomenological Models of Complex Biological Systems , 2015, PLoS Comput. Biol..

[95]  M. M. Mueller,et al.  Model structure and parameter identification of soil organic matter models , 2015 .

[96]  S. Billings,et al.  Model formulation of microbial CO2 production and efficiency can significantly influence short and long term soil C projections , 2018, The ISME Journal.

[97]  Pier Francesco Perri,et al.  Sensitivity analysis for environmental models , 2000 .

[98]  Model structures amplify uncertainty in predicted soil carbon responses to climate change , 2018, Nature Communications.

[99]  S. Blagodatsky,et al.  Microbial growth in soil and nitrogen turnover: Model calibration with laboratory data , 1998 .

[100]  Francesca Pianosi,et al.  A Matlab toolbox for Global Sensitivity Analysis , 2015, Environ. Model. Softw..

[101]  Hubert H. G. Savenije,et al.  Equifinality, a blessing in disguise? , 2001 .

[102]  Colin H. LaMont,et al.  Information-based inference in sloppy and singular models , 2015 .