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

En mettant en œuvre des simulations de convection anélastique en 3D de Jupiter qui isolent ou couplent la zone de convection profonde avec une couche météorologique peu profonde, l'étude révèle que les colonnes de Busse entraînent une superrotation équatoriale tandis que l'homogénéisation de la vorticité potentielle génère des jets aux hautes latitudes qui migrent lentement vers les pôles, la couche météorologique modifiant considérablement l'équilibre du vent thermique et l'alignement des jets.

L'essentiel

Imaginez Jupiter comme une gigantesque boule de gaz en rotation. Si vous l'observez à travers un télescope, vous voyez de magnifiques bandes : des courants de vent soufflant d'est en ouest, avec un puissant jet rapide juste à l'équateur. Les scientifiques débattent depuis longtemps de la manière dont ces vents sont créés. Sont-ils alimentés par la chaleur du Soleil frappant le sommet de l'atmosphère (un processus « superficiel »), ou par la chaleur remontant des profondeurs de la planète (un processus « profond ») ?

Cet article est comme une expérience virtuelle où les auteurs ont construit deux modèles numériques de Jupiter pour observer ce qui se produit lorsqu'ils activent simultanément les interrupteurs « profond » et « superficiel ».

Voici l'histoire de leurs découvertes, expliquée simplement :

Les deux moteurs du vent de Jupiter

Imaginez l'atmosphère de Jupiter comme ayant deux couches, telles qu'une maison à deux étages :

1. Le sous-sol profond (la zone de convection) : C'est l'intérieur chaud et tourbillonnant. Ici, la chaleur s'élève en d'immenses colonnes verticales de gaz qui se tordent à mesure que la planète tourne. Les auteurs les appellent des « colonnes de Busse ». Imaginez-les comme des tornades en rotation qui s'étendent du sol au plafond du sous-sol.
2. Le grenier (la couche météorologique) : C'est la couche supérieure fraîche et stable où nous voyons les nuages. Ici, l'air ne monte ni ne descend beaucoup ; il s'écoule simplement sur le côté dans des tourbillons plats et plats comme des crêpes.

La grande question était : sont-ce les colonnes du sous-sol ou les crêpes du grenier qui créent les bandes ?

L'expérience : deux simulations

L'équipe a exécuté deux simulations sur superordinateur :

• Simulation A : Juste le sous-sol (pas de grenier).
• Simulation B : Le sous-sol plus une fine couche de grenier stable au-dessus.

Qu'est-il arrivé ?

1. L'effet « escalier » (création des bandes)

Dans les deux simulations, le gaz en rotation s'est naturellement organisé en multiples bandes (jets).

• Comment cela fonctionne : Imaginez que le gaz essaie de se mélanger uniformément, comme si l'on remuait du sucre dans du café. Mais parce que la planète tourne si vite, elle ne peut pas tout mélanger de manière fluide. Au lieu de cela, elle crée des « marches » ou un « escalier » de différentes vitesses de vent.
• Le sous-sol : Les colonnes verticales créent des bandes alignées avec l'axe de la planète (comme les anneaux sur un tronc d'arbre).
• Le grenier : Les crêpes plates créent des bandes alignées avec la surface (comme les anneaux sur une sphère).
• Le résultat : Dans les premiers stades de la simulation, les deux couches ont réussi à créer plusieurs jets, tout comme nous les observons sur Jupiter réelle.

2. Le jet équatorial ultra-puissant

Les deux modèles ont produit un jet massif et rapide juste à l'équateur, tournant plus vite que la planète elle-même (appelé « superrotation »).

• Le rôle du sous-sol : Les auteurs ont constaté que les colonnes verticales du sous-sol agissent comme un tapis roulant. Parce que la planète est ronde, ces colonnes s'évasent légèrement en montant. Cet évasement pousse le moment angulaire (énergie de rotation) vers l'extérieur en direction de l'équateur, créant le jet ultra-rapide.
• Le rôle du grenier : Dans le modèle avec grenier, ce dernier n'a pas créé son propre super-jet. Au lieu de cela, il a simplement « attrapé » la rotation rapide du sous-sol en dessous, comme une personne sur un manège qui s'accroche à un poteau en rotation. Le vent du grenier n'était qu'un écho du vent du sous-sol.

3. La longue attente (le problème de la migration)

C'est la partie la plus surprenante.

• Les premiers jours : Au début de la simulation, les modèles semblaient parfaits. Ils possédaient de nombreuses bandes, tout comme Jupiter.
• La longue haleine : Les auteurs ont fait tourner les simulations pendant un temps très long (des milliers de fois plus longtemps que les études précédentes). Ils ont découvert que les bandes de haute latitude (celles près des pôles) ne sont pas stables.
• La dérive : Au fil du temps, ces petites bandes ont lentement dérivé vers les pôles et fusionné entre elles. C'est comme une foule de personnes marchant en cercle ; éventuellement, elles se heurtent les unes aux autres et fusionnent en moins de groupes, plus grands.
• L'état final : Après un temps très long, les modèles se sont stabilisés dans un état avec seulement trois jets par hémisphère : un rapide à l'équateur, et deux plus lents près des pôles.

La grande conclusion

L'article suggère que, bien que les couches « superficielles » (grenier) et « profondes » (sous-sol) puissent toutes deux créer des bandes de vent, la couche profonde est le véritable chef du jet équatorial super-rapide.

Cependant, il y a un mystère. Les auteurs ont constaté que dans leurs modèles 3D, les multiples bandes près des pôles finissent par disparaître et fusionner. Cela implique que Jupiter telle que nous la voyons aujourd'hui (avec ses nombreuses bandes) pourrait être dans un état temporaire, ou que nos modèles informatiques actuels manquent d'un « frein » ou d'une force de friction spécifique qui empêche les bandes de fusionner.

En bref : Les auteurs ont construit un Jupiter numérique pour voir comment ses vents se forment. Ils ont découvert que les colonnes profondes et les crêpes superficielles aident toutes deux à créer des bandes, mais que les colonnes profondes entraînent le vent équatorial ultra-rapide. Cependant, leurs modèles ont montré que les petites bandes près des pôles sont instables et ont tendance à fusionner au fil du temps, suggérant que le maintien des nombreuses bandes de Jupiter nécessite un équilibre délicat que nous essayons encore de comprendre.

Résumé technique

Résumé technique : Superrotation et migration des jets dans les simulations de la zone convective et de la couche météorologique de Jupiter

Énoncé du problème
L'écoulement zonal moyen observé sur Jupiter consiste en un motif complexe de jets alternés progrades et rétrogrades, surmonté d'un puissant jet équatorial superrotant (prograde). La compréhension théorique de ces écoulements a traditionnellement été divisée en deux cadres concurrents : la « conduite » peu profonde, où la turbulence barocline dans une couche météorologique (WL) stablement stratifiée génère des jets via une dynamique quasi bidimensionnelle (2D), et la « conduite » profonde, où la convection tridimensionnelle (3D) entraînée par la flottabilité et contrainte par la rotation dans la zone convective (CZ) génère des jets via des colonnes de Busse. Bien que les deux mécanismes soient considérés comme actifs, les simulations 3D précédentes ont souvent échoué à atteindre un état statistiquement stationnaire, laissant l'évolution à long terme des jets de haute latitude et l'influence spécifique d'une WL stable sur la convection profonde non quantifiées. Cette étude vise à unifier les concepts 2D et 3D en réalisant des simulations 3D entièrement non linéaires d'une planète de type jovien pour évaluer l'état stationnaire ultime des écoulements zonaux et les rôles distincts de l'homogénéisation de la vorticité potentielle (PV) dans différentes couches.

Méthodologie
Les auteurs mettent en œuvre deux simulations de convection anélastique en coquille sphérique, 3D et en rotation, utilisant le code open-source Rayleigh. Les simulations modélisent une planète en rotation rapide ($Ro \ll 1$) avec deux configurations distinctes :

1. Cas CZ uniquement : Une zone convective isolée avec des frontières imperméables.
2. Cas CZ–WL : Une zone convective surmontée d'une fine couche météorologique idéalisée stablement stratifiée, permettant un dépassement convectif.

Les modèles utilisent un nombre de Rayleigh basé sur le flux ($Ra_f \approx 10^7$) et un nombre de Rossby convectif ($Ro_c \approx 0.2$). Les simulations sont menées jusqu'à l'équilibre statistique, un processus nécessitant $\mathcal{O}(10)$ échelles de temps de diffusion visqueuse ($\sim 2.2 \times 10^5$ périodes de rotation), ce qui est significativement plus long que dans de nombreuses études antérieures. L'analyse se concentre sur la formation précoce des jets ($t \approx 5000$) et la migration à long terme ainsi que l'équilibration de la structure de l'écoulement. L'étude suit explicitement l'évolution de l'énergie cinétique zonale moyenne, de la circulation méridienne et de l'homogénéisation de formes distinctes de vorticité potentielle : $Q_z$ (PV axiale) dans la CZ et $Q_r$ (PV radiale) dans la WL.

Contributions et résultats clés

• Formes distinctes de vorticité potentielle et formation des jets :
  L'étude confirme que la forme appropriée de PV dépend de la symétrie des tourbillons. Dans la CZ, les colonnes de Busse (symétriques axialement) entraînent l'homogénéisation de la vorticité potentielle axiale $Q_z = \omega_z + 2\Omega_0/H$. Dans la WL, les tourbillons en forme de crêpes (symétriques radialement) entraînent l'homogénéisation de la vorticité potentielle radiale $Q_r = \omega_r + f$. Dans les deux couches, la formation initiale de multiples jets de haute latitude est cohérente avec le mécanisme de « l'escalier de PV », où les jets se forment à l'échelle de Rhines ($L_R$) séparant les régions de PV homogénéisée. La largeur locale des jets correspond bien à l'échelle locale de Rhines dans les régions de haute latitude des deux couches.

• Mécanismes de la superrotation :
  Il est constaté que la superrotation équatoriale est principalement entraînée par le transport de moment angulaire par les colonnes de Busse dans la CZ, plutôt que par l'homogénéisation de la PV seule. Les limites sphériques convexes provoquent l'évasement vers l'extérieur des colonnes de Busse, créant une corrélation entre les vitesses radiales et zonales fluctuantes qui transporte le moment angulaire vers l'extérieur (couple de contrainte de Reynolds). Ce mécanisme est robuste dans la CZ. Dans le cas CZ–WL, la superrotation dans la WL n'est pas entraînée par des contraintes de Reynolds locales ou l'homogénéisation de la PV, mais est plutôt une « empreinte » de la superrotation de la CZ, transmise vers le haut via la circulation méridienne et les contraintes visqueuses.

• Impact de la couche météorologique sur l'équilibre du vent thermique :
  La présence de la WL idéalisée modifie considérablement l'équilibre du vent thermique. Alors que les jets de la CZ restent largement alignés cylindriquement, la WL présente une circulation méridienne monocellulaire qui incline les isocontours de l'écoulement zonal par rapport à l'axe de rotation. Cela résulte en une grande baroclinicité à l'échelle dans la WL qui modifie l'équilibre du vent thermique, démontrant qu'une couche stable idéalisée peut même découpler l'alignement des écoulements profonds et peu profonds.

• Migration des jets et état stationnaire asymptotique :
  Une découverte critique est l'évolution à long terme de l'écoulement. Alors que la superrotation équatoriale s'équilibre relativement rapidement ($\mathcal{O}(1)$ temps de diffusion), les jets de haute latitude sont instables sur de longues échelles de temps. Sur $\mathcal{O}(10)$ temps de diffusion, les multiples jets de haute latitude migrent lentement vers les pôles et/ou fusionnent. L'état asymptotique final consiste en une structure en bandes « rapide-lent-rapide » : un fort jet superrotant à l'extérieur du cylindre tangent et seulement une paire de jets alternés par hémisphère à l'intérieur du cylindre tangent. Cela suggère que les multiples jets observés dans les simulations 3D antérieures, de durée plus courte, pourraient représenter des états transitoires plutôt que l'équilibre véritable.

Importance et affirmations
L'article affirme que le débat « conduite profonde versus conduite peu profonde » peut ne pas nécessiter une résolution unique, car les deux couches peuvent indépendamment générer de multiples jets via leurs mécanismes respectifs d'homogénéisation de la PV. Cependant, le couplage entre les couches est crucial pour déterminer l'alignement des écoulements zonaux. Les auteurs affirment que la superrotation équatoriale est fondamentalement un effet 3D entraîné par la géométrie des colonnes de Busse (évasement), qui ne peut être entièrement capturé par de simples arguments d'homogénéisation de la PV quasi 2D.

De plus, l'étude met en évidence une limitation potentielle dans les simulations 3D de pointe actuelles de Jupiter : les jets de haute latitude dans de nombreux modèles antérieurs étaient probablement encore en évolution. Les auteurs suggèrent que sans un mécanisme de dissipation physique (tel qu'un traînée indépendante de l'échelle) pour arrêter la cascade inverse d'énergie aux grandes échelles, les simulations 3D peuvent naturellement évoluer vers un état avec moins de jets de haute latitude que ceux actuellement observés sur Jupiter. L'article conclut que l'atteinte d'un état stationnaire avec de multiples jets de haute latitude dans les simulations 3D reste un défi, et que le véritable mécanisme soutenant le nombre de jets observé sur Jupiter peut impliquer une physique pas encore pleinement incorporée dans les modèles anélastiques globaux.