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: 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: We show the error in water-limited yields simulated by crop models which is associated with spatially aggregated soil and climate input data. Crop simulations at large scales (regional, national, continental) frequently use input data of low resolution. Therefore, climate and soil data are often generated via averaging and sampling by area majority. This may bias simulated yields at large scales, varying largely across models. Thus, we evaluated the error associated with spatially aggregated soil and climate data for 14 crop models. Yields of winter wheat and silage maize were simulated under water-limited production conditions. We calculated this error from crop yields simulated at spatial resolutions from 1 to 100 km for the state of North Rhine-Westphalia, Germany. Most models showed yields biased by <15% when aggregating only soil data. The relative mean absolute error (rMAE) of most models using aggregated soil data was in the range or larger than the inter-annual or inter-model variability in yields. This error increased further when both climate and soil data were aggregated. Distinct error patterns indicate that the rMAE may be estimated from few soil variables. Illustrating the range of these aggregation effects across models, this study is a first step towards an ex-ante assessment of aggregation errors in large-scale simulations.