Pan-basin warming now overshadows robust Pacific Decadal Oscillation

Received: 5 August 2024

Allison A. Cluett

1

, Steven J. Bograd

Mercedes Pozo Buil1 & Elliott L. Hazen

2,3, Michael G. Jacox

1,3,4

,

1,3

Accepted: 8 October 2025

Published online: 7 November 2025

Check for updates

Per quasi trent’anni la Pacific Decadal Oscillation (PDO) è stata trattata come una “variabile ponte” tra fisica del clima e risposta ecosistemica nel Pacifico settentrionale: un indice sintetico capace di descrivere, con un solo numero, una parte rilevante della variabilità decadale delle SST e dei suoi effetti a cascata su circolazione atmosferica, struttura della colonna d’acqua, produttività e risorse biologiche (il caso-scuola è quello delle salmon fisheries). L’impianto concettuale nasce con i lavori fondativi di Mantua et al. (1997), che legano i cambiamenti interdecadali nel Nord Pacifico a segnali coerenti negli stock salmonicoli e in più variabili fisiche e biologiche. 

Il punto chiave del lavoro di Cluett et al. (2025) è però che, nel clima attuale, la PDO non è “sparita” né si è indebolita nel senso statistico classico: ciò che è cambiato è il contesto di fondo (baseline) su cui quella variabilità si innesta. Gli autori mostrano che la configurazione spaziale della PDO e la sua forza relativa (quando la PDO è definita nel modo canonico, cioè come EOF principale delle anomalie di SST nel Nord Pacifico dopo rimozione del segnale di temperatura media globale) restano sostanzialmente stabili; tuttavia, se si guarda alle anomalie totali senza rimuovere il riscaldamento di fondo, emerge ormai un pattern dominante di segno uniforme su gran parte del bacino, che loro chiamano pan-basin pattern (PBP). In altre parole: la variabilità “interna” tipo PDO continua a esistere, ma viene sempre più sovrastata (in ampiezza e rilevanza pratica) dal segnale di riscaldamento pan-bacino. 

Il risultato è elegante anche dal punto di vista metodologico: usando finestre mobili di 30 anni, gli autori evidenziano che fino a circa il 2013 la rimozione della global-mean SST prima di calcolare la PDO è quasi “banale” (PDO calcolata con o senza rimozione del segnale globale è quasi identica), mentre includendo gli anni successivi il primo EOF delle SST totali tende a diventare unidirezionale (tutti i loadings dello stesso segno) e spiega una quota di varianza che cresce fino a valori molto elevati rispetto al passato (ordine di grandezza: dal ~20% fino a oltre il 50% nelle finestre più recenti, secondo quanto riportato nell’articolo). Contestualmente, la PDO “scende di rango” diventando un EOF successivo nelle SST totali, non perché si indebolisca fisicamente, ma perché la varianza del riscaldamento pan-bacino diventa dominante. 

Questa impostazione chiarisce anche un equivoco frequente: quando diciamo “PDO” possiamo riferirci a (i) un indice costruito per isolare la variabilità interna (detrending e/o rimozione della media globale) oppure a (ii) un PDO index operativo calcolato su dataset osservativi che, per costruzione, può inglobare parte della tendenza forzata se non si separano correttamente i contributi. Le pagine operative NOAA mostrano infatti che l’indice PDO (versioni NCEI/PSL) è derivato da ERSST e da una proiezione su una mappa di regressione, con scelte metodologiche precise che vanno sempre esplicitate quando si fanno attribuzioni o inferenze ecologiche. 

Da qui discende la parte più “applicativa” (e, se vuoi, più destabilizzante) del paper: molte correlazioni PDO–ecosistema non sono crollate perché la PDO abbia cambiato faccia, ma perché gli organismi e le comunità rispondono spesso a valori assoluti (o a soglie fisiologiche) più che a semplici anomalie standardizzate. Cluett et al. formalizzano questa idea con un modello concettuale in cui la risposta biologica è centrata su un optimum termico: a parità di oscillazione PDO, se la baseline si sposta verso temperature più alte per effetto del PBP, lo stesso “valore di PDO” può portarti prima verso l’optimum (decorrelazione), poi oltre l’optimum (inversione del segno della correlazione). È un modo molto pulito per spiegare perché una relazione storicamente positiva possa diventare nulla o addirittura negativa nell’ultimo decennio, senza invocare per forza una “mutazione” della PDO. 

Questo quadro si incastra bene con quanto sappiamo su altri modi interni del Nord Pacifico che continuano a essere fisicamente significativi ma che, nella metrica delle SST totali, possono perdere centralità apparente. Gli autori ricordano che l’EOF associato a NPGO/Victoria Mode resta riconoscibile e coerente nel tempo: ed è importante perché, da anni, la letteratura sottolinea come la NPGO catturi aspetti dinamici della variabilità del giro subtropicale e subpolare e della circolazione costiera che non sono riducibili alla sola PDO, con impatti distinti sulla biogeochimica e sugli ecosistemi del California Current e dell’Alaska Gyre. 

Il messaggio diventa ancora più concreto se lo connettiamo al tema delle marine heatwaves e degli eventi “no-analogue”. La definizione quantitativa oggi più usata per le marine heatwaves (durata, intensità, estensione, soglie percentile rispetto a una climatologia) deriva dal framework di Hobday et al. (2016), mentre l’episodio del 2013–2014 nel NE Pacific (“the Blob”) è uno degli esempi più studiati di anomalia persistente con impatti ecosistemici e atmosferici, descritto in modo classico da Bond et al. (2015). In un Nord Pacifico che si riscalda “di bacino”, la probabilità che oscillazioni interne (PDO/NPGO/ENSO-teleconnessioni) si sovrappongano a una baseline già elevata e producano estremi cresce, e questo è esattamente il tipo di regime non stazionario che Cluett et al. mettono al centro. 

Infine, sullo sfondo c’è il tema – più ampio e ancora molto dibattuto – del rapporto tra variabilità interna e forzante antropica nelle oscillazioni decadali del Pacifico. Il lavoro di Newman et al. (2016) rimane un riferimento perché chiarisce che la PDO non è un “unico meccanismo”, ma un indice emergente da più processi (teleconnessioni tropicali, accoppiamento oceano–atmosfera nel Nord Pacifico, rumore atmosferico integrato dall’oceano, ecc.). In parallelo, studi più recenti sostengono che una parte importante delle variazioni multidecadali dell’indice PDO osservato possa essere forzata (aerosol e gas serra) e che i modelli possano sottostimare l’ampiezza della variabilità forzata: è una linea interpretativa che, pur non coincidente con l’impostazione “PDO interna + PBP come baseline”, aiuta a capire perché il confine tra “modo interno” e “segnale di fondo” sia oggi un problema pratico, non solo teorico. 

In sintesi, Cluett et al. propongono un cambio di prospettiva molto operativo: la PDO resta un descrittore utile della dinamica interna, ma nel Nord Pacifico contemporaneo l’informazione che serve davvero per interpretare impatti e prevedibilità è spesso la sovrapposizione tra (i) oscillazioni interne e (ii) riscaldamento pan-bacino che rende la baseline non stazionaria. E questo è coerente con la cornice generale del riscaldamento oceanico e dei cambiamenti criosferici e di circolazione descritti nelle valutazioni IPCC: la “nuova normalità” è un sistema in cui trend e variabilità non si sommano più come componenti separabili in modo ingenuo, ma interagiscono nel determinare estremi e risposte ecosistemiche. 

Nel Pacifico settentrionale la variabilità climatica a bassa frequenza è stata per decenni interpretata attraverso la Pacific Decadal Oscillation (PDO), cioè il principale modo empirico della variabilità delle anomalie di temperatura superficiale del mare (SST) a nord di ~20°N: nella sua fase positiva tende a combinare SST relativamente più calde nel Pacifico orientale e lungo la costa nordamericana con SST relativamente più fredde nel Pacifico centrale-occidentale, accompagnate da riorganizzazioni coerenti della circolazione atmosferica e oceanica. Nella pratica operativa, però, “PDO” non è solo un concetto fisico ma anche una scelta metodologica: ad esempio, nella definizione NOAA/PSL la PDO è calcolata come EOF principale delle SSTA nel Nord Pacifico dopo rimozione della climatologia mensile e della media globale di SST, proprio per separare quel pattern dalla componente di riscaldamento globale; esistono poi indici affini mantenuti da altri centri (NCEI) con costruzioni leggermente diverse ma fortemente correlate. Questa distinzione è cruciale, perché la PDO è un modo statistico emergente (non un “ingranaggio unico”), che integra contributi da forzanti atmosferiche, teleconnessioni tropicali (ENSO) e dinamica extratropicale, con l’oceano che “arrossa” (reddening) il segnale atmosferico su scale interannuali-decadali. 

L’importanza storica della PDO è legata anche alla narrazione dei regime shifts: salti di segno o riorganizzazioni rapide del sistema (spesso datate attorno a 1925, 1947 e 1976–77) che hanno coinciso con cambiamenti marcati in molte variabili fisiche e biologiche. In questo quadro, il lavoro fondativo di Mantua e colleghi ha reso “gestibile” un problema complesso: collegare un pattern di SST su scala di bacino a risposte ecologiche, in particolare alle variazioni nella produzione e nelle catture di salmonidi nel Nord Pacifico. La transizione del 1976–77, spesso trattata come evento canonico, è stata associata a risposte opposte tra Alaska/Gulf of Alaska e costa occidentale USA, con l’idea (poi molto sfruttata) di “produzioni inverse” e di grandi riconfigurazioni nella distribuzione/abbondanza delle specie con impatti socioeconomici. Questa cornice ha alimentato l’uso della PDO come indicatore sintetico in contesti gestionali e previsionali, ma anche il rischio di interpretarla come causa unica, quando in realtà è un indice che riassume una parte della variabilità del sistema. 

Proprio perché la PDO è un indicatore composito, la letteratura ha da tempo sottolineato che nel Nord Pacifico convivono più “assi” di variabilità decadale con firme dinamiche differenti. Un esempio è la North Pacific Gyre Oscillation (NPGO), introdotta per catturare variazioni legate alla dinamica dei gyres e a segnali biogeochimici (salinità, nutrienti, clorofilla) che non sempre si allineano alla PDO: in molte applicazioni, PDO e NPGO forniscono informazioni complementari su struttura di circolazione e impatti ecosistemici. In parallelo, lavori di sintesi hanno chiarito che gran parte della variabilità decadale pacifica può essere letta come risultato dell’interazione tra teleconnessioni tropicali-extratropicali, risposta del Kuroshio-Oyashio Extension e integrazione dell’oceano del forcing atmosferico; quindi l’“indice giusto” dipende dalla domanda (fisica, ecologica, gestionale) e dalla metrica (anomalie relative vs valori assoluti). Anche la scelta tecnica di dataset e climatologia non è neutra: differenze tra prodotti SST e finestre di riferimento possono introdurre discrepanze tra indici PDO in monitoraggio reale, specialmente su scale stagionali-decadali. 

È qui che si inserisce il risultato centrale di Cluett et al. (2025): negli ultimi anni non osserviamo un “collasso” della PDO in senso stretto (né un cambiamento sostanziale del suo pattern spaziale o della sua forza), ma un cambio di regime nel rapporto tra variabilità interna e baseline termica. Gli autori mostrano che il riscaldamento su scala di bacino è diventato abbastanza grande da sovrastare la variabilità di SST tipicamente associata alla PDO, producendo periodi in cui le anomalie regionali si discostano da ciò che “ci aspetteremmo” guardando solo la fase PDO. Per descrivere questa baseline non stazionaria introducono un indice complementare, il pan-basin pattern (PBP), che rappresenta il segnale di riscaldamento diffuso del Nord Pacifico: le SST regionali risultano sempre più dalla sovrapposizione di PDO (variabilità interna) e PBP (spostamento della linea di base). In questo schema diventa naturale spiegare perché correlazioni storiche PDO–ecosistema possano indebolirsi o invertirsi: molti organismi rispondono a soglie e a valori assoluti (stress termico, mismatch trofico, habitat compression), per cui la stessa oscillazione “relativa” può produrre effetti diversi quando la baseline si è spostata verso stati più caldi. Il Research Briefing associato sintetizza bene la discontinuità: dal 2014 circa l’influenza della PDO sulle SST regionali risulta spesso superata da quella del riscaldamento di bacino, con impatti ecologici inattesi. 

Questa lettura è coerente con due filoni molto solidi della letteratura recente sugli estremi oceanici. Primo: l’aumento del contenuto di calore superficiale rende più probabile che la variabilità interna “spinga” il sistema oltre soglie critiche, favorendo condizioni no-analogue. Secondo: gli eventi di marine heatwave (definiti in modo standard come periodi prolungati in cui la SST supera una soglia stagionale locale, spesso il 90° percentile per almeno 5 giorni) sono aumentati in frequenza e durata su scala globale nel secolo scorso e sono destinati a intensificarsi in un oceano che continua a scaldarsi. Il caso del NE Pacific 2013–2014 (il “Blob”) è emblematico: un’anomalia calda persistente con impatti fisici ed ecosistemici ampiamente documentati, che ha contribuito a rendere intuitivo il concetto di “baseline spostata + variabilità interna = estremi persistenti”. In altre parole, se la PDO resta utile per descrivere una componente della dinamica interna, la gestione delle risorse marine (e la previsione di rischio) richiede sempre più di affiancare indici relativi con metriche legate agli assoluti termici e agli estremi, e di esplicitare sempre quale costruzione dell’indice si sta usando (ad esempio se si rimuove o meno la media globale, e con quale climatologia/dataset). Questo è esattamente il tipo di aggiornamento “operativo” che molte sintesi sulla climate-informed management considerano necessario per usare previsioni stagionali-decadali in modo robusto in un sistema non stazionario. 

Figura 1 — Pattern dominanti della variabilità delle SST nel Pacifico settentrionale: dalla PDO al riscaldamento pan-bacino

La Fig. 1 è costruita per rispondere a una domanda molto precisa: quando analizziamo la variabilità delle temperature superficiali del mare (SST) nel Nord Pacifico con una decomposizione EOF, i “modi dominanti” restano gli stessi nel tempo oppure cambiano perché la baseline termica si sta spostando? Gli autori calcolano le EOF su finestre mobili retrospettive di 30 anni (ogni riga), aggiornate ogni due anni, così da seguire l’evoluzione del sistema in modo quasi “cinematografico”. In ciascuna finestra mostrano i primi due modi (EOF1 ed EOF2), che per costruzione sono ortogonali e ordinati per varianza spiegata: EOF1 cattura la struttura spaziale che massimizza la varianza, EOF2 la successiva, e così via. È un punto cruciale: se nel dataset entra un segnale nuovo e molto energetico (ad esempio un riscaldamento coerente su larga scala), l’EOF1 tenderà naturalmente a “prenderselo”, anche se la dinamica interna sottostante non è cambiata.

Nelle prime due colonne (a sinistra) le EOF sono calcolate sulle anomalie totali di SST (“Total SST anomalies”), cioè senza rimuovere il contributo legato alla variazione della temperatura media globale del mare (GMSST). Nelle finestre più antiche (in alto) l’EOF1 mostra un pattern ben riconoscibile e fisicamente coerente con la PDO classica: una struttura dipolare su scala di bacino, con segni opposti tra settore orientale e settore centrale-occidentale del Pacifico settentrionale, la firma che storicamente è stata associata a cambiamenti atmosferici e oceanografici su vasta scala e a risposte ecosistemiche importanti (il caso emblematico è la riorganizzazione tardo-anni ’70 e i riflessi sulla pesca dei salmonidi). In altre parole, per una buona parte del periodo del tardo XX secolo la principale “geometria” della variabilità delle SST del Nord Pacifico, vista nelle anomalie totali, coincide con la geometria che chiamiamo PDO.

Scorrendo verso le finestre più recenti (righe basse), però, succede qualcosa di molto significativo: l’EOF1 delle anomalie totali perde progressivamente la struttura dipolare e tende a diventare quasi unipolare, con loadings dello stesso segno su gran parte del bacino. Nell’ultima finestra (1994–2023) questa configurazione è esplicitamente etichettata come PBP (Pan-Basin Pattern): non è più un “modo oscillatorio” che redistribuisce calore tra regioni, ma un segnale che rappresenta uno spostamento coerente della baseline termica del Nord Pacifico. Il senso fisico è semplice ma dirompente: quando il riscaldamento medio (o comunque una componente a grande scala, coerente) diventa abbastanza ampio, esso domina la varianza delle SST osservate e quindi “sale” al primo posto nella gerarchia EOF. Questo è esattamente il nucleo dell’articolo: la PDO può rimanere robusta, ma viene sovrastata in termini di varianza osservata dal riscaldamento pan-bacino. 

La prova più pulita arriva dalle due colonne di destra (“GMSST removed”), dove le EOF sono ricalcolate dopo aver rimosso le anomalie di GMSST. Questa procedura non è un artificio: è coerente con la costruzione canonica dell’indice PDO, che nasce proprio per isolare la variabilità interna del Nord Pacifico separandola dal segnale di riscaldamento su scala globale. Qui la figura mostra che il pattern della PDO non collassa e non “migra” in qualcosa di diverso: resta sostanzialmente stabile nelle varie finestre e continua a comparire come modo dominante o tra i primissimi modi. Questa è la distinzione concettuale decisiva proposta dagli autori: non è la PDO ad aver cambiato natura, è il contesto statistico-fisico in cui la osserviamo (SST totali vs SST depurate dal segnale globale) ad essere diventato non stazionario. 

Letta così, la Fig. 1 spiega in modo diretto perché negli ultimi anni alcune relazioni PDO-ecosistema abbiano mostrato indebolimenti o perfino inversioni di segno senza dover invocare “una PDO diversa”: se la SST regionale che interessa a un organismo (o a una catena trofica) è sempre più la somma di due contributi — una componente interna tipo PDO (dipolare) e una componente di baseline tipo PBP (unipolare) — allora può capitare che la fase PDO suggerisca un segnale “atteso”, ma la temperatura effettiva sia dominata dal riscaldamento di fondo, producendo condizioni termiche fuori scala rispetto al passato. In termini dinamici-statistici, la PDO resta un indice utile della variabilità interna (come discusso nelle sintesi classiche che sottolineano come la PDO sia un indice emergente da più processi, non un singolo meccanismo), ma la sua capacità di rappresentare da sola lo “stato termico” dell’oceano superficiale regionale diminuisce quando il trend/mean-state shift aumenta. 

Un ulteriore dettaglio importante della figura è che, mentre nelle SST totali l’EOF1 recente diventa PBP, la struttura dipolare “alla PDO” tende spesso a riapparire come EOF2 nelle finestre più recenti: questo è coerente con l’idea che la PDO continui a esistere ma venga declassata nella gerarchia della varianza perché la prima quota di varianza viene catturata dal riscaldamento coerente. Questo tipo di “riordinamento dei modi dominanti” sotto riscaldamento non è un risultato isolato nella letteratura sulla variabilità pacifica: diversi lavori recenti hanno discusso come la presenza di trend e cambiamenti del mean state possa indurre cambiamenti apparenti nelle modalità principali se non si separa correttamente la componente forzata dalla variabilità interna. 

Infine, la Fig. 1 si collega in modo naturale al concetto di regimi oceanici no-analogue e di crescita degli estremi: se la baseline del Nord Pacifico è più calda (PBP positivo) e sopra questa agisce ancora la variabilità interna (PDO e altri modi, come NPGO/Victoria Mode, spesso rilevanti per biogeochimica e dinamica dei gyre), allora la combinazione può produrre più facilmente eventi persistenti di caldo estremo e impatti ecosistemici non lineari. Il “Blob” del 2014 nel NE Pacific è un esempio noto di anomalia calda eccezionale con impatti fisici e biologici, utile qui come paradigma di cosa significa operare in un sistema dove variabilità e baseline si sommano su livelli termici già elevati. In sintesi, la figura formalizza un cambio di regime interpretativo: per leggere correttamente il Nord Pacifico contemporaneo (e usare indici climatici in gestione), serve distinguere con rigore la parte “PDO-like” (variabilità interna) dalla parte “PBP-like” (baseline in riscaldamento), perché nelle SST osservate entrambe concorrono a determinare lo stato termico regionale e quindi le risposte biologiche. 

Dominanza emergente del riscaldamento pan-bacino e “riordinamento” dei modi EOF nel Nord Pacifico (Fig. 1)

La sottosezione mette a fuoco un punto metodologico che, oggi, diventa anche un punto fisico: la PDO nasce come indice costruito per rappresentare variabilità interna del Nord Pacifico, e per farlo viene definita canonicamente come EOF1 delle anomalie di SST nel dominio 20–65°N, 105–260°E dopo rimozione del segnale di GMSST (temperatura media globale della superficie marina). Questa scelta non è un dettaglio “da statistici”: serve proprio a evitare che il cambiamento del mean state globale contamini il modo principale e a mantenere la PDO interpretativamente legata a meccanismi interni (forzante atmosferica extratropicale, teleconnessioni tropicali, integrazione oceanica del rumore atmosferico), come discusso nelle sintesi dinamiche che descrivono la PDO come pattern emergente da più processi e non come un singolo oscillatore fisico. 

Dentro questo quadro, il risultato “tranquillizzante” è che la PDO canonica risulta effettivamente stabile: quando si rimuove la GMSST e si calcolano le EOF su finestre mobili retrospettive di 30 anni, il pattern principale resta quello tipico PDO e la quota di varianza spiegata rimane in un intervallo relativamente stretto (circa 20–26% per EOF1). In parallelo, la EOF2 tende a riprodurre con continuità l’impronta della NPGO/Victoria Mode, anch’essa con varianza spiegata stabile (circa 12–16%), coerentemente con la letteratura che identifica la NPGO come modo complementare alla PDO, più direttamente connesso a dinamica dei gyres e variabili biogeochimiche/ecosistemiche. 

Il passaggio davvero cruciale, però, è ciò che accade quando non si rimuove la GMSST e si lavora sulle anomalie totali di SST (cioè includendo il riscaldamento di fondo): fino a prima del 2014 la differenza tra fare o non fare la rimozione è quasi “banale”, perché i pattern PDO ottenuti con e senza rimozione della GMSST risultano quasi identici (correlazioni spaziali molto alte, R ≈ 0,98–0,99) e la varianza spiegata cambia poco. Dopo il 2013, invece, emerge una discontinuità: l’EOF1 delle anomalie totali viene progressivamente occupata da un pattern unidirezionale (loadings dello stesso segno su quasi tutto il bacino), che gli autori definiscono PBP (Pan-Basin Pattern) nella finestra 1994–2023. In pratica, non è che la PDO “sparisca”: è che, nella metrica delle SST totali, la varianza associata al riscaldamento pan-bacino diventa così grande da sovrastare la varianza della PDO e quindi “sale” al primo posto nella gerarchia EOF. Questa sostituzione si riflette anche nei numeri: la varianza spiegata da EOF1 nelle anomalie totali cresce fino a valori dell’ordine 21–51%, superando qualunque intervallo storico trentennale in cui il primo modo fosse dominato dalla PDO. 

Qui vale la pena esplicitare il meccanismo statistico-fisico sottostante, perché è la chiave interpretativa della Fig. 1: le EOF non riconoscono “etichette climatiche”, ma catturano il modo che massimizza la varianza nel campione. Se nel campione entra un segnale coerente e crescente (il riscaldamento di fondo), l’EOF1 tenderà a rappresentarlo; la struttura dipolare PDO può allora riapparire come EOF2 (o comunque scendere di rango) senza che ciò implichi un cambiamento della dinamica interna. Questo tipo di lettura è coerente con la teoria degli errori campionari e con l’idea che, al variare della composizione del campione e della separazione tra autovalori, l’identificazione dei modi principali può “riordinarsi” anche a parità di processi fisici sottostanti—motivo per cui gli autori insistono su test di robustezza (stagione invernale soltanto, e finestre di 20–100 anni) mostrando che la transizione PDO→PBP resta confinata all’ultimo decennio anche cambiando ipotesi analitiche. 

Il fatto che l’indice PBP (PC1 delle anomalie totali 1994–2023) sia fortemente legato alla GMSST (correlazione lineare ~0,80) chiarisce la natura del segnale: non stiamo parlando di un “nuovo oscillatore decadale” indipendente, ma di un descrittore del baseline shift su scala di bacino, cioè della componente di riscaldamento che rende il sistema non stazionario rispetto al passato recente. Questo ha un riflesso diretto sulle applicazioni: la PDO, resa famosa anche dalla sua associazione a cambiamenti ecosistemici e a “regime shifts” (con la transizione 1976–77 come riferimento storico e impatti marcati su produzioni ittiche e catene trofiche), rimane utilissima come proxy di variabilità interna; ma se la temperatura regionale è sempre più la sovrapposizione tra PDO e PBP, allora le relazioni PDO-ecosistema possono indebolirsi o invertirsi semplicemente perché molte risposte biologiche dipendono da valori assoluti e soglie fisiologiche, non solo da anomalie relative. È esattamente il contesto in cui diventano più probabili condizioni “no-analogue” e periodi prolungati di caldo estremo, coerenti con la crescita e la formalizzazione del concetto di marine heatwave e con casi emblematici nel Nord-Est Pacifico (ad esempio l’anomalia calda del 2014) che hanno mostrato impatti fisici ed ecologici fuori scala rispetto alla climatologia recente. 

PDO robusta nonostante il riscaldamento predominante

In questa sottosezione gli autori fanno un’operazione concettualmente molto pulita: separano ciò che è davvero “cambiato” nel Nord Pacifico (la dominanza del riscaldamento secolare pan-bacino) da ciò che invece non è cambiato (la PDO come modo interno). Il punto di partenza è la definizione canonica: la PDO è il modo EOF principale delle anomalie di SST nel Nord Pacifico dopo aver rimosso la climatologia mensile e la media globale delle SST (GMSST), proprio per isolare la variabilità interna dal trend globale. 

Dentro questa cornice, il risultato chiave è che la PDO resta robusta “in senso dinamico-statistico”: nel periodo 1994–2023 il suo pattern non scompare, ma scende di rango quando si guardano le anomalie totali (senza rimozione della GMSST). In pratica, nelle SST totali la prima EOF viene occupata dal segnale pan-bacino (PBP), mentre il pattern tipico della PDO riemerge come EOF2; e la corrispondenza tra la PDO “EOF2 delle SST totali” e la PDO “EOF1 con GMSST rimossa” è quasi perfetta (correlazione spaziale R ≈ 0,99). Questo è il passaggio che taglia via molte interpretazioni fuorvianti: non c’è bisogno di ipotizzare una “nuova PDO” o una PDO indebolita; è la gerarchia della varianza a essere cambiata perché il segnale di riscaldamento comune è diventato dominante. 

La robustezza della PDO viene quantificata con misure che contano davvero in un’analisi EOF: gli autovalori associati alla PDO nel periodo più recente (EOF1 con GMSST rimossa e EOF2 nelle SST totali) risultano molto simili (circa 0,12–0,13) e rimangono nel range storico dell’EOF1 “de-trended” (circa 0,095–0,170), mentre le componenti principali (PC) corrispondenti restano fortemente correlate. In altre parole, l’energia statistica del modo PDO non collassa: semplicemente viene “messa in ombra” da un modo più grande, legato al riscaldamento del bacino. 

Il dettaglio più importante, però, è capire perché molti utenti “operativi” percepiscono una PDO più debole: non perché la PDO abbia perso ampiezza intrinseca, ma perché la quota percentuale di varianza che le viene attribuita nelle SST totali diminuisce (ad esempio dal ~23% al ~13% nell’ultimo decennio considerato dagli autori) mentre l’EOF1 (PBP) aumenta enormemente di importanza. È un effetto matematico inevitabile quando introduci nel dataset una componente coerente e crescente: le EOF massimizzano la varianza, quindi il riscaldamento pan-bacino tende a “prendersi” il primo modo e a comprimere in percentuale tutti gli altri, anche se in termini assoluti restano comparabili. Questo tipo di ri-organizzazione dei modi dominanti in presenza di un segnale di fondo persistente è coerente anche con risultati recenti che mostrano come periodi di marine heatwaves persistenti possano modificare la dominanza delle EOF di SST nel Nord Pacifico quando non si rimuove il trend/mean state. 

Un’altra conseguenza elegante del ragionamento è la “retrogradazione” del secondo grande modo interno: il footprint della NPGO/Victoria Mode, che storicamente compare spesso come EOF2, nel 1994–2023 passa a EOF3 nelle SST totali, pur spiegando una quantità di varianza assoluta simile a quella dei periodi precedenti (o quando si rimuove la GMSST). Anche qui il messaggio è identico: la NPGO non sparisce, ma viene spostata più in basso dalla comparsa di un modo di bacino molto energetico. Questo è perfettamente compatibile con la letteratura che tratta NPGO e PDO come modalità interne distinte e complementari, legate a dinamiche diverse (gyre circulation e variabili biogeochimiche nel caso NPGO, combinazione di forcing atmosferico e teleconnessioni nel caso PDO). 

In sintesi, la sottosezione formalizza un punto che per la gestione delle risorse e per l’interpretazione ecologica è decisivo: oggi il cambiamento dominante delle SST nel Nord Pacifico è il riscaldamento secolare, ma i due principali modi interni (PDO e NPGO/Victoria Mode) restano fisicamente presenti e statisticamente robusti quando li si isola correttamente dal segnale globale. Per questo, usare la PDO “come prima” senza esplicitare se si sta parlando di anomalie depurate dalla GMSST o di SST totali può portare a diagnosi sbagliate: la variabilità interna continua a modulare il sistema, ma lo fa sopra una baseline più calda che altera la lettura di segno/intensità e, soprattutto, la risposta biologica legata a soglie termiche.

La Fig. 2 mette “in numeri” ciò che la Fig. 1 mostrava in forma spaziale: quale modo EOF domina la varianza delle SST nel Nord Pacifico quando si usano finestre mobili di 30 anni, e come questa dominanza cambi a seconda che si lavori sulle anomalie totali (con dentro il riscaldamento globale) oppure sulle anomalie depurate dal segnale di GMSST. L’idea di fondo è che un’analisi EOF non identifica “meccanismi”, ma ordina i pattern in base alla varianza spiegata; se nel dataset entra un segnale sempre più grande e coerente (un riscaldamento diffuso su scala di bacino), quel segnale tende inevitabilmente a occupare l’EOF1, facendo scendere di rango i modi legati alla variabilità interna pur senza che essi si indeboliscano fisicamente. È esattamente il messaggio che Cluett et al. vogliono rendere evidente: la PDO resta robusta, ma nelle SST totali viene “oscurata” dalla crescita del Pan-Basin Pattern (PBP). 

Nel pannello superiore (“Total SST anomalies”) si vede la transizione in modo molto netto. Fino ai primi anni 2010, la EOF1 spiega grosso modo una quota stabile della varianza (ordine ~20–25%) ed è attribuita alla PDO: in quel periodo, l’organizzazione principale della variabilità delle SST osservate è ancora quella dipolare tipica del Nord Pacifico. Dopo la linea tratteggiata rossa (circa 2013–2015, indicata come momento di transizione), la varianza spiegata dalla EOF1 cresce rapidamente fino a valori molto più alti (nell’articolo viene riportato un intervallo che arriva a ~50% nelle finestre più recenti), e la modalità associata viene identificata come PBP, cioè un riscaldamento coerente su scala di bacino che domina la variabilità delle SST “così come sono”, senza rimozioni. In parallelo, nel pannello EOF2 delle anomalie totali, la PDO ricompare come secondo modo, ma la sua quota percentuale di varianza tende a ridursi nell’ultimo decennio, coerentemente con l’aumento schiacciante del PBP in EOF1: è la firma classica di un “riordinamento” dei modi EOF quando un segnale a grande scala aumenta la propria energia statistica. 

Il blocco inferiore (“GMSST removed”) è il controllo che chiude il cerchio: una volta rimossa la media globale delle SST, la gerarchia torna stabile. Qui la EOF1 rimane associata alla PDO con una varianza spiegata sostanzialmente costante nel tempo (nel testo: ~20–26%), mentre la EOF2 resta riconducibile al footprint della NPGO/Victoria Modecon una varianza anch’essa stabile (nel testo: ~12–16%). Questo comportamento è coerente con la definizione canonica e operativa della PDO, che nasce proprio come EOF principale delle SSTA nel Nord Pacifico dopo rimozione del contributo medio globale, così da descrivere variabilità interna distinta dai cambiamenti termici su scala planetaria. L’interpretazione fisica che ne deriva è molto precisa: non è la PDO ad essersi “rotta”, ma è la metrica delle SST totali ad essere diventata fortemente non stazionaria perché il riscaldamento di fondo è ormai il contributo dominante alla varianza; PDO e NPGO restano comunque i principali modi interni del sistema, in linea con la letteratura che descrive la PDO come pattern emergente da più processi (teleconnessioni tropicali, forcing atmosferico, integrazione oceanica del rumore) e non come un singolo oscillatore.

Rivalutare i cambiamenti climatici nel Pacifico settentrionale

La logica della sottosezione è costruire un ponte tra la “PDO canonica” (intesa come modo interno isolato rimuovendo la componente globale) e l’evoluzione storica completa delle SST osservate. Operativamente, gli autori prendono il pattern più recente della EOF1 con GMSST rimossa (1994–2023) e lo proiettano sull’intero archivio ERSST (1854–2024): così ottengono una ricostruzione coerente della componente PDO lungo tutta la serie strumentale, che riproduce le inversioni ben note dell’indice PDO utilizzato da NOAA/NCEI. Questo passaggio è importante perché mostra che la PDO “di oggi” resta pienamente riconoscibile e confrontabile con quella storica, in continuità con la formalizzazione classica del modo (regimi freddi/negativi 1890–1924 e 1947–1976; regimi caldi/positivi 1925–1946 e 1977–fine anni ’90), ma con una maggiore frequenza di inversioni sub-decadali dopo il 2000 (circa 2003, 2006, 2014, 2020), coerente con l’idea — ormai consolidata — che la PDO sia un pattern emergente dall’interazione tra teleconnessioni ENSO, “atmospheric bridge” e integrazione oceanica a bassa frequenza, più che un oscillatore periodico regolare. 

Dentro la stessa ricostruzione storica, il PBP (il segnale pan-bacino) cresce in modo pressoché monotono insieme alla GMSST, cioè lungo l’asse del riscaldamento globale. Qui la sottosezione introduce un punto delicato ma cruciale: la variabilità interna del Pacifico (inclusa la PDO) può modulare la temperatura media globale e, soprattutto, la sua pendenza su intervalli di 10–20 anni, attraverso riorganizzazioni della circolazione e della redistribuzione del calore in oceano (più accumulo subsuperficiale nel Pacifico tropicale durante fasi “La Niña-like”, con rallentamento apparente della crescita della GMST). Questo è esattamente il cuore della letteratura sul cosiddetto “hiatus” dei primi anni 2000: lavori osservazionali e modellistici hanno mostrato come un rafforzamento degli alisei e una maggiore efficienza di uptake/trasporto verso la termoclina e gli strati subsuperficiali possano ridurre temporaneamente il tasso di crescita della temperatura superficiale, pur in presenza di persistente squilibrio energetico del sistema climatico. In questa chiave, l’aumento “a gradini” del PBP attorno al 2014 può essere letto come combinazione di (i) fine/attenuazione del contributo interno che aveva smorzato la GMST e (ii) possibili cambiamenti del forcing aerosolico regionale, con evidenze recenti che collegano la riduzione degli aerosol in Asia orientale a un’accelerazione del riscaldamento post-2010, inclusi segnali nel Nord Pacifico. 

Sul piano spaziale e regionale, la sottosezione chiarisce perché PDO e PBP non siano intercambiabili: le anomalie “PDO-related” risultano più intense nel Pacifico orientale (carico maggiore nelle regioni nord-orientali), mentre le anomalie “PBP-related” sono più uniformi su scala di bacino. Anche se la somma lineare dei due contributi non deve coincidere perfettamente con le anomalie totali (perché un’analisi EOF produce modi ortogonali e non necessariamente fisicamente additivi in ogni punto e istante), il fatto che i primi due modi spieghino oltre il 60% della varianza recente implica che gran parte della variabilità osservata può essere interpretata come superposizione di: (a) variabilità interna tipo PDO e (b) traslazione della baseline termica tipo PBP. Questo si innesta direttamente sulla pratica operativa degli indici: l’indice PDO NOAA/NCEI, ad esempio, è calcolato proiettando le anomalie ERSST su una mappa di regressione legata al PDO “Mantua-like”, e nasce proprio per catturare il pattern interno dominante del Nord Pacifico; ma oggi, se non si separa esplicitamente il contributo di warming di fondo, lo stesso segno della PDO può non corrispondere più al segno dell’anomalia regionale assoluta di SST in aree chiave. 

La parte più applicativa arriva quando gli autori “calano” PDO e PBP nei Large Marine Ecosystems (LME). Nel Pacifico nord-orientale (California Current e Gulf of Alaska), prima del 2014 le anomalie di SST seguono prevalentemente la componente PDO: è coerente con il carico elevato della PDO in quell’area e con la tradizione di collegare le fasi PDO a riorganizzazioni ecosistemiche e a cambiamenti di produttività/aree di distribuzione (classicamente discussi già dalla letteratura sui regime shifts del 1976–77 e sugli impatti su salmonidi e reti trofiche). Nel Pacifico nord-occidentale (Kuroshio e Oyashio), invece, le SST storiche tendono a seguire di più il PBP, coerentemente con carichi PDO più deboli/variabili e con un segnale di warming più “coerente” su scala di bacino. Questo è un punto gestionale enorme: in aree dove la PDO “pesa” poco, aspettarsi risposte biologiche lineari a un indice PDO può essere fragile già in climatologia; con una baseline che si sposta rapidamente, diventa ancora più facile osservare correlazioni indebolite o invertite. 

Il passaggio finale, basato sui cross-plot PDO vs PBP, fornisce una griglia interpretativa molto concreta per gli ultimi anni: tra 2014 e 2016, PDO positiva e PBP elevato agiscono in modo costruttivo, producendo anomalie eccezionalmente calde e persistenti nel Nord-Est Pacifico, coerenti con la fase delle grandi marine heatwaves (incluso il “Blob”). Qui è utile ricordare che la comunità ha formalizzato metriche e definizioni operative di marine heatwave (durata, intensità, estensione) proprio per descrivere eventi persistenti sovrapposti a un trend di fondo; e il Blob 2013–2014/15 è diventato un caso di studio perché combina forzanti atmosferiche regionali, feedback oceano-atmosfera e contributi di teleconnessione tropicale che ne aumentano persistenza e impatti ecosistemici. Dopo il 2017 la PDO perde positività ma il PBP resta alto; e quando la PDO vira negativa attorno al 2021, il “raffreddamento PDO-related” nel Pacifico nord-orientale non basta più a compensare la baseline calda, generando per il 2021–2024 un regime nuovo nella storia recente: PDO negativa senza SST regionali negative nel Pacifico orientale. Questo è, in pratica, l’esempio più chiaro di “variabilità interna che continua a oscillare, ma sopra un pavimento che si alza”: la gamma della variabilità PDO resta, ma traslata lungo l’asse del cambiamento globale. 

Figura 3 — Serie temporali di PBP e PDO e loro impronte regionali: la SST come sovrapposizione tra variabilità interna e baseline in riscaldamento

La Fig. 3 è il passaggio in cui l’articolo rende “operativo” il nuovo schema concettuale: le anomalie regionali di temperatura superficiale del mare (SST) nel Pacifico settentrionale non sono più interpretabili come semplice funzione della fase PDO, perché lo stato termico osservato è sempre più la risultante della sovrapposizione tra un segnale di riscaldamento pan-bacino (PBP) e la variabilità interna rappresentata dalla PDO. Gli autori costruiscono entrambe le componenti con una procedura coerente: definiscono i pattern (EOF) nel periodo recente 1994–2023 e poi li proiettano sull’intero record ERSST (1854–2024), ottenendo due serie temporali “componibili” che permettono di confrontare storia lunga e periodo recente in modo omogeneo. 

Nei pannelli (a–c) si vedono i tre assi di riferimento. In (a) la GMSST (global mean sea-surface temperature) mostra la crescita di lungo periodo della temperatura superficiale oceanica media globale, che costituisce lo sfondo energetico su cui si innestano le oscillazioni decadali. In (b) la serie PBP (PC1) rappresenta il contributo del modo “unidirezionale” su scala di bacino: un segnale che cresce insieme alla GMSST ed esprime un cambiamento del mean state nel Nord Pacifico, non un’oscillazione redistributiva. In (c) la serie PDO ricostruita, ottenuta dalla EOF1 delle SST con GMSST rimossa (pattern 1994–2023 proiettato sul record completo), riproduce le note alternanze di regime e le inversioni sub-decadali nel XXI secolo; il legame con gli indici operativi è esplicitato perché la costruzione risulta coerente con l’indice PDO di NOAA/NCEI basato su ERSST, che viene calcolato tramite proiezione di anomalie ERSST su una mappa di regressione legata al PDO “Mantua-like”. 

Il pannello (c) va letto con un’attenzione particolare: qui non c’è alcuna evidenza di una “PDO scomparsa”; al contrario, la PDO continua a oscillare con ampiezza comparabile, in linea con la letteratura dinamica che descrive la PDO come pattern emergente dall’interazione tra forzante atmosferica extratropicale, teleconnessioni tropicali (ENSO/atmospheric bridge) e integrazione oceanica del forcing, più che come un oscillatore singolo e autonomo. Questa distinzione è decisiva perché spiega perché un indice PDO possa rimanere fisicamente robusto pur perdendo capacità di “predire” il segno delle SST assolute in alcune regioni: non cambia la variabilità interna, cambia il riferimento termico su cui essa si somma.

Nei pannelli (d–g) la figura entra nel dominio applicativo, confrontando per quattro Large Marine Ecosystems (LMEs) le anomalie totali osservate (linea nera continua) con le componenti previste associate a PBP (area gialla) e PDO (componenti rosso/blu), più la loro somma (linea nera tratteggiata). Il fatto che la curva tratteggiata segua spesso bene la curva continua è una conseguenza del risultato strutturale mostrato nelle figure precedenti: i primi due modi (PBP e PDO, nelle SST totali) catturano una frazione molto grande della varianza recente, quindi la decomposizione “spiega” gran parte della dinamica delle SST regionali anche se non può coincidere perfettamente punto per punto (perché le EOF sono ortogonali e il residuo include ulteriori modi e processi locali). 

Nel Pacifico orientale, GOA (d) e California Current (e), la PDO ha storicamente carichi elevati e quindi, fino a prima del 2014, le SST regionali tendono a seguire più da vicino la componente PDO. Dopo il 2014, però, la componente PBP diventa persistentemente positiva e “alza il pavimento”: quando PDO e PBP sono dello stesso segno, la somma produce anomalie eccezionalmente calde e durature. È in questo quadro che si collocano le grandi marine heatwaves del Nord-Est Pacifico, incluso il periodo 2014–2016 associato al “Blob”, che la letteratura ha ricondotto a una combinazione di forzante atmosferica anomala (ridotta perdita di calore e mixing) e feedback oceano-atmosfera, con impatti ecosistemici di ampia portata. La cornice metodologica con cui questi eventi vengono oggi definiti e confrontati è quella delle marine heatwaves come superamento del 90° percentile per almeno 5 giorni (su climatologia di riferimento), che rende particolarmente chiaro come un baseline shift aumenti la probabilità di superare soglie estreme anche a parità di variabilità interna. 

Il tratto più “nuovo” della figura, coerente con la tesi dell’articolo, si vede nel periodo 2021–2024: quando la PDO entra in fase negativa, ci si aspetterebbe un raffreddamento marcato nel Pacifico orientale; invece, la componente PBP resta così elevata che il contributo raffreddante associato alla PDO non basta più a portare le SST regionali sotto lo zero. In termini diagnostici, è il primo intervallo prolungato del record recente in cui una PDO “fredda” non coincide con SST negative nel Nord-Est Pacifico, cioè la relazione intuitiva “PDO− ⇒ est-Pacific cool” viene confusa da un riscaldamento di fondo ormai dominante. 

Nel Pacifico occidentale, Oyashio Current (f) e Kuroshio Current (g), la figura mostra un comportamento diverso: qui le anomalie PBP risultano più determinanti e la PDO tende a modulare in modo più debole/variabile, coerentemente con i carichi regionali discussi dagli autori. È per questo che possono comparire più facilmente periodi di disallineamento tra segno PDO e segno dell’anomalia totale: in aree dove la PDO non domina la varianza regionale, il segnale “assoluto” tende a seguire più da vicino la componente pan-bacino. Questa asimmetria est-ovest è esattamente ciò che motiva l’uso combinato di indici: PDO per la variabilità interna e PBP per la traslazione della baseline.

La figura collega poi questi risultati a un tema più ampio: la PDO non è solo un indice regionale, ma può contribuire a modulare la temperatura media globale su scale interdecadali tramite redistribuzione del calore nel Pacifico tropicale e variazioni dell’uptake oceanico subsuperficiale. Il “hiatus” dei primi anni 2000 è stato interpretato in questa chiave da più studi, che hanno collegato il rallentamento temporaneo del riscaldamento superficiale a condizioni tipo La Niña e a un rafforzamento degli alisei con maggiore accumulo di calore sotto la superficie. In parallelo, l’articolo richiama l’ipotesi che la recente accelerazione del riscaldamento e i segnali nel Pacifico possano essere stati amplificati anche da cambiamenti del forcing aerosolico, con lavori recenti che attribuiscono un ruolo importante alla riduzione degli aerosol in Asia orientale nel rafforzare lo squilibrio radiativo e i trend di riscaldamento, incluso il Nord Pacifico. 

Il pannello (h) completa la lettura localizzando i quattro LMEs e rendendo immediata la conclusione fisica della Fig. 3: la variabilità associata alla PDO conserva la propria escursione, ma “scorre” lungo un asse di cambiamento globale sempre più spostato verso il caldo. In altri termini, la PDO continua a organizzare una parte importante della variabilità interna (come da inquadramento dinamico classico), ma l’interpretazione delle SST regionali — e quindi di impatti biologici, stress termici e rischio di eventi estremi persistenti — richiede ormai di distinguere esplicitamente tra componente interna e baseline in riscaldamento, perché è la loro combinazione a determinare i nuovi regimi “no-analogue” osservati nell’ultimo decennio. 

Risposte ecologiche non stazionarie alla variabilità della PDO in un Nord Pacifico che si riscalda

Per molti decenni la Pacific Decadal Oscillation (PDO) è stata trattata come una sorta di “chiave di lettura” ecologica del Pacifico settentrionale, perché il suo pattern di SST e le teleconnessioni atmosferico-oceaniche associate si sono spesso accompagnati a cambiamenti coerenti in produttività, distribuzioni di specie e rese di pesca. Il caso storico più citato è quello dei salmonidi, dove le transizioni di regime del tardo XX secolo hanno contribuito a consolidare l’idea che un indice climatico di bacino potesse fungere da proxy gestionale per dinamiche ecosistemiche complesse. 

La sottosezione, però, spiega perché quelle correlazioni non sono “leggi di natura” ma relazioni contingenti: negli ultimi anni alcune associazioni PDO–ecosistema si sono indebolite o addirittura invertite, e il lavoro di Cluett et al. propone che la causa primaria non sia un collasso della PDO (che resta robusta), bensì la crescita di un segnale di riscaldamento pan-bacino (PBP) che sposta la baseline termica e rende non stazionario il rapporto tra indice e impatto biologico. In questo quadro, il caso del Gulf of Alaska (GOA) diventa esemplare: le catture di pink, sockeye e coho salmon risultavano fortemente correlate positivamente con la PDO nel 1965–1988, poi la relazione si attenua nel 1989–2013 e diventa significativamente negativa dopo il 2014, in sincronia con l’aumento della predominanza del PBP. 

Il cuore interpretativo è un modello concettuale volutamente semplice ma molto potente: una risposta biologica (abbondanza, sopravvivenza, crescita) può essere descritta come funzione non lineare della temperatura con un optimum termico, spesso rappresentabile tramite thermal performance curves (curve di prestazione) la cui forma è tipicamente unimodale, con un ramo crescente fino all’optimum e un ramo decrescente oltre soglie di stress. Questa idea è solidamente radicata nella fisiologia degli ectotermi e nella teoria ecologica della sensibilità termica: piccole variazioni di temperatura attorno all’optimum possono aumentare performance e fitness, mentre superare l’optimum può innescare cadute rapide per stress fisiologico e vincoli sistemici (inclusi limiti legati al trasporto di ossigeno e alla capacità aerobica nei pesci). 

Qui entra il passaggio chiave della non stazionarietà: se l’ampiezza della variabilità “PDO-related” resta sostanzialmente comparabile nel tempo, ma la baseline assoluta di SST sale per effetto del PBP, allora lo stesso valore dell’indice PDO corrisponde a temperature assolute diverse in epoche diverse. In termini geometrici, la variabilità PDO “scorre” lungo la curva di risposta mentre il riscaldamento secolare avanza. Prima del 1989, nel GOA, la banda termica esplorata dalla PDO avrebbe campionato soprattutto il ramo crescente (scenario A): PDO positiva → SST più alte ma ancora sotto l’optimum → risposta biologica positiva e quindi correlazione positiva. Con l’aumento della baseline (scenario B), le fasi PDO positive iniziano più spesso a superare l’optimum, producendo risposte negative o più deboli e quindi decorrelazione nel 1989–2013. Quando la baseline oltrepassa stabilmente l’optimum (scenario C), la variabilità PDO cade prevalentemente sul ramo discendente: PDO positiva → ancora più caldo → risposta negativa → inversione della correlazione dopo il 2014. 

Questo schema spiega bene anche perché possano indebolirsi correlazioni già negative: se il riscaldamento spinge il sistema verso la parte più “piatta” del ramo discendente, ulteriori aumenti di SST durante PDO positive producono incrementi marginali minori dell’impatto biologico, e la relazione statistica perde forza. Inoltre, il modello rende naturale la divergenza tra regioni: nel California Current (vicino al margine meridionale di distribuzione per diverse popolazioni), temperature di base più alte collocano la variabilità più “avanti” lungo la curva, rendendo più plausibili risposte storicamente diverse (o più deboli) rispetto al GOA. Il fatto che i salmonidi abbiano optimum e finestre termiche specie- e stadio-specifiche, con stress e mortalità crescenti oltre soglie, è ampiamente coerente con la letteratura recente su vulnerabilità e limiti termici dei salmoni in un clima che si riscalda. 

Detto questo, la sottosezione è molto chiara su un punto che spesso viene sottovalutato quando si parla di “temperatura”: non tutte le risposte biologiche alla PDO sono guidate primariamente dalle SST assolute. La PDO porta con sé anche variazioni di upwelling, stratificazione, profondità della termoclina, disponibilità trofica e tempistiche fenologiche; in questi casi le correlazioni PDO–biologia possono restare simili, cambiare in modo non monotono o diventare più complesse. In altre parole, la non stazionarietà può emergere sia per “traslazione termica” (baseline shift) sia per cambiamenti nella catena causale dominante (ad esempio da controllo termico diretto a controllo trofico o idrodinamico), e il modello proposto va letto come cornice diagnostica, non come spiegazione unica. 

Un aspetto che rende questa cornice particolarmente attuale è l’aumento di condizioni estreme persistenti, cioè le marine heatwaves, formalizzate operativamente come periodi in cui la SST supera una soglia stagionale locale (spesso il 90° percentile) per almeno 5 giorni: quando una baseline più calda (PBP) e una fase PDO favorevole al riscaldamento agiscono in modo costruttivo, la probabilità di superare soglie ecologicamente critiche cresce in modo non lineare. Il “Blob” del Nord-Est Pacifico (2013–2014 e persistenze successive) è un caso di studio classico di anomalia calda con forti impatti fisici ed ecosistemici, e nel framework PBP+PDO diventa un esempio quasi didattico di somma costruttiva su un livello termico già elevato. 

La conclusione operativa implicita (e la direzione di ricerca che gli autori suggeriscono) è che, dopo il 2014, usare la sola PDO come “indice ecologico” rischia di perdere informazione perché confonde variabilità interna e spostamento del mean state. Integrare PDO + PBP (o, in modo equivalente, combinare indici di variabilità interna con metriche di SST assoluta, stress termico e marine heatwaves) potrebbe recuperare parte della capacità inferenziale, soprattutto se ancorato a meccanismi: optimum termici specie-specifici, finestre di tolleranza, vincoli fisiologici e processi oceanografici (upwelling/termoclina). È un passaggio coerente con la traiettoria più ampia della letteratura su salmoni e cambiamento climatico, che insiste su risposte non lineari, differenze regionali e importanza di soglie e condizioni estreme.

Figura 4 — Un quadro geometrico per leggere insieme variabilità interna (PDO) e riscaldamento di fondo (PBP) nelle SST regionali

La Fig. 4 traduce in modo estremamente chiaro l’idea centrale dello studio: nel Nord Pacifico contemporaneo l’anomalia di SST osservata in una regione non è più interpretabile come “PDO pura”, ma come risultante della sovrapposizione tra un contributo di variabilità interna (PDO) e un contributo di cambiamento climatico di fondo(PBP, pan-basin warming). Nel pannello (a) questa sovrapposizione viene rappresentata come un piano cartesiano in cui l’asse x misura l’anomalia di SST attribuibile al PBP (spostamento della baseline lungo l’asse del riscaldamento globale), mentre l’asse y misura l’anomalia attribuibile alla PDO (modulazione interna). In questa geometria, i quadranti alto-destro e basso-sinistro corrispondono a combinazioni “costruttive” (PDO e PBP dello stesso segno: il segnale finale si amplifica), mentre alto-sinistro e basso-destro corrispondono a combinazioni “distruttive” (segni opposti: compensazione parziale). La figura aggiunge poi un elemento diagnostico decisivo: non basta sapere il segno, occorre capire chi domina in ampiezza. Per questo vengono evidenziate due regioni: la “clessidra” (hourglass), in cui la magnitudine del contributo PDO supera quella del PBP e quindi la variabilità interna controlla l’anomalia finale, e il “papillon” (bowtie), in cui il contributo PBP è maggiore e quindi il cambiamento di fondo sovrasta la PDO. Questo schema formalizza la non-stazionarietà discussa nell’articolo: la PDO resta un modo interno robusto, ma la crescente importanza del PBP cambia la gerarchia causale delle SST osservate e, di conseguenza, la leggibilità di molte relazioni PDO-ecosistema che erano state storicamente forti. 

I pannelli (b–e) applicano lo schema a quattro Large Marine Ecosystems (GOA, California Current, Oyashio Current e Kuroshio Current), riportando per ogni anno un punto la cui ascissa è l’anomalia “PBP-related” e l’ordinata è l’anomalia “PDO-related”; il colore del punto rappresenta l’anomalia di SST osservata. Questa scelta grafica è molto potente perché consente di vedere simultaneamente (i) la combinazione dei due contributi, (ii) il regime costruttivo/distruttivo e (iii) il passaggio da dominio PDO a dominio PBP. Nel Pacifico nord-orientale (GOA e California Current; pannelli b e c) la nube dei punti pre-2014 tende a concentrarsi lungo l’asse PDO (dinamica coerente con i carichi PDO storicamente più forti in quell’area), mentre dal 2014 in avanti i punti si spostano nettamente verso x positive, indicando che il PBP porta la regione dentro un contesto di riscaldamento persistente. Gli anni etichettati 2014–2016 cadono tipicamente nel quadrante alto-destro: è la somma costruttiva tra PDO positiva e PBP positivo che genera anomalie eccezionalmente calde e durature, coerenti con la fase delle grandi marine heatwaves del Nord-Est Pacifico, inclusa la sequenza associata al “Blob”. La letteratura sul Blob mostra infatti come un’anomalia calda persistente nel Nord-Est Pacifico possa emergere e mantenersi quando anomalie atmosferiche e condizioni oceaniche riducono la perdita di calore e il rimescolamento, e in un oceano con baseline più calda queste configurazioni diventano più efficaci nel produrre estremi. 

Il dettaglio più rivelatore della Fig. 4, però, è ciò che accade quando la PDO cambia segno: negli anni recenti (in particolare nel 2021–2024, come discusso nel testo) compaiono punti che pur trovandosi nel settore “distruttivo” (PDO negativa, quindi tendenza al raffreddamento regionale, ma PBP ancora positivo) mantengono colori compatibili con anomalie osservate spesso ancora positive. In termini fisici, questo significa che il contributo raffreddante della PDO non è più sufficiente a compensare un PBP elevato: è l’evidenza grafica del passaggio a un regime in cui la fase PDO “fredda” non garantisce più SST regionali fredde nel Pacifico orientale, cioè viene meno una regola empirica che aveva sostenuto molte inferenze ecosistemiche e gestionali. Questa dinamica è coerente con la definizione moderna di marine heatwave come evento anormalmente caldo che supera soglie percentile su una climatologia, perché un baseline shift verso temperature più alte aumenta la probabilità che anche una modulazione interna sfavorevole non riesca a riportare il sistema sotto soglie “fredde”. 

Nel Pacifico nord-occidentale (Oyashio e Kuroshio; pannelli d ed e) la nube dei punti è più compressa lungo l’asse PDO e più traslata lungo x, segnalando che qui il contributo PBP tende a dominare più spesso e che il ruolo della PDO è più debole/variabile. In pratica, lo schema spiega perché in queste regioni il segno della PDO possa disallinearsi più facilmente dal segno dell’anomalia totale osservata: quando i carichi PDO sono meno determinanti, la componente di riscaldamento pan-bacino diventa il driver principale del segno e dell’intensità della SST regionale. Ne esce un messaggio molto concreto: la variabilità interna (PDO) conserva una gamma di oscillazione che rimane riconoscibile, ma questa oscillazione si muove su un “piano inclinato” imposto dal riscaldamento di fondo (PBP), e la figura mostra letteralmente questa traslazione come migrazione della nube dei punti verso valori PBP più positivi. È esattamente questo meccanismo di sovrapposizione a rendere non stazionarie molte correlazioni clima-ecosistema storicamente associate alla PDO (celebrate già dagli studi fondativi sui legami tra oscillazioni interdecadali e produzione di salmone), perché la stessa fase PDO può corrispondere oggi a temperature assolute ed estremi ecologicamente molto diversi rispetto al passato. 

Discussione — Dal “PDO-centrismo” al quadro PDO+PBP: non-stazionarietà, estremi e implicazioni operative nel Nord Pacifico

Il punto di arrivo della discussione è netto: nel Nord Pacifico la modalità dominante delle variazioni interdecadali di SST non è più, in prima approssimazione, la classica impronta della PDO, ma un segnale a segno coerente su scala di bacino, definito dagli autori come Pan-Basin Pattern (PBP). In termini diagnostici questo significa che, lavorando sulle anomalie “totali” (senza sottrarre il contributo globale), l’EOF1 viene progressivamente occupata dal riscaldamento pan-bacino; e, soprattutto, per la prima volta nel record recente si osserva un periodo prolungato in cui una PDO negativa “fredda” non coincide con SST fredde nel Pacifico orientale, cioè viene meno un’associazione empirica storicamente molto robusta (si veda la divergenza mostrata in Fig. 3d,e e il quadro compositivo di Fig. 4). 

Un elemento cruciale, che evita letture fuorvianti, è che questa “perdita di dominanza” non equivale a un collasso fisico della PDO: la discussione sottolinea che la variabilità PDO resta robusta e che i driver che ne sostengono l’impronta di SST hanno mostrato resilienza almeno fino all’inizio del XXI secolo. Ciò è coerente con una visione ormai consolidata nella letteratura dinamica: la PDO non è un oscillatore singolo, ma un pattern statisticamente emergente che integra più processi, tra cui forcing atmosferico extratropicale (stocasticità associata all’Aleutian Low), teleconnessioni dal Pacifico tropicale (atmospheric bridge legato a ENSO) e risposta/integrale oceanico tramite onde di Rossby e aggiustamenti a bassa frequenza. In questo quadro, la “futura evoluzione” della PDO dipende da come il riscaldamento modificherà, cumulativamente e su scale diverse, ciascun anello della catena (forzante atmosferica, accoppiamento tropici-extratropici, tempi di aggiustamento oceanico). 

La discussione poi sposta l’attenzione dal “perché la PDO esiste” al “come gli ecosistemi rispondono”, mettendo in evidenza una asimmetria di conoscenza: i driver fisici della variabilità decadale nordpacifica sono relativamente ben inquadrati, mentre i meccanismi biofisici che trasformano un segnale climatico di bacino in risposta ecologica restano più incerti e, soprattutto, più eterogenei tra specie, regioni e livelli trofici. Storicamente la PDO ha funzionato come proxy affidabile di SST in ampie porzioni del bacino (ed è alla base di molta narrativa sui regime shifts e sulle riorganizzazioni delle risorse, salmoni inclusi), ma se oggi le SST regionali possono divergere dal segno “atteso” sulla base della PDO, allora diventa indispensabile separare gli effetti dovuti alle SST assolute (baseline più alta, soglie superate più spesso) dagli effetti dovuti alle dinamiche correlate alla PDO ma non riducibili alla temperatura (upwelling, trasporto, profondità della termoclina, masse d’acqua sorgente e cascata biogeochimica). In questo senso la PDO potrebbe rimanere ecologicamente utile non tanto come “termometro regionale”, quanto come indicatore di specifici processi dinamici che continuano a oscillare anche sopra un background più caldo. 

Un corollario diretto è l’aumento atteso della frequenza di condizioni che, nel passato, sarebbero state classificate come “eccezionali” (ad esempio >2 deviazioni standard) perché oggi la variabilità interna si somma a un baseline in crescita. Il periodo 2014–2020 è l’esempio didattico: quando PBP (riscaldamento di fondo) e contributo PDO si combinano in modo costruttivo, la probabilità di estremi caldi persistenti cresce rapidamente, e il sistema entra più facilmente nel regime delle marine heatwaves. Qui è utile ricordare che la definizione operativa di marine heatwave più adottata identifica eventi in cui la SST supera una soglia stagionale (spesso il 90° percentile) per almeno 5 giorni: con una baseline più calda, la stessa oscillazione interna supera più spesso le soglie percentile e aumenta durata/intensità degli eventi. L’episodio del “Blob” nel Nord-Est Pacifico (2013–2014 e persistenze successive) resta un riferimento fondamentale perché documenta sia le condizioni fisiche favorevoli (forzante atmosferica e flussi di calore/mixing) sia gli impatti ecosistemici a catena; e nel quadro PDO+PBP diventa un prototipo di estremo generato da superposizione tra variabilità e trend. 

Per rendere quantitativo il rischio di estremi derivante dall’interazione tra variabilità interna e riscaldamento antropogenico, la discussione invoca esplicitamente l’esigenza di “vincolare” ampiezza e stabilità della variabilità interna: è un tema metodologico centrale, perché senza una stima robusta della variabilità naturale (su scale decadali e multidecadali) è facile sovra- o sotto-attribuire singoli episodi estremi. Qui entrano in gioco due filoni complementari. Da un lato, gli archivi paleoclimatici (tree rings, coralli, ecc.) estendono la prospettiva oltre il breve record strumentale e aiutano a delimitare l’escursione naturale di indici tipo PDO su secoli, fornendo un contesto indispensabile per interpretare “quanto sia anomalo” il presente. Dall’altro, i large ensembles (ad esempio CESM-LE) e gli ensemble inizializzati permettono di separare forzante e rumore interno, campionando molte realizzazioni possibili del clima e quantificando quanta parte delle tendenze regionali e degli estremi possa emergere dalla variabilità interna a parità di forcing. 

La discussione tocca anche un punto che negli ultimi anni sta diventando sempre più rilevante per il Pacifico: la possibilità che cambiamenti nel forcing aerosolico e nelle interazioni aerosol-nubi contribuiscano ad accelerazioni regionali del riscaldamento oceanico e a condizioni favorevoli a grandi anomalie calde. Esistono lavori recenti che collegano riduzioni delle emissioni aerosoliche post-2011 a un rapido riscaldamento nel Pacifico nord-occidentale tramite variazioni di copertura nuvolosa e radiazione a onda corta superficiale; e altri che propongono un ruolo delle riduzioni degli effetti indiretti degli aerosol nel favorire eventi tipo “warm blob” nel Nord-Est Pacifico. Questi risultati non sostituiscono la dinamica interna (PDO/ENSO), ma rafforzano l’idea della discussione: i regimi “no-analogue” emergono spesso dall’interazione di più fattori (variabilità interna + forzanti antropogeniche), e richiedono un’interpretazione multicomponente. 

Il messaggio finale è operativo e, in un certo senso, “scomodo”: decisioni di gestione basate su baseline storiche e su correlazioni PDO-ecosistema calibrate nel XX secolo rischiano di fallire perché non catturano le interazioni emergenti tra warming e variabilità interna. Per conservare informazione predittiva, la discussione suggerisce di spostare l’enfasi da indicatori empirici usati come scorciatoie a approcci che incorporino modellistica dinamica e previsione del clima nordpacifico, capaci di distinguere contributi termici (baseline) da contributi dinamici (upwelling/trasporto/masse d’acqua). Ed è una conclusione che probabilmente non vale solo per il Pacifico: se perfino in un bacino con variabilità interna molto forte il riscaldamento pan-bacino ha iniziato a superare la PDO, è plausibile che una riconsiderazione simile degli indici di variabilità climatica diventi necessaria anche in altri oceani. 

Dalla PDO “classica” al Pacifico che cambia: perché le relazioni PDO–biologia possono indebolirsi, sparire o invertirsi con il riscaldamento secolare (spiegazione della Fig. 5 con riferimenti)

La figura 5 mette in scena un’idea semplice ma potentissima: il fatto che un indice climatico (la PDO) sia ancora “vivo” e variabile non garantisce affatto che mantenga nel tempo la stessa capacità di spiegare (o prevedere) una risposta biologica, perché l’ecosistema non risponde alle anomalie in astratto, ma allo stato termico assoluto e al contesto fisico-biogeochimico in cui quelle anomalie si innestano. La PDO, nella sua definizione canonica, è il principale modo empirico (EOF1) delle anomalie di SST nel Nord Pacifico, calcolato rimuovendo il segnale di temperatura media globale per isolare la variabilità interna; in fase positiva tende ad associare acque relativamente più calde nel Pacifico orientale e più fredde nel centrale/occidentale (e viceversa in fase negativa). In altre parole, la PDO nasce come descrizione di un “pattern” di variabilità interna del sistema oceano–atmosfera (forzante aleutinica, flussi aria–mare, teleconnessioni tropicali, onde di Rossby oceaniche), non come termometro del riscaldamento di lungo periodo. 

Il punto nuovo (e qui entra in gioco l’articolo da cui proviene la figura) è che, nel Pacifico contemporaneo, il segnale di riscaldamento “pan-bacino” può diventare così dominante da spostare la baseline su cui la PDO oscilla. Cluett et al. (2025) mostrano che, se guardi alle anomalie totali di SST senza sottrarre il contributo globale, il modo principale recente non assomiglia più alla PDO ma a un pattern quasi unidirezionale di riscaldamento su tutto il bacino nordpacifico, chiamato pan-basin pattern (PBP). Questa emergenza del PBP fa sì che, per la prima volta in modo prolungato, una PDO negativa (“fredda” secondo l’indice) non coincida necessariamente con SST realmente fresche nel Pacifico orientale: il segnale di fondo è salito e la variabilità interna si sovrappone a una piattaforma termica più calda. 

Dentro questo quadro, la figura 5 propone un modello concettuale che usa una curva di risposta biologica “a campana” rispetto alla SST (una versione semplificata delle thermal performance curves, TPC): esiste un intervallo di temperatura vicino a un optimum in cui crescita/sopravvivenza/reclutamento (o la produttività dell’ecosistema che sostiene i salmoni) è massima, mentre ai lati—troppo freddo o troppo caldo—la performance cala. Questo tipo di risposta non è un’invenzione grafica: è un pilastro della fisiologia ecologica e della teoria delle TPC (anche se nella realtà le curve sono spesso asimmetriche e dipendono dallo stadio vitale, dall’ossigeno e dall’alimentazione). 

A questo punto la figura fa un passo cruciale: assume che l’ampiezza della variabilità associata alla PDO (cioè quanto “oscilla” la SST per effetto PDO) resti più o meno comparabile, ma che il riscaldamento secolare (PBP) trasli l’intero intervallo di SST assolute verso destra (più caldo). Se le oscillazioni PDO si collocano sul ramo ascendente della curva (scenario A), allora quando la PDO aumenta la SST si avvicina all’optimum e la risposta biologica cresce: la correlazione PDO–biologia risulta positiva (m>0). Gli autori collegano questo caso a quanto osservato storicamente per le catture di salmoni nel Golfo dell’Alaska (GOA) prima della fine degli anni ’80, periodo in cui la PDO è stata effettivamente un indicatore “forte” di riorganizzazioni ecosistemiche e produttività alieutica nel Nord Pacifico (basti pensare ai lavori fondativi che hanno reso famosa la PDO proprio in relazione al salmone). 

Quando però la baseline si scalda e l’intervallo PDO finisce vicino al massimo della curva (scenario B), la variabilità della PDO continua a esserci, ma piccoli spostamenti di temperatura producono variazioni biologiche relativamente piccole e spesso dominate da altri fattori: ne risulta una decorrelazione (m≈0). Questo è coerente con l’idea di relazioni tempo-dipendenti e con quanto documentato per il GOA nel periodo 1989–2013: dopo il noto “climate shift” di fine anni ’80 (1988/89) molte associazioni tra PDO e variabili regionali/biologiche si indeboliscono o cambiano, perché il sistema atmosfera–oceano entra in un regime diverso di combinazioni tra SLP, mixing, trasporto costiero, ecc. 

Infine, se il riscaldamento prosegue e l’intervallo di SST associato alla PDO viene traslato oltre l’optimum (scenario C), la stessa variabilità PDO si colloca sul ramo discendente: PDO più alta equivale ora a condizioni più calde “di troppo” e dunque a una risposta biologica peggiore, invertendo il segno della correlazione (m<0). Ed è qui che la figura si incastra con l’evidenza degli anni recenti: nel GOA, durante la fase di riscaldamento estremo 2014–2019, è stato osservato un passaggio a correlazioni negative PDO–produzione/catture di salmone, interpretato come manifestazione di un “clima novel” (combinazioni di variabili fisiche senza analoghi nel record recente) in cui l’indice PDO non mappa più lo stesso pacchetto di forzanti atmosferiche e oceaniche di prima. Litzow et al. (2020) mostrano esplicitamente che in quel periodo il PDO si associa a una relazione più debole con la pressione/forzante aleutinica e, soprattutto, a condizioni regionali più calde, meno mixing e trasporti diversi, insieme a un’inversione della relazione PDO–salmone rispetto alle decadi precedenti. 

Un elemento chiave del “perché adesso” è che il 2014–2016 (e poi altri episodi) coincide con grandi anomalie persistenti nel Nord-Est Pacifico, inclusa la famosa marine heatwave “Blob”, che ha una letteratura specifica su cause e impatti: pattern di pressione anomali, riduzione degli scambi turbolenti e persistenza del calore superficiale. In termini del framework PBP+PDO, questi estremi diventano più probabili quando variabilità interna e riscaldamento di fondo si sommano in modo costruttivo, creando stati oceanici “fuori scala” rispetto al passato recente. 

Il valore aggiunto della figura 5, quindi, non è dire “la PDO non serve più”, ma l’esatto contrario: la PDO resta un descrittore importante della variabilità interna (e infatti la sua dinamica e la sua impronta statistica possono rimanere robuste), però da sola rischia di diventare un proxy ingannevole se la domanda che stai facendo è biologica o gestionale e dipende dalla temperatura assoluta e dalle combinazioni di forzanti. In questa prospettiva ha senso affiancare indici complementari (come PBP, o anche altri modi come NPGO/Victoria Mode che catturano aspetti diversi della variabilità del Nord Pacifico) e, soprattutto, lavorare più su variabili “primarie” (SST assoluta, stratificazione, mixing, trasporti, nutrienti, fenologia del plancton) che sugli indici come scorciatoia.

Metodi

Il cuore dei metodi sta nel separare (e confrontare) due cose che spesso, nella pratica, vengono “fuse” nello stesso indice: la variabilità interna del Nord Pacifico e lo spostamento della baseline dovuto al riscaldamento globale. Per farlo, gli autori partono da una serie storica molto lunga e omogenea di SST, l’ERSSTv5 di NOAA (griglia 2°×2°, mensile, 1854–oggi), un prodotto “reconstructed” costruito combinando osservazioni (soprattutto ICOADS) con tecniche statistiche che colmano i vuoti spaziali e temporali; è proprio questo tipo di dataset che consente di studiare trend secolari e pattern su bacino con una copertura ragionevole lungo il XX secolo, pur con incertezze maggiori nelle prime decadi. 

La strategia analitica è un’EOF analysis “a finestre” sul Nord Pacifico (105°E–80°W; 20–65°N), calcolando le EOF in intervalli mobili di 30 anni. La scelta della finestra è un compromesso: abbastanza lunga da includere informazione su scala interannuale–interdecadale, ma sufficientemente corta da far emergere eventuali cambiamenti nel tempo della struttura dominante della variabilità (cioè non assumere stazionarietà). L’aspetto più interessante è che le EOF vengono calcolate in due versioni parallele: (1) la versione “canonica” della PDO, in cui dalle anomalie locali di SST viene sottratta l’anomalia della temperatura media globale del mare (GMSST), così da ridurre la contaminazione del trend globale e mettere in evidenza la componente interna; e (2) la versione su anomalie totali, senza rimozione del GMSST, che lascia “entrare” esplicitamente lo spostamento della baseline. Questa distinzione è centrale nella letteratura: la PDO è spesso definita proprio come primo EOF delle SST nel Nord Pacifico dopo aver rimosso il segnale medio globale, mentre studi più recenti sottolineano come la PDO sia un pattern emergente e come la sua interpretazione (e i suoi legami con atmosfera/ecosistemi) possa cambiare nel tempo, soprattutto in un clima che si riscalda. 

Operativamente, le EOF sono ottenute con decomposizione ai valori singolari (SVD) usando pesi areali (tipicamente proporzionali a √cos(lat) per evitare che le alte latitudini “pesino” troppo solo perché hanno più punti griglia) e implementazione in Climate Data Operators (CDO); l’articolo specifica metodo Jacobi per la SVD, e questo dettaglio serve più a documentare la riproducibilità che a introdurre una scelta “esotica”: la sostanza è l’uso di una SVD pesata e robusta su campi grigliati. 

Un passaggio metodologico non banale è il modo in cui vengono definite le anomalie: mensili rispetto alla climatologia estesa 1854–2024 (non una climatologia “standard” tipo 1981–2010). La logica dichiarata è coerente con l’obiettivo: se vuoi descrivere come la variabilità interannuale–interdecadale “si muove” rispetto a una baseline che cambia, centrare tutto su una climatologia lunga evita di “inseguire” il riscaldamento con una climatologia che cambia a seconda del periodo e rende più leggibile lo spostamento del sistema nel tempo (pur accettando che le prime decadi abbiano maggiore incertezza osservativa). Questa scelta si collega direttamente al messaggio dell’articolo: quando il riscaldamento pan-bacino diventa dominante, il primo modo delle anomalie totali può non assomigliare più alla PDO classica, ma a un pattern di riscaldamento quasi unidirezionale sul bacino (il PBP), mentre la variabilità PDO può restare dinamicamente presente ma “secondaria” in termini di varianza spiegata e soprattutto di significato fisico-biologico. 

Per rendere confrontabili i pattern e costruire serie lunghe, gli autori fanno poi un’operazione tipica ma potente: prendono le EOF del periodo più recente (1994–2023), calcolate sia “con GMSST rimosso” sia “totali”, e le proiettanosull’intero archivio 1854–2024 per ottenere due serie di componenti principali (PC) estese nel tempo. La normalizzazione delle PC sulla deviazione standard dell’intero record serve a mettere su una scala comune i punteggi e a poter interpretare regressioni/ampiezze in modo coerente. A quel punto stimano “quanto vale” ciascun modo in termini di SST reali con una regressione delle PC sulle anomalie mensili di SST, ottenendo per ogni cella di griglia una funzione lineare PC→SST (di fatto, una mappa di regressione/transfer function). Questo consente di ricostruire, a ogni passo temporale, il campo di anomalie atteso associato a ciascun modo (PDO canonica vs modo “totale” tipo PBP) e di sintetizzare l’effetto medio su regioni marine chiave (LMEs come GOA, California Current, Kuroshio, Oyashio). In sostanza: non si limitano a dire “l’indice sale/scende”, ma traducono l’indice in un’anomalia fisica media regionale, che è ciò che poi ha più senso quando vuoi parlare di impatti (biologia, pesca, ondate di calore marine, ecc.). 

Due considerazioni “da lettore tecnico” aiutano a inquadrare la solidità (e i limiti) di questo impianto. Primo: i risultati dipendono dall’affidabilità delle SST ricostruite, e ERSSTv5 è stato aggiornato proprio per gestire meglio bias e cambi di strumentazione (ad esempio transizione nave–boa, uso di nuove sorgenti, validazioni e interconfronti), ma l’incertezza rimane più alta nel XIX secolo e nelle regioni/epoche con scarsa copertura. Secondo: l’EOF è uno strumento estremamente utile ma “statistico”: l’ordine dei modi e la loro forma possono cambiare con il periodo, la scelta del dominio, la rimozione o meno del segnale globale e perfino con la metrica/pesatura; in un clima non stazionario, non è un difetto, è spesso proprio il segnale che vuoi misurare, ma richiede cautela nell’interpretazione causale (cioè non confondere “modo principale” con “meccanismo unico”). Su questo, la letteratura moderna sulla PDO insiste: è un pattern emergente da più processi e non un singolo oscillatore, e le sue teleconnessioni possono mutare con il background state. 

Riferimenti essenziali collegati a questi metodi: Huang et al., 2017 (ERSSTv5); Newman et al., 2016 (PDO revisited e definizione/interpretazione); Climate Data Guide UCAR (definizione PDO e rimozione GMSST); Cluett et al., 2025 (PBP e non-stazionarietà recente); documentazione CDO (strumento e operatori).

Questo Extended Data Fig. 1 serve a mettere “nero su bianco” un punto che in climatologia statistica è cruciale: quando il clima cambia di stato (baseline che si scalda), anche i modi dominanti estratti con EOF possono riordinarsi, e un indice storico come la PDO può smettere di coincidere con il primo modo della variabilità delle SST calcolata sulle anomalie “totali”, pur restando una componente dinamicamente robusta del sistema. L’analisi è costruita su finestre mobili di 30 anni nel Nord Pacifico, usando ERSSTv5, un dataset ricostruito globale mensile a 2°×2° progettato proprio per descrivere trend e variabilità su scala secolare con un impianto omogeneo (pur con incertezze maggiori nelle prime decadi del record). 

Il pannello A confronta la somiglianza spaziale (pattern correlation arealmente pesata) tra le EOF ottenute in ciascuna finestra. La curva con i punti neri misura quanto l’EOF1 calcolata sulle anomalie totali di SST assomigli all’EOF1 calcolata dopo aver rimosso la GMSST (cioè la costruzione “canonica” della PDO). La curva con i punti bianchi misura invece quanto l’EOF2 delle anomalie totali assomigli alla stessa EOF1 “PDO” (GMSST-rimossa). Finché i punti neri stanno prossimi a 1, significa che il primo modo delle anomalie totali è praticamente la PDO classica: cioè la variabilità interna legata all’accoppiamento oceano–atmosfera nordpacifico (Aleutian Low, flussi aria–mare, teleconnessioni, onde di Rossby oceaniche) domina ancora il primo autovettore statistico, in linea con la visione “PDO come pattern emergente” descritta dalla letteratura moderna. Nelle finestre che terminano tra fine anni ’90 e primi 2000, infatti, la correlazione nera è quasi perfetta, mentre quella bianca è bassa: la PDO è “EOF1”, non “EOF2”. Poi però, entrando nelle finestre che terminano negli anni 2010, accade la cosa importante: la correlazione nera inizia a degradare e, quasi specularmente, la correlazione bianca cresce fino a valori molto alti. Il messaggio è: la mappa tipica della PDO non sparisce, ma migra dal primo al secondo modo delle anomalie totali. In pratica l’EOF1 delle SST totali non è più “la PDO”, bensì un pattern più coerente con un riscaldamento esteso sul bacino (il PBP) che, proprio perché è ampio e persistente, finisce per spiegare più varianza e quindi “si prende” il primo posto nella decomposizione. Questo è coerente con la tesi centrale di Cluett et al.: la modalità dominante recente delle SST nordpacifiche (sulle anomalie totali) è passata dalla PDO a un pattern pan-bacino, e ciò segnala una crescente non-stazionarietà delle relazioni PDO–SST regionali. 

Il pannello B ripete lo stesso test, ma non sulle mappe: lo fa sulle serie temporali (principal components, PC). Qui i punti neri sono la correlazione tra PC1 “totale” e PC1 “PDO” (GMSST rimossa), mentre i punti bianchi sono la correlazione tra PC2 “totale” e PC1 “PDO”. Anche qui la lettura è molto istruttiva: per molte finestre storiche la PDO è ben rappresentata dal primo componente principale delle anomalie totali (correlazione nera alta), quindi se tu avessi usato PC1-totale come “indice del Nord Pacifico” avresti, di fatto, seguito la PDO. Con l’avanzare degli anni 2010, però, la correlazione tra PC1-totale e PC1-PDO cala, mentre cresce quella tra PC2-totale e PC1-PDO: la variabilità PDO continua a essere presente e riconoscibile, ma viene “spinta” sul secondo componente principale perché il primo viene occupato da un segnale più grande, più coerente e più persistente, legato allo spostamento della baseline termica del bacino. È un risultato metodologicamente pulito perché mostra che non stiamo discutendo un’impressione qualitativa (“la PDO non funziona più”), ma un riassetto misurabile tra i primi due autospazi della variabilità.

Questa dinamica si collega bene anche a ciò che sappiamo sul sistema Pacifico–Alaska. Storicamente, la PDO è stata resa celebre proprio per i suoi legami con la produttività e le riorganizzazioni ecosistemiche (salmoni inclusi), e per l’idea di “regimi” o “shift” (1976/77 e 1988/89) che riorganizzano atmosfera, oceano e biologia su scale decadali. Il punto è che in un contesto recente caratterizzato da riscaldamenti molto marcati e persistenti nel Nord-Est Pacifico (es. l’anomalia calda 2013–2014 e l’evento associato al “Blob”), l’oceano può entrare in combinazioni di condizioni senza analoghi stretti nel record recente, e questo altera sia la relazione tra indice e variabili fisiche “primarie”, sia la relazione tra indice e risposta biologica. Non a caso, proprio nel Golfo dell’Alaska, durante il riscaldamento estremo 2014–2019 è stata documentata una transizione verso correlazioni PDO–clima e PDO–salmone più deboli o addirittura di segno opposto, interpretata come segnale di “novel climate”. In quel quadro, l’Extended Data Fig. 1 è quasi la cerniera matematica che rende coerente tutto: se il primo modo delle anomalie totali diventa un pattern pan-bacino legato al riscaldamento, allora è normale che molte vecchie scorciatoie basate su “PDO = primo modo = principale driver” perdano potere esplicativo; la PDO resta utile, ma va letta per quello che è (variabilità interna emergente) e soprattutto va affiancata da una descrizione esplicita dello stato di fondo, perché è lo stato di fondo che determina se una stessa oscillazione ti porta verso l’optimum o oltre l’optimum (come poi formalizza la figura concettuale sulle relazioni PDO–biologia). 

Questa Extended Data Fig. 2 è, di fatto, un “test di robustezza” stagionale: gli autori rifanno la stessa diagnostica a finestre mobili (30 anni) delle EOF delle anomalie di SST nel Nord Pacifico, ma usando solo l’inverno boreale NDJFM (novembre–marzo), cioè la stagione in cui il coupling oceano–atmosfera nel Pacifico settentrionale è più intenso (mixed layer più profondo, flussi turbolenti aria–mare più forti, risposta più diretta alla variabilità dell’Aleutian Low). In ogni riga vedi una finestra diversa (dalle più “vecchie” in alto alle più recenti in basso), mentre le colonne separano le EOF calcolate sulle anomalie totali di SST (“Total SSTa”, quindi includendo anche lo spostamento della baseline) da quelle calcolate dopo aver rimosso la GMSST (“GMSST removed”, cioè l’impostazione canonica usata per isolare la PDO). Ogni pannello è una mappa di “loading” dell’EOF (il segno è convenzionale: ciò che conta è la geometria del pattern e come si riorganizza nel tempo). 

Quello che salta all’occhio è che nelle finestre più datate l’EOF1 delle anomalie totali invernali assomiglia ancora molto alla PDO classica: un’impronta a larga scala con massimi nel Nord-Est Pacifico/Golfo dell’Alaska e contrasti con il Pacifico centrale, coerente con l’idea storica della PDO come pattern emergente della variabilità interna nordpacifica legata a circolazione atmosferica e feedback aria–mare. Scendendo verso le finestre che includono sempre più anni recenti, però, succede la transizione chiave: l’EOF1 calcolata sulle anomalie totali diventa progressivamente più “pan-bacino”, cioè tende a rappresentare un riscaldamento ampio e coerente su gran parte del Nord Pacifico, mentre l’impronta PDO-like si ritrova sempre più chiaramente nell’EOF2. È esattamente il riordino diagnostico mostrato anche dalle correlazioni dell’Extended Data Fig. 1: la PDO non sparisce, ma viene “retrocesso” a modo secondario quando nel calcolo lasci entrare il riscaldamento di fondo, perché quel segnale è spazialmente esteso, persistente e quindi capace di spiegare una quota crescente di varianza come primo autovettore. 

La parte più importante, qui, è che questo riordino si vede anche in inverno, e questo rafforza molto l’interpretazione fisica: non è un effetto confinato all’estate o a stagioni in cui la stratificazione e il forcing radiativo superficiale potrebbero “mascherare” la dinamica; emerge proprio nel semestre freddo, quando la SST risponde con maggiore efficienza alle anomalie di vento, mixing e flussi di calore associati alla variabilità del bacino e quando la PDO tende a mostrare forte persistenza/“reemergence” tra un inverno e l’altro. In altre parole, se perfino nella stagione in cui la PDO ha storicamente la massima impronta atmosferica ed oceanica l’EOF1 delle anomalie totali viene occupata da un pattern di riscaldamento diffuso, allora l’emergere del pan-basin warming non è un artefatto di stagionalità: è il segnale che la baseline termica del Pacifico nordico si è spostata a un punto tale da dominare la gerarchia dei modi della variabilità quando non la si rimuove esplicitamente. 

Il confronto con le mappe “GMSST removed” serve proprio a chiudere il cerchio: lì l’EOF1 rimane molto più stabile e riconducibile alla PDO lungo le finestre, perché la sottrazione del segnale medio globale impedisce al trend pan-bacino di “prendersi” il primo modo; questo è coerente con l’uso storico della PDO come indicatore di variabilità interna (e con il fatto che la PDO, come pattern, può rimanere dinamicamente robusta anche mentre cambia lo stato di fondo). Sullo sfondo c’è anche un messaggio operativo: negli ultimi decenni il Nord-Est Pacifico ha sperimentato episodi estremi di riscaldamento persistente (il caso emblematico è l’anomalia 2013–2014 e la sequenza associata al “Blob”), che mostrano come la combinazione tra forcing atmosferico stagionale e riscaldamento di lungo periodo possa produrre stati termici “nuovi” rispetto al passato recente; in un mondo così, aspettarsi che un indice costruito per descrivere la variabilità interna mantenga invariata la sua relazione con SST assolute, ecosistemi o teleconnessioni è sempre meno giustificato. 

La Extended Data Fig. 3 risponde a una domanda molto concreta: l’emergere del pan-basin warming (PBP) come modo dominante delle SST del Nord Pacifico dipende dal fatto che le EOF siano state calcolate su finestre “lunghe” (30 anni), oppure lo vedi comunque anche se riduci la finestra a 20 anni, quindi con un campionamento più breve e potenzialmente più “rumoroso”? L’impianto è lo stesso della figura principale: ogni riga è una finestra mobile (qui 20 anni; da 1988–2007 fino a 2004–2023), le due colonne a sinistra sono le EOF calcolate sulle anomalie totali di SST(“Total SSTa”, quindi includendo il trend/shift della baseline), mentre le due colonne a destra sono le EOF calcolate dopo aver rimosso la GMSST (“GMSST removed”), cioè la procedura tipica con cui si isola la PDO canonica come variabilità interna del Nord Pacifico. 

Il primo messaggio della figura è metodologico: passando a 20 anni, ti aspetti una maggiore sensibilità delle EOF a fluttuazioni campionarie, perché le EOF stimate da campioni finiti possono “oscillare” tra loro quando gli autovalori sono vicini (il classico problema della separazione tra EOF1 ed EOF2 e delle incertezze di campionamento; la regola pratica di North e lavori successivi servono proprio a capire quando due modi non sono ben separati e possono mescolarsi/ruotare tra finestre). E infatti, guardando le righe intermedie, si nota una transizione più nervosa e meno “liscia” rispetto alle finestre trentennali: è normale, perché in 20 anni hai meno realizzazioni indipendenti e la variabilità decadale stessa (che per definizione vive su scale comparabili o superiori) è più difficile da campionare stabilmente.

Detto questo, il secondo messaggio — quello scientificamente pesante — è che la storia non cambia: anche con finestre ventennali, nelle mappe “Total SSTa” l’EOF1 tende progressivamente a perdere l’impronta PDO-like e a diventare un pattern più coerente e diffuso di segno caldo su gran parte del bacino, mentre un pattern più simile alla PDO viene recuperato sempre più spesso nell’EOF2. In altre parole, il “riordino” tra i primi due modi non è un artefatto della finestra a 30 anni: lo vedi anche quando accetti un maggiore rumore statistico. Questa è esattamente l’idea del lavoro di Cluett et al.: nelle anomalie totali, il modo che spiega più varianza nel periodo recente non è più la PDO classica, ma un riscaldamento quasi unidirezionale a scala di bacino (PBP), mentre la PDO resta presente come variabilità interna robusta ma “scende di rango” nella decomposizione. 

Il confronto con “GMSST removed” serve a isolare la dinamica interna e a dimostrare che non stai “perdendo” la PDO: nelle colonne di destra l’EOF1 rimane molto più stabile e riconoscibile come PDO nel senso canonico (dominante pattern di variabilità mensile del Nord Pacifico una volta rimosso il contributo medio globale), coerentemente con la letteratura che interpreta la PDO come pattern emergente dall’integrazione di più processi (forcing aleutinico stocastico, teleconnessioni tropicali, onde di Rossby oceaniche) più che come un singolo oscillatore deterministico. Il punto quindi non è “la PDO è cambiata” (nel senso di scomparsa del pattern), ma che quando non togli il segnale globale, il riscaldamento di fondo diventa abbastanza grande e spazialmente coerente da dominare la varianza e conquistare l’EOF1.

C’è anche un’interpretazione fisica coerente con gli ultimi decenni osservati: il Nord-Est Pacifico e il Golfo dell’Alaska hanno sperimentato episodi di riscaldamento eccezionalmente persistenti (il caso scuola è la warm anomaly 2013–2014 e la sequenza della marine heatwave “Blob”), innescati e mantenuti da anomalie atmosferiche che riducono il raffreddamento invernale e modificano mixing e flussi aria–mare; quando questi episodi si sommano a una baseline che sale, diventa più facile ottenere campi SST con un segno caldo molto esteso, cioè esattamente ciò che un’EOF tende a catturare come primo modo nelle anomalie totali. E non è un dettaglio “solo fisico”: nello stesso periodo si osserva una non-stazionarietà delle relazioni tra PDO e variabili climatiche/ecosistemiche nel GOA, con correlazioni che si indeboliscono o si invertono durante il riscaldamento 2014–2019, proprio perché l’indice PDO non sta più campionando lo stesso spazio di stati oceanici/atmosferici di alcune decadi fa. 

In sintesi, questa figura dice: anche scegliendo una finestra più corta (e quindi più esposta a incertezze EOF), la transizione resta visibile e coerente — l’EOF1 delle anomalie totali tende a rappresentare il PBP nel periodo recente, mentre la PDO rimane più chiaramente rappresentata quando rimuovi la GMSST (o, nelle anomalie totali, quando guardi l’EOF2). Questo rafforza molto la conclusione che il risultato sia strutturale (cambiamento dello stato di fondo e della gerarchia dei modi), non una scelta arbitraria di timeframe. 

La Extended Data Fig. 4 ripete l’analisi delle EOF delle SST nel Nord Pacifico usando finestre mobili di 50 anni, quindi un campionamento molto più lungo rispetto alle versioni a 20 o 30 anni. Il motivo è semplice: se un risultato appare anche su scale cinquantennali, diventa difficile liquidarlo come fluttuazione campionaria o come “fase” temporanea della variabilità interna. Le mappe mostrano, per ciascuna finestra (righe, dal 1958–2007 fino al 1974–2023), i primi due modi (EOF1 ed EOF2) calcolati su anomalie totali di SST (“Total SSTa”, a sinistra) e su anomalie con GMSST rimossa (“GMSST removed”, a destra), che corrisponde alla costruzione canonica usata per isolare la PDO come variabilità interna del bacino. 

Nelle finestre più “antiche” (quelle che terminano prima della piena accelerazione del riscaldamento recente), l’EOF1 di Total SSTa conserva una struttura ampiamente riconducibile alla classica impronta PDO-like: un pattern a larga scala con massimi nel Nord-Est Pacifico/Golfo dell’Alaska e contrasti sul Pacifico centrale, cioè la firma statistica di una variabilità che integra forzante atmosferica (Aleutian Low e flussi aria-mare) e risposta oceanica. Questo è coerente con l’interpretazione moderna della PDO come pattern emergente (non un oscillatore “a sé”), risultante dalla combinazione di processi stocastici e dinamici che proiettano la loro varianza sulle SST del Nord Pacifico. 

Scendendo lungo le righe, quindi includendo un numero crescente di anni recenti nelle finestre cinquantennali, l’EOF1 di Total SSTa cambia natura: la mappa tende a diventare sempre più “monosimmetrica”, con loadings di segno caldo diffusi su gran parte del bacino, cioè un pattern compatibile con un riscaldamento pan-bacino (PBP) che domina la varianza quando non si rimuove esplicitamente il segnale globale. In parallelo, l’impronta PDO non scompare, ma si ritrova più chiaramente nell’EOF2 di Total SSTa. È esattamente la gerarchia descritta da Cluett et al.: la PDO resta robusta come variabilità interna, ma viene “sorpassata” dal riscaldamento di fondo nella competizione per il primo autovettore quando si analizzano le anomalie totali. 

Il confronto con le colonne GMSST removed è la chiave interpretativa: lì l’EOF1 rimane molto più stabile e persistentemente PDO-like lungo tutte le finestre di 50 anni. Questo dice che non stai osservando una “morte” della PDO né un collasso della sua fisica; stai osservando un cambio di regime statistico dovuto al fatto che, lasciando il trend/shift della baseline dentro il campo SST, quel segnale diventa abbastanza coerente e persistente da essere identificato come modo dominante dalla decomposizione. In altre parole, la variabilità interna rimane, ma cambia il peso relativo rispetto al segnale di riscaldamento del bacino. 

Il fatto che questo riordino si veda anche con finestre di 50 anni è particolarmente significativo dal punto di vista metodologico, perché le EOF sono note per essere sensibili al campionamento quando i modi sono poco separati; allungare la finestra aumenta il numero di realizzazioni effettive e riduce le fluttuazioni campionarie attese, secondo criteri classici sulle incertezze di stima delle EOF. Quindi la comparsa del PBP come EOF1 nelle finestre più recenti non è solo “rumore statistico”, ma un segnale robusto di non-stazionarietà della baseline termica del Nord Pacifico. 

Un ulteriore elemento di coerenza fisica è che la fase recente include episodi di riscaldamento persistente nel Nord-Est Pacifico (marine heatwaves) che, sovrapposti al trend, aumentano la probabilità di campi SST con segno caldo esteso e durata pluristagionale: l’evento 2013–2014 (il “Blob”) è un riferimento classico proprio perché mostra come anomalie atmosferiche e ridotto raffreddamento/mixing possano mantenere SST eccezionalmente elevate per mesi, e in un oceano già più caldo questi stati diventano più “competitivi” nel determinare la varianza totale. 

Infine, questa figura aiuta a capire perché molte relazioni storiche PDO-impatto (biologia, pesca, teleconnessioni regionali) possano degradarsi o invertire segno: se nelle anomalie totali la PDO non è più il primo modo, allora l’indice PDO “vede” una porzione diversa dello spazio degli stati climatici rispetto al passato, mentre la variabile che oggi domina la SST regionale è spesso la combinazione tra variabilità interna e riscaldamento pan-bacino. Questo quadro è compatibile con la letteratura sulle relazioni non stazionarie PDO-clima e PDO-salmoni nel Golfo dell’Alaska durante il riscaldamento estremo 2014–2019.

Questa Extended Data Fig. 5 è, in pratica, la “prova di resistenza” più severa dell’intero ragionamento: rifare l’analisi EOF con finestre mobili di 100 anni significa costringere la statistica a lavorare su campioni enormi, dove la variabilità interannuale e decadale viene fortemente mediata e dove, in teoria, i pattern davvero robusti dovrebbero emergere con chiarezza. Ogni riga rappresenta una finestra secolare che scorre in avanti (1908–2007, 1910–2009 … fino a 1924–2023). In ciascuna finestra gli autori calcolano i primi due modi (EOF1 e EOF2) in due modi diversi: a sinistra sulle anomalie totali di SST (“Total SSTa”, quindi includendo lo spostamento della baseline legato al riscaldamento globale), a destra sulle anomalie con GMSST rimossa (“GMSST removed”), cioè l’impostazione tipica con cui la PDO viene isolata come variabilità interna del Nord Pacifico (la definizione “canonica” è proprio: EOF1 delle SST nel Nord Pacifico dopo aver sottratto il segnale medio globale). Questo passaggio è in linea con la costruzione classica della PDO e con la letteratura che la interpreta come pattern emergente dall’integrazione di più processi (forzante aleutinica, feedback aria–mare, teleconnessioni tropicali, onde di Rossby oceaniche), non come un oscillatore unico. Qui i colori sono “loadings” dell’EOF: non sono gradi Celsius, ma coefficienti che dicono dove il modo ha massimi/minimi e quindi qual è la sua geometria.

La lettura fisica è: nelle finestre più “vecchie”, quando la parte recente e molto calda del record pesa meno nel campione, l’EOF1 delle anomalie totali mantiene un’impronta chiaramente riconoscibile come PDO-like: un centro d’azione marcato nel Nord-Est Pacifico/Golfo dell’Alaska e un contrasto su scala di bacino con il Pacifico centrale/occidentale. È esattamente il tipo di struttura che la comunità associa storicamente alla PDO (Mantua et al., 1997; Zhang et al., 1997) e che in molte applicazioni è stata usata come “bussola” per interpretare riorganizzazioni del Pacifico nordico e impatti ecosistemici. Ma scendendo lungo le righe, quindi includendo sempre più anni recenti nella finestra secolare, accade la transizione che lo studio vuole dimostrare: l’EOF1 delle anomalie totali diventa progressivamente più uniforme e “pan-bacino”, con segno caldo diffuso su larga parte del Nord Pacifico. In parallelo, la firma PDO non sparisce: tende a essere catturata sempre più spesso dall’EOF2 nel blocco “Total SSTa”. Questo è un punto chiave perché ti dice che la PDO continua a esistere come variabilità interna del sistema, ma non è più il modo che spiega la maggiore quota di varianza quando guardi le SST “reali” (totali); a dominare è lo spostamento della baseline termica, cioè quel contributo coerente e persistente che oggi è inevitabile in un oceano in riscaldamento (IPCC AR6; e, sul lato osservativo, dataset SST ricostruiti come ERSSTv5 descrivono chiaramente trend di lungo periodo).

Il confronto con la metà destra (“GMSST removed”) è ciò che rende l’interpretazione davvero pulita: lì l’EOF1 rimane molto più stabile e assomiglia in modo consistente alla PDO lungo tutte le finestre secolari. Questo implica che non stai assistendo a una “morte” dinamica della PDO; stai vedendo un riassetto gerarchico indotto dal fatto che, se non rimuovi il segnale globale, un pattern di riscaldamento ampio e coerente finisce per occupare il primo autovettore della decomposizione. È un comportamento atteso nelle EOF quando introduci nel campo una componente a grande scala, persistente e con forte proiezione spaziale: la SVD tende a catturare quella direzione come massima varianza, relegando la variabilità interna (che resta) a modi successivi. Questa logica è esattamente ciò che Newman et al. (2016) sottolineano quando discutono la PDO come pattern statistico emergente e potenzialmente sensibile al “background state”, e più in generale è coerente con l’idea di non-stazionarietà delle teleconnessioni in un clima che cambia.

Il fatto che la transizione si veda anche con finestre di 100 anni è metodologicamente rilevante per un altro motivo: le EOF, quando i primi autovalori sono relativamente vicini, possono “mescolarsi” o ruotare tra loro a seconda del campione (il classico problema di separazione dei modi e incertezza campionaria discusso, in forme diverse, già da North et al., 1982 e nella letteratura successiva sull’affidabilità delle EOF). Se un risultato fosse un artefatto di finestra, te lo aspetteresti fragile e incoerente quando allarghi il campione; qui invece l’emergere di un primo modo pan-bacino nelle anomalie totali diventa semmai più credibile, perché stai mostrando che non è una particolare fase decadale “presa al volo”, ma un cambiamento che lascia traccia persino quando medi su un secolo.

In termini concettuali, questa figura dice: nel XX secolo la PDO poteva apparire come “la” modalità dominante del Nord Pacifico perché il segnale di riscaldamento di fondo era più piccolo rispetto alla varianza interna e/o meno determinante nella gerarchia dei modi; nel clima recente, l’aumento termico persistente su scala di bacino diventa abbastanza grande da ridefinire lo spazio delle anomalie, e la PDO, pur restando fisicamente reale, viene statisticamente superata. È esattamente il passaggio che poi rende sensato tutto il discorso sugli impatti: se il principale modo delle SST totali non è più PDO-like, allora non è sorprendente che molte relazioni storiche PDO-ecosistemi o PDO-clima regionale mostrino decorrelazioni o inversioni di segno, perché l’indice PDO sta operando sopra una baseline spostata, e l’ecosistema risponde spesso alla temperatura assoluta e alle combinazioni di forzanti più che a un singolo indice (un quadro compatibile con i lavori su non-stazionarietà nel Golfo dell’Alaska durante il riscaldamento recente, e con la letteratura sulle marine heatwaves del Nord-Est Pacifico come amplificatori del segnale di fondo).

La Extended Data Fig. 6 è costruita per far vedere, in un colpo solo, la differenza tra (i) il segnale globale di riscaldamento che sposta la baseline termica degli oceani e (ii) la variabilità interna del Nord Pacifico descritta dagli indici “classici” (PDO, NPGO), e poi collegare tutto alle componenti principali (PC) ottenute dalle EOF calcolate nel periodo 1994–2023 e proiettate sul record lungo. Il pannello (a) mostra l’anomalia di GMSST derivata da NOAA ERSST: la tendenza positiva che accelera nel tardo XX–inizio XXI secolo è l’impronta del forcing di lungo periodo, ed è proprio la componente che, se lasciata “dentro” alle SST regionali, tende a dominare la varianza e quindi a emergere come modo principale nelle decomposizioni EOF. ERSSTv5 è un dataset ricostruito globale mensile a 2°×2° sviluppato per rappresentare in modo robusto trend e pattern di variabilità su lunghi periodi, con procedure di quality control, bias adjustment e ricostruzione statistica aggiornate rispetto alle versioni precedenti. 

Nel pannello (b) compare la serie dell’indice PDO fornito da NOAA PSL: è un segnale oscillatorio a scala interannuale–interdecadale, con alternanza di fasi positive e negative e senza un drift monotono paragonabile alla GMSST. Questo è coerente con la definizione operativa della PDO: il NOAA PSL la costruisce come EOF1 delle anomalie mensili di SST nel Nord Pacifico (a nord di ~20°N) dopo aver rimosso la climatologia mensile e, crucialmente, anche la SST media globale (proprio per ridurre la contaminazione del trend globale e isolare la variabilità interna). In altre parole, la PDO è pensata per descrivere un pattern “ENSO-like” a bassa frequenza nel Pacifico settentrionale che emerge dall’accoppiamento oceano–atmosfera e da teleconnessioni, e infatti la sua interpretazione storica come indice di regime e impatti ecosistemici risale ai lavori fondativi di fine anni ’90. 

Il pannello (c) aggiunge l’indice NPGO (NOAA PSL), introdotto in letteratura come un modo distinto della variabilità nordpacifica, spesso descritto come espressione oceanica (e in parte “giroscopica”) legata alla seconda EOF di altezza del livello del mare/sea surface height e con connessioni a trasporti, nutrienti e risposte ecosistemiche lungo i sistemi di correnti di margine. L’idea di includerlo qui è concettualmente importante: il Nord Pacifico non ha “un solo” modo; PDO e NPGO rappresentano famiglie diverse di variabilità e di meccanismi (atmosferici e oceanici) che possono proiettarsi in modo differente sulle SST e sugli impatti. 

La parte davvero diagnostica della figura sta nei pannelli (d–g), dove compaiono le componenti principali associate alle EOF calcolate nel periodo 1994–2023 e poi usate come “pattern di riferimento” proiettati sul record lungo. In (d) la PC1 delle anomalie totali di SST (Total SSTa PC1) è etichettata come PBP: la sua evoluzione è marcatamente non stazionaria e tende a crescere nettamente verso valori positivi nel periodo recente, comportamento compatibile con un segnale di riscaldamento diffuso a scala di bacino che “si prende” il primo modo quando non sottrai il contributo globale. Questo è esattamente il risultato centrale di Cluett et al.: nelle anomalie totali, il modo principale delle SST nordpacifiche è recentemente transitato da un’impronta PDO-like a un pattern di cambiamento quasi unidirezionale su tutto il bacino (pan-basin pattern, PBP), mentre la variabilità PDO rimane presente ma non domina più la varianza. 

Il pannello (e) completa la lettura: la PC2 delle anomalie totali (Total SSTa PC2) è etichettata come PDO e torna a comportarsi come un segnale oscillatorio attorno allo zero, molto più simile, per struttura temporale, alla serie dell’indice PDO mostrata in (b). Il messaggio è tecnico ma cruciale: quando lasci entrare la baseline globale, l’EOF1/PC1 tende a catturare la componente “di fondo” (PBP), e la PDO viene spesso ritrovata come secondo modo (EOF2/PC2) nel campo totale. È la stessa “migrazione gerarchica” che spiega perché, negli ultimi decenni, una PDO negativa possa non tradursi automaticamente in SST orientali fredde come accadeva storicamente: il pavimento termico si è alzato e la stessa variabilità interna opera su un livello medio diverso. 

Infine, (f) e (g) mostrano le PC quando le EOF sono calcolate su SST con GMSST rimossa. Qui il sistema “torna” a una descrizione dominata dalla variabilità interna: la PC1 GMSST-removed (g) è un segnale più simmetrico e oscillatorio, coerente con l’idea che, sottraendo la componente media globale, il primo modo recuperi la struttura PDO-like per costruzione e per fisica. La PC2 GMSST-removed (f) rappresenta ulteriore variabilità interna non riducibile alla PDO (e, a seconda delle finestre e dell’assetto spettrale, può includere contributi affini ad altri modi del Nord Pacifico), ma la cosa essenziale è che non mostra più la “deriva” tipica del PBP. In sintesi, la figura formalizza un punto operativo: parlare di “PDO che non funziona più” è fuorviante; ciò che cambia è il fatto che, nelle SST osservate totali, il riscaldamento pan-bacino può dominare il primo modo e spostare la PDO su un modo successivo, rendendo non stazionarie molte vecchie correlazioni indice–impatto se non si tiene conto della baseline termica.

Questa Extended Data Fig. 7 è pensata come un controllo “diagnostico” molto trasparente: prende le serie temporali mostrate nella Fig. 6 (PDO NOAA, GMSST e le componenti principali ottenute dalle EOF calcolate nel 1994–2023 e poi proiettate sul record lungo) e le mette a confronto con scatterplot mensili, in modo da chiarire quale componente sta catturando il riscaldamento di fondo e quale sta catturando la variabilità interna tipo PDO. Il fatto che i punti siano colorati in scala di grigi per il tempo (più chiaro = anni antichi, più scuro = anni recenti) serve a visualizzare immediatamente la non-stazionarietà: se un legame è dominato dal trend recente, i punti più scuri tendono a “spostarsi” lungo la retta; se invece il legame è interno/oscillatorio, la nube resta distribuita senza deriva sistematica. L’impianto è coerente con la definizione operativa della PDO usata da NOAA PSL (EOF1 delle SSTA nel Nord Pacifico con rimozione della media globale), proprio per separare la variabilità interna dal riscaldamento globale. 

I tre pannelli della riga superiore mettono un punto fermo: la PDO “fisica” è la stessa variabile, ma cambia “numero di PC” a seconda che tu rimuova o meno la GMSST. In alto a sinistra, PDO (asse x) contro GMSST-removed PC1 (asse y) dà R = 0.813, mostrando che quando togli il segnale globale la prima componente principale recupera molto bene l’indice PDO ufficiale. In alto al centro, PDO contro Total SSTa PC2 dà R = 0.792, cioè la stessa coerenza ma spostata sulla seconda componente quando lavori sulle anomalie totali (senza rimozione della GMSST). Il pannello in alto a destra chiude la dimostrazione: GMSST-removed PC1 contro Total SSTa PC2 ha R = 0.997, praticamente identità lineare, che è esattamente ciò che ti aspetti se la PDO non “sparisce” ma viene semplicemente retro-classificata da PC1 a PC2 quando nel campo SST lasci entrare un segnale a grande scala e persistente (il riscaldamento pan-bacino) capace di dominare la varianza. Questa lettura è perfettamente allineata sia con la definizione NOAA della PDO sia con la letteratura che la descrive come pattern emergente dalla dinamica oceano-atmosfera nordpacifica, non come un oscillatore isolato. 

La riga inferiore serve invece a identificare “chi è” la componente di riscaldamento. In basso a sinistra, GMSST (asse x) contro Total SSTa PC1 (asse y) mostra R = 0.801: una correlazione forte, con una nube che tende a organizzarsi lungo una relazione positiva e con i punti più scuri (periodo recente) che occupano più frequentemente la porzione alta della distribuzione. Questo pannello giustifica l’etichetta PBP/riscaldamento pan-bacino per la PC1 delle anomalie totali: non sta descrivendo un’oscillazione simmetrica attorno allo zero, ma una componente che segue da vicino lo spostamento della baseline termica globale. È un risultato coerente con il fatto che dataset come ERSSTv5 sono costruiti per rappresentare robustamente trend e variabilità su scale secolari, e con la conclusione del lavoro secondo cui, nel periodo recente, il primo modo delle SST totali nel Nord Pacifico tende a essere dominato da un segnale di riscaldamento diffuso. 

Gli ultimi due pannelli quantificano la separazione tra forzante di fondo e variabilità interna. In basso al centro, GMSST contro Total SSTa PC2 dà R = 0.150 (relazione debole): è esattamente ciò che ci si aspetta se PC2-total è PDO-like, cioè variabilità interna che oscilla senza essere “ancorata” linearmente al trend globale. In basso a destra, Total SSTa PC1 contro Total SSTa PC2 dà R = 0.127, quindi quasi indipendenza: statisticamente è un segnale che l’analisi sta separando bene due dimensioni diverse del sistema, una legata al riscaldamento di fondo e una legata alla variabilità interna del Nord Pacifico. Questa distinzione è concettualmente molto importante perché chiarisce anche un punto spesso sottovalutato nella diagnostica della variabilità decadale: i trend forzati possono proiettarsi su pattern simili a quelli della variabilità naturale e creare “cambiamenti apparenti” se non vengono rimossi o trattati correttamente; qui, invece, la bassa correlazione tra GMSST e PC2-total indica che la componente PDO-like rimane principalmente interna. 

Nel complesso, la figura formalizza con numeri (e non solo con mappe) una regola pratica che vale ben oltre questo studio: quando analizzi le SST del Nord Pacifico in un clima non stazionario, devi distinguere tra (i) la variabilità interna tipo PDO e (ii) lo spostamento della baseline globale. Qui la separazione emerge in modo pulito: Total PC1 ≈ GMSST (componente di fondo/PBP) e Total PC2 ≈ PDO ≈ GMSST-removed PC1 (variabilità interna). Il passaggio da PDO-come-PC1 a PDO-come-PC2 non è quindi una contraddizione, ma la firma statistica del fatto che il riscaldamento pan-bacino è diventato abbastanza grande da “occupare” il primo modo quando non lo si rimuove esplicitamente. 

Questa Extended Data Fig. 8 è, sostanzialmente, la versione “energetica” della storia: invece di mostrarti quantociascun modo spiega in termini percentuali (che è una misura relativa), ti mostra gli autovalori (eigenvalues) delle EOF calcolate su finestre mobili di 30 anni. Un autovalore, in questo contesto, è la varianza assoluta catturata da un modo (EOF/PC): se la varianza totale del campo SST aumenta perché entra un segnale spazialmente coerente e persistente (come il riscaldamento di fondo), allora può aumentare anche l’eigenvalue del modo che lo rappresenta, persino senza che la “forma” del pattern cambi drasticamente. È proprio per questo che gli autori lo usano qui: per dimostrare che l’emergere del pan-basin warming non è solo un riordino “geometrico” dei modi, ma anche un cambiamento nella quota assoluta di varianza associata al modo dominante. 

Nel pannello superiore (“Total SSTa”), cioè EOF calcolate sulle anomalie totali di SST (senza rimuovere la temperatura media globale del mare), l’EOF1 all’inizio delle finestre appare associata alla PDO (barre blu, autovalori relativamente stabili). Poi, nelle finestre che terminano negli anni più recenti, l’eigenvalue dell’EOF1 cresce in modo marcato e il modo viene etichettato come PBP (barre rosse): questo è il segnale quantitativo che il primo modo sta assorbendo una componente di varianza sempre più grande legata a un riscaldamento diffuso e coerente sul bacino. In parallelo, l’EOF2 nel caso “Total SSTa” tende a mostrare la PDO come secondo modo, coerentemente con l’idea centrale dell’articolo: quando non sottrai la GMSST, il riscaldamento pan-bacino diventa abbastanza dominante da occupare l’EOF1, mentre la variabilità PDO-like rimane robusta ma “scende” a EOF2. 

Il pannello inferiore (“GMSST removed”) è il controtest fisico: qui le EOF sono calcolate dopo aver sottratto la GMSST, cioè seguendo la costruzione canonica usata proprio per far sì che la PDO descriva la variabilità interna e non il segnale di global warming. In questa configurazione l’eigenvalue dell’EOF1 (PDO) rimane molto più stabilelungo le finestre, senza l’impennata vista sopra: è un indizio forte che l’esplosione dell’eigenvalue in “Total SSTa EOF1” non è una proprietà intrinseca della PDO, ma nasce dal fatto che il campo SST totale contiene un contributo a grande scala (il riscaldamento di fondo) che aumenta la varianza e si proietta con efficacia sul primo autovettore. Questo è perfettamente coerente con la definizione NOAA/PSL della PDO, che specifica esplicitamente la rimozione della media globale per separare questo pattern interno da un segnale di riscaldamento globale. 

Sempre nel blocco “GMSST removed”, l’EOF2 viene associata più spesso a un modo distinto (qui etichettato NPGO): anche questo è coerente con la letteratura che descrive l’NPGO come una modalità separata della variabilità nordpacifica (originariamente definita come secondo modo dominante della variabilità di sea surface height e collegata a dinamica dei gyre e dei trasporti), che può emergere come secondo EOF quando l’EOF1 è occupata dalla PDO interna. 

C’è poi un aspetto metodologico importante che questa figura “suggerisce” senza dirlo esplicitamente: nelle EOF a finestra mobile, quando due modi hanno autovalori vicini, possono mescolarsi o scambiarsi ordine tra una finestra e la successiva (problema classico di separazione dei modi e incertezza campionaria). Per questo, in un’analisi rigorosa, non basta guardare le mappe: serve verificare la continuità tramite correlazioni tra pattern e PC, e serve leggere gli autovalori proprio come qui, per capire quando un modo diventa davvero dominante. Questo principio è formalizzato nei lavori classici sugli errori di campionamento delle EOF (criterio di separazione degli autovalori). 

In sintesi, Extended Data Fig. 8 ci dice due cose insieme: (1) la PDO resta una variabilità interna robusta (lo vedi nella stabilità degli autovalori quando rimuovi la GMSST), ma (2) nelle SST totali la varianza associata a un modo di riscaldamento pan-bacino cresce fino a dominare l’EOF1 nelle finestre più recenti, rendendo inevitabile che molte relazioni “storiche” basate sulla PDO diventino non stazionarie se non si esplicita lo stato termico di fondo. 

Questa Extended Data Fig. 9 è molto utile perché “spacchetta” visivamente un punto che nelle altre figure viene dimostrato con correlazioni e autovalori: la fisica dei pattern (PDO e NPGO/Victoria Mode) resta riconoscibile, ma cambia l’ordine con cui emerge nelle EOF a seconda che tu lavori sulle SST totali oppure su SST con la GMSST rimossa. Nel pannello superiore (“Total SSTa”), cioè anomalie di SST che includono anche lo spostamento della baseline globale, l’EOF1 appare come un segnale ampio e coerente su gran parte del Nord Pacifico, con loadings prevalentemente dello stesso segno: questa è la firma di un contributo “a grande scala” che, in un clima in riscaldamento, tende a dominare la varianza e quindi ad occupare il primo autovettore. È esattamente il messaggio dello studio di Cluett et al.: il riscaldamento pan-bacino diventa statisticamente dominante nelle SST totali recenti pur senza che la PDO perda robustezza intrinseca. 

Sempre nella riga “Total SSTa”, l’EOF2 recupera invece una geometria tipicamente PDO-like, con opposizione di segno tra il settore nord-orientale (Golfo dell’Alaska / costa nordamericana) e porzioni più centrali–occidentali del bacino: è il tipo di impronta che storicamente ha reso la PDO un indice chiave per collegare variabilità del Pacifico e impatti su ecosistemi e risorse (per esempio salmonidi), a partire dai lavori classici di fine anni ’90. Il punto importante è che qui la PDO non “sparisce”: semplicemente non è più il primo modo quando nel campo SST lasci entrare il segnale globale, perché l’EOF è una decomposizione che ordina i modi per varianza spiegata; se aggiungi una componente spazialmente coerente e persistente (riscaldamento di fondo), quella componente tende a salire in cima alla gerarchia. Questo è anche il motivo per cui lo stesso studio insiste su controlli di robustezza (finestre mobili, correlazioni fra pattern/PC, autovalori): quando i modi sono vicini in varianza possono “mescolarsi” o scambiarsi ordine, un tema classico nella teoria degli errori di campionamento delle EOF. 

L’EOF3 della riga “Total SSTa” mostra un dipolo inclinato NE–SW nel Pacifico nord-orientale/subtropicale che gli autori associano a NPGO/Victoria Mode. Qui vale la pena essere precisi: in letteratura, la Victoria Mode (VM) è spesso definita come il secondo EOF delle SST extratropicali nordpacifiche, mentre la NPGO nasce come secondo EOF delle anomalie di sea surface height nel Pacifico centro-orientale; i due indici risultano però strettamente collegati e, in pratica, la NPGO “traccia” molto bene la seconda modalità SST (VM), con implicazioni note per circolazione di gyre, salinità, nutrienti e variabili ecosistemiche. Il fatto che, nelle SST totali, questa famiglia di variabilità appaia come EOF3 (e non EOF2) è coerente con l’idea che l’EOF1 venga “occupata” dal segnale pan-bacino e che l’EOF2 venga “presa” dalla PDO, spingendo la VM/NPGO un gradino più in basso.

La riga inferiore (“GMSST removed”) è la controprova concettuale: sottraendo la temperatura media globale del mare, l’EOF1 torna ad assomigliare chiaramente alla PDO classica, mentre la NPGO/Victoria Mode emerge come EOF2(come dichiara anche la didascalia). Questa configurazione è perfettamente coerente con la definizione operativa adottata da NOAA/PSL per la PDO, che rimuove esplicitamente la media globale proprio per isolare la variabilità interna del Nord Pacifico da qualunque segnale di “global warming” presente nei dati. In altre parole, la figura non sta solo dicendo “i pattern ci sono”: sta mostrando che la loro identità è stabile, ma la loro posizione nella classifica dei modi dipende dal fatto che tu stia analizzando un sistema “detrendizzato” (variabilità interna) oppure un sistema non stazionario in cui lo spostamento della baseline climatica è parte integrante del segnale osservato.