Search

Although exoplanetary science was not initially projected to be a substantial part of the Spitzer mission, its exoplanet observations set the stage for current and future surveys with JWST and Ariel. We present a comprehensive reduction and analysis of Spitzer's 4.5 micron phase curves of 29 hot Jupiters on low-eccentricity orbits. The analysis, performed with the Spitzer Phase Curve Analysis (SPCA) pipeline, confirms that BLISS mapping is the best detrending scheme for most, but not all, observations. Visual inspection remains necessary to ensure consistency across detrending methods due to the diversity of phase curve data and systematics. Regardless of the model selection scheme - whether using the lowest-BIC or a uniform detrending approach - we observe the same trends, or lack thereof. We explore phase curve trends as a function of irradiation temperature, orbital period, planetary radius, mass, and stellar effective temperature. We discuss the trends that are robustly detected and provide potential explanations for those that are not observed. While it is almost tautological that planets receiving greater instellation are hotter, we are still far from confirming dynamical theories of heat transport in hot Jupiter atmospheres due to the sample's diversity. Even among planets with similar temperatures, other factors like rotation and metallicity vary significantly. Larger, curated sample sizes and higher-fidelity phase curve measurements from JWST and Ariel are needed to firmly establish the parameters governing day-night heat transport on synchronously rotating planets., Comment: 16 pages, 6 figures, submitted to AAS journal

The geologically active south pole of Enceladus generates a plume of micron-sized particles, which likely form Saturn's tenuous E-ring extending from the orbit of Mimas to Titan. Interactions between these particles and satellites have been suggested, though only as very thin surficial phenomena. We scrutinize high-resolution images with a newly developed numerical shape model of Helene and find that the leading hemisphere of Helene is covered by thick deposits of E-ring particles, which occasionally collapse to form gully-like depressions. The depths of the resultant gullies and near-absence of small craters on the leading hemisphere indicate that the deposit is tens to hundreds of meters thick. The ages of the deposits are less than several tens of My, which coincides well with similar deposits found on Telesto and Calypso. Our findings as well as previous theoretical work collectively indicate that the cryovolcanic activity currently occurring on Enceladus is ephemeral., Comment: 23 pages 6 figures

With over 30 phase curves observed during the warm Spitzer mission, the complete data set provides a wealth of information relating to trends and three-dimensional properties of hot Jupiter atmospheres. In this work we present a comparative study of seven new Spitzer phase curves for four planets with equilibrium temperatures of T$_{eq}\sim$ 1300K: Qatar-2b, WASP-52b, WASP-34b, and WASP-140b, as well as the reanalysis of the 4.5 $\micron$ Qatar-1b phase curve due to the similar equilibrium temperature. In total, five 4.5 $\micron$ phase curves and three 3.6 $\micron$ phase curves are analyzed here with a uniform approach. Using these new results, in combination with literature values for the entire population of published Spitzer phase curves of hot Jupiters, we present evidence for a linear trend of increasing hot spot offset with increasing orbital period, as well as observational evidence for two classes of planets in apparent redistribution vs. equilibrium temperature parameter space, and tentative evidence for a dependence of hot spot offset on planetary surface gravity in our $\sim$ 1300 K sample. We do not find trends in apparent heat redistribution with orbital period or gravity. Non-uniformity in literature Spitzer data analysis techniques precludes a definitive determination of the sources or lack of trends., Comment: 22 pages, 13 figures, 6 tables. Accepted for publication in AAS journals

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: superrotation regime, midlatitude-jet regime, and 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., Comment: accepted by The Astrophysical Journal

We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $\mu$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $\mu$m and $1610 \pm 70$ ppm at 4.5 $\mu$m show no evidence of the previously reported dayside temperature inversion. We also empirically derive the mass and radius of XO-3b and its host star using Gaia DR3's parallax measurement and find a planetary mass $M_p=11.79 \pm 0.98 ~M_{\rm{Jup}}$ and radius $R_p=1.295 \pm 0.066 ~R_{\rm{Jup}}$. We compare our \textit{Spitzer} observations with multiple atmospheric models to constrain the radiative and advective properties of XO-3b. While the decorrelated 4.5 $\mu$m observations are pristine, the 3.6 $\mu$m phase curve remains polluted with detector systematics due to larger amplitude intrapixel sensitivity variations in this channel. We focus our analysis on the more reliable 4.5 $\mu$m phase curve and fit an energy balance model with solid body rotation to estimate the zonal wind speed and the pressure of the bottom of the mixed layer. Our energy balance model fit suggests an eastward equatorial wind speed of $3.13 ^{+0.26} _{-0.83}$ km/s, an atmospheric mixed layer down to $2.40 ^{+0.92} _{-0.16}$ bar, and Bond albedo of $0.106 ^{+0.008} _{-0.106}$. We assume that the wind speed and mixed layer depth are constant throughout the orbit. We compare our observations with a 1D planet-averaged model predictions at apoapse and periapse and 3D general circulation model (GCM) predictions for XO-3b. We also investigate the inflated radius of XO-3b and find that it would require an unusually large amount of internal heating to explain the observed planetary radius., Comment: 23 pages, 14 figures, accepted for publication in AJ

The large radii of many hot Jupiters can only be matched by models that have hot interior adiabats, and recent theoretical work has shown that the interior evolution of hot Jupiters has a significant impact on their atmospheric structure. Due to its inflated radius, low gravity, and ultra-hot equilibrium temperature, WASP-76b is an ideal case study for the impact of internal evolution on observable properties. Hot interiors should most strongly affect the non-irradiated side of the planet, and thus full phase curve observations are critical to ascertain the effect of the interior on the atmospheres of hot Jupiters. In this work, we present the first Spitzer phase curve observations of WASP-76b. We find that WASP-76b has an ultra-hot day side and relatively cold nightside with brightness temperatures of $2471 \pm 27~\mathrm{K}$/$1518 \pm 61~\mathrm{K}$ at $3.6~\micron$ and $2699 \pm 32~\mathrm{K}$/$1259 \pm 44~\mathrm{K}$ at $4.5~\micron$, respectively. These results provide evidence for a dayside thermal inversion. Both channels exhibit small phase offsets of $0.68 \pm 0.48^{\circ}$ at $3.6~\micron$ and $0.67 \pm 0.2^{\circ}$ at $4.5~\mu\mathrm{m}$. We compare our observations to a suite of general circulation models that consider two end-members of interior temperature along with a broad range of frictional drag strengths. Strong frictional drag is necessary to match the small phase offsets and cold nightside temperatures observed. From our suite of cloud-free GCMs, we find that only cases with a cold interior can reproduce the cold nightsides and large phase curve amplitude at $4.5~\micron$, hinting that the hot interior adiabat of WASP-76b does not significantly impact its atmospheric dynamics or that clouds blanket its nightside., Comment: 24 pages, 10 Figures, 5 Tables. Accepted to AJ. Co-First Authors

Brown dwarfs and directly imaged exoplanets exhibit observational evidence for active atmospheric circulation, raising critical questions about mechanisms driving the circulation, its fundamental nature, and time variability. Our previous work demonstrated the crucial role of cloud radiative feedback on driving a vigorous atmospheric circulation using local models that assume a Cartesian geometry and constant Coriolis parameters. In this study, we explore the properties of the global dynamics. We show that, under relatively strong dissipation in the bottom layers of the model, horizontally isotropic vortices are prevalent at mid-to-high latitudes while large-scale zonally propagating waves are dominant at low latitudes near the observable layers. The equatorial waves have both eastward and westward phase speeds, and the eastward components with typical speeds of a few hundred m/s usually dominate the equatorial time variability. Lightcurves of the global simulations show variability with amplitudes from 0.5 percent to a few percent depending on the rotation period and viewing angle. The time evolution of simulated lightcurves is critically affected by the equatorial waves, showing wave beating effects and differences in the lightcurve periodicity to the intrinsic rotation period. The vertical extent of clouds is the largest at the equator and decreases poleward due to the increasing influence of rotation with increasing latitude. Under weaker bottom dissipation, strong and broad zonal jets develop and modify wave propagation and lightcurve variability. Our modeling results help to explain the puzzling time evolution of observed lightcurves, a slightly shorter period of variability in IR than in radio wavelengths, and the viewing angle dependence of variability amplitude and IR colors., Comment: Accepted to MNRAS, 23 pages, 18 figures

Photochemical hazes have been suggested as candidate for the high-altitude aerosols observed in the transmission spectra of many hot Jupiters. We present 3D simulations of the hot Jupiter HD 189733b to study how photochemical hazes are transported by atmospheric circulation. The model includes spherical, constant-size hazes particles that gravitationally settle and are transported by the winds as passive tracers, with particle radii ranging from 1 nm to 300 $\mu$m. We identify two general types of haze distribution based on particle size: In the small-particle regime (<30 nm), gravitational settling is unimportant, and hazes accumulate in two large mid-latitude vortices centered on the night side that extend across the morning terminator. Therefore, small hazes are more concentrated at the morning terminator than at the evening terminator. In the large-particle regime (>30 nm), hazes settle out quickly on the nightside, resulting in more hazes at the evening terminator. For small particles, terminator differences in haze mass mixing ratio and temperature considered individually can result in significant differences in the transit spectra of the terminators. When combining both effects for HD189733b, however, they largely cancel out each other, resulting in very small terminator differences in the spectra. Transit spectra based on the GCM-derived haze distribution fail to reproduce the steep spectral slope at short wavelengths in the current transit observations of HD 189733b. Differing optical properties of hazes, hotter temperatures at low pressures because of heating by hazes, enhanced sub-grid-scale mixing, or star spots might explain the mismatch between the model and observations., Comment: This version incorporates the corrected versions of equations (7) and (8), as described in the correction published in July 2024, available at \url{https://academic.oup.com/mnras/article/531/3/3444/7691150}

Hot Jupiters have been predicted to have a strong day/night temperature contrast and a hot spot shifted eastward of the substellar point. This was confirmed by numerous phase curve observations probing the longitudinal brightness variation of the atmosphere. Global circulation models, however, systematically underestimate the phase curve amplitude and overestimate the shift of its maximum. We use a global circulation model including non-grey radiative transfer and realistic gas and cloud opacities to systematically investigate how the atmospheric circulation of hot Jupiters varies with equilibrium temperature from 1000 to 2200K. We show that the heat transport is very efficient for cloudless planets cooler than 1600K and becomes less efficient at higher temperatures. When nightside clouds are present, the day-to-night heat transport becomes extremely inefficient, leading to a good match to the observed low nightside temperatures. The constancy of this low temperature is, however, due to the strong dependence of the radiative timescale with temperature. We further show that nightside clouds increase the phase curve amplitude and decreases the phase curve offset at the same time. This change is very sensitive to the cloud chemical composition and particle size, meaning that the diversity in observed phase curves can be explained by a diversity of nightside cloud properties. Finally, we show that phase curve parameters do not necessarily track the day/night contrast nor the shift of the hot spot on isobars, and propose solutions to to recover the true hot-spot shift and day/night contrast., Comment: 32 pages, 21 figures, 4 years in the making. Accepted in MNRAS. See also the complementary study by Roman et al. today on ArXiv

Ground-based and spacecraft telescopic observations, combined with an intensive modeling effort, have greatly enhanced our understanding of hot giant planets and brown dwarfs over the past ten years. Although these objects are all fluid, hydrogen worlds with stratified atmospheres overlying convective interiors, they exhibit an impressive diversity of atmospheric behavior. Hot Jupiters are strongly irradiated, and a wealth of observations constrain the day-night temperature differences, circulation, and cloudiness. The intense stellar irradiation, presumed tidal locking and modest rotation leads to a novel regime of strong day-night radiative forcing. Circulation models predict large day-night temperature differences, global-scale eddies, patchy clouds, and, in most cases, a fast eastward jet at the equator-equatorial superrotation. The warm Jupiters may exhibit a wide range of rotation rates, obliquities, and orbital eccentricities, which, along with the weaker irradiation, leads to circulation patterns and observable signatures predicted to differ substantially from hot Jupiters. Brown dwarfs are typically isolated, rapidly rotating worlds; they radiate enormous energy fluxes into space and convect vigorously in their interiors. Their atmospheres exhibit patchiness in clouds and temperature on regional to global scales-the result of modulation by large-scale atmospheric circulation. Despite the lack of irradiation, such circulations can be driven by interaction of the interior convection with the overlying atmosphere, as well as self-organization of patchiness due to cloud-dynamical-radiative feedbacks. Finally, irradiated brown dwarfs help to bridge the gap between these classes of objects, experiencing intense external irradiation as well as vigorous interior convection. A hierarchy of modeling approaches have yielded major new insights into the dynamics governing these phenomena., Comment: Invited review for Space Science Reviews special issue on "Understanding the Diversity of Planetary Atmospheres". Accepted. 52 pages, 33 figures. This may be the last published work led by Adam Showman, whose sudden death during the revision process of this manuscript deprived the field of one of his giant. He will be missed by all

Observations of brown dwarfs (BDs), free-floating planetary-mass objects, and directly imaged extrasolar giant planets (EGPs) exhibit rich evidence of large-scale weather. Cloud radiative feedback has been proposed as a potential mechanism driving the vigorous atmospheric circulation on BDs and directly imaged EGPs, and yet it has not been demonstrated in three-dimensional dynamical models at relevant conditions. Here we present a series of atmospheric circulation models that self-consistently coupled dynamics with idealized cloud formation and its radiative effects. We demonstrate that vigorous atmospheric circulation can be triggered and self-maintained by cloud radiative feedback. Typical isobaric temperature variation could reach over 100 K and horizontally averaged wind speed could be several hundred m/s. The circulation is dominated by cloud-forming and clear-sky vortices that evolve over timescales from several to tens of hours. The typical horizontal lengthscale of dominant vortices is closed to the Rossby deformation radius, showing a linear dependence on the inverse of rotation rate. Stronger rotation tends to weaken the vertical transport of vapor and clouds, leading to overall thinner clouds. Domain-mean outgoing radiative flux exhibits variability over timescales of tens of hours due to the statistical evolution of storms. Different bottom boundary conditions in the models could lead to qualitatively different circulation near the observable layer. The circulation driven by cloud radiative feedback represents a robust mechanism generating significant surface inhomogeneity as well as irregular flux time variability. Our results have important implications for near-IR colors of dusty BDs and EGPs, including the scatter in the near-IR color-magnitude diagram and the viewing-geometry dependent near-IR colors., Comment: Accepted to MNRAS, 23 pages, 16 figures

We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~\mu$m and $0.30 \pm 0.02 \%$ at $4.5~\mu$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~\mu$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~\mu$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters., Comment: 14 pages, 8 figures. Accepted for publication in AJ

Tidally locked gas giants are typically in several-day orbits, implying a modest role for rotation in the atmospheric circulation. Nevertheless, there exist a class of gas-giant, highly irradiated objects---brown dwarfs orbiting white dwarfs in extremely tight orbits---whose orbital and hence rotation periods are as short as 1-2 hours. Phase curves and other observations have already been obtained for this class of objects, raising fundamental questions about the role of increasing planetary rotation rate in controlling the circulation. So far, most modeling studies have investigated rotation periods exceeding a day, as appropriate for typical hot Jupiters. Here we investigate atmospheric circulation of tidally locked atmospheres with decreasing rotation periods down to 2.5 hours. With decreasing rotation period, the width of the equatorial eastward jet decreases, consistent with the narrowing of the equatorial waveguide due to a decrease of the equatorial deformation radius. The eastward-shifted equatorial hot spot offset decreases accordingly, and the off-equatorial westward-shifted hot areas become increasingly distinctive. At high latitudes, winds become weaker and more rotationally dominated. The day-night temperature contrast becomes larger due to the stronger influence of rotation. Our simulated atmospheres exhibit variability, presumably caused by instabilities and wave interactions. Unlike typical hot Jupiter models, thermal phase curves of rapidly rotating models show a near alignment of peak flux to secondary eclipse. This result helps to explain why, unlike hot Jupiters, many brown dwarfs orbiting white dwarfs exhibit IR flux peaks aligned with secondary eclipse. Our results have important implications for understanding fast-rotating, tidally locked atmospheres., Comment: 22 pages, 13 figures, accepted to ApJ

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-grey radiative transfer. We find that hot Jupiter atmospheres can be time-variable at the $\sim 0.1-1\%$ level in globally averaged temperature and at the $\sim 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 $\lesssim 2\%$ level, the phase curve amplitude can change by $\lesssim 1\%$, the phase curve offset can shift by $\lesssim 5^{\circ}$, and terminator-averaged wind speeds can vary by $\lesssim 2~ \mathrm{km}~\mathrm{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 $\lesssim1500~\mathrm{K}$ and drag is weak. Lastly, we show that temporal variability may be observationally detectable in the infrared through secondary eclipse observations with JWST, phase curve observations with future space telescopes (e.g., ARIEL), and/or Doppler wind speed measurements with high-resolution spectrographs., Comment: 18 pages, 12 figures, updated to reflect published version

The nature and structure of the observed east-west flows on Jupiter and Saturn has been one of the longest-lasting mysteries in planetary science. This mystery has been recently unraveled due to the accurate gravity measurements provided by the Juno mission to Jupiter and the Grand Finale of the Cassini mission to Saturn. These two experiments, which coincidentally happened around the same time, allowed determination of the vertical and meridional profiles of the zonal flows on both planets. This paper reviews the topic of zonal jets on the gas giants in light of the new data from these two experiments. The gravity measurements not only allow the depth of the jets to be constrained, yielding the inference that the jets extend roughly 3000 and 9000 km below the observed clouds on Jupiter and Saturn, respectively, but also provide insights into the mechanisms controlling these zonal flows. Specifically, for both planets this depth corresponds to the depth where electrical conductivity is within an order of magnitude of 1 S/m, implying that the magnetic field likely plays a key role in damping the zonal flows., Comment: Submitted to Space Science Reviews. Part of ISSI special collection on Diversity of Atmospheres

The goal of Project GAUSS is to return samples from the dwarf planet Ceres. Ceres is the most accessible ocean world candidate and the largest reservoir of water in the inner solar system. It shows active cryovolcanism and hydrothermal activities in recent history that resulted in minerals not found in any other planets to date except for Earth's upper crust. The possible occurrence of recent subsurface ocean on Ceres and the complex geochemistry suggest possible past habitability and even the potential for ongoing habitability. Aiming to answer a broad spectrum of questions about the origin and evolution of Ceres and its potential habitability, GAUSS will return samples from this possible ocean world for the first time. The project will address the following top-level scientific questions: 1) What is the origin of Ceres and the origin and transfer of water and other volatiles in the inner solar system? 2) What are the physical properties and internal structure of Ceres? What do they tell us about the evolutionary and aqueous alteration history of icy dwarf planets? 3) What are the astrobiological implications of Ceres? Was it habitable in the past and is it still today? 4) What are the mineralogical connections between Ceres and our current collections of primitive meteorites? GAUSS will first perform a high-resolution global remote sensing investigation, characterizing the geophysical and geochemical properties of Ceres. Candidate sampling sites will then be identified, and observation campaigns will be run for an in-depth assessment of the candidate sites. Once the sampling site is selected, a lander will be deployed on the surface to collect samples and return them to Earth in cryogenic conditions that preserves the volatile and organic composition as well as the original physical status as much as possible., Comment: Section 3.4 revised; List of team members updated; Typos corrected

The atmospheres of the four giant planets of our Solar System share a common and well-observed characteristic: they each display patterns of planetary banding, with regions of different temperatures, composition, aerosol properties and dynamics separated by strong meridional and vertical gradients in the zonal (i.e., east-west) winds. On Jupiter, the reflective white bands of low temperatures, elevated aerosol opacities, and enhancements of quasi-conserved chemical tracers are referred to as 'zones.' Conversely, the darker bands of warmer temperatures, depleted aerosols, and reductions of chemical tracers are known as `belts.' On Saturn, we define cyclonic belts and anticyclonic zones via their temperature and wind characteristics, although their relation to Saturn's albedo is not as clear as on Jupiter. On distant Uranus and Neptune, the exact relationships between the banded albedo contrasts and the environmental properties is a topic of active study. This review is an attempt to reconcile the observed properties of belts and zones with (i) the meridional overturning inferred from the convergence of eddy angular momentum into the eastward zonal jets at the cloud level on Jupiter and Saturn and the prevalence of moist convective activity in belts; and (ii) the opposing meridional motions inferred from the upper tropospheric temperature structure, which implies decay and dissipation of the zonal jets with altitude above the clouds. These two scenarios suggest meridional circulations in opposing directions, the former suggesting upwelling in belts, the latter suggesting upwelling in zones. This presents an unresolved paradox for our current understanding of the banded structure of giant planet atmospheres, that could be addressed via a multi-tiered vertical structure of 'stacked circulation cells.' [Abridged], Comment: 25 pages, 6 figures, Space Science Reviews, accepted

Launching in 2028, ESA's Atmospheric Remote-sensing Exoplanet Large-survey (ARIEL) survey of $\sim$1000 transiting exoplanets will build on the legacies of Kepler and TESS and complement JWST by placing its high precision exoplanet observations into a large, statistically-significant planetary population context. With continuous 0.5--7.8~$\mu$m coverage from both FGS (0.50--0.55, 0.8--1.0, and 1.0--1.2~$\mu$m photometry; 1.25--1.95~$\mu$m spectroscopy) and AIRS (1.95--7.80~$\mu$m spectroscopy), ARIEL will determine atmospheric compositions and probe planetary formation histories during its 3.5-year mission. NASA's proposed Contribution to ARIEL Spectroscopy of Exoplanets (CASE) would be a subsystem of ARIEL's FGS instrument consisting of two visible-to-infrared detectors, associated readout electronics, and thermal control hardware. FGS, to be built by the Polish Academy of Sciences' Space Research Centre, will provide both fine guiding and visible to near-infrared photometry and spectroscopy, providing powerful diagnostics of atmospheric aerosol contribution and planetary albedo, which play a crucial role in establishing planetary energy balance. The CASE team presents here an independent study of the capabilities of ARIEL to measure exoplanetary metallicities, which probe the conditions of planet formation, and FGS to measure scattering spectral slopes, which indicate if an exoplanet has atmospheric aerosols (clouds and hazes), and geometric albedos, which help establish planetary climate. Our design reference mission simulations show that ARIEL could measure the mass-metallicity relationship of its 1000-planet single-visit sample to $>7.5\sigma$ and that FGS could distinguish between clear, cloudy, and hazy skies and constrain an exoplanet's atmospheric aerosol composition to $>5\sigma$ for hundreds of targets, providing statistically-transformative science for exoplanet atmospheres., Comment: accepted to PASP; 23 pages, 6 figures

We study the properties of photochemical hazes in super-Earths/mini-Neptunes atmospheres with particular focus on GJ1214b. We evaluate photochemical haze properties at different metallicities between solar and 10000$\times$solar. Within the four orders of magnitude change in metallicity, we find that the haze precursor mass fluxes change only by a factor of $\sim$3. This small diversity occurs with a non-monotonic manner among the different metallicity cases, reflecting the interaction of the main atmospheric gases with the radiation field. Comparison with relative haze yields at different metallicities from laboratory experiments reveals a qualitative similarity with our theoretical calculations and highlights the contributions of different gas precursors. Our haze simulations demonstrate that higher metallicity results into smaller average particle sizes. Metallicities at and above 100$\times$solar with haze formation yields of $\sim$10$\%$ provide enough haze opacity to satisfy transit observation at visible wavelengths and obscure sufficiently the H$_2$O molecular absorption features between 1.1 $\mu$m and 1.7 $\mu$m. However, only the highest metallicity case considered (10000$\times$solar) brings the simulated spectra into closer agreement with transit depths at 3.6 $\mu$m and 4.5 $\mu$m indicating a high contribution of CO/CO$_2$ in GJ1214b's atmosphere. We also evaluate the impact of aggregate growth in our simulations, in contrast to spherical growth, and find that the two growth modes provide similar transit signatures (for D$_f$=2), but with different particle size distributions. Finally, we conclude that the simulated haze particles should have major implications for the atmospheric thermal structure and for the properties of condensation clouds.

Aerosols appear to be ubiquitous in close-in gas giant atmospheres, and disequilibrium chemistry likely impacts the emergent spectra of these planets. Lofted aerosols and disequilibrium chemistry are caused by vigorous vertical transport in these heavily irradiated atmospheres. Here we numerically and analytically investigate how vertical transport should change over the parameter space of spin-synchronized gas giants. In order to understand how tracer transport depends on planetary parameters, we develop an analytic theory to predict vertical velocities and mixing rates ($K_\mathrm{zz}$) and compare the results to our numerical experiments. We find that both our theory and numerical simulations predict that, if the vertical mixing rate is described by an eddy diffusivity, then this eddy diffusivity $K_\mathrm{zz}$ should increase with increasing equilibrium temperature, decreasing frictional drag strength, and increasing chemical loss timescales. We find that the transition in our numerical simulations between circulation dominated by a superrotating jet and that with solely day-to-night flow causes a marked change in the vertical velocity structure and tracer distribution. The mixing ratio of passive tracers is greatest for intermediate drag strengths that corresponds to this transition between a superrotating jet with columnar vertical velocity structure and day-to-night flow with upwelling on the dayside and downwelling on the nightside. Lastly, we present analytic solutions for $K_\mathrm{zz}$ as a function of planetary effective temperature, chemical loss timescales, and other parameters, for use as input to one-dimensional chemistry models of spin-synchronized gas giant atmospheres., Comment: 25 pages, 12 figures, Accepted at ApJ

Here we present the analysis of multi-epoch secondary eclipse observations of HD 189733b and HD 209458b as a probe of temporal variability in the planetary climate using both Spitzer channels 1 and 2 (3.6 and 4.5 um). Constraining temporal variability will inform models and identify physical processes occurring at either length scales too small to directly observe or at pressure levels that are inaccessible to transit observations. We do not detect statistically significant variability and are able to place useful upper limits on the IR variability amplitudes in these atmospheres. There are very few planets with multi-epoch observations at the required precision to probe variability in dayside emission. The observations considered in this study span several years, providing insight into temporal variability at multiple timescales. In the case of HD 189733b, the best fit eclipse depths for the channel 2 observations exhibit a scatter of 102 ppm about a median depth of 1827 ppm and in channel 1 exhibit a scatter of 88 ppm about a median depth of 1481 ppm. For HD 209458b, the best fit eclipse depths for the channel 2 observations exhibit a scatter of 22 ppm about a median depth of 1406 ppm and in channel 1 exhibit a scatter of 131 ppm about a median depth of 1092 ppm. The precision and scatter in these observations allow us to constrain variability to less than (5.6% and 6.0%) and (12% and 1.6%) for channels (1,2) of HD 189733b and HD 209458b respectively. There is a difference in the best fit eclipse timing compared to the predicted time consistent with an offset hotspot as predicted by GCMs and confirmed in previous phase curve observations., Comment: Accepted for Publication in AJ 12/17/2019, 9 pages, 8 figures

We present the analysis of a full-orbit, spectroscopic phase curve of the ultra hot Jupiter WASP-18b, obtained with the Wide Field Camera 3 aboard the Hubble Space Telescope. We measure the planet's normalized day-night contrast as >0.96 in luminosity: the disk-integrated dayside emission from the planet is at 964+-25 ppm, corresponding to 2894+-30 K, and we place an upper limit on the nightside emission of <32ppm or 1430K at the 3-sigma level. We also find that the peak of the phase curve exhibits a small, but significant offset in brightness of 4.5+-0.5 degrees eastward. We compare the extracted phase curve and phase resolved spectra to 3D Global Circulation Models and find that broadly the data can be well reproduced by some of these models. We find from this comparison several constraints on the atmospheric properties of the planet. Firstly we find that we need efficient drag to explain the very inefficient day-night re-circulation observed. We demonstrate that this drag could be due to Lorentz-force drag by a magnetic field as weak as 10 Gauss. Secondly, we show that a high metallicity is not required to match the large day-night temperature contrast. In fact, the effect of metallicity on the phase curve is different from cooler gas-giant counterparts, due to the high-temperature chemistry in WASP-18b's atmosphere. Additionally, we compare the current UHJ spectroscopic phase curves, WASP-18b and WASP-103b, and show that these two planets provide a consistent picture with remarkable similarities in their measured and inferred properties. However, key differences in these properties, such as their brightness offsets and radius anomalies, suggest that UHJ could be used to separate between competing theories for the inflation of gas-giant planets., Comment: Accepted for publication in A&A

The L/T transition is an important evolutionary phase in brown dwarf atmospheres, providing us with a unique opportunity to explore the effects of clouds, convection, winds, gravity and metallicity across a very narrow temperature range. Understanding these physical processes is critical for understanding ultracool atmospheres. In the next decade, we will answer three key questions regarding L/T transition atmospheres: 1. What is the physical mechanism behind the L/T transition? 2. What is the spatial extent of atmospheric structures at the L/T transition? 3. How do gravity and metallicity affect the L/T transition? The theory and methods developed for brown dwarfs will be used in the 2030s and beyond for solar-system age giant exoplanets and eventually habitable zone earth analogues. Developing these techniques now are crucial.

Recent studies have shown that ocean dynamics can have a significant warming effect on the permanent night sides of 1 to 1 tidally locked terrestrial exoplanets with Earth-like atmospheres and oceans in the middle of the habitable zone. However, the impact of ocean dynamics on the habitable zone's boundaries (inner edge and outer edge) is still unknown and represents a major gap in our understanding of this type of planets. Here we use a coupled atmosphere-ocean global climate model to show that planetary heat transport from the day to night side is dominated by the ocean at lower stellar fluxes and by the atmosphere near the inner edge of the habitable zone. This decrease in oceanic heat transport (OHT) at high stellar fluxes is mainly due to weakening of surface wind stress and a decrease in surface shortwave energy deposition. We further show that ocean dynamics have almost no effect on the observational thermal phase curves of planets near the inner edge of the habitable zone. For planets in the habitable zone's middle range, ocean dynamics moves the hottest spot on the surface eastward from the substellar point. These results suggest that future studies of the inner edge may devote computational resources to atmosphere-only processes such as clouds and radiation. For studies of the middle range and outer edge of the habitable zone, however, fully coupled ocean-atmosphere modeling will be necessary. Note that due to computational resource limitations, only one rotation period (60 Earth days) has been systematically examined in this study; future work varying rotation period as well as other parameters such as atmospheric mass and composition is required., Comment: 34 pages, 13 figures, and 1 table

Growing observational evidence has suggested active meteorology in atmospheres of brown dwarfs (BDs) and directly imaged extrasolar giant planets (EGPs). In particular, a number of surveys have shown that near-IR brightness variability is common among L and T dwarfs. Despite initial understandings of atmospheric dynamics which is the major cause of the variability by previous studies, the detailed mechanism of variability remains elusive, and we need to seek a natural, self-consistent mechanism. Clouds are important in shaping the thermal structure and spectral properties of these atmospheres via large opacity, and we expect the same for inducing atmospheric variability. In this work, using a time-dependent one-dimensional model that incorporates a self-consistent coupling between the thermal structure, convective mixing, cloud radiative heating/cooling and condensation/evaporation of clouds, we show that radiative cloud feedback can drive spontaneous atmospheric variability in both temperature and cloud structure in conditions appropriate for BDs and directly imaged EGPs. The typical periods of variability are one to tens of hours with typical amplitude of the variability up to hundreds of Kelvins in effective temperature. The existence of variability is robust over a wide range of parameter space, but the detailed evolution of variability is sensitive to model parameters. Our novel, self-consistent mechanism has important implications for the observed flux variability of BDs and directly imaged EGPs, especially those evolve in short timescales. It is also a promising mechanism for cloud breaking, which has been proposed to explain the L/T transition of brown dwarfs., Comment: 19 pages, 14 figures, accepted in ApJ

We observed two full orbital phase curves of the transiting brown dwarf KELT-1b, at 3.6um and 4.5um, using the Spitzer Space Telescope. Combined with previous eclipse data from Beatty et al. (2014), we strongly detect KELT-1b's phase variation as a single sinusoid in both bands, with amplitudes of $964\pm36$ ppm at 3.6um and $979\pm54$ ppm at 4.5um, and confirm the secondary eclipse depths measured by Beatty et al. (2014). We also measure noticeable Eastward hotspot offsets of $28.4\pm3.5$ degrees at 3.6um and $18.6\pm5.2$ degrees at 4.5um. Both the day-night temperature contrasts and the hotspot offsets we measure are in line with the trends seen in hot Jupiters (e.g., Crossfield 2015), though we disagree with the recent suggestion of an offset trend by Zhang et al. (2018). Using an ensemble analysis of Spitzer phase curves, we argue that nightside clouds are playing a noticeable role in modulating the thermal emission from these objects, based on: 1) the lack of a clear trend in phase offsets with equilibrium temperature, 2) the sharp day-night transitions required to have non-negative intensity maps, which also resolves the inversion issues raised by Keating & Cowan (2017), 3) the fact that all the nightsides of these objects appear to be at roughly the same temperature of 1000K, while the dayside temperatures increase linearly with equilibrium temperature, and 4) the trajectories of these objects on a Spitzer color-magnitude diagram, which suggest colors only explainable via nightside clouds., Comment: AJ in press. Updated to reflect the accepted version

In Zhang $\&$ Showman (2018, hereafter Paper I), we developed an analytical theory of 1D eddy diffusivity $K_{zz}$ for global-mean vertical tracer transport in a 3D atmosphere. We also presented 2D numerical simulations on fast-rotating planets to validate our theory. On a slowly rotating planet such as Venus or a tidally locked planet (not necessarily a slow-rotator) such as a hot Jupiter, the tracer distribution could exhibit significant longitudinal inhomogeneity and tracer transport is intrinsically 3D. Here we study the global-mean vertical tracer transport on tidally locked planets using 3D tracer-transport simulations. We find that our analytical $K_{zz}$ theory in Paper I is validated on tidally locked planets over a wide parameter space. $K_{zz}$ strongly depends on the large-scale circulation strength, horizontal mixing due to eddies and waves and local tracer sources and sinks due to chemistry and microphysics. As our analytical theory predicted, $K_{zz}$ on tidally locked planets also exhibit three regimes In Regime I where the chemical and microphysical processes are uniformly distributed across the globe, different chemical species should be transported via different eddy diffusivity. In Regime II where the chemical and microphysical processes are non-uniform---for example, photochemistry or cloud formation that exhibits strong day-night contrast---the global-mean vertical tracer mixing does not always behave diffusively. In the third regime where the tracer is long-lived, non-diffusive effects are significant. Using species-dependent eddy diffusivity, we provide a new analytical theory of the dynamical quench points for disequilibrium tracers on tidally locked planets from first principles., Comment: Accepted at ApJ, 16 pages, 12 figures. This is the part II. Part I is "Global-mean Vertical Tracer Mixing in Planetary Atmospheres I: Theory and Fast-rotating Planets"

On hot Jupiter exoplanets, strong horizontal and vertical winds should homogenize the abundances of the important absorbers CH$_4$ and CO much faster than chemical reactions restore chemical equilibrium. This effect, typically neglected in general circulation models (GCMs), has been suggested as explanation for discrepancies between observed infrared lightcurves 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 $\mu$m band. We modified the SPARC/MITgcm to include disequilibrium abundances of CH$_4$, CO and H$_2$O by assuming that the CH$_4$/CO ratio is constant throughout the simulation domain. We ran simulations of hot Jupiter HD 189733b with 8 CH$_4$/CO ratios. In the more likely CO-dominated regime, we find temperature changes $\geq$50-100 K compared to the equilibrium chemistry case 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 600K and 1300K. We find that spectra in regions with strong methane absorption, including the Spitzer 3.6 and 8 $\mu$m bands, are strongly impacted by disequilibrium abundances. We expect chemical quenching to result in much larger nightside fluxes in the 3.6 $\mu$m band, in stark contrast to observations. Meanwhile, we find almost no effect on predicted observations in the 4.5 $\mu$m band, as the opacity changes due to CO and H$_2$O offset each other. We thus conclude that disequilibrium carbon chemistry cannot explain the observed low nightside fluxes in the 4.5 $\mu$m band., Comment: 28 pages, 12 figures; accepted for publication in ApJ. Key changes since v1: updated introduction and discussion to include recent literature more thoroughly; added paragraph and figure to clarify change in water abundance; ran simulations for additional 600 days simulation time (leading to small quantitative changes in some figures but no qualitative changes)

Brown dwarfs and directly imaged giant planets exhibit significant evidence for active atmospheric circulation, which induces a large-scale patchiness in the cloud structure that evolves significantly over time, as evidenced by infrared light curves and Doppler maps. These observations raise critical questions about the fundamental nature of the circulation, its time variability, and the overall relationship to the circulation on Jupiter and Saturn. Jupiter and Saturn themselves exhibit numerous robust zonal (east-west) jet streams at the cloud level; moreover, both planets exhibit long-term stratospheric oscillations involving perturbations of zonal wind and temperature that propagate downward over time on timescales of ~4 years (Jupiter) and ~15 years (Saturn). These oscillations, dubbed the Quasi Quadrennial Oscillation (QQO) for Jupiter and the Semi-Annual Oscillation (SAO) on Saturn, are thought to be analogous to the Quasi-Biennial Oscillation (QBO) on Earth, which is driven by upward propagation of equatorial waves from the troposphere. To investigate these issues, we here present global, three-dimensional, high-resolution numerical simulations of the flow in the stratified atmosphere--overlying the convective interior--of brown dwarfs and Jupiter-like planets. The effect of interior convection is parameterized by inducing small-scale, randomly varying perturbations in the radiative-convective boundary at the base of the model. In the simulations, the convective perturbations generate atmospheric waves and turbulence that interact with the rotation to produce numerous zonal jets. Moreover, the equatorial stratosphere exhibits stacked eastward and westward jets that migrate downward over time, exactly as occurs in the terrestrial QBO, Jovian QQO, and Saturnian SAO. This is the first demonstration of a QBO-like phenomenon in 3D numerical simulations of a giant planet., Comment: 27 pages, 15 figures, in press at ApJ; this is the revised (accepted) version, which includes a major new section providing detailed analysis of the types of wave modes present in the model, and characterizing the wave-mean-flow interactions by which they generate the QBO-like oscillations

A new class of exoplanets has emerged: the ultra hot Jupiters, the hottest close-in gas giants. Most of them have weaker than expected spectral features in the $1.1-1.7\mu m$ bandpass probed by HST/WFC3 but stronger spectral features at longer wavelengths probed by Spitzer. This led previous authors to puzzling conclusions about the thermal structures and chemical abundances of these planets. Using the SPARC/MITgcm, we investigate how thermal dissociation, ionization, H$^-$ opacity and clouds shape the thermal structures and spectral properties of ultra hot Jupiters with a special focus on WASP-121b. We expand our findings to the whole population of ultra hot Jupiters through analytical quantification of the thermal dissociation and its influence on the strength of spectral features. We predict that most molecules are thermally dissociated and alkalies are ionized in the dayside photospheres of ultra hot Jupiters. This includes H$_{\rm 2}$O, TiO, VO, and H$_{\rm 2}$ but not CO, which has a stronger molecular bond. The vertical molecular gradient created by the dissociation significantly weakens the spectral features from water while the $4.5\mu m$ CO feature remain unchanged. The water band in the HST/WFC3 bandpass is further weakened by H$^-$ continuum opacity. Molecules are expected to recombine before reaching the limb, leading to order of magnitude variations of the chemical composition and cloud coverage between the limb and the dayside. Overall, molecular dissociation provides a qualitative understanding of the lack of strong spectral feature of water in the $1-2\mu m$ bandpass observed in most ultra hot Jupiters. Quantitatively, however, our model does not provide a satisfactory match to the WASP-121b emission spectrum. Together with WASP-33b and Kepler-33Ab, they seem the outliers among the population of ultra hot Jupiters in need of a more thorough understanding., Comment: Accepted in A&A

We present thermal phase curve measurements for the hot Jupiter WASP-103b observed with Hubble/WFC3 and Spitzer/IRAC. The phase curves have large amplitudes and negligible hotspot offsets, indicative of poor heat redistribution to the nightside. We fit the phase variation with a range of climate maps and find that a spherical harmonics model generally provides the best fit. The phase-resolved spectra are consistent with blackbodies in the WFC3 bandpass, with brightness temperatures ranging from $1880\pm40$ K on the nightside to $2930 \pm 40$ K on the dayside. The dayside spectrum has a significantly higher brightness temperature in the Spitzer bands, likely due to CO emission and a thermal inversion. The inversion is not present on the nightside. We retrieved the atmospheric composition and found the composition is moderately metal-enriched ($\mathrm{[M/H]} = 23^{+29}_{-13}\times$ solar) and the carbon-to-oxygen ratio is below 0.9 at $3\,\sigma$ confidence. In contrast to cooler hot Jupiters, we do not detect spectral features from water, which we attribute to partial H$_2$O dissociation. We compare the phase curves to 3D general circulation models and find magnetic drag effects are needed to match the data. We also compare the WASP-103b spectra to brown dwarfs and young directly imaged companions and find these objects have significantly larger water features, indicating that surface gravity and irradiation environment play an important role in shaping the spectra of hot Jupiters. These results highlight the 3D structure of exoplanet atmospheres and illustrate the importance of phase curve observations for understanding their complex chemistry and physics., Comment: 25 pages, 17 figures, 7 tables; accepted to AJ

Most chemistry and cloud formation models for planetary atmospheres adopt a one-dimensional (1D) diffusion approach to approximate the global-mean vertical tracer transport. The physical underpinning of the key parameter in this framework, eddy diffusivity $K_{zz}$, is usually obscure. Here we analytically and numerically investigate vertical tracer transport in a 3D stratified atmosphere and predict $K_{zz}$ as a function of the large-scale circulation strength, horizontal mixing due to eddies and waves and local tracer sources and sinks. We find that $K_{zz}$ increases with tracer chemical lifetime and circulation strength but decreases with horizontal eddy mixing efficiency. We demarcated three $K_{zz}$ regimes in planetary atmospheres. In the first regime where the tracer lifetime is short compared with the transport timescale and horizontal tracer distribution under chemical equilibrium ($\chi_0$) is uniformly distributed across the globe, global-mean vertical tracer mixing behaves diffusively. But the traditional assumption in current 1D models that all chemical species are transported via the same eddy diffusivity generally breaks down. We show that different chemical species in a single atmosphere should in principle have different eddy diffusion profiles. In the second regime where tracer is short-lived but $\chi_0$ is non-uniformly distributed, a significant non-diffusive component might lead to a negative $K_{zz}$ under the diffusive assumption. In the third regime where the tracer is long-lived, global-mean vertical tracer transport is also largely influenced by non-diffusive effects. Numerical simulations of 2D tracer transport on fast-rotating zonally symmetric planets validate our analytical $K_{zz}$ theory over a wide parameter space., Comment: Accepted at ApJ, 19 pages, 11 figures. Paper was split into two papers during refereeing. This is the part I

We present one of the most precise emission spectra of an exoplanet observed so far. We combine five secondary eclipses of the hot Jupiter WASP-18 b (Tday=2900K) that we secured between 1.1 and 1.7 micron with the WFC3 instrument aboard the Hubble Space Telescope. Our extracted spectrum (S/N=50, R=40) does not exhibit clearly identifiable molecular features but is poorly matched by a blackbody spectrum. We complement this data with previously published Spitzer/IRAC observations of this target and interpret the combined spectrum by computing a grid of self-consistent, 1D forward models, varying the composition and energy budget. At these high temperatures, we find there are important contributions to the overall opacity from H- ions, as well as the removal of major molecules by thermal dissociation (including water), and thermal ionization of metals. These effects were omitted in previous spectral retrievals for very hot gas giants, and we argue that they must be included to properly interpret the spectra of these objects. We infer a new metallicity and C/O ratio for WASP-18 b, and find them well constrained to be solar ([M/H]=-0.01 (0.35), C/O<0.85 at 3 sigma confidence level), unlike previous work but in line with expectations for giant planets. The best fitting self-consistent temperature-pressure profiles are inverted, resulting in an emission feature at 4.5 micron seen in the Spitzer photometry. These results further strengthen the evidence that the family of very hot gas giant exoplanets commonly exhibit thermal inversions., Comment: Accepted on 26/02/2018

We present new 3.6 and 4.5 $\mu m$ Spitzer phase curves for the highly irradiated hot Jupiter WASP-33b and the unusually dense Saturn-mass planet HD 149026b. As part of this analysis, we develop a new variant of pixel level decorrelation that is effective at removing intrapixel sensitivity variations for long observations (>10 hours) where the position of the star can vary by a significant fraction of a pixel. Using this algorithm, we measure eclipse depths, phase amplitudes, and phase offsets for both planets at 3.6 $\mu m$ and 4.5 $\mu m$. We use a simple toy model to show that WASP-33b's phase offset, albedo, and heat recirculation efficiency are largely similar to those of other hot Jupiters despite its very high irradiation. On the other hand, our fits for HD 149026b prefer a very high albedo and an unusually high recirculation efficiency. We also compare our results to predictions from general circulation models, and find that while neither are a good match to the data, the discrepancies for HD 149026b are especially large. We speculate that this may be related to its high bulk metallicity, which could lead to enhanced atmospheric opacities and the formation of reflective cloud layers in localized regions of the atmosphere. We then place these two planets in a broader context by exploring relationships between the temperatures, albedos, heat transport efficiencies, and phase offsets of all planets with published thermal phase curves. We find a striking relationship between phase offset and irradiation temperature--the former drops with increasing temperature until around 3400 K, and rises thereafter. Although some aspects of this trend are mirrored in the circulation models, there are notable differences that provide important clues for future modeling efforts.

The NASA Exoplanet Exploration Program's SAG15 group has solicited, collected, and organized community input on high-level science questions that could be addressed with future direct imaging exoplanet missions and the type and quality of data answering these questions will require. Input was solicited through a variety of forums and the report draft was shared with the exoplanet community continuously during the period of the report development (Nov 2015 -- May 2017). The report benefitted from the input of over 50 exoplanet scientists and from multiple open-forum discussions at exoplanet and astrobiology meetings. The SAG15 team has identified three group of high-level questions, those that focus on the properties of planetary systems (Questions A1--A2), those that focus on the properties of individual planets (Questions B1--B4), and questions that relate to planetary processes (Questions C1--C4). The questions in categories A, B, and C require different target samples and often different observational approaches. For each questions we summarize the current body of knowledge, the available and future observational approaches that can directly or indirectly contribute to answering the question, and provide examples and general considerations for the target sample required., Comment: For the report with higher-resolution images and for updates and supporting documentations please visit http://eos-nexus.org/sag15/

Growing observations of brown dwarfs have provided evidence for strong atmospheric circulation on these objects. Directly imaged planets share similar observations, and can be viewed as low-gravity versions of brown dwarfs. Vigorous condensate cycles of chemical species in their atmospheres are inferred by observations and theoretical studies, and latent heating associated with condensation is expected to be important in shaping atmospheric circulation and influencing cloud patchiness. We present a qualitative description of the mechanisms by which condensational latent heating influence the circulation, and then illustrate them using an idealized general circulation model that includes a condensation cycle of silicates with latent heating and molecular weight effect due to rainout of condensate. Simulations with conditions appropriate for typical T dwarfs exhibit the development of localized storms and east-west jets. The storms are spatially inhomogeneous, evolving on timescale of hours to days and extending vertically from the condensation level to the tropopause. The fractional area of the brown dwarf covered by active storms is small. Based on a simple analytic model, we quantitatively explain the area fraction of moist plumes, and show its dependence on radiative timescale and convective available potential energy. We predict that, if latent heating dominates cloud formation processes, the fractional coverage area by clouds decreases as the spectral type goes through the L/T transition from high to lower effective temperature. This is a natural consequence of the variation of radiative timescale and convective available potential energy with spectral type., Comment: 13 pages, 8 figures, accepted for publication in ApJ

The full-phase infrared light curves of low-eccentricity hot Jupiters show a trend of increasing fractional dayside-nightside brightness temperature difference with increasing incident stellar flux, both averaged across the infrared and in each individual wavelength band. The analytic theory of Komacek & Showman (2016) shows that this trend is due to the decreasing ability with increasing incident stellar flux of waves to propagate from day to night and erase temperature differences. Here, we compare the predictions of this theory to observations, showing that it explains well the shape of the trend of increasing dayside-nightside temperature difference with increasing equilibrium temperature. Applied to individual planets, the theory matches well with observations at high equilibrium temperatures but, for a fixed photosphere pressure of $100 \ \mathrm{mbar}$, systematically under-predicts the dayside-nightside brightness temperature differences at equilibrium temperatures less than $2000 \ \mathrm{K}$. We interpret this as due to as the effects of a process that moves the infrared photospheres of these cooler hot Jupiters to lower pressures. We also utilize general circulation modeling with double-grey radiative transfer to explore how the circulation changes with equilibrium temperature and drag strengths. As expected from our theory, the dayside-nightside temperature differences from our numerical simulations increase with increasing incident stellar flux and drag strengths. We calculate model phase curves using our general circulation models, from which we compare the broadband infrared offset from the substellar point and dayside-nightside brightness temperature differences against observations, finding that strong drag or additional effects (e.g. clouds and/or supersolar metallicities) are necessary to explain many observed phase curves., Comment: Accepted at ApJ, 16 pages, 11 figures

Previous measurements of heat redistribution efficiency (the ability to transport energy from a planet's highly irradiated dayside to its eternally dark nightside) show considerable variation between exoplanets. Theoretical models predict a positive correlation between heat redistribution efficiency and temperature for tidally locked planets; however, recent HST WASP-43b spectroscopic phase curve results are inconsistent with current predictions. Using the Spitzer Space Telescope, we obtained a total of three phase curve observations of WASP-43b (P=0.813 days) at 3.6 and 4.5 microns. The first 3.6 micron visit exhibits spurious nightside emission that requires invoking unphysical conditions in our cloud-free atmospheric retrievals. The two other visits exhibit strong day-night contrasts that are consistent with the HST data. To reconcile the departure from theoretical predictions, WASP-43b would need to have a high-altitude, nightside cloud/haze layer blocking its thermal emission. Clouds/hazes could be produced within the planet's cool, nearly retrograde mid-latitude flows before dispersing across its nightside at high altitudes. Since mid-latitude flows only materialize in fast-rotating ($\lesssim1$ day) planets, this may explain an observed trend connecting measured day-night contrast with planet rotation rate that matches all current Spitzer phase curve results. Combining independent planetary emission measurements from multiple phases, we obtain a precise dayside hemisphere H2O abundance ($2.5\times 10^{-5} - 1.1\times 10^{-4}$ at 1$\sigma$ confidence) and, assuming chemical equilibrium and a scaled solar abundance pattern, we derive a corresponding metallicity estimate that is consistent with being solar (0.4 -- 1.7). Using the retrieved global CO+CO2 abundance under the same assumptions, we estimate a comparable metallicity of 0.3 - 1.7$\times$ solar. (Abridged), Comment: Accepted for publication in AJ

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., Comment: 25 pages, 22 figures

The nature of the vertical thermal structure of hot Jupiter atmospheres is one of the key questions raised by the characterization of transiting exoplanets over the last decade. There have been claims that many hot Jupiter's exhibit vertical profiles with increasing temperature with decreasing pressure in the infrared photosphere that leads to the reversal of molecular absorption bands into emission features (an inversion). However, these claims have been based on broadband photometry rather than the unambiguous identification of emission features with spectroscopy, and the chemical species that could cause the thermal inversions by absorbing stellar irradiation at high altitudes have not been identified despite extensive theoretical and observational effort. Here we present high precision HST WFC3 observations of the dayside emission spectrum of the hot Jupiter HD209458b; the first exoplanet suggested to have a thermal inversion. Our observations resolve a water band in absorption at 6.2 sigma confidence. When combined with Spitzer photometry the data are indicative of a monotonically decreasing temperature with pressure over the range 1-0.001 bar at 7.7 sigma confidence. We test the robustness of our results by exploring a variety of model assumptions including the temperature profile parameterization, presence of a cloud, and choice of Spitzer data reduction. We also introduce a new analysis method, "chemical retrieval-on-retrieval", to determine the elemental abundances from the spectrally retrieved mixing ratios with thermochemical self-consistency and find plausible abundances consistent with solar metallicity (0.06 - 10 x solar) and carbon-to-oxygen ratios less than unity. This work suggests that high-precision spectrophotometric results are required to robustly infer thermal structures and compositions of extra-solar planet atmospheres., Comment: submitted to ApJ

We present $Spitzer$/IRAC Ch1 and Ch2 monitoring of six brown dwarfs during 8 different epochs over the course of 20 months. For four brown dwarfs, we also obtained simulataneous $HST$/WFC3 G141 Grism spectra during two epochs and derived light curves in five narrow-band filters. Probing different pressure levels in the atmospheres, the multi-wavelength light curves of our six targets all exhibit variations, and the shape of the light curves evolves over the timescale of a rotation period, ranging from 1.4 h to 13 h. We compare the shapes of the light curves and estimate the phase shifts between the light curves observed at different wavelengths by comparing the phase of the primary Fourier components. We use state-of-the-art atmosphere models to determine the flux contribution of different pressure layers to the observed flux in each filter. We find that the light curves that probe higher pressures are similar and in phase, but are offset and often different from the light curves that probe lower pressures. The phase differences between the two groups of light curves suggest that the modulations seen at lower and higher pressures may be introduced by different cloud layers., Comment: 34 pages, 22 figures, accepted for publication in ApJ

We present results from an atmospheric circulation study of nine hot Jupiters that comprise a large transmission spectral survey using the Hubble and Spitzer Space Telescopes. These observations exhibit a range of spectral behavior over optical and infrared wavelengths which suggest diverse cloud and haze properties in their atmospheres. By utilizing the specific system parameters for each planet, we naturally probe a wide phase space in planet radius, gravity, orbital period, and equilibrium temperature. First, we show that our model "grid" recovers trends shown in traditional parametric studies of hot Jupiters, particularly equatorial superrotation and increased day-night temperature contrast with increasing equilibrium temperature. We show how spatial temperature variations, particularly between the dayside and nightside and west and east terminators, can vary by hundreds of K, which could imply large variations in Na, K, CO and CH4 abundances in those regions. These chemical variations can be large enough to be observed in transmission with high-resolution spectrographs, such as ESPRESSO on VLT, METIS on the E-ELT, or with MIRI and NIRSpec aboard JWST. We also compare theoretical emission spectra generated from our models to available Spitzer eclipse depths for each planet, and find that the outputs from our solar-metallicity, cloud-free models generally provide a good match to many of the datasets, even without additional model tuning. Although these models are cloud-free, we can use their results to understand the chemistry and dynamics that drive cloud formation in their atmospheres., Comment: 16 pages, 10 figures; accepted to the Astrophysical Journal

Over a large range of equilibrium temperatures, clouds shape the transmission spectrum of hot Jupiter atmospheres, yet their composition remains unknown. Recent observations show that the Kepler lightcurves of some hot Jupiters are asymmetric: for the hottest planets, the lightcurve peaks before secondary eclipse, whereas for planets cooler than $\sim1900\rm\,K$, it peaks after secondary eclipse. We use the thermal structure from 3D global circulation models to determine the expected cloud distribution and Kepler lightcurves of hot Jupiters. We demonstrate that the change from an optical lightcurve dominated by thermal emission to one dominated by scattering (reflection) naturally explains the observed trend from negative to positive offset. For the cool planets the presence of an asymmetry in the Kepler lightcurve is a telltale sign of the cloud composition, because each cloud species can produce an offset only over a narrow range of effective temperatures. By comparing our models and the observations, we show that the cloud composition of hot Jupiters likely varies with equilibrium temperature. We suggest that a transition occurs between silicate and manganese sulfide clouds at a temperature near $1600\rm\,K$, analogous to the L/T transition on brown dwarfs. The cold trapping of cloud species below the photosphere naturally produces such a transition and predicts similar transitions for other condensates, including TiO. We predict that most hot Jupiters should have cloudy nightsides, that partial cloudiness should be common at the limb and that the dayside hot spot should often be cloud-free., Comment: Accepted for publication in ApJ. Revised version include effect of metallicity, CaTiO3 clouds, and models with an updated Al2O3 condensation curve. Conclusions are unchanged

We present HST STIS optical transmission spectroscopy of the cool Saturn-mass exoplanet WASP-39b from 0.29-1.025 micron, along with complementary transit observations from Spitzer IRAC at 3.6 and 4.5 micron. The low density and large atmospheric pressure scale height ofWASP-39b make it particularly amenable to atmospheric characterization using this technique. We detect a Rayleigh scattering slope as well as sodium and potassium absorption features; this is the first exoplanet in which both alkali features are clearly detected with the extended wings predicted by cloud-free atmosphere models. The full transmission spectrum is well matched by a clear, H2-dominated atmosphere or one containing a weak contribution from haze, in good agreement with the preliminary reduction of these data presented in Sing et al. (2016). WASP-39b is predicted to have a pressure-temperature profile comparable to that of HD 189733b and WASP-6b, making it one of the coolest transiting gas giants observed in our HST STIS survey. Despite this similarity, WASP-39b appears to be largely cloud-free while the transmission spectra of HD 189733b and WASP-6b both indicate the presence of high altitude clouds or hazes. These observations further emphasize the surprising diversity of cloudy and cloud-free gas giant planets in short-period orbits and the corresponding challenges associated with developing predictive cloud models for these atmospheres.

The full-phase infrared light curves of low-eccentricity hot Jupiters show a trend of increasing dayside-to-nightside brightness temperature difference with increasing equilibrium temperature. Here we present a three-dimensional model that explains this relationship, in order to shed insight on the processes that control heat redistribution in tidally-locked planetary atmospheres. This three-dimensional model combines predictive analytic theory for the atmospheric circulation and dayside-nightside temperature differences over a range of equilibrium temperature, atmospheric composition, and potential frictional drag strengths with numerical solutions of the circulation that verify this analytic theory. This analytic theory shows that the longitudinal propagation of waves mediates dayside-nightside temperature differences in hot Jupiter atmospheres, analogous to the wave adjustment mechanism that regulates the thermal structure in Earth's tropics. These waves can be damped in hot Jupiter atmospheres by either radiative cooling or potential frictional drag. This frictional drag would likely be caused by Lorentz forces in a partially ionized atmosphere threaded by a background magnetic field, and would increase in strength with increasing temperature. Additionally, the amplitude of radiative heating and cooling increases with increasing temperature, and hence both radiative heating/cooling and frictional drag damp waves more efficiently with increasing equilibrium temperature. Radiative heating and cooling play the largest role in controlling dayside-nightside temperature temperature differences in both our analytic theory and numerical simulations, with frictional drag only important if it is stronger than the Coriolis force. As a result, dayside-nightside temperature differences in hot Jupiter atmospheres increase with increasing stellar irradiation and decrease with increasing pressure., Comment: 24 pages, 14 figures, Accepted to ApJ

We analyze full-orbit phase curve observations of the transiting hot Jupiters WASP-19b and HAT-P-7b at 3.6 and 4.5 $\mu$m obtained using the Spitzer Space Telescope. For WASP-19b, we measure secondary eclipse depths of $0.485\%\pm 0.024\%$ and $0.584\%\pm 0.029\%$ at 3.6 and 4.5 $\mu$m, which are consistent with a single blackbody with effective temperature $2372 \pm 60$ K. The measured 3.6 and 4.5 $\mu$m secondary eclipse depths for HAT-P-7b are $0.156\%\pm 0.009\%$ and $0.190\%\pm 0.006\%$, which are well-described by a single blackbody with effective temperature $2667\pm 57$ K. Comparing the phase curves to the predictions of one-dimensional and three-dimensional atmospheric models, we find that WASP-19b's dayside emission is consistent with a model atmosphere with no dayside thermal inversion and moderately efficient day-night circulation. We also detect an eastward-shifted hotspot, suggesting the presence of a superrotating equatorial jet. In contrast, HAT-P-7b's dayside emission suggests a dayside thermal inversion and relatively inefficient day-night circulation; no hotspot shift is detected. For both planets, these same models do not agree with the measured nightside emission. The discrepancies in the model-data comparisons for WASP-19b might be explained by high-altitude silicate clouds on the nightside and/or high atmospheric metallicity, while the very low 3.6 $\mu$m nightside planetary brightness for HAT-P-7b may be indicative of an enhanced global C/O ratio. We compute Bond albedos of 0 ($<0.08$ at $1\sigma$) and $0.38\pm 0.06$ for WASP-19b and HAT-P-7b, respectively. In the context of other planets with thermal phase curve measurements, we show that WASP-19b and HAT-P-7b fit the general trend of decreasing day-night heat recirculation with increasing irradiation., Comment: 22 pages, 29 figures, accepted by ApJ

Thousands of transiting exoplanets have been discovered, but spectral analysis of their atmospheres has so far been dominated by a small number of exoplanets and data spanning relatively narrow wavelength ranges (such as 1.1 to 1.7 {\mu}m). Recent studies show that some hot-Jupiter exoplanets have much weaker water absorption features in their near-infrared spectra than predicted. The low amplitude of water signatures could be explained by very low water abundances, which may be a sign that water was depleted in the protoplanetary disk at the planet's formation location, but it is unclear whether this level of depletion can actually occur. Alternatively, these weak signals could be the result of obscuration by clouds or hazes, as found in some optical spectra. Here we report results from a comparative study of ten hot Jupiters covering the wavelength range 0.3-5 micrometres, which allows us to resolve both the optical scattering and infrared molecular absorption spectroscopically. Our results reveal a diverse group of hot Jupiters that exhibit a continuum from clear to cloudy atmospheres. We find that the difference between the planetary radius measured at optical and infrared wavelengths is an effective metric for distinguishing different atmosphere types. The difference correlates with the spectral strength of water, so that strong water absorption lines are seen in clear-atmosphere planets and the weakest features are associated with clouds and hazes. This result strongly suggests that primordial water depletion during formation is unlikely and that clouds and hazes are the cause of weaker spectral signatures., Comment: This is the authors version of the manuscript, 18 pages including Methods. Published in Nature, available at http://nature.com/articles/doi:10.1038/nature16068 spectra also available at http://www.astro.ex.ac.uk/people/sing/David_Sing/Spectra.html

Rotational modulations of brown dwarfs have recently provided powerful constraints on the properties of ultra-cool atmospheres, including longitudinal and vertical cloud structures and cloud evolution. Furthermore, periodic light curves directly probe the rotational periods of ultra-cool objects. We present here, for the first time, time-resolved high-precision photometric measurements of a planetary-mass companion, 2M1207b. We observed the binary system with HST/WFC3 in two bands and with two spacecraft roll angles. Using point spread function-based photometry, we reach a nearly photon-noise limited accuracy for both the primary and the secondary. While the primary is consistent with a flat light curve, the secondary shows modulations that are clearly detected in the combined light curve as well as in different subsets of the data. The amplitudes are 1.36% in the F125W and 0.78% in the F160W filters, respectively. By fitting sine waves to the light curves, we find a consistent period of $10.7^{+1.2}_{-0.6}$ hours and similar phases in both bands. The J- and H-band amplitude ratio of 2M1207b is very similar to a field brown dwarf that has identical spectral type but different J-H color. Importantly, our study also measures, for the first time, the rotation period for a directly imaged extra-solar planetary-mass companion., Comment: Accepted for publication in ApJ. 13 pages, 5 figures, 2 tables