Cortisol levels in blood and hair of unanesthetized grizzly bears (Ursus arctos) following intravenous cosyntropin injection.

Hair cortisol concentration (HCC) is being used increasingly to evaluate long-term stress in many mammalian species. Most of the cortisol is assumed to passively diffuse from circulating blood into hair follicles and gradually accumulate in growing hair. However, our research with free-ranging grizzly bears (Ursus arctos) suggests HCC increases significantly within several hours following capture, a time too brief to be explained by this mechanism alone. In this study with captive grizzly bears, we sought to determine if a brief spike in blood cortisol concentration, thus mimicking a single stressful event, would cause an increase in HCC over a 7-day period. To do this, we administered a single intravenous dose (5 µg/kg) of cosyntropin to three captive unanaesthetised adult female grizzly bears on two occasions, during April when hair growth was arrested and during August when hair was growing. In both trials, the cosyntropin caused a two-fold or greater increase in serum cortisol levels within 1 hr but did not appear to influence HCC at 1, 48, and 168 hr following cosyntropin administration. We conclude the cosyntropin-induced cortisol spike was likely insignificant when compared to the adrenocortical response that occurs in free-ranging bears when captured. We suggest further study with a larger sample of captive bears to evaluate the combined effects of anaesthesia and multiple doses of cosyntropin administered over several hours would better simulate the adrenocortical response of free-ranging grizzly bears during capture.

Quantifying energetic costs and defining energy landscapes experienced by grizzly bears.

Animal movements are major determinants of energy expenditure and ultimately the cost-benefit of landscape use. Thus, we sought to understand those costs and how grizzly bears (Ursus arctos) move in mountainous landscapes. We trained captive grizzly bears to walk on a horizontal treadmill and up and down 10% and 20% slopes. The cost of moving upslope increased linearly with speed and slope angle, and this was more costly than moving horizontally. The cost of downslope travel at slower speeds was greater than the cost of traveling horizontally but appeared to decrease at higher speeds. The most efficient walking speed that minimized cost per unit distance was 1.19±0.11 m s-1 However, grizzly bears fitted with GPS collars in the Greater Yellowstone Ecosystem moved at an average velocity of 0.61±0.28 m s-1 and preferred to travel on near-horizontal slopes at twice their occurrence. When traveling uphill or downhill, grizzly bears chose paths across all slopes that were ∼54% less steep and costly than the maximum available slope. The net costs (J kg-1 m-1) of moving horizontally and uphill were the same for grizzly bears, humans and digitigrade carnivores, but those costs were 46% higher than movement costs for ungulates. These movement costs and characteristics of landscape use determined using captive and wild grizzly bears were used to understand the strategies that grizzly bears use for preying on large ungulates and the similarities in travel between people and grizzly bears that might affect the risk of encountering each other on shared landscapes.

Changing spring snow cover dynamics and early season forage availability affect the behavior of a large carnivore.

Changing climates are altering wildlife habitats and wildlife behavior in complex ways. Here, we examine how changing spring snow cover dynamics and early season forage availability are altering grizzly bear (Ursus arctos) behavior postden emergence. Telemetry data were used to identify spring activity dates for 48 individuals in the Yellowhead region of Alberta, Canada. Spring activity date was related to snow cover dynamics using a daily percent snow cover dataset. Snow melt end date, melt rate, and melt consistency explained 45% of the variation in spring activity date. We applied this activity date model across the entire Yellowhead region from 2000 to 2016 using simulated grizzly bear home ranges. Predicted spring activity date was then compared with a daily spring forage availability date dataset, resulting in "wait time" estimates for four key early season forage species. Temporal changes in both spring activity date and early season forage "wait times" were assessed using non-parametric regression. Grizzly bear activity date was found to have either remained constant (95%) or become earlier (5%) across the study area; virtually no areas with significantly later spring activity dates were detected. Similarly, the majority of "wait times" did not change (85%); however, the majority of significant changes in "wait times" for the four early season forage species indicated that "wait times" were lessening where changes were detected. Our results show that spring activity date is largely dictated by snow melt characteristics and that changing snow melt conditions may result in earlier spring activity. However, early season food stress conditions are likely to remain unchanged or improve as vegetation phenology also becomes earlier. Our findings extend the recent work examining animal movement in response to changing phenology from migratory birds and ungulates to an apex predator, further demonstrating the potential effects of changing climates on wildlife species.

Grizzly bear response to fine spatial and temporal scale spring snow cover in Western Alberta.

Snow dynamics influence seasonal behaviors of wildlife, such as denning patterns and habitat selection related to the availability of food resources. Under a changing climate, characteristics of the temporal and spatial patterns of snow are predicted to change, and as a result, there is a need to better understand how species interact with snow dynamics. This study examines grizzly bear (Ursus arctos) spring habitat selection and use across western Alberta, Canada. Made possible by newly available fine-scale snow cover data, this research tests a hypothesis that grizzly bears select for locations with less snow cover and areas where snow melts sooner during spring (den emergence to May 31st). Using Integrated Step Selection Analysis, a series of models were built to examine whether snow cover information such as fractional snow covered area and date of snow melt improved models constructed based on previous knowledge of grizzly bear selection during the spring. Comparing four different models fit to 62 individual bear-years, we found that the inclusion of fractional snow covered area improved model fit 60% of the time based on Akaike Information Criterion tallies. Probability of use was then used to evaluate grizzly bear habitat use in response to snow and environmental attributes, including fractional snow covered area, date since snow melt, elevation, and distance to road. Results indicate grizzly bears select for lower elevation, snow-free locations during spring, which has important implications for management of threatened grizzly bear populations in consideration of changing climatic conditions. This study is an example of how fine spatial and temporal scale remote sensing data can be used to improve our understanding of wildlife habitat selection and use in relation to key environmental attributes.

Natural regeneration on seismic lines influences movement behaviour of wolves and grizzly bears.

Across the boreal forest of Canada, habitat disturbance is the ultimate cause of caribou (Rangifer tarandus caribou) declines. Habitat restoration is a focus of caribou recovery efforts, with a goal to finding ways to reduce predator use of disturbances, and caribou-predator encounters. One of the most pervasive disturbances within caribou ranges in Alberta, Canada are seismic lines cleared for energy exploration. Seismic lines facilitate predator movement, and although vegetation on some seismic lines is regenerating, it remains unknown whether vegetation regrowth is sufficient to alter predator response. We used Light Detection and Ranging (LiDAR) data, and GPS locations, to understand how vegetation and other attributes of seismic lines influence movements of two predators, wolves (Canis lupus) and grizzly bears (Ursus arctos). During winter, wolves moved towards seismic lines regardless of vegetation height, while during spring wolves moved towards seismic lines with higher vegetation. During summer, wolves moved towards seismic lines with lower vegetation and also moved faster near seismic lines with vegetation <0.7 m. Seismic lines with lower vegetation height were preferred by grizzly bears during spring and summer, but there was no relationship between vegetation height and grizzly bear movement rates. These results suggest that wolves use seismic lines for travel during summer, but during winter wolf movements relative to seismic lines could be influenced by factors additional to movement efficiency; potentially enhanced access to areas frequented by ungulate prey. Grizzly bears may be using seismic lines for movement, but could also be using seismic lines as a source of vegetative food or ungulate prey. To reduce wolf movement rate, restoration could focus on seismic lines with vegetation <1 m in height. However our results revealed that seismic lines continue to influence wolf movement behaviour decades after they were built, and even at later stages of regeneration. Therefore it remains unknown at what stage of natural regeneration, if any, wolves cease to respond to seismic lines. To reduce wolf response to seismic lines, active restoration tactics like blocking seismic lines and tree planting, along with management of alternate prey, could be evaluated.

Using spatial mark-recapture for conservation monitoring of grizzly bear populations in Alberta.

One of the challenges in conservation is determining patterns and responses in population density and distribution as it relates to habitat and changes in anthropogenic activities. We applied spatially explicit capture recapture (SECR) methods, combined with density surface modelling from five grizzly bear (Ursus arctos) management areas (BMAs) in Alberta, Canada, to assess SECR methods and to explore factors influencing bear distribution. Here we used models of grizzly bear habitat and mortality risk to test local density associations using density surface modelling. Results demonstrated BMA-specific factors influenced density, as well as the effects of habitat and topography on detections and movements of bears. Estimates from SECR were similar to those from closed population models and telemetry data, but with similar or higher levels of precision. Habitat was most associated with areas of higher bear density in the north, whereas mortality risk was most associated (negatively) with density of bears in the south. Comparisons of the distribution of mortality risk and habitat revealed differences by BMA that in turn influenced local abundance of bears. Combining SECR methods with density surface modelling increases the resolution of mark-recapture methods by directly inferring the effect of spatial factors on regulating local densities of animals.

Functional macronutritional generalism in a large omnivore, the brown bear.

We combine a recently developed framework for describing dietary generalism with compositional data analysis to examine patterns of omnivory in a large widely distributed mammal. Using the brown bear (Ursus arctos) as a model species, we collected and analyzed data from the literature to estimate the proportions of macronutrients (protein, carbohydrate, and lipid) in the diets of bear populations. Across their range, bears consumed a diversity of foods that resulted in annual population diets that varied in macronutrient proportions, suggesting a wide fundamental macronutrient niche. The variance matrix of pairwise macronutrient log-ratios indicated that the most variable macronutrient among diets was carbohydrate, while protein and lipid were more proportional or codependent (i.e., relatively more constant log-ratios). Populations that consumed anthropogenic foods, such agricultural crops and supplementary feed (e.g., corn), had a higher geometric mean proportion of carbohydrate, and lower proportion of protein, in annual diets. Seasonally, mean diets were lower in protein and higher in carbohydrate, during autumn compared to spring. Populations with anthropogenic subsidies, however, had higher mean proportions of carbohydrate and lower protein, across seasons compared to populations with natural diets. Proportions of macronutrients similar to those selected in experiments by captive brown bears, and which optimized primarily fat mass gain, were observed among hyperphagic prehibernation autumn diets. However, the majority of these were from populations consuming anthropogenic foods, while diets of natural populations were more variable and typically higher in protein. Some anthropogenic diets were close to the proportions selected by captive bears during summer. Our results suggest that omnivory in brown bears is a functional adaptation enabling them to occupy a diverse range of habitats and tolerate variation in the nutritional composition and availability of food resources. Furthermore, we show that populations consuming human-sourced foods have different dietary macronutrient proportions relative to populations with natural diets.

Can concentrations of steroid hormones in brown bear hair reveal age class?

Although combining genetic and endocrine data from non-invasively collected hair samples has potential to improve the conservation of threatened mammals, few studies have evaluated this opportunity. In this study, we determined if steroid hormone (testosterone, progesterone, estradiol and cortisol) concentration profiles in 169 hair samples collected from free-ranging brown bears (Ursus arctos) could be used to accurately discriminate between immature and adult bears within each sex. Because hair samples were acquired opportunistically, we also needed to establish if interactions between hormones and several non-hormone factors (ordinal day, year, contact method, study area) were associated with age class. For each sex, we first compared a suite of candidate models by Akaike Information Criteria model selection, using different adult-age thresholds (3, 4 and 5 years), to determine the most supported adult age. Because hair hormone levels better reflect the endocrine state at an earlier time, possibly during the previous year, then at the time of sampling, we re-analysed the data, excluding the records for bears at the adult-age threshold, to establish if classification accuracy improved. For both sexes, candidate models were most supported based on a 3-year-old adult-age threshold. Classification accuracy did not improve with the 3-year-old bear data excluded. Male age class was predicted with a high degree of accuracy (88.4%) based on the concomitant concentrations of all four hormones. Female age class was predicted with less accuracy (77.1%) based only on testosterone and cortisol. Accuracy was reduced for females, primarily because we had poor success in correctly classifying immature bears (60%) whereas classification success for adult females was similar to that for males (84.5%). Given the small and unbalanced sample used in this study, our findings should be viewed as preliminary, but they should also provide a basis for more comprehensive future studies.

Space-time clusters for early detection of grizzly bear predation.

Accurate detection and classification of predation events is important to determine predation and consumption rates by predators. However, obtaining this information for large predators is constrained by the speed at which carcasses disappear and the cost of field data collection. To accurately detect predation events, researchers have used GPS collar technology combined with targeted site visits. However, kill sites are often investigated well after the predation event due to limited data retrieval options on GPS collars (VHF or UHF downloading) and to ensure crew safety when working with large predators. This can lead to missing information from small-prey (including young ungulates) kill sites due to scavenging and general site deterioration (e.g., vegetation growth). We used a space-time permutation scan statistic (STPSS) clustering method (SaTScan) to detect predation events of grizzly bears (Ursus arctos) fitted with satellite transmitting GPS collars. We used generalized linear mixed models to verify predation events and the size of carcasses using spatiotemporal characteristics as predictors. STPSS uses a probability model to compare expected cluster size (space and time) with the observed size. We applied this method retrospectively to data from 2006 to 2007 to compare our method to random GPS site selection. In 2013-2014, we applied our detection method to visit sites one week after their occupation. Both datasets were collected in the same study area. Our approach detected 23 of 27 predation sites verified by visiting 464 random grizzly bear locations in 2006-2007, 187 of which were within space-time clusters and 277 outside. Predation site detection increased by 2.75 times (54 predation events of 335 visited clusters) using 2013-2014 data. Our GLMMs showed that cluster size and duration predicted predation events and carcass size with high sensitivity (0.72 and 0.94, respectively). Coupling GPS satellite technology with clusters using a program based on space-time probability models allows for prompt visits to predation sites. This enables accurate identification of the carcass size and increases fieldwork efficiency in predation studies.

Recognizing the importance of an all-inclusive approach to brown bear conservation now and into the future.

In their critique of our recent article in Oecologia (Pigeon et al. Oecologia 181:1101-1116, 2016a) investigating the influence of ambient temperature on brown bear habitat selection, Ordiz et al. (2017, current issue) argue that we downplay the role of human disturbance on bear behavior, and that we wrongly report on the findings of Ordiz et al. (Oecologia 166:59-67, 2011). We argue that our previous article in Oecologia (Pigeon et al. 2016a) by no means downplays the influence of human factors on bear behavior, and that we correctly stated that Ordiz et al. (2011) did not adequately consider the potential influence of temperature on their findings. Finally, we stress the relevance of considering all-inclusive approaches to the common goal of successful wildlife conservation.

The quantification of reproductive hormones in the hair of captive adult brown bears and their application as indicators of sex and reproductive state.

Recognizing the potential value of steroid hormone measurements to augment non-invasive genetic sampling, we developed procedures based on enzyme-linked immunoassays to quantify reproductive steroid hormone concentrations in brown bear (Ursus arctos) hair. Then, using 94 hair samples collected from eight captive adult bears over a 2-year period, we evaluated (i) associations between hair concentrations of testosterone, progesterone, estradiol and cortisol; (ii) the effect of collecting by shaving vs. plucking; and (iii) the utility of reproductive hormone profiles to differentiate sex and reproductive state. Sample requirements (125 mg of guard hair) to assay all hormones exceeded amounts typically obtained by non-invasive sampling. Thus, broad application of this approach will require modification of non-invasive techniques to collect larger samples, use of mixed (guard and undercoat) hair samples and/or application of more sensitive laboratory procedures. Concentrations of hormones were highly correlated suggesting their sequestration in hair reflects underlying physiological processes. Marked changes in hair hormone levels during the quiescent phase of the hair cycle, coupled with the finding that progesterone concentrations, and their association with testosterone levels, differed markedly between plucked and shaved hair samples, suggests steroids sequestered in hair were likely derived from various sources, including skin. Changes in hair hormone concentrations over time, and in conjunction with key reproductive events, were similar to what has been reported concerning hormonal changes in the blood serum of brown bears. Thus, potential for the measurement of hair reproductive hormone levels to augment non-invasive genetic sampling appears compelling. Nonetheless, we are conducting additional validation studies on hair collected from free-ranging bears, representative of all sex, age and reproductive classes, to fully evaluate the utility of this approach for brown bear conservation and research.

Comparison of methanol and isopropanol as wash solvents for determination of hair cortisol concentration in grizzly bears and polar bears.

Methodological differences among laboratories are recognized as significant sources of variation in quantification of hair cortisol concentration (HCC). An important step in processing hair, particularly when collected from wildlife, is the choice of solvent used to remove or "wash" external hair shaft cortisol prior to quantification of HCC. The present study systematically compared methanol and isopropanol as wash solvents for their efficiency at removing external cortisol without extracting internal hair shaft cortisol in samples collected from free-ranging grizzly bears and polar bears. Cortisol concentrations in solvents and hair were determined in each of one to eight washes of hair with each solvent independently. •There were no significant decreases in internal hair shaft cortisol among all eight washes for either solvent, although methanol removed detectable hair surface cortisol after one wash in grizzly bear hair whereas hair surface cortisol was detected in all eight isopropanol washes.•There were no significant differences in polar bear HCC washed one to eight times with either solvent, but grizzly bear HCC was significantly greater in hair washed with isopropanol compared to methanol.•There were significant differences in HCC quantified using different commercial ELISA kits commonly used for HCC determinations.

A Double-Blinded, Randomized Comparison of Medetomidine-Tiletamine-Zolazepam and Dexmedetomidine-Tiletamine-Zolazepam Anesthesia in Free-Ranging Brown Bears (Ursus Arctos).

We compared anesthetic features, blood parameters, and physiological responses to either medetomidine-tiletamine-zolazepam or dexmedetomidine-tiletamine-zolazepam using a double-blinded, randomized experimental design during 40 anesthetic events of free-ranging brown bears (Ursus arctos) either captured by helicopter in Sweden or by culvert trap in Canada. Induction was smooth and predictable with both anesthetic protocols. Induction time, the need for supplemental drugs to sustain anesthesia, and capture-related stress were analyzed using generalized linear models, but anesthetic protocol did not differentially affect these variables. Arterial blood gases and acid-base status, and physiological responses were examined using linear mixed models. We documented acidemia (pH of arterial blood < 7.35), hypoxemia (partial pressure of arterial oxygen < 80 mmHg), and hypercapnia (partial pressure of arterial carbon dioxide ≥ 45 mmHg) with both protocols. Arterial pH and oxygen partial pressure were similar between groups with the latter improving markedly after oxygen supplementation (p < 0.001). We documented dose-dependent effects of both anesthetic protocols on induction time and arterial oxygen partial pressure. The partial pressure of arterial carbon dioxide increased as respiratory rate increased with medetomidine-tiletamine-zolazepam, but not with dexmedetomidine-tiletamine-zolazepam, demonstrating a differential drug effect. Differences in heart rate, respiratory rate, and rectal temperature among bears could not be attributed to the anesthetic protocol. Heart rate increased with increasing rectal temperature (p < 0.001) and ordinal day of capture (p = 0.002). Respiratory rate was significantly higher in bears captured by helicopter in Sweden than in bears captured by culvert trap in Canada (p < 0.001). Rectal temperature significantly decreased over time (p ≤ 0.05). Overall, we did not find any benefit of using dexmedetomidine-tiletamine-zolazepam instead of medetomidine-tiletamine-zolazepam in the anesthesia of brown bears. Both drug combinations appeared to be safe and reliable for the anesthesia of free-ranging brown bears captured by helicopter or by culvert trap.

Development and application of an antibody-based protein microarray to assess physiological stress in grizzly bears (Ursus arctos).

A novel antibody-based protein microarray was developed that simultaneously determines expression of 31 stress-associated proteins in skin samples collected from free-ranging grizzly bears (Ursus arctos) in Alberta, Canada. The microarray determines proteins belonging to four broad functional categories associated with stress physiology: hypothalamic-pituitary-adrenal axis proteins, apoptosis/cell cycle proteins, cellular stress/proteotoxicity proteins and oxidative stress/inflammation proteins. Small skin samples (50-100 mg) were collected from captured bears using biopsy punches. Proteins were isolated and labelled with fluorescent dyes, with labelled protein homogenates loaded onto microarrays to hybridize with antibodies. Relative protein expression was determined by comparison with a pooled standard skin sample. The assay was sensitive, requiring 80 µg of protein per sample to be run in triplicate on the microarray. Intra-array and inter-array coefficients of variation for individual proteins were generally <10 and <15%, respectively. With one exception, there were no significant differences in protein expression among skin samples collected from the neck, forelimb, hindlimb and ear in a subsample of n = 4 bears. This suggests that remotely delivered biopsy darts could be used in future sampling. Using generalized linear mixed models, certain proteins within each functional category demonstrated altered expression with respect to differences in year, season, geographical sampling location within Alberta and bear biological parameters, suggesting that these general variables may influence expression of specific proteins in the microarray. Our goal is to apply the protein microarray as a conservation physiology tool that can detect, evaluate and monitor physiological stress in grizzly bears and other species at risk over time in response to environmental change.

Staying cool in a changing landscape: the influence of maximum daily ambient temperature on grizzly bear habitat selection.

To fulfill their needs, animals are constantly making trade-offs among limiting factors. Although there is growing evidence about the impact of ambient temperature on habitat selection in mammals, the role of environmental conditions and thermoregulation on apex predators is poorly understood. Our objective was to investigate the influence of ambient temperature on habitat selection patterns of grizzly bears in the managed landscape of Alberta, Canada. Grizzly bear habitat selection followed a daily and seasonal pattern that was influenced by ambient temperature, with adult males showing stronger responses than females to warm temperatures. Cutblocks aged 0-20 years provided an abundance of forage but were on average 6 °C warmer than mature conifer stands and 21- to 40-year-old cutblocks. When ambient temperatures increased, the relative change (odds ratio) in the probability of selection for 0- to 20-year-old cutblocks decreased during the hottest part of the day and increased during cooler periods, especially for males. Concurrently, the probability of selection for 21- to 40-year-old cutblocks increased on warmer days. Following plant phenology, the odds of selecting 0- to 20-year-old cutblocks also increased from early to late summer while the odds of selecting 21- to 40-year-old cutblocks decreased. Our results demonstrate that ambient temperatures, and therefore thermal requirements, play a significant role in habitat selection patterns and behaviour of grizzly bears. In a changing climate, large mammals may increasingly need to adjust spatial and temporal selection patterns in response to thermal constraints.

Large Omnivore Movements in Response to Surface Mining and Mine Reclamation.

Increasing global demands have resulted in widespread proliferation of resource extraction. Scientists are challenged to develop environmental mitigation strategies that meet societal expectations of resource supply, while achieving minimal disruption to sensitive "wilderness" species. We used GPS collar data from a 9-year study on grizzly bears (Ursus arctos) (n = 18) in Alberta, Canada to assess movements and associated space use during versus after mining. Grizzly bear home range overlap with mined areas was lower during active mining except for females with cubs, that also had shortest movements on active mines. However, both females with cubs and males made shorter steps when on/close to mines following mine closure and reclamation. Our results show differences in bear movement and space-use strategies, with individuals from a key population segment (females with cubs) appearing most adaptable to mining disturbance. Preserving patches of original habitat, reclaiming the landscape and minimizing the risk of direct human-induced mortality during and after development can help conserve bears and other wildlife on industrially modified landscapes.

The impact of roads on the demography of grizzly bears in Alberta.

One of the principal factors that have reduced grizzly bear populations has been the creation of human access into grizzly bear habitat by roads built for resource extraction. Past studies have documented mortality and distributional changes of bears relative to roads but none have attempted to estimate the direct demographic impact of roads in terms of both survival rates, reproductive rates, and the interaction of reproductive state of female bears with survival rate. We applied a combination of survival and reproductive models to estimate demographic parameters for threatened grizzly bear populations in Alberta. Instead of attempting to estimate mean trend we explored factors which caused biological and spatial variation in population trend. We found that sex and age class survival was related to road density with subadult bears being most vulnerable to road-based mortality. A multi-state reproduction model found that females accompanied by cubs of the year and/or yearling cubs had lower survival rates compared to females with two year olds or no cubs. A demographic model found strong spatial gradients in population trend based upon road density. Threshold road densities needed to ensure population stability were estimated to further refine targets for population recovery of grizzly bears in Alberta. Models that considered lowered survival of females with dependant offspring resulted in lower road density thresholds to ensure stable bear populations. Our results demonstrate likely spatial variation in population trend and provide an example how demographic analysis can be used to refine and direct conservation measures for threatened species.

Idiosyncratic responses of grizzly bear habitat to climate change based on projected food resource changes.

Climate change vulnerability assessments for species of conservation concern often use species distribution and ecological niche modeling to project changes in habitat. One of many assumptions of these approaches is that food web dependencies are consistent in time and environmental space. Species at higher trophic levels that rely on the availability of species at lower trophic levels as food may be sensitive to extinction cascades initiated by changes in the habitat of key food resources. Here we assess climate change vulnerability for Ursus arctos (grizzly bears) in the southern Canadian Rocky Mountains using projected changes to 17 of the most commonly consumed plant food items. We used presence-absence information from 7088 field plots to estimate ecological niches and to project changes in future distributions of each species. Model projections indicated idiosyncratic responses among food items. Many food items persisted or even increased, although several species were found to be vulnerable based on declines or geographic shifts in suitable habitat. These included Hedysarum alpinum (alpine sweet vetch), a critical spring and autumn root-digging resource when little else is available. Potential habitat loss was also identified for three fruiting species of lower importance to bears: Empetrum nigrum (crowberry), Vaccinium scoparium (grouseberry), and Fragaria virginiana (strawberry). A general trend towards uphill migration of bear foods may result in higher vulnerability to bear populations at low elevations, which are also those that are most likely to have human-bear conflict problems. Regardless, a wide diet breadth of grizzly bears, as well as wide environmental niches of most food items, make climate change a much lower threat to grizzly bears than other bear species such as polar bears and panda bears. We cannot exclude, however, future alterations in human behavior and land use resulting from climate change that may reduce survival rates.

Macronutrient optimization and seasonal diet mixing in a large omnivore, the grizzly bear: a geometric analysis.

Nutrient balance is a strong determinant of animal fitness and demography. It is therefore important to understand how the compositions of available foods relate to required balance of nutrients and habitat suitability for animals in the wild. These relationships are, however, complex, particularly for omnivores that often need to compose balanced diets by combining their intake from diverse nutritionally complementary foods. Here we apply geometric models to understand how the nutritional compositions of foods available to an omnivorous member of the order Carnivora, the grizzly bear (Ursus arctos L.), relate to optimal macronutrient intake, and assess the seasonal nutritional constraints on the study population in west-central Alberta, Canada. The models examined the proportion of macronutrients that bears could consume by mixing their diet from food available in each season, and assessed the extent to which bears could consume the ratio of protein to non-protein energy previously demonstrated using captive bears to optimize mass gain. We found that non-selective feeding on ungulate carcasses provided a non-optimal macronutrient balance with surplus protein relative to fat and carbohydrate, reflecting adaptation to an omnivorous lifestyle, and that optimization through feeding selectively on different tissues of ungulate carcasses is unlikely. Bears were, however, able to dilute protein intake to an optimal ratio by mixing their otherwise high-protein diet with carbohydrate-rich fruit. Some individual food items were close to optimally balanced in protein to non-protein energy (e.g. Hedysarum alpinum roots), which may help explain their dietary prevalence. Ants may be consumed particularly as a source of lipids. Overall, our analysis showed that most food available to bears in the study area were high in protein relative to lipid or carbohydrate, suggesting the lack of non-protein energy limits the fitness (e.g. body size and reproduction) and population density of grizzly bears in this ecosystem.

Linking genotype, ecotype, and phenotype in an intensively managed large carnivore.

Numerous factors influence fitness of free-ranging animals, yet often these are uncharacterized. We integrated GPS habitat use data and genetic profiling to determine their influence on fitness proxies (mass, length, and body condition) in a threatened population of grizzly bears (Ursus arctos) in Alberta, Canada. We detected distinct genetic and habitat use (ecotype) clusters, with individual cluster assignments, or genotype/ecotype, being correlated (Pearson r = 0.34, P < 0.01). Related individuals showed evidence of similar habitat use patterns, irrespective of geographic distance and sex. Fitness proxies were influenced by sex, age, and habitat use, and homozygosity had a positive effect on these proxies that could be indicative of outbreeding depression. We further documented over 300 translocations occurring in the province since the 1970s, often to areas with significantly different habitat. We argue this could be unintentionally causing the pattern of outbreeding, although the heterozygosity correlation may instead be explained by the energetic costs associated with larger body size. The observed patterns, together with the unprecedented human-mediated migrations, make understanding the link between genotype, ecotype, and phenotype and mechanisms behind the negative heterozygosity-fitness correlations critical for management and conservation of this species.