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.

arXiv comment27 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 (36)

• Astropy: A Community Python Package for Astronomy
• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• 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
• Understanding the Lomb-Scargle Periodogram
• Kepler Presearch Data Conditioning II - A Bayesian Approach to Systematic Error Correction
• Kepler Presearch Data Conditioning I - Architecture and Algorithms for Error Correction in Kepler Light Curves
• Overview of the Kepler Science Processing Pipeline
• Large-scale magnetic topologies of mid-M dwarfs
• Assessing statistical significance of periodogram peaks
• Rotation and differential rotation of active Kepler stars
• Radial velocity planets de-aliased. A new, short period for Super-Earth 55 Cnc e
• Preparation of Kepler lightcurves for asteroseismic analyses
• Asteroseismology of red giants from the first four months of Kepler data: Fundamental parameters
• EVEREST: Pixel Level Decorrelation of K2 Light curves
• Detectability of Extrasolar Planets in Radial Velocity Surveys
• SOAP. A tool for the fast computation of photometry and radial velocity induced by stellar spots
• BGLS: A Bayesian formalism for the generalised Lomb-Scargle periodogram
• A candidate short-period sub-Earth orbiting Proxima Centauri
• The 55 Cnc system reassessed
• The Atacama Cosmology Telescope: A Measurement of the 600< ell <8000 Cosmic Microwave Background Power Spectrum at 148 GHz
• Magnetic activity and differential rotation in the very young star KIC 8429280
• 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
• Magnetic activity and differential rotation in the young Sun-like stars KIC 7985370 and KIC 7765135
• 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