Abstract
Sustainable food systems require malnutrition and climate change to be addressed in parallel. Here, we estimate the non-CO2 greenhouse gas emissions resulting from closing the world’s dietary nutrient gap—that between country-level nutrient supply and population requirements—for energy, protein, iron, zinc, vitamin A, vitamin B12 and folate under five climate-friendly intervention scenarios in 2030. We show that improving crop and livestock productivity and halving food loss and waste can close the nutrient gap with up to 42% lower emissions (3.03 Gt CO2eq yr−1) compared with business-as-usual supply patterns with a persistent nutrient gap (5.48 Gt CO2eq yr−1). Increased production and trade of vegetables, eggs, and roots and tubers can close the nutrient gap with the lowest emissions in most countries—with ≤23% increase in total caloric production required for 2030 relative to 2015. We conclude that the world’s nutrient gap could be closed without exceeding global climate targets and without drastic changes to national food baskets.
Main
The global syndemic—synchronous pandemics of malnutrition and climate change—poses a growing threat to humanity1, with the COVID-19 pandemic exacerbating these effects2. At the same time, food systems are responsible for one-third (14–22 Gt CO2eq yr−1 in 2015) of global greenhouse gas (GHG) emissions, ~33% of which are direct non-CO2 emissions (that is, CH4 and N2O) occurring on-farm3. The critical role of food systems in limiting mean temperature increase to 1.5 °C is now well established45, while nations have pledged to end all forms of malnutrition by 20306.
Despite a doubling of food production in caloric terms between 1995 and 2015, more than 40% of the global population continues to live in countries with inadequate micronutrient (for example, vitamins and minerals) supplies to meet population-level physiological requirements as a result of current food baskets that are largely dominated by cereals78. The shortfall between dietary nutrient (micro- and macronutrients) requirements and supply, at the country level, is termed the nutrient gap910, which implies that adequate nutrition is not possible even with equal distribution within countries. Those regions with higher nutrient gaps, such as sub-Saharan Africa and Southern Asia, tend to have much larger GHG emissions intensity per kilogram of animal protein due to low productivity11. They are also expected to have the highest population growth12 and may experience insufficient vegetable and fruit supply13. Hence, ensuring adequate nutrient supplies without exacerbating global warming requires carefully designed policies informed by appropriate indicators141516.
Recent assessments have largely focussed on income-driven demand1718 and wholesale dietary shifts (for example, towards flexitarian or vegetarian diets)1920 rather than physiological requirements and country-specific nutrient gaps. Some work has incorporated environmental boundaries in pursuit of optimal diets that provide recommended amounts of protein21, fat10 and other nutrients22. Production-based studies have incorporated composite productivity indicators linking nutrients to land1423 and water24. However, they have often been limited to specific regions and/or products. The ‘nutritional life cycle assessment’ approach has been applied to compare regional differences in environmental impacts of nutrient production, emphasizing the importance of a nutrition angle for better-informed comparisons25. Analyses of nutrients and emissions embedded in household food waste have suggested that global food waste is equivalent to 15% of recommended energy and vitamin A intake, and 6.6% of the food-related non-CO2 GHG limit to keep global warming below 2 °C26. Overall, there is now an urgent need to identify nutrient requirements lacking in national food supplies and to close these gaps with the lowest emissions27.
Here, we provide detailed estimates of the non-CO2 GHG emissions associated with closing the nutrient gap to address two dimensions of malnutrition, namely undernutrition and micronutrient deficiencies, under five climate-friendly intervention scenarios. We developed a composite indicator of emissions intensity of nutrient production to estimate non-CO2 emissions associated with closing energy, protein, iron, zinc, vitamin A, vitamin B12 and folate gaps (that is, meeting population-adjusted nutrient requirements) of populations in 2030, following Sustainable Development Goal (SDG) 2.2—‘By 2030, end all forms of malnutrition’. We used linear programming to optimize additional food production or trade to minimize the emissions from closing the nutrient gap. Given our focus on primary agricultural production and farm-level emissions, we did not consider the contribution of fortification to nutrient supplies. We optimized food supplies based on current highly disaggregated food baskets7 that reflect national food preferences to avoid fundamental changes in diets. Our optimization covers 156 crop and 40 (terrestrial) animal products for 128 countries. Four of the five climate-friendly scenarios involve increasing domestic production coupled with halving loss and waste, and improved crop and livestock productivity, while one assumes climate-friendly international food trade (Table 1).
Table 1 Scenario codes and brief descriptionFull size table
With a climate-friendly and nutrition-sensitive approach, we optimized food supply patterns to minimize the additional emissions while closing the nutrient gap (that is, all nutrient gaps) using linear programming. We used a range of data sources to calculate product-specific emissions. The scope of our analysis is limited to direct non-CO2 (CH4 and N2O) agricultural emissions following the classification by the Intergovernmental Panel on Climate Change (IPCC)28 due to a lack of data on product-specific emissions from land use and land-use change. Among direct agricultural emissions, we also excluded emissions from sewage applied to soils, liming and urea application, which do not have crop-specific attribution and represent less than 1% of agricultural emissions globally3. Our focus on non-CO2 emissions is in line with the existing literature on climate taxes29 and carbon footprints3031, as CO2 emissions from up- and downstream activities (for example, energy use) are allocated to other sectors including energy, processing and transport as per the IPCC classification.
We first assessed scenarios based on domestic production, which included productivity and food loss and waste interventions (Table 1). Next, a trade scenario was introduced to explore the potential of exploiting comparative advantage in minimizing the emissions while closing the nutrient gap. We presented our results based on default emissions factors from the IPCC Tier 1 method28, supplemented by lower and upper quantiles in parentheses (see Uncertainty estimates). To interpret our findings in the context of global climate targets, we compared them against what we call allowable food production emissions. The term refers to the non-CO2 agriculture, forestry and other land use (AFOLU) emissions in 2030 compatible with the Paris Agreement32 (see Paris Agreement and allowable food production emissions), which we downscaled according to the scope of this study (for example, population and emission sources).
Results
Emissions associated with closing the nutrient gap
Juxtaposition of country-level nutrient gaps and agricultural GHG emissions revealed that countries with large nutrient gaps, mainly concentrated in sub-Saharan Africa and Southern Asia, also tended to have a high emissions intensity of production (Fig. 1). Based on the United Nations medium variant population estimates for 2030, we estimated that total nutrient requirements would increase by ~21% for energy and ~29% for protein, vitamin A, vitamin B12, and folate (Source Data Fig. 1). In contrast, iron and zinc requirements would decrease by 36% due to decreasing birth rates and associated reductions in pregnancy in countries with large populations such as China, India and Indonesia.
Fig. 1: Global maps of dietary nutrient gaps under current food loss and waste rates in 2030 compared with emissions intensity of total nutrient production in 2015.
Countries are grouped into quartiles and coloured accordingly for ease of comparison. For example, light purple (‘Very low’) represents countries in the lowest nutrient gap quartile (see Source Data Fig. 1), while dark red (‘Very high’) represents countries in the highest quartile for emissions intensity of total production for a given nutrient. Higher emissions intensity of vitamin B12 production from animal sources, such as cow milk (see Source Data Fig. 1), particularly suggests low productivity in livestock production. Dominance of ruminant meat and dairy leads to high/very high emissions intensity of total nutrient production in countries such as Australia and Brazil, despite high livestock productivity in these countries, because livestock-related emissions represent the bulk of agricultural emissions and their dominance in national food baskets determines total emission volumes. Maps were drawn using the tmap R package79.
Total non-CO2 emissions from agricultural activities in 128 countries covering 89% of the global population reached 4.62 (4.27–6.26) Gt CO2eq yr−1 in 2013–2015 (Fig. 2). This represents ~77% of total AFOLU emissions, including CO2 emissions from drained organic soils, forestland and net forest conversion, in 229 countries. Assuming constant production-based emissions per capita and nutrient adequacy (that is, nutrient supply/population-level requirements) into the future (business-as-usual scenario; BaU) would result in 5.48 (4.76–7.02) Gt CO2eq yr−1 in 2030, exceeding 75th percentile of the allowable emission estimates compatible with the Paris Agreement. However, between-country differences in emissions intensity of nutrient production were as high as 200-fold for ruminant products such as cow’s milk (Source Data Fig. 1). Such heterogeneity in emissions intensity of production determined the effectiveness of productivity and trade scenarios.
Fig. 2: Total emissions from closing the nutrient gap against the allowable non-CO2 emissions for food production in 2030 compatible with the Paris Agreement.
Bars show the total emissions for 128 countries, that is, n = 128, with the top of the bar corresponding to the mode as the default measure of central tendency. Results are provided across five scenarios (Table 1) and current supply patterns extrapolated to 2030 populations (BaU). Hence, the nutrient gap persists under BaU. Error bars show the 25th and 75th percentiles and are negatively skewed (see Uncertainty estimates). The orange-shaded area represents the spread of allowable non-CO2 emissions for food production in 2030 (25th percentile (Q1): 4.33 Gt CO2eq yr−1 and 75th percentile (Q3): 5.31 Gt CO2eq yr−1). The solid red line represents the median (Q2) (4.67 Gt CO2eq yr−1). Allowable emissions are calculated based on an ensemble of models3277 as described in section Paris Agreement and allowable food production emissions. The values shown here are scaled based on the scope of emission sources and total population covered in this study. See Source Data Fig. 2 for full results.
Climate-friendly and nutritionally targeted increases in production closed the nutrient gap with lower emissions compared with the BaU for 2030 (Fig. 2). Under current productivity and loss and waste patterns (D-CP-FLW), emissions decreased by 11%, compared with the BaU, to 4.89 (4.52–6.70) Gt CO2eq yr−1 (Fig. 2). Closing the nutrient gap under the half loss and waste scenario (D-CP-HLW) resulted in a 22% reduction (compared with BaU) in emissions, with 4.28 (3.95–5.83) Gt CO2eq yr−1. Closing the crop yield gap increased baseline crop emissions by 6% due to increased fertilizer use, while enhancing livestock productivity decreased baseline livestock emissions by 28%. Overall, improving agricultural productivity (D-IP-FLW) reduced the emissions associated with closing the nutrient gap by 33% to 3.65 (3.35–5.00) Gt CO2eq yr−1.
When we combined half loss and waste with improved productivity (D-IP-HLW), the reduction in emissions was up to 42%, with 3.19 (2.93–4.34) Gt CO2eq yr−1 for 2030. Finally, closing the nutrient gap through increased imports of climate-friendly products (T-CP-FLW) showed a 14% decrease in the emissions, resulting in 4.70 (4.34–6.36) Gt CO2eq yr−1. Our findings under the domestic production (D-CP-FLW) and trade (T-CP-FLW) scenarios were similar because optimal food products were mostly plant-based sources, which had much smaller between-country differences in emissions intensity compared with livestock products. Overall, when compared against the allowable non-CO2 emissions for food production by 2030, closing the nutrient gap through improved productivity and half loss and waste scenarios (D-CP-HLW, D-IP-FLW and D-IP-HLW) helped keep median food system emissions below the 25th percentile of the allowable emissions compatible with the Paris Agreement.
Relative performance of climate-friendly scenarios varied slightly by national income level (Fig. 3). Owing to substantial differences in the nutrient gap, low- and lower-middle-income groups required higher production increases and together accounted for more than 80% of additional emissions. Halving food loss and waste (D-CP-HLW) and improved productivity (D-IP-FLW) mitigated emissions to a larger extent in the low-income group compared with other income groups. With current productivity and loss and waste patterns, increasing imports (T-CP-FLW) showed 13% and 6% lower emissions compared with increasing domestic production (D-CP-FLW) in the low- and lower-middle-income groups, respectively. On the other hand, emissions were similar under domestic production and trade scenarios in high- and upper-middle-income groups, which can be explained by the observation that they are the major trade partners exporting to low- and lower-middle-income groups.
Fig. 3: Breakdown of emissions by individual GHGs arising from closing the nutrient gap.
Bars show the total emissions for 128 countries based on default emissions factors, with the top of the bar corresponding to the mode as the default measure of central tendency. Results are provided by income level and across five climate-friendly scenarios and current supply patterns extrapolated to 2030 populations (BaU). Higher shares of livestock products in their food supply baskets result in larger contributions by CH4 to total GHG emission in the upper-middle- and high-income countries. Error bars show 25th and 75th percentiles for total GHG emissions. See Source Data Fig. 2 for full results. See Extended Data Fig. 1a for CH4 emissions (Mt CH4 yr−1) and Extended Data Fig. 1b for N2O emissions (Mt N2O yr−1).
Climate-friendly and nutrition-sensitive supply patterns
Of the individual nutrients lacking the most in national food supplies, vegetables (for example, carrots, spinach and tomatoes) and milk are major sources of vitamin A in several low-/lower-middle-income countries of Central, Southeastern and Southern Asia, while sweet potatoes have a larger contribution in sub-Saharan Africa. For vitamin B12, for which animal products are the only food sources, milk and seafood (particularly marine-sourced fish) are the major sources. Closing the nutrient gap without an optimization approach, that is, simply increasing total production altogether, would require doubling the global production because cereals dominate existing supply baskets (Fig. 5), and those rich in missing nutrients are underrepresented food groups such as vegetables.
Minimizing GHG emissions while closing the nutrient gap resulted in different food supply priorities at the country level, depending on which nutrients were lacking and specific scenario assumptions. Under all domestic-production-based scenarios, vegetables and vitamin-A-rich roots and tubers (for example, sweet potatoes, yams and cassava) were among the optimal food options that required the largest production increases (in caloric terms). Higher dominance of plant-based products in the optimal baskets of high-income countries (Fig. 4) was a result of an abundance of animal products in their food supplies, meaning that any vitamin B12 gap was effectively very small. In contrast, non-ruminant (for example, duck, rabbit and chicken) meat also featured in optimal solutions across low-, lower-middle- and upper-middle-income groups.
Fig. 4: Optimal food options to close the nutrient gap with the lowest emissions.
For food group compositions, please see Extended Data Table 3. The world map is coloured based on the food groups that require the largest additional supply, either domestically produced or imported, to close the nutrient gap with the lowest emissions. Grey-coloured countries are not included in the analysis. Bar graphs show the total number of countries that have the respective food group in their optimal solution, grouped by their income level. For example, under the first scenario (D-CP-FLW) and in the low-income group, there are 40 countries with roots and tubers in their optimal production baskets, while 5 countries have eggs. See Source Data Fig. 3 for full results and Table 1 for scenario descriptions. Maps were drawn using the tmap R package79.
In the trade scenario (T-CP-FLW), vegetables and other crops such as fruits replaced roots and tubers, and required the largest increases in a greater number of countries. This is because exporting countries in the optimal solution space were often higher-income countries with temperate climates whose production baskets did not include vitamin-A-rich roots and tubers that their partners needed, such as sweet potatoes, in their production and/or trade baskets. Therefore, other sources of such nutrients replaced roots and tubers. Additionally, non-ruminant meat, for example, chicken meat, replaced eggs in several low- and lower-middle-income countries because their exporting countries have more industrial livestock systems with much lower emissions intensities (Source Data Fig. 1).
Between 10% and 23% increase in global caloric production sufficed to close the world’s nutrient gap in 2030 with optimized supply patterns. Under domestic production scenarios, optimal production baskets involved up to 260% and 200% increases in global production of roots and tubers, and eggs, respectively (Fig. 5). Global vegetable production needed to increase by up to 116% under those scenarios. The largest production increase was observed for vegetables, by 48%, when countries resorted to increasing their imports (T-CP-FLW). Additionally, global production of non-ruminant meat showed a 37% increase under the trade scenario. Overall, all scenarios suggest some reduction (up to 17%) in the share of cereals in the global food basket.
Fig. 5: Agricultural production increase for closing the global nutrient gap in 2030 under five scenarios.
a–f, The percentage change is compared to baseline production levels in 2015 (a). Bar graphs show that global production of eggs, vegetables, and roots and tubers needed to increase by more than 100% under domestic-production-based scenarios (b–e), while the increase is between 10% and 23% for total calorie production from all food sources (including oilseeds and sugar crops in 2015) combined. Pie charts illustrate the relative contribution of each food group to global calorie production in 2015 (a) and under five climate-friendly scenarios (b–f). See Table 1 for scenario descriptions. Tool and crop icons from Flaticon.com.
At the country level, optimal food baskets entail production increases, which may be infeasible due to resource limitations. In these cases, it is more realistic to increase both production and imports. As vegetables and roots and tubers include several different individual products with distinct nutritional profiles, compositional change in national production baskets of these food groups can also reduce the need for large-scale increases. For instance, when the proportion of sweet potatoes, a source of vitamin A, is small compared with that of white potatoes, a high supply increase is needed to close the nutrient gap. In contrast, having a larger share of sweet potatoes would require smaller increases in the supply of roots and tubers.
Optimal food supply patterns revealed here help to explain why the climate mitigation potential of trade, compared with domestic production, was higher in the low- and lower-middle-income groups. This occurs because their trade partners included high-income countries where livestock emissions intensity is already much lower. Lastly, halving food loss and waste also performed better in the low-income group because products featuring in the optimal solutions, such as roots and tubers, vegetables and eggs, are subject to higher on-farm and post-harvest (including storage, distribution and processing/packaging) losses in those countries.
Intervention strategies and policy implications
Food loss and waste occur across food supply chains and current research emphasizes customized mitigation approaches3334. Preventing household waste of vegetables and fruits could be a priority in high-income countries where vitamin A is the most commonly lacking nutrient. Even though the emissions intensity per unit of nutrient production from livestock is often lower in developed countries (for example, New Zealand, the United States and France; Fig. 1), high production volume means that absolute levels are still substantial. Therefore, addressing household waste of animal products could also reduce the emissions associated with the nutrient gap closure in their trade partners. In this regard, several low-cost interventions such as smaller portion sizes and encouraging consumers to keep a diary of their waste can be effective34.
Targeting pre-/post-harvest (including processing and distribution) losses in fruits, vegetables, and roots and tubers could be prioritized in low/lower-middle-income countries. Similarly, low-cost and low-energy cold-chain solutions (for example, evaporative coolers) and improved preservation methods for animal products could offer remarkable contribution to closing the vitamin B12 gap. This would require investment in infrastructure, innovation, machinery, packaging and storage, as well as multi-stakeholder cooperation for awareness raising and the transfer of technology and knowledge35. Food shortages, extreme weather events and supply chain disruptions can force farmers to harvest too early or too late, leading to food losses. Hence, other effective measures include establishing standards, price guarantees for farmers, market information systems and public procurement schemes3336. Finally, such preventative measures can be combined with end-of-pipe solutions such as redistribution and donation of surplus from retail outlets and farms via appropriate regulations and financial incentives3336.
In countries with unattained potential yields, smallholder farmers may benefit from higher input use to increase yields3738. Farm diversification and integrating animals (fish and livestock) into crop production creates by-product circularity so that crop residues are used as animal feed, while animal manure is used as fertilizer38. This would also promote diversification of farmers’ incomes and avoid inflated feed costs while closing vitamin A and B12 gaps. For vegetables, access to markets and high-quality seeds are important in boosting productivity, as resource-constrained smallholders account for more than half of the global vegetable production3839. Biofortification programmes may also be effective to enhance nutrient supplies from primary production940. Since the Green Revolution, staples have received most agricultural subsidies, private sector investment and agricultural research focus4142, promoting substantial yield growth in major energy, protein and fat sources (for example, cereals and oilseeds)43. This resulted in a disconnect between the emerging challenges of malnutrition and food policies41. For long-term impacts, investment in infrastructure, innovation, capacity building, and research and development is needed.
Reducing emissions from livestock production has a substantial mitigation potential because animal products are essential to close the vitamin B12 gap. Improved livestock management by means of better feeding practices (for example, improving feed digestibility through lipid supplementation) and veterinary services, improved forage and grassland management, community-based rangeland management, modifying the proportion of the herd dedicated for reproduction, and younger age at first calving are among the least-cost interventions available4445. Access to low interest loans may facilitate such investment, but land tenure and other supporting public policies need to be secured for them to be beneficial for smallholders46. Nevertheless, sustainable livestock production in low-income countries needs more research that considers specific local/regional agri-environmental contexts1545. Extension and support services, production safety net programmes, tenure security, access to markets and affordable loans, and insurance programmes are all promising instruments to attain nutrition-sensitive and sustainable food systems4046.
Market forces have the potential to promote the best use of food-related resources if policies are tailored accordingly47. Currently, international food trade improves nutrient availability to a much lesser extent in low- and lower-middle-income countries with inadequate production7. Compounding this issue, emissions intensity of nutrient production is much higher in most of these countries due to low productivity (Fig. 1). Hence, current trade patterns are not optimal for either nutrition security or climate mitigation. Additionally, resource limitations, coupled with climate change, make several low-income countries increasingly dependent on imports48. Therefore, promoting trade between countries with production surpluses and lower emissions intensities and countries with production deficits and higher emissions intensity could help close nutrient gaps at lower emissions. Average tariffs applied on dairy, meat and seafood—high in most low-income countries7—could be selectively lowered to close vitamin B12 gaps49. In addition, internalizing the emissions of food (for example, border tax adjustments and climate measures in regional trade agreements) could facilitate trade patterns where nutrients flow from countries with lower emissions intensities to those with higher emissions intensities29. This would reverse the situation where less than half of the flows occurred from countries with lower emissions intensity to those with higher intensity, rendering the global nutrient trade ineffective as a mitigation mechanism50.
Discussion
There is growing consensus on the importance of food systems transformation, involving production- and consumption-based interventions, to tackle climate change and malnutrition simultaneously118. Targeted metrics, underpinned by comprehensive datasets and rigorous modelling, are crucial for effective policies1651. Here, we developed a multi-dimensional, high-product resolution, country-level dataset that combines dietary nutrient requirements, production, trade and resulting non-CO2 GHG emissions. With an optimization model, we evaluated the minimum emissions associated with closing the nutrient gap (for energy and six nutrients) across five climate-friendly intervention scenarios, 196 agri-food products and 128 countries. In contrast to wholesale dietary changes, a common approach in the literature831, we identified country-specific priority food sources that close the nutrient gap with the lowest emissions.
Interventions to address inefficiencies at home, for example, reducing loss and waste and livestock emissions, had a higher potential to mitigate the emissions compared with importing from the least emissions-intensive partner. Moreover, when emissions from transportation are incorporated, which may be up to 3 Gt CO2 yr−1 globally52, closing the nutrient gap via trade may offer lower mitigation potentials. Our results indicate that halving food loss and waste and improving agricultural productivity together can reduce the emissions by up to 42% while closing the nutrient gap at the same time. At-home interventions are also likely to be pro-poor because, otherwise, import substitution may harm producers in the lower-income countries18.
In terms of food sources, our high-product resolution, as opposed to the oft-cited8 aggregate food balance sheets, offers a more nuanced evaluation of not only the nutrient gap but also priority food sources by country. Increasing production of vegetables, roots and tubers, and non-ruminant meat would help close the world’s nutrient gap with the smallest emissions and result in 10–23% increase in global caloric production by 2030. This translates into a reduction in the share of cereals in food production/supply baskets, as is also suggested in the literature818313953. Our findings confirm the importance of addressing losses in and productivity of vegetable and fruit production131854.
Our model assumes that any increase in production is associated with an increase in total factor productivity and excludes emissions from land-use change due to a lack of high-resolution commodity-specific datasets that allow prospective mapping of different products with cropland and pasture change and associated changes in biomass—as is common practice in dynamic modelling approaches1931. Hence, the consequent emissions are potentially underestimated, particularly in regions nearing attainable yields. Nevertheless, abandoned farmland could also be re-cultivated in regions such as North America, Europe, Eastern and Central Asia, and Oceania. Our results can also be integrated with crop models to map where climate- and nutrition-sensitive cropping patterns are suitable.
Our model can be extended to include other environmental and nutritional constraints141516. Given our focus on the double burden of malnutrition, that is, undernutrition and micronutrient deficiencies, we neither set an upper limit on nor aimed to reduce calories to address overconsumption. However, we excluded energy-rich and micronutrient-poor sources, based on the way they are currently consumed, for example, sugar and oilseed crops, from optimization input because we aimed to address the double burden of malnutrition simultaneously and such food sources tend to encourage caloric overconsumption. This is similar to the ‘no sugar’ scenario in the literature8. Consequently, the required increase in global caloric production was between 10% and 23% due to selection of nutrient-rich foods such as vegetables by the optimization model. On the other hand, tackling overconsumption, particularly of animal products, in higher-income countries is paramount to address the global syndemic183155, with a caveat for low-income countries where undernutrition and micronutrient deficiencies necessitate increased consumption of animal products2131. Lastly, our scope is limited to availability, one of the four pillars (that is, availability, accessibility, utilization and stability) of nutrition security, and that adequate supply does not necessarily equate to adequate intake8. Nonetheless, without adequate supply, equal distribution alone does not suffice, with the caveat that complementary measures, for example, fortification, are also important to enhance nutrient supplies, such as iron and vitamin A, particularly from cereals and vegetable oils.
Climate change creates further challenges for nutrient availability and economic access to nutritious food, especially in low-income countries55. Economic responses and potential shifts in price, supply, demand and access should therefore be considered for integrated policies. For instance, productivity gains are found to encourage land expansion in non-food sectors18, which necessitates complementary instruments such as taxation1844. Despite the limitations regarding rebound effects in static models like ours, value-based economic models lack other dimensions such as biophysical limits and preservation of mass and energy balances56. By incorporating livestock production and crop models, our scenarios can directly inform the sustainable food systems debate. For instance, given the urgency of the global syndemic and the difficulties in achieving large-scale dietary shifts, which food supplies need to be scaled up to contribute to the food system transformation?
Without imposing drastic changes in the current food baskets (and hence, diets), increasing production of less emissions-intensive products that provide the set of nutrients lacking in the national supplies could close nutrient gaps and achieve substantial climate mitigation. This implies a different picture compared with that suggested by most global food demand models. While the literature suggests up to 30% increase in income-driven total food (calories) demand between 2010 and 203017, we offer an approach based on physiological needs47. Our findings suggest that less than 24% production increase from 2015 would suffice in closing the world’s nutrient gap by 2030, with 42% less GHG emissions compared with the BaU scenario, in which the nutrient gap remains.
Methods
Dietary nutrient production, supply and gap
Many countries do not produce and/or import enough nutrients to meet the recommended nutrient intake, an estimation of daily intake needed to provide the requirements of 97.5% of healthy individuals in a population group differentiated by age and sex, for their populations7. We obtained the dataset for country-level total nutrient production, trade, and supply from ref. 7 and methodology to estimate population-level nutrient requirements from ref. 57. These cover nutrient supply (energy, protein, iron, zinc, vitamin A, folate and vitamin B12) from crops, livestock and seafood. Nutrient supply estimates account for the share of crops that are fed to animals (see ref. 57 for more details). For better representation of nutrient supplies, we converted oilseeds into vegetable oils based on the share of oilseeds processed over their national supply by respective product (Extended Data Table 4). Finally, we subtracted the share of non-food uses, by using food balance sheets58, due to the use of large volumes of oilseeds and their oil derivatives (for example, linseed, soybean oil and rapeseed oil) for fuel purposes (that is, biodiesel).
We estimated future population-level nutrient requirements based on median variant population forecasts for 203012 and following an established nutrient adequacy approach that compares country-level nutrient supplies against age- and sex-adjusted requirements57. Because the requirements for zinc and iron almost double for zinc and quadruple for iron during pregnancy, lower fertility rates decreased the aggregate requirements for these nutrients.
The difference between nutrient supplies (that is, supply = domestic production + net trade − food loss and waste) and requirements provided country-level supply gaps for each nutrient, as shown in equation (1). We quantified the nutrient supply in 2015 under current food loss and waste patterns (that is, FLW) using regional and food group specific loss and waste rates provided by ref. 59 where p represent nutrients, a represents countries and c represents food products (C being the total number of food products). In equation (1), nutrient production, that is, \({\mathrm{NP}}_{a,c}^p\), refers to the weight before farm losses minus exports, which we adjusted by adding on values reported as excluding harvesting losses by the Food and Agriculture Organization of the United Nations (FAOSTAT).
$$\begin{array}{l}{\mathrm{Nutrient}}\,{\mathrm{gap}}_a^p = {\mathrm{Population}}\,{\mathrm{requirements}}_a^p-{\mathrm{Nutrient}}\,{\mathrm{supply}}_a^p\\ {\mathrm{where}}\\ {\mathrm{Nutrient}}\,{\mathrm{supply}}_a^p = \mathop {\sum}\nolimits_{c = 1}^C {\left[ {{\mathrm{NP}}_{a,c}^p \times {\mathrm{LW}}_{a,c}^{{\mathrm{FL}}} \times {\mathrm{LW}}_{a,c}^{{\mathrm{PHL}}} \times {\mathrm{LW}}_{a,c}^{{\mathrm{PPL}}} \times {\mathrm{LW}}_{a,c}^{{\mathrm{DL}}} \times {\mathrm{LW}}_{a,c}^{{\mathrm{HW}}}} \right.} \\ - \,{\mathrm{Exports}}_{a,c}^p \times {\mathrm{LW}}_{a,c}^{{\mathrm{PHL}}} \times {\mathrm{LW}}_{a,c}^{{\mathrm{PPL}}}\\ \left. { + \, {\mathrm{Imports}}_{a,c}^p \times{\mathrm{LW}}_{a,c}^{{\mathrm{DL}}} \times {\mathrm{LW}}_{a,c}^{{\mathrm{HW}}} } \right]\end{array}$$
(1)
Equation (1) considers the five different types of loss and waste (LW) as proportional values. For example, if farm loss is 5% of the production, then the LW for farm loss, \({\mathrm{LW}}_{a,c}^{{\mathrm{FL}}}\), becomes 0.95. Types of loss and waste include: farm loss (FL), post-harvest loss (PHL), processing and packaging loss (PPL), distribution loss (DL) and household waste (HW). We assumed constant supply, hence supply gaps, for each scenario. Consequently, the emissions associated with closing the nutrient gap varied by (1) how much production increase is needed to meet these gaps and/or (2) emissions intensity of production.
GHG emissions components
We developed a composite indicator to quantify farm-gate non-CO2 emissions intensity of national agricultural production to be used as input for optimization (equation (5)). Ref. 60 provides a comprehensive dataset of life cycle emissions for 188 commodities from 119 countries. However, most crops (excluding maize) and livestock products in that dataset originate from major producing countries, which does not allow reveal differences across countries from different income levels. Therefore, we chose to use process-based emissions estimates following the IPCC (2006)28 Tier I approach and account for the heterogeneity in emissions factors across different agro-ecological regions.
We obtained the data from FAOSTAT database (2020 revision)61. FAOSTAT emissions estimates follow the Tier I approach of the IPCC guidelines for national GHG inventories28. We focused on non-CO2 GHG emissions from agriculture, excluding emissions from land-use change (for example, forest conversion to cropland/grassland) due to the lack of product-specific attribution that would allow us to link production growth (for example, increase in vegetable production) with corresponding land-use change. In contrast, emissions from crop residues and fertilizer use are directly attributable to source products.
We included 128 countries (Extended Data Table 1) with comparable data for production, product-level bilateral trade and GHG emissions. We used the average GHG emissions values for the 2013–2015 period to smooth out yearly fluctuations61. We used the most widely used metric, 100-year global warming potential (GWP100), from the IPCC Sixth Assessment Report (AR6) in our calculations (CH4-non-fossil origin: 27; N2O: 273)62.
Our emissions scope covered rice cultivation (CH4 from decomposition of organic matter in paddy fields), synthetic fertilizers (direct emissions of N2O from nitrification and de-nitrification, and indirect emissions from volatilization/re-deposition and leaching processes), crop residues (N2O from decomposition of residue left on soils), burning crop residues (CH4 and N2O from the combustion of crop residues), enteric fermentation (CH4 from rumination in the digestive track), manure management (CH4 and N2O from manure decomposition), manure applied to soils (direct emissions of N2O from nitrification and de-nitrification, indirect emissions from volatilization/re-deposition and leaching processes from manure applied to cropland) and manure left on pasture (direct emissions of N2O from nitrification and de-nitrification, and indirect emissions from volatilization/re-deposition and leaching processes from manure left on pasture by grazing animals)61. Manure applied to soils is regarded as organic fertilizer and allocated to emissions from crop production.
Although recent assessments show that the share of total global food system emissions related to energy use is growing363, we excluded emissions from energy use in this study because product-specific global data were lacking. Interventions to mitigate emissions from energy use are similar to other economic sectors like industry and transport such that their emissions intensities are primarily determined by energy mix, which can be decarbonized via increasing the amount of renewable energy and technological innovation such as the electrification of transport.
Crops
Crop emissions cover the non-feed portion used for direct human consumption and are differentiated by product and country. The share of crops fed to livestock was obtained from food balance sheets of the FAOSTAT (see ref. 57 for more details). We also subtracted the quantity of crops used as feed in aquaculture (that is, aquafeed). As not all aquaculture is fed, we used several different sources to ensure reliable and comprehensive emissions accounting. Ref. 64 provides the share of aquaculture that is fed for 11 fish species: carps, tilapia, shrimp, catfishes, marine fish, salmon, freshwater crustaceans, other diadromous fishes, milkfish, trout and eel. The study also provides data on the share of fish sources (fishmeal and fish oil) in their compound feeds in addition to fish-in-fish-out ratios. Based on this information, we calculated the share coming from non-fish sources because we are interested in GHG emissions from growing crops used as feed. Ref. 65 estimated the share of individual crop sources (maize, soybeans, wheat, pulses and other oil crops) in aquafeed at the regional level. However, aquafeed does not entirely consist of these plant sources. To correct for this, we combined ratios65 of this regional crop-based feed with the share of plant sources in aquafeed64. Extended Data Table 8 shows the consequent aquafeed use ratios for products from maize, soybeans, wheat, pulses and other oil crops for fed aquaculture.
Finally, we estimated the primary crop equivalent of aquafeed sources from soybeans and other oil crops often used as by-products (for example, flour and oil cake) in aquafeed66. Extended Data Table 9 shows the conversion factors that are used in these calculations67. We used the weighted average of top producers (whose cumulative production accounts for ≥80% of global production). We assumed that other aquafeed uses of oil crops were largely rapeseed and sunflower seed. Consequently, we derived aquafeed-to-supply ratios per crop and country (βca), which we used in estimating the nutrient gap.
Total GHG emissions from crop residues (\({\mathrm{GHG}}_{s,a,c}^{{\mathrm{Crop}}\,{\mathrm{residues}}}\)), burning crop residues (\({\mathrm{GHG}}_{s,a,c}^{{\mathrm{Crop}}\,{\mathrm{residues}},\,{\mathrm{burning}}}\)) and rice cultivation (\({\mathrm{GHG}}_{s,a,c}^{{\mathrm{Paddy}}\,{\mathrm{rice}}}\)) were calculated for each crop product (c), country (a) and intervention scenario (s). We attributed emissions from synthetic (\({\mathrm{GHG}}_{s,\,a,c}^{{\mathrm{Synthetic}}\,{\mathrm{fertilizer}}}\)) and organic (\({\mathrm{GHG}}_{s,a,c}^{{\mathrm{Organic}}\,{\mathrm{fertilizer}}}\)) fertilizers to respective food sources based on data on fertilizer use by crops68 as detailed in Extended Data Table 2. Country-level total crop emissions were calculated as per equation (2) where θca represents livestock feed-to-supply ratio for a given country and product (calculated based on food balance sheets58 and βca shows aquafeed (that is, feed for aquaculture)-to-supply ratio (see Extended Data Table 7) for a given country and product6465. The ratio of crops that are not diverted to livestock or aquaculture is given by γ.
$$\begin{array}{*{20}{l}} {{\mathrm{GHG}}_{s,a}^{{\mathrm{Crops}}}} \hfill & = \hfill & {\mathop {\sum}\nolimits_{c = 1}^C {{\mathrm{GHG}}_{s,a,c}} \,{\mathrm{and}}} \hfill \\ {{\mathrm{GHG}}_{s,a,c}} \hfill & = \hfill & {\gamma _{c,a} \times \left( {{\mathrm{GHG}}_{s,a,c}^{{\mathrm{Crop}}\,{\mathrm{residues}}}} \right.} \hfill \\ {} \hfill & {} \hfill & { + \, {\mathrm{GHG}}_{s,a,c}^{{\mathrm{Crop}}\,{\mathrm{residues}},\,{\mathrm{burning}}}} \hfill \\ {} \hfill & {} \hfill & { + \, {\mathrm{GHG}}_{s,a,c}^{{\mathrm{Synthetic}}\,{\mathrm{fertilizer}}}} \hfill \\ {} \hfill & {} \hfill & { + \, {\mathrm{GHG}}_{s,a,c}^{{\mathrm{Organic}}\,{\mathrm{fertilizer}}}} \hfill \\ {} \hfill & {} \hfill & {\left. { + \, {\mathrm{GHG}}_{s,a,c}^{{\mathrm{Paddy}}\,{\mathrm{rice}}}} \right)} \hfill \\ {{\mathrm{where}}\,\gamma _{c,a}} \hfill & = \hfill & {1 - (\theta _{c,a} + \beta _{c,a})} \hfill \end{array}$$
(2)
Livestock
Livestock emissions included those from enteric fermentation (\({\mathrm{GHG}}_{s,a,c}^{{\mathrm{Enteric}}\,{\mathrm{fermentation}}}\)), manure management \(\left( {{\mathrm{GHG}}_{s,a,c}^{{\mathrm{Manure}}\,{\mathrm{management}}}} \right)\), manure left on pasture \(\left( {{\mathrm{GHG}}_{s,a,c}^{{\mathrm{Manure}}\,{\mathrm{left}}\,{\mathrm{on}}\,{\mathrm{pasture}}}} \right)\), fertilizers applied to grassland (\({\mathrm{GHG}}_{s,a,c}^{{\mathrm{Synthetic}}\,{\mathrm{fertilizer}}}\) and \({\mathrm{GHG}}_{s,a,c}^{{\mathrm{Organic}}\,{\mathrm{fertilizer}}}\)) and feed crops (\(\theta _{c,a} \times {\mathrm{GHG}}_{s,a,c}^{{\mathrm{Crops}}}\)) that are grown domestically (see Extended Data Table 6 for feed crop consumption by livestock group). We allocated the emissions from crops fed to livestock (that is, share of crops used as livestock and poultry feed according to food balance sheets of the FAOSTAT (see ref. 57) from crop emissions to the respective livestock type based on the relative amount of food crops estimated to be consumed by ruminants, pigs and poultry69 (consumption share by animal given in Extended Data Table 5). Fertilizer use by crops (Extended Data Table 2)68 also includes fertilizers applied to grassland. Hence, we linked the relevant share of emissions from synthetic (\({\mathrm{GHG}}_{s,a,c}^{{\mathrm{Synthetic}}\,{\mathrm{fertilizer}}}\)) fertilizers to respective animals based on the share of grass-fed animals by animal type, for example, ruminants, pigs and poultry69. Total country-level livestock emissions were calculated as per equation (3) for each livestock product (c), country (a) and intervention scenario (s):
$$\begin{array}{*{20}{l}} {{\mathrm{GHG}}_{s,a}^{{\mathrm{Livestock}}}} \hfill & = \hfill & {\mathop {\sum}\nolimits_{c = 1}^C {{\mathrm{GHG}}_{s,a,c}} ,} \hfill \\ {{\mathrm{GHG}}_{s,a,c}} \hfill & = \hfill & {{\mathrm{GHG}}_{s,a,c}^{{\mathrm{Enteric}}\,{\mathrm{fermentation}}}} \hfill \\ {} \hfill & {} \hfill & { + \, {\mathrm{GHG}}_{s,a,c}^{{\mathrm{Manure}}\,{\mathrm{management}}}} \hfill \\ {} \hfill & {} \hfill & { + \, {\mathrm{GHG}}_{s,a,c}^{{\mathrm{Manure}}\,{\mathrm{left}}\,{\mathrm{on}}\,{\mathrm{pasture}}}} \hfill \\ {} \hfill & {} \hfill & { + \, {\mathrm{GHG}}_{s,a,c}^{{\mathrm{Synthetic}}\,{\mathrm{fertilizer}}}} \hfill \\ {} \hfill & {} \hfill & { + \, \theta _{c,a} \times {\mathrm{GHG}}_{s,a,c}^{{\mathrm{Crops}}}} \hfill \end{array}$$
(3)
Emissions intensity of nutrient production
Emissions intensity of nutrient production, for visualization purposes in Fig. 1, is calculated by simply dividing total GHG emissions (\({\mathrm{GHG}}_{s,a,c}\)) by individual nutrient (for example, protein) production (see Dietary nutrient production, supply and gap). In order to construct our optimization model, we calculated GHG emissions intensity (that is, \(I_{s,a,c}^{p = {\mathrm{energy}}}\)) of unit caloric availability (that is, domestic production minus food loss and waste59):
$$I_{s,a,c}^{p = {\mathrm{energy}}} = {\mathrm{GHG}}_{s,a,c}/{\mathrm{Nutrient}}\,{\mathrm{availability}}_{s,a,c}^{p = {\mathrm{energy}}}$$
(4)
Scenario description
We introduced five climate-friendly intervention scenarios related to crop and livestock productivity, food loss and waste, and trade. Based on the assumptions imposed by each scenario, emissions intensity of energy availability (represented by \(I_{s,a,c}^{p = {\mathrm{energy}}}\)) changed and we calculated the agricultural production required to meet or exceed population-level requirements for energy, protein, iron, zinc, vitamin A, vitamin B12 and folate (see Dietary nutrient production, supply and gap) with the minimum emissions accordingly.
Current food loss and waste and productivity patterns
Under the current food loss and waste and productivity patterns (D-CP-FLW), emissions intensity (\(I_{s,a,c}^{p = {\mathrm{energy}}}\) where s = D-CP-FLW) was calculated based on current emissions, \({\mathrm{GHG}}_{s,a,c}\), and 2015 productivity patterns.
Halving food loss and waste
Losses arise before and after harvest, and during processing, packaging and distribution, while waste occurs at the household and retail level. Because there is no systematic evaluation of the extent of abatable food loss and waste in different regions (as there is for productivity), we assumed a 50% reduction in the half loss and waste scenario (HLW) in line with SDG 12 Target 12.3, which aims to halve food waste at retail and consumer levels70. For sugar crops, we assumed the same loss and waste rate as oil crops and pulses.
Depending on the stage at which loss/waste occurred as per the description in ref. 59, nutrient supply was calculated by halving loss and waste rates (equation (1)). As a result, for example, if farm loss is originally 5%, halving farm loss would result in \({\mathrm{LW}}_{a,c}^{{\mathrm{FL}}} = 0.975\). Because we assumed constant nutrient supply and gaps, halving loss and waste resulted in a reduction of baseline emissions and additional emissions through lower emissions intensity of energy availability \(\left( {I_{s,a,c}^{p = {\mathrm{energy}}}} \right)\).
Improved productivity
Under current productivity patterns (CP), we assumed the current (2013–2015) emissions intensity of production (that is, nutrient content of a given product/production-based GHG emissions) for every country (a). For increased productivity scenarios (D-IP-FLW and D-IP-HLW), we followed slightly different approaches for crops and livestock. For crops, we considered yield gap closure. We used the Global Agro-Ecological Zones (GAEZv3) model outputs, which include spatially resolved estimates of potential yields for dozens of individual crops under specific agro-climatic, soil, terrain and management conditions71. To quantify yield gaps, we compared historical crop yields with potential yields under high-input use. We then estimated the additional nitrogen (N) fertilizer requirement to achieve these high-input yields estimated at the regional level37. We assumed that yield increases were achieved by the increased use of synthetic fertilizers only and estimated the resulting emissions. Potential yields and fertilizer requirements were estimated at the regional level and downscaled to derive country-level estimates (see Supplementary Tables 1 and 2). Quantification of fertilizer GHG emissions followed the Tier 1 approach28, which assumes a default emissions factor of 0.01 kg N2O–N (kg N)−1.
For livestock productivity, we used potential mitigation in emissions intensities estimated by the Global Livestock Environmental Assessment Model (GLEAM)44. The model quantifies the environmental impacts of livestock production over its life cycle and draws on adaptation and mitigation scenarios for a more sustainable livestock sector. It provides the scope for mitigation in the livestock sector globally across five animal species and as case studies for five world regions that are applicable in the medium term (for example, up to two decades) (Extended Data Table 5). Beyond cattle, pigs and poultry covered by the GLEAM model, camelids are also a good source of nutrients in certain regions. They mostly occur in marginal lands of arid countries in Africa (for example, camels in North Africa and Sahelian countries), Asia (camels in West and Central Asia) and South America (alpacas and llamas in the Andean region), and are often kept for draught power. In this regard, any food product, such as milk and meat, provides additional income rather than being the main source of income for the farm. Hence, we did not treat them as regular farm animals and assumed no increase in their emissions intensity. The mitigation potentials were quantified based on constant output in the GLEAM model.
Feed quality and animal health and husbandry are critical factors for improving livestock productivity in low-/lower-middle-income countries72. We chose this approach to avoid imposing the same productivity patterns on developing countries as industrialized countries (that is, intensive production systems). The model also gives a breakdown of impacts from different intervention scenarios at the animal, herd, production unit and supply chain levels. These include optimized feed digestibility, animal health and mortality, genetics, grassland management, and manure management for ruminants and monogastrics, in addition to energy efficiency and anaerobic digestors for pig production44. We only considered interventions in feed quality (for example, digestibility), grazing management, manure management and reduced mortality. Feed conversion ratios are assumed to be constant and emissions from feed crops changed in accordance with crop-based emissions under the combined productivity scenarios.
The emissions intensity (\(I_{s,a,c}^{p = {\mathrm{energy}}}\), where s = D-IP-FLW) was calculated based on changes in crop- and livestock-related GHG emissions (\({\mathrm{GHG}}_{s,a,c}\)) and changes in crop production per unit emissions associated with crop yield gap closure.
Improved productivity and half food loss and waste
Under D-IP-HLW scenario, emissions intensity (\(I_{s,a,c}^{p = {\mathrm{energy}}}\), where s = D-IP-HLW) involved changes in crop- and livestock-related GHG emissions and lower production needs due to larger shares of production being available for consumption.
Domestic production versus imports
In contrast to the scenarios based on domestic production, the objective function included emissions intensity of production in the partners exporting to the given country a under T-CP-FLW. We constructed the objective function with the existing bilateral trade partnerships and baskets based on data provided by the FAO under the detailed trade matrix domain73. Any increase in the export volume of their trading partners was assumed to be met by corresponding increases in production. Hence, apparent consumption (for example, supply) in exporting countries remained unchanged. \(I_{s,a,c}^{p = {\mathrm{energy}}}\), where s = T-CP-FLW, was equal to \(I_{s,a,c}^{p = {\mathrm{energy}}}\), where s = D-CP-FLW, for all countries and products.
Optimization of food supply patterns to close the nutrient gap
We applied linear programming to identify the additional production (measured in caloric terms) required for a given country (a) and product (c) under a certain intervention scenario (s) to close the nutrient gap while minimizing food system non-CO2 GHG emissions. The objective function minimized GHG emissions from additional production such that the supply of all nutrients was adequate to meet national dietary requirements based on the existing production and bilateral trade baskets for each country. Therefore, additional production refers to domestic production of a given country a under domestic-production-based scenarios (D-CP-FLW, D-CP-HLW, D-IP-FLW and D-IP-HLW). In contrast, under the trade scenario (T-CP-FLW), it refers to additional production in partners that country a imports from. Similarly, \(I_{s,a,c}^{p = {\mathrm{energy}}}\) represents domestic emissions intensity under domestic-production-based scenarios, whereas it represents the vector of emissions intensity in partners exporting to the given country a. Composition of production and trade baskets (that is, the number of individual products) remained the same, although relative contribution by food source changed, with the assumption that diets do not observe radical changes in their composition (for example, complete elimination of certain food groups from diets and introduction of novel food products that are absent from current food baskets).
$$\begin{array}{l}{\mathrm{GHG}}_{s,a}^{{\mathrm{Min}}} = \\ {{{\mathrm{minimize}}}}:\mathop {\sum}\nolimits_{c = 1}^C {I_{s,a,c}^{p = {\mathrm{energy}}} \times } \Delta {\mathrm{NP}}_{s,a,c}^{p = {\mathrm{energy}}}\\ {{{\mathrm{subject}}}}\,{{{\mathrm{to}}}}\,{{{\mathrm{the}}}}\,{{{\mathrm{constraints}}}}:\\ {\mathrm{NER}}_{s,a,c}^p \times \Delta {\mathrm{NP}}_{s,a,c}^{p = {\mathrm{energy}}} \ge {\mathrm{NG}}_a^p\,{{{\mathrm{for}}}}\,{{{\mathrm{all}}}}\,{{{\mathrm{nutrients}}}}\,p\\ {\mathrm{and}}\\ \Delta {\mathrm{NP}}_{s,a,c}^{p = {\mathrm{energy}}} \ge 0\end{array}$$
(5)
where \({\mathrm{NP}}_{s,a,c}^{p = {\mathrm{energy}}}\) is the nutrient production of dietary energy (that is, calories). \(I_{s,a,c}^{p = {\mathrm{energy}}}\) was calculated as described in equation (4). \({\mathrm{NER}}_{s,a,c}^p\) is the nutrient-to-energy ratio for each nutrient p (that is, energy, protein, iron, zinc, vitamin A, vitamin B12 and folate), each country a and food product c, and \({\mathrm{NER}}_{s,a,c}^p\) = 1 for p = energy. \({\mathrm{NG}}_a^p\) is the nutrient gap for each nutrient p and country a. The general form of equation (5) was applied to each country for every scenario. We used the HiGHS solver from the linprog package from the SciPy library of Python74, which implements the interior-point method and features parallel programming.
Following equations (1)–(5), the emissions associated with closing the nutrient gap globally are the sum of baseline emissions in 2015 and additional emissions from the production increase required to close the nutrient gap (where A is the total number of countries):
$${\mathrm{GHG}}_s^{{\mathrm{Total}}} = \mathop {\sum}\nolimits_{a = 1}^A {{\mathrm{GHG}}_{s,a}^{{\mathrm{Crops}}}} + \mathop {\sum}\nolimits_{a = 1}^A {{\mathrm{GHG}}_{s,a}^{{\mathrm{Livestock}}}} + \mathop {\sum}\nolimits_{a = 1}^A {{\mathrm{GHG}}_{s,a}^{{\mathrm{Min}}}}$$
(6)
Uncertainty estimates
There are inherent uncertainties with our underlying data and approach. Uncertainty ranges are unknown for production/trade data that originate from FAOSTAT75, but the Tier 1 approach to estimate GHG emissions have known uncertainty ranges related to default emissions factors28. Therefore, in addition to the default factors, we included lower and upper bounds for emissions factors used in the Tier 1 approach (that is, Emissions = Activity data × Emissions Factor) to estimate GHG emissions28. The IPCC (2006)28 provides lower and upper bounds either as a percentage of deviation from the default value for some emissions sources (for example, enteric fermentation) or as an absolute value for others (for example, rice cultivation). Additionally, some emissions sources (for example, N2O from managed soils) have direct and indirect emissions. In that case, there is also uncertainty associated with the fraction of leaching and volatilization. The IPCC methodology suggests leaching only in regions where runoff occurs. However, the FAOSTAT assumes that leaching occurs in all regions due to absence of region-specific information61. Our estimates encompass uncertainty associated with both factors and converted absolute values to percentage change for ease of calculation (for example, percent deviation from the default value). Specific emissions factor ranges for each emission source are presented in Extended Data Table 9.
To construct uncertainty ranges, we assumed that our three-point estimates (results based on the default, lower and upper bound emissions factors) follow a program evaluation and review technique distribution. This distribution is defined by the most likely (that is, mode) and extreme (the minimum and the maximum) values a variable can take. We used the qpert function in the mc2d package in the R software to estimate the 25th and 75th percentiles76.
Paris Agreement and allowable food production emissions
To interpret our findings, we compared our results against the allowable emissions range of pathways compliant with the Paris Agreement32. Paris Agreement-compatible allowable emissions are estimated by selecting the Paris-compliant pathways from the full ensemble of pathways underlying the IPCC AR677. The ensemble was filtered by criteria for efforts to limit global warming to 1.5 °C (for example, <66% chance to overshoot of 1.5 °C), holding global warming well below 2 °C (for example, 90% chance) and achieving net zero emissions in the second half of the twenty-first century, in order to remain consistent with the Paris Agreement32.
In line with our scope, we considered only the CH4 and N2O emissions from the AFOLU sector in 2030. To derive the CO2eq emissions, we used the updated GWP100 factors from the IPCC AR662. To further ensure compatibility with our study scope, we adjusted the allowable AFOLU emissions based on the global population covered in this study (89% of the global population) under the assumption of fair-share GHG per capita. Furthermore, as we did not include N2O emissions from drained organic soils, which represents 2% of the global non-CO2 emissions, we rationalized the allowable non-CO2 emissions accordingly. In addition, our focus on food crops (that is, excluding fibre crops) corresponds to 99% of total agricultural production by weight75. Given that around 20% of agricultural non-CO2 emissions arose from crop production in 2013–2015, this amounted to ~0.3% downscaling in boundaries. Consequently, we scaled down the 25th (5.03 Gt CO2eq yr−1), 50th (5.43 Gt CO2eq yr−1) and 75th (6.17 Gt CO2eq yr−1) percentile values by 14% for comparison with our findings. The consequent allowable non-CO2 emissions range was 4.33–5.31 Gt CO2eq yr−1.
Assumptions and limitations
Recent research suggests that CO2 warming equivalents (CO2we), following the newly established GWP* model, may better account for the behaviour of the short-lived climate pollutants, such as CH4, in projecting temperature effects78. Given that the increase in CH4 emissions was smaller when optimization was introduced (Extended Data Fig. 1), compared with BaU, due to relatively smaller gaps in vitamin B12 supplies and associated increase in livestock production (the only source of vitamin B12 as an input to optimization), future research could enhance the understanding of temperature effects of decreasing growth in CH4 emissions by using GWP*. To provide an illustrative figure, we presented our optimization results based on CO2 warming equivalents using the GWP* model in the Supplementary Information. It suggests that despite smaller warming potentials suggested by GWP* because CH4 emissions either decrease (for example, with productivity improvements) or nearly stabilize (for example, with optimization only) compared with 2015, relative performance of our climate-friendly scenarios remains robust to the chosen equivalence method, that is, CO2we or CO2eq.
We acknowledge the complexity of adequate nutrition, which depends on a delicate balance of a diverse set of nutrients as well as other socioeconomic determinants and underlying health conditions that are not captured in this study. Similarly, fortification and supplementation, presented as food-based intervention options to fill the nutrient gap in diets, are not considered due to a lack of reliable production/trade data across all countries and products included in our study. More importantly, it is harder to estimate the contribution of fortification in countries where the nutrient gap is highest, such as low-income countries with a high share of rural population, because fortified foods may not be accessible in rural areas and implementation is difficult in small-scale mills9. See the Supplementary Information for a thorough discussion of our assumptions and their limitations.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
All input data are publicly available through online sources as given in the references. All other data supporting the findings of this study are available within the paper. Source data are provided with this paper.
Code availability
Codes related to optimization are publicly available via https://github.com/OzgeGe/opt.git. Further information is available upon request.
Change history
- 12 January 2023
A Correction to this paper has been published: https://doi.org/10.1038/s43016-023-00693-1
References
- Swinburn, B. A. et al. The global syndemic of obesity, undernutrition, and climate change: the Lancet Commission report. Lancet 393, 791–846 (2019).
- Popkin, B. M. et al. Individuals with obesity and COVID-19: a global perspective on the epidemiology and biological relationships. Obes. Rev. 21, e13128 (2020).
- Crippa, M. et al. Food systems are responsible for a third of global anthropogenic GHG emissions. Nat. Food 2, 198–209 (2021).
- Clark, M. A. et al. Global food system emissions could preclude achieving the 1.5° and 2 °C climate change targets. Science 370, 705–708 (2020).
- Rogelj, J. et al. in Global Warming of 1.5°C (eds Masson-Delmotte, V. et al) 93–174 (Cambridge University Press, 2022); https://doi.org/10.1017/9781009157940.004
- Zhang, P. et al. Global healthcare expenditure on diabetes for 2010 and 2030. Diabetes Res. Clin. Pract. 87, 293–301 (2010).
- Geyik, O., Hadjikakou, M., Karapinar, B. & Bryan, B. A. Does global food trade close the dietary nutrient gap for the world’s poorest nations? Glob. Food Sec. 28, 100490 (2021).
- Smith, N. W., Fletcher, A. J., Dave, L. A., Hill, J. P. & McNabb, W. C. Use of the DELTA model to understand the food system and global nutrition. J. Nutr. 151, 3253–3261 (2021).
- Bose, I., Baldi, G., Kiess, L. & Pee, S. The “fill the nutrient gap” analysis: an approach to strengthen nutrition situation analysis and decision making towards multisectoral policies and systems change. Matern. Child Nutr. 15, e12793 (2019).
- Bajželj, B., Laguzzi, F. & Röös, E. The role of fats in the transition to sustainable diets. Lancet Planet. Health 5, e644–e653 (2021).
- Herrero, M. et al. Biomass use, production, feed efficiencies, and greenhouse gas emissions from global livestock systems. Proc. Natl Acad. Sci. USA 110, 20888–20893 (2013).
- World Population Prospects 2019: Standard Projections (United Nations Department of Economic and Social Affairs, 2019); https://population.un.org/wpp/Download/Standard/Interpolated/
- Mason-D’Croz, D. et al. Gaps between fruit and vegetable production, demand, and recommended consumption at global and national levels: an integrated modelling study. Lancet Planet. Health 3, e318–e329 (2019).
- DeFries, R. et al. Metrics for land-scarce agriculture. Science 349, 238–240 (2015).
- Mehrabi, Z., Gill, M., Wijk, M., van, Herrero, M. & Ramankutty, N. Livestock policy for sustainable development. Nat. Food 1, 160–165 (2020).
- Fanzo, J. et al. Rigorous monitoring is necessary to guide food system transformation in the countdown to the 2030 global goals. Food Policy 104, 102163 (2021).
- van Dijk, M., Morley, T., Rau, M. L. & Saghai, Y. A meta-analysis of projected global food demand and population at risk of hunger for the period 2010–2050. Nat. Food 2, 494–501 (2021).
- Kuiper, M. & Cui, H. D. Using food loss reduction to reach food security and environmental objectives – a search for promising leverage points. Food Policy 98, 101915 (2021).
- Springmann, M. et al. Health and nutritional aspects of sustainable diet strategies and their association with environmental impacts: a global modelling analysis with country-level detail. Lancet Planet. Health 2, e451–e461 (2018).
- Springmann, M. et al. Options for keeping the food system within environmental limits. Nature 562, 519–525 (2018).
- Blakstad, M. M. et al. Life expectancy and agricultural environmental impacts in Addis Ababa can be improved through optimized plant and animal protein consumption. Nat. Food 2, 291–298 (2021).
- van Dooren, C. A review of the use of linear programming to optimize diets, nutritiously, economically and environmentally. Front. Nutr. 5, 48 (2018).
- Tuninetti, M., Ridolfi, L. & Laio, F. Ever-increasing agricultural land and water productivity: a global multi-crop analysis. Environ. Res. Lett. https://doi.org/10.1088/1748-9326/abacf8 (2020).
- Nyathi, M. K., Mabhaudhi, T., Van Halsema, G. E., Annandale, J. G. & Struik, P. C. Benchmarking nutritional water productivity of twenty vegetables - a review. Agric. Water Manag. 221, 248–259 (2019).
- Green, A., Nemecek, T., Smetana, S. & Mathys, A. Reconciling regionally-explicit nutritional needs with environmental protection by means of nutritional life cycle assessment. J. Clean. Prod. 312, 127696 (2021).
- Chen, C., Chaudhary, A. & Mathys, A. Nutritional and environmental losses embedded in global food waste. Resour. Conserv. Recycl. 160, 104912 (2020).
- Green, A., Nemecek, T., Chaudhary, A. & Mathys, A. Assessing nutritional, health, and environmental sustainability dimensions of agri-food production. Glob. Food Sec. 26, 100406 (2020).
- IPCC Guidelines for National Greenhouse Gas Inventories (IPCC, 2006).
- Blanford, D. Border and Related Measures in the Context of Adaptation and Mitigation to Climate Change. The State ofAgricultural Commodity Markets (SOCO): Background paper (FAO, 2018); https://www.fao.org/publications/card/es/c/CA2343EN/
- Springmann, M. & Freund, F. Options for reforming agricultural subsidies from health, climate, and economic perspectives. Nat. Commun. 13, 82 (2022).
- Springmann, M. et al. The healthiness and sustainability of national and global food based dietary guidelines: modelling study. Br. Med. J. 370, m2322 (2020).
- Schleussner, C., Ganti, G., Rogelj, J. & Gidden, M. J. An emission pathway classification reflecting the Paris Agreement climate objectives. Commun. Earth Environ. 3, 135 (2022).
- Diaz-Ruiz, R., Costa-Font, M., López-i-Gelats, F. & Gil, J. M. Food waste prevention along the food supply chain: a multi-actor approach to identify effective solutions. Resour. Conserv. Recycl. 149, 249–260 (2019).
- Reynolds, C. et al. Consumption-stage food waste reduction interventions – what works and how to design better interventions. Food Policy 83, 7–27 (2019).
- Stathers, T. et al. A scoping review of interventions for crop postharvest loss reduction in sub-Saharan Africa and south Asia. Nat. Sustain. 3, 821–835 (2020).
- Lipinski, B. et al. Reducing Food Loss and Waste Working Paper (World Resources Institute, 2013); http://www.worldresourcesreport.org
- Pradhan, P., Fischer, G., van Velthuizen, H., Reusser, D. E. & Kropp, J. P. Closing yield gaps: how sustainable can we be? PLoS ONE 10, e0129487 (2015).
- Stratton, A. E. et al. Mitigating sustainability tradeoffs as global fruit and vegetable systems expand to meet dietary recommendations. Environ. Res. Lett. 16, 055010 (2021).
- Schreinemachers, P., Simmons, E. B. & Wopereis, M. C. S. Tapping the economic and nutritional power of vegetables. Glob. Food Sec. 16, 36–45 (2018).
- Hawkes, C., Walton, S., Haddad, L. & Fanzo, J. 42 Policies and Actions to Orient Food Systems Towards Healthier Diets for All (Centre forFood Policy, City, Univ. London, 2020); https://researchcentres.city.ac.uk/food-policy#unit=publications
- Pingali, P. Agricultural policy and nutrition outcomes – getting beyond the preoccupation with staple grains. Food Secur. 7, 583–591 (2015).
- Arndt, C. et al. COVID-19 lockdowns, income distribution, and food security: an analysis for South Africa. Glob. Food Sec. 26, 100410 (2020).
- Hadjikakou, M. & Wiedmann, T. in Handbook on Growth and Sustainability (eds Victor, P. A. & Dolter, B.) 256–276 (Edward Elgar, 2017).
- Gerber, P. J. et al. Tackling Climate Change through Livestock—A Global Assessment of Emissions and Mitigation Opportunities (FAO, 2013); https://www.fao.org/publications/card/en/c/030a41a8-3e10-57d1-ae0c-86680a69ceea/
- Paul, B. K., Butterbach-Bahl, K., Notenbaert, A., Nduah Nderi, A. & Ericksen, P. Sustainable livestock development in low- and middle-income countries: shedding light on evidence-based solutions. Environ. Res. Lett. 16, 011001 (2021).
- Li, A. et al. Financial inclusion may limit sustainable development under economic globalization and climate change. Environ. Res. Lett. 16, 054049 (2021).
- Jackson, P. et al. Food as a commodity, human right or common. Nat. Food 2, 132–134 (2021).
- Fader, M. et al. Past and present biophysical redundancy of countries as a buffer to changes in food supply. Environ. Res. Lett. 11, 055008 (2016).
- Freund, F. & Springmann, M. Policy analysis indicates health-sensitive trade and subsidy reforms are needed in the UK to avoid adverse dietary health impacts post-Brexit. Nat. Food 2, 502–508 (2021).
- Bai, Z. et al. Food and feed trade has greatly impacted global land and nitrogen use efficiencies over 1961–2017. Nat. Food 2, 780–791 (2021).
- Hadjikakou, M., Ritchie, E. G., Watermeyer, K. E. & Bryan, B. A. Improving the assessment of food system sustainability. Lancet Planet. Health 3, e62–e63 (2019).
- Li, M. et al. Global food-miles account for nearly 20% of total food-systems emissions. Nat. Food 3, 445–453 (2022).
- Willett, W. et al. Food in the Anthropocene: the EAT–Lancet Commission on healthy diets from sustainable food systems. Lancet 393, 447–492 (2019).
- Anderson, J. R. & Birner, R.Fruits and vegetables in international agricultural research: a case of neglect? World Rev. Nutr. Diet. 121, 42–59 (2020).
- Nelson, G. et al. Income growth and climate change effects on global nutrition security to mid-century. Nat. Sustain. 1, 773–781 (2018).
- Delzeit, R. et al. Linking global CGE models with sectoral models to generate baseline scenarios: approaches, challenges, and opportunities. J. Glob. Econ. Anal. 5, 162–195 (2020).
- Geyik, O., Hadjikakou, M. & Bryan, B. A. Spatiotemporal trends in adequacy of dietary nutrient production and food sources. Glob. Food Sec. 24, 100355 (2020).
- Food Balances (-2013, old methodology and population) (FAOSTAT, 2017); https://www.fao.org/faostat/en/#data/FBSH
- Gustavsson, J., Cederberg, C., Sonesson, U., Van Otterdijk, R. & Meybeck, A. Global Food Losses and Food Waste – Extent, Causes and Prevention (FAO, 2011); https://www.fao.org/3/mb060e/mb060e00.htm
- Poore, J. & Nemecek, T. Reducing food’s environmental impacts through producers and consumers. Science 360, 987–992 (2018).
- Emissions - Agriculture (FAOSTAT, 2020); http://www.fao.org/faostat/en/#data/GT
- Forster, P. et al. in Climate Change 2021: The Physical Science Basis (eds Masson-Delmotte, V. et al.) 923–1054 (IPCC, Cambridge Univ. Press, 2021).
- Tubiello, F. N. et al. Greenhouse gas emissions from food systems: building the evidence base. Environ. Res. Lett. 16, 065007 (2021).
- Naylor, R. L. et al. A 20-year retrospective review of global aquaculture. Nature 591, 551–563 (2021).
- Tilman, D. & Clark, M. Global diets link environmental sustainability and human health. Nature 515, 518–522 (2014).
- Tacon, A. G. J., Hasan, M. R. & Metian, M. Demand and Supply of Feed Ingredients for Farmed Fish and Crustaceans: Trends and Prospects (FAO, 2011); http://www.fao.org/3/ba0002e/ba0002e.pdf
- Technical Conversion Factors for Agricultural Commodities (FAO, 2000); http://www.fao.org/economic/the-statistics-division-ess/methodology/methodology-systems/technical-conversion-factors-for-agricultural-commodities/en/
- Heffer, P., Gruère, A. & Roberts, T. Assessment of Fertilizer Use by Crop at the Global Level (International Fertilizer Association and International Plant Nutrition Institute, 2017); https://www.ifastat.org/plant-nutrition
- Bouwman, A. F., Van Der Hoek, K. W., Eickhout, B. & Soenario, I. Exploring changes in world ruminant production systems. Agric. Syst. 84, 121–153 (2005).
- Annex: Global Indicator Framework for the Sustainable Development Goals and Targets of the 2030 Agenda for Sustainable Development (United Nations, 2017).
- GAEZ v3.0: Global Agro-ecological Zones (IIASA/FAO, 2012); http://webarchive.iiasa.ac.at/Research/LUC/GAEZv3.0/
- Mottet, A. et al. Climate change mitigation and productivity gains in livestock supply chains: insights from regional case studies. Reg. Environ. Change 17, 129–141 (2017).
- Detailed Trade Matrix (FAOSTAT, 2019); http://www.fao.org/faostat/en/#data/TM
- Optimization (scipy.optimize) (SciPy, 2022); https://docs.scipy.org/doc/scipy/tutorial/optimize.html#id55
- Production (FAOSTAT, 2019); http://www.fao.org/faostat/en/#data
- Pouillot, R. & Delignette-Muller, M.-L. Evaluating variability and uncertainty in microbial quantitative risk assessment using two R packages. Int. J. Food Microbiol. 142, 330–340 (2010).
- Huppmann, D. et al. IAMC 1.5 °C scenario explorer and data hosted by IIASA. Zenodo https://doi.org/10.5281/zenodo.3363345 (2019).
- Smith, M. A., Cain, M. & Allen, M. R. Further improvement of warming-equivalent emissions calculation. npj Clim. Atmos. Sci. 4, 19 (2021).
- Tennekes, M. tmap: thematic maps in R. J. Stat. Softw. 84, 1–39 (2018).
- Data: World Bank Country and Lending Groups (World Bank, 2020); https://datahelpdesk.worldbank.org/knowledgebase/articles/906519-world-bank-country-and-lending-groups
Download references
Acknowledgements
The financial support for this research, from the German Research Foundation (DFG) through the Sustainable Food Systems Research Training Group (RTG 2654) and Deakin University, Australia (DUPR-STRATEGIC – 0000018831), both received by O.G., is gratefully acknowledged. We also thank P. Pradhan for kindly sharing preliminary R code used to estimate potential yields via the GAEZ model and G. Ganti for providing data used to determine non-CO2 emission pathways compatible with the Paris Agreement
Author information
Authors and Affiliations
- Centre for Integrative Ecology, School of Life and Environmental Sciences, Deakin University, Melbourne, Victoria, Australia
- Sustainable Food Systems Research Training Group (RTG 2654), University of Göttingen, Göttingen, Germany
Özge Geyik, Michalis Hadjikakou & Brett A. Bryan
Özge Geyik
Contributions
O.G., M.H. and B.A.B. designed the study. O.G. collated data and performed the analyses with the help of M.H., who also calculated crop yield gap estimates and climate boundaries. All authors discussed the methods and results, and helped shape the research, analysis and interpretation. O.G. took the lead in writing the manuscript with substantial contributions from all authors.
Corresponding author
Correspondence to Özge Geyik.
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Food thanks Laixiang Sun 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 Table 1 List of countriesFull size tableExtended Data Table 2 Relative proportion of N fertiliser use by crop groupFull size tableExtended Data Table 3 List of food productsFull size tableExtended Data Table 4 Technical conversion factors for oilseed processingFull size tableExtended Data Table 5 Livestock GHG mitigation potentialFull size tableExtended Data Table 6 Proportional consumption of crop feed by livestock groupFull size tableExtended Data Table 7 Average aquafeed use per unit productionFull size tableExtended Data Table 8 Technical conversion factorsFull size tableExtended Data Table 9 Uncertainty range of emission factorsFull size table
Extended Data Fig. 1 Breakdown of total GHG emissions by individual climate pollutant.
a) CH4 emissions results by scenario and income level. Results are presented in megatons CH4 yr-1. b) N2O emissions results by scenario and income level. Results are presented in megatons N2O yr-1. Bars show the total emissions for 128 countries, that is, n = 128, based on default emissions factors (corresponding to mode as the measure of center. Error bars show the 25th and 75th percentiles (see Uncertainty estimates).
Supplementary information
Supplementary Information
Supplementary Fig. 1 and Supplementary Discussion.
Reporting Summary
Supplementary Table
Supplementary tables.
Source data
Source Data Fig. 1
Statistical source data.
Source Data Figs. 2 and 3
Statistical source data.
Source Data Fig. 4
Statistical source data.
Source Data Fig. 5
Statistical source data.
Source Data Extended Data Fig. 1
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
Geyik, Ö., Hadjikakou, M. & Bryan, B.A. Climate-friendly and nutrition-sensitive interventions can close the global dietary nutrient gap while reducing GHG emissions. Nat Food 4, 61–73 (2023). https://doi.org/10.1038/s43016-022-00648-y