Using the RICE (Radio Ice Cherenkov Experiment) detector at the South Pole, we have estimated the variation in the index of refraction (n) of the firn, as a function of elevation (z) measured from the surface down to z =-150 m. Measurements were made by lowering a dipole transmitter into a dry (5 in (127 mm) caliber) borehole, originally drilled for the RICE experiment in 1998, and determining signal arrival times, as a function of transmitter depth, in the englacial RICE receiver array. We clearly confirm the expected correlation of n(z) with ice density. Our measurements are in fair agreement with previous laboratory characterizations of the dielectric properties of ice cores. These are the first such in situ measurements to be performed at the South Pole.

Observations between 1997 and 2001, of a 30% velocity increase and up to 60 m thinning of downstream parts of Jakobshavn Isbræ, Greenland, immediately following calving of about 4 km of its 15 km floating ice tongue, suggest that acceleration may have been initiated by the calving. Assuming that the force perturbation associated with such weakening is swiftly transmitted far up-glacier, I develop equations to estimate the perturbation. Initially, the observed changes are consistent with the comparatively small perturbation associated with the calving. Thereafter, it was probably sustained by thinning of the remaining ice tongue at rates of about 80 ma–1. Otherwise, the force perturbation would soon have been balanced by reduction in the hydrostatic driving force for longitudinal creep as the glacier thinned, with velocities dropping to their former values. The calculated force perturbation increases to a maximum about 10 km inland of the grounding line, consistent with decreasing weight forces as the glacier thins over bedrock that slopes uphill seawards. Further inland, it progressively decreases, probably because marginal drag increased as the glacier accelerated. Both here and on the floating tongue, marginal ice appears to have been softened by the influence of locally intense shear on ice temperature and/or fabric. More recent observations show continued acceleration and thinning, and most of the remaining ice tongue calved away in April 2003, so thinning is likely to continue.

We estimate Antarctic ice-flow balance velocities, which are the average speeds that ice must flow downslope through a volume assuming that there are equal amounts of ice entering and leaving the ice volume. We use the OSU (Ohio State University) digital elevation model of Antarctica, recent ice accumulation rate data and the BEDMAP ice-thickness data compilations for Antarctica to characterize the physical properties of the ice sheet that are included in the balance-velocity calculation. We adapt a flux algorithm from the hydrology literature that enables us to calculate the flux distribution from any cell in any order. Flux from one cell to its neighbors is partitioned as a simple function of surface slope direction. Digitized flow-stripe directions from satellite images minimize errors in flow direction where surface slopes are low or complex. We estimate errors in balance velocity arising from errors in the data and show semi-quantitatively how properties of the algorithm bias the balance-velocity result. We find a favorable comparison between our model and observed velocity data as well as the balance-velocity patterns reported by other researchers.

The net mass balance on Gulkana Glacier, Alaska, U.S.A., has been measured since 1966 by the glaciological method, in which seasonal balances are measured at three index sites and extrapolated over large areas of the glacier. Systematic errors can accumulate linearly with time in this method. Therefore, the geodetic balance, in which errors are less time-dependent, was calculated for comparison with the glaciological method. Digital elevation models of the glacier in 1974, 1993 and 1999 were prepared using aerial photographs, and geodetic balances were computed, giving – 6.0 ± 0.7 m w.e. from 1974 to 1993 and -11.8 ± 0.7 m w.e. from 1974 to 1999. These balances are compared with the glaciological balances over the same intervals, which were – 5.8 ± 0.9 and -11.2 ± 1.0 m w.e. respectively; both balances show that the thinning rate tripled in the 1990s. These cumulative balances differ by <6%. For this close agreement, the glaciologically measured mass balance of Gulkana Glacier must be largely free of systematic errors and be based on a time-variable area-altitude distribution, and the photography used in the geodetic method must have enough contrast to enable accurate photogrammetry.

On 20 September 2002, an enormous rock/ice slide and subsequent mud-flow occurred on the northern slope of the Kazbek massif, Northern Ossetia, Russian Caucasus. It started on the north- northeast wall of Dzhimarai-Khokh (4780 ma.s.l.) and seriously affected the valley of Genaldon/ Karmadon. Immediate governmental actions, available scientific information, first reconstructions, hazard assessments and monitoring activities as well as initial expert judgments/recommendations are documented in order to enable more detailed analyses and modelling of the event by the wider scientific community. Among the most remarkable aspects related to this event are (1) the relation between the recent event and somewhat smaller but quite similar events that occurred earlier in historical times (1835, 1902), (2) the interactions between unstable local geological structures and complex geothermal and hydraulic conditions in the starting zone with permafrost, cold to polythermal hanging glaciers and volcanic effects (hot springs) in close contact with each other, (3) the erosion and incorporation of a debris-covered valley glacier largely enhancing the sliding volume of rocks, ice, firn, snow, water and probably air to a total of about 100 × 106m3, and (4) the astonishingly high flow velocities (up to 300 km h-1) and enormous length of travel path (18 km plus 15 km of debris/mud-flow). This extraordinary case illustrates that large catastrophic events in high mountain regions typically involve a multitude of factors and require integrated consideration of complex chains of processes, a task which must be undertaken by qualified groups of experts.

We have investigated the formation of 10-50 mm long ice spikes that sometimes appear on the free surface of water when it solidifies. By freezing water under different conditions, we measured the probability of ice-spike formation as a function of: (1) the air temperature in the freezing chamber, (2) air motion in the freezing chamber (which promotes evaporative cooling), (3) the quantity of dissolved salts in the water, and (4) the size, shape and composing material of the freezing vessel. We found that the probability of ice-spike formation is greatest when the air temperature is near -7°C, the water is pure and the air in the freezing chamber is moving. Even small quantities of dissolved solids greatly reduce the probability of ice-spike formation. Under optimal conditions, approximately half the ice cubes in an ordinary ice-cube tray will form ice spikes. Guided by these observations, we have examined the Bally-Dorsey model for the formation of ice spikes. In this model, the density change during solidification forces super cooled water up through a hollow ice tube, where it freezes around the rim to lengthen the tube. We propose that any dissolved solids in the water will tend to concentrate at the tip of a growing ice spike and inhibit its growth. This can qualitatively explain the observation that ice spikes do not readily form using water containing even small quantities of dissolved solids.

A two-dimensional array of electrical conductivity measurements (ECM) has been used to image the acidity of the Siple Dome (Antarctica) ice core in a vertical plane. Annual layering and possible stratigraphic discontinuities are apparent. A brief disruption to the chemical stratigraphy is detected at 680 m. Below 800 m, weaker layering occurs and is interpreted as the result of post-depositional migration of chemical species. This technique provides a way to observe the horizontal continuity of chemical layers in an ice core and identify some types of flow irregularity.

Geophysical investigations on rock glaciers are often difficult because rock glaciers are covered by an unconsolidated debris mantle a few meters thick, are typically <50 m thick and are composed of an ice—rock mixture of unknown composition. Transient electromagnetics (TEM) is a method that allows some of these difficulties to be minimized, and data collection is relatively efficient. TEM, with calibration from terminus exposure, was used to determine the thickness (~60 m) of Fireweed rock glacier, Alaska, U.S.A., under complex valley geometry. A conductive layer beneath the rock glacier was identified, and its distribution is consistent with a till-like layer. Seismic refraction, used to resolve the debris-mantle thickness (2–4 m), suggests the presence of a discontinuity at 18–28 m depth within the rock glacier. The discontinuity is also indicated in the radio-echo sounding and the TEM data, but to a lesser extent. This discontinuity is important because the motion of the rock glacier may occur across this as a “shear plane”.

Ice-spike observations in nature have sparked much interest in the scientific and non-scientific communities alike, yet most research performed thus far has been largely qualitative. We have conducted a quantitative, systematic laboratory investigation in order to assess theories explaining ice-spike growth and to determine the conditions conducive to it. We observed ice-spike growth using time-lapse digital photography, using two water types in two different containers. We observed that ice spikes occurred much more frequently in distilled water than in tap water. Digital images were analyzed to determine the growth rate of the ice spikes. Water temperature was recorded throughout the freezing process, and the cooling rate was used to estimate a bulk heat transfer coefficient. Finally, a simple model, based on mass conservation, was derived and was found to give useful predictions of ice-spike growth rate.

Fireweed rock glacier is a large rock glacier in south central Alaska, U.S.A. It flows relatively fast, with velocities up to 3.5 ma–1, and exhibits both seasonal and annual velocity variations, some of which are related to periodic terminus calving and increased rainfall. Our analysis reveals that motion is likely concentrated in a pseudo-rectangular channel within the larger parabolic channel with a “shear plane” at ~27 m depth. There is likely motion along the shear plane as well as internal deformation above it. We estimate that the ice—rock mixture is up to seven times softer than clean glacier ice with a temperature of –2°C. Calving at the terminus is an important component of the mass balance of this rock glacier.

Hessbreen is a 5km2, polythermal glacier in southern Spitsbergen. It is believed to have surged at the end of the 19th century, and another surge started in 1969, not in 1972 as earlier studies suggested. The surge was probably triggered in the middle section of the glacier and affected the upper part approximately 4 years before including the whole glacier. Present velocities are low, increasing from 0.3m a-1 at the snout to 4.8 m a-1 above the equilibrium line. A distinct velocity reduction occurs where the glacier narrows. This suggests the velocities here are affected by the cold marginal and frontal ice. Large seasonal variations, with summer velocities of more than twice the winter velocities in the middle part of the glacier, indicate a seasonal change in the drainage system in this area. In the lower part there is no apparent sliding. The balance-flux estimate indicates that the actual ice flux is only about 50% of that required for steady state, so Hessbreen appears to be building up in the accumulation area. The surge duration is now estimated to be 8 years.

Geophysical inverse methods are used to calculate the basal motion of a glacier. They are applied to a one-dimensional forward model that can be linearized to make the analysis simpler. The inverse method finds a solution that fits the data within a given error. It selects for smooth solutions to discriminate against unrealistic oscillations. The method is applied to a simple model glacier of uniform shape and thickness to test how well a given basal motion field can be reconstructed. It shows, as expected, that optimizing for smoothness lowers maxima and increases minima of the solution. A step change in basal velocity is drawn out in the inversion over a distance that is given by the half-width of a resolving function. This is typically about three times the ice thickness, but is also affected by the sampling rate of the data. The method is then applied to two glaciers where suitable data are available: Brown Glacier on Heard Island, southern Indian Ocean, and McCall Glacier in the Brooks Range, Alaska, U.S.A. The McCall results agree well with earlier estimates of basal motion.

By using recently developed automatic instruments for fabric and texture measurements on ice, both the c-axis orientation and area of the individual crystals can be determined. Each grain can then be associated with its volume fraction, defined as a function of its measured cross-sectional area, to describe the microstructure of a poly-crystal. The relevance of this approach is studied using a three-dimensional microstructure obtained from the Potts model. In particular, the area weighting is compared to the classical implicit equal weighting used by glaciologists, which assumes that all the grains have the same volume fraction (discrete uniform distribution). Then, using the measurements of c-axis orientation and crystal size performed on the North Greenland Icecore Project (NorthGRIP) ice core, we compare area-weighted and equal-weighted fabrics. All these comparisons are made with respect to the orientation tensor. According to the ability of the Potts model to reproduce the ice microstructure, it is shown that using the grain cross-sectional area to infer its volume fraction improves the description of the actual polycrystal fabric.

Changes in the extent of the polar alpine glaciers within Taylor Valley, Antarctica, are important for understanding past climates and past changes in ice-dammed lakes. Comparison of ground-based photographs, taken over a 20 year period, shows glacier advances of 2-100 m. Over the past ~103 years the climate has warmed. We hypothesize that an increase in average air temperature alone can explain the observed glacier advance through ice softening. We test this hypothesis by using a flowband model that includes a temperature-dependent softness term. Results show that, for a 2˚C warming, a small glacier (50 km2) advances ~25m and the ablation zone thins, consistent with observations. A doubling of snow accumulation would also explain the glacial advance, but predicts ablation-zone thickening, rather than thinning as observed. Problems encountered in modeling glacier flow lead to two intriguing but unresolved issues. First, the current form of the shape factor, which distributes the stress in simple flow models, may need to be revised for polar glaciers. Second, the measured mass-balance gradient in Taylor Valley may be anomalously low, compared to past times, and a larger gradient is required to develop the glacier profiles observed today.

To better understand how internal radar echoes depend on ice-flow conditions and radar polarization, we surveyed two basins in East Antarctica using 179 MHz airborne radar. We compared radar echoes from three ice-flow conditions: parallel sheet flow in the main stream of a basin, convergent flow towards an ice stream, and longitudinal compression by nunataks. We detected a distinct zone of high radar scattering several hundred meters thick at middle depths in the latter two regions. This high-scattering zone was detected only when the radar polarization plane was parallel to the compression axis in ice. Such a high-scattering zone was not found in the parallel-flow region, regardless of the polarization. Using a recently developed theory of radar scattering in ice, we interpret the high-scattering zone as being caused by crystal-orientation-fabric alternations among adjacent ice layers due to difference in horizontal strain components. We argue that the spatial variation of the high-scattering zone is crucial for understanding past and present flow features.

The densification of dry polar snow and firn results in a continuous increase of density with depth accompanied by significant density fluctuations within seasonal layers. Density measurements of high spatial resolution reveal a persistent minimum of density fluctuations in the vicinity of the snow–firn transition (0.55–0.65 g cm-3) in firn-core records. In this study we give an explanation for the fluctuation minimum by applying a new method of X-ray microtomography to obtain three-dimensional (3-D) structural data of a Greenland firn core. At 13 different depths between 10 and 78 m a set of 16 samples of 40 cm total length for each depth interval was measured. A reconstructed firn segment of 40 cm covers 1–2 years of snow accumulation. Using digital image analysis techniques, different structural parameters are estimated including 3-D pore and particle sizes and specific surface areas. It is shown that the densification rates of snow and firn layers consisting of coarse particles are much higher than those consisting of fine particles within the same depth interval. This causes a density crossing of fine- and coarse-grained layers with a minimum of density variations at the crossover point. This crossing-over implies that formerly dense layers in the seasonal density signal are not of the same origin as dense layers in the deeper part of the firn column and that the seasonal density signal will totally change shape with depth. It is speculated that in coarse- and fine-grained firn the dominant mechanism of densification acts over different regimes of density.

We have estimated temperature gradients and melt rates at the bottom of the ice streams in West Antarctica. Measured velocities were used to include the effects of horizontal advection and strain heating in the temperature model and to determine shear heating at the bed. Our modeled temperatures agree well with measured temperatures from boreholes in regions of steady flow. We find that ice-stream tributaries and the inland ice account for about 87% of the total melt generated beneath the Ross ice streams and their catchments. Our estimates indicate that the ice plains of Whillans Ice Stream and Ice Stream C (even when active) have large areas subject to basal freezing, confirming earlier estimates that import of water from upstream is necessary to sustain motion. The relatively low melt rates on Whillans Ice Stream are consistent with observations of deceleration over the last few decades and suggest a shutdown may take place in the future, possibly within this century. While there are pockets of basal freezing beneath Ice Streams D and E, there are larger areas of basal melt that produce enough melt to more than offset the freezing, which is consistent with inferences of relatively steady flow for these ice streams over the last millennium.

We use data from four automatic weather stations (AWSs) in Dronning Maud Land, East Antarctica, to study the surface mass balance and its components. Distinct differences were found between the moisture climates of the high plateau, the katabatic wind zone and the coastal ice shelves: significant undersaturation occurs year-round in the katabatic wind zone, while on the high plateau and on the coastal ice shelf the air is usually close to saturation. In summer, absorption of shortwave radiation at the snow surface enhances surface sublimation at all sites, removing 3-9% of the annual solid precipitation. Significant summer melting is an equally important ablation term near the coast, but vanishes inland. Vertically integrated column drifting-snow sublimation was estimated using two different methods. This process appears to be similar to or greater in magnitude than surface sublimation. Because intervals between significant precipitation events may last as long as several months, sublimation and melt cause extended periods of surface ablation in summer. In summer, all ablation processes together remove 15-56% of the solid precipitation, or 6-27% on an annual basis.

Spatial and temporal variations of the flow of Hansbreen, a tidewater glacier in southern Spitsbergen, Svalbard, are investigated. During summer 1999, surface flow velocities were measured in the ablation zone of Hansbreen with a temporal resolution of 1—2 hours. Short events with strongly increased surface velocities and a typical duration of 1—2 days were observed. These “speed-up events” are related to periods of strongly increased water input to the glacier, due to rainfall or enhanced surface melt. A close relation is found between the surface velocities and water pressure recorded in a moulin. However, there are indications from a short time lag between velocity and water-pressure peak as well as from observed vertical surface uplifts that basal motion is related to basal water storage rather than directly to basal water pressure. The observed short-term velocity variations and associated processes on Hansbreen are very similar to those observed on land-based valley glaciers and suggest that the relevant mechanisms and physical processes that control the flow and its temporal variations are similar. In contrast to the flow of land-based glaciers, sliding velocities on Hansbreen are observed to be high all year round and velocities increase towards the calving front.

A vertical temperature profile through the West Antarctic ice sheet (WAIS) at the summit of Siple Dome reveals an elevated geothermal flux. This could be the root cause for the existence of a dynamic ice-stream system in the WAIS. Siple Dome is still frozen on its bed, but adjacent ice streams have bed temperatures at the pressure-melting point of ice. Although present-day temperature increases due to climatic change do not have an immediate effect on the basal conditions that control the velocity of the ice, indirect effects like a rapid disintegration of the floating ice shelves or additional melt-water input at the surface could give rise to speed-up of the ice streams with an ensuing rise in sea level. Ongoing melt at the base of the ice and changes at the margins will allow continued rapid flow of the ice streams with a possibility of disintegration, within a relatively short period of time, of at least part of the WAIS.