Superrotation and Jet Migration in Simulations of Jupiter's Convective Zone and Weather Layer

By implementing 3D anelastic convection simulations of Jupiter that isolate or couple the deep convective zone with a shallow weather layer, the study reveals that Busse columns drive equatorial superrotation while potential vorticity homogenization generates high-latitude jets that slowly migrate poleward, with the weather layer significantly altering thermal wind balance and jet alignment.

The Gist

Imagine Jupiter as a giant, spinning ball of gas. If you look at it through a telescope, you see beautiful stripes: bands of wind blowing east and west, with a massive, fast-moving jet right at the equator. Scientists have long debated how these winds are created. Is it driven by heat from the Sun hitting the top of the atmosphere (a "shallow" process), or is it driven by heat rising from deep inside the planet (a "deep" process)?

This paper is like a virtual experiment where the authors built two digital models of Jupiter to see what happens when they turn on both the "deep" and "shallow" switches at the same time.

Here is the story of what they found, explained simply:

The Two Engines of Jupiter's Wind

Think of Jupiter's atmosphere as having two layers, like a two-story house:

1. The Deep Basement (The Convective Zone): This is the hot, churning interior. Here, the heat rises in giant columns of gas (aligned with Jupiter's axis of rotation) that twist as the planet spins. The authors call these "Busse columns." Imagine them like spinning tornadoes that stretch from the floor to the ceiling of the basement.
2. The Attic (The Weather Layer): This is the cool, stable top layer where we see the clouds. Here, the gas doesn't move up and down much; it just flows sideways in flat, pancake-like swirls.

The big question was: Do the basement columns or the attic pancakes create the stripes?

The Experiment: Two Simulations

The team ran two super-computer simulations:

• Simulation A: Just the basement (no attic).
• Simulation B: The basement plus a thin, stable attic layer on top.

What Happened?

1. The "Staircase" Effect (Making the Stripes)

In both simulations, the spinning gas naturally organized itself into multiple stripes (jets).

• How it works: Imagine the gas is trying to mix itself out evenly, like stirring sugar into coffee. But because the planet is spinning so fast, it can't mix everything smoothly. Instead, it creates "steps" or a "staircase" of different wind speeds.
• The Basement: The vertical columns create stripes that are aligned with the planet's axis (like rings on a tree trunk).
• The Attic: The flat pancakes create stripes that are aligned with the spherical surface (like rings on a sphere).
• The Result: In the early stages of the simulation, both layers successfully created multiple jets, just like we see on real Jupiter.

2. The Super-Strong Equator Jet

Both models produced a massive, fast jet right at the equator that spins faster than the planet itself (called "superrotation").

• The Basement's Role: The authors found that the vertical columns in the basement act like a conveyor belt. Because the planet is round, these columns flare out slightly as they go up. This flare pushes angular momentum (spin energy) outward toward the equator, creating the super-fast jet.
• The Attic's Role: In the model with the attic, the attic didn't create its own superrotating jet. Instead, it just "caught" the fast spin from the basement below it, like a person on a merry-go-round grabbing onto a spinning pole. The attic's wind was just an echo of the basement's wind.

3. The Long Wait (The Migration Problem)

This is the most surprising part.

• The Early Days: At the start of the simulation, the models looked perfect. They had many stripes, just like Jupiter.
• The Long Haul: The authors ran the simulations for a very long time (substantially longer than previous studies, for thousands of rotation periods). They discovered that the high-latitude stripes (the ones near the poles) are not stable.
• The Drift: Over time, these smaller stripes slowly drifted toward the poles and merged with each other. It's like a crowd of people walking in a circle; eventually, they bump into each other and merge into fewer, larger groups.
• The Final State: After a very long time, the models settled into a state with only three jets per hemisphere: one fast one at the equator, and one slow/fast pair near the poles.

The Big Takeaway

The paper suggests that while the "shallow" (attic) and "deep" (basement) layers can both create wind stripes, the deep layer is the real boss of the equatorial superrotating jet.

However, there is a mystery. The authors found that in their 3D models, the multiple stripes near the poles eventually disappear and merge. This implies that the Jupiter we see today (with its many stripes) might be in a temporary state, or that our current computer models are missing a specific "brake" or friction force that stops the stripes from merging.

In short: The authors built a digital Jupiter to see how its winds form. They found that deep columns and shallow pancakes both help make stripes, but the deep columns drive the super-fast equator wind. However, their models showed that the smaller stripes near the poles are unstable and tend to merge over time, suggesting that keeping Jupiter's many stripes requires a delicate balance we are still trying to understand.

Technical Summary

Technical Summary: Superrotation and Jet Migration in Simulations of Jupiter's Convective Zone and Weather Layer

Problem Statement
The mean zonal flow observed on Jupiter consists of an intricate pattern of alternating prograde and retrograde jets, capped by a strong superrotating (prograde) equatorial jet. Theoretical understanding of these flows has traditionally been divided into two competing frameworks: "shallow" driving, where baroclinic turbulence in a stably stratified weather layer (WL) generates jets via quasi-two-dimensional (2D) dynamics, and "deep" driving, where buoyantly driven, rotationally constrained three-dimensional (3D) convection in the convective zone (CZ) generates jets via Busse columns. While both mechanisms are thought to be active, previous 3D simulations have often failed to reach a statistically steady state, leaving the long-term evolution of high-latitude jets and the specific influence of a stable WL on deep convection unquantified. This study aims to unify 2D and 3D concepts by performing fully nonlinear, 3D simulations of a Jovian-like planet to assess the ultimate steady state of zonal flows and the distinct roles of potential vorticity (PV) homogenization in different layers.

Methodology
The authors implement two rotating, 3D, spherical-shell, anelastic convection simulations using the open-source Rayleigh code. The simulations model a rapidly rotating planet ($Ro \ll 1$, where $Ro$, the Rossby number, measures the influence of rotation on the turbulent eddies) with two distinct configurations:

1. CZ-only case: An isolated convective zone with impenetrable boundaries.
2. CZ–WL case: A convective zone overlaid by a thin, stably stratified idealized weather layer, allowing convective overshoot.

The models utilize a flux-based Rayleigh number ($Ra_f \approx 10^7$) and a convective Rossby number ($Ro_c \approx 0.2$). The simulations are run to statistical equilibrium, a process requiring $\mathcal{O}(10)$ viscous diffusion timescales ($1.8 \times 10^4$ rotation periods), significantly longer than many prior studies. The analysis focuses on the early-time formation of jets ($t \approx 5000$ simulation time units, or $t \approx 400$ rotation periods) and the long-term migration and equilibration of the flow structure. The study explicitly tracks the evolution of mean zonal kinetic energy, meridional circulation, and the homogenization of distinct forms of potential vorticity: $Q_z$ (axial PV) in the CZ and $Q_r$ (radial PV) in the WL.

Key Contributions and Results

• Distinct Forms of Potential Vorticity and Jet Formation:
  The study confirms that the appropriate form of PV depends on the symmetry of the eddies. In the CZ, Busse columns (axially symmetric) drive the homogenization of axial potential vorticity $Q_z = (\omega_z + 2\Omega_0)/H$. In the WL, pancake eddies (radially symmetric) drive the homogenization of radial potential vorticity $Q_r = \omega_r + f$. In both layers, the initial formation of multiple high-latitude jets is consistent with the "PV staircase" mechanism, where jets form at the Rhines scale ($L_R$) separating regions of homogenized PV. The local width of the jets matches the local Rhines scale well in the high-latitude regions of both layers.

• Mechanisms of Superrotation:
  The equatorial superrotation is found to be driven primarily by the transport of angular momentum by Busse columns in the CZ, rather than solely by PV homogenization. The convex spherical boundaries cause Busse columns to flare outward, creating a correlation between fluctuating radial and zonal velocities that transports angular momentum outward (Reynolds stress torque). This mechanism is robust in the CZ. In the CZ–WL case, the superrotation in the WL is not driven by local Reynolds stresses or PV homogenization but is instead an "imprint" of the CZ superrotation, transmitted upward via meridional circulation and viscous stresses.

• Impact of the Weather Layer on Thermal Wind Balance:
  The presence of the idealized weather layer significantly alters the thermal wind balance. While the CZ jets remain largely cylindrically aligned, the WL exhibits a monocellular meridional circulation that tilts the zonal flow isocontours away from the rotation axis. This results in a large-scale baroclinicity in the WL that modifies the thermal wind balance, demonstrating that even an idealized stable layer (one with a single molecular species and no energy input from the Sun) can decouple the alignment of deep and shallow flows.

• Jet Migration and the Asymptotic Steady State:
  A critical finding is the long-term evolution of the flow. While the equatorial superrotation equilibrates relatively quickly ($\mathcal{O}(1)$ diffusion time), the high-latitude jets are unsteady on long timescales. Over $\mathcal{O}(10)$ diffusion times, the multiple high-latitude jets slowly migrate poleward and/or merge. The final asymptotic state consists of a "fast-slow-fast" banded structure: one strong superrotating jet outside the tangent cylinder and only one pair of alternating jets per hemisphere inside the tangent cylinder. This suggests that the multiple jets observed in earlier, shorter-duration 3D simulations may represent transient states rather than the true equilibrium.

Significance and Claims
The paper claims that the "deep versus shallow driving" debate may not require a single resolution, as both layers can independently generate multiple jets via their respective PV homogenization mechanisms. However, the coupling between the layers is crucial for determining the alignment of the zonal flows. The authors assert that the equatorial superrotation is fundamentally a 3D effect driven by the geometry of Busse columns (flaring), which cannot be fully captured by simple quasi-2D PV homogenization arguments alone.

Furthermore, the study highlights a potential limitation in current state-of-the-art 3D simulations of Jupiter: the high-latitude jets in many prior models were likely still evolving. The authors suggest that without a physical dissipation mechanism (such as a scale-independent drag) to halt the inverse energy cascade at large scales, 3D simulations may naturally evolve toward a state with fewer high-latitude jets than are currently observed on Jupiter. The paper concludes that achieving a steady state with multiple high-latitude jets in 3D simulations remains a challenge, and the true mechanism sustaining Jupiter's observed jet count may involve physics not yet fully incorporated into global simulations.