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.

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

We compare five general circulation models (GCMs) which have been recently used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and PEQMOD), under three test cases useful for assessing model convergence and accuracy. Such a broad, detailed intercomparison has not been performed thus far for extrasolar planets study. The models considered all solve the traditional primitive equations, but employ different numerical algorithms or grids (e.g., pseudospectral and finite volume, with the latter separately in longitude-latitude and `cubed-sphere' grids). The test cases are chosen to cleanly address specific aspects of the behaviors typically reported in hot extrasolar planet simulations: 1) steady-state, 2) nonlinearly evolving baroclinic wave, and 3) response to fast timescale thermal relaxation. When initialized with a steady jet, all models maintain the steadiness, as they should -- except MITgcm in cubed-sphere grid. A very good agreement is obtained for a baroclinic wave evolving from an initial instability in pseudospectral models (only). However, exact numerical convergence is still not achieved across the pseudospectral models: amplitudes and phases are observably different. When subject to a typical `hot-Jupiter'-like forcing, all five models show quantitatively different behavior -- although qualitatively similar, time-variable, quadrupole-dominated flows are produced. Hence, as have been advocated in several past studies, specific quantitative predictions (such as the location of large vortices and hot regions) by GCMs should be viewed with caution. Overall, in the tests considered here, pseudospectral models in pressure coordinate (PEBOB and PEQMOD) perform the best and MITgcm in cubed-sphere grid performs the worst.

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