Atmospheric Evolution on Low-gravity Waterworlds

TL;DR

This study investigates the atmospheric evolution of low-gravity waterworlds, showing they can maintain habitable conditions without runaway greenhouse effects, with habitability limited by atmospheric escape rates.

Contribution

The paper introduces a unified model combining energy balance, hydrodynamic escape, and radiative transfer to analyze habitability of low-gravity waterworlds, highlighting the physics of their atmospheric stability.

Findings

Low-gravity waterworlds do not undergo runaway greenhouse effects.

A long-lived habitable zone is possible for sufficiently massive low-gravity waterworlds.

Atmospheric escape rate limits the inner edge of habitability for these worlds.

Abstract

Low-gravity waterworlds ($M\lesssim 0.1 M_{\oplus}$) are of interest for their potential habitability. The weakly bound atmospheres of such worlds have proportionally larger radiative surfaces and are more susceptible to escape. We conduct a unified investigation into these phenomena, combining analytical energy balance and hydrodynamic escape with line-by-line radiative transfer calculations. Because outgoing radiation is forced to increase with surface temperature by the expansion of the radiative surface, we find that these worlds do not experience a runaway greenhouse. Furthermore, we show that a long-lived liquid water habitable zone is possible for low-gravity waterworlds of sufficient mass. Its inner edge is set by the rate of atmospheric escape, because a short-lived atmosphere limits the time available for life to evolve. In describing the physics of the parameter space transition from "planet-like" to "comet-like", our model produces a lower bound for habitability in terms of gravity. These results provide valuable insights in the context of the ongoing hunt for habitable exoplanets and exomoons.

Equations

$$\frac{R}{R_{\oplus }}=1.258{\left(\frac{M}{M_{\oplus }}\right)}^{0.302}$$

$${\left(\frac{r_{SW}}{r_s}\right)}^2\frac{1}{4}(1-A)F_{\text{stel}}-(g_s{r}_s+L)\Phi ={\left(\frac{r_{LW}}{r_s}\right)}^2{F}_{\text{out}}$$

$$\frac{dp}{dz}=-\rho g_s{\left(\frac{r_s}{r_s+z}\right)}^2$$

$$z_{\text{radiative surface}}=\frac{\int I(\nu )\times z_{\tau =1}(\nu )d\nu }{\int I(\nu )d\nu }$$

$$\frac{d}{dr}(r^2\rho w)=0$$

$$w\frac{dw}{dr}=-\frac{1}{\rho }\frac{dp}{dr}-\frac{GM}{r^2}$$

$$(w^2-w_c^2)\frac{1}{w}\frac{dw}{dr}=\frac{2w_c^2}{r}-\frac{GM}{r^2}$$

$$r_c=\frac{GM}{2w_c^2}$$

$$\Phi ={\rho }_s{w}_s={\rho }_s{w}_c{\left(\frac{r_c}{r_s}\right)}^2\text{exp}\left(-\frac{1}{2}+\frac{GM}{w_c^2}\left(\frac{1}{r_c}-\frac{1}{r_s}\right)\right)$$

$$\dot{M}=4\pi r_s^2\Phi$$

$${\tau }_{water}=\frac{M_{water}}{\dot{M}}$$

$$F_{\text{stel}}=\frac{(\frac{r_{LW}}{r_s}{)}^2{F}_{\text{out}}+(g_s{r}_s+L)\Phi }{(\frac{r_{SW}}{r_s}{)}^2\frac{1}{4}(1-A)}$$

$$A(T_s)=\{\begin{array}{ll}{\alpha }_I& T_s<273\text{ K}\\ {\alpha }_L& T_s\ge 273\text{ K}\end{array}$$

$$\rho w\frac{dw}{dr}+\frac{dp}{dz}=-\rho g$$

$$|w\frac{dw}{dr}|/|g(r)|$$

$$\ln (\frac{w}{w_c})-\frac{1}{2}(\frac{w}{w_c}{)}^2=-\frac{2r_c}{r}-2\ln (\frac{r}{r_c})+\frac{3}{2}$$

$$w(r)=w_c(\frac{r_c}{r}{)}^2\text{exp}(\frac{3}{2}-2\frac{r_c}{r})$$

$$w\frac{dw}{dr}=2w_c^2(\frac{r_c^5}{r^6}-\frac{r_c^4}{r^5})\text{exp}(3-4\frac{r_c}{r})$$

Full text

astrobiology – minor planets, asteroids: general – planets and satellites: atmospheres – planets and satellites: physical evolution

1 Introduction

Earth is the classic example of a terrestrial-mass, water-rich body with both an atmosphere and surface liquid H2O. Low-mass water-rich bodies (i.e. comets) have neither: the gravity is not strong enough to bind an atmosphere, and when they are heated their water ice is sublimated without any significant amounts of liquid water appearing on their surfaces. What happens at intermediate masses?

The evolution of low-gravity waterworlds differs from that of larger worlds primarily via two particular phenomena. First, the atmospheres of low-gravity bodies will expand outwards significantly as a function of surface temperature. The effective radiative surfaces of these atmospheres thus gain area relative to the surface, increasing both outgoing and absorbed radiation. Previously, it has been shown that this effect is already important for worlds of mass $\sim 0.1M_{\oplus}$ (Goldblatt, 2015). Second, the atmospheres of low-gravity worlds are extremely vulnerable to hydrodynamic escape via supersonic pressure-driven bulk outflow (Pierrehumbert, 2010; Catling and Kasting, 2017; Lehmer et al., 2017; Zahnle and Catling, 2017).

Kuramoto and Matsui (1994) studied these two phenomena together in the context of giant icy satellite accretion. That study used gray radiative transfer and focused on the implications for compositional evolution. Here, we present the first unified study of these phenomena and their implications for long-term atmospheric evolution and habitability over a wide range of starting parameters. To achieve this, we combine a line-by-line radiative-convective climate model (Section 2.3) with energy balance (Section 2.2) and hydrodynamic escape of the atmosphere (Section 2.4).

The habitability of low-gravity waterworlds has previously been discussed in a range of contexts. Examples of interest in our solar system include the icy Jovian moons Europa, Ganymede, and Callisto (Reynolds et al., 1983; Chyba and Phillips, 2001; Chela-Flores, 2010; Grasset et al., 2013), which all have masses below $0.03M_{\oplus}$. For terrestrial-mass worlds, habitability is characterized by a circumstellar “habitable zone”, commonly defined as the range in orbital distance within which the world could maintain surface liquid water. (Huang, 1959; Kasting et al., 1993; Kopparapu et al., 2013; Seager, 2013; Ramirez, 2018). The inner edge of the habitable zone is dependent on gravity even for near-terrestrial mass objects, albeit due to two competing mechanisms: as gravity decreases, a decreasing runaway greenhouse emission limit tends to move it further from the host star (Pierrehumbert, 2010; Kopparapu et al., 2014), and radiative surface expansion tends to move it closer (Goldblatt, 2015). However, because neither comets nor gas giants are habitable according to the criterion of surface liquid water (Kasting et al., 1993; Ramirez, 2018), there must ultimately also exist habitability boundaries in terms of gravity. The lower gravity boundary marks the transition between worlds that are effectively “planet-like” (habitable given the correct stellar flux), and “comet-like” (never habitable regardless of stellar flux). The model developed in this study provides a means to determine this boundary (Section 3.3), and shows what happens as it is approached.

In Section 3.1, we show that atmospheric expansion suppresses the runaway greenhouse and that escape rates determine the inner edge of the habitable zone for bodies of mass $\lesssim 0.1M_{\oplus}$. In Section 3.2, we show that a long-lived liquid water habitable zone is possible. In Section 3.4, we show that the ice-albedo feedback works against stellar flux-driven deglaciations to long-lived states. We discuss these results in Section 4 and conclude in Section 5.

2 Methods

2.1 Starting assumptions

We assume that the low-gravity waterworld has a pure water vapor atmosphere and a water reservoir fixed at $40\%$ of the planet’s total mass. For a given surface temperature, the atmospheric water content and temperature profile are set by the saturation vapor pressure. The water content of the atmosphere remains orders of magnitude below the world’s total water content in even the most extreme cases considered. We assume spherical symmetry; this is required by both the radiative transfer model and the hydrodynamic escape formulation. Our model, and therefore the “habitable zone” we consider, neglects the potential effects of CO2 cycling (Abbot et al., 2012; Ramirez and Levi, 2018), ocean chemistry (Wordsworth and Pierrehumbert, 2013; Kite and Ford, 2018), and three-dimensional dynamics (Pierrehumbert and Ding, 2016; Ding and Pierrehumbert, 2018).

It is convenient to parameterize the mass and radius of the waterworld in terms of a single quantity, surface gravity $g$. This is accomplished using the following scaling relation, accurate for waterworlds in the range $10^{-2}M_{\oplus}<M<M_{\oplus}$ (Sotin et al., 2007):

[TABLE]

2.2 Energy balance

A schematic of our model is shown in Figure 1. The waterworld gains energy from incoming stellar radiation, and loses it via longwave emission and atmospheric escape. We describe this by the following equation (2), which is similar to those of Pierrehumbert (2010) and Lehmer et al. (2017). However, we specifically account for radiative surface expansion in both the shortwave and longwave regimes:

[TABLE]

Here $A$ is the planetary albedo, $F_{stel}$ and $F_{out}$ are the incoming shortwave and outgoing longwave fluxes, respectively, $g_{s}$ and $r_{s}$ are the surface gravity and radius, respectively, $L$ is the latent heat of vaporization of water, $\Phi$ is the escaping mass flux, and $r_{SW}$ and $r_{LW}$ are, respectively the radii of the shortwave absorption and longwave emission surfaces. The critical radius $r_{c}$, where the outflow reaches the speed of sound $w_{c}$, is far enough above the radiative surfaces that the outflow near the radiative surfaces is negligible. Therefore, radiative transfer and hydrodynamic escape can be treated separately (Appendix A). In the low-gravity limit, the energy balance equation is consistent with expressions for volatile loss from comets (Lebofsky, 1975; Weissman, 1980): since there is no atmosphere, the scaling factors disappear, the outgoing flux reduces to the Stefan-Boltzmann law ($F_{out}=\sigma T_{s}^{4}$), and the gravity term ($g_{s}r_{s}$) becomes negligible. Conversely, in the high gravity limit, the scaling factors and escape terms are negligible, resulting in a standard runaway greenhouse scenario.

Here we obtain $r_{SW},r_{LW},F_{out}$ from line-by-line radiative transfer calculations (Section 2.3), and calculate $\Phi$ analytically, assuming hydrodynamic escape (Section 2.4). Then, (2) is solved numerically in instellation-gravity space.

2.3 Radiative transfer

To calculate outgoing longwave radiation and the radii of the radiative surfaces, we use the one-dimensional line-by-line radiative-convective model used previously by Schaefer et al. (2016) and Wordsworth et al. (2017). As in Schaefer et al. (2016), this model uses a truncated version of the HITEMP2010 line list (Rothman et al., 2010) accurate below 1000 K. We assume that the atmosphere is in phase equilibrium with the surface reservoir and therefore the temperature profile is given by the saturation vapor pressure curve of water. We limit ourselves to clear-sky calculations, accounting for cloud albedo forcing using the standard assumption (Goldblatt et al., 2013; Goldblatt, 2015) of an artificially large surface albedo (0.2). The potential effects of other cloud processes are discussed in Section 4.1.

The “radiative scaling factors” from (2), $(\frac{r_{SW}}{r_{s}})^{2}$ and $(\frac{r_{LW}}{r_{s}})^{2}$, are calculated using the method of Goldblatt (2015); we present a summary here. To a good approximation, emission of outgoing longwave radiation and absorption of incoming shortwave radiation occur at the altitude where optical depth $\tau=1$, measured from the top of the atmosphere. These altitudes can be calculated from the output of the radiative transfer model for a given surface gravity $g_{s}$ and temperature profile using the equation of hydrostatic balance (taking into account the decrease in gravity with altitude, which is non-negligible for small worlds):

[TABLE]

Here $r_{s}$ is the surface radius, $z$ the altitude, $p$ pressure, and $\rho$ mass density. Even for the most rapidly escaping atmospheres we consider, the velocity of the outflow near the radiative surface is small, such that hydrostatic balance can be assumed (Appendix A). Once the heights of unity optical depth as a function of wavenumber (spectral heights) $z_{\tau=1}(\nu)$ are obtained, the height of the radiative surface is calculated by a weighted average of $z_{\tau=1}(\nu)$ using the relevant irradiance $I(\nu)$.

[TABLE]

For the longwave emission surface, $I(\nu)$ is the outgoing longwave spectral irradiance given by the radiative transfer calculation. For the shortwave emission surface, it is the Planck blackbody spectral irradiance (hereafter, Planck curve) centered around the stellar effective emission temperature. The effect of deviations from blackbody behavior in real stellar spectra is neglected here.

The nature of the water vapor absorption spectrum is such that $z_{\tau=1}(\nu)$ shows a decreasing trend with increasing $\nu$. Therefore, this averaging process means that the emission surface is always above the absorption surface. Therefore, the longwave “scaling factor” is larger than its shortwave equivalent; indeed, their separation grows with temperature. This effect is critical to the rest of the study and is illustrated in Figure 2. The spectral heights of $\tau=1$ and the scaling factors produced by our model are in good agreement with those of Goldblatt (2015), who considered worlds with $0.12M_{\oplus}<M<10M_{\oplus}$.

2.4 Hydrodynamic escape

The problem of formulating the hydrodynamic escape from small bodies has been investigated previously. (Kuramoto and Matsui, 1994; Pierrehumbert, 2010; Catling and Kasting, 2017; Lehmer et al., 2017; Zahnle and Catling, 2017). Here we simplify the analysis by restricting ourselves to an isothermal escape formulation using the surface temperature $T_{s}$. This is because it is simple, easy to work with, and provides an upper bound on the escape rate. We are ultimately interested in water lifetime and habitability, and isothermal escape will provide a conservative estimate for both.

The isothermal sound speed is $\sqrt{\frac{kT_{s}}{m}}\equiv w_{c}$, where $k$ is Boltzmann’s constant and $m$ the molecular mass of water. Assuming transonic escape, there is a critical radius $r_{c}$ at which this speed is first reached. The outflow in spherically symmetric steady state is described by the continuity equation

[TABLE]

and the momentum conservation equation

[TABLE]

where $\rho$ is density, $w$ is vertical velocity, $p$ is pressure, $G$ is Newton’s gravitational constant, and $M$ is the mass of the body. These are combined (e.g. Catling and Kasting 2017, Lehmer et al. 2017) to form the isothermal planetary wind equation

[TABLE]

where the ideal gas law has also been used. At the critical radius $w=w_{c}$ and (7) yields

[TABLE]

Typically, $r_{c}\gtrsim 10r_{s}$. Following Lehmer et al. (2017), we integrate (7) from $r_{s}$ to $r_{c}$ and neglect the small kinetic energy term near the surface. This yields the mass flux, $\Phi$:

[TABLE]

$\Phi$ is a function of $T_{s}$ via $w_{c}$. The rate of mass loss is calculated as

[TABLE]

and the water lifetime $\tau_{water}$ is given by

[TABLE]

3 Results

From the radiative transfer calculations, we generate outgoing radiation and scaling factors over a $g-T_{s}$ grid. For a given stellar flux $F_{stel}$ and gravity $g$, we can now numerically solve (2) for $T_{s}$, the steady state surface temperature. The steady-state surface temperature further determines escape rates and water lifetimes, which form the basis of the results presented here.

3.1 No runaway greenhouse

Terrestrial-mass worlds can experience the phenomenon known as the “runaway greenhouse” (Komabayasi, 1967; Ingersoll, 1969; Kasting, 1988; Nakajima et al., 1992; Goldblatt and Watson, 2012): at high enough surface temperatures, total emission converges towards a limit. If the incoming stellar flux exceeds this threshold, there will be a runaway surface temperature increase, boiling off the oceans - making the world uninhabitable. Therefore, the runaway greenhouse is believed to set the ultimate inner edge of the habitable zone. (Kasting et al., 1993; Goldblatt and Watson, 2012; Kopparapu et al., 2013; Goldblatt et al., 2013)

Based on (2) and the preceding analysis, it is clear that low-gravity waterworlds should behave differently. Although the flux per unit area from the effective emission surface will still reach a limit as temperature increases, the emission surface itself can keep gaining area relative to the absorption surface (Figure 2). If the ratio between the longwave and shortwave scaling factors does not reach a limit before escape becomes important, a runaway greenhouse limit will not be reached.

We investigate this quantitatively by studying the steady-state instellation $F_{stel}$ as a function of surface temperature. (2) can be rewritten in the form

[TABLE]

$F_{stel}(T)$ is plotted in Figure 3, for a given $g=2.3$ ms*-2*, corresponding to $M\simeq 0.08M_{\oplus}$. This value is illustrative because the escape rate increases slowly enough with $T$ that the radiative surface expansion effect can be clearly observed. Water lifetimes are also added: they are a nonlinear function of temperature. For comparison, we also plot the results of two other radiation assumptions: the same line-by-line calculation with radiative surface expansion ignored, and blackbody emission from the surface ($\sigma T_{s}^{4}$) until the runaway greenhouse limit is reached. The latter approach was used by Lehmer et al. (2017), to study Europa- and Ganymede-like worlds. Hydrodynamic escape is assumed to be governed by (9) in all cases. This is responsible for the rapid increase in $F_{stel}(T)$ as higher temperatures are approached. We note that such states with high mass-flux are “quasi-steady”: the rapid expulsion of water into space is not sustainable indefinitely.

A given steady-state instellation will allow for a habitable (long-lived liquid water) state if the surface temperature is larger than 273 K and the water inventory is long-lived enough. In Figure 3 and the rest of this study we choose $10^{9}$ years as a cut-off; this choice is further discussed in Section 3.2. In our model, where radiative surface expansion is considered, $F_{stel}(T)$ never reaches a limit. Therefore, there is no runaway greenhouse. This leads to a relatively wide habitable zone. In contrast, neglecting radiative surface expansion but still using the line-by-line model would predict an extremely thin habitable zone, and assuming blackbody emission until the runaway greenhouse would predict no habitable zone at all.

3.2 Long-lived liquid water habitable zone

We now expand the picture of Section 3.1 and Figure 3 into gravity space, allowing the instellation to be set by the orbital distance. We choose a water lifetime of $10^{9}$ years as an inner habitability boundary. While the principal requirement for the water lifetime habitability boundary is that it feasibly allows the evolution of life, this particular choice of timescale also ensures that such worlds have a reasonable chance of being habitable at the time of observation - even when the system’s age is uncertain. With this definition, the “longer-lived states” in our model consistently have surface temperatures below proposed upper bounds for survival of life (Bains et al., 2015), suggesting that temperature is not additionally prohibitive.

As can be observed from Figure 2, stellar emission temperature plays a role in determining the size of the shortwave scaling factor. We thus consider two limiting cases: an M-star (modeled on AD Leonis) and a G-star (modeled on the Sun). The resulting steady-state temperatures and the $10^{9}$ year water lifetime contour are shown in Figure 4. We label five qualitatively different climate states: long-lived snowball, long-lived temperate, short-lived snowball, short-lived temperate and short-lived hot. The mass of Ganymede and a $10^{7}$ year water lifetime contour are plotted additionally for reference.

The “habitable zone” is defined by the area enclosed by the 273 K temperature contour and the $10^{9}$ year water lifetime contour. The conservative maximum width in the low-gravity regime is $\sim 0.08$ AU for G-stars and $\sim 0.008$ AU for M-stars. Relaxing the $10^{9}$ year habitability boundary would extend the habitable zone both spatially and to lower gravities. Worlds with appropriate geochemical cycling could of course have a wider habitable zone, because the outer edge would move outwards.

3.3 The planet-comet transition

The width of the habitable zone shown in Figure 4 decreases with gravity, ultimately disappearing. This marks the transition from “planet-like” (habitable at some stellar flux) to “comet-like” (never habitable regardless of stellar flux). For worlds with a mass-radius scaling accurately described by (1), our model conservatively gives this boundary at $g\simeq 1.48$ ms*-2*, or $M/M_{\oplus}\simeq 0.0268$. The boundary is essentially independent of stellar spectral type, because the differences in the shortwave scaling factor are negligible for waterworld surface temperatures near 273 K.

3.4 Stellar flux-driven deglaciation

Although our model predicts a habitable climate state, such a state is not necessarily easy to access. Figure 4 exhibits a range of possible climate states: short-lived hot, temperate or snowball, and long-lived temperate or snowball. The ice-albedo feedback can hamper transitions from snowball states to temperate states: this has already been demonstrated for terrestrial-mass worlds (Yang et al., 2017). We can incorporate the ice-albedo feedback into our model using the simple albedo step function

[TABLE]

where $\alpha_{I}$ is the albedo of the icy (snowball) state, and $\alpha_{L}$ is the albedo when there is surface liquid water. Hysteresis plots for different choices of $\alpha_{I}$ are shown in Figure 5. We observe that a snowball state experiencing a stellar flux-driven deglaciation generally bypasses the long-lived state entirely, except for very low $\alpha_{I}$ values. Although the mechanism setting the inner edge of the habitable zone is different, the conclusion of the habitable state likely being bypassed upon stellar flux-driven deglaciation is the same as that of Yang et al. (2017).

4 Discussion

4.1 Habitability: other considerations

Beyond those considered in this study, other factors play a role in determining the habitability of low-gravity waterworlds. For example, we have restricted ourselves to pure water vapor atmospheres; in reality, a variety of species could be present, depending on interior and atmospheric chemistry. We did not consider the potential effects of clouds except via an increased planetary albedo. Cloud greenhouse forcing could reduce outgoing radiation as it does for terrestrial-mass worlds, and clouds in an expanded low-gravity atmosphere could increase absorbed flux by scattering more stellar radiation towards the surface. Depending on whether cloud effects contribute more to cooling or warming, the size of the habitable zone presented in this study could either increase or decrease. The net effect should be carefully assessed in the future using three-dimensional models (e.g. Pierrehumbert and Ding 2016). Finally, as in our solar system, low-gravity waterworlds may often be found as exomoons around giant planets. Then, tidal forces are important for habitability (Heller and Barnes, 2013), becoming especially deleterious around M-stars (Zollinger et al., 2017).

Even if a habitable state exists, it is not necessarily “accessible”, as we showed in Section 3.4. Worlds with a warm history have access to the entire habitable zone derived here; the proto-atmosphere phenomenon discussed by Kuramoto and Matsui (1994) could make a warm start more likely. However, for worlds with a cold history, stellar-flux-driven deglaciation into a long-lived state only occurs for artificially low ice albedo values. Somewhat more favorable conditions might be found around M-stars, where the ice-albedo feedback is weakened (Joshi and Haberle, 2012; Shields et al., 2013). However, such systems experience another problem: the dramatic pre-main sequence evolution of the star may strip potentially habitable worlds of their water early on (Ramirez and Kaltenegger, 2014; Luger and Barnes, 2015; Tian and Ida, 2015). Further research is required to integrate all of these considerations.

In this study we have used the conventional but limited definition of habitability as the ability to maintain surface liquid water on sufficient timescales. While our model categorizes worlds of Europa’s or Ganymede’s size as “comet-like” (unable to maintain surface liquid water for Gyr timescales regardless of stellar flux), that does not necessarily make them uninhabitable. For example, there is a large body of work on how life could survive on Europa within a subsurface ocean, supported by energy from a geologically active interior and favorable chemistry (Reynolds et al., 1983; McCollom, 1999; Gaidos et al., 1999; Chyba, 2000). Detection of such life remotely would be extremely challenging. Thus the question of life on “comet-like” worlds ultimately remains open, if unanswerable at present.

4.2 Observing low-gravity waterworlds

Observation of low-gravity waterworlds outside of our solar system is possible in principle. The smallest exoplanet discovered to date is the rocky Kepler 37-b, with a radius of $0.303R_{\oplus}$ (Barclay et al., 2013). This is already smaller than all of the low-gravity waterworlds modeled in this study. Regarding exomoons, a previous search using the Kepler Space Telescope (Kipping et al., 2015) was sensitive to masses as low as $1.7$ $M_{\text{Ganymede}}\simeq 0.043$ $M_{\oplus}$, which includes most of the mass range for habitable low-gravity waterworlds shown in Figure 4. Unfortunately, due to the nature of transit photometry, low-gravity worlds can currently only be found in orbits with relatively small semi-major axes; this is likely to preclude long-lived water. Correctly classifying waterworlds presents another issue: under current observational uncertainties in this size regime (Jontof-Hutter et al., 2015), mass-density determinations alone are likely to be insufficient (Sotin et al., 2007).

The phenomena described in this study will have observable effects. Radiative surface expansion increases transit depths; planetary density will be underestimated by a factor of $r_{SW}^{3}/r^{3}$, a decrease of up to $15\%$ for temperate worlds. Clouds or escaping dust will amplify this: the latter already provides a means to observe the much smaller exo-comets (Kiefer et al, 2014; Rappaport et al., 2017), and has also been proposed to explain anomalously low-density sub-Neptunes (Wang and Dai, 2019). Escaping H atoms could produce a detectable Ly-$\alpha$ signal; for exomoons, these atoms would likely accumulate in a torus around the host planet (Lehmer et al., 2017). Transmission spectroscopy has already been successfully used to study Earth-mass planets (e.g. de Wit et al. 2016); for low-gravity waterworlds, the highly expanded atmosphere will amplify the variations in such a spectrum. Coupling the model developed here with a forward model for observables would be an interesting follow-up to this study, and would aid future detection efforts.

4.3 Future directions

Further work could consider more complicated models of hydrodynamic escape. Since the isothermal escape approximation represents an upper bound, such work is likely to widen the habitable zone presented in Figure 4.

The two key phenomena governing evolution on low-gravity waterworlds, radiative surface expansion and hydrodynamic escape, are treated separately in our model. This is possible while the outflow near the radiative surface is negligible, which is the case for the parameter space regime we consider. While our model is sufficient to estimate constraints on habitability, the scenario of non-negligible outflow near the radiative surface remains of physical interest, and should be explored further.

Finally, the planet-comet transition calculated in this study represents a hard lower gravity boundary for habitability in the conventional sense. An upper gravity boundary is likely given by the onset of runaway gas accretion, which distinguishes super-Earths from gas giants (e.g. Ginzburg et al. 2016). Such boundaries, when used in tandem, could be a useful tool when estimating a total number of potentially habitable worlds.

5 Conclusion

Low-gravity waterworlds experience radiative surface expansion (Goldblatt, 2015) and hydrodynamic escape (Kuramoto and Matsui, 1994; Pierrehumbert, 2010; Catling and Kasting, 2017; Zahnle and Catling, 2017; Lehmer et al., 2017). This study is the first to explore their combined effect on long-term atmospheric evolution and habitability. We showed that radiative surface expansion suppresses the runaway greenhouse (Section 3.1), the existence and location of a surface liquid water habitable zone whose inner edge is set by escape (Section 3.2), and that the ice-albedo feedback works against stellar-flux-driven deglaciations to long-lived states (Section 3.4). Combining these effects has allowed us for the first time to calculate the transition from “planet-like” to “comet-like” regimes for the evolution of a world’s H2O inventory (Section 3.3). Comet-like worlds cannot sustain surface liquid water on Gyr timescales. Therefore, by the criterion of long-lived surface liquid water, this provides a hard lower boundary on habitability in terms of gravity.

C. W. A. acknowledges support from the Harvard Origins of Life Initiative, and thanks E. Kite and P. Niraula for helpful discussions. This work was supported by NASA Habitable Worlds grant NNX16AR86G.

Appendix A Hydrostatic balance at the radiative surface

As stated in Section 2.2 and Section 2.3, our model can treat radiative surface expansion and hydrodynamic escape separately because the velocity of the outflow at and below the radiative surface is negligible. Specifically, it is small enough to allow a calculation of the radiative surface height using the hydrostatic relation (Section 2.3). Here, we make this argument quantitatively.

For spherically symmetric hydrodynamic escape, the most complete momentum equation is

[TABLE]

where the gravity $g$, pressure $p$, outflow $w$ are all functions of $r$. This is accurately approximated by a hydrostatic relation $\frac{dp}{dz}=-\rho g$ if the inertial term $\rho w\frac{dw}{dr}$ is sufficiently small. Thus it suffices to check whether the dimensionless number

[TABLE]

is much smaller than 1 at and below the (longwave) radiative surface. We can obtain an analytical expression for $w\frac{dw}{dr}$ by considering an integrated solar wind equation for isothermal escape (Parker, 1963; Catling and Kasting, 2017):

[TABLE]

For $w/w_{c}\lesssim 0.4$, it is a good approximation to drop the $w^{2}$ term and write

[TABLE]

Hence,

[TABLE]

Across the parameter range presented in Figure 4, using the steady state temperatures, we can thus calculate $w$ and $w\frac{dw}{dr}$ at the radiative surface using (A4) and (A5). In this regime, it is the case that $w\frac{dw}{dr}$ monotonically increases until the radiative surface. Therefore, when considering the accuracy of the hydrostatic approximation in calculating the radiative surface height, it is ultimately sufficient to evaluate $w\frac{dw}{dr}/g$ at the longwave radiative surface, and verify that it is much smaller than 1. In Figure 6, we plot logarithms of $w/w_{c}$ and $w\frac{dw}{dr}/g$ at the longwave radiative surface, over the same $g$-distance regime covered by Figure 4, finding that both of these terms are much smaller than 1. Since $w/w_{c}\ll 1$ at the longwave radiative surface, the use of (A5) to calculate $w\frac{dw}{dr}/g$ is justified. Since $w\frac{dw}{dr}/g\ll 1$ at the longwave radiative surface, the use of the hydrostatic approximation in Section 2.3 is justified.