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.
arXiv commentLanguage revised version accepted for publication in Astronomy & Astrophysics. 35 pages, 18 figures (including appendix)
References in corpus (56)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- 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
- Analytic Lightcurves for Planetary Transit Searches
- Detection of an Extrasolar Planet Atmosphere
- batman: BAsic Transit Model cAlculatioN in Python
- Fast and scalable Gaussian process modeling with applications to astronomical time series
- A continuum from clear to cloudy hot-Jupiter exoplanets without primordial water depletion
- Theoretical Transmission Spectra During Extrasolar Giant Planet Transits
- Transmission Spectra as Diagnostics of Extrasolar Giant Planet Atmospheres
- Rayleigh scattering in the transit spectrum of HD 189733b
- A Gaussian process framework for modelling stellar activity signals in radial velocity data
- The Transit Light Source Effect: False Spectral Features and Incorrect Densities for M-dwarf Transiting Planets
- The GAPS Programme with HARPS-N@TNG XIV. Investigating giant planet migration history via improved eccentricity and mass determination for 231 transiting planets
- Parameter Estimation from Time-Series Data with Correlated Errors: A Wavelet-Based Method and its Application to Transit Light Curves
- 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
- A Gaussian process framework for modelling instrumental systematics: application to transmission spectroscopy
- 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
- The Deep Blue Color of HD189733b: Albedo Measurements with Hubble Space Telescope/Space Telescope Imaging Spectrograph at Visible Wavelengths
- ACCESS I: An Optical Transmission Spectrum of GJ 1214b Reveals a Heterogeneous Stellar Photosphere
- ACCESS: A featureless optical transmission spectrum for WASP-19b from Magellan/IMACS
- Probing the haze in the atmosphere of HD 189733b with HST/WFC3 transmission spectroscopy
- VLT FORS2 comparative transmission spectroscopy: Detection of Na in the atmosphere of WASP-39b from the ground
- Avocado: Photometric Classification of Astronomical Transients with Gaussian Process Augmentation
- 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
- Reliable inference of exoplanet light curve parameters using deterministic and stochastic systematics models
- VLT/FORS2 comparative transmission spectroscopy II: confirmation of a cloud-deck and Rayleigh scattering in WASP-31b, but no potassium?
- Improving Cosmological Distance Measurements Using Twin Type Ia Supernovae
- The Very Low Albedo of WASP-12b From Spectral Eclipse Observations with
- WASP-31b: a low-density planet transiting a metal-poor, late-F-type dwarf star
- 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
- High-resolution atmospheric retrievals of WASP-121b transmission spectroscopy with ESPRESSO: Consistent relative abundance constraints across multiple epochs and instruments
- Gemini/GMOS Optical Transmission Spectroscopy of WASP-121b: signs of variability in an ultra-hot Jupiter?
- Revisiting the potassium feature of WASP-31b at high-resolution
- 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
- Simultaneous Optical Transmission Spectroscopy of a Terrestrial, Habitable-Zone Exoplanet with Two Ground-Based Multi-Object Spectrographs
- A Fast, 2D Gaussian Process Method Based on Celerite: Applications to Transiting Exoplanet Discovery and Characterization
- On the safe use of prior densities for Bayesian model selection
- Applications of a Gaussian Process Framework for Modelling of High-Resolution Exoplanet Spectra
- A new method to measure the spectra of transiting exoplanet atmospheres using multi-object spectroscopy
- ACCESS: Confirmation of no potassium in the atmosphere of WASP-31b
- 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