On Earth, methane is produced mainly by life, and it has been proposed that, under certain conditions, methane detected in an exoplanetary spectrum may be considered a biosignature. Here, we estimate how much methane may be produced in hydrothermal vent systems by serpentinization, its main geological source, using the kinetic properties of the main reactions involved in methane production by serpentinization. Hydrogen production by serpentinization was calculated as a function of the available FeO in the crust, given the current spreading rates. Carbon dioxide is the limiting reactant for methane formation because it is highly depleted in aqueous form in hydrothermal vent systems. We estimated maximum CH4 surface fluxes of 6.8×108 and 1.3×109 molecules cm−2 s−1 for rocky planets with 1 and 5 M⊕, respectively. Using a 1-D photochemical model, we simulated atmospheres with volume mixing ratios of 0.03 and 0.1 CO2 to calculate atmospheric methane concentrations for the maximum production of this compound by serpentinization. The resulting abundances were 2.5 and 2.1 ppmv for 1 M⊕ planets and 4.1 and 3.7 ppmv for 5 M⊕ planets. Therefore, low atmospheric concentrations of methane may be produced by serpentinization. For habitable planets around Sun-like stars with N2-CO2 atmospheres, methane concentrations larger than 10 ppmv may indicate the presence of life. Key Words: Serpentinization—Exoplanets—Biosignatures—Planetary atmospheres. Astrobiology 13, 550–559.
We study the influence of outer solar system architecture on the structural evolution of the Oort Cloud (OC) and the flux of Earth-crossing comets. In particular, we seek to quantify the role of the giant planets as “planetary protectors.” To do so, we have run simulations in each of four different planetary mass configurations to understand the significance of each of the giant planets. Because the outer planets modify the structure of the OC throughout its formation, we integrate each simulation over the full age of the solar system. Over this time, we follow the evolution of cometary orbits from their starting point in the protoplanetary disk to their injection into the OC to their possible re-entry into the inner planetary region. We find that the overall structure of the OC, including the location of boundaries and the relative number of comets in the inner and outer parts, does not change significantly between configurations; however, as planetary mass decreases, the trapping efficiency (TE) of comets into the OC and the flux of comets into the observable region increases. We determine that those comets that evolve onto Earth-crossing orbits come primarily from the inner OC but show no preference for initial protoplanetary disk location. We also find that systems that have at least a Saturn-mass object are effective at deflecting possible Earth-crossing comets but the difference in flux between systems with and without such a planet is less than an order of magnitude. We conclude by discussing the individual roles of the planets and the implications of incorporating more realistic planetary accretion and migration scenarios into simulations, particularly on existing discrepancies between low TE and the mass of the protoplanetary disk and on determining the structural boundaries of the OC.
The Kepler-11 planetary system contains six transiting planets ranging in size from 1.8 to 4.2 times the radius of Earth. Five of these planets orbit in a tightly packed configuration with periods between 10 and 47 days. We perform a dynamical analysis of the system based upon transit timing variations observed in more than three years of Kepler photometric data. Stellar parameters are derived using a combination of spectral classification and constraints on the star’s density derived from transit profiles together with planetary eccentricity vectors provided by our dynamical study. Combining masses of the planets relative to the star from our dynamical study and radii of the planets relative to the star from transit depths together with deduced stellar properties yields measurements of the radii of all six planets, masses of the five inner planets, and an upper bound to the mass of the outermost planet, whose orbital period is 118 days. We find mass-radius combinations for all six planets that imply that substantial fractions of their volumes are occupied by constituents that are less dense than rock. Moreover, we examine the stability of these envelopes against photoevaporation and find that the compositions of at least the inner two planets have likely been significantly sculpted by mass loss. The Kepler-11 system contains the lowest mass exoplanets for which both mass and radius have been measured.
The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed 10 transits of XO-2b over a period of 3 yr. We present measurements that confirm previous estimates of system parameters like the normalized semi-major axis (a/R sstarf), stellar density (ρsstarf), impact parameter (b), and orbital inclination (i orb). Our errors on system parameters like a/R sstarf and ρsstarf have improved by ~40% compared to previous best ground-based measurements. Our study of the transit times show no evidence for transit timing variations (TTVs) and we are able to rule out co-planar companions with masses ≥0.20 M ⊕ in low order mean motion resonance with XO-2b. We also explored the stability of the XO-2 system given various orbital configurations of a hypothetical planet near the 2:1 mean motion resonance. We find that a wide range of orbits (including Earth-mass perturbers) are both dynamically stable and produce observable TTVs. We find that up to 51% of our stable simulations show TTVs that are smaller than the typical transit timing errors (~20 s) measured for XO-2b, and hence remain undetectable.
Traditionally, stellar radiation has been the only heat source considered capable of determining global climate on long timescales. Here, we show that terrestrial exoplanets orbiting low-mass stars may be tidally heated at high-enough levels to induce a runaway greenhouse for a long-enough duration for all the hydrogen to escape. Without hydrogen, the planet no longer has water and cannot support life. We call these planets “Tidal Venuses” and the phenomenon a “tidal greenhouse.” Tidal effects also circularize the orbit, which decreases tidal heating. Hence, some planets may form with large eccentricity, with its accompanying large tidal heating, and lose their water, but eventually settle into nearly circular orbits (i.e., with negligible tidal heating) in the habitable zone (HZ).
Traditionally, stellar radiation has been the only heat source considered capable of determining global climate on long timescales. Here, we show that terrestrial exoplanets orbiting low-mass stars may be tidally heated at high-enough levels to induce a runaway greenhouse for a long-enough duration for all the hydrogen to escape. Without hydrogen, the planet no longer has water and cannot support life. We call these planets “Tidal Venuses” and the phenomenon a “tidal greenhouse.” Tidal effects also circularize the orbit, which decreases tidal heating. Hence, some planets may form with large eccentricity, with its accompanying large tidal heating, and lose their water, but eventually settle into nearly circular orbits (i.e., with negligible tidal heating) in the habitable zone (HZ).
Because of their large numbers, low-mass stars may be the most abundant planet hosts in our Galaxy. Furthermore, terrestrial planets in the habitable zones (HZs) around M-dwarfs can potentially be characterized in the near future and hence may be the first such planets to be studied. Recently, Dressing & Charbonneau used Kepler data and calculated the frequency of terrestrial planets in the HZ of cool stars to be $0.15^{+0.13}_{-0.06}$ per star for Earth-size planets (0.5-1.4 R ⊕). However, this estimate was derived using the Kasting et al. HZ limits, which were not valid for stars with effective temperatures lower than 3700 K. Here we update their result using new HZ limits from Kopparapu et al. for stars with effective temperatures between 2600 K and 7200 K, which includes the cool M stars in the Kepler target list.
White and brown dwarfs are astrophysical objects that are bright enough to support an insolation habitable zone (IHZ). Unlike hydrogen-burning stars, they cool and become less luminous with time; hence their IHZ moves in with time. The inner edge of the IHZ is defined as the orbital radius at which a planet may enter a moist or runaway greenhouse, phenomena that can remove a planet’s surface water forever. Thus, as the IHZ moves in, planets that enter it may no longer have any water and are still uninhabitable. Additionally, the close proximity of the IHZ to the primary leads to concern that tidal heating may also be strong enough to trigger a runaway greenhouse, even for orbital eccentricities as small as 10−6. Water loss occurs due to photolyzation by UV photons in the planetary stratosphere, followed by hydrogen escape. Young white dwarfs emit a large amount of these photons, as their surface temperatures are over 104 K. The situation is less clear for brown dwarfs, as observational data do not constrain their early activity and UV emission very well. Nonetheless, both types of planets are at risk of never achieving habitable conditions, but planets orbiting white dwarfs may be less likely to sustain life than those orbiting brown dwarfs. We consider the future habitability of the planet candidates KOI 55.01 and 55.02 in these terms and find they are unlikely to become habitable. Key Words: Extrasolar terrestrial planets—Habitability—Habitable zone—Tides—Exoplanets. Astrobiology 13, 279–291.
Identifying terrestrial planets in the habitable zones (HZs) of other stars is one of the primary goals of ongoing radial velocity (RV) and transit exoplanet surveys and proposed future space missions. Most current estimates of the boundaries of the HZ are based on one-dimensional (1D), cloud-free, climate model calculations by Kasting et al. However, this model used band models that were based on older HITRAN and HITEMP line-by-line databases. The inner edge of the HZ in the Kasting et al. model was determined by loss of water, and the outer edge was determined by the maximum greenhouse provided by a CO2 atmosphere. A conservative estimate for the width of the HZ from this model in our solar system is 0.95-1.67 AU. Here an updated 1D radiative-convective, cloud-free climate model is used to obtain new estimates for HZ widths around F, G, K, and M stars. New H2O and CO2 absorption coefficients, derived from the HITRAN 2008 and HITEMP 2010 line-by-line databases, are important improvements to the climate model. According to the new model, the water-loss (inner HZ) and maximum greenhouse (outer HZ) limits for our solar system are at 0.99 and 1.70 AU, respectively, suggesting that the present Earth lies near the inner edge. Additional calculations are performed for stars with effective temperatures between 2600 and 7200 K, and the results are presented in parametric form, making them easy to apply to actual stars. The new model indicates that, near the inner edge of the HZ, there is no clear distinction between runaway greenhouse and water-loss limits for stars with T eff lesssim 5000 K, which has implications for ongoing planet searches around K and M stars. To assess the potential habitability of extrasolar terrestrial planets, we propose using stellar flux incident on a planet rather than equilibrium temperature. This removes the dependence on planetary (Bond) albedo, which varies depending on the host star’s spectral type. We suggest that conservative estimates of the HZ (water-loss and maximum greenhouse limits) should be used for current RV surveys and Kepler mission to obtain a lower limit on η⊕, so that future flagship missions like TPF-C and Darwin are not undersized. Our model does not include the radiative effects of clouds; thus, the actual HZ boundaries may extend further in both directions than the estimates just given.
Earth-mass planets embedded in gaseous protoplanetary disks undergo Type I orbital migration. In radiative disks an additional component of the corotation torque scaling with the entropy gradient across the horseshoe region can counteract the general inward migration, Type I migration can then be directed either inward or outward depending on the local disk properties. Thus, special locations exist in the disk toward which planets migrate in a convergent way. Here we present N-body simulations of the convergent migration of systems of low-mass (M = 1−10M⊕) planets. We show that planets do not actually converge in convergence zones.
The detection of moons orbiting extrasolar planets (“exomoons”) has now become feasible. Once they are discovered in the circumstellar habitable zone, questions about their habitability will emerge. Exomoons are likely to be tidally locked to their planet and hence experience days much shorter than their orbital period around the star and have seasons, all of which works in favor of habitability. These satellites can receive more illumination per area than their host planets, as the planet reflects stellar light and emits thermal photons. On the contrary, eclipses can significantly alter local climates on exomoons by reducing stellar illumination. In addition to radiative heating, tidal heating can be very large on exomoons, possibly even large enough for sterilization. We identify combinations of physical and orbital parameters for which radiative and tidal heating are strong enough to trigger a runaway greenhouse. By analogy with the circumstellar habitable zone, these constraints define a circumplanetary “habitable edge.” We apply our model to hypothetical moons around the recently discovered exoplanet Kepler-22b and the giant planet candidate KOI211.01 and describe, for the first time, the orbits of habitable exomoons. If either planet hosted a satellite at a distance greater than 10 planetary radii, then this could indicate the presence of a habitable moon. Key Words: Astrobiology—Extrasolar planets—Habitability—Habitable zone—Tides. Astrobiology 13, 18–46.
Nearly half the exoplanets found within binary star systems reside1 in very wide binaries with average stellar separations greater than 1,000 astronomical units (one astronomical unit (AU) being the Earth–Sun distance), yet the influence of such distant binary companions on planetary evolution remains largely unstudied. Unlike their tighter counterparts, the stellar orbits of wide binaries continually change under the influence of the Milky Way’s tidal field and impulses from other passing stars. Here we report numerical simulations demonstrating that the variable nature of wide binary star orbits dramatically reshapes the planetary systems they host, typically billions of years after formation. Contrary to previous understanding2, wide binary companions may often strongly perturb planetary systems, triggering planetary ejections and increasing the orbital eccentricities of surviving planets. Although hitherto not recognized, orbits of giant exoplanets within wide binaries are statistically more eccentric than those around isolated stars. Both eccentricity distributions are well reproduced when we assume that isolated stars and wide binaries host similar planetary systems whose outermost giant planets are scattered beyond about 10 AU from their parent stars by early internal instabilities. Consequently, our results suggest that although wide binaries eventually remove the most distant planets from many planetary systems, most isolated giant exoplanet systems harbour additional distant, still undetected planets.
It is commonly assumed that the building blocks of the terrestrial planets were derived from a cosmochemical reservoir that is best represented by chondrites, the so-called chondritic Earth model. This view is possibly a good approximation for refractory elements (although it has been recently questioned; e.g., Caro et al. 2008), but for volatile elements, other cosmochemical reservoirs might have contributed to Earth, such as the solar nebula gas and/or cometary matter (Owen et al. 1992; Dauphas 2003; Pepin 2006). Hence, in order to get insights into the origin of the carbon in Earth, it is…
Starting from planetary systems with three giant planets and an outer disc of planetesimals, we use dynamical simulations to show how dynamical instabilities can transform planetesimal discs into 102–103 au-scale isotropic clouds. The instabilities involve a phase of planet–planet scattering that concludes with the ejection of one or more planets and the inward-scattering of the surviving gas giant(s) to remove them from direct dynamical contact with the planetesimals. ‘Mini-Oort clouds’ are thus formed from scattered planetesimals whose orbits are frozen by the abrupt disappearance of the perturbing giant planet. Although the planetesimal orbits are virtually isotropic, the surviving giant planets tend to have modest inclinations (typically ∼10°) with respect to the initial orbital plane. The collisional lifetimes of mini-Oort clouds are long (10 Myr to >10 Gyr) and there is a window of ∼100 Myr or longer during which they produce spherical clouds of potentially observable dust at 70 μm. If the formation channel for hot Jupiters commonly involves planetary close encounters, we predict a correlation between this subset of extrasolar planetary systems and mini-Oort clouds.
We present an analytic one-dimensional radiative-convective model of the thermal structure of planetary atmospheres. Our model assumes that thermal radiative transfer is gray and can be represented by the two-stream approximation. Model atmospheres are assumed to be in hydrostatic equilibrium, with a power-law scaling between the atmospheric pressure and the gray thermal optical depth. The convective portions of our models are taken to follow adiabats that account for condensation of volatiles through a scaling parameter to the dry adiabat. By combining these assumptions, we produce simple, analytic expressions that allow calculations of the atmospheric-pressure-temperature profile, as well as expressions for the profiles of thermal radiative flux and convective flux. We explore the general behaviors of our model. These investigations encompass (1) worlds where atmospheric attenuation of sunlight is weak, which we show tend to have relatively high radiative-convective boundaries; (2) worlds with some attenuation of sunlight throughout the atmosphere, which we show can produce either shallow or deep radiative-convective boundaries, depending on the strength of sunlight attenuation; and (3) strongly irradiated giant planets (including hot Jupiters), where we explore the conditions under which these worlds acquire detached convective regions in their mid-tropospheres. Finally, we validate our model and demonstrate its utility through comparisons to the average observed thermal structure of Venus, Jupiter, and Titan, and by comparing computed flux profiles to more complex models.
We present models for the formation of terrestrial planets, and the collisional evolution of debris disks, in planetary systems that contain multiple marginally unstable gas giants. We previously showed that in such systems, the dynamics of the giant planets introduces a correlation between the presence of terrestrial planets and cold dust, i.e., debris disks, which is particularly pronounced at λ ~ 70 μm. Here we present new simulations that show that this connection is qualitatively robust to a range of parameters: the mass distribution of the giant planets, the width and mass distribution of the outer planetesimal disk, and the presence of gas in the disk when the giant planets become unstable. We discuss how variations in these parameters affect the evolution. We find that systems with equal-mass giant planets undergo the most violent instabilities, and that these destroy both terrestrial planets and the outer planetesimal disks that produce debris disks.
Earth’s climate and biosphere have always shaped one another. James F. Kasting approves of an attempt to reveal the planet’s future by reading its past.
[1] The S‐wave velocity in the shallow subsurface within seismically active regions self‐organizes so that typical strong dynamic shear stresses marginally exceed the Coulomb elastic limit. The dynamic velocity from major strike‐slip faults yields simple dimensional relations. The near‐field velocity pulse is essentially a Love wave. The dynamic shear strain is the ratio of the measured particle velocity over the deep S‐wave velocity. The shallow dynamic shear stress is this quantity times the local shear modulus. The dynamic shear traction on fault parallel vertical planes is finite at the free surface. Coulomb failure occurs on favorably oriented fractures and internally in intact rock. I obtain the equilibrium shear modulus by starting a sequence of earthquakes with intact stiff rock extending all the way to the surface. The imposed dynamic shear strain in stiff rock causes Coulomb failure at shallow depths and leaves cracks in it wake. Cracked rock is more compliant than the original intact rock. Cracked rock is also weaker in friction, but shear modulus changes have a larger effect. Each subsequent event causes additional shallow cracking until the rock becomes compliant enough that it just reaches Coulomb failure over a shallow depth range of tens to hundreds of meters. Further events maintain the material at the shear modulus as a function where it just fails. The formalism provided in the paper yields reasonable representation of the S‐wave velocity in exhumed sediments near Cajon Pass and the San Fernando Valley of California. A general conclusion is that shallow rocks in seismically active areas just become nonlinear during typical shaking. This process causes transient changes in S‐wave velocity, but not strong nonlinear attenuation of seismic waves. Wave amplitudes significantly larger than typical ones would strongly attenuate and strongly damage the rock.
Submarine hydrothermal vents above serpentinite produce chemical potential gradients of aqueous and ionic hydrogen, thus providing a very attractive venue for the origin of life. This environment was most favourable before Earth’s massive CO2 atmosphere was subducted into the mantle, which occurred tens to approximately 100 Myr after the moon-forming impact; thermophile to clement conditions persisted for several million years while atmospheric pCO2 dropped from approximately 25 bar to below 1 bar. The ocean was weakly acid (pH ∼ 6), and a large pH gradient existed for nascent life with pH 9–11 fluids venting from serpentinite on the seafloor. Total CO2 in water was significant so the vent environment was not carbon limited. Biologically important phosphate and Fe(II) were somewhat soluble during this period, which occurred well before the earliest record of preserved surface rocks approximately 3.8 billion years ago (Ga) when photosynthetic life teemed on the Earth and the oceanic pH was the modern value of approximately 8. Serpentinite existed by 3.9 Ga, but older rocks that might retain evidence of its presence have not been found. Earth’s sequesters extensive evidence of Archaean and younger subducted biological material, but has yet to be exploited for the Hadean record.
Main sequence M stars are the most abundant stars in our galaxy, they have features that make them an attractive target for Astrobiology, but their strong chromospheric activity produces high energy radiation and charged particles that may be detrimental to life. We study the impact of a strong flare from the M dwarf, AD Leo, on the atmospheric chemistry of a hypothetical Earth-like planet located in the habitable zone. The simulations were performed using a 1-D photochemical model. We simulated
We used a standard equilibrium tidal model to compute the orbital evolution of a large ensemble of planet-BD systems. We tested the effect of numerous parameters such as the initial semi-major axis and eccentricity, the rotation period of the BD, the masses of both the BD and planet, and the tidal dissipation factors.
Recent numerical simulations have demonstrated that the Sun’s dynamical history within the Milky Way may be much more complex than that suggested by its current low peculiar velocity (Sellwood, J.A., Binney, J.J. [2002]. Mon. Not. R. Astron. Soc. 336, 785–796; Roškar, R., Debattista, V.P., Quinn, T.R., Stinson, G.S., Wadsley, J. [2008]. Astrophys. J. 684, L79–L82). In particular, the Sun may have radially migrated through the galactic disk by up to 5–6 kpc (Roškar, R., Debattista, V.P., Quinn, T.R., Stinson, G.S., Wadsley, J. [2008]. Astrophys. J. 684, L79–L82). This has important ramifications for the structure of the Oort Cloud, as it means that the Solar System may have experienced tidal and stellar perturbations that were significantly different from its current local galactic environment. To characterize the effects of solar migration within the Milky Way, we use direct numerical simulations to model the formation of an Oort Cloud around stars that end up on solar-type orbits in a galactic-scale simulation of a Milky Way-like disk formation.
The temperature profiles of micrometeors entering the atmospheres of Earth-like planets are calculated to determine the altitude at which exogenous organic compounds may be released. Previous experiments have shown that flash-heated micrometeorite analogs release organic compounds at temperatures from roughly 500 to 1000 K [1]. The altitude of release is of great importance because it determines the fate of the compound. Organic compounds that are released deeper in the atmosphere are more likely to rapidly mix to lower altitudes where they can accumulate to higher abundances or form more complex molecules and/or aerosols. Variables that are explored here are particle size, entry angle, atmospheric density profiles, spectral type of the parent star, and planet mass. The problem reduces to these questions: (1) How much atmosphere does the particle pass through by the time it is heated to 500 K? (2) Is the atmosphere above sufficient to attenuate stellar UV such that the mixing timescale is shorter than the photochemical timescale for a particular compound? We present preliminary results that the effect of the planetary and particle parameters have on the altitude of organic release
Most discussion of habitable planets has focused on Earth-like planets with globally abundant liquid water. For an “aqua planet” like Earth, the surface freezes if far from its sun, and the water vapor greenhouse effect runs away if too close. Here we show that “land planets” (desert worlds with limited surface water) have wider habitable zones than aqua planets. For planets at the inner edge of the habitable zone, a land planet has two advantages over an aqua planet: (i) the tropics can emit longwave radiation at rates above the traditional runaway limit because the air is unsaturated and (ii) the dry air creates a dry stratosphere that limits hydrogen escape. At the outer limits of the habitable zone, the land planet better resists global freezing because there is less water for clouds, snow, and ice. Here we describe a series of numerical experiments using a simple three-dimensional global climate model for Earth-sized planets.