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.

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.

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.

Negribreen, a tidewater glacier located in central eastern Svalbard, began actively surging after it experienced an initial collapse in summer 2016. The surge resulted in horizontal surface velocities of more than 25 m d−1, making it one of the fastest-flowing glaciers in the archipelago. The last surge of Negribreen likely occurred in the 1930s, but due to a long quiescent phase, investigations of this glacier have been limited. As Negribreen is part of the Negribreen Glacier System, one of the largest glacier systems in Svalbard, investigating its current surge event provides important information on surge behaviour among tidewater glaciers within the region. Here, we demonstrate the surge development and discuss triggering mechanisms using time series of digital elevation models (1969–2018), surface velocities (1995–2018), crevasse patterns and glacier extents from various data sources. We find that the active surge results from a four-stage process. Stage 1 (quiescent phase) involves a long-term, gradual geometry change due to high subglacial friction towards the terminus. These changes allow the onset of Stage 2, an accelerating frontal destabilization, which ultimately results in the collapse (Stage 3) and active surge (Stage 4).

Streaming ice accounts for a major fraction of global ice flux, yet we cannot yet fully explain the dominant controls on its kinematics. In this contribution, we use an anisotropic full-Stokes thermomechanical flow solver to characterize how mechanical anisotropy and temperature distribution affect ice flux. For the ice stream and glacier geometries we explored, we found that the ice flux increases 1–3% per °C temperature increase in the margin. Glaciers and ice streams with crystallographic fabric oriented approximately normal to the shear plane increase by comparable amounts: an otherwise isotropic ice stream containing a concentrated transverse single maximum fabric in the margin flows 15% faster than the reference case. Fabric and temperature variations independently impact ice flux, with slightly nonlinear interactions. We find that realistic variations in temperature and crystallographic fabric both affect ice flux to similar degrees, with the exact effect a function of the local fabric and temperature distributions. Given this sensitivity, direct field-based measurements and models incorporating additional factors, such as water content and temporal evolution, are essential for explaining and predicting streaming ice dynamics.

Computing predictions of future sea level that include well-defined uncertainty bounds requires models that are capable of robustly simulating the evolution of ice sheets and glaciers. Ice flow behaviour is known to be sensitive to the location and geometry of dynamic ice boundaries such as the grounding line (GRL), terminus position and ice surface elevation, so that any such model should track these interfaces with a high degree of accuracy. To address this challenge, we implement a numerical approach that uses the level-set method (LSM) that accurately models the evolution of the ice–air and ice–water interface as well as capturing topological changes in ice-sheet geometry. This approach is evaluated by comparing simulations of grounded and marine-terminating ice sheets to various analytical and numerical benchmark solutions. A particular advantage of the LSM is its ability to explicitly track the moving margin and GRL while using a fixed grid finite-difference scheme. Our results demonstrate that the LSM is an accurate and robust approach for tracking the ice surface interface and terminus for advancing and retreating ice sheets, including the transient marine ice-sheet interface and GRL positions.

Deglaciation of high mountain rockwalls alters slope stability as rockwalls become more sensitive to modifications in environmental factors (e.g. seasonal temperature variations). In the past decades, increasing efforts focused on studying deglaciated Alpine rockwalls. Yet, currently deglaciating rockfaces remain unstudied. Here, we quantify surface area variations of massive ice bodies lying on high mountain rockwalls (ice aprons) in the French sector of the Mont Blanc massif between the end of the Little Ice Age (LIA) and 2018. Surface area estimates are computed from terrestrial and aerial oblique photographs via photogrammetry. This technique allows using photographs taken without scientific intent, and to tap into diverse historical or recent photographic catalogs. We derive an ice apron surface area model from precipitation records and the positive degree-days. The studied ice aprons shrank from 1854 to the 1950s, before expanding until the end of the 1990s. The beginning of the 21st century shows a decrease in surface area, leading to the complete melt of one of the studied ice aprons in 2017. Observed variations correlate with modeled surface area, suggesting strong sensitivity of ice aprons to changes in climatic variables. By studying site-specific correlations, we explore the importance of local drivers over the balance of ice aprons.

Glacial environments exhibit temporally variable microseismicity. To investigate how microseismicity influences event detection, we implement two noise-adaptive digital power detectors to process seismic data from Taylor Glacier, Antarctica. We add scaled icequake waveforms to the original data stream, run detectors on the hybrid data stream to estimate reliable detection magnitudes and compare analytical magnitudes predicted from an ice crack source model. We find that detection capability is influenced by environmental microseismicity for seismic events with source size comparable to thermal penetration depths. When event counts and minimum detectable event sizes change in the same direction (i.e. increase in event counts and minimum detectable event size), we interpret measured seismicity changes as ‘true’ seismicity changes rather than as changes in detection. Generally, one detector (two degree of freedom (2dof)) outperforms the other: it identifies more events, a more prominent summertime diurnal signal and maintains a higher detection capability. We conclude that real physical processes are responsible for the summertime diurnal inter-detector difference. One detector (3dof) identifies this process as environmental microseismicity; the other detector (2dof) identifies it as elevated waveform activity. Our analysis provides an example for minimizing detection biases and estimating source sizes when interpreting temporal seismicity patterns to better infer glacial seismogenic processes.

The accuracy of sea-ice motion products provided by the National Snow and Ice Data Center (NSIDC) and the Ocean and Sea Ice Satellite Application Facility (OSI-SAF) was validated with data collected by ice drifters that were deployed in the western Arctic Ocean in 2014 and 2016. Data from both NSIDC and OSI-SAF products exhibited statistically significant (p < 0.001) correlation with drifter data. The OSI-SAF product tended to overestimate ice speed, while underestimation was demonstrated for the NSIDC product, especially for the melt season and the marginal ice zone. Monthly Lagrangian trajectories of ice floes were reconstructed using the products. Larger spatial variability in the deviation between NSIDC and drifter trajectories was observed than that of OSI-SAF, and seasonal variability in the deviation for NSIDC was observed. Furthermore, trajectories reconstructed using the NSIDC product were sensitive to variations in sea-ice concentration. The feasibility of using remote-sensing products to characterize sea-ice deformation was assessed by evaluating the distance between two arbitrary positions as estimated by the products. Compared with the OSI-SAF product, relative errors are lower (<11.6%), and spatial-temporal resolutions are higher in the NSIDC product, which makes it more suitable for estimating sea-ice deformation.

Drumlins form at the ice/bed interface through subglacial processes that are not directly observable. The internal stratigraphy of drumlins provides insight into how they developed and associated subglacial processes, but traditional stratigraphic logging techniques are limited to natural exposures and excavations. Using ground-penetrating radar, we imaged the internal stratigraphy of seven drumlins from a recently exposed drumlin field in the forefield of Múlajökull, Iceland. Data were collected with 100 and 200 MHz antennas with maximum resolvable depths of 8 and 4 m, respectively. Longitudinal echograms contained coherent down-ice dipping reflectors over the lengths of the drumlins. Near the drumlin heads (i.e., stoss sides), down-glacier dipping beds lie at high angles to the surface, whereas on the lee sides, the down-glacier dipping beds lie at low angles, or conform, to drumlin surfaces. Transverse echograms exhibited unconformities along the flanks of drumlin heads and conformable bedding across the lee side widths of the drumlins. These observations were ground-truthed with stratigraphic logs from a subset of drumlins and good agreement was found. The stratigraphic patterns support previous conclusions that drumlins at Múlajökull formed on a deformable bed through both depositional and erosional processes which may alternate between its surge and quiescent phases.

Climate-induced cryospheric changes can have a significant impact on the downstream water availability. In this study, the Open Global Glacier Model (OGGM) and the Glacio-hydrological Degree-day Model (GDM) are integrated to project the response of cryospheric and hydrological systems to climate change until 2100. The study area comprises six sub-basins of glacierized Koshi River basin covering Nepalese and Chinese territories. The output from OGGM is provided as input to GDM along with the spatial and hydro-meteorological data. The average glacier area change in all the sub-basins from 2021 to 2100 is estimated as 65 and 85% decrease and the average glacier volume change is estimated as 76 and 86% decrease for RCP 4.5 and 8.5 scenarios, respectively. The future simulated discharge shows an increasing trend in pre-monsoon and monsoon seasons and a decreasing trend in post-monsoon and winter seasons after 2060 in all the sub-basins, which can lead to wetter wet seasons and drier dry seasons in the far future. A shift in peak flow is observed from August to July in most of the sub-basins. The coupled modelling technique used in this study can largely improve our understanding of glacio-hydrological dynamics in the Himalayan region.

Using an ensemble of close- and long-range remote sensing, lake bathymetry and regional meteorological data, we present a detailed assessment of the geometric changes of El Morado Glacier in the Central Andes of Chile and its adjacent proglacial lake between 1932 and 2019. Overall, the results revealed a period of marked glacier down wasting, with a mean geodetic glacier mass balance of −0.39 ± 0.15 m w.e.a−1 observed for the entire glacier between 1955 and 2015 with an area loss of 40% between 1955 and 2019. We estimate an ice elevation change of −1.00 ± 0.17 m a−1 for the glacier tongue between 1932 and 2019. The increase in the ice thinning rates and area loss during the last decade is coincident with the severe drought in this region (2010–present), which our minimal surface mass-balance model is able to reproduce. As a result of the glacier changes observed, the proglacial lake increased in area substantially between 1955 and 2019, with bathymetry data suggesting a water volume of 3.6 million m3 in 2017. This study highlights the need for further monitoring of glacierised areas in the Central Andes. Such efforts would facilitate a better understanding of the downstream impacts of glacier downwasting.

Ice shelves play a critical role in modulating dynamic loss of ice from the grounded portion of the Antarctic Ice Sheet and its contribution to sea-level rise. Measurements of ice-shelf motion provide insights into processes modifying buttressing. Here we investigate the effect of seasonal variability of basal melting on ice flow of Ross Ice Shelf. Velocities were measured from November 2015 to December 2016 at 12 GPS stations deployed from the ice front to 430 km upstream. The flow-parallel velocity anomaly at each station, relative to the annual mean, was small during early austral summer (November–January), negative during February–April, and positive during austral winter (May–September). The maximum velocity anomaly reached several metres per year at most stations. We used a 2-D ice-sheet model of the RIS and its grounded tributaries to explore the seasonal response of the ice sheet to time-varying basal melt rates. We find that melt-rate response to changes in summer upper-ocean heating near the ice front will affect the future flow of RIS and its tributary glaciers. However, modelled seasonal flow variations from increased summer basal melting near the ice front are much smaller than observed, suggesting that other as-yet-unidentified seasonal processes are currently dominant.

The cross-sectional profile of a glacial valley can be obtained with a variational principle in which the friction against the valley walls and the glacier bed is extremized, subject to a Lagrangian constraint. We show that the actual valley profile maximizes the friction, thus settling an old debate.

A volume-area scaling relation is commonly used to estimate glacier volume or its future changes on a global scale. The presence of an insulating supraglacial debris cover alters the mass-balance profile of a glacier, potentially modifying the scaling relation. Here, the nature of scaling relations for extensively debris-covered glaciers is investigated. Theoretical arguments suggest that the volume-area scaling exponent for these glaciers is ~7% smaller than that for clean glaciers. This is consistent with the results from flowline-model simulations of idealised glaciers, and the available data from the Himalaya. The best-fit scale factor for debris-covered Himalayan glaciers is ~60% larger compared to that for the clean ones, implying a significantly larger stored ice volume in a debris-covered glacier compared to a clean one having the same area. These results may help improve scaling-based estimates of glacier volume and future glacier changes in regions where debris-covered glaciers are abundant.