Search

The global superrotation index S compares the integrated axial angular momentum of the atmosphere to that of a state of solid-body corotation with the underlying planet. The index S is similar to a zonal Rossby number, which suggests it may be a useful indicator of the circulation regime occupied by a planetary atmosphere. We investigate the utility of S for characterizing regimes of atmospheric circulation by running idealized Earthlike general circulation model experiments over a wide range of rotation rates Ω, 8ΩE to ΩE/512, where ΩE is Earth’s rotation rate, in both an axisymmetric and three-dimensional configuration. We compute S for each simulated circulation, and study the dependence of S on Ω. For all rotation rates considered, S is on the same order of magnitude in the 3D and axisymmetric experiments. For high rotation rates, S ≪ 1 and S ∝ Ω−2, while at low rotation rates S ≈ 1/2 = constant. By considering the limiting behavior of theoretical models for S, we show how the value of S and its local dependence on Ω can be related to the circulation regime occupied by a planetary atmosphere. Indices of S ≪ 1 and S ∝ Ω−2 define a regime dominated by geostrophic thermal wind balance, and S ≈ 1/2 = constant defines a regime where the dynamics are characterized by conservation of angular momentum within a planetary-scale Hadley circulation. Indices of S ≫ 1 and S ∝ Ω−2 define an additional regime dominated by cyclostrophic balance and strong equatorial superrotation that is not realized in our simulations.

Saturn's polar regions (polewards of ∼63° planetocentric latitude) are strongly dynamically active with zonal jets, polar cyclones and the intriguing north polar hexagon (NPH) wave. Here we analyze measurements of horizontal winds, previously obtained from Cassini images by Antuñano et al. (2015), https://doi.org/10.1002/2014je004709, to determine the spatial and spectral exchanges of kinetic energy (KE) between zonal mean zonal jets and nonaxisymmetric eddies in Saturn's polar regions. Eddies of most resolved scales generally feed KE into the eastward and westward zonal mean jets at rates between 4.3 × 10−5 and 1.4 × 10−4 W kg−1. In particular, the north polar jet (at 76°N) was being energized at a rate of ∼10−4 W kg−1, dominated by the contribution due to the zonal wavenumber m = 6 NPH wave itself. This implies that the hexagon was not being driven at this time through a barotropic instability of the north polar jet, but may suggest a significant role for baroclinic instabilities, convection or other internal energy sources for this feature. The south polar zonal mean jet KE was also being sustained by eddies in that latitude band across a wide range of m. In contrast, results indicate that the north polar vortex may have been weakly barotropically unstable at this time with eddies of low m gaining KE at the expense of the axisymmetric cyclone. However, the southern axisymmetric polar cyclone was gaining KE from non-axisymmetric components at this time, including m = 2 and its harmonics, as the elliptical distortion of the vortex may have been decaying.

We describe the life and scientific accomplishments of Professor Raymond Hide. He was a past President of the Royal Meteorological Society and a supreme example of a geophysicist much honoured in his lifetime. He covered a wide area of geophysics from geomagnetism, meteorology, geodesy, oceanography and related aspects of planetary physics. Raymond Hide was particularly known in meteorology as a founding father of geophysical fluid dynamics, especially for his experiments using a rotating cylindrical annulus to study atmospheric dynamics.

The weather layers of the gas giant planets, Jupiter and Saturn, comprise the shallow atmospheric layers that are influenced energetically by a combination of incoming solar radiation and localised latent heating of condensates, as well as by upwelling heat from their planetary interiors. They are also the most accessible regions of those planets to direct observations. Recent analyses in Oxford of cloud-tracked winds on Jupiter have demonstrated that kinetic energy is injected into the weather layer at scales comparable to the Rossby radius of deformation and cascades both upscale, mostly into the extra-tropical zonal jets, and downscale to the smallest resolvable scales in Cassini images. The large-scale flow on both Jupiter and Saturn appears to equilibrate towards a state which is close to marginal instability according to Arnol'd's 2nd stability theorem. This scenario is largely reproduced in a hierarchy of numerical models of giant planet weather layers, including relatively realistic models which seek to predict thermal and dynamical structures using a full set of parameterisations of radiative transfer, interior heat sources and even moist convection. Such models include the Jason GCM, developed in Oxford, which also represents the formation of (energetically passive) clouds of NH3, NH4SH and H2O condensates and the transport of condensable tracers. Recent results show some promise in comparison with observations from the Cassini and Juno missions, but some observed features (such as Jupiter's Great Red Spot and other compact ovals) are not yet captured spontaneously by any weather layer model. We review recent work in this vein and discuss a number of open questions for future study., 27 pages, including 7 figures. Submitted to Geoscience Letters as invited review from 2019 assembly of the Asia-Oceania Geosciences Society

Written by: Claire E. Newman (Aeolis Res.) Tanguy Bertrand (NASA Ames) Joseph Battalio (Yale Univ.) Mackenzie Day (UCLA) Manuel de la Torre Juarez (JPL) Meredith K. Elrod (NASA GSFC) Francesca Esposito (INAF-OAC) Lori Fenton (SETI Inst.) Claus Gebhardt (UAEU) Steven J. Greybush (Penn. State) Scott D. Guzewich (NASA GSFC) Henrik Kahanpaa (FMI) Melinda Kahre (NASA Ames) Ozgur Karatekin (Royal Obs. Belg.) Brian Jackson (Boise State Univ.) Mathieu Lapotre (Stanford Univ.) Christopher Lee (Aeolis Res.) Stephen R. Lewis (Open Univ.) Ralph D. Lorenz (APL) German Martinez Martinez (LPI) Javier Martin-Torres (Aberdeen U.) Michael A. Mischna (JPL) Luca Montabone (SSI) Lynn Neakrase (New Mexico State) Alexey Pankine (SSI) Jorge Pla-Garcia (CAB/SwRI/SSI) Peter L. Read (Univ. of Oxford) Isaac B. Smith (PSI/York Univ.) Michael D. Smith (NASA GSFC) Alejandro Soto (SwRI) Aymeric Spiga (Sorbonne Univ.) Christy Swann (NRL SSC) Leslie Tamppari (JPL) Orkun Temel (Royal Obs. Belgium) Daniel Viudez Moreiras (CAB) Danika Wellington (Ariz. State) Paulina Wolkenberg (INAF) Gerhard Wurm (Duisburg-Essen) Maria-Paz Zorzano (CAB)

Mars Polar Science is an integrated, compelling system that serves as a nearby analogue to numerous other planets, supports human exploration, and habitability. Mars possesses the closest and most easily accessible layered ice deposits outside of Earth, and accessing those layers to read the climate record would be a triumph for planetary science.

The regions of Saturn’s cloud-covered atmosphere polewards of 60o latitude are dominated in each hemisphere near the cloud tops by an intense, cyclonic polar vortex surrounded by a strong, high latitude eastward zonal jet. In the north, this high latitude jet takes the form of a remarkably regular zonal wavenumber m=6 hexagonal pattern that has been present at least since the Voyager spacecraft encounters with Saturn in 1980-81, and probably much longer. The origin of this feature, and the absence of a similar feature in the south, has remained poorly understood since its discovery. In this work, we present some new analyses of horizontal wind measurements at Saturn’s cloud tops polewards of 60 degrees in both the northern and southern hemispheres, previously published by Antuñano et al. (2015) using images from the Cassini mission, in which we compute kinetic energy spectra and the transfer rates of kinetic energy (KE) and enstrophy between different scales. 2D KE spectra are consistent with a zonostrophic regime, with a steep (~n-5) spectrum for the mean zonal flow (n is the total wavenumber) and a shallower Kolmogorov-like KE spectrum (~n-5/3) for the residual (eddy) flow, much as previously found for Jupiter’s atmosphere (Galperin et al. 2014; Young & Read 2017). Three different methods are used to compute the energy and enstrophy transfers, (a) as latitude-dependent zonal spectral fluxes, (b) as latitude-dependent structure functions and (c) as spatially filtered energy fluxes. The results of all three methods are largely in agreement in indicating a direct (forward) enstrophy cascade across most scales, averaged across the whole domain, an inverse kinetic energy cascade to large scales and a weak direct KE cascade at the smallest scales. The pattern of transfers has a more complex dependence on latitude, however. But it is clear that the m=6 North Polar Hexagon (NPH) wave was transferring KE into its zonal jet at 78o N (planetographic) at a rate of ∏E ≈ 1.8 x 10-4 W kg-1 at the time the Cassini images were acquired. This implies that the NPH was not maintained by a barotropic instability at this time, but may have been driven via a baroclinic instability or possibly from deep convection. Further implications of these results will be discussed. ReferencesAntuñano, A., T. del Río-Gaztelurrutia, A. Sánchez-Lavega, and R. Hueso (2015), Dynamics of Saturn’s polar regions, J. Geophys. Res. Planets, 120, 155–176, doi:10.1002/2014JE004709.Galperin, B., R. M.B. Young, S. Sukoriansky, N. Dikovskaya, P. L. Read, A. J. Lancaster & D. Armstrong (2014) Cassini observations reveal a regime of zonostrophic macroturbulence on Jupiter, Icarus, 229, 295–320.doi: 10.1016/j.icarus.2013.08.030Young, R. M. B. & Read, P. L. (2017) Forward and inverse kinetic energy cascades in Jupiter’s turbulent weather layer, Nature Phys., 13, 1135-1140. Doi:10.1038/NPHYS4227

Super-rotation is a phenomenon in atmospheric dynamics where the specific axial angular momentum of the wind (at some location) in an atmosphere exceeds that of the underlying planet at the equator. Hide's theorem states that in order for an atmosphere to super-rotate, non-axisymmetric disturbances (eddies) are required to induce transport of angular momentum up its local gradient. This raises a question as to the origin and nature of the disturbances that operate in super-rotating atmospheres to induce the required angular momentum transport. The primary technique employed to investigate this question has involved numerically modelling super-rotating atmospheres, and diagnosing the processes that give rise to super-rotation in the simulations. These modelling efforts can be separated into one of two approaches. The first approach utilises 'realistic', tailor-made models of Solar System atmospheres where super-rotation is present (e.g., Venus and Titan) to investigate the specific processes responsible for generating super-rotation on each planet. The second approach takes simple, 'Earth-like' models, typically dry dynamical cores with radiative transfer represented using a Newtonian cooling approach, and explores the effect of varying a single (or occasionally multiple) planetary parameters (e.g., the planetary radius or rotation rate) on the atmospheric dynamics. Notably, studies of this flavour have shown that super-rotation may emerge 'spontaneously' on planets with slow rotation rate or small radius (relative to the Earth's; Venus and Titan have these characteristics). However, the strength of super-rotation obtained in simulations of this type is far weaker than that observed in Venus' or Titan's atmospheres, or in tailored numerical models of either planet. In this work, our aim is to bridge the gap between these two modelling approaches. We will present results from a suite of simulations using an idealised general circulation model with a semi-grey representation of radiative transfer. Our experiments explore the effects of varying planetary size and rotation rate, atmospheric mass, and atmospheric absorption of shortwave radiation on the acceleration of super-rotation. A novel aspect of this work is that we vary multiple planetary properties away from their Earth-like 'defaults' in conjunction. This allows us to investigate how properties characteristic of the atmospheres of planets such as Venus and Titan combine to yield the strong super-rotation observed in their atmospheres (and realistic numerical models). We are also able to illustrate how features such as increased atmospheric mass and absorption of shortwave radiation modify the weakly super-rotating state obtained in simple, Earth-like models towards one more characteristic of Titan or Venus.

Super-rotation is a phenomenon in atmospheric dynamics where the axial angular momentum of an atmosphere in some way exceeds that of the underlying planet. In this presentation, we will discuss the dependency of both globally-integrated, and local metrics of super-rotation on planetary rotation rate, revealed through analysis of idealised General Circulation Model experiments. The model used here is based on the Held-Suarez benchmark for a dry, 'Earth-like' atmosphere, and results from both axisymmetric and three-dimensional experiments will be presented. Previous work has shown that the three-dimensional configuration used here will transition to a state of equatorial super-rotation if the rotation rate is reduced sufficiently from the Earth's. This motivates the question: How does super-rotation strength depend on rotation rate? We will use the term 'global super-rotation' to refer to an atmosphere with excess of globally-integrated axial angular momentum relative to that achieved by solid body co-rotation with the underlying planet, and 'local super-rotation' to refer to the existence of some region within the atmosphere where axial angular momentum exceeds that of the underlying planet at the equator. In an inviscid, axisymmetric atmosphere, the axial component of specific angular momentum is materially conserved. Consequently, in such a system local super-rotation is forbidden, although global super-rotation may still be achieved if a meridional circulation is able to transport fluid equilibrated with the equatorial surface poleward. If the restriction of axisymmetry is lifted, then local super-rotation may exist if non-axisymmetric disturbances that act to transport angular momentum up-gradient are present. The atmospheres of Venus, the Earth, Mars, and Titan may be considered to be globally super-rotating, however only Venus and Titan exhibit permanent local super-rotation at the equator. The results from axisymmetric experiments reveal that at high rotation rate (e.g., greater than 1/4 of the Earth's), the degree of global super-rotation scales inversely with the square of the rotation rate. In the low rotation rate limit, the degree of global super-rotation saturates, and becomes independent of rotation rate. We will show that the high, and low rotation rate dependencies can be predicted by a single analytic scaling for global super-rotation. Our three-dimensional experiments exhibit the same scaling behaviour for global super-rotation as observed in the axisymmetric experiments. The degree of global super-rotation achieved by the three-dimensional experiments is less than that of the axisymmetric experiments at high rotation rates, and greater at lower rotation rates, but in both limits the deviation from the axisymmetric 'base circulation' is small. In the low-rotation rate limit, local super-rotation is accelerated at the equator, which is consistent with the three-dimensional experiments obtaining a higher degree of global super-rotation than their axisymmetric counterparts. Estimates for global super-rotation strength on the Earth and Mars agree closely with the results of our three-dimensional numerical experiments, but Venus and Titan achieve substantially stronger global, and local super-rotation than found here. It appears that low rotation rate alone cannot induce substantial excess global super-rotation, relative to the axisymmetric base circulation we identify.

We present a series of experimental investigations in which a differentially-heated annulus was used to investigate the effects of topography on rotating, stratified flows. In particular, we investigate blocking effects via azimuthally varying differential-heating and compare them to previous experiments utilising partial mechanical barriers. The thermal topography used consisted of a flat patch of heating elements covering a small azimuthal extent of the base, forming an equivalent of a partial barrier, to study the difference between blocked and unblocked flow. These azimuthally-varying heating experiments produced results with many similarities to our previous experiments with a mechanical barrier, despite the lack of a physical obstacle or formation of bottom-trapped waves. In particular, a unique flow structure was found when the drifting flow and the topography interacted in the form of an “interference” regime at low Taylor number, but forming an erratic “irregular” regime at higher Taylor number. This suggests that blocking may be induced by either or both of a thermal or mechanical inhomogeneity. Evidence of coherent/persistent resonant wave triads was noted in both kinds of experiment, though the component wavenumbers of the wave-triads and their impact on the flow was found to depend on the topography in question.

We present a series of experimental investigations in which a differentially-heated annulus was used to investigate the effects of topography on rotating, stratified flows with similarities to the Earth’s atmospheric or oceanic circulation. In particular, we compare and investigate blocking effects via partial mechanical barriers to previous experiments by the authors utilising azimuthally-periodic topography. The mechanical obstacle used was an isolated ridge, forming a partial barrier, employed to study the difference between partially blocked and fully unblocked flow. The topography was found to lead to the formation of bottom-trapped waves, as well as impacting the circulation at a level much higher than the top of the ridge. This produced a unique flow structure when the drifting flow and the topography interacted in the form of an “interference” regime at low Taylor number, but forming an erratic “irregular” regime at higher Taylor number. The results also showed evidence of resonant wave-triads, similar to those noted with periodic wavenumber-3 topography by Marshall and Read (Geophys. Astrophys. Fluid Dyn., 2015, 109), though the component wavenumbers of the wave-triads and their impact on the flow were found to depend on the topography in question. With periodic topography, wave-triads were found to occur between both the baroclinic and barotropic components of the zonal wavenumber-3 mode and the wavenumber-6 baroclinic component, whereas with the partial barrier two nonlinear resonant wave-triads were noted, each sharing a common wavenumber-1 mode.

Jupiter’s turbulent weather layer contains phenomena of many different sizes, from local storms up to the Great Red Spot and banded jets. The global circulation is driven by complex interactions with (as yet uncertain) small scale processes. We have calculated structure functions and kinetic energy spectral fluxes from Cassini observations over a wide range of length scales in Jupiter’s atmosphere. We found evidence for an inverse cascade of kinetic energy from length scales comparable with the first baroclinic Rossby deformation radius to the global jet scale, but also a forward cascade of kinetic energy from the deformation radius to smaller scales. The latter disagrees with the traditional picture of Jupiter’s atmospheric dynamics, but has some similarities with mesoscale phenomena in the Earth’s atmosphere and oceans. We conclude that the inverse cascade driving Jupiter’s jets may have a dominant energy source at scales close to the deformation radius, such as baroclinic instability.

The atmospheric circulation in Venus is well known to exhibit strong super-rotation. However, the atmospheric mechanisms responsible for the formation of this super-rotation are still not fully understood. In this work, we developed a new Venus general circulation model to study the most likely mechanisms driving the atmosphere to the current observed circulation. Our model includes a new radiative transfer, convection and suitably adapted boundary layer schemes and a dynamical core that takes into account the dependence of the heat capacity at constant pressure with temperature. The new Venus model is able to simulate a super-rotation phenomenon in the cloud region quantitatively similar to the one observed. The mechanisms maintaining the strong winds in the cloud region were found in the model results to be a combination of zonal mean circulation, thermal tides and transient waves. In this process, the semi-diurnal tide excited in the upper clouds has a key contribution in transporting axial angular momentum mainly from the upper atmosphere towards the cloud region. The magnitude of the super-rotation in the cloud region is sensitive to various radiative parameters such as the amount of solar radiative energy absorbed by the surface, which controls the static stability near the surface. In this work, we also discuss the main difficulties in representing the flow below the cloud base in Venus atmospheric models. Our new radiative scheme is more suitable for 3D Venus climate models than those used in previous work due to its easy adaptability to different atmospheric conditions. This flexibility of the model was crucial to explore the uncertainties in the lower atmospheric conditions and may also be used in the future to explore, for example, dynamical-radiative-microphysical feedbacks., Accepted for publication in Planet. Space Sci

A Martian semiannual oscillation (SAO), similar to that in the Earth's tropical stratosphere, is evident in the Mars Analysis Correction Data Assimilation reanalysis dataset (MACDA) version 1.0, not only in the tropics, but also extending to higher latitudes. Unlike on Earth, the Martian SAO is found not always to reverse its zonal wind direction, but only manifests itself as a deceleration of the dominant wind at certain pressure levels and latitudes. Singular System Analysis (SSA) is further applied on the zonal-mean zonal wind in different latitude bands to reveal the characteristics of SAO phenomena at different latitudes. The second pair of principal components (PCs) is usually dominated by a SAO signal, though the SAO signal can be strong enough to manifest itself also in the first pair of PCs. An analysis of terms in the Transformed Eulerian Mean equation (TEM) is applied in the tropics to further elucidate the forcing processes driving the tendency of the zonal-mean zonal wind. The zonal-mean meridional advection is found to correlate strongly with the observed oscillations of zonal-mean zonal wind, and supplies the majority of the westward (retrograde) forcing in the SAO cycle. The forcing due to various non-zonal waves supplies forcing to the zonal-mean zonal wind that is nearly the opposite of the forcing due to meridional advection above ~3 Pa altitude, but it also partly supports the SAO between 40 Pa and 3 Pa. Some distinctive features occurring during the period of the Mars year (MY) 25 global-scale dust storm (GDS) are also notable in our diagnostic results with substantially stronger values of eastward and westward momentum in the second half of MY 25 and stronger forcing due to vertical advection, transient waves and thermal tides., Accepted for publication in Icarus. 20 pages, 5 figures

International audience

We review the recent progress in understanding the jet structures on exoplanets as well as on and inside the Sun. The emphasis is on the more robust aspects of observation and numerical modeling that relate directly to jets. For the exoplanets, the primary focus is on hot-Jupiters since many more observations are available for them presently than other types of exoplanets. Because not much is known about the morphology and strength of the jets on exoplanets, there is currently not much agreement. In contrast, the picture is very different for the Sun. In fact, the jet structure of the Sun is arguably one of the best known jet structures of all the planets and stars, due to the fact that Sun's disk is resolved and its interior can be probed with helioseismology. A discussion of several critical issues pertaining to the modeling of jets on exoplanets and the Sun is presented, along with a brief outlook on the subject.

The energetics of possible global atmospheric circulation patterns in an Earth-like atmosphere are explored using a simplified GCM based on the University of Hamburg's Portable University Model for the Atmosphere. Results from a series of simulations, obtained by varying planetary rotation rate {\Omega} with an imposed equator-to-pole temperature difference, were analysed to determine the heat transport and other contributions to the energy budget for the time-averaged, equilibrated flow. These show clear trends with {\Omega}, with the most intense Lorenz energy cycle for an Earth-sized planet occurring with a rotation rate around half that of the present day Earth. KE and APE spectra, E_K(n) and E_A(n) (where n is total spherical wavenumber), also show clear trends with \Omega, with n^{-3} enstrophy-dominated spectra around \Omega* = \Omega/\Omega_E = 1, where \Omega_E is the rotation rate of the Earth) and steeper (\sim n^{-5}) slopes in the zonal mean flow with little evidence for the n^{-5/3} spectrum anticipated for an inverse KE cascade. Instead, both KE and APE spectra become almost flat at scales larger than the internal Rossby radius, L_d, and exhibit near-equipartition at high wavenumbers. At \Omega* << 1, the spectrum becomes dominated by KE with E_K(n) \sim 2-3 E_A(n) at most wavenumbers and a slope \sim n^{-5/3} across most of the spectrum. Spectral flux calculations show that enstrophy and APE are almost always cascaded downscale, regardless of {\Omega}. KE cascades are more complicated, however, with downscale transfers across almost all wavenumbers, dominated by horizontally divergent modes, for \Omega* \lesssim 1/4. At higher rotation rates, transfers of KE become increasingly dominated by rotational components with strong upscale transfers (dominated by eddy-zonal flow interactions) for scales larger than L_d and weaker downscale transfers for scales smaller than L_d., Comment: 16 pages (19 as published), 9 figures

Scatter plots of Ertel potential vorticity, Q, versus Bernoulli streamfunction, B, on potential temperature surfaces, θ, are investigated for Mars using the global Mars Analysis Correction Data Assimilation (MACDA) reanalysis, which spans Mars Year (MY) 24.39 to 27.24. In midlatitudes, Mars exhibits monotonic, function-like Q(B) correlations on θ surfaces similar to those observed for Earth. We quantify this with linear regressions of Q versus B over the vertical range θ=400 to 900 K (∼30 to 60 km). In autumn, winter and spring, in both hemispheres, the non-dimensionalized correlation generally lies between zero and unity and gradually decreases with height, whereas in northern summer, it swings negative. These characteristics match Earth's lower mesosphere (θ= 2000 to 3000 K; z≈ 48 to 62 km) during the same seasons. The exception is southern summer, when the correlation on Mars nearly vanishes. In time series, the transition into and out of northern summer is sinuous and centred just after solar longitude Ls = 90°, whereas in southern summer it is abrupt and spans ΔLs≈120°, which is one third of a Mars year. A striking feature seen on Mars but not on Earth is a large range of Q over the narrow domain of B poleward of each winter polar jet, particularly in the north, which is consistent with the known annular structure of the Martian polar vortex. Froude number calculations suggest the existence of a planetary-scale hydraulic jump associated with the winter polar jet.

A two-dimensional regime diagram for geostrophic turbulence in the ocean is constructed by plotting observation-based estimates of the non-dimensional eddy length-scale against a nonlinearity parameter equal to the ratio of the root-mean-square eddy velocity and baroclinic Rossby phase speed. Two estimates of the eddy length-scale are compared: the equivalent eddy radius inferred from the area enclosed by contours of sea-surface height, and the ‘unsuppressed’ mixing length, based on an estimate of the eddy diffusivity with mean flow effects removed. For weak nonlinearity, as found in the Tropics, the mixing length mostly corresponds to the stability threshold for baroclinic instability whereas the eddy radius corresponds to the Rhines scale; it is suggested that this mismatch is indicative of the inverse energy cascade that occurs at low latitudes in the ocean and the zonal elongation of eddies. At larger values of nonlinearity, as found at mid- and high latitudes, the eddy length-scales are much shorter than the stability threshold, within a factor of 2.5 of the Rossby deformation radius.

Saturn's polar stratosphere exhibits the seasonal growth and dissipation of broad, warm, vortices poleward of $\sim75^\circ$ latitude, which are strongest in the summer and absent in winter. The longevity of the exploration of the Saturn system by Cassini allows the use of infrared spectroscopy to trace the formation of the North Polar Stratospheric Vortex (NPSV), a region of enhanced temperatures and elevated hydrocarbon abundances at millibar pressures. We constrain the timescales of stratospheric vortex formation and dissipation in both hemispheres. Although the NPSV formed during late northern spring, by the end of Cassini's reconnaissance (shortly after northern summer solstice), it still did not display the contrasts in temperature and composition that were evident at the south pole during southern summer. The newly-formed NPSV was bounded by a strengthening stratospheric thermal gradient near $78^\circ$N. The emergent boundary was hexagonal, suggesting that the Rossby wave responsible for Saturn's long-lived polar hexagon - which was previously expected to be trapped in the troposphere - can influence the stratospheric temperatures some 300 km above Saturn's clouds., 51 pages, 12 figures, published in Nature Communications

The regimes of possible global atmospheric circulation patterns in an Earth-like atmosphere are explored using a simplified GCM based on the University of Hamburg's Portable University Model for the Atmosphere with simplified (linear) boundary layer friction, a Newtonian cooling scheme and dry convective adjustment. A series of controlled experiments are conducted by varying planetary rotation rate and imposed equator-to-pole temperature difference. These defining parameters are cast into dimensionless forms to establish a parameter space, in which different circulation regimes are mapped and classified. Clear trends are found when varying planetary rotation rate and frictional and thermal relaxation timescales. The sequence of circulation regimes as a function of planetary rotation rate strongly resembles that obtained in laboratory experiments on rotating, stratified flows, especially if a topographic $\beta$-effect is included in those experiments to emulate the planetary vorticity gradients induced by the spherical curvature of the planet. A regular baroclinic wave regime is also obtained at intermediate values of thermal Rossby number and its characteristics and dominant zonal wavenumber depend strongly on the strength of radiative and frictional damping. These regular waves exhibit some strong similarities to baroclinic storms observed on Mars under some conditions. Multiple jets are found at the highest rotation rates, when the Rossby deformation radius and other eddy-related length scales are much smaller than the radius of the planet. These exhibit some similarity to the multiple zonal jets observed on gas giant planets. Jets form on a scale comparable to the most energetic eddies and the Rhines scale poleward of the supercritical latitude. The balance of heat transport varies strongly with {\Omega}* between eddies and zonally symmetric flows, becoming weak with fast rotation., Comment: 18 pages (21 pages as published), 11 figures

In this study, we consider how the wave number selection in spherical Couette flow, in the transition to azimuthal waves after the first instability, occurs in the presence of noise. The outer sphere was held stationary, while the inner sphere rotational speed was increased linearly from a subcritical flow to a supercritical one. In a supercritical flow, one of two possible flow states, each with different azimuthal wave numbers, can appear depending upon the initial and final Reynolds numbers and the acceleration value. Noise perturbations were added by introducing small disturbances into the rotational speed signal. With an increasing noise amplitude, a change in the dominant wave number from m to m +/- 1was found to occur at the same initial and final Reynolds numbers and acceleration values. The flow velocity measurements were conducted by using laser Doppler anemometry. Using these results, the role of noise as well as the behaviour of the amplitudes of the competing modes in their stages of damping and growth were determined. Published by AIP Publishing.

The superrotation of the atmospheres of Venus and Titan has puzzled dynamicists for many years and seems to put these planets in a very different dynamical regime from most other planets. In this review, we consider how to define superrotation objectively and explore the constraints that determine its occurrence. Atmospheric superrotation also occurs elsewhere in the Solar System and beyond, and we compare Venus and Titan with Earth and other planets for which wind estimates are available. The extreme superrotation on Venus and Titan poses some difficult challenges for numerical models of atmospheric circulation, much more difficult than for more rapidly rotating planets such as Earth or Mars. We consider mechanisms for generating and maintaining a superrotating state, all of which involve a global meridional overturning circulation. The role of nonaxisymmetric eddies is crucial, however, but the detailed mechanisms may differ between Venus, Titan, and other planets.

In this paper, we show that the behavior of weakly nonlinear waves in a 2-layer model of baroclinic instability on a b-plane with varying viscosity is determined by a single, degenerate codimension three bifurcation. In the process, we show how previous studies, using the method of multiple scales to derive evolution equations for the slowly varying amplitude of the growing wave, arise as special limits of the general evolution description.

This paper documents an experimental investigation in which a differentially-heated rotating annulus experiment was used to investigate the effects of topography on fluid flow under conditions similar to the atmospheric and oceanic circulation on Earth and other planets. In particular, the relationship between the effects of topographic resonance and the existence and mechanism for generation of low-frequency variability (LFV) were studied, motivated by outstanding questions in works such as Jin and Ghil (J. Atmos. Sci., 1990, 47) and Read and Risch (Geophys. Astrophys. Fluid Dyn., 2011, 105). Whilst employing sinusoidal wavenumber-3 topography a new regime was encountered within a region of stationary wavenumber-3 structural vacillation. Denoted as the “stationary-transition” regime, it featured periodic oscillations between a dominant stationary wavenumber-3 flow and axisymmetric or chaotic flow. Further investigation found that the “stationary-transition” regime appeared to be a near-resonant region wh...

A laboratory investigation is presented with the aim of studying systematically the occurrence and characteristics of low-frequency variability of flows resulting from the interaction of a baroclinic flow with periodic bottom topography. Low-frequency variability within the baroclinic wave regime occurred in two distinct forms in separate regions of parameter space. One corresponded to the transition region between the baroclinic travelling and stationary wave regimes. It involved primarily an interaction between the drifting baroclinic waves and stationary components of the topographically forced wave. The resulting flow had characteristics similar to amplitude vacillation and had a time-scale of 30–60 annulus revolutions (days), which also corresponded to the wave drift period. A new regime of low-frequency amplitude vacillation was discovered in the transition region with the axisymmetric flow regime. As the complexity of the flow increased the period of the vacillation cycles grew to ∼100–180 “days”. This slower vacillation seemed to involve a cyclic enabling and disabling of nonlinear interactions between the forced stationary wave and the growing and azimuthally drifting wave, which in turn was linked to a decrease in mean flow shear. Subsequent chains of wave-wave interactions characterised the complex but robust oscillation phenomenon. The resulting behaviour has several features in common with some recent models of intraseasonal oscillations in the mid-latitude troposphere and with sudden stratospheric warmings.

The non-axisymmetric structure of an unstable Stewartson shear layer generated via a differential rotation between flush disks and a cylindrical enclosure is investigated numerically using both three-dimensional direct numerical simulation and a quasi-two-dimensional model. Previous literature has only considered the depth-independent quasi-two-dimensional model due to its low computational cost. The three-dimensional model implemented here highlights the supercritical instability responsible for the polygonal deformation of the shear layer in the linear and nonlinear growth regimes and reveals that linear stability analysis is capable of accurately determining the preferred azimuthal wavenumber for flow conditions near the onset of instability. This agreement is lost for sufficiently forced flows where nonlinear effects encourage the coalescence of vortices towards lower-wavenumber structures. Time-dependent flows are found for large Reynolds numbers defined based on the Stewartson layer thickness and azimuthal velocity differential. However, this temporal behaviour is not solely characterized by Reynolds number but is rather a function of both the Rossby and Ekman numbers. At high Ekman and Rossby numbers, unsteady flow emerges through a small-scale azimuthal destabilization of the axial jets within the Stewartson layers; at low Ekman numbers, unsteady flow emerges through a modulation in the strength of one of the axial vortices rolled up by non-axisymmetric instability of the Stewartson layer.

We present a new radiation scheme for the Oxford Planetary Unified Model System for Venus, suitable for the solar and thermal bands. This new and fast radiative parameterization uses a different approach in the two main radiative wavelength bands: solar radiation ( 0.1 – 5.5 μ m ) and thermal radiation ( 1.7 – 260 μ m ) . The solar radiation calculation is based on the δ -Eddington approximation (two-stream-type) with an adding layer method. For the thermal radiation case, a code based on an absorptivity/emissivity formulation is used. The new radiative transfer formulation implemented is intended to be computationally light, to allow its incorporation in 3D global circulation models, but still allowing for the calculation of the effect of atmospheric conditions on radiative fluxes. This will allow us to investigate the dynamical-radiative-microphysical feedbacks. The model flexibility can be also used to explore the uncertainties in the Venus atmosphere such as the optical properties in the deep atmosphere or cloud amount. The results of radiative cooling and heating rates and the global-mean radiative–convective equilibrium temperature profiles for different atmospheric conditions are presented and discussed. This new scheme works in an atmospheric column and can be easily implemented in 3D Venus global circulation models.

We present a numerical study of axisymmetric flow in a rotating annulus in which local thermal forcing, via a heated annular ring on the outside of the base and a cooled circular disk in the centre of the top surface, drives convection. This new configuration is a variant of the classical thermally-driven annulus, where uniform heating and cooling are applied through the outer and inner sidewalls respectively. The annulus provides an analogue to a planetary circulation and the new configuration, with its more relaxed vertical thermal boundary conditions, is expected to better emulate vigorous convection in the tropics and polar regions as well as baroclinic instability in the mid-latitude baroclinic zone. Using the Met Office/Oxford Rotating Annulus Laboratory (MORALS) code, we have investigated a series of equilibrated, two dimensional axisymmetric flows across a large region of parameter space. These are characterized in terms of their velocity and temperature fields. When rotation is applied several distinct flow regimes may be identified for different rotation rates and strengths of differential heating. These regimes are defined as a function of the ratio of the horizontal Ekman layer thickness to the non-rotating thermal boundary layer thickness and are found to be similar to those identified in previous annulus experiments. Convection without rotation is also considered and the scaling of the heat transport with Rayleigh number is calculated. This is then compared with existing work on the classical annulus as well as horizontal and Rayleigh-Bénard convection. As with previous studies on both rotating and non-rotating convection the system’s behaviour is found to be aspect ratio dependent. This dependence is seen in the scaling of the non-rotating Nusselt number and in transitions between regimes in the rotating case although further investigation is required to fully explain these observations.

We review our current knowledge of the atmospheric dynamics of Venus prior to the Akatsuki mission, in the altitude range from the surface to approximately the cloud tops located at about 100 km altitude. The three-dimensional structure of the wind field in this region has been determined with a variety of techniques over a broad range of spatial and temporal scales (from the mesoscale to planetary, from days to years, in daytime and nighttime), spanning a period of about 50 years (from the 1960s to the present). The global panorama is that the mean atmospheric motions are essentially zonal, dominated by the so-called super-rotation (an atmospheric rotation that is 60 to 80 times faster than that of the planetary body). The zonal winds blow westward (in the same direction as the planet rotation) with a nearly constant speed of $\sim 100~\mathrm{m}\, \mathrm{s}^{-1}$ at the cloud tops (65–70 km altitude) from latitude 50°N to 50°S, then decreasing their speeds monotonically from these latitudes toward the poles. Vertically, the zonal winds decrease with decreasing altitude towards velocities $\sim 1\text{--}3~\mathrm{m}\,\mathrm{s}^{-1}$ in a layer of thickness $\sim 10~\text{km}$ close to the surface. Meridional motions with peak speeds of $\sim 15~\mathrm{m}\,\mathrm{s}^{-1}$ occur within the upper cloud at 65 km altitude and are related to a Hadley cell circulation and to the solar thermal tide. Vertical motions with speeds $\sim1\text{--}3~\mathrm{m}\, \mathrm{s}^{-1}$ occur in the statically unstable layer between altitudes of $\sim 50 \text{--} 55~\text{km}$ . All these motions are permanent with speed variations of the order of $\sim10\%$ . Various types of wave, from mesoscale gravity waves to Rossby-Kelvin planetary scale waves, have been detected at and above cloud heights, and are considered to be candidates as agents for carrying momentum that drives the super-rotation, although numerical models do not fully reproduce all the observed features. Momentum transport by atmospheric waves and the solar tide is thought to be an indispensable component of the general circulation of the Venus atmosphere. Another conspicuous feature of the atmospheric circulation is the presence of polar vortices. These are present in both hemispheres and are regions of warmer and lower clouds, seen prominently at infrared wavelengths, showing a highly variable morphology and motions. The vortices spin with a period of 2–3 days. The South polar vortex rotates around a geographical point which is itself displaced from the true pole of rotation by $\sim 3$ degrees. The polar vortex is surrounded and constrained by the cold collar, an infrared-dark region of lower temperatures. We still lack detailed models of the mechanisms underlying the dynamics of these features and how they couple (or not) to the super-rotation. The nature of the super-rotation relates to the angular momentum stored in the atmosphere and how it is transported between the tropics and higher latitudes, and between the deep atmosphere and upper levels. The role of eddy processes is crucial, but likely involves the complex interaction of a variety of different types of eddy, either forced directly by radiative heating and mechanical interactions with the surface or through various forms of instability. Numerical models have achieved some significant recent success in capturing some aspects of the observed super-rotation, consistent with the scenario discussed by Gierasch (J. Atmos. Sci. 32:1038–1044, 1975) and Rossow and Williams (J. Atmos. Sci. 36:377–389, 1979), but many uncertainties remain, especially in the deep atmosphere. The theoretical framework developed to explain the circulation in Venus’s atmosphere is reviewed, as well as the numerical models that have been built to elucidate the super-rotation mechanism. These tools are used to analyze the respective roles of the different waves in the processes driving the observed motions. Their limitations and suggested directions for improvements are discussed.

Mars is a relatively small, desert planet with only a thin atmosphere at present. These basic properties are all directly reflected in its global circulation, which is fundamentally similar to that on Earth in many ways, but very different in other ways. The investigation of the global circulation of the Martian atmosphere has a surprisingly long history, as a result of the basic terrestrial similarities and the very early spacecraft missions to Mars. The virtual explosion in spacecraft observations during the last two decades has made the global atmospheric circulation of Mars by far the best characterized and understood atmospheric circulation in the Solar System, other than that of Earth. By comparison with Earth our knowledge of the global atmospheric circulation of Mars remains quite limited: there are extremely few direct wind observations, very few surface pressure observations, and the local time coverage of the global observations is very limited. Our knowledge of the planetary boundary layer (which can be extremely deep by comparison to Earth) on Mars is limited at present. It is only with the aid of extensive numerical modeling that our fundamental understanding of the global circulation of the Martian atmosphere has been fleshed out. Data assimilation efforts for Mars are still at a relatively early stage but these certainly have great potential for the future – if observations of the global atmosphere continue to be obtained.

The Martian planetary boundary layer (PBL) consists of the layers of the atmosphere closest to the surface, within which interactions between the atmosphere and the surface itself are dominant. In general, this represents the lowest 1-10 km of the atmosphere, within which surface-driven intense convection may take place, with convective plumes and vortices rising to heights in excess of 5-10 km during the day (Thomas and Gierasch, 1985; Haberle et al., 1993b; Larsen et al., 2002; Balme and Greeley, 2006; Hinson et al., 2008). At night, convection is inhibited and radiative cooling produces a stably-stratified layer at the surface, and the PBL reduces to a shallow layer forced by mechanical turbulence at the bottom of the stable layer. It is therefore a highly dynamic and variable region of the atmosphere at virtually all locations on Mars, with additional variability induced by interactions with local surface topography.

We present an experimental study of flows in a cylindrical rotating annulus convectively forced by local heating in an annular ring at the bottom near the external wall and via a cooled circular disk near the axis at the top surface of the annulus. This new configuration is distinct from the classical thermally driven annulus analogue of the atmosphere circulation, in which thermal forcing is applied uniformly on the sidewalls, but with a similar aim to investigate the baroclinic instability of a rotating, stratified flow subjected to zonally symmetric forcing. Two vertically and horizontally displaced heat sources/sinks are arranged, so that in the absence of background rotation, statically unstable Rayleigh–Benard convection would be induced above the source and beneath the sink, thereby relaxing strong constraints placed on background temperature gradients in previous experimental configurations based on the conventional rotating annulus. This better emulates local vigorous convection in the tropics and polar regions of the atmosphere while also allowing stably-stratified baroclinic motion in the central zone of the annulus, as in mid-latitude regions in the Earth’s atmosphere. Regimes of flow are identified, depending mainly upon control parameters that in turn depend on rotation rate and the strength of differential heating. Several regimes exhibit baroclinically unstable flows which are qualitatively similar to those previously observed in the classical thermally driven annulus. However, in contrast to the classical configuration, they typically exhibit more spatio-temporal complexity. Thus, several regimes of flow demonstrate the equilibrated co-existence of, and interaction between, free convection and baroclinic wave modes. These new features were not previously observed in the classical annulus and validate the new setup as a tool for exploring fundamental atmosphere-like dynamics in a more realistic framework. Thermal structure in the fluid is investigated and found to be qualitatively consistent with previous numerical results, with nearly isothermal conditions, respectively, above and below the heat source and sink, and stably-stratified, sloping isotherms in the near-adiabatic interior.

Masoud Mohseni, Peter Read, Hartmut Neven and colleagues at Google's Quantum AI Laboratory set out investment opportunities on the road to the ultimate quantum machines.

The Mars Analysis Correction Data Assimilation (MACDA) dataset version 1.0 contains the reanalysis of fundamental atmospheric and surface variables for the planet Mars covering a period of about three Martian years (a Martian year is about 1.88 terrestrial years). This has been produced by data assimilation of observations from NASA's Mars Global Surveyor (MGS) spacecraft during its science mapping phase (February 1999–August 2004). In particular, we have used retrieved thermal profiles and total dust optical depths from the Thermal Emission Spectrometer (TES) on board MGS. Data have been assimilated into a Mars global climate model (MGCM) using the Analysis Correction scheme developed at the UK Meteorological Office. The MGCM used is the UK spectral version of the Laboratoire de Météorologie Dynamique (LMD, Paris, France) MGCM. MACDA is a joint project of the University of Oxford and The Open University in the UK.

The thermally driven rotating annulus is a laboratory experiment important for the study of the dynamics of planetary atmospheres under controllable and reproducible conditions. We use the analysis correction method to assimilate laboratory data into an annulus model. We analyze the 2S and 3AV regular flow regimes between rotation rates of 0.75 and 0.875 rad s−1 and the 3SV chaotic flow regime between rotation rates of 2.2 and 3.1 rad s−1. Our assimilated observations are irregularly distributed, which is more meteorologically realistic than gridded observations as used in recent applications of data assimilation to laboratory measurements. We demonstrate that data assimilation can be used successfully and accurately in this context. We examine a number of specific assimilation scenarios: a wave-number transition between two regimes, information propagation from data-rich to data-poor regions, the response of the assimilation to a strong disturbance to the flow, and a vortex-shedding instability phenomenon at high rotation rate. At the highest rotation rates we calculated the barotropic E-vectors using unobserved variables such as temperature and the vertical structure of the velocity field that are only available via the assimilation. These showed that the mean flow is weakened by the action of eddies, going some way towards explaining why vortices are shed at the very highest rotation rates but not at lower rotation. Rossby-wave stability theory suggests that the underlying instability leading to vortex shedding may be baroclinic in character.

The relative merits of three contrasting empirical orthogonal decomposition methods in common use (namely, Proper Orthogonal Decomposition, Biorthogonal Decomposition and Multivariate Singular Systems Analysis) are considered as applied to baroclinic flow data. The regimes analysed are a steady, drifting wave, a modulated amplitude vacillating wave flow and a neighbouring multi-mode state which exhibits intermittency. The results are used to make a qualitative comparison of the methods in terms of convergence properties, variance capture and eigenfunction structure. The feasibility of using the resulting empirical orthogonal functions to transform partial differential equations to ordinary differential equations by Galerkin projection is mentioned. © 1997 Elsevier Science B.V.

A new general circulation model (GCM) of Venus is being developed at Oxford. Venus presents unique numerical and physical challenges because of its thick atmosphere, slow underlying solid body rotation, and super-rotating atmosphere. Preliminary results from a GCM with simplified physical parameterizations are discussed. The current model uses linearized cooling and friction schemes, and spans five decades of pressure (0-90 km). The model is able to demonstrate significant global super-rotation, and although not yet fully realistic, future plans include more detailed representation of the Venusian atmosphere, such as the planetary boundary layer (PBL) scheme. The use of the model is discussed in supporting and interpreting data from future missions to Venus. (c) 2005 COSPAR. Published by Elsevier Ltd. All rights reserved.

Data assimilation is a technique for the analysis of atmospheric observations which combines current information with prior knowledge from previous observations, summarized and forecast in time via the use of a numerical model. A sequential data assimilation scheme has been implemented with a full general circulation model (GCM) of the martian atmosphere for the first time, and has been adapted for the types of atmospheric data which might be expected in the near future, e.g. remote-sensed temperature profiles from a polar orbiter mission such as Mars Surveyor '96 and '98. Tests demonstrate the performance of the scheme using artificial data generated from independent model experiments. (C) 1997 COSPAR. Published by Elsevier Science Ltd.