We present an analysis of seven primary transit observations of the hot Neptune GJ436b at 3.6, 4.5 and 8 µm obtained with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. After correcting for systematic effects, we fitted the light curves using the Markov Chain Monte Carlo technique. Combining these new data with the EPOXI, HST and ground-based V, I, H and K s published observations, the range 0.5 − 10 µm can be covered. Due to the low level of activity of GJ436, the effect of starspots on the combination of transits at different epochs is negligible at the accuracy of the dataset. Representative climate models were calculated by using a three-dimensional, pseudo-spectral general circulation model with idealised thermal forcing. Simulated transit spectra of GJ436b were generated using line-by-line radiative transfer models including the opacities of the molecular species expected to be present in such a planetary atmosphere. A new, ab-initio calculated, linelist for hot ammonia has been used for the first time. The photometric data observed at multiple wavelengths can be interpreted with methane being the dominant absorption after molecular hydrogen, possibly with minor contributions from ammonia, water and other molecules. No clear evidence of carbon monoxide and dioxide is found from transit photometry. We discuss this result in the light of a recent paper where photochemical disequilibrium is hypothesised to interpret secondary transit photometric data. We show that the emission photometric data are not incompatible with the presence

Capsule SummaryDuring austral spring 2019 the Antarctic stratosphere experienced record-breaking warming and a near-record polar vortex weakening, resulting in predictable extreme climate conditions throughout the Southern Hemisphere through December 2019.

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.

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.

Convectively generated gravity waves (CGGWs) are important for numerical weather prediction due to their effect on the quasi-biennial oscillation (QBO) in the stratosphere. Using global ECMWF IFS simulations at TCo7999 (or 1.25 km), TCo2559 (or 3.9 km) and TCo1279 (or 7.8 km) horizontal resolutions, sensitivity of resolved CGGWs to the horizontal resolution and to the explicit versus parametrized representation of deep convection is elucidated during the westerly shear phase of the QBO. Parametrized deep convection is found to inhibit CGGWs, resulting in a twofold reduction in CGGW forcing. When deep convection is explicitly resolved, the total CGGW forcing is almost unchanged across the horizontal resolutions. However, the contribution of long and mesoscale CGGWs (with horizontal wavelengths 100 km ≤ 𝜆 h < 1,900 km) to the total CGGW forcing decreases and the contribution of smaller-scale CGGWs (with 𝜆 h < 100 km) increases as the horizontal resolution increases. At the maximum CGGW forcing altitude, at TCo7999 resolution 43% of the total CGGW forcing is due to long and mesoscale waves, whereas at TCo2559 and TCo1279 resolutions their contribution is 70% and 90%, respectively. While CGGW forcing from long and mesoscale waves is similar at TCo7999 resolution with explicit deep convection and at TCo1279 resolution with parametrized deep convection, CGGW forcing from these waves is artificially enhanced at TCo1279 and TCo2559 resolutions with explicit deep convection. This is due to the explicit deep convection being too strong and having too much variance for 100 km ≤ 𝜆 h < 1,900 km. Therefore, parametrizations of deep convection and CGGWs (to account for forcing from waves with 𝜆 h < 100 km) are required even at TCo2559 resolution. Additionally, resolved CGGW forcing at TCo7999 resolution is examined for the easterly shear phase of the QBO; similar to the westerly shear phase, the smaller-scale waves contribute >55% to the total CGGW forcing at the maximum CGGW forcing altitude.