Effects of global core-mantle boundary topography on outer-core convection and topographic torques

Cette étude utilise des simulations numériques directes pour démontrer que la topographie globale de la limite noyau-manteau génère une nouvelle instabilité et renforce considérablement la convection du noyau externe ainsi que les couples topographiques, offrant ainsi un mécanisme cohérent avec les variations décennales observées de la durée du jour.

L'essentiel

Imaginez la Terre comme une gigantesque bille en rotation. À l'intérieur, un noyau de fer liquide tourbillonne et bouillonne, générant le champ magnétique de notre planète. Autour de ce noyau liquide se trouve le manteau solide, qui agit comme une coquille épaisse et rocheuse.

Habituellement, les scientifiques imaginent la frontière entre ce noyau liquide et cette coquille rocheuse (appelée la limite noyau-manteau, ou LNM) comme une sphère parfaitement lisse et ronde. Mais cet article soutient que cette frontière est en réalité bosselée et irrégulière, bien plus comme la surface d'une pomme de terre que comme celle d'une bille de billard. Ces bosses sont causées par d'énormes structures situées profondément dans le manteau rocheux, dont certaines mesurent des milliers de kilomètres de large.

Les chercheurs ont utilisé de puissants supercalculateurs pour simuler ce qui se produit lorsque ce « noyau liquide » tourbillonne contre une « coquille bosselée ». Voici ce qu'ils ont découvert, expliqué simplement :

1. La « Glissade Lisse » contre la « Piste Bosselée »

Dans une sphère parfaitement lisse, le liquide à l'intérieur souhaite s'écouler en anneaux circulaires réguliers autour de l'axe de rotation de la Terre. C'est comme un patineur tournant sur une patinoire parfaitement lisse ; il peut glisser sans effort en cercle.

Cependant, lorsque la frontière est bosselée, c'est comme placer une série de dos d'âne ou de collines sur cette patinoire. L'écoulement du liquide est forcé de changer de direction pour passer par-dessus ou autour de ces bosses. Les chercheurs ont découvert que ces bosses aident en réalité le liquide à se déplacer plus vite et à transporter la chaleur plus efficacement. C'est comme si les bosses agissaient comme un catalyseur, donnant au liquide une « poussée » qu'il n'obtiendrait pas sur une surface lisse. Dans leurs simulations, ces bosses ont augmenté la vitesse de l'écoulement et la quantité de chaleur se déplaçant du centre vers la périphérie jusqu'à 100 %.

2. La « Nouvelle Instabilité » (La Surprise Subcritique)

Il existe une règle en physique qui stipule que la convection liquide (comme l'eau qui bout) ne commence que lorsque la chaleur est suffisamment intense pour surmonter la résistance du fluide. Les chercheurs ont découvert quelque chose de surprenant : les bosses sur la frontière peuvent briser cette règle.

Même lorsque le noyau n'est pas assez chaud pour commencer à bouger de lui-même, les bosses peuvent créer un nouveau type d'instabilité qui met le liquide en mouvement de toute façon. Imaginez une balle reposant au fond d'une vallée profonde ; normalement, elle a besoin d'une grande poussée pour en sortir. Mais si la vallée a une forme étrange et bosselée, une toute petite pichenette peut suffire à faire rouler la balle. Cela signifie que le noyau de la Terre pourrait bouillonner et générer son champ magnétique même lorsqu'il est « plus frais » que nous ne le pensions auparavant.

3. Le « Couple » (La Toupie Qui Chancelle)

La Terre tourne comme une toupie. Parfois, la durée de notre journée change de fractions infimes de seconde (millisecondes) sur des périodes de 6 à 60 ans. Les scientifiques soupçonnent depuis longtemps que l'interaction entre le noyau liquide en rotation et le manteau solide est responsable de ces micro-chancellements.

Les chercheurs ont calculé le « couple » (la force de torsion) que le noyau liquide exerce sur la frontière bosselée. Ils ont découvert que les bosses créent une force de torsion significative.

• L'Analogie : Imaginez pousser un manège en rotation. Si vous le poussez sur un bord lisse, il est difficile de modifier sa vitesse. Mais si vous poussez contre un bord bosselé et irrégulier, vous pouvez vous accrocher aux bosses et faire tourner l'ensemble beaucoup plus efficacement.
• Le Résultat : Leurs calculs montrent que la force de torsion générée par ces bosses est suffisamment forte pour expliquer les changements observés dans la durée de notre journée.

4. L'Effet de « Verrouillage »

L'une des découvertes les plus intéressantes concerne la façon dont l'écoulement liquide interagit avec des formes spécifiques de bosses.

• L'Analogie : Imaginez un danseur essayant de bouger sur une musique. Si la musique (l'écoulement) et le motif du sol de danse (les bosses) correspondent parfaitement, le danseur pourrait se « verrouiller » dans un rythme spécifique.
• Le Résultat : Lorsque les bosses avaient une taille et une forme spécifiques correspondant au rythme naturel de l'écoulement liquide, celui-ci se « verrouillait » sur les bosses. Bien que cela rende l'écoulement très organisé, cela a en réalité réduit la force de torsion (couple) parce que le liquide ne se battait plus contre les bosses ; il se contentait de les suivre. Cela suggère que la forme des bosses compte tout autant que leur taille.

Résumé

Cet article utilise des modèles informatiques pour montrer que la frontière « bosselée » entre le noyau liquide et le manteau solide de la Terre n'est pas un simple mur passif. C'est un participant actif qui :

1. Accélère l'écoulement liquide et le transfert de chaleur.
2. Déclenche l'écoulement même lorsqu'il est trop froid pour bouger seul.
3. Tord la rotation de la Terre, expliquant pourquoi nos journées deviennent légèrement plus longues ou plus courtes au fil des décennies.

L'étude confirme que pour comprendre comment fonctionne le champ magnétique de la Terre et pourquoi la durée de nos journées change, nous ne pouvons pas traiter le noyau comme une sphère lisse et parfaite ; nous devons prendre en compte la réalité rugueuse et bosselée de la frontière.

Résumé technique

Résumé technique : Effets de la topographie globale de la limite noyau-manteau sur la convection du noyau externe et les couples topographiques

Énoncé du problème
La frontière entre le noyau externe liquide de la Terre et le manteau solide (la limite noyau-manteau, ou CMB) n'est pas parfaitement sphérique ou lisse. Les études sismiques et géodynamiques suggèrent la présence d'hétérogénéités topographiques, potentiellement liées aux provinces à faible vitesse de cisaillement de grande taille (LLSVP) dans le manteau profond. Bien qu'il soit hypothétisé que la topographie de la CMB couple le noyau et le manteau — potentiellement en entraînant les fluctuations décennales et sous-décennales observées dans la longueur du jour ($\Delta$LOD) et en influençant le champ géomagnétique — les mécanismes dynamiques spécifiques restent mal compris. Les investigations précédentes étaient limitées aux régimes linéaires, à des modèles bidimensionnels simplifiés ou à des caractéristiques topographiques isolées, échouant à capturer les effets d'une topographie globalement distribuée et spectralement large sur la convection pleinement turbulente.

Méthodologie
Les auteurs ont employé des simulations numériques directes (DNS) de convection dans une coquille en rotation pour étudier ces effets dynamiques. L'étude a utilisé le code spectral éléments Nek5000 pour résoudre les équations de Boussinesq adimensionnelles pour un fluide dans une coquille sphérique avec un rapport d'aspect $\chi = 0.35$.

• Géométrie : La frontière intérieure était sphérique, tandis que la frontière extérieure était déformée selon $r_{cmb} = r_o - \epsilon h(\theta, \phi)$, où $\epsilon$ est l'amplitude topographique adimensionnelle et $h$ la fonction de forme.
• Paramètres : Les simulations ont été menées à un nombre d'Ekman $Ek = 10^{-6}$ (et un sous-ensemble à $10^{-4}$) et un nombre de Rayleigh réduit $\widetilde{Ra} = 40$ pour le régime turbulent, avec un nombre de Prandtl $Pr = 1$.
• Modèles topographiques : Cinq formes topographiques distinctes ont été testées :
  1. Harmoniques sphériques uniques : $h_{2,1}$, $h_{8,4}$, $h_{32,16}$, et $h_{32,17}$.
  2. Un « modèle tomographique » ($h_t$) dérivé de cartes d'analyse de clusters sismiques des LLSVP, lissé pour représenter des intrusions continues.
• Analyse : L'étude a examiné le transport de chaleur (nombre de Nusselt, $Nu$), le transport de quantité de mouvement (nombre de Reynolds, $Re$), la morphologie de l'écoulement (spectres d'énergie cinétique, contours géostrophiques) et les couples topographiques axiaux résultants ($\Gamma_z$).

Contributions et résultats clés

1. Confirmation de l'instabilité topographique en géométrie globale :
   L'étude confirme l'existence d'une instabilité linéaire entraînée par la topographie à des nombres de Rayleigh inférieurs au seuil critique pour une coquille sphérique ($Ra < Ra_c$), un phénomène prédit précédemment par Bell & Soward (1996) dans un modèle d'anneau 2D. Bien que le mécanisme d'instabilité (la flottabilité entraînant l'écoulement le long des contours géostrophiques déformés) soit confirmé, l'échelle quantitative des vitesses d'écoulement diffère de la théorie de l'anneau en raison des différences géométriques dans le profil de température de fond et de la nature des contours géostrophiques dans une sphère.

2. Amélioration du transport par déformation des contours géostrophiques :
   La forme des contours géostrophiques (lignes de hauteur axiale constante) est centrale pour la dynamique. Dans une sphère parfaite, ce sont des cercles, et la flottabilité ne peut effectuer de travail sur l'écoulement géostrophique. La topographie déforme ces contours, permettant à la flottabilité d'effectuer un travail sur l'écoulement moyenné dans le temps.

   • Augmentation du transport : Pour les topographies présentant une symétrie équatoriale significative (par exemple, $h_{8,4}$, $h_{32,16}$), la déformation permet à la flottabilité d'entraîner directement l'écoulement géostrophique. Cela se traduit par des augmentations des nombres de Reynolds et de Nusselt allant jusqu'à $\sim 100\%$ par rapport au cas de contrôle sphérique.
   • Dépendance à la symétrie : Les topographies avec $\ell-m$ impair (par exemple, $h_{2,1}$, $h_{32,17}$) présentent une antisymétrie équatoriale, conduisant à des déformations d'amplitude plus faible des contours géostrophiques ($O(\epsilon^2)$) par rapport aux cas $\ell-m$ pairs ($O(\epsilon)$). Par conséquent, l'amélioration de l'écoulement géostrophique est supprimée dans les cas antisymétriques, bien que le transport puisse encore augmenter via une résonance avec les échelles convectives.

3. Échelle du couple topographique et verrouillage de l'écoulement :
   L'étude examine le couple axial exercé par le noyau sur la topographie de la CMB.

   • Loi d'échelle : Les résultats confirment et étendent la loi d'échelle $\Gamma'_z \sim \epsilon Re^2$, où le couple varie linéairement avec l'amplitude topographique et quadratiquement avec la vitesse d'écoulement.
   • Verrouillage de l'écoulement : Pour le cas $h_{32,16}$, où la longueur d'onde topographique est proche de la longueur d'onde convective critique, l'écoulement présente un « verrouillage » sur le motif topographique à fortes amplitudes. Cela entraîne une diminution de la composante de pression hors phase (qui génère le couple), conduisant à une saturation de l'amplitude du couple malgré l'augmentation des vitesses d'écoulement.
   • Théorie généralisée : Les auteurs proposent une procédure pour estimer les couples pour une topographie spectralement large (comme le modèle tomographique) en convoluant le spectre topographique (supposé suivre la règle de Kaula, $\eta_{\ell,m} \sim \ell^{-2}$) avec le spectre de pression fluctuante. Cela donne des estimations de couple cohérentes avec l'ampleur requise pour entraîner les variations observées de $\Delta$LOD.

4. Couples visqueux vs couples topographiques :
   Le rapport des fluctuations de couple visqueux sur le couple topographique ($\gamma$) a été trouvé inférieur à l'unité pour la plupart des paramètres, soutenant la conclusion que les couples topographiques, plutôt que le couplage visqueux, sont le mécanisme dominant pour les variations de $\Delta$LOD à l'échelle décennale.

Importance et affirmations
L'article revendique quatre résultats principaux :

1. Il valide le mécanisme d'instabilité de Bell & Soward (1996) dans une géométrie sphérique globale entièrement tridimensionnelle.
2. Il démontre que la topographie globale déforme les contours géostrophiques, permettant à la flottabilité d'effectuer un travail sur l'écoulement géostrophique, améliorant ainsi considérablement le transport de chaleur et de quantité de mouvement (jusqu'à des augmentations de 100 %).
3. Il montre que la topographie de grande amplitude peut remplacer le paradigme canonique visqueux $Ek^{1/3}$ pour le début du mouvement et l'écoulement géostrophique moyenné dans le temps, bien que la turbulence à petite échelle conserve son empreinte visqueuse.
4. Il établit la généralité de l'échelle de couple dynamique $\Gamma'_z \sim \epsilon Re^2$ à travers diverses formes topographiques globales. En développant une théorie d'échelle généralisée pour la topographie spectralement large, les auteurs estiment que les couples extrapolés atteignent des amplitudes de $\sim 10^{18}$ N m, ce qui est cohérent avec le couple requis pour entraîner les variations observées de $\Delta$LOD décennales et sous-décennales.

Les auteurs notent que, bien que leurs résultats fournissent un cadre dynamique robuste, la structure spécifique de la topographie de la CMB reste mal contrainte par la sismologie. Un travail futur est identifié comme nécessaire pour intégrer les champs magnétiques (magnétohydrodynamique), affiner les modèles topographiques avec des profils plus lisses, et explorer la transition vers des régimes où les échelles topographiques remplacent complètement les échelles visqueuses dans les conditions du noyau terrestre.