We use European Space Agency Envisat data to present the first observational evidence that a Northern Hemisphere tsunami triggered Antarctic ice-shelf calving more than 13 000 km away. The Honshu tsunami of 11 March 2011 traversed the Pacific Ocean in <18 hours where it impinged on the Sulzberger Ice Shelf, resulting in the calving of 125 km2 of ice from a shelf front that had previously been stable for >46 years. This event further illustrates the growing evidence of ocean-wave impact on Antarctic calving and emphasizes the teleconnection between the Antarctic ice sheet and events as far away as the Northern Hemisphere.

The role of basal melting within high-precipitation temperate glaciers represents a significant gap in understanding glacier melting processes. We use a basal melt equation to calculate geothermal and frictional heat-induced basal melt and develop an equation to calculate the rainfall-induced basal melt for Franz Josef Glacier, New Zealand, a high-precipitation, temperate glacier. Additionally, we calculate basal melt due to heat dissipation within water and ice. Data collated from published information on glacier dynamics and climate station readings show that total basal melt contributes on average ∼2.50 m a−1 over the lower to mid-ablation zone (300–1500 m a.s.l.), which is equivalent to >10% of the total ablation for the glacier. This indicates that basal melting is an important component of mass balance for high-precipitation, temperate glaciers.

The microstructures and microchemistry of four US International Trans-Antarctic Scientific Expedition ice cores were examined, at three depths (30, 60, 90 m) each, using scanning electron microscopy, including electron backscattered patterns and energy-dispersive spectroscopy (EDS), in order to assess the relationship between chemical and physical properties. The physical characteristics (grain size, porosity, density, internal surface volume, and crystallographic orientation) at the four sites were inhomogeneous, as expected on the basis of differences in the moisture content of deposited snow and the accumulation rate at their respective locations. Evidence of shallow subgrain boundary formation and trends in internal surface volume, having implications for the study of firn densification and ice-sheet modeling, were also found. Chemical characterization revealed that site-specific variations in particulate concentration and source could accurately be determined using EDS analysis. It was also found that the combination of elements predominant within the sample controls the morphology and microstructural location of the impurities.

A tomography-based methodology for the mass transport characterization of snow is presented. Five samples, characteristic for a wide range of seasonal snow, are considered. Their three-dimensional (3-D) geometrical representations are obtained by micro-computed tomography and used in direct pore-level simulations to numerically solve the governing mass and momentum conservation equations, allowing for the determination of their effective permeability and Dupuit–Forchheimer coefficient. The extension to the Dupuit–coefficient is useful near the snow surface, where Reynolds numbers higher than unity can appear. Simplified semi-empirical models of porous media are also examined. The methodology presented allows for the determination of snow’s effective mass transport properties, which are strongly dependent on the snow microstructure and morphology. These effective properties can, in turn, readily be used in snowpack volume-averaged (continuum) models such as strongly layered samples with macroscopically anisotropic properties.

The Altels hanging glacier in Switzerland broke off on 11 September 1895. The ice volume of this catastrophic rupture was estimated as 4 × 106 m3, the largest icefall event ever observed in the Alps. However, the causes of this collapse are not entirely clear. Based on previous studies, we reanalyzed this break-off event, with the help of a new numerical model, initially developed by Faillettaz and others (2010) for gravity-driven instabilities. The simulations indicate that a break-off event is only possible when the basal friction at the bedrock is reduced in a restricted area, possibly induced by the storage of infiltrated water within the glacier. Further, our simulations reveal a two-step behavior: (1) a first quiescent phase, without visible changes, with a duration depending on the rate of change in basal friction; (2) an active phase with a rapid increase of basal motion over a few days. The general lesson obtained from the comparison between the simulations and available observations is that detectable precursors (crevasse formation and velocity increase) of the destabilization process of a hanging glacier, resulting from a progressive warming of the ice/bed interface towards a temperate regime, will appear only a few days prior to the break-off.

While glacier fabric reflects the accumulated strain, detailed azimuthal information is required to link the microstructure to the flow, and this is not easily gathered at depth. Borehole logging provides a way to obtain a log of azimuthal orientation of tilted stratigraphic features that can be used to orient the core with respect to glacier flow. We demonstrate this using acoustic borehole logs and the ice core from a 162 m borehole in Upper Fremont Glacier, Wind River Range, Wyoming, USA. We measured the dip of tilted dust and bubble layers in the actual ice core, identified them on the borehole log, then used their strike to orient the core sections containing them. We examined the crystal orientation fabric of our samples, using electron backscatter diffraction in a scanning electron microscope. When we compared the orientation of the tilted layers in some samples with the fabric, we found that the normal to the foliation and the c-axes maxima both pointed in the direction of maximum shear stress. This illustrates the usefulness of borehole logs for orienting ice cores after removal from the borehole, and for developing a better understanding of fabric development.

Tidewater glaciers deposit sediment at their terminus, thereby reducing the relative water depth. Reduced water depth can lead to increased glacier stability through decreased rates of iceberg calving, glacier thinning and submarine melting. Here we investigate sedimentation processes at the termini of Kronebreen and Kongsvegen, Svalbard. We mapped the fjord floor bathymetry in August 2009 and calculate sedimentation rates based on our bathymetry and that from a similar study in 2005. A grounding-line fan is developing near the current position of the subglacial stream. An older, abandoned grounding-line fan that likely formed between ∼1987 and 2001 is degrading near the middle of the ice front. Our findings indicate that sediment gravity flows reduce the height of the sediment mound forming at the glacier terminus. The future impact of glacimarine sedimentation processes on glacier stability will depend on the net balance between the observed gravity flows and sediment deposition.

Surface ponds on Antarctic fast ice were examined by measuring temperature, salinity and concentrations of chlorophyll a (Chl-a), dissolved inorganic carbon (DIC) and nutrients (NO3 + NO2, PO4 and SiO2) in the surface pond water and under-ice water. Sea-ice cores were also collected from the bottom of a surface pond (pond-ice core) and from a site away from the pond (bare-ice core). Time-series measurements of surface pond water temperature showed that it varied with solar radiation rather than with air temperature. Comparison of water properties between surface pond water and under-ice water suggested that DIC and nutrients were consumed by biological productivity during pond formation. Depth profiles of nutrient concentrations in the pond-ice core suggested the remineralization of organic matter at the bottom of the surface pond. The Chl-a concentration was lower at the bottom of the pond-ice core than in the bare-ice core, suggesting that surface pond formation reduces ice algae abundance in sea ice because meltwater flushes algae from the porous sea ice into the under-ice water.

Snow avalanches possess two types of kinetic energy: the kinetic energy associated with the mean velocity in the downhill direction and the kinetic energy associated with individual particle velocities that vary from the mean. The mean kinetic energy is directional; the kinetic energy associated with the velocity fluctuations is non-directional in the sense that it is connected to random particle movements. However, the rigid, basal boundary directs the random fluctuation energy into the avalanche. Thus, the random energy flux is converted to free mechanical energy which lifts and dilates the avalanche flow mass, changing the flow density and increasing the normal (dispersive) pressure and, as a consequence, changing the flow resistance. In this paper we derive macroscopic relations that link the production of the random kinetic energy to the perpendicular acceleration of the avalanche’s center of mass. We show that a single burst of fluctuation energy will produce pressures that oscillate around the hydrostatic pressure. Because we do not include a damping process, the oscillations of the center of mass remain, even if the production of random kinetic energy stops. We formulate relationships that can be used within the framework of depth-averaged mass and momentum equations that are often used to simulate snow avalanches in realistic terrain.

Spectral measurements of solar radiation in the band 400–900 nm were performed above and inside the snowpack in two locations in Finland, using a spectroradiometer. The transmittance and extinction coefficient were estimated for different snow layers. Four small candle-shaped photosynthetically active radiation (PAR) sensors were also used to measure the transfer of PAR inside the snowpack. In addition to the light measurements, physical characterization of snow stratigraphy was done, including the thickness, density, hardness (hand test), salinity, and grain size and shape (photographs of crystals). The transmittance varied from <1% (0–12 cm layer) to 80% (0–4 cm layer), and the extinction coefficient was between 0.03 cm−1 (4–8 cm layer) and 0.8 cm −1 (0–4 cm layer). The physical properties of the snow varied considerably between locations and days. The density of the snow varied between 140 and 480 kg m−3.

The geometry of the ice sheets during the Pliocene to early Pleistocene is not well constrained. Here we apply an ice-flow model in the study of the Greenland ice sheet (GIS) during three extreme intervals of this period constrained by geological observations and climate reconstructions. We study the extent of the GIS during the Mid-Pliocene Warmth (3.3–3.0 Ma), its advance across the continental shelf during the late Pliocene to early Pleistocene glaciations (3.0–2.4 Ma) as implied by offshore geological studies, and the transition from glacial to interglacial conditions around 2.4 Ma as deduced from the deposits of the Kap København Formation, North Greenland. Our experiments show that no coherent ice sheet is likely to have existed in Greenland during the Mid-Pliocene Warmth and that only local ice caps may have been present in the coastal mountains of East Greenland. Our results illustrate the variability of the GIS during the Pliocene to early Pleistocene and underline the importance of including independent estimates of the GIS in studies of climate during this period. We conclude that the GIS did not exist throughout the Pliocene to early Pleistocene, and that it melted during interglacials even during the late Pliocene climate deterioration.

Ice cores (∼4 m long) obtained from areas of different surface velocities near the terminus of Sørsdal Glacier, East Antarctica, have been investigated using two versions of a fabric analyser (G50). In sections parallel to the flow plane, the microstructure is typically interlocking with elongate grains that parallel air-bubble elongation, X, reflecting their development in an earlier ductile regime. The c-axis fabric patterns vary with respect to X and vary from single–double maxima to asymmetric small-circle girdles oblique to the planar foliation, which can be attributed to a simple shear regime. The siteto-site variations in the c-axis patterns can be related to areas of different surface velocities, the asymmetry of fabrics correlating with localized strain variations and differences in the deformation path, but not to the current strain pattern recorded by the near-surface deformation conditions. Overprinting fractures have little effect on microstructure except for local dissolution and precipitation along stylolitic surfaces. Comparison of results from the two different fabric analysers reveals that with a higher pixel resolution the incorporation of additional monochromatic light-emitting diodes and repositioning of a retarder plate produce more reliable c-axis measurements.

Recent topographical and ground-penetrating radar (GPR) surveys of Kårsaglaciären, which is a small (<1 km2) and thin (<56 m) mountain glacier in Arctic Sweden, show that there are small areas of temperate ice in the lowermost part of the glacier. This is curious because we would expect such a small and thin glacier to have a fully cold ablation zone. Specifically, with our analyses of present glacier geometry and thickness and of the prevailing climate, we are unable to explain the presence of temperate ice within the snout of Kårsaglaciären using prevailing models of glacier thermal structure. This leads us to suggest that the presence of temperate ice within Kårsaglaciären is a remnant of a previous polythermal state that existed when the glacier was larger and thicker. Kårsaglaciären is thus out of synch with current geometry and climate and is exhibiting a ‘thermal lag’. We propose that, with time, Kårsaglaciären’s ablation zone and perhaps the entire glacier may well become fully cold as the temperate zone shrinks further. We anticipate that such a thermal lag is likely to be present within other Arctic glaciers. A thermal lag and an evolution to a fully cold thermal state have significant implications for the dynamic behaviour of small Arctic glaciers and for meltwater production from them.

Between 2003 and 2007, North Gasherbrum Glacier on the northeastern slope of the Karakoram mountains in Asia underwent a dramatic acceleration, during which a velocity wave propagated down the glacier. There was a significant transfer of ice from up-glacier downstream, which resulted in a strong surface elevation increase over the lower tongue, but only a moderate advance of the glacier snout. We interpret this behaviour as that of a glacier surge, and we explain the observations by means of a simple version of the Kamb drainage-switching theory.

As ice streams flow into the Ross Ice Shelf, West Antarctica, their bed coupling transitions from weak to transient to zero as the ice goes afloat. Here we explore the nature of the bed across these crucial grounding zones using ice-penetrating radar. We collected several ground-based 2 MHz radar transects across the grounding zones of Whillans and Kamb Ice Streams and inferred bed-reflectivity changes from in situ measurements of depth-averaged dielectric attenuation, made possible by the observation of both primary and multiple bed echoes. We find no significant change in the bed reflectivity across either grounding zone. Combined with reflectivity modeling, this observation suggests that a persistent layer of subglacial water (>∼0.2 m) is widespread several kilometers upstream of the grounding zone. Our results are consistent with previous inferences of gradual grounding zones across this sector of the Ross Ice Shelf from airborne radar and satellite altimetry. Separately, the only clear bed-reflectivity change that we observed occurs ∼40 km downstream of the Kamb Ice Stream grounding zone, which we attribute to the onset of marine ice accretion onto the base of the ice shelf. This onset is much nearer to the grounding zone than previously predicted.

The Greenland ice sheet (GrIS) subglacial hydrological system may undergo a seasonal evolution, with significant geophysical and biogeochemical implications. We present results from a new isotope-mixing model to quantify the relative contributions of surface snow, glacial ice and delayed flow to the bulk meltwater discharge from a small (∼5 km2) land-terminating GrIS outlet glacier during melt onset (May) and at peak melt (July). We use radioactive (222Rn) and stable isotopes (18O, deuterium) to differentiate the water source contributions. Atmospherically derived 7Be further constrains meltwater transit time from the glacier surface to the ice margin. We show that (1) 222Rn is a promising tracer for glacial waters stored at the bed and (2) a quantitative chemical mixing model can be constructed by combining 222Rn and the stable water isotopes. Applying this model to the bulk subglacial outflow from our study area, we find a constant delayed-flow (stored) component from melt onset through peak melt. This component is diluted first by snowmelt and then by increasing glacial ice melt as the season progresses. Results from this pilot study are consistent with the hypothesis that subglacial drainage beneath land-terminating sections of the GrIS undergoes a seasonal evolution from a distributed to a channelized system.

Static (or ‘normal’) grain growth, i.e. grain boundary migration driven solely by grain boundary energy, is considered to be an important process in polar ice. Many ice-core studies report a continual increase in average grain size with depth in the upper hundreds of metres of ice sheets, while at deeper levels grain size appears to reach a steady state as a consequence of a balance between grain growth and grain-size reduction by dynamic recrystallization. The growth factor k in the normal grain growth law is important for any process where grain growth plays a role, and it is normally assumed to be a temperature-dependent material property. Here we show, using numerical simulations with the program Elle, that the factor k also incorporates the effect of the microstructure on grain growth. For example, a change in grain-size distribution from normal to log-normal in a thin section is found to correspond to an increase in k by a factor of 3.5.

Sliding experiments were performed at −10°C on smooth surfaces of freshwater columnar-grained S2 ice sliding against itself at a velocity of 8 × 10−4 m s−1, with the purpose of examining whether normal stress parallel to the sliding plane affects frictional resistance. This component of the stress tensor was varied (0.20–1.83 MPa) using a loading system operated under biaxial compression, by orienting the sliding plane at two different angles, 26° and 64°, with respect to the principal loading direction. Under these conditions, no evidence was found to indicate that the normal stress in the direction of sliding affects the friction coefficient.

We provide an assessment of the supraglacial water budget of a moulin basin on the western margin of the Greenland ice sheet for 15 days in August 2009. Meltwater production, the dominant input term to the 1.14 ± 0.06 km2 basin, was determined from in situ ablation measurements. The dominant water-output terms from the basin, accounting for 52% and 48% of output, respectively, were moulin discharge and drainage into crevasses. Moulin discharge exhibits large diurnal variability (0.017–0.54 m3 s−1) with a distinct late-afternoon peak at 16:45 local time. This lags peak meltwater production by ∼2.8 ± 4.2 hours. An Extreme Ice Survey time-lapse photography sequence complements the observations of moulin discharge. We infer, from in situ observations of moulin geometry, previously published borehole water heights and estimates of the temporal lag between meltwater production and observed local ice surface uplift (‘jacking’), that the transfer of surface meltwater to the englacial water table via moulins is nearly instantaneous (<30 min). We employ a simple crevasse mass-balance model to demonstrate that crevasse drainage could significantly dampen the surface meltwater fluctuations reaching the englacial system in comparison to moulin discharge. Thus, unlike crevasses, moulins propagate meltwater pulses to the englacial system that are capable of overwhelming subglacial transmission capacity, resulting in enhanced basal sliding.

We use repeat-track laser altimeter data from the Ice, Cloud and land Elevation Satellite (ICESat) to map the grounding zone (GZ) of the Filchner–Ronne Ice Shelf, Antarctica. Ice flexure in the GZ occurs as the ice shelf responds to ocean-height changes due primarily to tides. We have identified three ‘ice plains’, regions of low surface slope near the GZ where the ice is close to hydrostatic equilibrium: one on Institute Ice Stream; another to its east; and another west of Foundation Ice Stream. The vertical information from repeated ICESat tracks enables us to study the topography, state of flotation and flexure characteristics across these features. In regions of ephemeral grounding, tidal migration of the grounding line allows us to estimate bed slope (∼1–2 × 10−3). From these studies we develop a classification scheme for ice plains, expressed in terms of the evolution, or ‘life cycle’, of these features. A lightly grounded ice plain progresses to a state of ephemeral grounding as the ice sheet thins near the GZ. Once sufficient thinning has occurred, the ice plain becomes a fully floating, relict ice plain with an undulated surface topography similar to that of lightly grounded ice; we expect viscous relaxation to a smooth ice-shelf surface to occur over a timescale of decades. Our improved insight into ice-plain evolution suggests added complexity in modeling ice in the vicinity of the GZ, and a role for ice-plain observations as a guide to relatively rapid changes in ice-sheet mass balance.