Multidecadal variability of the ENSO early-winter teleconnection to the North Atlantic in CMIP6 models: robustness and mechanisms

Technische Samenvatting: Multidecadale Variabiliteit van de ENSO Vroegwinter-teleconnectie naar de Noord-Atlantische Oceaan in CMIP6-modellen

Probleemstelling
El Niño-Southern Oscillation (ENSO) oefent een significante invloed uit op de Noord-Atlantische atmosferische circulatie tijdens het vroege winterseizoen (november–december), primair geprojecteerd op het East Atlantic Pattern (EAP). Hoewel operationele seizoensvoorspellingssystemen over het algemeen het ruimtelijke patroon van deze teleconnectie vastleggen, hebben recente observationele analyses aangegeven dat de sterkte van de ENSO-EAP-koppeling niet-stationair is en substantiële multidecadale variabiliteit vertoont. Zo is waargenomen dat de teleconnectie in verschillende tijdperken verschuift tussen een correlatie met de Noord-Atlantische Oscillatie (NAO) en de EAP. Het beoordelen van deze niet-stationariteit met behulp van observationele records wordt echter beperkt door de beperkte lengte van de instrumentele periode. De kleine steekproefomvang die inherent is aan sub-perioden van de observationele records maakt het moeilijk om statistisch te onderscheiden of de waargenomen veranderingen in de teleconnectie worden gedreven door genuwe dynamische mechanismen of dat ze louter artefactën zijn van steekproefvariabiliteit. Bov�elijk is het koppelen van deze veranderingen aan laagfrequente klimaatmodi, zoals de Pacific Decadal Oscillation (PDO) of de Atlantic Multidecadal Variability (AMV), uitdagend wanneer het observationele record slechts één of twee cycli van deze modi beslaat.

Methodologie
Om de beperkingen van de steekproefomvang te overwinnen en de stationariteit van de teleconnectie robuust te beoordelen, maakt deze studie gebruik van lange pre-industriële controle (piControl) simulaties van 10 Coupled Model Intercomparison Project Phase 6 (CMIP6) modellen. Deze simulaties, die meer dan 700 jaar beslaan met vaste antropogene forcings (1850-niveaus) en voorgeschreven natuurlijke forcings, bieden een grote steekproefomvang die het mogelijk maakt om de data te verdelen in afzonderlijke multidecadale sub-perioden zonder statistische degradatie.

De analyse hanteert het volgende methodologische kader:

1. Indexdefinitie: De Niño 3.4-index (N34) karakteriseert de ENSO-variabiliteit. De East Atlantic Index (EA) wordt afgeleid via Empirical Orthogonal Function (EOF) analyse van Mean Sea Level Pressure (MSLP) anomalieën in de Noord-Atlantische Oceaan, uitgelijnd met het ERA5-reanalysepatroon.
2. Beoordeling van Niet-Stationariteit: De studie berekent 21-jarige voortschrijdende correlaties tussen EA en N34. Om te bepalen of de waargenomen variabiliteit in deze correlaties de verwachting overstijgt die voortvloeit uit louter steekproefonzekerheid, wordt een lineair stochastisch model benadering (Rieke et al., 2021) toegepast. Dit omvat het genereren van 5.000 realisaties van een stationaire teleconnectie om een nulverdeling van de correlatievariabiliteit te creëren.
3. Periodeclassificatie: Perioden met statistisch significante positieve EA-N34 correlaties worden gelabeld als "POSCORR", terwijl perioden met zwakke correlaties worden gelabeld als "WEAKCORR".
4. Mechanisme-diagnose: De studie vergelijkt composiet-climatologieën en ENSO-regressiepatronen tussen POSCORR- en WEAKCORR-perioden. De onderzochte kernvariabelen zijn MSLP, 200 hPa geopotentiële hoogte (Z200), 200 hPa winden (u200), Outgoing Longwave Radiation (OLR) als proxy voor convectie, en Sea Surface Temperature (SST). Wave Activity Flux (WAF) wordt gebruikt om stationaire Rossby-golfpropagatie te diagnosticeren.

Belangrijkste Resultaten

1. Significante Multidecadale Variabiliteit: De analyse onthult dat de ENSO-EAP teleconnectie substantiële multidecadale variabiliteit vertoont in de CMIP6-modellen. Cruciaal is dat de vergelijking met het stochastische model aangeeft dat in 5 van de 10 modellen de variabiliteit van de 21-jarige voortschrijdende correlaties significant groter is dan wat op basis van steekproefonzekerheid alleen verwacht zou mogen worden. Dit ondersteunt de hypothese dat de teleconnectie niet-stationair is en wordt gedreven door dynamische mechanismen in plaats van willekeurige ruis.
2. Modelbiases: Hoewel de modellen in grote lijnen het ruimtelijke patroon van de ENSO-EAP teleconnectie vastleggen (negatieve MSLP-anomalieën ten westen van de Britse Eilanden), onderschatten zij de magnitude van de teleconnectie-sterkte aanzienlijk vergeleken met de ERA5-reanalyse. Deze onderschatting is consistent met bekende biases in GCM's en seizoensvoorspellingssystemen.
3. Dynamische Mechanismen: De studie identificeert drie primaire mechanismen die potentieel de multidecadale variabiliteit drijven:
   • Tropische Indische Oceaan Convectie: In POSCORR-perioden is er een sterkere ENSO-gerelateerde dipool in de tropische convectie boven de Indische Oceaan (verhoogde OLR-anomalieën). Dit versterkt de upper-level divergentie en fungeert als een krachtigerere bron van Rossby-golven, wat met name de wave activity flux vanuit de Indische Oceaan/Oost-Azië regio naar de Noord-Pacific versterkt.
   • Toestand van de Noord-Pacifische Jetstream: POSCORR-perioden zijn geassocieerd met een klimatologische verschuiving in de Noord-Pacifische jetstream (u200), gekenmerkt door een polaire verplaatsing en oostelijke elongatie. Deze configuratie faciliteert een efficiëntere propagatie van golfenergie van de tropen door de Noord-Pacific naar Noord-Amerika en vervolgens naar de Noord-Atlantische Oceaan.
   • Pacific Decadal Oscillation (PDO) Fase: De verschillen in de SST-climatologie tussen POSCORR- en WEAKCORR-perioden lijken op respectievelijk de negatieve en positieve fasen van de PDO. De studie vindt dat de ENSO-EAP teleconnectie sterker is tijdens perioden die lijken op de negatieve PDO-fase (geassocieerd met een sterkere extratropische meridionale SST-gradiënt in de Noord-Pacifische Oceaan) en zwakker tijdens positieve PDO-achtige fasen.
4. Golfpropagatie: Het verschil in Z200 regressiepatronen vertoont een golfachtige structuur over Oost-Azië en de Noord-Pacifische Oceaan die overeenkomt met de Indische Oceaan convectiedipool. Dit suggereert dat veranderingen in het traject van de Indische Oceaan de jetstream en de storm track respons in de Noord-Pacifische Oceaan moduleren, wat vervolgens de eddy-mean flow feedback beïnvloedt die de teleconnectie in de Noord-Atlantische Oceaan vormgeeft.

Betekenis en Claims
De auteurs claimen dat hun primaire bijdrage het bieden van een robuust kader is om de niet-stationariteit van de ENSO-EAP teleconnectie te beoordelen, waarbij de beperkingen van de steekproefomvang die inherent zijn aan observationele studies worden overwonnen. Door aan te tonen dat de multidecadale variabiliteit waarschijnlijk wordt gedreven door dynamische mechanismen in plaats van steekproefruis, valideert de studie de observationele bevindingen van verschuivende teleconnectiepatronen.

De auteurs stellen dat de geïdentificeerde mechanismen — specifelijk de modulatie van het Indische Oceaan convectietraject, de achtergrondtoestand van de Noord-Pacifische jetstream en de fase van de PDO — plausibele drijvers zijn van deze variabiliteit. Deze bevindingen verbeteren het begrip van de seizoensgebonden variabiliteit en voorspelbaarheid in de Noord-Atlantische Oceaan en Europa. De studie concludeert dat hoewel modelbiases bestaan, de consistentie van de signalen over verschillende modellen en hun uitlijning met eerdere observationele en modelleringsstudies (bijv. met betrekking tot de rol van de Indische Oceaan en de PDO) de hypothese van een betekenisvolle niet-stationariteit ondersteunt. De auteurs merken op dat hoewel hun resultaten robuust zijn tegen steekproefvariabiliteit, de aggregatie van modellen met verschillende biases het signaal in multi-model gemiddelden kan verzwakken, wat wijst op een behoefte aan gerichte sensitiviteitsexperimenten om deze mechanismen verder te isoleren.