Ice-marginal lakes impact glacier mass balance, water resources, and ecosystem dynamics and can produce catastrophic glacial lake outburst floods (GLOFs) via sudden drainage. Multitemporal inventories of ice-marginal lakes are a critical first step in understanding the drivers of historic change, predicting future lake evolution, and assessing GLOF hazards. Here, we use Landsat-era satellite imagery and supervised classification to semi-automatically delineate lake outlines for four ∼5-year time periods between 1984 and 2019 in Alaska and northwest Canada. Overall, ice-marginal lakes in the region have grown in total number (+183 lakes, 38 % increase) and area (+483 km², 59 % increase) between the time periods of 1984–1988 and 2016–2019. However, changes in lake numbers and area were notably unsteady and nonuniform. We demonstrate that lake area changes are connected to dam type (moraine, bedrock, ice, or supraglacial) and topological position (proglacial, detached, unconnected, ice, or supraglacial), with important differences in lake behavior between the sub-groups. In strong contrast to all other dam types, ice-dammed lakes decreased in number (six fewer, 9 % decrease) and area (−51 km², 40 % decrease), while moraine-dammed lakes increased (56 more, 26 % and +479 km², 87 % increase for number and area, respectively) at a faster rate than the average when considering all dam types together. Proglacial lakes experienced the largest area changes and rate of change out of any lake position throughout the period of study and moraine-dammed lakes which experienced the largest increases are associated with clean-ice glaciers (<19 % debris cover). By tracking individual lakes through time and categorizing lakes by dam type, subregion, and topological position, we are able to parse trends that would otherwise be aliased if these characteristics were not considered. This work highlights the importance of such lake characterization when performing ice-marginal lake inventories and provides insight into the physical processes driving recent ice-marginal lake evolution.

Glacial lakes can form and grow due to glacial retreat, and rapid lake drainage can produce destructive floods. Outburst flood compilations show a temporal increase in frequency; however, recent studies highlight the role of observational bias, creating uncertainty about current and future glacial-lake hazards. Here, we focus on the Alaska region, which generated a third of previously documented outbursts globally. Using multitemporal satellite imagery, we documented 1150 drainages from 106 ice-dammed lakes between 1985 and 2020. Documented events became more frequent over time, however, accounting for increasing image availability reveals no significant increase occurred. Most lakes decreased in area and volume, suggesting a reduction in regional flood hazard. Our satellite-based approach documented 60% more events in a 35-year period than had previously been documented over 100 years. This suggests that outburst floods have historically been underreported and warrants systematic study of other regions.

The Chaos Canyon landslide, which collapsed on the afternoon of 28 June 2022 in Rocky Mountain National Park, presents an opportunity to evaluate instabilities within alpine regions faced with a warming and dynamic climate. Video documentation of the landslide was captured by several eyewitnesses and motivated a rapid field campaign. Initial estimates put the failure area at 66 630 m², with an average elevation of 3555 m above sea level. We undertook an investigation of previous movement of this landslide, measured the volume of material involved, evaluated the potential presence of interstitial ice and snow within the failed deposit, and examined potential climatological impacts on the collapse of the slope. Satellite radar and optical measurements were used to calculate deformation of the landslide in the 5 years leading up to collapse. From 2017 to 2019, the landslide moved ∼5 m yr⁻¹, accelerating to 17 m yr⁻¹ in 2019. Movement took place through both internal deformation and basal sliding. Climate analysis reveals that the collapse took place during peak snowmelt, and 2022 followed 10 years of higher than average positive degree day sums. We also made use of slope stability modeling to test what factors controlled the stability of the area. Models indicate that even a small increase in the water table reduces the factor of safety to <1, leading to failure. We posit that a combination of permafrost thaw from increasing average temperatures, progressive weakening of the basal shear zone from several years of movement, and an increase in pore-fluid pressure from snowmelt led to the 28 June collapse. Material volumes were estimated using structure from motion (SfM) models incorporating photographs from two field expeditions on 8 July 2022 – 10 d after the slide. Detailed mapping and SfM models indicate that ∼1 258 000 ± 150 000 m³ of material was deposited at the slide toe and ∼1 340 000 ± 133 000 m³ of material was evacuated from the source area. The Chaos Canyon landslide may be representative of future dynamic alpine topography, wherein slope failures become more common in a warming climate.

Diminishing glaciers affect streamflow, and given the extent of glaciers in Alaska and adjacent Canada, continued glacier mass loss is likely to have profound effects on ecosystems sensitive to runoff. The effects of glacier mass loss on streamflow are likely to vary across the wide ranges of basin size, glacier cover, and precipitation in this region. In this study, we use U.S. Geological Survey (USGS) streamflow data with satellite‐based glacier volume change estimates to quantify how glacier mass loss subsidized streamflow over the 2000–2019 period for 116 glacierized basins. We examine interannual variability in that subsidy at three USGS‐monitored glaciers to explore the ability of the subsidy to buffer streamflow derived solely from precipitation. We found the relative importance of percent glacier cover on streamflow magnitude increases in drier basins. In the driest basins, glaciers produced 40 times greater percent glacier mass loss subsidies to streamflow for the percent glacier cover compared to the wettest basins. While the subsidy from glacier mass loss buffers interannual variability in streamflow to varying degrees, it can also increase streamflow variability. Smaller amounts of percent glacier cover are needed to produce summer‐melt‐dominated seasonal flow regimes in drier basins than in wetter basins. Decreasing glacier cover will eventually decrease summer streamflow, increasing spring streamflow in drier basins, and attenuating seasonality with increasing spring and autumnal streamflow in wetter basins. Quantifying the downstream effects of continued glacier mass loss without the computational expense of a hydrological model is broadly applicable in this changing climate.

Proglacial lakes at glacier termini have received widespread attention in the literature for their role in accelerating melt, velocity and contributing to cryospheric hazards. Although global and regional inventories for both glaciers and lakes exist, lake-terminating glaciers have not been consistently identified at the global scale. Based on the most recent global glacier inventory (RGI 7.0), which so far identifies marine-terminating glaciers in most regions, but not lake-terminating glaciers, we present a global inventory of lake-terminating glaciers, differentiating between three categories based on the degree of contact between the glacier terminus and any proglacial lakes. Contributors manually assigned categories to glaciers using satellite imagery from as close to the target date of 2000 as possible, aided by regional lake inventories where available. The resulting dataset corresponds to the year 2000 (±1.5), matching to the timestamp of RGI 7.0 outlines (2001 ± 6.2). We find that of 274 531 glaciers worldwide, 1.4 % terminate in lakes, with regional percentages varying between 0.5 % and 6.7 % across the 19 RGI regions. These glaciers account for 11.4 % of global glacier area (0.2 % to 41.8 % across regions). With multiple submissions available for 1260 individual glaciers, we find mapping conflicts between contributors to be low (6.7 %). The lake termini data set is available at https://doi.org/10.5281/zenodo.15524733 (Steiner et al., 2025) as well as at https://github.com/GLIMS-RGI/lake_terminating (last access: 2 February 2026). This dataset is integrated into the forthcoming update to the RGI, v7.1.

Lakes in direct contact with glaciers (ice-marginal lakes) are found across alpine and polar landscapes. Many studies characterize ice-marginal lake behavior over multi-decadal timescales using either episodic ~annual images or multi-year mosaics. However, ice-marginal lakes are dynamic features that experience short-term (i.e., day to year) variations in area and volume superimposed on longer-term trends. Through aliasing, this short-term variability could result in erroneous long-term estimates of lake change. We develop and implement an automated workflow in Google Earth Engine to quantify monthly behavior of ice-marginal lakes between 2013 and 2019 across south-central Alaska using Landsat 8 imagery. We employ a supervised Mahalanobis minimum-distance land cover classifier incorporating three datasets found to maximize classifier performance: shortwave infrared imagery, the normalized difference vegetation index (NDVI), and spatially filtered panchromatic reflectance. We observe physically-meaningful ice-marginal lake area variance on sub-annual timescales, with the median area fluctuation of an ice-marginal lake found to be 10.8% of its average area. The median signal (slow lake growth) to noise (physically-meaningful short-term area variability) ratio is 1.5:1, indicating that short-term variability is responsible for ~33% of observed area change in the median ice-marginal lake. The magnitude of short-term area variability is similar for ice-marginal and nonglacial lakes, suggesting that the cause of observed variations is not of glacial origin. These data provide a new context for interpreting behaviors observed in multi-decadal studies and encourage attention to sub-annual behavior of ice-marginal lakes even in long-term studies.