Observations of brown dwarfs (BDs), free-floating planetary-mass objects, and directly imaged extrasolar giant planets (EGPs) exhibit rich evidence of large-scale weather. Understanding the mechanisms driving the vigorous atmospheric circulation of BDs and directly imaged EGPs and its effects on their observed lightcurve variability and spectral properties is a pressing need. Our previous work has suggested a strong role of radiative cloud feedback on inducing a spontaneous time evolution in a simple one-dimensional framework. Yet the radiative cloud feedback in a three-dimensional (3D) dynamical framework remains unexplored for conditions relevant to these objects. Here we present a series of atmospheric circulation models that self-consistently couple 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. The circulation is dominated by cloud-forming and clear-sky vortices that evolve over timescales from several to tens of hours. The typical horizontal lengthscales of dominant vortices are strongly constrained by the rotation, showing a linear dependence on the inverse of rotation rate with stronger rotation leading to thinner clouds. Domain-mean outgoing radiative flux exhibits variability over timescales of tens of hours due to the statistical evolution of storms. The circulation driven by cloud radiative feedback represents a natural mechanism generating significant surface inhomogeneity as well as irregular flux variability. Our results also 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.

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.

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.

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.

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]

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Observations of the highly-eccentric (e~0.9) hot-Jupiter HD 80606b with Spitzer have provided some of best probes of the physics at work in exoplanet atmospheres. By observing HD 80606b during its periapse passage, atmospheric radiative, advective, and chemical timescales can be directly measured and used to constrain fundamental planetary properties such as rotation period, tidal dissipation rate, and atmospheric composition (including aerosols). Here we present three-dimensional general circulation models for HD 80606b that aim to further explore the atmospheric physics shaping HD 80606b's observed Spitzer phase curves. We find that our models that assume a planetary rotation period twice that of the pseudo-synchronous rotation period best reproduce the phase variations observed for HD~80606b near periapse passage with Spitzer. Additionally, we find that the rapid formation/dissipation and vertical transport of clouds in HD 80606b's atmosphere near periapse passage likely shapes its observed phase variations. We predict that observations near periapse passage at visible wavelengths could constrain the composition and formation/advection timescales of the dominant cloud species in HD 80606b's atmosphere. The time-variable forcing experienced by exoplanets on eccentric orbits provides a unique and important window on radiative, dynamical, and chemical processes in planetary atmospheres and an important link between exoplanet observations and theory.

Extrasolar planets on eccentric short-period orbits provide a laboratory in which to study radiative and tidal interactions between a planet and its host star under extreme forcing conditions. Studying such systems probes how the planet's atmosphere redistributes the time-varying heat flux from its host and how the host star responds to transient tidal distortion. Here, we report the insights into the planet-star interactions in HAT-P-2's eccentric planetary system gained from the analysis of 350 hr of 4.5 micron observations with the Spitzer Space Telescope. The observations show no sign of orbit-to-orbit variability nor of orbital evolution of the eccentric planetary companion, HAT-P-2 b. The extensive coverage allows us to better differentiate instrumental systematics from the transient heating of HAT-P-2 b's 4.5 micron photosphere and yields the detection of stellar pulsations with an amplitude of approximately 40 ppm. These pulsation modes correspond to exact harmonics of the planet's orbital frequency, indicative of a tidal origin. Transient tidal effects can excite pulsation modes in the envelope of a star, but, to date, such pulsations had only been detected in highly eccentric stellar binaries. Current stellar models are unable to reproduce HAT-P-2's pulsations, suggesting that our understanding of the interactions at play in this system is incomplete.

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.

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.

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)

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Detailed characterization of exoplanets has begun to yield measurements of their atmospheric properties that constrain the planets' origins and evolution. For example, past observations of the dayside emission spectrum of the hot Jupiter WASP-12b indicated that its atmosphere has a high carbon-to-oxygen ratio (C/O $>$ 1), suggesting it had a different formation pathway than is commonly assumed for giant planets. Here we report a precise near-infrared transmission spectrum for WASP-12b based on six transit observations with the Hubble Space Telescope/Wide Field Camera 3. We bin the data in 13 spectrophotometric light curves from 0.84 - 1.67 $\mu$m and measure the transit depths to a median precision of 51 ppm. We retrieve the atmospheric properties using the transmission spectrum and find strong evidence for water absorption (7$\sigma$ confidence). This detection marks the first high-confidence, spectroscopic identification of a molecule in the atmosphere of WASP-12b. The retrieved 1$\sigma$ water volume mixing ratio is between $10^{-5}-10^{-2}$, which is consistent with C/O $>$ 1 to within 2$\sigma$. However, we also introduce a new retrieval parameterization that fits for C/O and metallicity under the assumption of chemical equilibrium. With this approach, we constrain C/O to $0.5^{+0.2}_{-0.3}$ at $1\,\sigma$ and rule out a carbon-rich atmosphere composition (C/O$>1$) at $>3\sigma$ confidence. Further observations and modeling of the planet's global thermal structure and dynamics would aid in resolving the tension between our inferred C/O and previous constraints. Our findings highlight the importance of obtaining high-precision data with multiple observing techniques in order to obtain robust constraints on the chemistry and physics of exoplanet atmospheres.

We present full-orbit phase curve observations of the eccentric ($e\sim 0.08$) transiting hot Jupiter WASP-14b obtained in the 3.6 and 4.5 $\mu$m bands using the \textit{Spitzer Space Telescope}. We use two different methods for removing the intrapixel sensitivity effect and compare their efficacy in decoupling the instrumental noise. Our measured secondary eclipse depths of $0.1882\%\pm 0.0048\%$ and $0.2247\%\pm 0.0086\%$ at 3.6 and 4.5 $\mu$m, respectively, are both consistent with a blackbody temperature of $2402\pm 35$ K. We place a $2\sigma$ upper limit on the nightside flux at 3.6 $\mu$m and find it to be $9\%\pm 1\%$ of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 $\mu$m, the minimum planet flux is $30\%\pm 5\%$ of the maximum flux, corresponding to a brightness temperature of $1380\pm 65$ K. We compare our measured phase curves to the predictions of one-dimensional radiative transfer and three-dimensional general circulation models. We find that WASP-14b's measured dayside emission is consistent with a model atmosphere with equilibrium chemistry and a moderate temperature inversion. These same models tend to over-predict the nightside emission at 3.6 $\mu$m, while under-predicting the nightside emission at 4.5 $\mu$m. We propose that this discrepancy might be explained by an enhanced global C/O ratio. In addition, we find that the phase curves of WASP-14b ($7.8 M_{\mathrm{Jup}}$) are consistent with a much lower albedo than those of other Jovian mass planets with thermal phase curve measurements, suggesting that it may be emitting detectable heat from the deep atmosphere or interior processes.

In this work we present Spitzer 3.6 and 4.5 micron secondary eclipse observations of five new cool (<1200 K) transiting gas giant planets: HAT-P-19b, WASP-6b, WASP-10b, WASP-39b, and WASP-67b. We compare our measured eclipse depths to the predictions of a suite of atmosphere models and to eclipse depths for planets with previously published observations in order to constrain the temperature- and mass-dependent properties of gas giant planet atmospheres. We find that the dayside emission spectra of planets less massive than Jupiter require models with efficient circulation of energy to the night side and/or increased albedos, while those with masses greater than that of Jupiter are consistently best-matched by models with inefficient circulation and low albedos. At these relatively low temperatures we expect the atmospheric methane to CO ratio to vary as a function of metallicity, and we therefore use our observations of these planets to constrain their atmospheric metallicities. We find that the most massive planets have dayside emission spectra that are best-matched by solar metallicity atmosphere models, but we are not able to place strong constraints on metallicities of the smaller planets in our sample. Interestingly, we find that the ratio of the 3.6 and 4.5 micron brightness temperatures for these cool transiting planets is independent of planet temperature, and instead exhibits a tentative correlation with planet mass. If this trend can be confirmed, it would suggest that the shape of these planets' emission spectra depends primarily on their masses, consistent with the hypothesis that lower-mass planets are more likely to have metal-rich atmospheres.

HAT-P-20b is a giant exoplanet that orbits a metal-rich star. The planet itself has a high total density, suggesting that it may also have a high metallicity in its atmosphere. We analyze two eclipses of the planet in each of the 3.6- and 4.5 micron bands of Warm Spitzer. These data exhibit intra-pixel detector sensitivity fluctuations that were resistant to traditional decorrelation methods. We have developed a simple, powerful, and radically different method to correct the intra-pixel effect for Warm Spitzer data, which we call pixel-level decorrelation (PLD). PLD corrects the intra-pixel effect very effectively, but without explicitly using - or even measuring - the fluctuations in the apparent position of the stellar image. We illustrate and validate PLD using synthetic and real data, and comparing the results to previous analyses. PLD can significantly reduce or eliminate red noise in Spitzer secondary eclipse photometry, even for eclipses that have proven to be intractable using other methods. Our successful PLD analysis of four HAT-P-20b eclipses shows a best-fit blackbody temperature of 1134 +/-29K, indicating inefficient longitudinal transfer of heat, but lacking evidence for strong molecular absorption. We find sufficient evidence for variability in the 4.5 micron band that the eclipses should be monitored at that wavelength by Spitzer, and this planet should be a high priority for JWST spectroscopy. All four eclipses occur about 35 minutes after orbital phase 0.5, indicating a slightly eccentric orbit. A joint fit of the eclipse and transit times with extant RV data yields e(cos{omega}) = 0.01352 (+0.00054, -0.00057), and establishes the small eccentricity of the orbit to high statistical confidence. Given the existence of a bound stellar companion, HAT-P-20b is another excellent candidate for orbital evolution via Kozai migration or other three-body mechanism.

The recent discoveries of terrestrial exoplanets and super-Earths extending over a broad range of orbital and physical parameters suggest that these planets will span a wide range of climatic regimes. Characterization of the atmospheres of warm super-Earths has already begun and will be extended to smaller and more distant planets over the coming decade. The habitability of these worlds may be strongly affected by their three-dimensional atmospheric circulation regimes, since the global climate feedbacks that control the inner and outer edges of the habitable zone including transitions to Snowball-like states and runaway-greenhouse feedbacks depend on the equator-to-pole temperature differences, patterns of relative humidity, and other aspects of the dynamics. Here, using an idealized moist atmospheric general circulation model including a hydrological cycle, we study the dynamical principles governing the atmospheric dynamics on such planets. We show how the planetary rotation rate, stellar flux, atmospheric mass, surface gravity, optical thickness, and planetary radius affect the atmospheric circulation and temperature distribution on such planets. Our simulations demonstrate that equator-to-pole temperature differences, meridional heat transport rates, structure and strength of the winds, and the hydrological cycle vary strongly with these parameters, implying that the sensitivity of the planet to global climate feedbacks will depend significantly on the atmospheric circulation. We elucidate the possible climatic regimes and diagnose the mechanisms controlling the formation of atmospheric jet streams, Hadley and Ferrel cells, and latitudinal temperature differences. Finally, we discuss the implications for understanding how the atmospheric circulation influences the global climate.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Efforts to characterize extrasolar giant planet (EGP) atmospheres have so far emphasized planets within 0.05 AU of their stars. Despite this focus, known EGPs populate a continuum of orbital separations from canonical hot Jupiter values (0.03-0.05 AU) out to 1 AU and beyond. Unlike typical hot Jupiters, these more distant EGPs will not in general be synchronously rotating. In anticipation of observations of this population, we here present three-dimensional atmospheric circulation models exploring the dynamics that emerge over a broad range of rotation rates and incident stellar fluxes appropriate for warm and hot Jupiters. We find that the circulation resides in one of two basic regimes. On typical hot Jupiters, the strong day-night heating contrast leads to a broad, fast superrotating (eastward) equatorial jet and large day-night temperature differences. At faster rotation rates and lower incident fluxes, however, the day-night heating gradient becomes less important, and baroclinic instabilities emerge as a dominant player, leading to eastward jets in the midlatitudes, minimal temperature variations in longitude, and, in many cases, weak winds at the equator. Our most rapidly rotating and least irradiated models exhibit multiple eastward jets in each hemisphere--similar to the jets on Jupiter and Saturn--and illuminate the dynamical continuum between highly irradiated EGPs and the weakly irradiated giant planets of our own Solar System. We present infrared (IR) light curves and spectra of these models, which show that the amplitude and offset of the IR phase variation, as well as the shape of the spectra, depend significantly on incident flux and rotation rate. This provides a way to identify the regime transition in future observations and suggests that, in some cases, IR light curves can provide constraints on the rotation rate of non-synchronously rotating planets.

Kepler has detected numerous exoplanet transits by precise measurements of stellar light in a single visible-wavelength band. In addition to detection, the precise photometry provides phase curves of exoplanets, which can be used to study the dynamic processes on these planets. However, the interpretation of these observations can be complicated by the fact that visible-wavelength phase curves can represent both thermal emission and scattering from the planets. Here we present a semi-analytical model framework that can be applied to study Kepler and future visible-wavelength phase curve observations of exoplanets. The model efficiently computes reflection and thermal emission components for both rocky and gaseous planets, considering both homogeneous and inhomogeneous surfaces or atmospheres. We analyze the phase curves of the gaseous planet Kepler-7 b and the rocky planet Kepler-10 b using the model. In general, we find that a hot exoplanet's visible-wavelength phase curve having a significant phase offset can usually be explained by two classes of solutions: one class requires a thermal hot spot shifted to one side of the substellar point, and the other class requires reflective clouds concentrated on the same side of the substellar point. The two solutions would require very different Bond albedos to fit the same phase curve; atmospheric circulation models or eclipse observations at longer wavelengths can effectively rule out one class of solutions, and thus pinpoint the albedo of the planet, allowing decomposition of the reflection and the thermal emission components in the phase curve. Particularly for Kepler-7 b, reflective clouds located on the west side of the substellar point can best explain its phase curve. We further derive that the reflectivity of the clear part of the atmosphere should be less than 7% and that of the cloudy part should be greater than 80% (abridged)

The hot Jupiter WASP-43b has now joined the ranks of transiting hot Jupiters HD 189733b and HD 209458b as an exoplanet with a large array of observational constraints on its atmospheric properties. Because WASP-43b receives a similar stellar flux as HD 209458b but has a rotation rate 4 times faster and a much higher gravity, studying WASP-43b serves as a test of the effect of rotation rate and gravity on the circulation when stellar irradiation is held approximately constant. Here we present 3D atmospheric circulation models of WASP-43b using the SPARC/MITgcm, a coupled radiation and circulation model, exploring the effects of composition, metallicity, and frictional drag. We find that the circulation regime of WASP-43b is not unlike other hot Jupiters, with equatorial superrotation that yields an eastward-shifted hotspot and large day-night temperature variations (~600 K at photospheric pressures). We then compare our model results to observations from Stevenson et al. which utilize HST/WFC3 to collect spectrophotometric phase curve measurements of WASP-43b from 1.12-1.65 microns. Our results show the 5x solar model lightcurve provides a good match to the data, with a phase offset of peak flux and planet/star flux ratio that is similar to observations; however, the model nightside appears to be brighter. Nevertheless, our 5x solar model provides an excellent match to the WFC3 dayside emission spectrum. This is a major success, as the result is a natural outcome of the 3D dynamics with no model tuning, and differs significantly from 1D models that can generally only match observations when appropriately tuned. In sum, these results demonstrate that 3D circulation models can provide important insights in interpreting exoplanet atmospheric observations, even at high spectral resolution, and highlight the potential for future observations with HST, JWST and other next-generation telescopes.

We report Hubble Space Telescope (HST) optical to near-infrared transmission spectroscopy of the hot Jupiter WASP-6b, measured with the Space Telescope Imaging Spectrograph (STIS) and Spitzer's InfraRed Array Camera (IRAC). The resulting spectrum covers the range $0.29-4.5\,\mu$m. We find evidence for modest stellar activity of WASP-6b and take it into account in the transmission spectrum. The overall main characteristic of the spectrum is an increasing radius as a function of decreasing wavelength corresponding to a change of $\Delta (R_p/R_{\ast})=0.0071$ from 0.33 to $4.5\,\mu$m. The spectrum suggests an effective extinction cross-section with a power law of index consistent with Rayleigh scattering, with temperatures of $973\pm144$ K at the planetary terminator. We compare the transmission spectrum with hot-Jupiter atmospheric models including condensate-free and aerosol-dominated models incorporating Mie theory. While none of the clear-atmosphere models is found to be in good agreement with the data, we find that the complete spectrum can be described by models that include significant opacity from aerosols including Fe-poor Mg$_2$SiO$_4$, MgSiO$_3$, KCl and Na$_2$S dust condensates. WASP-6b is the second planet after HD189733b which has equilibrium temperatures near $\sim1200$ K and shows prominent atmospheric scattering in the optical.

We present time-resolved near-infrared spectroscopy of two L5 dwarfs, 2MASS J18212815+1414010 and 2MASS J15074759-1627386, observed with the Wide Field Camera 3 instrument on the Hubble Space Telescope (HST). We study the wavelength dependence of rotation-modulated flux variations between 1.1 $\mu$m and 1.7 $\mu$m. We find that the water absorption bands of the two L5 dwarfs at 1.15 $\mu$m and 1.4 $\mu$m vary at similar amplitudes as the adjacent continuum. This differs from the results of previous HST observations of L/T transition dwarfs, in which the water absorption at 1.4 $\mu$m displays variations of about half of the amplitude at other wavelengths. We find that the relative amplitude of flux variability out of the water band with respect to that in the water band shows a increasing trend from the L5 dwarfs toward the early T dwarfs. We utilize the models of Saumon & Marley (2008) and find that the observed variability of the L5 dwarfs can be explained by the presence of spatially varying high-altitude haze layers above the condensate clouds. Therefore, our observations show that the heterogeneity of haze layers - the driver of the variability - must be located at very low pressures, where even the water opacity is negligible. In the near future, the rotational spectral mapping technique could be utilized for other atomic and molecular species to probe different pressure levels in the atmospheres of brown dwarfs and exoplanets and uncover both horizontal and vertical cloud structures.

We present Hubble Space Telescope optical and near-IR transmission spectra of the transiting hot-Jupiter WASP-31b. The spectrum covers 0.3-1.7 $\mu$m at a resolution $R\sim$70, which we combine with Spitzer photometry to cover the full-optical to IR. The spectrum is dominated by a cloud-deck with a flat transmission spectrum which is apparent at wavelengths $>0.52\mu$m. The cloud deck is present at high altitudes and low pressures, as it covers the majority of the expected optical Na line and near-IR H$_2$O features. While Na I absorption is not clearly identified, the resulting spectrum does show a very strong potassium feature detected at the 4.2-$\sigma$ confidence level. Broadened alkali wings are not detected, indicating pressures below $\sim$10 mbar. The lack of Na and strong K is the first indication of a sub-solar Na/K abundance ratio in a planetary atmosphere (ln[Na/K]$=-3.3\pm2.8$), which could potentially be explained by Na condensation on the planet's night side, or primordial abundance variations. A strong Rayleigh scattering signature is detected at short wavelengths, with a 4-$\sigma$ significant slope. Two distinct aerosol size populations can explain the spectra, with a smaller sub-micron size grain population reaching high altitudes producing a blue Rayleigh scattering signature on top of a larger, lower-lying population responsible for the flat cloud deck at longer wavelengths. We estimate that the atmospheric circulation is sufficiently strong to mix micron size particles upward to the required 1-10 mbar pressures, necessary to explain the cloud deck. These results further confirm the importance of clouds in hot-Jupiters, which can potentially dominate the overall spectra and may alter the abundances of key gaseous species.

The characterization of exoplanetary atmospheres has come of age in the last decade, as astronomical techniques now allow for albedos, chemical abundances, temperature profiles and maps, rotation periods and even wind speeds to be measured. Atmospheric dynamics sets the background state of density, temperature and velocity that determines or influences the spectral and temporal appearance of an exoplanetary atmosphere. Hot exoplanets are most amenable to these characterization techniques; in the present review, we focus on highly-irradiated, large exoplanets (the "hot Jupiters"), as astronomical data begin to confront theoretical questions. We summarize the basic atmospheric quantities inferred from the astronomical observations. We review the state of the art by addressing a series of current questions and look towards the future by considering a separate set of exploratory questions. Attaining the next level of understanding will require a concerted effort of constructing multi-faceted, multi-wavelength datasets for benchmark objects. Understanding clouds presents a formidable obstacle, as they introduce degeneracies into the interpretation of spectra, yet their properties and existence are directly influenced by atmospheric dynamics. Confronting general circulation models with these multi-faceted, multi-wavelength datasets will help us understand these and other degeneracies. The coming decade will witness a decisive confrontation of theory and simulation by the next generation of astronomical data.

The water abundance in a planetary atmosphere provides a key constraint on the planet's primordial origins because water ice is expected to play an important role in the core accretion model of planet formation. However, the water content of the Solar System giant planets is not well known because water is sequestered in clouds deep in their atmospheres. By contrast, short-period exoplanets have such high temperatures that their atmospheres have water in the gas phase, making it possible to measure the water abundance for these objects. We present a precise determination of the water abundance in the atmosphere of the 2 $M_\mathrm{Jup}$ short-period exoplanet WASP-43b based on thermal emission and transmission spectroscopy measurements obtained with the Hubble Space Telescope. We find the water content is consistent with the value expected in a solar composition gas at planetary temperatures (0.4-3.5x solar at 1 $\sigma$ confidence). The metallicity of WASP-43b's atmosphere suggested by this result extends the trend observed in the Solar System of lower metal enrichment for higher planet masses.

Exoplanets that orbit close to their host stars are much more highly irradiated than their Solar System counterparts. Understanding the thermal structures and appearances of these planets requires investigating how their atmospheres respond to such extreme stellar forcing. We present spectroscopic thermal emission measurements as a function of orbital phase ("phase-curve observations") for the highly-irradiated exoplanet WASP-43b spanning three full planet rotations using the Hubble Space Telescope. With these data, we construct a map of the planet's atmospheric thermal structure, from which we find large day-night temperature variations at all measured altitudes and a monotonically decreasing temperature with pressure at all longitudes. We also derive a Bond albedo of 0.18 +0.07,-0.12 and an altitude dependence in the hot-spot offset relative to the substellar point.

We report secondary eclipse photometry of the hot Jupiter HAT-P-32Ab, taken with Hale/WIRC in H and Ks bands and with Spitzer/IRAC at 3.6 and 4.5 micron. We carried out adaptive optics imaging of the planet host star HAT-P-32A and its companion HAT-P-32B in the near-IR and the visible. We clearly resolve the two stars from each other and find a separation of 2.923" +/- 0. 004" and a position angle 110.64 deg +/- 0.12 deg. We measure the flux ratios of the binary in g' r' i' z' and H & Ks bands, and determine Teff = 3565 +/- 82 K for the companion star, corresponding to an M1.5 dwarf. We use PHOENIX stellar atmosphere models to correct the dilution of the secondary eclipse depths of the hot Jupiter due to the presence of the M1.5 companion. We also improve the secondary eclipse photometry by accounting for the non-classical, flux-dependent nonlinearity of the WIRC IR detector in the H band. We measure planet-to-star flux ratios of 0.090 +/- 0.033%, 0.178 +/- 0.057%, 0.364 +/- 0.016%, and 0.438 +/- 0.020% in the H, Ks, 3.6 and 4.5 micron bands, respectively. We compare these with planetary atmospheric models, and find they prefer an atmosphere with a temperature inversion and inefficient heat redistribution. However, we also find that the data are equally well-described by a blackbody model for the planet with Tp = 2042 +/- 50 K. Finally, we measure a secondary eclipse timing offset of 0.3 +/- 1.3 min from the predicted mid-eclipse time, which constrains e = 0.0072 +0.0700/-0.0064 when combined with RV data and is more consistent with a circular orbit.

We report secondary eclipse photometry of the hot Jupiter XO-3b in the 4.5~$\mu$m band taken with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We measure individual eclipse depths and center of eclipse times for a total of twelve secondary eclipses. We fit these data simultaneously with two transits observed in the same band in order to obtain a global best-fit secondary eclipse depth of $0.1580\pm 0.0036\%$ and a center of eclipse phase of $0.67004\pm 0.00013 $. We assess the relative magnitude of variations in the dayside brightness of the planet by measuring the size of the residuals during ingress and egress from fitting the combined eclipse light curve with a uniform disk model and place an upper limit of 0.05$\%$. The new secondary eclipse observations extend the total baseline from one and a half years to nearly three years, allowing us to place an upper limit on the periastron precession rate of $2.9\times 10^{-3}$ degrees/day the tightest constraint to date on the periastron precession rate of a hot Jupiter. We use the new transit observations to calculate improved estimates for the system properties, including an updated orbital ephemeris. We also use the large number of secondary eclipses to obtain the most stringent limits to date on the orbit-to-orbit variability of an eccentric hot Jupiter and demonstrate the consistency of multiple-epoch Spitzer observations.

The hot Jupiter HD 209458b is particularly amenable to detailed study as it is among the brightest transiting exoplanet systems currently known (V-mag = 7.65; K-mag = 6.308) and has a large planet-to-star contrast ratio. HD 209458b is predicted to be in synchronous rotation about its host star with a hot spot that is shifted eastward of the substellar point by superrotating equatorial winds. Here we present the first full-orbit observations of HD 209458b, in which its 4.5 $\mu$m emission was recorded with $Spitzer$/IRAC. Our study revises the previous 4.5 $\mu$m measurement of HD 209458b's secondary eclipse emission downward by $\sim$35% to $0.1391%^{+0.0072%}_{-0.0069%}$, changing our interpretation of the properties of its dayside atmosphere. We find that the hot spot on the planet's dayside is shifted eastward of the substellar point by $40.9^{\circ}\pm{6.0^{\circ}}$, in agreement with circulation models predicting equatorial superrotation. HD 209458b's dayside (T$_{bright}$ = 1499 $\pm$ 15 K) and nightside (T$_{bright}$ = 972 $\pm$ 44 K) emission indicates a day-to-night brightness temperature contrast smaller than that observed for more highly irradiated exoplanets, suggesting that the day-to-night temperature contrast may be partially a function of the incident stellar radiation. The observed phase curve shape deviates modestly from global circulation model predictions potentially due to disequilibrium chemistry or deficiencies in the current hot CH$_{4}$ line lists used in these models. Observations of the phase curve at additional wavelengths are needed in order to determine the possible presence and spatial extent of a dayside temperature inversion, as well as to improve our overall understanding of this planet's atmospheric circulation.

The hot-Jupiter HAT-P-2b has become a prime target for Spitzer Space Telescope observations aimed at understanding the atmospheric response of exoplanets on highly eccentric orbits. Here we present a suite of three-dimensional atmospheric circulation models for HAT-P-2b that investigate the effects of assumed atmospheric composition and rotation rate on global scale winds and thermal patterns. We compare and contrast atmospheric models for HAT-P-2b, which assume one and five times solar metallicity, both with and without TiO/VO as atmospheric constituents. Additionally we compare models that assume a rotation period of half, one, and two times the nominal pseudo-synchronous rotation period. We find that changes in assumed atmospheric metallicity and rotation rate do not significantly affect model predictions of the planetary flux as a function of orbital phase. However, models in which TiO/VO are present in the atmosphere develop a transient temperature inversion between the transit and secondary eclipse events that results in significant variations in the timing and magnitude of the peak of the planetary flux compared with models in which TiO/VO are omitted from the opacity tables. We find that no one single atmospheric model can reproduce the recently observed full orbit phase curves at 3.6, 4.5 and 8.0 microns, which is likely due to a chemical process not captured by our current atmospheric models for HAT-P-2b. Further modeling and observational efforts focused on understanding the chemistry of HAT-P-2b's atmosphere are needed and could provide key insights into the interplay between radiative, dynamical, and chemical processes in a wide range of exoplanet atmospheres.

More than a thousand exoplanets have been discovered over the last decade. Perhaps more excitingly, probing their atmospheres has become possible. With current data we have glimpsed the diversity of exoplanet atmospheres that will be revealed over the coming decade. However, numerous questions concerning their chemical composition, thermal structure, and atmospheric dynamics remain to be answered. More observations of higher quality are needed. In the next years, the selection of a space-based mission dedicated to the spectroscopic characterization of exoplanets would revolutionize our understanding of the physics of planetary atmospheres. Such a mission was proposed to the ESA cosmic vision program in 2014. Our paper is therefore based on the planned capabilities of the Exoplanet Characterization Observatory (EChO), but it should equally apply to any future mission with similar characteristics. With its large spectral coverage ($4-16\, \rm{\mu m}$), high spectral resolution ($\Delta\lambda/\lambda>300$ below $5\,\rm{\mu m}$ and $\Delta\lambda/\lambda>30$ above $5\,\rm{\mu m}$) and $1.5\rm{m}$ mirror, a future mission such as EChO will provide spectrally resolved transit lightcurves, secondary eclipses lightcurves, and full phase curves of numerous exoplanets with an unprecedented signal-to-noise ratio. In this paper, we review some of today's main scientific questions about gas giant exoplanets atmospheres, for which a future mission such as EChO will bring a decisive contribution.

A variety of observational evidence demonstrates that brown dwarfs exhibit active atmospheric circulations. In this study we use a shallow-water model to investigate the global atmospheric dynamics in the stratified layer overlying the convective zone on these rapidly rotating objects. We show that the existence and properties of the atmospheric circulation crucially depend on key parameters including the energy injection rate and radiative timescale. Under conditions of strong internal heat flux and weak radiative dissipation, a banded flow pattern comprising east-west jet streams spontaneously emerges from the interaction of atmospheric turbulence with the planetary rotation. In contrast, when the internal heat flux is weak and/or radiative dissipation is strong, turbulence injected into the atmosphere damps before it can self-organize into jets, leading to a flow dominated by transient eddies and isotropic turbulence instead. The simulation results are not very sensitive to the form of the forcing. Based on the location of the transition between jet-dominated and eddy-dominated regimes, we suggest that many brown dwarfs may exhibit atmospheric circulations dominated by eddies and turbulence (rather than jets) due to the strong radiative damping on these worlds, but a jet structure is also possible under some realistic conditions. Our simulated light curves capture important features from observed infrared lightcurves of brown dwarfs, including amplitude variations of a few percent and shapes that fluctuate between single-peak and multi-peak structures. More broadly, our work shows that the shallow-water system provides a useful tool to illuminate fundamental aspects of the dynamics on these worlds.

Kepler-13Ab (= KOI-13.01) is a unique transiting hot Jupiter. It is one of very few known short-period planets orbiting a hot A-type star, making it one of the hottest planets currently known. The availability of Kepler data allows us to measure the planet's occultation (secondary eclipse) and phase curve in the optical, which we combine with occultations observed by warm Spitzer at 4.5 mic and 3.6 mic and a ground-based occultation observation in the Ks band (2.1 mic). We derive a day-side hemisphere temperature of 2,750 +- 160 K as the effective temperature of a black body showing the same occultation depths. Comparing the occultation depths with one-dimensional planetary atmosphere models suggests the presence of an atmospheric temperature inversion. Our analysis shows evidence for a relatively high geometric albedo, Ag= 0.33 +0.04 -0.06. While measured with a simplistic method, a high Ag is supported also by the fact that the one-dimensional atmosphere models underestimate the occultation depth in the optical. We use stellar spectra to determine the dilution, in the four wide bands where occultation was measured, due to the visual stellar binary companion 1.15 +- 0.05" away. The revised stellar parameters measured using these spectra are combined with other measurements leading to revised planetary mass and radius estimates of Mp = 4.94 - 8.09 Mjup and Rp = 1.406 +- 0.038 Rjup. Finally, we measure a Kepler mid-occultation time that is 34.0 +- 6.9 s earlier than expected based on the mid-transit time and the delay due to light travel time, and discuss possible scenarios.

We present three-dimensional atmospheric circulation models of GJ 1214b, a 2.7 Earth-radius, 6.5 Earth-mass super Earth detected by the MEarth survey. Here we explore the planet's circulation as a function of atmospheric metallicity and atmospheric composition, modeling atmospheres with a low mean-molecular weight (i.e., H2-dominated) and a high mean-molecular weight (i.e. water- and CO2-dominated). We find that atmospheres with a low mean-molecular weight have strong day-night temperature variations at pressures above the infrared photosphere that lead to equatorial superrotation. For these atmospheres, the enhancement of atmospheric opacities with increasing metallicity lead to shallower atmospheric heating, larger day-night temperature variations and hence stronger superrotation. In comparison, atmospheres with a high mean-molecular weight have larger day-night and equator-to-pole temperature variations than low mean-molecular weight atmospheres, but differences in opacity structure and energy budget lead to differences in jet structure. The circulation of a water-dominated atmosphere is dominated by equatorial superrotation, while the circulation of a CO2-dominated atmosphere is instead dominated by high-latitude jets. By comparing emergent flux spectra and lightcurves for 50x solar and water-dominated compositions, we show that observations in emission can break the degeneracy in determining the atmospheric composition of GJ 1214b. The variation in opacity with wavelength for the water-dominated atmosphere leads to large phase variations within water bands and small phase variations outside of water bands. The 50x solar atmosphere, however, yields small variations within water bands and large phase variations at other characteristic wavelengths. These observations would be much less sensitive to clouds, condensates, and hazes than transit observations.

As a planet ages it cools and its radius shrinks, at a rate set by the efficiency with which heat is transported from the interior out to space. The bottleneck for this transport is at the boundary between the convective interior and the radiative atmosphere; the opacity there sets the global cooling rate. Models of planetary evolution are often one-dimensional, such that the radiative-convective boundary (RCB) is defined by a single temperature, pressure, and opacity. In reality the spatially inhomogenous stellar heating pattern and circulation in the atmosphere could deform the RCB, allowing heat from the interior to escape more efficiently through regions with lower opacity. We present an analysis of the degree to which the RCB could be deformed and the resultant change in the evolutionary cooling rate. In this initial work we calculate the upper limit for this effect by comparing an atmospheric structure in local radiative equilibrium to its 1D equivalent. We find that the cooling through an uneven RCB could be enhanced over cooling through a uniform RCB by as much as 10-50%. We also show that the deformation of the RCB (and the enhancement of the cooling rate) increases with a greater incident stellar flux or a lower inner entropy. Our results indicate that this mechanism could significantly change a planet's thermal evolution, causing it to cool and shrink more quickly than would otherwise be expected. This may exacerbate the well known difficulty in explaining the very large radii observed for some hot Jupiters.

We present secondary eclipse observations of the highly irradiated transiting brown dwarf KELT-1b. These observations represent the first constraints on the atmospheric dynamics of a highly irradiated brown dwarf, and the atmospheres of irradiated giant planets at high surface gravity. Using the Spitzer Space Telescope, we measure secondary eclipse depths of 0.195+/-0.010% at 3.6um and 0.200+/-0.012% at 4.5um. We also find tentative evidence for the secondary eclipse in the z' band with a depth of 0.049+/-0.023%. These measured eclipse depths are most consistent with an atmosphere model in which there is a strong substellar hotspot, implying that heat redistribution in the atmosphere of KELT-1b is low. While models with a more mild hotspot or even with dayside heat redistribution are only marginally disfavored, models with complete heat redistribution are strongly ruled out. The eclipse depths also prefer an atmosphere with no TiO inversion layer, although a model with TiO inversion is permitted in the dayside heat redistribution case, and we consider the possibility of a day-night TiO cold trap in this object. For the first time, we compare the IRAC colors of brown dwarfs and hot Jupiters as a function of effective temperature. Importantly, our measurements reveal that KELT-1b has a [3.6]-[4.5] color of 0.07+/-0.11, identical to that of isolated brown dwarfs of similarly high temperature. In contrast, hot Jupiters generally show redder [3.6]-[4.5] colors of ~0.4, with a very large range from ~0 to ~1. Evidently, despite being more similar to hot Jupiters than to isolated brown dwarfs in terms of external forcing of the atmosphere by stellar insolation, KELT-1b has an atmosphere most like that of other brown dwarfs. This suggests that surface gravity is very important in controlling the atmospheric systems of substellar mass bodies.

We present the first three-dimensional magnetohydrodynamic (MHD) simulations of the atmosphere of HD 209458b which self-consistently include reduction of winds due to the Lorentz force and Ohmic heating. We find overall wind structures similar to that seen in previous models of hot Jupiter atmospheres, with strong equatorial jets and meridional flows poleward near the day side and equatorward near the night side. Inclusion of magnetic fields slows those winds and leads to Ohmic dissipation. We find wind slowing ranging from 10%-40% for reasonable field strengths. We find Ohmic dissipation rates ~10^17 W at 100 bar, orders of magnitude too small to explain the inflated radius of this planet. Faster wind speeds, not achievable in these anelastic calculations, may be able to increase this value somewhat, but likely will not be able to close the gap necessary to explain the inflated radius. We demonstrate that the discrepancy between the simulations presented here and previous models is due to inadequate treatment of magnetic field geometry and evolution. Induced poloidal fields become much larger than those imposed, highlighting the need for a self-consistent MHD treatment of these hot atmospheres.

Unlike hot Jupiters or other gas giants, super-Earths are expected to have a wide variety of compositions, ranging from terrestrial bodies like our own to more gaseous planets like Neptune. Observations of transiting systems, which allow us to directly measure planet masses and radii and constrain atmospheric properties, are key to understanding the compositional diversity of the planets in this mass range. Although Kepler has discovered hundreds of transiting super-Earth candidates over the past four years, the majority of these planets orbit stars that are too far away and too faint to allow for detailed atmospheric characterization and reliable mass estimates. Ground-based transit surveys focus on much brighter stars, but most lack the sensitivity to detect planets in this size range. One way to get around the difficulty of finding these smaller planets in transit is to start by choosing targets that are already known to contain super-Earth sized bodies detected using the radial velocity technique. Here we present results from a Spitzer program to observe six of the most favorable RV detected super-Earth systems, including HD 1461, HD 7924, HD 156668, HIP 57274, and GJ 876. We find no evidence for transits in any of their 4.5 micron flux light curves, and place limits on the allowed transit depths and corresponding planet radii that rule out even the most dense and iron-rich compositions for these objects. We also observed HD 97658, but the observation window was based on a possible ground-based transit detection (Henry et al. 2011) that was later ruled out; thus the window did not include the predicted time for the transit detection recently made by MOST (Dragomir et al. 2013).

We report secondary eclipse photometry of two hot Jupiters, WASP-48b and HAT-P-23b, at 3.6 and 4.5 um taken with the InfraRed Array Camera aboard the Spitzer Space Telescope during the warm Spitzer mission and in the H and Ks bands with the Wide Field IR Camera at the Palomar 200-inch Hale Telescope. WASP-48b and HAT-P-23b are Jupiter-mass and twice Jupiter-mass objects orbiting an old, slightly evolved F star and an early G dwarf star, respectively. In the H, Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.047% +/- 0.016%, 0.109% +/- 0.027%, 0.176% +/- 0.013%, and 0.214% +/- 0.020% for WASP-48b. In the Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.234% +/- 0.046%, 0.248% +/- 0.019%, and 0.309% +/- 0.026% for HAT-P-23b. For WASP-48b and HAT-P-23b, respectively, we measure delays of 2.6 +/- 3.9 minutes and 4.0 +/- 2.4 minutes relative to the predicted times of secondary eclipse for circular orbits, placing 2-sigma upper limits on |e cos(omega)| of 0.0053 and 0.0080, both of which are consistent with circular orbits. The dayside emission spectra of these planets are well-described by blackbodies with effective temperatures of 2158 +/- 100 K (WASP-48b) and 2154 +/- 90 K (HAT-P-23b), corresponding to moderate recirculation in the zero albedo case. Our measured eclipse depths are also consistent with one-dimensional radiative transfer models featuring varying degrees of recirculation and weak thermal inversions or no inversions at all. We discuss how the absence of strong temperature inversions on these planets may be related to the activity levels and metallicities of their host stars.

We present an optical to near-infrared transmission spectrum of the hot Jupiter HAT-P-1b, based on HST observations, covering the spectral regime from 0.29 to 1.027{\mu}m with STIS, which is coupled with a recent WFC3 transit (1.087 to 1.687{\mu}m). We derive refined physical parameters of the HAT-P-1 system, including an improved orbital ephemeris. The transmission spectrum shows a strong absorption signature shortward of 0.55{\mu}m, with a strong blueward slope into the near-ultraviolet. We detect atmospheric sodium absorption at a 3.3{\sigma} significance level, but find no evidence for the potassium feature. The red data implies a marginally flat spectrum with a tentative absorption enhancement at wavelength longer than ~0.85{\mu}m. The STIS and WFC3 spectra differ significantly in absolute radius level (4.3 +/- 1.6 pressure scale heights), implying strong optical absorption in the atmosphere of HAT-P-1b. The optical to near-infrared difference cannot be explained by stellar activity, as simulta- neous stellar activity monitoring of the G0V HAT-P-1b host star and its identical companion show no significant activity that could explain the result. We compare the complete STIS and WFC3 transmission spectrum with theoretical atmospheric mod- els which include haze, sodium and an extra optical absorber. We find that both an optical absorber and a super-solar sodium to water abundance ratio might be a sce- nario explaining the HAT-P-1b observations. Our results suggest that strong optical absorbers may be a dominant atmospheric feature in some hot Jupiter exoplanets.

Infrared lightcurves of transiting hot Jupiters present a trend in which the atmospheres of the hottest planets are less efficient at redistributing the stellar energy absorbed on their daysides---and thus have a larger day-night temperature contrast---than colder planets. No predictive atmospheric model has been published that identifies which dynamical mechanisms determine the atmospheric heat redistribution efficiency on tidally locked exoplanets. Here we present a two-layer shallow water model of the atmospheric dynamics on synchronously rotating planets that explains the observed trend. Our model shows that planets with weak friction and weak irradiation exhibit a banded zonal flow with minimal day-night temperature differences, while models with strong irradiation and/or strong friction exhibit a day-night flow pattern with order-unity fractional day-night temperature differences. To interpret the model, we develop a scaling theory that shows that the timescale for gravity waves to propagate horizontally over planetary scales, t_wave, plays a dominant role in controlling the transition from small to large temperature contrasts. This implies that heat redistribution is governed by a wave-like process, similar to the one responsible for the weak temperature gradients in the Earth's tropics. When atmospheric drag can be neglected, the transition from small to large day-night temperature contrasts occurs when t_wave ~ sqrt(t_rad/Omega), where t_rad is the radiative relaxation time and Omega is the planetary rotation frequency. Alternatively, this transition criterion can be expressed as t_rad ~ t_vert, where t_vert is the timescale for a fluid parcel to move vertically over the difference in day-night thickness. These results subsume the commonly used timescale comparison for estimating heat redistribution efficiency between t_rad and the global horizontal advection timescale, t_adv.

Hot Jupiters exhibit atmospheric temperatures ranging from hundreds to thousands of Kelvin. Because of their large day-night temperature differences, condensable species that are stable in the gas phase on the dayside, such as TiO and silicates, may condense and gravitationally settle on the nightside. Atmospheric circulation may counterbalance this tendency to gravitationally settle. This three dimensional (3D) mixing of chemical species has not previously been studied for hot Jupiters, yet it is crucial to assess the existence and distribution of TiO and silicates in the atmospheres of these planets. We perform 3D global circulation models of HD209458b including passive tracers that advect with the 3D flow, including a source/sink on the nightside to represent condensation and gravitational settling of haze particles. We show that global advection patterns produce strong vertical mixing that can keep condensable species lofted as long as they are trapped in particles of sizes of a few microns or less on the night side. We show that vertical mixing results not from small-scale convection but from the large-scale circulation driven by the day-night heating contrast. Although this vertical mixing is not diffusive in any rigorous sense, a comparison of our results with idealized diffusion models allows a rough estimate of the vertical diffusion coefficient. Kzz=5x10**4/Sqrt(Pbar) m2/s can be used in 1D models of HD 209458b. Moreover, our models exhibit strong spatial and temporal variability in the tracer concentration that could result in observable variations during transit or secondary eclipse measurements. Finally, we apply our model to the case of TiO in HD209458b and show that the day-night cold trap would deplete TiO if it condenses into particles bigger than a few microns on the planet's night side, making it unable to create the observed stratosphere of the planet.

We present HST optical transmission spectra of the transiting hot Jupiter WASP-12b, taken with the STIS instrument. From the transmission spectra, we are able to decisively rule out prominent absorption by TiO in the exoplanet's atmosphere. Strong pressure-broadened Na and K absorption signatures are also excluded, as are significant metal-hydride features. We compare our combined broadband spectrum to a wide variety of existing aerosol-free atmospheric models, though none are satisfactory fits. However, we do find that the full transmission spectrum can be described by models which include significant opacity from aerosols: including Rayleigh scattering, Mie scattering, tholin haze, and settling dust profiles. The transmission spectrum follows an effective extinction cross section with a power-law of index alpha, with the slope of the transmission spectrum constraining the quantity alphaT = -3528+/-660 K, where T is the atmospheric temperature. Rayleigh scattering (alpha=-4) is among the best fitting models, though requires low terminator temperatures near 900 K. Sub-micron size aerosol particles can provide equally good fits to the entire transmission spectrum for a wide range of temperatures, and we explore corundum as a plausible dust aerosol. The presence of atmospheric aerosols also helps to explain the modestly bright albedo implied by Spitzer observations, as well as the near black body nature of the emission spectrum. Ti-bearing condensates on the cooler night-side is the most natural explanation for the overall lack of TiO signatures in WASP-12b, indicating the day/night cold-trap is an important effect for very hot Jupiters. These finding indicate that aerosols can play a significant atmospheric role for the entire wide range of hot-Jupiter atmospheres, potentially affecting their overall spectrum and energy balance.(abridged)

We present Hubble Space Telescope near-infrared transmission spectroscopy of the transiting hot-Jupiter HAT-P-1b. We observed one transit with Wide Field Camera 3 using the G141 low-resolution grism to cover the wavelength range 1.087- 1.678 {\mu}m. These time series observations were taken with the newly available spatial scan mode that increases the duty cycle by nearly a factor of two, thus improving the resulting photometric precision of the data. We measure a planet-to-star radius ratio of Rp/R*=0.11709+/-0.00038 in the white light curve with the centre of transit occurring at 2456114.345+/-0.000133 (JD). We achieve S/N levels per exposure of 1840 (0.061%) at a resolution of {\Delta\lambda}=19.2nm (R~70) in the 1.1173 - 1.6549{\mu}m spectral region, providing the precision necessary to probe the transmission spectrum of the planet at close to the resolution limit of the instrument. We compute the transmission spectrum using both single target and differential photometry with similar results. The resultant transmission spectrum shows a significant absorption above the 5-{\sigma} level matching the 1.4{\mu}m water absorption band. In solar composition models, the water absorption is sensitive to the ~1 mbar pressure levels at the terminator. The detected absorption agrees with that predicted by an 1000 K isothermal model, as well as with that predicted by a planetary-averaged temperature model.

A variety of observations provide evidence for vigorous motion in the atmospheres of brown dwarfs and directly imaged giant planets. Motivated by these observations, we examine the dynamical regime of the circulation in the atmospheres and interiors of these objects. Brown dwarfs rotate rapidly, and for plausible wind speeds, the flow at large scales will be rotationally dominated. We present 3D, global, numerical simulations of convection in the interior, which demonstrate that, at large scales, the convection aligns in the direction parallel to the rotation axis. Convection occurs more efficiently at high latitudes than low latitudes, leading to systematic equator-to-pole temperature differences that may reach ~1 K near the top of the convection zone. The interaction of convection with the overlying, stably stratified atmosphere will generate a wealth of atmospheric waves, and we argue that, as in the stratospheres of planets in the solar system, the interaction of these waves with the mean flow will cause a significant atmospheric circulation at regional to global scales. At large scales, this should consist of stratified turbulence (possibly organizing into coherent structures such as vortices and jets) and an accompanying overturning circulation. We present an approximate analytic theory of this circulation, which predicts characteristic horizontal temperature variations of several to ~50 K, horizontal wind speeds of ~10-300 m/sec, and vertical velocities that advect air over a scale height in ~10^5-10^6 sec. This vertical mixing may help to explain the chemical disequilibrium observed on some brown dwarfs. Moreover, the implied large-scale organization of temperature perturbations and vertical velocities suggests that, near the L/T transition, patchy clouds can form near the photosphere, helping to explain recent observations of brown-dwarf variability in the near-IR.

Exoplanetary transmission spectroscopy in the near-infrared using Hubble/NICMOS is currently ambiguous because different observational groups claim different results from the same data, depending on their analysis methodologies. Spatial scanning with Hubble/WFC3 provides an opportunity to resolve this ambiguity. We here report WFC3 spectroscopy of the giant planets HD209458b and XO-1b in transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decorrelating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels collected in spatial scan mode. Our errors are within 6-percent (XO-1) and 26-percent (HD209458b) of the photon-limit at a spectral resolving power of 70, and are better than 0.01-percent per spectral channel. Both planets exhibit water absorption of approximately 200 ppm at the water peak near 1.38 microns. Our result for XO-1b contradicts the much larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for HD209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of an exoplanet atmosphere by Charbonneau et al. (2002). Model atmospheres having uniformly-distributed extra opacity of 0.012 cm^2 per gram account approximately for both our water measurement and the sodium absorption in this planet. Our results for HD209458b support the picture advocated by Pont et al. (2013) in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by continuous opacity due to haze and/or dust. However, the extra opacity needed for HD209458b is grayer than for HD189733b, with a weaker Rayleigh component.

We measure the transmission spectrum of WASP-19b from 3 transits using low-resolution optical spectroscopy from the HST Space Telescope Imaging Spectrograph (STIS). The STIS spectra cover a wavelength range of 0.29-1.03 microns with resolving power R=500. The optical data are combined with archival near-IR data from the HST Wide Field Camera 3 (WFC3) G141 grism, covering the wavelength range 1.087-1.687 micron, with resolving power R=130. We obtain the transmission spectrum from 0.53-1.687 microns with S/N levels between 3000 and 11,000 in 0.1 micron bins. WASP-19 is a very active star, with optical stellar flux varying by a few per cent over time. We correct the transit light curves for the effects of stellar activity using ground-based activity monitoring with the Cerro Tololo Inter-American Observatory (CTIO). While we were not able to construct a transmission spectrum using the blue optical data due to the presence of large occulted star spots, we were able to use the spot crossings to help constrain the mean stellar spot temperature. To search for predicted features in the hot-Jupiter atmosphere, we also define spectral indices for differential radius measurements to specifically search for the presence of TiO and alkali line features. Our measurements rule out TiO features predicted for a planet of WASP-19b's equilibrium temperature (2050 K) in the transmission spectrum at the 2.7-2.9 sigma confidence level, depending on atmospheric model formalism. The WFC3 transmission spectrum shows strong absorption features due to the presence of H2O, which is detected at the 4 sigma confidence level between 1.1 and 1.4 microns. The results indicate that WASP-19b is a planet with no or low levels of TiO and without a high C/O ratio. The lack of observable TiO features are possibly due to rainout, breakdown from stellar activity or the presence of other absorbers in the optical.

We present photometry of the extrasolar planet WASP-5b in the 3.6 and 4.5 micron bands taken with the Spitzer Space Telescope's Infrared Array Camera as part of the extended warm mission. By examining the depth of the planet's secondary eclipse at these two wavelengths, we can place joint constraints on the planet's atmospheric pressure-temperature profile and chemistry. We measure secondary eclipse depths of 0.197% +/- 0.028% and 0.227% +/- 0.025% in the 3.6 micron and 4.5 micron bands, respectively. Our observations are best matched by models showing a hot dayside and, depending on our choice of model, a weak thermal inversion or no inversion at all. We measure a mean offset from the predicted center of eclipse of 0.078 +/- 0.032 hours, translating to ecos(omega) = 0.0031 +/- 0.0013 and consistent with a circular orbit. We see no evidence for any eclipse timing variations comparable to those reported in a previous transit study.

The meteorology of hot Jupiters has been characterized primarily with thermal measurements, but recent observations suggest the possibility of directly detecting the winds by observing the Doppler shift of spectral lines seen during transit. Motivated by these observations, we show how Doppler measurements can place powerful constraints on the meteorology. We show that the atmospheric circulation--and Doppler signature--of hot Jupiters splits into two regimes. Under weak stellar insolation, the day-night thermal forcing generates fast zonal jet streams from the interaction of atmospheric waves with the mean flow. In this regime, air along the terminator (as seen during transit) flows toward Earth in some regions and away from Earth in others, leading to a Doppler signature exhibiting superposed blueshifted and redshifted components. Under intense stellar insolation, however, the strong thermal forcing damps these planetary-scale waves, inhibiting their ability to generate jets. Strong frictional drag likewise damps these waves and inhibits jet formation. As a result, this second regime exhibits a circulation dominated by high-altitude, day-to-night airflow, leading to a predominantly blueshifted Doppler signature during transit. We present state-of-the-art circulation models including nongray radiative transfer to quantify this regime shift and the resulting Doppler signatures; these models suggest that cool planets like GJ 436b lie in the first regime, HD 189733b is transitional, while planets hotter than HD 209458b lie in the second regime. Moreover, we show how the amplitude of the Doppler shifts constrains the strength of frictional drag in the upper atmospheres of hot Jupiters. If due to winds, the ~2-km/sec blueshift inferred on HD 209458b may require drag time constants as short as 10^4-10^6 seconds, possibly the result of Lorentz-force braking on this planet's hot dayside.

The investigation of planets around other stars began with the study of gas giants, but is now extending to the discovery and characterization of super-Earths and terrestrial planets. Motivated by this observational tide, we survey the basic dynamical principles governing the atmospheric circulation of terrestrial exoplanets, and discuss the interaction of their circulation with the hydrological cycle and global-scale climate feedbacks. Terrestrial exoplanets occupy a wide range of physical and dynamical conditions, only a small fraction of which have yet been explored in detail. Our approach is to lay out the fundamental dynamical principles governing the atmospheric circulation on terrestrial planets--broadly defined--and show how they can provide a foundation for understanding the atmospheric behavior of these worlds. We first survey basic atmospheric dynamics, including the role of geostrophy, baroclinic instabilities, and jets in the strongly rotating regime (the "extratropics") and the role of the Hadley circulation, wave adjustment of the thermal structure, and the tendency toward equatorial superrotation in the slowly rotating regime (the "tropics"). We then survey key elements of the hydrological cycle, including the factors that control precipitation, humidity, and cloudiness. Next, we summarize key mechanisms by which the circulation affects the global-mean climate, and hence planetary habitability. In particular, we discuss the runaway greenhouse, transitions to snowball states, atmospheric collapse, and the links between atmospheric circulation and CO2 weathering rates. We finish by summarizing the key questions and challenges for this emerging field in the future.

The ongoing characterization of hot Jupiters has motivated a variety of circulation models of their atmospheres. Such models must be integrated starting from an assumed initial state, which is typically taken to be a wind-free, rest state. Here, we investigate the sensitivity of hot-Jupiter atmospheric circulation models to initial conditions. We consider two classes of models--shallow-water models, which have proven successful at illuminating the dynamical mechanisms at play on these planets, and full three-dimensional models similar to those being explored in the literature. Models are initialized with zonal jets, and we explore a variety of different initial jet profiles. We demonstrate that, in both classes of models, the final, equilibrated state is independent of initial condition--as long as frictional drag near the bottom of the domain and/or interaction with a specified planetary interior are included so that the atmosphere can adjust angular momentum over time relative to the interior. When such mechanisms are included, otherwise identical models initialized with vastly different initial conditions all converge to the same statistical steady state. In some cases, the models exhibit modest time variability; this variability results in random fluctuations about the statistical steady state, but we emphasize that, even in these cases, the statistical steady state itself does not depend on initial conditions. Although the outcome of hot-Jupiter circulation models depend on details of the radiative forcing and frictional drag, aspects of which remain uncertain, we conclude that the specification of initial conditions is not a source of uncertainty, at least over the parameter range explored in most current models.

We report near-infrared measurements of the terminator region transmission spectrum and dayside emission spectrum of the exoplanet WASP-12b obtained using the HST WFC3 instrument. The disk-average dayside brightness temperature averages about 2900 K, peaking to 3200 K around 1.46 microns. We modeled a range of atmospheric cases for both the emission and transmission spectrum and confirm the recent finding by Crossfield et al. (2012b) that there is no evidence for C/O >1 in the atmosphere of WASP-12b. Assuming a physically plausible atmosphere, we find evidence that the presence of a number of molecules is consistent with the data, but the justification for inclusion of these opacity sources based on the Bayesian Information Criterion (BIC) is marginal. We also find the near-infrared primary eclipse light curve is consistent with small amounts of prolate distortion. As part of the calibration effort for these data, we conducted a detailed study of instrument systematics using 65 orbits of WFC3-IR grims observations. The instrument systematics are dominated by detector-related affects, which vary significantly depending on the detector readout mode. The 256x256 subarray observations of WASP 12 produced spectral measurements within 15% of the photon-noise limit using a simple calibration approach. Residual systematics are estimated to be less than 70 parts per million.

We present Warm Spitzer/IRAC secondary eclipse time series photometry of three short-period transiting exoplanets, HAT-P-3b, HAT-P-4b and HAT-P-12b, in both the available 3.6 and 4.5 micron bands. HAT-P-3b and HAT-P-4b are Jupiter-mass, objects orbiting an early K and an early G dwarf stars, respectively. For HAT-P-3b we find eclipse depths of 0.112%+0.015%-0.030% (3.6 micron) and 0.094%+0.016%-0.009% (4.5 micron). The HAT-P-4b values are 0.142%+0.014%-0.016% (3.6 micron) and 0.122%+0.012%-0.014% (4.5micron). The two planets' photometry is consistent with inefficient heat redistribution from their day to night sides (and low albedos), but it is inconclusive about possible temperature inversions in their atmospheres. HAT-P-12b is a Saturn-mass planet and is one of the coolest planets ever observed during secondary eclipse, along with hot Neptune GJ 436b and hot Saturn WASP-29b. We are able to place 3$\sigma$ upper limits on the secondary eclipse depth of HAT-P-12b in both wavelengths: < 0.042% (3.6 micron) and <0.085% (4.5 micron). We discuss these results in the context of the {\it Spitzer} secondary eclipse measurements of GJ 436b and WASP-29b. It is possible that we do not detect the eclipses of HAT-P-12b due to high eccentricity, but find that weak planetary emission in these wavelengths is a more likely explanation. We place 3 sigma upper limits on the |e cos(omega)| quantity (where e is eccentricity and omega is the argument of periapsis) for HAT-P-3b (<0.0081) and HAT-P-4b (<0.0042), based on the secondary eclipse timings.

Of the over 800 exoplanets detected to date, over half are on non-circular orbits, with eccentricities as high as 0.93. Such orbits lead to time-variable stellar heating, which has implications for the planet's atmospheric dynamical regime. However, little is known about this dynamical regime, and how it may influence observations. Therefore, we present a systematic study of hot Jupiters on highly eccentric orbits using the SPARC/MITgcm, a model which couples a three-dimensional general circulation model with a plane-parallel, two-stream, non-grey radiative transfer model. In our study, we vary the eccentricity and orbit-average stellar flux over a wide range. We demonstrate that the eccentric hot Jupiter regime is qualitatively similar to that of planets on circular orbits; the planets possess a superrotating equatorial jet and exhibit large day-night temperature variations. We show that these day-night heating variations induce momentum fluxes equatorward to maintain the superrotating jet throughout its orbit. As the eccentricity and/or stellar flux is increased, the superrotating jet strengthens and narrows, due to a smaller Rossby deformation radius. For a select number of model integrations, we generate full-orbit lightcurves and find that the timing of transit and secondary eclipse viewed from Earth with respect to periapse and apoapse can greatly affect what we see in infrared (IR) lightcurves; the peak in IR flux can lead or lag secondary eclipse depending on the geometry. For those planets that have large day-night temperature variations and rapid rotation rates, we find that the lightcurves exhibit "ringing" as the planet's hottest region rotates in and out of view from Earth. These results can be used to explain future observations of eccentric transiting exoplanets.

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 \mu m bands of the Spitzer Space Telescope. The 3.6 and 4.5 \mu m data sets span an entire orbital period of HAT-P-2b, making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 \mu m that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39+/-0.28, 5.84+/-0.39, and 4.68+/-0.37 hours after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138%+/-0.0089%, 0.1162%+/-0.0080%, and 0.1888%+/-0.0072% in the 3.6, 4.5, and 8.0 \mu m bands respectively. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity and argument of periapse of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism.

We present an updated analysis of Jupiter's equatorial meteorology from Cassini observations. For two months preceding the spacecraft's closest approach, the Imaging Science Subsystem (ISS) onboard regularly imaged the atmosphere. We created time-lapse movies from this period in order to analyze the dynamics of equatorial hot spots and their interactions with adjacent latitudes. Hot spots are quasi-stable, rectangular dark areas on visible-wavelength images, with defined eastern edges that sharply contrast with surrounding clouds, but diffuse western edges serving as nebulous boundaries with adjacent equatorial plumes. Hot spots exhibit significant variations in size and shape over timescales of days and weeks. Some of these changes correspond with passing vortex systems from adjacent latitudes interacting with hot spots. Strong anticyclonic gyres present to the south and southeast of the dark areas appear to circulate into hot spots. Impressive, bright white plumes occupy spaces in between hot spots. Compact cirrus-like 'scooter' clouds flow rapidly through the plumes before disappearing within the dark areas. These clouds travel at 150-200 m/s, much faster than the 100 m/s hot spot and plume drift speed. This raises the possibility that the scooter clouds may be more illustrative of the actual jet stream speed at these latitudes. Most previously published zonal wind profiles represent the drift speed of the hot spots at their latitude from pattern matching of the entire longitudinal image strip. If a downward branch of an equatorially-trapped Rossby waves controls the overall appearance of hot spots, however, the westward phase velocity of the wave leads to underestimates of the true jet stream speed.

We present results regarding the dynamical meteorology of Jupiter's White Ovals at different points in their evolution. Starting from the era with three White Ovals FA, BC, and DE (Galileo), continuing to the post-merger epoch with only one Oval BA (Cassini), and finally to Oval BA's current reddened state (New Horizons), we demonstrate that the dynamics of their flow have similarly evolved along with their appearance. The northern periphery of Oval BA increased in speed by 20 m/s from Cassini to New Horizons, ending up at a speed nearly identical to that of the northern periphery of Oval DE during Galileo. However, the peak speeds along the southern rim of the newly formed Oval BA were consistently faster than the corresponding speeds in Oval DE, and they increased still further between Cassini and New Horizons, ending up at ~140 to 150 m/s. The modest strengthening of the winds in Oval BA, the appearance of red aerosols, and the appearance of a turbulent, cyclonic feature to Oval BA's northwest create a strong resemblance with the Great Red Spot from both a dynamical and morphological perspective. In addition to the White Ovals, we also measure the winds within two compact cyclonic regions, one in the Galileo data set and one in the Cassini data set. In the images, these cyclonic features appear turbulent and filamentary, but our wind field reveals that the flow manifests as a coherent high-speed collar surrounding relatively quiescent interiors. Although less studied than their anticyclonic cousins, these cyclones may offer crucial insights into the planet's cloud-level energetics and dynamical meteorology.

We have produced mosaics of the Great Red Spot (GRS) using images taken by the Galileo spacecraft in May 2000, and have measured the winds of the GRS using an automated algorithm that does not require manual cloud tracking. Our technique yields a high-density, regular grid of wind velocity vectors that is advantageous over a limited number of scattered wind vectors that result from manual cloud tracking. The high-velocity collar of the GRS is clearly seen from our velocity vector map, and highest wind velocities are measured to be around 170 m/s. The high resolution of the mosaics have also enabled us to map turbulent eddies inside the chaotic central region of the GRS, similar to those mapped by Sada et al. (1996) and Vasavada et al. (1998). Using the wind velocity measurements, we computed particle trajectories around the GRS as well as maps of relative and absolute vorticities. We have discovered a narrow ring of cyclonic vorticity that surrounds the main anti-cyclonic high-velocity collar. This narrow ring appears to correspond to a ring surrounding the GRS that is bright in 5-microns (Terrile and Beebe 1979). It appears that this cyclonic ring is not a transient feature of the GRS, as we have discovered it in a re-analysis of Galileo data taken in 1996 first analyzed by Vasavada et al. (1998). We also calculate how absolute vorticity changes as a function of latitude along a trajectory around the GRS and compare these measurements to similar ones performed by Dowling and Ingersoll (1988) using Voyager data. We show no dramatic evolution in the structure of the GRS since the Voyager era except for additional evidence for a counterrotating GRS core, an increase in velocity in the main velocity collar, and an overall decrease in the length of the GRS.

We present suggestive evidence for an inverse energy cascade within Jupiter's atmosphere through a calculation of the power spectrum of its kinetic energy and its cloud patterns. Using Cassini observations, we composed full-longitudinal mosaics of Jupiter's atmosphere at several wavelengths. We also utilized image pairs derived from these observations to generate full-longitudinal maps of wind vectors and atmospheric kinetic energy within Jupiter's troposphere. We computed power spectra of the image mosaics and kinetic energy maps using spherical harmonic analysis. Power spectra of Jupiter's cloud patterns imaged at certain wavelengths resemble theoretical spectra of two-dimensional turbulence, with power-law slopes near -5/3 and -3 at low and high wavenumbers, respectively. The slopes of the kinetic energy power spectrum are also near -5/3 at low wavenumbers. At high wavenumbers, however, the spectral slopes are relatively flatter than the theoretical prediction of -3. Our results also show the importance of calculating spectral slopes from full 2D velocity maps rather than 1D zonal mean velocity profiles, since at large wavenumbers the spectra differ significantly, though at low wavenumbers, the 1D zonal and full 2D kinetic energy spectra are practically indistinguishable. Furthermore, the difference between the image and kinetic energy spectra suggests some caution in the interpretation of power spectrum results solely from image mosaics and its significance for the underlying dynamics. Finally, we also report prominent variations in kinetic energy within the equatorial jet stream that appear to be associated with the 5 {\mu}m hotspots. Other eddies are present within the flow collar of the Great Red Spot, suggesting caution when interpreting snapshots of the flow inside these features as representative of a time-averaged state.

We present an analysis of data about Saturn's atmosphere from Cassini's Visual and Infrared Mapping Spectrometer (VIMS), focusing on the meteorology of the features seen in the 5-micron spectral window. We present VIMS mosaics and discuss the morphology and general characteristics of the features backlit by Saturn's thermal emission. We have also constructed a zonal wind profile from VIMS feature tracking observation sequences using an automated cloud feature tracker. Comparison with previously constructed profiles from Voyager and Cassini imaging data reveals broad similarities, suggesting minimal vertical shear of the zonal wind. However, areas of apparent wind shear are present in the VIMS zonal wind profile at jet stream cores. In particular, our analysis shows that the equatorial jet reaches speeds exceeding 450 m/s, similar to speeds obtained during the Voyager era. This suggests that recent inferences of relatively slower jet speeds of ~275-375 m/s are confined to the upper troposphere and that the deep (>1 bar) jet has not experienced a significant slowdown. Our measurements of the numerous dark, spotted features seen in the VIMS mosaics reveals that most of these features have diameters less than 1000 km and reside in confined zonal bands between jet stream cores. We propose that these spot features are vortices and that VIMS and ISS are sensing the same vortices at two different pressure levels. The local structure at the zonal jet streams remains complex, as VIMS may be sensing cloud features that are deeper than the NH3 cloud deck.

Heterogeneous clouds or temperature perturbations in rotating brown dwarfs produce variability in the observed flux. We report time-resolved simultaneous observations of the variable T6.5 brown dwarf 2MASSJ22282889-431026 over the wavelength ranges 1.1-1.7 microns and broadband 4.5 microns. Spectroscopic observations were taken with Wide Field Camera 3 on board the Hubble Space Telescope and photometry with the Spitzer Space Telescope. The object shows sinusoidal infrared variability with a period of 1.4 hr at most wavelengths with peak-to-peak amplitudes between 1.45% and 5.3% of the mean flux. While the light curve shapes are similar at all wavelengths, their phases differ from wavelength to wavelength with a maximum difference of more than half of a rotational period. We compare the spectra with atmospheric models of different cloud prescriptions, from which we determine the pressure levels probed at different wavelengths. We find that the phase lag increases with decreasing pressure level, or higher altitude. We discuss a number of plausible scenarios that could cause this trend of light curve phase with probed pressure level. These observations are the first to probe heterogeneity in an ultracool atmosphere in both horizontal and vertical directions, and thus are an ideal test case for realistic three dimensional simulations of the atmospheric structure with clouds in brown dwarfs and extrasolar planets.

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 micron, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733b's 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]

We report the results of an analysis of all Spitzer/MIPS 24 micron observations of HD 209458b, one of the touchstone objects in the study of irradiated giant planet atmospheres. Altogether we analyze 2.5 transits, 3 eclipses, and a 58-hour near-continuous observation designed to detect the planet's thermal phase curve. The results of our analysis are: (1) A mean transit depth of 1.484% +/- 0.035%, consistent with previous measurements and showing no evidence of variability in transit depth at the 3% level. (2) A mean eclipse depth of 0.332% +/- 0.026%, somewhat higher than that previously reported for this system; this new value brings observations into better agreement with models. The dayside flux shows no evidence of variability at the 12% level. (3) Eclipses in the system occur 32 s +/- 129 s earlier than would be expected from a circular orbit, which constrains the orbital quantity (e cos omega) to be 0.00004 +/- 0.00033. This result is fully consistent with a circular orbit and sets an upper limit of 140 m/s (3 sigma) on any eccentricity-induced velocity offset during transit. The phase curve observations (including one of the transits) exhibit an anomalous trend similar to the detector ramp seen in previous Spitzer/IRAC observations; by modeling this ramp we recover the system parameters. The photometry which follows the ramp and transit exhibits a gradual, ~0.2% decrease in flux, similar to that seen in pre-launch calibration data. The large uncertainties associated with this poorly-understood, likely instrumental effect prevent us from usefully constraining the planet's thermal phase curve. Our observations highlight the need for a thorough understanding of detector-related instrumental effects on long time scales when making the high-precision mid-infrared measurements planned for future missions such as EChO, SPICA, and JWST. [abridged]

We present a new model for Ellipsoidal Variations Induced by a Low-Mass Companion, the EVIL-MC model. We employ several approximations appropriate for planetary systems to substantially increase the computational efficiency of our model relative to more general ellipsoidal variation models and improve upon the accuracy of simpler models. This new approach gives us a unique ability to rapidly and accurately determine planetary system parameters. We use the EVIL-MC model to analyze Kepler Quarter 0-2 (Q0-2) observations of the HAT-P-7 system, an F-type star orbited by a nearly Jupiter-mass companion. Our analysis corroborates previous estimates of the planet-star mass ratio q = (1.10 +/- 0.06) x 10^(-3), and we have revised the planet's dayside brightness temperature to 2680 +10/-20 K. We also find a large difference between the day- and nightside planetary flux, with little nightside emission. Preliminary dynamical+radiative modeling of the atmosphere indicates this result is qualitatively consistent with high altitude absorption of stellar heating. Similar analyses of Kepler and CoRoT photometry of other planets using EVIL-MC will play a key role in providing constraints on the properties of many extrasolar systems, especially given the limited resources for follow-up and characterization of these systems. However, as we highlight, there are important degeneracies between the contributions from ellipsoidal variations and planetary emission and reflection. Consequently, for many of the hottest and brightest Kepler and CoRoT planets, accurate estimates of the planetary emission and reflection, diagnostic of atmospheric heat budgets, will require accurate modeling of the photometric contribution from the stellar ellipsoidal variation.

We have analyzed Warm Spitzer/IRAC observations of the secondary eclipses of three planets, XO-4b, HAT-P-6b and HAT-P-8b. We measure secondary eclipse amplitudes at 3.6{\mu}m and 4.5{\mu}m for each target. XO-4b exhibits a stronger eclipse depth at 4.5{\mu}m than at 3.6{\mu}m, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6{\mu}m, and is best-described by models without a temperature inversion. The eclipse depths of HAT-P-6b can be fitted with models with a small or no temperature inversion. We consider our results in the context of a postulated relationship between stellar activity and temperature inversions and a relationship between irradiation level and planet dayside temperature, as discussed by Knutson et al. (2010) and Cowan & Agol (2011), respectively. Our results are consistent with these hypotheses, but do not significantly strengthen them. To measure accurate secondary eclipse central phases, we require accurate ephemerides. We obtain primary transit observations and supplement them with publicly available observations to update the orbital ephemerides of the three planets. Based on the secondary eclipse timing, we set upper boundaries for e cos(\omega) for HAT-P-6b, HAT-P-8b and XO-4b and find that the values are consistent with circular orbits.

The increasing richness of exoplanet observations has motivated a variety of three-dimensional atmospheric circulation models of these planets. Under strongly irradiated conditions, models of tidally locked, short-period planets (both hot Jupiters and terrestrial planets) tend to exhibit a circulation dominated by a fast eastward, or "superrotating," jet stream at the equator. Under appropriate conditions, this phenomenon can cause the hottest regions to be displaced eastward from the substellar point by tens of degrees longitude. Such an offset has been observed on HD 189733b, supporting the possibility of equatorial jets on short-period exoplanets. Despite its relevance, however, the dynamical mechanisms responsible for generating the equatorial superrotation in such models have not been identified. Here, we show that the equatorial jet results from the interaction of the mean flow with standing, planetary-scale Rossby and Kelvin waves induced by the day-night thermal forcing. The Rossby waves develop phase tilts that pump eastward momentum from high latitudes to the equator, thereby inducing equatorial superrotation. We present an analytic theory demonstrating this mechanism and explore its properties in a hierarchy of one-layer (shallow-water) calculations and fully 3D models. The wave-mean-flow interaction produces an equatorial jet whose latitudinal width is comparable to that of the Rossby waves, namely the equatorial Rossby deformation radius modified by radiative and frictional effects. For conditions typical of synchronously rotating hot Jupiters, this length is comparable to a planetary radius, explaining the broad scale of the equatorial jet obtained in most hot Jupiter models. Our theory illuminates the dependence of the equatorial jet speed on forcing amplitude, strength of friction, and other parameters, as well as the conditions under which jets can form at all.

We have developed 1-D photochemical and thermochemical kinetics and diffusion models for the transiting exoplanets HD 189733b and HD 209458b to study the effects of disequilibrium chemistry on the atmospheric composition of "hot Jupiters." Here we investigate the coupled chemistry of neutral carbon, hydrogen, oxygen, and nitrogen species, and we compare the model results with existing transit and eclipse observations. We find that the vertical profiles of molecular constituents are significantly affected by transport-induced quenching and photochemistry, particularly on cooler HD 189733b; however, the warmer stratospheric temperatures on HD 209458b can help maintain thermochemical equilibrium and reduce the effects of disequilibrium chemistry. For both planets, the methane and ammonia mole fractions are found to be enhanced over their equilibrium values at pressures of a few bar to less than a mbar due to transport-induced quenching, but CH4 and NH3 are photochemically removed at higher altitudes. Atomic species, unsaturated hydrocarbons (particularly C2H2), some nitriles (particularly HCN), and radicals like OH, CH3, and NH2 are enhanced overequilibrium predictions because of quenching and photochemistry. In contrast, CO, H2O, N2, and CO2 more closely follow their equilibrium profiles, except at pressures < 1 microbar, where CO, H2O, and N2 are photochemically destroyed and CO2 is produced before its eventual high-altitude destruction. The enhanced abundances of HCN, CH4, and NH3 in particular are expected to affect the spectral signatures and thermal profiles HD 189733b and other, relatively cool, close-in transiting exoplanets. We examine the sensitivity of our results to the assumed temperature structure and eddy diffusion coefficientss and discuss further observational consequences of these models.

Theoretical arguments and observations suggest that the atmospheres of Brown Dwarfs and planets are very dynamic on chemical and on physical time scales. The modelling of such substellar atmospheres has, hence, been much more demanding than initially anticipated. This Splinter (http://star-www.st-and.ac.uk/~ch80/CS16/MultiDSplinter_CS16.html) has combined new developments in atmosphere modelling, with novel observational techniques, and new challenges arising from planetary and space weather observations.

We present photometry of the giant extrasolar planet WASP-4b at 3.6 and 4.5 micron taken with the Infrared Array Camera on board the Spitzer Space Telescope as part of Spitzer's extended warm mission. We find secondary eclipse depths of 0.319+/-0.031% and 0.343+/-0.027% for the 3.6 and 4.5 micron bands, respectively and show model emission spectra and pressure-temperature profiles for the planetary atmosphere. These eclipse depths are well fit by model emission spectra with water and other molecules in absorption, similar to those used for TrES-3 and HD 189733b. Depending on our choice of model, these results indicate that this planet has either a weak dayside temperature inversion or no inversion at all. The absence of a strong thermal inversion on this highly irradiated planet is contrary to the idea that highly irradiated planets are expected to have inversions, perhaps due the presence of an unknown absorber in the upper atmosphere. This result might be explained by the modestly enhanced activity level of WASP-4b's G7V host star, which could increase the amount of UV flux received by the planet, therefore reducing the abundance of the unknown stratospheric absorber in the planetary atmosphere as suggested in Knutson et al. (2010). We also find no evidence for an offset in the timing of the secondary eclipse and place a 2 sigma upper limit on |ecos(omega)| of 0.0024, which constrains the range of tidal heating models that could explain this planet's inflated radius.

Three-dimensional studies of convection in deep spherical shells have been used to test the hypothesis that the strong jet streams on Jupiter, Saturn, Uranus, and Neptune result from convection throughout the molecular envelopes. Due to computational limitations, these simulations must adopt viscosities and heat fluxes many orders of magnitude larger than the planetary values. Several numerical investigations have identified trends for how the mean jet speed varies with heat flux and viscosity, but no previous theories have been advanced to explain these trends. Here, we show using simple arguments that if convective release of potential energy pumps the jets and viscosity damps them, the mean jet speeds split into two regimes. When the convection is weakly nonlinear, the equilibrated jet speeds should scale approximately with F/nu, where F is the convective heat flux and nu is the viscosity. When the convection is strongly nonlinear, the jet speeds are faster and should scale approximately as (F/nu)^{1/2}. We demonstrate how this regime shift can naturally result from a shift in the behavior of the jet-pumping efficiency with heat flux and viscosity. Moreover, the simulations hint at a third regime where, at sufficiently small viscosities, the jet speed becomes independent of the viscosity. We show based on mixing-length estimates that if such a regime exists, mean jet speeds should scale as heat flux to the 1/4 power. Our scalings provide a good match to the mean jet speeds obtained in previous Boussinesq and anelastic, three-dimensional simulations of convection within giant planets over a broad range of parameters. When extrapolated to the real heat fluxes, these scalings suggest that the mass-weighted jet speeds in the molecular envelopes of the giant planets are much weaker--by an order of magnitude or more--than the speeds measured at cloud level.

We measure secondary eclipses of the hot giant exoplanets CoRoT-1 at 3.6 and 4.5 microns, and CoRoT-2 at 3.6 microns, both using Warm Spitzer. We find that the Warm Spitzer mission is working very well for exoplanet science. For consistency of our analysis we also re-analyze archival cryogenic Spitzer data for secondary eclipses of CoRoT-2 at 4.5 and 8 microns. We compare the total data for both planets, including optical eclipse measurements by the CoRoT mission, and ground-based eclipse measurements at 2 microns, to existing models. Both planets exhibit stronger eclipses at 4.5 than at 3.6 microns, which is often indicative of an atmospheric temperature inversion. The spectrum of CoRoT-1 is best reproduced by a 2460K blackbody, due either to a high altitude layer that strongly absorbs stellar irradiance, or an isothermal region in the planetary atmosphere. The spectrum of CoRoT-2 is unusual because the 8 micron contrast is anomalously low. Non-inverted atmospheres could potentially produce the CoRoT-2 spectrum if the planet exhibits line emission from CO at 4.5 microns, caused by tidal-induced mass loss. However, the viability of that hypothesis is questionable because the emitting region cannot be more than about 30-percent larger than the planet's transit radius, based on the ingress and egress times at eclipse. An alternative possibility to account for the spectrum of CoRoT-2 is an additional opacity source that acts strongly at wavelengths less than 5 microns, heating the upper atmosphere while allowing the deeper atmosphere seen at 8 microns to remain cooler. We obtain a similar result as Gillon et al. for the phase of the secondary eclipse of CoRoT-2, implying an eccentric orbit with e*cos(omega)=-0.0030 +/- 0.0004.

GJ436b is a unique member of the transiting extrasolar planet population being one of the smallest and least irradiated and possessing an eccentric orbit. Because of its size, mass and density, GJ436b could plausibly have an atmospheric metallicity similar to Neptune (20-60 times solar abundances), which makes it an ideal target to study the effects of atmospheric metallicity on dynamics and radiative transfer in an extrasolar planetary atmosphere. We present three-dimensional atmospheric circulation models that include realistic non-gray radiative transfer for 1, 3, 10, 30, and 50 times solar atmospheric metallicity cases of GJ436b. Low metallicity models (1 and 3 times solar) show little day/night temperature variation and strong high-latitude jets. In contrast, higher metallicity models (30 and 50 times solar) exhibit day/night temperature variations and a strong equatorial jet. Spectra and light curves produced from these simulations show strong orbital phase dependencies in the 50 times solar case and negligible variations with orbital phase in the 1 times solar case. Comparisons between the predicted planet/star flux ratio from these models and current secondary eclipse measurements support a high metallicity atmosphere (30-50 times solar abundances) with disequilibrium carbon chemistry at play for GJ436b. Regardless of the actual atmospheric composition of GJ436b, our models serve to illuminate how metallicity influences the atmospheric circulation for a broad range of warm extrasolar planets.

We compute models of the transmission spectra of planets HD 209458b, HD 189733b, and generic hot Jupiters. We examine the effects of temperature, surface gravity, and metallicity for the generic planets as a guide to understanding transmission spectra in general. We find that carbon dioxide absorption at 4.4 and 15 microns is prominent at high metallicity, and is a clear metallicity indicator. For HD 209458b and HD 189733b, we compute spectra for both one-dimensional and three-dimensional model atmospheres and examine the differences between them. The differences are usually small, but can be large if atmospheric temperatures are near important chemical abundance boundaries. The calculations for the 3D atmospheres, and their comparison with data, serve as constraints on these dynamical models that complement the secondary eclipse and light curve data sets. For HD 209458b, even if TiO and VO gases are abundant on the day side, their abundances can be considerably reduced on the cooler planetary limb. However, given the predicted limb temperatures and TiO abundances, the model's optical opacity is too high. For HD 189733b we find a good match with some infrared data sets and constrain the altitude of a postulated haze layer. For this planet, substantial differences can exist between the transmission spectra of the leading and trailing hemispheres, which is an excellent probe of carbon chemistry. In thermochemical equilibrium, the cooler leading hemisphere is methane-dominated, and the hotter trailing hemisphere is CO-dominated, but these differences may be eliminated by non-equilibrium chemistry due to vertical mixing. It may be possible to constrain the carbon chemistry of this planet, and its spatial variation, with JWST.

We survey the basic principles of atmospheric dynamics relevant to explaining existing and future observations of exoplanets, both gas giant and terrestrial. Given the paucity of data on exoplanet atmospheres, our approach is to emphasize fundamental principles and insights gained from Solar-System studies that are likely to be generalizable to exoplanets. We begin by presenting the hierarchy of basic equations used in atmospheric dynamics, including the Navier-Stokes, primitive, shallow-water, and two-dimensional nondivergent models. We then survey key concepts in atmospheric dynamics, including the importance of planetary rotation, the concept of balance, and scaling arguments to show how turbulent interactions generally produce large-scale east-west banding on rotating planets. We next turn to issues specific to giant planets, including their expected interior and atmospheric thermal structures, the implications for their wind patterns, and mechanisms to pump their east-west jets. Hot Jupiter atmospheric dynamics are given particular attention, as these close-in planets have been the subject of most of the concrete developments in the study of exoplanetary atmospheres. We then turn to the basic elements of circulation on terrestrial planets as inferred from Solar-System studies, including Hadley cells, jet streams, processes that govern the large-scale horizontal temperature contrasts, and climate, and we discuss how these insights may apply to terrestrial exoplanets. Although exoplanets surely possess a greater diversity of circulation regimes than seen on the planets in our Solar System, our guiding philosophy is that the multi-decade study of Solar-System planets reviewed here provides a foundation upon which our understanding of more exotic exoplanetary meteorology must build.

Three-dimensional numerical simulations show that large-scale latent heating resulting from condensation of water vapor can produce multiple zonal jets similar to those on the gas giants (Jupiter and Saturn) and ice giants (Uranus and Neptune). For plausible water abundances (3-5 times solar on Jupiter/Saturn and 30 times solar on Uranus/Neptune), our simulations produce ~20 zonal jets for Jupiter and Saturn and 3 zonal jets on Uranus and Neptune, similar to the number of jets observed on these planets. Moreover, these Jupiter/Saturn cases produce equatorial superrotation whereas the Uranus/Neptune cases produce equatorial subrotation,consistent with the observed equatorial jet direction on these planets.Sensitivity tests show that water abundance, planetary rotation rate,and planetary radius are all controlling factors, with water playing the most important role; modest water abundances, large planetary radii, and fast rotation rates favor equatorial superrotation, whereas large water abundances favor equatorial subrotation regardless of the planetary radius and rotation rate. Given the larger radii, faster rotation rates, and probable lower water abundances of Jupiter and Saturn relative to Uranus and Neptune, our simulations therefore provide a possible mechanism for the existence of equatorial superrotation on Jupiter and Saturn and the lack of superrotation on Uranus and Neptune. Nevertheless, Saturn poses a possible difficulty, as our simulations were unable to explain the unusually high speed (~400m/s) of that planet's superrotating jet. The zonal jets in our simulations exhibit modest violations of the barotropic and Charney-Stern stability criteria. Overall, our simulations, while idealized, support the idea that latent heating plays an important role in generating the jets on the giant planets.

Several processes can cause the shape of an extrasolar giant planet's shadow, as viewed in transit, to depart from circular. In addition to rotational effects, cloud formation, non-homogenous haze production and movement, and dynamical effects (winds) could also be important. When such a planet transits its host star as seen from Earth, the asphericity will introduce a deviation in the transit lightcurve relative to the transit of a perfectly spherical (or perfectly oblate) planet. We develop a theoretical framework to interpret planet shapes. We then generate predictions for transiting planet shapes based on a published theoretical dynamical model of HD189733b. Using these shape models we show that planet shapes are unlikely to introduce detectable lightcurve deviations (those >~1e-5 of the host star), but that the shapes may lead to astrophysical sources of systematic error when measuring planetary oblateness, transit time, and impact parameter.

Oxygen (O) and carbon (C) have been inferred recently to be subsolar in abundance from spectra of the atmosphere of the transiting hot Jupiter HD 189733b. Yet, the mass and radius of the planet coupled with structure models indicate a strongly supersolar abundance of heavy elements in the interior of this object. Here we explore the discrepancy between the large amount of heavy elements suspected in the planet's interior and the paucity of volatiles measured in its atmosphere. We describe the formation sequence of the icy planetesimals formed beyond the snow line of the protoplanetary disk and calculate the composition of ices ultimately accreted in the envelope of HD 189733b on its migration pathway. This allows us to reproduce the observed volatile abundances by adjusting the mass of ices vaporized in the envelope. The predicted elemental mixing ratios should be 0.15--0.3 times solar in the envelope of HD 189733b if they are fitted to the recent O and C determinations. However, our fit to the minimum mass of heavy elements predicted by internal structure models gives elemental abundances that are 1.2--2.4 times oversolar in the envelope of HD189733b. We propose that the most likely cause of this discrepancy is irradiation from the central star leading to development of a radiative zone in the planet's outer envelope which would induce gravitational settling of elements. Hence, all strongly irradiated extrasolar planets should present subsolar abundances of volatiles. We finally predict that the abundances of nitrogen (N), sulfur (S) and phosphorus (P) are of $\sim$ $2.8 \times 10^{-5}$, $5.3 \times 10^{-6}$ and $1.8 \times 10^{-7}$ relative to H$_2$, respectively in the atmosphere of HD 189733b.

We monitor the star HD 149026 and its Saturn-mass planet at 8.0 micron over slightly more than half an orbit using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We find an increase of 0.0227% +/- 0.0066% (3.4 sigma significance) in the combined planet-star flux during this interval. The minimum flux from the planet is 45% +/- 19% of the maximum planet flux, corresponding to a difference in brightness temperature of 480 +/- 140 K between the two hemispheres. We derive a new secondary eclipse depth of 0.0411% +/- 0.0076% in this band, corresponding to a dayside brightness temperature of 1440 +/- 150 K. Our new secondary eclipse depth is half that of a previous measurement (3.0 sigma difference) in this same bandpass by Harrington et al. (2007). We re-fit the Harrington et al. (2007) data and obtain a comparably good fit with a smaller eclipse depth that is consistent with our new value. In contrast to earlier claims, our new eclipse depth suggests that this planet's dayside emission spectrum is relatively cool, with an 8 micron brightness temperature that is less than the maximum planet-wide equilibrium temperature. We measure the interval between the transit and secondary eclipse and find that that the secondary eclipse occurs 20.9 +7.2 / -6.5 minutes earlier (2.9 sigma) than predicted for a circular orbit, a marginally significant result. This corresponds to e*cos(omega) = -0.0079 +0.0027 / -0.0025 where e is the planet's orbital eccentricity and omega is the argument of pericenter.

We present global, three-dimensional numerical simulations of HD 189733b and HD 209458b that couple the atmospheric dynamics to a realistic representation of non-gray cloud-free radiative transfer. The model, which we call the Substellar and Planetary Atmospheric Radiation and Circulation (SPARC) model, adopts the MITgcm for the dynamics and uses the radiative model of McKay, Marley, Fortney, and collaborators for the radiation. Like earlier work with simplified forcing, our simulations develop a broad eastward equatorial jet, mean westward flow at higher latitudes, and substantial flow over the poles at low pressure. For HD 189733b, our simulations without TiO and VO opacity can explain the broad features of the observed 8 and 24-micron light curves, including the modest day-night flux variation and the fact that the planet/star flux ratio peaks before the secondary eclipse. Our simulations also provide reasonable matches to the Spitzer secondary-eclipse depths at 4.5, 5.8, 8, 16, and 24 microns and the groundbased upper limit at 2.2 microns. However, we substantially underpredict the 3.6-micron secondary-eclipse depth, suggesting that our simulations are too cold in the 0.1-1 bar region. Predicted temporal variability in secondary-eclipse depths is ~1% at Spitzer bandpasses, consistent with recent observational upper limits at 8 microns. We also show that nonsynchronous rotation can significantly alter the jet structure. For HD 209458b, we include TiO and VO opacity; these simulations develop a hot (>2000 K) dayside stratosphere. Despite this stratosphere, we do not reproduce current Spitzer photometry of this planet. Light curves in Spitzer bandpasses show modest phase variation and satisfy the observational upper limit on day-night phase variation at 8 microns. (abridged)

We present new Spitzer observations of the phase variation of the hot Jupiter HD 189733b in the MIPS 24 micron bandpass, spanning the same part of the planet's orbit as our previous observations in the IRAC 8 micron bandpass (Knutson et al. 2007). We find that the minimum hemisphere-averaged flux from the planet in this bandpass is 76 +/- 3% of the maximum flux; this corresponds to minimum and maximum hemisphere-averaged brightness temperatures of 984 +/- 48 K and 1220 +/- 47 K, respectively. The planet reaches its maximum flux at an orbital phase of 0.396 +/- 0.022, corresponding to a hot region shifted 20-30 degrees east of the substellar point. Because tidally locked hot Jupiters would have enormous day-night temperature differences in the absence of winds, the small amplitude of the observed phase variation indicates that the planet's atmosphere efficiently transports thermal energy from the day side to the night side at the 24 micron photosphere, leading to modest day-night temperature differences. The similarities between the 8 and 24 micron phase curves for HD 189733b lead us to conclude that the circulation on this planet behaves in a fundamentally similar fashion across the range of pressures sensed by these two wavelengths. One-dimensional radiative transfer models indicate that the 8 micron band should probe pressures 2-3 times greater than at 24 micron, although the uncertain methane abundance complicates the interpretation. If these two bandpasses do probe different pressures, it would indicate that the temperature varies only weakly between the two sensed depths, and hence that the atmosphere is not convective at these altitudes. (abridged)

We present global, three-dimensional numerical simulations of the atmospheric circulation on HD 209458b and HD 189733b and calculate the infrared spectra and light curves predicted by these simulations, which we compare with available observations. Radiative heating/cooling is parameterized with a simplified Newtonian relaxation scheme. Our simulations develop day-night temperature contrasts that vary strongly with pressure. At low pressure (<10 mbar), air flows from the substellar point toward the antistellar point, both along the equator and over the poles. At deeper levels, the flow develops an eastward equatorial jet with speeds of 3-4 km/sec, with weaker westward flows at high latitudes. This basic flow pattern is robust to variations in model resolution, gravity, radiative time constant, and initial temperature structure. Nightside spectra show deep absorption bands of H2O, CO, and/or CH4, whereas on the dayside these absorption bands flatten out or even flip into emission. This results from the strong effect of dynamics on the vertical temperature-pressure structure; the temperature decreases strongly with altitude on the nightside but becomes almost isothermal on the dayside. In Spitzer bandpasses, our predicted planet-to-star flux ratios vary by a factor of 2-10 with orbital phase, depending on the wavelength and chemistry. For HD 189733b, where a detailed 8-micron light curve has been obtained, we correctly produce the observed phase offset of the flux maximum, but we do not explain the flux minimum and we overpredict the total flux variation. This discrepancy likely results from the simplifications inherent in the Newtonian relaxation scheme and provides motivation for incorporating realistic radiative transfer in future studies.

Hot Jupiters are new laboratories for the physics of giant planet atmospheres. Subject to unusual forcing conditions, the circulation regime on these planets may be unlike anything known in the Solar System. Characterizing the atmospheric circulation of hot Jupiters is necessary for reliable interpretation of the multifaceted data currently being collected on these planets. We discuss several fundamental concepts of atmospheric dynamics that are likely central to obtaining a solid understanding of these fascinating atmospheres. A particular effort is made to compare the various modeling approaches employed so far to address this challenging problem.

We present three-dimensional numerical simulations of the atmospheric dynamics of close-orbiting planets such as HD209458b. Our simulations show that winds of several km/sec and day-night temperature differences reaching 500-1000 K are possible at and above the photosphere. The circulation takes the form of a broad superrotating (eastward) equatorial jet. At 0.1-1 bar, the superrotation blows the hottest regions of the atmosphere downwind by 60 degrees of longitude, but at lower pressures the temperature pattern tracks the stellar illumination. We predict factors of several variation in the infrared flux received at Earth throughout an orbital cycle; if the photosphere is deep enough (>50-100 mbar pressure), the peak infrared emission should lead the time of secondary eclipse by 10 hours or more. Dynamics plays a key role in shaping the spectrum, clouds, chemistry, and long-term planetary evolution.

"Hot Jupiter" extrasolar planets are expected to be tidally locked because they are close (<0.05 astronomical units, where 1 AU is the average Sun-Earth distance) to their parent stars, resulting in permanent daysides and nightsides. By observing systems where the planet and star periodically eclipse each other, several groups have been able to estimate the temperatures of the daysides of these planets. A key question is whether the atmosphere is able to transport the energy incident upon the dayside to the nightside, which will determine the temperature at different points on the planet's surface. Here we report observations of HD 189733, the closest of these eclipsing planetary systems, over half an orbital period, from which we can construct a 'map' of the distribution of temperatures. We detected the increase in brightness as the dayside of the planet rotated into view. We estimate a minimum brightness temperature of 973 +/- 33 K and a maximum brightness temperature of 1212 +/- 11 K at a wavelength of 8 microns, indicating that energy from the irradiated dayside is efficiently redistributed throughout the atmosphere, in contrast to a recent claim for another hot Jupiter. Our data indicate that the peak hemisphere-integrated brightness occurs 16$\pm$6 degrees before opposition, corresponding to a hot spot shifted east of the substellar point. The secondary eclipse (when the planet moves behind the star) occurs 120 +/- 24 s later than predicted, which may indicate a slightly eccentric orbit.

We explore the infrared spectrum of a three-dimensional dynamical model of planet HD209458b as a function of orbital phase. The dynamical model predicts day-side atmospheric pressure-temperature profiles that are much more isothermal at pressures less than 1 bar than one-dimensional radiative-convective models have found. The resulting day-side thermal spectra are very similar to a blackbody, and only weak water absorption features are seen at short wavelengths. The dayside emission is consequently in significantly better agreement with ground-based and space-based secondary eclipse data than any previous models, which predict strong flux peaks and deep absorption features. At other orbital phases, absorption due to carbon monoxide and methane is also predicted. We compute the spectra under two treatments of atmospheric chemistry: one uses the predictions of equilibrium chemistry, and the other uses non-equilibrium chemistry, which ties the timescales of methane and carbon monoxide chemistry to dynamical timescales. As a function of orbital phase, we predict planet-to-star flux ratios for standard infrared bands and all Spitzer Space Telescope bands. In Spitzer bands, we predict 2-fold to 15-fold variation in planetary flux as a function of orbital phase with equilibrium chemistry, and 2-fold to 4-fold variation with non-equilibrium chemistry. Variation is generally more pronounced in bands from 3-10 $\mu$m than at longer wavelengths. The orbital phase of maximum thermal emission in infrared bands is 15--45 orbital degrees before the time of secondary eclipse. Changes in flux as a function of orbital phase for HD209458b should be observable with Spitzer, given the previously acheived observational error bars.

We present a method that employs the secondary eclipse light curves of transiting extrasolar planets to probe the spatial variation of their thermal emission. This technique permits an observer to resolve the surface of the planet without the need to spatially resolve its central star. We evaluate the feasibility of this technique for the HD 209458 system [..]. We consider two representations of the planetary thermal emission; a simple model parameterized by a sinusoidal dependence on longitude and latitude, as well as the results of a three-dimensional dynamical simulation of the planetary atmosphere previously published by Cooper & Showman. We find that observations of the secondary eclipse light curve are most sensitive to a longitudinal offset in the geometric and photometric centroids of the hemisphere of the planet visible near opposition. To quantify this signal, we define a new parameter, the ``uniform time offset,'' which measures the time lag between the observed secondary eclipse and that predicted by a planet with a uniform surface flux distribution. We compare the predicted amplitude of this parameter for HD 209458 with the precision with which it could be measured with IRAC. We find that IRAC observations at 3.6um a single secondary eclipse should permit sufficient precision to confirm or reject the Cooper & Showman model of the surface flux distribution for this planet. We quantify the signal-to-noise ratio for this offset in the remaining IRAC bands (4.5um, 5.8um, and 8.0um), and find that a modest improvement in photometric precision (as might be realized through observations of several eclipse events) should permit a similarly robust detection.

Chemical equilibrium considerations suggest that, assuming solar elemental abundances, carbon on HD 209458b is sequestered primarily as carbon monoxide (CO) and methane (CH4). The relative mole fractions of CO(g) and CH4(g) in chemical equilibrium are expected to vary greatly according to variations in local temperature and pressure. We show, however, that in the p = 1--1000 mbar range, chemical equilibrium does not hold. To explore disequilibrium effects, we couple the chemical kinetics of CO and CH4 to a three-dimensional numerical model of HD 209458b's atmospheric circulation. These simulations show that vigorous dynamics caused by uneven heating of this tidally locked planet homogenize the CO and CH4 concentrations at p < 1 bar, even in the presence of lateral temperature variations of ~500--1000 K. In the 1--1000 mbar pressure range, we find that over 98% of the carbon is in CO. This is true even in cool regions where CH4 is much more stable thermodynamically. Our work shows furthermore that planets 300--500 K cooler than HD 209458b can also have abundant CO in their upper layers due to disequilibrium effects. We demonstrate several interesting observational consequences of these results.

We calculate the meteorology of the close-in transiting extrasolar planet HD 209458b using a global, three-dimensional atmospheric circulation model. Dynamics are driven by perpetual irradiation of one hemisphere of this tidally locked planet. The simulation predicts global temperature contrasts of ~500 K at the photosphere and the development of a steady superrotating jet. The jet extends from the equator to mid-latitudes and from the top model layer at 1 mbar down to 10 bars at the base of the heated region. Wind velocities near the equator exceed 4 km/s at 300 mbar. The hottest regions of the atmosphere are blown downstream from the substellar point by 60 degrees of longitude. We predict from these results a factor of ~2 ratio between the maximum and minimum observed radiation from the planet over a full orbital period, with peak infrared emission preceding the time of the secondary eclipse by ~14 hours.

We examine the properties of the atmospheres of extrasolar giant planets at orbital distances smaller than 0.1 AU from their stars. We show that these ``51Peg b-like'' planets are rapidly synchronized by tidal interactions, but that small departures from synchronous rotation can occur because of fluid-dynamical torques within these planets. Previous radiative-transfer and evolution models of such planets assume a homogeneous atmosphere. Nevertheless, we show using simple arguments that, at the photosphere, the day-night temperature difference and characteristic wind speeds may reach ~500 K and ~2 km/s, respectively. Substantial departures from chemical equilibrium are expected. The cloud coverage depends sensitively on the dynamics; clouds could exist predominantly either on the dayside or nightside, depending on the circulation regime. Radiative-transfer models that assume homogeneous conditions are therefore inadequate in describing the atmospheric properties of 51Peg b-like planets. We present preliminary three-dimensional, nonlinear simulations of the atmospheric circulation of HD209458b that indicate plausible patterns for the circulation and generally agree with our simpler estimates. Furthermore, we show that kinetic energy production in the atmosphere can lead to the deposition of substantial energy in the interior, with crucial consequences for the evolution of these planets. Future measurements of reflected and thermally-emitted radiation from these planets will help test our ideas.

About one-quarter of the extrasolar giant planets discovered so far have orbital distances smaller than 0.1 AU. These ``51Peg b-like'' planets can now be directly characterized, as shown by the planet transiting in front the star HD209458. We review the processes that affect their evolution. We apply our work to the case of HD209458b, whose radius has been recently measured. We argue that its radius can be reproduced only when the deep atmosphere is assumed to be unrealistically hot. When using more realistic atmospheric temperatures, an energy source appears to be missing in order to explain HD209458b's large size. The most likely source of energy available is not in the planet's spin or orbit, but in the intense radiation received from the parent star. We show that the radius of HD209458b can be reproduced if a small fraction (~1%) of the stellar flux is transformed into kinetic energy in the planetary atmosphere and subsequently converted to thermal energy by dynamical processes at pressures of tens of bars.