Simulation of winter wheat response to variable sowing dates and densities in a high-yielding environment.

Crop multi-model ensembles (MME) have proven to be effective in increasing the accuracy of simulations in modelling experiments. However, the ability of MME to capture crop responses to changes in sowing dates and densities has not yet been investigated. These management interventions are some of the main levers for adapting cropping systems to climate change. Here, we explore the performance of a MME of 29 wheat crop models to predict the effect of changing sowing dates and rates on yield and yield components, on two sites located in a high-yielding environment in New Zealand. The experiment was conducted for 6 years and provided 50 combinations of sowing date, sowing density and growing season. We show that the MME simulates seasonal growth of wheat well under standard sowing conditions, but fails under early sowing and high sowing rates. The comparison between observed and simulated in-season fraction of intercepted photosynthetically active radiation (FIPAR) for early sown wheat shows that the MME does not capture the decrease of crop above ground biomass during winter months due to senescence. Models need to better account for tiller competition for light, nutrients, and water during vegetative growth, and early tiller senescence and tiller mortality, which are exacerbated by early sowing, high sowing densities, and warmer winter temperatures.

Crop loss identification at field parcel scale using satellite remote sensing and machine learning.

Identifying crop loss at field parcel scale using satellite images is challenging: first, crop loss is caused by many factors during the growing season; second, reliable reference data about crop loss are lacking; third, there are many ways to define crop loss. This study investigates the feasibility of using satellite images to train machine learning (ML) models to classify agricultural field parcels into those with and without crop loss. The reference data for this study was provided by Finnish Food Authority (FFA) containing crop loss information of approximately 1.4 million field parcels in Finland covering about 3.5 million ha from 2000 to 2015. This reference data was combined with Normalised Difference Vegetation Index (NDVI) derived from Landsat 7 images, in which more than 80% of the possible data are missing. Despite the hard problem with extremely noisy data, among the four ML models we tested, random forest (with mean imputation and missing value indicators) achieved the average AUC (area under the ROC curve) of 0.688±0.059 over all 16 years with the range [0.602, 0.795] in identifying new crop-loss fields based on reference fields of the same year. To our knowledge, this is one of the first large scale benchmark study of using machine learning for crop loss classification at field parcel scale. The classification setting and trained models have numerous potential applications, for example, allowing government agencies or insurance companies to verify crop-loss claims by farmers and realise efficient agricultural monitoring.

Comparative Analysis of Phenology Algorithms of the Spring Barley Model in APSIM 7.9 and APSIM Next Generation: A Case Study for High Latitudes.

Phenology algorithms in crop growth models have inevitable systematic errors and uncertainties. In this study, the phenology simulation algorithms in APSIM classical (APSIM 7.9) and APSIM next generation (APSIM-NG) were compared for spring barley models at high latitudes. Phenological data of twelve spring barley varieties were used for the 2014-2018 cropping seasons from northern Sweden and Finland. A factorial-based calibration approach provided within APSIM-NG was performed to calibrate both models. The models have different mechanisms to simulate days to anthesis. The calibration was performed separately for days to anthesis and physiological maturity, and evaluations for the calibrations were done with independent datasets. The calibration performance for both growth stages of APSIM-NG was better compared to APSIM 7.9. However, in the evaluation, APSIM-NG showed an inclination to overestimate days to physiological maturity. The differences between the models are possibly due to slower thermal time accumulation mechanism, with higher cardinal temperatures in APSIM-NG. For a robust phenology prediction at high latitudes with APSIM-NG, more research on the conception of thermal time computation and implementation is suggested.

Climate change impact and adaptation for wheat protein.

Wheat grain protein concentration is an important determinant of wheat quality for human nutrition that is often overlooked in efforts to improve crop production. We tested and applied a 32-multi-model ensemble to simulate global wheat yield and quality in a changing climate. Potential benefits of elevated atmospheric CO2 concentration by 2050 on global wheat grain and protein yield are likely to be negated by impacts from rising temperature and changes in rainfall, but with considerable disparities between regions. Grain and protein yields are expected to be lower and more variable in most low-rainfall regions, with nitrogen availability limiting growth stimulus from elevated CO2 . Introducing genotypes adapted to warmer temperatures (and also considering changes in CO2 and rainfall) could boost global wheat yield by 7% and protein yield by 2%, but grain protein concentration would be reduced by -1.1 percentage points, representing a relative change of -8.6%. Climate change adaptations that benefit grain yield are not always positive for grain quality, putting additional pressure on global wheat production.

Global wheat production with 1.5 and 2.0°C above pre-industrial warming.

Efforts to limit global warming to below 2°C in relation to the pre-industrial level are under way, in accordance with the 2015 Paris Agreement. However, most impact research on agriculture to date has focused on impacts of warming >2°C on mean crop yields, and many previous studies did not focus sufficiently on extreme events and yield interannual variability. Here, with the latest climate scenarios from the Half a degree Additional warming, Prognosis and Projected Impacts (HAPPI) project, we evaluated the impacts of the 2015 Paris Agreement range of global warming (1.5 and 2.0°C warming above the pre-industrial period) on global wheat production and local yield variability. A multi-crop and multi-climate model ensemble over a global network of sites developed by the Agricultural Model Intercomparison and Improvement Project (AgMIP) for Wheat was used to represent major rainfed and irrigated wheat cropping systems. Results show that projected global wheat production will change by -2.3% to 7.0% under the 1.5°C scenario and -2.4% to 10.5% under the 2.0°C scenario, compared to a baseline of 1980-2010, when considering changes in local temperature, rainfall, and global atmospheric CO2 concentration, but no changes in management or wheat cultivars. The projected impact on wheat production varies spatially; a larger increase is projected for temperate high rainfall regions than for moderate hot low rainfall and irrigated regions. Grain yields in warmer regions are more likely to be reduced than in cooler regions. Despite mostly positive impacts on global average grain yields, the frequency of extremely low yields (bottom 5 percentile of baseline distribution) and yield inter-annual variability will increase under both warming scenarios for some of the hot growing locations, including locations from the second largest global wheat producer-India, which supplies more than 14% of global wheat. The projected global impact of warming <2°C on wheat production is therefore not evenly distributed and will affect regional food security across the globe as well as food prices and trade.

Multimodel ensembles improve predictions of crop-environment-management interactions.

A recent innovation in assessment of climate change impact on agricultural production has been to use crop multimodel ensembles (MMEs). These studies usually find large variability between individual models but that the ensemble mean (e-mean) and median (e-median) often seem to predict quite well. However, few studies have specifically been concerned with the predictive quality of those ensemble predictors. We ask what is the predictive quality of e-mean and e-median, and how does that depend on the ensemble characteristics. Our empirical results are based on five MME studies applied to wheat, using different data sets but the same 25 crop models. We show that the ensemble predictors have quite high skill and are better than most and sometimes all individual models for most groups of environments and most response variables. Mean squared error of e-mean decreases monotonically with the size of the ensemble if models are added at random, but has a minimum at usually 2-6 models if best-fit models are added first. Our theoretical results describe the ensemble using four parameters: average bias, model effect variance, environment effect variance, and interaction variance. We show analytically that mean squared error of prediction (MSEP) of e-mean will always be smaller than MSEP averaged over models and will be less than MSEP of the best model if squared bias is less than the interaction variance. If models are added to the ensemble at random, MSEP of e-mean will decrease as the inverse of ensemble size, with a minimum equal to squared bias plus interaction variance. This minimum value is not necessarily small, and so it is important to evaluate the predictive quality of e-mean for each target population of environments. These results provide new information on the advantages of ensemble predictors, but also show their limitations.

Contribution of crop model structure, parameters and climate projections to uncertainty in climate change impact assessments.

Climate change impact assessments are plagued with uncertainties from many sources, such as climate projections or the inadequacies in structure and parameters of the impact model. Previous studies tried to account for the uncertainty from one or two of these. Here, we developed a triple-ensemble probabilistic assessment using seven crop models, multiple sets of model parameters and eight contrasting climate projections together to comprehensively account for uncertainties from these three important sources. We demonstrated the approach in assessing climate change impact on barley growth and yield at Jokioinen, Finland in the Boreal climatic zone and Lleida, Spain in the Mediterranean climatic zone, for the 2050s. We further quantified and compared the contribution of crop model structure, crop model parameters and climate projections to the total variance of ensemble output using Analysis of Variance (ANOVA). Based on the triple-ensemble probabilistic assessment, the median of simulated yield change was -4% and +16%, and the probability of decreasing yield was 63% and 31% in the 2050s, at Jokioinen and Lleida, respectively, relative to 1981-2010. The contribution of crop model structure to the total variance of ensemble output was larger than that from downscaled climate projections and model parameters. The relative contribution of crop model parameters and downscaled climate projections to the total variance of ensemble output varied greatly among the seven crop models and between the two sites. The contribution of downscaled climate projections was on average larger than that of crop model parameters. This information on the uncertainty from different sources can be quite useful for model users to decide where to put the most effort when preparing or choosing models or parameters for impact analyses. We concluded that the triple-ensemble probabilistic approach that accounts for the uncertainties from multiple important sources provide more comprehensive information for quantifying uncertainties in climate change impact assessments as compared to the conventional approaches that are deterministic or only account for the uncertainties from one or two of the uncertainty sources.

Author Correction: The uncertainty of crop yield projections is reduced by improved temperature response functions.

Nature Plants 3, 17102 (2017); published online 17 July 2017; corrected online 27 September 2017.

Erratum: The uncertainty of crop yield projections is reduced by improved temperature response functions.

This corrects the article DOI: 10.1038/nplants.2017.102.

Erratum: The uncertainty of crop yield projections is reduced by improved temperature response functions.

This corrects the article DOI: 10.1038/nplants.2017.102.

The uncertainty of crop yield projections is reduced by improved temperature response functions.

Increasing the accuracy of crop productivity estimates is a key element in planning adaptation strategies to ensure global food security under climate change. Process-based crop models are effective means to project climate impact on crop yield, but have large uncertainty in yield simulations. Here, we show that variations in the mathematical functions currently used to simulate temperature responses of physiological processes in 29 wheat models account for >50% of uncertainty in simulated grain yields for mean growing season temperatures from 14 °C to 33 °C. We derived a set of new temperature response functions that when substituted in four wheat models reduced the error in grain yield simulations across seven global sites with different temperature regimes by 19% to 50% (42% average). We anticipate the improved temperature responses to be a key step to improve modelling of crops under rising temperature and climate change, leading to higher skill of crop yield projections.

Multimodel ensembles of wheat growth: many models are better than one.

Crop models of crop growth are increasingly used to quantify the impact of global changes due to climate or crop management. Therefore, accuracy of simulation results is a major concern. Studies with ensembles of crop models can give valuable information about model accuracy and uncertainty, but such studies are difficult to organize and have only recently begun. We report on the largest ensemble study to date, of 27 wheat models tested in four contrasting locations for their accuracy in simulating multiple crop growth and yield variables. The relative error averaged over models was 24-38% for the different end-of-season variables including grain yield (GY) and grain protein concentration (GPC). There was little relation between error of a model for GY or GPC and error for in-season variables. Thus, most models did not arrive at accurate simulations of GY and GPC by accurately simulating preceding growth dynamics. Ensemble simulations, taking either the mean (e-mean) or median (e-median) of simulated values, gave better estimates than any individual model when all variables were considered. Compared to individual models, e-median ranked first in simulating measured GY and third in GPC. The error of e-mean and e-median declined with an increasing number of ensemble members, with little decrease beyond 10 models. We conclude that multimodel ensembles can be used to create new estimators with improved accuracy and consistency in simulating growth dynamics. We argue that these results are applicable to other crop species, and hypothesize that they apply more generally to ecological system models.

Greenhouse impact due to the use of combustible fuels: life cycle viewpoint and relative radiative forcing commitment.

Extensive information on the greenhouse impacts of various human actions is important in developing effective climate change mitigation strategies. The greenhouse impacts of combustible fuels consist not only of combustion emissions but also of emissions from the fuel production chain and possible effects on the ecosystem carbon storages. It is important to be able to assess the combined, total effect of these different emissions and to express the results in a comprehensive way. In this study, a new concept called relative radiative forcing commitment (RRFC) is presented and applied to depict the greenhouse impact of some combustible fuels currently used in Finland. RRFC is a ratio that accounts for the energy absorbed in the Earth system due to changes in greenhouse gas concentrations (production and combustion of fuel) compared to the energy released in the combustion of fuel. RRFC can also be expressed as a function of time in order to give a dynamic cumulative picture on the caused effect. Varying time horizons can be studied separately, as is the case when studying the effects of different climate policies on varying time scales. The RRFC for coal for 100 years is about 170, which means that in 100 years 170 times more energy is absorbed in the atmosphere due to the emissions of coal combustion activity than is released in combustion itself. RRFC values of the other studied fuel production chains varied from about 30 (forest residues fuel) to 190 (peat fuel) for the 100-year study period. The length of the studied time horizon had an impact on the RRFC values and, to some extent, on the relative positions of various fuels.