An official website of the United States government
Abstract Seismic tomography models indicate highly variable Earth structure beneath Antarctica with anomalously low shallow mantle viscosities below West Antarctica. An improved projection of the contribution of the Antarctic Ice Sheet to sea‐level change requires consideration of this complexity to precisely account for water expelled into the ocean from uplifting marine sectors. Here we build a high‐resolution 3‐D viscoelastic structure model based on recent inferences of seismic velocity heterogeneity below the continent. The model serves as input to a global‐scale sea‐level model that we use to investigate the influence of solid Earth deformation in Antarctica on future global mean sea‐level (GMSL) rise. Our calculations are based on a suite of ice mass projections generated with a range of climate forcings and suggest that water expulsion from the rebounding marine basins contributes 4%–16% and 7%–14% to the projected GMSL change at 2100 and 2500, respectively.
more »
« less
Three-dimensional glacial isostatic adjustment modeling reconciles conflicting geographic trends in North American marine isotope stage 5a relative sea level observations
Glacial isostatic adjustment (GIA) simulations using earth models that vary viscoelastic structure with depth alone cannot simultaneously fit geographic trends in the elevation of marine isotope stage (MIS) 5a relative sea level (RSL) indicators across continental North America and the Caribbean and yield conflicting estimates of global mean sea level (GMSL). We present simulations with a GIA model that incorporates three-dimensional (3-D) variation in North American viscoelastic earth structure constructed by combining high-resolution seismic tomographic imaging with a new method for mapping this imaging into lateral variations in lithospheric thickness and mantle viscosity. We pair this earth model with a global ice history based on updated constraints on ice volume and geometry. The GIA prediction provides the first simultaneous reconciliation of MIS 5a North American and Caribbean RSL highstands and strengthens arguments that MIS 5a peak GMSL reached values close to that of the Last Interglacial. This result highlights the necessity of incorporating realistic 3-D earth structure into GIA predictions with continent-scale RSL data sets.
more »
« less
- Award ID(s):
- 1927326
- PAR ID:
- 10434969
- Date Published:
- Journal Name:
- Geology
- ISSN:
- 0091-7613
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Recent ice-mass loss driven by warming along the Antarctic Peninsula has resulted in rapid changes in uplift rates across the region. Are such events only a function of recent warming? If not, does the Earth response to such events last long enough to be preserved in Holocene records of relative sea level (RSL), and thus have a bearing on global-scale glacial isostatic adjustment (GIA) models (e.g. ICE-6G)? Answering such questions in Antarctica is hindered by the scarcity of RSL reconstructions within the region. Here, a new RSL reconstruction for Antarctica is presented based on beach ridges from Joinville Island on the Antarctic Peninsula. We find that RSL has fallen 4.9 ± 0.58 m over the past 3100 yr, and that the island experienced a significant increase in the rate of RSL fall from 1540 ± 125 cal. (calibrated) yr B.P. to 1320 ± 125 cal. yr B.P. This increase in the rate of RSL fall is likely due to the viscoelastic response of the solid Earth to terrestrial ice-mass loss from the Antarctic Peninsula, similar to the Earth response experienced after ice-mass loss following acceleration of glaciers behind the collapsed Larsen B ice shelf in 2002 C.E. Additionally, slower rates of beach-ridge progradation from 695 ± 190 cal. yr B.P. to 235 ± 175 cal. yr B.P. potentially reflect erosion of beach ridges from a RSL rise induced by a local glacial advance. The rapid response of the Earth to minor ice-mass changes recorded in the RSL record further supports recent assertions of a more responsive Earth to glacial unloading and at time scales relevant for GIA of Holocene and Pleistocene sea levels. Thus, current continental and global GIA models may not accurately capture the ice-mass changes of the Antarctic ice sheets at decadal and centennial time scales.more » « less
-
Abstract. Predictions of future sea-level change and ice-sheet stability rely on accurate reconstructions of sea levels for past warm intervals, such as the mid-Pliocene Warm Period (MPWP; 3.264–3.025 Ma). The magnitude of MPWP glacial cycles and the relative contribution of meltwater sources remain uncertain. We explore this issue by modeling processes of glacial isostatic adjustment for a wide range of possible MPWP ice-sheet melt zones, including North America, Greenland, Eurasia, and West Antarctica, as well as the Wilkes Basin, the Aurora Basin, and the embayment of Prydz Bay in East Antarctica. As a case study, we use a series of ice histories together with a suite of viscoelastic Earth models to predict global changes in sea level from the Marine Isotope Stage (MIS) M2 glacial to the MIS KM3 interglacial. At the Whanganui Basin (New Zealand), a location with stratigraphic constraints on Pliocene glacial–interglacial sea-level amplitude, the calculated local-sea-level (LSL) rise is on average ∼ 15 % lower than the associated change in the global mean sea level (GMSL) in the ice-sheet scenarios explored here. In contrast, the calculated LSL rise over the deglaciation from MIS M2 to MIS KM3 at Enewetak Atoll is systematically larger than the GMSL change by 10 %. While no single LSL observation (field site) can provide a unique constraint on the sources of ice melt observed during this period, combinations of observations have the potential to yield a stronger constraint on GMSL change and to narrow the list of possible sources.more » « less
-
Abstract. Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological recordsand for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice lossin West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders ofmagnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreatand necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelasticresponse of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs toGIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, andhigh-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recentseismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the AmundsenSea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loadingchanges, we find that a grid resolution of ∼ 13 of the radius of the load or higher is required to accurately capture the elasticresponse of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records andfuture ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding lineconverge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2 % for grid resolutions of 3.75 km andhigher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscousdeformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21stcentury. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser gridin this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1130/G51257.1
