skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Continuous Simulations over the Last 34 Million Years with a Coupled Antarctic Ice Sheet and Sediment Model
Much of the knowledge of Antarctic Ice Sheet variations since its inception ~34 Ma derives from marine sediments on the continental shelf, deposited in glacimarine or sub-ice environments by advancing and retreating grounded ice, and observed today by seismic profiling and coring. If coupled ice-sheet and sediment models can simulate these deposits explicitly, direct comparisons with the sediment record would be valuable in linking it to Cenozoic ice and climate history. Here we apply an existing 3-D ice sheet and sediment model to the whole period of late Cenozoic Antarctic evolution. The ice-sheet model uses local parameterizations of grounding-line flux, ice-shelf hydrofracture and ice cliff failure. The sediment model includes quarrying of bedrock, sub-ice transport, and marine deposition. Atmospheric and oceanic forcing is determined by uniform shifts to modern climatology in proportion to records of atmospheric CO2, deep-sea-core d18O, and orbital insolation variations. Initial ice-free and sediment-free bedrock topography is prescribed from the 34 Ma reconstruction of Wilson et al., Palaeo3, 2011, and their estimated rate of tectonic subsidence is applied in West Antarctica. The model is run continuously from 34 Ma to the present, to capture the entire post-Eocene Antarctic landscape evolution and off-shore sediment packages in a single self-consistent simulation. In order to make these long simulations feasible, the model resolution is very coarse, 80 km. However the ice model's use of local parameterizations for fine-scale dynamical processes yields results that are not seriously degraded compared to finer resolutions in short tests. The primary goals are (1) to reproduce major recognized ice-sheet trends and fluctuations from the Eocene to today, and (2) to produce a 3-D model map of modern sediment deposits. "Strata" are tracked by recording times of deposition within the model sediment stacks. Unconformities in these strata occur in the model that can be compared with observed profiles. Initial results are presented, and preliminary overall comparisons are made with observed sediment packages, focusing on sensitivities to climate forcing, quarrying rates, and sediment parameters that stand in for alternate sediment rheologies.  more » « less
Award ID(s):
1664013
PAR ID:
10113501
Author(s) / Creator(s):
Date Published:
Journal Name:
AGU Fall Meeting
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Much of the knowledge of Antarctic Ice Sheet variations since its inception ∼34 Ma derives from marine sediments on the continental shelf, deposited in glacimarine or sub-ice environments by advancing and retreating grounded ice, and observed today by seismic profiling and coring. If coupled ice-sheet and sediment models can simulate these deposits explicitly, direct comparisons with the sediment record would help in linking it to Cenozoic ice and climate history. Here we apply an existing 3-D ice sheet and sediment model to the whole period of late Cenozoic Antarctic evolution. The ice-sheet model uses local parameterizations of grounding-line flux, ice-shelf hydrofracture and ice cliff failure. The sediment model includes quarrying of bedrock, sub-ice transport, and marine deposition. Atmospheric and oceanic forcing is determined by uniform shifts to modern climatology in proportion to records of atmospheric CO2, deep-sea-core d18O, and orbital insolation variations. Initial ice-free bedrock topography can either be prescribed from geologic reconstructions for ∼34 Ma (Wilson et al., Palaeo3, 2011) or deduced in an iterative procedure fitting to observed modern topography and total sediment amounts. The model is run continuously from 40 Ma to the present, capturing post-Eocene Antarctic landscape evolution and off-shore sediment packages in a single self-consistent simulation. In order to make these long simulations feasible, the model resolution is very coarse, 80 km. However the ice model’s use of local parameterizations for fine-scale dynamical processes yields results that are not seriously degraded compared to finer resolutions in short tests. The primary goals are (1) to reproduce major recognized ice-sheet trends and fluctuations from the Eocene to today, and (2) to produce a 3-D model map of modern sediment deposits. "Strata" are tracked by recording times of deposition within the model sediment stacks, which can be compared with observed seismic profiles. Initial results are presented, and preliminary overall comparisons are made with observed sediment packages and the modern ice and bedrock state. 
    more » « less
  2. null (Ed.)
    Much of the knowledge of Antarctic Ice Sheet variations since its inception ~34 Ma derives from marine sediments on the continental shelf, deposited in glacimarine or sub-ice environments by advancing and retreating grounded ice, and observed today by seismic profiling and coring. Here we apply a 3-D coupled ice sheet and sediment model from 40 Ma to the present, with the goal of directly linking ice-sheet variations with the sediment record. The ice-sheet model uses vertically averaged ice dynamics and parameterized grounding-line flux. The sediment model includes quarrying of bedrock, sub-ice transport, and marine deposition. Atmospheric and oceanic forcing are determined by uniform shifts to modern climatology in proportion to records of atmospheric CO2, deep-sea-core δ18O, and orbital insolation variations. The model is run continuously over the last 40 Myr at coarse resolution (80 or 160 km), modeling post-Eocene ice, landscape evolution and off-shore sediment packages in a single self-consistent simulation. Strata and unconformities are tracked by recording times of deposition within the model sediment stacks, which can be compared directly with observed seismic profiles. The initial bedrock topography is initialized to 34 Ma geologic reconstructions, or an iterative procedure is used that yields independent estimates of paleo bedrock topography. Preliminary results are compared with recognized Cenozoic ice-sheet variations, modern sediment distributions and seismic profiles, and modern and paleo bedrock topographies. 
    more » « less
  3. Abstract Drill cores from the Antarctic continental shelf are essential for directly constraining changes in past Antarctic Ice Sheet extent. Here, we provide a sedimentary facies analysis of drill cores from International Ocean Discovery Program (IODP) Site U1521 in the Ross Sea, which reveals a unique, detailed snapshot of Antarctic Ice Sheet evolution between ca. 18 Ma and 13 Ma. We identify distinct depositional packages, each of which contains facies successions that are reflective of past baseline shifts in the presence or absence of marine-terminating ice sheets on the outermost Ross Sea continental shelf. The oldest depositional package (>18 Ma) contains massive diamictites stacked through aggradation and deposited in a deep, actively subsiding basin that restricted marine ice sheet expansion on the outer continental shelf. A slowdown in tectonic subsidence after 17.8 Ma led to the deposition of progradational massive diamictites with thin mudstone beds/laminae, as several large marine-based ice sheet advances expanded onto the mid- to outer continental shelf between 17.8 Ma and 17.4 Ma. Between 17.2 Ma and 15.95 Ma, packages of interbedded diamictite and diatom-rich mudstone were deposited during a phase of highly variable Antarctic Ice Sheet extent and volume. This included periods of Antarctic Ice Sheet advance near the outer shelf during the early Miocene Climate Optimum (MCO)—despite this being a well-known period of peak global warmth between ca. 17.0 Ma and 14.6 Ma. Conversely, there were periods of peak warmth within the MCO during which diatom-rich mudstones with little to no ice-rafted debris were deposited, which indicates that the Antarctic Ice Sheet was greatly reduced in extent and had retreated to a smaller terrestrial-terminating ice sheet, most notably between 16.3 Ma and 15.95 Ma. Post-14.2 Ma, diamictites and diatomites contain unambiguous evidence of subglacial shearing in the core and provide the first direct, well-dated evidence of highly erosive marine ice sheets on the outermost continental shelf during the onset of the Middle Miocene Climate Transition (MMCT; 14.2–13.6 Ma). Although global climate forcings and feedbacks influenced Antarctic Ice Sheet advances and retreats during the MCO and MMCT, we propose that this response was nonlinear and heavily influenced by regional feedbacks related to the shoaling of the continental shelf due to reduced subsidence, sediment infilling, and local sea-level changes that directly influenced oceanic influences on melting at the Antarctic Ice Sheet margin. Although intervals of diatom-rich muds and diatomite indicating open-marine interglacial conditions still occurred during (and following) the MMCT, repeated advances of marine-based ice sheets since that time have resulted in widespread erosion and overdeepening in the inner Ross Sea, which has greatly enhanced sensitivity to marine ice sheet instability since 14.2 Ma. 
    more » « less
  4. Abstract Mass loss from polar ice sheets is poorly constrained in estimates of future global sea-level rise. Today, the marine-based West Antarctic Ice Sheet is losing mass at an accelerating rate, most notably in the Thwaites and Pine Island glacier drainage basins. Early Pliocene surface temperatures were about 4 °C warmer than preindustrial and maximum sea level stood ~20 m above present. Using data from a sediment archive on the Amundsen Sea continental rise, we investigate the impact of prolonged Pliocene ocean warmth on the ice-sheet−ocean system. We show that, in contrast to today, during peak ocean warming ~4.6 − 4.5 Ma, terrigenous muds accumulated rapidly under a weak bottom current regime after spill-over of dense shelf water with high suspended load down to the rise. From sediment provenance data we infer major retreat of the Thwaites Glacier system at ~4.4 Ma several hundreds of km inland from its present grounding line position, highlighting the potential for major Earth System changes under prolonged future warming. 
    more » « less
  5. Abstract The Antarctic ice sheet blankets >99% of the continent and limits our ability to study how subglacial geology and topography have evolved through time. Ice-rafted dropstones derived from the Antarctic subglacial continental interior at different times during the late Cenozoic provide valuable thermal history proxies to understand this geologic history. We applied multiple thermochronometers covering a range of closure temperatures (60–800 °C) to 10 dropstones collected during Integrated Ocean Drilling Program (IODP) Expedition 318 in order to explore the subglacial geology and thermal and exhumation history of the Wilkes Subglacial Basin. The Wilkes Subglacial Basin is a key target for study because ice-sheet models show it was an area of ice-sheet retreat that significantly contributed to sea-level rise during past warm periods. Depositional ages of dropstones range from early Oligocene to late Pleistocene and have zircon U-Pb or 40Ar/39Ar ages indicating sources from the Mertz shear zone, Adélie craton, Ferrar large igneous province, and Millen schist belt. Dropstones from the Mertz shear zone and Adélie craton experienced three cooling periods (1700–1500 Ma; 500–280 Ma; 34–0 Ma) and two periods of extremely slow cooling rates (1500–500 Ma; 280–34 Ma). Low-temperature thermochronometers from seven of the dropstones record cooling during the Paleozoic, potentially recording the Ross or Pan-African orogenies, and during the Mesozoic, potentially recording late Paleozoic to Mesozoic rifting. These dropstones then resided within ~500 m of the surface since the late Paleozoic and early Mesozoic. In contrast, two dropstones deposited during the mid-Pliocene, one from the Mertz shear zone and one from Adélie craton, show evidence for localized post-Eocene glacial erosion of ≥2 km. 
    more » « less