arXiv:1709.01219v1 [astro-ph.EP] 5 Sep 2017 [email protected] orsodn uhr hafiDong Chuanfei author: Corresponding ihrta h orsodn au f6.7 of value corresponding the than higher reyepoetermfiain rsn rmtels fwtrin,e significant. ions, water be of to loss accelerat likely the have are from effects arising could ramifications ejections) the mass explore briefly coronal frequent with (e.g., hr h o saertscnreach can rates escape ion the where vn.W eosrt htteinecp aefr(i,wt valu e a with one (ii), (iii) for and rate escape conditions, ion Sun-like wind the a solar that orbiting demonstrate normal We parameters ancient event. Earth-like (ii) with conditions, world magnitude water wind the a quantifying to accurately of model capable is that chemistry, 1 2 5 4 3 eateto srpyia cecs rneo Universi Princeton Sciences, Astrophysical of Department rneo etrfrHlohsc,PictnPam Physi Plasma Princeton Heliophysics, for Center Princeton onA alo colo niern n ple Sciences, Applied and Engineering of School Paulson A. John Cambridge, Astrophysics, Michi for of Center University Harvard-Smithsonian Modeling, Environment Space for Center epeetatreseismlifli H oe (H model MHD multi-fluid three-species a present We yee sn L 2017 using 6, Typeset September version Draft hafiDong, Chuanfei A T E X H EYRTO FWTRWRD I TOPEI LOSSES ATMOSPHERIC VIA WORLDS WATER OF DEHYDRATION THE twocolumn ,2 1, hnun Huang, Zhenguang tl nAASTeX61 in style ∼ 10 × 28 mtv Bhattacharjee Amitava 10 s − 25 1 3 r rca o nesadn o natv,eryslrtp st solar-type early active, an how understanding for crucial are , s aav Lingam, Manasvi − y rneo,N 84,USA 08544, NJ Princeton, ty, ABSTRACT 1 A018 USA 02138, MA sLbrtr,PictnUiest,Pictn J08544, NJ Princeton, University, Princeton Laboratory, cs a,AnAbr I419 USA 48109, MI Arbor, Ann gan, o i.Suiso o ossidcdb pc ete events, weather space by induced losses ion of Studies (i). for avr nvriy abig,M 23,USA 02138, MA Cambridge, University, Harvard + H , 2 O dtedpeino h xpae’ topee We atmosphere. exoplanet’s the of depletion the ed + fwtrinlse rmeolnt.W pl this apply We exoplanets. from losses ion water of n e and f6.0 of e pcal o lnt riigMdaf hr such where M-dwarfs orbiting planets for specially ,5 4, ,2 1, G trfrtrecss i urn omlsolar normal current (i) cases: three for star brT abor ´ − ,edwdwt h eust atmospheric requisite the with endowed ), × tee“arntntp”saeweather space “Carrington-type” xtreme 10 26 oth, ´ s − 1 3 saota re fmagnitude of order an about is , aa Gombosi, Tamas 3 USA and ar 2
1. INTRODUCTION Thus, when modeling atmospheric escape rates, it is The past decades have witnessed the discovery of necessary to take into account stellar wind-exoplanet in- thousands of exoplanets. Much effort has been devoted teractions, and quantify ion escape rates by using sophis- to locating rocky exoplanets in the habitable zone (HZ) ticated magnetohydrodynamic (MHD) models. Such of their host star, as these planets are theoretically ca- studies have been undertaken recently for Proxima Cen- pable of supporting liquid water - an essential ingredient tauri b (Airapetian et al. 2017; Dong et al. 2017b) and for life-as-we-know-it. However, a remarkable feature of the TRAPPIST-1 planets (Dong et al. 2017a). An im- exoplanetary science is that it has opened up the possi- portant to be noted is that these studies implicitly as- bility of discovering planets that are very dissimilar to sumed that the atmospheric composition was similar to those found in our Solar system. that of Mars, Venus or Earth. In this paper, we shall consider a scenario wherein the planet’s atmosphere is A noteworthy example in this category is ocean 1 planets, also referred to as water worlds. These primarily comprised of water vapor (Goldblatt 2015), planets are anticipated to be volatile-rich (Kuchner and determine the corresponding ion escape rates. 2003), and possess oceans that are conceivably hun- 2. MULTI-FLUID MHD MODEL dreds of kilometers deep (L´eger et al. 2004). No- In this Section, we describe the multi-fluid MHD table candidates for water worlds include Gliese 1214 model, endowed with the electron pressure equation, b (Charbonneau et al. 2009), the six planets in or- that is used to simulate the ion escape processes on wa- bit around Kepler-11 (D’Angelo & Bodenheimer 2016), ter worlds. We present the model equations and clarify Kepler-62e and Kepler-62f (Kaltenegger et al. 2013), the physics behind certain terms. and Kepler-22b (Borucki et al. 2012). Recent studies indicate that most super-Earths with radii & 1.6R⊕ are 2.1. The significance of multi-fluid MHD models not typically rocky since they are volatile-rich planets Before embarking on our description of the multi-fluid (Rogers 2015); a stringent upper bound on the radius MHD equations, a brief description of their utility and was also identified by Chen & Kipping (2017). The accuracy is warranted in the context of our Solar sys- habitability of ocean planets has been extensively in- tem. When we refer to a multi-fluid MHD model, it vestigated, such as their interior structure and dynam- must be understood that the continuity, momentum and ics (Nettelmann et al. 2011), mass-radius relationships energy equations for each species are prescribed, unlike (Benneke & Seager 2013), and atmospheric chemistry the better-known multi-species single-fluid MHD model (Kaltenegger et al. 2013; Goldblatt 2015). where only one equation is prescribed for the momentum We also emphasize that terrestrial planets in our Solar and for the energy (Ma et al. 2004). Although multi- system are believed to have possessed primordial oceans, fluid MHD models are computationally more expensive, although the evidence remains disputed. Most notably, they have the advantage of improved realism and accu- there have been several proposals that the Earth was racy (T´oth et al. 2012; Dong et al. 2014). initially an ocean-dominated world (Iizuka et al. 2010; Multi-fluid MHD models have been successfully ap- Arndt & Nisbet 2012), perhaps even up to the late- plied to analyze a wide range of planetary plasma en- Archean era at 2.5 Gya (Flament et al. 2008). In ad- vironments in the Solar system. A few notable ex- dition, many studies have suggested the existence of amples include Mars (Najib et al. 2011; Dong et al. a northern ocean of liquid water on ancient Mars; the 2014, 2015), Earth (Glocer et al. 2009; Brambles et al. reader may consult Wordsworth (2016) for more details. 2010), Europa (Rubin et al. 2015), Enceladus (Jia et al. Venus could also have had a shallow ocean (Way et al. 2012), Ganymede (Paty & Winglee 2004) and comets 2016) until about 1 Gya. Thus, on account of all these (Rubin et al. 2014; Huang et al. 2016a,b). In all of reasons, water-rich exoplanets merit further study. these instances, the obtained theoretical results from One of the most important factors studied in plan- the simulations have been shown to be in very good etary habitability is the long-term existence of an at- agreement with observations from missions such as mosphere. To this end, it is necessary to quantify the MAVEN, Galileo and Rosetta. For instance, in the extent of atmospheric losses from planets. When study- case of Mars, it was shown that the multi-fluid MHD ing the atmospheric escape rates of exoplanets, espe- equations characterized the ion escape rates from the cially in the context of water loss, it has become com- upper atmosphere more accurately as compared to the mon to rely upon models of thermal (Jeans or hydro- multi-species equations (Dong et al. 2014). dynamic) escape. However, as seen from NASA’s Mars Thus, in order to model the effects of water ion losses Atmosphere and Volatile EvolutioN (MAVEN) mission from ocean planets, we shall introduce a novel multi- (Jakosky et al. 2015; Brain et al. 2015), ion escape rates also play an important role; in fact, for planets larger than Earth, neutral losses are not necessarily the domi- 1 Although this assumption is not entirely realistic, it could nant mechanism (Dong et al. 2017b). serve as a reasonable approximation for planets like GJ 1214b that are potentially characterized by water-rich atmospheres (Miller-Ricci Kempton et al. 2012). 3
fluid MHD model which includes two ion fluids and one and η is the magnetic diffusivity, which is defined as electron fluid. The model incorporates all of the key 1 1 1 1 chemical reactions between neutral water molecules and η = = + (7) charged particles and is therefore well-suited for study- µ0σe µ0 σen σei ing the stellar wind-induced ion escape in ocean plan- where the electron conductivity, σe, comprises of contri- ets. The numerical code developed in this paper rep- 2 butions from electron-neutral (σen = e ne/Σn′ νen′ me) resents the first self-consistent model in the literature 2 that is capable of estimating the ion escape rates from and electron-ion (σei = e ne/Σs′ νes′ me) collisions. (exo)planets with water vapor atmospheres. In the above set of equations, note that the source (S) and loss (L) terms for species s associated with photoionization (ν ′ ), electron impact ionization 2.2. The multi-fluid MHD equations ph,s ←→ (νimp,s′ ), charge exchange (kis′ ), and recombination We use ρ, u, p, I , kB and γ to denote the mass (αR,s) are described below: density, velocity vector, pressure, the identity matrix, the Boltzmann constant and the specific heat ratio, re- S ′ ′ ′ ′ spectively. The multi-fluid MHD code simulates three s = msns νph,s + νimp,s + kis ni (8) i=ions ! fluids in total. Of this trio, two of them are ion fluids X + + - water group ions H2O and stellar wind protons H L ′ - that are distinguished via the subscript s. The third s = msns αR,sne + kst nt′ (9) t′=neutrals ! is the electron fluid that is denoted by the subscript e; X note that we also introduce the subscript n for neutral Se = me ns′ (νph,s′ + νimp,s′ ) (10) background species. The multi-fluid MHD equations are s′ X listed as follows (Najib et al. 2011; Rubin et al. 2014): Le = mene αR,sns (11) s=ions X ∂ρs Inelastic collisions between electrons and neutral wa- + ∇ · (ρsus)= Ss − Ls (1) ∂t ter molecules constitute an efficient way of cooling the ∂ (ρsus) ←→ electrons in the planetary lower ionosphere where colli- + ∇ · ρ usus + p I = n q (us − u+) × B ∂t s s s s sions are highly frequent. Hence, we include the overall qsns J B G (rotational, vibrational, and electronic) cooling rate co- + ( × − ∇pe )+ ρs inelastic 3 −1 ene efficient, Ren (in eV cm s ), in (4) by following the prescription provided in Gombosi (2015): +ρs νst(ut − us)+ Ssun − Lsus (2) t=all X inelastic −9 k(Te − Tn) ∂ps Ren =4 × 10 1 − exp − (12) + (us · ∇) p = −γ p (∇ · us) 0.033eV ∂t s s s kTe − 0.10eV ρsνst 2 2 + A 0.415 − exp − + 2kB (Tt − Ts)+ mt (ut − us) 0.10eV ms + mt 3 t=all X −9 SsTn − LsTs 1 2 where A =0 for kTe ≤ 0.188 eV and A =6.5 × 10 for un − us +kB + Ss ( ) (3) kTe > 0.188 eV. ms 3 In order to account for photoionization, we calculate ∂pe + (ue · ∇) pe = −γepe (∇ · ue) the optical depth of the neutral atmosphere by applying ∂t the Chapman functions based on the numerical evalua- ρeνet 2 2 + 2kB (Tt − Te)+ mt (ut − ue) tion given by Smith & Smith (1972). The photoelectron me + mt 3 Eexc t=s,n gains an excess energy ns through the photoionization X process (Huebner et al. 1992) while it loses the ioniza- LeTe 2 exc pot −kB + nn(νph,nE − νimp,nE ) tion energy of H2O during the electron impact ioniza- m 3 ns ns e tion process (Haynes 2014) as indicated in Eq. (4). The 2 inelastic 1 u u 2 electron number density and velocity can be obtained − nennRen + Se ( n − e) (4) 3 3 by assuming quasineutrality (T´oth et al. 2012), ∂B ∇× u × B − J = ( + η ) (5) 1 ∂t n = n q (13) e e s s s=ions where ν is the collision frequency between species, u+ is X the charge-averaged velocity, and in terms of the current, J B u qsnsus ∇× + = (6) u = u+ − = u+ − (14) en e s=ions e ene µ0ene X 4
Table 1. The elastic collision rates have been adopted from Kepler-22b; Earth may also have been an ocean planet Huang et al. (2016a). Zs, ms, ns and Ts denote the charge in the Hadean (or Archean) era, but its atmospheric 3 state, mass (in amu), number density (in cm− ) and tem- composition is expected to have involved other gases. m m perature (in K) of a given species. Here, mst = s t In accordance with our choice, we employ the atmo- ms+mt m T +m T and Tst = s t t s are the reduced mass and tempera- spheric temperature profile of Earth (Schunk & Nagy ms+mt ture. Csn represent numerical coefficients that have been 2009), and assume that its surface pressure is 1 bar. presented in Schunk & Nagy (2009). We also specify a fiducial value for the magnetic dipole moment that is equal to that of Earth, which is a rea- sonable assumption for Earth-like exoplanets. As noted 1 Elastic Collision Rates Value (s− ) earlier, we also assume that water vapor is the primary Z2Z2 m s t √ st nt neutral atmosphere component, and neglect other neu- ion-ion (νii) 1.27 × 3 2 ms T / st tral species for the purposes of this study. The neutral ion-neutral (νin) Csnnn atmospheric background adopted in our model incorpo- n Z2 rates the cold trap for water in the lower altitudes, but ν . s s electron-ion ( ei) 54 5 × 3/2 Te its effects on the ion losses are minimal due to the high- m n Z2 √ e s s ion-electron (νie) 1.27 × 3 2 altitude source of the escaping ions. In future investi- ms T / e gations, we will incorporate water vapor into primor- 5 0.62 electron-neutral (νen) 2.745 × 10− nnTe− dial or current Earth/Mars atmospheres, and study the (ionized) water vapor escape rates for a wider range of exoplanets endowed with other neutral and ion species. + The simulation domain starts at 100 km, where H2O where we have adopted Amp´ere’s law ∇×B = µ0J, with ion density satisfies the photochemical equilibrium con- µ0 denoting the vacuum permeability. dition and float boundary conditions for the velocity u By using (1) - (14), we are now in a position to sim- and the magnetic field B are applied. The simulation ulate the water ions, the stellar wind protons, and the box extends up to 200 planetary radii with a highly non- electrons in a fully self-consistent manner; both the in- uniform spherical grid. The lowest resolution is 10 km dividual behavior of each fluid and interactions between in the planetary lower ionosphere, as this value is small them have been taken into account. The array of pro- enough to capture the scale height variation in the lower cesses included in this model are summarized below. ionosphere and thermosphere. The horizontal (angular) • We have included the effects of both elastic and resolution is 3◦ × 3◦. The rest of the details concerning inelastic collisions. In Table 1, the elastic collision the simulation setup are similar to those described in rates between different fluid species are presented. Dong et al. (2017b). The above set of equations in Sec. 2.2 is solved by us- • The chemical reactions between species are man- ing an upwind finite-volume scheme based on an approx- ifested as source terms on the right hand side, imate Riemann solver to ensure the appropriate conser- including ionization (photoionization and electron vation of plasma variables. Moreover, hyperbolic clean- impact ionization), charge exchange between neu- ing is adopted for divergence control. In order to de- trals and ions, and recombination. The electron termine the steady-state solution, the simulation starts impact ionization rates have been adopted from with a two-stage local time-stepping scheme that enables Cravens et al. (1987). different grid cells to select varying time steps based on • Table 2 summarizes the chemical reactions and the the local stability condition, thereby accelerating the associated rates for inelastic collisions used in the convergence to the steady-state and saving computa- multi-fluid MHD calculations. tional resources. Because of the stiffness of the source terms, a point implicit scheme is used for handling them A detailed discussion of the source terms and the inher- (T´oth et al. 2012). ent chemistry can also be found in Rubin et al. (2014) We will investigate three different cases in this paper. and Huang et al. (2016a). The first utilizes the solar wind conditions at the present Earth, while the other two cases simulate the atmo- 3. SIMULATION SET-UP spheric loss arising from the ancient solar wind around 4 To this date, as noted in Sec. 1, no conclusive ocean Gya, and an extreme “Carrington-type” space weather planets have been detected, although several candidates event. The corresponding stellar wind parameters have have been identified. Consequently, there is a lack of been presented in Table 3. We have chosen to work concrete information regarding their atmospheric prop- with the Planet-Star-Orbital (PSO) coordinates, where erties such as the composition and mass. Thus, for the the X-axis points from the planet toward star, the Z-axis purposes of this preliminary study, we consider the sce- is perpendicular to the planetary orbital plane, and the nario where an Earth-sized ocean planet is orbiting a Y-axis completes the right-hand system. In this study, solar-type star. This situation roughly corresponds to 5
Table 2. Chemical reactions and associated rates for the ocean planets multi-fluid MHD code.
1 Chemical Reaction Rate (s− ) Primary Photolysis and Particle Impact
+ 7 1a H2O + hν → H2O + e− 8.28 × 10− s− Huebner & Mukherjee (2015) + e− + H2O → e− + H2O + e− see text Cravens et al. (1987)
3 1 Ion-Neutral Chemistry Rate (cm s− )
+ + 9 H2O + H2O → H2O+H2O 1.7 × 10− Gombosi et al. (1996) + + 9 H + H2O → H+H2O 1.7 × 10− Gombosi et al. (1996)
3 1 Electron Recombination Chemistry Rate (cm s− )
5 0.569 1.57 × 10− Te− , Te ≤ 800K + 5 0.74 H2O + e− → H2O 4.73 × 10− Te− , 800K
aThe photoionization rate represents the value at a distance of 1 AU from the Sun based on current solar cycle maximum conditions, i.e. the 1 EUV case in Table 3, and must be appropriately rescaled when dealing with early epochs or other stellar systems.
Table 3. Three hypothetical stellar wind input parameters for a solar-type star based on: (i) typical solar wind parameters at 1 AU at the current epoch (Schunk & Nagy 2009), (ii) ancient (4.02 Gyr) solar wind parameters at 1 AU (Boesswetter et al. 2010) and (iii) solar wind values at the maximum total pressure of an extreme “Carrington-type” space weather event (Ngwira et al. 2014). Note that 1 EUV (below) refers to the EUV flux received at Earth during the solar cycle maximum.
3 + 1 nsw (cm− ) vsw (km/s) IMF (nT) Radiation H2O loss rate (s− ) Current 8.7 (-468, 0, 0) (-4.4, 4.4, 0) 1 EUV 6.7×1025 Early 136.7 (-910, 0, 0) (-15.6, 30.2, 0) 12 EUV 6.0×1026 Carrington Event 424.5 (-1937.5, 6.7, -13.0) (0, 23.0, -194.3) 12 EUV 7.3×1027 we assume that the planetary dipole moment is along The ion escape rate reaches 7.3 × 1027 s−1, about one the Z-axis and has the same polarity as Earth. order of magnitude higher than the normal ancient case (6.0 × 1026 s−1). Our analysis suggests that such space 4. RESULTS weather events may prove to be a key driver of atmo- Our primary result concerns the magnitude of the ion spheric losses for (exo)planets orbiting an active young + escape rates for water ions, i.e. for H2O . We find that Sun-like star. + the escape rate for the current solar wind parameters is Fig. 1 presents the contour plots of the H2O ion den- 6.7 × 1025 s−1, while the corresponding value at 4 Gya sity, the sum of the ion and electron thermal pressures is 6.0 × 1026 s−1 as shown in Table 3. Hence, we find and the magnetic field lines for the three cases. The that the numerical escape rate 4 Gya is a factor of 9.0 planetary magnetosphere is significantly compressed higher than the current value. It is quite interesting from left to right due to the enhanced stellar wind dy- 2 that this factor appears to be in good agreement with namic pressure (Pdyn = nswmswvsw). An inspection of + the theoretical value of 8.7 for unmagnetized or weakly Fig. 1 reveals that H2O ions escape from the planet magnetized planets; the theoretical prediction follows along open geomagnetic field lines through the lobes ˙ ˙ −α + from N(t) ∝ M⋆ ∝ ts with α ∼ 1 (Johnstone et al. and magnetotail. The outflow of H2O ions (polar 2015; Lingam & Loeb 2017). wind) is mainly caused by their acceleration (Fig. 2, It is well-known that space weather events, such as first row) resulting from the electron pressure gradient, coronal mass ejections (CMEs), occur frequently on ∇pe (i.e., the ambipolar electric field), the J × B force, young Sun-like stars (Shibayama et al. 2013). Hence, and further pickup by the stellar wind at high altitudes we have also investigated an extreme “Carrington-type” due to the Lorentz force, nsqs (us − u+) × B, in the ion space weather event (Ngwira et al. 2014) to approxi- momentum equation (2); the latter exists only in the mately mimic a typical CME from a stormy young Sun. multi-fluid MHD model since it vanishes for the multi- 6
+ -3 + -3 + -3 n(H2O ) in cm (current) n(H2O ) in cm (early) n(H2O ) in cm (Carrington) 20 4 4 4 10 2 2 2
Z 0 0 0 0
-10 -2 -2 -2
-20 -4 -4 -4 -20 -10 0 10 -20 -10 0 10 -20 -10 0 10 X X X
Pion+Pelectron in nPa (current) Pion+Pelectron in nPa (early) Pion+Pelectron in nPa (Carrington) 20 3 3 3
2 2 2 10 1 1 1
Z 0 0 0 0 -10 -1 -1 -1
-20 -2 -2 -2 -20 -10 0 10 -20 -10 0 10 -20 -10 0 10 X X X
+ Figure 1. The logarithmic scale contour plots of the H2O ion density (first row) and the sum of the ion (Pion) and electron (Pelectron) pressure (second row) with magnetic field lines (in white) in the meridional plane for current solar wind conditions, early solar wind conditions and an extreme “Carrington-type” space weather event. species MHD model (Dong et al. 2014). As seen from planetary and stellar wind parameters equal to that of the second row of Fig. 1, a stronger stellar wind dynamic the current and ancient (4 Gya) Earth-Sun system. + pressure is commensurate with larger plasma thermal We found that the ion escape rate of H2O from this pressure in the magnetosheath, and is consistent with planet at the “current” epoch is 6.7 × 1025 s−1. This the magnetospheric shape. value is slightly higher than, or comparable to, the ion In Fig. 2, the smooth transition of the photoioniza- escape rates observed for the terrestrial planets in our tion rate around the terminator is a result of adopt- Solar system (Lammer 2013). We also demonstrated ing the Chapman function (Smith & Smith 1972). It that the escape rate at 4 Gya was about an order of + is noteworthy that the nightside H2O ions are mainly magnitude higher than the present-day value. For an produced by electron impact ionization, as seen from extreme “Carrington-type” space weather event that is equation (8) and Table 2. The contribution of charge ex- expected to have been frequent in the Hadean epoch, + 28 −1 change with stellar wind protons is concentrated around the H2O ion escape rate can reach ∼ 10 s , about the polar cap regions due to the proton precipitation two orders of magnitude higher than the current Earth along the open magnetic field lines; the charge exchange value. In other words, for active stars (of any spectral + with H2O has not been plotted explicitly. type), space weather events must be taken into account when studying the effects of ion atmospheric losses. At this stage, we wish to reiterate a few caveats. 5. CONCLUSIONS The planetary and stellar parameters chosen herein were We have developed a numerical code capable of based on fiducial values, and must be duly altered when accurately simulating the water ions escaping from studying a given exoplanet-star interaction. Second, the (exo)planets. Our code contains the relevant atmo- atmosphere was assumed to be consist mostly of water spheric chemistry, and evolves each species indepen- vapor, although we intend to generalize our model to in- dently in accordance with the multi-fluid MHD equa- clude other species in the future. Lastly, any statements tions. We applied this model to a water world with concerning atmospheric depletion must be interpreted 7
+ -3 + -3 n(H2O ) in cm (early) n(H2O ) in cm (Carrington) with caution since sources like outgassing and impacts 4 4 4 were not taken into account. 3 For all cases considered herein, we find that Earth-like 24 2 2 oceans (with a total mass of ∼ 10 g) will not be evap-
Z 2 orated over Gyr timescales as the ion escape rates are 0 0 far too low (by 3-4 orders of magnitude). In contrast, 1 for exoplanets in the HZ of M-dwarfs, the situation is -2 -2 0 very different because of the fact that the escape rate -2 -1 0 1 2 -2 -1 0 1 2 approximately scales with the inverse of the distance X X (Lingam & Loeb 2017). Hence, the escape rate (of wa- Phi -3 -1 Phi -3 -1 RH2O+ in cm s (early) RH2O+ in cm s (Carrington) ter ions) will be 2-3 orders of magnitude higher than 4 4 4 on Earth even without considering space weather ef- 3 3 3 fects. Thus, planets orbiting these stars could have their 2 2 oceans depleted over Gyr timescales, especially when ex-
Z 2 1 1 treme space weather events are taken into account. 0 0 Rapid desiccation would have important consequences 1 -1 -1 for the evolution of life on such planets, especially with 0 -2 -2 regards to oceanic and coastal biodiversity, productiv- -2 -1 0 1 2 -2 -1 0 1 2 ity, and food webs. Another important consequence of X X desiccation is that water has been identified as one of Imp -3 -1 Imp -3 -1 RH2O+ in cm s (early) RH2O+ in cm s (Carrington) the crucial ingredients for the functioning of plate tec- 4 4 4 tonics. If plate tectonics were to shut down, a wide 3 3 3 array of deleterious effects would follow (Lammer et al. 2 2 2009). One of the most notable among them is that
Z 2 1 1 the planetary dynamo could shut down, and unmagne- 0 0 tized planets are more vulnerable to sources of ionizing 1 -1 -1 radiation and atmospheric stripping by the stellar wind. 0 -2 -2 Hence, we conclude by observing that future studies -2 -1 0 1 2 -2 -1 0 1 2 of planetary habitability should take the effects of water X X ion losses into account. CX -3 -1 CX -3 -1 RH2O+ in cm s (early) RH2O+ in cm s (Carrington) 4 4 4 3 3 3 2 2
Z 2 1 1 0 0 1 -1 -1 0 -2 -2 The authors acknowledge fruitful discussions with -2 -1 0 1 2 -2 -1 0 1 2 X X Adam Burrows and Abraham Loeb. CD is sup- ported by the NASA Living With a Star Jack Eddy Figure 2. First row: The logarithmic scale contour plots Postdoctoral Fellowship Program, administered by + + of the H2O ion density with H2O ion velocity vectors the University Corporation for Atmospheric Research. (in black) and magnetic field lines (in white) in the merid- ML was supported by the Institute for Theory and ional plane for early solar wind conditions and an extreme Computation (ITC) at Harvard University. Re- “Carrington-type” space weather event. Note that the black sources for this work were provided by the NASA + arrows depict both the direction and the magnitude of H2O High-End Computing (HEC) Program through the + ion velocities. The acceleration of H2O ions is demon- NASA Advanced Supercomputing (NAS) Division strated by the variations in the vector magnitude with alti- at Ames Research Center. The Space Weather tudes. Second to fourth rows: The logarithmic scale contour Modeling Framework that contains the BATS-R-US Phi code used in this study is publicly available from plots of the photoionization rate RH2O+, electron impact ion- Imp CX http://csem.engin.umich.edu/tools/swmf. For dis- ization rate RH2O+, and charge exchange rate RH2O+ (with + tribution of the model results used in this study, please stellar wind protons) of H2O ions. contact the corresponding author.
REFERENCES Airapetian, V. S., Glocer, A., Khazanov, G. V., et al. 2017, Arndt, N. T., & Nisbet, E. G. 2012, AREPS, 40, 521 ApJL, 836, L3 8
Benneke, B., & Seager, S. 2013, ApJ, 778, 153 Jia, Y.-D., Ma, Y. J., Russell, C. T., et al. 2012, Icarus, Boesswetter, A., Lammer, H., Kulikov, Y., Motschmann, 218, 895 U., & Simon, S. 2010, Planet. Space Sci., 58, 2031 Johnstone, C. P., G¨udel, M., Brott, I., & L¨uftinger, T. Borucki, W. J., Koch, D. G., Batalha, N., et al. 2012, ApJ, 2015, A&A, 577, A28 745, 120 Kaltenegger, L., Sasselov, D., & Rugheimer, S. 2013, ApJL, Brain, D. A., McFadden, J. P., Halekas, J. S., et al. 2015, 775, L47 Geophys. Res. Lett., 42, 9142 Kuchner, M. J. 2003, ApJL, 596, L105 Brambles, O. J., Lotko, W., Damiano, P. A., et al. 2010, J. Lammer, H. 2013, Origin and Evolution of Planetary Geophys. Res. A, 115, A00J04 Atmospheres: Implications for Habitability, Springer Charbonneau, D., Berta, Z. K., Irwin, J., et al. 2009, Briefs in Astronomy (Springer), Nature, 462, 891 doi:10.1007/978-3-642-32087-3 Chen, J., & Kipping, D. 2017, ApJ, 834, 17 Lammer, H., Bredeh¨oft, J. H., Coustenis, A., et al. 2009, Cravens, T. E., Kozyra, J. U., Nagy, A. F., Gombosi, T. I., A&ARv, 17, 181 & Kurtz, M. 1987, J. Geophys. Res., 92, 7341 D’Angelo, G., & Bodenheimer, P. 2016, ApJ, 828, 33 L´eger, A., Selsis, F., Sotin, C., et al. 2004, Icarus, 169, 499 Dong, C., Bougher, S. W., Ma, Y., et al. 2014, GeoRL, 41, Lingam, M., & Loeb, A. 2017, Int. J. Astrobiol., 2708 arXiv:1707.02996 Dong, C., Jin, M., Lingam, M., et al. 2017a, Ma, Y., Nagy, A. F., Sokolov, I. V., & Hansen, K. C. 2004, arXiv:1705.05535, arXiv:1705.05535 JGRA, 109, A07211 Dong, C., Lingam, M., Ma, Y., & Cohen, O. 2017b, ApJL, Miller-Ricci Kempton, E., Zahnle, K., & Fortney, J. J. 837, L26 2012, ApJ, 745, 3 Dong, C., Ma, Y., Bougher, S. W., et al. 2015, Najib, D., Nagy, A. F., T´oth, G., & Ma, Y. 2011, J. Geophys. Res. Lett., 42, 9103 Geophys. Res. A, 116, A05204 Flament, N., Coltice, N., & Rey, P. F. 2008, E&PSL, 275, Nettelmann, N., Fortney, J. J., Kramm, U., & Redmer, R. 326 2011, ApJ, 733, 2 Glocer, A., T´oth, G., Ma, Y., et al. 2009, J. Geophys. Res. Ngwira, C. M., Pulkkinen, A., Kuznetsova, M. M., & A, 114, A12203 Glocer, A. 2014, J. Geophys. Res. A, 119, 4456 Goldblatt, C. 2015, Astrobiology, 15, 362 Paty, C., & Winglee, R. 2004, GeoRL, 31, L24806 Gombosi, T. I. 2015, in Washington DC American Rogers, L. A. 2015, ApJ, 801, 41 Geophysical Union Geophysical Monograph Series, Vol. Rubin, M., Combi, M. R., Daldorff, L. K. S., et al. 2014, 207, Magnetotails in the Solar System, ed. A. Keiling, ApJ, 781, 86 C. M. Jackman, & P. A. Delamere, 169–188 Gombosi, T. I., De Zeeuw, D. L., H¨aberli, R. M., & Powell, Rubin, M., Jia, X., Altwegg, K., et al. 2015, J. Geophys. K. G. 1996, J. Geophys. Res., 101, 15233 Res. A, 120, 3503 Haynes, W. M. 2014, CRC Handbook of Chemistry and Schunk, R., & Nagy, A. 2009, Ionospheres: Physics, Plasma Physics, 95th edn. (CRC Press) Physics, and Chemistry, Cambridge Atmospheric and Huang, Z., T´oth, G., Gombosi, T. I., et al. 2016a, J. Space Science Series (Cambridge Univ. Press), Geophys. Res. A, 121, 4247 doi:10.1017/CBO9780511635342 —. 2016b, MNRAS, 462, S468 Shibayama, T., Maehara, H., Notsu, S., et al. 2013, ApJS, Huebner, W. F., Keady, J. J., & Lyon, S. P. 1992, Ap&SS, 209, 5 195, 1 Smith, F. L., & Smith, C. 1972, J. Geophys. Res., 77, 3592 Huebner, W. F., & Mukherjee, J. 2015, Planet. Space Sci., T´oth, G., van der Holst, B., Sokolov, I. V., et al. 2012, J. 106, 11 Comp. Phys., 231, 870 Iizuka, T., Komiya, T., Rino, S., Maruyama, S., & Hirata, Way, M. J., Del Genio, A. D., Kiang, N. Y., et al. 2016, T. 2010, GeCoA, 74, 2450 GeoRL, 43, 8376 Jakosky, B. M., Grebowsky, J. M., Luhmann, J. G., et al. Wordsworth, R. D. 2016, AREPS, 44, 381 2015, Science, 350, 0210