skip to main content


Title: Modeling Northern Hemispheric Ice Sheet Dynamics, Sea Level Change, and Solid Earth Deformation Through the Last Glacial Cycle
Abstract

Retreat or advance of an ice sheet perturbs the Earth's solid surface, rotational vector, and the gravitational field, which in turn feeds back onto the evolution of the ice sheet over a range of timescales. Throughout the last glacial cycle, ice sheets over the Northern Hemisphere have gone through multiple growth and retreat phases, but the dynamics during these phases are not well understood. In this study, we apply a coupled ice sheet‐glacial isostatic adjustment model to simulate the Northern Hemisphere Ice Sheets over the last glacial cycle. We focus on understanding the influence of solid Earth deformation and gravitational field perturbations associated with surface (ice and water) loading changes on the dynamics of terrestrial and marine‐based ice sheets during different phases of the glacial cycle. Our results show that solid Earth deformation enhances glaciation during growth phases and melting during retreat phases in terrestrial regions through ice‐elevation feedback, and gravitational field perturbations have a stabilizing influence on marine‐based ice sheets in regions such as Hudson Bay in North America and Barents and Kara Seas in Eurasia during retreat phases through sea‐level feedback. Our results also indicate that solid Earth deformation influences the relative sensitivity of the North American and Eurasian ice sheets to climate and thus the timing and magnitude of their fluctuations throughout the last glacial cycle.

 
more » « less
Award ID(s):
1664013 1934477 1135417 2035080
NSF-PAR ID:
10450175
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
126
Issue:
4
ISSN:
2169-9003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Retreat and advance of ice sheets perturb the gravitational field, solidsurface and rotation of the Earth, leading to spatially variable sea-levelchanges over a range of timescales O(100−6 years), which in turn feedback onto ice-sheet dynamics. Coupled ice-sheet–sea-level models havebeen developed to capture the interactive processes between ice sheets, sealevel and the solid Earth, but it is computationally challenging to captureshort-term interactions O(100−2 years) precisely within longer O(103−6 years) simulations. The standard forward sea-level modelling algorithmassigns a uniform temporal resolution in the sea-level model, causing aquadratic increase in total CPU time with the total number of input icehistory steps, which increases with either the length or temporal resolutionof the simulation. In this study, we introduce a new “time window”algorithm for 1D pseudo-spectral sea-level models based on the normal modemethod that enables users to define the temporal resolution at which the iceloading history is captured during different time intervals before thecurrent simulation time. Utilizing the time window, we assign a finetemporal resolution O(100−2 years) for the period of ongoing andrecent history of surface ice and ocean loading changes and a coarsertemporal resolution O(103−6 years) for earlier periods in thesimulation. This reduces the total CPU time and memory required per modeltime step while maintaining the precision of the model results. We explorethe sensitivity of sea-level model results to the model temporal resolutionand show how this sensitivity feeds back onto ice-sheet dynamics in coupledmodelling. We apply the new algorithm to simulate sea-level changes inresponse to global ice-sheet evolution over two glacial cycles and the rapidcollapse of marine sectors of the West Antarctic Ice Sheet in the comingcenturies and provide appropriate time window profiles for each application.The time window algorithm reduces the total CPU time by ∼ 50 % in each of these examples and changes the trend of the total CPU timeincrease from quadratic to linear. This improvement would increase withlonger simulations than those considered here. Our algorithm also allows for couplingtime intervals of annual temporal scale for coupled ice-sheet–sea-levelmodelling of regions such as West Antarctica that are characterized byrapid solid Earth response to ice changes due to the thin lithosphere andlow mantle viscosities. 
    more » « less
  2. null (Ed.)
    International Ocean Discovery Program (IODP) Expedition 382, Iceberg Alley and Subantarctic Ice and Ocean Dynamics, investigated the long-term climate history of Antarctica, seeking to understand how polar ice sheets responded to changes in insolation and atmospheric CO2 in the past and how ice sheet evolution influenced global sea level and vice versa. Five sites (U1534–U1538) were drilled east of the Drake Passage: two sites at 53.2°S at the northern edge of the Scotia Sea and three sites at 57.4°–59.4°S in the southern Scotia Sea. We recovered continuously deposited late Neogene sediment to reconstruct the past history and variability in Antarctic Ice Sheet (AIS) mass loss and associated changes in oceanic and atmospheric circulation. The sites from the southern Scotia Sea (Sites U1536–U1538) will be used to study the Neogene flux of icebergs through “Iceberg Alley,” the main pathway along which icebergs calved from the margin of the AIS travel as they move equatorward into the warmer waters of the Antarctic Circumpolar Current (ACC). In particular, sediments from this area will allow us to assess the magnitude of iceberg flux during key times of AIS evolution, including the following: • The middle Miocene glacial intensification of the East Antarctic Ice Sheet, • The mid-Pliocene warm period, • The late Pliocene glacial expansion of the West Antarctic Ice Sheet, • The mid-Pleistocene transition (MPT), and • The “warm interglacials” and glacial terminations of the last 800 ky. We will use the geochemical provenance of iceberg-rafted detritus and other glacially eroded material to determine regional sources of AIS mass loss. We will also address interhemispheric phasing of ice sheet growth and decay, study the distribution and history of land-based versus marine-based ice sheets around the continent over time, and explore the links between AIS variability and global sea level. By comparing north–south variations across the Scotia Sea between the Pirie Basin (Site U1538) and the Dove Basin (Sites U1536 and U1537), Expedition 382 will also deliver critical information on how climate changes in the Southern Ocean affect ocean circulation through the Drake Passage, meridional overturning in the region, water mass production, ocean–atmosphere CO2 transfer by wind-induced upwelling, sea ice variability, bottom water outflow from the Weddell Sea, Antarctic weathering inputs, and changes in oceanic and atmospheric fronts in the vicinity of the ACC. Comparing changes in dust proxy records between the Scotia Sea and Antarctic ice cores will also provide a detailed reconstruction of changes in the Southern Hemisphere westerlies on millennial and orbital timescales for the last 800 ky. Extending the ocean dust record beyond the last 800 ky will help to evaluate dust-climate couplings since the Pliocene, the potential role of dust in iron fertilization and atmospheric CO2 drawdown during glacials, and whether dust input to Antarctica played a role in the MPT. The principal scientific objective of Subantarctic Front Sites U1534 and U1535 at the northern limit of the Scotia Sea is to reconstruct and understand how ocean circulation and intermediate water formation responds to changes in climate with a special focus on the connectivity between the Atlantic and Pacific basins, the “cold water route.” The Subantarctic Front contourite drift, deposited between 400 and 2000 m water depth on the northern flank of an east–west trending trough off the Chilean continental shelf, is ideally situated to monitor millennial- to orbital-scale variability in the export of Antarctic Intermediate Water beneath the Subantarctic Front. During Expedition 382, we recovered continuously deposited sediments from this drift spanning the late Pleistocene (from ~0.78 Ma to recent) and from the late Pliocene (~3.1–2.6 Ma). These sites are expected to yield a wide array of paleoceanographic records that can be used to interpret past changes in the density structure of the Atlantic sector of the Southern Ocean, track migrations of the Subantarctic Front, and give insights into the role and evolution of the cold water route over significant climate episodes, including the following: • The most recent warm interglacials of the late Pleistocene and • The intensification of Northern Hemisphere glaciation. 
    more » « less
  3. null (Ed.)
    International Ocean Discovery Program Expedition 382, Iceberg Alley and Subantarctic Ice and Ocean Dynamics, investigated the long-term climate history of Antarctica, seeking to understand how polar ice sheets responded to changes in insolation and atmospheric CO2 in the past and how ice sheet evolution influenced global sea level and vice versa. Five sites (U1534–U1538) were drilled east of the Drake Passage: two sites at 53.2°S at the northern edge of the Scotia Sea and three sites at 57.4°–59.4°S in the southern Scotia Sea. We recovered continuously deposited late Neogene sediments to reconstruct the past history and variability in Antarctic Ice Sheet (AIS) mass loss and associated changes in oceanic and atmospheric circulation. The sites from the southern Scotia Sea (Sites U1536–U1538) will be used to study the Neogene flux of icebergs through “Iceberg Alley,” the main pathway along which icebergs calved from the margin of the AIS travel as they move equatorward into the warmer waters of the Antarctic Circumpolar Current (ACC). In particular, sediments from this area will allow us to assess the magnitude of iceberg flux during key times of AIS evolution, including the following: • The middle Miocene glacial intensification of the East Antarctic Ice Sheet, • The mid-Pliocene warm period, • The late Pliocene glacial expansion of the West Antarctic Ice Sheet, • The mid-Pleistocene transition (MPT), and • The “warm interglacials” and glacial terminations of the last 800 ky. We will use the geochemical provenance of iceberg-rafted detritus and other glacially eroded material to determine regional sources of AIS mass loss. We will also address interhemispheric phasing of ice sheet growth and decay, study the distribution and history of land-based versus marine-based ice sheets around the continent over time, and explore the links between AIS variability and global sea level. By comparing north–south variations across the Scotia Sea between the Pirie Basin (Site U1538) and the Dove Basin (Sites U1536 and U1537), Expedition 382 will also deliver critical information on how climate changes in the Southern Ocean affect ocean circulation through the Drake Passage, meridional overturning in the region, water mass production, ocean–atmosphere CO2 transfer by wind-induced upwelling, sea ice variability, bottom water outflow from the Weddell Sea, Antarctic weathering inputs, and changes in oceanic and atmospheric fronts in the vicinity of the ACC. Comparing changes in dust proxy records between the Scotia Sea and Antarctic ice cores will also provide a detailed reconstruction of changes in the Southern Hemisphere westerlies on millennial and orbital timescales for the last 800 ky. Extending the ocean dust record beyond the last 800 ky will help to evaluate dust-climate couplings since the Pliocene, the potential role of dust in iron fertilization and atmospheric CO2 drawdown during glacials, and whether dust input to Antarctica played a role in the MPT. The principal scientific objective of Subantarctic Front Sites U1534 and U1535 at the northern limit of the Scotia Sea is to reconstruct and understand how intermediate water formation in the southwest Atlantic responds to changes in connectivity between the Atlantic and Pacific basins, the “cold water route.” The Subantarctic Front contourite drift, deposited between 400 and 2000 m water depth on the northern flank of an east–west trending trough off the Chilean continental shelf, is ideally situated to monitor millennial- to orbital-scale variability in the export of Antarctic Intermediate Water beneath the Subantarctic Front. During Expedition 382, we recovered continuously deposited sediments from this drift spanning the late Pleistocene (from ~0.78 Ma to recent) and from the late Pliocene (~3.1–2.6 Ma). These sites are expected to yield a wide array of paleoceanographic records that can be used to interpret past changes in the density structure of the Atlantic sector of the Southern Ocean, track migrations of the Subantarctic Front, and give insights into the role and evolution of the cold water route over significant climate episodes, including the following: • The most recent warm interglacials of the late Pleistocene and • The intensification of Northern Hemisphere glaciation. 
    more » « less
  4. Abstract The rate and consequences of future high latitude ice sheet retreat remain a major concern given ongoing anthropogenic warming. Here, new precisely dated stalagmite data from NW Iberia provide the first direct, high-resolution records of periods of rapid melting of Northern Hemisphere ice sheets during the penultimate deglaciation. These records reveal the penultimate deglaciation initiated with rapid century-scale meltwater pulses which subsequently trigger abrupt coolings of air temperature in NW Iberia consistent with freshwater-induced AMOC slowdowns. The first of these AMOC slowdowns, 600-year duration, was shorter than Heinrich 1 of the last deglaciation. Although similar insolation forcing initiated the last two deglaciations, the more rapid and sustained rate of freshening in the eastern North Atlantic penultimate deglaciation likely reflects a larger volume of ice stored in the marine-based Eurasian Ice sheet during the penultimate glacial in contrast to the land-based ice sheet on North America as during the last glacial. 
    more » « less
  5. null (Ed.)
    Abstract. It is widely accepted that orbital variations areresponsible for the generation of glacial cycles during the latePleistocene. However, the relative contributions of the orbital forcingcompared to CO2 variations and other feedback mechanisms causing thewaxing and waning of ice sheets have not been fully understood. Testingtheories of ice ages beyond statistical inferences, requires numericalmodeling experiments that capture key features of glacial transitions. Here,we focus on the glacial buildup from Marine Isotope Stage (MIS) 7 to 6covering the period from 240 to 170 ka (ka: thousand years before present). Thistransition from interglacial to glacial conditions includes one of thefastest Pleistocene glaciation–deglaciation events, which occurred during MIS 7e–7d–7c (236–218 ka). Using a newly developed three-dimensional coupledatmosphere–ocean–vegetation–ice sheet model (LOVECLIP), we simulate thetransient evolution of Northern Hemisphere and Southern Hemisphere ice sheets duringthe MIS 7–6 period in response to orbital and greenhouse gas forcing. For arange of model parameters, the simulations capture the evolution of globalice volume well within the range of reconstructions. Over the MIS 7–6period, it is demonstrated that glacial inceptions are more sensitive toorbital variations, whereas terminations from deep glacial conditions needboth orbital and greenhouse gas forcings to work in unison. For someparameter values, the coupled model also exhibits a critical North Americanice sheet configuration, beyond which a stationary-wave–ice-sheettopography feedback can trigger an unabated and unrealistic ice sheetgrowth. The strong parameter sensitivity found in this study originates fromthe fact that delicate mass imbalances, as well as errors, are integratedduring a transient simulation for thousands of years. This poses a generalchallenge for transient coupled climate–ice sheet modeling, with suchcoupled paleo-simulations providing opportunities to constrain suchparameters. 
    more » « less