Abstract: An environmental dynamic system is usually modeled as a nonlinear system described by a set of nonlinear ODEs. A central challenge in computational modeling of environmental systems is the determination of the model parameters. In these cases, estimating these variables or parameters from other easily obtained measurements can be extremely useful. This work addresses the problem of monitoring and modeling a leaf area index and soil moisture model (LSM) using state estimation. The performances of various conventional and state-of-the-art state estimation techniques are compared when they are utilized to achieve this objective. These techniques include the extended Kalman filter (EKF), particle filter (PF), and the more recently developed technique variational filter (VF). Specifically, two comparative studies are performed. In the first comparative study, the state variables (the leaf-area index LAI, the volumetric water content of the soil layer 1, HUR1 and the volumetric water content of the soil layer 2, HUR2) are estimated from noisy measurements of these variables, and the various estimation techniques are compared by computing the estimation root mean square error (RMSE) with respect to the noise-free data. In the second comparative study, the state variables as well as the model parameters are simultaneously estimated. In this case, in addition to comparing the performances of the various state estimation techniques, the effect of number of estimated model parameters on the accuracy and convergence of these techniques are also assessed. The results of both comparative studies show that the PF provides a higher accuracy than the EKF, which is due to the limited ability of the EKF to handle highly nonlinear processes. The results also show that the VF provides a significant improvement over the PF because, unlike the PF which depends on the choice of sampling distribution used to estimate the posterior distribution, the VF yields an optimum choice of the sampling distribution, which also accounts for the observed data. The results of the second comparative study show that, for all techniques, estimating more model parameters affects the estimation accuracy as well as the convergence of the estimated states and parameters. However, the VF can still provide both convergence as well as accuracy related advantages over other estimation methods. (C) 2013 Elsevier B.V. All rights reserved.

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: Increasing cereal yield is needed to meet the projected increased demand for world food supply of about 70% by 2050. Sirius, a process-based model for wheat, was used to estimate yield potential for wheat ideotypes optimized for future climatic projections for ten wheat growing areas of Europe. It was predicted that the detrimental effect of drought stress on yield would be decreased due to enhanced tailoring of phenology to future weather patterns, and due to genetic improvements in the response of photosynthesis and green leaf duration to water shortage. Yield advances could be made through extending maturation and thereby improve resource capture and partitioning. However the model predicted an increase in frequency of heat stress at meiosis and anthesis. Controlled environment experiments quantify the effects of heat and drought at booting and flowering on grain numbers and potential grain size. A current adaptation of wheat to areas of Europe with hotter and drier summers is a quicker maturation which helps to escape from excessive stress, but results in lower yields. To increase yield potential and to respond to climate change, increased tolerance to heat and drought stress should remain priorities for the genetic improvement of wheat.