Search

Direct geological information in Antarctica is limited to ice free regions along the coast, high mountain ranges, or isolated nunataks. Therefore, indirect methods are required to reveal subglacial geology and heterogeneities in crustal properties, which are critical steps toward interpreting geological history. We present a 3D crustal model of density and susceptibility distribution in the Wilkes Subglacial Basin (WSB) and the Transantarctic Mountains (TAM) based on joint inversion of airborne gravity and magnetic data. The applied “variation of information” technique enforces a coupling between inferred susceptibility and density, relating these quantities to the same gravity and magnetic sources to give an enhanced inversion result. Our model reveals a large body located in the interior of the WSB interpreted as a batholithic intrusive structure, as well as a linear dense body at the margin of the Terre Adélie Craton. Density and susceptibility relationships are used to inform the interpretation of petrophysical properties and the reconstruction of the origin of those crustal bodies. The petrophysical relationship indicates that the postulated batholitic intrusion is granitic, but independent from the Granite Harbor Igneous Complex described previously in the TAM area. Emplacement of a large volume of intrusive granites can potentially elevate local geothermal heat flow significantly. Finally, we present a new conceptual tectonic model based on the inversion results, which includes development of a passive continental margin with seaward dipping basalt horizons and magmatic underplating followed by two distinct intrusive events associated with the protracted Ross Orogen. Most rocks in Antarctica are hidden beneath a thick ice sheet. Therefore, indirect techniques are required to reveal rock provinces within Earth's crust below the ice. Rocks simultaneously influence the gravity and magnetic fields through their physical properties (density and susceptibility). Here we use both the gravity and magnetic fields to reveal rock provinces beneath the ice and use the relationship between density and susceptibility of the rocks to interpret the distribution of granitic rocks of the Transantarctic Mountains and the Wilkes Subglacial Basin region in East Antarctica. Granitic rocks can lead to elevated heat flow due to radiogenic decay of minerals within the rock and thus influence the overlying ice sheet. Based on our subsurface model of rock provinces we speculate on the tectonic evolution of the region. We present a new 3D crustal density and susceptibility distribution model based on joint inversion of gravity and magnetic dataDensity and susceptibility data are used to identify crustal level intrusions and the craton marginOur new conceptual tectonic model identifies a earlier intrusive event in the interior of the Wilkes Subglacial Basin, separate from those exposed in the Transantarctic Mountains We present a new 3D crustal density and susceptibility distribution model based on joint inversion of gravity and magnetic data Density and susceptibility data are used to identify crustal level intrusions and the craton margin Our new conceptual tectonic model identifies a earlier intrusive event in the interior of the Wilkes Subglacial Basin, separate from those exposed in the Transantarctic Mountains

Understanding the recent history of Thwaites Glacier, and the processes controlling its ongoing retreat, is key to projecting Antarctic contributions to future sea-level rise. Of particular concern is how the glacier grounding zone might evolve over coming decades where it is stabilized by sea-floor bathymetric highs. Here we use geophysical data from an autonomous underwater vehicle deployed at the Thwaites Glacier ice front, to document the ocean-floor imprint of past retreat from a sea-bed promontory. We show patterns of back-stepping sedimentary ridges formed daily by a mechanism of tidal lifting and settling at the grounding line at a time when Thwaites Glacier was more advanced than it is today. Over a duration of 5.5 months, Thwaites grounding zone retreated at a rate of >2.1 km per year—twice the rate observed by satellite at the fastest retreating part of the grounding zone between 2011 and 2019. Our results suggest that sustained pulses of rapid retreat have occurred at Thwaites Glacier in the past two centuries. Similar rapid retreat pulses are likely to occur in the near future when the grounding zone migrates back off stabilizing high points on the sea floor.

The mid-Cretaceous period was one of the warmest intervals of the past 140 million years1–5, driven by atmospheric carbon dioxide levels of around 1,000 parts per million by volume6. In the near absence of proximal geological records from south of the Antarctic Circle, it is disputed whether polar ice could exist under such environmental conditions. Here we use a sedimentary sequence recovered from the West Antarctic shelf—the southernmost Cretaceous record reported so far—and show that a temperate lowland rainforest environment existed at a palaeolatitude of about 82° S during the Turonian–Santonian age (92 to 83 million years ago). This record contains an intact 3-metre-long network of in situ fossil roots embedded in a mudstone matrix containing diverse pollen and spores. A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric carbon dioxide concentrations of 1,120–1,680 parts per million by volume and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo under high levels of atmospheric carbon dioxide.

Marine ice-cliff instability (MICI) processes could accelerate future retreat of the Antarctic Ice Sheet if ice shelves that buttress grounding lines more than 800 metres below sea level are lost. The present-day grounding zones of the Pine Island and Thwaites glaciers in West Antarctica need to retreat only short distances before they reach extensive retrograde slopes. When grounding zones of glaciers retreat onto such slopes, theoretical considerations and modelling results indicate that the retreat becomes unstable (marine ice-sheet instability) and thus accelerates. It is thought that MICI is triggered when this retreat produces ice cliffs above the water line with heights approaching about 90 metres. However, observational evidence confirming the action of MICI has not previously been reported. Here we present observational evidence that rapid deglacial ice-sheet retreat into Pine Island Bay proceeded in a similar manner to that simulated in a recent modelling study, driven by MICI. Iceberg-keel plough marks on the sea-floor provide geological evidence of past and present iceberg morphology, keel depth and drift direction. From the planform shape and cross-sectional morphologies of iceberg-keel plough marks, we find that iceberg calving during the most recent deglaciation was not characterized by small numbers of large, tabular icebergs as is observed today, which would produce wide, flat-based plough marks or toothcomb-like multi-keeled plough marks. Instead, it was characterized by large numbers of smaller icebergs with V-shaped keels. Geological evidence of the form and water-depth distribution of the plough marks indicates calving-margin thicknesses equivalent to the threshold that is predicted to trigger ice-cliff structural collapse as a result of MICI. We infer rapid and sustained ice-sheet retreat driven by MICI, commencing around 12,300 years ago and terminating before about 11,200 years ago, which produced large numbers of icebergs smaller than the typical tabular icebergs produced today. Our findings demonstrate the effective operation of MICI in the past, and highlight its potential contribution to accelerated future retreat of the Antarctic Ice Sheet.

Glaciological and oceanographic observations coupled with numerical models show that warm Circumpolar Deep Water (CDW) incursions onto the West Antarctic continental shelf cause melting of the undersides of floating ice shelves. Because these ice shelves buttress glaciers feeding into them, their ocean-induced thinning is driving Antarctic ice-sheet retreat today. Here we present a multi-proxy data based reconstruction of variability in CDW inflow to the Amundsen Sea sector, the most vulnerable part of the West Antarctic Ice Sheet, during the Holocene epoch (from 11.7 thousand years ago to the present). The chemical compositions of foraminifer shells and benthic foraminifer assemblages in marine sediments indicate that enhanced CDW upwelling, controlled by the latitudinal position of the Southern Hemisphere westerly winds, forced deglaciation of this sector from at least 10,400 years ago until 7,500 years ago—when an ice-shelf collapse may have caused rapid ice-sheet thinning further upstream—and since the 1940s. These results increase confidence in the predictive capability of current ice-sheet models.

Ice stream bed topography contains key evidence for the ways ice streams interact with, and are potentially controlled by, their beds. Here we present the first application of two‐dimensional Fourier analysis to 22 marine and terrestrial topographies from 5 regions in Antarctica and Canada, with and without mega‐scale glacial lineations (MSGLs). We find that the topography of MSGL‐rich ice stream sedimentary beds is characterized by multiple, periodic wavelengths between 300 and 1200 m and amplitudes from decimeters to a few meters. This periodic topography is consistent with the idea that instability is a key element to the formation of MSGL bedforms. Dominant wavelengths vary among locations and, on one paleo ice stream bed, increase along the direction of ice flow by 1.7 ± 0.52% km−1. We suggest that these changes are likely to reflect pattern evolution via downstream wavelength coarsening, even under potentially steady ice stream geometry and flow conditions. The amplitude of MSGLs is smaller than that of other fluvial and glacial topographies but within the same order of magnitude. However, MSGLs are a striking component of ice stream beds because the topographic amplitude of features not aligned with ice flow is reduced by an order of magnitude relative to those oriented with the flow direction. This study represents the first attempt to automatically derive the spectral signatures of MSGLs. It highlights the plausibility of identifying these landform assemblages using automated techniques and provides a benchmark for numerical models of ice stream flow and subglacial landscape evolution. MSGL topography consists of superimposed periodic wavelengths as expected in the instability theory of subglacial bedform formationMost of the dominant wavelengths present within one extensive MSGL field increase downstream, suggesting that MSGLs evolve via pattern coarseningFor MSGLs to generate periodic topography, sediment must be able to accumulate or erode (freely move) without fixed anchor points

The Transantarctic Mountains (TAMs) separate the warmer lithosphere of the Cretaceous‐Tertiary West Antarctic rift system and the colder and older provinces of East Antarctica. Low velocity zones beneath the TAM imaged in recent seismological studies have been interpreted as warm low‐density mantle material, suggesting a strong contribution of thermal support to the uplift of the TAM. We present new Curie Point Depth (CPD) and geothermal heat flow (GHF) maps of the northern TAM and adjacent Wilkes Subglacial Basin (WSB) based exclusively on high resolution magnetic airborne measurements. We find shallow CPD and high GHF beneath the northern TAM, reinforcing the hypothesis of thermal support of the topography of the mountain range. Additionally, this study demonstrates, that limiting spectral analysis to areas with a high density of aeromagnetic measurements increases the resolution of CPD estimates revealing localized shallow CPD and associated high heat flow in the Central Basin of the WSB and the Rennick Graben (RG). Across the study area the CPD ranges from 15 to 35 km and the GHF values range from 30 to 110 mW/m2. The recovered CPD range is compatible with recent Moho depth estimates, as the CPD predominantly lies within the crust, rather than in the magnetite‐poor mantle. GHF estimates, based on the CPD estimates, show a good agreement to sparse in situ GHF measurements and the location of active volcanoes. Comparison to existing continent‐wide GHF estimates shows strong differences from magnetically‐derived heat flow estimates, while seismologically‐derived heat flow estimates show the best agreement to our results. The Transantarctic Mountains (TAMs) separating the ancient East Antarctic region from the younger West Antarctic region. Recent studies found regions where acoustic waves move slowly beneath the TAM. Slow wave speeds are considered to be caused by warmer rocks at depth, which provide thermal support for the TAM. We present heat flow maps for this region. The estimated heat flow is based on “Curie Point Depths” (CPDs), which describes the depth at which magnetic rocks lose their magnetic properties due to increasing temperature. This temperature is called the Curie temperature (∼580°C for magnetite, the most common magnetic mineral in the Earth's crust). The heat flow at the surface is estimated from the depth at which this temperature is reached. We use for the first time for a CPD study in Antarctica exclusively high‐resolution magnetic data measured from aircraft. We demonstrated that limiting this method to airborne data improves the resolution of CPD variation, and thus heat flow estimates, substantially. We find shallow CPD and high heat flow beneath the TAM, which supports the suggestion that warmer rocks provide thermal support for the TAM. Additionally, we find high heat flow in the adjacent Wilkes Subglacial Basin and Rennick Graben. Spectral analysis of exclusively aeromagnetic data over East Antarctica provides new Curie Point Depth and geothermal heat flow estimatesHigh heat flow imaged beneath the Transantarctic Mountains reinforces the hypothesis that the topography is thermally supportedUsing exclusively airborne data improves the resolution of Curie depth estimation, revealing heat flow correlated to geological features Spectral analysis of exclusively aeromagnetic data over East Antarctica provides new Curie Point Depth and geothermal heat flow estimates High heat flow imaged beneath the Transantarctic Mountains reinforces the hypothesis that the topography is thermally supported Using exclusively airborne data improves the resolution of Curie depth estimation, revealing heat flow correlated to geological features

We present results of a multi-proxy study on marine sediment core JR179-PC466 recovered from the crest of a sediment drift off the West Antarctic Peninsula at approximately 2300 m water depth. The 10.45 m-long core consists dominantly of glaciomarine terrigenous sediments, with only traces of calcium carbonate (<1 wt%). Despite the very low abundance of calcareous foraminifera, planktonic shell numbers are sufficient for stable isotope analyses in two-thirds of the samples studied. The core chronology is based on oxygen isotope stratigraphy and correlation of its relative palaeomagnetic intensity (RPI) with a stacked reference curve. According to the age model, core PC466 spans the last 75 ka, with average sedimentation rates of between about 4 and 25 cm ka−1. Planktonic foraminifera abundances fluctuate between 0 and 30 individuals per gram throughout the core, with minima observed during Marine Isotope Stage (MIS) 2 (14–29 ka before present, BP) and MIS4 (57–71 ka BP). Planktonic foraminifera are present in the Holocene but more abundant in sediments deposited during MIS3 (29–57 ka BP), owing to less dilution by terrigenous detritus and/or better carbonate preservation. During MIS3, foraminifera maxima correlate with Antarctic warming events as recorded in the δ18O signal of the EPICA Dronning Maud Land (EDML) ice core. They indicate higher planktonic foraminifera production and better carbonate preservation west of the Antarctic Peninsula during that time. The abundance of ice-rafted detritus (IRD) in core PC466 increased during the last deglaciation between about 19 and 11 ka BP, when numerous icebergs drifted across the core site, thereby releasing IRD. During this time, sea-level rise destabilized the Antarctic Peninsula (APIS) and West Antarctic (WAIS) ice sheets that had advanced onto the shelf during the sea-level low-stand of the Last Glacial Maximum (LGM; c. 19–23 ka BP). Overall, our results demonstrate that it is possible to establish an age model and reconstruct palaeoceanographical and climatic changes at high temporal resolution from sedimentary sequences recovered at 2300 m water depth from a West Antarctic drift.

Ice shelves restrict outflow from many of the largest glaciers in Antarctica, thus limiting the Antarctic contribution to sea‐level rise. However, past ice‐shelf collapse events show they are highly vulnerable to surface and basal melting. Collapse of ice shelves in front of glaciers flowing on retrograde slopes could initiate runaway retreat processes. Difficulty in projecting how quickly these could play out makes dynamic ice loss from Antarctica the largest uncertainty in predicting future sea‐level rise. Basal melting can impact structural integrity of ice shelves in several ways. Results from analyses of variations in ice‐shelf roughness by Watkins et al. (2021; https://doi.org/10.1029/2021GL094743) raise the tantalizing prospect that this may provide a simple quantitative measure of how the structural integrity of an ice shelf has been impacted by basal melting. Applying the method to additional ice shelves would be useful to examine how other factors may contribute to roughness. In many places around Antarctica, and some around Greenland, glaciers flowing on beds that are hundreds of meters below sea level continue out over the sea for some distance, forming “ice shelves”. Over the past 40 years, break‐up of several ice shelves has been observed and shown to result from melting on their surface and/or at their base. Subsequent increases in the flow speed of glaciers that flowed into them confirmed that ice shelves restrict outflow of glaciers and thus limit their contribution to sea‐level rise. Furthermore, changes resulting from removal of ice shelves could lead to runaway retreat of glaciers flowing on beds that get deeper upstream, and lack of knowledge about whether and how quickly this will happen is the largest uncertainty in predicting future sea‐level rise. Basal melting weakens ice shelves in several ways. Results from analyses of variations in thickness, or “roughness”, of ice shelves by Watkins et al. (2021; https://doi.org/10.1029/2021GL094743) suggest that they may provide a measure of the extent to which basal melting has weakened an ice shelf. Conducting similar analyses on additional ice shelves would help show the extent to which other processes contribute to roughness. Ice shelves restrict glacier flow but are vulnerable to environmental changeBasal melting can impact the structural integrity of ice shelves in several waysA new study suggests that ice‐shelf roughness may provide a measure of the extent to which basal melting has impacted structural integrity Ice shelves restrict glacier flow but are vulnerable to environmental change Basal melting can impact the structural integrity of ice shelves in several ways A new study suggests that ice‐shelf roughness may provide a measure of the extent to which basal melting has impacted structural integrity

Multichannel seismic reflection profiles from the continental rise west of the Antarctic Peninsula between 63° and 69°S show the growth of eight very large mound-shaped sedimentary bodies. MCS profiles and long-range side-scan sonar (GLORIA) images show the sea floor between mounds is traversed by channels originating in a dendritic pattern near the base of the continental slope. The mounds are interpreted as sediment drifts, constructed mainly from the fine-grained components of turbidity currents originating on the continental slope, entrained in a nepheloid layer within the ambient southwesterly bottom currents and redeposited downcurrent.

Despite a spreading rate of 65–70 km Ma−1, the East Scotia Ridge has, along most of its length, a form typically associated with slower rates of sea floor spreading. This may be a consequence of cooler than normal mantle upwelling, which could be a feature of back-arc spreading. At the northern end of the ridge, recently acquired sonar data show a complex, rapidly evolving pattern of extension within 100 km of the South Sandwich Trench. New ridge segments appear to be nucleating at or near the boundary between the South American and Scotia Sea plates and propagating southwards, supplanting older segments. The most prominent of these, north of 56°30′S, has been propagating at a rate of approximately 60 km Ma−1 for at least 1 Ma, and displays a morphology unique on this plate boundary. A 40 km long axial high exists at the centre of this segment, forming one of the shallowest sections of the East Scotia Ridge. Beneath it, seismic reflection profiles reveal an axial magma chamber, or AMC, reflector, similar to those observed beneath the East Pacific Rise and Valu Fa Ridge. Simple calculations indicate the existence here of a narrow (<1 km wide) body of melt at a depth of approximately 3 km beneath the sea floor. From the topographic and seismic data, we deduce that a localised mantle melting anomaly lies beneath this segment. Rates of spreading in the east Scotia Sea show little variation along axis. Hence, the changes in melt supply are related to the unique tectonic setting, in which the South American plate is tearing to the east, perhaps allowing mantle flow around the end of the subducting slab. Volatiles released from the torn plate edge and entrained in the flow are a potential cause of the anomalous melting observed. A southward mantle flow may have existed beneath the axis of the East Scotia Ridge throughout its history.

Despite a spreading rate of 65–70 km Ma−1, the East Scotia Ridge has, along most of its length, a form typically associated with slower rates of sea floor spreading. This may be a consequence of cooler than normal mantle upwelling, which could be a feature of back-arc spreading. At the northern end of the ridge, recently acquired sonar data show a complex, rapidly evolving pattern of extension within 100 km of the South Sandwich Trench. New ridge segments appear to be nucleating at or near the boundary between the South American and Scotia Sea plates and propagating southwards, supplanting older segments. The most prominent of these, north of 56°30′S, has been propagating at a rate of approximately 60 km Ma−1for at least 1 Ma, and displays a morphology unique on this plate boundary. A 40 km long axial high exists at the centre of this segment, forming one of the shallowest sections of the East Scotia Ridge. Beneath it, seismic reflection profiles reveal an axial magma chamber, or AMC, reflector, similar to those observed beneath the East Pacific Rise and Valu Fa Ridge. Simple calculations indicate the existence here of a narrow (<1 km wide) body of melt at a depth of approximately 3 km beneath the sea floor. From the topographic and seismic data, we deduce that a localised mantle melting anomaly lies beneath this segment. Rates of spreading in the east Scotia Sea show little variation along axis. Hence, the changes in melt supply are related to the unique tectonic setting, in which the South American plate is tearing to the east, perhaps allowing mantle flow around the end of the subducting slab. Volatiles released from the torn plate edge and entrained in the flow are a potential cause of the anomalous melting observed. A southward mantle flow may have existed beneath the axis of the East Scotia Ridge throughout its history.