No Arabic abstract
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)
Atmospheric circulation on tidally-locked exoplanets is driven by the absorption and reradiation of heat from the host star. They are natural heat engines, converting heat into mechanical energy. A steady state is possible only if there is a mechanism to dissipate mechanical energy, or if the redistribution of heat is so effective that the Carnot efficiency is driven to zero. Simulations based on primitive, equivalent-barotropic, or shallow-water equations without explicit provision for dissipation of kinetic energy and for recovery of that energy as heat, violate energy conservation. More seriously perhaps, neglect of physical sources of drag may overestimate wind speeds and rates of advection of heat from the day to the night side.
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.
We present new Spitzer observations of the phase variation of the hot Jupiter HD 189733b in the MIPS 24 micron bandpass, spanning the same part of the planets orbit as our previous observations in the IRAC 8 micron bandpass (Knutson et al. 2007). We find that the minimum hemisphere-averaged flux from the planet in this bandpass is 76 +/- 3% of the maximum flux; this corresponds to minimum and maximum hemisphere-averaged brightness temperatures of 984 +/- 48 K and 1220 +/- 47 K, respectively. The planet reaches its maximum flux at an orbital phase of 0.396 +/- 0.022, corresponding to a hot region shifted 20-30 degrees east of the substellar point. Because tidally locked hot Jupiters would have enormous day-night temperature differences in the absence of winds, the small amplitude of the observed phase variation indicates that the planets atmosphere efficiently transports thermal energy from the day side to the night side at the 24 micron photosphere, leading to modest day-night temperature differences. The similarities between the 8 and 24 micron phase curves for HD 189733b lead us to conclude that the circulation on this planet behaves in a fundamentally similar fashion across the range of pressures sensed by these two wavelengths. One-dimensional radiative transfer models indicate that the 8 micron band should probe pressures 2-3 times greater than at 24 micron, although the uncertain methane abundance complicates the interpretation. If these two bandpasses do probe different pressures, it would indicate that the temperature varies only weakly between the two sensed depths, and hence that the atmosphere is not convective at these altitudes. (abridged)
We study the feasibility of observationally constraining the rotation rate of hot Jupiters, planets that are typically assumed to have been tidally locked into synchronous rotation. We use a three-dimensional General Circulation Model to solve for the atmospheric structure of two hot Jupiters (HD 189733b and HD 209458b), assuming rotation periods that are 0.5, 1, or 2 times their orbital periods (2.2 and 3.3 days, respectively), including the effect of variable stellar heating. We compare two observable properties: 1) the spatial variation of flux emitted by the planet, measurable in orbital phase curves, and 2) the net Doppler shift in transmission spectra of the atmosphere, which is tantalizingly close to being measurable in high-resolution transit spectra. Although we find little difference between the observable properties of the synchronous and non-synchronous models of HD 189733b, we see significant differences when we compare the models of HD 209458b. In particular, the slowly rotating model of HD 209458b has an atmospheric circulation pattern characterized by westward flow and an orbital phase curve that peaks after secondary eclipse (in contrast to all of our other models), while the quickly rotating model has a net Doppler shift that is more strongly blue-shifted than the other models. Our results demonstrate that the combined use of these two techniques may be a fruitful way to constrain the rotation rate of some planets, and motivate future work on this topic.
The measurement of the light scattered from extrasolar planets informs atmospheric and formation models. With the discovery of many hot Jupiter planets orbiting nearby stars, this motivates the development of robust methods of characterisation from follow up observations. In this paper we discuss two methods for determining the planetary albedo in transiting systems. First, the most widely used method for measuring the light scattered by hot Jupiters (Collier Cameron et al.) is investigated for application for typical echelle spectra of a transiting planet system, showing that detection requires high signal-to-noise ratio data of bright planets. Secondly a new Fourier analysis method is also presented, which is model-independent and utilises the benefits of the reduced number of unknown parameters in transiting systems. This approach involves solving for the planet and stellar spectra in Fourier space by least-squares. The sensitivities of the methods are determined via Monte Carlo simulations for a range of planet-to-star fluxes. We find the Fourier analysis method to be better suited to the ideal case of typical observations of a well constrained transiting system than the Collier Cameron et al. method. We apply the Fourier analysis method for extracting the light scattered by transiting hot Jupiters from high resolution spectra to echelle spectra of HD 209458 and HD 189733. Unfortunately we are unable to improve on the previous upper limit of the planet-to-star flux for HD 209458b set by space-based observations. A 1{sigma}upper limit on the planet-to-star flux of HD 189733b is measured in the wavelength range of 558.83-599.56 nm yielding {epsilon} < 4.5 times 10-4. Improvement in the measurement of the upper limit of the planet-to-star flux of this system, with ground-based capabilities, requires data with a higher signal-to-noise ratio, and increased stability of the telescope.