El Niño – Oscillation australe (ENSO) est le phénomène climatique le plus important affectant la météo à travers le monde. Ses phases, El Niño et La Niña, se caractérisent par des anomalies de température de surface de l'océan dans le Pacifique équatorial et ont un impact profond sur le climat global. Les observations montrent que ces phénomènes présentent un fort ancrage saisonnier : les pics surviennent généralement en décembre. Les explications précoces, comme la rétroaction positive de Bjerknes, décrivaient le mécanisme d'amplification des anomalies mais ne pouvaient expliquer leur alternance. Des modèles conceptuels, dont l'oscillateur à action retardée (DAO), ont été proposés pour formaliser l'interaction entre rétroactions positives et négatives. Ce travail se concentre sur le modèle DAO proposé par Suárez et Schopf, qui a récemment reçu une justification physique via la réduction de modèles spatialement distribués. Contrairement à des modèles plus complexes, il est suffisamment transparent pour une analyse mathématique détaillée, permettant d'étudier les phénomènes de résonance sous forçage périodique.
Les auteurs utilisent une combinaison de méthodes analytiques et numériques. Pour le modèle autonome de Suárez–Schopf (sans forçage externe), ils effectuent une analyse de stabilité linéaire des équilibres, obtenant des expressions analytiques pour les bifurcations fourche et de Hopf. Pour trouver les bifurcations non détectées par l'analyse linéaire (par exemple, les bifurcations nœud-col de cycles limites et les bifurcations homoclines), ils utilisent le logiciel de continuation DDE-BifTool. Pour le modèle périodiquement forcé, les auteurs développent une méthode de calcul du nombre de rotation sur les tores invariants attractifs. Cette méthode généralise une approche proposée antérieurement pour les équations différentielles ordinaires au cas de l'espace des phases infini-dimensionnel des systèmes à retard. L'idée clé consiste à itérer l'application stroboscopique jusqu'à ce que la trajectoire revienne dans un petit voisinage du point initial, puis le nombre de rotation est calculé comme le déplacement angulaire moyen. Pour surmonter les problèmes de projection dans l'espace infini-dimensionnel, deux projections différentes sur le plan sont utilisées, et on choisit celle qui donne une dépendance continue du nombre de rotation par rapport aux paramètres.
Pour le modèle autonome de Suárez–Schopf, un diagramme de bifurcation complet est construit dans le plan des paramètres (retard τ, force de rétroaction α). Il est montré que pour α > 1, le seul attracteur est un cycle limite symétrique issu d'une bifurcation de Hopf supercritique. Pour α < 1, on observe un tableau plus complexe : après la bifurcation fourche, deux équilibres symétriques apparaissent, qui perdent ensuite leur stabilité par une bifurcation de Hopf sous-critique ; entre les régions avec différents attracteurs se produisent une bifurcation nœud-col de cycles limites et une bifurcation homocline de collage. Pour le modèle périodiquement forcé (avec fréquence ω=2π), dans le plan (τ, c) (où c est l'amplitude du forçage), des régions avec un tore invariant attractif sont détectées. Les langues de résonance (régions de synchronisation de phase) s'avèrent très étroites et ne se chevauchent pas, ce qui indique une synchronisation faible. La structure des langues de résonance est organisée par les points critiques du graphe du nombre de rotation, à la fois sur les courbes de bifurcation du tore et à l'intérieur de la région d'existence du tore. Pour de grands retards, la structure de résonance se répète avec une période de 1 en τ, ce qui est décrit par l'application rationnelle Φ(p:q) = p:(q+2p). Pour α=0.75, on observe une bistabilité entre les orbites périodiques de période 1 et les tores invariants, incluant des bifurcations non classiques de « pli de tores » et de « collage de tores ».
Ce travail apporte une contribution significative à la compréhension de la dynamique des systèmes à retard périodiquement forcés, largement utilisés en modélisation climatique. La méthode développée pour calculer le nombre de rotation pour les systèmes à retard ouvre de nouvelles possibilités pour l'étude des phénomènes de résonance dans de tels systèmes. Le mécanisme découvert de répétition de la structure de résonance pour de grands retards, décrit par une transformation basée sur une matrice de SL(2,Z), est nouveau et pourrait être universel pour une large classe de modèles. Du point de vue climatologique, les résultats montrent que dans le cadre du modèle de Suárez–Schopf, le forçage saisonnier ne conduit qu'à une synchronisation de phase faible, ce qui est cohérent avec l'irrégularité observée de l'ENSO. Cependant, l'absence de chevauchement des langues de résonance signifie que le modèle ne reproduit pas la dynamique chaotique observée dans les données réelles, ce qui indique la nécessité de complexifier davantage le modèle.
Le sujet peut se développer dans plusieurs directions. Premièrement, il est intéressant d'étudier la dépendance de la structure de résonance par rapport à la fréquence de forçage ω, car des changements de ω peuvent entraîner des chevauchements de langues de résonance et une dynamique chaotique. Deuxièmement, on peut envisager une rupture de la symétrie Z2 du modèle pour tenir compte de l'asymétrie entre El Niño et La Niña. Troisièmement, l'inclusion d'un forçage stochastique pourrait renforcer l'irrégularité et rapprocher le modèle des données réelles. Enfin, l'étude d'autres modèles conceptuels de l'ENSO, comme le modèle de Gil–Zaliapin–Thompson, à l'aide des méthodes développées pourrait éclairer les schémas généraux des phénomènes de résonance dans les oscillateurs climatiques.
Les résultats pourraient influencer le développement de modèles conceptuels plus réalistes de l'ENSO, utilisés pour la prévision des anomalies climatiques. La méthode de calcul du nombre de rotation est applicable à un large éventail de problèmes avec forçage périodique en physique, biologie et ingénierie.
La prochaine étape immédiate est une étude systématique de l'influence de la fréquence de forçage ω sur la structure de résonance et l'apparition du chaos. Il est également important de comparer avec d'autres modèles de l'ENSO et avec les données d'observation pour évaluer l'adéquation du modèle.
Ce travail est directement lié au problème fondamental de la compréhension des mécanismes de l'irrégularité de l'ENSO. Les modèles conceptuels existants, y compris celui étudié, ne reproduisent pas le niveau de chaos observé, ce qui indique la nécessité de prendre en compte des facteurs supplémentaires tels que l'asymétrie des phases, le forçage stochastique ou des non-linéarités plus complexes. L'étude des phénomènes de résonance dans de tels modèles aide à identifier les conditions dans lesquelles une transition vers le chaos est possible, et ainsi à se rapprocher de l'explication de la variabilité réelle du système climatique.
Cette section s'adresse à l'auteur du travail et à ses collègues. Les travaux voisins ont été trouvés par une recherche sémantique dans notre archive, et non par correspondance de mots. Nous n'évaluons pas : nous montrons le champ autour.
Ce que fait ce travail
L'article porte sur la structure de résonance de l'équation de Suarez–Schopf avec délai sous forçage périodique. Les auteurs ont construit le diagramme de bifurcation du modèle autonome, développé une méthode de calcul du nombre de rotation pour les systèmes à délai et mis en évidence la répétition des langues de résonance pour de grands délais.
Où est sa force
Le point fort réside dans le développement et l'application d'une méthode de calcul du nombre de rotation sur des tores invariants pour des systèmes infini-dimensionnels à délai, ce qui permet d'étudier en détail les phénomènes de résonance.
Vers où avancer
À essayer
Pourquoi cela peut compter
Comprendre les phénomènes de résonance dans les modèles climatiques pourrait améliorer la prévision des phases d'El Niño et de leur calage saisonnier, ce qui est important pour l'évaluation des risques climatiques mondiaux.
Travaux voisins dans notre archive
L'archive est vide autour de ce travail : l'analyse la plus proche n'est qu'à 0.49, contre 0,64 d'ordinaire. Il semble qu'on ait peu foré à côté.
Le nombre de rotation sur un tore invariant peut être rationnel (la dynamique est alors périodique) ou irrationnel (la trajectoire remplit alors le tore de manière dense). Dans ce travail, le calcul du nombre de rotation a nécessité environ 768 000 évaluations, prenant environ 8 heures sur un ordinateur Apple M1.