Abstract: Soil water percolation is a key process in the life cycle of water in fields, watersheds, and river basins. The Environmental Policy Integrated Climate (EPIC) and the Agricultural Policy/Environmental eXtender (APEX) are continuous models developed for evaluating the environmental effects of agricultural management. Traditionally, these models have simulated soil water percolation processes using a tipping-bucket approach, with the rate of flow limited by the saturated hydraulic conductivity. This simple approach often leads to inaccuracy in simulating elevated soil water conditions where soil water content (SWC) levels may remain above field capacity under prolonged wet weather periods or limited drainage. To overcome this deficiency, a new sub-model, the variable saturation hydraulic conductivity (VSHC) method, was developed for simulating soil water percolation processes using a nonlinear equation to estimate the effective hydraulic conductivity as a function of the SWC and soil properties. The VSHC method was evaluated at three sites in the United States and two sites in Europe. In addition, a numerical solution of the Richards equation was used as a benchmark for SWC comparison. Results show that the VSHC method substantially improves the accuracy of the SWC simulation in long-term simulations, particularly during wet periods. At the watershed scale, results on the Riesel Y2 watershed indicate that the VSHC method enhances model performance in the high-flow regime of channel peak flows because of the improved estimation of SWC, which implies that the improved SWC simulation at the field scale is beneficial to hydrologic modeling at the watershed scale.

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: ► We compared nine crop simulation models for spring barley at seven sites in Europe. ► Applying crop models with restricted calibration leads to high uncertainties. ► Multi-crop model mean yield estimates were in good agreement with observations. ► The degree of uncertainty for simulated grain yield of barley was similar to winter wheat. ► We need more suitable data enabling us to verify different processes in the models. In this study, the performance of nine widely used and accessible crop growth simulation models (APES-ACE, CROPSYST, DAISY, DSSAT-CERES, FASSET, HERMES, MONICA, STICS and WOFOST) was compared during 44 growing seasons of spring barley (Hordeum vulgare L) at seven sites in Northern and Central Europe. The aims of this model comparison were to examine how different process-based crop models perform at multiple sites across Europe when applied with minimal information for model calibration of spring barley at field scale, whether individual models perform better than the multi-model mean, and what the uncertainty ranges are in simulated grain yields. The reasons for differences among the models and how results for barley compare to winter wheat are discussed. Regarding yield estimation, best performing based on the root mean square error (RMSE) were models HERMES, MONICA and WOFOST with lowest values of 1124, 1282 and 1325 (kg ha(-1)), respectively. Applying the index of agreement (IA), models WOFOST, DAISY and HERMES scored best having highest values (0.632, 0.631 and 0.585, respectively). Most models systematically underestimated yields, whereby CROPSYST showed the highest deviation as indicated by the mean bias error (MBE) (-1159 kg ha(-1)). While the wide range of simulated yields across all sites and years shows the high uncertainties in model estimates with only restricted calibration, mean predictions from the nine models agreed well with observations. Results of this paper also show that models that were more accurate in predicting phenology were not necessarily the ones better estimating grain yields. Total above-ground biomass estimates often did not follow the patterns of grain yield estimates and, thus, harvest indices were also different. Estimates of soil moisture dynamics varied greatly. In comparison, even though the growing cycle for winter wheat is several months longer than for spring barley, using RMSE and IA as indicators, models performed slightly, but not significantly, better in predicting wheat yields. Errors in reproducing crop phenology were similar, which in conjunction with the shorter growth cycle of barley has higher effects on accuracy in yield prediction.