Abstract: Ensembles of process-based crop models are increasingly used to simulate crop growth for scenarios of temperature and/or precipitation changes corresponding to different projections of atmospheric CO2 concentrations. This approach generates large datasets with thousands of simulated crop yield data. Such datasets potentially provide new information but it is difficult to summarize them in a useful way due to their structural complexities. An associated issue is that it is not straightforward to compare crops and to interpolate the results to alternative climate scenarios not initially included in the simulation protocols. Here we demonstrate that statistical models based on random-coefficient regressions are able to emulate ensembles of process-based crop models. An important advantage of the proposed statistical models is that they can interpolate between temperature levels and between CO2 concentration levels, and can thus be used to calculate temperature and [CO2] thresholds leading to yield loss or yield gain, without rerunning the original complex crop models. Our approach is illustrated with three yield datasets simulated by 19 maize models, 26 wheat models, and 13 rice models. Several statistical models are fitted to these datasets, and are then used to analyze the variability of the yield response to [CO2] and temperature. Based on our results, we show that, for wheat, a [CO2] increase is likely to outweigh the negative effect of a temperature increase of +2 degrees C in the considered sites. Compared to wheat, required levels of [CO2] increase are much higher for maize, and intermediate for rice. For all crops, uncertainties in simulating climate change impacts increase more with temperature than with elevated [CO2].

Abstract: We consider predictions of the impact of climate warming on rice development times in Sri Lanka. The major emphasis is on the uncertainty of the predictions, and in particular on the estimation of mean squared error of prediction. Three contributions to mean squared error are considered. The first is parameter uncertainty that results from model calibration. To take proper account of the complex data structure, generalized least squares is used to estimate the parameters and the variance-covariance matrix of the parameter estimators. The second contribution is model structure uncertainty, which we estimate using two different models. An ANOVA analysis is used to separate the contributions of parameter and model uncertainty to mean squared error. The third contribution is model error, which is estimated using hindcasts. Mean squared error of prediction of time from emergence to maturity, for baseline +2 °C, is estimated as 108 days2, with model error contributing 86 days2, followed by model structure uncertainty which contributes 15 days2 and parameter uncertainty which contributes 7 days2. We also show how prediction uncertainty is reduced if prediction concerns development time averaged over years, or the difference in development time between baseline and warmer temperatures.