Nutritional ecology of elk during summer and autumn in the Pacific Northwest

White-tailed deer population dynamics in a multipredator landscape shaped by humans

Large terrestrial mammals increasingly rely on human-modified landscapes as anthropogenic footprints expand. Land management activities such as timber harvest, agriculture, and roads can influence prey population dynamics by altering forage resources and predation risk via changes in habitat, but these effects are not well understood in regions with diverse and changing predator guilds. In northeastern Washington state, USA, white-tailed deer (Odocoileus virginianus) are vulnerable to multiple carnivores, including recently returned gray wolves (Canis lupus), within a highly human-modified landscape. To understand the factors governing predator-prey dynamics in a human context, we radio-collared 280 white-tailed deer, 33 bobcats (Lynx rufus), 50 cougars (Puma concolor), 28 coyotes (C. latrans), and 14 wolves between 2016 and 2021. We first estimated deer vital rates and used a stage-structured matrix model to estimate their population growth rate. During the study, we observed a stable to declining deer population (lambda = 0.97, 95% confidence interval: 0.88, 1.05), with 74% of Monte Carlo simulations indicating population decrease and 26% of simulations indicating population increase. We then fit Cox proportional hazard models to evaluate how predator exposure, use of human-modified landscapes, and winter severity influenced deer survival and used these relationships to evaluate impacts on overall population growth. We found that the population growth rate was dually influenced by a negative direct effect of apex predators and a positive effect of timber harvest and agricultural areas. Cougars had a stronger effect on deer population dynamics than wolves, and mesopredators had little influence on the deer population growth rate. Areas of recent timber harvest had 55% more forage biomass than older forests, but horizontal visibility did not differ, suggesting that timber harvest did not influence predation risk. Although proximity to roads did not affect the overall population growth rate, vehicle collisions caused a substantial proportion of deer mortalities, and reducing these collisions could be a win-win for deer and humans. The influence of apex predators and forage indicates a dual limitation by top-down and bottom-up factors in this highly human-modified system, suggesting that a reduction in apex predators would intensify density-dependent regulation of the deer population owing to limited forage availability.

Predation risk drives long-term shifts in migratory behaviour and demography in a large herbivore population

Migration is an adaptive life-history strategy across taxa that helps individuals maximise fitness by obtaining forage and avoiding predation risk. The mechanisms driving migratory changes are poorly understood, and links between migratory behaviour, space use, and demographic consequences are rare. Here, we use a nearly 20-year record of individual-based monitoring of a large herbivore, elk (Cervus canadensis) to test hypotheses for changing patterns of migration in and adjacent to a large protected area in Banff National Park (BNP), Canada. We test whether bottom-up (forage quality) or top-down (predation risk) factors explained trends in (i) the proportion of individuals using 5 different migratory tactics, (ii) differences in survival rates of migratory tactics during migration and whilst on summer ranges, (iii) cause-specific mortality by wolves and grizzly bears, and (iv) population abundance. We found dramatic shifts in migration consistent with behavioural plasticity in individual choice of annual migratory routes. Shifts were inconsistent with exposure to the bottom-up benefits of migration. Instead, exposure to landscape gradients in predation risk caused by exploitation outside the protected area drove migratory shifts. Carnivore exploitation outside the protected area led to higher survival rates for female elk remaining resident or migrating outside the protected area. Cause-specific mortality aligned with exposure to predation risk along migratory routes and summer ranges. Wolf predation risk was higher on migratory routes than summer ranges of montane-migrant tactics, but wolf predation risk traded-off with heightened risk from grizzly bears on summer ranges. A novel eastern migrant tactic emerged following a large forest fire that enhanced forage in an area with lower predation risk outside of the protected area. The changes in migratory behaviour translated to population abundance, where abundance of the montane-migratory tactics declined over time. The presence of diverse migratory life histories maintained a higher total population abundance than would have been the case with only one migratory tactic in the population. Our study demonstrates the complex ways in which migratory populations change over time through behavioural plasticity and associated demographic consequences because of individuals balancing predation risk and forage trade-offs.

Density-dependent changes in elk resource selection over successional time scales following forest disturbance

There is an increasing need to understand how animals respond to modifications of their habitat following landscape-scale disturbances such as wildfire or timber harvest. Such disturbances can promote increased use by herbivores due to changes in plant community structure that improve forage conditions, but can also cause avoidance if other habitat functions provided by cover are substantially reduced or eliminated. Quantifying the total effects of these disturbances, however, is challenging because they may not fully be apparent unless observed at successional timescales. Further, the effects of disturbances that improve habitat quality may be density dependent, such that the benefits are (1) less valuable to high-density populations because the per-capita benefits are reduced when shared among more users or, alternatively, (2) more valuable to animals living in high densities because resources may be more depleted from the greater intraspecific competition. We used 30 years of telemetry data on elk occurring at two distinct population densities to quantify changes in space use at diel, monthly, and successional timescales following timber harvest. Elk selected logged areas at night only, with selection strongest during midsummer, and peak selection occurring 14 years post harvest, but persisting for 26-33 years. This pattern of increased selection at night following a reduction in overhead canopy cover is consistent with elk exploiting improved nutritional conditions for foraging. The magnitude of selection for logged areas was 73% higher for elk at low population density, consistent with predictions from the ideal free distribution. Yet elk avoided these same areas during daytime for up to 28 years post logging and instead selected untreated forest, suggesting a role for cover to meet other life history requirements. Our results demonstrate that while landscape-scale disturbances can lead to increased selection by large herbivores and suggest that the improvement in foraging conditions can persist over short-term successional timescales, the magnitude of the benefits may not be equal across population densities. Further, the enduring avoidance of logging treatments during the daytime indicates a need for structurally intact forests and suggests that a mosaic of forest patches of varying successional stages and structural completeness is likely to be the most beneficial to large herbivores.

Evaluating the summer landscapes of predation risk and forage quality for elk (Cervus canadensis)

The recovery of carnivore populations in North American has consequences for trophic interactions and population dynamics of prey. In addition to direct effects on prey populations through killing, predators can influence prey behavior by imposing the risk of predation. The mechanisms through which patterns of space use by predators are linked to behavioral response by prey and nonconsumptive effects on prey population dynamics are poorly understood. Our goal was to characterize population- and individual-level patterns of resource selection by elk (Cervus canadensis) in response to risk of wolves (Canis lupus) and mountain lions (Puma concolor) and evaluate potential nonconsumptive effects of these behavioral patterns. We tested the hypothesis that individual elk risk-avoidance behavior during summer would result in exposure to lower-quality forage and reduced body fat and pregnancy rates. First, we evaluated individuals' second-order and third-order resource selection with a used-available sampling design. At the population level, we found evidence for a positive relationship between second- and third-order selection and forage, and an interaction between forage quality and mountain lion risk such that the relative probability of use at low mountain lion risk increased with forage quality but decreased at high risk at both orders of selection. We found no evidence of a population-level trade-off between forage quality and wolf risk. However, we found substantial among-individual heterogeneity in resource selection patterns such that population-level patterns were potentially misleading. We found no evidence that the diversity of individual resource selection patterns varied predictably with available resources, or that patterns of individual risk-related resource selection translated into biologically meaningful changes in body fat or pregnancy rates. Our work highlights the importance of evaluating individual responses to predation risk and predator hunting technique when assessing responses to predators and suggests nonconsumptive effects are not operating at a population scale in this system.

When can cropping rate compensate for increased vigilance?

Herbivores use vigilance to reduce predation risk and interact socially, yet it imposes a foraging efficiency cost. As individuals spend more time with their head up being vigilant, time available to search for and ingest food decreases. We explored whether ungulates can strategically modify behaviours to compensate for vigilance costs via increased cropping rate when food searching time was near-zero and bite sizes were small. We compared the proportion of time individuals had their head up to their cropping rate (bites/observation length) in 271 observations of Roosevelt elk (Cervus elaphus roosevelti). Using a linear mixed-effect model, we estimated the head up–cropping rate relationship and found that elk cropping rate was constant across varying lengths of time spent with their head up, indicating no vigilance compensation occurred via increased cropping rate. We discuss settings when cropping rate compensation is expected and other behaviours that might mitigate vigilance costs.

Influences on food supply from elk abundance and precipitation early in the growing season

Large grazing mammals should negatively affect forage biomass of their food supply, but documentation is lacking in free ranging populations. Furthermore, complications from factors such as weather patterns and spatial heterogeneity might obscure grazing effects on the food supply. We examined influences of Roosevelt elk (Cervus canadensis roosevelti (Merriam, 1897)) abundance and precipitation on forage biomass at two spatial scales; meadows that contained most of the food supply, and sectors nested in meadows. Spatial heterogeneity in forage biomass might also decline with increasing elk abundance. Elk abundance was estimated from population counts and varied 3.9-fold across the 15 years of study in northwestern California, USA. Each January, early in the growing season, we estimated forage biomass in the 50-ha meadow complex used by the elk population. Measures of palatable forage cover and height were taken in 270 ¼ m² plots dispersed throughout sectors. These measurements were then related to dried forage biomass. At both spatial scales, elk abundance was inversely, and precipitation was positively related to forage biomass. At the sector scale, analysis of a linear mixed effect model indicated heterogeneity. In some sectors both predictors were related to forage biomass and in other sectors they were not. Heterogeneity was not from uneven elk grazing as elk grazed sectors in proportion to forage biomass. The varied elk abundance–forage biomass relationships across sectors indicated that spatial heterogeneity declined with increasing elk abundance. Detecting relationships between free ranging ungulate populations and biomass of their food supply is not straightforward.

Non-invasive monitoring of multiple wildlife health factors by fecal microbiome analysis

Fecal microbial biomarkers represent a less invasive alternative for acquiring information on wildlife populations than many traditional sampling methodologies. Our goal was to evaluate linkages between fecal microbiome communities in Rocky Mountain elk (Cervus canadensis) and four host factors including sex, age, population, and physical condition (body-fat). We paired a feature-selection algorithm with an LDA-classifier trained on elk differential bacterial abundance (16S-rRNA amplicon survey) to predict host health factors from 104 elk microbiomes across four elk populations. We validated the accuracy of the various classifier predictions with leave-one-out cross-validation using known measurements. We demonstrate that the elk fecal microbiome can predict the four host factors tested. Our results show that elk microbiomes respond to both the strong extrinsic factor of biogeography and simultaneously occurring, but more subtle, intrinsic forces of individual body-fat, sex, and age-class. Thus, we have developed and described herein a generalizable approach to disentangle microbiome responses attributed to multiple host factors of varying strength from the same bacterial sequence data set. Wildlife conservation and management presents many challenges, but we demonstrate that non-invasive microbiome surveys from scat samples can provide alternative options for wildlife population monitoring. We believe that, with further validation, this method could be broadly applicable in other species and potentially predict other measurements. Our study can help guide the future development of microbiome-based monitoring of wildlife populations and supports hypothetical expectations found in host-microbiome theory.

Video recording and vegetation classification elucidate sheep foraging ecology in species-rich grassland

Factors influencing grazing behavior in species-rich grasslands have been little studied. Methodologies have mostly had a primary focus on grasslands with lower floristic diversity.We test the hypothesis that grazing behavior is influenced by both animal and plant factors and investigate the relative importance of these factors, using a novel combination of video technology and vegetation classification to analyze bite and step rates.In a semi-natural, partially wooded grassland in northern Estonia, images of the vegetation being grazed and records of steps and bites were obtained from four video cameras, each mounted on the sternum of a sheep, during 41 animal-hours of observation over five days. Plant species lists for the immediate field of view were compiled. Images were partnered by direct observation of the nearest-neighbor relationships of the sheep. TWINSPAN, a standard vegetation classification technique allocating species lists to objectively defined classes by a principal components procedure, was applied to the species lists and 25 vegetation classes (15 open pasture and 10 woodland) were identified from the images.Taking bite and step rates as dependent variables, relative importance of animal factors (sheep identity), relative importance of day, and relative importance of plant factors (vegetation class) were investigated. The strongest effect on bite rates was of vegetation class. Sheep identity was less influential. When the data from woodland were excluded, sheep identity was more important than vegetation class as a source of variability in bite rate on open pasture.The original hypothesis is therefore supported, and we further propose that, at least with sheep in species-rich open pastures, animal factors will be more important in determining grazing behavior than plant factors. We predict quantifiable within-breed and between-breed differences, which could be exploited to optimize conservation grazing practices and contribute to the sustainability of extensive grazing systems.

Seasonal patterns in nutritional condition of caribou (Rangifer tarandus) in the southern Northwest Territories and northeastern British Columbia, Canada

Evaluating nutritional condition provides insights of nutritional influences on wildlife populations. We sampled three measures of condition — body fat, body mass, and loin thickness — of adult female caribou (Rangifer tarandus (Linnaeus, 1758)) in boreal settings in the Northwest Territories (NT), Canada, in December and March, 2016–2018, and in mountain and boreal settings in British Columbia (BC), Canada, in December and February, 2014–2015. We evaluated the effect of calf-rearing on condition in December, compared influences of summer–autumn versus winter on condition over winter, and developed an annual profile of nutritional condition with estimates from caribou dying in summer. Mean December body fat was 8.4% in females with calves and 11.4% in females without calves, demonstrating the influence of lactation on condition. Over winter, nutritional condition did not decline in northeastern BC and it declined slightly in NT: body fat by 0.55 percentage points, mass by 2.8 kg, and loin thickness did not change. Body fat peaked in December, changed little over winter, but declined to a minimum by early summer, temporally coinciding with elevated rates of adult female mortality. Consistent with those of other ungulate studies worldwide, our findings suggest a need to focus on nutritional limitations operating in late spring through early autumn.

Evaluating indices of nutritional condition for caribou (Rangifer tarandus): which are the most valuable and why?

Body composition studies are critical for evaluating the accuracy of nutritional condition indices for predicting body components. We evaluated >40 indices of nutritional condition for caribou (Rangifer tarandus (Linnaeus, 1758)) using 29 female caribou captured from three populations in Alaska (USA) that ranged in condition from 2.3% to 11.2% ingesta-free body fat (IFBF) and 6 captive female caribou that ranged in condition from 8.1% to 26.0% IFBF. Estimates of body fat, protein, and gross energy were regressed against each index of nutritional condition. Generally, indices with linear or slightly curvilinear relations to body fat and those based on multiple fat depots were the most accurate in predicting nutritional condition and the most useful over the full range of nutritional condition. A scaledLIVINDEX (a combination of subcutaneous fat thickness and a condition score), CONINDEX (a combination of kidney fat and marrow fat), and a subset of the Kistner score (pericardium and kidneys only) had the strongest relationship with body fat (r2 > 0.86) and were useful over the entire range of nutritional condition. If used properly and with adequate training, indices of nutritional condition can be a critical tool for understanding the severity and seasonality of nutritional limitations in wild caribou populations.

Detecting Resource Limitation in a Large Herbivore Population Is Enhanced With Measures of Nutritional Condition

Resource limitation at the population level is a function of forage quality and its abundance relative to its per capita availability, which in turn, determines nutritional condition of individuals. Effects of resource limitation on population dynamics in ungulates often occur through predictable and sequential changes in vital rates, which can enable assessments of how resource limitation influences population growth. We tested theoretical predictions of bottom-up (i.e., resource limitation) forcing on moose (Alces alces) through the lens of vital rates by quantifying the relative influence of intrinsic measures of nutritional condition and extrinsic measures of remotely sensed environmental data on demographic rates. We measured rates of pregnancy, parturition, juvenile, and adult survival for 82 adult females in a population where predators largely were absent. Life stage simulation analyses (LSAs) indicated that interannual fluctuations in adult survival contributed to most of the variability in λ. We then extended the LSA to estimate vital rates as a function of bottom-up covariates to evaluate their influence on λ. We detected weak signatures of effects from environmental covariates that were remotely sensed and spatially explicit to each seasonal range. Instead, nutritional condition strongly influenced rates of pregnancy, parturition, and overwinter survival of adults, clearly implicating resource limitation on λ. Our findings depart from the classic life-history paradigm of population dynamics in ungulates in that adult survival was highly variable and generated most of the variability in population growth rates. At the surface, lack of variation explained by environmental covariates may suggest weak evidence of resource limitation in the population, when nutritional condition actually underpinned most demographics. We suggest that variability in vital rates and effects of resource limitation may depend on context more than previously appreciated, and density dependence can obfuscate the relationships between remotely sensed data and demographic rates.

Investigating the Dynamics of Elk Population Size and Body Mass in a Seasonal Environment Using a Mechanistic Integral Projection Model

Environmentally mediated changes in body size often underlie population responses to environmental change, yet this is not a universal phenomenon. Understanding when phenotypic change underlies population responses to environmental change is important for obtaining insights and robust predictions of population dynamics in a changing world. We develop a dynamic integral projection model that mechanistically links environmental conditions to demographic rates and phenotypic traits (body size) via changes in resource availability and individual energetics. We apply the model to the northern Yellowstone elk population and explore population responses to changing patterns of seasonality, incorporating the interdependence of growth, demography, and density-dependent processes operating through population feedback on available resources. We found that small changes in body size distributions can have large impacts on population dynamics but need not cause population responses to environmental change. Environmental changes that altered demographic rates directly, via increasing or decreasing resource availability, led to large population impacts in the absence of substantial changes to body size distributions. In contrast, environmentally driven shifts in body size distributions could occur with little consequence for population dynamics when the effect of environmental change on resource availability was small and seasonally restricted and when strong density-dependent processes counteracted expected population responses. These findings highlight that a robust understanding of how associations between body size and demography influence population responses to environmental change will require knowledge of the shape of the relationship between phenotypic distributions and vital rates, the population status with regard to its carrying capacity, and importantly the nature of the environmentally driven change in body size and carrying capacity.

State-dependent foraging by caribou with different nutritional requirements

Foraging by animals is hypothesized to be state-dependent, that is, varying with physiological condition of individuals. State often is defined by energy reserves, but state also can reflect differences in nutritional requirements (e.g., for reproduction, lactation, growth, etc.). Testing hypotheses about state-dependent foraging in ungulates is difficult because fine-scale data needed to evaluate these hypotheses generally are lacking. To evaluate whether foraging by caribou (Rangifer tarandus) was state-dependent, we compared bite and intake rates, travel rates, dietary quality, forage selection, daily foraging time, and foraging strategies of caribou with three levels of nutritional requirements (lactating adults, nonlactating adults, subadults 1-2 years old). Only daily foraging times and daily nutrient intakes differed among nutritional classes of caribou. Lactating caribou foraged longer per day than nonlactating caribou-a difference that was greatest at the highest rates of intake, but which persisted even when intake was below requirements. Further, at sites where caribou achieved high rates of intake, caribou in each nutritional class continued foraging even after satisfying daily nutritional requirements, which was consistent with a foraging strategy to maximize energy intake. Foraging time by caribou was partially state-dependent, highlighting the importance of accounting for physiological state in studies of animal behavior. Fine-scale foraging behaviors may influence larger-scale behavioral strategies, with potential implications for conservation and management.

Varied diets, including broadleaved forage, are important for a large herbivore species inhabiting highly modified landscapes

Diet quality is an important determinant of animal survival and reproduction, and can be described as the combination of different food items ingested, and their nutritional composition. For large herbivores, human landscape modifications to vegetation can limit such diet-mixing opportunities. Here we use southern Sweden’s modified landscapes to assess winter diet mixtures (as an indicator of quality) and food availability as drivers of body mass (BM) variation in wild moose (Alces alces). We identify plant species found in the rumen of 323 moose harvested in Oct-Feb, and link variation in average calf BM among populations to diets and food availability. Our results show that variation in calf BM correlates with variation in diet composition, diversity, and food availability. A varied diet relatively rich in broadleaves was associated with higher calf BM than a less variable diet dominated by conifers. A diet high in shrubs and sugar/starch rich agricultural crops was associated with intermediate BM. The proportion of young production forest (0–15 yrs) in the landscape, an indicator of food availability, significantly accounted for variation in calf BM. Our findings emphasize the importance of not only diet composition and forage quantity, but also variability in the diets of large free-ranging herbivores.

Animal responses to environmental variation: physiological mechanisms in ecological models of performance in deer (Cervidae)

Context Proper assessment of the consequences of environmental variation on animals depends on our ability to predict how they will perform under different circumstances. This requires two kinds of information. We need to know which environmental factors influence animal performance and their mode of action, i.e. whether a given factor acts alone or through interaction with other factors, directly or indirectly, instantaneously or after a delay and so on. This essentially correlative process falls within the domain of ecology. We also need to know what determines the direction, amplitude and limits of animal responses to environmental variation and change. This essentially experimental process falls within the domain of physiology. Physiological mechanisms are frequently poorly integrated within the correlative framework of ecological models. This is evident where programmed responses are attributed to environmental forcing and where the effect of environmental factors is evaluated without reference to the physiological state and regulatory capacity of the animal on which they act. Aims Here we examine ways in which the impacts of external (environmental) stimuli and constraints on performance are moderated by the animals (deer) on which they impinge. Key results The analysis shows (1) how trade-offs in foraging behaviour, illustrated by the timing of activity under the threat of predation, are modulated by integration of short-term metabolic feedback and animal emotions that influence the motivation to feed, (2) how the influence of thermal and nutritional challenges on performance, illustrated by the effect of weather conditions during gestation on the body mass of reindeer (Rangifer tarandus) calves at weaning, depends on the metabolic state of the female at the time the challenge occurs and (3) how annual cycles of growth, appetite and reproduction in seasonal species of deer are governed by innate circannual timers, such that their responses to seasonal changes in food supply are anticipatory and governed by rheostatic systems that adjust homeostatic set- points, rather than being purely reactive. Conclusions Concepts like ‘maintenance’ and ‘energy balance’, which were originally derived from non-seasonal domestic ruminants, are unable to account for annual cycles in metabolic and nutritional status in seasonal deer. Contrasting seasonal phenotypes (fat and anoestrous in summer, lean and oestrous in winter) represent adaptive solutions to the predictable challenges presented by contrasting seasonal environments, not failure of homeostasis in one season and its success in another. Implications The analysis and interpretation of responses to environment in terms of interaction between the external stimuli and the internal systems that govern them offer a more comprehensive, multifaceted understanding of the influence of environmental variation on performance in deer and open lines of ecological enquiry defined by non-intuitive aspects of animal function.

Herbicides and herbivory interact to drive plant community and crop-tree establishment

Land management practices often directly alter vegetation structure and composition, but the degree to which ecological processes such as herbivory interact with management to influence biodiversity is less well understood. We hypothesized that large herbivores compound the effects of intensive forest management on early seral plant communities and plantation establishment (i.e., tree survival and growth), and the degree of such effects is dependent on the intensity of management practices. We established 225-m² wild-ungulate (deer and elk) exclosures, nested within a manipulated gradient of management intensity (no-herbicide control, light herbicide, moderate herbicide, and intensive herbicide treatments), replicated at the scale of whole harvest units (10-19 ha). Vegetation structure, composition, and crop-tree responses to herbivory varied across the gradient of herbicide application during the first two years of stand establishment, with herbivory effects most evident at light and moderate herbicide treatments. In the moderate herbicide treatment, which approximates management applied to >2.5 million hectares in the Pacific Northwest, United States, foraging by deer and elk resulted in simplified, low-cover plant communities more closely resembling the intensive herbicide treatment. Herbivory further suppressed the growth of competing vegetation in the light herbicide treatment, improving crop-tree survival, and providing early evidence of an ecosystem service. By changing community composition and vegetation structure, intensive forest management alters foraging selectivity and subsequent plant-herbivore interactions; initial shifts in early seral communities are likely to influence understory plant communities and tree growth in later stages of forest development.

Straight from the caribou’s (Rangifer tarandus) mouth: detailed observations of tame caribou reveal new insights into summer–autumn diets

High-quality habitats for caribou (Rangifer tarandus (L., 1758)) are associated primarily with lichens, but lichens alone fail to satisfy summer nutritional requirements. To evaluate the summer forage value of plant communities across northeastern British Columbia (BC), where populations of northern and boreal ecotypes of caribou are declining, we observed foraging by tame, female caribou. We compared diet composition with forage abundance to determine forage selection and to quantify forage availability. Deciduous shrubs, not lichens, largely dominated summer diets. Caribou were highly selective foragers, with 28 species comprising 78% of diets. Caribou avoided ≥50% of understory vegetation in all communities, especially conifers, evergreen shrubs, mosses, and two genera of terrestrial lichens. Availability of accepted forage (species not avoided) was strongly heterogeneous across landscapes. Alpine shrub areas and mid-elevation spruce–fir stands in the mountains, as well as treed rich fens and white spruce communities in the boreal forests, provided the greatest quantities of accepted forage for caribou. Dry alpine sites and unproductive black spruce communities provided the least accepted forage. Our work has direct implications to caribou conservation by contributing to a greater understanding of the forage value of summer habitats, with implications to habitat selection, seasonal movements, and distribution ecology.