# Escape of the martian protoatmosphere and initial water inventory

N. V. Erkaev, H. Lammer, L. Elkins-Tanton, A. Stökl, P. Odert, Emmanuel

Marcq, E. A. Dorfi, K. G. Kislyakova, Yu. N. Kulikov, M. Leitzinger, et al.

► To cite this version:

N. V. Erkaev, H. Lammer, L. Elkins-Tanton, A. Stökl, P. Odert, et al. Escape of the martian protoatmosphere and initial water inventory. Planetary and Space Science, 2014, 98, pp.106-119. 10.1016/j.pss.2013.09.008. hal-00864122

HAL Id: hal-00864122

https://hal.science/hal-00864122

Submitted on 27 Sep 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Escape of the martian protoatmosphere and initial water inventory

<figure><img src="image_2.png" /></figure>

CrossMark

N.V. Erkaev<sup>a,b</sup>, H. Lammer<sup>c,*</sup>, L.T. Elkins-Tanton<sup>d</sup>, A. Stökl<sup>e</sup>, P. Odert<sup>c,h</sup>, E. Marcq<sup>f</sup>, E.
A. Dorfi<sup>e</sup>, K.G. Kislyakova<sup>c</sup>, Yu.N. Kulikov<sup>g</sup>, M. Leitzinger<sup>h</sup>, M. Güdel<sup>e</sup>

## ARTICLE INFO

**Article history:**

Received 28 December 2012

Received in revised form

20 July 2013

Accepted 6 September 2013

Available online 17 September 2013

**Keywords:**

Early mars

Protoatmospheres

Atmospheric escape

Evolution

## ABSTRACT

Latest research in planet formation indicates that Mars formed within a few million years (Myr) and remained as a planetary embryo that never grew to a more massive planet. It can also be expected from dynamical models that most of Mars' building blocks consisted of material that formed in orbital locations just beyond the ice line which could have contained ~0.1-0.2 wt.% of H₂O. By using these constraints, we estimate the nebula-captured and catastrophically outgassed volatile contents during the solidification of Mars' magma ocean and apply a hydrodynamic upper atmosphere model for the study of the soft X-ray and extreme ultraviolet (XUV) driven thermal escape of the martian protoatmosphere during the early active epoch of the young Sun. The amount of gas that has been captured from the protoplanetary disk into the planetary atmosphere is calculated by solving the hydrostatic structure equations in the protoplanetary nebula. Depending on nebular properties such as the dust grain depletion factor, planetesimal accretion rates and luminosities, hydrogen envelopes with masses $≥ 3 × 10^{19}$ g to $≤ 6.5 × 10^{22}$ g could have been captured from the nebula around early Mars. Depending on the before mentioned parameters, due to the planets low gravity and a solar XUV flux that was ~100 times stronger compared to the present value, our results indicate that early Mars would have lost its nebular captured hydrogen envelope after the nebula gas evaporated, during a fast period of ~0.1-7.5 Myr. After the solidification of early Mars' magma ocean, catastrophically outgassed volatiles with the amount of ~50-250 bar H₂O and ~10-55 bar CO₂ could have been lost during ~0.4-12 Myr, if the impact related energy flux of large planetesimals and small embryos to the planet's surface lasted long enough, that the steam atmosphere could have been prevented from condensing. If this was not the case, then our results suggest that the timescales for H₂O condensation and ocean formation may have been shorter compared to the atmosphere evaporation timescale, so that one can speculate that sporadically periods, where some amount of liquid water may have been present on the planet's surface. However, depending on the amount of the outgassed volatiles, because of impacts and the high XUV-driven atmospheric escape rates, such sporadically wet surface conditions may have also not lasted much longer than ~0.4-12 Myr. After the loss of the captured hydrogen envelope and outgassed volatiles during the first 100 Myr period of the young Sun, a warmer and probably wetter period may have evolved by a combination of volcanic outgassing and impact delivered volatiles ~4.0 ± 0.2 Gyr ago, when the solar XUV flux decreased to values that have been <10 times that of today's Sun.

© 2013 Elsevier Ltd. Open access under CC BY license.

## 1. Introduction

The formation of Mars' nebula-captured, catastrophically degassed and impact delivered protoatmosphere is directly connected to the planet's formation time scale, the nebula dissipation time, its orbital

location and the planet's small mass compared to Earth and Venus. Chassefière (1996a,b) investigated for the first time the hydrodynamic loss of oxygen from primitive atmospheres of Venus and Mars in detail. However, the pioneering studies of Chassefière (1996a,b) are based on meanwhile outdated terrestrial planet formation models in which the time of the final accretion for terrestrial planets occurred $≥ 100$ Myr after the formation of the Sun (Wetherill, 1986). Furthermore, in these pioneering studies by Chassefière (1996a,b) the cooling phase of the magma ocean was expected to occur after ~100 Myr,

<footnote><sup>a</sup> Institute for Computational Modelling, 660041 Krasnoyarsk 36, Russian Academy of Sciences, Russian Federation</footnote>

<footnote><sup>b</sup> Siberian Federal University, 660041 Krasnoyarsk, Russian Federation</footnote>

<footnote><sup>c</sup> Space Research Institute, Austrian Academy of Sciences, Schmiedlstrasse 6, A-8042 Graz, Austria</footnote>

<footnote><sup>d</sup> Department of Terrestrial Magnetism, Carnegie Institution for Science, Washington DC 20015, USA</footnote>

<footnote><sup>e</sup> Institute for Astronomy, University of Vienna, Türkenschanzstraße 17, 1180 Vienna, Austria</footnote>

<footnote><sup>f</sup> LATMOS, Université de Versailles Saint-Quentin-en-Yvelines, Guyancourt, France</footnote>

<footnote><sup>g</sup> Polar Geophysical Institute, Russian Academy of Sciences, Khalturina 15, 183010 Murmansk, Russian Federation</footnote>

<footnote><sup>h</sup> Institute of Physics, IGAM, University of Graz, Universitätsplatz 5, A-8010 Graz, Austria</footnote>

<footnote>* Corresponding author. Tel.: +43 31 6412 0641; fax: +43 31 6412 0690. E-mail address: helmut.lammer@oeaw.ac.at (H. Lammer).</footnote>
while more recent studies indicate that the solidification of magma oceans even with depths of up to ~2000 km is a fast process and mantle solidification of ~98% can be completed in ≤ 5 Myr (e.g. Elkins-Tanton, 2008, 2011; Marcq, 2012; Lebrun et al., 2013; Hamano et al., 2013). Moreover, it is also important to note that the assumption of several previous studies, that terrestrial planets, including early Mars finished their accretion late, resulted also in ages where the soft X-ray and extreme ultraviolet (XUV) flux of the young Sun was much lower compared to the high XUV flux values, which are now known from multi-wavelength observations of so-called young solar proxies (e.g., Güdel et al., 1997; Ribas et al., 2005; Güdel, 2007; Claire et al., 2012). Because of the lack of accurate data, Chassefière (1996a,b) applied an XUV enhancement factor as its highest value which was ~25 times higher than that of the present Sun.

In a recent review article on Mars' origin Brasser (2013) argued that Mars' small mass requires the terrestrial planets to have formed from a narrow annulus of material, rather than a disc extending to Jupiter. The truncation of the outer part of the disc was most likely related to migration of the gas giants, which kept the mass of Mars small. For the formation of the martian proto-atmosphere this evidence from planet formation and latest dynamical models has important implications, because it would mean that Mars formed within a few million years and can be considered as a planetary embryo that never grew to a "real" more massive planet. Moreover, from the latest martian formation modeling scenarios most likely related to migration of the giants (Walsh et al., 2011), it is expected that most of the planet's building blocks consists of material that formed in a region just behind the ice line, so that the materials were more water-rich than the materials that were involved in the accretion of Venus and Earth.

Brasser (2013) suggests that the building blocks of early Mars could have consisted of ~0.1–0.2 wt.% of H₂O. The results presented in Brasser (2013) which are based on studies by Walsh et al. (2011) agree in the amount of Mars' initial water inventory with Lunine et al. (2003) who applied also a dynamical model which yielded longer formation time scales. However, it should also be pointed out that model studies which consider different impact regimes than the before mentioned studies can also result in an early Mars which originated drier (Horner et al. 2009). Although, it is obvious that our current knowledge of terrestrial planet formation and its related hydration is presently insufficient there is geomorphological evidence for water on early Mars, where ~90% was most likely outgassed and/or delivered during the first Gyr (e.g., Chassefière, 1996b; Chassefière et al., 2013; Baker, 2001; Lammer et al., 2013a).

The main aim of the present study is to investigate in detail how long the before mentioned nebular captured and catastrophically outgassed protoatmospheres have been stable after Mars' origin, and to understand how long the early planet's protoatmosphere survived against thermal atmospheric escape. In Section 2 the formation of a nebula captured hydrogen envelope on early Mars and the expected catastrophically outgassed steam-type protoatmosphere based on materials which contain ~0.1–0.2 wt.% H₂O (Brasser, 2013) is described. In Section 3 we discuss the early XUV radiation environment of the young Sun and the life time of the nebula gas which determines the age when the planet's protoatmosphere was exposed freely to the high solar XUV radiation field. In Section 4 we study the upper atmosphere structure and the escape of the martian protoatmosphere by applying a time-dependent numerical algorithm, which is able to solve the system of 1-D fluid equations for mass, momentum, and energy conservation. Finally we describe the solar and atmospheric input parameters of the applied model and discuss the results.

## 2. Nebula-based and catastrophically outgassed protoatmospheres

For studying the potential habitability and atmosphere evolution of Mars, it is important to understand which sources and sinks contributed to the formation of the planet's initial atmosphere and water inventory. Furthermore, a detailed investigation on the escape-related evolution of the early martian protoatmosphere is important for understanding how long Mars may have had surface conditions that standing bodies of liquid water could have existed on the planet's surface. Generally four main processes are responsible for the formation of planetary atmospheres

* capture of hydrogen and other gases (He, noble gases, etc.) from the solar nebular,

* catastrophic outgassing of volatiles such as H₂O, CO₂, etc. and the formation of a steam atmosphere during and after the magma ocean solidification period,

* impact delivery of volatiles by asteroids and comets,

* degassing by volcanic processes during geological epochs.

Fig. 1 illustrates the expected atmosphere formation and loss scenarios for Mars during the planet's history. In the present work we focus on the origin and the evolution of the earliest martian protoatmosphere, consisting of hydrogen accumulated from the solar nebular and a catastrophically outgassed steam atmosphere after the planet finished its accretion and the magma ocean solidified.

### 2.1. Captured hydrogen envelope around early Mars

When proto-planets grow within the surrounding solar nebula by accretion of planetesimals, an extensive amount of gas will be attracted so that optically thick, dense hydrogen envelopes accumulate around a rocky core (e.g., Mizuno et al., 1978; Hayashi et al., 1979; Wuchterl, 1993; Ikoma et al., 2000; Ikoma and Genda, 2006; Rafikov, 2006). The structure of such nebular-based hydrogen atmospheres was investigated decades ago by Hayashi et al. (1979) and Nakazawa et al. (1985) for a wide range of planetary accretion rates, grain opacities, and gas disk densities. These pioneering studies obtained captured nebula gas around a Mars-mass body (i.e. ~0.1M⊕) of 8.4 × 10²⁴ g during nebular life times of ~1–10 Myr, equivalent to the hydrogen content of ~55 Earth oceans (1EOH ≈ 1.53 × 10²³ g). More recent studies improved on these earlier results by adoption of realistic gas and dust opacities as well as a realistic equation of state leading to significantly lower atmosphere masses around bodies with masses that are ~0.1M⊕ (Ikoma and Genda, 2006).

For the present investigation we computed a set of atmospheric models for Mars to obtain an estimate of the amount of gas collected from the protoplanetary disk into the planetary atmosphere. The hydrostatic structure equations have been solved by using the initial model integrator of the adaptive, implicit RHD-Code (TAPIR-code) the equation of state from Saumon et al. (1995), gas opacities from Freedman et al. (2008), and dust opacities by Semenov et al. (2003). Convective energy transport is included in TAPIR in the form of a turbulent convection model loosely based on the description by Kuhfuß (1987).

For the conditions of the solar nebula at the position of the Mars orbit we assumed a gas density of 5 × 10⁻¹⁰ g cm⁻³ and a temperature of 200 K. These values are in good agreement with restraints derived from the minimum-mass solar nebula (Hayashi, 1981). The minimum-mass solar nebula (MSN) is a protoplanetary disk that contains the minimum amount of solid material which is necessary to build the planets of the Solar system.
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Illustration of Mars' origin and protoatmosphere formation and evolution. The dotted lines correspond to the accumulation during the growth and escape of nebula-based hydrogen from proto-Mars. The onset of escape corresponds to the nebula dissipation time around ~3–10 Myr, which is also the expected time period when Mars finished its accretion (Brasser, 2013). The short dashed lines illustrate the catastrophically outgassed volatiles and their expected escape after the planet's magma ocean solidified. Later on when the solar activity decreased a secondary CO₂ atmosphere could have build up by volcanic activity (Grott et al., 2011; Lammer et al., 2013a) and the late heavy bombardment may also have delivered volatiles to Mars ~3.8 Gyr ago.</figcaption></figure>

The outer boundary conditions, i.e. nebula density and temperature, have been implemented at the Hill radius $r_{\text{Hill}}$ for all models as we consider $r_{\text{Hill}}$ to be a good approximation for the place where the essentially hydrostatic structure of the planetary atmosphere blends into the background disk structure. However, when calculating the captured atmospheric masses, i.e. the amount of gas in effect gravitationally bound to the planet, we used the minimum of $r_{\text{Hill}}$ and the Bondi radius $r_{\text{Bondi}}$, which turns out to be equal to the latter for all model runs by a margin of about a magnitude. The definition of the outer boundary condition seems to be, apart from the equation of state and nebular opacities, the main cause for the different captured atmospheric masses obtained by different authors. According to Ikoma (2012; private communication), the discrepancy between Ikoma and Genda (2006) and Hayashi et al. (1979) is a case in point. In general, as also described by Ikoma and Genda (2006), the atmospheres (and thus the atmospheric masses) of low-mass planets such as Mars are more dependent on outer boundary conditions than atmospheres of more massive Earth-like and super-Earth-type cores.

In order to get some measure of the uncertainties involved in our modeling, we covered a small parameter space by varying the most important atmospheric parameters: the planetary luminosity $L_{\text{pl}}$ and the dust grain depletion factor $f$. Table 1 and Fig. 2 summarize the results of our atmospheric calculations.

$L_{\text{pl}}$ is related to the rate of infalling planetesimals

$$L_{\text{pl}} \approx GM_{\text{pl}}\dot{M}_{\text{acc}}\left(\frac{1}{r_{\text{pl}}} - \frac{1}{r_{\text{Hill}}}\right), \quad (1)$$

with G the Newton gravitational constant, planetary mass $M_{\text{pl}}$, planetary radius $r_{\text{pl}}$ and planetesimal accretion rate $\dot{M}_{\text{acc}}$. Taking into account that according to Walsh et al. (2011) and Brasser (2013) Mars' formation was completed before or soon after the nebular gas disappeared at ~3–10 Myr, $M_{\text{pl}}/\dot{M}_{\text{acc}}$ should be several ~$10^6$ years or larger. On the other hand, according to Elkins-Tanton (2008) and Hamano et al. (2013) the cooling time scale of a Mars-size planet could be well above 1 Myr and thus it

<caption>**Table 1**</caption>

Integral parameters for Mars model atmospheres models with dust depletion factors $f$ of 0.1, 0.01, and 0.001 and for accretion rates $\dot{M}_{\text{acc}}$ between $1 \times 10^{-6}$ and $1 \times 10^{-9}$ Earth masses per year. $L$ is the luminosity resulting from the accretion of planetesimals; $M_{\text{atm}}$ is the atmospheric mass up to the Bondi radius; and surface pressure and temperature on the surface are denoted as $P_s$ and $T_s$, respectively.

<table><thead><tr><th>M<sub>acc</sub> (M<sub>Mars</sub>/yr)</th><th>f<sub>dust</sub></th><th>L (erg/s)</th><th>M<sub>atm</sub> (g)</th><th>P<sub>s</sub> (bar)</th><th>T<sub>s</sub> (K)</th></tr></thead><tbody><tr><td>9.35 &times; 10<sup>-9</sup></td><td>0.001</td><td>2.39 &times; 10<sup>22</sup></td><td>6.58 &times; 10<sup>22</sup></td><td>7.81</td><td>600</td></tr><tr><td>9.35 &times; 10<sup>-9</sup></td><td>0.01</td><td>2.39 &times; 10<sup>22</sup></td><td>3.21 &times; 10<sup>22</sup></td><td>3.38</td><td>639</td></tr><tr><td>9.35 &times; 10<sup>-9</sup></td><td>0.1</td><td>2.39 &times; 10<sup>22</sup></td><td>1.00 &times; 10<sup>22</sup></td><td>8.08</td><td>690</td></tr><tr><td>9.35 &times; 10<sup>-8</sup></td><td>0.001</td><td>2.38 &times; 10<sup>23</sup></td><td>2.66 &times; 10<sup>22</sup></td><td>2.58</td><td>693</td></tr><tr><td>9.35 &times; 10<sup>-8</sup></td><td>0.01</td><td>2.38 &times; 10<sup>23</sup></td><td>9.76 &times; 10<sup>21</sup></td><td>0.763</td><td>724</td></tr><tr><td>9.35 &times; 10<sup>-8</sup></td><td>0.1</td><td>2.38 &times; 10<sup>23</sup></td><td>2.86 &times; 10<sup>21</sup></td><td>0.155</td><td>754</td></tr><tr><td>9.35 &times; 10<sup>-7</sup></td><td>0.001</td><td>2.38 &times; 10<sup>24</sup></td><td>8.81 &times; 10<sup>21</sup></td><td>0.628</td><td>795</td></tr><tr><td>9.35 &times; 10<sup>-7</sup></td><td>0.01</td><td>2.38 &times; 10<sup>24</sup></td><td>2.84 &times; 10<sup>21</sup></td><td>0.151</td><td>784</td></tr><tr><td>9.35 &times; 10<sup>-7</sup></td><td>0.1</td><td>2.38 &times; 10<sup>24</sup></td><td>5.25 &times; 10<sup>20</sup></td><td>0.028</td><td>841</td></tr><tr><td>9.35 &times; 10<sup>-6</sup></td><td>0.001</td><td>2.38 &times; 10<sup>25</sup></td><td>2.70 &times; 10<sup>21</sup></td><td>0.132</td><td>885</td></tr><tr><td>9.35 &times; 10<sup>-6</sup></td><td>0.01</td><td>2.38 &times; 10<sup>25</sup></td><td>5.14 &times; 10<sup>20</sup></td><td>0.028</td><td>862</td></tr><tr><td>9.35 &times; 10<sup>-6</sup></td><td>0.1</td><td>2.38 &times; 10<sup>25</sup></td><td>3.21 &times; 10<sup>19</sup></td><td>0.005</td><td>960</td></tr></tbody></table>

seems plausible that during the nebula-gas accumulation phase the heat flux from the interior significantly adds to the planetary luminosity. The lower limit of $L_{\text{pl}}$ can be constrained by the radiogenic luminosity estimated to be ~$10^{20}$ erg s⁻¹ for Mars (Wänke and Dreibus, 1988).

For higher planetary luminosities $T_s$ almost reaches 1000 K and it is well likely that models with, e.g. other boundary conditions or different dust opacity data yield even higher surface temperatures. It is important to note that H₂O can also be produced on a planet if $T_s > 1500$ K. In such a case the planet's surface melts and atmospheric hydrogen can be oxidized by oxides such as wüstite, magnetite and fayalite, which are inside the planet to produce H₂O on the planet (Sasaki, 1990; Ikoma and Genda, 2006). However, the model results which yield high surface temperatures
<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Nebular-captured hydrogen envelopes for a Mars size and mass object at 1.5 AU, in units of Earth ocean equivalent amounts of hydrogen ($1\text{EO}_\text{H} = 1.53 \times 10^{23} \text{ g}$) as a function of luminosity for three different dust grain depletion factors $f=0.001$ (dashed-dotted line), $f=0.01$ (dashed line), and $f=0.1$ (dotted line).</figcaption></figure>

are also those with only comparatively thin hydrogen envelopes, which is reasonable as high luminosities and temperatures tend to inflate a planetary atmosphere. Therefore, one may speculate that Mars atmospheres with $T_s > 1500$ K will be too thin to allow for efficient H₂O production from a captured and oxidized hydrogen envelope.

Before we discuss the radiation environment of the young Sun during the first 100 Myr after Mars' origin and before we model the escape of the nebula-based hydrogen envelope we investigate the possible range of catastrophically outgassed steam atmospheres.

## 2.2. Magma ocean and outgassing of a steam atmosphere on early Mars

As discussed before the terrestrial planets are thought to have reached their final sizes by a series of giant accretionary impacts. These impacts were energetic enough to produce melting of some depth in the planet (e.g., Tonks and Melosh, 1993; Reese and Solomatov, 2006; Lebrun et al., 2013). This hypothesis is supported by the discovery of ¹⁴²Nd isotope anomalies in martian SNC meteoroids, which indicate that early Mars developed a magma ocean (Harper et al., 1995; Foley et al., 2005; Debaille et al., 2007). Therefore, the first major degassed volatile-rich atmospheres likely resulted from the solidification of these magma bodies, and their release into the growing atmosphere in excess of what can be held in crystallizing silicate minerals (Abe, 1993, 1997; Abe and Matsui, 1988; Matsui and Abe, 1986; Zahnle et al., 1988; Elkins-Tanton et al., 2005; Debaille et al., 2007; Elkins-Tanton, 2008, 2011; Hamano et al., 2013; Lebrun et al., 2013). In these models the magma ocean is expected to solidify from the bottom upward, because the slope of the adiabat is steeper than the slope of the solidus and thus they first intersect at depth. Because the energy and size of late accretionary impacts on early Mars are unknown, we consider a 500 km-deep magma ocean and, as an end-member, a 2000 km-deep or whole mantle magma ocean.

H₂O and CO₂ will be integrated in solidifying minerals in small quantities, will be enriched in solution in magma ocean liquids as solidification proceeds, and will degas into a growing steam atmosphere. At pressures and temperatures of magma ocean crystallization no hydrous or carbonate minerals will crystallize (Ohtani et al., 2004; Wyllie and Ryabchikov, 2000). Details of the solidification process, the mineral considered, their H₂O and carbon partitioning, and other methods can be found in Elkins-

<caption>**Table 2**</caption>

Modelled atmospheric partial surface pressures $P_{\text{H}_2\text{O}}$ and $P_{\text{CO}_2}$ in units of bar of catastrophically outgassed steam atmospheres dependent on initial H₂O and CO₂ contents in wt% inside a magma ocean with a minimum depth of 500 km and a maximum depth of 2000 km.

<table><thead><tr><th>Bulk magma ocean</th><th>Initial H<sub>2</sub>O (wt.%)</th><th>Initial CO<sub>2</sub> (wt.%)</th><th>P<sub>H<sub>2</sub>O</sub> (bar)</th><th>P<sub>CO<sub>2</sub></sub> (bar)</th></tr></thead><tbody><tr><td rowspan="2">500 km deep</td><td>0.1</td><td>0.02</td><td>52</td><td>11</td></tr><tr><td>0.2</td><td>0.04</td><td>108</td><td>22</td></tr><tr><td rowspan="2">2000 km deep</td><td>0.1</td><td>0.02</td><td>122</td><td>26</td></tr><tr><td>0.2</td><td>0.04</td><td>257</td><td>54</td></tr></tbody></table>

Tanton (2008). The quantity of water and carbon compounds available for degassing is dependent upon the bulk composition of the magma ocean. The terrestrial planets are likely to have been accreted from chondritic material and planetesimals built from chondrites.

Alexander et al. (2012) recently demonstrated that Earth's water, and therefore likely Mars' water, originated mainly from rocky meteoritic material. Wood (2005) reports up to 20 wt% of H₂O in primitive undifferentiated chondrites, and Jarosewich (1990) reports ~3 wt% H₂O in achondrites, though most are drier. Enstatite chondrites match the oxygen isotope composition of the Earth, but smaller fractions of the wide compositional range of other meteorite compositions (see also Alexander et al., 2012 and Drake and Righter, 2002; and references therein) though volatile-rich material from greater radii in the planetary disk may have been added later in planetary formation (e.g., Raymond et al., 2006; O'Brien et al., 2006). Here we assume that water and carbon are added to the growing rocky planets from rocky chondritic material.

Though the original quantity of water and carbon added during giant impacts remains unconstrained, we model two possible starting compositions, according to Brasser (2013) one with 1000 ppm H₂O, and one with 2000 ppm H₂O, each with one-fifth the CO₂ content. These initial compositions are conservatively supported by the data of Jarosewich (1990). For simplicity the carbon is assumed to be degassed as CO₂, though reducing conditions may have produced CO or even CH₄.

Elkins-Tanton (2008) showed that for a range of magma ocean bulk compositions with between ~500 and 5000 ppm H₂O, between ~70% and ~99% of the initial water and carbon is degassed into the planetary atmosphere. Magma ocean solidification is therefore the most significant degassing event in a planet's evolution; the remainder of the volatiles are stored in the interior, available for later degassing via volcanic processes (e.g., Grott et al., 2011).

Table 2 shows the partial surface pressures of catastrophically outgassed steam atmospheres, depending on the assumed bulk magma ocean depths and the initial H₂O and CO₂ contents in the magma ocean in wt.% according to the model of Elkins-Tanton (2008). One can see that a global magma ocean with the depth of ~500 km can produce a steam atmosphere with total surface pressures of ~60–130 bar. If the magma ocean contained the whole mantle, surface pressures between ~150 and 310 bar could have been outgassed.

## 3. Radiation environment during Mars' initial life time

The efficiency of thermal atmospheric escape is related to the planet's temperature at the base of the thermosphere which is located near the mesopause–homopause location in combination with the amount of the XUV flux that is absorbed in the upper atmosphere. The predicted evolution of the Sun's bolometric luminosity relative to its present value and the related equilibrium
<figure><img src="image_5.png" /><figcaption>Fig. 3. Evolution of the Sun's bolometric luminosity relative to its present value (upper panel) and the equilibrium temperature of Mars (lower panel). The solid line corresponds to an evolution track of Tognelli et al. (2011) and the dotted line to Baraffe et al. (1998), both for a star of solar mass and metallicity. The present-day values in both panels are indicated by diamonds. For the evolution of the T<sub>eq</sub>, a constant albedo was assumed.</figcaption></figure>

temperature T<sub>eq</sub> at Mars is shown in Fig. 3. We have chosen two stellar evolution tracks (Baraffe et al., 1998; Tognelli et al., 2011) which predict the lowest and highest luminosities, respectively, between 1 and 10 Myr compared to other authors (cf. Fig. 14 of Tognelli et al., 2011). From Baraffe et al. (1998), the track with parameters M = 1M<sub>☉</sub>, Y=0.282, Z=0.02 and mixing length parameter α = 1.9 was adopted, the track from Tognelli et al. (2011) has M = 1M<sub>☉</sub>, Y=0.288, Z=0.02 and α = 1.68. For planetary atmospheres that are in long-term radiative equilibrium the so-called planetary skin temperature is T<sub>eff</sub> ≈ T<sub>eq</sub>. The lower panel of Fig. 3 shows the corresponding evolution of the equilibrium temperature of Mars, which is about 200 K, 3–4 Myr after the Sun's origin. We assume a constant albedo over time and adopt a present-day value of T<sub>eq</sub> = 217 K. One should note that variations of the albedo due to the evolution of Mars' early atmospheric composition and the Sun's spectral energy distribution could alter the predicted evolution of T<sub>eq</sub> as shown in Fig. 3. Thermal escape of the martian protoatmosphere was driven by the XUV emission of the young Sun. The evolution of this high-energy emission of a solar-type star can be roughly divided into two regimes, the saturation phase and the post-saturation evolution. During the saturation phase the stellar X-ray flux does not scale with the stellar rotation period and is saturated about 0.1% of the bolometric luminosity L<sub>bol</sub> (Pizzolato et al., 2003; Jackson et al., 2012). After the Sun settled on the main sequence and began to spin down from a possibly shorter period to about 2 days due to spin-down via angular momentum loss by the solar wind, the post-saturation phase began. During this phase, the XUV emission of the Sun was determined by its rotation period. A reconstruction of the XUV-evolution during this time period was attempted in the “Sun in Time” program (Güdel, 2007 and references therein). By studying a sample of solar analogs of different ages Ribas et al. (2005) found that the Sun's XUV flux enhancement factor I<sub>XUV</sub> at Earth's orbit in the wavelength range

of 1–1200 Å can be calculated as

$$I_{XUV} = (t/4.56)^{-1.23} \qquad (2)$$

with the age t in Gyr. This relation was calibrated back to an age of 100 Myr corresponding to the youngest solar analog in their sample. However, deviations from this power law are possible during the first few hundred Myr because the stellar rotation periods, which determine the efficiency of the magnetic dynamo and, hence, the XUV emission during this phase, are not unique.

The XUV-evolution during the saturation phase was, as mentioned above, mainly determined by the evolution of L<sub>bol</sub>. Due to the difficulty of observing stars in the EUV because of the strong absorption by the interstellar medium, much of what is known about the high-energy emission of very young stars is extrapolated from X-ray observations. Between the zero-age main sequence (ZAMS), which the Sun reached at an age of about 50 Myr according to stellar evolution models (e.g. Baraffe et al., 1998; Siess et al., 2000), and the end of the saturation phase, the solar XUV flux should have been approximately constant because of the more or less constant bolometric luminosity. For pre-main sequence (PMS) stars, the observed X-ray luminosities are in the order of a few 10<sup>30</sup> erg s<sup>-1</sup> and show a large spread of more than an order of magnitude (Preibisch and Feigelson, 2005; Telleschi et al., 2007). These values are nevertheless consistent with the saturation level of main-sequence stars mentioned above because of the more luminous PMS-Sun and the observed evolution of the stellar X-ray emission during the first tens of Myr seems to be determined mainly by changes of L<sub>bol</sub> (Preibisch and Feigelson, 2005; Briggs et al., 2007).

The estimated past evolution of the Sun's XUV flux, scaled to the orbit of Mars and normalized to the average present solar value of 2 erg cm<sup>-2</sup> s<sup>-1</sup> (scaled from the present value at Earth of 4.64 erg cm<sup>-2</sup> s<sup>-1</sup>; Ribas et al., 2005), is shown in Fig. 4. The solid line indicates the post-saturation evolution after Eq. (2) and the symbols correspond to data from solar analog stars and the Sun. The dotted lines sketch a possible PMS-XUV evolution based on the evolution of L<sub>bol</sub> using theoretical evolutionary tracks for a solar mass star (Baraffe et al., 1998) and assuming that the XUV emission consists mainly of X-rays, so that L<sub>XUV</sub>/L<sub>bol</sub> ≈ L<sub>X</sub>/L<sub>bol</sub> ≈ 10<sup>-3.2 ± 0.3</sup>. The value of the saturation level is adopted from Pizzolato et al. (2003) for stars of about one solar mass. The

<figure><img src="image_6.png" /><figcaption>Fig. 4. Evolution of the Sun's XUV emission normalized to the present value and scaled to the present martian orbit at 1.52 AU. The solid line indicates the evolution during the post-saturation phase (Ribas et al., 2005) with data of solar analogs (black dots) and the Sun indicated. The dotted lines indicate the approximate evolution of the saturated XUV emission estimated by 10<sup>-3.2 ± 0.3</sup>L<sub>bol</sub> (Pizzolato et al., 2003), with the bolometric luminosity taken from stellar evolution tracks of a solar mass star (Baraffe et al., 1998). The shaded area indicates the expected formation time of Mars (Brasser, 2013). The dashed line shows our adopted average XUV value during the Sun's saturation phase.</figcaption></figure>
uncertainties of the Sun's XUV emission before the ZAMS are large because of the dependence of its activity level on the convection zone depth and the rotational history, which in turn depends on the disk locking history. Moreover, the contribution of EUV to the total XUV flux is observationally unconstrained because of strong absorption by the interstellar medium. Therefore we adopt a constant average XUV flux level of about 100 times the present value for our escape rate calculations.

The shaded area indicates the approximate formation time of Mars which occurred during the first few Myr (Brasser, 2013). The inner disk was still present after Mars formed, and the inner planets were still forming. An inner disk would have absorbed a significant fraction of the Sun's XUV radiation until it became optically thin so that the XUV flux actually received by Mars could have been lower than estimated in Fig. 4. Typically, inner disks disperse on timescales within a few Myr to 10 Myr (e.g. Mamajek et al., 2004; Najita et al., 2007; Hillenbrand, 2008).

Thus, if one compares the latest views of Mars' origin and age with that of the radiation history of the young Sun and the nebula dissipation time, Mars' nebula-based and/or outgassed steam atmosphere as well as volatiles which were delivered by frequent impacts were exposed to an XUV flux which was ~100 times stronger compared to that of the present Sun during ~95–100 Myr after the planet's origin. In the following section we investigate how long early Mars could have kept these hydrogen-rich protoatmospheres against XUV-driven thermal atmospheric escape.

## 4. Thermal escape of Mars' protoatmosphere

At present Mars the CO₂-rich thermosphere is in hydrostatic equilibrium, while a hydrogen-rich upper atmosphere of the protoatmosphere that is exposed to the high XUV flux of the young Sun will hydrodynamically expand and the bulk atmospheric particles can escape efficiently (e.g., Watson et al., 1981; Chassefiere, 1996a,b; Tian et al., 2009; Lammer, 2013; Lammer et al., 2013a, 2013b). For this reason we apply a 1-D hydrodynamic upper atmosphere model to the martian protoatmosphere and calculate the XUV-heated hydrogen-dominated dynamically expanding upper atmosphere structure and the thermal hydrogen escape rates, including dissociated and dragged heavier atmospheric main species.

### 4.1. Energy absorption and model description

The thermosphere is heated due to the absorption, excitation, dissociation and ionization of the gas by the incoming solar XUV radiation. By averaging the XUV volume heating rate over Mars' dayside the volume heat production rate $q_{XUV}$ due to the absorption of the solar radiation can then be written as (e.g., Erkaev et al., submitted for publication; Lammer et al., 2013b)

$$q(t, r) = \frac{\eta n \sigma_a}{2} \int_{0}^{\pi/2 + \arccos(1/r)} J(t, r, \Theta) \sin \Theta \, d\Theta, \quad (3)$$

with the polar angle $\Theta$ and $J(t, r, \Theta) = J_{XUV}e^{-\tau(t,r,\Theta)}$, where

$$\tau(t, r, \Theta) = \int_{r}^{\infty} \cos \Theta \, \sigma_a n \left( t, \sqrt{s^2 + r^2 \sin^2 \Theta} \right) ds, \quad (4)$$

$q$ is the volume heating rate depending on the radial distance, $n$ the atmospheric number density which is a function of time and spherical radius $r$, $\eta$ the heating efficiency which corresponds to the fraction of absorbed XUV radiation which is transformed into thermal energy. Depending on the availability of IR-cooling molecules such as $H_3^+$ or CO₂ it is known from various studies that $\eta \sim 15–60\%$ (Chassefiere, 1996a,b; Yelle, 2004; Lammer et al.,

2009; Leitzinger et al., 2011; Koskinen et al., in press). $\sigma_a$ is the absorption cross-section of hydrogen, and $J_{XUV}$ is the XUV energy flux of the young Sun outside the protoatmosphere.

For studying the XUV-exposed structure of the upper atmosphere we solve the system of the 1-D fluid equations for mass, momentum, and energy conservation in spherical coordinates by applying a non-stationary 1D hydrodynamic upper atmosphere model which is described in detail in Erkaev et al. (submitted for publication)

$$\frac{\partial \rho r^2}{\partial t} + \frac{\partial \rho v r^2}{\partial r} = 0, \quad (5)$$

$$\frac{\partial \rho v r^2}{\partial t} + \frac{\partial [r^2(\rho v^2 + P)]}{\partial r} = \rho g r^2 + 2Pr, \quad (6)$$

$$\frac{\partial r^2 \left[ \frac{\rho v^2}{2} + \frac{P}{\gamma - 1} \right]}{\partial t} + \frac{\partial v r^2 \left[ \frac{\rho v^2}{2} + \frac{\gamma P}{\gamma - 1} \right]}{\partial r} = \rho v r^2 g + q_{XUV} r^2, \quad (7)$$

with pressure

$$P = \frac{\rho}{m_H} kT, \quad (8)$$

and gravitational acceleration

$$g = -\nabla \Phi. \quad (9)$$

We note that we neglect the conduction term in the equations because as shown later the energy flux related to thermal conductivity is less important under these extreme conditions compared to the energy flux of the hydrodynamic flow. Here, $\rho$, $v$, $P$ and $T$ are the mass density, radial velocity, pressure and temperature of the atmosphere, $r$ is the radial distance from the center of the planet, $m_H$ is the mass of atomic hydrogen, $G$ is Newton's gravitational constant, $\gamma$ is the polytropic index or the ratio of the specific heats, and $k$ is the Boltzmann constant.

For computational convenience we introduce normalized parameters

$$\tilde{P} = P/(n_0 k T_0), \quad \tilde{\rho} = \rho/(n_0 m),$$

$$\tilde{v} = v/v_0, \quad v_0 = \sqrt{kT_0/m}, \quad \tilde{T} = T/T_0,$$

$$\tilde{q} = qr_0/(mn_0v_0^3), \quad \tilde{r} = r/r_0,$$

$$\tilde{t} = tv_0/r_0, \quad \beta = GmM_{pl}/(r_0kT_0). \quad (10)$$

Here $r_0$, $T_0$, $n_0$ and $v_0$ are the radius, temperature, number density and thermal velocity at the lower boundary of the simulation domain. $\beta$ is the so-called Jeans parameter (Chamberlain, 1963). For values of $\beta > 30$ the atmosphere can be considered as bound to the planet. For values which are lower classical Jeans escape happens. For $\beta$ values that are ~2–3.5 the thermal escape can be very high (Volkov et al., 2011; Volkov and Johnson, 2013) and for values $\le 1.5$ classical blow-off occurs and the atmosphere escapes uncontrolled. Using normalizations (10), we obtain the normalized XUV flux distribution in the planetary atmosphere

$$\tilde{J}(\tilde{r}, \Theta) = J/J_{XUV0} = \exp[-\tilde{\tau}(\tilde{r}, \Theta)], \quad (11)$$

where

$$\tilde{\tau}(\tilde{r}, \Theta) = \int_{\tilde{r}}^{\infty} \cos \Theta \, a\tilde{n} \left( \tilde{t}, \sqrt{s^2 + \tilde{r}^2 \sin^2 \Theta} \right) ds, \quad (12)$$
where $a = \sigma_a n_0 r_0$ is obtained due to the normalization of Eq. (4). The normalized heating rate is given by

$$ \tilde{q}(\tilde{r}) = A \tilde{n} \int_{0}^{\pi/2 + \arccos(1/\tilde{r})} \exp[-\tilde{r}(\tilde{r}, \Theta)] \sin \Theta d\Theta, \quad (13) $$

Integrating (13) over the whole domain we obtain the total energy absorption in the normalized units which is proportional to the incoming XUV flux.

$$ \int_{1}^{\infty} \tilde{q} 4\pi \tilde{r}^2 d\tilde{r} = \pi \frac{J_{XUV}}{mn_0 v_0^3} \frac{r_{XUV,eff}^2}{r_0^2}, \quad (14) $$

where $r_{XUV,eff}$ is the effective radius of the XUV energy absorption which is dependent on the density distribution. This effective radius can be determined from the following equation:

$$ r_{XUV,eff}^2/r_0^2 = 1 + 2 \int_{1}^{\infty} [1 - \tilde{J}(s, \pi/2)]s \, ds. \quad (15) $$

As shown by Watson et al. (1981) the effective radius can exceed the planetary radius quite substantially for a planetary body, which has a low gravity field and hence in low values of the $\beta$ parameter when its atmosphere is exposed by high XUV fluxes. We get the appropriate coefficient

$$ A = \frac{\eta \sigma_a r_0 J_{XUV}}{2 m v_0^3} \quad (16) $$

in formula (13) to satisfy Eq. (14) for a given value of $J_{XUV}$.

## 4.2. Boundary conditions at the lower thermosphere

The boundary conditions at the lower boundary of our simulation domain are the gas temperature $T_0$, number density $n_0$ and the corresponding thermal velocity $v_0$ near the mesopause-homopause level $r_0$, that is at present martian conditions located near the base of the thermosphere. The value of the number density $n_0$ at the base of the thermosphere can never be arbitrarily increased or decreased as much as by an order of magnitude, even if the surface pressure on a planet varies during its life time by many orders of magnitude. The reason for this is that the value of $n_0$ is strictly determined by the XUV absorption optical depth of the thermosphere. The temperature $T_0$ at the base of the thermosphere $z_0 = (r_0 - r_{pl})$ is determined only by the variation of the equilibrium or skin temperature of a planet, to which the base temperature $T_0$ is usually quite close. In a hotter environment corresponding to the catastrophically outgassed steam atmosphere, which is for instance strongly heated by frequent impacts, $z_0$ and the above estimated XUV effective radius $r_{XUV}$ simply rises to a higher altitude where the base pressure retains the same constant value as in a less dense atmosphere.

Marcq (2012) studied with a 1-D radiative-convective atmospheric model the coupling between magma oceans and outgassed steam atmospheres and found that for surface temperatures $T_s \ge 2350$ K, the radiative temperature of a planetary atmosphere $T_{eff}$ can rise from ~230 K to ~300–400 K, while $T_{eq}$ remains close to ~200 K. However, such extreme surface temperatures are only be reached during the totally and partially molten stage of the magma ocean, which last only for ≈0.1 Myr (Lebrun et al., 2013). For this reason we assume in the following thermal escape calculations similar as in Fig. 3 a temperature $T_0$ of 200 K at the base of the thermosphere which corresponds to the equilibrium $T_{eq}$, or skin temperature of Mars' orbit. We point out that an uncertainty of ±20 K as shown in the evolutionary path of $T_{eq}$ in Fig. 3 does not have a big influence in the modeled escape rates. We assume an atomic hydrogen density $n_0 = 10^{13} \text{ cm}^{-3}$ at the lower boundary of the hydrogen-rich protoatmosphere (e.g., Kasting and Pollack, 1983; Tian et al., 2005). According to Tian et al. (2005), similar number density values can be expected also

to H$_2$O mixing ratios ≥50% in a humid steam-like terrestrial planetary atmosphere.

The upper boundary of our simulation domain is chosen at $70r_{pl}$, but the results of our hydrodynamic model are considered as accurate only until the Knudsen number $Kn$, which is the ratio between the mean free path and the scale height, reaches 0.1 (Johnson et al., 2013). Because of the high XUV flux the whole bulk atmosphere reaches the martian escape velocity below or at this altitude level.

The high XUV flux of the young Sun will dissociate most H$_2$ and H$_2$O molecules in the thermosphere so that the upper part of the studied protoatmospheres should be mainly dominated by hydrogen atoms (Kasting and Pollack, 1983; Chassefière, 1996a; Yelle, 2004; Koskinen et al., 2010; Lammer, 2013). As it was shown by Marcq (2012), during periods of magma ocean related hot surface temperatures the tropopause location in an overlaying steam atmosphere can move at an Earth or Venus-like planet from its present altitude of ~30–40 km up to higher altitudes of ~300–550 km. Depending on the surface temperature and pressure of the steam atmosphere in such an environment the mesopause level would then also move to higher altitudes. By applying the model of Marcq (2012) to the outgassed steam atmospheres given in Table 2, we obtain mesopause altitudes of ~330–350 km, ~450–465 km, ~610–630 km and ~750–850 km for surface temperatures of ~1500 K, ~2000 K, ~2500 K and ~3000 K, respectively. This mesopause altitudes have been estimated by detailed modeling of the lowest 600 km of the steam atmospheres. The altitudes above 600 km are obtained from an extrapolation with a precision of ~20 km in the 600–700 km range and ~50 km above 700 km. The simulations used a grey approximation for the radiative transfer which can influence the profile by setting the mesospheric temperature and thus scale height to a slightly different value, but we do not expect that this uncertainty changes these altitudes by more than 20 km. As one can see, even in the most extreme case with a surface temperature of 3000 K, the mesopause altitude lies below 1000 km for a body with Earth's gravity. However, it will most likely be higher than 1000 km with a lower gravity such as Mars'. We plan to study the response to the mesopause location and its influence in the escape of outgassed steam atmospheres on Mars in detail in the near future.

However, for illustrating the importance and influence of the mesopause location in the escape efficiency we modeled also a case where we assumed that $z_0$ is located at 1000 km above the planet's surface. That hydrogen-dominated gas envelopes with hot surface temperatures will have larger radii compared to planets with present time atmospheres is also addressed in Mordasini et al. (2012). However, the planetary mass-radius relationship model results for small and low mass bodies remain highly uncertain.

# 5. Results

## 5.1. Thermospheric profiles and escape rates

By exposing the martian protoatmospheres with a 100 times higher XUV flux compared to today's solar value in martian orbit, we find that the convective thermal energy flux is less significant than the thermal energy flux related to the hydrodynamic flow. Fig. 5 compares the thermal energy flux due to the hydrodynamic flow (curves at the top: dotted lines: $\eta = 15\%$; dashed-lines: $\eta = 40\%$) per steradian of the atmospheric particles with the convective thermal energy flux (curves at the bottom: dotted lines: $\eta = 15\%$; dashed-lines: $\eta = 40\%$), obtained by our hydrodynamic model. The two sudden decreases in the convective thermal energy flux curves can be explained, because this flux is
proportional to the temperature gradient, and therefore it decreases in the vicinity of the temperature maximum and minimum. At first point we have a strong temperature maximum, and at the second point we have shallow temperature minimum. By comparing the two fluxes one can conclude that under such extreme conditions the influence of the thermal conduction on the atmospheric escape is expected to be rather small. Therefore we neglect the thermal conduction term in the energy equation.

Fig. 6 shows examples of the XUV volume heating rate and the corresponding upper atmosphere structure of a hydrogen dominated upper atmosphere of early Mars for a heating efficiency η of 15 and 40% with T₀ = 200 K, and n₀ = 10¹³ cm⁻³, which is exposed

to a XUV flux which is 100 times higher at the planet's orbit compared to that of the present Sun and assumed mesopause locations at 100 km and 1000 km. Under these assumptions the bulk atmosphere reaches the escape velocity v<sub>esc</sub> at about 35r₀ and 24r₀ for heating efficiencies η of 15% and 40%, respectively. One can see from the volume heating rate q<sub>XUV</sub> and the connected temperature profile that the XUV deposition peak occurs above 1.5r₀ for z₀ = 100 km and at ~2R₀ if z₀ = 1000 km. This can also be seen in the temperature profiles, which decrease first due to adiabatic cooling until the high XUV flux of the young Sun balances the cooling process due to the XUV heating, resulting in the more or less constant temperature profile between ~5 and 35r<sub>pl</sub> of ~50–70 K. One can also see that for a heating efficiency η of 40% the adiabatic cooling is stronger at distances that are ≤2.0r₀. The corresponding temperature drop is also larger for an η of 40% compared to that of 15%. For larger distances r > 2r₀, the energy absorption is larger and in the case of 40% efficiency, the additional heating exceeds the cooling. Therefore, the temperature decrease is less pronounced for large distances in the case of higher heating efficiencies compared to the lower value of η = 15%.

Table 3 shows the thermal hydrogen atom escape rates and relevant atmospheric parameters at the critical distance where the bulk atmosphere reaches sonic speed for a lower and higher heating efficiency η of 15% and 40% and for z₀ at 100 and 1000 km altitude. The temperature T₀ and the number density n₀ is assumed to be 200 K and 10¹³ cm⁻³ in all four cases. One can see from Table 3 that depending on z₀ the thermal hydrogen escape rates can reach values between ~2 × 10³² and ~10³³ H atoms per second. The present time thermal hydrogen atom escape from Mars by the classical Jeans escape is about ~1.5 × 10²⁶ s⁻¹ (e.g., Lammer et al., 2008), which indicates that the thermal escape of hydrogen from Mars' protoatmosphere could have been up to ~6–7 orders of magnitude higher.

<figure><img src="image_7.png" /><figcaption>Fig. 5. Comparison of the thermal energy flux per steradian of the hydrodynamical flow (upper dashed line: η = 40%; dotted line: η = 15%) with the thermal energy flux related only to the thermal conductivity (lower dashed: η = 40%; dotted line: η = 15%).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Fig. 6. Example of the XUV volume heating production rate (top left), temperature profile (top right), density profile (bottom left) and the velocity profile (bottom right) for a hydrogen-rich martian upper atmosphere with z₀=100 km by assuming a heating efficiency of 15% (dotted lines) and 40% (dashed lines) and a temperature T₀ at the base of the thermosphere of 200 K as a function of distance in planetary radii for a hydrogen-dominated upper atmosphere at Mars, that is exposed to a 100 times higher XUV flux compared to today's solar value. The solid line shown in the velocity profiles corresponds to the escape velocity v<sub>esc</sub> as a function of distance. The dashed-dotted lines (η = 15%) and the dashed-dotted-dotted-dotted lines (η = 40%) correspond to similar profiles but with z₀ = 1000 km. The hydrogen atoms reach the escape velocity below the theoretical exobase level at a location of ~35r₀ for η = 15% and at ~24 r<sub>pl</sub> for η = 40%.</figcaption></figure>
<caption>**Table 3**</caption>

Modeled atmospheric parameters and thermal hydrogen atom escape rates $L_{th}$ corresponding to a 100 times higher XUV flux compared to today's Sun at the critical distance $r_c \le r_{exo}$, where the dynamically outward flowing hydrogen dominated bulk atmosphere reaches (sonic speed) above the planetary surface and two heating efficiencies $\eta$ of 15% and 40%.

<table><thead><tr><th>CASES</th><th>η (%)</th><th>z<sub>0</sub> (km)</th><th>r<sub>XUV<sub>eff</sub></sub> (r<sub>0</sub>)</th><th>r<sub>c</sub> (r<sub>0</sub>)</th><th>n<sub>c</sub> (cm<sup>-3</sup>)</th><th>T<sub>c</sub> (K)</th><th>L<sub>th</sub> (s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>CI</td><td>15</td><td>100</td><td>3.4</td><td>32.5</td><td>2 × 10<sup>6</sup></td><td>40</td><td>1.8 × 10<sup>32</sup></td></tr><tr><td>CII</td><td>15</td><td>1000</td><td>4.5</td><td>30</td><td>6.2 × 10<sup>6</sup></td><td>36</td><td>7.0 × 10<sup>32</sup></td></tr><tr><td>CIII</td><td>40</td><td>100</td><td>3.2</td><td>21</td><td>6.5 × 10<sup>6</sup></td><td>60</td><td>3.0 × 10<sup>32</sup></td></tr><tr><td>CIIV</td><td>40</td><td>1000</td><td>4.2</td><td>20</td><td>1.7 × 10<sup>7</sup></td><td>50</td><td>1.0 × 10<sup>33</sup></td></tr></tbody></table>

## 5.2. Escape of the nebula captured hydrogen envelope

By knowing the escape rate of hydrogen atoms we can now estimate the loss of the expected nebula-based hydrogen envelope from proto-Mars. If we use the most massive captured hydrogen envelope shown in Table 1 of ~6.5 × 10²² g, corresponding to a luminosity of ~2.4 × 10²² erg s⁻¹ and a dust grain depletion factor f of 0.1, the envelope would be lost during ~1.3–7.5 Myr. The escape time span depends on the heating efficiency and the location distance of the lower thermosphere. A more realistic captured atmosphere with a mass of ~5 × 10²¹ g would be lost in ~0.1–0.5 Myr. From these escape estimates one can conclude that a captured nebular-based hydrogen envelope should have been lost very fast from the planet after the nebula dissipated. If the radius $r_0$ in the nebula captured hydrogen envelope was at further distances compared to our assumed values, then the escape rates would be higher.

## 5.3. Escape of the catastrophically outgassed steam atmosphere

According to the outgassing of the magma ocean depth dependent steam atmospheres shown in Table 2, even the deepest and most volatile-rich case completes solidification and degassing in ≤ 2 × 10⁵ years. The heat loss from the small planetary body is fast enough to allow rapid solidification in a convecting magma ocean. Theoretical studies by Elkins-Tanton (2008) showed that one can expect that the volatiles are likely to be released toward the end of solidification of the magma ocean in a “burst”.

The applied magma ocean model of Elkins-Tanton (2008) and the related results discussed in Section 2.2 predict a surface temperature of ≥ 800 K at the end of solidification, which lies above the condensation temperature for H₂O of ~645 K. If there is a solid-state mantle overturn, there will be a big temperature jump after ~2–4 Myr, when the hot mantle cumulates rise up in Mars because of their buoyancy, and advect their great heat with them. According to Brasser (2013) Mars' most likely finished its accretion or remained as a planetary embryo when the surrounding nebula was still present around the martian orbit location in its later stages. If this was the case the catastrophically outgassed volatiles could easily build up rapidly around the rocky embryo. As soon as this catastrophically outgassed steam atmosphere was released from the nebula, the efficient escape of the atmosphere which was driven by the high XUV flux of the young Sun began.

According to Lebrun et al. (2013), who studied the thermal evolution of an early martian magma ocean in interaction with a catastrophically outgassed ~43 bar H₂O and ~14 bar CO₂ steam atmosphere, water vapor would start to condense into liquid H₂O after ~0.1 Myr. On the other hand, such a fast cooling of the steam atmosphere contradicts the isotopic analysis of martian SNC meteorites by Debaille et al. (2007), where analyzed data can be

best explained by a progressive crystallization of a magma ocean with a duration of up to ~100 Myr. Therefore, Lebrun et al. (2013) suggest that frequent impacts of large planetesimals and small embryos, which have been not included in their study, could have kept the surface during longer times warmer. This suggestion is quite logical because one can also expect that during the first 100 Myr after the origin of the Solar System the young planets have been frequently hit by large impactors (e.g., Abe and Matsui, 1985, 1988; Genda and Abe, 2005; Lammer et al., 2013a), which may have enhanced the input energy flux above the value which is defined by the solar flux alone. In such a case one will obtain a hotter surface that prevent atmospheric H₂O vapor from condensing (e.g., Hayashi et al., 1979; Genda and Abe, 2005; Lammer et al., 2012; Lammer, 2013; Lebrun et al., 2013).

One should also note that for the surface temperatures of ~500 K, which are expected during the “Mush” stage (Lebrun et al., 2013), according to Kasting (1988) one can also expect water vapor mixing ratios at the mesopause level near to 1. For this reason H₂O will continue to escape effectively, even if there are periods of liquid water on the planet's surface. However, the mesopause level will be closer to the planet's surface and the escape rates will be reduced and may have values which correspond to case CI in Table 3.

In the outgassed steam atmosphere, the H₂O molecules in the upper atmosphere will be dissociated by the high XUV flux of the young Sun and by frequently occurring impacts in the lower thermosphere (e.g., Chassefière, 1996b; Lammer et al., 2012; Lammer, 2013). Tian et al. (2009) showed that for XUV fluxes which are > 10 times that of today's Sun, CO₂ and/or CH₄ molecules in the martian upper atmosphere will also be destroyed, so that C atoms can escape similar to O atoms with escape flux values which are ≥ 10¹¹ cm⁻² s⁻¹. From this study one can expect that for an XUV flux which is ~100 times stronger than the present solar value most CO₂ and/or CH₄ molecules will be dissociated as soon as they are exposed to the high XUV radiation. Therefore, one can assume that O and C atoms should also populate the lower hydrogen dominated thermosphere so that they can be dragged by the dynamically outward flowing hydrogen atom flux (Zahnle and Kasting, 1986; Zahnle et al., 1990; Chassefière, 1996a,b; Hunten et al., 1987; Lammer et al., 2012, 2013a).

As initial amount and composition of the outgassed atmosphere we adopt the four cases presented in Table 2. With the given partial surface pressures of H₂O and CO₂ and assuming that all molecules are dissociated under the high XUV flux of the young Sun, we calculate the initial inventories of atomic H, O, and C. For all four cases atomic hydrogen is the most abundant species ($N_H/N = 0.61$), followed by oxygen ($N_O/N = 0.36$), whereas C is just a minor constituent ($N_C/N = 0.03$). Hydrogen is assumed to escape at rates given in Table 3. The fractionation factors $x_i = L_i/(L_H f_i)$ for an escaping atmosphere composed of two major (here H, O) and several minor species (here only one, namely C) are given by Eqs. (19) and (35) of Zahnle et al. (1990), where $f_i = n_i/n_H = N_i/N_H$ is the mixing ratio with respect to H and $L_i$ are the escape fluxes of the heavy species i given in s⁻¹. Using the definition of x the escape fluxes of O and C can then be written as

$$L_O = L_H f_O x_O = L_H f_O \left( 1 - \frac{\mu_O - 1}{\mu_O \Phi_O} \frac{1}{1 + f_O} \right) \quad (17)$$

$$L_C = L_H f_C \frac{1 - \frac{\mu_C - 1}{\mu_C \Phi_C} + \frac{b_{HC}}{b_{OC}} f_O x_O + \frac{b_{HC}}{b_{HO}} f_O (1 - x_O)}{1 + \frac{b_{HC}}{b_{OC}} f_O} \quad (18)$$

with $\mu_i = m_i/m_H$, the binary diffusion parameters b, and the parameter
$$ \Phi_i = \frac{L_H k T}{3\pi G M m_i b_{Hi}} \qquad (19) $$

which represents approximately the ratio of drag to gravity (drag dominates if $\Phi_i > (\mu_{i-1}/\mu_i)$). The factor $3\pi$ stems from our adopted solid angle over which we assume that escape takes place and which is therefore included in the values of $L_H$. The binary diffusion parameter of O in H $b_{HO} = 4.8 \times 10^{17} T^{0.75} cm^{-1} s^{-1}$ was taken from Table 1 of Zahnle and Kasting (1986). $b_{HC}$ was assumed to be equal to $b_{HO}$, and $b_{OC}$ is roughly estimated as $2 \times 10^{17} T^{0.75} cm^{-1} s^{-1}$. However, we note that changing these parameters, as well as the adopted temperature, does not affect the results if the hydrogen escape rate is large.

Eqs. (17) and (18) were derived under the assumption that the flow is isothermal and subsonic (Zahnle et al., 1990), which is actually not valid during the phase of saturated solar XUV emission studied here. However, they showed that these simpler analytic approximations become comparable to the non-isothermal transonic solutions if $x_i \gg 1/\mu_i$ and $\Phi_i$ is large. These conditions are both fulfilled here because the masses of O and C are much larger than H and hydrogen escapes very efficiently (hence, $\Phi_i \gg$). It was also assumed that the mixing ratios $f_i$ are approximately constant with height. Expressions for $x_i$ without this constraint include terms with an exponential function that goes to zero for large $\Phi_i$ (Zahnle and Kasting, 1986; Zahnle et al., 1990) and would therefore vanish for the cases studied here.

Figs. 7–10 show the temporal evolution of the partial surface pressures of H, O, and C normalized to the initial total surface pressure for the four cases of the outgassed atmospheres given in Table 2. These results have been obtained by adopting the modeled hydrogen loss rates shown in the cases CI, CII, CIII and CIV in Table 3 corresponding to 100 times the present solar XUV flux and a lower boundary temperature $T_0$ of 200 K but low and high heating efficiencies $\eta$ of 15% and 40%, and $z_0$ at 100 km and

1000 km altitude. This temperature is also used for evaluating $\Phi_i$, but choosing a different value does not affect the results because the large $L_H$ dominates. The initial hydrogen inventory evolves with a constant escape rate, because the timescale for total hydrogen loss occurs during a time frame between ~0.4 and 12 Myr, well below the time it takes the Sun to drop out of its saturation phase. The evolution of the inventories, and hence partial surface pressures, of O and C are found numerically by integration of Eqs. (17) and (18).

From these figures one can see that the timescale for complete loss of H, O, and C, for a ~50 bar H₂O and ~10 bar CO₂ atmosphere for low $\eta$ and $z_0$ occurs in less than 2.5 Myr. If the base of the thermosphere would expand from 100 km to 1000 km and $\eta = 40\%$, such a steam atmosphere would be lost after ~0.4 Myr. Depending on the initial volatile content and assumed heating efficiencies and $z_0$, steam atmospheres with ~260 bar H₂O and ~55 bar CO₂ would be lost from early Mars between ~2.1 and 12 Myr. We also note that a magnetosphere would not protect the escape of the bulk atmosphere under these conditions because most of the atoms escape as neutrals until they become ionized due to the interaction with the early solar wind and plasma environment at large planetary distances (Kislyakova et al., submitted for publication; Lammer, 2013).

As discussed above, the time scale for cooling of the steam atmosphere to temperature-pressure values, that water can condense and build lakes or even oceans is very important and the influence of energy deposition on planetary surfaces by frequent impacts of large planetsimals or small embryos have to be studied in coupled magma ocean-protoatmosphere models in the future. Therefore, it is also possible that all of our studied steam atmosphere scenarios presented in Table 2 may have been lost within a few Myr, before the atmospheres cooled to temperatures that big lakes or oceans could have formed.

However, for outgassed steam atmospheres with surface pressures $\gg 50$ bar, the timescale for total escape compared to the

<figure><img src="image_9.png" /><figcaption>Fig. 7. Temporal evolution of the partial surface pressures $P_{surf}$ of H, O, and C normalized to the total initial surface pressure $P_{total}$ for the four compositions of outgassed atmospheres described in Table 2. The hydrogen inventory evolves assuming a constant escape rate and parameters according to CI in Table 3 valid for 100 XUV. Both O and C are dragged along with the escaping H.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 8.** Temporal evolution of the partial surface pressures *P*<sub>surf</sub> of H, O, and C normalized to the total initial surface pressure *P*<sub>total</sub> for the four compositions of outgassed atmospheres described in Table 2. The hydrogen inventory evolves assuming a constant escape rate and parameters according to CII in Table 3 valid for 100 XUV. Both O and C are dragged along with the escaping H.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 9.** Temporal evolution of the partial surface pressures *P*<sub>surf</sub> of H, O, and C normalized to the total initial surface pressure *P*<sub>total</sub> for the four compositions of outgassed atmospheres described in Table 2. The hydrogen inventory evolves assuming a constant escape rate and parameters according to CIII in Table 3 valid for 100 XUV. Both O and C are dragged along with the escaping H.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Fig. 10.** Temporal evolution of the partial surface pressures $P_{surf}$ of H, O, and C normalized to the total initial surface pressure $P_{total}$ for the four compositions of outgassed atmospheres described in Table 2. The hydrogen inventory evolves assuming a constant escape rate and parameters according to CIV in Table 3 valid for 100 XUV. Both O and C are dragged along with the escaping H.</figcaption></figure>

steam atmosphere cooling timescale could be larger, so that large lakes or water oceans could have been formed sporadically. In such scenarios water condensed and could have been present on the planet's surface for short time until the high thermal escape rates and impactors evaporated it again (Genda and Abe, 2005). During this time and also during later stages a fraction of condensed, or via later impacts delivered H$_2$O, may have been again incorporated by hydrothermal alteration processes such as serpentinization, so that remaining parts of it could be stored even today in subsurface serpentine (Chassefière et al., 2013).

We point out that a detailed photochemical study, which includes processes such as dissociation, ionization, etc. of the outgassed CO$_2$ molecules is beyond the scope of the present study. Our expectation that no dense CO$_2$ atmosphere has build up on early Mars during the first 100 Myr is also supported by a study of Tian et al. (2009), who showed that the thermal escape of C atoms was so efficient even during the early Noachian, > 4.1 Gyr ago, that a CO$_2$-dominated martian atmosphere could not have been maintained, and Mars most likely has begun its origin colder. In agreement with Lammer et al. (2013a) by the mid to late Noachian, as one can see from Fig. 4, the solar XUV flux would have become much weaker allowing the build up of a secondary CO$_2$ atmosphere by volcanic outgassing (Grott et al., 2011) and/or impact delivered volatiles.

Our results are also in agreement with the conclusions of Bibring et al. (2005), which are based on the so far not detected carbonates, that no major surface sink of CO$_2$ is present and the initial CO$_2$, if it was more abundant, should have been lost from Mars very early other than being stored in surface reservoirs after having been dissolved in liquid water at the surface. However, it should be noted that the accumulation of a secondary outgassed CO$_2$ atmosphere and volatiles, which could have been delivered by later impacts is highly dependent on less efficient atmospheric escape processes after the strong early hydrodynamic loss during

the XUV-saturation phase of the young Sun as well as by the efficiency of carbonate precipitation, and serpentinization during the Hesperian and Amazonian epochs (e.g., Chassefière and Leblanc, 2011a, 2011b; Lammer et al., 2013a; Niles et al., 2013).

Our result that Mars lost most likely the majority of its initial H$_2$O inventory very early is in support of the hypothesis presented by Albarède and Blichert-Toft (2007) that the planet could not develop an efficient plate tectonic regime due to the rapid removal of water by hydrodynamic escape. These authors suggest that the resulting low abundance of the remaining water in the martian mantle combined with weaker gravity than on Earth acted against the bending and foundering of lithospheric plates and the planet instead took the dynamic route of developing a thick stagnant lid. Because of the low size and gravity of Mars not enough water could be incorporated into the Martian mantle before it was lost to space so that plate tectonics never began.

## 6. Conclusions

The production and loss of the earliest martian atmosphere which consisted of captured nebula gas (H, He, etc.) and outgassed and impact delivered volatiles (e.g. H$_2$O, CO$_2$, and CH$_4$) have been studied. By using the latest knowledge of the origin of Mars summarized in Brasser (2013), we estimated the protoatmosphere masses and partial pressures and applied a 1-D hydrodynamic upper atmosphere model to the extreme XUV conditions of the young Sun. Depending on the amount of the outgassed volatiles, as well as the assumed heating efficiency and altitude location of the lower thermosphere, our results indicate that early Mars lost its nebular captured hydrogen envelope and catastrophically outgassed steam atmosphere most likely within ~0.4–12 Myr by hydrodynamic escape of atomic hydrogen. The main reasons for the fast escape of even a steam atmosphere with an amount of
~70% of an Earth ocean and ~50 bar CO₂ within <12 Myr are Mars' low gravity and the ~100 times higher XUV flux of the young Sun, which lasted ~100 Myr after the Solar Systems origin. The efficient escape of atomic hydrogen, drags heavier atoms within the escaping bulk atmosphere so that they can also be lost to space. Our results support the hypotheses of Tian et al. (2009) that early Mars could not build up a dense CO₂ atmosphere during the early Noachian. The results are also in agreement with the hypothesis presented in Lammer et al. (2013a) that after the planet lost its protoatmosphere the atmospheric escape rates were most likely balanced with the volatiles, which have been outgassed by volcanic activity and delivered by impacts until the activity of the young Sun decreased, so that the atmospheric sources could dominate over the losses ~4.2–3.8 Gyr ago.

## Acknowledgments

P. Odert, H. Lammer, K. G. Kislyakova and Yu. N. Kulikov acknowledge support from the Helmholtz Alliance project "Planetary Evolution and Life". E. Dorfi, M. Güdel, K. G. Kislyakova, H. Lammer, A. Stökl and E. A. Dorfi acknowledge the Austrian Science Fund (FWF) for supporting this study via the FWF NFN project S116 "Pathways to Habitability: From Disks to Active Stars, Planets and Life", and the related FWF NFN subprojects, S 116 02-N1 "Hydrodynamics in Young Star-Disk Systems", S116 604-N16 "Radiation & Wind Evolution from T Tauri Phase to ZAMS and Beyond", and S11607-N16 "Particle/Radiative Interactions with Upper Atmospheres of Planetary Bodies Under Extreme Stellar Conditions". M. Leitzinger and P. Odert also acknowledge the support from the FWF project P22950-N16. N. V. Erkaev acknowledges support by the RFBR Grant no 12-05-00152-a. Finally, H. Lammer thanks M. Ikoma from the Department of Earth and Planetary Science, of the University of Tokyo, Japan, for discussions related to the accumulation of nebular-based hydrogen envelopes around Mars-mass bodies. Finally the authors thank an anonymous referee for the interesting and important suggestions and recommendations that helped to improve the results of our study.

## References

* Abe, Y., 1993. Physical state of very early Earth. Lithosphere 30, 223–235.

* Abe, Y., 1997. Thermal and chemical evolution of the terrestrial magma ocean. Physics of the Earth and Planetary Interiors 100, 27–39.

* Abe, Y., Matsui, T., 1985. The formation of an impact-generated H₂O atmosphere and its implications for the early thermal history of the Earth. Journal of Geophysical Research: Atmospheres 90, 545–559.

* Abe, Y., Matsui, T., 1988. Evolution of an impact-generated H₂O–CO₂ atmosphere and formation of a hot proto-ocean on Earth. Journal of the Atmospheric Sciences 45, 3081–3101.

* Albarède, F., Blichert-Toft, J., 2007. The split fate of the early Earth, Mars, Venus, and Moon. Comptes Rendus Geoscience 339, 917–927.

* Alexander, C.M.O'D., Bowden, R., Fogel, M.L., Howard, K.T., Herd, C.D.K., Nittler, N.R., 2012. The provenances of asteroids, and their contributions to the volatile inventories of the terrestrial planets. Science 337 (6095), 721–723.

* Baker, V.R., 2001. Water and the martian landscape. Nature 412, 228–236.

* Baraffe, I., Chabrier, G., Allard, F., Hauschildt, P.H., 1998. Evolutionary models for solar metallicity low-mass stars: mass-magnitude relationships and color-magnitude diagrams. Astronomy and Astrophysics 337, 403–412.

* Bibring, J.-P., Langevin, Y., Gendrin, A., Gondet, B., Poulet, F., Berthé, M., Soufflot, A., Arvidson, R., Mangold, N., Mustard, J., Drossart, P., 2005. the OMEGA team, 2005. Mars surface diversity as revealed by the OMEGA/Mars Express observations. Science 307, 1576–1581.

* Brasser, R., 2013. The formation of Mars: building blocks and accretion time scale. Space Science Reviews 174, 11–25.

* Briggs, K.R., Güdel, M., Telleschi, A., Preibisch, T., Stelzer, B., Bouvier, J., Rebull, L., Audard, M., Scelsi, L., Micela, G., Grosso, N., Palla, F., 2007. The X-ray activity-rotation relation of T Tauri stars in Taurus-Auriga. Astronomy and Astrophysics 468, 413–424.

* Chamberlain, J.W., 1963. Planetary coronae and atmospheric evaporation. Planetary and Space Science 11, 901–996.

* Chassefière, E., 1996a. Hydrodynamic escape of hydrogen from a hot water-rich atmosphere: the case of Venus. Journal of Geophysical Research: Atmospheres 101, 26039–26056.

* Chassefière, E., 1996b. Hydrodynamic escape of oxygen from primitive atmospheres: application to the cases of Venus and Mars. Icarus 124, 537–552.

* Chassefière, E., Leblanc, F., 2011a. Constraining methane release due to serpentinization by the observed D/H ratio on Mars. Earth and Planetary Science Letters 310, 262–271.

* Chassefière, E., Leblanc, F., 2011b. Explaining the redox imbalance between the H and O escape fluxes at Mars by the oxidation of methane. Planetary and Space Science 59, 218–226.

* Chassefière, E., Langlais, B., Quesnel, Y., Leblanc, F., 2013. The fate of early Mars' lost water: The role of serpentinization. Journal of Geophysical Research: Atmospheres 118, 1123–1134.

* Claire, M.W., Sheets, J., Cohen, M., Ribas, I., Meadows, V.S., Catling, D.C., 2012. The evolution of solar flux from 0.1 nm to 160 µm: Quantitative estimates for planetary studies. Astrophysical Journal 757, art. id. 95.

* Debaille, V., Brandon, A.D., Yin, Q.Z., Jacobsen, B., 2007. Coupled $^{142}$Nd-$^{143}$Nd evidence for a protracted magma ocean in Mars. Nature 450, 525–528.

* Drake, M.J., Righter, K., 2002. Determining the composition of the Earth. Nature 416, 39–44.

* Elkins-Tanton, L.T., Zaranek, S., Parmentier, E.M., 2005. Early magnetic field and magmatic activity on Mars from magma ocean overturn. Earth and Planetary Science Letters 236, 1–12.

* Elkins-Tanton, L.T., 2008. Linked magma ocean solidification and atmospheric growth for Earth and Mars. Earth and Planetary Science Letters 271, 181–191.

* Elkins-Tanton, L.T., 2011. Formation of water ocean on rocky planets. Astrophysics and Space Science 332, 359–364.

* Erkaev, N.V., Lammer, H., Odert, P., Kulikov, Yu.N., Kislyakova, K.G., Khodachenko, M. L., Güdel, M., Hanslmeier, A., Biernat, H. XUV exposed non-hydrostatic hydrogen-rich upper atmospheres of Earth-size planets. Part I: atmospheric expansion and thermal escape. Astrobiology, submitted for publication. http://arxiv.org/abs/1212.4982.

* Freedman, R.S., Marley, M.S., Lodders, K., 2008. Line and mean opacities for ultracool dwarfs and extrasolar planets. Astrophysical Journal 174, 504–513.

* Foley, C.N., Wadhwa, M., Borg, L.E., Janney, P.E., Hines, R., Grove, T.L., 2005. The early differentiation history of Mars from new constraints on $^{182}$W-$^{142}$Nd isotope systematics in the SNC meteorites. Geochimica et Cosmochimica Acta 69, 4557–4571.

* Genda, H., Abe, Y., 2005. Enhanced atmospheric loss on protoplanets at the giant impact phase in the presence of oceans. Nature 433, 842–844.

* Grott, M., Morschhauser, A., Breuer, D., Hauber, E., 2011. Volcanic outgassing of CO₂ and H₂O on Mars. Earth and Planetary Science Letters 308, 391–400.

* Güdel, M., Guinan, E.F., Skinner, S.L., 1997. The X-ray Sun in time: a study of the long-term evolution of coronae of solar-type stars. Astrophysical Journal 483, 947–960.

* Güdel, M., 2007. The Sun in time: activity and environment. Living Reviews in Solar Physics 4 (3).

* Harper, C.L., Nyquist, L.E., Bansal, B., Wiesmann, H., Shih, C.-Y., 1995. Rapid accretion and early differentiation of Mars indicated by $^{142}$Nd/$^{144}$Nd in SNC meteorites. Science 267, 213–217.

* Hamano, K., Abe, Y., Genda, H., 2013. Emergence of two types of terrestrial planet on solidification of magma ocean. Nature 497, 607–610.

* Hayashi, C., Nakazawa, K., Mizuno, H., 1979. Earth's melting due to the blanketing effect of the primordial dense atmosphere. Earth and Planetary Science Letters 43, 22–28.

* Hayashi, C., 1981. Structure of the solar nebula. Growth and decay of magnetic fields, and effect of magnetic and turbulent viscosities on the nebula. Progress of Theoretical Physics Supplement 70, 35–53.

* Hillenbrand, L.A., 2008. Observational constraints on dust disk lifetimes: implications for planet formation. In: Livio, M., Sahu, K., Valenti, J. (Eds.), A Decade of Extrasolar Planets around Normal Stars. Space Telescope Science Institute Symposium Series, vol. 19. pp. 84–105.

* Hunten, D.M., Pepin, R.O., Walker, J.C.G., 1987. Mass fractionation in hydrodynamic escape. Icarus 69, 532–549.

* Horner, J., Mousis, O., Petit, J.-M., Jones, B.-W., 2009. Differences between the impact regimes of the terrestrial planets: implications for primordial D:H ratios. Planetary and Space Science 57, 1338–1345.

* Ikoma, M., Nakazawa, K., Emori, H., 2000. Formation of giant planets: dependences on core accretion rate and grain opacity. Astrophysical Journal 537, 1013–1025.

* Ikoma, M., Genda, H., 2006. Constraints on the mass of a habitable planet with water of nebular origin. Astrophysical Journal 648, 696–706.

* Jackson, A.P., Davis, T.A., Wheatley, P.J., 2012. The coronal X-ray-age relation and its implications for the evaporation of exoplanets. Monthly Notices of the Royal Astronomical Society 422, 2024–2043.

* Jaroszewich, E., 1990. Chemical analyses of meteorites: a compilation of stony and iron meteorite analyses. Meteoritics 25, 323–337.

* Johnson, R.E., Volkov, A.N., Erwin, J.T., 2013. Molecular-kinetic simulations of escape from the ex-planet and exoplanets: criterion for transonic flow. Astrophysical Journal Letters 768, L4.

* Kasting, J.F., Pollack, J.B., 1983. Loss of water from Venus. I. Hydrodynamic escape of hydrogen. Icarus 53, 479–509.

* Kasting, J.F., 1988. Runaway and moist greenhouse atmospheres and the evolution of Earth and Venus. Icarus 74, 472–494.

* Kislyakova, G.K., Lammer, H., Holmström, M., Panchenko, M., Khodachenko, M.L., Erkaev, N.V., Odert, P., Kulikov, Yu. N., Leitzinger, M., Güdel, M., Hanslmeier, A.
XUV exposed, non-hydrostatic hydrogen-rich upper atmospheres of Earth-size planets. II: hydrogen coronae and ion escape. Astrobiology, submitted for publication. http://arxiv.org/abs/1212.4710.

Koskinen, T.T., Yelle, R.V., Lavvas, P., Lewis, N.K., 2010. Characterizing the thermosphere of HD209458 b with UV tranist observations. Astrophysical Journal 723, 116–128.

Koskinen, T.T., Harris, M.J., Yelle, R.V., Lavvas, P. The escape of heavy atoms from the ionosphere of HD209458b. I. A photochemical-dynamical model of the thermosphere. Icarus, http://dx.doi.org/10.1016/j.icarus.2012.09.027, in press.

Kuhfuß, R., 1987. Ein Modell für zeitabhängige Konvektion in Sternen. Ph.D. Thesis, TU München.

Lammer, H., Kasting, J.F., Chassefière, E., Johnson, R.E., Kulikov, Yu. N., Tian, F., 2008. Atmospheric escape and evolution of terrestrial planets and satellites. Space Science Reviews 139, 399–436.

Lammer, H., Odert, P., Leitzinger, M., Khodachenko, M.L., Panchenko, M., Kulikov, Yu.N., Zhang, T.L., Lichtenegger, H.J.M., Erkaev, N.V., Wuchterl, G., Micela, G., Penz, T., Biernat, H.K., Weingrill, J., Steller, M., Ottacher, H., Hasiba, J., Hanslmeier, A., 2009. Determining the mass loss limit for close-in exoplanets: What can we learn from transit observations?. Astronomy and Astrophysics 506, 399–410.

Lammer, H., Kislyakova, K.G., Odert, P., Leitzinger, M., Schwarz, R., Pilat-Lohinger, E., Kulikov, Yu.N., Khodachenko, M.L., Güdel, M., Hanslmeier, A., 2012. Pathways to Earth-like atmospheres: extreme ultraviolet (EUV)-powered escape of hydrogen-rich protoatmospheres. Origins of Life and Evolution of Biospheres 41, 503–522.

Lammer, H., 2013. Origin and Evolution of Planetary Atmospheres: Implications for Habitability. Springer Briefs in Astronomy, Springer Publishing House, Heidelberg/New York.

Lammer, H., Chassefière, E., Karatekin, Ö., Morschhauser, A., Niles, P.B., Mousis, O., Odert, P., Möstl, U.V., Breuer, D., Dehant, V., Grott, M., Gröller, H., Hauber, E., San Pham, Lê Binh, 2013a. Outgassing history and escape of the martian atmosphere and water inventory. Space Science Reviews 174, 113–154.

Lammer, H., Erkaev, N.V., Odert, P., Kislyakova, K.G., Leitzinger, M., 2013b. Probing the blow-off criteria of hydrogen-rich "super-Earths". Monthly Notices of the Royal Astronomical Society 430, 1247–1256.

Lebrun, T., Massol, H., Chassefière, E., Davaille, A., Marcq, E., Sarda, P., Leblanc, F., Brandeis, G., 2013. Thermal evolution of an early magma ocean in interaction with the atmosphere. Journal of Geophysical Research: Atmospheres 118, 1–22.

Leitzinger, M., Odert, P., Kulikov, Yu. N., Lammer, H., Wuchterl, G., Penz, T., Guarcello, M.G., Micela, G., Khodachenko, M.L., Weingrill, J., Hanslmeier, A., Biernat, H.K., Schneider, J., 2011. Could CoRoT-7b and Kepler-10b be remnants of evaporated gas or ice giants? Planetary and Space Science 59, 1472–1481.

Lunine, J.I., Chambers, J., Morbidelli, A., Leshin, L.A., 2003. The origin of water on Mars. Icarus 165, 1–8.

Marcq, E., 2012. A simple 1-D radiative-convective atmospheric model designed for integration into coupled models of magma ocean planets. Journal of Geophysical Research: Atmospheres 117, E01001.

Matsui, T., Abe, Y., 1986. Evolution of an impact-induced atmosphere and magma ocean on the accreting Earth. Nature 319, 303–305.

Mamajek, E.E., Meyer, M.R., Hinz, P.M., Hoffmann, W.F., Cohen, M., Hora, J.L., 2004. Constraining the lifetime of circumstellar disks in the terrestrial planet zone: A mid-infrared survey of the 30 Myr old tucana-horologium association. Astrophysical Journal 612, 496–510.

Mizuno, H., Nakazawa, K., Hayashi, C., 1978. Instability of a gaseous envelope surrounding a planetary core and formation of giant planets. Progress of Theoretical Physics 60, 699–710.

Mordasini, C., Alibert, Y., Georgy, C., Dittkrist, K.-M., Henning, T., 2012. Characterization of exoplanets from their formation II. The planetary mass-radius relationship. Astronomy and Astrophysics 545, A112.

Nakazawa, K., Mizuno, H., Sekiya, M., Hayashi, C., 1985. Structure of the primordial atmosphere surrounding the early-Earth. Journal of Geomagnetism and Geoelectricity 37, 781–799.

Najita, J.R., Carr, J.S., Glassgold, A.E., Valenti, J.A., 2007. Gaseous inner disks. In: Reipurth, B., Jewitt, D., Keil, K. (Eds.), Protostars and Planets V, pp. 507–522.

Niles, P.B., Catling, D.C., Berger, G., Chassefière, E., Ehlmann, B.L., Michalski, J.R., Morris, R., Ruff, S.W., Sutter, B., 2013. Geochemistry of carbonates on Mars: implications for climate history and nature of aqueous environments. Space Science Reviews 174, 301–328.

O'Brien, D.P., Morbidelli, A., Levison, H.F., 2006. Terrestrial planet formation with strong dynamical friction. Icarus 184, 39–58.

Ohtani, E., Litasov, K., Hosoya, T., Kubo, T., Kondo, T., 2004. Water transport into the deep mantle and formation of a hydrous zone. Physics of the Earth and Planetary Interiors 143–144, 255–269.

Pizzolato, N., Maggio, A., Micela, G., Sciortino, S., Ventura, P., 2003. The stellar activity-rotation relationship revisited: Dependence of saturated and non-saturated X-ray emission regimes on stellar mass for late-type dwarfs. Astronomy and Astrophysics 397, 147–157.

Preibisch, T., Feigelson, E.D., 2005. The evolution of X-Ray emission in young stars. Astrophysical Journal Supplement Series 160, 390–400.

Rafikov, R.R., 2006. Atmospheres of protoplanetary cores: critical mass for nucleated instability. Astrophysical Journal 648, 666–682.

Raymond, S.N., Quinn, T., Lunine, J.I., 2006. High-resolution simulations of the final assembly of Earth-like planets I: Terrestrial accretion and dynamics. Icarus 183, 265–282.

Reese, C.C., Solomatov, V.S., 2006. Fluid dynamics of local martian magma oceans. Icarus 184 (1), 102–120.

Ribas, I., Guinan, E.F., Güdel, M., Audard, M., 2005. Evolution of the solar activity over time and effects on planetary atmospheres. I. High-energy irradiances (1–1700 Å). Astrophysical Journal 622, 680–694.

Sasaki, S., 1990. The primary solar-type atmosphere surrounding the accreting Earth: H₂O-induced high surface temperature. In: Newsom, H.E., Jones, J.H. (Eds.), Origin of the Earth. Oxford Univ. Press, New York, pp. 195–209.

Saumon, D., Chabrier, G., van Horn, H.M., 1995. An equation of state for low-mass stars and giant planets. Astrophysical Journal Supplement Series 99, 713.

Semenov, D., Henning, T., Helling, C., Ilgner, M., Sedlmayr, E., 2003. Rosseland and Planck mean opacities for protoplanetary discs. Astrophysics 410, 611–621.

Siess, L., Dufour, E., Forestini, M., 2000. An internet server for pre-main sequence tracks of low- and intermediate-mass stars. Astronomy and Astrophysics 358, 593–599.

Telleschi, A., Güdel, M., Briggs, K.R., Skinner, S.L., Audard, M., Franciosini, E., 2007. The first high-resolution X-ray spectrum of a Herbig star: AB Aurigae. Astronomy and Astrophysics 468, 425–442.

Tian, F., Toon, O.B., Pavlov, A.A., De Sterck, H., 2005. A hydrogen-rich early Earth atmosphere. Science 308, 1014–1017.

Tian, F., Kasting, J.F., Solomon, S.C., 2009. Thermal escape of carbon from the early Martian atmosphere. Geophysical Research Letters 36, L02205.

Tognelli, D., Prada Moroni, P.G., Degl'Innocenti, S., 2011. The Pisa pre-main sequence tracks and isochrones. A database covering a wide range of Z, Y, mass, and age values. Astronomy and Astrophysics 533, A109.

Tonks, Melosh, 1993. Magma ocean formation due to giant impacts. Journal of Geophysical Research: Atmospheres 98, 5319–5333.

Volkov, A.N., Johnson, R.E., Tucker, O.J., Erwin, J.T., 2011. Thermally driven atmospheric escape: transition from hydrodynamic to Jeans escape. Astrophysical Journal 729, L24.

Volkov, A.N., Johnson, R.E., 2013. Thermal escape in the hydrodynamic regime: Reconsideration of Parker's isentropic theory based on results of kinetic simulations. Astrophysical Journal 765, pp. 21, art. id. 90.

Wänke, H., Dreibus, G., 1988. Chemical composition and accretion history of terrestrial planets. Philosophical Transactions of the Royal Society A 325, 545–557.

Wash, K.J., Morbidelli, A., Raymond, S.N., O'Brien, D.P., Mandell, A.M., 2011. A low mass for mars from Jupiter's early gas-driven migration. Nature 475, 206–209.

Watson, A.J., Donahue, T.M., Walker, J.C.G., 1981. The dynamics of a rapidly escaping atmosphere: applications to the evolution of Earth and Venus. Icarus 48, 150–166.

Wetherill, G.W., 1986. Accumulation of terrestrial planets and implications concerning lunar origin. In: Hartmann, W.K., Phillips, R.J., Taylor, G.J. (Eds.), Origin of the Moon Lunar and Planet. University Arizona Press, Chicago, pp. 519–550.

Wood, J.A., 2005. The chondrite types and their origins. In: Krot, A.N., Scott, E.R.D., Reipurth, B. (Eds.), Chondrites and the Protoplanetary Disk. ASP Conference Series, vol. 341.

Wuchterl, G., 1993. The critical mass for protoplanets revised – Massive envelopes through convection. Icarus 106, 323–334.

Wyllie, P.J., Ryabchikov, I.D., 2000. Volatile components, magmas, and critical fluids in upwelling mantle. Journal of Petrology 41, 1195–1206.

Yelle, R.V., 2004. Aeronomy of extra-solar giant planets at small orbital distances. Icarus 170, 167–179.

Zahnle, K.J., Kasting, J.F., 1986. Mass fractionation during transonic escape and implications for loss of water from Mars and Venus. Icarus 68, 462–480.

Zahnle, K.J., Kasting, J.F., Pollack, J.B., 1988. Evolution of a steam atmosphere during Earth's accretion. Icarus 74, 62–97.

Zahnle, K.J., Kasting, J.F., Pollack, J.B., 1990. Mass fractionation of noble gases in diffusion-limited hydrodynamic hydrogen escape. Icarus 84, 502–527.