scholarly journals Surface Mass-Balance Gradients From Elevation and Ice Flux Data in the Columbia Basin, Canada

2021 ◽  
Vol 9 ◽  
Author(s):  
Ben M. Pelto ◽  
Brian Menounos

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 ± 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 ± 0.07 m w.e.) between observed and modeled mass balance over flux bins is relatively small, yet mean absolute error (0.55 ± 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.

2016 ◽  
Vol 62 (236) ◽  
pp. 1083-1092 ◽  
Author(s):  
SHUN TSUTAKI ◽  
SHIN SUGIYAMA ◽  
DAIKI SAKAKIBARA ◽  
TAKANOBU SAWAGAKI

ABSTRACTTo quantify recent thinning of marine-terminating outlet glaciers in northwestern Greenland, we carried out field and satellite observations near the terminus of Bowdoin Glacier. These data were used to compute the change in surface elevation from 2007 to 2013 and this rate of thinning was then compared with that of the adjacent land-terminating Tugto Glacier. Comparing DEMs of 2007 and 2010 shows that Bowdoin Glacier is thinning more rapidly (4.1 ± 0.3 m a−1) than Tugto Glacier (2.8 ± 0.3 m a−1). The observed negative surface mass-balance accounts for <40% of the elevation change of Bowdoin Glacier, meaning that the thinning of Bowdoin Glacier cannot be attributable to surface melting alone. The ice speed of Bowdoin Glacier increases down-glacier, reaching 457 m a−1 near the calving front. This flow regime causes longitudinal stretching and vertical compression at a rate of −0.04 a−1. It is likely that this dynamically-controlled thinning has been enhanced by the acceleration of the glacier since 2000. Our measurements indicate that ice dynamics indeed play a predominant role in the rapid thinning of Bowdoin Glacier.


2002 ◽  
Vol 35 ◽  
pp. 67-72 ◽  
Author(s):  
Edward Hanna ◽  
Philippe Huybrechts ◽  
Thomas L. Mote

AbstractWe used surface climate fields from high-resolution (~0.5660.56˚) European Centre for Medium-RangeWeather Forecasts (ECMWF) operational analyses (1992–98), together with meteorological and glaciological models of snow accumulation and surface meltwater runoff/retention, to produce novel maps of Greenland ice sheet (GIS) net accumulation, net runoff and surface mass balance (SMB). We compared our runoff maps with similar-scaled runoff (melt minus refreezing) maps based on passive-microwave satellite data. Our gross spatial/temporal patterns of runoff compared well with those from the satellite data, although amounts of modelled runoff are likely too low. Mean accumulation was 0.287 (0.307)ma–1, and mean runoff was 0.128 (0.151)ma–1, averaged across the W. Abdalati (T. L. Mote) GIS mask. Corresponding mean SMB was 0.159 (0.156)ma–1, with considerable interannual variability (standard deviation ~0.11ma–1) primarily due to variations in runoff. Considering best estimates of current iceberg calving, overall the GIS is probably currently losing mass. Our study shows great promise for meaningfully modelling SMB based on forthcoming ``second-generation’’ ECMWF re-analysis (ERA-40) data, and comparing the results with ongoing laser/radarmeasurements of surface elevation. This should help elucidate to what extent surface elevation changes are caused by short-term SMB variations or other factors (e.g. ice dynamics).


1977 ◽  
Vol 18 (80) ◽  
pp. 359-371 ◽  
Author(s):  
I. M. Whillans

Abstract The continuity relationship that is often used in the study of ice sheets and ice shelves is developed by integrating the equation of continuity through the ice thickness. This equation is then integrated again with respect to horizontal distance from an ice divide, showing that the difference between the true ice velocity and the balance velocity, which is defined, is a measure of the time chance of the mass of a column through the ice thickness. The relationship is applied using data from along the “Byrd” station strain network, Antarctica. This region is found to be thinning slowly (0.03 m a−1 of ice of mean density) and uniformly, but it is still close to steady-state. The calculations would show a larger thinning rate if bottom sliding contributed more to the ice movement and integral shear contributed less, but the “Byrd” station bore-hole tilting results of Garfield and Ueda (1975, 1976), together with surface velocity measurements at “Byrd” station, indicate that most of the ice flow is by deformation within the ice mass. This large amount of internal deformation is more than that predicted by most “flow laws”, probably because of the strongly oriented ice-crystal fabric in the ice sheet. The cause of ice thinning is probably decreased surface mass balance beginning before A.D. 1550. The consistent relationship between measured velocity and balance velocity indicates that the ice flow is simple and that flow lines are in the same direction at depth as at the surface when considered smoothed over a distance of 10 km. Because the ice sheet is at present thinning, the balance velocity, calculated only from flow line and surface mass-balance data, and the somewhat mistaken assumption of steady-state is 15% less than the true ice velocity. This rather small difference confirms the use of balance-velocity estimates where velocity measurements are not available.


2018 ◽  
Vol 64 (246) ◽  
pp. 568-582 ◽  
Author(s):  
GABRIELA COLLAO-BARRIOS ◽  
FABIEN GILLET-CHAULET ◽  
VINCENT FAVIER ◽  
GINO CASASSA ◽  
ETIENNE BERTHIER ◽  
...  

ABSTRACTWe simulate the ice dynamics of the San Rafael Glacier (SRG) in the Northern Patagonia Icefield (46.7°S, 73.5°W), using glacier geometry obtained by airborne gravity measurements. The full-Stokes ice flow model (Elmer/Ice) is initialized using an inverse method to infer the basal friction coefficient from a satellite-derived surface velocity mosaic. The high surface velocities (7.6 km a−1) near the glacier front are explained by low basal shear stresses (<25 kPa). The modelling results suggest that 98% of the surface velocities are due to basal sliding in the fast-flowing glacier tongue (>1 km a−1). We force the model using different surface mass-balance scenarios taken or adapted from previous studies and geodetic elevation changes between 2000 and 2012. Our results suggest that previous estimates of average surface mass balance over the entire glacier (Ḃ) were likely too high, mainly due to an overestimation in the accumulation area. We propose that most of SRG imbalance is due to the large ice discharge (−0.83 ± 0.08 Gt a−1) and a slightly positiveḂ(0.08 ± 0.06 Gt a−1). The committed mass-loss estimate over the next century is −0.34 ± 0.03 Gt a−1. This study demonstrates that surface mass-balance estimates and glacier wastage projections can be improved using a physically based ice flow model.


2015 ◽  
Vol 9 (6) ◽  
pp. 2009-2025 ◽  
Author(s):  
P. Kuipers Munneke ◽  
S. R. M. Ligtenberg ◽  
B. P. Y. Noël ◽  
I. M. Howat ◽  
J. E. Box ◽  
...  

Abstract. Observed changes in the surface elevation of the Greenland Ice Sheet are caused by ice dynamics, basal elevation change, basal melt, surface mass balance (SMB) variability, and by compaction of the overlying firn. The last two contributions are quantified here using a firn model that includes compaction, meltwater percolation, and refreezing. The model is forced with surface mass fluxes and temperature from a regional climate model for the period 1960–2014. The model results agree with observations of surface density, density profiles from 62 firn cores, and altimetric observations from regions where ice-dynamical surface height changes are likely small. In areas with strong surface melt, the firn model overestimates density. We find that the firn layer in the high interior is generally thickening slowly (1–5 cm yr−1). In the percolation and ablation areas, firn and SMB processes account for a surface elevation lowering of up to 20–50 cm yr−1. Most of this firn-induced marginal thinning is caused by an increase in melt since the mid-1990s and partly compensated by an increase in the accumulation of fresh snow around most of the ice sheet. The total firn and ice volume change between 1980 and 2014 is estimated at −3295 ± 1030 km3 due to firn and SMB changes, corresponding to an ice-sheet average thinning of 1.96 ± 0.61 m. Most of this volume decrease occurred after 1995. The computed changes in surface elevation can be used to partition altimetrically observed volume change into surface mass balance and ice-dynamically related mass changes.


2021 ◽  
Vol 2 (1) ◽  
Author(s):  
Youngmin Choi ◽  
Mathieu Morlighem ◽  
Eric Rignot ◽  
Michael Wood

AbstractThe mass loss of the Greenland Ice Sheet is nearly equally partitioned between a decrease in surface mass balance from enhanced surface melt and an increase in ice dynamics from the acceleration and retreat of its marine-terminating glaciers. Much uncertainty remains in the future mass loss of the Greenland Ice Sheet due to the challenges of capturing the ice dynamic response to climate change in numerical models. Here, we estimate the sea level contribution of the Greenland Ice Sheet over the 21st century using an ice-sheet wide, high-resolution, ice-ocean numerical model that includes surface mass balance forcing, thermal forcing from the ocean, and iceberg calving dynamics. The model is calibrated with ice front observations from the past eleven years to capture the recent evolution of marine-terminating glaciers. Under a business as usual scenario, we find that northwest and central west Greenland glaciers will contribute more mass loss than other regions due to ice front retreat and ice flow acceleration. By the end of century, ice discharge from marine-terminating glaciers will contribute 50 ± 20% of the total mass loss, or twice as much as previously estimated although the contribution from the surface mass balance increases towards the end of the century.


1977 ◽  
Vol 18 (80) ◽  
pp. 359-371 ◽  
Author(s):  
I. M. Whillans

AbstractThe continuity relationship that is often used in the study of ice sheets and ice shelves is developed by integrating the equation of continuity through the ice thickness. This equation is then integrated again with respect to horizontal distance from an ice divide, showing that the difference between the true ice velocity and the balance velocity, which is defined, is a measure of the time chance of the mass of a column through the ice thickness.The relationship is applied using data from along the “Byrd” station strain network, Antarctica. This region is found to be thinning slowly (0.03 m a−1 of ice of mean density) and uniformly, but it is still close to steady-state. The calculations would show a larger thinning rate if bottom sliding contributed more to the ice movement and integral shear contributed less, but the “Byrd” station bore-hole tilting results of Garfield and Ueda (1975, 1976), together with surface velocity measurements at “Byrd” station, indicate that most of the ice flow is by deformation within the ice mass. This large amount of internal deformation is more than that predicted by most “flow laws”, probably because of the strongly oriented ice-crystal fabric in the ice sheet. The cause of ice thinning is probably decreased surface mass balance beginning before A.D. 1550.The consistent relationship between measured velocity and balance velocity indicates that the ice flow is simple and that flow lines are in the same direction at depth as at the surface when considered smoothed over a distance of 10 km. Because the ice sheet is at present thinning, the balance velocity, calculated only from flow line and surface mass-balance data, and the somewhat mistaken assumption of steady-state is 15% less than the true ice velocity. This rather small difference confirms the use of balance-velocity estimates where velocity measurements are not available.


2015 ◽  
Vol 9 (3) ◽  
pp. 3541-3580 ◽  
Author(s):  
P. Kuipers Munneke ◽  
S. R. M. Ligtenberg ◽  
B. P. Y. Noël ◽  
I. M. Howat ◽  
J. E. Box ◽  
...  

Abstract. Observed changes in the surface elevation of the Greenland ice sheet are caused by ice dynamics, basal elevation change, surface mass balance (SMB) variability, and by compaction of the overlying firn. The latter two contributions are quantified here using a firn model that includes compaction, meltwater percolation, and refreezing. The model is forced with surface mass fluxes and temperature from a regional climate model for the period 1960–2013. The model results agree with observations of surface density, density profiles from 62 firn cores, and altimetric observations from regions where ice-dynamical surface height changes are likely small. We find that the firn layer in the high interior is generally thickening slowly (1–5 cm yr−1). In the percolation and ablation areas, firn and SMB processes account for a surface elevation lowering of up to 20–50 cm yr−1. Most of this firn-induced marginal thinning is caused by an increase in melt since the mid-1990s, and partly compensated by an increase in the accumulation of fresh snow around most of the ice sheet. The total firn and ice volume change between 1980 and 2013 is estimated at −3900 ± 1030 km3 due to firn and SMB, corresponding to an ice-sheet average thinning of 2.32 ± 0.61 m. Most of this volume decrease occurred after 1995. The computed changes in surface elevation can be used to partition altimetrically observed volume change into surface mass balance and ice-dynamically related mass changes.


2011 ◽  
Vol 5 (1) ◽  
pp. 299-313 ◽  
Author(s):  
G. E. Flowers ◽  
N. Roux ◽  
S. Pimentel ◽  
C. G. Schoof

Abstract. Glacier surges are a well-known example of an internal dynamic oscillation whose occurrence is not a direct response to the external climate forcing, but whose character (i.e. period, amplitude, mechanism) may depend on the glacier's environmental or climate setting. We examine the dynamics of a small (∼5 km2) valley glacier in Yukon, Canada, where two previous surges have been photographically documented and an unusually slow surge is currently underway. To characterize the dynamics of the present surge, and to speculate on the future of this glacier, we employ a higher-order flowband model of ice dynamics with a regularized Coulomb-friction sliding law in both diagnostic and prognostic simulations. Diagnostic (force balance) calculations capture the measured ice-surface velocity profile only when non-zero basal water pressures are prescribed over the central region of the glacier, coincident with where evidence of the surge has been identified. This leads to sliding accounting for 50–100% of the total surface motion in this region. Prognostic simulations, where the glacier geometry evolves in response to a prescribed surface mass balance, reveal a significant role played by a bedrock ridge beneath the current equilibrium line of the glacier. Ice thickening occurs above the ridge in our simulations, until the net mass balance reaches sufficiently negative values. We suggest that the bedrock ridge may contribute to the propensity for surges in this glacier by promoting the development of the reservoir area during quiescence, and may permit surges to occur under more negative balance conditions than would otherwise be possible. Collectively, these results corroborate our interpretation of the current glacier flow regime as indicative of a slow surge that has been ongoing for some time, and support a relationship between surge incidence or character and the net mass balance. Our results also highlight the importance of glacier bed topography in controlling ice dynamics, as observed in many other glacier systems.


2021 ◽  
Author(s):  
Lilian Schuster ◽  
David Rounce ◽  
Fabien Maussion

&lt;p&gt;A recent large model intercomparison study (GlacierMIP) showed that differences between the glacier models is a dominant source of uncertainty for future glacier change projections, in particular in the first half of the century.&amp;#160; Each glacier model has their own unique set of process representations and climate forcing methodology, which makes it impossible to determine the model components that contribute most to the projection uncertainty. This study aims to improve our understanding of the sources of large scale glacier model uncertainty using the Open Global Glacier Model (OGGM), focussing on the surface mass balance (SMB) in a first step. We calibrate and run a set of interchangeable SMB model parameterizations (e.g. monthly vs. daily, constant vs. variable lapse rates, albedo, snowpack evolution and refreezing) under controlled boundary conditions. Based on ensemble approaches, we explore the influence of (i) the parameter calibration strategy and (ii) SMB model complexity on regional to global glacier change. These uncertainties are then put in relation to a qualitative selection of other model design choices, such as the forcing climate dataset and ice dynamics model parameters.&amp;#160;&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document