Introduction
Tiger Panthera tigris populations continue to decline throughout their range (Dinerstein et al., Reference Dinerstein, Loucks, Heydlauff, Wikramanayake, Bryja and Forrest2006). This decline is because of the illegal trade in their body parts for traditional medicine, authorized removal of problem tigers following tiger-human conflicts (Tilson et al., Reference Tilson, Soemarna, Ramono, Lusli, Traylor-Holzer and Seal1994; Seidensticker et al., Reference Seidensticker, Christie, Jackson, Seidensticker, Christie and Jackson1999), depletion of prey species (Seidensticker, Reference Seidensticker, Miller and Everett1986; Karanth & Stith, Reference Karanth, Stith, Seidensticker, Christie and Jackson1999), and accelerated destruction of their natural habitat (Dinerstein et al., Reference Dinerstein, Loucks, Heydlauff, Wikramanayake, Bryja and Forrest2006).
Indonesia historically contained three subspecies of tigers but both the Javan subspecies P. t. sondaica and the Bali subspecies P. t. balica went extinct in the early 1980s and 1940s, respectively (Seidensticker, Reference Seidensticker, Miller and Everett1986; Seidensticker et al., Reference Seidensticker, Christie, Jackson, Seidensticker, Christie and Jackson1999). Currently, only the Sumatran subspecies P. t. sumatrae persists, in isolated populations throughout Sumatra (Dinerstein et al., Reference Dinerstein, Loucks, Heydlauff, Wikramanayake, Bryja and Forrest2006). In 1978 the Sumatran tiger population was estimated to be 1,000 (Borner, Reference Borner1978). In 1985 26 protected areas were found to contain a total of 800 tigers (Santiapillai & Ramono, Reference Santiapillai, Ramono and Tilson1987). In 1992, a Population and Habitat Viability Analysis estimated there were only 400–500 Sumatran tigers in the wild (Tilson et al., Reference Tilson, Soemarna, Ramono, Lusli, Traylor-Holzer and Seal1994). The aims of the present study were to make the first estimate of the density of tigers in the recently gazetted Batang Gadis National Park and to identify the environmental variables that best predict tiger occurrence in the Park.
Study area
Batang Gadis National Park in North Sumatra Province, Indonesia, is one of 12 forest blocks in Sumatra, including national parks, identified as a level 3 Tiger Conservation Landscape. Sanderson et al. (Reference Sanderson, Forrest, Loucks, Ginsberg, Dinerstein and Seidensticker2006) defined level 3 Tiger Conservation Landscapes as those that have habitat to support some tigers but that have moderate to high levels of threat and minimal conservation investment. Under a new legal framework that allows local governments to declare national parks, the governor of North Sumatra and the head of Mandailing Natal district declared 26% of the total forested area in the district as Batang Gadis National Park in December 2004. Covering a total of 1,080 km2 of tropical rainforest, the Park has an altitude range of 300–2,145 m (Sorik Merapi mountain; Fig. 1). Rainfall is seasonal, with an annual range of 1,900–2,800 mm over 2002–2006 (BMG, 2008).
Fig. 1 Batang Gadis National Park, Mandailing Natal District, North Sumatra, with the sampling block of 256 km2 (16 × 16 km divided into 4 × 4 km cells) where camera-trapping took place (see text for further details), and poor, medium and high quality tiger habitat as determined by a model of tiger presence (see Table 1 and text for further details). The rectangle on the inset indicates the location of the main map in Sumatra.
Methods
Large mammal surveys were conducted in eight periods between December 2005 and July 2006. Twenty-six passive infrared camera traps (Deercam, Park Falls, USA) were used, each with data packs that record time and date on the photographs. Cameras were placed in a 16 × 16 km2 sampling block in the centre of the Park (Fig. 1). The sampling block was divided into 16 equal-sized cells, with two cameras in each cell. Camera sites were chosen based on areas that had tiger or prey signs. Cameras were active 24 hours per day and operated for 20–30 days in each period. The number of trap days for each film was defined as from the activation of the camera until the film was retrieved, if the film had exposures remaining, or until the time and date stamp on the final exposure. Each animal photo was identified to species. We used the number of photographs of each species to construct a respective relative abundance index (O'Brien et al., Reference O'Brien, Kinnaird and Wibisono2003). We deployed the cameras at 32 different locations with a mean altitude of 983 ± SD 253 m. The average distance between cameras was 2,062 m (range 286–4,391 m). This prevented holes in the sampled area where the probability of capturing tigers was zero (Karanth et al., Reference Karanth, Kumar, Nichols, K.U. and J.D.2002).
Linkie et al. (Reference Linkie, Chapron, Martyr, Holden and Leader-Williams2006) suggested the following environmental variables as potential predictors for tiger occurrence: mean distance to public and logging roads, settlements and rivers, and altitude, slope and protected status. In addition, O'Brien et al. (Reference O'Brien, Kinnaird and Wibisono2003) and Linkie et al. (Reference Linkie, Haidir, Nugroho and Dinata2008) found a negative correlation between tiger occurrence and distance to forest edges. To test the relationship between tiger abundance and environmental variables a geographical information system (GIS) data set containing seven variables was constructed with ArcGIS v. 9.1 (ESRI, Redlands, USA): slope (degree); distance to rivers (50 m intervals); river length (km); distance from forest edge to its interior (50 m intervals); distance to district roads (km); distance to villages (km); and altitude (50 m intervals; Table 1). Altitude and slope data were obtained from Digital Elevation Map SRTM 07 (USGS, 2004) and data for the remaining variables from the National Coordinating Agency for Surveys and Mapping (BAKOSURTANAL) in 2005.
Table 1 The best multiple logistic regression model, obtained using stepwise selection procedures, for predicting tiger presence. A total of seven environmental variables were evaluated for inclusion (see text for further details).
We used the software CAPTURE (Patuxent Wildlife Research Center, Laurel, USA; Rexstad & Burnham, Reference Rexstad and Burnham1991) to estimate tiger capture probabilities and then abundance as recommended by Karanth & Nichols (Reference Karanth, Kumar, Nichols, K.U. and J.D.2002). In capture-recapture models abundance estimation requires that the population be closed i.e. no births, deaths, immigration or emigration during the sampling period (Otis et al., Reference Otis, Burnham, White and Anderson1978). The closure assumption was examined using the closure test in CAPTURE (Otis et al., Reference Otis, Burnham, White and Anderson1978). For each individual tiger, identified from its stripe pattern, we constructed unique capture histories (O'Brien et al., Reference O'Brien, Kinnaird and Wibisono2003). All analyses were confined to photographs of the right flank because this side had the largest number of identified and recaptured individuals. Each trapping day was treated as one sampling occasion, giving an average of 20 sampling occasions for the analysis. We used model Mh, which allows for a variable capture probability between individuals and is more reliable and robust to the violation of assumptions (Otis et al., Reference Otis, Burnham, White and Anderson1978). To calculate tiger density the estimated tiger abundance from model Mh was divided by the effective sampling area, defined as the outermost camera locations, plus an added buffer zone with equal width to the mean maximum distance moved (MMDM) of recaptured tigers (Nichols & Karanth, Reference Nichols, Karanth, K.U. and J.D.2002).
We used logistic regression to evaluate the probability of tiger presence based on the seven environmental variables (log transformed). Tiger occurrence was the binary response variable and the environmental factors the explanatory variables. We based our decisions on stepwise selection procedures. The fit of the data with the model was assessed using deviance and Pearson goodness-of-fit statistics. Potential collinearity of explanatory variables was assessed using a nonparametric Spearman correlation analysis. We used rs ± 0.33 as a cut off for strong collinearity between variables for n = 60 and P = 0.01 (Gibbons, Reference Gibbons1978). If two variables were highly correlated the one with the lower Wald χ2 in the full model was removed. Wald χ2 is the squared ratio between the point estimate and associated standard error. Thus we removed the variables with the lower Wald χ2 because they were weaker in predicting tiger occurrence.
A predictive model of tiger distribution was generated using distance from forest boundaries and elevation, the two environmental variables that best predicted tiger presence. A GIS analysis was used to create a predictive model of tiger habitat quality. Using the reclassify tool the altitude and distance from forest boundary layers were reclassified into 10 equal intervals and ranked between 1 and 10 with respect to the effect of each variable in predicting tiger occurrence. The two new layers were then combined using the map algebra tool of ArcGIS Spatial Analyst. The resulting layer was ranked into three categories of tiger habitat: poor, medium and high quality. The extent of this predictive map was based on forest boundaries.
Results
A total of 102 films were exposed at the 32 locations. We obtained a total of 364 animal photos, including tigers, representing at least 24 species in a total of 1,728 trap-days. A total of 10 tiger photographs of the right flank, including duplicates and unidentified individuals, were obtained. The quality of one photograph was too poor for identification. The number of unique individuals identified was six. Using CAPTURE model Mh we estimated an average capture probability of 0.08. The estimated tiger population size in the sampled area was 6 ± SE 2.4, with a 95% confidence interval (CI) of 6–21. The coefficient of variation (CV) was 40.0, indicating a lack of precision. The small sample size of only six unique tigers and three recaptures precluded more accurate estimates. Small sample sizes are common in tiger studies (Karanth & Nichols, Reference Karanth and Nichols2000; Lynam et al., Reference Lynam, Rabinowitz, Myint, Maung, Latt and Po2008).
A closure test indicated that our population did not differ significantly from a closed population (z = 0.126, P = 0.55). The MMDM between recaptured tiger individuals was 4.5 km, leading to an estimated total sampling coverage of 431 km2. We calculated the habitat actually sampled by subtracting the amount of non-forested areas from this total coverage, which gave an area of 327 km2. Using the estimated tiger population size in the sampled area with the effective sampling coverage of suitable habitat, we estimated the density of tigers to be 1.8 per 100 km2 (95% CI = 1.8–6.4). Assuming this estimate reflects tiger densities in suitable habitat in and adjacent to the Park of c. 1,607 km2, we estimated that 29–103 tigers live in this area.
We found significant collinearity between distance from rivers and river length (rs = -0.62, P < 0.001), distance from forest edges and villages (rs = 0.749, P < 0.001), distance from roads and villages (rs = 0.797, P < 0.001), distance from roads and distance from forest (rs = 0.650, P < 0.001), altitude and distance from rivers (rs = 0.357, P < 0.001), and altitude and river length (rs = 0.555, P < 0.01). Based on the Wald χ2 values of the full model we retained slope, altitude, river length and distance from forest edges in the logistic model. The stepwise selection procedure indicated that an increase in one unit (50 m) of altitude decreased the probability of tigers occurring (P < 0.019, odds ratio = 0.004, 95% CI = 0.001–0.403). Inversely, an increase in one unit (50 m) of distance from forest edges increased the probability of tigers occurring (P < 0.043, odds ratio = 5.646, 95% CI = 1.059–30.093). The Pearson (P > 0.194) and deviance (P > 0.6677) goodness-of-fit tests indicated a good fit because the observed values did not significantly differ from the values expected under the model in question.
The predictive model of tiger distribution was developed based on altitude and distance from forest boundaries, the two environmental variables included as the best predictors of tiger presence in the stepwise selection procedure. Based on the model we predicted 24.8% of the total landscape to be poor, 57.6% medium and 17.6% high quality habitat. Most of the high quality habitat was at the periphery of the Park (Fig. 1).
Discussion
This study is both the first to estimate tiger population size and density in Batang Gadis National Park. The results provide an estimate of the tiger population in the Park for the global tiger conservation assessment (Sanderson et al., Reference Sanderson, Forrest, Loucks, Ginsberg, Dinerstein and Seidensticker2006). The estimated range of population sizes from this study, together with the current development of infrastructure and conservation intervention in the Park, could potentially promote the Park's current status as Tiger Conservation Landscape level 3 and long-term priority (Sanderson et al., Reference Sanderson, Forrest, Loucks, Ginsberg, Dinerstein and Seidensticker2006) to Tiger Conservation Landscape level 2 and regional priority.
The Mh model estimated a mean of six tigers in the sampled area. The high CV of the model may be a result of low sample size and unequal sampling efforts between replicates. The mean density of 1.8 tigers per 100 km2 is comparable with other studies in high elevation tiger habitats in Sumatra. Linkie et al. (Reference Linkie, Martyr, Holden, Yanuar, Hartana, Sugardjito and Leader-Williams2003) found a relationship between altitude and relative abundance in Kerinci Seblat National Park in Central Sumatra, and estimated tiger densities of 2.0 and 1.5 per 100 km2 in hill and submontane forests, respectively (Linkie et al., Reference Linkie, Chapron, Martyr, Holden and Leader-Williams2006). Similarly, a tiger density of 1.7 per 100 km2 was estimated in the Bukit Barisan Selatan National Park in Southern Sumatra (O'Brien et al. Reference O'Brien, Kinnaird and Wibisono2003). Other studies in mainland Asia have also shown a wide variance in estimated tiger densities, with a range of 1.1–6.2 tigers per 100 km2 in Peninsular Malaysia (Kawanishi & Sunquist, Reference Kawanishi and Sunquist2004), 2.4–8.7 per 100 km2 in Thailand (Simcharoen et al., Reference Simcharoen, Pattanavibool, Karanth, Nichols and Kumar2007) and 1.98–7.31 per 100 km2 in Badhra, India (Karanth & Nichols, Reference Karanth and Nichols2000).
Tigers are considered to be habitat generalists (Schaller, Reference Schaller1967). However, this study found that the abundance of Sumatran tigers decreased with increasing altitude, as is the case in Bukit Barisan Selatan National Park (Wibisono, Reference Wibisono2006) and Kerinci Seblat National Park (Linkie et al., Reference Linkie, Chapron, Martyr, Holden and Leader-Williams2006). Higher elevations usually have more rugged terrain and thus a poorer prey base. Our predictive model did, however, indicate that good quality tiger habitat tended to be high in forest-town peripheries at some locations, in contrast to a recent study in selectively logged primary forest around Kerinci Seblat National Park (Linkie et al., Reference Linkie, Haidir, Nugroho and Dinata2008). Thus, our model could serve as a preliminary indicator of where conservation intervention should be concentrated with regard to tiger-human conflict mitigation and habitat protection.
Batang Gadis is a new National Park connecting two larger forest blocks of the Angkola and Barumun-Rokan ecosystems (Fig. 1). We detected the same tiger individual at two different camera locations at the northern and southern sections of the Park. It was unlikely that the tiger travelled through the villages dissecting the two sections. Rather, it probably moved through the connecting narrow isthmus of forest (Fig. 1). The Park therefore potentially serves as an important corridor for the Angkola and Barumun-Rokan ecosystems. The presence of tigers in such ecosystems has been confirmed by North Sumatra's Biodiversity Conservation Agency. Combining the Batang Gadis National Park with the Angkola and Barumun-Rokan ecosystems may increase the area's Tiger Conservation Landscape status from level 3 to level 1 (Sanderson et al., Reference Sanderson, Forrest, Loucks, Ginsberg, Dinerstein and Seidensticker2006). We estimate that Angkola, Batang Gadis National Park and Barumun-Rokan have a combined total area of c. 6,500 km2. Applying the estimated Batang Gadis National Park tiger density to this entire landscape gives an estimated tiger population of 117, representing a substantial contribution to the global Sumatran tiger population. Therefore, we believe that Batang Gadis National Park, together with the Angkola and Barumun–Rokan ecosystems could serve as a stronghold for tiger populations in northern Sumatra. The next step would be to develop an appropriate tiger management plan for this landscape, identify the conservation activities required, and outline the roles and responsibilities of the different agencies involved. The entire Batang Gadis, Angkola and Barumun-Rokan ecosystem should then be reevaluated and possibly upgraded to a higher level Tiger Conservation Landscape.
Acknowledgements
This study was part of a collaborative effort between Conservation International Indonesia and the Indonesian Ministry of Forestry's Department for Forest Protection and Nature Conservation, and was technically assisted by the Wildlife Conservation Society Indonesia Program. The work was funded by the US National Fish and Wildlife Foundation, Save the Tiger Fund, the Exxon Mobil Corporation, and Global Conservation Fund. We would like to thank Dr Jatna Supriatna, Dr Susie Ellis, Dr Noviar Andayani, Iwan Widjayanto, Dr Nick Brickle, Dr Mahendra Shrestha, Dr Didy Wurjanto, Maggie Meyer, Nenita de Guzman, Dr Jim Sanderson, Dr Erwin Widodo, Effendy Sumardja, Purbasari Surjadi, Kusuma Paramitha, Dian Melur, Madonna, Agus W. Boyce, Tengku Afriyenni, all CI staff, and Frida Mindasari Saanin from the Wildlife Conservation Society Indonesia Program. Dr Joel T. Heinen provided helpful advice and encouragement. We also thank Batang Gadis National Park staff and the many individuals living around the Park who have supported the project and Sumatran tiger conservation.
Biographical sketches
Hariyo T. Wibisono has worked as a field scientist and conservation manager since 1993. His research interest is in the population ecology of large mammals, especially Sumatran tigers. Joe J. Figel is currently working on jaguars in Latin America. Sugesti M. Arif has worked as a field technician since 1998. In this study he was responsible for the day-to-day monitoring of the camera traps. Anton Ario has worked as a conservation manager since 1998 and specializes in conservation education. Abu H. Lubis has also worked as a conservation manager since 1998 and specializes in primate ecology.
