Main

Human actions are rapidly transforming the planet, driving losses of nature at an unprecedented rate that negatively impacts societies and economies, from accelerating climate change to increasing zoonotic pandemic risk1,2. Recognizing the accelerating severity of the environmental crisis, the global community committed to Sustainable Development Goals and the Paris Agreement on climate change in 2015. In 2022, the UN Convention on Biological Diversity (CBD) will adopt new targets for conserving, restoring and sustainably managing multiple dimensions of biodiversity, including nature’s contributions to people (NCP)3. Collectively, these three policy frameworks will shape the sustainable development agenda for the next decade. All three depend heavily on safeguarding natural assets1,4, the living components of our lands and waters. For instance, restoring and ending conversion and degradation of forests, wetlands and peatlands could sequester 9 Gt CO2 per year by 2050 (ref. 5). While ambitious new targets to protect species and ecosystems have been proposed, including ‘half Earth’ (conserving half the Earth’s area for nature)6 and ‘30 by 30’ (30% protected by 2030) (ref. 7), these targets have been criticized for insufficiently accounting for the needs of people, including many Indigenous and local communities8. It is therefore essential to demonstrate how nature conservation contributes to human wellbeing. Yet despite the urgency of safeguarding natural assets around the world, we still have limited understanding of the spatial extent of ecosystems providing essential benefits to humanity9.

In this Article, leveraging recent advances in scientific understanding, data availability (Supplementary Tables 1 and 2) and computational power, we undertake a global analysis of 14 NCP (Extended Data Fig. 1), the most comprehensive set mapped globally so far10,11. Most of these NCP are considered ecosystem services, which, according to the IPBES Conceptual Framework3,12, are embedded within the larger, more inclusive concept of NCP, and therefore we use the latter term throughout to better align with international frameworks and policy13. Twelve of the NCP included in our analysis deliver primarily local benefits (though some subsequently enter global supply chains), including contributions to the provision of food, energy and raw materials; the regulation of water quality and disaster risk; and recreational activities (Fig. 1a). We prioritize these 12 ‘local’ NCP at the country level to identify areas needed to maintain their consistent provision within each country. In contrast, we prioritize at the global scale for two NCP related to climate regulation (terrestrial ecosystem carbon storage and vegetation-regulated atmospheric moisture recycling), whose benefits accrue at continental scales or globally.

Fig. 1: Critical natural assets, defined as the natural and semi-natural terrestrial and aquatic ecosystems required to maintain 12 of nature’s ‘local’ contributions to people (local NCP) on land (green) and in the ocean (blue).
figure 1

a, The 12 local NCP analysed (that is, not including global NCP, shown in Supplementary Fig. 4). b, The NCP accumulation curve, reflecting the total area required to maintain target levels of all NCP in every country, with dotted lines denoting the area of critical natural assets (90% of NCP in 30% of land area and 24% of EEZ area). Areas selected by optimization within each country are aggregated across all countries to create a single global accumulation curve; for area requirements in individual countries, see Supplementary Data 1. c, Map of critical natural assets, with darker shades connoting critical natural assets that are associated with higher levels of aggregated NCP. Grey areas show the extent of remaining natural assets not designated ‘critical’ but included in this analysis; white areas (cropland, urban and bare areas, ice and snow, and ocean areas outside the EEZ) were excluded from the optimization.

Through multi-objective spatial optimization at a resolution of 2 km, we map the location of the planet’s critical natural assets, defined as the natural and semi-natural ecosystems providing 90% of current levels of each NCP (that is, locations required for all NCP to meet or exceed a 90% target). Beyond this target there are diminishing returns, with disproportionately more natural area required to reach incrementally higher magnitudes of provision of NCP (Fig. 1b). We use the term ‘critical’ to convey that losing these natural assets would lead to disproportionately large losses in NCP. This usage is related to, though not the same as, the term ‘critical natural capital’, which has been formally defined in the ecological economics literature as ecosystems providing important functions that cannot be substituted14. However, as we do not explore substitutes to natural capital found in alternative forms of capital (that is, manufactured capital), we use the term ‘natural assets’ instead to avoid confusion. The assets we focus on are natural and semi-natural ecosystems (Supplementary Table 3), excluding developed lands (croplands and urban areas) and unvegetated areas, to provide insights for conservation priorities relevant to the CBD; global priorities for restoration15 or for management of developed areas16 are beyond the scope of this first effort to map critical natural assets. Our analysis reveals three key findings about critical natural assets: (1) their extent and location; (2) the number of people benefitting from, and living within, these critical areas; and (3) the degree of overlap between critical natural assets delivering local benefits and those delivering global benefits and between critical natural assets and other global priorities for the CBD (biodiversity and cultural diversity).

Results and discussion

Extent and location of critical natural assets

Critical natural assets providing the 12 local NCP (Fig. 1a) occupy only 30% (41 million km2) of total land area (excluding Antarctica) and 24% (34 million km2) of marine Exclusive Economic Zones (EEZs), reflecting the steep slope of the aggregate NCP accumulation curve (Fig. 1b). Despite this modest proportion of global land area, the shares of countries’ land areas that are designated as critical can vary substantially. The 20 largest countries require only 24% of their land area, on average, to maintain 90% of current levels of NCP, while smaller countries (10,000 to 1.5 million km2) require on average 40% of their land area (Supplementary Data 1). This high variability in the NCP–area relationship is primarily driven by the proportion of countries’ land areas made up by natural assets (that is, excluding barren, ice and snow, and developed lands), but even when this is accounted for, there are outliers (Extended Data Fig. 2). Outliers may be due to spatial patterns in human population density (for example, countries with dense population centres and vast expanses with few people, such as Canada and Russia, require far less area to achieve NCP targets) or large ecosystem heterogeneity (if greater ecosystem diversity yields higher levels of diverse NCP in a smaller proportion of area, which may explain patterns in Chile and Australia).

The highest-value critical natural assets (the locations delivering the highest magnitudes of NCP in the smallest area, denoted by the darkest blue or green shades in Fig. 1c) often coincide with diverse, relatively intact natural areas near or upstream from large numbers of people. Many of these high-value areas coincide with areas of greatest spatial congruence among multiple NCP (Extended Data Fig. 3). Spatially correlated pairs of local NCP (Supplementary Table 4) include those related to water (flood risk reduction with nitrogen retention and nitrogen with sediment retention); forest products (timber and fuelwood); and those occurring closer to human-modified habitats (pollination with nature access and with nitrogen retention). Coastal risk reduction, forage production for grazing, and riverine fish harvest are the most spatially distinct from other local NCP. In the marine realm, there is substantial overlap of fisheries with coastal risk reduction and reef tourism (though not between the latter two, which each have much smaller critical areas than exist for fisheries).

Number of people benefitting from critical natural assets

We estimate that ~87% of the world’s current population, 6.4 billion people, benefit directly from at least one of the 12 local NCP provided by critical natural assets, while only 16% live on the lands providing these benefits (and they may also benefit; Fig. 2a). To quantify the number of beneficiaries of critical natural assets, we spatially delineate their benefitting areas (which varies on the basis of NCP: for example, areas downstream, within the floodplain, in low-lying areas near the coast, or accessible by a short travel). While our optimization selects for the provision of 90% of the current value of each NCP, it is not guaranteed that 90% of the world’s population would benefit (since it does not include considerations for redundancy in adjacent pixels and therefore many of the areas selected benefit the same populations), so it is notable that an estimated 87% do. This estimate of ‘local’ beneficiaries probably underestimates the total number of people benefitting because it includes only NCP for which beneficiaries can be spatially delineated to avoid double-counting, yet it is striking that the vast majority, 6.1 billion people, live within 1 h travel (by road, rail, boat or foot, taking the fastest path17) of critical natural assets, and more than half of the world’s population lives downstream of these areas (Fig. 2b). Material NCP are often delivered locally, but many also enter global supply chains, making it difficult to delineate beneficiaries spatially for these NCP. However, past studies have calculated that globally more than 54 million people benefit directly from the timber industry18, 157 million from riverine fisheries19, 565 million from marine fisheries20 and 1.3 billion from livestock grazing21, and across the tropics alone 2.7 billion are estimated to be dependent on nature for one or more basic needs22.

Fig. 2: People benefitting from and living on critical natural assets (CNA).
figure 2

a,b, ‘Local’ beneficiaries were calculated through the intersection of areas benefitting from different NCP, to avoid double-counting people in areas of overlap; only those NCP for which beneficiaries could be spatially delineated were included (that is, not material NCP that enter global supply chains: fisheries, timber, livestock or crop pollination). Bars show percentages of total population globally and for large and small countries (a) or the percentage of relevant population globally (b). Numbers inset in bars show millions of people making up that percentage. Numbers to the right of bars in b show total relevant population (in millions of people, equivalent to total global population from Landscan 2017 for population within 1 h travel or downstream, but limited to the total population living within 10 km of floodplains or along coastlines <10 m above mean sea level for floodplain and coastal populations protected, respectively, and to rural poor populations for fuelwood).

Nearly all countries have a large percentage (>80%) of their populations benefitting from critical natural assets, but small countries have much larger proportions of their populations living within the footprint of critical natural assets than do large countries (Fig. 2a and Supplementary Data 2). When people live in these areas, and especially when current levels of use of natural assets are not sustainable, regulations or incentives may be needed to maintain the benefits these assets provide. While protected areas are an important conservation strategy, they represent only 15% of the critical natural assets for local NCP (Supplementary Table 5); additional areas should not necessarily be protected using designations that restrict human access and use, or they could cease to provide some of the diverse values that make them so critical23. Other area-based conservation measures, such as those based on Indigenous and local communities’ governance systems, Payments for Ecosystem Services programmes, and sustainable use of land- and seascapes, can all contribute to maintaining critical flows of NCP in natural and semi-natural ecosystems24.

Overlaps between local and global priorities

Unlike the 12 local NCP prioritized here at the national scale, certain benefits of natural assets accrue continentally or even globally. We therefore optimize two additional NCP at a global scale: vulnerable terrestrial ecosystem carbon storage (that is, the amount of total ecosystem carbon lost in a typical disturbance event25, hereafter ‘ecosystem carbon’) and vegetation-regulated atmospheric moisture recycling (the supply of atmospheric moisture and precipitation sustained by plant life26, hereafter ‘moisture recycling’). Over 80% of the natural asset locations identified as critical for the 12 local NCP are also critical for the two global NCP (Fig. 3). The spatial overlap between critical natural assets for local and global NCP accounts for 24% of land area, with an additional 14% of land area critical for global NCP that is not considered critical for local NCP (Extended Data Fig. 4). Together, critical natural assets for securing both local and global NCP require 44% of total global land area. When each NCP is optimized individually (carbon and moisture NCP at the global scale; the other 12 at the country scale), the overlap between carbon or moisture NCP and the other NCP exceeds 50% for all terrestrial (and freshwater) NCP except coastal risk reduction (which overlaps only 36% with ecosystem carbon, 5% with moisture recycling; Supplementary Table 4).

Fig. 3: Spatial overlaps between critical natural assets for local and global NCP.
figure 3

Red and teal denote where critical natural assets for global NCP (providing 90% of ecosystem carbon and moisture recycling globally) or for local NCP (providing 90% of the 12 NCP listed in Fig. 1), respectively, but not both, occur; gold shows areas where the two overlap (24% of the total area). Together, local and global critical natural assets account for 44% of total global land area (excluding Antarctica). Grey areas show natural assets not defined as ‘critical’ by this analysis, though still providing some values to certain populations. White areas were excluded from the optimization.

Synergies can also be found between NCP and biodiversity and cultural diversity. Critical natural assets for local NCP at national levels overlap with part or all of the area of habitat (AOH, mapped on the basis of species range maps, habitat preferences and elevation27) for 60% of 28,177 terrestrial vertebrates (Supplementary Data 3). Birds (73%) and mammals (66%) are better represented than reptiles and amphibians (44%). However, these critical natural assets represent only 34% of the area for endemic vertebrate species (with 100% of their AOH located within a given country; Supplementary Data 3) and 16% of the area for all vertebrates if using a more conservative representation target framework based on the IUCN Red List criteria (though, notably, achieving Red List representation targets is impossible for 24% of species without restoration or other expansion of existing AOH; Supplementary Data 4). Cultural diversity (proxied by linguistic diversity) has far higher overlaps with critical natural assets than does biodiversity; these areas intersect 96% of global Indigenous and non-migrant languages28 (Supplementary Data 5). The degree to which languages are represented in association with critical natural assets is consistent across most countries, even at the high end of language diversity (countries containing >100 Indigenous and non-migrant languages, such as Indonesia, Nigeria and India). This high correspondence provides further support for the importance of safeguarding rights to access critical natural assets, especially for Indigenous cultures that benefit from and help maintain them. Despite the larger land area required for maintaining the global NCP compared with local NCP, global NCP priority areas overlap with slightly fewer languages (92%) and with only 2% more species (60% of species AOH), although a substantially greater overlap is seen with global NCP if Red List criteria are considered (36% compared with 16% for local NCP; Supplementary Data 4). These results provide different insights than previous efforts at smaller scales, particularly a similar exercise in Europe that found less overlap with priority areas for biodiversity and NCP29. However, the 40% of all vertebrate species whose habitats did not overlap with critical natural assets could drive very different patterns if biodiversity were included in the optimization.

Although these 14 NCP are not comprehensive of the myriad ways that nature benefits and is valued by people23, they capture, spatially and thematically, many elements explicitly mentioned in the First Draft of the CBD’s post-2020 Global Biodiversity Framework13: food security, water security, protection from hazards and extreme events, livelihoods and access to green and blue spaces. Our emphasis here is to highlight the contributions of natural and semi-natural ecosystems to human wellbeing, specifically contributions that are often overlooked in mainstream conservation and development policies around the world. For example, considerations for global food security often include only crop production rather than nature’s contributions to it via pollination or vegetation-mediated precipitation, or livestock production without partitioning out the contribution of grasslands from more intensified feed production.

Gaps and next steps

Our synthesis of these 14 NCP represents a substantial advance beyond other global prioritizations that include NCP limited to ecosystem carbon stocks, fresh water and marine fisheries30,31,32, though still falls short of including all important contributions of nature such as its relational values33. Despite the omission of many NCP that were not able to be mapped, further analyses indicate that results are fairly robust to inclusion of additional NCP. Dropping one of the 12 local NCP at a time results in <1–3% change in the total global land area required to maintain 90% of current levels of these NCP (Supplementary Table 6) and a high degree of spatial agreement. In fact, 62% of the total area on land is shared by all optimization solutions, and 97% of the area is included in 11 of the 12 solutions; similar values are found across most countries (Supplementary Data 1). Nevertheless, this same multi-NCP optimization approach could accommodate additional NCP as spatial data become available at sufficient resolution and appropriate scales.

There also is uncertainty in the identification of critical natural assets related to model error in the individual NCP that we were able to include. We acknowledge that NCP models, like all models, have errors, and that consistent global-scale modelling will miss details important for certain specific locations. Validation of NCP is particularly difficult given that there are no direct measurements for many NCP with assessment reliant on remotely sensed proxies. We utilize the best available global modelling approaches and data, most of which have been validated in at least some locations19,25,34,35,36,37,38,39,40,41. Where uncertainty existed about what distance was most appropriate to model the delivery of NCP (for example, how far to model people downstream or how far people might travel to natural assets), we performed further sensitivity analysis and confirmed that the estimated land area of critical natural assets is robust to the distance chosen (impacting results by <5%; Supplementary Table 6). As availability of global models for many of these NCP increases, future work should move towards ensemble modelling approaches, which have been shown to increase accuracy and reduce uncertainty compared with individual model outputs42,43.

Data and modelling gaps prevented a broader exploration of issues relevant to the ecological supply side of NCP. Although results presented here suggest that nationally prioritized areas for local NCP can deliver on global priorities in many regions, they also highlight a need for integrated modelling to represent interactions between different NCP. For example, atmospheric moisture evapotranspired by Amazonian forests falls as rain in other parts of South America, supporting ecosystems that provide food, fuel and other benefits26. Further work is needed to move beyond the spatial overlaps explored here towards understanding functional inter-dependencies between NCP. We also acknowledge that urban and cropland systems were omitted from this analysis owing to data and modelling limitations that would fail to adequately capture the NCP supported by different land use types and land management practices within those systems. Likewise, arctic and desert ecosystems, owing to sparse vegetation and low human population densities, are not well represented in our NCP models and yet are very important to the people who live in and depend on them. As data and modelling gaps are filled, future assessment of critical natural assets should expand to recognize unique contributions of currently undervalued ecosystems and should include possible gains from restoring and sustainably managing human-dominated systems15,16 to consider how these different conservation strategies can complement one another.

There are also, perhaps even more pronounced, data and modelling gaps to fill on the social side of NCP. In the NCP modelled here, we represent realized benefits of natural assets—weighted by beneficiary population when feasible—but this understates the range of ways in which natural assets directly and indirectly contribute to people’s wellbeing. Limited socio-economic data and lack of reliable models linking NCP to wellbeing indicators preclude more precise valuation of most NCP at the global scale. Additional insight could be gained from mapping critical natural assets that support the most vulnerable or dependent22 people, including Indigenous peoples whose livelihoods and cultural identities directly depend on nature (and indeed overlap substantially with critical natural assets, on the basis of our estimates of Indigenous language diversity on these lands; Supplementary Data 5), and the poor who may lack access to anthropogenic substitutes for NCP (see also philosophical considerations in Supplementary Note 1). Recent progress in linking ecological modelling with integrated assessment modelling and general equilibrium economic modelling44 shows great promise for assessing the benefits of critical natural assets to society and the global economy. Such efforts could also reveal telecoupling of critical natural assets arising from transboundary flows between countries such as via international trade45.

Finally, further work is needed to investigate whether critical natural assets are necessary or sufficient for meeting humanity’s needs, by considering the availability of substitutes and what constitutes ‘enough’. Though motivated by the ecological economic concept of ‘critical natural capital’, we were not able to capture the degree to which anthropogenic assets could replace natural assets. However, given the high correspondence of many NCP (Supplementary Table 4), it seems unlikely that anthropogenic assets could substitute for all benefits provided by natural assets in a particular area (for example, green infrastructure delivers many co-benefits beyond the single benefit built infrastructure is designed to deliver46,47). It seems more likely that more than 44% of the planet may be required to secure the 14 NCP mapped here (let alone the other diverse values of nature we were not able to map), most importantly because it is unclear how much of the current need for nature is already unmet. In many parts of the world, natural ecosystems are already degraded or converted, so maintaining 90% of current levels of NCP may be far too little (for example, places prone to catastrophic flooding due to habitat conversion48, where grazing lands have been desertified49 and where fish populations have crashed50). Furthermore, what is critical now may change in the future owing to climate change, population growth or change in technology or consumption patterns1,11. Other areas may not directly benefit people but may be critical to avoid ecological tipping points and collapse of NCP in affected ecosystems51,52. We therefore propose the analysis conducted here as a first attempt to define a minimum set of critical natural assets, and suggest that the overall approach provides a useful framework for exploring such issues with future scenario modelling to develop more resilient conservation for nature and people.

We acknowledge that our approach is strongly anthropocentric (NCP, and all the concepts included in them, are by definition anthropocentric3). As such, it is not intended to capture the intrinsic values of nature, or the value of the ecosystems or species providing the NCP highlighted here to other non-human organisms. Our focus and findings should not be interpreted as dismissing those values, and further work could explore natural assets providing important contributions to biodiversity that are not captured by species maps alone. For example, regulating contributions, including water quality regulation, natural hazards resilience, pollination and atmospheric moisture recycling, maintain the conditions under which current biodiversity thrives. Delineating species or high-biodiversity areas as the ‘beneficiaries’ for many of these contributions may be an important step towards reflecting nature’s contributions to nature.

Conclusion

Identifying critical natural assets could enable national and global leaders to prioritize the conservation of a wide range of NCP. We find it encouraging that securing 90% of the NCP mapped here is feasible with an area comparable to other proposed conservation targets6,7,32,53. Global analyses such as this can set a broader context for local decisions including understanding of distant connections that extend outside a country’s borders, provide rapid information for global actors such as CBD signatories and non-governmental organizations working on conservation priorities across many countries, and supplement gaps in local information while it is still being generated54. National- and local-level use of spatial data is timely and relevant for countries seeking to access scarce international Sustainable Development Goal-related finance, as it can improve the ecological, social and economic aspects of project readiness55,56. However, we emphasize the value of the approach developed here more than the maps or data. This approach can be adapted and refined at the scales at which policy implementation occurs, with the best available data and complemented by input from local experts and diverse stakeholders, to improve accuracy and public legitimacy57,58 and to ensure that human rights and diverse human relationships with nature are safeguarded. Moreover, creation and use of spatially explicit information allows for a focus on ecosystem quality over quantity, helping to avoid potentially perverse outcomes of area targets for conservation59. This approach for identifying critical natural assets is a vital step forward in empowering actors at all levels to make decisions that benefit both nature and people.

Methods

Modelling NCP

The 14 NCP in this analysis (Extended Data Fig. 1) were chosen to span development and climate goals, and to be mappable with spatially explicit data representing the period 2000–2020. We use European Space Agency (ESA) 2015, for land cover at 10 arcsec (~300 m at the equator) resolution, Landscan 2017 for population60 at 30 arcsec (~1 km) resolution (which were the most current data available at the time we began our analysis). We focus on ‘nature’s contributions’ to key benefits of interest (for example, security in food, water, hazards, materials and culture), meaning we partition out the role of natural and semi-natural ecosystems in producing those benefits. For food security, we include the contributions of pollinator habitat sufficiency to pollination-dependent crop production, vegetation-mediated atmospheric moisture recycling to crop and livestock production (included as a global NCP that also links to climate security), grassland fodder production to livestock production, and wild riverine and marine fisheries. For water security, we include the contributions of water quality regulation, via sediment retention and nutrient retention, but not water yield since the role of ecosystems in determining the quantity of water is minimal (other than by evapotranspiration, which is already captured in the vegetation-mediated moisture recycling, and regulation of timing of flows which is captured in flood risk reduction). For security of protection from natural hazards, we include flood risk reduction and coastal risk reduction. For materials, we include timber production, fuelwood production and access to nature (which could be used for gathering, and also links to culture). For cultural benefits we include coral reef tourism (as the only globally mapped form of marine-based tourism) and access to nature again (which in addition to gathering also captures recreation or other uses of nearby greenspace). Finally, for climate security, in addition to moisture recycling mentioned above, we include total ecosystem carbon storage (as a global NCP). Below we briefly summarize the models used to map these local NCP (Supplementary Table 1) and global NCP (Supplementary Table 2); full information on each model is available in Supplementary Note 1.

Local NCP

  1. 1.

    Nitrogen retention to regulate water quality for downstream populations is modelled at native ESA resolution (at 10 arcsec or 300 m) using the InVEST61 Nutrient Delivery Ratio model, which is based on fertilizer application, precipitation, topography and the retention capacity of vegetation, and has been previously used in global applications11. The people benefitting from nitrogen retention are those who would otherwise be exposed to nitrogen contamination in their drinking water. In this analysis, the number of people downstream was calculated for every pixel of habitat, to provide a sense of which habitat potentially benefits the most people. Ideally, to map realized nitrogen retention, we would be able to convert biophysical service production into a measure of change in wellbeing, whether monetary, in health terms or otherwise. However, the state of the science and data available globally precludes this for most services, so our proxy was the number of people downstream who could potentially benefit from that retention. NCP for nitrogen retention is expressed as nitrogen retention on natural and semi-natural pixels multiplied by the number of people downstream of those pixels (for more details, see Supplementary Note 1 Section 1).

  2. 2.

    Sediment retention to regulate water quality for downstream populations is modelled at native ESA resolution (10 arcsec or 300 m) by adapting the InVEST Sediment Delivery Ratio model, which maps overland sediment generation and delivery to the stream using the Revised Universal Soil Loss Equation (RUSLE) based on climate, soil texture, topography and land cover, and a conductivity index based on the upslope and downslope areas of each pixel. Ideally sediment retention would be delineated for reservoirs, irrigation canals or other water delivery infrastructure that is most impacted by sedimentation, but lacking a comprehensive global dataset identifying all such infrastructure, we again use the proxy of number of people downstream (as described for nitrogen retention, above). NCP for sediment retention is expressed as sediment retention on natural and semi-natural pixels multiplied by the number of people downstream of those pixels (Supplementary Note 1 Section 2).

  3. 3.

    Pollinator habitat sufficiency for pollination-dependent crop production is modelled at native ESA resolution (10 arcsec or 300 m) with a simplified version of InVEST, improving upon previous global mapping of the potential contribution of wild pollinators to nutrition production by mapping the value from the farm field back to habitat11. This contribution is calculated as the sufficiency of habitat surrounding farmland multiplied by the pollination dependency of crops (the proportion of yields that are pollination dependent, based on the crop mixes grown there). While this does not capture other important features that affect pollinators, such as presence of specific species important for nesting or floral resources or the intensity or frequency of pesticide applications, it provides a proxy of the contribution of nearby natural habitat to pollination. This NCP is expressed in terms of the average equivalent number of people fed by pollination-dependent crop production, attributed to nearby ecosystems based on the area of pollinator habitat within pollinator flight distance of crops (Supplementary Note 1 Section 3).

  4. 4.

    Fodder production for livestock is modelled using version 3 of Co$ting Nature62 at a resolution of 5 arcmin (~10 km), and then masked to relevant ESA natural and semi-natural habitats (at 10 arcsec or 300 m resolution; Supplementary Table 3). Supply of fodder is calculated from remotely sensed dry matter productivity for the non-cropland cover fraction, and demand is estimated by the headcount of livestock in a grid cell multiplied by the average biomass intake requirements per animal. NCP for fodder production for livestock is expressed in terms of an index (0–1), re-scaled from the realized service, which is reported as the smaller of the supply or demand (if consumption exceeds productivity, the gap is assumed to be met with feed). The best available global inputs for dry matter productivity, livestock headcount, cropland and land cover are used as input (Supplementary Note 1 Section 4).

  5. 5.

    Timber production includes commercial (for example, for trade or export) and domestic (for example, for local use) timber, modelled using version 3 of Co$ting Nature as two spatially mutually exclusive layers, because they represent two different sets of beneficiaries. This model outputs data at a resolution of 5 arcmin (~10 km), and is then masked to relevant ESA natural and semi-natural habitats (at 10 arcsec or 300 m resolution; Supplementary Table 3). NCP for timber production is expressed as an index (0–1) based on remotely sensed forest productivity and accessibility for harvest. Total potential sustainable supply of timber is estimated from the best available global aboveground carbon stock map multiplied by fractional tree cover for rural areas only. The sustainable harvest is calculated as the reciprocal of the number of years taken to develop the stock at the annual sequestration rate, according to dry matter productivity data. Demand is calculated differently for commercial versus domestic timber based on different assumptions of accessibility. Commercial accessibility is defined as within 6 h travel time of a population centre of >50,000 people and on slope gradients <70%. Domestic accessibility is defined as areas inaccessible for commercial harvest, and harvest rates are based on per capita consumption multiplied by population within 10 km (Supplementary Note 1 Section 5).

  6. 6.

    Fuelwood production is calculated as a byproduct of the timber model from version 3 of Co$ting Nature. NCP for fuelwood production, like timber, is represented as an index (0–1) based on forest productivity and accessibility for harvest, but in this case specifically by rural people. Fuelwood can overlap spatially with domestic and commercial timber use, given that domestic and commercial timber harvest will not consume all sustainably available woody biomass in all places, owing to the slope gradient limit and/or in places where demand is less than supply, and fuelwood is often a byproduct of timber harvest. This model outputs data at a resolution of 5 arcmin (~10 km), and is then masked to relevant ESA natural and semi-natural habitats (at 10 arcsec or 300 m resolution; Supplementary Table 3) (Supplementary Note 1 Section 6).

  7. 7.

    Flood regulation is modelled using version 2 of the hydrologic model WaterWorld35. To map nature’s influence on flood risk reduction, we identify the upstream places where canopies, wetlands and soils (green storage) retain and slowly release rainfall, to the benefit of downstream communities. NCP for flood regulation is expressed as an index (0–1) of ‘green’ water storage (based on wetland storage, canopy storage and soil storage) multiplied by the number of people downstream. This model outputs data at a resolution of 5 arcmin (~10 km), and is then masked to ESA natural and semi-natural habitats (at 10 arcsec or 300 m resolution) (Supplementary Note 1 Section 7).

  8. 8.

    Access to nature is used as a proxy for numerous direct and indirect benefits to people, such as recreation, hunting and gathering, aesthetics, mental and physical health, and other cultural values that depend on the ability of people to access nature. This proxy NCP is expressed as the number of urban and rural63 people within 1 h (or 6 h, for sensitivity analysis) travel of natural and semi-natural habitat, taking the least-cost path (by foot, road, rail or boat) across a friction layer17 at a resolution of 2 km (equal-area projection) (Supplementary Note 1 Section 8).

  9. 9.

    Riverine fish catch is based on spatial disaggregation of nationally reported catch for 2007–2014 (ref. 19) and updated to include catch estimated by household consumption surveys in 32 countries with severe underreporting64. Catches from large lakes were excluded. To spatially disaggregate the global catch of 13.3 × 106 tonnes within country borders, a multiple linear regression model of total fish catch in river basins compiled from the literature was fitted with three predictor variables: population density, river discharge and percentage wetland cover (n = 40, R2adj = 0.69). NCP for riverine fish harvest is represented as metric tonnes of fish caught per km2 of land area per year, spatially allocated to the locations of the harvest. This model outputs data at a resolution of 5 arcmin (~10 km), and is then masked to ESA natural and semi-natural habitats (at 10 arcsec or 300 m resolution) (Supplementary Note 1 Section 9).

  10. 10.

    Marine fish catch is based on data from the Sea Around Us, Global Fishing Watch, and other sources to map fish catch for 2010–2014 within 30 min (~55 km) grid cells across the ocean65,66. NCP for marine fish harvest is represented as metric tonnes of fish caught per km2 per year, spatially disaggregated to the locations of the catch (Supplementary Note 1 Section 10).

  11. 11.

    Coral reef tourism is taken from the Mapping Ocean Wealth dataset34, which reports NCP for coral-reef-associated tourism as dollars of tourism expenditure (in deciles 1–10) at 15 arcsec (~500 m) resolution. National expenditure data are spatially distributed on the basis of three independent sources: hotel rooms from the commercial Global Accommodation Reference Database, dive shops and dive sites from Diveboard, and user-generated photos from the image-sharing website Flickr (Supplementary Note 1 Section 11).

  12. 12.

    Coastal risk reduction is modelled with InVEST for terrestrial and coastal/off-shore habitats67,68,69,70 at native ESA resolution (10 arcsec or 300 m) updating previous global modelling11 through the inclusion of new data and projecting the value from the shoreline back to the protective habitat (which may be off-shore). Coastal risk reduction depends on the physical exposure to coastal hazards (based on wind, waves, sea level rise, geomorphology and bathymetry), with and without natural habitat to attenuate storm surge, and the people exposed. NCP for coastal risk reduction is expressed as a unitless index of the coastal risk reduced by habitat multiplied by the number of people within the protective distance of the habitat (Supplementary Note 1 Section 12).

Global NCP

  1. 1.

    Vulnerable terrestrial ecosystem carbon storage is mapped at a resolution of 1 arcsec (~30 m), as the aboveground and belowground ecosystem carbon lost in a ‘typical’ disturbance event, rather than the total stock25. This includes terrestrial and coastal (mangrove, salt marsh and seagrass) ecosystem carbon pools (aboveground, belowground and soils), based on what carbon is likely to be released if the ecosystem were converted (Supplementary Note 1 Section 13).

  2. 2.

    Atmospheric moisture recycling is the process of water arising from the surface of the Earth as evaporation, flowing through the atmosphere as water vapour, and returning to the surface of the Earth as precipitation. Sources of evaporation include canopy interception, soil interception, soil evaporation, vegetation transpiration and open water evaporation. We employ an Eulerian moisture tracking model, WAM-2 layers26, to quantify the source of moisture, where it travels through the atmosphere and where it falls out downwind across a 1.5° grid (~166 km). NCP of moisture recycling, which is to say the moisture associated with vegetation, is expressed as the volume of water evapotranspired that falls on all rainfed productive land (cropland, rangelands and working forests71) (Supplementary Note 1 Section 14).

Attribution of value to natural assets

The first step in identifying critical natural assets, before the optimization can select the highest value areas across all NCP, is to attribute the magnitude of benefits and, where possible, the number of beneficiaries, to the ecosystems providing those benefits for each individual NCP on the basis of its unique characteristics (for example, attributing the value of pollination occurring on croplands to the nearby habitat supplying the pollinators, or the coastal risk reduced and number of people protected along the coastline to off-shore as well as on-shore habitats). We define natural assets as natural and semi-natural terrestrial ecosystems (including semi-natural lands like rangelands and production forests, but excluding cropland, urban areas, bare areas and permanent snow and ice; Supplementary Table 3) and inland and marine waters. Model outputs for pollination and coastal risk reduction are mapped back to habitat on the basis of pollinator flight distance (Supplementary Note 1 Section 3) and protective distance of coastal habitat (Supplementary Note 1 Section 12), respectively. For sediment and nitrogen retention, the count of people downstream of each habitat pixel was summed according to a hydrologic flow accumulation (Supplementary Note 1 Section 1), and for nature access, the count of people was calculated for each pixels by summing the population pixels within 1 h travel time according to a friction surface (Supplementary Note 1 Section 8). All other model outputs are coarser than ESA resolution and are masked to the land cover types defined as natural assets relevant to that NCP (for example, only forests for timber but excluding forests for grazing; Supplementary Table 3).

Optimization of NCP

Using integer linear programming (prioritizr, Supplementary Note 1 Section 22), we identify minimum areas required (1) within each country’s land borders and marine EEZs for the local NCP and (2) within all global land area (excluding Antarctica) and all countries’ combined EEZ area for the global NCP, to reach target levels (ranging from 5% to 100%) of every NCP. Before optimization the data were re-sampled to a 2 km equal-area projection for consistency (Supplementary Note 1 Section 24). The optimization selects for the highest values across all NCP, providing the most benefit and/or to the greatest number of people, but not accounting for complementarity or redundancies of adjacent pixels (that is, not dynamically optimizing after each pixel’s selection). Land and marine borders are based on Flanders Marine Institute (2020; Supplementary Table 2), and overlapping claims are excluded from the national analyses. The 12 local NCP are optimized for each country, then aggregated globally, while the two global NCP are optimized globally. In addition to these two main optimizations, we assess the sensitivity to scale by optimizing the 12 local NCP globally (instead of by each country; Extended Data Fig. 5), both with and without the two global NCP, and by substituting different scales of beneficiaries mapping for people downstream and for access to nature (Supplementary Table 6). We also assess the sensitivity of the area and location of critical natural assets (the optimization solution for the 90% target) to different NCP combinations. These variations include optimizing for each NCP individually, and optimizing for all NCP but dropping each local NCP from the set of 12 to evaluate its effect on the overall optimization (Supplementary Table 6). We also examine the correspondence between NCP and the robustness of these different solutions, by calculating the percentage of area shared by different pairs of services (Supplementary Table 4) or the percentage of area shared by all solutions (Supplementary Data 1). We summarize the land and ocean areas required by country in Supplementary Data 1.

Number of people benefitting from critical natural assets

We map the areas benefitting from critical natural assets to calculate the number of direct beneficiaries of these assets, and to compare the number of beneficiaries to the number of people living on the lands comprising these critical natural assets. For this analysis we are only able to include NCP for which the flow of the benefit can be spatially delineated: downstream water quality regulation (sediment retention and nitrogen retention), flood mitigation, nature access, fuelwood provision and coastal risk reduction. The benefitting areas of some of the material NCP that are traded (fish, timber, livestock and crops that are pollinated) or the location of people who buy those traded goods are not easily mapped, so people benefitting from these NCP are not included in this analysis of beneficiaries. However, people within 1 h of critical natural assets may provide a surrogate for many of the material NCP that are locally consumed. For water quality regulation, we take the population within the areas downstream (Supplementary Note 1 Section 1) of critical natural assets. For nature access, we take the population within 1 h travel (by foot, car, boat or rail; Supplementary Note 1 Section 8) of critical natural assets. Likewise, for the other NCP we take the relevant population downstream, within the protective distance, or a gathering distance of critical natural assets. The relevant population for each NCP is considered to be the total global population for nature access and water quality regulation, but is limited to the total population living within 10 km of floodplains for flood mitigation, population along coastlines (in exposed areas: <10 m above mean sea level) for coastal risk reduction, and rural poor populations for fuelwood. Total ‘local’ beneficiaries are calculated through the intersection of areas benefitting from different NCP, to avoid double-counting people in areas of overlap. We calculate the number of people and percent of relevant population benefitting globally for each NCP (Fig. 2b) and the total ‘local’ beneficiaries globally (Fig. 2a) and by country (Supplementary Data 2).

Overlap analysis

We evaluate how well local and global critical natural assets align spatially with each other, and with biodiversity (terrestrial vertebrate species AOH27; Supplementary Note 1 Section 15) and cultural diversity (proxied by the number of Indigenous and non-migrant languages28; Supplementary Note 1 Section 16), to identify synergies between these different potential priorities. To examine the level of overlap between areas identified as critical for the 12 local NCP versus the 2 global NCP, we calculate the area (globally and by country; Supplementary Data 1) where local NCP are selected by the optimization and global NCP are not, where global NCP are selected by the optimization and local NCP are not, and where both are selected by their respective optimizations (the overlap). To calculate the species and languages represented by critical natural assets, we count the number of species whose AOH area targets overlaps these areas, according to both log-linear representation targets (Supplementary Data 3) and targets based on IUCN Red List criteria (Supplementary Data 4), and the number of languages partially intersecting these areas (Supplementary Data 5) globally and within each country (for more detail, see Supplementary Note 1 Section 23).

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.