Abstract
Birds are among the best-studied animal groups, but their prehistoric diversity is poorly known due to low fossilization potential. Hence, while many human-driven bird extinctions (i.e., extinctions caused directly by human activities such as hunting, as well as indirectly through human-associated impacts such as land use change, fire, and the introduction of invasive species) have been recorded, the true number is likely much larger. Here, by combining recorded extinctions with model estimates based on the completeness of the fossil record, we suggest that at least ~1300–1500 bird species (~12% of the total) have gone extinct since the Late Pleistocene, with 55% of these extinctions undiscovered (not yet discovered or left no trace). We estimate that the Pacific accounts for 61% of total bird extinctions. Bird extinction rate varied through time with an intense episode ~1300 CE, which likely represents the largest human-driven vertebrate extinction wave ever, and a rate 80 (60–95) times the background extinction rate. Thus, humans have already driven more than one in nine bird species to extinction, with likely severe, and potentially irreversible, ecological and evolutionary consequences.
Similar content being viewed by others
Introduction
Humanity’s dispersal out of Africa and subsequent peopling of essentially all ice-free land across the globe has triggered waves of extinctions, which for many groups, including birds, have been particularly large across isolated archipelagos1,2,3 (Fig. 1). Even small human populations rapidly devastated island avifaunas as they introduced new threats outside the evolutionary experience of native species2. Drivers of human-driven bird extinctions include habitat loss associated with land clearance (cutting, burning) and the introduction of non-native plants and crops4, the introduction of alien species (domestic animals and/or human commensals)2,5,6 and the overexploitation of bird species via hunting and trapping (birds were hunted for their fat, protein, bones and colourful feathers)1,5. Yet, due to the incomplete avian fossil record3,5,7,8, many bird extinctions are likely to have gone undiscovered—i.e., not yet found due to a lack of research effort or because they left no trace in the discoverable fossil record4,7,9,10,11,12. Hence, the full extent of bird extinctions since the Late Pleistocene remains unknown (although see previous studies of the Pacific4,7,10,12). Previous analyses of bird extinction rate have therefore focused on well-recorded observed extinctions2,13,14, i.e., since 1500 CE (Common Era). However, ignoring fossil and undiscovered extinctions limits our understanding of avian extinction rates15 and can substantially underestimate the magnitude of human-driven biodiversity loss16, with implications for global biodiversity, evolutionary history, and conservation17,18,19,20,21.
Here, we go beyond previous studies of the Pacific4,7,10,12 to cover the globe, evaluating the spatial distribution of extinctions, and incorporating undiscovered extinctions, and their uncertainty, into temporal analyses of extinction rate. Hence, we quantify the total magnitude, distribution and rate of bird extinctions worldwide since the Late Pleistocene, including undiscovered extinctions. To estimate undiscovered extinctions, we first modelled fossil bird extinctions across 69 archipelagos (1488 islands; Supplementary Fig. 3) as explained by multiple environmental predictors and the completeness of the fossil record, as indicated by research effort (see Methods). We focus on archipelagos for the model of undiscovered extinctions because ~90% of recorded bird extinctions have occurred on islands2,22. We combined these estimates of undiscovered extinctions (archipelagos only) with estimates of fossil and observed extinctions across the globe. Subsequently, we estimated the extinction date for all bird species lost since the Late Pleistocene and inferred extinction rates of birds through time (see Methods). Our findings suggest that ~1300–1500 bird species (~12% of the total) have gone extinct since the Late Pleistocene, more than double the current estimate from the observed and fossil records alone. Moreover, we identify an intense human-driven extinction wave for birds (i.e., caused directly by human activities such as hunting, as well as indirectly through human-associated impacts such as deforestation, fire, and the introduction of invasive species), peaking ~1300 CE with a rate 80 (60–95) times the background extinction rate.
Results
Bird extinctions since the Late Pleistocene
Globally, we estimate 1430 (95% credible interval: 1327–1544) bird extinctions since the Late Pleistocene (Fig. 2a). Given current estimates of 10,865 extant bird species, this suggests that 11.6% (10.9–12.4) of all bird species (approximately one in nine) have gone extinct over the last ~126,000 years, with human activities likely contributing to almost all these extinctions. We estimate that 55% of these extinctions are currently undiscovered (788 species; 685–900). Thus, the total number of human-driven bird species extinctions is more than double the current estimate from the observed and fossil records alone (Fig. 2a). For the Pacific we estimate 875 (773–973) total bird extinctions with 554 (450–652) of these undiscovered, compared to 557 (508–605) total extinctions outside the Pacific with 235 (185–281) undiscovered. Hence, we estimate that the Pacific accounts for 61% of total bird extinctions and 70% of undiscovered bird extinctions.
Bird extinction rate since the Late Pleistocene
Bird extinction rate varies through time, with three major extinction waves since the Late Pleistocene (Fig. 3a and Supplementary Fig. 4). The first, peaking ~840 BCE (Before Common Era; Fig. 3a), was primarily driven by the arrival of people to islands across the Western Pacific (e.g., the Mariana, Tongan and Fijian Islands), as well as the Canary Islands (Fig. 3b). Known extinctions potentially associated with this wave include the Fiji Teal (Anas sp.) from Fiji, the New Caledonian Cockatoo (Cacatua sp.) from New Caledonia, and the Long-legged Bunting (Emberiza alcoveri) from the Canary Islands (Fig. 3b). The second, and most intense extinction wave, peaked ~1300 CE (Fig. 3a) with a maximum extinction rate of 1.6 (1.2–1.9) × 10−4/species/year, roughly equivalent to ~160 (120–190) extinctions per million species years (E/MSY). This is approximately 80 (60–95) times the background extinction rate, here set to 2 E/MSY following previous assessments18 (a relatively high estimate of background extinction rate). This wave was driven by the first arrival of people to islands across the Eastern Pacific, especially to Hawaii and the Marquesas Islands, as well as Aotearoa New Zealand (Fig. 3c). The arrival of humans, along with the domestic animals (pigs, dogs, chickens) and Polynesian Rats (Rattus exulans) that were translocated with them, brought habitat transformation and novel threats to the native avifauna, leading to rapid declines of naïve bird species2,5. Known extinctions potentially associated with this wave include the High-billed Crow (Corvus impluviatus) from Hawaii, Sinoto’s Lorikeet (Vini sinotoi) from the Marquesas Islands, and nine Moa (Dinornithiformes) species from Aotearoa New Zealand (Fig. 3c). By contrast, the ongoing extinction wave (Fig. 3a), includes extinctions across multiple disparate regions (Fig. 3d), and is driven by the global intensification of human threats, including habitat destruction, direct exploitation, pollution, and invasive non-native species2,5. Known extinctions associated with this ongoing wave include ʻĀmaui (Myadestes woahensis) from Hawaii, Lyall’s Wren (Traversia lyalli) from Aotearoa New Zealand, and the Colombian Grebe (Podiceps andinus) from continental South America (Fig. 3d).
Discussion
Previous assessments of bird extinction rates that only include observed extinctions (e.g., post-1500 CE2,13,14,18) clearly miss major extinction waves (Fig. 3a). Observed high extinction rates have therefore been considered unprecedented18. However, at least for birds, earlier extinction waves were greater in magnitude (Fig. 3), as has been hypothesized10 but not previously demonstrated. In addition, focusing on observed extinctions disregards the fluctuations in extinction rate between waves. Extinctions have typically occurred rapidly after human arrival2,5,23, but extinction rates have then slowed as vulnerable species were quickly lost and species more resilient to human impacts remained2 (reflected in Fig. 3a). Thus, the apparent steady increase in extinction rate observed since 1500 CE13,14,18 (Fig. 3a) may not be appropriate for estimating the true dynamics of the current biodiversity crisis. Instead, here we demonstrate the high spatial and temporal heterogeneity of human-driven extinctions (Fig. 3).
Nevertheless, extinction debts—the projected extinctions of numerous extant species that are yet to occur24,25—suggest that we may be moving towards a present-day bird extinction rate even greater than any since the Late Pleistocene13,26. Extinction debts can take centuries to be realized27 and thus today’s accumulating extinction debt is still being paid off19,25. A recent study estimated an effective (i.e., accounting for the accumulating extinction debt) extinction rate of 2.17 × 10−4/species/year for the present day13; compared to our maximum rate of 1.6 × 10−4/ species/year for ~1300 CE. Moreover, bird extinction rate over the next 100 years has been predicted26 to reach 6.98 × 10−4 /species/year, a rate that is unprecedented for birds (Fig. 3a). Together these studies predict an additional 226—738 bird extinctions over the next few hundred years13,26. If/when these extinctions are realized it would mean that, combined with our results, up to one in six bird species had been driven extinct. The extinction debt, as well as conservation successes2,13 and the delay in officially listing species as extinct28, could therefore explain the minor decline in extinction rate we detect since 1939 (Fig. 3a). Hence, in the future, extinctions may be identified and backdated, revising this small decline upwards to a likely increase.
Incomplete knowledge of bird extinctions necessarily underestimates the magnitude of species loss15,18. Still, any quantification of undiscovered extinct birds is unavoidably an estimate. Here, we used statistical models and a set of simple, conservative assumptions (in this case, assumptions which should lead to underestimates of the true number). Specifically, we assumed zero undiscovered extinctions since written observational records began and zero undiscovered continental extinctions. We also assumed complete knowledge of the well-studied Aotearoa New Zealand avifauna7,29, as a way of correcting for research effort (see Methods). Violations of these assumptions would result in higher overall extinction estimates. Overall, we provide a broad yet precise estimate of undiscovered bird extinctions.
Previous research has estimated the number of bird extinctions across the Pacific4,7,10,12, with the expectation that most extinctions would be located there (confirmed here as 61% of total bird extinctions). Specifically, these studies use faunal reconstruction4,10 and variants of mark-recapture7,12 to estimate bird extinctions. These estimates range from ~80012, to ~13007, to over 20004 total Pacific bird extinctions, although definitions of the Pacific vary (i.e., which islands are included/excluded), and which bird groups were studied also differ, making direct comparisons difficult. Still, for the Pacific we estimate 875 (95% credible interval: 773–973) total bird extinctions, which falls within the lower to mid-range of previous estimates4,7,10,12. Our Pacific estimate is similar to a previous study that suggested there were more than ~600 but fewer than ~1300 undiscovered extinctions10. By contrast, our estimate is lower than a recent estimate of ~1300 for all birds7, although this study also presents a median estimate for only non-passerine landbirds of 983 (95% credible interval: 731–1332). When considering the uncertainty of both estimates there is strong overlap in the credible intervals, highlighting their general agreement. Our estimate is similar to all previous estimates except for the estimate of over 2000 undiscovered extinctions4. This number has been previously described as a likely overestimate10,30, due to the occurrence of natural events (e.g., volcanism, tsunamis) across some islands assumed to host undiscovered endemic birds that might have prevented colonization or extirpated populations before speciation occurred10,30,31. In addition, our estimate of the ratio of undiscovered to discovered extinctions for the Pacific of 0.63 (95% credible interval: 0.59–0.67) is similar to previous estimates of 0.512 and 0.67 (95% credible interval: 0.46–0.84)7. Thus, our estimate is supported by a range of datasets and analytical techniques, both within our study (see Methods) and across previous studies4,7,10,12.
Overall, we estimate that at least ~1300–1500 bird species have gone extinct since the Late Pleistocene (~12% of the total). Moreover, we identify an intense extinction wave across the Eastern Pacific around the 14th century, which represents the largest human-driven vertebrate extinction wave reported to date. These extinctions have strong implications for our understanding of avian species richness2, ecological diversity19,21 and evolutionary history20. The loss of such a large proportion of bird species, with the identities (and therefore ecological roles) of 55% unknown, suggests underappreciated and potentially irreversible ecological impacts, and thus unknown debts in ecosystem functioning6,21,32,33, particularly across islands, where large native mammals are often absent, and instead birds have adapted to provide important ecosystem functions. Known examples include extinct megaherbviores such as the Elephant Birds (Aepyornithidae) of Madagascar, which influenced plant structure and diversity and ecosystem dynamics34, extinct aerial predators such as Haast’s Eagle (Hieraaetus moorei)35,36, and extinct seed-dispersers such as Seychelles Parakeet (Psittacula wardi)37. We estimate ~800 undiscovered bird extinctions, many of which could also have had key roles in their ecosystems. In addition, bird extinctions can have long evolutionary legacies20. These evolutionary impacts may be even greater for archipelagos with large proportions of undiscovered extinctions (e.g., Hawaii), and would amplify the already disproportionate contribution of islands to bird evolutionary diversity38.
Our study demonstrates the severe and long-term alteration of avifaunas globally1,3,7, representing a far higher human impact on avian diversity than previously recognized. Urgent protection of the remaining native biotas should constitute a high priority for conservation, to avoid a contemporary extinction wave of even greater magnitude than the prehistoric episodes revealed here.
Methods
In brief, we estimated the number and timing of bird extinctions since the Late Pleistocene across the globe. Although our outlook was global, we first modelled fossil bird extinctions across archipelagos, as the majority of known bird extinctions have occurred across islands2,22. We extrapolated our archipelago model to the research effort of Aotearoa New Zealand, constraining the extrapolations based on upper bounds to obtain estimates of undiscovered bird extinctions across archipelagos. Undiscovered extinctions are those that have not been scientifically recorded: some of these may not have left any physical trace, some may be described in the future9. We combined these estimated undiscovered bird extinctions across archipelagos with global fossil (extinct species only known from fossil evidence) and observed (extinct species with written observation records) extinction records to produce a total extinction estimate. We then estimated the extinction rate through time since the Late Pleistocene.
We used R version 4.0.4 (ref. 39), for all our analyses, except for pre-processing of shapefiles of the archipelagos, which we prepared in ArcGIS 10.4 (ref. 40) and rasterization of landmasses, which we ran in Python version 2.7.10 (ref. 41). See https://zenodo.org/records/10014585 for R code summarizing the major analytical steps. We used multiple R packages for data preparation, analysis, and visualization, including arm42 1.13-1, broom43 1.0.4, cowplot44 1.1.1, DHARMa45 0.4.6, dplyr46 1.1.2, gghalves47 0.1.4, ggplot248 3.4.2, HDInterval49 0.2.4, hydroGOF50 0.4-0, jtools51 2.2.1, letsR52 4.0, MASS53 7.3-53, purrr54 1.0.1, raster55 3.6-20, readr56 2.1.4, rsq57 2.5, scales58 1.2.1, scatterpie59 0.1.8, sf60 1.0-13, sp61 1.6-0, tibble62 3.2.1, tidyr63 1.3.0, and zoo64 1.8-12.
Archipelagos and islands
Although most bird species (>80%) live on continents65, the majority of recorded extinctions (~90%) have been on islands22. We therefore acknowledge that the number of continental undiscovered bird extinctions is minor relative to those on islands (i.e., for our analysis we assume zero continental undiscovered extinctions). In addition, post-1500 CE bird extinctions are generally well-documented2,22, and we therefore assume that the majority of post-1500 CE extinctions are recorded (but see ref. 9 for potential undiscovered extinctions post-1500 CE).
To estimate undiscovered bird extinctions across archipelagos we first extracted geographic and environmental data for 17,883 islands (ref. 66), we manually added Norfolk Island as it was missing from the dataset. We then identified 1488 islands from this dataset that were settled by humans before written observational records began (i.e., pre-1500 CE) and had the potential to support endemic bird species—islands with area >5 km2 (refs. 7, 10), isolated from the mainland66,67,68, and not glaciated69 during the Last Glacial Maximum (Supplementary Fig. 3). We compiled the date of first human arrival from the literature (Supplementary Data 2). We removed the islands of Cape Verde (humans first arrived 1462 CE70) and São Tomé and Príncipe (humans first arrived ca. 1470 CE71), as these were first settled close to 1500 CE and have zero fossil extinctions. We therefore conservatively assume that these islands have zero undiscovered extinctions. The area threshold was based on an earlier paper highlighting that only islands >5 km2 are likely to support endemic birds7,10. We only included islands isolated from the mainland, as islands connected to the mainland are likely to be dominated by continental processes and continental fauna. We compiled previous estimates of island area and connection to the mainland66, where connection to the mainland was based on global bathymetry data68 with a sea level of −122 m during the Last Glacial Maximum66,67. To identify islands not covered by glaciers during the Last Glacial Maximum we used the prehistoric distribution of ice sheets69.
We grouped these 1488 islands into 69 archipelagos (Supplementary Fig. 3) according to an existing classification66, with minor modifications (see isl_arch.csv at https://zenodo.org/records/10014585 for details on islands within archipelagos). Specifically, we aggregated archipelagos that were closely related (e.g., Canary Islands with Madeira), and we disaggregated large archipelagos, as informed by biogeographical (e.g., shared mainland) and ecological (e.g., shared species) similarities between islands, see also ref. 72. Archipelagos are considered the most-appropriate spatiotemporal insular unit for large-scale biodiversity analyses73,74,75,76. Here we refer to both true archipelagos (composed of multiple islands) and isolated insular units (single islands; e.g., Cyprus) as archipelagos. All fossil extinct birds were single-archipelago endemics except for nine species (2%): Alopecoenas nui, Caracara creightoni, Gallinago kakuki, Megapodius alimentum, Megapodius molistructor, Tyto noeli, Tyto pollens, Vini sinotoi, and Vini vidivici, which we assigned to their primary archipelago, based on island area. We therefore assume that all undiscovered extinctions were single-archipelago endemic species.
We constructed shapefiles in ArcGIS 10.4 (ref. 40) for each of the archipelagos from the database of global administrative areas (GADM), version 3.6 (gadm.org/data.html).
Data preparation
To estimate the number of undiscovered extinct species (not yet discovered or left no trace) we first modelled the number of fossil extinct bird species across the focal archipelagos. We calculated the number of fossil extinct birds for each of the archipelagos from a recently compiled database72.
We also compiled data for 14 predictors (aggregating island-scale measurements to the archipelago-scale where applicable; Supplementary Data 3), underpinned by biological a priori expectations. Although soil chemistry is known to affect fossil preservation77,78, suitable data at the global scale are too imprecise and the influence of soil properties on preservation is often local and complex78. We manually added data for Norfolk Island, which was missing from66, collating data from WorldClim79, SRTM30 (ref. 80), P. Weigelt (pers. comm., Feb 2020), and adjusting the Lord Howe Island estimate of archipelago plant richness based on the species-area relationship81,82. Here we outline the predictors and our respective a priori expectations.
-
1.
Research effort. We expect greater research effort, as approximated by the number of publications referring to biological fossils, to relate to greater probability of discovering bird fossils.
We estimated research effort as the number of eligible publications referring to biological fossils per archipelago. We followed the Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) guidelines83. To identify publications, we used Scopus (www.scopus.com; accessed 2 March 2020) with the search query string:
(TITLE-ABS-KEY (“Archipelago_names” OR “Island_names”)
AND TITLE-ABS-KEY (fossil)
AND NOT TITLE-ABS-KEY (“fossil fuel” OR “fossil reef” OR “fossil coral” OR “fossil beach*”))
The phrases “fossil fuel”, “fossil reef”, “fossil coral” or “fossil beach” were excluded to ensure the search was focused on terrestrial, biological references. We extracted archipelago names primarily from66, with manual additions based on the literature. We also included up to six alternative (sub-)archipelago names that might be referred to in the literature. We extracted island names from66 and the literature, and included alternative island names where applicable. We included both English and regional archipelago/island names where possible, e.g., Flores Island and Pulau Flores, and Scopus accounts for accented characters, e.g., Yucatan and Yucatán. We ran the search for each of the archipelagos. However, for Aotearoa New Zealand we removed all publications after 2009 (the date of discovery of the last fossil bird species). Thus, our measure of research effort reflects the research effort to date for all archipelagos (potentially incomplete avifaunas), except for Aotearoa New Zealand, where it reflects the research effort required to discover the complete avifauna29. Our search resulted in 4938 total publications (Supplementary Fig. 6). We removed 25 within-archipelago duplicates and screened the remaining publications for eligibility. We screened all publication titles, and when eligibility was uncertain, we read the abstracts and keywords. We excluded publications if the main study location was not the focal archipelago (this occurred where multiple archipelagos had the same island names, e.g. South Island), if the publication focused on marine fossil samples (e.g., benthic samples, fossil coral), if they did not focus on biological fossils (e.g., fossil fuels, volcanic fossils) and/or if they referred to living fossils—resulting in 2139 eligible publications (Supplementary Fig. 6).
-
2.
Isolation distance. Isolation is a major driver of island biodiversity, as isolation facilitates in-situ speciation, which leads to high levels of endemicity in islands (despite islands tending to have lower species richness overall)31,84. Isolation can increase extinction susceptibility, because the avifaunas of isolated archipelagos are likely to have evolved for longer in the absence of mainland predators, and so their species are more likely to react naïvely to both human hunters and exotic mammalian predators when they arrive3,85,86. Archipelagos that are more isolated are likely to have lower rates of colonization, leading to fewer species overall, but higher rates of endemism87, and thus greater probability of global, rather than regional, extinction. A classic example of this phenomenon is Hawaii88. Isolation distance also relates to the overwater distance from the mainland required to reach the archipelago, which acts as a dispersal filter, modifying community composition7.
-
3.
Surrounding landmass. As for isolation distance, surrounding landmass relates to species naïvety to predators. Unlike isolation distance, surrounding landmass accounts for coastline shape of large landmasses and therefore relates more closely to colonization pressure87, although surrounding landmass does not as closely represent the overwater distance required to reach the archipelago (i.e., the dispersal filter).
-
4.
Elevation. Mountainous regions generally have limited human accessibility, providing refuge from human hunters, but also impeding attempts to discover fossils7,89,90,91. Elevation also relates to topographic complexity and environmental heterogeneity, both of which can affect speciation and extinction processes.
-
5.
Temperature. Climatic factors, such as temperature, are key drivers of ecosystem processes and vegetation structure66, therefore influencing speciation, colonization, and extinction. Moreover, climatic factors can mediate human impacts on archipelagos; for example, dry islands are more likely than wet islands to end up deforested92.
-
6.
Precipitation. As for temperature.
-
7.
Temperature variability. As for temperature.
-
8.
Precipitation variability. As for temperature.
-
9.
Archipelago plant richness. Plant richness relates to the productivity and niche diversity within an archipelago and therefore the support for higher trophic levels such as birds93. Hence, archipelago plant richness relates to the number of bird, and therefore extinct bird, species supported within an archipelago.
-
10.
Total area. Area is a major driver of island biodiversity31,84. Archipelagos with greater total area support larger (meta-)populations of bird species, which are thus less prone to extinction3,31. Archipelagos with greater area are also more likely to provide refuge from extinction pressures3,7. Archipelagos with greater area likely require more sampling effort to discover fossil specimens. Thus, archipelago area can affect both the number, and detection, of extinct birds.
-
11.
Standard deviation (SD) in area. SD area relates to the disparity in island size within an archipelago and is therefore associated with meta-population processes, such as source-sink dynamics94. In addition, SD area captures differences in structure between archipelagos, e.g., archipelagos with a few large islands vs many small islands.
-
12.
Native rodents. The impact of (human-introduced) exotic predators, particularly rats and mice, is reduced on islands that possess native rodents, likely because these avifaunas are less naïve to the effects of nest predation95,96. Thus, insufficient eco-evolutionary exposure (i.e., absence of native rodents) likely drives avian naïvety towards exotic predators95, and therefore relates to extinction risk. Approximately a third of the focal archipelagos have native rodents97 and hence likely reduced extinction risk from exotic predators.
-
13.
Human arrival. Archipelagos colonized long ago have had more time for fossil specimens to degrade, deteriorate and disappear98. In addition, the timing of peopling of an archipelago could affect the magnitude and structure of human pressures6,99,100. Indeed, human arrival could reflect changes in agricultural practice and land modification (widespread land clearance and domesticated plants are more often associated with later human settlements)99, changes in hunting preferences, technologies and efficiency (e.g., use of bows, arrows, and spears; fishing capabilities)100,101, changes in the number of human introductions and commensals (e.g., pigs, rats)2,5,6, and changes in human behaviour and culture (e.g., communication/language, long-distance trade)99,100.
-
14.
Endemic birds. The number of (remaining) endemic bird species likely relates to evolutionary processes, such as dispersal and speciation, and therefore to general avian diversity. Hence, the number of endemic birds likely relates to the number of extinct birds.
We quantified the number of endemic bird species as the number of extant and observed extinct (i.e., those extant at 1500 CE) endemic birds per archipelago. To identify endemic birds we used expert-based species’ range maps for 11,165 species102. We restricted species ranges to only include areas where they were listed as native or reintroduced (i.e., we excluded introduced or vagrant areas, etc.) and where their seasonality was listed as resident or breeding season (i.e., we excluded passage areas, etc.). For 344 seabird species, we restricted their ranges further to only breeding season areas (i.e., removing areas where seabirds are resident at sea). We defined seabirds taxonomically as groups where all species feed at sea, either nearshore or offshore (Spheniscifromes, Procellariiformes, Phaethontiformes, Pelecaniformes [family Pelecanidae], Suliformes [families Fregatidae, Phalacrocoracidae and Sulidae] and Charadriiformes [families Stercorariidae, Laridae and Alcidae]103). We then constructed convex hulls for the archipelagos and buffered them by 4000 m (the shortest distance between any archipelago and the mainland = 4200 m) to account for any spatial mismatches between the species’ ranges and the archipelagos, and to ensure that small islands and islets (<5 km2; small islands and islets might be occupied by satellite bird populations but are unlikely to be large enough to support endemic species10) were included. We then calculated the number of species endemic to the archipelago by summing the number of species that were found within each buffered archipelago convex hull (st_covered_by function from the sf package60).
Total area was strongly correlated with SD area (Pearson r = 0.91, df = 67, P < 0.001), so we included only total area in the regression models, but total area can be interpreted in light of their covariation. Isolation distance and surrounding landmass were also strongly correlated (Pearson r = −0.84, df = 67, P < 0.001), so we included only isolation distance in the models, but again this means that the estimated effect can be interpreted as the effect of isolation distance or surrounding landmass. This resulted in 12 suitable predictors (Supplementary Data 3).
Model of fossil extinct birds
First, we fit a linear model (LM) of log(fossil extinct birds + 1) as explained by the main effects of the 12 suitable predictors (Supplementary Data 3) and the interaction between total area and research effort. We transformed the predictors where it improved the distribution of the residuals and/or the stability of the variance (Supplementary Data 3). All predictors were centred and scaled to zero mean and unit variance, to improve the interpretability of regression coefficients104. All variance inflation factors were <7 indicating acceptable levels of multicollinearity105. We found no evidence of spatial autocorrelation (Moran’s tests for one to 20 neighbours). To quantify model uncertainty we simulated >1000 posterior estimates106, via the sim function in the arm package42. For similar applications see refs. 107, 108.
Transformation of count data for modelling purposes is often not recommended, due to the expected poorer performance of LMs compared to generalized linear models (GLMs)109. Hence, we also fit a quasipoisson (the data were overdispersed) GLM. Although a GLM might seem preferable compared to LM109, we were concerned that the strongly asymmetrical errors of the GLM had the potential to generate upwards bias; hence the creation of the LM as well.
To test the performance of the LM vs GLM we used leave-one-out cross-validation. We evaluated the predictive performance of the models (comparing the observed to the predicted numbers of fossil bird extinctions) with three goodness of fit measures using the gof function (hydroGOF package50): Spearman’s rank correlation, cross-validated R2, and percent bias. Spearman’s rank correlation measured the model’s ability to correctly rank archipelagos according to their number of extinctions. Cross-validated R2 gave the proportion of the variance of the observed extinctions that is predictable from the estimated extinctions, highlighting the model’s efficiency. Percent bias measured the average tendency of the estimated extinctions to be larger or smaller than the observed extinctions. Overall, the LM had better predictive performance than the GLM, with similar Spearman’s rank correlation (LM = 0.73; GLM = 0.74), greater R2 (LM = 0.42; GLM = 0.14), and similar absolute percent bias (LM = -35; GLM = +29). Crucially, the GLM showed positive percent bias—the GLM overestimates bird extinctions—likely driven by the strongly asymmetrical errors of the GLM. Overestimation is particularly concerning for extrapolative approaches (see below), as extrapolation can amplify these overestimations. For instance, after extrapolation, 23% of the GLM estimates were greater than 1000 undiscovered extinctions for a single archipelago, a biologically implausible estimate. Hence, the LM shows both better performance and more conservative extinction estimates, so we chose to take the LM forward for further analysis.
The LM explained a substantial proportion of the variance in fossil bird extinctions (multiple R2 = 0.68; adjusted R2 = 0.60). The strongest predictor was research effort (partial R2 = 0.26), followed by human arrival (partial R2 = 0.13), and then isolation distance, elevation, temperature and native rodents (all partial R2 = 0.03) (Supplementary Data 4). See below for further validation, performance testing, and constraints applied to the LM.
Extrapolation
Next, we used the LM to extrapolate the number of extinct birds for each archipelago (Fig. 4). We extrapolated based on the research effort required to produce a complete inventory of extinct bird species for Aotearoa New Zealand. In other words, how many extinct birds would we expect if all archipelagos had the same research effort as Aotearoa New Zealand needed to discover their complete avifauna? An advantage of this extrapolative approach is that the total estimated number of extinct birds per archipelago must necessarily be greater than or equal to the observed number (biologically realistic), whereas other approaches (e.g., prediction-based approaches) can estimate fewer extinctions than observed (i.e., negative extinctions; biologically impossible). In addition, for Aotearoa New Zealand the estimated value is necessarily equal to the observed value. Thus, our estimates are benchmarked on Aotearoa New Zealand. We used Aotearoa New Zealand as our reference archipelago as extinction rates here are not subject to the same uncertainties associated with incomplete detection elsewhere7,29. Specifically, Aotearoa New Zealand was only recently colonized by humans (ca. 1250 CE110), contains abundant, well‐preserved remains of the Holocene bird assemblage (all modern native species have been found as remains111) and is the only place in the world where the pre-human avifauna is believed to have been comprehensively reconstructed29,85,112,113,114,115. Thus, Aotearoa New Zealand is the only reasonable point of reference for these analyses. In addition, Aotearoa New Zealand had a research effort of 343 eligible publications by 2009 when all extinct species were described, while the next most researched archipelago, Hispaniola*, had a research effort of 270. Thus, the analyses could also be interpreted as extrapolating to the maximum observed research effort (Fig. 4). For each archipelago we extrapolated the total number of extinct birds using the LM and the difference in research effort (and the scaling of increasing research effort with increasing total area, fit as an interaction term in the model) as compared to Aotearoa New Zealand, with all else being equal. We calculated total number of extinct birds for each of the >1000 posterior estimates from the LM. We back-transformed the extrapolations to the natural scale and then subtracted the number of known fossil extinct birds from the total to obtain estimates of the number of undiscovered extinct birds. We converted our extrapolations to species numbers (integers) using probability-based rounding. For instance, the value of 1.7 had a 0.7 probability of being 2 extinctions and a 0.3 probability of being 1 extinction.
We have described our estimates of undiscovered extinct birds as extrapolations, as we extrapolated outside the observed data for research effort. Yet, for the remaining predictors in the model, the extrapolations generally fell within the range of observed data (i.e., interpolations; Supplementary Fig. 7).
Archipelago upper bounds
To test and ensure our estimates of the number of undiscovered bird extinctions were biologically realistic we calculated archipelago-specific upper bounds (Fig. 4).
We calculated upper bounds to constrain extremely high values—maximum estimates, above which predictions would be biologically unrealistic. We based our upper bounds on island biogeography theory, which states that continental bird diversity is greater than island bird diversity per unit area, due to the smaller system size (higher extinction) and greater geographical isolation (reduced immigration) of island systems31,38,73,116. We therefore quantified the maximum bird diversity of the archipelagos as if they were continental, producing estimates of maximum potential bird diversity.
To quantify the upper bounds, we first calculated the number of extant and observed extinct (i.e., extant at 1500 CE) bird species at 1° resolution117 across continental areas. We used range maps for 11,165 bird species, prepared the same as for the endemic birds (see above; except instead of using the st_covered_by function [a species’ range is wholly within an archipelago, i.e., endemic] we used st_intersects [any of the range intersects with the archipelago]60). We converted the range maps into a global presence-absence raster using a polygon-to-grid procedure (lets.presab function; letsR package52) and summed the number of bird species per grid cell to obtain bird species richness.
We then compiled five environmental predictors at 1° resolution equal-area projection across continental areas: elevation, temperature, precipitation, temperature variability, and precipitation variability (transformed where applicable; Supplementary Data 5). We extracted elevation from SRTM3080, while the bioclimatic predictors were extracted from WorldClim79. We excluded grid cells with missing predictor data as well as grid cells within the arctic or Antarctic circles, due to the transient dynamics of birds in these regions118— resulting in 10,600 grid cells.
We modelled bird species richness as explained by the five environmental predictors with a quasipoisson GLM. To reduce residual spatial autocorrelation we also included a residuals autocovariate term119, built from the residuals of an initial non-spatial GLM, and an optimized spatial neighbourhood structure (first-order neighbourhood was optimal based on Moran’s I). To quantify uncertainty we simulated >1000 posterior estimates106 from the spatial GLM, using the sim function (arm package42).
We then predicted the species richness of the 68 focal archipelagos (excluding Aotearoa New Zealand) as if they were continental (median predicted continental richness across grid cells per archipelago). We scaled the grid-based continental richness estimates to the total area of the archipelagos using the canonical species-area relationship81,82. We then subtracted the number of fossil extinctions, observed extinctions, and extant bird species of an archipelago from the continental richness to obtain estimates of the potential additional species that an archipelago could support (Fig. 4). Where continental richness fell below recorded species richness, we set the number of potential additional species to zero (i.e., recorded species richness is at or close to the predicted continental maximum). We used these estimates of potential additional species as upper bounds, above which predictions would be biologically unrealistic.
The upper bound spatial model explained a substantial proportion of the variance in continental bird species richness (McFadden pseudo-R2 = 0.86; Supplementary Data 6). The upper bound model also showed good cross-validated predictive performance (Spearman’s rank correlation = 0.99, R2 = 0.93, percent bias = 0.00).
Global lower bound
Note, we also calculated a global lower bound of the number of describable species by fitting a Bayesian logistic function to a species-description curve120 using Markov Chain Monte Carlo (MCMC) estimation. Virtually all of our estimates of undiscovered extinctions were greater than the lower bound, so, for simplicity, we do not take the lower bound further, but methodological details are provided here.
We calculated a global lower bound, a minimum estimate, below which predictions would be biologically unrealistic. To characterize the lower bound, we estimated the number of describable (those species that have been described from fossils and those that are likely to be described in the future, based on current efforts) fossil extinct bird species. This is necessarily a lower bound, as it does not account for those species that are unlikely to ever be discovered, e.g., species that have gone extinct without fossilization or species that require intensive research effort to be discovered3,7,10.
To quantify the number of describable species we fit a logistic function to a species-description curve120, based on fossil bird extinctions and the date that they were described72. We then modelled the number of described species at time t as samples from a normal density:
with mean µt, following a logistic function and variance estimated from the data. We fit the Bayesian model using Markov Chain Monte Carlo (MCMC) estimation, with a sampling frequency of 5000, a burn-in period of 100,000 iterations and flat unbounded priors, except for the midpoint, which we forced to be in the observed time window (1843 CE–2017 CE). We therefore assume that the description of ≥50% of all prehistoric extinct birds occurs sometime after 1843 CE and before 2017 CE, which is a reasonable assumption given that there are fewer species being described per author than previously, suggesting it is already harder to discover new species121,122. To get the number of describable species, we obtained posterior estimates of the maximum of the logistic (excluding the burn-in period); we then subtracted the number of described species (283 species) to estimate the number of undescribed species (Fig. 5).
Rejection sampling
We used rejection sampling to integrate the estimates of undiscovered extinct birds from the LM with the archipelago upper bounds (Fig. 4). We rejected estimates of undiscovered bird extinctions where the archipelago had a value greater than their respective upper bound estimate. We rejected relatively few (median = 0.2% rejected; mean = 12% rejected) estimates per archipelago based on the upper bounds (Supplementary Figs. 8, 9). Archipelagos with high rejection rates generally had similar recorded bird species richness to continental richness, i.e., recorded bird richness is at or close to maximum, e.g., Sado*, Derawan*, Yucatan* (Supplementary Fig. 8).
Alternative model of fossil extinct birds
As an alternative modelling approach to the LM model described above, we also implemented a white noise model with multiple predictors using a Bayesian framework, where all parameters are sampled from their posterior distribution using a MCMC algorithm. The mean of the white noise process (μx) changes as a linear function of multiple predictors (P) and their estimated effect sizes (e), following the linear formula:
We used the same predictors (n = 12) in the white noise linear model as for the LM above (Supplementary Data 3). We also modelled explicitly the interaction between the two predictors total area and research effort with its own effect size parameter (e(n+1) × P(total area) × P(research effort)), which was added to the linear formula above. The variation (noise) around the mean of the white noise process is captured by a separate parameter (σ), describing the standard deviation of the normal distribution that is centred in μx. During the MCMC we sampled the parameter estimates for μ0, σ, and the effect sizes (ei) for all predictors, including the interaction term. The parameter μ0 can be interpreted as describing the expected baseline number of extinctions if all predictor values were 0. As the likelihood formula for the MCMC we applied the probability density function of the normal distribution. As priors we applied for μ0 a uniform prior with the limits 0–200, for σ a gamma prior with a shape of 1 and a rate of 0.01, and a normal prior on all effect size parameters with a mean of 0 and an exponential hyperprior on the standard-deviation with a rate of 1.
The MCMC was run for 300,000,000 iterations, sampling every 500 iterations. Although the MCMC converged already after 1,000,000 iterations, we continued the MCMC for this many iterations to ensure sufficient parameter samples for the following rejection sampling (after excluding 10% as burn-in).
Again, when comparing the predicted to the observed bird extinctions for the archipelagos, the white noise linear model (Spearman’s rank correlation = 0.65; R2 = 0.59; percent bias = +11) outperformed an exponential version of the white noise model (Spearman’s rank correlation = 0.60; R2 = 0.07; percent bias = +207; some extreme estimates, i.e., >1000 undiscovered extinctions for a single archipelago).
We applied the same extrapolative approach as above to obtain estimates of undiscovered extinctions for the white noise linear model. We also applied the same process of rejection sampling. Again, we rejected relatively few (median = 0% rejected; mean = 11% rejected) estimates per archipelago based on the upper bounds.
The results for the two modelling approaches (LM and white noise linear model) show moderate correlation (Spearman’s rank correlation = 0.5), although the total estimates are very similar (LM = 789 total undiscovered extinctions; white noise linear model = 935 total undiscovered extinctions). In the end, we decided to only take the LM forward for further analysis, as it was more conservative in terms of estimated extinctions, easier to interpret, and seemed to capture more of the differences between archipelagos (greater variance). Still, the similarity of the total estimates between the two modelling approaches gave us confidence in the validity of our extinction estimates.
Comparison to previous data
To further test the validity of our estimations we also compared them to previous estimates across the Pacific4,7,10,12. Although direct comparison is problematic, due to differing geographic and/or species inclusion criteria (e.g., different definitions of the Pacific, and/or different bird inclusion—single island endemics, only landbirds, etc.). For our estimate, we included archipelagos to match the coverage of ref. 10, representing Micronesia, Melanesia and Polynesia.
Extinction rate
After calculating the number of bird extinctions across the globe, we estimated the timing and rate of extinctions. To calculate the extinction rate, we first estimated extinction dates for all species.
Fossil
For all fossil bird extinctions we estimated extinction dates based on truncated exponential decays, following first human arrival. Hence we assume that the most extinction-prone species will be lost rapidly after human arrival and the rate of extinction will then slow8. For all archipelagos, except Madagascar, we modelled 75% of extinctions to occur within 200 years of first human arrival (half-life of 100 years), based on the well-studied extinction chronology of Aotearoa New Zealand23,123,124. We truncated the exponential distributions so that 100% of extinctions occurred at the 90th quantile, i.e., 100% of extinctions occurred within 332 years of first human arrival—to prevent predicting extinction dates far outside the window of initial human impact. We therefore assume that the majority of initial human disturbance and human-driven extinctions are likely to have occurred within a few centuries23,123,124. This estimate is likely conservative across the focal archipelagos, as most are smaller in area than Aotearoa New Zealand; therefore, it is likely that all else being equal, humans would have colonized—and impacted—more rapidly across these smaller archipelagos.
Madagascar, on the other hand, is much larger than Aotearoa New Zealand (2.2 times larger) and is often described as an island-continent, as it is influenced by both island-like and continent-like processes125,126. We therefore applied a longer extinction window to Madagascar compared to the other archipelagos. We modelled 75% of extinctions to occur within 2000 years of first human arrival (half-life of 1000 years), based on the well-studied extinction chronology of North American Pleistocene mammals127. In other words, we assume that the extinction chronology of Madagascar is more similar to the timing of continental extinctions than to the timing of island extinctions (where Aotearoa New Zealand and North American Pleistocene mammals are the only available, suitable extinction chronologies). For Madagascar, we truncated the exponential distribution so that no extinctions occurred post-1950 (i.e., 100% of extinctions occurred at the 75th quantile).
There were also 98 continental fossil bird extinctions, which we assigned to their respective biogeographic realms128. For the Australasia, Nearctic and Neotropic realms, which have better established and later human arrival dates129,130, compared to the Indo-Malay and Palearctic131, we estimated extinction dates based on truncated exponential decays after first human arrival. We used a Clovis-culture human arrival date for the Nearctic130, even though there is evidence of earlier human arrival132,133, as no consensus has been reached among archaeologists about the date of initial human arrival to the Americas (and evidence of human arrival south of the ice sheets is fragmented), but all agree that human populations were distributed across the North American continent by the Clovis period134, implying greater human influence at this point. We modelled 75% of extinctions to occur within 2000 years of first human arrival (half-life of 1000 years), based on the extinction chronology of North American Pleistocene mammals127, and truncated the exponential distributions so that 100% of extinctions occurred at the 90th quantile.
When human arrival date was uncertain for a region (archipelago or continent), we randomly sampled (for each of the 1000 extinction estimates) from a uniform distribution based on the arrival date bounds.
For the Indo-Malay and Palearctic realms we randomly estimated extinction dates per species from uniform distributions between the minimum age of the fossil and 1500 CE, i.e., the species went extinct sometime after the fossil was deposited and before preserved written observational records began. We used a different approach for these realms as they have much earlier and more complex human arrival sequences131.
Observed
For observed extinctions we used the date of the last reliable or confirmed record22,72,135, where available. For 17 species, date of the last reliable or confirmed record was unavailable. For these species, we estimated extinction dates from truncated exponential distributions (half-life of 100 years), following first human arrival (see above).
Observed possibly extinct
We also accounted for possibly extinct species22. Species are only classified as Extinct (or Extinct in the Wild) on the IUCN Red List if there is no reasonable doubt that the last individual has died136. This strict definition of extinction aims to reduce the Romeo error—when conservationists give up on a species prematurely137; but failing to recognize extinctions leads to underestimates of extinction rates22,28. Thus, we incorporated 46 possibly extinct species that have been previously identified22. To incorporate the possibly extinct species, we used the reported mean extinction probability22. Then, for each of the 1000 extinction estimates, we randomly sampled from a binomial distribution based on the species’ probability of extinction—resulting in ≤46 additional extinctions per estimate. For these species we used the date of the last reliable or confirmed record as the extinction date22,135.
Overall extinction rate
We generated 1000 potential extinction chronologies for each of the archipelagos, accounting for variation in the number of bird extinctions (variation in the number of undiscovered extinctions and possibly extinct species), timing of human arrival (uncertainty associated with human arrival dates) and exponential decay (randomization of exponential extinction sequences). We calculated extinction rate as the number of extinctions at time t divided by the number extant at time t (extinctions/species/year) for the period 124,050 BCE–2019 CE, with 10,865 extant birds at present. We then calculated a rolling mean of yearly extinction rate with a window of 100 years. We summarized extinction rate as the mean across the 1000 extinction estimates.
We report extinction rate on the natural scale (extinctions/species/year), although others have advocated for the scale extinctions per million species years (E/MSY)8. Our extinctions/species/year estimates were therefore very roughly converted to E/MSY by multiplying by one million13, however see13 for a discussion of the issues with E/MSY, including converting per-year estimates of extinction to E/MSY and comparing E/MSY across different numbers of species. Thus, following ref. 13 we focus on extinctions/species/year estimates, while E/MSY is only used to put our estimates in the context of background extinction rates8,18,138.
Although the exact causes of bird extinctions are debated1,6,23,139, here we assume humans contributed (directly or indirectly) to all bird extinctions, supported by multiple lines of evidence2,3,5,6,85,124,140,141. There are few meaningful exceptions to this assumption—the Maltese Giant Swan (Cygnus falconeri) is one likely natural extinction142. Excluding these species from our analyses would lead to incorrect extinction rate estimates, plus the impact of these natural extinctions on the total extinctions is minimal; thus, we included them in all analyses.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
All prepared and processed data are available at https://zenodo.org/records/10014585 (doi: 10.5281/zenodo.10014585). Raw data on island characteristics (e.g., island area, island precipitation)66 are available from https://doi.org/10.5061/dryad.fv94v. Shapefiles of the archipelagos are available from the database of global administrative areas (GADM) https://gadm.org/data.html. Information on fossil extinct birds72 is available from https://doi.org/10.5061/dryad.s1rn8pk66. Raw information on native rodents97 is available from https://doi.org/10.5281/zenodo.1250504. Shapefiles of bird distributions are available from http://datazone.birdlife.org/species/requestdis. Global elevation data80 are available from https://srtm.csi.cgiar.org/srtmdata/, while WorldClim79 data are available from https://www.worldclim.org/data/worldclim21.html or via the getData function from the raster R package55. Bird extinction probabilities and the date of recent bird extinctions22 are available from https://data.mendeley.com/datasets/vvjhpmyxb4/1. Source data are also provided with this paper. Specifically, the source data underlying Figs. 2, 3 and 5, and Supplementary Figs. 4, 5, 7–9 are provided. Source data are provided with this paper.
Code availability
R code to reproduce the analyses is available at https://zenodo.org/records/10014585 (https://doi.org/10.5281/zenodo.10014585).
References
Steadman, D. W. Extinction & Biogeography Of Tropical Pacific Birds (University of Chicago Press, 2006). https://doi.org/10.1080/15564890701615429.
Szabo, J. K., Khwaja, N., Garnett, S. T. & Butchart, S. H. M. Global patterns and drivers of avian extinctions at the species and subspecies level. PLoS ONE 7, e47080 (2012).
Blackburn, T. M., Cassey, P., Duncan, R. P., Evans, K. L. & Gaston, K. J. Avian extinction and mammalian introductions on oceanic islands. Science 305, 1955–1958 (2004).
Steadman, D. W. Prehistoric extinctions of Pacific island birds: biodiversity meets zooarchaeology. Science 267, 1123–1131 (1995).
Boyer, A. G. Extinction patterns in the avifauna of the Hawaiian islands. Divers. Distrib. 14, 509–517 (2008).
Wood, J. R. et al. Island extinctions: processes, patterns, and potential for ecosystem restoration. Environ. Conserv. 44, 348–358 (2017).
Duncan, R. P., Boyer, A. G. & Blackburn, T. M. Magnitude and variation of prehistoric bird extinctions in the Pacific. Proc. Natl Acad. Sci. USA 110, 6436–6441 (2013).
Pimm, S., Raven, P., Peterson, A., Şekercioǧlu, Ç. H. & Ehrlich, P. R. Human impacts on the rates of recent, present, and future bird extinctions. Proc. Natl Acad. Sci. USA 103, 10941–10946 (2006).
Boehm, M. M. A. & Cronk, Q. C. B. Dark extinction: the problem of unknown historical extinctions. Biol. Lett. 17, 20210007 (2021).
Curnutt, J. & Pimm, S. L. How many bird species in Hawaii and the Central Pacific before first contact? Stud. Avian Biol. 22, 15–30 (2001).
Pimm, S. L. et al. The biodiversity of species and their rates of extinction, distribution, and protection. Science 344, 1246752 (2014).
Pimm, S. L., Moulton, M. P. & Justice, L. J. Bird extinctions in the central Pacific. Philos. Trans. R. Soc. Lond. B: Biol. Sci. 344, 27–33 (1994).
Monroe, M. J., Butchart, S. H. M., Mooers, A. O. & Bokma, F. The dynamics underlying avian extinction trajectories forecast a wave of extinctions. Biol. Lett. 15, 20190633 (2019).
Loehle, C. & Eschenbach, W. Historical bird and terrestrial mammal extinction rates and causes. Divers. Distrib. 18, 84–91 (2012).
Hawksworth, D. L. & Cowie, R. H. The discovery of historically extinct, but hitherto undescribed, species: an under-appreciated element in extinction-rate assessments. Biodivers. Conserv. 22, 2429–2432 (2013).
Tedesco, P. A. et al. Estimating how many undescribed species have gone extinct. Conserv. Biol. 25, 1360–1370 (2014).
Faurby, S. & Svenning, J. C. Historic and prehistoric human-driven extinctions have reshaped global mammal diversity patterns. Divers. Distrib. 21, 1155–1166 (2015).
Ceballos, G. et al. Accelerated modern human-induced species losses: entering the sixth mass extinction. Sci. Adv. 1, e1400253 (2015).
Cooke, R. S. C., Eigenbrod, F. & Bates, A. E. Projected losses of global mammal and bird ecological strategies. Nat. Commun. 10, 2279 (2019).
Valente, L., Etienne, R. S. & Garcia-R, J. C. Deep macroevolutionary impact of humans on New Zealand’s unique avifauna. Curr. Biol. 29, 2563–2569 (2019).
Heinen, J. H., van Loon, E. E., Hansen, D. M. & Kissling, W. D. Extinction-driven changes in frugivore communities on oceanic islands. Ecography 41, 1245–1255 (2018).
Butchart, S. H. M. et al. Which bird species have gone extinct? A novel quantitative classification approach. Biol. Conserv. 227, 9–18 (2018).
Holdaway, R. N. & Jacomb, C. Rapid extinction of the moas (Aves: Dinornithiformes): model, test, and implications. Science 287, 2250–2254 (2000).
Tilman, D., May, R. M., Lehman, C. L. & Nowak, M. A. Habitat destruction and the extinction debt. Nature 371, 65–66 (1994).
Triantis, K. A. et al. Extinction debt on oceanic islands. Ecography 33, 285–294 (2010).
Andermann, T., Faurby, S., Cooke, R., Silvestro, D. & Antonelli, A. iucn_sim: a new program to simulate future extinctions based on IUCN threat status. Ecography 44, 162–176 (2021).
Halley, J. M., Monokrousos, N., Mazaris, A. D., Newmark, W. D. & Vokou, D. Dynamics of extinction debt across five taxonomic groups. Nat. Commun. 7, 12283 (2016).
Akçakaya, H. R. et al. Inferring extinctions III: a cost-benefit framework for listing extinct species. Biol. Conserv. 214 (2017).
Holdaway, R. N., Worthy, T. H. & Tennyson, A. J. D. A working list of breeding bird species of the New Zealand region at first human contact. N. Z. J. Zool. 28, 119–187 (2001).
Livezey, B. Evolution of flightlessness in rails (Gruiformes: Rallidae): phylogenetic, ecomorphological, and ontogenetic perspectives. Ornithol. Monogr. 53, iii–654 (2003).
MacArthur, R. H. & Wilson, E. O. The Theory of Island Biogeography (Princeton University Press, 1967).
Isbell, F., Tilman, D., Polasky, S. & Loreau, M. The biodiversity-dependent ecosystem service debt. Ecol. Lett. 18, 119–134 (2015).
Şekercioğlu, Ç. H., Daily, G. C. & Ehrlich, P. R. Ecosystem consequences of bird declines. Proc. Natl Acad. Sci. 101, 18042–18047 (2004).
Bond, W. J. & Silander, J. A. Springs and wire plants: anachronistic defences against Madagascar’s extinct elephant birds. Proc. R. Soc. B: Biol. Sci. 274, 1985–1992 (2007).
Holdaway, R. N. Systematics and Palaeobiology of Haast’s eagle (Harpagornis moorei Haast, 1872) (Aves: Accipitridae). https://doi.org/10.26021/6893 (1991).
Holdaway, R. N. New Zealand’s Pre-Human Avifauna and Its Vulnerability. N. Z. J. Ecol. 12, 11–25 (1989).
Forshaw, J. & Knight, F. Vanished and Vanishing Parrots. (CSIRO Publishing, 2017).
Jønsson, K. A. & Holt, B. G. Islands contribute disproportionately high amounts of evolutionary diversity in passerine birds. Nat. Commun. 6, 8538 (2015).
R. Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2023).
ESRI. ArcGIS 10.4 (Environmental Systems Research Institute, 2015).
Python Software Foundation. Python (Python Software Foundation, 2015).
Gelman, A. & Su, Y.-S. arm: Data Analysis Using Regression and Multilevel/Hierarchical Models. R package version 1.13-1 (2022).
Robinson, D., Hayes, A. & Couch, S. broom: Convert Statistical Objects into Tidy Tibbles. R package version 1.0.4. (2023).
Wilke, C. O. cowplot: Streamlined Plot Theme and Plot Annotations for ‘ggplot2’. R package version 1.1.1. (2020).
Hartig, F. DHARMa: Residual Diagnostics for Hierarchical (Multi-Level / Mixed) Regression Models. R package version 0.4.6. (2022).
Wickham, H., François, R., Henry, L. & Müller, K. dplyr: A Grammar of Data Manipulation. R package version 1.1.2. (2023).
Tiedemann, F. gghalves: Compose Half-Half Plots Using Your Favourite Geoms. R package version 0.1.4. (2022).
Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag, 2016).
Meredith, M. & Kruschke, J. HDInterval: Highest (Posterior) Density Intervals. R package version 0.2.0. (2018).
Zambrano-Bigiarini, M. hydroGOF: Goodness-of-fit functions for comparison of simulated and observed hydrological time series (2020).
Long, J. A. jtools: Analysis and Presentation of Social Scientific Data. R package version 2.2.1. (2022).
Vilela, B. & Villalobos, F. letsR: a new R package for data handling and analysis in macroecology. Methods Ecol. Evol. 6, 1229–1234 (2015).
Venables, W. N. & Ripley, B. D. Modern Applied Statistics with S (Springer, 2002).
Wickham, H. & Henry, L. purrr: Functional Programming Tools. R package version 1.0.1 (2023).
Hijmans, R. J. raster: Geographic Data Analysis and Modeling. R package version 3.6-20 (2023).
Wickham, H., Hester, J. & Bryan, J. readr: Read Rectangular Text Data. R package version 2.1.4 (2023).
Zhang, D. rsq: R-Squared and Related Measures. R package version 2.5 (2022).
Wickham, H. & Seidel, D. scales: Scale Functions for Visualization. R package version 1.2.1 (2022).
Yu, G. scatterpie: Scatter Pie Plot. R package version 0.1.8 (2022).
Pebesma, E. Simple features for R: standardized support for spatial vector data. R. J. 10, 439–446 (2018).
Pebesma, E. J. & Bivand, R. S. Classes and methods for spatial data in R. 5 (2005).
Müller, K. & Wickham, H. tibble: Simple Data Frames. R package version 3.2.1 (2023).
Wickham, H., Vaughan, D. & Girlich, M. H. K. tidyr: Tidy Messy Data. R package version 1.3.0 (2023).
Zeileis, A. & Grothendieck, G. zoo: S3 Infrastructure. Regul. Irregul. Time Ser. J. Stat. Softw. 14, 1–27 (2005).
Johnson, T. H. & Stattersfield, A. J. A global review of island endemic birds. Ibis 132, 167–180 (1990).
Weigelt, P., Jetz, W. & Kreft, H. Bioclimatic and physical characterization of the world’s islands. Proc. Natl Acad. Sci. USA 110, 15307–15312 (2013).
Miller, K. G. The Phanerozoic record of global sea-level change. Science 310, 1293–1298 (2005).
Amante, C. & Eakins, B. W. ETOPO1—1 arc-minute global relief model: procedures, data sources and analysis. in NOAA technical memorandum NESDIS NGDC-24 19 (National Geophysical Data Center, 2009). https://doi.org/10.1594/PANGAEA.769615.
Ehlers, J., Hughes, P. L. & Gibbard, P. D. Quaternary Glaciations—extent and Chronology: A Closer Look, Vol. 15 (Elsevier, 2011).
de Oliveira Boléo, J. O descobrimento e os descobridores das ilhas de Cabo Verde (1939).
Peres, D. História dos descobrimentos Portugueses (1960).
Sayol, F., Steinbauer, M. J., Blackburn, T. M., Antonelli, A. & Faurby, S. Anthropogenic extinctions conceal widespread evolution of flightlessness in birds. Sci. Adv. 6, eabb6095 (2020).
Valente, L. et al. A simple dynamic model explains the diversity of island birds worldwide. Nature 579, 92–96 (2020).
Triantis, K. A., Economo, E. P., Guilhaumon, F. & Ricklefs, R. E. Diversity regulation at macro-scales: species richness on oceanic archipelagos. Glob. Ecol. Biogeogr. 24, 594–605 (2015).
Triantis, K. A., Whittaker, R. J., Fernández-Palacios, J. M. & Geist, D. J. Oceanic archipelagos: a perspective on the geodynamics and biogeography of the World’s smallest biotic provinces. Front. Biogeogr. 8, e29605 (2016).
Santos, A. M. C. et al. Are species-area relationships from entire archipelagos congruent with those of their constituent islands? Glob. Ecol. Biogeogr. 19, 527–540 (2010).
Gardner, E. E., Walker, S. E. & Gardner, L. I. Palaeoclimate, environmental factors, and bird body size: A multivariable analysis of avian fossil preservation. Earth-Sci. Rev. 162, 177–197 (2016).
Dirrigl, F. J., Diaz, A. & Chavez, R. Effects of soil shrink-swell on bird bones: An experimental taphonomy pilot study. Quat. Int. 626–627, 6–13 (2022).
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978 (2005).
Farr, T. G. et al. The shuttle radar topography mission. Rev. Geophys. 45, 1–33 (2007).
Preston, F. W. The canonical distribution of commonness and rarity: part I. Ecology 43, 185–215 (1962).
Rosenzweig, M. L. Species Diversity in Space and Time (Cambridge University Press, 1995).
Moher, D., Liberati, A., Tetzlaff, J. & Altman, D. G. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. Ann. Intern. Med. 151, 264–269 (2009).
Whittaker, R. J. & Fernández-Palacios, J. M. Island Biogeography: Ecology, Evolution, and Conservation (Oxford University Press, 2007).
Duncan, R. P. & Blackburn, T. M. Extinction and endemism in the New Zealand avifauna. Glob. Ecol. Biogeogr. 13, 509–517 (2004).
Milberg, P. & Tyrberg, T. Naïve birds and noble savages—a review of man-caused prehistoric extinctions of island birds. Ecography 16, 229–250 (1993).
Weigelt, P. & Kreft, H. Quantifying island isolation—insights from global patterns of insular plant species richness. Ecography 36, 417–429 (2013).
Jacobi, J. D. & Atkinson, C. T. in Our Living Resources: A Report to the Nation on the Distribution, Abundance, and Health of U.S. Plants, Animals, and Ecosystems (eds. LaRoe, E. T., Farris, G. S., Puckett, C. E., Doran, P. D. & Mac, M. J.) 376–381 (National Biological Service, 1995).
Meyer, C., Kreft, H., Guralnick, R. & Jetz, W. Global priorities for an effective information basis of biodiversity distributions. Nat. Commun. 6, 8221 (2015).
Moura, M. R. et al. Geographical and socioeconomic determinants of species discovery trends in a biodiversity hotspot. Biol. Conserv. 220, 237–244 (2018).
Moura, M. R. & Jetz, W. Shortfalls and opportunities in terrestrial vertebrate species discovery. Nat. Ecol. Evol. 5, 631–639 (2021).
Rolett, B. & Diamond, J. Environmental predictors of pre-European deforestation on Pacific islands. Nature 431, 443–446 (2004).
Cappelli, S. L., Pichon, N. A., Mannall, T. & Allan, E. Partitioning the effects of plant diversity on ecosystem functions at different trophic levels. Ecol. Monogr. 92, e1521 (2022).
Matthews, T. J. On the biogeography of habitat islands: the importance of matrix effects, noncore species, and source-sink dynamics. Q. Rev. Biol. 96, 73–104 (2021).
Anton, A., Geraldi, N. R., Ricciardi, A. & Dick, J. T. A. Global determinants of prey naiveté to exotic predators. Proc. R. Soc. B: Biol. Sci. 287, 20192978 (2020).
Harper, G. A. & Bunbury, N. Invasive rats on tropical islands: their population biology and impacts on native species. Glob. Glob. Ecol. Conserv. 3, 607–627 (2015).
Faurby, S. et al. PHYLACINE 1.2: the phylogenetic atlas of mammal macroecology. Ecology 99, 2626–2626 (2018).
Davis, P. G. & Briggs, D. E. G. The impact of decay and disarticulation on the preservation of fossil birds. PALAIOS 13, 3–13 (1998).
Boivin, N., Crassard, R. & Petraglia, M. Human Dispersal and Species Movement: From Prehistory to the Present (Cambridge University Press, 2017).
Louys, J. et al. No evidence for widespread island extinctions after Pleistocene hominin arrival. Proc. Natl Acad. Sci. 118, e2023005118 (2021).
Darimont, C. T. et al. Humanity’s diverse predatory niche and its ecological consequences. Commun. Biol. 6, 1–10 (2023).
BirdLife International and Handbook of the Birds of the World. Bird species distribution maps of the world. Version 2018.1. http://datazone.birdlife.org/species/requestdis (2018).
Schreiber, E. A. & Burger, J. Biology of Marine Birds (CRC Press, 2001).
Schielzeth, H. Simple means to improve the interpretability of regression coefficients. Methods Ecol. Evol. 1, 103–113 (2010).
Vittinghoff, E., Glidden, D. V., Shiboski, S. C. & McCulloch, C. E. Regression Methods in Biostatistics. (Springer US, 2012). https://doi.org/10.1007/978-1-4614-1353-0.
Gelman, A. & Hill, J. Data Analysis Using Regression and Multilevel/hierarchical Models (Cambridge University Press, 2006). https://doi.org/10.1017/cbo9780511790942.
Nicolaus, M., Tinbergen, J. M., Ubels, R., Both, C. & Dingemanse, N. J. Density fluctuations represent a key process maintaining personality variation in a wild passerine bird. Ecol. Lett. 19, 478–486 (2016).
Sanders, D., Sutter, L. & van Veen, F. J. F. The loss of indirect interactions leads to cascading extinctions of carnivores. Ecol. Lett. 16, 664–669 (2013).
O’Hara, R. B. & Kotze, D. J. Do not log-transform count data. Methods Ecol. Evol. 1, 118–122 (2010).
Wilmshurst, J. M., Hunt, T. L., Lipo, C. P. & Anderson, A. J. High-precision radiocarbon dating shows recent and rapid initial human colonization of East Polynesia. Proc. Natl Acad. Sci. USA 108, 1815–1820 (2011).
Worthy, T. H. & Holdaway, R. N. The Lost World of the Moa: Prehistoric Life of New Zealand (Indiana University Press, 2002). https://doi.org/10.1671/0272-4634(2003)023[0475:TLWOTM]2.0.
Atkinson, I. A. E. & Millener, P. R. An ornithological glimpse into New Zealand’s pre-human past. Acta XX Congr. Int. Ornithol. 1, 129–191 (1991).
Worthy, T. H. & Holdaway, R. N. Quaternary fossil faunas from caves in the Punakaiki area, West Coast, South, New Zealand. J. R. Soc. N. Z. 23, 147–254 (1993).
Worthy, T. H. Quaternary fossil fauna of South Canterbury, South Island, New Zealand. J. R. Soc. N. Z. 27, 67–162 (1997).
Worthy, T. H. Quaternary fossil faunas of Otago, South Island, New Zealand. J. R. Soc. N. Z. 28, 421–521 (1998).
MacArthur, R. H. & Wilson, E. O. An equilibrium theory of insular zoogeography. Evolution 17, 373–387 (1963).
Hurlbert, A. H. & Jetz, W. Species richness, hotspots, and the scale dependence of range maps in ecology and conservation. Proc. Natl Acad. Sci. USA 104, 13384–9 (2007).
Kirby, J. S. et al. Key conservation issues for migratory land- and waterbird species on the world’s major flyways. Bird. Conserv. Int. 18, S49–S73 (2008).
Crase, B., Liedloff, A. C. & Wintle, B. A. A new method for dealing with residual spatial autocorrelation in species distribution models. Ecography 35, 879–888 (2012).
Benton, M. The completeness of the fossil record. Significance 6, 117–121 (2009).
Costello, M. J., Wilson, S. & Houlding, B. Predicting total global species richness using rates of species description and estimates of taxonomic effort. Syst. Biol. 61, 871 (2012).
Joppa, L. N., Roberts, D. L. & Pimm, S. L. The population ecology and social behaviour of taxonomists. Trends Ecol. Evol. 26, 551–553 (2011).
Tennyson, A. J. D. & Martinson, P. Extinct Birds of New Zealand (Te Papa Press, 2006).
Perry, G. L. W., Wheeler, A. B., Wood, J. R. & Wilmshurst, J. M. A high-precision chronology for the rapid extinction of New Zealand moa (Aves, Dinornithiformes). Quat. Sci. Rev. 105, 126–135 (2014).
Vences, M., Wollenberg, K. C., Vieites, D. R. & Lees, D. C. Madagascar as a model region of species diversification. Trends Ecol. Evol. 24, 456–465 (2009).
Wit, M. Jde Madagascar: heads it’s a continent, tails it’s an island. Annu. Rev. Earth Planet. Sci. 31, 213–248 (2003).
Faith, J. T. & Surovell, T. A. Synchronous extinction of North America’s Pleistocene mammals. Proc. Natl Acad. Sci. USA 106, 20641–20645 (2009).
Olson, D. M. et al. Terrestrial ecoregions of the worlds: a new map of life on Earth. BioScience 51, 933–938 (2001).
Clarkson, C. et al. Human occupation of northern Australia by 65,000 years ago. Nature 547, 306–310 (2017).
Waters, M. R. et al. The Buttermilk Creek complex and the origins of Clovis at the Debra L. Friedkin site, Texas. Science 331, 1599–1603 (2011).
Nielsen, R. et al. Tracing the peopling of the world through genomics. Nature 541, 302–310 (2017).
Shillito, L.-M. et al. Pre-Clovis occupation of the Americas identified by human fecal biomarkers in coprolites from Paisley Caves, Oregon. Sci. Adv. 6, eaba6404 (2020).
Davis, L. G. et al. Late Upper Paleolithic occupation at Cooper’s Ferry, Idaho, USA, ~16,000 years ago. Science 365, 891–897 (2019).
Surovell, T. A. et al. Late date of human arrival to North America: Continental scale differences in stratigraphic integrity of pre-13,000 BP archaeological sites. PLoS ONE 17, e0264092 (2022).
IUCN. The IUCN Red List of Threatened Species. Version 2019-3. https://www.iucnredlist.org (2020).
IUCN. IUCN Red List Categories and Criteria Version 3.1 (2001).
Collar, N. J. Extinction by assumption; or, the Romeo error on Cebu. Oryx 32, 239–244 (1998).
Barnosky, A. D. et al. Has the Earth’s sixth mass extinction already arrived? Nature 471, 51–57 (2011).
Grayson, D. K. The archaeological record of human impacts on animal populations. J. World Prehistory 15, 1–68 (2001).
Duncan, R. P. & Blackburn, T. M. Causes of extinction in island birds. Anim. Conserv. 10, 149–150 (2007).
Latham, A. D. M. et al. A refined model of body mass and population density in flightless birds reconciles extreme bimodal population estimates for extinct moa. Ecography 43, 353–364 (2020).
Thake, M. A. The biogeography of the Maltese Islands, illustrated by the Clausiliidae. J. Biogeogr. 12, 269–287 (1985).
Acknowledgements
The authors thank P. Weigelt for providing environmental data for Norfolk Island, as well as D. Silvestro and R. Duncan for help with the analyses. We also thank R. Smith and T. Matthews for comments on the manuscript. F.S. received support from “la Caixa” Foundation (ID 100010434) through a Junior Leader program (fellowship code LCF/BQ/PI23/11970019). T.A. received financial support from a SciLifeLab & Wallenberg Data Driven Life Science Program fellowship grant (KAW 2020.0239). M.J.S. acknowledges support by the Deutsche Forschungsgemeinschaft (STE 2360/2-1 embedded in FOR 2332) as well as ERC grant 741413 Humans on Planet Earth (HOPE). A.A. acknowledges financial support from the Swedish Research Council (2019-05191), the Swedish Foundation for Strategic Environmental Research MISTRA (Project BioPath), and the Kew Foundation. S.F. (and R.C. via S.F.) received support from the Swedish Research Council (2017-03862) and Carl Tryggers Stiftelse, project number CTS 18:105.
Author information
Authors and Affiliations
Contributions
The project was conceptualized by S.F., R.C., F.S., M.J.S. and T.M.B. Data were prepared by R.C. and F.S., and the data were analysed by R.C and T.A. The first version of the manuscript was written by R.C. and was constructively reviewed and edited by R.C., F.S., T.A., T.M.B., M.J.S., A.A. and S.F.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Julia Heinen, Stuart Pimm and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Cooke, R., Sayol, F., Andermann, T. et al. Undiscovered bird extinctions obscure the true magnitude of human-driven extinction waves. Nat Commun 14, 8116 (2023). https://doi.org/10.1038/s41467-023-43445-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-43445-2