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
-
-
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
-
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
-
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 The land surface beneath the Greenland and Antarctic Ice Sheets is isostatically suppressed by the mass of the overlying ice. Accurate computation of the land elevation in the absence of ice is important when considering, for example, regional geodynamics, geomorphology, and ice sheet behaviour. Here, we use contemporary compilations of ice thickness and lithospheric effective elastic thickness to calculate the fully re-equilibrated isostatic response of the solid Earth to the complete removal of the Greenland and Antarctic Ice Sheets. We use an elastic plate flexure model to compute the isostatic response to the unloading of the modern ice sheet loads, and a self-gravitating viscoelastic Earth model to make an adjustment for the remaining isostatic disequilibrium driven by ice mass loss since the Last Glacial Maximum. Feedbacks arising from water loading in areas situated below sea level after ice sheet removal are also taken into account. In addition, we quantify the uncertainties in the total isostatic response associated with a range of elastic and viscoelastic Earth properties. We find that the maximum change in bed elevation following full re-equilibration occurs over the centre of the landmasses and is +783 m in Greenland and +936 m in Antarctica. By contrast, areas around the ice margins experience up to 123 m of lowering due to a combination of sea level rise, peripheral bulge collapse, and water loading. The computed isostatic response fields are openly accessible and have a number of applications for studying regional geodynamics, landscape evolution, cryosphere dynamics, and relative sea level change.more » « less
