Optimal frequency-domain analysis for spacecraft time series: Introducing the missing-data multitaper power spectrum estimator

Abstract

While the Lomb-Scargle periodogram is foundational to astronomy, it has a significant shortcoming: the variance in the estimated power spectrum does not decrease as more data are acquired. Statisticians have a 60-year history of developing variance-suppressing power spectrum estimators, but most are not used in astronomy because they are formulated for time series with uniform observing cadence and without seasonal or daily gaps. Here we demonstrate how to apply the missing-data multitaper power spectrum estimator to spacecraft data with uniform time intervals between observations but missing data during thruster fires or momentum dumps. The F-test for harmonic components may be applied to multitaper power spectrum estimates to identify statistically significant oscillations that would not rise above a white noise-based false alarm probability. Multitapering improves the dynamic range of the power spectrum estimate and suppresses spectral window artifacts. We show that the multitaper - F-test combination applied to Kepler observations of KIC 6102338 detects differential rotation without requiring iterative sinusoid fitting and subtraction. Significant signals reside at harmonics of both fundamental rotation frequencies and suggest an antisolar rotation profile. Next we use the missing-data multitaper power spectrum estimator to identify the oscillation modes responsible for the complex "scallop shell" shape of the K2 light curve of EPIC 203354381. We argue that multitaper power spectrum estimators should be used for all time series with regular observing cadence.

27 pages, including 15 figures, published in AJ. While this version focuses exclusively on multitaper power spectrum estimation, v1 also demonstrated multitaper magnitude-squared coherence (MSC) and complex demodulation (CDM) for spacecraft data. Separate manuscripts documenting MSC and CDM are in preparation

References in corpus (20)

• The generalised Lomb-Scargle periodogram. A new formalism for the floating-mean and Keplerian periodograms
• Improved Age Estimation for Solar-Type Dwarfs Using Activity-Rotation Diagnostics
• A terrestrial planet candidate in a temperate orbit around Proxima Centauri
• Assessing statistical significance of periodogram peaks
• Radial velocity planets de-aliased. A new, short period for Super-Earth 55 Cnc e
• EVEREST: Pixel Level Decorrelation of K2 Light curves
• BGLS: A Bayesian formalism for the generalised Lomb-Scargle periodogram
• A candidate short-period sub-Earth orbiting Proxima Centauri
• Dynamics of the solar magnetic bright points derived from their horizontal motions
• Orbiting Clouds of Material at the Keplerian Co-Rotation Radius of Rapidly Rotating Low Mass WTTs in Upper Sco
• Efficient modeling of correlated noise II. A flexible noise model with fast and scalable methods
• TESS observations of Cepheid stars: first light results
• More Rapidly Rotating PMS M Dwarfs with Light Curves Suggestive of Orbiting Clouds of Material
• A method to discriminate solar and antisolar differential rotation in high-precision light curves
• Two super-Earths at the edge of the habitable zone of the nearby M dwarf TOI-2095
• Magnitude-squared coherence: A powerful tool for disentangling Doppler planet discoveries from stellar activity
• Characterizing the Motion of Solar Magnetic Bright Points at High Resolution
• Twinkle -- a small satellite spectroscopy mission for the next phase of exoplanet science
• A review of planetary systems around HD 99492, HD 147379 and HD 190007 with HARPS-N
• The Pointing Limits of Transiting Exoplanet Light Curve Characterization with Pixel Level De-correlation

Cited by in corpus (1)

• M${}^2$NuFFT: A Computationally Efficient Suboptimal Power Spectrum Estimator for Fast Exploration of Nonuniformly Sampled Time Series