To gain new insight into the mechanisms of basal motion, we have demonstrated the feasibility of an active seismic technique to measure temporal changes in basal conditions on sub-hourly time-scales. One region of the bed of Black Rapids Glacier, Alaska, U.S.A., was monitored for a period of 45 days using seismic reflections. The majority of these reflections were nearly identical. However, three significant anomalies were recorded several days apart. These corresponded with the englacial drainage of two ice-marginal lakes and one supraglacial pothole, each up-glacier of the study site, as well as dramatic increases in basal motion. Two of these seismic anomalies revealed identical changes over 1 km2 of the bed despite the fact that their drainage events occurred at different locations. Further, these two seismic anomalies were followed by records identical to the non-anomalous state, showing that the seismic changes were reversible. In one of these events, two records taken 36 min apart revealed that the transition between the anomalous and normal states occurred completely within this short interval.

Until now, an assumption of surface-parallel glacier flow has been used to express the vertical velocity component in terms of the horizontal velocity vector, permitting all three velocity components to be determined from synthetic aperture radar interferometry. We discuss this assumption, which neglects the influence of the local mass balance and a possible contribution to the vertical velocity arising if the glacier is not in steady state. We find that the mass-balance contribution to the vertical surface velocity is not always negligible as compared to the surface-slope contribution. Moreover, the vertical velocity contribution arising if the ice sheet is not in steady state can be significant. We apply the principle of mass conservation to derive an equation relating the vertical surface velocity to the horizontal velocity vector. This equation, valid for both steady-state and non-steady-state conditions, depends on the ice-thickness distribution. Replacing the surface-parallel-flow assumption with a correct relationship between the surface velocity components requires knowledge of additional quantities such as surface mass balance or ice thickness.

Using changes observed in daily seismic reflections, we have investigated the basal morphology of Black Rapids Glacier, Alaska, U.S.A. The englacial drainage of ice-marginal lakes caused significant changes in the daily reflections, as well as dramatic increases in basal motion. Changes in reflection arrival times and amplitudes indicate that there is a basal till layer at least 5 m thick at some locations beneath this surge-type glacier. Rapid changes in the observed reflection coefficients during the drainage events indicate that changes in till properties must occur throughout the entire 5 m thick layer, they must last for several days following the lake drainages and they must be completely reversible over as little as 36 min. Our seismic analysis shows that changes in effective pressure of the till are unlikely to cause the required changes in the reflection coefficients, but that a decrease in till saturation is likely. We therefore interpret the cause of the seismic anomalies as being a temporary decrease in saturation as water is input to the subglacial hydraulic system, and propose that such a change may occur quickly and reversibly by a redistribution of overburden pressure. Higher water pressures within the hydraulic system cause that region to support more of the glacier’s weight, leaving the remaining areas to support less. Any till within these areas of decreased normal stress would experience a consequent decrease in pore-water pressure, causing gas to exolve, thus decreasing saturation. This decrease in saturation would cause a change in the strength of the basal layer and may affect basal dynamics.

Repeat-pass L-band interferometric synthetic aperture radar (InSAR) data for part of Hielo Patagónico Sur, Chile, were collected by the space-shuttle-based Spaceborne Imaging Radar C (SIR-C) over a 4 day span in October 1994. Three co-registered complex SAR images are used to generate phase-coherence maps, a digital elevation model (DEM) and an ice-velocity map. The phase-coherence maps indicate low coherence in the 5–15 km approaching the termini due to large displacements, ice deformation and melting. However, the coherence is high over nearly all of the remaining imaged icefield. Ice-velocity precision is greater than 2 cm d−1, while the DEM is good to about 25 m. A flow divide between two of the glaciers is mapped by locating a narrow band of near-zero ice velocity. Horizontal ice-surface velocity profiles calculated along flowlines show there is a high degree of spatial variability reaching a peak value of 5.5 m d−1 located 3.5 km from the terminus of Glaciar Europa. Longitudinal strain rates along the center lines calculated from these velocities at the locations of the initiation of crevassing are used to compute the tensile strength of ice (169–224 kPa).

Snow glide is the translational slip of the entire snowpack over a sloping ground surface, and it is thought that rapid rates of snow glide precede full-depth avalanches. The nature of avalanches that release at the ground makes them difficult to predict and difficult to control using explosives.

On-slope instrumentation comprised of stainless-steel "glide shoes" was used to measure rates of snow glide for two winters on a bedrock slope adjacent to the Coquihalla Highway, Cascade Mountains, British Columbia, Canada. Climate data and avalanche occurrences were recorded by the British Columbia Ministry of Transportation and Highways.

Our results show that the supply of free water to the snow/ground interface by rain or snowmelt is the most important influence on full-depth avalanche release. Full-depth avalanche release responds to rainfall and snowmelt events within 12-24 hours. Occasionally, full-depth avalanches occur unexpectedly during clear, cold periods. Snowmelt by radiation is thought to contribute enough meltwater during these cold periods to induce higher rates of snow glide and full-depth avalanche release. The results also indicate that snow glide alone is not a reliable indicator for full-depth avalanche release.

Highly elongated bubbles are sometimes observed in ice-sheet ice. Elongation is favored by rapid ice deformation, and opposed by diffusive processes. We use simple models to show that vapor transport dominates diffusion except possibly very close to the melting point, and that latent-heat effects are insignificant. Elongation is favored by larger bubbles at pore close-off, but is nearly independent of bubble compression below close-off. The simple presence of highly elongated bubbles indicates only that a critical ice-strain rate has been exceeded for significant time, and provides no information on possible disruption of stratigraphic continuity by ice deformation.

A digital image-processing approach is proposed which allows the extraction of two-dimensional polycrystalline ice microstructure (grain boundaries) from thin sections observed between cross-polarisers. It is based on image segmentation of colour images. The method is applied to the preliminary analysis of the shallow ( Holocene) ice of the European Project for Ice Coring in Antarctica (EPICA) ice core at Dome Concordia. Structural parameters, such as the mean cross-sectional area, shape anisotropy and grain morphology, are obtained. The interest and limitations of this automatic procedure are discussed.

Ice at depth in ice sheets can be softer in bed-parallel shear than Glen’s flow law predicts. For example, at Dye 3, Greenland, enhancement factors of 3 4 are needed in order to explain the rate of borehole tilting Previous authors have identified crystal fabric as the dominant contributor, but the role of impurities and crystal size is still incompletely resolved. Here we use two formulations of anisotropic flow laws for ice (Azuma’s and Sachs’ models) to account for the effects of anisotropy, and show that the measured anisotropy of the ice at Dye 3 cannot explain all the detailed variations in the measured strain rates, the jump in enhancement across the Holocene–Wisconsin boundary is larger than expected from the measured fabrics alone. Dust and soluble-ion concentration divided by crystal size correlates well with the residual enhancement, indicating that most of the “excess deformation” may be due to impurities or crystal size. While the major features of the deformation at Dye 3 are explained by anisotropy and temperature, results also suggest that further research into the role of impurities and crystal size is warranted.

Measurement of glacier surface velocity provides some constraint on glacier flow models used to date ice cores recovered near the flow divide of remote high-altitude ice caps. The surface velocity is inferred from the change in position of a network of stakes estimated from the least-squares adjustment of geodetic observations – terrestrial and/or spaced-based – collected approximately 1 year apart. The lack of outliers in and the random distribution of the post-fit observation residuals are regarded as evidence that the observations contain no blunders. However, if the network lacks sufficient geometric redundancy, the estimated stake positions can shift to fit erroneous observations. To determine the maximum size of these potential undetected shifts, given the covariance of the observations and the approximate network geometry, expressions are developed to analyze a network for redundancy number and marginally detectable blunders (internal reliability), and the position shifts from marginally detectable blunders (external reliability). Two stake networks, one on the col of Huascar–n (9°–07'S, 77°–37'W; 6050 m a.s.l.) in the north-central Andes of Peru and one on the Guliya ice cap (35°–17'N, 81°–29'E; 6200 m a.s.l.) on the Qinghai–Tibetan Plateau in China, are examined for precision and internal and external reliability.

Thermal conductivity of snow has been investigated experimentally using the thermal-probe method, which is a transient method of measurement. The measurements have been made over a wide range of snow density (for fresh and dense snow), for varying temperatures and for different conditions of water content, snow-grain type, etc., both in the field and in the laboratory. The results are presented along with detailed sample descriptions. Thermal conductivity of snow increases with density and water content. It also increases with temperature, and the effect is more pronounced for temperatures between –15° and 0°C.

Slushflows were first recorded in the Iberian Peninsula on 18 December 1997. Three slushflows were released at the ski resort of El Port del Comte, in the Catalan Pyrenees, northeast Spain, during intense rainfall. Two of the slushflows originated on the pistes, and the third affected another piste. Three ski lifts were damaged. This paper analyzes the hydrogeological characteristics of the massif, the geomorphic features of the terrain and the meteorological and snowpack conditions that caused the release of the slushflows. Man's role in triggering the slushflows by compacting snow on the pistes is also considered. Drainage control for reducing the hazard is outlined, taking into consideration the low frequency of the phenomenon.

A mass-balance model based on the energy balance at the snow or ice surface is formulated, with particular attention paid to processes affecting absorption of radiation. The model is applied to a small glacier, Glacier AX010 in the Nepalese Himalaya, and tests of its mass-balance sensitivity to input and climatic parameters are carried out. Calculated and observed area-averaged mass balances of the glacier during summer 1978 (June-September) show good agreement, namely -0.44 and -0.46 m w.e., respectively.

Results show the mass balance is strongly sensitive to snow or ice albedo, to the effects of screening by surrounding mountain walls, to areal variations in multiple reflection between clouds and the glacier surface, and to thin snow covers which alter the surface albedo. In tests of the sensitivity of the mass balance to seasonal values of climatic parameters, the mass balance is found to be strongly sensitive to summer air temperature and precipitation but only weakly sensitive to relative humidity.

During the 1992 summer field season we installed arrays of “plough-meters” and water-pressure transducers beneath Trapridge Glacier. Yukon Territory, Canada, to study hydromechanical coupling at the ice–bed interface. Diurnal signals recorded with two of these ploughmeters appear to correlate with fluctuations in sub-glacial water pressure. These diurnal variations can be explained by changes in basal resistance to sliding as mechanical conditions at the bed vary temporally in response to changes in the subglacial hydrological system. We propose that a lubricating water film, associated with high water pressures, promotes glacier sliding, whereas low pressures cause increased basal drag resulting in “sticky” areas. Using a theoretical model, we analyze the sliding motion of glacier ice over a flat surface having variable basal drag and show that a consistent explanation can be developed. Results from our model calculations provide strong support for the existence of time-varying sticky spots which are associated with fluctuations in subglacial water pressure.

Measurements of ice thickness and surface elevation are prerequisite to many glaciological investigations. A variety of techniques has been developed for interpretation of these data, including means of constructing regularly gridded digital elevation models (DEMs) for use in numerical studies. Here we present a simple yet statistically sound method for processing ice-penetrating radar data and describe a technique for interpolating these data onto a regular grid. DEMs generated for Trapridge Glacier, Yukon Territory, Canada, are used to derive geometric quantities that give preliminary insights into the underlying basin-scale hydrological system. This simple geometric analysis suggests that at low water pressures a dendritic drainage network exists that evolves into a uniaxial morphology as water pressure approaches flotation. These predictions are compared to hydraulic connection probabilities based on borehole drilling.

Quadrature solutions for temperature fields in the central regions of ice sheets where flow is by internal deformation are presented and compared with numerical computations and the European Ice Sheet Modelling Initiative (EISMINT) Benchmark standards. The solution appears to be more accurate than the EISMINT Benchmark standards for thermomechanically uncoupled temperature distributions. The ability of finite-difference and pseudo-spectral methods to reproduce this solution is considered.

Errors in the EISMINT Benchmark are larger than expected. The possibility that they could arise from inaccurate evaluation of the vertical velocity is considered, and found to be unlikely. Formulae for computing the vertical velocity in Jenssen's σ coordinate are complicated and may lead to programming errors. A simple form of the heat-transport equation in the σ coordinate for flow by sliding and internal deformation is derived arid presented, and it is shown how this form is particularly suitable for finite-element computations. It is also shown that only two quadratures are necessary to compute the vertical and horizontal heat transport and the ice flux.

The rate of margin migration of an ice stream can be determined using repeat measurements of the surface velocity profile within the marginal shear zone. The method relies on the assumption that the velocity profile is a characteristic feature of the margin, and that this profile is shifted laterally as the margin migrates. Application of the method to Ice Stream B, Antarctica, indicates that the southern margin is moving outward into the inland ice at a rate of at least 9.7 ± 1.1 m a−1.

The kinematic boundary condition al the glacier surface can be used to give glacier mass balance at a point as a function of changes in the surface elevation, and of the horizontal and vertical velocities. Vertical velocity can in turn be estimated from basal slope, basal ice velocity and surface strain. In a pilot study on the tongue of Griesgletscher, Swiss Alps, the applicability of the relation for modelling area-wide ice flow and mass-balance distribution is tested. The key input of the calculations, i.e. the area-wide surface velocity field, is obtained using a newly developed photogrammetric technique. Ice thickness is derived from radar-echo soundings. Error estimates and comparisons with stake measurements show an average accuracy of approximately ±0.3 ma-1 for the calculated vertical ice velocity at the surface and ±0.7 ma-1 for the calculated mass balance. Due to photogrammetric restrictions and model-inherent sensitivities the applied model appeared to be most suitable for determining area-wide mass balance and ice flow on flat-lying ablation areas, but is so far not very well suited for steep ablation areas and most parts of accumulation areas. Nevertheless, the study on Griesgletscher opens a new and promising perspective for the monitoring of spatial and temporal glacier mass-balance variations.

Synthetic aperture radar- (SAR-)derived ice-motion vectors and SAR interferometry were used to study the sea-ice conditions in the region between the coast and 75° N (~ 560 km) in the East Siberian Sea in the vicinity of the Kolyma River. ERS-1 SAR data were acquired between 24 December 1993 and 30 March 1994 during the 3 day repeat Ice Phase of the satellite. The time series of the ice-motion vector fields revealed rapid (3 day) changes in the direction and displacement of the pack ice. Longer-term (≥ 1 month) trends also emerged which were related to changes in large-scale atmospheric circulation. On the basis of this time series, three sea-ice zones were identified: the near-shore, stationary-ice zone; a transitional-ice zone;and the pack-ice zone. Three 3 day interval and one 9 day interval interferometric sets (amplitude, correlation and phase diagrams) were generated for the end of December, the begining of February and mid-March. They revealed that the stationary-ice zone adjacent to the coast is in constant motion, primarily by lateral displacement, bending, tilting and rotation induced by atmospheric/oceanic forcing. The interferogram patterns change through time as the sea ice becomes thicker and a network of cracks becomes established in the ice cover. It was found that the major features in the interferograms were spatially correlated with sea-ice deformation features (cracks and ridges) and major discontinuities in ice thickness.

A detailed investigation was performed of ice-dynamic conditions across the southern grounding zone of Ekströmisen, at the mouth of a regular East Antarctic outlet glacier not characterized by ice-streaming. Accurate field measurements along a profile 20 km long served as input in a two-dimensional numerical ice-flow model in order to calculate the variation of stress, strain rate and velocity with depth. The model results point to a sharp transition between the mechanics of grounded and floating ice, with a transition zone of only a few km between the two. Vertical shear, most of it near the base, was found to be the dominant flow mechanism in grounded ice. This part is also characterized by a distinct succession of surface undulations which are in turn controlled by variations in resistive stresses at the bottom. The associated phase shift between driving stress and basal drag was found to be accommodated by differential longitudinal pushes and pulls at the base. The flow in the upper half of the profile, on the other hand, is extensive everywhere. The adjacent ice shelf is characterized by small stress and velocity gradients in both the vertical and horizontal directions of the vertical section, and therefore has little deformation. A derived longitudinal deviatoric normal stress of only one-tenth of the value required for freely floating ice shelves reflects a large back-stress originating from friction along the side-walls of the narrow embayment.