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 differential impact on soil organic carbon (SOC) of applying no tillage (NT) compared to conventional tillage (CT, i.e. mouldboard ploughing), along with three rates of nitrogen (N) fertilizer application (0,90 and 180 kg ha(-1) y(-1)), was studied under rain-fed Mediterranean conditions in a long-term experiment based on a durum wheat-maize rotation, in which crop residues were left on the soil (NT) or incorporated (CT). Observed SOC content following 8 and 12 years of continuous treatment application was significantly higher in the top 10 cm of the soil under NT than CT, but it was similar in the 10-40 cm layer. NT grain yields for both maize and durum wheat were below those attained under CT (on average 32% and 14% lower respectively) at a given rate of N fertilizer application. Soil, climate and crop data over 5 years were used to calibrate DSSAT model in order to simulate the impact of the different management practices over a 50-year period. Good agreement was obtained between observed and simulated values for crops grain yield, above-ground biomass and observed SOC values. Results from the simulations showed that under NT the weeds growing during the intercrop fallow period made a significant contribution to the observed SOC increase. When the contribution of the weed fallow was considered, NT significantly increased SOC in the top 40 cm of the soil at an average rate of 0.43, 0.31 and 0.03 t ha(-1) per year, respectively for 180,90 and 0 kg N ha(-1) year(-1), within the simulated 50 years. Under CT, a significant SOC increase was simulated under N180 and a significant decrease when no fertilizer was supplied.

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.