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.

The 2007–19 glaciological mass-balance series of Mera Glacier in the Everest Region, East Nepal, is reanalysed using the geodetic mass balance assessed by differencing two DEMs obtained from Pléiades stereo-images acquired in November 2012 and in October 2018. The glaciological glacier-wide annual mass balance of Mera Glacier has to be systematically decreased by 0.11 m w.e. a−1 to match the geodetic mass balance. We attribute part of the positive bias of the glaciological mass balance to an over-estimation of the accumulation above 5520 m a.s.l., likely due to a measurement network unable to capture its spatial variability. Over the period 2007–19, Mera Glacier has lost mass at a rate of −0.41 ± 0.20 m w.e. a−1, in general agreement with regional averages for the central Himalaya. We observe a succession of negative mass-balance years since 2013.

GRACE and ICESat Antarctic mass-balance differences are resolved utilizing their dependencies on corrections for changes in mass and volume of the same underlying mantle material forced by ice-loading changes. Modeled gravimetry corrections are 5.22 times altimetry corrections over East Antarctica (EA) and 4.51 times over West Antarctica (WA), with inferred mantle densities 4.75 and 4.11 g cm−3. Derived sensitivities (Sg, Sa) to bedrock motion enable calculation of motion (δB0) needed to equalize GRACE and ICESat mass changes during 2003–08. For EA, δB0 is −2.2 mm a−1 subsidence with mass matching at 150 Gt a−1, inland WA is −3.5 mm a−1 at 66 Gt a−1, and coastal WA is only −0.35 mm a−1 at −95 Gt a−1. WA subsidence is attributed to low mantle viscosity with faster responses to post-LGM deglaciation and to ice growth during Holocene grounding-line readvance. EA subsidence is attributed to Holocene dynamic thickening. With Antarctic Peninsula loss of −26 Gt a−1, the Antarctic total gain is 95 ± 25 Gt a−1 during 2003–08, compared to 144 ± 61 Gt a−1 from ERS1/2 during 1992–2001. Beginning in 2009, large increases in coastal WA dynamic losses overcame long-term EA and inland WA gains bringing Antarctica close to balance at −12 ± 64 Gt a−1 by 2012–16.

To assess the amount of ice volume stored in glaciers or ice caps, a method to estimate ice thickness distribution is required for glaciers where no direct observations are available. In this study, we use an existing inverse method to estimate the bedrock topography and ice thickness of the Renland Ice Cap, East Greenland, using satellite-based observations of the surface topography. The inverse approach involves a procedure in which an ice dynamical model is used to build-up an ice cap in steady state with climate forcing from a regional climate model, and the bedrock is iteratively adjusted until the modelled and observed surface topography match. We validate our model results against information from airborne radar data and satellite observed surface velocity, and we find that the inferred ice thickness and thereby the stored total volume of the ice cap is sensitive to the assumed ice softness and basal slipperiness. The best basal model parameters for the Renland Ice Cap are determined and the best estimated total ice volume of 384 km3 is found. The Renland Ice Cap is particularly interesting because of its location at a high elevation plateau and hence assumed low sensitivity to climate change.

Multi-decadal mass loss estimates are available for few glaciers of Central Asia. On Abramov Glacier (Pamir-Alay, Kyrgyzstan), comprehensive long-term glaciological measurements have been carried out from 1968 to 1999 and re-initiated in 2011. A climatological interpretation of this benchmark glacier in Central Asia requires bridging the gap between historical and renewed measurements. This is achieved here by computing the geodetic mass balance from 1975 to 2015 using previously unreleased Soviet aerial imagery and Pléaides stereo-imagery. During 1975–2015, Abramov Glacier lost 2.2 km2 (8.2%) of its area. The mean annual thickness change was − 0.43 ± 0.14 m a−1 for the period 1975–2015, corresponding to a volume change of − 0.45 ± 0.15 km3. The average specific geodetic mass balance amounts to − 0.38 ± 0.12 m w.e. a−1. The 1975–2015 glacier mass loss lies within the range of glaciological and geodetic mass-balance estimates that were previously published for disparate and shorter time intervals since 1968. This study covers a much longer time period than earlier geodetic estimates and demonstrates the capacity to geodetically constrain glacier change at high spatial resolution in Central Asia using historic aerial imagery and Structure from Motion techniques. Therefore, it could serve as a benchmark for future studies of regional mass change.

In the afternoon of 30 January 2017, a catastrophic outburst flood occurred in the Larsemann Hills (Princess Elizabeth Land, East Antarctica). The rapid drainage of both a thin supraglacial layer of water (near Boulder Lake) and Lake Ledyanoe into the englacial Lake Dålk provoked its overfill and outburst. As a result, a depression of 183 m × 220 m was formed in the place where Lake Dålk was located. This study summarises and clarifies the current state of knowledge on the flood that occurred in 2017. We present a phenomenological model of depression formation. We specify the reasons for the outburst of the system of lakes Boulder, Ledyanoe and Dålk. In addition, we carry out mathematical modelling of the outburst of each of the three lakes and estimate the flood severity. Outburst hydrographs, channel diameters, volume and duration of floods were calculated. Particular simulation results were validated with field data. In conclusion, we give an overview of the new outburst cycle of the lake system, which began in 2020 with the drainage of the lakes Boulder and Ledyanoe, and the new formation of Lake Dålk. Further research is required to improve our understanding of the lake system responses to changing external factors.

During the last few decades, bed-elevation profiles from radar sounders have been used to quantify bed roughness. Various methods have been employed, such as the ‘two-parameter’ technique that considers vertical and slope irregularities in topography, but they struggle to incorporate roughness at multiple spatial scales leading to a breakdown in their depiction of bed roughness where the relief is most complex. In this article, we describe a new algorithm, analogous to wavelet transformations, to quantify the bed roughness at multiple scales. The ‘Self-Adaptive Two-Parameter’ system calculates the roughness of a bed profile using a frequency-domain method, allowing the extraction of three characteristic factors: (1) slope, (2) skewness and (3) coefficient of variation. The multi-scale roughness is derived by weighted-summing of these frequency-related factors. We use idealized bed elevations to initially validate the algorithm, and then actual bed-elevation data are used to compare the new roughness index with other methods. We show the new technique is an effective tool for quantifying bed roughness from radar data, paving the way for improved continental-wide depictions of bed roughness and incorporation of this information into ice flow models.

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.

Due to climate change, worldwide glaciers are rapidly declining. The trend will continue into the future, with consequences for sea level, water availability and tourism. Here, we assess the future evolution of all glaciers in Scandinavia and Iceland until 2100 using the coupled surface mass-balance ice-flow model GloGEMflow. The model is initialised with three distinct past climate data products (E-OBS, ERA-I, ERA-5), while future climate is prescribed by both global and regional climate models (GCMs and RCMs), in order to analyze their impact on glacier evolution. By 2100, we project Scandinavian glaciers to lose between 67 ± 18% and 90 ± 7% of their present-day (2018) volume under a low (RCP2.6) and a high (RCP8.5) emission scenario, respectively. Over the same period, losses for Icelandic glaciers are projected to be between 43 ± 11% (RCP2.6) and 85 ± 7% (RCP8.5). The projected evolution is only little impacted by both the choice of climate data products used in the past and the spatial resolution of the future climate projections, with differences in the ice volume remaining by 2100 of 7 and 5%, respectively. This small sensitivity is attributed to our model calibration strategy that relies on observed glacier-specific mass balances and thus compensates for differences between climate forcing products.

Improving the knowledge on Himalayan glaciers mass balance is a key to understand the present and past annual atmospheric variations and future water availability in the region. Here, we present glaciological mass balance for Patsio Glacier, located in Himachal Pradesh (India), western Himalaya. Annual glacier-wide mass balance was measured for 7 consecutive years (2010/11 to 2016/17) and winter mass balance for 6 years (2011/12 to 2016/17). The cumulative mass balance over this period was −2.35 ± 0.37 m w.e. The corresponding mean mass balance was −0.34 m w.e. a−1. The mean annual ablation gradient excluding the debris-covered area was 0.47 m w.e. (100 m)−1. The annual ablation over the debris-covered area is reduced by an average of −1.0 m w.e. compared to the clean ice surface. Winter mass balance was consistently positive with a maximum of 1.34 m w.e. in 2014/15 and a minimum of 0.88 m w.e. in 2011/12. Multiple regression analysis between annual mass balance versus annual and winter precipitation of the Lahaul-Spiti region shows a significant positive correlation. Our results highlight the importance of monitoring seasonal mass balance and consideration of non-climatic parameters (debris and aspect) while estimating the glacier-wide mass balance.

The surface spectral albedo was measured over coastal landfast sea ice in Prydz Bay (off Zhongshan Station), East Antarctica from 5 October to 26 November of 2016. The mean albedo decreased from late-spring to early-summer, mainly responding to the change in surface conditions from dry (phase I) to wet (phase II). The evolution of the albedo was strongly influenced by the surface conditions, with alternation of frequent snowfall events and katabatic wind that induce snow blowing at the surface. The two phases and day-to-day albedo variability were more pronounced in the near-infrared albedo wavelengths than in the visible ones, as the near-infrared photons are more sensitive to snow metamorphism, and to changes in the uppermost millimeters and water content of the surface. The albedo diurnal cycle during clear sky conditions was asymmetric with respect to noon, decreasing from morning to evening over full and patchy snow cover, and decreasing more rapidly in the morning over bare ice. We conclude that snow and ice metamorphism and surface melting dominated over the solar elevation angle dependency in shaping the albedo evolution. However, we realize that more detailed surface observations are needed to clarify and quantify the role of the various surface processes.

Frontal ablation from tidewater glaciers is a major component of the total mass loss from the Greenland ice sheet. It remains unclear, however, how changes in atmospheric and oceanic temperatures translate into changes in frontal ablation, in part due to sparse observations at sufficiently high spatial and temporal resolution. We present high-frequency time-lapse imagery (photos every 30 min) of iceberg calving and meltwater plumes at Kangiata Nunaata Sermia (KNS), southwest Greenland, during June–October 2017, alongside satellite-derived ice velocities and modelled subglacial discharge. Early in the melt season, we infer a subglacial hydrological network with multiple outlets that would theoretically distribute discharge and enhance undercutting by submarine melt, an inference supported by our observations of terminus-wide calving during this period. During the melt season, we infer hydraulic evolution to a relatively more channelised subglacial drainage configuration, based on meltwater plume visibility indicating focused emergence of subglacial water; these observations coincide with a reduction in terminus-wide calving and transition to an incised planform terminus geometry. We suggest that temporal variations in subglacial discharge and near-terminus subglacial hydraulic efficiency exert considerable influence on calving and frontal ablation at KNS.

Bulk directional enhancement factors are determined for axisymmetric (girdle and single-maximum) orientation fabrics using a transversely isotropic grain rheology with an orientation-dependent non-linear grain fluidity. Compared to grain fluidities that are simplified as orientation independent, we find that bulk strain-rate enhancements for intermediate-to-strong axisymmetric fabrics can be up to a factor of ten larger, assuming stress homogenization over the polycrystal scale. Our work thus extends previous results based on simple basal slip (Schmid) grain rheologies to the transversely isotropic rheology, which has implications for large-scale anisotropic ice-flow modelling that relies on a transversely isotropic grain rheology. In order to derive bulk enhancement factors for arbitrary evolving fabrics, we expand the c-axis distribution in terms of a spherical harmonic series, which allows the rheology-required structure tensors through order eight to easily be calculated and provides an alternative to current structure-tensor-based modelling.

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.

Assessing the avalanche danger level requires snow stratigraphy and instability data. As such data are usually sparse, we investigated whether distributed snow cover modeling can be used to provide information on spatial instability patterns relevant for regional avalanche forecasting. Using Alpine3D, we performed spatially distributed simulations to evaluate snow instability for the winter season 2016–17 in the region of Davos, Switzerland. Meteorological data from automatic weather stations were interpolated to 100 m horizontal resolution and precipitation was scaled with snow depth measurements from airborne laser scanning. Modeled snow instability metrics assessed for two different weak layers suggested that the weak layer closer to the snow surface was more variable. Initially, it was less stable than the weak layer closer to the ground, yet it stabilized faster as the winter progressed. In spring, the simulated snowpack on south-facing slopes stabilized faster than on north-facing slopes, in line with the regional avalanche forecast. In the winter months January to March 2017, simulated instability metrics did not suggest that the snowpack on south-facing slopes was more stable, as reported in the regional avalanche forecast. Although a validation with field data is lacking, these model results still show the potential and challenges of distributed modeling for supporting operational avalanche forecasting.

With the assimilation of satellite-based sea-ice thickness (SIT) data, the new SIT reanalysis from the Towards an Operational Prediction system for the North Atlantic European coastal Zones (TOPAZ4) was released from 2014 to 2018. Apart from assimilating sea-ice concentration and oceanic variables, TOPAZ4 further assimilates CS2SMOS SIT. In this study, the 5-year reanalysis is compared with CS2SMOS, the Pan-Arctic Ice-Ocean Modeling and Assimilating System (PIOMAS) and the Combined Model and Satellite Thickness (CMST). Moreover, we evaluate TOPAZ4 SIT with field observations from upward-looking sonar (ULS), ice mass-balance buoys, Operation IceBridge Quicklook and Sea State Ship-borne Observations. The results indicate TOPAZ4 well reproduces the spatial characteristics of the Arctic SIT distributions, with large differences with CS2SMOS/PIOMAS/CMST mainly restricted to the Atlantic Sector and to the month of September. TOPAZ4 shows thinner ice in March and April, especially to the north of the Canadian Arctic Archipelago with a mean bias of −0.30 m when compared to IceBridge. Besides, TOPAZ4 simulates thicker ice in the Beaufort Sea when compared to ULS, with a mean bias of 0.11 m all year round. The benefit from assimilating SIT data in TOPAZ4 is reflected in a 34% improvement in root mean square deviation.

Spatio-temporal variability of the winter surface mass balance is a major uncertainty in the modelling of annual surface mass balance. Moreover, its measurement at high spatio-temporal resolution (sub-200 m) is very useful to force, calibrate or validate models. This study presents the results of year-round field campaigns to study the evolution of the surface mass balance in a ~2 km2 portion of the accumulation zone of the Mer de Glace (France). It is based on repeated LiDAR acquisitions, submergence-velocity measurements and meteorological records. The two methods used to quantify submergence velocities show good agreement. They present a linear temporal evolution without significant seasonal changes but display significant spatial variability. We conclude that a dense network of submergence velocity measurements is required to reduce the uncertainties when computing winter and annual surface mass balance from digital elevation model differencing. Finally, a hight spatio-temporal variability of the winter surface mass balance is highlighted (e.g., a std dev. of 0.92 m in April) even though the topography is homogeneous (std dev. of 25 m). Attempts to relate this variability to different morpho-topographic variables and wind-related indexes show the need for studies conducted at the snowfall event scale to obtain a better understanding of the variability in mass balance at the glacier scale.

Glacier surges are periodic episodes of mass redistribution characterized by dramatic increases in ice flow velocity and, sometimes, terminus advance. We use optical satellite imagery to document five previously unexamined surge events of Sít’ Kusá (Turner Glacier) in the St. Elias Mountains of Alaska from 1983 to 2013. Surge events had an average recurrence interval of ~5 years, making it the shortest known regular recurrence interval in the world. Surge events appear to initiate in the winter, with speeds reaching up to ~25 m d−1. The surges propagate down-glacier over ~2 years, resulting in maximum thinning of ~100 m in the reservoir zone and comparable thickening at the terminus. Collectively, the rapid recurrence interval, winter initiation and down-glacier propagation suggest Sít’ Kusá's surges are driven by periodic changes in subglacial hydrology and glacier sliding. Elevation change observations from the northern tributary show a kinematic disconnect above and below an icefall located 23 km from the terminus. We suggest the kinematic disconnect inhibits drawdown from the accumulation zone above the icefall, which leads to a steady flux of ice into the reservoir zone, and contributes to the glacier's exceptionally short recurrence interval.

Surface energy-balance models are commonly used in conjunction with satellite thermal imagery to estimate supraglacial debris thickness. Removing the need for local meteorological data in the debris thickness estimation workflow could improve the versatility and spatiotemporal application of debris thickness estimation. We evaluate the use of regional reanalysis data to derive debris thickness for two mountain glaciers using a surface energy-balance model. Results forced using ERA-5 agree with AWS-derived estimates to within 0.01 ± 0.05 m for Miage Glacier, Italy, and 0.01 ± 0.02 m for Khumbu Glacier, Nepal. ERA-5 data were then used to estimate spatiotemporal changes in debris thickness over a ~20-year period for Miage Glacier, Khumbu Glacier and Haut Glacier d'Arolla, Switzerland. We observe significant increases in debris thickness at the terminus for Haut Glacier d'Arolla and at the margins of the expanding debris cover at all glaciers. While simulated debris thickness was underestimated compared to point measurements in areas of thick debris, our approach can reconstruct glacier-scale debris thickness distribution and its temporal evolution over multiple decades. We find significant changes in debris thickness over areas of thin debris, areas susceptible to high ablation rates, where current knowledge of debris evolution is limited.