A Gaussian process framework for modelling instrumental systematics: application to transmission spectroscopy

Abstract

Transmission spectroscopy, which consists of measuring the wavelength-dependent absorption of starlight by a planet's atmosphere during a transit, is a powerful probe of atmospheric composition. However, the expected signal is typically orders of magnitude smaller than instrumental systematics, and the results are crucially dependent on the treatment of the latter. In this paper, we propose a new method to infer transit parameters in the presence of systematic noise using Gaussian processes, a technique widely used in the machine learning community for Bayesian regression and classification problems. Our method makes use of auxiliary information about the state of the instrument, but does so in a non-parametric manner, without imposing a specific dependence of the systematics on the instrumental parameters, and naturally allows for the correlated nature of the noise. We give an example application of the method to archival NICMOS transmission spectroscopy of the hot Jupiter HD 189733, which goes some way towards reconciling the controversy surrounding this dataset in the literature. Finally, we provide an appendix giving a general introduction to Gaussian processes for regression, in order to encourage their application to a wider range of problems.

6 figures, 1 table, accepted for publication in MNRAS

References in corpus (22)

• Analytic Lightcurves for Planetary Transit Searches
• Detection of an Extrasolar Planet Atmosphere
• Detection of Thermal Emission from an Extrasolar Planet
• Infrared radiation from an extrasolar planet
• HST Time-Series Photometry of the Transiting Planet of HD 209458
• The effect of red noise on planetary transit detection
• Detection of atmospheric haze on an extrasolar planet: The 0.55 - 1.05 micron transmission spectrum of HD189733b with the Hubble Space Telescope
• Transmission Spectra as Diagnostics of Extrasolar Giant Planet Atmospheres
• Efficient identification of exoplanetary transit candidates from SuperWASP light curves
• Parameter Estimation from Time-Series Data with Correlated Errors: A Wavelet-Based Method and its Application to Transit Light Curves
• The Transit Light Curve Project. IX. Evidence for a Smaller Radius of the Exoplanet XO-3b
• Stellar Limb-Darkening Coefficients for CoRot and Kepler
• Hubble Space Telescope times-series photometry of the planetary transit of HD189733: no moon, no rings, starspots
• The Transit Light Curve (TLC) Project. I. Four Consecutive Transits of the Exoplanet XO-1b
• A new look at NICMOS transmission spectroscopy of HD189733, GJ-436 and XO-1: no conclusive evidence for molecular features
• A transit timing analysis of nine RISE light curves of the exoplanet system TrES-3
• Re-assessing the radial-velocity evidence for planets around CoRoT-7
• Novel Methods for Predicting Photometric Redshifts from Broad Band Photometry using Virtual Sensors
• Transit infrared spectroscopy of the hot neptune around GJ 436 with the Hubble Space Telescope
• A transit timing analysis of seven RISE light curves of the exoplanet system HAT-P-3
• New Approaches To Photometric Redshift Prediction Via Gaussian Process Regression In The Sloan Digital Sky Survey
• Towards Real-time Classification of Astronomical Transients

Cited by in corpus (208)

• Fast and scalable Gaussian process modeling with applications to astronomical time series
• Infrared Transmission Spectroscopy of the Exoplanets HD209458b and XO-1b Using the Wide Field Camera-3 on the Hubble Space Telescope
• Planets and Stellar Activity: Hide and Seek in the CoRoT-7 system
• The prevalence of dust on the exoplanet HD 189733b from Hubble and Spitzer observations
• A Gaussian process framework for modelling stellar activity signals in radial velocity data
• Juliet: a versatile modelling tool for transiting and non-transiting exoplanetary systems
• 3.6 and 4.5 Micron Phase Curves and Evidence for Non-Equilibrium Chemistry in the Atmosphere of Extrasolar Planet HD 189733b
• Tau-REx I: A next generation retrieval code for exoplanetary atmospheres
• Early Release Science of the Exoplanet WASP-39b with JWST NIRSpec G395H
• Exoplanet population inference and the abundance of Earth analogs from noisy, incomplete catalogs
• Wotan: Comprehensive time-series de-trending in Python
• A rocky planet transiting a nearby low-mass star
• Inferring probabilistic stellar rotation periods using Gaussian processes
• Detection of H2O and evidence for TiO/VO in an ultra hot exoplanet atmosphere
• Detection of Fe I in the atmosphere of the ultra-hot Jupiter WASP-121b, and a new likelihood-based approach for Doppler-resolved spectroscopy
• A consistent retrieval analysis of 10 Hot Jupiters observed in transmission
• K2SC: Flexible systematics correction and detrending of K2 light curves using Gaussian Process regression
• Detection of titanium oxide in the atmosphere of a hot Jupiter
• An HST Optical to Near-IR Transmission Spectrum of the Hot Jupiter WASP-19b: Detection of Atmospheric Water and Likely Absence of TiO
• The Deep Blue Color of HD189733b: Albedo Measurements with Hubble Space Telescope/Space Telescope Imaging Spectrograph at Visible Wavelengths
• Transit Detection in the MEarth Survey of Nearby M Dwarfs: Bridging the Clean-First, Search-Later Divide
• No Thermal Inversion and a Solar Water Abundance for the Hot Jupiter HD209458b from HST WFC3 Emission Spectroscopy
• The Complete transmission spectrum of WASP-39b with a precise water constraint
• An optical transmission spectrum for the ultra-hot Jupiter WASP-121b measured with the Hubble Space Telescope
• An absolute sodium abundance for a cloud-free 'hot Saturn' exoplanet
• Determining the Mass of Kepler-78b With Nonparametric Gaussian Process Estimation
• Characterizing Transiting Planet Atmospheres through 2025
• H2O abundances in the atmospheres of three hot Jupiters
• A Gemini ground-based transmission spectrum of WASP-29b: a featureless spectrum from 515 to 720nm
• Probing the haze in the atmosphere of HD 189733b with HST/WFC3 transmission spectroscopy
• On Correlated-noise Analyses Applied To Exoplanet Light Curves
• The long egress of GJ~436b's giant exosphere
• Searching for Exoplanets Using Artificial Intelligence
• A uniform analysis of HD209458b Spitzer/IRAC lightcurves with Gaussian process models
• Precise time-series photometry for the Kepler-2.0 mission
• Using Machine Learning for Discovery in Synoptic Survey Imaging
• The spectrum of hot methane in astronomical objects using a comprehensive computed line list
• Large eccentricity, low mutual inclination: the three-dimensional architecture of a hierarchical system of giant planets
• A Ground-based Optical Transmission Spectrum of WASP-6b
• VLT FORS2 comparative transmission spectroscopy: Detection of Na in the atmosphere of WASP-39b from the ground
• The optical transmission spectrum of the hot Jupiter HAT-P-32b: clouds explain the absence of broad spectral features?
• Transit Timing Observations from Kepler: II. Confirmation of Two Multiplanet Systems via a Non-parametric Correlation Analysis
• The GTC exoplanet transit spectroscopy survey IX. Detection of Haze, Na, K, and Li in the super-Neptune WASP-127b
• Features in the Broadband Eclipse Spectra of Exoplanets: Signal or Noise?
• The Hunt for Exomoons with Kepler (HEK): V. A Survey of 41 Planetary Candidates for Exomoons
• Repeatability and Accuracy of Exoplanet Eclipse Depths Measured with Post-Cryogenic Spitzer
• Radial Velocity Observations and Light Curve Noise Modeling Confirm That Kepler-91b is a Giant Planet Orbiting a Giant Star
• Clouds on the hot Jupiter HD189733b: constraints from the reflection spectrum
• An emission spectrum for WASP-121b measured across the 0.8-1.1 micron wavelength range using the Hubble Space Telescope
• Let the Great World Spin: Revealing the Stormy, Turbulent Nature of Young Giant Exoplanet Analogs with the Spitzer Space Telescope
• Rayleigh scattering in the transmission spectrum of HAT-P-18b
• Reliable inference of exoplanet light curve parameters using deterministic and stochastic systematics models
• HAT-P-57b: A Short-Period Giant Planet Transiting A Bright Rapidly Rotating A8V Star Confirmed Via Doppler Tomography
• Measurement of the Nodal Precession of WASP-33 b via Doppler Tomography
• VLT/FORS2 comparative transmission spectroscopy II: confirmation of a cloud-deck and Rayleigh scattering in WASP-31b, but no potassium?
• Warm Ice Giant GJ 3470b. I. A Flat Transmission Spectrum Indicates a Hazy, Low-methane, and/or Metal-rich Atmosphere
• ACCESS & LRG-BEASTS: a precise new optical transmission spectrum of the ultrahot Jupiter WASP-103b
• Atmospheric Characterization of 5 Hot Jupiters with Wide Field Camera 3 on the Hubble Space Telescope
• Unveiling cloudy exoplanets: the influence of cloud model choices on retrieval solutions
• Near-infrared Thermal Emission Detections of a number of hot Jupiters and the Systematics of Ground-based Near-infrared Photometry
• Evidence for Atmospheric Cold-trap Processes in the Noninverted Emission Spectrum of Kepler-13Ab Using HST/WFC3
• A Comprehensive Reanalysis of $\text{Spitzer}$'s 4.5 $\mu$m Phase Curves, and the Phase Variations of the Ultra-hot Jupiters MASCARA-1b and KELT-16b
• Blind extraction of an exoplanetary spectrum through Independent Component Analysis
• The Very Low Albedo of WASP-12b From Spectral Eclipse Observations with $\text{Hubble}$
• Precision radial velocity measurements by the forward-modeling technique in the near-infrared
• MuSCAT2 multicolour validation of TESS candidates: an ultra-short-period substellar object around an M dwarf
• Physical Parameters of Late-type Contact Binaries in the Northern Catalina Sky Survey
• Ks band secondary eclipses of WASP-19b and WASP-43b with the Anglo-Australian Telescope
• LRG-BEASTS: transmission spectroscopy and retrieval analysis of the highly-inflated Saturn-mass planet WASP-39b
• Potassium detection in the clear atmosphere of a hot-Jupiter: FORS2 transmission spectroscopy of WASP-17b
• Kepler423b: a half-Jupiter mass planet transiting a very old solar-like star
• WASP-127b transmission spectrum]{Abundance measurements of H${}_2$O and carbon-bearing species in the atmosphere of WASP-127b confirm its super-solar metallicity
• Transmission spectroscopy of exoplanet XO-2b observed with HST NICMOS
• No Conclusive Evidence for Transits of Proxima b in MOST photometry
• CoRoT 223992193: A new, low-mass, pre-main sequence eclipsing binary with evidence of a circumbinary disk
• Evidence for disequilibrium chemistry from vertical mixing in hot Jupiter atmospheres. A comprehensive survey of transiting close-in gas giant exoplanets with warm-Spitzer/IRAC
• Improving transit characterisation with Gaussian process modelling of stellar variability
• Spi-OPS: Spitzer and CHEOPS confirm the near-polar orbit of MASCARA-1 b and reveal a hint of dayside reflection
• LRG-BEASTS III: Ground-based transmission spectrum of the gas giant orbiting the cool dwarf WASP-80
• Near-Infrared Emission Spectrum of WASP-103b using Hubble Space Telescope/Wide Field Camera 3
• Detection of Na, K and H${}_2$O in the hazy atmosphere of WASP-6b
• Very metal-poor stars observed by the RAVE survey
• A 0.8-2.4 Micron Transmission Spectrum of the Hot Jupiter CoRoT-1b
• SPIDERMAN: an open-source code to model phase curves and secondary eclipses
• Spin orbit alignment for KELT-7b and HAT-P-56b via Doppler tomography with TRES
• Gemini/GMOS Optical Transmission Spectroscopy of WASP-121b: signs of variability in an ultra-hot Jupiter?
• TOI-503: The first known brown dwarf-Am star binary from the TESS mission
• The Broadband and Spectrally-Resolved H-band Eclipse of KELT-1b and the Role of Surface Gravity in Stratospheric Inversions in Hot Jupiters
• A statistical approach to identify superluminous supernovae and probe their diversity
• On signals faint and sparse: the ACICA algorithm for blind de-trending of Exoplanetary transits with low signal-to-noise
• On the Detection of Molecules in the Atmosphere of HD189733b using HST NICMOS Transmission Spectroscopy
• Revisiting the exomoon candidate signal around Kepler-1625b
• The ground-based H, K, and L-band absolute emission spectra of HD 209458b
• HST PanCET Program: A Complete Near-UV to Infrared Transmission Spectrum for the Hot Jupiter WASP-79b
• The GTC exoplanet transit spectroscopy survey II: An overly-large Rayleigh-like feature for exoplanet TrES-3b
• TOI-674b: an oasis in the desert of exo-Neptunes transiting a nearby M dwarf
• The K2-ESPRINT Project III: A Close-in Super-Earth around a Metal-rich Mid-M Dwarf
• Revisiting the potassium feature of WASP-31b at high-resolution
• Exploring the GRB population: Robust afterglow modelling
• The GTC exoplanet transit spectroscopy survey XI. Possible detection of Rayleigh scattering in the atmosphere of the Saturn-mass planet WASP-69b
• A Featureless Infrared Transmission Spectrum for the Super-Puff Planet Kepler-79d
• An Unusual Transmission Spectrum for the Sub-Saturn KELT-11b Suggestive of a Sub-Solar Water Abundance
• AutoRegressive Planet Search: Application to the Kepler Mission
• The geometric albedo of the hot Jupiter HD 189733b measured with CHEOPS
• Spin-Orbit Misalignments of Three Jovian Planets via Doppler Tomography
• Regaining the FORS: optical ground-based transmission spectroscopy of the exoplanet WASP-19b with VLT+FORS2
• Multiple Clues for Dayside Aerosols and Temperature Gradients in WASP-69 b from a Panchromatic JWST Emission Spectrum
• Ground-Based Transmission Spectroscopy with FORS2: A featureless optical transmission spectrum and detection of H${}_2$O for the ultra-hot Jupiter WASP-103b
• Mitigating flicker noise in high-precision photometry. I -- Characterization of the noise structure, impact on the inferred transit parameters, and predictions for CHEOPS observations
• LRG-BEASTS: sodium absorption and Rayleigh scattering in the atmosphere of WASP-94A b using NTT/EFOSC2
• Detection of Na in WASP-21b's lower and upper atmosphere
• ACCESS: A Visual to Near-infrared Spectrum of the Hot Jupiter WASP-43b with Evidence of ${\mathrm{H}}_2\mathrm{O}$, but no evidence of Na or K
• Simultaneous Optical Transmission Spectroscopy of a Terrestrial, Habitable-Zone Exoplanet with Two Ground-Based Multi-Object Spectrographs
• The HARPS search for southern extra-solar planets. XXXIX. HD175607 b, the most metal-poor G dwarf with an orbiting sub-Neptune
• What can the 2008/10 broadband flare of PKS 1502+106 tell us? Nuclear opacity, magnetic fields, and the location of gamma rays
• Solar-to-supersolar sodium and oxygen absolute abundances for a "hot Saturn" orbiting a metal-rich star
• Benchmark Tests for Markov Chain Monte Carlo Fitting of Exoplanet Eclipse Observations
• Probing the atmosphere of WASP-69 b with low- and high-resolution transmission spectroscopy
• Exoplanet Transmission Spectroscopy using KMOS
• Spitzer transit follow-up of planet candidates from the K2 mission
• Modelling the Multiwavelength Variability of Mrk 335 using Gaussian Processes
• The phase curve and the geometric albedo of WASP-43b measured with CHEOPS, TESS and HST WFC3/UVIS
• Evidence of a Clear Atmosphere for WASP-62b: the Only Known Transiting Gas Giant in the JWST Continuous Viewing Zone
• Predicting binding energies of astrochemically relevant molecules via machine learning
• Probing the atmosphere of a sub-Jovian planet orbiting a cool dwarf
• The TESS Phase Curve of KELT-1b Suggests a High Dayside Albedo
• A Gaussian Process Regression Reveals No Evidence for Planets Orbiting Kapteyn's Star
• Using Gaussian processes to model light curves in the presence of flickering: the eclipsing cataclysmic variable ASASSN-14ag
• Searching for quasi-periodic oscillations in astrophysical transients using Gaussian processes
• LRG-BEASTS: Ground-based Detection of Sodium and a Steep Optical Slope in the Atmosphere of the Highly Inflated Hot-Saturn WASP-21b
• Demonstrating High-precision, Multiband Transit Photometry with MuSCAT: A Case for HAT-P-14b
• An enhanced slope in the transmission spectrum of the hot Jupiter WASP-104b
• Orbital alignment of the eccentric warm Jupiter TOI-677 b
• ACCESS: Confirmation of a Clear Atmosphere for WASP-96b and a Comparison of Light Curve Detrending Techniques
• The Dark World: A Tale of WASP-43b in Reflected Light with HST WFC3/UVIS
• Global Attenuation in Spiral Galaxies in Optical and Infrared Bands
• Trend Filtering -- I. A Modern Statistical Tool for Time-Domain Astronomy and Astronomical Spectroscopy
• Finding Black Holes with Black Boxes -- Using Machine Learning to Identify Globular Clusters with Black Hole Subsystems
• The GTC exoplanet transit spectroscopy survey VIII. Flat transmission spectrum for the warm gas giant WASP-80
• Hot Rocks Survey III: A deep eclipse for LHS 1140c and a new Gaussian process method to account for correlated noise in individual pixels
• CHEOPS finds KELT-1b darker than expected in visible light: Discrepancy between the CHEOPS and TESS eclipse depths
• A new method to measure the spectra of transiting exoplanet atmospheres using multi-object spectroscopy
• The GTC exoplanet transit spectroscopy survey X. Stellar spots versus Rayleigh scattering: the case of HAT-P-11b
• Evidence for TiO in the atmosphere of the hot Jupiter HAT-P-65 b
• Applications of a Gaussian Process Framework for Modelling of High-Resolution Exoplanet Spectra
• ACCESS: Confirmation of no potassium in the atmosphere of WASP-31b
• Detrending Exoplanetary Transit Light Curves with Long Short-Term Memory Networks
• HD 20329b: An ultra-short-period planet around a solar-type star found by TESS
• OGLE-2017-BLG-1186: first application of asteroseismology and Gaussian processes to microlensing
• Transmission spectroscopy with VLT FORS2: a featureless spectrum for the low-density transiting exoplanet WASP-88b
• Determining Empirical Stellar Masses and Radii Using Transits and Gaia Parallaxes as Illustrated by Spitzer Observations of KELT-11b
• Featureless transmission spectra of 12 giant exoplanets observed by GTC/OSIRIS
• A blind method to detrend instrumental systematics in exoplanetary light-curves
• Estimating dayside effective temperatures of hot Jupiters and associated uncertainties through Gaussian process regression
• Searching for the Transit of the Earth--mass exoplanet Proxima~Centauri~b in Antarctica: Preliminary Result
• Measuring the variability of directly imaged exoplanets using vector Apodizing Phase Plates combined with ground-based differential spectrophotometry
• Revisiting KELT-19Ab, WASP-156b and WASP-121b in the TESS Era
• Evidence for stellar contamination in the transmission spectra of HAT-P-12b
• Constraints on TESS albedos for five hot Jupiters
• A new approach to modelling gamma-ray burst afterglows: Using Gaussian processes to account for the systematics
• Allesfitter: Flexible Star and Exoplanet Inference From Photometry and Radial Velocity
• Validation and atmospheric exploration of the sub-Neptune TOI-2136b around a nearby M3 dwarf
• ACCESS: An optical transmission spectrum of the high-gravity, hot Jupiter HAT-P-23b
• TOI-1442 b and TOI-2445 b: two potentially rocky ultra-short period planets around M dwarfs
• Detection of faculae in the transit and transmission spectrum of WASP-69b
• Detection of Na and K in the atmosphere of the hot Jupiter HAT-P-1b with P200/DBSP
• A precise blue-optical transmission spectrum from the ground: Evidence for haze in the atmosphere of WASP-74b
• Tentative detection of titanium oxide in the atmosphere of WASP-69 b with a 4m ground-based telescope
• Absolute dimensions and apsidal motion of the eclipsing binaries V889 Aql and V402 Lac
• Exploring the directly imaged HD 1160 system through spectroscopic characterization and high-cadence variability monitoring
• WASP-35 and HAT-P-30/WASP-51: re-analysis using TESS and ground-based transit photometry
• An old warm Jupiter orbiting the metal-poor G-dwarf TOI-5542
• Ground-Based Transmission Spectroscopy with VLT FORS2: Evidence for faculae and clouds in the optical spectrum of the warm Saturn WASP-110b
• Knot a Bad Idea: Testing BLISS Mapping for Spitzer Space Telescope Photometry
• Ground-based Transit Observation of the Habitable-zone Super-Earth K2-3d
• Enabling Robust Exoplanet Atmospheric Retrievals with Gaussian Processes
• An Empirical Bayesian Approach to Limb-darkening in Modeling WASP-121b Transit Light Curves
• Benchmark tests of transmission spectroscopy using transiting white dwarfs
• Dynamical modelling of dwarf-spheroidal galaxies using Gaussian-process emulation
• The study on transiting systems HAT-P-13, HAT-P-16 and WASP-32 through combining ground-based and TESS photometry
• Discovery of a cold giant planet and mass measurement of a hot super-Earth in the multi-planetary system WASP-132
• How do wavelength correlations affect transmission spectra? Application of a new fast and flexible 2D Gaussian process framework to transiting exoplanet spectroscopy
• SSSpaNG! Stellar Spectra as Sparse, data-driven, Non-Gaussian processes
• Gaussian Process Modelling for Improved Resolution in Faraday Depth Reconstruction
• A statistical model of stellar variability. I. FENRIR: a physics-based model of stellar activity, and its fast Gaussian process approximation
• Simultaneous multicolour transit photometry of hot Jupiters HAT-P-19b, HAT-P-51b, HAT-P-55b, and HAT-P-65b
• Gaussian Processes and Nested Sampling Applied to Kepler's Small Long-period Exoplanet Candidates
• Bayesian Methods for Exoplanet Science
• Kepler-102: Masses and Compositions for a Super-Earth and Sub-Neptune Orbiting an Active Star
• Data-Driven Sensitivity Inference for Thomson Scattering Electron Density Measurement Systems
• A Transmission Spectrum Of HD 189733b From Multiple Broadband Filter Observations
• Spectrophotometric redshifts for $\mathrm{z}\sim 1$ galaxies and predictions for number densities with WFIRST and Euclid
• AutoRegressive Planet Search: Feasibility Study for Irregular Time Series
• Probabilistic orbits and dynamical masses of emission-line binaries
• Kernel-, mean- and noise-marginalised Gaussian processes for exoplanet transits and $H_0$ inference
• Detection of Na in the atmosphere of the hot Jupiter HAT-P-55b
• Temperate super-Earths/mini-Neptunes around M/K dwarfs Consist of 2 Populations Distinguished by Their Atmospheres
• An improved correction of radial-velocity systematic for the SOPHIE spectrograph
• A Two-limb Explanation for the Optical-to-infrared Transmission Spectrum of the Hot Jupiter HAT-P-32Ab
• Chasing the storm: Investigating the application of high-contrast imaging techniques in producing precise exoplanet light curves
• Evidence for stellar contamination and water absorption in NGTS-5b's transmission spectra with GTC/OSIRIS
• gallifrey: JAX-based Gaussian process structure learning for astronomical time series
• A new method to correct for host star variability in multi-epoch observations of exoplanet transmission spectra
• Gaussian processes for blazar variability studies
• Exposure-averaged Gaussian Processes for Combining Overlapping Datasets
• Using chromatic covariance to correct for scintillation noise in ground-based spectrophotometry
• Photometry on Structured Backgrounds: Local Pixelwise Infilling by Regression
• Panchromatic View of the Frigid Jovian Exoplanet COCONUTS-2 b