Earth-prints repository, logo   DSpace

About DSpace Software
|earth-prints home page | roma library | bologna library | catania library | milano library | napoli library | palermo library
Please use this identifier to cite or link to this item: http://hdl.handle.net/2122/1946

Authors: Cannelli, V.*
Melini, D.*
Piersanti, A.*
Spada, G.*
Title: Application of the Post-Widder Laplace inversion algorithm to postseismic rebound models
Issue Date: 11-Dec-2006
Keywords: postseismic deformation
Abstract: The computation of global postseismic rebound in a spherically symmetric, stratified, self-gravitating Earth with Maxwell viscoelastic rheology can be carried out semi-analytically with a normal-mode approach. The solution scheme usually involves the application of standard propagator techniques to the equivalent problem in the Laplace domain; to recover the temporal dependence a numerical Laplace anti-transform is required. This step involves the solution of the so-called “secular equation”, whose degree increases linearly with the detail of the stratification modeling, and whose coefficients become extremely ill-conditioned for high harmonic orders. As a result, the practically solvable models are limited to a few viscoelastic layers, and are anyway affected by severe numerical instabilities. To overcome these difficulties, alternative approaches have been explored by several authors, ranging from Runge-Kutta purely numerical integration to the evaluation of Laplace antitransform by a numerical discretization of the Bromwich integral. The Post-Widder algorithm allows the estimation of the Laplace antitransform by sampling numerically the transform on the positive real axis. This method, which has been recently applied to the computation of GIA viscoelastic Love numbers, allows to bypass completely the root-finding procedure while preserving at the same time the analytical normal-mode solution form. In this work, we apply the Post-Widder method to the computation of post-seismic rebound models. We perform a series of benchmarks to optimize the algorithm for speed while checking its stability against earlier results.
Appears in Collections:04.01.05. Rheology
Conference materials
04.03.08. Theory and Models
04.07.02. Geodynamics
Papers Published / Papers in press

Files in This Item:

File SizeFormatVisibility
AGU 2006b.pdf6.57 MBAdobe PDFView/Open

This item is licensed under a Creative Commons License
Creative Commons


Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.


Share this record
Del.icio.us

Citeulike

Connotea

Facebook

Stumble it!

reddit


 

Valid XHTML 1.0! ICT Support, development & maintenance are provided by CINECA. Powered on DSpace Software. CINECA