Robustly modeling the inner edge of the habitable zone is essential for determining the most promising potentially habitable exoplanets for atmospheric characterization. Global climate models (GCMs) have become the standard tool for calculating this boundary, but divergent results have emerged among the various GCMs. In this study, we perform an intercomparison of standard GCMs used in the field on a rapidly rotating planet receiving a G-star spectral energy distribution and on a tidally locked planet receiving an M-star spectral energy distribution. Experiments both with and without clouds are examined. We find relatively small difference (within 8 K) in global-mean surface temperature simulation among the models in the G-star case with clouds. In contrast, the global-mean surface temperature simulation in the M-star case is highly divergent (2030 K). Moreover, even differences in the simulated surface temperature when clouds are turned off are significant. These differences are caused by differences in cloud simulation and/or radiative transfer, as well as complex interactions between atmospheric dynamics and these two processes. For example we find that an increase in atmospheric absorption of shortwave radiation can lead to higher relative humidity at high altitudes globally and, therefore, a significant decrease in planetary radiation emitted to space. This study emphasizes the importance of basing conclusions about planetary climate on simulations from a variety of GCMs and motivates the eventual comparison of GCM results with terrestrial exoplanet observations to improve their performance.
One-dimensional (vertical) models of planetary atmospheres typically balance the net solar and internal energy fluxes against the net thermal radiative and convective heat fluxes to determine an equilibrium thermal structure. Thus, simple models of shortwave and longwave (optical and infrared) radiative transport can provide insight into key processes operating within planetary atmospheres. Here, we develop a simple, analytic expression for both the downwelling thermal and net thermal radiative fluxes in a planetary troposphere. For thermal wavelengths, we assume that the atmosphere is non-scattering and that opacities are grey (i.e., wavelength-independent). Additionally, we adopt an atmospheric thermal structure that follows a modified dry adiabat as well as a physically-motivated power-law relationship between grey thermal optical depth and atmospheric pressure. To verify aspects of our analytic treatment, we compare our model to more sophisticated full physics tools as applied to Venus, Earth, and a cloudfree Jupiter, thereby exploring a diversity of atmospheric conditions. Next, we seek to better understand our analytic model by exploring how thermal radiative flux profiles respond to variations in key physical parameters, such as the total grey thermal optical depth of the atmosphere. Using energy balance arguments, we derive convective flux profiles for the tropospheres of all Solar System worlds with thick atmospheres, and propose a scaling that enables inter-comparison of these profiles. Lastly, we use our analytic treatment to discuss the validity of other simple models of convective fluxes in planetary atmospheres. Our new expressions build on decades of analytic modeling exercises in planetary atmospheres, and further prove the utility of simple, generalized tools in comparative planetology studies.
Some atmospheric gases have been proposed as counter indicators to the presence of life on an exoplanet if remotely detectable at sufficient abundance (i.e., antibiosignatures), informing the search for biosignatures and potentially fingerprinting uninhabited habitats. However, the quantitative extent to which putative antibiosignatures could exist in the atmospheres of inhabited planets is not well understood. The most commonly referenced potential antibiosignature is CO, because it represents a source of free energy and reduced carbon that is readily exploited by life on Earth and is thus often assumed to accumulate only in the absence of life. Yet, biospheres actively produce CO through biomass burning, photooxidation processes, and release of gases that are photochemically converted into CO in the atmosphere. We demonstrate with a 1D ecosphere-atmosphere model that reducing biospheres can maintain CO levels of ~100 ppmv even at low H2 fluxes due to the impact of hybrid photosynthetic ecosystems. Additionally, we show that photochemistry around M dwarf stars is particularly favorable for the buildup of CO, with plausible concentrations for inhabited, oxygen-rich planets extending from hundreds of ppm to several percent. Since CH4 buildup is also favored on these worlds, and because O2 and O3 are likely not detectable with the James Webb Space Telescope, the presence of high CO (>100 ppmv) may discriminate between oxygen-rich and reducing biospheres with near-future transmission observations. These results suggest that spectroscopic detection of CO can be compatible with the presence of life and that a comprehensive contextual assessment is required to validate the significance of potential antibiosignatures.
Earthshine is the dominant source of natural illumination on the surface of the Moon during lunar night, and at some locations within permanently shadowed regions (PSRs) near the poles that never receive direct sunlight. As such, earthshine has the potential to enable the scientific investigation and exploration of conditions in areas of the Moon that are either temporarily or permanently hidden from the Sun. Earthshine has also been used to refer to Earthlight reflected from the lunar surface, but in this study we use it to refer specifically to Earthlight incident at the Moon. Under certain circumstances, the heat flux from earthshine could also influence the transport and cold-trapping of volatiles present in the very coldest areas within PSRs. In this study, Earth’s spectral irradiance, as it would appear at the Moon in the solar reflectance band (0.3–3.0 µm) and at thermal emission wavelengths (3–50 µm), is examined with a suite of model image cubes and whole-disk spectra created using the Virtual Planetary Laboratory (VPL) three-dimensional (latitude, longitude and altitude) modeling capability
Giant planets migrate though the protoplanetary disc as they grow their solid core and attract their gaseous envelope. Previously, we have studied the growth and migration of an isolated planet in an evolving disc. Here, we generalise such models to include the mutual gravitational interaction between a high number of growing planetary bodies. We have investigated how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 2040 AU when using nominal type-I and type-II migration rates and a pebble flux of approximately 100200 Earth masses per million years, enough to grow Jupiter within the lifetime of the solar nebula. The planetary embryos placed up to 30 AU migrate into the inner system (rP < 1AU). There they form super-Earths or hot and warm gas giants, producing systems that are inconsistent with the configuration of the solar system, but consistent with some exoplanetary systems. We also explored slower migration rates which allow the formation of gas giants from embryos originating from the 510 AU region, which are stranded exterior to 1 AU at the end of the gas-disc phase. These giant planets can also form in discs with lower pebbles fluxes (50100 Earth masses per Myr). We identify a pebble flux threshold below which migration dominates and moves the planetary core to the inner disc, where the pebble isolation mass is too low for the planet to accrete gas efficiently. In our model, giant planet growth requires a sufficiently high pebble flux to enable growth to out-compete migration. An even higher pebble flux produces systems with multiple gas giants. We show that planetary embryos starting interior to 5 AU do not grow into gas giants, even if migration is slow and the pebble flux is large. These embryos instead grow to just a few Earth masses, the mass regime of super-Earths. This stunted growth is caused by the low pebble isolation mass in the inner disc and is therefore independent of the pebble flux. Additionally, we show that the long-term evolution of our formed planetary systems can naturally produce systems with inner super-Earths and outer gas giants as well as systems of giant planets on very eccentric orbits.
Oxygen and methane are considered to be the canonical biosignatures of modern Earth, and the simultaneous detection of these gases in a planetary atmosphere is an especially strong biosignature. However, these gases may be challenging to detect together in the planetary atmospheres because photochemical oxygen radicals destroy methane. Previous work has shown that the photochemical lifetime of methane in oxygenated atmospheres is longer around M dwarfs, but M dwarf planet habitability may be hindered by extreme stellar activity and evolution. Here, we use a 1D photochemical-climate model to show that K dwarf stars also offer a longer photochemical lifetime of methane in the presence of oxygen compared to G dwarfs. For example, we show that a planet orbiting a K6V star can support about an order of magnitude more methane in its atmosphere compared to an equivalent planet orbiting a G2V star. In the reflected-light spectra of worlds orbiting K dwarf stars, strong oxygen and methane features could be observed at visible and near-infrared wavelengths. Because K dwarfs are dimmer than G dwarfs, they offer a better planet-star contrast ratio, enhancing the signal-to-noise ratio (S/N) possible in a given observation. For instance, a 50 hr observation of a planet at 7 pc with a 15 m telescope yields S/N = 9.2 near 1 ?m for a planet orbiting a solar-type G2V star, and S/N = 20 for the same planet orbiting a K6V star. In particular, nearby mid-late K dwarfs such as 61 Cyg A/B, Epsilon Indi, Groombridge 1618, and HD 156026 may be excellent targets for future biosignature searches.
At least 30% of main sequence stars host planets with sizes between 1 and 4 Earth radii and orbital periods of less than 100 days. We use N-body simulations including a model for gas-assisted pebble accretion and disk-planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-Neptunes (?15M?) and super-massive planetary cores potentially able to become gas giant planets (?15M?). Simulations with moderate pebble fluxes grow multiple super-Earth-mass planets that migrate inwards and pile up at the disk’s inner edge forming long resonant chains. We follow the long-term dynamical evolution of these systems and use the period ratio distribution of observed planet-pairs to constrain our model. Up to ?95% of resonant chains become dynamically unstable after the gas disk dispersal, leading to a phase of late collisions that breaks the resonant configuration. Our simulations match observations if we combine a dominant fraction (?95%) of unstable systems with a sprinkling (?5%) of stable resonant chains (the Trappist-1 system represents one such example). Our results demonstrate that super-Earth systems are inherently multiple (N?2) and that the observed excess of single-planet transits is a consequence of the mutual inclinations excited by the planet-planet instability. In simulations in which planetary seeds are initially distributed in the inner and outer disk, close-in super-Earths are systematically ice-rich. This contrasts with the interpretation that most super-Earths are rocky based on bulk density measurements of super-Earths and photo-evaporation modeling of their bimodal radius distribution. We investigate the conditions needed to form rocky super-Earths. The formation of rocky super-Earths (abridged)
Knowing the high-energy radiation environment of a star over a planet’s formation and evolutionary period is critical in determining if that planet is potentially habitable and if any biosignatures could be detected, as UV radiation can severely change or destroy a planet’s atmosphere. Current efforts for finding a potentially habitable planet are focused on M stars, yet K stars may offer more habitable conditions due to decreased stellar activity and more distant and wider habitable zones (HZs). While M star activity evolution has been observed photometrically and spectroscopically, there has been no dedicated investigation of K star UV evolution. We present the first comprehensive study of the near-UV, far-UV, and X-ray evolution of K stars. We used members of young moving groups and clusters ranging in age from 10 to 625 Myr combined with field stars and their archived GALEX UV and ROSAT X-ray data to determine how the UV and X-ray radiation evolve. We find that the UV and X-ray flux incident on an HZ planet is 550 times lower than that of HZ planets around early-M stars and 501000 times lower than those around late-M stars, due to both an intrinsic decrease in K dwarf stellar activity occurring earlier than for M dwarfs and the more distant location of the K dwarf HZ.
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.
We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon, or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar transit light curves, planetary secondary eclipse light curves, and planetplanet/planetmoon occultation light curves, as well as thermal (rotational) phase curves. In this paper, we also introduce starry, an open-source package written in C++ and wrapped in Python that computes these light curves. The algorithm in starry is six orders of magnitude faster than direct numerical integration and several orders of magnitude more precise. starry also computes analytic derivatives of the light curves with respect to all input parameters for use in gradient-based optimization and inference, such as Hamiltonian Monte Carlo (HMC), allowing users to quickly and efficiently fit observed light curves to infer properties of a celestial body’s surface map. (Please see https://github.com/rodluger/starry, https://rodluger.github.io/starry/, and https://doi.org/10.5281/zenodo.1312286).
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.