Introduction
The tiger Panthera tigris occurs in a wide range of climates and habitats, from the tropical forests of South Asia to the temperate forests of Far East Russia (Sunquist et al., Reference Sunquist, Karanth, Sunquist, Seidensticker and Jackson1999). However, the species is highly threatened across its entire range by poaching, habitat loss and prey depletion (Kenny et al., Reference Kenny, Smith, Starfield and McDougal1995; Wikramanayake et al., Reference Wikramanayake, Dinerstein, Robinson, Karanth, Rabinowitz and Olson1998; Karanth et al., Reference Karanth, Nichols, Kumar, Link and Hines2004). Its distribution is greatly restricted; populations that are viable, at least in the short term, occur only in Bangladesh, Bhutan, India, Indonesia, Malaysia, Nepal, Russia, and Thailand (Seidensticker, Reference Seidensticker2010; Walston et al., Reference Walston, Robinson, Bennett, Breitenmoser, Fonseca and Goodrich2010). The species is extinct in Bali, Java, southern China and central Asia, and populations in other range countries are now severely reduced and at imminent risk of extinction because of their small size (Carroll & Miquelle, Reference Carroll and Miquelle2006).
In response to these significant declines wildlife managers are attempting a number of strategies, alone or in combination, to increase the number of tigers. Efforts include restoring tiger habitat, re-establishing habitat connectivity, restoring prey populations, reducing poaching of tigers and their prey, and reintroducing animals to the wild (Seidensticker, Reference Seidensticker2010). All of these management options require an understanding of the relationship of tigers to their prey, and of the influence of ecological variables on both tigers and their prey. Although ecological variables can affect tigers directly, many are likely to have a pronounced influence on the density of the tiger's prey. Current models suggest that, in the absence of heavy poaching, the density of tigers is primarily determined by prey abundance (Karanth et al., Reference Karanth, Nichols, Kumar, Link and Hines2004; Wegge et al., Reference Wegge, Odden, Pokharel and Storaas2009; Harihar et al., Reference Harihar, Pandav and Gopal2009). Documentation of a positive relationship between tiger populations and their prey does not, however, explain the behavioural mechanisms underlying population growth or resilience. As reported for many solitary felids, breeding female tigers in Nepal and Russia occupy defended home ranges (Seidensticker, Reference Seidensticker1976; Sunquist, Reference Sunquist1981; Smith et al., Reference Smith, McDougal, Sunquist, Tilson and Seal1987; Goodrich et al., Reference Goodrich, Miquelle, Smirnov, Kerley, Quigley and Hornocker2010) and individual females have the sole responsibility for raising their young, using prey they acquire within their home range. Although prey abundance varies widely between these two countries, females in both localities are usually able to establish breeding territories large enough to raise young. Because previous studies have demonstrated that female home range size is a function of prey density (Smith et al., Reference Smith, McDougal, Sunquist, Tilson and Seal1987; Miquelle et al., Reference Miquelle, Goodrich, Smirnov, Stephens, Zaumyslova and Chapron2010) we hypothesize that this relationship holds across the tiger's range and that mean female home range size may be used to calculate local carrying capacity for tigers.
The number of female territories in a population determines that population's recruitment potential and ultimately its viability and resilience. Thus, a measurement of the number of breeding female tigers that can exist in a given area (i.e. carrying capacity) is a useful metric for managers. This measurement would be particularly useful if tiger carrying capacity could be approximated by estimating the abundance of prey in a unit of habitat. With this information, managers would be able to predict the response of breeding female home range sizes to changes in prey abundance, as well as to identify if tiger poaching, prey poaching or habitat quality were limiting tiger population size.
We used data from six resident breeding female tigers fitted with global positioning system (GPS) satellite collars in Huai Kha Khaeng Wildlife Sanctuary in western Thailand during 2005–2011 to explore the relationship between female home range size and prey density. We compared our results to those from studies in four other protected areas in which tiger home range size and prey abundance have been previously estimated. The objectives of this study are to (1) estimate home range size and home range overlap of tigers in western Thailand, (2) correlate female tiger home range size and prey density across the tiger's range, and (3) demonstrate the significance of female home range size and prey abundance as two important metrics for reserve managers.
Study area
The study was conducted in the Huai Kha Khaeng Wildlife Sanctuary (Fig. 1), which, in conjunction with the wildlife sanctuaries Thung Yai East and West, is designated as a UNESCO World Heritage Site. This protected area complex supports the largest tiger population in Thailand (DNP, 2010). In Huai Kha Khaeng Wildlife Sanctuary the estimated tiger density is 3.98 per 100 km2 (Simcharoen et al., Reference Simcharoen, Pattanavibool, Karanth, Nichols and Kumar2007). The sanctuary is 2,780 km2, altitudes are 200–1,600 m, annual temperature range is 8–38° C and mean annual rainfall is 1,375 mm (Khao Nang Rum Wildlife Research Station, unpubl. data). Normally the lowest temperatures occur in January and the highest in April. The dry season (November–April) has a mean rainfall of 298 mm and the wet season (May–October) a mean of 1,088 mm. The area has four main vegetation types, the occurrence of which depends on rainfall patterns and edaphic conditions: mixed deciduous forest (48%), dry evergreen forest (25%), hill evergreen forest (13%) and dry deciduous dipterocarp forest (7%; WEFCOM, 2004).
In addition to the tiger, other carnivores in the sanctuary include the leopard Panthera pardus, clouded leopard Neofelis nebulosa, marbled cat Pardofelis marmorata, leopard cat Prionailurus bengalensis, Asiatic jackal Canis aureus, wild dog Cuon alpinus, Asiatic black bear Ursus thibetanus and sun bear Helarctos malayanus. Four major tiger prey species occur in high numbers in the area: banteng Bos javanicus, gaur Bos gaurus, sambar Rusa unicolor and wild boar Sus scrofa (Srikosamatara, Reference Srikosamatara1993).
Methods
Capture and radio collaring
During June 2005–August 2011 we captured tigers in parts of Huai Kha Khaeng Wildlife Sanctuary that represent the range of habitat types and topographic features in the Sanctuary. Healthy live domestic cows were placed along trails where repeated recent tiger scent marks occurred. When a cow was killed three or four leg-hold snares were set near the cow to capture the tiger when it returned to feed the next night. The snares were designed by P. Pratumratanatarn, an experienced technician with the Department of National Parks, Wildlife and Plant Conservation. Captured tigers were anaesthetized with a mixture of tiletamine HCL and zolazepam HCL (Zoletil; Virbac Laboratories, Carros, France) at a dose of 4 mg kg−1 (Kreeger & Arnemo, Reference Kreeger and Arnemo2007). At the time of capture we recorded sex, age and reproductive status. The age of each tiger was estimated from tooth eruption, tooth wear and staining. Reproductive status of females was classified as nulliparous if nipples were pink; dark nipples indicated a female had produced young. Animal handling and immobilization were undertaken in accordance with the University of Minnesota IACUC protocol 0906A67489.
Breeding tigers were fitted with one of three radio collar models: Advanced Telemetry Systems (ATS) GPS collar model G2000 (Isante, Minnesota, USA), Telonics Argos Terrestrial Transmitter (Mesa, Arizona, USA), or Vectronic Aerospace GmbH GPS Plus (Berlin, Germany). The ATS collars were programmed to acquire locations every 2 hours. Telonics and Vectronic collars were programmed to acquire locations every hour. All radio collars were released with drop-off mechanisms. ATS and Telonics models had to be recovered to download data but the Vectronic model transmitted data via Iridium satellites and the data were available from the Vectronic website. Output of all collars included date, time, latitude, longitude, altitude and fix status. We used location data from both two- and three-dimensional fixes.
Home range size and overlap
We used 95 and 100% minimum convex polygons (MCP) to estimate home range size of tigers (Jennrich & Turner, Reference Jennrich and Turner1969). Because the 100% MCP can be influenced by outliers, the 95% MCP is often preferred, to avoid inflation of the estimate (Harris et al., Reference Harris, Cresswell, Forde, Trewhella, Woodlard and Wray1990; White & Garrott, Reference White and Garrott1990). The MCP is also the most widely used home range estimator (Harris et al., Reference Harris, Cresswell, Forde, Trewhella, Woodlard and Wray1990), which facilitates comparisons with other studies. To determine the number of locations needed to obtain stable home range areas we conducted a bootstrap simulation in which home range size was estimated from randomly selected locations. We began with 50 locations and increased the sample size by 50 until all radio locations were included in the calculations. We considered the home range established when the number of days or the number of fixes reached an asymptote (Harris et al., Reference Harris, Cresswell, Forde, Trewhella, Woodlard and Wray1990). The kernel method, in contrast to the MCP, is best as an estimator of the probability of use within a home range and thus can address questions related to third-order selection (sensu Johnson, Reference Johnson1980). But the kernel method can be problematic for estimating home range size, especially as the number of locations increases with the use of GPS collars (Downs & Horner, Reference Downs and Horner2008). Kernel estimates are particularly inflated when there is a large number of locations along a ‘hard’ boundary, which is often the case for animals, such as tigers, that defend home ranges (Downs & Horner, Reference Downs and Horner2008). Therefore, we believe that the 95% MCP is a suitable estimator for testing our hypothesis that prey density is an important factor determining the size of tiger home ranges. However, we include 100% MCP estimates here also to facilitate comparison with other studies.
To determine if a female home range was used exclusively by a single female to obtain food for herself and her young, we defined neighbouring females as those that had adjacent (with insufficient space or evidence for there to be another female living in between) and potentially overlapping home ranges during the same time period. Overlaps were based on the frequency a female used the same geographical locations as her female neighbours (at least once per month). As we did not have data for all female home ranges surrounding our study females we could only estimate the minimum percentage overlap (100% MCPs), defined here as HR ij = Aij/Ai , where HR ij is the proportion of individual i's home range that is overlapped by individual j's home range, A i is the area of individual i's home range and A ij is the area of overlap between the two individuals’ home ranges (Kernohan et al., Reference Kernohan, Gitzen, Millspaugh, Millspaugh and Marzluff2001). The home range sizes and degree of overlap of radio collared tigers were calculated using the ArcView v. 3.3 (ESRI, Redlands, USA) extension package Home Range (Rodgers & Carr, Reference Rodgers and Carr1998).
Prey abundance
We estimated prey abundance during the dry season (December–April) for 2009–2011. In addition, prey abundance in the entire Sanctuary has been monitored since 2006 as part of sanctuary management; these data suggested that prey density was approximately constant during 2006–2011 (Wildlife Conservation Society, unpubl. data).
We used faecal accumulation rate techniques (Bailey & Putman, Reference Bailey and Putman1981) to estimate prey abundance in the Sanctuary. Four c. 100 km2 areas that encompassed the six female tiger home ranges were selected to represent the range of habitat types and topographic features in the Sanctuary. Two sites were in the central valley of the Sanctuary along the lower and upper portions of the Huai Kha Khaeng River. The lower site was characterized by mixed deciduous forest. The upper site was a combination of mixed deciduous and dry evergreen forest and the habitat was more rugged than at the lower site. A third site, near the Khao Nang Rum Wildlife Research Station, was a drier area away from the Huai Kha Khaeng River. It was composed of a combination of dry dipterocarp and mixed deciduous habitat. The fourth site was in the northern part of the Sanctuary, and also away from the Huai Kha Khaeng River. It was characterized by mixed deciduous and dry evergreen habitat in rugged terrain.
A total of 90 square line transects, 200 × 200 m, were randomly located in each of the four sites (Fig. 1) to survey gaur and banteng, as dung of these species is relatively rare. In addition, 40 circular plots of 20 m2 were placed 20 m apart on each of these transects, to sample the more common ungulates (sambar, barking deer Muntiacus muntjac and wild boar). All ungulate pellet groups were counted and removed from all transects and circular plots. Following removal, new faecal accumulation was measured for a period of 30 days (over this time period no entire pellet groups were lost to decomposition or other natural processes during the dry season). A pellet group consisted of >15 pellets. Pellet groups that overlapped were distinguished by size and colour, which became lighter as pellets aged.
We searched for gaur and banteng dung piles up to 2 m either side of the transects. Observed piles were counted and marked, and the perpendicular distance from the transect was measured. We could not distinguish gaur and banteng dung and therefore combined the data, as ‘large bovid’. To determine absolute abundance of prey species we used published defecation rates of large bovids (9.5 dung piles per day) and barking deer (7.5 pellet groups per day; Srikosamatara, Reference Srikosamatara1993; Sukmasuang, Reference Sukmasuang2001). For sambar we recorded the defecation rate of a group of six animals of mixed sex and age (one male, two females and three young) at Khao Pratapchang Wildlife Breeding Centre. To control for possible seasonal differences in defecation rate, observations at the Centre were made in the same season as the field study. During a 5-day period deer grazed undisturbed; in a second 5-day period deer were stimulated so that they were more active than normal. Deer were fed a diet of natural browse that was cut near the breeding centre, to provide a diet similar in species and cellulose composition to their natural diet. The mean number of pellet groups per individual per day was 9.5 ± SD 1.41 (n = 6). We calculated deer and barking deer density as D = (P/A)/(Df Ds), where P = the number of pellet groups found on the second visit, A = the total sample area (km2), Df = the defecation rate and Ds = the number of days between the two visits. Distance v. 5.0 (Thomas et al., Reference Thomas, Laake, Strindberg, Marques, Buckland, Borchers and Marques2006) was used to estimate the density of dung, which was then used to estimate large bovid abundance. Data were right truncated at 1.2 m and grouped into 0.2 m intervals to improve model fit. Uniform and half-normal models for the detection function were fitted against the data, using cosine adjustments (Buckland et al., Reference Buckland, Anderson, Burnham, Laake, Borchers and Thomas2001). The best model selection and adjustment terms were based on Akaike's Information Criterion (AIC; Akaike, Reference Akaike1974). We calculated large bovid density BD = Dd/(DfDs), where Dd = bovid dung density, Df = defecation rate of bovids, and Ds = number of days between the two visits.
Home range size and prey abundance
We estimated prey biomass in six resident female tiger home ranges in the Sanctuary to determine the correlation between female home range size and prey abundance at the local level. Across the tiger's range we reviewed information on tiger home range size and prey abundance from five studies in Nepal, India, Bangladesh and Russia that had used radio telemetry to estimate home range size. In all studies home range size was estimated by either the 95 or 100% MCP. Using data from earlier studies we estimated prey biomass for each site as the product of prey density and mean body weight (Dhungel & O'Gara, 1991; Karanth & Sunquist, Reference Karanth and Sunquist1992). We then used simple linear regression to predict home range size in relation to prey abundance at the local level and Spearman's rank correlation to determine if the same relationship holds at a range-wide level (our results from Huai Kha Khaeng Wildlife Sanctuary combined with published data). Probabilities < 0.05 were considered significant.
Results
Home range size and overlap
A total of 11 tigers, seven resident breeding females and four resident breeding males, were fitted with GPS collars. One female and one male were collared twice, at separate times, and home range estimates were made for each period for each individual as size and shape of home ranges changed over time. The mean number of days used for home range analysis was 152 ± SD 72 days (range 77–358 days). The mean time to reach an asymptote in home range size was 77 ± SD 35 days and the mean number of fixes to reach an asymptote was 418 ± SD 160 locations. The mean number of locations used to estimate home range was 1,502 ± SD 1,215. Mean home range (95% MCP) of the resident females (n = 8) was 70.2 ± SD 33.2 km2 and that of the resident males was 267.6 ± SD 92.4 km2 (n = 5; Table 1). Mean home range estimated with the 100% MCP was 84.2 ± SD 40.8 km2 (n = 8) and 294.1 ± SD 100.3 km2 (n = 5) for the resident males and females, respectively (Table 1). The ratio of male to female tiger home range size was 3.8 : 1.
Home range overlap was examined from three pairs of adjacent females and one pair of adjacent males. Mean overlap was 4.5 ± SD 3% for the females and 18% for the males. We used five female home ranges within three male home ranges to examine overlap between male home ranges and sympatric female home ranges. Male home ranges overlapped nearly the entire sympatric home range of the females (90 ± SD 15%, n = 5).
Prey abundance
Mean densities of sambar and barking deer were 7.5 ± SD 6 and 3.5 ± SD 2.6 km−2 over the four sites. Our sample of wild boar dung (n = 77) was insufficient for precise density estimation and therefore we used the density of wild boar in the Sanctuary (2.4 ± SD 0.05 km−2) estimated by Sukmasuang (Reference Sukmasuang2001). To estimate large bovid density by distance sampling we first tested models of dung observability based on distance from the transects. During the surveys 105 bovid dung piles were detected. The best model for observability based on the lowest AIC score was a uniform cosine distribution. The bovid dung density was estimated to be 860 ± SD 126 km−2 and, from this, bovid density was estimated to be 3.0 ± SD 0.4 km−2.
Home range size and prey abundance
The simple linear regression indicated there was a significant negative correlation between size of female home range and prey biomass (r2 = 0.70, P < 0.05, n = 6). This relationship appeared to be log-linear as a log transformation of the home range size gave a better fit (r2 = 0.85, P < 0.05, n = 6; Fig. 2). At a regional scale (Table 2) we found the same relationship (Spearman rs = 0.88, P < 0.005, n = 6).
Discussion
Our results support the hypothesis that variation in the home range size of female tigers is partially explained by variation in prey density. The relationship appears to be log-linear, probably reflecting the biological upper and lower limits of home range sizes within the study area. This relationship indicates that home range size of female tigers is relatively constant when prey biomass is > 5,000 kg km−2, indicating minimum home ranges of 10–20 km2 (Fig. 2). In carnivores a direct relationship between higher prey abundance and smaller home range sizes has been reported in many species, including bobcat Felis rufus (Litvaitis et al., Reference Litvaitis, Sherburne and Bissonette1986) and Canadian lynx Lynx canadensis (Ward & Krebs, Reference Ward and Krebs1985). Schaller (Reference Schaller1972), Orsdol et al. (Reference Orsdol, Hanby and Bygott1985) and Hayward et al. (Reference Hayward, O'Brien and Kerley2007) also demonstrated this relationship for large African felids.
To test our hypothesis that a female's home range should be large enough to support the prey requirements for a female and her cubs from birth until they disperse at c. 18 months (Smith et al., Reference Smith, McDougal, Sunquist, Tilson and Seal1987), we needed a home range estimator that defined the area in which a female acquires food. We also needed to test the degree to which females have exclusive home ranges. Choice of home range estimators depends on the question being addressed (Kie et al., Reference Kie, Matthiopoulos, Fieberg, Powell, Cagnacci and Mitchell2010) and therefore we calculated home range as both the 95 and 100% MCPs (Smith et al., Reference Smith, McDougal, Sunquist, Tilson and Seal1987; Chundawat et al., Reference Chundawat, Gogate, Johnsingh, Seidensticker, Christie and Jackson1999; Karanth & Sunquist, Reference Karanth and Sunquist2000; Goodrich et al., Reference Goodrich, Miquelle, Smirnov, Kerley, Quigley and Hornocker2010; Barlow et al., Reference Barlow, Smith, Ahmad, Hossain, Rahman and Howlader2011). For female tigers home ranges are not only related to prey abundance but are exclusively used by an individual female. Tigers in Russia and Nepal exhibit territorial behaviour as a strategy to ensure access to adequate prey to raise their young to dispersal age (Smith et al., Reference Smith, McDougal, Sunquist, Tilson and Seal1987, Goodrich et al., Reference Goodrich, Miquelle, Smirnov, Kerley, Quigley and Hornocker2010). This behaviour is also common in other solitary female carnivores (Sandell, Reference Sandell and Gittleman1989). We did not document aggressive behaviour among females as an indicator of territorial behaviour but we did find that territorial overlap among females is small. For example, three pairs of adjacent females had an overlap of only 4%, which is similar to estimates from Nepal, where the overlap of females was 3.5–7% (Smith et al., Reference Smith, McDougal, Sunquist, Tilson and Seal1987), and Russia, where female tigers had a mean overlap of 9% (Goodrich et al., Reference Goodrich, Miquelle, Smirnov, Kerley, Quigley and Hornocker2010).
Factors other than prey, such as poaching of tigers, may reduce competition among females and result in home ranges larger than needed to supply food resources. These conditions may confound the negative relationship between female tiger home range size and prey abundance. When density of tigers was low in Nepal during periods when habitat was recovering from human disturbance, tiger home ranges were larger than needed and this allowed females to make room for their dispersing daughters (Smith et al., Reference Smith, McDougal, Sunquist, Tilson and Seal1987).
In contrast it is widely reported that home range size for male solitary carnivores is not influenced by food requirements, as males attempt to hold territories that encompass as many females as possible (Sandell, Reference Sandell and Gittleman1989). The home range size, and thus the number of females that a male may mate with, is constrained by the energy and time requirements expended visiting females to track their reproductive status, and patrolling boundaries to discourage intruding males that may kill cubs (Sandell, Reference Sandell and Gittleman1989). Mean male home range size in Huai Kha Khaeng Wildlife Sanctuary was 3.8 times larger than the mean female home range. This ratio is similar to the ratio of home range size between male and female tigers in Russia and Nepal (Smith et al., Reference Smith, McDougal, Sunquist, Tilson and Seal1987; Goodrich et al., Reference Goodrich, Miquelle, Smirnov, Kerley, Quigley and Hornocker2010).
Although our sample size for the correlation between home range and prey is small, it is currently the best information available. Tracking the relationship between female home range size and prey abundance is a useful tool for managers because it provides them with a metric for setting management goals. Hayward et al. (Reference Hayward, O'Brien and Kerley2007) suggested that the next step in refining the relationship between female home range size and prey abundance is to analyse the ecological factors that predict prey abundance. Detailed information on female home range size, prey abundance and the ecological factors that determine prey abundance will allow managers to develop appropriate management actions. For example, when female home range sizes are larger than needed to supply adequate prey for breeding females to raise their young, managers should investigate the possibility that tigers are being poached. Alternatively, when prey abundance is below the potential carrying capacity predicted by environmental variables, managers should examine whether prey are being poached. In contrast, if the mean size of home ranges is at the minimum size given the prey density, then the only management option to increase the viability of the tiger population is to enlarge the area available. If in the next few years female home range size increases in Huai Kha Khaeng Wildlife Sanctuary, managers should investigate if it is because of declining prey abundance or an increase in tiger poaching. Faced with continued threats to tigers, conservationists have focused on strategies that emphasize conserving source populations (Walston et al., Reference Walston, Robinson, Bennett, Breitenmoser, Fonseca and Goodrich2010) or entire landscapes (Sanderson, Reference Sanderson, Tilson and Nyhus2010). These efforts provide important guidelines, and strategic conservation planning is essential. Poaching, reduced carrying capacity and the loss of the land area available to tigers have been identified as the key threats to the species (Kenny et al., Reference Kenny, Smith, Starfield and McDougal1995; Chapron et al., Reference Chapron, Miquelle, Lambert, Goodrich, Legendre and Clobert2008; Seidensticker, Reference Seidensticker2010). Managers can use the relationship between female home range size and prey abundance, along with long-term camera-trap monitoring of tiger populations, to examine the threats specific to tigers in each protected area.
Acknowledgements
We acknowledge the Department of National Parks, Wildlife and Plant Conservation, and especially Chatchawan Pitdumkam and Theerapat Prayurasiddhi, for supporting this work. We also appreciate the help provided by the former chiefs of Huai Kha Khaeng Wildlife Sanctuary, Sunthorn Chaiwatana and Apicha Yusombune. Funding was provided by PTTEP (Thailand), the U.S. Fish and Wildlife Service, The Royal Golden Jubilee PhD Programme (Thailand) and the Biodiversity Research and Training Programme (Thailand). We are also grateful to research assistants Onsa Norasarn, Thavorn Thadvichit, Kerkpol Wongchu, Boonyang Srichan, Yingbune Chongsomchai, Kerati Phetong, Nattakarn Pengmark, Namkang Saelee and Parinyakorn Worawan for their help with the capture of tigers and field data collection. We are particularly grateful to Francesca J. Cuthbert, Erin Roche, Anak Pattanavibool, Vijak Chimchom, John Fieberg and Adam Barlow for their comments the article and for suggestions regarding data analysis.
Biographical sketches
Achara Simcharoen specializes in tiger ecology and tiger movement behaviour. Tommaso Savini's research includes the behavioural ecology of primates and pheasants. George A. Gale conducts research on the ecology and population dynamics of birds and mammals. Saksit Simcharoen studies the ecology and habitat use of tigers and leopards. Somphot Duangchantrasiri is in charge of tiger monitoring in the Western Forest Complex of Thailand. Somporn Pakpien undertakes research on tigers and their prey. James L. D. Smith investigates landscape-scale tiger ecology.