Abstract: This study presents a Bayesian approach for the parameters’ identification of the STICS crop model based on the recently developed Differential Evolution Adaptive Metropolis (DREAM) algorithm. The posterior distributions of nine specific crop parameters of the STICS model were sampled with the aim to improve the growth simulations of a winter wheat (Triticum aestivum L) culture. The results obtained with the DREAM algorithm were initially compared to those obtained with a Nelder-Mead Simplex algorithm embedded within the OptimiSTICS package. Then, three types of likelihood functions implemented within the DREAM algorithm were compared, namely the standard least square, the weighted least square, and a transformed likelihood function that makes explicit use of the coefficient of variation (CV). The results showed that the proposed CV likelihood function allowed taking into account both noise on measurements and heteroscedasticity which are regularly encountered in crop modelling. (C) 2013 Elsevier Ltd. All rights reserved.

Abstract: Reactive nitrogen (Nr) is an indispensable nutrient for agricultural production and human alimentation. Simultaneously, agriculture is the largest contributor to Nr pollution, causing severe damages to human health and ecosystem services. The trade-off between food availability and Nr pollution can be attenuated by several key mitigation options, including Nr efficiency improvements in crop and animal production systems, food waste reduction in households and lower consumption of Nr-intensive animal products. However, their quantitative mitigation potential remains unclear, especially under the added pressure of population growth and changes in food consumption. Here we show by model simulations, that under baseline conditions, Nr pollution in 2050 can be expected to rise to 102-156% of the 2010 value. Only under ambitious mitigation, does pollution possibly decrease to 36-76% of the 2010 value. Air, water and atmospheric Nr pollution go far beyond critical environmental thresholds without mitigation actions. Even under ambitious mitigation, the risk remains that thresholds are exceeded.

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.