Assessing Patterns of Human-Wildlife Conflicts and Compensation around a Central Indian Protected Area

Mitigating crop and livestock loss to wildlife and improving compensation distribution are important for conservation efforts in landscapes where people and wildlife co-occur outside protected areas. The lack of rigorously collected spatial data poses a challenge to management efforts to minimize loss and mitigate conflicts. We surveyed 735 households from 347 villages in a 5154 km2 area surrounding Kanha Tiger Reserve in India. We modeled self-reported household crop and livestock loss as a function of agricultural, demographic and environmental factors, and mitigation measures. We also modeled self-reported compensation received by households as a function of demographic factors, conflict type, reporting to authorities, and wildlife species involved. Seventy-three percent of households reported crop loss and 33% livestock loss in the previous year, but less than 8% reported human injury or death. Crop loss was associated with greater number of cropping months per year and proximity to the park. Livestock loss was associated with grazing animals inside the park and proximity to the park. Among mitigation measures only use of protective physical structures were associated with reduced livestock loss. Compensation distribution was more likely for tiger related incidents, and households reporting loss and located in the buffer. Average estimated probability of crop loss was 0.93 and livestock loss was 0.60 for surveyed households. Estimated crop and livestock loss and compensation distribution were higher for households located inside the buffer. Our approach modeled conflict data to aid managers in identifying potential conflict hotspots, influential factors, and spatially maps risk probability of crop and livestock loss. This approach could help focus allocation of conservation efforts and funds directed at conflict prevention and mitigation where high densities of people and wildlife co-occur.


Introduction
Reducing conflict between wildlife and people is considered a top conservation priority, particularly in landscapes where high densities of people and wildlife co-occur [1][2]. Increased visibility for conflict incidents may be attributed to actual increase in incidents taking place or just greater reporting by affected local people [3]. Dearth of knowledge about conflict loss and compensation distribution contributes to poor allocation of conservation resources [4][5]. Failure to address emerging issues with conflict losses and distribution of compensation may lead to escalation of tensions between people and wildlife, and promote retaliatory actions leading to extirpations of species [6][7]. Preventing conflict and improving distribution of compensation are important to fostering co-existence in landscapes that surround protected areas and function as critical buffers for wildlife [8][9].
Indian protected areas (PAs) support a huge array of wildlife that are prone to conflict with people. People tolerate some species such as Nilgai Boselaphus tragocamelus, Chinkara Gazzella bennetti and Blackbuck Antilope cervicapra but are less tolerant of other species such as wild pigs Sus scorfa and elephants Elephas maximus [10][11]. Crop loss is more common than livestock loss, human injury and death [12]. Local residents most often directly bear the costs of living alongside wildlife and may have limited ability to cope with losses [12]. Understanding the factors associated with conflict and where they are likely to occur is important for conservation management of conflicts [2,13].
In this study, we assess and map perceived conflict and compensation distribution to households around Kanha National Park in Central India. Loss refers to probabilities of surveyed households reporting crop raiding or livestock predation. By spatially modeling and mapping perceived conflict and selfreported compensation distribution, we potentially identify household characteristics and practices along with environmental factors that influence conflict and compensation. Data limitations on the timing and locations of actual conflicts and lack of access to detailed records of compensation restrict us to surveys of local people as the primary source of information (see Methods).
In this paper, we examined: a) What agricultural and demographic factors (e.g., land area, total number of crops grown in a year, cropping months in a year, household members, proportion of men), environmental factors (proximity to PA and water, elevation, forest cover), and mitigation measures (e.g. fencing, lighting, physical structures, night guards) are associated with households' self-reported conflict? b) How does loss vary between crop raiding versus livestock predation? c) Which demographic factors, conflict type, wildlife species, household location and reporting effort are associated with reported compensation received by households?
We expected households in closer proximity to the PA and water, surrounded by forest cover and in lower elevations with greater availability of forage from crops (numbers or types of crops grown in a year by a household and growing season) to be more prone to conflict [4,[14][15]. In contrast, we might expect demographic factors (more household members, more men), household characteristics (such as land size or agricultural income) and use of mitigation measures to lower household loss [16]. We also examined reported compensation distribution around Kanha. We expected households reporting conflict incidents, as well as those with more educated members, larger land size, and more livestock to receive better compensation [14]. We also expected compensation distribution to be influenced by species involved, perhaps better for tigers and leopards compared to ungulates and other carnivores [17][18][19].

Study Sites
Kanha National Park (22u 79-22u 279N, 80u 269-81u 39E), established in 1955, is one of India's most well known Tiger Reserves. Kanha (from here onwards) covers an area of 940 km 2 and is surrounded by an administrative buffer of 1005 km 2 (Fig. 1). Vegetation is comprised of Sal and mixed deciduous forests interspersed with grasslands, and supports carnivores such as tiger Panthera tigris, leopard Panthera pardus and wild dog Cuon alpinus, and herbivores including sambar Cervus unicolor, chital Cervus axis, barasingha Cervus duvaucelii and gaur Bos gaurus [9]. Local livelihood activities such as grazing and collection of resources are common occurrences in Kanha although all of these activities are legally prohibited inside Kanha [12]. Human population densities range between 182-195/km 2 , and livestock densities range between 65-79/km 2 in the three districts adjoining the PA [9]. Unlike most other Indian PAs, Kanha has an administratively designated buffer and is surrounded by forested patches that are interspersed with agricultural and barren land [9]. In the administrative buffer activities such as grazing, collection of forest resources and all activities are permitted. The presence of a buffer provides an opportunity to compare conflict experienced and compensation effectiveness among households within and outside the buffer.

Social Surveys and Interviews
Trained assistants conducted 735 structured and open-ended surveys with households selected from 347 villages surrounding Kanha (Table 1). Although there are households inside the park, we only surveyed households outside the park including within the administrative buffer. To enable systematic spatial sampling, we used a grid based approach to sample 5154 km 2 area by placing 218 grid cells (13 km 2 in size) in a 20 km radius outside the PA (Fig. 1). We selected a 20 km radius as a reasonable distance within which we would expect wildlife, particularly ungulates, from the PA to travel outside [20,Karanth et al. unpublished]. Due to the interest in a large geographic area, we sampled 50% of the grid cells in a checkerboard pattern (area of 2319 km 2 ), with 193 cells (we excluded 9 cells that were forested without villages, and 16 cells that were logistically inaccessible). Survey of India topographic maps, imagery and Google Earth were used to digitize 632 villages around Kanha (v6.1). Our goal was to sample at least 60% of villages in a cell and survey at least one household from each village. In each village we opportunistically selected and approached individual households. In some cells that had just 1 or 2 villages, we surveyed more than one household in such villages but ensured that these households were located as spatially apart as possible. We were able to opportunistically sample an average of 4 villages in each grid cell (ranging from 1-5 villages in a cell). We approached adult male and female respondents and they were questioned about household demographic and socio-economic characteristics and mitigation measures (use of fencing, lighting, guard animals etc., Table 1). Respondents were specifically questioned about recent (2010 and 2011) experience with all conflict incidents (crop raiding, livestock predation and human injury or death), as well as compensation reported and received. All surveys were conducted in October 2011 in Hindi and recorded responses were translated to English. All survey protocols were approved by Columbia University's IRB.
We also collected official year wise summary records of total compensation paid out to individual households by the authorities for crop and livestock loss for the years 2009-2011. These records only cover households located inside the PA and administrative buffer. Compensation by the park authorities is determined by the geographic location of households. Households apply to either territorial revenue or forest divisions in which they are located and the authorities compensate those households filing within their administrative division (Karanth et al. unpublished). Typically the amount of compensation provided is determined by officials assessing the extent of damage and do not differentiate among species (Shukla pers comm. 2012) and do not differentiate among species.

Variable and Model Selection
We used model selection to identify factors associated with selfreported household incidents of crop raiding, livestock predation and compensation distribution around Kanha. The variables chosen are in Table 2. To avoid colinearity of variables used in the model, we computed Pearsons' correlation coefficients for all pairs of variables. High correlation (.0.5) resulted in elimination of some variables [used distance to PA variable over distance to closest forest patch, (Table S1A-B) [21][22]. To permit direct comparison of estimated coefficient values, we scaled the regression inputs by dividing by two standard deviations [23].
Our models represented different hypotheses about factors potentially influencing reported household loss. We defined an a priori set of candidate models, assessed model fit and identified the variables significantly associated with crop and livestock loss using the Corrected Akaike's Information Criterion (AICc) as the model selection tool since the ratio of the number of data points to the number of structural parameters in our model set was, on average less than 40 [24] (Table S2). Additionally, the AICc criterion is also known to perform well under some degree of correlation among explanatory variables [24]. AICc weights represented the relative measure of appropriateness of a given model relative to the entire model set and rank the models (including a 'global model' that includes all explanatory variables of interest). The top model which provided the most optimal tradeoff between model parsimony and fit, in our goal to assess association, and closely competing top models (a subset of models accounting for .95% of the model uncertainty) was used to estimate probabilities of crop raiding, livestock predation and compensation distribution by multi-model inference via model averaging [24].

Modeling Household Risk of Crop-raiding and Livestock Predation
We defined household loss as a household reporting recent experiences with crop raiding or livestock predation (in 2010 and 2011). We chose to ask respondents only about the most recent year as we believed that respondents might forget or confuse incidents over multiple years. Since our responses were binary, we fit logistic regression models to model household conflict loss.
We modeled self-reported household crop raiding and livestock predation as a function of environmental factors. The environmental factors included distance to PA, distance to water, elevation and percentage of forest cover within 3 km of household for each household (data derived based on 9). The selection of 3 km buffer around a household represented a reasonable estimate for home ranges of crop raiding wildlife in tropical forests ( [20] and references therein). We also examined the influence of agricultural factors (number of crop types grown in a year, number of cropping months, land area owned by households reported in the surveys). We examined demographic characteristics (number of household members, proportion of men, gender of respondent). We also compared household use of different mitigation efforts: lighting versus fencing versus guard animals. We modeled reported compensation received by households as a function of individual respondents factors (gender, age), household characteristics (land size, total livestock owned), conflict type (crop loss, livestock loss, human injury, human death), reporting effort by households to authorities, wildlife species, proximity to the PA and location inside administrative buffer.
We constructed 31 a priori models representing our hypotheses about causes of 'crop-raiding' (Table S2), including a global model with 18 explanatory variables. We modeled household livestock predation as a function of 20 variables represented by 30 models (Table S2). Similarly, we proposed 36 models of compensation distribution using 22 variables (Table S3). We computed the weighted frequency of each variable among all the models to assess the relative importance of each variable. For each surveyed household, we estimated crop loss, livestock loss and compensation access probabilities for all the candidate models. The weighted average of these individual model estimates, weighted by the Akaike weights [24], was used to generate overall estimates for individual households crop loss, livestock loss and compensation access.

Spatial Mapping of Conflict Risk and Compensation Distribution
We derived estimated probabilities for crop loss, livestock loss and compensation access for every household from each individual model from the entire model set and these values are then weighted to derive model-averaged estimates for each household. We used ordinary kriging to estimate crop loss, livestock loss and compensation distribution probabilities in non-sampled areas. Kriging is a geostatistical interpolation technique that permits inferences about a parameter of interest (in our case, the probability of conflict or compensation) in non-sampled areas by using information available in sampled locations and accounting for uncertainty as the distance between spatial locations increases. The approach produces a semivariogram that describes the spatial correlation between the points. Several models are available for fitting a semivariogram, such as spherical, circular, exponential, Gaussian and linear [25][26]. Based on the Akaike's information criterion (AIC), we fit the spherical model to the observed data, using the Kriging Interpolator 3.2 extension in ArcView Spatial Analyst [27]. In our application, we treated the weighted estimate of conflict loss probability at each sampled location as ''truth'', assuming no variability around the estimate (i.e. no nugget effect) [28], to generate maps of conflict loss probabilities for crop loss and livestock predation, as well as distribution of compensation reported by households around Kanha.

Household Characteristics, Agricultural Practices and Reported Crop Loss
Household respondents were largely men (87%), with many (77%) having completed 8 th grade or less education. Majority of households (93%) were engaged in agriculture, 88% legally owned land with average land size of 4.3 acres. Households reported growing 12 different crops. The major crops grown were rice, millet and legumes, and the average length of time for crops in the field was 5 months (range 4-12 months, Table 1).
Households reporting crop loss listed 17 species as crop raiders including 10 herbivores, 4 carnivores, 2 primates and peacocks. Animals reported causing the most crop damage were wild pig Sus  Fig. 2). The highest numbers of raiding incidents were reported from September to December, peaking in October. Sixty-four percent of households reported experiencing more than five incidents per year and 32% households reported 2-5 incidents per year. From modeling factors associated with self-reported household crop loss, we found five top-ranked models with substantial weights (cumulative AICc .0.95, Table 3). As predicted, the total number of crops grown in a year was associated with increased households' loss (positive b coefficient in Table 3) but no other household characteristics were relevant. Although use of mitigation measures, in general, appeared important, no individual measure, such as fencing or lighting, appeared noticeably more important (Table 3). We suspect that the strong positive associations of any general mitigation measure (as indicated by positive beta coefficients in Table 3) are a result of these mitigation measures being put in place because households have historically experienced many incidents of crop raiding. Among the environmental factors distance to PA was associated with decreased crop loss (negative beta coefficient in Table 3). Distance to water was associated with increased Table 2. Details on variables collected from surveys and used in models.

Variable
Crop loss Livestock loss Compensation household loss (positive beta coefficient in Table 4) but high uncertainty in the estimates suggested that this parameter must be interpreted with caution. Contrary to our expectations other modeled factors did not appear to influence crop loss (Table 2  and Table S2), although we recognize that a much larger sample size might have assisted in providing greater support for other explanatory factors as well.
We mapped the probability of crop loss for households around Kanha (Fig. 3). Our modeling indicated that crop loss risk was high in most places, average estimated probabilities of crop loss was 0.93 (range 0.66-0.99). Probabilities of crop loss averaged for households inside the administrative buffer was 0.95 (S.E = 0.004, range 0.77-0.99) and for households outside the buffer was 0.92 (S.E = 0.003, range 0.66-0.99). Spatial modeling based on ordinary kriging (without the nugget effect) suggests that households located closer to the PA have higher risk regardless of location within or outside the administrative buffer. The spatial model locates hot spots of conflict (Fig. 3) and identified high risk particularly around some villages in the buffer (Ramepur, Moharai, Kirsari), and for other villages outside the buffer (Baila, Budhanwara, Maharajpur, Khorja).

Livestock Ownership and Perceived Predation Risk
Most (91%) of households owned livestock and ninety percent of animals were stall fed or locally grazed compared to 10% being park grazed. Noticeably households reported an increase in stallfeeding by 7.6% over 10 years and an increase in the number of days their animals were grazed in the last ten years (Table 1). This is perhaps indicative of declining forage areas in close proximity to people's households requiring them to graze livestock further away from home and well inside the PA.
Livestock losses to ten carnivores were reported and most troublesome species were jackal Canis aureus, wolf Canis lupus, tiger and leopard (Table 1). Livestock losses were distributed across the year but higher from April through July, which are the driest summer months in Kanha (9). Some (34%) households reported more than five incidents per year, 34% households reported 2-5 incidents per year and 31% reported one incident per year. Sevenpercent of households reported experiencing injury and 1% reported instances of human death (Table 1). From modeling factors associate with self-reported household risk to livestock loss outside Kanha, we identified nine top-ranked models (cumulative AICc .0.95, Table 4). These models indicate that as predicted grazing cows inside Kanha was associated with increased loss (positive beta coefficient in Table 4). Among the environmental factors, distance to PA was associated with decreased livestock loss (negative beta coefficient in Table 4). Distance to water was associated with increased household loss (positive beta coefficient in Table 4) but high standard errors suggest caution with interpreting this result. Two mitigation factors night watching and physical structures appeared in top ranked models (Table 3). Only use of physical structures (negative beta coefficient in Table 4) was associated with lower losses but high standard errors suggest that this must be interpreted with caution. Other modeled factors did not appear to be associated with household livestock loss (Table S2).
We mapped risk of livestock predation for households around Kanha (Fig. 4). Our modeling clearly indicates that livestock loss

Mitigation Measures
Households reported several mitigation measures used by them to protect property and lives. To protect crops the most common mitigation measures were night watching (67%) and scare devices (49%) and to protect livestock measures deployed by households were closer watch on animals (48%) and physical structures (47%, Table 1, Fig. 2). Our models suggest that potentially physical structures and guard animals appeared to be associated with lower crop loss, and no mitigation appeared to be associated lower livestock loss.

Compensation Distribution
Although 73% of those surveyed reported that they experienced crop loss, only 26% of these households reported losses to authorities and 22% of those reporting received compensation  Thirty percent of surveyed households reported to have experienced livestock loss, but only 34% of these affected households reported losses to authorities and 41% received compensation upon reporting. In Kanha official compensation paid from 2006-2011 to 524 households inside the PA was $84 and higher than the $63 compensation paid to 1644 households located in the buffer (Kanha Forest Department 2012). In our study, households reported receiving compensation averaging $22 (ranging $20-100). Although we cannot directly compare official compensation paid to individual households, the financial amounts of compensation from our surveyed households compare with official records. Surveys conducted in parks around Karnataka also find similar self-reporting levels of loss and compensation ( [12], Karanth et al. unpublished).
From modeling factors associated with reported household access to compensation we identify important factors associated with compensation receipt (Table 5). Households filing claims to authorities and located in the administrative buffer were more likely to receive compensation and this is supported by positive beta coefficients for reporting crop and livestock loss (Table 5). Species involved was also associated with compensation received by households, particularly reporting tigers had positive influence (positive beta coefficient in Table 5). Conflict incidents related to other species do not appear to have yielded compensatory benefits to households (negative beta coefficients in Table 5). Reporting incidents of human injury and death occur in fewer models but did not appear to be associated with improved compensation (positive beta coefficients in Table 5). Households' location in the legal buffer appeared to be associated with higher compensation (positive beta coefficients in Table 5). Other modeled factors did not appear to influence compensation (Table S3).
We compared losses experienced and reported to authorities, and compensation received by households inside and outside the buffer. Households reporting crop loss was similar inside and outside the buffer but households reporting livestock loss was higher inside the buffer (Table 6). We find that households inside the administrative buffer were more likely to report losses (46% for crop loss and 62% for livestock loss) compared to households outside the buffer (19% for crop loss and 18% for livestock loss). Of the households that reported loss to authorities, households inside the buffer were more likely to receive compensation (35% for crop loss and 48% for livestock loss) compared to households outside (11% for crop loss and 29% for livestock loss).
We mapped reported compensation distribution for households around Kanha. The estimated probabilities averaged 0.09 (range 0.01-0.72, Fig. 5). Households located within the administrative buffer averaged 0.21 (S. E = 0.02, range 0.03-0.72) compared to households outside 0.05 (S. E = 0.003, range 0.01-0.45). Compensation distribution for households located in the administrative buffer is higher than households located outside indicating some positive influence of management inside the buffer. Our spatial model suggests higher compensation around Moharai, Ramepur, Kirsari, Lamna, Umaerdehi and Chilpi inside the buffer (Fig. 5).

Discussion
Our results identified and quantified main factors associated with self-reported household crop raiding and livestock predation losses, as well as compensation distribution (Figures 3, 4, and 5), in 5145 km 2 area surrounding Kanha. Similar to other studies, we find spatial and temporal variation in conflict losses [4,14,18]. Overall, the probability of crop loss was high (average 0.93) for households and higher for households located within the administrative buffer. Crop loss was associated with greater cropping months and proximity to Kanha. No individual mitigation measure was associated with lowering crop loss. Overall, the probability of livestock loss was comparatively lower compared to crop loss (average 0.60) and higher for households located in the administrative buffer compared to households located outside. Livestock loss was associated with households grazing cows inside Kanha and located closer to the reserve. However, use of physical structures to protect livestock appeared to be the only mitigation measure associated with lowering of livestock loss. Interestingly, over the last ten years households reported a gradual increase in stall-feeding practices compared to free range grazing, suggesting the gradual modification of behavior. This result, however, may merely be a reflection of the inability of people to recall practices from the past. Our findings mirror other studies, where crop and livestock loss are a function of multiple factors, with individual choices and behavior influencing loss and often difficult to modify [29].
Considerable conservation attention and monies are currently being invested in mitigation factors that have rarely been evaluated for effectiveness [16,[30][31]. We find only use of physical structures was potentially associated with lowering Although 73% of households reported crop loss and 33% reported livestock loss, only 26% reported crop loss and 34% reported livestock loss to authorities. Among reporting households, 22% received compensation for crop loss and 41% received compensation for livestock loss.
Likelihood of compensation distribution received by households was higher for households that report incidents related to tigers (although most damage was reported by pig, jackal and wolf incidents, Table 1) and households located in the buffer. There is considerable debate about people's ability to accurately report losses to authorities, and other studies have found mismatches to official records [2,14,32]. These reported losses determined by surveying people have to be complemented with field-based monitoring of households over time. We do find among surveyed households in the administrative buffer were twice as likely to report crop loss and thrice as likely to report livestock loss compared to outside (Table 6). Households in the administrative buffer were more likely to receive compensation from authorities, suggesting that the legal designation and management of an administrative buffer improves reporting of losses and receipt of compensation by households (Table 6).
Overall, many (87%) households around Kanha report experiencing some kind of conflict incident with wildlife and this is comparatively higher than other places in the world [2,33].
However, there appears to be some degree of acceptance among people that there will be some losses to wildlife (Karanth pers. obs 2011, Karanth et al. unpublished), and perhaps households have evolved coping strategies [11][12]17,34].
Unlike most Indian PAs, Kanha has a legal administratively designated management buffer [9]. We estimated higher risk of crop and livestock loss inside the buffer, and find higher compensation in the buffer. This corroborates with other studies that suggest delineation of buffers or other management strategies should be based on ecological and economic realities [9]. Additionally, in this landscape we found household distance to the PA was more relevant for livestock and crop loss compared to distance to water. Our results potentially suggest that some crop losses to wildlife from the PA but we know that wildlife naturally persists outside the PA and have overlapping resource use with people and domestic animals (forage and water, Karanth et al. unpublished). So this raises important questions about who is responsible for compensating local people for wildlife damage outside the jurisdiction of PAs and park authorities [6,35]. Our spatial risk maps (Figs. 3 and 4) clearly identified locations that are hotspots for crop and livestock loss and can guide managers in improving existing conflict prevention and mitigation practices.

Management Implications
Most conflict studies are characterized by poor spatial sampling and modeling [36][37]. To improve efficacy and efficiency of conservation actions, managers require surveying and modeling approaches that are spatially explicit and rigorous [5,7]. Our approach identified factors associated with crop and livestock loss around Kanha National Park in India. We were able to link ground based survey data to landscape level data and spatially map conflict crop and livestock loss, as well as compensation distribution to identify conflict hotspots and focus ground based conservation efforts. For example, villages and households located in high-risk areas can be better educated about how to prevent and mitigate conflict and can report losses to authorities in an organized manner.
A monitoring system that systematically records and disperses information on conflict is much needed. Such an approach can guide the development of a risk database, and develop live warning and monitoring systems (such as those cell phone alert systems currently in use for Elephants in Valparai in southern India). Experiences from Uganda, Kenya and Sumatra suggest that establishing and long-term maintenance of monitoring systems in local communities remains a considerable challenge and is often difficult to sustain [15,30,38]. However, greater participation and involvement of multiple stakeholders (local households, PA authorities) can considerably decrease hostility towards wildlife species [29, 35,39].
Kanha is atypical of most Indian PAs because the administrative buffer gives it a lesser hard edge than other PAs and livelihoods of people living in the buffer fall within the purview of park management [9]. As consideration is currently under way for designating buffers around other protected areas in India, our results from Kanha suggest that compensation is more likely to be distributed to those who suffer losses if the administrative buffer is designated to include more susceptible locations for crop and livestock loss. This conclusion may apply to other parks located in human-dominated landscapes, though cultural and ecological factors will vary.
Our study has some limitations. The analysis relied primarily on individuals' recall and self-reporting obtained through surveys and did not observationally monitor actual incidents of conflict. This raises possibilities for over-reporting or exaggeration of losses, although there are no indications of a systematic bias across so many households. The number of surveys is large relative to most published studies and field-based human-wildlife conflict surveys that reduced the potential biases (although we recognize the limitations of self-reported data). To complement the surveys and provide objective data on human wildlife conflicts, active location based monitoring of households to measure conflict as it happens along with assessment of mitigation measures and compensation distribution process is much needed. Extension of surveys across multiple years is also needed to assess additional factors that might influence conflict risk that may not be identified in a single year assessment. Overall, our approach for identifying associated factors and spatially mapping risk is easily extendable to other landscapes where people and wildlife co-occur, so that allocation of physical efforts and funds are more effectively targeted at conflict prevention and mitigation.