Printer Friendly

EFFECTS OF SCAVENGING ON ASSUMPTIONS OF MORTALITY ANALYSES OF RADIO-MARKED GAMEBIRDS.

Survival is an important determinant of wildlife population biology. Over the past few decades, estimation of survival has improved in both insightfulness and rigor (Lebreton and others 1992; Murray 2006). In particular, the modeling of survival relative to intrinsic and environmental conditions has received significant attention (Murray and Patterson 2006). Recent development of time-dependent survival analyses has improved our ability to understand the demographic consequences of individual behaviors by establishing clear linkages between habitat use and survival (Johnson and others 2004; Halstead and others 2012). Use of models such as shared frailty models and the Andersen-Gill formulation of the Cox proportional hazards model has rapidly increased due to several advantages over other known-fate survival models, including the ability to incorporate spatial covariates on cause-specific mortality risk (Andersen and Gill 1982; Johnson and others 2004; Murray and Patterson 2006; Fleming and Harrington 2011; Halstead and others 2012; Coates and others 2017). These models have been used in studies of multiple taxa (Smith and others 2010; White and others 2010; Halstead and others 2012; Taylor and others 2016; Senner and others 2017), including gamebirds ([Angstrom]hlen and others 2013; Dinkins and others 2014; Winder and others 2014; Coates and others 2017).

Data for time-dependent survival models are structured so that a single animal has multiple observations, each of which represents a different interval of risk (Halstead and others 2012; Therneau and Grambsch 2013). For individuals that die, the final interval is represented by the mortality time and is compared to a series of intervals during which the live animal was tracked and which represent nonmortality locations. Mortality events are discovered or represented as a binomial status at the end of each interval, so the conditions influencing mortality risk are typically measured at that end location (Johnson and others 2004). Previous work found slight differences in coefficient estimates for habitat conditions affecting mortality risk depending on when the conditions for a given nonmortality interval were measured, but no systematic bias (Johnson and others 2004), suggesting that the habitat variables related to nonmortality locations adequately represent conditions influencing mortality risk. However, there have been no similar studies testing the assumption that habitat conditions at carcass locations represent the habitat conditions influencing mortality risk. If a carcass is moved after death, by either the original predator or a scavenger, then inferences based on the habitat conditions measured at the mortality location will be biased. Moreover, if scavenging activity confounds determination of cause-specific mortality, estimates of cause-specific mortality risks may result in incorrect inferences and inappropriate management recommendations.

We used radio-marked Chukar (Alectoris chukar) carcasses to test the common assumptions that scavenging activity does not confound field-based determination of cause-specific mortality, and mortality locations from telemetry studies of gamebirds accurately represent the habitat characteristics where individual animals die. Given relatively high but seasonally changing occurrence of both avian and mammalian scavengers at our study area (Vold 2018), we hypothesized that: (1) a non-zero proportion of carcasses would exhibit signs of scavenging within a week of being placed, thus resulting in incorrect classification of cause-specific mortality; (2) scavenged carcasses would be moved from their original locations; (3) the probability that a carcass is scavenged would be mediated by local habitat conditions; and (4) seasonality would influence scavenging activity and potential biases associated with survival estimation.

METHODS

This study was conducted in eastern Montana, 16 km south of Sidney, Montana (centered on UTM: Zone 13N, 0571465E, 5268907N), in April 2017 and July 2018, to assess scavenging activity in both cool (April, spring) and warm (July, summer) seasons. Our study area was located in the northern mixed-grass prairie, within an area typical for the region dominated by contiguous rangelands and managed for livestock grazing. The study area was primarily Great Plains mixed-grass prairie, interspersed with Great Plains badlands and wooded draws and ravines (LANDFIRE 2013). Average annual precipitation was 35.3 cm, and average high and low temperatures were 14.7 and -1.0[degrees]C in April and 29.4 and 12.7[degrees]C in July. Based on both point-count surveys and camera-trap data collected during 2016-2017, potential scavengers in our study area during the time of this experiment included Coyote (Canis latrans), Northern Raccoon (Procyon lotor), American Badger (Taxidea taxus), Striped Skunk (Mephitis mephitis), Shorttailed Weasel (Mustela ermine), Peregrine Falcon (Falco peregrinus), Prairie Falcon (Falco mexicanus), Northern Harrier (Circus cyaneus), Redtailed Hawk (Buteo jamaicensis), Golden Eagle, (Aquila chrisaetos), Great Horned Owl (Bubo virginianus), American Crow (Corvus brachynrhynchos), Black-billed Magpie (Pica hudsonia), and Turkey Vulture (Cathartes aura) (Vold 2018). All of the mammalian species are residents and, of the avian species, the Prairie Falcon, Golden Eagle, Great Horned Owl, and Black-billed Magpie are resident, while the remaining avian species are migratory, passing through the study area beginning in April, with some individuals remaining to breed.

Twenty-four Chukar carcasses were obtained in both April 2017 and July 2018. Carcasses used in 2017 were obtained from a National Shoot to Retrieve Trial in Savage, Montana; all birds were killed by hunters using shotguns. Hunter-killed carcasses were not available in 2018, so we obtained frozen carcasses directly from a private game bird production facility. Carcasses were outfitted with 18-g VHF radio-transmitters with an elastic necklace harness (model A4050; Advanced Telemetry Systems, Isanti, MN), and placed at randomly selected locations with a minimum distance of 350 m between points to satisfy the assumption of site independence (Lesmeister and others 2015). Carcasses were then checked at 3 and 7 d post-placement, and any scavenging activity and movement of either the carcass or the transmitter were recorded. Intervals of 3 and 7 d were chosen to represent intervals between consecutive locations used in many studies that were unable to monitor study animals daily (Murray 2006; Dinkins and others 2014; Winder and others 2018) and so represent relative time periods during which scavengers would have an opportunity to scavenge carcasses. Scavenging activity was classified as either mammal, avian, or unknown scavenger (Elbroch 2003). Scavenging activity from invertebrates was not recorded as it is not relevant to studies of gamebird mortality. Scavenging was classified as mammalian if bite marks, chewed feathers, or mammalian tracks or scat were present (Elbroch 2003). Scavenging was classified as avian if the carcass had been decapitated or cleaned of the breast muscle with no bite marks, or if the feathers had only been plucked (Elbroch 2003). Scavenging was classified as an unknown scavenger if the entire carcass was removed (only an unmarked transmitter remained). All Chukars used in this experiment were intact carcasses with no previous signs of damage other than gunshots, and all carcasses were handled similarly and placed in the field within 24 h of being obtained.

Habitat conditions were measured at both the initial locations where carcasses were placed and their locations at 3 d post-placement if the carcass had been moved >12 m. In our study, a distance of 12 m represents the minimum distance at which the subsequent local habitat conditions would differ from those measured at the initial carcass location (MC Milligan, unpubl. data). We recorded visual obstruction readings (VOR) at the carcass location and at 4 points 6 m from the carcass in each cardinal direction. At each point, VOR was measured in each cardinal direction from a distance of 2 m and a height of 0.5 m using a Robel pole (Robel and others 1970). We also estimated non-overlapping vegetation cover (percent new grass, residual grass, forbs, shrubs, bare ground, and litter) at 12 subsampling locations within 6 m of the carcass using a 20 x 50 cm sampling frame to characterize the vegetation composition (Daubenmire 1959). In addition to vegetation conditions, we measured habitat variables previously found to be important for prairie grouse (Tympanuchus spp.) survival, including distances to houses, oil pads, and roads, and 2 measures of habitat fragmentation: distance to grassland patch edge and density of edge habitat (Hovick and others 2014). Variables representing anthropogenic disturbances were based on road data-sets from Montana and North Dakota (Montana State Library, North Dakota GIS Hub Data Portal) and landcover analyses utilized the 30 m resolution LANDFIRE data depicting vegetation type (LANDFIRE 2013).

We 1st evaluated models examining the effect of season (April = cool-spring; July = warm-summer) on the probability that a carcass was scavenged during the 3-d interval. We assumed that the importance of season relative to scavenging activity could vary depending on the length of the relocation interval. We then separately evaluated competing models examining the effects of small-scale vegetation and habitat risk factors on the probability that a carcass was scavenged or moved during the 7-d interval. We also tested for additive effects with habitat risk factors and season. We assumed that habitat factors affecting the probability of scavenging would not differ between the 3-d and 7-d intervals. We used logistic regression to evaluate both the probability of scavenging and the probability that a carcass was moved [greater than or equal to] 12 m. Models were compared using Akaike's Information Criterion adjusted for small sample sizes ([AIC.sub.c]), and models with both large model weights ([w.sub.i] [greater than or equal to] 0.1) and [AIC.sub.c]- values [less than or equal to]2 from the best-fit model were considered equally parsimonious (Burnham and others 2011). Effects of variables were considered important if 85% confidence intervals of model-averaged slope coefficients did not overlap zero (Arnold 2010).

RESULTS

During the first 3-d post-placement, 7 of 24 carcasses (29%) were scavenged in April (6 by birds and 1 by a mammal), and 3 of 24 carcasses (12.5%) were scavenged in July (1 by birds and 2 by mammals). At 7-d post-placement, 12 of 24 carcasses (50%) were scavenged in April (10 by birds and 2 by mammals), and 19 of 24 carcasses (79.2%) were scavenged in July (2 by birds, 10 by mammals, and 7 by unknown scavengers). Only 2 carcasses (8%) were moved [greater than or equal to]12 m during the initial 3-d post-placement, with both being moved in July (Fig. 1). Two carcasses (8%) and 12 carcasses (50%) were moved [greater than or equal to]12 m during the 7-d post-placement in April and July, respectively. Both carcasses that were moved by 7-d post-placement in April were scavenged by birds and were moved 97 and 150 m, respectively. Half (50%) of the carcasses moved during the 7-d post-placement in July were scavenged by mammals, with the remainder scavenged by unknown scavengers. The carcasses placed in July were moved between 61 and 923 m, with a mean distance of 306 m (Fig.1). There was no significant difference in the probability of scavenging during the 3-d interval between April and July (null model: [DELTA][AIC.sub.c] = 0, [w.sub.i] = 0.51), but during the 7-d interval, the probability that a carcass was scavenged was higher during July than April (Table 1).

Carcasses were initially placed across a range of habitat conditions (Fig. 2), but local vegetation conditions did not influence the probability that a carcass was scavenged (Table 1). Rather, the probability that a carcass was scavenged during the 7-d interval was best predicted by season (cool-spring or warm-summer), which, in combination with other variables, accounted for 76% of model weight (Table 1; Fig. 3). The probability of scavenging was higher during July than April ([beta] = 1.37 [+ or -] 0.66). There was also evidence for an effect of anthropogenic structures on the probability that a carcass was scavenged (Table 1), with carcasses located closer to oil pads ([beta] = 0.03 [+ or -] 0.02; Fig. 3) and further from houses ([beta] = 0.06 [+ or -] 0.04) more likely to be scavenged. The odds of scavenging decreased 2.8% for every 100 m from an oil pad, and increased 6.6% with every 100 m from a house. The 85% confidence intervals on effects of all other variables overlapped zero and so were not considered important (Appendix Table S1). Four models describing the probability that a carcass was moved were supported by the data based on model weights >0.1 and [AIC.sub.c] values [less than or equal to]2 from the best-fit model (Table 2), but only the effects of distance to road, distance to oil pad, and distance to house were considered significant, with 85% confidence intervals that did not overlap zero (Appendix Table S2). The probability that a carcass was moved within 7 d was positively related to distance to the nearest road ([beta] = 0.08 [+ or -] 0.04) and distance to house ([beta] = 0.09 [+ or -] 0.05), and negatively related to distance to the nearest oil pad ([beta] = -0.04 [+ or -] 0.02). The odds of scavenging decreased 3.5% for every 100 m increase in distance from oil pads and increased 8.4% and 9.2% for every 100 m increase in distance from roads and houses, respectively. Although low sample sizes precluded statistical analyses, habitat conditions were similar between the original locations and the locations where transmitters were recovered (Fig. 4).

DISCUSSION

Incorrect predator classification can introduce bias into telemetry studies when scavengers obscure evidence from the original predator, and inflate estimates of mortality due to predation if animals that die due to factors other than predation, such as disease, are misclassified as a result of scavenging activity (Larsen and others 2008). Similar to previous studies, we observed high levels of scavenging activity on gamebird carcasses in eastern Montana that would result in incorrect classification of cause-specific mortality (Dumke and Pils 1973; Bumann and Stauffer 2002; Larsen and others 2008), with rates of scavenging that were significantly greater in summer than spring.

Scavenging activity did not significantly influence perceived mortality locations regardless of season or local habitat conditions with a relocation interval of 3 d, but mortality locations may be biased over longer periods, particularly in warm seasons or regions. Previous studies found that radio-tagged carcasses were moved minimal distances, even over longer time intervals (Bumann and Stauffer 2002; Larsen and others 2008). Furthermore, although sample sizes precluded analysis, habitat conditions were similar between the original locations and locations where transmitters were recovered. Therefore, for gamebird telemetry studies with a maximum relocation interval of 3 d, perceived habitat conditions at mortality locations should accurately represent those influencing mortality risk in the northern mixed-grass prairie, but caution should be exercised in studies with relocations intervals >3 d.

In our study, local vegetation conditions did not influence whether a carcass was scavenged. Previous work has found that factors including habitat, carcass size, temperature, and changes in the vertebrate scavenger community can all influence scavenging activity, although the importance of each factor varies by study (Selva and others 2005; Olson and others 2012; Moleon and others 2015; Turner and others 2017). While some studies have found specific habitat types to be important predictors of scavenging activity (Selva and others 2005; Turner and others 2017), other studies have found no influence of factors such as habitat structure, landscape connectivity or distance to edge habitat (Bumann and Stauffer 2002; Olson and others 2016), or have shown that other factors such as carcass type, temperature, or the local predator community were more important predictors (Olson and others 2012; Moleon and others 2015; Olson and others 2016). The habitat variables measured in our study were not categorical habitat types (for example see Selva and others 2005), but rather fine-scale continuous measures of vegetation structure, and differences in habitat measures could explain why our results differed from those of previous research.

In contrast, proximity to anthropogenic structures influenced both the probability of a carcass being scavenged and the probability that a carcass was moved. Carcasses located in closer proximity to oil pads were more likely to be scavenged and moved, while carcasses located close to roads and houses were less likely to be scavenged or moved. Anthropogenic structures related to energy development, including oil pads, have been associated with lower survival of grouse (Holloran 2005; Hovick and others 2014). Raptors are visual predators and utilize anthropogenic infrastructure as perch sites, which can facilitate the detection of both live prey and carcasses to scavenge (Prather and Messmer 2010; Slater and Smith 2010). More-over, energy development increases habitat fragmentation (Hovick and others 2014), and mesocarnivores often exhibit positive responses to habitat fragmentation (Andren 1995; Winter and others 2000). Carcasses were also less likely to be moved when in close proximity to roads and houses, potentially reflecting avoidance of human presence, a trend found in previous studies (Burr and others 2017; Vold 2018). In our study area, oil pads were not being actively drilled and so human presence was minimal.

Similar to other studies, (Bumann and Stauffer 2002; Olson and others 2016; Turner and others 2017), estimates of scavenging activity during the 7-d interval differed between seasons, with greater scavenging activity during the warmer summer season. This suggests that temperature could be a significant factor influencing whether a carcass was discovered and may reflect differences in the scavenger community between seasons. Both the abundance and composition of the local scavenger community can have an important influence on scavenging activity (Turner and others 2017), and the community of scavengers appeared to shift between the seasons, although sample sizes precluded statistical analyses. Seasonal changes in predator and scavenger communities have important implications for the determination of cause-specific mortality. In our study area, the majority of carcasses were scavenged by birds in April, while scavenging by mammals was more prevalent in July, which could be related to the difference in carcass movement between the 2 seasons. Avian scavengers in our study area, such as the Northern Harrier, typically pluck and consume prey on the ground, whereas it could be easier for larger-bodied mammals to move or cache a carcass (Harrington 1982; Thogmartin and Schaeffer 2000; Smith and others 2011). Misclassification of cause-specific mortality may be higher in areas where mammals are the primary scavengers (Bumann and Stauffer 2002). Therefore, both the scavenger assemblage and the timing of the study should be considered when applying our results to other locations. The biases associated with scavenging can vary throughout the year, with higher avian scavenging during migration and with greater scavenging by olfactory mammals during warmer seasons.

Our experiment only evaluates the effects of scavenging on the perceived mortality locations of gamebirds and does not include any movement of the carcass post-kill by the initial predator. The importance of this initial movement could vary depending on the predator assemblage in the area and, based on the scavenging activity in our study, could potentially be less pronounced in areas or for species where the dominant predators are avian. However, in areas with a higher proportion of mammalian scavengers, including well-known caching species such as Mountain Lions (Puma concolor) or woodrats (Neotoma spp.), this movement could be more significant (Bauer and others 2005; Larsen and others 2008).

Our study provides provisional evidence that scavenging can significantly confound field-based determination of cause-specific mortality, but, depending on the season and relocation interval of the study, may not have a significant influence on mortality locations of gamebirds in the northern mixed-grass prairie. However, both season and local vertebrate assemblages should be considered when determining the optimal tracking interval for telemetry studies if modeling the effect of habitat conditions on survival is an important priority. Biologists should be cautious when making inferences about cause-specific mortality when scavengers occupy their study areas.

ACKNOWLEDGMENTS

We would like to thank M Jepsen and the National Shoot to Retrieve Association for providing us with Chukar carcasses used in this study. Funding was provided by the Montana Agricultural Experiment Station, the Montana Department of Fish, Wildlife and Parks, and the US Fish and Wildlife Service Sportfish and Wildlife Restoration Program.

LITERATURE CITED

AHLEN PA, WILLEBRAND T, SJOBERG K, HORNELL-WILLE-BRAND M. 2013. Survival of female capercaillie Tetrao urogallus in northern Sweden. Wildlife Biology 19:368-373.

ANDERSEN PK, GILL RD. 1982. Cox's regression model for counting processes: A large sample study. Annals of Statistics 10:1100-1120.

ANDREN H. 1995. Effects of landscape composition on predation rates at habitat edges. In: Hansson L, Fahrig L, Merriam G, editors. Mosaic landscapes and ecological processes. New York, NY: Chapman and Hall. p 225-255.

ARNOLD TW. 2010. Uninformative parameters and model selection using Akaike's Information Criterion. Journal of Wildlife Management 74:1175-1178.

BAUER JW, LOGAN KA, SWEANOR LL, BOYCE WM. 2005. Scavenging behavior in Puma. Southwestern Naturalist 50:466-471.

BUMANN GB, STAUEEER DF. 2002. Scavenging of Ruffed Grouse in the Appalachians: Influences and implications. Wildlife Society Bulletin 30:853-860.

BURNHAM KP, ANDERSON DR, HUYVAERT KP. 2011. AIC model selection and multimodel inference in behavioral ecology: Some background, observations, and comparisons. Behavioral Ecology and Sociobiology 65:23-35.

BURR PC, ROBINSON AC, LARSEN RT, NEWMAN RA, ELLIS-FELEGE SN. 2017. Sharp-tailed Grouse nest survival and nest predator habitat use in North Dakota's Bakken oil field. PloS One 12:e0170177.

COATES PS, PROCHAZKA BG, RICCA MA, GUSTAESON KB, ZIEGLER P, CASAZZA ML. 2017. Pinyon and Juniper encroachment into sagebrush ecosystems impacts distribution and survival of Greater Sage-grouse. Rangeland Ecology and Management 70:25-38.

DAUBENMIRE RF. 1959. Canopy coverage method of vegetation analysis. Northwest Science 33:43-64.

DINKINS JB, CONOVER MR, KIROL CP, BECK JL, FREY SN. 2014. Greater Sage-grouse (Centrocercus urophasianus) hen survival: Effects of raptors, anthropogenic and landscape features, and hen behavior. Canadian Journal of Zoology 92:319-330.

DUMKE RT, PILS CM. 1973. Mortality of radio-tagged pheasants on the Waterloo Wildlife Area. Wisconsin Department of Natural Resources, Technical Bulletin 72.

ELBROCH M. 2003. Mammal tracks and sign: A guide to North American species. Mechanicsburg, PA: Stackpole.

FLEMING TR, HARRINGTON DP. 2011. Counting processes and survival analysis. Hoboken, NJ: John Wiley and Sons.

HALSTEAD B, WYLIE G, COATES P, VALCARCEL P, CASAZZA M. 2012. Bayesian shared frailty models for regional inference about wildlife survival. Animal Conservation 15:117-124.

HARRINGTON FH. 1982. Urine marking at food and caches in captive Coyotes. Canadian Journal of Zoology 60:776-782.

HOLLORAN MJ. 2005. Greater Sage-grouse (Centrocercus urophasianus) population response to natural gas field development in western Wyoming [dissertation]. Laramie, WY: University of Wyoming.

HOVICK TJ, ELMORE RD, DAHLGREN DK, FUHLENDORE SD, ENGLE DM. 2014. Evidence of negative effects of anthropogenic structures on wildlife: A review of grouse survival and behaviour. Journal of Applied Ecology 51:1680-1689.

JOHNSON CJ, BOYCE MS, SCHWARTZ CC, HAROLDSON MA. 2004. Modeling survival: Application of the Andersen--Gill model to Yellowstone Grizzly Bears. Journal of Wildlife Management 68:966-978.

LANDFIRE. 2013. LANDFIRE existing vegetation type layer. US Department of Interior, Geological Survey. Available: http://www.landfire.gov/index. php. Accessed 18 November 2016.

LARSEN RT, BENTLEY DF, FLINDERS JT. 2008. Implications of woodrats and other scavengers for avian telemetry studies. Journal of Wildlife Management 72:1152-1155.

LEBRETON J-D, BURNHAM KP, CLOBERT J, ANDERSON DR. 1992. Modeling survival and testing biological hypotheses using marked animals: A unified approach with case studies. Ecological Monographs 62:67-118.

LESMEISTER DB, NIELSEN CK, SCHAUBER EM, HELLGREN EC. 2015. Spatial and temporal structure of a mesocarnivore guild in midwestern North America. Wildlife Monographs 191:1-61.

MOLEON M, SANCHEZ-ZAPATA JA, SEBASTIAN-GONZALEZ E, OWEN-SMITH N. 2015. Carcass size shapes the structure and functioning of an African scavenging assemblage. Oikos 124:1391-1403.

MURRAY DL. 2006. On improving telemetry-based survival estimation. Journal of Wildlife Management 70:1530-1543.

MURRAY DL, PATTERSON BR. 2006. Wildlife survival estimation: Recent advances and future directions. Journal of Wildlife Management 70:1499-1503.

OLSON Z, BEASLEY J, DEVAULT TL, RHODES O. 2012. Scavenger community response to the removal of a dominant scavenger. Oikos 121:77-84.

OLSON ZH, BEASLEY JC, RHODES OE JR. 2016. Carcass type affects local scavenger guilds more than habitat connectivity. PloS One 11:e0147798.

PRATHER PR, MESSMER TA. 2010. Raptor and corvid response to power distribution line perch deterrents in Utah. Journal of Wildlife Management 74:796-800.

SELVA N, JEDRZEJEWSKA B, JEDRZEJEWSKI W, WAJRAK A. 2005. Factors affecting carcass use by a guild of scavengers in European temperate woodland. Canadian Journal of Zoology 83:1590-1601.

SENNER NR, STAGER M, SANDERCOCK BK. 2017. Ecological mismatches are moderated by local conditions for two populations of a long-distance migratory bird. Oikos 126:61-72.

SLATER SJ, SMITH JP. 2010. Effectiveness of raptor perch deterrents on an electrical transmission line in southwestern Wyoming. Journal of Wildlife Management 74:1080-1088.

SMITH DW, BANGS EE, OAKLEAE JK, MACK C, FONTAINE J, BOYD D, JIMENEZ M, PLETSCHER DH, NIEMEYER CC, MEIER TJ. 2010. Survival of colonizing Wolves in the northern Rocky Mountains of the United States, 1982-2004. Journal of Wildlife Management 74:620-634.

SMITH KG, WITTENBERG SR, MACWHIRTER RB, BILDSTEIN KL. 2011. Northern Harrier (Circus cyaneus), version 2.0. In: Rodewald PG, editor. The birds of North America. Ithaca, NY: Cornell Lab of Ornithology.

TAYLOR KL, BECK JL, HUZURBAZAR SV. 2016. Factors influencing winter mortality risk for Pronghorn exposed to wind energy development. Rangeland Ecology and Management 69:108-116.

THERNEAU TM, GRAMBSCH PM. 2013. Modeling survival data: Extending the Cox model. Berlin, Germany: Springer Science and Business Media.

THOGMARTIN WE, SCHAEEFER BA. 2000. Landscape attributes associated with mortality events of wild Turkeys in Arkansas. Wildlife Society Bulletin 28:865-874.

TURNER KL, ABERNETHY EF, MIKE CONNER L, RHODES OE, BEASLEY JC. 2017. Abiotic and biotic factors modulate carrion fate and vertebrate scavenging communities. Ecology 98:2413-2424.

VOLD ST. 2018. Effects of livestock grazing management on the ecology of grassland birds and their predators in a northern mixed-grass prairie ecosystem [thesis]. Bozeman, MT: Montana State University.

WHITE CG, ZAGER P, GRATSON MW. 2010. Influence of predator harvest, biological factors, and landscape on Elk calf survival in Idaho. Journal of Wildlife Management 74:355-369.

WINDER VL, MCNEW LB, GREGORY AJ, HUNT LM, WISELY SM, SANDERCOCK BK. 2014. Effects of wind energy development on survival of female Greater Prairie-Chickens. Journal of Applied Ecology 51:395-405.

WINDER VL, MCNEW LB, PITMAN JC, SANDERCOCK BK. 2018. Effects of rangeland management on survival of female Greater Prairie-Chickens. Journal of Wildlife Management 82:113-122.

WINTER M, JOHNSON DH, FAABORG J. 2000. Evidence for edge effects on multiple levels in tallgrass prairie. The Condor 102:256-266.

Submitted 25 February 2019, accepted 24 June 2019. Corresponding Editor: Robert Hoffman.
Appendix S1. Model averaged coefficient estimates for variables
relative to the probability of scavenging of Chukar carcasses
experimentally placed in eastern Montana in April 2017 and July 2018
during the 7-d interval post-placement. The estimate, standard error
(SE), and lower (LCI) and upper (UCI) 85% confidence intervals are
reported.

Variable                     Estimate   SE       LCI       UCI

LOCAL VEGETATION
VOR                          -0.1525    0.2432   -0.5025    0.1976
% New grass                   0.0144    0.0220   -0.0173    0.0460
% Residual                   -0.0267    0.0189   -0.0540    0.0006
% Forb                        0.0168    0.0442   -0.0467    0.0804
% Shrub                       0.0143    0.0203   -0.0149    0.0435
% Litter                     -0.0505    0.0445   -0.1145    0.0135
% Bare                        0.0019    0.0190   -0.0254    0.0292
HABITAT RISK FACTORS
Distance to house             0.0636    0.0439    0.0004    0.1268
Edge density                 -0.0049    0.0058   -0.0134    0.0035
Distance to grassland edge    0.3156    0.5226   -0.4367    1.0680
Distance to oil pad          -0.0285    0.0192   -0.0561   -0.0009
Distance to road              0.0090    0.0343   -0.0403    0.0583
Season                        1.3734    0.6603    0.4228    2.3240

Appendix S2. Model averaged coefficient estimates for variables
relative to the probability that Chukar carcasses experimentally placed
in eastern Montana in April 2017 and July 2018 were moved by
scavengers. The estimate, standard error (SE), and lower (LCI) and
upper (UCI) 85% confidence intervals are reported.

Variable                     Estimate   SE       LCI       UCI

LOCAL VEGETATION
VOR                          -0.3563    0.2569   -0.7262    0.0136
% New grass                  -0.0088    0.0218   -0.0402    0.0226
% Residual                    0.0004    0.0188   -0.0267    0.0276
% Forb                       -0.0229    0.0442   -0.0865    0.0407
% Shrub                      -0.0063    0.0176   -0.0316    0.0190
% Litter                      0.0661    0.0519   -0.0086    0.1409
% Bare                        0.0115    0.0214   -0.0193    0.0424
HABITAT RISK FACTORS
Distance to house             0.0881    0.0460    0.0219    0.1544
Edge density                  0.0011    0.0059   -0.0074    0.0096
Distance to grassland edge    0.3684    0.5336   -0.3998    1.1366
Distance to oil pad          -0.0356    0.0198   -0.0641   -0.0070
Distance to road              0.0807    0.0374    0.0268    0.1345
Season                       -0.6393    0.6617   -1.5919    0.3132


MEGAN C MILLIGAN AND LANCE B MCNEW

Department of Animal and Range Sciences, Montana State University, 103 Animal Biosciences Building, Bozeman, MT 59717-2900 USA; megan.milligan11@gmail.com
TABLE 1. Support for candidate models predicting the probability of
scavenging of Chukar carcasses experimentally placed in eastern Montana
in April 2017 and July 2018 during the 7-d interval post-placement.
Results for analyses examining effects of small-scale vegetation and
habitat risk factors are included. The number of parameters (K),
[AIC.sub.c] values, [DELTA][AIC.sub.c]- values, model weights
([w.sub.i]) and log-likelihoods are reported.

Model                  K   [AIC.sub.c]   [DELTA][AIC.sub.c]  [AIC.sub.c]
                                                             [w.sub.i]

LOCAL VEGETATION
Null                   1   63.60         0.00                0.22
% Residual             2   63.71         0.11                0.21
% Litter               2   64.54         0.94                0.14
% Shrub                2   65.23         1.63                0.10
% New grass            2   65.35         1.75                0.09
VOR                    2   65.39         1.79                0.09
% Forb                 2   65.64         2.04                0.08
Bare                   2   65.78         2.17                0.07
HABITAT RISK FACTORS
Season + Distance      3   61.96         0.00                0.19
to oil pad
Season + Distance to   3   62.09         0.13                0.18
house
Season                 2   62.10         0.14                0.18
Season + Edge density  3   63.64         1.68                0.08
Season + Distance to   3   63.99         2.03                0.07
grassland edge
Season + Distance to   3   64.31         2.35                0.06
road
Dist. To oil pad       2   64.48         2.52                0.06
Null                   1   64.49         2.53                0.05
Distance to house      2   64.60         2.64                0.05
Edge density           2   66.00         4.03                0.03
Distance to grassland  2   66.31         4.35                0.02
edge
Distance to road       2   66.60         4.64                0.02

Model                   Cum. [w.sub.i]   LL

LOCAL VEGETATION
Null                    0.22             -30.76
% Residual              0.43             -29.72
% Litter                0.57             -30.13
% Shrub                 0.66             -30.48
% New grass             0.76             -30.54
VOR                     0.85             -30.56
% Forb                  0.93             -30.68
Bare                    1.00             -30.75
HABITAT RISK FACTORS
Season + Distance       0.19             -27.71
to oil pad
Season + Distance to    0.38             -27.77
house
Season                  0.56             -28.92
Season + Edge density   0.64             -28.55
Season + Distance to    0.71             -28.72
grassland edge
Season + Distance to    0.77             -28.88
road
Dist. To oil pad        0.83             -30.11
Null                    0.88             -31.20
Distance to house       0.93             -30.17
Edge density            0.96             -30.86
Distance to grassland   0.98             -31.02
edge
Distance to road        1.00             -31.17

TABLE 2. Support for candidate models predicting the probability that
Chukar carcasses experimentally placed in eastern Montana in April 2017
and July 2018 were moved by scavengers. Results for analyses examining
effects of small-scale vegetation and habitat risk factors are
included. The number of parameters (K), [AIC.sub.c] values,
[DELTA][AIC.sub.c] values, model weights ([w.sub.i]) and
log-likelihoods are reported.

Model                    K  [AIC.sub.c]  [DELTA][AIC.sub.c]  [AIC.sub.c]
                                                             [w.sub.i]

LOCAL VEGETATION
Null                     1  60.95        0.00                0.22
VOR                      2  61.10        0.15                0.20
% Litter                 2  61.33        0.38                0.18
% Bare                   2  62.82        1.87                0.09
% Forb                   2  62.87        1.92                0.08
% New grass              2  62.98        2.02                0.08
% Shrub                  2  63.01        2.06                0.08
% Residual               2  63.14        2.18                0.07
HABITAT RISK FACTORS
Dist. To road            2  58.79        0.00                0.27
Dist. To house           2  59.81        1.02                0.16
Season + Dist.           3  60.09        1.30                0.14
To road
Dist. To oil pad         2  60.34        1.55                0.12
Season + Dist. To house  3  61.13        2.34                0.08
Season + Dist. To        3  61.67        2.88                0.06
oil pad
Null                     1  61.71        2.92                0.06
Season                   2  63.01        4.22                0.03
Dist. To grassland edge  2  63.39        4.60                0.03
Edge density             2  63.85        5.06                0.02
Season + Dist. To        3  64.78        5.99                0.01
grassland edge
Season + Edge density    3  65.26        6.46                0.01

Model                    Cum. [w.sub.i]   LL

LOCAL VEGETATION
Null                     0.22             -29.43
VOR                      0.42             -28.41
% Litter                 0.60             -28.53
% Bare                   0.69             -29.28
% Forb                   0.77             -29.30
% New grass              0.85             -29.35
% Shrub                  0.93             -29.37
% Residual               1.00             -29.43
HABITAT RISK FACTORS
Dist. To road            0.27             -27.26
Dist. To house           0.43             -27.77
Season + Dist.           0.57             -26.77
To road
Dist. To oil pad         0.69             -28.04
Season + Dist. To house  0.77             -27.29
Season + Dist. To        0.83             -27.56
oil pad
Null                     0.90             -29.81
Season                   0.93             -29.37
Dist. To grassland edge  0.95             -29.56
Edge density             0.98             -29.79
Season + Dist. To        0.99             -29.12
grassland edge
Season + Edge density    1.00             -29.35


[Please note: Some non-Latin characters were omitted from this article]
COPYRIGHT 2019 Society for Northwestern Vertebrate Biology
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2019 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Milligan, Megan C.; Mcnew, Lance B.
Publication:Northwestern Naturalist: A Journal of Vertebrate Biology
Geographic Code:1U8MT
Date:Dec 22, 2019
Words:5570
Previous Article:AMPHIBIAN MOVEMENT ACROSS A NEW RESIDENTIAL ROAD IN WESTERN WASHINGTON STATE.
Next Article:REOCCURRENCE OF THE UMPQUA CHUB (OREGONICHTHYS KALAWATSETI) IN THE NORTH UMPQUA AFTER ALMOST A CENTURY.
Topics:

Terms of use | Privacy policy | Copyright © 2020 Farlex, Inc. | Feedback | For webmasters