Observations from satellite and airborne platforms are combined with model calculations to infer the nature and efficiency of basal melting of the Pine Island Glacier ice shelf, West Antarctica, by ocean waters. Satellite imagery shows surface features that suggest ice-shelf-wide changes to the ocean’s influence on the ice shelf as the grounding line retreated. Longitudinal profiles of ice surface and bottom elevations are analyzed to reveal a spatially dependent pattern of basal melt with an annual melt flux of 40.5 Gt a−1. One profile captures a persistent set of surface waves that correlates with quasi-annual variations of atmospheric forcing of Amundsen Sea circulation patterns, establishing a direct connection between atmospheric variability and sub-ice-shelf melting. Ice surface troughs are hydrostatically compensated by ice-bottom voids up to 150 m deep. Voids form dynamically at the grounding line, triggered by enhanced melting when warmer-than-average water arrives. Subsequent enlargement of the voids is thermally inefficient (4% or less) compared with an overall melting efficiency beneath the ice shelf of 22%. Residual warm water is believed to cause three persistent polynyas at the ice-shelf front seen in Landsat imagery. Landsat thermal imagery confirms the occurrence of warm water at the same locations.

Acoustic sub-bottom profiler surveys on the northeast Antarctic Peninsula shelf indicate that parts of the seabed are underlain by an acoustically transparent layer that is thin on the inner shelf and becomes thicker and more extensive towards the outer shelf. Sedimentological and geophysical data are combined to construct a bed model where streaming ice flow, by both deformation and basal sliding, took place within cross-shelf troughs. The model suggests only limited deformation contributed to fast flow on the inner shelf, i.e. in the onset zone of ice streaming, where the bed was predominantly underlain by a stiff till. Thus, fast ice flow in this area might have been by basal sliding, with deformation confined to discontinuous patches of soft till <40 cm thick. Towards the middle and outer shelf, extensive, thick sequences of soft till suggest a change in the dominant subglacial process towards widespread deformation. This downstream change from basal sliding to subglacial deformation is manifest in the transition from stiff-till dominance to soft-till dominance, while a downstream increase in ice flow velocity is evident from the complex geomorphic imprint on the inner shelf evolving to the more restricted set of bedforms on the outer shelf.

The flow speed of Greenland outlet glaciers is governed by several factors, the relative importance of which is poorly understood. The delivery of surface-generated meltwater to the bed of alpine glaciers has been shown to influence glacier flow speed when the volume of water is sufficient to increase basal fluid pressure and hence basal lubrication. While this effect has also been demonstrated on the Greenland ice-sheet margin, little is known about the influence of surface melting on the large, marine-terminating outlet glaciers that drain the ice sheet. We use a validated model of meltwater input and GPS-derived surface velocities to quantify the sensitivity of glacier flow speed to changes in surface melt at Helheim Glacier during two summer seasons (2007–08). Our observations span ∼55 days near the middle of each melt season. We find that relative changes in glacier speed due to meltwater input are small, with variations of ∼45% in melt producing changes in velocity of ∼2–4%. These velocity variations are, however, of similar absolute magnitude to those observed at smaller glaciers and on the ice-sheet margin. We find that the glacier’s sensitivity to variations in meltwater input decreases approximately exponentially with distance from the calving front. Sensitivity to melt varies with time, but generally increases as the melt season progresses. We interpret the time-varying sensitivity of glacier flow to meltwater input as resulting from changes in subglacial hydraulic routing caused by the changing volume of meltwater input.

Observed meteorological data and a high-resolution (5 km) model were used to simulate Greenland ice sheet surface melt extent and trends before the satellite era (1960–79) and during the satellite era through 2010°. The model output was compared with passive microwave satellite observations of melt extent. For 1960–2010 the average simulated melt extent was 15 ± 5%. For the period 1960–72, simulated melt extent decreased by an average of 6%, whereas 1973–2010 had an average increase of 13%, with record melt extent in 2010. The trend in simulated melt extent since 1972 indicated that the melt extent in 2010 averaged twice that in the early 1970s. The maximum and mean melt extents for 2010 were 52% (∼9.5 × 105 km2) and 28% (∼5.2 × 105 km2), respectively, due to higher-than-average winter and summer temperatures and lower-than-average winter precipitation. For 2010, the southwest Greenland melt duration was 41–60 days longer than the 1960–2010 average, while the northeast Greenland melt duration was up to 20 days shorter. From 1960 to 1972 the melting period (with a >10% melt extent) decreased by an average of 3 days a−1. After 1972, the period increased by an average of 2 days a−1, indicating an extended melting period for the ice sheet of about 70 days: 40 and 30 days in spring and autumn, respectively.

A surface cooling of ∼1.7°C occurred over the ∼two millennia prior to ∼1700 CE at the West Antarctic ice sheet (WAIS) Divide site, based on trends in observed bubble number-density of samples from the WDC06A ice core, and on an independently constructed accumulation-rate history using annual-layer dating corrected for density variations and thinning from ice flow. Density increase and grain growth in polar firn are both controlled by temperature and accumulation rate, and the integrated effects are recorded in the number-density of bubbles as the firn changes to ice. Number-density is conserved in bubbly ice following pore close-off, allowing reconstruction of either paleotemperature or paleo-accumulation rate if the other is known. A quantitative late-Holocene paleoclimate reconstruction is presented for West Antarctica using data obtained from the WAIS Divide WDC06A ice core and a steady-state bubble number-density model. The resultant temperature history agrees closely with independent reconstructions based on stable-isotopic ratios of ice. The ∼1.7°C cooling trend observed is consistent with a decrease in Antarctic summer duration from changing orbital obliquity, although it remains possible that elevation change at the site contributed part of the signal. Accumulation rate and temperature dropped together, broadly consistent with control by saturation vapor pressure.

Abrupt changes in crystal-orientation fabric (COF), and therefore viscosity, are observed near the base of the ice sheet throughout West Antarctica. We report on active-source seismic observations from WAIS Divide, mid-stream and downstream on Thwaites Glacier, and the onset region of Bindschadler Ice Stream. These data reveal a prevalence of englacial seismic reflectivity in the bottom quarter of the ice sheet. The observed seismic reflectivity is complex but largely bed-conformable, with long-spatial-wavelength features observed in the flow direction and short-wavelength features observed across flow. A correspondence of englacial structures with bed features is also observed. We determine the origin of the reflectivity to be abrupt changes in the COF of ice, based on the following: (1) observations of englacial reflectivity are consistent with current knowledge of COF within ice sheets, (2) englacial reflectivity caused by COF contrasts requires the simplest genesis, especially at ice divides, and (3) amplitude analysis shows that the observed englacial reflectivity can be explained by contrasts in seismic velocity due to COF changes. We note that the downstream increase in the quantity and complexity of observations indicates that direct observations of COF at ice divides likely underestimate the role that fabric plays in ice-sheet dynamics.

There is growing recognition of the significance of biologically active supraglacial dust (cryoconite) for glacial mass balance and ecology. Nonetheless, the processes controlling the distribution, transport and fate of cryoconite particles in the glacial system remain somewhat poorly understood. Here, using a 216 hour time series of plot-scale (0.04 m2) images, we quantify the small-scale dynamics of cryoconite on Longyearbreen, Svalbard. We show significant fluctuations in the apparent cryoconite area and dispersion of cryoconite over the plot, within the 9 day period of observations. However, the net movement of cryoconite across the ice surface averaged only 5.3 mm d−1. High-resolution measurements of cryoconite granule motion showed constant, random motion but weak correlation with meteorological forcing factors and no directional trends for individual particle movement. The high-resolution time-series data suggest that there is no significant net transport of dispersed cryoconite material across glacier surfaces. The areal coverage and motion of particles within and between cryoconite holes appears to be a product of differential melting leading to changes in plot-scale microtopography, local meltwater flow dynamics and weather-dependent events. These subtle processes of cryoconite redistribution may be significant for supraglacial albedo and have bearing on the surface energy balance at the glacier scale.

Stakes at 2 km intervals were installed in January 1997 and remeasured in February 1998, January 1999, January 2005 and during the 2007/08 austral summer along a 1248 km traverse route from Zhongshan station to Dome A, East Antarctica. Based on topographical parameters, meteorological features and the records of ∼650 stakes and six stake arrays, the route is divided into five zones. We find that the snow accumulation rate decreases with increasing altitude as one progresses inland, except in the zone 800–1128 km from the coast, where the average annual accumulation rate is higher than in the zone 524–800 km from the coast. The Dome A zone (1128–1248 km) has the lowest accumulation rate (35 kg m−2 a−1, 2005–08) due to having the highest elevation and being furthest from the coast. The surface mass balance in the region 202–1128 km from the coast exhibits no temporal change from 1999–2005 to 2005–08, but there is a change in the accumulation distribution. The zone from 202 to 524 km shows a decrease in surface mass balance from 84 kg m−2 a−1 in 1999–2005 to 67 kg m−2 a−1 in 2005–08, while the zone between 800 and 1128 km shows an increase from 67 kg m−2 a−1 in 1999–2005 to 75 kg m−2 a−1 in 2005–08.

We have measured the glacier area changes in the central Southern Alps, New Zealand, between 1978 and 2002 and have compiled the 2002 glacier outlines using an image scene from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER). Three automated classification methods were tested: (1) band ratio, (2) normalized-difference snow index and (3) supervised classification. The results were compared with the glacier outlines photo-interpreted from the ASTER data, and were further validated using GPS-aided field mapping of selected test glaciers. The ASTER 3/4 band ratio provided the best results. However, all the classification methods failed to extract extensive debris-covered parts of the glaciers. Therefore, the photo-interpreted 2002 outlines were used when comparing with the existing 1978 glacier inventory derived from aerial photographs. Our results show a ∼17% reduction of glacier area, mainly driven by the retreat of the large valley glaciers. Despite the large climatic gradient from west to east, glaciers on both sides of the Main Divide lost similar percentages of area, except Franz Josef and Fox Glaciers which advanced. Smaller glaciers were found to have changed very little in the study period.

Taylor Glacier, an outlet of the East Antarctic ice sheet, flows through the Transantarctic Mountains and terminates in the Dry Valleys. Understanding how this glacier fluctuates is important for studies of glacial geology, paleoclimate, ice dynamics and ecology. Sublimation is the primary mass-loss process for most of the glacier. Four years of specific balance measurements from the ablation zone show sublimation rates up to 40 cm a−1. We used data from an array of weather stations as inputs to a model for latent heat flux and hence sublimation rate. Calculated and measured ablation rates agree to within uncertainties, indicating that wind speed and vapour pressure gradient (a function of temperature and humidity) are the governing variables, as expected from theory. Measurements and model results together allowed us to examine the spatial and temporal variations of sublimation on the glacier. On average, sublimation is about two times faster in summer than winter. Rapid sublimation occurs during storms and katabatic wind events, but such periods contribute less to the annual total than do slow, persistent losses. Spatially, sublimation reaches a maximum midway along the glacier, where descending surface air currents are focused by the topography of the aptly named tributary, Windy Gully.

We apply a novel one-dimensional glacier hydrology model that calculates hydraulic head to the tidewater-terminating Sermeq Avannarleq flowline of the Greenland ice sheet. Within a plausible parameter space, the model achieves a quasi-steady-state annual cycle in which hydraulic head oscillates close to flotation throughout the ablation zone. Flotation is briefly achieved during the summer melt season along a ∼17 km stretch of the ∼50 km of flowline within the ablation zone. Beneath the majority of the flowline, subglacial conduit storage ‘closes’ (i.e. obtains minimum radius) during the winter and ‘opens’ (i.e. obtains maximum radius) during the summer. Along certain stretches of the flowline, the model predicts that subglacial conduit storage remains open throughout the year. A calculated mean glacier water residence time of ∼2.2 years implies that significant amounts of water are stored in the glacier throughout the year. We interpret this residence time as being indicative of the timescale over which the glacier hydrologic system is capable of adjusting to external surface meltwater forcings. Based on in situ ice velocity observations, we suggest that the summer speed-up event generally corresponds to conditions of increasing hydraulic head during inefficient subglacial drainage. Conversely, the slowdown during fall generally corresponds to conditions of decreasing hydraulic head during efficient subglacial drainage.

Dispersed facies basal ice – massive (i.e. structureless) ice with dispersed debris aggregates – is present at the margins of many glaciers and, as a product of internal glacial processes, has the potential to provide important information about the mechanisms of glacier flow and the nature of the subglacial environment. The origin of dispersed facies is poorly understood, with several hypotheses having been advanced for its formation, and there is disagreement as to whether it is largely a sedimentary or a tectonic feature. We test these established hypotheses at the temperate glacier Svínafellsjökull, Iceland, and find that none fully account for dispersed facies characteristics at this location. Instead, dispersed facies physical, sedimentological and stable-isotope (δ18O, δD) characteristics favour a predominantly tectonic origin that we suggest comprises the regelation and strain-induced metamorphism of debris-rich basal ice that has been entrained into an englacial position by tectonic processes operating at the base of an icefall. Further thickening of the resultant dispersed facies may also occur tectonically as a result of ice flow against the reverse bed slope of a terminal overdeepening. Lack of efficient subglacial drainage in the region of the overdeepening may limit basal melting and thus favour basal ice preservation, including the preservation of dispersed facies. Despite the relatively low sediment content of dispersed facies ( ∼1.6% by volume), its thickness (up to 25 m) and ubiquity at Svínafellsjökull results in a significant contribution to annual sediment discharge (1635–3270 m3 a−1) that is ∼6.5 times that contributed by debris-rich stratified facies basal ice.

Using a snow/atmosphere coupled model, the evolution of the surface and near-surface snow temperature is modeled at Dome C, Antarctica, during the period 20–30 January 2010. Firstly, the detailed multilayer snow model Crocus is run in stand-alone mode, with meteorological input forcing data provided by local meteorological observations. The snow model is able to simulate the evolution of surface temperature with good accuracy. It reproduces the observed downward propagation of the diurnal heatwave into the upper 50 cm of the snowpack reasonably well. Secondly, a fully coupled 3-D snow/atmosphere simulation is performed with the AROME regional meteorological model, for which the standard single-layer snow parameterization is replaced by Crocus. In spite of a poor simulation of clouds, the surface and near-surface snow temperatures are correctly simulated, showing neither significant bias nor drifts during the simulation period. The model reproduces particularly well the average decrease of the diurnal amplitude of air temperature from the surface to the top of the 45 m instrumented tower. This study highlights the potential of snow/atmosphere coupled models over the Antarctic plateau and the need to improve cloud microphysics and data assimilation over polar regions.

We investigate the elevation and mass-balance response of tributary glaciers following the loss of the Larsen A and B ice shelves, Antarctic Peninsula (in 1995 and 2002 respectively). Our study uses MODIS imagery to track ice extent, and ASTER and SPOT5 digital elevation models (DEMs) plus ATM and ICESat laser altimetry to track elevation changes, spanning the period 2001–09. The measured Larsen B tributary glaciers (Hektoria, Green, Evans, Punchbowl, Jorum and Crane) lost up to 160 m in elevation during 2001–06, and thinning continued into 2009. Elevation changes were small for the more southerly Flask and Leppard Glaciers, which are still constrained by a Larsen B ice shelf remnant. In the northern embayment, continued thinning of >3 m a−1 on Drygalski Glacier, 14 years after the Larsen A ice shelf disintegrated, suggests that mass losses for the exposed Larsen B tributaries will continue for years into the future. Grounded ice volume losses exceed 13 km3 for Crane Glacier and 30 km3 for the Hektoria–Green–Evans glaciers. The combined mean loss rate for 2001–06 is at least 11.2 Gt a−1. Our values differ significantly from published mass-budget-based estimates for these embayments, but are a reasonable fraction of GRACE-derived rates for the region (∼40 Gt a−1).

We present a microstructural characterization of fine-grained layers from the top 90 m of firn from Summit, Greenland, performed using a combination of scanning electron microscopy techniques including secondary electron imaging, energy-dispersive spectroscopy and electron backscattered patterns, and X-ray microcomputed tomography. The impurities in the firn, both soluble impurities and dust particles, were found largely in the grain interiors. Both c- and a-axis pole figures do not show strong evidence of a preferred orientation of the grains even at the bottom of the firn column. The firn structure became increasingly anisotropic with vertical alignment in the top 3 m, probably due to vapor transport associated with dry-snow metamorphism. The anisotropy decreases below this level until at ∼50 m the average firn structure is close to isotropic. In the near surface, the level of anisotropy is weaker than at Hercules Dome, Antarctica, confirming that differences in accumulation rates and temperatures leave enduring evidence in the structure of the firn. The fraction of closed-off pores is relatively low until ∼65 m; below that it rises through the end of our sampling at 90 m. Our microstructure measurements on the microscale are consistent with in situ firn-air sampling measurements on a decimeter scale, both indicating the existence of the lock-in zone starting near 69 m depth, and pore close-off at 81 m at this site.

A novel ‘selection curve’ method is developed to interpret annual layers in the West Antarctic ice sheet (WAIS) Divide ice core based on dielectric properties (DEP). Because dielectric measurements are non-contact and represent the integrated response of the ice volume, they are particularly useful for the brittle zone of the core. Seasonal differences in ice chemistry create an annual signal in DEP, though multiple peaks of varying strength within a year may complicate the interpretation of annual layers. The selection curve algorithm uses a spline curve whose shape selects successive annual peaks in plots of DEP. This spline curve was scaled to the average annual-layer thickness at a given depth, where the layer thickness was best estimated using the fast Fourier transform (FFT) power spectrum within a sliding 10 m window. To explore the accuracy and stability of the method, several spline curves were generated from varying lengths of calibration data taken from multiple depths in the WAIS core. Using 50 m of manually interpreted calibration data, the selection curve method matched a manual interpretation throughout the entire 1200 m dataset to within 2% root-mean-square error (RMSE). This method is equally applicable to glaciochemical and other time/depth series measurements.

Jakobshavn Isbræ, West Greenland, underwent a large, rapid and not well understood change in flow dynamics in 1998, leading to a doubling of its ice discharge rates. We calculate the width-averaged forces controlling flow of Jakobshavn Isbræ in 1995, 2000 and 2005 to elucidate processes responsible for this change in flow speed. In contrast to earlier suggestions, we conclude that the observed acceleration was not caused by the loss of back-stress due to weakening and subsequent break-up of the floating ice tongue alone. Gradients in longitudinal stress are small at all times considered (∼3% of the driving stress) and basal and lateral drag provide resistance to flow. Over the 10 year period considered, the average driving stress increased by 20 kPa, which was balanced by a comparable increase in lateral drag. We surmise that the velocity changes resulted from weakening of the ice in the lateral shear margins and perhaps a change in properties at the bed. Possible mechanisms for weakening of ice in the lateral shear margins include cryo-hydrologic warming of subsurface ice in the ablation zone and hydraulic weakening due to higher water content of ice in the shear margins.