A fine-scale examination of parturition timing in temperate ungulates.

Parturition timing has long been a topic of interest in ungulate research. However, few studies have examined parturition timing at fine scale (e.g., <1 day). Predator activity and environmental conditions can vary considerably with diel timing, which may result in selective pressure for parturition to occur during diel times that maximize the likelihood of neonate survival. We monitored parturition events and early-life survival of elk (Cervus canadensis) and mule deer (Odocoileus hemionus) in Utah, USA to better understand diel timing of parturition in temperate ungulates. Diel timing of parturition was moderately synchronous among conspecifics and influenced by environmental variables on the date of parturition. For elk, parturition events were most common during the morning crepuscular period and generally occurred later (i.e., closer to 12:00) when a relatively large proportion of the moon was illuminated. For mule deer, parturition events were most common during the diurnal period and generally occurred later (i.e., closer to 15:00) on cold, wet dates. Diel timing of parturition did not influence neonate survival, but larger datasets may be required to verify the apparent lack of influence. Although additional work could evaluate alternative variables that might affect parturition timing, our data provide an improved and finer scale understanding of reproductive ecology and phenology in ungulates.

Spatial prey availability and pulsed reproductive tactics: Encounter risk in a canid-ungulate system.

Predation risk is a function of spatiotemporal overlap between predator and prey, as well as behavioural responses during encounters. Dynamic factors (e.g. group size, prey availability and animal movement or state) affect risk, but rarely are integrated in risk assessments. Our work targets a system where predation risk is fundamentally linked to temporal patterns in prey abundance and behaviour. For neonatal ungulate prey, risk is defined within a short temporal window during which the pulse in parturition, increasing movement capacity with age and antipredation tactics have the potential to mediate risk. In our coyote-mule deer (Canis latrans-Odocoileus hemionus) system, leveraging GPS data collected from both predator and prey, we tested expectations of shared enemy and reproductive risk hypotheses. We asked two questions regarding risk: (A) How does primary and alternative prey habitat, predator and prey activity, and reproductive tactics (e.g. birth synchrony and maternal defence) influence the vulnerability of a neonate encountering a predator? (B) How do the same factors affect behaviour by predators relative to the time before and after an encounter? Despite increased selection for mule deer and intensified search behaviour by coyotes during the peak in mule deer parturition, mule deer were afforded protection from predation via predator swamping, experiencing reduced per-capita encounter risk when most neonates were born. Mule deer occupying rabbit habitat (Sylvilagus spp.; coyote's primary prey) experienced the greatest risk of encounter but the availability of rabbit habitat did not affect predator behaviour during encounters. Encounter risk increased in areas with greater availability of mule deer habitat: coyotes shifted their behaviour relative to deer habitat, and the pulse in mule deer parturition and movement of neonatal deer during encounters elicited increased speed and tortuosity by coyotes. In addition to the spatial distribution of prey, temporal patterns in prey availability and animal behavioural state were fundamental in defining risk. Our work reveals the nuanced consequences of pulsed availability on predation risk for alternative prey, whereby responses by predators to sudden resource availability, the lasting effects of diversionary prey and inherent antipredation tactics ultimately dictate risk.

Experimental SARS-CoV-2 Infection of Elk and Mule Deer.

To assess the susceptibility of elk (Cervus canadensis) and mule deer (Odocoileus hemionus) to SARS-CoV-2, we performed experimental infections in both species. Elk did not shed infectious virus but mounted low-level serologic responses. Mule deer shed and transmitted virus and mounted pronounced serologic responses and thus could play a role in SARS-CoV-2 epidemiology.

Evaluating movement-based methods for estimating the frequency and timing of parturition in mule deer.

BACKGROUND: Information on reproduction of harvested species such as mule deer (Odocoileus hemionus) is vital for conservation and management. Furthermore, parturition in ungulates may be detected using patterns of movement logged by GPS transmitters. Several movement-based methods have been developed to detect parturition in ungulates including the Peterson method, behavioral change point analysis (BCPA), rolling minimum convex polygons (rMCP), individual-based method (IBM), and population-based method (PBM). Our objectives were to (1) test the accuracy and the precision of each previously described method and (2) develop an improved method optimized for mule deer that incorporated aspects of the other methods. METHODS: We determined parturition timing and status for female mule deer fitted with GPS collars and implanted with vaginal implant transmitters (VITs). We used movement patterns before and after parturition to set movement thresholds for each movement-based method. Following model training, we used location and birth date data from an external dataset to test the effectiveness of each movement-based method. Additionally, we developed a novel method for detecting parturition called the analysis of parturition indicators (API). We used two regression analyses to determine the accuracy and precision of estimates generated by each method. RESULTS: The six methods we employed varied in accuracy, with the API, rMCP, and BCPA being most accurate. Precision also varied among methods, with the API, rMCP, and PBM generating the most precise estimates of parturition dates. The API and the rMCP performed similarly and better overall than any of the other existing methods. CONCLUSIONS: We found that movement-based methods could be used to accurately and precisely detect parturition in mule deer. Further, we determined that the API and rMCP methods had the greatest overall success at detecting parturition in mule deer. The relative success of the API and rMCP may be attributed to the fact that both methods use home range size to detect parturition and are validated using known parturition dates of collared deer. We present the API as an efficient method of estimating birth status and timing of parturition of mule deer fitted with GPS transmitters, as well as affirm the effectiveness of a previously developed method, rMCP.

A test of the frost wave hypothesis in a temperate ungulate.

Growing evidence supports the hypothesis that temperate herbivores surf the green wave of emerging plants during spring migration. Despite the importance of autumn migration, few studies have conceptualized resource tracking of temperate herbivores during this critical season. We adapted the frost wave hypothesis (FWH), which posits that animals pace their autumn migration to reduce exposure to snow but increase acquisition of forage. We tested the FWH in a population of mule deer in Wyoming, USA by tracking the autumn migrations of n = 163 mule deer that moved 15-288 km from summer to winter range. Migrating deer experienced similar amounts of snow but 1.4-2.1 times more residual forage than if they had naïve knowledge of when or how fast to migrate. Importantly, deer balanced exposure to snow and forage in a spatial manner. At the fine scale, deer avoided snow near their mountainous summer ranges and became more risk prone to snow near winter range. Aligning with their higher tolerance of snow and lingering behavior to acquire residual forage, deer increased stopover use by 1 ± 1 day (95% CI) day for every 10% of their migration completed. Our findings support the prediction that mule deer pace their autumn migration with the onset of snow and residual forage, but refine the FWH to include movement behavior en route that is spatially dynamic.

Industrial development alters wolf spatial distribution mediated by prey availability.

Increasing resource extraction and human activity are reshaping species' spatial distributions in human-altered landscape and consequently shaping the dynamics of interspecific interactions, such as between predators and prey. To evaluate the effects of industrial features and human activity on the occurrence of wolves (Canis lupus), we used wildlife detection data collected in 2014 from an array of 122 remote wildlife camera traps in Alberta's Rocky Mountains and foothills near Hinton, Canada. Using generalized linear models, we compared the occurrence frequency of wolves at camera sites to natural land cover, industrial disturbance (forestry and oil/gas exploration), human activity (motorized and non-motorized), and prey availability (moose, Alces alces; elk, Cervus elaphus; mule deer, Odocoileus hemionus; and white-tailed deer, Odocoileus virginianus). Industrial block features (well sites and cutblocks) and prey (elk or mule deer) availability interacted to influence wolf occurrence, but models including motorized and non-motorized human activity were not strongly supported. Wolves occurred infrequently at sites with high densities of well sites and cutblocks, except when elk or mule deer were frequently detected. Our results suggest that wolves risk using industrial block features when prey occur frequently to increase predation opportunities, but otherwise avoid them due to risk of human encounters. Effective management of wolves in anthropogenically altered landscapes thus requires the simultaneous consideration of industrial block features and populations of elk and mule deer.

Hidden Markov movement models reveal diverse seasonal movement patterns in two North American ungulates.

Animal movement is the mechanism connecting landscapes to fitness, and understanding variation in seasonal animal movements has benefited from the analysis and categorization of animal displacement. However, seasonal movement patterns can defy classification when movements are highly variable. Hidden Markov movement models (HMMs) are a class of latent-state models well-suited to modeling movement data. Here, we used HMMs to assess seasonal patterns of variation in the movement of pronghorn (Antilocapra americana), a species known for variable seasonal movements that challenge analytical approaches, while using a population of mule deer (Odocoileus hemionus), for whom seasonal movements are well-documented, as a comparison. We used population-level HMMs in a Bayesian framework to estimate a seasonal trend in the daily probability of transitioning between a short-distance local movement state and a long-distance movement state. The estimated seasonal patterns of movements in mule deer closely aligned with prior work based on indices of animal displacement: a short period of long-distance movements in the fall season and again in the spring, consistent with migrations to and from seasonal ranges. We found seasonal movement patterns for pronghorn were more variable, as a period of long-distance movements in the fall was followed by a winter period in which pronghorn were much more likely to further initiate and remain in a long-distance movement pattern compared with the movement patterns of mule deer. Overall, pronghorn were simply more likely to be in a long-distance movement pattern throughout the year. Hidden Markov movement models provide inference on seasonal movements similar to other methods, while providing a robust framework to understand movement patterns on shorter timescales and for more challenging movement patterns. Hidden Markov movement models can allow a rigorous assessment of the drivers of changes in movement patterns such as extreme weather events and land development, important for management and conservation.

Mule deer fawn recruitment dynamics in an energy disturbed landscape.

Wildlife population dynamics are modulated by abiotic and biotic factors, typically climate, resource availability, density-dependent effects, and predator-prey interactions. Understanding whether and how human-caused disturbances shape these ecological processes is helpful for the conservation and management of wildlife and their habitats within increasingly human-dominated landscapes. However, many jurisdictions lack either long-term longitudinal data on wildlife populations or measures of the interplay between human-mediated disturbance, climate, and predator density. Here, we use a 50-year time series (1962-2012) on mule deer (Odocoileus hemionus) demographics, seasonal weather, predator density, and oil and gas development patterns from the North Dakota Badlands, USA, to investigate long-term effects of landscape-level disturbance on mule deer fawn fall recruitment, which has declined precipitously over the last number of decades. Mule deer fawn fall recruitment in this study represents the number of fawns per female (fawn:female ratio) that survive through the summer to October. We used this fawn recruitment index to evaluate the composite effects of interannual extreme weather conditions, energy development, and predator density. We found that density-dependent effects and harsh seasonal weather were the main drivers of fawn fall recruitment in the North Dakota Badlands. These effects were further shaped by the interaction between harsh seasonal weather and predator density (i.e., lower fawn fall recruitment when harsh weather was combined with higher predator density). Additionally, we found that fawn fall recruitment was modulated by interactions between seasonal weather and energy development (i.e., lower fawn fall recruitment when harsh weather was combined with higher density of active oil and gas wells). Interestingly, we found that the combined effect of predator density and energy development was not interactive but rather additive. Our analysis demonstrates how energy development may modulate fluctuations in mule deer fawn fall recruitment concurrent with biotic (density-dependency, habitat, predation, woody vegetation encroachment) and abiotic (harsh seasonal weather) drivers. Density-dependent patterns emerge, presumably due to limited quality habitat, being the primary factor influencing fall fawn recruitment in mule deer. Secondarily, stochastic weather events periodically cause dramatic declines in recruitment. And finally, the additive effects of human disturbance and predation can induce fluctuations in fawn fall recruitment. Here we make the case for using long-term datasets for setting long-term wildlife management goals that decision makers and the public can understand and support.

Fencing amplifies individual differences in movement with implications on survival for two migratory ungulates.

Fences have recently been recognized as one of the most prominent linear infrastructures on earth. As animals traverse fenced landscapes, they adjust movement behaviours to optimize resource access while minimizing energetic costs of coping with fences. Examining individual responses is key for connecting localized fence effects with population dynamics. We investigated the multi-scale effects of fencing on animal movements, space use and survival of 61 pronghorn and 96 mule deer on a gradient of fence density in Wyoming, USA. Taking advantage of the recently developed Barrier Behaviour Analysis, we classified individual movement responses upon encountering fences (i.e. barrier behaviours). We adopted the reaction norm framework to jointly quantify individual plasticity and behavioural types of barrier behaviours, as well as behaviour syndromes between barrier behaviours and animal space use. We also assessed whether barrier behaviours affect individual survival. Our results highlighted a high-level individual plasticity encompassing differences in the degree and direction of barrier behaviours for both pronghorn and mule deer. Additionally, these individual differences were greater at higher fence densities. For mule deer, fence density determined the correlation between barrier behaviours and space use and was negatively associated with individual survival. However, these relationships were not statistically significant for pronghorn. By integrating approaches from movement ecology and behavioural ecology with the emerging field of fence ecology, this study provides new evidence that an extraordinarily widespread linear infrastructure uniquely impacts animals at the individual level. Managing landscape for lower fence densities may help prevent irreversible behavioural shifts for wide-ranging animals in fenced landscapes.

Standardizing protocols for determining the cause of mortality in wildlife studies.

Mortality site investigations of telemetered wildlife are important for cause-specific survival analyses and understanding underlying causes of observed population dynamics. Yet, eroding ecoliteracy and a lack of quality control in data collection can lead researchers to make incorrect conclusions, which may negatively impact management decisions for wildlife populations. We reviewed a random sample of 50 peer-reviewed studies published between 2000 and 2019 on survival and cause-specific mortality of ungulates monitored with telemetry devices. This concise review revealed extensive variation in reporting of field procedures, with many studies omitting critical information for the cause of mortality inference. Field protocols used to investigate mortality sites and ascertain the cause of mortality are often minimally described and frequently fail to address how investigators dealt with uncertainty. We outline a step-by-step procedure for mortality site investigations of telemetered ungulates, including evidence that should be documented in the field. Specifically, we highlight data that can be useful to differentiate predation from scavenging and more conclusively identify the predator species that killed the ungulate. We also outline how uncertainty in identifying the cause of mortality could be acknowledged and reported. We demonstrate the importance of rigorous protocols and prompt site investigations using data from our 5-year study on survival and cause-specific mortality of telemetered mule deer (Odocoileus hemionus) in northern California. Over the course of our study, we visited mortality sites of neonates (n = 91) and adults (n = 23) to ascertain the cause of mortality. Rapid site visitations significantly improved the successful identification of the cause of mortality and confidence levels for neonates. We discuss the need for rigorous and standardized protocols that include measures of confidence for mortality site investigations. We invite reviewers and journal editors to encourage authors to provide supportive information associated with the identification of causes of mortality, including uncertainty.

Responses to natural gas development differ by season for two migratory ungulates.

While migrating, animals make directionally persistent movements and may only respond to human-induced rapid environmental change (HIREC), such as climate and land-use change, once a threshold of HIREC is surpassed. In contrast, animals on other seasonal ranges (e.g., winter range) make more localized and tortuous movements while foraging and may have the flexibility to adjust the location of their range and the intensity of use within it to minimize interactions with HIREC. Because of these seasonal differences in movement, animals on seasonal ranges should avoid areas that contain any level of HIREC, however, during migration, animals should use areas that contain low levels of HIREC, avoiding it only once a threshold of HIREC has been surpassed. We tested this hypothesis using a decade of GPS collar data collected from migratory mule deer (Odocoileus hemionus; n = 56 migration, 143 winter) and pronghorn (Antilocapra americana; n = 70 migration, 89 winter) that winter on and migrate through a natural gas field in western Wyoming. Using surface disturbance caused by well pads and roads as an index of HIREC, we evaluated behavioral responses across three spatial scales during winter and migration seasons. During migration, both species tolerated low levels of disturbance. Once a disturbance threshold was surpassed, however, they avoided HIREC. For mule deer, thresholds were consistently ~3%, whereas thresholds for pronghorn ranged from 1% to 9.25% surface disturbance. In contrast to migration, both species generally avoided all levels of HIREC while on winter range. Our study suggests that animal responses to HIREC are mediated by season-specific movement patterns. Our results provide further evidence of ungulates avoiding human disturbance on winter range and reveal disturbance thresholds that trigger mule deer and pronghorn responses during migration: information that managers can use to maintain the ecological function of migration routes and winter ranges.

Risky business: How an herbivore navigates spatiotemporal aspects of risk from competitors and predators.

Understanding factors that influence animal behavior is central to ecology. Basic principles of animal ecology imply that individuals should seek to maximize survival and reproduction, which means carefully weighing risk against reward. Decisions become increasingly complex and constrained, however, when risk is spatiotemporally variable. We advance a growing body of work in predator-prey behavior by evaluating novel questions where a prey species is confronted with multiple predators and a potential competitor. We tested how fine-scale behavior of female mule deer (Odocoileus hemionus) during the reproductive season shifted depending upon spatial and temporal variation in risk from predators and a potential competitor. We expected female deer to avoid areas of high risk when movement activity of predators and a competitor were high. We used GPS data collected from 76 adult female mule deer, 35 adult female elk, 33 adult coyotes, and six adult mountain lions. Counter to our expectations, female deer exhibited selection for multiple risk factors, however, selection for risk was dampened by the exposure to risk within home ranges of female deer, producing a functional response in habitat selection. Furthermore, temporal variation in movement activity of predators and elk across the diel cycle did not result in a shift in movement activity by female deer. Instead, the average level of risk within their home range was the predominant factor modulating the response to risk by female deer. Our results counter prevailing hypotheses of how large herbivores navigate risky landscapes and emphasize the importance of accounting for the local environment when identifying effects of risk on animal behavior. Moreover, our findings highlight additional behavioral mechanisms used by large herbivores to mitigate multiple sources of predation and potential competitive interactions.

Investigation into Causes of Antler Deformities in Mule Deer (Odocoileus hemionus) Bucks in Southern Utah, USA.

We investigated causes of antler deformities in mule deer (Odocoileus hemionus) bucks from the Paunsaugunt Plateau in southern Utah, US. A total of 10 hunter-harvested and nine live-captured bucks with antler deformities and six hunter-harvested and 43 live-captured bucks with normal antlers were included in the study. All were screened by serology for exposure to epizootic hemorrhagic disease virus (EHDV) types 1, 2, and 6, bluetongue virus, Brucella abortus, Brucella ovis, caprine arthritis, encephalitis virus, and bovine viral diarrhea virus (BVDV) types 1 and 2. Serum testosterone and trace minerals concentrations were measured, and whole-blood counts evaluated. Testicular tissue from the hunter-harvested bucks was tested by quantitative PCR for EHDV-1, -2, and -6. All bucks with antler deformities had low to non-detectable serum testosterone concentrations from end-stage fibrosing orchitis, and EHDV-2 was detected by quantitative PCR in one of the testicular tissues tested. All bucks with antler deformities were seropositive for EHDV-2. In comparison, only 53% of bucks with normal antlers were EHDV-2 seropositive. More than 67% of hunter-harvested and live-captured bucks had antibodies to BVDV-1 and BVDV-2, probably because of high cattle exposure in the area. Our results support previous research linking infection with EHDV-2 to development of antler deformities; however, it remains unclear why some infected bucks develop testicular fibrosis, and others recover from the infection.

MOUNTAIN LIONS (PUMA CONCOLOR) RESIST LONG-TERM DIETARY EXPOSURE TO CHRONIC WASTING DISEASE.

For nearly 18 yr, we evaluated susceptibility of captive mountain lions (Puma concolor) to chronic wasting disease (CWD) in the face of repeated exposure associated with consuming infected cervid carcasses. Three mountain lions with a monomorphic prion protein gene (PRNP) sequence identical to that described previously for the species had access to parts of ≥432 infected carcasses during ≥2,013 feeding occasions, conservatively representing >14,000 kg of infected feed material, during May 2002 to March 2020. The proportion of diet in infected carcass material averaged 43% overall but differed from year to year (minimally 11-74%). Most infected carcasses were mule deer (Odocoileus hemionus; ∼75%). We observed no clinical signs suggestive of progressive encephalopathy or other neurologic disease over the ∼14.5-17.9 yr between first known exposure and eventual death. Histopathology revealed no spongiform changes or immunostaining suggestive of prion infection in multiple sections of nervous and lymphoid tissue. Similarly, none of 133 free-ranging mountain lion carcasses sampled opportunistically during 2004-20 showed immunostaining consistent with prion infection in sections of brainstem or lymph node. These findings align with prior work suggesting that CWD-associated prions face strong barriers to natural transmission among species outside the family Cervidae.

THE RELATIONSHIP BETWEEN HARVEST MANAGEMENT AND CHRONIC WASTING DISEASE PREVALENCE TRENDS IN WESTERN MULE DEER (ODOCOILEUS HEMIONUS) HERDS.

We analyzed retrospective data on harvest management practices and corresponding chronic wasting disease (CWD) prevalence trends in 36 western US and Canadian mule deer (Odocoileus hemionus) management units (units). Our analyses employed logistic regression and model selection, exploiting variation in practices within and among jurisdictions to examine relationships between harvest management and apparent prevalence (the proportion of positive animals among those sampled). Despite notable differences in hunting practices among jurisdictions, our meta-analysis of combined data revealed strong evidence that the amount of harvest was related to CWD prevalence trends among adult male mule deer in the 32 units where prevalence at the start of the analysis period was ≤5%. All competitive models included the number of male deer harvested or number of hunters 1-2 yr prior as an explanatory variable, with increasing harvest leading to lower prevalence among males harvested in the following year. Competitive models also included harvest timing. Although less definitive than the number harvested, median harvest dates falling closer to breeding seasons were associated with lower prevalence in the following year. Our findings suggest harvest-when sufficient and sustained-can be an effective tool for attenuating CWD prevalence in adult male mule deer across western ranges, especially early in the course of an epidemic. Evidence of a broad relationship between the amount of harvest and subsequent changes in CWD prevalence among adult male mule deer provides an empirical basis for undertaking adaptive disease management experimentation aimed at suppressing or curtailing CWD epidemics.

Landscape Seroprevalence of Three Hemorrhagic Disease-Causing Viruses in a Wild Cervid.

Disease plays a major role in shaping wildlife populations worldwide, and changes in landscape conditions can significantly influence risk of pathogen exposure, a threat to vulnerable wild species. Three viruses that cause hemorrhagic disease affect cervid populations in the USA (Odocoileus hemionus adenovirus, bluetongue virus, and epizootic hemorrhagic disease virus), but little is known of their distribution and prevalence in wild populations. We explored the distribution and co-occurrence of seroprevalence of these three pathogens in southern mule deer (Odocoileus hemionus fuliginatus), a subspecies of conservation concern and a harvested species native to southern California, to evaluate the distribution of exposure to these pathogens relative to landscape attributes. We found that habitat type, level of development, and proximity to livestock may affect hemorrhagic disease seroprevalence in southern mule deer. Continued monitoring of hemorrhagic disease-causing viruses in areas where deer are in proximity to cattle and human development is needed to better understand the implications of future outbreaks in wild populations and to identify opportunities to mitigate disease impacts in southern mule deer and other cervid species.

Recombination Marks the Evolutionary Dynamics of a Recently Endogenized Retrovirus.

All vertebrate genomes have been colonized by retroviruses along their evolutionary trajectory. Although endogenous retroviruses (ERVs) can contribute important physiological functions to contemporary hosts, such benefits are attributed to long-term coevolution of ERV and host because germline infections are rare and expansion is slow, and because the host effectively silences them. The genomes of several outbred species including mule deer (Odocoileus hemionus) are currently being colonized by ERVs, which provides an opportunity to study ERV dynamics at a time when few are fixed. We previously established the locus-specific distribution of cervid ERV (CrERV) in populations of mule deer. In this study, we determine the molecular evolutionary processes acting on CrERV at each locus in the context of phylogenetic origin, genome location, and population prevalence. A mule deer genome was de novo assembled from short- and long-insert mate pair reads and CrERV sequence generated at each locus. We report that CrERV composition and diversity have recently measurably increased by horizontal acquisition of a new retrovirus lineage. This new lineage has further expanded CrERV burden and CrERV genomic diversity by activating and recombining with existing CrERV. Resulting interlineage recombinants then endogenize and subsequently expand. CrERV loci are significantly closer to genes than expected if integration were random and gene proximity might explain the recent expansion of one recombinant CrERV lineage. Thus, in mule deer, retroviral colonization is a dynamic period in the molecular evolution of CrERV that also provides a burst of genomic diversity to the host population.

Variation in movement patterns of mule deer: have we oversimplified migration?

BACKGROUND: Conservation and management of migratory animals has gained attention in recent years, but the majority of research has focused on stereotypical 'migrant' and 'resident' behaviors, often failing to incorporate any atypical behaviors or characterize migratory behaviors beyond distance and timing of the migration. With migration threatened by anthropogenic development and climate change, it is crucial that we understand the full range of migratory behaviors. Our objective was to demonstrate and characterize the variation in migration strategies, including typical and atypical migratory behaviors for mule deer (Odocoileus hemionus) in Utah, USA. METHODS: Because calculation of common metrics such as distance, timing, and use of stopovers during migration did not adequately describe the variation we observed in migratory behavior for this species-particularly when animals visited multiple (> 3) ranges for extended lengths of time-we developed additional methods and categories to describe observed variation in migratory behavior. We first categorized trajectories based on the number of discrete, separate ranges and range shifts between them. Then, we further characterized the variation in migration strategies by examining the timing, duration, and distance traveled within each of the categories. We also examined if and how frequently individual deer switched among categories from year to year. RESULTS: We classified 1218 movement trajectories from 722 adult female mule deer, and found that 54.4% were dual-range migrants, who made one round-trip to one distinct range. Multi-range migrants (23.6%) made one round-trip during which they stayed at multiple discrete ranges. Commuters (1.0%) traveled to the same range multiple times, and poly migrants (1.5%) made multiple round-trips to different ranges. Gradual movers (2.5%) did not show a discrete range shift but moved gradually between ranges, whereas residents (12.6%) never left their home ranges, and dispersers (4.4%) left but never returned. Of the deer that we monitored for multiple years, 51.2% switched among categories. CONCLUSION: We conclude that the substantial number of atypical migratory strategies, as well as the number of deer that switched categories, underlines the importance of studying these less-stereotyped behaviors that may be exhibited by large proportions of populations. Acknowledging and investigating the full complexity and diversity in migratory strategies might uncover unknowns with respect to underlying factors and drivers of migration, and can help shape effective conservation strategies.