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: • This study combines crop, livestock and economic models.
• Models interaction is through use of relevant input and output variables.
• Future climate change will result in increased grass and wheat dry matter yields.
• Changes in grass, wheat and milk yields in future reduce farm emissions intensity.
• Changes in future dry matter yields and emissions lead to increased profitability.
There is a scientific consensus that the future climate change will affect grass and crop dry matter (DM) yields. Such yield changes may entail alterations to farm management practices to fulfill the feed requirements and reduce the farm greenhouse gas (GHG) emissions from dairy farms. While a large number of studies have focused on the impacts of projected climate change on a single farm output (e.g. GHG emissions or economic performance), several attempts have been made to combine bio-economic systems models with GHG accounting frameworks. In this study, we aimed to determine the physical impacts of future climate scenarios on grass and wheat DM yields, and demonstrate the effects such changes in future feed supply may have on farm GHG emissions and decision-making processes. For this purpose, we combined four models: BASGRA and CSM-CERES-Wheat models for simulating forage grass DM and wheat DM grain yields respectively; HolosNor for estimating the farm GHG emissions; and JORDMOD for calculating the impacts of changes in the climate and management on land use and farm economics. Four locations, with varying climate and soil conditions were included in the study: south-east Norway, south-west Norway, central Norway and northern Norway. Simulations were carried out for baseline (1961–1990) and future (2046–2065) climate conditions (projections based on two global climate models and the Special Report on Emissions Scenarios (SRES) A1B GHG emission scenario), and for production conditions with and without a milk quota. The GHG emissions intensities (kilogram carbon dioxide equivalent: kgCO2e emissions per kg fat and protein corrected milk: FPCM) varied between 0.8 kg and 1.23 kg CO2e (kg FPCM)− 1, with the lowest and highest emissions found in central Norway and south-east Norway, respectively. Emission intensities were generally lower under future compared to baseline conditions due mainly to higher future milk yields and to some extent to higher crop yields. The median seasonal above-ground timothy grass yield varied between 11,000 kg and 16,000 kg DM ha− 1 and was higher in all projected future climate conditions than in the baseline. The spring wheat grain DM yields simulated for the same weather conditions within each climate projection varied between 2200 kg and 6800 kg DM ha− 1. Similarly, the farm profitability as expressed by total national land rents varied between 1900 million Norwegian krone (NOK) for median yields under baseline climate conditions up to 3900 million NOK for median yield under future projected climate conditions.