Glacier and ice sheet research is frequently justified on the basis of potential benefits to those communities that are most vulnerable to glacier change. In this glaciology research, funding priorities and communication to the broader public are strongly affected by the experiences and values of glaciology researchers. Using population data and newly available survey data from research organizations including glaciologists, we show that there is a substantial misalignment between the demographics of those who stand to benefit from glaciological research and those who produce glaciological knowledge. We discuss the potential negative consequences of this misalignment, which causes scientific research to be less effective, valuable and usable for communities. We conclude by outlining twenty evidence-based strategies that individuals and organizations can adopt to improve the recruitment and retention of a more representative group of scientists in glaciological research and encourage co-production with communities.

Calving glaciers respond quickly to atmospheric variability through ice dynamic adjustment. Particularly, single weather extremes may cause changes in ice-flow velocity and terminus position. Occasionally, this can lead to substantial event-driven mass loss at the ice front. We examine changes in terminus position, ice-flow velocity, and calving flux at the grounded lacustrine Schiaparelli Glacier in the Cordillera Darwin using geo-referenced time-lapse camera images and remote sensing data (Sentinel-1) from 2015 to 2022. Lake-level records, lake discharge measurements, and a coupled energy and mass balance model provide insight into the subglacial water discharge. We use downscaled reanalysis data (ERA5) to identify climate extremes and track land-falling atmospheric rivers to investigate the ice-dynamic response on possible atmospheric drivers.

Meltwater controls seasonal variations in ice-flow velocity, with an efficient subglacial drainage system developing during the warm season and propagating up-glacier. Calving accounts for 4.2% of the ice loss. Throughout the year, warm spells, wet spells, and landfalling atmospheric rivers promote calving. The progressive thinning of the ice destabilizes the terminus position, highlighting the positive feedback between glacier thinning, near-terminus ice-flow acceleration, and calving flux.

Citizen science and related engagement programmes have proliferated in recent years throughout the sciences but have been reasonably limited in the cryospheric sciences. In the semiarid Andes we at the Centro de Estudios Avanzados en Zonas Áridas have developed a range of initiatives together with the wider community and stakeholder institutions to improve our understanding of the role snow and ice play in headwater catchments. In this paper we reflect on ongoing engagement with communities living and working in and near study sites of cryospheric science in northern Chile as a strategy that can both strengthen the research being done and empower local communities.

Crevasses pose severe risks for mountaineers and field glaciologists. Smaller cracks between 0.5 and 2 m are still dangerous, but often not visible in medium resolution satellite imagery. If they are snow covered, they are completely undetectable by optical sensors. We set out to develop an approach to detect potentially crevassed areas by a minimum of geometric data, and to make the method generally applicable to glacier regions. On Vernagtferner, we compared a reference dataset of crevasses observed in high-resolution optical imagery with the curvature of the ice surface and the spatial gradients in driving stress. Both parameters can be derived from a digital surface model and a bedrock model, derived from ice thickness measurements. The correlation patterns show that crevasses preferably form in convex areas and in areas where the driving stress rapidly increases. This corresponds with the theory of crevasse formation. Although the method still misclassifies larger parts, the approach has the potential to define probable non-crevassed areas as well as to aid the planning of safe routes.

Raman spectra (532 nm diode laser) of gypsum, syngenite and görgeyite powders were studied from 273 to 83 K every −10 K. Although it was found that not all the Raman peaks are temperature-sensitive, the effect of temperature on the Raman signal of some bands was modelled. It was observed that the main bands of the three sulfates shift toward higher wavenumbers as temperature decreases. The strengthening of the bonds causes an increase in energy and, consequently, higher wavenumbers. However, –OH bands shift depending on the crystal structure they acquire as the temperature decreases, so not all –OH bands behave in the same way. Finally, the gypsum –OH bands change to more ordered crystalline structures with decreasing temperature, so their bandwidths become narrower and sharper.

Over the past 70 years, many different components of the cryosphere have been imaged with a variety of radar systems using increasingly sophisticated processing techniques. These systems use various pulse lengths, signal frequencies and, in some cases, modulated signals. The increasing diversity of radar systems has created the potential for confusion due to the use of non-consistent terminology. Here we provide an overview of state-of-the-science radar technologies and suggest a simplified and unified terminology for use by the cryosphere community. We recommend a terminology that is target independent but specifies the characteristics of the signal. Following this recommendation, commercial impulse systems that penetrate the subsurface should be referred to as ground-penetrating radar (GPR), and pulse radars as radio-echo sounding (RES). Continuous-wave (CW) radar systems should be referred to as ground-penetrating CW radars. We further suggest any additional characterisation of the system be expressed using descriptors that specify the platform it is mounted on (e.g. airborne) or the frequency range (e.g. HF (high frequency)) or modulation (e.g. FM (frequency modulated)).

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.

On 30 September 2017, an Air France Airbus A380-800 suffered a failure of its fourth engine while over Greenland. This failure resulted in the loss of the engine fan hub, fan blades and surrounding structure. An initial search recovered 30 pieces of light debris, but the primary part of interest, a ~220 kg titanium fan hub, was not recovered because it had a different fall trajectory than the light debris, impacted into the ice-sheet's snow surface, and was quickly covered by drifting snow. Here we describe the methods used for the detection of the fan hub and details of the field campaigns. The search area included two crevasse fields of at least 50 snow-covered crevasses 1 to ~30 m wide with similar snow bridge thicknesses. After 21 months and six campaigns, using airborne synthetic aperture radar, ground-penetrating radar, transient electromagnetics and an autonomous vehicle to survey the crevasse fields, the fan hub was found within ~1 m of a crevasse at a depth of ~3.3 to 4 m and was excavated with shovels, chain saws, an electric winch, sleds and a gasoline heater, by workers using fall-arrest systems.

We reanalyzed mass balance records at Taku and Lemon Creek Glaciers to better understand the relative roles of hypsometry, local climate and dynamics as mass balance drivers. Over the 1946–2018 period, the cumulative mass balances diverged. Tidewater Taku Glacier advanced and gained mass at an average rate of +0.25 ± 0.28 m w.e. a–1, contrasting with retreat and mass loss of −0.60 ± 0.15 m w.e. a−1 at land-terminating Lemon Creek Glacier. The uniform influence of regional climate is demonstrated by strong correlations among annual mass balance and climate data. Regional warming trends forced similar statistically significant decreases in surface mass balance after 1989: −0.83 m w.e. a–1 at Taku Glacier and −0.81 m w.e. a–1 at Lemon Creek Glacier. Divergence in cumulative mass balance arises from differences in glacier hypsometry and local climate. Since 2013 negative mass balance and glacier-wide thinning prevailed at Taku Glacier. These changes initiated terminus retreat, which could increase dramatically if calving begins. The future mass balance trajectory of Taku Glacier hinges on dynamics, likely ending the historic dichotomy between these glaciers.

The accumulation region of Fedchenko Glacier represents an extensive snow reservoir in the Pamir Mountains feeding the longest glacier in Central Asia. Observed elevation changes indicate a continuous ice loss in the ablation region of Fedchenko Glacier since 1928, while the mass balance of the accumulation region is largely unknown. In this study, we show that accumulation varies considerably in the main accumulation basin, with accumulation rates up to 2400 mm w.e. a−1 in the West, decreasing to <1000 mm w.e. a−1 in the center, although the elevation difference is <200 m. The combination of snow/firn samples and ground-penetrating radar profiles suggests that this accumulation pattern is persistent during the recent past. The recent accumulation history is reconstructed from internal radar reflectors using a firn densification model and shows strong interannual variations, but near constant mean values since 2002. Modeling of trajectories, based on accumulation and glacier geometry, results in an estimate of the depth/age relation close to the main divide. This region provides one of the most suitable locations for retrieving climate information with temporal high resolution for the last millennium, with a potential to cover most of the Holocene in less detail.

We compared elastic moduli in polar firn derived from diving wave refraction seismic velocity analysis, firn-core density measurements and microstructure modelling based on firn-core data. The seismic data were obtained with a small electrodynamic vibrator source near Kohnen Station, East Antarctica. The analysis of diving waves resulted in velocity–depth profiles for different wave types (P-, SH- and SV-waves). Dynamic elastic moduli of firn were derived by combining P- and S-wave velocities and densities obtained from firn-core measurements. The structural finite-element method (FEM) was used to calculate the components of the elastic tensor from firn microstructure derived from X-ray tomography of firn-core samples at depths of 10, 42, 71 and 99 m, providing static elastic moduli. Shear and bulk moduli range from 0.39 to 2.42 GPa and 0.68 to 2.42 GPa, respectively. The elastic moduli from seismic observations and the structural FEM agree within 8.5% for the deepest achieved values at a depth of 71 m, and are within the uncertainty range. Our observations demonstrate that the elastic moduli of the firn can be consistently obtained from two independent methods which are based on dynamic (seismic) and static (tomography and FEM) observations, respectively, for deeper layers in the firn below ~10 m depth.

Fedchenko Glacier experienced a large thickness loss since the first scientific investigations in 1928. As the largest glacier in the Pamir Mountains, this glacier plays an important role for the regional glacier mass budget. We use a series of Global Navigation Satellite Systems observations from 2009 to 2016 and TanDEM-X elevation models from 2011 to 2016 to investigate recent elevation changes. Accounting for radar wave penetration minimizes biases in elevation that can otherwise reach up to 6 m in dry snow on Fedchenko Glacier, with mean values of 3–4 m in the high accumulation regions. The seasonal elevation changes reach up to ±5 m. The glacier surface elevation decreased along its entire length over multi-year periods. Thinning rates increased between 2000 and 2016 by a factor of 1.8 compared with 1928–2000, resulting in peak values of 1.5 m a−1. Even the highest accumulation basins above 5000 m elevation have been affected by glacier thinning with change rates between −0.2 and −0.4 m a−1 from 2009 to 2016. The estimated glacier-wide mass-balance rates are −0.27 ± 0.05 m w.e. a−1 for 2000 to 2011 and −0.51 ± 0.04 m w.e. a−1 between 2011 and 2016.

The non-linearly viscous ice flow in the vicinity of a glacier terminus, an observation region, depends crucially on the upstream flow as well as on the local surface and bed conditions. The former requires a likely complex solution of the balance laws and boundary conditions for the complete glacier. However, if the profile and downstream surface tangential velocity in the observation region are measured at an observation time t = 0, and a two-dimensional flow approximation is satisfactory, the complete stress and velocity fields satisfying local reduced model equations in the observation region at time t = 0 can be determined by asymptotic expansions in upstream distance from the (moving) terminus. Thus the full strain-rate and stress tensors are determined without prescribing the basal conditions. The terminus velocity is determined in terms of the net accumulation or melt flux and surface velocity at the terminus, with bounds for advance or retreat. The analysis and illustration are presented for a plane flow approximation.

The flow of Glacier d’Argentière, French Alps, has been extensively studied by glaciologists for >100 years, but high-temporal- and high-spatial-resolution monitoring of its dynamics has rarely been carried out. In this paper, the terminal part of the main glacier was selected to undertake such a survey. The 2 month experiment (September–November 2013) included remote sensing to map the velocity field of the whole area, and a network of in situ measurements to record the ice velocity at the surface as well as at the base of the glacier with high precision and at high temporal resolution. These monitoring methods are combined to produce a comprehensive and reliable picture of the ice flow. Basal and surface in situ measurements are used to precisely monitor the temporal variability of the ice flow at a section of the glacier where crevassed areas and non-cracked ice coexist. The achieved results are an innovative and efficient tool to investigate the acceleration/deceleration pattern heterogeneities in ice flow that occur after heavy rainfall.

The cone penetration test has been successfully used to classify soil for ∼100 years. However, it has received only limited contemporary use in frozen geomaterials. The historical and contemporary use of the cone penetration test in various frozen geomaterials is considered here and contemporary data from recent work in polar snow are examined. It is probable that many material physical properties (e.g. density, strength and microstructure) can be obtained directly from cone penetration testing. It appears under-utilized as a contemporary scientific and engineering investigative tool in frozen geomaterials.

Despite their importance as freshwater reservoirs for downstream river systems, few glaciers in central Chile have been comprehensively surveyed. This study presents ground-penetrating radar (GPR) and field-based observations for characterizing the englacial and basal conditions of Glaciar Olivares Alfa (33°110 S, 70°130 W), central Chilean Andes. Using a 50 MHz radar mounted onto a helicopter platform, data were collected covering large portions of the glacier accumulation and ablation zones. The radar data revealed boundaries of a temperate-ice layer at the base of the eastern body of Glaciar Olivares Alfa which appears to be covered by colder ice that extends throughout large parts of the glacier. The thickness of the temperate ice layer is highly variable across the glacier, being on average 40% of the total ice thickness. Radar data analyses reveal regions of cold ice at the bottom/base of the glacier and also patterns of highly saturated sediments beneath the glacier. Using GPR data, this study represents the most exhaustive analysis of glacier ice structure performed in the central Chilean Andes. The results will enable improved estimations of the glacier’s mass balance and ice dynamics, helping us to understand its further development and its impact on water availability.

Marine ice is an important component of ice shelves in Antarctica. It accretes in substantial amounts at weak points and below ice shelves. It is likely to exhibit peculiar rheological properties, which are crucial to understanding its potential role in stabilizing ice-shelf flow. Due to its location and consolidation processes, marine ice can present a variety of textures which are likely to influence its rheological properties. We present a new dataset of unconfined uniaxial compression experiments on folded marine ice samples that have been cut at various angles to the folds. Texture and fabric analyses are described ‘before’ and ‘after’ the deformation experiment. It is shown that, in the given stress configuration, the geometry of the anisotropy controls the rheological behaviour of the marine ice. During secondary creep, folded marine ice is harder to deform than weakly textured ice when compressed parallel or perpendicular to the folds' hinges, while the reverse is true for ice compressed at 45°. The observed range of values for the n exponent in Glen's flow law is between 2.1 and 4.1. Surprisingly, we see that tertiary creep tends to develop at a higher total strain than for randomly oriented impurity-free meteoric ice.

Glacier dynamics are inextricably linked to the basal conditions of glaciers. Seismic reflection methods can image the glacier bed under certain conditions. However, where a seismically thin layer of material is present at the bed, traditional analyses may fail to fully characterize bed properties. We use a targeted full-waveform inversion algorithm to quantify the basal-layer parameters of a mountain glacier: thickness (d), P-wave velocity (α) and density (ρ). We simultaneously invert for the seismic quality factor (Q) of the bulk glacier ice. The inversion seeks to minimize the difference between the data and a one-dimensional reflectivity algorithm using a gradient-based search with starting values initialized from a Monte Carlo scheme. We test the inversion algorithm on four basal layer synthetic data traces with 5% added Gaussian noise. The inversion retrieved thin-layer parameters within 10% of synthetic test parameters with the exception of seismic Q. For the seismic dataset from Bench Glacier, Alaska, USA, inversion results indicate a thin basal layer of debris-rich ice within the study area having mean velocity 4000 ± 700 m s–1, density 1900 ± 200 kg m–3 and thickness 6 ± 1.5 m.

To better understand the risk of disasters due to glacial lake outburst floods (GLOFs), we synthetically analyze the spatial distribution and evolution of moraine-dammed lakes and potentially dangerous glacial lakes (PDGLs) in the Chinese Himalaya. Our county-based assessment of GLOF disaster risk combines PDGL outburst hazard, regional exposure, vulnerability of exposed elements and adaptation capability (risk management) using the analytic hierarchy process. We synthetically analyze the disaster risk using the weighted comprehensive method. Remote-sensing data show there are 329 moraine-dammed lakes (>0.02 km2; total area 125.43 km2) in the Chinese Himalaya, of which 116 (total area 49.49 km2) are identified as PDGLs. The zones at highest risk of GLOF disaster are mainly located in Nyalam, Tingri, Dinggyê, Lhozhag, Kangmar and Zhongba, in the mid-eastern Himalaya. Lowest-risk zones are located in the eastern Himalaya. On the county scale, Lhozhag and Lhunze have the highest hazard degrees and exposure, while Zhongba and Zando have the highest degree of vulnerability and lowest adaptation capacity. Our regionalization results for GLOF disaster risk are consistent with the distribution of historical disaster sites across the Chinese Himalaya.

The Wind River Range in Wyoming, USA, contains the largest concentration of glacial mass in the Rocky Mountains of the contiguous USA. Despite this distinction, only a few field or remotely sensed studies providing glacier volume changes have been published. The current study focuses on Continental Glacier located on the northern end of the range and uses two field datasets (high-accuracy GPS surface elevation points and ice-penetrating radar transects of the glacier bed) to create a three-dimensional model of glacier volume. Current surface elevations are compared with historical elevation data to calculate surface elevation change over time. An average thinning rate of 13.8 ± 7.8 m (0.30 ± 0.17 m a–1) between 1966 and 2012 was found. Surface elevation change rates varied across the glacier, ranging from +0.30 to –0.98 m a–1. Taking into account variable melt rates across the glacier, along with a glacial volume of 72.1 × 106 ± 10.8 × 106 m3, we estimate that Continental Glacier will be reduced in volume by 43% over the next 100 years and will disappear completely over the next 300–400 years, if current climatic conditions persist.