Hostname: page-component-78c5997874-g7gxr Total loading time: 0 Render date: 2024-11-05T04:01:54.025Z Has data issue: false hasContentIssue false

Priority areas for vulture conservation in the Horn of Africa largely fall outside the protected area network

Published online by Cambridge University Press:  21 July 2021

EVAN R. BUECHLEY*
Affiliation:
HawkWatch International, 2240 South 900 East, Salt Lake City, UT 84106, USA. School of Biological Sciences, The University of Utah, 257 South 1400 East, Salt Lake City, UT 84112, USA. Migratory Bird Center, Smithsonian Conservation Biology Institute, National Zoological Park, MRC 5503, Washington, DC 20013-7012, USA.
MARCO GIRARDELLO
Affiliation:
cE3c – Centre for Ecology, Evolution and Environmental Changes / Azorean Biodiversity Group and Universidade dos Açores – Faculty of Agriculture and Environment, Rua Capitão João d’Ávila, São Pedro, PT-9700-042 Angra do Heroísmo, Terceira, Açores, Portugal.
ANDREA SANTANGELI
Affiliation:
The Helsinki Lab of Ornithology, Finnish Museum of Natural History, FI-00014, University of Helsinki, Finland. FitzPatrick Institute of African Ornithology, DST-NRF Centre of Excellence, University of Cape Town, Cape Town, South Africa.
ALAZAR DAKA RUFFO
Affiliation:
Addis Ababa University, Faculty of Natural Science, Department of Zoological Sciences, Ethiopia.
GIRMA AYALEW
Affiliation:
Ethiopian Wildlife Conservation Authority, Addis Ababa, Ethiopia.
YILMA D. ABEBE
Affiliation:
P.O. Box 18112, Addis Ababa, Ethiopia.
DAVID R. BARBER
Affiliation:
Acopian Center for Conservation Learning, Hawk Mountain Sanctuary, 410 Summer Valley Road, Orwigsburg, PA 17961, USA.
RALPH BUIJ
Affiliation:
The Peregrine Fund, 5668 West Flying Hawk Lane Boise, ID 83709, USA. Wageningen University & Research, Droevendaalsesteeg 3A, 6708 PB Wageningen, The Netherlands.
KEITH BILDSTEIN
Affiliation:
116 Village Drive, Blandon, PA, 19510, USA.
BRUKTAWIT ABDU MAHAMUED
Affiliation:
Kotebe Metropolitan University, Addis Ababa, Ethiopia.
MONTAGUE H.C. NEATE-CLEGG
Affiliation:
School of Biological Sciences, The University of Utah, 257 South 1400 East, Salt Lake City, UT 84112, USA.
DARCY OGADA
Affiliation:
The Peregrine Fund, 5668 West Flying Hawk Lane Boise, ID 83709, USA. National Museums of Kenya, Ornithology Section, P.O. Box 40658- 00100, Nairobi, Kenya.
PETER P. MARRA
Affiliation:
Department of Biology and McCourt School of Public Policy, Georgetown University, 37th and O Streets NW, Washington, DC 20057, USA.
T. SCOTT SILLETT
Affiliation:
School of Biological Sciences, The University of Utah, 257 South 1400 East, Salt Lake City, UT 84112, USA.
JEAN-MARC THIOLLAY
Affiliation:
2 rue de la Rivière, 10220 Rouilly Sacey, France.
MARTIN WIKELSKI
Affiliation:
Department of Migration, Max Planck Institute of Animal Behavior, 78315 Radolfzell, Germany. Centre for the Advanced Study of Collective Behaviour, University of Konstanz, 78457 Konstanz, Germany.
PETER YAWORSKY
Affiliation:
Department of Anthropology, University of Utah, 260 Central Campus Drive, Salt Lake City, Utah, 84112, USA. Archaeological Center, University of Utah, 260 Central Campus Drive, Salt Lake City, Utah, 84112, USA.
ÇAĞAN H. ŞEKERCIOĞLU
Affiliation:
School of Biological Sciences, The University of Utah, 257 South 1400 East, Salt Lake City, UT 84112, USA. Department of Molecular Biology and Genetics, Koç University, Istanbul, Turkey. KuzeyDoğa Derneği, Ortakapı Mah. Şehit Yusuf Bey Cad. No:69 Kat:1 36100 Merkez, Kars, Turkey.
*
*Author for correspondence; email: [email protected]
Rights & Permissions [Opens in a new window]

Summary

Vulture populations are in severe decline across Africa and prioritization of geographic areas for their conservation is urgently needed. To do so, we compiled three independent datasets on vulture occurrence from road-surveys, GPS-tracking, and citizen science (eBird), and used maximum entropy to build ensemble species distribution models (SDMs). We then identified spatial vulture conservation priorities in Ethiopia, a stronghold for vultures in Africa, while accounting for uncertainty in our predictions. We were able to build robust distribution models for five vulture species across the entirety of Ethiopia, including three Critically Endangered, one Endangered, and one Near Threatened species. We show that priorities occur in the highlands of Ethiopia, which provide particularly important habitat for Bearded Gypaetus barbatus, Hooded Necrosyrtes monachus, Rüppell’s Gyps rüppelli and White-backed Gyps africanus Vultures, as well as the lowlands of north-eastern Ethiopia, which are particularly valuable for the Egyptian Vulture Neophron percnopterus. One-third of the core distribution of the Egyptian Vulture was protected, followed by the White-backed Vulture at one-sixth, and all other species at one-tenth. Overall, only about one-fifth of vulture priority areas were protected. Given that there is limited protection of priority areas and that vultures range widely, we argue that measures of broad spatial and legislative scope will be necessary to address drivers of vulture declines, including poisoning, energy infrastructure, and climate change, while considering the local social context and aiding sustainable development.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2021. Published by Cambridge University Press on behalf of BirdLife International

Introduction

Current species extinction rates are an estimated 1,000 times over the background rate and are expected to increase with growing anthropogenic pressures worldwide (Barnosky et al. Reference Barnosky, Matzke, Tomiya, Wogan, Swartz, Quental, Marshall, McGuire, Lindsey, Maguire, Mersey and Ferrer2011). Loss of species and wildlife abundance is compromising ecological processes, reducing ecosystem services, and directly affecting humans (Şekercioğlu Reference Şekercioğlu, Sodhi and Ehrlich2010, Cadotte et al. Reference Cadotte, Carscadden and Mirotchnick2011). Tropical biodiversity is generally less studied than temperate biodiversity (Trimble and van Aarde Reference Trimble and van Aarde2012) and the tropics are expected to experience the greatest biodiversity losses in the 21st century (Alroy Reference Alroy2017). Overcoming this trend will require international investment and local community support (Ghosh-Harihar et al. Reference Ghosh-Harihar, An, Athreya, Borthakur, Chanchani, Chetry, Datta, Harihar, Karanth, Mariyam, Mohan, Onial, Ramakrishnan, Robin, Saxena, Shahabuddin, Thatte, Vijay, Wacker, Mathur, Pimm and Price2019). As biodiversity and threats are unevenly distributed, limited research and conservation investments should be strategically targeted to maximize effectiveness (Brooks et al. Reference Brooks, Mittermeier, Fonseca, Gerlach, Hoffman, Lamoreux, Mittermeier, Pilgrim and Rodrigues2006) and to fill threat gaps (Joppa et al. Reference Joppa, O’Connor, Visconti, Smith, Geldmann, Hoffmann, Watson, Butchart, Virah-Sawmy, Halpern, Ahmed, Balmford, Sutherland, Harfoot, Hilton-Taylor, Foden, Minin, Pagad, Genovesi, Hutton and Burgess2016).

As scavengers, vultures constitute the most endangered functional guild of birds and their populations are in severe decline worldwide (Buechley and Şekercioğlu Reference Buechley and Şekercioğlu2016a,Reference Buechley and Şekercioğlub). They are a top conservation priority because they are highly threatened and provide critical ecological functions by quickly consuming carrion, contributing to nutrient cycling and regulation of problematic facultative scavengers (Ogada et al. 2012) and, potentially, controlling disease (Buechley and Şekercioğlu Reference Buechley and Şekercioğlu2016a, Devault et al. Reference Devault, Beasley, Olson, Moleón, Carrete, Margalida, Sánchez-zapata, Sekercioglu, Wenny and Whelan2016, Plaza et al. Reference Plaza, Blanco and Lambertucci2020). Vultures face a range of threats, including most notably, poisoning, toxic veterinary drugs, loss of food availability, and collision and electrocution on energy infrastructure (Buechley and Şekercioğlu Reference Buechley and Şekercioğlu2016a, Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017, Plaza et al. Reference Plaza, Martínez-López and Lambertucci2019). Given their broad distributions, large individual home ranges, and long-distance nomadic and migratory movements, conservation of vultures is challenging (e.g. Runge et al. Reference Runge, Martin, Possingham, Willis and Fuller2014, Reference Runge, Watson, Butchart, Hanson, Possingham and Fuller2015). Prioritizing geographic areas for conservation actions is needed (Mukherjee et al. Reference Mukherjee, Galligan, Prakash, Paudel, Khan, Prakash, Ranade, Shastri, Dave, Donald and Bowden2014, Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2019a). Indeed, the recently published Multi-Species Action Plan to Conserve African-Eurasian Vultures (Vulture MSAP), which is based on extensive expert input, provides a roadmap for vulture conservation across Africa and Eurasia, highlighting the need to identify threats and prioritize conservation actions at different spatial scales (Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017).

To identify spatial conservation priorities, it is essential to know the distribution and habitat use of species. Species distribution models (SDMs) are an important tool used to identify priority habitats and to forecast anthropogenic effects on species (Guisan and Thuiller Reference Guisan and Thuiller2005, Aryal et al. Reference Aryal, Shrestha, Ji, Ale, Shrestha, Ingty, Maraseni, Cockfield and Raubenheimer2016). Combined with information on protected areas and existing threats to species, SDMs can help identify conservation priorities (Ferraz et al. Reference Ferraz, De Barros, De Siqueira, Alexandrino, Da Luz and Do Couto2012, Evans et al. Reference Evans, Asner and Goossens2018). However, the accuracy of SDMs is limited by the availability of both occurrence (i.e. where a species occurs) and predictor variables (i.e. spatial layers relevant to model a species’ distribution) (Fletcher Jr et al. Reference Fletcher, Hefley, Robertson, Zuckerberg, McCleery and and Dorazio2018), both of which can be limited in the Global South. Targeted survey data are costly and time consuming, and thus generally scarce. Bio-logging (e.g. GPS tracking) can be used to model species’ distributions (e.g. Coxen et al. Reference Coxen, Frey, Carleton and Collins2017), but can suffer from small sample sizes. Citizen science data (e.g. eBird; Sullivan et al. Reference Sullivan, Aycrigg, Barry, Bonney, Bruns, Cooper, Damoulas, Dhondt, Dietterich, Farnsworth, Fink, Fitzpatrick, Fredericks, Gerbracht, Gomes, Hochachka, Iliff, Lagoze, La Sorte, Merrifield, Morris, Phillips, Reynolds, Rodewald, Rosenberg, Trautmann, Wiggins, Winkler, Wong, Wood, Yu and Kelling2014) are being increasingly used to fill information gaps and have been shown to be useful for modelling bird distributions and population trends (Horns et al. Reference Horns, Adler and Şekercioğlu2018, Fink et al. 2019), but have spatial and observer biases (Fletcher Jr et al. Reference Fletcher, Hefley, Robertson, Zuckerberg, McCleery and and Dorazio2018, Horns et al. Reference Horns, Adler and Şekercioğlu2018, Fink et al. 2019, Neate-Clegg et al. Reference Neate-Clegg, Horns, Adler, Kemahlı Aytekin and Şekercioğlu2020). Combining multiple data sources to model species’ distributions can help increase predictive accuracy and account for sampling biases (Fletcher Jr et al. Reference Fletcher, Hefley, Robertson, Zuckerberg, McCleery and and Dorazio2018, Miller et al. Reference Miller, Pacifici, Sanderlin and Reich2019).

Here we used an SDM framework to identify spatial conservation priorities for vultures in Ethiopia, a high-priority country for vulture conservation (Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2019a). Our aims were to: 1) model the distribution of each vulture species; 2) identify vulture conservation priority areas; and 3) assess how well protected vulture priority areas are. We incorporated three independent datasets on vulture occurrence (road-surveys, GPS-tracking, and citizen science) to build robust predictions on vulture space use. Our results provide a framework for where to prioritize conservation work on endangered vultures in the Horn of Africa.

Methods

Study location and species

The study took place in Ethiopia, a biologically diverse and little-studied country, which supports one of the most species rich and abundant vulture communities worldwide (Mundy et al. Reference Mundy, Bunchart, Ledger and Piper1993, Buechley et al. Reference Buechley, Santangeli, Girardello, Neate-Clegg, Oleyar, McClure and Şekercioğlu2019). It is also the second most populous country in Africa and has a high human population growth rate (2.46% per year) (World Bank 2019b). Despite being one of the poorest countries in the world, Ethiopia is developing rapidly (World Bank 2019a), and there is severe pressure on natural ecosystems. Seven vulture species reside in Ethiopia: Bearded Gypaetus barbatus (‘Near Threatened’), Egyptian Neophron percnopterus (‘Endangered’), Hooded Necrosyrtes monachus (‘Critically Endangered’), Lappet-faced Torgos tracheliotos (‘Endangered’), Rüppell’s Gyps rueppelli (‘Critically Endangered’), White-backed Gyps africanus (‘Critically Endangered’), and White-headed Trigonoceps occipitalis (‘Critically Endangered’) (Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017). The largest African populations of Bearded and Egyptian Vultures are thought to occur in Ethiopia (Arkumarev et al. Reference Arkumarev, Dobrev and Abebe2014, Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017).

Species occurrence data

Vulture occurrence data were acquired from three independent sources: road-surveys, satellite telemetry, and citizen science. Road-surveys took place in Ethiopia from 2010 to 2018, whereby two experienced raptor biologists completed 10,857 km of surveys, while driving at speeds <60 km per hour during daylight hours, identifying perched and flying vultures within 1 km of the road. Citizen science observations of vultures were downloaded via the comprehensive eBird Basic Dataset (Sullivan et al. Reference Sullivan, Aycrigg, Barry, Bonney, Bruns, Cooper, Damoulas, Dhondt, Dietterich, Farnsworth, Fink, Fitzpatrick, Fredericks, Gerbracht, Gomes, Hochachka, Iliff, Lagoze, La Sorte, Merrifield, Morris, Phillips, Reynolds, Rodewald, Rosenberg, Trautmann, Wiggins, Winkler, Wong, Wood, Yu and Kelling2014) for Ethiopia, and censored to minimise sampling and observer bias, following specific recommendations and code provided by Johnston et al. (Reference Johnston, Hochachka, Strimas-Mackey, Gutierrez, Robinson, Miller, Auer, Kelling and Fink2019). Satellite tracking data were collected by deploying 34 solar-powered GPS-transmitters on vultures in or ranging through Ethiopia from 2012 to 2018, including 15 Egyptian, eight White-backed, six Hooded, three Ruppell’s, and two Lappet-faced vultures. See Figure S1 in the online supplementary material for an overview of occurrence data. Further details on acquisition and treatment of occurrence data are in Appendix S1.

Each dataset had strengths and weaknesses. For example, road-surveys were designed to have broad geographic coverage and high data quality, but were limited by the extent of the road network and security concerns inhibiting access to some areas. Further, road-surveys took place only during the dry season, and seasonal differences in vulture ranging behaviour has been shown elsewhere in Africa (e.g. Kendall et al. Reference Kendall, Virani, Hopcraft, Bildstein and Rubenstein2014), indicating that these surveys may not have captured each species’ full environmental niche. In contrast, satellite-tracking provided a picture of the habitat use of individual vultures throughout the full annual cycle, but had a limited sample of individuals. Meanwhile, citizen science data had broad geographic and temporal coverage, but were spatially biased towards urban centres and birding “hotspots”, and may have suffered from lower data quality due to potential species misidentification (e.g. Fink et al. 2019). We recognize that each of the data types was therefore capturing different aspects of the species’ ecological niche and at somewhat different spatial and temporal scales, and we see that as a net benefit of our approach. By incorporating multiple data types, we maximized the number of species for which we had sufficient data to model distributions, while also incorporating data that captured different aspects of each species’ environmental niche.

All analyses were completed in R (Version 3.1.0; R Core Team 2019). For eBird data, occurrence points were observations for each species spatially rarefied to a 1-km distance, while the background sample was the location of all complete checklists (Coxen et al. Reference Coxen, Frey, Carleton and Collins2017). For road-survey data, we buffered road transects by 1 km on both sides of the road, overlaid a 1-km2 raster grid over the survey area, and aggregated observations for each species within each pixel, such that any pixel where a given vulture species was sighted was given a value of 1 (occurrence), whereas pixels where that vulture species was not encountered were given a value of 0 (background). For satellite tracking data, we calculated the 95% minimum convex polygon (MCP) cumulative home range (‘adehabitatHR’ package; Clemente Calenge Reference Calenge2011) of all tracked individuals of each species (with one exception: we calculated separate MCPs for two Hooded Vulture populations that did not overlap from transmitter deployments in northern and southern Ethiopia). Tracked vulture locations within each species’ MCP were used as occurrence points, while background points were drawn from a systematic sample of 10,000 evenly spaced points within the MCP for each species (Benson Reference Benson2013).

Environmental variables

We compiled environmental data that were expected to influence vulture use, including habitat type (European Space Agency 2017), elevation and ruggedness (Robinson et al. Reference Robinson, Regetz and Guralnick2014), human footprint (Venter et al. Reference Venter, Sanderson, Magrach, Allan, Beher, Jones, Possingham, Laurance, Wood, Fekete, Levy and Watson2016), climate (annually averaged temperature, precipitation, wind, solar radiation; Fick and Hijmans Reference Fick and Hijmans2017), latitude and longitude. Prior to modelling, we assessed correlation of predictor variables by creating a correlation matrix using ‘corrplot’ (Wei and Simko Reference Wei and Simko2013), and reviewed the variance inflation factors using the ‘usdm’ package (Naimi 2015). We used a cutoff of |r| = 0.60 as an indication of strong co-linearity (Crandall et al. 2015). Elevation and temperature were inversely correlated (|r| = -1) and we thus excluded elevation because temperature was expected to have a more direct biological effect. Wind, precipitation, and longitude were highly correlated with other variables and were removed, leaving us with 12 predictor variables: latitude, ruggedness, human footprint, distance to cropland, distance to desert, distance to forest, distance to grassland, distance to shrubland, distance to urban, categorical landcover class, temperature, and solar radiation (Figure S3). Further details on processing of predictor variables can be found in Appendix S1.

Species distribution models

We modelled the distributions of seven vulture species in Ethiopia using maximum entropy (Maxent; Phillips et al. Reference Phillips, Anderson and Schapire2006). Maxent is the most widely used species distribution modeling technique (Elith et al. Reference Elith, Phillips, Hastie, Dudík, Chee and Yates2011). To maximize the accuracy of Maxent models, it is important to account for sampling bias of occurrence points, to select the appropriate background sample, and to select the best regularization multiplier (Merow et al. Reference Merow, Smith and Silander2013). To address these concerns, we spatially thinned each dataset, carefully selected the background sample, and tested a range of regularization multipliers. Maxent was run using the package ‘ENMeval’ (Muscarella et al. Reference Muscarella, Galante, Soley-Guardia, Boria, Kass, Uriarte and Anderson2014), using “algorithm = ‘maxent.jar’”, which interfaces with the standalone Maxent program (Version 3.4.1; Phillips and Dudík Reference Phillips and Dudík2020). For each data type and species combination, we ran a suite of models with a range of regularization multiplier values (from 1 to 4, by 1) and feature class combinations (L, LQ, H, LQH, LQHP, LQHPT; where L = linear, Q = quadratic, H = hinge, P = product and T = threshold), or up to 24 models each. To assess the predictive capacity of each model, we used k-fold internal cross validation (Merow et al. Reference Merow, Smith and Silander2013) with four folds (we used 75% of the data to train the model and 25% of the data to test the model), using the ‘block’ data partitioning method, which spatially partitions data into four bins by the lines of latitude and longitude that divide occurrence localities as equally as possible. We selected the top model for each species and data type using Akaike information criteria corrected for small sample sizes (AICc) (Burnham and Anderson Reference Burnham and Anderson2004, Muscarella et al. Reference Muscarella, Galante, Soley-Guardia, Boria, Kass, Uriarte and Anderson2014). Predictive accuracy of top models was assessed from the AUCtest score, which measures the model’s ability to discriminate between conditions at occurrence versus background locations (Muscarella et al. Reference Muscarella, Galante, Soley-Guardia, Boria, Kass, Uriarte and Anderson2014). We also report ORMTP, which is the average omission rate of the occurrence records at the minimum training presence (MTP) threshold, where MTP represents an estimate of species habitat suitability (Muscarella et al. Reference Muscarella, Galante, Soley-Guardia, Boria, Kass, Uriarte and Anderson2014, Taylor et al. Reference Taylor, Hafen, Holley, González and Long2020). We deemed models with AUCtest values >0.60 to be informative (Randin et al. Reference Randin, Dirnböck, Dullinger, Zimmermann, Zappa and Guisan2006).

Next, we compared the top predictions for each species and data type, using Hellinger’s I statistic. Hellinger’s I values range from 0 (no overlap) to 1 (complete overlap), quantifying niche overlap over geographic space while making no assumptions about species density (Warren et al. Reference Warren, Glor and Turelli2008).

Then for each species we created an ensemble species distribution model. To do so, we weighted the top species-data type model by its respective AUC test score and then summed them. We mapped the final predicted ensemble distribution of each species across Ethiopia, interpreted as relative habitat suitability (Merow et al. Reference Merow, Smith and Silander2013).

To estimate uncertainty in our predictions, we took 10 bootstrap samples of the training datasets and ran predictions for each species and data type combination on each bootstrapped sample. We then calculated the standard deviation of the predicted values for each species and data type combination and derived an overall map of the uncertainty in the modelled distribution for each species by taking the average of the standard deviation values across each data type.

We identified the core distribution of each species within Ethiopia as the top 30% of its predicted distribution areas. While the 30% threshold is somewhat arbitrary, it follows other recent work prioritizing conservation areas for vultures across the Old World (Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2019a). We then calculated overlap of the core distribution with protected areas (PAs; considering IUCN protected area categories I to VI; IUCN and UNE P-WCMC 2018) and Important Bird and Biodiversity Areas (IBAs; BirdLife International 2020).

Identifying priority conservation areas

We combined the final ensemble species distributions for each vulture species to identify priority conservation areas for all vultures in Ethiopia using software Zonation Version 4.0 (Moilanen et al. Reference Moilanen, Pouzols, Meller, Veach, Arponen, Leppänen and Kujala2014). Zonation ranks areas of conservation importance by iteratively removing grid cells starting from those that have the lowest value for conservation. We used the core area cell removal method (CAZ), such that high priority areas include important habitat for each species, i.e. a ‘leave no species behind’ approach, at a 1-km2 pixel resolution (Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2019a). In order to assess the robustness of the main priority map, we incorporated the prediction uncertainty layers we created (see above) directly into the prioritization analysis using the distribution discounting tool in Zonation (Moilanen et al., Reference Moilanen, Wintle, Elith and Burgman2006), such that the highest-ranking sites have high conservation value and low uncertainty. We ran the zonation analyses three times, setting the uncertainty parameter to subtract 0, 0.5, and 1 SD, from the nominal estimates of the input layers of analysis. We report on the intermediary output (0.5 SD) in the main text and provide the other outputs in Figure S5. Through the prioritization, all species were given the same weight. An alternative approach could be to weight species according to their global conservation status (e.g. Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Di Minin and Moilanen2019). However, we decided to give species equal weight, provided that our exercise was national, and considering that Ethiopia may provide particularly important habitat for some of the less threatened species, particularly Bearded and Egyptian Vultures (Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017). Next, we identified “vulture priority areas” in Ethiopia as the top 30% of the zonation output and calculated the area of overlap and the average priority of each PA and IBA in the country for vulture conservation.

Results

Species distribution models

We produced an ensemble SDM for five of the seven vulture species considered. Ensemble models were based on three data types for Egyptian and Hooded, two data types for Bearded and Rüppell’s, and one data type for White-backed Vulture (Table 1). AUCtest values for models used ranged from 0.606 for the Rüppell’s Vulture road-survey model to 0.928 for the Egyptian Vulture road-survey model (Table 1). Predicted distributions for each species from different data types showed a high level of niche overlap (Hellinger’s I statistic values ranged from 0.759 to 0.973; Table S1), indicating that the different datasets had largely similar predictions of each species’ distribution.

Table 1. Features and evaluation metrics of Maxent distribution models for seven vulture species in Ethiopia based on up to three data types. AUCtest measures the model’s ability to discriminate between conditions at test versus background locations, while ORMTP measures model overfitting (see text for details). For eBird and road-survey models, n refers to the number of observations, whereas for telemetry models n refers to the number of tracked individuals. Models with AUCtest <0.6, shaded in grey, were deemed uninformative and excluded from further analyses. Models were not run for species lacking a sufficient sample size (<20 records for eBird or road-surveys, or <5 individuals for telemetry) and thus their values are NAs, but are included here to show sample size (n).

The leading predictor of Bearded Vulture occurrence was lower temperature, which was equivalent to higher elevations in Ethiopia, as well as more rugged areas. In contrast, Egyptian Vulture occurrence was best predicted by hotter temperatures, i.e. lowland desert areas. Hooded Vulture favoured human-dominated landscapes, while avoiding shrubland. Rüppell’s Vulture distribution was primarily predicted by ruggedness, as well as cooler, higher elevations. White-backed Vulture favoured cropland and forest habitats, proximity to urban centres, and cooler/higher elevations.

The core range of the Egyptian Vulture was best protected (PA coverage = 32.8%, IBA coverage = 7.8%), followed by White-backed (PA = 14.3%, IBA = 5.2%), Rüppell’s (PA = 11.9%, IBA = 9.6%), Bearded (PA = 10.5%, IBA = 13.1%), and Hooded (PA = 9.3%, IBA = 7.1%) Vultures (Table S2). The final ensemble species distribution models for each species are shown in Figure 1 and the average variable importance in the final ensemble model for each species is shown in Figure 2.

Figure 1. Ensemble species distribution models for five vulture species in Ethiopia based on up to three data sources per species (eBird, road-survey, satellite telemetry), weighted by their respective predictive accuracy, and summed. The predicted value, i.e. relative habitat suitability, ranges from 0 to 1 and is a measure of the relative probability that a species occurs on the landscape.

Figure 2. Average percent contribution of each predictor variable to the final ensemble distribution model for each species.

Priority conservation areas

Vulture conservation priority areas occurred throughout the highlands of Ethiopia, as well as the lowlands in the north-east of the country (Figure 3a). Southern and eastern Ethiopia had generally lower priority. The highest priority protected areas included Simien Mountains and Yangudi-Rassa National Parks, as well as several national forests (Figure 3b). These results were robust to uncertainty in predicted vulture distributions (Figure S5). Most of the highest priority IBAs were located in the vicinity of the capital city Addis Ababa in central Ethiopia (Figure 3c). The top 30% vulture priority areas encompassed an area of 338,388 km2, of which 19.4% was covered by protected areas and 8.6% was covered by IBAs. See Appendices S2 and S3 for a full list of protected areas and IBAs with their priority rank.

Figure 3. Spatial conservation priorities for vultures in Ethiopia. Panel a) shows the overall spatial priorities, which is the output of the zonation analysis. Panels b) and c) show the rank priority of the protected area (PAs) and Important Bird and Biodiversity Area (IBAs) networks, as well as the top-ten highest priority site. The colours of the filled polygons correspond with the scale, ranging from 0 to 1, which is a relative measure of spatial conservation importance for vultures.

Discussion

We provide the first quantitative analysis of vulture distributions in Ethiopia to identify spatial conservation priorities. We created ensemble models with up to three data sources per species to incorporate as much information as possible regarding each species’ ecological niche. Priority areas occurred broadly in the central highlands of Ethiopia, which provide particularly important habitat for Bearded, Hooded, Rüppell’s and White-backed Vultures, as well as the lowlands of north-eastern Ethiopia, which are particularly valuable for the Egyptian Vulture (Figure 1, Figure 3). One-fifth of vulture priority areas were covered by protected areas, while less than one-tenth were covered by IBAs. Species core ranges were variably protected, with the Egyptian Vulture faring best at one-third protected, followed by White-backed Vulture at one-sixth, and all other species at one-tenth. That the vast majority (i.e. approximately 80%) of priority areas for vultures in Ethiopia fall outside the protected area network is in contrast to other studies in Africa that have shown vultures persist primarily within protected areas across the continent (Thiollay Reference Thiollay2006, Virani et al. Reference Virani, Kendall, Njoroge and Thomsett2011, Pomeroy et al. Reference Pomeroy, Shaw, Opige, Kaphu, Ogada and Virani2014, Murn et al. Reference Murn, Mundy, Virani, Borello, Holloway and Thiollay2016), with exceptions for Egyptian and Hooded Vultures (Buechley et al. Reference Buechley, McGrady, Çoban and Şekercioğlu2018, Henriques et al. Reference Henriques, Granadeiro, Monteiro, Nuno, Lecoq, Cardoso, Regalla and Catry2018). Our results are strikingly similar to global Old World vulture conservation priority areas, for which one-fifth of the top-priority areas are protected and about one-tenth intersect IBAs (Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2019a).

The limited protection of vulture priority areas indicates that conservation actions must be effective outside the existing protected area network. As has been argued at the global level (Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017, Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2019a), measures of broad spatial and legislative scope are likely to be necessary to conserve vultures and their respective ecosystem services. Such measures could include legislation aimed at controlling the availability and use of veterinary drugs and poisons that are drivers of vulture declines across Africa (Ogada Reference Ogada2014). Further, comprehensive environmental impact assessments should occur prior to energy infrastructure developments to evaluate and mitigate threats to vultures and other biodiversity (Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2019a). It will also be imperative to work closely and collaboratively with stakeholders to identify and address local issues, such as human-wildlife conflict, which could contribute to vulture declines. Given that Ethiopia has a disproportionate role in supporting vulture populations (Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017, Buechley et al. Reference Buechley, Santangeli, Girardello, Neate-Clegg, Oleyar, McClure and Şekercioğlu2019, Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2019a), yet has among the least economic resources globally (World Bank 2019a), we call for increased international attention to and funding for vulture conservation efforts here.

An important next step to improve vulture conservation in the Horn of Africa will be a more detailed assessment of the spatial distribution of threats (Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017). The use of poisons to control feral dog populations by municipalities is apparently widespread and a major danger for vulture populations (Abebe Reference Abebe2013). Further, electrocution on powerlines has been shown to kill Egyptian Vultures in neighboring Sudan (Angelov et al. Reference Angelov, Hashim and Oppel2012), and, similarly concerning levels of electrocution of vultures have been recently reported in Ethiopia (Bakari et al. Reference Bakari, Mengistu, Tesfaye, Ruffo, Oppel, Arkumarev and Nikolov2020). Ethiopia also has some of the largest wind energy facilities in Africa, with ambitious development targets to grow this sector (Asress et al. Reference Asress, Simonovic, Komarov and Stupar2013). While renewable energy expansion is generally a positive strategy for the development of Ethiopia and global environmental sustainability, poorly situated turbines and high voltage transmission lines can be highly detrimental to soaring birds and other wildlife (Barrios and Rodríguez Reference Barrios and Rodríguez2004, Sánchez-Zapata et al. Reference Sánchez-Zapata, Clavero, Carrete, DeVault, Hermoso, Losada, Polo, Sánchez-Navarro, Pérez-García, Botella, Ibáñez and Donázar2016). Climate change is also expected to severely impact the region (Soultan et al. Reference Soultan, Wikelski and Safi2019), and, given that temperature was one of the most important predictors for the occurrence of Bearded, Rüppell’s, and White-backed Vultures (Figure 2), it could be contributing to range contractions (Simmons and Jenkins Reference Simmons and Jenkins2007, Phipps et al. Reference Phipps, Diekmann, MacTavish, Mendelsohn, Naidoo, Wolter and Yarnell2017). Further, we call for a review and update of the IBA network in Ethiopia, which is outdated and largely focused on wetland and riparian areas. The vulture conservation priority areas that were identified herein could be the basis for several additional IBAs, per the A1 Global IBA Criterion for globally threatened species (BirdLife International 2020).

The predicted distributions of vultures herein indicate that Egyptian, Hooded, and White-backed Vultures favour human dominated landscapes and proximity to urban areas. The Egyptian Vulture is known to regularly congregate at refuse dumps (Tauler-Ametller et al. Reference Tauler-Ametller, Hernández-Matías, Pretus and Real2017, Buechley et al. Reference Buechley, McGrady, Çoban and Şekercioğlu2018, McGrady et al. Reference McGrady, Karelus, Rayaleh, Sarrouf Willson, Meyburg, Oli and Bildstein2019), and Hooded Vultures are known to associate with humans, particularly in West Africa and Ethiopia (Mundy et al. Reference Mundy, Bunchart, Ledger and Piper1993, Mullié et al. Reference Mullié, Couzi, Diop, Piot, Peters and Reynaud2017, Henriques et al. Reference Henriques, Granadeiro, Monteiro, Nuno, Lecoq, Cardoso, Regalla and Catry2018, Thompson et al. Reference Thompson, Barber, Bechard, Botha, Wolter, Neser, Buechley, Reading, Garbett, Hancock, Maude, Virani, Thomsett, Lee, Ogada, Barlow and Bildstein2020). In contrast, the White-backed Vulture is more widely regarded as averse to human dominated landscapes, or having already declined precipitously in them, to the point that it occurs mainly within protected areas across much of Africa (e.g. Thiollay Reference Thiollay2006, Virani et al. Reference Virani, Kendall, Njoroge and Thomsett2011, Pomeroy et al. Reference Pomeroy, Shaw, Opige, Kaphu, Ogada and Virani2014; although note that tracked immatures in southern Africa spent the majority of their time outside of protected areas (Phipps et al. Reference Phipps, Willis, Wolter and Naidoo2013). Provided that we tracked Hooded and Egyptian Vultures that were caught in the vicinity of towns, there could be some sampling bias in our models. However, these results are fully congruent with our extensive observations throughout the country. Vulture coexistence with humans in Ethiopia may be regarded as a mutualistic relationship (Gangoso et al. Reference Gangoso, Agudo, Anadon, De la Riva, Suleyman, Porter and Donazar2013, Moleón et al. Reference Moleón, Sánchez-Zapata, Margalida, Carrete, Owen-Smith and Donázar2014). Vulture ecosystem services, while rarely robustly quantified, are now widely recognized (Şekercioğlu Reference Şekercioğlu2006, Buechley and Şekercioğlu Reference Buechley and Şekercioğlu2016a,b, Devault et al. Reference Devault, Beasley, Olson, Moleón, Carrete, Margalida, Sánchez-zapata, Sekercioglu, Wenny and Whelan2016), including efficient nutrient cycling and carrion removal, reduction in greenhouse gas emissions (Morales-Reyes et al. Reference Morales-Reyes, Pérez-García, Moleón, Botella, Carrete, Lazcano, Moreno-Opo, Margalida, Donázar and Sánchez-Zapata2015, Reference Morales-Reyes, Pérez-García, Moleón, Botella, Carrete, Donázar, Cortés-Avizanda, Arrondo, Moreno-Opo, Jiménez, Margalida and Sánchez-Zapata2017), controlling problematic facultative scavengers and insects (Ogada et al. Reference Ogada, Torchin, Kinnaird and Ezenwa2012b, Buechley and Şekercioğlu Reference Buechley and Şekercioğlu2016a), and potentially limiting the spread of disease (Markandya et al. Reference Markandya, Taylor, Longo, Murty, Murty and Dhavala2008, Buechley and Şekercioğlu Reference Buechley and Şekercioğlu2016a, Şekercioğlu et al. Reference Şekercioğlu, Whenny, Whelan, Şekercioğlu, Whenny and Whelan2016, Plaza et al. Reference Plaza, Blanco and Lambertucci2020). Nonetheless, the association of some species of vultures with humans, and their reliance on human waste may constitute an ecological trap for vultures in the 21st century (Ogada et al. Reference Ogada, Keesing and Virani2012a, Buechley and Şekercioğlu Reference Buechley and Şekercioğlu2016a).

Notably, we found that the Egyptian Vulture distribution was concentrated in north-eastern Ethiopia. Other research has shown the importance of this area for the species, particularly for overwintering migrant populations (Arkumarev et al. Reference Arkumarev, Dobrev and Abebe2014, Buechley et al. Reference Buechley, McGrady, Çoban and Şekercioğlu2018). However, surveys and observations of the species in Ethiopia in the 20th century (Ash and Atkins Reference Ash and Atkins2009) indicated a more widespread distribution, with resident populations throughout the country. Our sample of occurrence records for this species was likely biased towards migrants because 13 of 15 tracked Egyptian Vultures were migrants and road surveys were done when most migrants had already arrived in Ethiopia. Nonetheless, that we had few records across much of Ethiopia indicates potential decline in the resident population in Ethiopia, and merits further investigation.

Unfortunately, we were unable to produce informative models of the distributions of Lappet-faced and White-headed Vultures, due to low sample sizes. Surveys in the last century indicated widespread distributions of both of these species in Ethiopia (Ash and Atkins Reference Ash and Atkins2009), and recent work has predicted occurrence of the latter species in Ethiopia (Murn et al. Reference Murn, Mundy, Virani, Borello, Holloway and Thiollay2016). However, these species were the least reported from citizen science observations and road surveys (Table 1). This is not unexpected, as these species are more solitary and averse to human activity. We encourage further targeted survey work on these two species to better understand their distributions, status, and priority conservation areas in the Horn of Africa. Additionally, SDMs built on telemetry and citizen science data from White backed Vulture were uninformative (AUC <0.6; Table 1). We believe this occurred because the White-backed Vulture is a generalist species that can be found across all habitats and elevations, and in both human-dominated and wild environments in Ethiopia, which made it difficult to identify strong signals of habitat association in our models.

There are several potential limitations to our data sources and analyses. In addition to the paucity of data for certain species, as discussed above, we had sparse data from the Somali Region of south-eastern Ethiopia due to security issues inhibiting access. Thus, we caution that the modelled predictions for this region are extrapolatory and should be taken as a first indication of potential areas to target further investigations. Further, we recognize that spatial sampling biases across the three data types could have influenced our models. For example, road surveys were restricted to areas proximate to roads, citizen science data tend to be concentrated in human-populated areas, and tracking data may be influenced by individual behaviour and tagging locations. We worked to address these potential biases, following recommendations in Merow et al. (Reference Merow, Smith and Silander2013), by limiting the background sample in each species-data type model to areas that were surveyed or realistically accessible to tracked individuals. To assess how robust our predictions are, we also directly quantified niche overlap between different species-data-type model predictions (Table S1), quantified uncertainty in our predictions (Figure S4, Figure S5), and deducted areas with more uncertainty in our prioritization exercise. Overall results indicated that our models were largely robust to uncertainty. Further, our models lacked some predictor variables that could be important for vulture space use (for example, the location of carcasses). We expect that predictions could thus be further refined if predictive layers existed that more accurately captured important local aspects of vulture ecology (see relevant discussion here: Efrat et al. Reference Efrat, Hatzofe and Berger-tal2020, Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2020).

In summary, we provide here the first detailed assessment of vulture distributions and spatial conservation priorities for one of the most important countries for vulture conservation worldwide. The low coverage of vulture priority areas by protected areas, as well as the association between vultures and human-dominated landscapes, indicates a need for actions of wide temporal and spatial scope (Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017, Santangeli et al. Reference Santangeli, Girardello, Buechley, Botha, Minin and Moilanen2019a, Perrig et al. Reference Perrig, Lambertucci, Cruz, Alarcón, Plaza, Middleton, Blanco, Sánchez-Zapata, Donázar and Pauli2020). Enactment and enforcement of legislation will likely be essential to reduce the availability and use of veterinary drugs and poisons that threaten wildlife (Margalida et al. Reference Margalida, Bogliani, Bowden, Donazar, Genero, Gilbert, Karesh, Kock, Lubroth, Manteca, Naidoo, Neimanis, Sanchez-Zapata, Taggart, Vaarten, Yon, Kuiken and Green2014, Ogada Reference Ogada2014). Similarly, use of best practices and thorough implementation of environmental impact assessments prior to energy development projects are necessary. Several international conventions and frameworks exist to facilitate development and enforcement of policy relevant to vulture conservation. For example, Ethiopia is a signatory of the Convention on Biological Diversity (CBD), Convention of Migratory Species (CMS), and the Convention on the International Trade of Endangered Species of Wild Fauna and Flora (CITES), which provide guidelines to support the conservation of wildlife, including vultures. Encouragingly, in 2020 Ethiopia also became a signatory of the Raptors MoU, under CMS, which supports specific actions to reduce mortality of raptors and vultures from poisoning and energy infrastructure (Botha et al. Reference Botha, Andevski, Bowden, Gudka, Safford, Tavares and Williams2017). Such actions, if taken while considering the local social context and in a manner that aids sustainable development, would help to conserve threatened vultures and their critical ecological functions, while also benefiting many other species occurring in this biodiversity hotspot region of global importance. Further, we believe that it is imperative to reduce social and economic inequalities, promote peace, and invest in community-led initiatives in order to achieve biodiversity conservation in Ethiopia (Santangeli et al. Reference Santangeli, Girardello, Buechley, Eklund and Phipps2019b), especially provided the recent increase in social unrest and violence in the country.

Supplementary Materials

To view supplementary material for this article, please visit http://doi.org/10.1017/S0959270921000228.

Acknowledgements

We thank the Ethiopia Wildlife Conservation Authority for collaboration and research permits. We thank the Ethiopia Wildlife and Natural History Society and Addis Ababa University Department of Zoology for collaboration and technical support. Thanks to Sisay Seyfu for facilitating all work, and to Andres de la Cruz Muñoz, Juan Ramirez Roman, and Gabriel Caucal for participating in surveys. We thank Steffen Oppel, whose constructive comments helped to improve the manuscript. We are grateful to HawkWatch International, the National Geographic Society, the University of Utah, the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC 2117 – 422037984, the Wallace Research Foundation, The Peregrine Fund, Dutch zoos (notably Avifauna and Blijdorp), and WWF Netherlands for funding. The following HawkWatch International donors sponsored the vulture tracking project: Antczak Polich Law, Buddy Woodhouse, Circle of Life Fund, Doug and Tana Hunter, Eva Carlston Academy, Glen and Anneli Bowen, Jane Tatchell, Julia Shaw, Kirsten Collins, Nancy and John Matro, Scott and Amy Florell, SWCA Environmental Consultants, Team Kaddas, Valerie Walker, and Walter and Karen Loewenstern. AS was supported by a Finnish Academy fellowship (1307909).

References

Abebe, Y. D. (2013) Mass dog poisoning operation in Addis Ababa can have severe repercussions on vulture populations. Vulture News 64: 7476.Google Scholar
Alroy, J. (2017) Effects of habitat disturbance on tropical forest biodiversity. Proc. Natl. Acad. Sci. USA. 114: 60566061.10.1073/pnas.1611855114CrossRefGoogle ScholarPubMed
Angelov, I., Hashim, I. and Oppel, S. (2012) Persistent electrocution mortality of Egyptian Vultures Neophron percnopterus over 28 years in East Africa. Bird Conserv. Internatn. 23: 16.10.1017/S0959270912000123CrossRefGoogle Scholar
Arkumarev, V., Dobrev, V. and Abebe, Y. (2014) Congregations of wintering Egyptian Vultures Neophron percnopterus in Afar, Ethiopia: present status and implications for conservation. Ostrich 85: 139145.10.2989/00306525.2014.971450CrossRefGoogle Scholar
Aryal, A., Shrestha, U. B., Ji, W., Ale, S. B., Shrestha, S., Ingty, T., Maraseni, T., Cockfield, G. and Raubenheimer, D. (2016) Predicting the distributions of predator (snow leopard) and prey (blue sheep) under climate change in the Himalaya. Ecol. Evol. 6: 40654075.CrossRefGoogle ScholarPubMed
Ash, J. and Atkins, J. (2009) Birds of Ethiopia and Eritrea: an atlas of distribution. London, UK: Christopher Helm.Google Scholar
Asress, M. B., Simonovic, A., Komarov, D. and Stupar, S. (2013) Wind energy resource development in Ethiopia as an alternative energy future beyond the dominant hydropower. Renew. Sustain. Energy Rev. 23: 366378.10.1016/j.rser.2013.02.047CrossRefGoogle Scholar
Bakari, S., Mengistu, S., Tesfaye, M., Ruffo, A. D., Oppel, S., Arkumarev, V. and Nikolov, S. C. (2020) Technical report under action A3 of the “Egyptian Vulture New LIFE project”: Bird mortality due to hazardous powerlines in East Oromia and Afar regions, Ethiopia, 2019. Nairobi, Kenya: BirdLife Africa.Google Scholar
Barnosky, A. D., Matzke, N., Tomiya, S., Wogan, G. O. U., Swartz, B., Quental, T. B., Marshall, C., McGuire, J. L., Lindsey, E. L., Maguire, K. C., Mersey, B. and Ferrer, E. A. (2011). Has the Earth’s sixth mass extinction already arrived? Nature 471: 5157.10.1038/nature09678CrossRefGoogle ScholarPubMed
Barrios, L. and Rodríguez, A. (2004) Behavioural and environmental correlates of soaring-bird mortality at on-shore wind turbines. J. Appl. Ecol. 41: 7281.10.1111/j.1365-2664.2004.00876.xCrossRefGoogle Scholar
Benson, J. F. (2013) Improving rigour and efficiency of use-availability habitat selection analyses with systematic estimation of availability. Methods Ecol. Evol. 4: 244251.10.1111/2041-210x.12006CrossRefGoogle Scholar
BirdLife International (2020) BirdLife Data Zone: Global IBA criteria (consulted September 2020: http://datazone.birdlife.org/site/ibacritglob).Google Scholar
Botha, A. J., Andevski, J., Bowden, C. G., Gudka, M., Safford, R. J., Tavares, J. and Williams, N. P. (2017) Multi-species action plan to conserve African-Eurasian vultures (vulture MsAP). Abu Dhabi, United Arab Emirates: Coordinating Unit of the CMS Raptors MOU. (CMS Raptors MOU Technical Publication No. 5. CMS Technical Series No. 35).Google Scholar
Brooks, T. M., Mittermeier, R. A., Fonseca, G. A. B. da, Gerlach, J., Hoffman, M., Lamoreux, J. F., Mittermeier, C. G., Pilgrim, J. D. and Rodrigues, A. S. L. (2006) Global biodiversity conservation priorities. Science 313: 5862.10.1126/science.1127609CrossRefGoogle ScholarPubMed
Buechley, E. R., McGrady, M. J., Çoban, E. and Şekercioğlu, Ç. H. (2018) Satellite tracking of a wide ranging endangered vulture species to target conservation actions in the Middle East and Horn of Africa. Biodivers. Conserv. 27: 22932310.10.1007/s10531-018-1538-6CrossRefGoogle Scholar
Buechley, E. R., Santangeli, A., Girardello, M., Neate-Clegg, M. H., Oleyar, D. McClure, C. J. W. and Şekercioğlu, Ç. H. (2019) Global raptor research and conservation priorities: tropical raptors fall prey to knowledge gaps. Divers. Distrib. 25: 856869.10.1111/ddi.12901CrossRefGoogle Scholar
Buechley, E. R. and Şekercioğlu, Ç. H. (2016a) The avian scavenger crisis: looming extinctions, trophic cascades, and loss of critical ecosystem functions. Biol. Conserv. 198: 220228.10.1016/j.biocon.2016.04.001CrossRefGoogle Scholar
Buechley, E. R. and Şekercioğlu, Ç. H. (2016b) Vultures. Curr. Biol. 26: R560R561.10.1016/j.cub.2016.01.052CrossRefGoogle Scholar
Burnham, K. P. and Anderson, D. R. (2004) Multimodel inference: Understanding AIC and BIC in model selection. Sociol. Methods Res. 33: 261304.10.1177/0049124104268644CrossRefGoogle Scholar
Cadotte, M. W., Carscadden, K. and Mirotchnick, N. (2011) Beyond species: Functional diversity and the maintenance of ecological processes and services. J. Appl. Ecol. 48: 10791087.CrossRefGoogle Scholar
Calenge, C. (2011) Home range estimation in R: the adehabitatHR Package. 1-28. http://cran.r-project.org/web/packages/adehabitatHR/Google Scholar
Coxen, C. L., Frey, J. K., Carleton, S. A. and Collins, D. P. (2017) Species distribution models for a migratory bird based on citizen science and satellite tracking data. Glob. Ecol. Conserv. 11: 298311.CrossRefGoogle Scholar
Devault, T. L., Beasley, J. C., Olson, Z. H., Moleón, M., Carrete, M., Margalida, A. and Sánchez-zapata, J. A. (2016) Ecosystem services provided by avian scavengers. Pp. 235270 in Sekercioglu, C. H., Wenny, D. G., and Whelan, C. J., eds. Why birds matter: Avian ecological function and ecosystem services. Chicago, USA: The University of Chicago Press.Google Scholar
Efrat, R., Hatzofe, O. and Berger-tal, O. (2020) Translating large-scale prioritization models for vultures to local-scale decision making. Conserv. Biol. 5: 13051307.CrossRefGoogle Scholar
Elith, J., Phillips, S. J., Hastie, T., Dudík, M., Chee, Y. E. and Yates, C. J. (2011) A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 17: 4357.CrossRefGoogle Scholar
Evans, L. J., Asner, G. P. and Goossens, B. (2018) Protected area management priorities crucial for the future of Bornean elephants. Biol. Conserv. 221: 365373.CrossRefGoogle Scholar
Ferraz, K., De Barros, M. P. M., De Siqueira, M. F., Alexandrino, E. R., Da Luz, D. T. A. and Do Couto, H. T. Z. (2012) Environmental suitability of a highly fragmented and heterogeneous landscape for forest bird species in south-eastern Brazil. Environ. Conserv. 39: 316324.CrossRefGoogle Scholar
Fick, S. E. and Hijmans, R. J. (2017) Worldclim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37: 43025315.CrossRefGoogle Scholar
Fletcher, R. J. Jr., Hefley, T. J., Robertson, E. P., Zuckerberg, B., McCleery, R. A. and Dorazio, R. M. (2018) A practical guide for combining data to predict species distributions. Ecology 100: e02710.Google Scholar
Gangoso, L., Agudo, R., Anadon, J. D., De la Riva, M., Suleyman, A. S., Porter, R. and Donazar, J. A. (2013) Reinventing mutualism between humans and wild fauna: insights from vultures as ecosystem services providers. Conserv. Lett. 6: 172179.CrossRefGoogle Scholar
Ghosh-Harihar, M., An, R., Athreya, R., Borthakur, U., Chanchani, P., Chetry, D., Datta, A., Harihar, A., Karanth, K. K., Mariyam, D., Mohan, D., Onial, M., Ramakrishnan, U., Robin, V. V., Saxena, A., Shahabuddin, G., Thatte, P., Vijay, V., Wacker, K., Mathur, V. B., Pimm, S. L. and Price, T. D. (2019) Protected areas and biodiversity conservation in India. Biol. Conserv. 237: 114124.CrossRefGoogle Scholar
Guisan, A. and Thuiller, W. (2005) Predicting species distribution: Offering more than simple habitat models. Ecol. Lett. 8: 9931009.CrossRefGoogle ScholarPubMed
Henriques, M., Granadeiro, J. P., Monteiro, H., Nuno, A., Lecoq, M., Cardoso, P., Regalla, A. and Catry, P. (2018) Not in wilderness: African vulture strongholds remain in areas with high human density. PLoS One 13: e0190594.CrossRefGoogle ScholarPubMed
Horns, J. J., Adler, F. R. and Şekercioğlu, Ç. H. (2018) Using opportunistic citizen science data to estimate avian population trends. Biol. Conserv. 221: 151159.10.1016/j.biocon.2018.02.027CrossRefGoogle Scholar
IUCN and UNEP-WCMC (2018) The World Database on Protected Areas (WDPA). Gland, Switzerland: IUCN.Google Scholar
Johnston, A., Hochachka, W. M., Strimas-Mackey, M., Gutierrez, V. R., Robinson, O. J., Miller, E. T., Auer, T., Kelling, S. T. and Fink, D. (2019) Best practices for making reliable inferences from citizen science data: case study using eBird to estimate species distributions. bioRxiv 574392.Google Scholar
Joppa, L. N., O’Connor, B., Visconti, P., Smith, C., Geldmann, J., Hoffmann, M., Watson, J. E. M., Butchart, S. H. M., Virah-Sawmy, M., Halpern, B. S., Ahmed, S. E., Balmford, A., Sutherland, W. J., Harfoot, M., Hilton-Taylor, C., Foden, W., Minin, E. Di, Pagad, S., Genovesi, P., Hutton, J. and Burgess, N. D. (2016. Filling in biodiversity threat gaps. Science 352: 416418.CrossRefGoogle ScholarPubMed
Kendall, C. J., Virani, M. Z., Hopcraft, J. G. C., Bildstein, K. L. and Rubenstein, D. I. (2014) African vultures don’t follow migratory herds: scavenger habitat use is not mediated by prey abundance. PLoS One 9: 18.CrossRefGoogle Scholar
Margalida, A., Bogliani, G., Bowden, C. G. R., Donazar, J. A., Genero, F., Gilbert, M., Karesh, W. B., Kock, R., Lubroth, J., Manteca, X., Naidoo, V., Neimanis, A., Sanchez-Zapata, J. A., Taggart, M. A., Vaarten, J., Yon, L., Kuiken, T. and Green, R. E. (2014) Science and regulation. One Health approach to use of veterinary pharmaceuticals. Science 346: 12961298.CrossRefGoogle ScholarPubMed
Markandya, A., Taylor, T., Longo, A., Murty, M. N., Murty, S. and Dhavala, K. (2008) Counting the cost of vulture decline—an appraisal of the human health and other benefits of vultures in India. Ecol. Econ. 67: 194204.CrossRefGoogle Scholar
McGrady, M. J., Karelus, D. L., Rayaleh, H. A., Sarrouf Willson, M., Meyburg, B.-U., Oli, M. K. and Bildstein, K. (2019) Home ranges and movements of Egyptian Vultures Neophron percnopterus in relation to rubbish dumps in Oman and the Horn of Africa. Bird Study https://doi.org/10.1080/00063657.2018.1561648.CrossRefGoogle Scholar
Merow, C., Smith, M. J. and Silander, J. A. (2013) A practical guide to MaxEnt for modeling species’ distributions: What it does, and why inputs and settings matter. Ecography (Cop.). 36: 10581069.CrossRefGoogle Scholar
Miller, D. A. W., Pacifici, K., Sanderlin, J. S. and Reich, B. J. (2019) The recent past and promising future for data integration methods to estimate species’ distributions. Methods Ecol. Evol. 10: 2237.CrossRefGoogle Scholar
Moilanen, A., Pouzols, F. M., Meller, L., Veach, V., Arponen, A., Leppänen, J. and Kujala, H. (2014) Zonation- Spatial Conservation Planning Methods and Software Version 4 User Manual. https://researchportal.helsinki.fi/en/publications/zonation-spatial-conservation-planning-framework-and-software-verGoogle Scholar
Moilanen, A., Wintle, B. A., Elith, J. and Burgman, M. (2006) Uncertainty analysis for regional-scale reserve selection. Conserv. Biol. 20: 16881697.CrossRefGoogle ScholarPubMed
Moleón, M., Sánchez-Zapata, J. a., Margalida, A., Carrete, M., Owen-Smith, N. and Donázar, J. A. (2014) Humans and scavengers: The evolution of interactions and ecosystem services. Bioscience 64: 394403.CrossRefGoogle Scholar
Morales-Reyes, Z., Pérez-García, J. M., Moleón, M., Botella, F., Carrete, M., Donázar, J. A., Cortés-Avizanda, A., Arrondo, E., Moreno-Opo, R., Jiménez, J., Margalida, A. and Sánchez-Zapata, J. A. (2017) Evaluation of the network of protection areas for the feeding of scavengers in Spain: from biodiversity conservation to greenhouse gas emission savings. J. Appl. Ecol. 54: 11201129.CrossRefGoogle Scholar
Morales-Reyes, Z., Pérez-García, J. M., Moleón, M., Botella, F., Carrete, M., Lazcano, C., Moreno-Opo, R., Margalida, A., Donázar, J. A. and Sánchez-Zapata, J. A. (2015) Supplanting ecosystem services provided by scavengers raises greenhouse gas emissions. Sci. Rep. 5:.CrossRefGoogle ScholarPubMed
Mukherjee, A., Galligan, T. H., Prakash, V., Paudel, K., Khan, U., Prakash, S., Ranade, S., Shastri, K., Dave, R., Donald, P. and Bowden, C. (2014) Vulture Safe Zones to save Gyps Vultures in South Asia. Mistnet 53: 138.Google Scholar
Mullié, W. C., Couzi, F., Diop, M. S., Piot, B., Peters, T. and Reynaud, P. A. (2017) The decline of an urban Hooded Vulture Necrosyrtes monachus population in Dakar, Senegal, over 50 years. Ostrich 88: 131138.CrossRefGoogle Scholar
Mundy, P. P., Bunchart, D., Ledger, J. and Piper, S. (1993) The vultures of Africa. London, UK: Academic Press.Google Scholar
Murn, C., Mundy, P., Virani, M. Z., Borello, W. D., Holloway, G. J. and Thiollay, J.-M. (2016) Using Africa’s protected area network to estimate the global population of a threatened and declining species: a case study of the Critically Endangered White-headed Vulture Trigonoceps occipitalis. Ecol. Evol. 6: 10921103.CrossRefGoogle ScholarPubMed
Muscarella, R., Galante, P. J., Soley-Guardia, M., Boria, R. A., Kass, J. M., Uriarte, M. and Anderson, R. P. (2014) ENMeval: An R package for conducting spatially independent evaluations and estimating optimal model complexity for Maxent ecological niche models. Methods Ecol. Evol. 5: 11981205.CrossRefGoogle Scholar
Neate-Clegg, M.H.C., Horns, J.J., Adler, F.R., Kemahlı Aytekin, M.Ç., Şekercioğlu, Ç.H. 2020. Monitoring the world’s bird populations with community science data. Biological Conservation 248: 108653.CrossRefGoogle Scholar
Ogada, D. L. (2014) The power of poison: pesticide poisoning of Africa’s wildlife. Ann. N.Y. Acad. Sci. 1322: 120.CrossRefGoogle ScholarPubMed
Ogada, D. L., Keesing, F. and Virani, M. (2012a) Dropping dead: causes and consequences of vulture population declines worldwide. Ann. N.Y. Acad. Sci. 1249: 5771.CrossRefGoogle Scholar
Ogada, D., Torchin, M., Kinnaird, M. and Ezenwa, V. (2012b) Effects of vulture declines on facultative scavengers and potential implications for mammalian disease transmission. Conserv. Biol. 26: 453460.CrossRefGoogle Scholar
Perrig, P. L., Lambertucci, S. A., Cruz, J., Alarcón, P. A. E., Plaza, P. I., Middleton, A. D., Blanco, G., Sánchez-Zapata, J. A., Donázar, J. A. and Pauli, J. N. (2020) Identifying conservation priority areas for the Andean condor in southern South America. Biol. Conserv. 243: 108494.CrossRefGoogle Scholar
Phillips, S. J., Anderson, R. P. and Schapire, R. E. (2006) Maximum entropy modeling of species geographic distributions. 190: 231259.CrossRefGoogle Scholar
Phipps, W. L., Diekmann, M., MacTavish, L. M., Mendelsohn, J. M., Naidoo, V., Wolter, K. and Yarnell, R. W. (2017) Due south: A first assessment of the potential impacts of climate change on Cape vulture occurrence. Biol. Conserv. 210: 1625.CrossRefGoogle Scholar
Phipps, W. L., Willis, S. G., Wolter, K. and Naidoo, V. (2013) Foraging ranges of immature African White-Backed Vultures (Gyps africanus) and their use of protected areas in Southern Africa. PLoS One 8: e52813.CrossRefGoogle ScholarPubMed
Plaza, P. I., Blanco, G. and Lambertucci, S. A. (2020) Implications of bacterial, viral and mycotic microorganisms in vultures for wildlife conservation, ecosystem services and public health. Ibis 162: 11091124.CrossRefGoogle Scholar
Plaza, P. I., Martínez-López, E. and Lambertucci, S. A. (2019) The perfect threat: Pesticides and vultures. Sci. Total Environ. 687: 12071218.CrossRefGoogle ScholarPubMed
Pomeroy, D., Shaw, P., Opige, M., Kaphu, G., Ogada, D. L. and Virani, M. Z. (2014) Vulture populations in Uganda: using road survey data to measure both densities and encounter rates within protected and unprotected areas. Bird Conserv. Internatn. 25: 399414.CrossRefGoogle Scholar
R Core Team (2019) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Sytistical Computing.Google Scholar
Randin, C. F., Dirnböck, T., Dullinger, S., Zimmermann, N. E., Zappa, M. and Guisan, A. (2006) Are niche-based species distribution models transferable in space? J. Biogeogr. 33: 16891703.CrossRefGoogle Scholar
Robinson, N., Regetz, J. and Guralnick, R. P. (2014) EarthEnv-DEM90: A nearly-global, void-free, multi-scale smoothed, 90m digital elevation model from fused ASTER and SRTM data. ISPRS J. Photogramm. Remote Sens. 87: 5767.CrossRefGoogle Scholar
Runge, C. A., Martin, T. G., Possingham, H. P., Willis, S. G. and Fuller, R. A. (2014) Conserving mobile species. Front. Ecol. Environ. 12: 395402.CrossRefGoogle Scholar
Runge, C. A., Watson, J. E. M., Butchart, S. H. M., Hanson, J. O., Possingham, H. P. and Fuller, R. A. (2015) Protected areas and global conservation of migratory birds. Science 350: 12551258.CrossRefGoogle ScholarPubMed
Sánchez-Zapata, J. A., Clavero, M., Carrete, M., DeVault, T. L., Hermoso, V., Losada, M. A., Polo, M. J., Sánchez-Navarro, S., Pérez-García, J. M., Botella, F., Ibáñez, C. and Donázar, J. A. (2016) Effects of renewable energy production and infrastructure on wildlife. Wildl. Res. Monogr. Springer International Publishing, Switzerland, pp. 97123.Google Scholar
Santangeli, A., Girardello, M., Buechley, E., Botha, A, Di Minin, E., and Moilanen, A. (2019). Priority areas for conservation of Old World vultures. Conservation Biology 33:10561065.CrossRefGoogle ScholarPubMed
Santangeli, A., Girardello, M., Buechley, E., Botha, A., Minin, E. Di and Moilanen, A. (2019a) Priority areas for conservation of Old World vultures. Conserv. Biol. 5: 10561065.CrossRefGoogle Scholar
Santangeli, A., Girardello, M., Buechley, E., Botha, A., Minin, E. Di and Moilanen, A. (2020) Importance of complementary approaches for efficient vulture conservation: reply to Efrat et al. Conserv. Biol. 10.1111/cobi.13579.10.1111/cobi.13579CrossRefGoogle Scholar
Santangeli, A., Girardello, M., Buechley, E. R., Eklund, J. and Phipps, W. L. (2019b) Navigating spaces for implementing raptor research and conservation under varying levels of violence and governance in the Global South. Biol. Conserv. 239: 108212.CrossRefGoogle Scholar
Şekercioğlu, Ç. H. (2006) Ecological significance of bird populations. Handbook of the birds of the world. 11: 1551. Barcelona, Spain: Lynx Edicions.Google Scholar
Şekercioğlu, Ç. H. (2010) Ecosystem functions and services. Pp.4572 in Sodhi, N. S. and Ehrlich, P., eds. Conservation biology for all. Oxford, UK: Oxford University Press.CrossRefGoogle Scholar
Şekercioğlu, Ç. H., Whenny, D. and Whelan, C. J. (2016) Why birds matter: avian ecological function and ecosystem services, Şekercioğlu, Ç. H. Whenny, D. and Whelan, C. J. (eds.), Chicago, USA: University of Chicago Press.CrossRefGoogle Scholar
Simmons, R. E. and Jenkins, A. R. (2007) Is climate change influencing the decline of Cape and Bearded Vulures in southern Africa ? Vulture News 56: 4151.Google Scholar
Soultan, A., Wikelski, M. and Safi, K. (2019) Risk of biodiversity collapse under climate change in the Afro-Arabian region. Sci. Rep. 9: 112.CrossRefGoogle ScholarPubMed
Phillips, S. J. and Dudík, R. E. S. (2020) Maxent software for modeling species niches and distributions (consulted March 2020: http://biodiversityinformatics.amnh.org/open_source/maxent/).Google Scholar
Sullivan, B. L., Aycrigg, J. L., Barry, J. H., Bonney, R. E., Bruns, N., Cooper, C. B., Damoulas, T., Dhondt, A. A., Dietterich, T., Farnsworth, A., Fink, D., Fitzpatrick, J. W., Fredericks, T., Gerbracht, J., Gomes, C., Hochachka, W. M., Iliff, M. J., Lagoze, C., La Sorte, F. A., Merrifield, M., Morris, W., Phillips, T. B., Reynolds, M., Rodewald, A. D., Rosenberg, K. V., Trautmann, N. M., Wiggins, A., Winkler, D. W., Wong, W. K., Wood, C. L., Yu, J. and Kelling, S. (2014) The eBird enterprise: An integrated approach to development and application of citizen science. Biol. Conserv. 169: 3140.CrossRefGoogle Scholar
Tauler-Ametller, H., Hernández-Matías, A., Pretus, J. L. and Real, J. (2017) Landfills determine the distribution of an expanding breeding population of the endangered Egyptian Vulture Neophron percnopterus. Ibis 159: 757768.CrossRefGoogle Scholar
Taylor, A. T., Hafen, T., Holley, C. T., González, A. and Long, J. M. (2020) Spatial sampling bias and model complexity in stream-based species distribution models: A case study of Paddlefish (Polyodon spathula) in the Arkansas River basin, USA. Ecol. Evol. 10: 705717.CrossRefGoogle ScholarPubMed
Thiollay, J.-M. (2006) The decline of raptors in West Africa: long-term assesment and the role of protected areas. Ibis 148: 240254.CrossRefGoogle Scholar
Thompson, L. J., Barber, D., Bechard, M., Botha, A. J., Wolter, K., Neser, W., Buechley, E. R., Reading, R., Garbett, R., Hancock, P., Maude, G., Virani, M. Z., Thomsett, S., Lee, H., Ogada, D., Barlow, C. R. and Bildstein, K. L. (2020) Variation in monthly sizes of home-ranges of Hooded Vultures Necrosyrtes monachus in western, eastern, and southern Africa. Ibis 162: 13241338.CrossRefGoogle Scholar
Trimble, M. J. and van Aarde, R. J. (2012) Geographical and taxonomic biases in research on biodiversity in human-modified landscapes. Ecosphere 3: art119.CrossRefGoogle Scholar
Venter, O., Sanderson, E. W., Magrach, A., Allan, J. R., Beher, J., Jones, K. R., Possingham, H. P., Laurance, W. F., Wood, P., Fekete, B. M., Levy, M. A. and Watson, J. E. M. (2016) Global terrestrial human footprint maps for 1993 and 2009. Sci. Data 3: 110.CrossRefGoogle ScholarPubMed
Virani, M., Kendall, C., Njoroge, P. and Thomsett, S. (2011) Major declines in the abundance of vultures and other scavenging raptors in and around the Masai Mara ecosystem, Kenya. Biol. Conserv. 144: 746752.CrossRefGoogle Scholar
Warren, D. L., Glor, R. E. and Turelli, M. (2008) Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution. Evolution 62: 28682883.CrossRefGoogle ScholarPubMed
Wei, T., and Simko, V. (2013). corrplot: visualization of a correlation matriz. R package. <cran.us.r-project.org/web/packages/corrplot/corrplot.pdf>..>Google Scholar
World Bank (2019a) GDP per capita (consulted February 2019: https://data.worldbank.org/indicator/ny.gdp.pcap.cd?year_high_desc=false).Google Scholar
World Bank (2019b). Population growth (annual %). https://data.worldbank.org/indicator/SP.POP.GROW?view=chart). (consulted February 2019).Google Scholar
Figure 0

Table 1. Features and evaluation metrics of Maxent distribution models for seven vulture species in Ethiopia based on up to three data types. AUCtest measures the model’s ability to discriminate between conditions at test versus background locations, while ORMTP measures model overfitting (see text for details). For eBird and road-survey models, n refers to the number of observations, whereas for telemetry models n refers to the number of tracked individuals. Models with AUCtest <0.6, shaded in grey, were deemed uninformative and excluded from further analyses. Models were not run for species lacking a sufficient sample size (<20 records for eBird or road-surveys, or <5 individuals for telemetry) and thus their values are NAs, but are included here to show sample size (n).

Figure 1

Figure 1. Ensemble species distribution models for five vulture species in Ethiopia based on up to three data sources per species (eBird, road-survey, satellite telemetry), weighted by their respective predictive accuracy, and summed. The predicted value, i.e. relative habitat suitability, ranges from 0 to 1 and is a measure of the relative probability that a species occurs on the landscape.

Figure 2

Figure 2. Average percent contribution of each predictor variable to the final ensemble distribution model for each species.

Figure 3

Figure 3. Spatial conservation priorities for vultures in Ethiopia. Panel a) shows the overall spatial priorities, which is the output of the zonation analysis. Panels b) and c) show the rank priority of the protected area (PAs) and Important Bird and Biodiversity Area (IBAs) networks, as well as the top-ten highest priority site. The colours of the filled polygons correspond with the scale, ranging from 0 to 1, which is a relative measure of spatial conservation importance for vultures.

Supplementary material: File

Buechley et al. supplementary material

Buechley et al. supplementary material

Download Buechley et al. supplementary material(File)
File 1.5 MB