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: • ARMOSA model simulated a maize ideotype with drought adaptation under climate change. • The ideotype needs less water for higher yield compared to current hybrids. • Higher production involves more crop residues that enhance soil C sequestration. • Soil organic C may generally decrease and N leaching will increase in sandy soil. The expected climate change will affect the maize yields in view of air temperature increase and scarce water availability. The application of biophysical models offers the chance to design a drought-resistant ideotype and to assist plant breeders and agronomists in the assessment of its suitability in future scenarios. The aim of the present work was to perform a model-based estimation of the yields of two hybrids, current vs ideotype, under future climate scenarios (2030–2060 and 2070–2100) in Lombardy (northern Italy), testing two options of irrigation (small amount at fixed dates vs optimal water supply), nitrogen (N) fertilization (300 vs 400 kg N ha− 1), and crop cycle durations (current vs extended). For the designing of the ideotype we set several parameters of the ARMOSA process-based crop model: the root elongation rate and maximum depth, stomatal resistance, four stage-specific crop coefficients for the actual transpiration estimation, and drought tolerance factor. The work findings indicated that the current hybrid ensures good production only with high irrigation amount (245–565 mm y− 1). With respect to the current hybrid, the ideotype will require less irrigation water (− 13%, p < 0.01) and it resulted in significantly higher yield under water stress condition (+ 15%, p < 0.01) and optimal water supply (+ 2%, p < 0.05). The elongated cycle has a positive effect on yield under any combination of options. Moreover, higher yields projected for the ideotype implicate more crop residues to be incorporated into the soil, which are positively correlated with the SOC sequestration and negatively with N leaching. The crop N uptake is expected to be adequate in view of higher rate of soil mineralization; the N fertilization rate of 400 kg N ha− 1 will involve significant increasing of grain yield, and it is expected to involve a higher rate of SOC sequestration.

Abstract: Working with ensembles of crop models is a recent but important development in crop modeling which promises to lead to better uncertainty estimates for model projections and predictions, better predictions using the ensemble mean or median, and closer collaboration within the modeling community. There are numerous open questions about the best way to create and analyze such ensembles. Much can be learned from the field of climate modeling, given its much longer experience with ensembles. We draw on that experience to identify questions and make propositions that should help make ensemble modeling with crop models more rigorous and informative. The propositions include defining criteria for acceptance of models in a crop MME, exploring criteria for evaluating the degree of relatedness of models in a MME, studying the effect of number of models in the ensemble, development of a statistical model of model sampling, creation of a repository for MME results, studies of possible differential weighting of models in an ensemble, creation of single model ensembles based on sampling from the uncertainty distribution of parameter values or inputs specifically oriented toward uncertainty estimation, the creation of super ensembles that sample more than one source of uncertainty, the analysis of super ensemble results to obtain information on total uncertainty and the separate contributions of different sources of uncertainty and finally further investigation of the use of the multi-model mean or median as a predictor.