El Niño : une prédiction à 20 mois grâce à un réseau de neurones explicable ✛

Original: "Towards Long-Range ENSO Prediction with an Explainable Deep Learning Model"
arXiv:2503.19502 · 2025-03-25 · CC BY 4.0 · 6 min
Un hybride de réseau convolutif et de transformeur prédit El Niño à 20 mois et explique quels signaux de l'océan et de l'atmosphère déclenchent le phénomène.

Contexte

La prévision d'El Niño est cruciale pour l'agriculture, l'énergie et la gestion des risques d'inondations et de sécheresses dans le monde entier. Les modèles dynamiques classiques ne prédisent le phénomène avec confiance que 6 à 12 mois à l'avance, et la barrière de prévisibilité printanière réduit fortement la précision des prévisions faites au printemps. Les réseaux de neurones profonds ont déjà repoussé l'horizon à 17–18 mois, mais la plupart ne s'appuient que sur la température de surface de l'océan et ignorent les processus couplés océan-atmosphère. Ces travaux s'inscrivent dans les tentatives récentes d'utiliser des architectures hybrides et des données multidimensionnelles pour la prévision à long terme : des études voisines ont montré que les réseaux convolutifs capturent les motifs spatiaux et que les transformeurs saisissent les dépendances à longue portée, mais la combinaison de ces approches avec un ensemble complet de variables physiques restait un défi ouvert.

Méthodes

CTEFNet utilise une architecture hybride : trois blocs de couches convolutives extraient les caractéristiques spatiales locales, puis un transformeur avec encodeur et décodeur modélise les dépendances temporelles. Les données d'entrée sont 12 mois consécutifs de neuf variables : température de surface de l'océan, contenu thermique, profondeur de la couche de mélange, salinité, pression au niveau de la mer, composantes zonale et méridienne de la vitesse des courants et de la contrainte du vent. Tous les champs sont interpolés sur une grille de 1° × 2° et normalisés. Le modèle a été entraîné sur 18 modèles climatiques du projet CMIP6 selon le scénario SSP370 (2015–2100), validé sur les réanalyses ORAS5 et ERA5 (1958–1978), et testé sur GODAS et ERA5 (1980–2021). La variable cible était l'indice Niño 3.4 sur les 24 mois suivants. Pour l'interprétation, une analyse de sensibilité par gradient a été utilisée : la rétropropagation de l'erreur calculait la contribution de chaque point d'entrée à la prévision de l'indice, et la significativité statistique était sélectionnée par un test t au niveau de 95 %.

Résultats

CTEFNet maintient une corrélation de prévision supérieure à 0.7 pour toutes les saisons jusqu'à 12 mois, et supérieure à 0.6 jusqu'à 17 mois, tandis que l'horizon effectif (corrélation supérieure à 0.5) atteint 20 mois. Cela représente 3 à 5 mois de plus que les modèles dynamiques NMME, et 3 à 4 mois de plus que les meilleurs réseaux profonds (CNN, Geoformer, ResCNN, ResoNet, STPNet). L'analyse saisonnière montre que le modèle conserve une corrélation supérieure à 0.5 pendant plus de 20 mois pour les prévisions faites de juin à décembre, et 16 mois pour celles faites au printemps, ce qui indique un affaiblissement de la barrière printanière. L'analyse par gradient pour 11 événements El Niño a révélé une séquence de précurseurs : 11 mois avant le pic, une sensibilité positive du contenu thermique dans l'ouest du Pacifique indique l'accumulation d'eau chaude, et une sensibilité négative de la profondeur de la couche de mélange signale la suppression du mélange vertical. À 8 mois, les signaux se propagent vers le centre et l'est du Pacifique et l'est de l'océan Indien, ce qui correspond au déclenchement de la rétroaction de Bjerknes. Simultanément, le refroidissement de l'Atlantique tropical nord, via des ponts atmosphériques, affaiblit les alizés et favorise le développement d'El Niño. En octobre, les vents d'ouest et le transport d'eau chaude vers l'est s'intensifient, et la sensibilité négative de la température de surface dans l'océan Indien déclenche la réponse de Matsuno–Gill, qui alimente encore le phénomène.

Implications

Ce travail se distingue des études voisines par le fait qu'il combine pour la première fois des couches convolutives et un transformeur dans un encodeur-décodeur avec un ensemble complet de neuf variables océaniques et atmosphériques, et non seulement la température de surface. Cela permet au modèle de capturer les interactions inter-bassins couplées et d'atteindre un horizon de prévision de 20 mois, soit 2 à 3 mois de plus que les approches hybrides précédentes. De plus, l'analyse de sensibilité par gradient est utilisée ici comme substitut aux modèles adjoints, qui sont coûteux et linéaires : la rétropropagation fournit une estimation non linéaire et peu coûteuse en calcul de la contribution de chaque variable à chaque instant. Cela fait de l'interprétation non pas une simple démonstration, mais un outil de travail pour vérifier les mécanismes physiques.

Développements futurs

Le modèle peut être étendu pour distinguer les types d'El Niño : canonique (Pacifique est) et Modoki (Pacifique central), qui affectent différemment le climat mondial. Une autre direction est l'atténuation supplémentaire de la barrière de prévisibilité printanière, ce qui nécessiterait une modélisation plus fine des décalages saisonniers des vents et des courants, éventuellement avec des données à plus haute résolution ou des contraintes physiques supplémentaires dans la fonction de perte. Il serait également prometteur d'utiliser l'analyse par gradient proposée pour étudier l'asymétrie entre El Niño et La Niña, qui est déjà apparue dans les calculs préliminaires.

Impact

Une meilleure prévision d'El Niño affectera directement la planification saisonnière dans l'agriculture, la gestion des ressources en eau, l'assurance contre les catastrophes naturelles et la logistique énergétique dans les régions dépendantes des précipitations et de la température. La méthode de sensibilité par gradient pourrait être transférée à d'autres modes climatiques, comme le dipôle de l'océan Indien ou la circulation méridienne de retournement atlantique.

Prochaines étapes

La prochaine étape est l'adaptation de CTEFNet à la prévision d'indices distinguant les types d'El Niño, et la validation sur des données indépendantes au-delà de 1980–2021. Il conviendrait également de comparer la sensibilité du modèle avec les calculs directs de modèles adjoints pour des événements individuels, afin d'évaluer quantitativement la précision de l'approche par gradient.

Problèmes ouverts clés

En climatologie, la question des causes de la barrière de prévisibilité printanière reste ouverte : pourquoi les prévisions faites en mars–mai sont-elles systématiquement moins bonnes que celles des autres saisons ? CTEFNet atténue cette barrière sans l'éliminer, et son analyse par gradient pourrait indiquer quels processus deviennent instables pendant cette période. Un autre problème ouvert est le rôle relatif des trois bassins océaniques dans le déclenchement d'El Niño : la sensibilité du modèle à l'Atlantique et à l'océan Indien fournit une base quantitative pour tester les hypothèses sur les téléconnexions inter-bassins.

Le regard de la machine du savoir

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

Cet article présente CTEFNet, un réseau de neurones hybride combinant des couches convolutionnelles et un transformateur pour prévoir l'indice Niño 3.4 sur 24 mois. Le modèle a été entraîné sur l'ensemble CMIP6, testé sur les réanalyses GODAS et ERA5, et a montré une corrélation supérieure à 0.5 jusqu'à 20 mois, atténuant ainsi la barrière du printemps. L'analyse par gradient a révélé des précurseurs physiquement interprétables, tels que l'accumulation de chaleur dans l'ouest de l'océan Pacifique.

Où est sa force

La force de ce travail réside dans le fait qu'il ne se contente pas d'étendre l'horizon de prévision, mais utilise également un ensemble complet de variables physiques de l'océan et de l'atmosphère. De plus, l'interprétabilité via l'analyse par gradient relie la prévision à des mécanismes connus, ouvrant ainsi la voie à une analyse physique plus approfondie.

Vers où avancer

À essayer

  • Vérifier si le gain en horizon de prévision se maintient lorsque le modèle est entraîné uniquement sur la température de surface de l'océan et le contenu thermique, comme dans certains travaux voisins.
  • Utiliser la technique de lecture inverse des séries temporelles pour renforcer la sensibilité aux processus à évolution lente, tels que l'accumulation de chaleur dans l'ouest de l'océan Pacifique.
  • Essayer de remplacer une partie des blocs convolutionnels par des modules plus légers afin de réduire les coûts de calcul et de faciliter l'application du modèle dans les systèmes opérationnels.

Pourquoi cela peut compter

Si l'on parvient à maintenir la précision tout en simplifiant les données d'entrée ou l'architecture, cela rendrait la prévision à long terme d'El Niño plus accessible pour une utilisation opérationnelle, ce qui est important pour l'alerte précoce aux risques climatiques.

El Niño de 1997–1998 a été si fort que la température mondiale de surface de l'océan a augmenté de 0.2 °C par rapport à la moyenne décennale, et le niveau de la mer dans l'est du Pacifique a temporairement monté de 20 centimètres en raison de la dilatation thermique de l'eau.

Dans le roman de Kim Stanley Robinson « Quarante signes de pluie », les climatologues luttent contre les conséquences d'un changement soudain des courants océaniques provoqué par le réchauffement ; une prévision précise à long terme d'El Niño pourrait être l'outil que les héros utilisent pour prévenir la catastrophe.

Loss = \frac{1}{T_{out}} \sum_{t=1}^{T_{out}} \sqrt{(Nino3.4^{out}_t - Nino3.4^{tg}_t)^2}
Le modèle minimise l'écart entre l'indice Niño 3.4 prédit et l'indice réel sur un horizon de 24 mois.
Grad_t = \frac{\partial Nino3.4_t}{\partial Inputs_t}
Indique à quel point un changement de chaque variable d'entrée à un mois donné affecte la prévision de l'indice Niño 3.4.
AvgGrad = mean_t (Grad_t)
La sensibilité finale de chaque point d'entrée, moyennée sur tous les mois de prévision, est utilisée pour identifier les précurseurs.

Chiffres clés

  • horizon de prévision effectif: 20 mois
  • corrélation supérieure à 0.7: jusqu'à 12 mois
  • corrélation supérieure à 0.6: jusqu'à 17 mois
  • nombre de variables d'entrée: 9
  • nombre de modèles climatiques pour l'entraînement: 18
Notions
ENSO El Niño-Southern Oscillation spring predictability barrier Multivariate ENSO Index Jin-Timmermann model El Niño index equatorial wave propagation TAO/TRITON moorings La Niña ERSST Southern Oscillation NOAA PMEL / TAO NOAA CPC zebiak cane model North Equatorial Countercurrent prévision interannual variability Indian Ocean capacitor effect régression par réseau de neurones Warm water volume