Abstract: Grassland-based ruminant production systems are integral to sustainable food production in Europe, converting plant materials indigestible to humans into nutritious food, while providing a range of environmental and cultural benefits. Climate change poses significant challenges for such systems, their productivity and the wider benefits they supply. In this context, grassland models have an important role in predicting and understanding the impacts of climate change on grassland systems, and assessing the efficacy of potential adaptation and mitigation strategies. In order to identify the key challenges for European grassland modelling under climate change, modellers and researchers from across Europe were consulted via workshop and questionnaire. Participants identified fifteen challenges and considered the current state of modelling and priorities for future research in relation to each. A review of literature was undertaken to corroborate and enrich the information provided during the horizon scanning activities. Challenges were in four categories relating to: 1) the direct and indirect effects of climate change on the sward 2) climate change effects on grassland systems outputs 3) mediation of climate change impacts by site, system and management and 4) cross-cutting methodological issues. While research priorities differed between challenges, an underlying theme was the need for accessible, shared inventories of models, approaches and data, as a resource for stakeholders and to stimulate new research. Developing grassland models to effectively support efforts to tackle climate change impacts, while increasing productivity and enhancing ecosystem services, will require engagement with stakeholders and policy-makers, as well as modellers and experimental researchers across many disciplines. The challenges and priorities identified are intended to be a resource 1) for grassland modellers and experimental researchers, to stimulate the development of new research directions and collaborative opportunities, and 2) for policy-makers involved in shaping the research agenda for European grassland modelling under climate change.

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.