Ruza F. Ivanovic 1 , Lauren J. Gregoire 1 , Andrew D. Wickert2 , Paul J. Valdes 3 , and
Andrea Burke 4
1 School of Earth and Environment, University of Leeds, Leeds, UK, 2 Department of Earth Sciences, University of Minnesota,
Twin Cities, Minneapolis, Minnesota, USA, 3
School of Geographical Sciences, University of Bristol, Bristol, UK, 4 Department
of Earth and Environmental Sciences, University of St. Andrews, St. Andrews, UK
In questo studio, Ivanovic e colleghi interpretano il collasso della sella glaciale che univa le calotte Laurentide e Cordilleran come un evento chiave della deglaciazione tardoglaciale, capace di trasformare un rapido impulso di acqua di fusione in una perturbazione climatica emisferica: secondo le loro simulazioni, un rilascio equivalente a circa 7,26 m di innalzamento medio globale del livello del mare in circa 340 anni avrebbe indebolito l’AMOC di circa 6 Sverdrup, favorendo un raffreddamento diffuso dell’emisfero settentrionale compreso fra 1 e 5 °C, con massimi invernali nel settore nord-atlantico e possibili implicazioni per l’interruzione del riscaldamento del Bølling e per l’innesco o il rafforzamento dell’Older Dryas. Il lavoro si inserisce in una linea di ricerca ormai ben consolidata: Gregoire, Payne e Valdes avevano già mostrato che i collassi delle “ice saddles” possono produrre innalzamenti marini deglaciali molto rapidi; Gomez et al. hanno ulteriormente quantificato il contributo specifico della sella Laurentide-Cordilleran al Meltwater Pulse 1A; Gregoire et al. hanno poi evidenziato che il riscaldamento abrupto del Bølling poteva accelerare la fusione nordamericana fino a 3–4 m di sea-level equivalent, amplificati a 5–6 m quando veniva innescato il collasso della sella. In parallelo, Liu et al. e Obase et al. hanno sottolineato quanto la transizione Bølling-Allerød dipendesse criticamente dalla sensibilità dell’AMOC ai flussi di acqua dolce, mentre Ivanovic et al. (2018) e Matero et al. hanno rafforzato l’idea che fusioni accelerate delle calotte possano produrre rallentamenti della circolazione atlantica e raffreddamenti abrupti coerenti con i proxy paleoclimatici.
Introduzione
Nel quadro della deglaciazione che seguì il Last Glacial Maximum, quando le grandi calotte nordamericane iniziarono a perdere massa in risposta all’aumento dell’insolazione estiva boreale, all’incremento della CO₂ atmosferica e alla riorganizzazione della circolazione climatica, il lavoro di Ivanovic e colleghi interpreta il collasso della sella glaciale che univa i complessi Laurentide e Cordilleran come un passaggio dinamico cruciale, capace di convertire una fusione accelerata in una forzante climatica di scala emisferica; questa lettura si inserisce in un contesto paleoclimatico in cui il Bølling-Allerød rappresenta una delle più rapide transizioni calde dell’ultima deglaciazione, mentre il Meltwater Pulse 1A costituisce uno dei più rapidi episodi di innalzamento del livello del mare del tardo Quaternario. Gli autori sviluppano il problema in modo innovativo perché non assumono un semplice “impulso” idealizzato di acqua dolce, ma collegano la dinamica delle calotte a un modello di drenaggio ad alta risoluzione e poi a un modello climatico accoppiato, mostrando che la distribuzione geografica dell’acqua di fusione è decisiva per valutare la risposta dell’Atlantico settentrionale: proprio questo punto era già emerso in studi precedenti, nei quali la sensibilità dell’AMOC dipendeva fortemente dalla sorgente del forcing di acqua dolce, con scenari che ammettevano sia un contributo nordamericano dominante sia una partecipazione antartica o mista alla genesi del MWP-1A. In questa prospettiva, il collasso della sella Cordilleran-Laurentide non viene presentato come un semplice dettaglio geomorfologico, ma come un possibile meccanismo di soglia: il riscaldamento del Bølling avrebbe favorito l’abbassamento topografico e la destabilizzazione della sella, liberando grandi volumi di meltwater in tempi relativamente brevi; una volta immessa negli oceani, questa acqua dolce avrebbe ridotto la densità superficiale del Nord Atlantico, indebolito la formazione di acque profonde e quindi attenuato la circolazione meridionale di ribaltamento, producendo un raffreddamento diffuso nell’emisfero settentrionale coerente con l’idea che le variazioni dell’AMOC siano strettamente accoppiate alle oscillazioni climatiche deglaciali. In tal senso, il contributo principale dello studio è duplice: da un lato consolida l’ipotesi, già proposta da Gregoire et al. e rafforzata dai lavori di fingerprinting glaciale, che la sella nordamericana abbia fornito una parte sostanziale dell’acqua di fusione associata al MWP-1A; dall’altro mostra che tale apporto non sarebbe stato climaticamente neutro, ma avrebbe potuto contribuire alla fine del riscaldamento iniziale del Bølling o alla successiva fase fredda dell’Older Dryas, inserendo il collasso della sella in una catena causale che lega dinamica delle calotte, livello del mare, circolazione termoalina e brusche transizioni climatiche del tardo Pleistocene.
La figura 1 sintetizza in modo estremamente efficace l’idea centrale dello studio, cioè che il collasso della sella glaciale tra la Cordilleran Ice Sheet e la Laurentide Ice Sheet non rappresentò soltanto un episodio di rapida perdita di massa glaciale, ma un vero meccanismo di riorganizzazione idrologica e climatica durante la deglaciazione tardoglaciale: nel pannello (a) le tonalità rosse mostrano il forte assottigliamento del ghiaccio tra 14.5 e 14.3 ka proprio nella zona di contatto fra i due complessi glaciali, evidenziando la separazione delle due calotte e quindi il cosiddetto saddle collapse, mentre la sovrapposizione delle anomalie di salinità superficiale marina indica che questo rilascio di meltwater produce una marcata freshening del Nord Atlantico, con anomalie negative particolarmente estese nel bacino subpolare; si tratta di un risultato fisicamente cruciale, perché una diminuzione della salinità superficiale riduce la densità delle acque e tende a inibire la formazione di acque profonde, indebolendo l’AMOC. Questo schema è coerente sia con il lavoro originario di Ivanovic et al., che stima per il collasso della sella un apporto equivalente a 7.26 m di innalzamento medio globale del livello marino in circa 340 anni e un indebolimento dell’AMOC di circa 6 Sv, sia con gli studi precedenti di Gregoire et al., che avevano già identificato nel collasso della sella un meccanismo dinamico capace di produrre una rapida accelerazione della fusione attraverso un feedback quota–bilancio di massa, fino a generare diversi metri di sea-level equivalent in pochi secoli. I cerchi colorati e i bacini delimitati in nero mostrano inoltre che la questione non riguarda soltanto quanta acqua dolce venga rilasciata, ma soprattutto dove essa venga convogliata: i pannelli (b–e) confrontano infatti diversi scenari di drenaggio e mostrano che, quando il collasso della sella è incluso (SC_south e SC_east), compare un impulso molto più netto di scarico d’acqua dolce rispetto agli scenari senza collasso, confermando che la geografia del drenaggio è parte integrante del forcing climatico. Questo aspetto è in piena continuità con la letteratura paleoclimatica sul routing deglaciale: Wickert ha mostrato che la ricostruzione dei bacini e dei percorsi di drenaggio nordamericani varia sensibilmente nel tempo deglaciale, mentre Tarasov e Peltier avevano già evidenziato quanto una cronologia fisicamente coerente del drenaggio sia indispensabile per interpretare gli eventi abrupti; allo stesso modo, Condron e Winsor hanno dimostrato che acque dolci immesse attraverso il sistema artico–Mackenzie possono raggiungere molto più efficacemente le regioni di deep-water formation rispetto a un routing più meridionale via San Lorenzo, con un impatto molto maggiore sull’AMOC. In questo senso, la figura 1 non è una semplice carta descrittiva, ma la dimostrazione visiva che il collasso della sella Cordilleran-Laurentide costituì un evento di soglia capace di collegare in modo causale la dinamica delle calotte, il Meltwater Pulse 1A, la freshening del Nord Atlantico e il raffreddamento emisferico successivo, probabilmente contribuendo alla chiusura della fase di rapido riscaldamento del Bølling e alla transizione verso condizioni più fredde come l’Older Dryas.
2. Metodologia – 2.1 Storie di fusione delle calotte glaciali
In questa sottosezione, l’impianto metodologico dello studio si fonda su una scelta molto rilevante dal punto di vista paleoclimatico: invece di imporre al modello climatico un impulso di acqua dolce idealizzato, gli autori ricostruiscono la storia della fusione glaciale partendo da una simulazione dinamica delle calotte e la collegano a un modello di drenaggio ad alta risoluzione, così da trasferire all’oceano un forcing più fisicamente coerente nello spazio e nel tempo. Il nucleo di questa ricostruzione deriva da Gregoire et al. (2012), che usarono le uscite transitorie del modello climatico FAMOUS per gli ultimi 21 mila anni come forzante del modello glaciologico Glimmer-CISM, simulando l’evoluzione deglaciale delle calotte nordamericana e groenlandese e individuando nel collasso della sella Cordilleran-Laurentide un meccanismo capace di accelerare drasticamente la fusione; in quel quadro, l’abbassamento progressivo della sella tra due cupole glaciali favoriva un feedback altitudine–ablazione che poteva produrre rapidi impulsi di innalzamento del livello marino. Ivanovic et al. accettano questa ipotesi dinamica, ma compiono un passaggio metodologico decisivo: poiché nel modello originario il picco di fusione cadeva intorno a 11.6 ka, essi traslano l’evento di circa 2900 anni all’indietro per far coincidere il massimo della fusione nordamericana con il Meltwater Pulse 1A attorno a 14.5 ka, seguendo la possibilità che il collasso fosse stato innescato, o quantomeno fortemente favorito, dal rapido riscaldamento del Bølling; questa impostazione è coerente con Gregoire et al. (2016), che hanno mostrato come il brusco warming del Bølling potesse già di per sé generare 3–4 m di innalzamento medio globale del mare in 340 anni e predisporre il sistema glaciale nordamericano a una destabilizzazione della sella su scala plurisecolare. Il valore aggiunto del metodo sta inoltre nel fatto che la fusione non viene trattata come un’unica sorgente omogenea: i risultati glaciologici di Gregoire et al. sono combinati con il modello di drenaggio di Wickert, che ricostruisce i bacini nordamericani e i percorsi del deflusso in funzione di topografia, spessore del ghiaccio e aggiustamento isostatico, mentre per Eurasia e Antartide il forcing viene integrato attraverso la ricostruzione ICE-6G_C (VM5a), poi trasferito nel modello accoppiato HadCM3 come flusso globale e spazialmente distribuito di meltwater. In termini metodologici, questo passaggio è cruciale perché numerosi studi hanno mostrato che la risposta dell’AMOC e del clima superficiale dipende non solo dal volume totale di acqua dolce, ma anche dal suo punto di ingresso nei bacini oceanici; proprio per questo la ricostruzione dei percorsi di drenaggio è una variabile diagnostica, non un dettaglio secondario. La sottosezione evidenzia anche che gli autori non assumono un’origine esclusivamente nordamericana del MWP-1A in modo dogmatico, ma adottano una strategia intermedia: il Nord America resta la sorgente dominante del forcing, in accordo con la letteratura che attribuisce al collasso della sella circa 7 m di sea-level rise in circa 350 anni, mentre il fingerprinting del livello marino suggerisce che contributi addizionali da altre calotte siano verosimili; Gomez et al. hanno mostrato che uno scenario con collasso della sella più ~3–4 m equivalenti da altre masse glaciali soddisfa i vincoli dei record far-field, e Golledge et al. hanno indicato che l’Antartide potrebbe aver contribuito fino a circa 2 m in 340 anni durante MWP-1A. In coerenza con questa impostazione, Ivanovic et al. inseriscono nella forzante di fondo contributi di lungo periodo pari a 0.41 m dall’Antartide e 1.36 m dall’Eurasia in 340 anni, con portate massime inferiori rispettivamente a 0.02 e 0.05 Sv, ben al di sotto della scarica nordamericana, che nel loro schema resta sempre superiore a 0.1 Sv; ne emerge quindi una metodologia che non mira a chiudere definitivamente il dibattito sulle sorgenti del MWP-1A, ma a isolare in modo fisicamente plausibile l’effetto climatico specifico del collasso della sella nordamericana come primo passo per districare la sequenza di eventi avvenuti attorno a 14.5 ka.
2. Metodologia – 2.2 Calcolo del drenaggio dell’acqua di fusione
In questa sottosezione, il punto metodologico decisivo è che l’acqua di fusione non viene imposta al modello climatico come un forcing semplificato e puntuale, ma viene ricostruita come un flusso transitorio fisicamente distribuito, derivato dal bilancio di massa delle calotte e instradato verso l’oceano su una griglia di 30 arcsec (cioè inferiore a 1 km) in funzione della topografia superficiale, dello spessore del ghiaccio e dell’aggiustamento isostatico glaciale (GIA); per il Nord America e la Groenlandia gli autori usano le storie di fusione di Gregoire et al. con passo di 100 anni, mentre per Eurasia e Antartide incorporano la ricostruzione ICE-6G_C (VM5a) con passo di 500 anni, producendo così un forcing costiero globale a 100 anni che conserva la reale distribuzione geografica degli sbocchi di drenaggio. Il valore di questa scelta è elevato, perché la letteratura mostra da tempo che la risposta dell’Atlantico settentrionale dipende in modo critico non solo dal volume totale di acqua dolce, ma soprattutto dal routing, cioè dalla posizione in cui essa raggiunge il mare: Condron e Winsor, ad esempio, hanno mostrato che lo scarico attraverso i sistemi artici può perturbare la formazione di acque profonde molto più efficacemente rispetto a percorsi più meridionali, con conseguenze maggiori per l’AMOC e per il raffreddamento abrupto deglaciale. In questo quadro, Ivanovic e colleghi collegano direttamente il collasso della sella Cordilleran-Laurentide alla rete di drenaggio, integrando in un unico schema sia l’evento impulsivo sia la fusione di fondo, e trasformando quindi la dinamica glaciale in una forzante climatica geograficamente coerente; questa impostazione si basa sui risultati di Gregoire et al., secondo cui il collasso della sella poteva generare una forte accelerazione della perdita di massa glaciale, e sui successivi lavori che hanno proposto il Bølling warming come possibile innesco del collasso stesso. La sottosezione mette però in evidenza anche un nodo tecnico cruciale, cioè l’incertezza sul drenaggio dell’area dei Grandi Laghi: combinando la geometria glaciale di G12 con il modello della Terra solida VM2, Wickert ha ottenuto un GIA troppo forte, tale da creare un’insenatura marina artificiale nella regione laurenziana e quindi sbocchi continentali poco realistici per la fusione attorno e a ovest del Lago Superiore; per questo gli autori correggono manualmente il routing in modo da rispettare i vincoli geologici disponibili. La scelta di definire due scenari, SC_south e SC_east, nasce esattamente da questa incertezza: nel primo tutta l’acqua associata al settore del Lago Superiore viene inviata al Mississippi, nel secondo il 40% viene deviato verso il fiume Hudson, mentre gli altri bacini nordamericani risultano relativamente più stabili perché meno sensibili alla geometria del margine glaciale in quel momento. Questa impostazione è coerente con le ricostruzioni paleoidrografiche secondo cui, prima dell’apertura stabile del drenaggio verso il San Lorenzo intorno a 12.9 ka, gran parte del deflusso meridionale laurentidiano veniva convogliata verso il Mississippi e solo successivamente si osserva una riorganizzazione verso gli sbocchi orientali; tracce geochimiche e ricostruzioni di bacino sostengono infatti una modifica tardiva del routing al passaggio verso lo Younger Dryas. Nel complesso, questa sottosezione mostra bene che la metodologia dello studio non consiste semplicemente nel calcolare “quanta” acqua di fusione sia stata prodotta, ma nel determinare come, quando e soprattutto dove essa sia stata trasferita all’oceano, perché è proprio questa geografia del drenaggio a controllare la capacità del collasso della sella nordamericana di perturbare la salinità superficiale del Nord Atlantico, indebolire la circolazione di overturning e contribuire agli abrupt climate shifts di circa 14.5 ka; non a caso, anche lavori più recenti sulla tempistica della separazione Cordilleran-Laurentide confermano che il collasso della sella rappresentò un contributo sostanziale alla rapida immissione di meltwater durante il deglaciale, rafforzando la plausibilità dell’impostazione adottata qui.
2. Metodologia – 2.3 Il modello climatico
In questa sottosezione, la robustezza dell’impianto sperimentale dipende in larga misura dalla scelta di HadCM3, un modello climatico globale accoppiato atmosfera–oceano senza flux adjustment, a cui in questa configurazione viene associato anche un modulo di vegetazione interattiva; si tratta di una famiglia modellistica ampiamente impiegata negli studi paleoclimatici proprio perché combina una dinamica completa dei principali comparti del sistema climatico con un costo computazionale abbastanza contenuto da consentire esperimenti di lunga durata e numerose simulazioni di sensibilità. Nello specifico, il componente atmosferico opera alla risoluzione standard di circa 3.75° in longitudine × 2.5° in latitudine con 19 livelli verticali ibridi, mentre il componente oceanico ha una griglia più fine, pari a circa 1.25° × 1.25°, con 20 livelli verticali; questa asimmetria di risoluzione riflette una scelta tipica dei GCM della generazione HadCM3, nei quali la struttura oceanica deve descrivere in modo sufficientemente realistico il trasporto di calore, salinità e massa d’acqua, pur restando dentro un assetto numericamente stabile per integrazioni climatiche di lungo periodo. Accoppiato al sistema atmosfera–suolo è presente inoltre TRIFFID (Top-down Representation of Interactive Foliage and Flora Including Dynamics), cioè un modello dinamico della vegetazione che aggiorna copertura e struttura dei principali tipi funzionali vegetali in risposta ai bilanci di carbonio e alla competizione ecologica; la sua inclusione è metodologicamente rilevante perché consente di non trattare la superficie continentale come una semplice condizione al contorno fissa, ma come un elemento capace di retroagire sui flussi di energia, acqua e carbonio. Un aspetto cruciale, richiamato dagli autori, è però il fatto che l’oceano di HadCM3 utilizza un fixed lid, cioè un “coperchio rigido” che non permette la variazione del volume delle celle oceaniche e quindi del livello marino locale: per questo motivo i flussi idrologici, inclusi evaporazione, precipitazione, runoff fluviale e fusione glaciale, non vengono introdotti come veri apporti di massa ma come flussi virtuali di salinità, una scelta molto importante quando si simula l’immissione di acqua dolce deglaciale. Ivanovic e colleghi riconoscono che una rappresentazione con flussi reali potrebbe, in linea teorica, modificare leggermente i percorsi di propagazione del freshwater attraverso effetti sterici e di livello del mare, ma richiamano il fatto che per impulsi dell’ordine di ~0.1 Sv o inferiori tali differenze dovrebbero restare modeste; questa valutazione è coerente con l’analisi di Yin et al., che ha mostrato come il metodo del virtual salt flux possa introdurre alcune distorsioni, ma in molti casi senza alterare in modo sostanziale la risposta di larga scala della circolazione rispetto ai forcing di acqua dolce di entità moderata. Dal punto di vista della logica sperimentale, ciò significa che HadCM3 non viene usato qui come semplice “contenitore” numerico, ma come strumento capace di tradurre un forcing di meltwater geograficamente distribuito in una risposta coerente della salinità superficiale, della formazione di acque profonde e della circolazione atlantica; inoltre, il fatto che il modello disponga di bacini idrografici terrestri definiti sulla griglia, i quali trasferiscono istantaneamente l’acqua meteorica alla costa, mentre il routing dell’acqua di fusione glaciale viene fornito separatamente come input esterno variabile, permette agli autori di isolare in modo più pulito l’effetto specifico del collasso della sella Cordilleran-Laurentide rispetto al runoff continentale “ordinario”. In termini metodologici, è proprio questa combinazione fra GCM accoppiato, vegetazione interattiva, schema idrologico terrestre e forcing glaciale esterno distribuito a rendere la sottosezione particolarmente importante, perché mostra come lo studio cerchi di rappresentare il più realisticamente possibile non solo la quantità di acqua dolce immessa nell’oceano, ma anche la sua capacità di perturbare il sistema climatico emisferico attraverso i meccanismi di salinità e overturning oceanico.
2. Metodologia – 2.4 Progettazione degli esperimenti con il modello climatico
La progettazione sperimentale adottata in questa sottosezione è particolarmente solida perché mira a isolare nel modo più pulito possibile l’effetto climatico del collasso della sella glaciale Cordilleran-Laurentide dal rumore introdotto da altre forzanti deglaciali: gli autori partono infatti da una simulazione di controllo a 15 ka basata sul setup di Singarayer et al., con parametri orbitali derivati da Berger e Loutre, concentrazioni atmosferiche fissate a 225 ppmv per la CO₂, 473 ppbv per il CH₄ e 241 ppbv per il N₂O, e condizioni al contorno glacio-topografiche, batimetriche e costiere tratte dalla ricostruzione ICE-5G (VM2); tutte queste boundary conditions vengono mantenute costanti durante le integrazioni proprio per evitare che l’evoluzione indipendente di orbita, gas serra o geometria glaciale interferisca con l’attribuzione causale della risposta climatica al solo forcing di acqua di fusione. In termini metodologici, questa è una scelta classica ma molto efficace negli studi di sensibilità del sistema climatico, perché consente di trasformare l’esperimento in una sorta di “test controllato” sull’impatto del freshwater forcing, riducendo l’ambiguità interpretativa che spesso affligge le simulazioni deglaciali in cui più forzanti evolvono simultaneamente; inoltre, il ricorso a HadCM3, una famiglia modellistica nota per la sua stabilità nelle integrazioni paleoclimatiche di lunga durata e per la possibilità di lavorare senza flux correction, rende questo disegno particolarmente adatto a esplorare perturbazioni secolari e millenarie dell’AMOC. Il controllo viene fatto evolvere per 1000 anni con vegetazione prescritta e poi per altri 750 anni con TRIFFID attivo, per un totale di 1750 anni di spin-up, fino a raggiungere uno stato quasi stazionario di atmosfera, oceano e superficie continentale; la prosecuzione parallela della simulazione di controllo per altri 2000 anni durante gli esperimenti perturbati serve inoltre a verificare che il modello resti effettivamente vicino all’equilibrio per tutta la durata dell’analisi, rafforzando la credibilità diagnostica delle anomalie simulate. Su questa base vengono costruiti quattro esperimenti di 2000 anni nell’intervallo 15.5–13.5 ka: due simulazioni transitorie, SC_south e SC_east, includono l’impulso di meltwater legato al collasso della sella, mentre due simulazioni di riferimento, NoSC_south e NoSC_east, mantengono il flusso di fusione bloccato ai valori di 14.8 ka, così da rimuovere il pulse secolare e conservare soltanto la fusione di fondo su scala multimillenaria; il confronto fra queste coppie sperimentali è metodologicamente elegante perché permette di separare l’effetto del collasso rapido della sella dal background melt, cioè di distinguere una perturbazione impulsiva da una forzante lenta ma persistente. Un altro elemento importante è il modo in cui il forcing viene applicato al modello: i flussi di acqua di fusione, variabili nel tempo, sono interpolati sulla griglia climatica e distribuiti nelle celle oceaniche adiacenti alla costa con profondità superiore a 500 m, lungo il margine costiero e oltre la piattaforma continentale; questa scelta numerica evita il raggiungimento di salinità negative, che sarebbero fisicamente impossibili, e mantiene le proprietà del sistema entro il dominio valido dell’equazione di stato dell’acqua marina. Anche il fatto che i flussi non vengano concentrati in una sola cella, ma “spalmati” su più celle oceaniche, è metodologicamente cruciale, perché impedisce la generazione di artefatti numerici e rende più realistica la propagazione del segnale di freshening; si tratta di una cautela particolarmente rilevante nei modelli con fixed lid e flussi virtuali di salinità, dove la corretta rappresentazione del forcing di acqua dolce dipende fortemente dal modo in cui esso viene discretizzato nello spazio. Infine, gli autori scelgono di non simulare i flussi iperpicnali di acqua di fusione, giustificando questa omissione con la loro rarità e con il contributo generalmente trascurabile, su scala oceanica globale, dei piccoli bacini ripidi in cui essi tendono a svilupparsi; ciò conferma che l’esperimento non intende riprodurre ogni dettaglio sedimentologico locale, ma catturare con coerenza fisica il segnale climatico di larga scala associato al routing del meltwater. Nel complesso, questa sottosezione mostra una progettazione sperimentale molto ben costruita, in cui controllo equilibrato, boundary conditions fisse, coppie di scenari con e senza pulse e applicazione spazialmente distribuita del forcing convergono in un disegno capace di attribuire con buona robustezza il rallentamento dell’AMOC e le anomalie climatiche simulate al collasso della sella glaciale nordamericana.
3. Risultati – 3.1 Salinità superficiale del mare e circolazione oceanica
In questa sottosezione emerge con chiarezza che l’elemento decisivo non è soltanto la quantità totale di acqua di fusione prodotta dal collasso della sella Cordilleran-Laurentide, ma soprattutto il suo instradamento verso l’Artico e l’Atlantico: nello scenario SC_south circa il 50% del meltwater viene convogliato dal Mackenzie verso l’Oceano Artico e poco meno del 50% attraverso il Mississippi verso il Golfo del Messico e l’Atlantico, mentre nello scenario SC_east la quota diretta al Mississippi scende a circa il 30% e circa il 20% del flusso viene deviato verso la costa orientale nordamericana, soprattutto tramite San Lorenzo e Hudson; tuttavia, su scala secolare, gli autori mostrano che questa redistribuzione non produce differenze apprezzabili nella salinità superficiale né nella circolazione oceanica su larga scala, perché l’acqua dolce che entra nell’Atlantico occidentale viene rapidamente dispersa dai giri subtropicali e, in misura minore, subpolari. Questo risultato è metodologicamente importante, perché suggerisce una sensibilità relativamente bassa alla posizione esatta di immissione del meltwater nella regione del giro subtropicale atlantico, e giustifica il fatto che la discussione successiva si concentri soprattutto sullo scenario SC_south come riferimento principale.
Il cuore fisico del risultato è la formazione di un vero e proprio freshwater cap artico-nordatlantico: a circa 14.4 ka, cioè un secolo dopo il picco del collasso, l’Artico e il settore GIN raggiungono il minimo di salinità superficiale, con acque nei primi 200 m rese più fresche di circa 1–6 psu rispetto allo scenario senza impulso di collasso; questa anomalia si propaga poi attraverso lo Stretto di Fram verso il Nord Atlantico, dove si combina con circa 0.07 Sv di ulteriore fusione proveniente da Hudson Strait, costa orientale nordamericana e Mississippi. L’effetto dinamico è un aumento del gradiente verticale di densità, quindi una maggiore stratificazione della colonna d’acqua alle medie e alte latitudini, con conseguente riduzione della formazione di North Atlantic Deep Water: nel modello, l’AMOC diventa a 14.4 ka circa 6 Sv più debole — ossia quasi il 40% in meno — e 200–400 m più superficiale, per poi recuperare dopo 14.1 ka fino a mostrare un possibile overshoot, in linea con la dinamica di recupero già discussa da Liu et al. (2009)per il Bølling-Allerød. L’interpretazione proposta da Ivanovic et al. è coerente anche con Condron e Winsor (2012), che hanno mostrato come l’acqua dolce scaricata attraverso il sistema artico–Mackenzie raggiunga molto più efficacemente le regioni di formazione di acque profonde rispetto al drenaggio via San Lorenzo, indebolendo quindi l’AMOC in modo più netto; allo stesso tempo, il fatto che in questo studio la risposta del Pacifico settentrionale sia trascurabile non contraddice la letteratura, ma indica che il forcing qui simulato è probabilmente troppo modesto per attivare una vera riorganizzazione del Pacific overturning come quella discussa da Okazaki et al. (2010) e, in forma più debole o dipendente dal modello, da Chikamoto et al. (2012). Nel complesso, questa sottosezione documenta quindi una catena causale molto coerente: routing artico del meltwater → freshening superficiale e maggiore stratificazione → riduzione della convezione profonda nordatlantica → indebolimento multicentenario dell’AMOC, cioè esattamente il tipo di risposta oceanografica che rende plausibile un collegamento tra collasso della sella glaciale, fine del riscaldamento del Bølling e successivo raffreddamento dell’Older Dryas.
3. Risultati – 3.2 Clima superficiale
La sottosezione dedicata al clima superficiale mostra che l’indebolimento dell’AMOC simulato intorno a 14.4 ka non produce soltanto una risposta oceanica, ma si traduce in una riorganizzazione atmosferica emisferica ben strutturata: la riduzione del trasporto verso nord di calore oceanico superficiale e intermedio determina un raffreddamento diffuso dell’emisfero settentrionale e un lieve riscaldamento dell’emisfero meridionale, configurando un classico pattern bipolare che nel lavoro di Ivanovic et al. raggiunge in inverno boreale anomalie termiche generalmente comprese tra 1 e 5 °C, con punte ancora più marcate nei mari GIN e nel Labrador, proprio nelle aree maggiormente influenzate dalla freshening e dalla riduzione della convezione profonda. In questo assetto, il raffreddamento alle alte latitudini favorisce anche un aumento del ghiaccio marino nel Nord Atlantico e nei mari nordici, soprattutto nella stagione calda, e tale espansione agisce come feedback positivo tramite incremento dell’albedo e ulteriore isolamento termico tra oceano e atmosfera; ne consegue una risposta climatica coerente con molti esperimenti di freshwater hosing glaciali, nei quali l’indebolimento della circolazione atlantica amplifica il raffreddamento subpolare, modifica i gradienti meridionali di temperatura e riduce l’evaporazione e le precipitazioni alle alte latitudini dell’emisfero nord. Un aspetto particolarmente importante, e ben inserito nella letteratura dinamica sul clima tropicale, è lo spostamento verso sud dell’ITCZ: Ivanovic et al. lo interpretano come risposta alla nuova asimmetria energetica interemisferica, e questa lettura è pienamente compatibile con gli studi di Chiang e Bitz, che hanno mostrato come un raffreddamento imposto alle alte latitudini dell’emisfero nord tenda a far migrare la convergenza tropicale verso l’emisfero opposto, nel tentativo di riequilibrare il trasporto energetico atmosferico. Gli autori sottolineano inoltre che il pattern simulato è molto simile a quello ottenuto nei confronti multimodello di Kageyama et al., dove gli esperimenti LGM con hosing nordatlantico producono un segnale caratterizzato da raffreddamento marcato nel Nord Atlantico, riorganizzazione tropicale e teleconnessioni atmosferiche di ampia scala; ciò rafforza l’interpretazione secondo cui il collasso della sella Cordilleran-Laurentide possa aver avuto effetti climatici non solo regionali, ma emisferici, attraverso il nesso causale meltwater → AMOC debole → riorganizzazione termica e idrologica. Un risultato molto interessante è poi il riscaldamento di 1–3 °C sull’America nord-orientale, che contrasta con il raffreddamento generale dell’emisfero settentrionale e deriva, secondo il modello, da una riorganizzazione della circolazione atmosferica associata al forte raffreddamento dell’Atlantico: si instaura così un dipolo fra Nord Atlantico ad alte latitudini più freddo e settore continentale orientale nordamericano relativamente più mite, un segnale che compare anche in altri GCM quando il raffreddamento subpolare supera una certa soglia e si estende verso i tropici atlantici. Infine, la sottosezione mette in evidenza possibili effetti tropicali e continentali secondari, come il riscaldamento e l’essiccamento tra Panama e Brasile e una risposta della vegetazione amazzonica potenzialmente favorevole a un arretramento della foresta pluviale, un comportamento che trova analogie negli esperimenti di Bozbiyik et al. su AMOC indebolita e ciclo del carbonio, anche se gli stessi autori riconoscono che questa componente è meno robusta rispetto ai risultati principali su Atlantico nordico, ghiaccio marino e spostamento dell’ITCZ.
Collapse of the North American ice saddle 14,500 years ago caused widespread cooling and reduced ocean overturning circulation – Figura 2
La figura 2 rappresenta il nucleo dinamico dell’intero studio, perché mostra in forma sintetica la sequenza fisica che collega il collasso della sella glaciale Cordilleran–Laurentide alla risposta climatica emisferica: nel pannello (a) la salinità superficiale media dell’Artico e dei mari Groenlandia–Islanda–Norvegia cala rapidamente dopo l’avvio dell’evento, nel pannello (b) l’AMOC raggiunge il proprio minimo quasi in fase con questo intenso freshening, nel pannello (c) la Groenlandia centrale registra un marcato raffreddamento, mentre i pannelli (d) ed (e) mostrano che la risposta atmosferica è dominata da un’anomalia fredda sul settore nord-atlantico e artico, molto più intensa in inverno che in estate; nel lavoro di Ivanovic et al. questo insieme di segnali è interpretato come la risposta a un impulso di meltwater capace di indebolire l’AMOC di circa 6 Sv e di generare un raffreddamento diffuso dell’emisfero settentrionale dell’ordine di 1–5 °C. Il punto scientificamente più importante è che la figura non presenta anomalie indipendenti, ma una vera catena causale: l’acqua dolce riduce la salinità superficiale alle alte latitudini, aumenta la stratificazione della colonna d’acqua, inibisce la formazione di North Atlantic Deep Water e quindi riduce il trasporto meridionale di calore, con conseguente raffreddamento dell’Atlantico subpolare e delle terre circostanti; studi successivi degli stessi autori, riferiti ad altri intervalli deglaciali, hanno confermato lo stesso schema fisico, mostrando che la diminuzione della salinità alle alte latitudini può aumentare la stratificazione del 20–80% nei primi 200 m e indebolire in modo significativo la circolazione di overturning, segno che il comportamento illustrato nella figura 2 è dinamicamente robusto e non un semplice artefatto di configurazione modellistica. Dal punto di vista climatico, la maggiore intensità del segnale invernale rispetto a quello estivo è del tutto plausibile, perché l’espansione del ghiaccio marino rafforza l’albedo e soprattutto isola termicamente l’atmosfera dall’oceano, amplificando il raffreddamento sopra i mari GIN, il Labrador e la Groenlandia; nello stesso tempo, la figura evidenzia anche un’area relativamente più mite sull’America nord-orientale, cioè un dipolo già riconosciuto in diversi esperimenti di freshwater hosing, nei quali il forte raffreddamento del Nord Atlantico si accompagna a una riorganizzazione della circolazione atmosferica e a teleconnessioni continentali coerenti con quanto simulato qui. Su scala più ampia, il quadro termico mostrato dalla figura 2 è perfettamente compatibile con la letteratura classica sulla risposta interemisferica ai raffreddamenti dell’emisfero nord: quando il Nord Atlantico perde calore e l’AMOC si indebolisce, l’asimmetria energetica tende a spingere l’ITCZ verso sud, mentre il successivo recupero dell’AMOC, visibile nel pannello (b) e forse associato a un breve overshoot, richiama i meccanismi proposti da Liu et al. per la deglaciazione, in cui la diminuzione o cessazione del forcing di acqua dolce consente una rapida riattivazione della convezione nord-atlantica e un temporaneo superamento dei valori precedenti della circolazione.
4. Discussione e conclusioni
Nella parte conclusiva dello studio, Ivanovic e colleghi sostengono che il vero avanzamento non risiede soltanto nel risultato climatico finale, ma soprattutto nell’aver costruito una catena modellistica fisicamente coerente in cui una simulazione climatica transitoria forza un modello dinamico di calotta, il quale alimenta a sua volta un modello di drenaggio ad alta risoluzione; il forcing di acqua di fusione così ottenuto viene poi reintrodotto nel GCM, superando gli schemi troppo idealizzati di freshwater hosing che hanno caratterizzato molta modellistica precedente. In questo quadro, il collasso della sella Cordilleran–Laurentide emerge come un meccanismo credibile di instabilità deglaciale: già Gregoire et al. avevano mostrato che le “ice saddles” possono collassare rapidamente quando il riscaldamento abbassa la quota della sella e amplifica l’ablazione, producendo impulsi rapidi di innalzamento del livello marino, mentre Ivanovic et al. quantificano per l’evento intorno a 14.5 ka un indebolimento dell’AMOC di circa 6 Sv e un raffreddamento diffuso dell’emisfero settentrionale, con anomalie invernali fino a circa 6 °C sulla Groenlandia centrale. Un risultato particolarmente importante della discussione è che circa metà del meltwater del collasso viene instradato verso l’Oceano Artico e che la risposta climatica è molto più sensibile a questo routing artico che non alle differenze fra scarico sul Mississippi o sulla East Coast nordamericana; questo punto è in linea con Condron e Winsor, che hanno mostrato come il drenaggio artico possa raggiungere molto più efficacemente le regioni di formazione di acque profonde del Nord Atlantico e indebolire l’AMOC più del routing via San Lorenzo o subtropical gyre. La discussione insiste però anche sulle incertezze cronologiche: il collasso della sella è plausibile come risposta al rapido riscaldamento del Bølling, e Gregoire et al. hanno mostrato che il solo warming del primo Bølling può produrre 3–4 m di innalzamento medio globale del mare in 340 anni, amplificabili a circa 5–6 m se il riscaldamento innesca anche il collasso della sella; ciò rafforza l’idea che il meccanismo nordamericano possa aver contribuito in modo sostanziale al Meltwater Pulse 1A, pur senza esaurirne necessariamente l’intera entità. Proprio qui si inserisce il confronto con Gomez et al., che hanno mostrato come uno scenario con collasso della sella più circa 3 m di fusione esterna a Nord America e Groenlandia sia compatibile con i principali record far-field del MWP-1A, mentre Golledge et al. hanno fornito un supporto indipendente all’idea che anche l’Antartide possa aver contribuito all’episodio attraverso una fase di recessione accelerata legata alla riduzione dell’overturning nell’Oceano Meridionale. Ne deriva una conclusione equilibrata: il collasso della sella Cordilleran–Laurentide non va interpretato come unica causa di tutta la deglaciazione rapida di 14.5 ka, ma come un forcing sub-millenario molto robusto, capace di collegare riscaldamento del Bølling, rilascio accelerato di meltwater, indebolimento multicentenario dell’AMOC e successivo raffreddamento, plausibilmente associabile alla fine del warming iniziale del Bølling e/o all’Older Dryas.
Collapse of the North American ice saddle 14,500 years ago caused widespread cooling and reduced ocean overturning circulation – Figura 3
La figura 3 svolge una funzione interpretativa cruciale, perché non si limita a mostrare un risultato del modello, ma mette in relazione diretta il segnale simulato del collasso della sella glaciale Laurentide–Cordilleran con la variabilità termica ricostruita per la Groenlandia centrale durante la deglaciazione, utilizzando il record di temperatura di Buizert et al. come riferimento osservativo per valutare se il raffreddamento modellato abbia una controparte plausibile nei dati paleoclimatici. In questo senso, il pannello (a) va letto come una ricerca di “finestre” nel record groenlandese entro cui un impulso di meltwater nordamericano avrebbe potuto lasciare una firma climatica coerente con un temporaneo indebolimento dell’AMOC, mentre il pannello (b) traduce questa lettura in termini quantitativi, confrontando durata e massimo raffreddamento degli episodi identificati nel record con quelli dello scenario modellistico SC_south. Il quadro generale è coerente con la ricostruzione di Buizert et al., che mostra come la Groenlandia abbia attraversato nella deglaciazione una sequenza di cambiamenti termici estremamente rapidi, compreso il brusco riscaldamento del Bølling intorno a 14.7 ka, e con l’interpretazione di Ivanovic et al., secondo cui il collasso della sella nordamericana potrebbe aver introdotto, subito dopo quel warming, una perturbazione di segno opposto durata alcuni secoli.
Nel pannello (a), i riquadri colorati non rappresentano eventi “dimostrati” in senso stretto, ma episodi candidati di raffreddamento nel record della Groenlandia centrale che possono essere confrontati con il forcing simulato; la logica degli autori è che, se il collasso della sella ha davvero liberato grandi volumi di acqua dolce nell’Artico e nel Nord Atlantico, allora il risultato atteso è una riduzione multicentenaria del trasporto di calore oceanico verso nord, con un raffreddamento groenlandese di entità e durata paragonabili a quelle visibili nei box evidenziati. Il pannello (b) rende esplicito questo confronto: l’ellisse nera associata a SC_south si colloca nell’ordine di circa 600 anni di durata e circa 3 °C di massimo raffreddamento, cioè in un settore del diagramma che si sovrappone soprattutto agli eventi più lunghi e intensi del record, in particolare a quello attribuibile all’Older Dryas e a un episodio attorno a ~14.46 ka; al contrario, alcuni eventi più antichi, vicini a ~16.45 ka, risultano più brevi o meno intensi, ma restano comunque compatibili con uno scenario alternativo di collasso anticipato. Questa impostazione è perfettamente coerente con la conclusione del paper, secondo cui Scenario 1, collocato fra 14.5 e 14.0 ka, potrebbe spiegare il raffreddamento successivo al Bølling e/o l’Older Dryas, mentre Scenario 2, centrato attorno a ~16 ka, potrebbe corrispondere a un episodio freddo groenlandese più antico.
Dal punto di vista fisico, la figura 3 è importante perché suggerisce che il collasso della sella non vada interpretato solo come un contributo eustatico al Meltwater Pulse 1A, ma anche come un possibile modulatore della cronologia climatica deglaciale. Gregoire et al. hanno mostrato che il solo riscaldamento abrupto del Bølling poteva produrre 3–4 m di innalzamento medio globale del mare in 340 anni, ma che questo contributo poteva essere amplificato fino a 5–6 m se il warming innescava anche il collasso della sella glaciale; parallelamente, Gomez et al. hanno indicato che il collasso Laurentide–Cordilleran costituisce una sorgente plausibile e importante del MWP-1A, pur non esaurendone da solo tutta l’ampiezza. La figura 3 aggiunge a questa narrativa un tassello decisivo: se il collasso ha davvero contribuito in modo sostanziale al pulse di meltwater, allora deve aver lasciato non solo un segnale nel livello del mare, ma anche una firma termica transitoria nella Groenlandia e nel Nord Atlantico, ed è proprio questo che il confronto fra i box del pannello (a) e il diagramma durata–intensità del pannello (b) prova a dimostrare.
Il valore scientifico della figura, quindi, sta nel fatto che essa collega tre scale di evidenza che spesso restano separate: la dinamica delle calotte glaciali, la risposta della circolazione oceanica e il segnale paleotermico groenlandese. La conclusione implicita è che il collasso della sella Laurentide–Cordilleran costituisce un meccanismo plausibile per spiegare un raffreddamento multicentenario dopo il Bølling, soprattutto se si considera che Ivanovic et al. simulano per questo forcing un indebolimento dell’AMOC di circa 6 Sv, un raffreddamento diffuso dell’emisfero settentrionale di 1–5 °C e una risposta groenlandese abbastanza intensa da rientrare nel range degli episodi evidenziati nella figura 3. In altri termini, la figura non prova in modo definitivo quale box colorato corrisponda all’evento reale, ma mostra che la firma climatica attesa da un collasso della sella è compatibile con più di un episodio del record groenlandese, e che l’ipotesi più forte resta quella di un evento legato al post-Bølling e all’Older Dryas, senza escludere del tutto una collocazione più antica attorno a 16 ka.
https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1002/2016GL071849
0 commenti