An official website of the United States government
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
Inference of the Timescale‐Dependent Apparent Viscosity Structure in the Upper Mantle Beneath Greenland
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
- Award ID(s):
- 1928146
- PAR ID:
- 10398792
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- AGU Advances
- Volume:
- 4
- Issue:
- 2
- ISSN:
- 2576-604X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY Constitutive laws are a necessary ingredient in calculations of glacial isostatic adjustment (GIA) or other surface loading problems (e.g. loading by ocean tides). An idealized constitutive law governed by the Maxwell viscoelastic model is widely used but increasing attention is being directed towards more intricate constitutive laws that, in particular, include transient rheology. In this context, transient rheology collectively refers to dissipative mechanisms activated in addition to creep modelled by the Maxwell viscoelastic model. Consideration of such viscoelastic models in GIA is in its infancy and to encourage their wider use, we present constitutive laws for several experimentally derived transient rheologies and outline a flexible method in which to incorporate them into geophysical problems, such as the viscoelastic deformation of the Earth induced by surface loading. To further motivate this need, we demonstrate, via the Love number collocation method, how predictions of crustal displacement depart significantly between Earth models that adopt only Maxwell viscoelasticity and those with transient rheology. Throughout this paper, we highlight the differences in terminology and emphases between the rock mechanics, seismology and GIA communities, which have perhaps contributed towards the relative scarcity in integrating this broader and more realistic class of constitutive laws within GIA. We focus on transient rheology since the associated deformation has been demonstrated to operate on timescales that range from hours to decades. With ice mass loss enhanced at similar timescales as a consequence of anthropogenically caused climate change, the ability to model GIA with more accurate constitutive laws is an important tool to investigate such problems.more » « less
-
SUMMARY Earth structure beneath the Antarctic exerts an important control on the evolution of the ice sheet. A range of geological and geophysical data sets indicate that this structure is complex, with the western sector characterized by a lithosphere of thickness ∼50–100 km and viscosities within the upper mantle that vary by 2–3 orders of magnitude. Recent analyses of uplift rates estimated using Global Navigation Satellite System (GNSS) observations have inferred 1-D viscosity profiles below West Antarctica discretized into a small set of layers within the upper mantle using forward modelling of glacial isostatic adjustment (GIA). It remains unclear, however, what these 1-D viscosity models represent in an area with complex 3-D mantle structure, and over what geographic length-scale they are applicable. Here, we explore this issue by repeating the same modelling procedure but applied to synthetic uplift rates computed using a realistic model of 3-D viscoelastic Earth structure inferred from seismic tomographic imaging of the region, a finite volume treatment of GIA that captures this complexity, and a loading history of Antarctic ice mass changes inferred over the period 1992–2017. We find differences of up to an order of magnitude between the best-fitting 1-D inferences and regionally averaged depth profiles through the 3-D viscosity field used to generate the synthetics. Additional calculations suggest that this level of disagreement is not systematically improved if one increases the number of observation sites adopted in the analysis. Moreover, the 1-D models inferred from such a procedure are non-unique, that is a broad range of viscosity profiles fit the synthetic uplift rates equally well as a consequence, in part, of correlations between the viscosity values within each layer. While the uplift rate at each GNSS site is sensitive to a unique subspace of the complex, 3-D viscosity field, additional analyses based on rates from subsets of proximal sites showed no consistent improvement in the level of bias in the 1-D inference. We also conclude that the broad, regional-scale uplift field generated with the 3-D model is poorly represented by a prediction based on the best-fitting 1-D Earth model. Future work analysing GNSS data should be extended to include horizontal rates and move towards inversions for 3-D structure that reflect the intrinsic 3-D resolving power of the data.more » « less
-
Abstract The rheology of the crust and mantle and the interaction of viscoelastic flow with seismic/aseismic slip on faults control the state of stress in the lithosphere over multiple seismic cycles. The rheological behavior of rocks is well constrained in a laboratory setting, but thein situproperties of the lithosphere and its lateral variations remain poorly known. Here, we access the lower‐crustal rheology in Southern California by exploiting 8 years of geodetic postseismic deformation following the 2010 El Mayor‐Cucapah earthquake. The data illuminate viscoelastic flow in the lower crust with lateral variations of effective viscosity correlated with the geological province. We show that a Burgers assembly with dashpots following a nonlinear constitutive law can approximate the temporal evolution of stress and strain rate, indicating the activation of nonlinear transient creep before steady‐state dislocation creep. The transient and background viscosities in the lower crust of the Salton Trough are on the order of ~1018and ~1019 Pa s, respectively, about an order of magnitude lower than those in the surrounding regions. We highlight the importance of transient creep, nonlinear flow laws, and lateral variations of rheological properties to capture the entire history of postseismic relaxation following the El Mayor‐Cucapah earthquake.more » « less
-
Abstract Previous studies of flexure in continental settings assert that the elastic thickness of a multilayer lithosphere is controlled by the mechanically competent layer thicknesses only, following a cubic rule. More specifically, the cubic rule statesTe = (Te13+Te23+ … Ten3)1/3, whereTeis the total elastic thickness, andTeiis the elastic thickness of each competent layer. However, it is not necessarily clear thatTeshould be insensitive to the properties of intermediate weak layers (e.g., a weak lower crust) which may act to decouple the surface load from lower competent layers. To test this idea, we formulate 2D viscoelastic loading models with layered viscosity to compute the fully dynamic, time‐dependent response of a multilayer lithosphere with a weak lower crust. Results show that the flexural response of a multilayer lithosphere to a surface load is initially consistent with the cubic rule. However, this solution is transient because the stress associated with the load cannot be transmitted through the weak lower crust on long timescales. Stress in the mantle lithosphere relaxes with time and eventually does not support the load at all due to the decoupling effect of flow in the weak lower crust. The steady state flexure of a multilayer lithosphere is controlled solely by the mechanically competent upper crust such thatTe = Te1, and this new rule is the major finding of this study. Our new findings now explain small estimates ofTein continental settings with thick mantle lithosphere such as the Northern Tien Shan, which previously, were poorly explained by the cubic rule.more » « less
