We contrast the frontal dynamics of dilute powder snow avalanches with the behavior of their tail. While the former can be regarded as a fast-moving eruption current dominated by synergistic material injection into a short head, the latter behaves as a nearly arrested dilute cloud of particles expanding by progressive incorporation of ambient air, or by entrainment of snow-cover material by late avalanches trailing well behind the front.

We use time-lapse photography, MODIS satellite imagery, ocean wave measurements and regional broadband seismic data to demonstrate that icebergs that calve from Jakobshavn Isbræ, Greenland, can generate ocean waves that are detectable over 150 km from their source. The waves, which are recorded seismically, have distinct spectral peaks, are not dispersive and persist for several hours. On the basis of these observations, we suggest that calving events at Jakobshavn Isbræ can stimulate seiches, or basin eigenmodes, in both Ilulissat Icefjord and Disko Bay. Our observations furthermore indicate that coastal, land-based seismometers located near calving termini (e.g. as part of the new Greenland Ice Sheet Monitoring Network (GLISN)) can aid investigations into the largely unexplored, oceanographic consequences of iceberg calving.

We assess how the propagation of high-frequency elastic-flexural waves through an ice shelf is modified by the presence of spatially periodic crevasses. Analysis of the normal modes supported by the ice shelf with and without crevasses reveals that a periodic crevasse distribution qualitatively changes the mechanical response. The normal modes of an ice shelf free of crevasses are evenly distributed as a function of frequency. In contrast, the normal modes of a crevasse-ridden ice shelf are distributed unevenly. There are ‘band gaps’, frequency ranges over which no eigenmodes exist. A model ice shelf that is 50 km in lateral extent and 300 m thick with crevasses spaced 500 m apart has a band gap from 0.2 to 0.38 Hz. This is a frequency range relevant for ocean-wave/ice-shelf interactions. When the outermost edge of the crevassed ice shelf is oscillated at a frequency within the band gap, the ice shelf responds very differently from a crevasse-free ice shelf. The flexural motion of the crevassed ice shelf is confined to a small region near the outermost edge of the ice shelf and effectively ‘blocked’ from reaching the interior.

A new implementation of a calving model, using the finite-element code Elmer, is presented and used to investigate the effects of surface water within crevasses on calving rate. For this work, we use a two-dimensional flowline model of Columbia Glacier, Alaska. Using the glacier’s 1993 geometry as a starting point, we apply a crevasse-depth calving criterion, which predicts calving at the location where surface crevasses cross the waterline. Crevasse depth is calculated using the Nye formulation. We find that calving rate in such a regime is highly dependent on the depth of water in surface crevasses, with a change of just a few metres in water depth causing the glacier to change from advancing at a rate of 3.5 kma–1 to retreating at a rate of 1.9 km a–1. These results highlight the potential for atmospheric warming and surface meltwater to trigger glacier retreat, but also the difficulty of modelling calving rates, as crevasse water depth is difficult to determine either by measurement in situ or surface mass-balance modelling.

Simulations of grounding-line migration in ice-sheet models using a fixed grid have been shown to exhibit poor convergence at achievable resolutions. We present a series of ‘shelfy-stream’ flowline model experiments using an idealized set-up. We assess the performance of a range of grounding-line parameterizations (GLPs) over a large input space by varying bedrock gradient, rate factor, basal drag coefficient and net accumulation. The relative performance of GLPs is similar to Gladstone and others (2010a) except at low basal drag, in which case the grounding-line errors are very small for all GLPs. We find that grounding-line errors are far more sensitive to basal drag than to the other inputs or to choice of GLP. We then quantify grounding-line errors as a function of resolution while varying basal drag and channel width (using a parameterization to represent buttressing). Reducing either basal drag or channel width reduces the errors associated with the grounding line. Our results suggest that a structured fixed-grid shelfy-stream ice-sheet model would need to run at a horizontal resolution of ~1–2km to accurately simulate grounding-line positions of marine ice-sheet outlet glaciers such as Pine Island Glacier, Antarctica.

Full-glacier-thickness icebergs are frequently observed to capsize as they calve into the ocean. As they capsize they may collide with the glaciers’ termini; previous studies have hypothesized that such collisions are the source of teleseismic ‘glacial earthquakes’. We use laboratory-scale experiments, force-balance modeling and theoretical arguments to show that (1) the contact forces during these collisions are strongly influenced by hydrodynamic forces and (2) the associated glacial earthquake magnitudes (expressed as twice-integrated force histories) are related to the energy released by the capsizing icebergs plus a hydrodynamic term that is composed of drag forces and hydrodynamic pressure. Our experiments and first-order modeling efforts suggest that, due to hydrodynamic forces, both contact force and glacial earthquake magnitudes may not be directly proportional to the energy released by the capsizing icebergs (as might be expected). Most importantly, however, our results highlight the need to better understand the hydrodynamics of iceberg capsize prior to being able to accurately interpret seismic signals generated by iceberg collisions.

Direct measurements of processes occurring at the ice–ocean boundary are difficult to acquire because of the dangerous and dynamic nature of the boundary, yet these processes are among the least well understood in glaciology. Because sound travels well through water, passive underwater acoustics offers a method to remotely sense activity at this boundary. Here we present passive acoustic measurements and spectral analysis of the evolution of a subaerial calving event and the subsequent mini-tsunami and seiche at Meares Glacier, Alaska, USA. Using two hydrophones to record sound from 1 to 40 000 Hz, we find that each phase of a calving event has distinctive spectral characteristics. An event begins with an infrasound rumble (1–20 Hz), then the ice fractures (20–100 Hz), falls and impacts the water (200–600 Hz). High-frequency (>10 000 Hz) sound increases in intensity quickly as the iceberg oscillates, creating turbulence, spray and waves. Within 10 s, the low-frequency audible sound dissipates and the mini-tsunami and seiche sounds dominate (infrasound plus high frequencies) and continue for over 10 min. The specific frequencies and duration of each phase of a calving event depend on its size and location and the glacier and fjord characteristics.

Several hundred visible and thermal infrared satellite images of Antarctica’s southeast Amundsen Sea from 1986 to 2011, combined with aerial observations in 2009, show a strong inverse relation between prominent curvilinear surface depressions and the underlying basal morphology of the outer Pine Island Glacier ice shelf. Shipboard measurements near the calving front reveal positive temperature, salinity and current anomalies indicative of melt-laden, deep-water outflows near and above the larger channel termini. These buoyant plumes rise to the surface and are expressed as small polynyas in the sea ice and thermal signatures in the open water. The warm upwellings also trace the cyclonic surface circulation in Pine Island Bay. The satellite coverage suggests changing modes of ocean/ice interactions, dominated by leads along the ice shelf through 1999, fast ice and polynyas from 2000 to 2007, and larger areas of open water since 2008.

Ice flow from the ice sheets to the ocean contains the maximum potential contributing to future eustatic sea-level rise. In Antarctica most mass fluxes occur via the extended ice-shelf regions covering more than half the Antarctic coastline. The most extended ice shelves are the Filchner–Ronne and Ross Ice Shelves, which contribute ~30% to the total mass loss caused by basal melting. Basal melt rates here show small to moderate average amplitudes of <0.5ma–1. By comparison, the smaller but most vulnerable ice shelves in the Amundsen and Bellinghausen Seas show much higher melt rates (up to 30 ma–1), but overall basal mass loss is comparably small due to the small size of the ice shelves. The pivotal question for both characteristic ice-shelf regions, however, is the impact of ocean melting, and, coevally, change in ice-shelf thickness, on the flow dynamics of the hinterland ice masses. In theory, ice-shelf back-pressure acts to stabilize the ice sheet, and thus the ice volume stored above sea level. We use the three-dimensional (3-D) thermomechanical ice-flow model RIMBAY to investigate the ice flow in a regularly shaped model domain, including ice-sheet, ice-shelf and open-ocean regions. By using melting scenarios for perturbation studies, we find a hysteresis-like behaviour. The experiments show that the system regains its initial state when perturbations are switched off. Average basal melt rates of up to 2 ma–1 as well as spatially variable melting calculated by our 3-D ocean model ROMBAX act as basal boundary conditions in time-dependent model studies. Changes in ice volume and grounding-line position are monitored after 1000 years of modelling and reveal mass losses of up to 40 Gt a–1.

We have drilled 13 boreholes within and around a through-cutting rift on the (unofficially named) Roi Baudouin Ice Shelf, East Antarctica. Logging by optical televiewer (OPTV) combined with core inspection has resulted in the identification and characterization of several material facies. Outside the rift, OPTV-imaged annual layering indicates ~150 years of accumulation over the 66 m length of one of the boreholes. Luminosity analysis of this image also reveals the presence of numerous dark melt layers as well as a systematic decrease in background luminosity, interpreted in terms of a progressive increase in light transmission during firnification. We identify four material facies within the rift: snow, granular ice, marine ice and unconsolidated platelets. We interpret the granular ice facies as snow that has been saturated by percolating sea water, and the underlying marine ice as compacted buoyant platelets that have adhered to the rift base. Core sections reveal the presence of tubular channels within the marine ice, indicating that it is macroporous and permeable to sea water. The lower boundary of this facies merges into a mushy layer of unconsolidated platelets that were successfully imaged by OPTV, revealing irregular sub-horizontal layering similar to that reported previously on the basis of (directional) borehole video.

Tidewater calving glaciers can undergo large fluctuations not necessarily in direct response to climate, but rather owing to complex ice–water interactions at the glacier termini. One example of this process in Chilean Patagonia is Glaciar Jorge Montt, where two cameras were installed in February 2010, collecting up to four glacier photographs per day, until they were recovered on 22 January 2011. Ice velocities were derived from feature tracking of the geo-referenced photos, yielding a mean value of 13 ±4 md–1 for the whole lower part of the glacier. These velocities were compared to satellite-imagery-derived feature tracking obtained in February 2010, resulting in similar values. During the operational period of the cameras, the glacier continued to retreat (1 km), experiencing one of the highest calving fluxes ever recorded in Patagonia (2.4 km3 a–1). Comparison with previous data also revealed ice acceleration in recent years. These very high velocities are clearly a response to enhanced glacier calving activity into a deep water fjord.

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.