Pit-wall samples were collected from two sites about 2 km apart on Agassiz Ice Cap, Ellesmere Island, Northwest Territories, Canada, in 1992, 1993 and 1994, and from a site a further 1 km distant, in order to study spatial and seasonal variations in snow chemistry. Two of the pits were dug in wind-scoured zones and one in an unscoured zone. Although a large part of the winter snow is removed from the scoured zones (which do not show very negative δ18O values) the winter/spring anion peaks are still evident; this may be due to the predominance of dry deposition in mid-winter. The Cl− and SO42– ions peak in late winter/early spring, while NO3− peaks both in late winter/early spring and in summer. Vertical concentration profiles of all anions did not significantly alter over a 2 year period, indicating that there are no serious post-depositional changes due to evaporation, snow melting or photochemical reactions. However, comparisons between stake/board snow-accumulation measurements and those derived from the least scoured pit indicate that a single pit will represent annual accumulation rates for a local area only.

Pressure and tracer measurements in boreholes drilled to the bottom of Ice Stream B, West Antarctica, are used to obtain information about the basal water conduit system in which high water pressures are developed.These high pressures presumably make possible the rapid movement of the ice stream. Pressure in the system is indicated by the borehole water level once connection to the conduit system is made. On initial connection, here also called “breakthrough” to the basal water system, the water level drops in a few minutes to an initial depth in the range 96–117 m below the surface. These water levels are near but mostly somewhat deeper than the floation level of about 100 m depth (water level at which basal water pressure and ice overburden pressure are equal), which is calculated from depth-density profiles and is measured in one borehole. The conduit system can be modelled as a continuous or somewhat discontinuous gap between ice and bed; the thickness of the gap δ has to be about 2 mm to account for the water-level drop on breakthrough, and about 4 mm to fit the results of a salt-tracer experiment indicating downstream transport at a speed of 7.5 mm s−1. The above gap-conduit model is, however, ruled out by the way a pressure pulse injected into the basal water system at breakthrough propagates outward from the injection hole, and also by the large hole-to-hole variation in measured basal pressure, which if present in a gap-conduit system with δ = 2 or 4 mm would result in unacceptably large local water fluxes. An alternative model that avoids these objections, called the “gap opening” model, involves opening a gap as injection proceeds: starting with a thin film, the injection of water under pressure lifts the ice mass around the borehole, creating a gap 3 or 4mm wide at the ice/bed interface. Evaluated quantitatively, the gap-opening model accounts for the volume of water that the basal water system accepts on breakthrough, which obviates the gap-conduit model. In order to transport basal meltwater from upstream it is then necessary for the complete hydraulic model to contain also a network of relatively large conduits, of which the most promising type is the “canal” conduit proposed theoretically by Walder and Fowler (1994): flat, low conduits incised into the till, ∼0.1 m deep and perhaps ∼1 m wide, with a flat ice roof. The basal water-pressure data suggest that the canals are spaced ∼50–300 m apart, much closer than R-tunnels would be. The deepest observed water level, 117 m, is the most likely to reflect the actual water pressure in the canals, corresponding to a basal effective pressure of 1.6 bar. In this interpretation, the shallower water levels are affected by loss of hydraulic head in the narrow passageway (s) that connect along the bed from borehole to canal(s). Once a borehole has frozen up and any passageways connecting with canals have become closed, a pressure sensor in contact with the unfrozen till that underlies the ice will measure the pore pressure in the till, given enough time for pressure equilibration. This pressure varies considerably with time, over the equivalent water-level range from 100 to 113 m. Basal pressure sensors 500 m apart report uncorrelated variations, whereas sensors in boreholes 25 m араrt report mostly (but not entirely) well-correlated variations, of unknown origin. In part of the record, remarkable anticorrelated variations are interspersed with positively correlated ones, and there are rare, abrupt excursions to extreme water levels as deep as 125 m and as shallow as 74 m. A diurnal pressure fluctuation, intermittently observed, may possibly be caused by the ocean tide in the Ross Sea. The lack of any observed variation in ice-stream motion, when large percentagewise variations in basal effective pressure were occurring according to our data, suggests that the observed pressure variations are sufficiently local, and so randomly variable from place to place, that they are averaged out in the process by which the basal motion of the ice stream is determined by an integration over a large area of the bed.

The partitioning of resistive force between the bed and sides of Ice Stream B, Antarctica, is obtained for three large areas that have been measured using repeat aerial photogrammetry. Problems associated with data errors and local variations in ice strength and velocity are reduced by considering the areally averaged budget of forces for each photo block. Results indicate that the bed under Ice Stream B must be very weak and unable to provide much resistance. Mechanical control on this ice stream emanates almost entirely from the lateral margins.

Exposure of a full transverse cross-section of the terminus of Fireweed rock glacier, Wrangell Mountains, Alaska, revealed a thin layer of unconsolidated debris mantling a consolidated mélange of ice and rock. The main rock glacier is fed by three tributaries; at the terminus, contacts between the three are sharply defined. Ice content is >50% by volume. Bubble foliation and crystal morphologies of the ice matrix are similar to those reported from glacier ice. Folded ice-rich strata and lenses, foliation planes, and the long-intermediate axial planes of tabular-shaped englacial clasts dip sleeply toward the center line of the rock glacier. The planar structures generally parallel the steep walls of the gorge containing the trunk stream. These steeply dipping, longitudinal structures appear to result from transverse compression where the tributaries converge and the trunk stream narrows down-valley.

Bergschrund-like and moulin-like features at the heads of the middle and west tributaries, respectively, exposed admixed ice and rock similar to that of the terminus but lacking the well-developed flow texture. Talus from the unstable cirque headwall and polygenetic ice both accumulate at the base of the headwall and nourish the tributary rock glaciers.

Seismic reflection data from two lines on Rutford Ice Streem are presented and are compared with data already published from a third line on the ice stream. The amplitude and phase of the reflections have been used io investigate the properties ol the sub-ice material. Multiple reflections on long record-length data allowed calibration of the reflection coefficient at the ice-bed interlace and determination of the acoustic impedance of the bed material. The characteristics of the bed material vary both along and across the ice stream. The average acoustic impedance of the bed material across the glacier at the upstream line is 3.88 x 106 kg−2 s−1. This decreases to 3.19 x 10−6 kg m−2s−1 52 km further downstream. These values are within the rang which is typical of soft sediments. Using acoustic impedance as an indicator of subgiacial porosity, some areas of the ice-stream bed are interpreted as dilatant water-saturated sediments undergoing pervasive deformation. In other areas, the bed is not deforming and basal sliding may be a more important process. The proportional the ice-stream width over which bed deformation occurs increases downstream.

This paper examines the origin of moraine ridges in East Antarctic coastal oases and derives depositional models appropriate for the reconstruction of Quaternary history. On the basis of morphology, structure and sedimentology, four principal types of ridge may be identified: (1) type A moraines from when the basal debris zone crops out near an ice margin; (2) type B moraines form when large recumbent folds develop in the basal debris zone; (3) type C moraines are ice-contact screes and fans which form when debris accumulates at steep or cliffed ice margins; and (4) type D moraines are thrust-block moraines that form when unconsolidated sediment is entrained by freezing, shearing and thrusting of sediment blocks at the base of the glacier. Simple calculations of the rate of debris accumulation at ice margins suggest that type A, B and C moraines take thousands of years to form and record stable ice margins. Type D moraines are structural features that may form relatively quickly when ice margins override unconsolidated sediment. Constructing models to explain the origin of the moraines is an important part of reconstructing the Quaternary history of Antarctic coastal oases, because the models provide a basis for reconstructing the position and behaviour of the ice sheet during advance and retreat

The aim of this paper is to investigate the geophysical characteristics of the Antarctic ice sheet using radar altimetric observations. To do this, we use an altimetric waveform simulator, in situ observations, ERS-1 (European remote-sensing satellite) data and SPOT (Satellite pour l’observation de la terre) images. The small-scale study takes place at Dome C, Terre Adélie, which is a relatively flat region with gentle undulations and low wind speed. Despite this, the altimetric waveform parameters (height, energy, leading edge and trailing edge) are highly noisy. The effect of undulations on the waveform parameters is found to be dominant. The combination of a subsurface signal and a rough surface produces a linear effect on the altimetric backscattering or on the trailing edge of the waveform, but a strongly non-linear effect on the leading edge of the waveform or height estimation. As a consequence, the height measurement is very sensitive to the altimeter technical or orbital characteristics and is not reproducible from one mission to another. Observations show sastrugi fields that enhance the leading edge and affect the whole waveform. Observed local backscattering changes, probably due to local variations in surface microroughness, enhance the backscattered energy and may artificially create a topographic signal. The continental-scale study shows coherent patterns. Even if both surface and subsurface components affect the altimetric observation, the large-scale signal is mostly controlled by surface backscattering variations. The surface or near-subsurface characteristics of the snowpack may then be reached by altimetric observations.

Effects of the retreat of Gåsbreen, southwest Spitsbergen, Svalbard, on the evolution of the ice-dammed lake Goësvatnet are shown for the period 1899-1991. The retreat and lowering of the damming ice masses have changed not only the stored lake volume, the lake geometry and the elevation, slope and length of the subglacial outlet tunnel, but also the frequency and magnitude of outburst floods of Goësvatnet. For the estimation of peak discharges of outburst floods we computed an unbiased regression equation related to the progressive enlargement of subglacial tunnels using lake volume data and peak discharge data from the literature. The derived equation is very similar to the original form of the Clague-Mathews formula and answers the question why this formula has worked well in many cases. Peak discharges of Goësvatnet in various years were estimated by means of the derived equation. Effects of the changed lake geometry as well as the changed length and slope of the subglacial outlet tunnel on the discharge during outbursts will be discussed by means of the Nye-Clarke model. Observation of an outburst of Goësvatnet in summer 1991 indicates that the outbursts may have been triggered by pressure decrease in the subglacial outlet tunnel during increased discharge, whereas flotation of the ice dam can be excluded.

The evolution of ice-sheet configuration and flow pattern in Fennoscandia through the last glacial cycle was reconstructed using a glacial geological inversion model, i.e. a theoretical model that formalises the procedure of using the landform record to reconstruct ice sheets. The model uses mapped flow traces and deglacial melt-water landforms, as well as relative chronologies derived from cross-cutting striae and till lineations, as input data. Flow-trace systems were classified into four types: (i) time-transgressive wet-bed deglacial fans, (ii) time-transgressive frozen-bed deglacial fans, (iii) surge fans, and (iv) synchronous non-deglacial (event) fans. Using relative chronologies and aggregation of fans into glaciologically plausible patterns, a series of ice-sheet Configurations at different time slices was erected. A chronology was constructed through correlation with dated stratigraphical records and proxy data reflecting global ice volume. Geological evidence exists for several discrete ice-sheet configurations centred over the Scandinavian mountain range during the early Weichselian. The build-up of the main Weichselian Fennoscandian ice sheet started at approximately 70 Ka, and our results indicate that it was characterised by an ice sheet with a centre of mass located over southern Norway. This configuration had a flow pattern which is poorly reproduced by current numerical models of the Fennoscandian ice sheet. At the Last Glacial Maximum the main ice divide was located overthe Gulf of Bothnia. A major bend in the ice divide was caused by outflow of ice to the northwest over the lowest part of the Scandinavian mountain chain. Widespread areas of preserved pre-late-Weichselian landscapes indicate that the ice sheet had a frozen-bed core area, which was only partly diminished in size by inward-transgressive wet-bed zones during the decay phase.

A new deep ice-core drilling site has been identified in north Greenland at 75.12° N, 42.30° W, 316 km north-northwest (NNW) of the GRIР drill site on the summit of the ice sheet. The ice thickness here is 3085 m; the surface elevation is 2919 m.The North GRIP (NGRIP) site is identified so that ice of Eemian age (115–130 ka BP,calendar years before present) is located as far above bedrock as possible and so the thickness of the Eemian layer is as great as possible. An ice-flow model, similar to the one used to date the GRIP ice core, is used to simulate the flow along the NNW-trending ice ridge. Surface and bedrock elevations, surface accumulation-rate distribution and radio-echo sounding along the ridge have been used as model input.The surface accumulation rate drops from 0.23 m fee equivalent year−1 at GRIP to 0.19 m ice equivalent year−1 50 km from GRIP. Over the following 300km the accumulation is relatively constant, before it starts decreasing again further north. Ice thicknesses up to 3250 m bring the temperature of the basal ice up to the pressure-melting point 100–250 km from GRIP. The NGRIP site islocated 316 km from GRIP in a region where the bedrock is smooth and the accumulation rate is 0.19 m ice equivalent year−1. The modeled basal ice here has always been a few degrees below the pressure-melting point. Internal radio-echo sounding horizons can be traced between the GRIP and NGRIP sites, allowing us to date the ice down to 2300 m depth (52 ka BP). An ice-flow model predicts that the Eemian-age ice will be located in the depth range 2710–2800 m, which is 285 m above the bedrock. This is 120 m further above the bedrock, and the thickness of the Eemian layer of ice is 20 m thicker, than at the GRIP ice-core site.

The three-dimensional ice-sheet model SICOPOLIS is used to simulate the dynamic/thermody namic behaviour of the entire Greenland ice sheet from 250 000 a BP until today. External forcing consists of a surface-temperature history constructed from δ18O data of the GRIP core, a snowfall history coupled linearly to that of the surface temperature, a piecewise linear sea-level scenario and a constant geothermal heat flux. The simulated Greenland ice sheet is investigated in the vicinity of Summit, the position where the maximum elevation is taken, and where the two drill sites GRIP and GISP2 are situated 28km apart from each other. In this region, the agreement between modelled and observed topography and ice temperature turns out to be very good. Computed age-depth profiles for GRIP and GISP2 are presented, which can he used to complete the dating of these cores in the deeper regions where annual-layer counting is not possible. However, artificial diffusion influences the computed ages in a near-basal boundary layer of approximately 15% of the ice thickness, so that the age at the bottom of the cores cannot be predicted yet.

Observations are reported on cracks formed during compressive, unidirectional, constant-strain-rate deformation of columnar-grain ice. The axis of hexagonal crystallographic symmetry of each grain tended to be in the plane perpendicular to the long direction of the grains and to have a random orientation in that plane. For stress applied perpendicular to the long direction of the grains, the deformation was practically two-dimensional. It was found that the relative proportion of grain-boundary cracks increased with increasing strain rate, decreasing temperature and, for strain rate greater than 7 × 10−5 s−1, with decreasing grain-size. Almost all the grain-boundary cracks had at least one edge at a triple point. For each test, the grain-boundary and transcrystalline crack lengths tended to have a log-normal distribution. The logarithmic mean crack length (LMCL) decreased with increasing strain rate, decreasing grain-size and decreasing temperature and tended to a constant value of 0.75 mm at 10°C. For grain-size of 3 mm or greater, the LMCL had a maximum at a strain rate of 10−5 to 10−6 S−1 at −10°C. The LMCLs and the relative proportion of grain-boundary cracks tended to be normally distributed for given load conditions.

A degree-day glacier mass-balance model is coupled to a dynamic glacier model for temperate glaciers. The model is calibrated for two outlet glaciers from the Hofsjökull ice cap in central Iceland. It is forced with a climate scenario that has recently been defined for the Nordic countries for the purpose of outlining the hydrological consequences of future greenhouse warming. The scenario for Iceland specifies a warming rate of 0.25°C per decade in mid-summer and 0.35°C per decade in mid-winter with a sinusoidal variation through the year. The volume of the glaciers is predicted to decrease by approximately 40% over the next century, and the glaciers essentially disappear during the next 200 years. Runoff from the area that is presently covered by the glaciers is predicted to increase by approximately 0.5 m a−1 30 years from now due to the reduction in the volume of the glaciers. The runoff increase reaches a flat maximum of 1.5–2.0 m a−1 100–150 years from now and levels off after that. The predicted runoff increase leads to a significant increase in the discharge of rivers fed by meltwater from the outlet glaciers of the ice cap and may have important consequences for the operation and planning of hydroelectric power plants in Iceland.

Findelengletscher, Switzerland, advanced about 250 m between 1979 and 1985, and retreated thereafter. Subglacial water pressure, surface velocity and surface strain rate were determined at several sites. The measurements were made early in the melt seasons of 1980, 1982, 1985 and 1994 and in the autumn of 1983 and the winter of 1984. Changes of surface geometry were assessed from aerial photographs.

The estimated basal shear stress changed little between 1982 and 1994. Nevertheless, large changes in the relationship of subglacial water pressure and surface velocity were observed, which cannot be reconciled with the most commonly used sliding law unless it is modified substantially. Consideration of possible reasons indícales that a change in the subglacial drainage system occurred, probably involving a change in the degree of cavity interconnection. Isolated cavities damp the variations in sliding velocity that normally result from changes in water pressure, because the pressure in isolated cavities decreases as the sliding speed increases. In contrast, by transmitting water-pressure fluctuations to a larger area of the bed, interconnected cavities amplify the effect of water-pressure fluctuations on sliding speed. Thus, we suggest that an observed decrease in velocity (for a given water pressure) between 1982 and 1994 was a consequence of a decrease in the interconnectedness of the subglacial cavity system.

To study near-surface heat flow on the Antarctíc ice sheet, snow temperatures were measured at South Pole Station to a depth of 3 m at 15 min intervals during most of 1992. Solar heating and water-vapor transport were negligible during the 6 month Winter, as was inter-grain net thermal radiation, leaving conduction as the dominant heat-transport mechanism. The rate of temperature change at depth over 15 min intervals was smaller than that at the surface, by one order of magnitude at 20 cm depth and two orders of magnitude at 1 m depth. A finite-difference model, with conduction as the only heat-transfer mechanism and measured temperatures as the upper and lower boundary conditions, was applied to foursets of three thermistors each. The thermal conductivity was estimated as that which minimized the difference between modeled and measured 15 min changes in temperatures at the center thermistor. The thermal conductivity obtained at shallow depths (above 40 cm) was lower than that given by existing parameterizations based on density, probably because the snow grains were freshly deposited, cold and poorly bonded. A model using only vertical conduction explains on average 87% ofthe observed 15 min temperature changes at less than 60 cm depth and 92% below 60 cm. The difference between modeled andmeasured temperature changes decreased with depth. The discrepancies between model and observation correlated more strongly with the air-snow temperature difference than with the product of that difference with the square of the wind speed,suggesting that the residual errors are due more to non-vertical conduction and to sub-grid-scale variabilis of the conductivity than to windpumping. The residual heating rate not explained by the model of vertical conduction exceeds 0.2 W m−3 only in the top 60 cm of the near-surface snow.

Two-dimensional, isothermal, non-stationary, cold-glacier dynamics are analysed by perturbation methods, when the ice creep is described by Glen’s flow law. This approach allows us in model the ice-divide dynamics separately from that of the glacier and to pose a simplified problem for the ice divide. To order of the glacier-aspect ratio, all unknown characteristics near the ice divide can be searched for in an infinite layer with parallel boundaries, whose thickness coincides with that found by the shallow-ice approximation. The problem for the ice divide is stationary and does not depend on the ice and glacier characteristics, such as a flow-law constant, ice thickness and accumulation rate. At the ice divide, the ice upper surface curvature is finite and the shallow-ice approximation is inadequate.

In bedlock covered by continental ice, the depth to electrically conductive saline waters enriched below permafrost can be defined by electromagnetic measurements. Thus, if the ice thickness is also measured we obtain estimates of permafrost thickness variations in the survey area. Moreover, using the relation between ice and permafrost thicknesses and approximations of the melting temperature of saline waters combined with the average annual temperature at the surface, we can estimate temperature variations within both the bedrock and the ice. These can further be used for estimating the temperature values at the ice base.

The method was developed for interpreting two geophysical field profiles measured on continental ice in western Dronning Maud Land, Antarctica. Gravity, magnetic and electromagnetic profiles were measured between the nunataks Basen and Plogen (≈22 km) and Basen and Fossilryggen (≈38km). In the Basen–Fossilryggen profile the electromagnetic soundings revealed good conductors in the bedrock beneath the ice. The significant correlation between the ice thickness and the depth of the conductor in the bedrock strongly suggests that the conductors are not sedimentary rock layers but rather saline groundwater in which salinity has been enhanced by permafrost processes. From the regression curve we can conclude that the temperature vs depth curve in the ice can be represented by a first-degree line and that the “warm” ice base of about −2° C is expected to occur at depths of about 650 m, the base of permafrost in outcropping bedrock being at about 1 km. The mean temperature gradients in the ice and in the bedrock can be estimated as being about 20°Ckm−1 and 13°C km−1 respectively.

In the Basen–Plogen profile, conductors were detected only in the beginning of the profile, where ice thickness was below 650m. In the main part of the profile where the ice thickness exceeds 650 m, no conductors were observed. In those areas our interpretation is that the icebase is “warm” and no saline waters enriched by permafrost are present.

Electrical conductivity profiling allows definition of the hydrology of glacier boreholes occupied by contrasting water types. An inexpensive, robust conductivity profiling system is described that allows detailed description of borehole hydrostratigraphy with reasonable reproducibility. Field dato demonstrate the need to distinguish between bottom water in boreholes and basal water in drainage systems at the glacier bed.