skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

Title: Retreat and Regrowth of the Greenland Ice Sheet During the Last Interglacial as Simulated by the CESM2‐CISM2 Coupled Climate–Ice Sheet Model

During the Last Interglacial, approximately 129 to 116 ka (thousand years ago), the Arctic summer climate was warmer than the present, and the Greenland Ice Sheet retreated to a smaller extent than its current state. Previous model‐derived and geological reconstruction estimates of the sea‐level contribution of the Greenland Ice Sheet during the Last Interglacial vary widely. Here, we conduct a transient climate simulation from 127 to 119 ka using the Community Earth System Model (CESM2), which includes a dynamic ice sheet component (the Community Ice Sheet Model, CISM2) that is interactively coupled to the atmosphere, land, ocean, and sea ice components. Vegetation distribution is updated every 500 years based on biomes simulated using a monthly climatology to force the BIOME4 equilibrium vegetation model. Results show a substantial retreat of the Greenland Ice Sheet, reaching a minimum extent at 121.9 ka, equivalent to a 3.0 m rise in sea level relative to the present day, followed by gradual regrowth. In contrast, a companion simulation employing static vegetation based on pre‐industrial conditions shows a much smaller ice‐sheet retreat, highlighting the importance of the changes in high‐latitude vegetation distribution for amplifying the ice‐sheet response.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Paleoceanography and Paleoclimatology
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Sometime during the middle to late Holocene (8.2 ka to ∼ 1850–1900 CE), the Greenland Ice Sheet (GrIS) was smaller than its currentconfiguration. Determining the exact dimensions of the Holocene ice-sheetminimum and the duration that the ice margin rested inboard of its currentposition remains challenging. Contemporary retreat of the GrIS from itshistorical maximum extent in southwestern Greenland is exposing a landscapethat holds clues regarding the configuration and timing of past ice-sheetminima. To quantify the duration of the time the GrIS margin was near itsmodern extent we develop a new technique for Greenland that utilizes in situcosmogenic 10Be–14C–26Al in bedrock samples that have becomeice-free only in the last few decades due to the retreating ice-sheet margin atKangiata Nunaata Sermia (n=12 sites, 36 measurements; KNS), southwest Greenland. To maximizethe utility of this approach, we refine the deglaciation history of the regionwith stand-alone 10Be measurements (n=49) and traditional 14C agesfrom sedimentary deposits contained in proglacial–threshold lakes. We combineour reconstructed ice-margin history in the KNS region with additionalgeologic records from southwestern Greenland and recent model simulations ofGrIS change to constrain the timing of the GrIS minimum in southwestGreenland and the magnitude of Holocene inland GrIS retreat, as well as to explore theregional climate history influencing Holocene ice-sheet behavior. Our10Be–14C–26Al measurements reveal that (1) KNS retreated behindits modern margin just before 10 ka, but it likely stabilized near thepresent GrIS margin for several thousand years before retreating fartherinland, and (2) pre-Holocene 10Be detected in several of our sample sitesis most easily explained by several thousand years of surface exposure duringthe last interglaciation. Moreover, our new results indicate that the minimumextent of the GrIS likely occurred after ∼5 ka, and the GrISmargin may have approached its eventual historical maximum extent as early as∼2 ka. Recent simulations of GrIS change are able to match thegeologic record of ice-sheet change in regions dominated by surface massbalance, but they produce a poorer model–data fit in areas influenced by oceanicand dynamic processes. Simulations that achieve the best model–data fitsuggest that inland retreat of the ice margin driven by early to middleHolocene warmth may have been mitigated by increased precipitation. Triple10Be–14C–26Al measurements in recently deglaciated bedrockprovide a new tool to help decipher the duration of smaller-than-present iceover multiple timescales. Modern retreat of the GrIS margin in southwestGreenland is revealing a bedrock landscape that was also exposed during themigration of the GrIS margin towards its Holocene minimum extent, but it has yetto tap into a landscape that remained ice-covered throughout the entireHolocene. 
    more » « less
  2. This dataset contains the atmospheric river catalogues and the associated precipitation and temperature data for the Preindustrial and Last Glacial Maximum CESM2 simulations presented in the GRL manuscript:  Atmospheric river contributions to ice sheet hydro climate at the Last Glacial Maximum. The atmospheric river catalogue files (zipped) are in netcdf format and organized by year. There are 100 years of data for both simulations.  The Preindustrial simulation catalogue begins in model year 41 and ends in model year 140.  The LGM simulation catalogue begins in model year 1 and ends in year 100. Each yearly file has a temporal resolution of 6 hours (1460 time steps each file) and a spatial resolution of 0.9° x 1.25° (the native resolution of the CESM simulation). A variable in the file called "ar_binary_tag" indicates whether an atmospheric river is present at each grid cell and each tilmestep: 1 indicates an atmospheric river is present; 0 indicates an atmospheric river is not present.  The precipitation and temperature files are 100-year annual or 100-year seasonal averages of atmospheric river precipitation/temperature. See the Methods section of the article for more details on the atmospheric river detection algorithm and precipitation/temperature calculations.

    Associated article abstract:

    Atmospheric rivers (ARs) are an important driver of surface mass balance over today’s Greenland and Antarctic ice sheets. Using paleoclimate simulations with the Community Earth System Model, we find ARs also had a key influence on the extensive ice sheets of the Last Glacial Maximum (LGM). ARs provide up to 53% of total precipitation along the margins of the eastern Laurentide ice sheet and up to 22-27% of precipitation along the margins of the Patagonian, western Cordilleran, and western Fennoscandian ice sheets. Despite overall cold conditions at the LGM, surface temperatures during AR events are often above freezing, resulting in more rain than snow along ice sheet margins and conditions that promote surface melt. The results suggest  ARs may have had an important role in ice sheet growth and melt during previous glacial periods and may have accelerated ice sheet retreat following the LGM.

    more » « less
  3. Accurate reconstruction of Laurentide Ice Sheet volume changes following the Last Glacial Maximum is critical for understanding ice sheet contribution to sea-level rise, the resulting influence of meltwater on oceanic circulation, and the spatial and temporal patterns of deglaciation. Here, we provide empirical constraints on Laurentide Ice Sheet thinning during the last deglaciation by measuring in situ cosmogenic 10Be in 81 samples collected along vertical transects of nine mountains in the northeastern United States. In conjunction with 107 exposure age samples over five vertical transects from previous studies, we reconstruct ice sheet thinning history. At peripheral sites (within 200 km of the terminal moraine), we find evidence for ∼600 m of thinning between 19.5 ka and 17.5 ka, which is coincident with the slow initial margin retreat indicated by varve records. At locations >400 km north of the terminal moraine, exposure ages above and below 1200 m a.s.l. exhibit different patterns. Ages above this elevation are variable and older, while lower elevation ages are indistinguishable over 800−1000 m elevation ranges, a pattern that suggests a subglacial thermal boundary at ∼1200 m a.s.l. separating erosive, warm-based ice below and polythermal, minimally erosive ice above. Low-elevation ages from up-ice mountains are between 15 ka and 13 ka, which suggests rapid thinning of ∼1000 m coincident with Bølling-Allerød warming. These rates of rapid paleo-ice thinning are comparable to those of other vertical exposure age transects around the world and may have been faster than modern basin-wide thinning rates in Antarctica and Greenland, which suggests that the southeastern Laurentide Ice Sheet was highly sensitive to a warming climate. 
    more » « less
  4. null (Ed.)
    The West Antarctic Ice Sheet (WAIS) is largely marine based and thus highly sensitive to both climatic and oceanographic changes. Therefore, the WAIS has likely had a very dynamic history over the last several million years. A complete collapse of the WAIS would result in a global sea level rise of 3.3–4.3 m, yet the world’s scientific community is not able to predict its future behavior. Moreover, knowledge about past behavior of the WAIS is poor, in particular during geological times with climatic conditions similar to those expected for the near and distant future. Reconstructions and quantifications of partial or complete WAIS collapses in the past are urgently needed for constraining and testing ice sheet models that aim to predict future WAIS behavior and the potential contribution of the WAIS to global sea level rise. Large uncertainties exist regarding the chronology, extent, rates, and spatial and temporal variability of past advances and retreats of the WAIS across the continental shelves. These uncertainties largely result from the fundamental lack of data from drill cores recovered proximal to the WAIS. The continental shelf and rise of the Amundsen Sea are prime targets for drilling because the records are expected to yield archives of pure WAIS dynamics unaffected by other ice sheets and the WAIS sector draining into the Amundsen Sea Embayment (ASE) currently experiences the largest ice loss in Antarctica (Paolo et al., 2015). We propose a series of drill sites for the ASE shelf where seismic data reveal seaward-dipping sedimentary sequences that span from the preglacial depositional phase to the most recent glacial periods. Our strategy is to drill a transect from the oldest sequences close to the bedrock/basin boundary at the middle–inner shelf transition to the youngest sequences on the outer shelf in the eastern ASE. If the eastern ASE is inaccessible due to sea ice cover, a similar transect of sites can be drilled on the western ASE. The core transect will provide a detailed history of the glacial cycles in the Amundsen Sea region and allow comparison to the glacial history from the Ross Sea sector. In addition, deep-water sites on the continental rise of the Amundsen Sea are selected for recovering continuous records of glacially transported sediments and detailed archives of climatic and oceanographic changes throughout glacial–interglacial cycles. We will apply a broad suite of analytical techniques, including multiproxy analyses, to address our objectives of reconstructing the onset of glaciation in the greenhouse to icehouse transition, processes of dynamic ice sheet behavior during the Neogene and Quaternary, and ocean conditions associated with the glacial cycles. The five principal objectives of Expedition 379 are as follows: 1. To reconstruct the glacial history of West Antarctica from the Paleogene to recent times and the dynamic behavior of the WAIS during the Neogene and Quaternary, especially possible partial or full WAIS collapses, and the WAIS contribution to past sea level changes. Emphasis is placed in particular on studying the response of the WAIS at times when the pCO2 in Earth’s atmosphere exceeded 400 ppm and atmospheric and oceanic temperatures were higher than at present. 2. To correlate the WAIS-proximal records of ice sheet dynamics in the Amundsen Sea with global records of ice volume changes and proxy records for air and seawater temperatures. 3. To study the relationship between incursions of warm Circumpolar Deep Water (CDW) onto the continental shelf of the Amundsen Sea Embayment and the stability of marine-based ice sheet margins under warm water conditions. 4. To reconstruct the processes of major WAIS advances onto the middle and outer shelf that are likely to have occurred since the middle Miocene and compare their timing and processes to those of other Antarctic continental shelves. 5. To identify the timing of the first ice sheet expansion onto the continental shelf of the ASE and its possible relationship to the uplift of Marie Byrd Land. 
    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