Draft version March 20, 2018 Preprint typeset using LATEX style emulateapj v. 12/16/11 CHARACTERIZING EARTH ANALOGS IN REFLECTED LIGHT: ATMOSPHERIC RETRIEVAL STUDIES FOR FUTURE SPACE TELESCOPES Y. Katherina Feng1,2,3, Tyler D. Robinson1,4,5, Jonathan J. Fortney1,3, Roxana E. Lupu6,7, Mark S. Marley6, Nikole K. Lewis8,9, Bruce Macintosh10, Michael R. Line11 Draft version March 20, 2018 ABSTRACT Space-based high contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, an instrument noise model, and an ensemble Markov chain Monte Carlo sampler, we explore retrievals of atmospheric and planetary properties for Earth twins as a function of signal-to-noise ratio (SNR) and resolution (R). Our forward model includes Rayleigh scattering, single-layer water clouds with patchy coverage, and pressure-dependent absorption due to water vapor, oxygen, and ozone. We simulate data at R = 70 and R = 140 from 0.4{1.0 µm with SNR = 5; 10; 15; 20 at 550 nm (i.e., for HabEx/LUVOIR-type instruments). At these same SNR, we simulate data for WFIRST paired with a starshade, which includes two photometric points between 0.48{0.6 µm and R = 50 spectroscopy from 0.6{0.97 µm. Given our noise model for WFIRST-type detectors, we find that weak detections of water vapor, ozone, and oxygen can be achieved with observations with at least R = 70 / SNR = 15, or R = 140 / SNR = 10 for improved detections. Meaningful constraints are only achieved with R = 140 / SNR = 20 data. The WFIRST data offer limited diagnostic information, needing at least SNR = 20 to weakly detect gases. Most scenarios place limits on planetary radius, but cannot constrain surface gravity and, thus, planetary mass. 1. INTRODUCTION gets for these studies. Characterization of small, po- The scientific field of exoplanets has been rapidly ad- tentially rocky exoplanets is limited to worlds with cool vancing since the hallmark discovery of the first planet stellar hosts (K and M dwarfs), which offer favorable orbiting a Sun-like star (Mayor & Queloz 1995). Follow- planet-to-star size ratios. Recently, de Wit et al. (2016) ing the launch of NASA's Kepler mission (Borucki et al. studied the combined transmission spectra of two tran- 2003, 2011), the field has seen the discovery of thousands siting Earth-sized planets orbiting the ultracool dwarf of transiting exoplanets and the exciting result that plan- TRAPPIST-1 using the Hubble Space Telescope. While no gas absorption features were detected by de Wit et al. ets with radii between 0:75{2:5 R⊕ are common around solar-type stars (Burke et al. 2015). Only within the last (2016), this work highlights the improvements in signal decade have observational studies for exoplanet atmo- size when terrestrial-sized transiting planets are studied spheric characterization seen substantial development, around low-mass stars. Additionally, since the Habit- starting with the first detection of an exoplanet's atmo- able Zone (Kasting et al. 1993) around a low-mass star sphere by Charbonneau et al. (2002). is relatively close-in, characterization studies of poten- To date, the majority of exoplanet atmospheric char- tially habitable exoplanets around cool stars can benefit acterization investigations have focused on transiting from the frequency of transit events. However, for Sun- worlds. Hot Jupiters, owing to their large sizes and like hosts, the planet-to-star size ratio is much less fa- short orbital periods, are typically emphasized as tar- vorable and the Habitable Zone is located far from the star, thus severely limiting the potential for atmospheric 1 Department of Astronomy & Astrophysics, University of Cal- characterization. ifornia, Santa Cruz, CA 95064, USA Direct, high-contrast imaging has now emerged as an 2 NSF Graduate Research Fellow 3 essential technique for studying the atmospheres of plan- arXiv:1803.06403v1 [astro-ph.EP] 16 Mar 2018 University of California, Santa Cruz, Other Worlds Labora- tory ets at larger orbital separations from their host star (i.e., 4 now at Department of Physics & Astronomy, Northern Ari- at orbital distances & 1 au). Thus far, high-contrast zona University, Flagstaff, AZ 86011, USA 5 imaging has been proven successful in studying atmo- NASA Astrobiology Institute's Virtual Planetary Labora- spheres of young, self-luminous gas giants in the near- tory 6 NASA Ames Research Center, Moffett Field, CA 94035, USA infrared and mid-infrared (e.g., Barman et al. 2011; Ske- 7 Bay Area Environmental Research Institute, Petaluma, CA mer et al. 2014; Macintosh et al. 2015). These worlds, 94952, USA owing to their intrinsic brightness, have typical contrast 8 Space Telescope Science Institute, Baltimore, Maryland −4 21218, USA ratios of 10 with respect to their hosts. A true Jupiter 9 Department of Earth and Planetary Sciences, Johns Hopkins analog at visible wavelengths, by comparison, would have University, Baltimore, MD 21218, USA a contrast ratio of 10−9, while an Earth analog would 10 Department of Physics, Kavli Institute for Particle As- have a contrast ratio of order 10−10. Reflected light in the trophysics and Cosmology, Stanford University, Stanford, CA 94305, USA visible probes to atmospheric depths of up to ∼ 10 bar 11 School of Earth and Space Exploration, Arizona State Uni- for giant planets (Marley et al. 2014), which is compli- versity, Tempe, AZ 85287, USA mentary to the relatively low pressures probed in transit 2 Feng et al. observations (typically less than 10{100 mbar). Addi- dezvous concept. To accomplish this, we must perform tionally, the wavelength range of 0.4{1.0 µm holds rich atmospheric and instrument modeling to simulate the information about a planet's atmosphere, including sig- types of spectra we can expect to measure, and we must natures of methane, water vapor, and haze (Marley et al. develop tools to infer planetary properties from these 2014; Burrows 2014). simulated observations. In spite of the incredible technological challenges, there Traditionally, the comparison to a limited range of for- are multiple planned or in-development space-based mis- ward models has been used to infer atmospheric proper- sions that would be capable of high-contrast imaging ties (such as temperature structure and gas abundances) of exoplanets in reflected light. First among these from spectral observations. This involves iterating to a will be NASA's Wide-Field InfraRed Survey Telescope radiative-convective solution for a given set of planetary (WFIRST, Spergel et al. 2013), which was identified parameters (e.g., gravity, metallicity, equilibrium abun- as the top priority space mission in the 2010 National dances, incident flux), and can include detailed treatment Academy of Sciences Decadal Survey of Astronomy and of aerosols, chemistry, and dynamics within the model Astrophysics12. The WFIRST mission will carry a Coro- atmosphere (Marley & Robinson 2015). The goal is to nagraphic Instrument (CGI) with imaging capability and generate a spectrum that matches available data and, a visible-light integral field spectrograph of wavelength thus, offers one potential explanation for the world's at- resolution ∼ 50 (Noecker et al. 2016; Trauger et al. 2016; mospheric state (e.g., Konopacky et al. 2013; Macintosh Seo et al. 2016; Cady et al. 2016; Balasubramanian et al. et al. 2015; Barman et al. 2015). A more data-driven in- 2016; Groff et al. 2018). Although envisioned primarily terpretation of atmospheric observations is accomplished as a technology demonstrator, it may study the atmo- through inverse modeling, or retrievals. Developed for spheres of relatively cool gas giant exoplanets that have Solar System studies and remote sensing (e.g., Rodgers been previously detected using the radial velocity tech- 1976; Irwin et al. 2008), retrievals have become a valu- nique (Traub et al. 2016). able tool in constraining our understanding of the at- While WFIRST could also have some capability to mospheres of transiting exoplanets. Early exoplanet re- survey stars in the solar neighborhood for lower-mass trieval work invoked grid-based optimization schemes planetary companions (Burrows 2014; Greco & Burrows (Madhusudhan & Seager 2009), while subsequent works 2015; Spergel et al. 2015; Savransky & Garrett 2016; have taken advantage of Bayesian inference with meth- Robinson et al. 2016), it is anticipated that the core op- ods such as optimal estimation and Markov chain Monte tical throughput of the WFIRST CGI will be low for Carlo (MCMC) (e.g., Lee et al. 2012; Benneke & Seager planetary signals. This stems primarily from the com- 2012; Line et al. 2013). plexities of accommodating for WFIRST's on-axis sec- Several studies have examined the hypothetical yield ondary mirror and support structures within the high- from characterizing giant exoplanets observed with a contrast instruments (Traub et al. 2016; Krist et al. space-based coronagraph (such as WFIRST ) with re- 2016). Low throughput drives long requisite integration trieval techniques. Marley et al. (2014), for example, times, thereby likely making spectroscopic observations modeled spectra we could expect from known radial ve- of smaller, less-bright worlds (such as super-Earth exo- locity gas giants if observed by the WFIRST CGI. Given planets) unfeasible except around the very closest stars the diversity of cool giant planets, the model spectra have (Robinson et al. 2016). However, if the WFIRST space- a variety of input assumptions for clouds, surface gravity, craft were to be paired with an external starshade (Cash and atmospheric metallicity. Marley et al. (2014) then 2006; Kasdin et al. 2012), the CGI can be operated in applied retrieval methods to these synthetic spectra, en- a direct mode without coronagraphic masks, substan- abling the exploration of how well atmospheric parame- tially increasing throughput.