Two high-resolution (1 km grid) numerical model simulations of the Amundsen Sea, West Antarctica, are used to study the role of the ocean in the mass loss and grounding line retreat of Pine Island Glacier. The first simulation uses BEDMAP bathymetry under the Pine Island ice shelf, and the second simulation uses NASA IceBridge-derived bathymetry. The IceBridge data reveal the existence of a trough from the ice-shelf edge to the grounding line, enabling warm Circumpolar Deep Water to penetrate to the grounding line, leading to higher melt rates than previously estimated. The mean melt rate for the simulation with NASA IceBridge data is 28 ma–1, much higher than previous model estimates but closer to estimates from remote sensing. Although the mean melt rate is 25% higher than in the simulation with BEDMAP bathymetry, the temporal evolution remains unchanged between the two simulations. This indicates that temporal variability of melting is mostly driven by processes outside the cavity. Spatial melt rate patterns of BEDMAP and IceBridge simulations differ significantly, with the latter in closer agreement with satellite-derived melt rate estimates of ~50ma–1 near the grounding line. Our simulations confirm that knowledge of the cavity shape and its time evolution are essential to accurately capture basal mass loss of Antarctic ice shelves.

Comparison of the Ross Ice Shelf (RIS, Antarctica) response at near-front seismic station RIS2 with seismometer data collected on tabular iceberg B15A and with land-based seismic stations at Scott Base on Ross Island (SBA) and near Lake Vanda in the Dry Valleys (VNDA) allows identification of RIS-specific signals resulting from gravity-wave forcing that includes meteorologically driven wind waves and swell, infragravity (IG) waves and tsunami waves. The vibration response of the RIS varies with season and with the frequency and amplitude of the gravity-wave forcing. The response of the RIS to IG wave and swell impacts is much greater than that observed at SBA and VNDA. A spectral peak at near-ice-front seismic station RIS2 centered near 0.5 Hz, which persists during April when swell is damped by sea ice, may be a dominant resonance or eigenfrequency of the RIS. High-amplitude swell events excite relatively broadband signals that are likely fracture events (icequakes). Changes in coherence between the vertical and horizontal sensors in the 8–12 Hz band from February to April, combined with the appearance of a spectral peak near 10 Hz in April when sea ice damps swell, suggest that lower (higher) temperatures during austral winter (summer) months affect signal propagation characteristics and hence mechanical properties of the RIS.

Several hundred subglacial lakes have been identified beneath Antarctica so far. Their interaction with the overlying ice sheet and their influence on ice dynamics are still subjects of investigation. While it is known that lakes reduce the ice-sheet friction towards a free-slip basal boundary condition, little is known about how basal melting and freezing at the lake/ice interface modifies the ice dynamics, thermal regime and ice rheology. In this diagnostic study we simulate the Vostok Subglacial Lake area with a coupled full Stokes 3-D ice-flow model and a 3-D lake-circulation model. The exchange of energy (heat) and mass at the lake/ice interface increases (decreases) the temperature in the ice column above the lake by up to 10% in freezing (melting) areas, resulting in a significant modification of the highly nonlinear ice viscosity. We show that basal lubrication at the bottom of the ice sheet has a significant impact not only on the ice flow above the lake itself, but also on the vicinity and far field. While the ice flow crosses Vostok Subglacial Lake, flow divergence is observed and modelled. The heterogeneous basal-mass-balance pattern at the lake/ice interface intensifies this divergence. Instead of interactive coupling between the ice-flow model and the lake-flow model, only a single iteration is required for a realistic representation of the ice/water interaction. In addition, our study indicates that simplified parameterizations of the surface temperature boundary condition might lead to a velocity error of 20% for the area of investigation.

Dynamics of marine-terminating major outlet glaciers are of high interest because of their potential for drawing down large areas of the Greenland ice sheet. We quantify short-term changes in ice flow speed and calving at a major West Greenland glacier and examine their relationship to the presence of the sea-ice melange and tidal stage. A field campaign at the terminus of Store Gletscher (70.40˚N, 50.55˚W) spanning the spring and summer of 2008 included four broadband seismometers, three time-lapse cameras, a tide gauge, an automatic weather station and an on-ice continuous GPS station. Sub-daily fluctuations in speed coincide with two modes of oceanic forcing: (1) the removal of the ice melange from the terminus front and (2) tidal fluctuations contributing to speed increases following ice melange removal. Tidal fluctuations in ice flow speed were observed 16km from the terminus and possibly extend further. Seismic records suggest that periods of intensive calving activity coincide with ice-flow acceleration following breakup of the melange in spring. A synchronous increase in speed at the front and clearing of the melange suggests that the melange directly resists ice flow. We estimate a buttressing stress (~30–60 kPa) due to the presence of the ice melange that is greater than expected from the range of observed tides, though an order of magnitude less than the driving stress.

Ice sheets and smaller ice caps appear to behave in dynamically similar ways; both contain slow-moving ice that is probably frozen to the bed, interspersed with fast-flowing ice streams and outlet glaciers that terminate into the ocean. Academy of Sciences Ice Cap (Akademii Nauk ice cap; 5570 km2), Severnaya Zemlya, Russian High Arctic, provides a clear example of this varied flow regime. We have combined satellite measurements of elevation change and surface velocity to show that variable ice-stream dynamics dominate the mass balance of the ice cap. Since 1988, the ice cap has lost 58±16 Gt of ice, corresponding to ~3% of its mass or 0.16mm of sea-level rise. The climatic mass balance is estimated to be close to zero, and terminus positions have remained stable to within a few kilometers, implying that almost all mass loss has occurred through iceberg calving. The ice-cap calving rate increased from ~0.6 Gt a–1 in 1995 to ~3.0 Gt a–1 in 2000–02, but has recently decreased to ~1.4 Gt a–1 due to a likely slowdown of the largest ice stream. Such highly variable calving rates have not been reported before from High Arctic ice caps, suggesting that these ice masses may be less stable than previously thought.

Interaction of Greenland’s marine-terminating glaciers with the ocean has emerged as a key term in the ice-sheet mass balance and a plausible trigger for their recent acceleration. Our knowledge of the dynamics, however, is limited by scarcity of ocean measurements at the glacier/ocean boundary. Here data collected near six marine-terminating glaciers (79 North, Kangerdlugssuaq, Helheim and Petermann glaciers, Jakobshavn Isbræ, and the combined Sermeq Kujatdleq and Akangnardleq) are compared to investigate the water masses and the circulation at the ice/ocean boundary. Polar Water, of Arctic origin, and Atlantic Water, from the subtropical North Atlantic, are found near all the glaciers. Property analysis indicates melting by Atlantic Water (AW; found at the grounding line depth near all the glaciers) and the influence of subglacial discharge at depth in summer. AW temperatures near the glaciers range from 4.5˚C in the southeast, to 0.16˚C in northwest Greenland, consistent with the distance from the subtropical North Atlantic and cooling across the continental shelf. A review of its offshore variability suggests that AW temperature changes in the fjords will be largest in southern and smallest in northwest Greenland, consistent with the regional distribution of the recent glacier acceleration.

A combination of ground-penetrating radar surveys, physical sedimentology and ice composition measurements has been used to characterize ice and sediment accreted to the southern margin of the McMurdo Ice Shelf, Antarctica. The radar data and surface observations show that the ice-shelf margin consists of strongly layered debris-rich ice that contains marine sediment and fossils. A modified Rayleigh-based distillation system has been used to model the isotopic fractionation from sea water to ice in a closed system. The model of ice formation is consistent with formation during almost complete freezing of a sea-water reservoir. By contrast, ice on the upstream side of the grounding line has formed during the early stages of freezing in which a small fraction of the sea-water reservoir has frozen. The model results and the presence of delicate, well-preserved marine fossils are interpreted as evidence of anchor ice formation on the seabed, and rafting of glaciomarine sediment into the bottom of the ice shelf. We argue that repeated accretion of sea water and marine sediment has produced a stacked sequence of ice and glaciomarine debris that forms shore-parallel ice-cored moraines.

Ice-sheet modelling is an important tool for predicting the possible response of ice sheets to climate change in the past and future. An established ice-sheet model is SICOPOLIS (SImulation COde for POLythermal Ice Sheets), and for this study the previously grounded-ice-only model was complemented by an ice-shelf module. The new version of SICOPOLIS is applied to the Antarctic ice sheet, driven by standard forcings defined by the SeaRISE (Sea-level Response to Ice Sheet Evolution) community effort. A crucial point for simulations into the future is to obtain reasonable initial conditions by a palaeoclimatic spin-up, which we carry out over 125 000 years from the Eemian until today. We then carry out a set of experiments for 500 years into the future, in which the surface temperature and precipitation are kept at their present-day distributions, while sub-ice-shelf melting rates between 0 and 200 ma–1 are applied. These simulations show a significant, but not catastrophic, sensitivity of the ice sheet. Grounded-ice volumes decrease with increasing melting rates, and the spread of the results from the zero to the maximum melting case is ~0.65ms.l.e. (metres sea-level equivalent) after 100 years and ~2.25ms.l.e. after 500 years.

The largest dischargers of ice in Greenland are glaciers that terminate in the ocean and melt in contact with sea water. Studies of ice-sheet/ocean interactions have mostly focused on melting beneath near-horizontal floating ice shelves. For tidewater glaciers, melting instead takes place along the vertical face of the calving front. Here we modify the Massachusetts Institute of Technology general circulation model (MITgcm) to include ice melting from a calving face with the freshwater outflow at the glacier grounding line. We use the model to predict melt rates and their sensitivity to ocean thermal forcing and to subglacial discharge. We find that melt rates increase with approximately the one-third power of the subglacial water flux, and increase linearly with ocean thermal forcing. Our simulations indicate that, consistent with limited field data, melting ceases when subglacial discharge is shut off, and reaches several meters per day when subglacial discharge is high in the summer. These results are a first step toward a more realistic representation of subglacial discharge and of ocean thermal forcing on the subaqueous melting of tidewater glaciers in a numerical ocean model. Our results illustrate that the ice-front melting process is both complex and strongly time-dependent.

We have discovered a band of stones and coarse sand in the Ronne Ice Shelf, Antarctica, some 60 m above the ice shelf’s base, 40 km from its seaward edge and 420 km from the point where the ice originally went afloat. A study of ice-sounding radar data from across the Ronne Ice Shelf has revealed other areas likely to contain debris in significant quantities. It appears that basal debris at the margins of ice streams feeding the ice shelf can be buried in the ice shelf by sea water freezing-on at the ice-shelf base. These findings are evidence for a mechanism active in a present-day ice-sheet/shelf system, which enables icebergs to transport large volumes of ice-rafted debris, and which also provides a potential mechanism for the formation of ice rises near ice fronts. We anticipate that a seismics study of debris melted from the ice shelf and deposited beneath will provide a valuable control on the history of ice-shelf–ocean interactions.

Daugaard Jensen Gletscher, Greenland, is a large tidewater glacier terminating in the northwest corner of the Scoresby Sund fjords. We present a time series of surface flow speeds spanning 1985–2010 based on feature tracking of satellite images. The time series confirms that flow speeds remained stable and reveals a persistent summer acceleration of up to 10% over the lower 10 km of the glacier. The front of the 6 km floating tongue fluctuates by little more than the average size of calved icebergs, ~1 km. While we are unable to detect any imbalance between ice discharge and surface mass balance within our error estimates, observations suggest that the region is losing mass and experiencing decreases in surface elevation. We conclude that as flow speeds and surface mass balance have remained steady since 1985, the shift from balance to imbalance, leading to elevation decrease, must have occurred prior to this date. As for other stable Greenland tidewater glaciers, the seasonal melt cycle is the dominant influence on flow velocity variation but, if the apparent current thinning rates continue, there is potential for the grounding line to retreat, for calving rates to increase and for the glacier to accelerate.

We use the 7 km retreat of Kangerdlugssuaq Glacier (KG), East Greenland, to examine the mechanisms, interactions and relative significance of atmospheric forcing and ice/ocean interactions. Hydrographic data from 1991, 1993 and 2004 show that subtropical waters are common in Kangerdlugssuaq Fjord (KF), and that surface waters were warm in 2004 relative to 1991 and 1993. The main water column was nonetheless warmest in 1991. We contend that while flow of subtropical waters into fjords provides a setting in which rapid glacier retreat can occur, the triggering of retreat depends on additional environmental factors. The climatic variables standing out in our study of KG and KF are air temperature and katabatic winds. Both had strong positive anomalies during winter 2004/05, when KG retreated. We show that proglacial ice melange was absent and that fjord freeze-up did not occur until 11 April 2005, due to warm and windy conditions. We demonstrate that this setting is unusual and hypothesize that exposure to open water in winter months caused the retreat. Calculation of ice-front melt rates shows that discharge of basal meltwater, first from runoff and subsequently from frictional basal heating, should intensify the interaction between glacier and fjord.

We examine the pattern of spreading of warm subtropical-origin waters around Greenland for the years 1992–2009 using a high-resolution (4km horizontal grid) coupled ocean and sea-ice simulation. The simulation, provided by the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) project, qualitatively reproduces the observed warming of subsurface waters in the subpolar gyre associated with changes of the North Atlantic atmospheric state that occurred in the mid-1990s. The modeled subsurface ocean temperature warmed by 1.5˚C in southeast and southwest Greenland during 1994–2005 and subsequently cooled by 0.5˚C; modeled subsurface ocean temperature increased by 2–2.5˚C in central and then northwest Greenland during 1997–2005 and stabilized thereafter, while it increased after 2005 by <0.5˚C in north Greenland. Comparisons with in situ measurements off the continental shelf in the Labrador and Irminger Seas indicate that the model initial conditions were 0.4˚C too warm in the south but the simulated warming is correctly reproduced; while measurements from eastern Baffin Bay reveal that the model initial conditions were 1.0˚C too cold in the northwest but the simulated ocean warming brought modeled temperature closer to observations, i.e. the simulated warming is 1.0˚C too large. At several key locations, the modeled oceanic changes off the shelf and below the seasonal mixed layer were rapidly transmitted to the shelf within troughs towards (model-unresolved) fjords. Unless blocked in the fjords by shallow sills, these warm subsurface waters had potential to propagate down the fjords and melt the glacier fronts. Based on model sensitivity simulations from an independent study (Xu and others, 2012), we show that the oceanic changes have very likely increased the subaqueous melt rates of the glacier fronts, and in turn impacted the rates of glacier flow.

Recent satellite studies have shown that active subglacial lakes exist under the Antarctic ice streams and persist almost to their grounding lines. When the lowest-lying lakes flood, the water crosses the grounding line and enters the sub-ice-shelf cavity. Modeling results suggest that this additional freshwater influx may significantly enhance melting at the ice-shelf base. We examine the spatial and temporal variability in subglacial water supply to the grounding lines of the Siple Coast ice streams, by combining estimates for lake volume change derived from Ice, Cloud and land Elevation Satellite (ICESat) data with a model for subglacial water transport. Our results suggest that subglacial outflow tends to concentrate towards six embayments in the Siple Coast grounding line. Although mean grounding line outflow is ~60m3 s–1 for the entire Siple Coast, maximum local grounding line outflow may temporarily exceed 300 m3 s–1 during the synchronized flooding of multiple lakes in a hydrologic basin. Variability in subglacial outflow due to subglacial lake drainage may account for a substantial portion of the observed variability in freshwater flux out of the Ross Ice Shelf cavity. The temporal variability in grounding line outflow results in a net reduction in long-term average melt rate, but temporary peak melting rates may exceed the long-term average by a factor of three.

Prediction of future changes in dynamics of the Earth’s ice sheets, mass loss and resultant contribution to sea-level rise are the main objectives of ice-sheet modeling. Mass transfer from ice sheet to ocean is, in large part, through outlet glaciers. Subglacial topography plays an important role in ice dynamics; however, trough systems have not been included in bed digital elevation models (DEMS) used in modeling, because their size is close to the model resolution. Using recently collected CReSIS MCoRDs data of subglacial topography and an algorithm that allows topographically and morphologically correct integration of troughs and trough systems at any modeling scale (5 km resolution for SeaRISE), an improved Greenland bed DEM was developed that includes Jakobshavn Isbræ, Helheim, Kangerdlussuaq and Petermann glaciers (JakHelKanPet DEM). Contrasting the different responses of two Greenland ice-sheet models (UMISM and SICOPOLIS) to the more accurately represented bed shows significant differences in modeled surface velocity, basal water production and ice thickness. Consequently, modeled ice volumes for the Greenland ice sheet are significantly smaller using the JakHelKanPet DEM, and volume losses larger. More generally, the study demonstrates the role of spatial modeling of data specifically as input for dynamic ice-sheet models in assessments of future sea-level rise.

We 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.

The Finite Element Sea-ice Ocean Model (FESOM) has been augmented by an ice-shelf component with a three-equation system for diagnostic computation of boundary layer temperature and salinity. Ice-shelf geometry and global ocean bathymetry have been derived from the RTopo-1 dataset. A global domain with a triangular mesh and a hybrid vertical coordinate is used. To evaluate sub-ice-shelf circulation and melt rates for present-day climate, the model is forced with NCEP reanalysis data. Basal mass fluxes are mostly realistic, with maximum melt rates in the deepest parts near the grounding lines and marine ice formation in the northern sectors of the Ross and Filchner–Ronne Ice Shelves, Antarctica. Total basal mass loss for the ten largest ice shelves reflects the importance of the Amundsen Sea ice shelves; the Getz Ice Shelf is shown to be a major meltwater contributor to the Southern Ocean. Despite their modest melt rates, the ‘cold water’ ice shelves in the Weddell Sea are still substantial sinks of continental ice in Antarctica. Discrepancies between the model and observations can partly be attributed to deficiencies in the forcing data or to (sometimes unavoidable) smoothing of ice-shelf and bottom topographies.