Abstract: Highlights • The parameters to be fixed were consistent across sites. • Model calibration must be performed separately for each specific case. • Possible to reduce model parameters from 66 to 45. • Strong model reductions must be avoided. • The error term for the training data were characterised by timing (phase shift). Abstract Proper parameterisation and quantification of model uncertainty are two essential tasks in improvement and assessment of model performance. Bayesian calibration is a method that combines both tasks by quantifying probability distributions for model parameters and outputs. However, the method is rarely applied to complex models because of its high computational demand when used with high-dimensional parameter spaces. We therefore combined Bayesian calibration with sensitivity analysis, using the screening method by Morris (1991), in order to reduce model complexity by fixing parameters to which model output was only weakly sensitive to a nominal value. Further, the robustness of the model with respect to reduction in the number of free parameters were examined according to model discrepancy and output uncertainty. The process-based grassland model BASGRA was examined in the present study on two sites in Norway and in Germany, for two grass species (Phleum pratense and Arrhenatherum elatius). According to this study, a reduction of free model parameters from 66 to 45 was possible. The sensitivity analysis showed that the parameters to be fixed were consistent across sites (which differed in climate and soil conditions), while model calibration had to be performed separately for each combination of site and species. The output uncertainty decreased slightly, but still covered the field observations of aboveground biomass. Considering the training data, the mean square error for both the 66 and the 45 parameter model was dominated by errors in timing (phase shift), whereas no general pattern was found in errors when using the validation data. Stronger model reduction should be avoided, as the error term increased and output uncertainty was underestimated.

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.

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.