This study investigates the impact of sea-ice conditions on the production of iceberg plumes from two tidewater glaciers on Devon Ice Cap, Nunavut, Canada. These effects are quantified using a 12 year RADARSAT-1 satellite record from 1997–2008 that contains imagery from approximately every 1–2 weeks in the winter and every 1–4 days in the summer. Iceberg plumes identified in this record are verified against terrestrial time-lapse photography of Belcher Glacier from 2007–08. Results suggest a strong relationship between iceberg plumes and the retreat of sea ice from the glacier termini, with the plumes caused by both the release of previously calved icebergs (ice melange) and new glacier calving. Iceberg plumes are also sometimes observed at other times in the summer and in midwinter (occasionally on both glaciers simultaneously), with these events likely due to new glacier calving alone. Analysis of tides and air temperatures suggests that they provide a minor influence on the timing of iceberg plumes. Instead, it appears that changes in the presence of sea ice are dominant on seasonal timescales, although internal glacier dynamics likely play a significant role for winter plume events that occur when substantial thicknesses of landfast sea ice are present.

We identify a series of basal crevasses along a 31 km transect across the northern sector of the Larsen C ice shelf, Antarctica, using in situ ground-penetrating radar. The basal crevasses propagate from a region of multiple, shallow basal fractures to form widely spaced (0.5–2.0 km) but deeply incised (70–134 m) features. Surface troughs, observed in visible imagery, exist above the basal crevasses as the ice vertically shears to reach hydrostatic equilibrium, while widespread surface crevassing occurs along the crests and on the flanks of the undulations, primarily aligned with the topography. We suggest, based on the location of the surface crevasses and the along-flow evolution of the basal crevasses, that the former are induced by a bending stress created by gradients in hydrostatic forces. Using a linear elastic fracture mechanics model, we investigate the sensitivity of basal crevasse propagation to observed trends of ice-shelf thinning and acceleration. Basal crevasses are large-scale structural weaknesses that can both control meltwater ponding and induce surface crevassing. Together, these features may represent an important mechanism in both past and future ice-shelf disintegration events on the Antarctic Peninsula.

Outlet glaciers draining the Antarctic ice sheet into the Amundsen Sea Embayment (ASE) have accelerated in recent decades, most likely as a result of increased melting of their ice-shelf termini by warm Circumpolar Deep Water (CDW). An ocean model forced with climate reanalysis data shows that, beginning in the early 1990s, an increase in westerly wind stress near the continental shelf edge drove an increase in CDW inflow onto the shelf. The change in local wind stress occurred predominantly in fall and early winter, associated with anomalous high sea-level pressure (SLP) to the north of the ASE and an increase in sea surface temperature (SST) in the central tropical Pacific. The SLP change is associated with geopotential height anomalies in the middle and upper troposphere, characteristic of a stationary Rossby wave response to tropical SST forcing, rather than with changes in the zonally symmetric circulation. Tropical Pacific warming similar to that of the 1990s occurred in the 1940s, and thus is a candidate for initiating the current period of ASE glacier retreat.

Radar power returned from the basal interface along a 42 km long profile over an ice-rise promontory and the adjacent Roi Baudouin ice shelf, Dronning Maud Land, East Antarctica, is analyzed to infer spatial variations in basal reflectivity and hence the basal environment. Extracting basal reflectivity from basal returned power requires an englacial attenuation model. We estimate attenuation in two ways: (1) using a temperature-dependent model with input from thermomechanical ice-flow models; and (2) using a radar method that linearly approximates the geometrically corrected returned power with ice thickness. The two methods give different results. We argue that attenuation calculated using a modeled temperature profile is more robust than the widely used radar method, especially in locations where depth-averaged attenuation varies spatially or where the patterns of basal reflectivity correlate with the patterns of the ice thickness.

Summer 2000/01 ASTER and Landsat 7 scenes and semi-automated digitization were used to compile a glacier inventory for local glaciers of the Geikie Plateau region, central East Greenland. Of the 332 glaciers (41 591 km2), 120 are tidewater-terminating and drain 90% of the glacierized area. Differencing the 2000/01 tidewater margins from the 1980s GEUS map database ice polygons (113 glaciers) and from 2004/05 ASTER tidewater margins (78 glaciers) shows a cumulative tidewater terminus width decrease from 196 km to 183 km between the 1980s and 2000s, with a corresponding areal loss of ~31 km2 and an effective length change of –14.3 km. Between 2000/01 and 2004/05, areal loss was 26 km2. Average margin retreat rate increased two- to threefold, from 1.7–2.1km2 a–1 (1980s– 2000) to 3.9–5.7km2 a–1 (2000–05). Advances are negligible, apart from two surges, of which one was previously undetected. Coastal, ‘outer’ fjord-terminating, glaciers have two to three times larger areal and effective length retreat rates than ‘inner’ fjord-terminating glaciers. We investigate possible causes and hypothesize that, in addition to ocean temperature and sea ice, changes in sea fog may affect coastal-terminating more than inner fjord-terminating glaciers.

We compute the eigenmodes (seiches) of the barotropic and baroclinic hydrodynamic equations for an idealized fjord having length and depth scales similar to those of Ilulissat Icefjord, Greenland, into which Jakobshavn Isbræ (also known as Sermeq Kujalleq) discharges. The purpose of the computation is to determine the fjord’s seiche behavior when forced by iceberg calving, capsize and melange movement. Poorly constrained bathymetry and stratification details are an acknowledged obstacle. We are, nevertheless, able to make general statements about the spectra of external and internal seiches using numerical simulations of ideal one-dimensional channel geometry. Of particular significance in our computation is the role of weakly coupled ice melange, which we idealize as a simple array of 20 icebergs of uniform dimensions equally spaced within the fjord. We find that the presence of these icebergs acts to (1) slow down the propagation of both external and internal seiches and (2) introduce band gaps where energy propagation (group velocity) vanishes. If energy is introduced into the fjord within the period range covered by a band gap, it will remain trapped as an evanescent oscillatory mode near its source, thus contributing to localized energy dissipation and ice/melange fragmentation.

Submarine melting at the ice–ocean interface is a significant term in the mass balance of marine-terminating outlet glaciers. However, obtaining direct measurements of the submarine melt rate, or the ocean heat transport towards the glacier that drives this melting, has been difficult due to the scarcity of observations, as well as the complexity of oceanic flows. Here we present a method that uses synoptic velocity and temperature profiles, but accounts for the dominant mode of velocity variability, to obtain representative heat transport estimates. We apply this method to the Sermilik Fjord–Helheim Glacier system in southeastern Greenland. Using lowered acoustic Doppler current profiler (LADCP) and hydrographic data collected in summer 2009, we find a mean heat transport towards the glacier of 29 × 109W, implying a submarine melt rate at the glacier face of 650 ma–1. The resulting adjusted velocity profile is indicative of a multilayer residual circulation, where the meltwater mixture flows out of the fjord at the surface and at the stratification maximum.

We investigate the sensitivity of sub-ice-shelf melt rates under Pine Island Ice Shelf, West Antarctica, to changes in the oceanic state using an adjoint ocean model that is capable of representing the flow in sub-ice-shelf cavities. The adjoint code is based on algorithmic differentiation (AD) of the Massachusetts Institute of Technology’s ocean general circulation model (MITgcm). The adjoint model was extended by adding into the AD process the corresponding sub-ice-shelf cavity code, which implements a three-equation thermodynamic melt-rate parameterization to infer heat and freshwater fluxes at the ice-shelf/ocean boundary. The inferred sensitivities reveal dominant timescales of 30–60 days over which the shelf exit is connected to the deep interior via advective processes. They exhibit rich three-dimensional time-evolving patterns that can be understood in terms of a combination of the buoyancy forcing by inflowing water masses, the cavity geometry and the effect of rotation and topography in steering the flow in the presence of prominent features in the bedrock bathymetry. Dominant sensitivity pathways are found over a sill, as well as ‘shadow regions’ of very low sensitivities. To the extent that these transient patterns are robust they carry important information for decision-making in observation deployment and monitoring.

The mass balance of Antarctica is one of the crucial factors for determining sea-level change in a warming climate. The marginal zones of the continent, namely the ice shelves, are most sensitive to climate change. During the 2009/10 austral summer an extensive glaciological field campaign was carried out on Fimbulisen, an ice shelf in East Antarctica, to investigate its recent surface mass balance. Shallow (10–18 m) firn cores were drilled and accumulation rates and stable-isotope ratios determined. For firn-core dating, two different methods were compared: (1) seasonal variations of stable oxygen isotope ratios (δ18O), and (2) dielectric profiling, including using the volcanic eruptions of Pinatubo, Philippines (1991), and El Chichόn, Mexico (1982), as time markers. The mean annual accumulation for the period 1992–2009 ranges from 298 to 349 mmw.e. a–1. The interannual variability at the drilling sites is >30%. Accumulation rates show a weak decreasing trend during the past 20–30 years, which is statistically significant only for one of the cores. Stable-isotope ratios were compared to the snowfall temperature of Neumayer station. Neither the temperatures nor the δ18O values show any trend for the investigated time period.

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.