Search

The density structure of firn has implications for hydrological and climate modelling, and ice-shelf stability. The structure of firn can be evaluated from depth models of seismic velocity, widely obtained with Herglotz–Wiechert inversion (HWI), an approach that considers the slowness of refracted seismic arrivals. However, HWI is strictly appropriate only for steady-state firn profiles and the inversion accuracy can be compromised where firn contains ice layers. In these cases, full waveform inversion (FWI) may yield more success than HWI. FWI extends HWI capabilities by considering the full seismic waveform and incorporates reflected arrivals. Using synthetic firn density profiles, assuming both steady- and non-steady-state accumulation, we show that FWI outperforms HWI for detecting ice slab boundaries (5–80 m thick, 5–80 m deep) and velocity anomalies within firn. FWI can detect slabs thicker than one wavelength (here, 20 m, assuming a maximum frequency of 60 Hz) but requires the starting velocity model to be accurate to ±2.5%. We recommend for field practice that the shallowest layers of velocity models are constrained with ground-truth data. Nonetheless, FWI shows advantages over established methods, and should be considered when the characterisation of firn ice slabs is the goal of the seismic survey.

A precise knowledge of landfast sea-ice (fast-ice) thickness is relevant to many different disciplines. Sea Ice Monitoring Stations (SIMS) are used to measure time series of fast-ice thickness at a location. SIMS measure ice and ocean temperature via thermistor strings with many different methods for extracting sea-ice thickness from temperature existing. This study investigates: if thickness results from temperature recorded by SIMS of different designs, and analysed with different methods are comparable; which methods are recommended for their robustness, precision and accuracy and how they compare to independent thickness measurements; how otherwise unuseable data can be salvaged through specific SIMS designs. We present an analysis of fast-ice thickness calculated from SIMS deployed in McMurdo Sound, Antarctica and in the Chukchi Sea near Utqiaġvik, Alaska, over two decades. We find that median thicknesses derived by different methods agree within 1 ± 1.5 cm for McMurdo Sound and 2 ± 3 cm for Utqiaġvik. Thus, it is possible to confidently compare data collected from different stations and analysed with different methods. The vertical gradient of sea-ice temperature gives the best results for fast-ice thickness during the growth season and including standard resistors in a thermistor string can reduce potential data loss due to noise.

Links between proglacial lakes and glacier dynamics are poorly understood but are necessary to predict how mountain glaciers will react to a warmer, wetter climate, where such lakes are expected to increase both in number and volume. Here, we examine a long-term (~120 year) record of terminus retreat, thinning and surface velocities from in-situ and remote sensing observations at the terminus of Kaskawulsh Glacier, Yukon, Canada, and determine the impact of a local proglacial hydrological reorganisation on glacier dynamics. After an initial deceleration during the late 1990s, terminus velocities increased at a rate of 3 m a−2 from 2000–12, while proglacial Slims Lake area increased simultaneously. The rapid drainage of the lake in May 2016 substantially altered the velocity profile, decreasing annual velocities by 48% within 3 km of the terminus between 2015 and 2021, at an average rate of ~ 12.5 m a−2. A key cause of the rapid drop in glacier motion was a reduction in flotation of the lower part of the glacier terminus after lake drainage. This has important implications for glacier dynamics and provides one of the first assessments of the impacts of a rapid proglacial lake drainage event on local terminus velocities.

The Antarctic Peninsula Ice Sheet (APIS) has become a significant contributor to sea-level rise over recent decades. Accurately estimating the ice discharge from the outlet glaciers of the APIS is crucial to quantify the mass balance of the Antarctic Peninsula. We here compute the ice discharge from the outlet glaciers of the APIS north of 70${^\circ }$S for the five most widely used ice-thickness reconstructions, using a common surface velocity field and a common set of flux gates, so the differences in ice discharge can be solely attributed to the differences in ice thickness at the flux gates. The total volumetric ice discharge for 2015–2017 ranges within 45–141 km3 a−1, depending on the ice-thickness model, with a mean of 87 ± 44 km3 a−1. The substantial differences between the ice-discharge results, and a multi-model normalized root-mean-squared deviation of 0.91 for the whole data set, reveal large differences and inconsistencies between the ice-thickness models, giving an indication of the large uncertainty in the current ice-discharge estimates for the APIS. This manifests a fundamental problem of the region: the scarcity of appropriate ice-thickness measurements and the difficulty of the current models to reconstruct the ice-thickness distribution in this complex region.

Snow ice (SI) forms from freezing wet snow, known as slush, and contributes to the thickness of level and brash ice. However, the mechanism of snow-slush-snow ice transformation has not been extensively investigated to date, despite the difference in the freezing rate of slush in comparison with water is important for estimating the ice thickness. In this study, we examined the growth of initial congelation ice (CI) and snow ice (SI) in a fresh water tank exposed to outdoor weather conditions in Luleå, northern Sweden. The tank of size 1.8 × 0.65 × 1.2 m in length, width and height was divided into two compartments to facilitate the simultaneous growth of CI and SI. A total of 12 experiments were conducted in the years 2021 and 2022. The transformation from slush to snow ice was achieved by submerging various amounts of snow in the compartments. It was observed that approximately 35% of the initial snow transformed into SI. Snow ice grew 4 mm°C−0.5 d−0.5 faster than congelation ice. The CI growth under SI was 1 mm°C−0.5 d−0.5 slower than the CI growth under CI. This study provides valuable insights for modelling snow-slush-snow ice transformation and designing future laboratory-scale experiments.

Theory and observation show that glacier-flow regimes characterized by high basal slip enhance the projection of topographic detail to the surface, motivating this investigation into the efficacy of using glacier surges to improve bed estimation. Here we adapt a Bayesian inversion scheme and apply it to real and synthetic data as a proof of concept. Synthetic tests show a reduction in mean RMSE between true and inferred beds by more than half, and an increase in the mean correlation coefficient of ~0.5, when data from slip- versus deformation-dominated regimes are used. Multi-epoch inversions, which partition slip- and deformation-dominated regimes, are shown to outperform inversions that average over these flow regimes thereby squandering information. Tests with real data from a surging glacier in Yukon, Canada, corroborate these results, while highlighting the challenges of limited or inconsistent data. With the growing torrent of satellite-based observations, fast-flow events such as glacier surges offer potential to improve bed estimation for some of the world's most dynamic glaciers.

The density structure of firn has implications for hydrological and climate modelling, and ice-shelf stability. The structure of firn can be evaluated from depth models of seismic velocity, widely obtained with Herglotz–Wiechert inversion (HWI), an approach that considers the slowness of refracted seismic arrivals. However, HWI is strictly appropriate only for steady-state firn profiles and the inversion accuracy can be compromised where firn contains ice layers. In these cases, full waveform inversion (FWI) may yield more success than HWI. FWI extends HWI capabilities by considering the full seismic waveform and incorporates reflected arrivals. Using synthetic firn density profiles, assuming both steady- and non-steady-state accumulation, we show that FWI outperforms HWI for detecting ice slab boundaries (5–80 m thick, 5–80 m deep) and velocity anomalies within firn. FWI can detect slabs thicker than one wavelength (here, 20 m, assuming a maximum frequency of 60 Hz) but requires the starting velocity model to be accurate to ±2.5%. We recommend for field practice that the shallowest layers of velocity models are constrained with ground-truth data. Nonetheless, FWI shows advantages over established methods, and should be considered when the characterisation of firn ice slabs is the goal of the seismic survey. [ABSTRACT FROM AUTHOR]

During the second half of the 20th century, the Antarctic Peninsula region has undergone a long and sustained warming period, followed by a shorter but also sustained cooling period, and then a very recent return to warming conditions. All of these have profoundly impacted the glaciers peripheral to the Antarctic Peninsula. This paper focuses on the analysis of the surface mass balance monitoring of such glaciers by the glaciological method, complemented by the analysis of mass-balance estimates by geodetic methods, as well as frontal ablation estimates. We aim to summarize the current knowledge and outline the main challenges faced by investigating the mass balance of such peripheral glaciers and their current contribution to sea-level rise. [ABSTRACT FROM AUTHOR]

The debris that covers the ablation areas of high-elevation debris-covered glaciers contributes to the distinctive features and processes occurring both on and within such glaciers. Despite recent advances, knowledge of the subsurface environments of high-elevation debris-covered glaciers is still extremely limited. In particular, targeted field-based data are needed to parameterise and refine the projections of these glaciers in numerical models. Here, we outline the current understanding of the internal properties of high-elevation debris-covered glaciers based on direct field-based methods and suggest future research directions for field-based studies. [ABSTRACT FROM AUTHOR]

Links between proglacial lakes and glacier dynamics are poorly understood but are necessary to predict how mountain glaciers will react to a warmer, wetter climate, where such lakes are expected to increase both in number and volume. Here, we examine a long-term (~120 year) record of terminus retreat, thinning and surface velocities from in-situ and remote sensing observations at the terminus of Kaskawulsh Glacier, Yukon, Canada, and determine the impact of a local proglacial hydrological reorganisation on glacier dynamics. After an initial deceleration during the late 1990s, terminus velocities increased at a rate of 3 m a−2 from 2000–12, while proglacial Slims Lake area increased simultaneously. The rapid drainage of the lake in May 2016 substantially altered the velocity profile, decreasing annual velocities by 48% within 3 km of the terminus between 2015 and 2021, at an average rate of ~ 12.5 m a−2. A key cause of the rapid drop in glacier motion was a reduction in flotation of the lower part of the glacier terminus after lake drainage. This has important implications for glacier dynamics and provides one of the first assessments of the impacts of a rapid proglacial lake drainage event on local terminus velocities. [ABSTRACT FROM AUTHOR]

A precise knowledge of landfast sea-ice (fast-ice) thickness is relevant to many different disciplines. Sea Ice Monitoring Stations (SIMS) are used to measure time series of fast-ice thickness at a location. SIMS measure ice and ocean temperature via thermistor strings with many different methods for extracting sea-ice thickness from temperature existing. This study investigates: if thickness results from temperature recorded by SIMS of different designs, and analysed with different methods are comparable; which methods are recommended for their robustness, precision and accuracy and how they compare to independent thickness measurements; how otherwise unuseable data can be salvaged through specific SIMS designs. We present an analysis of fast-ice thickness calculated from SIMS deployed in McMurdo Sound, Antarctica and in the Chukchi Sea near Utqiaġvik, Alaska, over two decades. We find that median thicknesses derived by different methods agree within 1 ± 1.5 cm for McMurdo Sound and 2 ± 3 cm for Utqiaġvik. Thus, it is possible to confidently compare data collected from different stations and analysed with different methods. The vertical gradient of sea-ice temperature gives the best results for fast-ice thickness during the growth season and including standard resistors in a thermistor string can reduce potential data loss due to noise. [ABSTRACT FROM AUTHOR]

Snow ice (SI) forms from freezing wet snow, known as slush, and contributes to the thickness of level and brash ice. However, the mechanism of snow-slush-snow ice transformation has not been extensively investigated to date, despite the difference in the freezing rate of slush in comparison with water is important for estimating the ice thickness. In this study, we examined the growth of initial congelation ice (CI) and snow ice (SI) in a fresh water tank exposed to outdoor weather conditions in Luleå, northern Sweden. The tank of size 1.8 × 0.65 × 1.2 m in length, width and height was divided into two compartments to facilitate the simultaneous growth of CI and SI. A total of 12 experiments were conducted in the years 2021 and 2022. The transformation from slush to snow ice was achieved by submerging various amounts of snow in the compartments. It was observed that approximately 35% of the initial snow transformed into SI. Snow ice grew 4 mm°C−0.5 d−0.5 faster than congelation ice. The CI growth under SI was 1 mm°C−0.5 d−0.5 slower than the CI growth under CI. This study provides valuable insights for modelling snow-slush-snow ice transformation and designing future laboratory-scale experiments. [ABSTRACT FROM AUTHOR]

The Antarctic Peninsula Ice Sheet (APIS) has become a significant contributor to sea-level rise over recent decades. Accurately estimating the ice discharge from the outlet glaciers of the APIS is crucial to quantify the mass balance of the Antarctic Peninsula. We here compute the ice discharge from the outlet glaciers of the APIS north of 70 ${^\circ }$ S for the five most widely used ice-thickness reconstructions, using a common surface velocity field and a common set of flux gates, so the differences in ice discharge can be solely attributed to the differences in ice thickness at the flux gates. The total volumetric ice discharge for 2015–2017 ranges within 45–141 km3 a−1, depending on the ice-thickness model, with a mean of 87 ± 44 km3 a−1. The substantial differences between the ice-discharge results, and a multi-model normalized root-mean-squared deviation of 0.91 for the whole data set, reveal large differences and inconsistencies between the ice-thickness models, giving an indication of the large uncertainty in the current ice-discharge estimates for the APIS. This manifests a fundamental problem of the region: the scarcity of appropriate ice-thickness measurements and the difficulty of the current models to reconstruct the ice-thickness distribution in this complex region. [ABSTRACT FROM AUTHOR]

Theory and observation show that glacier-flow regimes characterized by high basal slip enhance the projection of topographic detail to the surface, motivating this investigation into the efficacy of using glacier surges to improve bed estimation. Here we adapt a Bayesian inversion scheme and apply it to real and synthetic data as a proof of concept. Synthetic tests show a reduction in mean RMSE between true and inferred beds by more than half, and an increase in the mean correlation coefficient of ~0.5, when data from slip- versus deformation-dominated regimes are used. Multi-epoch inversions, which partition slip- and deformation-dominated regimes, are shown to outperform inversions that average over these flow regimes thereby squandering information. Tests with real data from a surging glacier in Yukon, Canada, corroborate these results, while highlighting the challenges of limited or inconsistent data. With the growing torrent of satellite-based observations, fast-flow events such as glacier surges offer potential to improve bed estimation for some of the world's most dynamic glaciers. [ABSTRACT FROM AUTHOR]

Over the last 25 years, radiowave detection of neutrino-generated signals, using cold polar ice as the neutrino target, has emerged as perhaps the most promising technique for detection of extragalactic ultra-high energy neutrinos (corresponding to neutrino energies in excess of 0.01 Joules, or 1017 electron volts). During the summer of 2021 and in tandem with the initial deployment of the Radio Neutrino Observatory in Greenland (RNO-G), we conducted radioglaciological measurements at Summit Station, Greenland to refine our understanding of the ice target. We report the result of one such measurement, the radio-frequency electric field attenuation length $L_\alpha$. We find an approximately linear dependence of $L_\alpha$ on frequency with the best fit of the average field attenuation for the upper 1500 m of ice: $\langle L_\alpha \rangle = ( ( 1154 \pm 121) - ( 0.81 \pm 0.14) \, ( \nu /{\rm MHz}) ) \,{\rm m}$ for frequencies ν ∈ [145 − 350] MHz.

During the second half of the 20th century, the Antarctic Peninsula region has undergone a long and sustained warming period, followed by a shorter but also sustained cooling period, and then a very recent return to warming conditions. All of these have profoundly impacted the glaciers peripheral to the Antarctic Peninsula. This paper focuses on the analysis of the surface mass balance monitoring of such glaciers by the glaciological method, complemented by the analysis of mass-balance estimates by geodetic methods, as well as frontal ablation estimates. We aim to summarize the current knowledge and outline the main challenges faced by investigating the mass balance of such peripheral glaciers and their current contribution to sea-level rise.

The debris that covers the ablation areas of high-elevation debris-covered glaciers contributes to the distinctive features and processes occurring both on and within such glaciers. Despite recent advances, knowledge of the subsurface environments of high-elevation debris-covered glaciers is still extremely limited. In particular, targeted field-based data are needed to parameterise and refine the projections of these glaciers in numerical models. Here, we outline the current understanding of the internal properties of high-elevation debris-covered glaciers based on direct field-based methods and suggest future research directions for field-based studies.

The earliest airborne geophysical campaigns over Antarctica and Greenland in the 1960s and 1970s collected ice penetrating radar data on 35 mm optical film. Early subglacial topographic and englacial stratigraphic analyses of these data were foundational to the field of radioglaciology. Recent efforts to digitize and release these data have resulted in geometric and ice-thickness analysis that constrain subsurface change over multiple decades but stop short of radiometric interpretation. The primary challenge for radiometric analysis is the poorly-characterized compression applied to Z-scope records and the sparse sampling of A-scope records. Here, we demonstrate the information richness and radiometric interpretability of Z-scope records. Z-scope pixels have uncalibrated fast-time, slow-time, and intensity scales. We develop approaches for mapping each of these scales to physical units (microseconds, seconds, and signal to noise ratio). We then demonstrate the application of this calibration and analysis approach to a flight in the interior of East Antarctica with subglacial lakes and to a reflight of an East Antarctic ice shelf that was observed by both archival and modern radar. These results demonstrate the potential use of Z-scope signals to extend the baseline of radiometric observations of the subsurface by decades.

Over the last 25 years, radiowave detection of neutrino-generated signals, using cold polar ice as the neutrino target, has emerged as perhaps the most promising technique for detection of extragalactic ultra-high energy neutrinos (corresponding to neutrino energies in excess of 0.01 Joules, or 1017 electron volts). During the summer of 2021 and in tandem with the initial deployment of the Radio Neutrino Observatory in Greenland (RNO-G), we conducted radioglaciological measurements at Summit Station, Greenland to refine our understanding of the ice target. We report the result of one such measurement, the radio-frequency electric field attenuation length $L_\alpha$. We find an approximately linear dependence of $L_\alpha$ on frequency with the best fit of the average field attenuation for the upper 1500 m of ice: $\langle L_\alpha \rangle = ((1154 \pm 121) - (0.81 \pm 0.14) \, (\nu /{\rm MHz})) \,{\rm m}$ for frequencies ν ∈ [145 − 350] MHz. [ABSTRACT FROM AUTHOR]

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.

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.

During the second half of the 20th century, the Antarctic Peninsula region has undergone a long and sustained warming period, followed by a shorter but also sustained cooling period, and then a very recent return to warming conditions. All of these have profoundly impacted the glaciers peripheral to the Antarctic Peninsula. This paper focuses on the analysis of the surface mass balance monitoring of such glaciers by the glaciological method, complemented by the analysis of mass-balance estimates by geodetic methods, as well as frontal ablation estimates. We aim to summarize the current knowledge and outline the main challenges faced by investigating the mass balance of such peripheral glaciers and their current contribution to sea-level rise. [ABSTRACT FROM AUTHOR]

The debris that covers the ablation areas of high-elevation debris-covered glaciers contributes to the distinctive features and processes occurring both on and within such glaciers. Despite recent advances, knowledge of the subsurface environments of high-elevation debris-covered glaciers is still extremely limited. In particular, targeted field-based data are needed to parameterise and refine the projections of these glaciers in numerical models. Here, we outline the current understanding of the internal properties of high-elevation debris-covered glaciers based on direct field-based methods and suggest future research directions for field-based studies. [ABSTRACT FROM AUTHOR]

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

The earliest airborne geophysical campaigns over Antarctica and Greenland in the 1960s and 1970s collected ice penetrating radar data on 35 mm optical film. Early subglacial topographic and englacial stratigraphic analyses of these data were foundational to the field of radioglaciology. Recent efforts to digitize and release these data have resulted in geometric and ice-thickness analysis that constrain subsurface change over multiple decades but stop short of radiometric interpretation. The primary challenge for radiometric analysis is the poorly-characterized compression applied to Z-scope records and the sparse sampling of A-scope records. Here, we demonstrate the information richness and radiometric interpretability of Z-scope records. Z-scope pixels have uncalibrated fast-time, slow-time, and intensity scales. We develop approaches for mapping each of these scales to physical units (microseconds, seconds, and signal to noise ratio). We then demonstrate the application of this calibration and analysis approach to a flight in the interior of East Antarctica with subglacial lakes and to a reflight of an East Antarctic ice shelf that was observed by both archival and modern radar. These results demonstrate the potential use of Z-scope signals to extend the baseline of radiometric observations of the subsurface by decades. [ABSTRACT FROM AUTHOR]

Knowledge of glacier volume is crucial for ice flow modelling and predicting the impacts of climate change on glaciers. Rugged terrain, harsh weather conditions and logistic costs limit field-based ice thickness observations in the Himalaya. Remote-sensing applications, together with mathematical models, provide alternative techniques for glacier ice thickness and volume estimation. The objective of the present research is to assess the application of artificial neural network (ANN) modelling coupled with remote-sensing techniques to estimate ice thickness on individual glaciers with direct field measurements. We have developed two ANN models and estimated the ice thickness of Chhota Shigri Glacier (western Himalaya) on ten transverse cross sections and two longitudinal sections. The ANN model estimates agree well with ice thickness measurements from a ground-penetrating radar, available for five transverse cross sections on Chhota Shigri Glacier. The overall root mean square errors of the two ANN models are 24 and 13 m and the mean bias errors are ±13 and ±6 m, respectively, which are significantly lower than for other available models. The estimated mean ice thickness and volume for Chhota Shigri Glacier are 109 ± 17 m and 1.69 ± 0.26 km3, respectively.

Glaciers in the Tien Shan mountains contribute considerably to the fresh water used for irrigation, households and energy supply in the dry lowland areas of Kyrgyzstan and its neighbouring countries. To date, reconstructions of the current ice volume and ice thickness distribution remain scarce, and accurate data are largely lacking at the local scale. Here, we present a detailed ice thickness distribution of Ashu-Tor, Bordu, Golubin and Kara-Batkak glaciers derived from radio-echo sounding measurements and modelling. All the ice thickness measurements are used to calibrate three individual models to estimate the ice thickness in inaccessible areas. A cross-validation between modelled and measured ice thickness for a subset of the data is performed to attribute a weight to every model and to assemble a final composite ice thickness distribution for every glacier. Results reveal the thickest ice on Ashu-Tor glacier with values up to 201 ± 12 m. The ice thickness measurements and distributions are also compared with estimates composed without the use of in situ data. These estimates approach the total ice volume well, but local ice thicknesses vary substantially.

Airborne electromagnetic induction sensors have demonstrated their extensive capacities to measure sea-ice thickness distributions. However, biases can emerge when comparing these 1-D measurements to a broader 2-D regional scale due to the spatial anisotropy inherent to sea-ice cover. Automated processing of available sea-ice maps could significantly ease the decision on how to set up an optimised flight pattern, which would result in representative ice thickness numbers for the region. In this study, first we investigate the extent to which the sea-ice anisotropy can influence the representativeness of an airborne survey compared to the regional situation. Second, we propose a method to process sea-ice maps prior to flights to help preparing the most representative flight plan possible for the local area. The method is based on automated segmentation of radar satellite images and extensive simulation of flight transects over the image. The spatial analysis of these transects enables for the identification of the most representative survey trajectories for the area. The method was applied for seven different synthetic aperture radar satellite images over Arctic sea ice north of Svalbard. The results indicate that the proposed method improved the representativeness of the airborne survey by identifying the most suitable transect over the ice pack.

Several global datasets of glacier thickness exist, but the number of observations from western Canada are sparse and spatially biased. To supplement these limited observations, we measured ice thickness with ice penetrating radar on five glaciers in the Columbia Mountains, Canada. Our radar surveys, when combined with previous surveys for two glaciers in the Rocky Mountains, total 182 km of transects that represent 34 672 point measurements of ice thickness. Our measurements are, on average, 38% thicker than previous surface inversion model estimates of glacier thickness. Using our measurements within a cross-validation scheme, we model ice thickness with the Open Global Glacier Model (OGGM) driven with recent observations of surface mass balance and glacier elevation. We calibrated OGGM ice thickness by optimizing the ice creep parameter in the model. The optimized OGGM yields an ice volume for Columbia Basin of 122.5 ± 22.4 km3 for the year 2000, which is 23% greater than the range of previous estimates. At current rates of glacier mass loss for this region, glaciers would disappear from the basin in about 65–80 years. Disappearance of these glaciers will negatively affect the basin's surface hydrology, freshwater availability and aquatic ecosystems.

Glacier ice thickness is crucial to quantifying water resources in mountain regions, and is an essential input for ice-flow models. Using a surface velocity inversion method, we combine ice thickness measurements with detailed surface elevation and velocity data, and derive ice thickness and volume estimates for the Monte Tronador glaciers, North Patagonian Andes. We test the dependence of the inversion model on surface slope by resampling glacier slopes using variable smoothing filter sizes of 16–720 m. While total glacier volumes do not differ considerably, ice thickness estimates show higher variability depending on filter size. Smaller (larger) smoothing scales give thinner (thicker) ice and higher (lower) noise in ice thickness distribution. A filter size of 300 m, equivalent to four times the mean ice thickness, produces a noise-free thickness distribution with an accuracy of 35 m. We estimate the volume of the Monte Tronador glaciers at 4.8 ± 2 km3 with a mean ice thickness of 75 m. Comparison of our results with earlier regional and global assessments shows that the quality of glacier inventories is a significant source of discrepancy. We show that including surface slope as an input parameter increases the accuracy of ice thickness distribution estimates.

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

NASA's Operation IceBridge mission flew over the Ross Sea, Antarctica (20 and 27 November 2013) and collected data with Airborne Topographic Mapper (ATM) and Digital Mapping System (DMS). Using the DMS and reflectivity of ATM L1B, leads are detected to define local sea level height. The total freeboard is then obtained and converted to ice thickness. The estimated mean sea-ice thickness values are found to be in the 0.48–0.99 m range. Along the N-S track, sea ice was thinner southward rather than northward of the fluxgate, resulting in two peaks of modal thickness: 0.35 m (south) and 0.7 m (north). This supports that new ice produced in coastal polynyas is transported northward by katabatic winds off the ice-shelf. The lowest (2%) elevation method used for freeboard retrieval for ICESat is also tested for ATM data. It is found that the lowest elevation method tends to overestimate freeboard, but mean values are less affected than mode values. Using mean thickness values of ICESat and ATM along the ‘fluxgate’, separating the shelf from the deep ocean, the exported ice volume at this ‘fluxgate’ is found to be higher during the ICESat years (2003–2008) than during the IceBridge year (2013).

Direct measurements of spatially distributed vertical strain within ice masses are scientifically valuable but challenging to acquire. We use manual marker tracking and automatic cross correlation between two repeat optical televiewer (OPTV) images of an ~100 m-long borehole at Derwael Ice Rise (DIR), Antarctica, to reconstruct discretised, vertical strain rate and velocity at millimetre resolution. The resulting profiles decay with depth, from −0.07 a−1 at the surface to ~−0.002 a−1 towards the base in strain and from −1.3 m a−1 at the surface to ~−0.5 m a−1 towards the base in velocity. Both profiles also show substantial local variability. Three coffee-can markers installed at different depths into adjacent boreholes record consistent strain rates and velocities, although averaged over longer depth ranges and subject to greater uncertainty. Measured strain-rate profiles generally compare closely with output from a 2-D ice-flow model, while the former additionally reveal substantial high-resolution variability. We conclude that repeat OPTV borehole logging represents an effective means of measuring distributed vertical strain at millimetre scale, revealing high-resolution variability along the uppermost ~100 m of DIR, Antarctica.

This paper provides an update and overview of the Center for Remote Sensing of Ice Sheets (CReSIS) radars and platforms, including representative results from these systems. CReSIS radar systems operate over a frequency range of 14–38 GHz. Each radar system's specific frequency band is driven by the required depth of signal penetration, measurement resolution, allocated frequency spectra, and antenna operating frequencies (often influenced by aircraft integration). We also highlight recent system advancements and future work, including (1) increasing system bandwidth; (2) miniaturizing radar hardware; and (3) increasing sensitivity. For platform development, we are developing smaller, easier to operate and less expensive unmanned aerial systems. Next-generation platforms will further expand accessibility to scientists with vertical takeoff and landing capabilities.

Despite widespread use of radio-echo sounding (RES) in glaciology and broad distribution of processed radar products, the glaciological community has no standard software for processing impulse RES data. Dependable, fast and collection-system/platform-independent processing flows could facilitate comparison between datasets and allow full utilization of large impulse RES data archives and new data. Here, we present ImpDAR, an open-source, cross-platform, impulse radar processor and interpreter, written primarily in Python. The utility of this software lies in its collection of established tools into a single, open-source framework. ImpDAR aims to provide a versatile standard that is accessible to radar-processing novices and useful to specialists. It can read data from common commercial ground-penetrating radars (GPRs) and some custom-built RES systems. It performs all the standard processing steps, including bandpass and horizontal filtering, time correction for antenna spacing, geolocation and migration. After processing data, ImpDAR's interpreter includes several plotting functions, digitization of reflecting horizons, calculation of reflector strength and export of interpreted layers. We demonstrate these capabilities on two datasets: deep (~3000 m depth) data collected with a custom (3 MHz) system in northeast Greenland and shallow (

Basal motion of ice sheets depends in part on the roughness and material properties of the subglacial bed and the occurrence of water. To date, basal motion represents one of the largest uncertainties in ice-flow models. It is that component of the total flow velocity that can change most rapidly and can, therefore, facilitate rapid variations in dynamic behaviour. In this study, we investigate the subglacial properties of the East Antarctic Ice Sheet by statistically analysing the roughness of the bed topography, inferred from radio-echo sounding measurements. We analyse two sets of roughness parameters, one derived in the spatial and the other in the spectral domain, with two roughness parameters each. This enables us to compare the suitability of the four roughness parameters to classify the subglacial landscapes below the ice sheet. We further investigate the relationship of the roughness parameters with observed surface flow velocity and modelled basal temperatures of the ice sheet. We find that one of the roughness parameters, the Hurst exponent derived in the spatial domain, coincides with the thermal condition at the base of the ice sheet for slow flow velocities and varies with flow velocity.

Several autonomous phase-sensitive radio-echo sounders (ApRES) were deployed at Greenland glaciers to investigate ice deformation. Different attenuation settings were tested and it was observed that, in the presence of clipping of the deramped ApRES signal, each setting produced a different result. Specifically, higher levels of clipping associated with lower attenuation produced an apparent linear increase of diurnal vertical cumulative displacement with depth, and obscured the visibility of the basal reflector in the return amplitude. An example with a synthetic deramped signal confirmed that these types of artifacts result from the introduction of harmonics from square-wave-like features introduced by clipping. Apparent linear increase of vertical displacement with depth occurs when the vertical position of a near-surface internal reflector changes in time. Artifacts in the return amplitude may obscure returns from internal reflectors and the basal reflector, making it difficult to detect thickness evolution of the ice and to correctly estimate vertical velocities. Variations in surface melt during ApRES deployments can substantially modulate the received signal strength on short timescales, and we therefore recommend using higher attenuator settings for deployments in such locations.

Russian (former Soviet) systematic studies of Antarctica by radio-echo sounding (RES) and ground-penetrating radar technique (GPR) were commenced in 1964. Since that time airborne RES surveys have covered about 5.5 × 106 km2 of the icy continent discovering remarkable geographic objects such as Subglacial Gamburtsev Mountains, and allowed studies of Filchner-Ronne Ice Shelf, Amery Ice Shelf and Lambert Glacier. Ground-based investigations during the 1990s and 2000s revealed the structure of the Lake Vostok area and surveyed along the Mirny to Vostok and Progress to Vostok traverse routes. GPR studies during the 2010s were to select the site for a new snow-runway at Mirny Station, with the resumption of the aviation after a 25 year hiatus.

The thickness of glaciers in High-Mountain Asia (HMA) is critical in determining when the ice reserve will be lost as these glaciers thin but is remarkably poorly known because very few measurements have been made. Through a series of ground-based and airborne field tests, we have adapted a low-frequency ice-penetrating radar developed originally for Antarctic over-snow surveys, for deployment as a helicopter-borne system to increase the number of measurements. The manoeuvrability provided by helicopters and the ability of our system to detect glacier beds through thick, dirty, temperate ice makes it well suited to increase greatly the sample of measurements available for calibrating ice thickness models on the regional and global scale. The Bedmap Himalayas radar-survey system can reduce the uncertainty in present-day ice volumes and therefore in projections of when HMA's river catchments will lose this hydrological buffer against drought.

Passive seismology is becoming increasingly popular for glacier/ice-sheet structure investigations in Polar regions. Single-station passive seismic methods including P-wave receiver functions (PRFs), horizontal-to-vertical spectral ratio (HVSR) and a recently proposed autocorrelation method have been used to retrieve glacier/ice-sheet structures. Despite their successful applications, analysis regarding their detection abilities in different glaciological environments has not been reported. In this study, we compare ice thicknesses and vp/vs ratios obtained from the three methods using data collected at GAMSEIS and POLENET/ANET seismic arrays in Antarctica. Ice thickness estimates derived from the three methods are found to be consistent. Comparisons conducted under various model setups, including those involving tiled layers and sedimentary layers, show that the effectiveness of the autocorrelation method is not superior to the PRF method for retrieving ice-sheet structures. The autocorrelation method however can complement other methods as it only requires a single component seismic record.

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

Knowledge of glacier volume is crucial for ice flow modelling and predicting the impacts of climate change on glaciers. Rugged terrain, harsh weather conditions and logistic costs limit field-based ice thickness observations in the Himalaya. Remote-sensing applications, together with mathematical models, provide alternative techniques for glacier ice thickness and volume estimation. The objective of the present research is to assess the application of artificial neural network (ANN) modelling coupled with remote-sensing techniques to estimate ice thickness on individual glaciers with direct field measurements. We have developed two ANN models and estimated the ice thickness of Chhota Shigri Glacier (western Himalaya) on ten transverse cross sections and two longitudinal sections. The ANN model estimates agree well with ice thickness measurements from a ground-penetrating radar, available for five transverse cross sections on Chhota Shigri Glacier. The overall root mean square errors of the two ANN models are 24 and 13 m and the mean bias errors are ±13 and ±6 m, respectively, which are significantly lower than for other available models. The estimated mean ice thickness and volume for Chhota Shigri Glacier are 109 ± 17 m and 1.69 ± 0.26 km3, respectively. [ABSTRACT FROM AUTHOR]

Crevasses and englacial fracture networks route meltwater from a glacier's surface to the subglacial drainage system and thus influence glacial hydraulics. However, rapid fracture growth may also lead to sudden and potentially hazardous structural failure of unstable glaciers and ice dams, rifting of ice shelves, or iceberg calving. Here, we use passive seismic recordings to investigate the englacial fracture network on Glacier de la Plaine Morte, Switzerland. Glacier dynamics and the drainage of an ice-marginal lake give rise to numerous icequakes, the majority of which generate dispersed, high-frequency Rayleigh waves. A wide distribution of events allows us to study azimuthal anisotropy between 10 and 30 Hz in order to extract englacial seismic velocities in regions of preferentially oriented crevasses. Beamforming applied to a 100-m-aperture array reveals azimuthal anisotropy of Rayleigh-wave phase velocities reaching a strength of 8% at high frequencies. In addition, we find that the fast direction of wave propagation coincides with the observed surface strike of the narrow crevasses. Forward modeling and inversion of dispersion curves suggest that the azimuthal anisotropy is induced by a 40-m-thick crevassed layer at the surface of the glacier with 8% anisotropy in shear-wave velocity.

The long-term warming on the Antarctic Peninsula in the second half of the 20th century prompted rapid retreat of glaciers on the peninsula and surrounding islands. Retreat accelerated until the beginning of the new millennium when the regional warming trend significantly decreased. The response of glaciers to the change in temperature trend has been observed around the northern part of the Antarctic Peninsula but the timing of the shift from the surface lowering to mass gain remains unclear. Using historical aerial photographs, DEMs and satellite altimeter data from ICESat, we estimate areal and surface elevation changes of two small ice caps in the northern part of James Ross Island over the last 39 years. The glacierized area on Lachman Crags decreased from 4.337 ± 0.037 to 3.581 ± 0.014 km2 (−17.4%) between 1979 and 2006 and then increased to 3.597 ± 0.047 km2 (0.4%) until 2016. Surface lowering observed on ice caps after 1979 continued at least until 2008 as indicated by the ICESat data. The change from the lowering trend to increase in glacier surface elevation probably occurred after the ablation season 2008/09, which ranks among the warmest summers in the north-eastern Antarctic Peninsula since the mid-20th century.

Arctic sea ice is going through a dramatic change in its extent and volume at an unprecedented rate. Sea-ice thickness (SIT) is a controlling geophysical variable that needs to be understood with greater accuracy. For the first time, a SIT-retrieval method that exclusively uses only airborne SIT data for training the empirical algorithm to retrieve SIT from Soil Moisture Ocean Salinity (SMOS) brightness temperature (TB) at different polarization is presented. A large amount of airborne SIT data has been used from various field campaigns in the Arctic conducted by different countries during 2011–15. The algorithm attempts to circumvent the issue related to discrimination between TB signatures of thin SIT versus low sea-ice concentration. The computed SIT has a rms error of 0.10 m, which seems reasonably good (as compared to the existing algorithms) for analysis at the used 25 km grid. This new SIT retrieval product is designed for direct operational application in ice prediction/climate models.

Traditionally, helicopter-borne ground-penetrating radar (GPR) systems are operated with a single pair of bistatic dipole antennas to measure the thickness of glaciers. We demonstrate numerically that the directivity of the radiation pattern of single airborne dipoles do not correspond to an ideal full-space solution if the antennas are employed at typical flight heights. These directionality effects can degrade the quality of the subsurface images significantly, when the GPR antennas are orientated unfavorably. Since an adjustment of the antenna orientation is impractical during flight, we have developed a novel dual-polarization helicopter-borne GPR system consisting of two orthogonal pairs of commercial antennas in broadside configuration. To overcome the image quality deficits of the individual channels, we apply a pseudo-scalar approach in which we combine the data of both polarizations. Results of helicopter-borne GPR surveys on two alpine glaciers in Switzerland reveal more coherent bedrock reflections in the summed data compared with single dipole pair profiles. Generally, the dual-polarization setup is more suitable than a single antenna systems, because it is more versatile and less prone to directional effects caused by the placement of the dipole antennas in relation to undulating subsurface reflectors.

Despite their high value and importance for various glaciological applications, detailed ice thickness measurements of alpine glaciers are still very limited. Knowledge of bedrock topography is essential for paleoglaciological studies. The Guliya ice cap located on the Tibetan Plateau is one of the highest and largest ice caps in mid-low latitude regions. A detailed ground-penetrating radar (GPR) survey was conducted on the Guliya ice cap in 2015 using 20 and 40 MHz frequency antennas. An empirical Bayesian kriging method was used for ice thickness interpolation and uncertainty assessment. GPR measurements revealed complex basal topography of the Guliya glacier with a maximum thickness of 371.12 ± 13 m. The internal reflections caused by changes in the dielectric properties were registered on the 40 MHz radargrams at the summit and were attributed to density variations. As a result of this fieldwork, one of the largest ice thickness datasets in High Mountain Asia was obtained. Guliya glacier elevation changes were assessed by differencing digital elevation models. The glacier gained mass from 2000 to 2015 with an average rate of 0.270 ± 0.11 m w.e. a−1 at the summit and 0.279 ± 0.11 m w.e. a−1 at the lower elevations.

Glaciers in the Tien Shan mountains contribute considerably to the fresh water used for irrigation, households and energy supply in the dry lowland areas of Kyrgyzstan and its neighbouring countries. To date, reconstructions of the current ice volume and ice thickness distribution remain scarce, and accurate data are largely lacking at the local scale. Here, we present a detailed ice thickness distribution of Ashu-Tor, Bordu, Golubin and Kara-Batkak glaciers derived from radio-echo sounding measurements and modelling. All the ice thickness measurements are used to calibrate three individual models to estimate the ice thickness in inaccessible areas. A cross-validation between modelled and measured ice thickness for a subset of the data is performed to attribute a weight to every model and to assemble a final composite ice thickness distribution for every glacier. Results reveal the thickest ice on Ashu-Tor glacier with values up to 201 ± 12 m. The ice thickness measurements and distributions are also compared with estimates composed without the use of in situ data. These estimates approach the total ice volume well, but local ice thicknesses vary substantially. [ABSTRACT FROM AUTHOR]

Airborne electromagnetic induction sensors have demonstrated their extensive capacities to measure sea-ice thickness distributions. However, biases can emerge when comparing these 1-D measurements to a broader 2-D regional scale due to the spatial anisotropy inherent to sea-ice cover. Automated processing of available sea-ice maps could significantly ease the decision on how to set up an optimised flight pattern, which would result in representative ice thickness numbers for the region. In this study, first we investigate the extent to which the sea-ice anisotropy can influence the representativeness of an airborne survey compared to the regional situation. Second, we propose a method to process sea-ice maps prior to flights to help preparing the most representative flight plan possible for the local area. The method is based on automated segmentation of radar satellite images and extensive simulation of flight transects over the image. The spatial analysis of these transects enables for the identification of the most representative survey trajectories for the area. The method was applied for seven different synthetic aperture radar satellite images over Arctic sea ice north of Svalbard. The results indicate that the proposed method improved the representativeness of the airborne survey by identifying the most suitable transect over the ice pack. [ABSTRACT FROM AUTHOR]