The aim of deep ice drilling on Law Dome, Antarctica, has been to exploit the special characteristics of Law Dome summit, i.e. low temperature and high accumulation near an ice divide, to obtain a high-resolution ice core for climatic/environmental studies of the Holocene and the Last Glacial Maximum (LGM). Drilling was completed in February 1993, when basal ice containing small fragments of rock was reached at a depth of 1196 m. Accurate ice dating, obtained by counting annual layers revealed by fine-detail δ18О, peroxide and electrical-conductivity measurements, is continuous down to 399 m, corresponding to a date of AD 1304. Sulphate concentration measurements, made around depths where conductivity tracing indicates volcanic fallout, allow confirmation of the dating (for Agung in 1963 and Tambora in 1815) or estimates of the eruption date from the ice dating (for the Kuwae, Vanuatu, eruption ~1457). The lower part of the core is dated by extrapolating the layer-counting using a simple model of the ice flow. At the LGM, ice-fabric measurements show a large decrease (250 to 14 mm2) in crystal size and a narrow maximum in c-axis vertically. The main zone of strong single-pole fabrics however, is located higher up in a broad zone around 900 m. Oxygen-isotope (δ18O) measurements show Holocene ice down to 1113 m, the LGM at 1133 m and warm (δ18O) about the same as Holocene) ice near the base of the ice sheet. The LGM/Holocene δ18O shift of 7.0‰, only ~1‰ larger than for Vostok, indicates that Law Dome remained an independent ice cap and was not overridden by the inland ice sheet in the Glacial.

During ice-structure interaction, ice will fail in a brittle manner dominated by two processes. The first corresponds to the formation of macrocracks and the consequent spalling-off of large ice pieces. The second includes an intense shear-damage process in zones, termed critical zones, where high pressures are transmitted to the structure. The shear-damage process results in microstructural changes including microcrack formation and recrystallization. A range of tests on laboratory-prepared granular ice have been conducted to determine the fundamental behaviour of ice under various stress states and stress history, particularly as it relates to changes in microstructure. The test series was designed to study three aspects: the intrinsic creep properties of intact, undamaged ice; the enhancement of creep and changes in microstructure due to damage; and the effects of different stress paths. Tests on intact ice with triaxial confining pressures and low deviatoric stresses, aimed at defining the intrinsic creep response in the absence of microcracking, showed that an accelerated creep rate occurred at relatively low deviatoric stresses. Hence, a minimum Creep rate occurred under these conditions. Recrystallization to a smaller grain-size and void formation were observed. Ice damaged uniaxially and triaxially prior to testing showed enhancement of creep under both uniaxial and triaxial loading conditions Creep rates in triaxially damaged ice were found to be non-linear with high deviatoric stresses, corresponding to a power-law dependence of creep rate. Uniaxially damaged specimens contained microcracks parallel to the stressed direction which tended to close under triaxial confinement. Damage under triaxial conditions at low confining pressures produced small recrystallized grains near zones of microcracking. At high confining pressures, a fine-grained recrystallized structure with no apparent cracking was observed uniformly across the specimen. The recrystallization process contributes significantly to the enhanced creep rates found in damaged specimens.

Twenty-seven studies on the thermal conductivity of snow (Keff) have been published since 1886. Combined, they comprise 354 values of Keff, and have been used to derive over 13 regression equation and predicting Keff vs density. Due to large (and largely undocumented) differences in measurement methods and accuracy, sample temperature and snow type, it is not possible to know what part of the variability in this data set is the result of snow microstructure. We present a new data set containing 488 measurements for which the temperature, type and measurement accuracy are known. A quadratic equation,

where ρ is in g cm−3, and Keff is in W m−1K−1, can be fit to the new data (R2 = 0.79). A logarithmic expression,

can also be used. The first regression is better when estimating values beyond the limits of the data; the second when estimating values for low-density snow. Within the data set, snow types resulting from kinetic growth show density-independent behavior. Rounded-grain and wind-blown snow show strong density dependence. The new data set has a higher mean value of density but a lower mean value of thermal conductivity than the old set. This shift is attributed to differences in snow types and sample temperatures in the sets. Using both data sets, we show that there are well-defined limits to the geometric configurations that natural seasonal snow can take.

Measurements from the subarctic snowpack are used to explore the relationship between grain growth and vapor flow, the fundamental processes of dry-snow metamorphism. Due to extreme temperature gradients, the subarctic pack undergoes extensive depth-hoar metamorphism. By the end of the winter a five-layered structure with a pronounced weak layer near the base of the snow evolves. Grain-size increases by a factor of 2–3. while the number of grains per unit mass decreases by a factor of 10. Observed growth rates require significant net inter-particle vapor fluxes. Stable-isotope ratios show that there are also significant net layer-to-layer vapor fluxes. Soil moisture enters the base of the pack and mixes with the bottom 10 cm of snow, while isotopically light water vapor fractionates from the basal layer and is deposited up to 50 cm higher in the pack. End-of-winter density profiles for snow on the ground, compared with snow on tables, indicate the net layer-to-layer vapor flux averages 6 x 10−7 kg m−2 s−1, though detailed measurements show the net flux is episodic and varies with time and height in the pack, with peak net fluxes ten limes higher than average. A model, driven by observed temperature profiles, reproduces the layer-to-layer flux pattern and predicts the observed weak layer at the base of the snow. Calculated layer-to-layer vapor fluxes are ten times higher than inter-particle fluxes, which implies that depth-hoar grain growth is limited by factors other than the vapor supply. This finding suggests that gain and loss of water molecules due to sublimation from grains takes place at a rate many times higher than the rate at which grains grow, and it explains why grains can metamorphose into different forms so readily.

Layering affects the air flow through snow caused by surface pressure variations. The horizontal and total fluxes are high in hoar layers but the pressure perturbations and vertical components of the flow do not penetrate as deeply as in homogeneous snow. That is because the layers “pipe” the flow horizontally toward the area of low pressure. An ice layer at the surface reduces the total flow everywhere. The flow decreases as ice-layer thickness increases and, in general, flow changes with permeability. However, the magnitude of the effect is proportionately weaker when the ice layers are further from the surface. The residence time is reduced when hoar layers are present due to shorter flow paths, reduced penetration into the deeper snow and higher speeds.

Parameterizing the incoming longwave radiation L ↓ in terms of the fourth power of the absolute temperature at the reference height is used in glaciology for several purposes. In this paper, the validity of this kind of parameterization is investigated for the Greenland ice sheer, both by observations and by numerical simulation with a meso-scale model, It is found that such a parameterization severely underestimates the increase of L ↓ in response to large-scale warming in an area where surface melting is important. This is explained by the systematic influence that is exerted on the shape of the temperature profiles by surface melting.

The basal deformation of a gravity-driven linear creeping flow sliding frictionless over slowly varying bed undulations in two dimensions is analysed analytically, using results from second-order perturbation theory. One of the key results is that, close to sinusoidal bedrock undulations, up to two different spatial regions of local extrusion flow may arise. The offset and onset of extrusion flow is controlled primarily by the amplitude-to-wavelength ratio. Above the crest of a sinusoidal bed line, a local maximum of the surface-parallel velocity develops for ε : = ak < 0.138, where a is the amplitude and k is the wave number. As ε increases from zerо to this critical value, the vertical position of the velocity maximum moves from kz = 1 to kz ≈ 1.98, where z is the vertical distance above the mean bed line. Within and above the trough of a sinusoid, a region of local minimum of the surface-parallel velocity component develops, which shifts from kz = 1 towards the bed line as ε increases front zero to 1/2. Below this velocity minimum, and for some distance above the velocity maximum, the surface-parallel velocity increases with depth. This type of extrusion flow will cause a reversal of borehole-inclination profiles close to the bedrock.

The flow field of a medium sliding without friction over a strongly undulating surface is calculated numerically. The results are used to elucidate the basal-flow characteristics of glacier flow and they are discussed with reference to known analytical solutions. Extrusion flow is found to become increasingly pronounced as the value of n, where n is a parameter in Glen’s flow law, becomes larger. For sinusoidal bedrock undulations, a flow separation occurs if the amplitude-to-wavelength ratio exceeds a critical value of about 0.28. The main flow then sets up a secondary flow circulation within the trough, and the ice participating in this circular motion theoretically never leaves it. The sliding velocity is calculated numerically as a function of the mean basal shear stress, the amplitude-to-wavelength ratio and the flow parameter n. For moderate and high slope fluctuations, the sliding velocity is significantly different from what would be expected from results based on the small-slope approximation.

Poor consideration has been given in many Arctic circum-polar ice-core studies to the effect of summer snow melt on chemistry, stable-isotope concentrations and time-scales. Many of these corps are drilled close to the firn line where melt is intense. Some come from below the firn line where accumulation is solely in the form of super-imposed ice. In all cases, seasonal signals are reduced or removed and, in some, time gaps develop during periods of excessive melting which situate the drill site in the ablation zone. Consequently, cross correlations of assumed synchronous events among the cores are invalid, so that time-scales along the same cores differ between authors by factors of over 2. Many so-called climatic signals are imaginary rather than real. By reference to published analyses of cores from the superimposed ice zone on Devon Ice Cap (Koerner, 1970) and Meighen Ice Cap (Koerner and Paterson, 1974), it is shown how melt affects all the normally well-established ice-core proxies and leads to their misinterpretation. Despite these limitations, the cores can give valuable low-resolution records for all or part of the Holocene. They show that the thermal maximum in the circum-polar Arctic occurred in the early Holocene. This maximum, effected negative balances on all the ice caps and removed the smaller ones. Cooler conditions in the second half of the Holocene have caused the regrowth of these same ice caps.

Changes along the margin of the southern half of the 5900 km2 Barnes Ice Cap have been assessed using 1993 Landsat TM imagery in comparison with digitized 1:50 000 NTS maps based on 1961 photogrammetry. The average recession over the 183 km long southern perimeter was found to be at least 4 m a−1, with no significant difference between the southeast and southwest sectors. Viewed in conjunction with the sustained retreat previously reported for the northwest margin, these results indicate that a general reduction in the size of Barnes Ice Cap is occurring. The present retreat phase began under a regional climate warming in the late 19th to early 20th century period and continues, while the record of the ablation-season temperature since the mid-century has not shown any significant trend.

Response tests are widely used in ground-water studies to assess the hydraulic properties of sub-surface water-flow systems. The simplicity of such tests also makes them attractive for investigation of subglacial hydraulic conditions. This paper describes a systematic, quantitative approach to the analysis of borehole-response test data. The approach uses the theoretical model of Stone and Clarke (1993), which describes water motion in a coupled borehole—subglacial flow system; this framework provides the basis for an inversion scheme that is focused on quantifying physical properties of the basal-flow system, as it is characterized in the theoretical model. The inversion procedure was applied to response-test data from Trapridge Glacier, Yukon Territory, Canada. Results of the inversions suggest that the subglacial drainage network can be described as a confined layer comprising coarse-sand-to fine-gravel-sized sediments, having a thickness of 0.1 – 0.3 m, and a hydraulic conductivity of about 5 × 10−4ms−1. Based on the water-drainage rates from boreholes, as they connect with the subglacial water-flow system, specific storage of the sediment layer was estimated to be approximately 1 × 10−4m−1. Further consideration of subglacial water-flow conditions suggests that connection drainage test results may tend to underestimate specific storage of the overall glacier substrate.

Satellite imagery indicates that the floating terminus of Pine Island Glacier has changed little in extent over the past two decades. Data on the velocity and thickness of the glacier reveal that calving of 28 ± 4 Gta−1 accounts for only half of the ice input near the grounding line. The apparently steady configuration implies that the remainder of the input is lost by basal melting at a mean rate of 12 ± 3 ma−1. Ocean circulation in Pine Island Bay transports +1°C waters beneath the glacier and temperatures recorded in melt-laden outflows show that heat loss from the ocean is consistent with the requirements of the calculated melt rate. The combination of iceberg calving and basal melting lies at the lower end of estimates for the total accumulation over the catchment basin, drawing into question previous estimates of a significantly positive mass budget for this part of the ice sheet.

The surface energy budget was investigated during a cruise through the pack ice in the Southern Ocean. The time of observation was close to mid-summer. Some of the more important findings were:

Glacial mass-balance reconstruction for a long-term time-scale requires knowledge of the relation between climate change and mass-balance fluctuations. A large number of mass-balance reconstructions since the beginning of the century are based on statistical relations between monthly meteorological data and mass balance. The question examined in this paper is: are these relationships reliable enough for long-term time-scale extrapolation? From the glacier de Sarennes long mass-balance observations series, we were surprised to discover large discrepancies between relations resulting from different time periods. The importance of the albedo in relation to ablation and mass balance is highlighted, and it is shown that it is impossible to ignore glacier-surface conditions in establishing the empirical relation between mass-balance fluctuations and climatic variation; to omit this parameter leads to incorrect results for mass-balance reconstruction in the past based on meteorological data.

Sixty-three ice cores were collected in the Bellingshausen and Amundsen Seas in August and September 1993 during a cruise of the R.V. Nathaniel B. Palmer. The structure and stable-isotopic composition (18O/16O) of the cores were investigated in order to understand the growth conditions and to identify the key growth processes, particularly the contribution of snow to sea-ice formation. The structure and isotopic composition of a set of 12 cores that was collected for the same purpose in the Bellingshausen Sea in March 1992 are reassessed. Frazil ice and congelation ice contribute 44% and 26%, respectively, to the composition of both the winter and summer ice-core sets, evidence that the relatively calm conditions that favour congelation-ice formation are neither as common nor as prolonged as the more turbulent conditions that favour frazil-ice growth and pancake-ice formation. Both frazil- and congelation-ice layers have an av erage thickness of 0.12 m in winter, evidence that congelation ice and pancake ice thicken primarily by dynamic processes. The thermodynamic development of the ice cover relies heavily on the formation of snow ice at the surface of floes after sea water has flooded the snow cover. Snow-ice layers have a mean thickness of 0.20 and 0.28 m in the winter and summer cores, respectively, and the contribution of snow ice to the winter (24%) and summer (16%) core sets exceeds most quantities that have been reported previously in other Antarctic pack-ice zones. The thickness and quantity of snow ice may be due to a combination of high snow-accumulation rates and snow loads, environmental conditions that favour a warm ice cover in which brine convection between the bottom and top of the ice introduces sea water to the snow/ice interface, and bottom melting losses being compensated by snow-ice formation. Layers of superimposed ice at the top of each of the summer cores make up 4.6% of the ice that was examined and they increase by a factor of 3 the quantity of snow entrained in the ice. The accumulation of superimposed ice is evidence that melting in the snow cover on Antarctic sea-ice floes ran reach an advanced stage and contribute a significant amount of snow to the total ice mass.

The historical length variations in Unterer Grindelwaldgletscher have been simulated by coupling a numerical mass-balance model to a dynamic ice-flow model. As forcing functions, we used (partly reconstructed) local climatic records, which were transformed by the mass-balance model into a mass-balance history. The ice-flow model then computes the length variations that have occurred over the course of time.

In a model run from AD 1530 to the present, with both seasonal temperature and precipitation variations as forcing functions, the observed maximum length of the glacier around AD 1860 and the subsequent retreat are simulated. The observed AD 1600 maximum, however, does not show up in the simulation. This is probably due to an incorrect reconstruction of the mass balance for this period, as detailed climatic data are available only since 1865. The root-mean-square difference between the simulated and the observed front positions is 0.28 km. The simulated glacier geometry for 1987 fits the observed geometry for that year reasonably well.

Because of the success of the historical simulation, an attempt is made to predict future glacier retreat on the basis of two different greenhouse-gas scenarios. For a Business-as-Usual scenario, only 29% of the 1990 volume would remain in AD 2100.

The frequency of avalanches at a given location is the primary variable for calculating the risk as input to zoning applications and decisions about avalanche-control options. In this paper, we present an in-depth study of avalanche frequency using an extensive data base of avalanche-occurrence records from Rogers’ Pass, British Columbia (43 avalanche paths; 24 years of records). This study, the first of its kind for high-frequency avalanche paths, yields the result that the frequency of avalanches may be described by a Poisson distribution. Study of the relationship between terrain variables and precipitation estimates shows that avalanche frequency is significantly correlated with path roughness, 30 year maximum water equivalent, east-west location from Rogers’ Pass summit, wind exposure and run-out zone elevation and inclination. With the length of avalanche-occurrence records and quality of the data, we believe our study is the most comprehensive in existence about avalanche frequency and its relation to terrain variables.

In order to study, in situ, the rheology of a deforming subglacial till, various instruments were emplaced in till beneath Storglaciären, Sweden. Boreholes were used to gain access to the till beneath about 100 m of ice. Tiltmeters provided an estimate of the shear strain rate in the till. Two other instruments yielded measures of till strength. In addition, water pressures were recorded in boreholes and in the till, a computer-controlled distance meter provided an effectively continuous record of the surface velocity and data from frequent surveys of a stake network were used to estimate the mean basal drag, based on a force-balance calculation.

Tilt rates varied directly with effective pressure, so decreases in water pressure apparently increased the coupling between the glacier and the bed. Surface speed was either out of phase with tilt or varied independently of tilt. Thus, increases in speed were apparently a consequence either of longitudinal coupling or of reduced coupling between the glacier and the bed; they were not a result of till deformation! Till strength varied directly with effective pressure, which is consistent with it being a Mohr – Coulomb, or frictional material. The devices measuring till strength are presumed to have been pulled through the till at a speed that varied in phase with the surface speed but till strength did not vary systematically with surface speed. This implies that the residual strength of the till is insensitive to strain rate. Thus, the appropriate constitutive equation for till rheology may be of the form:

where k is a constant. This is consistent with experimental data reported in the geotechnical literature.

A physics-based snow model has been calibrated using data collected at Halley Bay, Antarctica, during the International Geophysical Year. Variations in snow temperature and density are well-simulated using values for the model parameters within the range reported from other polar field experiments. The effect of uncertainty in the parameter values on the accuracy of the predictions is no greater than the effect of instrumental error in the input data. Thus, this model can be used with parameters determined a priori rather than by optimization. The model has been validated using an independent data set from Halley Bay and then used to estimate 10 m temperatures on the Antarctic Peninsula plateau over the last half-century.