Un indice de mousson venu de la mer : comment l'apprentissage par renforcement améliore les prévisions de pluie en Thaïlande ✛

Original: "Reinforcement Learning to Discover a North-East Monsoon Index for Rainfall Prediction in Thailand"
· Kiattikun Chobtham
arXiv:2601.10181 · 2026-01-15 · CC BY 4.0 · 3 min
Un agent neuronal identifie des zones océaniques dont la température prédit les pluies en Thaïlande des mois à l'avance.

Contexte

Ce travail se distingue parmi les études qui relient les schémas climatiques à grande échelle aux précipitations régionales. Les articles voisins utilisent des indices globaux et l'apprentissage automatique pour prévoir les sécheresses, El Niño, les éruptions solaires et la météo spatiale. Contrairement à eux, cette étude se concentre sur un indice de mousson local, spécifique à la Thaïlande, et applique l'apprentissage par renforcement pour optimiser les zones géographiques, ce qui n'apparaît pas dans les travaux voisins.

Méthodes

L'auteur a utilisé les données de température de surface de la mer (OISST) de 1982 à 2024 et les précipitations de 74 stations du département météorologique thaïlandais (pour l'entraînement) et de 384 stations de l'Institut d'hydro-informatique (pour les tests). Les stations ont été regroupées en 12 clusters par classification hiérarchique avec ACP. Ensuite, un agent Deep Q-Network avec un espace d'actions discret (déplacement et modification de la taille des rectangles) a optimisé deux zones océaniques A et B, en maximisant une fonction objectif basée sur le carré de la corrélation entre l'indice (différence de SST) et les précipitations pendant les périodes d'installation (octobre à mars) et de retrait (avril à septembre) de la mousson du nord-est. La fraction océanique dans chaque zone devait être d'au moins 80 %. Les meilleures zones ont ensuite été utilisées pour calculer l'indice, qui a été ajouté à un modèle LSTM avec une fenêtre d'entrée de 24 mois et un horizon de prévision de 12 mois.

Résultats

L'indice optimisé a atteint une fonction objectif Q = 0.497 (contre 0.052 pour les zones d'origine), avec des corrélations allant jusqu'à -0.720 au début de la mousson dans le sud et jusqu'à -0.738 pendant le retrait dans le nord. En ajoutant l'indice au LSTM, l'erreur RMSE a diminué dans la plupart des clusters : par exemple, dans le cluster 1, elle est passée de 99.79 à 94.54 mm/mois, et dans le cluster 12, de 59.24 à 52.38 mm/mois. Seul le cluster 8 n'a montré aucune amélioration (59.81 contre 59.91).

Implications

Ce travail se distingue des autres en introduisant un indice local, optimisé spécifiquement pour la Thaïlande, plutôt que de s'appuyer sur des téléconnexions globales. C'est la première tentative d'utiliser l'apprentissage par renforcement pour trouver les zones optimales de calcul d'un indice de mousson. Le résultat montre que les schémas locaux peuvent améliorer considérablement les prévisions, même lorsque les indices globaux sont déjà inclus.

Développements futurs

À l'avenir, on pourrait utiliser des modèles spatio-temporels tels que Convolution-LSTM, les réseaux de neurones à graphes ou les opérateurs neuronaux pour tenir compte de la structure spatiale des précipitations et peut-être améliorer les prévisions. On pourrait également étendre l'approche à d'autres régions à climat de mousson.

Impact

Ce travail pourrait influencer la planification hydrologique, la gestion des réservoirs et la planification agricole en Thaïlande, où une prévision précise à long terme des précipitations est cruciale.

Prochaines étapes

La prochaine étape consistera à appliquer des modèles spatio-temporels avec le nouvel indice et à les tester sur des données indépendantes. Il faudrait également étudier la robustesse de l'indice face au changement climatique.

Problèmes ouverts clés

Ce travail est directement lié au problème de la représentation des processus climatiques locaux dans les modèles de prévision, qui reste ouvert en hydrométéorologie. Il aborde également la question de la détection automatique de prédicteurs significatifs dans de grands ensembles de données.

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

L'auteur propose un indice local de la mousson du nord-est de la Thaïlande, calculé comme la différence de SST entre deux régions optimisées, identifiées par un agent d'apprentissage profond par renforcement (Deep Q-Network). Cet indice est intégré à un modèle LSTM pour prévoir les précipitations à 12 mois sur 12 groupes de stations. L'indice optimisé atteint la fonction objectif Q=0,497, et son ajout au modèle réduit l'erreur RMSE dans la plupart des groupes, en particulier dans le sud de la Thaïlande.

Où est sa force

La force de ce travail réside dans l'introduction d'un indice local spécifique à la Thaïlande et dans la première application de l'apprentissage par renforcement pour déterminer les régions géographiques optimales de son calcul. Cela permet de prendre en compte les circulations de mousson locales qui pourraient être négligées si l'on n'utilisait que des indices globaux.

Vers où avancer

  • Il pourrait être intéressant d'examiner comment l'indice proposé interagit avec les indices globaux, comme El Niño, dans un modèle de prévision unifié. Dans les archives, une analyse voisine prédit El Niño avec un an et demi d'avance ; l'utilisation conjointe d'indices locaux et globaux pourrait améliorer les prévisions à différents horizons.appui : El Niño se laisse prédire dix-huit mois à l'avance arXiv:2601.12251
  • Une autre piste concerne la prévision de la sécheresse plusieurs mois à l'avance grâce à un réseau de neurones qui lit le temps à l'envers. Une approche similaire, avec une lecture rétrospective des séries temporelles de SST et de précipitations, pourrait améliorer la prévision ou révéler des prédicteurs supplémentaires pour l'indice de mousson.appui : Comment un réseau de neurones, en écoutant la météo à l'envers, annonce la sécheresse des mois à l'avance arXiv:2607.00331
  • Dans les archives, on trouve aussi un travail où la pluie cosmique (rayons cosmiques) est utilisée pour améliorer les prévisions d'ouragans. On pourrait vérifier si le flux de rayons cosmiques apporte des informations supplémentaires à l'indice de mousson pour la prévision des précipitations en Thaïlande, en particulier lors d'événements extrêmes.appui : Les muons, un thermomètre pour les cyclones arXiv:2509.04627

À essayer

  • Envisager d'inclure dans le modèle LSTM non seulement l'indice, mais aussi ses valeurs décalées ou des transformations non linéaires (par exemple, le carré ou la valeur absolue) pour tenir compte d'une éventuelle relation non linéaire avec les précipitations.
  • Tester la robustesse des régions identifiées en modifiant la période d'apprentissage, par exemple en entraînant l'agent sur la première moitié des données et en le validant sur la seconde, afin d'évaluer la stabilité de l'optimisation.
  • Comparer l'indice proposé avec d'autres indices locaux, s'ils existent, ou avec des alternatives simples (par exemple, la différence de SST dans des régions fixes choisies par des experts) pour démontrer l'avantage de l'optimisation.

Pourquoi cela peut compter

Si l'indice local optimisé par apprentissage par renforcement améliore effectivement les prévisions de précipitations, il pourrait servir de modèle pour d'autres régions à climat de mousson. Une telle approche pourrait contribuer à la mise en place de systèmes d'alerte précoce plus précis pour les sécheresses et les inondations, ce qui est crucial pour l'agriculture et la gestion des ressources en eau.

L'agent d'apprentissage par renforcement a exploré des milliers de combinaisons de rectangles dans l'océan pour trouver la paire de zones dont la différence de température prédit le mieux les pluies. Les zones optimales se sont avérées décalées de centaines de kilomètres par rapport à celles proposées dans les études précédentes.

Z(t) = \text{normalise}(SST_B(t) - SST_A(t))
Différence normalisée des températures moyennes de surface de la mer dans les zones B et A à l'instant t
Q(A,B) = \frac{1}{2}(R^2_{NEM onset} + R^2_{NEM retreat})
Moyenne des carrés des corrélations entre l'indice et les précipitations pendant les périodes d'installation et de retrait de la mousson

Chiffres clés

  • fonction objectif Q: 0.497 (contre 0.052 pour les zones d'origine)
  • corrélation maximale: -0.738 (cluster 5, retrait de la mousson)
  • réduction du RMSE dans le cluster 1: de 99.79 à 94.54 mm/mois
  • nombre de clusters de stations: 12
  • horizon de prévision: 12 mois
Notions
monsoon monsoon onset NOAA CPC ERA5 Oceanic Niño index forêt aléatoire IRI apprentissage par renforcement Multivariate ENSO Index interannual variability forecast skill relative ONI GPCP rmse ERSST precipitation anomaly Niño 4 apprentissage automatique probabilistic forecast tropospheric biennial oscillation