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 40 million years with a coupled Antarctic ice sheet-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. 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
Award ID(s):
1663693
PAR ID:
10283144
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Palaeogeography palaeoclimatology palaeoecology
Volume:
537
ISSN:
0031-0182
Page Range / eLocation ID:
109374
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. 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
  3. Abstract Oscillations in ice sheet extent during early and middle Miocene are intermittently preserved in the sedimentary record from the Antarctic continental shelf, with widespread erosion occurring during major ice sheet advances, and open marine deposition during times of ice sheet retreat. Data from seismic reflection surveys and drill sites from Deep Sea Drilling Project Leg 28 and International Ocean Discovery Program Expedition 374, located across the present-day middle continental shelf of the central Ross Sea (Antarctica), indicate the presence of expanded early to middle Miocene sedimentary sections. These include the Miocene climate optimum (MCO ca. 17–14.6 Ma) and the middle Miocene climate transition (MMCT ca. 14.6–13.9 Ma). Here, we correlate drill core records, wireline logs and reflection seismic data to elucidate the depositional architecture of the continental shelf and reconstruct the evolution and variability of dynamic ice sheets in the Ross Sea during the Miocene. Drill-site data are used to constrain seismic isopach maps that document the evolution of different ice sheets and ice caps which influenced sedimentary processes in the Ross Sea through the early to middle Miocene. In the early Miocene, periods of localized advance of the ice margin are revealed by the formation of thick sediment wedges prograding into the basins. At this time, morainal bank complexes are distinguished along the basin margins suggesting sediment supply derived from marine-terminating glaciers. During the MCO, biosiliceous-bearing sediments are regionally mapped within the depocenters of the major sedimentary basin across the Ross Sea, indicative of widespread open marine deposition with reduced glacimarine influence. At the MMCT, a distinct erosive surface is interpreted as representing large-scale marine-based ice sheet advance over most of the Ross Sea paleo-continental shelf. The regional mapping of the seismic stratigraphic architecture and its correlation to drilling data indicate a regional transition through the Miocene from growth of ice caps and inland ice sheets with marine-terminating margins, to widespread marine-based ice sheets extending across the outer continental shelf in the Ross Sea. 
    more » « less
  4. Ancient glaciated landscapes reveal interactions among ice dynamics, bed composition, and preglacial river networks. Subglacial landscapes, revealed in regions of recent ice-sheet retreat, provide a window into ice-sheet dynamics and interactions with evolving subglacial topography. Here, we document landscape evolution beneath the southern Laurentide Ice Sheet of North America since the end of the Pliocene, 2.6 million years (Ma) ago, by reconstructing the isostatically adjusted preglacial surface and modern bedrock topography at 250 m horizontal resolution. We use flow routing to reconstruct drainage networks and river longitudinal profiles, revealing the pattern and extent of their glacially forced reorganization. The overall mean Quaternary (2.6 Ma ago to present) erosion rate is 27 m/Ma, rising within ice-streaming corridors to 35 m/Ma (and locally reaching 400 m/Ma) and falling to 22 m/Ma in non–ice-streaming regions. Our results suggest that subglacial erosion was sufficient to lower the southern Laurentide Ice Sheet into warmer environments, thereby enhancing ablation and reducing ice-sheet extent over time. 
    more » « less
  5. Bedrock topography is a key boundary condition for ice sheet modeling, and determining changes in subglacial topography through time can provide insight into the timing of ice sheet development, the magnitude of glacial erosion, and the co-development of glaciers and glacial topography. West Antarctica hosts an unusually high geothermal gradient supported by hot, low-viscosity mantle which likely enhanced the lithospheric response to West Antarctic Ice Sheet (WAIS) cycles of growth and increased the sensitivity of thermochronometers to landscape evolution on million-year timescales. Thus, a valuable record of glacial landscape change might be recovered from apatite fission track [AFT 80-130°C range] and (U-Th)/He [AHe; 50-90°C] dating, provided that landscape evolution can be distinguished from tectonic signals, including the effects of faults. This study utilizes AFT-AHe thermochronology and thermo-kinematic Pecube modeling to investigate interactions between the hot geotherm, glacial erosion, and inferred crustal structures in the Ford Ranges and the DeVicq Glacier trough in western and central Marie Byrd Land (MBL), respectively. The Ford Ranges host glacial troughs (up to 3km relief) dissecting a low-relief erosional surface. Previous work suggests a majority of bedrock exhumation and cooling occurred at/by 80 Ma. However, new data hint at renewed exhumation linked to glacial incision since WAIS formation at 34 or 20 Ma. Prior (U-Th)/He zircon dates from exposures of crystalline bedrock span 90 – 67 Ma. New AHe bedrock dates are 41 to 26 Ma, while two glacial erratics (presumed to be eroded from walls or floor of glacial troughs) yielded AHe dates of 37 Ma and 16 Ma. Initial modeling results suggest a tectonic boundary between the Ford Ranges and Edward VII Peninsula separating regions with distinct exhumation histories. The boundary may cause differential WAIS incision at 34 or 20 Ma, a possibility being investigated with new models. The DeVicq Glacier trough (>3.5km relief) coincides with a prominent crustal lineament but lacks temperature-time information compared to other regions. The crustal structure may have accommodated motion between elevated central MBL and the subdued crust of the Ford Ranges. Here, owing to the lack of onshore non-volcanic bedrock exposure, we have employed AHe and AFT dating of glacial sediment marine core samples offshore of the DeVicq Glacier to investigate the timing and rates of exhumation of the bedrock carved by the DeVicq trough, with initial results revealing detrital AHe ages as young as 24 Ma. Our new Pecube models test a series of thermal, tectonic, and landscape evolution scenarios against a suite of thermochronologic data, allowing us to assess the timing of glacial incision and WAIS initiation in the Ford Ranges, and to seek evidence of an inferred tectonic boundary at DeVicq Trough. Modeling efforts will be aided by new AHe and AFT analyses from ongoing work. These models combine topographic, tectonic, thermal, and key thermochronologic datasets to produce new insight into the unique cryosphere-lithosphere interactions affecting landscape change in West Antarctica. 
    more » « less