Basal Mass Balance Research Articles

Compensating errors in inversions for subglacial bed roughness: same steady state, different dynamic response

Abstract. Subglacial bed roughness is one of the main factors controlling the rate of future Antarctic ice-sheet retreat and also one of the most uncertain. A common technique to constrain the bed roughness using ice-sheet models is basal inversion, tuning the roughness to reproduce the observed present-day ice-sheet geometry and/or surface velocity. However, many other factors affecting ice-sheet evolution, such as the englacial temperature and viscosity, the surface and basal mass balance, and the subglacial topography, also contain substantial uncertainties. Using a basal inversion technique intrinsically causes any errors in these other quantities to lead to compensating errors in the inverted bed roughness. Using a set of idealised-geometry experiments, we quantify these compensating errors and investigate their effect on the dynamic response of the ice sheet to a prescribed forcing. We find that relatively small errors in ice viscosity and subglacial topography require substantial compensating errors in the bed roughness in order to produce the same steady-state ice sheet, obscuring the realistic spatial variability in the bed roughness. When subjected to a retreat-inducing forcing, we find that these different parameter combinations, which per definition of the inversion procedure result in the same steady-state geometry, lead to a rate of ice volume loss that can differ by as much as a factor of 2. This implies that ice-sheet models that use basal inversion to initialise their model state can still display a substantial model bias despite having an initial state which is close to the observations.

Basal mass balance and prevalence of ice tongues in the Western ross sea

Ice tongues at the fringes of the Antarctic ice sheet lose mass primarily through both basal melting and calving. They are sensitive to ocean conditions which can weaken the ice both mechanically or through thinning. Ice tongues, which are laterally unconfined, are likely to be particularly sensitive to ocean-induced stresses. Here we examine ice tongues in the Western Ross Sea, by looking into the factors affecting their stability. We calculate the basal mass change of twelve Antarctic ice tongues using a flux gate approach, deriving thickness from ICESat-2 height measurements and ice surface velocities from Sentinel-1 feature-tracking over the same period (October 2018 to December 2021). The basal mass balance ranges between −0.14 ± 0.07 m yr−1 and −1.50 ± 1.2 m yr−1. The average basal mass change for all the ice tongues is −0.82 ± 0.68 m of ice yr−1. Low values of basal melt suggest a stable mass balance condition in this region, with low thermal ocean forcing, as other studies have shown. We found a heterogeneous basal melt pattern with no latitudinal gradient and no clear driver in basal melt indicating that local variables are important in the persistence of ice tongues in the absence of a strong oceanographic melting force. Moreover, thanks to the temporal resolution of the data we were able to resolve the seasonal variability of Drygalski and Aviator Ice Tongues, the two largest ice tongues studied.

Greenland-Ice-Sheet Surface Temperature and Melt Extent from 2000 to 2020 and Implications for Mass Balance

Accurate monitoring of surface temperature and melting on the Greenland Ice Sheet (GrIS) is important for tracking the ice sheet’s mass balance as well as global and Arctic climate change. Using a moderate-resolution-imaging-spectroradiometer (MODIS)-derived land-surface-temperature (LST) data product with a resolution of 1 km from 2000 to 2020, the temporal and spatial variations of annual and seasonal ‘clear-sky’ surface temperature were evaluated. We also monitored summer surface melting and studied the relationship between the mass balance of the ice sheet and changes in surface temperature and melting. The results show that the mean annual LST during the study period is −24.86 ± 5.46 °C, with the highest of −22.48 ± 5.61 °C in 2010 and the lowest temperature of −26.49 ± 5.30 °C in 2015. With the change of season, the spatial variation of the ice-sheet surface temperature changes greatly. 2012 and 2019 experienced the warmest summers (−5.92 ± 4.01 °C and −6.51 ± 3.93 °C), with extreme cumulative melting detected on the ice-sheet surface (89.9% and 89.7%, respectively), and 2002 also experienced a greater extent of melting. But short period of melt in 2002 and 2019 (30.6% and 31.4%, respectively), accounted for a larger proportion, with neither the duration nor intensity of the melt reaching that of 2012. There is a strong correlation between the GrIS surface temperature and its mass balance. By fitting the relationship between surface temperature and mass balance, it was found that 93.83% (6.17%) of the ice-sheet response to surface-temperature change was via surface-mass balance (discharge and basal-mass balance).

Thermal structure of the Amery Ice Shelf from borehole observations and simulations

Abstract. The Amery Ice Shelf (AIS), East Antarctica, has a layered structure, due to the presence of both meteoric and marine ice. In this study, the thermal structure of the AIS and its spatial pattern are evaluated and analysed through borehole observations and numerical simulations with Elmer/Ice, a full-Stokes ice sheet/shelf model. In the area with marine ice, a near-isothermal basal layer up to 120 m thick is observed, which closely conforms to the pressure-dependent freezing temperature of seawater. In the area experiencing basal melting, large temperature gradients, up to −0.36 ∘C m−1, are observed at the base. Three-dimensional (3-D) steady-state temperature simulations with four different basal mass balance (BMB) datasets for the AIS reveal a high sensitivity of ice shelf thermal structure to the distribution of BMB. We also construct a one-dimensional (1-D) transient temperature column model to simulate the process of an ice column moving along a flowline with corresponding boundary conditions, which achieves slightly better agreement with borehole observations than the 3-D simulations. Our simulations reveal internal cold ice advected from higher elevations by the AIS's main inlet glaciers, warming downstream along the ice flow, and we suggest the thermal structures dominated by these cold cores may commonly exist among Antarctic ice shelves. For the marine ice, the porous structure of its lower layer and interactions with ocean below determine the local thermal regime and give rise to the near-isothermal phenomenon. The limitations in our simulations identify the need for ice shelf–ocean coupled models with improved thermodynamics and more comprehensive boundary conditions. Given the temperature dependence of ice rheology, the depth-averaged ice stiffness factor B(T′)‾ derived from the most realistic simulated temperature field is presented to quantify the influence of the temperature distribution on ice shelf dynamics. The full 3-D temperature field provides a useful input to future modelling studies.

Greenland ice sheet mass balance from 1840 through next week

Abstract. The mass of the Greenland ice sheet is declining as mass gain from snow accumulation is exceeded by mass loss from surface meltwater runoff, marine-terminating glacier calving and submarine melting, and basal melting. Here we use the input–output (IO) method to estimate mass change from 1840 through next week. Surface mass balance (SMB) gains and losses come from a semi-empirical SMB model from 1840 through 1985 and three regional climate models (RCMs; HIRHAM/HARMONIE, Modèle Atmosphérique Régional – MAR, and RACMO – Regional Atmospheric Climate MOdel) from 1986 through next week. Additional non-SMB losses come from a marine-terminating glacier ice discharge product and a basal mass balance model. From these products we provide an annual estimate of Greenland ice sheet mass balance from 1840 through 1985 and a daily estimate at sector and region scale from 1986 through next week. This product updates daily and is the first IO product to include the basal mass balance which is a source of an additional ∼24 Gt yr−1 of mass loss. Our results demonstrate an accelerating ice-sheet-scale mass loss and general agreement (coefficient of determination, r2, ranges from 0.62 to 0.94) among six other products, including gravitational, volume, and other IO mass balance estimates. Results from this study are available at https://doi.org/10.22008/FK2/OHI23Z (Mankoff et al., 2021).

Coupling the U.K. Earth System Model to Dynamic Models of the Greenland and Antarctic Ice Sheets

AbstractThe physical interactions between ice sheets and the atmosphere and ocean around them are major factors in determining the state of the climate system, yet many current Earth System models omit them entirely or treat them very simply. In this work we describe how models of the Greenland and Antarctic ice sheets have been incorporated into the global U.K. Earth System model (UKESM1) via substantial technical developments with a two‐way coupling that passes fluxes of energy and water, and the topography of the ice sheet surface and ice shelf base, between the component models. File‐based coupling outside the running model executables is used throughout to pass information between the components, which we show is both physically appropriate and convenient within the UKESM1 structure. Ice sheet surface mass balance is computed in the land surface model using multi‐layer snowpacks in subgrid‐scale elevation ranges and compares well to the results of regional climate models. Ice shelf front discharge forms icebergs, which drift and melt in the ocean. Ice shelf basal mass balance is simulated using the full three‐dimensional ocean model representation of the circulation in ice‐shelf cavities. We show a range of example results, including from simulations with changes in ice sheet height and thickness of hundreds of meters, and changes in ice sheet grounding line and land‐terminating margin of many tens of kilometres, demonstrating that the coupled model is computationally stable when subject to significant changes in ice sheet geometry.

Mapping Basal Melt Under the Shackleton Ice Shelf, East Antarctica, From CryoSat-2 Radar Altimetry

Ice shelvesbuttress the fast-flowing glaciers draining ice sheets, so changes in the ice shelves may alter the buttressing effect and the discharge of grounded ice, therefore influencing sea level. Basal melting of the ice shelves would reduce their ability to restrain the discharge of grounded ice into the ocean. However, estimating the basal melt rate is challenging due to the large uncertainty in the calculation. In this study, we use a Lagrangian framework to improve the basal melt rate derivation and apply the method to the Shackleton ice shelf as a case study. We use the CryoSat-2 data to characterize the spatial distribution patterns of ice shelf surface elevation changes between 2010 and 2018. Combining these results with the ice surface velocity and output from the regional climate model, we obtain a map of the basal melt rate and calculate the total and average basal mass balance. The total basal meltwater production for the Shackleton ice shelf is 54.6 ± 7.2 Gt/yr. The highest melt rates, which exceed 50 m/yr, are found close to the grounding line and the main trunk of Denman glacier. Based on the analysis, we show that the Lagrangian method can provide more spatially coherent patterns of ice shelf surface elevation changes and reduce the uncertainty in basal melt rate calculation.

Variable Basal Melt Rates of Antarctic Peninsula Ice Shelves, 1994–2016

AbstractWe have constructed 23‐year (1994–2016) time series of Antarctic Peninsula (AP) ice‐shelf height change using data from four satellite radar altimeters (ERS‐1, ERS‐2, Envisat, and CryoSat‐2). Combining these time series with output from atmospheric and firn models, we partitioned the total height‐change signal into contributions from varying surface mass balance, firn state, ice dynamics, and basal mass balance. On the Bellingshausen coast of the AP, ice shelves lost 84 ± 34 Gt a−1 to basal melting, compared to contributions of 50 ± 7 Gt a−1 from surface mass balance and ice dynamics. Net basal melting on the Weddell coast was 51 ± 71 Gt a−1. Recent changes in ice‐shelf height include increases over major AP ice shelves driven by changes in firn state. Basal melt rates near Bawden Ice Rise, a major pinning point of Larsen C Ice Shelf, showed large increases, potentially leading to substantial loss of buttressing if sustained.

Continent-wide estimates of Antarctic strain rates from Landsat 8-derived velocity grids

ABSTRACTStrain rates are fundamental measures of ice flow and are used in a wide variety of glaciological applications including investigations of bed properties, calculations of basal mass balance on ice shelves, and constraints on ice rheological models. However, despite their extensive application, strain rates are calculated using a variety of methods and length scales and the details are often not specified. In this study, we compare the results of nominal and logarithmic strain-rate calculations based on a satellite-derived velocity field of the Antarctic ice sheet generated from Landsat 8 satellite data. Our comparison highlights the differences between the two common approaches in the glaciological literature. We evaluate the errors introduced by each approach and their impacts on the results. We also demonstrate the importance of choosing and specifying a length scale over which strain-rate calculations are made, which can strongly influence other derived quantities such as basal mass balance on ice shelves. Finally, we present strain-rate data products calculated using an approximate viscous length-scale with satellite observations of ice velocity for the Antarctic continent.

Detecting high spatial variability of ice shelf basal mass balance, Roi Baudouin Ice Shelf, Antarctica

Abstract. Ice shelves control the dynamic mass loss of ice sheets through buttressing and their integrity depends on the spatial variability of their basal mass balance (BMB), i.e. the difference between refreezing and melting. Here, we present an improved technique – based on satellite observations – to capture the small-scale variability in the BMB of ice shelves. As a case study, we apply the methodology to the Roi Baudouin Ice Shelf, Dronning Maud Land, East Antarctica, and derive its yearly averaged BMB at 10 m horizontal gridding. We use mass conservation in a Lagrangian framework based on high-resolution surface velocities, atmospheric-model surface mass balance and hydrostatic ice-thickness fields (derived from TanDEM-X surface elevation). Spatial derivatives are implemented using the total-variation differentiation, which preserves abrupt changes in flow velocities and their spatial gradients. Such changes may reflect a dynamic response to localized basal melting and should be included in the mass budget. Our BMB field exhibits much spatial detail and ranges from −14.7 to 8.6 m a−1 ice equivalent. Highest melt rates are found close to the grounding line where the pressure melting point is high, and the ice shelf slope is steep. The BMB field agrees well with on-site measurements from phase-sensitive radar, although independent radar profiling indicates unresolved spatial variations in firn density. We show that an elliptical surface depression (10 m deep and with an extent of 0.7 km × 1.3 km) lowers by 0.5 to 1.4 m a−1, which we tentatively attribute to a transient adaptation to hydrostatic equilibrium. We find evidence for elevated melting beneath ice shelf channels (with melting being concentrated on the channel's flanks). However, farther downstream from the grounding line, the majority of ice shelf channels advect passively (i.e. no melting nor refreezing) toward the ice shelf front. Although the absolute, satellite-based BMB values remain uncertain, we have high confidence in the spatial variability on sub-kilometre scales. This study highlights expected challenges for a full coupling between ice and ocean models.

GPS-derived estimates of surface mass balance and ocean-induced basal melt for Pine Island Glacier ice shelf, Antarctica

Abstract. In the last 2 decades, Pine Island Glacier (PIG) experienced marked speedup, thinning, and grounding-line retreat, likely due to marine ice-sheet instability and ice-shelf basal melt. To better understand these processes, we combined 2008–2010 and 2012–2014 GPS records with dynamic firn model output to constrain local surface and basal mass balance for PIG. We used GPS interferometric reflectometry to precisely measure absolute surface elevation (zsurf) and Lagrangian surface elevation change (Dzsurf∕ Dt). Observed surface elevation relative to a firn layer tracer for the initial surface (zsurf − zsurf0′) is consistent with model estimates of surface mass balance (SMB, primarily snow accumulation). A relatively abrupt ∼ 0.2–0.3 m surface elevation decrease, likely due to surface melt and increased compaction rates, is observed during a period of warm atmospheric temperatures from December 2012 to January 2013. Observed Dzsurf∕ Dt trends (−1 to −4 m yr−1) for the PIG shelf sites are all highly linear. Corresponding basal melt rate estimates range from ∼ 10 to 40 m yr−1, in good agreement with those derived from ice-bottom acoustic ranging, phase-sensitive ice-penetrating radar, and high-resolution stereo digital elevation model (DEM) records. The GPS and DEM records document higher melt rates within and near features associated with longitudinal extension (i.e., transverse surface depressions, rifts). Basal melt rates for the 2012–2014 period show limited temporal variability despite large changes in ocean temperature recorded by moorings in Pine Island Bay. Our results demonstrate the value of long-term GPS records for ice-shelf mass balance studies, with implications for the sensitivity of ice–ocean interaction at PIG.

Melting and freezing under Antarctic ice shelves from a combination of ice-sheet modelling and observations

ABSTRACTIce-shelf basal melting is the largest contributor to the negative mass balance of the Antarctic ice sheet. However, current implementations of ice/ocean interactions in ice-sheet models disagree with the distribution of sub-shelf melt and freezing rates revealed by recent observational studies. Here we present a novel combination of a continental-scale ice flow model and a calibration technique to derive the spatial distribution of basal melting and freezing rates for the whole Antarctic ice-shelf system. The modelled ice-sheet equilibrium state is evaluated against topographic and velocity observations. Our high-resolution (10-km spacing) simulation predicts an equilibrium ice-shelf basal mass balance of −1648.7 Gt a−1 that increases to −1917.0 Gt a−1 when the observed ice-shelf thinning rates are taken into account. Our estimates reproduce the complexity of the basal mass balance of Antarctic ice shelves, providing a reference for parameterisations of sub-shelf ocean/ice interactions in continental ice-sheet models. We perform a sensitivity analysis to assess the effects of variations in the model set-up, showing that the retrieved estimates of basal melting and freezing rates are largely insensitive to changes in the internal model parameters, but respond strongly to a reduction of model resolution and the uncertainty in the input datasets.

Modeling modified Circumpolar Deep Water intrusions onto the Prydz Bay continental shelf, East Antarctica

AbstractAn eddy‐resolving coupled regional ocean‐sea ice‐ice shelf model is employed to locate the hot spots where modified Circumpolar Deep Water (mCDW) intrudes onto the continental shelf within Prydz Bay, and locate the paths through which mCDW is transported to the Amery Ice Shelf (AIS) calving front. Evaluation of the model output is with satellite, hydrographic and borehole data. Two critical windows responsible for mCDW intrusions are identified. The first is the eastern branch of the cyclonic Prydz Bay gyre (PBG) that carries mCDW to the ice front line, accounting for an annual mean heat transport of ∼8.7 ×1011 J s−1. The second is located to the east of the Four Ladies Bank (FLB) where mCDW is channeled through submarine troughs, accounting for an annual mean heat transport of ∼16.2 ×1011 J s−1. The eddy‐induced heat transport accounts for ∼23% in the path of the PBG and ∼52% in the path of the eastern coastal current, with respect to their total onshore heat transport. The seasonal pulsing of mCDW intrusions is greatly dependent on the seasonal cycle of the Antarctic Slope Current (ASC) that peaks with a maximum of ∼29.3 Sv at 75°E in June. In austral winter, mCDW is allowed to access the eastern flank of the AIS calving front with potential consequences for the basal mass balance of the AIS. The dynamic effects of small‐scale troughs on the longshore ASC play an important role in the onshore mCDW transport.

Distribution of ice thickness and subglacial topography of the “Chinese Wall” around Kunlun Station, East Antarctica

As fundamental parameters of the Antarctic Ice Sheet, ice thickness and subglacial topography are critical factors for studying the basal conditions and mass balance in Antarctica. During CHINARE 24 (the 24th Chinese National Antarctic Research Expedition, 2007/08), the research team used a deep ice-penetrating radar system to measure the ice thickness and subglacial topography of the “Chinese Wall” around Kunlun Station, East Antarctica. Preliminary results show that the ice thickness varies mostly from 1600 m to 2800 m along the “Chinese Wall”, with the thickest ice being 3444 m, and the thinnest ice 1255 m. The average bedrock elevation is 1722 m, while the minimum is just 604 m. Compared with the northern side of the ice divide, the ice thickness is a little greater and the subglacial topography lower on the southern side, which is also characterized by four deep valleys. We found no basal freeze-on ice in the Gamburtsev Subglacial Mountains area, subglacial lakes, or water bodies along the “Chinese Wall”. Ice thickness and subglacial topography data extracted from the Bedmap 2 database along the “Chinese Wall” are consistent with our results, but their resolution and accuracy are very limited in areas where the bedrock fluctuates intensely. The distribution of ice thickness and subglacial topography detected by ice-penetrating radar clarifies the features of the ice sheet in this “inaccessible” region. These results will help to advance the study of ice sheet dynamics and the determination of future locations of the GSM’s geological and deep ice core drilling sites in the Dome A region.

Tide-corrected flow velocity and mass balance of Campbell Glacier Tongue, East Antarctica, derived from interferometric SAR

Accurate measurement of ice-flow velocity of floating glaciers is required to estimate ice mass balance from volume flux thinning/thickening, basal melting/freezing, and surface accumulation/ablation. We derived a tide-corrected ice velocity map and mass budget of Campbell Glacier Tongue (CGT) in East Antarctica by using 14 COSMO-SkyMed one-day tandem differential interferometric SAR (DInSAR) pairs obtained from January to November 2011. The vertical tidal deflection of CGT was estimated and removed from the DInSAR images by using a tide deflection ratio map generated by double-differential InSAR (DDInSAR) method. We then generated averaged ice-flow velocity (v map) and its standard deviation (σv map). Ice-flow velocity increased from the upper part of the grounding line (~0.20md−1) to the seaward edge of CGT (~0.67md−1) with σv less than ~0.04md−1 in the main stream of CGT. The eastern part of CGT flows slower than the western part because it is grounded along the flow line and thus experiences severe basal drag. Flux mass balance (FMB), i.e., ice thickness change by volume flux divergence, of CGT was obtained by combining the v map and an ice thickness value of 340±18m estimated by the ICESat GLAS data. We used mass conservation assumption in which total mass balance (TMB, −6.29±1.37ma−1 observed by ICESat GLAS data) is attributed to FMB, basal mass balance (BMB) and surface mass balance (SMB, 0.24±0.02ma−1). Mass loss in the freely floating zone of CGT is mainly caused by basal melting (BMB=−144.5±39.9Mta−1) while thinning by volume flux (FMB=−67.2±8.9Mta−1) is relatively small. In the hinge zone of CGT, mass change is contributed to FMB of −147.3±25.2Mta−1 and BMB of 15.3±39.8Mta−1. However, basal freezing derived for the hinge zone may be erroneous as a result of the constant ice thickness assumption extrapolated from the freely floating zone.

An evaluation of mass-balance methods applied to Castle creek Glacier, British Columbia, Canada

AbstractWe estimate the glacier mass balance of a 9.5 km2 mountain glacier using three approaches for balance years 2009, 2010 and 2011. The photogrammetric, GPS and glaciological methods yielded sampling densities of 100, 5 and 2 points km-2, with measurement precisions of ± 0.40, ± 0.10 and ± 0.10 m w.e. respectively. Our glaciological measurements likely include a positive bias, due to omission of internal and basal mass balance, and uncertainty in determining the interface between snow and firn with a probe (±0.10 m w.e.). Measurements from our photogrammetric method include a negative bias introduced by the manual operator and our temperature index model used to correct for different dates of imaging (0.15 m w.e.), whereas GPS measurements avoid these biases. The photogrammetric and GPS methods are suitable for estimating glacier-wide annual mass balance, and thus provide a valuable measure that complements the glaciological method. These approaches, however, cannot be used to estimate mass balance at a point or mass-balance profiles without a detailed understanding of the vertical component of ice velocity.

Calving fluxes and basal melt rates of Antarctic ice shelves

Iceberg calving has been assumed to be the dominant cause of mass loss for the Antarctic ice sheet, with previous estimates of the calving flux exceeding 2,000 gigatonnes per year. More recently, the importance of melting by the ocean has been demonstrated close to the grounding line and near the calving front. So far, however, no study has reliably quantified the calving flux and the basal mass balance (the balance between accretion and ablation at the ice-shelf base) for the whole of Antarctica. The distribution of fresh water in the Southern Ocean and its partitioning between the liquid and solid phases is therefore poorly constrained. Here we estimate the mass balance components for all ice shelves in Antarctica, using satellite measurements of calving flux and grounding-line flux, modelled ice-shelf snow accumulation rates and a regional scaling that accounts for unsurveyed areas. We obtain a total calving flux of 1,321 ± 144 gigatonnes per year and a total basal mass balance of -1,454 ± 174 gigatonnes per year. This means that about half of the ice-sheet surface mass gain is lost through oceanic erosion before reaching the ice front, and the calving flux is about 34 per cent less than previous estimates derived from iceberg tracking. In addition, the fraction of mass loss due to basal processes varies from about 10 to 90 per cent between ice shelves. We find a significant positive correlation between basal mass loss and surface elevation change for ice shelves experiencing surface lowering and enhanced discharge. We suggest that basal mass loss is a valuable metric for predicting future ice-shelf vulnerability to oceanic forcing.

Modeling the basal melting and marine ice accretion of the Amery Ice Shelf

The basal mass balance of the Amery Ice Shelf (AIS) in East Antarctica is investigated using a numerical ocean model. The main improvements of this model over previous studies are the inclusion of frazil formation and dynamics, tides and the use of the latest estimate of the sub‐ice‐shelf cavity geometry. The model produces a net basal melt rate of 45.6 Gt year−1 (0.74 m ice year−1) which is in good agreement with reviewed observations. The melting at the base of the ice shelf is primarily due to interaction with High Salinity Shelf Water created from the surface sea‐ice formation in winter. The temperature difference between the coldest waters created in the open ocean and thein situ freezing point of ocean water in contact with the deepest part of the AIS drives a melt rate that can exceed 30 m of ice year−1. The inclusion of frazil dynamics is shown to be important for both melting and marine ice accretion (refreezing). Frazil initially forms in the supercooled water layer adjacent to the base of the ice shelf. The net accretion of marine ice is 5.3 Gt year−1, comprised of 3.7 Gt year−1 of frazil accretion and 1.6 Gt year−1 of direct basal refreezing.

Basal melting at the Ekström Ice Shelf, Antarctica, estimated from mass flux divergence

AbstractWe characterize the basal mass balance of the Ekström Ice Shelf, Dronning Maud Land, Antarctica, using interferometrically derived surface velocities and ice thickness measurements from radio-echo sounding (RES). The surface velocities are based on data from European Remote-sensing Satellites-1 and -2 (ERS-1/2) during 1994–97. The ice thickness grid consists of 136 RES profiles acquired between 1996 and 2006. Mass fluxes are calculated along selected RES profiles where possible, to reduce uncertainties from ice thickness interpolation. Elsewhere large-scale mass fluxes are calculated using interpolated ice thickness data. Themass flux into the Ekström Ice Shelf from the main grounded drainage basins is estimated to be 3.19±0.4Gt a–1. The mass flux near the ice shelf front is 2.67±0.3Gt a–1. Assuming steady state, and based on the equation of continuity, we interpret the residual mass flux as a combined effect of snow accumulation and subglacial melting/refreezing. Using net snow accumulation rates from previous studies, we link the mass flux divergence in irregular-shaped polygons to processes beneath the ice shelf. The highest subglacial melt rates of ~1.1ma–1 are found near the grounding zone of two main inflow glaciers, and around the German station Neumayer III. The detection of unlikely refreezing in a small area ~15 km west of Neumayer III is attributed to both dataset inaccuracies and a (possibly past) violation of the steady-state assumption. In general, the method and input data allow mapping of the spatial distribution of basal melting and the results are in good agreement with several previous studies.

Ice flux divergence anomalies on 79north Glacier, Greenland

The ice flux divergence of a glacier is an important quantity to examine because it determines the rate of temporal change of its thickness. Here, we combine high‐resolution ice surface velocity observations of Nioghalvfjerdsfjorden (79north) Glacier, a major outlet glacier in north Greenland, with a dense grid of ice thickness data collected with an airborne radar sounder in 1998, to examine its ice flux divergence. We detect large variations, up to 100 m/yr, in flux divergence on grounded ice that are incompatible with what we know of the glacier surface mass balance, basal mass balance and thinning rate. We examine the hypothesis that these anomalies are due to the three‐dimensional flow of ice around and atop bumps and hollows in basal topography by comparing the flux divergence of three‐dimensional numerical models with its surface equivalent. We find that three‐dimensional effects have only a small contribution to the observed anomalies. On the other hand, if we degrade the spatial resolution of the data to 10 km the anomalies disappear. Further analysis shows that the source of the anomalies is not the ice velocity data but the interpolation of multiple tracks of ice thickness data onto a regular grid using a scheme (here block kriging) that does not conserve mass or ice flux. This problem is not unique to 79north Glacier but is common to all conventional ice thickness surveys of glaciers and ice sheets; and fundamentally limits the application of ice thickness grids to high‐resolution numerical modeling of glacier flow.