skip to main content


Search for: All records

Award ID contains: 1928146

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. 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
  2. Abstract

    The landscape hidden beneath the Greenland Ice Sheet remains one of the most sparsely mapped regions on Earth, but offers a unique record of environmental conditions prior to and during widespread glaciation, and of the ice sheet's response to changing climates. In particular, subglacial valleys observed across Greenland may preserve geomorphological information pertaining to landscape and ice sheet evolution. Here we analyze the morphology of a subglacial valley network in northern Greenland using bed elevation measurements acquired during multi‐year airborne radio‐echo sounding surveys. Channel profile morphologies are consistent with a primarily fluvial origin of the network, with evidence for localized modification by ice and/or meltwater. Gravity and magnetic anomalies suggest that the spatial organisation of the valley network is influenced by regional‐scale geological structure, implying a long‐lived and well‐established hydrological system. We also document two knickzones in the valley longitudinal profile and terraces above the channel floor in the lower course of the network. These observations, combined with stream power modeling, indicate that northern Greenland experienced two episodes of relative base level fall during the Neogene (∼150 m at ca. 12–3.7 Ma and ∼380 m at ca. 8.2–2.8 Ma) that resulted in channel profile adjustment via incision and knickzone retreat. The timing of the inferred base level fall correlates with other onshore and offshore records of uplift, denudation, and/or relative sea level change, and we suggest that tectonic and/or mantle‐driven uplift played an important role in the genesis of the modern landscape of northern Greenland.

     
    more » « less
  3. 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
  4. null (Ed.)