Gharial (Gavialis gangeticus) conservation in Bardia National Park, Nepal: Assessing population structure and habitat characteristics along the river channel amidst infrastructure development

Abstract Nepal initiated numerous hydropower and irrigation‐related infrastructure projects to enhance and promote green energy, water security, and agricultural productivity. However, these projects may pose risks to natural habitats and the well‐being of aquatic fauna, leading to significant effects on delicate ecosystems. To understand these potential impacts, it is crucial to gather reliable baseline data on the population status and habitat characteristics of species. This study specifically focuses on Gharials (Gavialis gangeticus), a critically endangered species. We recorded data on pre‐determined habitat variables at stations spaced 500 m apart along the two major river streams of Bardia National Park, as well as at locations where Gharials were sighted between February and March 2023. We used binary logistic regression with a logit link function to investigate the habitat characteristics related to the occurrence of Gharials. The presence/absence of Gharials at sampling points served as the dependent variable, while 10 other predetermined variables (ecological variables and disturbance variables) served as independent variables. Our study recorded 23 Gharials, comprising 14 adults, six sub‐adults, and three juveniles, with a sex ratio of 55.56 males per 100 females. Most individuals (83%) were found basking. Among the 10 habitat predictors, three variables (mid‐river depth, river width, and water temperature) were significantly correlated (p < .05) with the probability of Gharial occurrence. The model shows that Gharial detection probability increases with greater mid‐river depth and width and lower water temperature. This study establishes a population baseline for Gharials within the river system before the construction of large infrastructure projects, such as dams and irrigation canals. It also recommends continuous monitoring of Gharial populations after water release and/or diversion to evaluate the impact of large infrastructure projects on the population and their associated habitat characteristics. This will help enable more informed and targeted conservation efforts.

tion status and habitat characteristics of species.This study specifically focuses on Gharials (Gavialis gangeticus), a critically endangered species.We recorded data on pre-determined habitat variables at stations spaced 500 m apart along the two major river streams of Bardia National Park, as well as at locations where Gharials were sighted between February and March 2023.We used binary logistic regression with a logit link function to investigate the habitat characteristics related to the occurrence of Gharials.The presence/absence of Gharials at sampling points served as the dependent variable, while 10 other predetermined variables (ecological variables and disturbance variables) served as independent variables.Our study recorded 23 Gharials, comprising 14 adults, six sub-adults, and three juveniles, with a sex ratio of 55.56 males per 100 females.Most individuals (83%) were found basking.Among the 10 habitat predictors, three variables (mid-river depth, river width, and water temperature) were significantly correlated (p < .05)with the probability of Gharial occurrence.The model shows that Gharial detection probability increases with greater mid-river depth and width and lower water temperature.This study establishes a population baseline for Gharials within the river system before the construction of large infrastructure projects, such as dams and irrigation canals.It also recommends continuous monitoring of Gharial populations after water release and/or diversion to evaluate the impact of large infrastructure projects on the population and their associated habitat characteristics.This will help enable more informed and targeted conservation efforts.

| INTRODUC TI ON
Freshwater ecosystems, covering only 1% of the Earth's surface, host a third of all vertebrate species worldwide (Balian et al., 2008).However, the persistent of intense anthropogenic pressures on freshwater ecosystems worldwide have placed the rich biodiversity within them in jeopardy (Dudgeon et al., 2006;Reid et al., 2019), with a heightened focus on megafauna species.
These larger members of their taxonomic groups face a significantly greater risk of population declines and potential extinctions compared to their smaller counterparts (He et al., 2018).This situation is expected to persist or worsen in the coming years (He et al., 2019).Main threats encompass excessive exploitation, the erection of dams, habitat deterioration, pollution, and the incursion of invasive species (He et al., 2017).For instance, dams not only disrupt the connectivity of river channels but also exert significant impacts on the overall composition of riverine landscapes (Devkota et al., 2020;Nilsson et al., 2005;Nilsson & Berggren, 2000).River-dwelling habitat specialists encounter imminent and direct threats, including elevated harvesting pressures, constrained foraging opportunities, and the depletion of aquatic environments due to factors such as disrupted channel connectivity, altered flow regimes, and water extraction projects (Vashistha et al., 2021).If the impacts stemming from freshwater infrastructure developments, such as the construction of dams and irrigation canals, are left unattended, they have the potential to trigger declines in species populations, habitat fragmentation, and localized extinctions [e.g., Indus River dolphins (Braulik et al., 2014)] and regional extinction [e.g., Gharial in the Indus (Lang et al., 2019)].Thus, establishing baseline estimates for key biological parameters such as population size, distribution range, and habitat characteristics is however critical for directing conservation planning efforts towards threatened and rare taxa (Jeliazkov et al., 2022;Sutton et al., 2023), particularly those that are critically endangered (Dhami et al., 2023).Current conservation measures are inadequate in protecting freshwater habitats and biodiversity as these ecosystems are typically not a priority in conservation management plans and actions (Abell et al., 2017;Collen et al., 2014;Darwall et al., 2018;Harrison et al., 2018;Pant et al., 2021).
Nepal has experienced a greater loss of freshwater species compared to terrestrial species, with a nearly two-fold difference over the past five decades (WWF, 2016).Gharial (Gavialis gangeticus; see Figure 1) is a notable example among the freshwater reptiles at high risk of facing a rapid decline in Nepal's freshwater ecosystems.It is the only extant species of the family Gavialidae (Maskey, 1989) 1973(GoN, 1973) ) given its conservation status in the country.
Prior to the 1940s, Gharials thrived in significant numbers, estimated between 5000 and 10,000 individuals, and inhabited major river systems spanning from the Indus in Pakistan to the Gangetic floodplains of India, Nepal, Bangladesh, Bhutan, and Myanmar's Irrawaddy (Lang et al., 2019;Whitaker et al., 1974).By the early 1970s, they had vanished from about 80% of their historic range and experienced a staggering 95% decline in population size (Lang et al., 2019), leading to their recognition as a top priority EDGE (Evolutionary Distinct and Globally Endangered) reptile species (Gumbs et al., 2018).Presently, Gharials are limited to isolated areas in Nepal and India, with an estimated population of 300-900 breeding adults in the wild (Lang et al., 2019).
In Nepal, the once-thriving Gharial populations in the Koshi and Mahakali rivers became extinct, leaving only remnant populations in the Narayani-Rapti River system of Chitwan National Park (CNP) and the Karnali and Babai rivers within Bardia National Park (BNP) (Lang et al., 2019).The current population of the Rapti and Babai Rivers are comprised of species that were introduced from Breeding Center (GCBC) in Kasara, CNP in 1978(Ballouard et al., 2010).This effort was initiated during a period of rapid Gharial population decline when there were less than 81 individuals recorded, and hatchling survival rates in their natural habitat were less than 1% (DNPWC, 2018).
As per the latest census of 2016, the total population of Gharials in Nepal is estimated to be 198 individuals (Acharya et al., 2017).This marks a gradual increase in the population of Gharials since the previous surveys conducted in 2011 and 2013, which recorded 102 and 124 Gharials, respectively (Acharya et al., 2017).
The increase in the Gharial population has been credited to the implementation of ex-situ conservation measures and in-situ nesting success (Acharya et al., 2017;Bashyal et al., 2021).During an opportunistic visit in mid-June 2019, Bashyal et al. (2021) made a notable discovery at Dhanuse: three hatched nests with an estimated 100 hatchlings nearby.According to the most recent national survey, the Narayani River hosts a population of 84 Gharials, while the Rapti River is home to 82 individuals.In addition, the Babai River supports a population of 31 Gharials, while the Karnali River currently has only one individual remaining (DNPWC, 2018).
Despite the increase in the overall population, the number of adult Gharials has remained low and skewed towards females (Acharya et al., 2017).Due to their low population and the existing threats, it has become essential to conduct regular periodic monitoring of Gharials.The Government of Nepal has been systematically monitoring Gharials every 2 years since the initiation of the Gharial Conservation Action Plan in 2018, employing the direct count method (DNPWC, 2018).This enables conservation efforts to be implemented promptly to prevent any further decline in their population.
Further, it is important to recognize that infrastructure development is a significant cause of habitat disturbance and loss worldwide (Laurance et al., 2015); however, despite its far-reaching consequences, infrastructure development has not received as much attention as other drivers of biodiversity decline, such as agriculture (Pant et al., 2021;Simkin et al., 2022).The Government of Nepal (GoN) has initiated two large-scale projects in western Nepal.The Bheri Babai Diversion Multipurpose Project (BBDMP), which involves constructing an inter-basin water transfer system (Bheri River to the Babai River).This transfer of water, which amounts to 40 cubic meters per second, will enable irrigation of over 51,000 hectares of agricultural land in the Bardia and Banke districts.The project will also generate 46 MW of hydropower (GoN, 2023).Rani Jamara Kulariya Irrigation Scheme (RJKIS), which aims to develop canal system to improve irrigation water delivery to, and management, in the command area of 38,300 ha of Sudurpaschim Province, Nepal (GoN/DWRI, 2023).Both BBDMP and RJKIS have been labeled as "Projects of National Pride" by the GoN despite their adverse effects on aquatic ecosystems and their functions (Das, 2006;Lakra et al., 2011;Renöfält et al., 2010).These also include both upstream and downstream impacts on habitat for aquatic communities such as fish and reptiles, and water quantity and quality (Dudgeon, 2000;Pelicice et al., 2015).
The primary objective of this study is to meticulously establish a comprehensive baseline dataset encompassing critical information pertaining to the Gharial population status and the ecological characteristics of their habitats within the Karnali and Babai Rivers.This endeavor not only addresses the current knowledge gap but also lays a solid foundation for future research endeavors in the field of Gharial conservation and riverine ecosystem management, especially addressing potential impacts of large infrastructure development, such as irrigation projects.The insights from this study can aid in identifying areas vulnerable to infrastructure development and designing mitigation measures to minimize negative impacts on Gharials and other threatened species, such as freshwater dolphins, as well as their habitat.Ultimately, this study can contribute to shaping environmentally sustainable infrastructure plans that prioritize conservation objectives, facilitating more informed and responsible development decisions.
The Babai is a warm-water river system that enters BNP from its north-eastern boundary at Chepang and leaves at Babai Bridge (Parewa Odar), approximately 47 km downstream.The Karnali and Babai rivers, however, are not connected.Karnali is snow-fed and perennial, originating from the Tibetan Plateau that enters BNP at its north-western corner at the Chisapani Bridge and bifurcates into two major channels running downstream and south.We studied 30 km downstream that extends from Karnali bridge to Tallo Khaura Cross.The vegetation is sub-tropical, with floodplain vegetation on river sides, and Shorea robusta forest on surrounding uplands.BNP serves as the ecological lifeline for a diverse range of threatened species, including the Asian elephant (Elephas maximus), greater one-horned rhinoceros (Rhinoceros unicornis), tiger (Panthera tigris tigris), smooth-coated otter (Lutra perspicillata), and mugger crocodile (Crocodylus palustris) (Thapa et al., 2018).In this region, over 60% of the population is comprised of Indigenous Tharu communities residing along the rivers.While their primary sources of livelihood are agriculture and animal husbandry, they often face challenges in meeting their daily and incidental needs.
Consequently, they also rely on the rivers for resources, including fishing.This additional dependence on river resources can have negative impacts on Gharials and their habitat due to potential conflicts and resource depletion.

| Sampling design and data collection
This study was conducted during post-winter season in Nepal, between February and March 2023.This season is considered as favorable season for Gharial survey and monitoring because of low water level in the river channels compared to other seasons (Nishan et al., 2023).Further, the mild weather during this period allows Gharials to bask on the riverbank for a longer duration, thereby increasing the likelihood of sighting them (Adhikari et al., 2019).The data collection was done in two phases.

| Preliminary survey
In February 2023, preliminary surveys were conducted in two phases over a period of 6 days.During the first phase, first author conducted detailed informal interviews with five conservation officials of BNP and two personnel from the National Trust for Nature Conservation (NTNC), who had excellent knowledge of topography, ecological significance, human interactions, conservation efforts, and management considerations necessary for preserving these valuable natural resources of BNP, to understand the potential factors affecting the habitat use of Gharials such as presence of Invasive species and illegal fishing.These interviews assisted us in obtaining an overview of the current Gharial distribution within various segments of these rivers.In the second phase of our study, we used Google Earth Pro (Version: 7.3.4.8642) to digitize the mid-channel route of the river streams.We focused on the downstream midchannel route, covering a total distance of 47 km of the Babai River (Parewa Odar to Chepang) and 30 km of the Karnali River (Karnali bridge to Tallo Khaura Cross), where Gharials have been reported to inhabit.The extracted river channel routes were further divided into sections so that each sub-sectors can be surveyed in single day to minimize bias in habitat use measures due to the upstream and downstream movement of Gharial individuals within the river (Neupane et al., 2020;Poudyal et al., 2018) (Table 1).
Then, we monitored the Gharial population in the chosen river segments by wooden boat/walking downstream along the riverbanks.This monitoring was repeated three times to observe any differences in the number of Gharial individuals occupying the river segment.Once we confirmed that there were no significant differences in the number of observed individuals during three consecutive surveys (Figure S1), we proceeded with conducting a detailed habitat survey of the Gharial.During the preliminary field visit, an  , 2021).We observed that the data generated through virtual observation were close to the data generated through the GIS tool.Therefore, we decided to use the virtual estimation method for the major field survey.

| Major field survey and recording the habitat characteristics
During March of 2023, a detailed field survey was conducted, which involved generating sampling stations at 500 meter (m) intervals using ArcGIS 10.8.1.(ESRI, 2020).The rationale behind this strategy was that the habitat characteristics would differ significantly at 500 m intervals, and it was less likely that basking Gharials would move across consecutive stations during the survey.Several studies (Bhattarai et al., 2022;Lamichhane et al., 2022;Nishan et al., 2023) have utilized similar methods and assumptions for assessing habitat characteristics of various species of crocodiles, including the Gharials (Neupane et al., 2020).Habitat types were assessed on both sides of the river (left and right bank) at every survey station and at each Gharial sighting location.The survey was conducted during the mid-day between the hours 8:00-16:00 h as it is assumed that most Gharial individuals would be basking during that time (Adhikari et al., 2019).Further, breeding animals tend to bask in groups during late winter and early spring, and their courtship behavior provides a better opportunity to observe them (Choudhury & Rao, 1982).
Additionally, surveys were exclusively conducted during clear sunny weather conditions to mitigate potential biases in certain field measurements, such as water temperature, which were collected on different days.
A team of five individuals including two experienced local field guides aided with binoculars and a range finder moved downstream.
The team recorded the number of Gharials and their sex, and age class.
The sex was determined by observing distinct morphological feature called "Ghara" present in male adult only.Total length of the Gharial was estimated by calibrating with natural objects or landscape features using binoculars and was assigned into different age classes as hatchlings (<120 cm), juveniles (120-180 cm), sub-adults (180-270 cm) and adults (>270 cm) (Neupane et al., 2020;Rajbhandari & Acharya, 2015).
Based on preliminary field visit and published literature review (Bashyal et al., 2021;Neupane et al., 2020;Poudyal et al., 2018;Rajbhandari & Acharya, 2015), 10 pre-defined habitat characteristics were identified as potentially influencing Gharial occurrences; and those selected variables included: elevation, river width, mid-river depth, bank width, river substrate type, river bank slope, water temperature, presence/ absence of anthropogenic disturbance, invasive alien plant species and basking sites (see Table 2 for details).
The elevation at each sampling station and sighting location was recorded using a Garmin eTrex 30.River width and Bank width were estimated using the Vortex Optics Razor HD 4000 range finder.
The mid-river depth was estimated using a long-calibrated bamboo stick which was scaled at every 5 cm interval.River substrate type and presence/absence of [anthropogenic disturbance, invasive alien plant species (such as Lantana camara, Ageratum houstonianum, and Chromolaena odorata) and basking sites] were virtually estimated (Neupane et al., 2020;Nishan et al., 2023).To avoid the misidentification of the invasive alien species, primary author captured the high-resolution photograph of the recorded invasive species and conducted a taxonomic validation by cross-referencing it with the Book "Invasions of Alien Plant Species in Nepal" authored by Shrestha and Shrestha (2021).In our study, we rigorously assessed the presence or absence of Gharial basking sites by comparing riverbank width to the size of adult male Gharials.Sites with riverbanks equal to or larger than the dimensions of adult male Gharials were categorized as "basking site presence," while those with narrower riverbanks were labeled as "basking site absence."Further, the water temperature of the river was measured with temperature meter (Hanna HI991300).

| Data analysis
The present study utilized a binary logistic regression model to investigate the factors that influence the occurrence of Gharial at a selected site.The response variable was the binary presence or absence of Gharial at both sampling stations and sighting locations, whereas 10 pre-determined habitat variables (Table 2) were treated as independent (predictor) variables.The potential issue of multicollinearity was assessed before conducting regression analysis, and the variance inflation factor (VIF) was calculated using the "Faraway" R package (Boomsma, 2014).All independent variables were evaluated for VIF values less than 10, indicating no significant multicollinearity (Bowerman & O'connell, 1990).
Next, a global model logistic regression analysis was performed, with all independent variables included under the binomial family and logit link function.In order to rank the potential candidate models, we constructed multiple models and evaluated their performance using Akaike's Information Criterion with correction (AICc) for small sample sizes (Barton & Barton, 2020).The most promising candidate models with delta AIC ≤2 were selected and model averaging (covering 98% weightage (w)) was performed to obtain the final models, as suggested by Burnham and Anderson (2001) for minimizing uncertainty.
The "AICmodavg" package Mazerolle and Mazerolle (2017) was employed for this purpose.To assess the predictive performance of the best-fit model (models with ΔAIC ≤2), we calculated the Area Under Curve (AUC) of the Receiver Operator Characteristics (ROC) using the "ROCR" R package (Sing et al., 2005).The AUC values ranged from 1 to 0. The higher the AUC value, the better the model performance and vice-versa (Carter et al., 2016).According to Hosmer et al. (2000), receiver operator characteristic (ROC) values ranging from 0.7 to 0.8 are considered to exhibit acceptable discrimination, 0.8 to 0.9 are considered excellent, and values greater than 0.9 are considered superior.A p-value (Pr(>|z|)) less than .05 (at a significance level of 5%) was used to determine if the independent variables were significantly associated with the response variable in the model.Further, we utilized G-test (likelihood ratio test) (Agresti, 2012) to determine if there is a significant relationship between male and female-biased populations at 5% level of significance.Modeling analysis was conducted using the statistical package: R x 64 4.2.1 (R Core Team, 2020).

| Population abundance and distribution of Gharials in Bardia National Park
A total of 23 individuals of Gharials were recorded during the survey (Figure S1), with 22 individuals sighted in the Babai River and only one in the Karnali River.The total search effort was 1.93 km/h.
Among the recorded Gharials, 14 were adults (61%), six were subadults (26%), and three were juveniles (13%).Additionally, the study found that nine individuals could not be identified by sex, while five were identified as male and nine as female (Figure 3).Similarly, 83% of Gharials individuals were found basking, whereas the remainders were found swimming.The sex ratio was found to be 55.56 males per 100 females.However, result did not provide sufficient evidence to support the claim that the Gharial population in BNP has a biased sex ratio (χ 2 = 1.1429; p-value = .285).

| Influencing variables and probability of Gharial occurrence
Among the candidate models generated from binary logistic regression to elucidate the probability of Gharial occurrence, the model  4).Results of the model indicate that the probability of detecting Gharials in BNP increases with rise in mid-river depth, river width, and decreasing water temperature (Table 4 and Figure 5).

| DISCUSS ION
Given

| Population abundance of Gharials in Bardia National Park
Government records indicate that since the inception of the Gharial F I G U R E 4 ROC curve for optimal binary logistic regression model, displaying ACC (accuracy) and AUC (area under curve) metrics for performance evaluation.
River (35 individuals, 2.07%).These releases are an essential component of conservation efforts for the species, aimed at increasing their population and genetic diversity.However, the expected increase in the Gharial population has not been observed (Khadka et al., 2020;Neupane et al., 2020).The lower survival rate of the released individuals highlights the presence of multiple threats to the species, including monsoonal wash-offs into India, the presence of dams that hinder upstream passage from India to Nepal, natural mortality, and accidental deaths from entanglement in fishing nets (Acharya et al., 2017;Ballouard et al., 2010;Yadav et al., 2022).
The population status and habitat of Gharials in BNP have been the subject of numerous studies et al., 2019, 2021;DNPWC, 2018;Khadka et al., 2008;Thapa et al., 2018), which have provided valuable contributions to conservation efforts aimed at preserving this endangered species.Altogether a total of 23 Gharial individuals (five males, nine females, and nine unidentified) were observed in this study.Similarly, our findings are concurrent to the findings of Bashyal et al. (2021) where they recorded 19 individuals (18 Babai, one Karnali).However,  & Rao, 2017) in multiple sites in India and Nepal employed visual encounter surveys (Crump et al., 1994) while sailing random transects along riverine habitat to estimate the Gharial population size.
The direct count method does not account for imperfect detection (Barão-Nóbrega et al., 2022).Studies have indicated that imperfect detection is a prevalent issue in crocodilians, despite their substantial size (Balaguera-Reina et al., 2018;Barão-Nóbrega et al., 2022).Even Gharials, being among the largest species of crocodilians, may not be exempt from this problem during surveying.Therefore, we recommend combining modern methods such as Acoustic Monitoring, Drone Surveys, Remote Sensing and GIS and Citizen Science with traditional direct count methods to increase the accuracy of Gharial population count.Further, the majority of the Gharials were observed in the Babai River, with only a few individuals sighted in the Karnali river which is consistence with the findings of Bashyal et al. (2021).This may be due to lack of preferred sandbank habitat and high anthropogenic disturbance in Karnali river (Paudel et al., 2015).Further, the diversion of water from the Karnali River into irrigation canals for a commanded area The demographic composition of a species is a crucial element in population dynamics studies, as it can provide valuable insight into population health and reproductive success (di Sciara, 2018;Hoyle & Maunder, 2004).In our investigation, most Gharials observed were adults, while the number of sub-adults and juveniles was relatively small.These findings imply that recruitment levels within the population may not be sufficient to maintain a healthy population size.Moreover, the high proportion of Gharials with undetermined sex (39%), as observed in our study, raises significant concerns as it limits our ability to comprehend the population dynamics and the reproductive potential of the species.Similar findings were also reported by other studies that recorded a relatively higher number of unidentified sex Gharials (Bashyal et al., 2021;Khadka, 2011;Neupane et al., 2020;Poudyal et al., 2018).Use of genetic markers for sex identification is advisable yet multiple-year survey using current protocol is recommended.
A majority (83%) of the Gharials were found basking on sandy riverbanks, which is consistent with the findings of previous studies conducted by Neupane et al. (2020), Rajbhandari and Acharya (2013), Hussain (2009), andNair (2010).The underlying reason for the observed preference for sandy riverbanks among Gharials may be attributed to their greater ease in crawling on sandy surfaces, in contrast to rocky or clay ones (Neupane et al., 2020).
The sand found on riverbanks with lower moisture content acts as a moderator of extreme hot or cold environmental conditions, which ultimately lowers the likelihood of desiccation while basking in sun (Hussain, 2009;Neupane et al., 2020).Given the affinity of Gharials for sandy riverbanks as basking sites, it is advisable to prioritize the preservation and protection of such habitats in conservation initiatives.This may necessitate measures such as managing water flow and preventing soil erosion to ensure the continued availability of suitable basking areas.

| Influencing variables and probability of Gharial occurrence
The results of our binary logistic regression analysis indicate that three habitat predictors, namely mid-river depth, river width, and water temperature, were significantly associated with the likelihood of Gharial occurrence in BNP.However, the remaining variables did not demonstrate any significant predictive power in our analysis.
Our study's findings, which demonstrate a positive correlation between mid-river depth and the likelihood of Gharial occurrence, are consistent with earlier research (Hussain, 2009;Neupane et al., 2020;Rajbhandari & Acharya, 2015;Yadav et al., 2022) that highlight deep water areas as favorable habitats for Gharials.The availability of escape cover, such as deep water near basking sites, is a critical factor in their habitat selection since it allows them to retreat into water for safety when they feel threatened or disturbed by human activity (Hussain, 2009;Neupane et al., 2020).
Furthermore, female Gharials often choose nesting sites in areas near deeper water channels, which provide improved access to prey (Gupta et al., 2012), suggesting that the relationship between mid-river depth and Gharial occurrence may also be linked However, this hypothesis requires further investigation to establish its validity.
There exists a significant positive correlation between river width and Gharial occurrence, providing evidence that wider rivers may offer more suitable habitats for these reptiles.One possible explanation for this relationship is that wider rivers provide more space, diverse microhabitats, and better connectivity between different parts of the river ecosystem, which are important factors for Gharial survival and reproduction.It is important to note that our results differ from those of previous studies conducted by Neupane et al. (2020) andYadav et al. (2022), which found no significant impact of river width on the probability of Gharial occurrence.Further research is needed to reconcile these inconsistencies and gain a more comprehensive understanding of the ecological requirements and behaviors of Gharials.Further, due to the implementation of BBDMP, there will be a substantial increase in the flow of water in the Babai River (GoN, 2023), leading to alterations in the river's width.This change in flow may result in widening of the river in certain locations, while in other areas, it could cause erosion and lead to the narrowing of the river.This has both positive (more suitable habitat for feeding) and negative (reduce the availability of suitable nesting sites) implications for the conservation of the Gharial (Vashistha et al., 2021).Future studies should assess altered flow regime of Babai River due to the BBDMP and its impact on Gharial feeding and nesting ecology.
Our study revealed a significant negative correlation between water temperature and the probability of Gharial occurrence, which supports the findings of previous studies conducted by Shrestha (2000), and Adhikari et al. (2019).In winter, with air temperatures at around 18°C (64°F) and water temperatures near 17°C (63°F), Gharials basked throughout the day.Conversely, in summer, when air temperatures soared to about 35°C (95°F) and water temperatures reached approximately 25°C (77°F), gharials basked less and avoided midday sunbathing (Shrestha, 2000).Similarly, Adhikari et al. (2019) studied the effect of temperature fluctuation on the thermoregulatory behavior of captive Gharials and found that Gharials tended to remain in pools in the morning when the water temperature was warmer than the air temperature.However, more Gharials tended to move outside the water for basking when the water temperature dropped relative to the air temperature.These findings hold profound implications for the comprehensive conservation and sustainable management of Gharial populations in the wild.To bolster these initiatives, it is imperative to prioritize the development of environmentally conscious infrastructure that not only addresses economic demands but also harmonizes with the overarching objectives of conservation.This mandates a holistic strategy encompassing the preservation of critical habitats and meticulous consideration of the myriad factors influencing Gharial ecology, encompassing aspects like basking patterns, feeding behaviors, and nesting preferences.By highlighting the link between water temperature and Gharial behavior, our study also suggests that Gharials are sensitive to environmental conditions, which could make them particularly vulnerable to the effects of climate change.

| CON CLUS ION
This study contributes to the selection of interventions aimed at protecting Gharials and their associated habitats, specifically sandy banks in BNP.The Gharial population is comparable to previous studies and has remained stable over the years.This study establishes a population baseline for Gharials within the river system before the construction of large infrastructure projects, such as dams and irrigation canals.Continuous monitoring of Gharial populations along the river system after water release and/or diversion is recommended.This will help to evaluate the impact of large infrastructure projects, specifically the BBDMP and RJKIS, on the population and their associated habitat characteristics.Both projects are scheduled to be completed around 2025 or beyond.
The results of the best fit model indicate that the likelihood of detecting Gharials in BNP increases with higher mid-river depth and river width, as well as lower water temperature.These factors are crucial for developing sustainable infrastructure plans that ensure species conservation and maintain the ecological integrity of the river ecosystem.Considering the high demand for large infrastructure construction, such as dams in the river channel for electricity and irrigation, the assessment of environmental flow is crucial and likely to benefit Gharials.
Within BNP, riverbanks are protected yet habitat outside the core areas suffers from ongoing habitat degradation such as gravel extraction and mining.Use of citizen scientists (Bonney et al., 2014) to collect unbiased data to monitor both the population and habitat is recommended.There are analytical tools for robust analysis of unbiased data collected from the citizen scientist.Shreejan Gautam: Investigation (equal); resources (equal); writing -review and editing (equal).Santosh Ayer: Investigation (equal); and is listed as critically endangered on the Red List of the International Union for Conservation of Nature (Lang et al., 2019) and appended in Appendix I of the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES, 2021).Similarly, the Government of Nepal has listed in Schedule-I of the protected species under the National Parks and Wildlife Conservation Act regression, mid-river depth, occupancy, river width, sex biased, water temperature T A X O N O M Y C L A S S I F I C A T I O N Biodiversity ecology, Community ecology, Ecosystem ecology, Population ecology, Zoology F I G U R E 1 Gharial observed basking in the riverbank of Babai River stream.the Narayani River in 1985 and 1990, respectively.Recognizing the critical need for conservation, the government of Nepal, in collaboration with the Smithsonian Institute and the Frankfurt Zoological Society, established the Gharial Conservation and

F
Distribution map of gharial localities inside Bardia National Park is highlighted with light green circles.The sampling stations in Babai River are highlighted with yellow circles, and sampling stations in Karnali River are highlighted with orange circles.experienced local field guide from BNP, who had prior experience in studying crocodilian species, estimated the slopes of 40 riverbanks.The riverbanks were chosen to represent both sides of the river and varied topographic conditions.Simultaneously, a second experienced local field guide recorded the GPS coordinates of the visually estimated riverbank slopes, considering the absence of Gharials and other megafaunal species to avoid disturbance to species and any potential risks.Next, we validated the numerical values of the virtual slope estimation method using ArcGIS 10.8.1 (ESRI, 2020) by estimating the slope of each recorded point using DEM (Digital Elevation Model, 12.5 m resolution) (ASF Name and length (km)  of section of Babi and Karnali river of Bardia National Park.TA B L E 2Habitat variables used in the logistic regression analysis with their type, description, and reference.and sighting locations are characterized by the presence or absence of Gharials, where "1" denotes the presence and "0" denotes the absence of Gharials.Elevation [Elev]ContinuousElevation at mean sea level (m) at each sampling station and sighting location.Neupane et al. (2020) Mid-river depth [MRD] Continuous Mid-river depth (m) at each sampling station and sighting location.Neupane et al. (2020) River width [RW] Continuous River width (m) at each sampling station and sighting location.Neupane et al. (2020) Bank width [BW] Continuous Bank width (m) at each sampling station and sighting location.Neupane et al. (2020) Riverbank substrate type [RST] Categorical Sandy bank, sandy gravel, sandy grassy, sandy, grassy & gravel, gravel, gravel grassy and rocky Neupane et al. (2020) Riverbank slope [RBS] Categorical Gentle (<15° slope) Moderate (15-25° slope), Moderately steep (25-35° slope), Steep (35-55° slope) & very steep (>55° slope) He et al. (2020); Nishan et al. (2023); Thong et al. (2018) Water temperature [WT] Continuous Water temperature (°C) at each sampling station and sighting location.Adhikari et al. (2019); Hussain (1999) Presence/absence of basking sites [P.A.BS] Categorical Present = "1," absent = "0" of Gharial basking site at sampling station and sighting location.Bashyal et al. (2021) Presence/absence of anthropogenic disturbance [P.A.AD] Disturbance variable Categorical Evidence of any anthropogenic disturbance (washing clothes, cattle presence, grass collection, elephant rides, human walking trails, fishing, non-timber forest products (NTFPs) collection, swimming, gabion wall dikes) present = "1," absent = "0" at each sampling station and sighting location.Nishan et al. (2023) Presence/absence of alien invasive plant species [P.A.Inv]CategoricalEvidence of any alien invasive plant species present = "1," absent = "0" at each sampling station and sighting location.Nishan et al. (2023)

E 3
Bar graph showing the distribution of individuals based on age class (subplot a) and sex category (subplot b).Subplot "a" shows the number of juveniles, sub-adults, and adults, while subplot "b" shows the number of males, females, and unidentified individuals.withadditive effect comprising mid-river depth, availability of basking site, river width, and water temperature exhibited the lowest Akaike information criterion (AIC) value of 41.70 and the highest weight of 0.15.This model is considered the most parsimonious and provides the best explanation of the probability of Gharial occurrence (Table 3).The dominant model exhibited a remarkable level of proficiency, with an estimated area under the Receiver Operating Characteristic (ROC) curve of 0.99 and an accuracy score of 0.98 (98%) (Figure 4).Out of the 10 habitat predictors investigated in this study, three habitat predictors, mid-river depth (β = 3.958e+00, SE = 1.370e+00, p = .00386),river width (β = 5.598e-02, SE = 2.481e-02, p = .02404),and water temperature (β = −8.085e-01,SE = 3.368e-01, p = .01638)were significantly correlated with the probability of Gharial occurrence (Table the global decline of Gharial populations caused by anthropogenic disturbance, unique habitat requirements, low hatchling survival rates, and sex-biased population densities(DNPWC, 2018), urgent large-scale conservation efforts (ecofriendly infrastructure development, implementing captive breeding programs, promoting sustainable fisheries management, launching species translocation programs, and collaborative crossborder monitoring initiative) are necessary for their sustainable conservation.This study identifies the primary distribution locations of Gharials, assesses their current population status, and investigates influential variables affecting their distributions in two major river streams of western Nepal.The findings of this study are profoundly influential, as they provide the essential insights needed to shape infrastructure plans that prioritize environmental conservation, thus promoting robust and informed decision-making in this critical sphere.

F I G U R E 5
Figure "a," "b," and "c" showing the relation of occurrence of detected Gharial with mid-river depth, water temperature, and river width, respectively.
of 38,300 hectares (GoN/DWRI, 2023) is expected to result in a significant reduction in river discharge, which may have adverse effects on the population viability of Gharials inhabiting the Karnali River.Therefore, this anthropogenic water extraction for irrigation purposes requires careful consideration to ensure the conservation of the species and the maintenance of the ecological integrity of the river ecosystem.In addition,Vashistha et al. (2021) observed 72 individuals from the Katerniaghat Wildlife Sanctuary (KWS), which is located downstream of the Karnali river.It is possible that KWS Gharial population might have been formed by individuals released in Nepal who moved downstream (particularly hatchlings and juveniles) during periods of high-water levels or flood, but movement data to support this supposition is lacking.The Girijapuri barrage, located within the KWS is equipped with control gates that allow water to flow downstream but restrict the movement of Gharials upstream, once they cross the Girijapuri barrage gates and enter the Ghaghara River(Vashistha et al., 2021).Furthermore, evidence of breeding population of Gharial in Gaire khola -contributing river tributaries to Mohana River (5 km aerial stretch from western channel of Karnali River) -shows potential for possible movement of reptiles especially in high floods and their breeding success in an areas due low level of threats (especially illegal fishing and sand and gravel extraction) and availability of good nesting habitat.
to reproductive success.Significant population recovery, with evidence of breeding success, has been recorded through the strategic selection of sites with deep-water channels along the main tributaries of the Karnali River-Geruwa, flowing through protected areas, in the recent past (Personal communication: Bishnu Shrestha, Former Chief Warden-BNP).With the exception of the Babai River, the environmental flow of water discharge in the Karnali River is worth considering, given the changes in the level | 11 of 14 DHAMI et al.of water flow in the channel due to the construction or diversion of water in RJKIS, particularly during the dry season.Further, it is plausible that Gharials' preference for deep water segments is linked to their thermoregulation needs.Deeper water segments maintain warmer temperatures below the surface during winter and cooler temperatures in summer, providing a temperature gradient that Gharials can utilize for efficient temperature regulation.