Search

In order to characterize the global circulation of the subsurface ocean of Jovian and Saturnian moons, we analyze the properties of 21 three-dimensional simulations of Boussinesq thermal convection in a rapidly rotating spherical shell. Flow is driven by an adverse temperature contrast imposed across the domain, and is subjected to no-slip boundary conditions. We cover a region of parameter space previously unexplored by global simulations, both in terms of rapid rotation and vigor of convective forcing, closer to, yet still admittedly far from, the conditions appropriate for the subsurface ocean of Ganymede, Europa, Enceladus, and Titan. Our most extreme simulations exhibit a dynamic global circulation that combines powerful east-west zonal jets, planetary waves, and vortices. A spectral analysis of the kinetic energy distribution performed in cylindrical geometry reveals a high degree of anisotropy of the simulated flows. Specifically, the axisymmetric zonal energy spectra follow a steep $-5$ slope in wavenumber space, with the energy amplitude exclusively controlled by the rotation rate. In contrast, the non-axisymmetric residual spectra display a gentle $-5/3$ slope, with the energy amplitude controlled by the thermal buoyancy input power. This spectral behavior conforms with the theory of zonostrophic turbulence and allows us to propose tentative extrapolations of these findings to the more extreme conditions of icy satellites. By assuming that kinetic energy dissipates via Ekman friction we predict an upper bound for the zonal velocity ranging from a few centimeters per second for Enceladus to about one meter per second for Ganymede, with residual velocities smaller than the zonal velocity by an order of magnitude on each moon. These predictions yield typical jets size approaching the ocean depth of Titan, Ganymede and Europa and $10$ to $40\%$ of the ocean depth on Enceladus., Comment: 24 pages, 12 figures, 3 tables, accepted for publication in Icarus

We use 123 three dimensional direct numerical simulations to study fingering convection in non-rotating spherical shells. We investigate the scaling behaviour of the flow lengthscale, the non-dimensional heat and compositional fluxes $Nu$ and $Sh$ and the mean convective velocity over the fingering convection instability domain defined by $1 \leq R_\rho < Le$, $R_\rho$ being the ratio of density perturbations of thermal and compositional origins and $Le$ the Lewis number. We show that the chemical boundary layers are marginally unstable and adhere to the laminar Prandtl-Blasius model, hence explaining the asymmetry between the inner and outer spherical shell boundary layers. We develop scaling laws for two asymptotic regimes close to the two edges of the instability domain, namely $R_\rho \lesssim Le$ and $R_\rho \gtrsim 1$. For the former, we develop novel power laws of a small parameter $\epsilon$ measuring the distance to onset, which differ from theoretical laws published to date in Cartesian geometry. For the latter, we find that the Sherwood number $Sh$ gradually approaches a scaling $Sh\sim Ra_\xi^{1/3}$ when $Ra_\xi \gg 1$; and that the P\'eclet number accordingly follows $Pe \sim Ra_\xi^{2/3} |Ra_T|^{-1/4}$, $Ra_\xi$ being the chemical Rayleigh number. When the Reynolds number exceeds a few tens, we report on a secondary instability which takes the form of large-scale toroidal jets which span the entire spherical domain. Jets distort the fingers resulting in Reynolds stress correlations, which in turn feed the jet growth until saturation. This nonlinear phenomenon can yield relaxation oscillation cycles., Comment: 45 pages, 21 figures, 3 tables, accepted for publication in JFM

Convection occurs ubiquitously on and in rotating geophysical and astrophysical bodies. Prior spherical shell studies have shown that the convection dynamics in polar regions can differ significantly from the lower latitude, equatorial dynamics. Yet most spherical shell convective scaling laws use globally-averaged quantities that erase latitudinal differences in the physics. Here we quantify those latitudinal differences by analyzing spherical shell simulations in terms of their regionalized convective heat transfer properties. This is done by measuring local Nusselt numbers in two specific, latitudinally separate, portions of the shell, the polar and the equatorial regions, $Nu_p$ and $Nu_e$, respectively. In rotating spherical shells, convection first sets in outside the tangent cylinder such that equatorial heat transfer dominates at small and moderate supercriticalities. We show that the buoyancy forcing, parameterized by the Rayleigh number $Ra$, must exceed the critical equatorial forcing by a factor of $\approx 20$ to trigger polar convection within the tangent cylinder. Once triggered, $Nu_p$ increases with $Ra$ much faster than does $Nu_e$. The equatorial and polar heat fluxes then tend to become comparable at sufficiently high $Ra$. Comparisons between the polar convection data and Cartesian numerical simulations reveal quantitative agreement between the two geometries in terms of heat transfer and averaged bulk temperature gradient. This agreement indicates that spherical shell rotating convection dynamics are accessible both through spherical simulations and via reduced investigatory pathways, be they theoretical, numerical or experimental., Comment: 11 pages, 6 figures, accepted for publication in JFM

The ongoing NASA's Juno mission puts new constraints on the internal dynamics of Jupiter. Data gathered by its onboard magnetometer reveal a dipole-dominated surface magnetic field accompanied by strong localised magnetic flux patches. The gravity measurements indicate that the fierce surface zonal jets extend several thousands of kilometers below the cloud level before rapidly decaying below $0.94-0.96\,R_J$, $R_J$ being the mean Jovian radius at the one bar level. Several internal models suggest an intricate internal structure with a thin intermediate region in which helium would segregate from hydrogen, forming a compositionally-stratified layer. Here, we develop the first global Jovian dynamo which incorporates an intermediate stably-stratified layer between $0.82\,R_J$ and $0.86\,R_J$. Analysing the energy balance reveals that the magnetic energy is almost one order of magnitude larger than kinetic energy in the metallic region, while most of the kinetic energy is pumped into zonal motions in the molecular envelope. Those result from the different underlying force hierarchy with a triple balance between Lorentz, Archimedean and ageostrophic Coriolis forces in the metallic core and inertia, buoyancy and ageostrophic Coriolis forces controlling the external layers. The simulation presented here is the first to demonstrate that multiple zonal jets and dipole-dominated dynamo action can be consolidated in a global simulation. The inclusion of an stable layer is a necessary ingredient that allows zonal jets to develop in the outer envelope without contributing to the dynamo action in the deeper metallic region. Stable stratification however also smooths out the small-scale features of the magnetic field by skin effect. These constraints suggest that possible stable layers in Jupiter should be located much closer to the surface ($0.9-0.95\,R_J$)., Comment: 23 pages, 16 figures, 2 tables, accepted for publication in Icarus

We present a new pseudo-spectral open-source code nicknamed pizza. It is dedicated to the study of rapidly-rotating Boussinesq convection under the 2-D spherical quasi-geostrophic approximation, a physical hypothesis that is appropriate to model the turbulent convection that develops in planetary interiors. The code uses a Fourier decomposition in the azimuthal direction and supports both a Chebyshev collocation method and a sparse Chebyshev integration formulation in the cylindrically-radial direction. It supports several temporal discretisation schemes encompassing multi-step time steppers as well as diagonally-implicit Runge-Kutta schemes. The code has been tested and validated by comparing weakly-nonlinear convection with the eigenmodes from a linear solver. The comparison of the two radial discretisation schemes has revealed the superiority of the Chebyshev integration method over the classical collocation approach both in terms of memory requirements and operation counts. The good parallelisation efficiency enables the computation of large problem sizes with $\mathcal{O}(10^4\times 10^4)$ grid points using several thousands of ranks. This allows the computation of numerical models in the turbulent regime of quasi-geostrophic convection characterised by large Reynolds $Re$ and yet small Rossby numbers $Ro$. A preliminary result obtained for a strongly supercritical numerical model with a small Ekman number of $10^{-9}$ and a Prandtl number of unity yields $Re\simeq 10^5$ and $Ro \simeq 10^{-4}$. pizza is hence an efficient tool to study spherical quasi-geostrophic convection in a parameter regime inaccessible to current global 3-D spherical shell models., Comment: 20 pages, 10 figures, 4 tables, accepted for publication in GJI

Rayleigh-B\'enard convection in rotating spherical shells can be considered as a simplified analogue of many astrophysical and geophysical fluid flows. Here, we use three-dimensional direct numerical simulations to study this physical process. We construct a dataset of more than 200 numerical models that cover a broad parameter range with Ekman numbers spanning $3\times 10^{-7} \leq E \leq 10^{-1}$, Rayleigh numbers within the range $10^3 < Ra < 2\times 10^{10}$ and a Prandtl number unity. We investigate the scaling behaviours of both local (length scales, boundary layers) and global (Nusselt and Reynolds numbers) properties across various physical regimes from onset of rotating convection to weakly-rotating convection. Close to critical, the convective flow is dominated by a triple force balance between viscosity, Coriolis force and buoyancy. For larger supercriticalities, a subset of our numerical data approaches the asymptotic diffusivity-free scaling of rotating convection $Nu\sim Ra^{3/2}E^{2}$ in a narrow fraction of the parameter space delimited by $6\,Ra_c \leq Ra \leq 0.4\,E^{-8/5}$. Using a decomposition of the viscous dissipation rate into bulk and boundary layer contributions, we establish a theoretical scaling of the flow velocity that accurately describes the numerical data. In rapidly-rotating turbulent convection, the fluid bulk is controlled by a triple force balance between Coriolis, inertia and buoyancy, while the remaining fraction of the dissipation can be attributed to the viscous friction in the Ekman layers. Beyond $Ra \simeq E^{-8/5}$, the rotational constraint on the convective flow is gradually lost and the flow properties vary to match the regime changes between rotation-dominated and non-rotating convection. The quantity $Ra E^{12/7}$ provides an accurate transition parameter to separate rotating and non-rotating convection., Comment: 42 pages, 20 figures, 3 tables, accepted for publication in JFM

We simulate numerically Boussinesq convection in non-rotating spherical shells for a fluid with a unity Prandtl number and Rayleigh numbers up to $10^9$. In this geometry, curvature and radial variations of the gravitationnal acceleration yield asymmetric boundary layers. A systematic parameter study for various radius ratios (from $\eta=r_i/r_o=0.2$ to $\eta=0.95$) and gravity profiles allows us to explore the dependence of the asymmetry on these parameters. We find that the average plume spacing is comparable between the spherical inner and outer bounding surfaces. An estimate of the average plume separation allows us to accurately predict the boundary layer asymmetry for the various spherical shell configurations explored here. The mean temperature and horizontal velocity profiles are in good agreement with classical Prandtl-Blasius laminar boundary layer profiles, provided the boundary layers are analysed in a dynamical frame, that fluctuates with the local and instantaneous boundary layer thicknesses. The scaling properties of the Nusselt and Reynolds numbers are investigated by separating the bulk and boundary layer contributions to the thermal and viscous dissipation rates using numerical models with $\eta=0.6$ and a gravity proportional to $1/r^2$. We show that our spherical models are consistent with the predictions of Grossmann \& Lohse's (2000) theory and that $Nu(Ra)$ and $Re(Ra)$ scalings are in good agreement with plane layer results., Comment: 43 pages, 25 figures, 2 tables, accepted for publication in JFM

Spacecraft data reveal a very Earth-like Jovian magnetic field. This is surprising since numerical simulations have shown that the vastly different interiors of terrestrial and gas planets can strongly affect the internal dynamo process. Here we present the first numerical dynamo that manages to match the structure and strength of the observed magnetic field by embracing the newest models for Jupiter's interior. Simulated dynamo action primarily occurs in the deep high electrical conductivity region while zonal flows are dynamically constrained to a strong equatorial jet in the outer envelope of low conductivity. Our model reproduces the structure and strength of the observed global magnetic field and predicts that secondary dynamo action associated to the equatorial jet produces banded magnetic features likely observable by the Juno mission. Secular variation in our model scales to about 2000 nT per year and should also be observable during the one year nominal mission duration., Comment: 7 pages, 4 figures, accepted for publication in Geophysical Research Letters

The surface winds of Jupiter and Saturn are primarily zonal. Each planet exhibits strong prograde equatorial flow flanked by multiple alternating zonal winds at higher latitudes. The depth to which these flows penetrate has long been debated and is still an unsolved problem. Previous rotating convection models that obtained multiple high latitude zonal jets comparable to those on the giant planets assumed an incompressible (Boussinesq) fluid, which is unrealistic for gas giant planets. Later models of compressible rotating convection obtained only few high latitude jets which were not amenable to scaling analysis. Here we present 3-D numerical simulations of compressible convection in rapidly-rotating spherical shells. To explore the formation and scaling of high-latitude zonal jets, we consider models with a strong radial density variation and a range of Ekman numbers, while maintaining a zonal flow Rossby number characteristic of Saturn. All of our simulations show a strong prograde equatorial jet outside the tangent cylinder. At low Ekman numbers several alternating jets form in each hemisphere inside the tangent cylinder. To analyse jet scaling of our numerical models and of Jupiter and Saturn, we extend Rhines scaling based on a topographic $\beta$-parameter, which was previously applied to an incompressible fluid in a spherical shell, to compressible fluids. The jet-widths predicted by this modified Rhines length are found to be in relatively good agreement with our numerical model results and with cloud tracking observations of Jupiter and Saturn., Comment: 17 pages, 12 figures, 3 tables, accepted for publication in PEPI

The surface zonal winds observed in the giant planets form a complex jet pattern with alternating prograde and retrograde direction. While the main equatorial band is prograde on the gas giants, both ice giants have a pronounced retrograde equatorial jet. We use three-dimensional numerical models of compressible convection in rotating spherical shells to explore the properties of zonal flows in different regimes where either rotation or buoyancy dominates the force balance. We conduct a systematic parameter study to quantify the dependence of zonal flows on the background density stratification and the driving of convection. We find that the direction of the equatorial zonal wind is controlled by the ratio of buoyancy and Coriolis force. The prograde equatorial band maintained by Reynolds stresses is found in the rotation-dominated regime. In cases where buoyancy dominates Coriolis force, the angular momentum per unit mass is homogenised and the equatorial band is retrograde, reminiscent to those observed in the ice giants. In this regime, the amplitude of the zonal jets depends on the background density contrast with strongly stratified models producing stronger jets than comparable weakly stratified cases. Furthermore, our results can help to explain the transition between solar-like and "anti-solar" differential rotations found in anelastic models of stellar convection zones. In the strongly stratified cases, we find that the leading order force balance can significantly vary with depth (rotation-dominated inside and buoyancy-dominated in a thin surface layer). This so-called "transitional regime" has a visible signature in the main equatorial jet which shows a pronounced dimple where flow amplitudes notably decay towards the equator. A similar dimple is observed on Jupiter, which suggests that convection in the planet interior could possibly operate in this regime., Comment: 20 pages, 15 figures, 4 tables, accepted for publication in Icarus

Context: dynamo action in giant planets and rapidly rotating stars leads to a broad variety of magnetic field geometries including small scale multipolar and large scale dipole-dominated topologies. Previous dynamo models suggest that solutions become multipolar once inertia becomes influential. Being tailored for terrestrial planets, most of these models neglected the background density stratification. Aims: we investigate the influence of the density stratification on convection-driven dynamo models. Methods: three-dimensional nonlinear simulations of rapidly rotating spherical shells are employed using the anelastic approximation to incorporate density stratification. A systematic parametric study for various density stratifications and Rayleigh numbers allows to explore the dependence of the magnetic field topology on these parameters. Results: anelastic dynamo models tend to produce a broad range of magnetic field geometries that fall on two distinct branches with either strong dipole-dominated or weak multipolar fields. As long as inertia is weak, both branches can coexist but the dipolar branch vanishes once inertia becomes influential. The dipolar branch also vanishes for stronger density stratifications. The reason is the concentration of the convective columns in a narrow region close to the outer boundary equator, a configuration that favors non-axisymmetric solutions. In multipolar solutions, zonal flows can become significant and participate in the toroidal field generation. Parker dynamo waves may then play an important role close to onset of dynamo action leading to a cyclic magnetic field behavior. Conclusion: Our simulations also suggest that the fact that late M dwarfs have dipolar or multipolar magnetic fields can be explained in two ways. They may differ either by the relative influence of inertia or fall into the regime where both types of solutions coexist., Comment: 13 pages, 13 figures, 2 tables, accepted for publication in A&A

The banded structures observed on the surfaces of the gas giants are associated with strong zonal winds alternating in direction with latitude. We use three-dimensional numerical simulations of compressible convection in the anelastic approximation to explore the properties of zonal winds in rapidly rotating spherical shells. Since the model is restricted to the electrically insulating outer envelope, we therefore neglect magnetic effects. A systematic parametric study for various density scaleheights and Rayleigh numbers allows to explore the dependence of convection and zonal jets on these parameters and to derive scaling laws. While the density stratification affects the local flow amplitude and the convective scales, global quantities and zonal jets properties remain fairly independent of the density stratification. The zonal jets are maintained by Reynolds stresses, which rely on the correlation between zonal and cylindrically radial flow components. The gradual loss of this correlation with increasing supercriticality hampers all our simulations and explains why the additional compressional source of vorticity hardly affects zonal flows. All these common features may explain why previous Boussinesq models were already successful in reproducing the morphology of zonal jets in gas giants., Comment: 17 pages, 15 figures, 2 tables, accepted for publication in Icarus