Abstract: We compared the precision of simple random sampling (SimRS) and seven types of stratified random sampling (StrRS) schemes in estimating regional mean of water-limited yields for two crops (winter wheat and silage maize) that were simulated by fourteen crop models. We found that the precision gains of StrRS varied considerably across stratification methods and crop models. Precision gains for compact geographical stratification were positive, stable and consistent across crop models. Stratification with soil water holding capacity had very high precision gains for twelve models, but resulted in negative gains for two models. Increasing the sample size monotonously decreased the sampling errors for all the sampling schemes. We conclude that compact geographical stratification can modestly but consistently improve the precision in estimating regional mean yields. Using the most influential environmental variable for stratification can notably improve the sampling precision, especially when the sensitivity behavior of a crop model is known.

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.