Estimation of the significance of the Foster’s wavelet spectrum by means of a permutation test and its application for paleoclimate records

• J.M. Polanco-Martínez ^[1] ; S.H. Faria ^[1]
  1. [1] BC3 Basque Centre for Climate Change

     BC3 Basque Centre for Climate Change

     Leioa, España

     
• Localización: Boletín geológico y minero, ISSN 0366-0176, Vol. 129, Nº 3, 2018 (Ejemplar dedicado a: Análisis de series temporales en Ciencias de la Tierra ), págs. 549-564
• Idioma: inglés
• Títulos paralelos:
  • Estimación de la significación estadística del espectro wavelet de Foster mediante el test de permutación y aplicaciones en paleoclima
• Enlaces
  • Texto completo (pdf)
• Resumen
  • español

    En este trabajo se propone un test (no paramétrico pero de cálculo intensivo) de permutación para evaluar la significación estadística del espectro wavelet de Foster mediante simulaciones de Monte Carlo. Un procedimiento (algoritmo) se introduce para llevar a cabo esta tarea. Los wavelet de Foster son un método adecuado para analizar directamente series paleoclimáticas irregulares. Se simularon series temporales con conocida información espectral (eventos periódicos) con el objeto de evaluar el desempeño en la detección de eventos periódicos conocidos a priori, incrementando el porcentaje de datos faltantes de manera aleatoria (desde un 25% a un 75% de la cantidad total de datos de una serie de control equiespaciada). Se encontró que ante la presencia de una gran cantidad de datos faltantes (lo cuál es evidente desde un 50% y más pronunciado con un 75%), las señales periódicas no se detectan en el espectro wavelet. Esto es más evidente en las señales de corta frecuencia debido a la existencia de un sesgo en el espectro wavelet de Foster. Sin embargo, a pesar de lo anterior, el espectro wavelet de Foster en combinación con el test de permutación detectan adecuadamente señales periódicas en series sintéticas (simuladas) y reales con una moderada cantidad de datos faltantes. Finalmente, se aplicó este método a valores medios de δ18O de los testigos de hielo del GISP2 (Groenlandia) y se encontró la principal característica espectral de estos registros, esto es, la señal espectral ca. 1470 años.

  • English

    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 (δ18O) 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 δ18O record, i.e., the spectral peak around 1,470 years.