We report detailed chemical and isotopic data from a subglacial siliceous deposit on andesitic bedrock recently exposed by glacier retreat. Whereas a single, <1 μm, Si-rich layer covers the highly polished bedrock on the up-glacier (stoss) surfaces, distinct, lithified deposits commonly occur at the lee of small bedrock protuberances, on a scale <0.1 meter. The deposit is millimeters in thickness and consists of laminae tens to hundreds microns thick that differ from one another in color, rock-fragment abundance and chemical composition. Ca-rich laminae that are sufficiently enriched in uranium (~2–50 ppm) to permit U-series isotopic analysis suggest that the subglacial deposit formed 10–20 ka, much earlier than previously assumed. We conclude that (1) the siliceous deposit persisted for at least 10 000 years despite the intervening erosion and weathering, (2) distinct episodes of formation due to significant changes in hydrology and water chemistry are recorded in the deposit, and (3) a siliceous slurry may have existed at the ice-rock interface and influenced the local friction. This work reinforces earlier findings that subglacial chemical deposits can form and persist on geologic time scales and may have implications for the role of the cryosphere in the Earth's geochemical cycles and climate system.

Surface melting on Amery Ice Shelf (AIS), East Antarctica, produces an extensive supraglacial drainage system consisting of hundreds of lakes connected by surface channels. This drainage system forms most summers on the southern portion of AIS, transporting meltwater large distances northward, toward the ice front and terminating in lakes. Here we use satellite imagery, Landsat (1, 4 and 8), MODIS multispectral and Sentinel-1 synthetic aperture radar to examine the seasonal and interannual evolution of the drainage system over nearly five decades (1972–2019). We estimate seasonal meltwater input to one lake by integrating output from the regional climate model [Regional Atmospheric Climate Model (RACMO 2.3p2)] over its catchment defined using the Reference Elevation Model of Antarctica. We find only weak positive relationships between modeled seasonal meltwater input and lake area and between meltwater input and lake volume. Consecutive years of extensive melting lead to year-on-year expansion of the drainage system, potentially through a link between melt production, refreezing in firn and the maximum extent of the lakes at the downstream termini of drainage. These mechanisms are important when evaluating the potential of drainage systems to grow in response to increased melting, delivering meltwater to areas of ice shelves vulnerable to hydrofracture.

Englacial water transport is an integral part of the glacial hydrologic system, yet the geometry of englacial structures remains largely unknown. In this study, we explore the excitation of fluid resonance by small amplitude waves as a probe of englacial geometry. We model a hydraulic network consisting of one or more tabular cracks that intersect a cylindrical conduit, subject to oscillatory wave motion initiated at the water surface. Resulting resonant frequencies and quality factors are diagnostic of fluid properties and geometry of the englacial system. For a single crack–conduit system, the fundamental mode involves gravity-driven fluid sloshing between the conduit and the crack, at frequencies between 0.02 and 10 Hz for typical glacial parameters. Higher frequency modes include dispersive Krauklis waves generated within the crack and tube waves in the conduit. But we find that crack lengths are often well constrained by fundamental mode frequency and damping rate alone for settings that include alpine glaciers and ice sheets. Branching crack geometry and dip, ice thickness and source excitation function help define limits of crack detectability for this mode. In general, we suggest that identification of eigenmodes associated with wave motion in time series data may provide a pathway toward inferring englacial hydrologic structures.

The inaugural data from the first systematic program of sea-ice observations in Kotzebue Sound, Alaska, in 2018 coincided with the first winter in living memory when the Sound was not choked with ice. The following winter of 2018–19 was even warmer and characterized by even less ice. Here we discuss the mass balance of landfast ice near Kotzebue (Qikiqtaġruk) during these two anomalously warm winters. We use in situ observations and a 1-D thermodynamic model to address three research questions developed in partnership with an Indigenous Advisory Council. In doing so, we improve our understanding of connections between landfast ice mass balance, marine mammals and subsistence hunting. Specifically, we show: (i) ice growth stopped unusually early due to strong vertical ocean heat flux, which also likely contributed to early start to bearded seal hunting; (ii) unusually thin ice contributed to widespread surface flooding. The associated snow ice formation partly offset the reduced ice growth, but the flooding likely had a negative impact on ringed seal habitat; (iii) sea ice near Kotzebue during the winters of 2017–18 and 2018–19 was likely the thinnest since at least 1945, driven by a combination of warm air temperatures and a persistent ocean heat flux.

We present polar ice core organic matter (OM) fluorescence signatures to reconstruct ancient and modern atmospheric compositions and relate OM signals to past ecological changes. OM composition from three Arctic ice cores (Canada and Greenland) was characterized by fluorescence spectroscopy and compared to an Antarctic OM record. Diverse OM was measured in ancient and modern ice in both hemispheres and similarities existed across vast spatiotemporal scales. We determined three OM markers, indicating paleoclimate and modern carbon trends: (i) ‘humic-like’, detected in Holocene ice of more complex and aromatic character, supporting trends of higher plant influences in warmer climates, (ii) monolignol- and non-amino acid-like, describing simple, lignin-like OM precursors ubiquitous in the environment and the microbial degradation products of more complex materials from plants/soils, and (iii) amino acid- and tannin-like, indicating microbial degradation of simple OM chemical species, compared to the other markers. Concentration trends were inferred from fluorescence intensities of individual OM types and related to warmer temperatures. No indicators of freshly produced OM by microbes were detected; signals were interpreted as materials externally produced from the ice and transported to polar regions. This marks the first global comparison of atmospheric reconstructions from OM across vast spatiotemporal scales.

A variation in the α−β model which is a regression model that allows a deterministic prediction of the extreme runout to be expected in a given path, was applied for calculating avalanche runout in the Catalan Pyrenees. Present knowledge of major avalanche activity in this region and current mapping tools were used. The model was derived using a dataset of 97 ‘extreme’ avalanches that occurred from the end of 19th century to the beginning of 21st century. A multiple linear regression model was obtained using three independent variables: inclination of the avalanche path, horizontal length and area of the starting zone, with a good fit of the function (R2 = 0.81). A larger starting zone increases the runout and a larger length of the path reduces the runout. The new updated equation predicts avalanche runout for a return period of ~100 years. To study which terrain variables explain the extreme values of the avalanche dataset, a comparative analysis of variables that influence a longer or shorter runout was performed. The most extreme avalanches were treated. The size of the avalanche path and the aspect of the starting zone showed certain association between avalanches with longer or shorter runouts.

Southern Ocean sea ice can exhibit widespread flooding and subsequent snow-ice formation, due to relatively thick snow covers compared to the total ice thickness. Considerable subkilometer scale variability in snow and ice thickness causes poorly constrained uncertainties in determining the amount of flooding that occurs. Using datasets of snow depth and ice thickness acquired in the Weddell Sea during austral winter 2013 (AWECS campaign) from three floes, we demonstrate large spatial variability of a factor 10 and 5 for snow and combined snow and ice thickness, respectively. The temporal evolution after the floe visit was recorded by automatic weather station and ice mass balance buoys. Using a physics-based, multi-layer snow/sea ice model in a one-dimensional and distributed mode to simulate the thermodynamic processes, we show that the distributed simulations, modeling flooding across the entire heterogeneous floe, produced vastly different amounts of flooding than one-dimensional single point simulations. Three times the flooding is produced in the one-dimensional simulation for the buoy location than distributed (floe-averaged) simulations. The latter is in close agreement with buoy observations. The results suggest that using point observations or one-dimensional simulations to extrapolate processes on the floe-scale can overestimate the amount of flooding and snow-ice formation.

Accurate knowledge of the ice thickness distribution and glacier bed topography is essential for predicting dynamic glacier changes and the future developments of downstream hydrology, which are impacting the energy sector, tourism industry and natural hazard management. Using AIR-ETH, a new helicopter-borne ground-penetrating radar (GPR) platform, we measured the ice thickness of all large and most medium-sized glaciers in the Swiss Alps during the years 2016–20. Most of these had either never or only partially been surveyed before. With this new dataset, 251 glaciers – making up 81% of the glacierized area – are now covered by GPR surveys. For obtaining a comprehensive estimate of the overall glacier ice volume, ice thickness distribution and glacier bed topography, we combined this large amount of data with two independent modeling algorithms. This resulted in new maps of the glacier bed topography with unprecedented accuracy. The total glacier volume in the Swiss Alps was determined to be 58.7 ± 2.5 km3 in the year 2016. By projecting these results based on mass-balance data, we estimated a total ice volume of 52.9 ± 2.7 km3 for the year 2020. Data and modeling results are accessible in the form of the SwissGlacierThickness-R2020 data package.

To get a better overview of atmosphere-driven mass changes at the 79N Glacier (Nioghalvfjerdsfjorden Glacier), the largest outlet glacier of the northeast Greenland ice stream, the surface mass balance (SMB) is modeled by linking the COupled Snowpack and Ice surface energy and mass-balance model in PYthon (COSIPY) with the output of a regional atmospheric model (Polar WRF) for the years 2014–2018. After a manual model optimization, the model produces reliable results when compared to observations in the region and to values from the literature. High spatial resolution (1 km) simulations reveal strong interannual variability of the SMB. Stronger surface melting increased the ablation and runoff in years with high mass loss (2016 and 2017) whereas in other years (2015 and 2018) melting and refreezing inside the snowpack dominated the mass balance (MB). A cooler regional climate with higher snowfall-driven accumulation, higher albedo and reduced surface melt in the ablation period of 2018 resulted in a positive SMB in 2018, however, the annual total MB remained negative. The results suggest a promising new dataset for gaining more insights into mass-balance processes and their contribution to the acceleration of glacier retreat in northeast Greenland.

The thickness of a supraglacial layer is critical to the magnitude and time frame of glacier melt. Field-based, short pulse, ground-penetrating radar (GPR) has successfully measured debris thickness during a glacier's melt season, when there is a strong return from the ice–debris interface, but profiling with GPR in the absence of a highly reflective ice interface has not been explored. We investigated the performance of 960 MHz signals over 2 km of transects on Changri Nup Glacier, Nepal, during the post-monsoon. We also performed laboratory experiments to interpret the field data and investigate electromagnetic wave propagation into dry rocky debris. Laboratory tests confirmed wave penetration into the glacier ice and suggest that the ice–debris interface return was missing in field data because of a weak dielectric contrast between solid ice and porous dry debris. We developed a new method to estimate debris thicknesses by applying a statistical approach to volumetric backscatter, and our backscatter-based calculated thickness retrievals gave reasonable agreement with debris depths measured manually in the field (10–40 cm). We conclude that, when melt season profiling is not an option, a remote system near 1 GHz could allow dry debris thickness to be estimated based on volumetric backscatter.

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.

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.

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.

Microstructures, including crystallographic fabric, within the margin of streaming ice can exert strong control on flow dynamics. To characterize a natural setting, we retrieved three cores, two of which reached bed, from the flank of Jarvis Glacier, eastern Alaska Range, Alaska. The core sites lie ~1 km downstream of the source, with abundant water present in the extracted cores and at the base of the glacier. All cores exhibit dipping layers, a combination of debris bands and bubble-free domains. Grain sizes coarsen on average approaching the lateral margin. Crystallographic orientations are more clustered and with c-axes closer to horizontal nearer the lateral margin. The measured fabric is sufficiently weak to induce little mechanical anisotropy, but the data suggest that despite the challenging conditions of warm ice, abundant water and a short flow distance, many aspects of the microstructure, including measurable crystallographic fabric, evolved in systematic ways.

We calibrate the calving parameterisation implemented in the Open Global Glacier Model via two methods (velocity constraint and surface mass balance (SMB) constraint) and assess the impact of accounting for frontal ablation on the ice volume estimate of Greenland tidewater peripheral glaciers (PGs). We estimate an average regional frontal ablation flux of 7.38±3.45 Gta−1 after calibrating the model with two different satellite velocity products, and of 0.69±0.49 Gta−1 if the model is constrained using frontal ablation fluxes derived from independent modelled SMB averaged over an equilibrium reference period (1961–90). This second method makes the assumption that most PGs during that time have an equilibrium between mass gain via SMB and mass loss via frontal ablation. This assumption serves as a basis to assess the order of magnitude of dynamic mass loss of glaciers when compared to the SMB imbalance. The differences between results from both methods indicate how strong the dynamic imbalance might have been for PGs during that reference period. Including frontal ablation increases the estimated regional ice volume of PGs, from 14.47 to 14.64±0.12 mm sea level equivalent when using the SMB method and to 15.84±0.32 mm sea level equivalent when using the velocity method.

Hailstone structures have been studied for over a century, but so far mainly by manual optical means. This paper presents new texture and microstructure data (i.e. crystal lattice orientations, grain sizes and shapes) measured with an Automatic Ice Texture Analyzer, which gives access to high spatial and angular resolutions. The hailstones show two main characteristics: (1) they are structured with several concentric layers composed of alternating fine equiaxed grains and coarse elongated and radially oriented grains, and (2) they show two texture types with c-axes oriented either parallel or perpendicular to the radial direction. Such textures are compared with the ones observed in lake S1 and S2 ices, respectively. The S1 texture type (with c-axes parallel to the columnar crystals that grew in the radial direction) may result from epitaxial growth from a polycrystalline embryo, while the S2 texture (c-axes in the plane perpendicular to the column direction) may result from the growth from an embryo made of a few crystals with mainly one crystallographic orientation. Our novel high-resolution maps and measurements of both microstructure and texture may help to shed new light on the long-term discussion on the growth mechanisms of large hailstones.

The current paper studies the dynamics and age of the Triangle du Tacul (TDT) ice apron, a massive ice volume lying on a steep high-mountain rock wall in the French side of the Mont-Blanc massif at an altitude close to 3640 m a.s.l. Three 60 cm long ice cores were drilled to bedrock (i.e. the rock wall) in 2018 and 2019 at the TDT ice apron. Texture (microstructure and lattice-preferred orientation, LPO) analyses were performed on one core. The two remaining cores were used for radiocarbon dating of the particulate organic carbon fraction (three samples in total). Microstructure and LPO do not substantially vary with along the axis of the ice core. Throughout the core, irregularly shaped grains, associated with strain-induced grain boundary migration and strong single maximum LPO, were observed. Measurements indicate that at the TDT ice deforms under a low strain-rate simple shear regime, with a shear plane parallel to the surface slope of the ice apron. Dynamic recrystallization stands out as the major mechanism for grain growth. Micro-radiocarbon dating indicates that the TDT ice becomes older with depth perpendicular to the ice surface. We observed ice ages older than 600 year BP and at the base of the lowest 30 cm older than 3000 years.

The density–depth relationship of the Antarctic ice sheet is important for establishing a high-precision surface mass balance model and predicting future ice-sheet contributions to global sea levels. A new algorithm is used to reconstruct firn density and densification rate by inverting monostatic radio wave echoes from ground-operated frequency-modulated continuous wave radar data collected near four ice cores along the transect from Zhongshan Station to Dome A. The inverted density profile is consistent with the core data within 5.54% root mean square error. Due to snow redistribution, the densification rate within 88 km of ice core DT401 is correlated with the accumulation rate and varies greatly over horizontal distances of <5 km. That is, the depth at which a critical density of 830 kg m−3 is reached decreases and densification rate increases in high-accumulation regions but decreases in low-accumulation regions. This inversion technique can be used to analyse more Antarctic radar data and obtain the density distribution trend, which can improve the accuracy of mass-balance estimations.

Atmospheric ice accretion results from the exposure of technical equipment or facilities to cold and humid environments. Supercooled droplets in a cloud can impact an airplane's surface and quickly form an ice layer. The presence of air pockets in such a layer is well known and explains the white appearance of some of the accretions. However, estimation of its porosity values and studies on the pore formation mechanics remain limited. In this study, we performed tests in an icing wind tunnel and scans with micro-computed tomography to address these issues. Here, we show that the accretion has closed porosity below 1%, which is mostly produced by the interaction between a spray-like impact on the water surface. The insights we provide here are important to improve ice accretion modelling techniques and establish a different approach to address the interaction between the cloud and the surfaces exposed to atmospheric icing.

To improve Antarctic sea-ice simulations and estimations, an ensemble-based Data Assimilation System for the Southern Ocean (DASSO) was developed based on a regional sea ice–ocean coupled model, which assimilates sea-ice thickness (SIT) together with sea-ice concentration (SIC) derived from satellites. To validate the performance of DASSO, experiments were conducted from 15 April to 14 October 2016. Generally, assimilating SIC and SIT can suppress the overestimation of sea ice in the model-free run. Besides considering uncertainties in the operational atmospheric forcing data, a covariance inflation procedure in data assimilation further improves the simulation of Antarctic sea ice, especially SIT. The results demonstrate the effectiveness of assimilating sea-ice observations in reconstructing the state of Antarctic sea ice, but also highlight the necessity of more reasonable error estimation for the background as well as the observation.