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.

Abstract: As a part of the MACSUR task L2.4, we probabilistically calibrated the CARAIB dynamic vegetation model by Markov chain Monte Carlo (MCMC) simulation with the DREAMZS sampler. CARAIB is a mechanistic model that calculates the carbon assimilation of the vegetation as a function of the soil and climatic conditions, and can thus be used for simulating grassland production under cutting or grazing management. Bayesian model inversion was performed at 4 grassland sites across Europe: Oensingen, CH; Grillenburg, DE; Laqueuille, FR and Monte-Bodone, IT. Four daily measured variables from these sites: the Gross Primary Productivity (GPP), Net Ecosystem Exchange (NEE), Evapotranspiration (ET) and Soil Water Content (SWC) were used to sample 10 parameters related to rooting depth, stomatal conductance, specific leaf area, carbon-nitrogen ratio and water stresses. The maximized likelihood function therefore involved four objectives, whereas the applied Bayesian framework allowed for assessing the so called parameter posterior probability density function (pdf), which quantifies model parameter uncertainty caused by measurement and model errors. Sampling trials were performed using merged data from all sites (all-sites-sampling) and for each site (site-specific sampling) separately. The derived posterior parameter pdfs from the all-sites sampling and site-specific sampling runs showed differences in relation with the specificities of each site. Analysis of these distributions also revealed model sensitivity to parameters conditioned on the measured data, as well as parameter correlations.