James Y-K. Cho

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.

Data from the Galileo Probe, collected during its descent into Jupiter's atmosphere, is used to obtain a vertical profile of the zonal wind from ~0.5 bar (upper troposphere) to ~0.1 µbar (lower thermosphere) at the probe entry site. This is accomplished by constructing a map of gravity wave Lomb‐Scargle periodograms as a function of altitude. The profile obtained from the map indicates that the wind speed above the visible cloud deck increases with height to ~150 m s^‐1 and then levels off at this value over a broad altitude range. The location of the turbopause, as a region of wide wave spectrum, is also identified from the map. In addition, a cross‐equatorial oscillation of a jet, which has previously been linked to the QQO in the stratosphere, is suggested by the profile.

We investigate baroclinic instability in flow conditions relevant to hot extrasolar planets. The instability is important for transporting and mixing heat, as well as for influencing large-scale variability on the planets. Both linear normal mode analysis and non-linear initial value calculations are carried out -- focusing on the freely-evolving, adiabatic situation. Using a high-resolution general circulation model (GCM) which solves the traditional primitive equations, we show that large-scale jets similar to those observed in current GCM simulations of hot extrasolar giant planets are likely to be baroclinically unstable on a timescale of few to few tens of planetary rotations, generating cyclones and anticyclones that drive weather systems. The growth rate and scale of the most unstable mode obtained in the linear analysis are in qualitative, good agreement with the full non-linear calculations. In general, unstable jets evolve differently depending on their signs (eastward or westward), due to the change in sign of the jet curvature. For jets located at or near the equator, instability is strong at the flanks -- but not at the core. Crucially, the instability is either poorly or not at all captured in simulations with low resolution and/or high artificial viscosity. Hence, the instability has not been observed or emphasized in past circulation studies of hot extrasolar planets.

We elucidate the interplay between Newtonian thermal relaxation and numerical dissipation, of several different origins, in flow simulations of hot extrasolar planet atmospheres. Currently, a large range of Newtonian relaxation, or "cooling," times (~10 days to ~1 hr) is used among different models and within a single model over the model domain. In this study, we demonstrate that a short relaxation time (much less than the planetary rotation time) leads to a large amount of unphysical, grid-scale oscillations that contaminate the flow field. These oscillations force the use of an excessive amount of artificial viscosity to quench them and prevent the simulation from "blowing up." Even if the blow-up is prevented, such simulations can be highly inaccurate because they are either severely overdissipated or underdissipated, and are best discarded in these cases. Other numerical stability and time step size enhancers (e.g., Robert-Asselin filter or semi-implicit time-marching schemes) also produce similar, but less excessive, damping. We present diagnostics procedures to choose the "optimal" simulation and discuss implications of our findings for modeling hot extrasolar planet atmospheres.

We elucidate the interplay between Newtonian thermal relaxation and numerical dissipation, of several different origins, in flow simulations of hot extrasolar planet atmospheres. Currently, a large range of Newtonian relaxation, or "cooling", times (~10 days to ~1 hour) is used among different models and within a single model over the model domain. In this study we demonstrate that a short relaxation time (much less than the planetary rotation time) leads to a large amount of unphysical, grid-scale oscillations that contaminate the flow field. These oscillations force the use of an excessive amount of artificial viscosity to quench them and prevent the simulation from "blowing up". Even if the blow-up is prevented, such simulations can be highly inaccurate because they are either severely over-dissipated or under-dissipated, and are best discarded in these cases. Other numerical stability and timestep size enhancers (e.g., Robert-Asselin filter or semi-implicit time-marching schemes) also produce similar, but less excessive, damping. We present diagnostics procedures to choose the "optimal" simulation and discuss implications of our findings for modeling hot extrasolar planet atmospheres.

We report the detection of 24 micron variations from the planet-hosting upsilon Andromedae system consistent with the orbital periodicity of the system's innermost planet, upsilon And b. We find a peak-to-valley phase curve amplitude of 0.00130 times the mean system flux. Using a simple model with two hemispheres of constant surface brightness and assuming a planetary radius of 1.3 Jupiter radii gives a planetary temperature contrast of >900 K and an orbital inclination of >28 degrees. We further report the largest phase offset yet observed for an extrasolar planet: the flux maximum occurs ~80 degrees before phase 0.5. Such a large phase offset is difficult to reconcile with most current atmospheric circulation models. We improve on earlier observations of this system in several important ways: (1) observations of a flux calibrator star demonstrate the MIPS detector is stable to 10^-4 on long timescales, (2) we note that the background light varies systematically due to spacecraft operations, precluding use of this background as a flux calibrator (stellar flux measured above the background is not similarly affected), and (3) we calibrate for flux variability correlated with motion of the star on the MIPS detector. A reanalysis of our earlier observations of this system is consistent with our new result.

We use a general circulation model to study the three-dimensional (3D) flow and temperature distributions of atmospheres on tidally synchronized extrasolar planets. In this work, we focus on the sensitivity of the evolution to the initial flow state, which has not received much attention in 3D modeling studies. We find that different initial states lead to markedly different distributions-even under the application of strong forcing (large day-night temperature difference with a short "thermal drag time") that may be representative of close-in planets. This is in contrast with the results or assumptions of many published studies. In general, coherent jets and vortices (and their associated temperature distributions) characterize the flow, and they evolve differently in time, depending on the initial condition. If the coherent structures reach a quasi-stationary state, their spatial locations still vary. The result underlines the fact that circulation models are currently unsuitable for making quantitative predictions (e. g., location and size of a "hot spot") without better constrained, and well-posed, initial conditions.

We use a general circulation model to study the three-dimensional (3-D) flow and temperature distributions of atmospheres on tidally synchronized extrasolar planets. In this work, we focus on the sensitivity of the evolution to the initial flow state, which has not received much attention in 3-D modeling studies. We find that different initial states lead to markedly different distributions-even under the application of strong forcing (large day-night temperature difference with a short "thermal drag time") that may be representative of close-in planets. This is in contrast with the results or assumptions of many published studies. In general, coherent jets and vortices (and their associated temperature distributions) characterize the flow, and they evolve differently in time, depending on the initial condition. If the coherent structures reach a quasi- stationary state, their spatial locations still vary. The result underlines the fact that circulation models are currently unsuitable for making quantitative predictions (e.g., location and size of a "hot spot") without better constrained, and well posed, initial conditions.

We study the effects of gravity waves, or g-modes, on hot extrasolar planets. These planets are expected to possess stably-stratified atmospheres, which support gravity waves. In this paper, we review the derivation of the equation that governs the linear dynamics of gravity waves and describe its application to a hot extrasolar planet, using HD209458 b as a generic example. We find that gravity waves can exhibit a wide range of behaviors, even for a single atmospheric profile. The waves can significantly accelerate or decelerate the background mean flow, depending on the difference between the wave phase and mean flow speeds. In addition, the waves can provide significant heating (~100 to ~1000 K per planetary rotation), especially to the region of the atmosphere above about 10 scale heights from the excitation region. Furthermore, by propagating horizontally, gravity waves provide a mechanism for transporting momentum and heat from the dayside of a tidally locked planet to its nightside. We discuss work that needs to be undertaken to incorporate these effects in current atmosphere models of extrasolar planets.