|
Dumont, B., Leemans, V., Mansouri, M., Bodson, B., Destain, J. - P., & Destain, M. - F. (2014). Parameter identification of the STICS crop model, using an accelerated formal MCMC approach. Env. Model. Softw., 52, 121–135.
Abstract: This study presents a Bayesian approach for the parameters’ identification of the STICS crop model based on the recently developed Differential Evolution Adaptive Metropolis (DREAM) algorithm. The posterior distributions of nine specific crop parameters of the STICS model were sampled with the aim to improve the growth simulations of a winter wheat (Triticum aestivum L) culture. The results obtained with the DREAM algorithm were initially compared to those obtained with a Nelder-Mead Simplex algorithm embedded within the OptimiSTICS package. Then, three types of likelihood functions implemented within the DREAM algorithm were compared, namely the standard least square, the weighted least square, and a transformed likelihood function that makes explicit use of the coefficient of variation (CV). The results showed that the proposed CV likelihood function allowed taking into account both noise on measurements and heteroscedasticity which are regularly encountered in crop modelling. (C) 2013 Elsevier Ltd. All rights reserved.
|
|
|
Mansouri, M., & Destain, M. - F. (2015). Predicting biomass and grain protein content using Bayesian methods. Stoch. Environ. Res. Risk Assess., 29(4), 1167–1177.
Abstract: This paper deals with the problem of predicting biomass and grain protein content using improved particle filtering (IPF) based on minimizing the Kullback-Leibler divergence. The performances of IPF are compared with those of the conventional particle filtering (PF) in two comparative studies. In the first one, we apply IPF and PF at a simple dynamic crop model with the aim to predict a single state variable, namely the winter wheat biomass, and to estimate several model parameters. In the second study, the proposed IPF and the PF are applied to a complex crop model (AZODYN) to predict a winter-wheat quality criterion, namely the grain protein content. The results of both comparative studies reveal that the IPF method provides a better estimation accuracy than the PF method. The benefit of the IPF method lies in its ability to provide accuracy related advantages over the PF method since, unlike the PF which depends on the choice of the sampling distribution used to estimate the posterior distribution, the IPF yields an optimum choice of this sampling distribution, which also utilizes the observed data. The performance of the proposed method is evaluated in terms of estimation accuracy, root mean square error, mean absolute error and execution times.
|
|
|
Conradt, T., Gornott, C., & Wechsung, F. (2016). Extending and improving regionalized winter wheat and silage maize yield regression models for Germany: Enhancing the predictive skill by panel definition through cluster analysis. Agricultural and Forest Meteorology, 216, 68–81.
Abstract: Regional agricultural yield assessments allowing for weather effect quantifications are a valuable basis for deriving scenarios of climate change effects and developing adaptation strategies. Assessing weather effects by statistical methods is a classical approach, but for obtaining robust results many details deserve attention and require individual decisions as is demonstrated in this paper. We evaluated regression models for annual yield changes of winter wheat and silage maize in more than 300 German counties and revised them to increase their predictive power. A major effort of this study was, however, aggregating separately estimated time series models (STSM) into panel data models (PDM) based on cluster analyses. The cluster analyses were based on the per-county estimates of STSM parameters. The original STSM formulations (adopted from a parallel study) contained also the non-meteorological input variables acreage and fertilizer price. The models were revised to use only weather variables as estimation basis. These consisted of time aggregates of radiation, precipitation, temperature, and potential evapotranspiration. Altering the input variables generally increased the predictive power of the models as did their clustering into PDM. For each crop, five alternative clusterings were produced by three different methods, and similarities between their spatial structures seem to confirm the existence of objective clusters about common model parameters. Observed smooth transitions of STSM parameter values in space suggest, however, spatial autocorrelation effects that could also be modeled explicitly. Both clustering and autocorrelation approaches can effectively reduce the noise in parameter estimation through targeted aggregation of input data. (C) 2015 Elsevier B.V. All rights reserved.
|
|
|
Mansouri, M., Dumont, B., Leemans, V., & Destain, M. - F. (2014). Bayesian methods for predicting LAI and soil water content. Precision Agric., 15(2), 184–201.
Abstract: LAI of winter wheat (Triticum aestivum L.) and soil water content of the topsoil (200 mm) and of the subsoil (500 mm) were considered as state variables of a dynamic soil-crop system. This system was assumed to progress according to a Bayesian probabilistic state space model, in which real values of LAI and soil water content were daily introduced in order to correct the model trajectory and reach better future evolution. The chosen crop model was mini STICS which can reduce the computing and execution times while ensuring the robustness of data processing and estimation. To predict simultaneously state variables and model parameters in this non-linear environment, three techniques were used: extended Kalman filtering (EKF), particle filtering (PF), and variational filtering (VF). The significantly improved performance of the VF method when compared to EKF and PF is demonstrated. The variational filter has a low computational complexity and the convergence speed of states and parameters estimation can be adjusted independently. Detailed case studies demonstrated that the root mean square error of the three estimated states (LAI and soil water content of two soil layers) was smaller and that the convergence of all considered parameters was ensured when using VF. Assimilating measurements in a crop model allows accurate prediction of LAI and soil water content at a local scale. As these biophysical properties are key parameters in the crop-plant system characterization, the system has the potential to be used in precision farming to aid farmers and decision makers in developing strategies for site-specific management of inputs, such as fertilizers and water irrigation.
|
|
|
Minet, J., Laloy, E., Tychon, B., & François, L. (2015). Bayesian inversions of a dynamic vegetation model at four European grassland sites. Biogeosciences, 12(9), 2809–2829.
Abstract: Eddy covariance data from four European grassland sites are used to probabilistically invert the CARAIB (CARbon Assimilation In the Biosphere) dynamic vegetation model (DVM) with 10 unknown parameters, using the DREAM((ZS)) (DiffeRential Evolution Adaptive Metropolis) Markov chain Monte Carlo (MCMC) sampler. We focus on comparing model inversions, considering both homoscedastic and heteroscedastic eddy covariance residual errors, with variances either fixed a priori or jointly inferred together with the model parameters. Agreements between measured and simulated data during calibration are comparable with previous studies, with root mean square errors (RMSEs) of simulated daily gross primary productivity (GPP), ecosystem respiration (RECO) and evapotranspiration (ET) ranging from 1.73 to 2.19, 1.04 to 1.56 g C m(-2) day(-1) and 0.50 to 1.28 mm day(-1), respectively. For the calibration period, using a homoscedastic eddy covariance residual error model resulted in a better agreement between measured and modelled data than using a heteroscedastic residual error model. However, a model validation experiment showed that CARAIB models calibrated considering heteroscedastic residual errors perform better. Posterior parameter distributions derived from using a heteroscedastic model of the residuals thus appear to be more robust. This is the case even though the classical linear heteroscedastic error model assumed herein did not fully remove heteroscedasticity of the GPP residuals. Despite the fact that the calibrated model is generally capable of fitting the data within measurement errors, systematic bias in the model simulations are observed. These are likely due to model inadequacies such as shortcomings in the photosynthesis modelling. Besides the residual error treatment, differences between model parameter posterior distributions among the four grassland sites are also investigated. It is shown that the marginal distributions of the specific leaf area and characteristic mortality time parameters can be explained by site-specific ecophysiological characteristics.
|
|