Abstract
Globalization intensifies the demand for agricultural products from specific regions, resulting in intensive farming practices that can exacerbate local cropland soil phosphorus (P) depletion, thereby undermining long-term food security. By integrating global data on international trade and soil-P reserves and deficits from 1970 to 2017, we demonstrate that the contribution of trade to global soil-P deficits increased from 7% in 1970 to 18% in 2017, with 54% of this impact driven by non-food consumption. Over these 48 years, developing regions exported a net of 5.8 Mt P through agricultural trade, resulting in a net increase of 13 Mt soil-P deficits. These deficits are primarily concentrated in regions with low soil-P reserves, such as sub-Saharan Africa, Latin America and Southeast Asia, thereby heightening the risks of soil-P depletion in these areas and amplifying long-term concerns about food security. This insight underscores the imperative for a broader perspective on food security—prioritizing national soil productivity rather than merely boosting the availability of food in the global market when shaping global trade policies.
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 Springer Link
- 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 data supporting the findings of this study are available within the article and its Supplementary Information files. The data used in this study are mainly collected from the online statistical databases of the FAO (FAOSTAT; http://www.fao.org/statistics/databases/en/) and the International Fertilizer Association (IFASTAT; https://www.ifastat.org/nutrient-use-efficiency). The basemaps shown in Figs. 1 and 3 and Extended Data Figs. 2, 7 and 8 are sourced from Natural Earth (http://www.naturalearthdata.com), a public domain map dataset. The global soil-P deficit database that we established is available via Figshare at https://doi.org/10.6084/m9.figshare.25782579 (ref. 66). Source data are provided with this paper.
Code availability
The environmentally extended MRIO model analysis was conducted using the Australian IELab (https://ielab.info/) infrastructure and is accessible from the authors upon request.
References
Langhans, C., Beusen, A. H. W., Mogollón, J. M. & Bouwman, A. F. Phosphorus for sustainable development goal target of doubling smallholder productivity. Nat. Sustain. 5, 57–63 (2022).
MacDonald, G. K., Bennett, E. M., Potter, P. A. & Ramankutty, N. Agronomic phosphorus imbalances across the world’s croplands. Proc. Natl Acad. Sci. USA 108, 3086–3091 (2011).
Bekunda, M. & Manzi, G. Use of the partial nutrient budget as an indicator of nutrient depletion in the highlands of southwestern Uganda. Nutr. Cycl. Agroecosys. 67, 187–195 (2003).
Stoorvogel, J. J. & Smaling, E. Assessment of Soil Nutrient Depletion in Sub-Saharan Africa: 1983–2000 Report No. 0924-3062, Vol. 1 (SC-DLO, 1990).
de Astarloa, D. D. & Pengue, W. A. Nutrients metabolism of agricultural production in Argentina: NPK input and output flows from 1961 to 2015. Ecol. Econ. 147, 74–83 (2018).
Tan, Z. X., Lal, R. & Wiebe, K. D. Global soil nutrient depletion and yield reduction. J. Sustain. Agric. 26, 123–146 (2005).
Magnone, D. et al. The impact of phosphorus on projected Sub-Saharan Africa food security futures. Nat. Commun. 13, 6471 (2022).
Meijl, H. V., Tabeau, A., Stehfest, E., Doelman, J. & Lucas, P. How food secure are the green, rocky and middle roads: food security effects in different world development paths. Environ. Res. Commun. 2, 031002 (2020).
Chen, X. et al. Physical and virtual nutrient flows in global telecoupled agricultural trade networks. Nat. Commun. 14, 2391 (2023).
Cordell, D., Drangert, J. & White, S. The story of phosphorus: global food security and food for thought. Glob. Environ. Change 19, 292–305 (2009).
Grote, U., Craswell, E. & Vlek, P. Nutrient flows in international trade: ecology and policy issues. Environ. Sci. Policy 8, 439–451 (2005).
Nesme, T., Metson, G. S. & Bennett, E. M. Global phosphorus flows through agricultural trade. Glob. Environ. Change 50, 133–141 (2018).
Bai, Z. et al. Agricultural trade impacts global phosphorus use and partial productivity. Nat. Food 4, 762–773 (2023).
Hamilton, H. A. et al. Trade and the role of non-food commodities for global eutrophication. Nat. Sustain. 1, 314–321 (2018).
Wiedmann, T. Eutrophication’s neglected drivers. Nat. Sustain. 1, 273–274 (2018).
Li, M., Wiedmann, T. & Hadjikakou, M. Towards meaningful consumption-based planetary boundary indicators: the phosphorus exceedance footprint. Glob. Environ. Change 54, 227–238 (2019).
Zw, A. et al. International trade reduces global phosphorus demand but intensifies the imbalance in local consumption. Sci. Total Environ. 830, 154484 (2022).
Barbieri, P., MacDonald, G. K., Bernard De Raymond, A. & Nesme, T. Food system resilience to phosphorus shortages on a telecoupled planet. Nat. Sustain. 5, 114–122 (2022).
Lun, F. et al. Influences of international agricultural trade on the global phosphorus cycle and its associated issues. Glob. Environ. Change 69, 102282 (2021).
Zhang, J. et al. Spatiotemporal dynamics of soil phosphorus and crop uptake in global cropland during the 20th century. Biogeosciences 14, 2055–2068 (2017).
Mogollón, J. M., Beusen, A., Van Grinsven, H., Westhoek, H. & Bouwman, A. F. Future agricultural phosphorus demand according to the shared socioeconomic pathways. Glob. Environ. Change 50, 149–163 (2018).
Sattari, S. Z., Bouwman, A. F., Giller, K. E. & van Ittersum, M. K. Residual soil phosphorus as the missing piece in the global phosphorus crisis puzzle. Proc. Natl Acad. Sci. USA 109, 6348–6353 (2012).
Lenzen, M. et al. The Global MRIO Lab—charting the world economy. Econ. Syst. Res. 29, 158–186 (2017).
Schipanski, M. E. & Bennett, E. M. The influence of agricultural trade and livestock production on the global phosphorus cycle. Ecosystems 15, 256–268 (2012).
Yang, X., Post, W. M., Thornton, P. E. & Jain, A. K. Global Gridded Soil Phosphorus Distribution Maps at 0.5-Degree Resolution (Oak Ridge National Laboratory Distributed Active Archive Center, 2014); https://doi.org/10.3334/ORNLDAAC/1223
Batjes, N. H. Global Distribution of Soil Phosphorus Retention Potential No. 2011/06 (ISRIC—World Soil Information, 2011).
Zou, T., Zhang, X. & Davidson, E. A. Global trends of cropland phosphorus use and sustainability challenges. Nature 611, 81–87 (2022).
Ten Berge, H. F. et al. Maize crop nutrient input requirements for food security in sub-Saharan Africa. Glob. Food Sec. 23, 9–21 (2019).
Mineral Commodity Summaries 2020 (U.S. Geological Survey, 2020); https://doi.org/10.3133/mcs2020
Alewell, C. et al. Global phosphorus shortage will be aggravated by soil erosion. Nat. Commun. 11, 4546 (2020).
Awika, J. M., Piironen, V. & Bean, S. (eds) Advances in Cereal Science: Implications to Food Processing and Health Promotion (American Chemical Society, 2011).
Yang, X., Post, W. M., Thornton, P. E. & Jain, A. The distribution of soil phosphorus for global biogeochemical modeling. Biogeosciences 10, 2525��2537 (2013).
Torri, S. I., Urricariet, S. & Lavado, R. S. in Soil Nutrients (ed. Miransari, M.) 277–296 (Nova Science, 2011).
Álvarez, R., Steinbach, H. S. & de Paepe, J. L. Historical balance of nitrogen, phosphorus, and sulfur of the Argentine Pampas. Cienc. Suelo 34, 231–244 (2016).
Pengue, W. A. Transgenic crops in Argentina: the ecological and social debt. Bull. Sci. Technol. Soc. 25, 314–322 (2005).
Pavinato, P. S. et al. Revealing soil legacy phosphorus to promote sustainable agriculture in Brazil. Sci. Rep. 10, 15615 (2020).
Miettinen, J., Shi, C. & Liew, S. C. Land cover distribution in the peatlands of Peninsular Malaysia, Sumatra and Borneo in 2015 with changes since 1990. Glob. Ecol. Conserv. 6, 67–78 (2016).
Kassai, P. & Tóth, G. Agricultural soil phosphorus in Hungary: high resolution mapping and assessment of socioeconomic and pedological factors of spatiotemporal variability. Sustainability 12, 5311 (2020).
Wiedmann, T. & Lenzen, M. Environmental and social footprints of international trade. Nat. Geosci. 11, 314–321 (2018).
Lenzen, M. et al. The carbon footprint of global tourism. Nat. Clim. Change 8, 522–528 (2018).
Lenzen, M. et al. International trade drives biodiversity threats in developing nations. Nature 486, 109–112 (2012).
Dalin, C., Wada, Y., Kastner, T. & Puma, M. J. Groundwater depletion embedded in international food trade. Nature 543, 700–704 (2017).
Chaves, L. S. M. et al. Global consumption and international trade in deforestation-associated commodities could influence malaria risk. Nat. Commun. 11, 1258 (2020).
Zhang, Q. et al. Transboundary health impacts of transported global air pollution and international trade. Nature 543, 705–709 (2017).
World Food Summit: Rome Declaration for World Food Security (FAO, 1996); https://www.fao.org/3/w3613e/w3613e00.htm
Bouët, A., Odjo, S. P. & Zaki, C. Africa Agriculture Trade Monitor 2020 (International Food Policy Research Institute, 2020).
Cordell, D., Brownlie, W. J. & Esham, M. Time to take responsibility on phosphorus: towards circular food systems. Glob. Environ. Change 71, 102406 (2021).
Kahiluoto, H., Pickett, K. E. & Steffen, W. Global nutrient equity for people and the planet. Nat. Food 2, 857–861 (2021).
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. Global Biogeochem. Cycles 22, GB1022 (2008).
Robinson, T. P. et al. Mapping the global distribution of livestock. PLoS ONE 9, e96084 (2014).
Franke, N. A., Boyacioglu, H. & Hoekstra, A. Y. Grey Water Footprint Accounting: Tier 1 Supporting Guidelines Value of Water Research Report Series No. 65 (UNESCO-IHE, 2013).
Lenzen, M., Moran, D., Kanemoto, K. & Geschke, A. Building Eora: a global multi-region input–output database at high country and sector resolution. Econ. Syst. Res. 25, 20–49 (2013).
Lenzen, M., Kanemoto, K., Moran, D. & Geschke, A. Mapping the structure of the world economy. Environ. Sci. Technol. 46, 8374–8381 (2012).
Heckscher, E. F., Ohlin, B. G., Flam, H. & Flanders, M. J. Heckscher-Ohlin Trade Theory (MIT Press, 1991).
Prébisch, R. The Economic Development of Latin America and Its Principal Problems (Sede de la CEPAL en Santiago (Estudios e Investigaciones), Naciones Unidas Comisión Económica para América Latina y el Caribe (CEPAL), 1950).
Bibi, S. Prebisch and the terms of trade. Resour. Policy 90, 104813 (2024).
Duarte, R., Pinilla, V. & Serrano, A. Long term drivers of global virtual water trade: a trade gravity approach for 1965–2010. Ecol. Econ. 156, 318–326 (2019).
Frankel, J. A. & Rose, A. K. Is trade good or bad for the environment? Sorting out the causality. Rev. Econ. Stat. 87, 85–91 (2002).
Jelinek, M. & Porter, M. E. The competitive advantage of nations. Adm. Sci. Q. 37, 507 (1992).
Kaufmann, D., Kraay, A. & Mastruzzi, M. The worldwide governance indicators: methodology and analytical issues. Hague J. Rule Law 3, 220–246 (2011).
Lenzen, M. Aggregation versus disaggregation in input–output analysis of the environment. Econ. Syst. Res. 23, 73–89 (2011).
Bullard, C. W. & Sebald, A. V. Monte Carlo sensitivity analysis of input–output models. Rev. Econ. Stat. 70, 708–712 (1988).
Bullard, C. W. & Sebald, A. V. Effects of parametric uncertainty and technological change on input–output models. Rev. Econ. Stat. 59, 75–81 (1977).
Hondo, H., Sakai, S. & Tanno, S. Sensitivity analysis of total CO2 emission intensities estimated using an input–output table. Appl. Energy 72, 689–704 (2002).
Lenzen, M., Wood, R. & Wiedmann, T. Uncertainty analysis for multi-region input–output models—a case study of the UK’s carbon footprint. Econ. Syst. Res. 22, 43–63 (2010).
Niu, K., Li, M., Han, X. & Zhang, J. Global soil-P deficits database from 1970 to 2017. Figshare https://doi.org/10.6084/m9.figshare.25782579.v1 (2024).
Acknowledgements
This paper is supported by the National Natural Science Foundation of China (grant nos. 72104239 (K.N.), 42325707 (B.G.) and 42261144001 (B.G.)); the Agricultural Science and Technology Innovation Program (grant nos. 10-IAED-ZK-04-2024 (K.N.), 10-IAED-06-2024 (K.N.) and 10-IAED-SZD-01-2024 (K.N.)); the Australian Research Council through its Discovery Project grants (nos. DP0985522 (M. Lenzen), DP130101293 (M. Lenzen), DP190102277 (M. Lenzen), DP200103005 (A.M.) and DP200102585 (M. Lenzen and M. Li)), its Linkage Project (grant no. LP200100311 (A.M.)) and its Large Infrastructure, Equipment and Facilities grant (no. LE160100066 (T.W. and M. Lenzen)); the Australian Research Council Research Hub (grant no. IH190100009 (A.M.)); the Discovery Early Career Researcher Award (A.M.); and the National eResearch Collaboration Tools and Resources project (NeCTAR, T.W. and M. Lenzen) through its Industrial Ecology Virtual Laboratory VR201. NeCTAR projects are Australian government projects conducted as part of the Super Science initiative and financed by the Education Investment Fund. M. Lenzen acknowledges financial support from the Hanse-Wissenschaftskolleg in Delmenhorst, Germany, through its HWK Fellowships. M. Li acknowledges support from the University of Sydney Horizon Fellowship. A.M. acknowledges support from the University of Sydney SOAR Prize.
Author information
Authors and Affiliations
Contributions
K.N. and M. Li conceptualized and planned the paper. M. Li, M. Lenzen and A.M. constructed the MRIO matrix. K.N., M. Li and X.H. conducted the data collection, processing and chart creation. K.N., B.G., A.M., X.H., M. Li, M. Lenzen, T.W. and S.J. wrote the paper. All authors contributed to revising the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Sustainability thanks Daniel Magnone 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 Flow chart for this study.
We illustrate the process of calculating the soil-P budget at the production site (a) and constructing the soil-P deficits (SPD) database (b). The SPD is then integrated with the MRIO trade network (c) to assess the impact of food and non-food product consumption at the consumption site on the SPD of the production site (d). Subsequently, we combine this analysis with information on accumulated soil-P, natural soil-P level, soil properties and affected crops at the production site to determine the risk of soil-P depletion (e). The global agricultural trade network (ATN) centers around the agricultural supply chain indicated by the green arrows in the diagram. In contrast, the multi-regional input-output (MRIO) trade network employed in our analysis covers not only the impact of the agricultural supply chain (green arrows) but also that of the non-food supply chain, as illustrated by the red arrows. The cropland in the figure was drawn using Adobe illustrator.
Extended Data Fig. 2 Spatially explicit maps illustrating the soil-P budgets for various crop types in 1970 and 2017 across 44 nations exporting soil-P deficits (kg P/ha/yr).
a1-9 and b1-9 refer to the year 1970 and 2017, respectively. c1-c9 are the changes between 1970 and 2017, where blue signals rising soil-P deficits, and red signals declining soil-P deficits or increasing soil-P surplus. Basemaps from Natural Earth (https://www.naturalearthdata.com/).
Extended Data Fig. 3 Comparison of per capita consumption-driven soil-P deficits for 10 aggregated regions.
The size of the circle represents a region’s population. See Supplementary Table 7 for full names and abbreviations of aggregated regions.
Extended Data Fig. 4 Trade balance of soil-P deficits by region and sector over 48 years for top 20 net-importing and net-exporting countries/regions.
The green bars to the left of X = 0 signify the rise in soil-P deficits within a particular region’s agricultural system due to trade, whereas the varied bars on the right side of X = 0 indicate an escalation in soil-P deficits in foreign regions due to the consumption of products from a specific sector within that region. See Supplementary Table 6 for full names and abbreviations of 90 countries and regions. ‘Agri products’ refer to Agricultural products. ‘Manuf products’ refer to manufacturing products. ‘Acc & FS’ refer to Accommodation and food service. ‘Fin, Ret, & Whol’ refer to Finance, retail, and wholesale.
Extended Data Fig. 5 Modeling diagram of driving factors of net export of cropland soil-P deficits in each region.
The dashed downward curve in the figure represents the negative effect of population on the net export of soil-P deficits from a region, but it does not pass the statistical significance test.
Extended Data Fig. 6 Trade balance of soil-P deficits by region over 48 years.
A value less than 0 signals a rise in local soil-P deficits through trade, whereas a value greater than 0 indicates that imports of this region exacerbate the soil-P deficits in other foreign regions. The black bars represent low-income food-deficit nations recognized by the FAO, and it can be seen that most of them are net exporters of soil-P deficits. See Supplementary Table 6 for full names and abbreviations of 90 countries and regions. Unit, Mt.
Extended Data Fig. 7 Net-export density of soil-P deficits per unit area over 48 years.
a, countries with low level of food security; b, developing countries with large agricultural exports. c, sub-Saharan Africa’s net imports of P through food trade and soil-P deficits occurred due to trade from 1970 to 2017. Basemaps in a and b from Natural Earth (https://www.naturalearthdata.com/).
Extended Data Fig. 8 Aggregating 90 regions into 10 regions.
By considering information on economic development and geographic location, we aggregated 90 regions into 10 broad regions. See Supplementary Table 7 for full names and abbreviations of aggregated regions. Basemaps from Natural Earth (https://www.naturalearthdata.com/).
Supplementary information
Supplementary Information
Supplementary Tables 1–15 and text.
Supplementary Video 1
Top ten flows of soil-P deficits (in kt yr−1) embodied in traded commodities.
Source data
Source Data Figs. 1–5
Statistical source data.
Source Data Extended Data Figs. 2–4, 6 and 7
Statistical source data.
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
Niu, K., Li, M., Lenzen, M. et al. Impacts of global trade on cropland soil-phosphorus depletion and food security. Nat Sustain (2024). https://doi.org/10.1038/s41893-024-01385-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41893-024-01385-9