The effects of vigorous mixing in a convective model of zonal flow on the ice giants

Previous studies have used models of three-dimensional (3D) Boussinesq convection in a rotating spherical shell to explain the zonal flows on the gas giants, Jupiter and Saturn. In this paper we demonstrate that this approach can also generate flow patterns similar to those observed on the ice giants, Uranus and Neptune. The equatorial jets of Uranus and Neptune are often assumed to result from baroclinic cloud layer processes and have been simulated with shallow layer models. Here we show that vigorous, 3D convection in a spherical shell can produce the retrograde (westward) equatorial flows that occur on the ice giants as well as the prograde (eastward) equatorial flows of the gas giants. In our models, the direction of the equatorial jet depends on the ratio of buoyancy to Coriolis forces in the system. In cases where Coriolis forces dominate buoyancy, cylindrical Reynolds stresses drive prograde equatorial jets. However, as buoyancy forces approach and exceed Coriolis forces, the cylindrical nature of the flow is lost and 3D mixing homogenizes the fluid's angular momentum; the equatorial jet reverses direction, while strong prograde jets form in the polar regions. Although the results suggest that conditions involving strong atmospheric mixing are responsible for generating the zonal flows on the ice giants, our present models require roughly 100 and 10 times the internal heat fluxes observed on Uranus and Neptune, respectively.     

A 2D study of the effects of the size of a solid core on the equatorial flow in giant planets

Zonal winds simulated in two-dimensional computer models of turbulent convection in the equatorial plane of giant planets have greater surface amplitudes for cases with smaller solid cores, and therefore larger buoyancy driving, all other properties being equal. This differential rotation in radius is maintained by the convergence of angular momentum flux, which occurs because of the convective flow that develops due to the effects of planetary rotation and density stratification. The superposition of the convective flow and the stronger zonal flow produces wave-like, instead of cellular convection.     

Compressible convection in the deep atmospheres of giant planets

Fast rotating giant planets such as Jupiter and Saturn possess alternate prograde and retrograde zonal winds which are stable over long periods of time. We consider a compressible model of convection in a spherical shell with rapid rotation, using the anelastic approximation, to explore the parameter range for which such zonal flows can be produced.We consider models with a large variation in density across the layer. Our models are based only on the molecular H/He region above the metallic hydrogen transition at about 2 Mbar, and we do not include the hydromagnetic effects which may be important if the electrical conductivity is significant. We find that the convective velocities are significantly higher in the low density regions of the shell, but the zonal flow is almost independent of the z-coordinate parallel to the rotation axis. We analyse how this behaviour is consistent with the Proudman-Taylor theorem.We find that deep prograde zonal flow near the equator is a very robust feature of our models. Prograde and retrograde jets alternating in latitude can occur inside the tangent cylinder in compressible as well as Boussinesq models, particularly at lower Prandtl numbers. However, the zonal jets inside the tangent cylinder are suppressed if a no-slip condition is imposed at the inner boundary. This suggests that deep high latitude jets may be suppressed if there is significant magnetic dissipation.Our compressible calculations include the viscous dissipation in the entropy equation, and we find this is comparable to, and in some cases exceeds, the total heat flux emerging from the surface. For numerical reasons, these simulations cannot reach the extremely low Ekman number found in giant planets, and they necessarily also have a much larger heat flux than planets. We therefore discuss how our results might scale down to give solutions with lower dissipation and lower heat flux.     

Interaction of moist convection with zonal jets on Jupiter and Saturn

Observations suggest that moist convection plays an important role in the large-scale dynamics of Jupiter's and Saturn's atmospheres. Here we use a reduced-gravity quasigeostrophic model, with a parameterization of moist convection that is based on observations, to study the interaction between moist convection and zonal jets on Jupiter and Saturn. Stable jets with approximately the same width and strength as observations are generated in the model. The observed zonal jets violate the barotropic stability criterion but the modeled jets do so only if the flow in the deep underlying layer is westward. The model results suggest that a length scale and a velocity scale associated with moist convection control the width and strength of the jets. The length scale and velocity scale offer a possible explanation of why the jets of Saturn are stronger and wider than those of Jupiter.     

Generation of equatorial jets by large-scale latent heating on the giant planets

Three-dimensional numerical simulations show that large-scale latent heating resulting from condensation of water vapor can produce multiple zonal jets similar to those on the gas giants (Jupiter and Saturn) and ice giants (Uranus and Neptune). For plausible water abundances (3-5 times solar on Jupiter/Saturn and 30 times solar on Uranus/Neptune), our simulations produce ∼20 zonal jets for Jupiter and Saturn and 3 zonal jets on Uranus and Neptune, similar to the number of jets observed on these planets. Moreover, these Jupiter/Saturn cases produce equatorial superrotation whereas the Uranus/Neptune cases produce equatorial subrotation, consistent with the observed equatorial-jet direction on these planets. Sensitivity tests show that water abundance, planetary rotation rate, and planetary radius are all controlling factors, with water playing the most important role; modest water abundances, large planetary radii, and fast rotation rates favor equatorial superrotation, whereas large water abundances favor equatorial subrotation regardless of the planetary radius and rotation rate. Given the larger radii, faster rotation rates, and probable lower water abundances of Jupiter and Saturn relative to Uranus and Neptune, our simulations therefore provide a possible mechanism for the existence of equatorial superrotation on Jupiter and Saturn and the lack of superrotation on Uranus and Neptune. Nevertheless, Saturn poses a possible difficulty, as our simulations were unable to explain the unusually high speed (∼) of that planet’s superrotating jet. The zonal jets in our simulations exhibit modest violations of the barotropic and Charney-Stern stability criteria. Overall, our simulations, while idealized, support the idea that latent heating plays an important role in generating the jets on the giant planets.     

Constraints on deep-seated zonal winds inside Jupiter and Saturn

The atmospheres of Jupiter and Saturn exhibit strong and stable zonal winds. How deep the winds penetrate unabated into each planet is unknown. Our investigation favors shallow winds. It consists of two parts. The first part makes use of an Ohmic constraint; Ohmic dissipation associated with the planet's magnetic field cannot exceed the planet's net luminosity. Application to Jupiter (J) and Saturn (S) shows that the observed zonal winds cannot penetrate below a depth at which the electrical conductivity is about six orders of magnitude smaller than its value at the molecular-metallic transition. Measured values of the electrical conductivity of molecular hydrogen yield radii of maximum penetration of 0.96RJ and 0.86RS, with uncertainties of a few percent of R. At these radii, the magnetic Reynolds number based on the zonal wind velocity and the scale height of the magnetic diffusivity is of order unity. These limits are insensitive to difficulties in modeling turbulent convection. They permit complete penetration along cylinders of the equatorial jets observed in the atmospheres of Jupiter and Saturn. The second part investigates how deep the observed zonal winds actually do penetrate. As it applies heuristic models of turbulent convection, its conclusions must be regarded as tentative. Truncation of the winds in the planet's convective envelope would involve breaking the Taylor-Proudman constraint on cylindrical flow. This would require a suitable nonpotential acceleration which none of the obvious candidates appears able to provide. Accelerations arising from entropy gradients, magnetic stresses, and Reynolds stresses appear to be much too weak. These considerations suggest that strong zonal winds are confined to shallow, stably stratified layers, with equatorial jets being the possible exception.     

Turbulent convection in rapidly rotating spherical shells: A model for equatorial and high latitude jets on Jupiter and Saturn

The origin of zonal jets on the jovian planets has long been a topic of scientific debate. In this paper we show that deep convection in a spherical shell can generate zonal flow comparable to that observed on Jupiter and Saturn, including a broad prograde equatorial jet and multiple alternating jets at higher latitudes. We present fully turbulent, 3D spherical numerical simulations of rapidly rotating convection with different spherical shell geometries. The resulting global flow fields tend to be segregated into three regions (north, equatorial, and south), bounded by the tangent cylinder that circumscribes the inner boundary equator. In all of our simulations a strong prograde equatorial jet forms outside the tangent cylinder, whereas multiple jets form in the northern and southern hemispheres, inside the tangent cylinder. The jet scaling of our numerical models and of Jupiter and Saturn is consistent with the theory of geostrophic turbulence, which we extend to include the effect of spherical shell geometry. Zonal flow in a spherical shell is distinguished from that in a full sphere or a shallow layer by the effect of the tangent cylinder, which marks a reversal in the sign of the planetary β-parameter and a jump in the Rhines length. This jump is manifest in the numerical simulations as a sharp equatorward increase in jet widths—a transition that is also observed on Jupiter and Saturn. The location of this transition gives an estimate of the depth of zonal flow, which seems to be consistent with current models of the jovian and saturnian interiors.     

Deep jets on gas-giant planets

Three-dimensional numerical simulations of the atmospheric flow on giant planets using the primitive equations show that shallow thermal forcing confined to pressures near the cloud tops can produce deep zonal winds from the tropopause all the way down to the bottom of the atmosphere. These deep winds can attain speeds comparable to the zonal jet speeds within the shallow, forced layer; they are pumped by Coriolis acceleration acting on a deep meridional circulation driven by the shallow-layer eddies. In the forced layer, the flow reaches an approximate steady state where east-west eddy accelerations balance Coriolis accelerations acting on the meridional flow. Under Jupiter-like conditions, our simulations produce 25 to 30 zonal jets, similar to the number of jets observed on Jupiter and Saturn. The simulated jet widths correspond to the Rhines scale; this suggests that, despite the three-dimensional nature of the dynamics, the baroclinic eddies energize a quasi-two-dimensional inverse cascade modified by the β effect (where β is the gradient of the Coriolis parameter). In agreement with Jupiter, the jets can violate the barotropic and Charney-Stern stability criteria, achieving curvatures ^∂2u/∂y² of the zonal wind u with northward distance y up to 2β. The simulations exhibit a tendency toward neutral stability with respect to Arnol'd's second stability theorem in the upper troposphere, as has been suggested for Jupiter, although deviations from neutrality exist. When the temperature varies strongly with latitude near the equator, our simulations can also reproduce the stable equatorial superrotation with wind speeds greater than . Diagnostics show that barotropic eddies at low latitudes drive the equatorial superrotation. The simulations also broadly explain the distribution of jet-pumping eddies observed on Jupiter and Saturn. While idealized, these simulations therefore capture many aspects of the cloud-level flows on Jupiter and Saturn.     

Zonal jets in rotating convection with mixed mechanical boundary conditions

Large-scale zonal flows, as observed on the giant planets, can be driven by thermal convection in a rapidly rotating spherical shell. Most previous models of convectively-driven zonal flow generation have utilized stress-free mechanical boundary conditions (FBC) for both the inner and the outer surfaces of the convecting layer. Here, using 3D numerical models, we compare the FBC case to the case with a stress free outer boundary and a non-slip inner boundary, which we call the mixed case (MBC). We find significant differences in surface zonal flow profiles produced by the two cases. In low to moderate Rayleigh number FBC cases, the main equatorial jet is flanked by a strong, high-latitude retrograde jets in the northern and southern hemispheres. For the highest Rayleigh number FBC case, the equatorial jet is flanked by strong reversed jets as well as two additional large-scale alternating jets at higher latitudes. The MBC cases feature stronger equatorial jets but, much weaker, small-scale alternating zonal flows are found at higher latitudes. Our high Rayleigh number FBC results best compare with the zonal flow pattern observed on Jupiter, where the equatorial jet is flanked by strong retrograde jets as well as small-scale alternating jets at high latitude. In contrast, the MBC results compare better with the observed flow pattern on Saturn, which is characterized by a dominant prograde equatorial jet and a lack of strong high latitude retrograde flow. This may suggest that the mechanical coupling at the base of the jovian convection zone differs from that on Saturn.     

A dynamo driven by zonal jets at the upper surface: Applications to giant planets

We present a dynamo mechanism arising from the presence of barotropically unstable zonal jet currents in a rotating spherical shell. The shear instability of the zonal flow develops in the form of a global Rossby mode, whose azimuthal wavenumber depends on the width of the zonal jets. We obtain self-sustained magnetic fields at magnetic Reynolds numbers greater than 10³. We show that the propagation of the Rossby waves is crucial for dynamo action. The amplitude of the axisymmetric poloidal magnetic field depends on the wavenumber of the Rossby mode, and hence on the width of the zonal jets. We discuss the plausibility of this dynamo mechanism for generating the magnetic field of the giant planets. Our results suggest a possible link between the topology of the magnetic field and the profile of the zonal winds observed at the surface of the giant planets. For narrow Jupiter-like jets, the poloidal magnetic field is dominated by an axial dipole whereas for wide Neptune-like jets, the axisymmetric poloidal field is weak.     

2- and 3-D Simulations of Magnetocentrifugal Disk-Winds: Acceleration and Stability

Protostellar jets and winds are probably driven magnetocentrifugally from the surface of accretion disks close to the central stellar objects. The exact launching conditions on the disk, such as the distributions of magnetic flux and mass ejection rate, are poorly known. They could be constrained from observations at large distances, provided that a robust model is available to link the observable properties of the jets and winds at the large distances to the conditions at the base of the flow. We describe a set of 2D axisymmetric simulations that are able to follow the acceleration and propagation of the wind from the disk surface to arbitrarily large distances. After a typical 2D flow reaches the steady state, we impose on it nonaxisymmetric perturbations and follow numerically its 3D evolution. We find that the wind reverts quickly to its initial axisymmetric state, with no indication of rapid growth of instabilities leading to flow disruption. Our calculations strengthen the case for the magnetocentrifugal jet and wind launching.     

Decadal variations in a Venus general circulation model

The Community Atmosphere Model (CAM), a 3-dimensional Earth-based climate model, has been modified to simulate the dynamics of the Venus atmosphere. The most current finite volume version of CAM is used with Earth-related processes removed, parameters appropriate for Venus introduced, and some basic physics approximations adopted. A simplified Newtonian cooling approximation has been used for the radiation scheme. We use a high resolution (1° by 1° in latitude and longitude) to take account of small-scale dynamical processes that might be important on Venus. A Rayleigh friction approach is used at the lower boundary to represent surface drag, and a similar approach is implemented in the uppermost few model levels providing a ‘sponge layer’ to prevent wave reflection from the upper boundary. The simulations generate superrotation with wind velocities comparable to those measured in the Venus atmosphere by probes and around 50-60% of those measured by cloud tracking. At cloud heights and above the atmosphere is always superrotating with mid-latitude zonal jets that wax and wane on an approximate 10 year cycle. However, below the clouds, the zonal winds vary periodically on a decadal timescale between superrotation and subrotation. Both subrotating and superrotating mid-latitude jets are found in the approximate 40-60 km altitude range. The growth and decay of the sub-cloud level jets also occur on the decadal timescale. Though subrotating zonal winds are found below the clouds, the total angular momentum of the atmosphere is always in the sense of superrotation. The global relative angular momentum of the atmosphere oscillates with an amplitude of about 5% on the approximate 10 year timescale. Symmetric instability in the near surface equatorial atmosphere might be the source of the decadal oscillation in the atmospheric state. Analyses of angular momentum transport show that all the jets are built up by poleward transport by a meridional circulation while angular momentum is redistributed to lower latitudes primarily by transient eddies. Possible changes in the structure of Venus’ cloud level mid-latitude jets measured by Mariner 10, Pioneer Venus, and Venus Express suggest that a cyclic variation similar to that found in the model might occur in the real Venus atmosphere, although no subrotating winds below the cloud level have been observed to date. Venus’ atmosphere must be observed over multi-year timescales and below the clouds if we are to understand its dynamics.     

Numerical simulations of astrophysical jets from Keplerian discs – III. The effects of mass loading

We present 2.5D time-dependent simulations of the non-linear evolution of non-relativistic outflows from the surface of Keplerian accretion discs. The gas is accelerated from the surface of the disc (which is a fixed platform in these simulations) into a cold corona in stable hydrostatic equilibrium. We explore the dependence of the resulting jet characteristics upon the mass loading of the winds. Two initial configurations of the threading disc magnetic field are studied: a potential field and a uniform vertical field configuration.
We show that the nature of the resulting highly collimated, jet-like outflows (steady or episodic) is determined by the mass load of the disc wind. The mass load controls the interplay between the collimating effects of the toroidal field and the kinetic energy density in the outflow. In this regard, we demonstrate that the onset of episodic behaviour of jets appears to be determined by the quantity     which compares the speed for a toroidal Alfvén wave to cross the diameter of the jet, with the flow speed v _p along the jet. This quantity decreases with increasing load. For sufficiently large N (small mass loads), disturbances appear to grow leading to instabilities and shocks. Knots are then generated and the outflow becomes episodic. These effects are qualitatively independent of the initial magnetic configuration that we employed and are probably generic to a wide variety of magnetized accretion disc models.     

Deep zonal winds can result from shallow driving in a giant-planet atmosphere

A simple model shows that acceleration of Jupiter and Saturn's multiple jets at altitudes confined near the top of the adiabatic region (e.g., at a few bars pressure) can produce jets that penetrate deeply into the molecular envelope. This result disproves the common assertion that jet acceleration near the outer margin can only produce zonal winds that are confined to these outer layers.     

Constraints on Saturn’s tropospheric general circulation from Cassini ISS images

An automated cloud tracking algorithm is applied to Cassini Imaging Science Subsystem high-resolution apoapsis images of Saturn from 2005 and 2007 and moderate resolution images from 2011 and 2012 to define the near-global distribution of zonal winds and eddy momentum fluxes at the middle troposphere cloud level and in the upper troposphere haze. Improvements in the tracking algorithm combined with the greater feature contrast in the northern hemisphere during the approach to spring equinox allow for better rejection of erroneous wind vectors, a more objective assessment at any latitude of the quality of the mean zonal wind, and a population of winds comparable in size to that available for the much higher contrast atmosphere of Jupiter. Zonal winds at cloud level changed little between 2005 and 2007 at all latitudes sampled. Upper troposphere zonal winds derived from methane band images are ～10 m s^?1 weaker than cloud level winds in the cores of eastward jets and ～5 m s^?1 stronger on either side of the jet core, i.e., eastward jets appear to broaden with increasing altitude. In westward jet regions winds are approximately the same at both altitudes. Lateral eddy momentum fluxes are directed into eastward jet cores, including the strong equatorial jet, and away from westward jet cores and weaken with increasing altitude on the flanks of the eastward jets, consistent with the upward broadening of these jets. The conversion rate of eddy to mean zonal kinetic energy at the visible cloud level is larger in eastward jet regions (5.2 × 10^?5 m² s^?3) and smaller in westward jet regions (1.6 × 10^?5 m² s^?3) than the global mean value (4.1 × 10^?5 m² s^?3). Overall the results are consistent with theories that suggest that the jets and the overturning meridional circulation at cloud level on Saturn are maintained at least in part by eddies due to instabilities of the large-scale flow near and/or below the cloud level.     

Simulations of high-latitude spots on Jupiter: Constraints on vortex strength and the deep wind

We combine high-resolution observations of the dynamical behavior of small vortices (diameters ?5000 km) located at latitude 60°N on Jupiter with forward modeling, using the EPIC atmospheric model, to address two open questions: the dependence of the zonal winds with depth, and the strength of vortices that are too small to apply cloud tracking to their internal structure. The observed drift rates of the vortices can only be reproduced in the model when the zonal winds increase slightly with depth below the cloud tops, with a vertical shear that is less than was measured at 7°N at the southern rim of a 5-μm hotspot by the Galileo Probe Doppler Wind Experiment (DWE). This supports the idea that Jupiter's vertical shear may vary significantly with latitude. Our simulations suggest that the morphology of the mergers between vortices mainly depends on their maximum tangential velocities, the best results occurring when the tangential velocity is close to the velocity difference of the alternating jets constraining the zone in which the vortices are embedded. We use this correlation, together with the high-resolution data available for the White Ovals, to derive an empirical relationship between the maximum tangential velocity of a jovian vortex and its size, normalized by the strength and size of the encompassing shear zone. The Great Red Spot stands out as a significant anomaly to this relationship, but interestingly it is becoming less so with time.     

Stability of jets on Jupiter and Saturn

The linear stability of a zonal jet that decays with depth is investigated under the assumption that the thermal stratification is very small. A westerly cosine jet is found to be more stable than it is in a thin fluid shell with two-dimensional flow. This is in agreement with observations of Jupiter and Saturn, where jet curvature exceeds the barotropic stability criterion. This result constitutes an alternative hypothesis to that of Ingersoll and Pollard [Icarus 52 (1982) 62], who showed that deep jets extending through the interior are also more stable than thin shell jets. The flow regime assumed in the present work requires that a small stratification can exist and persist even in the presence of horizontal temperature gradients. Further work will be needed to test whether this is realistic.     

The deep wind structure of the giant planets: Results from an anelastic general circulation model

The giant gas planets have hot convective interiors, and therefore a common assumption is that these deep atmospheres are close to a barotropic state. Here we show using a new anelastic general circulation model that baroclinic vorticity contributions are not negligible, and drive the system away from an isentropic and therefore barotropic state. The motion is still aligned with the direction of the axis of rotation as in a barotropic rotating fluid, but the wind structure has a vertical shear with stronger winds in the atmosphere than in the interior. This shear is associated with baroclinic compressibility effects. Most previous convection models of giant planets have used the Boussinesq approximation, which assumes the density is constant in depth; however, Jupiter's actual density varies by four orders of magnitude through its deep molecular envelope. We therefore developed a new general circulation model (based on the MITgcm) that is anelastic and thereby incorporates this density variation. The model's geometry is a full 3D sphere down to a small inner core. It is nonhydrostatic, uses an equation of state suitable for hydrogen-helium mixtures (SCVH), and is driven by an internal heating profile. We demonstrate the effect of compressibility by comparing anelastic and Boussinesq cases. The simulations develop a mean state that is geostrophic and hydrostatic including the often neglected, but significant, vertical Coriolis contribution. This leads to modification of the standard thermal wind relation for a deep compressible atmosphere. The interior flow organizes in large cyclonically rotating columnar eddies parallel to the rotation axis, which drive upgradient angular momentum eddy fluxes, generating the observed equatorial superrotation. Heat fluxes align with the axis of rotation, and provide a mechanism for the transport of heat poleward, which can cause the observed flat meridional emission. We address the issue of over-forcing which is common in such convection models and analyze the dependence of our results on this; showing that the vertical wind structure is not very sensitive to the Rayleigh number. We also study the effect of rotation, showing how the transition from a rapidly to a slowly rotating system affects the dynamics.     

A theoretical and empirical study of the response of the high latitude thermosphere to the sense of the “Y” component of the interplanetary magnetic field

The strength and direction of the Interplanetary Magnetic Field (IMF) controls the transfer of solar wind momentum and energy to the high latitude thermosphere in a direct fashion. The sense of “ Y” component of the IMF (BY) creates a significant asymmetry of the magnetospheric convection pattern as mapped onto the high latitude thermosphere and ionosphere. The resulting response of the polar thermospheric winds during periods when BY is either positive or negative is quite distinct, with pronounced changes in the relative strength of thermospheric winds in the dusk-dawn parts of the polar cap and in the dawn part of the auroral oval. In a study of four periods when there was a clear signature of BY, observed by the ISEE-3 satellite, with observations of polar winds and electric fields from the Dynamics Explorer-2 satellite and with wind observations by a ground-based Fabry-Perot interferometer located in Kiruna, Northern Sweden, it is possible to explain features of the high latitude thermospheric circulation using three dimensional global models including BY dependent, asymmetric, polar convection fields. Ground-based Fabry-Perot interferometers often observe anomalously low zonal wind velocities in the (Northern) dawn auroral oval during periods of extremely high geomagnetic activity when BY is positive. Conversely, for BY negative, there is an early transition from westward to southward and eastward winds in the evening auroral oval (excluding the effects of auroral substorms), and extremely large eastward (sunward) winds may be driven in the auroral oval after magnetic midnight. These observations are matched by the observation of strong anti-sunward polar-cap wind jets from the DE-2 satellite, on the dusk side with BY negative, and on the dawn side with BY positive.     

Steady and Time-Dependent MHD Modelling of Jets

A brief review is given of some results of our work on the construction of (I) steady and (II) time-dependent MHD models for nonrelativistic and relativistic astrophysical outflows and jets, analytically and numerically. The only available exact solutions for MHD outflows are those in separable coordinates, i.e., with the symmetry of radial or meridional self-similarity. Physically accepted solutions pass from the fast magnetosonic separatrix surface in order to satisfy MHD causality. An energetic criterion is outlined for selecting radially expanding winds from cylindrically expanding jets. Numerical simulations of magnetic self-collimation verify the conclusions of analytical steady solutions. We also propose a two-component model consisting of a wind outflow from a central object and a faster rotating outflow launched from a surrounding accretion disk which plays the role of the flow collimator. We also discuss the problem of shock formation during the magnetic collimation of wind-type outflows into jets.