An official website of the United States government
Abstract Contemporary crustal uplift and relative sea level (RSL) change in Greenland is caused by the response of the solid Earth to ongoing and historical ice mass change. Glacial isostatic adjustment (GIA) models, which seek to match patterns of land surface displacement and RSL change, typically employ a linear Maxwell viscoelastic model for the Earth's mantle. In Greenland, however, upper mantle viscosities inferred from ice load changes and other geophysical phenomena occurring over a range of timescales vary by up to two orders of magnitude. Here, we use full‐spectrum rheological models to examine the influence of transient deformation within the Greenland upper mantle, which may account for these differing viscosity estimates. We use observations of shear wave velocity combined with constitutive rheological models to self‐consistently calculate mechanical properties including the apparent upper mantle viscosity and lithosphere thickness across a broad spectrum of frequencies. We find that the contribution of transient behavior is most significant over loading timescales of 102–103 years, which corresponds to the timeframe of ice mass loss over recent centuries. Predicted apparent lithosphere thicknesses are also in good agreement with inferences made across seismic, GIA, and flexural timescales. Our results indicate that full‐spectrum constitutive models that more fully capture broadband mantle relaxation provide a means of reconciling seemingly contradictory estimates of Greenland's upper mantle viscosity and lithosphere thickness made from observations spanning a range of timescales.
more »
« less
Glacial isostatic adjustment reduces past and future Arctic subsea permafrost
Abstract Sea-level rise submerges terrestrial permafrost in the Arctic, turning it into subsea permafrost. Subsea permafrost underlies ~ 1.8 million km2of Arctic continental shelf, with thicknesses in places exceeding 700 m. Sea-level variations over glacial-interglacial cycles control subsea permafrost distribution and thickness, yet no permafrost model has accounted for glacial isostatic adjustment (GIA), which deviates local sea level from the global mean due to changes in ice and ocean loading. Here we incorporate GIA into a pan-Arctic model of subsea permafrost over the last 400,000 years. Including GIA significantly reduces present-day subsea permafrost thickness, chiefly because of hydro-isostatic effects as well as deformation related to Northern Hemisphere ice sheets. Additionally, we extend the simulation 1000 years into the future for emissions scenarios outlined in the Intergovernmental Panel on Climate Change’s sixth assessment report. We find that subsea permafrost is preserved under a low emissions scenario but mostly disappears under a high emissions scenario.
more »
« less
- Award ID(s):
- 1841888
- PAR ID:
- 10500680
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 15
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract As ice sheets load Earth's surface, they produce ice‐marginal depressions which, when filled with meltwater, become proglacial lakes. We include self‐consistently evolving proglacial lakes in a glacial isostatic adjustment (GIA) model and apply it to the Laurentide ice sheet over the last glacial cycle. We find that the locations of modeled lakes and the timing of their disappearance is consistent with the geological record. Lake loads can deflect topography by >10 m, and volumes collectively approach 30–45 cm global mean sea‐level equivalent. GIA increases deglaciation‐phase lake volume up to five‐fold and average along‐ice‐margin depth ≤90 m compared to glaciation‐phase ice volume analogs—differences driven by changes in the position and size of the peripheral bulge. Since ice‐marginal lake depth affects grounding‐line outflow, GIA‐modulated proglacial lake depths could affect ice‐sheet mass loss. Indeed, we find that Laurentide ice‐margin retreat rate sometimes correlates with proglacial lake presence, indicating that proglacial lakes aid glacial collapse.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
-
ABSTRACT Raised shorelines provide important constraints on past sea levels, glacial isostatic adjustment (GIA), and rates and directions of vertical crustal motion. Although most raised shorelines across NW Scotland relate to post‐Last Glacial Maximum (LGM) glacial‐isostatic rebound, many undated shorelines lie above the marine limit established from isolation basins. Here, we present new optically stimulated luminescence (OSL) ages for a raised marine terrace at an elevation of 28 m in Slaggan Bay of NW Scotland. Four OSL ages suggest the feature is pre‐LGM, likely Marine Isotope Stage (MIS) 3. Global mean sea levels (GMSL) during MIS 3 are thought to have been ~40–60 m below present across most of the globe. We use a pair of GIA models to determine what ice sheet and sea‐level scenarios might provide an explanation for these anomalously high sea levels during MIS 3. Our results suggest that in the absence of tectonic activity, such high MIS 3 shorelines across NW Scotland require a MIS 4 ice sheet in Scotland, with postglacial rebound of the crustal depression following its demise during MIS 3 responsible for the elevated shoreline features at that time.more » « less
