Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Global crop production increase by soil organic carbon

Abstract

Soil organic carbon sequestration has been promoted to combat climate change while improving soil fertility. However, its quantitative contribution to crop productivity has proven elusive. Using data from 13,662 controlled field trials with 66,593 treatments across a broad range of soils, climates and management practices, we here show that yields increase with increased soil organic carbon, until no further increase (p < 0.05) occurs above mean optimum soil organic carbon of 43.2–43.9 g kg−1 for maize, 12.7–13.4 g kg−1 for wheat and 31.2–32.4 g kg−1 for rice. Sequestering soil organic carbon is one-fifth as effective (that is, 80% less) as nitrogen fertilization for improving crop yield where soil management is optimized. By increasing soil organic carbon beyond current technology to optimum levels, global production of the three most important staple crops increases by 4.3% (sufficient to provide calories for 640 million people). However, currently available management practices would increase crop production by only 0.7% once other production constraints have already been addressed. Therefore, yield improvements under currently available technologies are unlikely to drive adoption of soil organic carbon sequestration globally, except in hot-spot regions where crop production benefits most, or unless novel practices that allow greater soil organic carbon sequestration beyond current limitations can further increase yields cost-effectively.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Causal effect of SOC on crop yield and SOC yield gap.
Fig. 2: Drivers of crop yield.
Fig. 3: Drivers of SOC yield gap.
Fig. 4: Global potential to increase yields by increasing SOC.

Similar content being viewed by others

Data availability

Data that support these findings are available via figshare (https://doi.org/10.6084/m9.figshare.24137280).

Code availability

Codes for processing the data are available via figshare (https://doi.org/10.6084/m9.figshare.24137280).

References

  1. IPCC Climate Change 2021: The Physical Science Basis (eds Masson-Delmotte, V. et al.) (Cambridge Univ. Press, 2021).

  2. IPCC Special Report on Global Warming of 1.5°C (eds Masson-Delmotte, V. et al.) (WMO, 2018).

  3. Anderson, C. M. et al. Natural climate solutions are not enough. Science 363, 933–934 (2019).

    Article  Google Scholar 

  4. Field, C. B. & Mach, K. J. Rightsizing carbon dioxide removal. Science 356, 706–707 (2017).

    Article  Google Scholar 

  5. Batjes, N. H. Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sci. 47, 151–163 (1996).

    Article  Google Scholar 

  6. Minasny, B. et al. Soil carbon 4 per mille. Geoderma 292, 59–86 (2017).

    Article  Google Scholar 

  7. Amelung, W. et al. Towards a global-scale soil climate mitigation strategy. Nat. Commun. 11, 5427 (2020).

    Article  Google Scholar 

  8. Oldfield, E. E., Bradford, M. A. & Wood, S. A. Global meta-analysis of the relationship between soil organic matter and crop yields. Soil 5, 15–32 (2019).

    Article  Google Scholar 

  9. Lal, R. Soil carbon sequestration impacts on global climate change and food security. Science 304, 1623–1627 (2004).

    Article  Google Scholar 

  10. Lal, R. Soil organic matter content and crop yield. J. Soil Water Conserv. 75, 27A–32A (2020).

    Article  Google Scholar 

  11. Lal, R. Enhancing crop yields in the developing countries through restoration of the soil organic carbon pool in agricultural lands. Land Degrad. Dev. 17, 197–209 (2006).

    Article  Google Scholar 

  12. Pan, G., Smith, P. & Pan, W. The role of soil organic matter in maintaining the productivity and yield stability of cereals in China. Agric. Ecosyst. Environ. 129, 344–348 (2009).

    Article  Google Scholar 

  13. Yu, X. et al. Global meta-analysis of nitrogen fertilizer use efficiency in rice, wheat and maize. Agric. Ecosyst. Environ. 338, 108089 (2022).

    Article  Google Scholar 

  14. Oelofse, M. et al. Do soil organic carbon levels affect potential yields and nitrogen use efficiency? An analysis of winter wheat and spring barley field trials. Eur. J. Agron. 66, 62–73 (2015).

    Article  Google Scholar 

  15. Hudson, B. D. Soil organic matter and available water capacity. J. Soil Water Conserv. 49, 189–194 (1994).

    Google Scholar 

  16. Murphy, B. Impact of soil organic matter on soil properties—a review with emphasis on Australian soils. Soil Res. 53, 605–635 (2015).

    Article  Google Scholar 

  17. Kane, D. A., Bradford, M. A., Fuller, E., Oldfield, E. E. & Wood, S. A. Soil organic matter protects US maize yields and lowers crop insurance payouts under drought. Environ. Res. Lett. 16, 044018 (2021).

    Article  Google Scholar 

  18. Minasny, B. & McBratney, A. Limited effect of organic matter on soil available water capacity. Eur. J. Soil Sci. 69, 39–47 (2018).

    Article  Google Scholar 

  19. Paustian, K. et al. Climate-smart soils. Nature 532, 49–57 (2016).

    Article  Google Scholar 

  20. DeLisi, C. et al. The role of synthetic biology in atmospheric greenhouse gas reduction: prospects and challenges. Biodes. Res. 2020, 1016207 (2020).

    Article  Google Scholar 

  21. Ogle, S. M. et al. in 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories: Agriculture, Forestry and Other Land Use, Vol. 4 (eds Calvo Buendia, E. et al.) Ch. 2 (IPCC, 2019).

  22. Kimetu, J. M. et al. Reversibility of soil productivity decline with organic matter of differing quality along a degradation gradient. Ecosyst 11, 726–739 (2007).

    Article  Google Scholar 

  23. Oldfield, E. E. et al. Positive associations of soil organic matter and crop yields across a regional network of working farms. Soil Sci. Soc. Am. J. 86, 384–397 (2022).

    Article  Google Scholar 

  24. Lehmann, J. et al. Biochar in climate change mitigation. Nat. Geosci. 14, 883–892 (2021).

    Article  Google Scholar 

  25. South, P. F., Cavanagh, A. P., Lopez‐Calcagno, P. E., Raines, C. A. & Ort, D. R. Optimizing photorespiration for improved crop productivity. J. Integr. Plant Biol. 60, 1217–1230 (2018).

    Article  Google Scholar 

  26. Kroeger, M. E. et al. Microbial community composition controls carbon flux across litter types in early phase of litter decomposition. Environ. Microbiol. 23, 6676–6693 (2021).

    Article  Google Scholar 

  27. Qiao, L. et al. Soil quality both increases crop production and improves resilience to climate change. Nat. Clim. Change 12, 574–580 (2022).

    Article  Google Scholar 

  28. Zhang, F. et al. Integrated nutrient management for food security and environmental quality in China. Adv. Agron. 116, 1–40 (2012).

    Article  Google Scholar 

  29. Yang, H. S. et al. Hybrid-maize—a maize simulation model that combines two crop modeling approaches. Field Crops Res. 87, 131–154 (2004).

    Article  Google Scholar 

  30. Harris, I. C., Jones, P. D. & Osborn, T. CRU TS4.04: Climatic Research Unit (CRU) Time-Series (TS) Version 4.04 of High-Resolution Gridded Data of Month-by-Month Variation in Climate (Jan. 1901–Dec. 2019) (Univ. of East Anglia Climatic Research Unit, 2020); https://catalogue.ceda.ac.uk/uuid/89e1e34ec3554dc98594a5732622bce9

  31. Jiang, R. F. & Cui, J. Y. in Fertilization Technology Highlights (Ed. Zhang, F. S.) Ch. 5 (China Agricultural Univ. Press, 2006).

  32. National Soil Survey Office Chinese Soil (China Agricultural Press, 1998).

  33. Poggio, L. et al. SoilGrids 2.0: producing soil information for the globe with quantified spatial uncertainty. Soil 7, 217–240 (2021).

    Article  Google Scholar 

  34. Dobarco, M. R., Cousin, I., Le Bas, C. & Martin, M. P. Pedotransfer functions for predicting available water capacity in French soils, their applicability domain and associated uncertainty. Geoderma 336, 81–95 (2019).

    Article  Google Scholar 

  35. Breiman, L. Random forests. Mach. Learn. 45, 5–32 (2001).

    Article  Google Scholar 

  36. Altmann, A., Toloşi, L., Sander, O. & Lengauer, T. Permutation importance: a corrected feature importance measure. Bioinform 26, 1340–1347 (2010).

    Article  Google Scholar 

  37. Greenwell, B. M. pdp: an R package for constructing partial dependence plots. R J. 9, 421–436 (2017).

    Article  Google Scholar 

  38. Eilers, P. H. A perfect smoother. Anal. Chem. 75, 3631–3636 (2003).

    Article  Google Scholar 

  39. Altey, S. & Wager, S. Estimating treatment effects with causal forests: an application. Obs. Stud. 5, 37–51 (2019).

    Google Scholar 

  40. Reddy, S. et al. in 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories: Agriculture, Forestry and Other Land Use, Vol. 4 (eds Calvo Buendia, E. et al.) Ch. 3 (IPCC, 2019).

  41. Monfreda, C., Ramankutty, N. & Foley, J. A. Farming the planet: 2. geographic distribution of crop areas, yields, physiological types, and net primary production in the year 2000. Glob. Biogeochem. Cycles 22, GB1022 (2008).

    Article  Google Scholar 

  42. Roe, S. et al. Land-based measures to mitigate climate change: potential and feasibility by country. Glob. Change Biol. 27, 6025–6058 (2021).

    Article  Google Scholar 

  43. FAOSTAT: Crops and Livestock Products (FAO, 2021); https://www.fao.org/faostat/en/#data/QCL

  44. Alexandratos, N. & Mernies, J. in World Agriculture: Towards 2015/2030 An FAO Perspective: Prospects for Food and Nutrition (Ed. Bruinsma, J.) Ch. 2 (FAO, 2003).

  45. Nuss, E. T. & Tanumihardjo, S. A. Maize: a paramount staple crop in the context of global nutrition. Compr. Rev. Food Sci. Food Saf. 9, 417–436 (2010).

    Article  Google Scholar 

  46. Rome Declaration on World Food Security and World Food Summit Plan of Action (FAO, 1996).

  47. FAOSTAT: Suite of Food Security Indicators (FAO, 2022); https://www.fao.org/faostat/en/#data/FS

  48. Department of Economic and Social Affairs World Population Prospects (United Nations, 2019); https://population.un.org/wpp

  49. The State of Food Insecurity in the World (FAO, 2000).

  50. Global Food Security Index: Exploring Challenges and Developing Solutions (Economist Impact, 2021); https://impact.economist.com/sustainability/project/food-security-index/

Download references

Acknowledgements

We thank the National Key Research and Development Program of China (2017YFD0200108) and the National Natural Science Foundation of China (31972520) for providing financial support to Y.M., M.F. and L.Q. We thank the Cornell Institute for Digital Agriculture (CIDA RIF2019) for supporting D.W. and AI-CLIMATE (NIFA-2023-67021-39829) for supporting J.L. on this project. We thank the China Scholarship Council for providing funds to Y.M. to pursue his study at Cornell University. We thank L. M. Johnson from the Cornell Statistical Consulting Unit for supporting the data analysis work. We thank all agricultural scientists and extension personnel who provided local technical assistance and management during the project.

Author information

Authors and Affiliations

Authors

Contributions

Y.M., D.W., M.F. and J.L. designed the research; Y.M., D.W., M.F., L.Q., R.L. and J.L. collected field and other research data; Y.M., D.W., M.F. and J.L. contributed to methodology and data analysis; Y.M. and D.W. contributed to visualization; D.W., M.F. and J.L. supervised the work; Y.M. and D.W. wrote the first manuscript; Y.M., D.W., M.F. and J.L. edited the manuscript; all authors read and approved the final manuscript.

Corresponding authors

Correspondence to Mingsheng Fan or Johannes Lehmann.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Geoscience thanks Marcel van der Heijden and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editor: Xujia Jiang, in collaboration with the Nature Geoscience team.

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 Distribution of field trials for maize, wheat and rice.

ac, Maize (a), wheat (b) and rice (c). The round symbols on the three top panels represent the locations of field trials. The bottom panel indicates the eight IPCC climate zones (large fractions of climate zones where the three studied crops distribute are covered by field trials). Numbers in brackets under the climate name indicate the number of field trials for each climate zone of each crop. Numbers in brackets under the crop name indicate the total number of field trials for each crop. Basemap reproduced from Resource and Environment Science and Data Center, Institute of Geographic Sciences and Natural Resources, Chinese Academy of Sciences.

Extended Data Fig. 2 Distribution of climate-and-soil specific optimum SOC for maize, wheat and rice.

ac, Maize (a), wheat (b) and rice (c). The optimum SOC is the value above which no additional increase in crop yield can be achieved by increasing SOC. Basemap reproduced from ArcGIS/Esri.

Extended Data Fig. 3 Causal effect of soil organic carbon (SOC) on crop yield under different nitrogen input rates.

a, three crops combined. b, maize. c, wheat, d, rice. The solid lines represent mean partial dependence of crop yield on SOC, and the ribbons represent 95% confidence intervals. N0, N0.5, N1.0 and N1.5 indicate 0, 50%, 100% and 150% of optimum nitrogen input rates, respectively, with optimum phosphorus and potassium input rates. Control means no fertilizer inputs. For three crops combined, predicted yield was normalized as the ratio of on-farm measured yield in each trial to predicted minimum SOC yield gap for each crop under the same treatment. Predicted maximum SOC yield gap is the difference between predicted yield under optimum SOC and predicted minimum yield.

Extended Data Fig. 4 Causal effect of soil organic carbon (SOC) on crop yield across climate zones for maize, wheat and rice.

a-d, cool dry (a), cool moist (b), warm dry (c), and warm moist (d). The solid lines represent the mean partial dependence of normalized crop yield on SOC, and the ribbons represent 95% confidence intervals. Optimum SOC is the SOC level beyond which no significant (p < 0.05) yield increases were observed. The significance of the yield increase was estimated by two-sided t-tests. Normalized predicted yield was calculated as the ratio of predicted yield from partial dependence to predicted minimum mean yield for each crop under each climate zone. Normalized predicted maximum SOC yield gap is the difference between normalized predicted yield under optimum SOC and normalized predicted minimum yield.

Extended Data Fig. 5 Distribution of yield increases derived from soil organic carbon (SOC) for three studied crops.

Regions on the maps are derived from Fig. 4 and classified into quartiles. a, c, e, technical potential represents yield gains when SOC is increased to current technical potential SOC levels for maize (a), wheat (c), and rice (e). Current technical potential SOC levels can be achieved using a combination of high organic matter inputs and zero tillage; b, d, f, unconstrained represents yield gains when SOC increased to optimum levels for (b), wheat (d), and rice (f). Basemap reproduced from ArcGIS/Esri.

Extended Data Fig. 6 Relative importance (%) of variables for the global extrapolation of yield models and its predicted standard errors for each model.

a,b,c, bar charts indicate relative importance (permutation method; see Methods) of variables as predictors of random forest extrapolation yield models. d,e,f, line charts represent the predictions of standard errors for maize (d), wheat (e), and rice (f) estimated by the jackknife method (See methods; Uncertainty analysis of global extrapolations). NPK fertilization = nitrogen, phosphorus and potassium fertilizer rate; PET = potential evapotranspiration; PRE, Tmax and Tmin = precipitation, maximum mean temperature and minimum mean temperature; AWC = available water capacity.

Extended Data Fig. 7 Correspondence of distributions of climate zones where the majority of field trials are located and producing regions of the three studied crops.

a, maize, b, wheat, c, rice. Climate zones were divided based on IPCC criteria43. Ntreat is the number of treatments rather than the number of field trials in each climate zones. The distributions of maize, wheat and rice producing regions were obtained from EarthStat44. The shaded regions are the global producing regions for the three studied crops. Climate zones where only few field trials distribute are excluded from the map including maize tropical dry zone (Ntreat = 10), rice tropical dry zone (Ntreat = 15) and rice cool temperate dry zone (Ntreat = 15). Basemap reproduced from ArcGIS/Esri.

Extended Data Fig. 8 Global yield increases of maize, wheat and rice by increasing soil organic carbon (SOC) in each country.

Yield increases of the three crops by increasing SOC in each country was normalized to total yield increases of the three crops from SOC relative to annual current crop production of the three crops (a, b), to the yield increases of the three crops from SOC per unit harvest area (c, d), and to the total yield increases of the three crops from SOC relative to annual crop production of the three crops at best management (e, f). Total yield increase relative to annual current crop production was calculated by using yield increases of the three crops by increasing SOC in each country under technical potential and unconstrained increase in SOC divided by its current annual crop production of the three crops in each country. Yield increases of the three crops per unit harvested area was calculated by using the yield increase of the three crops by increasing SOC in each country under the technical potential and unconstrained SOC increase divided by the total harvested area of the three crops in each country. Total yield increases relative to annual crop production at best management was calculated by using yield increases of the three crops by increasing SOC in each country under technical potential and unconstrained increase in SOC divided by annual crop production of the three crops at best management in each country.

Extended Data Fig. 9 Depth of food deficit, prevalence of undernourishment and the proportions of food deficit addressed by improving SOC under the technical potential or unconstrained SOC increases.

Depth of food deficit (a) is the average calorie deficiency per capita per day below the minimum dietary daily energy requirement in each country. Prevalence of undernourishment (b) is the ratio of total undernourished population to the total population in each country that experiences national food deficits. Proportions of food deficit addressed were estimated by using the kilocalories provided by a daily increase of maize, wheat and rice per undernourished person resulting from a SOC increase under either the technical potential (c,e) or unconstrained SOC increases (d,f) for each country divided by the daily food deficit in kilocalories per undernourished person (Supplementary Table 6). The countries (c-f) were grouped based on the Global Food Security Index50 and classified as ‘weak and moderate food security environment’ or ‘good and very good food security environment’ (Supplementary Table 6). The technical potential represented yield gains when SOC increased to levels that are achieved using a combination of high organic matter inputs and zero tillage; an unconstrained increase in SOC represented yield gains when SOC increased to optimum levels.

Supplementary information

Supplementary Information

Supplementary Text 1–4, Figs. 1 and 2 and Tables 1–6.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ma, Y., Woolf, D., Fan, M. et al. Global crop production increase by soil organic carbon. Nat. Geosci. 16, 1159–1165 (2023). https://doi.org/10.1038/s41561-023-01302-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41561-023-01302-3

This article is cited by

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing