Assessing uncertainties in crop and pasture ensemble model simulations of productivity and N₂O emissions

Simulation models are extensively used to predict agricultural productivity and greenhouse gas (GHG) emissions. However, the uncertainties of (reduced) model ensemble simulations have not been assessed systematically for variables affecting food security and climate change mitigation, within multispecies agricultural contexts. We report an international model comparison and benchmarking exercise, showing the potential of multimodel ensembles to predict productivity and nitrous oxide (N2O) emissions for wheat, maize, rice and temperate grasslands. Using a multistage modelling protocol, from blind simulations (stage 1) to partial (stages 2-4) and full calibration (stage 5), 24 process-based biogeochemical models were assessed individually or as an ensemble against long-term experimental data from four temperate grassland and five arable crop rotation sites spanning four continents. Comparisons were performed by reference to the experimental uncertainties of observed yields and N2O emissions. Results showed that across sites and crop/grassland types, 23 to 40% of the uncalibrated individual models were within two standard deviations (s.d.) of observed yields, while 42 (rice) to 96% (grasslands) of the models were within one s.d. of observed N2O emissions. At stage 1, ensembles formed by the three lowest prediction model errors (RRMSE) predicted both yields and N2O emissions within experimental uncertainties for 44 and 33% of the crop and grassland growth cycles, respectively. Partial model calibration (stages 2 to 4) markedly reduced prediction errors of the full model ensemble E-median for crop grain yields (from 36% at stage 1 down to 4% on average) and grassland productivity (from 44 to 27%) and to a lesser and more variable extent for N2O emissions. Yield-scaled N2O emissions (N2O emissions divided by crop yields) were ranked accurately by 3-model ensembles across crop species and field sites. The potential of using process-based model ensembles to predict jointly productivity and N2O emissions at field scale is discussed