Abstract: Justifiable usage of large-scale crop model simulations requires transparent, comprehensive and spatially extensive evaluations of their performance and associated accuracy. Simulated crop yields of a Pan-European implementation of the Environmental Policy Integrated Climate (EPIC) crop model were satisfactorily evaluated with reported regional yield data from EUROSTAT for four major crops, including winter wheat, rainfed and irrigated maize, spring barley and winter rye. European-wide land use, elevation, soil and daily meteorological gridded data were integrated in GIS and coupled with EPIC. Default EPIC crop and biophysical process parameter values were used with some minor adjustments according to suggestions from scientific literature. The model performance was improved by spatial calculations of crop sowing densities, potential heat units, operation schedules, and nutrient application rates. EPIC performed reasonable in the simulation of regional crop yields, with long-term averages predicted better than inter-annual variability: linear regression R-2 ranged from 0.58 (maize) to 0.91 (spring barley) and relative estimation errors were between +/- 30% for most of the European regions. The modelled and reported crop yields demonstrated similar responses to driving meteorological variables. However, EPIC performed better in dry compared to wet years. A yield sensitivity analysis of crop nutrient and irrigation management factors and cultivar specific characteristics for contrasting regions in Europe revealed a range in model response and attainable yields. We also show that modelled crop yield is strongly dependent on the chosen PET method. The simulated crop yield variability was lower compared to reported crop yields. This assessment should contribute to the availability of harmonised and transparently evaluated agricultural modelling tools in the EU as well as the establishment of modelling benchmarks as a requirement for sound and ongoing policy evaluations in the agricultural and environmental domains. (C) 2013 The Authors. Published by Elsevier Ltd. All rights reserved.

Abstract: Eddy covariance data from four European grassland sites are used to probabilistically invert the CARAIB (CARbon Assimilation In the Biosphere) dynamic vegetation model (DVM) with 10 unknown parameters, using the DREAM((ZS)) (DiffeRential Evolution Adaptive Metropolis) Markov chain Monte Carlo (MCMC) sampler. We focus on comparing model inversions, considering both homoscedastic and heteroscedastic eddy covariance residual errors, with variances either fixed a priori or jointly inferred together with the model parameters. Agreements between measured and simulated data during calibration are comparable with previous studies, with root mean square errors (RMSEs) of simulated daily gross primary productivity (GPP), ecosystem respiration (RECO) and evapotranspiration (ET) ranging from 1.73 to 2.19, 1.04 to 1.56 g C m(-2) day(-1) and 0.50 to 1.28 mm day(-1), respectively. For the calibration period, using a homoscedastic eddy covariance residual error model resulted in a better agreement between measured and modelled data than using a heteroscedastic residual error model. However, a model validation experiment showed that CARAIB models calibrated considering heteroscedastic residual errors perform better. Posterior parameter distributions derived from using a heteroscedastic model of the residuals thus appear to be more robust. This is the case even though the classical linear heteroscedastic error model assumed herein did not fully remove heteroscedasticity of the GPP residuals. Despite the fact that the calibrated model is generally capable of fitting the data within measurement errors, systematic bias in the model simulations are observed. These are likely due to model inadequacies such as shortcomings in the photosynthesis modelling. Besides the residual error treatment, differences between model parameter posterior distributions among the four grassland sites are also investigated. It is shown that the marginal distributions of the specific leaf area and characteristic mortality time parameters can be explained by site-specific ecophysiological characteristics.

Abstract: We compare 27 wheat models’ yield responses to interannual climate variability, analyzed at locations in Argentina, Australia, India, and The Netherlands as part of the Agricultural Model Intercomparison and Improvement Project (AgMIP) Wheat Pilot. Each model simulated 1981-2010 grain yield, and we evaluate results against the interannual variability of growing season temperature, precipitation, and solar radiation. The amount of information used for calibration has only a minor effect on most models’ climate response, and even small multi-model ensembles prove beneficial. Wheat model clusters reveal common characteristics of yield response to climate; however models rarely share the same cluster at all four sites indicating substantial independence. Only a weak relationship (R-2 <= 0.24) was found between the models’ sensitivities to interannual temperature variability and their response to long-term warming, suggesting that additional processes differentiate climate change impacts from observed climate variability analogs and motivating continuing analysis and model development efforts. Published by Elsevier Ltd.