High-resolution CMIP6 climate projections for Ethiopia using the gridded statistical downscaling method

High-resolution climate model projections for a range of emission scenarios are needed for designing regional and local adaptation strategies and planning in the context of climate change. To this end, the future climate simulations of global circulation models (GCMs) are the main sources of critical information. However, these simulations are not only coarse in resolution but also associated with biases and high uncertainty. To make the simulations useful for impact modeling at regional and local level, we utilized the bias correction constructed analogues with quantile mapping reordering (BCCAQ) statistical downscaling technique to produce a 10 km spatial resolution climate change projections database based on 16 CMIP6 GCMs under three emission scenarios (SSP2-4.5, SSP3-7.0, and SSP5-8.5). The downscaling strategy was evaluated using a perfect sibling approach and detailed results are presented by taking two contrasting (the worst and best performing models) GCMs as a showcase. The evaluation results demonstrate that the downscaling approach substantially reduced model biases and generated higher resolution daily data compared to the original GCM outputs.

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.

Climate change impact on wheat and maize growth in Ethiopia: A multi-model uncertainty analysis

Ethiopia’s economy is dominated by agriculture which is mainly rain-fed and subsistence. Climate change is expected to have an adverse impact particularly on crop production. Previous studies have shown large discrepancies in the magnitude and sometimes in the direction of the impact on crop production. We assessed the impact of climate change on growth and yield of maize and wheat in Ethiopia using a multi-crop model ensemble. The multi-model ensemble (n = 48) was set up using the agroecosystem modelling framework Expert-N. The framework is modular which facilitates combining different submodels for plant growth and soil processes. The multi-model ensemble was driven by climate change projections representing the mid of the century (2021–2050) from ten contrasting climate models downscaled to finer resolution. The contributions of different sources of uncertainty in crop yield prediction were quantified. The sensitivity of crop yield to elevated CO₂, increased temperature, changes in precipitations and N fertilizer were also assessed. Our results indicate that grain yields were very sensitive to changes in [CO₂], temperature and N fertilizer amounts where the responses were higher for wheat than maize. The response to change in precipitation was weak, which we attribute to the high water holding capacity of the soils due to high organic carbon contents at the study sites. This may provide the sufficient buffering capacity for extended time periods with low amounts of precipitation. Under the changing climate, wheat productivity will be a major challenge with a 36 to 40% reduction in grain yield by 2050 while the impact on maize was modest. A major part of the uncertainty in the projected impact could be attributed to differences in the crop growth models. A considerable fraction of the uncertainty could also be traced back to different soil water dynamics modeling approaches in the model ensemble, which is often ignored. Uncertainties varied among the studied crop species and cultivars as well. The study highlights significant impacts of climate change on wheat yield in Ethiopia whereby differences in crop growth models causes the large part of the uncertainties.

Modelling climate change impacts on maize yields under low nitrogen input conditions in sub-Saharan Africa

Smallholder farmers in sub-Saharan Africa (SSA) currently grow rainfed maize with limited inputs including fertilizer. Climate change may exacerbate current production constraints. Crop models can help quantify the potential impact of climate change on maize yields, but a comprehensive multimodel assessment of simulation accuracy and uncertainty in these low-input systems is currently lacking. We evaluated the impact of varying [CO₂ ], temperature and rainfall conditions on maize yield, for different nitrogen (N) inputs (0, 80, 160 kg N/ha) for five environments in SSA, including cool subhumid Ethiopia, cool semi-arid Rwanda, hot subhumid Ghana and hot semi-arid Mali and Benin using an ensemble of 25 maize models. Models were calibrated with measured grain yield, plant biomass, plant N, leaf area index, harvest index and in-season soil water content from 2-year experiments in each country to assess their ability to simulate observed yield. Simulated responses to climate change factors were explored and compared between models. Calibrated models reproduced measured grain yield variations well with average relative root mean square error of 26%, although uncertainty in model prediction was substantial (CV = 28%). Model ensembles gave greater accuracy than any model taken at random. Nitrogen fertilization controlled the response to variations in [CO₂ ], temperature and rainfall. Without N fertilizer input, maize (a) benefited less from an increase in atmospheric [CO₂ ]; (b) was less affected by higher temperature or decreasing rainfall; and (c) was more affected by increased rainfall because N leaching was more critical. The model intercomparison revealed that simulation of daily soil N supply and N leaching plays a crucial role in simulating climate change impacts for low-input systems. Climate change and N input interactions have strong implications for the design of robust adaptation approaches across SSA, because the impact of climate change in low input systems will be modified if farmers intensify maize production with balanced nutrient management.

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.

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.

How do various maize crop models vary in their responses to climate change factors?

Potential consequences of climate change on crop production can be studied using mechanistic crop simulation models. While a broad variety of maize simulation models exist, it is not known whether different models diverge on grain yield responses to changes in climatic factors, or whether they agree in their general trends related to phenology, growth, and yield. With the goal of analyzing the sensitivity of simulated yields to changes in temperature and atmospheric carbon dioxide concentrations [CO2 ], we present the largest maize crop model intercomparison to date, including 23 different models. These models were evaluated for four locations representing a wide range of maize production conditions in the world: Lusignan (France), Ames (USA), Rio Verde (Brazil) and Morogoro (Tanzania). While individual models differed considerably in absolute yield simulation at the four sites, an ensemble of a minimum number of models was able to simulate absolute yields accurately at the four sites even with low data for calibration, thus suggesting that using an ensemble of models has merit. Temperature increase had strong negative influence on modeled yield response of roughly -0.5 Mg ha(-1) per °C. Doubling [CO2 ] from 360 to 720 μmol mol(-1) increased grain yield by 7.5% on average across models and the sites. That would therefore make temperature the main factor altering maize yields at the end of this century. Furthermore, there was a large uncertainty in the yield response to [CO2 ] among models. Model responses to temperature and [CO2 ] did not differ whether models were simulated with low calibration information or, simulated with high level of calibration information.

A dynamical model of environmental effects on allocation to carbon-based secondary compounds in juvenile trees

BACKGROUND AND AIMS

Patterns and variations in concentration of carbon-based secondary compounds in plant tissues have been explained by means of different complementary and, in some cases, contradictory plant defence hypotheses for more than 20 years. These hypotheses are conceptual models which consider environmental impacts on plant internal demands. In the present study, a mathematical model is presented, which converts and integrates the concepts of the 'Growth-Differentiation Balance' hypothesis and the 'Protein Competition' model into a dynamic plant growth model, that was tested with concentration data of polyphenols in leaves of juvenile apple, beech and spruce trees. The modelling approach is part of the plant growth model PLATHO that considers simultaneously different environmental impacts on the most important physiological processes of plants.

METHODS

The modelling approach for plant internal resource allocation is based on a priority scheme assuming that growth processes have priority over allocation to secondary compounds and that growth-related metabolism is more strongly affected by nitrogen deficiency than defence-related secondary metabolism.

KEY RESULTS

It is shown that the model can reproduce the effect of nitrogen fertilization on allocation patterns in apple trees and the effects of elevated CO(2) and competition in juvenile beech and spruce trees. The analysis of model behaviour reveals that large fluctuations in plant internal availability of carbon and nitrogen are possible within a single vegetation period. Furthermore, the model displays a non-linear allocation behaviour to carbon-based secondary compounds.

CONCLUSIONS

The simulation results corroborate the underlying assumptions of the presented modelling approach for resource partitioning between growth-related primary metabolism and defence-related secondary metabolism. Thus, the dynamical modelling approach, which considers variable source and sink strengths of plant internal resources within different phenological growth stages, presents a successful translation of existing concepts into a dynamic mathematical model.

Analysis of competition effects in mono- and mixed cultures of juvenile beech and spruce by means of the plant growth simulation model PLATHO

Inter- and intra-specific competition between plants for external resources is a critical process for plant growth in natural and managed ecosystems. We present a new approach to simulate competition for the resources light, water, and nitrogen between individual plants within a canopy. This approach was incorporated in a process-oriented plant growth simulation model. The concept of modelling competition is based on competition coefficients calculated from the overlap of occupied crown and soil volumes of each plant individual with the occupied volumes of its four nearest neighbours. The model was parameterised with data from a two-year phytotron experiment with juvenile beech and spruce trees growing in mono- and mixed cultures. For testing the model, an independent data set from this experiment and data from a second phytotron experiment with mixed cultures were used. The model was applied to analyse the consequences of start conditions and plant density on plant-plant competition. In both experiments, spruce dominated beech in mixed cultures. Based on model simulations, we postulate a large influence of start conditions and stand density on the outcome of the competition between the species. When both species have similar heights at the time of canopy closure, the model suggests a greater morphological plasticity of beech compared with spruce to be the crucial mechanism for competitiveness in mixed canopies. Similar to the experiment, in the model greater plasticity was a disadvantage for beech leading to it being outcompeted by the more persistent spruce.

Quantifying the effect of soil moisture on the aerobic microbial mineralization of selected pesticides in different soils

A standardized quantitative approach was developed to reliably elucidate the effect of increasing soil moisture on pesticide mineralization. The mineralization of three aerobically degradable and chemically different 14C-labeled pesticides (isoproturon, benazolin-ethyl, and glyphosate) was studied under controlled conditions in the laboratory at an identical soil density of 1.3 g cm(-3). The agricultural soils used are characterized by (i) large variations in soil texture (sand content 4-88%) and organic matter content (0.97-2.70% org. C), (ii) fairly diverse soil-water retention curves, and (iii) differing pH values. We quantified the effect of soil moisture on mineralization of pesticides and found that (i) at soil water potential < or = -20 MPa minimal pesticide mineralization occurred; (ii) a linear correlation (P < 0.0001) exists between increasing soil moisture (within a soil water potential range of -20 and -0.015 MPa), and increased relative pesticide mineralization; (iii) optimum pesticide mineralization was obtained at a soil water potential of -0.015 MPa, and (iv) when soil moisture approximated water holding capacity, pesticide mineralization was considerably reduced. As both selected pesticides and soils varied to a large degree, we propose that the correlation observed in this study may be also valid in the case of aerobic degradation of other native and artificial organic compounds in soils.

Contribution of current photosynthates to root respiration of non-nodulated Medicago sativa: effects of light and nitrogen supply

The effects of light (PFD) and nitrogen (N) supply on root respiration of new C (currently assimilated carbon, R(new)) and old C ( R(old)) were analysed in non-nodulated Medicago sativa. Plants were pre-treated with high/low PFD and high/low N supply with a regular 16/8 h light/dark cycle. Five to eight weeks after planting current photosynthates were labelled with (13)C and their contribution to root respiration was continuously measured during a 24 h day/night cycle. PFD conditions during labelling were either those of the pre-treatments (control, 25 or 6 mol m(-2) d(-1)) or, for high PFD plants, 6 mol m(-2) d(-1) by shortening the photoperiod or reducing irradiance. The fraction of new C in the respiratory CO2 increased during the light period, but remained constant in the dark period. In control plants, R(new) contributed 40 % to the daily root respiration in high PFD/high N conditions. Continuously low PFD increased (50 %) and low N decreased (26 %) the contribution of R(new). Exposing plants from high PFD pre-treatments to a short photoperiod or to low PFD stimulated R(old), indicating mobilisation of reserve C. This stimulation was more pronounced in plants with high N supply than in those with low N supply. Comparison with other legumes suggested that R(new) in root respiration was mainly defined by the ratio between the assimilatory capacity of the shoots and the maintenance costs of roots with a short-term capacity of buffering respiratory demand by mobilisation of reserves in situations of fluctuating PFD.

Uptake of terbuthylazine and its medium polar metabolites into maize plants

The uptake of terbuthylazine and its medium polar metabolites into maize plants under outdoor conditions is investigated. For this purpose, a dynamical fate model consisting of soil, plant and air is developed. The model calculations are compared with experimental results of outdoor lysimeter tests, carried out with(14)C-labelled herbicide applied to sandy agricultural soil at a single application rate of 890 g/ha. Approximately 0.3 % of the applied activity remains in all the plants after the vegetation period. The model predicts that about three times that amount is volatilized from the plants into the air. Activity uptaken from soil and volatilized from plant surface into air is predominately associated with metabolites. During the whole vegetation period the fraction of unchanged terbuthylazine in the plants is very small (less than 1 % of the extractable activity).