https://cpw.cvlcollections.org/files/original/26291721d2ca0b42e9087af70db8184d.pdf f886401fb409412bc388e05839f9f25c 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 �Herbivore Body Condition Response in Altered Environments: Mule Deer and Habitat Management Eric J. Bergman1*, Paul F. Doherty Jr.2, Chad J. Bishop1, Lisa L. Wolfe1, Bradley A. Banulis3 1 Mammals Research, Colorado Parks and Wildlife, Fort Collins, Colorado, United States of America, 2 Department of Fish, Wildlife, and Conservation Biology, Colorado State University, Fort Collins, Colorado, United States of America, 3 Terrestrial Programs, Colorado Parks and Wildlife, Montrose, Colorado, United States of America Abstract The relationships between habitat, body condition, life history characteristics, and fitness components of ungulates are interwoven and of interest to researchers as they strive to understand the impacts of a changing environment. With the increased availability of portable ultrasound machines and the refinement of hormonal assays, assessment of ungulate body condition has become an accessible monitoring strategy. We employed body condition scoring, estimation of % ingestafree body fat (%IFBF), assessment of free thyroid hormones (FT4 and FT3), and assessment of pregnancy, as metrics to determine if landscape-level habitat treatments affected body condition of adult ($1.5 years old) female mule deer (Odocoileus hemionus). All body condition related metrics were measured on 2 neighboring study areas — a reference area that had received no habitat treatments and a treatment study area that had received mechanical removal of pinyon pine (Pinyus edulis) - Utah juniper (Juniperus osteosperma) forest, chemical control of weeds, and reseeding with preferred mule �deer browse species. �A consistent trend of higher �%IFBF was observed �in the treatment study area � IFBF ~7:38,SD~1:31 than in the reference study area % � IFBF ~6:97,SD~2:16 , although variation of estimates % was larger than hypothesized. A similar pattern was observed with higher thyroid hormones concentrations being observed in the treatment study area, but large amounts of variation within concentration estimates were also observed. The consistent pattern of higher body condition related estimates in our treatment study area provides evidence that large mammalian species are sensitive to landscape change, although variation within estimates underlie the challenge in detecting population level impacts stemming from environmental change. Citation: Bergman EJ, Doherty PF Jr, Bishop CJ, Wolfe LL, Banulis BA (2014) Herbivore Body Condition Response in Altered Environments: Mule Deer and Habitat Management. PLoS ONE 9(9): e106374. doi:10.1371/journal.pone.0106374 Editor: Bernhard Kaltenboeck, Auburn University, United States of America Received April 7, 2014; Accepted August 5, 2014; Published September 3, 2014 Copyright: ß 2014 Bergman 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. Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the Supporting Information files. Funding: Funding for this research was provided by Colorado Parks and Wildlife Game Cash Funds, USFWS Federal Aid Research Grants, Colorado Parks and Wildlife Habitat Partnership Program, the Mule Deer Foundation, and Colorado Parks and Wildlife Big Game Auction and Raffle Grants. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * Email: eric.bergman@state.co.us the first are directly related to body condition of adult females ($1 year old). Thus, the relationships between habitat, body condition and life history characteristics are tightly interwoven [6,7]. Despite this broad history of investigation, the body condition and fitness of ungulates has not been used as a tool for evaluating habitat management. The need to evaluate individual-level metrics as a response to environmental change rests on the assumption that the effects of environmental change may be subtle. These subtle effects may have long-term fitness consequences that remain undetected at short time intervals when assessed with population level monitoring strategies. The increased availability of portable ultrasound machines, coupled with the development and validation of robust body condition estimation models [8–10] has made the assessment of ungulate body condition, and other fitness components, accessible monitoring strategies. Similarly, thyroid hormone concentrations reflect the metabolic condition of ungulates [11– 13], providing a window for assessing an individual’s ability to cope with current environmental conditions. Introduction Natural succession, climate mediated habitat change, deliberate habitat improvement, and direct habitat loss result in changing environments for wildlife populations. Due to the economic and social value of large ungulates, and in turn the datasets that management of these species foster, large ungulate populations are attractive to researchers hoping to elucidate the impacts of environmental change. Yet while wildlife professionals hasten to document the impacts of environmental change, the best barometer for measuring impacts to individuals and populations remains elusive. In general, the cascading effect of habitat quantity and quality on wildlife fitness has received attention for several decades [1–5]. Specifically, in bottom-up systems, the predicted sequence of density-dependent effects experienced by mammals as their populations saturate a landscape and approach the local carrying capacity have been succinctly predicted [1,2]: 1) reduced survival of juveniles, 2) delay in age of first pregnancy, 3) reduced neonatal and parturition rates of adults, and finally 4) reduced survival of adults. These predictions have subsequently been applied to large ungulate species [3,4]. Of these predictions, all but PLOS ONE | www.plosone.org 1 September 2014 | Volume 9 | Issue 9 | e106374 �Herbivore Body Condition and Habitat Management Total body fat and thyroid hormones can be viewed as metrics for the same general trait, overall deer health; however, they are parameters for different processes. Total body fat estimates, generated by merging ultrasonic rump fat measurements with body condition scoring [10], reflect the energetic reserve for an individual deer. Thyroid hormone concentrations reflect the ability of deer to utilize body fat reserves. Thyroxine (T4) hormone is a product of the thyroid gland and is a precursor to the triiodothyronine (T3) hormone [14,15]. The T3 hormone plays a direct role in regulating the basal metabolic rate and thermal regulation within animals [14,15]. Measurements of these 2 hormones typically occur in 2 forms, total hormone concentrations (T4 and T3) and free hormone concentrations (FT4 and FT3). Variation in hormone concentrations is indicative of physiological adjustment to changes in the environment. As managers implement habitat management actions, or as they consider alternative large scale changes to habitat (e.g., habitat response to wildfire or habitat alteration due to development), they often wish to know if ungulate populations have been affected. Experimental research has demonstrated a strong connection between maternal body condition (i.e, %IFBF and hormonal concentrations), pregnancy rates, as well as neonate and juvenile survival, when food was supplemented [7]. However, the study of Bishop et al [7] was designed to explore an ecological process, not to test a practical management scenario. In an attempt to evaluate those results [7,13] in the context of common habitat management techniques, we conducted a study that assessed late-winter body condition of adult female mule deer with respect to such management techniques. We employed body condition scoring, estimation of total body fat and assessment of thyroid hormones, as metrics to determine if landscape-level habitat manipulation affected body condition of adult ($1.5 years old) female mule deer. We hypothesized that estimates of these late winter condition metrics for adult females on the treatment study area would be consistent with animals in better overall condition, although we also hypothesized that our estimates would be lower than the experimentally elevated estimates reported in other research because increasing browse availability to similar ad libitum levels used by Bishop et al. was not a realistic expectation for our habitat management techniques [7,13]. Both study areas were located within Colorado Parks and Wildlife (CPW) Data Analysis Unit (DAU) 40. This 2,437 km2 DAU was managed for a post-hunt population size of 13,500– 15,000 mule deer. Each of these study areas was centered on public lands, although Buckhorn had private land at lower elevations. Likewise, both study areas declined in elevation from east to west. Mule deer arrival on each study area each winter was believed to have been heavily influenced by the building snowpack at higher elevations. Grazing pressure from domestic livestock was minimal on both study areas, with the majority of grazing occurring as livestock producers moved animals from summer range pastures to private pastures in the valley. Due to the proximity of the study areas, and to the overall topography, a high degree of spatial overlap on summer range occurred among deer that used these 2 distinct winter range segments (E. Bergman, Colorado Parks and Wildlife, unpublished data). Due to this mixing on summer range, we assumed that body condition was equalized among deer prior to their arrival on our winter range study areas. Habitat treatments For our research, habitat treatments occurred on BCSWA in 2 stages. The first stage occurred in 1998, during which 135.98 ha (,5%) of the 2,688 ha study area was exposed to mechanical roller-chopper treatments (see Fig. 1). Roller-chopper treatments consisted of a large drum, affixed with perpendicular blades, that was pulled behind a bulldozer [16]. Standing trees and taller vegetation were uprooted by the bulldozer and subsequently broken into smaller pieces by the drum. On BCSWA, rollerchopper treatments ranged in size between 6.8–24.7 ha with the objective of opening the forest canopy and increasing the edge/ area ratio (see Fig. 1). Treatments also created a mulched ground cover that was beneficial for holding moisture. Our second stage of habitat treatment efforts were reseeding and weed control that occurred concurrently with our study (2006–2008). Reseeding efforts targeted desirable browse species for mule deer by planting bitterbrush, cliffrose, sagebrush, serviceberry (Amelanchier alnifolia), and four-wing saltbush (Atriplex canescens). Weed control, via herbicide application, targeted cheatgrass (Bromus tectorum) and jointed goatgrass (Aegilops cylindrical). The second stage of habitat treatments occurred on 78.34 ha (,57%) of the original treatments. The delay between the first stage of habitat treatments and the initiation of mule deer body condition monitoring was a deliberate decision to accommodate temporal variation in vegetation response, post-treatment [17–20]. By allowing a time lag we afforded browse species ample opportunity to grow. Materials and Methods Study area We conducted this research on 2 study areas near the southeastern tip of the Uncompahgre Plateau, near the city of Montrose in southwest Colorado (38u 28.7269, 107u52.6249 – Fig. 1). One study area (Buckhorn) was maintained as a reference area, whereas the second study area (Billy Creek State Wildlife Area – BCSWA) was a treatment area (see Fig. 1). For both study areas, scientific collection was permitted by Colorado Parks and Wildlife (Licence No: CPW001 and CPW003). No component of this research involved endangered or protected species. The study areas were located in close proximity to one another to minimize spatial variation with Buckhorn being approximately 8.5 km north of BCSWA. Each study area was located on pinyon pine - Utah juniper forest winter range. These forests were lateseral stage, typified by open understory with occasional sagebrush (Artemisia spp.), cliffrose (Purshia Mexicana), antelope bitterbrush (Purshia tridentate), mountain mahogany (Cercocarpus spp.), or rabbitbrush (Ericameria spp.) plants. Grasses included western wheatgrass (Pascopyrum smithii), green needlegrass (Nassella viridula), Indian ricegrass (Achnatherum hymenoides) and bluegrass (Poa spp.). PLOS ONE | www.plosone.org Mule deer capture and handling During early March of each winter (2007–2009), 30 adult female deer were captured via helicopter net-gunning [21,22]. Upon capture, all deer were immediately blind-folded, hobbled and ferried to a central processing site (#3.2 km). At the field processing site, deer were weighed, age was estimated via tooth eruption and wear patterns [23–25], hind foot length was measured, and blood was drawn via jugular venipuncture. We measured the maximum subcutaneous fat thickness (cm) on the rump and the thickness of the longissimus dorsi muscle (cm) using a Sonovet 2000 (Universal Medical Systems, Bedford Hills, New York, USA) portable ultrasound machine and a 5-MHz linear transducer [8,9,26,27]. We also determined a body condition score for each animal by palpating the rump [9,10,28]. Capture, handling and radio-collaring procedures for all aspects of this study were approved by the Institutional Animal Care and Use 2 September 2014 | Volume 9 | Issue 9 | e106374 �Herbivore Body Condition and Habitat Management Figure 1. Location of 2 study areas in southwest Colorado. The Buckhorn study area was located in Montrose county, designated by hash marks from lower left to upper right. The Billy Creek study area was located in Ouray county, designated by hash marks from lower right to upper left. Habitat treatments are depicted by red polygons inside the Billy Creek study area boundaries. doi:10.1371/journal.pone.0106374.g001 Committees at CPW (protocol #10-2005) and Colorado State University (protocol #08-2006A). Body condition scores were combined with ultrasound measurements to generate a scaled estimate of the total percent of the body that was ingesta-free body fat (%IFBF) for each animal [10]. At the time of capture, pregnancy was determined via transabdominal ultrasonography [29–31] or via pregnancy-specific protein B concentrations [32] from blood serum samples (Biotracking, LLC, Moscow, Idaho, USA). Blood serum samples were also submitted to the Diagnostic Center for Population and Animal Health at Michigan State University (East Lansing, Michigan, USA) for estimation of T4, FT4, T3, and FT3 concentrations. 23 May 2014). While individual mass was collected for animals at the time of capture, these data were not directly used in the estimation process for %IFBF. For each of the 3 response variables, a total of 64 models were compared. These 64 models comprised a balanced model set in which all response variables were included in an equal number of models. Models containing multiplicative interactions were not included in these model sets. This modeling approach facilitated the computation and comparison of cumulative model weights for each predictor variable and response variable combination [34,35]. To assess the effect of habitat treatments, year, %IFBF, and age on pregnancy, we modeled the probability of an individual’s pregnancy status using logistic regression in the ‘‘Stats’’ package in R. To determine if there was evidence for a delay in age of first pregnancy, or senescence in pregnancy, second and third order polynomial models were also built. Finally, we conducted post hoc exploratory analyses to evaluate the conclusions and recommendations drawn by similar research [13] that regarded the utility of using blood serum thyroid concentrations to estimate %IFBF. This research [13] reported that the T4 and FT4 hormones were effective at predicting %IFBF (%IFBF = 24.8015–0.09466T4+ 0.0006036T42+0.14746FT4+0.14266chestgirth, R2 = 0.609). Following those methods [13], second and third order polynomials were allowed to occur in our later models. For all model sets, model fit was examined through residual and QQ plots. Assumptions of linear regression were upheld in all cases. Analytical methods Prior to building body condition models, we tested for correlation between response variables. Based on the results of correlation analyses, we modeled 3 of the 5 body condition measurements (%IFBF, FT4, and FT3) as a response to group covariates (study area and year) and to individual covariates (chest girth, hind foot length, pregnancy status, and age). For all analyses, model selection and evaluation was based on AICc [33]. Conditional model averaging of estimates was conducted such that average parameter estimates were generated using all models. For models in which individual parameters did not appear, b and standard error values of 0 were used. All possible combinations of additive multiple linear regression models were evaluated using the ‘‘MuMin’’, and ‘‘Stats’’ packages in R (R Foundation for Statistical Computing, version 3.1.0. www.r-project.org, accessed PLOS ONE | www.plosone.org 3 September 2014 | Volume 9 | Issue 9 | e106374 �Herbivore Body Condition and Habitat Management Our linear regression models for FT4 were similar to those of %IFBF. The best FT4 model was composed of study area, year, age, and pregnancy status (see Table 3), and estimated concentrations of FT4 for Buckhorn were consistently lower than those estimated for BCSWA (see Table 1 and Table 2). However, while the highest %IFBF was estimated during the second year of the study (2008), FT4 concentrations were highest during the first year of the study. The model-averaged parameter estimate for ^ = 21.739, pregnancy status in FT4 models was negative (b SE = 0.754), and a common factor in many of our top models (see Table 1). The morphometric measurements of chest girth and hind foot length explained little of the variation in our data and accounted for cumulative AICc weights were near 0.50 (see Table 1). The best linear regression model for FT3 deviated from the patterns established by %IFBF and FT4. The role of year and age had the greatest influence on model predictions for FT3, whereas study area did not (see Table 1). Concentration of FT3 was lower during 2008 and 2009, following the pattern observed for FT4. Based on AICc, as well as model-averaged parameter estimates, when pregnancy status was treated as a dependent variable there was no evidence that the probability of an adult female deer being pregnant varied between study areas or during years. Little difference in pregnancy rates was observed between BCSWA and Buckhorn during the 3 year period (BCSWA = 0.877 (SD = 0.329), Buckhorn = 0.862 (SD = 0.345)). When pooled across study areas, observed mean pregnancy was lower in 2008 than in 2007 and 2009 (2007: 0.896 (SD = 0.307), 2008: 0.833 (SD = 0.376), 2009: 0.883 (SD = 0.324)). Probability of being pregnant was best ^ = 3.354–0.27396age, with the effect of predicted by the model P ^ = 20.311, SE = 0.109). Our age on pregnancy being negative (b data did reflect some evidence for late age senescence (Fig. 3). Exploratory models that were structured with second and third order polynomial expressions in an attempt to accommodate delayed age of first pregnancy or late age senescence did not improve upon the simpler additive models. Results Estimated %IFBF was more correlated with T4 (0.25) and FT4 (0.18) than with T3 (0.07) and FT3 (0.09). However, the highest overall correlations were observed within categories of thyroid hormones. T4 and FT4 had the highest correlation (0.89), whereas the correlation between T3 and FT3 was slightly lower (0.70). Correlation of concentrations between the 2 T4 hormones and the 2 T3 hormones were consistently between 0.40–0.45. Correlations among predictor variables were low with the highest observed correlation occurring between individual chest girth and individual hind foot length (0.31). The pooled, mean estimate of %IFBF for deer during this study was 7.17% (SD = 1.79). The observed mean value for BCSWA � = 7.38, SD = 1.31) was higher than Buckhorn (% � = 6.97, (% SD = 2.16). Overall, the effect of year was an important component of model structures for all hormones (Table 1). When %IFBF was compared among years, the mean observed estimate � = 6.85, SD = 1.99) was less than 2008 (% � = 7.48, in 2007 (% � SD = 1.78) or 2009 (% = 7.19, SD = 1.56). The observed pattern of higher %IFBF and T4 in BCSWA was observed during all 3 years, although no pattern for T3 hormone levels was observed (Table 2). The difference in %IFBF between study areas and among years was subtle, with wide overlap in the estimates of variance (Fig. 2). The overall best model incorporated study area, year and individual chest girth (Table 3) and cumulative AICc weights from %IFBF models also reflected the importance of these 3 covariates (see Table 1). Estimated %IFBF was higher in BCSWA than in Buckhorn, and reflected a 1.086magnitude increase when pooled over 3 years. However, the performance of our best %IFBF model was not strong (see Table 3). Much of the variation within these data (.90%) remained unexplained. Based on cumulative AICc weights, the remaining covariates of interest (pregnancy status, hind foot length and age) accounted for less than 0.50 and contributed little to overall model predictions (see Table 1). P ^ ), and standard Table 1. Akaike’s Information Criterion cumulative model weights ( vi ), multiple linear regression coefficients (b i errors for body condition predictor variables for adult female mule deer. Predictor Response Variables Variable %IFBF FT4 0.722 FT3 Unit P 1.000 0.265 20.385 (0.191) 23.433 (0.614) 0.010 (0.028) Year ^ (SE) b i P vi 0.823 1.000 0.999 2008 ^ (SE) b i 0.745 (0.273) 23.447 (0.756) 20.108 (0.130) 2009 ^ (SE) b i P vi 0.267 (0.267) 25.115 (0.743) 20.872 (0.128) Chest 0.966 0.412 0.317 0.089 (0.028) 0.035 (0.028) 0.003 (0.003) Age ^ (SE) b i P vi 0.363 0.890 0.982 20.025 (0.024) 20.351 (0.138) 20.082 (0.026) Foot ^ (SE) b i P vi 0.293 0.533 0.504 0.017 (0.027) 20.170 (0.109) 20.026 (0.018) Pregnant ^ (SE) b i P vi 0.260 0.827 0.452 ^ (SE) b i 20.016 (0.103) 21.739 (0.754) 20.094 (0.071) vi Data were collected in southwest Colorado during early March, 2007–2009. A balanced model set was used such that all response variables appeared in an equal number of models, thus cumulative model weights .0.5 are attributed to variables that are most important. doi:10.1371/journal.pone.0106374.t001 PLOS ONE | www.plosone.org 4 September 2014 | Volume 9 | Issue 9 | e106374 �Herbivore Body Condition and Habitat Management Table 2. Observed mean estimates (with standard deviation) for 5 body condition variables from adult female mule deer in southwestern Colorado. Year Unit %IFBF T4 FT4 T3 FT3 2007 BCSWA 6.82 (1.51) 88.23 (15.53) 14.8 (3.98) 1.55 (0.53) 2.1 (0.7) Buckhorn Mountain 6.81 (2.36) 78.07 (22.34) 13.1 (4.66) 1.42 (0.31) 2.07 (0.56) 2008 BCSWA 7.91 (1.24) 94.3 (20.71) 13.37 (4.59) 1.17 (0.28) 1.98 (0.59) Buckhorn Mountain 7.05 (2.12) 56.17 (23.32) 8.37 (3.91) 1.17 (0.58) 2.13 (1.16) 2009 BCSWA 7.4 (0.94) 74.63 (14.61) 11.33 (3.46) 1.22 (0.32) 1.41 (0.52) Buckhorn Mountain 6.98 (1.99) 54.77 (19.34) 6.83 (3.17) 1.26 (0.35) 1.14 (0.44) Data were collected during early March in a treatment study area (Billy Creek State Wildlife Area – BCSWA) and a reference study area (Buckhorn Mountain). Variables include scaled percent ingesta-free body fat (%IFBF), as well as concentrations for the thyroid hormones: T4 (nanomole/L), T3 (nanomole/L), FT4 (picomole/L), and FT3 (picomole/L). doi:10.1371/journal.pone.0106374.t002 capacity to utilize those reserves (FT4) appeared to be higher in treatment deer than in reference deer. However, the considerable variation that occurred within those estimates tempers this conclusion. For both %IFBF and FT4 results, study area and year were consistently among the most important covariates. In the case of FT4, these covariates accounted for .99% of the cumulative model weight (see Table 1), demonstrating that these covariates were most useful in explaining variation in the data. In the case of %IFBF, the single covariate that explained most of the variation in the data was chest girth. Chest girth, a variable directly related to body size, helped distinguish between large bodied animals that had low %IFBF and small bodied animals that had high %IFBF. Annual variability in body condition among winters was expected to be an important factor in assessing late winter body condition, although its importance relative to habitat management efforts was difficult to predict prior to our study. This expectation was met as Results of our exploratory analysis in which %IFBF was modeled using thyroid hormones were not congruent with results from earlier research [7]. For our data, %IFBF was best predicted ^ IFBF = 1.911+0.18146TT420.0026TT42+ by the model % 0.0000076TT43. However, the predictive ability of our model was quite low (R2 = 0.106). When the model from similar research ^ IFBF = 21.359+0.0756TT42 [7] was fit to our data % 2 0.0036TT4 20.0506FT4+0.0586chestgirth, R2 = 0.120), the model only received 2.6% of the model weight and had low predictive ability. Discussion The patterns reflected by our data tend to support our hypothesis that late winter body condition of adult female mule deer was elevated in our treatment study area as compared to our reference area. Both total fat reserves (%IFBF) and the metabolic Figure 2. Estimates of winter body fat on treated and un-treated study areas. Scaled estimates of late winter percent ingest free body fat (%IFBF), with 95% prediction intervals, for adult female mule deer in southwest Colorado. Solid gray bars reflect estimates for our treatment study area (Billy Creek State Wildlife Area) and white bars reflect estimates for our reference study area (Buckhorn Mountain). Estimates and prediction intervals were generated according to the model %IFBF~{2:159{0:534|Buckhornz0:905|Year2008 z0:539|Year2009 z0:092|Chest in which chest girth was held constant at the observed mean of 95.476 cm and coefficient estimates have been model averaged based on model results. doi:10.1371/journal.pone.0106374.g002 PLOS ONE | www.plosone.org 5 September 2014 | Volume 9 | Issue 9 | e106374 �Herbivore Body Condition and Habitat Management Table 3. Best predictive multiple linear regression models, based on Akaike’s Information Criteria (AICc), for 3 different body condition related parameters in adult female mule deer in southwest Colorado. Model vi a Kb R2 %IFBF = Area + Year + Chest 0.222 6 0.096 FT4 = Area + Year + PRG + Age 0.233 7 0.360 FT3 = Year + Age 0.141 5 0.269 a AICc model weight. Number of estimated parameters. Body condition parameters included % ingesta-free body fat (%IFBF), concentration of FT4 thyroid hormones (FT4), and concentrations of FT3 thyroid hormones (FT3). Data were modeled using study area (Area), year (Year), individual pregnancy status (PRG), chest girth (Chest), age (Age) and hind foot length (Foot). doi:10.1371/journal.pone.0106374.t003 b set of treatments (i.e., no spatial replication) with limited temporal replication, thereby introducing the potential for confounding factors. By spatially pairing our study units we attempted to control for issues associated with environmental stochasticity (i.e., extreme weather events, geological and vegetational attributes associated at micro-scales, migratory behavioral characteristics displayed by deer within a single herd). Despite these efforts, these confounding factors likely influenced our results. Similarly, we believe the wide variation in our data would have been better explained by past reproductive history. The burden, or lack thereof, stemming from lactation and energetic transfer from a dam to 1–2 offspring could not be estimated as part of this study. Despite design limitation, one of our key objectives was to evaluate the utility of body condition metrics in the context of actual habitat treatments, as opposed to artificial feeding utilized by other studies [7]. In general, the mechanical habitat treatments utilized as part of our research did mirror the pattern stemming from the pelleted food ration [7], but as expected, the magnitude of our treatment effect was lower and more tenuous when variation in estimates was considered. The research of Bishop et al. [7] reported %IFBF estimates of 10.21%–13.90% in treatment units and 6.64%–7.60% in control units, reflecting a ,1.616 magnitude increase. We detected a 1.086 magnitude increase using common habitat management techniques. Our results do yearly variation appeared in most of the best models and never carried less than 82% of the cumulative AICc weight. We suspect that had we been able to increase the positive effects of habitat treatments, the importance of yearly variation may have been diminished. However, the treatments delivered as part of our research reflect those commonly utilized by land management agencies. We also note that .90% of the variation within our %IFBF data remained unexplained. Much of this variation was likely due to individual characteristics (i.e., past reproductive success or failure, energetic burdens due to lactation, and habitat selection behaviors at micro-scales). Our study did not evaluate these important sources of variation. Thus, we conclude that while the effect of habitat management techniques are positive, thereby elevating the late winter body condition of mule deer, the magnitude of those effects are subtle and not be strong enough to eliminate the roles of yearly or individual variation. While the variation surrounding our estimates limited our ability to make a robust conclusion, this variation also serves as a road map for future research. Specifically, exploration of the components of this variation is warranted. Spatial and temporal aspects of this variation, largely regulated by annual moisture and weather patterns, need to be better understood. We echo the sentiments of others [5,36] that long-term, large-scale, individual-based studies are needed. We also recognize that our research evaluated a single Figure 3. Probability of pregnancy for different age classes of adult female mule deer, with 95% confidence intervals, for adult mule deer in southwest Colorado. No discernible difference in probability of pregnancy between our treatment and reference study areas was observed, although evidence for senescence in older age classes was observed. doi:10.1371/journal.pone.0106374.g003 PLOS ONE | www.plosone.org 6 September 2014 | Volume 9 | Issue 9 | e106374 �Herbivore Body Condition and Habitat Management not support the recommendations of Bishop et al. [13] that thyroid hormones could be used to estimate %IFBF as even our best predictive model did not attain a satisfactory level of performance and our overall correlation between %IFBF and hormones were low. The fact that our results do not fully validate these earlier results [13] is noteworthy as it advances our knowledge about the power and utility of body condition as a metric for assessing the impacts of environmental change. Population-level impacts stemming from the differences in body condition on our study areas were likely nominal. For example, we did not detect a meaningful difference in pregnancy rates between study areas. Likewise, while our study did not assess neonatal rates we do not think the number of fetuses produced per adult female in the treated area was greater than that in the reference area. However, we note that we did not actively seek low quality habitat to serve as our reference area. Rather, the reference area was defined by pinyon-juniper winter range that had not received vegetation treatments. This allowed us to test the hypothesis that habitat manipulation and improvement could be used to improve winter range in terms of late winter body condition. The magnitude of improvement in body condition could be expected to be amplified, relative to pre-treatment levels, if habitat treatments were targeted to poor quality habitat. In parallel research [37], the overwinter survival of mule deer fawns was measured on these and multiple other study areas during the same time period as our study. This research [37] documented a 1.156 increase in fawn survival on treated study areas, including BCSWA. When considered in tandem with fawn survival estimates [37], our results can be used to evaluate the sequence of density dependent effects experienced by mammals as their populations approach carrying capacity [1,2]. Predictions state that survival of young is the first population parameter to reflect a response under habitat limited scenarios [1,2]. This prediction is supported by research paralleling our study [37]. The second prediction is a delay in the onset of first pregnancy [1,2]. While evidence for senescence in older age classes was observed, we did not detect a delay in age of first pregnancy. Despite a truncated evaluation, based on our data the sequence of densitydependent effects hypothesized by earlier research [1,2] were likely correct. Thus, while our study provides a small step in linking environmental change with the fitness components of ungulates, it also exemplifies the need for further evaluation of the variation that is inherent within fitness components. Supporting Information Appendix S1 Data used in our analyses are available in Appendix S1. (PDF) Acknowledgments Delivery of mechanical treatments were largely coordinated by the Uncompahgre Plateau Project, a non-profit group composed of many agencies, organizations and individuals who saw the value in working cooperatively to improve wildlife habitat. We are indebted to A. Cline, C. Harty, B. Lamont, R. Lockwood, D. Lucchesi, J. McMillan, C. Santana, C. Tucker, and K. Yeager for their assistance with field work. Fixed wing pilots S. Waters and D. Felix provided assistance with aerial observation and spotting during capture efforts. Helicopter pilots R. Swisher and M. Shelton of Quicksilver Air, Inc. provided safe and efficient capture services. The time and energy provided by CPW personnel were instrumental to this research. Valuable insight, discussion, assistance and training were also provided by C. Anderson, R. Cook, D. Freddy, A. Holland, and B. Watkins. Useful feedback on earlier drafts of this manuscript were provided by M. Alldredge, A. Apa, and H. Johnson. Author Contributions Conceived and designed the experiments: EJB PFD CJB. Performed the experiments: EJB CJB LLW BAB. Analyzed the data: EJB PFD. Contributed reagents/materials/analysis tools: EJB CJB LLW. Contributed to the writing of the manuscript: EJB PFD. References 13. Bishop CJ, Watkins BE, Wolfe LL, Freddy DJ, White GC (2009) Evaluating mule deer body condition using serum thyroid hormone concentrations. Journal of Wildlife Management 73(3):462–467. 14. Feldman EC, Nelson RW (2004) Canine and feline endocrinology and reproduction. Saint Louis: Elsevier Press. 1089p. 15. Meyer DJ, Harvey JW (2004) Veterinary laboratory medicine: interpretation and diagnosis. Saint Louis: Elsevier Press. 351p. 16. Watkins BE, Bishop CJ, Bergman EJ, Hale B, Wakeling BF, et al. (2007) Habitat guidelines for mule deer: Colorado Plateau shrubland and forest ecoregion. Mule Deer Working Group. Available: http://www.muledeerworkinggroup. com/Docs/CPE_Mule_Deer_Habitat_Guidelines.pdf. Accessed 11 August 2014. 17. Young J A, Evans RA, Rimbey C (1985) Weed control and revegetation following western juniper (Juniperus occidentalis) control. Weed Science 33: 513–517. 18. Bates JD, Miller RF, Svejcar T (1998) Understory patterns in cut western juniper (Juniper occidentalis spp. occidentalis hook.) woodlands. Great Basin Naturalist 58: 363–374. 19. Bates JD, Miller RF, Svejcar T (2000) Understory dynamics in cut and uncut western juniper woodlands. Journal of Range Management 53: 119–126. 20. Miller RF, Svejcar TJ, Rose JA (2000) Impacts of western juniper on plant community composition and structure. Journal of Range Management 53: 574– 585. 21. Webb SL, Lewis JS, Hewitt DG, Hellickson MW, Bryant FC (2008) Assessing the helicopter and net gun as a capture technique for white-tailed deer. Journal of Wildlife Management 72: 310–314. 22. Jacques CN, Jenks JA, Deperno CS, Sievers JD, Grovenburg TW, et al. (2009) Evaluating ungulate mortality associated with helicopter net-gun captures in the Northern Great Plains. Journal of Wildlife Management 73: 1282–1291. 23. Severinghaus CW (1949) Tooth development and war as criteria of age in whitetailed deer. Journal of Wildlife Management 13: 195–216. 24. Robinette WL, Jones DA, Rogers G, Gashwiler JS (1957) Notes of tooth development and wear for Rocky Mountain mule deer. Journal of Wildlife Management 21: 134–153. 1. Eberhardt LL (1977) ‘‘Optimal’’ management policies for marine mammals. Wildlife Society Bulletin 5: 162–169. 2. Eberhardt LL (1977) Optimal policies for conservation of large mammals with special reference to marine ecosystems. Environmental Conservation 4: 205– 212. 3. Gaillard JM, Festa-Bianchet M, Yoccoz NG (1998) Population dynamics of large herbivores: variable recruitment with constant adult survival. Trends in Ecology and Evolution 13: 58–63. 4. Gaillard JM, Festa-Bianchet M, Yoccoz NG, Loison A, Toigo C (2000) Temporal variation in fitness components and population dynamics of large herbivores. Annual Review of Ecology and Systematics 31: 367–393. 5. Massey JS, Cubaynes, T . Coulson (2013) Will central Wyoming elk stop migrating to Yellowstone, and should we care? Ecology 94: 1271–1274. 6. Keech MA, Bowyer RT, Ver Hoef JM, Boertje RD, Dale BW, et al.(2000) Lifehistory consequences of maternal condition in Alaskan moose. Journal of Wildlife Management 64: 450–462. 7. Bishop C.J, White GC, Freddy DJ, Watkins BE, Stephenson TR (2009) Effect of enhanced nutrition on mule deer population rate of change. Wildlife Monographs 172, 29p. 8. Stephenson TR, Hundertmark KJ, Schwartz CC, Van Ballenberge V (1998) Predicting body fat and body mass in moose with ultrasonography. Canadian Journal of Zoology 76: 717–722. 9. Cook RC, Cook JG, Murray DL, Zager P, Johnson BK, et al. (2001) Development of predictive models of nutritional condition for Rocky Mountain elk. Journal of Wildlife Management 65: 973–987. 10. Cook RC, Cook JG, Stephenson TR, Myers WL, 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. 11. Watkins BE, Nachreiner RF, Schmitt SM, Cooley TM, Ullrey DE (1982) Thyroid function in fed and fasted white-tailed deer fawns. Canadian Journal of Zoology 60: 1331–1338. 12. Watkins BE, Ullrey DE (1983) Thyroid iodine and serum thyroid-hormone levels in wild white-tailed deer (Odocoileus-virginianus) from central Michigan. Canadian Journal of Zoology. 61: 1116–1119. PLOS ONE | www.plosone.org 7 September 2014 | Volume 9 | Issue 9 | e106374 �Herbivore Body Condition and Habitat Management 31. Vahtiala SH, Dahl E, Eloranta E, Beckers JF, Ropstad E (2004) Ultrasonography in early pregnancy diagnosis and measurements of fetal size in reindeer (Rangifer tarandus tarandus). Theriogenology 61: 785–795. 32. Wood AK, Short RE, Darling A, Dusek GL, Sasser RG, et al. (1986) Serum assays for detecting pregnancy in mule and white-tailed deer. Journal of Wildlife Management 50: 684–687. 33. Burnham KP, Anderson DR (2002) Model Selection and multimodel inference: a practical information-theoretical approach. New York: Springer Press. 488p. 34. Doherty PF, White GC, Burnham KP (2012) Comparison of model building and selection strategies. Journal of Ornithology 152: S317–S323. 35. Barbieri MM, Berger JO (2004) Optimal predictive model selection. Annals of Statistics. 32: 870–897. 36. Gaillard JM (2013) Assessing fitness consequences of migratory tactics requires long-term individually based monitoring. Ecology 94: 1261–1264. 37. Bergman EJ, Bishop CJ, Freddy DJ, White GC, Doherty PF (2014) Habitat management influences overwinter survival of mule deer fawns in Colorado. Journal of Wildlife Management 78: 448–455. 25. Hamlin KL, Pac DF, Sime CA, DeSimone RM, Dusek GL (2000) Evaluating the accuracy of ages obtained by two methods for Montana ungulates. Journal of Wildlife Management 64: 441–449. 26. Stephenson TR, Hundertmark KJ, Schwartz CC, Van Ballenberge V (1993) Ultrasonic fat measurement of captive yearling bull moose. Alces 29: 115–123. 27. 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. 28. Cook RC, Stephenson TR, Myers WL, Cook JG, Shipley LA (2007) Validating predictive models of nutritional condition for mule deer. Journal of Wildlife Management 71: 1934–1943. 29. Smith RB, Lindzey FG (1982) Use of ultrasound for detecting pregnancy in mule deer. Journal of Wildlife Management 46: 1089–1092. 30. Houston EW, Hagberg PK, Fisher MT, Miller ME, Asa CS (2001) Monitoring pregnancy in babirusa (Babyrousa babyrussa) with transabdominal ultrasonography. Journal of Zoo and Wildlife Medicine 32: 366–372. PLOS ONE | www.plosone.org 8 September 2014 | Volume 9 | Issue 9 | e106374 � 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 Herbivore body condition response in altered environments: mule deer and habitat management Description An account of the resource <span>The relationships between habitat, body condition, life history characteristics, and fitness components of ungulates are interwoven and of interest to researchers as they strive to understand the impacts of a changing environment. With the increased availability of portable ultrasound machines and the refinement of hormonal assays, assessment of ungulate body condition has become an accessible monitoring strategy. We employed body condition scoring, estimation of % ingesta-free body fat (%IFBF), assessment of free thyroid hormones (FT4 and FT3), and assessment of pregnancy, as metrics to determine if landscape-level habitat treatments affected body condition of adult (≥1.5 years old) female mule deer (</span><em>Odocoileus hemionus</em><span>). All body condition related metrics were measured on 2 neighboring study areas — a reference area that had received no habitat treatments and a treatment study area that had received mechanical removal of pinyon pine (</span><em>Pinyus edulis</em><span>) - Utah juniper (</span><em>Juniperus osteosperma</em><span>) forest, chemical control of weeds, and reseeding with preferred mule deer browse species. A consistent trend of higher %IFBF was observed in the treatment study area </span><span class="inline-formula"></span><span> than in the reference study area </span><span class="inline-formula"></span><span>, although variation of estimates was larger than hypothesized. A similar pattern was observed with higher thyroid hormones concentrations being observed in the treatment study area, but large amounts of variation within concentration estimates were also observed. The consistent pattern of higher body condition related estimates in our treatment study area provides evidence that large mammalian species are sensitive to landscape change, although variation within estimates underlie the challenge in detecting population level impacts stemming from environmental change.</span> Bibliographic Citation A bibliographic reference for the resource. Recommended practice is to include sufficient bibliographic detail to identify the resource as unambiguously as possible. Bergman, E. J., P. F. Doherty, C. J. Bishop, L. L. Wolfe, and B. A. Banulis. 2014. Herbivore body condition response in altered environments: mule deer and habitat management. PLoS One 9(9):e106374. <a href="https://doi.org/10.1371/journal.pone.0106374" target="_blank" rel="noreferrer noopener">https://doi.org/10.1371/journal.pone.0106374</a> Creator An entity primarily responsible for making the resource Bergman, Eric J. Doherty Jr, Paul F. Bishop, Chad J. Wolfe, Lisa L. Banulis, Bradley A. Subject The topic of the resource Mule deer Habitat management Extent The size or duration of the resource. 8 pages Date Created Date of creation of the resource. 2015-09-03 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> Format The file format, physical medium, or dimensions of the resource application/pdf 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 Type The nature or genre of the resource Article