Patterns of strain rate and slope on the ice streams are unusual. They cannot be accounted for in the usual way as due to standing waves in ice flow over a basal obstruction to flow (such as a sticky spot). The features are studied using the force-budget technique. The conventional flow law is used, together with measurements of surface strain rate and shape of the glacier, to compute basal drag. The results for Ice Stream C are as expected, in that the drag varies from site to site but is directed inland, restraining the flow. The calculated drag at the base of Ice Stream B, on the other hand, is in places such that it acts to propel the glacier forward. This result is untenable. Either the conventional flow law is not applicable to Ice Stream Β or there are large spatial variations in ice stiffness, perhaps associated with foliation, or both.

The basal shear stress of an ice stream may be supported disproportionately on localized regions or “sticky spots”. The drag induced by large bedrock bumps sticking into the base of an ice stream is the most likely cause of sticky spots. Discontinuity of lubricating till can cause sticky spots, but they will collect lubricating water and therefore are unlikely to support a shear stress of more than a few tenths of a bar unless they contain abundant large bumps. Raised regions on the ice-air surface can also cause moderate increases in the shear stress supported on the bed beneath. Surveys of large-scale bed roughness would identify sticky spots caused by bedrock bumps, water-pressure measurements in regions of thin or zero till might reveal whether they were sticky spots, and strain grids across the margins of ice-surface highs would show whether the highs were causing sticky spots. Sticky spots probably are not dominant in controlling Ice Stream Β near the Upstream Β camp, West Antarctica.

Micro-earthquakes have been monitored at two locations on Ice Stream Β and one on Ice Stream C using a seismographic array built specifically for that purpose. Subglacial micro-earthquakes arc 20 times more abundant beneath Ice Stream C than beneath Ice Stream B, despite the 100 times more rapid movement of Ice Stream B. Triangulation shows the foci beneath Ice Stream C, like those beneath Ice Stream B, to be within a few meters of the base of the ice, presumably within the uppermost part of the bed, and fault-plane analysis indicates slips on horizontal planes at about a 30° angle to the presumed direction of formerly active flow. Source parameters, computed from spectra of the arrivals, confirmed that the speed of slip is three orders of magnitude faster beneath Ice Stream C than beneath Ice Stream Β which means that a five orders-of-magnitude greater fraction of the velocity of Ice Stream C is contributed by the faulting, although that fraction is still small. We attribute the difference in activity beneath the two ice streams to the loss of dilatancy in the till beneath Ice Stream C in the process that led to its stagnation.

The speed of Ice Stream B, Antarctica, was measured twice a day-over a 1 month study period, and found to be steady at about the ±3½% level, the sensitivity of the measurements. The vertical strain was measured at three sites over a 1 year period at 1 h intervals with sensitivities of 2 or 0.2 ppm. The strain rate varied on all time-scales. Events of high strain rate were observed, but never at more than one site at a time. They can probably be understood in terms of local modification of the strain field associated with crevassing. Diurnal variation in strain rate was observed at one and possibly two sites during two summers. The seismicity was measured at all three sites, and diurnal and seasonal variations were prominent at all, the seismicity being much more intense in winter. Several possible causes of the diurnal variations in strain and seismicity are considered: thermal and atmospheric effects, and the effects of tides in the Ross Sea.

Field data of ice thickness, velocity and accumulation of the ice plain region of Ice Stream B, Antarctica, are organized into a geographic information system (GIS) and used to calculate the spatial pattern of net mass balance. Overall, the ice plain is thickening at a rate of 0.13 ± 0.05 m a−1. Large uncertainties of the high-resolution calculations are reduced by spatial averaging, revealing a number of areas significantly out of balance. Ice in the broad diverging flow field is mostly thinning but a thinner region of the ice plain which includes ice raft “a” is thickening. Thickening is also indicated farther downstream on the ice plain, matching the thickening already calculated to be occurring upstream of Crary Ice Rise. The patterns of net mass balance for the two major tributaries of Ice Stream Β are asymmetric, demonstrating a real difference in the current dynamics of these two tributaries.

Measurements of velocity have been made on and next to Ice Streams Β and C, West Antarctica. The results are more precise than previous work and constitute a 93% increase in the number of values. These velocities are used to describe the confluence of flow into the ice streams and the development of fast ice-stream flow. The onset of fast-streaming flow occurs in many separate tributaries that coalesce down-glacier into the major ice streams. For those inter-stream ridges that have been studied, the flow is consistent with steady state. Along Ice Stream B, gradients in longitudinal stress offer little resistance to the ice flow. The transition from basal-drag control to ice-shelf flow is achieved through reduced drag at the glacier base and increased resistance associated with lateral drag. Velocities in the trunk of Ice Stream C are nearly zero but those at the up-glacial head are similar to those at the head of Ice Stream B.

By studying diatoms recovered from sediment taken from beneath the West Antarctic ice sheet (Ice Stream B), Scherer (1991) concluded that this feature must have collapsed at least once during the past 600 kyear and that the two likely candidates for time of collapse were oxygen-isotope stage 11 (362–423 kyear BP) and sub-stage 5e (110–128kyear BP). This conclusion does not stand up to critical examination of the data, however. Specifically, the diatom datum level, the last occurrence of Actinocyclus ingens, used to constrain the 600 kyear date is diachronous into higher latitudes and does not apply in sediments recovered from near Antarctica. Secondly, the additional diatoms used to constrain the time of collapse to the late Quaternary either range before the late Quaternary or have no published geological record. In spite of this, there are data to suggest that one or more late Quaternary interglacial intervals were as warm as, or warmer than, the present. As yet, however, no direct evidence exists to incorporate a West Antarctic ice sheet collapse into these scenarios. It is suggested that this debate is best resolved by study of deep-sea sediments of late Quaternary age recovered from around the Antarctic continent.

Digital airborne-radar data were collected during the 1988–89 Antarctic field season in six gridded blocks covering the upstream parts of Ice Streams A, Β and C. An automated processing procedure was developed for picking onset times, converting travel times, interpolating missing data, converting pressure-transducer readings, correcting navigational drift, performing cross-over analysis and zeroing remanent cross-over errors. Cross-over analysis was used to remove the effects of temporal variations in atmospheric pressure and to estimate errors. Interpolation between flight lines was carried out using the Kriging method. Surface elevation was referred to the Rapp Set A geoid by tying the gridded surface to satellite-surveyed ground stations, using a planar-model fit.

Maps of surface elevation, ice thickness and bottom topography with standard-error estimates of 4–9 m for surface elevation and 30–60 m for ice thickness and bottom topography were produced. These maps show that the locations of the ice streams are not clearly reflected in either the surface or basal topography, so are presumably controled by basal or internal conditions, that there is no clearly demarcated transition zone between sheet flow and streaming flow, that there is no clear cut evidence for the capture of the catchment of Ice Stream C by Ice Stream B, but that Ice Stream Β does drain virtually the entire region between the lateral boundaries of Ice Streams A and C.

The study of the phase of reflections of Ρ waves off the base of Ice Streams Β and C, and Ridge BC, indicates that acoustic impedances of the beds of both ice streams vary laterally. In some places, the impedance is higher than in the ice (a high-impedance bed) and in some places it is less (a low-impedance bed). The estimated impedances in a dilated bed (porosity 0.4) and in a model of the lowermost ice that takes into account the relatively low P-wave speed in ice at or very near the melting point are nearly the same. Whether the impedance in the bed is greater or less than in the ice could depend on minor changes in the nature of the sediments composing the bed, or the physical state of the bed (e.g. porosity) that could occur laterally. Lateral variations of this kind provide a ready explanation for the observations on Ice Stream B. The bed under a substantial part of Ice Stream C that exhibits a low-impedance bed also must be dilated. The evaluation of the state of the bed under the rest of Ice Stream C and on Ridge BC requires further analysis, which is in progress.

Numerous features of glaciological significance appear on two adjoining SPOT High Resolution Visible (HRV) images that cover the onset region of Ice Stream B, Antarctica. Many small-scale features, such as crevasses and drift plumes, have been observed in aerial photography. Subtle large-scale features, such as long flow traces that have not been mapped previously, are clear in the satellite imagery. Newly discovered features include ladder-like runners and rungs within certain shear margins, flow traces that are parallel to ice flow, unusual crevasse patterns and flow traces originating within shear margins.

Low-frequency surface-based radar-profiling experiments on Ice Streams Β and C, West Antarctica, have yielded high-resolution images which depict folding of the internal layers that can aid in the interpretation of ice-stream dynamics. Unlike folding seen in most earlier radar studies of ice sheets, the present structures have no relationship to bedrock topography and show tilting of their axial fold planes in the flow direction. Rather than being standing waves created by topography or local variations in basal shear stress, the data show that these folds originate upstream of the region of streaming flow and are advected into the ice streams. The mechanism for producing folds is hypothesized to be changes in the basal boundary conditions as the ice makes the transition from inland ice to ice-stream flow. Migration of this transition zone headward can then cause folds in the internal layering to be propagated down the ice streams.

Satellite imagery is used as a basis to review and critique the results of studies at the mouths of Ice Streams Β and C and Crary Ice Rise. In many cases, these past analyses are extended by taking advantage of the broad coverage within each image. New perspectives are provided by the image data and some longstanding controversies are resolved. The grounding line is easily delineated and mapped in areas covered by imagery. Extensive areas of grounded ice with complex patterns of flow stripes are identified on the flanks of Crary Ice Rise. The imagery also allows a corrected map of surface topography in the vicinity of the Downstream Β camp. New questions are posed by hitherto unseen features. Data from the IGY traverse of the Ross Ice Shelf in 1957 are included to demonstrate that large changes have occurred in the past almost 30 years in the area upstream of Crary Ice Rise. These changes include modifications in the surface topography, elimination of crevasses and increases in the ice thickness by approximately 60 m.

Five short-pulse radar profiles were run across the edge of inactive Ice Stream C, one of the “Ross” ice streams that flows from the West Antarctic inland ice sheet into the Ross Ice Shelf. Scatter from buried crevasses, which we presume were at the surface of the ice stream when it was active, creates hyperbolae on the radar records. A density-depth curve and local accumulation rates were used to convert the picked travel times of the apices of the hyperbolae into stagnation ages for the ice stream. Stagnation ages are 130 ± 25 year for the three profiles farthest downstream and marginally less (100 ± 30 year) for the fourth. The profile farthest upstream shows a stagnation age of only ~30 year. We believe that these results indicate a “wave” of stagnation propagating at a diminishing speed upstream from the mouth of the ice stream, and we suggest that the stagnation process involves a drop in water pressure at the bed due to a conversion from sheet flow to channelized water flow.

Repeat aerial photography is used to obtain closely spaced measurements of velocity and elevation over a complete transect of Ice Stream tributary B2, including the shear margins, the fast ice of the ice stream and several unusual features, as well as the UpB camp. Persistent features, mainly crevasses, are tracked to provide 1541 values of velocity and 1933 values of elevation. These are used to describe ice flow in the ice stream. Within the ice stream, the dominant velocity gradient is lateral shear. Crevasse patterns are studied in relation to measured velocity gradients. Crevasses intersect one another at acute angles, indicating that their origin is deeper than the depth to which crevasses penetrate. One feature within the ice stream seems to be a raft of stiff ice. Others are crevasse trains. Also, there are spreading ridges, perhaps due to upwelling ice. There is no evidence of large sticky spots within the studied transect, i.e. no steep surface slopes with associated surface stretching just up-glacier and surface compression down-glacier.

Solute acquisition by Alpine glacial meltwaters is the result of the coupling of different pairs of reactions, one of which usually involves dissolved gases. Hence, the availability of atmospheric gases to solution is an important control on the composition of glacial meltwaters. The chemical compositions of the two main components of the bulk meltwater, quick flow and delayed flow, are dominated by different geochemical processes. Delayed flow waters are solute-rich and exhibit high p(CO2) characteristics. The slow transit of these waters through a distributed drainage system and the predominance of relatively rapid reactions, such as sulphide oxidation and carbonate dissolution, in this environment maximize solute acquisition. Quick-flow waters are dilute, both because of their rapid transit through ice-walled conduits and open channels, and because the weathering reactions are fuelled by relatively slow gaseous diffusion of (CO2) into solution, despite solute acquisition being dominated by rapid surface exchange reactions. As a consequence, quick flow usually bears a low or open-system p(CO2) signature. Bulk meltwaters are more likely to exhibit low p(CO2) values when suspended-sediment concentrations are high, which promotes post-mixing reactions. This conceptual model suggests that the composition of both quick flow and delayed flow is likely to be temporally variable, since kinetic, rather than equilibrium, factors determine the composition.

The pattern of elevation change along the length of glaciers caused by retreat from Neoglacial maximum to present is investigated and described in terms of a profile-shape factor f defined as the ratio of the average thickness change to the local thickness change at the present terminus. The factor f is relevant to estimation of the change in volume of a glacier in relation to its change in length and the response time to a change in climate. Thickness-change profiles on a selection of 15 mountain valley glaciers give f in the range 0.1–0.4 with an average of 0.28. The elevation changes are in all cases largest near the terminus and decrease headward more rapidly than a linear variation with distance that would correspond to f of 0.50, a value assumed to be characteristic in the past. The smaller values of f indicate reduced estimates of changes in volume versus change in lengths and smaller response times for glaciers. Total volume change for glaciers and ice caps (excluding ice sheets in Greenland and Antarctica) as estimated roughly by these means, is consistent with previous estimates from mass-balance measurements and supports the view that melting of these ice masses has made a significant contribution to sea-level change over the last century. The observed values of f do not explain an apparent discrepancy between alternative theoretical methods for estimating the response time of glaciers based on length and terminus velocity or based on the thickness and terminus-ablation rate. The discrepancy may be caused in part by the relationship between terminus velocity and ablation rate associated with strongly non-steady-state conditions during advance or retreat. The patterns of elevation change measured for short time intervals of a few decades or less correspond to a wide spread of profile-shape factor f from less than 0.1 to more than 1.2 indicative of transients in the short time-scale response.

In order to estimate the accuracy of altimetric height measurements over ice sheets, an altimeter wave-form simulator has been developed, and different tracking methods have been tested. A large range of surface features, including large-scale and medium-scale features and micro-roughness have been taken into account for modeling of either surface-or volume-scattering. A large set of parameters affects the trailing edge of the radar wave form, so that re-tracking algorithms based on the detection of its leading edge provide better retrievals of the surface height than those based on the analysis of the whole wave form. A volume component is clearly present in the radar wave forms; its effect on the leading edge depends mostly on the snow grain-size (and therefore on the snow temperature) and on the pointing angle. However, on average, the induced error on the snow-surface height estimation should only be around 25 cm.

The importance of the atmospheric boundary layer for the coupling between the climate and an ice sheet is investigated using a slab model of the atmospheric boundary layer. The model is shown to give reasonable agreement with observations over Antarctica and it is used to look at the effect of different ice-sheet shapes on the boundary layer. The importance of entrainment in bringing heat to the surface is highlighted and is shown to be particularly significant when the ice profile becomes steeper. The model could be used as part of an energy-balance model of snow in order to incorporate the interplay of the boundary layer and ice-sheet shape in the ablation process. The slab model could also be used in a GCM as a parameterization of these sub-grid scale processes which are at present ignored in models on a global scale.

We present the results of tests done to determine the tensile behavior of first-year columnar sea ice over a range of temperatures from −20° to −3°C and strain rates of 10−5 and 10−3s−1. The temperature of a test specimen was dictated by its in-situ location within the sea-ice sheet; samples located near the top of the sea-ice sheet were tested at the lower temperatures. A tensile load was applied along the cylindrical axes of the test specimens, which were perpendicular to the growth direction of the ice. Results showed that the maximum stress reached during a test was most strongly influenced by temperature, while the failure strain and the modulus were principally affected by the loading rate. A model relating the tensile strength of the ice to its porosity based on temperature-dependent variations in the brine-pocket geometry is evaluated.

Grand Pacific and Melbern Glaciers, two of the largest valley glaciers in British Columbia, have decreased over 50% in volume in the last few hundred years (total ice loss = 250–300 km3). Melbern Glacier has thinned 300–600 m and retreated 15 km during this period; about 7 km of this retreat occurred between the mid-1970s and 1987, accompanied by the formation of one of the largest, presently existing, ice-dammed lakes on Earth. Grand Pacific Glacier, which terminates in Tarr Inlet at the British Columbia-Alaska boundary, retreated 24 km between 1879 and 1912. This rapid deglaciation has destabilized adjacent mountain slopes and produced spectacular ice-marginal land forms. The sediments and land forms produced by historic deglaciation in Melbern–Grand Pacific valley are comparable, both in style and scale, to those associated with the decay of the Cordilleran ice sheet at the end of the Pleistocene (c. 14–10 ka BP). Rates of historic and terminal Pleistocene deglaciation also may be comparable.