Tiger - Prey Report 2019
Tiger - Prey Report 2019
Tiger - Prey Report 2019
S TAT U S O F T I G E R S A N D P R E Y I N N E PA L 2 0 1 8
GOVERNMENT OF NEPAL
Ministry of Forests and Environment
Department of National Parks and Wildlife Conservation
and
Department of Forests and Soil Conservation
Kathmandu, Nepal
2018
TECHNICAL TASK FORCE
Laxman Prasad Poudyal, Ecologist, Department of National Parks and Wildlife Conservation
Bhupendra Prasad Yadav, Assistant Ecologist, Department of National Parks and Wildlife Conservation
Rishi Ranabhat, Assistant Ecologist, Department of National Parks and Wildlife Conservation
Sujan Maharjan, Assistant Planning Officer, Department of Forests and Soil Conservation
Sabita Malla, Wildlife Biologist-Senior Manager, WWF Nepal
Babu Ram Lamichhane, Wildlife Research Officer, National Trust for Nature Conservation
Samundra Ambuhang Subba, Research Officer, WWF Nepal
Saroj Koirala, GIS Officer, WWF Nepal
Sheren Shrestha, Senior Research Officer, WWF Nepal
Ashish Gurung, Conservation Officer, National Trust for Nature Conservation
Umesh Paudel, Natural Resource Conservation Assistant, National Trust for Nature Conservation
Tek Raj Bhatt, Senior Program Officer, ZSL-Nepal
Sailendra Raj Giri, Field Biologist, ZSL-Nepal
TECHNICAL TEAM
Gopal Prakash Bhattarai, Deputy Director General, Department of National Parks and Wildlife Conservation
Madhuri Karki (Thapa), Under Secretary, Department of Forests & Soil Conservation
Dhananjaya Lamichhane, Under Secretary, Ministry of Forests and Environment
Rajendra Dhungana, Assistant Planning Officer, Ministry of Forests and Environment
Hari Bhadra Acharya, Chief Conservation Officer, Parsa National Park
Ram Chandra Kandel, Chief Conservation Officer, Chitwan National Park
Dil Bahadur Purja Pun, Chief Conservation Officer, Banke National Park
Ramesh Kumar Thapa, Chief Conservation Officer, Bardia National Park
Bed Kumar Dhakal, Chief Conservation Officer, Shuklaphanta National Park
Shant Raj Jnawali, PhD, Chief of Party, Hariyo Ban Program, WWF Nepal
Shiv Raj Bhatta, Director of Programs, WWF Nepal
Naresh Subedi, PhD, Conservation Program Manager, National Trust for Nature Conservation
Kanchan Thapa, PhD, Technical Advisor, Hariyo Ban Program, WWF Nepal
Bhagawan Raj Dahal, PhD, Program Manager, ZSL-Nepal
Chiranjibi Prasad Pokheral, PhD, Program Manager, National Trust for Nature Conservation
ADVISORS
Man Bahadur Khadka, Director General, Department of National Parks and Wildlife Conservation
Krishna Prasad Acharya, Director General, Department of Forests and Soil Conservation
Maheshwar Dhakal, PhD, Chief, Biodiversity and Environment Division, Ministry of Forests and Environment
Govinda Gajurel, Member Secretary, National Trust for Nature Conservation
Ghana Shyam Gurung, PhD, Country Representative, WWF Nepal
Hem Sagar Baral, PhD, Country Representative, ZSL-Nepal
CITATION:
DNPWC and DFSC. (2018). Status of Tigers and Prey in Nepal. Department of National Parks and Wildlife Conservation
& Department of Forests and Soil Conservation. Ministry of Forests and Environment, Kathmandu, Nepal.
STATUS OF TIGERS AND PREY IN NEPAL 2018 iii
STATUS OF TIGERS AND PREY IN NEPAL 2018 v
ACKNOWLEDGEMENT
Nepal, among the 13 range states, is committed to Our special thanks go to all chief conservation officers
doubling tiger by 2022 as per St. Petersburg Declaration, and their staff of Parsa, Chitwan, Banke, Bardia, and
2010. Tiger is mainly found in the five protected areas Shuklaphanta National Parks; District Forest Officers of
and area outside the protected areas along the low land the (Dadeldhura, Kanchanpur, Kailali, Bardia, Surkhet,
area of Nepal. In line with Nepal Tiger Recovery Program Banke, Salyan, Dang, Arghakhanchi, Kapilvastu,
(NTRP) and Tiger Conservation Action Plan (2016-2021), Rupendehi, Palpa, Nawalparasi, Chitwan, Parsa,
third nationwide tiger and prey survey 2018 has been Makwanpur, Bara and Rautahat and their staff of tiger
conducted to update the national database on tiger and bearing forest districts, staff of NTNC, ZSL Nepal, WWF
their prey in Nepal. Nepal and TAL field offices for their participation and
support in the field survey implementations. Ms. Sabita
Third Nationwide Tiger and Prey Survey 2018 Nepal has Malla, Mr. Samundra Subba, Mr. Saroj Koirala, Mr.
shown 19% increase from 2013 estimate of 198 tigers and Shailendra Giri, Mr. Babu Ram Lamichhane deserves
results give us pride and encouragement to continue our special thanks for their outstanding contribution in
tireless efforts to save the magnificent animal which is designing the survey, field technical support, data
endangered globally. The survey was conducted by the analysis, report writing and publication. We would also
Government of Nepal in collaboration with National like to thank Dr. Rajan Amin, Dr. Marcella J. Kelly and Ms.
Trust for Nature Conservation, Zoological Society of Rebecca May for their thorough review and constructive
London-Nepal office, and World Wildlife Fund-Nepal comments in finalizing this report. Besides, Nepali Army,
office. We would like to thank Dr. Maheshwar Dhakal - law enforcement agencies, buffer zone user committees,
Joint Secretary, MoFE; Mr. Govinda Gajurel - Member community forest user groups, CBAPUs, citizen scientist,
Secretary, NTNC; Dr. Hem Sagar Baral - Country students and volunteers deserve special thanks for their
Representative, ZSL Nepal and Dr. Ghana Shyam Gurung role and participation in the survey.
- Country Representative, WWF Nepal for their support
and advice in making this national survey a success. The technical and financial support of WWF Nepal, ZSL
Nepal, NTNC, USAID-Hariyo Ban Program, Leonardo
We extend our thanks to the technical committee DiCaprio Foundation, KfW/IUCN, Panthera, WildCats
members: Mr. Gopal Prakash Bhattarai (DNPWC, Chair), Conservation Alliance for nationwide tiger survey
Mr. Dhananjaya Lamichhane (MoFE), Ms. Madhuri deserves special acknowledgement.
Karki (Thapa) (DFSC), Mr. Laxman Prasad Poudyal
(DNPWC), Dr. Shant Raj Jnawali (WWF Nepal), Mr. Shiv We hope this technical report will be useful to all policy
Raj Bhatta (WWF Nepal), Mr. Bhupendra Prasad Yadav makers, protected areas and divisional forest managers,
(DNPWC), Mr. Rishi Ranabhat (DNPWC), Mr. Sujan conservationist, academia, and general readers nationally
Maharjan (DFSC), Dr. Kanchan Thapa (WWF Nepal), Dr. and internationally. Finally, we reiterate our sincere
Naresh Subedi (NTNC), Dr. Chiranjibi Prasad Pokheral thanks to every individual and institutions who made
(NTNC), and Dr. Bhagawan Raj Dahal (ZSL Nepal) for the nationwide tiger and prey survey a huge success and
their untiring work throughout the survey. timely publication of this technical report.
EXECUTIVE SUMMARY
Tiger is an apex predator and an umbrella species that a total sampling effort of 2,485 km, provided the prey
ensures well-being of entire ecosystems across its habitat density estimates. An extensive effort of 53,843 person-
range in Asia. Its population and distribution range had days and 1,735 elephant-days was invested to complete
drastically declined in the last century, with the species the field work for the nationwide survey.
verging towards extinction. Over the past few decades,
however, implementation of strategic conservation Data analysis was done with established latest software.
interventions has aided in the revival of the tiger. Occupancy estimate was derived using PRESENCE,
tiger abundance and density was analyzed using ‘secr’
In 2010, Nepal along with other range countries and SPACECAP package in R environment, respectively.
endorsed the St. Petersburg declaration to double the DISTANCE software was used to estimate prey density.
tiger population by 2022. Since then, Nepal has been
conducting four-year periodic assessments to track the Tiger signs were detected in 12 districts (Bara, Parsa,
progress towards reaching the national target of 250 Makwanpur, Chitwan, Nawalparasi, Dang, Salyan,
tigers. The first and the second nationwide assessments Banke, Bardia, Surkhet, Kailali and Kanchanpur) out of
carried out in 2009 and 2013 estimated 121 and 198 18 districts surveyed across TAL-Nepal. Altogether, 494
tigers, respectively. unique tiger signs were detected, providing an occupancy
estimate of approximately 68% (11,057 km2 of the total
This report synthesizes the findings of the third 16,261 km2) across the landscape. Segregating further,
nationwide tiger and prey survey, led by the Department habitat occupied by tigers in PAs was found to be as high
of National Parks & Wildlife Conservation (DNPWC) and as 98% (6,828 km2), as against 60% (5,576 km2) outside
Department of Forests and Soil Conservation (DFSC) Protected Areas.
with the support from WWF Nepal, National Trust for
Nature Conservation (NTNC), and ZSL-Nepal. Citizen Tigers were captured in camera traps in 482 grids, or 29%
scientists and students of various institutions provided of the total 1,643 grids. Altogether, 4,388 photographs
for human resource needs to achieve this landscape-level of 209 individual tigers [Parsa National Park (PNP) -
exercise. 15, Chitwan National Park (CNP) - 85, Banke National
Park (BaNP) – 17, Bardia National Park (BNP) - 77 &
The survey was conducted in all potential tiger habitats Shuklaphanta National Park (ShNP) – 15] were obtained.
in the Terai Arc Landscape (TAL) of Nepal, following Independent detections (1,136) of identified individuals
nationally approved Tiger and Prey Base Monitoring were analyzed to estimate protected area-wise tiger
Protocol, 2017, between November 2017 to April 2018. populations. This revealed an estimate of 18 (16-24) tigers
The major objectives were to estimate: i) habitat occupied in Parsa National Park and adjoining forests, 93 (89-102)
by tigers, ii) tiger abundance and density, and iii) prey tigers in Chitwan National Park and adjoining forests,
density. This report also recommends appropriate 21 (18-30) tigers in Banke National Park and adjoining
measures to address challenges in tiger conservation and forests, 87 (82-97) tigers in Bardia National Park and
presents the road map towards achieving TX2 by 2022. adjoining forests, and 16 (15-21) tigers in Shuklaphanta
National Park and adjoining forests. This accounts for the
To implement this survey, various working committees total national estimate of 235 tigers. A naïve comparison
were formed at central to field levels. Orientation indicates an approximate increase in the national tiger
trainings were provided to field survey teams prior to population by 19%, within the four-year period. Notable
the field work. Habitat occupancy survey was carried increase was recorded in respective populations of BNP,
out in 112 (15 km×15 km) grid cells covering 16,261 km2 BaNP and PNP and their adjoining forests; population
of forested habitats to estimate tiger distribution across in ShNP remained stable, while a marginal decline was
TAL-Nepal. Camera-trap survey was carried out in 1,643 recorded in CNP.
(2 km×2 km) grid cells encompassing a) sampling area:
6,572 km2, b) effective sampling area (total sampling Tiger density (per 100 km2) in PAs and adjoining forests,
area plus buffer): 12,356.6 km2 (protected areas, buffer was estimated to be 0.92 (SD 0.15), 3.28 (SD 0.19),
zones, corridors, and adjoining forests) to estimate tiger 0.97 (SD 0.12), 4.74 (SD 0.28) and 0.96 (SD 0.14) in
population and density. Overall effective sampling effort PNP, CNP, BaNP, BNP, and ShNP, respectively. Habitat
was 27,829 trap days. Survey of 1,294 line transects, with occupancy and usage outside PAs, documented by the
viii STATUS OF TIGERS AND PREY IN NEPAL 2018
present survey, is well below the optimal potential; this better connectivity as well as greater support towards
provides opportunities for increase through management conservation by communities.
interventions that reduce human disturbances and
improve habitat quality outside PAs. Nepal has come a long way in its journey towards doubling
tiger numbers by 2022. Yet, strategic interventions over
Wild prey species detected during line transect survey the next four years will be critical to achieve this goal. This
included four deer species (spotted deer, sambar, hog report compiles necessary efforts needed at both national
deer, barking deer), two antelopes (blue bull and four- and site levels. These include policy initiatives, research
horned antelope), wild boar, gaur, and two primate priorities, further improvements in protection and
species (rhesus macaque and langur). Prey density (per management interventions as well as greater engagement
km2) in Protected Areas and adjoining forests were 22 with communities. Improving on current trends through
(SE 3.8), 70.7 (SE 7.5), 8.1 (SE 1.6), 77.5 (SE 6.6) and 68 these interventions and incorporation of new emerging
(SE 7) in PNP, CNP, BaNP, BNP, and ShNP, respectively. understanding, with sustained political commitment
The overall positive trends in habitat occupancy, tiger by the Government and consistent efforts of diverse
abundance in Nepal plausibly relates to positive outcomes stakeholders, Nepal may well become the first country to
of improved protection and management measures, achieve its commitment to global tiger conservation.
STATUS OF TIGERS AND PREY IN NEPAL 2018 ix
TABLE OF CONTENTS
Foreword....................................................................................................... iii
Acknowledgement............................................................................................v
Executive Summary.......................................................................................vii
1. Introduction.............................................................................................. 1
2. Objectives..................................................................................................3
3. Study Area.................................................................................................5
4.Materials and Methods..............................................................................7
4.1 Survey Organization............................................................................................... 7
4.2 Field Methods.........................................................................................................7
4.2.1 Field Training.............................................................................................. 7
4.2.2 Survey Time Frame and Human Resources............................................... 7
4.2.3 Tiger Habitat Occupancy Surveys............................................................... 7
4.2.4 Camera Trap Based Mark-Recapture of Tigers..........................................8
4.2.5. Line Transect Surveys for Prey Base Density Estimation..........................8
5. Data Analysis........................................................................................... 11
5.1 Tiger Habitat Occupancy Modeling......................................................................11
5.2 Tiger Population Abundance Estimation..............................................................11
5.3 Tiger Density Estimation.......................................................................................11
5.4 Tiger Prey Density Estimation............................................................................. 12
6 Results..................................................................................................... 13
6.1 Tiger Habitat Occupancy...................................................................................... 13
6.1.1 Sampling Effort and Sign Detection......................................................... 13
6.1.2 Tiger Habitat Occupancy and Detection Probability............................... 13
6.2 Tiger Abundance................................................................................................... 13
6.2.1 Sampling Effort and Tiger Captures........................................................ 13
6.2.2 Tiger Population Abundance Estimates.................................................. 15
6.3 Tiger Density Estimates........................................................................................ 16
6.4 Prey Density Estimates......................................................................................... 16
6.5 Habitat Use of Tigers Outside PAs....................................................................... 18
7.Discussion............................................................................................... 19
7.1 Tiger Distribution, Habitat Usage and Occupancy.............................................. 19
7.2 Tiger Abundance and Density ............................................................................. 19
7.2.1 Methods Used and the Extent of Areas Covered .................................... 19
7.2.2 Trend in Minimum Population Based on Individual Tiger Captures.....20
7.2.3 Tiger Abundance Estimates.....................................................................20
7.2.4 Tiger Density............................................................................................ 21
7.2.5 Factors Governing Tiger Population Abundance and Density................ 21
7.3 Prey Densities.......................................................................................................23
8. Management Implications and the Road Map Towards TX2....................25
9. References............................................................................................... 27
10. Annexure.................................................................................................29
x STATUS OF TIGERS AND PREY IN NEPAL 2018
1 . INTRODUCTION
The tiger (Panthera tigris) is an icon for wildlife Monitoring Protocol (DNPWC, 2017). Two nationwide
conservation. Presently, there are estimated to be around tiger assessments were carried out in 2009 and 2013,
3,900 tigers in the wild (WWF, 2016) and the tiger’s range respectively. The first assessment estimated 121 tigers in
has reduced by 95% since historical times (Dinerstein Nepal (Karki et al., 2009), and was crucial in bringing
et al 2007). Realizing this, the 13 tiger range countries major conservation policy changes in the country. The
led by Global Tiger Initiative (GTI) and supported by survey also established the tiger population baseline
global community joined hands to reverse the decline for the government’s commitment to double the tiger
of this iconic species (GTRP, 2010). This culminated in population from 121 to 250 tigers by 2022. The second
a global commitment to double the tiger population by assessment recorded a 63% increase in the country’s tiger
2022 (TX2) and the adoption of Global Tiger Recovery population from the 2009 baseline, with an estimated
Program (GTRP) in 2010 (GTRP, 2010). population of 198 tigers (Dhakal et al., 2014). It also
provided better insight of tigers along the transboundary
Tiger conservation in Nepal began with the launch of Terai Arc Landscape with empirical evidence of tiger
the tiger ecology project in 1972 in Chitwan (McDougal, movement across the borders (Chanchani et al., 2014).
1977; Smith, 1993) followed by the establishment of The study identified several site-specific management
Chitwan National Park (CNP), the first national park of and conservation gaps and recommended appropriate
Nepal. Since then, Nepal Government has established measures to address them.
an additional four protected areas (PAs); Parsa National
Park (PNP), Banke National Park (BaNP), Bardia National Accordingly, Nepal has invested intensive efforts over the
Park (BNP) and Shuklaphanta National Park (ShNP) for last four years to strive towards Tx2. Policy documents
the conservation of tigers. Nepal also gradually shifted such as National Biodiversity Strategy and Action Plan
its conservation focus from protecting isolated core PAs (2014-2020), Terai Arc Landscape Strategy and Action
to designing and managing conservation landscapes Plan (2015-2025), Forest Policy (2015), and President
focusing on providing connectivity between wildlife Chure Terai Madhesh Conservation and Management
populations (Wikramanayake et al., 1998). Terai Arc Master Plan (2017) and Tiger Conservation Action Plan
Landscape (TAL) was the first conservation landscape in (2016-2020) were developed and endorsed which are
Nepal that was designed based on tiger dispersal model the major guiding documents for tiger conservation in
with the larger goal of mainstreaming species and forest Nepal.
restoration into the rural development agenda. Almost
two decades of landscape level conservation coupled The third 2018 nationwide tiger and prey status
with law enforcement efforts has started to pay off with assessment was carried out by the Government of Nepal
increasing tiger numbers and extent illustrating tiger (GoN) led by Department of National Parks and Wildlife
dispersal is conceivable and breeding habitat can be Conservation (DNPWC) and Department of Forests and
restored and at the same time enhancing local livelihoods Soil Conservation (DFSC) in partnership with WWF-
(Chanchani et al., 2014; Thapa et al., 2018). Nepal, ZSL-Nepal and the National Trust for Nature
Conservation (NTNC). This report presents the findings
The Nepal Government conducts country-wide of the assessment along with appropriate measures to
assessment of the status of tiger and prey every four years, address identified challenges in tiger conservation and a
following the nationally approved Tiger and Prey Base road map for doubling tiger numbers by 2022.
2 STATUS OF TIGERS AND PREY IN NEPAL 2018
© DNPWC/WWF Nepal
STATUS OF TIGERS AND PREY IN NEPAL 2018 3
2 . OBJECTIVES
THE MAIN OBJECTIVES OF THE NATIONWIDE TIGER AND PREY SURVEY WERE:
© WWF Nepal
STATUS OF TIGERS AND PREY IN NEPAL 2018
Figure 1.Location of TAL – Nepal showing tiger-bearing PAs and corridors with camera trap locations (black dots) and habitat occupancy grid cells (15 km x15 km).
STATUS OF TIGERS AND PREY IN NEPAL 2018 5
3 . STUDY AREA
Tigers in Nepal are distributed across the Terai and species of herpeto-fauna and more than 125 species of
Churia habitats within TAL. The National Tiger and Prey fish (MoFSC, 2015). Other high profile threatened species
Survey - 2018 was conducted across TAL, Nepal (Figure include greater one-horned rhinoceros (Rhinoceros
1). The TAL is a global priority conservation landscape unicornis), swamp deer (Rucervus duvaucelii), Asian
for tigers that extends from the Bagmati river, Nepal in elephant (Elephas maximus), Gangetic dolphin
the east to the Yamuna river in Uttarakhand, India in the (Platanista gangetica), Bengal florican (Houbaropsis
west, with an area of 51,002 km2 (Wikramanayake et al., bengalensis), Gyps vulture (Gyps spp) and gharial
1998). TAL-Nepal is spread over 24,710 km2, covering (Gavialis gangeticus). This landscape is a mosaic of
18 districts; Dadeldhura, Kanchanpur, Kailali, Bardia, early successional tall grasslands established in the
Salyan, Surkhet, Banke, Dang, Arghakhanchi, Kapilvastu, alluvial floodplains to climax stage Sal forests at lower
Rupendehi, Palpa, Nawalparasi, Chitwan, Makwanpur, elevations, and broad-leaved forests in the Churia range.
Bara, Parsa and Rautahat (MoFSC, 2015). PNP, CNP, Major habitat types include Sal forests, riverine forests,
BaNP, BNP and ShNP are prime tiger habitats in TAL- mixed hardwood forests and grasslands (MoFSC, 2015).
Nepal; biological corridors (Someshwor, Barandabhar, TAL represents sub-tropical monsoonal climate with
Kamdi, Khata, Karnali, Basanta, Laljhadi, Brahmadev three distinct seasons: cool-dry (November-February),
and Jogbuda) provide habitat connectivity among these hot-dry (March-June) and monsoon (July-October).
PAs of Nepal and with the transboundary PAs in India The average temperature in the cool season drops to
(Chanchani et al., 2014). 5°C in January and rises to 40°C in the hot dry season
(MoFSC, 2015).
The highly productive alluvial grasslands and riverine
forests of TAL are the major habitats of tigers; these also The tiger protected areas within Nepal TAL are listed in
supports 85 species of mammals, 565 species of birds, 47 Table 1.
(** denotes the PAs with revised status of PA or changes in their former size). WR = Wildlife Reserve; NP = National Park.
6 STATUS OF TIGERS AND PREY IN NEPAL 2018
© WWF Nepal
STATUS OF TIGERS AND PREY IN NEPAL 2018 7
4.1 SURVEY ORGANIZATION 4.2.2 SURVEY TIME FRAME AND HUMAN RESOURCES
The field survey was carried out from December 2017
At the national level, an Advisory Committee was setup to April 2018 (Annex-9). It was initiated from PNP
under the chairmanship of the Director General, DNPWC following formal inauguration by Dr. Yubak Dhoj GC,
to provide overall guidance for the survey. Members of Secretary of MoFE. An extensive effort of 53,843 person-
this committee included Director General-DFSC, Member days and 1,735 elephant-days was invested to complete
Secretary - NTNC, Country Representative WWF-Nepal, the nationwide survey (Annex-9).
and Country Representative - ZSL-Nepal. A Technical
Committee chaired by Deputy Director General, DNPWC 4.2.3 TIGER HABITAT OCCUPANCY SURVEYS
was formed at the central level with representatives from Standardized method was followed for tiger occupancy
DNPWC, DFSC, WWF Nepal, ZSL-Nepal and NTNC for survey (DNPWC, 2012; DNPWC, 2017; Barber-Meyer et
overall coordination and supervision of the survey. A al., 2012). The occupancy survey was carried out in 112
Technical Task Force chaired by the Ecologist, DNPWC, grid cells (each measuring 15 km×15 km) that spanned
and comprising representatives of DNPWC and DFSC, across 16,261 km2 of potential tiger habitat (forested and
and wildlife biologists from WWF-Nepal, ZSL-Nepal and grassland) in TAL (Figure 1). Each grid cell was divided
NTNC, designed the survey, provided technical training into 16 sub-grid cells (3.75 km x 3.75 km). One sub-
and guidance for field work, analyzed data and produced grid cell was randomly selected to include an element
the report. of randomness in spatial distribution of survey routes
(Karanth et al., 2008; Barber-Meyer et al., 2012). The
Field Implementation Committees were formed at PA number of spatial replicates (i.e. km walked) per grid
level in PNP, CNP, BaNP, BNP, and ShNP under the cell was proportional to the percentage of tiger habitat
chairmanship of respective Chief Conservation Officers. (Karanth et al., 2008; Barber-Meyer et al., 2012). For
Members of these committees included respective District grid cells with 100% tiger habitat, 40 km was sampled
Forest Officers, Officers-in-Charge of NTNC field offices, by traversing random grid in every survey route. Grid
Managers of TAL-PABZ/CBRP, Field Officers from ZSL- cells with less than 10% habitat cover were discarded.
Nepal and other relevant stakeholders. The details of the Each contiguous 1 km segment was considered a ‘spatial
personnel involved is provided in Annex-10. replicate’ (Hines et al., 2010; Barber-Meyer et al., 2012).
Each replicate comprised 10 segments of 100 m each and
4.2 FIELD METHODS the data was recorded at every 100 m avoiding spatial
auto-correlation by accounting single records for each
4.2.1 FIELD TRAINING unique species per segment.
Training of survey field staff on occupancy surveys, camera
trap surveys, and line transect surveys were conducted in The trained personnel walked along high probability tiger
CNP (for Chitwan-Parsa complex), BNP (for Banke-Bardia sign areas such as forest trails, fire lines, ridge lines, river
Complex), and ShNP (for Shuklaphanta-Laljhadi-Jogbuda beds and streams searching for the signs of tiger (scats,
Complex). The trained personnel were deployed in groups pugmarks, scrapes, kills and urination), prey signs (dung,
of 6-8 at strategic locations across the study area. They footprints, sightings and calls) and human disturbances
were assigned to cover the allocated grid cells –to carry such as wood cutting, lopping, grazing and signs of
out tiger habitat occupancy surveys, to setup and monitor poaching.
camera traps, and to carry out line transect surveys.
8 STATUS OF TIGERS AND PREY IN NEPAL 2018
4.2.4 CAMERA TRAP SURVEYS FOR ABUNDANCE ESTIMATION OF TIGERS 4.2.5 LINE TRANSECT SURVEYS FOR TIGER PREY DENSITY ESTIMATION
A pair of camera traps was systematically placed in Distance sampling (Buckland et al., 2001) was used for
1,643 locations of the total 2,045 grid cells (2 km x 2 km) estimating tiger prey densities. Line transects of 1.5-2
covering the entire area of tiger-bearing PAs and adjoining km lengths were systematically placed on 2 km x 2 km
forests (6,572 km2, Figure 2). The camera trap locations camera trap grid cells; areas falling in hilly terrains were
were selected based on extensive field surveys for signs avoided to adhere to the straight-line assumption of
of tiger such as pugmarks, scats, scrape and urination. distance sampling (Figure 3). Global Positioning System
Camera trapping was carried out in shifting blocks in (GPS) locations of the start and end points of each of the
each survey site (Royle et al., 2009). Cuddeback (C1) and transects were uploaded onto GPS receivers prior to the
Panthera (V5 and V6) digital cameras were used to obtain survey and the straight line was navigated following the
high quality images for individual tiger identification. bearing using Suunto compass and GPS receiver.
Cameras were programmed to take 3 pictures per trigger
with no delay (FAP mode) using white flash. The camera Line transect survey was conducted either on foot or on
traps were deployed for 15-20 nights in each of the grid elephant back. Each transect was traversed by two people
cells. with no delay (FAP mode) using white flash. between 0630 hours and 0930 hours; each transect was
surveyed twice. Elephants were only used in tall flood radial distance to the animal/center of group and the
plain grasslands. The following data were recorded along bearing of the group and GPS locations of each sighting.
the transects - bearing, species sighted, group size, total Range finders and Suunto compass were used to measure
numbers of adult and young individuals in each group, radial distance and animal bearing.
10 STATUS OF TIGERS AND PREY IN NEPAL 2018
© DNPWC/WWF Nepal
© DNPWC/WWF Nepal
STATUS OF TIGERS AND PREY IN NEPAL 2018 11
5 . DATA ANALYSIS
5.1 TIGER HABITAT OCCUPANCY MODELLING Tiger abundance estimates were derived using Maximum
likelihood (ML) and Bayesian (B) based spatially-explicit
Standard occupancy modelling approach (MacKenzie et capture-recapture (SECR) (Royle et al., 2009; Efford and
al., 2002) was used to estimate occupancy (Psi), detection Fewster, 2012; Efford, 2018).
probabilities (p) and including modelling covariate
effects on detectability and occupancy. A detection/non- Input files i) a spatial capture history matrix, ii) a trap
detection history matrix was generated and imported into layout matrix and iii) a habitat mask excluding non-
the program PRESENCE 12.7 (Hines, 2013). Multiple habitat areas were prepared and analyzed using ‘secr’
season model (Hines et al., 2014) was used to estimate the package (version 3.1.6, Efford, 2018) in the R statistical
trend in occupancy dynamics across the landscape both environment (version 3.5.0, R Development Core Team,
spatially and temporally. For standard estimates, spatial 2018).
auto-correlation between sampled replicates was tested
using Hines et al., 2010 (for single season) and Hines et Range of models with biologically plausible covariates
al., 2014 (for multiple season) model respectively. on detection probability (g0) and space range (sigma)
were considered. The effects of time factor (t), time trend
5.2 TIGER POPULATION ABUNDANCE ESTIMATION (T), animal’s learned response (b), transient response
(B), animal x site learned response (bk), animal x site
Camera trap surveys are now a well-established transient response (Bk), and two-class mixtures (h2)
methodology for density and abundance estimation of were specified and modelled for both detection and
elusive carnivores (Karanth et al., 2008; Kelly and Holub distribution. All models were ranked based on Akaike’s
2008). Recent development of spatial capture-recapture Information Criterion (AICc) and model-averaging was
methods has led to greater clarity in abundance estimation done with models having delta AIC<2 to determine
by integrating spatial or “location” information of animal population estimates for each site.
captures. This involves identification of tigers based on
their unique stripe patterns, developing a capture history Since the PAs surveyed are contiguous (e.g. BaNP share
matrix detailing tiger ID, capture location and sampling its boundary with BNP, and PNP with CNP), taking
occasion over the sampling period (Karanth and Nichols population size (Ń) of PAs buffer would overestimate the
1998) and analysis of capture history data using maximum population. Therefore, SECR models were fitted using the
likelihood (Efford and Fewster, 2012; Efford, 2018), or stable buffer size first and then population estimates were
Bayesian framework (Royle et al., 2009; Gopalaswamy exclusively derived for the effectively sampled area or the
et al., 2012). The data is also amenable to analysis in a ellipse that contained all the detectors (camera traps).
non-spatial framework and can be used for conventional
mark-recapture analysis (White & Burnham, 1999). 5.3 TIGER DENSITY ESTIMATION
Individual tigers were visually identified by field SECR models under Bayesian framework using Markov-
technicians and trained biologists at three levels (i.e.- Chain Monte Carlo (MCMC) methods were used to
field technicians, field implementation committee and estimate site-specific tiger densities in SPACECAP
wildlife biologists) by thoroughly examining all the (version 1.1.0) (Gopalaswamy et al., 2012) in R 3.4.0 (R
images obtained. Only adult tigers (animals captured Core Team, 2017). Three input files - “animal capture” file
independently without mother) were used in the analysis detailing trap location, animal ID and sampling occasion,
(Karanth et al., 1998). Individual tigers were given a “trap flag” file, and, “habitat mask” were prepared. Trap
consistent ID based on the national tiger database. The flag was created and included in the model to specify
tigers were also segregated by sex where possible. Tiger active days of each camera trap station. This incorporated
images from protected areas with shared boundary the block sampling design and explicitly accounted
were also compared and common tigers were identified. for dysfunctional cameras on account of theft, wildlife
Common tigers were assigned to the protected area with damage or malfunction. Habitat mask was created for
maximum spatial coverage to estimate site level tiger area that included camera trap array (MCP: Minimum
population abundance. Convex Polygon) surrounded by a buffer of half mean
12 STATUS OF TIGERS AND PREY IN NEPAL 2018
maximum distance moved (1/2 MMDM) by the tigers camera locations of 2018 within the camera trap polygon
as range beyond this were all dominated by human of 2013 was used and tiger densities were estimated for
settlements. Pixelated habitat mesh size of 0.3364 all the study sites in TAL following the similar approach
km2 was used (Karanth et al., 2008). Models with four in SPACECAP.
different combinations - trap response present, trap
response absent, half normal and negative exponential 5.4 TIGER PREY DENSITY ESTIMATION
detection functions were used to fit the data.
Line transect data were analyzed using the program
MCMC simulations with over 1,00,000 iterations, burn- DISTANCE version 7.1 to obtain density estimates of
in of 15,000-25,000 and thinning rate of 1-5 and data prey species (Buckland et al., 2001; Thomas et al., 2010).
augmentation value of 5 times the number of animals Due to low detection of prey, multiple years’ data (2016-
captured was set for running the site-specific analysis. 2018) of similar season were used to estimate the prey
Geweke diagnostic scores (-1.64 to 1.64) was used to check density (Kumar et al., 2018). These yielded estimates
the convergence of chains and data fit (Gopalaswamy et of the density of principal prey species for each site.
al., 2012). Pixelated map showing the tiger density was Observation of all the species was pooled for fitting global
produced for each of the sites in ArcGIS (Ver. 10.1). detection function. For species with sufficient detections,
detection function was fitted at the species level. Chi
Two separate density estimates were derived by square goodness of fit test was used to assess the fit of the
accounting the area sampled in 2018 and 2013. For model, and the best model from the subset of models was
comparison with density estimates of 2013, site-specific selected using lowest AIC value.
© WWF Nepal
STATUS OF TIGERS AND PREY IN NEPAL 2018 13
6 . RESULTS
6.1 TIGER HABITAT OCCUPANCY (16,261 km2) in the landscape, an estimated 11,057.5 km2
(SE 663.5 km2) was occupied by tigers as of 2018.
6.1.1 SAMPLING EFFORT AND TIGER SIGN DETECTION
The team surveyed 112 grid cells in TAL (Figure 4) with a In addition, the modelled occupancy inside the PAs was
sampling effort of 2,838 km with an area of 16,261 km2. 0.98 (SE 0.06) covering 6,828 km2 and 0.60 (SE 0.11)
A total of 494 unique tiger signs were detected. Of these, 5,576 km2 outside PAs. The detection probability inside
409 (83%) signs were recorded within 45 grid cells located PAs was 0.79 (SE 0.03) and outside PAs 0.37 (SE 0.08)
inside the juristic boundary of PAs (Core area and Buffer respectively.
zone) and 85 (17%) signs were recorded outside PAs in
67 grid cells. The majority of the detections outside PAs 6.2 TIGER ABUNDANCE
(41 tiger signs) were recorded in 25 grid cells located in
corridors. Overall, tiger signs were recorded in 12 districts 6.2.1 SAMPLING EFFORT AND TIGER CAPTURES
(Bara, Parsa, Makwanpur, Chitwan, Nawalparasi, Dang, A pair of cameras were deployed in 1,643 grid cells (2 km
Salyan, Banke, Bardia, Surkhet, Kailali and Kanchanpur) x 2 km) across five tiger bearing protected areas and their
out of 18 districts surveyed. adjoining forests. Tigers were captured in 482 (29%) grid
cells (Table 2). Total effective sampling area (ESA) was
6.1.2 TIGER HABITAT OCCUPANCY AND DETECTION PROBABILITY estimated at 12,356.64 km2.
The naïve tiger occupancy was 0.6 where tiger signs were
detected in 68 out of 112 grid cells. The modelled occupancy Camera trapping effort of 27,829 days across all sites
(proportion of area occupied) in the landscape was 0.68 resulted in tiger trap rate of 0.04 per trap day (4.2%)
(SE 0.06) with an estimated detection probability of with 4,388 tiger images and 1,136 independent tiger
0.73 (SE 0.03). Out of the total potential habitat of tigers detections. Individual tigers were identified using stripe
Table 2. Number of grid cells surveyed and number of grid cells with tiger captures in each site.
Number of surveyed camera trap Number of grid cells with tiger
Site
grid cells captures
PNP and adjoining forests 305 49 (16%)
CNP and adjoining forests 509 199 (39%)
BaNP and adjoining Forests 254 38 (15%)
BNP and adjoining Forests 323 149 (46%)
ShNP and adjoining forests 252 47 (19%)
Total 1,643 482 (29%)
14 STATUS OF TIGERS AND PREY IN NEPAL 2018
patterns of all available “both flank” pictures and either 6.2.2 TIGER POPULATION ABUNDANCE ESTIMATES
“right or left flank pictures” for each of the study sites. The estimated abundance of tigers in PNP and adjoining
Juveniles and cubs (N=28) captured were not included in forests is 18 (16-24), CNP and adjoining forests is 93
the analysis (Table 3). (89-102), BaNP and adjoining forests is 21 (18-30),
BNP and adjoining forests is 87 (82-97) and ShNP and
These resulted into 209 individual tigers (PNP and adjoining forests is 16 (15-21). The SECR-ML based tiger
surrounding forests-15, CNP and surrounding forests-85, population estimate for each tiger bearing PAs and the
BaNP and surrounding forests-17, BNP and surrounding adjoining forests is provided in Table 4. Summing up the
forests-77, ShNP and surrounding forests-15) (Table 3) site-wise estimates, the forests of TAL-Nepal support 235
including 78 males, 118 females and 13 of unknown sex. tigers as of 2018. The details of the model used, and the
Site-wise breakdown of tiger numbers is provided in Table real parameters are provided in Annex-1.
3 with respective sampling efforts.
Table 3: Site wise sampling effort and the minimum tigers (Mt+1) identified
Number
Survey Effective
Number Number of of Adult
effort sampling Adults Adult
Site of tiger independent individual unknown Cubs
(trap area males females
photos detections tigers sex
days) (km2)
captured
PNP and
adjoining 4,810 3,634.50 294 102 15 5 10 - 3
forests
CNP and
adjoining 8,433 2,281.10 1,744 480 85 30 50 5 12
forests
BaNP and
adjoining 4,503 2,311.70 565 61 17 6 10 1 0
forests
BNP and
adjoining 5,479 2,832.80 1,554 404 77 28 42 7 11
forests
ShNP and
adjoining 4,604 2,154.30 231 89 15 9 6 - 2
forests
Total 27,829 12,356.6 4,388 1,136 209 78 118 13 28
The tiger population estimates generated using other 6.4 PREY DENSITY ESTIMATES
programs have been provided in Annex-2 for better
comparison with the earlier surveys. In total 1,294 transects were conducted covering 2,485
km. Tiger prey species recorded during the survey
6.3 TIGER DENSITY ESTIMATES included four deer species (spotted deer, sambar, hog
deer, barking deer), two antelope species (blue bull and
Data convergence was achieved for results of all the study four-horned antelope), wild boar, gaur, and two primate
sites accounting Geweke diagnostic score with other real species (rhesus macaque and langur).
parameters. The mean posterior density of tigers per 100
km2 in PNP and adjoining forests was 0.92 (SD 0.15), Combined density of prey (all prey per km2) varied
CNP and adjoining forests was 3.28 (SD 0.19), BaNP and between 8.1 and 77.51 animals per km2 across the sites.
adjoining forests was 0.97 (SD 0.12), BNP and adjoining The combined prey density per km2 in PNP and adjoining
forests was 4.74 (SD 0.28) and ShNP and adjoining forests is 22.03 (SE 3.8), CNP and adjoining forests is
forests was 0.96 (SD 0.14) respectively. 70.7 (SE 7.49), BaNP and adjoining forests is 8.1 (SE 1.6),
BNP and adjoining forests is 77.51 (SE 6.56) and ShNP
The density estimates with 95% confidence intervals are and adjoining forests is 68.04 (SE 6.95) respectively.
provided in Table 5. The summaries of real parameters
for each of the sites are provided in Annex-3 The site-wise sampling effort, number of observations
and prey density estimates are provided in Table 6
The pixelated tiger density map produced by combining and species-wise details and species-wise prey density
site-wise pixel values generated by program SPACECAP estimates are provided in Annex-7.
is provided in Figure 5.
Table 5: Tiger density estimates for the tiger-bearing protected areas including buffer zones, adjoining
forests and corridors.
Table 6.Overall prey density estimates in tiger bearing protected areas and adjoining forests.
Number
Effort No. of Density
Site of SE CV (%) 95% CI
(km) obs. (per km2)
transects
15.66 -
PNP and adjoining forests 482 248 194 22.02 3.8 17.48
30.96
57.49 -
CNP and adjoining forests 331.6 175 367 70.7 7.49 10.59
87.05
5.46 -
BaNP and adjoining forests 647 304 99 8.1 1.6 20.25
12.01
65.66 -
BNP and adjoining forests 745 414 776 77.51* 6.56 8.47
91.49
55.70 -
ShNP and adjoining forests 279 153 412 68.04* 6.95 10.22
83.10
Total 2,484.6 1,294 1,848
* denotes density estimates of all prey excluding swamp deer that was not detected during the line transect survey in both Bardia National Park and
Shuklaphanta National Park. SE: Standard Error of Mean, CI: Confidence Interval, CV: Coefficient of Variation (SD/Mean)
Figure 5. Tiger density within the Terai Arc Landscape, Nepal. The density map is composed of pixels (0.336 km2) representing potential activity centers of individual tigers.
STATUS OF TIGERS AND PREY IN NEPAL 2018 17
18 STATUS OF TIGERS AND PREY IN NEPAL 2018
6.5 HABITAT USE OF TIGERS OUTSIDE PAs 2 unknown) were observed in Someshwor hill forest
(CNP BZ), four tigers (2 male,2 female) were captured
In addition to the core tiger-bearing protected areas, the from Barandabhar corridor (CNP BZ) and two female
known tiger distribution range (buffer zones, corridors tigers were captured in the forests of Nawalparasi close to
and adjoining forests) within TAL were sampled using Binayee and Madhyabindu.
camera traps. Therefore, in the site-specific, spatially
explicit capture-recapture estimates of PAs and adjoining In western Terai, one tiger with a large spatial range
forests, the tigers with activity centers outside PAs have extending across BNP and BaNP was also observed in Dang
been included in the analysis. Altogether, twenty-two forest near Lauki guard post (Table 7; Annex-8: Figure
tigers were recorded outside PAs. However, only nine 8). Altogether, thirteen and three tigers were captured
tigers were exclusively captured outside PAs (Table 7). in Khata and Karnali river corridor respectively. Among
Of these, adult female captured in the district forest of the thirteen tigers captured in Khata corridor, four were
Nawalparasi was photographed with three healthy cubs. captured exclusively and nine tigers (6 male, 3 female) had
The rest of the fifteen tigers captured in district forests their home ranges extended within the habitats in BNP
had their home range extended to core and buffer zones and BNP BZ (Table 7; Annex-8: Figure 9).
(Table 7). Habitat use of tigers therefore was observed
to be minimal outside the juristic boundary of protected In far-west Terai, only one male tiger was recorded
areas during the study period (see maps in Annex-8). from outside ShNP in Laljhadi corridor. This male was
recorded from forests that spanned north-eastern part
In central Terai, a male and female tiger were recorded in of ShNP and fragmented forest patches in Laljhadi
the collaborative forests in Parsa, three tigers (1 female, corridor.
Table 7.Number of tigers captured within and outside the juristic boundary of protected areas and
exclusively outside protected areas.
S. N Site Male Female Unknown sex Total
1 Parsa collaborative forest* and PNP 1 1 - 2
2 Nawalparasi (exclusively) 0 2 - 2
3 Dang forest, BNP and BaNP 1 - - 1
4 Karnali river corridor (exclusively) - - 3 3
5 Khata corridor (exclusively) - 2 2 4
6 Khata, BNP BZ and BNP 6 3 - 9
7 ShNP core and Laljhadi corridor 1 - - 1
Total 9 8 5 22
*Collaborative forest: A forest jointly managed by national government, local government and local communities
STATUS OF TIGERS AND PREY IN NEPAL 2018 19
7 . DISCUSSION
7.1. TIGER DISTRIBUTION, HABITAT USAGE AND OCCUPANCY
KEY FINDINGS:
Altogether 494 unique tiger signs were recorded by
habitat occupancy survey carried out across TAL. About 68% of the total potential habitat was occupied
Maximum signs (83%) were recorded within PAs and 17% by tigers across TAL, covering 11,057 km2.
of tiger signs were detected outside PAs in TAL. Similarly,
tigers were captured in 482 out of 1,643 grid cells which Camera trap survey identified a minimum of 209
is only 29% of the total area surveyed. Combining both individual tigers.
tiger signs and tiger captured locations in camera traps,
tigers were recorded in 12 districts. However, majority PA wise tiger population was estimated at PNP-18,
(96%) of the photo-captured tigers (N=200) were mostly CNP-93, BaNP-21, BNP-87 & ShNP-16. This summed
confined within protected areas and 9 tigers (4%) were up to 235 tigers in Nepal. 88% of the estimated tiger
captured in forests outside the protected areas in TAL. population were photo-captured (Mt+1) in the survey.
Thirteen individual tigers that used core and buffer
zones also had their territorial range extended to forests Nine tigers were found exclusively in areas outside
outside PAs. Twenty-four tigers were captured from five juristic boundary of PAs with 13 tigers using forested
biological corridors viz: Laljhadi (1), Khata (13), Karnali habitats both within and outside PAs.
(3), Barandabhar (4) and Someshwor hill forest (3). Only
few tiger signs were recorded from Kamdi and Basanta Tiger density/100 km2 ranged from lowest 0.9 in BaNP
corridors. to highest 4.7 in BNP.
Modelled tiger occupancy in the landscape is showing an Prey density/km2 ranged from lowest 8.1 in BaNP to
increasing trend. Between 2009-2013, there was a 47% highest 77.5 in BNP.
(λ=1.47) increase in occupancy. Similarly, between 2013-
2018 there was a 12% (λ=1.12) increase in occupancy
across the landscape. The likely reason for the increase
in occupancy can be attributed to the increasing tiger PAs more conducive for tigers and to facilitate their
population that are expanding, establishing their safe dispersal, measures should be focused towards
territories in areas previously unoccupied [colonization improving habitat quality, increasing prey population
probability, γ- 0.45 (SE 0.06)] and ongoing restoration and minimizing human disturbances at the level that
efforts. comply with tolerance level of tigers. Replicating success
of community forestry as seen in Khata corridor could
In totality, tigers occupied 11,057 km2 (68%) of the provide a potential solution that benefits both tigers and
available habitat (16,261 km2) in the landscape. Within people.
PAs tigers occupied 6,828 km2 (98%) of the available
habitat (6,968 km2). However, tigers occupied only 5,576 7.2 TIGER ABUNDANCE AND DENSITY
km2 (60%) of the available habitat 9,293 km2 outside the
PAs. The large tracts of forest exist outside the protected 7.2.1 METHODS USED, AND THE EXTENT OF AREAS COVERED
areas in TAL, but majority of these forest patches face SECR-ML and SECR-B are commonly used techniques in
high anthropogenic pressure. The ground forest cover is deriving population and density estimates (Gopalaswamy
literally non-existent with high cattle grazing, and prey- et al. 2012; Royle et al. 2009; Elliot and Gopalaswamy
base is extremely low to support resident tigers. Unlike 2016). In the present survey, tiger population estimates
PAs, these forested habitats have minimal protection, and were derived using SECR ML while density estimates
therefore face risk of becoming a sink for tigers. Thus, the were derived using SECR-B. In 2013, both population
existing limited use of forests by tigers outside PAs can and density estimates were derived using SECR-B
be enhanced through protection and other management without considering the overlaps between the contiguous
interventions similar to PAs. To make habitat outside protected areas in Nepal and India (e.g. CNP, PNP and
20 STATUS OF TIGERS AND PREY IN NEPAL 2018
Valmiki Tiger Reserve (VTR)). This could have likely and adjoining forests-17, BNP and adjoining forests-77
overestimated the population abundance in the PAs that and ShNP and adjoining forests-15) compared to 142
shared contiguous habitats. The results of the present individuals (40 males, 102 females)
survey derived by both SECR-ML and SECR-B have
addressed the overlap issues and the estimates derived In PNP, BaNP and BNP where annual/bi-annual tiger
are similar with 95% CI overlaps (Table 8). In estimating monitoring was conducted, minimum population of
PA wise abundance coefficient of variation (CV %, tigers has increased against 2013 baseline. However,
Table 8) is relatively lower for SECR-ML (av. CV: 7.5%) since 2013, Mt+1 has remained stable in ShNP largely
as compared to SECR-B (av. CV: 11%). Therefore, PA because of the male biased sex ratio (2018 survey-1.5:1),
abundance estimate (∑ 235 tigers) have been reported that could be other reason impeding the growth in tiger
using SECR-ML. Similarly, as program CAPTURE and population.
Mark were used in 2009 and 2013 surveys, the results
obtained from these programs have also been reported In PNP, the highest number of individual tigers captured
for readers knowledge. For density estimates, the results was 19 in 2016. This has dropped down to 15 individuals
from both SECR-ML and SECR-B have been reported in the 2018 survey. Among the tigers not captured in
(Annex-6). PNP during the 2018 survey, two were earlier captured
in Someshwor hill corridor forest (a transboundary
The 2018 survey extensively covered potential tiger corridor linking CNP, PNP and VTR of India), indicating
habitat of TAL, Nepal. The sampling effort in this survey that these could be transients. It is only through regular/
was maximized by covering most of the known records annual surveys that the fate of individual tigers lost
of tiger distribution based on findings of annual tiger through natural death or dispersal and also recruitment
surveys since 2013 and thus, the extent of the area from births or immigration is understood, which can
covered by camera traps increased from 1,039 grid cells provide a rigorous audit for the successes or failures in
in 2013 survey to 1,643 grid cells in 2018 (58%). Tiger tiger conservation programs. Hence annual surveys are
captures were recorded in 29% (482) of the 1,643 grids. recommended to provide insights in tiger population
However, the 604 new grids added in 2018 accounted for dynamics.
only six tigers indicating minimal impact on the overall
status change. 7.2.3. TIGER ABUNDANCE ESTIMATES
Status of tigers in Nepal is estimated at 235 individuals
7.2.2 TREND IN MINIMUM POPULATION BASED ON INDIVIDUAL TIGER in 2018 which is simply a summation of PAs and their
CAPTURES adjoining forests-wise estimation without estimated
This study reports a minimum population of 209 variance. In 2013, a total of 198 tigers were estimated in
individual adult tigers; 78 males, 118 females and 13 Nepal (Dhakal et al., 2014). A naïve comparison indicates
unknow sex, from across the study sites (PNP and an increase of approximately 19% within the four-year
adjoining forests-15, CNP and adjoining forests-85, BaNP period. There has been notable increase in BNP, BaNP and
CV: Coefficient of Variation, Mt+1: Minimum individual identified, SE: Standard error of Mean, SECR: Spatially Explicit Capture Recapture, ML: Maximum
Likelihood, B: Bayesian
STATUS OF TIGERS AND PREY IN NEPAL 2018 21
120
80
60
40
20
0
2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019
-20
ShNP BNP BaNP PNP CNP
Linear (ShNP) Linear (BNP) Linear (BaNP) Linear (PNP) Linear (CNP)
PNP (including adjoining forests) whereas the population mostly confined within the boundaries of protected areas
has remained stable in ShNP and adjoining forests and and associated buffer zones. Therefore, to determine
declined in CNP and adjoining forests. Further, unlike the true change in density from 2013 to 2018, density
the other protected areas, CNP lacks annual monitoring estimates were generated for the effective sampled area
data to confirm the actual trend in population estimates. of 2013 for better comparison
In BaNP and adjoining forests, tiger population has Detailed outputs from SECR-B analysis for 2013 and
increased by more than five-fold i.e. from 4 (3-7) tigers in 2018 are provided Annex-4. The pixelated tiger density
2013 to 21 (18-30) tigers in 2018 (P =0.001). In Bardia, maps for 2013 and 2018 produced by combining site-
tiger population almost doubled from 50 (45-55) in wise pixel values generated by program SPACECAP are
2013 to 87 (82-97) in 2018 (P=0.002) Likewise, PNP provided in Annex-5.
and adjoining forests, the tiger population significantly
increased from 7 (4-7) in 2013 to 18 (16-24) tigers in 2018 Site-level tiger densities (individuals per 100 km2) from
(P=0.03). In contrast, the estimated population of CNP 2013 to 2018 increased from 0.65 (SD 0.28) to 1.49 (SD
has marginally declined from 120 (98-139) in 2013 to 93 0.23) in PNP (P =0.088), 0.16 (SD 0.05) to 1.38 (SD 0.17)
(89-102) tigers in 2018 (P = 0.06). The population has in BaNP (P =0.002) and 3.38 (SD 0.19) to 5.4 (SD 0.29)
remained stable in ShNP and adjoining forest with 17 (13- in BNP (P =0.03). The tiger density estimates for CNP
21) estimated tigers in 2013 to 16 (15-21) tigers in 2018 has marginally declined from 3.84 (SD 0.34) to 3.81 (SD
(P =0.86). 0.25) (P=0.97) whereas it has remained stable in ShNP at
3.4 (SD 0.65) (P=0.92).
7.2.4. TIGER DENSITY
Tiger density (no. of tigers per 100 km2) ranged from 0.9 7.2.5. FACTORS GOVERNING TIGER POPULATION ABUNDANCE AND DENSITY
in BaNP (lowest) to 4.7 in BNP (highest). Tiger density The conservation of tigers is dependent on appropriate
estimates for 2018 was lower compared to 2013. This was protection measures, prey densities, habitat connectivity,
because the additional habitats (~2400 km2) sampled habitat management of the critical habitats, park-people
in 2018 supported only a few tigers (n=6). Tigers were relationship and human-tiger interaction across TAL,
22 STATUS OF TIGERS AND PREY IN NEPAL 2018
2013 2018
6
5
5.4
3.81
3
3.84
3.38 3.48
2 3.4
1 1.49
1.38
0.65 0.6
0
PNP* CNP* BaNP* BNP* ShNP*
Figure 7. Site-level tiger density estimates with respective standard error bars for 2013 and 2018.
* denotes respective PAs and adjoining forests.
Nepal. These are discussed in detailed below. recorded from as early as 2014 and the cubs born in 2014,
2015 and 2016 have survived to adulthood. Therefore,
The significant increase in tiger population abundance a detailed study to understand the meta-population
in BNP and its adjoining forests can be attributed to dynamics of tigers is suggested as future priorities for
enhanced protection measures, increased support from managing the tigers in Banke-Bardia complex.
communities in buffer zone and corridors and regular
practice of habitat management by park authorities. The In PNP, study by Lamichhane et al. (2017) reported a total
park has been equipped with 10 additional protection of 25 individuals from the annual surveys carried out from
posts strategically placed at Lamidamar, Kalinara, 2013 to 2016. Of these, 10 tigers were earlier reported in
Ratamate, Lekhparajul, Taranga, Thuloshree, Sotkhola, CNP thus PNP has benefitted through tiger dispersal from
Banspani, Telpani and Okhariya in the last four years. the adjoining CNP. The tiger recovery can also be credited
Patrol coverage throughout the PA and buffer zone has to voluntary relocation of villages from the park; these areas
been achieved through the implementation of real-time now support higher numbers of ungulates and breeding
SMART covering 31 protection/army posts. BNP also tigers. In 2015, the core area of PNP was extended by 128
benefits from the adjoining Khata corridor that provides km2 to 627 km2. Park protection and management have been
direct habitat linkage to Katerniaghat Wildlife Sanctuary enhanced through construction of park posts in strategic
(KWS) in India through regular transboundary movement locations at northern part of the park including Annexed
of animals ensuring genetic exchange between the two area and increased number of Nepali Army personnel in the
populations (Chanchani et al.,2014; Thapa et al., 2018). park.
BNP is also acting as a source site for BaNP; three tigers (2 Chitwan-Valmiki-Parsa complex holds a significant
male, 1 female) born and raised in BNP have overlapped population of tigers within the eastern part of Terai Arc
territories with BaNP. Thus, the increase population in Landscape. CNP acts as source population replenishing
BaNP is partially contributed by the dispersal of tigers tigers to the adjoining and contiguous protected areas
from adjoining BNP. There has also been improvement in of PNP in the east and VTR in the south. Recent survey
park management and protection over the last four years. has showed a decline in tiger population in CNP and
Since 2013, total of 15 park posts were constructed along adjoining forests. This could be attributed to 1) dispersal
the northern border and in strategic sites in the south. In of tigers to neighboring sites (PNP and VTR), 2) increased
2013, real-time SMART was implemented starting with intraspecific competition, 3) human-wildlife conflicts
two park posts and extended to fourteen posts to date. (two were reported killed in retaliation), 4) poaching (one
Habitat improvement has also been undertaken including case), and 5) natural disasters such as the catastrophic
restoring and creating water holes (n=17) and managing flood of 2017.
grasslands (147 hectares). Signs of tiger breeding have been
STATUS OF TIGERS AND PREY IN NEPAL 2018 23
The loss of eight tigers due to territorial fights between note, existing method of line transect survey was heavily
2014-2018 suggests space could be a limiting factor concentrated in Terai, thus the species like sambar
within CNP. Two tigers were killed in retaliation (Rusa unicolor), four-horned antelope (Tetraceros
suggesting further interventions are needed in garnering quadricornis), ghoral (Naemorhedus goral) and
community support. One recorded case of poaching also Himalayan tahr (Hemitragus jemlahicus) that are well
captures the prevailing threat of poaching in the park. adapted to undulating Churias were not accounted for.
In addition, massive flood in August 2017 swept away
unaccounted numbers of wildlife in CNP. Park authorities Other possible reason for decline in prey density could
rescued a total of 10 rhinos from downstream riverine be the catastrophic flood in Chitwan in 2017 and Bardia
habitat in India. Furthermore, the flood could have in 2015. In one event alone, 28 hog deer carcasses were
affected significant number of prey species consequently found from a single site, after being washed away by the
affecting the tiger’s space use. However, a detailed study swelling Narayani and similar cases were reported in
(including carrying capacity) is needed to understand the Babai valley of Bardia in 2015.
tiger population in CNP.
Moreover, swamp deer which occurs in high density in
Tiger populations in PNP and VTR have increased ShNP has not been accounted as they were not sighted
substantially in recent years. The turnover rate (defined in the transect survey. This could have contributed
as loss and gain) of tigers in CNP based on camera trap to lower density estimates in ShNP. Recently, prime
data was found to be almost 70% in the last 4 years with wetlands (Rani taal, Salgaudi taal and Kalikich taal)
only 19 tigers recorded in 2013. Of the 31 individual tigers have deteriorated in the core area of ShNP primarily
captured in PNP between 2013-2018, 10 were earlier due to siltation, while grasslands have been encroached
captured in Chitwan and the rest 12 documented were not by woody perennials along Chaudhar river. This may
born and raised in PNP suggesting high dispersal from also have contributed in the decline of prey base when
CNP. In VTR, tiger population has increased from 10 compared to 2013 study.
individuals in 2008 to 31 individuals in 2016/17 (Source:
WWF India). Prey density in BaNP and PNP is notably lower than other
three tiger bearing PAs (CNP, BNP and ShNP) in Nepal.
The major reason for the stagnant population of ShNP is Larger parts of BaNP and PNP are covered by Churia
the male biased sex-ratio (1.5:1). The females are known forests, while the remaining lowland area is dominated
to give births occasionally but the survival of cubs by homogenous Sal forests that have lower potential to
to adulthood has been a major factor. Two cubs died support high prey density. CNP, BNP and ShNP have
because of possible infanticide in 2017 and one sub- extensive riparian forests and tall floodplain grasslands
adult tiger was killed by the dominant male in territorial regulated by annual floods where ungulates can reach
fight (ShNP, 2018). Furthermore, southern section of their highest densities. BaNP and PNP are also extremely
national park, along Lagga Bagga, offers an opportunity dry due to the physical characteristics of Bhabar region
to tigers to disperse into India. Security of dispersing and therefore water may be another limiting factor for
tigers and/or resident tigers along the transboundary prey species.
protected areas requires special attention.
Sightings of prey outside PAs were negligible owing to
7.3. PREY DENSITY large scale disturbances. Hence, prey recovery in the
surrounding forests of PAs with focus on regulating
Prey density estimates marginally declined across the and reducing disturbances should be emphasized in the
survey sites as compared to 2013. The prey densities management plan.
per km2 decreased from 25.33 to 22.02 (SE 3.8) in PNP
(P=0.6), 73.63 to 70.7 (SE 7.49) in CNP (P=0.8), 10.27 Based on documented records, the issue of prey poaching
to 8.1 (SE 1.6) in BaNP (P=0.6), 92.6 to 77.51 (SE 6.56) across the sites outside PAs cannot be overlooked
in BNP (P=0.3) and 78.62 to 68.04 (SE 3.8) in ShNP either, and therefore, it demands a timely investment to
(P=0.12) with no significant difference (Dhakal et al., elevate protection measures in all such vulnerable sites.
2014). Undoubtedly, low prey density will have an impact on
tiger populations as tiger densities are mediated mainly
The lower density estimates from the present survey by prey abundance (Karanth et al., 2004).
could be the results of difference in the study time frame.
Field survey for this study was conducted in mid-winter In addition, roads and highways bisect several forests
(Jan-Feb) when it was mostly misty until late morning, in TAL (tiger bearing PAs and critical corridors)
leading to poor visibility. The 2013 survey was conducted causing fatalities of predators and prey. On average, 26
in March-April when the weather was clear and the road kills in BNP and 55 road kills in BaNP occur per
grass much shorter (Dhakal et al., 2014). On another year, comprising 70-80% of prey species (chital, wild
24 STATUS OF TIGERS AND PREY IN NEPAL 2018
boar, sambar, barking deer and langur). The proposed the fragile corridors. Usage of underpasses by wildlife has
expansion of existing highways to four-lane and the been documented in Barandabhar corridor; placement of
construction of Hulaki road that will pass through crossing overs (under or over passes) along vulnerable
hundreds of kilometers along the Terai forests could sections could mitigate some losses.
further worsen this situation by severing connectivity in
STATUS OF TIGERS AND PREY IN NEPAL 2018 25
8.4 ADDRESSING THE IMPACT OF LINEAR INFRASTRUCTURES relief to victims of conflicts needs to be implemented.
Simplification of the claiming procedures for quick relief
The number and extent of infrastructure is likely to and enhancing transparency and efficiency at all levels is
increase within TAL, Nepal, further fragmenting and recommended.
severing the already fragile forested habitats. Efforts
need to continue to ensure the right balance between Focus must be made on understanding conflicts (both
conservation and development through i) engagement social and ecological dimensions), and on monitoring,
with policy makers to ensure that inviolate zones responding and preventing conflicts, as compensation
(critical tiger habitats) are avoided; ii) engagement measures in the long run will become a financial
with developmental agencies to promote smart green liability to the government. Promoting private
infrastructures (SGIs), and integration of effective insurance schemes (human, livestock, property and the
mitigation measures in the infrastructural planning crops) would help reduce the perpetual dependence on
process; and iii) developing appropriate mitigative the government. Likewise, the government’s Rapid
measures in existing linear infrastructures (for example- Response Team (RRT) in Chitwan and Bardia NPs
automated barriers to maintain speed limits, digital need to be strengthened and scaled up to additional
tracking of passing vehicles, construction of over/ sites to provide pro-active response to conflict cases.
underpasses or guiding fence in vulnerable sites to Rehabilitation of individual animals rescued from
provide safe passage for wildlife). conflicts should be carefully done and supported by
long-term monitoring.
8.5 ENHANCING TRANSBOUNDARY CO-OPERATION
8.7 STRENGTHENING ANNUAL TIGER, PREY AND HABITAT
Nepal shares approximately 800 km stretches of open MONITORING PROGRAMS
border with India, presenting opportunities for the two
countries to collaborate closely in conservation. The joint Data on survivorship, reproduction and social
tiger monitoring exercise carried out by the two countries structure in tiger populations is possible only through
in 2013 identified at least 10 tigers sharing the habitat standardized long-term monitoring programs. Annual
across borders, highlighting the need to manage tigers as monitoring would help keep pulse on core populations.
a transboundary metapopulation. However, the porous The management is recommended to establish long-
border between two countries adds challenges to the term monitoring programs in the respective PAs to keep
conservation of these endangered wildlife. Transboundary constant surveillance of the tiger population - increase
cooperation needs to be strengthened through intelligence (through new births or immigration) and losses (due to
sharing, joint patrols, standardized wildlife monitoring natural death, poaching or emigration). For Banke-Bardia
programs and data sharing, restoration and management and Chitwan-Parsa which serve as ecological units, it is
of the transboundary corridors, addressing the threats recommended to derive complex wise estimates for tiger
posed by infrastructure and knowledge sharing. population for greater ecological insights and to address
the issues of area overlap.
8.6 CREATING SAFE ENVIRONMENT FOR PEOPLE AND WILDLIFE
Similarly, prey monitoring is recommended during the
At the national level, tigers contribute only around time of the year when the visibility is optimal in stratified
2.78% of the total human-wildlife conflict compared to sampling blocks as per the habitat types for robust
67.52%, 12.85% and 10.65% caused by elephant, rhino estimates. This is also recommended in areas that are
and common leopard respectively (DNPWC, 2018). subject to annual habitat management for monitoring the
To manage conflicts, it is pertinent to understand the impact of habitat interventions.
overall conflict dynamics and consequently create safe
environment for both people and wildlife. The survival Increasing habitat potential of the existing habitats
of these species will depend upon tolerance level of the by increasing forage productivity would have direct
local communities which is determined by how well the relevance to increasing and sustaining growing tiger
overall conflict is managed, as well as ownership by the population. Therefore, research and monitoring
communities. programs such as estimating ecological carrying capacity
for the preparation of site specific management plans
Wildlife Damage Relief Guidelines (2069 BS, third and its implementation to provide scientific guidance for
amendment 2075) provisioned for providing the monetary habitat management in PAs and forests outside PAs.
STATUS OF TIGERS AND PREY IN NEPAL 2018 27
9 . REFERENCES
Barber-Meyer, S, Jnawali, S., Karki, J., Khanal, P., Raj., Subedi, N., Pradhan, N.M.B., Malla, S., Lamichhane,
Lohani, S., Long, B., MacKenzie, D., Pandav, B., B.R., Pokheral, C.P., Thapa, G.J., Oglethorpe, J., Subba,
Pradhan, N., Shrestha, R., Subedi, N., Thapa, K., and S.A., Bajracharya, P., Y, H. (2014). Status of Tigers and
Wikramanayake E. (2012). Influence of prey depletion Prey in Nepal. Department of National Parks and Wildlife
and human disturbance on tiger occupancy in Nepal. Conservation, Kathmandu, Nepal.
Journal of Zoology 289 (1): 10–18.
Dinerstein, E. (1979). An ecological survey of the
Buckland, S.T. (2001). Introduction to distance sampling: Royal Karnali-Bardia Wildlife Reserve, Nepal. Part
Estimating abundance of biological populations. Oxford I: Vegetation, modifying factors, and successional
University Press, USA. relationships. Biological Conservation, 15 (2): 127-150.
Burnham K.P and Anderson D.R (2004), Multi-model Efford, M.G. (2018). Package ‘secr’. University of Otago,
inference: understanding AIC and BIC in model selection. New Zealand.
Social Methods Res 33:261–304.
Efford, M.G., Fewster, R.M. (2012). Estimating
Burnham, K. P. and Anderson D. R. (1998). Model population size by spatially explicit capture-recapture.
selection and inference: A practical information- Oikos 122: 918-928.
theoretical approach. Springer-Verlag, New York, USA.
Elliot, N. B., & Gopalaswamy, A. M. (2017). Toward
Chanchani P., Lamichhane B. R., Malla S., Maurya accurate and precise estimates of lion density.
K., Bista A., Warrier R., Nair S., Almeida M., Ravi R., Conservation Biology, 31(4), 934–943.
Sharma R., Dhakal M., Yadav S. P., Thapa M., Jnawali
S. R., Pradhan N. M. B., Subedi N., Thapa G. J., Yadav Gopalaswamy, A.M., Royle, J.A., Hines, J.E., Singh, P.,
H., Jhala Y. V., Qureshi Q., Vattakaven J. and Borah J. Jathana, D., Kumar, N.S., Karanth, U.K. (2012). Program
(2014). Tigers of the Transboundary Terai Arc Landscape: SPACECAP: software for estimating animal density using
Status, distribution and movement in the Terai of India spatially explicit capture-recapture studies. Methods in
and Nepal. National Tiger Conservation Authority, Ecology and Evolution 3: 1067-1072.
Government of India, and Department of National Park
and Wildlife Conservation, Government of Nepal. GTRP. (2010). Global Tiger Recovery Program (Executive
Volume – September 20, 2010). St. Petersburg
Department of National Parks and Wildlife Conservation Declaration on Tiger Conservation at the International
(DNPWC). (2018). Annual report: Wildlife Crime. Tiger Forum (‘Tiger Summit’), held in St. Petersburg,
Department of National Parks and Wildlife Conservation, Russia, November 21–24, 2010.
Babarmahal, Kathmandu, Nepal.
Hines, J.E., Nichols, J.D., Collazo, J.A. (2014).
Department of National Parks and Wildlife Conservation Multiseason occupancy models for correlated replicate
(DNPWC). (2017). Tiger and Prey Base Monitoring surveys. Methods Ecol. Evol. 5, 583–591.
Protocol (Nepal). Government of Nepal, Ministry
of Forests and Soil Conservation, Department of Hines, J. E. (2013). PRESENCE 6.2: Software to estimate
National Parks and Wildlife Conservation, Babarmahal, patch occupancy and related parameters. United States
Kathmandu, Nepal. Geological Service, Patuxent Wildlife Research Center,
Laurel, Maryland, USA.
Department of National Parks and Wildlife Conservation
(DNPWC). (2016). Annual monitoring report, Western Hines, J.E., Nichols, J.D., Royle, J.A., MacKenzie, D.I.,
Terai. Department of National Parks and Wildlife Gopalaswamy, A.M., Kumar, N.S. and Karanth, K.U.
Conservation, Babarmahal, Kathmandu, Nepal. (2010). Tigers on trails: occupancy modeling for cluster
sampling. Ecol. Appl. 20: 1456–1466.
Dhakal, Maheshwar., Karki, Madhuri., Jnawali, Shant
28 STATUS OF TIGERS AND PREY IN NEPAL 2018
Karanth, K.U., Kumar, N.S., Srinivas, V. and Mishra, H.R. and Jefferies, M. (1991). Royal Chitwan
Gopalaswamy, A. (2008). Revised monitoring framework National Park: Wildlife Heritage of Nepal. The
for Tigers Forever: Panthera sites. Technical Support Mountaineers in association with David Bateman.
Team, Tigers Forever. Bangalore, India: WCS India.
Royle J.A., Nichols J.D., Karanth K.U. and Gopalaswamy,
Karanth, K.U., Nichols, J.D., Kumar, N.S., Link, W.A., M. (2009). A hierarchical model for estimating density
Hines, J.E., (2004). Tigers and their prey: Predicting in camera trap studies. Journal of Applied Ecology 46:
carnivore densities from prey abundance. Proc. Natl. 118-127.
Acad. Sci. U. S. A. 101: 4854–8.
Sanderson, J. (2010). Framing Tiger’s Troubles:
Karanth, K.U., Nichols, J.D. (1998). Estimation of Tiger Comparing Traditional and Social Media. International
Densities in India Using Photographic Captures and Journal of Sport Communication 3 (4): 438-453.
Recaptures. Ecology 79, 2852–2862.
Seidensticker, J. (2010). Saving wild tigers: A case study
Karki, J.B., Jnawali, S.R., Shrestha, R., Pandey, M.B., in biodiversity loss and challenges to be met for recovery
Gurung, G., Thapa, Karki, M. (2009). Tiger and their beyond 2010. Integrative Zoology 5 (4): 285-299.
prey base abundance in Terai Arc Landscape, Nepal.
Department of National Parks and Wildlife Conservation Smith, J.L.D. (1993). The Role of Dispersal in Structuring
and Department of Forests, Kathmandu, Nepal. the Chitwan Tiger Population. Behavior 124: 165-195.
Kelly, M.J., Holub, E.L. (2008). Camera trapping of Srivathsa, A., Parameshwaran, R., Sharma, S., Karanth,
carnivores: trap success among camera types and U.K. (2015). Estimating population sizes of leopard cats
across species and habitat selection by species, on salt in the Western Ghats using camera surveys. Journal of
pond mountains, Giles County, Virginia. Northeastern Mammalogy 96 (4): 742-750.
Naturalist, 15(2): 249-262.
Sunquist, M. E. (1981). The social organization of tigers
Kumar, C.V.P., Reddy, K.S., Srinivas, M.A.S.(2018). (Panthera tigris) in Royal Chitwan National Park, Nepal.
Dynamics of prey predator with Holling interactions and Smithsonian Contr. Zool. No. 336: 1-98.
stochastic influences. Alexandria Eng. J. 57: 1079–1086.
Thapa, K., Manandhar, S., Bista, M., Shakya, J., Sah,
Lamichhane, B.R., Pokharel, C.P., Poudel, S., Adhikari, G., Dhakal, M., Sharma, N., Llewellyn, B., Wultsch, C.,
D., Giri, S.R., Bhattarai, S., Bhatta, T.R., Pickles, R., Waits, L.P., Kelly, M.J., Hero, J., Huges, J., Karmacharya,
Amin, R., Acharya, K.P., Dhakal, M., Regmi, U.R., Ram, D. (2018). Assessment of genetic diversity, population
A.K., Subedi, N. (2017). Rapid recovery of tigers Panthera structure, and gene flow of tigers (Panthera tigris tigris)
tigris in Parsa Wildlife Reserve, Nepal. Oryx, 52(1): 16- across Nepal’s Terai Arc Landscape. PLoS ONE 13 (3):
24. 1-25.
MacKenzie, D.I. and Kendall, W.L. (2002). How should Thomas, L., Buckland, S.T., Rexstad, E.A., Laake, J.L.,
detection probability be incorporated into estimates of Strindberg, S., Hedley, S.L., Bishop, J.R.B. and Marques,
relative abundance. Ecology 83: 2387–2393. T.A. (2010). Distance software: design and analysis of
distance sampling surveys for estimating population size.
MacKenzie, D.I., Nichols, J.D., Lachman, G.B., Droege, Journal of Applied Ecology 47: 5-14.
S., Andrew R. J. and Langtimm, C.A. (2002). Estimating
site occupancy rates when detection probabilities are less White, G.C., Burnham, K.P.(1999). Program MARK:
than one. Ecology 83: 2248-2255. survival estimation from populations of marked animals.
Bird Study 46: S120–S139.
McDougal, C. (1977). The Face of the Tiger. Rivington
Books, London, UK. Wikramanayake, E.D., Dinerstein, E., Robinson, J.G.,
Karanth, U., Rabinowitz, A., Olson, D., Mathew, T.,
Ministry of Forests and Soil Conservation (MoFSC). Hedao, P., Conner, M., Hemley, G. and Bolze, D. (1998).
(2015). Strategy and Action Plan 2015-2025, Terai An ecology-based method for defining priorities for
Arc Landscape, Nepal. Ministry of Forests and Soil large mammal conservation: the tiger as case study.
Conservation, Singha Durbar, Kathmandu, Nepal. Conservation Biology 12 (4): 865-878.
STATUS OF TIGERS AND PREY IN NEPAL 2018 29
10 . ANNEXURE
Annex- 1. Summary of Spatially Explicit Capture Recapture (SECR) models for population estimation at site
level. Only the top models are presented.
PROTECTED
MODEL DETECTFN PAR LOGLIK AIC AICc AICc AICcWT g0 sigma
AREA
g0~bk
PNP hazard rate 5 -542.763 1095.526 1102.193 0 0.8517 0.02 3882.7
sigma~bk z~1
g0~bk
PNP hazard rate 4 -546.845 1101.689 1105.689 3.496 0.1483
sigma~1 z~1
g0~bk
CNP Exponential 3 -2565.33 5136.651 5136.947 0 1 0.03 2611
sigma~1
CNP g0~B sigma~1 Exponential 3 -2647.35 5300.708 5301.004 164.057 0
BaNP g0~B sigma~1 Exponential 3 -371.107 748.214 750.06 0 0.7197 0.03 2134.6
BaNP g0~B sigma~T Exponential 4 -370.999 749.999 753.332 3.272 0.1402
BaNP g0~T sigma~1 Exponential 3 -377.033 760.066 761.912 11.852 0
g0~bk
BNP Exponential 3 -2021.99 4049.987 4050.316 0 0.5544 0.04 2010.9
sigma~1
g0~bk
BNP Exponential 4 -2021.1 4050.197 4050.753 0.437 0.4456
sigma~T
BNP g0~B sigma~1 Exponential 3 -2091.47 4188.934 4189.262 138.946 0
g0~h2
ShNP sigma~1 Exponential 4 -512.691 1033.382 1037.382 0 0.7114 0.03 2369.2
pmix~h2
g0~h2
ShNP sigma~T Exponential 5 -511.26 1032.52 1039.186 1.804 0.2886
pmix~h2
g0: detection probability, Sigma:space range, T: time trend, B:transient response, bk: animal x site learned response, and h2: two-class mixtures
Annex- 2. Population abundance estimates with standard errors and 95% confidence intervals and detection probability estimates using program Capture and Mark
Annex- 3. Bayesian spatially explicit capture-recapture (SECR) analysis summary outputs from program SPACECAP. Density is presented as per 100 km2.
Annex- 4: Site-level tiger density estimates per 100 km2 based on the area surveyed in 2013 (for comparison).
Sites in bold show significant difference in density between 2013 and 2018.
Annex- 5. A comparative assessment of change in spatial density between 2013 and 2018 for each protected area.
Figure 1. Pixelated tiger density map of Parsa National Park and adjoining forests.
34 STATUS OF TIGERS AND PREY IN NEPAL 2018
Figure 2. Pixelated tiger density map of Chitwan National Park, adjoining forests and corridors.
STATUS OF TIGERS AND PREY IN NEPAL 2018 35
Figure 3. Pixelated tiger density map of Banke National Park, adjoining forests and corridors.
36 STATUS OF TIGERS AND PREY IN NEPAL 2018
Figure 4. Pixelated tiger density map of Bardia National Park, adjoining forests and corridors.
STATUS OF TIGERS AND PREY IN NEPAL 2018 37
Figure 5. Pixelated tiger density map of Shuklaphanta National Park, adjoining forests and
corridors.
38 STATUS OF TIGERS AND PREY IN NEPAL 2018
Annex- 6: A comparison of tiger density estimates derived by SECR-ML and SECR- Bayesian approaches.
(* denotes “and adjoining forest”) for effective sampling area of 2018 survey).
Annex-7: Prey density estimates and survey effort in each protected area and adjoining forests.
*Denotes the coverage of adjoining forests in line transect survey. ** denotes density estimates of all prey excluding swamp deer that was not detected during the line
transect survey in both Bardia National Park and Shuklaphanta National Park. (p denotes detection probability at defined area, ESW denotes Effective Strip Width of
detection, GOF- Chi P- denotes probability of chi square for goodness of fit test
STATUS OF TIGERS AND PREY IN NEPAL 2018 39
Annex- 8. Maps illustrating movement of tigers outside protected area boundary (the usage is illustrated by tiger
individuals having their MCP lying beyond the boundary)
Annex- 9. National Tiger Survey time frame and human resource involved.
No of No of
Person Elephant
SN Protected Area Training organized Field survey working personnel
days days
days involved
1 Chitwan NP, Parsa NP 28-29 November, 1 Dec, 2017- 3
63 411 25,893 1,575
and adjoining forests 2017 Feb 2018
2 Banke NP, Bardia NP 14-15 December, 17-Dec, 2017-
81 285 23,085 80
and adjoining forests 2018 11-March, 2018
3 Shuklaphanta NP, 23-March,
Laljhadi corridor and 3-4 February, 2018 2018-3-April, 20 81 1,620 80
Jogbuda forest 2018
4 1-Dec,
Habitat Occupancy Same date as to the
2017-3-April, 55 59 3,245
survey respective study sites
2018
Total 53,843 1,735
42 STATUS OF TIGERS AND PREY IN NEPAL 2018
PNP-FT01-RIGHT PNP-FT01-LEFT
PNP-FT02-RIGHT PNP-FT02-LEFT
PNP-FT07-RIGHT PNP-FT07-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 51
PNP-FT08-RIGHT PNP-FT08-LEFT
PNP-FT09-RIGHT PNP-FT09-LEFT
PNP-FT10-RIGHT PNP-FT10-LEFT
PNP-FT11-RIGHT PNP-FT11-LEFT
52 STATUS OF TIGERS AND PREY IN NEPAL 2018
PNP-FT12-RIGHT PNP-FT12-LEFT
PNP-FT15-RIGHT PNP-FT15-LEFT
PNP-FT16-RIGHT PNP-FT16-LEFT
PNP-MT01-RIGHT PNP-MT01-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 53
PNP-MT06-RIGHT PNP-MT06-LEFT
PNP-MT07-RIGHT PNP-MT07-LEFT
PNP-MT08-RIGHT PNP-MT08-LEFT
PNP-MT09-RIGHT PNP-MT09-LEFT
54 STATUS OF TIGERS AND PREY IN NEPAL 2018
CNP-FT41- LEFT
CNP-UT02- LEFT
CNP-UT04- LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 75
CNP-UT05 LEFT
BaNP-MT86-RIGHT BaNP-MT86-LEFT
76 STATUS OF TIGERS AND PREY IN NEPAL 2018
BaNP-FT87-RIGHT BaNP-FT87-LEFT
BaNP-FT88-RIGHT BaNP-FT88-LEFT
BaNP-MT92-RIGHT BaNP-MT92-LEFT
BaNP-FT94-RIGHT BaNP-FT94-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 77
BaNP-MT95-RIGHT BaNP-MT95-LEFT
BaNP-MT96-RIGHT BaNP-MT96-LEFT
BaNP-FT97-RIGHT BaNP-FT97-LEFT
BaNP-FT98-RIGHT BaNP-FT98-LEFT
78 STATUS OF TIGERS AND PREY IN NEPAL 2018
BaNP-FT99-RIGHT BaNP-FT99-LEFT
BaNP-FT100-RIGHT BaNP-FT100-LEFT
BaNP-MT101-RIGHT BaNP-MT101-LEFT
BaNP-FT115-RIGHT BaNP-FT115-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 79
BaNP-FT117-RIGHT BaNP-FT117-LEFT
BaNP-MT118-RIGHT BaNP-MT118-LEFT
BaNP-FT119-RIGHT BaNP-FT119-LEFT
BaNP-UT120-RIGHT BaNP-UT120-LEFT
80 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-MT01-RIGHT BNP-MT01-LEFT
BNP-MT02-RIGHT BNP-MT02-LEFT
BNP-FT03-RIGHT BNP-FT03-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 81
BNP-MT05-RIGHT BNP-MT05-LEFT
BNP-FT06-RIGHT BNP-FT06-LEFT
BNP-FT07-RIGHT BNP-FT07-LEFT
BNP-FT08-RIGHT BNP-FT08-LEFT
82 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-MT09-RIGHT BNP-MT09-LEFT
BNP-FT11-RIGHT BNP-FT11-LEFT
BNP-FT14-RIGHT BNP-FT14-LEFT
BNP-FT15-RIGHT BNP-FT15-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 83
BNP-MT16-RIGHT BNP-MT16-LEFT
BNP-FT17-RIGHT BNP-FT17-LEFT
BNP-FT18-RIGHT BNP-FT18-LEFT
BNP-FT19-RIGHT BNP-FT19-LEFT
84 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-MT20-RIGHT BNP-MT20-LEFT
BNP-FT21-RIGHT BNP-FT21-LEFT
BNP-MT22-RIGHT BNP-MT22-LEFT
BNP-FT23-RIGHT BNP-FT23-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 85
BNP-MT24-RIGHT BNP-MT24-LEFT
BNP-MT25-RIGHT BNP-MT25-LEFT
BNP-FT26-RIGHT BNP-FT26-LEFT
BNP-FT28-RIGHT BNP-FT28-LEFT
86 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-MT29-RIGHT BNP-MT29-LEFT
BNP-FT32-RIGHT BNP-FT32-LEFT
BNP-FT34-RIGHT BNP-FT34-LEFT
BNP-FT36-RIGHT BNP-FT36-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 87
BNP-FT37-RIGHT BNP-FT37-LEFT
BNP-MT38-RIGHT BNP-MT38-LEFT
BNP-FT39-RIGHT BNP-FT39-LEFT
BNP-FT40-RIGHT BNP-FT40-LEFT
88 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-MT41-RIGHT BNP-MT41-LEFT
BNP-MT42-RIGHT BNP-MT42-LEFT
BNP-FT43-RIGHT BNP-FT43-LEFT
BNP-MT44-RIGHT BNP-MT44-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 89
BNP-FT45-RIGHT BNP-FT45-LEFT
BNP-MT46-RIGHT BNP-MT46-LEFT
BNP-MT49-RIGHT BNP-MT49-LEFT
BNP-FT50-RIGHT BNP-FT50-LEFT
90 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-FT53-RIGHT BNP-FT53-LEFT
BNP-UT54-RIGHT BNP-UT54-LEFT
BNP-UT55-RIGHT BNP-UT55-LEFT
BNP-UT56-RIGHT BNP-UT56-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 91
BNP-UT57-LEFT
BNP-MT58-RIGHT BNP-MT58-LEFT
BNP-MT59-RIGHT BNP-MT59-LEFT
BNP-FT60-RIGHT BNP-FT60-LEFT
92 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-FT61-RIGHT BNP-FT61-LEFT
BNP-MT62-RIGHT BNP-MT62-LEFT
BNP-FT63-RIGHT BNP-FT63-LEFT
BNP-UT64-RIGHT BNP-UT64-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 93
BNP-UT65-RIGHT BNP-UT65-LEFT
BNP-FT66-RIGHT BNP-FT66-LEFT
BNP-MT67-LEFT
BNP-FT69-RIGHT BNP-FT69-LEFT
94 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-FT70-RIGHT BNP-FT70-LEFT
BNP-MT71-RIGHT BNP-MT71-LEFT
BNP-FT72-RIGHT BNP-FT72-LEFT
BNP-FT73-RIGHT BNP-FT73-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 95
BNP-MT74-RIGHT BNP-MT74-LEFT
BNP-FT75-RIGHT BNP-FT75-LEFT
BNP-FT76-RIGHT BNP-FT76-LEFT
BNP-MT77-RIGHT BNP-MT77-LEFT
96 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-FT78-RIGHT BNP-FT78-LEFT
BNP-MT82-RIGHT BNP-MT82-LEFT
BNP-FT84-RIGHT BNP-FT84-LEFT
BNP-MT85-RIGHT BNP-MT85-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 97
BNP-FT89-RIGHT BNP-FT89-LEFT
BNP-FT91-RIGHT BNP-FT91-LEFT
BNP-FT93-RIGHT BNP-FT93-LEFT
BNP-FT104-RIGHT BNP-FT104-LEFT
98 STATUS OF TIGERS AND PREY IN NEPAL 2018
BNP-MT105-RIGHT BNP-MT105-LEFT
BNP-MT106-RIGHT BNP-MT106-LEFT
BNP-FT109-RIGHT BNP-FT109-LEFT
BNP-FT110-RIGHT BNP-FT110-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 99
BNP-UT111-RIGHT BNP-UT111-LEFT
BNP-MT112-RIGHT BNP-MT112-LEFT
ShNP-FT02-RIGHT ShNP-FT02-LEFT
100 STATUS OF TIGERS AND PREY IN NEPAL 2018
ShNP-FT03-RIGHT ShNP-FT03-LEFT
ShNP-FT04-RIGHT ShNP-FT04-LEFT
ShNP-FT06-RIGHT ShNP-FT06-LEFT
ShNP-FT07-RIGHT ShNP-FT07-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 101
ShNP-FT08-RIGHT ShNP-FT08-LEFT
ShNP-MT01-RIGHT ShNP-MT01-LEFT
ShNP-MT02-RIGHT ShNP-MT02-LEFT
ShNP-MT03-RIGHT ShNP-MT03-LEFT
102 STATUS OF TIGERS AND PREY IN NEPAL 2018
ShNP-MT06-RIGHT ShNP-MT06-LEFT
ShNP-MT07-RIGHT ShNP-MT07-LEFT
ShNP-MT08-RIGHT ShNP-MT08-LEFT
ShNP-MT09-RIGHT ShNP-MT09-LEFT
STATUS OF TIGERS AND PREY IN NEPAL 2018 103
ShNP-MT10-RIGHT ShNP-MT10-LEFT
ShNP-MT11-RIGHT ShNP-MT11-LEFT
ACKNOWLEDGEMENT
The Government of Nepal acknowledges and appreciates the financial and technical support received
from the following organizations to carry out the Status of Tigers and Prey in Nepal, 2018:
STATUS OF TIGERS
AND PREY IN NEPAL
2018
S TAT U S O F T I G E R S A N D P R E Y I N N E PA L 2 0 1 8