Abstract Predator control and sport hunting are often used to reduce predator populations and livestock depredations, – but the efficacy of lethal control has rarely been tested. We assessed the effects of wolf mortality on reducing livestock depredations in Idaho, Montana and Wyoming from 1987–2012 using a 25 year time series. The number of livestock depredated, livestock populations, wolf population estimates, number of breeding pairs, and wolves killed were calculated for the wolf-occupied area of each state for each year. The data were then analyzed using a negative binomial generalized linear model to test for the expected negative relationship between the number of livestock depredated in the current year and the number of wolves controlled the previous year. We found that the number of livestock depredated was positively associated with the number of livestock and the number of breeding pairs. However, we also found that the number of livestock depredated the following year was positively, not negatively, associated with the number of wolves killed the previous year. The odds of livestock depredations increased 4% for sheep and 5–6% for cattle with increased wolf control - up until wolf mortality exceeded the mean intrinsic growth rate of wolves at 25%. Possible reasons for the increased livestock depredations at ≤25% mortality may be compensatory increased breeding pairs and numbers of wolves following increased mortality. After mortality exceeded 25%, the total number of breeding pairs, wolves, and livestock depredations declined. However, mortality rates exceeding 25% are unsustainable over the long term. Lethal control of individual depredating wolves may sometimes necessary to stop depredations in the near-term, but we recommend that non-lethal alternatives also be considered.

Citation: Wielgus RB, Peebles KA (2014) Effects of Wolf Mortality on Livestock Depredations. PLoS ONE 9(12): e113505. https://doi.org/10.1371/journal.pone.0113505 Editor: Joseph K. Bump, Michigan Technological University, United States of America Received: July 28, 2014; Accepted: October 24, 2014; Published: December 3, 2014 Copyright: © 2014 Wielgus, Peebles. 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. Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the Supporting Information files. Funding: Funding for this research was provided solely by a research grant from the Washington Department of Fish and Wildlife. 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 Predator control and sport hunting are often used to reduce predator populations, livestock depredations, and increase social acceptance of large carnivores such as brown bears (Ursus arctos) [1], wolves (Canis lupus) [2], cougars (Puma concolor) [3], jaguars (Panthera onca) [4], lions (Panthera leo) [5], leopards (Pantera pardus) [6], and others [7]. Gray wolves (our model animal) are currently being hunted in Idaho, Wyoming and Montana, in part, to reduce livestock depredations [8]. However, to our knowledge, the long-term effectiveness of lethal wolf control to reduce livestock depredations has not yet been rigorously tested. For example, Bradley and Pletscher [9] predicted that breeding pairs are responsible for most livestock depredations because they are bound to the den site, not natural prey distribution [10]. Brainerd et al. [11] predicted that increased wolf mortality could result in fracture of pack structure and increased breeding pairs. If this is the case, increased mortality of wolves could result in increased breeding pairs and livestock depredations following lethal control. In other species, Collins et al. [12] and Treves et al. [13] found increased damages by black bears (Ursus americanus) following high remedial mortality and Peebles et al. [14] found that increased mortality of cougars resulted in increased livestock depredations because of social disruption. In this paper we test the widely accepted, but untested, hypothesis that increased lethal control decreases wolf livestock depredations in a large scale (tri-state) long-term (25 year) quasi-experimental [15]. The “remedial control” hypothesis predicts that livestock depredations will decrease following increased lethal control.

Methods We obtained the confirmed number of cattle (Bos primigenius) and sheep (Ovis aries) depredated, wolf population estimates, number of breeding pairs, and the number of wolves killed in the wolf-occupied area of each state for each year between 1987–2012 from United States Fish and Wildlife Services Interagency Annual Wolf Reports [16] (Table S1). The number of wolves killed were wolves that were killed through control methods including wolves killed legally by livestock owners or through government control methods, these numbers do not include other sources, including natural mortality. Only the total numbers of livestock killed, not the number of confirmed livestock depredation incidents, were available from the USFW database. Numbers of livestock were similarly obtained from United States Department of Agriculture National Agricultural Statistics Service for counties where wolves are present [17] (Table S1). Livestock numbers for individual properties were not available so livestock tallies were made across the tri-state area using counts from wolf occupied counties. Following Peebles et al. [14], we used forward selection, negative binomial general linearized models to assess the relationship between livestock depredations and numbers of livestock, wolves, breeding pairs, and wolves killed. We used this method because depredations were over-dispersed and consisted of 0 to positive integer count data with a variance exceeding the mean [18]. The best statistical model was then selected using the lowest AIC (Akaike Information Criterion) and highest log-likelihood [19]. The rate ratio, analogous to odds-ratio, was computed from the coefficients to aid in interpreting the results [20]. For example, a rate ratio of 1.0 for any independent variable means the effect on the dependent variable is unchanged. A rate ratio of 1.5 means the odds are increased by 50%. To establish directionality, we analyzed the effect of independent variables in year 1 on number of livestock depredated in year 2. This one year time-lag between control kills and depredations removes the directional effect of depredations causing kills. After assessing the models, we plotted and interpreted the most important independent variables against depredations to provide a visual representation of model terms. Only the independent variables that had a rate ratio larger than 1.01or smaller than 0.99, meaning the change in the mean number of livestock depredated was increased or decreased by at least 1%, were examined further. We conducted our regressions on the entire tri-state area and the 3 separate states- but only report the larger tri-state area here because the results were basically the same in all cases (Figures S1 and S2).

Discussion Our results do not support the “remedial control” hypothesis of predator mortality on livestock depredations the following year. However, lethal control of wolves appears to be related to increased depredations in a larger area the following year. Our results are supported by the findings of Harper et al. (2008) in Minnesota where they found that across the state (large scale) none of their correlations supported the hypothesis that killing a high number of wolves reduced the following year's depredations. Harper et al also found that trapping and not catching wolves decreased depredations more than no trapping at all, suggesting that a mere increase in human activity at depredation sites reduced further depredations by those wolves in their study area. By contrast, Bjorge and Gunson (1985) found reducing the population from 40 to 3 wolves in 2 years in Alberta (a 10 fold reduction to near extirpation) resulted in a decline of livestock depredations for two years - followed by subsequent recolonization and increased depredations thereafter. Tompa (1983) also found that lethal control prevented conflict for more than a year in some areas of British Columbia. It should be noted that these 2 studies examined wolf control and livestock depredations at a fine scale (grazing allotment or wolf pack territory or management zone). They did not examine wolf control and livestock depredations at a larger scale (wolf occupied areas) as was done by Harper et al. (2008) and us (this study). It appears that wolf control is associated with reduced depredations at the local wolf pack scale but increased depredations at the larger wolf population scale. This appears consistent with Treves et al. (2005) prediction that the removal of carnivores generally only achieves a temporary reduction in livestock depredations locally when immigrants can rapidly fill the vacancies. There were several different factors that influenced the number of livestock depredated the following year by wolves. In order of importance, based on the values of the rate ratios, these include: the number of wolves removed through control methods, the number of breeding pairs, the minimum wolf population, and the number of livestock on the landscape. Consistent with expectations, each additional breeding pair on the landscape increased the expected mean number of cattle depredated by 8 to 9% and each additional wolf on the landscape increased the expected mean number of sheep depredated by 6%. Cattle were most affected by breeding pairs and sheep by wolves – perhaps because it takes more than one wolf (a pack) to kill a relatively larger cow and only one wolf to kill a smaller sheep. However, contrary to the “remedial control” hypothesis, each additional wolf killed increased the expected mean number of livestock depredated by 5–6% for cattle and 4% for sheep. It appears that lethal wolf control to reduce the number of livestock depredated is associated with increased, not decreased, depredations the following year, on a large scale – at least until wolf mortality exceeds 25%. Why 25%? The observed mean intrinsic growth rate of wolves in Idaho, Wyoming, and Montana is about 25% [21]. Therefore, once anthropogenic mortality exceeds 25%, the numbers of breeding pairs and wolves must decline – resulting in fewer livestock depredations. Below 25% mortality, lethal control may increase breeding pairs and wolves through social disruption and compensatory, density dependent effects. For example, wolf control efforts occur year round and often peak during grazing season in areas with livestock depredations [22], [23]. However, if control takes place during the breeding season and a member of the breeding pair is removed it may lead to pack instability and increased breeding pairs [24], [10]. Furthermore, loss of a breeder in a pack during or near breeding season can result in dissolution of territorial social groups, smaller pack sizes and compensatory density dependent effects – such as increased per-capita reproduction [11], [25], [26]. Culling of wolves may also cause frequent breeder turnover [11] and related social disruption – which can result in reduced effective prey use (through loss of knowledge of prey sources and ability to subdue prey) which may also result in increased livestock depredations [27], [28]. All of these effects could potentially result in increased livestock depredations. We would expect to see increased depredations, wolves killed, and breeding pairs as the wolf population grows and recolonizes the area - but our data suggest that lethal control exacerbates these increases. The secondary effects of time, wolf population growth rate, wolf occupied area, and wolf population size on depredations were already subsumed in the primary main effect terms of breeding pairs (cattle) and wolves (sheep), so those secondary effects cannot account for the positive effects of wolf kills on depredations. We do not yet know the exact mechanism of how increased wolf mortality up to ≤25% results in increased livestock depredations, but we do know that increased mortality is associated with compensatory increased breeding pairs, compensatory numbers of wolves, and depredations [24], [10], [27], [28], [11], [26]. Further research is needed to determine the exact causal mechanism(s). Annual mortality in excess of 25% will reduce future depredations, but that mortality rate is unsustainable and cannot be carried out indefinitely if federal relisting of wolves is to be avoided. Furthermore, a 5% (sheep) and 5% (cattle) kill rate of wolves yields the same number of cattle and sheep depredations as a 35% (cattle) and 30% (sheep) kill rate (Figures 3 & 6), but the 30% or 35% rate is unsustainable for wolf population persistence and the 5% rate is not. The worst possible case appears to be a high mortality rate at about 20–25%, since this corresponds to a “standing wave” of the highest livestock depredations. Further research is needed to test if this high level of anthropogenic wolf mortality (25%) is associated with high levels of predation on natural prey such as deer and elk. Further research is also needed to account for the limitations of our data set. The scale of our analysis was large (wolf occupied areas in each state in each year) and the scale of some other studies were small (wolf packs). Simultaneous, multi-scale analysis (individual wolf packs, wolf management zones, and wolf occupied areas) may yield further insights. Although lethal control is sometimes a necessary management tool in the near-term, we suggest that managers also consider testing non-lethal methods of wolf control [29] because these methods might not be associated with increased depredations in the long-term.

Acknowledgments This analysis and paper benefitted from the insights and comments of Hilary Cooley (U.S. Fish and Wildlife Service), and John Pierce, Donny Martorello, Brian Kertsen, Ben Maletzke, and Stephanie Simick (Washington Department of Fish and Wildlife).

Author Contributions Conceived and designed the experiments: RBW KAP. Performed the experiments: RBW KAP. Analyzed the data: RBW KAP. Contributed reagents/materials/analysis tools: RBW KAP. Wrote the paper: RBW KAP.