We present a new spectral analysis method for the identification of periodic signals in geophysical time series. We evaluate the power spectral density with the adaptive multitaper method, a non-parametric spectral analysis technique suitable for time series characterized by colored power spectral density. Our method provides a maximum likelihood estimation of the power spectral density background according to four different models. It includes the option for the models to be fitted on four smoothed versions of the power spectral density when there is a need to reduce the influence of power enhancements due to periodic signals. We use a statistical criterion to select the best background representation among the different smoothing+model pairs. Then, we define the confidence thresholds to identify the power spectral density enhancements related to the occurrence of periodic fluctuations (γ test). We combine the results with those obtained with the multitaper harmonic F test, an additional complex-valued regression analysis from which it is possible to estimate the amplitude and phase of the signals. We demonstrate the algorithm on Monte Carlo simulations of synthetic time series and a case study of magnetospheric field fluctuations directly driven by periodic density structures in the solar wind. The method is robust and flexible. Our procedure is freely available as a stand-alone IDL code at https://zenodo.org/record/3703168. The modular structure of our methodology allows the introduction of new smoothing methods and models to cover additional types of time series. The flexibility and extensibility of the technique makes it broadly suitable to any discipline.

Identifying the nature and source of Ultra Low Frequencies (ULF) waves (f ≤ 4 mHz) at discrete frequencies in the Earth’s magnetosphere is a complex task. The challenge comes from the simultaneous occurrence of externally and internally generated waves, and the ability to robustly identify such perturbations. Using a recently developed robust spectral analysis procedure, we study an interval that exhibited in magnetic field measurements at geosynchronous orbit and in ground magnetic observatories both internally supported and externally generated ULF waves. The event occurred on November 9, 2002 during the interaction of the magnetosphere with two interplanetary shocks that were followed by a train of 90 min solar wind periodic density structures. Using the Wang-Sheeley-Arge model, we mapped the source of this solar wind stream to an active region and a mid-latitude coronal hole just prior to crossing the Heliospheric current sheet. In both the solar wind density and magnetospheric field fluctuations, we separated broad power increases from enhancements at specific frequencies. For the waves at discrete frequencies, we used the combination of satellite and ground magnetometer observations to identify differences in frequency, polarization, and observed magnetospheric locations. The magnetospheric response was characterized by: (i) forced breathing by periodic solar wind dynamic pressure variations below ≈ 1 mHz; (ii) a combination of directly driven oscillations and wave modes triggered by additional mechanisms (e.g., shock and interplanetary magnetic field discontinuity impact, and substorm activity) between ≈ 1 and ≈ 4 mHz; and (iii) largely triggered modes above ≈ 4 mHz.