How do wavelength correlations affect transmission spectra? Application of a new fast and flexible 2D Gaussian process framework to transiting exoplanet spectroscopy
arXiv:2402.15204 · doi:10.1051/0004-6361/202347613
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 while avoiding the use of a 'common-mode' correction. 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 covariance matrix of the transmission spectrum. By 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 re-analysed, with our method strongly constraining the presence of wavelength-correlated noise in both datasets and recovering significantly different constraints on atmospheric features such as the scattering slope and strength of sodium and potassium features.
Language revised version accepted for publication in Astronomy & Astrophysics. 35 pages, 18 figures (including appendix)
References in corpus (32)
- Array Programming with NumPy
- The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
- Rayleigh scattering in the transit spectrum of HD 189733b
- Early Release Science of the exoplanet WASP-39b with JWST NIRSpec PRISM
- Water Vapor in the Spectrum of the Extrasolar Planet HD 189733b: 1. the Transit
- Early Release Science of the Exoplanet WASP-39b with JWST NIRSpec G395H
- Thermal emission from the Earth-sized exoplanet TRAPPIST-1 b using JWST
- HST hot-Jupiter transmission spectral survey: detection of potassium in WASP-31b along with a cloud deck and Rayleigh scattering
- exoplanet: Gradient-based probabilistic inference for exoplanet data & other astronomical time series
- No thick carbon dioxide atmosphere on the rocky exoplanet TRAPPIST-1 c
- Early Release Science of the exoplanet WASP-39b with JWST NIRCam
- Detection of titanium oxide in the atmosphere of a hot Jupiter
- ACCESS I: An Optical Transmission Spectrum of GJ 1214b Reveals a Heterogeneous Stellar Photosphere
- VLT FORS2 comparative transmission spectroscopy: Detection of Na in the atmosphere of WASP-39b from the ground
- Transit Signatures of Inhomogeneous Clouds on Hot Jupiters: Insights From Microphysical Cloud Modeling
- pyaneti II: A multidimensional Gaussian process approach to analysing spectroscopic time-series
- Empirical limb-darkening coefficients & transit parameters of known exoplanets from TESS
- VLT/FORS2 comparative transmission spectroscopy II: confirmation of a cloud-deck and Rayleigh scattering in WASP-31b, but no potassium?
- The Very Low Albedo of WASP-12b From Spectral Eclipse Observations with
- ExoClock project II: A large-scale integrated study with 180 updated exoplanet ephemerides
- Potassium detection in the clear atmosphere of a hot-Jupiter: FORS2 transmission spectroscopy of WASP-17b
- Efficient modeling of correlated noise. III. Scalable methods for jointly modeling several observables' time series with Gaussian processes
- A Spectral Survey of WASP-19b with ESPRESSO
- Gemini/GMOS Optical Transmission Spectroscopy of WASP-121b: signs of variability in an ultra-hot Jupiter?
- High-resolution atmospheric retrievals of WASP-121b transmission spectroscopy with ESPRESSO: Consistent relative abundance constraints across multiple epochs and instruments
- Exoplanet Spectroscopy with JWST NIRISS: Diagnostics and Case Studies
- LRG-BEASTS: sodium absorption and Rayleigh scattering in the atmosphere of WASP-94A b using NTT/EFOSC2
- On the safe use of prior densities for Bayesian model selection
- A new method to measure the spectra of transiting exoplanet atmospheres using multi-object spectroscopy
- Applications of a Gaussian Process Framework for Modelling of High-Resolution Exoplanet Spectra
- Impacts of high-contrast image processing on atmospheric retrievals
- A precise blue-optical transmission spectrum from the ground: Evidence for haze in the atmosphere of WASP-74b
Cited by in corpus (4)
- Hot Rocks Survey III: A deep eclipse for LHS 1140c and a new Gaussian process method to account for correlated noise in individual pixels
- Detection of faculae in the transit and transmission spectrum of WASP-69b
- Enabling Robust Exoplanet Atmospheric Retrievals with Gaussian Processes
- A statistical model of stellar variability. I. FENRIR: a physics-based model of stellar activity, and its fast Gaussian process approximation