Abstract
Increasing global food demand will require more food production1 without further exceeding the planetary boundaries2 while simultaneously adapting to climate change3. We used an ensemble of wheat simulation models with improved sink and source traits from the highest-yielding wheat genotypes4 to quantify potential yield gains and associated nitrogen requirements. This was explored for current and climate change scenarios across representative sites of major world wheat producing regions. The improved sink and source traits increased yield by 16% with current nitrogen fertilizer applications under both current climate and mid-century climate change scenarios. To achieve the full yield potential—a 52% increase in global average yield under a mid-century high warming climate scenario (RCP8.5), fertilizer use would need to increase fourfold over current use, which would unavoidably lead to higher environmental impacts from wheat production. Our results show the need to improve soil nitrogen availability and nitrogen use efficiency, along with yield potential.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 digital issues and online access to articles
$119.00 per year
only $9.92 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout



Similar content being viewed by others
Data availability
The measured and simulated crop data and model inputs for the New Zealand and South America field experiment data are available at https://doi.org/10.7910/DVN/XA4VA2 (ref. 54) and https://doi.org/10.7910/DVN/VKWKUP (ref. 55) respectively. The simulation protocols, model inputs and simulation results for the 34 global locations are available at https://doi.org/10.7910/DVN/6KBBI3 (ref. 56). National wheat production statistics and area for the period 2016 to 2019 were obtained from the FAO public database available at https://www.fao.org/faostat. Daily weather data were obtained from the AgMERRA (https://data.giss.nasa.gov/impacts/agmipcf/agmerra/) and NASA POWER (https://power.larc.nasa.gov/) climate datasets.
Code availability
The data analysis scripts were developed with R (v.4.1.3) and are available in GitHub at https://github.com/pmartre/AgMIPWheat4 (ref. 57). Documentation and codes of the crop models used in this study are available from the links or email addresses given in Supplementary Table 1.
References
Mottaleb, K. A., Kruseman, G., Frija, A., Sonder, K. & Lopez-Ridaura, S. Projecting wheat demand in China and India for 2030 and 2050: implications for food security. Front. Nutr. 9, 1077443 (2023).
de Vries, W., Kros, J., Kroeze, C. & Seitzinger, S. P. Assessing planetary and regional nitrogen boundaries related to food security and adverse environmental impacts. Curr. Opin. Environ. Sustain. 5, 392–402 (2013).
Asseng, S. et al. Climate change impact and adaptation for wheat protein. Glob. Change Biol. 25, 155–173 (2019).
Bustos, D. V., Hasan, A. K., Reynolds, M. P. & Calderini, D. F. Combining high grain number and weight through a DH-population to improve grain yield potential of wheat in high-yielding environments. Field Crops Res. 145, 106–115 (2013).
Erenstein, O. et al. in Wheat Improvement: Food Security in a Changing Climate (eds Reynolds, M. P. & Braun, H.-J.) Ch. 4 (Springer, 2022).
Godfray, H. C. J. et al. Food security: the challenge of feeding 9 billion people. Science 327, 812–818 (2010).
Chaudhary, A., Pfister, S. & Hellweg, S. Spatially explicit analysis of biodiversity loss due to global agriculture, pasture and forest land use from a producer and consumer perspective. Environ. Sci. Technol. 50, 3928–3936 (2016).
Breitburg, D. et al. Declining oxygen in the global ocean and coastal waters. Science 359, eaam7240 (2018).
Zhang, D., Shen, J., Zhang, F., Li, Y. & Zhang, W. Carbon footprint of grain production in China. Sci. Rep. 7, 4126 (2017).
Ladha, J. K., Pathak, H., Krupnik, T. J., Six, J. & van Kessel, C. Efficiency of fertilizer nitrogen in cereal production: retrospects and prospects. Adv. Agron. 87, 85–156 (2005).
Ladha, J. K. et al. Global nitrogen budgets in cereals: a 50-year assessment for maize, rice and wheat production systems. Sci. Rep. 6, 19355 (2016).
Zhang, X. et al. Sustainable nitrogen management index: definition, global assessment, and potential improvements. Front. Agric. Sci. Eng. 9, 356–365 (2022).
Nabuurs, G. J. et al. in Climate Change 2022: Mitigation of Climate Change (eds Shukla, P. R. et al.) Ch. 7 (IPCC, Cambridge Univ. Press, 2022).
Morris, M., Edmeades, G. & Pehu, E. The global need for plant breeding capacity: what roles for the public and private sectors? Hortscience 41, 30–39 (2006).
Reynolds, M. P. et al. A wiring diagram to integrate physiological traits of wheat yield potential. Nat. Food 3, 318–324 (2022).
Fischer, R. A. & Edmeades, G. O. Breeding and cereal yield progress. Crop Sci. 50, S-85–S-98 (2010).
van Ittersum, M. K. et al. Yield gap analysis with local to global relevance—a review. Field Crops Res. 143, 4–17 (2013).
Mueller, N. D. et al. Closing yield gaps through nutrient and water management. Nature 490, 254–257 (2012).
Senapati, N. et al. Global wheat production could benefit from closing the genetic yield gap. Nat. Food 3, 532–541 (2022).
Paleari, L. et al. A trait-based model ensemble approach to design rice plant types for future climate. Glob. Change Biol. 28, 2689–2710 (2022).
Molero, G. et al. Elucidating the genetic basis of biomass accumulation and radiation use efficiency in spring wheat and its role in yield potential. Plant Biotechnol. J. 17, 1276–1288 (2019).
van Grinsven, H. J. M. et al. Establishing long-term nitrogen response of global cereals to assess sustainable fertilizer rates. Nat. Food 3, 122–132 (2022).
Webber, H. et al. No perfect storm for crop yield failure in Germany. Environ. Res. Lett. 15, 104012 (2020).
Sapkota, T. B. et al. Identifying optimum rates of fertilizer nitrogen application to maximize economic return and minimize nitrous oxide emission from rice–wheat systems in the Indo-Gangetic Plains of India. Arch. Agron. Soil Sci. 66, 2039–2054 (2020).
Fast, A. et al. Integrating enhanced efficiency fertilizers and nitrogen rates to improve Canada Western Red Spring wheat. Can. J. Plant. Sci. 104, 144–160 (2023).
Pan, W. L., Kidwell, K. K., McCracken, V. A., Bolton, R. P. & Allen, M. Economically optimal wheat yield, protein and nitrogen use component responses to varying N supply and genotype. Front. Plant Sci. 10, 1790 (2020).
Ma, G. et al. Determining the optimal N input to improve grain yield and quality in winter wheat with reduced apparent N loss in the North China Plain. Front. Plant Sci. 10, 181 (2019).
Savin, R., Sadras, V. O. & Slafer, G. A. Benchmarking nitrogen utilisation efficiency in wheat for Mediterranean and non-Mediterranean European regions. Field Crops Res. 241, 107573 (2019).
Ludemann, C. I., Gruere, A., Heffer, P. & Dobermann, A. Global data on fertilizer use by crop and by country. Sci. Data 9, 501 (2022).
Bonilla-Cedrez, C., Chamberlin, J. & Hijmans, R. J. Fertilizer and grain prices constrain food production in sub-Saharan Africa. Nat. Food 2, 766–772 (2021).
Bouwman, A. F. et al. Lessons from temporal and spatial patterns in global use of N and P fertilizer on cropland. Sci. Rep. 7, 40366 (2017).
Tilman, D., Balzer, C., Hill, J. & Befort, B. L. Global food demand and the sustainable intensification of agriculture. Proc. Natl Acad. Sci. USA 108, 20260–20264 (2011).
Döring, T. F. & Neuhoff, D. Upper limits to sustainable organic wheat yields. Sci. Rep. 11, 12729 (2021).
Subbarao, G. V. et al. Enlisting wild grass genes to combat nitrification in wheat farming: a nature-based solution. Proc. Natl Acad. Sci. USA 118, e2106595118 (2021).
Kanter, D. R., Zhang, X. & Mauzerall, D. L. Reducing nitrogen pollution while decreasing farmers’ costs and increasing fertilizer industry profits. J. Environ. Qual. 44, 325–335 (2015).
Lemaire, G., Sinclair, T., Sadras, V. & Bélanger, G. Allometric approach to crop nutrition and implications for crop diagnosis and phenotyping. A review. Agron. Sustain. Dev. 39, 27 (2019).
Pingali, P. L. Green Revolution: impacts, limits, and the path ahead. Proc. Natl Acad. Sci. USA 109, 12302–12308 (2012).
Dueri, S. et al. Simulation of winter wheat response to variable sowing dates and densities in a high-yielding environment. J. Exp. Bot. 73, 5715–5729 (2022).
García, G. A. et al. Grain yield potential strategies in an elite wheat double-haploid population grown in contrasting environments. Crop Sci. 53, 2577–2587 (2013).
Liu, B. et al. Similar estimates of temperature impacts on global wheat yield by three independent methods. Nat. Clim. Change 6, 1130–1136 (2016).
Rattalino Edreira, J. I. et al. Spatial frameworks for robust estimation of yield gaps. Nat. Food 2, 773–779 (2021).
Ruane, A. C., Goldberg, R. & Chryssanthacopoulos, J. Climate forcing datasets for agricultural modeling: merged products for gap-filling and historical climate series estimation. Agric. Meteorol. 200, 233–248 (2015).
Ruane, A. C., Winter, J. M., Mcdermid, S. P. & Hudson, N. I. in Handbook of Climate Change and Agroecosystems: The Agricultural Model Intercomparison and Improvement Project (eds Rosenzweig, C. & Hillel, D.) 45–78 (Imperial College Press, 2015).
Ruane, A. C. & McDermid, S. P. Selection of a representative subset of global climate models that captures the profile of regional changes for integrated climate impacts assessment. Earth Perspect. 4, 1 (2017).
Taylor, K., Stouffer, R. & Meehl, G. An overview of CMIP5 and the experiment design. Bull. Am. Meteorol. Soc. 93, 485–498 (2012).
Müller, C. et al. Exploring uncertainties in global crop yield projections in a large ensemble of crop models and CMIP5 and CMIP6 climate scenarios. Environ. Res. Lett. https://doi.org/10.1088/1748-9326/abd8fc (2021).
Mossé, J., Huet, J. C. & Baudet, J. The amino acid composition of wheat grain as a function of nitrogen content. J. Cereal Sci. 3, 115–130 (1985).
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2022).
Lobell, D. B. Climate change adaptation in crop production: beware of illusions. Glob. Food Sec. 3, 72–76 (2014).
Heffer, P. Assessment of Fertilizer Use By Crop at the Global Level 2010–2010/11 (International Fertilizer Industry Association, 2013).
Heffer, P., Gruère, A. & Roberts, T. Assessment of Fertilizer Use By Crop at the Global Level 2014–2014/15 (International Fertilizer Association, International Plant Nutrition Institute, 2017).
Fertilizer Use by Crop 5th edn (FAO, IFDC, IPI, PPI, 2002).
Shahzad, A. N., Qureshi, M. K., Wakeel, A. & Misselbrook, T. Crop production in Pakistan and low nitrogen use efficiencies. Nat. Sustain. 2, 1106–1114 (2019).
Dueri, S. et al. Data from the winter wheat potential yield experiment in New Zealand and response to variable sowing dates and densities: field experiments and AgMIP-Wheat multi-model simulations. Harvard Dataverse https://doi.org/10.7910/DVN/XA4VA2 (2022).
Guarin, J. R. et al. Data from the AgMIP-Wheat high-yielding traits experiment for modeling potential production of wheat: field experiments and multi-model simulations. Harvard Dataverse https://doi.org/10.7910/DVN/VKWKUP (2022).
Martre, P. et al. Replication data for: global implications for nitrogen use of improved wheat yield under climate change. Harvard Dataverse https://doi.org/10.7910/DVN/6KBBI3 (2023).
Martre, P. AgMIPWheat4. GitHub https://github.com/pmartre/AgMIPWheat4 (2023).
Acknowledgements
This study was a part of the Agricultural Model Intercomparison and Improvement Project (AgMIP) Wheat Phase 4. The experimental work conducted at Valdivia, Chili by J. Herrera (UACh) is appreciated. P.M. and S.D. acknowledge support from the metaprogram Agriculture and forestry in the face of climate change: adaptation and mitigation (CLIMAE) of INRAE. This work was supported by the French National Research Institute for Agriculture, Food and Environment (INRAE); the International Maize and Wheat Improvement Center (CIMMYT) and the International Wheat Yield Partnership (IWYP, grant IWYP115 to P.M., S.A. and F.E.), CIMMYT and the Chilean Technical and Scientific Research Council (CONICYT-ANID) through FONDECYT (grant 1141048 to D. Calderini); the Foundation for Food and Agricultural Research (to M.R.); the German Federal Ministry of Education and Research (BMBF) through the BonaRes project ‘I4S’ (grant 031B0513I to K.C.K.); the Ministry of Education, Youth and Sports of Czech Republic through SustES - Adaption strategies for sustainable ecosystem services and food security under adverse environmental conditions (grant CZ.02.1.01/0.0/0.0/16_019/000797 to K.C.K. and C.N.); the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy (grant EXC 2070 – 390732324 to F.E. and T.G.) and the Collaborative Research Centre DETECT (grant No. SFB1502/1–2022 -450058266 to T.G.); the JPI-FACCE MACSUR2 project, funded by the Italian Ministry for Agricultural, Food and Forestry Policies (grant 24064/7303/15 to R.F. and G.P.) and the SYSTEMIC project funded by JPI-HDHL, JPI-OCEANS and FACCE-JPI under ERA-NET (grant 696295 to R.F. and G.P.); and BMBF in the framework of the funding measure ‘Soil as a Sustainable Resource for the Bioeconomy—BonaRes’, project BonaRes (Module A): BonaRes Center for Soil Research, subproject ‘Sustainable Subsoil Management—Soil3’ (grant 031B0151A to A.K.S.) and COINS (grant 01LL2204C to A.K.S.). A.C.R. received support from the National Aeronautics and Space Administration (NASA) Earth Science Division grant for the NASA Goddard Institute for Space Studies Climate Impacts Group. J.-P.C. and J.-C.D. received support from the CASDAR and Intercéréales funds.
Author information
Authors and Affiliations
Contributions
P.M., S.A., F.E. and H.W. conceived the research. P.M. designed the study. D. Calderini, G.M., M.R., D.M., G.G., H.B., M.G., R.C. J.-P.C. and J.-C.D. conducted the field experiments. P.M., S.D. and J.R.G. analysed the data. P.M. produced the figures and wrote the first draft of the manuscript. All authors contributed to the revision of the manuscript. Authors D. Cammarano, R.F., T.G., Y.G., Z.H., G.H., L.A.H., K.C.K., C.N., G.P., A.C.R., A.K.S., T.S., I.S., P.T., E.W., J.W., C.Z., M.B. and Z.Z. performed the crop model simulations and are listed in alphabetical order.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Plants thanks Ebrahim Jahanshiri and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Simulated global nitrogen (N) unlimited grain yield and protein concentration for cultivars with and without high-yielding traits.
Simulated global average N unlimited (a) grain yield and grain protein concentration (b) for locally adapted cultivars without (open boxes) and with (closed boxes) high-yielding traits for the 1981–2010 baseline period (yellow) and the 2040–2069 period under RCP4.5 (pink) and RCP8.5 (blue For each box plot, the vertical bars represent 1.5 times the interquartile range, the box represents the interquartile, and the horizontal lines inside the box represents the median. Data are the average of 30 years of yield simulated with 12 wheat crop growth models, and for future climate scenarios of five global climate models, using region-specific soils, cultivars, and sowing data and unlimited N supply.
Extended Data Fig. 2 Comparison of measured and simulated yield, yield components, and phenology of the spring wheat cultivar Bacanora and the best yielding line of the Bacanora x Weebil cross.
Measured and simulated (a) grain yield, (b) total above ground biomass, (c) harvest index, (d) grain number, (e) average grain dry mass, (f) days between sowing anthesis, and (g) grain filling duration for the spring wheat cultivar Bacanora and the best yielding doubled haploid line of the Bacanora x Weebil cross sown in the field in Valdivia, Chile (VA) in 2008 and 2009, in Buenos Aires, Argentina (BU) in 2009, and in Ciudad Obregon, Mexico (CO) in 2009, and for the average of the four experimental years (Average). Data are the average of 3 to 4 independent replicates for each year and location of measurements and the ensemble median of 12 wheat crop growth models.
Extended Data Fig. 3 Simulated relative climate change adaptation for global average wheat production, nitrogen (N) demand, and grain protein concentration for high-yielding traits at different rates of N fertilizer application.
Adaptation to climate change of (a) global production, (b) N demand, and (c) grain protein concentration for mid-century RCP4.5 and RCP8.5 climate scenarios relative to the baseline average for 1981–2010. Adaptation to climate change is calculated at the average N fertilization rate currently used by farmers in each country (yellow bars) and at N rates that maximize grain yield at each location considering (blue bars) or not (pink bars) a minimum grain protein concentration of 12%. Data are for the average of 30 years of yield simulated by 12 wheat crop growth models and five climate models. Relative adaptation values calculated at 34 representative global sites were weighed and aggregated to the globe based on national wheat production data. For each box plot, the vertical bars represent 1.5 times the interquartile range, the box represents the interquartile, and the horizontal lines inside the box represents the median.
Extended Data Fig. 4 Simulated grain yield response to nitrogen (N) fertilizer rate at 34 representative high rainfall or irrigated global sites.
Simulations are for the baseline period (1981–2010) for locally adapted cultivars without (black circles) and with high-yielding traits given in Supplementary Table 2 (red circles), and for mid-century RCP4.5 (blue circles) and RCP8.5 (green circles) climate scenarios with the high-yielding trait only. Data are the average of 30 years simulations with 12 wheat crop growth models, plus five global climate models for the climate change scenarios, using region-specific soils and sowing data. Red crosses show the current national average grain yield against the current national average N fertilizer application rate from reported national data, respectively.
Extended Data Fig. 5 Simulated grain nitrogen (N) use efficiency response to N fertilizer rate at 34 representative high rainfall or irrigated global sites.
Simulations are for the baseline period (1981–2010) for locally adapted cultivars without (black circles) and with high yielding traits given in Supplementary Table 2 (red circles), and for mid-century RCP4.5 (blue circles) and RCP8.5 (green circles) climate scenarios with the high-yielding trait only. Data are the average of 30 years simulations with 12 wheat crop growth models, plus five global climate models for the climate change scenarios, using region-specific soils and sowing data.
Extended Data Fig. 6 Simulated grain protein concentration response to nitrogen (N) fertilizer rate at 34 representative high rainfall or irrigated global sites.
Simulations are for the baseline period (1981–2010) for locally adapted cultivars without (black circles) and with high-yielding traits given in Supplementary Table 2 (red circles), and for mid-century RCP4.5 (blue circles) and RCP8.5 (green circles) climate scenarios with the high-yielding traits only. Data are the average of 30 years simulations with 12 wheat crop growth models, plus five global climate models for the climate change scenarios, using region-specific soils and sowing data.
Extended Data Fig. 7 Comparison of simulated and historical national average grain yield and N fertilizer rate for 23 countries.
In (a) Simulated national average grain yield at each studied location was interpolated on the relationship between simulated grain yield and N fertilizer application rate (Extended Data Fig. 3) at the N fertilizer application rate equal to the reported national application rate. In (b) at each studied location N fertilizer rate was interpolated on the relationship between simulated grain yield and N fertilizer application rate at the grain yield equal to the reported national average grain yield. Simulated data are means of 30 years (1981-20210) simulations with an ensemble of 12 wheat crop growth models. Error bars show 50% of the models (interquartile). Countries are represented with their two-letter code. Historical national yields are the average for the 2016 to 2019 harvests. Dotted line is the 1:1 relationship and solid line is linear regression. RRMSE, NU, LC, and SB are relative root mean squared error, non-unity slope, lack or correlation, and squared bias, respectively. Note that systematic biases (overestimation of national average yields and underestimation of the national average rate of nitrogen fertilization) are expected due to the impact of pests, diseases and weeds on yield, or to negative responses to certain extreme climatic impact factors23 that are not taken into account in the crop models used in this study.
Extended Data Fig. 8 Simulated impact of high-yielding traits on grain protein concentration with and without adaptation of N fertilizer application rate at 34 representative global sites.
Simulations are for the 1981–2010 baseline period (yellow) and for mid-century RCP4.5 (pink) and RCP8.5 (blue) climate scenarios. The relative change in N use efficiency are calculated at the average N fertilization rate currently used by farmers in each country (yellow bars) and at N rates that maximize grain yield at each location considering (blue bars) or not (pink bars) a minimum grain protein concentration of 12%. Data are the average of 30 years of yield simulated with 12 wheat crop growth models and five climate models using region-specific soils and sowing data. For each box plot, the vertical bars represent 1.5 times the interquartile range, the box represents the interquartile, and the horizontal lines inside the box represents the median.
Extended Data Fig. 9 Simulated impact of high-yielding traits on grain yield with and without adaptation of nitrogen (N) fertilizer application rate at 34 representative global sites.
Simulations are for the 1981–2010 baseline period (yellow) and for mid-century RCP4.5 (pink) and RCP8.5 (blue) climate scenarios. The relative change in grain yield are calculated at the average N fertilization rate currently used by farmers in each country (yellow bars) and at N rates that maximize grain yield at each location considering (blue bars) or not (pink bars) a minimum grain protein concentration of 12%. Data are the average of 30 years of yield simulated with 12 wheat crop growth models and five climate models using region-specific soils and sowing data. For each box plot, the vertical bars represent 1.5 times the interquartile range, the box represents the interquartile, and the horizontal lines inside the box represents the median.
Supplementary information
Supplementary Information
Supplementary Tables 1–8 and Figs. 1–8.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Martre, P., Dueri, S., Guarin, J.R. et al. Global needs for nitrogen fertilizer to improve wheat yield under climate change. Nat. Plants 10, 1081–1090 (2024). https://doi.org/10.1038/s41477-024-01739-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41477-024-01739-3