Draft version January 3, 2018 Preprint typeset using LATEX style emulateapj v. 12/16/11 HABITABILITY OF EXOPLANET WATERWORLDS Edwin S. Kite Department of Geophysical Sciences, University of Chicago, Chicago, IL 60637, USA [email protected] Eric B. Ford Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA Center for Exoplanets and Habitable Worlds, The Pennsylvania State University, University Park, PA 16802, USA Institute for CyberScience and Pennsylvania State Astrobiology Research Center Draft version January 3, 2018 ABSTRACT We model the evolution of ocean temperature and chemistry for rocky exoplanets with 10-1000× Earth's H2O but without H2, taking into account C partitioning, high-pressure ice phases, and atmosphere-lithosphere exchange. Within our model, for Sunlike stars, we find that: (1) habitability is strongly affected by ocean chemistry; (2) possible ocean pH spans a wide range; (3) exsolution-driven climate instabilities are possible; (4) surprisingly, many waterworlds stay habitable for >1 Gyr, and (contrary to previous claims) this longevity does not necessarily involve geochemical cycling. We also find, using an ensemble of N-body simulations that include volatile loss during giant impacts, that a substantial fraction of habitable-zone rocky planets emerge after the giant impact era with deep, ice-free water envelopes. This outcome is sensitive to our assumptions of low initial 26 60 abundances of Al and Fe in protoplanetary disks, plus H2-free accretion. We use the output of the N-body simulations as input to our waterworld evolution code. Thus, for the first time in an an end- to-end calculation, we show that chance variation of initial conditions, with no need for geochemical cycling, can yield multi-Gyr habitability on waterworlds. Keywords: planets and satellites: terrestrial planets | planets and satellites: physical evolution | planets and satellites: individual (Kepler-452b,Kepler-1638b, Kepler-1606b, Kepler-1090b, Kepler-22b, τ Ceti e, Proxima Cen b, TRAPPIST-1, GJ 667 C, LHS 1140 b, Ross 128 b, Kepler-62f, Kepler-186f, GJ 832 c, HD 40307 g, Kepler-442b, Kepler-1229 b) 1. INTRODUCTION. (Foley 2015; Abbot et al. 2012). Yet waterworlds should Habitable-zone small-radius exoplanets are common be common in the Galaxy (e.g. Mulders et al. 2015). Are (Burke et al. 2015; Dressing & Charbonneau 2015). waterworlds doomed? Not necessarily, because there is What fraction are habitable? For the purposes of this another (less-studied) road to long-term planetary hab- paper, an operational definition of a potentially habit- itability: able exoplanet is: maintains T < 450K liquid water on Cycle-independent planetary habitability: Planets its surface continuously for timescales that are relevant 7 where pCO2, ocean depth, and surface temperature Tsurf for biological macroevolution, 10 yr (Vermeij 2006; 7 1 remain within the habitable range for 10 yr without Carter 1983; Bains et al. 2015). With this definition, geochemical cycling. two pathways allow long-term planetary habitability: Habitability maintained by geochemical cy- In this paper, we first show that the long-term climate cles: Planets that maintain habitability despite large- evolution of waterworlds can be modeled independently of geochemical cycles (x2). Next, we set up (x3) and arXiv:1801.00748v1 [astro-ph.EP] 2 Jan 2018 amplitude perturbations (due to tectonics, stellar evolu- tion, e.t.c.) through a negative feedback. A proposed run (x4) such a waterworld evolution model. Using the mechanism for the negative feedback is a geochemical model, we estimate what combinations of water abun- cycle involving volcanic outgassing and the fixation of dance, initial carbon abundance, and geologic processes atmospheric gases into rocks (carbonate-silicate weath- allow habitable surface water to persist for >1 Gyr. We ering feedback; Walker et al. 1981; Kasting et al. 1993). focus on CO2+H2O(±N2) atmospheres (Wordsworth & Pierrehumbert 2013). We ignore H2 warming (Stevenson Climate-stabilizing geochemical cycles are a mainstay 1999). We emphasize long-term climate evolution, and of textbooks and review papers (Catling & Kasting defer discussion of nutrients to x6. In x5, we use an N- 2017; Knoll et al. 2012; Kaltenegger 2017, and refer- body model of planet assembly to demonstrate how wa- ences therein). However, such feedbacks probably do not terworlds form and migrate to the Habitable Zone (HZ). work on habitable-zone rocky planets with water mass- We discuss in x6 and conclude in x7. Figs. 15{17 show fractions 10× - 103× that of the Earth { \waterworlds" our main results. 1 Sub-ice oceans in extrasolar planetary systems may be habit- able, but this cannot be confirmed from Earth by remote sensing. 2 Kite & Ford 1.1. This paper in context Simulations of rocky-planet formation yield many rocky small-radius (Rpl <1.6 R⊕) planets with water 2 mass fractions fW = 10-1000 × fW;⊕ (Raymond et al. 2004, 2007; Ciesla et al. 2015; Mulders et al. 2015). −4 Earth's fW (∼ 10 ) appears to be the result of chance (e.g., Raymond et al. 2007; Morbidelli et al. 2016; Licht- enberg et al. 2016) (x5.2). Indeed, simulations suggest fW ≈ fW;⊕ may be uncommon on planets in general (e.g. Tian & Ida 2015). Among previous studies of waterworlds (e.g., Kuchner 2003; L´egeret al. 2004; Selsis et al. 2007; Fu et al. 2010; Abbot et al. 2012; Levi et al. 2017), the closest in in- tent to our own are those of Kitzmann et al. (2015) and Noack et al. (2016). Kitzmann et al. (2015) consider carbonate-system equilibria and find a \sweet-spot" in total planet C similar to our own, but they do not con- sider geological processes. Noack et al. (2016) empha- size the development of high-pressure H2O-ice (HP-ice) layers that isolate the liquid ocean from the nutrients supplied by silicate-rock leaching. Noack et al. (2016) also confirm the result of Kite et al. (2009) that deep oceans suppress silicate volcanism. We conservatively do not count planets with HP-ice as examples of habitability Figure 1. Sketch showing response of an Earth-mass planet to (see x6.4 for discussion). For Tsurf . 375K, this sets a perturbations of the geologic C cycle for (upper panel) an planet Tsurf -dependent upper limit on the depth of a habitable with shallow oceans but no geochemical feedbacks versus (lower ocean (Fig. 1.1), and we consider only <8 GPa oceans. panels) a waterworld. The two temperature-versus-time tracks for We go beyond Noack et al. (2016) by considering the each panel correspond to an increase (upper track) and a decrease effect of C on climate, by using an N-body code to cal- (lower track) in the CO2 flux from the lithosphere into the atmo- sphere+ocean. Color shows fraction of initial H2O lost to space, culate how volatile delivery and giant impacts \load the assuming XUV-energy-limited escape Lammer et al. (2009) when dice" by regulating the fraction of planets that can have the dimensionless moist number >1 (Wordsworth & Pierrehum- cycle-independent planetary habitability, and { most im- bert 2013). In the shallow-ocean case, the climate rapidly becomes uninhabitable (so in this toy model geochemical feedbacks are re- portantly { by tracking Tsurf and pCO2 for 10 Ga on the quired to restore and sustain habitability). In the deep-ocean case, habitable worlds we model. the planet remains habitable for many Gyr due to the subdued Most of the physical and chemical processes we dis- geologic C-cycle, plus ocean dilution. cuss have been investigated previously in an Earth con- text. The novel aspect of our paper is the application to exoplanet waterworlds (x2). 2. HOW TO MODEL HABITABILITY ON WATERWORLDS. Waterworlds are buffered against H2O loss (x2.1), but because CO2 also modulates planetary habitability (x2.2), we need to consider the chemistry that sets ocean pH and thus CO2 (x2.3). Fortunately, nature makes this task easier on waterworlds (x2.4), by chemically isolating the ocean+atmosphere from the lithosphere after 108 yr. This isolation permits cycle-independent planetary hab- itability. 2.1. Waterworlds are buffered against H2O loss For an initial water endowment of >50 Earth oceans in the habitable zone (HZ) of a Sunlike star, the moist greenhouse cannot cause complete ocean loss because XUV-limited water loss in the G-star habitable zone can- Figure 2. Cartoon showing the effect of partial pressure of not remove the bulk of the ocean (Kasting 1988; Lammer atmospheric CO2 and insolation L∗ on habitability, based on et al. 2009; Zahnle & Catling 2017) (Fig. 1). The XUV Wordsworth & Pierrehumbert (2013). The planet fates define a limit may be too pessimistic, in part because photolysis lozenge of habitability that is broadest for O(1) bar pCO2. Color can be limiting for high UV fluxes (e.g. Wordsworth ramp corresponds to Tsurf (K), and extends from 273K to 500K. The peak around 8 bars is an interpolation artifact. & Pierrehumbert 2013; Owen & Alvarez 2016. Also, 2 The H mass fraction of the atmosphere+ocean+crust+(silicate mantle), excluding any H in the core, multiplied by 9 to get H2O- equivalent mass fraction. Habitability of exoplanet waterworlds 3 atmosphere-interior exchange of H2O is not important for deep-ocean planets that have solid rock interiors, be- cause the H-storage capacity of solid Earth mantle rocks is only 1-10 Earth oceans (e.g. Hirschmann 2016) { this is swamped by fW on a waterworld (Tikoo & Elkins- Tanton 2017; for a contrary view, see Marty 2012). For such worlds, the ocean you have after cooling from the last giant impact is the ocean you keep. 2.2. Shallow-ocean worlds: vulnerable unless geochemical cycles maintain habitability. Planetary habitability is modulated by CO2 (Fig. 1). Very high partial pressures of atmospheric CO2 (pCO2) lead to temperatures too high for life, but intermediate pCO2 can stave off glaciation (up to an outer limit set by CO2 condensation and Rayleigh scattering) and thus ex- tend habitability (Fig.