Abstract: Process-based models (PBM) for simulation of weather dependent grass growth can assist farmers and plant breeders in addressing the challenges of climate change by simulating alternative roads of adaptation. They can also provide management decision support under current conditions. A drawback of existing grass models is that they do not take into account the effect of winter stresses, limiting their use for full-year simulations in areas where winter survival is a key factor for yield security. Here, we present a novel full-year PBM for grassland named BASGRA. It was developed by combining the LINGRA grassland model (Van Oijen et al., 2005a) with models for cold hardening and soil physical winter processes. We present the model and show how it was parameterized for timothy (Phleum pratense L.), the most important forage grass in Scandinavia and parts of North America and Asia. Uniquely, BASGRA simulates the processes taking place in the sward during the transition from summer to winter, including growth cessation and gradual cold hardening, and functions for simulating plant injury due to low temperatures, snow and ice affecting regrowth in spring. For the calibration, we used detailed data from five different locations in Norway, covering a wide range of agroclimatic regions, day lengths (latitudes from 59 degrees to 70 degrees N) and soil conditions. The total dataset included 11 variables, notably above-ground dry matter, leaf area index, tiller density, content of C reserves, and frost tolerance. All data were used in the calibration. When BASGRA was run with the maximum a-posteriori (MAP) parameter vector from the single, Bayesian calibration, nearly all measured variables were simulated to an overall normalized root mean squared error (NRMSE) <0.5. For many site x experiment combinations, NRMSE was <0.3. The temporal dynamics were captured well for most variables, as evaluated by comparing simulated time courses versus data for the individual sites. The results may suggest that BASGRA is a reasonably robust model, allowing for simulation of growth and several important underlying processes with acceptable accuracy for a range of agroclimatic conditions. However, the robustness of the model needs to be tested further using independent data from a wide range of growing conditions. Finally we show an example of application of the model, comparing overwintering risks in two climatically different sites, and discuss future model applications. Further development work should include improved simulation of the dynamics of C reserves, and validation of winter tiller dynamics against independent data. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: Rice models have been widely used in simulating and predicting rice phenology in contrasting climate zones, however the uncertainties from model structure (different equations or models) and/or model parameters were rarely investigated. Here, five rice phenological models/modules (Le., CERES-Rice, ORYZA2000, RCM, Beta Model and SIMRIW) were applied to simulate rice phenology at 23 experimental stations from 1992 to 2009 in two major rice cultivation regions of China: the northeastern China and the southwestern China. To investigate the uncertainties from model biophysical parameters, each model was run with randomly perturbed 50 sets of parameters. The results showed that the median of ensemble simulations were better than the simulation by most models. Models couldn’t simulate well in some specific years despite of parameters optimization, suggesting model structure limit model performance in some cases. The models adopting accumulative thermal time function (e.g., CERES-Rice and ORYZA2000) had better performance in the southwestern China, in contrast, those adopting exponential function (e.g., Beta model and RCM model) had better performance in the northeastern China. In northeastern China, the contribution of model structure and model parameters to model total variance was, respectively, about 55.90% and 44.10% in simulating heading date, and about 75.43% and 24.57% in simulating maturity date. In the southwestern China, the contribution of model structure and model parameters to model total variance was, respectively, about 79.97% and 27.03% in simulating heading date, about 92.15% and 7.85% in simulating maturity date. Uncertainty from model structure was the most relevant source. The results highlight that the temperature response functions of rice development rate under extreme climate conditions should be improved based on environment-controlled experimental data.

Abstract: • Systematic analysis of the influence of spatial soil data resolution on simulated regional yields and total growing season evapotranspiration. • The responses of four crop models of different complexity are compared. • Differences between models are larger than the effect of the chosen spatial soil data resolution. • Low influence of soil data resolution due to: high precipitation amount, methods for calculating water retention and method of data aggregation. The spatial variability of soil properties is an important driver of yield variability at both field and regional scale. Thus, when using crop growth simulation models, the choice of spatial resolution of soil input data might be key in order to accurately reproduce observed yield variability. In this study we used four crop models (SIMPLACE<LINTUL-SLIM>, DSSAT-CSM, EPIC and DAISY) differing in the detail of modeling above-ground biomass and yield as well as of modeling soil water dynamics, water uptake and drought effects on plants to simulate winter wheat in two (agro-climatologically and geo-morphologically) contrasting regions of the federal state of North-Rhine-Westphalia (Germany) for the period from 1995 to 2008. Three spatial resolutions of soil input data were taken into consideration, corresponding to the following map scales: 1:50 000, 1:300 000 and 1:1 000 000. The four crop models were run for water-limited production conditions and model results were evaluated in the form of frequency distributions, depicted by bean-plots. In both regions, soil data aggregation had very small influence on the shape and range of frequency distributions of simulated yield and simulated total growing season evapotranspiration for all models. Further analysis revealed that the small influence of spatial resolution of soil input data might be related to: (a) the high precipitation amount in the region which partly masked differences in soil characteristics for water holding capacity, (b) the loss of variability in hydraulic soil properties due to the methods applied to calculate water retention properties of the used soil profiles, and (c) the method of soil data aggregation. No characteristic “fingerprint” between sites, years and resolutions could be found for any of the models. Our results support earlier recommendation to evaluate model results on the basis of frequency distributions since these offer quick and better insight into the distribution of simulation results as compared to summary statistics only. Finally, our results support conclusions from other studies about the usefulness of considering a multi-model approach to quantify the uncertainty in simulated yields introduced by the crop growth simulation approach when exploring the effects of scaling for regional yield impact assessments.