Abstract
Vultures are among the most vulnerable birds in the world. The bearded vulture (Gypaetus barbatus) is among the threatened species of vultures and listed as Near Threatened. The species is widely distributed across the Palearctic, Afrotropical, and Indomalayan regions. The species faces several threats such as poisoning, direct persecution, habitat degradation, and collisions with powerlines and wind power farms. Thus, knowing the global habitat suitability of the species and environmental predictors of the species distribution can facilitate the species conservation. In this study, we applied a maximum entropy approach, 10,585 distribution records, and 10 environmental variables to model the bearded vulture's global habitat suitability at high spatial resolution [30-arc-second (1 km)]. We also estimated protected area coverage for the species' suitable habitats. We identified 8,117,231 km2 of suitable habitat for the species across its global range in Europe, Asia, and Africa. The results showed that topographic diversity is the most important predictor of the species distribution across its distribution range. Results of estimating the area of suitable habitats of the bearded vulture within protected areas revealed that only 16.26% of the species' suitable habitats are protected. The areas that were identified to have the highest suitability for the species have high priority for the conservation of this iconic species thus these areas should be included in the network of protected areas.
1. Introduction
Human activities over the past decades have led to habitat destruction and a decrease in numerous species (Hanski, 2011; Hansen et al., 2013; Haddad et al., 2015). Climate change, land use change, and road development are among the human activities that are causing the decline in biodiversity worldwide (Thomas et al., 2004; Hansen et al., 2013; Haddad et al., 2015; Newbold, 2018). Among the different groups of species, some are more sensitive to human activities such as raptors (hawks, harriers, kites, eagles, falcons, owls, and vultures). This group provides critical ecosystem services such as disease mitigation, agricultural production, and waste-disposal services (O'Bryan et al., 2018; Carucci et al., 2022). Among raptors, vultures are the most threatened and sensitive group (Botha et al., 2017). In fact, due to human persecution, inadvertent poisoning, and mortality associated with human infrastructures [i.e., roads, buildings, and wind power farms (Ogada et al., 2012; Safford et al., 2019)], vultures have become one of the most threatened taxa on the earth (Botha et al., 2017).
The bearded vulture (Gypaetus barbatus), a cliff-nesting bird of prey, is among the most threatened species of vultures and is listed as Near Threatened by the IUCN (BirdLife International, 2022). The species is widely distributed across the Palearctic, Afrotropical, and Indomalayan regions, with an area of ~61,700,000 km2 (BirdLife International, 2022). There has been a major decline in the distribution range of the habitat specialists' avian scavengers such as bearded vulture, Gypaetus barbatus (Acharya et al., 2010; Safford et al., 2019). As one of the most well-known vulture species, the bearded vulture has experienced a severe population decline in many of its distribution regions and may soon become endangered or reach the status of local extinction (Green et al., 2004; Sheykhi Ilanloo et al., 2020). However, population trends varied significantly throughout the species range; while the European population species has increased since 1980, some other local populations experienced a decline of 89.3% in Nepal (BirdLife International, 2022). At various phases of their life, bearded vulture chooses altitude grasslands and steep pastures as their primary habitat. When looking for food, bearded vultures often soar near mountain slopes or over the ground in their visual search for food, looking for bones and other remnants of carcasses (Margalida, 2008). Only free landscapes allow for this type of search for food (BirdLife International, 2022). Since the bearded vulture represents a significant umbrella species, any measures taken to ensure the survival of this species will also benefit other species that co-exist with it. As one of the most charismatic vulture types (Aguilera-Alcalà et al., 2020), the bearded vulture can only be found in the highest mountain ranges in Asia, Europe, and Africa (Orta et al., 2020). Our knowledge of the bearded vulture distribution varies from one mountain range to another one across the species distribution range (BirdLife International, 2022). While the species distribution is well documented in the Central European Highlands and Southeast European Highlands, it is not well known in other mountain ranges across its global distribution ranges, particularly in Africa (BirdLife International, 2022).
Species distribution models (SDMs) are used for many applications in ecology and conservation (Guisan and Thuiller, 2005; Austin and Van Niel, 2011; Araújo et al., 2019; Zurell et al., 2020). SDMs use species occurrence data along with environmental data characterizing climate, topography, and vegetation and predict spatial distribution/habitat suitability of target species (Guisan et al., 2017). These models are among the most frequently used methods in studying avian ecology, evolution, biogeography, and conservation (Engler et al., 2017; Fourcade et al., 2017; Brambilla et al., 2018, 2019, 2020; Ramellini et al., 2019; Burns et al., 2020; Ferrer-Paris and Sánchez-Mercado, 2021; Li et al., 2021; Lu et al., 2021; Mudereri et al., 2021; Avotins et al., 2022; Condro et al., 2022; Escobar-Luján et al., 2022). They have been applied to identifying suitable habitats of species (TellerÃa et al., 2019; Song et al., 2020; Bai et al., 2022; Buechley et al., 2022), predicting the impacts of climate change on species distribution (Escobar-Luján et al., 2022; Zhu et al., 2022), mapping avian hotspots and richness (Moradi et al., 2019), and effectiveness of protected areas for conservation of species (Buechley et al., 2022; Escobar-Luján et al., 2022; Tian et al., 2022). One particular application of SDMs is to model species' global habitat suitability, which is very critical for large-scale conservation planning of biodiversity (Panthi et al., 2021; Wang et al., 2022). For instance, Panthi et al. (2021) modeled the global habitat suitability of the Egyptian vulture (Neophron percnopterus) and identified environmental determinants of the species distribution. Thus, SDMs can be used in modeling the global habitat suitability of the bearded vulture at high spatial resolution.
This study aimed to model the global geographic distribution of the bearded vulture, to determine the most important environmental factors affecting the global habitat suitability, and to estimate protected area coverage for the suitable habitat of the species. It is known that at large scales, climate serves as the most important factor in shaping species distribution (Pearson and Dawson, 2003). Thus, we are expecting climatic factors to be the most influential predictor of the bearded vulture's global distribution. Human activities during recent decades are a further factor affecting the species distribution ranges (Newbold, 2018; Xu et al., 2019). For instance, environmental changes triggered by humans via their land usage and the resulting habitat degradation caused a severe decline in the distribution ranges of several species. Hence, we are expecting the human footprint index to play a key role in species distribution. The results of this study will facilitate large-scale conservation planning of this iconic and rapidly declining avian scavenger.
2. Materials and methods
2.1. Species global presence data
Global distribution records (Figure 1) of the bearded vulture (Figure 2) were collected from online databases such as the Global Biodiversity Information Facility (GBIF, 2022), eBird (https://ebird.org/home), and VertNet (http://vertnet.org/). In total, 43,670 distribution records were collected, and then duplicates were deleted. After removing duplicates, 17,506 presence points remained. We then assessed the representation of the species distribution records within major mountain ranges across the species' global distribution range (Snethlage et al., 2022).
Figure 1
Figure 2
2.2. Environmental variables
To map the global habitat suitability of the bearded vulture, we considered 14 environmental variables (Table 1) related to climate, topography, vegetation, and human presence. Climatic variables were obtained from the Wordclim at 1 km spatial resolution (Fick and Hijmans, 2017). To consider topography in our model, we estimated topographic heterogeneity based on the Shuttle Radar Topography Mission (SRTM) elevation model (Jarvis et al., 2008). For vegetation, we used the normalized difference vegetation index (NDVI). The human footprint index was used to quantify the anthropogenic impact on the habitat suitability of the species (Venter et al., 2016a,b). This index was created by combining data on the built environments, population density, electric infrastructure, crop lands, pasture lands, roads, railways, and navigable waterways (Venter et al., 2016a,b). We performed a VIF test and only included variables (10 variables) with VIF values less than 10 in the Maxent model (Guisan et al., 2017). Considering the spatial resolution of environmental variables (1 km resolution), we thinned the dataset to ensure that observations were at least 1 km apart. Thus, 10,585 distribution records with at least a 1 km distance were used in habitat suitability modeling.
Table 1
| Variables | Abbreviations |
|---|---|
| Mean diurnal range | Bio2 |
| Isothermality | Bio3 |
| Temperature seasonality | Bio4 |
| The mean temperature of the wettest quarter | Bio8 |
| The mean temperature of the warmest quarter | Bio10 |
| Precipitation of the driest month | Bio14 |
| Precipitation of the coldest quarter | Bio19 |
| Topographic heterogeneity | TH |
| Human footprint | HF |
| Normalized difference vegetation index | NDVI |
Environmental variables that were used in modeling the global habitat suitability of the bearded vulture.
2.3. Global habitat suitability modeling
Maximum entropy modeling (Maxent) is the most popular algorithm among many algorithms for modeling species habitat suitability. Like the other habitat suitability algorithms, Maxent uses presence data from the target species and predictors that present environmental conditions to find suitable habitats for the target species across the study area. In this study, the Kuenm R package (Cobos et al., 2019) was used for modeling the bearded vulture's global habitat suitability. We applied this package to create Maxent candidate models with multiple combinations of regularization multipliers, feature classes, and sets of variables. Then, the best parameters for modeling were selected based on the statistical significance, predictive power, and model complexity (Cobos et al., 2019). The area under the ROC curve (AUC) was used to evaluate the species distribution model's performance (Swets, 1988). To get an alternative estimate of variable importance for the global habitat suitability model of the bearded vulture, we conducted a jackknife test of variable importance (Phillips et al., 2006). The jackknife test excludes each variable, in turn, and creates a model with the remaining variables. Then, a model is created using each variable in isolation. In addition, a model is created using all variables (Phillips et al., 2006).
2.4. Protected area coverage for the bearded vulture
To estimate the protected area coverage of suitable habitats of the bearded vulture across its global range, the species continuous habitat suitability map was converted into a binary suitable–unsuitable map (Guisan et al., 2017). The 10-percentile training presence threshold was used to convert the continuous map into binary (Phillips et al., 2006). Then, the binary habitat suitability model was overlaid on the protected area layer. Protected areas' data were obtained from the protected planet (www.protectedplanet.net) (UNEP-WCMC and IUCN, 2022).
3. Results
3.1. Distribution in the major mountain range
We showed that the breaded vulture occurs in 26 major mountain ranges. Table 2 shows the number of distribution records in each major mountain range across the species distribution range. The highest number of distribution points was recorded in the Central European Highlands.
Table 2
| Regions | Mountain ranges | Number of records |
|---|---|---|
| Central Asia | Himalaya | 1,767 |
| Hindu Kush | 8 | |
| Karakoram | 240 | |
| Pamir mountains | 23 | |
| Tian Shan | 338 | |
| Tibetan Plateau | 498 | |
| East African highlands | Eastern Rift mountains | 31 |
| Ethiopian Highlands | 391 | |
| Horn of Africa Highlands | 1 | |
| East Asia | Central China Mountains | 6 |
| Europe | Central European highlands | 5,505 |
| East European highlands | 4 | |
| Mediterranean islands | 215 | |
| Northwest European highlands | 92 | |
| South European highlands | 5 | |
| Southeast European highlands | 4,877 | |
| North African highlands | Maghreb mountains | 19 |
| South Siberian mountains | Altai-Sayan region | 233 |
| Mongolian highlands | 42 | |
| Southern African ranges | Coastal Belt (South Africa) | 29 |
| Great Escarpment | 499 | |
| Western Asia | Anatolian/Armenian highlands | 580 |
| Arabian Peninsula | 1 | |
| Caucasus mountains | 937 | |
| Iranian plateau | 148 | |
| Levant ranges | 2 |
The number of distribution records in major mountain ranges in Europe, Asia, and Africa.
3.2. Global habitat suitability model
According to the AUS value (AUC = 0.924), the global habitat suitability model developed for the breaded vulture (Gypaetus barbatus) shows high performance. The Maxent model shows areas with high suitability for the species across Europe, Asia, and Africa. The model is created based on distribution records of the species collected from 2000 to 2022 across the species distribution range. We estimated 8,117,231 km2 of suitable habitat for the species across its global range in Europe, Asia, and Africa (Figure 3).
Figure 3
3.3. Variable importance
Results of estimating the importance of the environmental variables to the bearded vulture global habitat suitability model showed that topographic heterogeneity was the most important variable (with 46.2% contribution) in predicting suitable habitats of the species with a positive relationship (Table 3). Areas with high topographic heterogeneity have higher suitability for the species. Precipitation of the driest month (Bio14) was the second most important variable (with 21.6% contribution) in shaping the species' habitat suitability with a positive influence meaning that areas that receive more precipitation during the driest month have high suitability for the species.
Table 3
| Variable | Percent contribution |
|---|---|
| Topographic heterogeneity | 46.2 |
| Bio14 | 21.6 |
| Bio8 | 12.4 |
| Bio10 | 4.7 |
| Human footprint | 3.8 |
| Bio3 | 3.4 |
| Bio4 | 2.5 |
| Bio19 | 2.2 |
| Bio2 | 1.9 |
| NDVI | 1.3 |
Relative contributions of the environmental variables to the bearded vulture (Gypaetus barbatus) Maxent model.
Figure 4 shows the results of the jackknife test of variable importance. The environmental variable with the highest gain, when used in isolation, is topographic heterogeneity, which therefore appears to have the most useful information by itself. The environmental variable that decreases the gain the most when it is omitted is Bio14, which therefore appears to have the most useful information that is not present in the other variables (Figure 3).
Figure 4
3.4. Representation of suitable habitats within protected areas
Results of estimating the area of suitable habitats of the bearded vulture within protected areas showed that 1,320,533 km2 of suitable habitats of the species are protected. Considering that area of suitable habitat for the species is 8,117,231 km2, thus 16.26% of the suitable habitat of the bearded vulture is currently protected (Figure 5).
Figure 5
4. Discussion
Mountain ecosystems are hotspots of biodiversity and have high priority for conservation (Antonelli et al., 2018; Rahbek et al., 2019). These ecosystems are home to many iconic species in need of urgent conservation actions like the bearded vulture. This study presents a high-resolution global scale habitat suitability map for the bearded vulture to help large-scale conservation planning of the species. The overall extent of the occurrence of the species is estimated to be 61,700,000 km2 (BirdLife International, 2022). However, our model on the global distribution of the species showed that the area of suitable habitat for the species is much lower than what was previously thought (BirdLife International, 2022). This highlights the importance and usefulness of SDMs for precise estimation of the area of suitable habitat for species with conservation importance.
Despite we were expecting that climate will outperform other variables in predicting the species' global habitat suitability, the results showed that topography was the most influential factor. Topographic heterogeneity is an indicator of habitat diversity and species diversity (Stein et al., 2014; Badgley et al., 2017), which, in turn, influences food resources availability to bearded vultures (Sheykhi Ilanloo et al., 2020). Vignali et al. (2021) built several habitat suitability maps with respect to bearded vulture age class (immature and adult) and season (warm and cold) in the Swiss Alpine range. In most of the models, the environmental predictor that contributed most to explaining the habitat suitability of the bearded vulture was food availability. The precipitation of the driest month was the second most important predictor of the species distribution range showing the importance of precipitation for the species' habitat suitability. Sheykhi Ilanloo et al. (2020) modeled the habitat suitability of the bearded vulture in the Kopet Dagh Mountains in Asia and found that the precipitation was the most influential predictor of the species' suitable habitat in the area. Surprisingly, results of variable importance revealed that the human footprint index did not play an important role in shaping the species distribution showing that at least at a global scale, the species' suitable habitat is not affected by human activities. Confirming previous findings, our results show that species responses to environmental factors such as climate, topography, vegetation, and anthropogenic factors are species-specific and depend on the scale of studies (Roland and Taylor, 1997; Holland et al., 2004; McGarigal et al., 2016). We showed that while topography shaped the species' global distribution range, precipitation was the most important determinant of the species distribution in the Kopet Dagh Mountains in Asia (Sheykhi Ilanloo et al., 2020), and food availability was the most important variable in explaining habitat suitability of bearded vulture in the Swiss Alpine range. In comparison with other vulture species, Panthi et al. (2021) modeled the global habitat suitability of the Egyptian vulture (Neophron percnopterus) and found that climate is the most important determinant of the species' global habitat suitability.
Protected areas are important tools for the conservation of biodiversity (Gaston et al., 2008; Mawdsley et al., 2009; Watson et al., 2014; Gray et al., 2016). Threats that species generally face outside protected areas are absent or very low in protected areas (Leberger et al., 2020). Thus, knowing to what extent species' suitable habitats are located within protected areas can help species conservation and future selection of protected areas. In this regard, SDMs are very useful tools to investigate the effectiveness of protected areas in the conservation of species and identify gaps in protected area coverage for species' suitable habitats (Mudereri et al., 2021; Escobar-Luján et al., 2022). Here, we applied SDMs to estimate protected area coverage for the suitable habitat of the species and found that < 17% of the bearded vulture's suitable habitat is currently protected. This highlighted the necessity of new protected area establishment in mountain ecosystems in general and suitable habitats of this iconic and threatened avian species in particular.
Vultures are among the most threatened birds that are in need of urgent conservation actions (Botha et al., 2017). Like many other vultures, the bearded vulture has been listed as Near Threatened because it has undergone a moderately rapid population decline over the past three generations (BirdLife International, 2022). We found a high number of distribution records for the species in the following mountain ranges: Central European Highlands, Southeast European Highlands, Himalayas, Caucasus Mountains, and Anatolian/Armenian Highlands showing that more efforts have been granted for monitoring of this species or avian species in general in these mountains. While few observations have been reported in other mountain ranges such as Iranian Plateau, that the species is resident and common there (Sheykhi Ilanloo et al., 2020). Thus, more effort should be granted in data poor mountain ranges to monitor the species populations, particularly by using GPS transmitters and record their presence, which is a critical step in understanding the species' global distribution and conservation (Di Vittorio et al., 2018). Moreover, as suggested by Di Vittorio et al. (2018), the development of citizen science programs can significantly increase our knowledge about the species distribution and population in data poor regions. We particularly encourage awareness-raising programs to increase local people's awareness about the species and the species' highly suitable habitats identified in this study. The species population in Europe has increased due to conservation actions, thus aforementioned programs and conservation actions are necessary for the highly suitable habitat of the species in Asia and Africa in which the species experiencing dramatic population decline (BirdLife International, 2022). Species distribution models have been frequently used in studying avian ecology and conservation (Engler et al., 2017; Fourcade et al., 2017; Brambilla et al., 2018, 2019, 2020; Ramellini et al., 2019; Burns et al., 2020; Ferrer-Paris and Sánchez-Mercado, 2021; Li et al., 2021; Lu et al., 2021; Mudereri et al., 2021; Avotins et al., 2022; Condro et al., 2022; Escobar-Luján et al., 2022). However, studies that applied SDMs in studying avian global distribution are still limited (Engler et al., 2017; Ferrer-Paris and Sánchez-Mercado, 2021; Stiels et al., 2021). In this study, we modeled the global habitat suitability of the bearded vulture and showed that the species distribution range is very smaller than what was previously thought. In addition, we showed that a small proportion of the species' suitable habitats is located within protected areas. We believe that the establishment of new wind power farms should be avoided in highly suitable habitats or should be established with caution to reduce the species' collision with turbines (Carrete et al., 2012; Murgatroyd et al., 2021). Areas that are identified as the most suitable habitats for the species should be prioritized for the conservation of this charismatic vulture species.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
MY conceived the ideas, designed the methodology, collected the data, and wrote the original draft with inputs from SM (Section 1). AK conducted the analyses. MY, SM, and AK reviewed, edited the text, and gave their final approval for submission. All authors contributed to the article and approved the submitted version.
Acknowledgments
We thank Ali Khani for providing the photo of the bearded vulture. We would also like to thank two reviewers and editors of the Frontiers in Ecology and Evolution for their useful comments and suggestions.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AcharyaR.CuthbertR.BaralH. S.ChaudharyA. (2010). Rapid decline of the bearded vulture gypaetus barbatus in upper Mustang, Nepal. Forktail26, 117–120.
2
Aguilera-AlcalÃN.Morales-ReyesZ.MartÃn-LópezB.MoleónM.Sánchez-ZapataJ. A. (2020). Role of scavengers in providing non-material contributions to people. Ecol. Indic.117, 106643. 10.1016/j.ecolind.2020.106643
3
AntonelliA. W. DKisslingS. G. A.FlantuaBermúdez, A.MulchA. N. (2018). Geological and climatic influences on mountain biodiversity. Nat. Geosci. 11, 718–725. 10.1038/s41561-018-0236-z
4
AraújoM. B.AndersonR. P.BarbosaA. M.BealeC. M.DormannA. F.EarlyR.et al. (2019). Standards for distribution models in biodiversity assessments. Sci. Adv.5, 1–11. 10.1126/sciadv.aat4858
5
AustinM. P.Van NielK. P. (2011). Improving species distribution models for climate change studies: variable selection and scale. J. Biogeogr.38, 1–8. 10.1111/j.1365-2699.2010.02416.x
6
AvotinsA.KerusV. AAunins (2022). National scale habitat suitability analysis to evaluate and improve conservation areas for a mature forest specialist species. Glob. Ecol. Conserv.38, e02218. 10.1016/j.gecco.2022.e02218
7
BadgleyC.SmileyT. M.TerryR.DavisE. B.DeSantisL. R. G.FoxD. L.et al. (2017). Biodiversity and topographic complexity: modern and geohistorical perspectives. Trends Ecol. Evolut.32, 211–226. 10.1016/j.tree.2016.12.010
8
BaiJ.HouP.JinD.ZhaiJ.MaY.ZhaoJ. (2022). Habitat suitability assessment of black-necked crane (Grus nigricollis) in the Zoige Grassland Wetland ecological function zone on the Eastern tibetan plateau. Diversity14, 579. 10.3390/d14070579
9
BirdLife International (2022). Species factsheet: Gypaetus barbatus. Available online at: http://www.birdlife.org (accessed October 20, 2022).
10
BothaA. J.AndevskiJ.BowdenC. G. R. (2017). Multi-species action plan to conserve African-Eurasian vultures. CMS Raptors MOU Technical Publication No. 5. CMS technical series no. 3. Coordinating Unit of the CMS Raptors MOU, United Arab Emirates 2017.
11
BrambillaM.GustinM.CentoM.IlahianeL.CeladaC. (2019). Predicted effects of climate factors on mountain species are not uniform over different spatial scales. J. Avian Biol.50, 2162. 10.1111/jav.02162
12
BrambillaM.Resano-MayorJ.ArlettazR.BettegaC.BinggeliA.BoglianiG.et al. (2020). Potential distribution of a climate sensitive species, the White-winged Snowfinch Montifringilla nivalis in Europe. Bird Conserv. Int.30, 522–532. 10.1017/S0959270920000027
13
BrambillaM.Resano-MayorJ.ScridelD.AnderleM.BoglianiG.BraunischV.et al. (2018). Past and future impact of climate change on foraging habitat suitability in a high-alpine bird species: Management options to buffer against global warming effects. Biol. Conserv.221, 209–218. 10.1016/j.biocon.2018.03.008
14
BuechleyE. R.GirardelloM.SantangeliA.RuffoA. D.AyalewG.AbebeY. D.et al. (2022). Priority areas for vulture conservation in the Horn of Africa largely fall outside the protected area network. Bird Conserv. Int.32, 188–205. 10.1017/S0959270921000228
15
BurnsP.ClarkM.SalasL.HancockS.LelandD.JantzP.et al. (2020). Incorporating canopy structure from simulated GEDI lidar into bird species distribution models. Environ. Res. Lett.15, 095002. 10.1088/1748-9326/ab80ee
16
CarreteM.Sánchez-ZapataJ. A.BenÃtezJ. R.LobónM.MontoyaF.DonázarJ. A.et al. (2012). Mortality at wind-farms is positively related to large-scale distribution and aggregation in griffon vultures. Biol. Conserv.145, 102–108. 10.1016/j.biocon.2011.10.017
17
CarucciT.Whitehouse-TeddK.YarnellR. W.CollinsA.FitzpatrickF.BothaA.et al. (2022). Ecosystem services and disservices associated with vultures: a systematic review and evidence assessment. Ecosys. Serv.56, 101447. 10.1016/j.ecoser.2022.101447
18
CobosM. E.PetersonA. T.BarveN.Osorio-OlveraL. (2019). kuenm: an R package for detailed development of ecological niche models using Maxent. Peer. J.7, e6281. 10.7717/peerj.6281
19
CondroA. A.SyartiniliaH.HiguchiY. A.MulyaniR.RaffiudinL.RusniarsyahY.et al. (2022). Climate change leads to range contraction for Japanese population of the Oriental Honey-Buzzards: implications for future conservation strategies. Global Ecol. Conserv.34, e02044. 10.1016/j.gecco.2022.e02044
20
Di Vittorio HemaM.DendiE. M.AkaniD.CortoneG. C.Lopez-LopezG.et al. (2018). The conservation status of West African vultures: an updated review and a strategy for conservation. Vie Milieu.68, 33–43.
21
EnglerJ. O.StielsD.SchidelkoK.StrubbeD.QuillfeldtP.BrambillaM. (2017). Avian SDMs: current state, challenges, and opportunities. J. Avian. Biol.48, 1483–1504. 10.1111/jav.01248
22
Escobar-LujánJ.Castaño-QuinteroS. M.VillalobosF.Lira-NoriegaA.Chiappa-CarraraX.Yañez-Arenas (2022). Current and future geographic patterns of bird diversity dimensions of the Yucatan Peninsula and their representativeness in natural protected areas. Neotrop. Biodiv.8, 242–252. 10.1080/23766808.2022.2087282
23
Ferrer-ParisJ. R.Sánchez-MercadoA. (2021). Contributions of Distribution modelling to the ecological study of psittaciformes. Diversity:13, 611. 10.3390/d13120611
24
FickS. E.HijmansR. J. (2017). Worldclim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol.37, 4302–4315. 10.1002/joc.5086
25
FourcadeY.BesnardA. G.SecondiJ. (2017). Evaluating interspecific niche overlaps in environmental and geographic spaces to assess the value of umbrella species. J. Avian Biol.48, 1563–1574. 10.1111/jav.01153
26
GastonK. J.JacksonS. F.Cant'u-SalazarL.Cruz-PiñónG. (2008). The ecological performance of protected areas. Annu. Rev. Ecol. Evol. Syst. 3993–113. 10.1146/annurev.ecolsys.39.110707.173529
27
GBIF.org (2022). GBIF Occurrence. 10.15468/dl.zt3zp6. (accessed March 24, 2022).
28
GrayC. L.HillS. L. L.NewboldT.HudsonL. N.BörgerL.ContuS.et al. (2016). Local biodiversity is higher inside than outside terrestrial protected areas worldwide. Nat. Commun7, 12306. 10.1038/ncomms12306
29
GreenR. E.NewtonI.SchultzS. (2004). Diclofenac poisoning as a cause of vulture population declines across the Indian subcontinent. J. Appl. Ecol.41, 793–800. 10.1111/j.0021-8901.2004.00954.x
30
GuisanA.ThuillerW. (2005). Predicting species distribution: offering more than simple habitat models. Ecol. Lett.8, 993–1009. 10.1111/j.1461-0248.2005.00792.x
31
GuisanA.ThuillerW.ZimmermannN. E. (2017). Habitat Suitability and Distribution Models: With Applications in R. Cambridge: Cambridge University Press.
32
HaddadN. M.BrudvigL. A.ClobertJ.DaviesK. F.GonzalezA.HoltR. D.et al. (2015). Habitat fragmentation and its lasting impact on Earth's ecosystems. Sci. Adv.1, e1500052. 10.1126/sciadv.1500052
33
HansenM. C.PotapovP. V.MooreR.HancherM.TurubanovaS. A.TyukavinaA.et al. (2013). High-resolution global maps of 21st-century forest cover change. Science342, 850–853. 10.1126/science.1244693
34
HanskiI. (2011). Habitat loss, the dynamics of biodiversity, and a perspective on conservation. Ambio40, 248–255. 10.1007/s13280-011-0147-3
35
HollandJ. D.BertD. G.FahrigL. (2004). Determining the spatial scale of species' response to habitat. Bioscience54, 227–233.
36
JarvisA.ReuterH. I.NelsonA.GuevaraE. (2008). Hole-Filled SRTM for the Globe Version 4. Available online at: the CGIAR-CSI SRTM 90 m Database. Available online at: https://srtm.csi.cgiar.org (accessed April 15, 2015).
37
LebergerR.RosaI. M. D.GuerraC. A.WolfF.PereiraH. M. (2020). Global patterns of forest loss across IUCN categories of protected areas. Biol. Conserv. 241 108299. 10.1016/j.biocon.2019.108299
38
LiJ.ZhangY.ZhaoL.DengW.QianF.MaK.et al. (2021). Scale and landscape features matter for understanding waterbird habitat selection. Remote Sens.13, 4397. 10.3390/rs13214397
39
LuZ.ZhaiY.MengD.KouG.LiH.LiuC. (2021). Predicting the potential distribution of wintering Asian Great Bustard (Otis tarda dybowskii) in China: Conservation implications. Global Ecol. Conserv.31, e01817. 10.1016/j.gecco.2021.e01817
40
MargalidaA. (2008). Presence of bone remains in the ossuaries of Bearded Vultures (Gypaetus barbatus): storage or nutritive rejection?Auk125, 560–564. 10.1525/auk.2008.07124
41
MawdsleyJ. R. R.O'malleyD. S.Ojima (2009). A review of climate-change adaptation strategies for wildlife management and biodiversity conservation. Conserv. Biol. 23, 1080–1089. 10.1111/j.1523-1739.2009.01264.x
42
McGarigalK.WanH. Y.ZellerK. A.TimmB. C.CushmanS. A. (2016). Multi-scale habitat selection modeling: a review and outlook. Landscape Ecol.31, 1161–1175. 10.1007/s10980-016-0374-x
43
MoradiS.Sheykhi IlanlooS.KafashA.YousefiM. (2019). Identifying high-priority conservation areas for avian biodiversity using species distribution modeling. Ecol. Ind.97, 159–164. 10.1016/j.ecolind.2018.10.003
44
MudereriB. T.ChitataT.ChemuraA.MakaureJ.MukangaC.Abdel-RahmanE. C. (2021). Is the protected area coverage still relevant in protecting the Southern Ground-hornbill (Bucorvus leadbeateri) biological niche in Zimbabwe? Perspectives from ecological predictions. GIScie. Remote Sens.58, 405–424. 10.1080/15481603.2021.1883947
45
MurgatroydM.BoutenW.AmarA. A. (2021). predictive model for improving placement of wind turbines to minimise collision risk potential for a large soaring raptor. J. Appl. Ecol.58, 857–868. 10.1111/1365-2664.13799
46
NewboldT. (2018). Future effects of climate and land-use change on terrestrial vertebrate community diversity under different scenarios. Proceed. Royal Society B: Biol. Sci.285, 20180792. 10.1098/rspb.2018.0792
47
O'BryanC. J.BraczkowskiA. R.BeyerH. L.CarterN. H.WatsonJ. E. M.McDonald-MaddenE. (2018). The contribution of predators and scavengers to human wellbeing. Nat. Ecol. Evol2, 229–236. 10.1038/s41559-017-0421-2
48
OgadaD. L.KeesingF.ViraniM. Z. (2012). Dropping dead: causes and consequences of vulture population declines worldwide. Ann. N. Y. Acad. Sci.1249, 57–71. 10.1111/j.1749-6632.2011.06293.x
49
OrtaJ.de JuanaE.MarksJ. S.SharpeC. J.GarciaE. F. J. (2020). Bearded Vulture (Gypaetus barbatus). in Handbook of the Birds of the World Alive, eds J. del Hoyo, A. Elliott, J. Sargatal, D. A. Christie, and E. de Juana (Barcelona: Lynx Edicions) (accessed November 20, 2020).
50
PanthiS.PariyarS.LowM. (2021). Factors influencing the global distribution of the endangered Egyptian vulture. Sci. Rep.11, 21901. 10.1038/s41598-021-01504-y
51
PearsonR. G.DawsonT. P. (2003). Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful?. Glob. Ecol. Biogeogr12, 361–371. 10.1046/j.1466-822X.2003.00042.x
52
PhillipsS. J.AndersonR. P.SchapireR. E. (2006). Maximum entropy modeling of species geographic distributions. Ecol. Modell.190, 231–259. 10.1016/j.ecolmodel.2005.03.026
53
RahbekC.BorregaardM. K.ColwellR. K.DalsgaardB.HoltB. G.Morueta-HolmeN.et al. (2019). Humboldt's enigma: what causes global patterns of mountain biodiversity?Science365, 1108–1113. 10.1126/science.aax0149
54
RamelliniS.SimonciniA.FicetolaG. F.FalaschiM. (2019). Modelling the potential spread of the Red-billed Leiothrix lutea in Italy. Bird Study66, 550–560. 10.1080/00063657.2020.1732864
55
RolandJ.TaylorP. D. (1997). Insect parasitoid species respond to forest structure at different spatial scales. Nature386, 710–713. 10.1038/386710a0
56
SaffordR.AndevskiJ.BothaA.BowdenC. G. R.CrockfordN.GarbettR.et al. (2019). Vulture conservation: the case for urgent action. Bird Conserv. Int.29, 1–9. 10.1017/S0959270919000042
57
Sheykhi IlanlooS.KhaniA.KafashA.ValizadeganN.AshrafiS.LoercherF.et al. (2020). Applying opportunistic observations to model current and future suitability of the Kopet Dagh Mountains for a Near Threatened avian scavenger. Avian Biol. Res.14, 18–26. 10.1177/1758155920962750
58
SnethlageM. A.GeschkeJ.RanipetaA.JetzW.YoccozN. G.KörnerC.et al. (2022). A hierarchical inventory of the world's mountains for global comparative mountain science. Sci. Data9, 149. 10.1038/s41597-022-01256-y
59
SongK.MiC. R.YangN.SunL.SunY. H.XuJ. L. (2020). Improve the roles of nature reserves in conservation of endangered pheasant in a highly urbanized region. Sci. Reports10, 17673. 10.1038/s41598-020-74724-3
60
SteinA.GerstnerK.KreftH. (2014). Environmental heterogeneity as a universal driver of species richness across taxa, biomes and spatial scales. Ecol. Let.17, 866–880. 10.1111/ele.12277
61
StielsD.BastianH. V.BastianA.SchidelkoK.EnglerJ. O. (2021). An iconic messenger of climate change? Predicting the range dynamics of the European Bee-eater (Merops apiaster). J. Ornithol.162, 631–644. 10.1007/s10336-021-01867-z
62
SwetsJ. A. (1988). Measuring the accuracy of diagnostic systems. Science240, 1285–1293. 10.1126/science.3287615
63
TellerÃaJ. L.FandosG.TenaE.CarbonellR.OnrubiaA.QninbaA.et al. (2019). Constraints on raptor distribution at the southwestern boundary of the Palaearctic: implications for conservation. Biodiv. Conserv.28, 603–619. 10.1007/s10531-018-1677-9
64
ThomasC. D. ACameronR. E.GreenM.BakkenesL. J.BeaumontY. C.et al. (2004). Extinction risk from climate change. Nature427, 145–148. 10.1038/nature02121
65
TianS.LuS.HuaJ.ChangJ.LiJ.ZhangZ.et al. (2022). Integrating habitat suitability modelling and assessment of the conservation gaps of nature reserves for the threatened Reeves's Pheasant. Bird Conserv. Int.32, 384–397. 10.1017/S095927092100023X
66
UNEP-WCMC IUCN (2022) Protected Planet: The World Database on Protected Areas (WDPA) and World Database on Other Effective Area-based Conservation Measures (WD-OECM), July 2022, Cambridge, UK: UNEP-WCMC and IUCN. Available online at: www.protectedplanet.net (accessed July, 2022).
67
VenterO.SandersonE. W.MagrachA.AllanJ. R.BeherJ.JonesK. R.et al. (2016a). Sixteen years of change in the global terrestrial human footprint and implications for biodiversity conservation. Nat. Commun.7, 12558. 10.1038/ncomms12558
68
VenterO.SandersonE. W.MagrachA.AllanJ. R.BeherJ.JonesK. R.et al. (2016b). Global terrestrial Human Footprint maps for 1993 and 2009. Sci. Data3, 160067. 10.1038/sdata.2016.67
69
VignaliS.LörcherF.HegglinD.ArlettazR.BraunischV. (2021). Modelling the habitat selection of the bearded vulture to predict areas of potential conflict with wind energy development in the Swiss Alps. Glob. Ecol. Conserv. 25, e01405. 10.1016/j.gecco.2020.e01405
70
WangD.ShiC.AlamgirK.KwonS.PanL.ZhuY.et al. (2022). Global assessment of the distribution and conservation status of a key medicinal plant (Artemisia annua L.): the roles of climate and anthropogenic activities. Sci. Total Environ.821, 153378. 10.1016/j.scitotenv.2022.153378
71
WatsonJ. E. M.DudleyN.SeganD. B.HockingsM. (2014). The performance and potential of protected areas. Nature515, 67–73. 10.1038/nature13947
72
XuW. B.SvenningJ. C.ChenG. K.ZhangM. G.HuangJ. H.ChengB.et al. (2019). Human activities have opposing effects on distributions of narrow-ranged and widespread plant species in China. Proc. Natl. Acad. Sci. U.S.A.116, 26674–26681. 10.1073/pnas.1911851116
73
ZhuB. R.VerhoevenM. A.VelascoN.Sanchez-AguilarL.ZhangZ.PiersmaT.et al. (2022). Current breeding distributions and predicted range shifts under climate change in two subspecies of Black-tailed Godwits in Asia. Glob. Change Biol.28, 5416–5426. 10.1111/gcb.16308
74
ZurellD.FranklinJ.KönigC.BouchetP. J.DormannC. F.ElithJ.et al. (2020). A standard protocol for reporting species distribution models. Ecography43, 1261–1277. 10.1111/ecog.04960
Summary
Keywords
Gypaetus barbatus, avian conservation, avian ecology, conservation, global habitat suitability, Maxent
Citation
Yousefi M, Mohammadi S and Kafash A (2023) Modeling global habitat suitability and environmental predictor of distribution of a Near Threatened avian scavenger at a high spatial resolution. Front. Ecol. Evol. 11:1112962. doi: 10.3389/fevo.2023.1112962
Received
30 November 2022
Accepted
13 February 2023
Published
16 March 2023
Volume
11 - 2023
Edited by
Peter Convey, British Antarctic Survey (BAS), United Kingdom
Reviewed by
Daniel de Paiva Silva, Goiano Federal Institute (IFGOIANO), Brazil; Corrado Battisti, Roma Tre University, Italy
Updates
Copyright
© 2023 Yousefi, Mohammadi and Kafash.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Saeed Mohammadi smohammadi@uoz.ac.ir
This article was submitted to Biogeography and Macroecology, a section of the journal Frontiers in Ecology and Evolution
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.