Figures
Abstract
Wildlife-human conflicts occur wherever large carnivores overlap human inhabited areas. Conflict mitigation can be facilitated by understanding long-term dynamics and examining sex-structured conflict patterns. Predicting areas with high probability of conflict helps focus management strategies in order to proactively decrease carnivore mortality. We investigated the importance of cougar (Puma concolor) habitat, human landscape characteristics and the combination of habitat and human features on the temporal and spatial patterns of cougar-human conflicts in British Columbia. Conflicts (n = 1,727; 1978–2007) involved similar numbers of male and female cougars with conflict rate decreasing over the past decade. Conflicts were concentrated within the southern part of the province with the most conflicts per unit area occurring on Vancouver Island. For both sexes, the most supported spatial models for the most recent (1998–2007) conflicts contained both human and habitat variables. Conflicts were more likely to occur close to roads, at intermediate elevations and far from the northern edge of the cougar distribution range in British Columbia. Male cougar conflicts were more likely to occur in areas of intermediate human density. Unlike cougar conflicts in other regions, cattle density was not a significant predictor of conflict location. With human populations expanding, conflicts are expected to increase. Conservation tools, such as the maps predicting conflict hotspots from this study, can help focus management efforts to decrease carnivore-human conflict.
Citation: Teichman KJ, Cristescu B, Nielsen SE (2013) Does Sex Matter? Temporal and Spatial Patterns of Cougar-Human Conflict in British Columbia. PLoS ONE 8(9): e74663. https://doi.org/10.1371/journal.pone.0074663
Editor: Nei Moreira, Federal University of Parana (UFPR)) – Campus Palotina, Brazil
Received: April 25, 2013; Accepted: August 5, 2013; Published: September 11, 2013
Copyright: © 2013 Teichman 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 authors have no support or funding to report.
Competing interests: The authours have declared that no competing interests exist.
Introduction
Conflict between humans and wildlife is a widespread issue of concern in conservation of carnivores [1] [2]. Growing human populations and human encroachment into wildlife habitat, habitat loss and carnivore range expansion increase the likelihood of carnivore-human interactions. Understanding processes associated with temporal and spatial patterns of carnivore-human conflicts is essential to management planning [1].
Carnivore-human coexistence is dependent on maintaining low levels of conflicts [3], with high conflict levels potentially resulting in carnivore population sinks. An attractive sink is an area that animals may perceive as a favourable habitat but where risk of mortality is high [4]. Ideally, data on carnivore occurrence and human-caused mortality can be combined into habitat-based models that enable delineation of attractive sink-like habitats [4] [5]. However, in many cases conservation decisions need to be made in a timely fashion with relatively little information. A natural starting point is to simply describe and analyze spatial and temporal mortality patterns [6]. A more complex and rigurous approach is to develop predictive models for the spatial distribution of human-caused carnivore mortalities [7]. For example, models such as those developed in [8] are useful tools for identifying hotposts of livestock depredation by carnivores.
Loss of livestock to carnivores is a significant conservation challenge wherever livestock ranching practices overlap with carnivore range [9],[10]. In North America, all large carnivore species depredate on cattle (wolves [Canis lupus], [10]; grizzly bears [Ursus arctos], [11]; black bears [Ursus americanus], [12] and cougars [13]). Although radio telemetry-based monitoring of animals can reveal problem individuals such as cougars depredating on livestock [14], unmarked individuals, particularly sub-adult transients, may move into conflict hotspots. It may be these individuals that come into conflict with people [15] and livestock. Because depredation events are more likely to occur repeatedly in areas that have a history of conflict [9], predicting areas with high conflict potential would allow for more focused and preventative measures in conflict hotspots.
Several methods have been applied to reduce carnivore-human conflicts and these can be classified into lethal and non-lethal approaches. Historically, depredation events have been controlled with lethal management to reduce population densities, sometimes causing regional extirpations [16]. In addition, from the 1970’s to the 1990’s, translocation of ‘problem carnivores’ had been a standard management tool [17]. However, the loss or displacement of carnivores have poorly understood effects on population structure [18]. Removal of breeding adults may result in dependent young falling victim to starvation or killed by conspecifics moving into the vacant habitat [19]. Translocation may result in death by resident conspecifics, settlement and reproduction, or killing livestock in neighbouring areas [20]. Such conservation challenges can be minimized if a preventative approach is applied whereby management focuses on areas with high conflict probability [21].
Conflict mitigation is complicated by variability in behaviour between mammalian carnivore genders. For many solitary predators, males occupy larger home ranges [22] and have higher movement rates due to their evolutionary drive to search for females and defend their large territorial boundaries from intruding conspecifics [23],[24]. Also, males often disperse further distances compared to females [25],[26]. Seasonal differences in behaviour are also evident between males and females, for example, when females are raising young and have high energy demands thus spending more time searching and handling prey [27]. These differences may result in greater risk exposure of males in relation to humans such as when crossing roads or moving through ranchlands or wild-urban interfaces. Accounting for gender differences should therefore be considered when attempting to predict areas of high carnivore-human conflict.
Compared to other large carnivores in the Americas, cougars have lost the least of their historic range [28], with their persistence largely attributed to their secretive habits and solitary nature. Throughout their range, cougars come into conflict with people within areas dominated by livestock ranching [14],[29],[30],[31]. Furthermore, human encroachment in cougar habitat can result in cougar attacks on humans by dispersing sub-adult male cougars or emaciated, injured or diseased cougars that have difficulty capturing wild prey [15],[32]. The world’s greatest hotspot of cougar conflict is Vancouver Island located in British Columbia, Canada [15]. However, to date no study has quantified what specific habitat and human factors affect cougar conflicts for this area. It is suggested that British Columbia contains a stable cougar population estimated at 4,000–6,000 individuals, one of the highest of all states and provinces in North America [33]. Nonetheless, reliable estimation of cougar numbers is difficult and assigning ‘sustainable’ mortality rates that ensure cougar population persistence is not straighforward due to limited knowledge of cougar social systems and large-scale movements [18]. Thus, conservation and management decisions involving cougar populations will benefit from understanding cougar-human conflict patterns in order to minimize risk of conflict.
The objective of this study was to describe the temporal trend of cougar-human conflicts in British Columbia over the past thirty years and examine the spatial distribution of conflicts over a recent ten year period in relation to habitat and human features. We predicted that 1) most conflicts will occur in the summer which coincides with the greatest human presence on the landscape, 2) spatial models combining human and habitat features will be most supported and 3) male cougars will be more likely to get into conflict with people in association with human-related landscape features due to their wide ranging patterns and dispersal movements.
Methods
Study area
The study area encompasses cougar range in British Columbia, which coincides with the provincial extent excluding the extreme northern section of the province and most of the Pacific islands. Mean elevation is 1,184±499 m (mean ± sd) and 14.5% of the 765,703 km2 study area is protected under federal or provincial legislation. Mean human density is 4±68 individuals/km2, with the highest densities in the southern portion of the study area, which also has the highest road density. Cattle ranching is distributed mainly in the south-central and north-eastern parts of the province with mean cattle density at 0.3±10.8 cattle/km2.
Temporal patterns of cougar-human conflict
Our data set consisted of thirty years of conflict records collected by the British Columbia Ministry of Environment from 1978–2007, undifferentiated by type of conflict. We therefore define cougar-human 'conflict' as any incident of cougar road mortality, depredation on livestock or attack on humans. Each record included geographically referenced location (UTM coordinates), date, and the type of conflict. We assessed temporal dynamics of conflict occurrence for each sex by season, differentiating between spring (March 21-June 19), summer (June 20-September 21), fall (September 22-December 20), and winter (December 21-March 20). Chi-square goodness of fit tests were used to assess differences in the total number of conflicts between seasons with expected seasonal conflict rates set at 25% of total recorded conflicts. We also tested for differences in the number of conflicts between sexes by performing chi-square goodness of fit tests for each season.
Spatial patterns of cougar-human conflict
We used ten years of conflict data (1998–2007) and GIS spatial layers to estimate which variables adequately predict conflict and map the relative probability of cougar-human conflict. We did not use the entire thirty years of conflict data because accurate spatial information on conflicts was lacking for the 1978–1997 period and due to changes in landscape variables used to predict conflict hotspots. We carried out separate analyses for males and females without partitioning the data seasonally because of low sample sizes (female n = 222, male n = 222). To obtain a population level map of conflict potential, we computed the arithmetic mean between predicted values for the entire landscape, using the male and female predicted probabilites as inputs. Human and habitat characteristics at conflict locations were contrasted with those recorded at random locations to determine factors influencing conflict occurrence. Random locations were generated within British Columbia cougar range using Hawth's Tools [34] and a sampling intensity based on an estimated cougar home range size of 100 km2 [35] resulting in a total of 7,928 locations. Because we were interested in comparing conflict patterns between sexes, and conflicts involving cougars of unknown sex represented only 4.7% of total conflicts, these were excluded from the temporal and spatial analyses.
Predictors of cougar-human conflict.
GIS variables selected to predict conflict locations were used in three sets of models, which included human variables, habitat variables and combined human-habitat variables respectively. Human variables hypothesized to influence conflict occurrence included human population density (http://sedac.columbia.edu/gpw), density and distance to roads (http://geogratis.cgdi.gc.ca/) and cattle density (http://fao.org). Densities and distances to nearest paved and unpaved roads were not included in models because they were correlated with the combined paved-unpaved respective road variables. Habitat variables hypothesized to affect cougar and prey density and thus conflict levels included elevation and terrain ruggedness (http://cits.rncan.gc.ca), land cover (http://bioval.jrc.ec.europa.eu/) reclassified to represent natural habitat occupied by cougars (conifer, deciduous, mixed forest and shrubland), forest edge density, distance to water combined as the minimum of the distances to nearest river and lake, protected area (http://geogratis.cgdi.gc.ca/) and distance to northern edge of the cougar range (http://iucnredlist.org/) (Table 1).
Landscape characteristics in locations of cougar-human conflicts and at random locations were assessed using three moving windows in a GIS representing the average reported sizes of different cougar home ranges in British Columbia [35]. This approach followed Naves et al. [36] who applied it to their study of brown bears in Spain. The analysis resulted in three sets of landscape variables at the scales of 50 km2, 100 km2, and 200 km2 respectively. All layers were in raster format at 1 km2 grid cell size which also was employed in other studies [6],[8].
Modeling approach.
We formulated a priori candidate models using biologically meaningful variables and used identical variable combinations for male and female models respectively. Models were fitted using logistic regression with a 1 (conflict) or 0 (random) outcome using STATA/SE-64. Competing sets of models included human models (n = 7), habitat models (n = 12), and combined human-habitat models. Human-habitat models were further split into simple (without cattle as a predictor, n = 36; with cattle, n = 48) and complex (without cattle but with interaction terms, n = 30; with cattle and interaction terms, n = 40) forms. All variables were screened for correlations (r > |0.7|) and only uncorrelated variables were used within the same model. We computed robust standard errors for all models to account for potential autocorrelation and/or heteroskedasticity in the regression error estimation. Fit of individual models was assessed by carrying out goodness of fit tests and calculating percentage deviance explained.
To rank models we computedAkaike’s Information Criterion corrected for small sample sizes (AICc) for each model and used the difference from the null model (▵AICc) and Akaike weights (w) in all model ranking procedures described below. First, using ▵AICc, three univariate models (one for each scale: 50 km2, 100 km2 or 200 km2) were run for each variable to select the most appropriate scale for the respective covariate. Once the suitable scale was identified separately for each sex, we ran the within-sets candidate models and ranked them for each set. Models that received support in the within-sets ranking procedure (▵AICc < 10) [37] were included in a final (between-sets) AICc ranking protocol that competed human, habitat and combined human-habitat models.
Simpler models are often most appropriate for prediction and easiest to understand, therefore most useful for management decisions. Because information criteria can be sensitive to sample size, to obtain optimal parsimony we re-ranked models included in the final AICc ranking by using a % deviance rule-of-thumb, whereby models needed to improve the percentage deviance explained by at least 1% for each additional parameter included in the model. By this rational, a model with 5 parameters explaining 30.5% of deviance ranked higher than one with 6 parameters explaining 31% of deviance, even though the first model may have had a higher AICc than the second one. This approach resulted in two top models for each sex, one as selected by ▵AICc and a second one as obtained based on deviance explained. We used variance inflation factors to check for collinearity between linear predictor variables in all four top models, eliminating highly collinear variables from the same model when necessary [38]. Finally, we investigated the predictive power of all top models by using 5-fold Cross Validation [39] wherein we partitioned the data into 5 bins and tested the model iteratively on 20% witheld data. Area-adjusted frequencies with significant Spearman-rank correlation coefficients (rho) were considered indicative of good predictive power.
Results
A total of 1,727 conflicts involving cougars of known sex were recorded by the British Columbia Ministry of Environment over the thirty-year study period (female n = 847; male n = 880). The recent ten-year period (1998–2007) included a total of 444 conflicts with associated spatial information and known sex, resulting coincidently in a balanced data set for spatial analyses (female n = 222; male n = 222). Most conflicts occurred in central and southern British Columbia, with Vancouver Island having the largest number of conflicts of all regions within the province, a pattern which held during all three decades of data records (Figure 1).
(A) 1978–1987, (B) 1988–1997, (C) 1998–2007. Darkening shades of red indicate increase in number of conflicts. Numbers on maps represent mean conflict incidence per 10,000 km2 and are presented for each geographical region in the province.
Temporal patterns of cougar-human conflict
Minimum number of conflicts per unit area occurred during the first decade (1978–1987). Conflict numbers increased during the second decade (1988–1997) and finally decreased for the most recent decade (1998–2007), although this is still high relative to the first decade. Frequency of conflict differed between seasons for male (χ2 = 14.94, df = 3, P < 0.05) and female (χ2 = 16.66, df = 3, P < 0.05) cougars, with most conflicts recorded in summer and winter respectively (Figure 2). This pattern of conflict was also confirmed by contrasting frequencies of conflict between sexes on a seasonal basis. Conflict occurrence differed significantly between males and females for summer (χ2 = 5.34, df = 1, P < 0.05) and winter (χ2 = 5.67, df = 1, P < 0.05), but not for spring (χ2 = 0.20, df = 1, P = 0.66) and fall (χ2 = 1.56, df = 1, P = 0.21).
Graphs represent means and error bars represent standard errors.
Spatial patterns of cougar-human conflict
Models supported (▵AICc <10) in the final ranking procedure (between-sets) are reported in Table 2 (male) and Table 3 (female). All models had higher log pseudolikelihoods compared to the male and female null models respectively, including univariate models used to select the most appropriate scale for each variable. This finding confirmed that our a priori selected variables were adequate predictors of cougar-human conflict. For both sexes, models that incorporated human or habitat-only variables received no support (▵AICc ≥ 10) in the final (between-sets) ranking procedure. For both sexes, all human-habitat models competing in the final ranking received support, irrespective of whether interaction terms were present or not. The male and female top models selected via % deviance explained were from the combined human-habitat set, and so was the top ▵AICc - selected female model. The top ▵AICc - ranked male model, however, belonged to the set having human-habitat and interactions. As expected for both sexes, the top ▵AICc - selected model had a larger number of parameters than the top % deviance explained model. Models which included linear and squared terms for cattle density received less overall support than models with no cattle during initial ranking of the human model set, as well as during final ranking (between-sets).
Males.
Based on the top ▵AICc - selected model, male cougar conflicts were more likely to occur on/near roads, far from the northern cougar range in British Columbia, and at intermediate elevations in areas with high proportion of mixed forest (Figure 3a, Table 4). However, based on the interaction term between proportion mixed forest and distance to roads, mixed forest areas far from roads were less likely to have conflict. This model had 13 parameters, an AICc weight of 0.35 and explained 41.2% of the deviance. The model had good predictive power (mean ρ = 0.84, range 0.80−0.93) and its resulting predictions mapped in Figure 4a show conflict risk concentrated along roads and southern British Columbia, as well as high conflict for human populated areas in the central part of the province.
Predicted relative probability is based on variables from (A) Top male ▵AICc model, (B) Top % deviance explained male model, (C) Top ▵AICc female model, and (D) Top % deviance explained female model. Prediction for males are in blue and for females in red. Inset map illustrates conflict predictions for Vancouver Island, with elevation set as transparent in the background. For the inset only the top ▵AICc population-level model predictions are shown due to closely matching predictions with the correponding top % deviance explained model.
(A) Top ▵AICc male model, (B) Top % deviance explained male model, (C) Top ▵AICc female model, (D) Top % deviance explained female model, (E) Top ▵AICc population-level model, (F) Top % deviance explained population-level model. Predictions were based on conflict data for 1998–2007 (female n = 222; male n = 222).
Based on the top % deviance explained - selected model, conflicts were more likely to occur in areas of high human density, on/near roads, far from the northern cougar range edge and at intermediate elevations (Figure 3b, Table 4). The model included 6 parameters, and although it had an AICc weight of 0.01, it explained 39.8% of the deviance. This model also had good predictive power (mean ρ = 0.75, range 0.70−0.87) and its resulting predictions mapped in Figure 4b show conflict hotspots concentrated in southern British Columbia, along road corridors and less so in dense human areas compared to the previous model.
Females.
According to the top ▵AICc - selected model, female cougar conflicts were more likely on/near roads, far from the northern cougar range in British Columbia and at intermediate elevations (Figure 3c, Table 5). This model had 6 parameters, an AICc weight of 0.25 and explained 33.2% of the. The model also included human density but the confidence intervals of the coefficient estimate for this variable overlapped zero. The predictive performance of this model was good (mean ρ = 0.85, range 0.66−0.96) and its predictions in Figure 4c are similar to that of the % deviance model. The one exception is that areas with high human density are not predicted to have high conflict, a difference due to the quadratic form of the human density variable. When compared to the male ▵AICc - selected model, the corresponding female model predicted larger conflict hotspots on/near roads, but fewer conflicts in human populated areas in central British Columbia.
Based on the top % deviance explained-selected model, conflicts were more likely to occur on/near roads, far from the northern cougar range edge and at intermediary elevations (Figure 3d, Table 5). The model included 4 parameters, and although it had an AICc weight < 0.01, it explained 32.5% of the deviance. This model also had good predictive power (mean ρ = ρ0.85, range 0.66−0.93) and its predictions shown in Figure 4d illustrate conflict hotspots along road corridors, with more extensive conflict hotspots on/near roads than in the case of the top % deviance explained male model.
The ▵AICc - selected model had the highest mean Spearman correlation coefficient among the top two male models, which is indicative of slightly better predictive accuracy for the relative probability of occurence of a conflict location. In the case of the top two female models, the ▵AICc and % deviance explained selected models had similar predictive accuracies as illustrated by equal mean rho. Mean relative probabilities of conflict corresponding to population-level conflict risk are mapped in Figure 4e and Figure 4f.
Discussion
Conflicts were lowest in 1978–1987 potentially reflecting low cougar population densities as a result of a North American-wide attempt at reducing top predators [40]. However, cougar density estimates are hard to obtain and relating them to incidence of conflict can be questionable [41]. For example, high conflict incidences have been related to low or declining cougar populations [42],[43]. More recently, human encroachment and habitat fragmentation may have caused the increased conflict incidence observed during 1988–1997 and 1998–2007, although to a lesser extent in the latter decade. Possible causes for the recent downward trend in conflict include increased public awareness as a result of campaigns for mitigating carnivore-human conflict [44] or lowered cougar movement rates within their home ranges in response to increased deer populations, which is the opposite of increased carnivore movements at low prey densities [45].
The seasonal differences in conflict occurrence between male and female cougars likely are due to underlying biological processes characteristic of the two sexes. Conflicts involving males occurred more in the summer, which is the period with the highest human recreational use, although our data did not allow explicit testing of human recreation effects on conflict occurrence. Although females in this species are induced ovulators [46] and reproduction can be seasonal due to naturally fluctuating photoperiods, most births in North America occur during summer months [47],[48],[49]. Females have high energetic requirements during gestation and lactation [50],[51] and kill rates are highest for females with kittens than for any other cougar sex/reproductive class [52]. Accordingly, conflicts for females peaked in winter when cougars are often accompanied by fully dependent kittens. To fulfill their dietary requirements during this time, females may be active for longer periods, have different movement patterns and hunting behaviour, such as using risky areas near roads where ungulate prey may be concentrated for foraging [53]. In addition, females may be scavenging to supplementary feed [54], which would make them more susceptible to conflict if they exploit ungulate road kills.
Spatially, cougar-human conflicts were most likely to occur in southern and central British Columbia, far from the northern edge of the species’ distribution. Although prey mortality risk may increase with harsh snow conditions [55], cougars may not be effective predators at high snow cover [56], which may explain their northerly range boundary in British Columbia. Their distributional range follows closely the distribution of their preferred prey, which in Canada are whitetail deer (Odocoileus virginianus) and mule deer (Odocoileus hemionus) [52],[57]. Conflicts occurred mostly at intermediate elevations, as high mountain areas have low vegetative productivity and thus little availability of cougar prey [58]. With climate warming deer populations are expected to expand northwards [59] and upwards, which may in the future lead to cougar range expansion and increased conflict at higher latitudes and altitudes. Cougars have already been shown to persist at high elevations in the Andes [60].
Our top models for cougar-human conflict included distance to northern range edge and elevation habitat variables which are indicative of broad scale habitat preferences. Finer scale habitat information such as forest edge density or distance to nearest water, both good deer habitat predictors [61],[62], were not present in the best models. From the perspective of predictive value, this is excellent as broader scale environmental data are readily available for many regions. However, one of the top male models included mixed forest and shrub, which are preferred cougar habitat because of deer presence and stalking cover. These two habitat variables do not compromise a model’s predictive value because vegetative land cover data are now widely available. Conifer forest did not predict conflict well likely because needleleaf forest cover is extensive across Brisith Columbia and likely supports lower deer densities, especially relative to browse-rich mixed forest and shrub.
Although cougar distribution in human-modified landscapes is determined by habitat requirements as well as avoidance of people [63], encounter risk may be increased by human expansion into cougar habitat. All top models showed that conflicts were more likely to occur on/near roads irrespective of cougar sex. Roads facilitate human access into cougar areas, and if placed in good cougar habitat may increase vehicular collisions because of high incidence of cougar road crossings. Roads that cross through areas with visibility obstructions may be particularly risky for cougars and other wildlife because they decrease the drivers’ reaction time to a crossing animal. Roads thus represented ubiquous conflict hotspots [64] whereas human density was a good conflict predictor for males but not necessarily for females. Human density was present in both the top male models and in the ▵AICc – selected female model but was absent from the % deviance explained – selected model for females. During dispersal, male cougars encounter human inhabited areas where they are likely to get into conflict. These findings conform with worldwide patterns of widlife-human conflict in relation to roads [65],[66] and male-biased conflict in human-use areas [51],[67].
Cougars depredate on livestock [13],[14] and we expected an increased probability of conflict in British Columbia in areas with high cattle density. In particular, we suspected cattle density would be an important predictor of cougar conflict for males because transients are more likely to depredate livestock [68]. For both sexes conflicts were not predicted to occur in areas of high cattle density, possibly because of easy access to expanding white-tailed deer populations [69].
Carnivore conflict mitigation involves efforts at different scales and we recognize that more accurate conflict predictions could be obtained at the local scale (e.g., the extent of a township) if finer scale variables were utilized. For example, our forest edge density variable used as a surrogate for deer distribution was based on forest versus other land cover types and did not include edge density based on local disturbances such as cutblocks, roads, and seismic lines. Our mapped predictions are thus most informative for analyzing conflict probability at regional or continental levels, typically used in large scale conservation planning strategies. Future Resource Selection Function (RSF) modeling applications of conflict data could benefit from incorporating more refined information, for example records of animal age, and female lactation data gathered during carcass inspection. Reproductive status influences movement behaviour and risk taking in all wildlife species [70]. Such information could improve studies that did not consider demography in mortality predictions.
Conclusions
Carnivore-human conflict is a global conservation challenge that is expected to be on the rise with expanding human populations. In the case of cougars in British Columbia, distance to roads was the most important predictor of conflict locations out of all the potential human influences present in top models. This pattern is eloquently illustrated by the predicted conflict risk for Vancouver Island, the world’s largest cougar conflict hotpsot (Figure 3 inset). Fewer parameters were required to explain conflict patterns for females than for males, reflecting the higher variability of risk exposure by males during their wide ranging excursions. Decreasing current speed limits in predicted conflict hotpsots and introduction of cougar road crossing signs may help to reduce road mortality by creating more time for drivers to react to animals on roads, while simultaneously raising awareness of cougar presence. Because road mortality is common for many large carnivores [71],[72],[73],[74],[75],[76], when possible road access [74] and road density should be limited [77].
The predictive capacity for carnivore-human conflict may be increased by including reproductive status and age structure in conflict modeling. On the other hand, habitat data are increasingly available and we have shown that broad habitat patterns can be adequate predictors of cougar-human conflict. Predictions of potential areas that are likely to have cougar-human conflicts can help target management strategies that reduce cougar mortalities. With reduced occurrences of cougar-human conflicts, people are more likely to accept and thus coexist with cougars in a shared landscape such as southern British Columbia.
Acknowledgments
We thank Tony Hamilton for providing the British Columbia cougar-human conflict data, Charlene Nielsen for helping with GIS layer acquisition and processing and Evelyn Merrill for reviewing an earlier draft of this manuscript.
Author Contributions
Conceived and designed the experiments: SEN KJT. Analyzed the data: KJT BC. Wrote the paper: KJT BC. Performed multiple edits of manuscript: SEN.
References
- 1. Treves A, Wallace RB, Naughton-Treves L, Morales A (2006) Co-managing human-wildlife conflicts: a review. Hum Dimens Wildl 11: 383–396.
- 2. Inskip C, Zimmermann A (2009) Human-felid conflict: a review of patterns and priorities worldwide. Oryx 43: 18–34.
- 3. Loe J, Roskaft E (2004) Large carnivores and human safety: a review. Ambio 33: 283–288.
- 4. Nielsen SE, Herrero S, Boyce MS, Mace RD, Benn B, et al. (2004) Modelling the spatial distribution of human-caused grizzly bear mortalities in the Central Rockies ecosystem of Canada. Biol Conserv 120: 101–113.
- 5. Falcucci A, Ciucci P, Maiorano L, Gentile L, Boitani L (2009) Assessing habitat quality for conservation using an integrated occurrence-mortality model. J Appl Ecol 46: 600–609.
- 6. Rodriguez A, Delibes M (2004) Patterns and causes of non-natural mortality in the Iberian lynx during a 40-year period of range contraction. Biol Conserv 118: 151–161.
- 7. Nielsen SE, Stenhouse GB, Boyce MS (2006) A habitat-based framework for grizzly bear conservation in Alberta. Biol Conserv 130: 217–229.
- 8. Treves A, Naughton-Treves L, Harper EK, Mladenoff DJ, Rose RA, et al. (2004) Predicting human-carnivore conflict: a spatial model derived from 25 years of data on wolf predation on livestock. Conserv Biol 18: 114–125.
- 9. Karlsson J, Johansson O (2010) Predictability of repeated carnivore attacks on livestock favours reactive use of mitigation measures. J Appl Ecol 47: 166–171.
- 10. Morehouse AT, Boyce MS (2011) From venison to beef: seasonal changes in wolf diet composition in a livestock grazing landscape. Front Ecol Environ 9: 440–445.
- 11. Wilson SM, Madel MJ, Mattson DJ, Graham JM, Merrill T (2006) Landscape conditions predisposing grizzly bears to conflicts on private agricultural lands in the western USA. Biol Conserv 130: 47–59.
- 12. Horstman LP, Gunson JR (1982) Black bear Ursus americanus predation on livestock in Alberta Canada. Wildl Soc Bull 10: 34–39.
- 13. Cunningham SC, Gustavson CR, Ballard WB (1999) Diet selection of mountain lions in southeastern Arizona. J Range Manage 52: 202–207.
- 14. Polisar J, Maxit I, Scognamillo D, Farrell L, Sunquist ME, et al. (2003) Jaguars, pumas, their prey base, and cattle ranching: ecological interpretations of a management problem. Biol Conserv 109: 297–310.
- 15. Beier P (1991) Cougar attacks on humans in the United-States and Canada. Wildl Soc Bull 19: 403–412.
- 16. Treves A, Karanth KU (2003) Human-carnivore conflict and perspectives on carnivore management worldwide. Conserv Biol 17: 1491–1499.
- 17. Linnell JDC, Aanes R, Swenson JE, Odden J, Smith ME (1997) Translocation of carnivores as a method for managing problem animals: a review. Biodivers Conserv 6: 1245–1257.
- 18. Treves A (2009) Hunting for large carnivore conservation. J Appl Ecol 46: 1350–1356.
- 19. Czetwertynski SM, Boyce MS, Schmiegelow FK (2007) Effects of hunting on demographic parameters of American black bears. Ursus 18: 1–18.
- 20. Belden RC, Hagedorn BW (1993) Feasibility of translocating panthers into Northern Florida. J Wildl Manage 57: 388–397.
- 21. Treves A, Martin KA, Wydeven AP, Wiedenhoeft JE (2011) Forecasting environmental hazards and the application of risk maps to predator attacks on livestock. Biosci 61: 451–458.
- 22. Lindstedt SL, Miller BJ, Buskirk SW (1986) Home range, time, and body size in mammals. Ecol 67: 413–418.
- 23. Anderson AE, Bowden DC, Kattner DM (1992) The puma on Uncompahgre Plateau, Colorado. Colorado Division of Wildlife Technical Publication 40: 1–116.
- 24. Jedrzejewski W, Schmidt K, Okarma H, Kowalczyk R (2002) Movement pattern and home range use by the Eurasian Lynx in Bialowieza Primeval Forest (Poland). Ann Zool Fenn 39: 29–41.
- 25. Schwartz CC, Franzmann AW (1992) Dispersal and survival of subadult black bears from the Kenai Peninsula, Alaska. J Wildl Manage 56: 426–431.
- 26. Sweanor LL, Logan KA, Hornocker MG (2000) Cougar dispersal patterns, metapopulation dynamics, and conservation. Conserv Biol 14: 798–808.
- 27. Schmidt K, Jedrzejewski W, Okarma H (1997) Spatial organization and social relations in the Eurasian lynx population in Bialowieza Primeval Forest, Poland. Acta Theriol 42: 289–312.
- 28. Laliberte AS, Ripple WJ (2004) Range contractions of North American carnivores and ungulates. Biosci 54: 123–138.
- 29. Iriarte JA, Johnson WE, Franklin WL (1991) Feeding ecology of the Patagonia puma in southernmost Chile. Rev Chil Hist Nat 64: 145–156.
- 30. Mazzolli M, Graipel ME, Dunstone N (2002) Mountain lion depredation in southern Brazil. Biol Conserv 105: 43–51.
- 31.
McKinney T, Wakeling BF, O'Dell JC (2010) Mountain lion depredation harvests in Arizona, 1976–2005. Pages 303–313 in: C. van Riper III, B. F. Wakeling, and T. D. Sisk, editors. The Colorado Plateau IV: Shaping conservation through science and management. The University of Arizona Press, Arizona, USA.
- 32. Mattson D, Logan K, Sweanor L (2011) Factors governing risk of cougar attacks on humans. Hum-Wildl Interact 5: 135–158.
- 33.
Cougar Management Guidelines Working Group (2005) Cougar management guidelines. WildFutures, Bainbridge Island, Washington, USA.
- 34.
Beyer HL. (2007) Hawth's Analysis Tools for ArcGIS™ [online] URL: www.spatialecology.com/htools.
- 35. Spreadbury BR, Musil K, Musil J, Kaisner C, Kovak J (1996) Cougar population characteristics in southeastern British Columbia. J Wildl Manage 60: 962–969.
- 36. Naves J, Wiegand T, Revilla E, Delibes M (2003) Endangered species constrained by natural and human factors: The case of brown bears in northern Spain. Conserv Biol 17: 1276–1289.
- 37.
Burnham KP, Anderson DR (2002) Model selection and multimodel inference. A practical information-theoretic approach. 2nd ed. Springer, New York, USA.
- 38.
Chatterjee S, Hadi AS (2006) Regression analysis by example. John Wiley & Sons, Inc., Hoboke, New Jersey, USA.
- 39. Boyce MS, Vernier PR, Nielsen SE, Schmiegelow FKA (2002) Evaluating resource selection functions. Ecol Modell 157: 281–300.
- 40.
Gill RB (2009) To save a mountain lion: evolving philosophy of nature and cougars. Pages 5–16 in M. Hornocker, and S. Negri, editors. Cougar: ecology and conservation. University of Chicago Press, Chicago, Illinois, USA.
- 41.
Wainwright CJ, Darimont CT, Paquet PC (2010) British Columbia’s neglected carnivore: a conservation assessment and conservation planning guide for cougars. Version 01. Raincoast Conservation Foundation, Sidney, British Columbia, Canada.
- 42.
Lambert CS (2003) Dynamics and viability of a cougar population in the Pacific Northwest. MS Thesis. Washington State University, Pullman, Washington, USA.
- 43. Lambert CMS, Wielgus RB, Robinson HS, Katnik DD, Cruickshank HS, et al. (2006) Cougar population dynamics and viability in the Pacific Northwest. J Wildl Manage 70: 246–254.
- 44.
B.C. Wildlife–human conflict prevention strategy (2003) British Columbia Ministry of Water, Land and Air Protection.
- 45.
Fuller TK, Sievert PR (2001) Carnivore demography and the consequences of changes in prey availability. Pages 163–178 in J. Gittleman, D. Macdonald, S. Funk, and R. Wayne, editors. Carnivore conservation. Cambridge University Press, Cambridge, UK.
- 46. Brown JL (2006) Comparative endocrinology of domestic and nondemestic felids. Theriogenology 66: 25–36.
- 47. Jansen BD, Jenks JA (2012) Birth timing for mountain lions (Puma concolor): Testing the prey availability hypothesis. PLOS ONE 7: e44625.
- 48. Laundre JW, Hernandez L (2007) Do female pumas (Puma concolor) exhibit a birth pulse? J Mammal 88: 1300–1304.
- 49. Ross PI, Jalkotzy MG (1992) Characteristics of a hunted population of cougars in southwestern Alberta. J Wildl Manage. 56: 417–426.
- 50.
Ackerman BB, Lindzey FG, Hemker TP (1986) Predictive energetics model for cougars. Pages 333–352 in S.D. Miller, and D. D. Everett, editors. Cats of the world: biology, conservation, and management. Second International Cat Symposium, Kingsville, Texas, USA, October 4–6, 1982. National Wildlife Federation, Washington, D.C., USA.
- 51. Conde DA, Colchero F, Zarza H, Christensen Jr NL, Sexton JO, et al. (2010) Sex matters: Modeling male and female habitat differences for jaguar conservation. Biol Conserv 143: 1980–1988.
- 52. Knopff KH, Knopff AA, Kortello A, Boyce MS (2010b) Cougar kill rate and prey composition in a multiprey system. J Wildl Manage 74: 1435–1447.
- 53. Rea RV (2003) Modifying roadside vegetation management practices to reduce vehicular collisions with moose Alces alces. Wildl Biol 9: 81–91.
- 54. Knopff KH, Knopff AA, Boyce MS (2010a) Scavenging makes cougars susceptible to snaring at wolf bait stations. J Wildl Manage 74: 644–653.
- 55. Okarma H, Jedrzejewska B, Jedrzejewski W, Krasinski ZA, Milkowski L (1995) The roles of predation, snow cover, acorn crop, and man-related factors on ungulate mortality in Bialowieza Primeval Forest, Poland. Acta Theriol 40: 197–217.
- 56. Alexander SM, Logan TB, Paquet PC (2006) Spatio-temporal co-occurrence of cougars (Felis concolor), wolves (Canis lupus) and their prey during winter: a comparison of two analytical methods. J Biogeogr 33: 2001–2012.
- 57. Bacon MM, Becic GM, Epp MT, Boyce MS (2011) Do GPS clusters really work? Carnivore diet from scat analysis and GPS telemetry methods. Wildl Soc Bull 35: 409–415.
- 58. Pauley GR, Peek JM, Zager P (1993) Predicting white-tailed deer habitat use in Northern Idaho. J Wildl Manage 57: 904–913.
- 59.
Dawe KL (2011) Factors driving range expansion of white-tailed deer, Odocoileus virginianus, in the boreal forest of Northern Alberta, Canada. Ph.D. Thesis, University of Alberta, Edmonton, Alberta, Canada.
- 60. Lucherini M, Reppucci JI, Walker RS, Villalba ML, Wurstten A, et al. (2009) Activity pattern segregation of carnivores in the High Andes. J Mammal 90: 1404–1409.
- 61. Farrell MC, Tappe PA (2007) County-level factors contributing to deer-vehicle collisions in Arkansas. J Wildl Manage 71: 2727–2731.
- 62. Long RA, Muir JD, Rachlow RL, Kie JG (2009) A comparison of two modeling approaches for evaluating wildlife-habitat relationships. J Wildl Manage 73: 294–302.
- 63. Dickson BG, Beier P (2002) Home-range and habitat selection by adult cougars in southern California. J Wildl Manage 66: 1235–1245.
- 64.
Quigley H, Hornocker M (2010) Cougar population dynamics. Pages 59–75 in M. Hornocker and S. Negri, editors. Cougar: ecology and conservation. University of Chicago Press, Chicago, Illinois, USA.
- 65. Trombulak SC, Frissell CA (2000) Review of ecological effects of roads on terrestrial and aquatic communities. Conserv Biol 14: 18–30.
- 66.
Fahrig L, Rytwinski T (2009) Effects of roads on animal abundance: an empirical review and synthesis. Ecol Soc 14: 21 [online] URL: www.ecologyandsociety.org/vol14/iss1/art21/.
- 67. Sukumar R (1991) The management of large mammals in relation to male strategies and conflict with people. Biol Conserv 55: 93–102.
- 68. Linnell JDC, Odden J, Smith ME, Aanes R, Swenson JE (1999) Large carnivores that kill livestock: do "problem individuals" really exist? Wildl Soc Bull 27: 698–705.
- 69. Crete M, Daigle C (1999) Management of indigenous North American deer at the end of the 20th century in relation to large predators and primary production. Acta Vet Hung 47: 1–16.
- 70. Lima SL (1998) Stress and decision making under the risk of predation: Recent developments from behavioral, reproductive, and ecological perspectives. Stress Behav 27: 215–290.
- 71.
Gunther KA, Biel MJ, Robinson HL (1998) Factors Influencing the Frequency of Road-killed Wildlife in Yellowstone National Park. In the Proceedings of the International Conference on Wildlife Ecology and Transportation, Fort Myers, FL 395 – 405.
- 72. Dickson BG, Jenness JS, Beier P (2005) Influence of vegetation, topography, and roads on cougar movement in southern California. J Wildl Manage 69: 264–276.
- 73. McLellan BN (1989) Dynamics of a grizzly bear population during a period of industrial resource extraction. II. Mortality rates and causes of death. Can J Zool 67: 1861–1864.
- 74. Nielsen SE, Herrero S, Boyce MS, Mace RD, Benn B, et al. (2004) Modelling the spatial distribution of human-caused grizzly bear mortalities in the Central Rockies Ecosystem of Canada. Biol Conserv 120: 101–113.
- 75. Taylor SK, Buergelt CD, Roelke-Parker ME, Homer BL, Rotstein DS (2002) Causes of mortality of free-ranging Florida panthers. J Wildl Dis 38: 107–114.
- 76. Colino-Rabanal VJ, Lizana M, Peris SJ (2011) Factors influencing wolf Canis lupus roadkills in Northwest Spain. Eur J Wildl Res 57: 399–409.
- 77. Frair JL, Merrill EH, Beyer HL, Morales JM (2008) Thresholds in landscape connectivity and mortality risks in response to growing road networks. J Appl Ecol 45: 1504–1513.