Ice rises hold valuable records revealing the ice dynamics and climatic history of Antarctic coastal areas from the Last Glacial Maximum to today. This history is often reconstructed from isochrone radar stratigraphy and simulations focusing on Raymond arch evolution beneath the divides. However, this relies on complex ice-flow models where many parameters are unconstrained by observations. Our study explores quad-polarimetric, phase-coherent radar data to enhance understanding near ice divides and domes, using Hammarryggen Ice Rise (HIR) as a case study. Analysing a 5 km profile intersecting the dome, we derive vertical strain rates and ice-fabric properties. These align with ice core data near the summit, increasing confidence in tracing signatures from the dome to the flanks. The Raymond effect is evident, correlating with surface strain rates and radar stratigraphy. Stability is inferred over millennia for the saddle connecting HIR to the mainland, but dome ice-fabric appears relatively young compared to 2D model predictions. In a broader context, quad-polarimetric measurements provide valuable insights into ice-flow models, particularly for anisotropic rheology. Including quad-polarimetric data advances our ability to reconstruct past ice flow dynamics and climatic history in ice rises.

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.

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

During the past 20 years, multi-channel radar emerged as a key tool for deciphering an ice sheet's internal architecture. To assign ages to radar reflections and connect them over large areas in the ice sheet, the layer genesis has to be understood on a microphysical scale. Synthetic radar trace modelling based on the dielectric profile of ice cores allows for the assignation of observed physical properties’ variations on the decimetre scale to radar reflectors extending from the coring site to a regional or even whole-ice-sheet scale. In this paper we rely on the available dielectric profiling data of the northern Greenland deep ice cores: NGRIP, NEEM and EGRIP. The three records are well suited for assigning an age model to the stratigraphic radar-mapped layers, and linking up the reflector properties to observations in the cores. Our modelling results show that the internal reflections are mainly due to conductivity changes. Furthermore, we deduce fabric characteristics at the EGRIP drill site from two-way-travel-time differences of along and across-flow polarized radarwave reflections of selected horizons (below 980 m). These indicate in deeper parts of the ice column an across-flow concentrated c-axis fabric.

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.

The internal stratigraphy of snow and ice as imaged by ground-penetrating radar may serve as a source of information on past accumulation. This study presents results from two ground-based radar surveys conducted in Greenland in 2007 and 2015, respectively. The first survey was conducted during the traverse from the ice-core station NGRIP (North Greenland Ice Core Project) to the ice-core station NEEM (North Greenland Eemian Ice Drilling). The second survey was carried out during the traverse from NEEM to the ice-core station EGRIP (East Greenland Ice Core Project) and then onwards to Summit Station. The total length of the radar profiles is 1427 km. From the radar data, we retrieve the large-scale spatial variation of the accumulation rates in the interior of the ice sheet. The accumulation rates range from 0.11 to 0.26 m a−1 ice equivalent with the lowest values found in the northeastern sector towards EGRIP. We find no evidence of temporal or spatial changes in accumulation rates when comparing the 150-year average accumulation rates with the 321-year average accumulation rates. Comparisons with regional climate models reveal that the models underestimate accumulation rates by up to 35% in northeastern Greenland. Our results serve as a robust baseline to detect present changes in either surface accumulation rates or patterns.

The Northeast Greenland Ice Stream (NEGIS) is an important dynamic component for the total mass balance of the Greenland ice sheet, as it reaches up to the central divide and drains 12% of the ice sheet. The geometric boundary conditions and in particular the nature of the subglacial bed of the NEGIS are essential to understand its ice flow dynamics. We present a record of more than 8000 km of radar survey lines of multi-channel, ultra-wideband radio echo sounding data covering an area of 24 000 km2, centered on the drill site for the East Greenland Ice-core Project (EGRIP), in the upper part of the NEGIS catchment. Our data yield a new detailed model of ice-thickness distribution and basal topography in the region. The enhanced resolution of our bed topography model shows features which we interpret to be caused by erosional activity, potentially over several glacial–interglacial cycles. Off-nadir reflections from the ice–bed interface in the center of the ice stream indicate a streamlined bed with elongated subglacial landforms. Our new bed topography model will help to improve the basal boundary conditions of NEGIS prescribed for ice flow models and thus foster an improved understanding of the ice-dynamic setting.

Radar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of the Annals of Glaciology on ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciological Annals of Glaciology issue of 2014, the context of their history and future prospects.

The high-Alpine ice-core drilling site Colle Gnifetti (CG), Monte Rosa, Swiss/Italian Alps, provides climate records over the last millennium and beyond. However, the full exploitation of the oldest part of the existing ice cores requires complementary knowledge of the intricate glacio-meteorological settings, including glacier dynamics. Here, we present new ice-flow modeling studies of CG, focused on characterizing the flow at two neighboring drill sites in the eastern part of the glacier. The3-D full Stokes ice-flow model is thermo-mechanically coupled and includes firn rheology, firn densification and enthalpy transport, and is implemented using the finite element software Elmer/Ice. Measurements of surface velocities, accumulation, borehole inclination, density and englacial temperatures are used to validate the model output. We calculate backward trajectories and map the catchment areas. This constrains, for the first time at this site, the so-called upstream effects for the stable water isotope time series of the two ice cores drilled in 2005 and 2013. The model also provides a 3-D age field of the glacier and independent ice-core chronologies for five ice-core sites. Model results are a valuable addition to the existing glaciological and ice-core datasets. This especially concerns the quantitative estimate of upstream conditions affecting the interpretation of the deep ice-core layers.

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.

The detection and monitoring of meltwater within firn presents a significant monitoring challenge. We explore the potential of small wireless sensors (ETracer+, ET+) to measure temperature, pressure, electrical conductivity and thus the presence or absence of meltwater within firn, through tests in the dry snow zone at the East Greenland Ice Core Project site. The tested sensor platforms are small, robust and low cost, and communicate data via a VHF radio link to surface receivers. The sensors were deployed in low-temperature firn at the centre and shear margins of an ice stream for 4 weeks, and a ‘bucket experiment’ was used to test the detection of water within otherwise dry firn. The tests showed the ET+ could log subsurface temperatures and transmit the recorded data through up to 150 m dry firn. Two VHF receivers were tested: an autonomous phase-sensitive radio-echo sounder (ApRES) and a WinRadio. The ApRES can combine high-resolution imaging of the firn layers (by radio-echo sounding) with in situ measurements from the sensors, to build up a high spatial and temporal resolution picture of the subsurface. These results indicate that wireless sensors have great potential for long-term monitoring of firn processes.

The accumulation rate on Potsdam Glacier, East Antarctica, and its spatial and temporal variations are examined using ground-penetrating radar, snow samples and firn-core studies. Physical properties in snow samples and along firn cores provide distributions of density with depth, showing only small spatial variation. Counting of peaks in δ18O along the firn cores yields an age–depth distribution that is transferred to the stratigraphy of isochronal internal layers observed with radar. From two radar horizons we determine the spatial accumulation pattern, averaged over the periods 1970–80 and 1980–2004. The shape of internal layers indicates an ablation area at the eastern margin of the investigation area. Accumulation rates show a very high spatial variability, with a mean value of 141 kgm–2 a–1 for the period 1970–2004 and a standard deviation of almost 50%. Mean temporal variation of only a few per cent throughout the investigated area for the observed time interval is much less than the spatial variations. The mean accumulation values are somewhat less than values reported before from this region. Accumulation pattern and surface topography are linked in a way indicating that wind-borne redistribution of snow significantly contributes to the observed spatial variations of accumulation rates. The accumulation data and their variability complement and validate present and future satellite studies of Antarctica’s mass balance.

The distribution of surface mass balance on Amundsenisen, Dronning Maud Land, Antarctica, is investigated along a continous profile line. Ice-penetrating radar is used to map variations in ice-layer thickness within the upper 100 m of the ice sheet. The route passes several firn- and ice-core drilling sites over a distance of 320 km. Dielectric-profiling data of ice cores are used to calculate the depths of selected reflection horizons and the cumulative mass of the ice column. The local surface mass balance is determined as a temporal average, covering a time-span of almost two centuries. The findings indicate a complex accumulation pattern superimposed on a generally low surface mass balance, which is related to small-scale surface undulations. The results of the radar soundings are in general in good agreement with surface mass-balance data derived from firn-core studies. Discrepancies between these two datasets can be explained by spatial mismatch or by minor quality of either ice-core profiles or radar data. For regional comparison of radar-based accumulation data we use an accumulation distribution interpolated from point measurements. The surface mass balance varies up to 50% over short distances, with correlation lengths of <10 km. We conclude that the current utilization schemes of point sampling are only capable of reproducing local values and regional trends but provide no information on the small-scale variability of surface mass balance.

Isochronic continuous horizons between 20 and 90 m depth in a ground-penetrating radar (GPR) profile, recorded in Dronning Maud Land, Antarctica, are identified by comparison of synthetic and measured single radar traces. The measured radar-gram is derived from a stacked GPR profile; the synthetic radargram is computed by convolution of the complex reflection coefficient profile, based on dielectric profiling (DEP) data of a 150 m ice core, with a depth-invariant wavelet. It reproduces prominent reflections of the measured radargram to a considerable degree. Analyzing matching peaks in both radargrams enables us to identify isochronic reflections and transfer individual volcanic-event datings to the GPR profile. Reflections are primarily caused by changes in permittivity; changes in conductivity are of minor importance. However, several peaks in permittivity andconductivity show a good correlation and indicate that some reflections are related to acidic layers. The results demonstrate the possibility of reproducing radargrams from ice-core property profiles, a necessary step for the interpretation of remotely sensed radar data and the general significance of connecting ice-core and radar data for correct interpretations. Problems related to forward modeling, data gaps, origin of permittivity peaks, and GPR profiles used for comparison, are discussed.

During the summer of 2003, a ground-penetrating radar survey around the North Greenland Icecore Project (NorthGRIP) deep ice-core drilling site (75˚06’N, 42˚20’W; 2957ma.s.l.) was carried out using a shielded 250 MHz radar system. The drill site is located on an ice divide, roughly 300 km north-northwest of the summit of the Greenland ice sheet. More than 430 km of profiles were measured, covering a 10 km by 10 km area, with a grid centered on the drilling location, and eight profiles extending beyond this grid. Seven internal horizons within the upper 120 m of the ice sheet were continuously tracked, containing the last 400 years of accumulation history. Based on the age-depth and density-depth distribution of the deep core, the internal layers have been dated and the regional and temporal distribution of accumulation rate in the vicinity of NorthGRIP has been derived. The distribution of accumulation shows a relatively smoothly increasing trend from east to west from 145 kgm–2a–1 to 200 kg m–2 a -1 over a distance of 50 km across the ice divide. The general trend is overlain by small-scale variations on the order of 2.5 kgm–2a-1 km- 1 , i.e. around 1.5% of the accumulation mean. The temporal variations of the seven periods defined by the seven tracked isochrones are on the order of ± 4% of the mean of the last 400 years, i.e. at NorthGRIP ± 7 kg m–2 a-1. If the regional accumulation pattern has been stable for the last several thousand years during the Holocene, and ice flow has been comparable to today, advective effects along the particle trajectory upstream of NorthGRIP do not have a significant effect on the interpretation of climatically induced changes in accumulation rates derived from the deep ice core over the last 10 kyr.

Based on ground-penetrating radar profiles, we analyze area-wide spatial and temporal characteristics of accumulation rate in the vicinity of the EPICA (European Project for Ice Coring in Antarctica) deep-drilling site in Dronning Maud Land, Antarctica (EDML). An area of 1600km2 is covered by 500km of radar profiles, organized in a star-like pattern with eight 20–25km legs and a 10km grid with 1–3km spacing, each pattern being centred on the EDML drilling location. Distributions of density, cumulative mass, age and the electromagnetic wave speed with depth are available from physical ice-core records. Nine internal reflection horizons are continuously tracked within the upper 110m of ice over the whole area, yielding a spatial picture of accumulation rate history over >1000 years. The mean accumulation rate over the last 153 years varies between 50 and 75 kgm–2 a–1 over 50 km perpendicular to the ice divide; the spatial average is ~61 kgm–2 a–1. This general pattern is overlain by small-scale variations of accumulation rate on the order of 10% of the mean. Maximum local gradients in accumulation rate are ~2–3 kgm–2 a–1 km–1, about five times the regional accumulation rate gradient. Comparison of topography and accumulation rate along a 20 km profile in the direction of the mean winds indicates that variations in accumulation rate over short distances are linked to surface undulations. Our results show that advected spatial variations of accumulation rate are on the same order and even exceed temporal changes over the investigated periods. Ice flow and upstream effects therefore might influence accumulation rates reconstructed from the EDML ice core.

We investigate snowpack properties at a site in west-central Greenland with ground-penetrating radar (GPR), supplemented by stratigraphic records from snow pits and shallow firn cores. GPR data were collected at a validation test site for CryoSat (T05 on the Expéditions Glaciologiques Internationales au Groenland (EGIG) line) over a 100 m × 100 m grid and along 1 km sections at frequencies of 500 and 800 MHz. Several internal reflection horizons (IRHs) down to a depth of 10 m were tracked. IRHs are usually related to ice-layer clusters in vertically bounded sequences that obtain their initial characteristics near the surface during the melt season. Warm conditions in the following melt season can change these characteristics by percolating meltwater. In cold conditions, smaller melt volumes at the surface can lead to faint IRHs. The absence of simple mechanisms for internal layer origin emphasizes the need for independent dating to reliably interpret remotely sensed radar data. Our GPR-derived depth of the 2003 summer surface of 1.48 m (measured in 2004) is confirmed by snow-pit observations. The distribution of IRH depths on a 1 km scale reveals a gradient of increasing accumulation to the northeast of about 5 cm w.e. km−1. We find that point measurements of accumulation in this area are representative only over several hundred metres, with uncertainties of about 15% of the spatial mean.