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.