The mass-balance—elevation relation for a given glacier is required for many numerical models of ice flow. Direct measurements of this relation using remotely-sensed methods are complicated by ice dynamics, so observations are currently limited to glaciers where surface mass-balance measurements are routinely made. We test the viability of using the continuity equation to estimate annual surface mass balance between flux-gates in the absence of in situ measurements, on five glaciers in the Columbia Mountains of British Columbia, Canada. Repeat airborne laser scanning surveys of glacier surface elevation, ice penetrating radar surveys and publicly available maps of ice thickness are used to estimate changes in surface elevation and ice flux. We evaluate this approach by comparing modeled to observed mass balance. Modeled mass-balance gradients well-approximate those obtained from direct measurement of surface mass balance, with a mean difference of +6.6 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="m1"><mml:mo>±</mml:mo></mml:math> 4%. Regressing modeled mass balance, equilibrium line altitudes are on average 15 m higher than satellite-observations of the transient snow line. Estimates of mass balance over flux bins compare less favorably than the gradients. Average mean error (+0.03 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="m2"><mml:mo>±</mml:mo></mml:math> 0.07 m w.e.) between observed and modeled mass balance over flux bins is relatively small, yet mean absolute error (0.55 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="m3"><mml:mo>±</mml:mo></mml:math> 0.18 m w.e.) and average modeled mass-balance uncertainty (0.57 m w.e.) are large. Mass conservation, assessed with glaciological data, is respected (when estimates are within 1σ uncertainties) for 84% of flux bins representing 86% of total glacier area. Uncertainty on ice velocity, especially for areas where surface velocity is low (<10 m a −1 ) contributes the greatest error in estimating ice flux. We find that using modeled ice thicknesses yields comparable modeled mass-balance gradients relative to using observations of ice thickness, but we caution against over-interpreting individual flux-bin mass balances due to large mass-balance residuals. Given the performance of modeled ice thickness and the increasing availability of ice velocity and surface topography data, we suggest that similar efforts to produce mass-balance gradients using modern high-resolution datasets are feasible on larger scales.

The mass-balance—elevation relation for a given glacier is required for many numerical models of ice flow. Direct measurements of this relation using remotely-sensed methods are complicated by ice dynamics, so observations are currently limited to glaciers where surface mass-balance measurements are routinely made. We test the viability of using the continuity equation to estimate annual surface mass balance between flux-gates in the absence of in situ measurements, on five glaciers in the Columbia Mountains of British Columbia, Canada. Repeat airborne laser scanning surveys of glacier surface elevation, ice penetrating radar surveys and publicly available maps of ice thickness are used to estimate changes in surface elevation and ice flux. We evaluate this approach by comparing modeled to observed mass balance. Modeled mass-balance gradients well-approximate those obtained from direct measurement of surface mass balance, with a mean difference of +6.6 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="m1"><mml:mo>±</mml:mo></mml:math> 4%. Regressing modeled mass balance, equilibrium line altitudes are on average 15 m higher than satellite-observations of the transient snow line. Estimates of mass balance over flux bins compare less favorably than the gradients. Average mean error (+0.03 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="m2"><mml:mo>±</mml:mo></mml:math> 0.07 m w.e.) between observed and modeled mass balance over flux bins is relatively small, yet mean absolute error (0.55 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="m3"><mml:mo>±</mml:mo></mml:math> 0.18 m w.e.) and average modeled mass-balance uncertainty (0.57 m w.e.) are large. Mass conservation, assessed with glaciological data, is respected (when estimates are within 1σ uncertainties) for 84% of flux bins representing 86% of total glacier area. Uncertainty on ice velocity, especially for areas where surface velocity is low (<10 m a −1 ) contributes the greatest error in estimating ice flux. We find that using modeled ice thicknesses yields comparable modeled mass-balance gradients relative to using observations of ice thickness, but we caution against over-interpreting individual flux-bin mass balances due to large mass-balance residuals. Given the performance of modeled ice thickness and the increasing availability of ice velocity and surface topography data, we suggest that similar efforts to produce mass-balance gradients using modern high-resolution datasets are feasible on larger scales.

Western North American (WNA) glaciers outside of Alaska cover 14,384 km2 of mountainous terrain. No comprehensive analysis of recent mass change exists for this region. We generated over 15,000 multisensor digital elevation models from spaceborne optical imagery to provide an assessment of mass change for WNA over the period 2000–2018. These glaciers lost 117 ± 42 gigatons (Gt) of mass, which accounts for up to 0.32 ± 0.11 mm of sea level rise over the full period of study. We observe a fourfold increase in mass loss rates between 2000–2009 [−2.9 ± 3.1 Gt yr−1] and 2009–2018 [−12.3 ± 4.6 Gt yr−1], and we attribute this change to a shift in regional meteorological conditions driven by the location and strength of upper level zonal wind. Our results document decadal‐scale climate variability over WNA that will likely modulate glacier mass change in the future.

Abstract. Seasonal measurements of glacier mass balance provide insight into the relation between climate forcing and glacier change. To evaluate the feasibility of using remotely sensed methods to assess seasonal balance, we completed tandem airborne laser scanning (ALS) surveys and field-based glaciological measurements over a 4-year period for six alpine glaciers that lie in the Columbia and Rocky Mountains, near the headwaters of the Columbia River, British Columbia, Canada. We calculated annual geodetic balance using coregistered late summer digital elevation models (DEMs) and distributed estimates of density based on surface classification of ice, snow, and firn surfaces. Winter balance was derived using coregistered late summer and spring DEMs, as well as density measurements from regional snow survey observations and our glaciological measurements. Geodetic summer balance was calculated as the difference between winter and annual balance. Winter mass balance from our glaciological observations averaged 1.95±0.09 m w.e. (meter water equivalent), 4 % larger than those derived from geodetic surveys. Average glaciological summer and annual balance were 3 % smaller and 3 % larger, respectively, than our geodetic estimates. We find that distributing snow, firn, and ice density based on surface classification has a greater influence on geodetic annual mass change than the density values themselves. Our results demonstrate that accurate assessments of seasonal mass change can be produced using ALS over a series of glaciers spanning several mountain ranges. Such agreement over multiple seasons, years, and glaciers demonstrates the ability of high-resolution geodetic methods to increase the number of glaciers where seasonal mass balance can be reliably estimated.