Power spectral estimates using two-dimensional Morlet-fan wavelets with emphasis on the long wavelengths: jackknife errors, bandwidth resolution and orthogonality properties
MetadataShow full item record
We present a method for estimating the errors on local and global wavelet power spectra using the jackknife approach to error estimation, and compare results with jackknifed multitaper (MT) spectrum estimates. We test the methods on both synthetic and real data, the latter being free air gravity over the Congo basin. To satisfy the independence requirement of the jackknife we investigate the orthogonality properties of the two-dimensional Morlet wavelet. Although Morlet wavelets are non-orthogonal, we show that careful selection of parameters can yield approximate orthogonality in space and azimuth. We also find that, when computed via the Fourier transform, the continuous wavelet transform (CWT) contains errors at very long wavelengths due to the discretisation of large-scale wavelets in the Fourier domain. We hence recommend the use of convolution in the space-domain at these scales, even though this is computationally more expensive. Finally, in providing an investigation into the bandwidth resolution of CWT and MT spectra and errors at long wavelengths, we show that the Morlet wavelet is superior in this regard to Slepian tapers. Wavelets with higher bandwidth-resolution deliver smaller spectral error estimates, in contrast to the MT method, where tapers with higher bandwidth-resolution deliver larger errors. This results in the fan-WT having better spectral estimation properties at long wavelengths than Slepian multitapers.
This article has been accepted for publication in Geophysical Journal International ©: 2013 The Authors. Published by Oxford University Press on behalf of The Royal Astronomical Society. All rights reserved.
Showing items related by title, author, creator and subject.
Kirby, Jon; Swain, Christopher (2014)The strength of the cratonic lithosphere has been controversial. On the one hand, many estimates of effective elastic thickness (Te) greatly exceed the crustal thickness, but on the other the great majority of cratonic ...
Reducing the dimensionality of hyperspectral remotely sensed data with applications for maximum likelihood image classificationSantich, Norman Ty (2007)As well as the many benefits associated with the evolution of multispectral sensors into hyperspectral sensors there is also a considerable increase in storage space and the computational load to process the data. ...
A Synthetic Earth Gravity Model Designed Specifically for Testing Regional Gravimetric Geoid Determination AlgorithmsBaran, Ireneusz; Kuhn, Michael; Claessens, Sten; Featherstone, Will; Holmes, S.; Vanicek, P. (2006)A synthetic [simulated] Earth gravity model (SEGM) of the geoid, gravity and topography has been constructed over Australia specifically for validating regional gravimetric geoid determination theories, techniques and ...