No abstract

Unravelling the long‐term evolution of the subglacial landscape of Antarctica is vital for understanding past ice sheet dynamics and stability, particularly in marine‐based sectors of the ice sheet. Here we model the evolution of the bedrock topography beneath the Recovery catchment, a sector of the East Antarctic Ice Sheet characterized by fast‐flowing ice streams that occupy overdeepened subglacial troughs. We use 3‐D flexural models to quantify the effect of erosional unloading and mechanical unloading associated with motion on border faults in driving isostatic bedrock uplift of the Shackleton Range and Theron Mountains, which are flanked by the Recovery, Slessor, and Bailey ice streams. Inverse spectral (free‐air admittance) and forward modeling of topography and gravity anomaly data allow us to constrain the effective elastic thickness of the lithosphere (Te) in the Shackleton Range region to ~20 km. Our models indicate that glacial erosion, and the associated isostatic rebound, has driven 40–50% of total peak uplift in the Shackleton Range and Theron Mountains. A further 40–50% can be attributed to motion on normal fault systems of inferred Jurassic and Cretaceous age. Our results indicate that the flexural effects of glacial erosion play a key role in mountain uplift along the East Antarctic margin, augmenting previous findings in the Transantarctic Mountains. The results suggest that at 34 Ma, the mountains were lower and the bounding valley floors were close to sea level, which implies that the early ice sheet in this region may have been relatively stable.

The East Antarctic Ice Sheet (EAIS) is underlain by a series of low‐lying subglacial sedimentary basins. The extent, geology, and basal topography of these sedimentary basins are important boundary conditions governing the dynamics of the overlying ice sheet. This is particularly pertinent for basins close to the grounding line wherein the EAIS is grounded below sea level and therefore potentially vulnerable to rapid retreat. Here we analyze newly acquired airborne geophysical data over the Pensacola‐Pole Basin (PPB), a previously unexplored sector of the EAIS. Using a combination of gravity and magnetic and ice‐penetrating radar data, we present the first detailed subglacial sedimentary basin model for the PPB. Radar data reveal that the PPB is defined by a topographic depression situated ~500 m below sea level. Gravity and magnetic depth‐to‐source modeling indicate that the southern part of the basin is underlain by a sedimentary succession 2–3 km thick. This is interpreted as an equivalent of the Beacon Supergroup and associated Ferrar dolerites that are exposed along the margin of East Antarctica. However, we find that similar rocks appear to be largely absent from the northern part of the basin, close to the present‐day grounding line. In addition, the eastern margin of the basin is characterized by a major geological boundary and a system of overdeepened subglacial troughs. We suggest that these characteristics of the basin may reflect the behavior of past ice sheets and/or exert an influence on the present‐day dynamics of the overlying EAIS.

East Antarctica hosts large subglacial basins into which the East Antarctic Ice Sheet (EAIS) likely retreated during past warmer climates. However, the extent of retreat remains poorly constrained, making quantifying past and predicted future contributions to global sea level rise from these marine basins challenging. Geomorphological analysis and flexural modeling within the Wilkes Subglacial Basin are used to reconstruct the ice margin during warm intervals of the Oligocene-Miocene. Flat-lying bedrock plateaus are indicative of an ice sheet margin positioned >400-500 km inland of the modern grounding zone for extended periods of the Oligocene-Miocene, equivalent to a 2-m rise in global sea level. Our findings imply that if major EAIS retreat occurs in the future, isostatic rebound will enable the plateau surfaces to act as seeding points for extensive ice rises, thus limiting extensive ice margin retreat of the scale seen during the early EAIS. Plain Language SummaryThe Wilkes Subglacial Basin is a large, low-lying topographic depression situated beneath the Antarctic Ice Sheet. Because the land surface of the basin is currently situated below sea level, it is a potential site of ice sheet collapse and rapid retreat in a warming world. Understanding this landscape and how it has evolved through time in relation to past climate and sea level is therefore key to understanding the future dynamics of this part of the ice sheet. Here we report the discovery, using ice-penetrating radar data sets, of extensive subglacial bedrock plateaus within the Wilkes Subglacial Basin. We analyze the geomorphology of these plateau surfaces and reconstruct the evolution of the subglacial landscape through time. Our results indicate that this part of the Wilkes Subglacial Basin was free of ice for extensive and prolonged periods of time during the early stages of ice sheet development. These constraints on past ice sheet extent, together with our landscape reconstruction, can be used by the ice sheet modeling community to better understand the likely future dynamics of this part of the Antarctic Ice Sheet.

Preface: The East Antarctic Ice Sheet (EAIS) contains the vast majority of Earth's glacier ice (~52 metres sea-level equivalent), but is often viewed as less vulnerable to global warming than the West Antarctic or Greenland ice sheets. However, some regions of the EAIS have lost mass over recent decades, prompting the need to re-evaluate its sensitivity to climate change. Here we review the EAIS's response to past warm periods, synthesise current observations of change, and evaluate future projections. Some marine-based catchments that underwent significant mass loss during past warm periods are currently losing mass, but most projections indicate increased accumulation across the EAIS over the 21st Century, keeping the ice sheet broadly in balance. Beyond 2100, high emissions scenarios generate increased ice discharge and potentially several metres of sea-level rise within just a few centuries, but substantial mass loss could be averted if the Paris Agreement to limit warming below 2°C is satisfied.

Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. Alpine-style glacial erosion in the Gamburtsevs and outlet glacier-type selective linear erosion in the Lambert Rift, part of the East Antarctic Rift System. 3D flexural models indicate that valley incision has contributed ca. 500 m of peak uplift in the Gamburtsevs and up to 1.2 km in the Lambert Rift, which is consistent with the present-day elevation of Oligocene-Miocene glaciomarine sediments. Overall, we find that 17-25% of Gamburtsev peak uplift can be explained by erosional unloading. These relatively low values are typical of temperate mountain ranges, suggesting that most of the valley incision in the Gamburtsevs occurred prior to widespread glaciation at 34 Ma. The pre-incision topography of the Gamburtsevs lies at 2-2.5 km above sea-level, confirming that they were a key inception point for the development of the East Antarctic Ice Sheet. Tectonic and/or dynamic processes were therefore responsible for ca. 80% of the elevation of the modern Gamburtsev Subglacial Mountains.

Reconstructions of the bedrock topography of Antarctica since the Eocene‐Oligocene Boundary (approximately 34 Ma) provide important constraints for modeling Antarctic ice sheet evolution. This is particularly important in regions where the bedrock lies below sea level, since in these sectors the overlying ice sheet is thought to be most susceptible to past and future change. Here we use 3‐D flexural modeling to reconstruct the evolution of the topography of the Wilkes Subglacial Basin (WSB) and Transantarctic Mountains (TAM) in East Antarctica. We estimate the spatial distribution of glacial erosion beneath the East Antarctic Ice Sheet, and restore this material to the topography, which is also adjusted for associated flexural isostatic responses. We independently constrain our post‐34 Ma erosion estimates using offshore sediment stratigraphy interpretations. Our reconstructions provide a better‐defined topographic boundary condition for modeling early East Antarctic Ice Sheet history. We show that the majority of glacial erosion and landscape evolution occurred prior to 14 Ma, which we interpret to reflect more dynamic and erosive early ice sheet behavior. In addition, we use closely spaced 2‐D flexural models to test previously proposed hypotheses for a flexural origin of the TAM and WSB. The pre‐34 Ma topography shows lateral variations along the length of the TAM and WSB that cannot be explained by uniform flexure along the front of the TAM. We show that some of these variations may be explained by additional flexural uplift along the south‐western flank of the WSB and the Rennick Graben in northern Victoria Land.

Ice sheet behavior is strongly influenced by the bed topography. However, the effect of the progressive temporal evolution of Antarctica's subglacial landscape on the sensitivity of the Antarctic Ice Sheet (AIS) to climatic and oceanic change has yet to be fully quantified. Here we investigate the evolving sensitivity of the AIS using a series of data-constrained reconstructions of Antarctic paleotopography since glacial inception at the Eocene-Oligocene transition. We use a numerical ice sheet model to subject the AIS to schematic climate and ocean warming experiments and find that bed topographic evolution causes a doubling in ice volume loss and equivalent global sea level rise. Glacial erosion is primarily responsible for enhanced ice sheet retreat via the development of increasingly low-lying and reverse sloping beds over time, particularly within near-coastal subglacial basins. We conclude that AIS sensitivity to climate and ocean forcing has been substantially amplified by long-term landscape evolution. Plain Language Summary The Antarctic Ice Sheet is situated above a large landmass, the geometry of which is an important control on the behavior of the ice sheet and how it responds to climatic change. However, Antarctica's subglacial landscape has evolved significantly since the formation of the first ice sheets approximately 34 million years ago, which implies that the sensitivity of the ice sheet to climate change may also have changed over time. Our ice sheet model experiments show that the progressive evolution of Antarctica's bed topography has enhanced the sensitivity of the Antarctic Ice Sheet to climate and ocean warming, meaning that a greater volume of ice is lost for a given warming scenario when using the modern topography compared to past topographies. In particular, the lowering of bed elevations by glacial erosion has caused a notable increase in ice sheet sensitivity within subglacial basins close to the ice sheet margin.