No Arabic abstract
Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady state can be reached. Aims. A consistent physical understanding of superrotation is needed for arbitrary drag and radiative timescales, and the relevance of considering linear steady states needs to be assessed. Methods. We obtain an analytical expression for the structure, frequency and decay rate of propagating waves in hot Jupiter atmospheres around a state at rest in the 2D shallow water beta plane limit. We solve this expression numerically and confirm the robustness of our results with a 3D linear wave algorithm. We then compare with 3D simulations of hot Jupiter atmospheres and study the non linear momentum fluxes. Results. We show that under strong day night heating the dynamics does not transit through a linear steady state when starting from an initial atmosphere in solid body rotation. We further show that non linear effects favour the initial spin up of superrotation and that the acceleration due to the vertical component of the eddy momentum flux is critical to the initial development of superrotation. Conclusions. Overall, we describe the initial phases of the acceleration of superrotation, including consideration of differing radiative and drag timescales, and conclude that eddy-momentum driven superrotating equatorial jets are robust, physical phenomena in simulations of hot Jupiter atmospheres.
The macroturbulent atmospheric circulation of Earth-like planets mediates their equator-to-pole heat transport. For fast-rotating terrestrial planets, baroclinic instabilities in the mid-latitudes lead to turbulent eddies that act to transport heat poleward. In this work, we derive a scaling theory for the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanet atmospheres. This theory is built on the work of Jansen & Ferrari (2013), and determines how unstable the atmosphere is to baroclinic instability (the baroclinic criticality) through a balance between the baroclinic eddy heat flux and radiative heating/cooling. We compare our scaling theory to General Circulation Model (GCM) simulations and find that the theoretical predictions for equator-to-pole temperature contrast and bulk lapse rate broadly agree with GCM experiments with varying rotation rate and surface pressure throughout the baroclincally unstable regime. Our theoretical results show that baroclinic instabilities are a strong control of heat transport in the atmospheres of Earth-like exoplanets, and our scalings can be used to estimate the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanets. These scalings can be tested by spectroscopic retrievals and full-phase light curves of terrestrial exoplanets with future space telescopes.
Radiative transfer in planetary atmospheres is usually treated in the static limit, i.e., neglecting atmospheric motions. We argue that hot Jupiter atmospheres, with possibly fast (sonic) wind speeds, may require a more strongly coupled treatment, formally in the regime of radiation-hydrodynamics. To lowest order in v/c, relativistic Doppler shifts distort line profiles along optical paths with finite wind velocity gradients. This leads to flow-dependent deviations in the effective emission and absorption properties of the atmospheric medium. Evaluating the overall impact of these distortions on the radiative structure of a dynamic atmosphere is non-trivial. We present transmissivity and systematic equivalent width excess calculations which suggest possibly important consequences for radiation transport in hot Jupiter atmospheres. If winds are fast and bulk Doppler shifts are indeed important for the global radiative balance, accurate modeling and reliable data interpretation for hot Jupiter atmospheres may prove challenging: it would involve anisotropic and dynamic radiative transfer in a coupled radiation-hydrodynamical flow. On the bright side, it would also imply that the emergent properties of hot Jupiter atmospheres are more direct tracers of their atmospheric flows than is the case for Solar System planets. Radiation-hydrodynamics may also influence radiative transfer in other classes of hot exoplanetary atmospheres with fast winds.
We show that by changing a single non-dimensional number, the thermal Rossby number, global atmospheric simulations with only axisymmetric forcing pass from an Earth-like atmosphere to a superrotating atmosphere that more resembles the atmospheres of Venus or Titan. The transition to superrotation occurs under conditions in which equatorward-propagating Rossby waves generated by baroclinic instability at intermediate and high latitudes are suppressed, which will occur when the deformation radius exceeds the planetary radius. At large thermal Rossby numbers following an initial, nearly axisymmetric phase, a global baroclinic wave of zonal wavenumber one generated by mixed barotropic-baroclinic instability dominates the eddy flux of zonal momentum. The global wave converges eastward zonal momentum to the equator and deposits westward momentum at intermediate latitudes during spinup and before superrotation emerges, and the baroclinic instability ceases once superrotation is established. A global barotropic mode of zonal wavenumber one generated by a mix of high- and low-latitude barotropic instability is responsible for maintaining superrotation in the statistically steady state. At intermediate thermal Rossby numbers, momentum flux by the global baroclinic mode is subdominant relative to smaller baroclinic modes, and thus strong superrotation does not develop.
[Abridged] A key hypothesis in the field of exoplanet atmospheres is the trend of atmospheric thermal structure with planetary equilibrium temperature. We explore this trend and report here the first statistical detection of a transition in the near-infrared (NIR) atmospheric emission between hot and ultra-hot Jupiters. We measure this transition using secondary eclipse observations and interpret this phenomenon as changes in atmospheric properties, and more specifically in terms of transition from non-inverted to inverted thermal profiles. We examine a sample of 78 hot Jupiters with secondary eclipse measurements at 3.6 {mu}m and 4.5 {mu}m measured with Spitzer Infrared Array Camera (IRAC). We measure the deviation of the data from the blackbody, which we define as the difference between the observed 4.5 {mu}m eclipse depth and that expected at this wavelength based on the brightness temperature measured at 3.6 {mu}m. We study how the deviation between 3.6 and 4.5 {mu}m changes with theoretical predictions with equilibrium temperature and incoming stellar irradiation. We reveal a clear transition in the observed emission spectra of the hot Jupiter population at 1660 +/- 100 K in the zero albedo, full redistribution equilibrium temperature. We find the hotter exoplanets have even hotter daysides at 4.5 {mu}m compared to 3.6 {mu}m, which manifests as an exponential increase in the emitted power of the planets with stellar insolation. We propose that the measured transition is a result of seeing carbon monoxide in emission due to the formation of temperature
We revisit the issue of sensitivity to initial flow and intrinsic variability in hot-Jupiter atmospheric flow simulations, originally investigated by Cho et al. (2008) and Thrastarson & Cho (2010). The flow in the lower region (~1 to 20 MPa) `dragged to immobility and uniform temperature on a very short timescale, as in Liu & Showman (2013), leads to effectively a complete cessation of variability as well as sensitivity in three-dimensional (3D) simulations with traditional primitive equations. Such momentum (Rayleigh) and thermal (Newtonian) drags are, however, ad hoc for 3D giant planet simulations. For 3D hot-Jupiter simulations, which typically already employ strong Newtonian drag in the upper region, sensitivity is not quenched if only the Newtonian drag is applied in the lower region, without the strong Rayleigh drag: in general, both sensitivity and variability persist if the two drags are not applied concurrently in the lower region. However, even when the drags are applied concurrently, vertically-propagating planetary waves give rise to significant variability in the ~0.05 to 0.5 MPa region, if the vertical resolution of the lower region is increased (e.g. here with 1000 layers for the entire domain). New observations on the effects of the physical setup and model convergence in `deep atmosphere simulations are also presented.