An experimental study concerning the mass and heat balance of snowpatches was performed during the Finnish Antarctic Research Programme (FINNARP) 2004 and 2010 summer expeditions to Basen nunatak (73°03′ S, 13°25′ W). Data were collected from a snow stake line, snow pits and automated weather and snow recording systems. One 100 m perennial snowpatch and several smaller seasonal patches (<10 m) were monitored. Snow thickness decreased by 4.0–6.3 mm d−1 due to sublimation, compression and, close to lateral boundaries, meltwater runoff. The vertical mass loss was 1–2 mm snow water equivalent (SWE) d−1 and the lateral decay was −10 cm d−1. The net radiation was 20.2 W m−2 and the mean latent heat flux was −15.5 W m−2 .The mean surface energy flux was 4.9 W m −2 and the heat loss to the ground was 1.5 W m−2. Thin snow decayed faster due to surface thermomechanical erosion and melt from the bottom where the soil was heated by the solar radiation. Between the summers of 2004 and 2010, the thickness of the perennial snowpatch decreased by 230 mm.

Glacier dynamics are inextricably linked to the basal conditions of glaciers. Seismic reflection methods can image the glacier bed under certain conditions. However, where a seismically thin layer of material is present at the bed, traditional analyses may fail to fully characterize bed properties. We use a targeted full-waveform inversion algorithm to quantify the basal-layer parameters of a mountain glacier: thickness (d), P-wave velocity (α) and density (ρ). We simultaneously invert for the seismic quality factor (Q) of the bulk glacier ice. The inversion seeks to minimize the difference between the data and a one-dimensional reflectivity algorithm using a gradient-based search with starting values initialized from a Monte Carlo scheme. We test the inversion algorithm on four basal layer synthetic data traces with 5% added Gaussian noise. The inversion retrieved thin-layer parameters within 10% of synthetic test parameters with the exception of seismic Q. For the seismic dataset from Bench Glacier, Alaska, USA, inversion results indicate a thin basal layer of debris-rich ice within the study area having mean velocity 4000 ± 700 m s–1, density 1900 ± 200 kg m–3 and thickness 6 ± 1.5 m.

Flexural strength of iceberg and glacier ice was determined from Four-point beam-bending experiments. A large quantity of glacial ice was collected from four icebergs and one glacier, and a detailed ice-characterization program was performed on samples from the five sources. Beam-bending experiments were conducted at four temperatures in the range −1 ° to-16 °C and at strain rates of 10 −3 and 10−5 s−1. The flexural strength was found to increase with increasing strain rate (based on extreme fibre strain and decreasing temperature. The data suggest than air-bubble inclusions play an important role in determining the flexural strength of glacial ice and this can explain the significant differences in mean strength of the ice from the five sources. At a strain rate of 10−3 s−1 and temperature of —11 °C, the flexural strength was found to increase as the number of bubbles per unit volume increased. Reduction of crack-initiating stresses at grain boundaries by “softening” of grains due to intragranular air-bubble inclusions is thought to be the mechanism.

A theoretical density distribution of the Dumber of bonds per grain as a function of the mean number of bonds per grain is derived from the assumption of randomness and isotropy of grain and grain-bond location and orientation. The knowledge of the theoretical density distribution allows the determination of the effective distribution parameters from section planes or thin sections. The concept of the fundamental unit to describe the strength of snow is introduced. Structural parameters developed on the basis of the fundamental units show improved correlations with the tensile strength.

Further to systematic experiments on the flexural strength of laboratory-grown, fresh water ice loaded cyclically, this paper describes results from new experiments of the same kind on lake ice harvested in Svalbard. The experiments were conducted at −12 °C, 0.1 Hz frequency and outer-fiber stress in the range from ~ 0.1 to ~ 0.7 MPa. The results suggest that the flexural strength increases linearly with stress amplitude, similar to the behavior of laboratory-grown ice.

There are several cases in which large overthrusts and sub-horizontal faults appear to have occurred in temperate or cold glaciers. As a contribution to solving the problem of their origin, the stress field when there is an abrupt change in the bottom boundary conditions is determined, assuming ice to be isotropic, third-power-law viscous. Deviatoric stresses vary with the distance r to the singularity as r−1/4, and strain rates as r−3/4. They are computed numerically to a multiplicative factor, which is determined by the conditions at a large distance, but not computed here. Although the apparent viscosity varies as r1/2 times a function of the polar angle, the stress field around the singularity is not essentially different from that obtained assuming a constant viscosity. Some considerations on the apparition of faults follow, but at the scale of the microrelief the adopted model becomes oversimplified.

A diffusion equation is obtained that describes the mechanical dispersion of a dilute mixture of solid particles within an ice matrix that is undergoing deformation. It is shown that within the limits of time intervals and strain-rates appropriate to the movement of glaciers and ice sheets, the dispersal distance usually is no larger than a distance about one order of magnitude greater than the size of the particles themselves.

Two air-permeability devices were developed and measurements were made to study the relationship between air permeability and structure of different snow types. Permeabilities varied both with position and direction in the snow-pack, and changed with time as metamorphic changes occurred. A marked increase of air permeability was noted as faceted crystals grew due to radiation recrystallization. In other observations, as a snow-pack densified and crystals became more rounded (destructive metamorphism) the air permeability decreased. Measurements were made of air permeability of snows that had been subjected to a creeping tensile strain, and initial tests indicated that the changes in air permeability due to strain, even to rupture, were significant, but small compared with the intrinsic variability of snow.

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.

Mass loss at the Greenland Ice Sheet is influenced by atmospheric processes controlling its surface mass balance, and by submarine melt and calving where glaciers terminate in fjords. There, an ice mélange – a composite matrix of calved ice bergs and sea ice – may provide a buttressing force on a glacier terminus and control terminus dynamics. Kangerlussuaq Glacier is a major outlet of the Greenland Ice Sheet, for which recent major retreat events in 2004/2005 and 2016–2018 coincided with the absence of an ice mélange in Kangerlussuaq Fjord. To better understand the response of Kangerlussuaq Glacier to climatic and oceanic drivers, a 2D flowline model is employed. Results indicate that an ice mélange buttressing force exerts a major control on calving frequency and rapid retreat. When an ice mélange forms in Kangerlussuaq Fjord, it provides stabilising forces and conditions favourable for winter terminus re-advance. When it fails to form during consecutive years, model results indicate that Kangerlussuaq Glacier is primed to retreat into the large overdeepenings in Kangerlussuaq Fjord, and to terminus positions more than 30 km farther inland, implying that excessive mass loss from Kangerlussuaq Glacier by the year 2065 cannot be excluded.

Jakobshavns Isbræ (69 °10′ N, 49 °59′ W) drains about 6.5% of the Greenland ice sheet and is the fastest ice stream known. The Jakobshavns Isbræ basin of about 10 000 km2 was mapped photogrammetrically from four sets of aerial photography, two taken in July 1985 and two in July 1986. Positions and elevations of several hundred natural features on the ice surface were determined for each epoch by photogrammetric block aerial triangulation, and surface velocity vectors were computed from the positions. The two flights in 1985 yielded the best results and provided most common points (716) for velocity determinations and are therefore used in the modeling studies. The data from these irregularly spaced points were used to calculate ice elevations and velocity vectors at uniformly spaced grid points 3 km apart by interpolation. The field of surface strain rates was then calculated from these gridded data and used to compute the field of surface deviatoric stresses, using the flow law of ice, for rectilinear coordinates, X, Y pointing eastward and northward, and curvilinear coordinates. L, Τ pointing longitudinally and transversely to the changing ice-flow direction, Ice-surface elevations and slopes were then used to calculate ice thicknesses and the fraction of the ice velocity due to basal sliding. Our calculated ice thicknesses are in fair agreement with an ice-thickness map based on seismic sounding and supplied to us by K. Echelmeyer. Ice thicknesses were subtracted from measured ice-surface elevations to map bed topography. Our calculation shows that basal sliding is significant only in the 10–15 km before Jakobshavns Isbræ becomes afloat in Jakobshavns Isfjord.

Using an anti-icing coating to prevent ice accretion on the drill surface is a feasible solution to address the drilling difficulties in warm ice. In this study, four types of commercially available hydrophobic coating materials were tested to evaluate their water repellency and anti-icing properties, namely, a mixture of silica and fluorocarbon resin with polytrifluoroethylene, modified Teflon, silica-based emulsion and an acrylic-based copolymer. Their water contact angles are ~107°, 101°, 114° and 95°, respectively. All these hydrophobic coatings can significantly reduce the strength of the ice adhesion within a temperature range of −10 to −30°C on a planar or curved surface. The coating of an acrylic-based copolymer, in particular, can reduce the average tensile strength and the shear strength of the ice adhesion by 87.08 and 97.11% on planar surfaces at −30°C, and by 98.06 and 96.15% on a curved surface, respectively. The main challenge in the practical application of these coatings is their durability. An acrylic-based copolymer coating will lose its water repellency performance after 140 cycles of abrasion. The shear strength of ice adhered on curved surfaces coated with this material will approach that achieved on uncoated surfaces after 11 cycles of icing and de-icing tests.

Evolution of sedimentological characteristics for traction debris is investigated along a transport path down-glacier from a subglacial point source at the northern rim of the ice cap MỲrdalsjökull, Iceland. Losses and gains from this path are discussed. Close to mature characteristics as to particle size, shape, roundness, and fabric are developed after only short basal transport, about 500 m, while further transport only results in minor changes. This development is interpreted as the result of selective subglacial lodgement of large particles and/or particles with longitudinal α-axis orientation, while particles of median size and/or oblique orientation tend to stay in basal transport. Implications for basal sliding and subglacial sedimentation patterns are discussed.

Glacier de Saint-Sorlin, French Alps, left terminal moraines at 1.3, 2.9 and 3.7 km ahead of the present terminus. According to proxy data and to historical maps, these were formed in the 19th, 18th and 17th centuries, respectively. A plateau at 2700–2625 m was then surrounded by ice but never became an accumulation area. This fact shows that the equilibrium-line altitude (ELA) on the glacier never dropped below 2300 m. The following simple models apply sufficiently to yield reliable estimations of past ELA: (1) a uniform and constant vertical gradient of the mass balance, down to the terminus; and (2) a plane bed, with a slope of 8.5° and a uniform width. Then in a steady situation the accumulation–area ratio is 1/2. Compared to the mean for 1956–72, at the onset of the Little Ice Age the balances were higher by 3.75 m ice a−1, and the ELA was 400 m lower. Correlations between 1956–72 balances and meteorological data suggest that during the melting season the 0°C isotherm was about 800 m lower, while the winter precipitation at low altitudes did not change. These correlations may have been different in the past, but an equal lowering of the ELA and of the 0°C isotherm, as assumed by several authors, seems excluded.

Crary Ice Rise formed after the Ross Ice Shelf re-grounded ~1 kyr BP. We present new ice-penetrating radar data from two systems operating at center frequencies of 7 and 750 MHz that confirm the ice rise is composed of a former ice shelf buried by subsequent accumulation. Stacks of englacial diffraction hyperbolas are present almost everywhere across the central ice rise and extend up to ~350 m above the bed. In many cases, bed reflections beneath the diffraction hyperbolas are obscured for distances up to 1 km. Waveform modeling indicates that the diffraction hyperbolas are likely caused by marine ice deposits in former basal crevasses and rifts. The in-filling of rifts and basal crevasses may have strengthened the connection between the ice rise and the surrounding ice shelf, which could have influenced local and regional ice dynamics. Three internal reflection horizons mark the upper limit of disturbed ice and diffraction hyperbolas in different sections of the ice rise, indicating at least three stages of flow stabilization across the ice rise. A surface lineation visible in MODIS imagery corresponds spatially to deepening and strong deformation of these layers, consistent with the characteristics of former grounding lines observed elsewhere in Antarctica.

Liquid water can exist at temperatures well below freezing beneath glaciers and ice sheets, where subglacial water systems, fresh and saline, have been shown to host unique microbial ecosystems. Geophysical techniques sensitive to fluid-content contrasts, e.g. electromagnetics, can characterize subglacial water and its salinity. Here, we assess the ground-based transient electromagnetic (TEM) method for deriving the resistivity and salinity of subglacial water. We adapt an existing open-source Bayesian inversion algorithm, which uses independent depth constraints, to output posterior distributions of resistivity and pore fluid salinity with depth. A variety of synthetic models, including a thin (5 m), conductive (0.16 Ωm), hypersaline (147 psu) subglacial lake, are used to evaluate the TEM method for imaging under 800 m-thick ice. The study demonstrates that TEM methods can resolve conductive, saline bodies accurately using external depth constraints, for example, from radar or seismic data. The depth resolution of TEM can be limited beneath deep (>800 m), thick (>50 m) conductive, water bodies and additional constraints from passive electromagnetic (EM) methods could be used to reduce ambiguities in the TEM results. Subsequently, non-invasive active and passive EM methods could provide profound insights into remote aqueous systems under glaciers and ice sheets.

Under global warming, many glaciers worldwide are receding. However, recent studies have suggested the extension of the Karakoram Anomaly, a region of anomalous glacier mass gain, into the western Kunlun and eastern Pamir mountains. However, the eastern limit of this anomaly in the Kunlun Mountains is unclear. This study, using changes in glacier area and surface elevation, estimates the eastern limit of the Kunlun-Pamir-Karakoram anomaly at ~85°E. Over the past 50 years, glaciers west of 85°E in the Kunlun Mountains decreased in area from 8401 to 7945 km2 at a rate of −0.12 ± 0.07% a−1, showed a reduction in the rate of retreat through time and have recently gained mass, with surface elevation changes of 0.15 ± 0.35 m a−1 over the period of 2000–2013. Glaciers east of 85°E have experienced greater rates of area change (−61 ± 12 km2 and −0.43 ± 0.13% a−1) over the past 50 years, accelerated area loss in recent years and elevation change rate of −0.51 ± 0.18 m a−1 between 2000 and 2013. These patterns of elevation and area change are consistent with regional increases in summer temperature in the eastern Kunlun Mountains and slight cooling in the western Kunlun Mountains.

We use ice flow modelling to simulate the englacial stratigraphy of Blåskimen Island, an ice rise in Dronning Maud Land and elucidate the evolution of this data-sparse region. We apply a thermo-mechanically coupled Elmer/Ice model to a profile along flowlines and through the ice-rise summit, where surface mass balance (SMB), flow velocity and ice stratigraphy were recently measured. We conclude that: (i) the ice rise is presently thickening at a rate of 0.5~0.6 m ice equivalent per year (mieq a−1), which is twice an earlier estimate using the field data and the input–output method; (ii) present thickening started 20–40 years in the past, before which the ice rise was in a steady state; (iii) SMB contrast between the upwind and downwind slopes was stronger than the present value by ~23% (or 0.15 mieq a−1) prior to ~1100 years ago. Since then, this contrast has been decreasing overall. We surmise that these SMB changes are likely a result of synoptic-scale atmospheric changes, rather than local atmospheric changes controlled by local ice topography. Our technique effectively assimilates geophysical data, avoiding the complexity of ice flow beneath the ice divide. Thus, it could be applied to other ice rises to elucidate the recent glacial retreat.

The evolution of four dynamic radar glacier zones at the surface of an alpine icefield in British Columbia is monitored using a time series of 35 First European Remote Sensing Satellite (ERS-1) synthetic aperture radar (SAR) images acquired from 1992 to 1994. These zones result from changing wetness and textural properties, and appear to represent: (1) cold snow with no liquid water present; (2) an initial melt front with an upper boundary near the elevation of the 0° isotherm; (3) metamorphosed, rapidly melting first-year snow with a rough or pitted surface; and (4) bare ice. This interpretation is aided by temperature and runoff data, air photographs and field measurements of snowpack properties acquired with two ERS-1 SAR scenes, ice-surface elevations derived from 1:50 000 topographic maps and simulations of radar backscatter from a geometric optics model of surface scattering. Meltwater production is affected by the development of zones (2), (3) and (4), which form, migrate up-elevation and disappear each year between April and September.