Vulture distribution and people perception of vultures in Pokhara Valley, Nepal

Abstract Due to an abundance and diversity of vultures, Nepal is one of the most important countries for vulture conservation. Within Nepal, the Pokhara Valley is especially significant. We examine the distribution of vultures within the Pokhara Valley by conducting counts at 11 potential feeding or roosting sites using point count method. We further surveyed people of the valley regarding their perception of vulture ecology and conservation, knowledge of diclofenac use within the valley, and burial of livestock carcasses. We detected eight species of vultures, four of which are currently threatened with extinction. White‐rumped vulture Gyps bengalensis, Egyptian vulture Nephron percnopterus, and Himalayan vulture G. himalayensis were the most abundant. Almost all respondents (98%) had sighted the vultures in the wild. Formally educated respondents reported seeing vultures’ slightly more than nonformally educated respondents. Fifty‐eight percent respondents suspected habitat loss was the major threat for the vulture population decline in Pokhara Valley, and 97% respondents were not aware of any diclofenac use. The knowledge of vultures in people with different age groups suggests a more awareness programs are needed for local people, especially those who carry out animal husbandry and provide livestock to the vulture restaurant.

across South Asia were decimated by veterinary use of nonsteroidal anti-inflammatory drugs (NSAIDs; Margalida & Ogada, 2018;Ogada et al., 2012;Pain et al., 2008). In particular, diclofenac caused catastrophic declines of >95% in three species of Gyps vultures after being introduced in the 1990s (Oaks et al., 2004;Prakash, 1999;Prakash et al., 2003Prakash et al., , 2012. In response to these population declines, the governments of India, Pakistan, and Nepal banned the importation, production, and sale of veterinary diclofenac in 2006, with Bangladesh following in 2010 (Margalida & Ogada, 2018;Pain et al., 2008). Conservationists also promoted nontoxic meloxicam as an alternative to diclofenac (Swan et al., 2006;Swarup et al., 2007) and established Vulture Safe feeding sites (Vulture Restaurants) where diclofenac-free carcasses are provisioned as food for vultures .
Importantly, established Vulture Safe Zones also requires working with local communities to remove stocks of diclofenac, educate potential diclofenac users about its deleterious effects, and monitor the potential drug users in surrounding areas (Chaudhary et al., 2010. Following these conservation actions, some populations of Gyps vultures in South Asia seem to have stabilized from 2012 onwards (Chaudhry et al., 2012;Galligan et al., 2019;Prakash et al., 2012).
Despite some encouraging results, the threat of NSAIDs remains (Cuthbert, Dave, et al., 2011;Cuthbert, Taggart, et al., 2011;Margalida & Ogada, 2018). Illegal use of diclofenac in veterinary medicine still occurs to varying degrees across Africa and South Asia (Botha et al., 2017;Cuthbert, Dave, et al., 2011;Cuthbert, Taggart, et al., 2011), and several legal NSAIDs are toxic to vultures (Cuthbert et al., , 2016Naidoo et al., 2010). Less than one percent of livestock carcasses need to be contaminated with diclofenac to cause catastrophic losses observed during the Asian Vulture Crisis (Green et al., 2004). Actions including establishing Vulture Safe Zones and educating local people of the dangers of diclofenac therefore remain conservation priorities. In fact, according to IUCN Red List assessments, "education and awareness" and "law and policy" are the two most pressing conservation needs for Old World vultures (McClure et al., 2018)-highlighting that both legislative and educational actions are needed to assuage the threat of poisoning (Parvanov et al., 2018 Bhusal et al., 2019;Dhakal et al., 2019;DNPWC, 2015). Nepal is one of the most important areas on Earth for raptor conservation, especially in the context of vulture declines (Prakash et al., 2012). This relatively small nation contains six of the 10 most critical Important Bird and Biodiversity Areas (IBAs) for raptors and contains the 10th greatest number of declining raptor species of all nations (42 species;McClure et al., 2018). Indeed, the concept of Vulture Safe Zones was first successfully implemented in Nepal (Chaudhary et al., 2010). The Kaski District of Nepal might be particularly important-containing five ecoregions (Figure 1; Olson et al., 2001) and the Annapurna Conservation Area, which is the sixth most important IBA for raptors, globally (McClure et al., 2018). This district has one Vulture Safe Feeding Site at Ghachowk. Because of such diversity of vultures with this small area, it is important that conservationists understand the local people's perception toward vultures.
Here, we report the results of questionnaires completed by residents of the Kaski District of Nepal regarding their knowledge of vulture biology and ecosystem services, nearby use of diclofenac, and recent population status. We further report counts of eight vulture species, of which four are currently threatened with extinction (BirdLife International, 2019; Dhakal et al., 2019;Inskipp et al., 2016). We specifically emphasize local knowledge of the vulture (see Table 1 for scientific names) because these are rare , critically endangered This information is pre-requisite to update the existing conservation plan and to develop site-specific management plan at one of the IBA in Nepal. Therefore, we aimed to provide site-specific vulture abundance and people's knowledge on vulture conservation near to their settlement.

| Study area
We focused our study on the Pokhara Valley (28°7′N to 28°12′N,

| Vulture
We performed fixed-point surveys following Huff (2000) at 11 sites that we chose based on potential roosting and feeding habitat, and these sites cover the whole study area which comprises approximately 250 km 2 . The distance between first and last site was approximately 30 km apart and minimum distance between two sites was around 2 km. We recorded all individual vultures soaring, feeding, and roosting. We performed all surveys from vantage points that provided views of the entire sur- We communicated each observer if the soaring vulture seems potentially overlapping between nearest sites to avoid the double count of individual. We recorded the number of each vulture species at every 30-min interval, and we used the record of highest number of each vulture species from each station for data analysis. Surveys were conducted between 07h30 and 10h30 in the morning and 16h-18h in the evening during summer and 08h-11h30 in the morning and 15h-17h in the evening during winter.

| Questionnaire survey
Using a questionnaire (Appendix S1), we interviewed randomly selected respondents, who were living near to study area, that is, within 1 km periphery for questionnaire survey at each of the 11 vulture monitoring sites. We interviewed with 300 individuals, and the number of interviewee was varied according to sites (8-30 household) based on the household available within the 1 km periphery of study site. The interviewee was selected creating random number.
The selected sample size for interview at each site was based on the 5% margin of the error at 95% confidence interval (Kreb, 2014).
We surveyed only one adult (>16 years old) per household. We did not discriminate to the people based on education level, gender, ethnicity, or religion. We collected demographic data such as age, gender, and education. Furthermore, the questionnaire was designed to assess perception and knowledge of vultures as well as livestock holding practices including use of NSAIDS and methods of carcass disposal. In addition to the questionnaire, we showed pictures of different species of birds including all vulture species of Nepal, some eagle, and storks and asked them to identify which ones were vultures. The conduct of all aspects of this study was approved by the Department of National Parks and Wildlife Conservation, Nepal (permit DNPWC-074/75-2303).

| Analysis
We implemented state-space models in a Bayesian framework to estimate relative abundance of vulture species observed at each of our survey sites (Benson & McClure, 2019). A Bayesian framework allowed us to implement state-space models while accounting for effects from survey effort on vulture counts. State-space models can separate process from observation error in count data and can estimate imperfectly observed latent processes such as population  (Kéry & Schaub, 2012). These state-space models accounted for imperfect detection and misclassification of species, but assumed an equal probability of these errors. Therefore, we estimated an improved yet imperfect index of relative abundance rather than counts.
We modeled monthly relative abundance (N) as a Markovian process: log (N t+1,s ) = log (N t,s ) + r t,s (Kéry & Schaub, 2012), where r is the stochastic change in vulture counts and can be interpreted as the monthly change in abundance over time (t) at each survey site (Ss).
We included effort as a scaled and centered covariate (mean = 0 and SD = 1) to standardize counts to the mean time spent during sur- We implemented models using JAGS (Plummer, 2003) and the package jagsUI (Kellner, 2016) in R (R Core Team, 2017) and used three chains with 300,000 iterations, burn-in of 150,000, adaptation of 20,000, and we thinned one out of every 50 posterior draws resulting in three chains each having 3,000 posterior draws. We calculated the Gelman-Rubin statistic (R; Gelman & Rubin, 1992) for each parameter and determined convergence of chains when R < 1.1, and we created trace plots of parameter chains to assess convergence of chains. To further assess our models, we calculated the ratio of observation error to process error (i.e., the error ratio; Auger-Méthé et al., 2016) and considered our models to be adequate when the probability that the error ratio >10 was >0.90. We used vague priors for all parameters because these priors produce similar To compare species abundances, we calculated relative abundance across all months and study sites. To estimate an index of relative abundance that included all vulture species at a survey site, we summed each iteration of the posterior draws across all species. We present the median, 2.5th, and 97.5th percentiles for all reported estimates.
We examined the differences between education levels, age groups, and sexes regarding awareness of diclofenac, knowledge on the vulture species, and its ecosystem services. Our data were not normally distributed; therefore, we used Fisher's exact test and a Kruskal-Wallis tests for binary and numeric questionnaire responses, respectively. We performed all analyses in R program (R Core Team, 2017).

| RE SULTS
We detected eight species of vultures, of which Egyptian vulture and Bearded vulture were detected the most and least often, respectively ( Table 1). The species with the highest relative abundance per survey across sites were White-rumped vulture (median = 8), Egyptian vulture (median = 6), and Himalayan vulture (median = 3) ( Figure 2, Table 1, Table S1). Of the species analyzed, the three with the lowest relative abundance were Cinereous vulture, Griffon vulture, and Slender-billed Vulture (Figure 2, Figure   S1, Table S1). There was, however, substantial variance in relative abundance across sites for some species. For example, the Egyptian vulture was especially abundant at a few sites, leading to particularly skewed distributions of their survey-wide relative abundance ( Figure 2, Figure S1, Table 1). The three sites with the highest overall relative abundance were the Landfill, Dovilla, and Ghachowk,  Figure S1, Table   S1). No species had B values with positive 2.5th percentiles ( Table   S2). The probability that the error ratio >10 was <0.90 for all species (Table S3).
We interviewed to 300 respondents living near to study sites. In the interview, we did not find differences in age between male and female respondents (Kruskal-Wallis test = 1.269, df = 1, p = .26); however, there was variation in the age according to education (Kruskal-Wallis test = 47.962, df = 1, p = .001) in respondents (Table 2). In addition, there were differences in formally educated respondents between gender (Fisher's exact test, two-tailed = 0.015) ( Table 2). Almost all respondents (>95%) identified the vulture's pictures during this questionnaire survey.
Ninety-eight percent of respondents reported seeing wild vultures. Among these, formally educated respondents have reported seeing vultures slightly more than nonformally educated respondents (Table 2, Figure 3a,b). Only seven percent of respondents, among which 86% were formally educated, had seen a vulture nest in the wild (Figure 4a,b). Almost all formally educated and noneducated respondents (98%) confirmed the usefulness of vultures in an ecosystem services (Figure 5a,b). Among them, 86% of respondents recognized vultures as important carcass removal agents for environment cleaning, followed by 12% as preventing spread of disease, and two percent not sure. Only three percent of respondents were aware of the negative effects of diclofenac use on vulture populations (Figure 6a,b). Ninety-six percent of respondents believe vulture populations are declining in Pokhara, Kaski  success . We found one nest with a successfully reared nestling in 2017. However, population trends for long-lived species such as raptors are most sensitive to changes in adult survival (Clark & Martin, 2007;McClure et al., 2017;Saether & Bakke, 2000).

Comparison Statistics
These areas seem a good habitat for vulture occurrence, potentially due to food accessibility near to landfill sites, rivers, and vulture safe feeding site. However, the area is suffering from habitat loss and unsure about the illegal use of Diclofenac. In addition, vultures potentially threats from high collision risk with aircrafts. Pokhara Regional International Airport will start soon its fight from 2022, and vulture movements, nesting, and soaring sites are located within 2 km periphery of Airport. Every year >300 individuals of vultures were re- Knowledge and perception of a species often leads to support for conservation (Bhattarai & Fischer, 2014;Katuwal et al., 2021;Sharma et al., 2019). Therefore, if people understand the usefulness of vultures as an important scavenger (Ballejo et al., 2019;García-Alfonso et al., 2019), they might take initiation for conservation. Mainali for their support in every step of research.

CO N FLI C T O F I NTE R E S T
None.

F I G U R E 7
Respondents view whether vulture population in the wild in increasing or decreasing in the study area, (a) formally educated respondents and (b) non-formally educated respondents