Sea-level contribution from Thwaites Glacier doubles in an observation-constrained ice-sheet model

Sintesi Tecnica: Il contributo del ghiacciaio Thwaites al livello del mare raddoppia in un modello di calotta glaciale vincolato dalle osservazioni

Problematica
Il ghiacciaio Thwaites in l'Antartide Occidentale è un principale contributore alla perdita di massa antartica e una fonte critica di incertezza nelle proiezioni dell'innalzamento del livello del mare. Il suo rapido ritiro è guidato dalla fusione indotta dall'oceano alla linea di galleggiamento (grounding line), dove il ghiaccio ancorato transita in banchisa galleggiante. Sebbene studi precedenti abbiano evidenziato il ruolo delle acque profonde circolari (CDW) calde nell'erosione della base del ghiaccio, gli attuali modelli di calotta glaciale faticano a riprodurre simultaneamente i trend storici di ritiro, assottigliamento e accelerazione del ghiacciaio. Una carenza fondamentale identificata è il trattamento della linea di galleggiamento come un confine netto e discreto tra ghiaccio ancorato e galleggiante. Questo approccio trascura le intrusioni di acqua marina su scala chilometrica sotto il ghiaccio ancorato, recentemente osservate tramite satellite, che creano una "zona di galleggiamento" dinamica caratterizzata da fusione attiva. Inoltre, le parametrizzazioni standard della fusione spesso impongono tassi di fusione nulli alla linea di galleggiamento, contraddicendo le osservazioni che suggeriscono che la fusione sia sostanziale e potenzialmente dominante all'interno di questa zona.

Metodologia
Gli autori hanno utilizzato l'Ice-sheet and Sea-level System Model (ISSM) per simulare l'evoluzione del ghiacciaio Thwaites dal 1996 al 2023 e proiettarne il comportamento fino al 2100. Lo studio ha introdotto una nuova parametrizzazione per rappresentare esplicitamente una "zona di galleggiamento con fusione attiva" che si estende verso l'interno dalla linea di galleggiamento. In questo quadro, la zona di galleggiamento è definita come una regione di lunghezza $L$ in cui i tassi di fusione basale diminuiscono linearmente da un valore massimo alla linea di galleggiamento fino a zero alla distanza di intrusione.

Il modello è stato inizializzato per rappresentare lo stato di metà anni '90, incorporando una correzione localizzata di +100 m alla topografia del letto nelle regioni di ghiaccio galleggiante per garantire la coerenza con la posizione osservata della linea di galleggiamento del 1996. Lo studio ha impiegato un approccio di calibrazione bayesiana per valutare il modello rispetto a tre record osservativi indipendenti: velocità superficiale, spessore del ghiaccio e posizione della linea di galleggiamento. La calibrazione ha variato due parametri chiave: il tasso di fusione subacquea ($\dot{m}$) e la lunghezza della zona di galleggiamento ($L$). È stato eseguito un ensemble di hindcast di 105 simulazioni per identificare lo spazio dei parametri congiunto che minimizza il disaccordo con i trend storici. Dopo la calibrazione, le configurazioni con le migliori prestazioni sono state estese al 2100 sotto uno scenario di forcing climatico costante (mediato tra il 1996 e il 2023) per isolare l'impatto della fisica della zona di galleggiamento dalle incertezze del futuro riscaldamento oceanico.

Risultati Chiave

1. Ricostruzione dei Trend Storici: Il modello è riuscito a riprodurre il ritiro, l'assottigliamento e l'accelerazione osservati del ghiacciaio Thwaites nell'arco del periodo di 25 anni solo quando la zona di galleggiamento era rappresentata come una regione larga 2–5 km con tassi di fusione di circa 110 m l'anno. Le simulazioni che trattavano la linea di galleggiamento come un confine discreto richiedevano tassi di fusione irrealisticamente elevati (150 m l'anno) e non riuscivano comunque a catturare l'entità dell'assottigliamento e dell'accelerazione osservati a metà degli anni 2000.
2. Compensazioni tra i Parametri: L'analisi bayesiana ha rivelato una relazione compensativa tra la lunghezza della zona di galleggiamento e il tasso di fusione. Tuttavia, le simulazioni senza una zona di galleggiamento sono risultate quattro volte meno probabili nel riprodurre i trend storici rispetto a quelle con una zona di galleggiamento di diversi chilometri.
3. Proiezioni Future e Distacco: Sotto un forcing climatico costante, la lunghezza della zona di galleggiamento è emersa come il predittore dominante del fatto che il ghiacciaio si stacchi dalla stabilizzante cresta di Mouginot.
   • I modelli senza zona di galleggiamento o con una zona stretta di 1 km non si sono ritirati oltre la cresta entro il 2100.
   • I modelli con una zona di galleggiamento di 2–3 km hanno predetto una probabilità del 45% di distacco entro il 2100.
   • I modelli coerenti con le recenti osservazioni satellitari (zona di galleggiamento di 4–6 km) hanno mostrato una probabilità del ~60% di attraversare la cresta entro il 2100, con un potenziale distacco già nel 2038.
4. Contributo al Livello del Mare: Lo studio proietta che una lunghezza della zona di galleggiamento di 5±1 km (coerente con le recenti osservazioni ad alta risoluzione) risulta in un contributo al livello medio globale del mare di 30,1 ± 5,7 mm entro il 2100. Questo è circa il doppio della proiezione media dell'Ice Sheet Model Intercomparison Project per CMIP6 (ISMIP6) sotto lo stesso scenario di forcing costante (15,2 ± 21,8 mm) e rappresenta un incremento del 130% rispetto alla media ISMIP6.

Significato e Rivendicazioni
L'articolo sostiene che la risoluzione esplicita delle zone di galleggiamento su scala chilometrica e con fusione attiva altera fondamentalmente sia le capacità di hindcast che di previsione per il ghiacciaio Thwaites. Calibrando il modello rispetto a molteplici vincoli osservativi indipendenti, gli autori dimostrano che i processi della zona di galleggiamento governano il distacco dalle creste topografiche stabilizzanti, un meccanismo precedentemente sottorappresentato nei modelli.

Gli autori affermano che le loro scoperte stabiliscono un forte punto di riferimento per la capacità di proiezione, mostrando che l'estensione della zona di galleggiamento esercita un controllo più forte sul ritiro e sul distacco rispetto a variazioni moderate del forcing termico. Di conseguenza, lo studio sostiene che le zone di galleggiamento con fusione attiva devono essere sistematicamente incorporate nei modelli di calotta glaciale regionali e continentali, inclusi i quadri di confronto internazionali, per restringere l'incertezza sul contributo dell'Antartide al futuro innalzamento del livello del mare. Il documento sottolinea che, sebbene le simulazioni ghiaccio-oceano completamente accoppiate siano l'ideale, l'attuale approccio di calibrare i tassi di fusione efficaci rispetto ai cambiamenti di stato osservati fornisce una via affidabile per ridurre l'incertezza delle proiezioni.