https://cpw.cvlcollections.org/files/original/25699f5160237a87ec92d3c18c57421c.pdf 0352c334f3e492539c407151b6e6d9b4 PDF Text Text The research in this publication was partially or fully funded by Colorado Parks and Wildlife. Dan Prenzlow, Director, Colorado Parks and Wildlife • Parks and Wildlife Commission: Marvin McDaniel, Chair • Carrie Besnette Hauser, Vice-Chair Marie Haskett, Secretary • Taishya Adams • Betsy Blecha • Charles Garcia • Dallas May • Duke Phillips, IV • Luke B. Schafer • James Jay Tutchton • Eden Vardy �Migrating Mule Deer: Effects of Anthropogenically Altered Landscapes Patrick E. Lendrum1*, Charles R. Anderson Jr.2, Kevin L. Monteith1,3, Jonathan A. Jenks4, R. Terry Bowyer1 1 Department of Biological Sciences, Idaho State University, Pocatello, Idaho, United States of America, 2 Colorado Division of Parks and Wildlife, Grand Junction, Colorado, United States of America, 3 Wyoming Cooperative Fish and Wildlife Research Unit, University of Wyoming, Laramie, Wyoming, United States of America, 4 Department of Natural Resource Management, South Dakota State University, Brookings, South Dakota, United States of America Abstract Background: Migration is an adaptive strategy that enables animals to enhance resource availability and reduce risk of predation at a broad geographic scale. Ungulate migrations generally occur along traditional routes, many of which have been disrupted by anthropogenic disturbances. Spring migration in ungulates is of particular importance for conservation planning, because it is closely coupled with timing of parturition. The degree to which oil and gas development affects migratory patterns, and whether ungulate migration is sufficiently plastic to compensate for such changes, warrants additional study to better understand this critical conservation issue. Methodology/Principal Findings: We studied timing and synchrony of departure from winter range and arrival to summer range of female mule deer (Odocoileus hemionus) in northwestern Colorado, USA, which has one of the largest natural-gas reserves currently under development in North America. We hypothesized that in addition to local weather, plant phenology, and individual life-history characteristics, patterns of spring migration would be modified by disturbances associated with natural-gas extraction. We captured 205 adult female mule deer, equipped them with GPS collars, and observed patterns of spring migration during 2008–2010. Conclusions/Significance: Timing of spring migration was related to winter weather (particularly snow depth) and access to emerging vegetation, which varied among years, but was highly synchronous across study areas within years. Additionally, timing of migration was influenced by the collective effects of anthropogenic disturbance, rate of travel, distance traveled, and body condition of adult females. Rates of travel were more rapid over shorter migration distances in areas of high natural-gas development resulting in the delayed departure, but early arrival for females migrating in areas with high development compared with less-developed areas. Such shifts in behavior could have consequences for timing of arrival on birthing areas, especially where mule deer migrate over longer distances or for greater durations. Citation: Lendrum PE, Anderson CR Jr, Monteith KL, Jenks JA, Bowyer RT (2013) Migrating Mule Deer: Effects of Anthropogenically Altered Landscapes. PLoS ONE 8(5): e64548. doi:10.1371/journal.pone.0064548 Editor: Ofer Ovadia, Ben-Gurion University of the Negev, Israel Received January 13, 2013; Accepted April 17, 2013; Published May 14, 2013 Copyright: ß 2013 Lendrum et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: This project was funded and supported by the Colorado Parks and Wildlife (CPW). The authors thank C. Bishop, D. Freddy, and M. Michaels from CPW for helping administer the project. Additionally, the authors thank personnel at Little Hills State Wildlife Area for field support and coordination. The authors thank Quicksilver Air Inc. for assistance in capturing deer from helicopters, L. Gepfert and L. Coulter for fixed-wing aircraft support, and L. Wolfe, C. Bishop, and D. Finley of CPW for crucial assistance during capture efforts. Additional funding and support came from Federal Aid in Wildlife Restoration, Colorado Mule Deer Association, Colorado Mule Deer Foundation, Colorado Oil and Gas Conservation Commission, Williams Production LMT CO., EnCana Corp., Exxon Mobil Production Co., Shell Petroleum, Marathon Oil Corp., and Idaho State University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors acknowledge funding from multiple commercial sources (Williams Production LMT CO., EnCana Corp., Exxon Mobil Production Co., Shell Petroleum and Marathon Oil Corp) and declare that this does not alter the authors’ adherence to all the PLOS ONE policies on sharing data and materials. * E-mail: lendpatr@isu.edu migration is sufficiently plastic to compensate for such change is uncertain. Understanding effects of those disturbances on the ability of migratory ungulates to follow phenological gradients and thereby maximize energy intake, especially when parturition is looming, is of particular importance for conservation planning [7]. Indeed, the demise of migrating populations of large mammals is of increasing concern as extraction of non-renewable resources proceeds at unprecedented rates [4,6]. Migration is an adaptive strategy that is thought to allow animals to minimize resource shortages and reduce risk of predation [8] at broad geographic scales [9], both of which affect fitness [10]. In temperate and arctic environments, heterogeneous Introduction Ungulates generally migrate along traditional routes [1,2], and demonstrate high fidelity to seasonal ranges [3]. Recently, however, many of those migratory routes have been disrupted by increasing levels of anthropogenic disturbance [4,5]. Threats to remaining long-distance migration of ungulates include energy development, tourism, urban sprawl, highway mortality, and habitat fragmentation [6]. If traditional migration routes are blocked or impeded, individuals may not be able to modify their migratory behavior, which could compromise persistence of those populations [2]. How overlap of migration routes with oil and gas development will affect migratory patterns, and whether ungulate PLOS ONE | www.plosone.org 1 May 2013 | Volume 8 | Issue 5 | e64548 �Mule Deer Migration Patterns ing areas with higher levels of natural-gas extraction, compared with areas of less development. forages may be difficult to exploit by remaining sedentary, because forage quantity and quality are seasonally and spatially dynamic [3]. Migration thereby provides nutritional advantages to individuals that exploit seasonal peaks in resources at different locations, while avoiding inclement weather and reducing intraspecific competition [11,12]. Across western North America, mule deer (Odocoileus hemionus), North American elk (Cervus elaphus), pronghorn (Antilocapra americana), moose (Alces alces), and bighorn sheep (Ovis canadensis) may migrate 50–260 km between seasonal ranges [2]. Ungulate migrations in temperate environments are typified by movements between high elevations in summer and low elevations in winter [3,12,13]. Winter ranges are typically smaller and populations occur at relatively high density; therefore, movement to spring and summer ranges provides a release from a restricted food supply during a time when costs of reproduction are rising [8,12,14]. We evaluated factors that influenced patterns of spring migration for mule deer in the Piceance Basin in northwestern Colorado, USA, which included areas that were exposed to a wide range of natural-gas development. The deer in this study were the same population studied to examine resource selection along migration routes in Lendrum et al. [15]. Here we strive to determine how extrinsic and intrinsic factors influenced the timing of migration – we used different inferential techniques and a different focus, as well as a separate suite of hypotheses. Based on previous research [3,12], we expected timing of spring migration to be influenced by extrinsic factors associated with local weather and plant phenology. Onset of spring migration often is initiated by rising temperatures, decreasing snow cover, and increasing plant growth [3,12]. Accordingly, we tested whether spring migrations would be initiated by some combination of increased solar radiation and temperature, decreased snow depth, and advancing plant phenology. In addition, we predicted that deer residing at high elevations during summer, where snow depths were likely greater and green-up delayed [11,16], would either postpone initiation of, or exhibit a slower rate of migration than deer inhabiting lower-elevation ranges during summer. Although patterns of migration may be synchronized by environmental factors, migratory decisions can vary among individuals depending upon life-history characteristics and nutritional condition [12,17]. For example, Monteith et al. [12] observed that old female mule deer in the Sierra Nevada risked encountering severe weather by delaying autumn migration. Assessment of the relative vulnerability of migratory populations requires careful consideration of both extrinsic and intrinsic factors. Consequently, we predicted that older, more-experienced deer, and those in comparatively good condition would initiate spring migration earlier compared with young, inexperienced deer, and those in relatively poor condition [3,12]. Disturbances have the potential to override the contribution of other factors to ungulate migrations [5,7]. After accounting for effects of environmental and individual-based factors on patterns of migration, we further hypothesized that migration would be modified by the presence of natural-gas development. Lendrum et al. [15] documented longer step lengths by mule deer migrating through areas of greater development; however they did not examine what effects this pattern of movement may have had on overall timing of migration. Similarly, Sawyer et al. [18] observed that deer migrated faster through areas with a higher density of development compared with when the same deer were migrating through less-developed areas. Because mule deer sometimes avoid oil and gas developments [19,20], we predicted an early departure, faster movement rate, and early arrival for individuals experienc- PLOS ONE | www.plosone.org Methods Ethics Statement All aspects of animal handling and research complied with the methods adopted by the American Society of Mammalogists [21], and were approved by an Animal Care and Use Committee at Idaho State University (protocol # 670 0410). Permission to conduct all aspects of this field study was provided, in a collaborative effort, by Colorado Division of Parks and Wildlife and the U.S. Bureau of Land Management. Study Area The Piceance Basin supports one of the largest populations of migratory mule deer in North America, historically estimated at 21,000–27,000 animals [22]. This area also includes one of the largest natural-gas reserves in North America. Energy development throughout northwestern Colorado is projected to increase from approximately 500 to 17,000 wells over the next 20 years (U.S. Bureau of Land Management Executive Summary 2007). The Piceance Basin is topographically diverse and characterized by pinyon pine (Pinus edulis)-Utah juniper (Juniperus osteosperma) shrubland. Climate was typified by warm, dry summers (28uC high average) and cold winters (212uC low average), with most annual moisture coming from snow melt in spring (Western Regional Climate Center, 1893–2010). This area was dissected by numerous drainages with stands of big sagebrush (Artemisia tridentate), saltbrush (Atriplex spp.), black greasewood (Sarcobatus vermiculatus), and rabbitbrush (Crysothamnus spp.), with most of the primary drainage bottoms converted to fields of mixed-grass hay. Primary winter habitat for mule deer ranged from 1,675 to 2,285 m in elevation. Summer range for mule deer occurred at high elevations (2,100 to 2,700 m) with dominant vegetation communities of quaking aspen (Populus tremuloides)-Douglas-fir (Pseudotsuga menziesii) and Engelmann spruce (Picea engelmannii)subalpine fir (Abies lasiocarpa) forests [3]. The area contained additional large herbivores including North American elk (Cervus elaphus), wild horses (Equus caballus), and moose (Alces alces). Common species of predators included coyotes (Canis latrans), mountain lions (Puma concolor), bobcats (Lynx rufus), and black bears (Ursus americana). Lendrum et al. [15] provides a more complete description of the study area. We monitored four populations of mule deer that wintered in the Piceance Basin: 1) North Ridge (53 km2) in the northeastern portion of the Basin; 2) Ryan Gulch (141 km2) in the southwestern portion of the Basin; and 3) North Magnolia (79 km2); and 4) South Magnolia (83 km2) in the central portion of the Basin (Figure 1). During spring migration, mule deer from North Ridge and North Magnolia moved easterly to higher elevations across US Highway 13 towards the Flat Top Mountain Range. Mule deer from Ryan Gulch and South Magnolia migrated southerly through a fragmented landscape of well pads, compressor stations, pipelines, and roads to higher elevations along the Roan Plateau (Figure 1). Habitat characteristics (i.e, vegetation type) occurring along the migratory paths of mule deer was similar among study areas [15], which allowed us to control for such variation while examining what other variables may have influenced patterns of migration. Animal Capture We captured mule deer from a helicopter using net guns [23] to acquire a sample of adult (.1 years old) females (n = 205 deer; 2 May 2013 | Volume 8 | Issue 5 | e64548 �Mule Deer Migration Patterns Figure 1. Map of study area in the Piceance Basin, Colorado, USA. Seasonal ranges and general spring migration paths (arrows) of adult female mule deer for North Ridge (white with hash marks), North Magnolia (grey with hash marks), South Magnolia (white), and Ryan Gulch (grey). Active well pads ( ), Western Regional Climate Center (m) and SNOTEL weather station (¤). doi:10.1371/journal.pone.0064548.g001 N For each adult female captured during 2009 and 2010 (n = 160), we measured maximum thickness (cm) of subcutaneous fat deposition at the thickest point cranial to the cranial process of the tuber ischium with a portable ultrasound device (SONOVET 2000, Canmedical, Yanker, Ontario) and a 5-MHz transducer [25]. In addition, we used palpation to determine a bodycondition score, validated for mule deer, to provide an estimate of ingesta-free body fat (IFBFat), when subcutaneous fat reserves had been mobilized [26]. We used a combination of subcutaneous rump-fat thickness and body-condition score to estimate percent IFBFat [26]. We estimated deer age using tooth replacement and wear [27]. North Ridge = 60, North Magnolia = 43, South Magnolia = 42, and Ryan Gulch = 60) from 2008 to 2010. During 10–12 January 2008, we captured 45 adult female mule deer and fit them with Global Positioning Satellite (GPS) collars (14 with GPS-4400S; Lotek Wireless, Newmarket, Ontario, Canada; 31 with G2110D; Advanced Telemetry Systems, Isanti, Minnesota, USA). During late February – early March 2009 and 2010, we captured and radiocollared an additional 60 and 100 adult females, respectively (G2110D GPS collar). We programmed collars to obtain a locational fix every 5 h during the migration period, and only retained 3D fixes or 2D fixes with a positional dilution of precision ,10 [24]. All collars included mortality sensors and timed drop-off mechanisms that were scheduled to release during April of the year following deployment. We retrieved GPS collars, which also were equipped with Very High Frequency (VHF) transmitters, from the field as they became available following mortalities or collar drop each spring. Of the 205 adult females collared, we censored 16 animals from analyses because they did not provide complete data on migrations. PLOS ONE | www.plosone.org Characteristics of Spring Migration We used Hawth’s Analysis Tools in ArcGIS 9.3 (ESRI, Redlands, California, USA) to derive 95% kernel-density estimates of seasonal ranges for each individual. We determined the initiation of spring migration based on the day a particular deer left winter range on a trajectory path (i.e., three successive 3 May 2013 | Volume 8 | Issue 5 | e64548 �Mule Deer Migration Patterns We again used PCA, based on the correlation matrix, to derive independent composites that described daily weather [12]. The PCA included 14 variables representing absolute daily weather and a metric of change in daily weather for minimum, maximum, and average temperature (uC), relative humidity (%), solar radiation (Watts/m2), precipitation (cm), and snow depth (cm). If weather data were unavailable for a particular day (,1% were missing), we replaced the missing value with the average value of the previous and subsequent day. We selected five principal components that were both biologically relevant and explained .5% of the variation in daily weather [12]. Principle component 1 (WPC1) explained 47.3% of the variation in daily weather, and represented daily changes in temperature with an influence of humidity from cooling temperatures and moister days (negative loadings) to warming temperatures and dryer days (positive loadings). Principle component 2 (WPC2) explained 17.7% of the variation and reflected an absolute measure of daily snow depth with an influence of daily temperature from lower snow depths and warmer days (negative loadings) to higher snow depths and cooler days (positive loadings). Principle component 3 (WPC3) explained 9.9% of the variation and reflected daily changes in precipitation from wetter (negative loadings) to dryer (positive loadings) days. Principle component 4 (WPC4) explained 8.0% of the variation and was largely related to solar radiation from overcast (negative loadings) to sunny (positive loadings) days. Principle component 5 (WPC5) explained 6.5% of the variation and was almost entirely related to changes in snow depth from decreasing (negative loadings) to increasing (positive loadings) depths. locations leading away from winter range); arrival on summer range was determined as the first location inside the summer home range for that same deer [3]. We then calculated the distance and rate of travel (distance/days to complete migration) between winter and summer range along the migratory path using the Distance Between Points Tool in Hawth’s Analysis Tools. We also calculated elevation of summer range for each deer as the average elevation of all locations within their summer range. Levels of natural-gas development varied markedly among study areas within the Piceance Basin. Consequently, we calculated well-pad densities along spring migration routes for each radio-collared female. To define the width of migration corridors, we calculated the tortuosity of each movement path as log(N)/(log(N) + log(D/L)), where N is the number of line segments that make up the line, D is the distance between the start and end points of the line, and L is the cumulative length of all line segments. Tortuosity is a useful metric in linking fine-scale movements to habitat quality and human development [28]. We then used the number of well-pads within each buffered migratory path as a measure of well-pad density for individual deer. We tested for differences in characteristics among study areas and years using analysis of variance (ANOVA) with Bonferroni pairwise comparisons in Minitab 16.1.0 (State College, Pennsylvania, USA, 2010). Prior to interpretation of results, we examined residual plots of each dependent variable to test for compliance with assumptions of ANOVA. We hypothesized that characteristics of a migratory path would influence timing of migration including: elevation at summer range (m); distance migrated (km); well-pad density (wells/km2); and rate of travel to summer range (km/day). Because these intrinsic variables were interrelated, we used principal components analysis (PCA), based on the correlation matrix, to reduce dimensionality of those variables and derive independent composites that described characteristics of migratory routes for individual deer (Appendix S1). Principle component 1 (PC1) explained 43.8% of the variation of the intrinsic variables, and contrasted individuals migrating quickly through high well-pad densities for a short distance (positive loadings), with those traveling longer distances at a slower rate with no well pads (negative loadings). Principle component 2 (PC2) explained 26.8% of the variation in intrinsic variables, and was associated positively with elevation on summer range. Plant Phenology We used the Normalized Difference Vegetation Index (NDVI) to reflect primary productivity and greenness of vegetation [30,31], and thus, potential fluctuations in dietary quality and availability associated with spring green-up [12]. We obtained 7day composites of NDVI from MODIS (moderate-resolution imaging spectroradiometry; ftp://emodisftp.cr.usgs.gov/ eMODIS/CONUS/historical/TERRA/) with a 250-m2 spatial resolution. We extracted NDVI values associated with GPS locations of deer for each weekly composite. Once an individual departed winter range, we used the locations for that individual during the last week present on winter range to estimate phenological patterns for winter range for the remainder of the monitoring interval. Conversely, we used locations from the first week on summer range to extract NDVI values representative of phenological patterns on summer range, prior to arrival on summer range for each individual. We used this approach to obtain NDVI data, because simply following the path of an individual once it departed a seasonal range would inherently result in shifts in NDVI that would be caused by movement and be correlated with timing of migration, rather than allowing the assessment of the response of an individual to natural changes in greenness on a specific seasonal range. Local Weather We obtained data on daily weather including: minimum, maximum, and average temperature (uC); relative humidity (%); solar radiation (Watts/m2); and precipitation (cm) from a weather station located within winter range on the North Ridge study area (Western Regional Climate Center 2008–2010; Figure 1). Additionally, we obtained data on snow depth from a SNOTEL weather station located near the summer range of the North Magnolia subpopulation (2,865 m; Figure 1), which served as an index to snow depth for the entire study area. In addition to absolute measures of daily weather, we calculated a metric of change in weather based on the difference in the daily weather variables relative to the mean of that particular weather variable during the previous 2 weeks, because deer may respond to changes in weather rather than to absolute values [12,29]. We tested for effects of year on average temperature, humidity, solar radiation, and snow depth using multivariate analysis of variance (MANOVA). Following a significant main effect, we used canonical correlations to identify variables responsible for overall significance, and followed this with tests using separate ANOVAs for those variables. PLOS ONE | www.plosone.org Modeling of Migration We implemented time-to-event models in Program MARK to predict patterns of spring migration for mule deer, including departure from winter range and arrival on summer range. These methods were developed originally to estimate survival of marked animals [32], but they are a valuable tool for assessing factors that influence time to a specific event, and allow the incorporation of individual-based covariates [12,33]. We estimated daily probability of not migrating as a function of extrinsic and intrinsic factors using the known-fate option in Program MARK [32], and 4 May 2013 | Volume 8 | Issue 5 | e64548 �Mule Deer Migration Patterns arrival, absolute measures of daily weather also differed among years (Wilks’ l = 0.641, F8,112 = 3.47, p = 0.001). Canonical correlation analysis, however, indicated that overall significance was most influenced by mean snow depth (F2,61 = 5.55, p = 0.006; Appendix S2), which was greatest in 2008. Mean dates of departure were highly synchronous among three of four study areas (all p.0.10), with deer from North Ridge departing earliest (all p,0.02); dates of arrival were similar among study areas (F38,188 = 1.51, p = 0.212; Figure 3). Within study areas, measures of tortuosity varied little among individuals (all p.0.05); therefore, we used the average value for each study area as the buffer width for each migration path, the least tortuous being South Magnolia (ANOVA, F3,180 = 4.68, p = 0.004; North Ridge = 1.0860.060, North Magnolia = 1.0660.052, Ryan Gulch = 1.0760.076, South Magnolia = 1.0360.034). Mean well-pad density along migration routes (pads/km2) was lowest in North Ridge and North Magnolia, becoming progressively higher through Ryan Gulch and South Magnolia (ANOVA, F3,188 = 45.93, p,0.001; North Ridge ,0.0160.111, North Magnolia = 0.0560.014, Ryan Gulch = 0.1860.107, South Magnolia = 0.1960.132). Similarly, rates of movement by females during migration differed among study areas (F3,188 = 6.42, p,0.001), and increased with well-pad density. Movement rates (km/day) were highest for deer migrating through the most developed study site (South Magnolia correspondingly calculated daily probability of migrating as 1 minus the daily probability of not migrating. Beginning on 1 April, we used an 80-day interval to construct encounter histories for migration timing of each individual deer. We designated 1 April as the beginning of the interval that deer were available to migrate during each season, because that date was prior to any individual departing from winter range and the 80-day period encompassed the entire migratory period for mule deer in the Piceance Basin. We used a two-stage process to evaluate extrinsic and intrinsic factors that were related to patterns of migration, (sensu [12,34]). We examined extrinsic and intrinsic explanatory variables in a 2step process because we were interested in evaluating the effects environmental variables would have on patterns of migration at the population level, and then how those patterns might be modified by intrinsic variables, which were more specific to each individual. Our first step included examining all possible combinations of extrinsic variables that we predicted would influence timing of spring migration: daily weather variables and weather-change metrics from the PCA; weekly NDVI values on both winter and summer range; year; and study area. We included year and study area as grouping variables to account for variation that was not specifically addressed by our other environmental variables. To identify extrinsic variables that influenced timing of spring migration, we calculated Akaike’s Information Criterion adjusted for small sample size (AICc, DAICc, and Akaike weight wi; [35]) for each model. We determined model-averaged parameter estimates and unconditional standard errors (SE) for each predictor variable [35], and evaluated the significance of each variable by whether the corresponding 90% CI overlapped 0. We also calculated importance weights as the sum of wi for all models that contained a particular variable, and evaluated the relative importance of each variable based on those weights [12,35]. We considered variables to be biologically influential if they had an importance weight of .0.70. After identifying the extrinsic variables that influenced timing of migration, we evaluated the influence of intrinsic factors and lifehistory characteristics of individual mule deer. We combined the extrinsic factors identified as being influential in the first stage of the analysis with intrinsic covariates, PC1 (well-pad density, rate, and distance) and PC2 (elevation). We modeled all possible combinations of extrinsic variables that were previously determined to be important with the addition of the intrinsic variables. We again used model averaging, 90% CI, and importance weights to evaluate effects of individual characteristics on timing of spring migration [35]. In a final analysis using the same modeling approach, we evaluated the influence of age and nutritional condition of individual females, which were only collected during 2009 and 2010. In each stage of the analysis, we also investigated effects of interaction terms among extrinsic and intrinsic variables that we hypothesized would influence patterns of migration by examining DAICc and confidence intervals. We retained interaction terms in a model set if their inclusion resulted in an improvement of model fit (DAICc .2.0), and their parameter estimates differed from zero [35]. Results Figure 2. Yearly patterns of spring migration by mule deer. Model-averaged estimates of the cumulative proportion migrated (black) and the daily probability of migration (grey) relative to Julian date of departure from winter range (a) and arrival to summer range (b), for spring migration of adult female mule deer, Piceance Basin, Colorado, USA, 2008–2010. doi:10.1371/journal.pone.0064548.g002 During 2008–2010, we recorded 189 spring migrations by adult female mule deer. Mean date of departure (F38,188 = 2.12, p = 0.001) and arrival (F38,188 = 1.71, p = 0.013) differed among years, with the earliest migration occurring in 2009 (Appendix S2, Figure 2). Within 10 days of the average date of departure or PLOS ONE | www.plosone.org 5 May 2013 | Volume 8 | Issue 5 | e64548 �Mule Deer Migration Patterns Figure 3. Departure and arrival dates of spring migration by mule deer. Mean Julian date (695% CI) of departure from winter range and arrival to summer range for spring migration of adult female mule deer, Piceance Basin, Colorado, USA, 2008–2010. Mean values for study areas are of combined years and values above or below means are number of collared deer during each year. doi:10.1371/journal.pone.0064548.g003 = 11.666.38) compared with deer in the least developed areas (North Ridge = 7.664.33, North Magnolia = 7.563.55). Additionally, distance traveled differed among study areas (F3,188 = 6.30, p,0.001), with deer from North Ridge and North Magnolia traveling farther than deer from Ryan Gulch and South Magnolia (North Ridge = 53631.8 km, North Magnolia = 47620.7 km, Ryan Gulch = 38617.0 km, South Magnolia = 36610.2 km). Elevations of migratory paths varied among study areas (ANOVA, F3,188 = 16.57, p,0.001; North Ridge = 2,316622.6 m, North Magnolia = 2,417631.2 m, Ryan Gulch = 2,468612.9 m, South Magnolia = 2,497615.2 m), but remained relatively consistent within study areas among years (F2,188 = 2.81, p = 0.063; mean = 2,420 m, SE = 11.5). Median estimated age of female deer was 3.5 and ranged from 1.5 to .10.5 years-of-age. Age of females migrating was similar among years (ANOVA, F1,149 = 0.02, p = 0.89) and study areas (F3,149 = 1.57, p = 0.19). Mean IFBFat was 6.8% (SE = 0.12%) and ranged from 3.9 to 10.7%. Ingesta-free body fat of females was similar among years (F1,149 = 0.09, p = 0.76) and study areas (F3,149 = 2.06, p = 0.11). Figure 4. Probability of spring migration relative to environmental variables. Model-averaged estimates of the daily probability of migration (dark grey shaded region) and cumulative proportion migrated (light grey shaded region), average daily temperature (uC), daily snow depth (cm), and normalized difference vegetation index (NDVI) relative to Julian date of departure from winter range (a) and arrival on summer range (b) during spring migration of adult female mule deer, Piceance Basin, Colorado, USA, 2010. doi:10.1371/journal.pone.0064548.g004 probability of arrival on summer range increased as snow depths declined, temperatures warmed, and NDVI increased (Figs. 4, 5). Although NDVI on summer range had a positive effect on the probability of arrival on summer range, that relationship was dampened by increases in snow depth (Table 2). Individuals arrived 4.7 days earlier than average when their summer ranges displayed high values of NDVI (95% quantile) compared with a 4day delay from the mean for individuals with summer ranges with low NDVI levels (5% quantile; Figure 5). Intrinsic Variables. PC1, reflecting distance, rate of travel, and well-pad density, affected daily probability of departure from winter range (Table 1). Deer that traveled greater distances, at slower rates, through less natural-gas development (5% quantile) departed winter range 3.5 days earlier than average, whereas deer that traveled faster, over shorter distances, through greater levels of natural-gas development (95% quantile) remained on winter range 5.5 days longer than average (Figure 5). That same intrinsic variable (PC1) also affected timing of arrival on summer range; however, the relationship was the opposite of that for departure from winter range (Table 2). Deer that traveled slower for greater distances, through less development (5% quantile) arrived on summer range later than deer that traveled faster, over shorter distances, through greater development (95% quantile, Figure 5). PC1 had an effect similar in magnitude to that of NDVI on timing of arrival, with a deviation from the mean departure date of approximately 5.5 days (Figure 5). Elevation (PC2) of summer and winter ranges did not have a direct effect on probability of migration (Tables 1, 2). In a post-hoc analysis, however, we exchanged NDVI for elevation from the top model; elevation of Predictors of Spring Migration Extrinsic Variables. Departure from winter range was associated with daily absolute snow depth and temperature (WPC2), daily solar radiation (WPC4), daily change in snow depth (WPC5), and NDVI on winter range (WR-NDVI; importance weight .0.70; Table 1). Daily probability of leaving winter range was higher with reduced snow depth, warmer temperatures, cloudier days, and greater values of NDVI (Figs. 4, 5). Deer that wintered in areas with early increases in NDVI (95% quantile) left winter range 4.5 days prior to the mean departure day, whereas deer that experienced slower spring green-up (5% quantile) left winter range 6.5 days after mean departure day (Figure 5). Timing of arrival on summer range was affected by daily change in temperature (WPC1) and snow depth (WPC5), NDVI on summer range (SR-NDVI), and an interaction between absolute snow depth and temperature (WPC2) and NDVI of summer range (SR-NDVI; Table 2). Similar to departure from winter range, PLOS ONE | www.plosone.org 6 May 2013 | Volume 8 | Issue 5 | e64548 �Mule Deer Migration Patterns Initiation of spring migration was linked closely to patterns of local weather and plant phenology. Migration was delayed with increased winter severity associated with decreased temperature, increased humidity, and increased snow depth. Deep snow along migratory routes hampered travel of migrating deer by delaying their arrival to summer range, even if green-up was occurring at their destination. Indeed, migration strategies of ungulates are plastic and individuals delay migration in years with heavy snow pack and late green-up in spring, and migrate early in years with low snow pack and early emergence of vegetation [3,12,17]. Although weather and plant phenology helped synchronize spring migration, anthropogenic disturbance affected how individuals responded to such external stimuli. Female mule deer departed winter range and arrived on summer range following patterns of green-up (NDVI) on each respective range. Patterns of vegetation green-up on winter range affected timing of departure, whereas arrival on summer range was related to green-up on those areas. The ability of ungulates to follow the ‘‘green wave’’ along spatial and elevational gradients is a well-documented phenomenon, and is a critical aspect of their behavioral ecology, which allows individuals to enhance nutritional gain via access to high-quality forage during a crucial time of year [12,17,36]. Delayed migration to high-elevation ranges was more accurately described by patterns of NDVI. Indeed, when we exchanged NDVI for elevation from the top model, elevation of summer range was negatively related to timing of arrival, which supports the hypothesis that mule deer are following emerging vegetation along spatial and elevational gradients. Garrott et al. [3] hypothesized that deer must first improve their physiological condition prior to incurring the energetic costs associated with spring migration. In our study, female mule deer in better nutritional condition departed winter range earlier than females with lower ingesta-free body fat (IFBFat). Locomotive costs associated with migration are much less costly for large ungulates [37] compared with avian taxa, which may deplete substantial somatic reserves during migration [38]. Rather, we hypothesize that female mule deer with depleted nutritional reserves exhibited risk-averse behavior by remaining on winter range longer, where forage resources were likely less palatable and diverse, but more predictable (sensu [12]). Moreover, females in better nutritional condition that left winter range earlier risked encountering deep snow or late-winter storms at higher elevations, but could have been rewarded with newly emergent vegetation, thereby exhibiting risk-prone behavior [12]. Consequently, individual migration decisions can result in fitness inequalities, which may influence population dynamics [39]. Characteristics of migration routes of individuals affected how well their migration was synchronized with changes in weather and plant phenology; well-pad density and distance traveled likely affected timing of migration indirectly by influencing rate of travel during migration. Deer from the least-developed areas traveled slower over greater distances compared with deer that migrated through more developed areas over shorter distances – an effect that was comparable to green-up on timing of migration (Figure 5). As reported by Lendrum et al. [15], mule deer migrating through the study area of greatest well-pad density had longer step lengths compared with deer migrating through the least-developed areas. This outcome is consistent with the observations noted herein; that mule deer traveling through areas of high development traveled faster, and as a result, departed winter range later but still arrived to summer range earlier compared with deer migration through areas of lesser development. Unfortunately, we were not able to disentangle rate of travel and well-pad density from distance traveled, and had no pre-development data, which might have Figure 5. Predicted effect of migration under varying conditions. Relative difference in days of mean departure and arrival of adult female mule deer for variables identified to influence spring migration in the Piceance Basin Piceance Basin, Colorado, USA, during 2010. Estimates illustrate the predicted effects for the upper 95% and lower 5% quantiles of the normalized difference vegetation index (NDVI); distance traveled, rate of travel, well-pad density (PC1); and ingesta-free body fat (IFBFat), relative to the average date of migration. Predicted effects of each respective variable were determined by holding all other variables constant at their mean. doi:10.1371/journal.pone.0064548.g005 summer range was significant (90% CI not overlapping 0) and positively related (b = 0.001) to timing of arrival on summer range. Life-history Characteristics. Only IFBFat influenced probability of departure from winter range. Deer in better nutritional condition (95% quantile) were more likely to initiate spring migration earlier than deer with low body fat (5% quantile), which departed winter range 5 days later than deer in the upper 95% quantile (Figure 5). Timing of departure differed significantly among study areas, as did year for departure and arrival models (Tables 1, 2), indicating that some other factors among study areas and years still remained unexplained by our predictor variables. Discussion Migration by large herbivores is thought to be a strategy aimed at enhancing fitness, by increasing access to food, escape from predators, and avoidance of risky environmental conditions [8,14]. We evaluated whether anthropogenic disturbance modified patterns of spring migration of mule deer while accounting for other factors known to underpin migration of large herbivores. PLOS ONE | www.plosone.org 7 May 2013 | Volume 8 | Issue 5 | e64548 �Mule Deer Migration Patterns Table 1. Model-averaged parameter estimates, 90% CI, and Akaike importance weights for interval-censored models describing the relationship between the daily probability of departure from winter range for spring migration of mule deer, Piceance Basin, Colorado, USA, 2008–2010. 90% CI Model Group Parameter Extrinsic Intrinsic Life-history Estimate Lower Upper Importance Weight Significant WPC1 20.07 20.17 WPC2 23.57 24.11 0.03 0.52 No 23.02 1.00 WPC3 20.06 Yes 20.16 0.03 0.49 No WPC4 WPC5 20.64 20.90 20.37 1.00 Yes 20.30 20.49 20.10 0.91 Yes WR-NDVI 5.49 2.67 8.31 1.00 Yes SR-NDVI n/a n/a n/a 0.00 No SA n/a n/a n/a 1.00 Yes Yes Year n/a n/a n/a 1.00 PC1 20.28 20.44 20.12 0.96 Yes PC2 0.01 20.05 0.06 0.28 No IFBF 0.10 0.00 0.19 0.75 Yes Age 0.00 20.02 0.01 0.25 No Extrinsic variables included: change in temperature and humidity (WPC1), absolute snow depth and daily temperature (WPC2), precipitation (WPC3), solar radiation (WPC4), changes in snow depth (WPC5), normalized difference vegetation index on winter range (WR-NDVI), on summer range (SR-NDVI), and study area (SA). Intrinsic variables included: distance migrated, average daily distance traveled, well-pad density along migration routes (PC1), and average elevation on summer range (PC2). Individual covariates were ingesta-free body fat (IFBFat), age in years (Age). Year (Year) was included as a nuisance parameter. doi:10.1371/journal.pone.0064548.t001 further clarified this outcome. Nevertheless, ungulates commonly avoid areas of human disturbance, including roads and well pads, especially in late winter and at parturition [19,20,40], although well pads were not avoided by mule deer during this study [15]. We hypothesize that instead of avoiding disturbed areas during migration, mule deer altered timing and rate of movement to reduce exposure to disturbance, potentially reducing net energetic gain [37]. This display of behavioral plasticity may be an Table 2. Model-averaged parameter estimates, 90% CI, and Akaike importance weights for interval-censored models describing the relationship between the daily probability of arrival to summer range for spring migration of mule deer, Piceance Basin, Colorado, USA, 2008–2010. 90% CI Model Group Parameter Extrinsic WPC1 WPC3 WPC4 Intrinsic Life-history Estimate Upper Lower Importance Weight Significant 20.12 20.23 20.01 0.76 Yes 20.06 20.15 0.03 0.45 No 20.10 20.20 0.00 0.64 No WPC5 20.60 20.76 20.44 1.00 Yes SR-NDVIxWPC2 23.98 24.67 23.29 1.00 Yes SR-NDVI 2.16 0.68 3.64 1.00 Yes WR-NDVI n/a n/a n/a 0.00 No SA n/a n/a n/a 1.00 No Year n/a n/a n/a 1.00 Yes PC1 0.57 0.42 0.71 1.00 Yes PC2 20.09 20.19 0.01 0.64 No Age 0.00 20.02 0.02 0.34 No IFBF 20.02 20.08 0.05 0.63 No Extrinsic variables included: change in temperature and humidity (WPC1), absolute snow depth and daily temperature (WPC2), precipitation (WPC3), solar radiation (WPC4), changes in snow depth (WPC5), normalized difference vegetation index on winter range (WR-NDVI), on summer range (SR-NDVI), interaction between SR-NDVI and WPC2 (SR-NDVIxWPC2), and study area (SA). Intrinsic variables included: distance migrated, average daily distance traveled, well-pad density along migration routes (PC1), and average elevation on summer range (PC2). Individual covariates were ingesta-free body fat (IFBFat), age in years (Age). Year (Year) was included as a nuisance parameter. doi:10.1371/journal.pone.0064548.t002 PLOS ONE | www.plosone.org 8 May 2013 | Volume 8 | Issue 5 | e64548 �Mule Deer Migration Patterns important strategy where mule deer use traditional routes during migration that are affected by anthropogenic disturbances. To compensate for the increased rates of travel in areas of high disturbance, female deer that migrated faster departed winter range later, but still arrived on summer range earlier than deer that traveled through less-developed areas. This alteration in timing could disconnect timing of migration from phenological progression [41]. For migratory ungulates, failure to time lifehistory events in accordance with advances in plant phenology can have adverse effects on fitness by reducing net energetic gains [42]. Herbivores track phenological progression of forage plants by moving across landscapes to acquire newly emergent vegetation [30], a strategy of particular importance to pregnant females supporting the increasing demands of late gestation [17]. The relatively minor shift in departure and arrival dates as a function of disturbance levels that we observed (,8 days) may not contribute to a long term fitness consequence, but could do so if increased development activity intensified behavioral shifts in migratory patterns of mule deer [18]. Additionally, this potential disconnect between timing of migration and phenological progression could be exacerbated for deer that migrate greater distances or longer durations than those observed in the Piceance Basin. Migration is a critical life-history characteristic of ungulates that is at risk of disruption because of habitat loss and fragmentation, largely resulting from anthropogenic disturbances [4,6]. In some situations, the advantages acquired by migration could be outweighed by the risk, additional time, and energetic costs associated with avoidance of increased human development [5,7,18]. Despite high levels of energy development in the Piceance Basin, local weather patterns and plant phenology remained the predominant factors driving patterns of migration in mule deer; however, exposure to disturbance altered how individual deer responded to those environmental factors. As the level of natural-gas development expands across the Intermountain West, large areas of habitat for mule deer are being rapidly converted into gas fields consisting of networks of access roads, well pads, pipelines, and other infrastructure, which have potential to alter migratory behavior [18,20]. Mule deer in the Piceance Basin appear to avoid negative effects from development activity through behavioral shifts in timing and rate of migration. Continued monitoring of mule deer and energy-development interactions are necessary to identify potential development strategies that minimize behavioral shifts in traditional migratory patterns. Supporting Information Appendix S1 (DOC) Appendix S2 (DOC) Acknowledgments We thank C. Bishop, D. Freddy, and M. Michaels from Colorado Parks and Wildlife (CPW) for helping administer the project. Additionally, we thank personnel at Little Hills State Wildlife Area for field support and coordination. We thank Quicksilver Air Inc. for assistance in capturing deer from helicopters, L. Gepfert and L. Coulter for fixed-wing aircraft support, and L. Wolfe, C. Bishop, and D. Finley of CPW for crucial assistance during capture efforts. We also thank the White River Bureau of Land Management, U.S. Forest Service, and numerous private land owners for their cooperation. In addition, we thank J. Kie, J. Shallow, S. Collins, N. Guernsey, R. Long, and E. Bergman, who reviewed earlier versions of this manuscript and provided valuable comments. Author Contributions Conceived and designed the experiments: PEL CRA. Performed the experiments: PEL CRA. Analyzed the data: PEL KLM. Contributed reagents/materials/analysis tools: PEL CRA KLM JAJ RTB. Wrote the paper: PEL CRA KLM JAJ RTB. References 14. Mysterud A, Loe LE, Zimmerman B, Bischof R, Beiberg V, et al. (2011) Partial migration in expanding red deer populations at northern latitudes – a role for density dependence. Oikos 120: 1817–1825. 15. Lendrum PE, Anderson CR Jr, Long RA, Kie JG, Bowyer RT (2012) Habitat selection by mule deer during migration: effects of landscape structure and natural-gas development. Ecosphere 3: art82. 16. Mysterud A, Langvatn R, Yoccoz NG, Stenseth NC (2001) Plant phenology, migration and geographical variation in body weight of a large herbivore: the effect of a variable topography. Journal of Animal Ecology 70: 915–923. 17. Parker KL, Barboza PS, Gillingham MP (2009) Nutrition integrates environmental responses of ungulates. Functional Ecology 23: 57–69. 18. Sawyer H, Kauffman MJ, Middleton AD, Morrison TA, Nielson RM, et al. (2013) A framework for understanding semi-permeable barrier effects on migratory ungulates. Journal of Applied Ecology 50: 68–78. 19. Sawyer H, Nielson RM, Lindzey FG, McDonald LL (2006) Winter habitat selection of mule deer before and during development of a natural gas field. Journal of Wildlife Management 70: 396–403. 20. Sawyer H, Kauffman MJ, Nielson RM (2009) Influence of well pad activity on winter habitat selection patterns of mule deer. Journal of Wildlife Management 73: 1052–1061. 21. Sikes RS, Gannon WL, and the Animal Care and Use Committee of the American Society of Mammalogists (2011) Guidelines of the American Society of Mammalogists for the use of wild mammals in research. Journal of Mammalogy 92: 235–253. 22. White GC, Lubow BC (2002) Fitting population models to multiple sources of observed data. Journal of Wildlife Management 66: 300–309. 23. Krausman PR, Hervert JJ, Ordway LL (1985) Capturing deer and mountain sheep with a net-gun. Wildlife Society Bulletin 13: 71–73. 24. D’Eon RG, Delparte D (2005) Effects of radio-collar position and orientation on GPS radio-collar performance, and the implications of PDOP in data screening. Journal of Applied Ecology 42: 383–388. 25. Stephenson TR, Bleich VC, Pierce BM, Mulcahy GP (2002) Validation of mule deer body composition using in vivo and post-mortem indices of nutritional condition. Wildlife Society Bulletin 30: 557–564. 1. Andersen R (1991) Habitat deterioration and the migratory behavior of moose in Norway. Journal of Applied Ecology 28: 102–108. 2. Sawyer H, Kauffman MJ, Nielson RM, Horne JS (2009) Identifying and prioritizing ungulate migration routes for landscape-level conservation. Ecological Applications 19: 2016–2025. 3. Garrott RA, White GC, Bartmann RM, Carpenter LH, Alldredge AW (1987) Movements of female mule deer in northwest Colorado. Journal of Wildlife Management 51: 634–643. 4. Bolger DT, Newmark WD, Morrison TA, Doak DF (2008) The need for integrative approaches to understand and conserve migratory ungulates. Ecology Letters 11: 63–77. 5. Harris G, Thirgood S, Hopcraft G, Cromsigt JPGM, Berger J (2009) Global decline in aggregated migrations of large terrestrial mammals. Endangered Species Research 7: 55–76. 6. Berger J (2004) The last mile: how to sustain long-distance migration in mammals. Conservation Biology 18: 320–331. 7. Singh NJ, Milner-Gulland EJ (2011) Conserving a moving target: planning protection for a migratory species as its distribution changes. Journal of Applied Ecology 48: 35–46. 8. Fryxell JM, Sinclair ARE (1988) Causes and consequences of migration by large herbivores. Trends in Ecology and Evolution 3: 237–241. 9. Boyce MS (1991) Migratory behavior and management of elk (Cervus elaphus). Applied Animal Behavior Science 29: 239–250. 10. Nicholson MC, Bowyer RT, Kie JG (1997) Habitat selection and survival of mule deer: tradeoffs associated with migration. Journal of Mammalogy 78: 483– 504. 11. Albon SD, Langvatn R (1992) Plant phenology and the benefits of migration in a temperate ungulate. Oikos 65: 502–513. 12. Monteith KL, Bleich VC, Stephenson TR, Pierce BM, Conner MM, et al. (2011) Timing of seasonal migration in mule deer: effects of climate, plant phenology, and life-history characteristics. Ecosphere 2: art47. 13. Mysterud A (1999) Seasonal migration pattern and home range of roe deer (Capreolus capreolus) in an altitudinal gradient in southern Norway. Journal of Zoology 247: 479–486. PLOS ONE | www.plosone.org 9 May 2013 | Volume 8 | Issue 5 | e64548 �Mule Deer Migration Patterns 26. Cook RC, Cook JC, Stephenson TR, Myers WI, McCorquodale SM, et al. (2010) Revisions of rump fat and body scoring indices for deer, elk, and moose. Journal of Wildlife Management 74: 880–896. 27. Robinette WL, Jones DA, Rogers G, Gashwiler JS (1957) Notes on tooth development and wear for Rocky Mountain mule deer. Journal of Wildlife Management 21: 134–153. 28. Whittington J, St. Clair CC, Mercer G (2004) Path tortuosity and the permeability of roads and trails to wolf movement. Ecology and Society 9: art4. 29. Grovenburg TW, Jenks JA, Klaver RW, Swanson CC, Jacques CN, et al. (2009) Seasonal movements and home ranges of white-tailed deer in north-central South Dakota. Canadian Journal of Zoology 87: 876–885. 30. Hebblewhite M, Merrill E, McDermid G (2008) A multi-scale test of the forage maturation hypothesis in a partially migratory ungulate population. Ecological Monographs 78: 141–166. 31. Pettorelli N, Vik JO, Mysterud A, Gaillard JM, Tucker CJ, et al. (2005) Using satellite-derived NDVI to assess ecological responses to environmental change. Trends in Ecology and Evolution 20: 503–510. 32. White GC, Burnham KP (1999) Program MARK: survival estimation from populations of marked animals. Bird Study 46: 120–139. 33. Fieberg J, DelGiudice G (2008) Exploring migration data using interval-censored time-to event models. Journal of Wildlife Management 72: 1211–1219. 34. St-Louis A, Hamel S, Mainguy J, Côté SD (2013) Factors influencing the reaction of mountain goats towards all-terrain vehicles. Journal of Wildlife Management 77: 599–605. PLOS ONE | www.plosone.org 35. Burnham KP, Anderson DR (2002) Model selection and multimodel inference: a practical information theoretic approach. Second edition. Springer-Verlag: New York. 36. Sawyer H, Kauffman MJ (2011) Stopover ecology of a migratory ungulate. Journal of Animal Ecology 80: 1078–1087. 37. Parker KL, Robbins CT, Hanley TA (1984) Energy expenditures for locomotion by mule deer and elk. Journal of Wildlife management 48: 474–488. 38. Pulido F (2007) Phenotypic changes in spring arrival: evolution, phenotypic plasticity, effects of weather and condition. Climate Research 35: 5–23. 39. Alves JA, Gunnarsson TG, Potts PM, Gélinaud G, Sutherland WJ, et al. (2012) Overtaking on migration: does longer distance migration always inqure a penalty? Oikos 121: 464–470. 40. Cameron RD, Smith WT, White RG, Griffith B (2005) Central arctic caribou and petroleum development: distributional, nutritional, and reproductive implications. Arctic 58: 1–9. 41. Bishof R, Loe LE, Meisingset EL, Zimmermann B, Moorter BV, et al. (2012) A migratory northern ungulate in the pursuit of spring: jumping of surfing the green wave. American Naturalist 180: 407–424. 42. Post E, Pederson C, Wilmers CC, Forchhammer MC (2008) Warming, plant phenology and the spatial dimension of trophic mismatch for large herbivores. Proceedings of the Royal Society of London B 275: 2005–213. 10 May 2013 | Volume 8 | Issue 5 | e64548 � https://cpw.cvlcollections.org/files/original/b19016fb443a5d8cbc9e7cbf714c0b67.pdf 74f1bc13113d53b16a0ec76e209ae37b PDF Text Text Appendix I. Component loadings for principal component analysis for extrinsic and intrinsic variables used to model timing of spring migration for mule deer in the Piceance Basin, Colorado, USA, 2008–2010. Variables Extrinsic Snow depth ∆ Snow Temperature ∆ Temperature Humidity ∆ Humidity Solar radiation ∆ Solar ∆ Precipitation Precipitation Intrinsic Distance Rate Well density Elevation 1 51.49 6.60 -4.00 -0.32 3.40 0.95 -0.51 0.11 0.07 -0.19 Component loadings 2 3 4 2.33 -1.88 0.37 -1.86 14.89 0.75 3.79 0.27 3.29 3.64 0.26 0.23 -20.96 0.41 -4.62 -20.97 -1.81 5.22 1.38 0.02 0.34 0.32 0.12 -0.02 -1.33 -0.23 0.38 -1.32 -0.17 0.08 -0.81 0.75 0.73 0.03 0.14 -0.20 0.33 0.95 5 0.19 -0.19 3.19 2.04 1.60 -0.68 -0.15 -0.11 0.21 0.45 � https://cpw.cvlcollections.org/files/original/965ab3f8019d9ec02e7754f27b9d3ddf.pdf 86d52c1f0fad1259fdfb07a97358f002 PDF Text Text Appendix II. Mean date of departure and arrival, and weather variables for 10 days prior to and post mean departure and arrival dates of spring migration of mule deer, Piceance Basin, Colorado, USA, 2008–2010. Snow depth was collected from SNOTEL weather station; temperature, humidity, and solar radiation, was collected from Western Regional Climate Center. Year Departure Day Mean NDVI SD Mean SE Snow Depth (cm) Temperature (°C) Mean Mean SE SE Humidity (%) Mean Solar Radiation (Watts/m2) SE Mean SE 2008 134.30 10.62 0.34 0.001 77.05 3.91 10.11 1.12 50.24 4.27 6.69 0.41 2009 123.90 5.84 0.31 0.001 57.69 5.13 10.37 0.82 41.57 4.28 7.26 0.43 2010 132.50 7.75 0.34 0.001 61.21 4.02 8.38 0.87 50.70 4.11 6.88 0.40 2008 142.40 9.81 0.42 0.002 59.75 5.55 11.25 1.11 45.86 4.41 7.16 0.46 2009 131.10 6.16 0.38 0.002 20.68 4.79 12.67 0.81 38.00 3.44 7.51 0.42 2010 139.50 7.72 0.40 0.001 37.62 5.19 10.35 1.07 48.33 4.47 6.86 0.39 Arrival � Dublin Core The Dublin Core metadata element set is common to all Omeka records, including items, files, and collections. For more information see, http://dublincore.org/documents/dces/. Title A name given to the resource Journal Articles Description An account of the resource CPW peer-reviewed journal publications Text A resource consisting primarily of words for reading. Examples include books, letters, dissertations, poems, newspapers, articles, archives of mailing lists. Note that facsimiles or images of texts are still of the genre Text. Dublin Core The Dublin Core metadata element set is common to all Omeka records, including items, files, and collections. For more information see, http://dublincore.org/documents/dces/. Title A name given to the resource Migrating mule deer: effects of anthropogenically altered landscapes Rights Information about rights held in and over the resource <a href="http://rightsstatements.org/vocab/InC-NC/1.0/" target="_blank" rel="noreferrer noopener">In Copyright - Non-Commercial Use Permitted</a> Date Created Date of creation of the resource. 2013-05-14 Subject The topic of the resource Migration <em>Odocoileus hemionus</em> Northwestern Colorado Oil and gas development Description An account of the resource <em>Background</em>: Migration is an adaptive strategy that enables animals to enhance resource availability and reduce risk of predation at a broad geographic scale. Ungulate migrations generally occur along traditional routes, many of which have been disrupted by anthropogenic disturbances. Spring migration in ungulates is of particular importance for conservation planning, because it is closely coupled with timing of parturition. The degree to which oil and gas development affects migratory patterns, and whether ungulate migration is sufficiently plastic to compensate for such changes, warrants additional study to better understand this critical conservation issue.<br /><br /><em>Methodology/Principal Findings</em>: We studied timing and synchrony of departure from winter range and arrival to summer range of female mule deer (Odocoileus hemionus) in northwestern Colorado, USA, which has one of the largest natural-gas reserves currently under development in North America. We hypothesized that in addition to local weather, plant phenology, and individual life-history characteristics, patterns of spring migration would be modified by disturbances associated with natural-gas extraction. We captured 205 adult female mule deer, equipped them with GPS collars, and observed patterns of spring migration during 2008–2010.<br /><br /><em>Conclusions/Significance</em>: Timing of spring migration was related to winter weather (particularly snow depth) and access to emerging vegetation, which varied among years, but was highly synchronous across study areas within years. Additionally, timing of migration was influenced by the collective effects of anthropogenic disturbance, rate of travel, distance traveled, and body condition of adult females. Rates of travel were more rapid over shorter migration distances in areas of high natural-gas development resulting in the delayed departure, but early arrival for females migrating in areas with high development compared with less-developed areas. Such shifts in behavior could have consequences for timing of arrival on birthing areas, especially where mule deer migrate over longer distances or for greater durations. Creator An entity primarily responsible for making the resource Lendrum, Patrick E. Anderson Jr, Charles R. Monteith, Kevin L. Jenks, Jonathan A. Bowyer, R. Terry Language A language of the resource English Is Part Of A related resource in which the described resource is physically or logically included. PLoS One Format The file format, physical medium, or dimensions of the resource application/pdf Extent The size or duration of the resource. 10 pages Bibliographic Citation A bibliographic reference for the resource. Recommended practice is to include sufficient bibliographic detail to identify the resource as unambiguously as possible. Lendrum, P. E., C. R. Anderson, Jr., K. L. Monteith, J. A. Jenks and R. T. Bowyer. 2013. Migrating mule deer: effects of anthropogenically altered landscapes. PLoS One 8(5): e64548. <a href="https://doi.org/10.1371/journal.pone.0064548" target="_blank" rel="noreferrer noopener">https://doi.org/10.1371/journal.pone.0064548</a> Type The nature or genre of the resource Article