Hydrogen-dominated atmospheres of hot exoplanets expand and escape due to the intense heating by the X-ray and extreme ultraviolet (XUV) irradiation of their host stars. Excess absorption of neutral hydrogen has been observed in the Ly$\alpha$ line during transits of several close-in exoplanets, indicating such extended atmospheres. For the hot Jupiter HD 189733b, this absorption shows temporal variability. Variations in stellar XUV emission and/or stellar wind conditions have been invoked to explain this effect. We apply a 1D hydrodynamic upper atmosphere model and a 3D MHD stellar wind flow model to study the effect of variations of the stellar XUV and wind conditions on the neutral hydrogen distribution, including the production of energetic neutral atoms (ENAs), and the related Ly$\alpha$ transit signature. We obtain comparable, albeit slightly higher Ly$\alpha$ absorption as observed in 2011 with a stellar XUV flux of $1.8\times10^4$ erg cm$^{-2}$ s$^{-1}$, rather typical activity conditions for this star. Flares similar to the one observed 8 h before the transit are unlikely to have caused a significant modulation of the transit signature. The resulting Ly$\alpha$ absorption is dominated by atmospheric broadening, whereas the contribution of ENAs is negligible, as they are formed inside the bow shock from decelerated wind ions that are heated to high temperatures. Thus, within our modeling framework and assumptions, we find an insignificant dependence on the stellar wind parameters. Since the transit absorption can be modeled with typical stellar XUV and wind conditions, it is possible that the non-detection of the absorption in 2010 was affected by less typical stellar activity conditions, such as a very different magnitude and/or shape of the star's spectral XUV emission, or temporal/spatial variations in Ly$\alpha$ affecting the determination of the transit absorption.

In the present chapter we present the results of evolutionary studies of exoplanetary atmospheres. We mostly focus on the sub- to super-Earth domain, although these methods are applicable to all types of exoplanets. We consider both thermal and nonthermal loss processes. The type of thermal loss mechanism depends on so-called escape parameter $\beta$, which is the ratio of the gravitational energy of a particle to its thermal energy. While $\beta$ is decreasing, an exoplanet switches from classical Jeans to modified Jeans and finally to blow-off escape mechanisms. During blow-off the majority of the atmospheric particles dispose of enough energy to escape the planet's gravity field. This leads to extreme gas losses. Although nonthermal losses never exceed blow-off escape, they are of significant importance for planets with relatively weak Jeans-type escape. From the diversity of nonthermal escape mechanisms, in the present chapter we focus on ion pickup and discuss the importance of other loss mechanisms. The general conclusion of the chapter is, that escape processes strongly shape the evolution of the exoplanets and determine, if the planet loses its atmosphere due to erosion processes or, on the contrary, stays as mini-Neptune type body, which can probably not be considered as a potential habitat as we know it.

The influence of the hydrogen hydrodynamic upper atmosphere escape, driven by the solar soft X-ray and extreme ultraviolet radiation (XUV) flux, on an expected magma ocean outgassed steam atmosphere of early Venus is studied. By assuming that the young Sun was either a weak or moderate active young G star, we estimated the water loss from a hydrogen dominated thermosphere due to the absorption of the solar XUV flux and the precipitation of solar wind produced energetic hydrogen atoms (ENAs). The production of ENAs and their interaction with the hydrodynamic extended upper atmosphere, including collision-related feedback processes, have been calculated by means of Monte Carlo models. ENAs that collide in the upper atmosphere deposit their energy and heat the surrounding gas mainly above the main XUV energy deposition layer. It is shown that precipitating ENAs modify the thermal structure of the upper atmosphere, but the enhancement of the thermal escape rates caused by these energetic hydrogen atoms is negligible. Our results also indicate that the majority of oxygen arising from dissociated H$_2$O molecules is left behind during the first 100 Myr. It is thus suggested that the main part of the remaining oxygen has been absorbed by crustal oxidation.

The atmospheres of close-in planets are strongly influenced by mass loss driven by the high-energy (X-ray and extreme ultraviolet, EUV) irradiation of the host star, particularly during the early stages of evolution. We recently developed a framework to exploit this connection and enable us to recover the past evolution of the stellar high-energy emission from the present-day properties of its planets, if the latter retains some remnants of their primordial hydrogen-dominated atmospheres. Furthermore, the framework can also provide constraints on planetary initial atmospheric mass fractions. The constraints on the output parameters improve when more planets can be simultaneously analysed. This makes the Kepler-11 system, which hosts six planets with bulk densities between 0.66 and 2.45g cm^{-3}, an ideal target. Our results indicate that the star has likely evolved as a slow rotator (slower than 85\% of the stars with similar masses), corresponding to a high-energy emission at 150 Myr of between 1-10 times that of the current Sun. We also constrain the initial atmospheric mass fractions for the planets, obtaining a lower limit of 4.1% for planet c, a range of 3.7-5.3% for planet d, a range of 11.1-14% for planet e, a range of 1-15.6% for planet f, and a range of 4.7-8.7% for planet g assuming a disc dispersal time of 1 Myr. For planet b, the range remains poorly constrained. Our framework also suggests slightly higher masses for planets b, c, and f than have been suggested based on transit timing variation measurements. We coupled our results with published planet atmosphere accretion models to obtain a temperature (at 0.25 AU, the location of planet f) and dispersal time of the protoplanetary disc of 550 K and 1 Myr, although these results may be affected by inconsistencies in the adopted system parameters.

Planet atmospheric escape induced by high-energy stellar irradiation is a key phenomenon shaping the structure and evolution of planetary atmospheres. Therefore, the present-day properties of a planetary atmosphere are intimately connected with the amount of stellar flux received by a planet during its lifetime, thus with the evolutionary path of its host star. Using a recently developed analytic approximation based on hydrodynamic simulations for atmospheric escape rates, we track within a Bayesian framework the evolution of a planet as a function of stellar flux evolution history, constrained by the measured planetary radius, with the other system parameters as priors. We find that the ideal objects for this type of study are close-in sub-Neptune-like planets, as they are highly affected by atmospheric escape, and yet retain a significant fraction of their primordial hydrogen-dominated atmospheres. Furthermore, we apply this analysis to the HD3167 and K2-32 planetary systems. For HD3167, we find that the most probable irradiation level at 150 Myr was between 40 and 130 times solar, corresponding to a rotation period of 1.78^{+2.69}_{-1.23} days. For K2-32, we find a surprisingly low irradiation level ranging between half and four times solar at 150 Myr. Finally, we show that for multi-planet systems, our framework enables one to constrain poorly known properties of individual planets.

We modeled the transit signatures in the Lya line of a putative Earth-sized planet orbiting in the HZ of the M dwarf GJ436. We estimated the transit depth in the Lya line for an exo-Earth with three types of atmospheres: a hydrogen-dominated atmosphere, a nitrogen-dominated atmosphere, and a nitrogen-dominated atmosphere with an amount of hydrogen equal to that of the Earth. We calculated the in-transit absorption they would produce in the Lya line. We applied it to the out-of-transit Lya observations of GJ 436 obtained by the HST and compared the calculated in-transit absorption with observational uncertainties to determine if it would be detectable. To validate the model, we also used our method to simulate the deep absorption signature observed during the transit of GJ 436b and showed that our model is capable of reproducing the observations. We used a DSMC code to model the planetary exospheres. The code includes several species and traces neutral particles and ions. At the lower boundary of the DSMC model we assumed an atmosphere density, temperature, and velocity obtained with a hydrodynamic model for the lower atmosphere. We showed that for a small rocky Earth-like planet orbiting in the HZ of GJ436 only the hydrogen-dominated atmosphere is marginally detectable with the STIS/HST. Neither a pure nitrogen atmosphere nor a nitrogen-dominated atmosphere with an Earth-like hydrogen concentration in the upper atmosphere are detectable. We also showed that the Lya observations of GJ436b can be reproduced reasonably well assuming a hydrogen-dominated atmosphere, both in the blue and red wings of the Lya line, which indicates that warm Neptune-like planets are a suitable target for Lya observations. Terrestrial planets can be observed in the Lya line if they orbit very nearby stars, or if several observational visits are available.

Studies of planetary atmospheric composition, variability, and evolution require appropriate theoretical and numerical tools to estimate key atmospheric parameters, among which the mass-loss rate is often the most important. In evolutionary studies, it is common to use the energy-limited formula, which is attractive for its simplicity but ignores important physical effects and can be inaccurate in many cases. To overcome this problem, we consider a recently developed grid of about 7000 one-dimensional upper-atmosphere hydrodynamic models computed for a wide range of planets with hydrogen-dominated atmospheres from which we extract the mass-loss rates. The grid boundaries are [1:39] MEARTH in planetary mass, [1:10] REARTH in planetary radius, [300:2000] K in equilibrium temperature, [0.4:1.3] MSUN in host star's mass, [0.002:1.3] au in orbital separation, and about [10^{26}:5*10^{30}] erg/s in stellar X-ray and extreme ultraviolet luminosity. We then derive an analytical expression for the atmospheric mass-loss rates based on a fit to the values obtained from the grid. The expression provides the mass-loss rates as a function of planetary mass, planetary radius, orbital separation, and incident stellar high-energy flux. We show that this expression is a significant improvement to the energy-limited approximation for a wide range of planets. The analytical expression presented here enables significantly more accurate planetary evolution computations without increasing computing time.

There is growing observational and theoretical evidence suggesting that atmospheric escape is a key driver of planetary evolution. Commonly, planetary evolution models employ simple analytic formulae (e.g., energy limited escape) that are often inaccurate, and more detailed physical models of atmospheric loss usually only give snapshots of an atmosphere's structure and are difficult to use for evolutionary studies. To overcome this problem, we upgrade and employ an already existing upper atmosphere hydrodynamic code to produce a large grid of about 7000 models covering planets with masses 1 - 39 Earth mass with hydrogen-dominated atmospheres and orbiting late-type stars. The modeled planets have equilibrium temperatures ranging between 300 and 2000 K. For each considered stellar mass, we account for three different values of the high-energy stellar flux (i.e., low, moderate, and high activity). For each computed model, we derive the atmospheric temperature, number density, bulk velocity, X-ray and EUV (XUV) volume heating rates, and abundance of the considered species as a function of distance from the planetary center. From these quantities, we estimate the positions of the maximum dissociation and ionisation, the mass-loss rate, and the effective radius of the XUV absorption. We show that our results are in good agreement with previously published studies employing similar codes. We further present an interpolation routine capable to extract the modelling output parameters for any planet lying within the grid boundaries. We use the grid to identify the connection between the system parameters and the resulting atmospheric properties. We finally apply the grid and the interpolation routine to estimate atmospheric evolutionary tracks for the close-in, high-density planets CoRoT-7 b and HD219134 b,c...

We investigate under which conditions supermassive hot Jupiters can sustain source regions for radio emission, and whether this emission could propagate to an observer outside the system. We study Tau Bootis b-like planets (a supermassive hot Jupiter with 5.84 Jupiter masses and 1.06 Jupiter radii), but located at different orbital distances (between its actual orbit of 0.046 AU and 0.2 AU). Due to the strong gravity of such planets and efficient radiative cooling, the upper atmosphere is (almost) hydrostatic and the exobase remains very close to the planet, which makes it a good candidate for radio observations. We expect similar conditions as for Jupiter, i.e. a region between the exobase and the magnetopause that is filled with a depleted plasma density compared with cases where the whole magnetosphere cavity is filled with hydrodynamically outward flowing ionospheric plasma. Thus, unlike classical hot Jupiters like the previously studied planets HD 209458b and HD 189733b, supermassive hot Jupiters should be in general better targets for radio observations.

The K2-33 planetary system hosts one transiting ~5 R_E planet orbiting the young M-type host star. The planet's mass is still unknown, with an estimated upper limit of 5.4 M_J. The extreme youth of the system (<20 Myr) gives the unprecedented opportunity to study the earliest phases of planetary evolution, at a stage when the planet is exposed to an extremely high level of high-energy radiation emitted by the host star. We perform a series of 1D hydrodynamic simulations of the planet's upper atmosphere considering a range of possible planetary masses, from 2 to 40 M_E, and equilibrium temperatures, from 850 to 1300 K, to account for internal heating as a result of contraction. We obtain temperature profiles mostly controlled by the planet's mass, while the equilibrium temperature has a secondary effect. For planetary masses below 7-10 M_E, the atmosphere is subject to extremely high escape rates, driven by the planet's weak gravity and high thermal energy, which increase with decreasing mass and/or increasing temperature. For higher masses, the escape is instead driven by the absorption of the high-energy stellar radiation. A rough comparison of the timescales for complete atmospheric escape and age of the system indicates that the planet is more massive than 10 M_E.

Planetary embryos form protoplanets via mutual collisions, which can lead to the development of magma oceans. During their solidification, large amounts of the mantles' volatile contents may be outgassed. The resulting H$_2$O/CO$_2$ dominated steam atmospheres may be lost efficiently via hydrodynamic escape due to the low gravity and the high stellar EUV luminosities. Protoplanets forming later from such degassed building blocks could therefore be drier than previously expected. We model the outgassing and subsequent hydrodynamic escape of steam atmospheres from such embryos. The efficient outflow of H drags along heavier species (O, CO$_2$, noble gases). The full range of possible EUV evolution tracks of a solar-mass star is taken into account to investigate the escape from Mars-sized embryos at different orbital distances. The envelopes are typically lost within a few to a few tens of Myr. Furthermore, we study the influence on protoplanetary evolution, exemplified by Venus. We investigate different early evolution scenarios and constrain realistic cases by comparing modeled noble gas isotope ratios with observations. Starting from solar values, consistent isotope ratios (Ne, Ar) can be found for different solar EUV histories, as well as assumptions about the initial atmosphere (either pure steam or a mixture with accreted H). Our results generally favor an early accretion scenario with a small amount of accreted H and a low-activity Sun, because in other cases too much CO$_2$ is lost during evolution, which is inconsistent with Venus' present atmosphere. Important issues are likely the time at which the initial steam atmosphere is outgassed and/or the amount of CO$_2$ which may still be delivered at later evolutionary stages. A late accretion scenario can only reproduce present isotope ratios for a highly active young Sun, but then very massive steam atmospheres would be required.

The strong, nearly wavelength-independent absorption cross section of aerosols produces featureless exoplanet transmission spectra, limiting our ability to characterize their atmospheres. Here we show that even in the presence of featureless spectra, we can still characterize certain atmospheric properties. Specifically, we constrain the upper and lower pressure boundaries of aerosol layers, and present plausible composition candidates. We study the case of the bloated Saturn-mass planet WASP-49b, where near-infrared observations reveal a flat transmission spectrum between 0.7 and 1.0 {\microns}. First, we use a hydrodynamic upper-atmosphere code to estimate the pressure reached by the ionizing stellar high-energy photons at $10^{-8}$ bar, setting the upper pressure boundary where aerosols could exist. Then, we combine HELIOS and Pyrat Bay radiative-transfer models to constrain the temperature and photospheric pressure of atmospheric aerosols, in a Bayesian framework. For WASP-49b, we constrain the transmission photosphere (hence, the aerosol deck boundaries) to pressures above $10^{-5}$ bar (100$\times$ solar metallicity), $10^{-4}$ bar (solar), and $10^{-3}$ bar (0.1$\times$ solar) as lower boundary, and below $10^{-7}$ bar as upper boundary. Lastly, we compare condensation curves of aerosol compounds with the planet's pressure-temperature profile to identify plausible condensates responsible for the absorption. Under these circumstances, we find as candidates: Na$_{2}$S (at 100$\times$ solar metallicity); Cr and MnS (at solar and 0.1$\times$ solar); and forsterite, enstatite, and alabandite (at 0.1$\times$ solar).

We investigate the interaction between the magnetized stellar wind plasma and the partially ionized hydrodynamic hydrogen outflow from the escaping upper atmosphere of non- or weakly magnetized hot Jupiters. We use the well-studied hot Jupiter HD 209458b as an example for similar exoplanets, assuming a negligible intrinsic magnetic moment. For this planet, the stellar wind plasma interaction forms an obstacle in the planet's upper atmosphere, in which the position of the magnetopause is determined by the condition of pressure balance between the stellar wind and the expanded atmosphere, heated by the stellar extreme ultraviolet (EUV) radiation. We show that the neutral atmospheric atoms penetrate into the region dominated by the stellar wind, where they are ionized by photo-ionization and charge exchange, and then mixed with the stellar wind flow. Using a 3D magnetohydrodynamic (MHD) model, we show that an induced magnetic field forms in front of the planetary obstacle, which appears to be much stronger compared to those produced by the solar wind interaction with Venus and Mars. Depending on the stellar wind parameters, because of the induced magnetic field, the planetary obstacle can move up to ~0.5-1 planetary radii closer to the planet. Finally, we discuss how estimations of the intrinsic magnetic moment of hot Jupiters can be inferred by coupling hydrodynamic upper planetary atmosphere and MHD stellar wind interaction models together with UV observations. In particular, we find that HD 209458b should likely have an intrinsic magnetic moment of 10-20% that of Jupiter.

Stimulated by the discovery of a number of close-in low-density planets, we generalise the Jeans escape parameter taking hydrodynamic and Roche lobe effects into account. We furthermore define $\Lambda$ as the value of the Jeans escape parameter calculated at the observed planetary radius and mass for the planet's equilibrium temperature and considering atomic hydrogen, independently of the atmospheric temperature profile. We consider 5 and 10 $M_{\oplus}$ planets with an equilibrium temperature of 500 and 1000 K, orbiting early G-, K-, and M-type stars. Assuming a clear atmosphere and by comparing escape rates obtained from the energy-limited formula, which only accounts for the heating induced by the absorption of the high-energy stellar radiation, and from a hydrodynamic atmosphere code, which also accounts for the bolometric heating, we find that planets whose $\Lambda$ is smaller than 15-35 lie in the "boil-off" regime, where the escape is driven by the atmospheric thermal energy and low planetary gravity. We find that the atmosphere of hot (i.e. $T_{\rm eq}\gtrapprox$ 1000 K) low-mass ($M_{\rm pl}\lessapprox$ 5 $M_{\oplus}$) planets with $\Lambda$ < 15-35 shrinks to smaller radii so that their $\Lambda$ evolves to values higher than 15-35, hence out of the boil-off regime, in less than $\approx$500 Myr. Because of their small Roche lobe radius, we find the same result also for hot (i.e. $T_{\rm eq}\gtrapprox$ 1000 K) higher mass ($M_{\rm pl}\lessapprox$ 10 $M_{\oplus}$) planets with $\Lambda$ < 15-35, when they orbit M-dwarfs. For old, hydrogen-dominated planets in this range of parameters, $\Lambda$ should therefore be $\geq$15-35, which provides a strong constraint on the planetary minimum mass and maximum radius and can be used to predict the presence of aerosols and/or constrain planetary masses, for example.

We present a uniform analysis of the atmospheric escape rate of Neptune-like planets with estimated radius and mass (restricted to $M_{\rm p}<30\,M_{\oplus}$). For each planet we compute the restricted Jeans escape parameter, $\Lambda$, for a hydrogen atom evaluated at the planetary mass, radius, and equilibrium temperature. Values of $\Lambda\lesssim20$ suggest extremely high mass-loss rates. We identify 27 planets (out of 167) that are simultaneously consistent with hydrogen-dominated atmospheres and are expected to exhibit extreme mass-loss rates. We further estimate the mass-loss rates ($L_{\rm hy}$) of these planets with tailored atmospheric hydrodynamic models. We compare $L_{\rm hy}$ to the energy-limited (maximum-possible high-energy driven) mass-loss rates. We confirm that 25 planets (15\% of the sample) exhibit extremely high mass-loss rates ($L_{\rm hy}>0.1\,M_{\oplus}{\rm Gyr}^{-1}$), well in excess of the energy-limited mass-loss rates. This constitutes a contradiction, since the hydrogen envelopes cannot be retained given the high mass-loss rates. We hypothesize that these planets are not truly under such high mass-loss rates. Instead, either hydrodynamic models overestimate the mass-loss rates, transit-timing-variation measurements underestimate the planetary masses, optical transit observations overestimate the planetary radii (due to high-altitude clouds), or Neptunes have consistently higher albedos than Jupiter planets. We conclude that at least one of these established estimations/techniques is consistently producing biased values for Neptune planets. Such an important fraction of exoplanets with misinterpreted parameters can significantly bias our view of populations studies, like the observed mass--radius distribution of exoplanets for example.

We suggest an additional possible plasma source to explain part of the phenomena observed for the transiting hot Jupiters WASP-12b and HD 189733b in their ultraviolet (UV) light curves. In the proposed scenario, material outgasses from the molten surface of Trojan satellites on tadpole orbits near the Lagrange points L$_4$ and L$_5$. We show that the temperature at the orbital location of WASP-12b is high enough to melt the surface of rocky bodies and to form shallow lava oceans on them. In case of WASP-12b, this leads to the release of elements such as Mg and Ca, which are expected to surround the system. The predicted Mg and Ca outgassing rates from two Io-sized WASP-12b Trojans are $\approx 2.2 \times 10^{27}$ s$^{-1}$ and $\approx 2.2 \times 10^{26}$ s$^{-1}$, respectively. Trojan outgassing can lead to the apparent lack of emission in Mg{\sc ii}\,h\&k and Ca{\sc ii}\,H\&K line cores of WASP-12. For HD 189733b, the mechanism is only marginally possible due to the lower temperature. This may be one of the reasons that couldn't explain the early ingress of HD 189733b observed in the far-UV (FUV) C{\sc ii} doublet due to absence of carbon within elements outgassed by molten lava. We investigate the long-term stability region of WASP-12b and HD 189733b in case of planar and inclined motion of these satellites and show that unlike the classical exomoons orbiting the planet, Io-sized Trojans can be stable for the whole systems life time.

For the hot exoplanets CoRoT-24b and CoRoT-24c, observations have provided transit radii R$_{\rm T}$ of 3.7$\pm$0.4 R$_{\oplus}$ and 4.9$\pm$0.5 R$_{\oplus}$, and masses of $\le$5.7 M$_{\oplus}$ and 28$\pm$11 M$_{\oplus}$, respectively. We study their upper atmosphere structure and escape applying an hydrodynamic model. Assuming R$_{\rm T} \approx$ R$_{\rm PL}$, where R$_{\rm PL}$ is the planetary radius at the pressure of 100 mbar, we obtained for CoRoT-24b unrealistically high thermally-driven hydrodynamic escape rates. This is due to the planet's high temperature and low gravity, independent of the stellar EUV flux. Such high escape rates could last only for $<$100 Myr, while R$_{\rm PL}$ shrinks till the escape rate becomes less than or equal to the maximum possible EUV-driven escape rate. For CoRoT-24b, R$_{\rm PL}$ must be therefore located at $\approx 1.9-2.2$ R$_{\oplus}$ and high altitude hazes/clouds possibly extinct the light at R$_{\rm T}$. Our analysis constraints also the planet's mass to be 5$-$5.7 M$_{\oplus}$. For CoRoT-24c, R$_{\rm PL}$ and R$_{\rm T}$ lie too close together to be distinguished in the same way. Similar differences between R$_{\rm PL}$ and R$_{\rm T}$ may be present also for other hot, low-density sub-Neptunes.

We investigate the loss rates of the hydrogen atmospheres of terrestrial planets with a range of masses and orbital distances by assuming a stellar extreme ultraviolet (EUV) luminosity that is 100 times stronger than that of the current Sun. We apply a 1D upper atmosphere radiation absorption and hydrodynamic escape model that takes into account ionization, dissociation and recombination to calculate hydrogen mass loss rates. We study the effects of the ionization, dissociation and recombination on the thermal mass loss rates of hydrogen-dominated super-Earths and compare the results to those obtained by the energy-limited escape formula which is widely used for mass loss evolution studies. Our results indicate that the energy-limited formula can to a great extent over- or underestimate the hydrogen mass loss rates by amounts that depend on the stellar EUV flux and planetary parameters such as mass, size, effective temperature, and EUV absorption radius.

Terrestrial planets formed within gaseous protoplanetary disks can accumulate significant hydrogen envelopes. The evolution of such an atmosphere due to XUV driven evaporation depends on the activity evolution of the host star, which itself depends sensitively on its rotational evolution, and therefore on its initial rotation rate. In this letter, we derive an easily applicable method for calculating planetary atmosphere evaporation that combines models for a hydrostatic lower atmosphere and a hydrodynamic upper atmosphere. We show that the initial rotation rate of the central star is of critical importance for the evolution of planetary atmospheres and can determine if a planet keeps or loses its primordial hydrogen envelope. Our results highlight the need for a detailed treatment of stellar activity evolution when studying the evolution of planetary atmospheres.

We study the origin and escape of catastrophically outgassed volatiles (H$_2$O, CO$_2$) from exomoons with Earth-like densities and masses of $0.1M_{\oplus}$, $0.5M_{\oplus}$ and $1M_{\oplus}$ orbiting an extra-solar gas giant inside the habitable zone of a young active solar-like star. We apply a radiation absorption and hydrodynamic upper atmosphere model to the three studied exomoon cases. We model the escape of hydrogen and dragged dissociation products O and C during the activity saturation phase of the young host star. Because the soft X-ray and EUV radiation of the young host star may be up to $\sim$100 times higher compared to today's solar value during the first 100 Myr after the system's origin, an exomoon with a mass $ < 0.25M_{\oplus}$ located in the HZ may not be able to keep an atmosphere because of its low gravity. Depending on the spectral type and XUV activity evolution of the host star, exomoons with masses between $\sim0.25-0.5M_{\oplus}$ may evolve to Mars-like habitats. More massive bodies with masses $ > 0.5M_{\oplus}$, however, may evolve to habitats that are a mixture of Mars-like and Earth-analogue habitats, so that life may originate and evolve at the exomoon's surface.

By considering martian-like planetary embryos inside the habitable zone of solar-like stars we study the behavior of the hydrodynamic atmospheric escape of hydrogen for small values of the Jeans escape parameter $\beta < 3$, near the base of the thermosphere, that is defined as a ratio of the gravitational and thermal energy. Our study is based on a 1-D hydrodynamic upper atmosphere model that calculates the volume heating rate in a hydrogen dominated thermosphere due to the absorption of the stellar soft X-ray and extreme ultraviolet (XUV) flux. We find that when the $\beta$ value near the mesopause/homopause level exceeds a critical value of $\sim$2.5, there exists a steady hydrodynamic solution with a smooth transition from subsonic to supersonic flow. For a fixed XUV flux, the escape rate of the upper atmosphere is an increasing function of the temperature at the lower boundary. Our model results indicate a crucial enhancement of the atmospheric escape rate, when the Jeans escape parameter $\beta$ decreases to this critical value. When $\beta$ becomes $\leq$2.5, there is no stationary hydrodynamic transition from subsonic to supersonic flow. This is the case of a fast non-stationary atmospheric expansion that results in extreme thermal atmospheric escape rates.

In the present series of papers we propose a consistent description of the mass loss process. To study the effects of intrinsic magnetic field of a close-orbit giant exoplanet (so-called Hot Jupiter) on the atmospheric material escape and formation of planetary inner magnetosphere in a comprehensive way, we start with a hydrodynamic model of an upper atmosphere expansion presented in this paper. While considering a simple hydrogen atmosphere model, we focus on selfconsistent inclusion of the effects of radiative heating and ionization of the atmospheric gas with its consequent expansion in the outer space. Primary attention is paid to investigation of the role of specific conditions at the inner and outer boundaries of the simulation domain, under which different regimes of material escape (free- and restricted- flow) are formed. Comparative study of different processes, such as XUV heating, material ionization and recombination, H3+ cooling, adiabatic and Lyman-alpha cooling, Lyman-alpha reabsorption is performed. We confirm basic consistence of the outcomes of our modeling with the results of other hydrodynamic models of expanding planetary atmospheres. In particular, we obtain that under the typical conditions of an orbital distance 0.05 AU around a Sun-type star a Hot Jupiter plasma envelope may reach maximum temperatures up to ~9000K with a hydrodynamic escape speed ~9 km/s resulting in the mass loss rates ~(4-7)*10^10 g*s . In the range of considered stellar-planetary parameters and XUV fluxes that is close to mass loss in the energy limited case. The inclusion of planetary intrinsic magnetic fields in the model is a subject of the following up paper (Paper II).

With the discovery of hundreds of exoplanets and a potentially huge number of Earth-like planets waiting to be discovered, the conditions for their habitability have become a focal point in exoplanetary research. The classical picture of habitable zones primarily relies on the stellar flux allowing liquid water to exist on the surface of an Earth-like planet with a suitable atmosphere. However, numerous further stellar and planetary properties constrain habitability. Apart from "geophysical" processes depending on the internal structure and composition of a planet, a complex array of astrophysical factors additionally determine habitability. Among these, variable stellar UV, EUV, and X-ray radiation, stellar and interplanetary magnetic fields, ionized winds, and energetic particles control the constitution of upper planetary atmospheres and their physical and chemical evolution. Short- and long-term stellar variability necessitates full time-dependent studies to understand planetary habitability at any point in time. Furthermore, dynamical effects in planetary systems and transport of water to Earth-like planets set fundamentally important constraints. We will review these astrophysical conditions for habitability under the crucial aspects of the long-term evolution of stellar properties, the consequent extreme conditions in the early evolutionary phase of planetary systems, and the important interplay between properties of the host star and its planets.

We investigate the influence of impacts of large planetesimals and small planetary embryos on the early Martian surface on the hydrodynamic escape of an early steam atmosphere that is exposed to the high soft X-ray and EUV flux of the young Sun. Impact statistics in terms of number, masses, velocities, and angles of asteroid impacts onto the early Mars are determined via n-body integrations. Based on these statistics, smoothed particle hydrodynamics (SPH) simulations result in estimates of energy transfer into the planetary surface material and according surface heating. For the estimation of the atmospheric escape rates we applied a soft X-ray and EUV absorption model and a 1-D upper atmosphere hydrodynamic model to a magma ocean-related catastrophically outgassed steam atmosphere with surface pressure values of 52 bar H2O and 11 bar CO2. The estimated impact rates and energy deposition onto an early Martian surface can account for substantial heating. The energy influx and conversion rate into internal energy is most likely sufficient to keep a shallow magma ocean liquid for an extended period of time. Higher surface temperatures keep the outgassed steam atmosphere longer in vapor form and therefore enhance its escape to space within about 0.6 Myr after its formation.

We study the interactions between stellar wind and the extended hydrogen-dominated upper atmospheres of planets and the resulting escape of planetary pick-up ions from the 5 "super-Earths" in the compact Kepler-11 system and compare the escape rates with the efficiency of the thermal escape of neutral hydrogen atoms. Assuming the stellar wind of Kepler-11 is similar to the solar wind, we use a polytropic 1D hydrodynamic wind model to estimate the wind properties at the planetary orbits. We apply a Direct Simulation Monte Carlo Model to model the hydrogen coronae and the stellar wind plasma interaction around Kepler-11b-f within a realistic expected heating efficiency range of 15-40%. The same model is used to estimate the ion pick-up escape from the XUV heated and hydrodynamically extended upper atmospheres of Kepler-11b-f. From the interaction model we study the influence of possible magnetic moments, calculate the charge exchange and photoionization production rates of planetary ions and estimate the loss rates of pick-up H+ ions for all five planets. We compare the results between the five "super-Earths" and in a more general sense also with the thermal escape rates of the neutral planetary hydrogen atoms. Our results show that for all Kepler-11b-f exoplanets, a huge neutral hydrogen corona is formed around the planet. The non-symmetric form of the corona changes from planet to planet and is defined mostly by radiation pressure and gravitational effects. Non-thermal escape rates of pick-up ionized hydrogen atoms for Kepler-11 "super-Earths" vary between approximately 6.4e30 1/s and 4.1e31 1/s depending on the planet's orbital location and assumed heating efficiency. These values correspond to non-thermal mass loss rates of approximately 1.07e7 g/s and 6.8e7 g/s respectively, which is a few percent of the thermal escape rates.

We investigate the origin and loss of captured hydrogen envelopes from protoplanets between `sub-Earth'-like bodies of 0.1$M_{\oplus}$ up to `super-Earths' with 5$M_{\oplus}$ in the HZ of a Sun like G star, assuming their rocky cores had formed before the nebula dissipated. We model the gravitational accumulation of nebula gas around a core as a function of protoplanetary luminosity during accretion and calculate the resulting surface temperature by solving the hydrostatic structure equations for the protoplanetary nebula. Depending on nebular properties and resulting luminosities, for planetary bodies of 0.1--1$M_{\oplus}$ we obtain hydrogen envelopes with masses between $\sim 2.5\times 10^{19}$--$1.5\times 10^{26}$ g. For `super-Earths' with masses between 2--5$M_{\oplus}$ hydrogen envelopes within the mass range of $\sim 7.5\times 10^{23}$--$1.5\times 10^{28}$ g can be captured. To study the escape of these hydrogen-dominated protoatmospheres, we apply a hydrodynamic upper atmosphere model and calculate the loss rates due to the heating by the high XUV flux of the young star. Our results indicate that under most nebula conditions `sub-Earth' and Earth-mass planets can lose their envelopes by thermal escape during the first $100$ Myr after the disk dissipated. However, if a nebula has a low dust depletion factor or low accretion rates resulting in low protoplanetary luminosities, it is possible that even protoplanets with Earth-mass cores may keep their hydrogen envelopes during their whole lifetime. In contrast to lower mass protoplanets, `super-Earths' accumulate a huge amount of nebula gas and lose only tiny fractions of their primordial envelopes. Our results agree with the fact that Venus, Earth, and Mars are not surrounded by dense hydrogen envelopes, as well as with the recent discoveries of low density `super-Earths' that most likely could not get rid of their protoatmospheres.

Latest research in planet formation indicate that Mars formed within a few million years (Myr) and remained 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 H2O. 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 >=3x10^{19} g to <=6.5x10^{22} g could have been captured from the nebula around early Mars. Depending of 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 H2O and ~10-55 bar CO2 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... (continued)

The recently discovered low-density "super-Earths" Kepler-11b, Kepler-11f, Kepler-11d, Kepler-11e, and planets such as GJ 1214b represent most likely planets which are surrounded by dense H/He envelopes or contain deep H2O oceans also surrounded by dense hydrogen envelopes. Although these "super-Earths" are orbiting relatively close to their host stars, they have not lost their captured nebula-based hydrogen-rich or degassed volatile-rich steam protoatmospheres. Thus it is interesting to estimate the maximum possible amount of atmospheric hydrogen loss from a terrestrial planet orbiting within the habitable zone of late main sequence host stars. For studying the thermosphere structure and escape we apply a 1-D hydrodynamic upper atmosphere model which solves the equations of mass, momentum and energy conservation for a planet with the mass and size of the Earth and for a "super-Earth" with a size of 2 R_Earth and a mass of 10 M_Earth. We calculate volume heating rates by the stellar soft X-ray and EUV radiation and expansion of the upper atmosphere, its temperature, density and velocity structure and related thermal escape rates during planet's life time. Moreover, we investigate under which conditions both planets enter the blow-off escape regime and may therefore experience loss rates which are close to the energy-limited escape. Finally we discuss the results in the context of atmospheric evolution and implications for habitability of terrestrial planets in general.

We study the interactions between the stellar wind plasma flow of a typical M star, such as GJ 436, and hydrogen-rich upper atmospheres of an Earth-like planet and a "super-Earth" with the radius of 2 R_Earth and a mass of 10 M_Earth, located within the habitable zone at ~0.24 AU. We investigate the formation of extended atomic hydrogen coronae under the influences of the stellar XUV flux (soft X-rays and EUV), stellar wind density and velocity, shape of a planetary obstacle (e.g., magnetosphere, ionopause), and the loss of planetary pick-up ions on the evolution of hydrogen-dominated upper atmospheres. Stellar XUV fluxes which are 1, 10, 50 and 100 times higher compared to that of the present-day Sun are considered and the formation of high-energy neutral hydrogen clouds around the planets due to the charge-exchange reaction under various stellar conditions have been modeled. Charge-exchange between stellar wind protons with planetary hydrogen atoms, and photoionization, leads to the production of initially cold ions of planetary origin. We found that the ion production rates for the studied planets can vary over a wide range, from ~1.0x10^{25} s^{-1} to ~5.3x10^{30} s^{-1}, depending on the stellar wind conditions and the assumed XUV exposure of the upper atmosphere. Our findings indicate that most likely the majority of these planetary ions are picked up by the stellar wind and lost from the planet. Finally, we estimate the long-time non-thermal ion pick-up escape for the studied planets and compare them with the thermal escape. According to our estimates, non-thermal escape of picked up ionized hydrogen atoms over a planet's lifetime varies between ~0.4 Earth ocean equivalent amounts of hydrogen (EO_H) to <3 EO_H and usually is several times smaller in comparison to the thermal atmospheric escape rates.

The discovery of transiting "super-Earths" with inflated radii and known masses such as Kepler-11b-f, GJ 1214b and 55 Cnc e, indicates that these exoplanets did not lose their nebula-captured hydrogen-rich, degassed or impact-delivered protoatmospheres by atmospheric escape processes. Because hydrodynamic blow-off of atmospheric hydrogen atoms is the most efficient atmospheric escape process we apply a time-dependent numerical algorithm which is able to solve the system of 1-D fluid equations for mass, momentum, and energy conservation to investigate the criteria under which "super-Earths" with hydrogen-dominated upper atmospheres can experience hydrodynamic expansion by heating of the stellar XUV (soft X-rays and extreme ultraviolet) radiation and thermal escape via blow-off. Depending on orbit location, XUV flux, heating efficiency and the planet's mean density our results indicate that the upper atmospheres of all "super-Earths" can expand to large distances, so that besides of Kepler-11c all of them experience atmospheric mass-loss due to Roche lobe overflow. The atmospheric mass-loss of the studied "super-Earths" is one to two orders of magnitude lower compared to that of "hot Jupiters" such as HD 209458b, so that one can expect that these exoplanets cannot lose their hydrogen-envelopes during their remaining lifetimes.

Observational evidence of a hydrodynamically evaporating upper atmosphere of HD209458b (Vidal-Madjar et al. 2003; 2004) and recent theoretical studies on evaporation scenarios of ``Hot Jupiters'' in orbits around solar-like stars with the age of the Sun indicate that the upper atmospheres of short-periodic exoplanets experience hydrodynamic blow-off conditions resulting in loss rates of the order of about 10^10 - 10^12 g s^-1 (Lammer et al. 2003; Yelle 2004; Baraffe et al. 2004; Lecavlier des Etangs et al. 2004; Jaritz et al. 2005, Tian et al. 2005; Penz et al. 2007). By studying the effect of the Roche lobe on the atmospheric loss from short-periodic gas giants we found, that the effect of the Roche lobe can enhance the hydrodynamic evaporation from HD209458b by about 2 and from OGLE-TR-56b by about 2.5 times. For similar exoplanets which are closer to their host star than OGLE-TR-56b, the enhancement of the mass loss can be even larger. Moreover, we show that the effect of the Roche lobe raises the possibility that ``Hot Jupiters'' can reach blow-off conditions at temperatures which are less than expected (< 10000 K) due to the stellar X-ray and EUV (XUV) heating.