MNRAS 000,1–9 (2020) Preprint 6 July 2020 Compiled using MNRAS LATEX style file v3.0
Losing Oceans: The Effects of Composition on the Thermal Component of Impact-driven Atmospheric Loss
John B. Biersteker1? and Hilke E. Schlichting1,2 1Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139-4307, USA 2UCLA, 595 Charles E. Young Drive East, Los Angeles, CA 90095, USA
Accepted XXX. Received YYY; in original form ZZZ
ABSTRACT The formation of the solar system’s terrestrial planets concluded with a period of gi- ant impacts. Previous works examining the volatile loss caused by the impact shock in the moon-forming impact find atmospheric losses of at most 20–30 per cent and essentially no loss of oceans. However, giant impacts also result in thermal heating, which can lead to significant atmospheric escape via a Parker-type wind. Here we show that H2O and other high-mean molecular weight outgassed species can be effi- ciently lost through this thermal wind if present in a hydrogen-dominated atmosphere, substantially altering the final volatile inventory of terrestrial planets. Specifically, we demonstrate that a giant impact of a Mars-sized embryo with a proto-Earth can re- move several Earth oceans’ worth of H2O, and other heavier volatile species, together with a primordial hydrogen-dominated atmosphere. These results may offer an expla- nation for the observed depletion in Earth’s light noble gas budget and for its depleted xenon inventory, which suggest that Earth underwent significant atmospheric loss by the end of its accretion. Because planetary embryos are massive enough to accrete primordial hydrogen envelopes and because giant impacts are stochastic and occur concurrently with other early atmospheric evolutionary processes, our results suggest a wide diversity in terrestrial planet volatile budgets. Key words: planets and satellites: atmospheres – planets and satellites: formation – planets and satellites: terrestrial planets
1 INTRODUCTION .60 Myr after the formation of the solar system (Canup & Asphaug 2001; Asphaug 2014; Barboni et al. 2017). The growth of rocky planets in the solar system is believed to have culminated in a series of giant impacts where plan- The atmospheres of young terrestrial planets evolve sig- etary embryos, bodies which are thousands of kilometers in nificantly over their first 100–200 Myr. Isotopic evidence size, collide and merge over a timescale of ∼10–100 million indicates that the formation timescale for Mars is compa- years to form the terrestrial planets (Morbidelli et al. 2012). rable to the lifetime of the gas disk (Dauphas & Pour- This giant impact phase is a natural outcome of the growth mand 2011), suggesting that embryos that formed the ter- of planetary embryos in the stabilizing presence of a gas restrial planets grew to at least the mass of Mars (∼0.1 M⊕) arXiv:2007.01446v1 [astro-ph.EP] 3 Jul 2020 disk. While the gas is present, interactions between the em- prior to the dissipation of nebular gas. Mars-sized and bryos and the gas damp the embryosˆaA˘Z´ orbital eccentric- larger embryos are massive enough to directly accrete a ities. When the gas is removed, gravitational interactions primary atmosphere of hydrogen and helium (H/He) from between the planetary embryos lead to excitation of their the disk (Hayashi et al. 1979; Lammer et al. 2014; Inam- orbital eccentricities, resulting in orbit crossing and colli- dar & Schlichting 2015). The noble gas composition in the sions (e.g., Chambers & Wetherill 1998; Walsh & Levison Earth’s mantle may reflect the dissolution of nebular gas 2019). This period of instability has an estimated duration from such a primordial atmosphere into the proto-Earth’s of ∼30−100 Myr at ∼1 au (Morbidelli et al. 2012; Izidoro et al. magma ocean (Harper & Jacobsen 1996; Mukhopadhyay 2015; Walsh & Levison 2019). In the case of Earth, bulk ac- 2012; Williams & Mukhopadhyay 2019). In addition, the cretion was likely completed with the moon-forming impact ubiquity of sub-Neptunes–exoplanets residing within 1 au of their host stars that are shrouded in hydrogen-dominated envelopes containing roughly 1 per cent of the planet’s to- tal mass–suggests that rapid core formation and gas accre- ? Email: [email protected] tion from the primordial disk is commonplace in planet for-
© 2020 The Authors 2 J. B. Biersteker and H. E. Schlichting mation (Fressin et al. 2013; Ginzburg et al. 2016). In some mixed atmospheres, consisting of primordial and outgassed planet formation models (O’Brien et al. 2014; Johansen & material, and on secondary atmospheres, which are com- Lambrechts 2017), the embryos that eventually form the ter- posed of heavier volatile species. In Section2 we review the restrial planets can grow to &50 per cent of their final size model of atmospheric evolution and loss following an im- before the dissipation of the gas disk. It is therefore pos- pact and present resulting analytical estimates for condi- sible that planetary embryos could have not only accreted tions where substantial atmospheric loss can be expected. but also retained significant H/He envelopes over 10–100 Myr In Section3, we present numerical simulations for four case timescales, well into the giant impact phase (Lammer et al. studies spanning a range from pure H/He envelopes to com- 2014). pletely secondary steam atmospheres. We close in Section4 In addition to directly accreting an atmosphere, as a with a discussion. planet cools and its magma ocean solidifies, volatiles that were dissolved in the magma ocean can be outgassed. The magma oceans of terrestrial planets without substantial pri- 2 ATMOSPHERIC EVOLUTION MODEL mary atmospheres can cool quickly (1–10 Myr), likely pro- ducing an episode of catastrophic outgassing near the end of We calculate the extent of atmospheric loss from planetary the magma ocean’s crystallization that releases a secondary winds following an impact using the model developed in atmosphere dominated by water vapor (H2O) and carbon Biersteker & Schlichting(2019), adapted to handle the di- dioxide (CO2)(Elkins-Tanton 2008, 2012). This steam atmo- verse atmospheric compositions and the range of impact sce- sphere can persist for ∼10 Myr before cooling and condensing narios expected in terrestrial planet formation. This model (Elkins-Tanton 2008, 2011). If instead the magma ocean is assumes a two-layer atmosphere: the lower layer is opti- blanketed by an optically thick hydrogen-dominated enve- cally thick and convecting, as expected for envelopes with lope, the added insulation will keep surface temperatures the masses and compositions considered here, and the up- high, allowing a magma ocean to persist through the giant per layer is radiative and approximately isothermal (Rafikov impact phase. Chemical equilibration between the long-lived 2006; Piso & Youdin 2014; Ginzburg et al. 2016). The magma ocean and the atmosphere can produce mixed atmo- boundary between these atmospheric layers is defined as the spheres of predominantly H/He and H2O, with the relative radiative-convective boundary (RCB) where r = Rrcb denotes abundances determined by the redox state of the magma the radial distance from centre of the planet to the radiative- ocean (Kite et al. 2020). convective boundary. Because these planets are young and The young atmospheres of terrestrial planets are sub- possess insulating envelopes, we assume a molten surface jected to loss processes beyond those associated with im- that efficiently exchanges heat with the atmosphere. Atmo- pacts. In particular, stellar radiation-powered hydrodynamic spheric mass loss occurs through hydrodynamic outflow be- escape of atmospheric hydrogen, either nebular in origin or yond the outer radius (Rout), which is the smaller of the 1/3 produced through photodissociation of H2O, is thought to Hill radius (RH = a(Mp/3Mstar) ) and the Bondi radius 2 partially explain the observed variation in volatile abun- (RB = 2GMp/cs ), where Mp is the planet mass, Mstar is the dance and noble gas mass fractionation in the solar system’s mass of the host star, a is the planet’s orbital radius, cs is the terrestrial planets (Zahnle & Kasting 1986; Odert et al. 2018; isothermal sound speed of the gas, and G is the gravitational Lammer et al. 2018). Depending on the assumed high energy constant. The mass loss rate is limited by either the density flux from the host star, primary and secondary atmospheres of gas at the outer radius, or the rate at which gas can be can be significantly eroded or even lost over 10–100 Myr delivered from the convective region to the outer radius. Fol- timescales (Johnstone et al. 2015; Lammer et al. 2014; Odert lowing a giant impact, some of the energy deposited by the et al. 2018). Given the timing and duration of these key pro- impactor heats the planet, causing the envelope to thermally cesses driving atmospheric evolution, the giant impact phase expand, leading to accelerated atmospheric escape. likely spans a wide diversity of atmospheric conditions on young terrestrial planets; an early impact may occur on a planet with a thick H/He envelope, while an impactor arriv- 2.1 Atmospheric structure ing near the end of the giant impact phase may encounter a We model the inner convective region of the envelope with much less massive secondary atmosphere. an adiabatic profile defined by Giant impacts have the potential to dramatically alter 1 the atmosphere of an accreting planet. Past work examin- ρ γ − 1 Rp γ−1 ing volatile loss in giant impacts focused on the loss trig- = Λ − 1 + 1 , (1) ρp γ r gered directly by the impact shock (Genda & Abe 2003, γ 2005; Schlichting et al. 2015). Here, we determine the at- P γ − 1 Rp γ−1 = Λ − 1 + 1 , (2) mospheric loss as a result of the thermal heating caused by Pp γ r the giant impact, which can drive atmospheric escape via a T γ − 1 Rp Parker-type wind (Pepin 1997). It has already been demon- = Λ − 1 + 1, (3) T γ r strated that this thermal loss can lead to complete deple- p tion of primordial H/He envelopes during the formation of where Rp is the planet’s radius, γ is the adiabatic index of super-Earths and sub-Neptunes and that, for such planets, the atmosphere, and ρp, Pp, and Tp are the density, pressure, this thermal loss component exceeds that due to the impact and temperature at the surface of the planet. We define Λ ≡ shock by an order of magnitude (Biersteker & Schlichting GMp µmp/(Rp kBTp), where kB is the Boltzmann constant, 2019). Here, we examine the effect of atmospheric composi- mp is the proton mass, and µ is the mean molecular weight of tion on this impact-triggered thermal loss, with a focus on the atmosphere. Beyond the RCB, we model the atmosphere
MNRAS 000,1–9 (2020) Effects of composition on impact-driven atmospheric loss 3 as isothermal, with a temperature given by the equilibrium p 10 temperature, Teq ' Tstar Rstar/2a, where Rstar and Tstar are the radius and effective temperature of the host star. The density profile in this region is exponential, 0Myr R R 5 ρ = ρ exp rcb rcb − 1 , (4) rcb h r RCB p 2 / where h = kBTeqRrcb GMp µ is the scale height. For the atmo- 3
r/R 1% H/He spheres considered here, the majority of the mass and energy Pre-impact 10 Myr is contained in the convective region of the envelope. These quantities can be obtained by integrating over the adiabatic 2 profile:
¹ Rrcb 3 2 0.1% H2O (steam) Menv = 4πRp ρp ρr dr (5) R 1 p 12 10 8 6 4 2 10 10 10 10 10 10 ¹ Rrcb 2 ⇢ [g cm 3] Eg = −4πGMp ρp Rp ρrdr (6) Rp ¹ Rrcb 3 kBTp 2 Eth = 4πρp Rp Tr dr. (7) Figure 1. Examples of the evolution of atmospheric structure µmp(γ − 1) R p following an impact. Atmospheric density (ρ) at different times Because we are interested in young (.100 Myr) planets is shown as a function of radius (r/Rp ) for a H/He atmosphere that are insulated by optically thick envelopes, we make the with a mass fraction of 1 per cent (blue) and an H2O (steam) simplifying assumption that the solid planet can be treated atmosphere with a mass fraction of 0.1 per cent (orange). Three as fully molten and isothermal. We additionally assume that times are shown: pre-impact (solid), soon after the impact and thermal re-equilibration of the planet (dashed), and 10 Myr af- the molten surface efficiently exchanges heat with the over- ter impact (dash-dotted). Circles mark the radiative-convective lying envelope so that the planet’s temperature is roughly boundary (RCB), where the density profile transitions from adia- equal to the base temperature of the envelope, Tp. The ther- batic to exponential. The thermal energy deposited by the impact mal energy of the planet is then E ∼ cV,p MpTp where we causes the envelope to expand, increasing the density at high al- 6 −1 −1 adopt cV,p ∼ 7.5 × 10 erg g K as the approximate heat titudes by orders of magnitude and promoting atmospheric loss. capacity of the solid planet (Alf`eet al. 2001; Guillot et al. As mass is lost and the envelope cools, the envelope contracts and 1995). mass loss is quenched. The steam atmosphere, due to its higher The atmospheric mass loss rate is determined by the mean molecular weight, experiences only minor inflation and no atmospheric profile and is limited either by the density of gas appreciable loss, while the hydrogen-dominated envelope expands dramatically and is reduced to ∼10 per cent of its original mass. at the outer radius or the rate at which gas can be resupplied to this radius from the RCB. The energy required to lift gas from Rrcb to Rout is provided by a combination of the Û Û Û then given by E = Eenv,m + Eenv,L. We numerically integrate cooling luminosity of the planet (Lrcb) and incoming stellar these equations over 2 Gyr using the odeint routine provided radiation, yielding an additional limit on the mass loss rate. by scipy to determine the final atmospheric structure and Combining these conditions yields cumulative atmospheric mass loss for a given planet. The Û 2 LrcbRrcb model described above assumes that the atmosphere remains M = min 4πRout ρoutcs, , (8) GMp collisional at altitudes below Rout. We check this condition is satisfied at each time step. We do this by comparing the p where ρout is the density at Rout, cs = γkBT/µmp is the altitude of the exobase Rexo to Rout, where the exobase is sound speed, and the luminosity at the RCB is conservatively defined as the location where h (from Equa- 3 tion (4)) is equal to the mean free path of a gas molecule. γ − 1 64πσTrcbGMp µmp Lrcb = , (9) We find that halting atmospheric mass loss once the flow is γ 3κR ρrcbkB no longer collisional has only a minor effect on the overall where κR is the Rosseland mean opacity. We adopt a Rosse- results. This is because the flow typically becomes collision- 2 −1 land mean opacity of κR = 0.1 cm g for pure H/He and less when either the atmosphere has already been nearly 2 −1 mixed H/He+H2O atmospheres, and κR = 0.05 cm g when completely lost, or when our model already predicts limited considering H2O atmospheres as appropriate approximations atmospheric mass loss. for the conditions at the RCB (Freedman et al. 2008; Valen- The initial conditions for our calculations of atmo- cia et al. 2013). The energy lost through atmospheric mass spheric evolution are determined by the pre-impact state of Û Û loss is then Eenv,m ≈ GMp M/Rrcb, while the radiative cooling the planet and the chosen impactor mass and velocity. We Û is given by Eenv,L = −Lrcb. set an initial temperature for the base of the atmosphere of We combine these models of the planet and its atmo- Tp,0 = 2000 K, consistent with our assumption of a molten sphere to numerically calculate the thermal evolution and planetary surface and expectations for young planets hosting mass loss of the planetary atmosphere following an impact substantial insulating atmospheres. We assume the impact 0 (see Figure1). The total energy, E = Ep + Eg + Eth uniquely results in perfect accretion, so that Mp = Mp + Mimp, where determines the atmospheric structure, which in turn defines the prime indicates the post-impact value. The energy added ' 2 / the mass loss rate (Equation (8)). The change in energy is by the impact is Eimp ηMimpvimp 2, where Mimp and vimp
MNRAS 000,1–9 (2020) 4 J. B. Biersteker and H. E. Schlichting
are the impactor mass and velocity, and η describes the frac- tion of the impact energy that is used to heat the envelope 16 and the core (where the core refers to the iron and silicate portion of the planet). In our calculations we assume η ∼ 0.5 =4/3 √
2006, 2014; Carter et al. 2020). max 14 µ =7/5
MNRAS 000,1–9 (2020) Effects of composition on impact-driven atmospheric loss 5 mosphere with mass fractions of fH/He = MH/He/Mp ∼ 1 per for the fractionation factor (Equation (14)) and for mass Û Û Û cent and fH2O ∼ 0.1 per cent would contain the equivalent of conservation (M = Mi + Mj ). In the supply-limited regime, ∼4 Earth oceans of water, but have a mean molecular weight the hydrogen flux, taken to be species i, in the fractionation Û 2 of µ ' 2.5. If such an atmosphere is well-mixed to the outer factor is determined by Mi = 4πRout ρi,out. The density of radius, the above calculations would indicate it could be lost hydrogen is calculated from in a giant impact even though a pure steam atmosphere with Rrcb Rrcb the same water mass could not. ρi = ρi,rcb exp − 1 , (15) hi r where, assuming R = R so that atmospheric species 2.3 Compositional differentiation and homo rcb are diffusively separated, the scale height of the hydrogen hydrodynamic drag 2 /( ) is hi = kBTeqRrcb GMp µimp . The total mass loss rate is If an atmosphere is well-mixed at all altitudes, the relative then determined by combining the hydrogen flux with the abundances in the escaping gas are the same as those in flux of the heavier, dragged species determined by Equation bulk atmosphere. If this is not the case, then the prefer- (13). We typically find that the high mass loss rates follow- ential loss of one species over another leads to fractiona- ing an impact lead to fractionation factors xj ∼ 1, so that tion in the remaining atmosphere. This effect during the the relative flux rates of each species are proportional to hydrodynamic escape of hydrogen has been previously stud- their abundances in the total atmosphere. In other words, ied and proposed as a possible explanation for the relative the mass loss proceeds similarly to the case where the at- abundances of noble gases on the terrestrial planets (Hunten mosphere is assumed to be well-mixed to the outer radius 1973; Zahnle et al. 1990; Pepin 1991). The extent of this ef- (Rhomo = Rout). fect depends on the efficiency of vertical mixing and the rate This can be understood as follows. Approximating Fi = Û 2 Û 2 of atmospheric loss. Mi/(4πr µimp)' M/(4πr µmp) in a hydrogen-dominated at- An atmosphere is vertically well-mixed when its rate of mosphere, the fractionation factor can be rewritten as turbulent mixing is more rapid than molecular diffusion. The 4πµ µ − µ m2 b altitude where these processes become equally important is GMp i j i p i, j xj = 1 − (16) defined as the homopause. Above the homopause molecular MÛ kBTeq 1 + nj /ni diffusion dominates and atmospheric species may segregate 0.01 τ ⇒ x ∼ 1 − loss and (17) themselves according to their respective molecular weights. H2O f 7 Gyr Therefore, the density of a heavier atmospheric component 0.01 τ may be substantially depleted at high altitudes compared loss xCO2 ∼ 1 − , (18) to its bulk composition in the atmosphere. We assume the f 3 Gyr atmosphere is well-mixed to the RCB, Rhomo = Rrcb. This is 17 0.75 −1 −1 where we have used bH/He,H O ' 2.7 × 10 T cm s both because vigorous convection below the RCB promotes 2 and b ' 2.3 × 1017T0.75 cm−1 s−1 (Zahnle & Kast- turbulent mixing and because higher densities in this region H/He,CO2 ing 1986), and made the substitution MÛ ∼ M /τ ∼ slow the rate of molecular diffusion, which is ∝n−1 where n is env loss f M /τ . Following giant impacts, typical mass loss the atmospheric number density. In a static atmosphere, we p loss timescales for H/He envelopes with a mass fraction f ∼ 1 would therefore expect much lower concentrations of heavier per cent are 100 Myr (Biersteker & Schlichting 2019). species at altitudes much higher than the RCB. . Over these timescales, the fractionation factors for H O and In cases of rapid atmospheric escape, however, the 2 CO are ∼1, with H O marginally more easily entrained escaping flow of light molecules can drag along heavier 2 2 in the hydrogen outflow. Because the source of the out- molecules, leading to their escape. For a flux of a light species flow is the well-mixed convective region, this implies that F , the flux of a heavier species can be written (Zahnle et al. i F /n ' F /n , implying that the total escape 1990) H2O H2O H/He H/He flux is similar to what is expected from a well-mixed enve- Fj = Fi xj nj /ni, (13) lope (Rhomo = Rout). where nj and ni are the number densities of species j and i and xj ∈ [0, 1] is the fractionation factor, which determines the efficiency with which heavier molecules are dragged by 3 NUMERICAL RESULTS the flow. The fractionation factor is defined as To capture the full range of possible outcomes during the g µ − µ m b j i p i, j giant impact phase, we calculate the post-impact evolution xj = 1 − , (14) Fi kBT 1 + nj /ni of four example atmospheres, representing plausible atmo- where g is the gravitational acceleration and bi, j is the bi- spheric states during terrestrial planet formation, with the nary interaction cross-section. Radially varying quantities model described above. Specifically, we examine the loss of a are evaluated at the base of the outflow region, which we pure H/He atmosphere, a steam atmosphere, and two mixed take to be the RCB. H/He + H2O atmospheres. In each case, we assume that the We include the effects of compositional differentiation initial planet has Mp = 1M⊕, Tc,0 = 2000 K, and an orbital and hydrodynamic drag in our calculations and track the radius of a = 1.0 au around a sun-like star. We neglect any behavior of individual species as the atmosphere evolves. atmospheric mass contributed by the impactor and assume The mass loss rate is still governed by Equation (8). In the that the impact results in perfect accretion so that the new core mass is M 0 = M + M . Finally, we adopt η = 0.5 and energy-limited regime, the mass loss rate of each species is √ p p imp determined by simultaneously solving the above equation vimp = 2vesc (see Section4).
MNRAS 000,1–9 (2020) 6 J. B. Biersteker and H. E. Schlichting
increase the water vapor fraction of the atmosphere to the 1.0 1% H/He equivalent of 20 Earth oceans. Even then, a typical giant 1% H/He impact results in efficient loss of the atmosphere (red line +0.5% H2O in Figure3). These results can be understood by examining 0.8 1% H/He the mean molecular weights of the mixed atmospheres. The +0.1% H O 2 mean molecular weights of the atmospheres with ∼4 and 20 0.6 oceans’ worth of water vapor are µ ' 2.5 and 3.3, respec-
X tively. In both cases, this is significantly less than the max- imum mean molecular weight atmosphere that can be lost 0.4 Impact shocks given by Equation (12). As expected, the higher mean molec- ular weight of the mixed atmospheres does require slightly 0.2 larger impactors to achieve the same loss fraction as pure H/He atmospheres (see Equation (11)). Additionally, we find 0.1% H O (steam) 0.0 2 that simulations using a well-mixed atmosphere model re- quire slightly larger impactors to achieve atmospheric loss 0.01 0.1 1.0 comparable to that in simulations using the fractionation M /M imp p model (dashed and solid green lines, respectively, in Figure 3). This is because the atmospheric hydrogen is more ex- tended in the case where the radiative region is un-mixed Figure 3. Mass fraction of envelope lost (X) as a function of (hH/He > h), lowering the impactor mass needed to initiate impactor mass (Mimp/Mp ) for planets with Mp = 1M⊕, Rp = 1R⊕, atmospheric loss. Because the loss of heavier species is driven T = 2000 K, a = 1 au, and different envelope compositions and p by the hydrogen flux, heavier molecules are also lost more mass fractions (as indicated by their respective labels). The two easily in the un-mixed case, despite their much smaller scale green and red lines for H/He+H2O envelope compositions represent the two different mixing models: the bold solid line corresponds heights. to hydrodynamic escape with mass fractionation, while the thin Finally, comparison with the estimated shock-induced dashed line is for a well-mixed envelope. The black dashed line losses from Schlichting et al.(2015) indicates that atmo- indicates the atmospheric envelope loss from the impact shock for spheric loss from thermal expansion dominates the primor- vimp ' vesc (Schlichting et al. 2015, Equation (33)). dial and mixed envelopes for most of the impactor mass range considered, with the two processes becoming compa- rable only when Mimp ∼ Mp. For the secondary atmosphere, For the primordial atmosphere case, we assume an ini- however, shocks are the more important process at all im- tial H/He envelope with a mass fraction f = Menv/Mp = 1 pactor masses. per cent. This approximately corresponds to the envelope mass that an Earth-mass planet could accrete if it grew to 50 per cent of its current mass prior to the dissipation of & 4 DISCUSSION AND CONCLUSIONS the gas disk (Lammer et al. 2014; Lee & Chiang 2015; Inam- dar & Schlichting 2015). We find, consistent with previous The above case studies allow us to draw several broad con- work (Biersteker & Schlichting 2019), that the H/He enve- clusions about the potential role of giant impacts in the at- lope is readily lost during giant impacts for the impactor mospheric evolution of terrestrial planets. Planets with light mass fractions Mimp/Mp & 0.05 (Figure3), implying that hydrogen-dominated envelopes can lose the bulk of their at- a Mars-sized impactor hitting the proto-Earth would have mospheres through a planetary wind following a large im- ejected the entire H/He envelope. pact (Mimp/Mp & 0.05). Thermal loss following an impact The next case we consider is a secondary atmosphere can also deplete the inventory of high mean molecular weight entirely consisting of outgassed H2O (orange line in Figure volatiles, provided they are sufficiently diluted by H/He in 3). Here we assume an atmospheric mass of f = 0.1 per cent, the bulk atmosphere. These results are remarkable in that which is equivalent in mass to about four Earth oceans, a they provide an efficient way to lose several oceans’ worth moderate estimate for the outgassed water content of the of water and other volatile species that are predicted to be magma ocean (Elkins-Tanton 2008; Hamano et al. 2013). present in the atmospheres of terrestrial planets as a result In contrast to a primordial H/He atmosphere, the pure sec- of outgassing. Thermal atmospheric loss triggered by giant ondary steam atmosphere undergoes no appreciable atmo- impacts can therefore radically alter the volatile budget of spheric loss due to the heating and thermal expansion of terrestrial planets. Total atmospheric loss in these scenarios the envelope following the impact in any of the scenarios would likely be enhanced by impact shocks and high-energy considered here. stellar irradiation, so that the cumulative loss after an im- Most notable are the atmospheric mass-loss results for pact may be even greater (e.g., Lopez et al. 2012; Inamdar secondary atmospheres that are outgassed into a hydrogen- & Schlichting 2015; Liu et al. 2015; Erkaev et al. 2016). dominated envelope (green and red lines in Figure3). To Atmospheric loss following a giant impact is highly de- model this scenario, we consider two atmospheric composi- pendent on the atmosphere’s composition. For impactors tions. First, we combine the primary and secondary atmo- with vimp ∼ vesc, we find that thermal escape of a purely spheres described above for a total mass fraction f = 1.1 outgassed secondary atmosphere following an impact is neg- per cent. For such a mixed atmosphere, we find complete ligible, and that the total impact-driven loss is dominated by atmospheric loss for impactor masses exceeding about 5 per shocks. We do not model impact shocks in this work, but pre- cent of the target mass (green line in Figure3). Second, we vious studies suggest that, for impactors with vimp ∼ vesc, re-
MNRAS 000,1–9 (2020) Effects of composition on impact-driven atmospheric loss 7
Initial planet (A) Mixed atmosphere lost in impact
Early impact Final states Late impact
Primary atmosphere lost (B) Outgassing
Figure 4. Example of divergent atmospheric outcomes following a collision at different times during the giant impact phase. Planets A and B (top and bottom rows, respectively) begin with identical primary atmospheres (shown in blue) and outgas identical secondary atmospheres (depicted in green). However, planet A loses its mixed primordial-secondary atmosphere during an impact and is left with a negligible atmosphere after the impact and with a significantly depleted volatile budget. Planet B suffers an impact after the primary atmosphere is lost by either photoevaporation or an earlier impact such that it has a pure secondary atmosphere at the time of the giant impact depicted in the figure. In this case, it retains most of its secondary atmosphere and hence most of its volatiles during subsequent giant impacts.
moval of the entire envelope requires Mimp/Mp ∼ 1 (Schlicht- planetary embryos, Walsh & Levison(2019) find a delayed ing et al. 2015; Inamdar & Schlichting 2015). In the canonical onset for the dynamical instability with impacts beginning moon-forming impact (Mimp/Mp ∼ 0.1), atmospheric losses at ∼10–20 Myr. For an estimated value of η ' 0.5 (Carter of ∼20 per cent are expected due to shock, though this can be et al. 2020) most of these impacts can be expected to en- enhanced for some choices of the proto-Earth’s initial con- tirely remove H/He envelopes. The only exceptions are the ditions (Genda & Abe 2003, 2005) and scenarios with larger impactors at the low end of the mass range in the “Grand impactor masses (Mimp/Mp ∼ 1) have been proposed (Canup Tack” scenario. Therefore, while similar planets may expe- 2012). Therefore, while some scenarios may produce signif- rience a diversity of impact scenarios, it appears that the icant erosion of secondary atmospheres, these atmospheres most plausible range of impactor masses and velocities will are generally difficult to dislodge in a single giant impact efficiently remove primordial H/He envelopes as well as en- event. velopes that consist of a mix of primordial hydrogen and secondary outgassed species. In addition to atmospheric composition, the outcome of a giant impact depends on the details of the impact The variable timing of impacts, however, may result scenario, in particular the impactor size (Mimp/Mp), im- in significantly different evolutionary pathways for planets. pact velocity (vimp), and timing. Terrestrial planet forma- This is because the giant impact phase occurs concurrently tion likely encompasses a broad array of plausible scenar- with multiple other evolutionary processes for terrestrial ios and the details for any particular planet are stochas- planet atmospheres. Specifically, a proto-planet with a large tic, potentially producing dramatically different results. In fraction of its volatiles mixed into a hydrogen-dominated simulations adopting either the current orbits of the giant atmosphere may lose its entire atmosphere in a giant im- planets or a more compact and circular configuration, im- pact, resulting in a substantially depleted volatile inventory. pactor mass ratios range from Mimp/Mp ∼ 0.1–0.3, impact On the other hand, an identical proto-planet that lost its velocities are ∼1–1.3vesc, and impact times are ∼10–100 Myr primordial atmosphere either by an early giant impact or after the formation of the solar system (e.g., O’Brien et al. through high-energy stellar radiation (Johnstone et al. 2015) 2006). In the “Grand Tack” scenario, characterized by in- prior to outgassing of a secondary atmosphere will be able ward then outward migration of Jupiter in the presence of to retain most of its secondary atmosphere and, hence, its the gas disk (Walsh et al. 2011), a wider range of impact volatile budget, for the duration of the giant impact phase scenarios occur: Mimp/Mp ∼ 0.01–0.3, vimp ∼ 1–2vesc, and (see Figure4). Early removal of a hydrogen envelope by an impact times of ∼0.1–100 Myr (O’Brien et al. 2014). In a impact may allow a planet’s magma ocean to cool and out- model including both the dissipating disk and the growth of gas earlier than it otherwise would. The components of this
MNRAS 000,1–9 (2020) 8 J. B. Biersteker and H. E. Schlichting outgassed atmosphere might then face greater erosion from Barboni M., Boehnke P., Keller B., Kohl I. E., Schoene B., Young increased cumulative exposure to planetesimal impacts and E. D., McKeegan K. D., 2017, Science Advances, 3, e1602365 high-energy stellar radiation (Schlichting et al. 2015). Biersteker J. B., Schlichting H. E., 2019, MNRAS, 485, 4454 In the specific case of the solar system, evidence for the Canup R. M., 2012, Science, 338, 1052 accretion of Mars on a timescale comparable to the gas disk Canup R. M., Asphaug E., 2001, Nature, 412, 708 Carter P. J., Lock S. J., Stewart S. T., 2020, Journal of Geophys- lifetime (Dauphas & Pourmand 2011) suggests the plane- ical Research (Planets), 125, e06042 tary embryos that formed the inner planets grew quickly Chambers J. E., Wetherill G. W., 1998, Icarus, 136, 304 enough to accrete a primordial atmosphere. The composi- Dauphas N., Pourmand A., 2011, Nature, 473, 489 tion of noble gases in the Earth’s mantle is consistent with Elkins-Tanton L. T., 2008, Earth and Planetary Science Letters, the dissolution of gas from such an atmosphere into the early 271, 181 Earth’s magma ocean (Harper & Jacobsen 1996; Mukhopad- Elkins-Tanton L. T., 2011, Ap&SS, 332, 359 hyay 2012; Williams & Mukhopadhyay 2019). If this is the Elkins-Tanton L. T., 2012, Annual Review of Earth and Planetary case, the atmospheric evolution of these planets may have Sciences, 40, 113 been significantly affected by giant impacts. This scenario is Erkaev N. V., Lammer H., Odert P., Kislyakova K. G., Johnstone consistent with xenon abundance measurements that sug- C. P., Gudel¨ M., Khodachenko M. L., 2016, MNRAS, 460, 1300 gest the proto-Earth once hosted roughly four times its Freedman R. S., Marley M. S., Lodders K., 2008, ApJS, 174, 504 presently observed volatile inventory (Avice & Marty 2014) Fressin F., et al., 2013, ApJ, 766, 81 and therefore could offer an explanation for the depletion Genda H., Abe Y., 2003, Icarus, 164, 149 of Earth’s and Mars’s light noble gas budgets compared to Genda H., Abe Y., 2005, Nature, 433, 842 Venus (Pepin 1997; Halliday 2013). Ginzburg S., Schlichting H. E., Sari R., 2016, ApJ, 825, 29 Giant impacts are an expected component of terres- Guillot T., Chabrier G., Gautier D., Morel P., 1995, ApJ, 450, trial planet formation, and have the potential to signifi- 463 cantly shape the volatile inventories of accreting terrestrial Halliday A. N., 2013, Geochim. Cosmochim. Acta, 105, 146 planets. We find that the outcome of an impact depends Hamano K., Abe Y., Genda H., 2013, Nature, 497, 607 strongly on the atmospheric composition. Low mean molec- Harper Charles L. J., Jacobsen S. B., 1996, Science, 273, 1814 ular weight envelopes can be efficiently removed by a plane- Hayashi C., Nakazawa K., Mizuno H., 1979, Earth and Planetary Science Letters, 43, 22 tary wind following these events, and loss by this mechanism Hunten D. M., 1973, Journal of Atmospheric Sciences, 30, 1481 significantly exceeds that caused by the initial impact shock. Hunter J. D., 2007, Computing In Science & Engineering, 9, 90 Additionally, the rate of hydrogen escape following an im- Inamdar N. K., Schlichting H. E., 2015, MNRAS, 448, 1751 pact is so rapid that large reservoirs of high mean molecular Izidoro A., Raymond S. N., Morbidelli A. r., Winter O. C., 2015, weight volatiles can be lost if they are sufficiently diluted in MNRAS, 453, 3619 a hydrogen-dominated atmosphere. In contrast, purely sec- Johansen A., Lambrechts M., 2017, Annual Review of Earth and ondary atmospheres (e.g., steam envelopes) are only mod- Planetary Sciences, 45, 359 erately diminished during giant impacts, unless Mimp ∼ Mp. Johnstone C. P., et al., 2015, ApJ, 815, L12 For these atmospheres, loss is due to the impact shock and Kite E. S., Fegley Bruce J., Schaefer L., Ford E. B., 2020, ApJ, thermal loss following the impact is negligible. Because gi- 891, 111 Lammer H., et al., 2014, MNRAS, 439, 3225 ant impacts occur stochastically during a time when nascent Lammer H., et al., 2018, A&A Rev., 26, 2 planetary atmospheres are being altered by other processes, Lee E. J., Chiang E., 2015, ApJ, 811, 41 including outgassing and photoevaporation, these impacts Liu S.-F., Hori Y., Lin D. N. C., Asphaug E., 2015, ApJ, 812, 164 likely play a key role in determining the ultimate volatile Lopez E. D., Fortney J. J., Miller N., 2012, ApJ, 761, 59 inventories of terrestrial planets, leading to a large diversity Morbidelli A., Lunine J. I., O’Brien D. P., Raymond S. N., Walsh in composition among otherwise similar worlds. K. J., 2012, Annual Review of Earth and Planetary Sciences, 40, 251 Mukhopadhyay S., 2012, Nature, 486, 101 O’Brien D. P., Morbidelli A., Levison H. F., 2006, Icarus, 184, 39 ACKNOWLEDGEMENTS O’Brien D. P., Walsh K. J., Morbidelli A., Raymond S. N., Man- This research has been supported by the National Aeronau- dell A. M., 2014, Icarus, 239, 74 tics and Space Administration under grant No. 17 XRP17 2− Odert P., et al., 2018, Icarus, 307, 327 0055 issued through the Exoplanet Research Program. This Pepin R. O., 1991, Icarus, 92, 2 Pepin R. O., 1997, Icarus, 126, 148 research made use of the software packages NumPy (Van Piso A.-M. A., Youdin A. N., 2014, ApJ, 786, 21 Der Walt et al. 2011), SciPy (Virtanen et al. 2020), Astropy Rafikov R. R., 2006, ApJ, 648, 666 (Astropy Collaboration et al. 2013), and matplotlib (Hunter Schlichting H. E., Sari R., Yalinewich A., 2015, Icarus, 247, 81 2007). Valencia D., O’Connell R. J., Sasselov D., 2006, Icarus, 181, 545 Valencia D., Guillot T., Parmentier V., Freedman R. S., 2013, ApJ, 775, 10 REFERENCES Van Der Walt S., Colbert S. C., Varoquaux G., 2011, Computing in Science & Engineering, 13, 22 Alf`eD., Price G. D., Gillan M. J., 2001, Phys. Rev. B, 64, 045123 Virtanen P., et al., 2020, Nature Methods, 17, 261 Asphaug E., 2014, Annual Review of Earth and Planetary Sci- Walsh K. J., Levison H. F., 2019, Icarus, 329, 88 ences, 42, 551 Walsh K. J., Morbidelli A., Raymond S. N., O’Brien D. P., Man- Astropy Collaboration et al., 2013, A&A, 558, A33 dell A. M., 2011, Nature, 475, 206 Avice G., Marty B., 2014, Philosophical Transactions of the Royal Williams C. D., Mukhopadhyay S., 2019, Nature, 565, 78 Society of London Series A, 372, 20130260 Zahnle K. J., Kasting J. F., 1986, Icarus, 68, 462
MNRAS 000,1–9 (2020) Effects of composition on impact-driven atmospheric loss 9
Zahnle K., Kasting J. F., Pollack J. B., 1990, Icarus, 84, 502
This paper has been typeset from a TEX/LATEX file prepared by the author.
MNRAS 000,1–9 (2020)