Source author record

Adam P. Showman

Adam P. Showman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

75works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

75 published item(s)

preprint2022arXiv

Influences of internal forcing on atmospheric circulations of irradiated giant planets

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

preprint2022arXiv

Particle deposition on the saturnian satellites from ephemeral cryovolcanism on Enceladus

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

preprint2021arXiv

Atmospheric circulation of brown dwarfs and directly imaged exoplanets driven by cloud radiative feedback: effects of rotation

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

preprint2021arXiv

Atmospheric circulation of brown dwarfs and directly imaged exoplanets driven by cloud radiative feedback: global and equatorial dynamics

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

preprint2021arXiv

Thermal Phase Curves of XO-3b: an Eccentric Hot Jupiter at the Deuterium Burning Limit

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

preprint2020arXiv

How well do we understand the belt/zone circulation of Giant Planet 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]

preprint2020arXiv

Smaller than expected bright-spot offsets in Spitzer phase curves of the hot Jupiter Qatar-1b

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~μ$m and $0.30 \pm 0.02 \%$ at $4.5~μ$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~μ$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~μ$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.

preprint2019arXiv

Comparison of the deep atmospheric dynamics of Jupiter and Saturn in light of the Juno and Cassini gravity measurements

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.

preprint2019arXiv

Evaluating Climate Variability Of The Canonical Hot Jupiters Hd 189733b & Hd 209458b Through Multi-epoch Eclipse Observations

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.

preprint2019arXiv

Temporal Variability in Hot Jupiter Atmospheres

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.

preprint2016arXiv

$Extrasolar~Storms$: Pressure-dependent Changes In Light Curve Phase In Brown Dwarfs From Simultaneous $Hubble$ and $Spitzer$ Observations

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.

preprint2016arXiv

3.6 and 4.5 $μ$m ${\it Spitzer}$ Phase Curves of the Highly-Irradiated Hot Jupiters WASP-19b and HAT-P-7b

We analyze full-orbit phase curve observations of the transiting hot Jupiters WASP-19b and HAT-P-7b at 3.6 and 4.5 $μ$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 $μ$m, which are consistent with a single blackbody with effective temperature $2372 \pm 60$ K. The measured 3.6 and 4.5 $μ$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 $μ$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σ$) 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.

preprint2016arXiv

Atmospheric Circulation of Hot Jupiters: Dayside-Nightside Temperature Differences

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.

preprint2016arXiv

HST hot-Jupiter transmission spectral survey: Clear skies for cool Saturn WASP-39b

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.

preprint2016arXiv

No Thermal Inversion and a Solar Water Abundance for the Hot Jupiter HD209458b from HST WFC3 Emission Spectroscopy

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.

preprint2016arXiv

The atmospheric circulation of a nine-hot Jupiter sample: Probing circulation and chemistry over a wide phase space

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.

preprint2016arXiv

Transitions in the cloud composition of hot Jupiters

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.

preprint2015arXiv

3.6 and 4.5 $μ$m Phase Curves of the Highly-Irradiated Eccentric Hot Jupiter WASP-14b

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 $μ$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 $μ$m, respectively, are both consistent with a blackbody temperature of $2402\pm 35$ K. We place a $2σ$ upper limit on the nightside flux at 3.6 $μ$m and find it to be $9\%\pm 1\%$ of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 $μ$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 $μ$m, while under-predicting the nightside emission at 4.5 $μ$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.

preprint2015arXiv

A continuum from clear to cloudy hot-Jupiter exoplanets without primordial water depletion

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 μ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.

preprint2015arXiv

A Detection of Water in the Transmission Spectrum of the Hot Jupiter WASP-12b and Implications for its Atmospheric Composition

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 $μ$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$σ$ confidence). This detection marks the first high-confidence, spectroscopic identification of a molecule in the atmosphere of WASP-12b. The retrieved 1$σ$ water volume mixing ratio is between $10^{-5}-10^{-2}$, which is consistent with C/O $>$ 1 to within 2$σ$. 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\,σ$ and rule out a carbon-rich atmosphere composition (C/O$>1$) at $>3σ$ 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.

preprint2015arXiv

A Semi-Analytical Model of Visible-Wavelength Phase Curves of Exoplanets and Applications to Kepler-7 b and Kepler-10 b

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)

preprint2015arXiv

Atmospheric dynamics of terrestrial exoplanets over a wide range of orbital and atmospheric parameters

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.

preprint2015arXiv

Discovery of Rotational Modulations in the Planetary-Mass Companion 2M1207b: Intermediate Rotation Period and Heterogeneous Clouds in a Low Gravity Atmosphere

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.

preprint2015arXiv

Spitzer Secondary Eclipse Observations of Five Cool Gas Giant Planets and Empirical Trends in Cool Planet Emission Spectra

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.

preprint2015arXiv

Spitzer Secondary Eclipses of the Dense, Modestly-irradiated, Giant Exoplanet HAT-P-20b Using Pixel-Level Decorrelation

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.

preprint2015arXiv

The Atmospheric Circulation of the Hot Jupiter WASP-43b: Comparing Three-Dimensional Models to Spectrophotometric Data

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.

preprint2015arXiv

Three-dimensional Atmospheric Circulation of Warm and Hot Jupiters: Effects of Orbital Distance, Rotation Period, and Non-Synchronous Rotation

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.

preprint2014arXiv

A Precise Water Abundance Measurement for the Hot Jupiter WASP-43b

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 $σ$ 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.

preprint2014arXiv

Atmospheric Characterization of the Hot Jupiter Kepler-13Ab

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.

preprint2014arXiv

Atmospheric Circulation of Brown Dwarfs: Jets, Vortices, and Time Variability

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.

preprint2014arXiv

Atmospheric Circulation of Eccentric Hot Jupiter HAT-P-2b

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.

preprint2014arXiv

Atmospheric Dynamics of Hot Exoplanets

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.

preprint2014arXiv

Characterization of the atmosphere of the hot Jupiter HAT-P-32Ab and the M-dwarf companion HAT-P-32B

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.

preprint2014arXiv

Constraints on the Atmospheric Circulation and Variability of the Eccentric Hot Jupiter XO-3b

We report secondary eclipse photometry of the hot Jupiter XO-3b in the 4.5~$μ$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.

preprint2014arXiv

HST Rotational Spectral Mapping of Two L-Type Brown Dwarfs: Variability In and Out of Water Bands Indicates High-Altitude Haze Layers

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 $μ$m and 1.7 $μ$m. We find that the water absorption bands of the two L5 dwarfs at 1.15 $μ$m and 1.4 $μ$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 $μ$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.

preprint2014arXiv

Spitzer and z' Secondary Eclipse Observations of the Highly Irradiated Transiting Brown Dwarf KELT-1b

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.

preprint2014arXiv

The 4.5 $μ$m full-orbit phase curve of the hot Jupiter HD 209458b

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 $μ$m emission was recorded with $Spitzer$/IRAC. Our study revises the previous 4.5 $μ$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.

preprint2014arXiv

The atmospheric circulation of the super Earth GJ 1214b: Dependence on composition and metallicity

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.

preprint2014arXiv

Thermal structure of an exoplanet atmosphere from phase-resolved emission spectroscopy

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.

preprint2014arXiv

Unveiling the atmospheres of giant exoplanets with an EChO-class mission

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{μm}$), high spectral resolution ($Δλ/λ>300$ below $5\,\rm{μm}$ and $Δλ/λ>30$ above $5\,\rm{μ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.

preprint2013arXiv

3D mixing in hot Jupiter atmospheres. I. application to the day/night cold trap in HD 209458b

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.

preprint2013arXiv

Atmospheric circulation of hot Jupiters: insensitivity to initial conditions

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.

preprint2013arXiv

Atmospheric Circulation of Terrestrial Exoplanets

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.

preprint2013arXiv

Atmospheric dynamics of brown dwarfs and directly imaged giant planets

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.

preprint2013arXiv

Atmospheric Heat Redistribution on Hot Jupiters

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.

preprint2013arXiv

Cloud Features and Zonal Wind Measurements of Saturn's Atmosphere as Observed by Cassini/VIMS

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.

preprint2013arXiv

Doppler Signatures of the Atmospheric Circulation on Hot Jupiters

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.

preprint2013arXiv

Infrared Transmission Spectroscopy of the Exoplanets HD209458b and XO-1b Using the Wide Field Camera-3 on the Hubble Space Telescope

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.

preprint2013arXiv

Meteorology of Jupiter's Equatorial Hot Spots and Plumes from Cassini

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.

preprint2013arXiv

Orbital Phase Variations of the Eccentric Giant Planet HAT-P-2b

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 μm bands of the Spitzer Space Telescope. The 3.6 and 4.5 μ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 μ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 μ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.

preprint2013arXiv

Power Spectral Analysis of Jupiter's Clouds and Kinetic Energy from Cassini

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 μ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.

preprint2013arXiv

Secondary Eclipse Photometry of the Exoplanet WASP-5b with Warm Spitzer

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.

preprint2013arXiv

The Evolving Flow of Jupiter's White Ovals and Adjacent Cyclones

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.

preprint2013arXiv

Three-dimensional atmospheric circulation of hot Jupiters on highly eccentric orbits

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.

preprint2013arXiv

Velocity and Vorticity Measurements of Jupiter's Great Red Spot Using Automated Cloud Feature Tracking

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.

preprint2013arXiv

Warm Spitzer and Palomar Near-IR Secondary Eclipse Photometry of Two Hot Jupiters: WASP-48b and HAT-P-23b

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.

preprint2013arXiv

Warm Spitzer Photometry of Three Hot Jupiters: HAT-P-3b, HAT-P-4b and HAT-P-12b

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$σ$ 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.

preprint2012arXiv

3.6 and 4.5 Micron Phase Curves and Evidence for Non-Equilibrium Chemistry in the Atmosphere of Extrasolar Planet HD 189733b

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]

preprint2012arXiv

Spitzer/MIPS 24 micron Observations of HD 209458b: 3 eclipses, 2.5 transits, and a Phase Curve Corrupted by Instrumental Sensitivity Variations

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]

preprint2012arXiv

The EVIL-MC Model for Ellipsoidal Variations of Planet-Hosting Stars and Applications to the HAT-P-7 System

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.

preprint2012arXiv

Vertical Atmospheric Structure in a Variable Brown Dwarf: Pressure-dependent Phase Shifts in Simultaneous Hubble Space Telescope-Spitzer Light Curves

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.

preprint2011arXiv

Disequilibrium Carbon, Oxygen, and Nitrogen Chemistry in the Atmospheres of HD 189733b and HD 209458b

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.

preprint2011arXiv

Equatorial superrotation on tidally locked exoplanets

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.

preprint2011arXiv

Warm Spitzer Photometry of XO-4b, HAT-P-6b and HAT-P-8b

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μm and 4.5μm for each target. XO-4b exhibits a stronger eclipse depth at 4.5μm than at 3.6μm, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6μ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(ω) for HAT-P-6b, HAT-P-8b and XO-4b and find that the values are consistent with circular orbits.

preprint2010arXiv

Atmospheric Circulation of Eccentric Hot Neptune GJ436b

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.

preprint2010arXiv

Scaling laws for convection and jet speeds in the giant planets

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.

preprint2010arXiv

Secondary Eclipse Photometry of WASP-4b with Warm Spitzer

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.

preprint2010arXiv

Warm Spitzer Photometry of the Transiting Exoplanets CoRoT-1 and CoRoT-2 at Secondary Eclipse

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.

preprint2009arXiv

Atmospheric circulation of hot Jupiters: Coupled radiative-dynamical general circulation model simulations of HD 189733b and HD 209458b

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)

preprint2009arXiv

Detecting the Wind-Driven Shapes of Extrasolar Giant Planets from Transit Photometry

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.

preprint2009arXiv

Elemental abundances and minimum mass of heavy elements in the envelope of HD 189733b

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.

preprint2009arXiv

Generation of equatorial jets by large-scale latent heating on the giant planets

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.

preprint2009arXiv

The 8 Micron Phase Variation of the Hot Saturn HD 149026b

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.

preprint2007arXiv

A map of the day-night contrast of the extrasolar planet HD 189733b

"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.

preprint2002arXiv

Atmospheric Circulation and Tides of "51Peg b-like" Planets

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.