Abstract
Abstract. Geophysical time series are sometimes sampled irregularly along the time axis. The situation is particularly frequent in palaeoclimatology. Yet, there is so far no general framework for handling the continuous wavelet transform when the time sampling is irregular. Here we provide such a framework. To this end, we define the scalogram as the continuous-wavelet-transform equivalent of the extended Lomb–Scargle periodogram defined in Part 1 of this study (Lenoir and Crucifix, 2018). The signal being analysed is modelled as the sum of a locally periodic component in the time–frequency plane, a polynomial trend, and a background noise. The mother wavelet adopted here is the Morlet wavelet classically used in geophysical applications. The background noise model is a stationary Gaussian continuous autoregressive-moving-average (CARMA) process, which is more general than the traditional Gaussian white and red noise processes. The scalogram is smoothed by averaging over neighbouring times in order to reduce its variance. The Shannon–Nyquist exclusion zone is however defined as the area corrupted by local aliasing issues. The local amplitude in the time–frequency plane is then estimated with least-squares methods. We also derive an approximate formula linking the squared amplitude and the scalogram. Based on this property, we define a new analysis tool: the weighted smoothed scalogram, which we recommend for most analyses. The estimated signal amplitude also gives access to band and ridge filtering. Finally, we design a test of significance for the weighted smoothed scalogram against the stationary Gaussian CARMA background noise, and provide algorithms for computing confidence levels, either analytically or with Monte Carlo Markov chain methods. All the analysis tools presented in this article are available to the reader in the Python package WAVEPAL.
Reference33 articles.
1. Berger, A. L.: Long-Term Variations of Caloric Insolation Resulting from the Earth's Orbital Elements, Quaternary Res., 9, 139–167, https://doi.org/10.1016/0033-5894(78)90064-9, 1978.
2. Berger, A. L., Loutre, M. F., and Mélice, J. L.: Instability of the atsronomical periods from 1.5 Myr BP to 0.5 Myr AP, Paleoclimates, 2, 239–280, 1998.
3. Brockwell, P. and Davis, R.: Time Series: Theory and Methods, Springer Series in Statistics, Second edn., Springer, New York, USA, 1991.
4. Carmona, R., Hwang, W., and Torresani, B.: Characterization of signals by the ridges of their wavelet transforms, IEEE T. Signal Proces., 45, 2586–2590, https://doi.org/10.1109/78.640725, 1997.
5. Carmona, R., Hwang, W., and Torresani, B.: Multiridge detection and time-frequency reconstruction, IEEE T. Signal Proces., 47, 480–492, https://doi.org/10.1109/78.740131, 1999.
Cited by
6 articles.
订阅此论文施引文献
订阅此论文施引文献，注册后可以免费订阅5篇论文的施引文献，订阅后可以查看论文全部施引文献