Global Gridded Crop Model Intercomparison Research Articles

The statistical emulators of GGCMI phase 2: responses of year-to-year variation of crop yield to CO2, temperature, water, and nitrogen perturbations

Abstract. Understanding the impact of climate change on year-to-year variation of crop yield is critical to global food stability and security. While crop model emulators are believed to be lightweight tools to replace the models, few emulators have been developed to capture such interannual variation of crop yield in response to climate variability. In this study, we developed a statistical emulator with a machine learning algorithm to reproduce the response of year-to-year variation of four crop yields to CO2 (C), temperature (T), water (W), and nitrogen (N) perturbations defined in the Global Gridded Crop Model Intercomparison Project (GGCMI) phase 2. The emulators were able to explain more than 52 % of the variance of simulated yield and performed well in capturing the year-to-year variation of global average and gridded crop yield over current croplands in the baseline. With the changes in CO2–temperature–water–nitrogen (CTWN) perturbations, the emulators could reproduce the year-to-year variation of crop yield well over most current cropland. The variation of R and the mean absolute error was small under the single CTWN perturbations and dual-factor perturbations. These emulators thus provide statistical response surfaces of yield, including both its mean and interannual variability, to climate factors. They could facilitate spatiotemporal downscaling of crop model simulation, projecting the changes in crop yield variability in the future and serving as a lightweight tool for multi-model ensemble simulation. The emulators enhanced the flexibility of crop yield estimates and expanded the application of large-ensemble simulations of crop yield under climate change.

Evaluation of global gridded crop models (GGCMs) for the simulation of major grain crop yields in China

Abstract Multimodel ensembles are powerful tools for evaluating agricultural production. Multimodel simulation results provided by the Global Gridded Crop Model Intercomparison (GGCMI) facilitate the evaluation of the grain production situation in China. With census crop yield data, the performance of nine global gridded crop models (GGCMs) in China was evaluated, and the yield gaps of four crops (maize, rice, soybean, and wheat) were estimated. The results showed that GGCMs better simulated maize yields than those of other crops in the northeast, north, northwest, east, and center. GEPIC (CLM-CROP) performed best in simulating maize (wheat) yield in the north, east, and northwest (southwest and south), due to reasonable parameter (cultivar and phenology parameters) settings. Because the rice phenology parameters were calibrated against phenological observation networks and a simple nitrogen limitation index was introduced, ORCHIDEE-CROP performed well in rice yield simulation and soybean yield simulation (center and southwest). Among four crops, wheat has the largest yield gap (7.3–14.1%), in which the poor soil of northwest (14.1%) exposes wheat to relatively high nutritional stress. Thus, in northwest China, optimizing nitrogen management in wheat production can effectively mitigate the negative impact of climate change on crop production.

Climate impacts on global agriculture emerge earlier in new generation of climate and crop models

Potential climate-related impacts on future crop yield are a major societal concern. Previous projections of the Agricultural Model Intercomparison and Improvement Project's Global Gridded Crop Model Intercomparison based on the Coupled Model Intercomparison Project Phase 5 identified substantial climate impacts on all major crops, but associated uncertainties were substantial. Here we report new twenty-first-century projections using ensembles of latest-generation crop and climate models. Results suggest markedly more pessimistic yield responses for maize, soybean and rice compared to the original ensemble. Mean end-of-century maize productivity is shifted from +5% to -6% (SSP126) and from +1% to -24% (SSP585)-explained by warmer climate projections and improved crop model sensitivities. In contrast, wheat shows stronger gains (+9% shifted to +18%, SSP585), linked to higher CO2 concentrations and expanded high-latitude gains. The 'emergence' of climate impacts consistently occurs earlier in the new projections-before 2040 for several main producing regions. While future yield estimates remain uncertain, these results suggest that major breadbasket regions will face distinct anthropogenic climatic risks sooner than previously anticipated.

Agricultural breadbaskets shift poleward given adaptive farmer behavior under climate change.

Modern food production is spatially concentrated in global "breadbaskets." A major unresolved question is whether these peak production regions will shift poleward as the climate warms, allowing some recovery of potential climate-related losses. While agricultural impacts studies to date have focused on currently cultivated land, the Global Gridded Crop Model Intercomparison Project (GGCMI) Phase 2 experiment allows us to assess changes in both yields and the location of peak productivity regions under warming. We examine crop responses under projected end of century warming using seven process-based models simulating five major crops (maize, rice, soybeans, and spring and winter wheat) with a variety of adaptation strategies. We find that in no-adaptation cases, when planting date and cultivar choices are held fixed, regions of peak production remain stationary and yield losses can be severe, since growing seasons contract strongly with warming. When adaptations in management practices are allowed (cultivars that retain growing season length under warming and modified planting dates), peak productivity zones shift poleward and yield losses are largely recovered. While most growing-zone shifts are ultimately limited by geography, breadbaskets studied here move poleward over 600km on average by end of the century under RCP 8.5. These results suggest that agricultural impacts assessments can be strongly biased if restricted in spatial area or in the scope of adaptive behavior considered. Accurate evaluation of food security under climate change requires global modeling and careful treatment of adaptation strategies.

Projections of thermal growing season indices over China under global warming of 1.5 °C and 2.0 °C

Global warming may prolong and intensify the thermal growing season of vegetation. It is not yet clear how the Paris Agreement's long-term temperature goals will affect the thermal growing season of vegetation, particularly crops, in China. Based on the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) datasets and the Global Gridded Crop Model Intercomparison (GGCMI) phase 1 dataset, we have investigated changes in spatiotemporal patterns of four thermal growing season indices (Growing Degree Days, GDD; Length of Growing Season, GSL; the Start of Growing Season, GSS; the End of Growing Season, GSE) over China under global warming scenarios of 1.5 °C and 2.0 °C with four representative concentration pathway (RCP) scenarios. Our results indicate that during the periods which achieve the global warming of 1.5 °C and 2.0 °C, only 3.82% and 29.15% of the total areas in China have higher warming levels beyond the global warming targets. For warmer RCP scenarios (except RCP2.6), there was a rising trend for GSE, GDD and GSL and a decreasing trend for GSS in China. Many crop regions in China have also shown an advance of GSS, an extension of GSL and an earlier end of GSE under the global warming of 1.5 °C and 2.0 °C, suggesting that crop planting and harvesting dates need to be adjusted accordingly in order to capture appropriate timing for crop maturity and to achieve a maximum yield.

Strong regional influence of climatic forcing datasets on global crop model ensembles

We present results from the Agricultural Model Intercomparison and Improvement Project (AgMIP) Global Gridded Crop Model Intercomparison (GGCMI) Phase I, which aligned 14 global gridded crop models (GGCMs) and 11 climatic forcing datasets (CFDs) in order to understand how the selection of climate data affects simulated historical crop productivity of maize, wheat, rice and soybean. Results show that CFDs demonstrate mean biases and differences in the probability of extreme events, with larger uncertainty around extreme precipitation and in regions where observational data for climate and crop systems are scarce. Countries where simulations correlate highly with reported FAO national production anomalies tend to have high correlations across most CFDs, whose influence we isolate using multi-GGCM ensembles for each CFD. Correlations compare favorably with the climate signal detected in other studies, although production in many countries is not primarily climate-limited (particularly for rice). Bias-adjusted CFDs most often were among the highest model-observation correlations, although all CFDs produced the highest correlation in at least one top-producing country. Analysis of larger multi-CFD-multi-GGCM ensembles (up to 91 members) shows benefits over the use of smaller subset of models in some regions and farming systems, although bigger is not always better. Our analysis suggests that global assessments should prioritize ensembles based on multiple crop models over multiple CFDs as long as a top-performing CFD is utilized for the focus region.

The GGCMI Phase 2 experiment: global gridded crop model simulations under uniform changes in CO<sub>2</sub>, temperature, water, and nitrogen levels (protocol version 1.0)

Abstract. Concerns about food security under climate change motivate efforts to better understand future changes in crop yields. Process-based crop models, which represent plant physiological and soil processes, are necessary tools for this purpose since they allow representing future climate and management conditions not sampled in the historical record and new locations to which cultivation may shift. However, process-based crop models differ in many critical details, and their responses to different interacting factors remain only poorly understood. The Global Gridded Crop Model Intercomparison (GGCMI) Phase 2 experiment, an activity of the Agricultural Model Intercomparison and Improvement Project (AgMIP), is designed to provide a systematic parameter sweep focused on climate change factors and their interaction with overall soil fertility, to allow both evaluating model behavior and emulating model responses in impact assessment tools. In this paper we describe the GGCMI Phase 2 experimental protocol and its simulation data archive. A total of 12 crop models simulate five crops with systematic uniform perturbations of historical climate, varying CO2, temperature, water supply, and applied nitrogen (“CTWN”) for rainfed and irrigated agriculture, and a second set of simulations represents a type of adaptation by allowing the adjustment of growing season length. We present some crop yield results to illustrate general characteristics of the simulations and potential uses of the GGCMI Phase 2 archive. For example, in cases without adaptation, modeled yields show robust decreases to warmer temperatures in almost all regions, with a nonlinear dependence that means yields in warmer baseline locations have greater temperature sensitivity. Inter-model uncertainty is qualitatively similar across all the four input dimensions but is largest in high-latitude regions where crops may be grown in the future.

Parameterization-induced uncertainties and impacts of crop management harmonization in a global gridded crop model ensemble

Global gridded crop models (GGCMs) combine agronomic or plant growth models with gridded spatial input data to estimate spatially explicit crop yields and agricultural externalities at the global scale. Differences in GGCM outputs arise from the use of different biophysical models, setups, and input data. GGCM ensembles are frequently employed to bracket uncertainties in impact studies without investigating the causes of divergence in outputs. This study explores differences in maize yield estimates from five GGCMs based on the public domain field-scale model Environmental Policy Integrated Climate (EPIC) that participate in the AgMIP Global Gridded Crop Model Intercomparison initiative. Albeit using the same crop model, the GGCMs differ in model version, input data, management assumptions, parameterization, and selection of subroutines affecting crop yield estimates via cultivar distributions, soil attributes, and hydrology among others. The analyses reveal inter-annual yield variability and absolute yield levels in the EPIC-based GGCMs to be highly sensitive to soil parameterization and crop management. All GGCMs show an intermediate performance in reproducing reported yields with a higher skill if a static soil profile is assumed or sufficient plant nutrients are supplied. An in-depth comparison of setup domains for two EPIC-based GGCMs shows that GGCM performance and plant stress responses depend substantially on soil parameters and soil process parameterization, i.e. hydrology and nutrient turnover, indicating that these often neglected domains deserve more scrutiny. For agricultural impact assessments, employing a GGCM ensemble with its widely varying assumptions in setups appears the best solution for coping with uncertainties from lack of comprehensive global data on crop management, cultivar distributions and coefficients for agro-environmental processes. However, the underlying assumptions require systematic specifications to cover representative agricultural systems and environmental conditions. Furthermore, the interlinkage of parameter sensitivity from various domains such as soil parameters, nutrient turnover coefficients, and cultivar specifications highlights that global sensitivity analyses and calibration need to be performed in an integrated manner to avoid bias resulting from disregarded core model domains. Finally, relating evaluations of the EPIC-based GGCMs to a wider ensemble based on individual core models shows that structural differences outweigh in general differences in configurations of GGCMs based on the same model, and that the ensemble mean gains higher skill from the inclusion of structurally different GGCMs. Although the members of the wider ensemble herein do not consider crop-soil-management interactions, their sensitivity to nutrient supply indicates that findings for the EPIC-based sub-ensemble will likely become relevant for other GGCMs with the progressing inclusion of such processes.

The Global Gridded Crop Model Intercomparison phase 1 simulation dataset

The Global Gridded Crop Model Intercomparison (GGCMI) phase 1 dataset of the Agricultural Model Intercomparison and Improvement Project (AgMIP) provides an unprecedentedly large dataset of crop model simulations covering the global ice-free land surface. The dataset consists of annual data fields at a spatial resolution of 0.5 arc-degree longitude and latitude. Fourteen crop modeling groups provided output for up to 11 historical input datasets spanning 1901 to 2012, and for up to three different management harmonization levels. Each group submitted data for up to 15 different crops and for up to 14 output variables. All simulations were conducted for purely rainfed and near-perfectly irrigated conditions on all land areas irrespective of whether the crop or irrigation system is currently used there. With the publication of the GGCMI phase 1 dataset we aim to promote further analyses and understanding of crop model performance, potential relationships between productivity and environmental impacts, and insights on how to further improve global gridded crop model frameworks. We describe dataset characteristics and individual model setup narratives.

Global patterns of crop yield stability under additional nutrient and water inputs.

Agricultural production must increase to feed a growing and wealthier population, as well as to satisfy increasing demands for biomaterials and biomass-based energy. At the same time, deforestation and land-use change need to be minimized in order to preserve biodiversity and maintain carbon stores in vegetation and soils. Consequently, agricultural land use needs to be intensified in order to increase food production per unit area of land. Here we use simulations of AgMIP’s Global Gridded Crop Model Intercomparison (GGCMI) phase 1 to assess implications of input-driven intensification (water, nutrients) on crop yield and yield stability, which is an important aspect in food security. We find region- and crop-specific responses for the simulated period 1980–2009 with broadly increasing yield variability under additional nitrogen inputs and stabilizing yields under additional water inputs (irrigation), reflecting current patterns of water and nutrient limitation. The different models of the GGCMI ensemble show similar response patterns, but model differences warrant further research on management assumptions, such as variety selection and soil management, and inputs as well as on model implementation of different soil and plant processes, such as on heat stress, and parameters. Higher variability in crop productivity under higher fertilizer input will require adequate buffer mechanisms in trade and distribution/storage networks to avoid food price volatility.

Yield risks in global maize markets: Historical evidence and projections in key regions of the world

Simultaneous worldwide crop failures stemming from a more unstable climate may reduce the scope for international trade to compensate food shortages and stabilize food prices across the various regions of the world. Understanding the effects of changes in crop productivity on global markets requires knowledge about the extent to which crop yields may be systematically related across producing and consuming centers. This short communication contributes to this knowledge by investigating the potential changes in the strength of two key sources of supply risks in global maize markets: yield variance and cross-country yield correlation. We focus on the largest producing and consuming countries of the world. We capitalize on yield projections from the Global Gridded Crop Model Intercomparison project. Exploratory analysis of the skill of the underlying GGCMI models in reproducing key moments of the distribution of observed yields reveals that they overstate observed variances but faithfully reproduce observed patterns of cross-country correlations. We find no evidence of an increase in the degree of cross-country dependency of maize yields. We also find a higher incidence of what would be considered extremely low maize yields by present-time standards stemming from the projected downward trend in yields levels toward mid-century. The weak dependency of maize yields across countries, an the possibility of reducing the higher incidence of extremes through policies aimed to reverse the climate-induced downward trends in yields, suggest that international trade can become a valuable tool to ameliorate the effects of more unstable crop yields.

Global gridded crop model evaluation: benchmarking, skills, deficiencies and implications

Abstract. Crop models are increasingly used to simulate crop yields at the global scale, but so far there is no general framework on how to assess model performance. Here we evaluate the simulation results of 14 global gridded crop modeling groups that have contributed historic crop yield simulations for maize, wheat, rice and soybean to the Global Gridded Crop Model Intercomparison (GGCMI) of the Agricultural Model Intercomparison and Improvement Project (AgMIP). Simulation results are compared to reference data at global, national and grid cell scales and we evaluate model performance with respect to time series correlation, spatial correlation and mean bias. We find that global gridded crop models (GGCMs) show mixed skill in reproducing time series correlations or spatial patterns at the different spatial scales. Generally, maize, wheat and soybean simulations of many GGCMs are capable of reproducing larger parts of observed temporal variability (time series correlation coefficients (r) of up to 0.888 for maize, 0.673 for wheat and 0.643 for soybean at the global scale) but rice yield variability cannot be well reproduced by most models. Yield variability can be well reproduced for most major producing countries by many GGCMs and for all countries by at least some. A comparison with gridded yield data and a statistical analysis of the effects of weather variability on yield variability shows that the ensemble of GGCMs can explain more of the yield variability than an ensemble of regression models for maize and soybean, but not for wheat and rice. We identify future research needs in global gridded crop modeling and for all individual crop modeling groups. In the absence of a purely observation-based benchmark for model evaluation, we propose that the best performing crop model per crop and region establishes the benchmark for all others, and modelers are encouraged to investigate how crop model performance can be increased. We make our evaluation system accessible to all crop modelers so that other modeling groups can also test their model performance against the reference data and the GGCMI benchmark.

Spatial and temporal uncertainty of crop yield aggregations

Abstract The aggregation of simulated gridded crop yields to national or regional scale requires information on temporal and spatial patterns of crop-specific harvested areas. This analysis estimates the uncertainty of simulated gridded yield time series related to the aggregation with four different harvested area data sets. We compare aggregated yield time series from the Global Gridded Crop Model Intercomparison project for four crop types from 14 models at global, national, and regional scale to determine aggregation-driven differences in mean yields and temporal patterns as measures of uncertainty. The quantity and spatial patterns of harvested areas differ for individual crops among the four data sets applied for the aggregation. Also simulated spatial yield patterns differ among the 14 models. These differences in harvested areas and simulated yield patterns lead to differences in aggregated productivity estimates, both in mean yield and in the temporal dynamics. Among the four investigated crops, wheat yield (17% relative difference) is most affected by the uncertainty introduced by the aggregation at the global scale. The correlation of temporal patterns of global aggregated yield time series can be as low as for soybean (r = 0.28). For the majority of countries, mean relative differences of nationally aggregated yields account for 10% or less. The spatial and temporal difference can be substantial higher for individual countries. Of the top-10 crop producers, aggregated national multi-annual mean relative difference of yields can be up to 67% (maize, South Africa), 43% (wheat, Pakistan), 51% (rice, Japan), and 427% (soybean, Bolivia). Correlations of differently aggregated yield time series can be as low as r = 0.56 (maize, India), r = 0.05 (wheat, Russia), r = 0.13 (rice, Vietnam), and r = −0.01 (soybean, Uruguay). The aggregation to sub-national scale in comparison to country scale shows that spatial uncertainties can cancel out in countries with large harvested areas per crop type. We conclude that the aggregation uncertainty can be substantial for crop productivity and production estimations in the context of food security, impact assessment, and model evaluation exercises.

The Global Gridded Crop Model Intercomparison: data and modeling protocols for Phase 1 (v1.0)

Abstract. We present protocols and input data for Phase 1 of the Global Gridded Crop Model Intercomparison, a project of the Agricultural Model Intercomparison and Improvement Project (AgMIP). The project includes global simulations of yields, phenologies, and many land-surface fluxes using 12–15 modeling groups for many crops, climate forcing data sets, and scenarios over the historical period from 1948 to 2012. The primary outcomes of the project include (1) a detailed comparison of the major differences and similarities among global models commonly used for large-scale climate impact assessment, (2) an evaluation of model and ensemble hindcasting skill, (3) quantification of key uncertainties from climate input data, model choice, and other sources, and (4) a multi-model analysis of the agricultural impacts of large-scale climate extremes from the historical record.

Assessing agricultural risks of climate change in the 21st century in a global gridded crop model intercomparison

Here we present the results from an intercomparison of multiple global gridded crop models (GGCMs) within the framework of the Agricultural Model Intercomparison and Improvement Project and the Inter-Sectoral Impacts Model Intercomparison Project. Results indicate strong negative effects of climate change, especially at higher levels of warming and at low latitudes; models that include explicit nitrogen stress project more severe impacts. Across seven GGCMs, five global climate models, and four representative concentration pathways, model agreement on direction of yield changes is found in many major agricultural regions at both low and high latitudes; however, reducing uncertainty in sign of response in mid-latitude regions remains a challenge. Uncertainties related to the representation of carbon dioxide, nitrogen, and high temperature effects demonstrated here show that further research is urgently needed to better understand effects of climate change on agricultural production and to devise targeted adaptation strategies.