Abstract: Livestock farming is the world’s largest land use sector and utilizes around 60% of the global biomass harvest. Over the coming decades, climate change will affect the natural resource base of livestock production, especially the productivity of rangeland and feed crops. Based on a comprehensive impact modeling chain, we assess implications of different climate projections for agricultural production costs and land use change and explore the effectiveness of livestock system transitions as an adaptation strategy. Simulated climate impacts on crop yields and rangeland productivity generate adaptation costs amounting to 3% of total agricultural production costs in 2045 (i.e. 145 billion US$). Shifts in livestock production towards mixed crop-livestock systems represent a resource-and cost-efficient adaptation option, reducing agricultural adaptation costs to 0.3% of total production costs and simultaneously abating deforestation by about 76 million ha globally. The relatively positive climate impacts on grass yields compared with crop yields favor grazing systems inter alia in South Asia and North America. Incomplete transitions in production systems already have a strong adaptive and cost reducing effect: a 50% shift to mixed systems lowers agricultural adaptation costs to 0.8%. General responses of production costs to system transitions are robust across different global climate and crop models as well as regarding assumptions on CO2 fertilization, but simulated values show a large variation. In the face of these uncertainties, public policy support for transforming livestock production systems provides an important lever to improve agricultural resource management and lower adaptation costs, possibly even contributing to emission reduction.

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.