Geology Reports⌕ Search

Geology topics

Louis C. Sass

Publications and source records attributed to Louis C. Sass.

At least 19 recordsLinked to original sources
USGS2025

Automated snow cover detection on mountain glaciers usingspaceborne imagery and machine learning

Tracking the extent of seasonal snow on glaciers over time is critical for assessing glacier vulnerability and the response of glacierized watersheds to climate change. Existing snow cover products do not reliably distinguish seasonal snow from glacier ice and firn, preventing their use for glacier snow cover detection. Despite previous efforts to classify glacier surface facies using machine learning on local scales, currently there is no published comparison of machine learning models for classifying glacier snow cover across different satellite image products. We present an automated snow detection workflow for mountain glaciers using supervised machine-learning-based image classifiers and Landsat 8 and 9, Sentinel-2, and PlanetScope satellite imagery. We develop the image classifiers by testing numerous machine learning algorithms with training and validation data from the U.S. Geological Survey Benchmark Glacier Project glaciers. The workflow produces daily to twice monthly time series of several glacier mass balance and snowmelt indicators (snow-covered area, accumulation area ratio, and seasonal snow line) from 2013 to present. Workflow performance is assessed by comparing automatically classified images and snow lines to manual interpretations at each glacier site. The image classifiers exhibit overall accuracies of 92%–98%, K scores of 84%–96%, and F scores of 93%–98% for all image products. The median difference between automatically and manually delineated median snow line altitudes is 31m (IQR of 73to0m)across all image products. The Sentinel-2 classifier (support vector machine) produces the most accurate glacier mass balance and snowmelt indicators and distinguishes snow from ice and f irn the most reliably. Although they are less accurate, the Landsat- and PlanetScope-derived estimates greatly enhance the temporal coverage of observations. The transient accumulation area ratio produces the least noisy time series, making it the most reliable indicator for characterizing seasonal snow trends. The temporally detailed accumulation area ratio time series reveal that the timing of minimum snow cover conditions varies by up to a month between Arctic (63°N) and midlatitude (48°N) sites, underscoring the potential for bias when estimating glacier minimum snow cover conditions from a single late-summer image. Widespread application of our automated snow detection workflow has the potential to improve regional assessments of glacier mass balance, land ice representations within Earth system models, water resources, and the impacts of climate change on snow cover across broad spatial scales.

The Cryosphere↗
USGS2025

Equilibrium line altitudes, accumulation areas, and the vulnerability of glaciers in Alaska

The accumulation area ratio (AAR) of a glacier reflects its current state of equilibrium, or disequilibrium, with climate and its vulnerability to future climate change. Here, we present an inventory of glacier-specific annual accumulation areas and equilibrium line altitudes (ELAs) for over 3000 glaciers in Alaska and northwest Canada (88% of the regional glacier area) from 2018 to 2022 derived from Sentinel-2 imagery. We find that the 5 year average AAR of the entire study area is 0.41, with an inter-annual range of 0.25–0.49. More than 1000 glaciers, representing 8% of the investigated glacier area, were found to have effectively no accumulation area. Summer temperature and winter precipitation from ERA5-Land explained nearly 50% of the inter-annual ELA variability across the entire study region ( R 2 = 0.47 ). An analysis of future climate scenarios (SSP2-4.5) projects that ELAs will rise by ∼170 m on average by the end of the 21st century. Such changes would result in a loss of 25% of the modern accumulation area, leaving a total of 1900 glaciers (22% of the investigated area) with no accumulation area. These results highlight the current state of glacier disequilibrium with modern climate, as well as glacier vulnerability to projected future warming.

Alaska↗
USGS2025

Streamflow response to glacier mass loss varies with basin precipitation across Alaska

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.

Alaska↗
USGS2024

GNSS reflectometry from low-cost sensors for continuous in situ contemporaneous glacier mass balance and flux divergence

Recent advances in remote sensing have produced global glacier surface elevation change data. Parsing these elevation change signals into contributions from the climate (i.e. climatic mass balance) and glacier dynamics (i.e. flux divergence) is critical to enhance our process-based understanding of glacier change. In this study, we evaluate three approaches for direct, continuous measurements of the climatic mass balance, flux divergence, and elevation change at a site on Gulkana Glacier in Alaska using low-cost GNSS sensors, GNSS interferometric reflectometry (GNSS-IR), banded ablation stakes with time-lapse cameras, and combinations thereof. Cumulative climatic mass balance over the season was 4.85 m and the three approaches were within 0.08 m through early July before the snowpack melted, and within 0.28 m through mid-August. The flux divergence increased from 0.52 ± 0.03 cm d -1 before June 3 to roughly 0.73 cm d -1 after June 27. We demonstrate a single GNSS system fixed atop an ablation stake can measure contemporaneous climatic mass balance, flux divergence, and elevation change based on the antenna’s position and GNSS-IR techniques. The ability of these systems to measure glacier mass balance and flux divergence offers unique opportunities for year-round observations on mountain glaciers in the future.

Alaska↗
USGS2024

Direct measurements of firn-density evolution from 2016 to 2022 at Wolverine Glacier, Alaska

Knowledge of snow and firn-density change is needed to use elevation-change measurements to estimate glacier mass change. Additionally, firn-density evolution on glaciers is closely connected to meltwater percolation, refreezing and runoff, which are key processes for glacier mass balance and hydrology. Since 2016, the U.S. Geological Survey Benchmark Glacier Project has recovered firn cores from a site on Wolverine Glacier in Alaska's Kenai Mountains. We use annual horizons in repeat cores to track firn densification and meltwater retention over seasonal and interannual timescales, and we use density measurements to quantify how the firn air content (FAC) changes through time. The results suggest the firn is densifying due primarily to compaction rather than refreezing. Liquid-water retention in the firn is transient, likely due to gravity-fed drainage and irreducible-water-content decreases that accompany decreasing porosity. We show that the uncertainty (±60 kg m −3 ) in the commonly used volume-to-mass conversion factor of 850 kg m −3 is an underestimation when glacier-wide FAC variability exceeds 12% of the glacier-averaged height change. Our results demonstrate how direct measurements of firn properties on mountain glaciers can be used to better quantify the uncertainty in geodetic volume-to-mass conversions.

Alaska↗
USGS2024

Observing glacier elevation changes from spaceborne optical and radar sensors – an inter-comparison experiment using ASTER and TanDEM-X data

Observations of glacier mass changes are key to understanding the response of glaciers to climate change and related impacts, such as regional runoff, ecosystem changes, and global sea level rise. Spaceborne optical and radar sensors make it possible to quantify glacier elevation changes, and thus multi-annual mass changes, on a regional and global scale. However, estimates from a growing number of studies show a wide range of results with differences often beyond uncertainty bounds. Here, we present the outcome of a community-based inter-comparison experiment using spaceborne optical stereo (ASTER) and synthetic aperture radar interferometry (TanDEM-X) data to estimate elevation changes for defined glaciers and target periods that pose different assessment challenges. Using provided or self-processed digital elevation models (DEMs) for five test sites, 12 research groups provided a total of 97 spaceborne elevation-change datasets using various processing approaches. Validation with airborne data showed that using an ensemble estimate is promising to reduce random errors from different instruments and processing methods but still requires a more comprehensive investigation and correction of systematic errors. We found that scene selection, DEM processing, and co-registration have the biggest impact on the results. Other processing steps, such as treating spatial data voids, differences in survey periods, or radar penetration, can still be important for individual cases. Future research should focus on testing different implementations of individual processing steps (e.g. co-registration) and addressing issues related to temporal corrections, radar penetration, glacier area changes, and density conversion. Finally, there is a clear need for our community to develop best practices, use open, reproducible software, and assess overall uncertainty to enhance inter-comparison and empower physical process insights across glacier elevation-change studies.

The Cryosphere↗
USGS2023

How to handle glacier area change in geodetic mass balance

Innovations in geodesy enable widespread analysis of glacier surface elevation change and geodetic mass balance. However, coincident glacier area data are less widely available, causing inconsistent handling of glacier area change. Here we quantify the bias introduced into meters water equivalent (m w.e.) specific geodetic mass balance results when using a fixed, maximum glacier area, and illustrate the bias for five North American glaciers. Sites span latitudes from the northern U.S. Rocky Mountains (48°N) to the Central Alaska Range (63°N) between 1948 and 2021. Results show that fixed (maximum) area treatment subdues the m w.e. mass change signal, underestimating mass balance by up to 19% in our test cases. This bias scales with relative glacier area change and the mass balance magnitude. Thus, the bias for specific geodetic mass balances will be most pronounced across rapidly deglaciating regions. Our analysis underscores the need for temporally resolved glacier area in geodetic mass balance studies.

Journal of Glaciology↗
USGS2022

Beyond glacier-wide mass balances: Parsing seasonal elevation change into spatially resolved patterns of accumulation and ablation at Wolverine Glacier, Alaska

We present spatially distributed seasonal and annual surface mass balances of Wolverine Glacier, Alaska, from 2016 to 2020. Our approach accounts for the effects of ice emergence and firn compaction on surface elevation changes to resolve the spatial patterns in mass balance at 10 m scale. We present and compare three methods for estimating emergence velocities. Firn compaction was constrained by optimizing a firn model to fit three firn cores. Distributed mass balances showed good agreement with mass-balance stakes (RMSE = 0.67 m w.e., r = 0.99, n = 41) and ground-penetrating radar surveys (RMSE = 0.36 m w.e., r = 0.85, n = 9024). Fundamental differences in the distributions of seasonal balances highlight the importance of disparate physical processes, with anomalously high ablation rates observed in icefalls. Winter balances were found to be positively skewed when controlling for elevation, while summer and annual balances were negatively skewed. We show that only a small percent of the glacier surface represents ideal locations for mass-balance stake placement. Importantly, no suitable areas are found near the terminus or in elevation bands dominated by icefalls. These findings offer explanations for the often-needed geodetic calibrations of glaciological time series.

Alaska↗
USGS2022

Uncertainty of ICESat-2 ATL06- and ATL08-derived snow depths for glacierized and vegetated mountain regions

Seasonal snow melt dominates the hydrologic budget across a large portion of the globe. Snow accumulation and melt vary over a broad range of spatial scales, preventing accurate extrapolation of sparse in situ observations to watershed scales. The lidar onboard the Ice, Cloud, and land Elevation, Satellite (ICESat-2) was designed for precise mapping of ice sheets and sea ice, and here we assess the feasibility of snow depth-mapping using ICESat-2 data in more complex and rugged mountain landscapes. We explore the utility of ATL08 Land and Vegetation Height and ATL06 Land Ice Height differencing from reference elevation datasets in two end member study sites. We analyze ∼3 years of data for Reynolds Creek Experimental Watershed in Idaho's Owyhee Mountains and Wolverine Glacier in southcentral Alaska's Kenai Mountains. Our analysis reveals decimeter-scale uncertainties in derived snow depth and glacier mass balance at the watershed scale. Both accuracy and precision decrease as slope increases: the magnitudes of the median and median of the absolute deviation of elevation errors (MAD) vary from ∼0.2 m for slopes <5° to >1 m for slopes >20°. For glacierized regions, failure to account for intra- and inter-annual evolution of glacier surface elevations can strongly bias ATL06 elevations, resulting in under-estimation of the mass balance gradient with elevation. Based on these results, we conclude that ATL08 and ATL06 observations are best suited for characterization of watershed-scale snow depth and mass balance gradients over relatively shallow slopes with thick snowpacks . In these regions, ICESat-2 elevation residual-derived snow depth and mass balance transects can provide valuable watershed scale constraints on terrain parameter- and model-derived estimates of snow accumulation and melt.

Alaska, Idaho↗
USGS2021

The imminent calving retreat of Taku Glacier

Along the rugged Southeast Alaska coast, 30 kilometers northeast of the state capital Juneau, a tidewater glacier has largely defied global trends by steadily advancing for most of the past century while most glaciers on Earth retreated. This 55-kilometer-long and nearly 1,500-meter-thick tidewater glacier, named Taku Glacier, or T'aaḵú Ḵwáan Sít'i in the language of the Indigenous Tlingit people, has been the focus of continuous scientific study for more than 70 years. Some records even extend back to the mid-18th century. With this long observation record and the glacier’s year-round accessibility and proximity to Juneau and adjacent research facilities, Taku provides an unparalleled locale to study tidewater glaciers and their response to Earth’s rapidly changing climate.

Eos, American Geophysical Union↗
USGS2020

Explaining mass balance and retreat dichotomies at Taku and Lemon Creek Glaciers, Alaska

We reanalyzed mass balance records at Taku and Lemon Creek Glaciers to better understand the relative roles of hypsometry, local climate and dynamics as mass balance drivers. Over the 1946–2018 period, the cumulative mass balances diverged. Tidewater Taku Glacier advanced and gained mass at an average rate of +0.25±0.28 m w.e. a–1, contrasting with retreat and mass loss of –0.60±0.15 m w.e. a-1 at land-terminating Lemon Creek Glacier. The uniform influence of regional climate is demonstrated by strong correlations among annual mass balance and climate data. Regional warming trends forced similar statistically significant decreases in surface mass balance after 1989: –0.83 m w.e. a–1 at Taku Glacier and –0.81 m w.e. a–1 at Lemon Creek Glacier. Divergence in cumulative mass balance arises from differences in glacier hypsometry and local climate. Since 2013 negative mass balance and glacier-wide thinning prevailed at Taku Glacier. These changes initiated terminus retreat, which could increase dramatically if calving begins. The future mass balance trajectory of Taku Glacier hinges on dynamics, likely ending the historic dichotomy between these glaciers.

Alaska↗
USGS2019

Reanalysis of the U.S. Geological Survey Benchmark Glaciers: Long-term insight into climate forcing of glacier mass balance

Mountain glaciers integrate climate processes to provide an unmatched signal of regional climate forcing. However, extracting the climate signal via intercomparison of regional glacier mass balance records can be problematic when methods for extrapolating and calibrating direct glaciological measurements are mixed or inconsistent. To address this problem, we reanalyzed and compared long-term mass balance records from the U.S. Geological Survey Benchmark Glaciers. These five glaciers span maritime and continental climate regimes of the western United States and Alaska. Each glacier exhibits cumulative mass loss since the mid-20th century, with average rates ranging from –0.58 to –0.30 m water equivalent (w.e.) a-1. We produced a set of solutions using different extrapolation and calibration methods to inform uncertainty estimates, which range from 0.22–0.44 m w.e. a-1. Mass losses are primarily driven by increasing summer warming. Continentality exerts a stronger control on mass loss than latitude. Similar to elevation, topographic shading, snow redistribution, and glacier surface features often exert first-order control on mass balance. The reanalysis underscores the value of geodetic calibration to resolve mass balance magnitude, as well as the irreplaceable value of direct measurements in contributing to process-based understanding of glacier mass balance.

Journal of Glaciology↗
USGS2017

Hypsometric control on glacier mass balance sensitivity in Alaska and northwest Canada

Glacier hypsometry provides a first‐order approach for assessing a glacier's response to climate forcings. We couple the Randolph Glacier Inventory to a suite of in situ observations and climate model output to examine potential change for the ∼27,000 glaciers in Alaska and northwest Canada through the end of the 21st century. By 2100, based on Representative Concentration Pathways (RCPs) 4.5–8.5 forcings, summer temperatures are predicted to increase between +2.1 and +4.6°C, while solid precipitation (snow) is predicted to decrease by −6 to −11%, despite a +9 to +21% increase in total precipitation. Snow is predicted to undergo a pronounced decrease in the fall, shifting the start of the accumulation season back by ∼1 month. In response to these forcings, the regional equilibrium line altitude (ELA) may increase by +105 to +225 m by 2100. The mass balance sensitivity to this increase is highly variable, with the most substantive impact for glaciers with either limited elevation ranges (often small (<1 km 2 ) glaciers, which account for 80% of glaciers in the region) or those with top‐heavy geometries, like icefields. For more than 20% of glaciers, future ELAs, given RCP 6.0 forcings, will exceed the maximum elevation of the glacier, resulting in their eventual demise, while for others, accumulation area ratios will decrease by >60%. Our results highlight the first‐order control of hypsometry on individual glacier response to climate change, and the variability that hypsometry introduces to a regional response to a coherent climate perturbation.

Alaska, British Columbia, Yukon↗
USGS2017

Geometry, mass balance and thinning at Eklutna Glacier, Alaska: an altitude-mass-balance feedback with implications for water resources

We analyzed glacier surface elevations (1957, 2010 and 2015) and surface mass-balance measurements (2008–2015) on the 30 km 2 Eklutna Glacier, in the Chugach Mountains of southcentral Alaska. The geodetic mass balances from 1957 to 2010 and 2010 to 2015 are −0.52 ± 0.46 and −0.74 ± 0.10 m w.e. a −1 , respectively. The glaciological mass balance of −0.73 m w.e. a −1 from 2010 to 2015 is indistinguishable from the geodetic value. Even after accounting for loss of firn in the accumulation zone, we found most of the mass loss over both time periods was from a broad, low-slope basin that includes much of the accumulation zone of the main branch. Ice-equivalent surface elevation changes in the basin were −1.0 ± 0.8 m a −1 from 1957 to 2010, and −0.6 ± 0.1 m a −1 from 2010 to 2015, shifting the glacier hypsometry downward and resulting in more negative mass balances: an altitude-mass-balance feedback. Net mass loss from Eklutna Glacier accounts for 7 ± 1% of the average inflow to Eklutna Reservoir, which is entirely used for water and power by Anchorage, Alaska's largest city. If the altitude-mass-balance feedback continues, this ‘deglaciation discharge dividend’ is likely to increase over the short-term before it eventually decreases due to diminishing glacier area.

Alaska↗
USGS2017

Snow and ice

Temperature and precipitation are key determinants of snowpack levels. Therefore, climate change is likely to affect the role of snow and ice in the landscapes and hydrology of the Chugach National Forest region. Downscaled climate projections developed by Scenarios Network for Alaska and Arctic Planning (SNAP) are useful for examining projected changes in snow at relatively fine resolution using a variable called “snowday fraction (SDF),” the percentage of days with precipitation falling as snow. We summarized SNAP monthly SDF from five different global climate models for the Chugach region by 500 m elevation bands, and compared historical (1971–2000) and future (2030–2059) SDF. We found that: Snow-day fraction and snow-water equivalent (SWE) are projected to decline most in late autumn (October to November) and at lower elevations. Snow-day fraction is projected to decrease 23 percent (averaged across five climate models) from October to March, between sea level and 500 m. Between sea level and 1000 m, SDF is projected to decrease by 17 percent between October and March. Snow-water equivalent is projected to decrease most in autumn (October and November) and at lower elevations (below 1500 m), an average of -26 percent for the 2030–2059 period compared to 1971– 2000. Averaged across the cool season and the entire domain, SWE is projected to decrease at elevations below 1000 m because of increased temperature, but increase at higher elevations because of increased precipitation. Compared to 1971–2000, the percentage of the landscape that is snowdominant in 2030–2059 is projected to decrease, and the percentage in which rain and snow are co-dominant (transient hydrology) is projected to increase from 27 to 37 percent. Most of this change is at lower elevations. Glaciers on the Chugach National Forest are currently losing about 6 km3 of ice per year; half of this loss comes from Columbia Glacier (Berthier et al. 2010). Over the past decade, almost all glaciers surveyed within the Chugach have lost mass (with one exception), including glaciers that have advancing termini (Larsen et al. 2015). Glaciers that are not calving into the ocean are typically thinning by 3 m/year at their termini (Larsen et al. 2015). In the future, glaciers not calving into the ocean will retreat and shrink at rates equivalent to or higher than current rates of ice loss (Larsen et al. 2015). Columbia Glacier will likely retreat another 15 km and break into multiple tributaries over the next 20 years before stabilizing. Other tidewater glaciers have uncertain futures, but likely will not advance significantly in coming decades. These impacts will likely affect recreation and tourism through changes in reliable snowpack and access to recreation and viewsheds.

Alaska↗
USGS2015

End-of-winter snow depth variability on glaciers in Alaska

A quantitative understanding of snow thickness and snow water equivalent (SWE) on glaciers is essential to a wide range of scientific and resource management topics. However, robust SWE estimates are observationally challenging, in part because SWE can vary abruptly over short distances in complex terrain due to interactions between topography and meteorological processes. In spring 2013, we measured snow accumulation on several glaciers around the Gulf of Alaska using both ground- and helicopter-based ground-penetrating radar surveys, complemented by extensive ground truth observations. We found that SWE can be highly variable (40% difference) over short spatial scales (tens to hundreds of meters), especially in the ablation zone where the underlying ice surfaces are typically rough. Elevation provides the dominant basin-scale influence on SWE, with gradients ranging from 115 to 400 mm/100 m. Regionally, total accumulation and the accumulation gradient are strongly controlled by a glacier's distance from the coastal moisture source. Multiple linear regressions, used to calculate distributed SWE fields, show that robust results require adequate sampling of the true distribution of multiple terrain parameters. Final SWE estimates (comparable to winter balances) show reasonable agreement with both the Parameter-elevation Relationships on Independent Slopes Model climate data set (9–36% difference) and the U.S. Geological Survey Alaska Benchmark Glaciers (6–36% difference). All the glaciers in our study exhibit substantial sensitivity to changing snow-rain fractions, regardless of their location in a coastal or continental climate. While process-based SWE projections remain elusive, the collection of ground-penetrating radar (GPR)-derived data sets provides a greatly enhanced perspective on the spatial distribution of SWE and will pave the way for future work that may eventually allow such projections.

Alaska↗
USGS2014

Assessing streamflow sensitivity to variations in glacier mass balance

The mountains ringing the Gulf of Alaska (GOA) receive upwards of 4–8 m yr −1 of precipitation (Simpson et al.2005; Weingartner et al. 2005; O’Neel 2012), much of which runs off into productive coastal waters. The alpine landscape is heavily glacierized, and storage and turnover of water by glaciers substantially influences the regional surface water balance (Neal et al. 2010). In turn, the land-to-ocean flux of freshwater impacts the biogeochemistry, physical oceanography, freshwater and marine ecology of the downstream components of the GOA ecosystem (e.g., Royer et al. 2001; Hood and Scott 2008). In this way, the links between terrestrial and ocean ecosystems along the GOA have widespread impacts on regional socioeconomic issues including water and hydropower resources, fish populations, and sea level change (Dorava and Milner 2000; Royer and Grosch 2006; Cherry et al. 2010; Gardner et al. 2013). Moreover, predicting future changes in physical, chemical and biological processes in near-shore ecosystems along the GOA hinges, in part, on developing a robust understanding of water storage and transfer by glaciers through streams to the ocean. Glacierized basins (i.e. presently ice covered as opposed to glaciated, or historically ice covered) are very efficient producers of runoff, yielding 2–10 times greater runoff than similarly sized, non-glacierized basins (Mayo 1984). The unique energy balance that characterizes these basins (Jansson et al. 2003; Hock 2005) results in substantial alterations to streamflow, even when fractional ice coverage is very small (Stahl and Moore 2006). Consistent and precise treatment of glacier runoff is key to accurate assessment of hydrologic, ecological and socioeconomic impacts, but previously used definitions for glacier runoff are variable. They include: 1) meltwater produced as negative annual mass balance (e.g., Fountain and Tangborn 1985); 2) storage changes in the monthly water budget, where solid precipitation is balanced by melt and evaporation (Huss 2011, concept #2); 3) meltwater derived from melting ice only (irrespective of melting snow or mass balance) (Nolin et al. 2010; Huss 2011, concept #1); 4) all meltwater derived from the glacier surface (Cogley et al. 2011, meltwater runoff); 5) total runoff from the glacier surface (meltwater runoff plus rain on the glacier) (Neal et al. 2010). Total glacier runoff (Definitions 4 and 5 above) includes a contribution from annual mass balance, i.e. the sum of accumulation and ablation through a mass balance year (Definition 1), or what has historically been referred to as the “net” balance (Cogley et al. 2011). Indeed, annual balance has been shown to be an important driver of streamflow trends in glacierized basins, with periods of persistent negative annual balance resulting in statistically significant increases in streamflow (e.g., Pellicciotti et al. 2010). However, in maritime climates, anomalies in glacier runoff can be disconnected from annual balance because of the high variability in winter precipitation. For example, positive anomalies in winter accumulation can result in elevated levels of glacier runoff in times of positive annual mass balance (Thayyen and Gergan 2010). Quantifying the impacts of changing glacier geometries (annual balance) on glacier runoff is essential for predicting future changes in streamflow in glacierized basins. However, determining the role that this component plays in total glacier runoff (Definition 5) requires consistent measurements of seasonal (or shorter period) mass balances, measurements of precipitation at multiple locations within a basin, and streamflow measurements in close proximity to a glacier’s terminus. Practical and logistical challenges associated with assembling such data sets typically preclude such partitioning. As a result, most analyses of the relationship between annual mass balance and streamflow rely on some component of model output to compute glacier runoff (e.g. Huss et al. 2008; Kaser et al. 2010). Ultimately, developing an understanding of how total glacier runoff will change in the future is critical for predicting downstream ecological impacts associated with changes in riverine fluxes of water, sediment, and solutes (e.g., metals and nutrients) to near-shore coastal ecosystems. The purpose of this paper is to evaluate relationships among seasonal and annual glacier mass balances, glacier runoff and streamflow in two glacierized basins in different climate settings. We use long-term glacier mass balance and streamflow datasets from the United States Geological Survey (USGS) Alaska Benchmark Glacier Program to compare and contrast glacier-streamflow interactions in a maritime climate (Wolverine Glacier) with those in a continental climate (Gulkana Glacier). Our overall goal is to improve our understanding of how glacier mass balance processes impact streamflow, ultimately improving our conceptual understanding of the future evolution of glacier runoff in continental and maritime climates.

Alaska↗