Abstract: We show the error in water-limited yields simulated by crop models which is associated with spatially aggregated soil and climate input data. Crop simulations at large scales (regional, national, continental) frequently use input data of low resolution. Therefore, climate and soil data are often generated via averaging and sampling by area majority. This may bias simulated yields at large scales, varying largely across models. Thus, we evaluated the error associated with spatially aggregated soil and climate data for 14 crop models. Yields of winter wheat and silage maize were simulated under water-limited production conditions. We calculated this error from crop yields simulated at spatial resolutions from 1 to 100 km for the state of North Rhine-Westphalia, Germany. Most models showed yields biased by <15% when aggregating only soil data. The relative mean absolute error (rMAE) of most models using aggregated soil data was in the range or larger than the inter-annual or inter-model variability in yields. This error increased further when both climate and soil data were aggregated. Distinct error patterns indicate that the rMAE may be estimated from few soil variables. Illustrating the range of these aggregation effects across models, this study is a first step towards an ex-ante assessment of aggregation errors in large-scale simulations.

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: The European Union Effort Sharing Regulation (ESR) will require a 30% reduction in greenhouse gas (GHG) emissions by 2030 compared with 2005 from the sectors not included in the European Emissions Trading Scheme, including agriculture. This will require the estimation of current and future emissions from agriculture, including dairy cattle production systems. Using a farm-scale model as part of a Tier 3 method for farm to national scales provides a more holistic and informative approach than IPCC (2006) Tier 2 but requires independent quality control. Comparing the results of using models to simulate a range of scenarios that explore an appropriate range of biophysical and management situations can support this process by providing a framework for placing model results in context. To assess the variation between models and the process of understanding differences, estimates of GHG emissions from four farm-scale models (DailyWise, FarmAC, HolosNor and SFARMMOD) were calculated for eight dairy farming scenarios within a factorial design consisting of two climates (cool/dry and warm/wet) x two soil types (sandy and clayey) x two feeding systems (grass only and grass/maize). The milk yield per cow, follower cow ratio, manure management system, nitrogen (N) fertilisation and land area were standardised for all scenarios in order to associate the differences in the results with the model structure and function. Potential yield and application of available N in fertiliser and manure were specified separately for grass and maize. Significant differences between models were found in GHG emissions at the farm-scale and for most contributory sources, although there was no difference in the ranking of source magnitudes. The farm-scale GHG emissions, averaged over the four models, was 10.6 t carbon dioxide equivalents (CO(2)e)/ha per year, with a range of 1.9 t CO(2)e/ha per year. Even though key production characteristics were specified in the scenarios, there were still significant differences between models in the annual milk production per ha and the amounts of N fertiliser and concentrate feed imported. This was because the models differed in their description of biophysical responses and feedback mechanisms, and in the extent to which management functions were internalised. We conclude that comparing the results of different farm-scale models when applied to a range of scenarios would build confidence in their use in achieving ESR targets, justifying further investment in the development of a wider range of scenarios and software tools.