Information on bed topography and basal conditions is essential to developing the next-generation ice-sheet models needed to generate a more accurate estimate of ice-sheet contribution to sea-level rise. Synthetic aperture radar (SAR) images of the ice–bed can be analyzed to obtain information on bed topography and basal conditions. We developed a wideband SAR, which was used during July 2005 to perform measurements over a series of tracks between the GISP2 and GRIP cores near Summit Camp, Greenland. The wideband SAR included an eight-element receive-antenna array with multiple-phase centers. We applied the MUltiple SIgnal Classification (MUSIC) algorithm, which estimates direction of arrival signals, to single-pass multichannel data collected as part of this experiment to obtain fine-resolution bed topography. This information is used for producing fine-resolution estimates of bed topography over a large swath of 1600 m, with a 25 m posting and a relative accuracy of approximately 10 m. The algorithm-derived estimate of ice thickness is within 10 m of the GRIP ice-core length. Data collected on two parallel tracks separated by 500 m and a perpendicular track are compared and found to have difference standard deviations of 9.1 and 10.3 m for the parallel and perpendicular tracks, respectively.

Firn microstructure is accurately characterized using images obtained from scanning electron microscopy (SEM). Visibly etched grain boundaries within images are used to create a skeleton outline of the microstructure. A pixel-counting utility is applied to the outline to determine grain area. Firn grain sizes calculated using the technique described here are compared to those calculated using the techniques of Gow (1969) and Gay and Weiss (1999) on samples of the same material, and are found to be substantially smaller. The differences in grain size between the techniques are attributed to sampling deficiencies (e.g. the inclusion of pore filler in the grain area) in earlier methods. The new technique offers the advantages of greater accuracy and the ability to determine individual components of the microstructure (grain and pore), which have important applications in ice-core analyses. The new method is validated by calculating activation energies of grain boundary diffusion using predicted values based on the ratio of grain-size measurements between the new and existing techniques. The resulting activation energy falls within the range of values previously reported for firn/ice.

The temporal variability of surface snow and glacier melt flux and runoff are investigated for the ablation area of Jakobshavn Isbræ, West Greenland. High-resolution meteorological observations both on and outside the Greenland ice sheet were used as model input. SnowModel, a physically based spatially distributed meteorological and snow evolution modeling system, is used to simulate the temporal variability of Jakobshavn Isbræ accumulation and ablation processes for 2000/01–2006/07. Winter snow depth observations and MODIS satellite-derived summer melt observations are used for model validation of accumulation and ablation. The modeled interannual runoff variability varied from 1.81 × 109 m3 (2001/02) to 5.21 × 109 m3 (2004/05), yielding a cumulative runoff at the Jakobshavn Glacier terminus of ∼2.25 to ∼4.5 m w.e. The average modeled Jakobshavn runoff of ∼3.4 km3 a−1 was merged with previous estimates of Jakobshavn ice discharge to quantify the freshwater flux to Illulissat Icefjord. For both runoff and ice discharge the average trends are similar, indicating increasing (insignificant) influx of fresh water to Ilulissat Icefjord for the period 2000/01–2006/07. This study suggests that surface runoff forms a minor part of the overall Jakobshavn freshwater flux to the fjord: about 7% (∼3.4 km3 a−1) of the average annual freshwater flux of ∼51.0 km3 a−1 originates from the surface runoff.

The volume- and area-change evolution of glaciers can be obtained by employing the volume–area scaling approach during mass-balance modelling. This method usually requires information on the initial surface area and ice volume to adjust the volume–area relation to the specific ice body. However, absolute volumetric data on glaciers are very rare, so the applicability of volume–area scaling is limited. In order to use volume–area scaling on glaciers for which only limited information is available, a new method is presented to calibrate the volume–area relation without prior knowledge of this relation by using glacier extent information from different times. To validate the method and illustrate its practicability, we model the range of probable future changes in ice volume and surface area of ‘Glaciar Noroeste’, an outlet glacier of Gran Campo Nevado ice cap, southern Chilean Patagonia, during the 21st century, based on IPCC SRES scenarios B1 and A2.

A simple circuit has been developed to allow measurement of brief subglacial water-pressure pulses. This circuit continuously powers a pressure transducer and captures the peak output of the transducer during each measurement interval, thus allowing determination of the maximum pressure attained during the interval. This circuit provides an alternative to setting a data logger to perform rapid repeated measurements, and overcomes some key limitations imposed by rapid measurement. Benefits include significantly lower demands on the data-logger microprocessor, which allows additional instruments to be monitored simultaneously, reduced memory usage and moderately lower power consumption. The reduced microprocessor and memory loads allow older and slower logger models, many of which are still in common use, to be used to obtain data that compare favourably with high-frequency data obtained using newer data loggers.

We present results from ground-penetrating radar (GPR) and seismic surveys for the Mullins Valley debris-covered glacier (Mullins Glacier), Antarctica, that yield local ice-thickness estimates of 80–110 m in upper Mullins Valley and 150 m in upper Beacon Valley. Englacial debris in upper Mullins Glacier occurs as scattered cobbles and as discrete layers. One extensive englacial debris layer, which appears as a coherent reflector dipping 40–45° up-valley, intersects the ground surface within an ∼8 m high ice-cored ridge, the largest of several ridges that mark the glacier surface. Field excavations reveal that this englacial interface consists of multiple debris bands that can be directly correlated with ridge microtopography. Englacial debris layers most probably originate as concentrated rockfall in ice accumulation zones and/or as surface lags that form as dirty ice sublimes during periods of negative mass balance. A similar pattern of surface ridges on Friedman Glacier (∼2.5 km west of Mullins Glacier) suggests regional environmental changes are involved in ridge formation. These observations carry implications for evaluating debris entrainment and surface ridge formation mechanisms in cold-based, debris-covered glaciers and provide a glaciological framework for evaluating and interpreting paleoclimate records from Mullins Glacier.

Digital elevation models (DEMs) of the ablation area of Hailuogou glacier, China, produced from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data obtained in 2009, differential GPS (DGPS) data surveyed in 2008 and aerial photographs acquired in 1966 and 1989 are differenced to estimate long- and short-term glacier surface elevation change (dh/dt). The mean dh/dt of the ablation area over 43 years (1966–2009) is −1.1 ± 0.4 m a−1. Since 1989 the thinning has accelerated significantly. Ice velocities measured by DGPS at 28 fixed stakes implanted in the ablation area increase with distance from the glacier terminus, ranging from 41.0 m a−1 approaching the glacier terminus to a maximum of 205.0 m a−1 at the base of an icefall. Our results reveal that the overall average ice velocity in the ablation area has undergone significant temporal variability over the past several decades. Changes in glacier surface elevation in the ablation area result from the combined effects of climate change and glacier dynamics, which are driven by different factors for different regions and periods.

Deposition of snow from precipitation and wind events creates layering within seasonal snowpacks. The thickness and horizontal continuity of layers within seasonal snowpacks can be highly variable, due to snow blowing around topography and vegetation, and this has important implications for hydrology, remote sensing and avalanche forecasting. In this paper, we present practical field and post-processing protocols for recording lateral variations in snow stratigraphy using near-infrared (NIR) photography. A Fuji S9100 digital camera, modified to be sensitive to NIR wavelengths, was mounted on a rail system that allowed for rapid imaging of a 10 m long snow trench excavated on the north side of Toolik Lake, Alaska (68°38′ N, 149°36′ W). Post-processing of the images included removal of lens distortion and vignetting. A tape measure running along the base of the trench provided known locations (control points) that permitted scaling and georeferencing. Snow layer heights estimated from the NIR images compared well with manual stratigraphic measurements made at 0.2 m intervals along the trench (n = 357, R2 = 0.97). Considerably greater stratigraphic detail was captured by the NIR images than in the manually recorded profiles. NIR imaging of snow trenches using the described protocols is an efficient tool for quantifying continuous microscale variations in snow layers and associated properties.

The basal melting and freezing rates under the Amery Ice Shelf, East Antarctica, are evaluated, and their spatial distributions mapped. Ice velocity, surface elevation and accumulation rate datasets are employed in the analysis, along with a column-averaged ice density model. Our analysis shows that the total area of basal melting is 34 700 km2, with a total annual melt of 62.5 ± 9.3 Gt and an average melting rate of 1.8 ± 0.3 m a−1. Basal freezing mainly occurs in the northwestern part of the ice shelf, over a total area of 26 100 km2 and with a maximum freezing rate of 2.4 ± 0.4 m a−1. The total marine ice that accretes to the ice-shelf base is estimated to be 16.2 ± 2.4 Gt a−1. Using a redefined grounding line and geometry of the Amery Ice Shelf, we estimate the net melt over the ice-shelf base is about 46.4 ± 6.9 G ta−1, which is higher than previous modeling and oceanographic estimates. Net basal melting accounts for about half of the total ice-shelf mass loss, with the rest being from iceberg discharge. Our basal melting and freezing distribution map provides a scientific basis for quantitative analysis of ice–ocean interaction at the ice-shelf–ocean interface.

Supraglacial Lac de Rochemelon was formed 50 years ago behind an ice dam and grew steadily until 2004. In October 2004, the volume of the lake reached 650 000 m3, bringing its surface within 0.2 m of the top of the ice dam. To eliminate the threat to towns located below in the event of an overflow, the lake was drained artificially in October 2004 and during the summer of 2005. Once the volume had been reduced to 250 000 m3 by siphoning, a channel was dug with explosives and the remaining water overflowed naturally. This offered a very good opportunity to investigate the breaching of an ice dam accompanied by thermal erosion of the drainage channel. Extensive field measurements were carried out during drainage. Analysis of the energy dissipated in the channel reveals that only half of the available energy was used for breach erosion. A numerical model was used to simulate the evolution of a number of variables during drainage and to study the sensitivity of discharge and ice erosion to different parameters, revealing a high sensitivity to water temperature. Model simulations indicate that natural drainage of this lake at the beginning of October 2004 would have led to a peak discharge of <6 m3 s−1.

Ice types, albedos and impurity content are characterized for the ablation zone of the Greenland ice sheet in Kronprinz Christians Land (80° N, 24° W). Along this ice margin the width of the ablation zone is only about 8 km. The emergence and melting of old ice in the ablation zone creates a surface layer of dust that was originally deposited with snowfall high on the ice sheet. This debris cover is augmented by locally derived wind-blown sediment. Subsequently, the surface dust particles often aggregate together to form centimetre-scale clumps that melt into the ice, creating cryoconite holes. The debris in the cryoconite holes becomes hidden from sunlight, raising the area-averaged albedo relative to surfaces with uniform debris cover. Spectral and broadband albedos were obtained for snow, ice hummocks, debris-covered ice, cryoconite-studded ice and barren tundra surfaces. Broadband ice albedos varied from 0.2 (for ice with heavy loading of uniform debris) to 0.6 (for ice hummocks with cryoconite holes). The cryoconite material itself has albedo 0.1 when wet. Areal distribution of the major surface types was estimated visually from a transect video as a function of distance from the ice edge (330 m a.s.l.). Ablation rates were measured along a transect from the ice margin to the slush zone 8 km from the margin (550 m a.s.l.), traversing both Pleistocene and Holocene ice. Ablation rates in early August averaged 2 cm d−1. Impurity concentrations were typically 4.3 mg L−1 in the subsurface ice. Surface concentrations were about 16 g m−2 on surfaces with low impurity loading, and heavily loaded surfaces had concentrations as high as 1.4 kg m−2. The mineralogical composition of the cryoconite material is comparable with that of the surrounding soils and with dust on a snowdrift in front of the ice margin, implying that much of the material is derived from local sources. A fine mode (clay) is present in the oldest ice but not in the nearby soil, suggesting that its origin is from wind deposition during Pleistocene glaciation.

The sensitivity of glaciers to climatic change is key information in assessing the response and sea-level implications of projected future warming. New Zealand glaciers are important globally as an example of how maritime glaciers will contribute to sea-level rise. A spatially distributed energy-balance model is applied to Brewster Glacier, New Zealand, in order to calculate glacier mass balance, run-off and sensitivity to climate change. The model successfully simulates four annual mass-balance cycles. Close to half (52%) of the energy available for melt on the glacier is supplied by turbulent heat fluxes, with radiation less important, except during the winter. Model sensitivity to temperature change is one of the largest reported on Earth, at −2.0 m w.e. a−1 °C−1. In contrast, a 50% change in precipitation is required to offset the mass-balance change resulting from a 1 °C temperature change. Meltwater runoff sensitivity is also very high, increasing 60% with a 1°C warming. The extreme sensitivity of mass balance to temperature change suggests that significant ice loss will occur with even moderate climate warming.

The southernmost glaciers in Europe are located on the Iberian, Apennine and Balkan Peninsulas in mid-latitudes between 41° N and 44° N at altitudes ranging from 2000 to 3000 m a.s.l. All these glaciers are a legacy of the Little Ice Age (LIA). They survive in a relatively warm environment (mean annual temperature 0°C to +1°C) due to local topographic controls and high levels of accumulation as a result of avalanche and wind-blown snow. In the Pirin Mountains, Bulgaria, Snezhnika glacieret has been cored, providing an archive of recent climate change. Small glaciers such as this respond quickly to climatic extremes. Since the LIA maximum during the 19th century, all southern European glaciers have retreated, losing 30–100% of their volume. However, despite the trend towards warmer years since the late 1970s, some glaciers still survive, even after some of the hottest summers on record. Predicted future warming, especially in summer, and drier conditions in the Mediterranean basin may result in the disappearance of all glacier features at these latitudes in Europe within the next few decades.

Digital optical televiewing (OPTV) of hot-water-drilled boreholes is evaluated as a technique for the investigation of englacial ice and debris structures on the basis of six boreholes drilled in the terminus region of midre Lovénbreen, Svalbard. The resulting OPTV logs successfully reveal several visually distinctive englacial ice properties and deformation structures (e.g. oblique englacial fractures imaged here for the first time). Combining these OPTV logs with surface mapping has resulted in the identification of eight separate structural elements, several of which can be interpolated onto 3-D grids at a node spacing of 1 m vertically and 10 m horizontally. Basally derived englacial sediment layers are also found to be intercalated with primary stratification, elevated into near-vertical planes around a central fold axis by large-scale lateral folding. The analysis also allows supraglacial longitudinal debris ridges to be subclassified into two types: a previously described (type-I) form, which are the exposed fold axes of large-scale lateral folds, and a new (type-II) form experiencing secondary deformation by small-scale horizontal folding in association with vertical displacements across arcuate shear planes in response to longitudinally compressive stresses near the glacier terminus. Although using boreholes to investigate glacier structure is limited (e.g. by parallelism with vertical planes), applying OPTV to multiple boreholes at midre Lovénbreen has successfully revealed a range of 3-D structural elements at high spatial resolution. As such, interpolating between multiple OPTV logs overcomes many of the problems associated with interpretations made solely on the basis of surface-based structural mapping, and combining the two techniques represents a powerful glaciological tool.

The ice streams feeding the Ross Ice Shelf, Antarctica, have large tidally modulated (sinusoidal and stick–slip) flow, but the interaction with the ice shelf is poorly understood. We show that the flow of the Ross Ice Shelf front, up to ∼650 km from the ice streams, exhibits smooth, sinusoidal motions corresponding to tidal modulation. These observations suggest a possible linking of the ice shelf with the ice streams to form a unified system that responds to small perturbations in stresses associated with ocean tides. If this is the case, the presence of the sinusoidal motion but the absence of stick–slip motion suggests there is damping of very high-frequency signals. The dissimilar signatures of the motions observed in the ice streams and at the front of the ice shelf present challenges to model development aimed at understanding the dynamics of coupled ice-stream/ice-shelf flow and the movement of ice across grounding lines.

We present a newly developed high-resolution frequency-modulated continuous-wave (FMCW) radar system for sounding the sub-ice topography. The system, working in the frequency range from 200 to 400 MHz, was developed to measure thickness and internal layers in cold ice with a resolution better than 1 m. This system has the potential to measure accumulation rates, an important input for improved knowledge of the mass balance of polar ice sheets. First measurements for the test, calibration and optimization of the new ice sounder were made in December 2005 near the Ellsworth Mountains, specifically in Patriot Hills, West Antarctica, at 80°18′ S, 81°22′ W. The complete radar system was installed on a ski-mounted support frame, and towed by hand across the ice surface. The measurement results show the capability of this system to measure ice thickness up to 1000 m and to define internal layers within the ice body.

Modulation of the flow of Rutford Ice Stream, Antarctica, has been reported previously at semi-diurnal, diurnal, fortnightly and semi-annual periods. A model that includes non-linear response to tidal forcing has been shown to fit closely observations at fortnightly frequencies. Here we examine that model further and test its fit at a larger set of observed frequencies, including the large semi-annual displacement. We show analytically that, when forced by major tidal terms, the model (using a basal shear stress exponent m = 3) predicts several discrete response periods from 4 hours to 0.5 years. We examine a 1.5 year GPS record from Rutford Ice Stream and find that the model, when forced by a numerical tide model, is able to reproduce the reported semi-annual signal. We confirm that about 5% of the mean flow is due solely to the (m = 3) non-linear response to tidally varying basal shear stress. Our best-fitting set of model parameters is similar to those originally reported using a much shorter data record, although with noticeably improved fit, suggesting these parameters are robust. We find that m ≈ 3 fits the data well, but that m ≈ 2 does not. Furthermore, we find that a small variation in flow over the 18.6 year lunar node tide cycle is expected. Fits to semi-diurnal and diurnal terms remain relatively poor, possibly due to viscoelastic effects that are not included in the model and reduced GPS data quality at some discrete periods. For comparison, we predict the response of Bindschadler Ice Stream and Lambert Glacier and show, given identical model parameters, a similar response pattern but with ∼1–2 orders of magnitude smaller variability; these may still be measurable and hence useful in testing the applicability of this model to other locations.