The macroturbulent atmospheric circulation of Earth-like planets mediates their equator-to-pole heat transport. For fast-rotating terrestrial planets, baroclinic instabilities in the mid-latitudes lead to turbulent eddies that act to transport heat poleward. In this work, we derive a scaling theory for the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanet atmospheres. This theory is built on the work of Jansen & Ferrari and determines how unstable the atmosphere is to baroclinic instability (the baroclinic “criticality”) through a balance between the baroclinic eddy heat flux and radiative heating/cooling. We compare our scaling theory to General Circulation Model (GCM) simulations and find that the theoretical predictions for equator-to-pole temperature contrast and bulk lapse rate broadly agree with GCM experiments with varying rotation rate and surface pressure throughout the baroclincally unstable regime. Our theoretical results show that baroclinic instabilities are a strong control of heat transport in the atmospheres of Earth-like exoplanets, and our scalings can be used to estimate the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanets. These scalings can be tested by spectroscopic retrievals and full-phase light curves of terrestrial exoplanets with future space telescopes.
We implement a search for exoplanets in campaigns zero through eight (C08) of the K2 extension of the Kepler spacecraft. We apply a modified version of the QATS planet search algorithm to K2 light curves produced by the EVEREST pipeline, carrying out the C08 search on 1.5 × 105 target stars with magnitudes in the range of Kp = 9?15. We detect 818 transiting planet candidates, of which 374 were undiscovered by prior searches, with {64, 15, 5, 2, 1} in {2, 3, 4, 5, 6}-planet multiplanet candidate systems, respectively. Of the new planets detected, 100 orbit M dwarfs, including one that is potentially rocky and in the habitable zone. A total of 154 of our candidates reciprocally transit with our solar system: they are geometrically aligned to see at least one solar system planet transit. We find candidates that display transit timing variations and dozens of candidates on both period extremes with single transits or ultrashort periods. We point to evidence that our candidates display similar patterns in frequency and sizeperiod relation to confirmed planets, such as tentative evidence for the radius gap. Confirmation of these planet candidates with follow-up studies will increase the number of K2 planets by up to 50%, and characterization of their host stars will improve statistical studies of planet properties. Our sample includes many planets orbiting bright stars amenable for radial velocity follow-up and future characterization with JWST. We also list the 579 eclipsing binary systems detected as part of this search.
The membranes of the first protocells on the early Earth were likely self-assembled from fatty acids. A major challenge in understanding how protocells could have arisen and withstood changes in their environment is that fatty acid membranes are unstable in solutions containing high concentrations of salt (such as would have been prevalent in early oceans) or divalent cations (which would have been required for RNA catalysis). To test whether the inclusion of amino acids addresses this problem, we coupled direct techniques of cryoelectron microscopy and fluorescence microscopy with techniques of NMR spectroscopy, centrifuge filtration assays, and turbidity measurements. We find that a set of unmodified, prebiotic amino acids binds to prebiotic fatty acid membranes and that a subset stabilizes membranes in the presence of salt and Mg2+. Furthermore, we find that final concentrations of the amino acids need not be high to cause these effects; membrane stabilization persists after dilution as would have occurred during the rehydration of dried or partially dried pools. In addition to providing a means to stabilize protocell membranes, our results address the challenge of explaining how proteins could have become colocalized with membranes. Amino acids are the building blocks of proteins, and our results are consistent with a positive feedback loop in which amino acids bound to self-assembled fatty acid membranes, resulting in membrane stabilization and leading to more binding in turn. High local concentrations of molecular building blocks at the surface of fatty acid membranes may have aided the eventual formation of proteins.
coronagraph is an open-source Python package for generalized telescope noise modeling for extrasolar planet (exoplanet) science. This package is based on Interactive Data Language (IDL) code originally developed by T. Robinson (Robinson, 2018), and described in detail with science applications in (Robinson, Stapelfeldt, & Marley, 2016).
This article presents the proportional evolutionary time (PET) hypothesis, which posits that the mean time required for the evolution of complex life is a function of stellar mass. The biological available window is defined as the region of a stellar spectrum between 200 and 1200?nm that generates free energy for life. Over the ?4 Gyr history of Earth, the total energy incident at the top of the atmosphere and within the biological available window is ?1034 J. The hypothesis assumes that the rate of evolution from the origin of life to complex life is proportional to this total energy, which would suggest that planets orbiting other stars should not show signs of complex life if the total energy incident on the planet is below this energy threshold. The PET hypothesis predicts that late K- and M-dwarf stars (M??0.7 ) represent the best targets for the next generation of space telescopes to search for spectroscopic biosignatures indicative of complex life.
We investigate planetesimal accretion via a direct N-body simulation of an annulus at 1 au orbiting a 1 M? star. The planetesimal ring, which initially contains N = 106 bodies is evolved into the oligarchic growth phase. Unlike previous lower resolution studies, we find that the mass distribution of planetesimals develops a bump at intermediate mass after the oligarchs form. This feature marks a boundary between growth modes. The smallest planetesimals are packed tightly enough together to populate mean motion resonances with the oligarchs, which heats the small bodies, enhancing their growth. If we depopulate most of the resonances by decreasing the width of the annulus, this effect becomes weaker. To clearly demonstrate the dynamics driving these growth modes, we also examine the evolution of a planetary embryo embedded in an annulus of collisionless planetesimals. In this case, we find that the resonances push planetesimals away from the embryo, decreasing the surface density of the bodies adjacent to the embryo. This effect only occurs when the annulus is wide enough and the mass resolution of the planetesimals is fine enough to populate the resonances. The bump we observe in the mass distribution resembles the 100 km power-law break seen in the size distribution of asteroid belt objects. Although the bump produced in our simulations occurs at a size larger than 100 km, we show that the bump location is sensitive to the initial planetesimal mass, which implies that this feature is potentially useful for constraining planetesimal formation models.
Using Gaussian process regression to analyze the Martian surface methane Tunable Laser Spectrometer (TLS) data reported by Webster et al. (2018), we find that the TLS data, taken as a whole, do not indicate seasonal variability. Enrichment protocol CH4 data are consistent with either stochastic variation or a spread of periods without seasonal preference.
Thermal history models, historically used to understand Earth’s geologic history, are being coupled to climate models to map conditions that allow planets to maintain life. However, the lack of structural uncertainty assessment has blurred guidelines for how thermal history models can be used toward this end. Structural uncertainty is intrinsic to the modeling process. Model structure refers to the cause and effect relations that define a model and are assumed to adequately represent a particular real world system. Intrinsic/structural uncertainty is different from input and parameter uncertainties (which are often evaluated for thermal history models). A full uncertainty assessment requires that input/parametric and intrinsic/structural uncertainty be evaluated (one is not a substitute for the other). We quantify the intrinsic uncertainty for several parameterized thermal history models (a subclass of planetary models). We use single perturbation analysis to determine the reactance time of different models. This provides a metric for how long it takes low?amplitude, unmodeled effects to decay or grow. Reactance time is shown to scale inversely with the strength of the dominant model feedback (negative or positive). A perturbed physics analysis is then used to determine uncertainty shadows for model outputs. This provides probability distributions for model predictions. It also tests the structural stability of a model (do model predictions remain qualitatively similar, and within assumed model limits, in the face of intrinsic uncertainty?). Once intrinsic uncertainty is accounted for, model outputs/predictions and comparisons to observational data should be treated in a probabilistic way.
The search for an inhabited planet, beyond our own, is a driver of planetary exploration in our solar system and beyond. Using information from our own planet to inform search strategies allows for a targeted search. It is, however, worth considering some span in the strategy and in a priori expectation. An inhabited, Earth-like planet is one that would be similar to Earth in ways that extend beyond having biota. To facilitate a comparative cost/risk/benefit analysis of different potential search strategies, we use a metric akin to the Earth-similarity index. The metric extends from zero, for an inhabited planet that is like Earth in all other regards (i.e., zero differences), toward end-member values for planets that differ from Earth but maintain life potential. The analysis shows how finding inhabited planets that do not share other Earth characteristics could improve our ability to assess galactic life potential without a large increase in time-commitment costs. Search strategies that acknowledge the possibility of such planets can minimize the potential of exploration losses (e.g., searching for long durations to reach conclusions that are biased). Discovering such planets could additionally provide a test of the Gaia hypothesisa test that has proven difficult when using only Earth as a laboratory. Finally, we discuss how an Earth2.0 narrative that has been presented to the public as a search strategy comes with nostalgia-laden baggage that does not best serve exploration.
Context. The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised for follow-up observations of transiting planet candidate Kepler objects of interest with large transit timing variations (TTVs). The main goal of KOINet is the completion of their TTV curves as the Kepler telescope stopped observing the original Kepler field in 2013.
Aims. We ensure a comprehensive characterisation of the investigated systems by analysing Kepler data combined with new ground-based transit data using a photodynamical model. This method is applied to the Kepler-82 system leading to its first dynamic analysis.
Methods. In order to provide a coherent description of all observations simultaneously, we combine the numerical integration of the gravitational dynamics of a system over the time span of observations with a transit light curve model. To explore the model parameter space, this photodynamical model is coupled with a Markov chain Monte Carlo algorithm.
Results. The Kepler-82b/c system shows sinusoidal TTVs due to their near 2:1 resonance dynamical interaction. An additional chopping effect in the TTVs of Kepler-82c hints to a further planet near the 3:2 or 3:1 resonance. We photodynamically analysed Kepler long- and short-cadence data and three new transit observations obtained by KOINet between 2014 and 2018. Our result reveals a non-transiting outer planet with a mass of mf = 20.9 ± 1.0 M? near the 3:2 resonance to the outermost known planet, Kepler-82c. Furthermore, we determined the densities of planets b and c to the significantly more precise values ?b = 0.98?0.14+0.10 g cm?3 and ?c = 0.494?0.077+0.066 g cm?3.
The current goals of the astrobiology community are focused on developing a framework for the detection of biosignatures, or evidence thereof, on objects inside and outside of our solar system. A fundamental aspect of understanding the limits of habitable environments (surface liquid water) and detectable signatures thereof is the study of where the boundaries of such environments can occur. Such studies provide the basis for understanding how a once inhabitable planet might come to be uninhabitable
The discovery of planets orbiting stars other than the Sun has accelerated over the past decade, and this trend will continue as new space- and ground-based observatories employ next-generation instrumentation to search the skies for habitable worlds. However, many factors and processes can affect planetary habitability and must be understood to accurately determine a planet’s habitability potential. While climate models have long been used to understand and predict climate and weather patterns on the Earth, a growing community of researchers has begun to apply these models to extrasolar planets. This work has provided a better understanding of how orbital, surface, and atmospheric properties affect planetary climate and habitability; how these climatic effects might change for different stellar and planetary environments; and how the habitability and observational signatures of newly discovered planets might be influenced by these climatic factors. This review summarizes the origins and evolution of the burgeoning field of exoplanet climatology, discusses recent work using a hierarchy of computer models to identify those planets most capable of supporting life, and offers a glimpse into future directions of this quickly evolving subfield of exoplanet science.
The partial pressure of atmospheric hydrogen (pH2) on the early Earth is important because it has been proposed that high pH2 warmed the planet or allowed prebiotic chemistry in the early atmosphere. However, such hypotheses lack observational constraints on pH2 . Here, we use the existence of detrital magnetites in (? 3.0 Ga) Archean riverbeds to constrain pH2 . Under the condition of high pH2 , magnetite should disappear via reductive dissolution. We investigated the timescale for a magnetite particle in a river to dissolve, which depends on pH2 and pCO2 . Using published estimates of Archean pCO2 and assuming the presence of Fe(III)-reducing microbes, the survival timescale is ? 1 kyr when pH2 is ?10-2 bar , and decreases as pH2 increases. Considering that the residence time of a particle in a short river (< 1000 km) is ? 1 kyr , the existence of detrital magnetite particles in Archean riverbeds likely indicates that pH2 was below ?10-2 bar . Such a level would preclude H2 as a greenhouse gas or a strongly reducing Archean atmosphere. It is also consistent with limits imposed on H2 by consumption by methanogens because conversion to CH4 is thermodynamically favored.
Super-Earths planets with sizes between the Earth and Neptune are found in tighter orbits than that of the Earth around more than one third of main sequence stars. It has been proposed that super-Earths are scaled-up terrestrial planets that also formed similarly, through mutual accretion of planetary embryos, but in discs much denser than the solar protoplanetary disc. We argue instead that terrestrial planets and super-Earths have two clearly distinct formation pathways that are regulated by the pebble reservoir of the disc. Through numerical integrations, which combine pebble accretion and N-body gravity between embryos, we show that a difference of a factor of two in the pebble mass flux is enough to change the evolution from the terrestrial to the super-Earth growth mode. If the pebble mass flux is small, then the initial embryos within the ice line grow slowly and do not migrate substantially, resulting in a widely spaced population of approximately Mars-mass embryos when the gas disc dissipates. Subsequently, without gas being present, the embryos become unstable due to mutual gravitational interactions and a small number of terrestrial planets are formed by mutual collisions. The final terrestrial planets are at most five Earth masses. Instead, if the pebble mass flux is high, then the initial embryos within the ice line rapidly become sufficiently massive to migrate through the gas disc. Embryos concentrate at the inner edge of the disc and growth accelerates through mutual merging. This leads to the formation of a system of closely spaced super-Earths in the five to twenty Earth-mass range, bounded by the pebble isolation mass. Generally, instabilities of these super-Earth systems after the disappearance of the gas disc trigger additional merging events and dislodge the system from resonant chains. Therefore, the key difference between the two growth modes is whether embryos grow fast enough to undergo significant migration. The terrestrial growth mode produces small rocky planets on wider orbits like those in the solar system whereas the super-Earth growth mode produces planets in short-period orbits inside 1 AU, with masses larger than the Earth that should be surrounded by a primordial H/He atmosphere, unless subsequently lost by stellar irradiation. The pebble flux which controls the transition between the two growth modes may be regulated by the initial reservoir of solids in the disc or the presence of more distant giant planets that can halt the radial flow of pebbles.
The discovery of the first interstellar object passing through the Solar System, 1I/2017 U1 (Oumuamua), provoked intense and continuing interest from the scientific community and the general public. The faintness of Oumuamua, together with the limited time window within which observations were possible, constrained the information available on its dynamics and physical state. Here we review our knowledge and find that in all cases, the observations are consistent with a purely natural origin for Oumuamua. We discuss how the observed characteristics of Oumuamua are explained by our extensive knowledge of natural minor bodies in our Solar System and our current knowledge of the evolution of planetary systems. We highlight several areas requiring further investigation.
Cyanide plays a critical role in origin of life hypotheses that have received strong experimental support from cyanide-driven synthesis of amino acids, nucleotides, and lipid precursors. However, relatively high cyanide concentrations are needed. Such cyanide could have been supplied by reaction networks in which hydrogen cyanide in early Earths atmosphere reacted with iron to form ferrocyanide salts, followed by thermal decomposition of ferrocyanide salts to cyanide. Using an aqueous model supported by new experimental data, we show that sodium ferrocyanide salts precipitate in closed-basin, alkaline lakes over a wide range of plausible early Earth conditions. Such lakes were likely common on the early Earth because of chemical weathering of mafic or ultramafic rocks and evaporative concentration. Subsequent thermal decomposition of sedimentary sodium ferrocyanide yields sodium cyanide (NaCN), which dissolves in water to form NaCN-rich solutions. Thus, geochemical considerations newly identify a particular geological setting and NaCN feedstock nucleophile for prebiotic chemistry.
Terrestrial planets orbiting M dwarfs may soon be observed with the James Webb Space Telescope (JWST) to characterize their atmospheric composition and search for signs of habitability or life. These planets may undergo significant atmospheric and ocean loss due to the superluminous pre-main-sequence phase of their host stars, which may leave behind abiotically generated oxygen, a false positive for the detection of life. Determining if ocean loss has occurred will help assess potential habitability and whether or not any O2 detected is biogenic. In the solar system, differences in isotopic abundances have been used to infer the history of ocean loss and atmospheric escape (e.g., Venus, Mars). We find that isotopologue measurements using transit transmission spectra of terrestrial planets around late-type M dwarfs like TRAPPIST-1 may be possible with JWST, if the escape mechanisms and resulting isotopic fractionation were similar to Venus. We present analyses of post-ocean-loss O2- and CO2-dominated atmospheres containing a range of trace gas abundances. Isotopologue bands are likely detectable throughout the near-infrared (18 ?m), especially 34 ?m, although not in CO2-dominated atmospheres. For Venus-like D/H ratios 100 times that of Earth, TRAPPIST-1b transit signals of up to 79 ppm are possible by observing HDO. Similarly, 18O/16O ratios 100 times that of Earth produce signals at up to 94 ppm. Detection at signal-to-noise ratio = 5 may be attained on these bands with as few as four to 11 transits, with optimal use of JWST’s NIRSpec Prism. Consequently, H2O and CO2 isotopologues could be considered as indicators of past ocean loss and atmospheric escape for JWST observations of terrestrial planets around M dwarfs.
The James Webb Space Telescope (JWST) will offer the first opportunity to characterize terrestrial exoplanets with sufficient precision to identify high mean molecular weight atmospheres, and TRAPPIST-1’s seven known transiting Earth-sized planets are particularly favorable targets. To assist community preparations for JWST observations, we use simulations of plausible post-ocean-loss and habitable environments for the TRAPPIST-1 exoplanets, and test simulations of all bright object time-series spectroscopy modes and all Mid-Infrared Instrument photometry filters to determine optimal observing strategies for atmospheric detection and characterization using both transmission and emission observations. We find that transmission spectroscopy with the Near-Infrared Spectrograph Prism is optimal for detecting terrestrial, CO2-containing atmospheres, potentially in fewer than 10 transits for all seven TRAPPIST-1 planets, if they lack high-altitude aerosols. If the TRAPPIST-1 planets possess Venus-like H2SO4 aerosols, up to 12 times more transits may be required to detect an atmosphere. We present optimal instruments and observing modes for the detection of individual molecular species in a given terrestrial atmosphere and an observational strategy for discriminating between evolutionary states. We find that water may be prohibitively difficult to detect in both Venus-like and habitable atmospheres, due to its presence lower in the atmosphere where transmission spectra are less sensitive. Although the presence of biogenic O2 and O3 will be extremely challenging to detect, abiotically produced oxygen from past ocean loss may be detectable for all seven TRAPPIST-1 planets via O2O2 collisionally induced absorption at 1.06 and 1.27 ?m, or via NIR O3 features for the outer three planets. Our results constitute a suite of hypotheses on the nature and detectability of highly evolved terrestrial exoplanet atmospheres that may be tested with JWST.
The habitable zone (HZ) is commonly defined as the range of distances from a host star within which liquid water, a key requirement for life, may exist on a planet’s surface. Substantially more CO2 than present in Earth’s modern atmosphere is required to maintain clement temperatures for most of the HZ, with several bars required at the outer edge. However, most complex aerobic life on Earth is limited by CO2 concentrations of just fractions of a bar. At the same time, most exoplanets in the traditional HZ reside in proximity to M dwarfs, which are more numerous than Sun-like G dwarfs but are predicted to promote greater abundances of gases that can be toxic in the atmospheres of orbiting planets, such as carbon monoxide (CO). Here we show that the HZ for complex aerobic life is likely limited relative to that for microbial life. We use a 1D radiative-convective climate and photochemical models to circumscribe a Habitable Zone for Complex Life (HZCL) based on known toxicity limits for a range of organisms as a proof of concept. We find that for CO2 tolerances of 0.01, 0.1, and 1 bar, the HZCL is only 21%, 32%, and 50% as wide as the conventional HZ for a Sun-like star, and that CO concentrations may limit some complex life throughout the entire HZ of the coolest M dwarfs. These results cast new light on the likely distribution of complex life in the universe and have important ramifications for the search for exoplanet biosignatures and technosignatures.
Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location of water vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.
Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location of water vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.
We describe the circumstances that led to the discovery of Kepler-36b, and the subsequent characterization of its host planetary system. The Kepler-36 system is remarkable for its physical properties: the close separation of the planets, the contrasting densities of the planets despite their proximity, and the short chaotic timescale. Its discovery and characterization was also remarkable for the novelty of the detection technique and for the precise characterization due to the large transit-timing variations caused by the close proximity of the planets, as well as the precise stellar parameters due to asteroseismology. This was the first multi-planet system whose transit data was processed using a fully consistent photometric-dynamical model, using population Markov Chain Monte Carlo techniques to precisely constrain system parameters. Amongst those parameters, the stellar density was found to be consistent with a complementary, concurrent asteroseismic analysis. In a first, the 3D orientation of the planets was constrained from the lack of transit-duration variations. The system yielded insights into the composition and evolution of short-period planet systems. The denser planet appears to have an Earth-like composition, with uncertainties comparable to the highest precision rocky exoplanet measurements, and the planet densities foreshadowed the rocky/gaseous boundary. The formation of this system remains a mystery, but should yield insights into the migration and evolution of compact exoplanet systems.
Kepler-62f is the first exoplanet small enough to plausibly have a rocky composition orbiting within the habitable zone (HZ) discovered by the Kepler Mission. The planet is 1.4 times the size of the Earth and has an orbital period of 267 days. At the time of its discovery, it had the longest period of any small planet in the habitable zone of a multi-planet system. Because of its long period, only four transits were observed during Kepler’s interval of observations. It was initially missed by the Kepler pipeline, but the first three transits were identified by an independent search by Eric Agol, and it was identified as a planet candidate in subsequent Kepler catalogs. However in the latest catalog of exoplanets (Thompson et al., 2018), it is labeled as a false positive. Recent exoplanet catalogues have evolved from subjective classification to automatic classifications of planet candidates by algorithms (such as Robovetter). While exceptionally useful for producing a uniform catalogue, these algorithms sometimes misclassify planet candidates as a false positive, as is the case of Kepler-62f. In particularly valuable cases, i.e., when a small planet has been found orbiting in the habitable zone (HZ), it is important to conduct comprehensive analyses of the data and classification protocols to provide the best estimate of the true status of the detection. In this paper we conduct such analyses and show that Kepler-62f is a true planet and not a false positive. The table of stellar and planet properties has been updated based on GAIA results.
A rigorous definition of the habitable zone and its dependence on planetary properties is part of the search for habitable exoplanets. In this work, we use the general circulation model ExoCAM to determine how the inner edge of the habitable zone of tidally locked planets orbiting M dwarf stars depends on planetary radius, surface gravity, and surface pressure. We find that the inner edge of the habitable zone for more massive planets occurs at higher stellar irradiation, as found in previous 1D simulations. We also determine the relative effects of varying planetary radius and surface gravity. Increasing the planetary radius leads to a lower planetary albedo and warmer climate, pushing the inner edge of the habitable zone to lower stellar irradiation. This results from a change in circulation regime that leads to the disruption of the thick, reflective cloud deck around the substellar point. Increasing gravity increases the outgoing longwave radiation, which moves the inner edge of the habitable zone to higher stellar irradiation. This is because the column mass of water vapor decreases with increasing gravity, leading to a reduction in the greenhouse effect. The effect of gravity on the outgoing longwave radiation is stronger than the effect of radius on the planetary albedo, so that increasing gravity and radius together causes the inner edge of the habitable zone to move to higher stellar irradiation. Our results show that the inner edge of the habitable zone for more massive terrestrial planets occurs at a larger stellar irradiation.
Recent observations have found a valley in the size distribution of close-in super-Earths that is interpreted as a signpost that close-in super-Earths are mostly rocky in composition. However, new models predict that planetesimals should first form at the water ice line such that close-in planets are expected to have a significant water ice component. Here we investigate the water contents of super-Earths by studying the interplay between pebble accretion, planet migration and disc evolution. Planets compositions are determined by their position relative to different condensation fronts (ice lines) throughout their growth. Migration plays a key role. Assuming that planetesimals start at or exterior to the water ice line (r > rH2O), inward migration causes planets to leave the source region of icy pebbles and therefore to have lower final water contents than in discs with either outward migration or no migration. The water ice line itself moves inward as the disc evolves, and delivers water as it sweeps across planets that formed dry. The relative speed and direction of planet migration and inward drift of the water ice line is thus central in determining planets water contents. If planet formation starts at the water ice line, this implies that hot close-in super-Earths (r < 0.3 AU) with water contents of a few percent are a signpost of inward planet migration during the early gas phase. Hot super-Earths with larger water ice contents on the other hand, experienced outward migration at the water ice line and only migrated inwards after their formation was complete either because they become too massive to be contained in the region of outward migration or in chains of resonant planets. Measuring the water ice content of hot super-Earths may thus constrain their migration history.