A Mechanism for the Poleward Propagation of Zonal Mean Flow Anomalies- Un meccanismo per la propagazione verso i poli delle anomalie del flusso zonale medio

SUKYOUNG LEE , SEOK -WOO SON , AND KEVIN GRISE
Department of Meteorology, The Pennsylvania State University, University Park, Pennsylvania
STEVEN B. FELDSTEIN
Earth and Environmental Systems Institute, The Pennsylvania State University, University Park, Pennsylvania
(Manuscript received 30 November 2005, in final form 30 June 2006)

Nella variabilità intrastagionale della circolazione extratropicale emerge spesso un comportamento “a migrazione”: le anomalie del vento zonale medio (cioè la componente del vento da Ovest verso Est mediata lungo i paralleli) tendono a spostarsi progressivamente verso latitudini più alte con una certa quasi-periodicità. Questo segnale, già documentato in analisi osservazionali e reanalisi (ad esempio nella propagazione polare delle anomalie del momento angolare relativo zonal-mean) e collegato concettualmente alla famiglia dei pattern annulari/zonal index, è interessante perché suggerisce che la circolazione non fluttui solo “in ampiezza”, ma anche “in posizione” lungo il meridiano. 

Il lavoro di Lee, Son, Grise e Feldstein (pubblicato su Journal of the Atmospheric Sciences) affronta proprio la domanda fisica cruciale: che cosa “muove” l’anomalia del flusso zonale medio verso i poli? L’idea centrale è che la propagazione non sia un semplice effetto lineare, ma il prodotto di una sequenza autosostenuta in cui interagiscono (i) propagazione meridionale delle onde di Rossby, (ii) breaking non lineare con mescolamento/omogeneizzazione della vorticità potenziale (PV), e (iii) rilassamento radiativo che ristabilizza la fascia tropicale. In questo quadro, l’innesco è segnato dalla comparsa di anomalie negative del vento zonale nei Tropici: esse non nascono “localmente” per forzante tropicale, ma come conseguenza del fatto che onde generate alle medie latitudini riescono a propagare meridionalmente e poi “rompono” (overturning/breaking) a latitudini più basse, riorganizzando PV e vento medio. 

Il passaggio dinamico chiave è il ruolo del breaking come meccanismo irreversibile di rimescolamento della PV: quando le onde rompono, stirano e ripiegano i contorni di PV fino a produrre filamenti e regioni a PV quasi omogenea (PV mixing), e questa ristrutturazione modifica direttamente il gradiente meridionale di PV e quindi la “guida d’onda” (waveguide) per le onde successive. In termini qualitativi, si forma una “fascia” in cui la PV è più omogenea associata alle anomalie di vento (negative ai Tropici e positive in fascia subtropicale nel loro scenario): tale fascia diventa una barriera dinamica che riduce o interrompe la successiva propagazione equatoriale delle onde di media latitudine. Il risultato è che l’attività ondosa, non potendo più scaricare momento/vorticità verso l’equatore, tende a concentrarsi e a rompere più a nord, sul margine polare della regione a PV omogeneizzata. È proprio questo spostamento del luogo di assorbimento/breaking e dei flussi di vorticità degli eddies che genera “nuove” coppie di anomalie (westerly/easterly) a latitudini leggermente più elevate, permettendo alla firma zonal-mean di avanzare verso il polo a passi successivi. 

Un elemento spesso sottovalutato, ma determinante nella loro interpretazione, è il rilassamento radiativo: una volta che la regione tropicale è “schermata” dall’invasione di onde extratropicali, il modello tende a riportare i venti tropicali verso lo stato indisturbato (westerlies “di base”) tramite la termodinamica semplificata di Newtonian cooling/relaxation. In parallelo, alle latitudini extratropicali le anomalie di westerlies vengono invece sostenute dalla convergenza dei flussi di quantità di moto/vorticità associati agli eddies (eddy–mean flow interaction). Questa divisione dei ruoli è importante perché suggerisce che la propagazione sia, in sostanza, una “staffetta”: il breaking ridisegna PV e waveguide, i flussi degli eddies costruiscono l’anomalia sul bordo polare, e il rilassamento radiativo ripulisce i Tropici, preparando il campo per il passo successivo. 

Questa cornice si innesta naturalmente nella letteratura sugli annular modes e sul feedback eddy–jet. I pattern annulari (NAM/SAM) possono essere visti come la principale modalità di variabilità del vento zonale medio extratropicale; la loro persistenza e la loro evoluzione dipendono in modo sostanziale dai flussi di quantità di moto degli eddies e dal modo in cui il jet si sposta/si riorganizza (feedback positivo o negativo tra anomalie del getto e storm track). Studi classici e successivi hanno quantificato come le anomalie del vento associate a queste modalità siano in larga parte “guidate” dagli eddies e come l’interazione eddy–mean flow possa introdurre una scala temporale intrastagionale distinta, potenziando la persistenza e, in certe condizioni, favorendo componenti propaganti (non puramente stazionarie) nelle principali EOF del vento zonale medio. 

Non a caso, lavori più recenti hanno riformulato la propagazione polare anche in termini di “annular modes propaganti”, usando strumenti come EOF coniugate e modelli ridotti eddy–zonal flow feedback, mostrando che la propagazione non è necessariamente ubiquitaria: emerge o si indebolisce a seconda del regime dinamico di fondo (forza del waveguide/gradiente di PV, struttura del jet, dissipazione effettiva, ecc.). Questo è coerente con un punto molto concreto del paper: la propagazione compare soprattutto quando il gradiente meridionale di PV ha intensità moderata. Se è troppo forte, la guida d’onda risulta molto efficiente e l’evoluzione tende a rimanere più lineare e meno “traslante”; se è troppo debole, il breaking e il rimescolamento possono diventare più diffusi e meno capaci di produrre una sequenza ordinata di schermatura–riattivazione. In altre parole, serve un equilibrio fine tra propagazione lineare e non linearità del breaking perché la migrazione diventi un ciclo coerente. 

Dal punto di vista operativo (prevedibilità), l’implicazione più interessante è che la propagazione non è solo “rumore”: se il sistema entra in un regime in cui si è già formata una regione a PV omogeneizzata e la propagazione equatoriale delle onde è parzialmente bloccata, allora la posizione del successivo breaking e quindi la tendenza del vento zonale medio hanno una memoria intrastagionale legata alla geometria del waveguide e alla distribuzione dell’attività d’onda. In alcuni contesti osservativi regionali (ad esempio sul Nord Atlantico) sono state cercate tracce dirette del meccanismo proposto, proprio per capire quando la catena propagazione → breaking → riorganizzazione della PV → nuova anomalia sia effettivamente diagnosticabile nei dati. Questo filone è utile perché collega un comportamento “zonal-mean” a processi sinottici molto concreti (breaking ciclonico/anticiclonico, spostamenti del getto, variazioni della storm track), cioè esattamente ciò che i modelli provano a prevedere su scala sub-stagionale. 

1. Introduzione

Negli ultimi decenni la dinamica del flusso zonale medio (cioè la componente del vento zonale mediata lungo i paralleli) è tornata al centro dell’interesse perché rappresenta una sorta di “scheletro” della circolazione extratropicale: quando cambia l’intensità o la latitudine dei westerlies medi, cambiano anche storm track, scambi meridiani e regimi di blocco. Storicamente, questo filone nasce con il concetto di zonal index e del suo ciclo: Rossby (1939) e, soprattutto, Namias (1950) cercavano di descrivere la variabilità della circolazione come alternanza tra stati più zonali e stati più ondulati/“bloccati”, con implicazioni dirette per la prevedibilità oltre la scala sinottica. 

Dentro questa tradizione, l’Introduzione del lavoro mette in evidenza un aspetto che, pur comparendo in letteratura da molti anni, è stato relativamente poco “internalizzato” come fenomeno sistematico dalla comunità: la propagazione verso i poli delle anomalie del vento zonale medio, dall’equatore fino a ~70° in entrambi gli emisferi, con una velocità quasi costante. L’idea è forte perché non descrive solo un’anomalia che cresce o decade “in loco”, ma una struttura che migra meridionalmente in modo ordinato, come se la circolazione avesse un grado di organizzazione intrinseco su scala intrastagionale. Gli autori sottolineano che questo comportamento emerge sia in modelli idealizzati (citano James et al. 1994) sia nelle osservazioni (Riehl et al. 1950; Feldstein 1998), e che esempi di questo tipo di migrazione risultano chiaramente diagnosticabili in analisi su reanalisi tramite grandezze come il momento angolare relativo zonal-mean (MR) e i campi di vento zonale mediato. 

Un punto concettuale molto importante, che l’Introduzione chiarisce bene, è in che senso questa propagazione sia “diversa” da altre quasi-periodicità celebri. Fenomeni come la Madden–Julian Oscillation (intrastagionale, tropicale, a propagazione prevalentemente zonale verso Est) o la QBO (stratosferica, con propagazione verticale di shear alternati) hanno firme spazio-temporali relativamente “pulite”, ma non sono tipicamente caratterizzate dal mantenimento coerente della struttura nel cuore della troposfera extratropicale, dove la turbolenza e l’attività d’onda rendono il flusso molto energetico e caotico. La propagazione polare delle anomalie zonal-mean è invece notevole proprio perché conserva una geometria riconoscibile anche alle medie latitudini, suggerendo un legame stretto con i processi di interazione eddy–mean flow (trasporto e convergenza di quantità di moto e vorticità da parte degli eddies). 

Gli autori insistono poi sulle implicazioni previsionali: già nel primo studio citato sul tema, Riehl e colleghi (1950) collegarono una fase specifica di questa evoluzione propagante a un aumento dell’attività ciclonica tropicale, cioè a un possibile “ponte” tra variabilità zonal-mean e risposta dei tropici (o quantomeno una covariazione diagnostica utile). Ma il punto operativo più intrigante è un altro: se un osservatore a latitudine fissa vede passare queste anomalie come una sequenza abbastanza regolare, allora può percepire un’alternanza quasi periodica di westerlies più forti e più deboli; in quel caso la prevedibilità dell’indice zonale può aumentare sensibilmente rispetto a un’atmosfera dove le anomalie nascono e muoiono localmente senza una “traiettoria” coerente. È esattamente il motivo per cui, storicamente, il zonal index sembrava promettere molto per l’extended range, salvo poi deludere quando ci si accorse che la periodicità non era robusta in modo generale; l’Introduzione qui recupera quel dibattito in modo moderno: non “promette” periodicità universale, ma suggerisce che la propagazione persistente possa creare finestre dinamiche con memoria più lunga. 

Sul piano dei meccanismi, l’Introduzione passa in rassegna due idee guida che convergono su un messaggio: il motore non è un forcing esterno semplice, ma un circuito interno di feedback tra onde e flusso medio. Da un lato, James e Dodd (1996) propongono che un’anomalia di vento zonale “confinata” in latitudine modifichi ampiezza e numero d’onda meridionale degli eddies, riorganizzando la convergenza del flusso di quantità di moto degli eddies (eddy momentum flux convergence) in modo tale da spingere l’anomalia verso la regione sorgente delle onde; se la sorgente è a latitudini più alte, la deriva risultante è poleward. Dall’altro lato, Robinson (2000) enfatizza un quadro “baroclinico”: la struttura a bassa troposfera (baroclinità, gradienti termici e quindi generazione di eddies) può retroagire sistematicamente sul getto, alimentando un eddy feedback che organizza la variabilità del zonal index. Entrambe le visioni sono coerenti con l’evidenza osservativa di Feldstein (1998), che mette al centro il ruolo essenziale dei transient eddy momentum fluxesnel guidare la propagazione verso i poli delle anomalie zonal-mean. 

A questo punto l’Introduzione prepara il terreno per la scelta metodologica del paper: usare un modello numerico idealizzato non lineare. La logica è molto solida: le osservazioni sono ricche ma “confondenti” (mille forzanti, mille retroazioni), mentre un modello idealizzato permette di testare ipotesi causali e, soprattutto, di verificare quanto la propagazione sia robusta al cambiare del regime di fondo (ad esempio, variando parametri che controllano intensità del getto, gradiente di PV e dissipazione). Questa scelta si inserisce bene nella letteratura sugli annular modes e sul feedback eddy–getto: lì si è visto che molte proprietà (persistenza, varianza a bassa frequenza, struttura spaziale) dipendono fortemente da come gli eddies rinforzano o smorzano le anomalie zonal-mean. In altre parole, l’Introduzione sta dicendo: se capiamo perché e quando le anomalie si mettono in marcia verso il polo, capiamo anche quando il sistema atmosferico guadagna memoria e quindi prevedibilità su scala estesa—ed è esattamente ciò che le sezioni successive proveranno a dimostrare in modo meccanicistico.

Figura 1 – Evoluzione delle anomalie del momento angolare relativo (MR) e propagazione verso i poli

La Figura 1 è un diagramma latitudine–tempo (un Hovmöller meridiano) che mostra come evolvono, lungo la stagione, le anomalie del momento angolare relativo atmosferico calcolate come media zonale: nel pannello (a) per l’inverno boreale 1991–1992, nel pannello (b) per l’estate australe 1990–1991. Le isolinee misurano l’anomalia di MR (intervallo 4×10²³ kg m² s⁻¹): linee continue positive e tratteggiate negative; lo zero è omesso. In termini fisici, MR è una misura “integrata” di quanta circolazione zonale (vento medio da ovest verso est) è presente attorno all’asse terrestre: quando MR è anormalmente alto, tipicamente c’è un surplus di westerlies medi (o una loro traslazione/rafforzamento) nelle fasce latitudinali dove l’anomalia è massima; quando MR è anormalmente basso, vale l’opposto. Proprio per questo Feldstein usa MR come proxy robusto della variabilità del zonal index/annular modes: invece di guardare un singolo livello o una singola regione, si segue una grandezza dinamicamente coerente con l’equazione del momento zonale medio e con il ruolo dei flussi di quantità di moto degli eddies. 

Il messaggio visivo più importante è l’inclinazione “diagonale” delle strutture: le bande di anomalia non restano ferme a una latitudine fissa, ma migrano sistematicamente verso i poli con una velocità quasi costante, dall’area tropicale/subtropicale fino a circa 60–70° in entrambi gli emisferi. Se le anomalie fossero puramente locali e stazionarie, vedresti massimi/minimi che si accendono e spengono in colonne verticali; qui invece la pendenza delle isolinee indica una vera propagazione meridionale su scala intrastagionale (settimane). Feldstein mostra che questo comportamento è riscontrabile in entrambi gli emisferi e a diverse latitudini; inoltre, dalla diagnostica sul reanalysis emerge che la velocità di propagazione tende a essere più rapida in fascia subtropicale rispetto alle alte latitudini, coerentemente con differenze di guida d’onda e di interazione eddy–flusso medio lungo il meridiano. 

Dinamicamente, la figura suggerisce che la variabilità del flusso zonale medio non sia solo “rafforza/indebolisci il jet”, ma un ciclo organizzato in cui il luogo di massima convergenza dei flussi di quantità di moto (e quindi di accelerazione/decelerazione dei westerlies medi) si sposta progressivamente verso nord/sud. Questa lettura si integra bene con la teoria moderna degli annular modes: una parte rilevante della persistenza e dell’evoluzione del vento zonale medio è spiegata dal feedback tra anomalie del getto e risposta degli eddies (storm track), che tende a rinforzare o rimodellare la struttura zonal-mean nel tempo. In reanalisi, Lorenz e Hartmann hanno mostrato come l’interazione eddy–zonal flow sia capace di sostenere e strutturare la variabilità del vento zonale medio, fornendo “memoria” oltre la scala sinottica; Thompson e Wallace hanno inquadrato la geometria anulare di queste anomalie come modalità dominante della circolazione extratropicale. La Figura 1, in questo senso, può essere letta come una manifestazione propagante (non puramente stazionaria) della stessa famiglia di dinamiche. 

Il passo in più — ed è ciò che motiva il paper di Lee et al. — è capire perché la struttura propaghi invece di restare ancorata: il meccanismo proposto è una sequenza orchestrata di (i) propagazione meridionale delle onde di Rossby, (ii) breaking non lineare con mescolamento/omogeneizzazione della PV che “riorganizza” la waveguide e blocca la radiazione equatoriale delle onde, e (iii) rilassamento radiativo che ristabilisce i westerlies tropicali mentre le anomalie extratropicali vengono rigenerate dai flussi di vorticità/momento degli eddies più a nord (o più a sud). In altre parole, la Figura 1 non è solo un tracciato statistico: è la firma osservativa di un ciclo in cui l’assorbimento/breaking dell’attività ondosa e la risposta del flusso medio si spostano di latitudine in latitudine, producendo la “marcia” verso i poli. Esistono anche lavori successivi che cercano evidenza osservativa specifica del meccanismo di Lee et al. nei dati atmosferici, proprio partendo da strutture tipo quelle in Figura 1. 

Infine, la figura chiarisce bene l’aspetto applicativo: se a una latitudine fissa passano in successione anomalie positive e negative (perché la struttura migra), un osservatore locale percepisce un’alternanza quasi periodica di stati con westerlies più forti e più deboli — esattamente la promessa storica dello zonal index per la previsione estesa. Il punto moderno non è “c’è sempre una periodicità”, ma “quando la propagazione è persistente e coerente, il sistema acquista una componente prevedibile intrastagionale”, perché l’anomalia non evolve solo per crescita/decadimento locale ma segue una traiettoria dinamicamente vincolata. Questo è il ponte concettuale tra Figura 1, la diagnostica eddy–mean flow e l’interesse per la predicibilità su scala sub-stagionale. 

Figura 2 – Segmento di 1600 giorni del vento zonale medio a 250 hPa: regime propagante vs regime “zonal index”

La Figura 2 confronta due comportamenti distinti della variabilità del vento zonale medio a 250 hPa (cioè nella fascia tipica del getto troposferico): nel pannello (a) un caso in cui le anomalie del vento zonale medio propagano verso i poli, nel pannello (b) un caso in cui la variabilità è invece dominata dal cosiddetto zonal index, ossia un’alternanza tra stati con westerlies più forti e più deboli che rimane sostanzialmente ancorata alla stessa latitudine. Il cuore della figura è un diagramma latitudine–tempo (tipo Hovmöller meridiano) che mostra dove e quando compaiono anomalie positive e negative: le aree ombreggiate indicano valori positivi (rafforzamento dei westerlies), mentre le isolinee tratteggiate indicano valori negativi (indebolimento dei westerlies); la linea di zero non è disegnata per rendere più leggibile la struttura. Questa scelta diagnostica è coerente con l’impostazione osservativa classica, in cui la variabilità del flusso zonale medio viene spesso interpretata tramite grandezze integrate e robuste, come il momento angolare relativo o i campi di vento zonale mediato lungo i paralleli, perché sono direttamente legati al bilancio del momento e all’azione degli eddies.

Nel pannello (a) la firma della propagazione è evidente: le bande di anomalia appaiono inclinate, segno che i massimi e i minimi si spostano progressivamente dalle basse latitudini verso le medio-alte latitudini con il passare del tempo. Questa inclinazione è, di fatto, la “traccia” grafica di una propagazione meridionale quasi regolare su scala intrastagionale (settimane). La cosa dinamicamente notevole è che non si tratta di un semplice rinforzo del getto “in loco”, ma di un processo in cui la regione in cui i westerlies risultano anomali tende a rigenerarsi progressivamente più a nord, suggerendo un’interazione organizzata tra onde di Rossby, breaking ondoso e risposta del flusso medio. È proprio questa tipologia di comportamento che, nelle analisi osservazionali, Feldstein aveva evidenziato come propagazione poleward delle anomalie legate al zonal index/momento angolare: la Figura 2a mostra che un modello idealizzato può riprodurre un pattern con caratteristiche analoghe.

Nel pannello (b), invece, l’evoluzione è diversa: le anomalie tendono a disporsi in strutture più “verticali”, cioè compaiono e si attenuano senza un chiaro spostamento latitudinale coerente. Questo è il comportamento tipico della variabilità da zonal index: il getto alterna fasi più intense e più deboli, ma la latitudine caratteristica del massimo rimane relativamente stabile. In termini di dinamica, qui l’interazione con gli eddies agisce soprattutto come modulazione dell’intensità del flusso medio, senza innescare quel ciclo di traslazione meridionale che produce la propagazione verso i poli. Questo tipo di distinzione è del tutto coerente con la letteratura sugli annular modes, dove la variabilità principale del vento zonale medio può essere scomposta in pochi modi dominanti legati a cambi di forza e/o posizione del getto, e dove la persistenza e la forma del segnale dipendono dal feedback fra getto e storm track (cioè dall’organizzazione dei flussi di quantità di moto degli eddies).

La parte destra di ciascun pannello serve proprio a “diagnosticare” questa differenza. La linea spessa rappresenta il profilo medio del vento zonale (la climatologia del getto), mentre le curve sottili mostrano i primi due modi principali di variabilità ricavati statisticamente (EOF). Nel caso propagante (a) i primi due modi hanno importanza comparabile, indicando che il sistema combina in modo significativo variazioni di intensità e spostamenti latitudinali: questa combinazione rende possibile una sequenza spazio-temporale organizzata e quindi la propagazione. Nel caso zonal-index (b), invece, il primo modo domina molto più nettamente del secondo: ciò implica che la variabilità è più “monomodale” e tende a manifestarsi come pulsazione dell’intensità alla stessa latitudine, con una componente di migrazione molto più debole. Questo quadro è in linea con gli studi sul feedback eddy–getto che mostrano come la risposta degli eddies possa selezionare e rinforzare alcuni tipi di anomalie del flusso medio, determinando se il sistema preferisce oscillazioni locali o evoluzioni più traslanti.

Un dettaglio spesso trascurato ma molto utile per la lettura sono anche le linee rette che, come indicato in didascalia, marcano i massimi locali del vento zonale medio: nel pannello (a) queste linee aiutano a visualizzare una traiettoria più ordinata, coerente con una migrazione del core del getto (o, più precisamente, della regione di anomalia zonal-mean associata al getto); nel pannello (b) risultano più disperse e meno coerenti in senso propagante. Questo conferma che la differenza tra i due regimi non è un artefatto della scelta grafica, ma riflette un comportamento fisico distinto della circolazione.

In sintesi, la Figura 2 dimostra che la propagazione verso i poli delle anomalie del vento zonale medio non è un dettaglio secondario, ma un vero e proprio regime dinamico che emerge solo in una parte dello spazio dei parametri del modello. Quando emerge, il sistema mostra un’evoluzione organizzata su scala intrastagionale, potenzialmente associata a maggiore “memoria” e dunque a un incremento della prevedibilità rispetto a un’atmosfera in cui il getto oscilla soprattutto in intensità senza migrare. È esattamente questa distinzione che prepara il terreno alle sezioni successive del lavoro, in cui gli autori attribuiscono la propagazione a un ciclo combinato di propagazione delle onde di Rossby, breaking e rilassamento radiativo.

2. Modello

In questa sottosezione gli autori chiariscono che l’intero studio poggia su un “dynamical core” secco a equazioni primitive (primitive equations), impostato apposta per isolare i meccanismi fondamentali dell’interazione onde–flusso medio senza l’interferenza di processi complessi come condensazione, convezione parametrizzata realistica, orografia e relative onde stazionarie. La scelta è molto in linea con una tradizione consolidata di esperimenti idealizzati: configurazioni con rilassamento termico (Newtonian cooling) verso un profilo di equilibrio prescritto e con smorzamenti/frizioni semplificate sono usate da decenni come banco di prova per studiare la circolazione generale e la fisica dei getti in un contesto controllato. 

Il modello è un globale spettrale con 10 livelli sigma (coordinate verticali di tipo sigma), e una risoluzione orizzontale tipicamente troncata a T30: non è una scelta “povera”, è una scelta funzionale. A questa risoluzione il modello rappresenta bene la dinamica a grande scala e l’attività baroclinica essenziale (storm track, trasporto di quantità di moto degli eddies, modifiche del getto), ma resta abbastanza economico da consentire integrazioni lunghissime (migliaia di giorni), indispensabili quando l’obiettivo è caratterizzare regimi statistici e transizioni in uno spazio dei parametri. 

Il forcing termico è interamente “chiuso” in una temperatura di equilibrio radiativo idealizzata che varia con latitudine e quota sigma. Questa temperatura di riferimento non è un singolo profilo liscio: è costruita come somma di una componente a grande scala (che impone un gradiente meridionale di temperatura e quindi una baroclinicità di base) più due correzioni mirate: un riscaldamento tropicale (H) e un raffreddamento alle alte latitudini (C). In questa famiglia di modelli, H tende a controllare soprattutto l’intensità/struttura del getto subtropicale, mentre C modula l’ampiezza latitudinale e l’efficienza della zona baroclina extratropicale (e quindi la natura dei getti e della waveguide). È esattamente per questo che gli autori esplorano sistematicamente lo spazio dei parametri C/H: variando il bilancio tra “spinta” tropicale e raffreddamento polare si ottengono stati medi con getti differenti e, di conseguenza, differenti modalità di variabilità (propagazione verso i poli vs zonal-index). 

La termodinamica viene implementata tramite un rilassamento della temperatura del modello verso l’equilibrio radiativo con una scala temporale di 30 giorni. Questa scelta è cruciale: da un lato fornisce una sorgente/sink energetica semplice e regolare (quindi “pulita” dal punto di vista interpretativo), dall’altro stabilisce una scala temporale relativamente lenta che permette alla circolazione e agli eddies di organizzarsi e retroagire sul flusso medio, anziché essere continuamente “schiacciati” da un forcing troppo aggressivo. È la stessa filosofia che rende il paradigma Held–Suarez così utile come benchmark dinamico: forcing termico semplice + drag semplice, per far emergere una circolazione statisticamente realistica senza complicare la fisica. 

Sul lato dissipativo, il modello include ingredienti standard ma importanti: attrito superficiale (qui con drag non lineare e scala tipica di circa un giorno), diffusione verticale, e una diffusione orizzontale selettiva di ordine elevato (iper-diffusione) che agisce soprattutto sulle piccole scale, stabilizzando l’integrazione senza smorzare eccessivamente la dinamica sinottica. In un modello spettrale, questo pacchetto è quasi “obbligatorio” per controllare l’accumulo di energia alle alte onde e per ottenere climatologie stabili su integrazioni molto lunghe. 

Infine, dal punto di vista metodologico, la decisione di basare i risultati sugli ultimi 4500 giorni di integrazioni da 5000 giorni è un segnale di rigore: gli autori vogliono lavorare in regime statisticamente stazionario, riducendo il peso dello spin-up e delle transizioni iniziali. Questo è particolarmente rilevante perché l’oggetto dello studio non è un singolo evento, ma l’esistenza (o meno) di un comportamento robusto — la propagazione verso i poli delle anomalie del flusso zonale medio — che emerge solo in certe condizioni di fondo e richiede campioni lunghi per essere diagnosticato con affidabilità. 

3. Robustezza della propagazione verso i poli

In questa parte gli autori fanno un passaggio metodologico decisivo: prima di discutere come avvenga la propagazione verso i poli delle anomalie del flusso zonale medio, verificano quando essa esista davvero e quanto sia sensibile allo stato medio della circolazione nel modello. L’idea è semplice ma potente: se la propagazione è un meccanismo fisico reale (e non un artefatto di una singola simulazione), allora deve emergere in modo coerente in una regione identificabile dello spazio dei parametri che controlla la climatologia del modello. È esattamente ciò che viene fatto esplorando sistematicamente il rapporto tra raffreddamento alle alte latitudini e riscaldamento tropicale, cioè la “manopola” che nel loro framework regola l’assetto dei getti e della baroclinicità di fondo. 

All’interno di questo spazio dei parametri, la variabilità del vento zonale medio in alta troposfera si organizza in due classi. La prima è quella che interessa al paper: un comportamento quasi-periodico in cui le anomalie “nascono” a latitudini più basse e poi migrano ordinatamente verso le alte latitudini (il caso mostrato nella loro Fig. 2a). La seconda classe è la più “classica” nella letteratura sullo zonal index: fluttuazioni del getto e dell’indice zonale che appaiono molto più stocastiche, con anomalie che si intensificano e decadono rimanendo più ancorate alla stessa fascia latitudinale (Fig. 2b). Questa distinzione non è cosmetica: segnala due stati dinamici diversi dell’interazione tra eddies e flusso medio. 

Il risultato chiave sulla “robustezza” è che la propagazione verso i poli non si manifesta ovunque, ma preferenzialmente nella regione in cui il raffreddamento polare è relativamente forte e il riscaldamento tropicale relativamente debole. In quel regime, i westerlies medi sono più deboli e, soprattutto, la struttura meridionale della vorticità potenziale di fondo risulta più ampia e liscia (meno “tagliente”), cioè con un gradiente meno concentrato in una fascia stretta. Al contrario, quando il riscaldamento tropicale è più intenso e il raffreddamento polare più debole, il modello tende a produrre westerlies medi più forti e un gradiente di vorticità potenziale più marcato e affilato; ed è lì che domina la variabilità da zonal index “rumorosa”. In pratica, l’esistenza stessa della propagazione dipende dalle proprietà del campo di fondo che funge da waveguide: se quel waveguide è troppo rigido e concentrato, l’evoluzione tende a rimanere più locale e meno “traslante”; se è più largo e moderato, diventano possibili sequenze spazio-temporali organizzate. 

Questo punto collega direttamente i risultati del paper a una linea di ricerca più ampia. Da un lato, Feldstein aveva mostrato nelle reanalisi che anomalie legate al momento angolare/vento zonale medio possono esibire una propagazione intrastagionale verso i poli, suggerendo che non si tratti di un dettaglio regionale ma di una modalità reale della circolazione. Dall’altro lato, teorie e modelli ridotti avevano già proposto che la direzione della “deriva” delle anomalie dipendesse da come gli eddies riorganizzano la convergenza dei flussi di quantità di moto rispetto alla sorgente d’onda (James e Dodd) oppure da feedback sistematici legati alla baroclinicità e alla storm track (Robinson). La sezione sulla robustezza dà un contenuto operativo a queste idee: la propagazione emerge solo quando lo stato medio permette una particolare configurazione della guida d’onda e quindi del feedback eddy–getto. 

Un altro messaggio forte è di natura “storica”: gli autori sottolineano che il comportamento sorprendentemente periodico visto in esperimenti idealizzati precedenti (come quelli richiamati da James et al.) non è necessariamente un tratto generico dei modelli, ma può essere il prodotto di trovarsi in una specifica porzione dello spazio dei parametri. In altre parole, non bisogna aspettarsi periodicità ovunque: essa appare come proprietà emergente di un regime ben definito del flusso medio. E, cosa ancora più interessante, Son e Lee avevano argomentato che questa regione “propagante” potrebbe non essere lontanissima da condizioni realistiche della circolazione generale, il che rende sensato studiarla in dettaglio anche in chiave di prevedibilità intrastagionale. 

Con questo quadro, il lavoro restringe poi il focus: seleziona una specifica integrazione nel regime propagante e costruisce anche una versione “ridotta” (un modello settoriale) in cui lo spettro dei numeri d’onda zonali viene drasticamente semplificato, mantenendo una componente fondamentale e poche armoniche. Lo scopo non è estetico ma causale: se la propagazione sopravvive anche quando si limitano i gradi di libertà ondosi, allora il fenomeno può essere attribuito a meccanismi essenziali (propagazione/assorbimento d’onda e risposta del flusso medio) più che a una miscela complicata di interazioni multi-scala. Questa scelta prepara il terreno alle sezioni successive, dove la propagazione viene interpretata come sequenza auto-consistente di propagazione delle onde di Rossby, breaking e “ripristino” del flusso medio nelle basse latitudini. 

Figura 3 – Sintesi delle statistiche del flusso nello spazio dei parametri C/H

La Figura 3 è, di fatto, la “mappa di regime” dell’articolo: mostra come cambiano struttura del jet, tipo di variabilità del vento zonale medio e configurazione del gradiente di vorticità potenziale (PV) in alta troposfera al variare dei due ingredienti termici imposti nel modello, cioè il raffreddamento alle alte latitudini (C) sull’asse orizzontale e il riscaldamento tropicale (H) sull’asse verticale. L’idea fisica è che, cambiando C e H, cambi l’assetto medio dei westerlies e della baroclinicità, e quindi cambi anche “come” gli eddies riescono a forzare e organizzare le anomalie del flusso zonale medio. Questa impostazione deriva direttamente dal lavoro di Son e Lee sulla risposta dei getti al forcing termico in un modello a equazioni primitive, dove la combinazione tra forcing tropicale e polare determina la comparsa di stati a getto singolo o doppio getto e ne regola intensità e latitudine. 

Nel pannello (a) (“Jet Structure”) ogni punto dello spazio C–H riassume com’è fatto il getto medio: i cerchi indicano la posizione del massimo del vento zonale medio (la latitudine è scritta sotto ciascun cerchio) e i numeri dentro i cerchi indicano la velocità massima del jet. Il simbolo è informativo anche sulla morfologia: cerchio singolo = stato a getto singolo, cerchio doppio = stato a doppio getto (tipicamente una componente subtropicale e una extratropicale separate). Le due simulazioni evidenziate (cerchi scuri) sono quelle usate nella Figura 2 per mostrare, rispettivamente, un caso “propagante” e un caso “zonal index”. Questo pannello va letto come una diagnosi del “campo di base” su cui poi agisce l’interazione onde–flusso medio: un jet più intenso e più concentrato implica in genere una waveguide più rigida; un jet più debole e/o più largo implica una waveguide meno “tagliente” e potenzialmente più favorevole a riassorbimento/breaking e ri-organizzazione meridionale dell’attività d’onda. 

Il pannello (b) (“Summary of variability”) è la conseguenza dinamica: in alto a sinistra compare il regime in cui domina lo zonal index / jet meander, cioè una variabilità principalmente stocastica con oscillazioni del getto che restano in gran parte ancorate alla stessa fascia latitudinale; in basso a destra emerge invece la regione in cui prevale la propagazione verso i poli delle anomalie del vento zonale medio, con un comportamento più quasi-periodico (come in Figura 2a). Tra i due c’è una fascia di transizione. Il punto cruciale è che la propagazione non è “sempre presente”: è un regime che richiede un certo background. Questo messaggio è coerente con quanto si osserva anche nei dati: Feldstein aveva mostrato che la propagazione verso i poli delle anomalie zonal-mean (diagnosticata tramite momento angolare relativo) è un fenomeno reale ma non necessariamente onnipresente, e la letteratura successiva sottolinea proprio la dipendenza dal regime di fondo e dalla struttura dei modi dominanti. 

Il pannello (c) spiega il “perché” in termini di PV/waveguide: non sta dicendo solo che “cambia la PV”, ma che cambia la larghezza meridionale della regione in cui il gradiente di PV in alta troposfera funge da guida per le onde di Rossby. La zona ombreggiata (larghezze grandi) corrisponde, in modo molto suggestivo, alla regione del pannello (b) in cui compare la propagazione verso i poli. Tradotto: quando il gradiente di PV è relativamente ampio e più graduale, il sistema favorisce una sequenza meridionale organizzata (le anomalie possono “spostare” il luogo di assorbimento/breaking e rigenerarsi più a nord); quando invece il gradiente è più stretto e più intenso, le onde tendono a essere più intrappolate e la variabilità resta più locale, producendo il comportamento da zonal index. È esattamente la logica del paper: la propagazione emerge dall’interazione tra propagazione d’onda, breaking e risposta del flusso medio, ma può farlo solo se il background non rende la waveguide troppo rigida. 

Nel complesso, la Figura 3 “chiude il cerchio” tra climatologia e variabilità: mostra che il modello passa da un mondo in cui il getto è forte e la variabilità è rumorosa (zonal index) a un mondo in cui il getto è più debole/meno affilato e la variabilità diventa organizzata e propagante, e indica chiaramente che la discriminante fisica non è una singola variabile locale, ma la struttura meridionale di fondo che controlla la propagazione e l’assorbimento delle onde. Questo tipo di lettura si incastra bene anche con la visione più ampia della variabilità del momento angolare atmosferico e dei modi annulari, dove la “mobilità” del sistema dipende in modo critico da come eddies e flusso medio si accoppiano lungo il meridiano.

4. Procedura di analisi dell’output del modello

In questa fase gli autori definiscono un protocollo di analisi pensato per “isolare” in modo oggettivo gli intervalli in cui la propagazione verso i poli è più netta e poi descriverne la struttura media. L’idea di fondo è tipica degli studi sulla variabilità intrastagionale: quando il segnale è intermittente e immerso in rumore sinottico, conviene identificare episodi robusti con un indice semplice e poi usare tecniche compositive per aumentare il rapporto segnale/rumore, evitando che la media cancelli strutture coerenti. Questo approccio è ampiamente discusso nella manualistica statistica per le scienze atmosferiche, dove compositi e selezioni “event-based” vengono presentati come strumenti standard per evidenziare pattern ricorrenti e le loro evoluzioni temporali. 

Poiché nel loro caso la propagazione è visivamente più chiara tra equatore e 50° di latitudine, scelgono come indice operativo l’anomalia del vento zonale medio a 25°, cioè la latitudine intermedia del “corridoio” in cui il segnale è più leggibile (come mostrato nella loro Figura 2a). La scelta è metodologicamente sensata perché un indice in posizione centrale tende a catturare sia la fase di innesco alle basse latitudini sia la fase in cui la risposta massima si organizza sul getto guidato dagli eddies, che nel loro stato medio si colloca intorno a 50°N. La Figura 2a suggerisce inoltre una relazione di fase coerente: quando l’indice a 25° raggiunge i valori più negativi, le anomalie a circa 50°N tendono a risultare massime e positive, segnalando una struttura meridionalmente organizzata (non un semplice rafforzamento locale del getto). Questo tipo di lettura è in continuità con la diagnostica osservativa della propagazione intrastagionale del flusso zonale medio e del momento angolare relativo, dove la migrazione meridionale delle anomalie viene riconosciuta proprio da relazioni di fase tra fasce latitudinali diverse. 

La selezione degli episodi avviene con una regola chiara e replicabile: si cercano i minimi locali dell’indice che scendono sotto una soglia pari a una deviazione standard negativa (quindi si prendono solo gli eventi “forti”), e si impone che due episodi consecutivi siano separati da almeno 100 giorni, così da escludere sequenze anomale troppo ravvicinate rispetto al periodo tipico del fenomeno. Questa doppia condizione ha due implicazioni importanti: da un lato riduce il rischio di includere fluttuazioni brevi e non rappresentative; dall’altro rende più probabile che gli episodi selezionati siano statisticamente indipendenti o quantomeno non semplici “doppioni” della stessa fase. Nel loro set di simulazioni ciò porta a identificare 21 episodi nell’emisfero nord del modello e 22 nell’emisfero sud; considerando che il periodo tipico della propagazione nel modello è dell’ordine di qualche mese (circa 160 giorni), questo conteggio copre una porzione molto ampia della serie analizzata (circa tre quarti). Il punto, qui, è cruciale: non stanno inseguendo pochi casi rari, ma un comportamento che nel regime propagante diventa quasi la normalità del sistema. 

Una volta definiti gli episodi, costruiscono la “firma media” del fenomeno con un composito centrato sul giorno mediano di ciascun episodio, definito come giorno a ritardo zero. È la classica logica della composizione a ritardi: si allineano gli eventi su un riferimento comune e si mediano i campi a vari ritardi temporali prima e dopo, ottenendo una sequenza evolutiva tipica. Il vantaggio è che, se esiste un’evoluzione coerente (ad esempio l’innesco ai tropici e la successiva traslazione verso latitudini più alte), questa emerge con molta più chiarezza rispetto a una singola realizzazione dominata dalla variabilità sinottica. 

Per i campi che variano anche lungo la longitudine (quindi non solo grandezze zonal-mean), gli autori aggiungono un passaggio fondamentale: l’allineamento di fase in longitudine, proprio per evitare che onde simili ma “posizionate” in longitudes diverse finiscano per annullarsi nella media (il classico problema dell’interferenza distruttiva nei compositi). Operativamente, per ogni episodio identificano, nel giorno a ritardo zero e lungo la latitudine di riferimento (25°N), la longitudine in cui l’indice raggiunge il minimo; quella longitudine diventa il riferimento dell’episodio. Poi traslano longitudinalmente tutti i campi di quell’episodio (per tutti i ritardi e tutte le latitudini) in modo che il riferimento coincida sempre con una longitudine comune scelta per il composito. In questo modo le strutture ondose associate alla propagazione vengono “messe in fase” prima della media, e ciò consente di ricostruire un quadro composito che preserva geometria e segnale dinamico, invece di ridurli a un campo smorzato. Questa è una scelta metodologica molto coerente con l’obiettivo del lavoro, perché la tesi degli autori è che la propagazione verso i poli derivi da una sequenza organizzata di propagazione d’onda, assorbimento/breaking e risposta del flusso medio: se non allinei correttamente la fase, rischi di perdere proprio la parte ondosa che alimenta il meccanismo. 

Infine, presentano i compositi costruiti sugli episodi dell’emisfero nord, specificando che quelli dell’emisfero sud forniscono sostanzialmente lo stesso quadro. In un modello idealizzato, secco e senza orografia, questa quasi-simmetria emisferica è attesa e rafforza l’interpretazione “meccanicistica” del fenomeno: la propagazione non dipende da dettagli regionali, ma da proprietà generali dello stato medio e dell’interazione tra eddies e flusso zonale medio, in linea con la letteratura che inquadra la variabilità del momento angolare e dei modi annulari come risultato di accoppiamenti dinamici su scala planetaria.

5. Risultati dal modello completo

In questa sezione gli autori presentano una serie di compositi (costruiti sugli episodi di propagazione verso i poli selezionati in precedenza) e usano il campo delle anomalie del vento zonale medio a 250 hPa come “sfondo” di riferimento nelle Figure 4a, 4b, 4d e 4e. In tal modo, la lettura degli altri campi compositi avviene sempre in relazione alla fase della propagazione mostrata dal vento zonale medio. Gli autori specificano inoltre che i campi in Figg. 4a–e presentano ampie regioni statisticamente significative al 95% secondo un t test bilaterale, informazione utile perché segnala che le strutture discusse non sono limitate a poche aree puntuali ma risultano estese nei compositi. 

Un primo risultato descritto riguarda la relazione temporale tra le anomalie del vento zonale medio e la convergenza del flusso di quantità di moto degli eddies a 250 hPa: nel composito, la convergenza del flusso di momento tende a precedere leggermente l’anomalia del vento zonale medio. Gli autori notano che questa sfasatura è coerente con l’evidenza osservativa riportata da Feldstein sull’evoluzione intrastagionale del flusso zonale medio e del momento angolare relativo zonal-mean, dove i termini eddy (in particolare i flussi di momento) mostrano un ruolo temporale anticipatore rispetto alla risposta del vento zonale medio. Questa impostazione (diagnostica del bilancio del momento e confronto lead–lag) è anche in linea con la letteratura sugli annular modes e sul feedback eddy–getto, che utilizza proprio la convergenza dei flussi di momento per descrivere la covariazione tra storm track e vento zonale medio. 

Per caratterizzare le sorgenti di attività ondosa e i processi energetici associati, gli autori analizzano anche le anomalie del flusso di calore degli eddies a 850 hPa, livello al quale (nel loro modello) questo flusso è massimo. Nel composito, il flusso di calore mostra un segnale di propagazione verso i poli, ma in modo chiaramente confinato a una fascia più ristretta, grosso modo tra 20° e 50° di latitudine, quindi più limitata rispetto all’estensione meridionale con cui compaiono le anomalie del vento zonale medio e della convergenza del flusso di momento. Nella descrizione della Fig. 4c gli autori evidenziano anche che, oltre alla componente che varia con la stessa periodicità del ciclo di propagazione, compaiono fluttuazioni a scala temporale più breve: tra lag −40 e lag 0 giorni e intorno a 40° di latitudine si riconoscono tre massimi distinti nel flusso di calore, ciascuno con una scala dell’ordine di ~10 giorni, sovrapposti alla modulazione a più lunga scala temporale. 

Un ulteriore elemento riportato nasce dalla sovrapposizione (nella Fig. 4c) dei contorni della convergenza del flusso di momento sopra il campo (ombreggiato) del flusso di calore: durante le fasi di flusso di calore aumentato, viene descritta l’emersione di una nuova banda di anomalia negativa della convergenza del flusso di momento in area tropicale, mentre un’altra banda negativa già presente si sposta verso la latitudine tipica del getto del fronte polare nel modello (nell’intervallo di lag compreso tra −40 e 0 giorni). Durante le fasi di flusso di calore ridotto viene descritto il comportamento opposto. Nello stesso passaggio gli autori riportano che la relazione di fase tra i due flussi non è uniforme con la latitudine: nei subtropici la convergenza del flusso di momento tende a precedere il flusso di calore, mentre alle medie latitudini accade l’inverso, soprattutto per i segnali a scala di ~10 giorni identificati nel periodo tra lag −40 e 0. 

Infine, viene indicato un dettaglio temporale aggiuntivo: all’inizio della fase di propagazione (in area tropicale e subtropicale) le anomalie della convergenza del flusso di momento compaiono molto prima delle anomalie del flusso di calore, mentre per le strutture a scala di ~10 giorni il flusso di calore tende a precedere di alcuni giorni la convergenza del flusso di momento. Gli autori segnalano che queste relazioni tra flussi eddy verranno riprese più approfonditamente nella sezione successiva, mantenendo qui il focus sulla descrizione dei pattern compositi nel modello completo. 

gli autori mettono in fila due “strade” possibili per l’innesco della propagazione verso i poli delle anomalie del vento zonale medio e chiariscono, con diagnostiche mirate, quale di queste strade sia compatibile con il loro modello a spettro completo. Da un lato, riconoscono che in atmosfera reale non si può escludere un avvio “tropicale”: un’anomalia di riscaldamento convettivo ai tropici può intensificare (o indebolire) la circolazione di Hadley e produrre, alle basse latitudini, un’anomalia coerente del vento zonale medio. A supporto di questa possibilità richiamano risultati osservativi che legano la propagazione meridionale del momento angolare atmosferico (e quindi del flusso zonale medio) a variabilità tropicale organizzata: Feldstein (1999) mostra una propagazione verso i poli della convergenza dei flussi di momento angolare quando i campi vengono regressi rispetto alla tendenza della durata del giorno (LOD), grandezza che risulta fortemente connessa alla variabilità intrastagionale tropicale associata alla Madden–Julian Oscillation. Su scale più lunghe, Dickey et al. (1992) documentano una propagazione coerente verso i poli delle fluttuazioni interannuali del momento angolare atmosferico, in fase con la variabilità ENSO, suggerendo un possibile ruolo di anomalie convettive tropicali come “trigger” su quelle frequenze. 

Dall’altro lato, gli autori verificano se nel loro sistema idealizzato la sequenza sia effettivamente quella “Hadley-driven” che ci si aspetterebbe in un innesco tropicale: in quel caso, la banda di vento zonale medio positivo alle basse latitudini dovrebbe collocarsi al bordo polare di un’anomalia di flusso divergente meridionale coerente con una Hadley anormalmente intensa. Nelle loro figure (in particolare il composito del vento divergente e quello del riscaldamento diabatico), il quadro descritto è invece compatibile con una circolazione meridionale che emerge come risposta ai flussi di calore e di quantità di moto degli eddies, e con anomalie di riscaldamento/raffreddamento che risultano allineate alle regioni di moto verticale implicate da tale circolazione. Questo comportamento, nel contesto del loro modello, è anche “atteso” perché il riscaldamento diabatico non è parametrizzato tramite convezione realistica, ma nasce dal rilassamento newtoniano verso un equilibrio radiativo, che finisce per riflettere soprattutto la dinamica eddy-driven. In termini operativi, ciò che gli autori evidenziano è che anche se nella realtà la forzante tropicale può avviare episodi di propagazione, nel loro setup la firma tropicale (vento divergente e riscaldamento diabatico) può essere prodotta dal lato extratropicale attraverso la catena eddies → circolazioni indotte → risposta termica, e questa stessa firma diventa quindi un possibile “marcatore” diagnostico da cercare in analisi future su dati osservativi. 

Per risalire a quali processi generino la convergenza dei flussi di quantità di moto degli eddies, gli autori spostano poi l’attenzione sulla struttura orizzontale della vorticità potenziale, usando la PV di Ertel su 250 hPa e sovrapponendo le anomalie locali del vento zonale. Qui la logica è quella classica della diagnostica PV: anche se la PV “più elegante” sarebbe su superfici isentropiche, le mappe PV restano uno strumento estremamente informativo per riconoscere strutture d’onda, filamenti, overturning e mixing e collegarli a segnali cinematici e dinamici nella troposfera in coordinate di pressione. La descrizione che forniscono è che il contributo al segnale zonal-mean (cioè all’anomalia del vento zonale medio) proviene tipicamente da disturbi baroclinici localizzati longitudinalmente, non da un pattern perfettamente zonalmente uniforme: è un dettaglio importante perché spiega perché, nella media zonale, emergano strutture meridionali a tripolo o quadrupolo. In questa lettura, le anomalie positive e negative del vento zonale medio si associano alle regioni in cui i contorni di PV vengono intensificati, ribaltati o “spezzati” (wave breaking), cioè a quella fase non lineare del ciclo d’onda in cui il gradiente di PV viene localmente deformato e si innescano processi di mescolamento quasi-irreversibile. Questa cornice concettuale è coerente con la letteratura fondativa sul Rossby wave breaking e sul mescolamento della PV, dove l’overturning del gradiente meridionale di PV e la successiva fase di mixing sono visti come ingredienti essenziali per la riorganizzazione del flusso medio e per la creazione di regioni tipo “surf zone” o strati critici non lineari. 

Nella sequenza descritta per i tropici (attorno alla latitudine di riferimento), l’anomalia negativa del vento zonale medio viene fatta corrispondere a una fase in cui la struttura d’onda passa da “ripida” a “overturning” fino al massimo ribaltamento, e poi decade quando il mescolamento prende il sopravvento. Poiché su queste scale temporali la PV è quasi conservata lungo le traiettorie (al netto dei termini diabatici/frizionali), gli autori richiamano implicitamente un vincolo geometrico: se la PV viene omogeneizzata dentro la regione di breaking, i bordi nord e sud di quella regione tendono a presentare un gradiente di PV più marcato, il che è compatibile con la formazione di anomalie di segno opposto ai lati dell’anomalia centrale e quindi con la comparsa di tripoli meridionali. Questo tipo di “impronta” (centro di segno opposto flanqueato da due lobi del segno contrario) è una firma ricorrente quando il flusso medio risponde a un’area di mescolamento PV indotta dal breaking, come discusso in molti contesti troposferici e stratosferici. 

Infine, sulla base del fatto che nei tropici profondi la struttura dei flussi eddy (in particolare la divergenza/convergenza del flusso di quantità di moto) viene associata a onde che arrivano dalle medie latitudini e poi vanno incontro a steepening e breaking, gli autori introducono l’idea che la dinamica di “latitudine critica” possa essere rilevante: nella teoria delle onde di Rossby, l’avvicinamento a regioni in cui il flusso medio e la velocità di fase si eguagliano favorisce assorbimento e rottura d’onda, spesso in forma non lineare, con forte mixing PV e trasferimenti di quantità di moto. Tuttavia evidenziano anche un limite pratico: con uno spettro completo di onde zonali, stimare in modo pulito velocità di fase e latitudini critiche (e attribuire il segnale a specifiche componenti d’onda) non è banale; per questo motivano il passaggio a un modello settoriale a spettro ridotto, che può rendere più trasparente la diagnosi dei meccanismi senza perdere, idealmente, i processi dinamici portanti del caso completo. 

La Figura 4 è costruita come un grande diagramma “tempo–latitudine” (asse orizzontale = lag in giorni, asse verticale = latitudine) e serve a visualizzare, nello stesso riferimento composito, la sequenza tipica di un episodio di propagazione verso i poli. Nei pannelli (a), (b), (d) ed (e) lo sfondo in ombreggiatura è sempre l’anomalia del vento zonale medio a 250 hPa: le bande inclinate che risalgono da basse verso medio-alte latitudini sono la firma grafica della propagazione meridionale del segnale nel modello completo. Gli autori indicano anche che molte strutture in questi compositi superano il 95% di significatività (t test bilaterale), cioè che non si tratta di dettagli puntuali ma di pattern estesi e coerenti nel campione di episodi. 

Sopra questa “traccia” del vento zonale medio, nel pannello (a) compaiono i contorni della convergenza del flusso di quantità di moto degli eddy in alta troposfera (250 hPa). La caratteristica chiave è la relazione di fase: le massime aree di convergenza/divergenza tendono a comparire poco prima che l’anomalia del vento zonale medio raggiunga la massima ampiezza alla stessa latitudine. È esattamente il tipo di lead-lag atteso quando la variabilità del vento zonale medio è legata al bilancio del momento zonale: gli eddy trasferiscono quantità di moto meridionalmente e, dove il trasporto converge, il flusso medio viene accelerato; dove diverge, viene decelerato. Questa lettura è coerente sia con le diagnosi classiche basate su flussi alla Eliassen–Palm/eddy-mean flow, sia con la letteratura sui modi annulari e sul feedback eddy-getto, dove la convergenza dei flussi di momento è un indicatore diretto del forcing eddy sul getto. 

Nel pannello (b) i contorni rappresentano il flusso di calore degli eddy a 850 hPa (il livello in cui nel loro modello questo flusso è massimo), mentre l’ombreggiatura resta l’anomalia del vento zonale medio a 250 hPa. Qui si vede un dettaglio importante: il segnale “propagante” del flusso di calore è più evidente soprattutto tra ~20° e 50°, quindi in una fascia che ricorda la zona di attività baroclinica/storm track. Nello stesso tempo, la struttura del flusso di calore mostra una componente a scala più lunga (in fase con il ciclo di propagazione) su cui si innestano pulsazioni più brevi: gli autori notano tre massimi distinti tra lag −40 e 0 attorno a ~40°, con tempi tipici dell’ordine di ~10 giorni, cioè una scala sinottica. In termini di diagnostica energetica, il flusso di calore meridionale è spesso usato come proxy dell’attività baroclinica e della conversione di energia potenziale disponibile dal flusso medio agli eddy, per cui una modulazione a più scale (intrastagionale + sinottica) è fisicamente plausibile in un regime dominato da instabilità baroclina. 

Il pannello (c) rende più leggibile il rapporto tra i due flussi: qui l’ombreggiatura è il flusso di calore a 850 hPa, mentre i contorni sono di nuovo la convergenza del flusso di quantità di moto in alta troposfera. La sovrapposizione mette in evidenza che le due grandezze non sono perfettamente in fase ovunque: in alcuni tratti (soprattutto ai tropici/subtropici) la convergenza del flusso di momento tende a comparire prima del segnale termico, mentre alle medie latitudini — in particolare per le pulsazioni più rapide — il flusso di calore può precedere di qualche giorno la risposta nel flusso di momento. Questo tipo di sfasamento latitudine-dipendente è tipico quando coesistono processi diversi (sviluppo baroclino, crescita e propagazione di onde, aggiustamento del getto) che non hanno la stessa cronologia lungo il meridiano. 

I pannelli (d) ed (e) aggiungono la parte “secondaria” della risposta media. Nel (d) i contorni descrivono il vento meridionale divergente a 250 hPa: è un modo compatto per visualizzare la componente divergente della circolazione meridionale associata all’episodio (quindi una cella/overturning composita). Nel (e) i contorni mostrano le anomalie di riscaldamento diabatico, che nel loro modello derivano dal rilassamento newtoniano verso l’equilibrio radiativo. La co-locazione tra aree di riscaldamento/raffreddamento e le regioni di moto verticale implicate dal campo divergente è un controllo di coerenza fisica: indica che la parte termica del segnale è “organizzata” insieme alla circolazione meridionale composita e non appare come rumore scollegato. In un modello in cui il riscaldamento diabatico non è guidato da convezione interattiva realistica ma da rilassamento, è naturale che le anomalie termiche risultino fortemente legate alle circolazioni indotte dagli eddy e ai flussi che le sostengono. 

Letta nel suo insieme, la Figura 4 non è solo una conferma visiva della propagazione del vento zonale medio: è una diagnosi “multicampo” che permette di seguire, lungo la stessa coordinata tempo-latitudine, come si dispongono (e con quale sfasamento) i termini chiave dell’interazione eddy–flusso medio (flussi di momento e di calore), insieme alla componente divergente della circolazione e alla risposta termica. Questo tipo di rappresentazione è strettamente in continuità con l’impostazione osservativa che aveva motivato il problema (propagazione intrastagionale del flusso zonale medio/momento angolare) e con le cornici teoriche usate per interpretare la variabilità dei getti come esito di forcing eddy e feedback eddy-getto. 

La Figura 5 è costruita per far vedere che ciò che, nelle figure tempo–latitudine, appare come un’anomalia “zonal-mean” che si organizza e migra verso i poli, in realtà nasce da disturbi ondosi baroclinici longitudinalmente localizzati e dalla loro evoluzione non lineare in termini di vorticità potenziale (PV) in alta troposfera. I quattro pannelli mostrano compositi della PV di Ertel a 250 hPa ai lag −16, −8, 0 e +8 giorni; nel pannello centrale di ciascun lag le linee sottili sono i contorni di PV (intervallo più largo), mentre i contorni spessi rappresentano le anomalie locali del vento zonale a 250 hPa (linea continua positiva, tratteggiata negativa). Il riquadro evidenzia l’area in cui il contributo alle anomalie zonal-mean è massimo, e l’inserto a sinistra è un ingrandimento della stessa regione con contorni PV più fitti; sul margine destro compare, come riferimento, il profilo latitudinale dell’anomalia del vento zonale mediata lungo la longitudine (cioè il “pezzo” che entra nel segnale zonal-mean). Questa scelta diagnostica (PV su una superficie isobarica) è deliberata: anche se la PV è concettualmente più “naturale” su superfici isentropiche, le mappe PV restano estremamente efficaci per riconoscere strutture d’onda, filamenti, ribaltamenti e mescolamento, e qui servono soprattutto come guida cinematica per interpretare i flussi eddy che forzano il vento zonale medio studiato a 250 hPa. 

Guardando la sequenza, a lag −16 si osserva una perturbazione ancora relativamente “ordinata”: i contorni di PV mostrano un’ondulazione ben definita e una fase di steepening (addensamento e inclinazione dei contorni) nella regione evidenziata, coerente con uno stadio di crescita baroclinica in cui il gradiente di PV associato al getto viene deformato e intensificato localmente; le anomalie di vento zonale (contorni spessi) sono già organizzate in lobi di segno alterno, e il profilo a destra suggerisce una struttura meridionale che può proiettare su configurazioni tipo tripolo/quadrupolo quando si fa la media zonale. A lag −8 l’onda entra più chiaramente nello stadio non lineare: nel riquadro (e ancor meglio nell’inserto) i contorni di PV cominciano a curvarsi e ad avvicinarsi fino a configurazioni che preludono al ribaltamento (overturning), cioè alla fase in cui l’onda non si limita più a propagare ma inizia a riorganizzare in modo irreversibile il gradiente meridionale di PV. Questo tipo di transizione è un tratto classico dei cicli di vita delle onde barocliniche, dove l’evoluzione tarda può culminare in forme di wave breaking (con geometrie cicloniche o anticicloniche) e in un forte rimaneggiamento della PV in alta troposfera. 

Il pannello a lag 0 è il punto di massima “firma” dinamica: nel riquadro i contorni di PV risultano fortemente ribaltati e avvolti, con un aspetto tipico del Rossby wave breaking, cioè della fase in cui il gradiente di PV viene localmente rovesciato e si attivano processi di mescolamento che tendono a omogeneizzare la PV all’interno della regione di breaking (“surf zone”), lasciando però gradienti più intensi ai bordi della zona mescolata. Questa è esattamente la fenomenologia descritta nella letteratura classica di McIntyre e Palmer sul breaking e sul mescolamento di PV: l’overturning produce filamenti e intrusione reciproca di PV, e il mixing conseguente modifica in modo netto la struttura del flusso medio perché la PV, su scale temporali sinottiche, è quasi conservata e quindi il rimaneggiamento deve “chiudere i conti” rinforzando i gradienti ai margini della regione rimescolata. In Figura 5 questo si riflette nel fatto che l’anomalia di vento zonale associata al breaking non è un singolo nucleo: i lobi di vento (continui e tratteggiati) si dispongono in modo da produrre, nella media zonale mostrata a destra, una sequenza di segni alterni coerente con la proiezione tripolare/quadrupolare discussa nel testo. 

A lag +8 la struttura nel riquadro appare più “rilassata”: i contorni di PV sono meno arrotolati e più regolari, segno che la fase più intensa di overturning è passata e che l’anomalia si sta riorganizzando/decadendo; i contorni spessi del vento zonale mostrano ancora lobi coerenti ma con un assetto latitudinale spostato rispetto ai pannelli precedenti, e il profilo a destra restituisce la corrispondente variazione del contributo zonal-mean. Il messaggio che la figura rende molto concreto è che, in ogni istante, la componente zonal-mean dell’anomalia (quella che poi “sembra propagare” nelle figure lag–latitudine) deriva soprattutto da una porzione longitudinale dove l’onda è in fase di forte deformazione e/o breaking: fuori da quella finestra, le anomalie di PV sono molto più deboli e contribuirebbero poco alla media zonale. Questo è coerente con l’impostazione stessa dell’articolo, che usa la Figura 5 per legare i termini di convergenza dei flussi eddy (mostrati in Fig. 4) a una dinamica di crescita–steepening–overturning–mixing di onde che, pur originate alle medie latitudini, possono produrre firme marcate anche alle basse latitudini quando incontrano condizioni favorevoli al breaking e all’assorbimento. 

La Figura 6 ripropone la stessa diagnostica della Figura 4, ma applicata al sector model (modello “a settore”), cioè una versione ridotta in cui lo spettro delle onde zonali è limitato (nel lavoro, alla fondamentale e a poche armoniche), così da rendere più “leggibile” la dinamica ondosa che alimenta la propagazione verso i poli. 

La lettura è la seguente: asse orizzontale = lag (giorni) rispetto al giorno centrale dell’episodio, asse verticale = latitudine. Nei pannelli (a), (b), (d) ed (e) l’ombreggiatura è sempre l’anomalia del vento zonale medio a 250 hPa, quindi le bande inclinate sono la traccia spazio-temporale della propagazione meridionale del segnale; i contorni, invece, cambiano da pannello a pannello e permettono di confrontare la fase della propagazione con i termini eddy e con la risposta “secondaria” della circolazione. 

Nel pannello (a) i contorni rappresentano la convergenza del flusso di quantità di moto degli eddies in alta troposfera (250 hPa). La caratteristica che la figura rende evidente è la coerenza tra le strutture di convergenza/divergenza del flusso di momento e le bande di vento zonale medio: le anomalie dei flussi eddy si dispongono lungo traiettorie inclinate simili a quelle del vento, e il loro massimo tende a collocarsi leggermente “a monte” in fase rispetto al picco del vento zonale medio alla stessa latitudine, come ci si aspetta in un quadro di interazione onde–flusso medio in cui i flussi di momento forzano l’accelerazione/decelerazione del getto. Questa impostazione diagnostica è direttamente in continuità con la letteratura su flussi eddy, feedback eddy-getto e formulazioni tipo Eliassen–Palm per la retroazione delle onde sul flusso medio. 

Nel pannello (b) i contorni mostrano il flusso di calore degli eddies a 850 hPa, che nel modello è un proxy dell’attività baroclina a bassa troposfera (la componente sinottica/storm-track). Qui il segnale appare più confinato, soprattutto tra fascia subtropicale e medie latitudini, e mostra pulsazioni più “granulari” nel tempo rispetto al pannello (a), coerenti con la natura episodica dell’attività baroclina. La figura consente quindi di vedere, nello stesso riferimento lag-latitudine, come la modulazione baroclina (flusso di calore) e la riorganizzazione del flusso zonale medio (bande ombreggiate) coesistano ma non abbiano necessariamente la stessa estensione meridionale né la stessa regolarità temporale. 

Il pannello (c) è costruito proprio per chiarire la relazione tra i due “motori” eddy: l’ombreggiatura rappresenta il flusso di calore a 850 hPa, mentre i contorni tornano a essere la convergenza del flusso di quantità di moto in alta troposfera. La sovrapposizione fa emergere che i massimi di flusso di calore e quelli di convergenza del flusso di momento non risultano ovunque simultanei: in alcune fasce latitudinali la firma del trasferimento di momento in quota è più precoce rispetto al segnale termico, mentre altrove le pulsazioni del flusso di calore risultano più avanzate rispetto alla risposta nel flusso di momento. Questa asimmetria di fase è compatibile con l’idea (ben nota nelle diagnosi eddy-mean flow) che i termini associati alla crescita baroclina e quelli associati al trasporto/assorbimento del momento possano organizzarsi con tempistiche diverse lungo il meridiano e lungo il ciclo di vita delle onde. 

Nel pannello (d) i contorni mostrano il vento meridionale divergente a 250 hPa, cioè una misura della circolazione meridionale “secondaria” associata agli episodi compositati. La struttura segue la stessa geometria inclinata del segnale di vento zonale medio, indicando che, nel composito del sector model, esiste una risposta organizzata della componente divergente che accompagna la migrazione meridionale delle anomalie del getto. In chiave interpretativa generale (senza entrare qui nei dettagli del bilancio), questo tipo di campo viene spesso usato per distinguere una risposta meridionale coerente indotta dai flussi eddy rispetto a un forcing puramente imposto dall’esterno, ed è esattamente il motivo per cui nel lavoro viene affiancato al pannello (e). 

Nel pannello (e), infine, i contorni rappresentano il riscaldamento diabatico composito. Nel contesto del modello utilizzato dagli autori, questa variabile è strettamente legata allo schema di rilassamento radiativo (Newtonian relaxation) e quindi tende a organizzarsi in modo consistente con le anomalie di circolazione divergente e con la struttura meridionale del segnale di vento zonale medio. La figura mostra dunque, nello stesso diagramma lag-latitudine, l’allineamento tra la componente termica (riscaldamento/raffreddamento) e la componente dinamica (circolazione divergente) durante la sequenza propagante nel sector model, fornendo una lettura integrata della parte “termodinamica” che accompagna i flussi eddy e l’anomalia del getto. 

La propagazione verso il polo nel modello a settore: fedeltà al modello completo, vincoli spettrali e implicazioni dinamiche

Per valutare quanto la propagazione meridiana verso il polo (poleward propagation) sia un risultato “robusto” e legato a pochi ingredienti essenziali della dinamica, viene costruita una versione a settore del modello a spettro completo, scegliendo come onda fondamentale il numero d’onda zonale più energetico nella simulazione full-spectrum. In questo caso emerge che k = 5 domina il contenuto energetico; per non sterilizzare il trasferimento spettrale e consentire una cascata di enstrofia, si mantengono inoltre alcune armoniche (un punto cruciale nei modelli ideali, dove la dinamica turbolenta tende a redistribuire varianza e vorticità su più scale). L’integrazione prolungata (migliaia di giorni) e l’analisi su un sottoinsieme “stazionario” dell’output servono a isolare il segnale intrinseco del sistema, minimizzando la dipendenza da transitori iniziali.

Il risultato centrale è che la propagazione verso il polo del run di controllo viene riprodotta molto bene dal modello a settore: i campi compositi associati all’evento (quelli che, nel lavoro originale, sono mostrati nelle Figg. 6a–e) risultano molto simili ai corrispettivi del modello completo, suggerendo che il settore cattura i processi fisici di base della propagazione. Questo è un passaggio concettualmente importante: se una versione drasticamente semplificata sul piano zonale conserva la fenomenologia principale, allora la propagazione non dipende in modo determinante da dettagli “accessori” della variabilità tridimensionale, ma da un meccanismo più fondamentale di interazione onda–flusso medio e di feedback tra flusso zonale e forzanti turbolente. In termini più generali, questo tipo di comportamento si collega alla letteratura sulla propagazione meridiana delle anomalie zonali e sul ruolo dei feedback eddy-driven nella modulazione e migrazione dei getti. 

Un elemento particolarmente diagnostico è la relazione di fase tra flusso di calore turbolento (eddy heat flux, legato alla conversione baroclina e al trasporto meridiano di calore) e convergenza del flusso di quantità di moto turbolento(eddy momentum flux convergence, legata alla forzante barotropa sul getto). Il fatto che, alle medie latitudini e su scale relativamente brevi, il flusso di calore tenda a precedere di qualche giorno la convergenza del flusso di quantità di moto è coerente con l’idea classica del ciclo di vita delle onde barocline: crescita alimentata dalla baroclinicità e successiva fase in cui i flussi di quantità di moto rimodellano il getto e la struttura del vento medio. Studi sul life cycle baroclino e sul ruolo dei flussi di quantità di moto mostrano infatti come la fase “barotropa” e i processi di breaking possano diventare dominanti dopo la crescita iniziale. 

La coerenza tra modello completo e a settore vale anche per un dettaglio non banale: le oscillazioni quasi periodiche su ~10 giorni (tra lag −40 e 0) nei campi di flusso, e soprattutto l’inversione della sequenza ai tropici/subtropici (dove la convergenza del flusso di quantità di moto può anticipare il flusso di calore). Questo pattern “a due regimi latitudinali” è un indizio che la propagazione non è un semplice segnale locale, ma emerge dalla struttura latitudinale del feedback eddy–mean flow e dalla diversa natura dei transitori (baroclini extratropicali vs contributi più barotropi/di riflessione/riorganizzazione del getto a latitudini inferiori), un tema che ricorre anche nelle interpretazioni moderne della propagazione delle anomalie zonali e dei modi annulari propaganti. 

Il punto forse più “forte” della comparazione è però un altro: la somiglianza tra i campi compositi dei due modelli suggerisce che la localizzazione zonale delle onde e dei pacchetti d’onda (cioè il fatto che nel modello completo si possano vedere wave packets più o meno confinati longitudinalmente) non è essenziale per la propagazione verso il polo. In altre parole, la propagazione sembra essere governata principalmente da quantità zonali-medie e da un meccanismo che sopravvive anche quando si impone una struttura zonale più “rigida” (settoriale). Questo tipo di evidenza è coerente con l’idea che la propagazione sia una proprietà emergente del sistema eddy-driven, sostenuta da rottura d’onda, mixing di PV e migrazione delle latitudini critiche, più che da dettagli di pacchetti localizzati. 

A fronte delle forti somiglianze, emergono due differenze dinamicamente istruttive. La prima riguarda la scala temporale: nel modello a settore il periodo della propagazione risulta più corto (circa 125 giorni) rispetto alla banda ~120–160 giorni tipica del modello completo e di esperimenti affini, ma resta comunque molto più lungo del periodo osservato di circa 60 giorni, già discusso in ambito osservativo sin dagli studi storici e poi quantificato in modo più sistematico nelle reanalisi. Il fatto che un modello ideale produca periodi più lunghi non è sorprendente: semplificazioni fisiche (forzanti rilassate, assenza di ciclo stagionale, mancanza di eterogeneità superficie–orografia–contrasti termici reali) possono alterare la forza dei feedback, la persistenza del getto e i tempi di riassetto della baroclinicità. Ma la discrepanza con l’osservato resta un “vincolo” importante, perché implica che il modello cattura il meccanismo, ma non necessariamente i tempi caratteristici del sistema reale. 

La seconda differenza è spaziale: nel settore la propagazione tende a terminare attorno a 50°, mentre nel modello completo prosegue a latitudini più alte. Qui il motivo proposto è molto fisico e legato alla geometria sferica e alla scelta spettrale: se si fissa il numero d’onda fondamentale a 5, alle alte latitudini la lunghezza d’onda fisica associata a quel numero d’onda si riduce perché la circonferenza zonale decresce con cosφ. In un sistema reale (o in un modello con spettro più ricco), una struttura che “sale” di latitudine può efficacemente traslare verso numeri d’onda più bassi per conservare scale orizzontali comparabili; nel settore, invece, le onde più lunghe (k < 5) non esistono, e questo può distorcere l’interazione onda–flusso medio proprio dove la propagazione dovrebbe continuare. In termini di meccanismi, si tratta di un limite che può modificare breaking, mixing di PV e la posizione delle latitudini critiche, “disinnescando” il feedback necessario a sostenere l’avanzamento verso il polo. 

Nel complesso, il confronto full-spectrum vs settore fornisce quindi un messaggio molto pulito: la propagazione verso il polo è un fenomeno che sembra poggiare su pochi processi cardine (feedback tra flusso zonale ed eddies, sequenze coerenti tra flussi di calore e di quantità di moto, e riorganizzazione della struttura del getto), tanto da sopravvivere a una forte semplificazione della variabilità zonale. Al tempo stesso, la rappresentazione spettrale impone vincoli reali su periodo e raggiungimento delle alte latitudini: ed è proprio lì che si vede quanto la “gerarchia di modelli” sia utile, perché separa ciò che è essenziale per il meccanismo da ciò che controlla tempi e dettagli spaziali, che nel mondo osservato restano più rapidi e più estesi in latitudine. 

Evoluzione temporale della propagazione verso il polo nel modello a settore: PV di Ertel, diagnostica EP e risposta della circolazione media

Una volta verificato che il modello a settore riesce a riprodurre la propagazione verso il polo del modello completo, l’analisi dell’evoluzione temporale del campo di moto si basa su una lettura “coerente” di tre pezzi che si parlano tra loro: la sequenza composita del potenziale vorticoso di Ertel (PV) a 250 hPa, la comparsa e la migrazione di una banda di anomalia negativa del vento zonale medio alle basse latitudini, e i segnali diagnostici associati ai flussi di Eliassen–Palm (EP) e alla circolazione meridionale media (MMC). L’uso del PV come traccia principale non è solo una scelta grafica: nella dinamica quasi-adiabatica e a grande scala, le mappe di PV sono uno strumento molto potente per riconoscere steepening, overturning e mixing legati al breaking, perché rendono visibile l’erosione del gradiente meridiano e la riorganizzazione delle strutture d’onda. 

Nel composito, la finestra temporale considerata copre circa un ciclo: dall’avvio della banda negativa alle bassissime latitudini tropicali fino al “riavvio” del ciclo con una nuova banda che ricompare vicino all’equatore. In parallelo, la presenza (a margine dei pannelli) della componente divergente del vento meridionale mediata zonalmente consente un confronto diretto tra quanto velocemente la banda migra in latitudine e quanto rapidamente scorre il flusso meridionale divergente di fondo.

All’inizio del ciclo, quando la banda negativa comincia a formarsi equatorward di circa 10°, il PV mostra solo un accenno di ripidimento dell’onda: è la fase in cui l’anomalia negativa del vento zonale medio “si aggancia” ai primi segnali di deformazione dell’onda. In questa fase i segnali nella diagnostica EP risultano deboli ma organizzati: compare una convergenza alle basse latitudini dominata dalla componente orizzontale del flusso EP, un comportamento che si allinea bene con la lettura classica in termini di dinamica di latitudine critica, cioè con l’idea che l’interazione onda–flusso medio tenda a intensificarsi dove la velocità di fase dell’onda e il vento medio favoriscono assorbimento e deposito di quantità di moto. 

Proseguendo verso metà ciclo, la firma del wave breaking diventa molto più netta nel campo di PV: le creste d’onda mostrano una marcata inclinazione nord-est/sud-ovest, una geometria che, nella lettura standard delle onde di Rossby, è indicativa di una connessione con sorgenti e propagazione provenienti dalle medie latitudini e con trasporti di quantità di moto coerenti con quella struttura inclinata. Nello stesso momento, l’orientamento dei vettori EP è compatibile con una propagazione dell’attività d’onda che “arriva” dalle medie latitudini e poi interagisce in modo più intenso alle basse latitudini. 

Durante la fase di overturning (il ribaltamento dell’onda), il centro della banda negativa avanza verso il polo in modo relativamente rapido; segue poi una fase più lunga in cui domina il mescolamento del PV, cioè la parte irreversibile del ciclo in cui le strutture di PV vengono rimescolate e i gradienti si attenuano. Questo quadro è pienamente coerente con il paradigma del “surf zone” e, più in generale, con la descrizione del breaking come processo che trasforma un disturbo d’onda organizzato in una regione di mixing capace di ridistribuire PV e, tramite i flussi turbolenti, di forzare il flusso medio. 

Un aspetto quantitativo interessante della sequenza composita è che la velocità di migrazione della banda, stimata in modo semplice dalle sue posizioni a lag successivi, risulta molto vicina alla velocità del flusso meridionale divergente poleward mostrato nei pannelli. Questo aggancio tra “deriva della banda” e flusso divergente è reso dinamicamente interpretabile guardando la MMC: la cella anomala alle basse latitudini appare legata a un forcing degli eddies, nel senso che la convergenza/divergenza dei flussi diagnostici (riassunta dalla divergenza del flusso EP) può sostenere una circolazione anomala termicamente diretta in regione tropicale-subtropicale. In termini TEM, è proprio qui che la diagnostica EP diventa preziosa: traduce l’azione collettiva degli eddies sulla circolazione zonalmente media in un segnale che si può seguire nel tempo e in quota. 

La stessa sequenza si ripete con un secondo episodio di breaking: nei primi giorni l’avanzamento latitudinale del centro della banda è più rapido, seguito poi da una deriva più lenta verso il polo con velocità di nuovo molto simile a quella del flusso divergente ambientale. In questa fase, nel flusso divergente totale compaiono due massimi positivi: quello più vicino alla latitudine della banda coincide bene con la struttura della MMC anomala; quello più equatoriale è ricondotto alla componente diabatico-riscaldante del modello, un punto che torna spesso quando si interpretano segnali tropicali in sistemi idealizzati dove la rappresentazione del riscaldamento può imprimere una firma netta sulla circolazione meridionale. 

Verso la fine del secondo episodio, l’attività d’onda torna a rinvigorirsi: nella diagnostica EP questo si manifesta con un rafforzamento della componente verticale del trasporto di attività d’onda, e nel PV con una nuova intensificazione della deformazione fino a un breaking molto marcato al lag 0, che nel composito appare come un evento di ristrutturazione ampia e intensa del campo di PV a 250 hPa. 

Nel terzo episodio, rispetto a quanto visto nei compositi al lag −32 giorni (primo episodio) e al lag −18 giorni (secondo episodio), la regione di wave breaking mostra un’estensione meridionale molto più ampia: è un segnale dinamicamente coerente con un breaking più fortemente non lineare, perché l’overturning del campo di PV tende a coinvolgere una fascia latitudinale più larga quando la distorsione delle superfici a PV quasi-conservato e il mescolamento diventano più intensi. L’interpretazione di questi passaggi tramite mappe di PV è in linea con la letteratura “classica” sull’uso diagnostico del PV isentropico per riconoscere intrusione, filamentazione, overturning e mixing associati a breaking e riorganizzazioni del flusso medio. 

Dopo questo terzo evento segue una fase relativamente quieta, dal lag 0 a circa il lag +50 giorni, in cui l’attività d’onda che raggiunge le basse latitudini risulta attenuata (coerentemente con l’indebolimento dei segnali nei flussi turbolenti e con vettori EP anomali prevalentemente orientati verso il basso lungo la colonna). In un contesto simile, un ruolo crescente viene assunto dai processi di rilassamento radiativo del modello, che possono ricostruire progressivamente i venti occidentali nella troposfera tropicale superiore quando la forzante ondosa dalle medie latitudini non “disturba” in modo efficace la regione tropicale. La rapida scomparsa della struttura ondulata interna alla banda tropicale del vento zonale medio anomalo tra il lag 0 e il lag +8 giorni, e la sua persistenza fino a circa il lag +33 giorni, è compatibile con un regime in cui la componente ondosa perde organizzazione e la parte zonalmente media torna a essere controllata soprattutto da bilanci lenti (radiativi e di circolazione media), più che da deposizione episodica di quantità di moto da parte degli eddies.

In questa fase quieta, la struttura zonalmente media evidenzia un massimo positivo del vento zonale medio anomalo centrato intorno a 15°, collocato sul bordo polare di una regione relativamente ampia associata alla componente divergente del vento meridionale. Questa configurazione richiama il legame, ben noto nella teoria delle circolazioni assialsimmetriche, tra ramo discendente della cella di Hadley e formazione/mantenimento del getto subtropicale in prossimità del bordo della circolazione tropicale: un accoppiamento che, nei modelli ideali, emerge come conseguenza della redistribuzione di momento angolare e dei vincoli di equilibrio termico-dinamico. 

Un tassello diagnostico aggiuntivo viene dalla lettura congiunta della latitudine critica e del gradiente meridiano medio di PV a 250 hPa, sovrapposti alle anomalie del vento zonale medio. Nei tre episodi di wave breaking, la latitudine critica tende a sovrapporsi sia alle anomalie negative del vento zonale medio sia alle regioni con gradiente meridiano di PV ridotto, e la sua posizione coincide con le aree in cui nei campi compositi si osservano overturning e breaking. Questo comportamento è coerente con la teoria delle latitudini critiche e dei critical layers, dove l’avvicinamento a condizioni critiche favorisce assorbimento, forte interazione onda–flusso medio e transizione a regime non lineare; in tali regioni la propagazione “lineare” delle onde si interrompe e la risposta del flusso medio può essere forzata in modo marcato dalla convergenza/divergenza dei flussi associati alle onde. 

La stessa figura diagnostica mostra anche che una banda positiva del vento zonale medio anomalo tende a comparire già attorno al lag −20 giorni sul lato equatoriale della regione a gradiente meridiano di PV debole. Se quella fascia a gradiente ridotto agisce come “barriera” alla propagazione delle onde di Rossby (nel senso che ostacola l’invasione di attività d’onda dalle medie latitudini verso i Tropici profondi), allora sul lato equatoriale della barriera i processi di rilassamento radiativo possono dominare più facilmente, perché la circolazione locale risulta relativamente schermata dall’impatto dei flussi d’onda extratropicali. L’idea che il breaking vada inteso come un processo generale di perdita di regolarità della propagazione ondosa e di transizione verso mixing e deposizione di quantità di moto è coerente con le definizioni estese di wave breaking proposte per onde di Rossby e di gravità, e con il quadro in cui i critical layers possono comportarsi come regioni di assorbimento (o, in certi regimi, di riflessione/over-reflection). 

Il ruolo del rilassamento radiativo risulta supportato dalla coerenza tra tempi: durante il periodo quieto la latitudine critica arretra verso l’equatore fino a valori intorno a 10° in circa 30 giorni (dal lag −20 al lag +10), lo stesso ordine di grandezza della scala temporale di rilassamento radiativo impostata nel modello, proprio mentre si consolidano le anomalie occidentali tropicali. In seguito, attorno al lag +30 giorni, quando la barriera legata al gradiente di PV debole inizia a cedere, il gradiente meridiano di PV alle medie latitudini torna ad aumentare e ricompare una nuova banda negativa alle basse latitudini, predisponendo il sistema alla ripartenza del ciclo di propagazione verso il polo, con un ritorno a episodi di attività d’onda extratropicale e breaking il cui posizionamento latitudinale tende a rimanere strettamente connesso alla migrazione della latitudine critica. 

Test di sensibilità: risoluzione orizzontale, diffusione sottogriglia e ruolo relativo tra assorbimento alle latitudini critiche e forcing turbolento

Il quadro interpretativo delineato in precedenza attribuisce la propagazione verso il polo a una sequenza ricorrente di Rossby wave breaking e di apparente assorbimento dell’attività d’onda in prossimità delle latitudini critiche, cioè nelle regioni in cui la propagazione ondosa tende a interrompersi e il momento (e più in generale l’attività d’onda) viene trasferito al flusso medio. Questa lettura è coerente con la teoria dei critical layers, che mostra come tali regioni possano comportarsi, a seconda del regime e dei gradienti di vorticità potenziale, come zone di assorbimento, ma anche (in circostanze diverse) di riflessione o addirittura di over-reflection; in altre parole, l’esito non è “garantito” a priori e dipende in modo sensibile dalla non linearità e dai dettagli dissipativi. 

Proprio per questo, è legittimo chiedersi se in un modello con risoluzione orizzontale relativamente grossolana e diffusione sottogriglia intensa l’assorbimento non risulti eccessivo: la diffusione numerica (o parametrizzata) può infatti aumentare l’efficienza con cui l’attività d’onda viene “smorzata” nel critical layer, limitando la possibilità di riflessione e rendendo più netta la firma di assorbimento. Nella letteratura sui critical layers di onde di Rossby, lavori storici come quelli di Stewartson (e affini) hanno mostrato come la viscosità/“viscosità efficace” condizioni l’evoluzione del critical layer e quindi l’equilibrio tra assorbimento e riflessione; in regimi meno dissipativi, la riflessione diventa una possibilità concreta e, nel quadro proposto dagli autori, ciò dovrebbe ostacolare la propagazione verso il polo. 

I test numerici riportati mirano quindi a stressare proprio questi ingredienti “numerici”: viene aumentata la risoluzione (configurazione a risoluzione più alta rispetto al controllo) e si varia in modo sostanziale il coefficiente di diffusione orizzontale. Il risultato è che la propagazione verso il polo rimane una caratteristica evidente nelle diverse configurazioni, suggerendo una robustezza del fenomeno rispetto a variazioni ragionevoli di risoluzione e dissipazione. Questo tipo di messaggio è rilevante perché, in generale, la climatologia e l’intensità degli eddies possono rispondere sensibilmente allo schema e al livello di diffusione: già in studi “classici” sulla sensibilità dei GCM, modificare la diffusione orizzontale può alterare energia cinetica turbolenta e conversioni energetiche, con ricadute sul vento zonale medio. Allo stesso tempo, lavori più moderni sulla progettazione della dissipazione (ad esempio iperdiffusione con migliori proprietà di conservazione energetica) mostrano che una parte della sensibilità alla dissipazione può essere mitigata se lo smorzamento sottogriglia è formulato in modo più fisicamente coerente, migliorando anche la convergenza con la risoluzione. In questo contesto, il fatto che la propagazione persista al variare della diffusione va letto come indicazione che il “cuore” del meccanismo non è un artefatto stretto di un singolo settaggio dissipativo, pur restando vero che i dettagli del breaking e del deposito di quantità di moto possono dipendere dalla formulazione numerica.

Un secondo blocco di sensibilità riguarda i parametri che, nei modelli idealizzati, rappresentano i “tempi lenti” della fisica: il coefficiente di drag superficiale (Rayleigh drag nel boundary layer) e la scala temporale del rilassamento radiativo (Newtonian cooling verso un profilo di equilibrio). Questi ingredienti sono centrali nei benchmark di circolazione generale secca, dove l’obiettivo è ottenere un’atmosfera baroclinicamente attiva con forzanti minime ma controllabili; in tale cornice, drag e rilassamento fissano vincoli fondamentali, ma non determinano da soli la risposta del vento zonale medio, che spesso emerge dall’equilibrio tra forzante termica e trasporti turbolenti. 

Nei test descritti, raddoppiare o dimezzare il drag, così come accorciare o allungare la scala di rilassamento radiativo (ad esempio da valori dell’ordine di due settimane a valori dell’ordine di due mesi), non produce cambiamenti apprezzabili nella velocità di propagazione. Questo punto è dinamicamente informativo: indica che la propagazione (e, in particolare, il suo ritmo) non è controllata direttamente dal tempo imposto del rilassamento radiativo, e quindi non è semplicemente un “orologio termico” del modello. Una lettura coerente, già implicita nella diagnostica EP/TEM, è che l’aggiustamento del flusso zonale medio sia dominato soprattutto dai flussi turbolenti e dalla circolazione meridionale media indotta dagli eddies, cioè da quella parte della risposta media che nasce dall’interazione onda–flusso medio più che dal ritorno radiativo all’equilibrio. Il linguaggio operativo che lega in modo formale flussi d’onda, forzante del vento zonale medio e circolazione residua è proprio quello della teoria TEM e delle sezioni EP, sviluppata e resa strumentale per l’interpretazione dei dati da Andrews & McIntyre e da Edmon, Hoskins & McIntyre. 

La figura 7 mostra, per il modello a settore, l’evoluzione composita del potenziale vorticoso (PV) a 250 hPa in una serie di lag (ad es. −60, −32, −30, −20, −18, −12 giorni). Il PV a questa quota è un diagnostico particolarmente “parlante” perché rende visibile la geometria delle onde di Rossby e, soprattutto, i passaggi verso la non linearità: steepening, overturning e mixing (mescolamento) emergono come piegamenti e arrotolamenti delle strutture di PV, cioè come perdita di monotonia latitudinale e sviluppo di filamenti e “lingue” che segnalano breaking e trasporto irreversibile di PV. Questa lettura è coerente con la tradizione del “PV-thinking” e con l’uso delle mappe isentropiche/isentropico-equivalenti per interpretare i processi dinamici a grande scala. 

Le ombreggiature (scuro/chiaro/assenza) non sono una classificazione fisica “rigida” ma un espediente grafico: le soglie di PV cambiano da colonna a colonna (come indicato in didascalia) proprio per mettere in risalto le strutture ondose che gli autori vogliono seguire nel tempo. In altre parole, l’informazione chiave non è tanto “quanto PV” in assoluto, ma come si deforma il campo: quando compaiono intrusioni meridiane di PV e quando le strutture diventano sufficientemente inclinate e arrotolate da indicare overturning e mixing. Un’idea utile qui è quella di “surf zone”, introdotta per descrivere regioni dove il PV viene stirato e rimescolato in modo irreversibile durante il breaking delle onde di Rossby. 

Sopra il PV sono tracciati i contorni spessi delle anomalie del vento zonale a 250 hPa: linee continue per valori positivi, tratteggiate per valori negativi (la linea di zero è omessa). Questo overlay serve a collegare visivamente le fasi di breaking/mixing con la nascita e l’evoluzione della banda di anomalia negativa del vento zonale medio (il segnale che, nel lavoro, viene poi seguito nella sua propagazione verso il polo). A lato di ogni pannello compare il profilo zonalmente medio: la linea spessa è l’anomalia zonal-mean del vento zonale, mentre la linea sottile è il vento meridionale divergente medio (moltiplicato per 30 solo per renderlo confrontabile sulla stessa scala grafica). Questo margine destro è cruciale perché consente di confrontare la deriva latitudinale delle anomalie con la ventilazione meridiana divergente che, nel quadro interpretativo del paper, contribuisce a trasportare lentamente la banda dopo i “salti” rapidi associati ai singoli breaking.

Al lag −60, il campo di PV suggerisce un regime ancora relativamente “lineare”: la struttura ondosa è presente ma non mostra ancora un overturning marcato alle basse latitudini. In questa fase iniziale, la banda negativa nelle anomalie di vento è solo abbozzata in prossimità della fascia tropicale più bassa; dinamicamente, è il tipo di configurazione in cui ci si aspetta un primo segnale di steepening (ripidimento) senza ancora un mixing esteso del PV. Nel linguaggio del wave–mean flow interaction, è lo stadio in cui l’attività d’onda può iniziare a convergere verso una regione favorevole all’interazione, ma senza aver ancora “scaricato” in modo evidente il suo effetto sul flusso medio.

Passando al lag −32 e soprattutto al lag −30, la deformazione del PV diventa più organizzata e compaiono segnali più chiari di overturning alle basse latitudini: la geometria delle creste e la loro inclinazione indicano un trasporto meridiano associato alle onde, e la co-localizzazione con la banda negativa delle anomalie di vento segnala che la regione sta entrando in un regime di interazione più non lineare. Questo è precisamente il tipo di firma che, nella letteratura sul Rossby wave breaking, viene associato alla fase terminale del ciclo di vita di un’onda: da una perturbazione che propaga e cresce a una struttura che si ripiega e rimescola PV, con effetti persistenti sul vento zonale medio. 

I pannelli attorno a lag −20 e lag −18 evidenziano un secondo episodio, in cui le strutture di PV alle basse latitudini appaiono più vigorose e più estese: si riconosce meglio la logica “intrusiva” del breaking, con tongue/streamer che indicano scambi meridiani più forti e un mixing più efficace. In questa fase, la banda negativa delle anomalie di vento si consolida e tende a collocarsi a latitudini un po’ più alte rispetto alla fase iniziale, coerentemente con l’idea che ogni episodio di breaking produca un riassetto del flusso medio e predisponga la successiva traslazione verso il polo.

Il pannello a lag −12 mostra uno stadio più avanzato del ciclo: rispetto ai lag precedenti, l’impressione è di un campo di PV più “rimaneggiato”, cioè con una traccia più netta del mixing precedente e con strutture ondose che prefigurano un ulteriore rinvigorimento (nel paper, questo viene poi collegato ai burst successivi e alla preparazione di breaking più intenso). La sequenza complessiva è compatibile con un’interpretazione in termini di critical-layer dynamics: il breaking tende a concentrarsi dove le condizioni favoriscono l’assorbimento/trasformazione dell’attività d’onda e dove il sistema passa dalla propagazione quasi lineare alla non linearità con mixing di PV. È importante ricordare che, in generale, i critical layers possono comportarsi come regioni di assorbimento ma, in regimi meno dissipativi o con gradienti di PV particolari, anche la riflessione (o l’over-reflection) può diventare rilevante: questo è uno dei motivi per cui gli autori discutono poi test di sensibilità su risoluzione e diffusione. 

Sebbene la figura 7 non mostri direttamente i flussi di Eliassen–Palm, la lettura “PV + anomalie di vento” è perfettamente coerente con la diagnostica TEM: quando il breaking e il mixing diventano evidenti, ci si aspetta un concomitante trasferimento di quantità di moto dalle onde al flusso medio, formalizzabile come convergenza/divergenza dei flussi d’onda e risposta della circolazione residua. In questo senso, la figura 7 è la controparte planimetrica (a 250 hPa) di ciò che le EP cross-sections rappresentano in sezione: dove l’onda si organizza, dove si deforma, e dove lascia un’impronta persistente sul vento zonale medio. 

Infine, il margine destro suggerisce anche un collegamento “di equilibrio” tra la struttura del vento zonale tropicale/subtropicale e la circolazione meridionale divergente: la collocazione dei massimi del vento zonale medio anomalo in prossimità del bordo della ventilazione divergente richiama il legame tra cella di Hadley e getto subtropicale nel quadro delle circolazioni assialsimmetriche, in cui il trasporto di momento angolare e la dinamica del ramo discendente contribuiscono a fissare la fascia di massimi occidentali subtropicali. Questo riferimento è utile perché, nella fase in cui l’attività d’onda si attenua, i processi “lenti” di aggiustamento radiativo e di circolazione media possono riemergere nel controllo della struttura del vento tropicale superiore. 

La figura 8 sintetizza, in un’unica diagnostica, tre aspetti che in un’analisi “onda–flusso medio” sono inseparabili: (i) dove e in che direzione viene trasportata l’attività d’onda (vettori del flusso di Eliassen–Palm anomalo), (ii) dove tale attività viene effettivamente trasferita al flusso medio (divergenza del flusso EP, rappresentata dai contorni sottili e dall’ombreggiatura per i valori negativi), e (iii) quale risposta produce la circolazione media meridiana(contorni della funzione di corrente di massa anomala, cioè la MMC anomala). Questo impianto interpretativo è quello del formalismo TEM e della diagnostica EP resa operativa in ambito atmosferico da lavori fondamentali come Andrews e McIntyre (1976, 1978) ed Edmon, Hoskins e McIntyre (1980), dove il flusso EP fornisce una misura compatta della propagazione dell’attività d’onda e la sua convergenza/divergenza quantifica la forzante netta degli eddies sul vento zonale medio e sulla circolazione residua; il quadro è poi divenuto standard anche nella lettura troposferica/stratosferica dei processi di forcing ondoso (Holton et al., 1995).

In ogni pannello la coordinata verticale è la pressione (quindi si leggono le strutture dall’alta troposfera verso il basso), mentre l’asse orizzontale è la latitudine. I vettori mostrano il flusso EP anomalo: quando prevale una componente meridiana, la figura suggerisce un trasporto anomalo dell’attività d’onda lungo latitudine; quando la componente verticale è marcata, indica un’anomalia di trasferimento verticale dell’attività d’onda, spesso legata a fasi di intensificazione della propagazione e a scambi più efficaci tra livelli. I contorni sottili mostrano la divergenza del flusso EP: l’ombreggiatura evidenzia i valori negativi (cioè la convergenza del flusso EP anomalo). Nel linguaggio TEM, una convergenza del flusso EP corrisponde a un torque ondoso sul flusso zonale medio: è il segnale che l’attività d’onda non sta solo “passando”, ma sta venendo assorbita/trasformata in una regione, producendo un cambiamento persistente del vento zonale medio. Questo è concettualmente lo stesso tipo di firma che, nella teoria dei critical layers e nelle interpretazioni del breaking di onde di Rossby, accompagna la transizione da propagazione quasi lineare a regime non lineare con mixing di vorticità potenziale e deposito di quantità di moto (McIntyre e Palmer, 1983; Killworth e McIntyre, 1985).

I contorni mediamente spessi (con la linea di zero enfatizzata) rappresentano la MMC anomala: in pratica descrivono le celle anomale di circolazione nel piano meridiano. Qui la figura fa vedere un punto chiave: la circolazione media non è un “ornamento” diagnostico, ma parte della risposta forzata dagli eddies. In un quadro TEM, quando la convergenza/divergenza del flusso EP si organizza su scala substagionale, la circolazione meridiana media tende a riorganizzarsi in modo coerente, perché deve bilanciare la forzante ondosa e i vincoli di equilibrio dinamico. Questa coerenza è esattamente ciò che permette di collegare i segnali in planimetria (PV e bande di vento zonale) alla dinamica in sezione: dove le onde depositano quantità di moto, lì si struttura una risposta della MMC, e quella risposta contribuisce a trasportare e mantenere le anomalie del vento zonale medio.

Guardando la sequenza dei lag, il pannello a lag −60 mostra un segnale ancora relativamente debole: i vettori EP anomali sono meno organizzati e le aree ombreggiate (convergenza) sono limitate. Questa fase corrisponde a un contesto in cui l’interazione onda–flusso medio sta iniziando a “precondizionare” le basse latitudini, ma il deposito netto di attività d’onda è ancora modesto. Nonostante ciò, la MMC anomala è già riconoscibile, il che è coerente con l’idea che anche segnali ondosi piccoli ma persistenti possano avviare una risposta media che poi si amplifica quando l’attività d’onda aumenta (un tema comune negli studi sulla forzante eddy-driven della circolazione media).

Nei pannelli a lag −32 e −30, che corrispondono alla prima fase di overturning/breaking evidenziata nelle mappe di PV, la figura mostra un salto di qualità: aumentano sia l’organizzazione dei vettori EP sia l’estensione delle regioni ombreggiate, indicando che il trasporto anomalo dell’attività d’onda non solo è presente, ma tende a convergere in una fascia latitudinale ben definita. Questo è il tipo di segnale che ci si aspetta quando onde provenienti dalle medie latitudini entrano in un regime in cui l’assorbimento (o la trasformazione non lineare) diventa efficiente, con conseguente trasferimento di quantità di moto al flusso medio. La MMC anomala risponde in modo più netto, con celle più marcate: questa co-variazione tra convergenza EP e intensificazione della MMC è una firma tipica del coupling eddy–mean flow e, nella pratica, è ciò che rende plausibile l’interpretazione secondo cui una parte del flusso meridionale divergente “ambientale” che trasporta lentamente la banda di vento zonale è in realtà eddy-driven (Andrews e McIntyre, 1976; Edmon et al., 1980).

Nei pannelli a lag −20 e −18 (secondo episodio), la struttura appare ancora più robusta e coerente: le aree di convergenza del flusso EP sono più evidenti e la MMC anomala risulta più energica e ben separata in celle. Questa fase, nel linguaggio del lavoro, coincide con un nuovo burst di attività d’onda e un nuovo episodio di breaking: la figura 8 mostra come quell’episodio abbia una controparte dinamica “integrata”, cioè un incremento del trasporto/propagazione dell’attività d’onda e, soprattutto, un incremento della sua conversione in forzante sul flusso medio (convergenza). In termini più generali, è la stessa logica che emerge nei cicli di vita delle onde barocline e nelle classificazioni del breaking (ad esempio gli archetipi LC1/LC2), dove l’evoluzione non lineare determina pattern diversi di trasporto e deposito di quantità di moto, con conseguenze differenti sul vento zonale medio (Thorncroft et al., 1993).

Il pannello a lag −12 mostra un assetto in cui la forzante ondosa anomala rimane organizzata e la MMC continua a rispondere, suggerendo una fase di transizione verso l’ulteriore rinvigorimento dell’attività d’onda descritto nel testo (preludio al breaking successivo). Anche senza entrare nei dettagli numerici, l’informazione fisica importante è che l’intero ciclo non è governato solo da “onde che si vedono nel PV”, ma da una catena coerente: incremento dell’attività d’onda, trasporto anomalo rappresentato dal flusso EP, convergenza come impronta del deposito/assorbimento (spesso in prossimità di latitudini critiche), e risposta della circolazione media meridiana che contribuisce a riorganizzare il vento zonale medio. È proprio questa coerenza tra vettori EP, loro convergenza e MMC che rende la figura 8 il ponte diagnostico tra la fenomenologia planimetrica del breaking e la dinamica media che sostiene la propagazione verso il polo, nel solco della letteratura su interazione onda–flusso medio e critical-layer dynamics (Edmon et al., 1980; Andrews e McIntyre, 1978; Killworth e McIntyre, 1985; Holton et al., 1995).

La figura 9 è pensata per rendere “visibile” il legame tra tre ingredienti che, nella dinamica atmosferica a grande scala, governano l’interazione onda–flusso medio: la propagazione dell’attività d’onda (rappresentata dai vettori del flusso di Eliassen–Palm totale), la struttura del vento zonale medio (mostrata dai contorni) e le latitudini critiche (evidenziate con l’ombreggiatura chiara). In questo tipo di diagnostica, diventata standard a partire dai lavori fondativi su Trasformed Eulerian Mean e flusso EP (Andrews e McIntyre; Edmon, Hoskins e McIntyre), le frecce non indicano “vento”, ma la direzione lungo cui l’attività d’onda tende a propagare e a trasferire la sua influenza sul flusso medio. La presenza del vento zonale medio sullo stesso pannello è essenziale perché la propagazione delle onde di Rossby non è libera: dipende fortemente dal profilo dei venti occidentali e dalla presenza di regioni con venti orientali, che possono agire come barriere dinamiche o deviare i percorsi di propagazione, un concetto coerente anche con la teoria classica della propagazione verticale delle onde planetarie (Charney e Drazin).

I contorni mostrano quindi dove i venti occidentali sono più intensi e come cambiano con latitudine e quota; la campitura scura mette in evidenza le aree con venti orientali, che in genere risultano poco favorevoli alla propagazione “classica” delle onde di Rossby verso l’alto o verso latitudini più basse. La campitura chiara, invece, identifica i settori in cui il vento zonale medio è molto vicino alla velocità di fase dell’onda fondamentale: sono le cosiddette latitudini critiche. Dal punto di vista teorico, le latitudini critiche sono il luogo in cui la descrizione lineare della propagazione tende a rompersi e dove diventano probabili assorbimento, forte interazione non lineare e deposito di quantità di moto sul flusso medio; la letteratura sui critical layers delle onde di Rossby (in particolare Killworth e McIntyre) sottolinea che lì si gioca l’equilibrio tra assorbimento e, in alcuni regimi meno dissipativi, possibile riflessione o over-reflection. In altre parole: la figura non sta solo “marcando un punto”, ma sta evidenziando la regione dove la teoria si aspetta che l’onda smetta di comportarsi come una perturbazione che attraversa il dominio e inizi invece a modificare attivamente la circolazione media.

La sequenza temporale dei pannelli (lag −60, −32, −30, −20, −18, −12) mostra che l’attività d’onda tende a organizzarsi e a dirigersi ripetutamente verso la fascia subtropicale/tropicale superiore, e che questa direzione di propagazione intercetta spesso la regione identificata come critica. Ai lag più precoci (−60) la zona critica è già presente e le frecce suggeriscono un trasporto di attività d’onda che la raggiunge; nei lag associati agli episodi più attivi (−32/−30 e poi −20/−18) la coerenza fra “dove puntano” le frecce e “dove sta” la latitudine critica diventa più evidente. Questo è il tipo di firma che ci si aspetta quando burst di attività d’onda provenienti dalle medie latitudini vengono guidati in una fascia dove la propagazione si arresta e l’interazione con il flusso medio diventa efficiente: è esattamente il ponte diagnostico che collega la fenomenologia del breaking/mixing osservabile nelle mappe di PV (figura 7) con la forzante ondosa e la risposta media (figura 8), ma qui espresso in termini di vincolo imposto dal vento medio e dalla latitudine critica.

Un punto importante che la figura 9 aiuta a interiorizzare è che la latitudine critica non è un “luogo fisso”: si muove nel tempo perché si muove il profilo del vento zonale medio. Basta che la fascia di venti occidentali si rafforzi, si indebolisca o trasli in latitudine perché anche la regione in cui il vento medio e la velocità di fase dell’onda si equivalgono cambi posizione. Nel quadro del lavoro, questa mobilità è fondamentale perché rende plausibile che anche la sede preferenziale di assorbimento/breaking e, di conseguenza, la banda di anomalia del vento zonale che si organizza nel ciclo, possano migrare con un certo ordine temporale. In sintesi, la figura 9 non è “una figura in più”: è la prova geometrico-dinamica che il percorso dell’attività d’onda è guidato dalla struttura dei venti medi e che le regioni critiche – dove la teoria prevede massima interazione onda–flusso medio – sono effettivamente collocate lungo quel percorso e variano con esso.

La figura 10 condensa in un solo diagramma tempo–latitudine ciò che, nelle figure precedenti, era distribuito tra mappe planari (PV e breaking) e sezioni meridiane (flusso EP): qui l’obiettivo è mostrare come l’evoluzione del vento zonale medio in alta troposfera, la posizione delle latitudini critiche e l’emergere di una fascia a gradiente meridiano di vorticità potenziale quasi-geostrofica molto debole siano parte dello stesso ciclo fisico. Sull’asse orizzontale c’è il lag (giorni), sull’asse verticale la latitudine; i contorni rappresentano l’anomalia del vento zonale a 250 hPa (con un passo fine, e con alcune isolinee aggiuntive per evidenziare segnali deboli), mentre i puntini quadrati tracciano la latitudine critica definita in figura 9, cioè la fascia in cui il vento medio “visto” dall’onda si avvicina alla sua velocità di fase. In questa cornice, la latitudine critica non è un dettaglio grafico: è il luogo in cui la propagazione ondosa di tipo rossbiano tende a interrompersi e dove l’interazione onda–flusso medio può diventare massimamente efficace, favorendo assorbimento e processi non lineari. Questa idea discende direttamente dalla formulazione TEM e dall’uso diagnostico dei flussi di Eliassen–Palm, che legano propagazione dell’attività d’onda e forcing sullo stato zonalmente medio. 

L’ombreggiatura, invece, evidenzia le regioni in cui il gradiente meridiano della vorticità potenziale quasi-geostrofica è sotto una soglia (quindi “insolitamente debole”). La ragione fisica per cui questo è importante è che le onde di Rossby, per propagare e mantenere una struttura coerente, hanno bisogno di un ambiente con gradienti di PV sufficientemente marcati; quando il gradiente si appiattisce, la propagazione diventa inefficiente e la regione può comportarsi come una barriera dinamica, riducendo l’invasione dell’attività d’onda verso le basse latitudini. In termini più generali, questa dipendenza dalla struttura del flusso medio e dei suoi gradienti è coerente con la teoria classica della propagazione delle onde planetarie, che mostra come i “corridoi” di propagazione siano controllati dallo sfondo dei venti zonali e dalle proprietà di rifrazione del mezzo. 

Il punto chiave della figura 10 è la co-localizzazione sistematica, durante gli episodi attivi del ciclo, tra (i) latitudine critica (puntini), (ii) anomalia negativa del vento zonale e (iii) fascia a gradiente di PV debole. Questa sovrapposizione è esattamente ciò che ci si aspetta se l’evoluzione sia governata da dinamica di critical layer: quando l’attività d’onda proveniente da latitudini più alte raggiunge una regione prossima alla latitudine critica, la descrizione lineare perde validità, l’onda tende a essere assorbita o trasformata, e diventa probabile l’innesco di breaking con mescolamento del PV e con un’impronta persistente sul vento medio. La teoria sottolinea anche che l’esito (assorbimento puro oppure riflessione/over-reflection) dipende dalla dissipazione efficace e dalla struttura locale dei gradienti di PV; il fatto che, nel composito, i puntini critici “stiano dentro” la fascia di gradiente debole mentre si organizzano le anomalie negative del vento è quindi una firma qualitativa molto coerente con un regime dominato da assorbimento e mixing. 

La figura 10 aggiunge inoltre un’informazione temporale che nelle mappe non è immediata: la latitudine critica migra nel tempo perché dipende dal profilo del vento zonale medio. Quando l’anomalia del vento si intensifica o si sposta, anche la fascia in cui il vento medio “eguaglia” la velocità di fase dell’onda trasla, e con essa si sposta la regione preferenziale di interazione. Questa mobilità fornisce un vincolo geometrico-dinamico al ciclo: non è l’onda a scegliere liberamente dove rompersi, ma è il flusso medio che, evolvendo, “sposta” il punto in cui la propagazione diventa critica. La coerenza tra questa migrazione e la sequenza degli episodi attivi è precisamente il tipo di evidenza che ha reso i diagrammi EP/TEM uno strumento così potente per collegare fenomenologia (breaking, mixing) e bilanci medi (forzante eddy e risposta del vento zonale). 

Un ulteriore elemento diagnostico è che la fascia ombreggiata (gradiente di PV debole) può essere interpretata come una “zona schermante” che tende a ridurre l’arrivo di flussi d’onda dalle medie latitudini verso i Tropici profondi. Questo concetto è in linea con la visione del breaking come processo che appiattisce gradienti e genera regioni di mixing (in stratosfera è stato storicamente descritto come “surf zone”, ma la logica di fondo—erosione del gradiente e aumento della miscelazione—è trasferibile anche a contesti troposferici idealizzati). In tale configurazione, quando la barriera è ben formata, diventa più facile che i processi lenti del modello (in particolare il rilassamento radiativo del vento medio verso un equilibrio) ricostruiscano i westerlies tropicali superiori, proprio perché la forzante ondosa in ingresso risulta attenuata; quando invece la barriera si indebolisce e i gradienti recuperano, le onde possono nuovamente penetrare più in profondità e il ciclo riparte con la formazione di una nuova anomalia negativa alle basse latitudini. 

La figura 11 è uno schema concettuale che distilla, in forma “meccanicistica”, ciò che nelle figure diagnostiche precedenti viene mostrato con campi e compositi: la propagazione verso il polo non è un moto continuo e uniforme della banda di vento zonale medio anomalo, ma un ciclo a impulsi in cui l’attività d’onda extratropicale, l’interazione con la latitudine critica e la risposta della circolazione media si agganciano tra loro in modo ricorrente. Nel pannello superiore, la curva di EKE (eddy kinetic energy) riassume proprio questa natura intermittente: la fase attiva è scandita da più massimi successivi, interpretabili come “burst” di attività ondosa alle medie latitudini, seguiti da una fase più quieta in cui l’attività d’onda si attenua nettamente. Questa idea di ciclicità a episodi, con crescita e decadimento della turbolenza sinottica e della conversione energetica, è coerente con la letteratura sui cicli di vita delle onde barocline e sulle transizioni tra regimi di propagazione e breaking, in cui i flussi turbolenti diventano temporaneamente molto efficienti nel forzare il flusso medio (ad esempio, la famiglia di lavori su life cycle e breaking come Thorncroft, Hoskins e McIntyre 1993).

La riga centrale sintetizza, a quattro tempi ideali, l’evoluzione zonal-mean del vento zonale e della circolazione media meridiana associata. Le linee orizzontali marcate come c.l. rappresentano le latitudini critiche, cioè quelle fasce in cui il vento medio “visto” dall’onda si avvicina alla velocità di fase: è il punto in cui la propagazione lineare tende a perdere validità e dove la teoria dei critical layers prevede assorbimento e forte interazione onda–flusso medio, con elevata probabilità di processi non lineari e deposizione netta di quantità di moto sul flusso medio (Killworth e McIntyre 1985; e, più in generale, il quadro TEM/EP sviluppato da Andrews e McIntyre 1976, 1978 e reso operativo in diagnostica atmosferica da Edmon, Hoskins e McIntyre 1980). I simboli “+” e “−” non sono decorativi: indicano, rispettivamente, convergenza e divergenza dei flussi di quantità di moto turbolenti, cioè i luoghi in cui gli eddies accelerano o decelerano localmente il vento zonale medio. In questo schema il segnale centrale è la comparsa di una banda anomala negativa di vento zonale alle basse latitudini in coincidenza con l’interazione ondosa vicino alla latitudine critica: l’interpretazione è che le onde, originate alle medie latitudini, propagano verso sud fino a incontrare la regione critica, dove vengono assorbite/rompono, trasferendo quantità di moto e rimodellando il profilo del vento. È precisamente il tipo di meccanismo che, nel linguaggio del PV-thinking, corrisponde a breaking e mixing del PV con formazione di una “zona di surf” e appiattimento dei gradienti, concetti introdotti e discussi ampiamente per l’atmosfera da McIntyre e Palmer (1983, 1984).

Il passaggio successivo, mostrato nei pannelli centrali intermedi, è altrettanto importante: la propagazione verso il polo del segnale zonal-mean non è attribuita solo al momento dell’overturning/breaking (fase “rapida”), ma anche a una deriva più lenta associata alla circolazione meridiana media (indicata nello schema come poleward MMC). In termini TEM, quando i flussi d’onda convergono/divergono, la circolazione residua tende a rispondere in modo coerente, contribuendo a redistribuire momento e a trasportare le anomalie; questo accoppiamento tra forcing turbolento e risposta della circolazione media è un punto cardine del formalismo di Andrews e McIntyre e della diagnostica EP. La figura 11 enfatizza infatti che la MMC non è un semplice “contorno” diagnostico, ma una parte attiva della catena causale: una componente del flusso meridionale divergente che accompagna la banda anomala può essere interpretata come eddy-driven, cioè sostenuta dall’azione netta dei flussi turbolenti (Edmon et al. 1980).

Un altro ingrediente chiave, rappresentato esplicitamente nello schema con la dicitura “small” o “moderate” associata al gradiente meridiano di PV, è la formazione di una regione a gradiente debole che agisce da barriera alla propagazione delle onde. Quando i gradienti di PV si appiattiscono per effetto del mixing, l’ambiente diventa meno favorevole alla propagazione rossbiana verso le basse latitudini: in modo qualitativo, è come se venisse meno la “guida” dinamica che consente alle onde di mantenere una propagazione coerente. Questa idea si raccorda al principio, ben noto già dalla teoria della propagazione delle onde planetarie, secondo cui la struttura del vento medio e dei gradienti associati controlla corridoi e barriere alla propagazione (Charney e Drazin 1961). Nello schema, tale barriera è ciò che inaugura la fase quieta: riducendo l’invasione di attività d’onda verso i Tropici profondi, consente ai processi lenti del modello (in particolare il rilassamento radiativo del vento zonale medio verso un profilo di equilibrio) di ricostruire progressivamente westerlies tropicali in alta troposfera. Il collegamento concettuale con la dinamica assialsimmetrica tropicale è implicito: quando l’interferenza ondosa extratropicale è attenuata, la circolazione tropicale tende a riemergere come attore dominante nella strutturazione dei venti subtropicali, in linea con il quadro classico del legame tra cella di Hadley e getto subtropicale (Held e Hou 1980; Schneider 1977).

La riga inferiore, infine, traduce la stessa sequenza in termini di struttura orizzontale dell’onda: a ciascuno dei quattro tempi della riga centrale corrisponde un diverso assetto delle creste e dei flussi, e soprattutto un diverso posizionamento rispetto alla latitudine critica. Il messaggio è che il wave packet/onda, quando raggiunge la fascia critica, entra in una fase di forte deformazione e breaking che produce un pattern coerente di convergenza e divergenza dei flussi di quantità di moto; quel pattern è ciò che imprime la risposta zonal-mean (banda anomala) e ne innesca lo spostamento latitudinale. In questo senso, la figura 11 funziona da “ponte” tra fenomenologia (burst di EKE, breaking, mixing) e bilanci medi (forzante turbolenta e risposta della MMC), collocando il tutto nel quadro teorico consolidato dell’interazione onda–flusso medio e dei critical layers (Andrews e McIntyre 1976, 1978; Edmon et al. 1980; Killworth e McIntyre 1985; McIntyre e Palmer 1983, 1984).

7. Conclusioni

Come discusso nella sezione 6b, la propagazione verso il polo delle anomalie del flusso zonale medio nel modello a settore risulta riconducibile alla cooperazione tra tre ingredienti dinamici: la propagazione (in prima approssimazione lineare) delle onde di Rossby generate alle medie latitudini, la loro transizione non lineare verso il wave breaking con conseguente riorganizzazione del flusso medio, e i processi di rilassamento radiativo che ripristinano gradualmente lo stato di fondo, soprattutto nella fascia tropicale e subtropicale superiore. L’eventuale obiezione metodologica è che, per costruzione, un modello a settore possa favorire un certo tipo di rottura d’onda e quindi amplificare artificialmente il peso del wave breaking rispetto a un modello a spettro completo; tuttavia, la forte somiglianza tra i due esperimenti (settore e completo) rende poco persuasiva questa lettura riduttiva, perché entrambi mostrano un’organizzazione temporale molto simile dell’attività ondosa e della risposta del flusso medio, inclusa la ricorrenza di tre episodi distinti di crescita ondosa e breaking. Una coerenza di questo tipo è tipicamente interpretata come evidenza che il meccanismo fondamentale non dipenda da dettagli “accidentali” della configurazione, ma da un accoppiamento robusto tra pacchetti d’onda, deposito di quantità di moto e aggiustamento della circolazione media. 

Nel quadro teorico dell’interazione onda–flusso medio, la presenza di breaking di tipo “cat’s eye” (anche se osservato su campi compositi) è particolarmente informativa perché richiama la dinamica dei critical layers: regioni in cui le onde, avvicinandosi a una latitudine in cui il flusso medio è “in risonanza” con la velocità di fase, tendono a perdere la capacità di propagazione quasi-lineare e a trasferire in modo efficiente quantità di moto e attività d’onda al flusso medio, favorendo mixing e irreversibilità. La letteratura sui Rossby-wave critical layers sottolinea inoltre che l’esito locale può includere assorbimento e, in altri regimi, anche riflessione o over-reflection, con un ruolo sensibile della dissipazione efficace e dei gradienti di vorticità potenziale; ciò rende ancora più significativo il fatto che, nelle diagnostiche del lavoro, breaking e massimi cambiamenti del vento zonale medio si organizzino proprio attorno alla regione critica. 

Poiché il flusso di quantità di moto associato alle onde di Rossby appare il cardine della catena causale, questa interpretazione è in continuità con l’idea di James e Dodd, che attribuiva la propagazione verso il polo a onde generate alle medie latitudini e alla successiva interazione con il flusso zonale medio. Nello stesso tempo, il meccanismo evidenziato qui si discosta in punti cruciali: nel modello a settore emergono segnali forti che la dinamica di latitudine critica non sia un dettaglio marginale, ma un perno che organizza dove avvengono il deposito di quantità di moto e i principali cambiamenti del vento zonale medio; inoltre, nelle fasce tropicali e subtropicali, la traslazione verso il polo non sembra spiegabile soltanto come risposta locale ai flussi di quantità di moto ondosi, ma risulta coerente anche con un contributo di avvezione da parte del flusso meridionale divergente associato alla circolazione media (un accoppiamento ben rappresentato nella diagnostica TEM/EP). Infine, sempre in Tropici e subtropici, la ricostruzione delle anomalie positive del vento zonale durante le fasi “quiet” viene ricondotta in modo preferenziale al rilassamento radiativo del modello, più che a una convergenza persistente del flusso di quantità di moto ondoso: questo punto è concettualmente coerente con il fatto che, quando l’invasione ondosa è schermata o indebolita, i processi lenti di aggiustamento termodinamico e di circolazione media riacquistano peso nel definire la struttura del vento subtropicale superiore. 

Resta aperta, sul piano logico, la possibilità che un contributo “alla James–Dodd” coesista con il quadro basato sul wave breaking, ma le evidenze diagnostiche spingono a collocare l’azione più intensa proprio nelle vicinanze della latitudine critica, dove breaking e massimi cambiamenti del vento zonale medio tendono a sovrapporsi; in un caso del genere, la distinzione operativa tra i due processi diventerebbe intrinsecamente difficile perché entrambi si manifesterebbero nello stesso spazio-fase della circolazione. Inoltre, poiché la formulazione di James e Dodd nasce in un contesto debolmente non lineare e non incorpora esplicitamente il breaking, un’eventuale componente di quel tipo risulterebbe più naturalmente interpretabile come contributo secondario rispetto ai processi non lineari di rottura e mixing, che qui compaiono come elemento organizzatore della variabilità e del deposito di quantità di moto sul flusso medio. Nel quadro dell’interazione onde–flusso zonalmente medio, è utile mettere a fuoco l’alternativa proposta da Robinson: la propagazione verso il polo non sarebbe primariamente “trascinata” dal breaking e dall’assorbimento in prossimità della latitudine critica, ma dalla deriva poleward della zona baroclina, realizzata da una circolazione meridiana media termicamente indiretta forzata dalla convergenza/divergenza del flusso di quantità di moto degli eddies. In questa famiglia di idee, lo spostamento della baroclinicità (e quindi della storm track) diventa il motore del segnale zonal-mean: la circolazione media, indotta dalla forzante di quantità di moto, ristruttura il gradiente termico e fa migrare il massimo di instabilità, che a sua volta riorganizza la generazione d’onda. Il punto interessante, però, è che nel caso che stai traducendo questo “quadro Robinson” non emerge in modo pulito: la regione di forte propagazione verticale dell’attività d’onda (la firma classica di rigenerazione baroclina e sviluppo sinottico) mostra sì una traslazione poleward, ma tende a collocarsi sul bordo poleward della cella termicamente indiretta invece che nel suo cuore, e questo rende meno naturale leggere la deriva della baroclinicità come causa prima della propagazione. 

Una chiave interpretativa alternativa è collegare quel comportamento a ciò che spesso si osserva negli studi sui baroclinic life cycles multipli: dopo un primo ciclo, l’evoluzione non “si spegne” semplicemente, ma può favorire un secondo ciclo in cui la zona baroclina e la massima attività d’onda risultano traslate di diversi gradi verso il polo, come esito del modo in cui il primo ciclo ha precondizionato PV, shear e struttura del getto. In questa prospettiva, lo spostamento poleward del massimo flusso di calore turbolento e della rigenerazione baroclina può essere un effettodella propagazione (e del riassetto zonal-mean indotto dagli eddies), più che la sua causa: le onde in quota già presenti possono indurre perturbazioni al suolo e riattivare la crescita baroclina “da sopra”, in un modo concettualmente affine alla ciclogenesi di tipo B (sviluppo fortemente condizionato dal forcing in quota e dall’accoppiamento con la struttura frontale al basso livello). Questa lettura è coerente con l’idea generale che il breaking e la riorganizzazione del PV non siano solo un “finale” del ciclo sinottico, ma possano modulare la variabilità interna della circolazione troposferica e predisporre nuovi episodi di crescita/rottura in posizioni latitudinali diverse. 

Sul piano diagnostico, il fatto che sulle scale temporali lunghe della propagazione (intraseasonali) la convergenza del flusso di quantità di moto emerga come segnale guida è in linea con la lettura TEM/EP: il deposito di quantità di moto e la risposta della circolazione media sono il canale privilegiato con cui l’attività d’onda “scrive” un’anomalia persistente nel vento zonale medio. La sfumatura importante che il testo mette in evidenza è quella di scala: su finestre più brevi, durante la sequenza degli episodi, compaiono prima flussi di calore coerenti e solo dopo flussi di quantità di moto ben organizzati, suggerendo una sorta di handover tra crescita baroclina (che proietta subito sul calore) e maturazione/rottura (che proietta più nettamente sulla quantità di moto). Questo è perfettamente compatibile con la fenomenologia dei life cycles (LC1/LC2 e varianti) in cui la cronologia tra conversioni energetiche, tilt delle onde e breaking controlla quando e dove si massimizza il forcing sul flusso medio. 

Resta poi il nodo della velocità di propagazione, che nelle simulazioni risulta troppo bassa rispetto alle stime osservative citate dagli autori: anche avendo un quadro meccanicistico robusto, capire cosa “setta” la velocità richiede un controllo più diretto sui parametri effettivi della riorganizzazione zonal-mean. Un’ipotesi fisicamente sensata è che il modello generi onde che rompono con una scala meridiana troppo piccola: se il mixing e la modifica del flusso medio avvengono su un intervallo latitudinale ristretto, anche lo spostamento poleward della regione critica e della zona di interazione tende a essere limitato per episodio, e la propagazione composta su molti episodi risulta lenta. Un secondo controllo plausibile è l’ampiezza degli eddies: onde più energetiche mescolano più rapidamente e possono traslare più velocemente le anomalie del vento zonale medio (e quindi la geometria della successiva interazione), un’idea che si raccorda sia alla teoria dei critical layers sia alla visione “PV-mixing” della rottura d’onda. Proprio per isolare questi controlli, la proposta di un esperimento futuro con un modello idealizzato (per esempio barotropico) forzato da un wave maker disegnato ad hoc è metodologicamente centrata: permette di variare in modo pulito scala e ampiezza dell’onda, e misurare come cambiano mixing, spostamento della regione critica e velocità di migrazione del segnale zonal-mean. 

La figura 12 è un test di robustezza “da dinamical core” che serve a verificare se la propagazione verso il polo delle anomalie del vento zonale medio in alta troposfera sia un prodotto fragile della dissipazione numerica, oppure un comportamento dinamico che sopravvive a variazioni ragionevoli di risoluzione e diffusione. In ciascun pannello viene mostrato un segmento di circa 2000 giorni della serie temporale del vento zonalmente medio a 250 hPa, rappresentato in un diagramma tempo–latitudine: l’inclinazione delle bande (aree ombreggiate per anomalie positive, aree non ombreggiate per anomalie negative) è la firma visiva della migrazione meridiana delle anomalie. Le strutture diagonali che “risalgono” in latitudine con il tempo corrispondono proprio alla propagazione poleward discussa nel lavoro; se la propagazione fosse un artefatto, ci si aspetterebbe che cambi drasticamente (o scompaia) al variare della dissipazione.

Il confronto è organizzato in due dimensioni. La prima riguarda il tipo di modello: nei pannelli (a) e (b) è mostrato il modello completo, nei pannelli (c) e (d) il modello a settore. La seconda riguarda l’intensità della diffusione orizzontale: (a) e (c) utilizzano la configurazione di controllo, mentre (b) e (d) riducono la diffusione (indicata come “weak diffusion”). Questo è un punto metodologicamente delicato perché, nei modelli spettrali e nei modelli idealizzati, la diffusione (spesso iperdiffusione) non è un dettaglio marginale: controlla quanto rapidamente l’enstrofia e le piccole scale vengono smorzate, e quindi può alterare l’equilibrio tra propagazione, assorbimento e possibile riflessione delle onde. Proprio per questo, la sensibilità alla dissipazione è un controllo standard negli esperimenti di interconfronto dei “dynamical cores”, come sottolineato nella proposta di benchmark Held–Suarez, che ha reso esplicito quanto le statistiche della circolazione possano dipendere dalle scelte numeriche e dai termini dissipativi, pur mantenendo forzanti idealizzate. 

Il messaggio principale della figura 12 è che la propagazione poleward rimane chiaramente presente in tutti e quattro i casi. Cambiano alcuni dettagli (regolarità delle bande, intermittente intensificazione locale, grado di “rumorosità” alle alte latitudini), ma la geometria fondamentale—bande coerenti che migrano verso latitudini più alte su scale substagionali—persiste sia nel modello completo sia nel modello a settore, e soprattutto persiste quando si indebolisce la diffusione. Questo risultato è importante perché riduce la probabilità che la propagazione dipenda in modo critico da un eccesso di smorzamento numerico: se fosse la diffusione a “fabbricare” l’assorbimento ondoso necessario alla migrazione, ci si aspetterebbe che riducendola emergano comportamenti qualitativamente diversi (ad esempio maggiore riflessione, maggiore stazionarietà o disorganizzazione del segnale). Invece, la propagazione resta un tratto robusto della dinamica simulata.

Il test è anche fisicamente motivato da ciò che la teoria e gli studi sui critical layers suggeriscono: la transizione tra assorbimento e riflessione delle onde di Rossby può dipendere in modo sensibile da viscosità/dissipazione e dalla struttura dei gradienti di vorticità potenziale, perché la regione critica è un luogo in cui la descrizione lineare si rompe e i processi non lineari (mixing, filamentazione, trasferimento verso armoniche) diventano centrali. Per questo, l’idea che una diffusione più debole possa aumentare la riflessione e ostacolare la propagazione è perfettamente plausibile in astratto; lavori classici sulla dinamica dei critical layers discutono proprio come la presenza o assenza di dissipazione efficace condizioni l’evoluzione della regione critica e la risposta a lungo termine. La figura 12 mostra che, almeno nel range esplorato, la propagazione non collassa quando si riduce lo smorzamento, il che è coerente con un meccanismo dominato da interazione onda–flusso medio e mixing non lineare che non dipende “solo” da un parametro diffusivo.

Un ulteriore valore della figura 12 è che avvicina questo studio alla discussione moderna sulla dissipazione sub-grid nei modelli atmosferici: non conta soltanto “quanta” dissipazione si applica, ma anche “come” essa interagisce con il bilancio energetico e con la cascata verso le piccole scale. Studi più recenti hanno mostrato che schemi dissipativi energeticamente più coerenti possono modificare in modo non banale la qualità delle simulazioni a bassa risoluzione e il comportamento delle strutture eddy-driven, proprio perché la dissipazione controlla quanto rapidamente l’energia viene rimossa e come ciò retroagisce sulla macroturbolenza. In questo contesto, il fatto che la propagazione poleward compaia sia in controllo sia in weak diffusion, e in entrambi i “dynamical cores” (full e sector), suggerisce che il segnale sia legato a un’organizzazione a grande scala del sistema (burst di attività ondosa, deposito di quantità di moto, risposta del flusso medio) più che a un delicato equilibrio numerico.

In sintesi, la figura 12 funziona come una prova empirica che la propagazione verso il polo delle anomalie del vento zonale medio a 250 hPa è un comportamento strutturale del modello: compare nel modello completo e nel modello a settore, e resta evidente anche quando si riduce la diffusione orizzontale. Questo non elimina ogni sensibilità quantitativa (periodo, ampiezza e regolarità possono variare), ma rafforza l’interpretazione che il fenomeno sia sostenuto da processi dinamici interni—interazione onda–flusso medio e transizioni verso il breaking—piuttosto che da un artefatto legato a una specifica scelta di dissipazione numerica.

8. Discussione e implicazioni

Sebbene Riehl et al. (1950) non abbiano formulato un meccanismo fisico esplicito per la propagazione verso il polo, il loro lavoro contiene diversi indizi osservativi che, riletti con l’impianto interpretativo di questo studio, acquisiscono una coerenza dinamica notevole. In particolare, l’associazione tra le anomalie del vento zonale e configurazioni riconducibili a wave breaking è compatibile con l’idea che la propagazione delle anomalie zonal-mean sia legata a episodi di crescita d’onda, interazione con il flusso medio e successiva riorganizzazione non lineare del campo di PV, come suggerito dalla moderna “PV-thinking” e dalla letteratura sul breaking rossbiano (McIntyre e Palmer, 1983, 1984). Il dettaglio riportato da Riehl et al. (1950)—il massimo dell’anomalia occidentale che coincide con il bordo meridionale di una saccatura in quota e con una configurazione locale di breaking anticiclonico—è esattamente il tipo di firma che ci si aspetta quando la dinamica si concentra in prossimità di una regione critica e produce un deposito netto di quantità di moto sul flusso medio. Ancora più significativo è l’altro elemento descritto nello studio del 1950: quando il massimo dell’anomalia occidentale attraversa la latitudine del getto medio, l’ampiezza delle onde alle medie latitudini risulta piccola. Questa relazione è coerente con un quadro in cui, nelle fasi di attraversamento del getto climatologico, l’attività turbolenta misurabile tramite flusso di calore ed energia cinetica turbolenta tende a un minimo, come spesso avviene quando il sistema si colloca in un “intervallo di transizione” tra episodi di crescita/rottura e periodi più quieti; un comportamento di questo tipo è anche in linea con la visione TEM/EP dell’interazione onda–flusso medio, in cui la struttura e la variabilità del vento zonale medio riflettono l’alternanza tra fasi di forcing ondoso efficiente e fasi di ridotta penetrazione dell’attività d’onda (Andrews e McIntyre, 1976, 1978; Edmon, Hoskins e McIntyre, 1980).

Un passaggio concettualmente importante riguarda la possibile connessione con il “zonal index cycle” classico di Rossby (1939) e Namias (1950). Il rallentamento marcato della propagazione alle latitudini superiori a circa 40°—comportamento riscontrato anche in analisi osservative moderne (Feldstein, 1998)—suggerisce che, oltre una certa soglia latitudinale, la variabilità del vento zonale medio possa essere descritta in modo comparabile come una transizione tra stati a indice zonale alto e basso, più che come un’unica anomalia che trasla in modo regolare. In termini contemporanei, questa descrizione richiama l’idea di “regimi” o stati preferenziali della circolazione extratropicale, in cui la dinamica alterna configurazioni relativamente persistenti e transizioni rapide, con implicazioni dirette sulla prevedibilità substagionale. Il punto proposto dagli autori—che l’impressionante propagazione descritta da Riehl et al. e l’apparente successo delle previsioni basate sull’indice zonale di Namias possano essere due facce dello stesso fenomeno, anche perché i periodi di dati si sovrappongono—è plausibile: se una componente quasi periodica o quasi regolare del vento zonale medio emerge in certe finestre temporali, allora sia l’approccio “indice zonale” sia la lettura in termini di propagazione poleward possono catturare aspetti reali della stessa dinamica, con la differenza che il secondo quadro tenta di legare tale regolarità a un ciclo fisico di crescita d’onda, breaking e riaggiustamento radiativo.

L’implicazione più ampia è legata alla prevedibilità e a una possibile gerarchia di strutture prevedibili in un sistema atmosferico non lineare. In un limite idealizzato lineare, soluzioni ondose perfettamente periodiche renderebbero la previsione del “tempo” concettualmente banale; nella pratica, la non linearità distrugge la periodicità delle singole onde e accelera la perdita di memoria iniziale. Tuttavia, la non linearità può anche favorire l’emergere di nuove strutture coerenti su scale spazio-temporali maggiori, come i pacchetti d’onda e la loro propagazione lungo le storm tracks, concetto formalizzato in modo generale dalla teoria delle onde modulanti e dei pacchetti (Whitham, 1974) e documentato in atmosfera da lavori classici sui treni d’onda e sulla propagazione dell’attività sinottica (Namias e Clapp, 1944; Trenberth, 1986; Chang, 1993). Anche la modellistica idealizzata ha mostrato che i pacchetti d’onda possono essere riprodotti in modo robusto e che la loro coerenza dipende dal regime dinamico e dai parametri che controllano la dispersione e la non linearità (Lee e Held, 1993). In questo contesto, l’osservazione che i pacchetti in propagazione verso est risultino meno coerenti proprio nelle regioni dello spazio dei parametri in cui la propagazione poleward del vento zonale medio è più evidente suggerisce un trade-off dinamico: quando aumentano gli spostamenti meridiani e il mixing (quindi l’irreversibilità associata al breaking), si perde parte della coerenza dei pacchetti, ma può emergere una nuova “firma” prevedibile a scala più lenta, cioè la migrazione organizzata delle anomalie zonal-mean. Ne deriva un’ipotesi di lavoro stimolante: aumentando la non linearità, alcune forme di prevedibilità “fine” (legate alle singole onde) decadono, mentre si rafforzano forme di prevedibilità “grossolana” legate all’attività d’onda nel suo insieme (pacchetti) e, su scale ancora più lente, alle anomalie del flusso medio che propagano verso il polo. Questa prospettiva, pur preliminare, è metodologicamente utile perché sposta l’attenzione dalla sola previsione deterministica degli stati sinottici alla previsione di proprietà organizzate della circolazione (attività d’onda, regimi, anomalie zonal-mean), cioè esattamente il territorio concettuale in cui si collocano molte idee moderne sulla prevedibilità substagionale e sui limiti dinamici imposti dall’interazione onda–flusso medio.

https://t.co/BK92kii58b