Previous observations have shown that rift propagation on the Amery Ice Shelf (AIS), East Antarctica, is episodic, occurring in bursts of several hours with typical recurrence times of several weeks. Propagation events were deduced from seismic swarms (detected with seismometers) concurrent with rapid rift widening (detected with GPS receivers). In this study, we extend these results by deploying seismometers and GPS receivers in a dense network around the tip of a propagating rift on the AIS over three field seasons (2002/03, 2004/05 and 2005/06). The pattern of seismic event locations shows that icequakes cluster along the rift axis, extending several kilometers back from where the rift tip was visible in the field. Patterns of icequake event locations also appear aligned with the ice-shelf flow direction, along transverse-to-rift crevasses. However, we found some key differences in the seismicity between field seasons. Both the number of swarms and the number of events within each swarm decreased during the final field season. The timing of the slowdown closely corresponds to the rift tip entering a suture zone, formed where two ice streams merge upstream. Beneath the suture zone lies a thick band of marine ice. We propose two hypotheses for the observed slowdown: (1) defects within the ice in the suture zone cause a reduction in stress concentration ahead of the rift tip; (2) increased marine ice thickness in the rift path slows propagation. We show that the size–frequency distribution of icequakes approximately follows a power law, similar to the well-known Gutenberg–Richter law for earthquakes. However, large icequakes are not preceded by foreshocks nor are they followed by aftershocks. Thus rift-related seismicity differs from the classic foreshock and aftershock distribution that is characteristic of large earth quakes.

We studied the diffusion rate of stable isotopes of water ice in a controlled laboratory experiment, where isotopically different layers of fabricated firn were stacked. This experiment was done in order to study the diffusion rate between the isotopically different layers and to compare the measurements with forward-calculated diffusion rates of the fabricated firn stack, using three different analytical and numerical methods. The layers varied in thickness within the stack, and the studied period includes variations in firn temperature. We sampled the firn stack on five different dates along the 144 day long experiment. We obtained a quantitatively good match between the calculated and measured diffusion rates, with a root-mean-square similarity of 67–80% in the depth-averaged diffusion rates. Despite the good average match, we found some qualitative discrepancy between our calculations and our measurements. The calculated diffusion rates predict >5% faster isotopic decay of thinner layers and 13% slower decay of thicker layers than our measurements show. We speculate that grain size and grain fabric may influence the tortuosity of the firn pack, and that the current generalization of tortuosity, being a function of firn density as the only free parameter, is an oversimplification.

Analysis of aerial photographs, GPS mapping and comparison of digital elevation models have been used to quantify the losses in extent and volume observed (1981–2005) in the glaciers of the Maladeta massif (Spanish Pyrenees). The data are examined in relation to different climatic (temperature, precipitation) and topographic factors that control glacial retreat both at the general and local scales. The evolution observed in the massif is characterized by the remarkable surface and volumetric shrinkage registered in all the glaciers: (1) glacial ice decreased in extent 35.7%, reducing from 240.9 ha to 155.0 ha; (2) losses in total ice volume reached 0.0137 km3 (75.6 m w.e.); (3) mean altitude of the studied glaciers increased 43.5 m. These changes seem to have been forced by climatic change (reduction in the snowfall contributions and increase in the maximum temperatures) during the past few decades in this Pyrenean region. In addition, local variables such as the orientation of each glacier, their altitude and their initial size seem to have induced significant spatial differences in the magnitude of the losses.

Global positioning system (GPS) campaigns were conducted during the 2003 and 2004 austral summer seasons to obtain insight into the velocity and strain-rate distribution on Schirmacher Glacier, central Dronning Maud Land, East Antarctica. GPS data were collected at 21 sites and analyzed to estimate the site coordinates, baselines and velocities. The short-term precision of the base station, MAIT, is estimated from the daily coordinate repeatability solutions during the two years. All GPS points on the glacier were constrained with respect to MAIT and nearby International GPS Service stations. Horizontal velocities of the glacier sites lie between 1.89 ± 0.01 and 10.88 ± 0.01 ma−1 to the north-northeast, with an average velocity of 6.21 ± 0.01 m a−1. The principal strain rates provide a quantitative measurement of extension rates, which range from (0.11 ± 0.01) × 10−3 to (1.48 ± 0.85) × 10−3a−1, and shortening rates, which range from (0.04 ± 0.02) × 10−3 to (0.96 ± 0.16) × 10−3a−1. The velocity and strain-rate distributions across the GPS network in Schirmacher Glacier are spatially correlated with topography, subsurface undulations, fracture zones/crevasses and the partial blockage of the flow by nunataks and the Schirmacher Oasis.

The deformation, resulting from a surge in 1985, of Glaciar Horcones Inferior is analyzed using structural geological models. During the surge, previously continuous debris cover was deformed by the formation of regularly separated and rotated ice blocks, suggesting a system of linked rotational extensional faults. Block tilting was measured from photographs taken shortly after the surge, showing rotation of the debris-covered surface. Fault inclination was assumed to be coincident with the debris-free side of the block. Glacier advance during the surge was obtained by comparing pre-surge aerial photographs with the position of maximum advance after the surge. Glacier thinning was estimated from the debris surface average lowering (relief generated at lateral scarps coincident with shear zones) and ice thickness measurements after surge termination. Three independent sets of information, geometry of the deformation (i.e. depth of detachment, fault traces, fault spacing, block rotation), glacier thinning and net advance, limit possible interpretations. Surface geometry suggests a domino-style or a linked planar rotational extensional fault system. In the observed configuration, however, these models can only explain a 12–13% extension. Glacier thinning suggests 30% local extension, and total glacier advance implies 16% minimum extension, which does not account for some frontal compression, as observed. A linked curved rotational extensional fault model fits the data well, implying a significant degree of internal deformation within each block. This model satisfactorily explains the observed deformation produced by the surge. It may also explain some modes of fast glacier flow, since the observed style of block tilting is present in other glaciers with high relief.

The quality of the ice core from Siple Dome, West Antarctica, varied widely, with significant fracturing below 400 m. Bubbly ice persisted to the ice–rock interface at 1004 m and constituted the brittle zone. The core has undergone minimal relaxation and has remained brittle and prone to fracturing more than 5 years after it was drilled. This behavior is attributed to unrelieved stresses from Kamb and Bindschadler Ice Streams (former Ice Streams C and D) flanking the dome. Melt layers were identified sporadically throughout the core, as were inclined layers tilted at angles that occasionally exceeded 10°. Structurally, the ice was characterized by extensive recrystallization including grain-size changes from 0.074 cm2 at 59 m to >50 cm2 at 992 m, and major transitions in c-axis fabrics. Unusual fabrics included vertical c-axis clusters superimposed on vertical girdles that may reflect vertical compression acting in conjunction with horizontal tension. The sudden appearance of a shear-type fabric at 700–800 m appears closely linked to the occurrence of abundant tephra particles embedded in the ice. The occurrence of dispersed sediment in the bottom 2 m is attributed to freeze-on of basal meltwater.

Approximately 300 volcanic ash and dust layers were observed in the Siple Dome (Antarctica) ice core. Most of this tephra, deposited between 700 and 800 m depth, consisted primarily of glass shards with varying amounts of crystalline material and groundmass fragments. The pattern of distribution of tephra fallout closely replicates that found in the Byrd ice core, indicative of contemporaneous deposition at both locations. Peak fallout occurred approximately 19 500 years ago, based on methane tie points in the Siple Dome and Greenland Ice Sheet Project 2 (GISP2) ice cores. Mount Berlin was identified as a potential source of tephra, although other volcanoes in West and East Antarctica appear to have contributed ash and dust. Ice between 697 and 730 m, in which fine-grained tephra is concentrated, has undergone enhanced thinning compared to ice with a similar concentration of tephra deposited contemporaneously between 1300 and 1540 m at Byrd. It is speculated that this thinning has occurred in response to dynamic interaction between ice at Siple Dome and the two ice streams flanking it. A dramatic change to a shear fabric appears to be directly related to the higher concentration of volcanic particles in the ice between 700 and 800 m.

Ice-sheet flow solutions commonly adopt a simple isotropic viscous law, in which the deviatoric stress is coaxial with the strain rate, and the single response function depends on only one invariant, the second principal invariant of the deviatoric stress. This can be correlated with single-stress component tests, which cannot, however, verify the validity of the simplification. Morland (2007) has shown how combined compression and shear tests can verify and correlate a general quadratic isotropic viscous relation, or simply third invariant dependence, and there is evidence that at least third principal invariant dependence is required. Morland (2007) showed that a significant quadratic term changes the relative stress magnitudes in the reduced model for ice-sheet flow, that crucial simplifications are not achieved and formally noted that dependence on the third invariant in the coaxial relation also prevented the simplification. It is now shown that, provided the third invariant dependence does not dominate the isotropic relation, further asymptotic expansion does yield the reduced model simplification, and the influence of different weightings of third invariant dependence is illustrated by a comparison of steady radial flow solutions over a flat bed. It is found that the third invariant dependence does not have a large influence in these examples.

Little is known about the Himalayan glaciers, although they are of particular interest in terms of future water supply, regional climate change and sea-level rise. In 2002, a long-term monitoring programme was started on Chhota Shigri Glacier (32.2° N, 77.5° E; 15.7 km2, 6263–4050 ma.s.l., 9 km long) located in Lahaul and Spiti Valley, Himachal Pradesh, India. This glacier lies in the monsoon–arid transition zone (western Himalaya) which is alternately influenced by Asian monsoon in summer and the mid-latitude westerlies in winter. Here we present the results of a 4 year study of mass balance and surface velocity. Overall specific mass balances are mostly negative during the study period and vary from a minimum value of –1.4 m w.e. in 2002/03 and 2005/06 (equilibrium-line altitude (ELA) ∼5180 m a.s.l.) to a maximum value of +0.1 m w.e. in 2004/05 (ELA 4855 m a.s.l.). Chhota Shigri Glacier seems similar to mid-latitude glaciers, with an ablation season limited to the summer months and a mean vertical gradient of mass balance in the ablation zone (debris-free part) of 0.7mw.e.(100 m)–1, similar to those reported in the Alps. Mass balance is strongly dependent on debris cover, exposure and the shading effect of surrounding steep slopes.

The relationship between velocity field and isochrone geometry along a steady flowline of an ice sheet is examined. The method is analytical and based upon the stream function and its vertically normalized form, the normalized stream function (NSF). We show that the slope of the isochrones is the slope of the iso-NSF lines, plus a path term which is the cumulative result of the past trajectory of the ice particles. We illustrate this path term in three different examples: varying basal melting, varying basal sliding (Weertman effect) and varying velocity profile around a divide (Raymond effect). The path term generally counteracts the slope of the iso-NSF lines. In the case of the Raymond effect, it can even lead to depressions surrounding the bumps if the transition from dome to flank velocity profile is sufficiently abrupt.

Snow pits sampled during two consecutive years (2001, 2002) at the summit of Lomonosovfonna ice cap in central Spitsbergen, Svalbard, showed that ion concentrations were spatially homogeneous. The snowpack on Lomonosovfonna shows no evidence of aerosol deposition from Arctic haze, in contrast to Holtedahlfonna (a glacier at a similar altitude in northern Spitsbergen) where there is a clear signature. In common with many other ice caps in the Arctic, Lomonosovfonna experiences periodic melting, and the deepest of the snow pits contained a record of one exceptionally warm (2001) and one long summer (2000). The most easily eluted species are nitrate and the divalent ions. Very low ion concentrations and high values of a melt indicator log([Na+]/[Mg2+]) were a result of either deep percolation or runoff of ions during melting. Comparing the snow-pit record with the ion record of more than 800 years from an ice core drilled on Lomonosovfonna in 1997 reveals some layers with similar composition to those that suffered significant melting in the snowpack: a few years in the 20th century and around AD 1750, and all of the core from before AD 1200 show unusually heavy melting.

A record of accumulation and ablation from a network of 47 stakes at a coastal blue-ice area in Terre Adélie, Antarctica, is presented and analyzed. The record covers early 2004 to early 2006, from 25 field surveys including some in austral winter. The two years are very different, with a virtually null surface mass balance during the 2004 winter but large accumulation during the 2005 winter. A snow/ice energy- and mass-balance model is used to reproduce the accumulation and ablation record. A parameterization for snow erosion by wind is included. Input meteorology is from the European Centre for Medium-Range Weather Forecasts (ECMWF) analyses and forecasts, corrected using 1 year of local meteorological observations from an automatic weather station. Model results agree reasonably well with the observations. Wind erosion is the largest contributor to ablation, removing much of the precipitation. Sublimation and, to a lesser extent, melt/runoff together account for >60 cm w.e. of ablation in 2 years, mainly in summer. Although the record is short, it confirms high interannual variability and thus high sensitivity to meteorology and climate. Monitoring and understanding the mass balance of such coastal blue-ice areas may help monitor and detect climate change in the Antarctic coastal regions.

The University of Maine Ice Sheet Model was used to study basal conditions during retreat of the Laurentide ice sheet in Maine. Within 150 km of the margin, basal melt rates average ∼5 mm a−1 during retreat. They decline over the next 100 km, so areas of frozen bed develop in northern Maine during retreat. By integrating the melt rate over the drainage area typically subtended by an esker, we obtained a discharge at the margin of ~1.2 m3 s-1. While such a discharge could have moved the material in the Katahdin esker, it was likely too low to build the esker in the time available. Additional water from the glacier surface was required. Temperature gradients in the basal ice increase rapidly with distance from the margin. By conducting upward into the ice all of the additional viscous heat produced by any perturbation that increases the depth of flow in a flat conduit in a distributed drainage system, these gradients inhibit the formation of sharply arched conduits in which an esker can form. This may explain why eskers commonly seem to form near the margin and are typically segmented, with later segments overlapping onto earlier ones.

A diagnostic, dynamic/thermodynamic ice-shelf model is applied to the George VI Ice Shelf, situated in the Bellinghausen Sea, Antarctica. The George VI Ice Shelf has a peculiar flow geometry which sets it apart from other ice shelves. Inflow occurs along the two longest, and almost parallel, sides, whereas outflow occurs on the two ice fronts that are relatively short and situated at opposite ends of the ice shelf. Two data sources were used to derive the ice thickness distribution: conventional radioecho sounding from the British Antarctic Survey was combined with thickness inferred from surface elevation obtained by the NASA GLAS satellite system assuming hydrostatic equilibrium. We simulate the present ice flow over the ice shelf that results from the ice thickness distribution, the inflow at the grounding line and the flow rate factor. The high spatial resolution of the ice thickness distribution leads to very detailed simulations. The flow field has some extraordinary elements (e.g. the stagnation point characteristics resulting from the unusual ice-shelf geometry).

We present seismic and radar data from the onset region of Rutford Ice Stream, West Antarctica, which show the form and internal structure of a variety of bed forms beneath an active ice stream. The ice flow in the area of our survey accelerates from 72 to >200 m a–1, the ice is 2200–3200 m thick, and the bed of the ice stream lies up to 2000 m below present sea level. We have imaged the internal structure of the bed forms with seismic reflection techniques and also observed radar reflections from below the bed in some circumstances. We observed a transverse moraine 2 km wide and 1.5 km long beneath the slower-flowing part of the ice stream, which we interpret to be composed of unconsolidated sediment undergoing active deformation near the ice–sediment interface. We observed drumlins of classical form with elongation ratios of between 1:1.5 and 1:4.0 where the surface flow speed exceeded 95 m a–1. The conformity of the internal structure of the bed forms with the ice base suggests that the bed forms are active depositional features in congruence with the observation of a contemporary drumlin-forming episode in the distal part of the same ice stream. These observations provide the first direct evidence of the association between ice-stream flow speed and bed-form shape.

Glacier and lake variations in the Yamzhog Yumco basin in southern Tibet were studied by integrating series of spatial data from topographic maps and Landsat images at three different times: 1980, 1988/90 and 2000. The results indicate that the total glacier area has decreased from 218 km2 in 1980 to 215 km2 in 2000, a total reduction of 3 km2 (i.e. a 1.5% decrease). Glacier recession rates were clearly larger in the 1990s than the 1980s due to the warmer climate. The total lake area decreased by about 67 km2 during 1980–90 and increased by 32 km2 during 1990–2000. It is suggested that change of lake area in the basin was rapid and most likely caused primarily by the change in precipitation and evaporation in the basin, and secondarily by the increased water supply from melting glaciers.

A measure of snow density is required to estimate water equivalent ice-sheet surface mass balance (SMB) from stake measurements. Previous studies have utilized the snow density at different depths within the snow. By considering the snow densification process in the time interval between stake height measurements, we find that use of the snow density at the base of the stake is more appropriate. We assume the stakes are firmly anchored at the bottom and that Sorge’s law holds, i.e. the density–depth profile does not change with time. Applying this method to the data for 36 snow stakes on Dome Fuji, the SMB in 2003 was 36.5 kg m-2 a−1, 27% larger than the previous estimate, which used surface snow density. Correct selection of the snow density for SMB estimations is important, especially for Antarctic inland areas where accumulation is low (e.g. Dome Fuji, Vostok, Dome C and South Pole) and where the snow density near the surface varies markedly.

Samples of snow and firn from accumulation zones on Clark, Commonwealth, Blue and Victoria Upper Glaciers in the McMurdo Dry Valleys (∼77–78º S, 161–164º E), Antarctica, are evaluated chemically and isotopically to determine the relative importance of local (site-specific) factors vs regional-scale influences in defining glaciochemistry. Spatial variation in snow and firn chemistry confirms documented trends within individual valleys regarding major-ion deposition relative to elevation and to distance from the coast. Sodium and methylsulfonate (MS–), for example, follow a decreasing gradient with distance from the coast along the axis of Victoria Valley (350–119 µg L−1 for Na+; 33–14 µg L−1 for MS–); a similar pattern exists between Commonwealth and Newall Glaciers in the Asgaard Range. When comparing major-ion concentrations (e.g. Na+, MS−, Ca2+) or trace metals (e.g. Al, Fe) among different valleys, however, site-specific exposures to marine and local terrestrial chemical sources play a dominant role. Because chemical signals at all sites respond to particulates with varying mixtures of marine and terrestrial sources, each of these influences on site glaciochemistry must be considered when drawing temporal climate inferences on regional scales.

The spatial pattern of accumulation rate can be inferred from internal layers in glaciers and ice sheets. Non-dimensional analysis determines where finite strain can be neglected (‘shallow-layer approximation’) or approximated with a local one-dimensional flow model (‘local-layer approximation’), and where gradients in strain rate along particle paths must be included (‘deep layers’). We develop a general geophysical inverse procedure to infer the spatial pattern of accumulation rate along a steady-state flowband, using measured topography of the ice-sheet surface, bed and a ‘deep layer’. A variety of thermomechanical ice-flow models can be used in the forward problem to calculate surface topography and ice velocity, which are used to calculate particle paths and internal-layer shapes. An objective tolerance criterion prevents over-fitting the data. After making site-specific simplifications in the thermomechanical flow algorithm, we find the accumulation rate along a flowband through Taylor Mouth, a flank site on Taylor Dome, Antarctica, using a layer at approximately 100 m depth, or 20% of the ice thickness. Accumulation rate correlates with ice-surface curvature. At this site, gradients along flow paths critically impact inference of both the accumulation pattern, and the depth-age relation in a 100 m core.

A number of improvements have been made to an existing two-dimensional ice-flow model applied to an alpine glacier. Analysis of the results of the existing model revealed several shortcomings. The first concerns the lack of mass conservation of the applied alternating-direction-implicit (ADI) scheme. A semi-implicit (SI) scheme is therefore proposed and the effects on mass conservation assessed by a comparison with the ADI scheme. The comparison is first carried out with a simple theoretical glacier for which the improvement is significant. Concerning the real case of Glacier de Saint-Sorlin, France, the initial deviation in mass conservation was much less pronounced such that the new scheme, although improving mass conservation, does not significantly change the modelled dynamics. However, other shortcomings that have a more profound impact on the modelling of glacier behaviour have been identified. The ice thickness may become negative over some gridpoints, leading to an inconsistency. The problem is partly resolved by incorporating extra checks on critical gridpoints at the glacier border. Finally, with the help of ice particle tracking, unrealistic ice settlement above the bergschrund has been identified as the main reason for spurious dynamic effects and has been corrected.