A conspicuous precursor of catastrophic ice-shelf break-up along the Antarctic Peninsula, reported widely in the literature, is the gradual increase in surface melting and consequent proliferation of supraglacial lakes and dolines. Here we present analytical and numerical solutions for the flexure stresses within an ice shelf covered by lakes and dolines, both isolated and arrayed. We conclude that surface water promotes ice-shelf instability in two ways: (1) by water-assisted crevasse penetration, as previously noted, and (2) by the inducement of strong tensile flexure stresses (exceeding background spreading stress by 10–100 times) in response to surface water mass loads and ‘hydrostatic rebound’ occurring when meltwater lakes drain.

A large set of ice ablation data from a glacier in the French Alps is used to investigate the sensitivity of ice melting to solar radiation and temperature. The data come from 7 years of observations on a small network of 16 stakes set up on Glacier de Saint-Sorlin. The high spatial variability of ice ablation is shown to be strongly dependent on potential solar radiation. On the other hand, temporal variations are highly correlated with air temperatures measured both at a nearby and at a remote meteorological station, but poorly correlated with incoming shortwave radiation measured at an automatic weather station located close to the glacier terminus. Spatial variations of ice ablation therefore appear to be mainly driven by potential solar radiation, while temporal variations are driven by temperature. This result suggests that minimizing the influence of temperature in an enhanced temperature-index melt model may compromise the model’s ability to simulate interannual changes in melt. These new results may help to improve the calibration, and thus the performance, of these empirical melt models used for long-term simulations of glacier mass balances.

An interpolated bedrock topography is presented for Flask Glacier, one of the tributaries of the remnant part of the Larsen B ice shelf, Antarctic Peninsula. The ice thickness distribution is derived by combining direct but sparse measurements from airborne radio-echo soundings with indirect estimates obtained from ice-flow modelling. The ice-flow model is applied to a series of transverse profiles, and a first estimate of the bedrock is iteratively adjusted until agreement between modelled and measured surface velocities is achieved. The adjusted bedrock is then used to reinterpret the radio-echo soundings, and the recovered information used to further improve the estimate of the bedrock itself. The ice flux along the glacier center line provides an additional and independent constraint on the ice thickness. The resulting bedrock topography reveals a glacier bed situated mainly below sea level with sections having retrograde slope. The total ice volume of 120 ±15 km3 for the considered area of 215 km2 corresponds to an average ice thickness of 560 ± 70 m.

We explore a mathematical model that couples together a thermomechanically evolving subglacial channel, distributed cavity drainage, and basal sliding along a subglacial flood path fed by a jökulhlaup lake. It allows water transfer between channel and cavities and a migrating subglacial water divide or ‘seal’ to form between floods. Notably, it accounts for full coupling between the lake and subglacial drainage in terms of both discharge and pressure, unlike models that neglect the pressure coupling by imposing a known history of lake discharge at the channel inlet. This means that flood hydrographic evolution and its impact on glacier motion are consistently determined by our model. Numerical simulations for a model alpine lake yield stable limit cycles simulating repeating jökulhlaups, with the channel drawing water from the cavities at a varying rate that modulates basal sliding during each flood. A wave of fast sliding propagates down-glacier at flood initiation, followed by deceleration as the growing channel sucks water from the cavities. These behaviours cannot be correctly simulated without the full coupling. We show that the flood’s peak discharge, its initiation threshold and the magnitude of the ‘fast sliding’ wave decrease with the background water supply to the cavities.

Glacier mass balance in Norway is only observed over a small portion (<15%) of the glacierized surface and only for short time periods (<10 years) for most sites. To provide a comprehensive overview of the temporal mass-balance evolution, we modeled surface mass balance for the glacierized area of mainland Norway from 1961 to 2010. The model is forced by operationally gridded daily temperature and precipitation fields which are available at 1 km horizontal resolution from 1957 until the present. The applied mass-balance model accounts for melting of snow and ice by using a distributed temperature-index approach. The precipitation input is corrected to obtain agreement between modeled and observed winter mass balance, and a melt factor and two radiation coefficients are optimized to the corresponding summer balance. The model results show positive trends of winter balance between 1961 and 2000 followed by a remarkable decrease in both summer and winter balances which resulted in an average annual balance of –0.86 ± 0.15 m w.e. a-1 between 2000 and 2010 after four decades of zero to slightly positive annual mass balances.

The turbulent sensible and latent heat fluxes are important components of the surface energy balance over glaciers in the Southern Alps of New Zealand, contributing over half the energy available for ablation during large melt events. To calculate these terms confidently in glacier mass-balance models it is essential to use appropriate parameterizations for surface roughness and atmospheric stability. Eddy covariance measurements at Brewster Glacier were obtained over an ice surface to help facilitate an assessment of the calculation of the turbulent heat fluxes. The roughness length for momentum was found to be 3.6 x 10−3m, while the roughness lengths for temperature and humidity were two orders of magnitude smaller, in agreement with surface renewal theory. A Monte Carlo approach was used to assess the uncertainty in turbulent heat fluxes calculated using the bulk aerodynamic method. It was found that input-data and roughness-length uncertainty could not explain underestimates of observed sensible heat fluxes during periods with low wind speed and large temperature gradients. During these periods a katabatic wind speed maximum alters the formulation of the turbulent exchange coefficient to that typically observed in a neutral atmosphere and this has implications for glacier mass-balance sensitivity.

Surface melt on a glacier can perturb the glaciochemical record beyond the natural variability. While the centre of the Greenland ice sheet is usually devoid of surface melt, many high-Arctic and alpine ice cores document frequent summertime melt events. Current hypotheses interpreting melt-affected ice-core chemistry rely on preferential elution of certain major ions. However, the precise nature of chemistry alteration is unknown because it is difficult to distinguish natural variability from melt effects in a perennially melt-affected site. We use eight trace-element snow chemistry records recovered from Summit, Greenland, to study spatial variability and melt effects on insoluble trace chemistry and physical stratigraphy due to artificially introduced meltwater. Differences between non-melt and melt-affected chemistry were significantly greater than the spatial variability in chemistry represented by nearest-neighbour pairs. Melt-perturbed trace elements, particularly rare earth elements, retained their seasonal stratigraphies, suggesting that trace elements may serve as robust chemical indicators for annual layers even in melt-affected study areas. Results suggest trace-element transport via meltwater percolation will deposit eluted material down-pit in refrozen areas below the nearest-surface chemistry peak. In our experiments, snow chemistry analyses are more sensitive to melt perturbations than density changes or unprocessed near-infrared digital imagery.

Icelandic glaciers cover ∼11 000 km2 in area and store ∼3600 km3 of ice. Starting in 2008 during the International Polar Year, accurate digital elevation models (DEMs) of the glaciers are being produced with airborne lidar. More than 90% of the glaciers have been surveyed in this effort, including Vatnajökull, Hofsjökull, Myrdalsjökull, Drangajökull, Eyjafjallajökull and several smaller glaciers. The publicly available DEMs are useful for glaciological and geological research, including studies of ice-volume changes, estimation of bias in mass-balance measurements, studies of jökulhlaups and subglacial lakes formed by subglacial geothermal areas, and for mapping of crevasses. The lidar mapping includes a 500-1000 m wide ice-free buffer zone around the ice margins which contains many glacio-geomorphological features, and therefore the new DEMs have proved useful in geological investigations of proglacial areas. Comparison of the lidar DEMs with older maps confirms the rapid ongoing volume changes of the Icelandic ice caps which have been shown by mass-balance measurements since 1995/96. In some cases, ice-volume changes derived by comparing the lidar measurements with older DEMs are in good agreement with accumulated ice-volume changes derived from traditional mass-balance measurements, but in other cases such a comparison indicates substantial biases in the traditional mass-balance records.

This study presents a method that allows continuous monitoring of mass balance for remote or inaccessible glaciers, based on repeated oblique photography. Hourly to daily pictures from two automatic cameras overlooking two large valley glaciers in the Swiss Alps are available for eight ablation seasons (2004–11) in total. We determine the fraction of snow-covered glacier surface from orthorectified and georeferenced images and combine this information with simple accumulation and melt modelling using meteorological data. By applying this approach, the evolution of glacier-wide mass balance throughout the ablation period can be directly calculated, based on terrestrial remote-sensing data. Validation against independent in situ mass-balance observations indicates good agreement. Our methodology has considerable potential for the remote determination of mountain glacier mass balance at high temporal resolution and could be applied using both repeated terrestrial and air-/spaceborne observations.

Disintegration of several ice shelves along the Antarctic Peninsula demonstrates a mechanism that involves the conversion of a contiguous ice shelf into an expanding plume of ice-shelf fragments that spreads rapidly across the ocean surface. The growth of surface area and energetic expansion are hypothesized to be driven by gravitational potential energy release associated with iceberg capsize and break-up. Here we investigate this process using a water tank filled with plastic icebergs scaled to represent a laboratory analogue of an expanding plume of ice-shelf fragments (icebergs). Our experiments suggest that hydrodynamic pressure within the water separating neighbouring icebergs is sufficient to couple the motion when their separation is comparable to the iceberg size. This allows one iceberg’s capsize to initiate a ‘domino-like’ effect, where the entire array will subsequently capsize in the same direction and expand across the water surface. Our experimental results motivate the suggestion that cooperative iceberg hydrodynamics is a process that enhances the expansion of ice-shelf fragment plumes during ice-shelf disintegration.

We used measurements of radar-detected stratigraphy, surface ice-flow velocities and accumulation rates to investigate relationships between local valley-glacier and regional ice-sheet dynamics in and around the Schmidt Hills, Pensacola Mountains, Antarctica. Ground-penetrating radar profiles were collected perpendicular to the long axis of the Schmidt Hills and the margin of Foundation Ice Stream (FIS). Within the valley confines, the glacier consists of blue ice, and profiles show internal stratigraphy dipping steeply toward the nunataks and truncated at the present-day ablation surface. Below the valley confines, the blue ice is overlain by firn. Data show that upward-progressing overlap of actively accumulating firn onto valley-glacier ice is slightly less than ice flow out of the valleys over the past ∼1200 years. The apparent slightly negative mass balance (-0.25 cm a-1) suggests that ice-margin elevations in the Schmidt Hills may have lowered over this time period, even without a change in the surface elevation of FIS. Results suggest that (1) mass-balance gradients between local valley glaciers and regional ice sheets should be considered when using local information to estimate regional ice surface elevation changes; and (2) interpretation of shallow ice structures imaged with radar can provide information about local ice elevation changes and stability.

Using June–August surface meteorological measurements for 2004–10 we found that the duration of the summer melt season at elevations of 1800, 1300 and 1000 m a.s.l. on Devon Ice Cap, Canada, increased at a rate of 3.4 (p < 0.1; weak significance), 6.1 (p< 0.01) and 8.8 (p < 0.05) days a−1, respectively. The calculated surface melt rate at 1300m a.s.l. increased from 74 to 133 cm w.e. a−1 between 2007 and 2010. These changes are linked to two types of change in the Arctic atmospheric circulation. Strengthening of the 500 hPa ridge over the Arctic in June–July resulted in increases in both the advection of warm air into the region and the occurrence of cloud-free conditions over the ice cap, causing the available melt energy to increase by 4–24% relative to June–July 2007–10 mean conditions. More frequent southwesterly low-pressure systems in August after 2004 accounted for a 12–38% increase in available melt energy relative to the August 2007–10 daily mean due to advection of warm air into the Arctic, reduced incoming shortwave radiation and increased net longwave radiation. This delayed the timing of freeze-up by an average of 5.5 days a–1 (p < 0.05) at the three sites over 2004–10.

Glaciers on King George Island, Antarctica, have shown retreat and surface lowering in recent decades, concurrent with increasing air temperatures. A large portion of the glacier perimeter is ocean-terminating, suggesting possible large mass losses due to calving and submarine melting. Here we estimate the ice discharge into the ocean for the King George Island ice cap. L-band synthetic aperture radar images covering the time-span January 2008 to January 2011 over King George Island are processed using an intensity-tracking algorithm to obtain surface velocity measurements. Pixel offsets from 40 pairs of radar images are analysed and inverted to estimate a weighted average surface velocity field. Ice thicknesses are derived from simple principles of ice flow mechanics using the computed surface velocity fields and in situ thickness data. The maximum ice surface speeds reach >225 m a-1, and the total ice discharge for the analysed flux gates of King George Island is estimated to be 0.720 ± 0.428 Gt a−1, corresponding to a specific mass loss of 0.64 ± 0.38 m w.e. a-1 over the area of the entire ice cap (1127 km2).

Near-surface air temperature, typically measured at a height of 2 m, is the most important control on the energy exchange and the melt rate at a snow or ice surface. It is distributed in a simplistic manner in most glacier melt models by using constant linear lapse rates, which poorly represent the actual spatial and temporal variability of air temperature. In this paper, we test a simple thermodynamic model proposed by Greuell and Böhm in 1998 as an alternative, using a new dataset of air temperature measurements from along the flowline of Haut Glacier d’Arolla, Switzerland. The unmodified model performs little better than assuming a constant linear lapse rate. When modified to allow the ratio of the boundary layer height to the bulk heat transfer coefficient to vary along the flowline, the model matches measured air temperatures better, and a further reduction of the root-mean-square error is obtained, although there is still considerable scope for improvement. The modified model is shown to perform best under conditions favourable to the development of katabatic winds – few clouds, positive ambient air temperature, limited influence of synoptic or valley winds and a long fetch – but its performance is poor under cloudy conditions.

The Patagonia Icefields are characterized by a large number of outlet glaciers calving into lakes and the ocean. In contrast to the recent intensive research activities on tidewater glaciers in other regions, very few observations have been made on calving glaciers in Patagonia. We analysed satellite images of Glaciar Upsala, the third largest freshwater calving glacier in the Southern Patagonia Icefield, to investigate changes in its front position, ice velocity and surface elevation from 2000 to 2011. Our analyses revealed a clear transition from a relatively stable phase to a rapidly retreating and fast-flowing condition in 2008. The glacier front receded by 2.9 km, and the ice velocity increased by 20–50%, over the 2008–11 period. We also found that the ice surface lowered at a rate of up to 39 m a−1 from 2006 to 2010. This magnitude and the rate of changes in the glacier front position, ice velocity and surface elevation are greater than previously reported for Glaciar Upsala, and comparable to recent observations of large tidewater glaciers in Greenland. Our data illustrate details of a rapidly retreating calving glacier in Patagonia that have been scarcely reported despite their importance to the mass budget of the Patagonia Icefields.

We estimate temperature sensitivity of observed mass balance for glaciers in the European Alps to compare with our earlier estimates from modelling. We treat 1961-90 as our reference period and compare mass balance for 1991-2010 with this reference period. There are eight Alpine glaciers with long enough records. The mean mass balance for the eight glaciers is close to zero for the reference period, very negative for 1991-2010 and highly negative in the exceptionally warm year 2003. Mass-balance changes are compared with 1991-2010 temperature anomalies to give an average temperature sensitivity of approximately -0.7 ± 0.2 m w.e. a-1 K-1 for the eight glaciers. There is considerable variation between individual glaciers, with higher sensitivity for western glaciers and lower sensitivity for eastern glaciers. Temperature sensitivity for the eight glaciers is strongly correlated with the standard deviation of observed mass balance, which from earlier work is known to be correlated with winter balance. The mean temperature sensitivity for observed mass balance is in reasonable agreement with values from earlier models, but changes in glacier hypsometry during 1991-2010 may have further increased observed mass loss.

We report the results of a comparative analysis focusing on grain size, mineralogical composition and spectral reflectance values (400-2500 nm) of cryoconite samples collected from Jakobshavn Isbræ, West Greenland, and Canada Glacier, McMurdo Dry Valleys, Antarctica. The samples from the Greenland site were composed of small particles clumped into larger rounded agglomerates, while those from the site in Antarctica contained fragments of different sizes and shapes. Mineralogical analysis indicates that the samples from Jakobshavn Isbræ contained a higher percentage of quartz and albite, whereas those from Canada Glacier contained a higher percentage of amphibole, augite and biotite. Spectral measurements confirmed the primary role of organic material in reducing the reflectance over the measured spectrum. The reflectance of the samples from the Antarctic site remained low after the removal of organic matter because of the higher concentration of minerals with low reflectance. The reflectance of dried cryoconite samples in the visible region was relatively low (e.g. between ∼0.1 and ∼0.4) favouring increased absorbed solar radiation. Despite high reflectance values in the shortwave infrared region, the effect of the presence of cryoconite is negligible at infrared wavelengths where ice reflectance is low.

The Bering Glacier–Bagley Icefield system in Alaska is currently surging (2011). Large-scale elevation changes and small-scale elevation-change characteristics are investigated to understand surge progression, especially mass transport from the pre-surge reservoir area to the receiving area and propagation of the kinematic surge wave as manifested in heavy crevassing characteristic of rapid, brittle deformation. This analysis is based on airborne laser altimeter data collected over Bering Glacier in September 2011. Results include the following: (1) Maximal crevasse depth is 60 m, reached in a rift that separates two deformation domains, indicative of two different flow regimes. Otherwise surge crevasse depth reaches 20–30 m. (2) Characteristic parameters of structural provinces are derived by application of geostatistical classification. Parameters include significance and spacing of crevasses, surface roughness and crevasse-edge curvature (indicative of crevasse age). A classification based on these parameters serves to objectively discriminate structural provinces, indicative of surge progression down-glacier and up-glacier. (3) Elevation changes from 2011 and 2010 altimetry show 40–70 m surface lowering in the reservoir area in lower central Bering Glacier and 20–40m thickening near the front in Tashalich arm. Combining elevation changes with results of crevasse profilometry and pattern analysis, the rapid progression of the surge can be mathematically–physically reconstructed.

Deriving glacier outlines from satellite data has become increasingly popular in the past decade. In particular when glacier outlines are used as a base for change assessment, it is important to know how accurate they are. Calculating the accuracy correctly is challenging, as appropriate reference data (e.g. from higher-resolution sensors) are seldom available. Moreover, after the required manual correction of the raw outlines (e.g. for debris cover), such a comparison would only reveal the accuracy of the analyst rather than of the algorithm applied. Here we compare outlines for clean and debris-covered glaciers, as derived from single and multiple digitizing by different or the same analysts on very high- (1 m) and medium-resolution (30 m) remote-sensing data, against each other and to glacier outlines derived from automated classification of Landsat Thematic Mapper data. Results show a high variability in the interpretation of debris-covered glacier parts, largely independent of the spatial resolution (area differences were up to 30%), and an overall good agreement for clean ice with sufficient contrast to the surrounding terrain (differences ∼5%). The differences of the automatically derived outlines from a reference value are as small as the standard deviation of the manual digitizations from several analysts. Based on these results, we conclude that automated mapping of clean ice is preferable to manual digitization and recommend using the latter method only for required corrections of incorrectly mapped glacier parts (e.g. debris cover, shadow).