Abstract: Croplands and grasslands are agricultural systems that contribute to land–atmosphere exchanges of carbon (C). We evaluated and compared gross primary production (GPP), ecosystem respiration (RECO), net ecosystem exchange (NEE) of CO2, and two derived outputs – C use efficiency (CUE=-NEE/GPP) and C emission intensity (IntC= -NEE/Offtake [grazed or harvested biomass]). The outputs came from 23 models (11 crop-specific, eight grassland-specific, and four models covering both systems) at three cropping sites over several rotations with spring and winter cereals, soybean and rapeseed in Canada, France and India, and two temperate permanent grasslands in France and the United Kingdom. The models were run independently over multi-year simulation periods in five stages (S), either blind with no calibration and initialization data (S1), using historical management and climate for initialization (S2), calibrated against plant data (S3), plant and soil data together (S4), or with the addition of C and N fluxes (S5). Here, we provide a framework to address methodological uncertainties and contextualize results. Most of the models overestimated or underestimated the C fluxes observed during the growing seasons (or the whole years for grasslands), with substantial differences between models. For each simulated variable, changes in the multi-model median (MMM) from S1 to S5 was used as a descriptor of the ensemble performance. Overall, the greatest improvements (MMM approaching the mean of observations) were achieved at S3 or higher calibration stages. For instance, grassland GPP MMM was equal to 1632 g C m−2 yr-1 (S5) while the observed mean was equal to 1763 m-2 yr-1 (average for two sites). Nash-Sutcliffe modelling efficiency coefficients indicated that MMM outperformed individual models in 92.3 % of cases. Our study suggests a cautious use of large-scale, multi-model ensembles to estimate C fluxes in agricultural sites if some site-specific plant and soil observations are available for model calibration. The further development of crop/grassland ensemble modelling will hinge upon the interpretation of results in light of the way models represent the processes underlying C fluxes in complex agricultural systems (grassland and crop rotations including fallow periods).

Abstract: • We simulate biomass, soil water content (SWC) and temperature (ST) in grasslands. • We compare nine models to the multi-model median (MMM) at nine sites. • With model calibration, we obtain satisfactory estimates of ST, less of SWC and biomass. • We observe discrepancies across models in the simulation of grassland processes. • We improve performance with multi-model approach. This study presents results from a major grassland model intercomparison exercise, and highlights the main challenges faced in the implementation of a multi-model ensemble prediction system in grasslands. Nine, independently developed simulation models linking climate, soil, vegetation and management to grassland biogeochemical cycles and production were compared in a simulation of soil water content (SWC) and soil temperature (ST) in the topsoil, and of biomass production. The results were assessed against SWC and ST data from five observational grassland sites representing a range of conditions – Grillenburg in Germany, Laqueuille in France with both extensive and intensive management, Monte Bondone in Italy and Oensingen in Switzerland – and against yield measurements from the same sites and other experimental grassland sites in Europe and Israel. We present a comparison of model estimates from individual models to the multi-model ensemble (represented by multi-model median: MMM). With calibration (seven out of nine models), the performances were acceptable for weekly-aggregated ST (R² > 0.7 with individual models and >0.8–0.9 with MMM), but less satisfactory with SWC (R² < 0.6 with individual models and < ∼ 0.5 with MMM) and biomass (R² < ∼0.3 with both individual models and MMM). With individual models, maximum biases of about −5 °C for ST, −0.3 m3 m−3 for SWC and 360 g DM m−2 for yield, as well as negative modelling efficiencies and some high relative root mean square errors indicate low model performance, especially for biomass. We also found substantial discrepancies across different models, indicating considerable uncertainties regarding the simulation of grassland processes. The multi-model approach allowed for improved performance, but further progress is strongly needed in the way models represent processes in managed grassland systems.

Abstract: Predicting regional and global carbon (C) and water dynamics on grasslands has become of major interest, as grasslands are one of the most widespread vegetation types worldwide, providing a number of ecosystem services (such as forage production and C storage). The present study is a contribution to a regional-scale analysis of the C and water cycles on managed grasslands. The mechanistic biogeochemical model PaSim (Pasture Simulation model) was evaluated at 12 grassland sites in Europe. A new parameterization was obtained on a common set of eco-physiological parameters, which represented an improvement of previous parameterization schemes (essentially obtained via calibration at specific sites). We found that C and water fluxes estimated with the parameter set are in good agreement with observations. The model with the new parameters estimated that European grassland are a sink of C with 213 g C m(-2) yr(-1), which is close to the observed net ecosystem exchange (NEE) flux of the studied sites (185 g C m(-2) yr(-1) on average). The estimated yearly average gross primary productivity (GPP) and ecosystem respiration (RECO) for all of the study sites are 1220 and 1006 g C m(-2) yr(-1), respectively, in agreement with observed average GPP (1230 g C m(-2) yr(-1)) and RECO (1046 g C m(-2) yr(-1)). For both variables aggregated on a weekly basis, the root mean square error (RMSE) was similar to 5-16 g C week(-1) across the study sites, while the goodness of fit (R-2) was similar to 0.4-0.9. For evapotranspiration (ET), the average value of simulated ET (415 mmyr(-1)) for all sites and years is close to the average value of the observed ET (451 mm yr(-1)) by flux towers (on a weekly basis, RMSE similar to 2-8 mm week(-1); R-2 = 0.3-0.9). However, further model development is needed to better represent soil water dynamics under dry conditions and soil temperature in winter. A quantification of the uncertainties introduced by spatially generalized parameter values in C and water exchange estimates is also necessary. In addition, some uncertainties in the input management data call for the need to improve the quality of the observational system.