This paper reports detailed textural and gas measurements conducted in cold basal ice (–17°C from the margin of Taylor Glacier, an outlet glacier of the East Antarctic ice sheet. The analyzed samples were retrieved from a basal ice sequence excavated at the end of a subglacial tunnel dug near the glacier snout. The basal sequence exhibits two contrasting ice facies, defined as the englacial and stratified facies. On the one hand, analysis of ice crystal textures from the basal ice sequence provides evidence for localized ductile deformation, especially within the stratified facies where significant dynamic recrystallization was detected. On the other hand, high-resolution gas analyses reveal that strong changes in gas composition occurred at the structural interfaces of the stratified facies. These gas composition changes are typical of melting–refreezing processes but are not associated with any significant loss of gas volume. Given the specific subglacial thermal conditions at the margin of Taylor Glacier, we interpret this phenomenon as resulting from microscopic phase changes involving selective gas redistribution through the pre-melt phase. It is argued that such processes may play an important role in the post-genetic geochemical evolution of cold debris-laden ice and may be enhanced through intense strain conditions.

A numerical simulation of the Ross Ice Shelf, Antarctica, in which discharge from Ice Streams A-E is changed suddenly between extreme states, is used to investigate ice-shelf thickness and flow anomalies generated by ice-stream transience. At one extreme, ice-stream discharge rates specified as model boundary conditions are balanced individually with snow accumulation in the ice-stream catchment areas. At the other, discharge rates are fixed at current observed values which widely depart from mass balance. The simulated ice-shelf evolution between initial and final steady states suggests that ice-thickness and velocity fields adjust to new ice-stream conditions over a relatively short time span (approximately 500 years). In contrast, transitory geometrics of medial moraines and relict-crevasse bands persist over a longer time span (up to 2000 years). Contortions of medial moraines and relict-crevasse bands thus may provide a useful long-term history of past ice-stream activity. The past stoppage of Ice Stream C, for example, should be evident today in some medial moraine trajectories even if the stoppage occurred over 1000 years ago. Ice-shelf thickness fluctuations induced by ice-stream activity are generally restricted to the neighborhood of the grounding line. These fluctuations may constitute a trigger for ice-rise formation near ice-stream outlets.

We recorded electrical resistivity data at the base of four boreholes drilled through Haut Glacier d’Arolla, Switzerland. The data were acquired repetitively every hour over two diurnal hydrological cycles in the late melt season, separated by 10 days. Constrained three-dimensional (3-D) data inversion allowed reconstruction of hourly variations in bulk resistivity in the subglacial sediment layer. Inverted resistivity models reflect the establishment of channelized subglacial drainage in the study area between the two hydrological cycles, in agreement with previous work. Daily variations in bulk and water resistivity are in phase, and bulk resistivity amplitudes decrease away from the subglacial channel. Using selected electrical–hydraulic relationships, we estimate metre-scale changes in the hydraulic conductivity and porosity of the subglacial sediment layer, accounting for increasing clay content and decreasing median grain radius with distance from the channel. Hydraulic conductivity and porosity were respectively calculated to decrease from (6.4 ± 2.1) × 10–2ms–1 and 0.34 ± 0.01 at the channel to (3.3 ± 2.2) × 10–2ms–1 and 0.26 ± 0.01 at a distance of 5m from it. The hydraulic conductivity estimates are in agreement with previously inferred values, and the porosity estimates fall within the expected range for unlithified subglacial sediments. We conclude that collection and inversion of repeat 3-D subglacial resistivity data is feasible and has the capacity to generate multidimensional images of subglacial hydraulic processes and properties.

Coastal ice cores provide an opportunity to investigate regional climate and sea-ice variability in the past to complement hemispheric-scale climate reconstructions from ice-sheet-interior ice cores. Here we describe robust proxies of Baffin Bay temperature and sea-ice concentration from the coastal 2Barrel ice core collected in the Thule region of northwest Greenland. Over the 1990–2010 record, 2Barrel annually averaged methanesulfonic acid (MSA) concentrations are significantly correlated with May–June Baffin Bay sea-ice concentrations and summer temperatures. Higher MSA is observed during warmer years with less sea ice, indicative of enhanced primary productivity in Baffin Bay. Similarly, 2Barrel annually averaged deuterium excess (d-excess) values are significantly correlated with annual Baffin Bay sea-ice concentrations and summer and annual temperatures. Warm (cool) years with anomalously low (high) sea-ice concentration are associated with proportionally more (less) low-d-excess Baffin Bay moisture at the ice-core site. Multilinear regression models incorporating 2Barrel MSA, d-excess and snow accumulation account for 38–51% of the Baffin Bay sea-ice and temperature variance. The annual temperature model is significantly correlated with temperatures throughout most of Greenland and eastern Arctic Canada due to the strong influence of the North Atlantic Oscillation and Atlantic Multidecadal Oscillation.

Conventional avalanche forecasting is practiced as a mix of deterministic treatment for snow and weather parameters and inductive logic to reach actual forecast decisions. Inductive logic of the scientific method dominates, making frequent use of iteration and redundancy to minimize decision uncertainties. The mental processes involved are holistic rather than analytical. Elementary information theory can be used rationally to sort data categories for minimum entropy and optimize inductive reasoning. Recognizing these principles affords a chance to improve the practice and teaching of conventional forecasting techniques.

We derive the surface and basal radar reflectance and backscatter coefficients of the southern McMurdo Ice Shelf (SMIS) and part of the nearby Ross Ice Shelf (RIS), Antarctica, from radar statistical reconnaissance using a 60-MHZ airborne survey. The surface coefficients are further inverted in terms of snow density and roughness, providing a spatial distribution of the processes contributing to the surface boundary conditions. We disentangle the basal coefficients from surface transmission losses, and we provide the basal coherent content, an indicator of the boundary geometric disorder that is also self-corrected from englacial attenuation. The basal radar properties exhibit sharp gradients along specific iso-depths, suggesting an abrupt modification of the ice composition and geometric structure. We interpret this behavior as locations where the pressure-melting point is reached, outlining fields of freezing and melting ice. Basal steps are observed at both SMIS and RIS, suggesting a common geometric expression of widespread basal processes. This technique offers a simultaneous view of both the surface and basal boundary conditions to help investigate the ice-shelf stability, while its application to airborne data significantly improves coverage of the difficult-to-observe ice–ocean boundary. It also provides constraints on thermohaline circulation in ice shelves cavities, which are analogs for ice-covered ocean worlds.

Cryoconite is a dark-coloured granular sediment that contains biological and mineralogical components, and it plays a pivotal role in geochemistry, carbon cycling and glacier mass balance. In this work, we collected cryoconite samples from Laohugou Glacier No. 12 (LHG) on the north-eastern Tibetan Plateau during the summer of 2015 and measured the spectral albedo. To explore the impacts of this sediment on surface ablation, the ice melting differences between the cryoconite-free (removed) ice and the intact layers were compared. The results showed that the mean concentrations of black carbon (BC), organic carbon (OC) and total iron (Fe) in the LHG cryoconite were 1.28, 11.18 and 39.94 mg g−1, respectively. BC was found to play a stronger role in solar light adsorption than OC and free Fe. In addition, ice covered by cryoconite exhibited the lowest mean reflectance (i.e., <0.1). Compared with the cryoconite-free ice surface, cryoconite effectively absorbed solar energy and enhanced glacial melting at a rate of 2.27–3.28 cm d−1, and free Fe, BC and OC were estimated to contribute 1.01, 0.99 and 0.76 cm d−1, respectively. This study provides important insights for understanding the role of cryoconite in the glacier mass balance of the northern Tibetan Plateau.

In a temperate glacier ice the water content results from the past history of the ice, and the temperature adjusts itself to the appropriate value. Results obtained by the Laboratoire de Glaciologie are summarized. Several facts seem to prove some migration of the liquid phase and a coalescence of liquid inclusions. The enthalpy flux is calculated, but it cannot be proved that a spontaneous gathering together of liquid inclusions could occur. The water-rich, small grain-size layers should come from the lowering of the viscosity when the water content increases, as discovered by Duval. The form of Weertman’s sliding law without cavitation can be derived from simple dimensional considerations. Numerically, it cannot account for high sliding velocities, for which subglacial cavitation must be introduced. The cavities are infilled with more stagnant regelation ice than water. Two subglacial hydraulic regimes are distinguished: autonomous and interconnected. The effect which the salt content may have on cavitation is estimated.

The weakening and/or removal of floating ice shelves in Antarctica can induce inland ice flow acceleration. Numerical modelling suggests these processes will play an important role in Antarctica's future sea-level contribution, but our understanding of the mechanisms that lead to ice tongue/shelf collapse is incomplete and largely based on observations from the Antarctic Peninsula and West Antarctica. Here, we use remote sensing of structural glaciology and ice velocity from 2001 to 2020 and analyse potential ocean-climate forcings to identify mechanisms that triggered the rapid disintegration of ~2445 km2 of ice mélange and part of the Voyeykov Ice Shelf in Wilkes Land, East Antarctica between 27 March and 28 May 2007. Results show disaggregation was pre-conditioned by weakening of the ice tongue's structural integrity and was triggered by mélange removal driven by a regional atmospheric circulation anomaly and a less extensive latent-heat polynya. Disaggregation did not induce inland ice flow acceleration, but our observations highlight an important mechanism through which floating termini can be removed, whereby the break-out of mélange and multiyear landfast sea ice triggers disaggregation of a structurally-weak ice shelf. These observations highlight the need for numerical ice-sheet models to account for interactions between sea-ice, mélange and ice shelves.

An end of summer snowline (EOSS) photographic dataset for Aotearoa New Zealand contains over four decades of equilibrium line altitude (ELA) observations for more than 50 index glaciers. This dataset provides an opportunity to create a climatological ELA reference series that has several applications. Our work screened out EOSS sites that had low temporal coverage and also removed limited observations when the official survey did not take place. Snowline data from 41 of 50 glaciers in the EOSS dataset were retained and included in a normalised master snowline series that spans 1977–2020. Application of the regionally representative normalised master snowline series in monthly and seasonally resolved climate response function analyses showed consistently strong relationships with austral warm-season temperatures for land-based stations west of the Southern Alps and the central Tasman Sea. There is a trend towards higher regional snowlines since the 1990s that has been steepening in recent decades. If contemporary decadal normalised master snowline series trends are maintained, the average Southern Alps snowline elevation will be displaced at least 200 m higher than normal by the 2025–2034 decade. More frequent extremely high snowlines are expected to drive more extreme cumulative mass-balance losses that will reduce the glacierised area of Aotearoa New Zealand.

To predict the distribution of motion beneath glaciers on soft beds, the strength of the coupling between the ice and the bed and its variation with effective pressure must be known. A record of shear strain, acquired with a tiltmeter emplaced in till beneath Storglaciären, Sweden, indicates that fluctuations in water pressure cause variations in the local shear stress on the bed and that the bed deforms elastically in response to these variations. To estimate the shear stress from the elastic component of the total shear strain, the shear modulus of the till was measured in relaxation tests conducted in the laboratory with a ring-shear device. After accounting for the elastic compliance of the device, these tests yielded shear moduli of about 1000 and 1800 kPa at confining pressures of 85 and 280 kPa, respectively. These values are comparable to those of other granular materials undergoing recoverable shear strains of the same magnitude. The local shear stress on the till, calculated by applying the measured shear moduli to the tilt record, scales with Pe1.7, where Pe is the effective pressure. This relation implies that as Pe decreases at the ice/till interface, shear stresses on the till are reduced and concentrated elsewhere on the bed, perhaps where the till is absent or the glacier is frozen to the bed. When compared with the strength of the till determined from ring-shear tests, this relation also accounts for the lack of permanent deformation at depth in the bed during periods of low Pe and indicates that most basal motion was by sliding or ploughing.

A mass-balance model based on the energy balance at the snow or ice surface is formulated, with particular attention paid to processes affecting absorption of radiation. The model is applied to a small glacier, Glacier AX010 in the Nepalese Himalaya, and tests of its mass-balance sensitivity to input and climatic parameters are carried out. Calculated and observed area-averaged mass balances of the glacier during summer 1978 (June-September) show good agreement, namely -0.44 and -0.46 m w.e., respectively.

Results show the mass balance is strongly sensitive to snow or ice albedo, to the effects of screening by surrounding mountain walls, to areal variations in multiple reflection between clouds and the glacier surface, and to thin snow covers which alter the surface albedo. In tests of the sensitivity of the mass balance to seasonal values of climatic parameters, the mass balance is found to be strongly sensitive to summer air temperature and precipitation but only weakly sensitive to relative humidity.

Despite their high value and importance for various glaciological applications, detailed ice thickness measurements of alpine glaciers are still very limited. Knowledge of bedrock topography is essential for paleoglaciological studies. The Guliya ice cap located on the Tibetan Plateau is one of the highest and largest ice caps in mid-low latitude regions. A detailed ground-penetrating radar (GPR) survey was conducted on the Guliya ice cap in 2015 using 20 and 40 MHz frequency antennas. An empirical Bayesian kriging method was used for ice thickness interpolation and uncertainty assessment. GPR measurements revealed complex basal topography of the Guliya glacier with a maximum thickness of 371.12 ± 13 m. The internal reflections caused by changes in the dielectric properties were registered on the 40 MHz radargrams at the summit and were attributed to density variations. As a result of this fieldwork, one of the largest ice thickness datasets in High Mountain Asia was obtained. Guliya glacier elevation changes were assessed by differencing digital elevation models. The glacier gained mass from 2000 to 2015 with an average rate of 0.270 ± 0.11 m w.e. a−1 at the summit and 0.279 ± 0.11 m w.e. a−1 at the lower elevations.

Major ice loss has recently been observed along coastal outlet glaciers of the West Antarctic ice sheet, mainly due to increased melting below the ice shelves. However, the behavior of this marine ice sheet is poorly understood, leading to significant shortcomings in ice-sheet models attempting to predict future sea-level rise. The stability of a marine ice sheet is controlled by the dynamics at the grounding line, the boundary between the grounded ice stream and the floating ice shelf. One of the largest contributors to current sea-level rise is the fast-flowing Thwaites Glacier, which flows into the Amundsen Sea. Here we use an ice-stream/ice-shelf model and perform a number of experiments along a central flowline to analyze the sensitivity of its grounding line on centennial timescales. In the absence of width and buttressing effects, we find that the grounding line retreats by ˜300 km in 200 years from the present day (rate of 1.5 km a–1). With variable glacier width implemented in the model, flow convergence slows the retreat of Thwaites grounding line at 0.3–1.2 km a–1. The parameterization of ice-shelf buttressing according to different observed scenarios further reduces the glacier retreat and can even lead to a slight advance in the most buttressed case.

A flow-line model is presented for calculating the surface profile and the velocity, strain-rate, and stress fields in an ice sheet with given base-elevation profile, ice thickness at the dome (divide), flow-law parameters, mass-balance distribution, and convergence/divergence conditions along the flow line. The model, which is based on a “quasi-similarity” hypothesis as regards the horizontal velocity-depth profiles, accounts for changes along the flow line in the depth distributions of temperature, normal stress deviators, and possible enhanced flow of deep ice of Wisconsin origin. A curvilinear coordinate system is applied with horizontal axes along flow lines and surface-elevation contours, respectively. The flow equations are reduced to two differential equations, one for the surface-elevation profile, and the other for a profile function that determines the depth distributions of velocities and strain-rates. The two equations are coupled through a profile parameter that communicates the influence of velocity-profile changes to the surface-profile equation. It is shown that the variation along the flow line of this parameter should also be considered when deriving flow-law parameters from ice-sheet flow-line data. For a symmetric dome, explicit expressions are derived for the depth distributions of the vertical velocity, strain-rates, and stresses. The strain-rate profiles display an inflection about half-way down the ice sheet, and, in the case of isothermal ice, have surface values 2.2 times their depth-averaged values. The depth distribution of the vertical velocity indicates that a relatively thick layer of almost stagnant ice is present at the ice-sheet base below a dome.

Remote sensing of ice motion by tracking displacement of surface features is a valuable tool in glaciology. Efficient image feature-tracking programs, such as IMCORR, based on fast Fourier transform methods can produce misleading correlations if there are data gaps in either or both of the reference and search images. This is particularly problematic if the data gaps are regular in character, such as for Landsat7 images collected after the failure of the Scan Line Corrector (SLC-off images). We demonstrate that this situation can be alleviated by filling the data gaps with suitably chosen random data. We modified IMCORR to achieve this automatically (source code is included), but generic image-processing software could be used to modify inputs for other correlation packages. We test our method using images of Pine Island Ice Shelf, Antarctica, and document the acceleration of the velocity field for the floating extension of Pine Island Glacier over the decade 2001–11. We also combine our velocities with recent NASA Operation IceBridge ice thickness data from CReSIS to estimate the basal melt rates.

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.

Time-dependent ice-sheet modelling of a 176 000 km2 area in western Dronning Maud Land, East Antarctica, provided information on the ice sheet’s response to six climate-change scenarios. Another experiment was done to study changes in ice thickness, flow and basal temperature conditions between the present ice configuration and a simulated maximum palaeo-ice sheet. The input to the model included new datasets of bed and surface topography compiled for this study. The results of the six climate-change experiments, including a 0.5°C per century global-warming scenario, show that the ice sheet has a robust behaviour with respect to the different climate changes. The maximum change in ice volume was <5% of the initial volume in all climate runs. This is for only relatively short-term climate changes without major changes in global sea level, and also a simulated ice sheet without an ice shelf. The modelled long-term response time of the ice sheet, 20 kyr or more, indicates that the ice sheet may still be adjusting to the climate change that ended the Last Glacial Maximum. In the maximum palaeo-ice-sheet simulation, with a 5°C climate cooling and the grounding line located at the continental-shelf margin, ice thickness increased drastically downstream from the Heimefrontfjella mountain range but remained basically unaffected on the upstream polar plateau. Compared to present conditions, complex changes in basal temperatures were observed. The extent of areas with basal melting increased, for example in the deep trough of the Veststraumen ice stream. Areas at intermediate elevations in the landscape also experienced increased basal temperatures, with significant areas reaching the melting point. In contrast, high-altitude areas that today are clearly cold-based, such as around Heimefrontfjella and Vestfjella and the Högisen dome, experienced a 5–10°C decrease in basal temperatures in the palaeo-ice-sheet reconstruction. The results suggest that the alpine landscape within these mountain regions was formed by wet-based local glaciers and ice sheets prior to the late Cenozoic.

Modelling the hydrology of the Greenland ice sheet, including the filling and drainage of supraglacial lakes, requires melt inputs generated at high spatial and temporal resolution. Here we apply a high spatial (100 m) and temporal (1 hour) mass-balance model to a 450 km2 subset of the Paakitsoq region, West Greenland. The model is calibrated by adjusting the values for parameters of fresh snow density, threshold temperature for solid/liquid precipitation and elevation-dependent precipitation gradient to minimize the error between modelled output and surface height and albedo measurements from three Greenland Climate Network stations for the mass-balance years 2000/01 and 2004/05. Bestfit parameter values are consistent between the two years at 400 kg m-3, 2°C and +14% (100 m)-1, respectively. Model performance is evaluated, first, by comparing modelled snow and ice distribution with that derived from Landsat-7 ETM+ satellite imagery using normalized-difference snow index classification and supervised image thresholding; and second, by comparing modelled albedo with that retrieved from the MODIS sensor M0D10A1 product. Calculation of mass-balance components indicates that 6% of surface meltwater and rainwater refreezes in the snowpack and does not become runoff, such that refreezing accounts for 31% of the net accumulation.