Abstract: Crop growth simulation models are increasingly used for regionally assessing the effects of climate change and variability on crop yields. These models require spatially and temporally detailed, location-specific, environmental (weather and soil) and management data as inputs, which are often difficult to obtain consistently for larger regions. Aggregating the resolution of input data for crop model applications may increase the uncertainty of simulations to an extent that is not well understood. The present study aims to systematically analyse the effect of changes in the spatial resolution of weather input data on yields simulated by four crop models (LINTUL-SLIM, DSSAT-CSM, EPIC and WOFOST) which were utilized to test possible interactions between weather input data resolution and specific modelling approaches representing different degrees of complexity. The models were applied to simulate grain yield of spring barley in Finland for 12 years between 1994 and 2005 considering five spatial resolutions of daily weather data: weather station (point) and grid-based interpolated data at resolutions of 10 km x 10 km; 20 km x 20 km; 50 km x 50 km and 100 km x 100 km. Our results show that the differences between models were larger than the effect of the chosen spatial resolution of weather data for the considered years and region. When displaying model results graphically, each model exhibits a characteristic ‘fingerprint’ of simulated yield frequency distributions. These characteristic distributions in response to the inter-annual weather variability were independent of the spatial resolution of weather input data. Using one model (LINTUL-SLIM), we analysed how the aggregation strategy, i.e. aggregating model input versus model output data, influences the simulated yield frequency distribution. Results show that aggregating weather data has a smaller effect on the yield distribution than aggregating simulated yields which causes a deformation of the model fingerprint. We conclude that changes in the spatial resolution of weather input data introduce less uncertainty to the simulations than the use of different crop models but that more evaluation is required for other regions with a higher spatial heterogeneity in weather conditions, and for other input data related to soil and crop management to substantiate our findings. Our results provide further evidence to support other studies stressing the importance of using not just one, but different crop models in climate assessment studies. (C) 2013 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.

Abstract: In this study, the performance of nine widely used and accessible crop growth simulation models (APES-ACE, CROPSYST, DAISY, DSSAT-CERES, FASSET, HERMES, MONICA, STICS and WOFOST) was compared during 44 growing seasons of spring barley (Hordeum vulgare L) at seven sites in Northern and Central Europe. The aims of this model comparison were to examine how different process-based crop models perform at multiple sites across Europe when applied with minimal information for model calibration of spring barley at field scale, whether individual models perform better than the multi-model mean, and what the uncertainty ranges are in simulated grain yields. The reasons for differences among the models and how results for barley compare to winter wheat are discussed. Regarding yield estimation, best performing based on the root mean square error (RMSE) were models HERMES, MONICA and WOFOST with lowest values of 1124, 1282 and 1325 (kg ha(-1)), respectively. Applying the index of agreement (IA), models WOFOST, DAISY and HERMES scored best having highest values (0.632, 0.631 and 0.585, respectively). Most models systematically underestimated yields, whereby CROPSYST showed the highest deviation as indicated by the mean bias error (MBE) (-1159 kg ha(-1)). While the wide range of simulated yields across all sites and years shows the high uncertainties in model estimates with only restricted calibration, mean predictions from the nine models agreed well with observations. Results of this paper also show that models that were more accurate in predicting phenology were not necessarily the ones better estimating grain yields. Total above-ground biomass estimates often did not follow the patterns of grain yield estimates and, thus, harvest indices were also different. Estimates of soil moisture dynamics varied greatly. In comparison, even though the growing cycle for winter wheat is several months longer than for spring barley, using RMSE and IA as indicators, models performed slightly, but not significantly, better in predicting wheat yields. Errors in reproducing crop phenology were similar, which in conjunction with the shorter growth cycle of barley has higher effects on accuracy in yield prediction. (C) 2012 Elsevier B.V. All rights reserved.