Loading [MathJax]/jax/output/HTML-CSS/fonts/TeX/fontdata.js
Articles | Volume 19, issue 3
https://doi.org/10.5194/tc-19-1153-2025
https://doi.org/10.5194/tc-19-1153-2025
Research article
 | 
13 Mar 2025
Research article |  | 13 Mar 2025

Age–depth distribution in western Dronning Maud Land, East Antarctica, and Antarctic-wide comparisons of internal reflection horizons

Steven Franke, Daniel Steinhage, Veit Helm, Alexandra M. Zuhr, Julien A. Bodart, Olaf Eisen, and Paul Bons
Abstract

Radio-echo sounding provides the opportunity to study the internal architecture of ice sheets through imaging stratified englacial reflections, known as internal reflection horizons (IRHs). They represent consistent time horizons formed at the former ice-sheet surface and buried over time, thus reflecting the ice sheet's age–depth architecture. Their analysis allows crucial insights into past and present glaciological conditions, e.g. bed topography, surface and basal mass balance, and physical properties and ice dynamics. This study presents a comprehensive data set of IRHs and insight into the age–depth distribution in western Dronning Maud Land (DML), East Antarctica, spanning the Holocene to the Last Glacial Period (4.8–91.0 ka). Using data from various radar systems deployed by the Alfred Wegener Institute between 1996 and 2023, we traced and dated nine IRHs over an area of 450 000 km2. A precise age could be assigned to the IRHs by two-way travel time to depth conversion and employing radar forward modelling based on conductivity peaks of the EPICA DML ice core. Six IRHs correlate with the timing of past volcanic eruptions, and our findings suggest that most IRHs correspond to IRHs of similar age in other regions of East and West Antarctica, thus likely originating from the same physical reflectors at depth, although some could not be physically connected. This work enhances understanding of the englacial architecture and relationships with snow accumulation and ice-dynamic processes of this sector of the Antarctic ice sheet and provides boundary conditions for numerical ice flow models and paleoclimatic studies.

Share
1 Introduction

Studying the dynamics of the Antarctic ice sheet (AIS) through a combination of geophysical observations and ice-sheet modelling is crucial for understanding its response to climate change and better projecting future sea-level rise (IPCC2023). Deducing its past evolution helps us to understand the processes driving ice flow, accumulation, and basal melting. As one unique observable property, the ice sheet's internal stratigraphy represents a valuable record of past kinematics and dynamics (Siegert et al.2004; Bons et al.2016; Leysinger Vieli et al.2018; Jansen et al.2024), offering insights into atmosphere–ice–ocean interactions (Drews et al.2020; Višnjević et al.2022) and how the ice-sheet system responded to natural climatic variations over a broad range of timescales ranging from hundreds to hundreds of thousands of years (e.g. Leysinger Vieli et al.2004; Sutter et al.2021; Bodart et al.2023).

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f01

Figure 1Radar profiles used for IRH tracing. (a) Radar profiles are sorted by system: EMR short pulse (white lines), UWB narrowband (orange lines), and UWB wideband (blue lines). (b) Radar data sorted by acquisition season. Although some profiles extend down to the ice shelf, the IRH tracing region is restricted upstream of the Maud Belt. The background map shows ice surface flow velocity (Mouginot et al.2019a) superimposed on the hillshaded surface elevation (Howat et al.2019). The IMBIE (Ice sheet Mass Balance Inter-comparison Exercise) drainage basins (Rignot et al.2019) are marked by thin black lines, and deep ice core sites are marked with a white circle.

Radar data provide the opportunity to study the internal architecture of ice sheets (i.e. their stratigraphy) via continuous stratified englacial reflections (Bingham et al.2024). These internal reflection horizons (IRHs) occur when radar waves are reflected at boundary layers where the dielectric properties of the ice change (Glen and Paren1975; Robin et al.1977). Dielectric contrasts in ice sheets are primarily caused by variations in density (Robin et al.1969), electrical conductivity (Paren and Robin1975), and changes in the preferred orientation of ice crystals (Harrison1973; Fujita et al.1999). IRHs that are located sufficiently deep below the ice surface, typically within a few hundred metres in Antarctica, are primarily linked to conductivity and density contrasts and can sometimes be detected thousands of kilometres across the ice sheet (Millar1981; MacGregor et al.2015; Winter et al.2019). Strong conductivity contrasts are often linked to deposits of acidic materials from major volcanic eruptions (Hammer1980; Millar1981). Moreover, IRHs are often considered isochronous because they form simultaneously at the ice-sheet surface and are then buried under subsequent snow accumulation, carrying with them changes in dielectric properties which radar systems are sensitive to (Bingham et al.2024).

Information about the internal stratigraphy and age–depth architecture of the Antarctic ice sheet is crucial as it serves as a climatic and ice-dynamic record. For instance, it provides insights into past accumulation rates and basal melting, informing us about past surface and basal mass balance (Leysinger Vieli et al.2004; Eisen et al.2005; Leysinger Vieli et al.2011; Cavitte et al.2018; Bodart et al.2023; Koch et al.2023). Since many IRHs can be dated with ice cores, they can also serve as calibration points for ice-sheet models (Sutter et al.2021; Višnjević et al.2022) and to model basal conditions (Leysinger Vieli et al.2018; Fudge et al.2023; Chung et al.2023b; Wang et al.2023a). Additionally, the geometry, depth, and continuity of IRHs contain information on the cumulative deformation due to ice flow, making them suitable passive markers for ice-dynamic activity and its changes over time, e.g. deciphering folding processes and ice stream activity (Siegert et al.2004; Bons et al.2016; Franke et al.2022a; Jansen et al.2024), as well as ice-dynamic processes causing IRH discontinuity (Panton and Karlsson2015; Sanderson et al.2023b).

The generation and compilation of IRHs using radar data in Antarctica and Greenland remain largely a manual task (Bingham et al.2024). While automated methods can successfully trace IRHs to a certain extent (see e.g. review by Moqadam and Eisen2024), these algorithms often fail to consistently follow the same reflection across multiple intersecting radar profiles and different radar systems, where reflections have a different radar signature due to their different wavelengths. Nevertheless, IRHs have been traced and dated across large parts of East Antarctica (Leysinger Vieli et al.2011; Winter et al.2019; Cavitte et al.2016, 2021; Beem et al.2021b; Wang et al.2023a; Sanderson et al.2024), West Antarctica (Siegert et al.2005; Jacobel and Welch2005; Ashmore et al.2020a; Bodart et al.2021b; Muldoon et al.2018), and Greenland (Karlsson et al.2013; MacGregor et al.2015; Franke et al.2023a; MacGregor et al.2025). For Antarctica, in particular, these studies have been motivated by the AntArchitecture initiative, a Scientific Committee on Antarctic Research Action Group which aims to build a 3D age–depth model of Antarctica from IRHs (Bingham et al.2024).

Table 1Alfred Wegener Institute (AWI) radar system specifications. The range resolution applies for a dielectric permittivity of εr=3.15 in ice.

a Hamburg University of Technology. b Center for Remote Sensing and Integrated Systems, formerly known as Center for Remote Sensing of Ice Sheets.

Download Print Version | Download XLSX

In this study, we present a comprehensive and detailed insight into the dated radiostratigraphy of western Dronning Maud Land (DML; East Antarctica), spanning from the Holocene to the Last Glacial Period, using radar data collected over the past 3 decades by various recording systems of the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI; Fig. 1). Through radar forward modelling based on conductivity data from the EPICA DML (EDML) ice core, we can accurately date IRHs and, in most cases, link them to deposits from past volcanic eruptions. We discuss the age–depth distribution and links to glaciological parameters, such as ice flow velocity, accumulation rates, and bed topography. Moreover, we infer the extent to which our dated IRHs can be correlated with those from other Antarctic-wide studies and demonstrate the potential for extrapolating some of these IRHs to much larger areas where an AWI radar system with coarser vertical resolution was used.

2 Data and methods

2.1 Radar data

In this study, we utilize three types of radar products that provide sufficiently high range resolution (0.35–5 m) and penetration depth (Table 1). These radar products were obtained during campaigns conducted by AWI between 1996 and 2023 with two different radar systems (Fig. 1). The radar data were collected using AWI's polar aircraft Polar 2, Polar 5, and Polar 6 (Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung2016a), operating between the stations Neumayer and Kohnen (Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung2016b), Novolazarevskaya, Troll, and Dome Fuji.

2.1.1 AWI EMR short-pulse data

The largest pool of radar data is provided by AWI's EMR (Electromagnetic Reflection System), which has been operational in Antarctica and Greenland since 1994 (Nixdorf et al.1999). For IRH tracing, we utilize EMR radar data from nine Antarctic seasons spanning from 1996/97 to 2016/17 (Fig. 1 and Table 1). The EMR system consists of two dipole antennas mounted underneath the wings of AWI's polar aircraft. The transmission signal is a burst with pulse lengths of 60 or 600 ns at a frequency of 150 MHz, with a pulse repetition frequency (PRF) of 20 kHz. For this study, we exclusively utilize EMR 60 ns pulse data (short-pulse data) due to their higher range resolution to detect IRHs compared to the 600 ns pulse data (see Wang et al.2023a). The short-pulse data have a range resolution of approximately 5 m and a vertical sample interval of 13.33 ns. EMR radar data processing comprises filtering and along-track stacking with a factor of 7 or 10, which results in a mean trace spacing of ∼35 and ∼50 m, respectively (Steinhage et al.2001; Steinhage2001; Wang et al.2023a).

EMR survey lines cover the entire western DML (Fig. 1). Surveys between 1997 and 1999 focused on the pre-site survey of the EDML ice core site (Steinhage et al.2001; Steinhage2001), while surveys between 2001 and 2011 were used to investigate the ice internal structure (Eisen et al.2006, 2007; Drews et al.2009; Steinhage et al.2013) as well as ice thickness for solid earth geophysical studies in DML (Riedel et al.2012; Eisermann et al.2020). The survey in 2016/17 extensively surveyed ice thickness at the Dome Fuji ice core site (Karlsson et al.2018; Wang et al.2023a).

2.1.2 AWI ultra-wideband data

In addition, we use two data sets acquired with AWI's ultra-wideband (AWI UWB) radar system. The system is an improved version of the Multichannel Coherent Radar Depth Sounder (MCoRDS, version 5), which was developed at the Center for Remote Sensing and Integrated Systems (CReSIS) at the University of Kansas (Rodriguez-Morales et al.2013; Hale et al.2016) and has been in operation since 2016 (Kjær et al.2018). The radar configuration for both campaigns in 2018/19 and 2023/24 consisted of an eight-element antenna array mounted under AWI's Polar 5 or Polar 6 Basler BT-67 aircraft's fuselage, which serves as a transmit and receive antenna array. The transmission signal is composed of several staged modulated chirp signals, which provide high resolution at different depths (Franke et al.2022b). Standard processing techniques were performed with the OPR Toolbox (Open Polar Radar Toolbox, formerly termed CReSIS Toolbox; Open Polar Radar2023). The main steps include motion compensation, pulse compression, synthetic aperture radar (SAR) focusing in the fk (frequency–wavenumber) domain, and array processing. Aeroplane location and orientation are from Global Positioning System (GPS) precise point positioning (PPP) post-processed with a final estimated accuracy (commercial software package Waypoint 8.4) of better than 3 cm for latitude and longitude and better than 10 cm for altitude.

The two AWI UWB campaigns were flown in two different acquisition settings in the 2018/19 and 2023/24 seasons, respectively. In the 2018/19 season, data were acquired with a frequency range between 180–210 MHz (narrowband), which corresponds to a range resolution of 4.3 m. The acquisition settings were optimized to sound deep englacial reflections as well as the bed topography at the onset of Jutulstraumen Glacier (Franke et al.2021, 2024b). AWI UWB data in the 2023/24 season were acquired with a frequency range between 150–520 MHz (wideband), with a range resolution of  35 cm. The acquisition settings were chosen to map near-surface IRHs at high spatial resolution (Koch et al.2023) in the upstream part of the Jutulstraumen catchment as well as the central plateau south and west of the EDML ice core.

2.2 Internal reflection horizon tracing

For tracing IRHs, we utilize the commercial software Paradigm 19 from AspenTech Subsurface Science & Engineering. To enhance the visibility of isochrones at greater depths, we applied an automatic gain control (AGC) algorithm to the radargrams. IRHs were traced with the following rationale: (1) existing IRHs, such as the 38 and 74 ka isochrones previously traced by Winter et al. (2019) on segments of the same data we used here; (2) the shallowest and deepest clearly identifiable reflection in the EMR short-pulse data was traced; and (3) additional isochrones were traced at regular intervals between those defined in (1) and (2), which were clearly discernible in all radar products.

All isochrones were manually traced with the assistance of a semi-automatic picker, which searches for the local maximum in return power within a window width of 5 ns along two picks. Based on the intersection points of the radar profiles, IRHs were transferred between the different radar products. In cases where the assignment at the intersection point was not clear or the radar reflection was no longer traceable within a profile, no picks were created. We clarify that we use the existing IRHs from Winter et al. (2019) in DML and integrate these data for our representation, dating, and analyses.

2.3 Ice base determination

The determination of the ice base reflection for this study is composed of existing and newly traced data. We used existing ice base pick data from the AWI UWB 2018/19 campaign from Franke et al. (2021). Ice base determination from the AWI UWB 2023/2024 season was generated using the OPR Toolbox (Open Polar Radar2023) and Paradigm. For the EMR data, we were able to partially rely on existing data (Steinhage et al.2001; Riedel et al.2012), which were either directly traced in the 60 ns short-pulse data or, if flown in toggle mode, traced in the 600 ns long-pulse data and projected to the 60 ns data. For EMR short-pulse profiles used in this study that did not have existing ice base picks, the ice base reflection was traced at locations where IRHs were present and the base reflection was clearly visible.

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f02

Figure 2Comparison between simulated radar data for 150 and 335 MHz (150 MHz for the EMR profile and 335 MHz for the AWI UWB wideband profile) based on measured conductivity at the EDML ice core and radar data collected near EDML. (a) Input parameters for dielectric conductivity and relative permittivity. (b) Comparison between the gprMax simulated radar trace at 150 MHz with overlaid conductivity peaks (yellow) and the EMR profile 20023154 (transmit frequency is 150 MHz). (c) Comparison between the gprMax simulated radar trace at 335 MHz with overlaid conductivity peaks and the UWB wideband profile 20231211_01_024 (centre frequency is 335 MHz). The correspondence between simulated reflections, conductivity peaks, and the reflections of our IRHs is marked with coloured arrows. Note that the simulated radargram is the 2-D representation of a single trace (without AGC), and the measured radargram has an along-track range of 20 km (with AGC) to enable a better comparison.

Download

2.4 IRH depth and normalized depth

We derived the depth of all IRHs by calculating the two-way travel time (TWT), t, between the ice surface and the IRH. Then, we converted the TWT to depth (e.g. Steinhage et al.2001; MacGregor et al.2015; Cavitte et al.2016; Winter et al.2017; Ashmore et al.2020a; Bodart et al.2021b) using

Zirh=tcair2εr+Zfirn,

where cair is the electromagnetic wave speed in air (cair=2.9971×108 m s−1). For the calculation of IRH depths, we use a constant dielectric permittivity of εr=3.15. This value deviates slightly from the empirically determined value at EDML (Eisen et al.2006), but it is a realistic estimate for the spatial extrapolation within our study area and has been widely used in other studies. Additionally, we apply a firn correction of 13 m (Steinhage2001), which is the best-known value in this region. However, even though a 13 m firn correction was applied uniformly, snowfall and accumulation rates vary across western Dronning Maud Land (e.g. Rotschky et al.2007), potentially affecting depth calculations of IRHs slightly, i.e. on the order of a few metres.

Additionally, we calculated the normalized depth of the IRHs within the ice column, which provides a better comparison to determine whether IRHs lie deeper or shallower within the ice column in specific areas. For this, we divided the IRH depth below the ice surface by the respective ice thickness. Here, we only use ice thickness values at locations where ice base picks are available, as differences in ice thickness compared to gridded products like BedMachine Antarctica (Morlighem et al.2020) can result in significant discrepancies, which can strongly influence the relative position of IRHs within the ice column. These discrepancies are due to interpolation methods and the fact that not all ice thickness data used here are yet included in BedMachine Antarctica.

2.5 Internal reflection horizon dating

To determine the age of the IRHs, we used two different methods: (1) age dating by converting the IRH radar wave two-way travel time (TWT) from the ice surface to depth and correlating it with the age–depth scale of the EDML ice core (e.g. Cavitte et al.2016; Bodart et al.2021b) and (2) matching reflections of our IRHs with simulated radar reflections based on conductivity peaks derived from measured DEP (dielectric profiling) data (e.g. Eisen et al.2003, 2006; Winter et al.2017). Method (1) has the advantage of allowing the determination of an age for each IRH, albeit with the acceptance of a larger error in depth and age assignment. Method (2) offers greater accuracy and linkage to volcanic events because the conductivity peaks represent the physical origin of reflections (Eisen et al.2003); however, this approach is not applicable for all IRHs due to some reflections not rising above background noise in the radar data.

Table 2Overview of IRH key characteristics, such as the mean IRH TWT at EDML, mean IRH depth at EDML derived from TWT-to-depth conversion, IRH depths at EDML derived from matching conductivity peaks, IRH ages based on the TWT-to-depth converted depths, IRH ages based on the depth from matched conductivity peaks, cumulative length of IRHs (line kilometres), and link to possible volcanic eruption ages. The IRH name contains the region (DML) and age in ka (e.g. 4_8_ka for 4.8 ka). All TWT-to-depth converted depth values are relative to the ice surface reflection at first break. The conductivity peak depths are from the DEP data set and refer to the depth of the maximum value of the respective conductivity peak. The IRH ages highlighted in bold italics in the TWT-to-depth-based age and conductivity-peak-based age columns are the ones we have selected as the final IRH ages.

a Large Northern Hemisphere higher-latitude eruption (rank 18 in Lin et al.2022).
b Large lower-latitude or Southern Hemisphere eruption (rank 32 in Lin et al.2022).
c Taupo Oruanui eruption (Dunbar et al.2017).
d Large lower-latitude eruption (rank 3 in Lin et al.2022).
e Large bipolar eruption (rank 5 in Lin et al.2022).
f Toba eruption dated to 75.0 ± 0.9 ka (see references in Svensson et al.2013).
g Toba eruption dated to 73.88 ± 0.32 ka (see references in Svensson et al.2013).

Download Print Version | Download XLSX

2.5.1 IRH dating by TWT-to-depth conversion

For the IRH dating using the AICC2023 age–depth scale of the EDML ice core (Bouchet et al.2023), we used three radar profiles from the seasons 1997/98, 2001/02, and 2023/24, which are located close to the drilling site. Based on the radar traces closest to the ice core for each radar profile (25 m for profile 19983101, 65 m for profile 20023154, and 280 m for profile 20231211_01_024), we determined the travel time between the ice surface reflection and the IRHs. For the depth conversion at EDML, we used the constant radar wave speed in ice suggested by Eisen et al. (2006) of  1.69 × 108 m s−1, corresponding to a dielectric permittivity of εr=3.145. Additionally, we added a value of 13 m to the IRH depths as a firn correction, calculated from measurements of the complex permittivity of shallow firn cores in the region (Steinhage2001; Steinhage et al.2001). To account for the variation in acquisition years and radar products, we averaged the depth values for all three radar profiles. We disregarded the snow accumulation between the recording times of the radar profiles, as the average accumulation rate at EDML is approximately 5 cm ice equivalent per year, which amounts to about 1.25 m over 25 years and is significantly below the range resolution of 5 m of the EMR system.

For the depth determination error, we included the following parameters: (1) the standard deviation of the depth differences of the respective IRHs from the three radar profiles; (2) the smallest range resolution of 5 m from the EMR system; (3) the span of the tracing window of 5 ns, which corresponds to approximately 0.5 m; and (4) an error in the dielectric permittivity of 1 % (Bohleber et al.2012). In addition to the age error resulting from the depth error range, the uncertainty of the age in the AICC2023 chronology is also included.

2.5.2 IRH dating by DEP-based radar forward modelling

For age assignment based on simulated radar reflections derived from the measured conductivity of the EDML ice core (Eisen et al.2006; Winter et al.2017; Mojtabavi et al.2022), we used the open-source software gprMax (Giannopoulos2005; Warren et al.2016) to simulate electromagnetic waves in three-dimensional space. This software solves Maxwell's equations using the finite-difference time-domain (FDTD) method (Taflove and Hagness2005) and Yee cells (Yee1966), which has proven effective for 1-D and 2-D simulations in ice sheets and glaciers (Franke et al.2023b; Santin et al.2023). gprMax uses the same algorithm as already successfully employed for ice-core-based forward modelling by Eisen et al. (2006) but includes an efficient solution for model boundary conditions using perfectly matched layers (PMLs; Giannopoulos2012) while still offering reasonable computation times for 2-D models. All simulations presented here were performed with gprMax version 3.1.5 (Big Smoke).

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f03

Figure 3Representation of a subset of radargrams (AGC version) of the three different radar products used in this study and traced IRHs in the western DML. (a) Visualization of selected radar data in a 3D canvas with the bed topography (Morlighem et al.2020) and traced IRHs. Radar profiles with a black line at the top of the fence radargrams are EMR short-pulse profiles, those with a red line are narrow-bandwidth profiles, and those with a yellow line are wide-bandwidth profiles. (b) Representation of an EMR short-pulse radargram (profile ID 20023154), (c) a narrow-bandwidth radargram (profile ID 20181227_03_001), and (d) a wide-bandwidth radargram (profile ID 20231211_01_024) with IRHs. The radargrams in panels (b)(d) are located close to the EDML ice core and are aligned with the surface reflection.

Download

For our simulations, we defined a 2-D model domain extending 2.4 m in the x direction (width) and 2774 m in the z direction (depth) to represent the full depth of the ice core. To reduce the model size, the y-direction extension was limited to a single cell size, rendering the model in 2-D to effectively operate in the transverse-magnetic mode. The cell size in all dimensions is 0.02 m, which enables high resolution and numerical stability and avoids numerically induced dispersion by ensuring the cell size is at least 10 times smaller than the smallest wavelength in our model. Consequently, the time step Δt is 0.047 ns for our cell size. The transmit and receive antennas are placed 1 m below the upper limit of the model domain, in the centre of the x domain (zero offset). The simulation window for all simulations is 30 µs, covering the entire ice thickness at EDML. The antenna represents a Hertzian dipole transmitting a Ricker wavelet polarized in the y direction. For all simulations, we defined 30 PML cells at the outer boundary of the x and z directions.

We conducted simulations using two fixed transmission frequencies: 150 and 335 MHz. These frequencies represent the centre frequencies of the two radar profiles located close to the ice core (150 MHz for the EMR profile and 335 MHz for the UWB wideband profile). As input data, we used the conductivity measured by dielectric profiling (DEP; Moore1993; Wilhelms et al.1998) along the EDML ice core. We used a constant dielectric permittivity value of εr=3.145 for the depth range of solid ice. This corresponds to an EM wave propagation velocity of  1.69 × 108 m s−1, which is the average value determined by Eisen et al. (2006) at EDML. For the upper  180 m, we fitted a function to the permittivity measured with the DEP until it reached a constant value. The relative magnetic permeability was set to a constant value of μr=1, and the magnetic loss factor was set to σ*=0.

We applied the following signal-processing steps to the simulation output. We derived the reflected energy P of the synthetic radar data following the methods described by Winter et al. (2017) and Franke et al. (2023b). This involved calculating the envelope of the dielectric field strength (y component of the electric field, Ey) using the Hilbert magnitude transform. Additionally, the radar trace was smoothed using a 1-D Gaussian filter. This approach enhances the representation of the impulse response of a transmitted radar signal and reflections from a smooth reflector. The final signal strength was obtained by converting the electric field envelope (Ey) to decibel (dB) (20log10(|Ey|)). We did not apply corrections for geometric spreading or englacial attenuation.

For the age assignment, we examined which of the simulated reflections originated from conductivity peaks and matched the IRHs (Eisen et al.2006). We then determined the depth of the corresponding conductivity peaks (if present and unambiguous to determine) and used this depth to assign the age via the AICC2023 age–depth chronology (Bouchet et al.2023). This method is less dependent on the choice of radar wave propagation velocity in ice, firn correction factor, and ice surface determination, relying solely on the assignment of conductivity peaks, and provides a clearer age determination with a smaller associated error. Hence, for the error estimation, we only consider the dating error provided in the AICC2023 age–depth chronology. In addition, it can be used to better estimate the average relative permittivity of the ice at the ice core site.

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f04

Figure 4Location, vertical overlap, and spatial distribution of IRHs. (a) Spatial distribution of IRHs from this study (dark blue) in DML and of IRHs from other studies in this region (Winter et al.2019; Wang et al.2023a, orange and purple, respectively). (b) Colour-coded overlap of IRHs. (c–e) Magnified views showing the spatial coverage of IRHs in key regions: (c) Jutulstraumen Glacier onset, (d) EDML ice core site, and (e) ice divide southwest of EDML (locations indicated by white boxes in panel b).

2.6 IRH depth and age uncertainty quantification

We follow Cavitte et al. (2016, 2021) and Winter et al. (2019) to quantify the depth and age uncertainty for all IRH depths and consider the following sources of depth uncertainties: (i) the range resolution of the radar system, (ii) the uncertainty associated with the value for the firn correction, and (iii) the uncertainty in the chosen dielectric permittivity. From the obtained total depth uncertainties we derive an age uncertainty by determining the lower and upper age bound of the respective IRH age at the EDML ice core chronology.

The range resolution of the EMR system in the short-pulse mode is approximately 5 m (Nixdorf et al.1999). For the UWB, the range resolution is determined by the bandwidth and the windowing factor in processing what results in a range resolution for the narrowband product of 4.3 m (Franke et al.2022b) and for the broadband product of 0.35 m (Koch et al.2023). For determining the uncertainty in the firn correction, we use the typically applied values in the range of 10–15 m. The values are uncertain due to the paucity of measurements and can vary spatially due to different rates of firn densifications over large regional scales. Previous studies used uncertainty values ranging from 0.5–2 m in East Antarctica (Cavitte et al.2016; Winter et al.2019; Wang et al.2023a) and 3 m in West Antarctica (Ashmore et al.2020a; Bodart et al.2021b). Therefore, we estimate an uncertainty of ± 2.5 m for our firn correction. For the uncertainty in the dielectric permittivity of dry ice, we apply an error of 1 % (Bohleber et al.2012). The total IRH depth uncertainty is the root-sum-square error of the three uncertainties and is calculated for each IRH data point.

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f05

Figure 5Colour-coded absolute depth representation of the nine IRHs from this study. The background map displays the ice surface elevation Reference Elevation Model of Antarctica (REMA) (Howat et al.2019) in hillshade. The fine black lines delineate the drainage basins according to Rignot et al. (2019).

2.7 Data validation

We validated our traced IRHs using a cross-point analysis for each IRH. To achieve this, we calculated crossover points along the traced horizons and determined the depth difference between two intersecting IRHs (see Appendix A for details). This ensured that the same IRHs were traced across intersecting transects across our study area. In total, we analysed more than 2000 crossover points, and the differences for all IRHs in the crossover analysis range between 0 and 74 m, with a mean standard deviation of 3.75 m. Since we traced the ice surface at the highest gradient of surface return power, while the IRHs were traced at their maximum return power, a certain fluctuation in depth at the cross-points is expected. The individual results of the cross-point analysis and details on the method are presented in Appendix A.

3 Results

3.1 IRHs ages in EMR short-pulse and UWB data

Two different dating methods provided IRH ages ranging from 4.8 to 91.0 ka (Table 2) covering the Holocene and Last Glacial Period. We use the age derived from the depth of conductivity peaks to date the IRHs due to their greater accuracy and only use the age from TWT–depth conversion in the radar profiles when no conductivity peak could be assigned. The radar reflections of seven of nine IRHs could be clearly associated with conductivity peaks in the EDML core using radar forward modelling (Fig. 2).

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f06

Figure 6Colour-coded normalized depth representation of the nine IRHs from this study. The background map displays the ice surface elevation REMA (Howat et al.2019) in hillshade. The fine black lines delineate the IMBIE drainage basins (Rignot et al.2019).

3.2 Spatial extent of IRHs in DML

The overall coverage of traceable IRHs in this study spans the entire western DML, covering an area of approximately 450 000 km2. IRHs extend over 700 km south of the EDML ice core,  200 km north, and  400 km southwest and northeast. IRHs also extend along the divide between EDML and Dome Fuji ice core sites into a region where IRHs were traced extensively in AWI's EMR long-pulse data (Fig. 4a; Wang et al.2023a). The horizontal density of traced IRHs varies significantly (Fig. 4c–e) and is highest within a radius of about 30 km around the EDML ice core (Fig. 4d). Additionally, there is dense coverage in the Jutulstraumen drainage basin south of the Maud Belt, particularly at the onset of the Jutulstraumen ice stream (Fig. 4c). Furthermore, there is a high density of IRHs southwest of the EDML ice core at the divide between the Slessor, Bailey, and Recovery drainage basins and the drainage basin feeding the Riiser-Larsen and Brunt Ice Shelf (Fig. 4e). Generally, IRHs cover mostly slow-moving regions, in particular at the ice divides. Regions of higher flow velocities are covered at the onset of Slessor ice stream ( 35 m yr−1) and Jutulstraumen Glacier (<100 m yr−1).

IRHs could be consistently traced in most radar profiles of the different radar products (Fig. 3). However, there are gaps, particularly in the area of the Jutulstraumen ice stream, which increase downstream as the internal layers become more folded due to ice flow convergence and the deep canyon system (Franke et al.2021). Strongly dipping internal layers show a loss in return power depending on their dip angle at depth due to off-nadir ray path losses and destructive interference (Holschuh et al.2014). Moreover, gaps in the 4.8 and sometimes in the 7.2 ka IRH occur because the surface multiple reflection overlays these IRHs, which makes continuous tracing difficult.

The IRH overlap is greatest in the vicinity of the EDML ice core, where all nine IRHs could be traced nearly continuously (Fig. 4b). Here, the layer continuity is least disturbed by ice flow, and both shallow and deep IRHs are well resolved in all radar systems. There is also high IRH overlap at the onset of the Jutulstraumen ice stream and along the ice divide towards Dome Fuji. The lowest overlap is generally found in the part of radar profiles closest to the Maud Belt, where deeper and older IRHs become absent.

In our region, IRHs with ages of 38.1 and 73.4 ka were previously traced and published by Winter et al. (2019) and Steinhage et al. (2013). The most extensive of the two IRH archives, published by Winter et al. (2019), is now significantly expanded by several hundreds of kilometres around EDML. The 38.1 and 73.4 ka IRHs are extended by  450 % and  480 %, respectively. Additionally, seven more IRHs were added with comparable coverage (Fig. 4b). Based on the existing data set of IRHs from Winter et al. (2019), this represents a total increase in IRH data coverage by a factor of 20.

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f07

Figure 7Depth distribution of the 38.1 ka IRH in the catchment area of the Jutulstraumen (a–c) and upstream of the Slessor ice stream (d–e). Background maps in panels (a) and (d) show ice flow velocity (Mouginot et al.2019a), panels (b) and (e) show bed elevation (Morlighem et al.2020), and panels (c) and (f) show ice thickness (Morlighem et al.2020) with RACMO (Regional Atmospheric Climate Model) snow accumulation represented as contour lines (Lenaerts et al.2013).

3.3 Spatial variation in IRH depth distribution

Our IRHs, which range between 4.8 and 91.0 ka in age, are found at minimum and maximum depths of 200 and 2200 m below the ice surface, respectively (Fig. 5). Both the absolute depth below the ice surface and normalized depth distribution of the IRHs show significant variation. We observe a general pattern where IRHs around the ice divide near EDML and to the west are deeper compared to IRHs further south and towards Dome Fuji (Figs. 5 and 6). This pattern is even more pronounced in the normalized IRH depths along the ice divide between EDML and Dome Fuji. Additionally, we note a contrasting trend when comparing absolute and relative IRH depths at the edges of the Maud Belt and on the plateau. For the 73.4 ka IRH, for example, it is evident that northeast of EDML, this IRH lies at depths of approximately 1000 to 1200 m, similar to the southeastern part near Dome Fuji (comparison of Feature A and Feature B in Fig. 5). However, the relative depth reveals that the 73.4 ka IRH is about 65 % below the surface northeast of EDML but only about 45 % near Dome Fuji. This pattern is also evident in the 11.1, 25.6, and 38.1 ka IRHs but with smaller differences as the IRHs become younger (comparison of Feature A and Feature B in Fig. 6).

The IRH depths within the Jutulstraumen drainage basin are particularly variable and show the following pattern in their relative depths. In the southwestern part of the Jutulstraumen onset region, the IRHs are significantly deeper (absolute depth) than in the northeastern area. The IRHs in this region (denoted as Feature A in Figs. 5 and 6) are significantly shallower compared to the southwest and east but show similar relative depths to EDML. This general pattern is evident in all IRHs that cover this region. A different pattern emerges when considering absolute depths. We find a greater variation with very shallow IRHs in the northeast and very deep IRHs in the Jutulstraumen Trough area. In the 38.1 ka IRH, these depth differences amount to up to 1500 m.

Additionally, the variation in depth of the IRHs is striking in two profiles extending from EDML towards the south-southeast. Around 79.9° S, there is a significant increase in absolute depth followed by a decrease to more shallow depths, observable in all IRHs from 7.2 ka onward (Feature C in Fig. 5). This pattern is again more pronounced in older IRHs compared to younger IRHs. However, there is no information on the relative IRH depth at this location since the ice base reflection is not visible in this section of the radar profile. Moreover, we observe an increase in absolute IRH depth at the southwestern edge of our IRH data, which is particularly pronounced in the 25.6 and 38.1 ka IRHs (Feature D in Fig. 5). Notably, older IRHs of this particular region could not be traced in the radar data.

4 Discussion

4.1 Interpretation of IRH depth distribution variations

Bed topography undulations, the overall accumulation patterns and complex ice flow dynamics are among the most likely reasons why IRH depths in DML vary widely across different subregions. The overall snow accumulation pattern in East Antarctica (Lenaerts et al.2013), velocity field (Mouginot et al.2019a), and ice thickness distribution (Morlighem et al.2020) agree on the broader scale with IRH depths, e.g. deep IRHs on the plateau and shallower IRHs towards the coast. Similar to the behaviour seen, for example, over the Gamburtsev Mountains (Sanderson et al.2024) or near Jutulstraumen Glacier (Franke et al.2021), IRHs in DML likely show variability in depth where bedrock is uneven, particularly near subglacial ridges and troughs.

Feature A (Figs. 5, 6, and 7) represents a notable anomaly in a region characterized by particularly shallow IRHs in both absolute and relative depths. This region deviates from the general trend of having a gradient from shallower IRHs on the plateau to deeper IRHs (especially in relative depth) near the Maud Belt. This general trend is likely linked to the overall snow accumulation pattern, with a higher accumulation rate towards the Antarctic coast and a lower accumulation rate on the plateau, as evidenced by the RACMO2 (Regional Atmospheric Climate Model; Lenaerts et al.2013) snow accumulation data (Fig. 7c) and by firn core analysis (Rotschky et al.2007). The fact that this pattern is also evident in older and thus deeper IRHs suggests that this low-accumulation anomaly may have persisted in this region not only during the Holocene but also during the Last Glacial Period, given that the ice in this region shows very low flow velocities and relatively undisturbed layering throughout the ice column.

Additionally, IRHs are particularly deep (both absolute and relative) in the southwestern area of the Jutulstraumen onset (Figs. 5 and 6), which could be linked to the ice stream activity and bed topography (Fig. 7a and b). Besides dynamic thinning due to acceleration of ice flow, increased friction at the ice base could lead to basal melting, causing the IRHs to subside. However, our work also shows that ice older than 73.4 ka is present in the main trunk of Jutulstraumen, at least at its onset. Further northeast, the shallow depths of the IRHs might be explained not only by the low-accumulation anomaly but also by reduced basal sliding, less basal melting, and even a positive basal mass balance due to basal freeze-on of subglacial water (Franke et al.2021, 2024b).

A comparison of IRH depths between Feature A (south of the Maud Belt) and Feature B (near Dome Fuji) highlights the importance of considering the relative depth of IRHs in relation to ice thickness. The fact that all IRHs traced near Dome Fuji are significantly shallower relative to ice thickness than those further northwest can be explained by the ice sheet's geometry and surface mass balance. Ice and firn cores indicate that both regions are characterized by low accumulation rates (Rotschky et al.2007; Oyabu et al.2023). This low-accumulation pattern is confirmed in the wider Dome Fuji area via shallow-sounding radar measurements in combination with firn cores and snow pits (Fujita et al.2011; Van Liefferinge et al.2021; Cavitte et al.2023). However, the ice at Dome Fuji is much thicker and significantly older at greater depths compared to further north. This results in lower relative depths for the IRHs. The low-accumulation pattern at Dome Fuji likely persisted over time (hundreds of thousands of years), impacting the long-term stratigraphy.

The sudden IRH depth increase of all ages along the southward-trending radar profiles (Feature C) and some IRHs in the southwest (Feature D) could be explained by changes in bed elevation and ice dynamics, respectively. Although there are no visible bed reflections in the radar profiles where the IRHs were traced, ice thickness data from nearby regions suggest a connection with the underlying topography. The bed topography data from BedMachine (Morlighem et al.2020) indicate that the IRH depth increase in Feature C is mostly associated with a subglacial valley with a width between 15–25 km and not with snow accumulation or flow velocity anomalies (Fig. 7d–f). In contrast, the IRHs at Feature D are located at the onset of one of the tributaries of the Slessor ice stream, where faster flow dominates. Higher flow velocity at this point would dynamically thin the ice, which must be compensated by lateral mass influx or increased snow accumulation, which would increase IRH depths. Hence, the relative drop in elevations of the IRHs in this area could indicate the long-term activity of this tributary (Fig. 7d), as previously suggested by Rippin et al. (2006). The fact that older IRHs, such as the 52.2 and 73.4 ka, could not be traced in this area further supports the hypothesis of a significant dynamic component causing IRH subsidence or layer disruption, which could suggest that this ice stream system has been stable at least for the past  70 ka.

4.2 Link between IRHs and volcanic eruptions

The dominant origin of IRHs was attributed to changes in conductivity due to variations in acidity originating from deposits of large volcanic events (Millar1981). Six of the IRHs that we dated using conductivity peaks and one IRH dated using the EDML ice core stratigraphy can be associated with major bipolar volcanic eruptions during the Holocene and Last Glacial Period (Lin et al.2022; Svensson et al.2013). Our 11.1 ka IRH, associated with a moderate conductivity peak, can be correlated with the deposits of a bipolar eruption event dated to 11.3 ka. These deposits have been identified in both Antarctic and Greenland ice cores (Lin et al.2022). Additionally, we consider the possibility that our 15.4 ka IRH could be linked to an event dated to 15.6 ka. Although there is no clear and strong conductivity peak in the EDML ice core for this event, one explanation could be that this event is classified as a relatively weaker bipolar eruption compared to others (rank 32; Lin et al.2022).

The radar reflection of the 25.6 ka IRH can be linked to sulfate deposits from the New Zealand Taupo Oruanui eruption, 25.32 ka, which were detected in two Greenland (although undetected in the GISP2 ice core) and three Antarctic ice cores (Lin et al.2022) and confirmed by tephra deposits in Antarctica (Dunbar et al.2017). The IRH we dated to 38.1 ka is very likely associated with the third strongest bipolar eruption in the Last Glacial Period known to date, which occurred at 38.13 ka. The exact eruption location is unknown, but it is likely a low-latitude or Southern Hemisphere eruption with an estimated average climate forcing twice as strong as the Tambora eruption in 1815 CE (Lin et al.2022). Our 52.2 ka IRH could correspond to an eruption that occurred at 52.3 ka and is listed as the fifth strongest eruption in Lin et al. (2022). The deposits from this eruption can be traced in three Greenland and three Antarctic ice cores (Lin et al.2022), and the possible eruption site is estimated to be in the high-latitude Northern Hemisphere (above 40° N).

From tracing IRHs in the radar data in this region, we find that the reflection of the 73.4 ka IRH is the most discernible and extensive in the deeper third of the ice sheet. As per Winter et al. (2019), we also speculate that this reflection most likely corresponds to the deposits of the Toba eruption in northwest Indonesia, dated using tephra and Ar–Ar dating to 73.88 ± 0.32 ka and 75.0 ± 0.9 ka, depending on the location of the deposits (see references in Svensson et al.2013). The Toba eruption is considered the largest known supereruption of the last 2.5 million years (Chesner2012) and likely occurred during the cooling transition from Greenland Interstadial 20 to Greenland Stadial 20 (Lin et al.2023, and references therein).

Although our IRH ages align closely with major volcanic events, minor age discrepancies are inherent to the dating methods for the IRHs and volcanic events, respectively. Inconsistencies can arise depending on the specific dating method for volcanic events (e.g. the different ages for the Toba eruption; Svensson et al.2013) as well as due to atmospheric transport and deposition processes of the volcanic ash and sulfate deposits. Additionally, ages assigned to IRHs are core-specific and vary between ice cores (Steinhage et al.2013; Winter et al.2019).

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f08

Figure 8IRHs from published data sets whose ages match with the (a) DML 4.8 ka horizon, (b) DML 38.1 ka horizon, (c) DML 73.4 ka horizon, and (d) DML 91.0 ka horizon. References to the scientific studies and data publications are provided in Table 3.

Beem et al. (2021b)Beem et al. (2021a)Steinhage et al. (2013)Ashmore et al. (2020a)Ashmore et al. (2020b)Bodart et al. (2021b)Bodart et al. (2021a)Muldoon et al. (2018)Muldoon et al. (2023)Steinhage et al. (2013)Ashmore et al. (2020a)Ashmore et al. (2020b)Bodart et al. (2021b)Bodart et al. (2021a)Siegert and Payne (2004)Beem et al. (2021b)Beem et al. (2021a)Chung et al. (2023b)Mulvaney et al. (2023)Muldoon et al. (2018)Steinhage et al. (2013)Beem et al. (2021b)Beem et al. (2021a)Steinhage et al. (2013)Cavitte et al. (2021)Cavitte et al. (2020)Chung et al. (2023b)Mulvaney et al. (2023)Leysinger Vieli et al. (2011)Sanderson et al. (2024)Sanderson et al. (2023a)Wang et al. (2023a)Wang et al. (2023b)Winter et al. (2019)Winter et al. (2018)Beem et al. (2021b)Beem et al. (2021a)Beem et al. (2021b)Beem et al. (2021a)Steinhage et al. (2013)Cavitte et al. (2021)Cavitte et al. (2020)Chung et al. (2023b)Chung et al. (2023a)Chung et al. (2023b)Mulvaney et al. (2023)Wang et al. (2023a)Wang et al. (2023b)Leysinger Vieli et al. (2011)Beem et al. (2021b)Beem et al. (2021a)Chung et al. (2023b)Chung et al. (2023a)Sanderson et al. (2024)Sanderson et al. (2023a)Winter et al. (2019)Winter et al. (2018)

Table 3IRHs from published data sets whose ages match with the IRHs of this study in western DML and likely represent the same physical reflector. The following abbreviations apply to the “Region” column: East Antarctica (EA), Central East Antarctica (CEA), Central West Antarctica (CWA), Dronning Maud Land (DML), Institute and Möller (I&M), and South Pole (SP).

a The age and error for the respective IRHs from Ashmore et al. (2020a) are based on the dating at the WAIS Divide ice core by Bodart et al. (2021b), with vertical conformity between the two data sets shown in Bodart et al. (2021b). b The age and error for the respective IRH from Winter et al. (2019) are based on the dating at the Dome C ice core.

Download Print Version | Download XLSX

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f09

Figure 9Potential to extrapolate IRHs to AWI EMR long-pulse data. (a) Overview of EMR short-pulse data (red) superimposed with the IRHs traced in this study (black) and in the background showing the EMR long-pulse data (dark blue). The IRHs from Wang et al. (2023a) are depicted in purple. The lines A–A (black) and B–B (red) in panel (a) indicate a long-pulse and short-pulse radargram at Dome Fuji (Fig. 10) where IRHs from this study and IRHs from Wang et al. (2023a) intersect. (b) Composite radargram (automatic gain control; AGC version) from EMR short-pulse and long-pulse data (profile ID 19992116, long-pulse, right; 19993116, short-pulse, left). At the intersection, transitions where IRHs from this study are also represented in the long-pulse data are marked with arrows.

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f10

Figure 10Illustration of two radargrams (automatic gain control; AGC version) and the IRHs intersecting at the Dome Fuji ice core. Panel (a) shows the radargrams with IRHs, and panel (b) shows the radargrams without IRHs but where the Wang et al. (2023a) IRHs intersect this study and Winter et al.'s 2019 IRHs at the Dome Fuji ice core site. (a) The radargram on the left is an EMR long-pulse profile with IRHs from Wang et al. (2023a), while the radargram on the right is an EMR short-pulse profile with IRHs from this study together with the IRHs from Winter et al. (2019). (b) The location at the Dome Fuji core, where the 36.3 and 75.3 ka IRH from Wang et al. (2023a) intersects with the 38.1 and 73.4 ka IRH from this study and Winter et al. (2019) are marked with arrows. The colours of the IRHs in the left radargram correspond to those in Fig. 2 in Wang et al. (2023a). The location of the two profiles is shown in Fig. 9a.

Download

4.3 Comparison with IRHs from other studies

Steinhage et al. (2013) traced eight IRHs along the ice divide connecting EMR short-pulse profiles between EDML and Dome Fuji and dated them via TWT-to-depth conversion using the EDML AICC2012 age–depth chronology (Bazin et al.2013). Several of these isochrones correspond with those we have traced over a much larger spatial extent in this study (corresponding IRH ages dated at the EDML core by Steinhage et al.2013, shown in parentheses): 4.8 ka (4.7 ka), 7.2 ka (7.4 ka), 25.6 ka (25.1 ka), 38.1 ka (38.1 ka), and 73.4 ka (72.4 ka) (Table 3). The discrepancies arise due to a combination of the choice of radar profiles, different dating methods, the type of surface pick, and the selected propagation speed for radar waves in ice, all of which are within the stated error margins. Moreover, the IRHs in Steinhage et al. (2013) were also dated using the Dome Fuji age–depth chronology, resulting in different ages. In addition to the study from Steinhage et al. (2013), the 38.1 and 73.4 ka IRHs had also already been traced extensively in our study area by Winter et al. (2019) (Fig. 8 and Table 3), and their paper also discusses briefly the spatially limited presence of IRHs with ages of 4.8, 7.6, 15.4, and 25.0 ka. Together, these data provide an essential foundation for this study and enable a direct and independent comparison of both tracing and dating methods.

A comparison with traced IRHs in East Antarctica, in the region of the Vostok ice core and EPICA Dome C ice core (Leysinger Vieli et al.2011; Cavitte et al.2016, 2021; Winter et al.2019; Chung et al.2023b), suggests that some of these IRHs likely correspond to those identified in our study (Table 3). Most of these isochrones extend far beyond an age of 100 ka, reaching up to approximately 700 ka, with the dating of the oldest IRHs being particularly challenging due to the inherent uncertainties associated with the 1-D age–depth model used to date (Cavitte et al.2021). The 38.1 and 73.4 ka IRHs (dated to the same ages by Cavitte et al.2021) are almost certainly the same reflections in this region, as they appear as very strong reflections in our radar data and correlate with significant past volcanic eruptions in ice cores (Lin et al.2022; Svensson et al.2013). Furthermore our 11.1, 38.1 and 74.3 ka IRHs also likely correspond to the 10.5, 39.2, and 75.2 ka IRHs from Chung et al. (2023b) at Dome C traced in ground-based DELORES (Deep Looking Radio Echo Sounder) data acquired in the 2016/17 and 2017/18 seasons as well as two IRHs (dated to 73.7 and 90.3 ka) traced in radar data acquired with the UA LDC-VHF system (Lilien et al.2021) at Little Dome C (Chung et al.2023b, Table 3).

Our isochrones also extend into the area around Dome Fuji, for which there are dated IRHs from Wang et al. (2023a) (Fig. 8). However, these were traced using EMR long-pulse data, which have a vertical resolution of  50 m, making the direct comparison more difficult. Nonetheless, they also identified IRHs dated to 36.3 ± 3.6 and 75.3 ± 7.0 ka, which, when examined at the intersection points with our EMR short-pulse radar data, most likely represent the same IRHs (Figs. 9 and 10). A large amount of our IRH ages correlate with IRHs from Beem et al. (2021b) at Titan Dome, particularly their 4.7, 10.7, 37.8, 51.4, 72.5, and 93.9 ka IRHs. Moreover, we see a connection of our 38.1 and 91.0 ka IRHs to the 38.2 ± 2.0 and 90.4 ± 3.6 ka IRHs observed in radar data between South Pole and Dome A (Sanderson et al.2024) acquired with the British Antarctic Survey's Polarimetric radar Airborne Science Instrument (PASIN and PASIN 2) during the AGAP and PolarGAP surveys (Fig. 8 and Table 3).

Dated IRHs with similar ages to ours are not confined to East Antarctica. Several studies have traced, dated, and published IRHs in West Antarctica (Muldoon et al.2018; Ashmore et al.2020a; Bodart et al.2021b, 2023). In particular, all of these studies identified a significantly widespread IRH dated at 4.7 ka that is present across most of West Antarctica (e.g. Bodart et al.2023) and which matches our shallowest IRH over DML (Fig. 8 and Table 3). For instance, Muldoon et al. (2018) identified a 4.7 ka IRH across West Antarctica's ice divide and throughout the Thwaites Glacier catchment, and Bodart et al. (2021b, 2023) identified a 4.72 ka IRH across the divide and throughout the Pine Island Glacier catchment, with direct links showing vertical conformity with another IRH widely traced across the Institute and Möller ice stream catchment by Ashmore et al. (2020a). Additionally, there may be further links between our 7.2 ka IRH and the 6.94 ka IRH from Ashmore et al. (2020a) and Bodart et al. (2021b), respectively, and our 25.6 ka IRH and the 24.9 ka IRH from Muldoon et al. (2018) (Fig. 8 and Table 3). However, to accurately determine if these dated isochrones indeed originate from the same reflector Antarctic-wide, it would be necessary to link not only the IRH ages (which can vary depending on the ice core, method, and determined TWT between the ice surface and the IRH) but also the backscatter patterns in the radar data. Here, radar lines connecting data sets between East and West Antarctica would be beneficial (Fig. 8). However, this presents challenges, such as long flights over the Antarctic Plateau, and geographical difficulties, such as the loss of layer continuity across the Transantarctic Mountains. For a definitive connection of IRHs that relate to the same reflectors and, therefore, the same deposits, geochemical analyses of ice cores could provide additional constraints.

The AntArchitecture initiative (Bingham et al.2024) has been aiming to compile all such IRH data sets into a continent-wide radiostratigraphic database, which will be used to create a three-dimensional age–depth model of the ice sheet by interpolating between known IRH along-track data sets similar to already existing gridded topographic models such as Bedmap (Lythe and Vaughan2001; Fretwell et al.2013; Frémand et al.2022; Pritchard et al.2025) and BedMachine Antarctica (Morlighem et al.2020). For such a task to be successful in the future, additional IRH tracing will need to take place in poorly surveyed areas of the ice sheet or away from ice divides and ice core locations using existing RES coverage to avoid over-fitting the interpolation algorithms. By utilizing 3 decades of AWI RES data from multiple radar systems in the area, this study connects a significant gap in the East Antarctic radiostratigraphic record for IRHs older than 7.2 ka and provides additional evidence of the presence of West Antarctica's 4.7–4.8 ka IRH deep in the East Antarctic ice sheet. The success of identifying similar IRHs across large swaths of both East and West Antarctica as demonstrated here (Fig. 8) invites further opportunities to make use of the available historical record of yet unused radar data across Antarctica to accelerate the retrieval of a continent-wide age–depth model of the ice sheet, as envisaged by AntArchitecture (Bingham et al.2024).

4.4 Potential to expand IRH coverage in and beyond DML

In our study, we have not fully utilized the complete archive of AWI radar data. We focused on EMR short-pulse data, UWB narrow data, and UWB wideband data because these have similar vertical resolutions (5 m or better), enabling the capture of a higher number of IRHs as well as shallow IRHs that are comparably represented across these systems. However, there are additional areas in DML covered exclusively by EMR long-pulse data (approximately 50 m vertical resolution but deep sounding). At intersections between EMR short-pulse and long-pulse data where we were able to trace IRHs, it becomes apparent that the 25.6, 38.1, 52.2, and 73.4 ka IRHs are also represented in the long-pulse data (Fig. 9), which aligns with the IRHs of similar age in Wang et al. (2023a) (Fig. 8 and Table 3). However, the range resolution of the long-pulse system used in Wang et al. (2023a) is significantly lower, resulting in a larger uncertainty in age dating compared to our IRHs. A direct comparison of two radargrams intersecting at Dome Fuji confirms that particularly the 36.3 and 75.3 ka IRHs from Wang et al. (2023a) correspond to the same reflections identified in our study, which we date as 38.1 and 73.4 ka, respectively (Fig. 10). Additionally, it is likely that the reflection of the 95.5 ka IRH from Wang et al. (2023a) corresponds to our 91.0 ka IRH (Fig. 10). However, the offset in IRH depth and age is relatively large.

The representation of the same reflectors in EMR long-pulse and short-pulse data does not only account for reflections caused by changes in electrical conductivity but also for reflections caused by changes in ice crystal lattice orientation (Eisen et al.2007). Considering the AWI EMR radar line coverage in Fig. 9a, there is additional potential, particularly for regions upstream of Princess Ragnhild and Princess Astrid coasts, as well as the onsets of the Slessor and Recovery ice streams. Additionally, the spatial coverage of the East Antarctic Plateau between EDML and Dome Fuji could be significantly improved. It may even be possible to establish connections between our IRHs and the western Coats Land at the ice divide to the Bailey and Slessor basins using AWI's and the British Antarctic Survey's radar data (Frémand et al.2022), which are of similar range resolution.

4.5 Significance of dated IRHs for reconstructing the ice-sheet history in western DML

Reconstructing past ice-sheet configurations is crucial for understanding ice-sheet processes and their impacts on the Earth system, as well as providing essential constraints for ice-sheet models. The spatial depth distribution of our nine IRHs in western DML can significantly contribute to understanding ice-sheet evolution. By comparing IRH depths and spatial anomalies, we can test hypotheses about underlying mechanisms, such as surface or basal mass balance and ice-dynamic processes, as well as how they may have changed over time. Our dated IRHs also provide an opportunity to validate or expand existing findings and hypotheses on ice surface changes based on cosmogenic nuclides (Andersen et al.2020) about past ice dynamics in DML (Andersen et al.2023; Braga et al.2023). In addition, the absolute and relative depth of our IRHs offer the opportunity to complement studies on, for example, preserved paleo-geomorphological structures (e.g. Näslund1997; Rose et al.2015; Franke et al.2021; Carter et al.2024) or regions of basal freeze-on of subglacial water (e.g. Bell et al.2011; Leysinger Vieli et al.2018; Franke et al.2024b) providing further insights into past flow behaviour and ice-sheet configurations.

Furthermore, we see significant potential for improving regional and Antarctic-wide ice-sheet models to provide better projections for future sea-level rise (Sutter et al.2021). Indeed, IRHs have previously been used as a tool for calibrating results in model simulations (Sutter et al.2021; Višnjević et al.2022) and would be well suited to testing hypotheses of the stability of the dynamic Slessor and Jutulstraumen ice stream onset regions over the last interglacial period (Rippin et al.2003, 2006). Conversely, IRH data from less dynamically active regions, such as the ice divide southwest from EDML, along the divide between EDML and Dome Fuji, the region south of EDML, and the shallow-IRH anomaly region (Feature A), are well suited for inferring the spatial distribution of past accumulation rates (Eisen et al.2005; Huybrechts et al.2009; Leysinger Vieli et al.2011; Bodart et al.2023).

An additional approach to decipher the past ice dynamics of a region using IRHs lies in the structural analysis of three-dimensional englacial structures, particularly the geometries of folds and fold axial planes, as well as the development of fold amplitudes with depth (Bons et al.2016; Franke et al.2023a; Jansen et al.2024). Especially for Greenland, it has been shown in various dynamic ice regions that a three-dimensional representation aiming to depict deformation structures like folds accurately can contribute to understanding ice mechanical properties (Bons et al.2016; Zhang et al.2024) and past flow patterns (Franke et al.2022a; Jansen et al.2024). Notably, the profiles in the upstream region of the Jutulstraumen ice stream, as well as the region around EDML, where the high density of IRHs allows resolving small-scale englacial structures, and the profiles at the ice divides have the potential to provide insights into past ice dynamics, such as ice flow direction changes or shifts in the location of ice divides.

5 Conclusions

We have mapped nine internal reflection horizons (IRHs) using radar data from various systems of the Alfred Wegener Institute, collected over the last 3 decades. Our IRHs cover the western DML, south of the Maud Belt, and range from 4.8 to 91.0 ka, covering the Holocene and Last Glacial Period. Accurate dating of the IRHs was achieved by combining the age–depth chronology of the EPICA DML ice core and DEP-based radar forward modelling. Additionally, six of the new IRHs could be linked to deposits from significant historic bipolar volcanic eruptions, facilitating synchronization of these reflectors across Antarctica and potentially Greenland. Many of the IRHs mapped in this study likely correspond to the same englacial reflectors found in extensive regions over East and West Antarctica. A comparison with AWI EMR long-pulse data suggests that some of the IRHs identified here could be extended to even larger sectors of East Antarctica.

Our results show that the influence of bed topography, accumulation rates, and ice dynamics is engraved in the age–depth stratigraphy of the ice in the western DML. Thus, our study significantly contributes to a broader understanding of the englacial age architecture of the DML region in East Antarctica and highlights the potential for linking individual IRHs to other regions in Antarctica. Furthermore, the findings presented here are fundamental to enhancing our comprehension of past ice-sheet processes and are crucial for numerical ice flow models aimed at improving our understanding of the paleoclimatic history of Antarctica.

Appendix A: Cross-point analysis

We performed an analysis of IRH depth differences at crossover points of intersecting radar lines to validate that we assigned the same reflections for the respective IRHs. For the analysis of crossover differences in the respective IRHs, we calculate the exact intersection points between the geometric lines formed by connecting the individual IRH picks. We create line subsets of one radar profile where IRH gaps exist to avoid creating lines, and therefore potential fake intersection points, where no IRHs are traced within a radar profile. Using these lines, we determine the location of the resulting intersection points. We then create a circular buffer around these intersection points and capture all IRH picks within this buffer. We select a buffer radius of 50 m and average the IRH depths of all picks within a profile to minimize small-scale variations introduced by the semi-automatic tracker. The difference is derived from the mean depths of picks from both intersecting lines within the buffer.

Examining the histograms in Fig. A1, we observe that most intersection point differences are below 10 m depth difference, and the depth difference generally decreases for higher values. For the 4.8 and 7.2 ka IRHs, we see a small local maximum around a 5 m intersection point difference. Moreover, the differences become larger for older IRHs, leading to a higher number of larger cross-point errors compared to younger IRHs.

(based on REMA; Howat et al.2019)

Table A1IRH data set description for column-separated text files.

 The ice base reflection is not present in all radar profiles.

Download Print Version | Download XLSX

https://tc.copernicus.org/articles/19/1153/2025/tc-19-1153-2025-f11

Figure A1Differences of IRH depths up to 50 m at crossing points for the respective IRHs. The histograms in the lower-left corner of each panel show the total count of differences of IRH depths at crossing points up to 50 m (x axis). The y axis shows the number of crossing points, and the bin width is 1 m.

Download

Appendix B: IRH data set description

We publish our dated IRHs on the PANGAEA data publisher (Felden et al.2023) in a column-separated text file format, and the key information is summarized in Table A1. It is worth noting that the vertical position of the IRHs is provided in different formats: (1) as absolute TWT in the radargram, representing the direct reference to the radargram, and (2) as TWT (in seconds) and depth (in metres) beneath our ice surface pick. To ensure this reference is clear for future users, we also provide the absolute TWT of the ice surface pick, the elevation of the ice surface (based on the REMA surface DEM version 2; Howat et al.2019) at each location, and the depth of the ice base if the bed reflection is visible in the radargram. For conversion between TWT to depth, we used a dielectric permittivity of εr=3.15 and a firn correction term of Zfirn=13 m (see Eq. 1). We also provide the depth uncertainty of each IRH depth as well as an estimate for the resulting IRH age error (see “Data and methods” section). The “profile_id” is the ID of the radargram, which also appears in the filename. The “paradigm_id” is an internal AWI ID for referencing the AWI radar database and radar processing software where the picks are created and archived.

Data availability

Our IRHs are available on PANGAEA (Franke et al.2024a): https://doi.org/10.1594/PANGAEA.973266. The AICC2023 EDML ice core age–depth chronology (Bouchet et al.2023) is available on PANGAEA: https://doi.org/10.1594/PANGAEA.961019. Ice surface velocities from Mouginot et al. (2019a) are available from the National Snow and Ice Data Center (NSIDC), https://doi.org/10.7280/D10D4Z (Mouginot et al.2019b). The IMBIE drainage system boundaries (Rignot et al.2019) can be obtained from http://imbie.org/imbie-3/drainage-basins/ (last access: 17 January 2025). The BedMachine Antarctica v03 bed topography data from Morlighem et al. (2020) are available at https://doi.org/10.5067/FPSU0V1MWUB6 (Morlighem2022). The Reference Elevation Model of Antarctica (REMA) is available from the Polar Geospatial Center at https://www.pgc.umn.edu/data/rema/ (last access: 17 January 2025) and https://doi.org/10.7910/DVN/EBW8UC (Howat et al.2022). RACMO2.3 snow accumulation data can be obtained at the Institute for Marine and Atmospheric research Utrecht at Utrecht University: https://www.projects.science.uu.nl/iceclimate/models/racmo-data.php (last access: 17 January 2025) and https://doi.org/10.5281/zenodo.10854319 (Van Dalum et al.2024). For the data availability of IRHs used for comparison, we refer to Table 3.

Author contributions

SF conceptualized the study, wrote the manuscript, traced and dated all IRHs, performed all data analysis, created all figures, and validated the data. AMZ, DS, and SF acquired radar data from the 2023/24 AWI UWB campaign, while AMZ processed the data from this campaign with contributions of SF and VH. AMZ, VH, and SF contributed to tracing the ice–bed interface. JAB provided context of comparable studies and for the discussion about connecting IRHs in this study over Antarctica. VH developed an algorithm for standard ice surface tracking and processed GPS data. OE, DS, SF, and VH implemented access to and usability of AWI's radar data archive. All authors jointly revised and edited the manuscript.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

Acknowledgements

We thank the AWI polar aircraft technicians Eduard Gebhard and Christoph Petersen for their support in the field during the 2023/24 radar campaign, the Kenn Borek crew of AWI's polar research aircraft, and all supporters of previous radar campaigns mentioned in acknowledgements of related publications. We acknowledge support via airborne radar campaign funding grants AWI_PA_02004, AWI_PA_02020, AWI_PA_02084, AWI_PA_02105, and AWI_PA_02146. Logistical support in the field over the past 3 decades has been provided by Neumayer III station (Germany), Troll station (Norway), Kohnen station (Germany), Princess Elisabeth station (Belgium), and Novolazarevskaya airbase (Russia). We acknowledge the use of software from Open Polar Radar generated with support from the University of Kansas; NASA grants 80NSSC20K1242 and 80NSSC21K0753; and NSF grants OPP-2027615, OPP-2019719, OPP-1739003, IIS-1838230, RISE-2126503, RISE-2127606, and RISE-2126468. The authors would like to thank Aspen Technology, Inc. for providing software licenses and support. This study was motivated by the AntArchitecture Scientific Committee on Antarctic Research Action Group (https://scar.org/science/cross/antarchitecture, last access: 17 January 2025).

Steven Franke was funded by the Walter Benjamin Programme of the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation; project no. 506043073). Alexandra M. Zuhr was funded by the DFG in the framework of the priority programme SPP 1158 “Antarctic Research with Comparative Investigations in Arctic Ice Areas” (grant no. 522419679). Julien A. Bodart acknowledges funding from the Swiss National Science Foundation (grant no. 211542).

Financial support

This research has been supported by the Deutsche Forschungsgemeinschaft (grant nos. 506043073 and 522419679) and the Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (grant no. 211542).

This open-access publication was funded by the Open Access Publication Fund of the University of Tübingen.

Review statement

This paper was edited by Joseph MacGregor and reviewed by Rebecca Sanderson and Marie G. P. Cavitte.

References

Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung: Polar aircraft Polar5 and Polar6 operated by the Alfred Wegener Institute, Journal of large-scale research facilities JLSRF, 2, 87, https://doi.org/10.17815/jlsrf-2-153, 2016a. a

Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung: Neumayer III and Kohnen Station in Antarctica operated by the Alfred Wegener Institute, Journal of large-scale research facilities JLSRF, 2, 85, https://doi.org/10.17815/jlsrf-2-152, 2016b. a

Andersen, J., Newall, J., Blomdin, R., Sams, S., Fabel, D., Koester, A., Lifton, N., Fredin, O., Caffee, M., Glasser, N. F., Rogozhina, I., Suganuma, Y., Harbor, J., and Stroeven, A.: Ice surface changes during recent glacial cycles along the Jutulstraumen and Penck Trough ice streams in western Dronning Maud Land, East Antarctica, Quaternary Sci. Rev., 249, 106636, https://doi.org/10.1016/j.quascirev.2020.106636, 2020. a

Andersen, J. L., Newall, J. C., Fredin, O., Glasser, N. F., Lifton, N. A., Stuart, F. M., Fabel, D., Caffee, M., Pedersen, V. K., Koester, A. J., Suganuma, Y., Harbor, J. M., and Stroeven, A. P.: A topographic hinge-zone divides coastal and inland ice dynamic regimes in East Antarctica, Commun. Earth Environ., 4, 9, https://doi.org/10.1038/s43247-022-00673-6, 2023. a

Ashmore, D. W., Bingham, R. G., Ross, N., Siegert, M. J., Jordan, T. A., and Mair, D. W. F.: Englacial Architecture and Age‐Depth Constraints Across the West Antarctic Ice Sheet, Geophys. Res. Lett., 47, e2019GL086663, https://doi.org/10.1029/2019gl086663, 2020a. a, b, c, d, e, f, g, h, i

Ashmore, D. W., Bingham, R. G., Ross, N., Siegert, M. J., Jordan, T. A., and Mair, D. W. F.: Radiostratigraphy of the Weddell Sea Sector of West Antarctica, Zenodo [data set], https://doi.org/10.5281/zenodo.3635940, 2020b. a, b

Bazin, L., Landais, A., Lemieux-Dudon, B., Toyé Mahamadou Kele, H., Veres, D., Parrenin, F., Martinerie, P., Ritz, C., Capron, E., Lipenkov, V., Loutre, M.-F., Raynaud, D., Vinther, B., Svensson, A., Rasmussen, S. O., Severi, M., Blunier, T., Leuenberger, M., Fischer, H., Masson-Delmotte, V., Chappellaz, J., and Wolff, E.: An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120–800 ka, Clim. Past, 9, 1715–1731, https://doi.org/10.5194/cp-9-1715-2013, 2013. a

Beem, L. H., Young, D. A., Greenbaum, J., Ng, G., Blankenship, D. D., Cavitte, M. G. P., Jingxue, G., and Bo, S.: Titan Dome, East Antarctica, Aerogeophysical Survey [data set], https://doi.org/10.15784/601437, 2021a. a, b, c, d, e, f

Beem, L. H., Young, D. A., Greenbaum, J. S., Blankenship, D. D., Cavitte, M. G. P., Guo, J., and Bo, S.: Aerogeophysical characterization of Titan Dome, East Antarctica, and potential as an ice core target, The Cryosphere, 15, 1719–1730, https://doi.org/10.5194/tc-15-1719-2021, 2021b. a, b, c, d, e, f, g, h

Bell, R. E., Ferraccioli, F., Creyts, T. T., Braaten, D., Corr, H., Das, I., Damaske, D., Frearson, N., Jordan, T., Rose, K., Studinger, M., and Wolovick, M.: Widespread Persistent Thickening of the East Antarctic Ice Sheet by Freezing from the Base, Science, 331, 1592–1595, https://doi.org/10.1126/science.1200109, 2011. a

Bingham, R. G., Bodart, J. A., Cavitte, M. G. P., Chung, A., Sanderson, R. J., Sutter, J. C. R., Eisen, O., Karlsson, N. B., MacGregor, J. A., Ross, N., Young, D. A., Ashmore, D. W., Born, A., Chu, W., Cui, X., Drews, R., Franke, S., Goel, V., Goodge, J. W., Henry, A. C. J., Hermant, A., Hills, B. H., Holschuh, N., Koutnik, M. R., Leysinger Vieli, G. J.-M. C., Mackie, E. J., Mantelli, E., Martín, C., Ng, F. S. L., Oraschewski, F. M., Napoleoni, F., Parrenin, F., Popov, S. V., Rieckh, T., Schlegel, R., Schroeder, D. M., Siegert, M. J., Tang, X., Teisberg, T. O., Winter, K., Yan, S., Davis, H., Dow, C. F., Fudge, T. J., Jordan, T. A., Kulessa, B., Matsuoka, K., Nyqvist, C. J., Rahnemoonfar, M., Siegfried, M. R., Singh, S., Višnjević, V., Zamora, R., and Zuhr, A.: Review Article: Antarctica's internal architecture: Towards a radiostratigraphically-informed age–depth model of the Antarctic ice sheets, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2024-2593, 2024. a, b, c, d, e, f

Bodart, J., Bingham, R., Ashmore, D., Karlsson, N., Hein, A., and Vaughan, D.: Dated radar stratigraphy of the Pine Island Glacier catchment (West Antarctica) derived from BBAS-PASIN (2004-05) and OIB-MCoRDS2 (2016/2018) surveys (Version 1.0), UK Polar Data Centre, Natural Environment Research Council, UK Research & Innovation [data set], https://doi.org/10.5285/f2de31af-9f83-44f8-9584-f0190a2cc3eb, 2021a. a, b

Bodart, J. A., Bingham, R. G., Ashmore, D. W., Karlsson, N. B., Hein, A. S., and Vaughan, D. G.: Age‐Depth Stratigraphy of Pine Island Glacier Inferred From Airborne Radar and Ice‐Core Chronology, J. Geophys. Res.-Earth Surf., 126, e2020JF005927, https://doi.org/10.1029/2020jf005927, 2021b. a, b, c, d, e, f, g, h, i, j, k

Bodart, J. A., Bingham, R. G., Young, D. A., MacGregor, J. A., Ashmore, D. W., Quartini, E., Hein, A. S., Vaughan, D. G., and Blankenship, D. D.: High mid-Holocene accumulation rates over West Antarctica inferred from a pervasive ice-penetrating radar reflector, The Cryosphere, 17, 1497–1512, https://doi.org/10.5194/tc-17-1497-2023, 2023. a, b, c, d, e, f

Bohleber, P., Wagner, N., and Eisen, O.: Permittivity of ice at radio frequencies: Part I. Coaxial transmission line cell, Cold Reg. Sci. Technol., 82, 56–67, https://doi.org/10.1016/j.coldregions.2012.05.011, 2012. a, b

Bons, P. D., Jansen, D., Mundel, F., Bauer, C. C., Binder, T., Eisen, O., Jessell, M. W., Llorens, M.-G., Steinbach, F., Steinhage, D., and Weikusat, I.: Converging flow and anisotropy cause large-scale folding in Greenland's ice sheet, Nat. Commun., 7, 11427, https://doi.org/10.1038/ncomms11427, 2016. a, b, c, d

Bouchet, M., Landais, A., Grisart, A., Parrenin, F., Prié, F., Jacob, R., Fourré, E., Capron, E., Raynaud, D., Lipenkov, V. Y., Loutre, M.-F., Extier, T., Svensson, A. M., Martinerie, P., Leuenberger, M. C., Jiang, W., Ritterbusch, F., Lu, Z.-T., and Yang, G.-M.: The Antarctic ice core chronology (AICC2023), PANGAEA [data set bundled publication], https://doi.org/10.1594/PANGAEA.961019, 2023. a

Bouchet, M., Landais, A., Grisart, A., Parrenin, F., Prié, F., Jacob, R., Fourré, E., Capron, E., Raynaud, D., Lipenkov, V. Y., Loutre, M.-F., Extier, T., Svensson, A., Legrain, E., Martinerie, P., Leuenberger, M., Jiang, W., Ritterbusch, F., Lu, Z.-T., and Yang, G.-M.: The Antarctic Ice Core Chronology 2023 (AICC2023) chronological framework and associated timescale for the European Project for Ice Coring in Antarctica (EPICA) Dome C ice core, Clim. Past, 19, 2257–2286, https://doi.org/10.5194/cp-19-2257-2023, 2023. a, b

Braga, M. M. e., Jones, R. S., Bernales, J., Andersen, J. L., Fredin, O., Morlighem, M., Koester, A. J., Lifton, N. A., Harbor, J. M., Suganuma, Y., Glasser, N. F., Rogozhina, I., and Stroeven, A. P.: A thicker Antarctic ice stream during the mid-Pliocene warm period, Commun. Earth Environ., 4, 321, https://doi.org/10.1038/s43247-023-00983-3, 2023. a

Carter, C. M., Bentley, M. J., Jamieson, S. S. R., Paxman, G. J. G., Jordan, T. A., Bodart, J. A., Ross, N., and Napoleoni, F.: Extensive palaeo-surfaces beneath the Evans–Rutford region of the West Antarctic Ice Sheet control modern and past ice flow, The Cryosphere, 18, 2277–2296, https://doi.org/10.5194/tc-18-2277-2024, 2024. a

Cavitte, M. G. P., Blankenship, D. D., Young, D. A., Schroeder, D. M., Parrenin, F., Lemuer, E., MacGregor, J. A., and Siegert, M. J.: Deep radiostratigraphy of the East Antarctic plateau: connecting the Dome C and Vostok ice core sites, J. Glaciol., 62, 323–334, https://doi.org/10.1017/jog.2016.11, 2016. a, b, c, d, e, f

Cavitte, M. G. P., Parrenin, F., Ritz, C., Young, D. A., Van Liefferinge, B., Blankenship, D. D., Frezzotti, M., and Roberts, J. L.: Accumulation patterns around Dome C, East Antarctica, in the last 73 kyr, The Cryosphere, 12, 1401–1414, https://doi.org/10.5194/tc-12-1401-2018, 2018. a

Cavitte, M. G. P., Young, D. A., Mulvaney, R., Ritz, C., Greenbaum, J. S., Ng, G., Kempf, S. D., Quartini, E., Muldoon, G. R., Paden, J., Frezzotti, M., Roberts, J. L., Tozer, C. R., Schroeder, D. M., and Blankenship, D. D.: A detailed radiostratigraphic data set for the central East Antarctic Plateau spanning from the Holocene to the mid-Pleistocene, Earth Syst. Sci. Data, 13, 4759–4777, https://doi.org/10.5194/essd-13-4759-2021, 2021. a, b, c, d, e, f, g

Cavitte, M. G. P., Goosse, H., Matsuoka, K., Wauthy, S., Goel, V., Dey, R., Pratap, B., Van Liefferinge, B., Meloth, T., and Tison, J.-L.: Investigating the spatial representativeness of East Antarctic ice cores: a comparison of ice core and radar-derived surface mass balance over coastal ice rises and Dome Fuji, The Cryosphere, 17, 4779–4795, https://doi.org/10.5194/tc-17-4779-2023, 2023. a

Cavitte, M. M., Blankenship, D. D., Frezzotti, M., Kempf, S. D., and Muldoon, G.: Ice-penetrating radar internal stratigraphy over Dome C and the wider East Antarctic Plateau, U.S. Antarctic Program (USAP) Data Center [data set], https://doi.org/10.15784/601411, 2020. a, b

Chesner, C. A.: The Toba Caldera Complex, Quaternary Int., 258, 5–18, https://doi.org/10.1016/j.quaint.2011.09.025, 2012. a

Chung, A., Parrenin, F., Steinhage, D., Mulvaney, R., Martín, C., Cavitte, M. G. P., Lilien, D. A., Helm, V., Taylor, D., Gogineni, P., Ritz, C., Frezzotti, M., O′Neill, C., Miller, H., Dahl-Jensen, D., and Eisen, O.: Internal reflecting horizons at Little Dome C, Antarctica, PANGAEA [data set], https://doi.org/10.1594/PANGAEA.957176, 2023a. a, b

Chung, A., Parrenin, F., Steinhage, D., Mulvaney, R., Martín, C., Cavitte, M. G. P., Lilien, D. A., Helm, V., Taylor, D., Gogineni, P., Ritz, C., Frezzotti, M., O'Neill, C., Miller, H., Dahl-Jensen, D., and Eisen, O.: Stagnant ice and age modelling in the Dome C region, Antarctica, The Cryosphere, 17, 3461–3483, https://doi.org/10.5194/tc-17-3461-2023, 2023b. a, b, c, d, e, f, g, h, i

Drews, R., Eisen, O., Weikusat, I., Kipfstuhl, S., Lambrecht, A., Steinhage, D., Wilhelms, F., and Miller, H.: Layer disturbances and the radio-echo free zone in ice sheets, The Cryosphere, 3, 195–203, https://doi.org/10.5194/tc-3-195-2009, 2009. a

Drews, R., Schannwell, C., Ehlers, T. A., Gladstone, R., Pattyn, F., and Matsuoka, K.: Atmospheric and Oceanographic Signatures in the Ice Shelf Channel Morphology of Roi Baudouin Ice Shelf, East Antarctica, Inferred From Radar Data, J. Geophys. Res.-Earth Surf., 125, e2020JF005587, https://doi.org/10.1029/2020jf005587, 2020. a

Dunbar, N. W., Iverson, N. A., Eaton, A. R. V., Sigl, M., Alloway, B. V., Kurbatov, A. V., Mastin, L. G., McConnell, J. R., and Wilson, C. J. N.: New Zealand supereruption provides time marker for the Last Glacial Maximum in Antarctica, Sci. Rep., 7, 12238, https://doi.org/10.1038/s41598-017-11758-0, 2017. a, b

Eisen, O., Wilhelms, F., Nixdorf, U., and Miller, H.: Revealing the nature of radar reflections in ice: DEP‐based FDTD forward modeling, Geophys. Res. Lett., 30, 1218, https://doi.org/10.1029/2002gl016403, 2003. a, b

Eisen, O., Rack, W., Nixdorf, U., and Wilhelms, F.: Characteristics of accumulation around the EPICA deep-drilling site in Dronning Maud Land, Antarctica, Ann. Glaciol., 41, 41–46, https://doi.org/10.3189/172756405781813276, 2005. a, b

Eisen, O., Wilhelms, F., Steinhage, D., and Schwander, J.: Improved method to determine radio-echo sounding reflector depths from ice-core profiles of permittivity and conductivity, J. Glaciol., 52, 299–310, https://doi.org/10.3189/172756506781828674, 2006. a, b, c, d, e, f, g, h

Eisen, O., Hamann, I., Kipfstuhl, S., Steinhage, D., and Wilhelms, F.: Direct evidence for continuous radar reflector originating from changes in crystal-orientation fabric, The Cryosphere, 1, 1–10, https://doi.org/10.5194/tc-1-1-2007, 2007. a, b

Eisermann, H., Eagles, G., Ruppel, A., Smith, E. C., and Jokat, W.: Bathymetry Beneath Ice Shelves of Western Dronning Maud Land, East Antarctica, and Implications on Ice Shelf Stability, Geophys. Res. Lett., 47, e2019GL086724, https://doi.org/10.1029/2019gl086724, 2020. a

Felden, J., Möller, L., Schindler, U., Huber, R., Schumacher, S., Koppe, R., Diepenbroek, M., and Glöckner, F. O.: PANGAEA – Data Publisher for Earth & Environmental Science, Sci. Data, 10, 347, https://doi.org/10.1038/s41597-023-02269-x, 2023. a

Franke, S., Eisermann, H., Jokat, W., Eagles, G., Asseng, J., Miller, H., Steinhage, D., Helm, V., Eisen, O., and Jansen, D.: Preserved landscapes underneath the Antarctic Ice Sheet reveal the geomorphological history of Jutulstraumen Basin, Earth Surf. Process. Landf., 46, 2728–2745, https://doi.org/10.1002/esp.5203, 2021. a, b, c, d, e, f

Franke, S., Bons, P. D., Westhoff, J., Weikusat, I., Binder, T., Streng, K., Steinhage, D., Helm, V., Eisen, O., Paden, J. D., Eagles, G., and Jansen, D.: Holocene ice-stream shutdown and drainage basin reconfiguration in northeast Greenland, Nat. Geosci., 15, 995–1001, https://doi.org/10.1038/s41561-022-01082-2, 2022a. a, b

Franke, S., Jansen, D., Binder, T., Paden, J. D., Dörr, N., Gerber, T. A., Miller, H., Dahl-Jensen, D., Helm, V., Steinhage, D., Weikusat, I., Wilhelms, F., and Eisen, O.: Airborne ultra-wideband radar sounding over the shear margins and along flow lines at the onset region of the Northeast Greenland Ice Stream, Earth Syst. Sci. Data, 14, 763–779, https://doi.org/10.5194/essd-14-763-2022, 2022b. a, b

Franke, S., Bons, P. D., Streng, K., Mundel, F., Binder, T., Weikusat, I., Bauer, C. C., Paden, J. D., Dörr, N., Helm, V., Steinhage, D., Eisen, O., and Jansen, D.: Three-dimensional topology dataset of folded radar stratigraphy in northern Greenland, Sci. Data, 10, 525, https://doi.org/10.1038/s41597-023-02339-0, 2023a. a, b

Franke, S., Gerber, T., Warren, C., Jansen, D., Eisen, O., and Dahl-Jensen, D.: Investigating the Radar Response of Englacial Debris Entrained Basal Ice Units in East Antarctica Using Electromagnetic Forward Modeling, IEEE T. Geosci. Remote, 61, 1–16, https://doi.org/10.1109/tgrs.2023.3277874, 2023b. a, b

Franke, S., Steinhage, D., Helm, V., Zuhr, A., and Eisen, O.: Internal Reflecting Horizons (IRHs) in western Dronning Maud Land, East Antarctica from airborne radar surveys (4.8–91.0 ka), PANGAEA [data set], https://doi.org/10.1594/PANGAEA.973266, 2024a. a

Franke, S., Wolovick, M., Drews, R., Jansen, D., Matsuoka, K., and Bons, P. D.: Sediment Freeze‐On and Transport Near the Onset of a Fast‐Flowing Glacier in East Antarctica, Geophys. Res. Lett., 51, e2023GL107164, https://doi.org/10.1029/2023gl107164, 2024b. a, b, c

Fretwell, P., Pritchard, H. D., Vaughan, D. G., Bamber, J. L., Barrand, N. E., Bell, R., Bianchi, C., Bingham, R. G., Blankenship, D. D., Casassa, G., Catania, G., Callens, D., Conway, H., Cook, A. J., Corr, H. F. J., Damaske, D., Damm, V., Ferraccioli, F., Forsberg, R., Fujita, S., Gim, Y., Gogineni, P., Griggs, J. A., Hindmarsh, R. C. A., Holmlund, P., Holt, J. W., Jacobel, R. W., Jenkins, A., Jokat, W., Jordan, T., King, E. C., Kohler, J., Krabill, W., Riger-Kusk, M., Langley, K. A., Leitchenkov, G., Leuschen, C., Luyendyk, B. P., Matsuoka, K., Mouginot, J., Nitsche, F. O., Nogi, Y., Nost, O. A., Popov, S. V., Rignot, E., Rippin, D. M., Rivera, A., Roberts, J., Ross, N., Siegert, M. J., Smith, A. M., Steinhage, D., Studinger, M., Sun, B., Tinto, B. K., Welch, B. C., Wilson, D., Young, D. A., Xiangbin, C., and Zirizzotti, A.: Bedmap2: improved ice bed, surface and thickness datasets for Antarctica, The Cryosphere, 7, 375–393, https://doi.org/10.5194/tc-7-375-2013, 2013. a

Frémand, A. C., Bodart, J. A., Jordan, T. A., Ferraccioli, F., Robinson, C., Corr, H. F. J., Peat, H. J., Bingham, R. G., and Vaughan, D. G.: British Antarctic Survey's aerogeophysical data: releasing 25 years of airborne gravity, magnetic, and radar datasets over Antarctica, Earth Syst. Sci. Data, 14, 3379–3410, https://doi.org/10.5194/essd-14-3379-2022, 2022. a, b

Fudge, T. J., Hills, B. H., Horlings, A. N., Holschuh, N., Christian, J. E., Davidge, L., Hoffman, A., O'Connor, G. K., Christianson, K., and Steig, E. J.: A site for deep ice coring at West Hercules Dome: results from ground-based geophysics and modeling, J. Glaciol., 69, 538–550, https://doi.org/10.1017/jog.2022.80, 2023. a

Fujita, S., Maeno, H., Uratsuka, S., Furukawa, T., Mae, S., Fujii, Y., and Watanabe, O.: Nature of radio echo layering in the Antarctic Ice Sheet detected by a two‐frequency experiment, J. Geophys. Res.-Sol. Ea., 104, 13013–13024, https://doi.org/10.1029/1999jb900034, 1999. a

Fujita, S., Holmlund, P., Andersson, I., Brown, I., Enomoto, H., Fujii, Y., Fujita, K., Fukui, K., Furukawa, T., Hansson, M., Hara, K., Hoshina, Y., Igarashi, M., Iizuka, Y., Imura, S., Ingvander, S., Karlin, T., Motoyama, H., Nakazawa, F., Oerter, H., Sjöberg, L. E., Sugiyama, S., Surdyk, S., Ström, J., Uemura, R., and Wilhelms, F.: Spatial and temporal variability of snow accumulation rate on the East Antarctic ice divide between Dome Fuji and EPICA DML, The Cryosphere, 5, 1057–1081, https://doi.org/10.5194/tc-5-1057-2011, 2011. a

Giannopoulos, A.: Modelling ground penetrating radar by GprMax, Construct. Build. Mater., 19, 755–762, https://doi.org/10.1016/j.conbuildmat.2005.06.007, 2005. a

Giannopoulos, A.: Unsplit Implementation of Higher Order PMLs, IEEE T. Antennas Prop., 60, 1479–1485, https://doi.org/10.1109/TAP.2011.2180344, 2012. a

Glen, J. W. and Paren, J. G.: The Electrical Properties of Snow and Ice, J. Glaciol., 15, 15–38, https://doi.org/10.3189/s0022143000034249, 1975. a

Hale, R., Miller, H., Gogineni, S., Yan, J. B., Rodriguez-Morales, F., Leuschen, C., Paden, J., Li, J., Binder, T., Steinhage, D., Gehrmann, M., and Braaten, D.: Multi-Channel Ultra-Wideband Radar Sounder and Imager, 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), 2112–2115, https://doi.org/10.1109/igarss.2016.7729545, 2016. a

Hammer, C. U.: Acidity of Polar Ice Cores in Relation to Absolute Dating, Past Volcanism, and Radio–Echoes, J. Glaciol., 25, 359–372, https://doi.org/10.3189/s0022143000015227, 1980. a

Harrison, C. H.: Radio Echo Sounding of Horizontal Layers in Ice, J. Glaciol., 12, 383–397, https://doi.org/10.3189/s0022143000031804, 1973. a

Holschuh, N., Christianson, K., and Anandakrishnan, S.: Power loss in dipping internal reflectors, imaged using ice-penetrating radar, Ann. Glaciol., 55, 49–56, https://doi.org/10.3189/2014aog67a005, 2014. a

Howat, I. M., Porter, C., Smith, B. E., Noh, M.-J., and Morin, P.: The Reference Elevation Model of Antarctica, The Cryosphere, 13, 665–674, https://doi.org/10.5194/tc-13-665-2019, 2019. a, b, c, d, e

Howat, I., Porter, C., Noh, M., Husby, E., Khuvis, S., Danish, E., Tomko, K., Gardiner, J., Negrete, A., Yadav, B., Klassen, J., Kelleher, C., Cloutier, M., Bakker, J., Enos, J., Arnold, G., Bauer, G., and Morin, P.: The Reference Elevation Model of Antarctica – Mosaics, Version 2, Harvard Dataverse, V1 [data set], https://doi.org/10.7910/DVN/EBW8UC, 2022. a

Huybrechts, P., Rybak, O., Steinhage, D., and Pattyn, F.: Past and present accumulation rate reconstruction along the Dome Fuji–Kohnen radio-echo sounding profile, Dronning Maud Land, East Antarctica, Ann. Glaciol., 50, 112–120, https://doi.org/10.3189/172756409789097513, 2009. a

IPCC: Ocean, Cryosphere and Sea Level Change, Cambridge University Press, 1211–1362, https://doi.org/10.1017/9781009157896.011, 2023. a

Jacobel, R. W. and Welch, B. C.: A time marker at 17.5 kyr BP detected throughout West Antarctica, Ann. Glaciol., 41, 47–51, https://doi.org/10.3189/172756405781813348, 2005. a

Jansen, D., Franke, S., Bauer, C. C., Binder, T., Dahl-Jensen, D., Eichler, J., Eisen, O., Hu, Y., Kerch, J., Llorens, M.-G., Miller, H., Neckel, N., Paden, J., Riese, T. d., Sachau, T., Stoll, N., Weikusat, I., Wilhelms, F., Zhang, Y., and Bons, P. D.: Shear margins in upper half of Northeast Greenland Ice Stream were established two millennia ago, Nat. Commun., 15, 1193, https://doi.org/10.1038/s41467-024-45021-8, 2024. a, b, c, d

Karlsson, N. B., Dahl-Jensen, D., Gogineni, S. P., and Paden, J. D.: Tracing the depth of the Holocene ice in North Greenland from radio-echo sounding data, Ann. Glaciol., 54, 44–50, https://doi.org/10.3189/2013aog64a057, 2013. a

Karlsson, N. B., Binder, T., Eagles, G., Helm, V., Pattyn, F., Van Liefferinge, B., and Eisen, O.: Glaciological characteristics in the Dome Fuji region and new assessment for “Oldest Ice”, The Cryosphere, 12, 2413–2424, https://doi.org/10.5194/tc-12-2413-2018, 2018. a

Kjær, K. H., Larsen, N. K., Binder, T., Bjørk, A. A., Eisen, O., Fahnestock, M. A., Funder, S., Garde, A. A., Haack, H., Helm, V., Houmark-Nielsen, M., Kjeldsen, K. K., Khan, S. A., Machguth, H., McDonald, I., Morlighem, M., Mouginot, J., Paden, J. D., Waight, T. E., Weikusat, C., Willerslev, E., and MacGregor, J. A.: A large impact crater beneath Hiawatha Glacier in northwest Greenland, Sci. Adv., 4, eaar8173, https://doi.org/10.1126/sciadv.aar8173, 2018. a

Koch, I., Drews, R., Franke, S., Jansen, D., Oraschewski, F. M., Muhle, L. S., Višnjević, V., Matsuoka, K., Pattyn, F., and Eisen, O.: Radar internal reflection horizons from multisystem data reflect ice dynamic and surface accumulation history along the Princess Ragnhild Coast, Dronning Maud Land, East Antarctica, J. Glaciol., 70, e18, https://doi.org/10.1017/jog.2023.93, 2023. a, b, c

Lenaerts, J. T. M., Meijgaard, E. v., Broeke, M. R. v. d., Ligtenberg, S. R. M., Horwath, M., and Isaksson, E.: Recent snowfall anomalies in Dronning Maud Land, East Antarctica, in a historical and future climate perspective, Geophys. Res. Lett., 40, 2684–2688, https://doi.org/10.1002/grl.50559, 2013. a, b, c

Leysinger Vieli, G. J. C., Hindmarsh, R. C. A., Siegert, M. J., and Bo, S.: Time‐dependence of the spatial pattern of accumulation rate in East Antarctica deduced from isochronic radar layers using a 3‐D numerical ice flow model, J. Geophys. Res.-Earth Surf., 116, F02028, https://doi.org/10.1029/2010jf001785, 2011. a, b, c, d, e, f

Leysinger Vieli, G. J.-M. C., Siegert, M. J., and Payne, A. J.: Reconstructing ice-sheet accumulation rates at ridge B, East Antarctica, Ann. Glaciol., 39, 326–330, https://doi.org/10.3189/172756404781814519, 2004. a, b

Leysinger Vieli, G. J.-M. C., Martín, C., Hindmarsh, R. C. A., and Lüthi, M. P.: Basal freeze-on generates complex ice-sheet stratigraphy, Nat. Commun., 9, 4669, https://doi.org/10.1038/s41467-018-07083-3, 2018. a, b, c

Lilien, D. A., Steinhage, D., Taylor, D., Parrenin, F., Ritz, C., Mulvaney, R., Martín, C., Yan, J.-B., O'Neill, C., Frezzotti, M., Miller, H., Gogineni, P., Dahl-Jensen, D., and Eisen, O.: Brief communication: New radar constraints support presence of ice older than 1.5 Myr at Little Dome C, The Cryosphere, 15, 1881–1888, https://doi.org/10.5194/tc-15-1881-2021, 2021. a

Lin, J., Svensson, A., Hvidberg, C. S., Lohmann, J., Kristiansen, S., Dahl-Jensen, D., Steffensen, J. P., Rasmussen, S. O., Cook, E., Kjær, H. A., Vinther, B. M., Fischer, H., Stocker, T., Sigl, M., Bigler, M., Severi, M., Traversi, R., and Mulvaney, R.: Magnitude, frequency and climate forcing of global volcanism during the last glacial period as seen in Greenland and Antarctic ice cores (60–9 ka), Clim. Past, 18, 485–506, https://doi.org/10.5194/cp-18-485-2022, 2022. a, b, c, d, e, f, g, h, i, j, k, l

Lin, J., Abbott, P. M., Sigl, M., Steffensen, J. P., Mulvaney, R., Severi, M., and Svensson, A.: Bipolar ice-core records constrain possible dates and global radiative forcing following the ∼74 ka Toba eruption, Quaternary Sci. Rev., 312, 108 162, https://doi.org/10.1016/j.quascirev.2023.108162, 2023. a

Lythe, M. B. and Vaughan, D. G.: BEDMAP: A new ice thickness and subglacial topographic model of Antarctica, J. Geophys. Res.-Sol. Ea., 106, 11335–11351, https://doi.org/10.1029/2000jb900449, 2001. a

MacGregor, J. A., Fahnestock, M. A., Catania, G. A., Paden, J. D., Gogineni, S. P., Young, S. K., Rybarski, S. C., Mabrey, A. N., Wagman, B. M., and Morlighem, M.: Radiostratigraphy and age structure of the Greenland Ice Sheet, J. Geophys. Res.-Earth Surf., 120, 212–241, https://doi.org/10.1002/2014jf003215, 2015. a, b, c

MacGregor, J. A., Fahnestock, M. A., Paden, J. D., Li, J., Harbeck, J. P., and Aschwanden, A.: A revised and expanded deep radiostratigraphy of the Greenland Ice Sheet from airborne radar sounding surveys between 1993–2019, Earth Syst. Sci. Data Discuss. [preprint], https://doi.org/10.5194/essd-2024-578, in review, 2025. a

Millar, D. H. M.: Radio-echo layering in polar ice sheets and past volcanic activity, Nature, 292, 441–443, https://doi.org/10.1038/292441a0, 1981. a, b, c

Mojtabavi, S., Eisen, O., Franke, S., Jansen, D., Steinhage, D., Paden, J., Dahl-Jensen, D., Weikusat, I., Eichler, J., and Wilhelms, F.: Origin of englacial stratigraphy at three deep ice core sites of the Greenland Ice Sheet by synthetic radar modelling, J. Glaciol., 68, 799–811, https://doi.org/10.1017/jog.2021.137, 2022. a

Moore, J. C.: High-resolution dielectric profiling of ice cores, J. Glaciol., 39, 245–248, https://doi.org/10.3189/s0022143000015902, 1993. a

Moqadam, H. and Eisen, O.: Review article: Feature tracing in radio-echo sounding products of terrestrial ice sheets and planetary bodies, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2024-1674, 2024. a

Morlighem, M.: MEaSUREs BedMachine Antarctica. (NSIDC-0756, Version 3), Boulder, Colorado USA, NASA National Snow and Ice Data Center Distributed Active Archive Center [data Set], https://doi.org/10.5067/FPSU0V1MWUB6, 2022. a

Morlighem, M., Rignot, E., Binder, T., Blankenship, D., Drews, R., Eagles, G., Eisen, O., Ferraccioli, F., Forsberg, R., Fretwell, P., Goel, V., Greenbaum, J. S., Gudmundsson, H., Guo, J., Helm, V., Hofstede, C., Howat, I., Humbert, A., Jokat, W., Karlsson, N. B., Lee, W. S., Matsuoka, K., Millan, R., Mouginot, J., Paden, J., Pattyn, F., Roberts, J., Rosier, S., Ruppel, A., Seroussi, H., Smith, E. C., Steinhage, D., Sun, B., Broeke, M. R. v. d., Ommen, T. D. v., Wessem, M. v., and Young, D. A.: Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet, Nat. Geosci., 13, 132–137, https://doi.org/10.1038/s41561-019-0510-8, 2020. a, b, c, d, e, f, g, h

Mouginot, J., Rignot, E., and Scheuchl, B.: Continent‐Wide, Interferometric SAR Phase, Mapping of Antarctic Ice Velocity, Geophys. Res. Lett., 46, 9710–9718, https://doi.org/10.1029/2019gl083826, 2019a. a, b, c, d

Mouginot, J., Rignot, E., Scheuchl, B.: Phase-Based Antarctica Ice Velocity Map, Dryad [data set], https://doi.org/10.7280/D10D4Z, 2019b. a

Muldoon, G., Blankenship, D. D., Jackson, C., and Young, D. A.: AGASEA 4.7 ka Englacial Isochron over the Thwaites Glacier Catchment, USAP-DC [data set], https://doi.org/10.15784/601673, 2023. a

Muldoon, G. R., Jackson, C. S., Young, D. A., and Blankenship, D. D.: Bayesian estimation of englacial radar chronology in Central West Antarctica, Dynamics and Statistics of the Climate System, 3, dzy004, https://doi.org/10.1093/climatesystem/dzy004, 2018. a, b, c, d, e, f

Mulvaney, R., Frezzotti, M., Ritz, C., Martín, C., Rix, J., and Chung, A.: Internal reflecting horizons in the Dome C region of Antarctica, PANGAEA [data set], https://doi.org/10.1594/PANGAEA.963470, 2023. a, b, c

Nixdorf, U., Steinhage, D., Meyer, U., Hempel, L., Jenett, M., Wachs, P., and Miller, H.: The newly developed airborne radio-echo sounding system of the AWI as a glaciological tool, Ann. Glaciol., 29, 231–238, https://doi.org/10.3189/172756499781821346, 1999. a, b

Näslund, J.: Subglacial Preservation of Valley Morphology at Amundsenisen, Western Dronning Maud Land, Antarctica, Earth Surf. Process. Landf., 22, 441–455, https://doi.org/10.1002/(sici)1096-9837(199705)22:5<441::aid-esp696>3.0.co;2-4, 1997. a

Open Polar Radar: Open Polar Radar. opr (Version 3.0.1), Zenodo [software], https://doi.org/10.5281/zenodo.5683959, 2023. a, b

Oyabu, I., Kawamura, K., Fujita, S., Inoue, R., Motoyama, H., Fukui, K., Hirabayashi, M., Hoshina, Y., Kurita, N., Nakazawa, F., Ohno, H., Sugiura, K., Suzuki, T., Tsutaki, S., Abe-Ouchi, A., Niwano, M., Parrenin, F., Saito, F., and Yoshimori, M.: Temporal variations of surface mass balance over the last 5000 years around Dome Fuji, Dronning Maud Land, East Antarctica, Clim. Past, 19, 293–321, https://doi.org/10.5194/cp-19-293-2023, 2023. a

Panton, C. and Karlsson, N. B.: Automated mapping of near bed radio-echo layer disruptions in the Greenland Ice Sheet, Earth Planet. Sc. Lett., 432, 323–331, https://doi.org/10.1016/j.epsl.2015.10.024, 2015. a

Paren, J. G. and Robin, G. D. Q.: Internal Reflections in Polar Ice Sheets, J. Glaciol., 14, 251–259, https://doi.org/10.3189/s0022143000021730, 1975. a

Pritchard, H. D., Fretwell, P. T., Fremand, A. C., Bodart, J. A., Kirkham, J. D., Aitken, A., Bamber, J., Bell, R., Bianchi, C., Bingham, R. G., Blankenship, D. D., Casassa, G., Christianson, K., Conway, H., Corr, H. F. J., Cui, X., Damaske, D., Damm, V., Dorschel, B., Drews, R., Eagles, G., Eisen, O., Eisermann, H., Ferraccioli, F., Field, E., Forsberg, R., Franke, S., Goel, V., Gogineni, S. P., Greenbaum, J., Hills, B., Hindmarsh, R. C. A., Hoffman, A. O., Holschuh, N., Holt, J. W., Humbert, A., Jacobel, R. W., Jansen, D., Jenkins, A., Jokat, W., Jong, L., Jordan, T. A., King, E. C., Kohler, J., Krabill, W., Maton, J., Gillespie, M. K., Langley, K., Lee, J., Leitchenkov, G., Leuschen, C., Luyendyk, B., MacGregor, J. A., MacKie, E., Moholdt, G., Matsuoka, K., Morlighem, M., Mouginot, J., Nitsche, F. O., Nost, O. A., Paden, J., Pattyn, F., Popov, S., Rignot, E., Rippin, D. M., Rivera, A., Roberts, J. L., Ross, N., Ruppel, A., Schroeder, D. M., Siegert, M. J., Smith, A. M., Steinhage, D., Studinger, M., Sun, B., Tabacco, I., Tinto, K. J., Urbini, S., Vaughan, D. G., Wilson, D. S., Young, D. A., and Zirizzotti, A.: Bedmap3 updated ice bed, surface and thickness gridded datasets for Antarctica, Sci. Data, 12, 414, https://doi.org/10.1038/s41597-025-04672-y, 2025. a

Riedel, S., Jokat, W., and Steinhage, D.: Mapping tectonic provinces with airborne gravity and radar data in Dronning Maud Land, East Antarctica, Geophys. J. Int., 189, 414–427, https://doi.org/10.1111/j.1365-246x.2012.05363.x, 2012. a, b

Rignot, E., Mouginot, J., Scheuchl, B., Broeke, M. V. D., Wessem, M. J. v., and Morlighem, M.: Four decades of Antarctic Ice Sheet mass balance from 1979–2017, P. Natl. Acad. Sci. USA, 116, 1095–1103, https://doi.org/10.1073/pnas.1812883116, 2019 (data available at: http://imbie.org/imbie-3/drainage-basins/, last access: 17 January 2025). a, b, c, d

Rippin, D. M., Bamber, J. L., Siegert, M. J., Vaughan, D. G., and Corr, H. F. J.: Basal topography and ice flow in the Bailey/Slessor region of East Antarctica, J. Geophys. Res.-Earth Surf., 108, 6008, https://doi.org/10.1029/2003jf000039, 2003. a

Rippin, D. M., Siegert, M. J., Bamber, J. L., Vaughan, D. G., and Corr, H. F. J.: Switch‐off of a major enhanced ice flow unit in East Antarctica, Geophys. Res. Lett., 33, L15501, https://doi.org/10.1029/2006gl026648, 2006. a, b

Robin, G. D. Q., Evans, S., and Bailey, J. T.: Interpretation of radio echo sounding in polar ice sheets, Philos. T. Roy. Soc. Lond. A, 265, 437–505, https://doi.org/10.1098/rsta.1969.0063, 1969. a

Robin, G. D. Q., Drewry, D. J., and Meldrum, D. T.: International studies of ice sheet and bedrock, Philos. T. Roy. Soc. Lond. B, 279, 185–196, https://doi.org/10.1098/rstb.1977.0081, 1977. a

Rodriguez-Morales, F., Byers, K., Crowe, R., Player, K., Hale, R. D., Arnold, E. J., Smith, L., Gifford, C. M., Braaten, D., Panton, C., Gogineni, S., Leuschen, C. J., Paden, J. D., Li, J., Lewis, C. C., Panzer, B., Alvestegui, D. G.-G., and Patel, A.: Advanced Multifrequency Radar Instrumentation for Polar Research, IEEE T. Geosci. Remote, 52, 2824–2842, https://doi.org/10.1109/tgrs.2013.2266415, 2013. a

Rose, K. C., Ross, N., Jordan, T. A., Bingham, R. G., Corr, H. F. J., Ferraccioli, F., Le Brocq, A. M., Rippin, D. M., and Siegert, M. J.: Ancient pre-glacial erosion surfaces preserved beneath the West Antarctic Ice Sheet, Earth Surf. Dynam., 3, 139–152, https://doi.org/10.5194/esurf-3-139-2015, 2015. a

Rotschky, G., Holmlund, P., Isaksson, E., Mulvaney, R., Oerter, H., Broeke, M. R. V. D., and Winther, J.-G.: A new surface accumulation map for western Dronning Maud Land, Antarctica, from interpolation of point measurements, J. Glaciol., 53, 385–398, https://doi.org/10.3189/002214307783258459, 2007. a, b, c

Sanderson, R., Ross, N., Winter, K., Bingham, R., Callard, L., Jordan, T., and Young, D.: Dated radar stratigraphy between Dome A and South Pole, East Antarctica, derived from AGAP North PASIN (2008-2009) and PolarGAP PASIN2 (2015–2016) surveys (Version 1.0), NERC EDS UK Polar Data Centre [data set], https://doi.org/10.5285/cfafb639-991a-422f-9caa-7793c195d316, 2023a. a, b

Sanderson, R. J., Winter, K., Callard, S. L., Napoleoni, F., Ross, N., Jordan, T. A., and Bingham, R. G.: Englacial architecture of Lambert Glacier, East Antarctica, The Cryosphere, 17, 4853–4871, https://doi.org/10.5194/tc-17-4853-2023, 2023b. a

Sanderson, R. J., Ross, N., Winter, K., Bingham, R. G., Callard, S. L., Jordan, T. A., and Young, D. A.: Dated radar-stratigraphy between Dome A and South Pole, East Antarctica: old ice potential and ice sheet history, J. Glaciol., 70, e74, https://doi.org/10.1017/jog.2024.60, 2024. a, b, c, d, e

Santin, I., Roncoroni, G., Forte, E., Gutgesell, P., and Pipan, M.: GPR modelling and inversion to quantify the debris content within ice, Near Surf. Geophys., 22, 220–234, https://doi.org/10.1002/nsg.12274, 2023. a

Siegert, M. J. and Payne, A. J.: Past rates of accumulation in central West Antarctica, Geophys. Res. Lett., 31, L12403, https://doi.org/10.1029/2004gl020290, 2004. a

Siegert, M. J., Welch, B., Morse, D., Vieli, A., Blankenship, D. D., Joughin, I., King, E. C., Vieli, G. J.-M. C. L., Payne, A. J., and Jacobel, R.: Ice Flow Direction Change in Interior West Antarctica, Science, 305, 1948–1951, https://doi.org/10.1126/science.1101072, 2004. a, b

Siegert, M. J., Pokar, M., Dowdeswell, J. A., and Benham, T.: Radio‐echo layering in West Antarctica: a spreadsheet dataset, Earth Surf. Process. Landf., 30, 1583–1591, https://doi.org/10.1002/esp.1238, 2005. a

Steinhage, D.: Beitrage aus geophysikalischen Messungen in Dronning Maud Land, Antarktis, zur Auffindung eines optimalen Bohrpunktes für eine Eiskerntiefbohrung, Berichte zur Polar und Meeresforschung, 2001. a, b, c, d

Steinhage, D., Nixdorf, U., Meyer, U., and Miller, H.: Subglacial topography and internal structure of central and western Dronning Maud Land, Antarctica, determined from airborne radio echo sounding, J. Appl. Geophys., 47, 183–189, https://doi.org/10.1016/s0926-9851(01)00063-5, 2001. a, b, c, d, e

Steinhage, D., Kipfstuhl, S., Nixdorf, U., and Miller, H.: Internal structure of the ice sheet between Kohnen station and Dome Fuji, Antarctica, revealed by airborne radio-echo sounding, Ann. Glaciol., 54, 163–167, https://doi.org/10.3189/2013aog64a113, 2013. a, b, c, d, e, f, g, h, i, j, k, l

Sutter, J., Fischer, H., and Eisen, O.: Investigating the internal structure of the Antarctic ice sheet: the utility of isochrones for spatiotemporal ice-sheet model calibration, The Cryosphere, 15, 3839–3860, https://doi.org/10.5194/tc-15-3839-2021, 2021. a, b, c, d

Svensson, A., Bigler, M., Blunier, T., Clausen, H. B., Dahl-Jensen, D., Fischer, H., Fujita, S., Goto-Azuma, K., Johnsen, S. J., Kawamura, K., Kipfstuhl, S., Kohno, M., Parrenin, F., Popp, T., Rasmussen, S. O., Schwander, J., Seierstad, I., Severi, M., Steffensen, J. P., Udisti, R., Uemura, R., Vallelonga, P., Vinther, B. M., Wegner, A., Wilhelms, F., and Winstrup, M.: Direct linking of Greenland and Antarctic ice cores at the Toba eruption (74 ka BP), Clim. Past, 9, 749–766, https://doi.org/10.5194/cp-9-749-2013, 2013. a, b, c, d, e, f

Taflove, A. and Hagness, S. C.: Computational electrodynamics: the finite-difference time-domain method, in: The Electrical Engineering Handbook, 3rd edn., edited by: Chen, W.-K., Artech House, ISBN 0-12-170960-4, 2005. a

Van Dalum, C., Van de Berg, W. J., and Van den Broeke, M.: Monthly RACMO2.4p1 data for Greenland (11 km) and Antarctica (27 km) for SMB, SEB, near-surface temperature and wind speed (2006–2015), Zenodo [data set], https://doi.org/10.5281/zenodo.10854319, 2024. a

Van Liefferinge, B., Taylor, D., Tsutaki, S., Fujita, S., Gogineni, P., Kawamura, K., Matsuoka, K., Moholdt, G., Oyabu, I., Abe-Ouchi, A., Awasthi, A., Buizert, C., Gallet, J.-C., Isaksson, E., Motoyama, H., Nakazawa, F., Ohno, H., O’Neill, C., Pattyn, F., and Sugiura, K.: Surface Mass Balance Controlled by Local Surface Slope in Inland Antarctica: Implications for Ice-Sheet Mass Balance and Oldest Ice Delineation in Dome Fuji, Geophys. Res. Lett., 48, e2021GL094966, https://doi.org/10.1029/2021GL094966, 2021.  a

Višnjević, V., Drews, R., Schannwell, C., Koch, I., Franke, S., Jansen, D., and Eisen, O.: Predicting the steady-state isochronal stratigraphy of ice shelves using observations and modeling, The Cryosphere, 16, 4763–4777, https://doi.org/10.5194/tc-16-4763-2022, 2022. a, b, c

Wang, Z., Chung, A., Steinhage, D., Parrenin, F., Freitag, J., and Eisen, O.: Mapping age and basal conditions of ice in the Dome Fuji region, Antarctica, by combining radar internal layer stratigraphy and flow modeling, The Cryosphere, 17, 4297–4314, https://doi.org/10.5194/tc-17-4297-2023, 2023a. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u

Wang, Z., Chung, A., Steinhage, D., Parrenin, F., Freitag, J., and Eisen, O.: Radar internal layer stratigraphy in the Dome Fuji region, Antarctica, PANGAEA [data set], https://doi.org/10.1594/PANGAEA.958462, 2023b. a, b

Warren, C., Giannopoulos, A., and Giannakis, I.: gprMax: Open source software to simulate electromagnetic wave propagation for Ground Penetrating Radar, Comput. Phys. Commun., 209, 163–170, https://doi.org/10.1016/j.cpc.2016.08.020, 2016. a

Wilhelms, F., Kipfstuhl, J., Miller, H., Heinloth, K., and Firestone, J.: Precise dielectric profiling of ice cores: a new device with improved guarding and its theory, J. Glaciol., 44, 171–174, https://doi.org/10.3189/s002214300000246x, 1998. a

Winter, A., Steinhage, D., Arnold, E. J., Blankenship, D. D., Cavitte, M. G. P., Corr, H. F. J., Paden, J. D., Urbini, S., Young, D. A., and Eisen, O.: Comparison of measurements from different radio-echo sounding systems and synchronization with the ice core at Dome C, Antarctica, The Cryosphere, 11, 653–668, https://doi.org/10.5194/tc-11-653-2017, 2017. a, b, c, d

Winter, A., Steinhage, D., Creyts, T. T., and Eisen, O.: Radio-echo sounding isochrone depths in the East Antarctic Ice Sheet, PANGAEA [data set], https://doi.org/10.1594/PANGAEA.895528, 2018. a, b

Winter, A., Steinhage, D., Creyts, T. T., Kleiner, T., and Eisen, O.: Age stratigraphy in the East Antarctic Ice Sheet inferred from radio-echo sounding horizons, Earth Syst. Sci. Data, 11, 1069–1081, https://doi.org/10.5194/essd-11-1069-2019, 2019. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u

Yee, K.: Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media, IEEE T. Antenn. Propag., 14, 302–307, https://doi.org/10.1109/TAP.1966.1138693, 1966. a

Zhang, Y., Sachau, T., Franke, S., Yang, H., Li, D., Weikusat, I., and Bons, P. D.: Formation Mechanisms of Large‐Scale Folding in Greenland's Ice Sheet, Geophys. Res. Lett., 51, e2024GL109492, https://doi.org/10.1029/2024gl109492, 2024. a

Download
Short summary
The study presents internal reflection horizons (IRHs) over an area of 450 000 km² from western...
Share