Astronomy and Astrophysics – Astronomy
Scientific paper
May 2005
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2005geoji.161..421l&link_type=abstract
Geophysical Journal International, Volume 161, Issue 5, pp. 421-444.
Astronomy and Astrophysics
Astronomy
29
Crustal Motions, Finite Volumes, Glacial Isostatic Adjustment, 3-D Viscoelastic Earth Models
Scientific paper
We describe and present results from a finite-volume (FV) parallel computer code for forward modelling the Maxwell viscoelastic response of a 3-D, self-gravitating, elastically compressible Earth to an arbitrary surface load. We implement a conservative, control volume discretization of the governing equations using a tetrahedral grid in Cartesian geometry and a low-order, linear interpolation. The basic starting grid honours all major radial discontinuities in the Preliminary Reference Earth Model (PREM), and the models are permitted arbitrary spatial variations in viscosity and elastic parameters. These variations may be either continuous or discontinuous at a set of grid nodes forming a 3-D surface within the (regional or global) modelling domain. In the second part of the paper, we adopt the FV methodology and a spherically symmetric Earth model to generate a suite of predictions sampling a broad class of glacial isostatic adjustment (GIA) data types (3-D crustal motions, long-wavelength gravity anomalies). These calculations, based on either a simple disc load history or a global Late Pleistocene ice load reconstruction (ICE-3G), are benchmarked against predictions generated using the traditional normal-mode approach to GIA. The detailed comparison provides a guide for future analyses (e.g. what grid resolution is required to obtain a specific accuracy?) and it indicates that discrepancies in predictions of 3-D crustal velocities less than 0.1 mm yr-1 are generally obtainable for global grids with ~3 × 106 nodes; however, grids of higher resolution are required to predict large-amplitude (>1 cm yr-1) radial velocities in zones of peak post-glacial uplift (e.g. James bay) to the same level of absolute accuracy. We conclude the paper with a first application of the new formulation to a 3-D problem. Specifically, we consider the impact of mantle viscosity heterogeneity on predictions of present-day 3-D crustal motions in North America. In these tests, the lateral viscosity variation is constructed, with suitable scaling, from tomographic images of seismic S-wave heterogeneity, and it is characterized by approximately 2 orders of magnitude (peak-to-peak) lateral variations within the lower mantle and 1 order of magnitude variations in the bulk of the upper mantle (below the asthenosphere). We find that the introduction of 3-D viscosity structure has a profound impact on horizontal velocities; indeed, the magnitude of the perturbation (of order 1 mm yr-1) is as large as the prediction generated from the underlying (1-D) radial reference model and it far exceeds observational uncertainties currently being obtained from space-geodetic surveying. The relative impact of lateral viscosity variations on predicted radial motions is significantly smaller.
Christara Christina C.
Komatitsch Dimitri
Latychev Konstantin
Mitrovica Jerry X.
Tamisiea Mark E.
No associations
LandOfFree
Glacial isostatic adjustment on 3-D Earth models: a finite-volume formulation does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Glacial isostatic adjustment on 3-D Earth models: a finite-volume formulation, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Glacial isostatic adjustment on 3-D Earth models: a finite-volume formulation will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1434839