The recent detections of temperate terrestrial planets orbiting nearby stars and the promise of characterizing their atmospheres motivate a need to understand how the diversity of possible planetary parameters affects the climate of terrestrial planets. In this work, we investigate the atmospheric circulation and climate of terrestrial exoplanets orbiting both Sun-like and M dwarf stars over a wide swath of possible planetary parameters, including the planetary rotation period, surface pressure, incident stellar flux, surface gravity, planetary radius, and cloud particle size. We do so using a general circulation model (GCM) that includes nongray radiative transfer and the effects of clouds. The results from this suite of simulations generally show qualitatively similar dependencies of circulation and climate on planetary parameters to idealized GCMs, with quantitative differences due to the inclusion of additional model physics. Notably, we find that the effective cloud particle size is a key unknown parameter that can greatly affect the climate of terrestrial exoplanets. We confirm a transition between low and high dayside cloud coverage of synchronously rotating terrestrial planets with increasing rotation period. We determine that this cloud transition is due to eddy-driven convergence near the substellar point and should not be parameterization dependent. Finally, we compute full-phase light curves from our simulations of planets orbiting M dwarf stars, finding that changing incident stellar flux and rotation period affect observable properties of terrestrial exoplanets. Our GCM results can guide expectations for planetary climate over the broad range of possible terrestrial exoplanets that will be observed with future space telescopes.
Containing only a few percentages of the mass of the moon, the current asteroid belt is around three to four orders of magnitude smaller than its primordial mass inferred from disk models. Yet dynamical studies have shown that the asteroid belt could not have been depleted by more than about an order of magnitude over the past ~4 Gyr. The remainder of the mass loss must have taken place during an earlier phase of the solar system’s evolution. An orbital instability in the outer solar system occurring during the process of terrestrial planet formation can reproduce the broad characteristics of the inner solar system. Here, we test the viability of this model within the constraints of the main belt’s low present-day mass and orbital structure. Although previous studies modeled asteroids as massless test particles because of limited computing power, our work uses graphics processing unit acceleration to model a fully self-gravitating asteroid belt. We find that depletion in the main belt is related to the giant planets’ exact evolution within the orbital instability. Simulations that produce the closest matches to the giant planets’ current orbits deplete the main belt by two to three orders of magnitude. These simulated asteroid belts are also good matches to the actual asteroid belt in terms of their radial mixing and broad orbital structure.
The solar systems dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations of the Early Instability scenario. We use an N-body integration scheme that accounts for collisional fragmentation, and also perform a large set of control simulations that do not include an early giant planet instability. Since the total particle number decreases slower when collisional fragmentation is accounted for, the growing planets orbits are damped more strongly via dynamical friction and encounters with small bodies that dissipate angular momentum (eg: hit-and-run impacts). Compared with simulations without collisional fragmentation, our fully evolved systems provide better matches to the solar systems terrestrial planets in terms of their compact mass distribution and dynamically cold orbits. Collisional processes also tend to lengthen the dynamical accretion timescales of Earth analogs, and shorten those of Mars analogs. This yields systems with relative growth timescales more consistent with those inferred from isotopic dating. Accounting for fragmentation is thus supremely important for any successful evolutionary model of the inner solar system.
Each of the giant planets within the Solar system has large moons but none of these moons have their own moons (which we call submoons). By analogy with studies of moons around short-period exoplanets, we investigate the tidal-dynamical stability of submoons. We find that 10?km-scale submoons can only survive around large (1000?km-scale) moons on wide-separation orbits. Tidal dissipation destabilizes the orbits of submoons around moons that are small or too close to their host planet; this is the case for most of the Solar systems moons. A handful of known moons are, however, capable of hosting long-lived submoons: Saturns moons Titan and Iapetus, Jupiters moon Callisto, and Earths Moon. Based on its inferred mass and orbital separation, the newly discovered exomoon candidate Kepler-1625b-I can in principle host a large submoon, although its stability depends on a number of unknown parameters. We discuss the possible habitability of submoons and the potential for subsubmoons. The existence, or lack thereof, of submoons may yield important constraints on satellite formation and evolution in planetary systems.
The TRAPPIST-1 planetary system provides an unprecedented opportunity to study terrestrial exoplanet evolution with the James Webb Space Telescope (JWST) and ground-based observatories. Since M dwarf planets likely experience extreme volatile loss, the TRAPPIST-1 planets may have highly evolved, possibly uninhabitable atmospheres. We used a versatile, 1D terrestrial planet climate model with line-by-line radiative transfer and mixing length convection (VPL Climate) coupled to a terrestrial photochemistry model to simulate environmental states for the TRAPPIST-1 planets. We present equilibrium climates with self-consistent atmospheric compositions and observational discriminants of postrunaway, desiccated, 10–100 bar O2- and CO2-dominated atmospheres, including interior outgassing, as well as for water-rich compositions
(The Astrophysical Journal, 2018)
CO2?driven changes to climate have occurred during many epochs of Earth’s history when the solar insolation, atmospheric CO2 concentration, and surface temperature of the planet were all significantly different than today. Each of these aspects affects the implied radiative forcings, climate feedbacks, and resultant changes in global mean surface temperature. Here we use a three?dimensional climate system model to study the effects of increasing CO2 on Earth’s climate, across many orders of magnitude of variation, and under solar inputs relevant for paleo, present, and future Earth scenarios. We find that the change in global mean surface temperature from doubling CO2 (i.e., the equilibrium climate sensitivity) may vary between 2.6 and 21.6 K over the course of Earth’s history. In agreement with previous studies, we find that the adjusted radiative forcing from doubling CO2 increases at high concentrations up to about 1.5 bars partial pressure, generally resulting in larger changes in the surface temperature. We also find that the cloud albedo feedback causes an abrupt transition in climate for warming atmospheres that depends both on the mean surface temperature and the total solar insolation. Climate sensitivity to atmospheric CO2 has probably varied considerably across Earth’s history.
We present empirical evidence, supported by a planet formation model, to show that the curve $R/{R}_{oplus }=1.05{(F/{F}_{oplus })}^{0.11}$ approximates the location of the so-called photo-evaporation valley. Planets below that curve are likely to have experienced complete photo-evaporation, and planets just above it appear to have inflated radii; thus we identify a new population of inflated super-Earths and mini-Neptunes. Our N-body simulations are set within an evolving protoplanetary disk and include prescriptions for orbital migration, gas accretion, and atmospheric loss due to giant impacts. Our simulated systems broadly match the sizes and periods of super-Earths in the Kepler catalog. They also reproduce the relative sizes of adjacent planets in the same system, with the exception of planet pairs that straddle the photo-evaporation valley. This latter group is populated by planet pairs with either very large or very small size ratios (R out/R in Gt 1 or R out/R in Lt 1) and a dearth of size ratios near unity. It appears that this feature could be reproduced if the planet outside the photo-evaporation valley (typically the outer planet, but sometimes not) has its atmosphere significantly expanded by stellar irradiation. This new population of planets may be ideal targets for future transit spectroscopy observations with the upcoming James Webb Space Telescope.
A leading model for the origin of super-Earths proposes that planetary embryos migrate inward and pile up on close-in orbits. As large embryos are thought to preferentially form beyond the snowline, this naively predicts that most super-Earths should be very water-rich. Here we show that the shortest period planets formed in the migration model are often purely rocky. The inward migration of icy embryos through the terrestrial zone accelerates the growth of rocky planets via resonant shepherding. We illustrate this process with a simulation that provided a match to the Kepler-36 system of two planets on close orbits with very different densities. In the simulation, two super-Earths formed in a Kepler-36-like configuration; the inner planet was pure rock while the outer one was ice-rich. We conclude from a suite of simulations that the feeding zones of close-in super-Earths are likely to be broad and disconnected from their final orbital radii.
The main asteroid belt (MB) is low in mass but dynamically excited. Here we propose a new mechanism to excite the MB during the giant planet (the “Nice model”) instability, which is expected to feature repeated close encounters between Jupiter and one or more ice giants (“jumping Jupiter” or JJ). We show that, when Jupiter temporarily reaches a high-enough level of excitation, both in eccentricity and inclination, it induces strong forced vectors of eccentricity and inclination across the MB region. Because during the JJ instability Jupiter’s orbit “jumps” around, the forced vectors keep changing both in magnitude and phase throughout the whole MB region. The entire cold primordial MB is thus excited as a natural outcome of the JJ instability. The level of such an excitation, however, is typically larger than the current orbital excitation observed in the MB. We show that the subsequent evolution of the solar system is capable of reshaping the resultant overexcited MB to its present-day orbital state, and that a strong mass depletion (~90%) is associated with the JJ instability phase and its subsequent evolution throughout the age of the solar system.
Aims. Evidence of mutually inclined planetary orbits has been reported for giant planets in recent years. Here we aim to study the impact of eccentric and inclined massive giant planets on the terrestrial planet formation process, and investigate whether it can possibly lead to the formation of inclined terrestrial planets.
Methods. We performed 126 simulations of the late-stage planetary accretion in eccentric and inclined giant planet systems. The physical and orbital parameters of the giant planet systems result from n-body simulations of three giant planets in the late stage of the gas disc, under the combined action of Type II migration and planet-planet scattering. Fourteen two- and three-planet configurations were selected, with diversified masses, semi-major axes (resonant configurations or not), eccentricities, and inclinations (including coplanar systems) at the dispersal of the gas disc. We then followed the gravitational interactions of these systems with an inner disc of planetesimals and embryos (nine runs per system), studying in detail the final configurations of the formed terrestrial planets.
Results. In addition to the well-known secular and resonant interactions between the giant planets and the outer part of the disc, giant planets on inclined orbits also strongly excite the planetesimals and embryos in the inner part of the disc through the combined action of nodal resonance and the LidovKozai mechanism. This has deep consequences on the formation of terrestrial planets. While coplanar giant systems harbour several terrestrial planets, generally as massive as the Earth and mainly on low-eccentric and low-inclined orbits, terrestrial planets formed in systems with mutually inclined giant planets are usually fewer, less massive (<0.5 M?), and with higher eccentricities and inclinations. This work shows that terrestrial planets can form on stable inclined orbits through the classical accretion theory, even in coplanar giant planet systems emerging from the disc phase.
Solar and thermal radiation are critical aspects of planetary climate, with gradients in radiative energy fluxes driving heating and cooling. Climate models require that radiative transfer tools be versatile, computationally efficient, and accurate. Here, we describe a technique that uses an accurate full-physics radiative transfer model to generate a set of atmospheric radiative quantities which can be used to linearly adapt radiative flux profiles to changes in the atmospheric and surface state—the Linearized Flux Evolution (LiFE) approach. These radiative quantities describe how each model layer in a plane-parallel atmosphere reflects and transmits light, as well as how the layer generates diffuse radiation by thermal emission and by scattering light from the direct solar beam. By computing derivatives of these layer radiative properties with respect to dynamic elements of the atmospheric state, we can then efficiently adapt the flux profiles computed by the full-physics model to new atmospheric states. We validate the LiFE approach, and then apply this approach to Mars, Earth, and Venus, demonstrating the information contained in the layer radiative properties and their derivatives, as well as how the LiFE approach can be used to determine the thermal structure of radiative and radiative-convective equilibrium states in one-dimensional atmospheric models.
Many dynamical aspects of the solar system can be explained by the outer planets experiencing a period of orbital instability sometimes called the Nice Model. Though often correlated with a perceived delayed spike in the lunar cratering record known as the Late Heavy Bombardment (LHB), recent work suggests that this event may have occurred much earlier; perhaps during the epoch of terrestrial planet formation. While current simulations of terrestrial accretion can reproduce many observed qualities of the solar system, replicating the small mass of Mars requires modification to standard planet formation models. Here we use 800 dynamical simulations to show that an early instability in the outer solar system strongly influences terrestrial planet formation and regularly yields properly sized Mars analogs. Our most successful outcomes occur when the terrestrial planets evolve an additional 110 million years (Myr) following the dispersal of the gas disk, before the onset of the giant planet instability. In these simulations, accretion has begun in the Mars region before the instability, but the dynamical perturbation induced by the giant planets scattering removes large embryos from Mars vicinity. Large embryos are either ejected or scattered inward toward Earth and Venus (in some cases to deliver water), and Mars is left behind as a stranded embryo. An early giant planet instability can thus replicate both the inner and outer solar system in a single model.
VPL Scientist Dr. Shawn Domagal-Goldman (NASA GSFC) was recently featured on the Ask an Astrobiologist program! Dr. Domagal-Goldman is a…
Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star’s habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here, we use 1-D coupled climate-photochemical models to generate self-consistent atmospheres for several evolutionary scenarios, including high-O2, high-CO2, and more Earth-like atmospheres, with both oxic and anoxic compositions. We show that these modeled environments can be habitable or uninhabitable at Proxima Cen b’s position in the habitable zone. We use radiative transfer models to generate synthetic spectra and thermal phase curves for these simulated environments, and use instrument models to explore our ability to discriminate between possible planetary states.
VPL Scientist Dr. Sean Raymond (Université de Bordeaux) was recently featured on the Ask an Astrobiologist program! Dr. Raymond’s research…
The obliquity of the Earth, which controls our seasons, varies by only ~2fdg5 over ~40,000 years, and its eccentricity varies by only ~0.05 over 100,000 years. Nonetheless, these small variations influence Earth’s ice ages. For exoplanets, however, variations can be significantly larger. Previous studies of the habitability of moonless Earth-like exoplanets have found that high obliquities, high eccentricities, and dynamical variations can extend the outer edge of the habitable zone by preventing runaway glaciation (snowball states). We expand upon these studies by exploring the orbital dynamics with a semianalytic model that allows us to map broad regions of parameter space. We find that, in general, the largest drivers of obliquity variations are secular spin–orbit resonances. We show how the obliquity varies in several test cases, including Kepler-62 f, across a wide range of orbital and spin parameters. These obliquity variations, alongside orbital variations, will have a dramatic impact on the climates of such planets.
The TRAPPIST-1, Proxima Centauri, and LHS 1140 systems are the most exciting prospects for future follow-up observations of potentially inhabited planets. All of the planets orbit nearby M-stars and are likely tidally locked in 1:1 spinorbit states, which motivates the consideration of the effects that tidal locking might have on planetary habitability. On Earth, periods of global glaciation (snowballs) may have been essential for habitability and remote signs of life (biosignatures) because they are correlated with increases in the complexity of life and in the atmospheric oxygen concentration. In this paper, we investigate the snowball bifurcation (sudden onset of global glaciation) on tidally locked planets using both an energy balance model and an intermediate-complexity global climate model. We show that tidally locked planets are unlikely to exhibit a snowball bifurcation as a direct result of the spatial pattern of insolation they receive. Instead, they will smoothly transition from partial to complete ice coverage and back. A major implication of this work is that tidally locked planets with an active carbon cycle should not be found in a snowball state. Moreover, this work implies that tidally locked planets near the outer edge of the habitable zone with low CO2 outgassing fluxes will equilibrate with a small unglaciated substellar region rather than cycling between warm and snowball states. More work is needed to determine how the lack of a snowball bifurcation might affect the development of life on a tidally locked planet.
Read More >>
In his article “Getting Real about the Oxygen Biosignature”, Marc Kaufman discusses Dr. Victoria Meadows‘ most recent paper “Reflections on O2…
VPL researchers Jacob Haqq-Misra, Ravi Kopparapu, and Eric Wolf are recent published authors of the International Journal of Astrobiology. In…
David Crisp was awarded an NASA Exceptional Public Service Medal for contributions to radiative transfer advances used across planetary and earth…
As lower-mass stars often host multiple rocky planets, gravitational interactions among planets can have significant effects on climate and habitability over long timescales. Here we explore a specific case, Kepler-62f, a potentially habitable planet in a five-planet system with a K2V host star. N-body integrations reveal the stable range of initial eccentricities for Kepler-62f is 0.00⩽e⩽0.32, absent the effect of additional, undetected planets. We simulate the tidal evolution of Kepler-62f in this range and find that, for certain assumptions, the planet can be locked in a synchronous rotation state. Simulations using LMD Generic GCM indicate that with 3 bars of CO2 in its atmosphere, Kepler-62f would only be warm enough for surface liquid water at the upper limit of this eccentricity range, providing it has a high planetary obliquity (between 60∘ and 90∘). A climate similar to modern-day Earth is possible for the entire range of stable eccentricities if atmospheric CO2 is increased to 5-bar levels. In a low-CO2 case, simulations with CCSM4 and LMD Generic GCM indicate that increases in planetary obliquity and orbital eccentricity coupled with an orbital configuration that places the summer solstice at or near pericenter permit regions of the planet with above-freezing surface temperatures. This may melt ice sheets formed during colder seasons. If Kepler-62f is synchronously rotating and has an ocean, CO2 levels above 3 bars would be required to distribute enough heat to the night side of the planet to avoid atmospheric freeze-out and permit a large enough region of open water at the planet’s substellar point to remain stable. Overall, we find multiple plausible combinations of orbital and atmospheric properties that permit surface liquid water on Kepler-62f.
Jupiter’s core is generally assumed to have formed beyond the snow line. Here we consider an alternative scenario, that Jupiter’s core may have accumulated in the innermost parts of the protoplanetary disk. A growing body of research suggests that small particles (“pebbles”) continually drift inward through the disk. If a fraction of drifting pebbles is trapped at the inner edge of the disk a several Earth-mass core can quickly grow. Subsequently, the core may migrate outward beyond the snow line via planet-disk interactions. Of course, to reach the outer Solar System Jupiter’s core must traverse the terrestrial planet-forming region.
Reproducing Uranus and Neptune remains a challenge for simulations of solar system formation. The ice giants’ peculiar obliquities suggest that they both suffered giant collisions during their formation. Thus, there must have been an epoch of accretion dominated by collisions among large planetary embryos in the primordial outer solar system. We test this idea using N-body numerical simulations including the effects of a gaseous protoplanetary disk. One strong constraint is that the masses of the ice giants are very similar – the Neptune and Uranus mass ratio is ~1.18. We show that similar-sized ice giants do indeed form by collisions between planetary embryos beyond Saturn. The fraction of successful simulations varies depending on the initial number of planetary embryos in the system, their individual and total masses. Similar-sized ice giants are consistently reproduced in simulations starting with five to ten planetary embryos with initial masses of ~3–6 M⊕. We conclude that accretion from a population of planetary embryos is a plausible scenario for the origin of Uranus and Neptune.