Guido Vettoretti1*, Ruei- Jia Hu2, Ingo Bethke3, Kirstin Krüger4,5, Michael Sigl4,6, Stephen Outten7
,
Jaimei Lin1, Roman Nuterman1, Anders Svensson1, Peter Ditlevsen1, Markus Jochum1
Il lavoro di Vettoretti e colleghi, pubblicato su Science Advances nel febbraio 2026, affronta una questione di grande rilievo per la paleoclimatologia dinamica: se le grandi eruzioni vulcaniche possano non soltanto perturbare il clima su scala annuale o decadale, ma anche contribuire all’innesco di transizioni climatiche persistenti su scala millenaria. L’elemento innovativo dello studio consiste nell’integrare un forcing vulcanico realistico, ricostruito e vincolato dai segnali di solfato conservati nelle carote di ghiaccio, in un ampio ensemble di simulazioni accoppiate atmosfera-oceano riferite al clima glaciale. In questo quadro gli autori mostrano che eruzioni equatoriali eccezionalmente intense possono alterare profondamente la AMOC, cioè la circolazione meridionale atlantica di ribaltamento, attraverso una catena di processi che coinvolge la circolazione atmosferica, la circolazione oceanica e soprattutto i flussi di galleggiamento aria-mare, con la possibilità di spingere il sistema climatico da uno stato relativamente mite a uno freddo, o viceversa, per tempi dell’ordine dei millenni. Il punto cruciale è che il vulcanismo non viene interpretato come causa unica della transizione, ma come perturbazione esterna capace di “attivare” una risposta non lineare di un sistema già vicino a una soglia critica.
Questo risultato si inserisce in una cornice teorica consolidata, secondo cui molte delle oscillazioni abrupthe dell’ultima glaciazione, in particolare gli eventi di Dansgaard-Oeschger, riflettono passaggi rapidi tra stati multipli della AMOC. Lo stesso articolo ricorda che durante il MIS 3 tali eventi produssero in Groenlandia aumenti termici dell’ordine di 10–15 °C nel giro di pochi decenni, ed è proprio questa rapidità a suggerire un comportamento di tipo “tipping” del sistema oceanico-atmosferico nordatlantico. Le ricostruzioni paleoceanografiche hanno progressivamente rafforzato questa interpretazione: la review di Lynch-Stieglitz evidenzia che il segnale globale dei cambiamenti abrupti dell’ultimo ciclo glaciale è coerente con forti variazioni della AMOC, mentre Böhm et al. mostrano che una modalità profonda e vigorosa della circolazione atlantica ha dominato gran parte dell’ultimo ciclo glaciale, con indebolimenti severi soprattutto durante gli Heinrich stadials e in corrispondenza di condizioni particolarmente favorevoli alla stratificazione e alla perturbazione della convezione nordatlantica. Ne emerge dunque un quadro in cui la circolazione overturning glaciale non è un sistema statico, ma una struttura metastabile, capace di permanere a lungo in uno stato robusto e tuttavia vulnerabile a forzanti esterne quando si avvicina a soglie dinamiche critiche.
L’aspetto forse più interessante del lavoro di Vettoretti et al. è che fornisce un meccanismo fisico plausibile per collegare la forzante vulcanica a tali transizioni. Lo studio indica che l’eruzione modifica dapprima i campi termici e dinamici dell’atmosfera e dell’oceano, quindi altera densità superficiale, salinità, copertura di ghiaccio marino e convezione profonda nelle aree chiave del Nord Atlantico. In altre parole, il segnale radiativo vulcanico viene amplificato e riorganizzato dal sistema oceano-ghiaccio marino, fino a produrre una risposta della AMOC molto più duratura del forcing iniziale. Questa interpretazione è coerente con il lavoro di Lohmann e Svensson, che usando archivi bipolari da carote di ghiaccio hanno mostrato che le grandi eruzioni avvengono con frequenza significativamente superiore al caso poco prima dell’inizio di vari eventi di Dansgaard-Oeschger; su 20 warmings abrupti tra 12 e 60 ka, 5 si verificano entro 20 anni e 7 entro 50 anni da una grande eruzione bipolare. A loro volta, studi più recenti sul carattere “state-dependent” della risposta vulcanica nel Nord Atlantico indicano che in clima glaciale il raffreddamento successivo a cluster eruttivi può risultare più intenso proprio perché dipende dalla stratificazione preesistente dell’oceano superiore e quindi dallo stato di fondo del sistema. In termini dinamici, ciò significa che la stessa eruzione non produce necessariamente lo stesso esito: la risposta dipende in modo decisivo da quanto il sistema sia già vicino a un cambio di regime.
Il contributo dello studio, quindi, non sta soltanto nell’attribuire un ruolo al vulcanismo, ma nel ridefinire il rapporto tra forcing esterno e variabilità interna durante le grandi oscillazioni glaciali. L’eruzione agisce come impulso relativamente breve, ma il suo effetto climatico di lunga durata emerge soltanto perché intercetta una AMOC dotata di multistabilità e di forte sensibilità ai flussi di galleggiamento e alla convezione profonda. Gli autori sottolineano inoltre che la variabilità naturale non forzata può, in alcune circostanze, esercitare un effetto stabilizzante nelle decadi successive all’eruzione, rallentando o impedendo il superamento della soglia critica. Questo punto è molto importante perché sposta l’attenzione da una visione deterministica a una probabilistica: il sistema glaciale non reagisce in modo lineare, ma in funzione della sua distanza dal tipping point. Anche per questo il lavoro dialoga direttamente con il dibattito contemporaneo sulla stabilità della AMOC, oggi considerata uno dei principali tipping elements del sistema climatico. Studi recenti basati su indicatori osservativi e teorici hanno rilanciato l’attenzione sul rischio di indebolimento critico della circolazione atlantica, mostrando come la nozione di soglia, isteresi e perdita di resilienza non appartenga solo al passato glaciale, ma costituisca un tema centrale anche per il clima attuale e futuro. In questa prospettiva, il lavoro di Vettoretti et al. è particolarmente prezioso perché suggerisce che gli archivi paleoclimatici non descrivono soltanto eventi remoti, ma rappresentano veri e propri laboratori naturali per comprendere come un sistema climatico vicino all’instabilità possa reagire a una perturbazione rapida e intensa.
INTRODUZIONE
Durante l’ultimo periodo glaciale, compreso approssimativamente tra 115.000 e 11.700 anni fa, il sistema climatico terrestre fu caratterizzato da una marcata instabilità millenaria, espressa dall’alternanza tra condizioni stadiali fredde e interstadiali relativamente più miti. La manifestazione più spettacolare di questa variabilità è rappresentata dagli eventi di Dansgaard–Oeschger, particolarmente frequenti e intensi durante il Marine Isotope Stage 3, tra circa 59 e 27 ka, quando le carote di ghiaccio groenlandesi documentano aumenti termici dell’ordine di 10–15 °C sviluppatisi nell’arco di pochi decenni. La letteratura più recente interpreta questi eventi non come semplici oscillazioni locali dell’Atlantico settentrionale, ma come transizioni rapide di un sistema climatico accoppiato, con propagazione di segnali atmosferici e idrologici su scala quasi globale; una sintesi moderna di questa prospettiva è offerta dal quadro multiscala discusso da Riechers et al., mentre la sincronizzazione di molte risposte climatiche in Asia monsonica, Sud America ed Europa-Mediterraneo entro scale decadali rispetto ai warming groenlandesi è stata evidenziata dalla compilazione speleotemica di Corrick et al.
All’interno di questo contesto, la AMOC è considerata il cardine dinamico principale della variabilità abrupta glaciale. Le review paleoceanografiche mostrano infatti che il paradigma dominante collega i passaggi rapidi tra stadiali e interstadiali a riorganizzazioni della circolazione overturning atlantica, cioè del trasporto verso nord di calore e salinità e della formazione di acque profonde nel Nord Atlantico. La cornice teorica che sostiene questa interpretazione deriva sia dall’evidenza paleoceanografica sia dai modelli di multistabilità: il sistema oceanico nordatlantico può esistere in stati relativamente robusti ma distinti, e il passaggio da uno stato all’altro può avvenire rapidamente quando il sistema viene spinto oltre una soglia critica. In questa prospettiva, gli eventi D-O costituiscono un esempio paleoclimatico di dinamica da tipping point, con l’oceano atlantico settentrionale che alterna fasi di convezione profonda relativamente vigorosa e fasi di circolazione più debole, fortemente influenzate da ghiaccio marino, stratificazione e forzanti di galleggiamento. La stessa distinzione tra condizioni interglaciali, pienamente glaciali e glaciali intermedie è rilevante, perché studi di stabilità indicano che la AMOC del MIS 3 era probabilmente più vulnerabile alle transizioni rispetto a quella del Last Glacial Maximum, offrendo una spiegazione fisicamente plausibile del fatto che la variabilità D-O sia stata più pronunciata in condizioni glaciali intermedie che in pieno glaciale.
Dal punto di vista dei meccanismi, durante gli stadiali freddi il Nord Atlantico ad alte latitudini era coperto da un’estesa banchisa invernale, in alcuni casi spinta molto a sud, e la presenza di una forte stratificazione alina nei mari del Labrador, di Groenlandia, di Irminger e di Norvegia favoriva il confinamento del calore atlantico negli strati di sottosuperficie. Questo accumulo graduale di calore nei primi centinaia di metri della colonna d’acqua costituiva una riserva energetica potenzialmente destabilizzante: una volta superata una soglia critica, il sistema poteva liberare rapidamente parte di quel contenuto termico verso l’atmosfera, innescando una transizione interstadiale con riattivazione della convezione profonda e rafforzamento della AMOC. A quel punto entravano in gioco feedback amplificanti, primo fra tutti il feedback ghiaccio marino–albedo, insieme all’intensificazione del trasporto meridiano di calore. Successivamente, l’aumento del deflusso di acqua dolce verso l’Artico e l’export di ghiaccio marino verso i mari GIN tendevano a ripristinare una forte stratificazione, sopprimendo la formazione di acque profonde e favorendo il ritorno a uno stato stadiale. Questa interpretazione ciclica, già presente nella letteratura precedente e ripresa nello studio di Vettoretti et al., è coerente con una visione non lineare del clima glaciale, in cui accumulo di calore oceanico, copertura di ghiaccio marino e idrologia superficiale modulano l’accesso del sistema a stati alternativi della circolazione nordatlantica.
In tale quadro, il vulcanismo emerge come possibile fattore di innesco, più che come causa unica, di alcune transizioni abrupthe. Il recente studio di Vettoretti et al. mostra, mediante simulazioni accoppiate atmosfera-oceano in condizioni glaciali forzate con una ricostruzione realistica delle grandi eruzioni basata sui solfati delle carote di ghiaccio, che eruzioni equatoriali molto grandi possono perturbare in modo sostanziale la AMOC attraverso modifiche della circolazione atmosferica e oceanica e dei flussi di galleggiamento aria-mare, fino a spingere il sistema tra stati persistenti freddi e caldi della durata di millenni. Questo risultato si inserisce bene in un filone già sviluppato da Lohmann e Svensson, i quali hanno mostrato che le grandi eruzioni bipolari compaiono con frequenza significativamente superiore al caso immediatamente prima dell’avvio di vari eventi D-O: su 20 warming abrupti tra 12 e 60 ka, 5 si collocano entro 20 anni e 7 entro 50 anni da un’eruzione bipolare. Lo stesso quadro è rafforzato dal catalogo glaciologico di Lin et al., che documenta 85 grandi eruzioni bipolari tra 60 e 9 ka, di cui 25 più grandi di qualsiasi eruzione degli ultimi 2500 anni e 69 con emissioni di zolfo superiori a quelle del Tambora del 1815. Ne deriva l’idea che il clima glaciale fosse soggetto non solo a oscillazioni interne, ma anche a impulsi esterni capaci di attivare risposte sproporzionate quando lo stato di fondo era già vicino a una soglia dinamica.
Questa lettura ha anche un valore concettuale per il clima contemporaneo. Il nesso tra eventi D-O, soglie della AMOC e risposta del sistema a perturbazioni relativamente brevi rende infatti il record glaciale un laboratorio naturale per comprendere la sensibilità dell’overturning atlantico. Gli studi recenti sulla AMOC moderna, inclusi quelli che propongono indicatori fisici precoci di perdita di stabilità, ribadiscono che la circolazione atlantica è oggi considerata uno dei principali tipping elements del sistema climatico. Per questo motivo l’analisi del ruolo del vulcanismo nelle dinamiche glaciali non riguarda soltanto il passato remoto: aiuta anche a chiarire in che modo un sistema climatico già predisposto all’instabilità possa reagire a forzanti esterne o a fluttuazioni interne, e perché la distanza dalla soglia critica conti tanto quanto l’intensità del disturbo stesso. Il possibile ruolo delle grandi eruzioni vulcaniche nell’avvio o nella modulazione della variabilità climatica su scala millenaria è discusso nella letteratura paleoclimatica da diversi decenni, ma solo recentemente il tema ha iniziato a essere affrontato con strumenti cronologici e modellistici sufficientemente robusti da permettere valutazioni meno speculative. Per lungo tempo, infatti, molte ipotesi di collegamento tra vulcanismo e transizioni abrupthe del clima glaciale si sono basate soprattutto su coincidenze temporali tra proxy climatici e segnali di attività eruttiva, senza che fosse disponibile una solida ricostruzione dei tempi, della magnitudo delle emissioni solforate e, soprattutto, dei passaggi fisici capaci di tradurre un forcing radiativo di durata relativamente breve in una risposta climatica persistente. Il nuovo studio di Vettoretti et al. si inserisce proprio in questo vuoto interpretativo, mostrando che il problema non è soltanto stabilire se esista una correlazione tra eruzioni e variabilità tipo Dansgaard–Oeschger, ma comprendere in quali condizioni dinamiche un’eruzione possa agire da innesco o da amplificatore di un sistema climatico glaciale già prossimo a una soglia critica della AMOC. In questo senso, il vulcanismo non viene trattato come causa esclusiva delle transizioni millenarie, bensì come perturbazione esterna capace di spingere un sistema multistabile verso un cambio di regime quando la circolazione oceanica, la copertura di ghiaccio marino e i flussi di galleggiamento aria-mare si trovano in una configurazione particolarmente sensibile.
Uno dei motivi per cui questa ipotesi è stata a lungo difficile da valutare riguarda la qualità e la completezza degli archivi vulcanici del passato. Le stime petrologiche del rilascio di zolfo stratosferico per le grandi eruzioni quaternarie possono divergere enormemente, e l’attribuzione delle anomalie di solfato osservate nelle carote di ghiaccio a specifici vulcani sorgente resta in molti casi incerta. A ciò si aggiunge il fatto che il trasporto atmosferico degli aerosol verso Groenlandia e Antartide, nonché la loro efficienza di scavenging lungo il tragitto, dipendono dallo stato del clima di fondo: durante gli stadiali, più secchi, la ridotta precipitazione tende a limitare la rimozione degli aerosol, mentre durante gli interstadiali più umidi la maggiore disponibilità di vapore acqueo può favorire una scavenging più efficiente. Proprio per questo le carote di ghiaccio ad alta risoluzione costituiscono oggi una risorsa eccezionale: nello stesso archivio preservano infatti sia il segnale vulcanico, tramite solfati e cenere, sia variabili climatiche cruciali come temperatura, polvere e metano, consentendo confronti diretti tra forcing e risposta climatica con una precisione temporale molto superiore a quella di gran parte degli altri archivi paleoclimatici. I lavori di Lohmann e Svensson e, più recentemente, di Lohmann et al. hanno mostrato che i record bipolari sincronizzati permettono non solo di individuare grandi eruzioni durante gli ultimi 60 mila anni, ma anche di testare la dipendenza dallo stato climatico della risposta registrata nelle serie isotopiche polari, un passaggio essenziale per separare semplici coincidenze da legami dinamicamente plausibili.
Dal punto di vista fisico, il forcing vulcanico agisce innanzitutto attraverso l’iniezione nella stratosfera di anidride solforosa, cenere e particolato, con successiva formazione di aerosol solfatici capaci di aumentare la riflessione della radiazione solare incidente e quindi di raffreddare la superficie terrestre per alcuni anni. Tuttavia, il punto decisivo emerso dalla letteratura è che gli effetti del vulcanismo non si esauriscono nel raffreddamento radiativo immediato. Studi modellistici e osservativi indicano che grandi eruzioni possono influenzare il contenuto di calore oceanico, l’estensione del ghiaccio marino, i gradienti termici tropicali e subtropicali e, tramite questi, l’intera dinamica accoppiata oceano-atmosfera. Nel Nord Atlantico, per esempio, il raffreddamento superficiale, la riduzione delle precipitazioni e del runoff e l’aumento della salinità superficiale possono favorire una maggiore densità delle acque, destabilizzare la colonna d’acqua e intensificare la convezione profonda, con un rafforzamento transitorio della AMOC in alcuni contesti modellistici. Van Dijk et al. mostrano che sequenze eruttive ravvicinate possono sostenere anomalie di raffreddamento attraverso feedback oceano–ghiaccio marino su tempi molto più lunghi del forcing iniziale, mentre studi sul clima storico e olocenico hanno evidenziato come cluster di eruzioni possano modulare fasi fredde multidecadali o multicentenarie nell’Atlantico settentrionale. In questo quadro, l’impatto di una singola eruzione o di una sequenza eruttiva dipende in modo cruciale dalla struttura preesistente del sistema, in particolare dalla stratificazione del Nord Atlantico e dall’estensione del ghiaccio marino, due elementi centrali anche nella dinamica degli eventi D-O.
Un altro aspetto di rilievo riguarda il ruolo delle teleconnessioni atmosferiche e oceaniche. La risposta della AMOC al forcing vulcanico non sembra dipendere soltanto dai processi locali di convezione e densificazione, ma anche da aggiustamenti su larga scala che coinvolgono la North Atlantic Oscillation, il Pacifico tropicale, i gradienti di temperatura superficiale nel tropical Atlantic e la posizione della ITCZ. Nello stesso studio di Vettoretti et al. si sottolinea che grandi eruzioni equatoriali possono modificare la circolazione atmosferica e oceanica e i flussi di galleggiamento aria-mare in modo tale da spingere il sistema climatico glaciale tra stati persistenti caldi e freddi. Questo punto è particolarmente importante perché suggerisce che l’effetto di un’eruzione non debba essere pensato come un semplice raffreddamento globale uniforme, ma come una perturbazione in grado di riorganizzare i trasporti energetici tra tropici e alte latitudini, anche attraverso meccanismi riconducibili alla compensazione di Bjerknes. In un sistema già caratterizzato da multistabilità, come quello dell’Atlantico glaciale, una simile riorganizzazione può produrre conseguenze assai più durevoli della vita media degli aerosol stratosferici. Non sorprende quindi che gli autori insistano sulla necessità di integrare i dati delle carote di ghiaccio con simulazioni accoppiate ad alta complessità: i soli archivi proxy, pur fondamentali, non possono descrivere pienamente la propagazione spaziale del segnale né chiarire in modo univoco perché alcune eruzioni restino climaticamente “inermi” mentre altre riescano a innescare o amplificare una transizione abrupta.
Alla luce di questi risultati, la continuazione dello studio si colloca in una linea di ricerca molto promettente, nella quale il record vulcanico delle carote di ghiaccio non è più soltanto un catalogo di eventi atmosferici, ma diventa una chiave per analizzare la resilienza del sistema climatico glaciale. Il contributo più innovativo non sta soltanto nell’ipotizzare che il vulcanismo possa precedere alcuni warming Dansgaard–Oeschger, cosa già suggerita dal lavoro statistico di Lohmann e Svensson, ma nel tentativo di spiegare come la risposta dipenda dal background climatico e dalla vicinanza del sistema a una soglia della AMOC. In questo senso, il segmento successivo del paper si inserisce coerentemente in una prospettiva process-based: chiarire se le eruzioni abbiano agito da semplice forcing radiativo transitorio, da fattore di amplificazione di oscillazioni interne già in atto, oppure da vero e proprio trigger capace di convertire una perturbazione breve in una riorganizzazione millenaria della circolazione atlantica.
RISULTATI
Nel segmento dei risultati, Vettoretti e colleghi costruiscono anzitutto un ponte quantitativo tra archivio paleovulcanico e sperimentazione modellistica, mostrando che i record di solfato delle carote di ghiaccio per Olocene e MIS 3 possono essere trattati con strumenti propri della teoria dei valori estremi. L’uso di una distribuzione generalized extreme value (GEV) è metodologicamente coerente, perché questo tipo di funzione è ampiamente impiegato per descrivere la coda degli eventi rari ma ad alto impatto e per stimarne i livelli di ritorno; nel caso specifico, gli autori trovano che la frequenza delle grandi e grandissime eruzioni si adatta bene a una GEV, pur con una sottostima nella coda più estrema, cioè per gli intervalli di ricorrenza più lunghi. Il risultato più rilevante è che un evento equivalente a 230 Tg di aerosol vulcanico, scelto come forcing di riferimento nelle simulazioni, presenta un tempo di ritorno dell’ordine di circa 1500 anni nell’Olocene e superiore a 3000 anni nel MIS 3, mentre se si considera l’equivalente deposizione di circa 217 Tg di SO₄ osservabile nelle carote di ghiaccio la ricorrenza effettiva si avvicina a ~1000 anni. Questo scarto tra Olocene e MIS 3 non va però letto in modo ingenuamente climatico, perché riflette anche differenze metodologiche nelle ricostruzioni: nel record glaciale il maggiore assottigliamento degli strati annuali con la profondità peggiora il rapporto segnale-rumore e rende plausibile che le stime di ricorrenza per il MIS 3 siano conservative. In questo senso, la figura 1 non si limita a descrivere quante eruzioni estreme si osservano nei proxy, ma definisce il quadro probabilistico entro cui valutare quanto realistico sia un forcing “Samalas-like” applicato a simulazioni glaciali di lunga durata.
Su questa base statistica, il lavoro passa alla dinamica climatica vera e propria usando il CCSM4 in configurazione a bassa risoluzione, scelta necessaria per integrare simulazioni accoppiate atmosfera-oceano su tempi plurimillenari in condizioni glaciali. Un elemento importante è che il massimo raffreddamento globale indotto da una grande eruzione equatoriale di tipo Samalas resta sostanzialmente simile indipendentemente dallo stato di fondo, moderno o glaciale, stadiale o interstadiale: l’anomalia media globale raggiunge circa −2,5 °C, con una profondità ottica media globale degli aerosol prossima a 1,75 a 550 nm pochi mesi dopo l’evento. Questo implica che la differenza tra i diversi esperimenti non nasce tanto dall’intensità radiativa iniziale del forcing, quanto dalla maniera in cui quello stesso impulso viene filtrato dal background climatico, dalla stratificazione oceanica e dalla sensibilità della AMOC. È un passaggio cruciale, perché suggerisce che il problema scientifico non sia soltanto “quanto raffredda” una super-eruzione, ma “quanto efficacemente” quel raffreddamento viene convertito in una perturbazione dell’oceano nordatlantico e del ghiaccio marino. La letteratura recente sulle ricostruzioni del forcing vulcanico stratosferico conferma del resto che il segnale dei grandi eventi è ben registrato dagli archivi glaciali e olocenici, mentre i lavori dedicati al Samalas del 1257 mostrano perché tale eruzione costituisca un riferimento plausibile per un esperimento numerico estremo ma fisicamente credibile.
Un ulteriore aspetto di notevole interesse riguarda il modo in cui gli autori contestualizzano cronologicamente il MIS 3 attraverso il matching bipolare delle deposizioni di solfato tra Groenlandia e Antartide. Questo approccio, ormai centrale nella paleoclimatologia del Quaternario recente, permette di vincolare con precisione decadale la sequenza degli eventi abrupti e di collocare le anomalie vulcaniche all’interno del ritmo dei Dansgaard–Oeschger. La robustezza di tale strategia è stata rafforzata dai lavori di Svensson e collaboratori, che hanno identificato numerose eruzioni bipolari nel periodo glaciale, e da ricostruzioni più recenti che hanno ampliato ulteriormente la sincronizzazione interemisferica. In questo quadro, il richiamo a eruzioni come Mazama/Crater Lake o alla super-eruzione Oruanui del Taupō non è meramente illustrativo: serve a mostrare che anche eventi extratropicali di grandissima magnitudo possono lasciare firme solfatiche e, in alcuni casi, tefritiche su scala emisferica o interemisferica. Ciò rafforza l’idea che il record delle carote di ghiaccio non sia solo un catalogo di aerosol atmosferici, ma un archivio fisicamente integrato in cui il forcing vulcanico può essere confrontato direttamente con temperatura, polveri e altri indicatori climatici. In termini interpretativi, questa parte dei risultati prepara quindi il terreno a una domanda più profonda, che il segmento successivo affronta in chiave dinamica: non soltanto se le grandi eruzioni fossero abbastanza intense da lasciare una traccia nei ghiacci, ma se fossero abbastanza rare, abbastanza grandi e collocate in una fase abbastanza sensibile del sistema glaciale da alterare il ciclo stesso della AMOC. In questa parte dei risultati emerge con particolare chiarezza che il ruolo del forcing vulcanico nel modulare la variabilità di tipo Dansgaard–Oeschger non può essere interpretato in modo lineare o universalmente deterministico. Gli autori sottolineano anzitutto che, pur avendo adottato una ricostruzione della sorgente vulcanica coerente con i segnali registrati nelle carote di ghiaccio, una trattazione più completa delle risposte climatiche glaciali dovrebbe includere in futuro la stagionalità dell’iniezione aerosolica, la diversa distribuzione emisferica del forcing e il trasporto dinamico degli aerosol nella stratosfera. Questo punto è cruciale perché, nel sistema glaciale, gli effetti di un’eruzione non dipendono solo dall’intensità radiativa globale, ma anche da come il segnale viene ridistribuito tra emisfero nord ed emisfero sud. In tale quadro si inserisce il cosiddetto bipolar seesaw, cioè il ben noto accoppiamento interemisferico per cui un raffreddamento nel Nord Atlantico associato a un indebolimento della AMOC tende a corrispondere, con proprie scale temporali e modalità fisiche, a un riscaldamento relativo dell’emisfero australe. Le review più autorevoli su questo tema mostrano infatti che il seesaw termico non è un semplice schema descrittivo, ma il risultato di una riorganizzazione dei trasporti di calore oceanici e atmosferici, con compensazioni energetiche che coinvolgono anche l’oceano globale e la criosfera australe. Per questo motivo, simulare eruzioni con forcing confinato emisfericamente, o valutare eventi extratropicali di grande entità, richiede modelli in cui il trasporto degli aerosol e la risposta dinamica dell’atmosfera siano trattati in modo esplicito, aspetto che nel lavoro di Vettoretti et al. viene correttamente riconosciuto come limite metodologico e come naturale sviluppo successivo dell’analisi.
Il cuore del segmento, tuttavia, è lo studio di sensibilità con forcing vulcanico equatoriale periodico, pensato per verificare se e in quali condizioni una grande eruzione possa perturbare il ciclo D-O. Gli autori impongono un’eruzione da 230 Tg di aerosol vulcanico ogni 1000 anni e confrontano poi questa configurazione con scenari ridotti a 115 Tg e 58 Tg. Il risultato è molto istruttivo: la risposta del sistema climatico glaciale si mostra nettamente più evidente per il forcing estremo, mentre al di sotto di circa 115 Tg, quindi per intensità paragonabili a un evento di ordine Tambora, l’effetto sul ciclo D-O diventa fortemente attenuato; con 58 Tg, cioè per una magnitudo assimilabile a Pinatubo, la risposta non è più distinguibile dalle oscillazioni spontanee della variabilità interna. Questo implica che, nel modello, non tutte le grandi eruzioni possiedono la capacità di alterare sensibilmente la traiettoria della AMOC o di innescare un ritorno rapido verso uno stato stadiale. Il dato importante non è soltanto la presenza di una soglia di efficacia del forcing, ma il fatto che tale soglia emerga in un sistema già caratterizzato da multistabilità, isteresi e forte dipendenza dallo stato di fondo. In letteratura, gli eventi D-O sono sempre più spesso interpretati proprio come manifestazioni altamente non periodiche di un sistema dinamico non lineare, in cui durate molto diverse degli interstadiali e la presenza di eventi precursori brevi possono verificarsi anche sotto condizioni esterne simili. In questo senso, il risultato numerico di Vettoretti et al. si accorda molto bene con i più recenti modelli concettuali e con gli studi sul noise-induced tipping della AMOC, secondo i quali la vicinanza a una soglia critica rende il sistema estremamente sensibile non solo al forcing esterno, ma anche alla struttura temporale della variabilità interna.
Da qui deriva forse il messaggio scientificamente più interessante di questa continuazione: la risposta alle eruzioni appare probabilistica, non deterministica. Durante la lenta fase di decadimento di un interstadiale D-O, la variabilità interna può produrre finestre decadali in cui la AMOC è insolitamente robusta; se un’eruzione interviene in una di queste finestre, il sistema può mostrare una risposta minima o addirittura nulla. In altri momenti dello stesso ciclo, invece, una perturbazione analoga può favorire un comportamento di tipping. Gli autori mostrano infatti che due eruzioni identiche inserite in fasi diverse dell’interstadiale non generano necessariamente lo stesso esito, a riprova del fatto che il sistema non reagisce meccanicamente all’ampiezza del forcing, ma alla combinazione tra forcing, stato medio e rumore climatico interno. Ancora più significativa è l’osservazione che la fase iniziale dell’interstadiale, caratterizzata da un overshoot della AMOC durante la transizione rapida dal regime stadiale freddo a quello interstadiale caldo, rende il sistema più resiliente verso perturbazioni vulcaniche che in altre fasi potrebbero favorire il ritorno al freddo. Questo comportamento è coerente con la moderna teoria dei tipping elements: la resilienza di un sistema vicino a una biforcazione non è costante nel tempo, ma oscilla con la traiettoria interna del sistema stesso. Per questo il paper insiste sulla necessità di usare ensemble molto ampi: solo un grande numero di simulazioni permette infatti di stimare in modo statisticamente robusto la probabilità che una stessa classe di eruzione produca o meno una transizione della AMOC in un clima glaciale dominato da forte variabilità interna e da stati multipli della circolazione overturning.
Figura 1. Record del vulcanismo esplosivo e quadro probabilistico degli eventi estremi
La figura 1 svolge una funzione metodologica centrale nello studio, perché non si limita a mostrare una successione di eruzioni, ma costruisce il legame tra archivio paleovulcanico, variabilità climatica glaciale e scelta del forcing applicato nelle simulazioni del CCSM4. Nei pannelli A e B sono riportate le serie temporali del carico stratosferico di solfati ricostruito dalle carote di ghiaccio rispettivamente per l’Olocene e per il MIS 3, mentre nel pannello B compare anche la temperatura NGRIP della Groenlandia in verde, utile a collocare i picchi eruttivi rispetto alla successione di stadiali e interstadiali del ciclo Dansgaard–Oeschger. Il messaggio principale è che il record glaciale e quello olocenico contengono una popolazione di eventi rari ma molto intensi, con contributi sia equatoriali sia extratropicali, e che proprio le carote di ghiaccio restano oggi uno degli archivi più robusti per ricostruire il forcing aerosolico pre-strumentale, perché consentono di confrontare nello stesso supporto naturale deposizione vulcanica e risposta climatica. Questo approccio si innesta direttamente sui progressi ottenuti dalle ricostruzioni oloceniche ad alta risoluzione e dalla sincronizzazione bipolare Groenlandia–Antartide, che hanno migliorato notevolmente la cronologia degli eventi eruttivi e delle oscillazioni abrupthe del clima glaciale.
Dal punto di vista fisico, il pannello B è probabilmente il più importante, perché mette in relazione il catalogo del vulcanismo del MIS 3 con una delle più note espressioni della variabilità millenaria del Quaternario recente, cioè i bruschi riscaldamenti groenlandesi. La linea verde mostra infatti la forte oscillazione della temperatura NGRIP associata agli eventi D-O, che nei record groenlandesi corrispondono a transizioni molto rapide verso condizioni più miti all’interno di un clima glaciale fortemente instabile; la sovrapposizione con i picchi di solfato non implica automaticamente causalità, ma serve a visualizzare l’ipotesi che alcune grandi eruzioni possano aver agito come perturbazioni capaci di modificare la probabilità di una transizione quando la AMOC era prossima a una soglia di stabilità. In questa prospettiva, la figura 1 non va letta come una prova diretta che ogni grande eruzione inneschi un evento D-O, bensì come il fondamento empirico di una domanda dinamica più sottile: se e quando un forcing vulcanico eccezionale possa interferire con un sistema già caratterizzato da multistabilità, feedback di ghiaccio marino e forte sensibilità della circolazione overturning nordatlantica. Questa interpretazione è coerente con il lavoro di Lohmann e Svensson, secondo cui le grandi eruzioni bipolari compaiono più frequentemente del caso immediatamente prima dell’inizio di diversi eventi D-O, pur senza esaurire da sole la spiegazione dell’intera variabilità glaciale.
I pannelli C e D traducono poi il record osservativo in termini statistici, mostrando che la distribuzione delle grandi eruzioni può essere descritta con una generalized extreme value distribution, cioè una distribuzione dei valori estremi adatta a rappresentare eventi rari ma ad alto impatto. Gli autori osservano che la frequenza degli eventi vulcanici si adatta bene a una GEV, sebbene con una certa sottostima nella parte più estrema della coda, e usano questo schema per derivare i tempi di ritorno. È qui che la figura acquista un valore decisivo per l’impostazione degli esperimenti: un’eruzione equivalente a 230 Tg di aerosol vulcanico, cioè il forcing più intenso usato nelle simulazioni, corrisponde a un tempo di ritorno dell’ordine di circa 1500 anni nell’Olocene e superiore a 3000 anni nel MIS 3; tuttavia, tenendo conto che il carico modellistico equivale a una deposizione nelle carote di ghiaccio di circa 217 Tg di SO₄, il periodo di ritorno osservazionale associato scende verso ~1000 anni. Gli stessi autori precisano però che le stime per il MIS 3 sono probabilmente conservative, perché l’assottigliamento degli strati annuali con la profondità degrada il rapporto segnale-rumore e rende più difficile ricostruire con continuità il forcing glaciale. Ne risulta che la figura 1 non serve soltanto a descrivere “quanto vulcanismo” ci sia stato, ma soprattutto a mostrare che il forcing scelto nel modello rappresenta un evento eccezionale ma fisicamente plausibile su scala millenaria, abbastanza raro da essere climatically special e abbastanza realistico da costituire un test credibile della sensibilità della AMOC glaciale.
Un ultimo elemento importante è dato dalle tre linee orizzontali tratteggiate del pannello B, che indicano i livelli di forcing usati nel CCSM4: 230, 115 e 58 Tg. Questi livelli corrispondono, nello schema sperimentale del paper, a un’eruzione molto grande di tipo Samalas, a un evento di ordine Tambora e a un forcing di ordine Pinatubo. La figura anticipa quindi già la logica del segmento successivo dei risultati: verificare se il sistema glaciale risponda solo agli eventi veramente estremi oppure anche a eruzioni più “ordinarie”. Il fatto che gli autori abbiano ancorato tali esperimenti a distribuzioni di ritorno ricavate dai record di ghiaccio rende la figura 1 non una semplice introduzione grafica, ma la base quantitativa dell’intero impianto numerico dello studio.
La figura 2 è una delle parti più istruttive del lavoro, perché mostra in modo diretto che il ciclo Dansgaard–Oeschger simulato dal modello non risponde al forcing vulcanico in modo lineare, ma secondo una dinamica di soglia fortemente dipendente dallo stato della AMOC. In tutti e tre i pannelli la curva blu rappresenta il controllo senza vulcanismo, mentre le curve colorate mostrano esperimenti in cui il modello viene perturbato ogni 1000 anni da un’eruzione equatoriale a partire dall’anno 2601. Il parametro seguito è l’AMOC max, espresso in Sverdrup, cioè l’intensità massima della circolazione di overturning atlantica. Nel pannello A, con forcing di 230 Tg di aerosol vulcanico, la traiettoria rossa si discosta più volte in modo marcato dal controllo, mostrando indebolimenti pronunciati della AMOC e successive fasi di recupero; è proprio questo comportamento che dà senso al titolo del lavoro, perché suggerisce che eruzioni equatoriali molto grandi possano occasionalmente spingere il sistema glaciale da uno stato interstadiale relativamente forte verso uno stato più debole e stadiale, prima di un successivo riassetto. Questa lettura è coerente con il quadro paleoclimatico generale secondo cui gli eventi D-O riflettono riorganizzazioni rapide della AMOC all’interno di un sistema nordatlantico multistabile, sensibile a ghiaccio marino, stratificazione e flussi di galleggiamento superficiali.
Il confronto con i pannelli B e C chiarisce però che non conta solo la presenza del forcing, ma soprattutto la sua magnitudo. Quando l’ampiezza dell’impulso viene ridotta a 115 Tg, cioè a un ordine di grandezza assimilabile a un’eruzione tipo Tambora, la risposta della curva verde diventa molto più attenuata: in alcuni intervalli il sistema devia dal controllo, ma in altri resta vicino alla sua traiettoria spontanea. Quando poi il forcing scende a 58 Tg, comparabile a un evento tipo Pinatubo, la curva viola è sostanzialmente indistinguibile dalla variabilità interna del controllo. Il messaggio fisico è molto importante: nel modello esiste una soglia oltre la quale il forcing vulcanico può modificare in modo apprezzabile la probabilità di una transizione del ciclo D-O, mentre al di sotto di quella soglia il segnale viene assorbito dal rumore climatico interno. Questo risultato si inserisce bene nella letteratura più recente, che considera gli eventi Dansgaard–Oeschger come manifestazioni di un sistema climatico non lineare, nel quale la risposta a una perturbazione esterna dipende dalla vicinanza a una soglia critica e non soltanto dall’ampiezza del forcing stesso.
La figura 2 è particolarmente convincente anche perché mostra che persino nel caso estremo da 230 Tg la risposta non è perfettamente ripetibile. Eruzioni identiche, inserite in momenti diversi del ciclo interstadiale, non producono necessariamente lo stesso esito sulla AMOC: talvolta il sistema cede verso uno stato più debole, talvolta oppone una notevole resistenza. Questo significa che la risposta è probabilistica e state-dependent, come sottolineano esplicitamente gli autori. In termini dinamici, durante la fase iniziale dell’interstadiale, quando la AMOC presenta un overshoot dopo la transizione rapida dal freddo al caldo, il sistema appare più resiliente; durante la fase di decadimento lento dell’interstadiale, invece, il background può diventare più favorevole a un tipping indotto da una perturbazione vulcanica. È esattamente questo il punto che collega la figura 2 ai lavori di Lohmann e Svensson, i quali avevano già mostrato che grandi eruzioni bipolari compaiono più spesso del caso poco prima dell’avvio di alcuni eventi D-O: il vulcanismo non genera automaticamente il cambiamento abrupto, ma può agire come trigger quando il sistema è già vicino a una soglia di instabilità.
In definitiva, la figura 2 non dimostra soltanto che le grandi eruzioni possono perturbare il clima glaciale, ma precisa qualcosa di più sottile e più utile dal punto di vista dinamico: l’efficacia del forcing vulcanico dipende dall’interazione tra magnitudo dell’eruzione, fase del ciclo D-O e variabilità interna della AMOC. È per questo che gli autori insistono sulla necessità di ensemble molto ampi: in un sistema governato da multistabilità, isteresi e rumore interno, una singola simulazione non basta a distinguere un effetto causale robusto da una semplice coincidenza di fase. La figura, quindi, rappresenta un passaggio chiave dell’intero articolo, perché trasforma l’idea generale di “vulcanismo come possibile innesco” in una dimostrazione modellistica di sensibilità di soglia, con implicazioni molto rilevanti per l’interpretazione degli eventi abrupti del Quaternario.
Impatti vulcanici sul clima interstadiale e stadiale
Nel quadro interpretativo proposto da Vettoretti e colleghi, il clima glaciale del MIS 3 non viene trattato come un sistema che risponde in modo lineare a una singola forzante esterna, ma come un sistema dinamico complesso dotato di soglie, isteresi e transizioni probabilistiche tra stati quasi stabili. La simulazione di riferimento, impostata con valori di CO₂ coerenti con il primo MIS 3, circa 230 ppm, produce una lunga fase interstadiale e una fase stadiale più breve, fornendo così il contesto da cui vengono fatti partire grandi ensemble perturbati. In questo impianto sperimentale, eruzioni equatoriali molto intense vengono introdotte in diversi momenti del ciclo Dansgaard–Oeschger per verificare se un forcing impulsivo di breve durata possa alterare la stabilità della circolazione meridionale atlantica e favorire un passaggio verso uno stato climatico persistente alternativo. Il risultato più importante è che l’eruzione non agisce come “causa automatica” del cambiamento abrupto: essa aumenta piuttosto la probabilità di transizione quando il sistema è già vicino a una soglia critica. Questa lettura è coerente con la visione moderna dell’AMOC come componente del sistema climatico capace di esibire forte non linearità, ampia sensibilità alle condizioni di fondo e possibili passaggi rapidi tra configurazioni di intensità molto diversa.
Dal punto di vista dinamico, il cuore della sottosezione è l’ipotesi che durante le fasi glaciali il sistema nord-atlantico possa trovarsi in un regime di bistabilità o quasi-bistabilità, nel quale una AMOC vigorosa e una più debole rappresentano due configurazioni climatiche relativamente persistenti. Studi precedenti avevano già indicato che le oscillazioni D-O sono strettamente legate a riorganizzazioni della circolazione atlantica e che la stabilità di tale circolazione cambia profondamente tra condizioni glaciali e interglaciali. In particolare, Sévellec e Fedorov mostrano che durante i glacials l’AMOC tende a diventare intrinsecamente più instabile, e che questa instabilità è fortemente controllata dalla posizione media del margine del ghiaccio marino, la quale modifica il rapporto tra sito di formazione delle acque profonde, bilancio di precipitazione-evaporazione e feedback oceanici. Anche il lavoro di Buizert e colleghi sottolinea che la stabilità dell’AMOC glaciale dipende da controlli remoti e dal background oceanico globale, compreso il ruolo dell’Oceano Meridionale, il che rafforza l’idea che il Nord Atlantico glaciale non fosse un sistema semplice, ma una struttura prossima a biforcazioni e suscettibile a forzanti relativamente brevi purché dinamicamente ben collocate nel tempo.
In questa prospettiva, il forcing vulcanico diventa scientificamente rilevante non tanto per la sua durata radiativa diretta, limitata a pochi anni, quanto per la sua capacità di innescare una cascata di feedback aria-mare-ghiaccio. Il raffreddamento indotto dagli aerosol stratosferici aumenta l’albedo planetaria, favorisce un’estensione del ghiaccio marino e modifica i flussi di galleggiamento superficiale e la convezione nel Nord Atlantico; da qui può originarsi una risposta dell’AMOC che inizialmente può anche rafforzarsi, per poi indebolirsi su scala decadale fino a spingere il sistema oltre una soglia dinamica. Questo elemento è particolarmente raffinato, perché mostra che l’effetto dell’eruzione non coincide banalmente con un semplice raffreddamento superficiale, ma con una riorganizzazione transitoria della circolazione accoppiata atmosfera-oceano. La plausibilità di un ruolo vulcanico come fattore di innesco riceve inoltre un supporto indipendente dai dati paleoclimatici: Lohmann e Svensson hanno mostrato che, tra 12 e 60 ka, grandi eruzioni bipolari compaiono più frequentemente del caso immediatamente prima dell’inizio di alcuni eventi di riscaldamento D-O, con 5 eventi su 20 entro 20 anni e 7 entro 50 anni da una grande eruzione. Nel complesso, quindi, la sottosezione suggerisce che i vulcani possano aver agito come perturbazioni brevi ma dinamicamente efficaci, capaci di anticipare una transizione già resa possibile dalla vulnerabilità interna del sistema glaciale.
Un altro contributo importante del brano è l’insistenza sulla natura probabilistica della transizione. Le simulazioni senza forcing vulcanico mostrano infatti che, quando il sistema si avvicina al tipping point, cresce già da solo il numero di transizioni indotte dalla variabilità interna; le simulazioni con forcing vulcanico aumentano ulteriormente questa probabilità perché la perturbazione è molto più intensa del rumore climatico di fondo. Questo quadro richiama la tradizione interpretativa della stochastic resonance, secondo cui il rumore interno può cooperare con forcing esterni deboli o intermittenti nel produrre salti tra stati climatici distinti, e si collega anche alla più ampia letteratura sui fenomeni abrupti del Pleistocene come espressione di dinamiche soglia, overshoot e rilassamento. Tuttavia, la letteratura più recente mostra che il dibattito è ancora aperto: accanto ai modelli bistabili esistono proposte alternative in cui i cicli D-O emergono da una dinamica monostabile ma eccitabile, nella quale grandi perturbazioni o rumore intermittente generano escursioni di stato simili a interstadiali senza richiedere necessariamente due attrattori stabili pienamente distinti. Questo non indebolisce il valore della sottosezione; al contrario, la rende ancora più interessante, perché il risultato di Vettoretti et al. può essere letto come evidenza che il forcing vulcanico interagisce con un sistema dotato di forte non linearità, indipendentemente dal fatto che tale non linearità venga formalizzata come bistabilità classica o come eccitabilità multiscala.
Nel suo significato più profondo, dunque, questa sottosezione suggerisce che gli eventi abrupti del tardo Pleistocene non vadano interpretati come fenomeni rigidamente periodici né come semplici risposte deterministicamente forzate dall’esterno, ma come il prodotto dell’interazione tra forzanti impulsive, rumore interno, feedback del ghiaccio marino, convezione oceanica e stabilità variabile dell’AMOC. Le grandi eruzioni vulcaniche, in questo schema, non costruiscono da sole il ciclo Dansgaard–Oeschger, ma possono modificarne il timing, amplificarne l’espressione e, in circostanze favorevoli, spingere il sistema oltre una soglia che conduce a un nuovo regime climatico persistente per secoli o millenni. È proprio questo passaggio, dalla semplice nozione di “raffreddamento da aerosol” alla concezione di “forzante impulsiva capace di deviare la traiettoria di un sistema vicino al tipping point”, a rendere il testo particolarmente rilevante anche per la climatologia contemporanea, perché fornisce un esempio paleoclimatico concreto di come perturbazioni relativamente brevi possano produrre conseguenze di lunga durata quando agiscono su un sistema climatico già predisposto all’instabilità
Figura 3 – Probabilità di transizione dell’AMOC durante fasi interstadiali e stadiali sotto forcing vulcanico
La figura 3 mostra con grande chiarezza che, nel modello CCSM4, l’effetto di una grande eruzione vulcanica equatoriale non è quello di imporre automaticamente una transizione climatica, ma di modificarne la probabilità in funzione dello stato dinamico in cui si trova il sistema. I pannelli A ed E rappresentano la simulazione di controllo: nella prima parte domina una fase interstadiale lunga, con AMOC forte, mentre successivamente compare una fase stadiale più breve, con AMOC debole. Su questa traiettoria di fondo gli autori costruiscono, in sei momenti diversi del ciclo D-O, due ensemble paralleli da 50 membri ciascuno: uno con eruzione molto intensa e uno senza eruzione, così da separare il contributo del forcing vulcanico da quello della sola variabilità interna. Gli istogrammi e le curve di densità a destra di ciascun pannello sintetizzano la distribuzione dell’intensità finale dell’AMOC 300 anni dopo la perturbazione, cioè mostrano se il sistema sia rimasto nello stesso attrattore o se sia transitato verso l’altro stato quasi stabile. In sostanza, la figura non descrive soltanto una risposta media, ma la struttura probabilistica della stabilitàdel clima glaciale.
Nei pannelli B, C e D viene analizzata la fase interstadiale, cioè lo stato caldo con AMOC relativamente intensa. Qui la domanda è se l’eruzione possa favorire un collasso dell’AMOC verso lo stato freddo. La risposta della figura è nettamente stato-dipendente: nel pannello B, quando il sistema è ancora relativamente lontano dalla soglia critica, le transizioni vulcaniche sono limitate, circa 8%, mentre quelle dovute al solo rumore interno sono praticamente assenti; nel pannello C, quando il sistema si è avvicinato alla regione di instabilità, la probabilità di collasso sale a circa 60%nel gruppo forzato e a circa 6% nel gruppo non forzato; nel pannello D, in prossimità del tipping point, si arriva a circa 72% di transizioni vulcaniche contro circa 28% di transizioni indotte dal solo rumore. Il messaggio fisico è molto importante: la barriera dinamica che separa lo stato forte e quello debole dell’AMOC si abbassa progressivamente, per cui anche piccole fluttuazioni interne diventano più efficaci, mentre una perturbazione intensa come quella vulcanica risulta ancora più capace di spingere il sistema oltre soglia. Questa lettura è coerente con la letteratura che interpreta gli eventi abrupti glaciali come manifestazioni di un’AMOC soggetta a instabilità, soglie e possibili equilibri multipli; inoltre è in linea con le evidenze paleoceanografiche che collegano i grandi eventi stadiali a un indebolimento persistente della circolazione atlantica e con studi che attribuiscono un ruolo importante al ghiaccio marino e al controllo dell’Oceano Meridionale sulla stabilità dell’AMOC glaciale.
I pannelli F, G e H mostrano il problema opposto, cioè la fase stadiale, quando il sistema si trova nello stato freddo con AMOC debole e l’interesse è capire se l’eruzione possa favorire un recupero verso lo stato caldo. Anche qui la probabilità di transizione dipende in modo marcato dalla posizione del sistema rispetto alla soglia. Nel pannello F, all’inizio della fase stadiale analizzata, il passaggio verso AMOC forte resta raro, circa 8% sia nell’ensemble vulcanico sia in quello non forzato. Nel pannello G, quando il sistema si trova in una regione più vulnerabile, le transizioni salgono a circa 50% nel gruppo vulcanico e 46% in quello con solo rumore. Nel pannello H, infine, quando il sistema è ormai molto vicino al punto di recupero, le transizioni arrivano a circa 94% con eruzione e a circa 60% senza eruzione. Questo comportamento indica che anche la fase stadiale possiede una propria dinamica di soglia, ma suggerisce anche che, in certi momenti, il solo rumore interno possa quasi uguagliare l’efficacia della perturbazione vulcanica. Gli autori collegano questo comportamento a una struttura a due stati quasi stabili e sottolineano inoltre che, durante il passaggio dal freddo al caldo, compare un marcato overshoot dell’AMOC, cioè un temporaneo superamento dei valori medi del ramo forte, elemento che caratterizza la dinamica della ripresa.
Dal punto di vista paleoclimatico, la figura 3 è molto rilevante perché offre un meccanismo fisico e statistico per interpretare risultati emersi anche dalle carote di ghiaccio. Lohmann e Svensson hanno mostrato che, nel periodo 12–60 ka, le grandi eruzioni bipolari compaiono più spesso del caso immediatamente prima dell’inizio di alcuni eventi Dansgaard–Oeschger: 5 eventi su 20 ricadono entro 20 anni da una grande eruzione e 7 su 20 entro 50 anni, suggerendo che il vulcanismo possa agire come fattore di innesco a breve termine quando il sistema è già vicino a una soglia. Un successivo studio del 2024 ha inoltre evidenziato che l’impatto climatico delle grandi eruzioni glaciali è stato-dipendente, cioè varia in funzione del contesto climatico di fondo, un risultato che si adatta molto bene alla logica mostrata in figura 3. Per questo la figura non va letta come dimostrazione che i vulcani “causino” tutti gli eventi D-O, ma come evidenza che il forcing vulcanico può modulare il timing delle transizioni in un sistema climatico non lineare già predisposto all’instabilità. Va anche ricordato che la letteratura non converge su un’unica formalizzazione teorica: accanto ai modelli bistabili esistono lavori recenti che interpretano i cambiamenti abrupti glaciali come fenomeni multiscala di tipo eccitabile. Tuttavia, anche in questa prospettiva alternativa, il risultato centrale della figura resta robusto: la risposta del clima glaciale dipende dall’interazione tra forcing impulsivo, rumore interno, ghiaccio marino, convezione oceanica e vicinanza a un tipping point.
Teleconnessioni post-vulcaniche e rafforzamento dell’AMOC
In questa sottosezione il punto decisivo è che la risposta post-vulcanica dell’AMOC viene interpretata dentro un contesto glaciale e non come una semplice estensione dei meccanismi osservati nel clima moderno. Le eruzioni impiegate da Vettoretti et al. sono infatti molto più intense di quelle documentate nell’epoca strumentale, e agiscono in un sistema in cui grandi calotte continentali, albedo più elevata, ghiaccio marino più esteso e diversa struttura della circolazione atmosferica modificano profondamente la trasmissione del segnale vulcanico. La base fisica resta il classico deficit radiativo da aerosol solfatici stratosferici, che induce un raffreddamento globale iniziale e altera il ciclo idrologico, ma sia i dati da carote di ghiaccio sia i modelli recenti mostrano che l’impatto delle grandi eruzioni durante il clima glaciale fu fortemente stato-dipendente, con raffreddamenti ad alte latitudini dell’emisfero nord più pronunciati rispetto a quelli del periodo recente e con una maggiore capacità del sistema nord-atlantico di amplificare il segnale tramite feedback tra oceano e ghiaccio marino. In questo senso, la sottosezione si colloca bene accanto a lavori che hanno mostrato come le eruzioni maggiori del tardo Pleistocene compaiano più spesso del caso poco prima dell’avvio di alcuni eventi Dansgaard–Oeschger e come, in generale, i cluster eruttivi possano prolungare anomalie fredde attraverso feedback graduali oceano–ghiaccio marino.
L’elemento più raffinato dell’interpretazione proposta dagli autori è la distinzione tra rafforzamento iniziale e indebolimento successivo dell’AMOC. Il testo afferma che subito dopo l’eruzione la circolazione meridionale atlantica si intensifica, e collega questa risposta a una perdita di galleggiamento invernale nelle alte latitudini del Nord Atlantico, associata a anomalie positive di salinità superficiale: in termini fisici, acque superficiali più saline e più dense favoriscono convezione profonda e ventilazione delle regioni di formazione della NADW. Questa idea è coerente con una parte importante della letteratura sul post-vulcanico nord-atlantico: già Mignot et al. avevano mostrato che, nel modello del clima dell’ultimo millennio, dopo alcune grandi eruzioni l’AMOC può intensificarsi dopo circa 5–10 anni prima di entrare in una fase di lieve indebolimento nel decennio seguente, mentre studi più recenti indicano che un rafforzamento iniziale della circolazione del subpolar gyre, l’aumento del trasporto di calore verso nord e l’incremento dei flussi di galleggiamento nelle zone convettive possono precedere una risposta più lenta e complessa dell’AMOC. Nella fase successiva, però, il segnale cambia segno: la crescita del ghiaccio marino, lo sviluppo di una aloclina più marcata, la limitazione degli scambi aria-mare e le modifiche dello stress del vento sul vortice subpolare aumentano il galleggiamento superficiale nelle zone di formazione delle acque profonde e favoriscono la tendenza opposta, cioè l’indebolimento della convezione e quindi dell’AMOC su scala decadale. Proprio questa sequenza bifasica, prima densificazione e poi stratificazione, è ciò che rende il meccanismo descritto nella sottosezione molto più interessante di un semplice “raffreddamento da vulcano”.
Il richiamo alla NAO negativa va poi letto con cautela, come fanno gli stessi autori: non come catena causale lineare e necessaria, ma come possibile espressione della riorganizzazione dei modi interni di variabilità atmosferica dopo il raffreddamento vulcanico. La letteratura moderna mostra che la relazione tra vulcani e NAO non è univoca, perché dipende da latitudine dell’eruzione, stagione, stato del vortice polare, condizioni iniziali di ENSO e struttura media del Nord Atlantico; lavori recenti trovano ad esempio risposte invernali atlantiche differenti tra eruzioni tropicali ed extratropicali dell’emisfero nord, mentre gli studi sul legame NAO–AMOC mostrano che la NAO è effettivamente in grado di modulare la circolazione overturning attraverso anomalie di flusso di calore, acqua dolce e momento. Inoltre, nei modelli CMIP6 il segnale AMOC associato alla NAO risulta molto sensibile allo stato medio subpolare: i modelli più caldo-salini generano variazioni di densità subsuperficiale più forti e dunque una risposta AMOC ritardata più intensa. In altre parole, quando Vettoretti et al. collegano l’anomalia iniziale di salinità superficiale e il riassetto atmosferico a una NAO negativa, stanno inserendo il loro risultato dentro una cornice fisica plausibile ma deliberatamente non deterministica, nella quale il forcing vulcanico agisce soprattutto come perturbazione capace di riorganizzare un sistema già fortemente non lineare. Su questo sfondo diventano allora centrali, nelle pagine che seguono, la scomposizione dei flussi superficiali termici e salini e il ruolo del wind stress, perché è lì che il forcing radiativo viene trasformato in un’anomalia di galleggiamento e quindi in una risposta della circolazione termoalina. Con il proseguimento del testo, il quadro fisico delineato da Vettoretti et al. acquista una struttura molto più chiara: il forcing vulcanico non viene più letto soltanto come un impulso radiativo negativo, ma come un disturbo capace di riorganizzare in modo stagionalmente selettivo i flussi di galleggiamento superficiale, la circolazione atmosferica nord-atlantica, il trasporto salino del subpolar gyre e, infine, la stabilità della convezione profonda. L’aspetto più importante è che il rafforzamento iniziale dell’AMOC coincide con una perdita di galleggiamento alle alte latitudini del Nord Atlantico, soprattutto in inverno, cioè proprio nella stagione in cui la formazione di acque profonde contribuisce in modo decisivo alla circolazione overturning. In questa fase precoce, il termine termico domina la risposta invernale, mentre in estate la perdita di galleggiamento appare più estesa e più fortemente controllata dal contributo alino. Ne emerge una risposta bifasica: dapprima il sistema viene spinto verso una maggiore densificazione superficiale e un rafforzamento della convezione, poi, nel giro di alcuni anni, subentrano processi di segno opposto che favoriscono una progressiva stabilizzazione della colonna d’acqua e l’indebolimento dell’AMOC. Questa lettura è coerente con il risultato generale dello studio, secondo cui grandi eruzioni equatoriali possono modificare l’AMOC tramite cambiamenti accoppiati della circolazione atmosferica e oceanica e dei flussi aria-mare di galleggiamento, spingendo il sistema glaciale verso stati persistenti diversi quando esso è già vicino a una soglia dinamica.
Il passaggio successivo del meccanismo riguarda la riorganizzazione delle teleconnessioni atmosferiche. Gli autori mostrano che nei primi inverni post-eruttivi compare una configurazione assimilabile a El Niño sul Pacifico tropicale, innestata però su un raffreddamento globale di fondo, e che nell’emisfero nord emerge una risposta di tipo Pacific–North American accompagnata, sul settore atlantico, da un indebolimento dell’Icelandic Low e dell’Azores High, quindi da una configurazione coerente con NAO negativa. Il significato dinamico di questa risposta è notevole: la riduzione del gradiente meridionale di pressione attenua i venti occidentali e riduce il trasporto di umidità verso le alte latitudini del Nord Atlantico, con una conseguente diminuzione dell’apporto di acqua dolce nelle aree chiave di formazione della NADW. In questo modo aumenta la salinità superficiale del mare nei mari d’Irminger e di Norvegia, cioè proprio nelle zone convettive più sensibili, e l’anomalia salina positiva favorisce la perdita di galleggiamento e il rafforzamento iniziale dell’AMOC. Questo schema è in buona continuità con studi precedenti sul forcing vulcanico nel Nord Atlantico, che avevano già evidenziato come la risposta oceanica dipenda da aggiustamenti atmosferici rapidi, da anomalie di vento e da modifiche del bilancio di acqua dolce, pur mostrando una forte dipendenza dal contesto climatico e dal modello utilizzato.
Un ulteriore elemento decisivo è il ruolo del subpolar gyre, che in questa interpretazione non rappresenta un semplice dettaglio regionale, ma uno dei principali mediatori tra forcing vulcanico e circolazione overturning. Le anomalie concomitanti di pressione al livello del mare e di stress del vento rafforzano infatti la circolazione del vortice subpolare proprio nella fase in cui l’AMOC si intensifica. Ciò implica un maggiore trasporto verso nord di acque subtropicali relativamente più saline, convogliate dal sistema della Gulf Stream verso le regioni convettive prossime alla Groenlandia, con un ulteriore rinforzo della salinità superficiale e quindi della densità. In termini dinamici, questo configura un feedback positivo temporaneo tra subpolar gyre e AMOC. La letteratura recente supporta fortemente questo punto: lavori dedicati all’accoppiamento tra gyre subpolare e overturning mostrano che cambiamenti autosostenuti nella circolazione del Nord Atlantico possono produrre transizioni abrupthe nella temperatura groenlandese, mentre studi più recenti sul vulcanismo in climi glaciali sottolineano che la risposta del Nord Atlantico dipende in modo cruciale dallo stato medio dell’oceano superiore, dall’estensione del ghiaccio marino e dalla localizzazione delle aree convettive. In altre parole, il forcing vulcanico non opera su un fondo neutro, ma su una struttura oceanica già predisposta ad amplificare oppure smorzare la perturbazione.
La parte forse più interessante del brano è però quella che descrive il cambio di segno della risposta nel decennio successivo. Dopo la fase iniziale di salinizzazione e densificazione, compaiono anomalie di diminuzione della salinità superficiale e un aumento del galleggiamento superficiale nelle aree di formazione delle acque profonde, con contributi termici invernali e alini estivi che agiscono nella direzione opposta rispetto all’avvio del processo. Gli autori legano questo ribaltamento alla crescita del ghiaccio marino alle alte latitudini, indotta dal forte raffreddamento post-vulcanico. L’espansione del ghiaccio sopra le principali regioni convettive favorisce una stratificazione più stabile, mentre la fusione stagionale estiva tende a rinfrescare lo strato superficiale e a consolidare una aloclina persistente. Il risultato è la soppressione della convezione profonda e l’avvio di un feedback ghiaccio marino–NADW che spinge l’AMOC al di sotto dei livelli dell’ensemble di controllo non vulcanico. Questo aspetto coincide molto bene con quanto emerso in studi recenti sul carattere stato-dipendente della risposta vulcanica nel Nord Atlantico: Dutta et al. mostrano che in condizioni glaciali il raffreddamento nord-atlantico può essere più intenso e più persistente proprio perché l’aumento del galleggiamento superficiale e dell’estensione del ghiaccio marino indebolisce l’AMOC, ma anche che l’ampiezza della risposta dipende fortemente dalla stabilità dell’oceano superiore e dai bias con cui i modelli rappresentano convezione e criosfera.
Nel complesso, il proseguimento del testo consente di vedere con grande precisione la logica completa della sequenza post-vulcanica. Prima, il raffreddamento globale riorganizza le teleconnessioni tropicali ed extratropicali e modifica i campi di pressione nord-atlantici; poi, la riduzione del trasporto di umidità verso le alte latitudini e il rafforzamento del subpolar gyre aumentano la salinità nelle regioni di formazione delle acque profonde e rafforzano temporaneamente l’AMOC; successivamente, l’espansione del ghiaccio marino, la freschizzazione stagionale e la formazione di una forte aloclina inibiscono la convezione e avviano l’indebolimento della circolazione overturning. Questa sequenza è particolarmente convincente perché collega in modo continuo forcing radiativo, dinamica atmosferica, bilancio di acqua dolce, trasporto salino, feedback del ghiaccio marino e stabilità della NADW. Inoltre si inserisce bene nel quadro paleoclimatico più ampio: le evidenze da carote di ghiaccio indicano che grandi eruzioni bipolari si verificano più spesso del caso prima dell’avvio di alcuni eventi Dansgaard–Oeschger, mentre la modellistica suggerisce che il vulcanismo, da solo, non sempre basta a produrre un raffreddamento millenario stabile, ma può agire come trigger efficace se il sistema è già vicino a una soglia di instabilità. È proprio questa interazione tra perturbazione impulsiva e vulnerabilità interna del sistema glaciale a rendere il passaggio particolarmente rilevante anche dal punto di vista teorico.
Figura 4
La figura 4 è centrale perché traduce il segnale post-vulcanico in termini di flussi di galleggiamento superficiale invernali sul Nord Atlantico, cioè nella grandezza fisica che più direttamente controlla la densità delle acque superficiali, la convezione profonda e quindi l’intensità della formazione di North Atlantic Deep Water e dell’AMOC. I tre pannelli di colonna separano il contributo totale, quello alino e quello termico, mentre le tre righe mostrano l’evoluzione temporale dopo l’eruzione: anni 7701–7705, 7706–7715 e 7716–7720. Nella convenzione adottata dagli autori, anomalie positive indicano un oceano superficiale più galleggiante e quindi più resistente alla formazione di acque profonde, mentre anomalie negative indicano rimozione di galleggiamento, cioè una superficie più densa e più favorevole alla convezione. Letta in questo modo, la figura mostra una risposta chiaramente bifasica: una prima fase in cui il Nord Atlantico perde galleggiamento e tende a rafforzare la convezione, seguita da una fase in cui il galleggiamento aumenta e la colonna d’acqua diventa progressivamente più stabile. Questa sequenza è esattamente coerente con il risultato generale di Vettoretti et al., secondo cui le grandi eruzioni equatoriali, in clima glaciale, possono prima rafforzare temporaneamente l’AMOC e poi favorirne l’indebolimento nel giro di uno o due decenni.
Nella prima riga, relativa agli anni immediatamente successivi all’eruzione, il segnale dominante è una anomalia negativa del flusso totale di galleggiamento sulle alte latitudini del Nord Atlantico, in particolare tra Mare di Norvegia, settore groenlandese e aree convettive limitrofe. Il fatto più importante è che questa risposta iniziale è dominata soprattutto dal termine termico del pannello C, mentre il contributo alino del pannello B appare più localizzato e meno energico. In termini fisici, ciò indica che nei primi inverni post-eruttivi il raffreddamento e l’aumento della perdita di calore oceano-atmosfera rendono le acque superficiali più dense, favorendo la convezione profonda proprio nella stagione in cui l’AMOC è più sensibile ai flussi di galleggiamento. Questo quadro è in linea con studi precedenti sul forcing vulcanico nord-atlantico, che avevano già evidenziato come alcune grandi eruzioni possano produrre, su scala pluriennale, una riorganizzazione dinamica dell’oceano con temporaneo rafforzamento della circolazione overturning o delle sue componenti subpolari, anche se con forte dipendenza dal contesto climatico e dal modello impiegato.
La seconda e la terza riga mostrano però che questa risposta non persiste. Tra il periodo 7706–7715 e quello 7716–7720 emerge infatti una vasta area di anomalie positive del flusso di galleggiamento sulle regioni convettive nord-atlantiche e nordiche, soprattutto nel pannello totale e in quello termico, segnalando che la superficie oceanica diventa più leggera e più resistente al rimescolamento verticale. In parallelo, il contributo alino indica una crescente importanza dei processi di freschizzazione e di stratificazione. L’interpretazione proposta dagli autori è che, dopo la fase iniziale di densificazione, il forte raffreddamento ad alte latitudini favorisca l’espansione del ghiaccio marino, la successiva fusione stagionale superficiale e la formazione di una aloclina persistente, tutti processi che riducono l’efficienza della convezione profonda. È un risultato molto robusto anche nel confronto con la letteratura recente: Dutta et al. hanno mostrato che, in condizioni glaciali, l’aumento graduale del flusso di galleggiamento nella zona convettiva porta a un indebolimento di lungo periodo dell’AMOC, e lavori precedenti di Vettoretti e di altri autori avevano già attribuito al feedback tra ghiaccio marino, stratificazione e formazione di NADW un ruolo chiave nelle transizioni abrupthe del Nord Atlantico glaciale.
In una lettura più ampia, la figura 4 suggerisce che il vulcano non agisca come semplice agente di raffreddamento uniforme, ma come perturbazione capace di riorganizzare la dinamica dell’Atlantico subpolare. Nella fase iniziale il sistema perde galleggiamento e la convezione si intensifica; successivamente, il rafforzamento del ghiaccio marino e della stratificazione superficiale inverte il segno della risposta. Questo è esattamente il tipo di comportamento atteso in un sistema vicino a una soglia, dove piccoli cambiamenti nei flussi superficiali possono alterare il bilancio tra gyre subpolare, trasporto salino e overturning. Studi recenti sul legame tra subpolar gyre e AMOC mostrano infatti che variazioni nella circolazione subpolare possono contribuire a transizioni rapide del clima groenlandese, mentre le evidenze da carote di ghiaccio indicano che le grandi eruzioni bipolari si verificano più spesso del caso poco prima di alcuni eventi Dansgaard–Oeschger. Per questo la figura 4 non va vista come una semplice mappa di anomalie di flusso, ma come una rappresentazione molto efficace del passaggio da una fase iniziale di rafforzamento convettivo a una successiva fase di stabilizzazione stratificata, cioè del meccanismo fisico che collega il forcing vulcanico al possibile indebolimento multidecadale dell’AMOC in condizioni glaciali.
Figura 5
La figura 5 mostra le anomalie estive del flusso di galleggiamento superficiale nel Nord Atlantico dopo l’eruzione, separate in contributo totale, alino e termico, e va letta come il complemento stagionale della figura 4. Qui il punto chiave è che, in estate, il sistema non risponde semplicemente con un raffreddamento uniforme, ma con una riorganizzazione dei processi che controllano densità, stratificazione e predisposizione alla convezione profonda. Nei pannelli superiori, relativi agli anni 7701–7705, domina un segnale diffusamente negativo del flusso totale di galleggiamento sul Nord Atlantico subpolare, cioè una condizione che rende temporaneamente la superficie oceanica meno galleggiante e quindi più favorevole alla densificazione; tuttavia, gli autori sottolineano che nella stagione estiva questa perdita iniziale di galleggiamento è legata soprattutto alla componente alina, a differenza della risposta invernale che è invece più nettamente controllata dal termine termico. Questo passaggio è importante perché indica che l’eruzione modifica subito anche il bilancio di acqua dolce e la distribuzione superficiale della salinità, non soltanto gli scambi di calore aria-mare. In termini dinamici, la figura suggerisce quindi che il primo effetto estivo del forcing vulcanico è ancora compatibile con il temporaneo rafforzamento dell’AMOC mostrato nella figura 3, ma attraverso una firma stagionale diversa da quella invernale.
La parte più istruttiva della figura emerge però nelle due righe successive, corrispondenti agli intervalli 7706–7710 e 7716–7720, quando nelle regioni convettive del Nord Atlantico e soprattutto nel settore del subpolar gyre compaiono anomalie positive di galleggiamento via via più evidenti. In pratica, la superficie oceanica diventa più leggera e più resistente alla formazione di acque profonde, e questo segnale appare dominato soprattutto dal contributo alino del pannello B, mentre quello termico del pannello C resta più debole e spazialmente secondario. La lettura fisica proposta da Vettoretti et al. è che il forte raffreddamento post-vulcanico favorisca l’espansione del ghiaccio marino alle alte latitudini; successivamente, la fusione stagionale e la riorganizzazione dei flussi di acqua dolce freschizzano lo strato superficiale e contribuiscono alla formazione di una aloclina persistente, che stabilizza la colonna d’acqua e inibisce la convezione. È esattamente questo tipo di meccanismo che la letteratura recente considera cruciale per la stabilità dell’AMOC glaciale: il controllo dei flussi di galleggiamento superficiale nelle regioni subpolari, soprattutto quando ghiaccio marino e freschizzazione modificano la trasformazione di massa d’acqua, può spostare il sistema da una fase più convettiva a una fase più stratificata e meno favorevole alla NADW.
Rispetto alla figura 4, la figura 5 aggiunge quindi un’informazione molto preziosa: mostra che il ramo estivo della risposta post-vulcanica è il luogo in cui si costruisce progressivamente la memoria superficiale del sistema, cioè quella stabilizzazione salina che poi prepara l’indebolimento dell’AMOC su scala decadale. In inverno la dinamica è più direttamente connessa alla perdita di calore e alla convezione immediata; in estate, invece, il sistema accumula gli effetti di ghiaccio marino, fusione, salinità superficiale e redistribuzione dell’acqua dolce. Questa distinzione stagionale è molto coerente con gli studi che attribuiscono ai feedback tra ghiaccio marino e overturning un ruolo fondamentale nelle transizioni del Nord Atlantico glaciale: Ando e colleghi hanno mostrato che l’isteresi dell’AMOC glaciale è fortemente controllata da feedback termici associati al raffreddamento atmosferico e alle transizioni del ghiaccio marino, mentre Prange et al. hanno evidenziato che la variabilità multicentenaria nord-atlantica in clima glaciale viene amplificata proprio dai feedback del ghiaccio marino e dalla distribuzione della salinità. In questa prospettiva, la figura 5 non è soltanto una mappa estiva di anomalie, ma la rappresentazione del passaggio da una risposta inizialmente favorevole alla densificazione a una configurazione sempre più dominata dalla freschizzazione superficiale e dalla stratificazione, cioè dal meccanismo che prepara il successivo declino dell’AMOC.
Nel complesso, la figura 5 conferma che il forcing vulcanico agisce in un sistema glaciale altamente non lineare e fortemente stato-dipendente. La risposta estiva non replica semplicemente quella invernale, ma la completa: dapprima contribuisce alla fase di aggiustamento rapido, poi favorisce l’emergere di un oceano superficiale più galleggiante nelle aree convettive, più stabile e meno capace di sostenere la ventilazione profonda. Questo risultato è in piena sintonia con gli studi più recenti sulle risposte vulcaniche del Nord Atlantico in clima glaciale, che mostrano come l’effetto delle eruzioni dipenda non solo dall’intensità del forcing, ma anche dalla struttura termica dell’oceano superiore, dall’estensione del ghiaccio marino e dalla sensibilità del bilancio di acqua dolce nelle zone di formazione delle acque profonde. Per questo la figura 5 è essenziale: rende visibile il modo in cui la stagione estiva trasforma una perturbazione impulsiva in una anomalia di galleggiamento persistente, capace di modulare la traiettoria successiva dell’AMOC.
Prospettive dei sistemi dinamici sulla variabilità climatica decadale post-vulcanica
In questa parte del lavoro, Vettoretti e colleghi compiono un passaggio concettuale molto importante: dopo aver descritto con un approccio process-based la risposta iniziale del sistema climatico glaciale all’eruzione, spostano l’attenzione sulla scala multidecadale, dove il comportamento dell’ensemble diventa rapidamente divergente e non può più essere interpretato in modo convincente seguendo, membro per membro, una catena causale lineare di processi fisici. Per affrontare questo problema, gli autori adottano il linguaggio dei sistemi dinamici non lineari e introducono un modello ridotto, stocastico, costruito per imitare qualitativamente la dinamica del modello climatico complesso. In pratica, il forcing vulcanico viene rappresentato come una perturbazione breve sovrapposta a equazioni differenziali stocastiche integrate con schema di Euler–Maruyama, mentre, una volta terminato l’impulso, il sistema viene lasciato evolvere liberamente, così da isolare la risposta post-vulcanica dell’AMOC. Il punto metodologico di fondo è molto robusto e si inserisce bene nella tradizione dei modelli concettuali del clima: la letteratura recente sottolinea infatti che, per i problemi di tipping dell’AMOC, i modelli ridotti non sono meri esercizi astratti, ma strumenti interpretativi capaci di mettere a fuoco i feedback dominanti e di costruire un ponte teorico tra box models, modelli di complessità intermedia e GCM accoppiati.
Il risultato più rilevante è che, regolando tre ingredienti essenziali — intensità del rumore interno, ampiezza della perturbazione vulcanica e posizione del punto fisso del sistema, controllata nel modello semplice dal parametro associato alla CO₂ di fondo — gli autori riescono a riprodurre in modo qualitativo molte delle caratteristiche già osservate nel CCSM4: la presenza di due stati quasi stabili dell’AMOC, la possibilità di transizioni indotte dal rumore e il fatto che un impulso vulcanico possa spingere il sistema oltre il confine del bacino di attrazione, favorendo un passaggio temporaneo verso uno stato freddo a bassa AMOC. Questa impostazione è pienamente coerente con una parte importante della letteratura sui cicli Dansgaard–Oeschger, che da tempo interpreta tali eventi come espressione di una dinamica interna del sistema atmosfera–ghiaccio–oceano del Nord Atlantico e dell’Artico, piuttosto che come risposta rigidamente periodica a un forcing esterno. In particolare, Ganopolski e Rahmstorf mostrarono già nel 2002 che la risonanza stocastica poteva costituire un meccanismo plausibile per la variabilità millenaria glaciale, mentre la review di Li e colleghi del 2019 ha ribadito che gli eventi D-O possono essere letti come oscillazioni spontanee del sistema accoppiato delle Northern Seas. Il modello semplice usato qui si colloca esattamente in questa linea di pensiero: il vulcano non “crea” la dinamica, ma la devia quando il sistema è già predisposto a transizioni tra stati persistenti.
Da questo punto di vista, la sottosezione è particolarmente interessante perché suggerisce che, alle condizioni di fondo imposte nelle simulazioni glaciali con CO₂ fissa a 230 ppm, le transizioni dell’AMOC siano soprattutto noise-induced, cioè indotte dalla variabilità interna, più che da un’instabilità deterministica pienamente sviluppata. Gli autori osservano infatti che il miglior accordo con il modello complesso si ottiene quando il sistema ridotto possiede un punto fisso stabile, corrispondente a uno stato caldo ad AMOC forte, collocato tuttavia vicino a una soglia dinamica. In questo quadro, il sistema può apparire come un regime stabile ma vulnerabile, in cui il rumore interno e l’impulso vulcanico modulano la probabilità di attraversare una frontiera di bacino; tuttavia, poiché nel GCM la variabilità interna non può essere rimossa, resta difficile distinguere in modo netto questo scenario da uno in cui il sistema sia invece già entrato in un regime oscillatorio vicino a una biforcazione di Hopf. Questa prudenza interpretativa è ben fondata. Lavori recenti sul CESM sostengono che il collasso dell’AMOC in quel modello sia compatibile con il quadro classico di una biforcazione saddle-node, mentre la teoria degli early warning signals collega l’avvicinamento a un tipping point alla perdita di resilienza e al critical slowing down; gli stessi autori, però, sottolineano che segnali di questo tipo non sono ancora stati dimostrati in modo convincente nelle loro simulazioni glaciali. Il testo, quindi, non pretende di chiudere il dibattito sulla geometria esatta della biforcazione, ma usa il formalismo dinamico per dare una lettura fisicamente coerente dell’elevata sensibilità del sistema in prossimità della soglia.
Molto raffinata è poi l’analisi del ruolo del rumore. Nel modello semplice senza rumore, la perturbazione vulcanica può spingere con relativa efficienza il sistema oltre il confine tra i bacini di attrazione e indurre la transizione verso il manifold inferiore, cioè verso lo stato freddo con AMOC debole, dal quale il sistema recupera successivamente lo stato caldo perché l’attrattore superiore resta stabile. Quando però si aumenta progressivamente il livello di rumore, la dinamica cambia: il comportamento diventa meno deterministico, la variabilità interna pesa sempre di più e, a rumore elevato, la perturbazione vulcanica può persino diventare meno efficace nel produrre la transizione, perché il sistema tende a restare nel manifold superiore. Questo è un risultato molto interessante, poiché spiega in modo elegante perché nel CCSM4 alcune realizzazioni attraversino effettivamente la soglia mentre altre, pur sottoposte allo stesso forcing vulcanico, recuperino oppure non completino il passaggio stadiale. In altre parole, il rumore non è soltanto il fattore che rende possibili le transizioni: può anche renderle meno prevedibili e perfino ostacolarle. Tale visione è compatibile sia con gli studi che attribuiscono alla posizione media del margine del ghiaccio marino un controllo decisivo sulla stabilità dell’AMOC glaciale, sia con sviluppi più recenti che propongono per gli eventi D-O una dinamica multiscala di tipo eccitabile, non riducibile in modo univoco a una semplice bistabilità classica. Per questo il valore della sottosezione non sta solo nel mostrare che un modello ridotto può imitare il GCM, ma nel chiarire che la risposta climatica post-vulcanica va interpretata come il prodotto congiunto di forzante impulsiva, vicinanza a una biforcazione, struttura del bacino di attrazione e variabilità interna.
Figura 6
La figura 6 rappresenta uno dei passaggi più importanti dell’intero studio, perché ricostruisce in modo organico la sequenza delle risposte climatiche del Nord Atlantico dopo l’eruzione vulcanica, mettendo in relazione atmosfera, oceano superficiale, ghiaccio marino, circolazione subpolare e intensità dell’AMOC. Le tre righe mostrano l’evoluzione temporale della risposta nei diversi intervalli successivi all’eruzione, mentre le colonne illustrano, rispettivamente, le anomalie della temperatura dell’aria e della pressione al livello del mare, della salinità superficiale e della profondità dello strato rimescolato, del bilancio superficiale dell’acqua dolce e dell’estensione del ghiaccio marino, della circolazione barotropica e dello stress del vento, e infine della stessa AMOC. Considerata nel suo insieme, la figura mostra con notevole chiarezza che il forcing vulcanico non produce soltanto un raffreddamento generalizzato, ma innesca una vera riorganizzazione del sistema nord-atlantico, nella quale i cambiamenti atmosferici iniziali vengono rapidamente trasferiti all’oceano attraverso modifiche della salinità, del trasporto di umidità, del vento e della stabilità della colonna d’acqua.
Nel primo intervallo temporale successivo all’eruzione emerge un raffreddamento diffuso sul Nord Atlantico e sulle regioni circostanti, accompagnato da una configurazione della pressione atmosferica coerente con una fase di NAO negativa. Questo assetto implica un indebolimento del gradiente barico meridionale e una riduzione dei venti occidentali che normalmente convogliano umidità verso le alte latitudini atlantiche. La conseguenza è una diminuzione dell’apporto di acqua dolce nelle principali regioni di formazione delle acque profonde, in particolare nel settore compreso tra il Mare d’Irminger, il Mare di Norvegia e le aree limitrofe. La risposta oceanica immediata è un aumento della salinità superficiale proprio in queste zone, accompagnato da un approfondimento dello strato rimescolato. Questo significa che le acque superficiali diventano più dense e più favorevoli alla convezione profonda. Nello stesso tempo, la figura mostra anche un rafforzamento del vortice subpolare, sostenuto dalle anomalie di vento e di circolazione barotropica, che tende a incrementare il trasporto verso nord di acque relativamente più saline provenienti da latitudini inferiori. L’effetto combinato di questi processi si riflette nell’ultima colonna, dove compare una temporanea intensificazione dell’AMOC. In sostanza, la prima fase post-vulcanica è dominata da una risposta che favorisce la densificazione superficiale e la ventilazione profonda del Nord Atlantico, un risultato coerente con altri studi che hanno evidenziato come grandi eruzioni possano produrre, almeno inizialmente, un rafforzamento della circolazione overturning attraverso modifiche accoppiate della circolazione atmosferica e dei campi di salinità. (Vettoretti et al., 2026; Lohmann e Svensson, 2022)
La parte più interessante della figura emerge però negli intervalli successivi, quando il segnale cambia progressivamente segno. Nel secondo periodo, e ancora più chiaramente nel terzo, il raffreddamento alle alte latitudini favorisce un’espansione del ghiaccio marino nelle regioni nord-atlantiche e nordiche. Questo passaggio è cruciale, perché il ghiaccio marino non agisce solo come indicatore termico, ma come elemento dinamico capace di modificare in profondità gli scambi tra oceano e atmosfera e di alterare la struttura verticale dell’oceano superiore. Con il passare degli anni, nelle aree convettive principali si osserva infatti una riduzione della salinità superficiale, associata a una minore profondità dello strato rimescolato e a una crescente stabilizzazione della colonna d’acqua. In altre parole, il sistema evolve da una fase iniziale favorevole alla convezione verso una configurazione in cui l’oceano superficiale diventa più stratificato e meno predisposto al rimescolamento verticale. Una parte di questo cambiamento è legata al fatto che la crescita del ghiaccio marino e la successiva fusione stagionale contribuiscono a rendere più fresco e meno salino lo strato superficiale, favorendo la formazione di una struttura aloclina più robusta. Questo processo riduce l’efficienza della formazione di acque profonde e indebolisce progressivamente il legame tra superficie e oceano intermedio-profondo.
La quarta colonna della figura mostra che anche il contributo dinamico del vortice subpolare si attenua nelle fasi successive rispetto alla risposta iniziale. Il rafforzamento precoce della circolazione girometrica, che nei primi anni aveva contribuito a trasportare verso nord acque più saline, lascia il posto a un assetto meno efficace nel sostenere la densificazione convettiva. Di conseguenza, l’ultima colonna evidenzia il passaggio fondamentale dell’intera sequenza post-vulcanica: dopo un primo rafforzamento, l’AMOC inizia a indebolirsi e, nel periodo più avanzato mostrato nella figura, presenta anomalie negative nelle fasce latitudinali chiave. Questo comportamento è pienamente coerente con la dinamica descritta nel resto dell’articolo e con una parte importante della letteratura paleoclimatica e modellistica, secondo cui il Nord Atlantico glaciale risponde alle perturbazioni impulsive attraverso una sequenza bifasica: prima una fase di maggiore densificazione e rafforzamento convettivo, poi una fase dominata da stratificazione, crescita del ghiaccio marino e progressiva soppressione della convezione. Studi recenti hanno mostrato infatti che, in condizioni glaciali, il ruolo del ghiaccio marino e del bilancio superficiale dell’acqua dolce è determinante per la stabilità dell’AMOC e per la possibilità che una perturbazione relativamente breve produca conseguenze climatiche di durata molto più lunga. (Dutta et al., 2025; Buizert e Schmittner, 2015)
Nel complesso, la figura 6 ha un valore particolarmente elevato perché non si limita a mostrare una serie di anomalie spaziali, ma rende visibile il meccanismo fisico che collega il forcing vulcanico alla risposta decadale della circolazione atlantica. La sequenza che emerge è chiara: inizialmente il raffreddamento globale riorganizza la circolazione atmosferica e riduce il trasporto di umidità verso le alte latitudini; ciò aumenta la salinità superficiale nelle aree di formazione delle acque profonde e rafforza temporaneamente sia il vortice subpolare sia l’AMOC. In una fase successiva, però, il persistente raffreddamento alle alte latitudini promuove la crescita del ghiaccio marino, la stabilizzazione della colonna d’acqua e la riduzione della convezione, aprendo la strada a un indebolimento della circolazione overturning. La figura 6 rappresenta quindi il ponte interpretativo tra le anomalie di galleggiamento mostrate nelle figure precedenti e le probabilità di transizione dell’AMOC discusse in figura 3, perché mostra in modo esplicito come una perturbazione vulcanica possa spingere il sistema nord-atlantico glaciale prima verso una maggiore instabilità convettiva e poi verso una configurazione più stratificata e vulnerabile a una transizione di regime.
Figura 7
La figura 7 è fondamentale perché mostra come gli autori passino da una descrizione process-based della risposta post-vulcanica a una lettura esplicitamente dinamica del sistema climatico glaciale. Invece di seguire, per ogni membro dell’ensemble del modello complesso, l’intera catena di cause fisiche che porta o meno alla transizione dell’AMOC, Vettoretti et al. costruiscono un modello dinamico semplificato capace di riprodurre qualitativamente la stessa logica: uno stato superiore con AMOC forte, uno stato inferiore con AMOC debole, una regione di soglia che separa i due regimi, e un forcing vulcanico breve che può spingere il sistema oltre quel limite. Nel pannello A, la traiettoria blu rappresenta il controllo senza vulcanismo; le linee tratteggiate indicano le soglie dinamiche, mentre la fascia verde identifica il confine del bacino di attrazione, cioè la regione oltre la quale il sistema tende a cambiare ramo evolutivo. I pannelli B, C e D mostrano poi cosa accade quando l’eruzione viene applicata in tre momenti diversi dell’interstadiale, sempre più vicini al tipping point: la probabilità di transizione vulcanica cresce dal 12% nel pannello B al 66% nel pannello C fino all’88% nel pannello D, mentre anche le transizioni indotte dal solo rumore aumentano da 0% a 8% e poi a 18%. Il significato fisico è molto chiaro: più il sistema si avvicina alla soglia, più diventa vulnerabile sia alla perturbazione vulcanica sia alla variabilità interna. Questo risultato è del tutto coerente con l’idea, già classica nella letteratura sui cicli Dansgaard–Oeschger, che il Nord Atlantico glaciale possa funzionare come un sistema non lineare con stati persistenti multipli o quasi-multipli, nel quale forzanti impulsive e rumore interno modulano la probabilità di salto tra regimi climatici.
La metà destra della figura rende il messaggio ancora più raffinato, perché isola il ruolo del rumore interno. Nel pannello E viene fissato un unico intervallo temporale e, nei pannelli F, G e H, si aumenta progressivamente il livello di rumore del modello. Quando il rumore è nullo, nel pannello F, la risposta è quasi perfettamente deterministica: il forcing vulcanico induce la transizione nel 100% dei casi, mentre in assenza di eruzione non si verifica alcuna transizione. Quando il rumore cresce, nel pannello G, la probabilità di transizione vulcanica scende all’80%; nel pannello H, con rumore ancora più elevato, scende al 74%, mentre compaiono anche poche transizioni spontanee nel gruppo non vulcanico, pari al 4%. Questo è uno dei risultati più interessanti dell’intera figura, perché mostra che il rumore non agisce soltanto come fattore che “aggiunge instabilità”, ma può anche ridurre l’efficacia della perturbazione vulcanica, deviando le traiettorie e mantenendo una parte dell’ensemble nel ramo superiore con AMOC forte. In altre parole, la variabilità interna non è un semplice contorno statistico: è una parte costitutiva della dinamica, capace sia di favorire sia di ostacolare il tipping. Questa interpretazione si inserisce molto bene nella tradizione della stochastic resonance e, più in generale, negli approcci che leggono gli eventi abrupti glaciali come emergenti dall’interazione tra forcing esterni deboli o intermittenti e variabilità interna del sistema oceano-atmosfera-ghiaccio marino.
Dal punto di vista teorico, la figura 7 suggerisce che il modello complesso e quello ridotto condividano una stessa geometria di fondo: un sistema con un attrattore caldo relativamente stabile, ma collocato abbastanza vicino a una soglia da rendere possibili transizioni premature verso uno stato freddo. Vettoretti et al. spiegano infatti che il miglior accordo con il modello complesso si ottiene quando il modello semplice viene configurato con un punto fisso stabilecorrispondente a un’AMOC forte, mentre le transizioni risultano soprattutto noise-induced alle condizioni di anidride carbonica di fondo usate negli esperimenti. La figura, quindi, non dimostra solo che un modello semplice “imita” il GCM, ma mostra qualcosa di più importante: che il comportamento post-vulcanico dell’AMOC può essere compreso come il risultato di una competizione tra distanza dalla biforcazione, ampiezza del forcing e livello di rumore interno. Questo quadro è in linea sia con la review di Li e Born, che interpreta gli eventi D-O come oscillazioni spontanee del sistema accoppiato delle Northern Seas condizionate da ghiaccio marino, venti e subpolar gyre, sia con studi più recenti che identificano vere e proprie soglie termiche dell’AMOC glaciale e valorizzano l’uso di emulatori e modelli ridotti per interpretare i tipping points della circolazione atlantica. Perciò la figura 7 va letta come il ponte concettuale tra la fisica dettagliata mostrata nelle figure precedenti e la teoria della variabilità abrupta: il vulcano fornisce l’impulso, ma l’esito finale dipende dalla struttura del paesaggio dinamico in cui quel sistema si trova in quel momento.
Figura 8
La figura 8 costituisce la sintesi concettuale dell’intero lavoro, perché traduce i risultati del modello complesso e del modello ridotto in un unico schema dinamico del ciclo Dansgaard–Oeschger. L’asse verticale rappresenta la risposta rapida del sistema, identificata con la forza dell’AMOC, mentre l’asse orizzontale rappresenta un forcing lento legato all’Oceano Meridionale; su questo piano gli autori collocano un ramo superiore associato a un clima interstadiale caldo con AMOC forte, un ramo inferiore associato a un clima stadiale freddo con AMOC debole e una struttura intermedia instabile che separa i due regimi. Nel quadro proposto per condizioni glaciali con anidride carbonica attorno a 230 ppm, il sistema si troverebbe, in assenza di rumore interno, vicino a un punto fisso stabile sul ramo superiore, cioè in uno stato caldo e non oscillatorio. Questa lettura è coerente con quanto Vettoretti et al. scrivono esplicitamente nella descrizione della figura e si inserisce bene nella letteratura che interpreta gli eventi D-O come espressione di una dinamica interna del sistema accoppiato atmosfera–ghiaccio–oceano del Nord Atlantico, piuttosto che come semplice risposta lineare a un forcing esterno. Anche la review di Li e Born sottolinea infatti che tali oscillazioni possono emergere da fisica nota, purché il clima di fondo sia sufficientemente diverso da quello interglaciale, mentre gli esperimenti di isteresi di Ando e Oka mostrano che l’AMOC glaciale può esibire equilibri multipli e transizioni di regime controllate da feedback termici e dal ghiaccio marino.
Il nucleo teorico della figura è il rapporto tra punti sella, biforcazioni di Hopf e perdita di stabilità del regime forte. Gli autori propongono che, se il parametro di controllo associato alla concentrazione di anidride carbonica si sposta dalla configurazione “alta” verso quella intermedia, il sistema attraversi una biforcazione di Hopf in prossimità del punto sella e possa allora entrare in un regime oscillatorio, cioè in un vero e proprio oscillatore D-O tra stato forte e stato debole dell’AMOC. È importante notare che nel testo viene anche precisato che, nelle simulazioni glaciali considerate, segnali classici di critical slowing down o altri precursori di instabilità non sono ancora stati dimostrati in modo netto, per cui la figura 8 va letta come uno schema interpretativo forte ma non come una dimostrazione definitiva della geometria della biforcazione. Proprio per questo il suo valore è notevole: essa mostra come un sistema climatico ad altissima dimensionalità possa comunque essere interpretato attraverso una struttura a bassa dimensionalità, in cui pochi feedback dominanti organizzano la dinamica di transizione. Questo è perfettamente in linea con lavori recenti sui tipping dell’AMOC, nei quali si sostiene che anche modelli climatici molto complessi possano mostrare una struttura riconducibile a una biforcazione saddle-node, e con riflessioni metodologiche più ampie secondo cui i modelli concettuali restano strumenti essenziali per comprendere la dinamica dei tipping elements climatici.
La traiettoria azzurra della figura rappresenta la possibilità che la variabilità interna del sistema generi da sola un ciclo tra i due rami, cioè una sequenza di transizioni noise-induced tra stato caldo e stato freddo. Qui la figura dialoga direttamente con la tradizione della stochastic resonance: già Ganopolski e Rahmstorf mostrarono che la variabilità millenaria glaciale può emergere dall’interazione tra rumore interno e struttura non lineare del sistema climatico, e Vettoretti et al. riprendono questa idea sostenendo che, nelle loro simulazioni, il clima glaciale può essere pensato come un sistema a due stati quasi stabili nel quale il rumore rende possibili transizioni che, in assenza di variabilità interna, sarebbero molto più rare o del tutto assenti. La traiettoria rossa, invece, rappresenta la transizione indotta dal vulcano: una grande eruzione non crea da zero la dinamica del ciclo D-O, ma può spingere il sistema oltre soglia prima che esso raggiunga naturalmente la regione del tipping point, anticipando così il passaggio verso il ramo inferiore a bassa AMOC. Questo è uno dei messaggi più forti della figura: il vulcano agisce come un impulso breve che deforma il timing della transizione in un sistema già predisposto alla vulnerabilità dinamica.
L’aspetto forse più raffinato è però l’ultima implicazione della figura 8: la variabilità interna non solo può favorire la transizione, ma può anche impedirla. Gli autori scrivono infatti che il rumore atmosferico e oceanico può in alcuni casi ostacolare il passaggio verso lo stato debole, mantenendo il sistema sul ramo superiore nonostante l’impulso vulcanico. Questo significa che la risposta post-vulcanica dell’AMOC non è deterministica, bensì probabilistica, e dipende dalla posizione del sistema nel paesaggio dinamico, dall’ampiezza della perturbazione e dall’intensità del rumore interno. In questo senso la figura 8 non è soltanto una spiegazione grafica del paper, ma un vero quadro teorico unificante: collega la bistabilità dell’AMOC, il ruolo dell’Oceano Meridionale come forcing lento, i feedback di ghiaccio marino e convezione, l’azione del vulcanismo come trigger e la variabilità interna come fattore che modula sia l’innesco sia la resilienza del sistema. È anche per questo che essa risulta pienamente coerente con studi recenti sulle isteresi rumorosedell’AMOC, secondo i quali le transizioni diventano più probabili vicino a una biforcazione ma la stessa presenza del rumore può produrre traiettorie molto diverse a parità di forcing esterno.
Discussione
La discussione proposta da Vettoretti et al. è particolarmente rilevante perché sposta il fuoco interpretativo dall’idea di un vulcanismo capace di “produrre” direttamente gli eventi Dansgaard–Oeschger a una visione più sofisticata, nella quale le grandi eruzioni agiscono soprattutto come inneschi condizionati dallo stato di fondo del sistema climatico. In questa prospettiva, configurazione orbitale, concentrazione dei gas serra, volume delle calotte continentali, stato del ghiaccio marino e stabilità dell’AMOC definiscono una sorta di finestra dinamica favorevole, un vero “sweet spot”, entro il quale il sistema glaciale può diventare particolarmente sensibile a perturbazioni impulsive. Questa impostazione è coerente sia con il lavoro stesso, che mostra come il forcing vulcanico aumenti la probabilità di transizione solo quando il sistema è già vicino a una soglia, sia con studi che attribuiscono alla stabilità dell’AMOC glaciale un forte controllo da parte dell’Oceano Meridionale, del margine del ghiaccio marino e della struttura delle acque profonde. Anche la recente letteratura di sintesi sui modelli di MIS 3 sottolinea che il comportamento tipo D-O non emerge in condizioni preindustriali, ma richiede boundary conditions glaciali specifiche, mentre le evidenze da carote di ghiaccio mostrano che l’impatto delle grandi eruzioni nel glaciale fu nettamente state-dependent.
In questo quadro, il richiamo alla maggiore resilienza del clima olocenico rispetto a eruzioni isolate molto grandi è plausibile e scientificamente ben fondato. La semplice assenza, nell’Olocene iniziale e medio, di salti climatici comparabili agli eventi D-O, pur in presenza di forti eruzioni registrate nelle carote di ghiaccio, suggerisce infatti che il sistema fosse collocato al di fuori della finestra di massima suscettibilità che caratterizzò parte del MIS 3. La stessa logica può estendersi al MIS 5a, che rappresentava un contesto meno glaciale del MIS 3, con livello marino più alto e dunque calotte continentali ridotte, cioè con condizioni di fondo verosimilmente meno favorevoli a transizioni abrupthe guidate da un AMOC prossima alla soglia. In questo senso la discussione degli autori si inserisce bene nella più ampia idea che la variabilità millenaria nord-atlantica non dipenda solo dalla perturbazione esterna, ma dal fatto che il sistema si trovi o meno in una geometria dinamica predisposta all’amplificazione.
Molto convincente è anche la proposta di spostare l’attenzione futura dagli eventi estremi isolati ai cluster di eruzioni, perché è proprio la persistenza del forcing, più che il solo picco iniziale, a poter modificare in modo duraturo il bilancio tra oceano, ghiaccio marino e atmosfera. Studi recenti sul clima olocenico mostrano che gruppi ravvicinati di eruzioni possono produrre raffreddamenti prolungati attraverso feedback lenti oceano–ghiaccio marino, mentre simulazioni dedicate al Nord Atlantico suggeriscono che cluster vulcanici possono generare raffreddamenti intensi ma non sempre sufficienti, da soli, a sostenere uno stato di tipo Younger Dryas. Questo punto rafforza la prudenza della discussione: il vulcanismo può essere un attivatore importante della variabilità millenaria, ma il suo effetto dipende dalla durata del forcing, dalla frequenza degli eventi e soprattutto dalla struttura termica e salina dell’oceano superiore al momento dell’impatto. Ne deriva una visione meno semplicistica e molto più realistica, in cui il vulcanismo viene trattato come parte di una dinamica accoppiata e non come una causa unica.
La parte finale della discussione, dedicata alle implicazioni per il clima moderno, è forse la più delicata ma anche la più stimolante. Gli autori suggeriscono che una maggiore variabilità interna possa, in alcuni contesti, aumentare la resilienza del sistema a un collasso abrupto dell’AMOC, perché il rumore non favorisce soltanto il tipping, ma può anche ostacolarlo; tuttavia, trasferire direttamente questa conclusione dal mondo glaciale a quello contemporaneo richiede molta cautela. Sul futuro dell’AMOC, infatti, la letteratura recente non converge su un’unica diagnosi: da un lato, studi basati su segnali precoci e indicatori statistici suggeriscono che il sistema possa avvicinarsi a una transizione critica, con finestre temporali potenzialmente anche entro questo secolo; dall’altro, analisi multi-modello più recenti indicano che, pur in presenza di un marcato indebolimento, un collasso completo entro il XXI secolo potrebbe essere meno probabile di quanto ipotizzato in altri lavori, perché l’upwelling dell’Oceano Meridionale tende a sostenere una AMOC indebolita ma persistente. Resta però un punto fermo: sia l’IPCC sia la gran parte delle simulazioni CMIP6 indicano un indebolimento dell’AMOC nel corso del XXI secolo, e proprio per questo gli analoghi paleoclimatici, pur non essendo trasferibili in modo diretto, restano essenziali per capire tempi, soglie e feedback della circolazione overturning. La forza della discussione sta dunque nell’aver mostrato che il passato glaciale non fornisce una previsione semplice del futuro, ma offre un laboratorio fisico prezioso per comprendere come impulsi brevi, variabilità interna e stati di fondo diversi possano interagire nel portare il sistema verso o lontano da una soglia critica.
0 commenti