Intra-seasonal variation in wolf Canis lupus kill rates.
In our study, we observed how kill rates for wolves killing white-tailed deer Odocoileus virginianus varied from week to week throughout a four-month period (December-March) for each of five packs observed in each of four winters. We estimated kill rates by following wolf tracks in the snow. This is distinctive because most researchers used aerial observations to estimate kill rate for wolves. For example, in a review of kill rate estimates for wolves, researchers used ground tracking for only two of 124 estimates (Schmidt & Mech 1997; the ground-based estimates were from Stenlund (1955) and Kolensky (1972)).
Our ground-tracking experiences and the tremendously variable kill rates which we observed (see section Results) inspired us to hypothesize that many published estimates of wolf kill rates were affected by two important, but unappreciated, methodological factors: 1) duration of observation and 2) whether estimates were based on ground or aerial observations. Our observations suggested that wolves only sustained high kill rates for a short time. For this reason, we thought that the duration of study could affect published estimates of kill rate. Here, we aimed to assess the extent of bias in higher estimates of kill rates, especially those associated with shorter study durations.
We also observed many instances in which wolves either killed multiple deer in one predation event or bedded far from their kill site. These patterns provided a basis for better estimating the extent to which aerial observations lead to underestimates of kill rates for deer-eating wolves. We used published data to assess these hypotheses.
Our study area encompassed 1,940 [km.sup.2] in the western portion of Michigan's Upper Peninsula, USA (Fig. 1). Wolves recolonized the study area in the early 1990s (Beyer et al. 2009). Aerial telemetry and ground-based track surveys conducted each year of this study indicated the presence of 25-35 wolves in six to nine packs in the study area. The track survey methods we used are described in Potvin et al. (2005), and we used aerial telemetry methods similar to those described in Fuller (1989).
Wolf density varies nearly five-fold across the upper Michigan counties, with the greatest densities (26 wolves/1,000 [km.sup.2]) occurring in Upper Michigan's two southwestern-most counties (7,900 [km.sup.2]), and the lowest densities (5 wolves/1,000 [km.sup.2]) occurring in the seven eastern-most counties (26,400 [km.sup.2]). Compared to these values, wolf densities in the study area (13-18 wolves/1,000 [km.sup.2]) were intermediate (D. Beyer, unpubl. data).
The study area represents high-quality wolf habitat insomuch as deer densities were high (> 12 deer/ [km.sup.2]; Hill 1999), and human density (< 7 people/[km.sup.2]) and road density (< 0.45 km/[km.sup.2]) were low (Potvin et al. 2005). Habitat analysis and recent patterns of population growth suggest that the wolf population in this area had not yet reached equilibrium with the deer population during our study (Potvin et al. 2005).
Northern hardwoods and transition boreal forest cover almost all of the pack territories that we studied. Elevations range from 184 to 604 m. Small lakes and rivers are numerous. Most of the study area is located within the Ottawa National Forest. Commercial timber harvesting (selective and clear-cut) and recreational deer hunting are important extractive activities in the area. Mean monthly temperatures are -12.2[degrees]C in January and 18.6[degrees]C in July. Mean annual snowfall in the study area varies between 300 and 450 cm (National Weather Service data collected between 1951 and 1989).
We estimated daily per capita and per pack kill rates for five different packs (2-4 packs/year) on as many days as possible during the winters (1 December-1 May) of 2000/01, 2001/02, 2002/03 and 2003/04 in the study area. To do this, we followed wolf tracks in the snow by foot and snow machine. We used radio telemetry to aid in finding wolves, to determine when tracks or kills were made and to keep us from disturbing pack movements (i.e. to maintain a minimum distance of approximately one kilometer from the wolves). From the tracks, we determined the number of wolves in the hunting group each day.
Because packs occasionally split into smaller hunting groups for several days at a time, the daily hunting group was not always the same as the pack size (group size using one territory). From tracks, we also estimated the number of deer killed each day. Kills from one day were distinguished from kills of the next or previous day by telemetry and track-based knowledge of the pack's movements and location on each day, and by condition of the kill (e.g. freshness of kill, how frozen the carcass was and the amount of snow covering the kill). We only included days when we were confident that we accounted for all the activities of the pack.
Our analysis was based on 359 reliable estimates of daily kill rate throughout the four winters (see Table 1 and Fig. 2). Reliable estimates were not obtainable for every day during this interval because at times fresh snowfall, deer tracks or snowmobile traffic would cover up the wolves' tracks before we could follow them. Our assessment of temporal trends in kill rates are reliable if days missed did not tend to be days with greater than average or less than average kill rates. We are not aware of any possibility of systematic bias in our data based on days with missing data.
The reliability of the 359 daily estimates also depends on the certainty of a deer's cause of death. We based the 359 estimates of daily kill rate on necropsies of 206 deer. Of these deer, we judged wolf predation to be the cause of death for 82% of the specimens, based on tracks and other signs that wolves chased and killed the deer. In the remaining 18% of cases, the cause of death was not certain but evidence strongly suggested that wolves were responsible. We refer to these as probable wolf kills. Probable wolf kills were typically associated with tracks of multiple predators including bobcats Lynx rufus or coyotes Canis latrans. In five additional cases of possible wolf predation, there was substantial uncertainty in the cause of death, and these were not included in this analysis. This occurred, for example, when we observed wolf tracks at a deer carcass, but it was unclear whether the wolves had killed it or had merely scavenged the deer after death from some other cause (i.e. starvation, malnutrition, disease, human/hunter, car collision, bobcat or coyote). Because we were conservative in our classification, the best estimate of kills would include definite kills and most of the probable wolf kills. Nevertheless, we performed all subsequent analysis on two data sets, one comprised of only definite kills and a second comprised of both definite and probable kills. Because the results for both analyses were qualitatively identical, we only report the results of the analysis for the combined data set of definite and probable kills. We estimated consumption rates (kg/ wolf/day) by assuming that the consumable mass from each deer was 54.5 kg (Schmidt & Mech 1997).
We analyzed these data to assess: 1) overall variability of intra-seasonal kill rates, 2) influence of snow depth on kill rates, and 3) the extent to which kill rates tended to increase throughout winter. We conducted analyses using Sigmastat and SPSS. Poisson regression is appropriate for these data because the response variable is count data (number of kills/pack/day) with a range of 0-9 (see legend for Fig. 4). We assessed the 'goodness-of-fit' of Poisson models, by the residual deviance divided by the degrees of freedom. If this quotient is less than approximately four, a model's goodness-of-fit may bejudged reasonable (Eberhardt 1978; see also Burnham & Anderson 1998:53). We also used P-values and [AIC.sub.c] statistics (Anderson et al. 2000) to judge the appropriateness of each model. For models assessing the extent to which snow depth accounts for intra-seasonal variation in kill rate, we estimated the proportion of observed variation that a Poisson regression can potentially explain (Fridstrom et al. 1995):
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].
Wolf radio collaring
The wolves in our study had been fitted with very high frequency (VHF) radio-collars (Telonics Inc., Mesa, Arizona, USA) as part of a larger study conducted by the Michigan Department of Natural Resources, whose aim had been to assess population abundance and survival rate in Upper Michigan (Beyer et al. 2009). We captured wolves to attach radio-collars in spring and summer using methods similar to Mech (1974) and Kuehn et al. (1986). We live-captured wolves with foot-hold traps modified to reduce injury (Minnesota Trapline 760, also New-house Modified 14 and McBride No. 7), and we chemically anesthetized (using ketamine hydrochloride and xylazine, both at 100 mg/ml) wolves at doses of 0.11 mg/kg and 2 mg/kg, respectively. We also administered penicillin, and until 2003, we administered vaccinations against sarcoptic mange, canine distemper and canine parvovirus. Each pack that we studied included at least one radio-collared wolf for each year of observation.
Duration of observation
To assess the influence of study duration on estimates of kill rate (i.e. kg/wolf/day), we used 71 published estimates of kill rate compiled from 11 studies by Schmidt & Mech (1997). Their analyses indicate that per capita kill rate tends to decline with pack size and is substantially lower for wolves that feed primarily on white-tailed deer than for wolves that feed primarily on moose Alces alces. Although Schmidt & Mech (1997) collected information on the duration over which each kill rate was estimated, they did not assess the relationship between kill rate and duration of study. In their data set, the duration of study varied from six to 179 days (inner quartile range = 19, 27.5; N = 71) for wolves that prey on moose. Because the duration of study varied little for deer-eating wolves (inner quartile range = 120, 120; N = 31), we restricted our analysis to moose-eating wolves. We quantified the relationship between study duration and kill rates using Poisson regression models.
Of the 701 carcasses that we observed wolves to have fed on, 91% were white-tailed deer, 4% were snowshoe hare Lepus americanus, 3% were grouse Bonasa umbellus and Falcipennis canadensis and 2% were beaver Castor canadensis. Of the 638 deer carcasses on which wolves had fed, wolves had killed 77%, 10% were scavenged from bobcat or coyote kills, 4% had starved, 7% were killed by hunters and 2% had been killed by vehicles.
Inter-annual variation in winter kill rates
We calculated 12 estimates of winter kill rates for five different packs in four different winters (Table 1). The average kill rate, among these packs and among these winters, was 0.68 kills/pack/day or 7.7 kg/wolf/ day. Accounting for differences in pack size, the kill rate for an average wolf, during an average year, was 6.7 kg/wolf/day (i.e. this value is a weighted average of data in Table 1, where the weights varied according to pack size). The coefficient of variation (CV) and interquartile range among season-long kill rates for different years was 55% and 0.33, 0.97 (see Table 1).
The maximum kill rate (based on 38 days of observation) was 2.2 times the average, and the minimum kill rate (based on 27 days of observation) was only 30% of the average. Variation in pack size accounts for 29% of the variation in kill rate (P = 0.07, linear regression). Of the remaining variation, 43% is attributable to differences between years (P = 0.09), 29% to unidentified differences between packs (perhaps attributable to differences in habitat quality or hunting skill; P = 0.01), and 27% is attributable to unidentified sources that include sampling error and random fortune (ANOVA on residuals of the pack-size/ kill-rate regression).
Intra-annual variation in winter kill rates
Daily kill rate per pack was substantially lower in early winter than in late winter (Fig. 2). On average (i.e. all packs and all years combined), kill rates in late winter (mid-February-April) were twice the kill rates observed during early winter (December-mid-February; 0.88 vs 0.44 kills/pack/day, P = 1.5 x [10.sup.-3], [N.sub.late] = 182, [N.sub.early] = 105; t-test). A moving average of the temporal pattern in daily kill rate reinforced this interpretation (see Fig. 2). As winter progressed, a slight decrease occurred in the proportion of days with zero kills, but a substantial increase occurred in the proportion of days with multiple kills per day (Fig. 3). Importantly, as average kill rate increased so did the variation in kill rate (Fig. 4A). For each of the years between 2001 and 2004, CV in weekly kill rates were 90, 830, 48 and 890% (see Fig. 2). The interquartile ranges for kills/pack/day were 0.3, 1.0, 0.4, 1.3, 0.3, 0.7 and 0.4, 1.2 for each of these years, respectively.
To assess the extent to which kill rates tend to increase throughout the winter, we fit two sets of Poisson regression models to our data. One set of regression models included Julian day (JD and [JD.sup.2]) and pack size (PS and [PS.sup.2]) as candidate predictor variables. The most parsimonious of these models is shown in Figure 4A and was:
kills/pack/day = exp(-1.8 + 1.4 x [10.sup.-2]JD + 0.15PS) (1).
Kill rate increases significantly with JD (P = 1.5 x [10.sup.-6]) and PS (P = 1.9 x [10.sup.-6]). The goodness-of-fit for this model was reasonable (1.30). Although we constructed a model that included a quadratic term for Julian day (which would allow for the detection of declining kill rate at the end of winter), equation (1) slightly outperformed the quadratic model ([AIC.sub.c] = 0.55 for the model with a quadratic term for JD, and zero for equation (1)). However, examination of the moving averages suggested that kill rate commonly declined toward the end of March (see Fig. 2).
The regression analysis represented by equation (1) ignores inter-annual variation. Thus, we conducted similar regression analyses to assess the extent to which intra-annual variation differs among years (see Fig. 4B). During a mild winter (2003), the expected daily kill rate in late winter (March) was 30% greater than it was in early winter (early January). During a severe winter (2002), the expected daily kill rate in late winter (March) was 7.5 times greater than it was in early winter (December-January). Even in the mild winter, intra-annual variation in kill rate seems biologically significant from the perspective of both wolf nutrition and deer demography.
To assess the extent to which snow depth accounts for intra-seasonal variation in kill rate, we analyzed kill rate in response to PS and daily snow depth (SD). We estimated SD by taking an average of multiple snow depths recorded along the travel routes made by wolves. Specifically, for a set of Poisson models with SD, [SD.sup.2], PS and [PS.sup.2] as candidate predictor variables, the most parsimonious of these models was (see Fig. 4C):
kills/pack/day = exp(-3.5 + 6.8 x [10.sup.-2]SD - 3.9 x [10.sup.4][SD.sup.2] + 0.16PS) (2).
Kill rate was significantly related to snow depth (P = 3.0 x [10.sup.-3] for SD, and P = 3.9 x [10.sup.-2] for [SD.sup.2]) and pack size (P = 4.0 x [10.sup.-6]). The goodness-of-fit for this model was reasonable (i.e. 1.2). Equation (2) is, according to [AIC.sub.c] weights, 3.5 times more likely than a model that includes only linear terms for SD and PS. Equation (2) also outperformed equation (1) (i.e. [AIC.sub.c] = 16.2 for equation (1) compared to equation (2)).
To better understand the effect of pack size and snow depth on kill rate, we calculated the proportion of observed variation ([R.sup.2.sub.p]) potentially explained by the three different models (i.e. equation (1), a model with only PS and a model with only SD and [SD.sup.2]). The values of [R.sup.2.sub.p] were 0.08 for the pack size only model, 0.18 for the snow depth model and 0.30 for equation (1).
Duration of observation
We assessed the influence of study duration on kill rate by building the following Poisson regression model with PS and study duration (ln(D)) as predictor variables (Fig. 5):
kg/wolf/day = exp(3.85 - 0.065PS - 0.315 ln(D)) (3).
Per capita kill rate declined significantly with pack size (P < [10.sup.-4]) and duration of study (P < [10.sup.-4]). To be sure that the influence of pack size is isolated on the right side of the equation, we also estimated a regression model with kg/pack/day as the response variable and with PS and ln(D) as predictor variables, and the best fitting model for observed kg/ pack/day was:
kg/pack/day = exp(4.41 + 0.094PS - 0.301 ln(D)) (4).
Per pack kill rate again declined with duration of study (P < 7.0 x [10.sup.-4]) and increased with pack size (P < [10.sup.-4]). Although models based on D, rather than ln(D), were also statistically significant, they exhibited a non-linearity (when plotted on a log-log scale) that was not apparent in the observations (see Fig. 5).
The goodness-of-fit for equations (3) and (4) were 4.3 and 27.5, respectively, indicating that equation (3) had a better fit. These goodness-of-fits are likely low because equations (3) and (4) do not account for many factors known to affect kill rate, such as prey density and winter severity (Post et al. 1999, Vucetich et al. 2002). However, it is unlikely that failure to account for these factors led to the spurious appearance of a relationship with D.
Nevertheless, to further increase confidence in a belief that study duration affects kill rate, we estimated models represented by equations (3) and (4) using median regression, which makes no assumptions concerning the distribution of the response variable (Koenker & Bassett 1978). PS and ln(D) were statistically significant (Ps < [10.sup.-4]) for the median regression models.
To highlight the biological significance of these relationships, consider that estimates of kill rate based on 30 days of monitoring tend to be 50% greater than estimates based on 120 days (Table 2).
Ground-based vs aerial estimates of kill rate
Mean kill rate for aerial-based observations of deer-eating wolves was 2.9 [+ or -] 0.8 (SE) kg/wolf/day and 7.7 [+ or -] 1.5 kg/wolf/day for ground-based estimates (Fig. 6). Regression analyses (see Fig. 6) indicate that ground-based estimates of kill rate (data from Table 1), for any given pack size, tend to be greater than aerial-based estimates (from Schmidt & Mech 1997). More specifically, for a typical pack size of four, ground-based estimates were 3.3 times greater than aerial-based estimates (2.4 vs 7.9 kg/wolf/day; see Fig. 6).
We observed two kinds of wolf behaviour that may account for this previously underappreciated discrepancy. First, more often than previously appreciated, we observed wolves killing more than one deer in a single location during a single night (see Fig. 3). Multiple kills occurred on 16% of the observed days, but represented 54% of the total number of observed kills. The frequency of multiple kills is likely related to winter severity (see Fig. 3). Second, for 51% of the 274 observed predation events, wolves bedded > 500 meters from the kill site during the day that followed a kill (Table 3). When wolves bedded > 500 meters from a kill site the mean distance was 1.9 km (the interquartile range was (0.5 km, 2.0 km), and the range of distances was (0.5 km, 12.0 km)). When we could make precise estimates of distance between kill and bed sites (N = 67), > 70% of the bed sites were > 0.5 km from the kill site (Fig. 7).
Kill rate variability
Previous work has indicated that major sources of variability in kill rate are attributable to pack size, prey density, winter severity and unidentified sources (Schmidt & Mech 1997, Post et al. 1999, Vucetich et al. 2002). Our results are consistent with these and add additional insights. Previous research has only shown us that wolf kill rates for deer tend to increase with deeper snow (Pimlott et al. 1969, Mech & Frenzel 1971, Kolensky 1972, Nelson & Mech 1986, Huggard 1993, DelGiudice 1998). The synchrony between kill rate and snow depth (see Fig. 2), together with our field observations of conditions when packs killed > 1 deer/day, suggest how kill rates are associated with specific winter events, such as storms and crust formation, the timing of which are unpredictable. Snow depth explained about twice the intra-seasonal variation in kill rate than pack size did (see [R.sup.2.sub.p] statistics in the section Results). In addition, predicted kill rates seem to asymptote at snow depths of about 80 cm (see Fig. 4C), suggesting this may be the snow depth where snow conditions are so difficult for deer that wolves become satiated or kill rates become constrained by the time it takes to capture and consume prey.
In addition to snow depth, other difficult-to-measure aspects of snow also likely have an important influence on wolf-deer interactions (e.g. snow density, day-night freeze-thaw phenomena and variability of snow crusts in relation to coniferous tree cover). For example, in our study area, snow becomes crusted (increased hardness) from mid-winter thaws and rain. These conditions favour wolves because they have lower foot loadings than their ungulate prey (Mech & Peterson 2003). Important details may be unpredictable; i.e. the three highest daily kill rates that we measured were associated with heavy snowfall accompanied by high winds.
Our analyses also add insights about temporal variation of kill rates. The CV and interquartile ranges for weekly kill rates were much greater than for annual kill rates suggesting that the variability in rates of food acquisition from week to week far exceeds the variability observed from year to year.
Among wolves in Upper Michigan, important sources of variation included pack size, year effects (which includes temporal variability in climate and prey density), pack effects (which includes variation in pack killing ability, spatial variation in prey density and habitat quality), sampling error and random fortune. The fundamental importance of diverse sources of variation in determining kill rate is not always appreciated (e.g. Messier 1994).
Our results indicate that late winter (March) kill rates were between 40% and 12 times greater than early winter (December) kill rates (see Fig. 4B). These patterns are comparable to those for elk-killing wolves in Yellowstone National Park, where later winter kill rates tend to be between 15 and 61% (interquartile range for N = 12 years) greater than early winter kill rates (Smith et al. 2004). In Yellowstone, early winter kill rates are only roughly correlated with late winter kill rates (r = 0.63, P = 0.03, N = 10; Smith et al. 2004). Though only four years of data (i.e. N=4) are available, it also seems as though early and late winter kill rates are not strongly correlated for Michigan's deer-killing wolves (r = -0.38, P = 0.62).
For wolves, and most predators, annual kill rates are often extrapolated from measurements taken from a period of a few weeks (e.g. Messier 1994, Varley & Boyce 2006). The extreme variability of kill rates within a winter suggests that such an extrapolation is unreliable. Rather than short-term estimates of kill rate, the best way to understand the annual impact of wolf predation is likely the assessment of how predation rate affects total annual mortality rate of prey (Marshal & Boutin 1999, Vucetich et al. 2011).
For a typical pack size of four wolves, ground-based estimates were 3.3 times greater than aerial-based estimates, even when the latter were inflated by 150% (see Fig. 6). By comparison, Fuller (1989) presented indirect evidence suggesting that actual kill rates for deer-killing wolves could be two times the estimates based on aerial observations.
The reasons for these differences are likely to include: 1) the difficulty of seeing, from an aircraft, deer carcasses, which are smaller and consumed more quickly, 2) the tendency for wolves to kill multiple deer in one predation event, but not multiple moose, and 3) wolves not bedding near their kill sites.
Multiple kills seem associated with: 1) wolves feeding extensively on one kill and returning (sometimes days) later to consume the other carcasses, and 2) a first kill being made early in the evening (i.e. within 1-2 km of the previous day's resting site). Ground-tracked wolves typically traveled 10-20 km/ night.
More often than not, wolves do not bed near their kill sites during the day (see Table 3 and Fig. 7). In cases where wolves bedded far from their kill, it appeared that they did so for a variety of reasons: to bed at an area farther from human activity or higher up on a hilltop or ridge (46%), to be near a past den site (20%), they had completely consumed the carcass (14%), and multiple kills had been made and the pack could only be near one kill site at a time (19%). On one occasion, a pack that was chasing a trespassing wolf interrupted the chase to kill a deer, and then resumed pursuit after feeding briefly on their kill.
These behaviours, in addition to other difficulties of detecting kills from aircraft, may explain much of the discrepancy between ground-based and aerial-based estimates of kill rate. In mid-winter (February), counting multiple kills as single kills will result in an estimated kill rate that is only 70% the actual kill rate (0.66 and 0.46 deer/pack/day, N = 146; Table 4). In later winter (March-April), the same error results in an estimate that is only 55% of the actual kill rate (0.95 and 0.52 deer/pack/day, N = 159). Other kills detected from the ground were associated with very few remains; these kills would also likely have remained undetected from the air. It is also plausible that high deer density or deer:wolf ratios explain, in part, the higher kill rates that we observed (Vucetich et al. 2002). These observations may explain, at least in part, why deer-eating wolves have appeared to acquire food at slower rates (kg/wolf/ day) than moose-eating wolves (Schmidt & Mech 1997).
Although the kill rates that we observed are the highest average kill rates observed in a population of deer-eating wolves (see Fig. 6), they likely reflect that other reported values had been underestimates. Moreover, the kill rates that we observed were comparable to those observed among moose-eating wolves (compare Table 1 and Fig. 5). The tendency to underestimate kill rates for deer-eating wolves is associated with this observation: Although physiological considerations indicate that free-ranging wolves require approximately 3.0-3.25 kg/wolf/day (Peterson & Ciucci 2003), wolves lose 20-50%, depending largely on pack size, of their captured prey to scavengers (Fig. 4.10 of Peterson & Ciucci 2003, Vucetich et al. 2004). Even in our study, we observed two extreme cases of kleptoparasitism. In one case, dozens of ravens Corvus spp. and eagles Haliaeetus alone removed ~20 kg of meat in 4.5 hours from a wolf-killed deer. In a second case, ravens and eagles consumed or cached ~80 kg from a large adult doe in two days. According to aerially-derived estimates, wolves acquire on average between 1 and 3.4 kg/wolf/day (see Fig. 6). These observations are inconsistent because after scavenging losses are taken into account, aerial-based estimates of kill rate suggest that wolves acquire, on average, less prey than is required for survival. In some cases, the apparent discrepancy might be attributable to many wolves being chronically malnourished and often dying from malnutrition. However, biased estimates of kill rate may represent a better explanation of these discrepancies.
Our analyses also show that kill rate estimates are sensitive to duration of observation. Although most estimates of per capita kill rate are based on < 30 days of observation (see Fig. 5), season-long estimates of winter kill rate tend to be only 2/3 as great as those based on 30 days of observation (see Table 2).
Because wolf-prey ecology depends critically on interpreting and comparing kill rate estimates, the lack of appreciation of these important methodological issues is a cause for concern (see also Hebblewhite et al. 2003). Moreover, while most estimates are based on aerial observations, recent endeavours based on ground-tracking has produced other insight that would not have been possible from a light aircraft (Jedrzejewski et al. 2002, Hebblewhite et al. 2003).
Acknowledgements--we thank D. Lonsway, B. Johnson, M. Haen, K. Sitar, D. Wagner, T. Gouza, J. Lukowski, B. Roell, R. Atkinson, M. Joseph and T. Weise for their contributions to field work, and T. Cooley for performing wolf necropsies. We thank R. Loo, D. Minet and N. Harri for their piloting skills. We also thank P. Lederle for administrative support. Our work was supported in part by Federal Aid in Wildlife Restoration, Michigan Project W-147-R.
Anderson, D.R., Burnham, K.P. & Thompson, W.L. 2000: Null hypothesis testing: Problems, prevalence, and an alternative. --Journal of Wildlife Management 64: 912-923.
Beyer, D.E., Peterson, R.O., Vucetich, J.A. & Hammill, J.H. 2009: Wolf Population Changes in Michigan.--In: Wydeven, A.P., VanDeelen, T.R. & Heske, E.J. (Eds.); Recovery of Gray wolves in the Great Lakes Region of the United States: An Endangered Species Success Story. Springer Press, New York, New York, USA, pp. 65-86.
Burnham, K.P. & Anderson, D.R. 1998: Model Selection and Inference.--Springer-Verlag, New York, New York, USA, 353 pp.
DelGiudice, G.D. 1998: Surplus killing of white-tailed deer by wolves in northcentral Minnesota.--Journal of Mammalogy 79: 227-235.
DelGiudice, G.D., Mech, L.D. & Seal, U.S. 1990: Effects of winter undernutrition on body composition and physiological profiles of white-tailed deer.--Journal of Wildlife Management 54: 539-550.
Eberhardt, L.L. 1978: Appraising variability in population studies.--Journal of Wildlife Management 42: 207-238.
Fridstrom, L., Ifver, J., Ingebrigsten, S., Kulmala, R. & Thomsen, L.K. 1995: Measuring the contribution of randomness, exposure, weather, and daylight to the variation in the road accident counts.--Accident Analysis and Prevention 27: 1-20.
Fritts, S.H. & Mech, L.D. 1981: Dynamics, movements and feeding ecology of a newly protected wolf population in northwestern Minnesota.--Wildlife Monographs 80, pp. 3-79.
Fuller, T.K. 1989: Population dynamics of wolves in north-central Minnesota.--Wildlife Monographs 105, pp. 3-41.
Hebblewhite, M., Paquet, P.C., Pletscher, D.H., Lessard, R.B. & Callighan, C.J. 2003: Development and application of a ratio estimator to estimate wolf kill rates and variance in a multiple-prey system.--Wildlife Society Bulletin 31: 933-946.
Hill, H.R. 1999: The 1999 deer pellet group surveys.--Michigan Department of Natural Resources, Lansing, Michigan, USA, 15 pp.
Huggard, D.J. 1993: Effect of snow depth on predation and scavenging by gray wolves.--Journal of Wildlife Management 57: 382-388.
Jedrzejewski, W., Schmidt, K., Theuerkauf, J., Jedrzejewska, B., Selva, N., Zub, K. & Szymura, L. 2002: Kill rates and predation by wolves on ungulate populations in BialowieZa Primeval Forest (Poland).--Ecology 83: 13411356.
Koenker, R. & Bassett, G. 1978: Regression quantiles.--Econometrica 46: 33-50.
Kolensky, G.B. 1972: Wolf predation on wintering deer in east-central Ontario.--Journal of Wildlife Management 36: 357-369.
Kuehn, D.W., Fuller, T.K., Mech, L.D., Paul, W.J., Fritts, S.H. & Berg, W.E. 1986: Trap-related injuries to gray wolves in Minnesota.--Journal of Wildlife Management 50: 90-91.
Marshal, J.P. & Boutin, S. 1999: Power analysis of wolf-moose functional responses.--Journal of Wildlife Management 63: 396-402.
Mech, L.D. 1974: Current techniques in the study of elusive wilderness carnivores.--Transactions of the 11th Congress of the International Union of Game Biologists 11: 315-322.
Mech, L.D. & Frenzel, L.D. 1971: Ecological studies of the timber wolf in northeastern Minnesota.--U.S. Department of Agriculture Forest Service, North Central Forest Experiment Station Research Paper NC-52, St. Paul, Minnesota, USA, 62 pp.
Mech, L.D. & Peterson, R.O. 2003: Wolf-Prey Relations.--In: Mech, L.D. & Boitani, L. (Eds.); Wolves: Behavior, Ecology and Conservation. University of Chicago Press, Chicago, Illinois, USA, pp. 131-160.
Messier, F. 1994: Ungulate population-models with predation --a case-study with the North American moose. Ecology 75: 478-488.
Metz, M., Smith, D., Vucetich, J.A., Stahler, D. & Peterson, R.O. 2012: Seasonal patterns of predation for gray wolves in the multi-prey system of Yellowstone National Park.--Journal of Animal Ecology 81(3): 553-563.
Nelson, M.E. & Mech, L.D. 1986: Relationship between snow depth and gray wolf predation on white-tailed deer.--Journal of Wildlife Management 50: 471-474.
Peterson, R.O. & Ciucci, P. 2003: The wolf as a carnivore.--In Mech, L.D. & Boitani, L. (Eds.); Wolves: Behavior, Ecology and Conservation. University of Chicago Press, Chicago, Illinois, USA, pp. 104-130.
Pimlott, D.H., Shannon, J.A. & Kolenosky, G.B. 1969: The ecology of the timber wolf in Algonquin Provincial Park.--Ontario Department of Forest and Lands, Research Report (Wildlife) No. 87, Toronto, Ontario, Canada, 92 pp.
Post, E., Peterson, R.O., Stenseth, N.C. & McLaren, B.E. 1999: Ecosystem consequences of wolf behavioural response to climate.--Nature 401: 905-907.
Potvin, M.J., Drummer, T., Vucetich, J.A., Beyer, D.E., Peterson, R.O. & Hammill, J.H. 2005: Monitoring and habitat analysis for wolves in Upper Michigan.--Journal of Wildlife Management 69: 1660-1669.
Sand, H., Wabakken, P., Zimmermann, B., Johansson, O., Pedersen, H.C. & Liberg, O. 2008: Summer kill rates and predation pattern in a wolf-moose system: can we rely on winter estimates?--Oecologia 156: 53-64.
Schmidt, P.A. & Mech, L.D. 1997: Wolf pack size and food acquisition.--American Naturalist 150: 513-517.
Smith, D.W., Drummer, T.D., Murphy, K.M., Guernsey, D.S. & Evans, S.B. 2004: Winter prey selection and estimation of wolfkill rates in Yellowstone National Park, 1995-2000.--Journal of Wildlife Management 68:153-166.
Stenlund, M.H. 1955: A field study of the timber wolf (Canis lupus) on the Superior National Forest, Minnesota.--Minnesota Department of Conservation Technical Bulletin 4, Minnesota, USA, 55 pp.
Varley, N. & Boyce, M.S. 2006: Adaptive management for reintroductions: updating a wolf recovery model for Yellowstone National Park.--Ecological Modelling 193: 315-339.
Vucetich J.A., Hebblewhite, M., Smith, D.W. & Peterson, R.O. 2011: Predicting prey population dynamics from kill rate, predation rate and predator-prey ratios in three wolf-ungulate systems.--Journal of Animal Ecology 80: 1236-1245.
Vucetich, J.A., Peterson, R.O. & Schaefer, C.L. 2002: The effect of prey and predator densities on wolf predation.--Ecology 83: 3003-3013.
Vucetich, J.A., Peterson, R.O. & Waite, T.A. 2004: Raven scavenging favours group foraging in wolves.--Animal Behaviour 67: 1117-1126.
John A. Vucetich, Brett A. Huntzinger, Rolf O. Peterson & Leah M. Vucetich, School of Forest Resources and Environmental Sciences, Michigan Technological University, Houghton, Michigan 49931, USA--e-mail addresses: email@example.com (John A. Vucetich); firstname.lastname@example.org (Brett A. Huntzinger); email@example.com (Rolf O. Peterson); lmvuceti@mtu. edu (Leah M. Vucetich)
James H. Hammill, Michigan Department of Natural Resources, Crystal Falls, Michigan 49920, USA--e-mail: firstname.lastname@example.org
Dean E. Beyer, Jr., Michigan Department of Natural Resources, Northern Michigan University, Marquette, Michigan 49855, USA--e-mail: email@example.com
Corresponding author: John A. Vucetich
Received 15 July 2011, accepted 6 April 2012
Associate Editor: Stefano Focardi
Table 1. Upper Michigan kill rate data summarized for each wolf pack (1-5) and winter (2001-2004) with [greater than or equal to] 4 days of observation. 'Days' refers to number of days available from which to estimate kill rate. Kg/wolf/day is calculated based on the assumption that the consumable mass from each deer was 54.5 kg (Schmidt & Mech 1997). Year Pack ID# Pack size Days Kills/pack/day Kg/wolf/day 2001 2 5 14 0.29 3.16 2001 1 5 42 1.05 11.45 2002 3 4 21 0.24 3.27 2002 1 3 38 0.94 17.08 2002 2 8 43 1.16 7.90 2003 4 3 31 0.45 8.18 2003 2 8 32 0.34 2.32 2003 1 3 36 0.89 16.17 2004 2 9 18 0.50 3.03 2004 5 6 26 1.19 10.81 2004 3 5 27 0.22 2.40 2004 4 7 31 0.84 6.54 Table 2. Estimates of daily kill rates (in kg/wolf/day) for wolf packs of different size (3-9) which primarily eat moose. Values were determined by best-fit Poisson regression analysis of data previously compiled by Schmidt & Mech (1997) (i.e. equation (4)). Pack size Duration of study in days 10 30 60 120 3 18.6 13.2 10.6 8.5 5 16.4 11.6 9.3 7.5 9 12.7 9.0 7.2 5.8 Table 3. Daily kill rate of wolf packs, total number of daily kill rates and of deer killed, and the number of kills from which wolves did not bed within 500 m during daylight hours during the four winters (2001-2004). Daily kill rate 0 1 2+ Total Number of daily kill rates 198 115 58 371 Number of deer killed 0 115 159 274 Kills from which wolves 0 38 103 141 bedded > 500 m away Percentage of kills wolves 0 33 65 51 did not bed near Table 4. Bias in estimates of kill rate when days with multiple kills are recorded as days with only a single kill. As winter (Julian day) progresses, so does the percentage of the kills occurring on multiple kill days. Average kill rate (with pooled pack sizes) is the estimate that accounts for days with multiple kills, and bias in kill rate estimation is the proportion by which kill rate would have been underestimated if multiple kills had been mistaken for only single kills. Numbers in parentheses are sample sizes (i.e. number of days) associated with each estimate. Period Average Kill rate if Bias in kill kill rate multiple kill rate estimation (deer/pack/day) days were recorded as singles December-January 0.32 0.32 0% (N = 91) February 0.66 0.46 30% (N = 146) March-April 0.95 0.52 45% (N = 159)
Please note: Some tables or figures were omitted from this article.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Original article|
|Author:||Vucetich, John A.; Huntzinger, Brett A.; Peterson, Rolf O.; Vucetich, Leah M.; Hammill, James H.; Be|
|Date:||Sep 1, 2012|
|Previous Article:||Model evaluation in statistical population reconstruction.|
|Next Article:||Impacts of rural development on Yellowstone wildlife: linking grizzly bear Ursus arctos demographics with projected residential growth.|