Abstract Protected areas are extremely important for the long term viability of biodiversity in a densely populated country like India where land is a scarce resource. However, protected areas cover only 5% of the land area in India and in the case of large carnivores that range widely, human use landscapes will function as important habitats required for gene flow to occur between protected areas. In this study, we used photographic capture recapture analysis to assess the density of large carnivores in a human-dominated agricultural landscape with density >300 people/km2 in western Maharashtra, India. We found evidence of a wide suite of wild carnivores inhabiting a cropland landscape devoid of wilderness and wild herbivore prey. Furthermore, the large carnivores; leopard (Panthera pardus) and striped hyaena (Hyaena hyaena) occurred at relatively high density of 4.8±1.2 (sd) adults/100 km2 and 5.03±1.3 (sd) adults/100 km2 respectively. This situation has never been reported before where 10 large carnivores/100 km2 are sharing space with dense human populations in a completely modified landscape. Human attacks by leopards were rare despite a potentially volatile situation considering that the leopard has been involved in serious conflict, including human deaths in adjoining areas. The results of our work push the frontiers of our understanding of the adaptability of both, humans and wildlife to each other’s presence. The results also highlight the urgent need to shift from a PA centric to a landscape level conservation approach, where issues are more complex, and the potential for conflict is also very high. It also highlights the need for a serious rethink of conservation policy, law and practice where the current management focus is restricted to wildlife inside Protected Areas.

Citation: Athreya V, Odden M, Linnell JDC, Krishnaswamy J, Karanth U (2013) Big Cats in Our Backyards: Persistence of Large Carnivores in a Human Dominated Landscape in India. PLoS ONE 8(3): e57872. https://doi.org/10.1371/journal.pone.0057872 Editor: Matt Hayward, Australian Wildlife Conservancy, Australia Received: September 4, 2012; Accepted: January 30, 2013; Published: March 6, 2013 Copyright: © 2013 Athreya et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: The project was funded by a Kaplan Graduate Award (http://www.panthera.org/); a Norwegian Embassy and Norwegian Research Council grant (http://www.nina.no/ninaenglish/Researchprojects/HumanWildlifeConflicts.aspx). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.

Introduction Charismatic predatory species have long held a central place in global conservation strategies, both in terms of attracting the public’s attention and serving as a focus for research and conservation effort [1]. Accordingly, they are frequently cast in the roles of flagships and umbrellas [2]–[4]. This is based on symbolic (their flagship role) and functional (their umbrella role as mediators of top-down cascade effects) perception of their role in ecosystem processes and their presumed dependence on wild nature and wilderness. A consequence of this is that strategic planning for their conservation is often based around protected areas that by design largely exclude or minimize inclusion of populated areas with human settlements and agricultural land-use. The protected area focus for large predator conservation has been a powerful argument for justifying the setting aside of wilderness areas, especially in tropical developing countries [5] where there has been a long standing skepticism about the ability of wildlife to persist in unprotected landscapes with moderate to high human densities [6]. India is no exception to this pattern, and most of its conservation focus, in terms of conservation actions, research, and legislation, is focused on protected areas which aim to minimize human settlements and agro-pastoral land-use. In some cases, such as tiger (Panthera tigris) conservation, this priority given to protected areas may well be justified [7], [8]. Some ecological studies, albeit in areas with sparse human density, have documented that large predators in general [9]–[11] benefit from areas with high densities of their natural, wild prey. Although there is also some evidence of carnivores that adapt to a wide variety of habitats even with a long history of human impacts and influences [12]. What is important for a pragmatic approach to the conservation of these wide ranging carnivores is to identify the limits of their tolerance, as much as the nature of their preferences, and these can be assessed only when the large carnivores occur outside protected areas, where their interface with humans is high. This knowledge is crucial when planning for long term conservation objectives which will need to integrate these species into the wider landscapes matrices where protected areas [8], [13] can be connected via human use landscapes. During the last twenty years, there has been an increasing awareness of the ability of some large predators like wolves (Canis lupus) and pumas (Puma concolor) to live in very human-dominated, even sub-urban, environments in the developed world [14], [15]. Accordingly there has been a great deal of research aimed at these situations and appropriate management responses have been developed [15]. Most examples of large carnivores in urban landscapes are from countries with a low human population density although in the case of mountain lions in Southern California, it has been seen that part of their home ranges overlap with densely populated urban landscapes [15]. Recent results from Africa show comparable densities for lions, cheetahs and jackals inside and outside the protected areas, although the human population density there was low [12]. Very little is known about the ecology of ‘urban’ carnivores in densely populated countries where the potential for conflict can be very high. For instance, India has a high diversity of large carnivores, many of which share spaces with one of the highest human and livestock populations in the world [16]–[18]. The most common is the leopard (Panthera pardus), which is frequently reported from many human dominated landscapes across India where it is involved in a wide range of conflict situations, often with fatal outcomes for humans. The main management response, in the absence of robust information on large carnivore ecology in human dominated areas has been their translocation to nearby protected areas [19] based on a belief that these leopards are “stray” individuals that have dispersed from protected forest areas and need to be “helped” back to the forests. Recently, this has been documented to worsen the situation, leading to increased attacks on people near the sites of release [18]. Clearly there is a need to document the status of leopards living in these human-dominated landscapes, as well as to understand the community structure of wild carnivores in highly modified ecosystems that are also home to high density of humans. In this paper, we present evidence of the presence of an entire community of wild carnivores that share space with very high densities of humans. Although the work focused on estimating the density of the leopard populations, we also provide a density estimate of the striped hyaena (Hyaena hyaena) and an overview of the occurrence of other carnivore species. Using the results of this work, we discuss the potential importance of broadening the current protected area conservation focus to include a consideration of the value of human-dominated landscapes [20].

Methods Study Area The study was conducted in a densely populated, irrigated valley in Akole Tahsil located in the Ahmednagar district (19.576959 N 73.937123 E to 19.460715 N 74.089954 E) of western Maharashtra, India. An Indian district is administratively comparable to a ‘county’ and the Tahsil is a sub-unit of it. Approximately 80% of the human population is rural with farming of sugar cane, millets, and vegetables being the major source of livelihood. Rainfall varies from 1000 to 2000 mm per year. Akole Tahasil contains 191 villages (as per the 2001 census; http://ahmednagar.gov.in/html_docs/GEO-Main.htm accessed 26th June 2012) with a human population density of 177 km−2. The intensive study area covered a 179 km2 area in the irrigated valley around Akole town. Digitized maps of all the households in the study area, combined with household interviews were used to determine the average number of people per household (Athreya et al. unpublished) which was 357 humans/km2 in the intensively cultivated area. Only 15% (1515 km2) of Akole Tahasil is protected (legally designated as government-owned Reserved Forest; [21]), and mainly occurs in its western part or in small, dispersed patches unsuitable for cultivation (see Figure 1). The only wildlife protected area is the Kalsubai Harishchandragarh Wildlife Sanctuary (299 km2), which is located on the western boundary of Akole Tahasil, 18 km beyond the edge of our intensive study area (Figure 2). PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Figure 1. Map of the study area which consisted of an irrigated valley around the town of Akole in the Ahmednagar District, Maharashtra, India. https://doi.org/10.1371/journal.pone.0057872.g001 PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Figure 2. Map of the Ahmednagar district with the study area polygon and the nearest protected area of Kalsubai Harishchandragarh Wildlife Sanctuary. https://doi.org/10.1371/journal.pone.0057872.g002 Ethics Statement Relevant permissions to carry out the ecological research were obtained from the Office of the Chief Wildlife Warden, Maharashtra Forest Department. Camera Trap Surveys Forty camera trap locations were selected following study design approaches prescribed for large felids [22] at sites that had intensive signs of their usage. Most trap sites were on human trails. The duration of the camera trap survey was fixed at 30 days between November 2008 and December 2008 to meet the underlying assumption of population closure [23]. The closure test in program CAPTURE was used to test the assumption [24]–[26]. The sampling was carried out in two trapping sessions of 15 days each covering two blocks, each with 20 trap locations (Figure 1). The capture history data from the two blocks were combined and analyzed using program CAPTURE. Of the 40 trap locations, three became non-operational due to theft of cameras. Based on typical minimum home range sizes of leopards [27], [28], an average trap spacing of about 1.5 km was maintained to ensure that all individual leopards were potentially exposed to trapping. Each camera trap had two camera units facing each other 5 m apart, fixed at locations judged to be optimal leopard travel routes. Deercam DC 300 passive detection camera traps protected by metal shells were used and set at a height of 40–45 cm above the ground. Camera traps were checked twice daily to record the number of exposures, to change the films or batteries if needed, and to switch them on for night operation (around 6 pm to 7 am). The cameras were turned off during daytime because of the large amount of human and livestock traffic. Individual leopards and hyaenas were identified from photos based on pelage marking patterns. However, in some cases we could not obtain photographs of both flanks of hyaenas, and therefore we used only the right flank to identify them. We provide an index of trap visitation rates for other wild carnivores based on the maximum number of images obtained per species per trap site in all 37 traps during the survey period. In the absence of an appropriately designed survey for all the other species, and adequate samples of detections we could not use abundance models of occupancy [29]–[31]. Estimating the Density of Leopards and Hyenas using Capture-recapture Sampling Images of leopards more than a year old were used for estimating abundances [32] because the probability of capturing younger animals is very low (due to the time delay between first and second images if a family group pass a camera), and because the high mortality rate of this age class could lead to violation of population closure assumptions. There are two different approaches to estimation of animal density from photographic capture-recapture data. The earlier approach, hereafter referred to as conventional capture recapture (CR) model [33], [34] relies on capture-frequency data derived from photographed animals to estimate population size (Appendix 1), and thereafter uses distances between locations of recaptures to obtain the effective area sampled for estimating density. However, recently this approach has been complimented by spatially explicit capture-recapture (SECR) models that incorporates information on capture location directly into the modeling and the estimation of capture probabilities [35], [36]. However, for the purpose of comparing our results with earlier studies (e.g. [27], [37], [38]), we present estimates based on both CR and SECR models. For estimating abundance with the CR models, we used the program CAPTURE [33], [39]. In CAPTURE, different models are compared based on expected effects of individual heterogeneity in capture probability (Model M h ), behavioural responses to trapping (Model M b ), and changes in capture probability among sampling occasions (Model M t ) [34]. The null model assumes that capture probability is not affected by any of these factors. As suggested by Karanth & Nichols [40], we rejected the null model because the model assumption is likely to be violated due to territorial behaviour and trap response. We estimated the effective sampling area (for density estimates) by adding a circular buffer around each camera trap location equal to half the average mean maximum distance between locations of individuals photo-captured more than once [27], [37]. For more reliable estimation of leopard densities we preferred the recently developed spatially explicit capture-recapture models (SECR). In addition to individual animal capture histories, SECR models also use the spatial information from capture locations in the density estimation process (see [36] for a current summary of SECR models). Of the two types of SECR models available, we preferred the Bayesian models [35], [41] over the likelihood based models [42] because the former does not rest on asymptotic assumptions, and, can potentially be extended to open model scenarios to estimate survival and other parameters with multi-year data. For leopard density estimation we used program SPACECAP [36], which reports interval estimates of density as direct probabilities without asymptotic assumptions for our relatively small sample of captures [35]. Initially a buffer of 15 km was generated around a grid of equally spaced points (4757), each 0.336 km2 in size, to represent all probable leopard activity centers (see [35], [41] for a full description of the methods). An area of 1589 km2 leopard habitat over which these activity centers S i were uniformly distributed was estimated after deducting an area of 2.67 km2 of water bodies. We used three input data files consisting of animal capture details (which individual was photographed at which site and on which sampling occasion), trap deployment details (which traps were active when and where) and the potential home range or activity centers. We used 52000 iterations, of which the initial 2000 were discarded, a thinning rate was set at 50, and augmentation of 110 individuals were used (see [35], [41] for details of estimation methodology). Information on leopard deaths or captures that occurred at the study site, in the camera trapping period were obtained from the Forest Officials. Unlike most other parts of the world, government sanctioned lethal control is not the norm in India and leopards not killed legally.

Discussion Photographic capture recapture has been used for estimating densities of many secretive large carnivore species, including tigers [34], [43], [44], leopards [27], [38], [44], snow leopards (Panthera uncia) [45], jaguars (Panthera onca) [46] and hyaenas [47]–[49]. The previous conventional analyses of obtaining population density have given way to more robust methods that use the spatial information of location of traps rather than methods that estimate the size of the effective area heurestically [35], [41]. The spatial capture recapture model is also less affected by the small sample sizes often associated with camera trap data from large felids [35], [41]. A range of studies have estimated leopard and hyaena densities using photographic capture-recapture. Balme et al. [27] obtained a density of 7.17 per 100 km−2 in a protected buffer area in South Africa while the non-protected farmlands had leopards at 2.49/100 km−2. A recent study from a protected area in Cambodia, devoid of human habitations, obtained leopard density estimates of 3.6/100 km2 and 3.8/100 km2 using SECR (Spatially Explicit Capture Recapture) and conventional CR (Capture Recapture) methods [50]. Estimate of leopard density in India is available from only within protected areas with 15/100 km2 [38] but this area is largely devoid of people and agricultural land-use. Harihar et al. [44] found that the density of leopards in the Rajaji National Park decreased from 9.76/100 km2 to 2.07/100 km2 with a concurrent increase of tiger populations following relocation of people from within the Park. Recent estimates from camera trapping studies on hyaenas outside protected areas in India report density of 3.67–6.5/100 km2 [48] whereas higher densities of 15.1/100 km2 [47] and 3.9–5.67/100 km2 [49] were obtained from within the protected areas of Sariska Tiger Reserve and Rajaji National Park, India. Prey biomass is seen to strongly influence tiger [51] and leopard [10], [52] density although other factors such as interspecific competition [44] and disease [53] also affect carnivore densities. Leopards in protected areas in India feed on small to medium sized wild prey such as cheetal (Axis axis), sambar (Rusa unicolor) and langur (Semnopithecus spp.) [54], [55]. Our study site contains no other apex predator and no wild ungulate prey species suitable for leopards. Data (unpublished results) indicates that leopards in our study area primarily subsist on a diet of domestic dogs and livestock, which are abundant. Thus in our case, models of predator-prey density need to include estimates of the density of domestic prey species as well as wild prey. We also obtained clear evidence in the camera-trap photos, and from cubs that were found dead or rescued from wells in the study area, that the study population consisted of resident and reproductive individuals and was not made up of only sub-adult animals which had potentially dispersed from some distant patch of forest habitat. These lines of evidence along with information from collared leopards (unpublished results) combine to overturn the popular view that leopards in such human-dominated habitats represent a few occasional “stray” individuals. It is clear that the study area contained a dense, established, breeding population of leopards. The results are also interesting for evolving theories on the ways in which community structure changes following human interventions [4], [56]. Conventional theory would not have predicted the persistence of such a carnivore biased community with eight species of carnivore persisting in a high density human use landscape devoid of wild herbivore species. The leopard is the apex predator here with the tiger and the wolf being absent from the system, although wolves have been reported in the dry hills surrounding the irrigated valley where this study was conducted. This merely underlines how aspects of carnivore behavior that permit adaptation to diverse habitats are often just as crucial as the more conventional life-history metrics such as body size and reproductive rates that are often used to predict species and community persistence [4], [56]. Our results show that large predators, like leopards and striped hyaena, are probably not very suitable as either flagships, umbrellas or indicator species for wild nature in India. It appears that sugarcane and other tall crops, domestic dogs and livestock are sufficient as habitat and prey, respectively, for the leopards outside designated protected areas. In these contexts, leopards can serve as important flagships under an alternative approach in the conservationist’s toolkit. This is the philosophy based on sharing space and integrating wildlife into human-modified landscapes where the focus is as much on knowing the “social carrying capacity”, which is defined by the tolerance of humans towards predators [57], in addition to the ecological carrying capacity. This will require an acceptance of situations where humans and wildlife share multi-use landscapes to the extent possible, outside of protected areas. To achieve this, the research focus must include areas outside protected areas as much as those inside, and social science research as much as ecological research. It also requires enabling flexible and pragmatic legislation acceptable to rural people because it needs to take into account their concerns and interests also. This study shows that leopards can persist in the human-modified landscapes and is possibly dependent mainly on the social tolerance. Despite their subsistence on a diet of domestic animals the levels of conflict in the study area are quite low (unpublished results). Although leopards are implicated in the highest number of fatal attacks by a large felid on humans in some other parts of India, no fatal attacks are known from the study area, even though more than 300 people/km2 share the same space with 5–6 leopards/100 km2. There are indications from other parts of India that this phenomena of carnivores surviving in highly modified landscapes is not confined to leopards. The Asiatic lion (Panthera leo) has recently extended its range and is known to now use even the areas outside the protected area [17] as well as some tiger populations which have recently been reported to use sugarcane areas [58] in northern India. Our study documents for the first time that a whole guild of predators can persist in totally human dominated landscape in India. This probably has a lot to do with India's laws which makes it illegal to kill any wildlife for sport or for consumption. In the case of large cats, even killing man-eaters requires permission from the state authorities, unlike most other countries where even livestock killers are often removed immediately. India has also been known for its tolerance towards other life forms, even the large, potentially dangerous species. It also implies that a much greater area of potential "tolerance habitat" is available outside protected areas, and potentially a far greater degree of connectivity between protected areas than generally expected. The ability of conflict causing species to persist in close proximity to humans greatly expands the spatial extent of human – wildlife interfaces beyond the narrow “zone of influence” [20] that surrounds protected areas. This poses many challenges for India’s legislation and wildlife management structures which are heavily focused on protected areas and are very wildlife-centric. There is a clear need to recognize that these potentially conflict causing species [59] can, and will, colonise many areas and that their management cannot only be based on a hands-off policy. That being said there is a clear need to ensure that management interventions do not make the situations any worse [18]. The results of this study add to an emerging body of empirical results that demonstrate the conservation value of unprotected, human-dominated landscapes for large carnivores [60]–[64] as an important supplement to protected areas. While this approach has long been recognized in temperate areas (e.g. [65], [66]), it is only recently being demonstrated in tropical areas [67], although its general applicability as a conservation model may be highly species and context specific [62] depending on a range of ecological, social, cultural and economic factors.

Acknowledgments We would like to thank the Maharashtra Forest Department for their support. We would also like to thank the Centre for Wildlife Studies and the Asian Nature Conservation Foundation, Bangalore for institutional support and N. Samba Kumar and D. Jathanna for providing support for the camera trapping study and for useful discussions and advice.

Author Contributions Conceived and designed the experiments: VA MO JDCL JK UK. Performed the experiments: VA. Analyzed the data: VA MO. Contributed reagents/materials/analysis tools: JDCL JK UK. Wrote the paper: VA.