Digital elevation models (DEMs) from the spaceborne interferometric radar mission TanDEM-X hold a large potential for glacier change assessments. However, a bias is potentially introduced through the penetration of the X-band signal into snow and firn. To improve our understanding of radar penetration on glaciers, we compare DEMs derived from the almost synchronous acquisition of TanDEM-X and Pléiades optical stereo-images of Grosser Aletschgletscher in March 2021. We found that the elevation bias – averaged per elevation bin – can reach up to 4–8 m in the accumulation area, depending on post co-registration corrections. Concurrent in situ measurements (ground-penetrating radar, snow cores, snow pits) reveal that the signal is not obstructed by the last summer horizon but reaches into perennial firn. Because of volume scattering, the TanDEM-X surface is determined by the scattering phase centre and does not coincide with a specific firn layer. We show that the bias corresponds to more than half of the decadal ice loss rate. To minimize the radar penetration bias, we recommend to select DEMs from the same time of the year and over long observation periods. A correction of the radar penetration bias is recommended, especially when combining optical and TanDEM-X DEMs.

We present a quantitative ground-penetrating radar (GPR) analysis of the Lower Calderone Glacieret to highlight the recent evolution of one of the southernmost glacial systems in Europe. The Upper and Lower Calderone Glacierets are the last two perennial ice bodies in the Apennine Mountains (Italy), and their continuous monitoring is important for glaciology, hydrology and climate science. We applied a proprietary auto-picking algorithm to track reflections accurately and objectively within three pseudo-3-D GPR data sets that were acquired over the Lower Glacieret in different years. After the time-to-depth conversion, the basal reflections were projected onto the normal versors of a 3-D topographic model of the survey area, at the different GPR trace positions. We then applied an Adjusted Inverse-Distance Weighted Spatial Interpolation method to extrapolate the ice-bedrock interface within the areas not directly covered by the GPR profiles and compare it with the topographic surface to recover the glacieret volume. In this paper, we critically examine the accuracy of the reconstructed models, including possible issues related to the challenging survey areas, such as local artifacts in the interpolated interface caused by irregular GPR coverage. We further discuss the various advantages of the implemented procedure with respect to more traditional glacier monitoring techniques.

This article describes a series of steps to integrate multiple modes of archaeological mapping in arid and agricultural settings. We use the coastal region of Peru as a case study and share our recent field experience at Cerro San Isidro, a multicomponent hill site located in the agriculture-intensive and mid-elevation (about 500 m asl) Moro region of the Nepeña Valley. In June and July 2022, we spent eight weeks deploying a combination of drone aerial imagery, pedestrian GPS reconnaissance, and GPR survey to map the surface and subsurface features at the site and in the adjacent agricultural fields. Our efforts suggest that the ancient settlement extended over an area of at least 50 ha, well beyond the visible surface architecture. Using a multimodal approach to confirming the partial destruction of archaeological vestiges by modern agricultural encroachment is both time-effective and noninvasive. The article offers insights from our experience, including the sequence of field operations, technical troubleshooting, and the collection and integration of datasets. We discuss the methodological potential and implications of this combination of multimodal mapping and its deployment in coastal Peru, a region that, like many others in the world, is increasingly subject to rapid agricultural expansion and other anthropogenic developments.

We review successes and challenges from five recent subglacial bedrock drilling campaigns intended to find evidence for Antarctic Ice Sheet retreat during warm periods in the geologic past. Insights into times when the polar ice sheets were smaller than present serve as guiding information for modeling efforts that aim to predict the rate and magnitude of future sea level rise that would accompany major retreat of the Antarctic Ice Sheet. One method to provide direct evidence for the timing of deglaciations and minimum extent of prior ice sheets is to extract subglacial bedrock cores for cosmogenic nuclide analysis from beneath the modern ice sheet surface. Here we summarize the lessons learned from five field seasons tasked with obtaining bedrock cores from shallow depths (<120 m beneath ice surface) across West Antarctica since 2016. We focus our findings on drilling efforts and technology and geophysical surveys with ground-penetrating radar. Shallow subglacial drilling provides a high risk, high reward means to test for past instabilities of the Antarctic Ice Sheet, and we highlight key challenges and solutions to increase the likelihood of success for future subglacial drilling efforts in polar regions.

Ice shelves, which regulate ice flow from the Antarctic ice sheet towards the ocean, are shaped by spatiotemporal patterns of surface accumulation, surface/basal melt and ice dynamics. Therefore, an ice dynamic and accumulation history are imprinted in the internal ice stratigraphy, which can be imaged by radar in the form of internal reflection horizons (IRHs). Here, IRHs were derived from radar data combined across radar platforms (airborne and ground-based) in coastal eastern Dronning Maud Land (East Antarctica), comprising three ice rises and adjacent two ice shelves. To facilitate interpretation of dominant spatiotemporal patterns of processes shaping the local IRH geometry, traced IRHs are classified into three different types (laterally continuous, discontinuous or absent/IRH-free). Near-surface laterally continuous IRHs reveal local accumulation patterns, reflecting the mean easterly wind direction, and correlate with surface slopes. Areas of current and past increased ice flow and internal deformation are marked by discontinuous or IRH-free zones, and can inform about paleo ice-stream dynamics. The established IRH datasets extend continent-wide mapping efforts of IRHs to an important and climatically sensitive ice marginal region of Antarctica and are ready for integration into ice-flow models to improve predictions of Antarctic ice drainage.

We analyze the internal structure of two polythermal glaciers, Hurd and Johnsons, located on Livingston Island, Antarctica, using 200 and 750 MHz GPR data collected in 2003/04, 2008/09 and 2016/17 field campaigns. Based on the different permittivities of snow and ice, we determined the thickness distribution of the end-of winter snow cover and of the cold ice layer. Their knowledge is fundamental for mass balance and glacier dynamics studies due to the different densities and rheological properties of such media. The average measured thicknesses for the snow and cold ice layers (the latter including the snow layer) were of 1.44 ± 0.09 and 29.1 ± 1.5 m, and their corresponding maxima were of 2.45 ± 0.21 and 80.8 ± 2.5 m. GPR snow profiling allowed for extension of the coverage of the snow thickness survey, but added little information to that supplied by snow pits, stake readings and manual snow probing, because of the multiplicity of reflections within the seasonal snowpack caused by internal ice layers and lenses. The polythermal structure determined for Hurd Glacier fits into the so-called Scandinavian type, seldom reported for the Antarctic region.

Observations remain sparse for peripheral glaciers and ice caps in Greenland. Here, we present the results of a multi-frequency radar survey of Lyngmarksbræen Ice Cap in West Greenland conducted in April 2017. Radar measurements show thick ice of up to ~120 m in subglacial valleys associated with the largest outlet glaciers, while relatively thin ice cover the upper plateau ice divides, suggesting future vulnerability to ice cap fragmentation. At the time of the radar survey, Lyngmarksbræen Ice Cap had a total volume of 0.82 ± 0.1 km3. Measurements show a 1.5–2 m thick end-of-winter snowpack, and that firn is largely absent, signifying a prolonged period of negative mass balance for most of the ice cap. The thermal regime of Lyngmarksbræen Ice Cap is investigated through analysis of scattering observed along radar profiles. Results show that the ice cap is largely below the pressure melting point, but that temperate ice exists both in deep basal pockets and in shallow zones that some places extend from ~15 m depth and to the ice base. The distribution of shallow temperate ice appears unrelated to variations in ice thickness; instead we find a strong correlation to the presence of nearby surface crevasses.

Ground penetrating radar (GPR) has been extensively used in glaciology to infer glacier's ice thickness, liquid water content, water drainage pathways, and other properties. The interpretation of such GPR data is not always straightforward and for temperate glaciers, the signal is often affected by strong scattering and attenuation. It has often been suggested that such effects originate from englacial water inclusions, since water and ice have a large contrast in their di-electric permittivity. To investigate such effects quantitatively, we perform an extensive numerical modeling study of GPR signals. By exploring how different liquid water contents (LWC) and water-inclusions size affect the GPR signal, we show that their effects are much larger than the potential presence of a wet snowpack or a heterogeneous distribution of ice permittivity. In particularly, we show that the presence of such water inclusions is a necessary and sufficient condition for reproducing the typical characteristics of GPR data acquired in the field. Further, we find that for 25 MHz GPR antennas, a bulk LWC $\gtrsim$ 0.2%, associated with decimeters-scale water inclusions already limits bedrock detectability for ice thicknesses $\gtrsim 100$ m. Since these values are typical for Alpine glaciers, they clarify why the quality of GPR data is often poor in such environments.

Measuring surface mass-balance in the accumulation areas of glaciers is challenging because of the high spatial variability of snow accumulation and the difficulty of conducting annual field glaciological measurements. Here, we propose a method that can solve both these problems for many locations. Ground-penetrating radar measurements and firn cores extracted from a site in the French Alps were first used to reconstruct the topography of a buried end-of-summer snow horizon from a past year. Using these data and surface elevation observations from LiDAR and Global Navigation Satellite System instruments, we calculated the submergence velocities over the period between the buried horizon and more recent surface elevation observations. The differences between the changes in surface elevation and the submergence velocities were then used to calculate the annual surface mass-balances with an accuracy of ±0.34 m w.e. Assuming that the submergence velocities remain stable over several years, the surface mass-balance can be reconstructed for subsequent years from the differences in surface elevation alone. As opposed to the glaciological method that requires substantial fieldwork year after year to provide only point observations, this method, once submergence velocities have been calculated, requires only remote-sensing data to provide spatially distributed annual mass-balances in accumulation areas.

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)).

We assess the composition and geometry of four individual rock glaciers in Alaska, Wyoming and Colorado by measuring their radio wave speed and applying these results to ground-penetrating radar depth corrections and dielectric mixing models. Our method includes a correction for subsurface reflector dip angle, which we show can lead to an incorrect determination of wave speeds using common midpoint configurations. By observing the radar properties of the rock glaciers and their supraglacial debris, we find that some of the sites exhibit nearly pure ice cores, and all of the sites indicate volumetric ice fractions >50%. These results have implications for terrestrial glaciology and hydrology because the present ice volume is connected to past ice accumulation and subsurface ice preservation, which may affect the future availability of alpine water resources. An understanding of the processes that govern rock glacier evolution over a wide range of latitudes and elevations will also contribute to the exploration of planetary surfaces such as Mars, which hosts a significant population of debris-covered glaciers. Our subsurface composition and geometry estimates will inform simulations of rock glacier formation and evolution to test hypothesized ice origin mechanisms along with the preservation of climate signals.

Data assimilation in high-order ice flow modeling is a challenging and computationally costly task, yet crucial to find ice thickness and ice flow parameter distributions that are consistent with ice flow mechanics and mass balance while best matching observations. Failing to find these distributions that are required as initial conditions leads to a disequilibrium between mass balance and ice flow, resulting in nonphysical transient effects in the prognostic model. Here we tackle this problem by inverting an emulator of the Stokes ice flow model based on deep learning. By substituting the ice flow equations using a convolutional neural network emulator, we simplify, make more robust and dramatically speed up the solving of the underlying optimization problem thanks to automatic differentiation, stochastic gradient methods and implementation of graphics processing unit (GPU). We demonstrate this process by simultaneously inferring the ice thickness distribution, ice flow parametrization and ice surface of ten of the largest glaciers in Switzerland. As a result, we obtain a high degree of assimilation while guaranteeing an equilibrium between mass-balance and ice flow mechanics. The code runs very efficiently (optimizing one large-size glacier at 100 m takes < 1 min on a laptop) while it is open-source and publicly available.

The future contributions of the Antarctic Ice Sheet to sea level rise will depend on the evolution of its surface mass balance (SMB), which could amplify/dampen mass losses increasingly observed at the ice sheet's edge. In situ constraints of SMB over annual-to-decadal timescales consist mostly of firn/ice cores that have a surface footprint $\sim$cm$^{2}$. SMB constraints also come from climate models, which have a higher temporal resolution but a larger surface footprint of several km$^{2}$. We use ice-penetrating radar data to obtain an intermediate spatial and temporal resolution SMB record over three ice rises along the Princess Ragnhild Coast. The co-located ice cores allow us to obtain absolute radar-derived SMB rates at a multi-annual-to-decadal temporal resolution. By comparing the ice core SMB measurements and the radar-derived SMB records, we determine that pointwise measurements of SMB are representative of a small surface area, $\sim 200-500$ m radius extending from the ice core drill site for the ice rises studied here, and that the pointwise measurements are systematically 7–15 cm w.e. a$^{-1}$ lower than the mean SMB value calculated for the whole ice rises. However, ice core records are representative of an entire ice rise's temporal variability at the temporal resolution examined.

Constructional sedimentary history of coastal barriers can provide insights regarding meteorological and oceanographic processes, and relative sea-level changes. We investigated the Holocene evolution of a segment of the Rio Grande do Sul Coastal Plain in southernmost Brazil. Data were obtained from ground-penetrating radargrams, optically stimulated luminescence dating, altimetry measurements, Google Earth imagery, and aerial photographs. These data allowed a qualitative and quantitative analysis of the beach profiles recorded in the radargrams. From which eolian, backshore/foreshore, breaker bars, and upper and lower shoreface radar facies were identified. The beach-related radar facies are recorded in twenty-eight progradational sand units (sigmoidal bodies). These units record increments of relatively steady deposition bounded by erosional surfaces produced by anomalous, high-magnitude storm events taking place about every 250 years. The upper shoreface strata include two to three breaker bars. Several pieces of evidence (number of breaker bars, upper shoreface gradient and progradation rate) suggest an alternation between dissipative and intermediate stages of barrier morphodynamics and a decrease of wave energy from 2 ka onwards. The barrier prograded during the last 7.1 ka, and initially, barrier progradation occurred because of a normal regression during a relative sea-level rise followed by a stillstand. Later, barrier progradation took place as a forced regression driven by a relative sea-level fall.

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.

Understanding glacier drainage system behaviour and its response to increased meltwater production faces several challenges in the High Arctic because many glaciers are transitioning from polythermal to almost entirely cold thermal structures. We, therefore, used ground-penetrating radar data to investigate the thermal structure and drainage system of Waldemarbreen in Svalbard: a small High Arctic glacier believed to be undergoing thermal change. We found that Waldemarbreen retains up to 80 m of temperate ice in its upper reaches, but this thickness most likely is a relict from the Little Ice Age when greater ice volumes were insulated from winter cooling and caused greater driving stresses. Since then, negative mass balance and firn loss have prevented latent heat release and allowed near-surface ice temperatures to cool in winter, thus reducing the thickness of the temperate ice. Numerous reflectors that can be traced up-glacier are interpreted as englacial channels formed by hydrofracturing in the crevassed upper region of the glacier. The alternative cut and closure mechanism of conduit initiation only forms conduits in parts of the lower ablation area. Consequently, Waldemarbreen provides evidence that hydrofracturing at higher elevations can play a major role in englacial water drainage through cold ice.

The coastal Droning Maud Land in East Antarctica is characterized by small ice shelves with numbers of promontories and locally grounded isles, both called ice rises. These ice rises are typically dome-shaped and surface elevations are hundreds of meters above the surrounding ice shelves, which cause strong orographic effects on surface mass balance (SMB). We conducted shallow ice-penetrating radar sounding to visualize firn stratigraphy in the top 35 m over ~400 km of profiles across the Nivlisen Ice Shelf, and in a grid pattern over two adjacent ice rises (Djupranen and Leningradkollen). We tracked six reflectors (isochrones) and dated them using two ice cores taken at the ice rise summits, from which SMB over six periods in the past three decades was retrieved. The overall SMB pattern across the ice shelf remained similar for all periods; however, the eastwest contrast in SMB varies by a factor of 1.5–2 between the Leningradkollen and Djupranen grounding lines. The SMB patterns over the ice rises are more varied owing to complex interactions between topography, snowfall and wind. We use our results to evaluate the regional climate model RACMO2.3p2 in terms of the spatial SMB distribution and temporal changes over the ice shelf and ice rises at regional scale.

For tribes whose preservation values and mitigation strategies for managing cultural heritage are built on an ethic of avoidance and minimal disturbance, geophysical technologies can be key components of the research design. These technologies, most notably ground-penetrating radar, have been used with great success in identifying and evaluating the depth, extent, and composition of some of those resources for heritage research and management purposes, easing tensions when working with sensitive ancestral places. Additionally, research in archaeological geophysics has shifted from feature finding in order to excavate targets of interest to the recognition that geophysical survey can provide data and interpretations for whole sites and landscapes complementary to or beyond that of excavation, especially regarding the intactness and sensitivity of cultural heritage sites. This use of geophysics as a primary method for research rather than a precursor to archaeological research has empowered tribes with another tool to advocate for low-impact investigation of ancestral sites and landscapes that position tribes as pro-science. Geophysical technologies provide scientifically rigorous yet minimally impactful strategies for investigating heritage while satisfying the requirements of academic and compliance archaeology in ways that can also be culturally appropriate for a much broader spectrum of tribal cultural heritage under consideration.

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.

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.