LSTperiod software: spectral analysis of multiple irregularly sampled time series

Author(s)

Caminha-Maciel, George ^{[1, 2]} ; Ernesto, Marcia ^[3]

Affiliation

^[1] Univ Hawaii Manoa, SOEST, HIGP, Petrofabr & Paleomagnet Lab, 1680 East West Rd, Honolulu, HI 96822 - USA

^[2] Univ Fed Santa Catarina, Dept Geociencias, COMPUTAGEO Lab Geofis Computac, Campus Univ Trindade, BR-88040900 Florianopolis, SC - Brazil

^[3] Univ Sao Paulo, Inst Astron Geofis & Ciencias Atmosfer, Dept Geofis, Rua Matao, 1226, Cidade Univ, BR-05508090 Sao Paulo - Brazil

• Document type: Journal article
• Source: ANNALS OF GEOPHYSICS; v. 62, 2019.
• Web of Science Citations: 0

Abstract

Irregularly sampled time series are common in several different areas, such as astronomy, meteorology, biology, oceanography, cyclostratigraphy, and others. The periodogram is a primary tool to extract meaningful information from irregularly spaced and noisy time series. It is an element of decision theory, meaning the periodogram usually transforms the data, and its ordinates are subsequently submitted to a statistical test compared to a population originating from a known stochastic model (white Gaussian noise). If some ordinate f(0) (usually a local maximum, a peak) fails in this test, we declare that it is a `periodicity' at a frequency f(0). Besides its full usage, this method until now suffer from numerous theoretical difficulties in adapting to real case situations and shows lack of usefulness for very poorly sampled and high noise cases. All of it implies low usefulness for applying in most sedimentary sequences at our disposal nowadays. The LSTperiod is an application, written in Matlab, conceived to perform spectral analysis of multiple irregularly sampled time series. It combines information from Lomb-Scargle periodogram estimates over different time series sampling the same phenomenon, enabling the recovering of signals from very poorly sampled and noisy time series. The software comprises a set of four Graphical User Interfaces (GUIs) that allow the user to: (1) Have broad choices of the frequency-domain range and density for spectral estimation; (2) Select possible spectral features (i.e., pick ``T{''}) for testing as a model {[} A{*}sin(2 pi/T t - theta)] through the visualization of several goodness-of-fit statistics; (3) Visualize the fitting residuals in the time domain, for each time series, for the chosen sinusoidal model. These tools help the user to identify and analyze any suspected feature in the estimated spectra through its related linear system responses. All estimated parameter can be saved on worksheets and the visualizations in several different figure formats. We illustrate the use of the software with a set of Ocean Drilling Program (ODP) data series that show long-period Milankovitch-related spectral features and demonstrate its performance using synthetic time series. (AU)