Peering into the black box: forward-modeling the uncertainty budget of high-resolution spectroscopy of exoplanet atmospheres

TL;DR

This paper develops a forward-modeling simulation tool to analyze and understand the uncertainty sources and biases in high-resolution spectroscopy of exoplanet atmospheres, enhancing confidence in the method's results.

Contribution

It introduces 'scope', a comprehensive HRCCS observation simulator that models spectral contributions and tests the impact of various systematic effects on detection and measurement accuracy.

Findings

Principal component analysis does not bias signals with few components.

Mild telluric variations only slightly reduce detection significance.

Strong telluric changes can bias velocities and gas abundance estimates.

Abstract

Ground-based high-resolution cross-correlation spectroscopy (HRCCS; R >~ 15,000) is a powerful complement to space-based studies of exoplanet atmospheres. By resolving individual spectral lines, HRCCS can precisely measure chemical abundance ratios, directly constrain atmospheric dynamics, and robustly probe multidimensional physics. But the subtleties of HRCCS datasets -- e.g., the lack of exoplanetary spectra visible by eye and the statistically complex process of telluric removal -- can make interpreting them difficult. In this work, we seek to clarify the uncertainty budget of HRCCS with a forward-modeling approach. We present a HRCCS observation simulator, scope (https://github.com/arjunsavel/scope), that incorporates spectral contributions from the exoplanet, star, tellurics, and instrument. This tool allows us to control the underlying dataset, enabling controlled experimentation with complex HRCCS methods. Simulating a fiducial hot Jupiter dataset (WASP-77Ab emission with IGRINS), we first confirm via multiple tests that the commonly used principal components analysis does not bias the planetary signal when few components are used. Furthermore, we demonstrate that mildly varying tellurics and moderate wavelength solution errors induce only mild decreases in HRCCS detection significance. However, limiting-case, strongly varying tellurics can bias the retrieved velocities and gas abundances. Additionally, in the low-SNR limit, constraints on gas abundances become highly non-Gaussian. Our investigation of the uncertainties and potential biases inherent in HRCCS data analysis enables greater confidence in scientific results from this maturing method.