The Vavilov ice cap was perforated in 1988 by a drilling which reached the underlying frozen sediments. In contrast to the overlying glacier ice, the basal ice is composed of different ice layers with a variable debris load. The stable-isotope composition of these layers shows δ values much lower than everywhere else in the core or in the Vavilov ice cap. This is most probably the signature of a remnant of Pleistocene ice which, for the first time, is shown to occur in the Russian Arctic.

We use the Wetherald and Manabe climate model to predict the response of mountain glaciers to a doubling of atmospheric carbon dioxide. The response is measured in terms of a change in the equilibrium-line altitude (ELA) and the glacier terminus altitude (GTA), net accumulation–ablation on these altitudes and the melt runoff for 12 mountain-glacier regions in southeastern Alaska, U.S.A. The methods we use involve extrapolating climate-model temperature fields to a glacier’s location, and empirical–statistical relationships between air temperature and percentage of solid precipitation, and between summer air temperature and ablation and melt runoff. Our study shows that, under global warming, glaciation in southeastern Alaska will not disappear, but mass exchange of glaciers will be more intensive and the ELA value will increase by 300–760 m, depending on the glacier’s distance from the ocean.

Stresses and velocities at depth are calculated across Jutulstraumen, an ice stream in Dronning Maud Land, draining about 1% of the Antarctic ice sheet. The force-balance study is based on data from kinematic GPS measurements on three strain nets, each consisting of 3 × 3 stakes. The maximum measured velocity is 443 m a−1 and the velocity variation over short distances is large compared with studied ice streams in West Antarctica. The surface topography together with the measured velocities across the profile indicate that the bottom topography has a great influence on the flow direction, even where the ice thickness is more than 2000 m. The basal shear stresses are calculated as 180, 227 and 146 kPa in the three Strain nets, while the corresponding driving stresses are 180, 122 and 111 kPa (±5%). The heat produced by sliding and internal deformation is sufficient to keep the base at the pressure-melting point. The annual basal melting is estimated to be about 60 mm. Investigations on the effect of temperature softening show that the flow parameter’s influence on the effective strain rate is more important than the flow parameter’s direct softening in the flow low alone. The mass flow calculated by the force-balance method is between 87 and 96% of pure plug flow and total discharge is calculated to be 13.3 ± 10 km3a-1.

Following construction of a glacial ice runway on the Ross Ice Shelf, Antarctica, and prior to flight operations, the runway was proof-rolled. The proof exercise was designed to simulate typical heavy aircraft. Initial testing produced numerous brittle surface failures in the runway ice. Thin sections of ice cores taken from the failed areas showed large crystals (с axis vertical) of clear, blue ice with long, vertical bubbles, indicative of ice formed directly from meltwater. Uniaxial unconfined compression tests on core samples were used to compare runway ice strength with published data for polycrystalline laboratory ice. Since the frequent failure of surface ice had not been expected, it was critical to understand the formation and mechanical properties of the weak ice to prevent its occurrence in the future and to strengthen the existing problem areas. We discuss the likely scenarios for development of weak ice on the airstrip and the physical properties of this type of ice. Also, the procedure used to repair successfully the runway surface is described, which culminated in test flights, followed by full flight operations.

40Ar/39Ar ages of most single ice-ratted amphiboles from Heinrich layer 2 (H2) from a core in the Labrador Sea, a core in the eastern North Atlantic and a core in the western North Atlantic range from 1600 to 2000 Ma. This range is identical to that for K/Ar ages from the Churchill province of the Canadian Shield that outcrops at Hudson Strait and forms the basement of the northern part of Hudson Bay. The ambient glacial sediment includes some younger and older grains derived from Paleozoic, Mesoproterozoic and Archean sources, but still the majority of the amphiboles have ages in the 1600–2000 Ma interval. The Ca/K ratios of these 1600–2000 Ma old amphiboles, however, have a bimodal distribution in contrast with the uniformity of the Ca/K ratios of H2 amphiboles. This indicates that 1600–2000 Ma old amphiboles of the ambient sediment were derived from an additional Early Proterozoic source besides Churchill province. In H2, Churchill-derived grains constitute 20–40% of the ice-rafted debris (IRD). The fraction in the ambient glacial sediment is 65–80%. Results presented here are consistent with the hypothesis that Heinrich events were produced by a sudden intensification of the iceberg discharge through Hudson Strait that mixed, in the North Atlantic, with icebergs that continued to calve from other ice sheets. The shift from mixed sources in the background sediment to a large dominance of Churchill province grains in H2 indicates that, even if calving of other ice sheets intensified during the Heinrich episode, the increase in the iceberg discharge via Hudson Strait from the Hudson Bay drainage basin of the Laurentide ice sheet was by far the largest.

Two large glacier tongues, which extend substantially across the coastline of King George V Land in East Antarctica, have been studied by remote sensing (synthetic aperture radar, JERS-1). The tongue of Mertz Glacier is in a state of advance, while the Ninnis Glacier tongue is retreating. Some more specific points are:

The relations between mass balance and meltwater refreezing were examined on the basis of glaciological observations carried out in summer 1993 on Xiao Dongkemadi Glacier, Tanggula Mountains, central Tibetan Plateau. On this glacier, a part of meltwaler refreezes at the snow/ice interface as superimposed ice. The amount of superimposed ice formation was determined by both meltwater supply and temperature condition of the glacier. Snow-layer thickness on the glacier ice body is less than 2 m, even in the higher accumulation zone. About 60% of meltwaler generated in the accumulation zone for the period May–September was trapped at the snow/ice interface by refreezing, and was not discharged out of the glacier. About 26% of accumulated snow to the glacier surface was replaced on the snow/ice interface by refreezing in the accumulation zone. These facts indicate that superimposed ice formation is quite significant for water retention in glaciers under low-precipitation conditions.

Over 75 000 surface-velocity measurements are extracted from sequential satellite imagery of Ice Streams D and E to reveal a complex pattern of flow not apparent from previous measurements. Horizontal and vertical strain rates, calculated from surface velocity, indicate that the bed experiences larger basal shear where the surface of these ice streams is rougher. Ten airborne-radar profiles and one surface-based radar profile of ice thickness make possible the calculation of mass balance for longitudinal sections of each ice stream. Improved data-collection methods increase data density, substantially reducing random errors in velocity. However, systematic errors continue to limit the ability of the flux-differencing technique used here to resolve local variations in mass balance. Nevertheless, significant local variations in mass balance are revealed, while, overall, Ice Streams D and E are in approximate equilibrium. An earlier estimate of the net mass balance for Ice Stream D is improved.

Over a floating glacier ice tongue or an ice shelf, the glacier motion measured by a single, repeat-pass, radar interferogram is difficult to analyze, because the long-term, steady motion of the ice is intermixed with its cyclic, downward motion induced by tidal forcing. Multiple interferograms and a quadruple-difference technique are necessary to separate the tidal signal from the long-term, steady motion of the ice. An example of application of this technique is given here using ERS-1 radar images of Petermann Gletscher, a major outlet glacier of northern Greenland. Tidal displacements are measured with < 5 mm statistical noise. The long-term ice Velocity is measured with a precision of 1 ma−1. The inferred tidal displacements agree well with model predictions from a fixed elastic beam with an elastic damping factor of 0.47 ± 0.01 km3. The hinge line is mapped with a precision of 20-80 m.

Combining the interferometric ice velocities with ice thickness data, the glacier ice discharge is calculated at and below the hinge line. At the hinge line, the ice flux is 12.1 ± 1 km3 a−1. At the ice front, calf-ice production is only 0.59 km3−1 a−1, meaning that 95% of the ice that crosses the grounding line melts before it reaches the calving front. Assuming steady-state conditions, the melt rate of the glacier tongue averages 12 ± 1 m a−1, with peak values exceeding 20 ma−1 near the hinge line. This high melt rate cannot be accommodated by surface ablation alone (only about 23 ma−1) and is attributed to pronounced basal melting of the ice tongue. Basal melting, often assumed to be negligible in Greenland, is the dominant process of mass release from the floating section of Petermann Gletscher.

We present new interpretations of deglaciation in McMurdo Sound and the western Ross Sea, with observationally based reconstructions of interactions between East and West Antarctic ice at the last glacial maximum (LGM), 16000, 12000, 8000 and 4000 BP. At the LGM, East Antarctic ice from Mulock Glacier split; one branch turned westward south of Ross Island but the other branch rounded Ross Island before flowing southwest into McMurdo Sound. This flow regime, constrained by an ice saddle north of Ross Island, is consistent with the reconstruction of Stuiver and others (1981a). After the LGM, grounding-line retreat was most rapid in areas with greatest water depth, especially along the Victoria Land coast. By 12000 BP, the ice-now regime in McMurdo Sound changed to through-flowing Mulock Glacier ice, with lesser contributions from Koettlitz, Blue and Ferrar Glaciers, because the former ice saddle north of Ross Island was replaced by a dome. The modern flew regime was established ∼4000 BP. Ice derived from high elevations on the Polar Plateau but now stranded on the McMurdo Ice Shelf, and the pattern of the Transantarctic Mountains erratics support our reconstructions of Mulock Glacier ice rounding Minna Bluff but with all ice from Skelton Glacier ablating south of the bluff. They are inconsistent with Drewry’s (1979) LGM reconstruction that includes Skelton Glacier ice in the McMurdo-Sound through-flow. Drewry’s (1979) model closely approximates our results for 12000-4000 BP. Ice-sheet modeling holds promise for determining whether deglaciation proceeded by grounding-line retreat of an ice sheet that was largely stagnant, because it never approached equilibrium flowline profiles after the Ross Ice Shelf grounded, or of a dynamic ice sheet with flowline profiles kept low by active ice streams that extended northward from present-day outlet glaciers after the Ross Ice Shelf grounded.

Antarctic subglacial lakes provide аn important boundary condition for thermal analysis of the ice sheet in that the basal ice temperature over lakes may be assumed to be at the pressure-melting point. We have used a one-dimensional vertical heat-transfer equation to determine theoretical temperature values for the ice-sheet base above 77 subglacial lakes identified from airborne radio-echo-sounding data covering 50% of Antarctica. Variations in our temperature results to below the pressure-melting temperature over lakes are due to either our estimate of the geothermal heat flux or a neglect of heat derived from (a) internal ice deformation and (b) basal sliding, in the thermal model. Our results indicate that, when the geothermal heat flux is set at 54 m W m−2, the ice-sheet base above 70% of the known Antarctic subglacial lakes is calculated to be at the pressure-melting value. These lakes are located mainly around Dome C, Ridge B and Vostok station. For the ice sheet above subglacial lakes located hundreds of kilometres from the ice divide, using the same thermal model, loss of heat due to vertical advection is calculated to be relatively high. In such regions, if the ice-sheet base is at the pressure-melting point, heat lost due to vertical advection must be supplemented by heat from other sources. For the three lakes beneath Terre Adélie and George V Land, for instance, the basal thermal gradient calculated to produce pressure melting at the ice-sheet base is equivalent to 1.5–2 times the value obtained when 54 m W m−2 of geothermal heat is used as the sole basal thermal component. We suggest that, as distance from the ice divide increases, so too does the amount of heat due to internal ice deformation and basal sliding. Moreover, by considering the ice-sheet basal thermal characteristics above subglacial lakes which lie on the same ice flowline, we demonstrate empirically that the heat due to these horizontal ice-motion terms varies pseudo-exponentially with distance from the ice divide. The location along a flowline where a rapid increase in the basal heat gradient is calculated may correspond to the onset of large-scale basal sliding.

Geodetic surveying and ground-based radar profiling were used to determine geometry and surface motion of the ice sheet on the Dyer Plateau, Antarctica, in the vicinity of an ice-core site on a local dome. Vertical strain measurements in the core hole constrain the depth profile of vertical velocity. These geophysical measurements are used to analyze the profiles of density and annual layer thickness measured on the ice core to estimate the current mass balance of the ice column and the past history of accumulation rate. Consideration of horizontal and vertical mass-flow divergence shows that the profiles of density and vertical velocity are not fully consistent with steady state. Mean density of the firn layer appears to be increasing, which leads to the deduction of a small rate of mass increase (≈ 0.02 m a− 1 ice-equivalent thickness). Over the last 200a there has been a gradual increase in accumulation rate from about 0.46 m a− 1 to 0.54 m a− 1 ice-equivalent thickness in recent time.

To obtain a better understanding of the compliance of saline ice, we developed a simple conceptual model of a process that presumably takes place in saline ice. Using elastic models, we assessed the role that phase transformations might play when brine is sealed into small cells during crystal growth. Cooling of a sealed brine cell, leading to the precipitation of ice, provides a mechanism for the accumulation of “large” localized stresses. Based on our analysis, this mechanism has the potential to be a significant source for the nucleation of dislocations, and can conceivably make a significant contribution to the greater compliance (softness) of saline ice relative to non-saline ice.

The results of this model are consistent with the observation that laboratory-grown saline crystals sometimes display extensive differences in mechanical behavior that appear to be due to variations in the growth and storage conditions experienced by the crystals.

Radar-sounding data taken at Uvérsbreen glacier, Svalbard, are used to estimate the relative water content of the temperate part of the glacier. The relative change in water content of the temperate part of the glacier. The relative change in water content is estimated from the strength of the backscttered radar energy without knowing the size of the scatterers. The main findings are that the water content is on average changing by a factor of five in the ice, with the exception of some bright spots believed to stem fro larger water moulins. The general trend is also that the water content is increasing with lower altitude. The applied technique may be a useful method for large-scale mapping of the relative water-content variation in temperate and polythermal glaciers.

A relatively lightweight and simple airborne system for surface elevation profiling of glaciers in narrow mountain valleys has been developed and tested. The aircraft position is determined by kinematic global positioning system (GPS) methods. The distance to the glacier surface is determined with a laser ranger. The accuracy is about 0.3 m, sufficient to permit future changes to be observed over short time intervals. Long-term changes can be estimated by comparison of profiles with existing maps. Elevation profiles obtained in 1993–94 from three glaciers in central and south-central Alaska are compared with maps made about 1950. The resulting area-averaged, seasonally corrected thickness changes during the interval are: Gulkana Glacier (central Alaska Range)–11 m, Worthington Glacier (central Chugach Mountains) +7 m, and Bear Lake Glacier (Kenai Mountains) −12 m. All three glaciers retreated during the interval of comparison. The estimated uncertainty in the average thickness change is ±5 m. which is mainly due to errors in the existing maps. Constraints on the accuracy of the maps are obtained by profiling in proglacial areas.

White Glacier is a valley glacier at 79.5°N with an area of 38.7 km2. Its mass balance has been measured, over 32 years with a 3 year gap, by standard techniques using the stratigraphic system with a stake density of the order of one stake per km2. Errors in stake mass balance are about ±(200–250) mm, due largely to the local unrepresentativeness of measurements. Errors in the whole-glacier mass balance B are of the same order as single-slake errors. However, the lag-1 autocorrelation in the time series of B is effectively zero, so it consists of independent random samples, and the error in the long-term “balance normal” 〈B〉 is noticeably less. 〈B〉 is −100 ± 48 mm. The equilibrium-line altitude (ELA) averages 970 m. with a range of 470–1400 m. Mass balance is well correlated with ELA, but detailed modelling shows that the equilibrium line is undetectable on visible-band satellite images. A reduced network of a few stakes could give acceptable but less accurate estimates of the mass balance, as could estimates based on data from a weather station 120 km away. There is no evidence of a trend in the mass balance of White Glacier. To detect a climatologically plausible trend will require a ten-fold reduction of measurement error, a conclusion which may well apply to most estimates of mass balance based on similar stake densities.

Satellite radar interferometry provides glaciologists with an important new tool for determining the motion and topography of large ice sheets. We examine the sources of error in interferometrically derived ice-motion measurements, including those errors due to inaccurate estimates of the interfero-metric baseline. Several simulations are used to assess baseline accuracy in terms of tie-point error and the number and distribution of tie points. These results give insight into how best to select tie points, and also demonstrate the level of accuracy that can be achieved. Examination of two representative cases likely to occur in mapping ice-sheet motion leads to the conclusion that with adequate tie-point information ice velocity can be measured accurately to within a few meters per year. A method to correct horizontal velocity estimates for the effect of vertical displacement using surface slopes is also developed. Finally, we estimate the single-component velocity field for an area on Humboldt Glacier, northern Greenland, using interferograms formed from ERS-1 SAR image. We estimatе that these velocity measurements are accurate to within 2.3 m year1.

We describe a high-precision (0.1–1.0 mK) borehole-temperature (BT) logging system developed at the United States Geological Survey (USGS) for use in remote polar regions. We discuss calibration, operational and data-processing procedures, and present an analysis of the measurement errors. The system is modular to facilitate calibration procedures and field repairs. By interchanging logging cables and temperature sensors, measurements can be made in either shallow air-filled boreholes or liquid-filled holes up to 7 km deep. Data can be acquired in either incremental or continuous-logging modes. The precision of data collected by the new logging system is high enough to detect and quantify various thermal effects at the milli-Kelvin level. To illustrate this capability, we present sample data from the 3 km deep borehole at GISP2, Greenland, and from a 130 m deep air-filled hole at Taylor Dome, Antarctica. The precision of the processed GISP2 continuous temperature logs is 0.25–0.34 mK, while the accuracy is estimated to be 4.5 mK. The effects of fluid convection and the dissipation of the thermal disturbance caused by drilling the borehole are clearly visible in the data. The precision of the incremental Taylor Dome measurements varies from 0.11 to 0.32 mK. depending on the wind strength during the experiments. With this precision, we found that temperature fluctuations and multi-hour trends in the BT measurements correlate well with atmospheric-pressure changes.

Large, single ice crystals containing nо air bubbles and free of both small-angle grain boundaries and visible stress birefringence can be grown using a very simple growth chamber within a temperature-controlled, outer enclosure. The method relies upon the spontaneous formation of an ice crystal with its с axis accurately normal to a free, slightly supercooled, water surface.