An ancient river landscape preserved beneath the East Antarctic Ice Sheet | Nature Communications

• Explore content
• About the journal
• Publish with us
• Sign up for alerts
• RSS feed

1. nature
2. nature communications
3. articles
4. article

• Article
• Open access
• Stewart S. R. Jamieson
• Neil Ross
• Guy J. G. Paxman
• Fiona J. Clubb
• Duncan A. Young
• Shuai Yan
• Jamin Greenbaum
• Donald D. Blankenship &
• Martin J. Siegert

Nature Communications 14, Article number: 6507 (2023) Cite this article

49k Accesses

7 Citations

2014 Altmetric

Metrics

Abstract

The East Antarctic Ice Sheet (EAIS) has its origins ca. 34 million years ago. Since then, the impact of climate change and past fluctuations in the EAIS margin has been reflected in periods of extensive vs. restricted ice cover and the modification of much of the Antarctic landscape. Resolving processes of landscape evolution is therefore critical for establishing ice sheet history, but it is rare to find unmodified landscapes that record past ice conditions. Here, we discover an extensive relic pre-glacial landscape preserved beneath the central EAIS despite millions of years of ice cover. The landscape was formed by rivers prior to ice sheet build-up but later modified by local glaciation before being dissected by outlet glaciers at the margin of a restricted ice sheet. Preservation of the relic surfaces indicates an absence of significant warm-based ice throughout their history, suggesting any transitions between restricted and expanded ice were rapid.

Similar content being viewed by others

Mid-Pleistocene aridity and landscape shifts promoted Palearctic hominin dispersals

Article Open access 27 November 2024

Lateral river erosion impacts the preservation of Neolithic enclosures in alluvial plains

Article Open access 02 October 2023

Late Cenozoic climate change paces landscape adjustments to Yukon River capture

Article 20 July 2020

Introduction

The glaciation of Antarctica was triggered by global climatic cooling over the Cenozoic Era (Fig. 1). During the Eocene, glaciation was likely restricted to ephemeral ice masses and small-scale mountain glaciers in regions of high topography12. However, a step-change in ice extent and volume occurred at the Eocene-Oligocene transition (EOT; ca. 34.0–33.5 Ma; Fig. 1) when the first widespread Antarctic glaciation was recorded in marine sediment records34. This transition to a glaciated Antarctica was potentially caused by a combination of CO2 dropping below a key threshold5, associated feedbacks within the carbon cycle67, and the opening/deepening of circum-Antarctic ocean gateways8. In East Antarctica, the ice sheet likely nucleated on the high topography of the Gamburtsev Subglacial Mountains, Transantarctic Mountains, and Dronning Maud Land259 (Fig. 1). Expansion and coalescence of independent ice masses on these highlands led to the growth of the continental-scale EAIS510.

Fig. 1: Location and climate evolution.

a

Surface elevation of the Antarctic Ice Sheet from the Reference Elevation Model of Antarctica

68

. The black line marks the present-day grounding line

74

; pale areas are floating ice shelves.

b

Bed elevation of Antarctica beneath the grounded ice sheet from BedMachine v.3

64

. Surface and bed elevations are relative to the global mean sea level (the EIGEN-6C4 geoid). The black dashed box shows the extent of the maps in Figs.

2

and

5

, and lines of latitude are separated by 5°. ASB Aurora Subglacial Basin, DML Dronning Maud Land, GSM Gamburtsev Subglacial Mountains, TAM Transantarctic Mountains, WSB Wilkes Subglacial Basin.

c

Composite global oxygen isotope record from benthic foraminifera

75

. Lower (more negative)

δ

18

O values indicate lower ice volume and/or warmer ocean temperatures; higher (more positive)

δ

18

O values indicate greater ice volume and/or cooler ocean temperatures. Labels show the timings of major climate intervals or transitions (red), glacial expansions (black), and prevailing glacial conditions (blue). Plio. = Pliocene; Plt. = Pleistocene; note that the Holocene epoch (0.012–0 Ma) is not shown on the timescale.

Full size image

The EAIS has fluctuated on a range of spatial and temporal scales since the EOT, with warm-based glaciers reaching the coast and advancing and retreating at orbitally-paced timescales during the Oligocene1112, and major variations in ice extent and volume occurring during the early- to mid-Miocene13141516. Following the Mid-Miocene Climatic Optimum (ca. 17.0–14.8 Ma; Fig. 1), further atmospheric cooling led to the establishment of an arid polar climate, and by ca. 14 Ma, a continental-scale ice sheet became persistent, with evidence for widespread cold-based ice and advances and retreats of the ice-sheet margin to and from the continental shelf edge17. Offshore sedimentary records suggest the EAIS may have experienced retreat during subsequent periods of warmer climate, for example, the mid-Piacenzian warm period (mPWP; ca. 3.26–3.02 Ma)181920 and potentially the Pleistocene interglacials of MIS 11c (ca. 430–400 ka)21 and 5e (ca. 125 ka)222324. Numerical ice sheet models indicate that retreat may have been focussed in low-lying marine-based sectors of the EAIS, such as the Wilkes and Aurora Subglacial Basins (Fig. 1) during Plio-Pleistocene warm intervals, although the extent of this retreat remains uncertain252627.

Much of the evidence for past ice margin fluctuations and dynamic change of the EAIS is derived from global records of sea-level change and/or ice volume28 or from local records of glacier fluctuation on the exposed land around the coast and in the offshore sediments on the continental shelf12181929. Although we have been able to infer from these lines of evidence that significant past changes in the EAIS have occurred, there is limited direct evidence that records the degree to which the ice margin retreated inland of its present-day position during past warm periods. This means the extent and possible rates of past EAIS retreats are poorly resolved. However, as ice sheets fluctuate, they modify the landscape upon which they rest, leaving a fingerprint representing the long-term average glacial conditions to which it has been subjected30. Interpretation of subglacial topography, measured using radio-echo sounding (RES) surveying, can therefore be used to infer past processes of landscape evolution and former ice dynamics93132333435. Processes that may be reflected in the subglacial landscape include but are not limited to, glacial erosion, solid earth adjustments to erosional unloading, and fluvial incision under scenarios of more restricted EAIS extent36. Accordingly, interpreting the morphology of the East Antarctic subglacial landscape and its long-term evolution offers a significant and hitherto relatively little exploited opportunity for understanding the past extent, fluctuations, flow, and thermal regime of the EAIS.

Here we use geophysical data to analyse an area in East Antarctica with a subglacial landscape that appears to be at odds with the modern ice sheet. It lies adjacent to the low-lying Aurora and Schmidt Subglacial Basins and inland of the Denman and Totten glaciers; a sector of the EAIS known to be sensitive to past and potentially future climate and ocean warming (Fig. 2)1333. Young et al.37 identified palaeo fjord landscapes cutting the highlands that bound these basins (Fig. 2a), interpreting them as evidence of Oligocene and/or early Miocene incision (i.e., between ca. 34 and 14 Ma) under an ice sheet significantly smaller than today37. Our geophysical study focuses on the higher elevation land (now bed) surfaces between the fjords, which may be more ancient features preserved beneath the EAIS. Our overarching aim is to evaluate the evolution and age of these land surfaces and to assess their implications for EAIS history.

Fig. 2: Characteristics of the subglacial highlands and troughs of ‘Highland A’37 .

a

BedMachine version 3 bed elevation

64

, contours (black) in 250 m increments. Bed dataset is transparent and underlain by a hillshade of the same dataset. Extent of Fig.

3a

shown by red box;

b

BedMachine version 3 ice thickness

64

. Black lines are radar profiles, with red parts showing where bed was picked. Ice thickness contours (black) are in 125 m increments. Ice thickness dataset is transparent and underlain by a hillshade of the same dataset;

c

Ice surface elevation

68

, contours (black) in 50 m increments. Ice surface dataset is transparent and underlain by a hillshade of the same dataset;

d

RADARSAT image mosaic of Antarctic ice surface morphology

42

;

e

ice sheet velocity

41

;

f

Basal melting

47

identified as a percentage of 500 ensemble members recording basal melt. Backdrop imagery for (

e

) and (

f

) is a RADARSAT image mosaic of Antarctica

42

.

Full size image

To test the hypothesis of an ancient land surface preserved beneath the modern EAIS, we quantify landscape characteristics from RES data collected as part of the International Collaborative Exploration of the Cryosphere through Airborne Profiling (ICECAP) airborne geophysical survey of the Aurora Subglacial Basin and Totten Glacier between 2008 and 201133373839 (Fig. 2b). To support this hypothesis, we would expect a specific set of criteria to be met, which would be relevant anywhere on the continent. These are: (1) topographic features should be similar to each other but not consistent with continental-scale ice flow; (2) the elevations of the surfaces should be coherent when corrected for more recent glacial trough incision and ice loading; (3) the basal thermal regime must be cold-based to preserve the landscape. To evaluate these criteria, we: (1) map topographic features using satellite data that record changes in ice surface slope, which have previously been shown to reflect the large-scale subglacial topography3240 (Fig. 2d); (2) use along-track RES data (Fig. 2b) to ground-truth the geomorphology and quantify the subglacial topographic relief; (3) apply flexural modelling to understand whether the highland blocks are consistent with having once been a single land surface that was subsequently incised and uplifted by selective linear erosion37; (4) evaluate our findings in the context of modern and past ice flow and basal thermal conditions.

Results

We identify three distinct landscape blocks that constitute ‘Highland A’37. The blocks range between 121 and 173 km long and 73 and 85 km wide; each covers an area between ~7000 and 10,000 km2 (Fig. 3). The blocks are separated by two large-scale troughs, which are up to ~40 km wide and thought to have been generated by selective glacial erosion, hence previously referred to as ‘fjords’37. The southern trough is significantly deeper than the northern trough, with minimum elevations of 1480 m vs 620 m below global mean sea level, respectively. On a regional scale, the three highland landscape blocks appear to have similar morphological characteristics, each comprising rough terrain containing an intricate network of ridges and valleys (Fig. 3). If we were to assume these landscape blocks were once a coherent single land surface, it would have covered an area of more than 32,000 km2 (~1.5 x the area of Wales) with a long-axis of ~300 km. RES data confirm that across the three highland blocks, the buried topographic peaks have relatively uniform heights of 660–850 m above global mean sea level, with a subtle ‘domed’ pattern of tilting away from the southern and central blocks, which exhibit the highest elevations (Fig. 3).

Fig. 3: Mapping the subglacial landscape blocks in the Highland A region reveals a coherent valley network spanning three blocks.

a

Valley and ridge planforms (lines) identified in ice surface imagery compared with locations of topographic ‘highs’ (ridges) and ‘lows’ (valleys) identified from RES survey data (points). The background shows the RADARSAT ice surface imagery

42

. Dashed lines denote inferred former valley connections across south (S), central (C) and northern (N) topographic blocks. Dashed line A–A’ shows the location of the radargram in panel (

b

). The map is projected in EPSG:3031 and rotated by 25 degrees and its location is shown in Fig.

2a

.

b

Depth corrected 2D focused radargram from the January 2008 ICECAP campaign confirms surface features correspond to peaks and valleys. Depths relative to WGS84 ellipsoid.

c

Kernel density estimation plot of the ridge-to-valley spacing based on mapping in A. Vertical dashed lines denote means. Valley widths would be double the ridge-valley spacing. Source data are provided as a Source Data File.

Full size image

The first criterion, that the landscape fragments should be similar to each other but should not be consistent with modern continental-scale ice flow, was assessed by quantifying the valley morphology and comparing it to modern ice velocity data (Fig. 2e). Currently, modern ice flows broadly northwards towards the Wilkes Land coast and is very slow (<5 m yr−1) over the blocks and slightly faster (~25 m yr−1) as it flows over the deeper basins surrounding our study area41. RADARSAT ice surface imagery42, calibrated with airborne RES data3839, shows that the landscape blocks host a complex network of trunk and tributary valleys separated by sinuous mountain ridges, with valleys ranging between approx. 3 and 75 km in length (Fig. 3a).

Spatial analysis shows that the valley networks on all three highland blocks have a similar structure, with two dominant, near-orthogonal trends in valley orientation (approximately E–W (~70–90° / ~250–270°) and approximately N–S (~330–20° / ~150–200°); Fig. 4). The branching, dendritic structure of valleys (Fig. 3a) indicates that the landscape may have a fluvial origin. Purely fluvial systems tend to be characterised by relatively intricate valley networks, with high stream orders, numerous tributaries, and acute tributary junction angles, whereas glacial valley networks generally have lower stream orders, fewer contributing channels, and near-orthogonal junction angles43.

Fig. 4: Highland A blocks have similar hypsometries.

a

Combined hypsometry for all three highland blocks under in situ (i.e., ice loaded; yellow) and isostatically rebounded (i.e., ice-free; red) conditions using a laterally heterogeneous Te model

45

. Lower panels show rebounded hypsometry for the individual landscape blocks;

b

south block,

c

central block,

d

north block. Blocks are delineated by the outlines shown in Fig.

3a

. All elevations were obtained from the airborne RES data crossing the landscape blocks and are relative to present-day global mean sea level. N is the number of RES-derived bed elevation measurements within each block. Insets show polar histograms (sector diagrams) for the orientations of the valleys in each highland block. Valleys were assigned a geographic orientation but not directionality, hence the sector diagrams are symmetrical. The sector diagrams have been rotated to align with the approximate orientation of the map view in Fig.

3a

. Source data are provided as a Source Data file.

Full size image

Figure 3a shows evidence for both of these types of networks, with some valleys containing complex branching structures with multiple tributaries that are suggestive of fluvial incision, while others have a simpler network structure of tributaries entering the main stem at a relatively high junction angle, more indicative of glacial incision. The valley cross profiles tend to be U-shaped, separated by sharp ridges (Fig. 3b), indicative of erosive modification by valley glaciers. It is unlikely that these valleys were incised by subglacial meltwater, as in such cases, networks tend to be non-dendritic and oriented parallel to the overall direction of ice flow44.

The RADARSAT-derived valley networks also provide evidence for the potential coherency of the three blocks, as a number of N–S-trending valleys suggest a previous connection between the blocks (Fig. 3a). These connections have now been erased by the insertion of the troughs, so the valleys must predate the trough incision. There are also similar valley widths across each block, with median ridge-to-valley spacings of 2.1–2.7 km (Fig. 3c). Analysis of RES data (Fig. 3b) demonstrates that the mean relief (difference in elevation between ridge crests and valley floors) of the buried topography is approx. 800 m (South: 823 m; Central: 814 m; North: 757 m) with maximum relief of 1197 m (1021, 1197, and 1058 m). On average, the valley floors currently occupy similar elevations across the three blocks, with mean thalweg elevations of 238 m below sea level (South), 125 m below sea level (Central), and 206 m below sea level (North). However, valley floors are located over a range of elevations, with standard deviations of 427 m (South), 209 m (Central), and 241 m (North). Drainage density, a measure of valley length per unit area, is similar across all three blocks (0.13, 0.09, and 0.09 km per km2). The approximate uniformity of these metrics across the three blocks is indicative of a coherent single land surface.

Our analysis indicates that valley structure therefore meets criterion #1 due to its similarity across the three blocks and because it is not consistent with modern continental-scale ice flow. The dendritic valley structure, potential connectivity of valleys across the blocks, valley-ridge relief/morphology, and valley cross-sectional profiles are all consistent with a fluvial landscape subsequently modified by local-scale glaciation.

The second criterion, that the elevation of the ancient surface should be coherent when corrected for more recent glacial trough incision and ice loading, was assessed using flexural modelling. Our calculations show that the displacement induced by erosional unloading within the large-scale troughs and smaller-scale valleys has caused land surface uplift of up to 500 m (Fig. 5). This value is likely to be an upper limit because the two large troughs may have been partially developed prior to glaciation, for example, due to inherited tectonic structure and/or fluvial incision103335. The spatial pattern of uplift is dependent on the chosen value of the effective elastic thickness of the lithosphere (Te; Fig. 5). For low Te values (e.g., 5–10 km), the flexure exhibits short wavelength warping close to the margins of the large troughs that is not consistent with the observed elevations of the intervening highlands. By contrast, the flexure computed for higher Te values of 20–40 km predicts a pattern of uplift characterised by gentle domal tilting away from the deep trough between the southern and central blocks, within which the largest eroded thicknesses are predicted (Fig. 5). This pattern of tilting is more consistent with the observed peak elevations. Higher values of Te up to 80 km or more45 cannot be precluded, given that the sensitivity of the flexure to Te at the wavelengths of interest (up to 300 km) drops off at high Te values46.

Fig. 5: Flexural modelling of erosion-driven uplift indicates the coherency of the three highland blocks.

a

BedMachine v.3 bed topography

64

isostatically adjusted for the removal of the modern ice sheet load using a heterogeneous

Te

model

45 ,70

. White and blue circles mark topographic peaks used to estimate the distribution of valley and trough erosion.

b

Estimated erosion distribution, calculated by differencing the ice-free topography from a smooth accordant surface interpolated between the peaks shown in panel (

a

).

c

Erosion estimate along profile A–A’. The thickness of eroded material is shaded in yellow. Note that dashed lines indicate regions of the bed that are not imaged by RES data.

d

Same as panel (

c

) but for profile B–B’. Profile locations are shown in panels (

a

) and (

b

).

e

Computed flexural response to erosional unloading along profile A–A’ for five

Te

values. To aid visual comparison between the calculated flexure patterns and the tilt of the highland surfaces, the flexure has been vertically shifted upwards by 1200 m. The white line shows in situ bed elevation beneath the modern East Antarctic Ice Sheet.

f

Same as panel (

e

) but for profile B–B’. All elevations are relative to present-day global mean sea level.

Full size image

The flexural modelling results are consistent with the interpretation that a coherent highland landscape existed prior to the incision of the two deep troughs/fjords and that the three blocks were separated, uplifted, and subtly upwarped/domed as a result of this incision (thereby meeting criterion 2). The computed maximum amplitude of the erosion-driven flexure is ~500 m, whereas the elevations of the highland peaks are up to 1500–1700 m above sea level under ice-free conditions (Figs. 4 and 5). This implies that the highland peaks were situated at least 1–1.2 km above sea level prior to glacial incision, assuming a lack of vertical displacement driven by dynamic topography processes, which remain unconstrained in this region. The similarity of the hypsometries of the three blocks when isostatically corrected for ice loading, with consistently located modal peaks for the valley floors and ridge tops (Fig. 4), also offers further support for criterion 1.

The third criterion, that the basal thermal regime must have been cold-based to preserve the landscape, was tested by comparing our landscape to outcomes from previous numerical models of basal melting and with modern ice flow patterns. Modern ice flow is slow (<5 m yr−1), flows in a northwards direction across the blocks and does not appear to be topographically steered by the local valleys in our mapped landscape (Fig. 2e). Over the wider region, where the ice thickness is greater (~4000 m) (Fig. 2b), elevated ice flow (10–50 ma−1) at the onset zones of the Denman and Totten Glaciers is ‘steered’ around the highland blocks, with basal melting potentially occurring in the floors of the dissecting fjords47 (Fig. 2f). Therefore, the scale, orientation and local relief of the upland valleys are consistent with a cold basal ice temperature at the present day, and indeed since the incision of the local valleys within each block. A final piece of potential evidence for a cold basal thermal regime at present is the indication of internal reflection horizons that may indicate older ice is present over the southern block (Fig. 3b) and suggesting that basal melting must therefore be minimal otherwise, this ice would not have survived.

Discussion

The three criteria tested by our analysis have been met, confirming the existence of a relic landscape in central East Antarctica. We now consider the sequence of events and possible age of this landscape before addressing the implications for ice sheet history.

We propose that the following relative sequence of events is required to modify an original large-scale fluvial landscape to become a set of three remnant upland blocks separated by deep troughs. This landscape evolution scenario (Fig. 6) requires interactions between tectonic structures and climate in order to dissect the original fluvial drainage network, modify it with local ice, and subsequently preserve it beneath the continental-scale EAIS.

The original landscape (Fig. 6a) likely developed under pre-glacial conditions before, during and after the break-up of Gondwana, when fluvial incision would have been the main agent of landscape evolution on a continental scale, with deep weathering occurring in warmer, more humid areas48. It is possible that the preserved dendritic valleys originally formed part of a palaeofluvial drainage pathway that radiated from the interior highlands of East Antarctica (e.g., the Gamburtsev Mountains) to the Wilkes Land coast and potentially across to Australia prior to continental separation49. We expect this fluvial landscape to have been characterised by relatively low local relief (<800 m) and peak elevations approximately 1 km above sea level.

Jurassic–Cretaceous Gondwana break-up likely caused the reactivation of pre-existing (Late Palaeozoic–Triassic) rift structures in East Antarctica, as well as uplift of the Australo-Antarctic passive margin4950. We suggest that faults between the three highland blocks may have been activated at this time (Fig. 6b), with associated extensional/transtensional motion separating the highlands and thus breaking the inferred cross-block valley connections (Fig. 3a). The similarity in the way the nearby Highland B region is separated into multiple blocks37 suggests it may also have been affected by tectonic activity at this time. We infer that separation of the blocks must have occurred prior to local glaciation since otherwise, local ice would likely not have remained topographically steered51 due to the broad uniformity in elevation/relief across the region and would instead likely have formed a broader domed icecap.

Initial glaciation of the landscape blocks (Fig. 6c) likely saw localised warm-based valley glaciers that were steered by the inherited dendritic topography1051 to deepen the valleys to ~800 m of relief and generate U-shaped cross-sectional profiles. Peaks were likely preserved beneath cold-based ice and experienced isostatic uplift due to erosional unloading within the valleys (Fig. 5). This local glaciation, which may have been cyclical and paced by orbital forcing, likely modified the landscape over timescales on the order of 106 years. However, the duration and/or erosivity of this glaciation were unable to cause incision below sea level (Fig. 4).

The ice sheet rapidly expanded to a continental scale (Fig. 6d, e). During growth, the ice sheet likely exploited the existing tectonic structures and topographic lows separating the three highlands (Fig. 6d), further reinforcing (via glacial erosion) the isolation of the three pieces of the original land surface. Preferential steering of warm-based ice through the troughs enabled a cold-based thermal regime to become established across the entirety of the adjacent highland land surfaces5253 (i.e., not just on the peaks but also in the upland valley floors), thus enabling widespread preservation of the landscape whilst the ice in the deeper troughs and adjacent basins was contemporaneously warm-based.

Subsequent ice margin fluctuations may have occurred but none of these was sufficient to significantly alter the basal thermal regime or ice flow pattern over the highland blocks. Glaciological conditions amenable to landscape preservation therefore persisted.

Fig. 6: Inferred evolution of the ancient landscape.

a

Original fluvial land surface;

b

Gondwana break-up activates pre-existing tectonic structures, causing drainage re-organisation;

c

Local warm-based ice forms on the uplands and incises local-scale glacial valleys;

d

Ice expansion towards continental-scale ice sheet with cold-based ice rapidly developed over highlands, and warm-based ice steering through and eroding the intervening structurally-controlled valleys so they become fjords;

e

Continental-scale ice with cold-based ice on the highlands and continental-scale ice flowing slowly over the region, largely ignoring topography.

Full size image

Our proposed landscape evolution scenario indicates that the most recent erosive modification of the three highland blocks was by local-scale glaciation. Since this modification must have occurred at a time of much-reduced EAIS extent relative to the present day, a key question for EAIS history is: when did this local-scale glacial erosion last occur? Here, we evaluate the likelihood of the landscape of local glaciation being older than four key ‘thresholds’ of geological age pertinent to EAIS evolution. We consider that the land surface has a likely age of at least 14 Ma based on the following logic:

First, local-scale glacial modification is not occurring contemporarily because present-day continental-scale ice flow speed, direction41 (Fig. 2e), and locally frozen basal conditions47 (Fig. 2f) are inconsistent with the development of the valley-scale glacial geomorphology (Fig. 3).

Second, the landscape is most likely older than 1.5 Ma. This is because numerical modelling indicates the present-day cold basal state of this part of the EAIS is likely to have persisted since at least 1.5 Ma, even allowing for greater ice thickness at the Last Glacial Maximum and a range of geothermal heat fluxes54. Indeed, van Leifferinge et al.54 find that the only way to induce melting over the past 1.5 Ma would be to increase the geothermal heat flux to above 75 mW/m2, a significant shift from predicted modern rates of between 30 and 60 mW/m255. Furthermore, modelling of Pleistocene interglacials does not indicate a significant response of this sector of the EAIS during those times565758. The landscape therefore most likely pre-dates significant Pleistocene interglacials such as MIS 5e, 11c, and 31.

Third, the landscape is likely older than ca. 14 Ma. The mPWP (ca. 3.26–3.02 Ma) is likely the warmest (highest CO2 and with orbital parameters favourable to deglaciation59) period since ca. 14 Ma and is therefore the most likely candidate for significant retreat of the EAIS via the Denman and Totten Glaciers, which would be necessary for modification of the landscape by local ice flow. However, ice sheet model intercomparison indicates that most models do not simulate such retreat in the Aurora and Schmidt Subglacial Basins under peak mPWP climate conditions60 unless marine ice cliff instability is invoked25. We therefore suggest that these surfaces predate the EAIS switch to cold polar conditions at ca. 14 Ma17.

Finally, there remains a possibility that the landscape is older than ca. 34 Ma. This is because after East Antarctica was initially glaciated following the EOT, the extent and chronology of fluctuations of the ice margin are uncertain1114. However, there is evidence from the Aurora subglacial basin that ice was reaching the continental shelf prior to the EOT29. Thus, while the landscape may have formed under restricted EAIS conditions between 34 and 14 Ma, it may alternatively (or additionally) reflect glacial conditions immediately prior to the establishment of the first continental-scale ice sheet at the EOT2910. Therefore, on the basis of this evaluation of the existing understanding of Antarctic Cenozoic glacial history, we conclude that the Highland A land surface was last incised by local valley glaciers at least prior to ca. 14 Ma and possibly immediately prior to the initial growth of continental-scale ice at 34 Ma.

The presence of this preserved landscape has implications for understanding long-term EAIS behaviour. As previously shown in other regions of Antarctica93132343561, ancient landscapes can survive beneath the EAIS for periods of at least several million years and likely much longer. In order for this to be possible, the transition between two glacial states—(1) local-scale glaciation with warm-based erosive glaciers cutting U-shaped valleys and (2) a continental-scale ice sheet which is locally cold-based and preserves the landscape of local glaciation—must occur rapidly. Previous modelling of ice sheet growth has demonstrated rapid ice expansion and a snap-change in basal thermal regime from warm- to cold-based occurred in the Gamburtsev Mountains region at the EOT10, meaning they were eroded very little62 and we propose a similar mechanism to have occurred across Highland A. Such rapid changes in ice extent and dynamics are supported by recent ice sheet-climate modelling63, where over a period of a few million years, either side of the EOT simulated ice volumes oscillate between local- and continental-scale in a quasi-binary and geologically near-instantaneous manner. During the growth of the EAIS, it is therefore possible that ice existed in a relatively small local form capable of eroding the Highland A land surface or at a continental scale that preserved the landscape of local glaciation, with little time in between.

The preserved landscape we identify is somewhat unique given that it lies on the periphery of one of the ‘great basins’ of East Antarctica, which are of substantial interest as potential locations of past/future margin retreat via marine ice sheet instability5664. The survival of the relic landscape implies that there is long-term thermal stability over these blocks of ancient terrain, with cold-based ice being the average glacial condition since landscape formation. If the ice margin were to have retreated substantially into the Aurora and Schmidt Subglacial Basins during significant past warm intervals such as the mPWP, independent ice masses remaining on the highland blocks would likely have been drained radially by outlet glaciers. We do not see evidence of radial drainage, but instead a more complex dendritic pattern and therefore deem it likely that the ice margin has never retreated far enough to enable this region to become an independent ice centre. If our interpretation is correct, this landscape may have survived because any retreat of this sector of the EAIS during past warm periods did not extend as far inland as these highland blocks, which are ~300 km inland of the modern-day grounding line. We note, however, that this landscape is likely not sensitive to, and therefore does not preclude, fluctuations within 300 km of the modern ice margin33.

Our analysis shows that local-scale glaciation modified a pre-existing fluvial landscape. However, it is not clear whether this would have occurred under the single initial state of early ice growth at or before the EOT or whether this reflects average glacial conditions throughout the Oligocene and early Miocene prior to the onset of full polar ice sheet conditions at ca. 14 Ma17. Sampling of subglacial bedrock and potentially subglacial sediment would help better constrain the landscape age and we propose that given this region may not have been subject to significant erosion or basal transport for millions of years, it may be a promising area to recover a deep ice and subglacial sediment core in order to understand glacial conditions during those times. In addition, a more detailed local RES survey would help resolve the elevation and geometry of potential cirque floors231, which could then be used to constrain the climate that supported local-scale glaciers on the landscape. Such a survey could be directed by our ice surface feature mapping (Fig. 3a) and the hints of cirques visible in the existing RES survey data (Fig. 3b).

Finally, given this discovery of an ancient landscape hidden in plain sight, and that of others9324061, we propose that there will be other similar, as yet undiscovered, ancient landscapes beneath the EAIS. Systematic mapping of these landscapes, along with sampling via direct bed access, may provide significant new insights into the growth and past fluctuations of the EAIS. They may also facilitate a greater understanding of the long-term landscape evolution and tectonic history of interior East Antarctica, which is largely inferred from potential field data with little direct supporting evidence.

In conclusion, we use satellite remote sensing of ice surface morphology and imaging of the ice-bed interface from RES surveying to map a previously undiscovered landscape in the Aurora-Schmidt basins in the upper Denman and Totten glacier catchments of East Antarctica. We find that the landscape contains a valley network with a dendritic (i.e., fluvial) planform and was once a continuous land surface that has since been dissected by deep glacial troughs that may inherit earlier tectonic weaknesses. Cross profiles of the terrain from RES data indicate that this fluvial landscape has been modified by local-scale glaciers and has not been modified by ice flowing at a continental scale. The age of the land surface is uncertain, but it is likely to predate the mid-Miocene (ca. 14 Ma) and perhaps dates from the transition from warm to glacial conditions in Antarctica following the EOT at ca. 34 Ma. The survival of the landscape implies a long-term stable basal thermal regime and that the ice margin is unlikely to have retreated as far inland as this locality during warm periods of the last 14 Ma. The landscape is more consistent with the idea that this region of East Antarctica requires greater warming than seen in (e.g.,) the mPWP for retreat to reach this far into the Aurora-Schmidt subglacial basins.

Given that modern atmospheric CO2 and temperature conditions have reached levels unprecedented since the Pliocene6566, we are now on course to develop atmospheric conditions similar to those that prevailed between 34 and 14 Ma (ca. 3–7 °C warmer than present), with CO2 conditions reaching above 500 ppm between now and 2100 under continued fossil fuel burning67. Under such anthropogenic forcing and the associated global heating, our work suggests the EAIS may eventually retreat enough that local ice caps would once again exist on ancient preserved landscapes at the margin of the Aurora and Schmidt basins.

Methods

Geomorphological mapping

Following the approach of Ross et al.32 and Jamieson et al.40, we use satellite remote sensing to map the likely ridge and valley structure of the Highland A region. We use contrasts in the RADARSAT42 and REMA68 ice surface datasets to infer minor changes in ice surface slope, which are a consequence of buried topography. RADARSAT imagery records the brightness or intensity of the reflection of the satellite-emitted radar signal, which depends on the angle of the surface slope. Typically, subglacial valleys are positioned below darker areas (lower reflection intensity) and subglacial ridges below lighter areas (higher reflection intensity) in the RADARSAT image.

REMA is an ice surface DEM; we calculate the plan and profile curvature of the DEM to identify breaks in slope that are indicative of subglacial valleys and ridges. In some locations, the change in ice surface brightness/curvature is abrupt, for example, where potentially sharp ridges or peaks are present beneath the ice. In other areas, the contrast is significantly muted, making mapping more challenging and less certain, particularly when attempting to determine valley/ridge orientations.

Using the processed REMA and RADARSAT images, we digitise the locations of valleys and ridges where clear breaks in slope or changes in radar reflection intensity occur. Connections between valleys are inferred such that if a valley is continuous across an entire block, it is then connected to the next block if this block hosts an apparently contiguous valley on the other side of the intervening trough. Within each topographic block, the spatial statistics for the mapped valleys are calculated, including their length, the ridge-valley spacing, and the drainage density (length of valley per unit area).

The location and form of the valleys are then confirmed using the available HiCARS RES survey data373839. Data were focused following the approaches described in Peters et al.69. Using the focussed RES lines, we automatically extract ridge and valley locations along each RES line as the points of maximum and minimum elevation within a 10 km-wide moving window, respectively. The ridge and valley locations are used to ground-truth the locations of the features mapped from the REMA and RADARSAT data. We also determine the mean relief (i.e., the vertical distance between individual ridges and valleys) within the three landscape blocks. In addition to identifying the locations of ridges and valleys, we qualitatively examine the RES lines to identify whether the valleys are U-shaped or V-shaped and thus assess whether the landscape records the erosive imprint of either local valley glaciers or rivers.

Flexural modelling

We use flexural modelling to quantify the effect of large-scale valley incision on the elevation of the highland blocks. After isostatically adjusting bed elevations for the removal of the modern ice load70, we construct an accordant surface between the highland peaks. Peaks are identified as high points in the landscape within a circular moving window of diameter 10 km, and a smooth surface representing the land surface in the absence of valley incision is interpolated between these peaks (Fig. 5). This summit accordance surface also extends across the two wide troughs that separate the three highland blocks and across the surrounding low-lying terrain. The distribution of erosion is, in turn, inferred to be the difference between the accordant surface and the (ice-free) bed topography7172.

We use a 2D elastic plate model to compute the flexure, w, induced by the removal of the eroded material of thickness h. The general 2D flexure equation is expressed as

$${\nabla }^{2}\left[D\left(x,y\right){\nabla }^{2}w\left(x,y\right)\right]+\left({\rho }_{{mantle}}-{\rho }_{{infill}}\right)\,g\,w\left(x,y\right)=\left({\rho }_{{load}}-{\rho }_{{displace}}\right)\,g\,h\left(x,y\right)$$

(1)

where

$$D\left(x,y\right)=\frac{{{ET}}_{e}{\left(x,y\right)}^{3}}{12\left(1-{\nu }^{2}\right)}$$

(2)

The parameter Te represents the effective elastic thicknesses (a proxy for the flexural rigidity, D) of the East Antarctic lithosphere. For simplicity, we assume a spatially uniform Te, allowing Eq. (1) to be solved analytically using a fast Fourier transform of the erosional unload and convolution with a 2D flexural isostatic response function46. We assume a Young’s modulus (E) of 100 GPa, a Poisson ratio (ν) of 0.25, a gravitational acceleration (g) of 9.81 m s−2 and densities of 2670 and 3330 kg m−3 for the eroded material (ρload) and mantle (ρmantle), respectively. Note that for erosional unloading, the material infilling and displaced by the flexure is air, so ρinfill and ρdisplace are set to 0 kg m−3. We perform the calculation for a range of constant Te values to test the sensitivity of the flexural response to lithospheric rigidity and compare the observed land surface elevations to the predictions of our elastic plate model for different Te values.

Data availability

The data used for this work is the radio-echo sounding data from the ICECAP project, which is openly accessible via the Blankenship 2017 references3839 (HICARS1: https://doi.org/10.5067/F5FGUT9F5089; HICARS2: https://doi.org/10.5067/9EBR2T0VXUDG). The mapping data generated in this study (Fig. 3a) are openly available as GIS shapefiles at https://doi.org/10.5281/zenodo.815922373. Source data are provided with this paper—these relate to the data that underlies Figs. 3c and 4. Source data are provided with this paper.

References

1. Carter, A., Riley, T. R., Hillenbrand, C.-D. & Rittner, M. Widespread Antarctic glaciation during the Late Eocene. Earth Planet. Sci. Lett. 458, 49–57 (2017).

Google Scholar

2. Barr, I. D. et al. 60 million years of glaciation in the Transantarctic Mountains. Nat. Commun. 13, 5526 (2022).

PubMed PubMed Central Google Scholar

3. Coxall, H. K., Wilson, P. A., Palike, H., Lear, C. H. & Backman, J. Rapid stepwise onset of Antarctic glaciation and deeper calcite compensation in the Pacific Ocean. Nature 433, 53–57 (2005).

PubMed Google Scholar

4. Zachos, J. C., Rea, D. K., Seto, K., Nomura, R. & Niitsuma, N. Paleogene and Early Neogene deep water paleoceanography of the Indian Ocean as determined from benthic foraminifer stable carbon and oxygen isotope records. Synthesis of Results from Scientific Drilling in the Indian Ocean (AGU, 1992).
5. DeConto, R. M. & Pollard, D. Rapid Cenozoic glaciation of Antarctica induced by declining atmospheric CO2. Nature 421, 245–249 (2003).

PubMed Google Scholar

6. Tripati, A., Backman, J., Elderfield, H. & Ferretti, P. Eocene bipolar glaciation associated with global carbon cycle changes. Nature 436, 341–346 (2005).

PubMed Google Scholar

7. Zachos, J. C. & Kump, L. R. Carbon cycle feedbacks and the initiation of Antarctic glaciation in the earliest Oligocene. Glob. Planet. Change 47, 51–66 (2005).

Google Scholar

8. Sauermilch, I. et al. Gateway-driven weakening of ocean gyres leads to Southern Ocean cooling. Nat. Commun. 12, 6465 (2021).

PubMed PubMed Central Google Scholar

9. Rose, K. C. et al. Early East Antarctic Ice Sheet growth recorded in the landscape of the Gamburtsev Subglacial Mountains. Earth Planet. Sci. Lett. 375, 1–12 (2013).

Google Scholar

10. Jamieson, S. S. R., Sugden, D. E. & Hulton, N. R. J. The evolution of the subglacial landscape of Antarctica. Earth Planet. Sci. Lett. 293, 1–27 (2010).

Google Scholar

11. Naish, T. R. et al. Orbitally induced oscillations in the East Antarctic ice sheet at the Oligocene/Miocene boundary. Nature 413, 719–723 (2001).

PubMed Google Scholar

12. Naish, T. et al. Obliquity-paced Pliocene West Antarctic ice sheet oscillations. Nature 458, 322–328 (2009).

PubMed Google Scholar

13. Gasson, E., DeConto, R. M., Pollard, D. & Levy, R. H. Dynamic Antarctic ice sheet during the early to mid-Miocene. Proc. Natl Acad. Sci. USA 113, 3459–3464 (2016).

PubMed PubMed Central Google Scholar

14. Levy, R. et al. Antarctic ice sheet sensitivity to atmospheric CO2 variations in the early to mid-Miocene. Proc. Natl Acad. Sci. USA 113, 3453 (2016).

PubMed PubMed Central Google Scholar

15. Passchier, S. et al. Early and middle Miocene Antarctic glacial history from the sedimentary facies distribution in the AND-2A drill hole, Ross Sea, Antarctica. GSA Bull. 123, 2352–2365 (2011).

Google Scholar

16. Pérez, L. F. et al. Early and middle Miocene ice sheet dynamics in the Ross Sea: results from integrated core-log-seismic interpretation. GSA Bull. 134, 348–370 (2021).

Google Scholar

17. Sugden, D. E. & Denton, G. Cenozoic landscape evolution of the Convoy Range to Mackay Glacier area, Transantartic Mountains: onshore to offshore synthesis. Bull. Geol. Soc. Am. 116, 840–857 (2004).

Google Scholar

18. Williams, T. et al. Evidence for iceberg armadas from East Antarctica in the Southern Ocean during the late Miocene and early Pliocene. Earth Planet. Sci. Lett. 290, 351–361 (2010).

Google Scholar

19. Cook, C. P. et al. Dynamic behaviour of the East Antarctic ice sheet during Pliocene warmth. Nat. Geosci. 6, 765–769 (2013).

Google Scholar

20. Bertram, R. A. et al. Pliocene deglacial event timelines and the biogeochemical response offshore Wilkes Subglacial Basin, East Antarctica. Earth Planet. Sci. Lett. 494, 109–116 (2018).

Google Scholar

21. Blackburn, T. et al. Ice retreat in Wilkes Basin of East Antarctica during a warm interglacial. Nature 583, 554–559 (2020).

PubMed Google Scholar

22. Wilson, D. J. et al. Ice loss from the East Antarctic Ice Sheet during late Pleistocene interglacials. Nature 561, 383–386 (2018).

PubMed Google Scholar

23. Crotti, I. et al. Wilkes subglacial basin ice sheet response to Southern Ocean warming during late Pleistocene interglacials. Nat. Commun. 13, 5328 (2022).

PubMed PubMed Central Google Scholar

24. Iizuka, M. et al. Multiple episodes of ice loss from the Wilkes Subglacial Basin during the Last Interglacial. Nat. Commun. 14, 2129 (2023).

PubMed PubMed Central Google Scholar

25. DeConto, R. M. et al. The Paris Climate Agreement and future sea-level rise from Antarctica. Nature 593, 83–89 (2021).

PubMed Google Scholar

26. Edwards, T. L. et al. Revisiting Antarctic ice loss due to marine ice-cliff instability. Nature 566, 58–64 (2019).

PubMed Google Scholar

27. Stokes, C. R. et al. Response of the East Antarctic Ice Sheet to past and future climate change. Nature 608, 275–286 (2022).

PubMed Google Scholar

28. Miller, K. G. et al. Cenozoic sea-level and cryospheric evolution from deep-sea geochemical and continental margin records. Sci. Adv. 6, eaaz1346 (2020).

PubMed PubMed Central Google Scholar

29. Gulick, S. P. S. et al. Initiation and long-term instability of the East Antarctic Ice Sheet. Nature 552, 225–229 (2017).

PubMed Google Scholar

30. Porter, S. C. Some geological implications of average Quaternary glacial conditions. Quat. Res. 32, 245–261 (1989).

Google Scholar

31. Bo, S. et al. The Gamburtsev mountains and the origin and early evolution of the Antarctic Ice Sheet. Nature 459, 690–693 (2009).

PubMed Google Scholar

32. Ross, N. et al. The Ellsworth Subglacial Highlands: inception and retreat of the West Antarctic Ice Sheet. Geol. Soc. Am. Bull. 126, 3–15 (2014).

Google Scholar

33. Aitken, A. R. A. et al. Repeated large-scale retreat and advance of Totten Glacier indicated by inland bed erosion. Nature 533, 385–389 (2016).

PubMed Google Scholar

34. Paxman, G. J. G. et al. Bedrock erosion surfaces record former East Antarctic Ice Sheet extent. Geophys. Res. Lett. 45, 4114–4123 (2018).

Google Scholar

35. Franke, S. et al. Preserved landscapes underneath the Antarctic Ice Sheet reveal the geomorphological history of Jutulstraumen Basin. Earth Surf. Process. Landf. 46, 2728–2745 (2021).

Google Scholar

36. Paxman, G. J. G. et al. Reconstructions of Antarctic topography since the Eocene–Oligocene boundary. Palaeogeogr. Palaeoclimatol. Palaeoecol. 535, 109346 (2019).

Google Scholar

37. Young, D. A. et al. A dynamic early East Antarctic Ice Sheet suggested by ice-covered fjord landscapes. Nature 474, 72–75 (2011).

PubMed Google Scholar

38. Blankenship, D. D. et al. IceBridge HiCARS 1 L1B Time-Tagged Echo Strength Profiles (NASA National Snow and Ice Data Center, 2017).
39. Blankenship, D. D. et al. IceBridge HiCARS 2 L1B Time-Tagged Echo Strength Profiles (NASA National Snow and Ice Data Center, 2017).
40. Jamieson, S. S. R. et al. An extensive subglacial lake and canyon system in Princess Elizabeth Land, East Antarctica. Geology 44, 87–90 (2016).
41. Mouginot, J., Rignot, E. & Scheuchl, B. Continent-wide, interferometric SAR phase, mapping of Antarctic Ice velocity. Geophys. Res. Lett. 46, 9710–9718 (2019).

Google Scholar

42. Jezek, K., Curlander, J. C., Carsey, F., Wales, C. & Barry, R. RAMP AMM-1 SAR Image Mosaic of Antarctica, Version 2 (NASA National Snow and Ice Data Center, 2002).
43. Grau Galofre, A. & Jellinek, A. M. The geometry and complexity of spatial patterns of terrestrial channel networks: distinctive fingerprints of erosional regimes. J. Geophys. Res. Earth Surf. 122, 1037–1059 (2017).

Google Scholar

44. Sugden, D. E., Denton, G. H. & Marchant, D. R. Subglacial meltwater channel systems and ice sheet overriding, Asgard Range, Antarctica. Geogr. Ann. A 73, 109–121 (1991).

Google Scholar

45. Swain, C. J. & Kirby, J. F. Effective elastic thickness map reveals subglacial structure of East Antarctica. Geophys. Res. Lett. 48, e2020GL091576 (2021).

Google Scholar

46. Watts, A. B. Isostasy and Flexure of the Lithosphere (Cambridge University Press, 2001).
47. McCormack, F. S. et al. Fine-scale geothermal heat flow in Antarctica can increase simulated subglacial melt estimates. Geophys. Res. Lett. 49, e2022GL098539 (2022).

Google Scholar

48. Sugden, D. E. & Jamieson, S. S. R. The pre-glacial landscape of Antarctica. Scott. Geogr. J. 134, 203–223 (2018).

Google Scholar

49. Maritati, A., Halpin, J. A., Whittaker, J. M., Daczko, N. R. & Wainman, C. C. Provenance of Upper Jurassic–Lower Cretaceous strata in the Mentelle Basin, southwestern Australia, reveals a trans-Gondwanan fluvial pathway. Gondwana Res. 93, 128–141 (2021).

Google Scholar

50. Aitken, A. R. A. et al. The subglacial geology of Wilkes Land, East Antarctica. Geophys. Res. Lett. 41, 2014GL059405 (2014).

Google Scholar

51. Kessler, M. A., Anderson, R. S. & Briner, J. P. Fjord insertion into continental margins driven by topographic steering of ice. Nat. Geosci. 1, 365–369 (2008).

Google Scholar

52. Stern, T. A., Baxter, A. K. & Barrett, P. J. Isostatic rebound due to glacial erosion within the Transantarctic Mountains. Geology 33, 221 (2005).

Google Scholar

53. Sugden, D. E. et al. Emergence of the Shackleton Range from beneath the Antarctic Ice Sheet due to glacial erosion. Geomorphology 208, 190–199 (2014).

Google Scholar

54. Van Leifferinge, B. et al. Promising Oldest Ice sites in East Antarctica based on thermodynamical modelling. Cryosphere 12, 2773–2787 (2018).

Google Scholar

55. Lösing, M. & Ebbing, J. Predicting geothermal heat flow in Antarctica with a machine learning approach. J. Geophys. Res. Solid. Earth 126, e2020JB021499 (2021).

Google Scholar

56. DeConto, R. M. & Pollard, D. Contribution of Antarctica to past and future sea-level rise. Nature 531, 591–597 (2016).

PubMed Google Scholar

57. DeConto, R. M., Pollard, D. & Kowalewski, D. Modeling Antarctic ice sheet and climate variations during Marine Isotope Stage 31. Glob. Planet. Change 88–89, 45–52 (2012).

Google Scholar

58. Sutter, J. et al. Modelling the Antarctic Ice Sheet across the mid-Pleistocene transition—implications for Oldest Ice. Cryosphere 13, 2023–2041 (2019).

Google Scholar

59. Rae, J. W. B. et al. Atmospheric CO2 over the past 66 million years from marine archives. Annu. Rev. Earth Planet. Sci. 49, 609–641 (2021).

Google Scholar

60. de Boer, B., Lourens, L. J. & van de Wal, R. S. W. Persistent 400,000-year variability of Antarctic ice volume and the carbon cycle is revealed throughout the Plio-Pleistocene. Nat. Commun. 5, 2999 (2014).

PubMed Google Scholar

61. Rose, K. C. et al. Ancient pre-glacial erosion surfaces preserved beneath the West Antarctic Ice Sheet. Earth Surf. Dyn. 3, 139–152 (2015).

Google Scholar

62. Cox, S. E., Thomson, S. N., Reiners, P. W., Hemming, S. R. & van de Flierdt, T. Extremely low long-term erosion rates around the Gamburtsev Mountains in interior East Antarctica. Geophys. Res. Lett. 37, L22307 (2010).
63. Van Breedam, J., Huybrechts, P. & Crucifix, M. Modelling evidence for late Eocene Antarctic glaciations. Earth Planet. Sci. Lett. 586, 117532 (2022).

Google Scholar

64. Morlighem, M. et al. Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet. Nat. Geosci. 13, 132–137 (2020).

Google Scholar

65. Haywood, A. M., Dowsett, H. J. & Dolan, A. M. Integrating geological archives and climate models for the mid-Pliocene warm period. Nat. Commun. 7, 10646 (2016).

PubMed PubMed Central Google Scholar

66. McClymont, E. L. et al. Lessons from a high-CO2 world: an ocean view from ~3 million years ago. Clim. Past 16, 1599–1615 (2020).

Google Scholar

67. IPCC. Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, 2021).
68. Howat, I. M., Porter, C., Smith, B. E., Noh, M.-J. & Morin, P. The reference elevation model of Antarctica. Cryosphere 13, 665–674 (2019).

Google Scholar

69. Peters, M. E. et al. Along-track focusing of airborne radar sounding data from West Antarctica for improving basal reflection analysis and layer detection. IEEE Trans. Geosci. Remote. Sens. 45, 2725–2736 (2007).

Google Scholar

70. Paxman, G. J. G., Austermann, J. & Hollyday, A. Total isostatic response to the complete unloading of the Greenland and Antarctic Ice Sheets. Sci. Rep. 12, 11399 (2022).

PubMed PubMed Central Google Scholar

71. Small, E. E. & Anderson, R. S. Pleistocene relief production in Laramide mountain ranges, western United States. Geology 26, 123–126 (1998).

Google Scholar

72. Champagnac, J. D., Molnar, P., Anderson, R. S., Sue, C. & Delacou, B. Quaternary erosion-induced isostatic rebound in the western Alps. Geology 35, 195–198 (2007).

Google Scholar

73. Jamieson, S. S. R. et al. Subglacial valleys and ridges of Highland A, East Antarctica (1.0.0). Zenodo https://doi.org/10.5281/zenodo.8159223 (2023).
74. Rignot, E., Mouginot, J. & Scheuchl, B. Antarctic grounding line mapping from differential satellite radar interferometry. Geophys. Res. Lett. 38, L10504 (2011).

Google Scholar

75. Westerhold, T. et al. An astronomically dated record of Earth’s climate and its predictability over the last 66 million years. Science 369, 1383–1387 (2020).

PubMed Google Scholar

Download references

Acknowledgements

We thank Felicity McCormack for supplying modelled basal melt data. Funding for geophysical data acquisition was provided by NSF grants ANT-0733025 (D.D.B., D.A.Y.), ANT-1443690 (D.D.B., D.A.Y.), NASA grant NNX09AR52G (D.D.B., D.A.Y.), the G. Unger Vetlesen Foundation (D.D.B., D.A.Y.) and the UK NERC grant NE/D003733/1 (M.J.S.). G.J.G.P. was funded by UK NERC grant NE/L002590/1. For the purpose of open access, the author has applied a Creative Commons Attribution (CC BY) license to any Author Accepted Manuscript version arising from this submission. This is University of Texas Institute of Geophysics contribution 3963.

Author information

1. These authors contributed equally: Stewart S. R. Jamieson, Neil Ross.

Authors and Affiliations

1. Department of Geography, Durham University, Durham, DH1 3LE, UK

Stewart S. R. Jamieson, Guy J. G. Paxman & Fiona J. Clubb

2. School of Geography, Politics and Sociology, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK

Neil Ross

3. University of Texas Institute for Geophysics, Jackson School of Geosciences, University of Texas at Austin, Austin, USA

Duncan A. Young, Shuai Yan & Donald D. Blankenship

4. Department of Geosciences, Jackson School of Geosciences, University of Texas at Austin, Austin, USA

Shuai Yan

5. Scripps Institute for Oceanography, University of California at San Diego, San Diego, USA

Jamin Greenbaum

6. Tremough House, University of Exeter, Penryn Campus, Penryn, Cornwall, TR10 9FE, UK

Martin J. Siegert

Contributions

S.S.R.J. and N.R. conceived the work. D.A.Y., D.D.B., M.J.S. and J.G. provided geophysical survey data. S.S.R.J., G.J.G.P., F.J.C., D.A.Y. and S.Y. conducted the analysis. S.S.R.J., N.R., G.J.G.P., F.J.C. and M.J.S. drafted the manuscript and S.S.R.J., N.R., G.J.G.P., F.J.C., D.A.Y., S.Y., D.D.B., J.G. and M.J.S. provided revisions to the manuscript.

Corresponding author

Correspondence to Stewart S. R. Jamieson.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks Paul Hoffman, Robert Larter and the other anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Peer Review File

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Jamieson, S.S.R., Ross, N., Paxman, G.J.G. et al. An ancient river landscape preserved beneath the East Antarctic Ice Sheet. Nat Commun 14, 6507 (2023). https://doi.org/10.1038/s41467-023-42152-2

Download citation

• Received10 May 2023
• Accepted02 October 2023
• Published24 October 2023
• DOIhttps://doi.org/10.1038/s41467-023-42152-2

Share this article

Anyone you share the following link with will be able to read this content:

Provided by the Springer Nature SharedIt content-sharing initiative

Subjects

• Geomorphology
• Geophysics

This article is cited by

The role of mountains in shaping the global meridional overturning circulation

1. Haijun Yang
2. Rui Jiang
3. Jianping Huang

Nature Communications (2024)

Download PDF Associated content Collection The Top 25 Earth, Environmental, and Planetary Sciences Articles of 2023 • Abstract • Introduction • Results • Discussion • Methods • Data availability • References • Acknowledgements • Author information • Ethics declarations • Peer review • Additional information • Supplementary information • Source data • Rights and permissions • About this article

• • Abstract
• • Introduction
• • Results
• • Discussion
• • Methods
• • Data availability
• • References
• • Acknowledgements
• • Author information
• • Ethics declarations
• • Peer review
• • Additional information
• • Supplementary information
• • Source data
• • Rights and permissions
• • About this article