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

Predicting future grizzly bear habitat use in the Bitterroot Ecosystem under recolonization and reintroduction scenarios

Many conservation actions must be implemented with limited data. This is especially true when planning recovery efforts for extirpated populations, such as grizzly bears (Ursus arctos) within the Bitterroot Ecosystem (BE), where strategies for reestablishing a resident population are being evaluated. Here, we applied individual-based movement models developed for a nearby grizzly bear population to predict habitat use in and near the BE, under scenarios of natural recolonization, reintroduction, and a combination. All simulations predicted that habitat use by grizzly bears would be higher in the northern half of the study area. Under the natural recolonization scenario, use was concentrated in Montana, but became more uniform across the northern BE in Idaho over time. Use was more concentrated in east-central Idaho under the reintroduction scenario. Assuming that natural recolonization continues even if bears are reintroduced, use remained widespread across the northern half of the BE and surrounding areas. Predicted habitat maps for the natural recolonization scenario aligned well with outlier and GPS collar data available for grizzly bears in the study area, with Spearman rank correlations of ≥0.93 and mean class values of ≥9.1 (where class 10 was the highest relative predicted use; each class 1-10 represented 10% of the landscape). In total, 52.4% of outlier locations and 79% of GPS collar locations were in class 10 in our predicted habitat maps for natural recolonization. Simulated grizzly bears selected habitats over a much larger landscape than the BE itself under all scenarios, including multiple-use and private lands, similar to existing populations that have expanded beyond recovery zones. This highlights the importance of recognizing and planning for the role of private lands in recovery efforts, including understanding resources needed to prevent and respond to human-grizzly bear conflict and maintain public acceptance of grizzly bears over a large landscape.

Long Arm of Motorway-The Impact of Fenced Road on the Mortality of European Badgers

We studied the impact of the new fenced and accident-safe motorway on the mortality of European badgers Meles meles on local roads in western Poland in 2010-2015. We monitored the badgers mortality on local roads of three categories: main roads, secondary roads and county roads. The study was conducted before and after the opening of the motorway in 2012. We hypothesized that the mortality of badgers is lower due to traffic concentration on motorway. Ninety two badgers were killed in collisions with vehicles on all monitored roads. Mean number of killed badgers was lowest in 2010 before the motorway opening and the highest in 2012. The mortality of badgers on regional roads was highest after the opening of the motorway due to the changes in traffic on the access roads. Within the road network, the mortality of badgers was 5.8 individuals/10 km of road per whole study period with the highest rate on main roads 8.5 individuals/10 km. The badgers mortality was highest on county roads but it was lower than expected in relation to the road network density. The highest vehicle collision risk for badgers of both sexes occurred in June. Distance to human settlements was the only environmental factor that was positively related to badger mortality on roads. We conclude that the new motorway did not reduce the mortality rate of badgers on the adjacent roads because the status of local roads has changed and now they mainly function as access roads to the motorway.

Working Together for Grizzly Bears: A Collaborative Approach to Estimate Population Abundance in Northwest Alberta, Canada

Grizzly bears are a threatened species in Alberta, Canada, and their conservation and management is guided by a provincial recovery plan. While empirical abundance and densities estimates have been completed for much of the province, empirical data are lacking for the northwest region of Alberta, a 2.8 million hectare area called Bear Management Area 1 (BMA 1). In part, this is due to limited staff capacity and funding to cover a vast geographic area, and a boreal landscape that is difficult to navigate. Using a collaborative approach, a multi-stakeholder working group called the Northwest Grizzly Bear Team (NGBT) was established to represent land use and grizzly bear interests across BMA 1. Collectively, we identified our project objectives using a Theory of Change approach, to articulate our interests and needs, and develop common ground to ultimately leverage human, social, financial and policy resources to implement the project. This included establishing 254 non-invasive genetic hair corral sampling sites across BMA 1, and using spatially explicit capture-recapture models to estimate grizzly bear density. Our results are two-fold: first we describe the process of developing and then operating within a collaborative, multi-stakeholder governance arrangement, and demonstrate how our approach was key to both improving relationships across stakeholders but also delivering on our grizzly bear project objectives; and, secondly we present the first-ever grizzly bear population estimate for BMA 1, including identifying 16 individual bears and estimating density at 0.70 grizzly bears/1,000 km2-the lowest recorded density of an established grizzly bear population in Alberta. Our results are not only necessary for taking action on one of Alberta's iconic species at risk, but also demonstrate the value and power of collaboration to achieve a conservation goal.

Landscape condition influences energetics, reproduction, and stress biomarkers in grizzly bears

Environmental change has been shown to influence mammalian distribution, habitat use, and behavior; however, few studies have investigated the impact on physiological function. This study aimed to determine the influence of landscape condition on the expression of target proteins related to energetics, reproduction, and stress in grizzly bears. We hypothesized that changes in landscape condition explains protein expression. Skin biopsies were collected from free-ranging grizzly bears in Alberta, Canada from 2013-2019 (n = 86 individuals). We used an information theoretic approach to develop 11 a priori candidate generalized linear mixed models to explain protein expression. We compared models using Akaike Information Criteria (AICc) weights and averaged models with ΔAICc < 2 for each protein. Food resources, represented by increased distance to coal mines and decreased crown closure, positively influenced energetic proteins (adiponectin and alpha-1-acid glycoprotein). Proteins related to reproduction (ceruloplasmin and serpin B5) were positively associated with increased wetland and upland food resources in addition to movement, but negatively associated with increased distance to roads. One stress related protein, complement C3, was positively influenced by increased percent conifer. Given the need to detect emerging threats to wildlife, we suggest the assessment of physiological function will lead to improved monitoring of species in rapidly changing landscapes.

Response of grizzly bears (Ursus arctos) to pipelines in Alberta

This research provides the first in-depth analysis of fine-scale grizzly bear habitat selection and movement patterns in response to the linear footprints cleared for below-ground pipelines in Alberta. Using an extensive set of GPS location data from collared grizzly bears, we were able to determine that grizzly bears selected for younger pipelines (mean age since last construction~6.5 years), which are known to have a greater abundance of important bear foods. Bears also selected for wider corridors that were disturbed for construction more than once. During the spring season, sex/age class was an important predictor of grizzly bear use of pipelines, with adult female bears more likely to use these features than other sex/age classes. Examination of movement patterns revealed that pipeline density influenced grizzly bears' movement rates and path straightness, particularly in the spring, when bears moved more slowly and movement paths were more tortuous in areas with higher pipeline densities. These movement patterns are consistent with foraging behavior and further indicate that bears are not exhibiting avoidance behaviors or displacement by pipeline features, and pipelines may be functioning as seasonally important foraging areas for grizzly bears in Alberta.

Road visibility influences habitat selection by grizzly bears (Ursus arctos horribilis)

Anthropogenic disturbances, including roads, are known to influence animal habitat selection and mortality. In this study, we consider the role of sensory perception in understanding why and how animals respond to disturbances. Our goal was to investigate the effect of visual perception (visibility) around roads on grizzly bear (Ursus arctos horribilis Ord, 1815) habitat selection and mortality in Alberta, Canada. We used detailed topographic and vegetation data from airborne light detection and ranging (lidar) to estimate visibility around roads. We modelled habitat selection as a function of road visibility and environmental variables using GPS telemetry data from 39 grizzly bears and integrated step selection analysis (iSSA). Finally, we assessed mortality risk in visible areas by comparing habitat selection between grizzly bears that died and those that survived. We found that grizzly bears were less likely to select visible areas when moving slowly or resting, but were more likely to select visible areas when travelling. We found that grizzly bears that survived selected for areas farther from roads than grizzly bears that died. However, no difference in selection for visible areas was observed. An exploratory analysis showed that grizzly bear mortalities commonly occurred in visible areas. Our findings highlight the importance of sensory perception in understanding animal behaviour.

Development and validation of protein biomarkers of health in grizzly bears

Large carnivores play critical roles in the maintenance and function of natural ecosystems; however, the populations of many of these species are in decline across the globe. Therefore, there is an urgent need to develop novel techniques that can be used as sensitive conservation tools to detect new threats to the health of individual animals well in advance of population-level effects. Our study aimed to determine the expression of proteins related to energetics, reproduction and stress in the skin of grizzly bears (Ursus arctos) using a liquid chromatography and multiple reaction monitoring mass spectrometry assay. We hypothesized that a suite of target proteins could be measured using this technique and that the expression of these proteins would be associated with biological (sex, age, sample location on body) and environmental (geographic area, season, sample year) variables. Small skin biopsies were collected from free-ranging grizzly bears in Alberta, Canada, from 2013 to 2019 (n = 136 samples from 111 individuals). Over 700 proteins were detected in the skin of grizzly bears, 19 of which were chosen as targets because of their established roles in physiological function. Generalized linear mixed model analysis was used for each target protein. Results indicate that sample year influenced the majority of proteins, suggesting that physiological changes may be driven in part by responses to changes in the environment. Season influenced the expression of proteins related to energetics, reproduction and stress, all of which were lower during fall compared to early spring. The expression of proteins related to energetics and stress varied by geographic area, while the majority of proteins that were affected by biological attributes (age class, sex and age class by sex interaction) were related to reproduction and stress. This study provides a novel method by which scientists and managers can further assess and monitor physiological function in wildlife.

Genetic tagging in the Anthropocene: scaling ecology from alleles to ecosystems

The Anthropocene is an era of marked human impact on the world. Quantifying these impacts has become central to understanding the dynamics of coupled human-natural systems, resource-dependent livelihoods, and biodiversity conservation. Ecologists are facing growing pressure to quantify the size, distribution, and trajectory of wild populations in a cost-effective and socially acceptable manner. Genetic tagging, combined with modern computational and genetic analyses, is an under-utilized tool to meet this demand, especially for wide-ranging, elusive, sensitive, and low-density species. Genetic tagging studies are now revealing unprecedented insight into the mechanisms that control the density, trajectory, connectivity, and patterns of human-wildlife interaction for populations over vast spatial extents. Here, we outline the application of, and ecological inferences from, new analytical techniques applied to genetically tagged individuals, contrast this approach with conventional methods, and describe how genetic tagging can be better applied to address outstanding questions in ecology. We provide example analyses using a long-term genetic tagging dataset of grizzly bears in the Canadian Rockies. The genetic tagging toolbox is a powerful and overlooked ensemble that ecologists and conservation biologists can leverage to generate evidence and meet the challenges of the Anthropocene.

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 31^st). 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.

EcoAnthromes of Alberta: An example of disturbance-informed ecological regionalization using remote sensing

Humans influence ecosystems on magnitudes that often exceed that of natural forces such as climate and geology; however, frameworks rarely include anthropogenic disturbance when delineating unique ecological regions. A critical step toward understanding, managing and monitoring human-altered ecosystems is to incorporate disturbance into ecological regionalizations. Furthermore, quantitative regionalization approaches are desirable to provide cost-effective, repeatable and statistically sound stratification for environmental monitoring. We applied a two-stage multivariate clustering technique to identify 'EcoAnthromes' across a large area - the province of Alberta, Canada - at 30 m spatial resolution, and using primarily remotely sensed inputs. The EcoAnthrome clusters represent regions with unique ecological characteristics based on a combination of natural ecological potential (e.g., climatic and edaphic factors) and disturbance, both natural and anthropogenic. Compared to existing expert-derived Natural Subregions in Alberta, the model-based EcoAnthromes showed greater class separation and explained more variance for an assortment of variables related to land cover, disturbance and species intactness. The EcoAnthromes successfully separated important ecological regions that are defined by complex assemblages of topography, climate and disturbance, such as gravel-bed river valleys, boreal forests, grasslands, post-fire recovery areas and highly disturbed agricultural, industrial and urban landscapes. In addition to presenting a flexible method for EcoAnthrome regionalization, we group and describe the EcoAnthromes created for Alberta and discuss how they can complement expert-derived regionalizations to aid in environmental management efforts, such as species recovery planning and monitoring for threatened species.

Identifying attributes associated with brown bear (Ursus arctos) road-crossing and roadkill sites

Habitat fragmentation caused by transportation infrastructure is an issue of growing concern worldwide. We show how secondary roads may affect landscape permeability for brown bears (Ursus arctos Linnaeus, 1758). We focused on identifying environmental variables that govern the selection of road-crossing zones by bears (crossing model). We also investigated whether variables that characterize road-crossing zones differ from those that are typical for bear–vehicle collision sites (collision model). The study area was located in north-central Slovakia. To identify road-crossing sites, we used the GPS fixes of 27 bears and identified 35 bear–vehicle collision sites from a different data set. We used mixed-effects logistic regression to model resource selection at road-crossing sites and to compare bear-crossing sites with bear-kill sites. The crossing model showed that the traffic volume with distance to forest and grassland were the most influential factors in bear selection of road-crossing sites. Results of the collision model indicated that successful road crossings by bears were located at different road sections from vehicle collisions, which differed by a traffic volume of 5000 vehicles/24 h. The outcomes of this study can facilitate improved mitigation measures on secondary roads.

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.

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.

Individuals matter: predicting koala road crossing behaviour in south-east Queensland

The ability to predict the frequency with which koalas (Phascolarctos cinereus) cross roads is an important step in developing mitigation strategies aimed at maintaining metapopulation viability for this species. In south-east Queensland, Australia, rapid urban development, including fragmentation and vehicle-related mortality, has contributed to a dramatic decline in local koala populations. Assessment of wildlife crossing structures that mitigates the impact of roads demonstrates these to be beneficial for many species. However, building enough structures to enable animals to move freely across impacted landscapes may be prohibitive in complex urban landscapes. The focus for this study was to consider the role of differing characteristics for species at risk. We examined the frequency of road crossings by individual koalas among six subpopulations in south-east Queensland. We explored the influence of potential predictor variables including: age, sex, and distance from roads. We established that road crossings were limited to a subset of individuals, with only 18 of 51 (~35%) koalas studied ever crossing a road. We found that koalas were disproportionately more likely to cross a road if they had initially been found within a distance of 100 m of that road. Males were more likely to cross roads than females, as were koalas aged less than five years. We suggest that successful wildlife road mitigation can be improved by better understanding the road-crossing behaviour of species such as the koala.

Forbidden fruit: human settlement and abundant fruit create an ecological trap for an apex omnivore

Habitat choice is an evolutionary product of animals experiencing increased fitness when preferentially occupying high-quality habitat. However, an ecological trap (ET) can occur when an animal is presented with novel conditions and the animal's assessment of habitat quality is poorly matched to its resulting fitness. We tested for an ET for grizzly (brown) bears using demographic and movement data collected in an area with rich food resources and concentrated human settlement. We derived measures of habitat attractiveness from occurrence models of bear food resources and estimated demographic parameters using DNA mark-recapture information collected over 8 years (2006-2013). We then paired this information with grizzly bear mortality records to investigate kill and movement rates. Our results demonstrate that a valley high in both berry resources and human density was more attractive than surrounding areas, and bears occupying this region faced 17% lower apparent survival. Despite lower fitness, we detected a net flow of bears into the ET, which contributed to a study-wide population decline. This work highlights the presence and pervasiveness of an ET for an apex omnivore that lacks the evolutionary cues, under human-induced rapid ecological change, to assess trade-offs between food resources and human-caused mortality, which results in maladaptive habitat selection.

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.

Effects of habitat quality and anthropogenic disturbance on grizzly bear (Ursus arctos horribilis) home-range fidelity

In the Rocky Mountain eastern slopes of Alberta, Canada, grizzly bears (Ursus arctos horribilis Ord, 1815) live in a landscape heavily impacted by industrial development and human disturbance. We characterized the role of changing habitat quality and new disturbance features on patterns of grizzly bear seasonal home-range fidelity and drift by comparing consecutive-year seasonal home ranges. We relied on the geographic technique “spatial–temporal analysis of moving polygons” (STAMP) to examine changes in habitat quality and new development between zones of home-range fidelity, expansion, and contraction. Areas considered to be high-quality habitat were selected at a greater frequency than available and retained in zones of home-range fidelity, but also vacated during home-range contraction. Areas of decreasing habitat quality were equally present in zones of contraction, expansion, and stability. The proportion of new forest harvest areas and roads developed within the past year did not differ between zones of home-range change, but the proportion of new well sites was higher in contraction zones than in stability zones. Our results showed that while considerable drift occurs, changes in habitat quality and recent anthropogenic disturbances cannot account for annual variation in home ranges, suggesting other important factors influencing behaviour and movement.