Les El Niño extrêmes, comme celui de 1997-1998, provoquent sécheresses, inondations et mortalité des récifs coralliens dans le monde entier. Selon les estimations, leur fréquence pourrait doubler au XXIe siècle, mais les détails précis de la distribution des événements les plus forts restent flous, car la simulation directe nécessite des centaines de milliers d'années de simulation. Ce travail s'inscrit dans le cadre des recherches sur l'échantillonnage des événements rares dans les modèles climatiques : des méthodes de Monte Carlo par diffusion ont déjà été appliquées aux saisons chaudes et aux cyclones tropicaux, ainsi que le splitting adaptatif multi-niveaux pour des modèles simples. Cet article teste pour la première fois systématiquement l'algorithme TEAMS sur un modèle intermédiaire d'El Niño, où une comparaison avec une longue référence est possible.
Nous avons utilisé le modèle Zebiak-Cane, un modèle intermédiaire de l'ENSO avec un bloc atmosphérique réagissant au chauffage et un bloc océanique en eaux peu profondes. Nous avons lancé un ensemble de 16 ancêtres de 25 ans chacun, sélectionnés à partir d'une courte simulation directe. À chaque étape, nous avons sélectionné 25 % des membres ayant le plus petit maximum de l'indice NINO3, et pour les autres, nous avons généré des descendants avec une petite perturbation aléatoire de la température de surface de l'océan (amplitude 0,1 °C) 7,5 ans avant le franchissement du seuil. Nous avons répété 50 exécutions indépendantes avec un budget de 2000 années-modèle chacune. Le temps de retour a été estimé par la méthode des maxima par blocs avec une approximation de Poisson, et l'incertitude par bootstrap sur 5000 échantillons.
TEAMS a reproduit les estimations du temps de retour de la simulation directe jusqu'à 60 000 ans pour un coût de 100 000 années-modèle contre 500 000 pour la référence. La distribution des maxima par blocs, pondérée par les probabilités de branchement, a bien coïncidé avec la référence, en particulier dans la queue. Le seuil NINO3 a augmenté rapidement dans les premières générations, puis a atteint un plateau autour de 4,2 °C, ce qui indique une possible limite interne de l'extrémité dans le modèle. Le taux de succès de génération de descendants dépassant le seuil était d'environ 77 %. L'analyse des trajectoires a montré qu'une petite perturbation peut d'abord provoquer une Nina, puis renforcer l'El Niño suivant par une croissance non normale des perturbations.
Ce travail montre que l'échantillonnage des événements rares peut fournir des estimations statistiquement fiables de la queue de la distribution d'El Niño à un coût nettement inférieur à celui de la simulation directe. Contrairement aux travaux voisins, où les méthodes ont été appliquées à des systèmes plus simples ou différents, ici un accord quantitatif avec la référence est atteint pour la première fois pour un modèle ENSO de complexité intermédiaire. Cela renforce la confiance dans la transférabilité de l'approche aux modèles climatiques complets, où les simulations directes de la durée nécessaire sont impossibles.
Une optimisation des hyperparamètres de TEAMS (taux d'élimination, nombre d'ancêtres, temps d'avance) est possible par recherche bayésienne ou analyse par grille. Il est prometteur d'étudier des climats non stationnaires, par exemple sous réchauffement anthropique, où la statistique des El Niño extrêmes peut changer. Il convient également d'explorer des stratégies alternatives de gestion du temps de splitting pour éviter un biais vers des pics plus précoces et améliorer l'efficacité du boost.
La méthode est applicable aux modèles climatiques globaux coûteux pour évaluer les risques d'El Niño extrêmes.
Appliquer TEAMS à un modèle climatique complet, par exemple l'un des modèles CMIP, et comparer avec les longues simulations disponibles.
En climatologie, la forme de la queue de la distribution des El Niño extrêmes et sa sensibilité au changement climatique restent des questions ouvertes. Ce travail fournit un outil pour obtenir des estimations fiables dans des conditions de ressources de calcul limitées, ce qui est directement lié au problème de la quantification des risques climatiques rares.
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
Les auteurs ont appliqué la méthode d'échantillonnage des événements rares TEAMS au modèle de Zebiak-Cane pour estimer les temps de retour des El Niño extrêmes. La comparaison avec une simulation directe de référence de 500 000 ans a montré que TEAMS reproduit les estimations pour un coût d'environ 20 % de celui de la simulation directe.
Où est sa force
Ce travail vérifie pour la première fois de manière systématique l'algorithme TEAMS sur un modèle intermédiaire de l'ENSO, obtenant un accord quantitatif avec la référence pour la queue de la distribution. Cela renforce la confiance dans la transférabilité de l'approche à des modèles climatiques complets.
Vers où avancer
À essayer
Pourquoi cela peut compter
Le développement de méthodes d'échantillonnage des événements rares pour les modèles climatiques pourrait améliorer l'évaluation des risques d'événements extrêmes tels que les sécheresses et les inondations associés à El Niño, et aider à la planification de l'adaptation au changement climatique.
Travaux voisins dans notre archive
L'archive est vide autour de ce travail : l'analyse la plus proche n'est qu'à 0.47, contre 0,64 d'ordinaire. Il semble qu'on ait peu foré à côté.
Le coût de TEAMS de 100 000 années-modèle équivaut à simuler le climat de 1997 à environ 102 000 après J.-C., mais il est réalisé grâce au branchement de seulement 50 courtes exécutions.
L'idée de brancher des trajectoires pour explorer des issues rares fait écho au multivers dans la nouvelle de Ted Chiang « L'histoire de ta vie », où différentes branches du futur sont envisagées simultanément.