Search

Sub-Neptunes (planets with radii between 2 and 4 R ⊕) are abundant around M-dwarf stars, yet the atmospheric dynamics of these planets is relatively unexplored. In this paper, we aim to provide a basic underpinning of the dry dynamics of general low-mean-molecular-weight, temperate sub-Neptune atmospheres. We use the ExoFMS general circulation model (GCM) with an idealized gray-gas radiation scheme to simulate planetary atmospheres with different levels of instellation and rotation rates, using the atmosphere of K2-18b as our control. We find that the atmospheres of tidally locked (TL), temperate sub-Neptunes have weak horizontal temperature gradients owing to their slow rotation rates and hydrogen-dominated composition. The zonal wind structure is dominated by high-latitude cyclostrophic jets driven by the conservation of angular momentum. At low pressures we observe superrotating equatorial jets, which we propose are driven by a Rossbyâ€"Kelvin instability similar to the type seen in simulations of idealized atmospheres with axisymmetric forcing. By viewing the flow in TL coordinates, we find the predominant overturning circulation to be between the day side and night side, and we derive scaling relations linking the TL stream function and vertical velocities to instellation. Comparing our results to the only other GCM study of K2-18b, we find significant qualitative differences in dynamics, highlighting the need for further collaboration and investigation into the effects of different dynamical cores and physical parameterizations. This paper provides a baseline for studying the dry dynamics of temperate sub-Neptunes, which will be built on in part II with the introduction of moist effects. [ABSTRACT FROM AUTHOR]

Titan’s stratosphere has been observed in a superrotation state, where the atmosphere rotates many times faster than the surface does. Another characteristic of Titan’s atmosphere is the presence of a thick haze layer. In this paper, we performed numerical experiments using a general circulation model to explore the effects of the haze layer on the stratospheric superrotation. We employed a semigray radiation model of Titan’s atmosphere following McKay et al., which takes account of sunlight absorption by the haze particles. The phase change of methane or seasonal changes were not taken into account. Our model with radiation parameters tuned for Titan yielded a global eastward wind around the equator with larger velocities at higher altitudes, except at around 70 km, after 105 Earth days. Although the atmosphere is not in an equilibrium state, the zonal wind profiles are approximately consistent with the observed one. By changing the parameters of the radiation model, we found that the intensity and the location of the maximum zonal wind velocity highly depended on the optical thickness and the altitude of the haze layer, respectively. Analysis of our experiments suggests that the quasi-stationary stratospheric superrotation is maintained by the balance between the meridional circulation decoupled from the surface and the eddies that transport angular momentum equatorward. This is different from, but similar to, the so-called Gierasch mechanism, in which momentum is supplied from the surface. This structure may explain the no-wind region at about 80 km in altitude. [ABSTRACT FROM AUTHOR]

Tidally locked terrestrial planets around low-mass stars are the prime targets of finding potentially habitable exoplanets. Several atmospheric general circulation models have been employed to simulate their possible climates; however, model intercomparisons showed that there are large differences in the results of the models even when they are forced with the same boundary conditions. In this paper, we examine whether model resolution contributes to the differences. Using the atmospheric general circulation model ExoCAM coupled to a 50 m slab ocean, we examine three different horizontal resolutions (440 km × 550 km, 210 km × 280 km, and 50 km × 70 km in latitude and longitude) and three different vertical resolutions (26, 51, and 74 levels) under the same dynamical core and the same schemes of radiation, convection, and clouds. Among the experiments, the differences are within 5 K in global-mean surface temperature and within 0.007 in planetary albedo. These differences are from cloud feedback, water vapor feedback, and the decreasing trend of relative humidity with increasing resolution. Relatively small-scale downdrafts between upwelling columns over the substellar region are better resolved and the mixing between dry and wet air parcels and between anvil clouds and their environment are enhanced as the resolution is increased. These reduce atmospheric relative humidity and high-level cloud fraction, causing a lower clear-sky greenhouse effect, a weaker cloud longwave radiation effect, and subsequently a cooler climate with increasing model resolution. Overall, the sensitivity of the simulated climate of tidally locked aquaplanets to model resolution is small. [ABSTRACT FROM AUTHOR]

The runaway greenhouse represents the ultimate climate catastrophe for rocky, Earth-like worlds: when the incoming stellar flux cannot be balanced by radiation to space, the oceans evaporate and exacerbate heating, turning the planet into a hot wasteland with a steam atmosphere overlying a possibly molten magma surface. The equilibrium state beyond the runaway greenhouse instellation limit depends on the radiative properties of the atmosphere and its temperature structure. Here, we use 1D radiative-convective models of steam atmospheres to explore the transition from the tropospheric radiation limit to the post-runaway climate state. To facilitate eventual simulations with 3D global circulation models, a computationally efficient band-gray model is developed, which is capable of reproducing the key features of the more comprehensive calculations. We analyze two factors that determine the equilibrated surface temperature of post-runaway planets. The infrared cooling of the planet is strongly enhanced by the penetration of the dry adiabat into the optically thin upper regions of the atmosphere. In addition, thermal emission of both shortwave and near-IR fluxes from the hot lower atmospheric layers, which can radiate through window regions of the spectrum, is quantified. Astronomical surveys of rocky exoplanets in the runaway greenhouse state may discriminate these features using multiwavelength observations. [ABSTRACT FROM AUTHOR]

The observation of radio, X-ray, and Hα emission from substellar objects indicates the presence of plasma regions and associated high-energy processes in their surrounding envelopes. This paper numerically simulates and characterizes critical velocity ionization (CVI), a potential ionization process, that can efficiently generate plasma as a result of neutral gas flows interacting with seed magnetized plasmas. By coupling a gas–magnetohydrodynamic (MHD) interactions code (to simulate the ionization mechanism) with a substellar global circulation model (to provide the required gas flows), we quantify the spatial extent of the resulting plasma regions, their degree of ionization, and their lifetime for a typical substellar atmosphere. It is found that the typical average ionization fraction reached at equilibrium (where the ionization and recombination rates are equal and opposite) ranges from 10−5 to 10−8, at pressures between 10−1 and 10−3 bar, with a trend of increasing ionization fraction with decreasing atmospheric pressure. The ionization fractions reached as a result of CVI are sufficient to allow magnetic fields to couple to gas flows in the atmosphere. [ABSTRACT FROM AUTHOR]

Atmospheric retrievals of exoplanet transmission spectra allow constraints on the composition and structure of the dayâ€"night terminator region. Such retrievals in the past have typically assumed one-dimensional (1D) temperature structures which were adequate to explain extant observations. However, the increasing data quality expected from exoplanet spectroscopy with the James Webb Space Telescope (JWST) motivates considerations of multidimensional atmospheric retrievals. We present Aura-3D, a three-dimensional atmospheric retrieval framework for exoplanet transmission spectra. Aura-3D includes a forward model that enables rapid computation of transmission spectra in 3D geometry for a given atmospheric structure and can, therefore, be used for atmospheric retrievals as well as for computing spectra from general circulation models (GCMs). In order to efficiently explore the space of possible 3D temperature structures in retrievals, we develop a parametric 3D pressureâ€"temperature profile which can accurately represent azimuthally averaged temperature structures of a range of hot Jupiter GCMs. We apply our retrieval framework to simulated JWST observations of hot Jupiter transmission spectra, obtaining accurate estimates of the dayâ€"night temperature variation across the terminator as well as the abundances of chemical species. We demonstrate an example of a model hot Jupiter transmission spectrum for which a traditional 1D retrieval of JWST-quality data returns biased abundance estimates, whereas a retrieval including a dayâ€"night temperature gradient can accurately retrieve the true abundances. Our forward model also has the capability to include inhomogeneous chemistry as well as variable clouds/hazes. This new retrieval framework opens the field to detailed multidimensional atmospheric characterization using transmission spectra of exoplanets in the JWST era. [ABSTRACT FROM AUTHOR]

Cloud is critical for planetary climate and habitability, but it is also one of the most challenging aspects of studying planets in and beyond the solar system. Here we use a cloud-resolving model (CRM) with high resolution (2 km) in a 2D configuration to simulate the clouds and circulation on tidally locked aquaplanets. We find that the substellar area is covered by deep convective clouds, the nightside is dominated by low-level clouds, and the two are linked by a global-scale Walker circulation. We further find that uniform surface warming causes the substellar cloud width to decrease, but a reduction in the dayâ€"night surface temperature contrast or an increase in the longwave radiative cooling rate causes the substellar cloud width to increase. These relationships can be roughly interpreted in accordance with simple thermodynamic theories. Comparing the results between the CRM and the global 3D general circulation model (GCM), we find that they are qualitatively consistent, including the Walker circulation, the substellar clouds, and the responses of the substellar ascending area and strength to changes in the surface temperature or in its zonal contrast. But large quantitative differences exist, such as the magnitude of the cloud water path, the cloud width, and their responses to external forcings. These results increase our confidence in using GCMs to model exoplanetary climates, although large quantitative uncertainties are always likely to exist. Future work is required to use 3D CRMs with realistic radiative transfer and the Coriolis force to examine the clouds and climates of tidally locked planets. [ABSTRACT FROM AUTHOR]

Close-in giant planets with strong stellar irradiation show atmospheric circulation patterns with strong equatorial jets and global-scale stationary waves. So far, almost all modeling works on atmospheric circulations of such giant planets have mainly considered external radiation alone, without taking into account the role of internal heat fluxes or just treating it in very simplified ways. Here, we study atmospheric circulations of strongly irradiated giant planets by considering the effect of internal forcing, which is characterized by small-scale stochastic interior thermal perturbations, using a three-dimensional atmospheric general circulation model. We show that the perturbation-excited waves can largely modify atmospheric circulation patterns in the presence of relatively strong internal forcing. Specifically, our simulations demonstrate three circulation regimes: a superrotation regime, a midlatitude-jet regime, and a quasi-periodic oscillation regime, depending on the relative importance of external and internal forcings. It is also found that strong internal forcing can cause noticeable modifications of the thermal phase curves. [ABSTRACT FROM AUTHOR]

Optical secondary eclipse measurements made by Kepler reveal a diverse set of geometric albedos for hot Jupiters with equilibrium temperatures between 1550 and 1700 K. The presence or absence of high-altitude condensates, such as Mg2SiO4, Fe, Al2O3, and TiO2, can significantly alter optical albedos, but these clouds are expected to be confined to localized regions in the atmospheres of these tidally locked planets. Here, we present 3D general circulation models and corresponding cloud and albedo maps for six hot Jupiters with measured optical albedos in this temperature range. We find that the observed optical albedos of K2-31b and K2-107b are best matched by either cloud-free models or models with relatively compact cloud layers, while Kepler-8b's and Kepler-17b's optical albedos can be matched by moderately extended (fsed = 0.1) parametric cloud models. HATS-11b has a high optical albedo, corresponding to models with bright Mg2SiO4 clouds extending to very low pressures (fsed = 0.03). We are unable to reproduce Kepler-7b's high albedo, as our models predict that the dayside will be dominated by dark Al2O3 clouds at most longitudes. We compare our parametric cloud model with a microphysical cloud model. We find that even after accounting for the 3D thermal structure, no single cloud model can explain the full range of observed albedos within the sample. We conclude that a better knowledge of the vertical mixing profiles, cloud radiative feedback, cloud condensate properties, and atmospheric metallicities is needed in order to explain the unexpected diversity of albedos in this temperature range. [ABSTRACT FROM AUTHOR]

Observations of scattered light and thermal emission from hot Jupiter exoplanets have suggested the presence of inhomogeneous aerosols in their atmospheres. 3D general circulation models (GCMs) that attempt to model the effects of aerosols have been developed to understand the physical processes that underlie their dynamical structures. In this work, we investigate how different approaches to aerosol modeling in GCMs of hot Jupiters affect high-resolution thermal emission spectra throughout the duration of the planet's orbit. Using results from a GCM with temperature-dependent cloud formation, we calculate spectra of a representative hot Jupiter with different assumptions regarding the vertical extent and thickness of clouds. We then compare these spectra to models in which clouds are absent or simply post-processed (i.e., added subsequently to the completed clear model). We show that the temperature-dependent treatment of clouds in the GCM produces high-resolution emission spectra that are markedly different from the clear and post-processed cases—both in the continuum flux levels and line profiles—and that increasing the vertical extent and thickness of clouds leads to bigger changes in these features. We evaluate the net Doppler shifts of the spectra induced by global winds and the planet's rotation and show that they are strongly phase dependent, especially for models with thicker and more extended clouds. This work further demonstrates the importance of radiative feedback in cloudy atmospheric models of hot Jupiters, as this can have a significant impact on interpreting spectroscopic observations of exoplanet atmospheres. [ABSTRACT FROM AUTHOR]

Using a general circulation model (GCM), we investigate trends in simulated hot Jupiter atmospheres for a range of irradiation temperatures (1500–4000 K), surface gravities (10 and 40 m s−2), and cloud conditions. Our models include simplified temperature-dependent clouds with radiative feedback and show how different cloud compositions, vertical thicknesses, and opacities shape hot Jupiter atmospheres by potentially increasing planetary albedos, decreasing photospheric pressures and nightside temperatures, and in some cases producing strong dayside thermal inversions. With decreasing irradiation, clouds progressively form on the nightside and cooler western limb, followed by the eastern limb and central dayside. We find that clouds significantly modify the radiative transport and affect the observable properties of planets colder than Tirr ≈ 3000 K (Teq ≈ 2100 K) depending on the clouds' vertical extent. The precise strength of expected effects depends on the assumed parameters, but trends in predicted phase curves emerge from an ensemble of simulations. Clouds lead to larger phase-curve amplitudes and smaller phase-curve offsets at IR wavelengths, compared to cloud-free models. At optical wavelengths, we predict mostly westward phase-curve offsets at intermediate temperatures (Tirr ≈ 2000–3500 K) with clouds confined to the nightside and western limb. If clouds are vertically compact (i.e., on the order of a pressure scale height in thickness), their distributions and effects become more complicated as different condensates form at different heights—some too deep to significantly affect the observable atmosphere. Our results have implications for interpreting the diversity of phase-curve observations of planets with Tirr ≲ 3000 K. [ABSTRACT FROM AUTHOR]

Brown dwarfs were recently found to display rotational modulations, commonly attributed to cloud cover of varying thickness, possibly modulated by planetary-scale waves. However, the long-term, continuous, high-precision monitoring data to test this hypothesis for more objects is lacking. By applying our novel photometric approach to Transiting Exoplanet Survey Satellite data, we extract a high-precision lightcurve of the closest brown dwarfs, which form the binary system Luhman 16 AB. Our observations, which cover about 100 rotations of Luhman 16 B, display continuous lightcurve evolution. The periodogram analysis shows that the rotational period of the component that dominates the lightcurve is 5.28 hr. We also find evidence for periods of 2.5, 6.94, and 90.8 hr. We show that the 2.5 and 5.28 hr periods emerge from Luhman 16 B and that they consist of multiple, slightly shifted peaks, revealing the presence of high-speed jets and zonal circulation in this object. We find that the lightcurve evolution is well fit by the planetary-scale waves model, further supporting this interpretation. We argue that the 6.94 hr peak is likely the rotation period of Luhman 16 A. By comparing the rotational periods to observed v sin(i) measurements, we show that the two brown dwarfs are viewed at angles close to their equatorial planes. We also describe a long-period (P ∼ 91 hr) evolution in the lightcurve, which we propose emerges from the vortex-dominated polar regions. Our study paves the way toward direct comparisons of the predictions of global circulation models to observations via periodogram analysis. [ABSTRACT FROM AUTHOR]

The atmospheric circulation of tidally locked planets is dominated by a superrotating eastward equatorial jet. We develop a predictive theory for the formation of this jet, proposing a mechanism in which the three-dimensional stationary waves induced by the day–night forcing gradient produce an equatorial acceleration. This is balanced in equilibrium by an interaction between the resulting jet and the vertical motion of the atmosphere. The three-dimensional structure of the zonal acceleration is vital to this mechanism. We demonstrate this mechanism in a hierarchy of models. We calculate the three-dimensional stationary waves induced by the forcing on these planets and show the vertical structure of the zonal acceleration produced by these waves, which we use to suggest a mechanism for how the jet forms. General circulation model simulations are used to confirm the equilibrium state predicted by this mechanism, where the acceleration from these waves is balanced by an interaction between the zonal-mean vertical velocity and the jet. We derive a simple model of this using the "Weak Temperature Gradient" approximation, which gives an estimate of the jet speed on a terrestrial tidally locked planet. We conclude that the proposed mechanism is a good description of the formation of an equatorial jet on a terrestrial tidally locked planet and should be useful for interpreting observations and simulations of these planets. The mechanism requires assumptions such as a large equatorial Rossby radius and weak acceleration due to transient waves, and a different mechanism may produce the equatorial jets on gaseous tidally locked planets. [ABSTRACT FROM AUTHOR]

Direct-imaging techniques of exoplanets have made significant progress recently and will eventually enable monitoring of photometric and spectroscopic signals of Earth-like habitable planets. The presence of clouds, however, would remain as one of the most uncertain components in deciphering such direct-imaged signals of planets. We attempt to examine how the planetary obliquity produces different cloud patterns by performing a series of general circulation model simulation runs using a set of parameters relevant for our Earth. Then we use the simulated photometric lightcurves to compute their frequency modulation that is due to the planetary spin–orbit coupling over an entire orbital period, and we attempt to see to what extent one can estimate the obliquity of an Earth twin. We find that it is possible to estimate the obliquity of an Earth twin within the uncertainty of several degrees with a dedicated 4 m space telescope at 10 pc away from the system if the stellar flux is completely blocked. While our conclusion is based on several idealized assumptions, a frequency modulation of a directly imaged Earth-like planet offers a unique methodology to determine its obliquity. [ABSTRACT FROM AUTHOR]

Using a 3D general circulation model (GCM), we investigate the sensitivity of the climate of tidally locked Earthlike exoplanets, Trappist-1e and Proxima Centauri b, to the choice of a convection parameterization. Compared to a mass-flux convection parameterization, a simplified convection adjustment parameterization leads to a >60% decrease of the cloud albedo, increasing the mean dayside temperature by. The representation of convection also affects the atmospheric conditions of the night side, via a change in planetary-scale wave patterns. As a result, using the convection adjustment scheme makes the nightside cold traps warmer by 17–36 K for the planets in our simulations. The day–night thermal contrast is sensitive to the representation of convection in 3D GCM simulations, so caution should be taken when interpreting emission phase curves. The choice of convection treatment, however, does not alter the simulated climate enough to result in a departure from habitable conditions, at least for the atmospheric composition and planetary parameters used in our study. The near-surface conditions both in the Trappist-1e and Proxima b cases remain temperate, allowing for an active water cycle. We further advance our analysis using high-resolution model experiments, in which atmospheric convection is simulated explicitly. Our results suggest that in a hypothetical global convection-permitting simulation, the surface temperature contrast would be higher than in the coarse-resolution simulations with parameterized convection. In other words, models with parameterized convection may overestimate the inter-hemispheric heat redistribution efficiency. [ABSTRACT FROM AUTHOR]

It is well known that planets with radii between that of Earth and Neptune have been the most commonly detected to date. To classify these planets as either terrestrial or gaseous, typically we turn to mass–radius relations and composition curves to determine the likelihood of such a planet being rocky or gaseous. While these methods have set a likely transition radius of approximately 1.5 R⊕, we cannot expect that any change between terrestrial and gaseous compositions will be a sharp cutoff, and composition curve predictions result in ambiguous designations for planets right near this transition radius. In this work, we present 3D general circulation models (GCMs) of transition planets, wherein we study the effects of a surface on observable quantities such as the latitudinal variations and eclipse depths. We present our updated GCM, validated on the circulation of Earth, before discussing our modeling choices for this transition planet. Finally, we discuss the results of this study and explore the prospects of detecting the presence of a surface through observations of secondary eclipses in the future. [ABSTRACT FROM AUTHOR]

Robust atmospheric and radiative transfer modeling will be required to properly interpret reflected-light and thermal emission spectra of terrestrial exoplanets. This will help break observational degeneracies between the numerous atmospheric, planetary, and stellar factors that drive planetary climate. Here, we simulate the climates of earthlike worlds around the Sun with increasingly slow rotation periods, from earthlike to fully Sun-synchronous, using the ROCKE-3D general circulation model. We then provide these results as input to the Spectral Planet Model, which employs the Spectral Mapping Atmospheric Radiative Transfer model to simulate the spectra of a planet as it would be observed from a future space-based telescope. We find that the primary observable effects of slowing planetary rotation rate are the altered cloud distributions, altitudes, and opacities that subsequently drive many changes to the spectra by altering the absorption band depths of biologically relevant gas species (e.g., , , and). We also identify a potentially diagnostic feature of synchronously rotating worlds in mid-infrared absorption/emission lines. [ABSTRACT FROM AUTHOR]

Planets residing in circumstellar habitable zones offer us the best opportunities to test hypotheses of life's potential pervasiveness and complexity. Constraining the precise boundaries of habitability and its observational discriminants is critical to maximizing our chances at remote life detection with future instruments. Conventionally, calculations of the inner edge of the habitable zone (IHZ) have been performed using both 1D radiative-convective and 3D general circulation models. However, these models lack interactive 3D chemistry and do not resolve the mesosphere and lower thermosphere region of the upper atmosphere. Here, we employ a 3D high-top chemistry-climate model (CCM) to simulate the atmospheres of synchronously rotating planets orbiting at the inner edge of habitable zones of K- and M-dwarf stars (between Teff = 2600 and 4000 K). While our IHZ climate predictions are in good agreement with general circulation model studies, we find noteworthy departures in simulated ozone and HOx photochemistry. For instance, climates around inactive stars do not typically enter the classical moist greenhouse regime even with high (≳10−3 mol mol−1) stratospheric water vapor mixing ratios, which suggests that planets around inactive M-stars may only experience minor water-loss over geologically significant timescales. In addition, we find much thinner ozone layers on potentially habitable moist greenhouse atmospheres, as ozone experiences rapid destruction via reaction with hydrogen oxide radicals. Using our CCM results as inputs, our simulated transmission spectra show that both water vapor and ozone features could be detectable by instruments NIRSpec and MIRI LRS on board the James Webb Space Telescope. [ABSTRACT FROM AUTHOR]

The search for water-rich Earth-sized exoplanets around low-mass stars is rapidly gaining attention because they represent the best opportunity to characterize habitable planets in the near future. Understanding the atmospheres of these planets and determining the optimal strategy for characterizing them through transmission spectroscopy with our upcoming instrumentation is essential in order to constrain their environments. For this study, we present simulated transmission spectra of tidally locked Earth-sized ocean-covered planets around late-M to mid-K stellar spectral types, utilizing the results of general circulation models previously published by Kopparapu et al. as inputs for our radiative transfer calculations performed using NASA’s Planetary Spectrum Generator (psg.gsfc.nasa.gov). We identify trends in the depth of H2O spectral features as a function of planet surface temperature and rotation rate. These trends allow us to calculate the exposure times necessary to detect water vapor in the atmospheres of aquaplanets through transmission spectroscopy with the upcoming James Webb Space Telescope as well as several future flagship space telescope concepts under consideration (the Large UV Optical Infrared Surveyor and the Origins Space Telescope) for a target list constructed from the Transiting Exoplanet Survey Satellite (TESS) Input Catalog (TIC). Our calculations reveal that transmission spectra for water-rich Earth-sized planets around low-mass stars will be dominated by clouds, with spectral features <20 ppm, and only a small subset of TIC stars would allow for the characterization of an ocean planet in the habitable zone. We thus present a careful prioritization of targets that are most amenable to follow-up characterizations with next-generation instrumentation, in order to assist the community in efficiently utilizing precious telescope time. [ABSTRACT FROM AUTHOR]

We present gray gas general circulation model (GCM) simulations of the tidally locked mini-Neptune GJ 1214b. On timescales of 1000–10,000 Earth days, our results are comparable to previous studies of the same planet, in the sense that they all exhibit two off-equatorial eastward jets. Over much longer integration times (50,000–250,000 Earth days) we find a significantly different circulation and observational features. The zonal-mean flow transitions from two off-equatorial jets to a single wide equatorial jet that has higher velocity and extends deeper. The hot spot location also shifts eastward over the integration time. Our results imply a convergence time far longer than the typical integration time used in previous studies. We demonstrate that this long convergence time is related to the long radiative timescale of the deep atmosphere and can be understood through a series of simple arguments. Our results indicate that particular attention must be paid to model convergence time in exoplanet GCM simulations, and that other results on the circulation of tidally locked exoplanets with thick atmospheres may need to be revisited. [ABSTRACT FROM AUTHOR]

The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day–night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5–12 μm range with JWST/MIRI. To prepare for these observations, we performed extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 μm spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the nightside emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, eventually providing the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition. [ABSTRACT FROM AUTHOR]

Hot Jupiters receive intense incident stellar light on their daysides, which drives vigorous atmospheric circulation that attempts to erase their large dayside-to-nightside flux contrasts. Propagating waves and instabilities in hot Jupiter atmospheres can cause emergent properties of the atmosphere to be time-variable. In this work, we study such weather in hot Jupiter atmospheres using idealized cloud-free general circulation models with double-gray radiative transfer. We find that hot Jupiter atmospheres can be time-variable at the ∼0.1%–1% level in globally averaged temperature and at the ∼1%–10% level in globally averaged wind speeds. As a result, we find that observable quantities are also time-variable: the secondary eclipse depth can be variable at the ≲2% level, the phase-curve amplitude can change by ≲1%, the phase-curve offset can shift by ≲5°, and terminator-averaged wind speeds can vary by ≲2 km s−1. Additionally, we calculate how the eastern and western limb-averaged wind speeds vary with incident stellar flux and the strength of an imposed drag that parameterizes Lorentz forces in partially ionized atmospheres. We find that the eastern limb is blueshifted in models over a wide range of equilibrium temperature and drag strength, while the western limb is only redshifted if equilibrium temperatures are ≲1500 K and drag is weak. Lastly, we show that temporal variability may be observationally detectable in the infrared through secondary eclipse observations with the James Webb Space Telescope, phase-curve observations with future space telescopes (e.g., ARIEL), and/or Doppler wind speed measurements with high-resolution spectrographs. [ABSTRACT FROM AUTHOR]

The habitable fraction of a planet's surface is important for the detectability of surface biosignatures. The extent and distribution of habitable areas are influenced by external parameters that control the planet's climate, atmospheric circulation, and hydrological cycle. We explore these issues using the ROCKE-3D general circulation model, focusing on terrestrial water fluxes and thus the potential for the existence of complex life on land. Habitability is examined as a function of insolation and planet rotation for an Earth-like world with zero obliquity and eccentricity orbiting the Sun. We assess fractional habitability using an aridity index that measures the net supply of water to the land. Earth-like planets become "superhabitable" (a larger habitable surface area than Earth) as insolation and day-length increase because their climates become more equable, reminiscent of past warm periods on Earth when complex life was abundant and widespread. The most slowly rotating, most highly irradiated planets, though, occupy a hydrological regime unlike any on Earth, with extremely warm, humid conditions at high latitudes but little rain and subsurface water storage. Clouds increasingly obscure the surface as insolation increases, but visibility improves for modest increases in rotation period. Thus, moderately slowly rotating rocky planets with insolation near or somewhat greater than modern Earth's appear to be promising targets for surface characterization by a future direct imaging mission. [ABSTRACT FROM AUTHOR]

The TRAPPIST-1 system, consisting of an ultracool host star having seven known Earth-sized planets, will be a prime target for atmospheric characterization with the James Webb Space Telescope (JWST). However, the detectability of atmospheric molecular species may be severely impacted by the presence of clouds and/or hazes. In this work, we perform 3D general circulation model (GCM) simulations with the LMD-G model supplemented by 1D photochemistry simulations at the terminator with the Atmos model to simulate several possible atmospheres for TRAPPIST-1e, 1f, and 1g: (1) modern Earth, (2) Archean Earth, and (3) CO2-rich atmospheres. The JWST synthetic transit spectra were computed using the GSFC Planetary Spectrum Generator. We find that the TRAPPIST-1e, 1f, and 1g atmospheres, with clouds and/or hazes, could be detected using JWST's NIRSpec Prism from the CO2 absorption line at 4.3 μm in less than 15 transits at 3σ or less than 35 transits at 5σ. However, our analysis suggests that other gases would require hundreds (or thousands) of transits to be detectable. We also find that H2O, mostly confined in the lower atmosphere, is very challenging to detect for these planets or similar systems if the planets' atmospheres are not in a moist greenhouse state. This result demonstrates that the use of GCMs, self-consistently taking into account the effect of clouds and subsaturation, is crucial to evaluate the detectability of atmospheric molecules of interest, as well as for interpreting future detections in a more global (and thus robust and relevant) approach. [ABSTRACT FROM AUTHOR]

We investigate how the meridional circulation and baroclinic eddies change with insolation and rotation rate, under high and zero-obliquity setups, using a general circulation model. The total circulation is considered as superposition of circulations driven by different physics processes, such as diabatic and adiabatic processes. We decompose the meridional circulation into diabatic and adiabatic components, in order to understand their different responses to changes of insolation and rotation rate. As insolation or rotation period increases, the meridional circulation tends to become more diabatically dominant, regardless of the obliquity. The low obliquity circulation is always dominated by diabatic processes, while the high obliquity configuration has two circulation regimes: an adiabatic-dominant regime in the limit of low insolation and fast rotation, and a diabatic-dominant regime in the opposite limit. This regime transition may be observable via its signature on the upper atmospheric zonal wind and the column cloud cover. The momentum-driven circulation, the dominant circulation component in the weak-insolation and fast-rotating regimes is found to resemble that in a dry dynamic model forced by a reversed meridional temperature gradient, indicating the relevance of using a dry dynamic model to understand planetary general circulations under high obliquity. [ABSTRACT FROM AUTHOR]

The potential habitability of known exoplanets is often categorized by a nominal equilibrium temperature assuming a Bond albedo of either ∼0.3, similar to Earth, or 0. As an indicator of habitability, this leaves much to be desired, because albedos of other planets can be very different, and because surface temperature exceeds equilibrium temperature due to the atmospheric greenhouse effect. We use an ensemble of general circulation model simulations to show that for a range of habitable planets, much of the variability of Bond albedo, equilibrium temperature and even surface temperature can be predicted with useful accuracy from incident stellar flux and stellar temperature, two known parameters for every confirmed exoplanet. Earth’s Bond albedo is near the minimum possible for habitable planets orbiting G stars, because of increasing contributions from clouds and sea ice/snow at higher and lower instellations, respectively. For habitable M star planets, Bond albedo is usually lower than Earth’s because of near-IR H2O absorption, except at high instellation where clouds are important. We apply relationships derived from this behavior to several known exoplanets to derive zeroth-order estimates of their potential habitability. More expansive multivariate statistical models that include currently non-observable parameters show that greenhouse gas variations produce significant variance in albedo and surface temperature, while increasing length of day and land fraction decrease surface temperature; insights for other parameters are limited by our sampling. We discuss how emerging information from global climate models might resolve some degeneracies and help focus scarce observing resources on the most promising planets. [ABSTRACT FROM AUTHOR]

The macroturbulent atmospheric circulation of Earth-like planets mediates their equator-to-pole heat transport. For fast-rotating terrestrial planets, baroclinic instabilities in the mid-latitudes lead to turbulent eddies that act to transport heat poleward. In this work, we derive a scaling theory for the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanet atmospheres. This theory is built on the work of Jansen & Ferrari and determines how unstable the atmosphere is to baroclinic instability (the baroclinic “criticality”) through a balance between the baroclinic eddy heat flux and radiative heating/cooling. We compare our scaling theory to General Circulation Model (GCM) simulations and find that the theoretical predictions for equator-to-pole temperature contrast and bulk lapse rate broadly agree with GCM experiments with varying rotation rate and surface pressure throughout the baroclincally unstable regime. Our theoretical results show that baroclinic instabilities are a strong control of heat transport in the atmospheres of Earth-like exoplanets, and our scalings can be used to estimate the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanets. These scalings can be tested by spectroscopic retrievals and full-phase light curves of terrestrial exoplanets with future space telescopes. [ABSTRACT FROM AUTHOR]

On hot Jupiter exoplanets, strong horizontal and vertical winds should homogenize the abundances of the important absorbers CH4 and CO much faster than chemical reactions restore chemical equilibrium. This effect, typically neglected in general circulation models (GCMs), has been suggested to explain discrepancies between observed infrared light curves and those predicted by GCMs. On the nightsides of several hot Jupiters, GCMs predict outgoing fluxes that are too large, especially in the Spitzer 4.5 μm band. We modified the SPARC/MITgcm to include disequilibrium abundances of CH4, CO, and H2O by assuming that the CH4/CO ratio is constant throughout the simulation domain. We ran simulations of hot Jupiter HD 189733b with eight CH4/CO ratios. In the more likely CO-dominated regime, we find temperature changes ≥50–100 K compared to the simulation for equilibrium chemistry across large regions. This effect is large enough to affect predicted emission spectra and should thus be included in GCMs of hot Jupiters with equilibrium temperatures between 600 and 1300 K. We find that spectra in regions with strong methane absorption, including the Spitzer 3.6 and 8 μm bands, are strongly impacted by disequilibrium abundances. We expect chemical quenching to result in much larger nightside fluxes in the 3.6 μm band, in stark contrast to observations. Meanwhile, we find almost no effect on predicted observations in the 4.5 μm band, because the changes in opacity due to CO and H2O offset each other. We thus conclude that disequilibrium carbon chemistry cannot explain the observed low nightside fluxes in the 4.5 μm band. [ABSTRACT FROM AUTHOR]

Exploring diverse planetary atmospheres requires modeling tools that are both accurate and flexible. Here, we develop a three-dimensional general circulation model (3D GCM) that, for the first time, uses a line-by-line approach to describe the radiative transfer. We validate our GCM by comparing with published results done by different 1D and 3D models. To demonstrate the versatility of the model, we apply the GCM to the hot Earth-sized exoplanet GJ 1132b and study its climate and circulation assuming an atmosphere dominated by abiotic oxygen (O2). Our simulations show that a minor CO2 composition can change the circulation pattern substantially, intensifying the equatorial superrotation in particular. Computation of the phase-resolved spectroscopy indicates that the vertical profile of the superrotating jet could be inferred in future spectrophotometric observations by the phase shift of the hotspot in the CO2 principle absorption band centered at 667 cm−1. We also show that atmospheric mass could potentially be constrained by the phase amplitude in the O2 vibrational fundamental band for planets with O2-rich atmospheres, although further experimental and/or theoretical O2–O2 collision-induced absorption data at high temperatures is needed to confirm this. More physical schemes, such as moist dynamics, will be implemented in the GCM in the future so that it can be used to tackle a wide variety of planetary climate problems. [ABSTRACT FROM AUTHOR]

Spectral features in the observed spectra of exoplanets depend on the composition of their atmospheres. A good knowledge of the main atmospheric processes that drive the chemical distribution is therefore essential to interpret exoplanetary spectra. An atmosphere reaches chemical equilibrium if the rates of the forward and backward chemical reactions converge to the same value. However, there are atmospheric processes, such as atmospheric transport, that destabilize this equilibrium. In this work we study the changes in composition driven by a 3D wind field in WASP-43b using our Global Circulation Model, THOR. Our model uses validated temperature- and pressure-dependent chemical timescales that allow us to explore the disequilibrium chemistry of CO, CO2, H2O, and CH4. In WASP-43b the formation of the equatorial jet has an important impact on the chemical distribution of the different species across the atmosphere. At low latitudes the chemistry is longitudinally quenched, except for CO2 at solar abundances. The polar vortexes have a distinct chemical distribution since these are regions with lower temperature and atmospheric mixing. Vertical and latitudinal mixing have a secondary impact on the chemical transport. We determine graphically the effect of disequilibrium on the observed emission spectra. Our results do not show any significant differences in the emission spectra between the equilibrium and disequilibrium solutions for C/O = 0.5. However, if C/O is increased to 2.0, differences in the spectra due to the disequilibrium chemistry of CH4 become non-negligible. In some spectral ranges the emission spectra can have more than 15% departure from the equilibrium solution. [ABSTRACT FROM AUTHOR]

It has recently been proposed that Earth-like planets in the outer regions of the habitable zone experience unstable climates, repeatedly cycling between glaciated and deglaciated climatic states. While this result has been confirmed and also extended to explain early Mars climate records, all existing work relies on highly idealized low-dimensional climate models. Here, we confirm that the phenomenology of climate cycles remains in 3D Earth climate models with considerably more degrees of freedom. To circumvent the computational barrier of integrating climate on Gyr timescales, we implement a hybrid 0D-3D integrator that uses a general circulation model (GCM) as a short relaxation step along a long evolutionary climate sequence. We find that GCM climate cycles are qualitatively consistent with reported low-dimensional results. This establishes on a firmer ground the notion that outer habitable zone planets may be preferentially found in transiently glaciated states. [ABSTRACT FROM AUTHOR]

Conventional definitions of habitability require abundant liquid surface water to exist continuously over geologic timescales. Water in each of its thermodynamic phases interacts with solar and thermal radiation and is the cause for strong climatic feedbacks. Thus, assessments of the habitable zone require models to include a complete treatment of the hydrological cycle over geologic time. Here, we use the Community Atmosphere Model from the National Center for Atmospheric Research to study the evolution of climate for an Earth-like planet at constant CO2, under a wide range of stellar fluxes from F-, G-, and K-dwarf main sequence stars. Around each star we find four stable climate states defined by mutually exclusive global mean surface temperatures (Ts); snowball (Ts ≤ 235 K), waterbelt (235 K ≤ Ts ≤ 250 K), temperate (275 K ≤ Ts ≤ 315 K), and moist greenhouse (Ts ≥ 330 K). Each is separated by abrupt climatic transitions. Waterbelt, temperate, and cooler moist greenhouse climates can maintain open-ocean against both sea ice albedo and hydrogen escape processes respectively, and thus constitute habitable worlds. We consider the warmest possible habitable planet as having Ts ∼ 355 K, at which point diffusion limited water-loss could remove an Earth ocean in ∼1 Gyr. Without long timescale regulation of non-condensable greenhouse species at Earth-like temperatures and pressures, such as CO2, habitability can be maintained for an upper limit of ∼2.2, ∼2.4, and ∼4.7 Gyr around F-, G-, and K-dwarf stars respectively, due to main sequence brightening. [ABSTRACT FROM AUTHOR]

Super Earths and mini Neptunes likely have a wide range of atmospheric compositions, ranging from low molecular mass atmospheres of H2 to higher molecular atmospheres of water, CO2, N2, or other species. Here we systematically investigate the effects of atmospheric bulk compositions on temperature and wind distributions for tidally locked sub-Jupiter-sized planets, using an idealized 3D general circulation model (GCM). The bulk composition effects are characterized in the framework of two independent variables: molecular weight and molar heat capacity. The effect of molecular weight dominates. As the molecular weight increases, the atmosphere tends to have a larger day–night temperature contrast, a smaller eastward phase shift in the thermal phase curve, and a smaller zonal wind speed. The width of the equatorial super-rotating jet also becomes narrower, and the “jet core” region, where the zonal-mean jet speed maximizes, moves to a greater pressure level. The zonal-mean zonal wind is more prone to exhibit a latitudinally alternating pattern in a higher molecular weight atmosphere. We also present analytical theories that quantitatively explain the above trends and shed light on the underlying dynamical mechanisms. Those trends might be used to indirectly determine the atmospheric compositions on tidally locked sub-Jupiter-sized planets. The effects of the molar heat capacity are generally small. But if the vertical temperature profile is close to adiabatic, molar heat capacity will play a significant role in controlling the transition from a divergent flow in the upper atmosphere to a jet-dominated flow in the lower atmosphere. [ABSTRACT FROM AUTHOR]

We have designed and developed, from scratch, a global circulation model (GCM) named THOR that solves the three-dimensional nonhydrostatic Euler equations. Our general approach lifts the commonly used assumptions of a shallow atmosphere and hydrostatic equilibrium. We solve the “pole problem” (where converging meridians on a sphere lead to increasingly smaller time steps near the poles) by implementing an icosahedral grid. Irregularities in the grid, which lead to grid imprinting, are smoothed using the “spring dynamics” technique. We validate our implementation of spring dynamics by examining calculations of the divergence and gradient of test functions. To prevent the computational time step from being bottlenecked by having to resolve sound waves, we implement a split-explicit method together with a horizontally explicit and vertically implicit integration. We validate our GCM by reproducing the Earth and hot-Jupiter-like benchmark tests. THOR was designed to run on graphics processing units (GPUs), which allows for physics modules (radiative transfer, clouds, chemistry) to be added in the future, and is part of the open-source Exoclimes Simulation Platform (www.exoclime.org). [ABSTRACT FROM AUTHOR]

Circumbinary planets experience a time-varying irradiation pattern as they orbit their two host stars. In this work, we present the first detailed study of the atmospheric effects of this irradiation pattern on known and hypothetical gaseous circumbinary planets. Using both a one-dimensional energy balance model (EBM) and a three-dimensional general circulation model (GCM), we look at the temperature differences between circumbinary planets and their equivalent single-star cases in order to determine the nature of the atmospheres of these planets. We find that for circumbinary planets on stable orbits around their host stars, temperature differences are on average no more than 1.0% in the most extreme cases. Based on detailed modeling with the GCM, we find that these temperature differences are not large enough to excite circulation differences between the two cases. We conclude that gaseous circumbinary planets can be treated as their equivalent single-star case in future atmospheric modeling efforts. [ABSTRACT FROM AUTHOR]