Abstract: Working with ensembles of crop models is a recent but important development in crop modeling which promises to lead to better uncertainty estimates for model projections and predictions, better predictions using the ensemble mean or median, and closer collaboration within the modeling community. There are numerous open questions about the best way to create and analyze such ensembles. Much can be learned from the field of climate modeling, given its much longer experience with ensembles. We draw on that experience to identify questions and make propositions that should help make ensemble modeling with crop models more rigorous and informative. The propositions include defining criteria for acceptance of models in a crop MME, exploring criteria for evaluating the degree of relatedness of models in a MME, studying the effect of number of models in the ensemble, development of a statistical model of model sampling, creation of a repository for MME results, studies of possible differential weighting of models in an ensemble, creation of single model ensembles based on sampling from the uncertainty distribution of parameter values or inputs specifically oriented toward uncertainty estimation, the creation of super ensembles that sample more than one source of uncertainty, the analysis of super ensemble results to obtain information on total uncertainty and the separate contributions of different sources of uncertainty and finally further investigation of the use of the multi-model mean or median as a predictor.

Abstract: Interactions between soil properties and climate affect forage grass productivity. Dynamic models, simulating crop performance as a function of environmental conditions, are valid for a specific location with given soil and weather conditions. Extrapolations of local soil properties to larger regions can help assess the requirement for soil input in regional yield estimations. Using the LINGRA model, we simulated the regional yield level and variability of timothy, a forage grass, in Akershus and Ostfold counties, Norway. Soils were grouped according to physical similarities according to 4 sets of criteria. This resulted in 66, 15, 5 and 1 groups of soils. The properties of the soil with the largest area was extrapolated to the other soils within each group and input to the simulations. All analyses were conducted for 100 yr of generated weather representing the period 1961-1990, and climate projections for the period 2046-2065, the Intergovernmental Panel on Climate Change greenhouse gas emission scenario A1B, and 4 global climate models. The simulated regional seasonal timothy yields were 5-13% lower on average and had higher inter-annual variability for the least detailed soil extrapolation than for the other soil extrapolations, across climates. There were up to 20% spatial intra-regional differences in simulated yield between soil extrapolations. The results indicate that, for conditions similar to these studied here, a few representative profiles are sufficient for simulations of average regional seasonal timothy yield. More spatially detailed yield analyses would benefit from more detailed soil input.

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.