skip to main content

Explore Research Products in the NSF-PAR

Title: Sea-level trends across The Bahamas constrain peak last interglacial ice melt

During the last interglacial (LIG) period, global mean sea level (GMSL) was higher than at present, likely driven by greater high-latitude insolation. Past sea-level estimates require elevation measurements and age determination of marine sediments that formed at or near sea level, and those elevations must be corrected for glacial isostatic adjustment (GIA). However, this GIA correction is subject to uncertainties in the GIA model inputs, namely, Earth’s rheology and past ice history, which reduces precision and accuracy in estimates of past GMSL. To better constrain the GIA process, we compare our data and existing LIG sea-level data across the Bahamian archipelago with a suite of 576 GIA model predictions. We calculated weights for each GIA model based on how well the model fits spatial trends in the regional sea-level data and then used the weighted GIA corrections to revise estimates of GMSL during the LIG. During the LIG, we find a 95% probability that global sea level peaked at least 1.2 m higher than today, and it is very unlikely (5% probability) to have exceeded 5.3 m. Estimates increase by up to 30% (decrease by up to 20%) for portions of melt that originate from the Greenland ice sheet (West Antarctic ice sheet). Altogether, this work suggests that LIG GMSL may be lower than previously assumed.

  more » « less
Award ID(s):
1841888
NSF-PAR ID:
10285778
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
Proceedings of the National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
118
Issue:
33
ISSN:
0027-8424
Page Range / eLocation ID:
Article No. e2026839118
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ( , Geophysical Journal International)
    null (Ed.)
    It is generally agreed that the Last Interglacial (LIG; ∼130 – 115 ka) was a time when global average temperatures and global mean sea level were higher than they are today. However, the exact timing, magnitude, and spatial pattern of ice melt is much debated. One difficulty in extracting past global mean sea level from local observations is that their elevations need to be corrected for glacial isostatic adjustment (GIA), which requires knowledge of Earth’s internal viscoelastic structure. While this structure is generally assumed to be radially symmetric, evidence from seismology, geodynamics, and mineral physics indicates that large lateral variations in viscosity exist within the mantle. In this study, we construct a new model of Earth’s internal structure by converting shear wave speed into viscosity using parameterisations from mineral physics experiments and geodynamical constraints on Earth’s thermal structure. We use this 3D Earth structure, which includes both variations in lithospheric thickness and lateral variations in viscosity, to calculate the first 3D GIA prediction for LIG sea level. We find that the difference between predictions with and without lateral Earth structure can be meters to 10s of meters in the near field of former ice sheets, and up to a few meters in their far field. We demonstrate how forebulge dynamics and continental levering are affected by laterally varying Earth structure, with a particular focus on those sites with prominent LIG sea level records. Results from four 3D GIA calculations show that accounting for lateral structure can act to increase local sea level by up to ∼1.5 m at the Seychelles and minimally decrease it in Western Australia. We acknowledge that this result is only based on a few simulations, but if robust, this shift brings estimates of global mean sea level from these two sites into closer agreement with each other. We further demonstrate that simulations with a suitable radial viscosity profile can be used to locally approximate the 3D GIA result, but that these radial profiles cannot be found by simply averaging viscosity below the sea level indicator site. 
    more » « less
  2. ; ; ; ; ; ( , Journal of Geophysical Research: Earth Surface)

    Previous studies have interpreted Last Interglacial (LIG;∼129–116 ka) sea‐level estimates in multiple different ways to calibrate projections of future Antarctic ice‐sheet (AIS) mass loss and associated sea‐level rise. This study systematically explores the extent to which LIG constraints could inform future Antarctic contributions to sea‐level rise. We develop a Gaussian process emulator of an ice‐sheet model to produce continuous probabilistic projections of Antarctic sea‐level contributions over the LIG and a future high‐emissions scenario. We use a Bayesian approach conditioning emulator projections on a set of LIG constraints to find associated likelihoods of model parameterizations. LIG estimates inform both the probability of past and future ice‐sheet instabilities and projections of future sea‐level rise through 2150. Although best‐available LIG estimates do not meaningfully constrain Antarctic mass loss projections or physical processes until 2060, they become increasingly informative over the next 130 years. Uncertainties of up to 50 cm remain in future projections even if LIG Antarctic mass loss is precisely known (±5 cm), indicating that there is a limit to how informative the LIG could be for ice‐sheet model future projections. The efficacy of LIG constraints on Antarctic mass loss also depends on assumptions about the Greenland ice sheet and LIG sea‐level chronology. However, improved field measurements and understanding of LIG sea levels still have potential to improve future sea‐level projections, highlighting the importance of continued observational efforts.

     
    more » « less
  3. ; ; ; ; ; ; ; ( , Scientific Reports)
    Abstract

    Sea-level reconstructions are important for understanding past ice sheet variability and its response to past and future warming. Here we present Neogene and Quaternary sea-level snapshots using phreatic overgrowths on speleothems (POS) from caves on Mallorca, Spain. POS are excellent sea level index points because of their clear relationship to sea level and precise U–Pb chronology. We find that local sea-level before and at the onset of the Messinian Salinity Crisis was at 33.3 ± 0.25 m (6.54 ± 0.37 Ma) and 31.8 ± 0.25 m (5.86 ± 0.60 Ma) above present levels, respectively. We further present global mean sea level (GMSL) estimates, i.e. local sea level corrected for glacial isostatic adjustment and long-term uplift, for three other POS. The results show that GMSL during the Pliocene–Pleistocene Transition was 6.4 m (− 2.0–8.8 m) at 2.63 ± 0.11 Ma and during the beginning and the end of the Mid-Pleistocene Transition was − 1.1 m (− 5.6–2.4 m) and 5 m (1.5–8.1 m), respectively. These estimates provide important constraints for the past evolution of sea level and show that local sea level prior to the MSC was similar to the highest stand during the Pliocene, with markedly lower position afterwards.

     
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al ( , Geoscientific Model Development)
    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
  5. ; ; ; ; ( , Geophysical Journal International)
    SUMMARY

    Fossil corals are commonly used to reconstruct Last Interglacial (∼125 ka, LIG) sea level. Sea level reconstructions assume the water depth at which the coral lived, called the ‘relative water depth’. However, relative water depth varies in time and space due to coral reef growth in response to relative sea level (RSL) changes. RSL changes can also erode coral reefs, exposing older reef surfaces with different relative water depths. We use a simplified numerical model of coral evolution to investigate how sea level history systematically influences the preservation of corals in the Bahamas and western Australia, regions which house >100 LIG coral fossils. We construct global ice histories spanning the uncertainty of LIG global mean sea level (GMSL) and predict RSL with a glacial isostatic adjustment model. We then simulate coral evolution since 132 ka. We show that preserved elevations and relative water depths of modelled LIG corals are sensitive to the magnitude, timing and number of GMSL highstand(s). In our simulations, the influence of coral growth and erosion (i.e. the ‘growth effect’) can have an impact on RSL reconstructions that is comparable to glacial isostatic adjustment. Thus, without explicitly accounting for the growth effect, additional uncertainty is introduced into sea level reconstructions. Our results suggest the growth effect is most pronounced in western Australia due to Holocene erosion, but also plays a role in the Bahamas, where LIG RSL rose rapidly due to the collapsing peripheral bulge associated with Laurentide Ice Sheet retreat. Despite the coral model's simplicity, our study highlights the utility of process-based RSL reconstructions.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email