Next Article in Journal
Metabolic Responses and Resilience to Environmental Challenges in the Sedentary Batrachoid Halobatrachus didactylus (Bloch & Schneider, 1801)
Previous Article in Journal
Pet Ownership in Aotearoa New Zealand: A National Survey of Cat and Dog Owner Practices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Factors Influencing the Potential Distribution of Globally Endangered Egyptian Vulture Nesting Habitat in Nepal

1
Central Department of Zoology, Institute of Science and Technology, Tribhuvan University, Kirtipur, Kathmandu 44618, Nepal
2
Bird Conservation Nepal, Kathmandu 44600, Nepal
3
Department of Zoology, Prithvi Narayan Campus, Tribhuvan University, Pokhara 33700, Nepal
4
Pokhara Bird Society, Pokhara, Kaski 33700, Nepal
5
Department of National Parks and Wildlife Conservation, Babarmahal, Kathmandu 44600, Nepal
6
Department of Fisheries and Wildlife, Michigan State University, East Lansing, MI 44824, USA
*
Author to whom correspondence should be addressed.
Animals 2023, 13(4), 633; https://doi.org/10.3390/ani13040633
Submission received: 28 November 2022 / Revised: 1 February 2023 / Accepted: 3 February 2023 / Published: 11 February 2023
(This article belongs to the Section Birds)

Abstract

:

Simple Summary

Potentially suitable habitat for nest-site selection of the globally endangered Egyptian vulture in Nepal was influenced primarily by precipitation and forested areas containing cliffs. Sites near to forests and human settlements were most suited habitat for Egyptian vultures in Nepal. Land-use change could negatively impact the nesting and distribution of Egyptian vultures and threaten their persistence in some areas. Potentially suitable nesting areas generally occurred outside of protected areas Therefore, to protect the potential nesting habitat of Egyptian vultures, forested areas with large numbers of cliffs should be conserved.

Abstract

Habitat suitability is crucial to ensure the long-term persistence of species and can be identified based on relationships between species occurrences and underlying abiotic and biotic factors. We identified potential nesting habitat for the Egyptian vulture (Neophron percnopterus) in Nepal using ecological niche modeling with climatic variables. We estimated the currently suitable nesting habitat for Egyptian vulture in Nepal at 38,204 km2. We found a high probability of suitable nesting habitat on east-facing aspects, and the probability of a suitable nesting habitat was greater in more mountainous areas, particularly in central and western regions of Nepal. Precipitation was a major factor for predicting probability of the presence of nest sites for Egyptian vultures. After identifying potentially suitable habitat, we identified environmental factors affecting landscape-level suitable nesting habitat for Egyptian vultures using generalized linear models. For Egyptian vultures, sites near forests and human settlements were most suitable for nesting, roosting, and foraging, especially in central and western Nepal. Based on potentially suitable nesting habitat and previous work on Egyptian vulture foraging and roosting habitat, we recommend protecting forests near water sources and open areas for their long-term conservation.

1. Introduction

The Egyptian vulture (Neophron percnopterus) is a widely distributed species native to Asia, Africa, and Europe [1]. Egyptian vultures are categorized as endangered by the IUCN Red List of Threatened Species [1]. The estimated global population of Egyptian vultures is 12,400 to 36,000 individuals in the wild [1] and in Nepal; the species is non migratory, a widespread resident with an estimated population of 1000 individuals [2]. The global population of Egyptian vultures is decreasing rapidly (about 50–79% in three generations) including a 10% decline in Europe, 91% in Africa, 40–50% in Israel, and >40% in Greece in last three decades [1]. This species is already extinct from South Africa. Since 1999 this species has declined more than 35% annually in India [1,2]. Major threats include ingestion of the veterinary drug diclofenac, along with human disturbances, lead poisoning, electrocution, collisions, reduced food availability, and habitat change [1,3,4,5]. In Nepal, the Egyptian vulture has been assessed as vulnerable because of its population decline from diclofenac poisoning and other threats, including the destruction of nesting habitat, persecution, and other disturbances [2,6].
Vultures perform numerous beneficial activities including scavenging, which removes carcasses and can control the spread of disease [7]. However, the survival of many vulture species is threatened globally [1]. The causes of population decline vary geographically and include declines in carcass availability, habitat alteration, and poisoning [1]. In Nepal, vultures are currently threatened from forest degradation due to forest resource extraction [2], and disturbances to breeding, roosting, and feeding habitat [1,2]. Consequently, the government of Nepal developed a Vulture Conservation Action Plan in 2009 to prevent the potential loss of vultures, including Egyptian vultures, through the establishment of vulture safe-feeding sites and strengthening vulture-safe zones for providing diclofenac- and poison-free carcasses [8]. However, little information is available on suitable habitats for Egyptian vultures in Nepal, including nesting habitat.
Detailed information on the suitable habitat of Egyptian vultures is limited, and their nests are not uniformly distributed in Nepal [9]. Understanding the factors shaping species distributions and habitat suitability is critical to ensure their long-term persistence [10]. Species distribution models (SDMs) predict species distributions based on relationships between species occurrences and underlying abiotic and biotic factors [11,12]. Maps of species potential distributions generated from SDMs can contribute to conservation planning. For example, species potential distribution maps can help to identify habitat patches that may be connected to create functional meta-populations [11]. Identified suitable areas can also be used for the translocation of threatened species (e.g., the greater one-horned rhino (Rhinoceros unicornis)) [13,14] and to facilitate reserve designs (e.g., the important bird areas proposed for the Iberian Peninsula to conserve Golden Eagles (Aquila chrysaetos)) [15]. In addition, SDMs can be used to predict invasive species occurrences, nesting habitat, and threatened species distributions, the future prediction of species distribution based on climate change [16,17,18,19,20].
Nepal is a mountainous country containing many cliffs that provide a habitat for Egyptian vultures; therefore, we hypothesized that suitable nesting habitat for Egyptian vultures would be widely distributed throughout Nepal. We aimed to estimate the potential distribution of Egyptian vulture nesting habitat in Nepal. These baseline data can be used in support of Egyptian vulture management and conservation, particularly the biodiversity-conservation-related targets in the Vulture Conservation Action Plan (2015–2019), National Biodiversity Strategy and Action Plan of Nepal, and the National Parks and Wildlife Conservation Act 1973.

2. Materials and Methods

2.1. Study Area

Nepal comprises 147,516 km2 (26°22′ to 30°27′ N and 80°04′ to 88°12′ E) and is bordered by China to the north and India to the south, east, and west. Elevations range from 60 to 8850 m above sea level. The climate varies among physiographic regions: the Tarai is considered tropical, the Siwalik subtropical, the Mid Hills and mountains are temperate, and the high mountain consists of subalpine and alpine climates [21]. Varied climatic zones of Nepal support the occurrence of more than 17,097 faunal species, including 886 bird species, of which nine are vulture species [22]. The Madi River corridor occurs in Tanahun, Kaski, and Lumjung districts in the Gandaki Province of Nepal (Figure 1) between the Tarai and High Mountains at an elevation range of 300–8155 m above the sea level. This corridor contains many steppe areas; dominant tree species include needlewood (Schima wallichii), sal (Shorea robusta), chestnut (Castanopsis indica), pine (Pinus wallichiana), banyan (Ficus benghalensis), pipal tree (Ficus religiosa), red cedar (Toona ciliate), gular (Ficus racemosa), red silk cotton (Bombax ceiba), and Indian rosewood (Dalbergia sissoo). In addition to Egyptian vultures, wildlife occurring in the Madi River corridor includes leopards (Panthera pardus), Himalayan black bears (Ursus thibrtanus), Assamese monkeys (Macaca assamensis), golden jackals (Canis aurens), jungle cats (Felis chaus), yellow-throated martens (Martes flavigula), and crab-eating mongooses (Herpestes urva) [23,24]. Overall, 992,110 people inhabit the Tanahun, Kaski, and Lumjung districts of Nepal [25].

2.2. Study Species

The Egyptian vulture is an old-world vulture and a long-lived scavenger. The body measures 47–65 cm from beak tip to tail feathers, with a wingspan about 2.7 times longer than its body length [26]. Adult plumage is white with black flight feathers on the wings and yellow faces; however, sub-adults possess black plumage distinct from adults and other birds. The beak is slender and hook-shaped, adapted for removing meat from carcasses. Males are slightly smaller than females, with body mass ranging from 2.0 to 2.5 kg [27]. Egyptian vulture diets are diverse and include carrion, tortoises, organic waste, insects, young vertebrates, eggs, and feces [28]. Egyptian vultures are solitary nesters and monogamous, with a breeding season from February through May. They generally use the same nest each year, constructed on cliffs and rocky outcrops, but occasionally nest in large trees and utility poles [29,30].

2.3. Data Collection and Management

We obtained 60 locations of active nests of Egyptian vultures throughout Nepal 2006–2021. The nearest distance between nests was 3 km. We obtained nest location data mainly during the breeding season (late February–May). Breeding location of the Egyptian Vultures was checked during field surveys 2020–2021. We used a nest with either eggs or nestlings as an indicator of confirmed breeding. These data were collected during regular vulture-nest monitoring conducted by biologists of the Bird Conservation Nepal, Pokhara Bird Society and supplemented with other research on this species (e.g., [6,31,32] and under the project of Division Forest Kaski, Gandaki Province, Nepal.

2.4. Spatial Modeling

We used species distribution models which envisage potentially suitable habitat for species based on variables supportive to species occurrences [33]. Numerous model types are available for developing SDMs; however, maximum entropy (e.g., Maxent) is one of most effective tools for species having limited spatial presence-only data [34,35,36]. Maximum entropy forecasts the species probability distribution of random events with uniform and high stability [37,38].
We determined potential habitat layers by first selecting 19 bioclimatic variables at 30 arc-second (about 1 km2) resolution from Worldclim [39] and extracted aspect from a digital elevation model (about 1-km2 resolution; Shuttle Radar Topography Mission (SRTM); https://ita.cr.usgs.gov, accessed on 28 November 2022) (Supplementary Table S1). We selected eight environmental variables including mean annual temperature (Bio_1), mean diurnal range (Bio_02), iso-thermality (Bio_3), annual precipitation (Bio_12), precipitation of driest month (Bio_14), precipitation seasonality (Bio_15), precipitation of coldest quarter (Bio_19), and aspect (see Figure 2) based on jackknife analysis, which uses a leave-one-out approach to estimate variable importance (Supplementary Materials Table S1; Figure S1). We excluded highly correlated (|r| > 0.70) variables to reduce the potential for model overfitting [40]. We used maximum entropy modeling software (Maxent version 3.3.3 K) [36] to estimate the potential nesting distribution of Egyptian vultures.
We ran Maxent models using 75% presence data for calibration and the remaining 25% for model validation [41] using the bootstrapping procedure with 10 replications and a maximum of 500 iterations. We used the default settings for the number of background samples and auto features for model construction. We evaluated model validation and accuracy using the area under the curve (AUC) of the receiver operating characteristics [42]. However, because of criticism in the use of AUC for model evaluation [43], we also evaluated model performance using the true skill statistic (TSS), sensitivity, and specificity [44]. The TSS values range from −1 to 1, with values approaching 1 indicating good model performance and values <0 indicating performance no better than random. We converted habitat into binary form (suitable or unsuitable) using the 10th-percentile training presence logistic threshold [45] and used the corresponding raster layer to identify suitable potential nesting habitat for Egyptian vultures.
We used predicted suitable nesting habitat to identify whether the model correctly predicted nesting habitat as well as existing factors influencing the occurrence of Egyptian vulture nests along the Madi River corridor from Damauli, Tanahu to Tanting, Kaski. This area supports abundant Egyptian vultures and we assumed this abundance reflected their residence in this area, where cliffs and forests occur near human settlements. We established 35 stations at 3-km intervals along 135 km of the Madi River. Each station was 500 m × 500 m and within each we established five 100 m × 100 m plots, four at the corners and one at the center. During October 2020–February 2021, we searched each plot for 3 h, recording the occurrence of Egyptian vulture nests and the presence of Egyptian vultures roosting or feeding. We recorded the occurrence of feeding vultures while sitting quietly at the corner of each plot, and for observations of nesting and roosting we searched throughout the plots. From the center of each plot we recorded the elevation using GPS, and measured the distance to the nearest forest, area of agricultural land, water, and human settlement. We measured the nearest distance to human settlement using a GIS if the area had more than one house and were used for living.

2.5. Data Analysis

We used generalized linear models to estimate the effects of distance to forest, distance to agricultural land, distance to nearest water, and distance to human settlement on the presence of nesting and feeding/roosting vultures. We ranked models using the Akaike information criterion adjusted for small samples [46] and Akaike model weights to estimate the relative strength of evidence for each model. We considered models with AICc scores within 4 of the most parsimonious model to have support [46]. We conducted model averaging using all models within 4 AICc of the top model to estimate 95% confidence intervals for each variable and accepted statistical significance at α = 0.05. All analyses were performed in the program R [47].

3. Results

We obtained the locations of 60 nest sites for use in the SDM. Maxent had good performance as indicated by the high AUC (0.952), TSS (0.600), sensitivity (0.869), and specificity (0.730). Annual precipitation (Bio_12) contributed 54%, followed by mean diurnal range (Bio_02: 27%), and mean annual temperature (Bio_1: 8%) (Figure 3).
The greatest probability of Egyptian vulture potential nest habitat suitability for annual precipitation was within the range of 268–3008 mm and annual temperatures of 20–33 °C. (Figure 4). The probability of suitable nesting habitat was greater on east-facing aspects. The probability of suitable Egyptian vulture nesting habitat was greater in central and western regions than other parts of the country and in more mountainous areas, and overall comprised 38,204 km2 (Figure 5A). Suitable areas mostly occurred outside protected areas (Figure 5B).
Among 175 sampling plots in the predicted suitable habitat, Egyptian vultures were sighted at 56 plots (32%). We observed Egyptian vultures feeding (n = 17), soaring (n = 33), and roosting (n = 36), but no nests were detected. The occurrence of Egyptian vultures in their natural habitat was greater near forests (Table 1 and Table 2). The best model included areas near forests and water. Proximity to human settlement also positively influenced the occurrence of Egyptian vultures. The probability of Egyptian vultures’ presence increased with increasing distances to agricultural land and slightly with increasing distance to water.

4. Discussion

The Egyptian vulture’s nesting distribution appears to be influenced by climatic, environmental, and anthropogenic factors. Among climatic factors, annual precipitation contributed the most, which could be a result of abundant precipitation delaying rapid carcass depletion that benefits vulture scavenging [48]. The nests of Egyptian vultures are more common in the central and western mid-hill regions [49,50,51] and high annual precipitation occurs these areas [52]. In these areas, rubbish dumps are common in most municipalities, which undoubtedly attract Egyptian vultures for feeding (Sharma, H.P. personal observation). East-facing slopes also were more likely to have Egyptian vulture nest sites, and increased temperatures from earlier sun exposure could promote thermals to aid flight earlier in the day to facilitate foraging [53].
We were unable to detect nests while searching in the predicted suitable nest area; however, we were able to find foraging and roosting Egyptian vultures in the study area. The majority of suitable habitats were located in the central and western regions, probably due to vulture safe-feeding sites in these regions, which provide diclofenac- and poison-free food to the vultures [8]. The distribution and population of Egyptian vultures inside the protected areas of Nepal is little known. The majority of the protected areas seem less supportive for vulture conservation, so that the government of Nepal declared a vulture safe-feeding zone throughout Nepal for vulture conservation [8]. We found Egyptian vultures selected nest sites close to the forest, likely in part to facilitate nest construction. The Madi River corridor is in the Mid Hills region with a steppe habitat, which might provide visibility for foraging [54]. There is a greater prevalence of steep cliffs in forests of the Mid Hills region (e.g., Kali Gandaki and Madi river corridors) than other regions of Nepal, where vultures construct nests on cliffs and in nearby trees [55,56]. Cliffs can be advantageous to Egyptian vultures due to the greater survival of young in cliff nests and the greater visibility of carrion at lower elevations (Bhusal, K. personal observation). The nesting and roosting of Egyptian vultures on cliffs and trees high above the ground has been previously reported [57]. In addition, a higher probability of the presence of nests near the forest edge could be due to greater diversity and availability of food, including small mammals, young birds, fish, eggs, and decomposed fruits, all of which are consumed by Egyptian vultures [58]. The probability of Egyptian vultures’ nest habitat was not influenced by the distance to water. Forests and cliffs are more prevalent near water, especially rivers (Sharma, H. P., personal observation), and water may be important for vultures to clean the waste from feeding on carcasses and human refuse [59,60]. It is possible that the unimportance of water was an artifact of our study area occurring within a river corridor in which water was not limited.
The distance to agricultural land did not influence the occurrence of Egyptian vultures. In Nepal, people currently bury dead livestock; therefore, carcasses are not available in agricultural land in comparison to forests. That Egyptian vultures were more likely to be present near human settlements is also likely due to greater food availability, including slaughterhouses and dumping sites near settlements in the Mid Hills region [6]. However, human disturbances including roads and associated vehicle traffic can reduce vulture occurrence [49,61].
We acknowledge several limitations in our study. First, we were uncertain whether we detected every vulture or nest that was present or whether multiple observations represented the same individual vulture. However, that we detected multiple nests used in our initial model but none during our assessment for model validation suggests that few, if any, nests were undetected. Further, that an individual bird was observed in multiple plots would not influence our results, as we were addressing population-level responses to various factors that influence vulture presence. Finally, we note that our use of aspect as a continuous variable may have been less desirable than use of alternatives (e.g., cardinal directions), and may require a larger sample size than that which we used.
We developed the first national map for Nepal that estimates areas of suitable nesting habitat for Egyptian vultures. This research identifies areas of potential nesting sites and identifies factors influencing the presence of the species. Our results support the post-1990 distribution of Egyptian vultures, which was estimated to be most limited in east-central and eastern Nepal [2]. Our study also further informs and supports the provisional vulture safe zone, a pioneering landscape conservation approach for vultures in Nepal [62].

5. Conclusions

The nest-site selection of Egyptian vultures was influenced by nearest forests for nest construction and also likely near to human settlement. Based on these factors acting on the occurrence of Egyptian vultures, we have provided the first estimate of its potential persistent nesting distribution in Nepal. Climate, but more likely land-use change, could limit nesting habitat and the overall distribution of this species and threaten its persistence in some areas. Importantly, potential suitable nesting areas occurred predominantly outside of protected areas and we recommend further safeguards to protect the nesting habitat of Egyptian vultures, particularly maintaining forested areas with nearby cliffs, a recommendation within the Vulture Conservation Action Plan of Nepal (2015–2019).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani13040633/s1, Figure S1: Jackknife analysis used to select predictive climactic variables based on their importance to model potential suitable nest habitat of Egyptian vultures, Nepal; Table S1: Predictive variables considered in model to estimate potential distribution of suitable nest habitat for Egyptian vultures, Nepal.

Author Contributions

H.P.S.: designed the study, analyzed the data, and wrote the original draft; S.D.: collected field data and reviewed and edited the draft; K.P.B.: collected field data and reviewed and edited the draft; H.D.: collected field data and reviewed and edited the draft; R.G.: collected field data, and reviewed and edited the draft; A.B.J.: collected field data and reviewed and edited the draft; D.B.R.: collected field data and reviewed and edited the draft; M.G.: collected field data and reviewed and edited the draft; S.G.: collected field data and reviewed and edited the draft; J.L.B.: designed the study and reviewed and edited the draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Division Forest Office, Kaski- Research Permission ref. 836, 31 May 2020.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the Central Department of Zoology, Bird Conservation Nepal, and Pokhara Bird Society for logistic support. We thank the Ishwari Prasad Chaudhary and Hiru Lal Dangaura for field support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. BirdLife International. Neophron percnopterus. The IUCN Red List of Threatened Species; Birdlife Internatiuonal: Cambridge, UK, 2022; e.T22695180A154895845; Available online: https://www.iucnredlist.org/species/22695180/205187871 (accessed on 7 August 2022).
  2. Inskipp, C.; Baral, H.S.; Phuyal, S.; Bhatt, T.R.; Khatiwada, M.; Inskipp, T.; Khatiwada, A.; Gurung, S.; Singh, P.B.; Murray, L.; et al. The Status of Nepal’s Birds: The National Red List Series; Zoological Society of London: London, UK, 2016. [Google Scholar]
  3. Velevski, M.; Nikolov, S.C.; Hallmann, B.; Dobrev, V.; Sidiropoulos, L.; Saravia, V.; Tsiakiris, R.; Arkumarev, V.; Galanaki, A.; Kominos, T.; et al. Population decline and range contraction of the Egyptian Vulture Neophron percnopterus in the Balkan Peninsula. Bird Conserv. Int. 2015, 25, 440–450. [Google Scholar] [CrossRef]
  4. Bounas, A.; Ganoti, M.; Giannakaki, E.; Akrivos, A.; Vavylis, D.; Zorrilla, I.; Saravia, V. First confirmed case of lead poisoning in the endangered Egyptian Vulture (Neophron percnopterus) in the Balkans. Vulture News 2016, 70, 22–29. [Google Scholar] [CrossRef]
  5. Donázar, J.A.; Cortés-Avizanda, A.; Fargallo, J.A.; Margalida, A.; Moleón, M.; Morales-Reyes, Z.; Moreno-Opo, R.; Pérez-García, J.M.; Sánchez-Zapata, J.A.; Zuberogoitia, I.; et al. Roles of raptors in a changing world: From flagships to providers of key ecosystem services. Ardeola 2016, 63, 181–234. [Google Scholar] [CrossRef]
  6. Dhakal, H.; Bhusal, P.K.; Ghimire, S.M. Vultures of Nepal (Nepali Version); Pokhara Bird Society: Pokhara, Nepal, 2019. [Google Scholar]
  7. Ogada, D.L.; Keesing, F.; Virani, M.Z. Dropping dead: Causes and consequences of vulture population declines worldwide. Ann. N. Y. Acad. Sci. USA 2012, 1249, 57–71. [Google Scholar] [CrossRef]
  8. DNPWC. Vulture Conservation Action Plan for Nepal (2015—2019); Department of National Parks and Wildlife Conservation, Ministry of Forests and Soil Conservation, Government of Nepal: Kathmandu, Nepal, 2015; p. 40.
  9. Gurung, S.; Subedi, T.R.; Baral, R.; Pérez-García, J.M.; Ghimire, M.; Baral, H.S.; Virani, M.; Buij, R. Breeding Habitat and Factors Affecting the Cliff Selection by Egyptian Vultures in Central-West Nepal. J. Rapt. Res. 2022, 57, 1–11. [Google Scholar] [CrossRef]
  10. Thuiller, W.; Münkemüller, T.; Lavergne, S.; Mouillot, D.; Mouquet, N.; Schiffers, K.; Gravel, D. A road map for integrating eco-evolutionary processes into biodiversity models. Ecol. Lett. 2013, 16, 94–105. [Google Scholar] [CrossRef]
  11. Guisan, A.; Thuiller, W. Predicting species distribution: Offering more than simple habitat models. Ecol. Lett. 2005, 8, 993–1009. [Google Scholar] [CrossRef]
  12. Guisan, A.; Zimmermann, N.E. Predictive habitat distribution models in ecology. Ecol. Modell. 2000, 135, 147–186. [Google Scholar] [CrossRef]
  13. Dinerstein, E.; Price, L. Demography and Habitat Use by Greater One-Horned Rhinoceros in Nepal. J. Wildl. Manag. 1991, 55, 401–411. [Google Scholar] [CrossRef]
  14. Thapa, K.; Williams, A.C.; Khaling, S.; Bajimaya, S. Observations on habitat preference of translocated rhinos in Bardia National Park and Suklaphanta Wildlife Reserve, Nepal. Pachyderm 2008, 45, 108–114. [Google Scholar]
  15. López-López, P.; García-Ripollés, C.; Soutullo, Á.; Cadahía, L.; Urios, V. Identifying potentially suitable nesting habitat for golden eagles applied to “important bird areas” design. Anim. Conserv. 2007, 10, 208–218. [Google Scholar] [CrossRef]
  16. Kandel, K.; Huettmann, F.; Suwal, M.K.; Ram Regmi, G.; Nijman, V.; Nekaris, K.A.I.; Lama, S.T.; Thapa, A.; Sharma, H.P.; Subedi, T.R. Rapid multi-nation distribution assessment of a charismatic conservation species using open access ensemble model GIS predictions: Red panda (Ailurus fulgens) in the Hindu-Kush Himalaya region. Biol. Conserv. 2015, 181, 150–161. [Google Scholar] [CrossRef]
  17. Adhikari, S.; Sharma, H.P.; Rimal, B.; Belant, J.L.; Katuwal, H.B. Road as a major driver for the potential distribution of the invasive giant African land snail in Nepal. Trop. Ecol. 2020, 61, 583–588. [Google Scholar] [CrossRef]
  18. Sharma, H.P.; Rimal, B.; Zhang, M.; Sharma, S.; Poudyal, L.P.; Maharjan, S.; Kunwar, R.; Kaspal, P.; Bhandari, N.; Baral, L.; et al. Potential distribution of the critically endangered Chinese pangolin (Manis pentadactyla) in different land covers of Nepal: Implications for conservation. Sustainibility 2020, 12, 1282. [Google Scholar] [CrossRef]
  19. Hong, S.H.; Lee, Y.H.; Lee, G.; Lee, D.-H.; Adhikari, P. Predicting Impacts of Climate Change on Northward Range Expansion of Invasive Weeds in South Korea. Plants 2021, 10, 1604. [Google Scholar] [CrossRef] [PubMed]
  20. Baral, K.; Adhikari, B.; Bhandari, S.; Kunwar, R.M.; Sharma, H.P.; Aryal, A.; Ji, W. 2023. Impact of climate change on distribution of common leopard (Panthera pardus) and its implication on conservation and conflict in Nepal. Heliyon 2023, 9, e12807. [Google Scholar] [CrossRef] [PubMed]
  21. LRMP. Land Utilization Report; Land Resource Mapping Project, Kenting Earth Science Canada and Department of Topography, Government of Nepal: Kathmandu, Nepal, 1986; p. 112.
  22. GoN. Nepal’s Sixth National Report to the Convention on Biological Diversity; Ministry of Forests and Environment: Kathmandu, Nepal, 2018.
  23. Jnawali, S.R.; Baral, H.; Lee, S.; Acharya, K.; Upadhyay, G.; Pandey, M.; Griffiths, J. The Status of Nepal Mammals: The National Red List Series; Department of National Parks and Wildlife Conservation Kathmandu: Kathmandu, Nepal, 2011.
  24. Sharma, H.P.; Adhikari, B.; Bhandari, S.; Baral, K.; Kunwar, R.M. Crab-Eating Mongoose Herpestes urva: Occurrence and its Activity in Mid-Hills of Nepal. J. Inst. Sci. Technol. 2022, 26, 53–60. [Google Scholar] [CrossRef]
  25. CBS. National Population and Housing Census 2011 (National Report); Government of Nepal, National Planning Commission Secretariat Central Bureau of Statistics: Kathmandu, Nepal, 2012.
  26. Ferguson-Lees, J.; Christie, D.A. Egyptian Vulure. Raptors of the World; Christopher Helm A and C Black Publications Ltd.: London, UK, 2001; pp. 417–424. [Google Scholar]
  27. Donázar, J.A.; Palacios, C.J.; Gangoso, L.; Ceballos, O.; González, M.J.; Hiraldo, F. Conservation status and limiting factors in the endangered population of Egyptian vulture (Neophron percnopterus) in the Canary Islands. Biol. Conserv. 2002, 107, 89–97. [Google Scholar] [CrossRef]
  28. Margalida, A.; Benitez, J.R.; Sanchez-Zapata, J.A.; Ávila, E.; Arenas, R.; Donázar, J.A. Long-term relationship between diet breadth and breeding success in a declining population of Egyptian Vultures Neophron percnopterus. Ibis 2012, 154, 184–188. [Google Scholar] [CrossRef]
  29. Sarà, M.; Di Vittorio, M. Factors influencing the distribution, abundance and nest-site selection of an endangered Egyptian Vulture (Neophron percnopterus) population in Sicily. Anim. Conserv. 2003, 6, 317–328. [Google Scholar] [CrossRef]
  30. Naoroji, R. Birds of Prey of the Indian Subcontinent; Christopher Helm: London, UK, 2006. [Google Scholar]
  31. Dhakal, H.; Sharma, H.P.; McClure, C.J.; Virani, M.; Rolek, B.W.; Pradhan, N.M.; Bhusal, K.P. Vulture distribution and people perception of vultures in Pokhara Valley, Nepal. Ecol. Evol. 2022, 12, e8528. [Google Scholar] [CrossRef]
  32. Gautam, R.; Baral, N.; Sharma, H.P. Are white-rumped vultures (Gyps bengalensis) scavengers or predators at a vulture safe feeding site of Nepal? Biodiv. J. Biol. Diver. 2022, 23, d230757. [Google Scholar] [CrossRef]
  33. Hirzel, A.H.; Hausser, J.; Chessel, D.; Perrin, N. Ecological-niche factor analysis: How to compute habitat-suitability maps without absence data? Ecology 2002, 83, 2027–2036. [Google Scholar] [CrossRef]
  34. Phillips, S.J.; Dudík, M.; Schapire, R.E. A maximum entropy approach to species distribution modeling. In Proceedings of the Twenty-First International Conference on Machine Learning, Banff, AB, Canada, 4–8 July 2004; p. 83. [Google Scholar]
  35. Phillips, S.J.; Dudík, M. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography 2008, 31, 161–175. [Google Scholar] [CrossRef]
  36. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 2011, 17, 43–57. [Google Scholar] [CrossRef]
  37. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 2006, 190, 231–259. [Google Scholar] [CrossRef]
  38. Duan, R.Y.; Kong, X.Q.; Huang, M.Y.; Fan, W.Y.; Wang, Z.G. The predictive performance and stability of six species distribution models. PLoS ONE 2014, 9, e112764. [Google Scholar] [CrossRef]
  39. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  40. Peterson, A.T. Ecological niche conservatism: A time-structured review of evidence. J. Biogeogr. 2011, 38, 817–827. [Google Scholar] [CrossRef]
  41. Fielding, A.H.; Bell, J.F. A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ. Conserv. 1997, 24, 38–49. [Google Scholar] [CrossRef]
  42. Hajian-Tilaki, K. Receiver operating characteristic (ROC) curve analysis for medical diagnostic test evaluation. Casp. J. Intern. Med. 2013, 4, 627–635. [Google Scholar]
  43. Lobo, J.M.; Jiménez-valverde, A.; Real, R. AUC: A misleading measure of the performance of predictive distribution models. Glob. Ecol. Biogeog. 2008, 17, 145–151. [Google Scholar] [CrossRef]
  44. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  45. Liu, C.; Newell, G.; White, M. On the selection of thresholds for predicting species occurrence with presence-only data. Ecol. Evol. 2016, 6, 337–348. [Google Scholar] [CrossRef]
  46. Anderson, D.R.; Burnham, K.P. Avoiding pitfalls when using information-theoretic methods. J. Wildl. Manag. 2002, 66, 912–918. [Google Scholar] [CrossRef]
  47. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022; Available online: https://www.R-project.org/ (accessed on 7 August 2022).
  48. Selva, N.; Jedrzejewska, B.; Jedrzejewski, W.; Wajrak, A. Scavenging on European bison carcasses in Bialowieza primeval forest (eastern Poland). Ecoscience 2003, 10, 303–311. [Google Scholar] [CrossRef]
  49. Ghimire, P.; Dhakal, H.; Sharma, B.; Ghimire, M.; Baral, M.; Bhusal, K.P. Status of Egyptian Vulture Neophron percnopterus in Pokhara Valley, Nepal. Vulture Bull. 2020, 9, 1–4. [Google Scholar]
  50. Bhusal, K.P.; Pandey, M.; Magar, T.G. Diversity and Status of Birds in Argha Important Bird and Biodiversity Area, Western Midhill of Nepal. Danphe 2020, 29, 1–11. [Google Scholar]
  51. KC, K.B.; Koju, N.P.; Bhusal, K.P.; Low, M.; Ghimire, S.K.; Ranabhat, R.; Panthi, S. Factors influencing the presence of the endangered Egyptian Vulture Neophron percnopterus in Rukum, Nepal. Glob. Ecol. Conserv. 2019, 20, e00727. [Google Scholar] [CrossRef]
  52. Kansakar, S.R.; Hannah, D.M.; Gerrard, J.; Rees, G. Spatial pattern in the precipitation regime of Nepal. Int. J. Climatol. J. Royal Meteorol. Soc. 2004, 24, 1645–1659. [Google Scholar] [CrossRef]
  53. Sen, B.; Tavares, J.P.; BİLGİN, C.C. Nest site selection patterns of a local Egyptian Vulture Neophron percnopterus population in Turkey. Bird Conserv. Int. 2017, 27, 568–581. [Google Scholar] [CrossRef]
  54. Sen, B. Breeding Ecology of the Egyptian Vulture (Neophron percnopterus) Population Beypazarı. Master’s Thesis, Middle East Technical University, Ankara, Turkey, 2012. [Google Scholar]
  55. Ceballos, O.L.; Donázar, J.A. Roost-tree characteristics, food habits and seasonal abundance of roosting Egyptian Vultures in northern Spain. J. Rapt. Res. 1990, 24, 19–25. [Google Scholar]
  56. Ceballos, O.; Donázar, J.A. Factors influencing the breeding density and nest-site selection of the Egyptian Vulture (Neophron percnopterus). J. Ornithol. 1989, 130, 353–359. [Google Scholar] [CrossRef]
  57. Subedi, T.R.; DeCandido, R. Population and breeding success of Red-headed Vulture Sarcogyps calvus and Egyptian Vulture Neophron percnopterus in central west Nepal. Vulture News 2014, 67, 21–32. [Google Scholar] [CrossRef]
  58. Milchev, B.; Spassov, N.; Popov, V. Diet of the Egyptian Vulture (Neophron percnopterus) after livestock reduction in eastern Bulgaria. North-West. J. Zool. 2012, 8, 315–323. [Google Scholar]
  59. Kanaujia, A.; Kushwaha, S. Vulnerable Vultures of India: Population, Ecology and Conservation; Rare Animals of India, Bentham Science Publishers: Sharjah, UAE, 2013; pp. 113–144. [Google Scholar]
  60. Buechley, E.R.; Sekercioglu, C.H. Vultures. Cur. Biol. 2016, 26, R560–R561. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Margalida, A.; García, D.; Cortés-Avizanda, A. Factors influencing the breeding density of Bearded vultures, Egyptian Vultures and Eurasian griffon vultures in Catalonia (NE Spain): Management implications. Anim. Biodiv. Conserv. 2007, 30, 189–200. [Google Scholar]
  62. Bhusal, K.P. Vulture Safe Zone: A landscape level approach to save the threatened vultures in Nepal. Himal. Nat. 2018, 1, 25–26. [Google Scholar]
Figure 1. Madi River corridor and tributaries where field study was carried out after identifying potential suitable nesting habitat of Egyptian vultures in Tanahun, Kaski, and Lamjung districts, Nepal.
Figure 1. Madi River corridor and tributaries where field study was carried out after identifying potential suitable nesting habitat of Egyptian vultures in Tanahun, Kaski, and Lamjung districts, Nepal.
Animals 13 00633 g001
Figure 2. Spearman pairwise correlation coefficients between predictive variables aspect, mean annual temperature (bio_01), mean diurnal range (bio_02), iso-thermality (bio_03), temperature seasonality (bio_04), maximum temperature of the warmest month (bio_05), minimum temperature of coldest month (bio_06), temperature annual range (bio_07), mean temperature of wettest quarter (bio_08), mean temperature of driest quarter (bio_09), mean temperature of the warmest quarter (bio_10), mean temperature of coldest quarter (bio_11), annual precipitation (bio_12), precipitation of wettest month (bio_13), precipitation of driest month (bio_14), precipitation seasonality (bio_15), precipitation of driest quarter (bio_16), precipitation of wettest quarter (bio_17), precipitation of the warmest quarter (bio_18) and precipitation of the coldest quarter (bio_19) considered in the Egyptian vulture nesting habitat distribution model, Nepal. Color variation indicates correlation values.
Figure 2. Spearman pairwise correlation coefficients between predictive variables aspect, mean annual temperature (bio_01), mean diurnal range (bio_02), iso-thermality (bio_03), temperature seasonality (bio_04), maximum temperature of the warmest month (bio_05), minimum temperature of coldest month (bio_06), temperature annual range (bio_07), mean temperature of wettest quarter (bio_08), mean temperature of driest quarter (bio_09), mean temperature of the warmest quarter (bio_10), mean temperature of coldest quarter (bio_11), annual precipitation (bio_12), precipitation of wettest month (bio_13), precipitation of driest month (bio_14), precipitation seasonality (bio_15), precipitation of driest quarter (bio_16), precipitation of wettest quarter (bio_17), precipitation of the warmest quarter (bio_18) and precipitation of the coldest quarter (bio_19) considered in the Egyptian vulture nesting habitat distribution model, Nepal. Color variation indicates correlation values.
Animals 13 00633 g002
Figure 3. Percent contribution of climactic variables aspect, mean annual temperature (bio_01), mean diurnal range (bio_02), iso-thermality (bio_03), annual precipitation (bio_12), precipitation of wettest month (bio_13), precipitation of driest month (bio_14), precipitation seasonality (bio_15) and precipitation of the coldest quarter (bio_19) for model construction for the potential distribution of suitable nest habitat for Egyptian vultures, Nepal.
Figure 3. Percent contribution of climactic variables aspect, mean annual temperature (bio_01), mean diurnal range (bio_02), iso-thermality (bio_03), annual precipitation (bio_12), precipitation of wettest month (bio_13), precipitation of driest month (bio_14), precipitation seasonality (bio_15) and precipitation of the coldest quarter (bio_19) for model construction for the potential distribution of suitable nest habitat for Egyptian vultures, Nepal.
Animals 13 00633 g003
Figure 4. Response curves of aspect, mean annual temperature (bio_1), mean diurnal range (bio_02), iso-thermality (bio_3), annual precipitation (bio_12), precipitation of driest month (bio_14), precipitation seasonality (bio_15) and precipitation of coldest quarter (bio_19) on the probability of potential suitable nest habitat for Egyptian vultures, Nepal.
Figure 4. Response curves of aspect, mean annual temperature (bio_1), mean diurnal range (bio_02), iso-thermality (bio_3), annual precipitation (bio_12), precipitation of driest month (bio_14), precipitation seasonality (bio_15) and precipitation of coldest quarter (bio_19) on the probability of potential suitable nest habitat for Egyptian vultures, Nepal.
Animals 13 00633 g004
Figure 5. (A) Potential suitable nest habitat of Egyptian vultures, Nepal. (B) Potential suitable nest habitat of Egyptian vultures and protected areas, Nepal.
Figure 5. (A) Potential suitable nest habitat of Egyptian vultures, Nepal. (B) Potential suitable nest habitat of Egyptian vultures and protected areas, Nepal.
Animals 13 00633 g005
Table 1. Outputs from generalized linear models to identify factors related to the probability of presence of Egyptian vultures, Madi River corridor, Nepal. Egyptian vulture (presence) used as response variable; elevation (m); distance to forest (m), distance to agricultural land (m), distance to human settlement (m), distance to water (m) are predictive variables. K is the number of parameters, ∆AICc is the difference between the AICc value of the best-supported model and successive models, and wi is the Akaike model weight. Only models with ∆AICc < 4 are reported.
Table 1. Outputs from generalized linear models to identify factors related to the probability of presence of Egyptian vultures, Madi River corridor, Nepal. Egyptian vulture (presence) used as response variable; elevation (m); distance to forest (m), distance to agricultural land (m), distance to human settlement (m), distance to water (m) are predictive variables. K is the number of parameters, ∆AICc is the difference between the AICc value of the best-supported model and successive models, and wi is the Akaike model weight. Only models with ∆AICc < 4 are reported.
ParametersKlogLikAICc∆AICcwi
Distance to forest + Distance to water3−37.0880.200.000.12
Distance to forest + Distance to water + Distance to settlement4−36.2680.500.360.10
Distance to forest + Distance to water + Distance to settlement + Elevation5−35.3280.600.460.09
Distance to forest + Distance to water + Elevation4−36.5281.000.870.08
Distance to agricultural land + Distance to forest + Distance to water + Distance to settlement5−35.9781.901.760.05
Distance to forest + Elevation3−37.9882.001.790.05
Distance to water2−39.0382.001.860.05
Distance to agricultural land + Distance to forest + Distance to water4−37.0782.101.970.04
Distance to forest + Distance to settlement + Elevation4−37.1782.302.180.04
Distance to agricultural land + Distance to forest + Distance to water + Distance to settlement + Elevation6−35.2682.502.360.04
Distance to water + Elevation3−38.3482.702.510.03
Distance to agricultural land + Distance to forest + Distance to water + Elevation5−36.3682.702.540.03
Distance to water + Distance to settlement + Elevation4−37.4883.002.780.03
Distance to water + Distance to settlement3−38.5183.002.850.03
Distance to forest2−39.6783.303.160.02
Distance to agricultural land + Distance to forest + Elevation4−37.7383.503.290.02
Elevation2−39.7883.603.390.02
Distance to agricultural land + Distance to water3−39.0184.003.840.02
Distance to agricultural land + Distance to water + Distance to settlement4−38.0384.103.890.02
Table 2. Model-averaged parameter estimates and lower and upper 95% confidence limits describing the occurrence of Egyptian vultures in Nepal. Parameter estimates were averaged from all models reported in Table 1.
Table 2. Model-averaged parameter estimates and lower and upper 95% confidence limits describing the occurrence of Egyptian vultures in Nepal. Parameter estimates were averaged from all models reported in Table 1.
VariablesEstimateSELower CLUpper CLzp
(Intercept)1.3670.794−0.1902.100861.7210.085
Distance to forest−0.0050.002−0.0090.000041.9760.048
Distance to water0.0040.003−0.0010.01011.5330.125
Distance to settlement−0.0020.002−0.0050.00101.2950.195
Elevation0.0010.001−0.0010.00291.3440.179
Distance to agricultural land0.00030.002−0.0040.00490.1180.906
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sharma, H.P.; Dhakal, S.; Bhusal, K.P.; Dhakal, H.; Gautam, R.; Joshi, A.B.; Rana, D.B.; Ghimire, M.; Ghimire, S.; Belant, J.L. Factors Influencing the Potential Distribution of Globally Endangered Egyptian Vulture Nesting Habitat in Nepal. Animals 2023, 13, 633. https://doi.org/10.3390/ani13040633

AMA Style

Sharma HP, Dhakal S, Bhusal KP, Dhakal H, Gautam R, Joshi AB, Rana DB, Ghimire M, Ghimire S, Belant JL. Factors Influencing the Potential Distribution of Globally Endangered Egyptian Vulture Nesting Habitat in Nepal. Animals. 2023; 13(4):633. https://doi.org/10.3390/ani13040633

Chicago/Turabian Style

Sharma, Hari Prasad, Santosh Dhakal, Krishna Prasad Bhusal, Hemanta Dhakal, Ramji Gautam, Ankit Bilash Joshi, Deu Bahadur Rana, Manshanta Ghimire, Suman Ghimire, and Jerrold L. Belant. 2023. "Factors Influencing the Potential Distribution of Globally Endangered Egyptian Vulture Nesting Habitat in Nepal" Animals 13, no. 4: 633. https://doi.org/10.3390/ani13040633

APA Style

Sharma, H. P., Dhakal, S., Bhusal, K. P., Dhakal, H., Gautam, R., Joshi, A. B., Rana, D. B., Ghimire, M., Ghimire, S., & Belant, J. L. (2023). Factors Influencing the Potential Distribution of Globally Endangered Egyptian Vulture Nesting Habitat in Nepal. Animals, 13(4), 633. https://doi.org/10.3390/ani13040633

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop