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.

Drought risk for agricultural systems in South Africa: Drivers, spatial patterns, and implications for drought risk management

The regular drought episodes in South Africa highlight the need to reduce drought risk by both policy and local community actions. Environmental and socioeconomic factors in South Africa's agricultural system have been affected by drought in the past, creating cascading pressures on the nation's agro-economic and water supply systems. Therefore, understanding the key drivers of all risk components through a comprehensive risk assessment must be undertaken in order to inform proactive drought risk management. This paper presents, for the first time, a national drought risk assessment for irrigated and rainfed systems, that takes into account the complex interaction between different risk components. We use modeling and remote sensing approaches and involve national experts in selecting vulnerability indicators and providing information on human and natural drivers. Our results show that all municipalities have been affected by drought in the last 30 years. The years 1981-1982, 1992, 2016 and 2018 were marked as the driest years during the study period (1981-2018) compared to the reference period (1986-2015). In general, the irrigated systems are remarkably less often affected by drought than rainfed systems; however, most farmers on irrigated land are smallholders for whom drought impacts can be significant. The drought risk of rainfed agricultural systems is exceptionally high in the north, central and west of the country, while for irrigated systems, there are more separate high-risk hotspots across the country. The vulnerability assessment identified potential entry points for disaster risk reduction at the local municipality level, such as increasing environmental awareness, reducing land degradation and increasing total dam and irrigation capacity.

The use of remote sensing to derive maize sowing dates for large-scale crop yield simulations

One of the major sources of uncertainty in large-scale crop modeling is the lack of information capturing the spatiotemporal variability of crop sowing dates. Remote sensing can contribute to reducing such uncertainties by providing essential spatial and temporal information to crop models and improving the accuracy of yield predictions. However, little is known about the impacts of the differences in crop sowing dates estimated by using remote sensing (RS) and other established methods, the uncertainties introduced by the thresholds used in these methods, and the sensitivity of simulated crop yields to these uncertainties in crop sowing dates. In the present study, we performed a systematic sensitivity analysis using various scenarios. The LINTUL-5 crop model implemented in the SIMPLACE modeling platform was applied during the period 2001-2016 to simulate maize yields across four provinces in South Africa using previously defined scenarios of sowing dates. As expected, the selected methodology and the selected threshold considerably influenced the estimated sowing dates (up to 51 days) and resulted in differences in the long-term mean maize yield reaching up to 1.7 t ha^-1 (48% of the mean yield) at the province level. Using RS-derived sowing date estimations resulted in a better representation of the yield variability in space and time since the use of RS information not only relies on precipitation but also captures the impacts of socioeconomic factors on the sowing decision, particularly for smallholder farmers. The model was not able to reproduce the observed yield anomalies in Free State (Pearson correlation coefficient: 0.16 to 0.23) and Mpumalanga (Pearson correlation coefficient: 0.11 to 0.18) in South Africa when using fixed and precipitation rule-based sowing date estimations. Further research with high-resolution climate and soil data and ground-based observations is required to better understand the sources of the uncertainties in RS information and to test whether the results presented herein can be generalized among crop models with different levels of complexity and across distinct field crops.

Effects of soil- and climate data aggregation on simulated potato yield and irrigation water requirement

Input data aggregation affects crop model estimates at the regional level. Previous studies have focused on the impact of aggregating climate data used to compute crop yields. However, little is known about the combined data aggregation effect of climate (DAE_c) and soil (DAE_s) on irrigation water requirement (IWR) in cool-temperate and spatially heterogeneous environments. The aims of this study were to quantify DAE_c and DAE_s of model input data and their combined impacts for simulated irrigated and rainfed yield and IWR. The Agricultural Production Systems sIMulator Next Generation model was applied for the period 1998-2017 across areas suitable for potato (Solanum tuberosum L.) in Tasmania, Australia, using data at 5, 15, 25 and 40 km resolution. Spatial variances of inputs and outputs were evaluated by the relative absolute difference (rAD¯) between the aggregated grids and the 5 km grids. Climate data aggregation resulted in a rAD¯ of 0.7-12.1%, with high values especially for areas with pronounced differences in elevation. The rAD¯ of soil data was higher (5.6-26.3%) than rAD¯ of climate data and was mainly affected by aggregation of organic carbon and maximum plant available water capacity (i.e. the difference between field capacity and wilting point in the effective root zone). For yield estimates, the difference among resolutions (5 km vs. 40 km) was more pronounced for rainfed (rAD¯ = 14.5%) than irrigated conditions (rAD¯ = 3.0%). The rAD¯ of IWR was 15.7% when using input data at 40 km resolution. Therefore, reliable simulations of rainfed yield require a higher spatial resolution than simulation of irrigated yields. This needs to be considered when conducting regional modelling studies across Tasmania. This study also highlights the need to separately quantify the impact of input data aggregation on model outputs to inform about data aggregation errors and identify those variables that explain these errors.

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 CO₂ 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.

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 CO₂ 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 CO₂ . Introducing genotypes adapted to warmer temperatures (and also considering changes in CO₂ 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.

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.

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.

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.