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.

Application of the Post-Widder Laplace inversion algorithm to postseismic rebound models

SPADA, GIORGIO
2009

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.
Calcolo Scientifico nella Fisica Italiana - CSFI 2008
197
200
Cannelli, V.; Melini, D.; Piersanti, A.; SPADA, GIORGIO
File in questo prodotto:
Eventuali allegati, non sono esposti

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/11585/771913
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 1
  • ???jsp.display-item.citation.isi??? 1
social impact