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 km² (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.

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 km² (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.

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

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.