Geology topics

Research about Glacier National Park

Source-linked reports with geographic coverage including Glacier National Park.

At least 19 recordsLinked to original sources

USGS2026

Climate-driven shifts in snowpack conditions favorable for wet snow avalanches: Insights from snow cover model simulations

Forecasting wet snow avalanches remains a major challenge due to their temporal and spatial variability. Long-term climate trends are producing more variable and extreme weather, creating conditions outside historical experience, and increasing forecasting uncertainties. To assess how climate change affects the frequency and timing of wet snow avalanches, we examined the long-term conditions that support wet snow formation. Specifically, we ask: Will the onset of wet snow avalanche conditions or season end shift toward the mid century? And will maximum snow height decrease in avalanche starting zones? To address these questions, we analyze long-term snow cover simulations across our study site in Glacier National Park, Montana, located in an intermountain snow climate. We examined trends in simulated snow depth, and timing of liquid water infiltration over the past four decades (1979–2023) using a physically based snow cover model. To force the model with meteorological data, we use the long-term regional hydroclimate reanalysis dataset U.S. Geological Survey CONUS404 with 4 km spatial resolution and 1h temporal resolution. The use of gridded meteorological inputs to drive a one-dimensional point model introduces inherent scale mismatches. To address this limitation, we use multiple grid points within the study area to simulate a more representative, generalized snowpack, and filter the grid points based on elevation to distinguish between upper- and mid-elevation starting zones. Finally, we analyzed shifts in snow depth, melt-out-day, and LWCindex (bulk liquid water content (LWC) normalized by the residual water content). We use snow depth as a proxy for changes in slab thickness, the LWCindexas an indicator of wet snow avalanche conditions, and melt-out date as a measure of shifts in season length. We ran simulations on flat terrain and on north- and south-facing virtual slopes (38°), across mid-elevation (1,600–1,900 m) and upper-elevation (≥1,900 m) sites for the historical dataset as baseline and a Pseudo-Global-Warming (PGW) dataset to identify a warming driven shift. Across all aspects and both elevation bands, PGW forcing produces a consistent reduction in annual maximum snow depth (median shift −0.15 to −0.26 m, p<0.001) over 44 seasons, with the largest reductions on north-facing slopes and a modest but significant strengthening of this reduction over the 44-year simulated record. Melt-out starts 10–14 days earlier under PGW across all aspect–elevation combinations (p<0.001), again most pronounced on north-facing, higher-elevation terrain. The LWCindex threshold-crossing dates likewise shift earlier under PGW, indicating an earlier onset of wet snow avalanche conditions across the domain. Together, these results point to a snowpack that is both thinner and wetter earlier in the season under future warming, with the strongest signal consistently on north-facing slopes at higher elevation. This work lays the groundwork for climate-informed forecasting and mitigation and offers historical context for avalanche operations in a changing climate. Understanding long-term trends can help avalanche professionals interpret current observations, improve wet snow avalanche forecasting, and develop more effective mitigation strategies in a warming climate.

Montana↗

USGS2026

Mountain goat declines in a protected, interior, native population

A shifting climate poses threats to alpine-adapted species including mountain goats. We used long-term (12 years) citizen science monitoring data and Bayesian N-mixture modeling to estimate population trends and drivers of population metrics among mountain goats in Glacier National Park (GNP). Median goats per site ( n = 37 sites) declined by 45% (95% credible interval [CRI] = 32%, 57%) from 77.8 (95% CRI = 64.4, 95.1) in 2008 to 42.3 (95% CRI = 34.3, 52.2) in 2019, with consistent declines from 2008 until 2015, when the number of estimated goats stabilized. The decline exceeds IUCN criteria for classifying a population as vulnerable, >30% declines over only two generations. Across years, relatively few goats occupied northwestern GNP. Goat numbers declined the most at northeastern sites, trended toward decline in most southern sites, and increased at only two west-central sites. The proportion of permanent snow and glaciers, the presence of natural mineral licks, and habituation strongly increased the initial abundance of goats in the area. Weather variables had the greatest influence on population growth rates, particularly precipitation between May 15 and June 15 of the previous summer, the neonatal period. Lower growth occurred with less snow water equivalent and lower mean winter temperature, early summer temperature, and early summer precipitation. Projected reductions of permanent snow, increasing spring and summer temperatures, and insufficient and variable spring precipitation raise concerns for the future of native goats in this region. Our analyses reveal ways to improve detection rates of goats during surveys, which is important for optimizing the precision of estimates and the power to detect future trends. Detection increased with goat habituation, retention of observers with experience, use of binoculars, and conducting surveys at lower temperatures and earlier dates. Improving detection will be particularly important given the lower number of goats currently observed in the park. Research to estimate park-wide population size, evaluate genetic structure and diversity, assess changing habitat, human recreation levels and forage, and forward-project climate effects on persistence will be crucial to understanding the context of these results and conserving this iconic, metapopulation at the southern edge of the distribution of native mountain goats.

Montana↗

USGS2025

Ecosystem-engineered infections: Beaver-modified wetlands are associated with conflicting drivers of amphibian pathogen prevalence

Beavers are ecosystem engineers and keystone species that protect freshwater resources and increase biodiversity. Beaver reintroductions are promoted for amphibian conservation, yet their impact on Batrachochytrium dendrobatidis (Bd), a pathogen linked with amphibian population declines worldwide, remains unclear. We investigated the abiotic and biotic drivers of Bd prevalence in Columbia spotted frogs ( Rana luteiventris ) and western toads ( Anaxyrus boreas ) in 20 beaver-modified and 23 non-beaver wetlands in Glacier National Park, USA. We found that beavers increased wetland hydroperiod, which was associated with higher Bd prevalence. However, beavers also reduced wetland canopy cover, which was associated with lower Bd prevalence. Our models also predicted higher Bd prevalence associated with higher adult density of both species of amphibians, although species’ densities were similar in beaver-modified and non-beaver wetlands. These results suggest that beavers have a cumulatively negligible net effect on Bd prevalence owing to their effects on both hydroperiod and canopy cover, which is encouraging for amphibian conservation. Our findings also suggest that decreasing canopy cover may be a potential management option to reduce Bd prevalence. In addition, these findings indicate that beaver-mimicking restoration projects may harm amphibian populations if they increase wetland hydroperiods without reducing canopy cover.

Montana↗

USGS2025

Documenting, quantifying, and modeling a large glide avalanche in Glacier National Park, Montana, USA

Glide avalanches present a significant and repetitive challenge to many operational forecasting programs, and they are likely to become more frequent. While the spatial location of glide release areas is extremely consistent, the onset of glide avalanche release is notoriously difficult to forecast, and their destructive potential can be immense. Thus, the timing and dynamics of glide avalanches is an important area of study. To better understand these processes, and to improve assessments of risk to transportation corridors and infrastructure, event documentation is key. Here, we survey a large glide avalanche event along the Going-to-the-Sun Road in Glacier National Park, Montana, USA, during road opening operations in the spring of 2022. Using three sets of terrestrial lidar data (pre-event, post-event, and snow-off), we quantified key aspects of the avalanche and created powerful visualizations for analysis. Further, we evaluated meteorological data from automated weather stations between the onset of glide cracking and avalanche release. Last, we synthesized lidar data with a numerical dynamics model to replicate the event in a simulated environment. Using the tuned model, we determined the critical mean snow depth in the release area necessary for an avalanche to reach the road (4.2 m). Our method may be of particular use for glide avalanches, which tend to release in roughly the same place and time each year at a known interface. This could make the calculated critical depths more consistently reliable and preclude the need for additional tuning in dynamics models. As 1) lidar technology continues to improve and reduce in cost, 2) transportation corridors continue to extend into avalanche terrain, and 3) glide avalanches potentially become increasingly frequent, the synthesis outlined here provides a valuable tool for operational forecasters considering infrastructure threatened by glide events.

Montana↗

USGS2025

Myiasis infection by the toad fly (Lucilia bufonivora; Calliphoidae) in amphibians in Montana, USA

Toad flies in the genus Lucilia (previously referred to as Bufolucilia spp.) parasitize and cause myiasis in several amphibian species in North America. From 2019 to 2022, we documented Lucilia bufonivora infections in post-metamorphic western toads (Anaxyrus boreas ) during amphibian surveys in four wetlands in Glacier National Park, Montana, US. We found nine infected adult toads in 2019, seven infected adults in 2020, one infected juvenile in 2021, and five infected adults plus one infected juvenile in 2022. We also captured Columbia spotted frogs ( Rana luteiventris ) during these same surveys but detected no infections. Only one of the four wetlands had infected toads in 2019, despite their proximity and hydrologic connectivity, but two of these wetlands had infections in 2020, and a third had a single infection in 2021. The same three of four wetlands had infections in 2022. In 2008, a similar parasitic infection in one western toad had been noted at the same wetland as in 2019. That toad had been captured again two years later without signs of infection.

Montana↗

USGS2025

Assessing diet and genotyping success of goat pellet surveys from 2019 in Glacier National Park

Fecal pellets contain genetic information and can be used to identify individuals, their diet, and more. Individual identification can be useful in understanding movements of individuals, developing population estimates, assessing vital rates, genetic diversity and structure, and evaluating trends over time (e.g., Epps et al 2024). Successful genotyping depends on the quality of the sample, which can be influenced by many things including the initial state of the sample, environmental conditions that the pellets were exposed to, and the method of storage. The diet of mountain goats ( Oreamnos americanus ) is hard to study as they are habitat specialists in some of the most remote and rugged alpine environments (Festa-Bianchet & Côté 2012). These iconic ungulates epitomize the renowned wilderness character of Glacier National Park (GNP), yet the population has declined by 45% since 2008 likely due at least in part to changing climate (Graves et al. 2025). This is well above the International Union for the Conservation of Nature’s population decline criterion used to designate a species as vulnerable (IUCN 2012). Mountain goats are known as dietary generalists employing highly variable diets across systems, but diets have not been identified for GNP. Identifying goat diet composition in GNP can provide the information needed to determine whether this population decline could be related to climate-change induced shifts in forage phenology, diversity, productivity, and plant community structure. Other anthropogenic changes in GNP, such as increases in exotic and invasive plant species, will continue to affect plant community composition (Lesica et al. 1993). These alterations may reduce optimal forage availability and quality, which could have a deleterious effect on mountain goat survival. To examine these possibilities, we must first ask- What precisely are the mountain goats of Glacier National Park eating? DNA metabarcoding of fecal pellets provides a cutting-edge, non-invasive method for assessing herbivore diet. Previous studies of mountain goat diet have used methods such as fecal crude protein analysis (Festa-Bianchet & Côté 2012), rumen analysis (Saunders 1955), fecal microhistology (Cobb et al. 2012), and bite-for-bite feeding observations (Daily et al. 1984). Many of these studies have only been able to identify diet components to functional plant type and the methods used are biased towards less digestible diet components such as grasses (McInnis et al. 1983). Metabarcoding of fecal pellets can provide a higher taxonomic resolution than previous methods (Scasta et al. 2019, Stapleton et al. 2022) and open the field for increased participation of volunteers, also called community or citizen scientists. Here, we use mountain goat fecal pellets to conduct 1) an evaluation of which conditions might influence genotyping success and 2) a preliminary assessment of diet using metabarcoding. This pilot study allowed us to evaluate which pellets to focus on collecting, identify the best storage methods, explore approaches for analyzing metabarcoding data, examine the utility of metabarcoding for acquiring taxonomically specific diet information, and to identify genetic sequences that could not be resolved to genus or species levels using these methods. We examined the frequency of occurrence and relative diet composition of identified forbs, shrubs, graminoids, mosses, and trees in mountain goat diets.

Montana↗

USGS2024

Evaluating density-weighted connectivity of black bears (Ursus americanus) in Glacier National Park with spatial capture–recapture models

Background Improved understanding of wildlife population connectivity among protected area networks can support effective planning for the persistence of wildlife populations in the face of land use and climate change. Common approaches to estimating connectivity often rely on small samples of individuals without considering the spatial structure of populations, leading to limited understanding of how individual movement links to demography and population connectivity. Recently developed spatial capture-recapture (SCR) models provide a framework to formally connect inference about individual movement, connectivity, and population density, but few studies have applied this approach to empirical data to support connectivity planning. Methods We used mark-recapture data collected from 924 genetic detections of 598 American black bears ( Ursus americanus ) in 2004 with SCR ecological distance models to simultaneously estimate density, landscape resistance to movement, and population connectivity in Glacier National Park northwest Montana, USA. We estimated density and movement parameters separately for males and females and used model estimates to calculate predicted density-weighted connectivity surfaces. Results Model results indicated that landscape structure influences black bear density and space use in Glacier. The mean density estimate was 16.08 bears/100 km 2 (95% CI 12.52–20.6) for females and 9.27 bears/100 km 2 (95% CI 7.70–11.14) for males. Density increased with forest cover for both sexes. For male black bears, density decreased at higher grizzly bear ( Ursus arctos ) densities. Drainages, valley bottoms, and riparian vegetation decreased estimates of landscape resistance to movement for male and female bears. For males, forest cover also decreased estimated resistance to movement, but a transportation corridor bisecting the study area strongly increased resistance to movement presenting a barrier to connectivity. Conclusions Density-weighed connectivity surfaces highlighted areas important for population connectivity that were distinct from areas with high potential connectivity. For black bears in Glacier and surrounding landscapes, consideration of both vegetation and valley topography could inform the placement of underpasses along the transportation corridor in areas characterized by both high population density and potential connectivity. Our study demonstrates that the SCR ecological distance model can provide biologically realistic, spatially explicit predictions to support movement connectivity planning across large landscapes.

Montana↗

USGS2024

Characterizing vegetation and return periods in avalanche paths using lidar and aerial imagery

Snow avalanches are a hazard and ecological disturbance across mountain landscapes worldwide. Understanding how avalanche frequency affects forests and vegetation improves infrastructure planning, risk management, and avalanche forecasting. We implemented a novel approach using lidar, aerial imagery, and a random forest model to classify imagery-observed vegetation within avalanche paths in southern Glacier National Park, Montana, USA. We calculated spatially explicit avalanche return periods using a physically based spatial interpolation method and characterized the vegetation within those return period zones. The automated vegetation classification model differed slightly between avalanche paths, but the combination of lidar and spectral signature metrics provided the best accuracy (88–92 percent) for predicting vegetation classes within complex avalanche terrain rather than lidar or spectral signature metrics alone. The highest frequency avalanche return periods were broadly characterized by grassland and shrubland, but the influence of topography greatly influences the vegetation classes as well as the return periods. Furthermore, statistically significant differences in lidar-derived vegetation canopy height exist between categorical return periods. The ability to characterize vegetation within various avalanche return periods using remote sensing data provides land use planners and avalanche forecasters a tool for assessing the spatial extent of large-magnitude avalanches in individual avalanche paths.

Montana↗

USGS2024

Glacier National Park bumble bee survey report 2023

Glacier National Park's grasslands provide important contributions to the character and ecology of the park. In 1999-2001, Glacier National Park (hereafter, the park) established 155 permanent vegetation monitoring plots to inventory grassland vegetation communities east of the Continental Divide. In June of 2023, the Blackfeet Nation (Amskapi Piikuni) re-introduced a herd of free roaming bison on adjacent tribal land (the linnii Initiative), and those bison may enter the park. Because grasslands are important areas for bison grazing, and because biological communities have likely changed since the initial inventories, the park’s vegetation program began revisiting grassland monitoring plots in 2018. New data will provide insights into how these grasslands, and the wildlife that depend on them, have changed over the past two decades. Data on pollinators in these grasslands, however, are limited. Bison have been shown to influence plant communities, for example, by increasing herbaceous forb diversity and drought resilience (Elson and Hartnett 2017; Ratajczak et al. 2022). These changes directly influence resources available to foraging pollinators. To provide understanding of pollinators prior to bison reintroduction, bumble bee surveying at these grassland plots began in 2020 and has continued through 2023. Target sites for 2020-2022 focused on collecting data at plots where new vegetation surveys were recently completed to match current vegetation data with the bumble bee data.

Montana↗

USGS2023

Mapping a glide avalanche with terrestrial lidar in Glacier National Park, USA

Thorough documentation of large avalanche events is important for forecasting efforts, infrastructure planning, and investigating the processes involved in avalanche formation and release. However, due in part to the isolated and dangerous nature of avalanche terrain, collecting in-situ, spatially continuous, and quantitative information surrounding avalanches remains difficult. Advances in remote sensing continue to address this knowledge gap. For example, terrestrial laser scanners (TLSs) can produce snow depth measurements at fine spatial resolutions over large areas. Repeat data acquisitions between precipitation events also allow for depth quantification atop an interface, as well as precise estimations of release volume and runout area after avalanche failure. Here, we explore the benefits of TLS-derived documentation from a large avalanche event by examining the development and release of a glide avalanche that occurred in Glacier National Park, Montana, USA, during the spring of 2022. Three sets of lidar point cloud data were acquired in the Haystack Creek drainage, focused on a well-known glide avalanche site. Lidar scans were collected after glide cracks emerged but prior to glide failure, and shortly (~ 1.5 days) after avalanche occurrence, in addition to a snow-free scan later in the year. With this temporal dataset, we were able to account for and visualize the spatial variability of snow depth across the avalanche start zone, such that we could precisely calculate the release volume (18674 m3) and average start zone depth (3.3 m) of the avalanche. Furthermore, TLS data were used to map the extent of the runout area and entrainment zone.

Montana↗

USGS2020

Specialized meltwater biodiversity persists despite widespread deglaciation

Glaciers are important drivers of environmental heterogeneity and biological diversity across mountain landscapes. Worldwide, glaciers are receding rapidly due to climate change, with important consequences for biodiversity in mountain ecosystems. However, the effects of glacier loss on biodiversity have never been quantified across a mountainous region, primarily due to a lack of adequate data at large spatial and temporal scales. Here, we combine high-resolution biological and glacier change (ca. 1850–2015) datasets for Glacier National Park, USA, to test the prediction that glacier retreat reduces biodiversity in mountain ecosystems through the loss of uniquely adapted meltwater stream species. We identified a specialized cold-water invertebrate community restricted to the highest elevation streams primarily below glaciers, but also snowfields and groundwater springs. We show that this community and endemic species have unexpectedly persisted in cold, high-elevation sites, even in catchments that have not been glaciated in ∼170 y. Future projections suggest substantial declines in suitable habitat, but not necessarily loss of this community with the complete disappearance of glaciers. Our findings demonstrate that high-elevation streams fed by snow and other cold-water sources continue to serve as critical climate refugia for mountain biodiversity even after glaciers disappear.

Montana↗

USGS2020

Modelling pinniped abundance and distribution by combining counts at terrestrial sites and in-water sightings

Pinnipeds are commonly monitored using aerial photographic surveys at land- or ice-based sites, where animals come ashore for resting, pupping, molting, and to avoid predators. Although these counts form the basis for monitoring population change over time, they do not provide information regarding where animals occur in the water, which is often of management and conservation interest. In this study, we developed a hierarchical model that links counts of pinnipeds at terrestrial sites to sightings-at-sea and estimates abundance, spatial distribution, and the proportion of time spent on land (attendance probability). The structure of the model also allows for the inclusion of predictors that may explain variation in ecological and observation processes. We applied the model to Steller sea lions ( Eumetopias jubatus ) in Glacier Bay, Alaska using counts of sea lions from aerial photographic surveys and opportunistic in-water sightings from vessel surveys. Glacier Bay provided an ideal test and application of the model because data are available on attendance probability based on long-term monitoring. We found that occurrence in the water was positively related to proximity to terrestrial sites, as would be expected for a species that engages in central-place foraging. The proportion of sea lions in attendance at terrestrial sites and overall abundance estimates were consistent with reports from the literature and monitoring programs. The model we describe has benefit and utility for park managers who wish to better understand the overlap between pinnipeds and visitors, and the framework that we present has potential for application across a variety of study systems and taxa.

Alaska↗

USGS2020

Parsing complex terrain controls on mountain glacier response to climate forcing

Glaciers are a key indicator of changing climate in the high mountain landscape. Glacier variations across a mountain range are ultimately driven by regional climate forcing. However, changes also reflect local, topographically driven processes such as snow avalanching, snow wind-drifting, and radiation shading as well as the initial glacier conditions such as hypsometry and ice thickness. Here we assess the role of these various terrain influences on change to Little Ice Age (LIA) glaciers in Glacier National Park, U.S.A . With available data for LIA and modern glacier areas, we estimate glacier volumes using simple ice flow assumptions, and topographically driven processes using terrain proxies. At the LIA glacial maxima there were 82 glaciers larger than 0.1 km 2 ranging from 0.11 to 4.97 km 2 . Over the course of the 20 th century, every single LIA glacier decreased in area and 60% (49 glaciers) diminished to below the 0.1 km 2 threshold. Glaciers with large initial area (>1.5 km 2 ) at the end of LIA persisted. Within the intermediate size class (0.5 km 2 < area < 1.5 km 2 ), LIA glacier persistence is poorly explained by initial glacier volume, ice thickness, or elevation. Instead, wind exposure is an important explanatory factor. Our analysis demonstrates the complex response of cirque glaciers to post-LIA climate change in this region: individual glaciers have not necessarily undergone equivalent and synchronous change. Nevertheless, that all glaciers in this mountain range experienced retreat demonstrates that local processes mediated adjustments of some glaciers, but completely decoupled none from the regional climate forcing.

Montana↗

USGS2019

Glacier retreat in Glacier National Park, Montana

Currently, the volume of land ice on Earth is decreasing, driving consequential changes to global sea level and local stream habitat. Glacier retreat in Glacier National Park, Montana, U.S.A., is one example of land ice loss and glacier change. The U.S. Geological Survey Benchmark Glacier Project conducts glaciological research and collects field measurements across select North American glaciers, including Sperry Glacier located in Glacier National Park. These records provide direct evidence of glacier change. The amount and timing of Earth’s land ice that will be lost in the future depends on the future trajectory of greenhouse gas emissions. Ongoing glaciological research will enhance physical understanding and forecasting capabilities; however, the continued retreat and ice mass loss of small glaciers in Glacier National Park is certain.

Montana↗

USGS2019

Glacier recession since the Little Ice Age: Implications for water storage in a Rocky Mountain landscape

Glacial ice is a significant influence on local climate, hydrology, vegetation, and wildlife. We mapped a complete set of glacier areas from the Little Ice Age (LIA) using very high-resolution satellite imagery (30-cm) within Glacier National Park, a region that encompasses over 400,000 hectares. We measured glacier change across the park using LIA glacier area as a baseline and used this to estimate change in glacier area and volume over time. An estimated 146 glaciers existed within the current boundaries of Glacier National Park during the LIA. By 2005, only 51 (35%) persisted. Nearly 90% of LIA glaciers had lost 50% of their area by 2005. This decrease in glacier area equates to an estimated loss of ice volume of 1.52 km3, or 1.37 km3 of water storage, roughly equivalent to 40% of Lake McDonald, the largest lake in the park. Understanding rates of deglaciation and implications for water storage and use can assist local resource managers and downstream communities in planning for change.

Montana↗

USGS2019

Microbial assemblages reflect environmental heterogeneity in alpine streams

Alpine streams are dynamic habitats harboring substantial biodiversity across small spatial extents. The diversity of alpine stream biota is largely reflective of environmental heterogeneity stemming from varying hydrological sources. Globally, alpine stream diversity is under threat as meltwater sources recede and stream conditions become increasingly homogeneous. Much attention has been devoted to macroinvertebrate diversity in alpine headwaters, yet to fully understand the breadth of climate change threats, a more thorough accounting of microbial diversity is needed. We characterized microbial diversity (specifically Bacteria and Archaea) of 13 streams in two disjunct Rocky Mountain subranges through 16S rRNA gene sequencing. Our study encompassed the spectrum of alpine stream sources (glaciers, snowfields, subterranean ice, and groundwater) and three microhabitats (ice, biofilms, and streamwater). We observed no difference in regional (γ) diversity between subranges but substantial differences in diversity among (β) stream types and microhabitats. Within‐stream (α) diversity was highest in groundwater‐fed springs, lowest in glacier‐fed streams, and positively correlated with water temperature for both streamwater and biofilm assemblages. We identified an underappreciated alpine stream type—the icy seep—that are fed by subterranean ice, exhibit cold temperatures (summer mean <2°C), moderate bed stability, and relatively high conductivity. Icy seeps will likely be important for combatting biodiversity losses as they contain similar microbial assemblages to streams fed by surface ice yet may be buffered against climate change by insulating debris cover. Our results show that the patterns of microbial diversity support an ominous trend for alpine stream biodiversity; as meltwater sources decline, stream communities will become more diverse locally, but regional diversity will be lost. Icy seeps, however, represent a source of optimism for the future of biodiversity in these imperiled ecosystems.

Montana↗

USGS2019

Mapping a keystone shrub species, huckleberry (Vaccinium membranaceum), using seasonal colour change in the Rocky Mountains

Black huckleberries ( Vaccinium membranaceum ) provide a critical food resource to many wildlife species, including apex omnivores such as the grizzly bear ( Ursus arctos ), and play an important socioeconomic role for many communities in western North America, especially indigenous peoples. Remote sensing imagery offers the potential for accurate landscape-level mapping of huckleberries because the shrub changes colour seasonally. We developed two methods, for local and regional scales, to map a shrub species using leaf seasonal colour change from remote sensing imagery. We assessed accuracy with ground-based vegetation surveys. The high-resolution supervised random forest classification from one-meter resolution National Agricultural Imagery Program (NAIP) imagery achieved an overall accuracy of 75.31% (kappa = 0.26). The approach using multi-temporal 30-meter Landsat imagery similarly had an overall accuracy of 79.19% (kappa = .31). We found underprediction error was related to higher forest cover and a lack of visible colour change on the ground in some plots. Where forest cover was low, both models performed better. In areas with <10% forest cover, the high-resolution classification achieved an accuracy of 80.73% (kappa = 0.48), while the Landsat model had an accuracy of 82.55% (kappa = 0.47). Based on the fine-scale predictions, we found that 94% of huckleberry shrubs identified in our study area of Glacier National Park, Montana, USA are over 100 meters from human recreation trails. This information could be combined with productivity and phenology information to estimate the timing and availability of food resources for wildlife and to provide managers with a tool to identify and manage huckleberries. The development of the multi-temporal Landsat models sets the stage for assessment of impacts of disturbance at regional scales on this ecologically, culturally, and economically important shrub species. Our approach to map huckleberries is straightforward, efficient and accessible to wildlife and environmental managers and researchers in diverse fields.

Montana↗

USGS2018

Using spatially‐explicit capture–recapture models to explain variation in seasonal density patterns of sympatric ursids

Understanding how environmental factors interact to determine the abundance and distribution of animals is a primary goal of ecology, and fundamental to the conservation of wildlife populations. Studies of these relationships, however, often assume static environmental conditions, and rarely consider effects of competition with ecologically similar species. In many parts of their shared ranges, grizzly bears Ursus arctos and American black bears U. americanus have nearly complete dietary overlap and share similar life history traits. We therefore tested the hypothesis that density patterns of both bear species would reflect seasonal variation in available resources, with areas of higher primary productivity supporting higher densities of both species. We also hypothesized that interspecific competition would influence seasonal density patterns. Specifically, we predicted that grizzly bear density would be locally reduced due to the ability of black bears to more efficiently exploit patchy food resources such as seasonally abundant fruits. To test our hypotheses, we used detections of 309 grizzly and 597 black bears from two independent genetic sampling methods in spatially‐explicit capture–recapture (SECR) models. Our results suggest grizzly bear density was lower in areas of high black bear density during spring and summer, although intraspecific densities were also important, particularly during the breeding season. Black bears had lower densities in areas of high grizzly bear density in spring; however, density of black bears in early and late summer was best explained by primary productivity. Our results are consistent with the hypothesis that smaller‐bodied, more abundant black bears may influence the density patterns of behaviorally‐dominant grizzly bears through exploitative competition. We also suggest that seasonal variation in resource availability be considered in efforts to relate environmental conditions to animal density.

Montana↗