skip to main content


Title: Simulating Marine Isotope Stage 7 with a coupled climate–ice sheet model
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
Award ID(s):
1903197
NSF-PAR ID:
10232854
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Climate of the Past
Volume:
16
Issue:
6
ISSN:
1814-9332
Page Range / eLocation ID:
2183 to 2201
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. The penultimate deglaciation (PDG, ∼138–128 thousand years before present, hereafter ka) is the transition fromthe penultimate glacial maximum (PGM)to the Last Interglacial (LIG, ∼129–116 ka).The LIG stands out as one of the warmest interglacials of the last 800 000 years (hereafter kyr),with high-latitude temperature warmer than today and global sea level likely higher by at least 6 m.Considering the transient nature of the Earth system,the LIG climate and ice-sheet evolution were certainly influenced by the changesoccurring during the penultimate deglaciation.It is thus importantto investigate, with coupled atmosphere–ocean general circulation models (AOGCMs),the climate and environmental response to the large changesin boundary conditions(i.e. orbital configuration, atmospheric greenhouse gas concentrations, ice-sheet geometry and associated meltwater fluxes) occurring during the penultimate deglaciation. A deglaciation working group has recently been set up as part of the Paleoclimate Modelling Intercomparison Project (PMIP) phase 4, with a protocolto perform transient simulations of the last deglaciation (19–11 ka; although the protocol covers 26–0 ka).Similar to the last deglaciation, the disintegration of continental ice sheets during the penultimate deglaciation led to significant changesin the oceanic circulation during Heinrich Stadial 11 (∼136–129 ka).However, the two deglaciations bear significant differences in magnitude and temporal evolution of climate and environmental changes. Here, as part of the Past Global Changes (PAGES)-PMIP working group on Quaternary interglacials (QUIGS), we propose a protocol to perform transient simulations of the penultimate deglaciationunder the auspices of PMIP4.This design includes time-varying changes in orbital forcing, greenhouse gas concentrations, continental ice sheets as well as freshwater input from the disintegration ofcontinental ice sheets.This experiment is designed for AOGCMs to assessthe coupled response of the climate system to all forcings.Additional sensitivity experiments are proposed to evaluate the response to each forcing.Finally, a selection of paleo-records representing different parts of the climate system is presented, providing an appropriatebenchmark for upcoming model–data comparisons across the penultimate deglaciation. 
    more » « less
  2. A climate/vegetation model simulates episodic wetter and drier periods at the 21,000-y precession period in eastern North Africa, the Arabian Peninsula, and the Levant over the past 140,000 y. Large orbitally forced wet/dry extremes occur during interglacial time, ∼130 to 80 ka, and conditions between these two extremes prevail during glacial time, ∼70 to 15 ka. Orbital precession causes high seasonality in Northern Hemisphere (NH) insolation at ∼125, 105, and 83 ka, with stronger and northward extended summer monsoon rains in North Africa and the Arabian Peninsula and increased winter rains in the Mediterranean Basin. The combined effects of these two seasonally distinct rainfall regimes increase vegetation and narrow the width of the Saharan–Arabian desert and semidesert zones. During the opposite phase of the precession cycle (∼115, 95, and 73 ka), NH seasonality is low, and decreased summer insolation and increased winter insolation cause monsoon and storm track rains to decrease and the width of the desert zone to increase. During glacial time (∼70 to 15 ka), forcing from large ice sheets and lowered greenhouse gas concentrations combine to increase winter Mediterranean storm track precipitation; the southward retreat of the northern limit of summer monsoon rains is relatively small, thereby limiting the expansion of deserts. The lowered greenhouse gas concentrations cause the near-equatorial zone to cool and reduce convection, causing drier climate with reduced forest cover. At most locations and times, the simulations agree with environmental observations. These changing regional patterns of climate/vegetation could have influenced the dispersal of early humans through expansions and contractions of well-watered corridors.

     
    more » « less
  3. 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
  4. 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
  5. 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