We present an analysis of the active hydrologic system of MacAyeal Ice Stream (MacIS), West Antarctica, from a synthesis of multiple remote-sensing techniques: satellite laser altimetry; satellite image differencing; and hydrologic potential mapping (using a satellite-derived DEM and a bedrock DEM from airborne radio-echo sounding). Combining these techniques augments the information provided by each one individually, and allows us to develop a protocol for studying subglacial hydrologic systems in a holistic manner. Our study reveals five large active subglacial lakes under MacIS, the largest of which undergoes volume changes of at least 1.0 km3. We discuss the hydrologic properties of this system and present evidence for links between the lakes. At least three of the lakes are co-located with sticky spots, i.e. regions of high local basal shear stress. We also find evidence for surface elevation changes due to ice-dynamic effects (not just water movement) caused by changes in basal resistance. Lastly, we show that satellite radar altimetry is of limited use for monitoring lake activity on fast-flowing ice streams with surfaces that undulate on ∼10 km length scales.

Calving from tidewater glaciers and ice shelves accounts for around half the mass loss from both polar ice sheets, yet the process is not well represented in prognostic models of ice dynamics. Benn and others proposed a calving criterion appropriate for both grounded and floating glacier tongues or ice shelves, based on the penetration depth of transverse crevasses near the calving front, computed using Nye’s formula. The criterion is readily incorporated into glacier and ice-sheet models, but has not been fully validated with observations. We apply a three-dimensional extension of Benn and others’ criterion, incorporated into a full-Stokes model of glacier dynamics, to estimate the current position of the calving front of Johnsons Glacier, Antarctica. We find that two improvements to the original model are necessary to accurately reproduce the observed calving front: (1) computation of the tensile deviatoric stress opening the crevasse using the full-stress solution and (2) consideration of such a tensile stress as a function of depth. Our modelling results also suggest that Johnsons Glacier has a polythermal structure, rather than the temperate structure suggested by earlier studies.

Temperate glaciers are more sensitive to climate changes than polar or continental glaciers, and can drive remarkable runoff variation in local water catchments. Here we present recent glacier shrinkage and runoff change for Hailuogou glacier, a typical monsoon temperate glacier on the east slope of Mount Gongga (Minya Konga), China. The surface area of Hailuogou glacier has decreased by 3.5% (0.92 km2) between 1966 (aerial photographs) and 2007 (ASTER images). Flow measurements at a stream gauge about 500 m down-glacier commencing in 1994 display a remarkable increase in annual runoff (mostly during July–September) since 1999. Annual runoff over the same period in a nonglacierized but forested subcatchment (9.17 km2) did not experience significant change. By separating the daily rainfall component from the daily total discharge, monthly catchment water-balance series were calculated for the period 1994–2005, which shows an increasing trend of glacier storage loss. We concluded that air-temperature rise (with a trend of +0.2°C (10 a)−1 between 1988 and 2005, recorded at nearby weather stations) has had an increased effect on glacier mass loss and river runoff change during the past 20 years.

For the first time, to our knowledge, a scientific study of the formation and evolution of waterfall ice, the ephemeral ice structures that form from the freezing of liquid water seeping on steep rock, was performed. We surveyed and analysed three waterfall ice structures near Glacier d’Argentière, Mont Blanc massif, France, between winter 2007 and spring 2009. We reconstruct the global evolution of two vertical ice structures using automatic digital cameras, while the internal ice microstructure was analysed using ice coring and sampling. Macro- and microstructural observations are considered, along with temperature conditions recorded at a nearby meteorological station and directly within the ice structure. They reveal that vertical structures initially grow rapidly from the aggregation of stalactites with microstructures indicative of temperature conditions during their crystallization. After this initial stage, the volume of the ice structure reaches an asymptotic value, as water continues to flow inside the structure, isolated from the outside cold ice; the outer surface remains dry. At the end of the season, the collapse of the free-standing structure does not occur by progressive melting, but is initiated by a horizontal crack propagation at the top. The initiation of this crack seems to be triggered by a drastic temperature decrease.

Due to a lack of in situ measurements, model-based studies of glacier mass balance in the Tibetan Plateau are very limited. An energy-balance model is applied to analyse the mass-balance sensitivity of Xibu glacier, in the Nyainqêntanglha mountain range, to climatic change. A sensitivity calculation shows that a temperature change of ±1°C or a precipitation change of ±35% changes the equilibrium-line altitude (ELA) by 140 ± 125 m. We use a clustering method to link local weather parameters, including wet-season temperatures, to weather types observed over the period 1955–2006. Modelled variability of the Xibu glacier mass balance seems to be controlled by air-temperature variations during the wet season (May–September) and by a long-term warming trend that is unrelated to weather type. The observed wet-season temperature trend of 0.23°C (10 a)−1 leads to an estimated lengthening of the ablation season by 8 days at the glacier terminus (5000 m a.s.l.) and by 23 days at the ELA (5590 m a.s.l.) over the period 1966–2005. The calculated rise in the ELA was 49 m (10 a)−1.

We conducted seismic refraction surveys in the upper ablation area of Storglaciären, a small valley glacier located in Swedish Lapland. We estimated seismic-wave attenuation using the spectral-ratio method on the energy travelling in the uppermost ice with an average temperature of approximately −1 °C. Attenuation values were derived between 100 and 300 Hz using the P-wave quality factor, QP, the inverse of the internal friction. By assuming constant attenuation along the seismic line we obtained mean QP = 6 ± 1. We also observed that QP varies from 8 ± 1 to 5 ± 1 from the near-offset to the far-offset region of the line, respectively. Since the wave propagates deeper at far offsets, this variation is interpreted by considering the temperature profile of the study area; far-offset arrivals sampled warmer and thus more-attenuative ice. Our estimates are considerably lower than those reported for field studies in polar ice (∼500–1700 at −28°C and 50–160 at −10°C) and, hence, are supportive of laboratory experiments that show attenuation increases with rising ice temperature. Our results provide new in situ estimates of QP for glacier ice and demonstrate a valuable method for future investigations in both alpine and polar ice.

We address the usefulness of the unstable manifold correction (UMC) in a Picard iteration for the solution of the velocity field in higher-order ice-flow models. We explain under- and overshooting and how one can remedy them. We then discuss the rationale behind the UMC, initially developed to remedy overshooting, and how it was previously introduced in a Picard iteration to calculate the velocity field in higher-order models. Using a laminar-flow experiment with two higher-order model implementations, we demonstrate that it is not overshooting, but undershooting, that is the main problem when using a proper implementation for the calculation of the velocity field in higher-order models. We also consider a variant of the original UMC algorithm that often enables a relatively fast solution, but is theoretically less sound. Therefore, neither the variant nor the original algorithm is suited for these problems. We present a more appropriate, stable and simple relaxed Picard algorithm and demonstrate that, compared to the true Picard iteration and the variant of the original UMC algorithm, it results in a faster solution of the velocity field in higher-order models for problems with real data.

The air–sea-ice CO2 flux was measured in the ice-covered Saroma-ko, a lagoon on the northeastern coast of Hokkaido, Japan, using a chamber technique. The air–sea-ice CO2 flux ranged from −1.8 to +0.5 mg C m−2 h−1 (where negative values indicate a sink for atmospheric CO2). The partial pressure of CO2 (pCO2) in the brine of sea ice was substantially lower than that of the atmosphere, primarily because of the influence of the under-ice plume from the Saromabetsu river located in the southeastern part of the lagoon. This suggests that the brine had the ability to take up atmospheric CO2 into the sea ice. However, the snow deposited over the sea ice and the superimposed ice that formed from snowmelting and refreezing partially blocked CO2 diffusion, acting as an impermeable medium for CO2 transfer. Our results suggest that the air–sea-ice CO2 flux was dependent not only on the difference in pCO2 between the brine and the overlying air, but also on the status of the ice surface. These results provide the necessary evidence for evaluation of the gas exchange processes in ice-covered seas.

We consider an idealized problem of a sphere of ice growing symmetrically in a spherical cavity within a porous rock in order to identify and quantify different physical mechanisms that can result in fracturing the rock. We show that if the permeability of the rock is very small then high pressures can develop in the cavity as the water inside it expands on freezing. However, given typical permeabilities of most rocks, the pressure is relieved by flow out of the cavity through the rock pores. When ice fills the cavity, there remains a microscopic film of water separating the ice from the rock, owing to disjoining forces, and these forces can stress the rock and have the potential to fracture it. The elastic pressure in the rock depresses the freezing temperature, which can limit the potential for fracturing. This simple example reveals the important interactions between disjoining forces, elasticity and fluid flow in determining the pressure exerted during freezing of water-saturated cavities in rocks.

Gornergletscher, Switzerland, is located adjacent to the marginal lake Gornersee, which periodically drains. We measured glacier ice-flow velocities during two drainage events of the lake, in 2004 and 2006. The common feature of these events was that, during both, Gornersee filled to its maximum level and then overflowed. The events differed in that in 2004 Gornersee rapidly drained via a sudden subglacial connection, whereas in 2006 the lake water continued to overflow and slowly discharged into a nearby moulin. We analysed the changes in ice-flow velocities in the vicinity of Gornersee during the two drainage events, using a three-dimensional ice-flow model which is able to (1) simulate locally variable enhanced basal motion of the ice and (2) account for the load and release of water pressure exerted on the ice margin. We demonstrate that the key features of the observed flow changes can be reproduced adequately in the numerical model by considering these two effects as the main mechanisms. We interpret the 2006 flow changes to be dominated by the release of lake water pressure acting on the ice during the lake drainage. The 2004 ice-flow changes can be explained by enhanced basal motion, and the impact of the lake water pressure provides certain clues to some observations insufficiently explained by enhanced basal motion.

A physically based energy-balance model with improved parameterization of solar radiation for a sloped ice surface has been developed to estimate the backwasting rate of an ice cliff in a debris-covered area. The model has been tested against observations between 5 August and 5 September 2008 on 38 ice cliffs in the debris-covered area of Koxkar glacier, Tuomuer mountain, China. We calculated that the energy-balance model gave a good estimate of the backwasting rates, with errors in the range ±1.96 cm d−1 and root-mean-square errors of 0.99 cm d−1. Errors arising from setting of surface albedo and turbulent flux parameterization were limited. We found that shortwave radiation is the most important heat source for ice-cliff ablation, contributing about 76% of the total heat available for ice melt, while the sensible heat flux provides nearly 24% of the total heat for ice-cliff wastage. The latent heat flux and net longwave radiation are comparatively small according to the model calculation. The mean backwasting rate of ice cliffs in the debris-covered area of Koxkar glacier is estimated at 7.64 m a−1 when the winter ablation is neglected. With this annual backwasting rate and given a mean slope angle of 46.4°, the backwasting of ice cliffs produces about 1.60 × 106 m3 of meltwater, accounting for about 7.3% of the total melt runoff from the debris-covered area.

Supraglacial dust (cryoconite) is an important but poorly understood component of the glacial system. There is a lack of primary data on cryoconite form, extent and dynamics. Here we present a suite of rapid, low-cost methodologies for quantification of granule geometry and supraglacial cryoconite coverage using image data captured by commercially available digital cameras. We develop robust, transferable protocols for analysis of (1) cryoconite granule geometry (major axis, Feret diameter, circularity); (2) centimetre–metre scale supraglacial extent (m2cryoconite m−2surface); and (3) temporal change in supraglacial extent at hourly intervals over several days. Image-processing methodologies were developed using the public domain software ImageJ. Manual (supervised) controls were used to estimate sources of error, and measurements then automated using simple scripting tools (macros). Fully automated processing successfully identified ∼90% of a sample of isolated granules ranging between 2.5 and 39.2 mm, with uncertainties of <20%. Particle sphericity (inferred from circularity) decreased as particle size increased. Supraglacial cryoconite extent was obtained with a mean uncertainty of 37% and 22% for data from field sites in Greenland and Svalbard, respectively. These methods will facilitate acquisition and analysis of datasets for cryoconite across a range of spatial scales, supporting research into cryoconite impacts on supraglacial hydrological connections, nutrient and carbon cycling, and initiation of primary succession in deglaciating environments.

We use topographic maps, historical data, multispectral satellite data and real-time kinematic GPS data to analyze glacier area, length and ice-elevation changes of two glaciers in the central Qilian Shan, China, between 1956 and 2007. We find that the fronts of Yanglonghe glacier No. 1 (5Y432A1) and Yanglonghe glacier No. 5 (5Y432A5) have retreated by 266.5 ± 37.1 m (5.2 ± 0.73 m a−1) and 181.4 ± 37.1 m (3.6 ± 0.73 m a−1) respectively, and that this retreat accelerated after 1999. During the study period, the glacier areas decreased by ∼4.1% and 15.9% respectively. In addition, spatially non-uniform thinning, which averaged 20.2 ± 11 m (0.4 ± 0.22 m a−1) and 16.9 ± 11 m (0.33 ± 0.22 m a−1) in the ablation areas of 5Y432A1 and 5Y432A5 respectively, is observed using digital elevation models constructed using data from 1956, 1977 and 2007. The ice-volume depletion from 5Y432A1 (2.91 × 107 m3) was 2.7 times greater than from the smaller 5Y432A5 (1.08 × 107 m3). Based on records from nearby Tuole weather station, increasing annual temperatures are principally responsible for the observed glacier thinning and retreat.

Using a Monte Carlo (MC) method, we determine the accumulation-rate profile along a flowband, the influx of ice into the upstream end of the flowband and the age of an internal layer. The data comprise the depth profile of the internal layer, a few velocity measurements at the surface and the average accumulation at one location. The data in our example were collected at Taylor Mouth, a flank site off Taylor Dome, Antarctica. We present three alternative formulations of this inverse problem. Depending on the formulation used, this particular inverse problem can have up to four solutions, each corresponding to a different spatial accumulation-rate pattern. This study demonstrates the ability of a MC method to find several solutions to this inverse problem, and how to use a Metropolis algorithm to determine the probability distribution of each of these different solutions. The only disadvantage of the MC method is that it is computationally more expensive than other inverse methods, such as the Gradient method.

A new vibration apparatus for measuring the shear strength of snow has been designed and fabricated. The force applied to a snow block is calculated using Newton’s second law. Results from this apparatus concerning the dependence of the shear strength on snow density, overburden load and strain rate are in reasonable agreement with those obtained from the work of previous researchers. Snow densities ranged from 160 to 320 kg m−3. The overburden load and strain rate ranged from 1.95 × 10−1 to 7.79 × 10−1 kPa and 2.9 × 10−4 to 9.1 × 10−3 s−1 respectively.

Air bubbles in ice cores play an essential role in climate research, not only because they contain samples of the palaeoatmosphere, but also because their shape, size and distribution provide information about the past firn structure and the embedding of climate records into deep ice cores. In this context, we present profiles of average bubble size and bubble number for the entire EDML (Antarctica) core and the top 600 m of the EDC (Antarctica) core, and distributions of bubble sizes from selected depths. The data are generated with an image-processing framework which automatically extracts position, orientation, size and shape of an elliptical approximation of each bubble from thick-section micrographs, without user interaction. The presented software framework allows for registration of overlapping photomicrographs to yield accurate locations of bubble-like features. A comparison is made between the bubble parameterizations in the EDML and EDC cores and data published on the Vostok (Antarctica) ice core. The porosity at the firn/ice transition is inferred to lie between 8.62% and 10.48% for the EDC core and between 10.56% and 12.61 % for the EDML core.

Glacier surfaces support unique microbial food webs dominated by organic and inorganic debris called ‘cryoconite’. Observations from Longyearbreen, Spitsbergen, show how these aggregate particles can develop an internal structure following the cementation of mineral grains (mostly quartz and dolomite) by filamentous microorganisms. Measurements of carbon and dissolved O2 show that these microorganisms, mostly cyanobacteria, promote significant rates of photosynthesis (average 17 μgC g−1 d−1) which assist aggregate growth by increasing the biomass and producing glue-like extracellular polymeric substances. The primary production takes place not only upon the surface of the aggregates but also just beneath, due to the translucence of the quartz particles. However, since total photosynthesis is matched by respiration (average 19 μgC g−1 d−1), primary production does not contribute directly to cryoconite accumulation upon the glacier surface. The microorganisms therefore influence the surface albedo most by cementing dark particles and organic debris together, rather than simply growing over it. Time-lapse photographs show that cryoconite is likely to reside upon the glacier for years as a result of this aggregation. These observations therefore show that a better understanding of the relationship between supraglacial debris and ablation upon glaciers requires an appreciation of the biological processes that take place during summer.