Exploring uncertainties in global crop yield projections in a large ensemble of crop models and CMIP5 and CMIP6 climate scenarios

Concerns over climate change are motivated in large part because of their impact on human society. Assessing the effect of that uncertainty on specific potential impacts is demanding, since it requires a systematic survey over both climate and impacts models. We provide a comprehensive evaluation of uncertainty in projected crop yields for maize, spring and winter wheat, rice, and soybean, using a suite of nine crop models and up to 45 CMIP5 and 34 CMIP6 climate projections for three different forcing scenarios. To make this task computationally tractable, we use a new set of statistical crop model emulators. We find that climate and crop models contribute about equally to overall uncertainty. While the ranges of yield uncertainties under CMIP5 and CMIP6 projections are similar, median impact in aggregate total caloric production is typically more negative for the CMIP6 projections (+1% to −19%) than for CMIP5 (+5% to −13%). In the first half of the 21st century and for individual crops is the spread across crop models typically wider than that across climate models, but we find distinct differences between crops: globally, wheat and maize uncertainties are dominated by the crop models, but soybean and rice are more sensitive to the climate projections. Climate models with very similar global mean warming can lead to very different aggregate impacts so that climate model uncertainties remain a significant contributor to agricultural impacts uncertainty. These results show the utility of large-ensemble methods that allow comprehensively evaluating factors affecting crop yields or other impacts under climate change. The crop model ensemble used here is unbalanced and pulls the assumption that all projections are equally plausible into question. Better methods for consistent model testing, also at the level of individual processes, will have to be developed and applied by the crop modeling community

Direct soil moisture controls of future global soil carbon changes: An important source of uncertainty

The nature of the climate–carbon cycle feedback depends critically on the response of soil carbon to climate, including changes in moisture. However, soil moisture–carbon feedback responses have not been investigated thoroughly. Uncertainty in the response of soil carbon to soil moisture changes could arise from uncertainty in the relationship between soil moisture and heterotrophic respiration. We used twelve soil moisture–respiration functions (SMRFs) with a soil carbon model (RothC) and data from a coupled climate–carbon cycle general circulation model to investigate the impact of direct heterotrophic respiration dependence on soil moisture on the climate carbon cycle feedback. Global changes in soil moisture acted to oppose temperature‐driven decreases in soil carbon and hence tended to increase soil carbon storage. We found considerable uncertainty in soil carbon changes due to the response of soil respiration to soil moisture. The use of different SMRFs resulted in both large losses and small gains in future global soil carbon stocks, whether considering all climate forcings or only moisture changes. Regionally, the greatest range in soil carbon changes across SMRFs was found where the largest soil carbon changes occurred. Further research is needed to constrain the soil moisture–respiration relationship and thus reduce uncertainty in climate–carbon cycle feedbacks. There may also be considerable uncertainty in the regional responses of soil carbon to soil moisture changes since climate model predictions of regional soil moisture changes are less coherent than temperature changes

How can the First ISLSCP Field Experiment contribute to present-day efforts to evaluate water stress in JULESv5.0?

The First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE), Kansas, US, 1987–1989, made important contributions to the understanding of energy and CO2 exchanges between the land surface and the atmosphere, which heavily influenced the development of numerical land-surface modelling. Now, 30 years on, we demonstrate how the wealth of data collected during FIFE and its subsequent in-depth analysis in the literature continue to be a valuable resource for the current generation of land-surface models. To illustrate, we use the FIFE dataset to evaluate the representation of water stress on tallgrass prairie vegetation in the Joint UK Land Environment Simulator (JULES) and highlight areas for future development. We show that, while JULES is able to simulate a decrease in net carbon assimilation and evapotranspiration during a dry spell, the shape of the diurnal cycle is not well captured. Evaluating the model parameters and results against this dataset provides a case study on the assumptions in calibrating “unstressed” vegetation parameters and thresholds for water stress. In particular, the responses to low water availability and high temperatures are calibrated separately. We also illustrate the effect of inherent uncertainties in key observables, such as leaf area index, soil moisture and soil properties. Given these valuable lessons, simulations for this site will be a key addition to a compilation of simulations covering a wide range of vegetation types and climate regimes, which will be used to improve the way that water stress is represented within JULES

Carbon residence time dominates uncertainty in terrestrial vegetation responses to future climate and atmospheric CO2.

Future climate change and increasing atmospheric CO2 are expected to cause major changes in vegetation structure and function over large fractions of the global land surface. Seven global vegetation models are used to analyze possible responses to future climate simulated by a range of general circulation models run under all four representative concentration pathway scenarios of changing concentrations of greenhouse gases. All 110 simulations predict an increase in global vegetation carbon to 2100, but with substantial variation between vegetation models. For example, at 4 °C of global land surface warming (510-758 ppm of CO2), vegetation carbon increases by 52-477 Pg C (224 Pg C mean), mainly due to CO2 fertilization of photosynthesis. Simulations agree on large regional increases across much of the boreal forest, western Amazonia, central Africa, western China, and southeast Asia, with reductions across southwestern North America, central South America, southern Mediterranean areas, southwestern Africa, and southwestern Australia. Four vegetation models display discontinuities across 4 °C of warming, indicating global thresholds in the balance of positive and negative influences on productivity and biomass. In contrast to previous global vegetation model studies, we emphasize the importance of uncertainties in projected changes in carbon residence times. We find, when all seven models are considered for one representative concentration pathway × general circulation model combination, such uncertainties explain 30% more variation in modeled vegetation carbon change than responses of net primary productivity alone, increasing to 151% for non-HYBRID4 models. A change in research priorities away from production and toward structural dynamics and demographic processes is recommended.The research leading to these results has received funding from the European Community’s Seventh Framework Programme (FP7 2007-2013) under Grant 238366. R.B., R.K., R.D., A.W., and P.D.F. were supported by the Joint Department of Energy and Climate Change/Department for Environment, Food and Rural Affairs Met Office Hadley Centre Climate Programme (GA01101). A.I. and K.N. were supported by the Environment Research and Technology Development Fund (S-10) of the Ministry of the Environment, Japan. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for the Coupled Model Intercomparison Project (CMIP), and we thank the climate modeling groups responsible for the GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, and NorESM1-M models for producing and making available their model output. For CMIP, the US Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. This work has been conducted under the framework of the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP). The ISI-MIP Fast Track project was funded by the German Federal Ministry of Education and Research (BMBF) with project funding Reference 01LS1201A.This is the author accepted manuscript. The final version is available from PNAS via http://dx.doi.org/10.1073/pnas.122247711

First look at changes in flood hazard in the Inter-Sectoral Impact Model Intercomparison Project ensemble

Climate change due to anthropogenic greenhouse gas emissions is expected to increase the frequency and intensity of precipitation events, which is likely to affect the probability of flooding into the future. In this paper we use river flow simulations from nine global hydrology and land surface models to explore uncertainties in the potential impacts of climate change on flood hazard at global scale. As an indicator of flood hazard we looked at changes in the 30-y return level of 5-d average peak flows under representative concentration pathway RCP8.5 at the end of this century. Not everywhere does climate change result in an increase in flood hazard: decreases in the magnitude and frequency of the 30-y return level of river flow occur at roughly one-third (20–45%) of the global land grid points, particularly in areas where the hydrograph is dominated by the snowmelt flood peak in spring. In most model experiments, however, an increase in flooding frequency was found in more than half of the grid points. The current 30-y flood peak is projected to occur in more than 1 in 5 y across 5–30% of land grid points. The large-scale patterns of change are remarkably consistent among impact models and even the driving climate models, but at local scale and in individual river basins there can be disagreement even on the sign of change, indicating large modeling uncertainty which needs to be taken into account in local adaptation studies

Climate change impacts on global agriculture

Based on predicted changes in the magnitude and distribution of global precipitation, temperature and river flow under the IPCC SRES A1B and A2 scenarios, this study assesses the potential impacts of climate change and CO2 fertilization on global agriculture. The analysis uses the new version of the GTAP-W model, which distinguishes between rainfed and irrigated agriculture and implements water as an explicit factor of production for irrigated agriculture. Future climate change is likely to modify regional water endowments and soil moisture. As a consequence, the distribution of harvested land will change, modifying production and international trade patterns. The results suggest that a partial analysis of the main factors through which climate change will affect agricultural productivity provide a false appreciation of the nature of changes likely to occur. Our results show that global food production, welfare and GDP fall in the two time periods and SRES scenarios. Higher food prices are expected. No matter which SRES scenario is preferred, we find that the expected losses in welfare are significant. These losses are slightly larger under the SRES A2 scenario for the 2020s and under the SRES A1B scenario for the 2050s. The results show that national welfare is influenced both by regional climate change and climate-induced changes in competitiveness

The relevance of uncertainty in future crop production for mitigation strategy planning

In order to achieve climate change mitigation, long-term decisions are required that must be reconciled with other societal goals that draw on the same resources. For example, ensuring food security for a growing population may require an expansion of crop land, thereby reducing natural carbon sinks or the area available for bio-energy production. Here, we show that current impact-model uncertainties pose an important challenge to long-term mitigation planning and propose a new risk-assessment and decision framework that accounts for competing interests. Based on cross-sectorally consistent simulations generated within the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP) we discuss potential gains and limitations of additional irrigation and trade-offs of the expansion of agricultural land as two possible response measures to climate change and growing food demand. We describe an illustrative example in which the combination of both measures may close the supply demand gap while leading to a loss of approximately half of all natural carbon sinks. We highlight current limitations of available simulations and additional steps required for a comprehensive risk assessment

Simulation of salinity effects on past, present, and future soil organic carbon stocks

Soil organic carbon (SOC) models are used to predict changes in SOC stocks and carbon dioxide (CO(2)) emissions from soils, and have been successfully validated for non-saline soils. However, SOC models have not been developed to simulate SOC turnover in saline soils. Due to the large extent of salt-affected areas in the world, it is important to correctly predict SOC dynamics in salt-affected soils. To close this knowledge gap, we modified the Rothamsted Carbon Model (RothC) to simulate SOC turnover in salt-affected soils, using data from non-salt-affected and salt-affected soils in two agricultural regions in India (120 soils) and in Australia (160 soils). Recently we developed a decomposition rate modifier based on an incubation study of a subset of these soils. In the present study, we introduce a new method to estimate the past losses of SOC due to salinity and show how salinity affects future SOC stocks on a regional scale. Because salinity decreases decomposition rates, simulations using the decomposition rate modifier for salinity suggest an accumulation of SOC. However, if the plant inputs are also adjusted to reflect reduced plant growth under saline conditions, the simulations show a significant loss of soil carbon in the past due to salinization, with a higher average loss of SOC in Australian soils (55 t C ha(-1)) than in Indian soils (31 t C ha(-1)). There was a significant negative correlation (p < 0.05) between SOC loss and osmotic potential. Simulations of future SOC stocks with the decomposition rate modifier and the plant input modifier indicate a greater decrease in SOC in saline than in non-saline soils under future climate. The simulations of past losses of SOC due to salinity were repeated using either measured charcoal-C or the inert organic matter predicted by the Falloon et al. equation to determine how much deviation from the Falloon et al. equation affects the amount of plant inputs generated by the model for the soils used in this study. Both sets of results suggest that saline soils have lost carbon and will continue to lose carbon under future climate. This demonstrates the importance of both reduced decomposition and reduced plant input in simulations of future changes in SOC stocks in saline soils.Raj Setia, Pete Smith, Petra Marschner, Pia Gottschalk, Jeff Baldock, Vipan Verma, Deepika Setia and Jo Smit