Mapping of the structural glaciology of Kongsvegen, Svalbard, reveals evidence for four main deformational structures. These are stratification, longitudinal foliation, thrusts and crevasse traces. These structures are considered in terms of their contribution to debris entrainment, transport and subsequent landform development. Stratification is associated with small amounts of supraglacial debris that has been folded with flow-parallel axes; longitudinal foliation in places incorporates basal glacial sediments along folds with flow-parallel axes; and thrusts transport basal debris to the glacier surface. Crevasse traces are not significant in terms of debris entrainment. The entrainment of basal debris along longitudinal foliation is not a universally recognised process. At Kongsvegen this process is attributed to the development of a transposition foliation, in combination with incorporation of debris-rich basal ice or soft basal sediment in the fold complex. Mapping of the landforms in the proglacial area shows that debris incorporated along longitudinal foliation is released as “foliation-parallel ridges” and that transverse ridges mark debris-bearing thrusts. The role of longitudinal foliation in landform development has never been documented in this manner. Although the preservation potential of such ridges may be limited, recognition of foliation-parallel ridges in the Pleistocene landform record has important implications for the interpretation of the dynamics of former ire masses.

Airborne surface elevation profiles of the Harding Ice field, south-central Alaska, were made in 1991 and 1996. Thirteen glaciers were profiled, along with the tipper region of the icefield. The profiles were compared to U.S. Geological Survey topographie maps made in the 1950s, to obtain elevation and volume changes. Comparison of the changes for the different glaciers shows no significant correlation between volume change and the type of glacier or characteristics such as location, aspect, size, slope or terminus changes. Estimated total volume change tor this ~43 year period is about -34 km3, which corresponds to an area-average elevation change of-21 m. The estimated error in this elevation change of 5 m is mainly due to errors in the maps at higher elevations. Our measurements provide an accurate baseline against which future determinations of volume change can be made.

Digital airborne radar data were collected during the 1987-88 Antarctic field season in nine gridded blocks covering the downstream portions of Ice Stream B (6km spacing) and Ice Stream C (11 km spacing), together with a portion of ridge BC between them. An automated processing procedure was used for picking onset times of the reflected radar pulses, converting travel times to distances, interpolating missing data, converting pressure transducer readings, correcting navigational drift, performing crossover analysis, and zeroing rémanent crossover errors. Interpolation between flight-lines was carried out using the minimum curvature method.

Maps of ice thickness (estimated accuracy 20 m) and basal-reflection strength (estimated accuracy 1 dB) were produced. The ice-thickness map confirms the characteristics of previous reconnaissance maps and reveals no new features. The reflection-strength map shows pronounced contrasts between the ice streams and ridge BC and between the two ice streams themselves. We interpret the reflection strengths to mean that the bed of Ice Stream C, as well as that of Ice Stream B, is unfrozen, that the bed of ridge BC is frozen and that the boundary between the frozen bed of ridge BC and the unfrozen bed of Ice Stream C lies precisely below the former shear margin of the ice stream.

Based on glacier mass-balance records collected, compiled and published by the World Glacier Monitoring Service (WGMS), an attempt is made here to determine whether there were abrupt decadal-scale changes in annual balance on glaciers in the Tien Shan Mountains. There appears to have been an abrupt change in the mid-1970s; two-thirds of the change was caused by an increase in annual ablation and the other one-third by a decrease in the annual accumulation. Thus, the change was caused mainly by a summer-temperature increase, supplemented by a decrease in annual snowfall. Climate records suggest that this event, in effect, is in response to a shift hi the background state of the coupled Ocean -atmosphere system over the Tropical Pacific during the mid-1970s.

Ablation climate studies were made at two locations in northern Greenland in the summers of 1993 and 1994, respectively. Daily ablation was measured at ten stakes within a small area, and the data were compared with each other to detect gross errors. For example, high standard deviations for data taken on the same day, or low correlations between data series at different stakes, indicate erroneous data. After discarding data for one stake in 1993 and two stakes in 1994, random errors in daily ablation data for individual stakes are ± 5 kg m 2 d−1, which is further reduced to only about db 2 kg m−2d−1 by averaging over eight or nine stakes. Random errors in calculated energy balances using the present ablation data are much lower than found in earlier stuthes in West Greenland where ablation was only measured on three stakes without any attempt to detect gross errors. Aside from day-to-day errors, there are ±10% differences in mean ablation at different stakes, which are probably caused by small-scale variations in surface albedo. Such interstate differences give a ± 10% uncertainty in positive degree-day factors, which are 9.8 ± 0.9 and 5.9 ± 0.6 kg m 2 d −1 deg −1 for the two sites.

Will worldwide sea level soon rise rapidly because of a shrinkage of the West Antarctic ice sheet (WAIS)? Here I give a personal perspective of that probability. The crucial question is not whether large changes in ice mass can occur, but how likely it is that a large, rapid change, say a several-fold increase in the 20th-century rate of about 2 mm a-1, will occur in the next century or two from a West Antarctic cause.

Twenty years ago Weertman proposed that a marine ice sheet is inherently unstable. But Weertman’s analysis was based on a simple model of a marine ice sheet that did not include fast-flowing, wet-based ice streams, which are now known to dominate the grounded ice sheet. Modern analyses do not definitively determine just how ice streams affect the stability of the WAIS, but it can at least be said that there is no compelling theoretical reason to expect a rapid rise in sea level from the WAIS triggered by ice-shelf thinning.

Of the three main ice-drainage systems in the WAIS, the one that flows into Pine Island Bay might be a particularly likely site for accelerated flow since there is no ice shelf to restrain the inflowing ice streams, yet measurements show that this system is not significantly out of mass balance. If the “Ross Embayment” system, which has undergone several sudden glacial reorganizations in the last thousand years, were unstable one might expect a history of large changes in the total outflow of ice into the Ross Ice Shelf, yet the total outflow in the “Ross Embayment” has remained relatively unchanged despite the large internal perturbations, a fact that, points to a stable, not an unstable, system. Study of the third major drainage from the WAIS, into the Ronne Ice Shelf, also suggests that there is no gross discordance between the present velocity vectors and flow tracers in the ice shelf, although the evidence is limited.

In the light of the evidence for recent stability, it is difficult to see how climate warming (whether anthropogenic or natural) could trigger a collapse of the WAIS in the next century or two. Thus, I believe that a rapid rise in sea level in the next century or two from a West Antarctic cause could only occur if a natural (not induced) collapse of the WAIS were imminent. Based on a concept of pseudo-random collapse once per major glacial cycle, I estimate the chances of that to be on the order of one in a thousand.

During the 1991 94 summer held seasons, time-correlated measure-merits of water pressure and surface speed were made over and down-glacier from a major riegel on Storglaciären, Sweden. Measurements were made at sub-hourly time-scales in order to discern details in the diurnal cycle. Large water-input events, typically associated with rain storms, produced coherent, lagged surface-velocity responses that could be understood in terms of till deformation or decoupling, and these have been discussed elsewhere. The consequences of smaller diurnal water-pressure events w ere more enigmatic, in that acceleration of ice flow generally preceded the onset of the local water-pressure rise. From consideration of these data and other work done on the hydrology of Storglaciären, we infer that the ice in this area is generally pushed from behind via a relaxation in extensional strain across the riegel. Hence, accelerations occur in response to increases in water pressure that occur up-glacier and that precede local water-pressure rises. In addition, following a period of large storm events, surface speeds became more spatially coherent and were in phase with the diurnal water-pressure cycle. This suggests that the large water-pressure events lead to a spatially more homogeneous subglacial drainage system. Sliding laws need to take into account such temporal changes in spatial coherence of the subglacial drainage system.

This paper addresses the coupling of flows of ice, till and water, and the issue of whether such coupling provides mechanisms for meso-scale (kilometres to tens of kilometres) variability in ice-sheet flow and texture. The question of whether effective pressures at the ice-bed interface are statically or hydraulically controlled is examined in this paper. The answer is scale dependent, and has a significant effect on the relationship between ice surface and basal topography.

The consequences of these considerations on till flow, coupled ice–till flow and coupled ice, till and water flow are examined. An analysis of till-flow kinematics and shock formation is carried out. The linear stability of coupled long-wavelength ice-till flow is analysed, and regions in parameter space where this flow is unstable, with rather small rate constants are found. Upstream-moving ice surface waves are predicted. The linear stablity of coupled ice–till–water flow is examined, where water flow is modelled using a basal flow system with effective-pressure-dependent properties. Again, regions in parameter space where the system is linearly unstable are found, this time with relatively rapid rate constants. The water pressure exhibits “breather” modes.

These analyses assume that there is a substantial basal traction. A problem with models of ice streams wholly restrained at the side is identified: they seem to predict erosion rates which are unfeasibly large.

There appears to be sufficient variability in the ice–till– water system to potentially explain texture in ice-stream surfaces, variations in ice-stream thickness of tens of metres not directly relatable to topography, and waves moving upstream or downstream. Most importantly, the ice-stream–bed system is shown to exhibit meso-scale variability simply by coupling ice flow according to the shallow-ice approximation, till flow according to the hydrostatic thin-till approximation and water flow according to an effective-pressure-dependent hydraulics.

High-resolution maps of the glacier bed are developed through a pseudo-three-dimensional migration of a dense array of radio-echo sounding profiles. Resolution of three-dimensional maps of subglacial surfaces is determined by the radio-echo sounding wavelength, data sparing in the field, and migration. Based on synthetic radio-echo sounding profile experiments, the maximum resolution of the final map cannot exceed one half-wavelength. A methodology of field and processing techniques is outlined to develop a maximum-resolution map of the glacier bed. The field and processing techniques are used to develop a map of the glacier bed below part of Worthington Glacier, a temperate valley glacier in south-central Alaska. The field techniques and the processing steps used on the glacier result in a map of 20 m x 20 m resolution.

A three-dimensional, finite-difference model based on a first-order solution of the ice-flow equations is applied to Haut Glacier d’Arolla, Switzerland. The numerical model successfully converges at horizontal resolutions down to 70 m, so a number of detailed comparisons with field data can be made. Modelled surface velocities with no basal sliding component are compared with surface velocities observed on the glacier over four different time periods. The best fit is achieved with over-winter surface velocities (R2 = 0.75) using a rate factor, A, in Glen’s flow law of 0.063 a−1 bar−3. Surface zones of maximum computed effective stress display a high level of coincidence with observed crevassing, the orientation of which is successfully predicted by the direction of the tensile component of the computed principal surface stress. Comparison of the relative magnitude and direction of computed principal stresses with principal strains measured at the ice surface also correspond closely. In an attempt to simulate the observed annual velocity distribution within a cross-section of the glacier tongue, we incorporate two basal-motion patterns into the model. By treating net annual ice motion as a time-weighted composite of three separate flow situations: normal sliding, enhanced sliding and no sliding, we are able to reproduce the key features of the observed cross-sectional ice and basal slip velocity distributions. These experiments indicate there may be substantial decoupling taking place along an elongated narrow zone at the bed of Haut Glacier d’Arolla and that this decoupling interacts in a complex manner with the englacial stress and strain field.

Ice temperature was measured in and around the chaotically crevassed south margin of Ice Stream B, Antarctica, from 1992 to 1994. The temperatures at 30 m depth in the chaotic zone are about 12 K lower than in the adjacent uncrevassed ice, due to the ponding of cold winter air. At depths greater than 150 m, there is clear evidence of internal heating of the ice due to the large shear déformation rate in the marginal zone. Analysis of the depth of cooling below the crevasses and of the internal heating gives two pieces of information. First, over the last half century the lateral shear stress averaged 2.0 x 105 Pa in the top third of the margin and, second, the margin moved outward at an average rate of 7.3 m a−1. These values do not involve any assumptions about the How law of ice. The uncertainties are roughly 20%. The value of lateral shear stress indicates that the most of the drag on the ice stream is along its sides.

Using an active transponder with the ERS-I and ERS-2 radar altimeters, the distance to the satellite was measured at a location close to the GRIP site, Greenland, at an altitude of 3.2 km. The measurement was executed while the transponder was in the “ice-tracking mode”. It includes a bias due to the propagation delay. The location of the transponder was determined using the global positioning system.

The transponder signal was modeled and the distance from the altimeter to the effective reflection point of the transponder was determined. Since the transponder was located within 1 km of the ground tracks, the measurement was corrected for this offset. A correction was also done for the surface slope, resulting in the distance (plus bias) to the closest sub-satellite point on the surface of the (compact) snow.

The transponder signal was then removed from the radar altimeter waveform, enabling the determination of the distance (plus bias from the altimeter to the first reflective surface within the snow. The différence between this distance and that obtained using the transponder was < 2 m. This shows that the surface which gives rise to the first return of the reflection agrees with the surface of the (compact, dry) snow at this high-altitude location. This is an important result to be used when studying ice-cap topography using satellite radar altimetry.

Dielectric profiling (DEP), a capacitive method for measuring an ice core’s electrical conductivity and permittivity, has been in use for several years. We have built an improved DEP instrument with a modified guarded-electrode geometry whose capacitance can be determined from potential theory. The new instrument can measure an ice core’s conductivity and permittivity to within a few per cent error, a high enough accuracy to infer its density using well-known models.

Under idealized conditions, when pressurized water has access to all low-pressure areas at the glacier bed, a sliding instability exists at a critical pressure, pc, well below the overburden pressure, p0. The critical pressure is given by , where l is the wave length and a is the amplitude of a sinusoidal bedrock, and T is the basal shear stress. When the subglacial water pressure, pw, approaches this critical value, the area of ice-bed contact, △l, becomes very small and the pressure on the contact area becomes very large. This pressure is calculated from a force balance and the corresponding rate of compression is obtained using Glen’s flow law for ice. On the assumption that compression in the vicinity of the contact area occurs over a distance of the order of the size of this area, Δl, a deformational velocity is estimated. The resultant sliding velocity shows the expected instability at the critical water pressure. The dependency on other parameters, such as wavelength l and roughness a/l, was found to be the same as for sliding without bed separation.

A ring-shear device was used to study the factors that control the ultimate(steady) strength of till at high shear strains.Tests at a steady strain rate and at different stresses normal to the shearing direction yielded ultimate friction angles of 26.3° and 18.6° for tills containing 4% and 30% clay-sized particles, respectively Other tests at steady normal stresses and variable shear-strain rates indicated a tendency for both tills to weaken slightly with increasing strain rate. This weakening may be due to small increases in till porosity.

These results provide no evidence of viscous behavior and suggest that a Coulomb-plastic idealization is reasonable for till deformation. However, viscous behavior has often been suggested on the basis of distributed shear strain observed in subglacial till. We hypothesize that deformation may become distributed in till that is deformed cyclically in response to fluctuations in basal water pressure. During a deformation event, transient dilation of discrete shear zones should cause a reduction in internal pore-water pressure that should strengthen these zones relative to the surrounding till, a process called dilatant hardening. Consequent changes in shear-zone position, when integrated over time, may yield the observed distributed strain.

A derivative of the Danish tipping-lower intermediate-depth ice-coring drill system has been developed for use in remote areas. The system as described is capable of drilling to depths of 250 m, although the drill has been proven to 350 m. The drill system is designed to be transported by hand and consumes 700 W maximum power from a 24 V Ac. battery bank charged by solar panels or a small generator. A field test in the St Elias Mountains, Yukon, produced good-quality core to 161.6 m depth; drilling at 6500 m a.s.l. on Far East Rongbuk Glacier, Mount Everest, produced good core to 42 m depth.

A one-dimensional time-dependent flowline model of Rhonegletscher, Switzerland, has been used to test the glacier’s response to climatic warming. Mass-balance variations over the last 100 years are obtained from observations of the equilibrium-line altitude (ELA) and a reconstruction of the ELA based on a statistical correlation between temperature and ELA. For the period prior to AD 1882, for which no reliable climate data exist, we chose equilibrium-line altitudes that enabled us to simulate accurately the glacier length from AD 1602.

The model simulates the historical glacier length almost perfectly and glacier geometry very well. It underestimates glacier-surface velocities by 1-18%. Following these reference experiments, we investigated the response of Rhonegletscher to a number of climate-change scenarios for the period AD 1990-2100. For a constant climate equal to the 1961-90 mean, the model predicts a 6% decrease in glacier volume by AD 2100. Rhonegletscher will retreat by almost 1 km over the next 100 years at this scenario. At a warming rate of 0.04 K a-1, only 4% of the glacier volume will be left by AD 2100.

A new portable device for ice-core drilling, specially designed for thin tephra-layer sampling, was tested on the South Shetland glaciers during the 1994-95 Antarctic summer. The machine is based on a combination of the standard paleomagnetism drilling machine and a specially built drill-bit designed for ice-coring.

The non-linearity of the ice-flow law or a local accumulation low over an ice divide can cause isochrones (internal layers) to be shallower under the divide relative to the flanks, forming a “divide bump” in the internal layer pattern. This divide signature is analyzed using ice-flow models and inverse techniques to detect and quantify motion of the Siple Dome ice divide, West Antarctica. The principal feature indicating that migration has occurred is a distinct tilt of the axis of the peaks of the warped internal layers beneath the divide. The calculated migration rate is 0.05-0.50 m a−1 toward Ice Stream D and depends slightly on whether the divide bump is caused by the non-linearity of ice flow or by a local accumulation low. Our calculations also suggest a strong south-north accumulation gradient of 5-10 x 10−6 a−1 in a narrow zone north of the divide. A consequence of divide migration is that pre-Holocene ice is thickest about 0.5 km south of the present divide position. Divide motion indicates that non-steady processes, possibly associated with activity of the bounding ice streams, are affecting the geometry of Siple Dome. The migration rate is sufficiently slow that the divide bump is maintained in the internal layer pattern at all observable depths. This suggests that major asynchronous changes in the elevation or position of the bounding ice streams are unlikely over at least the past 103-104 years.

Kongsvegen is a 102 km2 sub-polar (polythermal) surge-type glacier in northwest Spitsbergen, Svalbard. It surged just before 1948 and is currently in its quiescent phase. Measurements of surface geometry since 1966 show a retreat of the front and strong thinning of up to 75 m in the ablation area, and a build-up of up to 32 m in the accumulation area. Present-day annual velocities along the glacier are low, from 1.4 up to 3.6 m a-1. The measured mean net balance for the period 1987-94 and the balance reconstructed back to 1967 show a weak positive balance of about 0.1 m w.e. The measured actual ice flux is low and the mass transfer down-glacier at the ELA is only about 3-20% of that required for steady state. Thus, the glacier is building up towards a new surge. The total thickening rate on Kongsvegen is somewhat higher than in other cases from Svalbard, but it is small compared with other well-studied surge-type glaciers in Alaska and the Pamirs. This relatively low rate of change is a function of the low accumulation rate and the relatively cold climate compared to other areas and is common for surge-type glaciers in Svalbard.