Abstract:The use of Gaussian processes (GPs) is a common approach to account for correlated noise in exoplanet time series, particularly for transmission and emission spectroscopy. This analysis has typically been performed for each wavelength channel separately, with the retrieved uncertainties in the transmission spectrum assumed to be independent. However, the presence of noise correlated in wavelength could cause these uncertainties to be correlated, which could significantly affect the results of atmospheric retrievals. We present a method that uses a GP to model noise correlated in both wavelength and time simultaneously for the full spectroscopic dataset. To make this analysis computationally tractable, we introduce a new fast and flexible GP method that can analyse 2D datasets when the input points lie on a (potentially non-uniform) 2D grid - in our case a time by wavelength grid - and the kernel function has a Kronecker product structure. This simultaneously fits all light curves and enables the retrieval of the full covariance matrix of the transmission spectrum. Our new method can avoid the use of a `common-mode' correction, which is known to produce an offset to the transmission spectrum. Through testing on synthetic datasets, we demonstrate that our new approach can reliably recover atmospheric features contaminated by noise correlated in time and wavelength. In contrast, fitting each spectroscopic light curve separately performed poorly when wavelength-correlated noise was present. It frequently underestimated the uncertainty of the scattering slope and overestimated the uncertainty in the strength of sharp absorption peaks in transmission spectra. Two archival VLT/FORS2 transit observations of WASP-31b were used to compare these approaches on real observations. Our method strongly constrained the presence of wavelength-correlated noise in both datasets, and significantly different constraints on atmospheric features such as the scattering slope and strength of sodium and potassium features were recovered.

Robust Detrending of Spatially Correlated Systematics in Kepler Light Curves Using Low-Rank Methods

Robust Detrending of Spatially Correlated Systematics in Kepler Light Curves Using Low-rank Methods

Principal Component Analysis to correct data systematics. Case study: K2 light curves

The SARS algorithm: detrending CoRoT light curves with Sysrem using simultaneous external parameters

Identification and Removal of Noise Modes in Kepler Photometry

Bayesian Methods for Joint Exoplanet Transit Detection and Systematic Noise Characterization

The democratic detrender: Ensemble-Based Removal of the Nuisance Signal in Stellar Time-Series Photometry

Jitter Across 15 Years: Leveraging Precise Photometry from Kepler and TESS to Extract Exoplanets from Radial Velocity Time Series

Exposure-based Algorithm for Removing Systematics out of the CoRoT Light Curves

Utilizing Photometry from Multiple Sources to Mitigate Stellar Variability in Precise Radial Velocities: A Case Study of Kepler-21

Detrending Exoplanetary Transit Light Curves with Long Short-Term Memory Networks

Single Transit Detection In Kepler With Machine Learning And Onboard Spacecraft Diagnostics

Demystifying Kepler Data: A Primer for Systematic Artifact Mitigation

nuance: Efficient detection of planets transiting active stars

Gaussian Processes and Nested Sampling Applied to Kepler's Small Long-period Exoplanet Candidates

EVEREST: Pixel Level Decorrelation of K2 Light curves

A study of two periodogram algorithms for improving the detection of small transiting planets

How do wavelength correlations affect transmission spectra? Application of a new fast and flexible 2D Gaussian process framework to transiting exoplanet spectroscopy