Where Might We Find Ecologically Intact Communities?

Conservation efforts should target the few remaining areas of the world that represent outstanding examples of ecological integrity and aim to restore ecological integrity to a much broader area of the world with intact habitat and minimal species loss while this is still possible. There have been many assessments of “intactness” in recent years but most of these use measures of anthropogenic impact at a site, rather than faunal intactness or ecological integrity. This paper makes the first assessment of faunal intactness for the global terrestrial land surface and assesses how many ecoregions have sites that could qualify as Key Biodiversity Areas (KBAs – sites contributing significantly to the global persistence of biodiversity) based on their outstanding ecological integrity (under KBA Criterion C). Three datasets are combined on species loss at sites to create a new spatially explicit map of numbers of species extirpated. Based on this map it is estimated that no more than 2.9% of the land surface can be considered to be faunally intact. Additionally, using habitat/density distribution data for 15 large mammals we also make an initial assessment of areas where mammal densities are reduced, showing a further decrease in surface area to 2.8% of the land surface that could be considered functionally intact. Only 11% of the functionally intact areas that were identified are included within existing protected areas, and only 4% within existing KBAs triggered by other criteria. Our findings show that the number of ecoregions that could qualify as Criterion C KBAs could potentially increase land area up to 20% if their faunal composition was restored with the reintroduction of 1–5 species. Hence, if all necessary requirements are met in order to reintroduce species and regain faunal integrity, this will increase ecological integrity across much of the area where human impacts are low (human footprint ≤4). Focusing restoration efforts in these areas could significantly increase the area of the planet with full ecological integrity.

Conservation efforts should target the few remaining areas of the world that represent outstanding examples of ecological integrity and aim to restore ecological integrity to a much broader area of the world with intact habitat and minimal species loss while this is still possible. There have been many assessments of "intactness" in recent years but most of these use measures of anthropogenic impact at a site, rather than faunal intactness or ecological integrity. This paper makes the first assessment of faunal intactness for the global terrestrial land surface and assesses how many ecoregions have sites that could qualify as Key Biodiversity Areas (KBAs -sites contributing significantly to the global persistence of biodiversity) based on their outstanding ecological integrity (under KBA Criterion C). Three datasets are combined on species loss at sites to create a new spatially explicit map of numbers of species extirpated. Based on this map it is estimated that no more than 2.9% of the land surface can be considered to be faunally intact. Additionally, using habitat/density distribution data for 15 large mammals we also make an initial assessment of areas where mammal densities are reduced, showing a further decrease in surface area to 2.8% of the land surface that could be considered functionally intact. Only 11% of the functionally intact areas that were identified are included within existing protected areas, and only 4% within existing KBAs triggered by other criteria. Our findings show that the number of ecoregions that could qualify as Criterion C KBAs could potentially increase land area up to 20% if their faunal composition was restored with the reintroduction of 1-5 species. Hence,

INTRODUCTION
Ecologically intact ecosystems have long been recognized as an important conservation objective for protection. Wilderness areas were identified as important priorities over 30 years ago (McCloskey and Spalding, 1989;Mittermeier et al., 1998), but more recently there have been attempts to be more explicit about what is being measured when referring to wilderness, with a focus on quantifying intact habitat (Potapov et al., 2008(Potapov et al., , 2017Watson et al., 2016;Kennedy et al., 2019;Grantham et al., 2020) and the values of intact habitat . Many assessments have focused on mapping anthropogenic influence (Human Footprint: Sanderson et al., 2002a;Venter et al., 2016b;Anthromes: Goldewijk et al., 2017; Global Human Modification: Kennedy et al., 2018;Low Impact Areas: Jacobson et al., 2019;Intact Forest: Potapov et al., 2008 and have generated several maps of anthropogenic impact that independently estimate that 20-40% of the earth's terrestrial surface is under low human influence (Riggio et al., 2020). Recently, the integrity of natural ecosystems has been recognized as an important goal in the post 2020 global biodiversity framework (CBD, 2020), and the value of intact forest also been recognized both for biodiversity conservation and for ecosystem services . Systematic identification of wilderness areas with outstanding ecological integrity requires a clear definition of ecological integrity and indicators of intactness that can be used to assess areas with potential ecological integrity at a global scale.
There is no clear definition of what is meant by intactness and the term is used loosely in the scientific literature. However, the KBA Standard, published in 2016, following an extensive global consultation process, provides a global standard for identifying "sites contributing significantly to the global persistence of biodiversity" (IUCN, 2016, p. 9). Five over-arching criteria are used to identify KBAs, with Criterion C identifying sites of globally outstanding ecological integrity. The KBA Standard defines Criterion C KBAs in terms of ecological intactness as: A site is one of ≤ 2 per ecoregion characterized by wholly intact ecological communities, comprising the composition and abundance of native species and their interactions. The KBA Standard further defines "intact ecological community" as: "An ecological community having the complete complement of species known or expected to occur in a particular site or ecosystem, relative to a regionally appropriate historical benchmark, which will often correspond to pre-industrial times." The KBA Standard provides some guidance on how ecological integrity should be measured and states that it "should be observed or inferred from both direct measures of species composition and abundance/biomass across taxonomic groups (particularly for species indicative of long-term structural stability and functionality or those known to be highly sensitive to human impact) and absence (or very low levels) of direct industrial human impact (as quantified by appropriate indices at the scale of interest and verified on the ground or in the water)." There is also a recommendation that Criterion C KBAs "should ideally be delineated to be at least 10,000 km 2 in size, within the confines of manageability" and "based on historical bounds of variation using a regionally appropriate benchmark (e.g., the past 500 years)." Ecological integrity needs to be assessed against some regionally appropriate historical benchmark. Species distributions have been changing for millennia because of human activity [e.g., the loss of large mammals across continents (Sandom et al., 2014)] and we chose to focus on changes that have occurred since the year 1500 AD, because this is the baseline date for assessing species extinctions within the IUCN Red List of Threatened species (IUCN, 2020).
Building on these definitions, we distinguish three types of intactness which together make up ecological integrity. Habitat intactness occurs where there is no anthropogenic habitat conversion/transformation. Habitat intactness can be estimated and monitored using remote sensing techniques and measures of the impact of the human footprint on ecosystem structure, although it is easier to estimate for forests than some other habitat types. However, habitat intactness does not necessarily equate to faunal intactness (Benítez-López et al., 2019;Plumptre et al., 2019). Several studies have shown that habitats that appear intact are often missing large, wide-ranging species that play important ecological roles in the ecosystem (Morrison et al., 2007;Plumptre et al., 2019). Medium-size and even small species also play key ecological and functional roles in ecosystems (e.g., small carnivores, bats), significantly contributing to integrity and intactness (Camargo-Sanabria et al., 2015;Brodie et al., 2018). Faunal intactness is therefore differentiated here and occurs where there is no loss of faunal species. We also include the term functional intactness to refer to where there is no reduction in faunal densities below ecologically functional densities to differentiate this intactness at sites from those where species are just present but at diminished density. Identifying sites with outstanding ecological integrity requires an assessment of habitat, faunal and functional intactness.
Efforts to measure intactness to date have mostly mapped human influence that is visible from remotely sensed data (e.g., roads, human settlement, night lights, forest canopy intactness, etc.) or where global-scale data are available from other sources of human activity (e.g., human population density; livestock density, etc.), and thus provide measures of the human footprint on habitat/land cover change. Habitat loss and resulting ecosystem degradation is one of the key drivers of species loss (Chase et al., 2020) but it is not the only driver. While measurements of the Human Footprint (HF) on habitat are useful, they do not fully capture the effects of human activities on biodiversity. Overexploitation, effects of invasive species, and disease, for example, are more difficult to assess using HF-type datasets. In particular, hunting is one threat that is much more difficult to assess because hunters can penetrate large areas of intact habitat, spending days by moving camps, and this activity is impossible to detect remotely (Redford, 1992;Peres et al., 2006). Benítez-López et al. (2019) provided the first attempt to map the impact of hunting across the tropical forests of the world, developing a model of hunting impacts using data from many local hunting studies from South America, Africa, and Asia, and predicting mammal defaunation (both extirpation and reduction in density) across the tropical forested region. This analysis showed that at least 50% of the tropical forests likely have partial defaunation of their mammal biota, with 52% of Intact Forests and 62% of Wilderness Areas (where the Human Footprint ≤ 2) being partially depleted of large mammals. Clearly habitat loss and hunting have compounded impacts on mammal populations (Romero-Muñoz et al., 2020); reducing the area of occupied habitat by more than half, with 29% of the total area attributed to hunting activities (Gallego-Zamorano et al., 2020).
Where species have fallen below a "functional density" at a site, it would not be considered functionally intact because important biotic interactions are likely no longer fulfilling their ecological role (Soulé et al., 2003). Overexploitation, invasive species and disease may reduce species below functional densities leading to loss of functional intactness even if species are not completely extirpated. Functional densities are particularly important for species that play key ecological roles in an ecosystem, such as seed dispersal, nutrient cycling or top-down regulation (Estes et al., 2011;Camargo-Sanabria et al., 2013;Peres et al., 2016). Where mesopredators or competitive species have been introduced by humans, these can have an impact on both species composition and on ecosystem function (Gordon and Letnic, 2016). Assessment of functional intactness should therefore assess both loss of species at a site (faunal intactness) as well as decline in the functional role of individual species as a result of human influence.
To support identification of sites with outstanding ecological integrity, we provide an initial scoping analysis and explore how measuring habitat, faunal and functional intactness narrow down the areas where sites might be identified as having outstanding ecological integrity under KBA Criterion C. More detailed analyses specific to single ecoregions would be needed for a complete assessment of Criterion C sites. We assess how the percentage of land that would qualify for Criterion C KBA would change if we consider just habitat intactness or if we include faunal and functional intactness, which is arguably more difficult to measure. We compile data across multiple taxonomic groups, but the focus is primarily on mammals, particularly large-and medium-sized mammals, because these are often the first species to be affected by hunting (Redford, 1992;Atwood et al., 2020;Gallego-Zamorano et al., 2020). We then assess what proportion of the terrestrial surface of the world could be considered faunally and functionally intact at a minimum. We compiled three sets of data layers, representing areas that have low human footprint, are faunally intact, and functionally intact. We applied these layers cumulatively to identify areas that are ecologically intact (i.e., low human footprint and faunally intact and functionally intact). We then calculated the percentage of the earth's land surface and the number/percent of ecoregion with at least one polygon remaining in the different grades of intactness, and identified areas with potential for restoration of ecological integrity through reintroduction of extirpated species.

MATERIALS AND METHODS
To undertake a global scoping of ecological integrity we used three main data layers and an analytical process to: (1) Identify areas where human impact is low; (2) Assess where species have been historically extirpated; and (3) Assess where species densities are low where this is likely to have been caused by human influence.

Areas of Low Human Impact
We used the Human Footprint (HF) map (Venter et al., 2016a) as our index of human impact. The HF is a composite index of eight global pressures and ranges between 0 and 50. HF scores of four or more are considered modified habitat; scores below this are considered as low impact (Venter et al., 2016b;Mokany et al., 2020;Williams et al., 2020), although, in some recent papers, a value of below three is considered low impact (Di Marco et al., 2018;O'Bryan et al., 2020). The HF mainly maps human infrastructure, population or agriculture, but does attempt to predict impacts up to 15 km on either side of linear features (rivers and roads). For our scoping we selected a threshold value for the HF of four or lower, to ensure that sites with scores between three and four that might still meet faunal intactness were considered. Following the methodology in Plumptre et al. (2019), we smoothed the global HF scores before identifying where the HF score is equal to or lower than four. This was done to eliminate isolated micro-patches of high or low HF scores that would add to analysis time while providing no substantial benefit to global findings. We then converted the result to a vector file. This process developed our baseline map of estimated habitat intactness.
We calculated the percentage loss of area for each ecoregion (Dinerstein et al., 2017) and the average percentage loss across all ecoregions where there was habitat with a HF score ≤ 4. Criterion C sites are suggested to be large (10,000 km 2 or larger). However, our initial results showed that few ecoregions can meet this size with complete functional intactness. We therefore also assessed how the number of ecoregions and the area of land surface would change when selecting sites with minimum areas of 1,000, 5,000, and 10,000 km 2 . We excluded Antarctica as it is not designated as an ecoregion but included within a category of rock and ice (Dinerstein et al., 2017).

Areas of Where Fauna Have Not Been Extirpated
We used three data layers that provided measures of species extirpation at a site: two global, and one focused on tropical forested areas. The first layer was an estimate of the number of historical mammal extinctions across the world (Mammal extirpations). The data for this were taken from Faurby et al. (2018), with the original data mainly coming from Faurby and Svenning (2015) and IUCN 2016-3. We note that these data are not explicitly linked to any specific time point but rather compare current ranges with estimated non-human modified ones. However, we omitted any global or continental extinctions happening before 1500 AD. Whilst a few species have had prolonged range contractions within the continents they still occupy, the data are generally an adequate approximation of the mammalian community if 1500 AD is used as a baseline and have previously been used as such (Belote et al., 2020). The resolution of this map was at ∼ 100 × 100 km and we downscaled it to 10 × 10 km for this study using a bilinear interpolation in ArcGIS 10.5.
The second layer was a map of extinct and possibly extinct range areas compiled from the IUCN Red List of Threatened Species for all taxonomic groups assessed by June 2020 (IUCN 2020 -henceforth IUCN Red List). The IUCN Red List assesses where species have been extirpated (extinct) since the year 1500 AD and this is recorded in the "Presence" field of the range polygons (Plumptre et al., 2019). We updated the map of all extinct ranges provided by Plumptre et al. (2019) for the IUCN Red List 2020 v2. Range mapping on the IUCN Red List is not consistent and has changed over time, so that not all species have their extinct range mapped, but it provides a measure across many taxonomic groups. For this analysis we excluded the extinct ranges of mammal species as these were included already in the first layer (Mammal extirpations). We refer to this layer as non-mammal extirpations (Non-mammal extirpations). Because IUCN Red List Assessors do not consistently map extirpated range for species when making assessments, we preferred to use the first layer for mammals as the mammal layer uses the same approach for nearly all mammal species.
To account for the localized effect of hunting, we used a third layer, a binary raster map of defaunation (i.e., defaunated or not) within tropical forest areas which represented where tropical mammals are predicted to be locally extirpated due to hunting pressure (Gallego-Zamorano et al., 2020). This dataset contains tropical mammal species with at least 95% of their range within the tropical biome. In its preparation, each species was modeled separately at 1 km resolution using presence information in hunted and unhunted areas. A binomial mixed-effects model that quantified the species-specific probability of persistence under hunting pressure was used to identify local extirpations per species (Gallego-Zamorano et al., 2020). We calculated this layer by aggregating all species to create a map of predicted numbers of large and medium mammals extirpated within the areas identified as "defaunated" in the binary map (Tropical forest hunting extirpations).
From these three extirpation maps we created two outputs: (1) A binary map where at least one species has been extirpated by combining all three maps and selecting all cells where extirpated species exceed zero; and (2) A raster of numbers of species extirpated at a scale of 10 × 10 km which was derived by summing the numbers of mammal extirpations (Belote et al., 2020) and the non-mammal extirpations generated for this analysis, scaled to the same 10 × 10 km resolution.

Assess Areas Where Species Are Reduced Below Functional Density
We compiled maps for widespread species for which predicted density or probability of suitable habitat has been modeled across their global range. We targeted wide-ranging large mammals (including large ungulates, great apes, and carnivores), and landscape species that are often the first species to be lost from an ecological community because of anthropogenic impacts (Sanderson et al., 2002b). While this is not a comprehensive mapping of potential reduction in densities of species across the world, these species may be considered indicators of ecological integrity because they are sensitive to human impacts. We assembled density data for the following species or subspecies: African forest elephant (Loxodonta africana cyclotis) (Maisels et al., 2013), western lowland gorilla (Gorilla gorilla gorilla) and central chimpanzee (Pan troglodytes troglodytes) (Strindberg et al., 2018), western chimpanzee (P.t.verus) (Heinicke et al., 2019), Grauer's gorilla (Gorilla beringei graueri) (Plumptre et al., 2016), Bornean Orangutan (Pongo pygmaeus) (Voigt et al., 2018), Sumatran Orangutan (P. abelii) (Wich et al., 2016), Tapanuli Orangutan (P. tapanuliensis) (Nater et al., 2017), jaguar (Panthera onca) (Jȩdrzejewski et al., 2018), lion (Panthera leo) (Riggio et al., 2013), and African buffalo (Syncerus caffer) (Wint and Kiara, 2017). In addition, species distribution models showing probability of suitable habitat were compiled for eastern chimpanzee (P. t. schweinfurthii) (Plumptre et al., 2010), and four bear species: polar bear (Ursus maritimus), brown bear (Ursus arctos), American black bear (Ursus americanus) and Asian black bear (Ursus thibetanus) (Luna-Aranguré et al., 2020). Threshold densities or probabilities were identified for each species and for this study were assumed to predict where a species is at low abundance and will not be playing its functional role in the ecosystem (Supplementary Table 1). We compiled these into an overall map showing where densities are likely to be low (Selected large mammal reduced density).
As a second input layer to measure reduced density, Benítez-López et al. (2019) mapped areas of predicted defaunation (defaunation map) resulting from hunting across tropical forests. Their maps show the areas where mammal species are predicted to have declined in abundance, but do not identify areas of local extinction. However, the same data were used to predict where species have been extirpated by hunting (Gallego-Zamorano et al., 2020) -our Tropical Forest Hunting Extirpation map above -and we used this as a mask to remove areas from the defaunation map which have lost species, thereby leaving sites where species have been reduced in absolute abundance, but not lost. We used a threshold of <0.9 (i.e., more than 10% average decline in abundance across mammal species) as recommended by Benítez-López et al. (2019) on their defaunation score to represent areas where species were below functional density (Tropical forest hunting reduced density).

Analyses
We assessed the cumulative effects of these layers for potential Criterion C sites by sequentially applying the following filters to Frontiers in Forests and Global Change | www.frontiersin.org vector files that we had derived from the various maps described above (Figure 1):

Habitat intactness: Select all areas where Human
Footprint ≤ 4; 2. Faunally intact: Subtract areas from 1 (Habitat intactness), where one or more species have been extirpated; 3. Functionally intact: Subtract areas from 2 (Faunally intact) where selected large mammal density is low.
We assessed the effect of limiting sites to a minimum polygon area of 1,000, 5,000, and 10,000 km 2 for each ecoregion. We calculated changes in numbers and percentage area of ecoregions and calculated median polygon area for each ecoregion that had at least one polygon for each minimum area class. We then calculated the total percentage of land that would meet a definition of faunal and functional intactness for each minimum area class.
We assessed which biomes of the world would contain areas of faunal intactness that are at least 10,000 km 2 by intersecting the polygons generated above with the map of terrestrial biomes generated by Dinerstein et al. (2017) as attributes in their ecoregions shapefile. For the ecoregions that were identified as having functionally intact areas, we also assessed what percentage of these areas was covered by protected areas using the World Database of Protected Areas (UNEP-WCMC and IUCN, 2020) and percentage covered by existing KBAs using the World Database of KBAs (BirdLife International, 2020). We measured this coverage using the same threshold areas above to evaluate changes with minimum polygon area.
Finally, having mapped the number of species extirpated in each 10 km × 10 km cell, we also mapped the extent of land where Human Footprint ≤ 4 and where it might be possible to restore ecological integrity by assessing where few species (1-5) have been extirpated and might be reintroduced if they have not gone extinct globally. The number of species extirpated did not correlate well with the Human Footprint score ( Supplementary  Figure 1) but the requirement of Criterion C sites to be in areas with low human impact meant we restricted our analyses to the area of low human impact (Human footprint ≤ 4). We first removed the number that were globally extinct on the IUCN Red List from the map of number of species extirpated, so that we could assess areas where restoration might be possible in land where HF ≤ 4. Applying the three minimum area classes, the total area of a site was computed for sites incorporating 0-5 species extirpated (for example a 10,000 km 2 site assessed for 2 species extirpated could contain areas within the 10,000 km 2 where 0, 1 or 2 species were extirpated).

Areas Under the Types of Intactness
Habitat intactness (as measured by HF ≤ 4), covered 43.4% of the terrestrial surface of the planet. When selecting areas of land exclusively based on the HF ≤ 4 score, only 214 of the 847 ecoregions (25.3%) qualify as having polygons of 10,000 km 2 or larger. This rises to 464 ecoregions (54.8%) if areas as small as 1,000 km 2 are considered. Removing each of the individual faunal layers from this baseline habitat intactness map for a minimum area of 10,000 km 2 showed that most of the reduction in percentage of ecoregions or terrestrial land surface occurs as a result of mammal and non-mammal extirpations and the reduced densities of selected species (Figure 2).
Combining the three layers measuring extirpation (Tropical forest hunting extirpation, Non-mammal extirpation, and Mammal extirpation) into a single mask layer enabled us to assess faunal intactness. Removing regions where species have been lost from the HF ≤ 4 map has a huge impact, reducing the number of ecoregions to 42 at 10,000 km 2 or 124 at 1,000 km 2 (Figure 3A; Supplementary Figure 2). Only 2.9% of land surface ( Figure 3B) could be considered faunally intact for a minimum area of 10,000 km 2 (3.4% at 1,000 km 2 ). Finally, we combined the two layers measuring low density (tropical forest hunting reduced density and selected large mammal density) to assess functional intactness (i.e., faunally intact areas which also have functional densities of species). This reduced the number of ecoregions that would qualify further from 38 at 10,000 km 2 to 119 at 1,000 km 2 and only 2.8% of the land surface can be considered Frontiers in Forests and Global Change | www.frontiersin.org functionally intact for a minimum area of 10,000 km 2 (3.3% at 1,000 km 2 ).
Approximately 25% of the world's terrestrial ecoregions have at least one habitat intact polygon >10,000 km 2 remaining; whereas only 5.0% have at least one faunally intact polygon remaining; and 4.5% have at least one functionally intact polygon remaining. The median percentage area of ecoregions that remain in the analysis also drops as we sequentially filter for HF ≤ 4, faunal intactness, and functional intactness (Figure 3C). This shows that not only do the numbers of ecoregions decrease, but also the area remaining within the ecoregions that can be considered intact decreases greatly.

Biome Representation
A few biomes -tundra, deserts, tropical and subtropical grasslands -have an average of more than 40% of each ecoregion in faunally and ecologically intact areas of >10,000 km 2 . Other biomes have a much lower percentage (Figure 4 and Supplementary Figure 2).

Proportion of Intact Land Covered by Protected Areas or KBAs
Only about 11% of the area identified as functionally intact is included in existing protected areas (categories I-IV). The median percentage of functionally intact land covered by protected areas of all types is larger but still only 34% (Figure 5). The percentage of functionally intact land area covered by KBAs identified with other KBA criteria is only 4% at 10,000 km 2 but rises to 7% at 1,000 km 2 .

Potential for Recovery of Functional Intactness
The very low percentage of land that might qualify as functionally intact led us to assess what remaining land exists where only five or fewer species have been extirpated (but are not globally extinct) within the area of habitat intactness (Figure 6). We calculated changes in land surface and number of ecoregions that could in future become faunally intact (and eventually functionally intact) if restoration occurred by sequentially adding 1, 2, up to 5 extant species where these have been extirpated (Table 1). Restoring 1-5 species would increase the percentage of the earth's land surface that is faunally intact from approximately 3% at present up to 19% (10,000 km 2 minimum area) to 24% (1,000 km 2 ).

DISCUSSION
While wilderness areas are increasingly recognized as important for biodiversity conservation, few areas of the world remain that can be characterized as having outstanding ecological integrity, (i.e., retaining intact species assemblages at ecologically functional densities). We found that only 2.8% of the terrestrial surface of the planet is represented in areas of 10,000 km 2 or larger with low human footprint, no known species loss and no species known to be reduced below functional densities. This compares with estimates of 20-40% from mapped habitat intactness in the literature (Riggio et al., 2020).
The datasets we have applied in our analysis are the best we could compile at a global scale; we recognize that at an ecoregional scale it will be possible to make more detailed analyses of human footprint, identification of areas which have experienced species loss, and mapping of population densities of species likely impacted by humans to assess whether they are likely functional. In particular, we selected a subset of large mammals to assess functional density, choosing species that were relatively widespread and likely hunted, as surrogates for other species. Few species have their density mapped across their global range and we compiled what we could find that was relevant to the analysis. We strongly encourage more detailed analyses to be made to start scoping where ecologically functional KBAs might be identified and that sensitivity analyses are made using variation in the datasets. Detailed assessments will likely reduce the area that can be described as ecologically intact still further. However, detailed assessments may also reveal that smaller sites are sufficient to maintain populations and ecological processes in a few ecoregions, especially on islands or island-like ecosystems with few area-demanding species. For example, the East African montane moorlands ecoregion is a high altitude ecoregion where few large-bodied fauna occur naturally. The total area of this ecoregion is less than 5,000 km 2 but it could potentially have areas that are considered functionally intact. We therefore explored the effect of relaxing the area requirement, but emphasize that sites that are smaller than 10,000 km 2 would still need to maintain functionally intact populations of all species for that ecoregion to qualify as KBAs under Criterion C.
A surprisingly low percentage of the area identified as ecologically intact in our analysis is covered by protected areas (11% covered by PAs in categories I-IV), indicating that Other Effective Conservation Measures (OECMs) have an important role to play in maintaining ecological integrity. Many of the areas identified as ecologically intact coincide with territories managed by indigenous communities, who have played a vital role in maintaining the ecological integrity of these areas (Figure 6). We show that there is potential to restore ecological integrity over a much broader area. We found that up to 20% of the earth's land surface could be restored to faunal intactness through the reintroduction of up to five species. Clearly restoring faunal intactness will involve meeting the requirements for successful species' reintroduction and the ability to build up population numbers so that they can reach functional densities. Applying the IUCN/SSC reintroduction and translocation guidelines (IUCN/SSC, 2013) will be necessary to minimize disease risks and possible conflicts with people, but our findings show there is potential to increase faunal intactness, and ultimately functional intactness over time, if species can be reintroduced to sites with low human footprint. In this "decade of restoration" (Aronson et al., 2020), we believe that, as well as focusing on restoring habitat, a proactive approach should assess which areas would be amenable to restore faunal intactness with the aim of regaining full functional intactness, by tackling the threats that led to the loss in the first place. These areas will likely be located in eastern Russia, northern Canada and Alaska, the Amazon Basin, parts of the Sahara, and the Congo Basin. Planning for ecological restoration to  include the reintroduction of extirpated but globally extant species, as has occurred in several conservation projects around the world (Hale and Koprowski, 2018), could significantly increase the number of ecoregions that could achieve sites with outstanding ecological integrity. Examples would include reintroducing forest elephants in areas of the Congo Basin where they have been extirpated, or reintroducing some of the large ungulates that have been lost from much of Africa's woodlands and savannas because of overhunting (e.g., buffalo, giraffe, zebras etc.), as long as overhunting has ceased. We acknowledge that other key factors need be taken into account for successful species reintroduction projects, including but not limited to adequate source populations, pedigree and genetics, and disease assessments (Dinca et al., 2018;Hunter-Ayad et al., 2020). Restoration projects should also consider assessing which ecoregions currently cannot meet the requirement of faunal intactness and determine where they could best support restoration in these ecoregions to rebuild faunal intactness and, over time, functional intactness as densities rebuild. The method identified sites with the minimum area class for polygons with zero up to five species extirpated somewhere within the polygon.
Surprisingly, human footprint indices did not correlate well with faunal intactness (Supplementary Figure 1), because while threats such as habitat loss may correlate well, several major threats to faunal intactness (e.g., hunting, invasive species) are not properly accounted for by the human footprint maps (Benítez-López et al., 2019;Belote et al., 2020). There is considerable geographic variation in patterns of hunting pressure, but many tropical forest regions characterized by relatively low human densities and land cover change have also suffered from extensive and unsustainable hunting pressure (Redford, 1992;Tilker et al., 2019); invasive predators have encroached into wilderness areas far from infrastructure networks (Doherty et al., 2016); and diseases such as ebola can spread deep into remote forests (Rizkalla et al., 2007). None of these threats are picked up by remote sensing tools. Threats such as hunting primarily affect large-and medium-sized mammals (Benítez-López et al., 2019; Gallego-Zamorano et al., 2020), although regional variations exist because of hunters' preferences, cultural beliefs, and taboos regarding the consumption of some species (Ngoufo et al., 2014;Bobo et al., 2015). For instance, in Uganda ungulates but not primates are targeted by hunters, unlike just over the border in the Democratic Republic of Congo where primates are widely hunted as well (A. Plumptre pers. obs.). Along the Gabon and Congo coastline, chimpanzees but not gorillas are considered taboo, and this is reversed in some areas of Cameroon; a factor "eat gorillas" or "eat chimpanzees" was included in the spatial model that predicted the density of both taxa in western Equatorial Africa (Strindberg et al., 2018). Predicting the impacts of hunting needs much more detailed analysis (Benítez-López et al., 2019;Bogoni et al., 2020;Romero-Muñoz et al., 2020) than making assumptions about impacts from access distance to road and river networks. Similarly, the impacts of introduced species cannot be assessed from maps such as the Human Footprint but have had a major impact on extirpated species numbers in places such as Australia (Woinarski et al., 2015).
Nevertheless, human footprint indices can play an essential role in identifying areas with potential for ecological restoration through reintroduction of extirpated species, as there will typically be fewer barriers to successful restoration in areas with low human footprint. As the world develops the next goals and targets for the Convention on Biological Diversity under the post 2020 global biodiversity framework, "intact habitat" has been recognized as an important target. We would argue here that "intact habitat" is clearly an important conservation target, and that we must also consider faunal and functional intactness, and direct resources to areas that still retain a fully functioning ecological community as this is clearly very rare. Retaining large areas will give these sites a better chance of retaining their functional intactness. Working to restore ecological integrity to those areas of the world that are close to ecological integrity should be a focus of the current decade of restoration also. Moreover, identifying Criterion C KBAs could help focus attention on these sites for conservation action and investment.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
AP and DB conceived the original idea, compiled layers, made analyses, and drafted manuscript. RB, EV-D, SF, WJ, HKi, HKü, AB-L, CL-A, MV, SW, WW, and JG-Z provided data layers and contributed to the manuscript. CB contributed ideas on manuscript structure and analyses. All authors contributed to the article and approved the submitted version. Supplementary Table 1 | The species name, threshold density and source of data for rangewide species distribution maps for the species that were used to assess impacts on functional intactness.