“Fresh-looking” moraines lying within a mile of glaciers in the lowlands of the McMurdo Sound area, Victoria Land, Antarctica have been interpreted by some to be only 50 or 100 yr. old and to have been left by recent glacier retreats. The present author believes that these moraines are considerably older because a comparison of the position of existing glacier fronts in the area with positions on photographs of 50 yr. ago shows no movement of the front or appreciable thickening or thinning of the glaciers during the last 50 yr.

Radiocarbon dating of mummified seal carcasses lying near glacier fronts indicates that the glaciers have not been more extensive for at least 1,000 yr. The youngest moraines in the area are ice-cored moraines to the existing glaciers. Radiocarbon dates of algae from extinct ephemeral ponds in the ice-cored moraines indicate the moraines to be at least 6,000 yr. old.

Additional evidence for the antiquity of the near-glacier moraines is the presence of a rather uniformly well-developed micro-relief pattern of sand-wedge polygons on both the end moraines and the ground moraines extending up to the existing glacier front. Studies show that it requires hundreds if not thousands of years for such well-developed polygonal patterns to form.

Bottom melting rates have been calculated for a large number of sites on George VI Ice Shelf from measurements of its kinematic behaviour. No simple explanation for the melt-rate pattern was found in terms of ice-shelf parameters, assuming steady-state conditions. Values of apparent melt rates varied from 1 to 8 m a−1 of ice. Along different flow lines the melt rate would sometimes increase with distance from the grounding line and sometimes the melt rate would decrease with distance. Large melt rates were found both where ice flowed off Palmer Land and where the ice shelf butted against Alexander Island. Although oceanographic conditions probably control bottom melting rates the complex pattern with large spatial variation seems to indicate that some areas of ice shelf are changing in thickness.

Glaciers in the Russian High Arctic have undergone accelerated mass loss due to atmospheric and oceanic warming in the Barents–Kara Sea region. Most studies have concentrated on the western Barents–Kara sector, despite evidence of accelerating mass loss as far east as Severnaya Zemlya. However, long-term trends in glacier change on Severnaya Zemlya are largely unknown and this record may be complicated by surge-type glaciers. Here, we present a long-term assessment of glacier change (1965–2021) on Severnaya Zemlya and a new inventory of surge-type glaciers using declassified spy-satellite photography (KH-7/9 Hexagon) and optical satellite imagery (ASTER, Sentinel-2A, Landsat-4/5 TM and 8 OLI). Glacier area reduced from 17 053 km2 in 1965 to 16 275 in 2021 (−5%; mean: −18%, max: −100%), with areal shrinkage most pronounced at land-terminating glaciers on southern Severnaya Zemlya, where there is a recent (post-2010s) increase in summer atmospheric temperatures. We find that surging may be more widespread than previously thought, with three glaciers classified confirmed as surge-type, eight as likely to have surged and nine as possible, comprising 11% of Severnaya Zemlya's 190 glaciers (37% by area). Under continued warming, we anticipate accelerated retreat and increased likelihood of surging as basal thermal regimes shift.

From data taken on the British Trans-Arctic Expedition it is calculated that 9% of the Arctic Ocean surface between the North Pole and Spitsbergen was hummocked or ridged ice, 17% was unridged ice less than a year old, 73% was unridged old ice and 0.6% was ice-free. The mode of 250 thickness measurements taken through level areas of old floes along the entire traverse lies between 2.25 and 2.75 m. The mean end-of-winter thickness of the ice is calculated to be 4.6 m in the Pacific Gyral and 3.9 m in the Trans-Polar Drift Stream. From measurements of the percentage coverage and thickness of the various ice forms, it is calculated that the total annual ice accumulation in the Arctic Ocean is equivalent to a continuous layer of ice 1.1 m thick. 47% of this accumulation occurs in ice-free areas and under ice less than 1 year old. 20% of the total ice production is either directly or indirectly related to ridging or hummocking. An ice-ablation rate of 500 kg m−2 measured on a level area of a multi-year floe is compared with the rate on deformed and ponded ice. Greatest melting occurs on new hummocks and least on old smooth hummocks. The annual balance of ice older than 1 year but younger than multi-year ice is calculated from a knowledge of ice-drift patterns and the percentage coverage of first-year ice. The same calculations give a mean-maximum drift period of 5 years for ice in the Trans-Polar Drift Stream and 16 years in the Pacific Gyral. It is calculated that for the period February 1968 to May 1969 the annual ice export was 5 580 km3.

Seismic-reflection methods were used to determine the ice thickness and basal topography of Jakobshavns Isbræ, a large, fast-moving ice stream/outlet glacier in West Greenland. A method of data analysis was developed which involves the pointwise migration of data from a linear seismic array and a single explosive source; the method yields the depth, horizontal position and slope of the basal reflector. A deep U-shaped subglacial trough was found beneath the entire length of the well-defined ice stream. The trough is incised up to 1500 m into bedrock, and its base lies 1200–1500 m below sea level for at least 70 km inland. Center-line ice thickness along most of the channel is about 2500 m, or about 2.5 times that of the surrounding ice sheet. This prominent bedrock trough was not apparent in existing radio-echo-sounding data. Reflection coefficients indicate that much of the basal interface is probably underlain by compacted, non-deforming sediment. The large ice thickness, coupled with relatively steep surface slopes, leads to high basal shear stresses (200–300 k Pa) along the ice stream. The large shear stresses and lack of a deformable bed imply that internal deformation plays a dominant role in the dynamics of Jakobshavns Isbræ.

Two–dimensional hydrodynamic equations for laminar, viscous flow, and admitting a frictional slip-plane lower boundary are applied to the modeling of snow-avalanche impact on rigid wall structures. Predicted maximum pressures and pressures versus time are compared with published experimental results, and general correspondence is established. Impact pressure versus time is found to depend upon the shape of the avalanche leading edge, for which general information is lacking. Computer modeling of more complex structural configurations is feasible using the methodology reported.

Arctic ice shelves have declined over the past several decades, one of many indications of a rapidly changing cryosphere. Here we use a collection of off-nadir Landsat 8 images, a 1978 digital orthophotograph and photogrammetrically derived DEM, satellite altimetry and other data to examine the causes of an Arctic ice-shelf retreat in northernmost Greenland, the Hunt Fjord Ice Shelf (HFIS). HFIS has several distinct provenance regions comprised of glacier-derived ice and corrugated multi-decadal fast ice, with varying ice thicknesses (5–64 m). Available imagery shows little change in HFIS between 1978 and 2012, after which several midsummer calving events occurred (2012, 2016 and 2019) that reduced the HFIS by 42.5 km2 (~56%). Shelf area losses began as the number of surface melt days on the adjacent ice sheet more than doubled relative to the 1980s. Recent calving events also occurred during open-water periods at the ice-shelf front. Prior to mid-2012, there were no calving events during similar open-water periods. HFIS tributary glaciers have thinned by 3–20 m near their grounding zones, and may have accelerated since the 1980s, likely due to increased basal melting from contact with warm Atlantic Water.

A glacier-caused slide occurred in July 1975 within the Devastation Glacier basin, British Columbia. An estimated 2.5 × 106 m3 of glacier ice, representing about 9% of the total mass, was incorporated into displaced debris. Four human lives were lost.

Processes driving changes in sea-ice seasonality and sea-ice thickness were explored for a ‘warm-shelf’ region along the West Antarctic Peninsula using vertically coupled sea-ice-ocean thermodynamic simulations, with and without assimilated satellite sea-ice observations and moored ocean temperature observations. Simulations with assimilated sea-ice observations permitted investigation of surface [thermodynamic and dynamic (e.g., wind-driven)] processes affecting sea-ice thickness and seasonality. Assimilation of quasi-weekly variability in the depth and temperature of the deep warm pycnocline permitted examination of subsurface processes affecting sea-ice. Simulations using assimilated sea-ice observations (and implied motion) always produced greater surface heat fluxes and overall thinner sea ice. Assimilating seasonal and quasi-weekly variability in the depth and temperature of the pycnocline modified the start of the sea-ice season by −23 to +1 d, and also modified the sea ice thickness/seasonality to be thinner/shorter or thicker/longer at sub-seasonal and seasonal timescales, highlighting a mechanism where a shoaling pycnocline enhanced upward deep-water heat fluxes as transient surface-induced turbulence had a greater effect on a reduced mixed layer volume. The observed interplay of surface, subsurface, and sea-ice modulation of ocean-atmosphere heat transfer underscores the importance of representing the interaction between sea-ice concentration and upper ocean variability in climate projections.

Air temperature inversions, a situation in which atmospheric temperature increases with height, are key components of the Arctic planetary boundary layer. The present study investigates the spatial and temporal variations of temperature inversions over different surface types (rock, gravel, snow, ice) along the Mittivakkat valley (southeast Greenland). For this purpose, 113 vertical profiles with high spatio-temporal resolution of air temperature and relative humidity were collected with unoccupied aerial vehicles (UAVs) during a 13-day field campaign in summer 2019. Air temperature inversions were present in 83% of the profiles, of which 24% were surface-based inversions and 76% were elevated inversions. The proglacial area covered with bare rock and gravel induces surface heating and convection during the day and, through interaction with local circulation patterns, leads to the frequent formation of elevated inversions. In contrast, the glacier surface itself acts as a persistent cooling surface and leads to the formation of surface-based inversions. A low-level fog layer that forms under the inversion layer may be causing non-linear vertical ablation gradients on Mittivakkat Gletsjer. Furthermore, we demonstrate that atmospheric measurements using UAVs can better capture small-scale processes than other products like radiosonde or modeled reanalysis data.

We use a mathematical model to investigate the effect of basal topography and ice surface slope on transport and deposition of sediment within a water-filled subglacial channel. In our model, three zones of different behaviour occur. In the zone furthest upstream, variations in basal topography lead to sediment deposition under a wide range of conditions. In this first zone, even very small and gradually varying basal undulations (~5 m amplitude) can lead to the deposition of sediment within a modelled channel. Deposition is concentrated on the downstream gradient of subglacial ridges, and on the upstream gradient of subglacial troughs. The thickness and steepness of the ice sheet has a substantial impact on deposition rates, with shallow ice profiles strongly promoting both the magnitude and extent of sediment deposition. In a second zone, all sediment is transported downstream. Finally, a third zone close to the ice margin is characterised by high rates of sediment deposition. The existence of these zones has implications for esker formation and the dynamics of the subglacial environment.

Measurements have been made on Austerdalsbreen, Norway, to test the theory that wave ogives at the foot of an ice fall are formed by pressure. The pattern of deformation over the first three waves from the Odinsbre ice fall was studied by measuring the absolute and relative motions of 35 stakes during August 1956. There are large variations of longitudinal compression through the wave pattern, but no simple correlation exists between the compression and the crests; nor are the crests consistently rising with respect to the average surface. We conclude that the waves are not forming by pressure in this area. A detailed quantitative explanation of the observed pattern of deformation is found by improving the existing theory of glacier flow. The deformations taking place are explained as essentially independent of the presence of the waves; they are primarily due (a) to the compression that must occur in the ice as its slope decreases and it becomes thicker, and (b) to the bending and unbending to which the ice is subjected as it passes over a bed of changing curvature. Additional contributions to the deformation arise from the widening of the glacier channel, which is important at the immediate foot of the ice fall, and from the annual ablation.

The theory explains quantitatively the rotation and bending of the tunnel excavated in 1955; it was also used to make a prediction, which was successfully verified in a further experiment on the glacier in July and August 1957.

A pressure mechanism of wave formation can be reconciled with the observations if it is supposed that the compression in the lower two-thirds of the ice fall, which was measured photogrammetrically in August 1956, varies significantly with the seasons. However, the positions and amplitudes of the observed waves are fully accounted for by the combined plastic deformation and ablation mechanism recently described in another paper. There is therefore no longer any need for the pressure hypothesis as a primary cause of the waves.

The age hardening of artificially and naturally compacted snow has been investigated at the South Pole. Results show that the age-hardening process is greatly retarded at low temperatures. Artificially compacted samples of density 0.55 g./cm.3 attained a compressive strength of less than 3.0 kg./cm.2 after one year’s ageing at −49° C. Exposure to solar radiation accelerated the age hardening. Irradiated samples attained a strength of 6.0 kg./cm.2 after 100 hr., increasing to a virtual maximum of 8.0 kg./cm.2 at the end of 600 hr. Compressive strengths increased with decrease in snow-particle size and with increasing angularity of the particles. Below 3 m. the strength of naturally compacted snow was found to increase rapidly with increase in density. Naturally compacted snow of density 0.55 g./cm.3 possessed considerably greater strength than any of the age-hardened samples of artificially compacted snow of the same density. Thin-section studies show that age hardening can be correlated with the formation and growth of intergranular bonds, and that bond growth falls off rapidly with decreasing temperature. In view of the low strengths found in both naturally compacted snows near the surface and in artificially compacted snow at the South Pole, “cut-and-cover” under-snow camp construction may not prove too practical at the South Pole.

A procedure is described for testing cored ice specimens under uniaxial tension. The procedure incorporates an end-capping system described by Cole and others (1985), two carpeted Synthane end caps, and extensometers mounted directly on the ice. Use of the procedure has resulted in a successful test rate of 75%

Torsion creep tests were performed on glacier ice at temperatures above -12°C. The polycrystalline ice, when unloaded, exhibits creep recovery. The time-dependent recoverable component of deformation (or anelastic strain) ϵa was found to be adequately described by a relationship of the form: ϵa = Δτ log(r + αt)/h, where Δτ is the stress decrement, α a constant, and t the time. The anelastic modulus h defined for times t in excess of 3 h is always smaller than the dynamic elastic modulus. The movement of dislocations composing the sub-boundaries or in dislocation pile-ups may produce this important reversible deformation. The time-dependent recovery is explained in a similar way to the transient creep behaviour observed at low temperatures for metals. The small temperature dependence of creep recovery would arise from the existence of a distribution of internal stresses values.

The modeling of deformation in snow-pack of realistic geometry is reported. Experimentally measured deformation rates are compared with corresponding rates predicted from a finite-element model developed using local stratigraphy in the experimental area. Specifically designed deformation gages placed on an avalanche-sensitive snow slope at Berthoud Pass, Colorado, during the winter of 1972-73, were used to measure the experimental values. The finite-element code, based upon a plane-strain idealization, was used to model the steady-state flow of the slope. The program incorporates visco-elastic properties of the snow-pack, snow-layer stratigraphy, and the orthotropic properties of snow when both tension and compression components of stress act at a point. The degree of correlation between experimental and computed values of deformation rates is reported, and limitations of the analysis in conjunction with the type of experimental data taken are discussed.

Since their discovery, Antarctic subglacial lakes have become of great interest to the science community. It is hypothesized that they may hold unique forms of biological life and that they hold detailed sedimentary records of past climate change. According to the latest inventory, a total of 387 subglacial lakes have been identified in Antarctica (Wright and Siegert, 2011). However, exploration using scientific probes has yet to be performed. We propose a generic, formal approach to manage the operational risk of deploying probes during clean access to subglacial lake exploration. A representation of the entire probe deployment process is captured in a Markov chain. The transition from one state to the next depends on several factors, including reliability of components and processes. We use fault trees to quantify the probability of failure of the complex processes that must take place to facilitate the transition from one state to another. Therefore, the formal framework consists of integrating a Markov chain, fault trees, component and subsystem reliability data and expert judgment. To illustrate its application we describe how the approach can be used to address a series of what-if scenarios, using the intended Ellsworth Subglacial Lake probe deployment as a case study.

Robin (1967) and Budd (1968, unpublished) have succeeded in connecting the variations in surface slope of an ice sheet with variations in the gradient of the longitudinal strain-rate. This paper tries to improve the theoretical basis of their work. By choice of a suitable coordinate system and suitable redefinition of the variables, Budd’s formula for the basal shear stress is derived with a minimum of restrictive assumptions. The resulting formula, containing the gradient of a longitudinal stress, is thought to be of high accuracy for the two-dimensional flow of cold ice sheets, and is valid for slopes of any magnitude.

Temperate glacier ice is neither dry nor impermeable, as the standard theory of glacier sliding assumes. This fact leads to the already published concept of locally stress-controlled temperatures. Why the temperature is determined by the highest principal pressure, why the microscopic stress equals more or less the macroscopic one, and why water may flow in the capillary network even when water lenses at grain boundaries are freezing is explained. The new concept is applied to ice sliding on a hard bed having a sine profile, without cavitation. First, the stress field for a Newtonian viscous material and a vanishing roughness are used; next, an improved one, that takes into account the non-linear viscosity of ice and the finite amplitude of the micro-relief. It appears that water migrates from the stoss sides of the bumps to the lee sides within a bottom layer of thickness hw. Moreover, there is less ice melting at the sole on the former ones than ice accretion on the latter, a fact that yields a trend of ice accretion at the glacier sole. It is balanced by internal melting near the bed and water oozing at the interface from the soaked ice. Consequently, a thin layer of accreted regelation ice with a constant mean thickness hi should exist at the interface. Modelling realistically mountain glaciers, hw ∼ 20 cm and hi ∼ 3.5 cm are found.

Walder recently analyzed the initial instability of water flow in a sheet under a glacier that is produced by greater heat production in the flowing water and thus the larger rate of ice melting where a perturbation has increased the water sheet thickness. We have looked at the problem from the other presumed final state. We assume that instabilities have finally caused all the water to flow in channels (R-channels) at the bed. We investigated whether these channels can collect enough of the water that is produced by the geothermal heat and the heat sliding to remain in existence. When a basal shear stress is present, the distance out to which a channel can collect water is not that much greater than the channel radius itself. It is concluded that it is not likely for the channel to be able to collect appreciable amounts of water that is produced at the bed. Hence despite the indication that an instability initially might grow in a water sheet, this instability either cannot develop to the point when the water flow is primarily in channels or else the water flow alternates cyclically between a state primarily of sheet flow and a state of channel flow. It would appear that the channels that do exist under a glacier have their primary origin in the melt water from the upper surface that pours down moulins and thus is already channelized by the time it reaches the bed.