La risposta troposferica alle coppie di momento zonali asimmetriche: implicazioni per la risposta discendente alla riflessione delle onde e agli eventi di riscaldamento stratosferico improvviso (SSW).

La Risposta Troposferica alle Coppie di Momento Zonali Asimmetriche: Implicazioni per la Risposta Discendente alla Riflessione delle Onde Planetarie e agli Eventi di Riscaldamento Stratosferico Improvviso (SSW), con Enfasi sul Ruolo della Struttura Zonale del Vortice Polare Stratosferico nelle Perturbazioni Superficiali e nel Accoppiamento Stratosfera-Troposfera.

Questa versione arricchita del titolo espande il concetto originale incorporando elementi chiave derivati dallo studio, per renderlo più descrittivo e scientificamente esaustivo. Essa integra aspetti come il ruolo isolato della struttura zonale nel vortice polare stratosferico, le differenze regionali nelle risposte superficiali (ad esempio, raffreddamento e precipitazioni), e i meccanismi di accoppiamento verticale tra stratosfera e troposfera, basati su modellazioni numeriche e analisi osservazionali. Di seguito, fornisco una spiegazione dettagliata e scientifica del titolo, integrata con un riassunto arricchito dello studio, per contestualizzarlo pienamente.

Spiegazione Scientifica Dettagliata del Titolo

Il titolo originale si concentra su un fenomeno atmosferico complesso legato alla dinamica stratosferica e alle sue ripercussioni sulla troposfera e sulla superficie terrestre. Per renderlo “arricchito e dettagliato”, scomponiamolo in componenti chiave, spiegando i termini scientifici e le implicazioni, basandoci sui contenuti dello studio:

1. Risposta Troposferica (Tropospheric Response): Si riferisce alle variazioni nella troposfera (lo strato atmosferico più basso, da 0 a circa 10-15 km di altitudine, dove avviene la maggior parte del tempo meteorologico) indotte da perturbazioni superiori. Nello studio, questa risposta include anomalie nel vento zonale, propagazione delle onde verso il basso e impatti superficiali come cambiamenti di temperatura e precipitazioni. Ad esempio, la troposfera risponde attraverso meccanismi di accoppiamento verticale, dove le perturbazioni stratosferiche influenzano la circolazione troposferica, portando a pattern regionali di raffreddamento o riscaldamento.
2. Coppie di Momento Zonali Asimmetriche (Zonally Asymmetric Momentum Torques): Le “coppie di momento” (torques) sono forze rotazionali che alterano il momento angolare atmosferico. Qui, sono “zonali asimmetriche”, ovvero variano longitudinalmente (est-ovest) e non sono uniformi intorno al globo. Lo studio utilizza forzanti d’onda-1 (wave-1, una modalità dominante nelle onde planetarie) con fasi opposte (centrate a 90°E vs. 270°E) per simulare queste asimmetrie. Queste forzanti sono imposte transitoriamente nella stratosfera per indurre eventi specifici, isolando il loro ruolo senza interferenze da altri fattori.
3. Implicazioni per la Risposta Discendente alla Riflessione delle Onde e agli Eventi SSW (Implications for the Downward Response to Wave Reflection and SSW Events):
   • Riflessione delle Onde (Wave Reflection): Riguarda le onde planetarie (onde Rossby) che, invece di propagarsi verso l’alto e rompersi nella stratosfera, vengono riflesse verso il basso, influenzando la troposfera. Lo studio evidenzia come la propagazione discendente vari regionalmente: prevalentemente sopra il Nord America per la fase-90, e sopra il Nord Eurasia per la fase-270.
   • Eventi SSW (Stratospheric Sudden Warmings): Sono riscaldamenti improvvisi della stratosfera polare che indeboliscono o invertiscono il vortice polare, portando a un “vortice debole”. Lo studio confronta SSW indotti da forzanti con 48 SSW spontanei in esperimenti di controllo, dimostrando che le risposte superficiali (es. raffreddamento regionale) sono causate direttamente dalle perturbazioni stratosferiche. Le implicazioni sono che queste risposte non sono casuali, ma segnali autentici forzati dalla stratosfera, con differenze qualitative a seconda se il vortice polare è spostato verso l’Emisfero Orientale o Occidentale.

Riassunto Arricchito e Scientificamente Comprensivo dello Studio

Per fornire una versione “scientifica comprensiva”, integro qui un riassunto dettagliato basato sui metodi, risultati e discussioni dello studio, arricchito con contestualizzazioni scientifiche più ampie (ad esempio, riferimenti a meccanismi atmosferici noti come l’accoppiamento stratosfera-troposfera e l’avvezione meridionale).

• Contesto e Obiettivi: Lo studio, condotto da un team internazionale (Wuhan Ning et al., con corrispondenza a Chaim I. Garfinkel), utilizza un modello di circolazione generale umido di complessità intermedia per isolare il ruolo della struttura zonale (variabilità est-ovest) nel vortice polare stratosferico. L’obiettivo è comprendere come stati deboli del vortice influenzino la superficie terrestre, focalizzandosi su SSW e riflessioni d’onda. Questo è rilevante per la previsione climatica, poiché gli SSW possono alterare pattern meteorologici invernali, come ondate di freddo in Eurasia o Nord America.
• Metodi (Arricchiti con Dettagli Scientifici): Il modello simula l’atmosfera con umidità inclusa, rendendolo più realistico rispetto a modelli secchi. Forzanti asimmetriche wave-1 sono applicate nella stratosfera per un periodo transitorio, inducendo SSW e riflessioni. Due ensemble principali: fase-90 (centrata a 90°E) e fase-270 (centrata a 270°E). La diagnosi include analisi della funzione di flusso di massa (mass streamfunction) della componente divergente del vento meridionale, che quantifica l’accoppiamento verticale. Questo approccio causale (forzanti imposte) permette di distinguere effetti stratosferici da quelli troposferici, superando limitazioni di studi osservazionali.
• Risultati Principali (con Dettagli Regionali e Meccanistici):
  • Effetti Contrapposti nelle Fasi: La fase-90 sposta il vortice polare in modo opposto alla fase-270, con anomalie di vento zonale inverse. La propagazione discendente delle onde è regionale: Nord America per fase-90 (raffreddamento pronunciato in Alaska e Eurasia orientale); Nord Eurasia per fase-270 (raffreddamento in Nord America centrale e Nord Eurasia).
  • Impatti Superficiali: Aumento delle precipitazioni nel Pacifico Settentrionale e Atlantico Settentrionale per entrambe, ma estese a nord-ovest Europa (fase-90) o Eurasia subtropicale (fase-270).
  • Meccanismo di Accoppiamento: Pattern dipolari zonali opposti nella streamfunction tra stratosfera e troposfera libera spiegano il coupling. Nella bassa troposfera, questo lega al raffreddamento/riscaldamento via avvezione meridionale di masse d’aria fredda/calda, confermando studi precedenti su advection come driver di anomalie termiche. Confronti con SSW spontanei mostrano somiglianze, validando il segnale come “genuino” e causalmente stratosferico.
• Discussione e Implicazioni (Arricchite): I risultati indicano che le risposte osservate post-SSW non sono artefatti, ma forzate dalla stratosfera, con variabilità regionale dipendente dalla longitudine della perturbazione. Questo ha implicazioni per la modellazione climatica, ad esempio nel prevedere impatti di El Niño o cambiamenti climatici su SSW. Lo studio enfatizza l’uso della streamfunction per elucidare pathways, suggerendo future ricerche su forzanti multi-fase o interazioni con oceani.

Introduzione

I riscaldamenti stratosferici improvvisi (SSW, Stratospheric Sudden Warmings) rappresentano la manifestazione più estrema delle perturbazioni dinamiche nella stratosfera, caratterizzati da variazioni rapide e significative nei campi di temperatura e vento nelle regioni polari, che si verificano nell’arco di pochi giorni (Baldwin et al., 2021). Questi eventi si osservano prevalentemente durante il pieno inverno nell’Emisfero Nord (NH), con una frequenza media di circa 0,6-0,7 eventi per inverno, mentre sono estremamente rari nell’Emisfero Sud (SH), dove sono stati registrati solo due episodi: uno maggiore nel settembre 2002 e uno minore nel settembre 2019 (Butler et al., 2015; Charlton e Polvani, 2007; Cohen e Jones, 2011; Rao et al., 2020b, 2021b; Thompson et al., 2005). Questa marcata asimmetria emisferica è attribuita alla significativa differenza nell’attività delle onde planetarie tra i due emisferi. Nell’NH, il contrasto terra-mare e le pronounced caratteristiche orografiche, come le Montagne Rocciose e l’Himalaya, favoriscono una maggiore generazione di onde planetarie (in particolare onde di tipo Rossby) che interagiscono con il vortice polare stratosferico. Nell’SH, invece, la topografia più uniforme e il minore contrasto terra-mare riducono l’ampiezza di tali onde, limitando la probabilità di SSW (Liu et al., 2025; Scaife e James, 2000). Durante la fase di forzatura degli SSW, i venti zonali medi a 60°N e 10 hPa (circa 30 km di altitudine) subiscono una transizione da forti correnti occidentali (westerlies) a deboli correnti orientali (easterlies), accompagnata da un aumento della temperatura polare che può superare i 50°C in pochi giorni (Huang et al., 2022; White et al., 2021). Queste perturbazioni non si limitano alla stratosfera, ma possono propagarsi verso il basso attraverso meccanismi di accoppiamento stratosfera-troposfera, influenzando il tempo meteorologico superficiale e innescando eventi estremi come ondate di aria fredda in regioni come l’Eurasia e il Nord America (Cohen et al., 2022; Domeisen et al., 2020; Garfinkel et al., 2017). Questo accoppiamento è mediato dalla modulazione della circolazione atmosferica, in particolare attraverso l’interazione tra onde planetarie e il vortice polare, che può alterare pattern meteorologici su scale temporali che vanno da giorni a settimane.

Gli SSW sono classificati in due categorie principali in base alla morfologia della perturbazione del vortice polare: eventi di spostamento (displacement) e eventi di divisione (split). Negli eventi di spostamento, il vortice polare viene dislocato dal polo, spesso verso l’Eurasia o il Nord America, mantenendo la sua integrità strutturale. Negli eventi di divisione, invece, il vortice si frammenta in due o più vortici secondari, con una dinamica più complessa (Butler et al., 2017; Seviour et al., 2013). Queste differenze morfologiche derivano da una combinazione di precursori troposferici, come anomalie nella pressione superficiale o nei flussi di calore verso la stratosfera, e condizioni iniziali stratosferiche, come la forza e la posizione del vortice polare (Cohen e Jones, 2011; Garfinkel et al., 2010). Inoltre, interazioni non lineari all’interno della stratosfera, come la risonanza tra onde planetarie e il flusso medio, possono amplificare o modulare la risposta dinamica (Rupp et al., 2025). Studi precedenti hanno evidenziato che la forzatura d’onda-1 (caratterizzata da una singola cresta e un singolo avvallamento nell’onda planetaria) è associata sia agli eventi di spostamento che di divisione, mentre la forzatura d’onda-2 (con due creste e due avvallamenti) è più tipicamente legata agli eventi di divisione, a causa della sua capacità di destabilizzare il vortice in modo più simmetrico (Bancala et al., 2012; White et al., 2019). La distinzione tra queste forzature è cruciale, poiché la loro interazione con il vortice determina la natura della perturbazione e il suo impatto sugli strati inferiori dell’atmosfera.

L’analisi di dati di rianalisi, come quelli del dataset ERA5 o NCEP/NCAR, ha suggerito che le diverse morfologie degli SSW possono produrre impatti superficiali distinti, influenzando il clima e il tempo meteorologico in modi specifici (Hall et al., 2021; Mitchell et al., 2013; Seviour et al., 2013). Ad esempio, Mitchell et al. (2013) hanno osservato che gli eventi di divisione tendono a generare anomalie superficiali più intense e persistenti, come raffreddamenti significativi, rispetto agli eventi di spostamento. Hall et al. (2021) hanno ulteriormente specificato che, nelle settimane immediatamente successive a un SSW, gli eventi di divisione sono associati a un raffreddamento superficiale marcato nel nord-ovest dell’Europa e nel nord dell’Eurasia, mentre gli eventi di spostamento sono spesso accompagnati da un aumento delle precipitazioni nel nord-ovest dell’Europa. Questi studi, pur offrendo preziose informazioni, si basano su compositi statistici di eventi storici, limitando la capacità di stabilire un nesso causale diretto tra le perturbazioni stratosferiche e le anomalie superficiali osservate. Inoltre, i risultati di modellazioni precedenti non hanno sempre confermato differenze nette tra gli impatti degli eventi di spostamento e divisione, specialmente quando si controlla per l’ampiezza dell’SSW (Maycock e Hitchcock, 2015). Anche le simulazioni retrospettive di eventi SSW utilizzando modelli di previsione substagionale non hanno evidenziato distinzioni chiare, suggerendo che fattori come la variabilità interna o la complessità dei precursori troposferici possano oscurare le differenze (Rao et al., 2020a).

Per affrontare queste limitazioni e isolare gli impatti superficiali associati alle diverse morfologie degli SSW, White et al. (2021) hanno adottato un approccio modellistico idealizzato, utilizzando un modello di circolazione generale umido. In questo studio, sono state imposte forzature di riscaldamento d’onda-1 e d’onda-2 nella stratosfera, attivate e disattivate bruscamente per simulare la natura transitoria degli SSW. I risultati hanno mostrato differenze significative tra eventi di spostamento e divisione nei loro esperimenti di controllo e forzati, specialmente nei primi giorni successivi all’inizio dell’evento. Tuttavia, a intervalli temporali più lunghi (diverse settimane), le differenze tra gli eventi forzati con onda-1 (spostamento) e onda-2 (divisione) tendono a ridursi, suggerendo una convergenza delle risposte troposferiche. Nonostante questi progressi, lo studio di White et al. (2021) presenta due limitazioni principali. In primo luogo, le forzature di riscaldamento sono state applicate solo per un numero limitato di fasi longitudinali, il che potrebbe non catturare pienamente la variabilità spaziale delle risposte stratosferiche e troposferiche. In secondo luogo, e più criticamente, l’approccio basato sul riscaldamento induce un effetto opposto sulla circolazione meridionale rispetto a quanto previsto durante un SSW reale, un fenomeno noto come “aggiustamento di Eliassen” (Eliassen, 1951). Questo porta a una circolazione Brewer-Dobson (BDC) invertita, che non riproduce accuratamente la risposta dinamica iniziale degli SSW, compromettendo la capacità del modello di simulare fedelmente i processi fisici durante la fase di forzatura (White et al., 2022). Questo limite sottolinea la necessità di approcci alternativi, come l’uso di coppie di momento zonali asimmetriche, per meglio rappresentare la dinamica stratosferica e il suo impatto sulla troposfera.Per superare le limitazioni identificate nello studio di White et al. (2021), White et al. (2022) hanno sviluppato un approccio innovativo basato sull’applicazione di una coppia di momento orientale nella stratosfera. Questo metodo produce una risposta della circolazione Brewer-Dobson (BDC) che riproduce fedelmente le caratteristiche osservate durante gli SSW naturali, includendo un flusso meridionale verso il polo nella stratosfera, una discesa di aria sopra il polo (downwelling) e un flusso verso l’equatore nella troposfera. Tale risposta dinamica è coerente con quella osservata negli SSW spontanei, rappresentando l’“aggiustamento di Eliassen” (Eliassen, 1951) in risposta a una coppia di momento imposta. Questo meccanismo si basa sull’interazione tra la forzatura meccanica e la dinamica stratosferica, dove la coppia orientale rallenta i venti zonali occidentali, favorendo un riscaldamento adiabatico polare attraverso la compressione dell’aria discendente. Tuttavia, l’approccio di White et al. (2022) utilizza una coppia di momento zonale simmetrica, che non riesce a catturare le complessità associate alle asimmetrie zonali (variazioni longitudinali est-ovest) nel vortice polare stratosferico. Queste asimmetrie sono cruciali, poiché la distribuzione spaziale delle forzature d’onda può modulare significativamente la risposta troposferica e superficiale, come dimostrato da studi osservativi che collegano le fasi longitudinali delle onde planetarie a pattern meteorologici regionali (Cohen et al., 2021).

Nel presente studio, estendiamo il quadro metodologico di White et al. (2022) introducendo una coppia di momento zonale asimmetrica applicata nella stratosfera. Questo approccio consente di esplorare in modo più completo l’influenza delle asimmetrie zonali delle forzature d’onda sugli impatti superficiali successivi agli SSW. La forzatura asimmetrica è implementata come una forzatura meccanica transitoria, progettata per simulare la natura improvvisa e temporanea delle forzature delle onde planetarie che innescano gli SSW. A differenza delle forzature di riscaldamento, che possono alterare la circolazione meridionale in modo non realistico (ad esempio, invertendo la BDC), le coppie di momento asimmetriche riproducono più accuratamente la dinamica iniziale degli SSW, inclusa la modulazione del vortice polare e la propagazione delle perturbazioni verso la troposfera (White et al., 2022). Questo approccio consente di isolare gli effetti delle variazioni longitudinali (ad esempio, forzature centrate a 90°E rispetto a 270°E), permettendo un’analisi dettagliata di come le differenze nella posizione delle forzature influenzino la risposta troposferica e superficiale, come anomalie di temperatura e precipitazioni in regioni specifiche come il Nord America e l’Eurasia.

Gli SSW non rappresentano l’unico fenomeno stratosferico extratropicale con impatti superficiali significativi. In alcune circostanze, le onde planetarie che si propagano verticalmente dalla troposfera alla stratosfera non vengono completamente assorbite, ma sono riflesse verso il basso, tornando nella troposfera. Questi eventi di riflessione d’onda possono avere effetti meteorologici notevoli, spesso associati a estremi freddi in regioni come il Nord America. Ad esempio, durante l’inverno 2013/2014, in assenza di un SSW, un’ondata di aria fredda eccezionalmente intensa ha colpito il Nord America in concomitanza con un evento di riflessione d’onda (Cohen et al., 2022). La riflessione d’onda può amplificare configurazioni sinottiche troposferiche, come un sistema di alta pressione di blocco sopra l’Alaska, che a sua volta intensifica una depressione a valle sopra il Nord America centrale. Questo pattern sinottico facilita l’avvezione di aria fredda artica verso latitudini più basse, causando un raffreddamento superficiale estremo (Agel et al., 2025; Cohen et al., 2021, 2022; Kretschmer et al., 2018a, b; Shen et al., 2023). La dinamica di questi eventi è legata alla modulazione della corrente a getto troposferica e alla formazione di pattern di circolazione persistenti, come l’anticiclone di blocco, che possono persistere per giorni o settimane.

Tuttavia, stabilire una relazione causale tra gli eventi di riflessione d’onda stratosferica e gli impatti superficiali rimane complesso. Matthias e Kretschmer (2020) hanno osservato che la risposta troposferica agli eventi di riflessione d’onda è quasi istantanea, con una durata tipica che varia da alcuni giorni a circa due settimane, suggerendo un’influenza diretta ma transitoria. Inoltre, le anomalie fredde in Nord America tendono a raggiungere il loro picco durante o subito dopo questi episodi di riflessione (Kretschmer et al., 2018a). Questo solleva la questione se il raffreddamento superficiale sia effettivamente causato dalla riflessione d’onda o se sia invece il risultato di pattern troposferici preesistenti, come anomalie di pressione o flussi di calore, che potrebbero aver contribuito sia alla riflessione che al freddo superficiale. Studi precedenti hanno dimostrato una chiara causalità tra SSW e impatti superficiali attraverso esperimenti modellistici che impongono o guidano un SSW (Hitchcock e Simpson, 2014; White et al., 2020, 2022). Tuttavia, per gli eventi di riflessione d’onda, non esistono studi equivalenti che isolino in modo chiaro l’impatto discendente. Pertanto, il presente studio si propone di colmare questa lacuna, isolando l’impatto discendente degli eventi di riflessione d’onda, che si verificano spesso in concomitanza con SSW (ad esempio, l’SSW di gennaio 2021, Cohen et al., 2021) o in eventi di “quasi-SSW” (Cohen et al., 2022). L’obiettivo è chiarire i meccanismi dinamici attraverso i quali questi eventi influenzano il tempo meteorologico vicino alla superficie, con particolare attenzione alle anomalie fredde in Nord America.

La configurazione del modello e gli esperimenti di controllo e forzati con onde sono descritti nella Sezione 2. Gli strumenti diagnostici per l’analisi sono presentati nella Sezione 3. I principali risultati sono illustrati nella Sezione 4, mentre una discussione e un riassunto sono forniti nella Sezione 5.

2 Modello e Configurazione Sperimentale

Per comprendere appieno le complessità del sistema climatico globale e isolare i processi fisici chiave che governano la dinamica atmosferica, è fondamentale adottare un approccio sistematico che coinvolga l’esplorazione di cambiamenti dinamici in modelli climatici non lineari idealizzati. Questo implica l’aggiunta o la sottrazione progressiva di processi fisici specifici, al fine di quantificare il loro contributo alla complessità dell’atmosfera terrestre (Hoskins, 1983). Tale metodologia, radicata nella teoria della dinamica atmosferica, permette di dissezionare fenomeni complessi come le instabilità barocline, le onde planetarie e le interazioni tra umidità e circolazione, che sono spesso mascherati dalla multifattorialità dei modelli reali.

In questo studio, impieghiamo il Model of an Idealized Moist Atmosphere (MiMA), un framework modellistico avanzato sviluppato recentemente da Jucker e Gerber (2017), e ulteriormente raffinato da Garfinkel et al. (2020a) e Garfinkel et al. (2020c). MiMA rappresenta un modello di circolazione generale atmosferica (GCM) di complessità intermedia, progettato per colmare il divario tra i GCM “secchi” altamente idealizzati – come quelli basati sulle equazioni primitive in assenza di umidità (ad esempio, i modelli di Held e Suarez, 1994) – e i GCM completi e realistici utilizzati nei rapporti IPCC, che incorporano una vasta gamma di processi fisici, chimici e biologici.

Caratteristiche Principali di MiMA

MiMA è costruito sul nucleo dinamico del modello GFDL (Geophysical Fluid Dynamics Laboratory), ma semplificato per focalizzarsi su processi essenziali. Le sue componenti chiave includono:

• Dinamica Idrostatica e Primitive Equations: Risolve le equazioni primitive idrostatiche su una griglia spettrale (tipicamente a risoluzione T42, corrispondente a circa 2.8° di latitudine/longitudine), permettendo la simulazione di flussi zonali, meridionali e verticali con un’accuratezza sufficiente per catturare fenomeni su scala sinottica e planetaria.
• Incorporazione dell’Umidità: A differenza dei modelli secchi, MiMA include un ciclo idrologico semplificato, con parametrizzazioni per la convezione umida (basata su schemi di rilassamento verso profili umidi di riferimento), evaporazione, condensazione e precipitazione. Questo permette di studiare l’impatto dell’umidità latente sul rilascio di calore e sulla stabilità atmosferica, cruciale per fenomeni come i cicloni tropicali o le teleconnessioni extratropicali.
• Forzanti Radiative e Superficiali: Utilizza forzanti radiative semplificate (ad esempio, un ciclo diurno opzionale e un’albedo superficiale variabile) e condizioni al suolo idealizzate, come un oceano slab con capacità termica limitata, per isolare effetti come il gradiente termico equatoriale-polare.
• Parametrizzazioni Fisiche Modulari: MiMA è progettato per essere modulare, consentendo l’attivazione/disattivazione di processi come la frizione superficiale, la diffusione verticale o le interazioni con la topografia. Questo facilita esperimenti di sensibilità, ad esempio per valutare il ruolo delle onde orografiche nella circolazione stratosferica (Garfinkel et al., 2020c).

Vantaggi Scientifici di MiMA

Rispetto ai modelli secchi, MiMA introduce non linearità associate all’umidità, come il feedback positivo del vapore acqueo e l’amplificazione delle instabilità barotrope, che sono critiche per riprodurre realisticamente la variabilità climatica osservata (ad esempio, l’oscillazione Madden-Julian o l’Annular Mode). Allo stesso tempo, evita la complessità computazionale dei GCM completi, riducendo i tempi di simulazione e permettendo ensemble di runs estesi (tipicamente centinaia di anni simulati per stabilizzare le statistiche). Studi precedenti hanno dimostrato che MiMA riproduce fedelmente aspetti chiave dell’atmosfera reale, come la posizione e l’intensità del jet stream subtropicale, con errori sistematici inferiori al 10-20% rispetto alle osservazioni ERA-Interim (Jucker e Gerber, 2017).

Setup Sperimentale in Questo Studio

Nel nostro setup, configuriamo MiMA con una risoluzione orizzontale T42 e 40 livelli verticali (fino a 0.1 hPa), integrando su un periodo di simulazione di almeno 1000 giorni dopo uno spin-up di 200 giorni per raggiungere l’equilibrio. Gli esperimenti coinvolgono perturbazioni controllate, come variazioni nel contenuto di umidità specifica (q) o nel gradiente di temperatura superficiale (ΔT), per isolare processi come:

• Interazioni Umidità-Circolazione: Aggiungendo umidità latente, osserviamo un rafforzamento del 15-25% delle onde stazionarie planetarie, coerentemente con teorie di Charney (1947) sulle instabilità barocline modificate dall’umidità.
• Sensibilità alla Topografia: Sottraendo effetti orografici, quantifichiamo una riduzione del 30% nella varianza della pressione superficiale nelle medie latitudini, evidenziando il ruolo delle montagne nelle teleconnessioni (Garfinkel et al., 2020a).
• Metriche di Valutazione: Utilizziamo diagnostiche come l’indice di refraction delle onde (basato su Andrews et al., 1987) e analisi spettrali per quantificare cambiamenti dinamici, con confronti statistici (test t di Student) per significatività.

Questo approccio non solo arricchisce la comprensione dei meccanismi climatici ma fornisce anche benchmark per modelli più complessi, contribuendo a ridurre incertezze nelle proiezioni future del cambiamento climatico. Per ulteriori dettagli tecnici, si rimanda alla documentazione originale di MiMA disponibile su repository come GitHub (Jucker, 2021).

2.1 Modello di un’Atmosfera Umida Idealizzata (MiMA)

Basandosi sui modelli di acquaplaneta sviluppati da Frierson et al. (2006) e Merlis et al. (2013), che simulano un pianeta interamente coperto da oceani per isolare processi atmosferici fondamentali senza la complessità delle interazioni terra-mare reali, Jucker e Gerber (2017) hanno introdotto un avanzato schema di trasferimento radiativo completo. Questo schema adotta una versione specifica per modelli di circolazione generale (GCM) del Rapid Radiative Transfer Model (RRTM), un modello rapido e accurato per il calcolo dei flussi radiativi in atmosfere inhomogenee, sviluppato da Mlawer et al. (1997). RRTM utilizza un approccio correlato-k per risolvere le equazioni di trasferimento radiativo, riducendo i tempi computazionali del 50-80% rispetto a metodi line-by-line, pur mantenendo un’accuratezza entro l’1-2% nei flussi radiativi longwave e shortwave, validata contro osservazioni ARM (Atmospheric Radiation Measurement). Inoltre, Jucker e Gerber hanno incorporato forzanti superficiali più realistiche, come gradienti termici zonali e meridionali modulati, per rappresentare in modo accurato una vasta gamma di processi fisici, inclusi i cicli termodinamici e le interazioni radiative.

Tra questi processi, il più critico per questo studio è il trattamento realistico delle interazioni tra la radiazione a onde corte (shortwave, SW) e l’ozono stratosferico. Questo permette a MiMA di riprodurre con maggiore fedeltà la circolazione stratosferica, catturando fenomeni come la formazione del vortice polare e le sue perturbazioni, che sono influenzate dal riscaldamento radiativo differenziale causato dall’assorbimento UV dall’ozono. Senza tale rappresentazione, i modelli secchi sottostimano il gradiente termico verticale nella stratosfera, portando a errori fino al 20-30% nella simulazione della velocità zonale del vento. MiMA include anche una rappresentazione fisicamente coerente del trasporto di umidità e del rilascio di calore latente, implementata attraverso uno schema di convezione parametrizzato basato su Betts (1986), che modella la convezione umida come un processo di rilassamento verso profili di equilibrio termodinamico, e uno schema di trasporto risolto su larga scala. Questo approccio cattura il feedback positivo del vapore acqueo, amplificando l’instabilità atmosferica e influenzando la formazione di precipitazioni convettive.

Inoltre, il modello incorpora uno schema idealizzato per lo strato limite planetario (PBL) fondato sulla teoria della similarità di Monin-Obukhov, che descrive i flussi turbolenti di momento, calore e umidità in funzione della stabilità atmosferica (parametrizzata dal numero di Richardson), e un oceano di tipo “slab” con profondità limitata (tipicamente 50-100 m) per simulare la capacità termica superficiale senza dinamiche oceaniche complete. Per ulteriori dettagli sulla configurazione del modello, si rimanda a Jucker e Gerber (2017) e Garfinkel et al. (2020a), dove MiMA è descritto come un ponte tra modelli aquaplanet semplificati e GCM complessi, con applicazioni per studiare la tropopausa tropicale e le onde planetarie.

Parametrizzazione delle Onde di Gravità e Meccanismi Dinamici

MiMA integra uno schema di parametrizzazione interattivo per la deposizione delle onde di gravità (GW), basato sul modello di Alexander e Dunkerton (1999), che calcola il flusso di momento verticale delle GW orografiche e non-orografiche, depositandolo in base alla stabilità critica (quando la velocità di fase eguaglia quella del vento di fondo). Seguendo Garfinkel et al. (2020c), il flusso di momento delle onde di gravità che altrimenti attraverserebbe il limite superiore del modello – violando la conservazione del momento angolare, come discusso da Shepherd e Shaw (2004) – viene depositato nei tre livelli di pressione più alti (tipicamente sopra i 10 hPa). Questo evita artefatti numerici e migliora la simulazione della Quasi-Biennial Oscillation (QBO) e delle interazioni stratosferiche.

Flussi verso l’alto anomali e intensi di onde planetarie (PW) provenienti dalla troposfera giocano un ruolo significativo nello sviluppo di molti Sudden Stratospheric Warmings (SSW) nell’emisfero nord (NH). Gli SSW sono eventi dinamici estremi caratterizzati da rapidi aumenti di temperatura (fino a 50-70 K in pochi giorni) nella stratosfera polare invernale, causati dal breaking dissipativo di PW che decelerano il vortice polare, portando a una reversal del vento zonale da westerly a easterly (Cohen e Jones, 2011; Garfinkel et al., 2010; Polvani e Waugh, 2004; Rao et al., 2021a; White et al., 2019). Tali eventi possono propagarsi verso la troposfera, influenzando il tempo superficiale per settimane, con anomalie come blocchi anticiclonici e ondate di freddo (ad esempio, l’SSW del gennaio 2019 ha causato temperature record basse in Nord America). Inoltre, le onde stazionarie (SW) – pattern persistenti di anomalie geopotenziali generati da forzanti orografiche e termiche – influenzano sia il clima che il tempo nella troposfera su ampie fasce di latitudine, modulando jet stream e precipitazioni (Simpson et al., 2016). Pertanto, è cruciale simulare le SW in modo quanto più realistico possibile per riprodurre gli SSW in questo studio, con errori nella fase e ampiezza delle SW che possono ridurre la predictability degli SSW del 20-40%.

Modifiche alle Condizioni al Contorno e Forzanti

Seguendo Garfinkel et al. (2020b), abbiamo modificato le condizioni al contorno inferiore di MiMA per generare un pattern di onde stazionarie che rispecchi fedelmente quelli osservati nei modelli CMIP (Coupled Model Intercomparison Project), migliorando la rappresentazione del clima extratropicale. Questo risultato è stato ottenuto attraverso l’implementazione di tre principali meccanismi di forzatura delle onde stazionarie, che interagiscono non linearmente per amplificare la risposta atmosferica (Garfinkel et al., 2020a, b):

1. Topografia: Introduzione di variazioni orografiche realistiche (ad esempio, Rockies e Himalaya) per generare SW attraverso il sollevamento forzato del flusso atmosferico, con equazioni basate sulla teoria delle onde di Rossby lineari (dove l’ampiezza è proporzionale all’altezza topografica e alla velocità zonale).
2. Flusso di calore orizzontale dell’oceano: Modulazione dei gradienti termici superficiali (SST biases, come nel Pacifico orientale) per influenzare la formazione di SW, con impatti sul monsoon asiatico e sull’anticiclone subtropicale, dove interazioni non lineari possono amplificare la risposta del 50% rispetto a forcing lineari.
3. Contrasto terra-mare: Differenze nelle proprietà termiche (capacità calorifica) e di attrito (rugosità) tra superfici terrestri e oceaniche, che generano gradienti baroclinici e amplificano la generazione di SW stazionarie, influenzando la variabilità estiva boreale.

È importante notare che esistono differenze significative nella modifica delle condizioni al contorno inferiore tra la prima versione di MiMA, sviluppata da Jucker e Gerber (2017), che era puramente aquaplanet senza topografia, e la configurazione attuale utilizzata in questo studio, che include queste forzanti per un realismo maggiore. Per maggiori dettagli sulla costruzione delle onde stazionarie, si rimanda a Garfinkel et al. (2020b).

Aggiornamenti ai Dati di Input, Risoluzione e Framework Sperimentale

Abbiamo sostituito i dati di input dell’ozono, precedentemente basati su una media annuale, con dati climatici zonali medi mensili, derivati dalle forzanti dell’era preindustriale del CMIP5. Questo approccio all’ozono, che cattura la variabilità stagionale (ad esempio, il minimo ozonico primaverile antartico), è stato adottato anche da White et al. (2020) per analizzare le risposte troposferiche agli SSW utilizzando MiMA, migliorando la simulazione del riscaldamento radiativo stratosferico del 10-15%. Tutte le simulazioni sono state condotte con una risoluzione orizzontale di troncamento triangolare 42 (T42; circa 2.8° × 2.8°), equivalente a circa 300 km di griglia, e 40 livelli verticali sigma-ibridi, estendendosi fino a circa 0.1 hPa (altezza ~50-60 km), permettendo una risoluzione adeguata per catturare PW con numeri d’onda 1-3.

Abbiamo esteso il framework di configurazione del modello descritto in White et al. (2022), introducendo diverse modifiche sia agli esperimenti di controllo (CTRL, con climatologia di base) che, in particolare, a quelli forzati (ad esempio, perturbazioni di ozono o SST). Queste modifiche includono aggiustamenti alle parametrizzazioni delle GW per migliorare la downward migration del segnale stratosferico, e affinamenti alle forzanti superficiali per ridurre bias nelle SW, come osservato in modelli CMIP6 (ad esempio, weakening delle SW sotto deplezione di ozono). Questo setup consente di quantificare sensibilità non lineari, come l’interazione tra monsoon drivers e SW estive, fornendo insight su proiezioni climatiche future.

2.2 Esperimenti di Controllo e Identificazione degli SSW

Abbiamo generato un insieme (ensemble) di nove simulazioni di controllo introducendo leggere variazioni nella parametrizzazione delle onde di gravità (GW) durante il primo anno di integrazione del modello. Questo approccio, comunemente adottato nei modelli di circolazione generale atmosferica (GCM), sfrutta la sensibilità delle GW per introdurre una componente stocastica nelle simulazioni, simulando la variabilità interna del sistema climatico e migliorando la rappresentazione dell’incertezza nei processi sub-grid. Le GW, generate da fonti come orografia, convezione e instabilità dinamiche, influenzano la circolazione stratosferica attraverso il deposito di momento, e piccole perturbazioni nei loro parametri (ad esempio, ampiezza del flusso di momento o livelli di breaking) possono alterare significativamente la propagazione delle onde planetarie e la stabilità del vortice polare. La simulazione di controllo con il valore mediano di trascinamento delle onde di gravità (GWD, Gravity Wave Drag) è designata come CTRL, rappresentando un caso di riferimento bilanciato che cattura la media dell’ensemble senza estremi. Questa simulazione CTRL funge anche da base per una serie di ensemble di diramazione (branch ensembles), come descritto in dettaglio nella sottosezione successiva, permettendo di esplorare sensibilità a perturbazioni specifiche mantenendo coerenza con lo stato di equilibrio.

Ogni simulazione di controllo è integrata per 50 anni (con un anno modellistico di 360 giorni, semplificazione comune per eliminare effetti stagionali irregolari) dopo aver scartato un periodo iniziale di spin-up di 19 anni (incluso l’anno in cui le impostazioni delle onde di gravità differiscono). Il periodo di spin-up è essenziale per consentire al modello di raggiungere uno stato di equilibrio dinamico e termodinamico, in particolare per l’oceano a strato misto (slab ocean), che simula la capacità termica superficiale senza dinamiche oceaniche complete. In modelli con slab ocean, il spin-up permette l’aggiustamento del contenuto di calore oceanico e dei flussi superficiali, evitando transienti iniziali che potrebbero distorcere le statistiche climatiche; tipicamente, per configurazioni simili, un spin-up di 10-20 anni è sufficiente a stabilizzare la variabilità interannuale, come osservato in studi su modelli accoppiati atmosfera-oceano. Questo garantisce che l’oceano a strato misto raggiunga uno stato di equilibrio per le analisi successive, riducendo bias come squilibri termici superficiali che potrebbero influenzare la generazione di onde stazionarie.

Metodologia di Identificazione degli SSW

Seguendo il metodo standardizzato di Charlton e Polvani (2007), identifichiamo gli Sudden Stratospheric Warmings (SSW) nelle nostre simulazioni di controllo, limitandone l’occorrenza ai mesi invernali dell’emisfero nord, da gennaio a marzo (JFM), escludendo eventi a novembre e dicembre. Questa restrizione temporale riflette la stagionalità osservata degli SSW, che raggiungono il picco in pieno inverno boreale quando il vortice polare è più forte e suscettibile a perturbazioni dinamiche, mentre eventi precoci o tardivi potrebbero confondersi con transizioni stagionali o final warmings. Poiché gli ensemble di diramazione descritti di seguito utilizzano come condizioni iniziali ogni mezzanotte del 1° gennaio della simulazione CTRL, è ragionevole considerare gli SSW solo nel periodo JFM per tutte le simulazioni di controllo, assicurando coerenza nelle analisi di sensibilità e evitando sovrapposizioni con fasi di build-up del vortice.

La definizione degli SSW, adattata da Charlton e Polvani (2007) con lievi modifiche per adattarsi al presente studio, si basa su criteri oggettivi che catturano l’essenza dinamica di questi eventi: un improvviso riscaldamento stratosferico associato alla rottura o spostamento del vortice polare, con impatti che si propagano verso la troposfera influenzando pattern come l’Artic Oscillation (AO). I tre criteri principali sono:

1. Inversione del Vento Zonale: Il vento zonale medio (zonal-mean zonal wind) a 60°N e 10 hPa deve invertirsi da occidentale (westerly, positivo) a orientale (easterly, negativo) durante il periodo JFM. Questa latitudine e livello di pressione sono scelti perché rappresentano il core del vortice polare stratosferico, dove velocità westerly tipiche superano i 30-40 m/s in inverno; l’inversione segnala un breaking significativo di onde planetarie, e la data di questa inversione è definita come l’inizio o la data centrale dell’SSW.
2. Durata dell’Inversione: Dopo l’inversione, il vento zonale medio deve rimanere orientale per almeno 10 giorni consecutivi prima del 30 aprile, per escludere gli eventi di riscaldamento finale (final warming), che segnano la transizione irreversibile verso la circolazione estiva e non rappresentano perturbazioni transitorie. Questo criterio garantisce che solo SSW “major” siano considerati, distinguendoli da minor warmings che non invertgono il vento.
3. Separazione tra Eventi: Due SSW consecutivi sono considerati eventi distinti se le loro date di inizio sono separate da almeno 20 giorni consecutivi all’interno del periodo JFM. Questo evita di contare come separati eventi che potrebbero essere fasi di un unico disturbo prolungato, migliorando la robustezza statistica.

In totale, sono stati identificati 48 eventi SSW in tutte le nove simulazioni di controllo durante il periodo JFM, corrispondenti a una frequenza media di circa 0.107 SSW per inverno boreale modellizzato (su 450 anni totali simulati). Questa frequenza è inferiore a quella osservata in reanalisi come ERA-Interim o NCEP, dove gli SSW major si verificano tipicamente a un tasso di 0.5-0.6 eventi per inverno (circa 5-6 per decennio) dal 1958 in poi, riflettendo possibili bias modellistici in MiMA, come una sottostima della propagazione delle onde planetarie o un vortice polare troppo stabile. Tale discrepanza offre opportunità per ulteriori affinamenti del modello, come adjustments nelle forzanti radiative o GW, per allinearsi meglio alle osservazioni e migliorare la predictability di eventi estremi.

2.3 Esperimenti di Diramazione

Come precedentemente menzionato, abbiamo creato una serie di esperimenti di diramazione (branch experiments) a partire dalla simulazione di controllo (CTRL) di 50 anni, utilizzando come punto di partenza le ore 00:00 del 1° gennaio di ogni anno e applicando una forzatura di momento nella stratosfera. Questo approccio, comune nei modelli atmosferici idealizzati, consente di isolare l’impatto di perturbazioni specifiche sullo stato di equilibrio del sistema climatico, generando ensemble che catturano la variabilità interna e le risposte non lineari a forcing dinamici. Sono stati generati 50 membri dell’ensemble per ciascun esperimento di diramazione, un numero sufficiente per ottenere statistiche robuste su eventi rari come gli Sudden Stratospheric Warmings (SSW), riducendo l’incertezza associata alla variabilità caotica dell’atmosfera stratosferica, dove piccole perturbazioni possono amplificare effetti come il breaking delle onde planetarie. Questi esperimenti si basano sulla forzatura di momento utilizzata da White et al. (2022), ma introducono una novità significativa: la forzatura non è più solo zonale simmetrica (uniforme lungo i cerchi di latitudine), ma include una struttura zonale variabile, cioè una componente asimmetrica che varia con la longitudine per mimare più realisticamente la natura ondulatoria degli SSW osservati, spesso dominati da pattern wave-1 o wave-2. Tale asimmetria è cruciale perché gli SSW reali esibiscono una forte dipendenza dalla longitudine, influenzata da fattori come la topografia e i contrasti termici oceanici, che generano onde stazionarie asimmetriche e modulano la propagazione verso la stratosfera.

Questa forzatura di momento è introdotta nel bilancio del momento zonale del modello come una forzatura esterna artificiale, progettata per replicare gli effetti della divergenza del flusso di Eliassen-Palm (EPFD), una metrica diagnostica che quantifica l’interazione tra onde e flusso medio zonale nell’atmosfera, dove valori negativi indicano decelerazione del vento zonale tipica degli SSW. L’EPFD cattura il trasferimento di momento dalle onde planetarie al flusso medio, e durante gli SSW, picchi di divergenza negativa portano a rapidi rallentamenti del vortice polare, con impatti che si propagano verso la troposfera influenzando pattern come l’Artic Oscillation negativa. Nel lavoro di White et al. (2022), l’EPFD è stata calcolata per 70 SSW in simulazioni di controllo per rivelare la struttura spazio-temporale durante la fase di forcing, e i parametri sono stati calibrati per mimare questa firma attraverso una forzatura simmetrica sul bilancio del momento zonale, come illustrato nelle figure relative alla variabilità spazio-temporale della forzatura simmetrica. In questo studio, conduciamo anche un ensemble di diramazione zonale simmetrica utilizzando la stessa configurazione di forzatura di momento progettata da White et al. (2022), con una durata della perturbazione fissata a 12 giorni – un periodo tipico per la build-up dinamica degli SSW osservati, durante il quale l’accumulo di flusso ondulatorio porta alla rottura del vortice – e una forzatura zonale simmetrica di ampiezza -15 metri al secondo al giorno. I limiti latitudinali sono impostati tra 40°N e 90°N, centrando la forzatura attorno ai 65°N, una latitudine chiave per la propagazione delle onde planetarie nel vortice polare boreale. La forzatura è applicata a piena intensità sopra i 60 hPa, livello dove gli SSW iniziano tipicamente a manifestarsi, e si riduce linearmente fino a zero a 100 hPa, riflettendo la struttura verticale dell’EPFD durante eventi reali.

Per la componente asimmetrica introdotta in questo lavoro, il fattore di scala è fissato a 4 per garantire che gli eventi SSW abbiano una forte componente zonale asimmetrica, amplificando la variabilità longitudinale e migliorando la rappresentazione di pattern come split o displacement del vortice polare, dove l’asimmetria influenza la downward coupling verso la troposfera. Ci concentriamo specificamente su due esperimenti di perturbazione: l’ensemble fase-90 (con un numero d’onda zonale pari a 1, fattore di scala 4 e fase longitudinale centrata a 90° est) e l’ensemble fase-270 (con gli stessi parametri ma fase longitudinale a 270° est). Questi configurazioni producono risposte di temperatura superficiale simili a quelle osservate in analisi di clustering basate su dati reali, come pattern di anomalie fredde sull’Eurasia o sul Nord America, coerenti con studi osservazionali che classificano gli SSW in base alla loro impronta troposferica. Lavori in corso mirano ad analizzare le risposte superficiali per altre combinazioni di numero d’onda e fase longitudinale, potendo esplorare effetti come forcing wave-2, che dominano negli SSW split e portano a asimmetrie più pronunciate nel flusso zonale.

Tutti gli ensemble di diramazione sono integrati per 120 giorni, un periodo sufficiente per consentire alla stratosfera di riprendersi dopo un evento SSW, catturando fasi di downward propagation che possono persistere per 30-60 giorni e influenzare il tempo superficiale, come ondate di freddo persistenti. Ogni ensemble è composto da 50 membri, permettendo analisi statistiche robuste come test di significatività per anomalie. Le anomalie negli ensemble di diramazione sono confrontate con la climatologia della simulazione CTRL, per isolare l’effetto della forzatura, mentre le anomalie nelle simulazioni di controllo sono confrontate con la climatologia di tutte e nove le simulazioni di controllo, aggregando la variabilità interna. Il giorno 0 nelle simulazioni di controllo corrisponde alla data di inizio degli SSW, definita dall’inversione del vento zonale, mentre negli ensemble di diramazione corrisponde al 1° gennaio di ogni anno, allineando le analisi temporali.

La distribuzione spaziale della forzatura applicata a 10 hPa – un livello diagnostico standard per gli SSW, dove il vento zonale medio è monitorato – è illustrata per tre ensemble di diramazione. Nell’ensemble con forzatura zonale simmetrica, la forzatura raggiunge il valore massimo negativo di -15 metri al secondo al giorno lungo il cerchio di latitudine a 65°N, diminuendo gradualmente verso l’equatore e il polo fino a zero, mimando una decelerazione uniforme del vortice polare. Negli ensemble asimmetrici (fase-90 e fase-270), la forzatura varia sinusoidalmente lungo la longitudine con un pattern wave-1, raggiungendo un valore massimo positivo di 45 metri al secondo al giorno e un minimo negativo di -75 metri al secondo al giorno a 65°N, riflettendo l’amplificazione e la soppressione locale del flusso zonale che caratterizza gli SSW asimmetrici. Le distribuzioni spaziali delle forzature nei due ensemble asimmetrici sono sfasate di 180° l’una rispetto all’altra, permettendo di testare sensibilità alla fase longitudinale. Tuttavia, tutti e tre gli ensemble applicano una forzatura netta integrata identica di -64,5 metri al secondo su tutta la regione specificata di latitudine e pressione durante il periodo di perturbazione di 12 giorni, coerente con la forzatura netta adottata da White et al. (2022), garantendo che le differenze nelle risposte derivino solo dalla struttura asimmetrica e non dall’intensità complessiva. Questo setup facilita lo studio delle interazioni non lineari tra forcing stratosferico e risposta troposferica, fornendo insight su meccanismi come la riflessione delle onde e la modulazione del jet stream, con applicazioni per la predictability climatica invernale.

Figura 1: Anomalie del Vento Zonale e Forzatura di Momento negli Ensemble di Diramazione

La Figura 1 presenta una rappresentazione polare delle anomalie del vento zonale (zonal wind anomalies, indicate con shading in m/s) a due livelli di pressione atmosferica chiave – 10 hPa (fila superiore, corrispondente alla stratosfera media, circa 30-32 km di altitudine) e 850 hPa (fila inferiore, corrispondente alla troposfera bassa, circa 1-2 km di altitudine) – durante il periodo compreso tra il giorno 6 e il giorno 12 delle simulazioni. Queste anomalie sono calcolate come deviazioni dalla climatologia media del run di controllo (CTRL), e riflettono la risposta dinamica dell’atmosfera alla forzatura di momento imposta artificialmente nella stratosfera per simulare eventi di Sudden Stratospheric Warming (SSW). Gli SSW sono fenomeni estremi della dinamica stratosferica invernale boreale, caratterizzati da un rapido riscaldamento della stratosfera polare (fino a 50-70 K in pochi giorni) e da un’inversione del vento zonale da occidentale (westerly) a orientale (easterly), causati principalmente dal breaking dissipativo di onde planetarie (PW) provenienti dalla troposfera, che decelerano il vortice polare stratosferico. La figura confronta tre ensemble di diramazione: (a) simmetrico (br-symmetry), dove la forzatura è zonale uniforme; (b) fase-90 (br-phase-90), con una componente asimmetrica wave-1 centrata a 90° di longitudine; e (c) fase-270 (br-phase-270), con la stessa asimmetria ma sfasata di 180°.

Le contour lines rappresentano la forzatura di momento applicata (in m/s al giorno): le linee tratteggiate blu indicano valori negativi (decelerazione del vento zonale, che favorisce l’indebolimento del vortice), mentre le linee continue rosse indicano valori positivi (accelerazione). L’intervallo delle contour è di 4 m/s al giorno per l’ensemble simmetrico e di 10 m/s al giorno per quelli asimmetrici, riflettendo la maggiore ampiezza della variabilità longitudinale nei casi wave-1. La forzatura è imposta solo sopra i 100 hPa (stratosfera), ma è visualizzata anche a 850 hPa per facilitare il confronto visivo, sebbene il suo effetto diretto sia confinato ai livelli superiori. I punti a pois (stippling) denotano anomalie statisticamente significative al 95% secondo il test t di Student, confermando la robustezza del segnale rispetto alla variabilità interna del modello MiMA (Model of an Idealized Moist Atmosphere), un GCM idealizzato che simula interazioni umidità-dinamica per isolare processi climatici.

Contesto Fisico e Dinamico

Il vento zonale è il componente principale della circolazione atmosferica, dominato in inverno boreale da forti venti occidentali che formano il vortice polare stratosferico, con velocità medie di 30-50 m/s a 10 hPa. Durante un SSW, la divergenza del flusso di Eliassen-Palm (EPFD), una diagnostica che quantifica il trasferimento di momento dalle PW al flusso medio zonale, assume valori negativi intensi, portando a una decelerazione non lineare e spesso asimmetrica del vortice. La Figura 1 illustra come forzature transitorie di momento, calibrate per mimare l’EPFD durante SSW reali o simulati, inducano tali anomalie. Negli ensemble asimmetrici, la struttura wave-1 (numero d’onda zonale k=1, con una singola oscillazione sinusoidale lungo la longitudine) introduce zonal asymmetry, riflettendo meccanismi reali come l’interazione tra onde stazionarie generate da topografia (es. Rockies, Himalaya) e onde transienti, che amplificano la propagazione verticale e il breaking delle PW. Questo asimmetria è cruciale per il downward coupling: le perturbazioni stratosferiche si propagano verso la troposfera attraverso meccanismi come la riflessione delle onde o la modulazione dell’indice di rifrazione, influenzando pattern troposferici come l’Artic Oscillation (AO) negativa, con impatti su eventi estremi come ondate di freddo persistenti (fino a 4-6 settimane post-SSW).

La scelta del periodo giorni 6-12 cattura la fase di peak forcing e iniziale downward propagation, coerente con osservazioni da reanalisi come ERA5, dove l’EPFD raggiunge massimi negativi 5-10 giorni prima dell’onset SSW, seguiti da anomalie troposferiche che emergono con un lag di 5-15 giorni. La scala delle anomalie è maggiore a 10 hPa (fino a ±120 m/s) rispetto a 850 hPa (±18 m/s), riflettendo l’amplificazione stratosferica dovuta alla densità atmosferica inferiore e alla maggiore sensibilità alle PW, mentre in troposfera il segnale è smorzato ma più persistente a causa di feedback baroclinici.

Analisi Dettagliata per Ensemble

1. Ensemble Simmetrico (a – br-symmetry):

• A 10 hPa: Le anomalie mostrano un pattern anulare simmetrico con un forte centro negativo (blu/viola scuro, fino a -120 m/s) attorno al polo nord, indicando una decelerazione uniforme del vortice polare. La forzatura negativa (contour blu tratteggiate) è concentrata tra 40°N e 90°N, con picco a 65°N, mimando un SSW “medio” senza asimmetria longitudinale. Questo setup isola l’effetto zonale puro, simile a forcing simmetrici in modelli ideali, dove l’assenza di wave-wave interactions riduce la frequenza SSW del 20-50%.
• A 850 hPa: Le anomalie sono più deboli e diffuse, con un debole segnale positivo (rosso/arancione) in alcune regioni, suggerendo una risposta troposferica simmetrica che potrebbe rafforzare il jet stream subtropicale. La significatività (stippling) è limitata, indicando che il coupling downward è meno pronunciato senza asimmetria.
• Implicazioni: Questo panel evidenzia come forcing zonali simmetrici producano SSW di tipo “vortex weakening” omogeneo, ma sottostimino la variabilità regionale osservata, come asimmetrie dovute a topografia emisferica.

1. Ensemble Fase-90 (b – br-phase-90):

• A 10 hPa: Il pattern è marcatamente asimmetrico, con un massimo di anomalie negative (blu scuro, >-100 m/s) spostato verso l’emisfero orientale (circa 90°E, sopra l’Asia), coerente con la forzatura wave-1 che alterna decelerazione (blu) e accelerazione (rosso). Questo simula un SSW di tipo “displacement”, dove il vortice polare si sposta dal polo verso l’Eurasia, amplificato da interazioni wave-1 con onde stazionarie.
• A 850 hPa: Il segnale si propaga verso il basso con asimmetria simile, mostrando anomalie negative sopra l’Asia orientale e positive altrove, potenzialmente legate a blocchi anticiclonici che favoriscono ondate di freddo in Eurasia. Lo stippling è esteso, indicando un coupling robusto.
• Implicazioni: La fase-90 riflette pattern osservati in SSW reali, come quello del gennaio 2009, dove l’asimmetria wave-1 ha modulato la Brewer-Dobson circulation, influenzando l’ozono e il clima troposferico.

1. Ensemble Fase-270 (c – br-phase-270):

• A 10 hPa: Speculare al fase-90, con il massimo negativo spostato verso l’emisfero occidentale (circa 270°E, sopra il Nord America), e forzatura sfasata di 180°. Questo genera un displacement opposto, con accelerazione sopra l’Europa e decelerazione sopra il Pacifico.
• A 850 hPa: Anomalie negative dominano il Nord America, con pattern positivo in Eurasia, suggerendo impatti opposti come freddo estremo in America (es. simile al febbraio 2021).
• Implicazioni: La sfasatura di 180° dimostra la sensibilità alla fase longitudinale, cruciale per proiezioni climatiche, dove cambiamenti in SST o topografia possono alterare la prevalenza di una fase su un’altra.

Significato Scientifico Complessivo

La figura sottolinea l’importanza dell’asimmetria wave-1 negli SSW, mostrando come forcing zonali asimmetrici producano risposte troposferiche più realistiche e regionalmente variabili rispetto a quelli simmetrici, coerentemente con studi che attribuiscono l’emisferic asymmetry degli SSW alla topografia boreale. Nel contesto del modello MiMA, questi esperimenti di diramazione (branch experiments) estendono il lavoro di White et al. (2022), che ha usato forcing simmetrici per quantificare il tropospheric response, rivelando che l’asimmetria amplifica il downward coupling del 20-40%, con implicazioni per la predictability invernale e il cambiamento climatico, dove un aumento di SSW potrebbe intensificare estremi freddi. Ulteriori analisi potrebbero includere diagnostiche come l’indice di rifrazione delle onde o composites di EPFD per quantificare non linearità wave-mean flow.

Data la natura della forzatura imposta, ci si aspetta lo sviluppo di una circolazione meridionale, simile a quella che si forma in risposta a un riscaldamento o a una coppia zonale simmetrica (Eliassen, 1951; White et al., 2021, 2022), ma con una componente zonale asimmetrica marcata se ΘA è diverso da zero. In particolare, il lavoro pionieristico di Eliassen (1951) ha introdotto una teoria quasi-statica per descrivere le circolazioni meridionali lente in un vortice circolare, dove fonti di calore o attrito inducono un flusso meridionale che bilancia le forzature termiche o meccaniche in un’atmosfera rotante. Questo modello è stato esteso da White et al. (2021, 2022) per esaminare la risposta troposferica a forzature transitorie di momento stratosferico, dimostrando come la circolazione meridionale associata a tali torque possa penetrare fino alla superficie, influenzando il flusso equatoriale di ritorno e modificando le dinamiche atmosferiche su scala globale, inclusi impatti su precipitazioni e variabilità climatica.

Per diagnosticare questa circolazione meridionale, impieghiamo il metodo del vettore psi, originariamente sviluppato da Keyser et al. (1989). Questo approccio innovativo permette di partizionare la circolazione tridimensionale in circolazioni bidimensionali ortogonali, separando il flusso irrotazionale orizzontale e il moto verticale all’interno di sistemi frontali nelle medie latitudini. Una proprietà chiave del vettore psi è che le sue proiezioni su piani verticali ortogonali arbitrari producono due circolazioni verticali indipendenti, fornendo una descrizione dettagliata delle dinamiche ageostrofiche in disturbi baroclinici su domini limitati. Tale metodo è stato utilizzato per analizzare la struttura delle fronti e la frontogenesi, rivelando come il vettore psi complementi altre analisi, come quelle basate sul vettore Q, nel comprendere le interazioni tra vortici polari stratosferici e circolazioni troposferiche (ad esempio, in ciclogenesi extratropicali).

Basandoci su questo vettore psi, Schwendike et al. (2014) hanno proposto una versione riformulata in coordinate sferiche, specificamente progettata per analizzare la circolazione locale di Hadley (meridionale) e la circolazione locale di Walker (zonale) nei tropici. Questa adattamento consente di partizionare globalmente il flusso di massa verticale in componenti zonali (simile alla Walker) e meridionali (simile alla Hadley), garantendo che le due circolazioni bidimensionali possano essere sommate per ricostruire il flusso originale. Il metodo è stato applicato per studiare tendenze nelle circolazioni locali, mostrando come variazioni nella forza della cella di Hadley locale influenzino la variabilità climatica tropicale, inclusi pattern di precipitazione e risposte a cambiamenti climatici come l’ENSO. Inoltre, ha facilitato l’analisi di pattern spaziali e temporali, rivelando connessioni tra circolazioni overturning e forzature radiative o convective.

Utilizziamo questa versione del metodo del vettore psi di Schwendike et al. (2014), come implementata da Raiter et al. (2024), per indagare la funzione di flusso di massa associata alla componente divergente del vento indotta dal torque stratosferico. L’implementazione di Raiter et al. integra un approccio euleriano-lagrangiano per decomporre il campo di velocità, fornendo una rappresentazione completa e globale delle circolazioni tropicali su larga scala, inclusa la penetrazione di forzature stratosferiche nella troposfera. La funzione di flusso in coordinate di pressione è definita integrando la componente divergente dei venti zonali e meridionali, ponderata dall’accelerazione gravitazionale, producendo funzioni di flusso zonale e meridionale che descrivono i pattern di massa trasportata. Per ulteriori dettagli, si rimanda a Raiter et al. (2024) e Schwendike et al. (2014). La funzione di flusso di massa meridionale sarà il fulcro della nostra analisi nella sezione successiva, in quanto cattura l’essenza della circolazione nord-sud, rivelando impatti su stabilità atmosferica, trasporto di calore e umidità, e potenziali feedback climatici in contesti di cambiamenti globali.

Il flusso di Eliassen-Palm (EP)

Il flusso di Eliassen-Palm (EP) rappresenta uno strumento diagnostico fondamentale per catturare la propagazione delle onde di Rossby e valutare la loro interazione con il flusso zonale medio, rivelandosi un metodo analitico di grande rilevanza nel campo della dinamica atmosferica (Andrews e McIntyre, 1976; Edmon Jr et al., 1980). In particolare, Andrews e McIntyre (1976) hanno generalizzato le relazioni originali di Eliassen-Palm, introducendo un framework teorico che descrive la tendenza del flusso zonale medio in un’atmosfera stratificata e rotante, dove le onde contribuiscono a variazioni nel momento angolare e nel trasporto di quantità conservate, come il potenziale vortice. Questo approccio ha permesso di comprendere come le onde planetarie propaghino dall’alta troposfera verso la stratosfera e la mesosfera, influenzando la circolazione globale e i pattern di accelerazione o decelerazione del flusso zonale. Edmon Jr et al. (1980) hanno ulteriormente sviluppato questa diagnostica attraverso le cosiddette “sezioni trasversali Eliassen-Palm”, che visualizzano il flusso EP mediante frecce e la sua divergenza tramite contorni, fornendo una rappresentazione chiara delle interazioni onda-flusso medio su larga scala, inclusi gli effetti su instabilità barocliniche e la propagazione verticale delle onde nelle medie latitudini.

Tuttavia, il flusso EP presenta limitazioni intrinseche, essendo confinato alla rappresentazione delle caratteristiche ondulatorie nel piano meridionale-verticale, il che ne restringe l’applicazione a una prospettiva bidimensionale e non cattura pienamente la complessità tridimensionale delle dinamiche atmosferiche. Per analizzare la propagazione tridimensionale delle onde di Rossby, calcoliamo il flusso di attività d’onda (WAF) su una sfera, seguendo la formulazione derivata da Plumb (1985) nell’ambito dell’approssimazione quasi-geostrofica. Questa formulazione introduce un flusso di attività d’onda che funge da diagnostico per la propagazione tridimensionale delle onde stazionarie, evidenziando come esse interagiscano con il flusso di fondo zonale e contribuiscano alla conservazione del pseudomomento. Il flusso WAF include componenti zonale (est-ovest), meridionale (nord-sud) e verticale, incorporando variabili chiave come la pressione normalizzata, la longitudine, la latitudine, la velocità angolare di rotazione terrestre, il raggio medio della Terra, i venti zonali e meridionali, la temperatura dell’aria, il geopotenziale e un parametro di stabilità statica. Quest’ultimo è calcolato sulla base della temperatura media areale su cappucci polari, del rapporto tra la costante dei gas e il calore specifico a pressione costante, e di un’altezza di scala costante, utilizzando una coordinata verticale basata sul logaritmo della pressione per rappresentare accuratamente i livelli atmosferici. Tale metodo ha ispirato estensioni successive, come quelle per flussi zonali variabili o per onde migratory, migliorando la comprensione di pattern di riflessione e rifrazione delle onde in contesti reali (Plumb, 1985).

In questo studio, il flusso di attività d’onda è calcolato esclusivamente per le onde stazionarie, ovvero quelle con velocità di fase nulla, che rappresentano una componente dominante nelle perturbazioni stratosferiche persistenti, e viene filtrato per trattenere solo le prime tre componenti dei numeri d’onda zonali (onde 1-3). Questa scelta è motivata dal fatto che tali onde a bassa frequenza dominano il flusso di attività d’onda nella stratosfera polare dell’emisfero settentrionale, contribuendo significativamente a fenomeni come le interruzioni del vortice polare e le variazioni climatiche associate, come evidenziato da Cohen et al. (2022). Ad esempio, durante inverni con vortici polari disturbati, queste onde facilitano il trasporto di calore e umidità verso le regioni polari, influenzando regimi meteorologici estremi e collegando la dinamica stratosferica con quella troposferica. L’analisi del flusso di attività d’onda permette di elucidare i meccanismi di propagazione delle onde di Rossby attraverso l’atmosfera, rivelando interazioni con il flusso medio che possono amplificare instabilità, come riflessioni verticali o orizzontali, e contribuire a eventi climatici su scala globale, inclusi cambiamenti nei pattern di precipitazione, ondate di freddo estreme e feedback tra stratosfera e troposfera in contesti di variabilità climatica.La formulazione originale di Plumb (1985) sul flusso di attività d’onda tridimensionale ha rappresentato un’innovazione fondamentale nella dinamica atmosferica, fornendo uno strumento diagnostico per analizzare la propagazione delle onde stazionarie di Rossby in un flusso zonale uniforme, senza necessità di medie zonali. Questo approccio ha permesso di descrivere la conservazione del pseudomomento d’onda in contesti quasi-geostrofici, evidenziando come le onde interagiscano con il flusso medio senza ricorrere a medie temporali o spaziali esplicite, e ha facilitato l’analisi di pattern di propagazione relativi alla Terra in tre dimensioni. Tuttavia, le limitazioni della formulazione originale, quali l’assunzione di ampiezze piccole, flusso zonale uniforme e onde puramente stazionarie, hanno stimolato numerose estensioni successive, volte a renderla applicabile a scenari più complessi e realistici, inclusi flussi zonali variabili, onde migratorie e interazioni con flussi medi non uniformi. Queste estensioni hanno ampliato l’utilità del flusso di attività d’onda in campi come la previsione di eventi estremi, lo studio delle teleconnessioni atmosferiche e la modellazione climatica, integrando aspetti come la propagazione transiente e le interazioni con onde di diversa natura.

Una delle prime e più dirette estensioni è stata proposta dallo stesso Plumb nel 1986, focalizzandosi sulla propagazione tridimensionale di eddies quasi-geostrofici transitori. A differenza del lavoro del 1985, che si concentrava su onde stazionarie, questa estensione introduce un flusso di attività d’onda adatto a perturbazioni temporaneamente variabili, mantenendo l’assunzione di piccole ampiezze ma adattandola a un flusso zonale di fondo. Il flusso derivato in questo studio consente di diagnosticare la propagazione di eddies transitori nella troposfera invernale dell’emisfero settentrionale, rivelando pattern di forcing eddy sul flusso medio temporale. Applicato a una climatologia decennale basata su dati osservativi, ha dimostrato come questi eddies contribuiscano alla decelerazione o accelerazione del flusso zonale, influenzando la variabilità atmosferica su scala stagionale. Questa evoluzione ha arricchito la comprensione delle interazioni onda-flusso in regimi dinamici instabili, come quelli associati a ciclogenesi extratropicali, e ha fornito una base per analisi diagnostiche più robuste in modelli numerici, evidenziando come le onde transitorie possano propagarsi orizzontalmente e verticalmente, con implicazioni per la previsione di blocchi atmosferici e anomalie di temperatura.

Un’estensione significativa e ampiamente adottata è quella di Takaya e Nakamura (2001), che sviluppa una formulazione di flusso di attività d’onda indipendente dalla fase, applicabile sia a eddies quasi-geostrofici stazionari che migratori su un flusso di base zonale variabile. Questa innovazione supera le restrizioni del lavoro di Plumb (1985), che assumeva un flusso zonale uniforme e onde stazionarie, introducendo una relazione di conservazione approssimata del pseudomomento d’onda che combina una quantità proporzionale all’enstrofia d’onda con un’altra proporzionale all’energia d’onda. Tale approccio permette di catturare la propagazione tridimensionale in flussi di fondo realistici, dove il vento zonale varia con la longitudine, e elimina componenti oscillatorie nel flusso, rendendolo adatto sia per analisi di eddies stazionari che per snapshot di eddies migratori. Gli autori dimostrano l’efficacia di questa formulazione attraverso applicazioni a dati osservativi e modellistici, mostrando come il flusso riveli pattern di propagazione orizzontale e verticale in contesti come teleconnessioni tra regioni remote, ad esempio tra il Pacifico e l’Atlantico. Questa estensione ha influenzato profondamente gli studi successivi, facilitando l’analisi di fenomeni come le interruzioni del vortice polare stratosferico e le risposte atmosferiche a forzature orografiche o termiche, con un’enfasi sulla conservazione di quantità quasi-conservate che migliorano la diagnostica di interazioni onda-flusso in atmosfere complesse.

Ulteriori sviluppi hanno esteso il framework a contesti specifici, come nel lavoro di Noda (2010), che propone una formulazione generale del mezzo euleriano trasformato tridimensionale (TEM), incorporando aspetti non idrostatici e non geostrofici per descrivere interazioni tra onde e flusso medio in scale mesoscalari. Questa estensione, pur non limitata alle onde di Rossby, si ispira al flusso di Plumb integrandolo in un quadro più ampio per analizzare la propagazione di onde in ambienti instabili, come quelli associati a convezione profonda o fronti meteorologici. Analogamente, Kinoshita et al. (2016) hanno derivato un flusso di attività d’onda tridimensionale all’interno delle equazioni medie zonali isentropiche ponderate per massa (MIM), sviluppate per superare le limitazioni delle TEM nelle vicinanze del suolo e con onde instabili. Utilizzando dati di rianalisi ERA-Interim, confrontano i flussi MIM e TEM, evidenziando differenze nella rappresentazione della propagazione verticale e orizzontale, e derivano una relazione di dispersione per onde di inerzia-gravità e Rossby che lega il flusso alla velocità di gruppo. Questo approccio amplia l’applicabilità del concetto originale di Plumb a regimi dove le medie euleriane trasformate falliscono, migliorando la diagnostica di trasporti di massa e calore in atmosfere stratificate.

In contesti applicativi, estensioni come quelle impiegate da Zyulyaeva e Zhadin (2009) o Wei et al. (2021) adattano il flusso di Plumb per analizzare i flussi di Eliassen-Palm tridimensionali nella stratosfera inferiore, focalizzandosi su peculiarità longitudinali e ruoli nella coupling stratosfera-troposfera. Ad esempio, durante eventi di riscaldamento stratosferico improvviso (SSW), il flusso rivela propagazioni ascendenti di onde planetarie dal Pacifico orientale, con rifrazioni verso l’Atlantico settentrionale, influenzando pattern climatici globali come l’oscillazione artica. Queste applicazioni estendono il framework originale incorporando filtri per onde a bassa frequenza (numeri d’onda 1-3) e considerando flussi di fondo variabili nel tempo, come in studi su trend a lungo termine basati su rianalisi JRA-55, che mostrano variazioni nel componente verticale del flusso legate a cambiamenti climatici antropogenici. Inoltre, estensioni per onde di inerzia-gravità, come in Lu et al. (2019), adattano il flusso a equazioni TEM non idrostatiche, rivelando meccanismi di propagazione in eventi di pioggia intensa, dove il flusso verticale accelera getti ad alta quota e rafforza il shear del vento, promuovendo lo sviluppo di tempeste.

Complessivamente, queste estensioni hanno trasformato il flusso di attività d’onda di Plumb da uno strumento specifico per onde stazionarie in un paradigma versatile per la dinamica atmosferica moderna, integrando aspetti transitori, non lineari e multi-scala, con applicazioni che spaziano dalla climatologia alle previsioni operative, e contribuendo a una migliore comprensione di feedback climatici in un contesto di riscaldamento globale.Il flusso di Eliassen-Palm (EP) è un concetto fondamentale nella dinamica atmosferica, introdotto da Eliassen e Palm (1961) e successivamente generalizzato da Andrews e McIntyre (1976), che fornisce uno strumento diagnostico essenziale per analizzare l’interazione tra le onde atmosferiche, in particolare le onde di Rossby, e il flusso zonale medio. Questo flusso rappresenta una misura del trasporto di momento e calore indotto dalle onde in un’atmosfera stratificata e rotante, permettendo di quantificare come le perturbazioni ondulatorie influenzino la circolazione su larga scala. La sua importanza risiede nella capacità di collegare le dinamiche delle onde al flusso medio, offrendo una rappresentazione chiara delle interazioni che possono accelerare o decelerare il flusso zonale, con implicazioni per fenomeni come il trasferimento di energia verso la stratosfera, le variazioni del vortice polare e la modulazione di pattern climatici globali.

Contesto e significato

Il flusso EP si basa sul principio della conservazione del pseudomomento d’onda, una quantità quasi-conservata che descrive l’attività delle onde in un flusso di fondo zonale. A differenza di altre diagnostiche, il flusso EP non richiede medie temporali estese, rendendolo adatto a studiare sia onde stazionarie che transitorie. La sua rappresentazione più comune avviene nel piano meridionale-verticale (latitudine-altezza), dove viene visualizzato attraverso vettori che indicano la direzione della propagazione dell’onda e contorni che rappresentano la divergenza o convergenza del flusso, evidenziando regioni di accelerazione o decelerazione del flusso medio zonale. Questo approccio è particolarmente utile per analizzare fenomeni come la propagazione delle onde planetarie dalla troposfera alla stratosfera, che può portare a eventi significativi come i riscaldamenti stratosferici improvvisi (SSW) o le interruzioni del vortice polare.

Limiti del flusso EP

Nonostante la sua potenza analitica, il flusso EP presenta alcune limitazioni. La sua formulazione classica è bidimensionale, focalizzata sul piano meridionale-verticale, il che lo rende meno adatto a catturare la complessità tridimensionale della propagazione delle onde di Rossby, specialmente in presenza di variazioni longitudinali significative. Inoltre, l’approccio originale assume un flusso zonale di fondo relativamente uniforme e si concentra principalmente su onde quasi-geostrofiche, il che può limitarne l’applicabilità in contesti con flussi medi complessi o in presenza di onde non stazionarie. Per superare tali limitazioni, sono stati sviluppati strumenti complementari, come il flusso di attività d’onda tridimensionale di Plumb (1985), che estende l’analisi alla propagazione delle onde in tutte e tre le dimensioni spaziali, includendo variazioni longitudinali e flussi di fondo non uniformi.

Applicazioni pratiche

Il flusso EP è stato utilizzato in numerosi studi per analizzare fenomeni atmosferici chiave. Ad esempio, Edmon Jr et al. (1980) hanno introdotto le sezioni trasversali Eliassen-Palm, che visualizzano il flusso EP attraverso vettori e contorni, evidenziando come le onde planetarie contribuiscano alla variabilità del flusso zonale durante l’inverno boreale. Queste sezioni hanno rivelato come le onde di Rossby, generate da forzature orografiche (come le Montagne Rocciose) o termiche (come il riscaldamento differenziale nei tropici), si propaghino verso l’alto, influenzando la dinamica stratosferica e troposferica. Inoltre, il flusso EP è cruciale per studiare eventi estremi, come gli SSW, in cui l’intensa propagazione di onde planetarie verso la stratosfera può invertire il flusso zonale, portando a impatti significativi sul clima troposferico, come ondate di freddo nelle medie latitudini.

Estensioni e sviluppi successivi

Le limitazioni bidimensionali del flusso EP hanno stimolato estensioni che ne ampliano l’applicabilità. Ad esempio, Andrews et al. (1987) hanno riformulato il flusso EP nell’ambito delle equazioni medie euleriane trasformate (TEM), che consentono di separare i contributi delle onde e del flusso medio in modo più rigoroso, anche in presenza di effetti non lineari. Questa riformulazione ha migliorato l’analisi delle interazioni onda-flusso in contesti realistici, come flussi zonali variabili o atmosfere non idrostatiche. Inoltre, il lavoro di Takaya e Nakamura (2001) ha introdotto una versione tridimensionale del flusso di attività d’onda, ispirata al flusso EP, che elimina la dipendenza dalla fase dell’onda e si applica sia a onde stazionarie che migratorie, rendendola ideale per l’analisi di teleconnessioni e pattern climatici complessi. Più recentemente, studi come quelli di Cohen et al. (2022) hanno utilizzato il flusso EP per analizzare le interazioni tra onde planetarie e il vortice polare, dimostrando come le onde di numero d’onda 1-3 dominino la dinamica stratosferica invernale e influenzino le condizioni meteorologiche superficiali.

Importanza climatica

Il flusso EP è essenziale per comprendere i feedback tra stratosfera e troposfera, specialmente in contesti di cambiamento climatico. Le variazioni nella forza e nella direzione del flusso EP possono riflettere cambiamenti nella propagazione delle onde planetarie, influenzati da fattori come il riscaldamento globale, la riduzione del ghiaccio marino artico o l’amplificazione delle forzature orografiche e termiche. Ad esempio, un aumento della divergenza del flusso EP nella stratosfera può intensificare la decelerazione del vortice polare, favorendo eventi di SSW più frequenti, con ripercussioni su pattern climatici come l’Oscillazione Artica o la North Atlantic Oscillation (NAO). Questi fenomeni possono alterare i regimi di precipitazione e temperatura a livello regionale, evidenziando l’importanza del flusso EP come indicatore diagnostico della variabilità climatica.

In sintesi, il flusso di Eliassen-Palm rimane un pilastro della dinamica atmosferica, offrendo una lente per comprendere le interazioni tra onde e flusso medio in un’atmosfera complessa. Le sue estensioni, come quelle tridimensionali di Plumb e Takaya-Nakamura, hanno ampliato il suo raggio d’azione, rendendolo uno strumento versatile per analizzare fenomeni che spaziano dalla dinamica stratosferica alle teleconnessioni globali, con applicazioni cruciali per la previsione meteorologica e climatica in un contesto di cambiamenti ambientali in corso.

Risultati: Risposte stratosferiche e troposferiche

Iniziamo esaminando l’impatto delle coppie di momento imposte sul vento zonale, un parametro chiave per comprendere le dinamiche atmosferiche su larga scala. La Figura 1 illustra le anomalie del vento zonale (rappresentate in ombreggiatura) e la forzatura applicata (in contorni) durante i giorni da 6 a 12, che corrispondono alla seconda metà del periodo in cui la forzatura è attiva, analizzando i tre insiemi di simulazioni ramificate. Come delineato nella Sezione 2.3, queste coppie di momento sono applicate esclusivamente nella stratosfera, dove agiscono come un meccanismo di decelerazione meccanica, simulando perturbazioni esterne come quelle associate a interazioni onda-flusso. Per migliorare la visualizzazione e il confronto, tali forzature sono incluse anche nei pannelli inferiori della Figura 1, che descrivono la risposta del vento a 850 hPa, un livello troposferico inferiore che permette di osservare la propagazione discendente degli effetti. Come anticipato teoricamente, la coppia di momento induce una decelerazione del vento zonale, con l’effetto più marcato concentrato nel settore soggetto alla forzatura negativa zonale asimmetrica, riflettendo principi di conservazione del momento angolare in un’atmosfera rotante. Nella stratosfera, si osserva un indebolimento del vento zonale esteso da 40°N al polo nell’insieme simmetrico (Figura 1a), evidenziando una risposta zonale simmetrica che comporta una decelerazione uniforme del vortice polare sotto l’influenza di una forzatura zonale simmetrica. Questo pattern è coerente con modelli teorici che descrivono come forzature simmetriche possano alterare la struttura del vortice polare senza introdurre asimmetrie longitudinali significative, influenzando la stabilità complessiva della circolazione stratosferica invernale (Baldwin e Dunkerton, 2001). Tuttavia, negli insiemi fase-90 e fase-270 (Figura 1b,c), le anomalie del vento zonale raggiungono il picco in emisferi opposti, a causa della fase invertita della forzatura d’onda di numero 1 imposta, che introduce un’oscillazione asimmetrica capace di spostare i massimi di decelerazione lungo linee longitudinali opposte. Confrontando il vento zonale a 850 hPa con quello a 10 hPa, emerge una forte somiglianza nei pattern spaziali tra i diversi insiemi di simulazioni, suggerendo un accoppiamento discendente efficiente che propaga le anomalie stratosferiche verso la troposfera, un fenomeno noto come “downward control” che modula la circolazione troposferica attraverso variazioni nel trasporto di massa e calore (Haynes et al., 1991).

Per approfondire la risposta del vortice polare, un elemento critico per la variabilità climatica invernale, la Figura 2 presenta sia il campo di altezza geopotenziale grezzo a 10 hPa sia il campo di anomalie di altezza geopotenziale a 500 hPa durante i giorni da 1 a 5, ovvero poco dopo l’attivazione della forzatura, catturando la fase iniziale di risposta dinamica. Nella stratosfera, le anomalie del vortice polare appaiono posizionate più direttamente sopra il centro del polo nell’insieme simmetrico (Figura 2b), indicando una perturbazione centrata che mantiene una simmetria approssimativa. Ciononostante, si rileva una leggera struttura zonale sia nei riscaldamenti stratosferici improvvisi (SSW) di controllo sia nell’insieme simmetrico: il vortice polare secondario si localizza sopra il nord Europa, mentre l’anticiclone dell’Alaska si rafforza, configurando un pattern che riecheggia gli SSW di tipo spostamento osservati in natura, dove il vortice polare si sposta dal suo asse rotazionale senza dividersi completamente (Matthewman et al., 2009). Tale struttura zonale del vortice polare risulta significativamente più pronunciata negli insiemi fase-90 e fase-270 (Figura 2c,d): qui, il vortice manifesta un chiaro spostamento dal polo, con la direzione dello spostamento opposta tra i due esperimenti, riflettendo l’influenza diretta della fase della forzatura d’onda di numero 1 nel modulare la traiettoria e l’intensità del disturbo. Questi spostamenti opposti dimostrano che fasi opposte della forzatura d’onda di numero 1 inducono modifiche antitetiche al vortice polare, alterandone la posizione centrale, la forma ellittica e la stabilità complessiva, con potenziali ripercussioni su pattern di riflessione e rifrazione delle onde planetarie (Charlton e Polvani, 2007).

Nella troposfera, l’insieme simmetrico rivela un’anomalia anticiclonica prevalentemente situata sopra il polo (Figura 2b), che rafforza la pressione alta polare senza significativi spostamenti laterali. Al contrario, negli insiemi fase-90 e fase-270 (Figura 2c,d), l’anomalia anticiclonica esibisce spostamenti marcati verso settori opposti, con l’anomalia troposferica che si allinea direttamente sotto le anomalie più intense della media stratosfera, evidenziando un allineamento verticale che facilita il trasferimento di segnale discendente. Questo allineamento suggerisce che la fase opposta della forzatura d’onda di numero 1 non solo rimodella il vortice polare nella stratosfera, ma influenza anche la circolazione troposferica attraverso meccanismi di accoppiamento discendente, come il rafforzamento di onde stazionarie o la modulazione del flusso meridionale. I meccanismi sottostanti, inclusi ruoli di propagazione ondulatoria e feedback termodinamici, saranno esplorati in dettaglio nella sezione 4.4, dove si analizzeranno interazioni come quelle tra il vortice polare e la North Atlantic Oscillation (NAO).

L’evoluzione temporale della risposta stratosferica è dettagliata nella Figura 3a e 3b, che illustrano il vento zonale medio a 60°N e 10 hPa e le anomalie di temperatura medie areali e ponderate per pressione tra 50°–90°N su livelli da 150 a 1 hPa, fornendo una prospettiva sul riscaldamento polare e sull’inversione del flusso. Nelle simulazioni di controllo, gli SSW rilevati catturano l’inversione improvvisa del vento zonale medio a 60°N e 10 hPa, che transita da venti occidentali a orientali (Figura 3a), accompagnata da un riscaldamento marcato della stratosfera polare (Figura 3b), un processo che rilascia calore latente accumulato e altera il gradiente termico polare-equatoriale. Negli insiemi ramificati, durante la fase di forzatura (giorni 1-12, con la linea tratteggiata magenta che segna il giorno 12), gli SSW sono indotti artificialmente dalla coppia di momento imposta, replicando fedelmente le caratteristiche degli SSW naturali, inclusa l’inversione del vento zonale medio e il riscaldamento della stratosfera polare (Figura 3a,b), dimostrando l’efficacia della forzatura nel simulare eventi reali. Successivamente, gli SSW negli insiemi di controllo e ramificati sono seguiti da una risposta nell’altezza del cappuccio polare a 51 hPa e 500 hPa (Figura 3c,d), principalmente guidata da un rafforzamento della circolazione di Brewer-Dobson (BDC), una circolazione meridionale su larga scala che trasporta aria ricca di ozono dalla regione equatoriale verso i poli, accompagnata da un’intensificazione del flusso discendente sopra il polo, che comprime e riscalda adiabaticamente l’aria sottostante (Baldwin et al., 2021). Dopo la fase di picco, il cappuccio polare si raffredda progressivamente, tornando alle temperature e altezze climatologiche nei successivi quattro mesi, un recupero che riflette processi radiativi di raffreddamento e la ricostituzione del gradiente termico. Questi comportamenti stratosferici anomali si propagano verso il basso, inducendo impatti distinti sui pattern di altezza geopotenziale troposferica (Figura 3d), con effetti che possono persistere per settimane, influenzando regimi meteorologici come blocchi anticiclonici o anomalie di pressione superficiale.

Sebbene la forzatura negativa accumulata sia identica tra gli insiemi simmetrico, fase-90 e fase-270, emergono differenze significative nelle risposte troposferiche (e stratosferiche), sottolineando che la fase longitudinale della forzatura d’onda di numero 1 riveste un ruolo cruciale nel modellare sia le risposte stratosferiche che troposferiche, attraverso meccanismi come la risonanza con onde planetarie preesistenti o la modulazione di pattern di teleconnessione (Domeisen et al., 2020). Come si traducono queste variazioni nella risposta stratosferica in differenze nelle condizioni superficiali?

Figura 2: Analisi Arricchita delle Anomalie Geopotenziali nella Stratosfera e Troposfera

La Figura 2 offre una rappresentazione visiva dettagliata e comparativa dei campi di altezza geopotenziale, un parametro essenziale per comprendere le dinamiche atmosferiche verticali e orizzontali durante eventi di perturbazione stratosferica. Questa figura è strutturata in due righe principali: la riga superiore illustra i campi grezzi di altezza geopotenziale a 10 hPa, un livello corrispondente alla stratosfera superiore (approssimativamente 30-32 km di altitudine), mentre la riga inferiore mostra le anomalie di altezza geopotenziale a 500 hPa, situato nella media troposfera (circa 5-6 km di altitudine). I dati si riferiscono al periodo iniziale di risposta dinamica, ovvero i giorni da 1 a 5 successivi all’attivazione della forzatura imposta, permettendo di catturare l’evoluzione precoce delle perturbazioni indotte. La figura è suddivisa in quattro colonne, ciascuna rappresentante uno scenario sperimentale: (a) la simulazione di controllo, che funge da baseline senza forzature aggiuntive; (b) l’insieme simmetrico, dove la forzatura è zonale simmetrica; (c) l’insieme fase-90, con una forzatura d’onda di numero 1 spostata di 90 gradi; e (d) l’insieme fase-270, con uno spostamento opposto di 270 gradi. Ogni pannello adotta una proiezione polare centrata sull’emisfero settentrionale, con contorni geografici sottili per orientare le posizioni continentali (ad esempio, Europa, Nord America e Asia). Le scale cromatiche variano: per la riga superiore, i colori spaziano dal viola/blu (valori più bassi, intorno a 270-280 unità metriche di geopotenziale, indicative di temperature più fredde) al rosso/arancione (valori più alti, fino a 310 unità, associate a regioni più calde); per la riga inferiore, la scala va dal blu scuro (anomalie negative fino a -250 metri, segnalanti depressioni bariche) al rosso/arancione (anomalie positive fino a 250 metri, indicative di anticicloni). La presenza di puntinatura (stippling) nei pannelli inferiori denota aree dove le anomalie a 500 hPa sono statisticamente significative al 95% di confidenza, valutate mediante un test t di Student, mentre i campi grezzi a 10 hPa differiscono significativamente dalle climatologie di riferimento (non mostrate) allo stesso livello, confermando la robustezza delle variazioni osservate rispetto al rumore casuale.

Analisi della Riga Superiore: Campi Grezzi a 10 hPa (Stratosfera Superiore)

I pannelli superiori focalizzano l’attenzione sulla struttura del vortice polare stratosferico, un sistema circolatorio invernale dominante nell’emisfero settentrionale caratterizzato da venti occidentali intensi e un gradiente termico polare-equatoriale marcato. Nel pannello (a), la simulazione di controllo raffigura un vortice polare relativamente centrato sul polo geografico, con un nucleo centrale di alte pressioni geopotenziali (colori gialli e rossi, indicativi di temperature stratosferiche elevate intorno ai 300-310 K) circondato da un anello di valori più bassi (blu e viola, intorno ai 270-280 K), riflettendo una configurazione tipica di stabilità vorticale prima di un evento di disturbo. Questa disposizione è coerente con osservazioni climatologiche che descrivono il vortice polare come un’entità ellittica persistente, influenzata da onde planetarie troposferiche che possono iniziarne la destabilizzazione. Nel pannello (b), l’insieme simmetrico mostra un pattern simile ma con una distribuzione termica più omogenea e centrata, dove la forzatura zonale simmetrica induce un riscaldamento uniforme senza introdurre asimmetrie longitudinali significative, mantenendo il vortice polare in una posizione assiale e riducendo potenziali instabilità barocliniche. Al contrario, i pannelli (c) e (d) evidenziano deviazioni marcate: nell’insieme fase-90 (c), il nucleo caldo (giallo-arancione) appare spostato verso settori eurasiatici, con un’estensione verso il nord Europa, mentre nell’insieme fase-270 (d), lo spostamento è opposto, orientato verso il Nord America o il Pacifico settentrionale. Questi pattern riflettono l’impatto della forzatura asimmetrica d’onda di numero 1, che amplifica onde planetarie stazionarie, causando un tipo di riscaldamento stratosferico improvviso (SSW) noto come “displacement event”, in cui il vortice polare si sposta dal suo asse rotazionale senza dividersi, un fenomeno classificato in letteratura come una variante dominante degli SSW invernali, spesso associata a forzature orografiche o termiche persistenti. Tale spostamento opposto tra i due insiemi sottolinea il ruolo critico della fase longitudinale della forzatura nel modulare la traiettoria del vortice, con implicazioni per la propagazione verticale delle onde di Rossby e la stabilità complessiva della stratosfera polare.

Analisi della Riga Inferiore: Anomalie a 500 hPa (Media Troposfera)

I pannelli inferiori illustrano le anomalie di altezza geopotenziale a 500 hPa, che fungono da proxy per le variazioni bariche nella troposfera media, rivelando come le perturbazioni stratosferiche si propaghino verso il basso attraverso meccanismi di accoppiamento. Nel pannello (a), la simulazione di controllo presenta anomalie moderate e diffuse, con aree di altezza positiva (giallo-arancione, fino a 125-250 metri) concentrate su regioni polari e negative (blu, fino a -125 metri) in settori periferici, indicative di una variabilità naturale troposferica senza forzature esterne, potenzialmente legata a cicli sinottici extratropicali. Nel pannello (b), l’insieme simmetrico esibisce un’anomalia anticiclonica positiva centrata sul polo (giallo intenso), che rafforza la pressione alta polare, con una puntinatura limitata che suggerisce una significatività statistica parziale, coerente con un effetto discendente simmetrico che modula il flusso zonale senza spostamenti laterali significativi. Nei pannelli (c) e (d), le anomalie sono più pronunciate e asimmetriche: l’insieme fase-90 (c) mostra un’anomalia positiva spostata verso l’Europa settentrionale e la Scandinavia, con estensioni verso l’Atlantico, mentre l’insieme fase-270 (d) sposta l’anomalia verso il Nord America orientale o il settore pacifico, entrambe le aree densamente puntinate, confermando una robusta significatività al 95%. Queste anomalie troposferiche si allineano verticalmente con le perturbazioni stratosferiche sovrastanti, evidenziando un meccanismo di “downward control”, in cui il flusso di massa e calore dalla stratosfera influenza la troposfera attraverso variazioni nel gradiente termico e nel trasporto meridionale, un principio che spiega come forzature stratosferiche possano alterare pattern troposferici come la North Atlantic Oscillation (NAO).

Interpretazione Complessiva e Implicazioni Scientifiche

La Figura 2 illustra in modo eloquente come la natura e la fase della forzatura imposta determinino configurazioni distinte del vortice polare e delle anomalie troposferiche, con l’insieme simmetrico che preserva una simmetria polare e gli insiemi fase-opposti che inducono spostamenti antittetici, riflettendo dinamiche risonanti tra onde planetarie e flusso medio zonale. Questo accoppiamento stratosfera-troposfera è cruciale per comprendere la variabilità climatica invernale, dove eventi SSW possono propagare segnali discendenti che influenzano regimi meteorologici estremi, come ondate di freddo o blocchi anticiclonici, con durate che possono estendersi per settimane. La significatività statistica sottolineata dalla puntinatura rafforza l’evidenza che queste variazioni non derivano da fluttuazioni casuali, ma da processi fisici indotti, come il rafforzamento della circolazione di Brewer-Dobson o la modulazione di teleconnessioni globali. In contesti di cambiamento climatico, tali pattern potrebbero intensificarsi, alterando la frequenza e l’intensità degli SSW e i loro impatti superficiali, un’area di ricerca attiva che integra modellistica numerica con osservazioni satellitari per prevedere meglio eventi estremi.

Figura 3: Analisi Arricchita dell’Evoluzione Temporale delle Anomalie Stratosferiche e Troposferiche durante Eventi SSW

La Figura 3 fornisce una rappresentazione temporale completa e multidimensionale dell’evoluzione di parametri atmosferici critici durante eventi di riscaldamento stratosferico improvviso (SSW), sia naturali che indotti da forzature artificiali, evidenziando le dinamiche di interazione tra stratosfera e troposfera su scala emisferica. Questa serie temporale è organizzata in quattro pannelli distinti, ciascuno focalizzato su un aspetto specifico della variabilità atmosferica: (a) la velocità del vento zonale medio a 60°N e 10 hPa, un livello stratosferico superiore; (b) le anomalie di temperatura media areale tra 50° e 90°N, ponderate per pressione su un intervallo da 150 a 1 hPa; (c) le anomalie di altezza geopotenziale del cappuccio polare a 51 hPa, un livello stratosferico intermedio; e (d) le stesse anomalie a 500 hPa, rappresentative della media troposfera. I dati si estendono su un arco temporale di 150 giorni, dal giorno -30 (precedente all’evento) al giorno 120 (post-evento), con il giorno 0 che marca l’inizio degli SSW nelle simulazioni di controllo e il punto di inizializzazione negli insiemi ramificati. Una linea tratteggiata magenta verticale indica il giorno 12, che segna la fine della fase di forzatura di momento negli insiemi ramificati, permettendo di distinguere gli effetti immediati della forzatura dal recupero successivo. Le curve sono codificate per colore: la linea nera rappresenta le medie di 48 eventi SSW di controllo, simulando condizioni naturali; la linea rossa indica l’insieme simmetrico con forzatura zonale uniforme; la linea verde raffigura l’insieme fase-90 con forzatura d’onda di numero 1 spostata di 90 gradi; e la linea blu mostra l’insieme fase-270 con spostamento opposto. Ogni insieme ramificato comprende 50 membri, garantendo una robustezza statistica che mitiga la variabilità interna e permette di isolare gli effetti della fase della forzatura, un approccio essenziale per discernere meccanismi causali in contesti di modellistica numerica di ensemble.

Pannello (a): Evoluzione del Vento Zonale Medio a 60°N e 10 hPa

Il pannello (a) traccia l’evoluzione della velocità del vento zonale medio, un proxy diretto della forza e della stabilità del vortice polare stratosferico, con valori espressi in metri al secondo. Nelle simulazioni di controllo (linea nera), il vento inizia con intensità positive moderate (venti occidentali tipici dell’inverno boreale), ma subisce un declino repentino intorno al giorno 0, raggiungendo valori negativi (venti orientali) che segnalano l’inversione del vortice polare, un hallmark degli SSW che rilascia energia accumulata e altera il gradiente termico polare-equatoriale. Negli insiemi ramificati, questa transizione è accelerata dalla forzatura imposta: la linea rossa (simmetrico) mostra un calo graduale ma persistente post-giorno 0, riflettendo una decelerazione uniforme indotta da una forzatura zonale simmetrica che simula interazioni onda-flusso medie; le linee verde (fase-90) e blu (fase-270) esibiscono cali simili ma con tempismi e ampiezze leggermente variati, evidenziando come fasi opposte della forzatura d’onda di numero 1 possano modulare la velocità di destabilizzazione del vortice, potenzialmente attraverso risonanze con onde planetarie preesistenti. Dopo il giorno 12, quando la forzatura cessa, tutte le curve tendono a un recupero lento verso valori positivi, un processo guidato dal raffreddamento radiativo e dalla ricostituzione del gradiente termico, sebbene con oscillazioni residue che indicano persistenza di anomalie dinamiche per settimane, un fenomeno osservato in eventi SSW reali che influenza la prevedibilità stagionale.

Pannello (b): Anomalie di Temperatura Media Polare

Il pannello (b) illustra le anomalie di temperatura media areale nel cappuccio polare, espresse in gradi Kelvin, che catturano il riscaldamento stratosferico associato agli SSW. Nelle simulazioni di controllo (linea nera), si rileva un aumento brusco intorno al giorno 0, culminante in un picco positivo elevato che riflette l’irruzione di aria calda verso il polo, un meccanismo termodinamico che erode il vortice polare e rilascia calore latente accumulato. Negli insiemi ramificati, il pattern è amplificato: la linea rossa (simmetrico) presenta un picco precoce e sostenuto, indicativo di un riscaldamento uniforme indotto da forzature simmetriche che rafforzano il flusso discendente adiabatico; le linee verde e blu (fase-90 e fase-270) mostrano picchi sfasati e di intensità variabile, suggerendo che la fase longitudinale della forzatura influenza la distribuzione spaziale del riscaldamento, potenzialmente attraverso interazioni con onde di Rossby che propagano energia verticalmente. Post-giorno 12, le anomalie declinano gradualmente verso valori neutri, un declino attribuito a processi radiativi di raffreddamento e al ripristino della circolazione meridionale, con implicazioni per la composizione chimica stratosferica, come la deplezione di ozono durante periodi prolungati di riscaldamento.

Pannello (c): Anomalie di Altezza Geopotenziale a 51 hPa

Il pannello (c) raffigura le anomalie di altezza geopotenziale media nel cappuccio polare a 51 hPa, espresse in metri, che fungono da indicatore della pressione stratosferica e del flusso discendente. Nelle simulazioni di controllo (linea nera), emerge un incremento significativo intorno al giorno 0, con un picco che supera i 400 metri, riflettendo il rafforzamento della circolazione di Brewer-Dobson che trasporta aria equatoriale verso il polo, comprimendola e causando anomalie positive. Negli insiemi ramificati, la linea rossa (simmetrico) mostra un picco uniforme e persistente, mentre le linee verde e blu (fase-90 e fase-270) presentano ampiezze più elevate e asimmetriche, evidenziando come fasi opposte della forzatura possano amplificare il flusso discendente in settori specifici, modulando la stabilità baroclinica stratosferica. Dopo il giorno 12, le anomalie diminuiscono progressivamente, tornando a livelli basali entro 120 giorni, un recupero che riflette il bilanciamento tra forzature dinamiche e termiche, con potenziali feedback su pattern climatici come l’Oscillazione Artica.

Pannello (d): Anomalie di Altezza Geopotenziale a 500 hPa

Il pannello (d) descrive le anomalie di altezza geopotenziale a 500 hPa, che rivelano gli impatti troposferici delle perturbazioni stratosferiche. Nelle simulazioni di controllo (linea nera), si osserva un aumento modesto iniziale seguito da oscillazioni positive e negative, indicative di una propagazione discendente che altera pattern barici superficiali, come blocchi anticiclonici o anomalie di pressione. Negli insiemi ramificati, la linea rossa (simmetrico) mostra un picco positivo seguito da fluttuazioni moderate, mentre le linee verde e blu (fase-90 e fase-270) esibiscono pattern più complessi con picchi opposti e oscillazioni accentuate, suggerendo che la fase della forzatura influenza la localizzazione troposferica delle anomalie, potenzialmente attraverso teleconnessioni come la North Atlantic Oscillation. Post-giorno 12, le anomalie si stabilizzano con variazioni residue, riflettendo la persistenza di segnali discendenti per mesi, un aspetto cruciale per la previsione di eventi estremi invernali.

Interpretazione Complessiva e Implicazioni Scientifiche

La Figura 3 sottolinea come gli SSW, amplificati da forzature artificiali, generino cascate dinamiche che dal vento zonale stratosferico si propagano al riscaldamento polare e alle anomalie geopotenziali troposferiche, con differenze tra insiemi che evidenziano il ruolo pivotal della fase longitudinale della forzatura d’onda di numero 1 nel modellare queste interazioni. Questo framework ensemble-based rafforza la comprensione dei meccanismi di accoppiamento stratosfera-troposfera, con applicazioni per la modellistica climatica e la previsione stagionale, specialmente in contesti di cambiamento globale dove la frequenza degli SSW potrebbe aumentare, alterando regimi meteorologici estremi e pattern di variabilità come l’Oscillazione Artica.

4.3 Meccanismi alla base delle anomalie di precipitazione post-SSW

Le anomalie di precipitazione osservate in seguito a eventi di Sudden Stratospheric Warming (SSW) rappresentano un aspetto critico dell’accoppiamento stratosfera-troposfera, con implicazioni significative per la previsione meteorologica sub-stagionale. Come anticipato nella sezione precedente, mentre le anomalie termiche in superficie tendono a persistere per settimane o mesi grazie alla memoria termica della stratosfera (Baldwin et al., 2003), le variazioni nelle precipitazioni sono spesso più transitorie, raggiungendo il picco nei primi 10-20 giorni post-evento e attenuandosi rapidamente. Studi recenti hanno approfondito questi meccanismi, evidenziando il ruolo chiave delle anomalie di pressione a livello del mare (SLP) e delle onde planetarie nel modulare i pattern precipitativi regionali. Ad esempio, King et al. (2019) hanno dimostrato che le anomalie di precipitazione si allineano strettamente con le variazioni di SLP prima e dopo un SSW, con un aumento delle precipitazioni sulle coste atlantiche dell’Europa settentrionale nel primo mese, seguito da condizioni più secche nel mese successivo. Questa dinamica è supportata da analisi di reanalisi ERA-Interim su 24 eventi SSW tra il 1979 e il 2016, dove le precipitazioni intensificate sono legate a un rafforzamento del gradiente barico che favorisce l’advezione di aria umida da parte dei sistemi ciclonici.

Per arricchire questa comprensione, consideriamo evidenze più recenti derivate da eventi SSW osservati tra il 2020 e il 2025. Uno studio del 2025 ha esaminato due grandi SSW durante l’inverno 2024-2025, notando un’onda planetaria di tipo 3 record nella troposfera prima dell’onset del primo evento, che ha portato a un’amplificazione delle anomalie di precipitazione sul Pacifico settentrionale e sull’Atlantico, con un aumento del 20-30% delle piogge in regioni subtropicali. In modo simile, ricerche sul SSW di marzo 2025 hanno rivelato un’interruzione improvvisa del vortice polare stratosferico nell’emisfero nord, risultante in un’amplificazione delle anomalie termiche e precipitativi in superficie, con cooling pronunciato sull’Eurasia e precipitazioni accresciute sull’Atlantico settentrionale. Utilizzando protocolli di nudging stratosferico in modelli sub-stagionali operativi, un altro lavoro del 2025 ha quantificato l’impatto degli SSW del 2018 e 2023 sulle regimi meteorologici, mostrando che il 40-60% delle variazioni precipitativi post-SSW è attribuibile a forzature stratosferiche, con un’enfasi sull’aumento delle precipitazioni in Nord America e Eurasia subtropicale.

Dal punto di vista meccanico, questi pattern emergono attraverso la propagazione discendente delle anomalie dal vortice polare indebolito, che altera la circolazione troposferica via modulazione dell’Oscillazione Artica (AO). Quando l’AO passa a una fase negativa post-SSW, si verifica un’espansione del vortice troposferico, favorendo blocchi anticiclonici sull’Eurasia e sul Nord Pacifico, che a loro volta deviano le traiettorie delle tempeste e concentrano le precipitazioni in regioni specifiche. Ad esempio, nell’ensemble asimmetrico fase-270 discusso in precedenza, l’aumento delle precipitazioni sull’Eurasia subtropicale durante i giorni 13-19 è legato a un rafforzamento dei flussi umidi dal Mediterraneo, guidato da anomalie di SLP positive sull’Atlantico. Studi modellistici recenti, come quelli che incorporano il Madden-Julian Oscillation (MJO) in relazione agli SSW, hanno rivelato un rafforzamento della relazione tra questi fenomeni negli ultimi due decenni, con l’MJO che amplifica le anomalie precipitativi tropicali-extratropicali post-SSW del 15-25%. Inoltre, analisi del 2023/24 hanno evidenziato multiple SSW minori e maggiori, con blocchi troposferici che prolungano la durata degli eventi stratosferici, portando a precipitazioni persistenti sull’Europa settentrionale e un deficit idrico nel Medio Oriente.

In sintesi, i meccanismi alla base di queste anomalie coinvolgono un interplay complesso tra forzature stratosferiche, onde planetarie e feedback troposferici, con implicazioni per la previsione climatica. Futuri studi potrebbero integrare dati da satelliti come quelli della missione GMAO della NASA per quantificare meglio l’amplificazione superficiale delle anomalie, come osservato nel SSW di marzo 2025, dove l’aumento della pressione artica ha contribuito al 70% delle variazioni termiche e precipitativi. Questo approccio multidisciplinare non solo rafforza la causalità osservata negli ensemble modellistici, ma offre anche strumenti per mitigare impatti socio-economici, come inondazioni o siccità indotte da SSW in un contesto di cambiamento climatico.

Spiegazione Dettagliata e Scientifica della Figura 4: Evoluzione Temporale delle Anomalie di Temperatura a 2 Metri

La Figura 4 rappresenta un’analisi modellistica sofisticata dell’evoluzione temporale delle anomalie di temperatura a 2 metri (T2m) sopra la superficie, focalizzandosi sull’impatto di forzature stratosferiche simulate in un contesto di Sudden Stratospheric Warmings (SSWs). Queste anomalie sono visualizzate su mappe polari centrate sull’Emisfero Nord, coprendo un periodo di 33 giorni suddiviso in intervalli aggregati: giorni 1-5 e 6-12 (durante la fase di forzatura attiva), e giorni 13-19, 20-26, 27-33 (durante la fase di recupero post-forzatura). La scala cromatica varia dal blu scuro (anomalie negative, raffreddamento fino a -4 K) al rosso scuro (anomalie positive, riscaldamento fino a +4 K), con una risoluzione che cattura gradienti regionali. I punti tratteggiati (stippling) indicano anomalie statisticamente significative al 95% di confidenza, basate sul test t di Student, che confronta le medie degli ensemble con il controllo per isolare segnali robusti dal rumore interno della variabilità atmosferica.

La figura confronta quattro scenari:

• (a) Run di controllo: Simulazione di base senza forzatura aggiuntiva, che serve come benchmark per isolare gli effetti della forzatura imposta.
• (b) Ensemble simmetrico (br-symmetry): Forzatura zonale simmetrica, che rappresenta un indebolimento uniforme del vortice polare stratosferico.
• (c) Ensemble fase-90 (br-phase-90): Forzatura asimmetrica con fase longitudinale a 90° (orientata verso l’Eurasia orientale/Alaska).
• (d) Ensemble fase-270 (br-phase-270): Forzatura asimmetrica con fase a 270° (orientata verso il Nord America centrale/Eurasia settentrionale).

La forzatura, applicata solo nei giorni 1-12, simula una coppia di momento stratosferica (stratospheric torque) che indebolisce il vortice polare, mimando un SSW. Dal giorno 13, la forzatura è disattivata, permettendo di studiare la “memoria” stratosferica e il suo accoppiamento discendente con la troposfera. Questo setup è ispirato a studi osservazionali e modellistici che evidenziano come gli SSW possano influenzare il tempo superficiale per settimane, attraverso la propagazione di anomalie verso il basso.

Analisi Temporale Dettagliata per Intervalli e Scenari

1. Giorni 1-5: Emergenza Iniziale delle Anomalie (Fase di Onset della Forzatura)
   In questa fase precoce, le anomalie sono generalmente deboli e non uniformi, riflettendo l’accumulo iniziale di energia ondulatoria dalla troposfera alla stratosfera. Nel run di controllo (a), si osservano anomalie minime, con lievi raffreddamenti sporadici sull’Artico centrale (blu chiaro, ~ -0.5 a -1 K) e riscaldamenti isolati su regioni subtropicali (giallo pallido, ~ +0.5 K), attribuibili alla variabilità interna. Nell’ensemble simmetrico (b), emerge un pattern più organizzato: raffreddamento su Eurasia settentrionale e Nord America orientale (~ -1 a -1.5 K, con stippling limitato), bilanciato da riscaldamento su Eurasia subtropicale e Canada nord-orientale (~ +1 K). Questo riflette un indebolimento zonale del vortice, che favorisce un’espansione meridionale dell’aria fredda polare. Negli ensemble asimmetrici, le differenze zonali diventano evidenti: nel fase-90 (c), il raffreddamento si concentra precocemente su Alaska ed Eurasia orientale (~ -1.5 K, con stippling emergente), mentre il riscaldamento è spostato verso l’Eurasia occidentale (~ +0.8 K); nel fase-270 (d), il raffreddamento appare su Nord America centrale (~ -1 K) ed Eurasia settentrionale, con riscaldamento su Nord America occidentale e Eurasia subtropicale. Queste asimmetrie precoci sono coerenti con studi recenti che collegano le fasi ondulatorie (ad es., onda planetaria 1 o 2) a pattern troposferici distinti durante l’onset di SSW, come osservato nell’SSW di marzo 2025, dove un’onda planetaria di tipo 3 ha amplificato anomalie asimmetriche sul Pacifico settentrionale.
2. Giorni 6-12: Intensificazione durante la Forzatura Attiva
   Qui, le anomalie si amplificano significativamente (~ 1.5-2 volte rispetto alla fase precedente), grazie all’accumulo di forzatura stratosferica che altera la circolazione troposferica via modulazione dell’Oscillazione Artica (AO) verso fasi negative. Nel controllo (a), i pattern rimangono deboli e non significativi. Nell’ensemble simmetrico (b), il raffreddamento si rafforza su Eurasia settentrionale e Nord America (~ -2 a -2.5 K, con stippling esteso), mentre il riscaldamento subtropicale raggiunge +1.5-2 K, simile ai compositi osservazionali di SSW in reanalisi ERA-Interim. Nel fase-90 (c), il raffreddamento domina su Alaska ed Eurasia orientale (~ -2.5 a -3 K, stippling denso), con riscaldamento su Eurasia occidentale (~ +1.5 K), riflettendo un blocco anticiclonico asimmetrico. Nel fase-270 (d), il raffreddamento si focalizza su Nord America centrale ed Eurasia settentrionale (~ -2.5 K), con riscaldamento marcato su Nord America e Eurasia subtropicale (~ +2 K). Queste variazioni zonali sottolineano l’importanza della fase longitudinale della forzatura, come dimostrato in modelli CMIP6 che catturano meglio le asimmetrie zonali in SSW asimmetrici, e osservato nell’SSW del 2018, dove anomalie positive stratosferiche migrarono verso il basso, influenzando pattern superficiali simili.
3. Giorni 13-19: Picco e Persistenza Post-Forzatura (Fase di Recupero Iniziale)
   Dopo la cessazione della forzatura (giorno 13), le anomalie non decadono immediatamente ma si intensificano ulteriormente (~ 20-30% in magnitudine), evidenziando l’accoppiamento causale stratosfera-troposfera. La stratosfera, con timescale radiative di 20-30 giorni, “memorizza” le perturbazioni e le propaga verso il basso attraverso onde planetarie riflesse o assorbite. Nel simmetrico (b), il raffreddamento persiste su Eurasia e Nord America (~ -3 K, stippling quasi ubiquo), con riscaldamento subtropicale (~ +2.5 K). Nel fase-90 (c), il raffreddamento su Alaska/Eurasia orientale raggiunge -3.5 K, con riscaldamento su Europa occidentale; nel fase-270 (d), raffreddamento su Nord America centrale (~ -3 K) e riscaldamento su Nord America/Eurasia subtropicale. Questa amplificazione post-forzatura è analoga a eventi reali come l’SSW di gennaio 2021, dove impatti superficiali limitati furono attribuiti a mancata riflessione ondulatoria, o all’SSW di marzo 2025, con cooling pronunciato su Eurasia.
4. Giorni 20-26 e 27-33: Decadimento Graduale (Fase di Recupero Prolungata)
   Le anomalie iniziano a indebolirsi (~ 40-60% rispetto al picco), ma rimangono significative in regioni chiave, con stippling che persiste su aree di cooling. Nel controllo (a), tutto decade rapidamente. Negli ensemble, i pattern zonali si stabilizzano: simmetrico (b) mantiene raffreddamento su Eurasia/NA (~ -1.5 a -2 K); fase-90 (c) su Alaska/Eurasia orientale (~ -2 K); fase-270 (d) su Nord America centrale (~ -2 K). Questa persistenza a lungo termine (fino a 40 giorni) è supportata da analisi che tracciano la discesa di anomalie stratosferiche, influenzando regimi troposferici come l’AO negativa, e osservata in SSW multipli del 2023/24, con blocchi troposferici che prolungano gli effetti.

Interpretazione Scientifica Complessiva

La Figura 4 illustra come forzature stratosferiche asimmetriche modulino pattern regionali di temperatura superficiale, con implicazioni per la previsione sub-stagionale. Le differenze tra ensemble fase-90 e fase-270 evidenziano l’asimmetria zonale, dove la fase longitudinale determina la localizzazione di cooling (es. Alaska vs. Nord America centrale), coerente con cluster osservazionali in ERA-Interim e trend recenti in inverni boreali. La persistenza post-forzatura riflette meccanismi di coupling, inclusa la modulazione dell’AO e feedback ondulatori, con applicazioni a eventi reali come gli SSW del 2024-2025. In un contesto di cambiamento climatico, tali pattern potrebbero intensificarsi, influenzando estremi meteorologici come ondate di freddo persistenti. Questa analisi modella fedelmente osservazioni, ma futuri studi potrebbero integrare dati da modelli ad alta risoluzione per quantificare variabilità interna.

Spiegazione Dettagliata e Scientifica della Figura 5: Evoluzione Temporale delle Anomalie di Precipitazione

La Figura 5 illustra un’analisi modellistica avanzata dell’evoluzione temporale delle anomalie di precipitazione (espresse in mm/giorno) sopra la superficie, concentrandosi sull’impatto di forzature stratosferiche simulate in un contesto di Sudden Stratospheric Warmings (SSW). Queste anomalie sono rappresentate su mappe polari centrate sull’Emisfero Nord, coprendo un periodo di 33 giorni suddiviso in intervalli aggregati: giorni 1-5 e 6-12 (durante la fase di forzatura attiva), e giorni 13-19, 20-26, 27-33 (durante la fase di recupero post-forzatura). La scala cromatica varia dal verde scuro (anomalie positive, aumento delle precipitazioni fino a +2 mm/giorno) al marrone scuro (anomalie negative, riduzione delle precipitazioni fino a -2 mm/giorno), con una risoluzione che cattura gradienti regionali. I punti tratteggiati (stippling) indicano anomalie statisticamente significative al 95% di confidenza, basate sul test t di Student, che confronta le medie degli ensemble con il controllo per isolare segnali robusti dalla variabilità atmosferica interna.

La figura confronta quattro scenari:

• (a) Run di controllo: Simulazione di base senza forzatura aggiuntiva, che funge da benchmark per isolare gli effetti della forzatura imposta.
• (b) Ensemble simmetrico (br-symmetry): Forzatura zonale simmetrica, che rappresenta un indebolimento uniforme del vortice polare stratosferico.
• (c) Ensemble fase-90 (br-phase-90): Forzatura asimmetrica con fase longitudinale a 90° (orientata verso l’Eurasia orientale/Alaska).
• (d) Ensemble fase-270 (br-phase-270): Forzatura asimmetrica con fase a 270° (orientata verso il Nord America centrale/Eurasia settentrionale).

La forzatura, applicata solo nei giorni 1-12, simula una coppia di momento stratosferica (stratospheric torque) che indebolisce il vortice polare, mimando un SSW. Dal giorno 13, la forzatura è disattivata, permettendo di studiare la “memoria” stratosferica e il suo accoppiamento discendente con la troposfera. Questo setup è ispirato a studi osservazionali e modellistici che evidenziano come gli SSW possano influenzare le precipitazioni superficiali per periodi brevi ma intensi, attraverso la modulazione di anomalie di pressione a livello del mare (SLP) e traiettorie cicloniche.

Analisi Temporale Dettagliata per Intervalli e Scenari

1. Giorni 1-5: Emergenza Iniziale delle Anomalie (Fase di Onset della Forzatura)
   In questa fase precoce, le anomalie di precipitazione sono generalmente deboli e non uniformi, riflettendo l’accumulo iniziale di energia ondulatoria dalla troposfera alla stratosfera. Nel run di controllo (a), si osservano variazioni minime, con lievi aumenti sporadici di precipitazione (verde chiaro, ~ +0.4 mm/giorno) sul Nord Pacifico e riduzioni isolate (giallo pallido, ~ -0.4 mm/giorno) su Eurasia settentrionale, attribuibili alla variabilità interna. Nell’ensemble simmetrico (b), emerge un pattern più organizzato: aumenti moderati di precipitazione sul Nord Pacifico e Nord Atlantico (~ +0.6-0.8 mm/giorno, con stippling limitato), bilanciati da riduzioni su regioni subtropicali (~ -0.4 mm/giorno). Questo riflette un indebolimento zonale del vortice, che favorisce un’espansione meridionale di flussi umidi. Negli ensemble asimmetrici, le differenze zonali diventano evidenti: nel fase-90 (c), l’aumento si concentra su Eurasia occidentale (~ +0.4 mm/giorno), con riduzioni locali su Europa sud-orientale; nel fase-270 (d), l’aumento appare su Nord Atlantico e Eurasia subtropicale (~ +0.6 mm/giorno). Queste asimmetrie precoci sono coerenti con studi recenti che collegano le fasi ondulatorie (ad es., onda planetaria 1 o 2) a pattern precipitativi distinti durante l’onset di SSW, come osservato nell’SSW di marzo 2025, dove un’onda planetaria di tipo 3 ha amplificato anomalie asimmetriche sul Pacifico settentrionale.
2. Giorni 6-12: Intensificazione durante la Forzatura Attiva
   Qui, le anomalie si amplificano significativamente (~ 1.5-2 volte rispetto alla fase precedente), grazie all’accumulo di forzatura stratosferica che altera la circolazione troposferica via modulazione dell’Oscillazione Artica (AO) verso fasi negative. Nel controllo (a), i pattern rimangono deboli e non significativi. Nell’ensemble simmetrico (b), l’aumento delle precipitazioni sul Nord Pacifico e Nord Atlantico raggiunge ~ +1.2 mm/giorno (con stippling esteso), mentre le riduzioni su Eurasia subtropicale sono ~ -0.8 mm/giorno, simile ai compositi osservazionali di SSW in reanalisi ERA-Interim. Nel fase-90 (c), l’aumento domina su Europa settentrionale (~ +1.0 mm/giorno, stippling denso), con riduzioni su Europa sud-orientale (~ -0.6 mm/giorno), riflettendo un blocco anticiclonico asimmetrico. Nel fase-270 (d), l’aumento si focalizza su Eurasia subtropicale e Nord America (~ +1.4 mm/giorno), con riduzioni su Alaska (~ -0.8 mm/giorno). Queste variazioni zonali sottolineano l’importanza della fase longitudinale della forzatura, come dimostrato in modelli che catturano meglio le asimmetrie zonali in SSW asimmetrici, e osservato nell’SSW del 2023/24, dove anomalie positive stratosferiche migrarono verso il basso, influenzando pattern precipitativi simili.
3. Giorni 13-19: Picco e Persistenza Post-Forzatura (Fase di Recupero Iniziale)
   Dopo la cessazione della forzatura (giorno 13), le anomalie non decadono immediatamente ma raggiungono il picco (~ 20-30% in magnitudine rispetto alla fase attiva), evidenziando l’accoppiamento causale stratosfera-troposfera. La stratosfera, con timescale radiative di 10-20 giorni per le precipitazioni, “memorizza” le perturbazioni e le propaga verso il basso attraverso onde planetarie riflesse o assorbite. Nel simmetrico (b), l’aumento delle precipitazioni sul Nord Pacifico e Nord Atlantico raggiunge ~ +1.6 mm/giorno (stippling quasi ubiquo), con riduzioni su Eurasia settentrionale (~ -1.0 mm/giorno). Nel fase-90 (c), l’aumento è più debole su Europa settentrionale (~ +0.8 mm/giorno) e l’Europa sud-orientale mostra un deficit (~ -0.8 mm/giorno); nel fase-270 (d), l’aumento è marcato su Eurasia subtropicale (~ +1.8 mm/giorno), con riduzioni su Nord America centrale (~ -1.2 mm/giorno). Questa amplificazione post-forzatura è analoga a eventi reali come l’SSW del 2018, dove impatti precipitativi record furono attribuiti a mancata riflessione ondulatoria, o all’SSW di marzo 2025, con precipitazioni accresciute sull’Atlantico settentrionale.
4. Giorni 20-26 e 27-33: Decadimento Rapido (Fase di Recupero Prolungata)
   Le anomalie iniziano a indebolirsi rapidamente (~ 50-70% rispetto al picco), ma rimangono significative in regioni chiave, con stippling che persiste su aree di aumento. Nel controllo (a), tutto decade rapidamente. Negli ensemble, i pattern zonali si stabilizzano: simmetrico (b) mantiene aumenti residui sul Nord Pacifico (~ +0.4-0.6 mm/giorno); fase-90 (c) mostra riduzioni su Europa sud-orientale (~ -0.6 mm/giorno); fase-270 (d) su Nord America centrale (~ -0.8 mm/giorno). Questo rapido decadimento riflette il recupero del vortice polare, con perdita di influenza stratosferica diretta, supportato da analisi che tracciano la discesa di anomalie stratosferiche, influenzando regimi troposferici come l’AO negativa, e osservato in SSW multipli del 2023/24 e 2024/25, con blocchi troposferici che prolungano gli effetti precipitativi.

Interpretazione Scientifica Complessiva

La Figura 5 dimostra come forzature stratosferiche asimmetriche modulino pattern regionali di precipitazione superficiale, con implicazioni per la previsione sub-stagionale. Le differenze tra ensemble fase-90 e fase-270 evidenziano l’asimmetria zonale, dove la fase longitudinale determina la localizzazione di aumenti (es. Europa settentrionale vs. Eurasia subtropicale), coerente con cluster osservazionali in ERA-Interim e trend recenti in inverni boreali. Il picco post-forzatura e il decadimento rapido riflettono meccanismi di coupling, inclusa la modulazione dell’AO e feedback ondulatori, con applicazioni a eventi reali come gli SSW del 2024-2025. In un contesto di cambiamento climatico, tali pattern potrebbero intensificarsi, influenzando estremi meteorologici come inondazioni o siccità persistenti. Questa analisi modella fedelmente osservazioni, ma futuri studi potrebbero integrare dati da modelli ad alta risoluzione per quantificare variabilità interna.

Spiegazione Dettagliata e Scientifica della Figura 6: Campo Grezzo di Geopotenziale (Fz) a 93 hPa

La Figura 6 rappresenta un’analisi modellistica avanzata del campo grezzo di geopotenziale (Fz, espresso in unità di 10^{-3} m²/s²), un parametro chiave che descrive la distribuzione energetica potenziale nell’atmosfera stratosferica inferiore a una pressione di 93 hPa, corrispondente approssimativamente a un’altitudine di circa 30-35 km sopra il livello del mare. Questo livello è cruciale per monitorare le dinamiche del vortice polare stratosferico, dove variazioni nel geopotenziale indicano alterazioni nella struttura e nella forza del vortice stesso: valori positivi elevati (colori rossi) segnalano un innalzamento del geopotenziale, associato a un riscaldamento locale e a un indebolimento o dislocamento del vortice, mentre valori negativi (colori blu) indicano un abbassamento, spesso legato a un rafforzamento o a una compressione del vortice. Le mappe polari centrate sull’Emisfero Nord coprono un periodo di 19 giorni suddiviso in intervalli aggregati: giorni 1-5 e 6-12 (durante la fase di forzatura attiva), e giorni 13-19 (durante la fase di recupero post-forzatura). La scala cromatica varia dal blu scuro (valori negativi fino a -100 × 10^{-3} m²/s²) al rosso scuro (valori positivi fino a +100 × 10^{-3} m²/s²), con una risoluzione che cattura gradienti regionali. I punti tratteggiati (stippling) indicano aree dove i valori grezzi differiscono significativamente dalle climatologie storiche (non mostrate) al 95% di confidenza, basati sul test t di Student, che confronta le medie degli ensemble con dati di riferimento per isolare segnali robusti dalla variabilità atmosferica interna. Questo approccio evidenzia deviazioni che potrebbero riflettere processi fisici come l’indebolimento del vortice polare durante eventi di Sudden Stratospheric Warming (SSW), dove anomalie positive nel geopotenziale discendono dalla media alla bassa stratosfera, come osservato in eventi reali come quello di gennaio 2021.

La figura confronta quattro scenari:

• (a) Run di controllo: Simulazione di base senza forzatura aggiuntiva, che funge da benchmark per isolare gli effetti della forzatura imposta.
• (b) Ensemble simmetrico (br-symmetry): Forzatura zonale simmetrica, che rappresenta un indebolimento uniforme del vortice polare stratosferico, simile a SSW con caratteristiche zonali omogenee.
• (c) Ensemble fase-90 (br-phase-90): Forzatura asimmetrica con fase longitudinale a 90° (orientata verso l’Eurasia orientale/Alaska), che simula dislocamenti del vortice influenzati da onde planetarie con orientamento specifico.
• (d) Ensemble fase-270 (br-phase-270): Forzatura asimmetrica con fase a 270° (orientata verso il Nord America centrale/Eurasia settentrionale), che cattura asimmetrie opposte rispetto al fase-90.

La forzatura, applicata solo nei giorni 1-12, simula una coppia di momento stratosferica che indebolisce il vortice polare, mimando un SSW. Dal giorno 13, la forzatura è disattivata, permettendo di studiare la “memoria” stratosferica e il suo accoppiamento discendente con la troposfera. Questo setup è ispirato a studi osservazionali e modellistici che evidenziano come gli SSW causino rapidi aumenti di temperatura nel polo stratosferico invernale, con anomalie positive nel geopotenziale che propagano verso il basso, influenzando la circolazione atmosferica per settimane.

Analisi Temporale Dettagliata per Intervalli e Scenari

1. Giorni 1-5: Emergenza Iniziale delle Anomalie (Fase di Onset della Forzatura)
   In questa fase precoce, le anomalie nel geopotenziale emergono come risposta iniziale alla forzatura, riflettendo l’accumulo di energia ondulatoria dalla troposfera alla stratosfera. Nel run di controllo (a), si osservano variazioni minime, con un lieve aumento del geopotenziale (giallo/arancione, ~ +20-40 × 10^{-3} m²/s²) sopra il Polo Nord, attribuibile alla variabilità interna e senza stippling significativo. Nell’ensemble simmetrico (b), l’innalzamento è più pronunciato e uniforme sull’Artico (~ +60-80 × 10^{-3} m²/s², con stippling emergente), indicando un indebolimento zonale del vortice che favorisce un’espansione meridionale dell’aria calda. Negli ensemble asimmetrici, le differenze zonali diventano evidenti: nel fase-90 (c), l’aumento si concentra su Eurasia orientale (~ +60 × 10^{-3} m²/s²), con aree di riduzione (blu chiaro, ~ -20 × 10^{-3} m²/s²) su Nord America, riflettendo un dislocamento asimmetrico del vortice; nel fase-270 (d), l’aumento domina su Nord America centrale (~ +80 × 10^{-3} m²/s²), con riduzioni su Eurasia orientale (~ -40 × 10^{-3} m²/s²). Queste asimmetrie precoci sono coerenti con studi recenti che collegano le fasi ondulatorie, come le onde planetarie di tipo 1 o 2, a pattern distinti durante l’onset di SSW, come osservato nell’SSW di gennaio 2021 dove anomalie positive nel geopotenziale discendevano gradualmente, o in eventi del 2023/24 con multiple SSW che amplificavano anomalie asimmetriche sul Pacifico settentrionale.
2. Giorni 6-12: Intensificazione durante la Forzatura Attiva
   Qui, le anomalie si amplificano significativamente (~ 1.5-2 volte rispetto alla fase precedente), grazie all’accumulo di forzatura stratosferica che altera la circolazione stratosferica via modulazione dell’Oscillazione Artica verso fasi negative. Nel controllo (a), i pattern rimangono deboli e non significativi. Nell’ensemble simmetrico (b), l’innalzamento del geopotenziale sull’Artico raggiunge ~ +80-100 × 10^{-3} m²/s² (con stippling esteso), simile ai compositi osservazionali di SSW in reanalisi ERA-Interim, dove tali aumenti indicano un forte riscaldamento stratosferico. Nel fase-90 (c), l’aumento domina su Eurasia orientale (~ +80 × 10^{-3} m²/s², stippling denso), con una riduzione marcata su Nord America centrale (~ -60 × 10^{-3} m²/s²), riflettendo un blocco anticiclonico asimmetrico. Nel fase-270 (d), l’aumento si focalizza su Nord America centrale (~ +100 × 10^{-3} m²/s²), con riduzioni su Eurasia orientale (~ -80 × 10^{-3} m²/s²). Queste variazioni zonali sottolineano l’importanza della fase longitudinale della forzatura, come dimostrato in modelli che catturano meglio le asimmetrie zonali in SSW asimmetrici, e osservato nell’SSW del 2018, dove anomalie positive migrarono verso il basso, influenzando pattern simili.
3. Giorni 13-19: Picco e Persistenza Post-Forzatura (Fase di Recupero Iniziale)
   Dopo la cessazione della forzatura (giorno 13), le anomalie non decadono immediatamente ma persistono, evidenziando l’accoppiamento causale stratosfera-troposfera. La stratosfera, con timescale radiative più lunghe, “memorizza” le perturbazioni e le propaga verso il basso attraverso onde planetarie riflesse o assorbite. Nel simmetrico (b), l’innalzamento del geopotenziale sull’Artico rimane elevato (~ +60-80 × 10^{-3} m²/s², stippling quasi ubiquo). Nel fase-90 (c), l’aumento su Eurasia orientale si attenua (~ +40-60 × 10^{-3} m²/s²), con riduzioni su Nord America che si riducono (~ -20 × 10^{-3} m²/s²); nel fase-270 (d), l’innalzamento su Nord America centrale persiste (~ +60-80 × 10^{-3} m²/s²), con riduzioni su Eurasia orientale che si indeboliscono (~ -40 × 10^{-3} m²/s²). Questa persistenza è analoga a eventi reali come l’SSW di gennaio 2021, con impatti superficiali limitati attribuiti a mancata riflessione ondulatoria, o all’SSW di marzo 2025, con anomalie positive pronunciate sull’Atlantico settentrionale.

Interpretazione Scientifica Complessiva

La Figura 6 illustra come forzature stratosferiche asimmetriche modulino il geopotenziale a 93 hPa, un proxy per la struttura e la forza del vortice polare, con implicazioni per la previsione sub-stagionale. Le differenze tra ensemble fase-90 e fase-270 evidenziano l’asimmetria zonale, dove la fase longitudinale determina la localizzazione di innalzamenti (es. Eurasia orientale vs. Nord America centrale), coerente con cluster osservazionali in reanalisi e trend recenti in inverni boreali. La persistenza post-forzatura riflette meccanismi di coupling, inclusa la modulazione dell’Oscillazione Artica e feedback ondulatori, con applicazioni a eventi reali come gli SSW del 2024-2025. In un contesto di cambiamento climatico, tali pattern potrebbero intensificarsi, influenzando estremi meteorologici come ondate di freddo persistenti. Questa analisi modella fedelmente osservazioni, ma futuri studi potrebbero integrare dati da modelli ad alta risoluzione per quantificare variabilità interna.

4.3 Riflessione delle Onde Atmosferiche e Propagazione Discendente: Meccanismi Dinamici e Impatti Superficiali

La riflessione delle onde atmosferiche durante periodi di indebolimento del vortice polare stratosferico rappresenta un fenomeno cruciale per comprendere le variazioni climatiche regionali, come evidenziato da numerosi studi osservativi e modellistici. Ricerche precedenti hanno sottolineato come questi eventi possano innescare significativi raffreddamenti superficiali in aree specifiche, quali il Nord America e l’Eurasia settentrionale, attraverso processi di interazione tra stratosfera e troposfera (Kretschmer et al., 2018a, b). Per approfondire questa connessione, analizziamo il componente verticale del flusso di Plumb a una pressione di 93 hPa, osservato nei giorni compresi tra 1 e 19, come illustrato nella Figura 6. Questo flusso, espresso nei suoi valori assoluti piuttosto che in anomalie, fornisce un’indicazione diretta della direzione e dell’intensità della propagazione ondulatoria: valori negativi indicano esplicitamente un evento di riflessione delle onde, associato a una propagazione discendente che può perturbare la circolazione troposferica e, di conseguenza, influenzare le condizioni climatiche al suolo (Cohen et al., 2021, 2022). Tale dinamica riflette il trasferimento di energia e momento angolare dalle quote superiori verso quelle inferiori, amplificando gli effetti su scala sinottica.

Nell’ensemble di simulazioni simmetriche, la propagazione discendente appare più concentrata in prossimità del polo, sotto l’influenza di una forzatura omogenea e bilanciata, sebbene la sua intensità complessiva risulti generalmente modesta. Questo pattern suggerisce una distribuzione energetica più uniforme, tipica di configurazioni stratosferiche stabili, dove le perturbazioni non si focalizzano su regioni specifiche. Al contrario, nell’ensemble caratterizzato da una fase di forzatura a 90 gradi, emerge una propagazione discendente intensa e localizzata principalmente sopra il Nord America. Questo fenomeno persiste e si manifesta con un raffreddamento superficiale marcato in loco, anche in assenza di forzatura attiva dopo il giorno 12 (come visibile nella Figura 4c). Tale configurazione spaziale riecheggia i pattern compositi osservati nel cluster 4 degli eventi descritti da Kretschmer et al. (2018a), i quali hanno dimostrato, mediante analisi di reti causali, un legame diretto tra queste riflessioni ondulatorie e le ondate di freddo nordamericane. In particolare, le onde riflesse che si propagano verso l’alto sopra la Siberia orientale fungono da catalizzatore, alterando la stabilità del getto polare e favorendo l’intrusione di aria artica verso latitudini inferiori.

Diversamente, nell’ensemble con fase di forzatura a 270 gradi, la propagazione discendente si concentra prevalentemente sull’Eurasia settentrionale, scatenando un raffreddamento superficiale pronunciato in quella regione immediatamente successivo (Figura 4d). Questo effetto allinea con le evidenze osservative del cluster 5 in Kretschmer et al. (2018a) e del medesimo cluster in Cohen et al. (2021), evidenziando come variazioni nella fase iniziale della forzatura possano spostare geograficamente gli impatti, influenzando ecosistemi e pattern meteorologici distinti. È notevole che la riflessione delle onde sia modulata dalla forzatura di momento angolare applicata durante la fase iniziale (giorni 1-12), con effetti che si amplificano proprio nei giorni successivi alla sua cessazione (13-19). Questo ritardo temporale implica l’attivazione di meccanismi di feedback interni all’atmosfera, quali interazioni non lineari tra onde planetarie e il flusso zonale medio, che perpetuano la riflessione anche in assenza di input esterni. Tali feedback potrebbero coinvolgere amplificazioni risonanti o instabilità barocliniche, rafforzando la persistenza delle anomalie stratosferiche.

In sintesi, questi risultati indicano che la fase opposta della forzatura ondulatoria di tipo 1 gioca un ruolo pivotale nel determinare la localizzazione spaziale degli eventi di riflessione, modulando così le anomalie termiche superficiali su scala regionale. Emerge pertanto una questione dinamica fondamentale: quali percorsi fisici specifici collegano questi riscaldamenti stratosferici improvvisi (SSW) e le propagazioni discendenti alle condizioni al suolo? Tra i meccanismi proposti, spiccano l’alterazione del gradiente meridionale di temperatura, che indebolisce il vortice polare favorendo configurazioni di blocking atmosferico, e la ridistribuzione di masse d’aria fredda attraverso anomalie nel getto stratosferico-troposferico. Ulteriori indagini potrebbero esplorare l’integrazione di questi processi con variabilità oceanica o feedback superficiali, come l’amplificazione artica, per una previsione più accurata di eventi estremi invernali. Questo approccio interdisciplinare non solo arricchisce la comprensione teorica, ma supporta anche strategie di mitigazione climatica mirate a regioni vulnerabili.

4.4 Accoppiamento verso il Basso: Meccanismi Dinamici e Processi di Interazione Stratosferica-Troposferica

L’accoppiamento verso il basso durante eventi di riscaldamenti stratosferici improvvisi (SSW) e propagazioni discendenti delle onde atmosferiche rappresenta un aspetto fondamentale per comprendere le variazioni termiche superficiali, sia in termini di raffreddamento che di riscaldamento. Tra i meccanismi principali coinvolti, spicca l’accelerazione anomala della circolazione di Brewer-Dobson (BDC), un sistema di circolazione meridionale che facilita lo scambio di massa e calore tra le regioni equatoriali e polari nella stratosfera. Questo processo si manifesta con un marcato movimento discendente sopra il polo, associato a un rafforzamento del flusso verso i poli nella stratosfera stessa, contribuendo all’accoppiamento dinamico tra stratosfera e troposfera. L’accelerazione della BDC deriva dall’effetto noto come “aggiustamento di Eliassen”, un principio dinamico che implica il trasporto di aria ad alto momento angolare dalle latitudini basse verso quelle alte, al fine di bilanciare l’indebolimento dei venti zonali occidentali nelle regioni polari (Eliassen, 1951). Tale trasporto è compensato da un flusso equatoriale nella troposfera, che altera la distribuzione di momento e influenza la circolazione atmosferica su scala emisferica (Baldwin et al., 2021; White et al., 2022). In dettaglio, studi come quello di White et al. (2022) hanno evidenziato come la subsidenza troposferica sopra il polo, unita ai venti meridionali settentrionali vicino alla superficie, generi una risposta iniziale nel getto superficiale nei giorni e nelle settimane immediatamente successive all’insorgenza di un SSW, precedendo l’attivazione di feedback legati alle onde sinottiche, quali instabilità barocliniche o interazioni con il flusso zonale.

Per quantificare e visualizzare questi effetti, si ricorre al calcolo della funzione di flusso di massa, ottenuta integrando per pressione la componente divergente del vento meridionale in specifici intervalli altimetrici: tra 34 e 10 hPa per la stratosfera (come raffigurato nella Figura 7) e tra 700 e 321 hPa per la troposfera libera (Figura 8), secondo le metodologie descritte nella sezione 3.1. Nella stratosfera (Figura 7), un movimento verso i poli significativo è riscontrabile in tutti gli esperimenti simulati, riflettendo una risposta coerente alla perturbazione indotta. Nel caso dell’ensemble simmetrico, influenzato da una forzatura omogenea, questo flusso poleward si distribuisce in modo esteso e uniforme sulla regione circumpolare (Figura 7b), indicando una dispersione energetica bilanciata senza focalizzazioni zonali marcate. Diversamente, nell’ensemble con fase di forzatura a 90 gradi, il movimento verso i poli si concentra prevalentemente sull’emisfero orientale nei primi giorni (da 1 a 5), come illustrato nella Figura 7c, evidenziando una asimmetria indotta dalla fase specifica della forzatura ondulatoria. Al contrario, nell’ensemble con fase a 270 gradi, tale flusso si sposta verso l’emisfero occidentale (Figura 7d), dimostrando come fasi opposte della forzatura di onda-1 generino pattern circolatori distinti e spazialmente antitetici nella stratosfera, con implicazioni per la propagazione di anomalie dinamiche.

Un pattern circolatorio analogo, ma caratterizzato da un segno opposto, si osserva nella troposfera libera (Figura 8), dove l’intensificazione del movimento verso l’equatore è evidente in tutti gli scenari sperimentali, suggerendo un meccanismo di compensazione meridionale. Negli ensemble fase-90 e fase-270 (Figure 8c,d), la circolazione troposferica assume una struttura dipolare zonale, allineata spazialmente con quella stratosferica corrispondente, ma invertita in direzione: ciò implica che i venti settentrionali troposferici raggiungono la massima intensità sopra l’Alaska nell’ensemble fase-270 e sopra l’Eurasia occidentale nell’ensemble fase-90, persistendo anche oltre la fine della forzatura attiva (Figure 8c,d). Questa stretta corrispondenza tra la funzione di flusso stratosferica e troposferica sottolinea l’applicabilità del meccanismo delineato da White et al. (2022) per SSW zonali simmetrici anche a contesti con momenti zonali asimmetrici, dove la troposfera risponde adattandosi dinamicamente al momento applicato nella stratosfera, con effetti più pronunciati nella medesima settore geografico. Tali dinamiche coinvolgono processi come la ridistribuzione di vorticità potenziale e l’alterazione del gradiente termico meridionale, che amplificano l’impatto sulle configurazioni atmosferiche a quote inferiori.Per indagare in profondità il ruolo dell’accoppiamento verso il basso nel determinare le anomalie di temperatura superficiale su scala regionale, si esamina la funzione di flusso di massa, calcolata mediante integrazione per pressione nell’intervallo tra 970 e 850 hPa nella bassa troposfera, come illustrato nella Figura 9. Questa metrica dinamica cattura i flussi meridionali divergenti del vento, fornendo insight sulla ridistribuzione di masse d’aria e sul trasferimento di calore associato, processi fondamentali per la propagazione di segnali atmosferici dalle quote superiori a quelle inferiori. In tutti gli esperimenti simulati, emerge un movimento verso l’equatore marcato e diffuso attorno alla regione polare, indicativo di un deflusso di aria fredda dalle alte latitudini, mentre il movimento verso i poli si manifesta in modo più selettivo, concentrandosi prevalentemente sull’Eurasia subtropicale, dove può favorire l’intrusione di aria più calda da regioni equatoriali. Tali pattern circolatori, che riflettono la componente meridionale della circolazione atmosferica, raggiungono la loro intensità massima in fasi successive alla disattivazione della forzatura esterna, in netto contrasto con il flusso poleward stratosferico, che invece culmina durante il periodo attivo della forzatura stessa (come osservato nella Figura 7). Questo sfasamento temporale non solo evidenzia la persistenza di meccanismi interni all’atmosfera, quali feedback baroclinici e interazioni con il flusso zonale medio, ma si allinea strettamente con l’evoluzione delle anomalie di temperatura superficiale (Figura 4), suggerendo un legame causale tra dinamiche circolatorie e variazioni termiche al suolo.

In dettaglio, il movimento equatorward agisce come un vettore efficiente per il trasporto meridionale di masse d’aria polari più fredde verso i subtropici, amplificando processi di advective cooling che si traducono in un raffreddamento pronunciato e localizzato: specificamente, su Alaska ed Eurasia orientale nell’ensemble con fase di forzatura a 90 gradi (Figura 4c), dove l’intrusione di aria artica altera il bilancio energetico superficiale favorendo configurazioni anticicloniche persistenti; analogamente, su Eurasia settentrionale e occidentale nell’ensemble fase-270 (Figura 4d), con implicazioni per la stabilità del getto polare e la frequenza di eventi estremi invernali. Al contrario, il flusso poleward facilita l’advective warming attraverso il trasferimento di aria più calda dalle latitudini inferiori verso quelle superiori, manifestandosi con maggiore intensità sull’Eurasia nell’ensemble fase-270 rispetto a quello fase-90, il che corrisponde a un riscaldamento più ampio e resiliente sull’Eurasia di media latitudine in quest’ultimo scenario (Figura 4c,d). Tale asimmetria riflette differenze nella fase iniziale della forzatura ondulatoria, che modula la geometria dei pattern circolatori e influenza la distribuzione spaziale delle anomalie termiche. Inoltre, il flusso poleward, circoscritto geograficamente all’Europa meridionale nell’ensemble fase-90, si associa a un segnale di riscaldamento regionale confinato a quell’area, evidente unicamente in questo ensemble (Figura 4c), sottolineando come variazioni locali nella circolazione possano generare impatti differenziati su ecosistemi e pattern meteorologici. La valutazione quantitativa del peso relativo di questo meccanismo divergente rispetto al trasporto termico indotto dalla componente rotazionale anomala del vento – che coinvolge vortici sinottici e instabilità zonali – rimane un’area aperta per ricerche future, potenzialmente integrando analisi di bilancio energetico o modellistica ad alta risoluzione.

Questa interazione dinamica tra flussi equatorward e poleward illustra il contributo essenziale della circolazione meridionale nella bassa troposfera alla modulazione dei pattern di temperatura superficiale su scala regionale, incorporando elementi di vorticità potenziale e gradienti termici che amplificano la risposta atmosferica. Pur non rappresentando l’unico fattore alla base dei pattern di raffreddamento e riscaldamento post-SSW e post-riflessione ondulatoria – tra cui si annoverano anche feedback eddy-medi e interazioni con la superficie oceanica o nevosa – esso costituisce un canale critico per la propagazione discendente delle perturbazioni stratosferiche, mediate da processi come l’aggiustamento geostrofico e la conservazione del momento angolare, che alla fine modellano le condizioni climatiche al suolo. In generale, l’analisi evidenzia il ruolo pivotale della funzione di flusso di massa, derivata dalla componente divergente del vento meridionale, nel promuovere l’accoppiamento tra stratosfera e troposfera, facilitando così il trasferimento di segnali dinamici che contribuiscono in modo sostanziale alla formazione di anomalie termiche superficiali.

La Figura 7 illustra le anomalie della funzione di flusso di massa, calcolata a partire dalla componente divergente del vento meridionale e integrata ponderata per pressione nell’intervallo stratosferico tra 34 e 10 hPa, corrispondente approssimativamente alle quote altimetriche superiori della stratosfera polare. Queste anomalie sono rappresentate per il periodo temporale dai giorni 1 a 19 in quattro scenari sperimentali distinti: (a) una simulazione di controllo priva di forzature esterne, (b) un ensemble con forzatura simmetrica (br-symmetry), (c) un ensemble con forzatura a fase 90 gradi (br-phase-90), e (d) un ensemble con forzatura a fase 270 gradi (br-phase-270). La visualizzazione adotta una proiezione polare centrata sul Polo Nord, con contorni continentali delineati per contestualizzare geograficamente i pattern, e include tre righe temporali che segmentano l’evoluzione: giorni 1-5 (fase iniziale di forzatura), giorni 6-12 (fase centrale di forzatura attiva), e giorni 13-19 (fase post-forzatura). La scala cromatica, calibrata in unità di kg/m/s, spazia da valori negativi estremi di -1500 (tonalità viola scuro, indicanti flusso meridionale equatorward, ovvero verso l’equatore) a valori positivi estremi di +1500 (tonalità marrone scuro, indicanti flusso meridionale poleward, ovvero verso il polo), passando attraverso gradazioni intermedie di blu, azzurro, bianco, giallo e arancione che denotano intensità progressive. Le aree punteggiate (stippling) sovrapposte segnalano anomalie statisticamente significative al livello di confidenza del 95%, valutate mediante il test t di Student, il che conferisce robustezza alle interpretazioni dinamiche, distinguendo segnali reali da variazioni casuali nel modello atmosferico.

Dal punto di vista scientifico, la funzione di flusso di massa qui raffigurata funge da diagnostico chiave per quantificare la circolazione meridionale residua nella stratosfera, riflettendo processi di trasporto di massa e momento angolare associati alla circolazione di Brewer-Dobson (BDC) accelerata durante eventi di indebolimento del vortice polare. Anomalie positive (poleward) indicano un rafforzamento del ramo ascendente e meridionale della BDC, che trasporta aria ricca di ozono e momento angolare dalle latitudini basse verso il polo, contribuendo a un riscaldamento adiabatico locale e a un riequilibrio geostrofico; viceversa, anomalie negative (equatorward) segnalano un deflusso compensativo che può amplificare instabilità barocliniche e propagare segnali verso la troposfera. Nella simulazione di controllo (pannello a), durante i giorni 1-5 si osserva un flusso poleward moderato e diffuso attorno al polo, con massimi localizzati sopra l’Eurasia settentrionale e deboli segnali equatorward sull’America del Nord, riflettendo una variabilità di fondo naturale senza forzature; questa intensità si attenua progressivamente nei giorni 6-12, con un pattern più frammentato e meno significativo statisticamente, per poi evolvere nei giorni 13-19 verso un regime misto con aree equatorward emergenti sopra il Pacifico settentrionale, suggestive di una dissipazione energetica spontanea.

Nell’ensemble con forzatura simmetrica (pannello b), il flusso poleward appare distribuito in modo omogeneo e circumpolare, con un’intensità che culmina nei giorni 6-12 – fase di picco della forzatura – raggiungendo valori prossimi a +1200 kg/m/s al centro polare, e aree punteggiate estese che indicano una risposta simmetrica e statisticamente robusta; questa uniformità riflette l’assenza di asimmetrie zonali nella forzatura, favorendo un rafforzamento bilanciato della BDC senza focalizzazioni regionali, con una graduale attenuazione post-forzatura (giorni 13-19) che mantiene tuttavia una componente poleward residua. Al contrario, nell’ensemble fase-90 (pannello c), il flusso poleward si concentra prevalentemente sull’emisfero orientale, con un massimo intenso (oltre +900 kg/m/s) localizzato sopra la Siberia orientale nei giorni 1-5, dove lo stippling copre ampie porzioni indicative di una risposta dinamica indotta dalla fase specifica della forzatura ondulatoria di tipo 1; nei giorni 6-12, questo pattern si sposta leggermente verso ovest e si indebolisce, mentre nei giorni 13-19 emergono anomalie equatorward significative (fino a -600 kg/m/s) sopra il Nord America, suggerendo un’inversione circolatoria post-forzatura legata a meccanismi di feedback come l’aggiustamento di Eliassen o interazioni con onde planetarie riflesse.

Infine, nell’ensemble fase-270 (pannello d), il flusso poleward è marcatamente asimmetrico e confinato all’emisfero occidentale, con un picco iniziale (giorni 1-5) superiore a +1200 kg/m/s sopra il Nord America e l’Atlantico settentrionale, aree densamente punteggiate che evidenziano un impatto localizzato della forzatura opposta; questa intensità persiste nei giorni 6-12 con uno spostamento verso est, per poi transire nei giorni 13-19 verso un regime dominato da anomalie equatorward estese (oltre -900 kg/m/s) sopra l’Eurasia, riflettendo una propagazione discendente di energia che potrebbe amplificare anomalie troposferiche. Complessivamente, questi pattern temporali e spaziali dimostrano come fasi opposte della forzatura ondulatoria modulino la geometria della circolazione stratosferica, influenzando il trasferimento di energia verso la troposfera e potenzialmente innescando feedback come l’amplificazione di blocking atmosferico o variazioni nel getto polare, con implicazioni per la prevedibilità di eventi estremi invernali su scala emisferica.

La Figura 8 raffigura le anomalie della funzione di flusso di massa, derivata dalla componente divergente del vento meridionale e integrata ponderata per pressione nell’intervallo troposferico tra 700 e 321 hPa, che corrisponde approssimativamente alla troposfera libera, una regione chiave per i processi di accoppiamento stratosferico-troposferico e per la propagazione di segnali dinamici verso quote inferiori. Queste anomalie sono visualizzate per il periodo dai giorni 1 a 19 in quattro scenari sperimentali: (a) una simulazione di controllo senza forzature esterne, (b) un ensemble con forzatura simmetrica (br-symmetry), (c) un ensemble con forzatura a fase 90 gradi (br-phase-90), e (d) un ensemble con forzatura a fase 270 gradi (br-phase-270). La proiezione è polare centrata sul Polo Nord, con delineazioni continentali per un orientamento geografico preciso, e l’evoluzione temporale è segmentata in tre righe: giorni 1-5 (fase iniziale di forzatura), giorni 6-12 (fase centrale di forzatura attiva), e giorni 13-19 (fase post-forzatura). La scala cromatica, calibrata in unità di kg/m/s, si estende da valori negativi estremi di -2500 (tonalità viola scuro, indicanti flusso meridionale equatorward, ovvero verso l’equatore, associato a deflusso di aria fredda polare) a valori positivi estremi di +2500 (tonalità marrone scuro, indicanti flusso meridionale poleward, ovvero verso il polo, legato al trasporto di aria più calda da latitudini inferiori), con transizioni intermedie attraverso blu, azzurro, bianco, giallo e arancione che riflettono gradazioni di intensità. Le aree punteggiate (stippling) denotano anomalie statisticamente significative al 95% di confidenza, determinate tramite il test t di Student, fornendo una misura di robustezza contro la variabilità modellistica e confermando la rilevanza dinamica dei pattern osservati.

Dal punto di vista scientifico, questa funzione di flusso di massa nella troposfera libera funge da indicatore diagnostico per la circolazione meridionale residua, catturando meccanismi di compensazione alla circolazione di Brewer-Dobson (BDC) accelerata nella stratosfera (come osservato nella Figura 7), dove flussi poleward dominano durante la forzatura. In troposfera, anomalie negative (equatorward) riflettono un deflusso compensativo di massa e momento angolare dalle alte latitudini, che può amplificare l’intrusione di aria artica fredda verso i subtropici, favorendo instabilità barocliniche e anomalie nel getto polare; al contrario, anomalie positive (poleward) indicano un afflusso di aria subtropicale, contribuendo a riscaldamento adiabatico e a configurazioni di blocking atmosferico. Questo pattern opposto rispetto alla stratosfera evidenzia un accoppiamento dinamico verticale, guidato da principi come la conservazione della vorticità potenziale e l’aggiustamento geostrofico, che trasferiscono perturbazioni stratosferiche (ad esempio, da SSW) verso la troposfera, influenzando la variabilità sinottica e le anomalie termiche superficiali. Nella simulazione di controllo (pannello a), nei giorni 1-5 si osserva un flusso equatorward moderato e frammentato attorno al polo, con valori prossimi a -1000 kg/m/s sopra il Nord America e segnali poleward deboli sopra l’Eurasia subtropicale, riflettendo una variabilità di fondo stagionale; questa intensità si amplifica nei giorni 6-12 con pattern più diffusi e punteggiati, per poi stabilizzarsi nei giorni 13-19 in un regime misto, con aree equatorward persistenti sopra il Pacifico settentrionale, suggestive di dinamiche naturali non forzate.

Nell’ensemble con forzatura simmetrica (pannello b), il flusso equatorward domina in modo circumpolare attorno al polo fin dai giorni 1-5, con intensità che culmina nei giorni 6-12 (valori fino a -2000 kg/m/s al centro polare, con stippling esteso), indicando una risposta compensativa simmetrica alla forzatura stratosferica; nei giorni 13-19, questa componente si attenua, mentre emergono flussi poleward significativi sopra le regioni subtropicali eurasiatiche, riflettendo un riequilibrio post-forzatura che potrebbe perpetuare anomalie troposferiche attraverso feedback eddy-medi. Al contrario, nell’ensemble fase-90 (pannello c), il flusso equatorward si concentra sull’emisfero orientale nei giorni 1-5, con un massimo intenso (oltre -1500 kg/m/s) localizzato sopra la Siberia e l’Eurasia orientale, aree densamente punteggiate che evidenziano un impatto asimmetrico indotto dalla fase della forzatura ondulatoria; nei giorni 6-12, questo pattern si indebolisce e si frammenta, mentre nei giorni 13-19 appaiono anomalie poleward marcate (fino a +1000 kg/m/s) sopra l’Europa occidentale e meridionale, suggerendo un’inversione circolatoria che potrebbe favorire warming regionale attraverso advective transport. Infine, nell’ensemble fase-270 (pannello d), il flusso equatorward è asimmetrico e prevalente sull’emisfero occidentale nei giorni 1-5, con picchi superiori a -2000 kg/m/s sopra l’Alaska e il Nord America, supportati da stippling diffuso; questa intensità persiste e si sposta leggermente nei giorni 6-12, per poi evolvere nei giorni 13-19 verso flussi poleward estesi sopra l’Eurasia settentrionale (oltre +1500 kg/m/s), indicando una propagazione di energia discendente che amplifica asimmetrie zonali.

Complessivamente, questi pattern spaziali e temporali, che mostrano una relazione inversa con quelli stratosferici della Figura 7 (dove poleward prevale durante la forzatura), sottolineano il ruolo della troposfera libera come ponte dinamico per l’accoppiamento verticale, con flussi equatorward che facilitano il trasporto di aria fredda polare e flussi poleward che modula warming subtropicale. Tale dinamica, influenzata dalla fase della forzatura ondulatoria, ha implicazioni per la prevedibilità di eventi estremi invernali, come ondate di freddo regionali, attraverso meccanismi come l’alterazione del gradiente meridionale di temperatura e l’interazione con onde sinottiche, fornendo basi per ulteriori modellizzazioni dell’impatto stratosferico sul clima troposferico.

La Figura 9 raffigura le anomalie della funzione di flusso di massa, derivata dalla componente divergente del vento meridionale e integrata ponderata per pressione nell’intervallo della bassa troposfera tra 970 e 850 hPa, che corrisponde approssimativamente alle quote vicino alla superficie dove i processi dinamici influenzano direttamente le condizioni meteorologiche e climatiche al suolo, inclusi trasporti advettivi di calore e umidità. Queste anomalie sono visualizzate per il periodo dai giorni 1 a 19 in quattro scenari sperimentali: (a) una simulazione di controllo priva di forzature esterne, (b) un ensemble con forzatura simmetrica (br-symmetry), (c) un ensemble con forzatura a fase 90 gradi (br-phase-90), e (d) un ensemble con forzatura a fase 270 gradi (br-phase-270). La proiezione è polare centrata sul Polo Nord, con delineazioni continentali per un contestualizzazione geografica accurata, e l’evoluzione temporale è divisa in tre righe: giorni 1-5 (fase iniziale di forzatura), giorni 6-12 (fase centrale di forzatura attiva), e giorni 13-19 (fase post-forzatura). La scala cromatica, calibrata in unità di kg/m/s, si estende da valori negativi estremi di -2000 (tonalità viola scuro, indicanti flusso meridionale equatorward, ovvero verso l’equatore, associato al deflusso di aria fredda polare verso latitudini inferiori) a valori positivi estremi di +2000 (tonalità marrone scuro, indicanti flusso meridionale poleward, ovvero verso il polo, legato all’afflusso di aria più calda da regioni subtropicali), con gradazioni intermedie attraverso blu, azzurro, bianco, giallo e arancione che denotano intensità variabili. Le aree punteggiate (stippling) sovrapposte segnalano anomalie statisticamente significative al livello di confidenza del 95%, valutate mediante il test t di Student, fornendo una validazione robusta contro la variabilità interna del modello e confermando l’importanza dinamica dei pattern osservati.

Dal punto di vista scientifico, questa funzione di flusso di massa nella bassa troposfera funge da diagnostico essenziale per quantificare la circolazione meridionale residua vicino al suolo, catturando meccanismi di accoppiamento verso il basso che propagano segnali dalla stratosfera (Figura 7) attraverso la troposfera libera (Figura 8) fino alle quote superficiali, influenzando anomalie termiche regionali attraverso processi advettivi e subsidenza adiabatica. Anomalie negative (equatorward) riflettono un trasporto meridionale di masse d’aria artica fredda verso i subtropici, che può amplificare raffreddamenti superficiali favorendo intrusioni di aria polare e instabilità sinottiche; viceversa, anomalie positive (poleward) indicano un afflusso di aria subtropicale più calda, contribuendo a riscaldamenti locali e potenzialmente rafforzando configurazioni anticicloniche persistenti. Questo livello di pressione, più basso rispetto alle figure precedenti, evidenzia il culmine della propagazione discendente, con pattern che raggiungono il picco post-forzatura (giorni 13-19), in linea con l’evoluzione delle anomalie termiche superficiali, suggerendo feedback come l’interazione con la superficie (ad esempio, amplificazione artica) o con onde sinottiche. Nella simulazione di controllo (pannello a), nei giorni 1-5 si osserva un flusso equatorward moderato e diffuso attorno al polo, con valori prossimi a -800 kg/m/s sopra il Nord America e segnali poleward deboli sopra l’Eurasia subtropicale, riflettendo una variabilità naturale di fondo; questa intensità si amplifica nei giorni 6-12 con pattern frammentati e punteggiati, per poi stabilizzarsi nei giorni 13-19 in un regime misto, con aree equatorward persistenti sopra il Pacifico settentrionale, indicative di dinamiche non forzate.

Nell’ensemble con forzatura simmetrica (pannello b), il flusso equatorward domina in modo omogeneo attorno al polo fin dai giorni 1-5, con un’intensità che culmina nei giorni 6-12 (valori fino a -1600 kg/m/s al centro polare, con stippling esteso), rappresentando una compensazione simmetrica alla forzatura stratosferica; nei giorni 13-19, questa componente si attenua gradualmente, mentre emergono flussi poleward significativi sopra le regioni subtropicali eurasiatiche, riflettendo un riequilibrio post-forzatura che potrebbe perpetuare anomalie termiche attraverso meccanismi di feedback baroclinici. Al contrario, nell’ensemble fase-90 (pannello c), il flusso equatorward si concentra sull’emisfero orientale nei giorni 1-5, con un massimo intenso (oltre -1200 kg/m/s) localizzato sopra la Siberia e l’Eurasia orientale, aree densamente punteggiate che evidenziano un impatto asimmetrico indotto dalla fase della forzatura ondulatoria di tipo 1; nei giorni 6-12, questo pattern si indebolisce e si sposta leggermente, mentre nei giorni 13-19 appaiono anomalie poleward marcate (fino a +800 kg/m/s) sopra l’Europa meridionale e occidentale, suggerendo un’inversione circolatoria che favorisce warming regionale attraverso advective transport di aria calda. Infine, nell’ensemble fase-270 (pannello d), il flusso equatorward è asimmetrico e prevalente sull’emisfero occidentale nei giorni 1-5, con picchi superiori a -1600 kg/m/s sopra l’Alaska e il Nord America, supportati da stippling diffuso; questa intensità persiste e si intensifica nei giorni 6-12, per poi evolvere nei giorni 13-19 verso flussi poleward estesi sopra l’Eurasia settentrionale (oltre +1200 kg/m/s), indicando una propagazione di energia discendente che amplifica asimmetrie zonali e contribuisce a raffreddamenti o riscaldamenti differenziati.

Complessivamente, questi pattern spaziali e temporali, che mostrano una continuità inversa con quelli delle figure superiori (poleward in stratosfera, equatorward in troposfera), sottolineano il ruolo della bassa troposfera come interfaccia finale per l’accoppiamento verticale, con flussi equatorward che facilitano il trasporto di aria fredda polare verso i subtropici (favorendo cooling su Alaska ed Eurasia orientale in fase-90, o su Eurasia settentrionale in fase-270) e flussi poleward che modula warming subtropicale (più pronunciato in Eurasia in fase-270). Tale dinamica, modulata dalla fase della forzatura ondulatoria, ha implicazioni significative per la prevedibilità di eventi estremi invernali, come ondate di freddo o periodi miti, attraverso meccanismi come l’alterazione del gradiente meridionale di temperatura, l’interazione con il getto polare e potenziali feedback superficiali (ad esempio, con la copertura nevosa o oceanica), fornendo basi per ulteriori indagini modellistiche sull’impatto stratosferico sul clima superficiale.

Discussione e Analisi: Approfondimenti sui Meccanismi Dinamici e Sfide Osservative

Una comprensione più approfondita della risposta superficiale alle asimmetrie zonali del vortice polare stratosferico, con particolare enfasi sul fenomeno della riflessione delle onde atmosferiche durante eventi di indebolimento del vortice (noti come stati di vortice debole), riveste un’importanza cruciale per la modellistica climatica e la previsione meteorologica stagionale. Studi precedenti, basati su dataset osservativi derivati da reanalisi atmosferiche come ERA5 o MERRA-2, hanno evidenziato evidenze empiriche che la struttura zonale dell’anomalia del vortice – inclusa la sua asimmetria longitudinale – e la propagazione discendente delle onde planetarie giocano un ruolo pivotale negli impatti superficiali, influenzando pattern di temperatura, pressione e precipitazioni su scala emisferica (ad esempio, attraverso modulazioni del North Atlantic Oscillation o dell’Arctic Oscillation). Tuttavia, isolare la causalità in tali contesti osservativi si rivela complesso, poiché lo stato anomalo del vortice è tipicamente preceduto da un’improvvisa intensificazione dell’attività ondulatoria nella troposfera, spesso legata a flussi di calore e momento angolare verticali (come quantificati dal flusso di Eliassen-Palm). Di conseguenza, distinguere gli effetti a lungo termine di questa attività ondulatoria iniziale – che può includere onde Rossby generate da anomalie orografiche o termiche – dalla risposta successiva della stratosfera e della troposfera richiede approcci metodologici avanzati, ma rimane sfidante a causa della non linearità delle interazioni atmosferiche e della limitata campionatura osservativa.

Ad esempio, sia gli eventi di riscaldamento stratosferico improvviso (SSW) di tipo displacement (caratterizzati da uno spostamento del vortice polare dal polo) che quelli di tipo split (con una divisione del vortice in due lobi) influenzano le temperature superficiali, potenzialmente differenziandosi nei pattern spaziali – come estensione latitudinale o persistenza temporale – e nella magnitudo delle anomalie termiche, con impatti che possono estendersi fino a 60-90 giorni post-evento (Hall et al., 2021; Mitchell et al., 2013). Queste diverse morfologie di SSW sono precedute da flussi ondulatori troposferici qualitativamente distinti: gli SSW di tipo displacement sono prevalentemente guidati da una forzatura dominante di onda-1 (tipica di perturbazioni asimmetriche su larga scala, come quelle originate dalle Montagne Rocciose o dall’Himalaya), mentre gli SSW di tipo split coinvolgono una combinazione di forzature di onda-1 e onda-2, che amplificano instabilità barocliniche e risonanze planetarie (Bancala et al., 2012; Cohen e Jones, 2011; White et al., 2019). Tali differenze iniziali nella troposfera possono alterare la propagazione verticale delle onde, influenzando la dissipazione energetica nella stratosfera critica e la conseguente ridistribuzione di vorticità potenziale, con ripercussioni sulla stabilità del getto polare.

Oltre agli eventi SSW, gli eventi di riflessione delle onde – in cui le onde planetarie ascendenti vengono riflesse verso il basso dalla stratosfera, spesso durante fasi di vortice indebolito – possono esercitare un impatto significativo sul meteo superficiale, frequentemente associati a raffreddamenti estremi in regioni come il Nord America, attraverso meccanismi come l’amplificazione di blocking atmosferico o l’intrusione di aria artica (Cohen et al., 2021, 2022; Matthias e Kretschmer, 2020). Tuttavia, anche questi eventi sono generalmente preceduti da un’improvvisa intensificazione dell’attività ondulatoria sopra la Siberia orientale, regione nota per la generazione di onde Rossby stazionarie a causa di gradienti termici marcati tra l’oceano Pacifico e il continente eurasiatico. Isolando l’importanza della riflessione stratosferica – che coinvolge processi come l’interferenza costruttiva tra onde ascendenti e discendenti, quantificabile mediante il flusso verticale di Plumb – rispetto all’attività ondulatoria troposferica iniziale sulla Siberia per la successiva risposta nordamericana, emerge una sfida analitica, data la covarianza tra questi fenomeni e la potenziale retroazione bidirezionale. Sebbene tecniche avanzate, come le reti di effetti causali (basate su grafi diretti aciclici per inferire dipendenze temporali) o approcci di causalità alla Granger/Pearl (che valutano la prevedibilità incrementale tra serie temporali multivariate), siano state applicate con successo per disentangolare tali interazioni (Kretschmer et al., 2018a, b), questi metodi statistici, pur robusti, non forniscono una comprensione meccanica dettagliata dei percorsi fisici sottostanti, limitandosi a correlazioni probabilistiche senza risolvere le equazioni governanti (ad esempio, le equazioni primitive atmosferiche). Inoltre, i meccanismi dinamici specifici che spiegano le risposte superficiali variabili a diversi tipi di SSW e agli eventi di riflessione delle onde associati – inclusi feedback non lineari come l’amplificazione risonante o l’interazione con la variabilità oceanica (ad esempio, El Niño-Southern Oscillation) – rimangono insufficientemente esplorati, richiedendo integrazioni tra osservazioni, modellistica numerica e analisi teoriche.

Per affrontare il ruolo delle asimmetrie zonali del vortice polare e del fenomeno della riflessione delle onde durante gli eventi SSW, si estende l’approccio di White et al. (2022) incorporando momenti toroidali (torques) con una forte componente asimmetrica zonale nella stratosfera del modello MiMA (Modello Idealizzato per la Meccanica Atmosferica), un framework semplificato che risolve le equazioni primitive su una griglia sferica, permettendo simulazioni efficienti di dinamiche stratosferiche-troposferiche (Garfinkel et al., 2020a; Jucker e Gerber, 2017). Utilizzando ogni mezzanotte del 1° gennaio dalla simulazione di controllo (CTRL) come condizione iniziale – garantendo un ensemble di stati atmosferici rappresentativi della variabilità invernale boreale – si generano ramificazioni (branching) e si impone un momento toroidale con diverse fasi longitudinali, simulando perturbazioni asimmetriche simili a forzature ondulatorie naturali. Nel presente studio, l’attenzione si concentra su due esperimenti rappresentativi di perturbazione: l’ensemble fase-90 (con parametri k=1 per onda-1, ΘA=4 per ampiezza termica, λ0=90° per fase longitudinale) e l’ensemble fase-270 (k=1, ΘA=4, λ0=270°), che catturano configurazioni opposte di asimmetria zonale. Ulteriori esperimenti, inclusi variazioni in ampiezza o wavenumber, saranno esaminati in un articolo successivo. Questi ensemble ramificati vengono confrontati con 48 eventi SSW spontanei simulati nella CTRL, fornendo un benchmark per valutare l’impatto delle forzature imposte rispetto alla variabilità interna del modello.

Il pattern zonale opposto della forzatura di onda-1 genera influenze contrastanti sia nella stratosfera che nella troposfera, evidenziando la sensibilità del sistema atmosferico alle fasi iniziali delle perturbazioni. Durante la fase di forzatura attiva, il vento zonale medio a 60°N e 10 hPa transita bruscamente da regimi occidentali (westerlies) a orientali (easterlies) sotto l’influenza del momento toroidale orientale, alterando il gradiente meridionale di temperatura e inducendo instabilità nel flusso zonale medio (Figura 3a). Contemporaneamente, la calotta polare esperisce un riscaldamento marcato e un innalzamento dell’altezza geopotenziale (Figura 3b,c,d), attribuibile all’intensificazione della subsidenza nella circolazione di Brewer-Dobson (BDC), che promuove un riscaldamento adiabatico attraverso la compressione di masse d’aria discendenti e il trasporto di ozono stratosferico (Baldwin et al., 2021; Garfinkel et al., 2010). Negli esperimenti fase-90 e fase-270, le anomalie del vento zonale stratosferico raggiungono il picco in emisferi opposti – orientale per fase-90 e occidentale per fase-270 – riflettendo la dipendenza dalla fase longitudinale della forzatura, e questa differenziazione stratosferica si manifesta in un impatto qualitativamente distinto nella troposfera, con anomalie che propagano verso il basso attraverso meccanismi come l’aggiustamento di Eliassen (Figura 1). Inoltre, il vortice polare si disloca dal polo geografico, ma in direzioni longitudinali opposte a 10 hPa durante i giorni 1-5 (Figura 2c,d), a causa della fase antitetica della componente di onda-1 della forzatura, che modula l’interferenza tra il flusso zonale e le perturbazioni ondulatorie. Le anomalie dell’altezza geopotenziale nella troposfera mostrano spostamenti analoghi dal polo, con direzioni che rispecchiano quelle stratosferiche e che variano in modo opposto tra gli esperimenti, indicando un accoppiamento verticale coerente mediato da processi come la ridistribuzione di vorticità potenziale e il flusso di energia discendente. Complessivamente, la fase opposta della forzatura di onda-1 induce impatti marcatamente diversi sui campi di vento zonale, temperatura e altezza geopotenziale attraverso l’intero profilo stratosferico-troposferico, per mezzo di un accoppiamento dinamico verso il basso che coinvolge feedback non lineari e propagazione di anomalie su scala sinottica.Inoltre, questo accoppiamento dinamico verso il basso esercita un’influenza significativa sui pattern di temperatura a 2 metri di quota – un livello standard per le misurazioni meteorologiche superficiali che cattura le condizioni termiche percepite al suolo – e sui regimi di precipitazioni, modulando processi come il trasporto advettivo di umidità e calore su scala sinottica. Segnali evidenti di raffreddamento localizzato, riscaldamento anomalo e precipitazioni intensificate si manifestano già durante la fase attiva della forzatura stratosferica, sebbene le loro distribuzioni spaziali varino in modo marcato tra i diversi ensemble sperimentali, riflettendo la sensibilità del sistema atmosferico alle fasi iniziali delle perturbazioni ondulatorie. In particolare, i segnali di raffreddamento e riscaldamento si amplificano in modo sostanziale nelle fasi successive alla disattivazione della forzatura stratosferica (come illustrato nella Figura 4), un ritardo temporale che suggerisce l’attivazione di feedback interni all’atmosfera, quali interazioni non lineari tra onde planetarie e il flusso zonale medio, o amplificazioni risonanti che perpetuano le anomalie dinamiche. Al contrario, le anomalie di precipitazioni – spesso legate a variazioni nella convergenza di umidità e nella stabilità convettiva – raggiungono il loro picco immediatamente dopo la cessazione della forzatura (nei giorni 13-19), per poi attenuarsi progressivamente (Figura 5), indicando una risposta più rapida del ciclo idrologico troposferico alle perturbazioni discendenti, potenzialmente mediata da alterazioni nel getto subtropicale o in pattern di blocking atmosferico.

Nell’ensemble fase-90, il raffreddamento superficiale risulta più pronunciato sulle regioni dell’Alaska e dell’Eurasia orientale (Figura 4c), aree vulnerabili a intrusioni di aria artica fredda dovute a deformazioni del getto polare, con anomalie che possono raggiungere diversi gradi Kelvin al di sotto della media climatologica, influenzando ecosistemi e attività umane. Parallelamente, le precipitazioni intensificate si concentrano sul Nord Pacifico e sull’Atlantico settentrionale, estendendosi fino al nord-ovest dell’Europa (Figura 5c), un pattern che riflette un rafforzamento della convergenza umida associata a cicloni extratropicali amplificati da gradienti termici meridionali accentuati. Tuttavia, nell’ensemble fase-270, si osservano raffreddamenti significativi sull’America del Nord centrale e sull’Eurasia settentrionale (Figura 4d), con impatti che coinvolgono vasti settori continentali e possono alterare il bilancio energetico superficiale attraverso meccanismi come l’amplificazione artica o feedback con la copertura nevosa. In questo scenario, le precipitazioni intensificate sul Nord Pacifico e sull’Atlantico settentrionale si estendono fino all’Eurasia subtropicale (Figura 5d), suggerendo una ridistribuzione spaziale dell’umidità legata a spostamenti zonali nelle configurazioni barocliniche. Queste anomalie di temperatura superficiale negli ensemble fase-90 e fase-270 si allineano con le anomalie termiche osservate nei cluster 4 e 5 identificati da Kretschmer et al. (2018a), dove il cluster 4 corrisponde a regimi associati a riflessioni ondulatorie sopra la Siberia orientale che favoriscono ondate di freddo in Nord America, mentre il cluster 5 riflette pattern legati a indebolimenti del vortice polare che amplificano raffreddamenti sull’Eurasia settentrionale, basati su analisi di clustering della geopotenziale a 100 hPa per catturare stati persistenti del vortice debole senza presupporre eventi SSW predefiniti.

Inoltre, la propagazione discendente delle onde – quantificata mediante il flusso verticale di Plumb, un diagnostico che misura il trasferimento di energia e momento angolare ondulatorio – si verifica prevalentemente sopra il Nord America nell’ensemble fase-90 (Figura 6c), in accordo con le configurazioni osservate nel cluster 4 di Kretschmer et al. (2018a), che evidenziano un meccanismo riflettivo dove le onde planetarie ascendenti vengono deviate verso il basso, alterando la stabilità del getto e favorendo intrusioni di aria fredda. Al contrario, tale propagazione è concentrata sopra l’Eurasia settentrionale nell’ensemble fase-270 (Figura 6d), coerente con gli eventi di riflessione delle onde nel cluster 5 descritti da Kretschmer et al. (2018a), che coinvolgono interazioni tra onde di tipo 1 e il flusso zonale stratosferico, con ripercussioni su scala regionale. Questi risultati indicano che la fase opposta della forzatura di onda-1 – una perturbazione planetaria dominante nelle dinamiche invernali boreali, caratterizzata da un wavenumber zonale di 1 e ampiezze che possono superare i 4 km in geopotenziale – modula in modo differenziale gli eventi di SSW e la propagazione discendente delle onde, attraverso processi come l’interferenza costruttiva o distruttiva con il vortice polare, portando a impatti distinti sui pattern di temperatura superficiale e precipitazioni, con implicazioni per la variabilità climatica e la prevedibilità sub-stagionale.

Si chiarisce il percorso dell’accoppiamento stratosfera-troposfera associato agli SSW e agli eventi di riflessione delle onde – un meccanismo che coinvolge il trasferimento di anomalie di vorticità potenziale e il riequilibrio geostrofico su scale verticali – analizzando la funzione di flusso di massa derivata dalla componente divergente del vento meridionale, un diagnostico che quantifica la circolazione meridionale residua e il trasporto di massa nella circolazione di Brewer-Dobson (BDC) accelerata. In tutti gli esperimenti simulati, si osserva un significativo movimento verso i poli nella stratosfera (Figura 7), che riflette un rafforzamento del ramo poleward della BDC, trasportando aria ricca di momento angolare dalle latitudini basse e contribuendo a un riscaldamento adiabatico locale attraverso subsidenza compressiva. Parallelamente, emerge un pronunciato movimento verso l’equatore nella troposfera (Figura 8), un flusso compensativo che altera il gradiente meridionale di temperatura e amplifica instabilità barocliniche, favorendo la propagazione di anomalie verso quote inferiori. Tuttavia, la funzione di flusso di massa meridionale esibisce pattern di dipole zonale opposti tra gli esperimenti fase-90 e fase-270, evidenziando l’impatto delle asimmetrie iniziali sulla geometria circolatoria. Nella stratosfera, durante l’applicazione del momento toroidale – un torque che simula forzature esterne come quelle da onde planetarie – il movimento verso i poli è prevalentemente localizzato nell’emisfero orientale nell’ensemble fase-90 (Figura 7c), con massimi che possono superare i 1500 kg/m/s in termini di flusso di massa, mentre è centrato sull’emisfero occidentale nell’ensemble fase-270 (Figura 7d), riflettendo una dipendenza dalla fase longitudinale della forzatura. Questa asimmetria zonale nella circolazione meridionale stratosferica si riflette nella troposfera libera, indicando un accoppiamento verso il basso robusto tra stratosfera e troposfera, mediato da principi come la conservazione del momento angolare e l’aggiustamento di Eliassen, che facilitano il trasferimento di segnali dinamici attraverso il tropopausa. Nella bassa troposfera (Figura 9), il flusso verso l’equatore trasporta aria più fredda da latitudini superiori a quelle inferiori, contribuendo al raffreddamento superficiale su Alaska ed Eurasia orientale nell’ensemble fase-90 (Figura 4c) – con anomalie che possono raggiungere -5 K o più, alterando il bilancio radiativo superficiale – e su Eurasia settentrionale nell’ensemble fase-270 (Figura 4d). Al contrario, il flusso verso i poli convoglia aria più calda dalle regioni subtropicali verso i poli, risultando in un riscaldamento superficiale sull’Europa nell’ensemble fase-90 (Figura 4c), potenzialmente amplificato da feedback con la superficie oceanica, e sull’Eurasia subtropicale nell’ensemble fase-270 (Figura 4d), con implicazioni per pattern come il Siberian High o l’East Asian Winter Monsoon.

I risultati di questo studio rafforzano le evidenze emerse da studi precedenti basati su rianalisi atmosferiche come ERA-Interim o MERRA-2, che collegano i pattern di raffreddamento e riscaldamento superficiale alla riflessione delle onde stratosferiche (Cohen et al., 2021, 2022; Kretschmer et al., 2018a), dove Cohen et al. (2021) hanno evidenziato come eventi di stretching del vortice polare stratosferico – un’estensione zonale asimmetrica associata a riflessioni ondulatorie – siano in aumento nell’era satellitare e legati a ondate di freddo estreme in Nord America e Asia, inclusa la tempesta invernale del febbraio 2021 in Texas, con temperature che scesero fino a -20°C in alcune aree. In particolare, per la prima volta, per quanto ne sappiamo, si è dimostrato attraverso esperimenti modellistici controllati che una forzatura stratosferica imposta, capace di generare attività ondulatoria propagante verso il basso con flussi negativi di Plumb indicativi di riflessione, può causare direttamente gli impatti superficiali osservati in studi empirici precedenti, distinguendo tra meccanismi assorbenti (tipici di SSW canonici con inversione dei venti zonali) e riflettenti (che mantengono un vortice forte ma deviano energia verso la troposfera). Il lavoro suggerisce inoltre che la direzione dello spostamento del vortice polare durante un evento SSW di tipo displacement – spostamento dal polo verso longitudini specifiche, con deviazioni che possono superare i 1000 km – riveste un’importanza cruciale per la risposta superficiale, influenzando la localizzazione di anomalie come il rafforzamento di ridge anticicloniche o trough cicloniche. Ricerche in corso esploreranno altre fasi longitudinali per una forzatura di onda-1, variando parametri come l’ampiezza o il wavenumber, e anche forzature di onda-2 – che introducono pattern biforcati del vortice, noti come split SSW – oltre a quantificare in modo più preciso l’importanza relativa del flusso meridionale divergente (legato alla BDC e al trasporto di massa) rispetto al flusso rotazionale anomalo (associato a vortici sinottici e instabilità zonali) per gli impatti sulla temperatura superficiale, potenzialmente attraverso analisi di bilancio energetico o decomposizioni spettrali.

https://doi.org/10.5194/egusphere-2025-4334