We assess the evolution of glaciological structures during the 2003–05 surge in the Paulabreen glacier system, Svalbard. Glaciological structures on the glacier surface were mapped using aerial photographs captured in the early stages of the surge (2003) and 5 years after surge termination (2011). Three-dimensional measurements of glaciological structures were collected at the tidewater front in 2013. These datasets document the physical changes during (1) the late quiescent phase; (2) the early phase of the surge as the surge front propagated down Skobreen and advanced into Paulabreen and (3) the final stages of the surge following the surge front reaching the glacier terminus. Crevasse patterns and clusters of arcuate shear planes record zones of compressive and extensional flow associated with the downglacier progression of the surge front. The transfer of surging ice from Skobreen into Paulabreen caused lateral displacement of the medial moraines to the northeast. At the ice front, this movement tilted glaciological structures in the same direction. Structures at the southwest margin record strike–slip faulting and the elevation of debris into the ice in a zone of compression and transpression. We summarise these observations in a schematic reconstruction of structural evolution during the surge.

In situ observations of wave attenuation by sea ice are required to develop and validate wave–ice interaction parameterizations in coupled wave models. To estimate ice-induced wave attenuation in the field, the wave field is typically assumed to be stationary. In this study we investigate the validity of this assumption by creating a synthetic wave field in sea ice for different attenuation rates. We observe that errors in estimates of the wave attenuation rates are largest when attenuation rates are small or temporal averaging periods are short. Moreover, the adoption of the wave stationarity assumption can lead to negative estimates of the instantaneous wave attenuation rate. These apparent negative values should therefore not be attributed to wave growth or erroneous measurements a priori. Surprisingly, we observe that the validity of the wave stationarity assumption is irrelevant to the accuracy of estimates of wave attenuation rates as long as the temporal averaging period is taken sufficiently long. This may provide opportunities in using satellite-derived products to estimate wave attenuation rates in sea ice at global scales.

One of the largest contributors to uncertainty in predictions of sea-level rise from ice-sheet models is a lack of knowledge about the bed topography beneath ice sheets. Bed topography maps are normally made by interpolating between linear radar surveys using methods that include kriging, mass conservation and flowline diffusion, all of which may miss influential mesoscale (2–30 km) bedforms. Previous works have explored an Ice-Flow Perturbation Analysis (IFPA) approach for estimating bed topography using the surface expression of these mesoscale bedforms. Using regions of Pine Island Glacier that have been intensively surveyed by ice-penetrating radar as test sites, and a refined IFPA methodology, we find that IFPA detects bedforms capable of influencing ice flow which are not represented in Bedmachine Antarctica and other interpolated bed products. We further explore the ability of IFPA to estimate relative bed slipperiness, finding higher slipperiness in the main trunk and tributaries. Alongside other methods which estimate ice thickness, bed topography maps from IFPA have the potential to constrain projections of future sea-level rise, especially where radar data are sparse.

Thwaites Glacier (TG) plays an important role in future sea-level rise (SLR) contribution from the West Antarctic Ice Sheet. Recent observations show that TG is losing mass, and its grounding zone is retreating. Previous modeling has produced a wide range of results concerning whether, when, and how rapidly further retreat will occur under continued warming. These differences arise at least in part from ill-constrained processes, including friction from the bed, and future atmosphere and ocean forcing affecting ice-shelf and grounding-zone buttressing. Here, we apply the Ice Sheet and Sea-level System Model (ISSM) with a range of specifications of basal sliding behavior in response to varying ocean forcing. We find that basin-wide bed character strongly affects TG's response to sub-shelf melt by modulating how changes in driving stress are balanced by the bed as the glacier responds to external forcing. Resulting differences in dynamic thinning patterns alter modeled grounding-line retreat across Thwaites' catchment, affecting both modeled rates and magnitudes of SLR contribution from this critical sector of the ice sheet. Bed character introduces large uncertainties in projections of TG under equal external forcing, pointing to this as a crucial constraint needed in predictive models of West Antarctica.

Debris-covered glaciers are an important component of the mountain cryosphere and influence the hydrological contribution of glacierized basins to downstream rivers. This study examines the potential to make estimates of debris thickness, a critical variable to calculate the sub-debris melt, using ground-based thermal infrared radiometry (TIR) images. Over four days in August 2019, a ground-based, time-lapse TIR digital imaging radiometer recorded sequential thermal imagery of a debris-covered region of Peyto Glacier, Canadian Rockies, in conjunction with 44 manual excavations of debris thickness ranging from 10 to 110 cm, and concurrent meteorological observations. Inferring the correlation between measured debris thickness and TIR surface temperature as a base, the effectiveness of linear and exponential regression models for debris thickness estimation from surface temperature was explored. Optimal model performance (R2 of 0.7, RMSE of 10.3 cm) was obtained with a linear model applied to measurements taken on clear nights just before sunrise, but strong model performances were also obtained under complete cloud cover during daytime or nighttime with an exponential model. This work presents insights into the use of surface temperature and TIR observations to estimate debris thickness and gain knowledge of the state of debris-covered glacial ice and its potential hydrological contribution.

Ice streams are warmed by shear strain, both vertical shear near the bed and lateral shear at the margins. Warm ice deforms more easily, establishing a positive feedback loop in an ice stream where fast flow leads to warm ice and then to even faster flow. Here, we use radar attenuation measurements to show that the Siple Coast ice streams are colder than previously thought, which we hypothesize is due to along-flow advection of cold ice from upstream. We interpret the attenuation results within the context of previous ice-temperature measurements from nearby sites where hot-water boreholes were drilled. These in-situ temperatures are notably colder than model predictions, both in the ice streams and in an ice-stream shear margin. We then model ice temperature using a 1.5-dimensional numerical model which includes a parameterization for along-flow advection. Compared to analytical solutions, we find depth-averaged temperatures that are colder by 0.7°C in the Bindschadler Ice Stream, 2.7°C in the Kamb Ice Stream and 6.2–8.2°C in the Dragon Shear Margin of Whillans Ice Stream, closer to the borehole measurements at all locations. Modelled cooling corresponds to shear-margin thermal strengthening by 3–3.5 times compared to the warm-ice case, which must be compensated by some other weakening mechanism such as material damage or ice-crystal fabric anisotropy.

Ice records provide a qualitative rather than a quantitative indication of the trend of climate change. Using the bulk aerodynamic method and degree day model, this study quantified ice mass loss attributable to sublimation/evaporation (S/E) and meltwater on the basis of integrated observations (1960–2006) of glacier-related and atmospheric variables in the northeastern Tibetan Plateau. During 1961–2005, the average annual mass loss in the ice core was 95.33 ± 20.56 mm w.e. (minimum: 78.97 mm w.e. in 1967, maximum: 146.67 mm w.e. in 2001), while the average ratio of the revised annual ice accumulation was 21.2 ± 7.7% (minimum: 11.0% in 1992, maximum 44.8% in 2000). A quantitative formula expressing the relationship between S/E and air temperature at the monthly scale was established, which could be extended to estimation of S/E changes of other glaciers in other regions. The elevation effect on alpine precipitation determined using revised ice accumulation and instrumental data was found remarkable. This work established a method for quantitative assessment of the temporal variation in ice core mass loss, and advanced the reconstruction of long-term precipitation at high elevations. Importantly, the formula established for reconstruction of S/E from temperature time series data could be used in other regions.

Glacial hydraulic tremor (GHT) can be monitored to observe changes in location and distribution of water flow beneath glacial ice, allowing the spatiotemporal evolution of subglacial hydrology to be studied continuously and remotely. We use frequency-dependent polarization analysis (FDPA) to classify types of GHT and assess its spatio-temporal extent beneath Rhonegletscher, Switzerland, in a continuous seismic record through 2018 and 2019 at three ice-proximal bedrock-based seismometers. We determine the frequency bands composing the GHT and calculate back azimuth angles pointing to a previously known subglacial channel. We also inspect the relationship between GHT seismic power and water discharge from the glacier to compare daily and seasonal shifts with the observed GHT. We observed the seasonal shift from a distributed system to a channelized system, and our dataset allowed comparison of channel locations within and across seasons, with implications for sediment evacuation and bed erosion. The successful use of this method to assess GHT previously on Taku glacier (the methods of which this project follows) and now Rhonegletscher shows that existing ice-proximal passive seismic installations can be used to easily and continuously monitor subglacial hydraulic activity.

To examine the influence of global warming including increased heavy snowfall frequency on the potential of natural dry snow avalanche frequency and magnitude, we estimated the frequency of weak layer formation and the associated slab overload above the weak layer over northern Japan. The estimation was numerically performed using climate models' output for 1800 winter simulations in each of the historical (1951–2010) and +4°C experiments by forcing a physical-based snowpack model with the result of the climate models. Here the +4°C experiment was defined as a climate when the global mean air temperature had increased by 4°C from the preindustrial level. The estimation results showed that the probability of weak layer formation, identified by the natural stability index, would decrease all over the area because of the shorter age of the weak layers caused by a warmer climate, indirectly indicating a potential decrease in avalanche frequency. However, because of increased heavy snowfall frequency, slab overload would increase by 10–15% in inland areas for weak layers of decomposing fragments/precipitation particles and the mountainous area facing the Sea of Japan for weak layers of facets/depth hoar, thereby potentially indicating an increased magnitude of avalanches.

Debris-covered ice is widespread in mountain regions with debris an important control on surface ice melt and glacier retreat. Quantifying debris cover extent and its evolution through time over large regions remains a challenge. This study develops two Normalized Difference Supraglacial Debris Indices for mapping debris-covered ice based on thermal and near Infrared Landsat 8 bands. They were calibrated with field data. Validation suggests that they have a high level of accuracy. They are then applied to Landsat data for 2016 to produce the first detailed glacier inventory of the Afghanistan Hindu Kush Himalaya that includes debris cover. 3408 glaciers were identified which, for those ⩾0.01 km2 in area, gives an ice cover of 2,222 ± 11 km2 and a debris cover of 619 ± 40 km2. Principal components analysis was used to identify the most influential drivers of debris-covered ice extent. Lower proportions of debris cover were associated with glaciers with a higher elevation range, that were larger, longer and wider. These relations were statistically clearer when the dataset was broken down into climate and geological zones. A glaciers continue to shrink, the proportion of debris cover will become higher, making it more important to map debris cover reliably.

Characterising the structures within glaciers can give unique insight into ice motion processes. On debris-covered glaciers, traditional structural glaciological mapping is challenging because the lower glacier is hidden by the supraglacial debris layer. Here, we use high-resolution optical televiewer (OPTV) image logs from four boreholes drilled into Khumbu Glacier, Nepal, to overcome this limitation and investigate englacial structural features within a Himalayan debris-covered glacier. The OPTV logs show structural features that are up to an order of magnitude thinner than those observed at the glacier surface and reveal five structural units: (I) primary stratification of ice; (II) debris-rich planes that conform with the primary stratification; (III) water-healed crevasse traces; (IV) healed crevasse traces; and (V) steeply dipping planes of basally derived fine sediment near the glacier terminus. The OPTV logs also reveal that the primary stratification both decreases in dip with depth (by up to 56° over 20 m) and rotates with depth (by up to 100° over 20 m) towards parallelism with the proximal lateral moraine. This transformation and the presence of relict layers of basally derived sediment raised into an englacial position – possibly involving thrusting – near the glacier's now stagnant terminus reveal a previously more dynamic glacier regime.

Global glacier mass balance decreased rapidly over the last two decades, exceeding mass loss from the Greenland and Antarctic Ice Sheets. In Greenland, peripheral glaciers and ice caps (GICs) cover only ~5% of Greenland's area but contributed ~20% of the island's ice mass loss between 2000 and 2018. Although Greenland GIC mass loss due to surface meltwater runoff has been estimated using atmospheric models, mass lost to changes in ice discharge into oceans (i.e., dynamic mass loss) remains unquantified. We use the flux gate method to estimate discharge from Greenland's 585 marine-terminating peripheral glaciers between 1985 and 2018, and compute dynamic mass loss as the discharge anomaly relative to the 1985–98 period. Greenland GICs discharged between 2.94 ± 0.23 and 4.03 ± 0.23 Gt a−1 from 1985 to 1998, depending on the gap-filling method, and abruptly increased to 5.10 ± 0.21 Gt a−1 from 1999 to 2018. The resultant ~1–2 Gt a−1 dynamic mass loss was driven by synchronous widespread acceleration around Greenland. The mass loss came predominantly from the southeast region, which contains 39% of the glaciers. Although changes in discharge over time were small relative to surface mass-balance changes, our speed and discharge time series suggest these glaciers may quickly accelerate in response to changes in climate.

Increased surface melt in the percolation zone of the Greenland ice sheet causes significant changes in the firn structure, directly affecting the amount and timing of meltwater runoff. Here we force an energy-balance model with automatic weather stations data at two sites in the percolation zone of southwest Greenland ($2040$ and 2360 m a.s.l.) between spring $2017$ and fall $2019$. Extensive model validation and sensitivity analysis reveal that the skin layer formulation used to compute the surface temperature by closing the energy balance leads to a consistent overestimation of melt by more than a factor of two or three depending on the site. In contrast, model results match the observations well when the model is forced by observed surface temperatures; however, unexplained residuals in the energy balance occur. The sensible and ground heat flux differ markedly in the two simulations accounting largely for the difference in modeled melt amounts. This indicates that the energy available for melt is highly sensitive to small changes in surface temperature. Thus, regional climate models that also use the skin layer formulation may have a bias in surface temperature and melt energy in the percolation zone of the ice sheet.

Ice aprons are defined as very small ice bodies covering steep rock slopes. They have only been the subject of increased scientific interest for a few years, despite the fact that they are a condition for mountaineering and obvious elements in the high-alpine landscapes. However, very little is known about their distribution, evolution and physical characteristics. In this paper, we review the existing knowledge on ice aprons, which have almost exclusively been investigated in the Mont-Blanc massif, Western Alps. We supplement this review with novel results from recent surveys of ice aprons. We used a wide array of methodologies, from remote sensing (multi-source imagery) to in situ (stakes and thermal monitoring) and laboratory (radiocarbon dating and texture analysis) glaciological investigations. In the Mont-Blanc massif, ice aprons occupy 4.2 km2 within the alpine permafrost zone. Temperature measured at the ice–rock interface is indeed largely negative. Thinness of ice aprons coupled with the cold context implies a quasi-stationary shear regime without basal Sliding. Only ice at the surface can possibly melt in warm periods. After a shrinking period from the end of the Little Ice Age to the mid-to-late-1960s, ice aprons experienced a short period of expansion, followed by an accelerated shrinkage since the beginning of the 21st century. This shrinkage now favours rockfall triggering and poses a serious threat to a glaciological heritage since ice aprons host several-thousand-year-old ice. Finally, we synthesize this information to assess the existing definition of ice aprons, and propose some future research directions.

We perform Bayesian inference of the parameters of a time-dependent model of ice flow and calving at Helheim Glacier, East Greenland. We find that, while a time-independent calving parameterization can recover the mean observed terminus position, such a model is unable to recover sub-annual variability, even when forced with seasonally varying climate. To address this, we develop a simple stochastic model relating surface runoff rates and calving threshold. Again inferring model parameters from observations, we find that this parameterization is able to reproduce observations with respect to both mean position and characteristic temporal variability. This result demonstrates the importance of considering potential sub-annual controls on calving rates in numerical models, which may include variable undercutting rates or surface runoff-dependent surface crevasse propagation.

We developed a novel laser melting sampler (LMS) for ice cores to measure the stable water isotope ratios (δ18O and δD) as temperature proxies at sub-centimeter depth resolutions. In this LMS system, a 2 mm diameter movable evacuation nozzle holds an optical fiber through which a laser beam irradiates the ice core. The movable nozzle intrudes into the ice core, the laser radiation meanwhile melts the ice cylindrically, and the meltwater is pumped away simultaneously through the same nozzle and transferred to a vial for analysis. To avoid isotopic fractionation of the ice through vaporization, the laser power is adjusted to ensure that the temperature of the meltwater is always kept well below its boiling point. A segment of a Dome Fuji shallow ice core (Antarctica), using the LMS, was then demonstrated to have been discretely sampled with a depth resolution as small as 3 mm: subsequent analysis of δ18O, δD, and deuterium excess (d) was consistent with results obtained by hand segmentation within measurement uncertainties. With system software to control sampling resolution, the LMS will enable us to identify temperature variations that may be detectable only at sub-centimeter resolutions in ice cores.

Himalayan glaciers have been shrinking and losing mass rapidly since 1970s with an enhanced rate after 2000. The shrinkage is, however, quite heterogeneous and it is important to document individual glacier characteristics and their changes at the basin scale. We present an updated glacier inventory of the Upper Alaknanda Basin (UAB), Central Himalaya for the year 2020 and report area, debris cover and length changes for the periods 1994–2006 and 2006–2020 based on remote-sensing data. We identified 198 glaciers, comprising an area of 354.6 ± 8.5 km2, and classified them according to their size and morphology. The glaciers of the basin lost 4.2 ± 2.9% (0.16 ± 0.11% a−1) of their frontal area (from 368.6 ± 9.2 to 353.0 ± 5.3 km2) from 1994 to 2020. The average retreat rate was higher in the period 2006–2020 (13.3 ± 1.8 m a−1) in comparison to 1994–2006 (9.3 ± 1.9 m a−1). However, the area change rate was similar for the two periods (0.14 ± 0.27% a−1 for 1994–2006 and 0.16 ± 0.19% a−1 for 2006–2020). The debris-covered area has increased by 13.4 ± 4.4% from 1994 to 2020. A comparison with previous studies in UAB indicates consistent area loss of ~0.15% a−1 since the 1960s.

Owing to drifting snow processes, snow accumulation and surface density in polar environments are variable in space and time. We present new field data of manual measurements, repeat terrestrial laser scanning and snow micro-penetrometry from Dronning Maud Land, Antarctica, showing the density of new snow accumulations. We combine these data with published drifting snow mass flux observations, to evaluate the performance of the 1-D, detailed, physics-based snow cover model SNOWPACK in representing drifting snow and surface density. For two sites in East Antarctica with multiple years of data, we found a coefficient of determination for the simulated drifting snow of r2 = 0.42 and r2 = 0.50, respectively. The field observations show the existence of low-density snow accumulations during low wind conditions. Successive high wind speed events generally erode these low-density layers while producing spatially variable erosion/deposition patterns with typical length scales of a few metres. We found that a model setup that is able to represent low-density snow accumulating during low wind speed conditions, as well as subsequent snow erosion and redeposition at higher densities during drifting snow events was mostly able to describe the observed temporal variability of surface density in the field.

We present the surface mass balance (SMB) dataset from Vostok Station's accumulation stake farms which provide the longest instrumental record of its kind obtained with a uniform technique in central Antarctica over the last 53 years. The snow build-up values at individual stakes demonstrate a strong random scatter related to the interaction of wind-driven snow with snow micro-relief. Because of this depositional noise, the signal-to-noise ratio (SNR) in individual SMB time series derived at single points (from stakes, snow pits or firn cores) is as low as 0.045. Averaging the data over the whole stake farm increases the SNR to 2.3 and thus allows us to investigate reliably the climatic variability of the SMB. Since 1970, the average snow accumulation rate at Vostok has been 22.5 ± 1.3 kg m−2 yr−1. Our data suggest an overall increase of the SMB during the observation period accompanied by a significant decadal variability. The main driver of this variability is local air temperature with an SMB temperature sensitivity of 2.4 ± 0.2 kg m−2 yr−1 K−1 (11 ± 2% K−1). A covariation between the Vostok SMB and the Southern Oscillation Index is also observed.

A conceptual model, based on field observations and assumed physics of a perennial firn aquifer near Helheim Glacier (southeast Greenland), is evaluated via steady-state 2-D simulation of liquid water flow and energy transport with phase change. The simulation approach allows natural representation of flow and energy advection and conduction that occur in vertical meltwater recharge through the unsaturated zone and in lateral flow within the saturated aquifer. Agreement between measured and simulated aquifer geometry, temperature, and recharge and discharge rates confirms that the conceptual field-data-based description of the aquifer is consistent with the primary physical processes of groundwater flow, energy transport and phase change. Factors that are found to control simulated aquifer configuration include surface temperature, meltwater recharge rate, residual total-water saturation and capillary fringe thickness. Simulation analyses indicate that the size of perennial firn aquifers depends primarily on recharge rates from surface snowmelt. Results also imply that the recent aquifer expansion, likely due to a warming climate, may eventually produce lakes on the ice-sheet surface that would affect the surface energy balance.