Estimation of the significance of the Foster's wavelet spectrum by means of the permutation test, and applications to paleoclimate records

Here we propose a permutation test-a non-parametric computing-intensive test-to evaluate the statistical significance of the Foster s wavelet spectrum by means of Monte Carlo simulations. A procedure (algorithm) is introduced in order to carry out this aim. The Foster s wavelet is an adequate method...

Full description

Bibliographic Details
Published in:BOLETÍN GEOLÓGICO Y MINERO
Main Authors: Polanco-Martínez, J.M., Faria, S.H.
Format: Article in Journal/Newspaper
Language:Spanish
Published: Instituto Geologico y Minero de Espana 2018
Subjects:
Online Access:http://hdl.handle.net/10810/44134
https://doi.org/10.21701/bolgeomin.129.3.006
Description
Summary:Here we propose a permutation test-a non-parametric computing-intensive test-to evaluate the statistical significance of the Foster s wavelet spectrum by means of Monte Carlo simulations. A procedure (algorithm) is introduced in order to carry out this aim. The Foster s wavelet is an adequate method to cope directly with unevenly spaced paleoclimatic time series. We have conducted time series simulations to study the performance of the Foster s wavelet spectrum and applied the permutation test to localize periodic signals known a priori by randomly increasing the fraction of missing data (from 25% to 75% of the total amount of data of an evenly spaced time series). We found that the periodic signals are progressively lost as larger amounts of data are removed. This loss becomes noticeable at circa 50% of the missing data and gets more evident when 75% of the data is removed. The signal loss is more pronounced in the high-frequency range, due to an intrinsic bias in the wavelet spectrum. Notwithstanding, the Foster s wavelet spectrum and the permutation test succeeded in locating the periodic signals, even with a moderate amount of missing data. Finally, we applied our procedure to the oxygen isotope ratio (d 18 O) data of the GISP2 deep ice core (Greenland) and we were able to detect the main spectral signature of the unevenly spaced time series of the GISP2 d 18 O record, i.e., the spectral peak around 1,470 years. © 2018, Instituto Geologico y Minero de Espana. All rights reserved. JMPM was funded by a Basque Government post-doctoral fellowship (Ref. No. POS_2015_1_0006). SHF acknowledges financial support from the Ramón y Cajal Research Fellowship of the Ministry of Economy and Competitiveness of Spain (Ref. No. RYC-2012-12167).