The ice-cored Galena Creek Rock Glacier, Wyoming, USA, has been the subject of a number of studies that sought to determine the origin of its ice. We present new observations of the rock glacier's internal structure from ground-penetrating radar to constrain ice and debris distribution and accumulation. We imaged dipping reflectors in the center of the glacier that are weak and discontinuous, in contrast to strong reflectors toward the edge of the cirque beneath large debris-avalanche chutes. These reflectors form a network of concave-up, up-glacier dipping layers. We interpret these as englacial debris bands formed by large debris falls buried by subsequent ice and snow accumulation. They are discontinuous where ice outpaces debris accumulation, but with sufficient debris accumulation an interleaved pattern of ice and debris layers can form. We propose a model in which the ice in these interleaved layers is snowfall preserved by debris-facilitated accumulation. Large debris falls that occur in early spring bury sections of the snowpack, which are then preserved through summer and incorporated into the rock glacier body over time. This study highlights the importance of sequential accumulation of ice and debris for understanding the dynamics of rock glaciers and debris-covered glaciers.

The stress balance within an ice shelf is key to the resistance, or buttressing, it can provide and in part controls the rate of ice discharge from the upstream ice sheet. Unconfined ice shelves are widely assumed to provide no buttressing. However, theory and laboratory-scale analogue experiments have shown that unconfined, floating viscous flows generate buttressing via hoop stresses. Hoop stress results from the viscous resistance to spreading perpendicular to the flow direction in a diverging flow. We build on theoretical work to explore the controls on the magnitude of hoop-stress buttressing, deducing that buttressing increases with increasing effective viscosity and increasing divergence. We use an idealised model calibrated to unconfined sections of Antarctic ice shelves and find that many shelves have low effective viscosity, most likely due to extensive damage resulting from high extensional stresses. Therefore, they are unable to sustain the large hoop stresses required to resist flow. Some ice shelves that are surrounded by sea ice year-round have a greater effective viscosity and can provide buttressing, suggesting that sea ice reduces fracturing. However, we find that most unconfined ice shelves provide insignificant buttressing today, even when hoop stresses are considered in the stress balance.

Antarctica's ice shelves modulate the grounded ice flow, and weakening of ice shelves due to climate forcing will decrease their ‘buttressing’ effect, causing a response in the grounded ice. While the processes governing ice-shelf weakening are complex, uncertainties in the response of the grounded ice sheet are also difficult to assess. The Antarctic BUttressing Model Intercomparison Project (ABUMIP) compares ice-sheet model responses to decrease in buttressing by investigating the ‘end-member’ scenario of total and sustained loss of ice shelves. Although unrealistic, this scenario enables gauging the sensitivity of an ensemble of 15 ice-sheet models to a total loss of buttressing, hence exhibiting the full potential of marine ice-sheet instability. All models predict that this scenario leads to multi-metre (1–12 m) sea-level rise over 500 years from present day. West Antarctic ice sheet collapse alone leads to a 1.91–5.08 m sea-level rise due to the marine ice-sheet instability. Mass loss rates are a strong function of the sliding/friction law, with plastic laws cause a further destabilization of the Aurora and Wilkes Subglacial Basins, East Antarctica. Improvements to marine ice-sheet models have greatly reduced variability between modelled ice-sheet responses to extreme ice-shelf loss, e.g. compared to the SeaRISE assessments.

Glacier basal motion is responsible for the majority of ice flux on fast-flowing glaciers, enables rapid changes in glacier motion and provides the means by which glaciers shape alpine landscapes. In an effort to enhance our understanding of basal motion, we investigate the evolution of glacier velocity and ice-marginal lake stage on Kennicott Glacier, Alaska, during the spring–summer transition, a time when subglacial drainage is undergoing rapid change. A complicated record of > 50 m fill-and-drain sequences on a hydraulically-connected ice-marginal lake likely reflects the punctuated establishment of efficient subglacial drainage as the melt season begins. The rate of change of lake stage generally correlates with diurnal velocity maxima, both in timing and magnitude. At the seasonal scale, the up-glacier progression of enhanced summer basal motion promotes uniformity of daily glacier velocity fluctuations throughout the 10 km study reach, and results in diurnal velocity patterns suggesting increasingly efficient meltwater delivery to and drainage from the subglacial channel system. Our findings suggest the potential of using an ice-marginal lake as a proxy for subglacial water pressure, and show how widespread basal motion affects bulk glacier behavior.

Direct measurements of spatially distributed vertical strain within ice masses are scientifically valuable but challenging to acquire. We use manual marker tracking and automatic cross correlation between two repeat optical televiewer (OPTV) images of an ~100 m-long borehole at Derwael Ice Rise (DIR), Antarctica, to reconstruct discretised, vertical strain rate and velocity at millimetre resolution. The resulting profiles decay with depth, from −0.07 a−1 at the surface to ~−0.002 a−1 towards the base in strain and from −1.3 m a−1 at the surface to ~−0.5 m a−1 towards the base in velocity. Both profiles also show substantial local variability. Three coffee-can markers installed at different depths into adjacent boreholes record consistent strain rates and velocities, although averaged over longer depth ranges and subject to greater uncertainty. Measured strain-rate profiles generally compare closely with output from a 2-D ice-flow model, while the former additionally reveal substantial high-resolution variability. We conclude that repeat OPTV borehole logging represents an effective means of measuring distributed vertical strain at millimetre scale, revealing high-resolution variability along the uppermost ~100 m of DIR, Antarctica.

The response of glaciers to climate change has major implications for sea-level change and water resources around the globe. Large-scale glacier evolution models are used to project glacier runoff and mass loss, but are constrained by limited observations, which result in models being over-parameterized. Recent systematic geodetic mass-balance observations provide an opportunity to improve the calibration of glacier evolution models. In this study, we develop a calibration scheme for a glacier evolution model using a Bayesian inverse model and geodetic mass-balance observations, which enable us to quantify model parameter uncertainty. The Bayesian model is applied to each glacier in High Mountain Asia using Markov chain Monte Carlo methods. After 10,000 steps, the chains generate a sufficient number of independent samples to estimate the properties of the model parameters from the joint posterior distribution. Their spatial distribution shows a clear orographic effect indicating the resolution of climate data is too coarse to resolve temperature and precipitation at high altitudes. Given the glacier evolution model is over-parameterized, particular attention is given to identifiability and the need for future work to integrate additional observations in order to better constrain the plausible sets of model parameters.

Antarctic lakes with perennial ice covers provide the opportunity to investigate in-lake processes without direct atmospheric interaction, and to study their ice-cover sensitivity to climate conditions. In this study, a numerical model – driven by radiative, atmospheric and turbulent heat fluxes from the water body beneath the ice cover – was implemented to investigate the impact of climate change on the ice covers from two Antarctic lakes: west lobe of Lake Bonney (WLB) and Crooked Lake. Model results agreed well with measured ice thicknesses of both lakes (WLB – RMSE= 0.11 m over 16 years of data; Crooked Lake – RMSE= 0.07 m over 1 year of data), and had acceptable results with measured ablation data at WLB (RMSE= 0.28 m over 6 years). The differences between measured and modeled ablation occurred because the model does not consider interannual variability of the ice optical properties and seasonal changes of the lake's thermal structure. Results indicate that projected summer air temperatures will increase the ice-cover annual melting in WLB by 2050, but that the ice cover will remain perennial through the end of this century. Contrarily, at Crooked Lake the ice cover becomes ephemeral most likely due to the increase in air temperatures.

Glacial extent and mass balance are sensitive climate proxies providing solid information on past climatic conditions. However, series of annual mass-balance measurements of more than 60 years are scarce. To our knowledge, this is the first time the latewood density data (MXD) of the Swiss stone pine (Pinus cembra L.) have been used to reconstruct the summer mass balance (Bs) of an Alpine glacier. The MXD-based Bs well correlates with a Bs reconstruction based on the May to September temperature. Winter precipitation has been used as an independent proxy to infer the winter mass balance and to obtain an annual mass-balance (Bn) estimate dating back to the glaciological year 1811/12. The reconstructed MXD/precipitation-based Bn well correlates with the data both of the Careser and of other Alpine glaciers measured by the glaciological method. A number of critical issues should be considered in both proxies, including non-linear response of glacial mass balance to temperature, bedrock topography, ice thinning and fragmentation, MXD acquisition and standardization methods, and finally the ‘divergence problem’ responsible for the recently reduced sensitivity of the dendrochronological data. Nevertheless, our results highlight the possibility of performing MXD-based dendroglaciological reconstructions using this stable and reliable proxy.

We reanalyzed mass balance records at Taku and Lemon Creek Glaciers to better understand the relative roles of hypsometry, local climate and dynamics as mass balance drivers. Over the 1946–2018 period, the cumulative mass balances diverged. Tidewater Taku Glacier advanced and gained mass at an average rate of +0.25 ± 0.28 m w.e. a–1, contrasting with retreat and mass loss of −0.60 ± 0.15 m w.e. a−1 at land-terminating Lemon Creek Glacier. The uniform influence of regional climate is demonstrated by strong correlations among annual mass balance and climate data. Regional warming trends forced similar statistically significant decreases in surface mass balance after 1989: −0.83 m w.e. a–1 at Taku Glacier and −0.81 m w.e. a–1 at Lemon Creek Glacier. Divergence in cumulative mass balance arises from differences in glacier hypsometry and local climate. Since 2013 negative mass balance and glacier-wide thinning prevailed at Taku Glacier. These changes initiated terminus retreat, which could increase dramatically if calving begins. The future mass balance trajectory of Taku Glacier hinges on dynamics, likely ending the historic dichotomy between these glaciers.

Herein, we present results obtained from time-lapse imagery acquired by a digital single-lens reflex camera during 2014–18 used to monitor the Planpincieux Glacier on the Italian side of the Grandes Jorasses (Mont Blanc massif). We processed the images using image cross-correlation to measure the surface kinematics of the most active lobe of the glacier that presents a high velocity and frequent ice detachments. During the monitoring, we observed two or three periods of sharp acceleration per year that culminated with large break-offs followed by analogous decelerations. Overall, we registered more than 350 failures with a volume >100 m3, of which, 14 events had volumes larger than 10 000 m3. The study identified a monotonic positive relationship between the velocity and failure volume that may be used to estimate the volume of the collapses before an event. We identified the thresholds of velocity and acceleration that characterise the activation of the speed-up periods. The study allowed the characterisation of three different instability processes that lead to the break-off of ice chunks from the glacier terminus: (i) disaggregation, (ii) slab fracture and (iii) water tunnelling failure which can be differentiated based on the rheology, the volume involved and the trigger process.

Velocities within an icefall are typically the fastest within a glacier system and experience complex flow. The combination of convergent and fast flow, and steep slope generate a quickly changing and intensely fractured surface. This complicates velocity extraction from repeat satellite images, especially when common pattern matching procedures are used. In this study, we exploit the high temporal revisit of medium-resolution satellite images using a novel image matching technique, ensemble matching, making it possible to generate a high-resolution (30 m) velocity field from high-repeat image sequences despite challenging image conditions. We demonstrate this technique for the first time in the glaciology domain using repeat Sentinel-2 optical data over the famous Khumbu icefall, situated on the southern slopes of Mount Everest. Estimates of velocity go just over 1 m d−1, which is slower than summer velocities from noisy single pair image matching. This icefall is frequently crossed by high-altitude mountaineers who use a route confined by fixed ropes and ladders set out every season. The mountain climbers typically record their trajectory on their personal satellite navigation device. We use such volunteered geographic information to verify our velocity estimates, confirming our underestimation with ensemble matching. Besides unprecedented remotely sensed surface velocities over the icefall, we also note that the generated velocity field can aid with the planning of a safe passage through this icefall.

We present a newly developed 1-D numerical energy-balance and phase transition supraglacial lake model: GlacierLake. GlacierLake incorporates snowfall, in situ snow and ice melt, incoming water from the surrounding catchment, ice lid formation, basal freeze-up and thermal stratification. Snow cover and temperature are varied to test lake development through winter and the maximum lid thickness is recorded. Average wintertime temperatures of −2 to $-30^{\circ }{\rm C}$ and total snowfall of 0 to 3.45 m lead to a range of the maximum lid thickness from 1.2 to 2.8 m after ${\sim }250$ days, with snow cover exerting the dominant control. An initial ice temperature of $-15^{\circ }{\rm C}$ with simulated advection of cold ice from upstream results in 0.6 m of basal freeze-up. This suggests that lakes with water depths above 1.3 to 3.4 m (dependent on winter snowfall and temperature) upon lid formation will persist through winter. These buried lakes can provide a sizeable water store at the start of the melt season, expedite future lake formation and warm underlying ice even in winter.

Under the potential to reconstruct the past climatic and atmospheric conditions from a deep ice core in the coastal Antarctic site (Styx Glacier), an 8.84 m long firn core (73°50.975′ S, 163°41.640′ E; 1623 m a.s.l.) was initially studied to propose a reliable age scale for the local estimation of snow accumulation rate. The seasonal variations of δ18O, methanesulfonic acid (MSA) and non-sea-salt sulfate (nssSO42–) were used for the firn core dating and revealed 25 annual peaks (from 1990 to 2014) with volcanic sulfate signal. The observed declining trend in annual accumulation rate with a mean value of 146 ± 60 kg m–2 a–1 is likely to be linked to the changes of sea-ice extent in the Ross Sea region. Moreover, the temporal variation of the annual mean δ18O, an annual flux of MSA and nssSO42– also likely to be under the influence of ice-covered and open water area. This study suggests a potential to recover past changes in an oceanic environment and will be useful for the interpretation of the long ice core drilled at the same site.

On 29 April 2015, four beacons were deployed onto an ice island in the Strait of Belle Isle to record positional data. The ice island later broke up into many fragments, four of which were tracked by the beacons. The relative influences of wind drag, current drag, Coriolis force, sea surface height gradient and sea-ice force on the drift of the tracked ice island fragments were analyzed. Using atmospheric and oceanic model outputs, the sea-ice force was calculated as the residual of the fragments' net forces and the sum of all other forces. This was compared against the force obtained through ice concentration-dependent relationships when sea ice was present. The sea-ice forces calculated from the residual approach and concentration-dependent relationships were significant only when sea ice was present at medium-high concentrations in the vicinity of the ice island fragments. The forces from ocean currents and sea surface tilt contributed the most to the drift of the ice island fragments. Wind, however, played a minimal role in the total force governing the drift of the four ice island fragments, and Coriolis force was significant when the fragments were drifting at higher speeds.

We present a 3-year record of seasonal variations in ice speed and frontal ablation of 10 marine-terminating outlet glaciers along the coast of Prudhoe Land in northwestern Greenland. The glaciers showed seasonal speedup initiated between late May and early June, and terminated between late June and early July. Ice speed subsequently decreased from July to September. The timing of the speedup coincided with the onset of the air temperature rise to above freezing, suggesting an influence of meltwater availability on the glacier dynamics. No clear relationship was found between the speedup and the terminus position or the sea-ice/ice-mélange conditions. These results suggest that the meltwater input to the glacier bed triggered the summer speedup. The excess of summer speed (June–August) over the mean for the rest of the year accounted for 0.5–13% of the annual ice motion. Several glaciers showed seasonal frontal variations, i.e. retreat in summer and advance in winter. This was not due to ice-speed variations, but was driven by seasonal variations in frontal ablation. The results demonstrate the dominant effect of glacier surface melting on the seasonal speedup, and the importance of seasonal speed patterns on longer-term ice motion of marine-terminating outlet glaciers in Greenland.

Physically-based avalanche propagation models must still be locally calibrated to provide robust predictions, e.g. in long-term forecasting and subsequent risk assessment. Friction parameters cannot be measured directly and need to be estimated from observations. Rich and diverse data are now increasingly available from test-sites, but for measurements made along flow propagation, potential autocorrelation should be explicitly accounted for. To this aim, this work proposes a comprehensive Bayesian calibration and statistical model selection framework. As a proof of concept, the framework was applied to an avalanche sliding block model with the standard Voellmy friction law and high rate photogrammetric images. An avalanche released at the Lautaret test-site and a synthetic data set based on the avalanche are used to test the approach and to illustrate its benefits. Results demonstrate (1) the efficiency of the proposed calibration scheme, and (2) that including autocorrelation in the statistical modelling definitely improves the accuracy of both parameter estimation and velocity predictions. Our approach could be extended without loss of generality to the calibration of any avalanche dynamics model from any type of measurement stemming from the same avalanche flow.

We present the first direct comparison of turbulence conditions measured simultaneously over exposed ice and a 0.08 m thick supraglacial debris cover on Suldenferner, a small glacier in the Italian Alps. Surface roughness, sensible heat fluxes (~20–50 W m−2), latent heat fluxes (~2–10 W m−2), topology and scale of turbulence are similar over both glacier surface types during katabatic and synoptically disturbed conditions. Exceptions are sunny days when buoyant convection becomes significant over debris-covered ice (sensible heat flux ~ −100 W m−2; latent heat flux ~ −30 W m−2) and prevailing katabatic conditions are rapidly broken down even over this thin debris cover. The similarity in turbulent properties implies that both surface types can be treated the same in terms of boundary layer similarity theory. The differences in turbulence between the two surface types on this glacier are dominated by the radiative and thermal contrasts, thus during sunny days debris cover alters both the local surface turbulent energy fluxes and the glacier component of valley circulation. These variations under different flow conditions should be accounted for when distributing temperature fields for modeling applications over partially debris-covered glaciers.

The Himalayan glaciers contribute significantly to regional water resources. However, limited field observations restrict our understanding of glacier dynamics and behaviour. Here, we investigated the long-term in situ mass balance, meteorology, ice velocity and discharge of the Chhota Shigri Glacier. The mean annual glacier-wide mass balance was negative, −0.46 ± 0.40 m w.e. a−1 for the period 2002–2019 corresponding to a cumulative wastage of −7.87 m w.e. Winter mass balance was 1.15 m w.e. a−1 and summer mass balance was −1.35 m w.e. a−1 over 2009–2019. Surface ice velocity has decreased on average by 25–42% in the lower and middle ablation zone (below 4700 m a.s.l.) since 2003; however, no substantial change was observed at higher altitudes. The decrease in velocity suggests that the glacier is adjusting its flow in response to negative mass balance. The summer discharge begins to rise from May and peaks in July, with a contribution of 43%, followed by 38% and 19% in August and September, respectively. The discharge pattern closely follows the air temperature. The long-term observation on the ‘Chhota Shigri – a benchmark glacier’, shows a mass wastage which corresponds to the slowdown of the glacier in the past two decades.

Glaciers in the eastern Pamir have reportedly been gaining mass during recent decades, even though glaciers in most other regions in High Mountain Asia have been in recession. Questions still remain about whether the trend is strengthening or weakening, and how far the positive balances extend into the eastern Pamir. To address these gaps, we use three different digital elevation models to reconstruct glacier surface elevation changes over two periods (2000–09 and 2000–15/16). We characterize the eastern Pamir as a zone of transition from positive to negative mass balance with the boundary lying at the northern end of Kongur Tagh, and find that glaciers situated at higher elevations are those with the most positive balances. Most (67% of 55) glaciers displayed a net mass gain since the 21st century. This led to an increasing regional geodetic glacier mass balance from −0.06 ± 0.16 m w.e. a−1 in 2000–09 to 0.06 ± 0.04 m w.e. a−1 in 2000–15/16. Surge-type glaciers, which are prevalent in the eastern Pamir, showed fluctuations in mass balance on an individual scale during and after surges, but no statistical difference compared to non-surge-type glaciers when aggregated across the region.

Systematic experiments reveal that the flexural strength of freshwater S2 columnar-grained ice loaded normal to the columns increases upon cyclic loading. Specifically, over the range of stress amplitudes 0.1–2.6 MPa the flexural strength increases linearly with increasing stress amplitude. The experiments were conducted upon both reversed and non-reversed cyclic loading over ranges of frequencies from 0.03 to 2 Hz and temperatures from −25 to −3°C. Strengthening can also be imparted through bending-induced creep. The fundamental requirement for strengthening is that the surface that undergoes maximum tensile stress during failure must have been pre-stressed in tension. Flexural strength is governed by crack nucleation. We suggest that the process is resisted by an internal back-stress that opposes the applied stress and builds up through either crystal dislocations piling up or grain boundaries sliding.