Abstract: During the past few years, several studies have compared the performance of crop simulation models to assess the uncertainties in model-based climate change impact assessments and other modelling studies. Many of these studies have concentrated on cereal crops, while fewer model comparisons have been conducted for grasses. We compared the predictions for timothy grass (Phleum pratertse L.) yields for first and second cuts along with the dynamics of above-ground biomass for the grass simulation models BASGRA and CATIMO, and the soil -crop model STICS. The models were calibrated and evaluated using field data from seven sites across Northern Europe and Canada with different climates, soil conditions and management practices. Altogether the models were compared using data on timothy grass from 33 combinations of sites, cultivars and management regimes. Model performances with two calibration approaches, cultivar-specific and generic calibrations, were compared. All the models studied estimated the dynamics of above-ground biomass and the leaf area index satisfactorily, but tended to underestimate the first cut yield. Cultivar-specific calibration resulted in more accurate first cut yield predictions than the generic calibration achieving root mean square errors approximately one third lower for the cultivar-specific calibration. For the second cut, the difference between the calibration methods was small. The results indicate that detailed soil process descriptions improved the overall model performance and the model responses to management, such as nitrogen applications. The results also suggest that taking the genetic variability into account between cultivars of timothy grass also improves the yield estimates. Calibrations using both spring and summer growth data simultaneously revealed that processes determining the growth in these two periods require further attention in model development.

Abstract: Eleven widely used crop simulation models (APSIM, CERES, CROPSYST, COUP, DAISY, EPIC, FASSET, HERMES, MONICA, STICS and WOFOST) were tested using spring barley (Hordeum vulgare L.) data set under varying nitrogen (N) fertilizer rates from three experimental years in the boreal climate of Jokioinen, Finland. This is the largest standardized crop model inter-comparison under different levels of N supply to date. The models were calibrated using data from 2002 and 2008, of which 2008 included six N rates ranging from 0 to 150 kg N/ha. Calibration data consisted of weather, soil, phenology, leaf area index (LAI) and yield observations. The models were then tested against new data for 2009 and their performance was assessed and compared with both the two calibration years and the test year. For the calibration period, root mean square error between measurements and simulated grain dry matter yields ranged from 170 to 870 kg/ha. During the test year 2009, most models failed to accurately reproduce the observed low yield without N fertilizer as well as the steep yield response to N applications. The multi-model predictions were closer to observations than most single-model predictions, but multi-model mean could not correct systematic errors in model simulations. Variation in soil N mineralization and LAI development due to differences in weather not captured by the models most likely was the main reason for their unsatisfactory performance. This suggests the need for model improvement in soil N mineralization as a function of soil temperature and moisture. Furthermore, specific weather event impacts such as low temperatures after emergence in 2009, tending to enhance tillering, and a high precipitation event just before harvest in 2008, causing possible yield penalties, were not captured by any of the models compared in the current study.

Abstract: Field-scale crop models are often applied at spatial resolutions coarser than that of the arable field. However, little is known about the response of the models to spatially aggregated climate input data and why these responses can differ across models. Depending on the model, regional yield estimates from large-scale simulations may be biased, compared to simulations with high-resolution input data. We evaluated this so-called aggregation effect for 13 crop models for the region of North Rhine-Westphalia in Germany. The models were supplied with climate data of 1 km resolution and spatial aggregates of up to 100 km resolution raster. The models were used with 2 crops (winter wheat and silage maize) and 3 production situations (potential, water-limited and nitrogen-water-limited growth) to improve the understanding of errors in model simulations related to data aggregation and possible interactions with the model structure. The most important climate variables identified in determining the model-specific input data aggregation on simulated yields were mainly related to changes in radiation (wheat) and temperature (maize). Additionally, aggregation effects were systematic, regardless of the extent of the effect. Climate input data aggregation changed the mean simulated regional yield by up to 0.2 t ha(-1), whereas simulated yields from single years and models differed considerably, depending on the data aggregation. This implies that large-scale crop yield simulations are robust against climate data aggregation. However, large-scale simulations can be systematically biased when being evaluated at higher temporal or spatial resolution depending on the model and its parameterization.