Please use this identifier to cite or link to this item:
|Authors: ||Cannelli, V.*|
|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:||Papers Published / Papers in press|
04.03.08. Theory and Models
Files in This Item:
|AGU 2006b.pdf||6.57 MB||Adobe PDF||View/Open
This item is licensed under a Creative Commons License
Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.