In the search for life on Earth-like planets around other stars, the first (and likely only) information will come from the spectroscopic characterization of the planet’s atmosphere. Of the countless number of chemical species terrestrial life produces, only a few have the distinct spectral features and the necessary atmospheric abundance to be detectable. The easiest of these species to observe in Earth’s atmosphere is O2 (and its photochemical byproduct, O3).
Atmospheric collapse is likely to be of fundamental importance to tidally locked rocky exoplanets but remains understudied. Here, general results on the heat transport and stability of tidally locked terrestrial-type atmospheres are reported. First, the problem is modeled with an idealized 3D general circulation model (GCM) with gray gas radiative transfer. It is shown that over a wide range of parameters the atmospheric boundary layer, rather than the large-scale circulation, is the key to understanding the planetary energy balance. Through a scaling analysis of the interhemispheric energy transfer, theoretical expressions for the day-night temperature difference and surface wind speed are created that reproduce the GCM results without tuning. Next, the GCM is used with correlated-k radiative transfer to study heat transport for two real gases (CO2 and CO). For CO2, empirical formulae for the collapse pressure as a function of planetary mass and stellar flux are produced, and critical pressures for atmospheric collapse at Earth’s stellar flux are obtained that are around five times higher (0.14 bar) than previous gray gas estimates. These results provide constraints on atmospheric stability that will aid in future interpretation of observations and exoplanet habitability modeling.
We use a 3‐D general circulation model to compare the primitive Martian hydrological cycle in “warm and wet” and “cold and icy” scenarios. In the warm and wet scenario, an anomalously high solar flux or intense greenhouse warming artificially added to the climate model are required to maintain warm conditions and an ice‐free northern ocean. Precipitation shows strong surface variations, with high rates around Hellas basin and west of Tharsis but low rates around Margaritifer Sinus (where the observed valley network drainage density is nonetheless high). In the cold and icy scenario, snow migration is a function of both obliquity and surface pressure, and limited episodic melting is possible through combinations of seasonal, volcanic, and impact forcing. At surface pressures above those required to avoid atmospheric collapse (∼0.5 bar) and moderate to high obliquity, snow is transported to the equatorial highland regions where the concentration of valley networks is highest. Snow accumulation in the Aeolis quadrangle is high, indicating an ice‐free northern ocean is not required to supply water to Gale crater. At lower surface pressures and obliquities, both H2O and CO2 are trapped as ice at the poles and the equatorial regions become extremely dry. The valley network distribution is positively correlated with snow accumulation produced by the cold and icy simulation at 41.8∘ obliquity but uncorrelated with precipitation produced by the warm and wet simulation. Because our simulations make specific predictions for precipitation patterns under different climate scenarios, they motivate future targeted geological studies.
Most of the properties of the Earth–Moon system can be explained by a collision between a planetary embryo (giant impactor) and the growing Earth late in the accretion process1,2,3. Simulations show that most of the material that eventually aggregates to form the Moon originates from the impactor1,4,5. However, analysis of the terrestrial and lunar isotopic compositions show them to be highly similar6,7,8,9,10,11. In contrast, the compositions of other Solar System bodies are significantly different from those of the Earth and Moon12,13,14, suggesting that different Solar System bodies have distinct compositions. This challenges the giant impact scenario, because the Moon-forming impactor must then also be thought to have a composition different from that of the proto-Earth. Here we track the feeding zones of growing planets in a suite of simulations of planetary accretion15, to measure the composition of Moon-forming impactors. We find that different planets formed in the same simulation have distinct compositions, but the compositions of giant impactors are statistically more similar to the planets they impact. A large fraction of planet–impactor pairs have almost identical compositions. Thus, the similarity in composition between the Earth and Moon could be a natural consequence of a late giant impact.
We show that terrestrial planets in the habitable zones of M dwarfs older than ∼1 Gyr could have been in runaway greenhouses for several hundred million years following their formation due to the star’s extended pre-main sequence phase, provided they form with abundant surface water. Such prolonged runaway greenhouses can lead to planetary evolution divergent from that of Earth. During this early runaway phase, photolysis of water vapor and hydrogen/oxygen escape to space can lead to the loss of several Earth oceans of water from planets throughout the habitable zone, regardless of whether the escape is energy-limited or diffusion-limited.
What fascinates people about astrobiology is that it seeks answers to long-standing unsolved questions: How quickly did life evolve on Earth and why did life persist here? Is there life elsewhere in the Solar System or beyond? Astrobiology: A Very Short Introduction explores some of the big unanswered questions about the universe, considers the origins of life on Earth and its evolution, and brings together the ideas of microbiologists, astronomers, planetary scientists, and geologists. It introduces the origins of astrobiology and demonstrates its impact on current astronomical research and…
Earth-like planets within the liquid water habitable zone of M-type stars may evolve into synchronous rotators. On these planets, the substellar hemisphere experiences perpetual daylight while the opposing antistellar hemisphere experiences perpetual darkness. Because the night-side hemisphere has no direct source of energy, the air over this side of the planet is prone to freeze out and deposit on the surface, which could result in atmospheric collapse. However, general circulation models (GCMs) have shown that atmospheric dynamics can counteract this problem and provide sufficient energy transport to the antistellar side. Here, we use an idealized GCM to consider the impact of geothermal heating on the habitability of synchronously rotating planets.
Embedded in the gaseous protoplanetary disk, Jupiter and Saturn naturally become trapped in 3:2 resonance and migrate outward. This serves as the basis of the Grand Tack model. However, previous hydrodynamical simulations were restricted to isothermal disks, with moderate aspect ratio and viscosity. Here we simulate the orbital evolution of the gas giants in disks with viscous heating and radiative cooling. We find that Jupiter and Saturn migrate outward in 3:2 resonance in modest-mass (M disk ≈ M MMSN, where MMSN is the “minimum-mass solar nebula”) disks with viscous stress parameter α between 10–3 and 10–2. In disks with relatively low-mass (M disk lesssim M MMSN), Jupiter and Saturn get captured in 2:1 resonance and can even migrate outward in low-viscosity disks (α ≤ 10–4). Such disks have a very small aspect ratio (h ~ 0.02-0.03) that favors outward migration after capture in 2:1 resonance, as confirmed by isothermal runs which resulted in a similar outcome for h ~ 0.02 and α ≤ 10–4. We also performed N-body runs of the outer solar system starting from the results of our hydrodynamical simulations and including 2-3 ice giants. After dispersal of the gaseous disk, a Nice model instability starting with Jupiter and Saturn in 2:1 resonance results in good solar systems analogs. We conclude that in a cold solar nebula, the 2:1 resonance between Jupiter and Saturn can lead to outward migration of the system, and this may represent an alternative scenario for the evolution of the solar system.
National Aeronautics and Space Administration’s (NASA) Kepler Space Telescope has detected over 3,000 planet candidates, about a dozen of which are probably rocky planets within the liquid-water habitable zones of their parent stars. Climate-modeling calculations discussed here shed light on the width of that zone. Within the next several years, NASA may obtain spectra of nearby transiting Earth-sized planets around M stars, using its James Webb Space Telescope. NASA hopes to build an even more capable space telescope to perform direct imaging of Earth-like exoplanets and take spectra of their atmospheres. Once data are obtained from either of these missions, correct interpretation of possible biomarker gases will become critical. We discuss here how those interpretations might be made.
We explore the impact of obliquity variations on planetary habitability in hypothetical systems with high mutual inclination. We show that large-amplitude, high-frequency obliquity oscillations on Earth-like exoplanets can suppress the ice-albedo feedback, increasing the outer edge of the habitable zone. We restricted our exploration to hypothetical systems consisting of a solar-mass star, an Earth-mass planet at 1 AU, and 1 or 2 larger planets. We verified that these systems are stable for 108 years with N-body simulations and calculated the obliquity variations induced by the orbital evolution of the Earth-mass planet and a torque from the host star. We ran a simplified energy balance model on the terrestrial planet to assess surface temperature and ice coverage on the planet’s surface, and we calculated differences in the outer edge of the habitable zone for planets with rapid obliquity variations. For each hypothetical system, we calculated the outer edge of habitability for two conditions: (1) the full evolution of the planetary spin and orbit and (2) the eccentricity and obliquity fixed at their average values. We recovered previous results that higher values of fixed obliquity and eccentricity expand the habitable zone, but we also found that obliquity oscillations further expand habitable orbits in all cases. Terrestrial planets near the outer edge of the habitable zone may be more likely to support life in systems that induce rapid obliquity oscillations as opposed to fixed-spin planets. Such planets may be the easiest to directly characterize with space-borne telescopes. Key Words: Exoplanets—Habitable zone—Energy balance models. Astrobiology 14, 277–291.
A minimum atmospheric temperature, or tropopause, occurs at a pressure of around 0.1 bar in the atmospheres of Earth1, Titan2, Jupiter3, Saturn4, Uranus and Neptune4, despite great differences in atmospheric composition, gravity, internal heat and sunlight. In all of these bodies, the tropopause separates a stratosphere with a temperature profile that is controlled by the absorption of short-wave solar radiation, from a region below characterized by convection, weather and clouds5,6. However, it is not obvious why the tropopause occurs at the specific pressure near 0.1 bar. Here we use a simple, physically based model7 to demonstrate that, at atmospheric pressures lower than 0.1 bar, transparency to thermal radiation allows short-wave heating to dominate, creating a stratosphere. At higher pressures, atmospheres become opaque to thermal radiation, causing temperatures to increase with depth and convection to ensue. A common dependence of infrared opacity on pressure, arising from the shared physics of molecular absorption, sets the 0.1 bar tropopause. We reason that a tropopause at a pressure of approximately 0.1 bar is characteristic of many thick atmospheres, including exoplanets and exomoons in our galaxy and beyond. Judicious use of this rule could help constrain the atmospheric structure, and thus the surface environments and habitability, of exoplanets.
The presence of valleys on ancient terrains of Mars suggests that liquid water flowed on the martian surface 3.8 Gyr ago or before. The above-freezing temperatures required to explain valley formation could have been transient, in response to the frequent large meteorite impacts on early Mars, or they could have been caused by long-lived greenhouse warming. Climate models that consider only the greenhouse gases carbon dioxide and water have been unable to recreate warm surface conditions, given the lower solar luminosity at that time. Here we use a one-dimensional climate model to demonstrate that an atmosphere containing 1.3–4 bar of CO2 and water, in addition to 5–20% H2, could have raised the mean surface temperature of early Mars above the freezing point of water. Vigorous volcanic outgassing from a highly reduced early martian mantle is expected to provide sufficient atmospheric H2 and CO2—the latter from the photochemical oxidation of outgassed CH4 and CO—to form a CO2 and H2 greenhouse. Such a dense early martian atmosphere is consistent with independent estimates of surface pressure based on cratering data.
The Great Oxidation Event (GOE) was an increase in atmospheric oxygen levels from less than 1 ppm to 0.2–2% by volume 2.4–2.2 billion years ago. In atmospheric chemistry, hydrogen-bearing reduced gases, such as methane and hydrogen, are adversaries of O2. When the concentration of hydrogen-bearing reduced gases goes up, O2 declines, and vice versa. Thus, the pre-GOE atmospheric redox chemistry should have been dominated by methane and hydrogen. Before the GOE, oxygen was driven to trace levels by reactions with volcanic and metamorphic reductants, including dissolved cations (e.g., Fe2 +) in surface waters and reducing gases (H2, CH4, CO, SO2, and H2S). Rapid escape of hydrogen to space from such an atmosphere would have slowly oxidized the Earth. A ‘tipping point’ was reached when the flux of O2 associated with the burial of organic carbon exceeded O2 losses. Oxidative weathering then became significant. Models suggest that methane level fell before the GOE and such loss of greenhouse gases plausibly caused global cooling. Multiple glaciations during 2.45–2.22 Ga hint that the climate and atmospheric composition oscillated until permanent oxygenation was established. Subsequent levels of O2 were sufficient to protect the Earth’s surface from harmful ultraviolet with a stratospheric ozone layer.
Observations by the SPICAV/SOIR instruments aboard Venus Express have revealed that the upper haze (UH) of Venus, between 70 and 90 km, is variable on the order of days and that it is populated by two particle modes. We use a one-dimensional microphysics and vertical transport model based on the Community Aerosol and Radiation Model for Atmospheres to evaluate whether interaction of upwelled cloud particles and sulfuric acid particles nucleated in situ on meteoric dust are able to generate the two observed modes, and whether their observed variability are due in part to the action of vertical transient winds at the cloud tops. Nucleation of photochemically produced sulfuric acid onto polysulfur condensation nuclei generates mode 1 cloud droplets, which then diffuse upwards into the UH. Droplets generated in the UH from nucleation of sulfuric acid onto meteoric dust coagulate with the upwelled cloud particles and therefore cannot reproduce the observed bimodal size distribution. By comparison, the mass transport enabled by transient winds at the cloud tops, possibly caused by sustained subsolar cloud top convection, are able to generate a bimodal size distribution in a time scale consistent with Venus Express observations.
We study the evolution of multiple protoplanets of a few Earth masses embedded in a non-isothermal protoplanetary disc. The protoplanets are located in the vicinity of a convergence zone located at the transition between two different opacity regimes. Inside the convergence zone, Type I migration is directed outward and outside the zone migration is directed inward.
In a multiplanet system, tides acting on the inner planet can significantly affect the orbital evolution of the entire system. While tides usually damp eccentricities, a novel mechanism identified by Correia et al. (Astrophys J Lett 744, article id. L23, 2012) tends to raise eccentricities as a result of the tides’ effect on the inner planet’s rotation. Our analytical description of this spin-driven tidal (SDT) effect shows that, while the inner planet’s eccentricity undergoes pumping, the process is more completely described by an exchange of strength between the two eigenmodes of the dynamical system. Our analysis allows derivation of criteria for two-planet coplanar systems where the SDT effect can reverse tidal damping, and may preclude the effect’s being significant for realistic systems. For the specific case quantified by Correia et al., the effect is strong because of the large adopted tidal time lag, which may not be appropriate for the assumed Saturn-like inner planet. On the other hand, the effective
The field of exoplanetary science has experienced a recent surge of new systems that is largely due to the precision photometry provided by the Kepler mission. The latest discoveries have included compact planetary systems in which the orbits of the planets all lie relatively close to the host star, which presents interesting challenges in terms of formation and dynamical evolution. The compact exoplanetary systems are analogous to the moons orbiting the giant planets in our solar system, in terms of their relative sizes and semimajor axes. We present a study that quantifies the scaled sizes and separations of the solar system moons with respect to their hosts. We perform a similar study for a large sample of confirmed Kepler planets in multi-planet systems. We show that a comparison between the two samples leads to a similar correlation between their scaled sizes and separation distributions. The different gradients of the correlations may be indicative of differences in the formation and/or long-term dynamics of moon and planetary systems.
The potential habitability of newly discovered exoplanets is initially assessed by determining whether their orbits fall within the circumstellar habitable zone of their star. However, the habitable zone (HZ) is not static in time or space, and its boundaries migrate outward at a rate proportional to the increase in luminosity of a star undergoing stellar evolution, possibly including or excluding planets over the course of the star’s main sequence lifetime. We describe the time that a planet spends within the HZ as its “habitable zone lifetime.” The HZ lifetime of a planet has strong astrobiological implications and is especially important when considering the evolution of complex life, which is likely to require a longer residence time within the HZ. Here, we present results from a simple model built to investigate the evolution of the “classic” HZ over time, while also providing estimates for the evolution of stellar luminosity over time in order to develop a “hybrid” HZ model. These models return estimates for the HZ lifetimes of Earth and 7 confirmed HZ exoplanets and 27 unconfirmed Kepler candidates. The HZ lifetime for Earth ranges between 6.29 and 7.79×109 years (Gyr). The 7 exoplanets fall in a range between ∼1 and 54.72 Gyr, while the 27 Kepler candidate planets’ HZ lifetimes range between 0.43 and 18.8 Gyr. Our results show that exoplanet HD 85512b is no longer within the HZ, assuming it has an Earth analog atmosphere. The HZ lifetime should be considered in future models of planetary habitability as setting an upper limit on the lifetime of any potential exoplanetary biosphere, and also for identifying planets of high astrobiological potential for continued observational or modeling campaigns. Key Words: Exoplanet habitability metrics—Continuously habitable zone—Stellar evolution—Planetary habitability. Astrobiology 13, 833–849.
Planetary climate can be affected by the interaction of the host star spectral energy distribution with the wavelength-dependent reflectivity of ice and snow. In this study, we explored this effect with a one-dimensional (1-D), line-by-line, radiative transfer model to calculate broadband planetary albedos as input to a seasonally varying, 1-D energy balance climate model. A three-dimensional (3-D) general circulation model was also used to explore the atmosphere’s response to changes in incoming stellar radiation, or instellation, and surface albedo. Using this hierarchy of models, we simulated planets covered by ocean, land, and water-ice of varying grain size, with incident radiation from stars of different spectral types. Terrestrial planets orbiting stars with higher near-UV radiation exhibited a stronger ice-albedo feedback. We found that ice extent was much greater on a planet orbiting an F-dwarf star than on a planet orbiting a G-dwarf star at an equivalent flux distance, and that ice-covered conditions occurred on an F-dwarf planet with only a 2% reduction in instellation relative to the present instellation on Earth, assuming fixed CO2 (present atmospheric level on Earth).
If mutual gravitational scattering among exoplanets occurs, then it may produce unique orbital properties. For example, two-planet systems that lie near the boundary between circulation and libration of their periapses could result if planet-planet scattering ejected a former third planet quickly, leaving one planet on an eccentric orbit and the other on a circular orbit. We first improve upon previous work that examined the apsidal behavior of known multiplanet systems by doubling the sample size and including observational uncertainties. This analysis recovers previous results that demonstrated that many systems lay on the apsidal boundary between libration and circulation. We then performed over 12,000 three-dimensional N-body simulations of hypothetical three-body systems that are unstable, but stabilize to two-body systems after an ejection. Using these synthetic two-planet systems, we test the planet-planet scattering hypothesis by comparing their apsidal behavior, over a range of viewing angles, to that of the observed systems and find that they are statistically consistent regardless of the multiplicity of the observed systems. Finally, we combine our results with previous studies to show that, from the sampled cases, the most likely planetary mass function prior to planet-planet scattering follows a power law with index -1.1. We find that this pre-scattering mass function predicts a mutual inclination frequency distribution that follows an exponential function with an index between -0.06 and -0.1.
The atmospheres of terrestrial planets are expected to be in long-term radiation balance: an increase in the absorption of solar radiation warms the surface and troposphere, which leads to a matching increase in the emission of thermal radiation. Warming a wet planet such as Earth would make the atmosphere moist and optically thick such that only thermal radiation emitted from the upper troposphere can escape to space. Hence, for a hot moist atmosphere, there is an upper limit on the thermal emission that is unrelated to surface temperature. If the solar radiation absorbed exceeds this limit, the planet will heat uncontrollably and the entire ocean will evaporate—the so-called runaway greenhouse. Here we model the solar and thermal radiative transfer in incipient and complete runaway greenhouse atmospheres at line-by-line spectral resolution using a modern spectral database. We find a thermal radiation limit of 282 W m−2 (lower than previously reported) and that 294 W m−2 of solar radiation is absorbed (higher than previously reported). Therefore, a steam atmosphere induced by such a runaway greenhouse may be a stable state for a planet receiving a similar amount of solar radiation as Earth today. Avoiding a runaway greenhouse on Earth requires that the atmosphere is subsaturated with water, and that the albedo effect of clouds exceeds their greenhouse effect. A runaway greenhouse could in theory be triggered by increased greenhouse forcing, but anthropogenic emissions are probably insufficient.
Excess emission, associated with warm, dust belts, commonly known as exozodis, has been observed around a third of nearby stars. The high levels of dust required to explain the observations are not generally consistent with steady-state evolution. A common suggestion is that the dust results from the aftermath of a dynamical instability, an event akin to the Solar system’s Late Heavy Bombardment. In this work, we use a data base of N-body simulations to investigate the aftermath of dynamical instabilities between giant planets in systems with outer planetesimal belts. We find that, whilst there is a significant increase in the mass of material scattered into the inner regions of the planetary system following an instability, this is a short-lived effect. Using the maximum lifetime of this material, we determine that even if every star has a planetary system that goes unstable, there is a very low probability that we observe more than a maximum of 1 per cent of sun-like stars in the aftermath of an instability, and that the fraction of planetary systems currently in the aftermath of an instability is more likely to be limited to ≤0.06 per cent. This probability increases marginally for younger or higher mass stars. We conclude that the production of warm dust in the aftermath of dynamical instabilities is too short-lived to be the dominant source of the abundantly observed exozodiacal dust.
Given their tendency to be incorporated into the core during differentiation, the highly-siderophile elements (HSEs) in Earth’s mantle are thought to have been accreted as a “late veneer” after the end of the giant impact phase. Bottke et al. (Bottke, W.F., Walker, R.J., Day, J.M.D., Nesvorny, D., Elkins-Tanton, L. [2010]. Science 330, 1527) proposed that the large Earth-to-Moon HSE abundance ratio can be explained if the late veneer was characterized by large (D = 1000–4000 km) impactors. Here we simulate the evolution of the terrestrial planets during a stochastic late veneer phase from the end of accretion until the start of the late heavy bombardment ∼500 Myr later. We show that a late veneer population of 0.05M⊕ dominated by large (D > 1000 km) bodies naturally delivers a ∼0.01M⊕ veneer to Earth, consistent with geochemical constraints.