Options
Application of the Post-Widder Laplace inversion algorithm to postseismic rebound models
Author(s)
Spada, Giorgio
Istituto di Fisica, Università di Urbino "Carlo Bo" - Urbino, Italy
Type
Conference paper
Language
English
Obiettivo Specifico
3.3. Geodinamica e struttura dell'interno della Terra
Status
Published
Conference Name
Calcolo Scientifico nella fisica italiana - CSFI 2008
Issued date
May 2008
Conference Location
Rimini, 27-30 Maggio
Senigallia, 31 Maggio
Senigallia, 31 Maggio
Abstract
The postseismic response of a viscoelastic Earth can be computed
analytically with a normal-mode approach, based on the application of propagator
methods. This framework suffers from many limitations, mostly connected with the
solution of the secular equation, whose degree scales with the number of viscoelastic
layers so that only low-resolution models can be practically solved. Recently, a
viable alternative to the normal-mode approach has been proposed, based on the
Post-Widder inversion formula. This method allows to overcome some of the intrinsic
limitations of the normal-mode approach, so that Earth models with arbitrary
radial resolution can be employed and general linear non-Maxwell rheologies can be
implemented. In this work, we test the robustness of the method against a standard
normal-mode approach in order to optimize computation performance while
ensuring the solution stability. As an application, we address the issue of finding
the minimum number of layers with distinct elastic properties needed to accurately
describe the postseismic relaxation of a realistic Earth model.
analytically with a normal-mode approach, based on the application of propagator
methods. This framework suffers from many limitations, mostly connected with the
solution of the secular equation, whose degree scales with the number of viscoelastic
layers so that only low-resolution models can be practically solved. Recently, a
viable alternative to the normal-mode approach has been proposed, based on the
Post-Widder inversion formula. This method allows to overcome some of the intrinsic
limitations of the normal-mode approach, so that Earth models with arbitrary
radial resolution can be employed and general linear non-Maxwell rheologies can be
implemented. In this work, we test the robustness of the method against a standard
normal-mode approach in order to optimize computation performance while
ensuring the solution stability. As an application, we address the issue of finding
the minimum number of layers with distinct elastic properties needed to accurately
describe the postseismic relaxation of a realistic Earth model.
File(s)
Loading...
Name
Cannelli et al, NCC 2009.pdf
Description
Nuovo Cimento C, Vol. 32 C, N. 2 Marzo-Aprile 2009
Size
86.24 KB
Format
Adobe PDF
Checksum (MD5)
df7e06cb87b2de93eeac2654975c6310