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.

Abstract: Climate change scenarios suggest that long periods without rainfall will occur in the future often causing instability of the agricultural products market. The aim of our research was to build a model describing the amount of precipitation and droughts for forecasting crop yields in the future. In this study, we analysed a non-standard mixture of gamma and one point distributions as the model of rainfall. On the basis of the rainfall data, one can estimate parameters of the distribution. Parameter estimators were constructed using a method of maximum likelihood. The obtained rainfall data allow confirming the hypothesis of the adequacy of the proposed rainfall models. Long series of droughts allow one to determine the probabilities of adverse phenomena in agriculture. Based on the model, yields of barley in the years 2030 and 2050 were forecasted which can be used for the assessment of other crops productivity. The results obtained with this approach can be used to predict decreases in agricultural production caused by prospective rainfall shortages. This will enable decision makers to shape effective agricultural policies in order to learn how to balance the food supplies and demands through an appropriate management of stored raw food materials and import/export policies.

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.