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.

Abstract: Heat and drought stress can reduce crop yields considerably which is increasingly assessed with crop models for larger areas. Applying these models originally developed for the field scale at large spatial extent typically implies the use of input data with coarse resolution. Little is known about the effect of data resolution on the simulated impact of extreme events like heat and drought on crops. Hence, in this study the effect of input and output data aggregation on simulated heat and drought stress and their impact on yield of winter wheat is systematically analyzed. The crop model SIMPLACE was applied for the period 1980-2011 across Germany at a resolution of 1 km x 1 km. Weather and soil input data and model output data were then aggregated to 10 km x 10 km, 25 km x 25 km, 50 km x 50 km and 100 km x 100 km resolution to analyze the aggregation effect on heat and drought stress and crop yield. We found that aggregation of model input and output data barely influenced the mean and median of heat and drought stress reduction factors and crop yields simulated across Germany. However, data aggregation resulted in less spatial variability of model results and a reduced severity of simulated stress events, particularly for regions with high heterogeneity in weather and soil conditions. Comparisons of simulations at coarse resolution with those at high resolution showed distinct patterns of positive and negative deviations which compensated each other so that aggregation effects for large regions were small for mean or median yields. Therefore, modelling at a resolution of 100 km x 100 km was sufficient to determine mean wheat yield as affected by heat and drought stress for Germany. Further research is required to clarify whether the results can be generalized across crop models differing in structure and detail. Attention should also be given to better understand the effect of data resolution on interactions between heat and drought impacts. (C) 2015 Elsevier B.V. All rights reserved.

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.