This paper contains supplementary material that can be found online at Http://journals.cambridge.org
Recently, the revision of the taxonomy of clouded leopards using molecular and morphological data has shown that Bornean and Sumatran clouded leopards form a distinct species, Neofelis diardi, and that each island is occupied by separate subspecies, N. diardi diardi in Sumatra and N. diardi borneensis in Borneo (Wilting et al., Reference Wilting, Buckley-Beason, Feldhaar, Gadau, O'Brien and Linsenmair2007, 2011). This revision was adopted by the IUCN Red List, with N. diardi categorized as Vulnerable and both subspecies as Endangered (Hearn et al., Reference Hearn, Sanderson, Ross, Wilting and Sunarto2011).
Clouded leopards are the least known of the pantherine felids. Earlier density estimates were based on extremely scarce data (Davies & Payne, Reference Davies and Payne1982; Wilting et al., Reference Wilting, Fischer, Abu Bakar and Linsenmair2006), raising concerns regarding reliability. Recently several intensive camera-trapping studies have been conducted in Sabah, Malyasian Borneo. Here we present rigorous density estimates for the Sunda clouded leopard from two commercially used forest reserves in this region.
We carried out camera-trap surveys in the 500 km2 Tangkulap-Pinangah Forest Reserve and 572 km2 Segaliud Lokan Forest Reserve, both in the lowlands of central Sabah. During the last 50 years both areas were repeatedly logged by conventional selective logging techniques. Consequently they are covered today with degraded secondary lowland dipterocarp forests. Logging ceased in Tangkulap-Pinangah Forest Reserve in 2001, and in June 2011 the Reserve received certification from the Forest Stewardship Council. In Segaliud Lokan Forest Reserve, which is privately managed by KTS-Plantations Sdn. Bhd, logging practices changed in 1998 to reduced impact logging, and the Reserve is certified by the Malaysian Timber Certification Scheme.
We set up 64 and 55 camera-trap stations, covering areas of 122 and 114 km2 in Tangkulap-Pinangah and Segaliud Lokan Forest Reserves, respectively (Fig. 1). Stations were spaced c. 1.7 km apart, each comprising two camera traps (Cuddeback Expert & Capture, Non Typical Inc., Green Bay, USA) facing each other to capture both flanks of any passing animal. This permitted identification of individuals based on their unique spot patterns.
Because of our limited number of camera traps we divided Tangkulap-Pinangah Forest Reserve into three blocks and Segaliud Lokan Forest Reserve into two blocks and sampled blocks separately. Each block was sampled for 42 (Tangkulap-Pinangah Forest Reserve) or 48 (Segaliud Lokan Forest Reserve) consecutive days. Total sampling duration was approximately 4 months, a period within the range of durations used in similar studies to approximate a closed population (e.g. Kawanishi & Sunquist, Reference Kawanishi and Sunquist2004; Royle et al., Reference Royle, Magoun, Gardner, Valkenburg and Lowell2011). Tangkulap-Pinangah Forest Reserve was sampled during April–September 2009 and Segaliud Lokan Forest Reserve during January–April 2010.
For each study site we compiled the number of photographs of each individual at each camera-trap station (Appendix 1) and analysed the data using spatial capture–recapture models (Efford, Reference Efford2004; Royle & Young, Reference Royle and Young2008). These models combine a model of individual movements with a model describing detection by camera traps. We adopted a model with a Poisson trap encounter process as described by Royle & Gardner (Reference Royle, Gardner, O'Connell, Nichols and Karanth2011).
Because data at each site were sparse we coupled models from both sites and assumed that the parameter describing individual movements, σ, was constant across sites (i.e. σ is estimated from data generated at both sites). This reflects the assumption that individuals in the two reserves have similar-sized home ranges. We believe this is justified because the habitats in both study sites are similar (AW & AM, unpubl. data), suggesting that clouded leopards should have similar area requirements. Detection and density were estimated separately for each site.
In spatial capture–recapture models density is estimated as the number of individuals, N, occurring within the state space, an area encompassing the trapping grid large enough to contain all individuals that could have been exposed to trapping (Royle & Gardner, Reference Royle, Gardner, O'Connell, Nichols and Karanth2011). We defined the state space for each site by adding a 20-km buffer to the outermost coordinates of the two trapping grids.
We implemented the model in a Bayesian framework using WinBUGS (Gilks et al., Reference Gilks, Thomas and Spiegelhalter1994), accessed through R v. 2.13.1 (R Development Core Team, 2011) using the package R2WinBUGS (Sturtz et al., Reference Sturtz, Ligges and Gelman2005; see Appendix 2 for the WinBUGs code). We present results as means±SE and the 2.5 and 97.5 percentiles of the posterior distributions. The latter characterize the 95% Bayesian credibility interval, analogous to a 95% confidence interval.
We obtained 29 (Tangkulap-Pinangah Forest Reserve) and 15 (Segaliud Lokan Forest Reserve) photographs of Sunda clouded leopards and in both areas we identified five individuals. The number of captures per individual was 1–17 and 1–9 in Tangkulap-Pinangah and Segaliud Lokan Forest Reserves, respectively. In both areas three individuals were males and two (Tangkulap-Pinangah Forest Reserve) and one (Segaliud Lokan Forest Reserve) were females; the sex of the remaining individual could not be determined. Baseline trap encounter rates λ 0 were similar in both reserves and Sunda clouded leopards showed large movements (σ = 6.5 km). With 0.84±SE 0.42 and 1.04±SE 0.58 individuals per 100 km2, density was similar in both sites (Table 1).
* σ=movement parameter, related to home range radius; λ 0=baseline trap encounter rate, the detection parameter of the spatial capture–recapture model; N = number of individuals in state space; D=density
As these are among the first density estimates for the Sunda clouded leopard using spatial–recapture modelling comparison with most existing estimates is difficult. The provisional estimates of 25 individuals (Davies & Payne, Reference Davies and Payne1982) and 8 individuals per 100 km2 (Wilting et al., Reference Wilting, Fischer, Abu Bakar and Linsenmair2006) were based on scarce data. Densities obtained with non-spatial capture–recapture models in two study areas in Sabah (both previously selectively logged but with different forestry management strategies and histories) were 4.8–7.3 individuals per 100 km2 (A.J. Hearn, J. Ross, D. Pamin, H. Bernard, L. Hunter & D.W. Macdonald, unpubl. data). Although these estimates could indicate higher numbers of Sunda clouded leopards in these areas, the differences in densities could also be a methodological artifact, as non-spatial capture–recapture modelling often leads to higher density estimates (e.g. Sharma et al., Reference Sharma, Jhala, Qureshi, Vattakaven, Gopal and Nayak2010). In contrast, Brodie & Giordano (Reference Brodie and Giordano2012) estimated the density of the Sunda clouded leopard in a protected area in Sabah, using spatial capture–recapture models, to be similar to those we report.
We recorded more males than females and individual males were photographed almost four times more often. Similarly, in other pantherine species male cats move over larger areas than females and are more readily detected by camera traps (e.g. Sollmann et al., Reference Sollmann, Furtado, Gardner, Hofer, Jácomo, Tôrres and Silveira2011). Furthermore, female clouded leopards may also spend more time in trees than males: because of their lower weight they may be more agile and better climbers. Thus females may be detected less often in ground-based camera traps. Unfortunately, data were too sparse to model movement and detection separately for both sexes. Assuming a lower probability of detection and smaller movements for females, the density estimates could increase if differences in movements by sex were accounted for (Sollmann et al., Reference Sollmann, Furtado, Gardner, Hofer, Jácomo, Tôrres and Silveira2011).
Even for a large cat a density of c. 1 individual per 100 km2 is low. This raises several conservation concerns. Firstly, studying clouded leopards is difficult. Even our large and narrow-spaced trapping grid did not result in sufficient data to incorporate potentially important covariates in the model (sex, on/off road camera placement), and density estimates had large standard errors, as did the estimates of Brodie & Giordano (Reference Brodie and Giordano2012). To obtain more precise estimates and to monitor populations over time, substantial field monitoring will be required. Secondly, areas to protect clouded leopards need to be very large to ensure stable populations, probably several thousand km2. Consequently, few existing protected areas are large enough to ensure the long-term viability of any Sunda clouded leopard population. Thirdly, therefore, forests beyond the borders of protected areas, such as timber concessions, are important for the survival of the species. This requires sustainable management of these areas because it is unknown to what extent Sunda clouded leopards can survive in highly degraded forests. The danger is that overexploited non-certified forests are often transformed into oil palm plantations, which are unlikely to support resident clouded leopards. To address these concerns further research should focus on timber concessions, to foster our understanding of habitat and space requirements of the Sunda clouded leopard in these altered environments.
Acknowledgements
We thank KTS-Plantations Sdn. Bhd., Dr Henry Bernard (University Malaysia Sabah) and Raymond Alfred (formerly WWF–Malaysia) for their support and the Sabah Economic Planning Unit for issuing a research permit. Financial support was provided by Point Defiance Zoo and Aquarium, Clouded Leopard Project, WWF–Germany, WWF–Malaysia, Cleveland Metroparks, Minnesota Zoo, Houston Zoo, Nashville Zoo, Panthera Foundation, Zoological Society for the Conservation of Species and Populations and the Leibniz Institute for Zoo and Wildlife Research.
Appendices 1–2
The appendices for this article are available online at Http://journals.cambridge.org
Biographical sketches
Andreas Wilting focuses on the ecology and evolution of mammals, particularly carnivores, in the Sunda Shelf. He integrates results from field research, distribution modelling and phylogenetic approaches to gain insights into the past and current distribution of these species. He and all the other authors have a strong interest in the study and conservation of tropical carnivores, their prey, and their interactions with humans. Azlan Mohamed uses camera-trapping to study mammals in Peninsular Malaysia. Laurentius N. Ambu's main interest lies in sustainable management of wildlife populations, and Peter Lagan and Sam Mannan focus on sustainable forest management in Sabah, Malaysian Borneo. Heribert Hofer's focus is on integrating the behavioural ecology, conservation biology and evolutionary epidemiology of predominantly African carnivores and their pathogens. Rahel Sollmann studies carnivores in Brazil and is currently working on the application of non-invasive methods and related statistical modelling techniques to this group of species.