Demography of a small, isolated tiger (Panthera tigris tigris) population in a semi-arid region of western India
BMC Zoology volume 2, Article number: 16 (2017)
Tiger populations have declined globally due to poaching, prey depletion, and habitat loss. The westernmost tiger population of Ranthambhore in India is typified by bottlenecks, small size, and isolation; problems that plague many large carnivore populations worldwide. Such populations are likely to have depressed demographic parameters and are vulnerable to extinction due to demographic and environmental stochasticity. We used a combination of techniques that included radio telemetry, camera traps, direct observations, and photo documentation to obtain 3492 observations on 97 individually known tigers in Ranthambhore between 2006 and 2014 to estimate demographic parameters. We estimated tiger density from systematic camera trap sampling using spatially explicit capture-recapture (SECR) framework and subsequently compared model inferred density with near actual density.
SECR tiger density was same as actual density and recovered from 4.6 (SE 1.19) to 7.5 (SE 1.25) tigers/100km2 over the years. Male: female ratio was 0.76 (SE 0.07), and cub: adult tigress ratio at 0.48 (SE 0.12). Average litter size was estimated at 2.24 (SE 0.14). Male recruitment from cub to sub-adult stage (77.8%, SE 2.2) was higher than that of females (62.5%, SE 2.4). But male recruitment rate as breeding adults from the sub-adult stage (72.6%, SE 2.0) was lower than females (86.7%, SE 1.3). Annual survival rates, estimated by known-fate models, of cubs (85.4%, CI95% 80.3–90.5%) were lower than that of juvenile (97.0%, CI95% 95.4–98.7%) and sub-adult (96.4%, CI95% 94.0–98.9%) tigers. Adult male (84.8%, CI95% 80.6–89.2%) and female (88.7%, CI95% 85.3–92.2%) annual survival rates were similar. Human-caused mortality was 47% in cubs and 38% in adults. Mean dispersal age was 33.9 months (SE 0.8), males dispersed further (61 Km, SE 2) than females (12 Km, SE 1.3). Higher age of first reproduction (54.5 months, SE 3.7) with longer inter-birth intervals (29.6 months, SE 3.15) was likely to be an effect of high tiger density.
Demographic parameters of Ranthambhore tigers were similar to other tiger populations. With no signs of inbreeding depression there seems to be no eminent need for genetic rescue. The best long-term conservation strategy would be to establish and manage a metapopulation in the Ranthambhore landscape.
At the onset of the nineteenth century, India was home to nearly 40,000 tigers (Panthera tigris tigris, Linnaeus) , while currently there are around 2200 left . The decline in tigers was primarily due to hunting, prey depletion, followed by habitat loss . A timely and proactive conservation measure, in the form of Project Tiger initiated in 1973 by the Indian Government , initially halted the rapid decline caused by trophy hunting. But an increased demand for tiger body parts in China and Southeast Asia in the past 25 years has severely impacted wild tiger populations. Demand driven poaching resulted in the local extinction of tigers in Sariska and Panna Tiger Reserves in India [5, 6]. Most tiger populations currently are small, isolated, and highly structured [3, 7, 8]. Such populations are vulnerable to extinction events caused by environmental and demographic stochasticity [9, 10]. Ranthambhore was a famous hunting reserve for the Maharajas of Jaipur, and numerous shikar (hunting) camps were organized in pre and post-independence era . Subsequent to India’s independence, intensity of hunting increased since tiger shikar was considered a social status symbol. This unregulated hunting caused a severe decline in Ranthambhore tiger population, and before the onset of Project Tiger (1973), there were around 14 tigers left in Ranthambhore . After an initial recovery in 1980’s, rampant poaching in 1992 and 2005, caused Ranthambhore tiger population to decline below 15 individuals from about 40 [13, 14]. Local extinctions in the last five decades suggest that tigers of the semi-arid region of western India are most vulnerable . The tiger population of Ranthambhore is the only population that survives in western India. It typifies the problems many large carnivore populations face globally i.e. small founder population and lack of connectivity with other source populations. Small isolated populations like Ranthambhore are susceptible to loss of genetic variability caused by genetic drift and inbreeding depression [16, 17]. Such populations often lose their ability to adept in response to environmental changes and some manifest deleterious effects in the form of morphological abnormalities and depressed population vigour [18, 19]. Hence, understanding demographic parameters of a potentially genetically compromised population to determine the need for genetic rescue is important for developing appropriate conservation strategies . Quantifying demographic parameters needs long-term data over multiple generations and for long-lived carnivores such datasets are rare . Till date most population studies conducted on tigers aim at estimating abundance [22,23,24,25,26,27,28,29], while studies on demographic parameters (like survival rates, litter size, sexual maturity, reproductive success) have been sparse (but see [30,31,32] for P. t. tigris and, [33, 34] for P. t. altaica).
Herein, we report demographic parameters of free ranging tigers from Ranthambhore Tiger Reserve from a nine-year study where 97 individually known tigers were monitored, and annual density estimated by spatially explicit capture-recapture using camera traps. By 2012 due to intensive monitoring, we had photo-captured almost all tigers of Ranthambhore and had developed a catalogue for individually identifying them. We use this information to compare snapshot density estimated by model-based inference with near reality. We also compare demographic parameters of Ranthambhore tigers with those of other tiger populations and conclude that though Ranthambhore tigers have undergone population bottlenecks, with limited gene flow and small population size, their demographic parameters do not seem to be compromised.
The study was conducted from 2006 to 2014 in Ranthambhore Tiger Reserve (hereafter RTR, latitudes 25041′ N to 26022′ N and longitudes 76016′ E to 77014′ E) which is situated at the junction of two ancient mountain ranges, the Aravalli and the Vindhya. RTR is part of the western Indian landscape that has Sariska Tiger Reserve in the north, Kuno Wildlife Sanctuary and Madhav National Park in the east, Ramgarh Visdhari Wildlife Sanctuary and Mukundara Hills Tiger Reserve in the south-western part (Fig. 1). The core area of RTR was composed of Ranthambhore National Park (392 km2), Sawai Mansingh Sanctuary (290 km2) while Kailadevi Wildlife Sanctuary (630 km2) was designated as the buffer zone of RTR. Within this western Indian landscape, tigers were only present in Ranthambhore NP during the commencement of this study and subsequently colonised Sawai Mansingh Sanctuary in 2008–09. These together comprise the only source population of tigers in the landscape. During this study, tigers from this population were reintroduced in Sariska  and six tigers dispersed into northern as well as south-eastern and eastern part of the landscape (Fig. 1).
The sub-tropical dry climate of RTR experiences three distinct seasons: mostly dry winters (October–February, minimum average temperature 5 °C, relative humidity ~10%), hot summers (March–June, mean maximum temperature 45 °C, relative humidity 10–15%), and humid monsoons (July–September, average rainfall 700 mm, relative humidity >60%). RTR primarily comprises of steep hills, gentle slopes, plateaus, and narrow valleys dotted with shallow man-made perennial lakes. The area is representative of dry deciduous Anogeissus pendula forests in association with Acacia, Butea, Capparis, Zizyphus and Prosopis species (5B/C2 - Northern Tropical Dry Deciduous forests. 6B/DS1 - Zizyphus scrub, DS1 - Dry deciduous scrub and 5/DS4 - Dry Grasslands, of Champion & Seth  classification). A diverse assemblage of carnivore species (17 species from 7 different families) were recorded during the course of the study, which include tiger, leopard (Panthera pardus, Linnaeus), sloth bear (Melursus ursinus, Shaw), striped hyena (Hyaena hyaena, Linnaeus), caracal (Caracal caracal, Schreber), fishing cat (Prionailurus viverrinus, Bennett), jungle cat (Felis chaus, Schreber), desert cat (Felis silvistris, Schreber), rusty-spotted cat (Prionailurus rubuginosa, I. Geoffroy Saint-Hilaire), golden jackal (Canis aureus, Linnaeus), Indian fox (Vulpes bengalensis, Shaw), honey badger (Mellivora capensis, Schreber), common palm civet (Paradoxurus hermaphorditus, Pallas), small Indian civet (Viverricula indica, É. Geoffroy Saint-Hilaire), Indian gray mongoose (Herpestes edwardsii, É. Geoffroy Saint-Hilaire), small Indian mongoose (Herpestes auropunctatus, Hodgson), and ruddy mongoose (Herpestes smithi, Gray). Tiger prey species present in the study area were chital (Axis axis, Erxleben), sambar (Rusa unicolor, Kerr), nilgai (Boselaphus tragocamelus, Pallas), chinkara (Gazella bennetti, Skyes), wild pig (Sus scrofa, Linnaeus), common langur (Semnopithecus entellus, Dufresne) and rhesus macaque (Macaca mulata, Zimmermann).
Field methods: Monitoring of tigers
We monitored 97 individual tigers during the study (2006 to 2014) through camera traps, radio-telemetry and routine patrolling (for direct sightings and photo-documentation) by researchers and forest staff. We developed criteria for classifying tigers into age groups by observing known-age individuals and use teeth eruption, ware, and body characteristics (see Additional file 1 for age estimation of tigers) similar to that of lions [37, 38]. We classified tigers into six age classes, namely, cubs (< 12 months), juveniles (12-24 months), sub-adults (2–3 years), young adults (4–5 years), prime adults (6–10 years) and old adults (> 10 years). Of the 97 tigers, 74 were known since cub stage, and their age was known to the exact month. Sex of the individuals was ascertained by the time cubs were 6 months old.
Camera traps were used a) in a systematic grid-based design (4km2 in 2006, 2009, 2012 and 2013, and 2km2 in 2014) for a short duration to estimate tiger abundance, and b) to target specific areas so as to determine the presence and use of the area by particular individual tigers throughout the study period.
a) The entire tiger occupied part of RTR was sampled in a systematic manner during the study by placing a pair of camera traps at each selected location within a grid. After conducting a reconnaissance survey, camera traps were placed on dirt roads, animal trails, fire lines and dry river beds at locations that maximized the chances of photo-capturing tigers. A pair of camera traps (TrailMaster® Lenexa KS USA, Cuddeback™ Green Bay USA, MOULTRIE® Alabama USA, or Stealth Cam® LLC Grand Prairie USA), facing each other, were placed at each location to get both flank photos of tigers at the same time. Each Camera was programmed with unique trap ID, time and date stamp on each photograph. Location of each camera trap was recorded by handheld GPS unit (Garmin 72™ and Garmin Etrex® 10, Kansas, USA) and plotted on a digitized map of RTR in GIS domain to ensure no sampling holes were present. Cameras were checked every 2–3 days to ensure proper functioning and recovery of data. Each photo captured tiger was identified to individuals by comparing their stripe patterns and given a unique id (e.g. T1, T2, and so on). The entire study area was sampled simultaneously within a period of 28 to 51 days so as to adhere to the assumption of population closure . The camera trapped area sampled each year as estimated by joining the outer most camera traps ranged between 139 to 492 km2 (Table 1).
b) Besides systematic camera trapping conducted once in a year, cameras were also used in a need based manner to record photographs of specific tigers. Areas of tigers that were not seen for over 60 days were specifically targeted. These special efforts were carried out by placing multiple camera traps and conducting an intensive ground search in most probable locations of that tiger. The effort continued till the fate of that tiger was ascertained by locating it or confirming its death or dispersal. Since the landscape outside of the tiger reserve was human-dominated, the presence of dispersing tigers was quickly detected by reports of livestock kills and sighting of tigers or their signs by villagers. The research team along with forest department staff subsequently tried to locate each such tiger and ascertain its identity through sighting, camera traps or hand held photography.
Eight tigers (three adult males, one adult female, three sub-adult males and one sub-adult female) were radio-collared between April 2007 to May 2009. Tigers were anesthetized with ketamine hydrochloride in combination with medetomidine injected intramuscularly using a gas-powered projectile dart delivery system . Tigers were collared with a Very High-Frequency transmitter (Telonics, Arizona, USA) and in most cases with a Global Positioning System with ground download facility (HABIT, British Columbia, Canada). Soon after the radio-collaring operation, Atipamezole was administered for reversing the effect of the tranquilizer . Animals were left after natural reflexes and behaviour returned and subsequently monitored through telemetry. Radio-collared tigers were tracked and regularly monitored (>three times a week) from a vehicle or on foot throughout the functional period of those collars. Tigers were tracked with the help of a hand-held directional 3- element Yagi antenna and Telonics and HABIT receivers. After the battery life of the radio-collars, the surviving collared tigers were monitored through camera traps and visual sightings.
Due to intensive camera trapping over the years and photo documentation by researchers and park officers, we were reasonably certain that almost all tigers of RTR were photo-captured by 2012 and were individually known. Our claim of knowing almost all RTR tigers was substantiated by the fact that in subsequent years no unknown adult tiger was recorded either by camera traps or by any other means. All additions to RTR population from 2012 onward till date were from known cubs. A photo album was developed and shared with the staff and officers of RTR in 2012 with additions and deletions done every 6 months. Since adult tiger numbers ranged between 20 to 40 individuals, it was possible to identify each tiger on most occasions it was sighted and often photographed with digital cameras by forest department staff and researchers. When in doubt these photographs were compared with the photo-album or on the computer to ascertain the tiger’s identity. A daily record on all tiger observations was maintained in a register which was subsequently transferred to a database.
Analytical methods: Estimating demographic parameters of tigers
We used likelihood based spatially explicit capture-recapture (SECR, [41, 42] in package ‘secr’ on R platform [43, 44] to estimate tiger density from the systematically sampled camera trap data over the years (2006, 2009, 2012, 2013, and 2014; Table 1). SECR consists of two sub-models. The distribution sub-model depicts the spatial distribution of detectors and animal captures in the landscape. The detection sub-model (g(x)) declines with increasing distance between the animal’s activity centre, and this spatial scale of detection is parameterized by sigma (σ). A spatial capture history matrix, a trap layout matrix, and a habitat mask that excluded non-habitat areas from the SECR model space were prepared and used in secr. Home-range size (as indexed by σ) is often correlated with density . Since the tiger population of RTR increased during our study period, we, parameterized σ and g0 (capture probability at the activity centre) separately for each year. Male and female tigers were likely to differ in their ranging patterns. Hence we used gender as a covariate to model heterogeneity in movement parameter (σ). Half-normal detection function was used to model σ. We used AICc (Akaike Information Criterion corrected for sample size)  to compare models with the null model (where g0 and σ were constant) and amongst themselves. Models with less than five delta AICc values were considered probable and the parameter estimates were obtained by AICc weighted model averages . Cubs (<12 months) were excluded from density estimation since this cohort is underrepresented in camera traps and has relatively high mortality [22, 48]. Since by 2012 we were reasonably certain that almost all tigers of RTR were individually known to us, we use this information to compare estimates obtained by SECR with actual known density.
Sex ratio and female reproductive parameters
Since sex of all tigers was known, we have calculated sex ratio (adult males: Adult females) for each year (sampling without replacement) by counting the total number of males and females in the population, and an estimate of its versatility between years .
Before giving birth, females restricted themselves to a small area of their territory [30, 50] and were often detected through frequent photo-captures and higher intensity of use at a particular site. Births were confirmed from photo-captures or direct sighting of lactating females (Additional file 1), while most cubs were recorded (photo-captured or sighted) with their mother only after they were about 2 months old. Since litter size at birth was rarely known, our reported litter size could be an underestimate as mortality before 2 months’ age was not known. The ratio of cub: adult tigress was computed. Age at first reproduction was determined by recording first birth of tigresses that were monitored since they were cubs. Inter-birth intervals were recorded from intensively monitored tigresses that littered more than once during the study period. We recorded intervals between two successive litters when all cubs of the previous litter died before reaching dispersal age (more than 2 years) and compare these with intervals between two successive litters when cubs of the previous litter survived beyond dispersal age.
Reproductive success and recruitment
We recorded reproductive success of tigresses as the number of cubs that survived to the age of independence (24 months). Individually identified tiger cubs that were monitored up to young adult stage allowed us to determine the age at which these tigers acquired territories either by displacing established tigers or in vacant habitats through dispersal. We calculated recruitment as the proportion of cubs that survived to the sub-adult stage (≥ 24 month). We also computed recruitment of sub- adults to successful breeders as the proportion of sub-adults that subsequently established territories. Territoriality was inferred when a young-adult exclusively used an area that was earlier used by other adult tigers of the same gender or dispersed to a vacant habitat and lived there for several months.
We estimated stage specific (cub, juvenile, sub-adult, young adult, prime adult, and old adult) survival probabilities of tigers by using known-fate model  in program MARK (ver.7, ). Known-fate models use Kaplan-Meier estimator [53, 54] to estimate survival which requires that fate of the individual is known with certainty during the period a particular individual is monitored. All individuals do not enter the study simultaneously but are added in a staggered manner known as the staggered entry design . Radio-telemetered individuals provide the best data for such analyses . However, in our case, besides the eight collared tigers that were located several times a week, each known individual tiger was observed (camera trap photo-captured, directly observed or photographed from hand-held cameras) at least once every month, and its fate recorded (see Additional file 2). In rare cases when an individual was not observed during an interval the known-fate model allows its observation to be ‘censored’ from the analysis . Data on survival/death of tigers was compiled on a monthly basis. We subsequently pooled this data for a three-month period to have a good number of observations on each tiger as well as to have an interval that was meaningful for survivorship analysis of tigers that are reasonably long-lived. Live-dead encounter history matrix for tigers was made by pooling encounter histories of 3 months into one interval where 10 represented survival of the individual throughout the interval, 11 represented mortality of the individual during the interval, and 00 represented censoring the individual when that individual was not observed during the three-month interval. Only 59 observations on seven tigers out of 3492 observations from 97 tigers were censored during the study (see Additional file 2), forming a very small proportion (1.6%) of observation where fate could not be ascertained for that interval. All of these seven tigers were observed in subsequent periods, but by censoring them for intervals they were not recorded, we add to the uncertainty in our estimates . Seven tigers (two juveniles, one sub-adult male, one young adult male, and three old adult females) out of 97 that we monitored went ‘missing’ during the study period. We considered two extreme scenarios for these individuals in our survival analysis: a) a highly likely conservative scenario, that these tigers were poached, and b) a more optimistic but less likely scenario that these tigers dispersed and we were unable to trace them. In case of scenario ‘a’, we modelled these tigers as dead, and in case of scenario ‘b’, we censored them from our analysis. The resultant estimates of survival would encompass the true estimate. For our analysis, the matrix was right censored for all the surviving tigers at the end of December 2014. Tigers that survived and were monitored from cub to adulthood were included in subsequent stages with the assumption that survival rates were independent for different stages. As part of tiger reintroduction program, seven tigers were translocated from RTR to Sariska Tiger Reserve during the study period . These individuals were censored in our survival analysis at the time of their translocation since these tigers would be exposed to a different set of environmental factors in Sariska that determine their survival probability. We formulated and run different candidate models which were ecologically plausible and compare them using AICc values. Model averaged survival estimates were obtained for models with less than five delta AICc values .
Mortality events were confirmed when carcasses were recovered, and causes of mortality were recorded on the basis of post-mortem report compiled by experienced veterinary personnel and/or by questioning eyewitnesses (if any). Cub mortalities were confirmed when found dead. However, not all cub carcasses could be recovered. A cub was considered dead if it was not detected (not photo-captured or sighted) with its mother for more than a month . There were no incidents where a cub that was not recorded for over a month (considered dead) was ever seen again. We considered seven tiger cubs (ranging from 2 to 10 months) that were provisioned by park managers after the death of their mother, as dead since they were unlikely to have survived in the wild on their own.
The small and isolated nature of RTR made it possible for us to follow and ascertain dispersal events of tigers. Tigers that left the National Park had to traverse human-dominated areas where their presence was detected readily by local communities who were always on a high vigil for large carnivores. The presence of tigers was detected from signs and examination of livestock kills (which were compensated by the Government). Our research team and park managers ascertained the identity of these tigers by targeted effort of camera trapping and visual sighting. Due to the high profile nature of RTR as well as the small tiger population that was vulnerable to poaching, the park management along with the research team made a concentrated effort to locate each tiger individually to ascertain its wellbeing. Mortality events were categorized into natural mortality (death of the mother, infanticide, old age, disease, intra-specific strife), human-caused mortality (poaching, poisoning, accidents due to human causes), and unknown.
Dispersal age and distances
We considered tigers to have reached dispersal age once they became ≥24 month old (considering the lowest age of dispersal). At this age most tigers no longer moved with their mothers and were capable of hunting on their own. Dispersals from the natal area were confirmed through telemetry, camera trap photographs, direct observations and hand-held photography of tigers. Dispersal distances (Euclidean distance) of these tigers were measured from the centre of their natal area to the most extreme location of that tiger.
On the average an effort of 3715 (SE 1324) camera trap nights were invested each year (Table 1). In most years the model having g0 as constant (.) and movement parameter (σ) having sex-specific responses was selected as the best-fit model. In year 2009 the null model (g0(.), σ(.)) was selected as the best model (Additional file 3: Table S1). Annual model averaged density estimates varied from 4.6 to 8.7 tigers per 100 km2 (Table 1). The movement parameter (σ) was consistently larger for males compard to females and was found to decline asymptotically with an increase in density for males (Fig. 2). The density estimate obtained from SECR did not differ from near actual density (Table 1).
Sex ratio and female reproductive parameters
Adult sex ratio (male: female) in the initial years (2006–07) was female biased (0.38, SE 0.04) and became marginally female biased (0.91, SE 0.04, 2008–14) in subsequent years. Overall male: female ratio was 0.76 (SE 0.07) during the study period. The total number of breeding females in RTR ranged between 12 to 15 (see Additional file 4: Figure S1). We recorded litter size from 33 litters of 18 females, and the mean litter size was 2.24 (SE 0.14; range = 1–4 cubs). Most of the litters were of two (50%) or three (31%) cubs (Fig. 3). Overall cub: adult female ratio (n = 9 years) was 0.48 (0.12 SE). We could ascertain the age at first reproduction for 11 tigresses that were monitored since they were cubs as 54.5 months (3.7 SE; range 33–68 months). Eight out of these 11 tigresses produced their first litter after 4 years of age, while only two tigresses gave birth before 3 years (see Additional file 5: Figure S2). The average interval between two successive litters (inter-birth interval) was 29.6 months (SE 3.15; range = 7–51 months, n = 14 intervals from 8 tigresses). Intervals were shorter when all cubs of previous litter died before reaching independence (24 months, 15.0, SE 4.04 months, n = 3) than when cubs of previous litters survived till independence (33.64, SE 2.83 months, n = 11, see Additional file 6: Figure S3).
Reproductive success and recruitment
More than 50% of the intensively monitored breeding females (n = 18) successfully raised all their litters to independence, while around 10% of the females failed to raise any of their cubs during the study period (Fig. 4). Data on 51 individuals indicated that male recruitment rate from cub to the sub-adult stage was higher (77.8%, SE 2.2) than females (62.5%, SE 2.4). But male recruitment rate as breeding adults in the population from the sub-adult stage was lower (72.6%, SE 2.0) than females (86.7%, SE 1.3, Table 2).
The average annual survival rate of cubs (85.4%) were comparatively lower than that of other age classes (Table 3). Annual survival rate increased after cub stage, remained constant till prime adult stage after which it declined (Table 3). Survival rates for males and females were similar for younger stages while females had marginally higher survival in older stages (Table 3). The survival estimates with our conservative approach of considering seven missing tigers as dead did not statistically differ from survival estimates when these tigers were censored from the analysis (Additional file 7: Table S2).
We recorded 25 mortality events (17cubs and juveniles, and 8 adults) during our study. Amongst all cub and juvenile mortality, 41% were natural (infanticide and death of mother due to natural causes), 47% were human caused (mother poached, accidents due to human causes and poisoning) and 12% of these could not be ascertained (Fig. 5). Amongst all adult mortality, 50% were natural (old age, disease, and intra-specific strife), and 38% were human-caused (poisoning and poaching, Fig. 5), the cause of 12% adult mortality could not be determined.
Dispersal age and distances
We recorded dispersal of 29 tigers, of these six were long distance dispersal out of RTR (Fig. 1). Mean dispersal age of tigers in RTR was 33.9 months (SE 0.8, range = 24–42 months). Mean male dispersal distances (60.6, SE 2.1 Km) were larger than that of females (12.1, SE 1.3 Km, Fig. 6). Most of the females established their territory near their natal area while the majority of the males dispersed further from their natal areas. Five out of 16 males and one out of 13 females dispersed outside RTR and settled in forest patches within the larger landscape. These long distance dispersal movements ranged from 56 to 220 km Euclidean distance from their natal areas.
Tigers are conservation dependent species and require substantial investments in terms of management and protection for their long-term persistence . Information on survival, recruitment, litter size, reproductive parameters, dispersal, sex ratio and density that we provide in this paper are the basis of estimating population viability and planning management interventions.
We mostly use standard methodology, but adept some and develop a few approaches required for procuring data from endangered carnivores. The subsequent analyses of these data do not violate any underlying assumptions of the analytical procedures. We believe that we were justified in using Known-Fate model for estimating survival due to the high frequency (every month) of observations of 97 individually known tigers for determining their fate (see Additional file 2). As explained in our methods, we ‘censored’ surviving individuals when the study was completed and for the few occasions in-between when we could not determine the fate of individuals . Open population capture-mark-recapture models cannot distinguish between emigration and mortality, and are confounded by the nuisance parameter of imperfect detection [51, 57]. Therefore, known-fate models though being data intensive, provide more precise and more informative parameter estimates.
During the initial phase of the study (i.e. 2006) tiger density in RTR was low (Table 1) as the population was recovering from a recent decline caused by poaching. In subsequent years, density increased with good protection and fluctuated between 5.6 to 8.7 tigers/ 100km2 with a mean density of 7.5 (SE 2.7) tigers/ 100km2 (Table 1). This fluctuation in tiger density was likely due to synchrony in breeding by several females and recruitment of a large cohort of sub-adults that became available for camera trap sampling at an approximate interval of 2 years. As tiger density increases within a limited area, we would expect home-range to either decrease and/or show an increase in overlap. The movement parameter σ is an index of home-range radius over short time duration . Efford et al.  show an asymptotic decline in σ with increase in tiger density from data across India. Herein, we demonstrate a similar relationship between σ and tiger density within a single population which we believe is ecologically more meaningful (Fig. 2). The asymptotic nature of σ for males and a relatively constant σ for females suggests that male tigers’ home range decline to some extent with increase in tiger density, while home range size in females, which is based primarily on food availability , does not change with density. This could be interpreted to suggest that either tigresses do not invest energy in acquiring a home range larger than required to rear cubs or that Ranthambhore National Park was already near carrying capacity density with little scope of reduction in home-range sizes for tigresses. The increase in tiger numbers in RTR was accompanied with an increase in tiger occupancy while the number of breeding females remained relatively constant (12 to 15) during our study period (see Additional file 4: Figure S1). Camera trap data suggests that the entire core area of RTR (Ranthambhore National Park and Sawai Mansingh Sanctuary) was occupied by tigers by 2014. We also witnessed intense competition between mothers and daughters for breeding territories. This was another indication that tiger density was at or near carrying capacity within the core of RTR. For high density populations near carrying capacity regulatory mechanisms may include delayed age of first reproduction, increased inter-birth intervals, smaller litter sizes and depressed survival . Delayed female age of first reproduction and longer intervals between two successive litters observed in RTR tigers were likely an effect of high tiger density at or nearing carrying capacity inside the core of RTR. Parameters that could potentially be depressed by inbreeding depression like litter size, cub survival, and disease caused mortality were either similar or better in comparison to other tiger populations (Table 4).
In highly inbred populations morphological abnormalities are often observed and resistance to disease is often compromised [18, 59]. During our study we did not encounter cubs with abnormalities or skeletal defects and mortality of tigers attributed to disease was only one. The morphometric measurement of RTR tigers captured for radio-collaring were among those of the largest recorded for tigers in India (YV Jhala, unpublished data). These observations along with comparable demographic parameters to other tiger populations, suggest that there were no deleterious effects that had as yet manifested among RTR tigers due to population bottlenecks, small size and isolation . This seems likely, since, till recent times the semiarid zone tiger population of central India was large and well connected  and the population bottlenecks that RTR tigers passed through did not remain very small for long periods of time (less than one generation time) . However, we lack data on early infant mortality and foetal loss during pregnancy, which are some of the parameters that would be influenced by inbreeding depression. Reddy et al.  conclude that RTR tigers had reasonable genetic diversity comparable to other central Indian tiger populations. Studies that link genetic variability with population size, connectivity, and ultimately with demography are required for conservation management of large carnivores.
Females that had lost cubs before their recruitment age gave birth within a smaller time interval as also reported in tigers  and in Asiatic lions . Due to a high level of protection, poaching was rare inside RTR during the study period. Therefore, long and stable tenures of territorial males were recorded (mean 5.6 years, range 4 to 11 years, n = 13). Only three cases of infanticide were observed during this study; this contrasts with studies on tigers in Chitwan National Park in Nepal  and Asiatic lions in Gir , where infanticide was a major cause of cub mortality. Mean dispersal age of RTR tigers (~33 months) was higher than reported for Chitwan (~23 months) , and Amur tigers (~19 months, Kerley et al. 2003). Young tigers were not compelled to leave their natal territories due to new male takeovers . This stability in territorial male tenures allowed young tigers, especially males, to continue to live longer within their natal area enhancing their survival and achieving rapid growth to breeding body size. On four occasions we observed sub-adult male tigers occupy small ranges in peripheral areas of resident male territories. Once such tigers become sufficiently large so as to challenge resident males, they expand their range into that of the residential males’ territory. By following this strategy relatively smaller and inexperienced sub-adult males avoid lethal encounters with larger and experienced males during their initial dispersal stage .
Ten out of 13 sub-adult females established territories near their natal areas, two females established their territory inside RTR but ~20 km away from their natal area, and one female dispersed outside the reserve (~90 km). We observed that females either established their territory beside their mothers (n = 3) or occupied a part of their mother’s territory and gradually pushed their mother off (n = 5). One female that initially occupied her mother’s territory by displacing her mother produced her first litter in this territory, but subsequently shifted her territory with her 12-month-old cubs (of her second litter) ~10 km away from her natal territory. This shift coincided with the takeover of her natal territory by a new male tiger. Her natal territory was subsequently occupied by her daughter from her first litter.
The dispersing sex in tigers is known to be males . Amongst RTR tigers, males did disperse larger distances than females, but opportunities to disperse were restricted in this landscape. Large dispersal distances travelled by tigers (once out of RTR), and older age of dispersal reflects the difficulty a tiger faces in locating appropriate vacant habitat to settle. The small reserve size combined with very little and disjunct tiger habitat available outside the reserve system within a hostile human-dominated habitat matrix restricted dispersal. Males that managed to disperse and locate habitat patches tried to settle there, but due to lack of female tigers in these patches, were unable to breed. If areas of Kailadevi Wildlife Sanctuary (part of RTR buffer), Kuno Wildlife Sanctuary, Ramgarh Visdhari Wildlife Sanctuary and Mukundara Hills Tiger Reserve (Fig. 1) are made free of human settlements through incentivised voluntary relocation scheme (WPA 1972; 2006 amendment), these areas could harbour breeding population of tigers. Simultaneously dispersal corridors between these tiger habitats and with core of RTR need to be secured and restored to promote a metapopulation in this landscape. Initially, once these areas have been appropriately restored and have sufficient prey base, this tiger dispersal could be aided by translocating tigers. Establishing and managing the tiger population in the larger Ranthambhore landscape as a metapopulation  would be desirable for long-term conservation.
Park authorities often intervene by treating injured tigers in-situ or supplementing food resources for orphan/sick cubs. These activities, although well intentioned, hinder the natural process of selection and the social dynamics of the species [38, 64]. Such interventions also reduce the genetic fitness of the population over time by ensuring the survival of unfit individuals especially if sick animals are treated, saved and allowed to breed . Interfering with the natural process may be necessary for highly endangered populations where every living individual counts. However, RTR is now a high tiger density area , and therefore, management interference of health care should be extremely selective, if any.
Our study did not find any evidence of detrimental effects resulting from a small population that could potentially be inbred and there seems to be no current need for genetic rescue  of RTR tigers. Currently RTR has about 15 breeding units and adult female survival of about 88%, these are bare minimal requirements to ensure long-term persistence . Managing the Ranthambhore landscape by restoring habitat patches and connectivity with RTR so as to promote a metapopulation of tigers would enhance the potential of long-term persistence of this last remaining semi-arid tiger population in western India.
- AICc :
Akaike Information Criterion, corrected for sample size
- CI95% :
95% Confidence Interval
Capture probability at activity centre
Ranthambhore Tiger Reserve
Spatially Explicit Capture Recapture
Sigma, movement parameter
Gee EP. The wildlife of India. London: UK. Collins; 1964.
Jhala YV, Qureshi Q, Gopal R. The status of tigers in India 2014. New Delh and Dehradun: National Tiger Conservation Authority and The Wildlife Institute of India; 2015.
Dinerstein E, Loucks C, Wikramanayake E, Ginsberg J, Sanderson E, Seidensticker J, Forrest J, Bryja G, Heydlauff A, Klenzendorf S, Leimgruber P. The fate of wild tigers. Bioscience. 2007;57:508–14.
Panwar HS. Project Tiger: the reserves, the tigers and their future. In: Tilson RL, Seal US, editors. Tigers of the world: the biology, biopolitics, management and conservation of an endangered species. Park Ridge: Noyes Publications; 1987. p. 110–7.
Check E. The tiger’s retreat. Nature. 2006;441:927–30.
Gopal R, Qureshi Q, Bhardwaj M, Singh RKJ, Jhala YV. Evaluating the status of the endangered tiger Panthera Tigris and its prey in Panna Tiger Reserve, Madhya Pradesh, India. Fauna and Flora International, Oryx. 2010;44:383–9.
Ranganathan J, Chan KMA, Karanth KU, Smith JLD. Where can tiger persist in the future? A landscape-scale, density-based population model for the Indian subcontinent. Biol Conserv. 2008;141:67–77.
Yumnam B, Jhala YV, Qureshi Q, Maldonado JE, Gopal R, Saini S, Srinivas Y, Fleischer RC. Prioritizing Tiger conservation through landscape genetics and habitat linkages. PLoS One. 2014;9(11):e111207.
Caughley GC. Directions in conservation biology. J Anim Ecol. 1994;63:215–44.
Purvis A, Gittleman JL, Cowlishaw G, Mace GM. Predicting extinction risk in declining species. Proc Biol Sci. 2000;267:1947–52.
Singh K. Shikar Camps. In: Rangarajan M, editor. The Oxford anthology of Indian wildlife volume I: hunting and shooting. New Delhi: Oxford University Press; 1999. p. 35–50.
IBWL (Indian Board for Wild Life). Project Tiger, a planning proposal for the preservation of the Tiger (Panthera tigris tigris Linn.) in India. Dehradun: F.R.I. Press; 1972.
Jackson P. Fifty years in the tiger world: an introduction. In: Tilson RL, Seal US, editors. Tigers of the world: the biology, politics, management and conservation of an endangered species. Park Ridge: Noyes Publications; 2010. p. 1–15.
Sharma S, Wright B. Monitoring tigers in Ranthambhore using digital pugmark technique: Technical Report. New Delhi: Wildlife Protection Society of India; 2005.
Chundawat RS, Sharma K, Gogate N, Malik PK, Vanak AT. Size matters: scale mismatch between space use patterns of tigers and protected area size in a tropical dry Forest. Biol Conserv. 2016;197:146–53.
Allendorf FW, Leary RF. Heterozygosity and fitness in natural populations of animals. In: Conservation biology: the science of scarcity and diversity. M. E. Smith, editor, Sinauer, Sunderland, MA; 1986. p. 57–76.
Frankham R, Briscoe DA, Ballou JD. Introduction to conservation genetics. Cambridge: Cambridge University Press; 2002.
Keller LK, Waller DM. Inbreeding effects in wild populations. Trends Ecol Evol. 2002;17:230–41.
O'brien SJ. Tears of the cheetah: the genetic secrets of our animal ancestors. St Martin's Press New York: Thomas Dunne Books; 2003.
Pimm SL, Dollar L, Bass OL. The genetic rescue of the Florida Panther. Anim Conserv. 2006;9:115–22.
Balme GA, Batchelor A, Britz NDEW, Seymour G, Grover M, Hes L, Macdonald DW, Hunter LTB. Reproductive success of female leopards Panthera pardus: the importance of top-down processes. Mammal Rev. 2012;43:221–37.
Karanth KU. Estimating tiger (Panthera tigris) populations from camera-trap data using capture-recapture models. Biol Conserv. 1995;71(3):333–8.
Karanth KU, Nichols JD. Estimation of tiger densities in India using photographic captures and recaptures. Ecology. 1998;79:2852–62.
Karanth KU, Nichols JD, Kumar NS, Link WA, Hines JE. Tigers and their prey: predicting carnivore densities from prey abundance. PNAS. 2004;101:4854–8.
Kawanishi K, Sunquist ME. Conservation status of tigers in a primary rainforest of peninsular Malaysia. Biol Conserv. 2004;120:329–44.
Barlow ACD, Ahmed MIU, Rahman MM, Howlader A, Smith AC, Smith JLD. Linking monitoring and intervention for improved management of tigers in the Sundarbans of Bangladesh. Biol Conserv. 2008;141:2031–40.
Jhala YV, Gopal R, Qureshi Q. Status of tigers, co-predators, and prey in India. New Delhi and Dehradun: National Tiger Conservation Authority, Govt. of India and Wildlife Institute of India; 2008.
Jhala YV, Qureshi Q, Gopal R. Can the abundance of tigers be assessed from their signs? J Appl Ecol. 2011;48:14–24.
Harihar A, Prasad DL, Ri C, Pandav B, Goyal SP. Losing ground: tigers Panthera Tigris in the north-western Shivalik landscape of India. Oryx. 2009;43:35–43.
Smith JLD, McDougal C. The contribution of variance in lifetime reproduction to effective population size in tigers. Conserv Biol. 1991;5:484–90.
Smith JLD. The role of dispersal in structuring the Chitwan tiger population. Behaviour. 1993;124:165–95.
Singh R, Mazumdar A, Sankar K, Qureshi Q, Goyal SP, Nigam P. Interbirth interval and litter size of free-ranging Bengal tiger (Panthera tigris tigris) in dry tropical deciduous forests of India. Eur J Wildlife Res. 2013;59:629–36.
Kerley LL, Goodrich JM, Miquelle DG, Smirnov EN, Quigley H, Hornocker MG. Reproductive parameters of wild female Amur (Siberian) tigers (Panthera tigris altaica). J Mammal. 2003;84:288–98.
Goodrich JM, Kerley LL, Smirnov EN, Miquelle DG, McDonald L, Quigley HB, Hornocker MG, McDonald T. Survival rates and causes of mortality of Amur tigers on and near the Sikhote-Alin biosphere Zapovednik. J Zool. 2008;276:323–9.
Sankar K, Goyal SP, Qureshi Q. Assessment of status of tiger (Panthera tigris) in Sariska Tiger Reserve, Rajasthan. Dehra Dun: A Report submitted to the Project Tiger, Ministry of Environment & Forests, Govt. of India, New Delhi Wildlife Institute of India; 2005. p. 1–26.
Champion HG, Seth SK. A revised survey of the Forest types of India. Dehli: Manager of Publications; 1968.
Schaller GB. The Serengeti lion. Chicago: University of Chicago Press; 1972.
Banerjee K, Jhala YV. Demographic parameters of endangered Asiatic lions (Panthera leo persica) in Gir forests. India Journal of Mammalogy. 2012;93:1420–30.
Chao A, Huggins RM. Classical closed-population capture-recapture models. In: Amstrup SC, TL MD, Manly BF, editors. Handbook of capture-recapture analysis. Princeton: Princeton University Press; 2010. p. 22–35.
Kreeger TJ. Handbook of wildlife chemical immobilization. 1st ed. Laramie: International Wildlife Veterinary Services Inc; 1996. p. 175.
Efford MG. Estimation of population density by spatially explicit capture–recapture analysis of data from area searches. Ecology. 2011;92:2202–7.
Borchers DL, Efford MG. Spatially explicit maximum likelihood methods for capture–recapture studies. Biometrics. 2008;64:377–85.
Efford MG. Secr: spatially explicit capture-recapture models. R package version 2.9.3. 2015. http://CRAN.R-project.org/package=secr.
R Core Team. R: a language and environment for statistical computing. Vienna: R foundation for Statistical Computing. http://www.R-project.org/. Accessed 15 Mar 2017.
Efford MG, Dawson DK, Jhala YV, Qureshi Q. Density-dependent home-range size revealed by spatially explicit capture–recapture. Ecography. 2015;39:676–88.
Akaike H. A new look at the statistical model identification. IEEE Trans Autom Control. 1974;19:716–23.
Burnham KP, Anderson DR. Model selection and multimodel inference- a practical information-theoretic approach. 2nd ed. New York: Springer; 2002.
Sharma RK, Jhala Y, Qureshi Q, Vattakaven J, Gopal R, Nayak K. Evaluating capture–recapture population and density estimation of tigers in a population with known parameters. Anim Conserv. 2010;13(1):94–103.
Skalski JR, Ryding KE, Millspaugh JJ. Wildlife demography: analysis of sex, age, and count data. Burlington: Elsevier Academic Press; 2005.
Sunquist ME. The social organization of tigers (Panthera tigris) in Royal Chitawan National Park, Nepal. Washington: Smithsonian Institution Press; 1981.
Williams BK, Nichols JD, Conroy MJ. Analysis and management of animal populations. Oxford: Academic Press; 2002.
Cooch W, White GC. Program MARK: a gentle introduction. 17th ed; 2017. p. 16-1–16-25. http://www.phidot.org/software/mark/docs/book/.
Kaplan EL, Meier P. Non-parametric estimation from incomplete observations. J Am Stat Assoc. 1958;53:457–81.
Pollock KH, Winterstein SR, Bunck CM, Curtis PD. Survival analysis in telemetry studies: the staggered entry design. J Wildlife Manage. 1989;53:7–15.
Sankar K, Nigam P, Malik PK, Qureshi Q, Bhattarcharjee S. Monitoring of reintroduced tigers (Panthera tigris tigris) in Sariska Tiger Reserve, Rajasthan. Technical report −1. Dehradun: Wildlife Institute of India; 2013.
Sanderson EW, Forrest J, Loucks C, Ginsberg J, Dinerstein E, Seidensticker J, Leimgruber P, Songer M, Heydlauff A, O’Brien T, Bryja G. Setting priorities for tiger conservation: 2005–2015. In: Tilson RL, Seal US, editors. Tigers of the world: the science, politics, and conservation of Panthera tigris. Park Ridge, New Jersey: Noyes Publications; 2010. p. 143–61.
Nichols JD. Modern open-population capture-recapture models. In: Amstrup SC, TL MD, Manly BF, editors. Handbook of capture-recapture analysis. Princeton: Princeton University Press; 2010. p. 88–123.
Derocher AE, Stirling I. The population dynamics of polar bears in western Hudson Bay. In: McCullough D, Barrett R, editors. Wildlife 2001: populations. England: H. Elsevier Science Publishers Ltd; 1992. p. 1150–9.
Räikkönen J, Vucetich JA, Peterson RO, Nelson MP. Congenital bone deformities and the inbred wolves (Canis Lupus) of isle Royale. Biol Conserv. 2009;142:1025–31.
Reddy GV, Tyagi RK, Bhatnagar D, Soni RG, Daima ML, Sen A. Management plan of Ranthambhore Tiger Reserve (2002–2003 to 2011–2012). Rajasthan, India: Forest Department; 2002. p. 363–6.
Reddy PA, Gour DS, Bhavanishankar M, Jaggi K, Hussain SM, Harika K, Shivaji S. Genetic evidence of Tiger population structure and migration within an isolated and fragmented landscape in Northwest India. PLoS One. 2012;7:e 29827.
Karanth KU. Tiger ecology and conservation in the Indian subcontinent. J Bombay Nat Hist Soc. 2003;100:169–89.
Hanski I. Metapopulation dynamics. Nature. 1998;396:41–9.
Packer C, Brink H, Kissui BM, Maliti H, Kushnir H, Caro T. Effects of trophy hunting on lion and leopard populations in Tanzania. Conserv Biol. 2011;25(1):142–53.
Sadhu A, Gupta D, Latafat K, George S, Jhala YV, Qureshi Q. Ranthambhore Tiger Reserve. In: Jhala YV, Qureshi Q, Gopal R, editors. The status of tigers in India 2014. New Delhi and Dehradun: National Tiger Conservation Authority and The Wildlife Institute of India; 2015. p. 167–9.
Chapron G, Miquelle DG, Lambert A, Goodrich JM, Legendre S, Clobert J. The impact on tigers of poaching versus prey depletion. J Appl Ecol. 2008;45:1667–74.
Sankhala K. Tiger! The story of the Indian tiger. London: Collins; 1978.
Chundawat RS, Gogate N, Malik PK. Understanding tiger ecology in the tropical dry deciduous forests of Panna Tiger Reserve. Final report. Dehradun: Wildlife Institute of India; 2002.
Rajesh Gopal is specially thanked for encouragement and facilitation. Director, Dean and Research Coordinator of the Wildlife Institute of India are acknowledged for logistic support. We thank G. V. Reddy, Rajesh Gupta, and Y. K. Sahu for support. We thank all frontline staff of Ranthambhore Tiger Reserve and our field assistants Ram Prasad, Mujahid, Kannaiya, and Javed for their sincere effort. S. Dutta, U. Kumar, and K. Banerjee are acknowledged for assistance with data analysis. G. S. Bhardwaj and D. S. Shaktawat are acknowledged for providing photographs.
This work was funded by the National Tiger Conservation Authority for the research project “Monitoring Source Population of Tiger in Ranthambhore Tiger Reserve”.
Availability of data and materials
All data are either reported in tables in the text or in the supplementary material online.
Permissions for capture and collaring tigers were obtained under the Wildlife (Protection) Act 1972 from the Ministry of Environment and Forests, Government of India and the Chief Wildlife Warden, Government of Rajasthan. The technical committee of the National Tiger Conservation Authority which also considers the well-being of animals and ethics of research approved the research project. Tigers were anesthetized using standard drugs under supervision by qualified veterinarians. All tigers were observed from a safe distance till they fully recovered from the anaesthesia and walked away into the forest.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Field Guide for Aging Tigers. (PDF 8068 kb)
Individual tiger monitoring data used for “Known Fate” survival analysis. (PDF 98 kb)
Sampling details and parameters estimates of annual tiger density from camera trap based spatial capture-recapture analysis in Ranthambhore Tiger Reserve. (PDF 84 kb)
Number of breeding tigresses in each year observed in Ranthambhore Tiger Reserve during the study period. (TIFF 258 kb)
Age at first reproduction of tigresses (n = 11) observed in Ranthambhore Tiger Reserve. (TIFF 515 kb)
Inter-birth intervals of tigresses in Ranthambhore Tiger Reserve when all cubs of previous litter died before reaching independence (n = 3 litters) and when cubs of previous litters survived beyond 24 months (n = 11). (TIFF 344 kb)
Survival rates of tigers (n = 97) estimated in an optimistic (where we have censored seven ‘missing’ tigers) as well as a conservative (a more likely scenario where we have considered seven ‘missing’ tigers as dead) approach in Ranthambhore between 2006 to 2014. (PDF 8 kb)
About this article
Cite this article
Sadhu, A., Jayam, P.P.C., Qureshi, Q. et al. Demography of a small, isolated tiger (Panthera tigris tigris) population in a semi-arid region of western India. BMC Zool 2, 16 (2017). https://doi.org/10.1186/s40850-017-0025-y
- Camera traps
- Inter-birth interval
- Known fate
- Litter size
- Spatially explicit capture-recapture