USGS Science⌕ Search

Geology topics

Karl W. Birkeland

Publications and source records attributed to Karl W. Birkeland.

15 recordsLinked to original sources

Stability changes after weak layer collapse: Insights from an unexpected dataset

While conducting avalanche mitigation or traveling through the backcountry, avalanche professionals sometimes experience a "whumpf": a collapse-mode crack accompanied by surface cracking that does not produce an avalanche. These whumpfs are a clear warning sign of instability in the surrounding terrain. But assessing the collapsed slope itself, if it is steep enough to avalanche, carries much more uncertainty. Can slopes that whumpf and are steep enough to slide still avalanche? The answer to this question supports avalanche mitigators and ski guides in assessing whether to open such slopes. Previous research on collapsed slopes is limited to one study conducted before tests that estimate crack propagation propensity, such as the Propagation Saw Test (PST), were developed. This work showed that weak layer shear strength decreased immediately after collapsing and then increased over time. However, no research has examined crack propagation propensity after a whumpf. Studying propagation propensity is key for assessing how whumpfs immediately affect slope stability. To improve our understanding of changes in propagation propensity following weak layer collapse, we conducted PSTs on a test slope in the hours and days after an unexpected collapse occurred on a layer of surface hoar. Two PSTs at our study site three days prior to this collapse provided an estimate of weak layer propagation propensity. PST cut lengths decreased immediately following the collapse, suggesting that the collapse increased the propagation propensity. However, over the next 4-5 hours, PST cut lengths steadily increased at our site, likely due to ongoing sintering and strengthening in the collapsed weak layer. By the end of day two, we were no longer able to initiate propagation, suggesting substantial stabilization of the weak layer. Our results support prior findings and observations suggesting that slopes may remain capable of producing avalanches immediately after a collapse but likely stabilize over the next day or two. Although this research does not provide answers about precisely when a slope becomes stable after a collapse, it does highlight uncertainties in our understanding and interpretation of the stability of slopes that have whumped.

Montana↗

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↗

Changing drivers of regional large magnitude avalanche frequency throughout Colorado, USA

Large magnitude snow avalanches (destructive size ≥ D3) impact settlements, transportation corridors, and public safety worldwide. In Colorado, United States, avalanches have killed more people than any other natural hazard since 1950. In March 2019, a large magnitude avalanche cycle occurred throughout the entire mountainous portion of Colorado resulting in more than 1000 reported avalanches during a two-week period. Nearly 200 of these avalanches were size D4 or larger with at least three D5 avalanches. However, placing this 2019 large magnitude avalanche cycle in historic context requires data prior to the instrumental record. Here, we paired tree disturbance data from dendrochronology (1698 to 2020) with meteorological data from the modeled and instrumental record (1901 to 2020) to understand the frequency and climate drivers of large magnitude snow avalanche cycles. The extensive number of downed trees from the 2019 avalanche cycle allowed us to collect 1,188 cross-sections and cores from 1023 individual trees within 24 avalanche paths across the state. From these samples we identified 4135 avalanche-related growth disturbances. We employed a strategic nested sampling design to spatially aggregate avalanche frequency from individual avalanche paths, to counties, to three major sub-regions (i.e., north, central, and south), and across the entire region (i.e., state of Colorado). Over a period spanning more than three centuries (1698 to 2020), we identified 76 avalanche years within 24 individual avalanche paths. Large magnitude avalanche event frequency varied across paths and sub-regions with several notable region-wide avalanche cycles. Both tree-ring and historical written records highlighted 1899 as a year with widespread and large magnitude avalanche activity similar to the March 2019 avalanche cycle. Since the early-20th century (1900 to 2020) regional avalanche probability declined significantly in parallel with decreasing snowpack throughout Colorado. Similarly, dominant avalanche regimes shifted from large magnitude regional cycles driven by above average snowfall years over most of the record, to regional avalanche cycles occurring more commonly in average to low snow years since 1988. In recent decades, a lack of December precipitation and above average March precipitation characterized years with regional large magnitude avalanche activity. Even with declining snow water equivalent, truly extreme regional large magnitude avalanche cycles remain possible – as demonstrated by the 2019 cycle. This underscores that rare but high-impact events are not eliminated by long-term trends. Understanding the changing snow and weather drivers and subsequent behavior of large magnitude avalanche cycles across multiple spatial scales may improve avalanche forecasting and the products and mitigations strategies developed by structural engineers to mitigate avalanche danger. This can decrease the avalanche risk to the public and improve infrastructure design in avalanche terrain.

Colorado↗

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↗

The relationship between whumpf observations and avalanche activity in Colorado, USA

Triggering whumpfs is a primary indicator of unstable snowpack conditions. Although backcountry travelers and avalanche forecasters rely on whumpfs as a warning sign of potential avalanches, there is little formal research to confirm this relationship. This study investigated the temporal correlation between whumpfs and avalanche activity in data from Colorado's Front Range and southern San Juan Mountains between the winters of 2010/11 and 2022/23. To assess changing conditions over a variety of seasons, we compared the timing of whumpfs and avalanches to the total snow depth at a representative site. We used a 13-inch (33 cm) rolling-window average snow depth versus the median for observed whumpfs, and small avalanches (D1 to D1.5), and large to very large avalanches (D2 and greater). Our results support informal observations that whumpfs are important indicators of avalanche activity, especially at shallower snow depths. Later in the season, when snow depths are deeper and basal weak layers become more difficult to trigger, whumpfs become less common even during periods of increasing avalanche activity. Some of our results may be due to the thin, weak, and wind-affected snow in the Colorado Front Range, where whumpfing typically occurs due to collapsing basal depth hoar. Our findings are important for backcountry travelers assessing stability and for backcountry avalanche forecasters communicating conditions to the public. Our data show that although whumpfs generally indicate unstable conditions and correlate with avalanche activity, the largest avalanches of the winter may not always be preceded by whumpfing.

Colorado↗

Using tree rings to compare Colorado’s 2019 avalanche cycle to previous large avalanche cycles

Large magnitude avalanches (size ≥D3) impact settlements, transportation corridors, and public safety worldwide. In Colorado, United States, avalanches have killed more people than any other natural hazard since 1950. In March 2019, a historically large magnitude avalanche cycle occurred throughout the entire mountainous portion of Colorado resulting in more than 1000 reported avalanches during a 2-week period. Nearly 200 of these avalanches were size D4 or larger with at least three D5 avalanches. The extensive number of downed trees from this avalanche cycle allowed us to collect 1188 cross-sections and cores from 1023 unique trees within 24 avalanche paths across the state. We recorded 4135 growth disturbances in these samples. These data comprise the largest known avalanche tree-ring dataset in the world. We employed a strategic nested sampling design to account for scale by including several individual avalanche paths within a given drainage to create sub-regions and then sampled six major sub-regions (counties) throughout the greater region (state). We identified 76 avalanche years within 24 individual avalanche paths from 1698 to 2020. Large magnitude empirical avalanche event frequency varied across paths and sub-regions. Our results indicate the most widespread avalanche cycle in our study area prior to 2019 occurred in 1899, where 12 avalanche paths show evidence of large magnitude avalanche activity. Historical records also highlight 1899 as a year with widespread and large magnitude avalanche activity. These results indicate the avalanche cycle of March 2019 was of similar magnitude. Understanding the spatial extent and return frequency of large magnitude avalanche cycles across multiple spatial scales, from individual paths to an entire state, helps avalanche forecasters improve their products and mitigation strategies and assists infrastructure planners when designing and planning in avalanche terrain.

Colorado↗

Big avalanches in a changing climate: Using tree-ring derived avalanche chronologies to examine avalanche frequency across multiple climate types

Large-magnitude snow avalanches pose a hazard to humans and infrastructure worldwide. Analyzing the spatiotemporal behavior of avalanches and the contributory climate factors is important for understanding historical variability in climate-avalanche relationships as well as improving avalanche forecasting. This study uses established dendrochronological methods to develop long-term regional avalanche chronologies for three different climate types: high-latitude maritime climate of southeast Alaska, intermountain climate of the northern Rocky Mountains, and continental climate of Colorado. In the maritime study area, we collected 434 cross sections throughout six avalanche paths near Juneau, Alaska. This resulted in 2706 identified avalanche growth disturbances between year 1720 and 2018 Common Era (CE), which allowed us to reconstruct 82 years with large magnitude avalanche activity across three sub-regions. By combining this tree-ring derived avalanche dataset with a suite of climate and atmospheric variables and applying a generalized linear model to fit a binomial regression, we found February and March precipitation and the Oceanic Niño Index (ONI) were significant predictors of large magnitude avalanche activity in the southeast Alaska study area. In the intermountain climate study area, tree-rings from 647 trees exhibited 2134 avalanche-related growth disturbances in the northern Rocky Mountains of northwest Montana from 1867 to 2019. The data show that the amount of snowpack across the northern Rocky Mountain region is directly related to avalanche probability. Coincident with warming and regional snowpack reductions, a decline of ~ 14% (~ 2% per decade) in overall large magnitude avalanche probability is apparent through the period 1950–2017 CE. In the continental climate of Colorado, we sampled 24 avalanche paths throughout the state and collected 1188 total samples with 4135 identified growth disturbances from 1698 to 2019. Preliminary results suggest years with large magnitude avalanche activity across the sub-regions of this study area are generally characterized by stormy winters with above average snowpack development but that early and late winter temperature and precipitation also play an important role in large avalanche activity. Characterizing historical climate-avalanche relationships across different climate types provides a broad baseline for understanding potential future changes in avalanche activity. Overall, this work helps forecasters and planners better understand the influence of climate on large magnitude avalanche frequency, and how potential changes in avalanche character and occurrence will affect their operations in the context of a warming climate.

Alaska, Colorado, Montana↗

Assessing the seasonal evolution of snow depth spatial variability and scaling in complex mountain terrain

Dynamic natural processes govern snow distribution in mountainous environments throughout the world. Interactions between these different processes create spatially variable patterns of snow depth across a landscape. Variations in accumulation and redistribution occur at a variety of spatial scales, which are well established for moderate mountain terrain. However, spatial patterns of snow depth variability in steep, complex mountain terrain have not been fully explored due to insufficient spatial resolutions of snow depth measurement. Recent advances in uncrewed aerial systems (UASs) and structure from motion (SfM) photogrammetry provide an opportunity to map spatially continuous snow depths at high resolutions in these environments. Using UASs and SfM photogrammetry, we produced 11 snow depth maps at a steep couloir site in the Bridger Range of Montana, USA, during the 2019–2020 winter. We quantified the spatial scales of snow depth variability in this complex mountain terrain at a variety of resolutions over 2 orders of magnitude (0.02 to 20 m) and time steps (4 to 58 d) using variogram analysis in a high-performance computing environment. We found that spatial resolutions greater than 0.5 m do not capture the complete patterns of snow depth spatial variability within complex mountain terrain and that snow depths are autocorrelated within horizontal distances of 15 m at our study site. The results of this research have the potential to reduce uncertainty currently associated with snowpack and snow water resource analysis by documenting and quantifying snow depth variability and snowpack evolution on relatively inaccessible slopes in complex terrain at high spatial and temporal resolutions.

Montana↗

Climate drivers of large magnitude snow avalanche years in the U.S. northern Rocky Mountains

Large magnitude snow avalanches pose a hazard to humans and infrastructure worldwide. Analyzing the spatiotemporal behavior of avalanches and the contributory climate factors is important for understanding historical variability in climate-avalanche relationships as well as improving avalanche forecasting. We used established dendrochronological methods to develop a long-term (1867–2019) regional avalanche chronology for the Rocky Mountains of northwest Montana using tree-rings from 647 trees exhibiting 2134 avalanche-related growth disturbances. We then used principal component analysis and a generalized linear autoregressive moving average model to examine avalanche-climate relationships. Historically, large magnitude regional avalanche years were characterized by stormy winters with positive snowpack anomalies, with avalanche years over recent decades increasingly influenced by warmer temperatures and a shallow snowpack. The amount of snowpack across the region, represented by the first principal component, is shown to be directly related to avalanche probability. Coincident with warming and regional snowpack reductions, a decline of ~ 14% (~ 2% per decade) in overall large magnitude avalanche probability is apparent through the period 1950–2017. As continued climate warming drives further regional snowpack reductions in the study region our results suggest a decreased probability of regional large magnitude avalanche frequency associated with winters characterized by large snowpacks and a potential increase in large magnitude events driven by warming temperatures and spring precipitation.

Montana↗

A regional spatio-temporal analysis of large magnitude snow avalanches using tree rings

Snow avalanches affect transportation corridors and settlements worldwide. In many mountainous regions, robust records of avalanche frequency and magnitude are sparse or non-existent. However, dendrochronological methods can be used to fill this gap and infer historical avalanche patterns. In this study, we developed a tree-ring-based avalanche chronology for large magnitude avalanche events (size &#x2265; &#x223C; D 3 "> ≥ ∼ D 3 ) using dendrochronological techniques for a portion of the US northern Rocky Mountains. We used a strategic sampling design to examine avalanche activity through time and across nested spatial scales (i.e., from individual paths, four distinct subregions, and the region). We analyzed 673 samples in total from 647 suitable trees collected from 12 avalanche paths from which 2134 growth disturbances were identified over the years 1636 to 2017 CE. Using existing indexing approaches, we developed a regional avalanche activity index to discriminate avalanche events from noise in the tree-ring record. Large magnitude avalanches, common across the region, occurred in 30 individual years and exhibited a median return interval of approximately 3 years (mean = 5.21 years). The median large magnitude avalanche return interval (3–8 years) and the total number of avalanche years (12–18) varies throughout the four subregions, suggesting the important influence of local terrain and weather factors. We tested subsampling routines for regional representation, finding that sampling 8 random paths out of a total of 12 avalanche paths in the region captures up to 83 % of the regional chronology, whereas four paths capture only 43 % to 73 %. The greatest value probability of detection for any given path in our dataset is 40 %, suggesting that sampling a single path would capture no more than 40 % of the regional avalanche activity. Results emphasize the importance of sample size, scale, and spatial extent when attempting to derive a regional large magnitude avalanche event chronology from tree-ring records.

Montana↗

Research Note: How old are the people who die in avalanches? A look into the ages of avalanche victims in the United States (1950-2018)

Since the winter of 1950-1951, 1084 individuals perished in snow avalanches in the United States. In this study, we analyze the ages of those killed (n=900) by applying non-parametric methods to annual median ages and for age groups and primary activity groups. Change point detection results suggest a significant change in 1990 in the median age of avalanche fatalities. Significant positive trends exist for both the 1950 to 2018 and 1990 to 2018 median age of victims. The median age of victims from 1950 to 1989 is 27 and 33 from 1990 to 2018. Since 1990, the 30-39 and 40-49 age groups are the only age category to exhibit a positive trend in the number of fatalities. There is no significant difference in median ages between snowmobilers and other categories. These results can be used to enhance avalanche education and forecasting efforts in the United States.

Journal of Outdoor Recreation and Tourism↗

Detecting snow depth change in avalanche path starting zones using uninhabited aerial systems and structure from motion photogrammetry

Understanding snow depth distribution and change is useful for avalanche forecasting and mitigation, runoff forecasting, and infrastructure planning. Advances in remote sensing are improving the ability to collect snow depth measurements. The development of structure from motion (SfM), a photogrammetry technique, combined with the use of uninhabited aerial systems (UASs) allows for high resolution mapping of snow depth over complex terrain. The primary objective of this study was to determine the feasibility and efficacy of SfM to examine snow depth distribution and variability in complex terrain such as avalanche path starting zones at multiple times during the season. We used a 3DR Solo quadcopter UAS equipped with a Ricoh GR II camera at 90 m above ground level to acquire images of one avalanche starting zone in northwest Montana, USA. We also placed 4 to 13 ground control points (GCPs) around the area of interest to avoid traveling in steep, avalanche terrain. Ground control measurements resulted in 5 to10 cm horizontal accuracy and 5 to 15 cm vertical accuracy for 90 to 95 % of the collected points (a minimum of 100 points collected at each GCP). In-situ measurements of snow depth difference between sampling days ranged from 20 to 60 cm. We processed the images to create point clouds and digital surface models (DSMs). The resolution of the resultant DSMs was approximately 5 cm. Preliminary DSM and point cloud differencing efforts suggest relative change detection of snow depth at 5 to 15 cm resolution. The use of these relatively low cost and easily accessible methods of snow depth data collection will enhance accuracy of snow depth change estimates in starting zones and can be used to inform avalanche forecasting and mitigation efforts.

Conference Paper↗

Identifying major avalanche years from a regional tree-ring based avalanche chronology for the U.S. Northern Rocky Mountains

Avalanches not only pose a major hazard to people and infrastructure, but also act as an important ecological disturbance. In many mountainous regions in North America, including areas with existing transportation corridors, reliable and consistent avalanche records are sparse or non-existent. Thus, inferring long-term avalanche patterns and associated contributory climate and weather factors requires the use of dendrochronological methods...

Conference Paper↗

Meteorological variables to aid forecasting deep slab avalanches on persistent weak layers

Deep slab avalanches are particularly challenging to forecast. These avalanches are difficult to trigger, yet when they release they tend to propagate far and can result in large and destructive avalanches. We utilized a 44-year record of avalanche control and meteorological data from Bridger Bowl ski area in southwest Montana to test the usefulness of meteorological variables for predicting seasons and days with deep slab avalanches. We defined deep slab avalanches as those that failed on persistent weak layers deeper than 0.9 m, and that occurred after February 1st. Previous studies often used meteorological variables from days prior to avalanches, but we also considered meteorological variables over the early months of the season. We used classification trees and random forests for our analyses. Our results showed seasons with either dry or wet deep slabs on persistent weak layers typically had less precipitation from November through January than seasons without deep slabs on persistent weak layers. Days with deep slab avalanches on persistent weak layers often had warmer minimum 24-hour air temperatures, and more precipitation over the prior seven days, than days without deep slabs on persistent weak layers. Days with deep wet slab avalanches on persistent weak layers were typically preceded by three days of above freezing air temperatures. Seasonal and daily meteorological variables were found useful to aid forecasting dry and wet deep slab avalanches on persistent weak layers, and should be used in combination with continuous observation of the snowpack and avalanche activity.

Montana↗

Meteorological variables associated with deep slab avalanches on persistent weak layers

Deep slab avalanches are a particularly challenging avalanche forecasting problem. These avalanches are typically difficult to trigger, yet when they are triggered they tend to propagate far and result in large and destructive avalanches. For this work we define deep slab avalanches as those that fail on persistent weak layers deeper than 0.9m (3 feet), and that occur after February 1st. We utilized a 44-year record of avalanche control and meteorological data from Bridger Bowl Ski Area to test the usefulness of meteorological variables for predicting deep slab avalanches. As in previous studies, we used data from the days preceding deep slab cycles, but we also considered meteorological metrics over the early months of the season. We utilized classification trees for our analyses. Our results showed warmer temperatures in the prior twenty-four hours and more loading over the seven days before days with deep slab avalanches on persistent weak layers. In line with previous research, extended periods of above freezing temperatures led to days with deep wet slab avalanches on persistent weak layers. Seasons with either dry or wet avalanches on deep persistent weak layers typically had drier early months, and often had some significant snow depth prior to those dry months. This paper provides insights for ski patrollers, guides, and avalanche forecasters who struggle to forecast deep slab avalanches on persistent weak layers late in the season.

Conference Paper↗