Multivariate global sensitivity analysis for dynamic crop models

Abstract Dynamic crop models are frequently used in ecology, agronomy and environmental sciences for simulating crop and environmental variables at a discrete time step. They often include a large number of parameters whose values are uncertain, and it is often impossible to estimate all these parameters accurately. A common practice consists in selecting a subset of parameters by global sensitivity analysis, estimating the selected parameters from data, and setting the others to some nominal values. For a discrete-time model, global sensitivity analyses can be applied sequentially at each simulation date. In the case of dynamic crop models, simulations are usually computed at a daily time step and the sequential implementation of global sensitivity analysis at each simulation date can result in several hundreds of sensitivity indices, with one index per parameter per simulation date. It is not easy to identify the most important parameters based on such a large number of values. In this paper, an alternative method called multivariate global sensitivity analysis was investigated. More precisely, the purposes of this paper are (i) to compare the sensitivity indices and associated parameter rankings computed by the sequential and the multivariate global sensitivity analyses, (ii) to assess the value of multivariate sensitivity analysis for selecting the model parameters to estimate from data. Sequential and multivariate sensitivity analyses were compared by using two dynamic models: a model simulating wheat biomass named WWDM and a model simulating N 2 O gaseous emission in crop fields named CERES-EGC. N 2 O measurements collected in several experimental plots were used to evaluate how parameter selection based on multivariate sensitivity analysis can improve the CERES-EGC predictions. The results showed that sequential and multivariate sensitivity analyses provide modellers with different types of information for models which exhibit a high variability of sensitivity index values over time. Conversely, when the parameter influence is quite constant over time, the two methods give more similar results. The results also showed that the estimation of the parameters with the highest sensitivity indices led to a strong reduction of the prediction errors of the model CERES-EGC.

[1]  Gilbert Saporta,et al.  Plans d'expériences : applications à l'entreprise , 1997 .

[2]  S. Pacala,et al.  Forest models defined by field measurements : Estimation, error analysis and dynamics , 1996 .

[3]  Stefano Tarantola,et al.  Sensitivity analysis practices: Strategies for model-based inference , 2006, Reliab. Eng. Syst. Saf..

[4]  B. Silverman,et al.  Functional Data Analysis , 1997 .

[5]  Andrea Saltelli,et al.  An effective screening design for sensitivity analysis of large models , 2007, Environ. Model. Softw..

[6]  J. Germon,et al.  NEMIS, a predictive model of denitrification on the field scale , 2000 .

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

[8]  David Makowski,et al.  Global sensitivity analysis for calculating the contribution of genetic parameters to the variance of crop model prediction , 2006, Reliab. Eng. Syst. Saf..

[9]  David Makowski,et al.  Multivariate global sensitivity analysis for discre te-time models. , 2008 .

[10]  James W. Jones,et al.  POTENTIAL USES AND LIMITATIONS OF CROP MODELS , 1996 .

[11]  Hiromi Seno,et al.  Ecological balance in the native population dynamics may cause the paradox of pest control with harvesting. , 2008, Journal of theoretical biology.

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

[13]  T. W. Anderson.  An Introduction to Multivariate Statistical Analysis , 1959 .

[14]  Daniel Wallach,et al.  Parameter Estimation for Crop Models , 2001 .

[15]  Ludek Berec,et al.  Modeling population dynamics of two cockroach species: effects of the circadian clock, interspecific competition and pest control. , 2007, Journal of theoretical biology.

[16]  Luca Bechini,et al.  Parameterization of a crop growth and development simulation model at sub-model components level. An example for winter wheat (Triticum aestivum L.) , 2006, Environ. Model. Softw..

[17]  A. Saltelli,et al.  The role of sensitivity analysis in ecological modelling , 2007 .

[18]  Sabrina Gaba,et al.  Combined use of local and ANOVA-based global sensitivity analyses for the investigation of a stochastic dynamic model: Application to the case study of an individual-based model of a fish population , 2006 .

[19]  William N. Venables,et al.  Modern Applied Statistics with S , 2010 .

[20]  Stefano Tarantola,et al.  Sensitivity Analysis as an Ingredient of Modeling , 2000 .

[21]  Stefano Tarantola,et al.  Sensitivity Analysis in Practice , 2002 .

[22]  K. Butterbach-Bahl,et al.  Quantifying the regional source strength of N-trace gases across agricultural and forest ecosystems with process based models , 2004, Plant and Soil.

[23]  George E. P. Box,et al.  Empirical Model‐Building and Response Surfaces , 1988 .

[24]  Brian J. Williams,et al.  Sensitivity analysis when model outputs are functions , 2006, Reliab. Eng. Syst. Saf..

[25]  D. Wallach,et al.  Comparison of parameter estimation methods for crop models , 2004 .

[26]  Benoit Gabrielle,et al.  Process‐based modeling of nitrous oxide emissions from wheat‐cropped soils at the subregional scale , 2006 .

[27]  J. Passioura.  Simulation Models: Science, Snake Oil, Education, or Engineering? , 1996 .

[28]  David Makowski,et al.  Bayesian calibration of the nitrous oxide emission module of an agro-ecosystem model , 2009 .

[29]  A. Saltelli,et al.  Sensitivity Anaysis as an Ingredient of Modeling , 2000 .

[30]  B. Gabrielle,et al.  Simulation of Nitrous Oxide Emissions from Wheat-cropped Soils using CERES , 2006, Nutrient Cycling in Agroecosystems.

[31]  Benoit Gabrielle,et al.  Predicting in situ soil N2O emission using NOE algorithm and soil database , 2005 .

[32]  Mimmo Iannelli,et al.  Population dynamics and conservation biology of the over-exploited Mediterranean red coral. , 2007, Journal of theoretical biology.

[33]  Carlos Castillo-Chavez,et al.  Interplay between local dynamics and dispersal in discrete-time metapopulation models. , 2002, Journal of theoretical biology.