We use seismic refraction data to investigate the firn structure across a suture zone on the Amery Ice Shelf, East Antarctica, and the possible role of glacier dynamics in firn evolution. In the downstream direction, the data reveal decreasing compressional-wave velocities and increasing penetration depth of the propagating wave in the firn layer, consistent with $\sim$1 m firn thickening every 6 km. The boundary between the Lambert Glacier unit to the west and a major suture zone and the Mawson Escarpment Ice Stream unit to the east, is marked by differences in firn thicknesses, compressional-wave velocities and seismic anisotropy in the across-flow direction. The latter does not contradict the presence of a single-maximum crystal orientation fabric oriented 45–$90^{\circ }$ away from the flow direction. This is consistent with the presence of transverse simple shear governing the region's underlying ice flow regime, in association with elevated strain along the suture zone. The confirmation and quantification of the implied dynamic coupling between firn and the underlying ice requires integration of future seismic refraction, coring and modelling studies. Because firn is estimated to cover $\sim$98% of the Antarctic continent any such coupling may have widespread relevance to ice-sheet evolution and flow.

Measurements of elastic wave velocities enable non-destructive estimation of the mechanical properties, elastic moduli and density of snow and firn. The variation of elastic moduli with porosity in dry snow and firn is modeled using a differential effective medium scheme modified to account for the critical porosity above which the bulk and shear moduli of the ice frame vanish. A comparison of predicted and measured elastic moduli indicates that the shear modulus of ice in snow is lower than that computed from single crystal elastic stiffnesses of ice. This may indicate that the bonds between snow particles are more deformable under shear than under compression. A partial alignment of ice crystals also may contribute. Good agreement between elastic stiffnesses of the ice frame obtained from elastic wave velocity measurements and the predictions of the theory is observed. The approach is simple and compact, and does not require the use of empirical fits to the data. Owing to its simplicity, this model may prove useful in a variety of potential applications such as construction on snow, interpretation of seismic measurements to monitor and locate avalanches and estimation of density within compacting snow deposited on glaciers and ice sheets.

Up to 30% of the current tidewater mass loss in Svalbard corresponds to frontal ablation through submarine melting and calving. We developed two-dimensional (2-D) glacier–line–plume and glacier–fjord circulation coupled models, both including subglacial discharge, submarine melting and iceberg calving, to simulate Hansbreen–Hansbukta system, SW Svalbard. We ran both models for 20 weeks, throughout April–August 2010, using different scenarios of subglacial discharge and crevasse water depth. Both models showed large seasonal variations of submarine melting in response to transient fjord temperatures and subglacial discharges. Subglacial discharge intensity and crevasse water depth influenced calving rates. Using the best-fit configuration for both parameters our two coupled models predicted observed front positions reasonably well (±10 m). Although the two models showed different melt-undercutting front shapes, which affected the net-stress fields near the glacier front, no significant effects on the simulated glacier front positions were found. Cumulative calving (91 and 94 m) and submarine melting (108 and 118 m) along the simulated period showed in both models (glacier–plume and glacier–fjord) a 1:1.2 ratio of linear frontal ablation between the two mechanisms. Overall, both models performed well on predicting observed front positions when best-fit subglacial discharges were imposed, the glacier–plume model being 50 times computationally faster.

Distant glacial areas are interconnected by a complex system of fractures and water channels which run in the glacier interior and characterize the englacial realm. Water can slowly freeze in these channels where the slow freezing excludes air bubbles giving the ice a clear aspect. This ice is uplifted to the surface ablation zone by glacial movements and can therefore be observed in the form of clear surface ice bands. We employed an indirect method to sample englacial water by coring these ice bands. We were able, for the first time, to compare microbial communities sampled from clear (i.e. frozen englacial water bands) and cloudy ice (i.e. meteoric ice) through 16S rRNA gene sequencing. Although microbial communities were primarily shaped and structured by their spatial distribution on the glacier, ice type was a clear secondary factor. One area of the glacier, in particular, presented significant microbial community clear/cloudy ice differences. Although the clear ice and supraglacial communities showed typical cold-adapted glacial communities, the cloudy ice had a less defined glacial community and ubiquitous environmental organisms. These results highlight the role of englacial channels in the microbial dispersion within the glacier and, possibly, in the shaping of glacial microbial communities.

Lemon Creek Glacier, a temperate valley glacier in the Juneau Icefield of Southeast Alaska, is the site of long running (>60 years) glaciological studies. However, the most recent published estimates of its thickness and subglacial topography come from two ~50 years old sources that are not in agreement and do not account for the effects of years of negative mass balance. We collected a 1-km long active-source seismic line on the upper section of the glacier parallel and near to the centerline of the glacier, roughly straddling the equilibrium-line altitude. We used these data to perform joint reflection-refraction velocity modeling and reflection imaging of the glacier bed. We find that this upper section of Lemon Creek Glacier is as much as 150 m (~65%) thicker than previously suggested with a large overdeepening in an area previously believed to have a uniform thickness. Our results lead us to reinterpret the impact of basal motion on ice flow and have a significant impact on expectations of subglacial hydrology. We suggest that further efforts to develop a whole-glacier model of subglacial topography are necessary to support studies that require accurate models of ice thickness and subglacial topography.

Greenland's future contribution to sea-level rise is strongly dependent on the extent to which dynamic perturbations, originating at the margin, can drive increased ice flow within the ice-sheet interior. However, reported observations of ice dynamical change at distances >~50 km from the margin have a very low spatial and temporal resolution. Consequently, the likely response of the ice-sheet's interior to future oceanic and atmospheric warming is poorly constrained. Through combining GPS and satellite-image-derived ice velocity measurements, we measure multi-decadal (1993–1997 to 2014–2018) velocity change at 45 inland sites, encompassing all regions of the ice sheet. We observe an almost ubiquitous acceleration inland of tidewater glaciers in west Greenland, consistent with acceleration and retreat at glacier termini, suggesting that terminus perturbations have propagated considerable distances (>100 km) inland. In contrast, outside of Kangerlussuaq, we observe no acceleration inland of tidewater glaciers in east Greenland despite terminus retreat and near-terminus acceleration, and suggest propagation may be limited by the influence of basal topography and ice geometry. This pattern of inland dynamical change indicates that Greenland's future contribution to sea-level will be spatially complex and will depend on the capacity for dynamic changes at individual outlet glacier termini to propagate inland.

Casting snow is necessary to prevent metamorphism and deformation prior to X-ray micro-computed tomography (μCT) imaging. Current methods are insufficient for large-scale field sampling of snow due to safety considerations associated with the casting medium and/or lengthy sample preparation times. Here, a casting method using contrast-enhanced diethylphthalate (DEP) for μCT of snow is presented. The X-ray contrast of DEP is enhanced with barium titanate nanoparticles (BaTiO3) and iodine (I2). A partially unsupervised, three-phase segmentation method utilizing traditional Gaussian smoothing followed by a three-step process to address transition voxels is also presented. Synthetic images derived from real snow samples are used to evaluate the segmentation method with various configurations of trapped air bubbles. Real snow samples spanning a range of specific surface areas (SSAs) (8–28 m2 kg−1) and densities (135–463 kg m−3) are used to assess the performance of the segmentation method on real, cast samples. The method yields SSA, density and correlation length errors of less than 10% for synthetic images with air bubble surface areas less than 333 m−1 per sample volume for eight of the nine snow samples. For eight of the nine cast samples, the method yields errors of less than 10% for all three parameters.

In recent years, researchers have focused on the applications of uncrewed aerial vehicles (UAVs) in environmental remote sensing tasks. However, studies on glacier monitoring using UAV technology are relatively scarce, especially for high mountain glacier monitoring. To explore the feasibility of UAV technology for high mountain glaciers, four UAV surveys were deployed on two glaciers of the central Tibetan Plateau. Based on the images retrieved by UAV in 2017 and 2019, orthomosaics and digital elevation models were produced to quantify the length, area and elevation changes in the ablation zone of these two glaciers at different times. Additionally, we utilized several Landsat scenes to derive glacier changes over the last 30 years and combined these with the UAV data to assess the advantages and disadvantages of UAV technology in mountain glacier monitoring.

In this study, we explore small-scale (~1 to 20 km) thermal-refractive effects on basal geothermal heat flux (BGHF) at subglacial boundaries resulting from lateral thermal conductivity contrasts associated with subglacial topography and geologic contacts. We construct a series of two-dimensional, conductive, steady-state models that exclude many of the complexities of ice sheets in order to demonstrate the effect of thermal refraction. We show that heat can preferentially flow into or around a subglacial valley depending on the thermal conductivity contrast with underlying bedrock, with anomalies of local BGHF at the ice–bedrock interface between 80 and 120% of regional BGHF and temperature anomalies on the order of ±15% for the typical range of bedrock conductivities. In the absence of bed topography, subglacial contacts can produce significant heat flux and temperature anomalies that are locally extensive (>10 km). Thermal refraction can result in either an increase or decrease in the likelihood of melting and ice-sheet stability depending on the conductivity contrast and bed topography. While our models exclude many of the physical complexities of ice behavior, they illustrate the need to include refractive effects created by realistic geology into future glacial models to improve the prediction of subglacial melting and ice viscosity.

We develop a two-dimensional, plan-view formulation of ice-shelf flow and viscoelastic ice-shelf flexure. This formulation combines, for the first time, the shallow-shelf approximation for horizontal ice-shelf flow (and shallow-stream approximation for flow on lubricated beds such as where ice rises and rumples form), with the treatment of a thin-plate flexure. We demonstrate the treatment by performing two finite-element simulations: one of the relict pedestalled lake features that exist on some debris-covered ice shelves due to strong heterogeneity in surface ablation, and the other of ice rumpling in the grounding zone of an ice rise. The proposed treatment opens new venues to investigate physical processes that require coupling between the longitudinal deformation and vertical flexure, for instance, the effects of surface melting and supraglacial lakes on ice shelves, interactions with the sea swell, and many others.

Earlier gas measurements of firn air (atmosphere in open pore channels) at polar sites have revealed the occurrence of gas fractionation phenomena during bubble close-off, in addition to well-known thermal and gravitational gas separation. Nevertheless, because of difficulties posed by measurement, little is known about the distribution of air constituents in already closed pores (bubbles) in firn. Herein, we describe the application of high-sensitivity pure rotational Raman spectroscopy, combined with sample immersion in the fluorocarbon-based inert fluid for removing the optical disturbance by diffused reflection. That application efficiently elicits information about nitrogen and oxygen composition ratios (N2/O2 or O2/N2) for each air bubble in firn. The developed methodology presents important implications for elucidating how gas records are formed and modified in the course of pore close-off in polar firn.

Alaska's largest city, Anchorage, depends on Eklutna Glacier meltwater for drinking water and hydropower generation; however, the 29 km2 glacier is rapidly retreating. We used a temperature-index model forced with local weather station data to reconstruct the glacier's mass balance for the period 1985–2019 and quantify the impacts of glacier change on discharge. Model calibration involved a novel combination of in situ, geodetic mass-balance measurements and observed snowlines from satellite imagery. A resulting ensemble of 250 best-fitting model parameters was used to model mass balance and discharge. Eklutna Glacier experienced a significant negative trend (−0.31 m w.e. decade−1) in annual mean surface mass balance (mean: −0.62 ± 0.06 m w.e.). The day of the year when 95% of annual melt occurs was five days later in 2011–19 than in 1985–93, demonstrating a prolongation of melt season (May–September). Modeled mean specific discharge increased at 0.14 m decade−1, indicating peak water, the year when annual discharge reaches a maximum due to glacier retreat, has not been reached. Four of the five highest discharge years occurred since 2000. Increases in discharge quantity and melt season length require water resource managers consider future decreased discharge as the glacier continues to shrink.

The potential of capillary forces to retain water in pores is an important property of snow and firn at glaciers. Meltwater suspended in pores does not contribute to runoff and may refreeze during winter, which can affect the climatic mass balance and the subsurface density and temperature. However, measurement of firn water content is challenging and few values have been reported in the literature. Here, we use subsurface temperature and density measured at the accumulation zone of Lomonosovfonna (1200 m a.s.l.), Svalbard, to derive water content of the firn profiles after the 2014 and 2015 melt seasons. We do this by comparing measured and simulated rates of freezing front propagation. The calculated volumetric water content of firn is ~1.0–2.5 vol.% above the depth of 5 m and <0.5 vol.% below. Results derived using different thermistor strings suggest a prominent lateral variability in firn water content. Reported values are considerably lower than those commonly used in snow/firn models. This is interpreted as a result of preferential water flow in firn leaving dry volumes within wetted firn. This suggests that the implementation of irreducible water content values below 0.5 vol.% within snow/firn models should be considered at the initial phase of water infiltration.

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.

Processes governing meltwater penetration into cold firn remain poorly constrained. Here, in situ experiments are used to develop a grain-scale model to investigate physical limitations on meltwater infiltration in firn. At two sites in Greenland, drilling pumped water into cold firn to >75 m depth, and the thermo-hydrologic evolution of the firn column was measured. Rather than filling all available pore space, the water formed perched aquifers with downward penetration halted by thermal and density conditions. The two sites formed deep aquifers at ~40 m depth and at densities considerably less than the air pore close-off density (~725 kg m−3 at −18°C, and ~750 kg m−3 at −14°C), demonstrating that some pore space at depth remains inaccessible. A geometric grain-scale model of firn is constructed to quantify the limits of a descending fully saturated wetting front in cold firn. Agreement between the model and field data implies the model includes the first-order effects of water and heat flow in a firn lattice. The model constrains the relative importance of firn density, temperature and grain/pore size in inhibiting wetting front migration. Results imply that deep infiltration, including that which leads to firn aquifer formation, does not have access to all available firn pore space.

The dynamics of seasonal snow ablation on six glaciers in central Spitsbergen (Dicksonland) were assessed by examining a set of Sentinel-2 satellite images covering the summer ablation season for the period 2016–19. All glaciers lost 80% or more of their surface snow cover during the studied ablation seasons. This bolsters the recently observed trend of local glacier thinning, even at higher altitudes. Snow ablation dynamics are highly dependent on the glaciers altitudes, their position relative to the prevailing wind direction and the exposure to insolation. The accumulation areas of the studied glaciers were delimited based on the overlap of the minimum extent of snow-covered areas in the four consecutive studied summer seasons. The high temporal and spatial resolutions of available images enabled a detailed description of the seasonal snow ablation dynamics. Moreover, an estimate of the average number of days with below threshold glacier snow cover was made. This study contributes to our understanding of recent processes and might further support the modelling of glacier melt and subsequent runoff.

Few surges on the Central Tibetan Plateau have been reported. Here, we report observations of a recent surging event of the Gangjiaquba Glacier in the Geladandong Peak region using surface velocity and morphology changes that were extracted from Landsat MSS/TM/ETM+/OLI images obtained from 1973 to 2019. The results reveal that the active surge of this glacier initiated at the end of summer in 2014 and terminated in 2016. The surge resulted in a total advance of 500 ± 11.2 m and many fresh crevasses in the surging zone. The maximum velocity was 1100 m a−1 during the active surge phase, which is much smaller than those observed in Karakoram but similar to observations in West Kunlun.