Il quadro che emerge dal lavoro di Välisuo e colleghi è molto netto: quando si passa dalla “semplice” presenza di fusione in un giorno qualunque alla misura dell’impatto stagionale (quanti giorni fonde, quanto a lungo e su quanta superficie si accumula la fusione), l’atmosfera che conta di più non è solo “più calda”, ma soprattutto “più umida” e organizzata in configurazioni sinottiche capaci di convogliare calore e vapore verso la calotta. La distinzione tra estensione massima della fusione e indicatori integrati stagionalmente (numero di giorni di fusione, estensione cumulativa, quantità di fusione modellata) è fisicamente cruciale: l’estensione massima può essere dominata da eventi intensi ma brevi, su scale di pochi giorni, mentre le metriche cumulative “sommano” l’effetto persistente di circolazioni favorevoli lungo l’intera stagione. Questo è coerente con ciò che molti studi hanno evidenziato: gli estremi di fusione su vasta scala in Groenlandia sono spesso legati a episodi sinottici molto organizzati (creste/ridging, blocchi, trasporti meridionali di umidità), ma l’entità complessiva stagionale dipende dalla frequenza e persistenza di tali episodi.
Nel periodo 2000–2014, il fatto che la total column water (contenuto totale d’acqua nella colonna atmosferica) risulti fortemente correlata con NMD, estensione cumulativa e fusione modellata va letto come una firma dinamico-termodinamica: più vapore in colonna, in Groenlandia, raramente è “locale”, ma tipicamente segnala trasporto avvettivo di masse d’aria più miti e umide dall’Atlantico verso la calotta, spesso in presenza di un pattern di alta pressione/blocco che guida e sostiene l’afflusso. In termini energetici, l’umidità agisce su più canali: aumenta l’effetto serra locale (incremento della radiazione infrarossa discendente), favorisce nubi basse liquide con forte impatto radiativo in estate artica, e può intensificare i flussi turbolenti (sensibile/latente) quando l’aria mite scorre su superfici fredde o su aree prossime a 0 °C. Studi “iconici” sull’evento del luglio 2012 mostrano proprio questo meccanismo: nubi basse liquide e un’atmosfera insolitamente umida hanno contribuito in modo determinante all’anomalia di fusione, amplificando la radiazione infrarossa verso la superficie.
Il legame (più debole ma significativo) con indici di grande scala come NAO e soprattutto con il Greenland Blocking Index (GBI) è l’altra faccia della medaglia: la fusione estiva su Groenlandia tende ad aumentare quando si instaurano configurazioni di blocking che promuovono cieli più stabili, avvezione meridionale sui settori chiave e un’anomalia di geopotenziale a 500 hPa persistente sull’area groenlandese. Il GBI, definito come media regionale del geopotenziale a 500 hPa su Groenlandia e dintorni, è diventato uno strumento operativo proprio perché “riassume” l’intensità/persistenza del blocco e la sua efficacia nel predisporre condizioni di fusione. E la NAO, come modalità dominante di variabilità del Nord Atlantico, si intreccia spesso con questi assetti: fasi di NAO che favoriscono ondulazioni e creste sul settore groenlandese aumentano la probabilità di afflussi miti/umidi e quindi di fusione.
A scala più “meso-sinottica”, l’analisi per coste e per direzione del vento chiarisce perché la Groenlandia non risponde in modo uniforme. Sulle aree costiere, l’anomalia di total column water è fortemente legata ai giorni di fusione perché lì l’atmosfera marittima può penetrare più facilmente e il bilancio energetico superficiale è spesso vicino alla soglia di fusione: piccoli incrementi di temperatura/umidità e di radiazione infrarossa possono “spostare” molte giornate da non-fusione a fusione. Inoltre, in estati specifiche, anomalie positive dei venti meridionali coincidono con anomalie positive di NMD: è il segnale classico dell’advezione calda e umida lungo i margini della calotta. Questo schema è molto coerente con la letteratura sugli atmospheric rivers (AR), che rappresentano i casi più organizzati e intensi di trasporto di vapore: in vari lavori si mostra che un trasporto anomalo di umidità, spesso in connessione con un blocco, aumenta la probabilità di eventi di fusione estesi e di impatti sul bilancio di massa superficiale.
La dicotomia 2012–2013 è un laboratorio naturale particolarmente istruttivo. Il 2012 è stato caratterizzato da condizioni atmosferiche favorevoli a trasporti caldi/umidi e a configurazioni sinottiche persistenti che hanno facilitato la fusione su larghissima scala; il 2013, al contrario, ha mostrato una stagione relativamente sfavorevole alla fusione diffusa, con una marcata differenza nella frequenza/origine delle masse d’aria che raggiungono aree sensibili (come la Groenlandia nord-orientale). Studi specifici sul contrasto tra 2012 e 2013 confermano quanto sia determinante la circolazione atmosferica nel modulare l’estensione della fusione osservata da satellite.
Un passaggio molto importante del lavoro è quello relativo alla costa orientale: qui anomalie positive di NMD compaiono con downslope winds più forti, suggerendo un ruolo dominante del riscaldamento per subsidenza/foehnrispetto alla semplice avvezione fredda dall’interno della calotta. In pratica, quando il flusso scavalca l’orografia e discende verso la costa, l’aria può riscaldarsi adiabaticamente, ridurre l’umidità relativa e intensificare i flussi turbolenti: una combinazione capace di produrre fusione anche in un contesto in cui l’aria “di origine” potrebbe non sembrare estremamente calda. Questo filone è oggi molto consolidato, soprattutto per la Groenlandia nord-orientale, dove si è mostrato che gli eventi più intensi possono nascere da una “catena” dinamica: AR e trasporto di vapore su un settore (spesso nord-occidentale) che, con un certo lag temporale, innesca condizioni favorevoli a foehn e fusione sul nord-est.
Infine, l’aspetto forse più elegante sul piano diagnostico è l’associazione tra la massima estensione di fusione e la frequenza di masse d’aria di origine nord-orientale a Summit, che nel lavoro risulta la variabile con correlazione negativa più forte con la massima estensione. È una lettura coerente con la fisica: un contributo più frequente da nord-est tende a implicare aria relativamente più fredda e/o meno ricca di vapore (o comunque meno “melt-efficient”), riducendo la probabilità di episodi sinottici capaci di portare la calotta verso la fusione estesa. In termini di dinamica atmosferica, molti studi recenti sulle grandi ondate di fusione ad alta quota mostrano che la firma ricorrente è una cresta in quota e un pattern di alta pressione che favorisce traiettorie meridionali e “intrusioni” di aria mite/umida: quando questo pattern è assente e dominano apporti più settentrionali/orientali, la fusione massima tende a essere meno probabile.
1. Introduzione
La calotta glaciale della Groenlandia (GrIS) è la maggiore riserva di ghiaccio e neve dell’emisfero nord e, se fusa completamente, corrisponde a un potenziale innalzamento del livello medio globale del mare di circa 7,4 m: un ordine di grandezza che da solo basta a spiegare perché ogni variazione nella fusione superficiale sia un segnale climatico ad alta “sensibilità” e ad alta rilevanza socio-ambientale. In effetti, osservazioni e ricostruzioni basate su sensori microonde e modellistica regionale mostrano un aumento statisticamente significativo della fusione superficiale dalla fine del XX secolo, con un’accelerazione evidente negli anni 2000 e con estati “spartiacque” come il 2007, quando l’estensione/quantità di fusione raggiunse valori eccezionali rispetto alle decadi precedenti.
La letteratura degli ultimi due decenni converge su un punto: gli anni di fusione intensa sulla GrIS non sono spiegati soltanto da un riscaldamento “di fondo”, ma soprattutto dalla configurazione della circolazione estiva nord-atlantica che predispone advezione calda e umida verso la Groenlandia e condizioni radiative favorevoli alla fusione. In termini di teleconnessioni, ricorre spesso una combinazione tra NAO estiva negativa e blocking groenlandese, cioè una persistente anomalia anticiclonica (e quindi di geopotenziale) centrata o prossima alla calotta, capace di deviare e ondulare il getto polare e di canalizzare flussi meridionali lungo il margine occidentale. Tuttavia, la risposta della GrIS non è “monotona”: in certe estati il segnale NAO/GBI viene modulato da altri pattern euro-atlantici (per esempio il ramo EA/WR e la componente scandinava), che riorganizzano i centri d’azione della pressione e quindi cambiano traiettorie e caratteristiche termo-igrometriche delle masse d’aria in ingresso.
Dentro questo quadro dinamico entra con forza il tema delle nubi e del bilancio energetico superficiale (SEB). Sulla Groenlandia, le nubi non sono un semplice “disturbo”: a seconda di quota, fase microfisica e spessore ottico, possono raffreddare schermando la radiazione solare oppure scaldare aumentando la radiazione infrarossa discendente (downwelling longwave). Gli eventi estremi mostrano quanto questa competizione sia decisiva: nel luglio 2012, l’estensione eccezionale della fusione fu amplificata da nubi basse contenenti acqua liquida che incrementarono la componente infrarossa verso il suolo e favorirono temperature prossime alla superficie più elevate, in un contesto già predisposto da blocking e avvezione mite. Studi successivi hanno generalizzato il concetto evidenziando che regimi nuvolosi e interazioni nube-strato limite-superficie possono modulare i flussi radiativi e turbolenti e, quindi, la “traduzione” della circolazione sinottica in energia effettivamente disponibile per fondere.
Un secondo aspetto chiave è come misuriamo la fusione. Molti lavori storici si sono basati sull’estensione della fusione (melt extent) derivata da microonde passive: è un indicatore robusto, ottimo per identificare episodi rapidi e diffusi, ma non sempre descrive bene la persistenza locale della fusione lungo la stagione. Per questo, l’uso del numero di giorni di fusione (NMD) è concettualmente potente: integra la componente di durata senza dover inseguire definizioni talvolta ambigue di inizio/fine stagione (quando congelamento e fusione si alternano) e, soprattutto, permette di risolvere meglio la variabilità spaziale, distinguendo “dove” e “quanto spesso” si fonde, non solo “quanto lontano” arriva la fusione in un singolo picco. In parallelo, la quantificazione della quantità di fusione o del runoff richiede modelli (PDD o SEB completi): strumenti indispensabili, ma più esposti a incertezze legate a parametrizzazioni e forcing, in particolare sui flussi turbolenti e sulla fisica dello strato limite in ambiente polare.
Da qui nasce la logica della comparazione multi-metrica adottata nello studio: affiancare l’NMD (stimabile in modo relativamente diretto tramite temperatura superficiale) a indicatori satellitari di estensione massima e cumulativa della fusione, usando il dataset MEaSUREs Greenland Surface Melt Daily basato su SSM/I-SSMIS, e confrontare il tutto con una stima modellistica della fusione stagionale prodotta da un RCM (HIRHAM5) impiegato nelle simulazioni di bilancio di massa superficiale. Questo approccio è particolarmente utile perché molte conclusioni “classiche” sugli anni estremi derivano da un singolo descrittore (spesso l’estensione massima), mentre la fisica sottostante può differire tra un evento di pochi giorni e una stagione complessivamente molto fusiva: la prima è più sensibile a picchi sinottici e a soglie termiche, la seconda alla frequenza/persistenza dei regimi di trasporto caldo-umido e alla modulazione radiativa (inclusi i regimi nuvolosi). Non a caso, lavori recenti mostrano anche il ruolo di trasporti organizzati di vapore (atmospheric rivers) nel forzare la fusione attraverso combinazioni spazialmente variabili di avvezione, nubi e termini del SEB.
In questa cornice, l’obiettivo diventa quindi leggere quindici estati (2000–2014) non come una sequenza di “casi studio” isolati, ma come un campione in cui far emergere relazioni statistiche tra variabili di fusione (NMD, estensioni massima/cumulativa, fusione modellata) e forzanti atmosferiche (circolazione su larga scala, origine delle masse d’aria, contenuto d’acqua in colonna, temperatura troposferica). È un passaggio metodologico importante perché aiuta a capire quando e perché un pattern sinottico si traduce in fusione diffusa: per esempio, mettendo in evidenza dove la correlazione tra contenuto d’acqua in colonna e NMD è massima (fasce costiere specifiche) e come la frequenza di masse d’aria da nord-est a Summit possa associarsi a una riduzione della massima estensione di fusione, con un contrasto particolarmente istruttivo tra l’estate estrema 2012 e la stagione molto più debole del 2013.
Sintesi interpretativa dei risultati riportati in Tabella 1 (studi sugli estremi di fusione in Groenlandia, 1998–2012)
La tabella mette a confronto quindici filoni di analisi che, tra fine anni ’90 e il 2012, hanno definito e spiegato gli “estremi” di fusione della calotta groenlandese usando metriche e strumenti diversi, mostrando in modo molto chiaro che la narrativa scientifica sugli eventi di melt dipende in modo sensibile da (i) che cosa si intende per fusione (estensione, durata, frequenza, runoff, esposizione di ghiaccio nudo, ecc.), (ii) come la si misura (microonde passive, MODIS, scatterometro, carote di firn, modelli regionali), e (iii) quale catena causale si assume dominante (dinamica atmosferica, feedback radiativi/di albedo, precondizionamento nivologico-firnico, trasporti di umidità organizzati). Nelle prime annate considerate, come il 1998, la fusione viene spesso descritta tramite melt extent ricavata da combinazioni di modello climatico regionale e microonde satellitari (Fettweis, Tedesco et al., 2011), ma in alcuni casi l’interpretazione fisica non viene esplicitata (“not given”), un dettaglio non banale perché segnala come, in quella fase, la priorità fosse consolidare serie osservazionali e coerenza tra prodotti più che attribuire in modo robusto la causa sinottica.
Con il 2002, la tabella mostra un salto di qualità nell’attribuzione: l’estremo viene letto attraverso MODIS IST(temperatura superficiale da sensore ottico/IR) e associato a avvezione calda da sud-est e cieli sereni (Hall et al., 2013), oppure a blocking atmosferico sulla Groenlandia con anomalie termiche positive alla superficie e in libera atmosfera usando MODIS e rianalisi NCEP/NCAR (Häkkinen et al., 2014). Nello stesso contesto compaiono studi che introducono una dimensione “metodologica” cruciale: Mernild et al. (2011) combinano microonde (SMMR/SSM/I), SnowModel e rianalisi per descrivere non solo l’estensione massima ma anche frequenza, durata, runoff e ablazione, evidenziando anche la possibilità che estensioni osservate molto elevate possano risentire di aspetti algoritmici del rilevamento microonde (melt extent osservata > melt extent modellata). Steffen et al. (2004), passando da “quanto” fonde a “quando e per quanto”, collegano invece l’estremo a inizio precoce e stagione di fusione lunga, con conseguente riduzione dell’albedo: qui la causalità non è solo dinamica (pattern atmosferico), ma anche di memoria della superficie, perché la perdita di riflettività rende più efficiente l’assorbimento radiativo e amplifica la fusione successiva.
Per gli anni successivi (2003, 2005, 2007) la tabella evidenzia che molti lavori restano centrati su melt extent e su metodi misti (microonde + modelli regionali), ma con attribuzioni talvolta non dichiarate (Fettweis, Tedesco et al., 2011) o più esplicite quando si introducono temperature prossime al suolo: McGrath et al. (2013) attribuiscono per esempio l’estremo 2005 a riscaldamento delle temperature near-surface. In parallelo compaiono interpretazioni che legano la fusione a condizioni pregresse invernali: Frauenfeld et al. (2011) richiamano una NAO invernale positivacapace di ridurre l’accumulo nevoso (meno neve = albedo più basso e transizione più rapida verso ghiaccio nudo) e, insieme a temperature costiere estive elevate, predisporre condizioni favorevoli alla fusione. In questa parte della tabella emerge un messaggio trasversale: il melt estremo può essere innescato da un episodio sinottico caldo, ma spesso diventa “estremo” perché la superficie è già stata precondizionata (nevicate ridotte, neve più sottile, firn più caldo, albedo compromesso).
Un elemento di raccordo temporale importante è lo studio di Mote (2007), che non si focalizza su un singolo anno ma segnala un rafforzamento della fusione nel periodo 2002–2006 durante giugno–agosto, interpretandolo come riduzione del cold content della neve (minor “capacità frigorifera” del manto) e incremento delle temperature del firn: è un’indicazione forte di come, oltre alla variabilità interannuale, esista un cambiamento di fondo nelle condizioni termo-fisiche che rende la calotta più pronta a fondere quando arrivano configurazioni atmosferiche favorevoli.
Nel 2010 la tabella diventa molto più ricca perché entrano in gioco metriche più “di impatto” (runoff, esposizione di ghiaccio nudo, bilancio di massa superficiale) e metodi integrati (MODIS albedo + stazioni automatiche + modello regionale come MAR): Tedesco et al. (2011) indicano temperature near-surface sopra la norma e feedback di albedo come meccanismi principali. Mernild et al. (2011), usando ancora microonde + SnowModel + rianalisi, interpretano l’estremo come risultato di avvezione calda da sud in primavera-estate, con precondizionamento dovuto a un inverno caldo e secco (basse precipitazioni) e un possibile legame statistico con la variabilità multidecadale atlantica. Hanna et al. (2013) e McGrath et al. (2013), introducendo con più decisione anche il runoff e le temperature misurate, sottolineano un punto metodologico cruciale: l’estensione della fusione non sempre è la variabile più sensibile alla circolazione rispetto a runoff/quantità; in altre parole, può esistere una stagione in cui l’area che “va in melt” aumenta moderatamente ma l’acqua prodotta/defluita aumenta molto (o viceversa), a seconda del bilancio energetico e dello stato della superficie.
Nel 2011 e soprattutto nel 2012, la tabella concentra numerosi studi perché il 2012 è l’evento emblematico di fusione eccezionalmente estesa. McGrath et al. (2013) rimarcano l’importanza del riscaldamento near-surface e, per il 2012, richiamano un ruolo di precondizionamento pluriennale, con migrazione verso l’alto della linea di equilibrio e della dry-snow line (segnali coerenti con una maggiore predisposizione della calotta a estendere la fusione verso quote più elevate). Hall et al. (2013) evidenziano l’albedo basso come fattore amplificante. Häkkinen et al. (2014) insistono sul blocking groenlandese e sulle anomalie termiche positive su più livelli. Hanna et al. (2014) ricostruiscono la catena sinottica classica dell’evento: configurazione NAO/GBI, getto polare e avvezione calda meridionale lungo la costa occidentale. Bennartz et al. (2013) aggiungono il tassello radiativo-microfisico decisivo: nubi basse con gocce liquideche aumentano la radiazione infrarossa discendente e quindi la temperatura vicino alla superficie, trasformando un pattern dinamico favorevole in un bilancio energetico nettamente positivo per la fusione. Nghiem et al. (2012) sottolineano il carattere “precursore” di un GBI molto forte prima dell’evento, cioè una persistenza del blocco tale da rendere probabile una risposta di fusione su larga scala. Neff et al. (2014) chiamano in causa il trasporto di aria calda e umida dall’Atlantico verso la costa occidentale favorito dalla formazione di atmospheric rivers, cioè corridoi di trasporto di vapore altamente organizzati che possono aumentare simultaneamente temperatura, contenuto d’acqua in colonna e copertura nuvolosa “riscaldante”. Infine, Keegan et al. (2014), con una prospettiva paleoclimatica e deposizionale (carote di firn + scatterometro + MODIS), propongono un ponte concettuale fra 2012 e 1889 e introducono con forza il ruolo di incendi e black carbon: l’aumento delle temperature favorirebbe più incendi, maggiore deposizione di particolato assorbente e quindi un feedback positivo di albedo (superficie più scura → più assorbimento → più fusione), ampliando il ventaglio di meccanismi oltre la sola dinamica atmosferica.
Letta nel suo insieme, la tabella racconta che gli estremi di fusione groenlandesi non sono spiegati da un singolo “interruttore”, ma da combinazioni ricorrenti: (1) un pattern di circolazione che favorisce advezione calda/umida e spesso blocking, (2) processi radiativi che amplificano l’energia disponibile alla superficie, in particolare attraverso nubi basse liquide e incremento della radiazione infrarossa, (3) feedback di albedo (ghiaccio nudo, neve più sottile, impurità come black carbon) e (4) precondizionamento del manto nevoso/firn che determina quanta energia serve per portare la superficie a 0 °C e mantenerla in fusione. Il valore aggiunto dei “dati della tabella” è proprio questo: mostrare come, a parità di anno estremo, cambino diagnosi e attribuzione a seconda che si guardi all’estensione massima, alla durata, al runoff o ai controlli radiativi, e come l’evento del 2012 rappresenti la convergenza più completa di dinamica, umidità, nubi e feedback superficiali in un unico episodio.
2.1. Prodotto MODIS della temperatura superficiale del ghiaccio (IST)
In questo studio l’identificazione della fusione viene costruita a partire dalla serie di Ice Surface Temperature (IST)derivata da MODIS e disponibile per il periodo 2000–2014 (Hall et al., 2012, 2013). La scelta non è casuale: l’IST MODIS rappresenta una misura della skin temperature della superficie glaciale, cioè della temperatura radiativa degli strati più esterni del manto nevoso/ghiacciato, ed è quindi direttamente legata alla condizione fisica più rilevante per la diagnosi di melt, ovvero il raggiungimento (o la prossimità) del punto di fusione. L’analisi utilizza esclusivamente i pixel che ricadono sulla calotta groenlandese secondo la maschera ghiaccio–terra inclusa nel prodotto, in modo da evitare contaminazioni da superfici rocciose, aree costiere miste o porzioni marginali non glaciali; sia la maschera sia l’IST operano a risoluzione 1,56 × 1,56 km, un dettaglio metodologicamente importante perché consente di descrivere gradienti termici e transizioni spaziali della fusione con una granularità molto più fine rispetto ai prodotti microonde, tipicamente più grossolani ma più robusti in presenza di nubi.
Come spesso accade per le temperature superficiali da sensori in infrarosso termico, anche l’IST MODIS tende a presentare un bias freddo rispetto a misure in situ, aspetto già discusso nelle validazioni precedenti (Hall et al., 2008). Un confronto più recente con osservazioni al suolo (Shuman et al., 2014) conferma però un punto cruciale per l’impianto di questo lavoro: vicino a 0 °C (cioè nella fascia termica che decide “melt sì/melt no”) il prodotto risulta strettamente compatibile con le misure ground-based. Il bias diventa via via più pronunciato al diminuire della temperatura, con una sottostima che parte da circa −0,5 °C per condizioni prossime allo zero e può arrivare fino a circa −5 °C intorno a −60 °C. Dal punto di vista fisico-strumentale, questa non linearità è coerente con limiti ben noti della termometria satellitare in ambiente criosferico: l’accuratezza dipende dalla stima dell’emissività, dalla correzione atmosferica (in particolare in condizioni di aria molto secca e fredda), dalla sensibilità ai residui di nube o foschia glaciale e, soprattutto, dal fatto che MODIS osserva la temperatura radiativa dello strato superficiale, che può disaccoppiarsi dalla temperatura dell’aria a 2 m in presenza di forti inversioni o stratificazione stabile. Proprio perché l’obiettivo qui è la diagnosi di melt e non la ricostruzione accurata della climatologia del freddo estremo, gli autori considerano tale bias un limite secondario: la metrica di interesse è centrata sulle condizioni prossime alla fusione, dove la qualità del prodotto è migliore e dove un criterio basato sull’IST è più informativo rispetto a variabili indirette.
Un ulteriore vantaggio operativo è legato al timing di osservazione: l’IST viene fornita alle 17:00 UTC, che corrispondono approssimativamente alle 13:00–16:00 di tempo solare locale in Groenlandia. Questa finestra cade in prossimità delle ore mediamente più calde del ciclo diurno, quando la probabilità di superare la soglia di fusione (o di raggiungere valori prossimi allo zero) è massima; di conseguenza, il prodotto è particolarmente adatto a intercettare i giorni in cui la superficie entra effettivamente in melt o vi si avvicina in modo critico. In termini di diagnostica, ciò rende l’IST MODIS un indicatore efficace per costruire conteggi stagionali di giorni di fusione e per descrivere la variabilità spaziale degli episodi di melt, mantenendo una coerenza con la letteratura che ha già impiegato lo stesso record e ne ha discusso in modo dettagliato prestazioni e limiti (Hall et al., 2008, 2012, 2013; Shuman et al., 2014).
Figura 1 | Disponibilità del MODIS IST e climatologia dei giorni di fusione estivi (JJA, 2000–2014)
La Figura 1 serve, di fatto, a “mettere in sicurezza” l’intera analisi: prima mostra quanto il dataset MODIS di temperatura superficiale del ghiaccio (IST) sia effettivamente disponibile durante l’estate groenlandese, poi fa vedere quale struttura spaziale assume, in media, la fusione estiva quando la si diagnostica come numero di giorni di fusione (NMD). Nel pannello (a) la scala è una percentuale: indica per ogni pixel della calotta quante giornate di giugno–agosto, mediamente tra 2000 e 2014, presentano una stima valida di IST. Questa informazione è cruciale perché MODIS lavora nel canale infrarosso termico: la misura è molto informativa vicino al punto di fusione, ma è intrinsecamente vulnerabile alla copertura nuvolosa (e, più in generale, a condizioni atmosferiche che impediscono una “vista” chiara della superficie). Il disegno che emerge è tipico della Groenlandia estiva: un’ampia parte della calotta mostra disponibilità elevata, mentre alcune aree – soprattutto verso il sud e lungo settori costieri più esposti a ingressi umidi atlantici – evidenziano disponibilità più bassa, coerente con una maggiore frequenza di nubi e sistemi perturbati. In altre parole, la mappa non sta dicendo “qui la temperatura è più bassa”, ma “qui il satellite riesce meno spesso a misurare la temperatura della superficie”, e questo è un passaggio metodologico fondamentale perché, dove l’atmosfera è spesso nuvolosa, i giorni mancanti possono introdurre una tendenza a sottostimare (o quantomeno rendere più incerta) la frequenza della fusione se non si gestiscono bene i vuoti osservativi. È anche il motivo per cui gli autori insistono sulla bontà del prodotto proprio vicino a 0 °C: le validazioni precedenti (Hall et al., 2008; Hall et al., 2012, 2013; Shuman et al., 2014) mostrano che l’eventuale bias freddo del satellite diventa davvero problematico alle temperature molto basse, mentre nel “range decisionale” della fusione la coerenza con misure in situ è molto migliore, rendendo l’IST un proxy robusto per contare i giorni in cui la superficie entra (o sfiora criticamente) la condizione di melt.
Il pannello (b) traduce questa base dati in climatologia fisica: il numero medio di giorni di fusione in JJA disegna una geometria molto riconoscibile, con un “anello” periferico di NMD elevati e un interno della calotta dove la fusione estiva è rara o praticamente assente su base media. È la firma combinata di quota e bilancio energetico: lungo i margini, a basse quote, la superficie è più spesso prossima a 0 °C e basta un forcing relativamente modesto (più radiazione netta o un episodio di avvezione mite) per produrre fusione frequente; sull’altopiano centrale, invece, le temperature medie estive restano in genere sufficientemente basse da impedire che il bilancio energetico superi spesso la soglia di fusione, quindi l’NMD collassa verso valori molto piccoli. Questa distribuzione è perfettamente in linea con ciò che sappiamo dai modelli di bilancio energetico e dagli studi osservativi: nelle zone di ablazione costiera diventano dominanti, a seconda delle situazioni sinottiche, l’aumento dell’assorbimento di radiazione solare quando affiora ghiaccio nudo e si attiva il feedback dell’albedo (Tedesco et al., 2011; Box et al., 2012), ma anche i contributi di radiazione infrarossa e flussi turbolenti quando arrivano masse d’aria più calde e umide. Il fatto che gli NMD massimi si addensino soprattutto nei settori meridionali e lungo porzioni della fascia costiera è quindi il risultato integrato di (i) condizioni termiche di base più favorevoli alla fusione, (ii) maggiore probabilità di avvezioni miti/umide dal Nord Atlantico e (iii) una superficie che, una volta “aperta” (neve assottigliata o ghiaccio esposto), risponde con maggiore efficienza radiativa.
Il valore scientifico della figura sta soprattutto nella lettura “a due livelli” che propone. Da un lato, (a) ti dice dove il prodotto MODIS è più affidabile dal punto di vista della continuità campionaria e dove invece la nuvolosità può creare buchi; dall’altro, (b) ti mostra che la fusione estiva media non è un fenomeno diffuso uniformemente, ma un processo fortemente confinato ai margini, con un gradiente netto verso l’interno. È un punto chiave anche per interpretare gli anni estremi (come il 2012): quando la circolazione estiva favorisce blocking groenlandese e/o una fase NAO estiva negativa, il getto e i campi di pressione tendono a predisporre avvezioni meridionali lungo la Groenlandia occidentale e a modificare nubi e flussi radiativi, così che la “fascia di NMD elevati” può espandersi verso quote più alte e regioni più interne (Fettweis et al., 2013; Hanna et al., 2014). In quel tipo di configurazione, la microfisica delle nubi può diventare decisiva: l’aumento della radiazione infrarossa discendente associata a nubi basse con acqua liquida è stato identificato come un amplificatore importante della fusione del 2012 (Bennartz et al., 2013), e proprio per questo la figura 1, pur essendo “solo” climatologica, prepara il terreno: mostra qual è la struttura media su cui poi agiscono, anno per anno, le anomalie dinamiche e radiative. In sintesi, (a) ti dice quanto è solida la base osservativa estiva, (b) ti fa vedere la geografia climatologica della fusione: un fenomeno marginale per definizione, che diventa “centrale” solo quando la circolazione atmosferica riesce a spingere calore e umidità (e quindi energia disponibile) verso l’interno della calotta.
2.2. Integrazione del MODIS IST con la Skin Temperature di ERA-Interim
Un limite strutturale dei prodotti MODIS di temperatura superficiale del ghiaccio (IST) è la presenza di lacune giornaliere locali dovute alla copertura nuvolosa: quando le nubi schermano la superficie, il sensore satellitare non riesce a stimare la temperatura radiativa e, in quel punto, la classificazione “fusione / non fusione” diventa impossibile. La Figura 1a quantifica bene questo problema nel trimestre estivo (giugno–agosto) 2000–2014: nel Nord-Est groenlandese l’IST risulta disponibile per oltre il 70% dei giorni estivi, mentre nel Sud la disponibilità scende tipicamente a circa 40–50%, coerentemente con una climatologia più spesso nuvolosa e con maggiore influenza delle perturbazioni atlantiche. Questa disomogeneità non è un dettaglio tecnico: se manca la misura, non solo non si può “vedere” la fusione, ma non si può nemmeno calcolare in modo univoco il numero esatto di giorni di fusione, perché un giorno classificato come “assenza di fusione” potrebbe in realtà essere un giorno “non osservato” (cioè nuvoloso) in cui la fusione è avvenuta ma non è stata rilevata. Proprio su questo punto, Hall et al. (2013) avevano proposto una soluzione più prudente ma anche più grossolana: una stima categoriale dei giorni di fusione, costruendo mappe annuali in quattro classi (nessuna fusione, 1 giorno, 2 giorni, più di 2 giorni/anno). Tuttavia, come già discusso in Hall et al. (2012, 2013), la classe “no melt” può inglobare anche giornate nuvolose senza dato IST, quindi con potenziale sottostima degli episodi di fusione.
Per ridurre questo problema di “missing data” in condizioni nuvolose, gli autori integrano l’IST con la skin temperature (SKT) fornita dalla rianalisi ERA-Interim di ECMWF (Dee et al., 2011). L’idea è metodologicamente chiara: l’IST MODIS è un prodotto pensato per la qualità climatica della stima radiativa della superficie, ma soffre di buchi spaziali/temporali; la SKT di ERA-Interim, al contrario, garantisce copertura completa perché deriva dal bilancio energetico superficiale modellato (quindi può essere calcolata anche “sotto le nubi”), pur portandosi dietro le incertezze tipiche della modellistica della superficie in ambiente polare (parametrizzazioni dei flussi turbolenti, nubi, radiazione, conduzione in neve/firn). Dal punto di vista operativo, ERA-Interim ha una risoluzione orizzontale nativa di circa 0,73° in lat-lon; qui si usano valori di SKT interpolati a 0,125°, un passaggio utile per la co-locazione con la griglia di lavoro, pur sapendo che l’interpolazione non crea informazione fisica nuova ma riduce soprattutto i disallineamenti geometrici. Per massimizzare la coerenza temporale con MODIS (che campiona nel primo pomeriggio locale), la SKT viene presa alle 18 UTC, e la regola è semplice: la fusione viene diagnosticata primariamente con IST, e solo quando IST manca si sostituisce con la SKT allo stesso orario e nella stessa posizione.
Proprio perché le due grandezze nascono da catene di produzione diverse (radiometria satellitare con schermatura nuvolosa e correzioni atmosferiche, contro rianalisi con modello di superficie e forcing assimilati), IST e SKT non coincidono e la loro differenza va valutata. Per questo gli autori confrontano sistematicamente i due campi durante marzo–agosto 2000, portandoli sulla stessa griglia e usando solo i punti in cui entrambi sono disponibili, calcolando un bias definito come SKT − IST. Nel periodo complessivo il bias può arrivare fino a 3,4 °C, ma soprattutto emerge una dipendenza fisicamente plausibile: più la temperatura è alta, più il bias tende a ridursi, segno che l’accordo migliora proprio nel range termico più critico per la fusione. Quando si restringe l’analisi ai casi in cui la fusione è “potenzialmente” in gioco, cioè con IST sufficientemente vicino alla soglia, il quadro cambia sensibilmente: all’inizio della stagione di fusione (aprile–maggio) il bias medio nelle aree che mostrano condizioni compatibili con fusione è intorno a ~2 °C, mentre in giugno–agosto il bias sui punti in fusione (diagnosticati con IST) resta in genere positivo ma per lo più inferiore a 2 °C. Questo comportamento è importante perché delimita l’errore atteso quando si usa SKT come “riempitivo” dei buchi MODIS: l’aggiunta della rianalisi riduce il problema dei giorni non osservati, ma introduce un potenziale spostamento termico che va gestito con soglie coerenti.
La scelta delle soglie è infatti un passaggio concettuale centrale. In linea con Hall et al. (2012, 2013), la fusione viene identificata per ciascun pixel MODIS quando IST supera −1 °C (e non 0 °C), perché in presenza di neve/ghiaccio è possibile avere fusione interna/subsuperficiale anche con temperatura superficiale “vera” leggermente sotto lo zero: è una conseguenza della struttura termica del manto e del fatto che la diagnostica radiativa può non catturare perfettamente l’inizio dei processi di melt all’interno del firn (Kuipers Munneke et al., 2009; Van den Broeke et al., 2008; Hall et al., 2012). Per la SKT, invece, viene adottata una soglia di 0,0 °C, motivata dal fatto che in piena estate la SKT tende a risultare più alta dell’IST: quindi, quando IST supera la soglia indicativa di possibile fusione interna (−1 °C), è ragionevole attendersi che SKT superi 0 °C. Gli autori segnalano però un punto sottile ma operativo: a inizio stagione la SKT può rimanere chiaramente sotto 0 °C mentre l’IST indica già condizioni di melt nella stessa area, evidenziando una divergenza che verrà discussa più avanti (sezione 2.7) e che richiama la maggiore difficoltà dei modelli di superficie nel catturare transizioni precoci legate a stratificazione stabile, metamorfismo della neve e rapidi cambiamenti di albedo/irraggiamento.
A questo punto il calcolo dell’NMD (Number of Melt Days) è costruito in modo diretto: per ogni giorno e per ogni punto grigliato si verifica se l’IST (oppure, quando mancante, la SKT) supera la rispettiva soglia di fusione; sommando nel tempo si ottiene il numero di giorni di fusione. La climatologia media estiva 2000–2014 (Figura 1b) mostra il risultato atteso e fisicamente coerente: i valori massimi di NMD si trovano sul margine meridionale della calotta, con valori prossimi a 90, indicando aree dove la fusione si verifica praticamente ogni giorno d’estate. In sintesi, questa sezione formalizza una strategia “ibrida” molto pragmatica: preservare l’informazione ad alta risoluzione e prossima alla soglia di MODIS quando disponibile, ma ridurre la sottostima indotta dalle nubi usando una rianalisi (Dee et al., 2011) come supporto, dopo averne quantificato il bias nel range termico più rilevante per la diagnostica della fusione.
2.3. Variabili atmosferiche basate su rianalisi
Per descrivere in modo coerente le condizioni meteorologiche sopra la Groenlandia e nelle aree circostanti, gli autori si appoggiano alla rianalisi ERA-Interim, scelta perché ampiamente valutata alle alte latitudini e considerata, nel complesso, affidabile nel riprodurre la variabilità sinottica e molte variabili near-surface su ghiaccio marino (pur con bias non trascurabili nello strato limite). In studi di validazione sul Mar Glaciale Artico, per esempio, ERA-Interim risulta spesso tra i prodotti meglio posizionati rispetto ad altre rianalisi per temperatura e vento, anche se può mostrare bias caldo e/o umido nel boundary layer (ordine di 1–2 K o più in specifiche condizioni), un punto che impone prudenza soprattutto quando si interpretano variabili prossime alla superficie. In parallelo, viene sottolineato un elemento metodologico importante: tra le rianalisi globali “classiche”, ERA-Interim e JRA-55 sono (storicamente) quelle che impiegano l’assimilazione variazionale 4D-Var, una caratteristica considerata vantaggiosa in regioni povere di osservazioni, perché consente di ottimizzare l’analisi nello spazio e nel tempo in modo più consistente con la dinamica del modello. Nel contesto artico, inoltre, la capacità di rappresentare trasporto e bilanci di umidità è particolarmente critica, perché nubi e contenuto di vapore modulano direttamente il bilancio radiativo e quindi la probabilità di fusione; su questo aspetto, lavori di confronto fra reanalisi evidenziano differenze sostanziali nella chiusura e nella scomposizione del ciclo idrologico artico, confermando che la scelta del prodotto di rianalisi non è un dettaglio “tecnico” ma una componente del segnale fisico che si vuole interpretare.
Su questa base, lo studio utilizza medie mensili (e non il dato istantaneo “grezzo”) di un set di variabili atmosferiche che coprono sia la forzante sinottica sia i termini più direttamente collegati a avvezione, stabilità e umidità: componenti del vento a 10 m e a 500 hPa (per distinguere il trasporto nel basso strato e la circolazione di guida in media troposfera), altezze di geopotenziale a 500 e 850 hPa (Z500 e Z850, utili per diagnosticare creste, depressioni e pattern di blocking e per legare la fusione ai regimi di circolazione), temperatura dell’aria a 2 m e ai livelli 850/500 hPa (per separare il segnale termico near-surface da quello troposferico), Total Column Water (TCW) integrato verticalmente (includendo vapore, acqua e ghiaccio di nube e contributi legati a precipitazione), e la velocità verticale a 500 hPa(proxy dinamico di subsidenza/ascesa sinottica, quindi utile per collegare ad esempio condizioni anticicloniche a compressione adiabatica, cielo più stabile e/o soppressione della convezione nuvolosa, oppure fasi cicloniche a sollevamento e maggiore copertura nuvolosa). Tutti questi campi vengono usati alla risoluzione di uscita interpolata a 0,125°, in continuità con quanto fatto per la SKT: l’idea non è “inventare” dettaglio fine che la rianalisi non possiede, ma garantire una griglia di lavoro comune e una co-locazione stabile tra variabili atmosferiche e diagnostiche di fusione.
Un passaggio chiave della sottosezione è la scelta di non usare direttamente, nell’analisi finale, il prodotto ERA-Interim di copertura delle low clouds (LCC). Gli autori lo valutano, ma riscontrano che spiega poco le anomalie delle variabili di melt; la spiegazione proposta è coerente con un problema ben noto: la frazione di nube “diagnostica” in molte rianalisi è una variabile fortemente dipendente da parametrizzazioni e post-processing, e può presentare errori sistematici quando confrontata con osservazioni. Inoltre, anche sul piano fisico, l’effetto radiativo delle nubi basse non è uniforme: nelle aree elevate dell’interno groenlandese, dove il contenuto di acqua liquida può essere ridotto, l’impatto radiativo netto tende a essere più debole rispetto alle fasce costiere, mentre lungo i margini e in specifiche configurazioni sinottiche nubi basse liquide possono invece aumentare sensibilmente la radiazione infrarossa discendente e favorire fusione (come mostrato in modo molto chiaro per il 2012). Per questi motivi, invece di “inseguire” una variabile di cloud cover potenzialmente rumorosa, lo studio preferisce usare TCW come descrittore più integrato del ruolo combinato di vapore + nubi nel predisporre un ambiente radiativamente e termodinamicamente favorevole alla fusione: non è una sostituzione perfetta (TCW non equivale a cloud radiative forcing), ma è spesso più robusta per correlazioni interannuali perché cattura la componente di umidità e condensato che accompagna le intrusioni miti/umide e i regimi sinottici di trasporto.
Infine, per collegare la variabilità della fusione non solo agli indici e ai campi medi, ma anche alla “storia” delle masse d’aria, gli autori introducono un tassello lagrangiano: calcolano traiettorie retrograde di 5 giorni per la stazione Summit (72,58°N; 38,46°W) a 100 m sopra la superficie nevosa, ripetendo il calcolo per ogni giorno estivo (JJA) nel 2000–2014. Le traiettorie sono ottenute tramite il servizio METEX (Meteorological Data Explorer), che consente calcoli sistematici basati su rianalisi/analisi globali; qui vengono utilizzati i prodotti NCEP CFSR e la sua estensione/upgrade legata a CFSv2, così da avere un forcing coerente per la diagnostica dei percorsi e dell’origine delle masse d’aria. Quando serve mettere a fuoco in modo più mirato la Groenlandia nord-orientale, lo stesso schema viene replicato su un punto rappresentativo (77°N; 30°W). Il senso fisico di questa scelta è molto pratico: le correlazioni tra melt e campi medi (Z500, TCW, T850, vento) dicono “quale configurazione” favorisce fusione, mentre le traiettorie dicono “da dove arriva” l’aria associata a quegli anni/episodi, permettendo confronti diagnostici robusti (ad esempio tra estati con masse d’aria più spesso sud-occidentali e altre con contributo più frequente nord-orientale) senza ridurre tutto a un singolo indice di teleconnessione.
2.4. Quantità di fusione modellata con HIRHAM5
Per disporre di una stima “process-based” della produzione di acqua di fusione sulla Groenlandia, lo studio utilizza l’output di una simulazione di bilancio di massa superficiale ottenuta con il modello climatico regionale HIRHAM5(Langen et al., 2017). HIRHAM5 nasce come modello ibrido: combina la dinamica del sistema di previsione numerica HIRLAM (qui nella famiglia HIRLAM7) con schemi fisici derivati dal modello di circolazione generale ECHAM5, un’impostazione pensata per mantenere una descrizione dinamica robusta su area limitata e, al tempo stesso, una fisica atmosferica completa e coerente con la modellistica climatica. Nel caso specifico delle applicazioni sulla GrIS, questo impianto è importante perché consente di forzare ai bordi il dominio con una rianalisi globale (qui ERA-Interim) e di sviluppare internamente, con risoluzione chilometrica, la struttura sinottica e mesoscala che controlla trasporti di calore/umidità, nubi e flussi energetici alla superficie, elementi chiave per la fusione e per l’evoluzione del firn.
Un punto centrale della sottosezione è che la stima della “melt amount” non deriva solo da un bilancio energetico superficiale “semplice”, ma da una versione aggiornata dello schema di superficie/sottosuperficie che incorpora processi idrologici e strutturali del manto nevoso: densificazione della neve/firn, conducibilità idraulica variabile in funzione dello stato della neve, e soprattutto la saturazione irreducibile (la quota d’acqua liquida trattenuta per capillarità nei pori, che non drena facilmente), oltre a scelte che influenzano percolazione e ritenzione dell’acqua liquida e la formazione di ghiaccio sovrapposto/strati impermeabili. La saturazione irreducibile richiamata nel MOD-ref è collegata alla classica parametrizzazione sperimentale di Coléou & Lesaffre (1998), che ha quantificato in laboratorio quanta acqua può rimanere “intrappolata” in neve bagnata a seconda delle proprietà del mezzo poroso; nella pratica modellistica, questa quantità controlla quanta acqua resta nel firn (potenzialmente rifrezza, rilasciando calore latente) e quanta invece prosegue verso il basso o contribuisce a runoff.
Il setup numerico richiamato dagli autori fornisce anche l’ordine di scala della rappresentazione: HIRHAM5 è qui impiegato a circa 5–6 km (0,05° su griglia a polo ruotato), con 31 livelli verticali in atmosfera, e con condizioni ai bordi derivate da ERA-Interim, così da assicurare coerenza con la circolazione osservata nel periodo 1980–2014 e, nello specifico di questo lavoro, nel 2000–2014. In letteratura, l’uso di RCM a risoluzione dell’ordine di 5 km è considerato un compromesso efficace per Groenlandia perché migliora la rappresentazione dell’orografia e dei gradienti costieri rispetto a griglie più grossolane, con ricadute dirette su vento catabatico, advezioni lungo costa e distribuzione spaziale dei flussi energetici—tutti fattori che governano la geografia della fusione.
Nel presente studio viene selezionata la variabile di fusione proveniente dall’esperimento MOD-ref, che include tre scelte fisiche chiave: (i) la percolazione dell’acqua liquida nella neve descritta secondo la legge di Darcy (quindi con un flusso dipendente dalla conducibilità idraulica e dai gradienti idraulici), (ii) la parametrizzazione della saturazione irreducibile secondo Coléou & Lesaffre, e (iii) un’albedo della neve derivata da MODIS, che migliora la coerenza tra stato radiativo della superficie e forzante energetica della fusione. Inoltre, MOD-ref assume che, quando in uno strato nevoso è presente ghiaccio, esso occupi metà dell’area orizzontale dell’unità di griglia, influenzando la conducibilità idraulica effettiva: è un modo pragmatico per rappresentare l’eterogeneità sub-grid (neve/firn con inclusioni di ghiaccio e quindi vie preferenziali o barriere al flusso) senza risolvere esplicitamente la micro-struttura.
È altrettanto importante chiarire cosa rappresenta esattamente la variabile usata qui: gli autori impiegano la fusione modellata di neve e/o ghiaccio glaciale come produzione di melt nel punto di griglia, ma la variabile non include alcuni passaggi idrologici successivi (ad esempio instradamento del runoff tra celle adiacenti o contributi di rifreezing contabilizzati come “destino” dell’acqua nel bilancio idrologico della colonna), perché l’obiettivo è confrontare un indicatore di potenziale di fusione con metriche osservate (NMD ed estensione di fusione) che descrivono soprattutto la presenza e la persistenza del melt, non necessariamente la redistribuzione laterale dell’acqua. In questo senso, il dato HIRHAM5 usato nello studio va interpretato come una misura modellistica coerente della tendenza a produrre acqua di fusione nelle diverse estati, mentre la trasformazione di quella produzione in runoff effettivo dipende da ritenzione, refreezing, stratificazione del firn e dalla presenza di strati impermeabili, processi che Langen et al. discutono proprio come determinanti per l’evoluzione termica e idrologica del manto.
La Figura 2 mette in fila, per ogni estate giugno–agosto (JJA) dal 2000 al 2014, due indici che in letteratura sono diventati quasi “diagnostici” quando si vuole capire perché alcune stagioni groenlandesi producano molta più fusione di altre: il North Atlantic Oscillation Index (NAOI, curva arancione) e il Greenland Blocking Index (GBI, curva blu). Il messaggio fisico è che, quando il sistema nord-atlantico si dispone in modo da favorire un’anomalia anticiclonica persistente sulla Groenlandia (GBI alto), la NAO estiva tende spesso a collocarsi in fase negativa(NAOI basso), cioè con una circolazione meno zonale e più ondulata, un getto polare più “meandriforme” e un assetto sinottico favorevole a subsidenza/ridging sul settore groenlandese e, contemporaneamente, a trasporti meridionali di aria più mite lungo i margini della calotta. Questo legame non è una regola matematica perfetta, ma è un comportamento ricorrente già discusso in diversi lavori che collegano NAO estiva, blocking groenlandese e variabili di melt (ad esempio Box et al., 2012; Fettweis et al., 2013; Hanna et al., 2013, 2014, 2016).
Dal punto di vista dinamico, il GBI è concettualmente una misura dell’intensità del “blocco” in media troposfera (tipicamente costruito a partire dall’altezza di geopotenziale a 500 hPa mediata su un dominio groenlandese): valori più alti implicano un’anomalia di cresta più robusta e spesso più persistente. In estate, una cresta/blocco sopra la Groenlandia tende a favorire (i) subsidenza e stabilizzazione della colonna, (ii) una maggiore probabilità di giornate con bilancio radiativo più favorevole (più radiazione solare netta se la nuvolosità diminuisce e/o cambia regime), e (iii) un assetto dei venti tale da facilitare episodi di warm-air advection sul lato occidentale e sud-occidentale della calotta, cioè proprio le regioni che “rispondono” più rapidamente con aumento dei giorni di fusione e con ampliamento dell’area in melt. Questa catena fisica è coerente con l’interpretazione degli eventi estremi: l’estate 2012, che nella figura spicca per GBI molto elevato e NAOI fortemente negativo, è l’esempio canonico di fusione eccezionale in cui il forcing principale è atmosferico (configurazione di blocking/NAO-) e agisce sia via dinamica (advezioni miti) sia via termini radiativi e di umidità (Hanna et al., 2014; Fettweis et al., 2013). Proprio per il 2012, studi mirati hanno mostrato anche come la microfisica nuvolosa possa amplificare l’energia disponibile alla superficie: nubi basse con contenuto di acqua liquida aumentano la radiazione infrarossa discendente e spingono la temperatura prossima alla superficie verso condizioni più favorevoli al melt, rendendo l’evento non soltanto “caldo”, ma anche radiativamente efficiente (Bennartz et al., 2013). Altri lavori hanno enfatizzato la componente di trasporto di umidità, includendo casi di corridoi di vapore (atmospheric rivers) diretti verso la Groenlandia occidentale che, in presenza di un regime bloccato, possono rafforzare sia l’avvezione calda sia i termini radiativi legati a vapore e nubi (Neff et al., 2014).
La sequenza temporale 2000–2014 mostrata in Figura 2 è utile anche perché evidenzia che la Groenlandia non vive una “modalità” unica: si osserva una marcata variabilità interannuale con passaggi repentini tra estati più bloccate e estati più zonali. In modo qualitativo, la figura suggerisce una fase in cui, dal 2007 al 2012, il GBI tende spesso a mantenersi su valori più alti e il NAOI più spesso su valori negativi: questo è coerente con la narrativa di diversi studi che descrivono, in quel periodo, una maggiore frequenza di configurazioni favorevoli a fusione elevata e a record di melt extent/SMB negativo (Box et al., 2012; Hanna et al., 2013). Il 2013, al contrario, emerge come un anno “opposto” rispetto al 2012, con una brusca inversione del segnale: NAOI positivo e GBI molto più basso, quindi una circolazione estiva meno favorevole a blocchi persistenti sulla Groenlandia e, in media, più sfavorevole a una fusione diffusa e duratura. Questo contrasto 2012–2013 è particolarmente prezioso perché dimostra che, a parità di trend climatico di fondo, la risposta della calotta su scala annuale può essere dominata dalla configurazione sinottica stagionale, cioè dalla posizione/intensità della cresta groenlandese, dall’ondulazione del getto e dalle traiettorie delle masse d’aria (tema che nello studio viene poi affrontato anche con le backward trajectories).
Infine, la figura aiuta a capire un punto metodologico spesso sottovalutato: NAO e blocking sono “parenti stretti”, ma non sono sinonimi. Ci sono anni in cui l’anticorrelazione è evidente (GBI alto ↔ NAOI negativo) e anni in cui la relazione è più sfumata, segnalando che altri pattern euro-atlantici (come East Atlantic pattern o Scandinavian pattern) e la posizione longitudinale delle anomalie possono modulare l’efficacia del blocco nel produrre avvezione calda/umida o nel cambiare il regime nuvoloso sopra la calotta (Lim, 2015). In altre parole, la Figura 2 non è solo un “contorno” descrittivo: è una chiave interpretativa che colloca le estati 2000–2014 lungo un asse dinamico (NAO/GBI) ben documentato in letteratura e che, per la Groenlandia, si traduce direttamente in differenze sostanziali nei giorni di fusione, nell’estensione del melt e nella quantità di acqua prodotta, soprattutto quando la circolazione entra in regimi di blocking persistente come nel 2012.
2.5. Estensione della fusione dal dataset MEaSUREs Greenland Surface Melt
In questa parte gli autori introducono una seconda famiglia di metriche di fusione, complementare ai “giorni di fusione” derivati da MODIS: l’estensione della fusione (melt extent) ricavata dal dataset MEaSUREs Greenland Surface Melt, che fornisce un’indicazione giornaliera della presenza di fusione superficiale sulla Groenlandia con una risoluzione orizzontale di circa 25 km (Mote, 2014). Il principio osservativo è quello classico della criosfera a microonde passive: la fusione viene diagnosticata a partire dalla temperatura di brillanza misurata da radiometri satellitari, sfruttando il fatto che l’emissione a microonde del manto nevoso cambia in modo rapido e marcato quando compare acqua liquida nei pori della neve (Mote & Andersson, 1995). In termini fisici, la presenza anche di piccole frazioni di liquid water modifica le proprietà dielettriche del mezzo, riduce la scattering “da neve asciutta” e tende ad aumentare l’emissività effettiva, producendo un segnale radiometrico facilmente distinguibile rispetto alle condizioni asciutte; questo rende i sensori a microonde particolarmente adatti a rilevare la fusione perché operano indipendentemente dall’illuminazione solare e con minore sensibilità alla copertura nuvolosa rispetto ai sensori infrarossi/ottici. Nella costruzione del record MEaSUREs, le osservazioni provengono da SSM/I fino al 31 dicembre 2007 e da SSMIS successivamente; nello studio vengono usati tre descrittori: (i) l’estensione giornaliera e la sua cumulata lungo l’intera stagione di fusione, e (ii) la massima estensione giornaliera di ciascuna stagione. Questa distinzione è importante perché le tre metriche non raccontano la stessa storia: la massima estensione è sensibile a episodi brevi ma molto intensi (pochi giorni “esplosivi” possono dominare il massimo stagionale), mentre l’estensione cumulativa integra la persistenza e la frequenza della fusione durante tutta la stagione, risultando più affine a variabili che “accumulano” nel tempo (durata della stagione di melt, energia netta disponibile, predisposizione della superficie). Non a caso, la letteratura mostra spesso che le conclusioni sugli “anni estremi” cambiano a seconda che si guardi al massimo giornaliero oppure a una misura integrata: un evento sinottico eccezionale può far impennare il massimo, ma una stagione complessivamente favorevole alla fusione (più giorni con condizioni prossime allo zero, più trasporto caldo-umido, più feedback radiativi/di albedo attivati) tende a emergere meglio nelle metriche cumulative.
Gli autori ricordano anche due aspetti metodologici che, nella pratica, determinano la robustezza delle comparazioni tra studi. Primo: il dataset MEaSUREs è reso disponibile (almeno nella versione utilizzata qui) solo fino al 2012, per cui l’analisi dell’estensione di fusione basata su questo prodotto non può coprire l’intero 2000–2014 come altre variabili dello studio. Secondo: benché la fusione rilevata con SSM/I e SSMIS sia stata confrontata più volte con osservazioni in situ o con fusione modellata (Bhattacharya et al., 2009; Fettweis et al., 2006; Fettweis, Tedesco, et al., 2011; Nghiem et al., 2012; Mote & Andersson, 1995), i risultati non sono sempre direttamente sovrapponibili perché cambiano gli algoritmi di detezione (soglie, frequenze, criteri di classificazione) e perché esiste un problema strutturale di rappresentatività di scala: un sensore a 25 km “vede” una media spaziale che può includere mosaici di neve asciutta, neve bagnata, ghiaccio nudo e zone di rifreezing, mentre una misura al suolo o una stazione automatica è puntuale e può non descrivere la cella satellitare, soprattutto in aree costiere con forti gradienti altimetrici e termici. Proprio per questo gli autori non insistono tanto su una corrispondenza uno-a-uno, quanto su un consenso più generale, ben consolidato nella comunità, secondo cui le microonde passive sono complessivamente molto efficaci nel rilevare la comparsa di acqua liquida nel manto nevoso e quindi nel mappare la fusione su larga scala (Mote & Andersson, 1995), pur con limiti interpretativi noti: la risoluzione relativamente grossolana attenua i dettagli locali, e alcune condizioni (transizioni rapide, pioggia su neve, superfici miste, metamorfismo e refreezing) possono complicare l’attribuzione “melt/no melt” se si pretendono diagnosi puntuali. In questo studio, l’uso dell’estensione MEaSUREs accanto a NMD e alla fusione modellata HIRHAM5 ha quindi un ruolo preciso: fornire una metrica satellitare indipendente, robusta su scala sinottica, che cattura bene la geografia e l’ampiezza areale degli episodi di fusione e che permette di distinguere, in modo più chiaro, tra stagioni dominate da un picco estremo (massimo giornaliero) e stagioni caratterizzate da fusione più persistente (cumulata).
2.6. Indici di circolazione su larga scala
Per verificare se la variabilità del numero di giorni di fusione (NMD) sia legata alla circolazione atmosferica emisferica e, in particolare, al settore nord-atlantico, lo studio adotta due indici “chiave” e li mette in relazione con le serie di melt: il North Atlantic Oscillation Index (NAOI) e il Greenland Blocking Index (GBI) (richiamati in Figura 2). I valori mensili del NAOI provengono dal NOAA Climate Prediction Center e sono costruiti con una procedura standard di monitoraggio teleconnettivo basata su Rotated Principal Component Analysis (RPCA) applicata alle anomalie standardizzate dell’altezza di geopotenziale a 500 hPa sull’emisfero nord (dominio a nord di 20°N), con standardizzazione riferita al periodo base 1950–2000; questa impostazione segue la tradizione metodologica di Barnston & Livezey (1987) ed è esplicitata nella documentazione CPC sulle modalità di calcolo degli indici teleconnettivi. In termini fisici, usare un indice definito su Z500 è particolarmente sensato per la Groenlandia perché la fusione estiva risponde in modo marcato alla presenza/posizione delle creste anticicloniche in media troposfera, che modulano getto, avvezioni meridionali e stabilità della colonna, con ricadute dirette su temperatura, umidità e bilancio radiativo alla superficie. Parallelamente, il GBI viene trattato come un indicatore più “locale” e direttamente mirato alla Groenlandia: è definito come la media del geopotenziale a 500 hPa su un riquadro centrato sulla calotta (60–80°N, 20–80°W), concetto introdotto da Fang (2004) e poi ampiamente utilizzato nella letteratura recente sul blocking groenlandese (ad esempio da Hanna e collaboratori). In pratica, mentre il NAOI descrive la modalità dominante della variabilità nord-atlantica su scala più ampia (e quindi il “contesto” emisferico), il GBI quantifica in modo più diretto l’intensità del ridging/blocco sopra la Groenlandia, cioè una delle configurazioni più efficienti nel creare estati predisposte a fusione diffusa e persistente. Coerentemente con questa definizione, nello studio il GBI viene calcolato a partire dalle rianalisi NCEP-CFSR e NCEP-CFSv2, in modo da disporre di campi consistenti di Z500 e quindi di una misura omogenea del blocking nel periodo analizzato.
2.7. Stima degli errori
Gli autori valutano che l’incertezza principale sul numero di giorni di fusione (NMD) derivi da tre famiglie di fattori: (i) l’accuratezza delle temperature impiegate per diagnosticare la fusione (MODIS IST e, quando IST manca, ERA-Interim SKT), (ii) la scelta della soglia termica con cui si decide “melt/no-melt”, e (iii) la qualità della maschera ghiaccio–terra utilizzata per selezionare i soli punti della calotta. La soglia adottata per l’IST è −1,0 °C (coerente con l’impostazione classica dei prodotti MODIS per la Groenlandia), scelta che non coincide con 0 °C perché in ambiente nivale/glaciale la fusione può iniziare e/o persistere nel sottosuolo (firn/snowpack) anche quando la temperatura radiativa della pelle superficiale resta leggermente sottozero: è un effetto legato alla penetrazione della radiazione a onde corte e alla conseguente produzione di calore negli strati sub-superficiali, oltre ai processi di percolazione e rifreezing che rilasciano calore latente e modificano la termica del manto. In quest’ottica, la soglia a −1 °C è un compromesso: riduce il rischio di “perdere” giorni in cui è plausibile melt interno, ma resta abbastanza selettiva da non trasformare in melt un semplice riscaldamento radiativo superficiale senza reale presenza di acqua liquida.
Per quantificare quanto questa scelta conti davvero sul risultato finale, viene eseguita un’analisi di sensibilità ricalcolando l’NMD (per marzo–agosto 2000) con soglie alternative −0,5 °C e −1,5 °C (cioè ±0,5 °C rispetto alla soglia di riferimento). L’effetto emerge soprattutto ai margini della calotta, dove le temperature estive oscillano più spesso attorno alla soglia: nei punti effettivamente sensibili al cambio di soglia, una soglia più permissiva (−1,5 °C) tende ad aumentare l’NMD e una più restrittiva (−0,5 °C) a ridurlo, con un impatto medio dell’ordine di 1 giorno in marzo–maggio e 2 giorni in giugno–agosto. Se però si considera l’intera GrIS (quindi includendo vaste aree interne sempre ben sotto soglia), l’effetto medio “areale” si comprime drasticamente: circa 0,01 giorni in marzo–maggio e 0,5 giorniin giugno–agosto, a conferma che l’incertezza di soglia è un problema soprattutto costiero/marginale, cioè proprio dove la fusione è climatologicamente più frequente e dove i gradienti termici spaziali sono più ripidi.
Resta poi il tema dell’integrazione con ERA-Interim SKT nei giorni in cui IST non è disponibile per nuvolosità. Da un lato, la rianalisi garantisce copertura completa perché la SKT deriva dal bilancio energetico superficiale modellato; dall’altro, introduce un’incertezza legata a errori nelle componenti del SEB (radiazione, turbolenza, conduzione nella neve/firn, rappresentazione delle nubi e dell’umidità), che in ambiente polare possono essere non trascurabili. Gli autori notano che la stagionalità del bias SKT–IST è difficile da correggere in modo esplicito, ma sottolineano anche che questa incertezza “entra in gioco” soprattutto sotto cieli nuvolosi a inizio stagione di fusione, quando il melt è ancora limitato e quindi pesa meno sul conteggio stagionale; per questo motivo la considerano una sorgente d’errore secondaria per l’NMD marzo–agosto e, ancora di più, per l’NMD estivo (JJA), che è la metrica centrale dello studio.
Il terzo pilastro dell’errore riguarda la maschera ghiaccio–terra del prodotto MODIS, a 1,56 × 1,56 km, che viene usata per definire quali pixel appartengano alla calotta. La logica dell’errore è semplice ma importante: se alcuni punti di griglia vengono classificati erroneamente come “ghiaccio” quando in realtà sono “terra”, l’NMD può venire sovrastimato, perché fuori calotta la neve stagionale tende a scomparire in estate lasciando suolo nudo capace di scaldarsi molto più del ghiaccio, innalzando la temperatura superficiale e falsando la diagnosi di melt. La buona risoluzione della maschera riduce questo rischio, ma non lo annulla del tutto nelle zone marginali più complesse (fjord, nunatak, aree miste). Un controllo empirico sull’estate 2000 mostra infatti che in una piccola quota di punti “di calotta” (tipicamente 2–6% dei pixel) l’IST risulta sopra 0 °C senza superare circa 2,9 °C; gli autori interpretano questi casi come aree dove la superficie glaciale è localmente influenzata da rocce e detriti (o affioramenti rocciosi/nunatak prossimi), che assorbono più radiazione solare e possono far salire la temperatura radiativa oltre lo zero pur in un contesto glaciale. Un test di sensibilità sul luglio 2000 evidenzia che l’effetto di questi punti si concentra soprattutto sul margine orientale, proprio dove l’NMD è già elevato: quindi, anche quando c’è un impatto assoluto, l’errore relativo tende a rimanere piccolo e difficilmente altera in modo sostanziale la variabilità interannuale. (È un punto coerente anche con il fatto che, in generale, le metriche di “melt days” basate su soglie sono più vulnerabili alle ambiguità di superficie nelle fasce marginali, e proprio per questo vengono spesso validate o discusse insieme a osservazioni in situ e ad analisi di coerenza del prodotto.)
Tirando le fila, la conclusione operativa è che le fonti dominanti di incertezza sull’NMD, per come è costruito qui, restano l’accuratezza di IST/SKT e la soglia di classificazione, mentre gli effetti della maschera ghiaccio–terra e delle superfici “miste” incidono soprattutto localmente ai margini. In termini quantitativi, gli autori stimano un effetto medio dell’ordine di ±2 giorni in giugno–agosto nelle aree costiere: una forchetta che va letta come incertezza tipica sul conteggio locale dei melt days, ma che, proprio perché relativamente contenuta e confinata soprattutto alle fasce marginali, è ritenuta poco capace di distorcere il segnale delle variazioni interannuali su scala di calotta o su grandi regioni.
Serie temporali delle variabili di fusione in Groenlandia (2000–2014): cosa misurano davvero NMD, estensione di fusione e “melt amount” modellato
Nel periodo 2000–2014 i principali descrittori della fusione superficiale sulla calotta groenlandese (GrIS) non oscillano all’unisono, perché non “vedono” la stessa cosa: il numero di giorni di fusione (NMD) è una misura di persistenza temporale (quante giornate superano una soglia di fusione), l’estensione cumulativa integra quanto a lungo e quanto estensivamente la fusione interessa la superficie nel corso della stagione, l’estensione massima fotografa invece un picco episodico (spesso di pochi giorni), mentre la quantità di fusione modellata traduce l’energia disponibile in un flusso/volume di melt, quindi è più vicina alla “intensità” fisica del processo rispetto ai soli indicatori binari “melt/non-melt”. Questa distinzione concettuale è cruciale per interpretare perché alcune annate risultino coerenti tra metriche cumulative, ma divergano quando si guarda ai massimi giornalieri.
La scelta di normalizzare le serie al 2012 è metodologicamente sensata perché quell’anno rappresenta un riferimento estremo comune a tutte le metriche: nella prima metà di luglio 2012 una configurazione di blocco/alta pressione persistente e anomala favorì cieli più sereni, avvezione calda e un forte forcing radiativo, culminando in un episodio di fusione quasi generalizzata, documentato da più sensori e da diverse catene algoritmiche. Studi osservativi hanno mostrato che l’evento del 2012 fu associato a un’anomala cresta stagnante e a condizioni atmosferiche eccezionali capaci di spingere la fusione fino alle alte quote normalmente “protette”.
Dentro questo quadro, l’NMD calcolato tra marzo e agosto serve a rendere confrontabile una metrica derivata da temperatura superficiale/condizione di melt con grandezze calcolate sull’intera stagione di fusione. Il ciclo medio stagionale evidenzia che marzo–maggio contribuisce poco all’NMD annuale: non è un dettaglio banale, perché implica che gran parte dell’informazione interannuale è controllata dalla finestra estiva (soprattutto giugno–luglio), quando la radiazione netta e i feedback albedo-melt diventano dominanti nelle aree marginali e a quote più basse. In pratica, la variabilità tra anni dipende molto più da quanto a lungo e dove si mantiene la fusione estiva che non da un’eventuale “anticipazione” primaverile.
Quando si confrontano NMD medio su GrIS ed estensione cumulativa della fusione, emerge una coerenza attesa: entrambe sono metriche cumulative e tendono a premiare annate con fusione persistente lungo i margini e nelle fasce a bassa elevazione. Non sorprende quindi che 2012, 2010 e 2002 risultino elevate e che 2000 e 2009 risultino depresse anche nella stima da microonde (con il caveat della disponibilità dati). La microonde passiva, infatti, è storicamente molto adatta a mappare la presenza di acqua liquida nella neve/firn attraverso variazioni della brillanza, ed è stata alla base di molte ricostruzioni dell’estensione di fusione in Groenlandia già dagli anni ’90–2000, con algoritmi come XPGR e successive evoluzioni.
La divergenza più istruttiva arriva però quando si considera l’estensione massima di fusione. Per definizione, un massimo può essere guidato da un singolo episodio sinottico “perfetto” (ad esempio un breve ma intenso blocco con avvezione calda e cieli sereni) anche se la stagione complessiva non è eccezionalmente persistente. È per questo che alcune annate possono mostrare massimi relativamente piccoli ma NMD/estensione cumulativa medio-alti: significa che la fusione è rimasta confinata (spesso ai margini) ma tenace per settimane. Il 2010 è emblematico di questo comportamento: una fusione non necessariamente “esplosiva” in termini di picco spaziale, ma capace di accumularsi nel tempo. Al contrario, 2002 e 2005 (oltre al 2012) superano soglie areali molto ampie e si collocano tra gli anni con massimi più estesi, coerentemente con quanto discusso anche nella letteratura sul ruolo di eventi estivi anomali.
Questa distinzione tra episodico e stagionale non è soltanto “statistica”: riflette processi fisici diversi. Un massimo areale molto grande richiede che per alcuni giorni condizioni prossime alla fusione si estendano verso l’interno e alle alte quote (talvolta fino a regioni come Summit), cosa che dipende in modo sensibile dalla circolazione atmosferica (blocchi groenlandesi, ondulazioni del getto, pattern NAO estivo) e dalla struttura termodinamica in colonna, inclusa la presenza di nubi e il bilancio radiativo. In particolare, la letteratura sul 2012 collega l’evento a configurazioni di blocco e a condizioni capaci di sostenere fusione estesa e, in alcuni casi, prolungata.
Passando al confronto con la fusione modellata da HIRHAM5, la somiglianza dell’andamento interannuale con NMD ed estensione cumulativa indica che il modello intercetta bene il “segnale” di stagione (anni più o meno favorevoli alla fusione persistente). Tuttavia, il fatto che nelle fasi iniziali (2000–2006) la fusione modellata risulti più bassa, mentre dal 2007 in poi diventi comparabile o maggiore, è coerente con due famiglie di spiegazioni che spesso coesistono: (i) questioni di forcing/assimilazione e condizioni al contorno (che modulano l’energia disponibile e la frequenza dei regimi sinottici pro-melt), e (ii) rappresentazione dei processi di superficie e nel firn (albedo, metamorfismo della neve, ritenzione/ricongelamento, drenaggio), che determinano quanta parte dell’energia si traduce in “melt amount” e quanta acqua liquida rimane intrappolata o ricongela. In generale, gli studi su HIRHAM5 e su altri RCM per la Groenlandia mostrano che la componente di fusione/runoff è quella che ha intensificato maggiormente il declino della surface mass balance negli ultimi decenni, con forte sensibilità ai processi di superficie e alle anomalie di circolazione estiva.
L’anomalia del 2002 (melt elevata in NMD/estensione cumulativa ma relativamente bassa nel “melt amount” modellato, una volta normalizzato al 2012) è particolarmente diagnostica: suggerisce una stagione in cui la presenzadella fusione (giorni/mappa) è stata diffusa o persistente, ma l’intensità energetica o la conversione in quantità di fusione/ruscellamento è stata inferiore rispetto a quanto accadrebbe in annate “più radiativamente efficienti” o con albedo più degradato (ad esempio per impurità, metamorfismo e ampliamento delle aree di ghiaccio scoperto). In altre parole, due stagioni possono avere simili “giorni di fusione” ma produrre quantità diverse di melt se cambiano radiazione netta, copertura nuvolosa, temperatura dell’aria, vento (turbulenza) e condizioni iniziali del manto nevoso/firn. Questo è anche uno dei motivi per cui la comunità sta investendo molto nel confrontare e validare algoritmi di rilevamento melt (microonde, L-band, combinazioni dinamiche) e nel collegare meglio le osservazioni di “melt extent” con variabili fisiche come contenuto d’acqua liquida e runoff.
Infine, il punto più “climatico” del passaggio: la tendenza crescente e statisticamente significativa nella quantità di fusione modellata lungo il periodo di studio, a differenza delle metriche puramente descrittive/areali che oscillano fortemente di anno in anno. Questo tipo di risultato è coerente con l’idea che la variabilità interannuale (dominata da regimi sinottici e blocchi estivi) si sovrapponga a un forzante di fondo che aumenta la probabilità e l’efficienza della fusione e del runoff: più energia disponibile, maggiore espansione delle facies di ablazione e maggiore predisposizione del firn a saturarsi e a trasformare fusione in deflusso. Studi recenti mostrano, ad esempio, segnali di aumento del runoff anche da aree di firn che storicamente funzionavano da “spugna” (ritenzione/ricongelamento), con implicazioni dirette per la sensibilità idrologica della calotta.
La Figura 3 è costruita per farci vedere, nello stesso colpo d’occhio, come cambiano da un anno all’altro (2000–2014) quattro “facce” diverse della fusione sulla calotta groenlandese (GrIS) e, in parallelo, quando durante la stagione (marzo–agosto) questi segnali si concentrano. Il punto chiave è che le quattro curve non misurano la stessa cosa: alcune descrivono soprattutto persistenza e accumulo stagionale, altre catturano picchi episodici, altre ancora provano a tradurre la fusione in una quantità fisica (melt amount) coerente con un bilancio energetico.
Nel pannello (a) tutte le serie sono normalizzate al 2012 (cioè il 2012 vale 1 per ciascuna metrica), così i confronti diventano “relativi”: quanto vale un certo anno rispetto all’anno estremo. Il 2012 è un riferimento naturale perché fu caratterizzato da un evento di fusione eccezionale: osservazioni satellitari indicarono fusione “a” o “vicino” alla superficie su circa 98.6% dell’estensione della calotta il 12 luglio 2012, associata a una cresta di aria calda anomala e stagnante sulla Groenlandia. Questa eccezionalità emerge nel grafico perché tutte le metriche (giorni di fusione, estensione cumulativa, massimo areale, fusione modellata) convergono sul picco 2012: è il segnale tipico di un anno in cui l’anomalia è stata sia spazialmente vasta sia temporaneamente persistente.
Le quattro curve del pannello (a) derivano anche da filiere osservative/modellistiche diverse. Il numero medio di giorni di fusione (NMD) qui è legato a prodotti di temperatura superficiale (MODIS IST, ice-surface temperature) e alla skin temperature di una rianalisi (ERA-Interim). L’IST MODIS per la Groenlandia è stato sviluppato proprio come record “climate-quality”, con controlli di qualità e omogeneità temporale, per analisi di variabilità e trend. L’estensione di fusione (cumulativa e massima) proviene invece da microonde passive (dataset in ambito MEaSUREs/NSIDC): è una famiglia di prodotti molto usata perché la microonda è sensibile alla presenza di acqua liquida nella neve/firn e non soffre la copertura nuvolosa come l’ottico/IR. La logica di rilevamento melt in microonde (ad esempio tramite XPGR e soglie collegate alla wet snow) è un caposaldo storico per Groenlandia e mostra anche un risultato “climatologico” coerente con il pannello (b): la fusione tende a piccare a fine luglio. Infine, la curva di melt amount modellato(HIRHAM5) non è un “sì/no” di fusione, ma una quantità legata alla fisica del modello regionale e ai parametri di superficie/firn; HIRHAM5 è ampiamente usato per la surface mass balance groenlandese e gli studi sottolineano che il declino recente della SMB è guidato soprattutto da melt e runoff in aumento.
Una volta chiarita la natura delle metriche, diventa più leggibile perché alcune annate si assomiglino e altre no. In generale, NMD medio ed estensione cumulativa tendono a covariare perché sono entrambi indicatori “integrati” sull’intera stagione: se la fusione resta attiva a lungo (anche solo sui margini a bassa quota) accumuli molti giorni di melt e, sommando giorno dopo giorno, accumuli anche molta estensione cumulativa. Questo spiega perché anni come 2010 e 2002 risultino elevati in entrambe. La massima estensione, invece, è spesso più “capricciosa”: è sensibile a un singolo episodio sinottico favorevole (un’ondata calda breve ma molto vasta) e quindi può risultare relativamente bassa in anni in cui la fusione è stata sì persistente, ma confinata alle fasce marginali. È esattamente il tipo di situazione che gli autori richiamano quando osservano anni con NMD/cumulativa alti ma massimo più modesto: significa “estate lunga ai bordi” più che “grande botta areale” nell’interno.
Questa distinzione ha un significato fisico: un massimo areale molto grande richiede che, anche alle quote più elevate e nelle facies più fredde, il bilancio energetico superficiale vada almeno per qualche giorno vicino alla fusione (ruolo combinato di avvezione calda, cieli sereni o nubi “efficaci” radiativamente, e feedback dell’albedo). Nel 2012, ad esempio, oltre alla configurazione di blocco, è stato mostrato un contributo importante delle nubi basse liquide: possono ridurre un po’ la shortwave entrante ma aumentare la longwave verso il basso, e in determinate condizioni il risultato netto è un potenziamento della fusione superficiale su aree molto estese. Al contrario, una stagione con fusione persistente ma periferica può gonfiare NMD e cumulativa senza produrre un massimo “da record”, perché l’interno resta sotto soglia per quasi tutto il tempo.
Il pannello (b) mette ordine sul timing: mostra il numero medio mensile di giorni di fusione (media su tutta la GrIS) da marzo ad agosto, includendo la media 2000–2014 e alcuni anni rappresentativi. Qui si vede chiaramente che primavera (marzo–aprile) pesa pochissimo: la fusione media sulla calotta è quasi nulla, poi compare debolmente a maggio, accelera in giugno e raggiunge il picco in luglio, per poi calare in agosto. Questo profilo è coerente con la climatologia della fusione da microonde (picco a fine luglio) e, in termini energetici, con la combinazione di massima insolazione stagionale, progressiva riduzione dell’albedo (metamorfismo della neve, affioramento di ghiaccio più scuro nelle zone di ablazione) e maggiore efficienza dei flussi turbolenti nelle aree marginali.
La cosa più interessante del pannello (b) non è solo la forma “a campana”, ma dove gli anni estremi si separano dalla media. Il 2012 non è soltanto “più alto in luglio”: mostra un surplus già marcato tra giugno e luglio, cioè nel cuore della stagione in cui la calotta è più sensibile ai feedback (albedo-melt e, in alcuni casi, anche effetti radiativi delle nubi). Il 2013, al contrario, resta sotto la media soprattutto nel periodo di massimo potenziale di fusione, segnalando che la differenza tra anni non dipende tanto da “un po’ di melt primaverile”, quanto dal modo in cui si configura l’estate: circolazione atmosferica, persistenza dei pattern sinottici e condizioni radiative che determinano se la fusione rimane episodica o diventa strutturale su settimane.
In sintesi, la figura ci sta dicendo tre cose operative: (1) 2012 è un outlier robusto perché è estremo per tutte le metriche (spazio, tempo e intensità). (2) Se vuoi descrivere “quanto è stata lunga la stagione di melt”, NMD e cumulativa sono i descrittori più coerenti; se vuoi cogliere “l’evento shock” areale, la massima estensione è la metrica più sensibile ma anche più episodica. (3) La metrica modellata (HIRHAM5) aggiunge una dimensione fisica (quantità di fusione) e può mostrare trend o scarti rispetto agli indicatori binari perché dipende dalla rappresentazione di albedo, firn e ritenzione/ricongelamento; ed è in linea con la letteratura che attribuisce al melt/runoff un ruolo centrale nel declino recente della surface mass balance groenlandese.
Stagioni di fusione e condizioni atmosferiche (2000–2014): perché i descrittori di melt possono “raccontare storie diverse”
Per interpretare le discrepanze tra i descrittori di fusione mostrati in Figura 3a (numero medio di giorni di fusione, estensione cumulativa, estensione massima e melt amount modellato), il passaggio chiave è ragionare in termini di scala temporale e meccanismo fisico: alcune metriche rispondono soprattutto alla persistenza stagionale (quante giornate sopra soglia, per quanto tempo rimane attiva la fusione ai margini), altre sono dominate da episodi sinottici brevi ma intensi (picchi areali), mentre la fusione modellata tende a pesare maggiormente l’energia disponibile nelle zone a bassa quota dove si concentra la maggior parte del melt e dove la stagione è più lunga. In questo quadro, indici come NAOI e Greenland Blocking Index (GBI) sono strumenti diagnostici utili perché riassumono lo “stato” della circolazione estiva sull’Atlantico settentrionale e sulla Groenlandia: estati con blocco più marcato su Groenlandia (GBI alto) sono spesso associate a pattern che favoriscono advezione calda e/o cieli più sereni e quindi condizioni più favorevoli alla fusione, anche se il dettaglio dipende da posizione e persistenza del blocco e dal ruolo delle nubi.
L’estate 2002 è un esempio didattico di come un anno possa produrre un segnale molto forte nell’estensione massima(secondo valore più alto, ~86% del record 2012) pur mostrando, nella metrica modellata, un’intensità complessiva relativamente più contenuta (circa metà del 2012). In quel caso la fusione estesa è descritta come episodica, con un picco concentrato in pochi giorni di inizio luglio: è esattamente la firma di una dinamica dominata da intrusioni di aria calda da sud che spingono temporaneamente la fusione in regioni (come Groenlandia nord-centrale e nord-orientale) dove normalmente non è frequente su vasta scala. In altre parole, un singolo “colpo” sinottico ben posizionato può gonfiare l’estensione massima senza che ciò implichi, automaticamente, una stagione uniformemente estrema in termini di energia integrata o produzione totale di melt.
Sul piano delle condizioni medie JJA, l’interpretazione proposta (pressione anomala più bassa su gran parte della calotta) è coerente con uno schema che favorisce flusso più meridionale e una riduzione della componente occidentale su Groenlandia centrale e orientale: ciò si riflette in anomalia positiva della componente meridiana del vento e negativa della zonale, con conseguente incremento delle temperature in quota fino a circa 500 hPa sul settore nord-orientale, dove infatti l’anomalia positiva di NMD risulta più marcata. È importante anche il punto “negativo”: il fatto che le anomalie di velocità verticale a 500 hPa non si allineino chiaramente al pattern di NMD suggerisce che, su scala stagionale, il contributo di riscaldamento/raffreddamento adiabatico non sia stato il controllore principale del melt, mentre risultano più determinanti l’advezione e il bilancio radiativo/turbolento al suolo.
La discrepanza tra segnali elevati in NMD/estensione cumulativa e un melt amount modellato relativamente ridottonel 2002 si capisce se si considera dove “pesa” la fusione totale: i modelli di SMB (inclusi RCM come HIRHAM5) sono fortemente sensibili alle condizioni estive medie nelle aree di ablazione a bassa quota, soprattutto nel settore sud-occidentale, perché lì la stagione è più lunga e la produzione di melt può trasformarsi più efficacemente in runoff. Se in quelle regioni prevalgono anomalie termiche negative o condizioni meno favorevoli, la quantità totale modellata può rimanere moderata anche in presenza di un episodio areale eccezionale altrove.
Nel 2008 (per come è descritto nel brano) emergono due ingredienti interessanti. Il primo è il possibile contributo di un’anomalia positiva di u10 (vento a 10 m) alla fusione nel nord-est e in una stretta fascia del sud-est, dove le anomalie di vento e NMD risultano co-localizzate: questo richiama meccanismi in cui il vento amplifica la fusione attraverso flussi turbolenti (sensibile e latente), rimescolamento dello strato limite e/o episodi di venti discendenti caldi e secchi (foehn) che, nelle regioni orientali della Groenlandia, sono oggi riconosciuti come driver efficaci di eventi di melt intenso e rapido. Il secondo è l’idea che una stagione di fusione più lunga lungo la costa nord-orientale sia associata a una lieve anomalia positiva di TCW: maggiore disponibilità di vapore acqueo e umidità colonnare può sostenere un forcing radiativo aggiuntivo (specie via longwave in presenza di nubi basse liquide) e/o segnalare fasi di trasporto umido (in alcuni casi legate ad atmospheric rivers) che modulano fortemente il bilancio energetico superficiale. Allo stesso tempo, il fatto che nel sud-ovest compaiano anomalie negative di NMD spiega perché il massimo NMD annuale resti vicino alla media: è un promemoria che “il cuore” della produzione di melt totale non coincide sempre con le aree dove si osservano anomalie positive più appariscenti in altre regioni.
L’estate 2010 (che nella Figura 3a risulta alta per NMD/estensione cumulativa ma con estensione massima relativamente più contenuta rispetto agli anni record) è coerente con uno scenario di fusione persistente ma più periferica: tante giornate di melt e accumulo stagionale lungo le fasce marginali, senza necessariamente un singolo episodio capace di espandere la fusione su porzioni enormi dell’interno. Questo tipo di comportamento si colloca bene nella letteratura che collega la variabilità della fusione groenlandese alla frequenza/intensità degli episodi di blocco e alle anomalie di geopotenziale su Groenlandia, più che a un singolo indice “semplice” preso isolatamente.
Il 2012 rappresenta invece il caso in cui le diverse metriche si “allineano” perché l’anomalia è stata contemporaneamente vasta, persistente e fisicamente efficiente: record di estensione di fusione nella serie satellitare, durata eccezionale della stagione, e un contesto di circolazione anomala con blocco e intrusioni calde che hanno favorito condizioni radiative estreme. Qui è fondamentale non cadere in una lettura monocausale “solo cieli sereni”: per l’evento 2012 è stato mostrato anche un ruolo chiave delle nubi basse liquide nell’aumentare la radiazione infrarossa verso il basso e quindi la temperatura prossima alla superficie in fasi critiche, mentre su scale più ampie la copertura nuvolosa può tanto amplificare quanto attenuare il melt a seconda del regime radiativo e della stagione.
Il 2013, infine, è il contraltare: la Figura 3 lo colloca tra le estati a fusione più contenuta nel periodo, e osservazioni indipendenti mostrano un massimo di area in fusione molto più basso rispetto al 2012 (ordine del 44% contro ~97% nell’estate precedente). Sul piano dinamico-radiativo, diversi lavori interpretano la fase post-2013 come caratterizzata più spesso da condizioni coerenti con NAO estiva positiva, maggiore nuvolosità totale e riduzione della shortwave al suolo, fattori che contribuiscono a “raffreddare” il bilancio energetico superficiale e quindi a ridurre la fusione.
Mettendo insieme questi casi, la lezione della figura (e della discussione associata) è che la fusione groenlandese va letta come risultante di un equilibrio tra forzanti sinottiche episodiche (che dominano l’estensione massima), persistenza stagionale e localizzazione altimetrica (che dominano NMD ed estensione cumulativa), e bilancio energetico nelle aree di ablazione (che domina la quantità modellata e, in ultima analisi, il contributo al runoff). In pratica, un anno può “sembrare enorme” se guardi solo al picco areale, ma non esserlo altrettanto in termini di melt totale; oppure può non mostrare un picco spettacolare e tuttavia produrre molta fusione perché mantiene condizioni favorevoli per settimane lungo i margini. Nel triennio 2010–2013 si vede benissimo perché, in Figura 3a, i descrittori di fusione possano divergere: il sistema può produrre (i) stagioni lunghe e “marginali” (molti giorni di fusione lungo le coste e a bassa quota) senza necessariamente generare un picco areale massimo molto elevato, oppure (ii) episodi sinottici rari capaci di spingere il melt su porzioni vaste dell’interno, facendo esplodere l’estensione massima e rendendo coerenti tutte le metriche. In questo quadro, gli indici di circolazione (NAO/GBI), i campi di Z500 e le anomalie di umidità colonnare (TCW)aiutano a collegare la dinamica atmosferica alla risposta criosferica, perché riassumono quanto spesso e quanto efficacemente arrivino sulla calotta masse d’aria calde/umide e quanto persistano configurazioni di blocking su Groenlandia.
Nel 2010 la fusione risulta sopra media in quasi tutte le metriche, tranne l’estensione massima, ed è un segnale molto coerente con una stagione che inizia presto, dura a lungo e presenta un picco netto in luglio, ma senza “colpi” di avvezione calda abbastanza forti da portare la fusione in modo esteso verso le parti centrali della calotta (da qui l’assenza di grandi picchi giornalieri nell’area in fusione). Il punto fisico è che puoi avere un NMD elevato perché la fusione è persistente lungo i margini (dove la stagione è naturalmente più lunga), e al contempo avere un massimo areale relativamente contenuto se l’interno rimane spesso sotto soglia. Il risultato di Mernild e colleghi, che distingue una durata di fusione molto più lunga nel sud-ovest e più breve nel nord-est nel 2010, va esattamente in questa direzione: è la geografia (quote basse vs alte) a determinare quanta fusione “si accumula” nel tempo.
Sempre nel 2010, l’assetto atmosferico descritto (temperature troposferiche positive fino a 500 hPa, anomalie di TCW positive quasi ovunque eccetto il nord-est, NAO negativa east-based e GBI leggermente sopra media) suggerisce una stagione in cui il contributo dominante alla fusione lungo le coste sia arrivato da maggiori apporti di calore/umidità e dalla persistenza stagionale del melt, più che da un singolo episodio “da record” nell’interno. Qui si innesta anche un meccanismo cruciale e molto documentato: il feedback dell’albedo. Se l’innevamento è relativamente ridotto e il ghiaccio nudo affiora prima, la superficie scurisce, assorbe più radiazione e “auto-rafforza” la fusione, amplificando la risposta soprattutto nelle aree di ablazione del versante occidentale, che sono anche quelle che pesano di più sulla quantità totale di melt/runoff.
Il 2012 è l’estremo opposto (e infatti in Figura 3a tutte le curve vanno insieme): non serve un inizio stagione eccezionalmente anticipato per fare un anno record; basta che, già da giugno, si instaurino condizioni capaci di imporre riscaldamento diffuso dell’intera colonna troposferica e anomalie positive di TCW su larga scala, mantenendo poi valori estremi anche in agosto. La letteratura sull’evento 2012 converge su una catena causale robusta: una circolazione anomala con blocking/alte altezze geopotenziali su Groenlandia (GBI alto, NAO fortemente negativa west-based) ha favorito un trasporto eccezionalmente efficace di calore e umidità e un riscaldamento profondo, in alcuni casi legato a un Atmospheric River che convoglia aria calda (anche continentale) dal Nord America verso la Groenlandia e poi sopra la calotta. Un dettaglio importante, spesso sottovalutato quando si parla solo di “cieli sereni”, è che durante la fase critica del melt 2012 le nubi basse liquide possono aumentare la radiazione infrarossa verso il basso e spingere la superficie oltre soglia anche in condizioni non necessariamente ottimali per la shortwave: Bennartz et al. mostrano proprio il ruolo radiativo di queste nubi nel potenziare l’estensione della fusione nel luglio 2012.
Il 2013, infine, si comporta come un “anti-2012”: tutte le metriche si avvicinano ai minimi, NMD è sotto media quasi ovunque e il quadro di circolazione medio stagionale torna coerente con una NAO positiva e un GBI molto basso, cioè una configurazione meno favorevole al trasporto da sud-ovest e più incline a mantenere masse d’aria fredde sul dominio groenlandese. Häkkinen et al. collegano esplicitamente la bassa fusione del 2013 a una NAO positiva, mentre la reportistica climatica artica di quell’anno documenta temperature estive più basse su Groenlandia e una conseguente riduzione di estensione/durata del melt. In pratica: nel 2013 non basta avere localmente segnali di subsidenza o ω “favorevoli” in qualche settore, perché la disponibilità di aria calda/umida e il bilancio termico della colonna sono complessivamente sfavorevoli; e quando TCW e temperature troposferiche risultano negative su vasta scala, la fusione “stagionale” si spegne.
In questo passaggio si capisce in modo molto “pulito” che la differenza fra 2012 e 2013 non è solo una questione di “quanto” ha fuso, ma di quali traiettorie atmosferiche hanno rifornito (o privato) la calotta di calore e umidità, soprattutto nelle aree interne dove il sistema è spesso vicino alla soglia di fusione. Nel 2013 le anomalie delle componenti del vento a 10 m risultano in generale deboli; il segnale più incisivo, coerente con le anomalie negative di NMD, è l’indebolimento dei venti meridionali nella metà meridionale della Groenlandia centrale, cioè una riduzione del contributo advettivo “mite” verso l’interno. Il confronto fra i pattern spaziali di NMD del 2012 e del 2013 mette in evidenza che la regione con la maggiore divergenza è l’alto plateau del nord-est: un’area che, per quota e climatologia, tende a essere estremamente sensibile a piccoli cambiamenti nella massa d’aria in ingresso (temperatura, contenuto di vapore, nubi e bilancio radiativo). Nel lavoro a cui stai facendo riferimento, le traiettorie retrospettive in quel punto (77°N, 30°W) mostrano che nel 2013 diminuisce la quota di masse d’aria di origine sud-occidentale rispetto al 2012, mentre aumentano le provenienze sud-orientali e nord-orientali; e nei casi nord-orientali la temperatura media a Summit risulta addirittura di ~7 °C più bassa nel 2013 rispetto al 2012. Questo è un dettaglio enorme dal punto di vista fisico: sul plateau, dove spesso “manca poco” per raggiungere la fusione, un raffreddamento di quella entità nella massa d’aria associata a determinate traiettorie può annullare il melt anche se altri indicatori dinamici (ad esempio segnali locali di ω) suggerirebbero condizioni favorevoli. È coerente con la letteratura che mette al centro, per la Groenlandia, il ruolo del trasporto di umidità e calore (TCW e temperatura in colonna) e dei corridoi di flusso legati a blocking/NAO, inclusi i casi in cui gli ingressi umidi e caldi assumono la forma di eventi tipo “atmospheric river” capaci di modificare rapidamente il bilancio energetico superficiale.
C’è poi un punto molto istruttivo, quasi da manuale di circolazione artica: nel 2013 lo stesso assetto su larga scala che limita il melt groenlandese può contemporaneamente favorire condizioni opposte su altri settori dell’Artico atlantico. Il fatto che venga citato un melt record alle Svalbard associato a venti sud–sudoccidentali proprio nel 2013 è un esempio perfetto di “redistribuzione” dei flussi: non è che il sistema emisferico diventi ovunque più freddo o più caldo, ma cambia la geometria delle traiettorie e dei canali di trasporto, penalizzando la Groenlandia e alimentando invece l’advezione mite verso Svalbard. Lang et al. documentano per il 2013 un cambio di circolazione verso sud–sudovest su Svalbard e un melt/SMB record nella serie 1979–2013.
Quando il testo passa a riassumere le altre estati, la trama comune è sempre la stessa: il “background” stagionale(TCW, temperature troposferiche, regime di Z500/GBI/NAO) decide se la stagione è predisposta a una fusione diffusa e persistente, mentre gli eventi sinottici determinano se compariranno “spike” tali da gonfiare l’estensione massima. L’estate 2000 incarna la dicotomia: anomalie negative in tutte le variabili di melt per condizioni complessivamente fredde e secche, ma con un’estensione massima non così bassa perché due eventi sinottici riescono comunque a produrre episodi relativamente estesi. Questo tipo di contrasto tra cumulato basso e massimo relativamente alto non è un dettaglio del singolo dataset: è un comportamento che emerge anche quando l’estensione di melt viene ricostruita da microonde passive (XPGR e affini) e confrontata con simulazioni RCM (MAR), cioè lungo filiere indipendenti.
Nel 2001 la combinazione è quasi “speculare”: tutte le variabili hanno anomalie negative, ma con valori molto piccoli sia per melt amount modellato sia per estensione massima e, al tempo stesso, NMD ed estensione cumulativa relativamente più alte. La lettura fisica è coerente con quanto discusso dagli autori: un’estate mediamente fredda e senza episodi sinottici capaci di espandere la fusione sull’altopiano (quindi niente picchi nell’estensione giornaliera), ma con una certa persistenza di melt marginale sufficiente a far crescere le metriche cumulative rispetto a quanto ci si aspetterebbe guardando solo il massimo.
Il 2003 è un caso didattico per capire perché maximum melt extent e modeled melt amount possano “litigare”: l’estensione massima resta tra le più basse, mentre il melt amount modellato risulta relativamente alto (addirittura superiore al 2002, una volta normalizzato al 2012), verosimilmente sostenuto da TCW un po’ più elevato e da condizioni favorevoli nelle aree a bassa quota dove si concentra la produzione di melt che domina la metrica modellata. In altre parole: puoi avere una stagione energeticamente efficace lungo i margini senza avere il singolo evento sinottico capace di portare melt esteso verso l’interno. Questo tipo di separazione tra controlli “stagionali” e controlli “event-scale” è esattamente uno dei messaggi centrali dell’analisi 2000–2014.
Nel 2004 l’NMD positivo nella metà meridionale viene ricondotto a un assetto di Z500 (anomalia negativa su Baffin Bay) che favorisce venti occidentali/sudoccidentali al suolo; eppure il melt amount modellato resta basso. Anche qui la lezione è che la circolazione favorevole non basta: serve che il bilancio energetico superficiale (radiazione, flussi turbolenti, albedo, stato del firn) converta quell’assetto in fusione “integrata” nelle zone chiave di ablazione. Nel 2005la vicinanza alla media nasce da compensazioni spaziali (anomalie negative al sud e centro, positive lungo un’ampia fascia marginale), con un ruolo plausibile delle anomalie positive di TCW nel sostenere la fusione lungo i bordi.
Per 2006–2007 le anomalie restano deboli; ma il 2007 è particolarmente interessante perché TCW e GBI risultano elevati e “intuitivamente” ci si aspetterebbe più fusione. Häkkinen et al. spiegano però un punto molto robusto: anche con molti giorni di blocking, se le anomalie di temperatura atmosferica sono troppo piccole, non si innescano advezioni calde sufficienti a produrre estremi di melt. È un promemoria: il blocking è una condizione di contorno favorevole, ma la risposta dipende dalla termodinamica della massa d’aria e dalla sua traiettoria effettiva verso la Groenlandia.
Il confronto 2009 vs 2011 mette in guardia dall’usare NAO/GBI come “interruttori” deterministici: Lim et al. discutono come, in entrambi gli anni, un blocking locale (più debole che nel 2010 e 2012) possa aver inibito il riscaldamento della GrIS nonostante una NAO molto negativa; tuttavia, nelle metriche analizzate, il 2009 resta limitato (anche perché il campo medio Z500 favorisce componenti più orientali, meno efficaci nel trasporto sud-occidentale), mentre il 2011 emerge tra le estati più intense del periodo, con segnali più favorevoli a nord legati a Z500 positivo (circolazione anticiclonica e trasporto di aria più mite), ma con freni regionali lungo costa est e sud per venti settentrionali più forti. È il classico caso in cui lo stesso “segno” dell’indice grande scala può produrre risposte diverse a seconda della posizione del massimo di geopotenziale, dei gradienti e dei percorsi delle masse d’aria.
Infine, il 2014 mostra in modo quasi “chirurgico” la fragilità concettuale della sola estensione massima: NMD e melt amount modellato sono vicini alla media, ma l’estensione massima è la più bassa dell’intera serie perché il massimo stagionale viene raggiunto molto presto (17 giugno) e poi la stagione procede solo con eventi sinottici di estensione leggermente inferiore. Qui la metrica del massimo non descrive “una stagione scarsa” in senso generale: descrive l’assenza del singolo giorno-evento capace di espandere la fusione su un’area più ampia.
Nel complesso, questa continuazione rafforza il messaggio metodologico della figura e del testo: per capire davvero la fusione groenlandese bisogna tenere insieme (i) la persistente forzante stagionale (TCW e temperatura in colonna, spesso modulati da blocking/NAO e da trasporti umidi organizzati), (ii) la geografia del melt (bassa quota/margini vs plateau interno), e (iii) la componente episodica che governa l’estensione massima. Se manca uno di questi tre piani, il rischio è interpretare male perché un anno appare “grande” in una metrica e “normale” in un’altra.
La Figura 4 rappresenta l’andamento giornaliero dell’estensione della fusione superficiale (melt extent) sulla calotta groenlandese per il periodo 2000–2012, ricavato dal prodotto MEaSUREs “Greenland Surface Melt Daily” basato su osservazioni a microonde passive (griglia EASE-Grid 2.0). In termini fisici, questa variabile non misura “quanta” acqua di fusione viene prodotta, ma quanta superficie della calotta presenta condizioni compatibili con melt (o presenza di acqua liquida nella neve/firn) in ciascun giorno: è dunque un indicatore estremamente sensibile alla geografia e alla sinottica del riscaldamento estivo, e risulta particolarmente utile per confrontare stagioni dominate da persistenza rispetto a stagioni dominate da episodi brevi ma intensi.
Nel grafico, la linea nera è la media climatologica giornaliera 2000–2012; l’area grigia rappresenta la forbice interannuale (il range di variabilità) per ciascun giorno dell’anno; le linee blu (2002) e rossa (2012) evidenziano due estati con comportamento “estremo”, ma per ragioni in parte diverse. La climatologia (linea nera) descrive un ciclo stagionale netto: l’estensione di fusione resta prossima allo zero per gran parte dell’anno, poi cresce rapidamente tra tarda primavera e inizio estate, raggiunge il massimo nel cuore dell’estate boreale (intorno ai giorni dell’anno tipici di inizio–metà luglio) e decresce verso fine estate. Questo profilo è coerente con il controllo dominante del bilancio energetico superficiale (radiazione netta, flussi turbolenti) e con il fatto che, in Groenlandia, la fusione più estesa si concentra nella finestra in cui insolazione e temperatura dell’aria rendono più facile superare la soglia di melt anche oltre le sole fasce costiere.
La fascia grigia è forse l’elemento più informativo: mostra quanto il sistema sia “guidato dagli eventi”. A parità di giorno dell’anno, l’estensione di fusione può oscillare enormemente tra anni diversi, segnalando che la fusione groenlandese non è soltanto un fenomeno radiativo stagionale “regolare”, ma risponde in modo amplificato ai regimi di circolazione estiva (blocking, traiettorie delle masse d’aria, trasporto di calore/umidità). Questa variabilità interannuale è ampiamente documentata negli studi che collegano le estati più intense a configurazioni anticicloniche persistenti su Groenlandia (GBI elevato) e/o a NAO estiva negativa, che favoriscono advezioni calde verso la calotta e condizioni radiative/termiche più favorevoli al melt.
Il confronto tra 2002 (blu) e 2012 (rosso) chiarisce bene la differenza tra un’estate con picchi molto acuti e un’estate con persistenza + picchi record. Nel 2002 la curva mostra un’impennata molto marcata in corrispondenza di pochi giorni di luglio: è la firma tipica di episodi sinottici di advezione calda che, per breve tempo, riescono a spingere la fusione su aree molto ampie, includendo settori interni o nord-orientali dove normalmente il melt è raro su grande scala. Questo tipo di dinamica è coerente con interpretazioni che attribuiscono al 2002 alcuni casi “concentrati” di trasporto caldo verso Groenlandia, capaci di gonfiare l’estensione massima senza necessariamente implicare una stagione uniformemente estrema in termini di durata.
Il 2012 si colloca invece in una classe a parte: la curva rossa non solo raggiunge i valori più elevati (nel grafico si arriva a circa 1.6 milioni km², cioè quasi l’intera calotta), ma mostra anche una fase estiva con estensione molto alta non limitata a un singolo spike. Il massimo più spettacolare coincide con l’episodio di metà luglio 2012: osservazioni satellitari indipendenti indicano che il 12 luglio 2012 condizioni di fusione erano presenti “a” o “vicino” alla superficie su ~98.6% dell’estensione della calotta, un evento eccezionale nella serie moderna. La letteratura attribuisce questo comportamento a una combinazione di fattori atmosferici ad ampia scala: una cresta di alta pressione persistente con circolazione anticiclonica su Groenlandia, capace di favorire (a seconda della fase) cieli più sereni e/o traiettorie di aria più mite verso la calotta, oltre a un marcato riscaldamento della colonna troposferica. Un aspetto molto robusto (e spesso sottovalutato) riguarda inoltre il ruolo delle nubi basse liquide: nel luglio 2012 è stato mostrato che nubi ricche di acqua liquida, attraverso l’aumento della radiazione infrarossa verso il basso, possono incrementare la temperatura prossima alla superficie e amplificare l’estensione del melt durante le fasi critiche, contribuendo a spiegare perché l’evento sia stato così diffuso anche in aree non “facili” da far fondere.
Letta in questa chiave, la figura diventa quasi una “radiografia” dei meccanismi: la linea nera ti dice qual è la finestra stagionale in cui la Groenlandia è predisposta alla fusione estesa; la banda grigia quantifica quanto contino le configurazioni sinottiche nel modulare giorno per giorno l’area in melt; le curve 2002 e 2012 mostrano rispettivamente un caso dominato da episodi intensi ma relativamente concentrati e un caso in cui forzante atmosferica + feedback radiativi portano a un regime davvero estremo e più persistente. Questa distinzione è cruciale anche per collegare Figura 4 a Figura 3: un anno può avere un’estensione massima molto alta (per pochi giorni) senza necessariamente presentare lo stesso grado di eccezionalità nelle metriche cumulative, mentre un anno come il 2012 tende ad allineare massimi e cumulati perché l’anomalia è sia vasta sia sostenuta.
La Figura 5 è pensata per ricostruire, in modo “diagnostico”, perché l’estate 2002 abbia prodotto un segnale di fusione rilevante pur mostrando, su scala stagionale, un quadro atmosferico non necessariamente estremo come quello del 2012. Il punto metodologico è che qui si mettono insieme un indicatore di risposta della calotta (NMD, numero di giorni di fusione) e una batteria di variabili atmosferiche di ERA-Interim (Z500, TCW, temperatura a 500 hPa, vento a 10 m, ω a 500 hPa), tutte espresse come anomalie rispetto alla climatologia JJA 2000–2014. In questo modo si distinguono i meccanismi che favoriscono la fusione per persistenza stagionale (anomalie termodinamiche e di trasporto) da quelli che la amplificano per episodi sinottici brevi (intrusioni calde, picchi areali), che possono essere “smussati” dalle medie JJA.
Nel pannello (a) le anomalie di NMD (colori) sono sovrapposte alle anomalie di Z500 (contorni). La distribuzione suggerisce che nel 2002 l’incremento dei giorni di fusione non è uniforme: il segnale positivo è più evidente in settori specifici (in particolare verso il nord-est e lungo porzioni marginali), mentre altre aree mostrano anomalie più deboli o negative. La co-occorrenza con anomalie di Z500 indica che la circolazione in media troposfera ha modulato la probabilità di condizioni pro-melt, ma soprattutto suggerisce un aspetto chiave evidenziato da studi sul melt groenlandese: non basta “un” indice di circolazione, conta dove si colloca l’anomalia di geopotenziale e come essa ridisegna i flussi di massa d’aria verso la calotta, in particolare l’efficacia delle intrusioni calde e la loro capacità di raggiungere settori interni/settentrionali.
Il pannello (b) (NMD totale in JJA 2002) mostra un pattern molto coerente con la climatologia fisica della Groenlandia: i valori più alti si concentrano nelle fasce costiere e a bassa quota, mentre l’altopiano interno resta su valori bassi. Questo è importante perché spiega perché metriche “cumulative” (giorni di fusione, estensione cumulativa) possono risultare elevate anche senza un’estensione massima eccezionalmente persistente nell’interno: la fusione che “fa massa” (e spesso anche runoff) si gioca soprattutto nelle aree marginali, dove la stagione è più lunga e il bilancio energetico è più spesso positivo.
Il pannello (c) (anomalia di TCW) mostra, nel 2002, un segnale complessivamente modesto: alcune aree presentano lievi anomalie positive (più marcate sul lato orientale), altre lievi negative. Questo è un indizio cruciale: l’umidità colonnare è spesso un amplificatore del melt perché aumenta la radiazione infrarossa verso il basso e perché segnala trasporto di masse d’aria potenzialmente più miti. Tuttavia, nel 2002 la “firma” media stagionale di TCW non è così forte, e ciò è perfettamente compatibile con l’interpretazione che il grande melt del 2002 abbia avuto una componente importante episodica, concentrata in finestre di pochi giorni. In effetti, analisi basate su MODIS indicano che una parte consistente della fusione 2002 fu di tipo effimero (1–2 giorni) e concentrata a cavallo tra fine giugno e inizio luglio, associata a configurazioni sinottiche favorevoli.
Il pannello (d) (anomalia di temperatura a 500 hPa) aiuta a capire dove la colonna troposferica fosse più favorevole al melt: nel 2002 emerge un contrasto tra un settore con anomalie più fredde (in prossimità del lato occidentale) e un settore con anomalie più miti verso nord-est. Questo è coerente con un’idea ben documentata: gli anni con maggiore melt non dipendono solo dal blocking in sé, ma dalla combinazione tra pattern di circolazione e anomalia termica associata alle intrusioni; Häkkinen et al. mostrano che la temperatura anomala legata alle intrusioni subtropicali è un ingrediente determinante per spingere la frazione di area in melt verso valori elevati.
I pannelli (e–h) (vento a 10 m) servono a separare due contributi: la componente zonale u10 (e, f) e la componente meridiana v10 (g, h). In media JJA, il campo di u10 mostra la struttura attesa della circolazione vicino-superficie attorno alla calotta, mentre l’anomalia (f) suggerisce nel 2002 un indebolimento/localizzazione dei flussi zonali in alcune aree (riduzione dei westerlies o rafforzamento relativo di componenti easterly). Ancora più informativo è v10: la mappa media (g) e la sua anomalia (h) indicano che, in porzioni importanti della Groenlandia centrale/orientale, nel 2002 si è favorito un contributo più “meridiano” (componente da sud relativamente più efficace o indebolimento dei flussi da nord), coerente con l’idea che la fusione aumenti quando la calotta viene raggiunta più spesso da masse d’aria trasportate lungo traiettorie che favoriscono advezione calda e, talvolta, umida. Questa lettura è coerente anche con la letteratura sul 2002 basata su microonde passive, che descrive quell’estate come un’anomalia di melt significativa nella Groenlandia recente e discute il ruolo delle condizioni atmosferiche nel modulare estensione e durata della fusione.
Il pannello (i) (anomalia di ω a 500 hPa) è quasi “negativo” per design: mostra che nel 2002 non c’è una corrispondenza netta tra le zone con ω anomala (sollevamento/subsidenza) e quelle con NMD anomalo. Questo è un risultato concettualmente importante perché ridimensiona, su scala stagionale, il ruolo diretto del riscaldamento/raffreddamento adiabatico come controllore principale del melt: la fusione sembra invece rispondere più robustamente alla combinazione di advezione, temperatura in colonna e condizioni radiative/umidità, soprattutto nei settori dove il sistema è già predisposto (margini). È anche il motivo per cui gli autori insistono sul fatto che il 2002, pur avendo avuto una massima estensione di melt molto grande, sia spiegabile come somma di episodi caldi (a volte brevissimi) che spingono la fusione in aree insolite, mentre le medie JJA di alcune variabili (TCW, GBI/Z500) possono apparire relativamente “poco eclatanti” perché diluite nel trimestre.
Un ultimo tassello metodologico riguarda come si osserva la fusione: l’estensione di melt in molti dataset groenlandesi (inclusi prodotti MEaSUREs) deriva da microonde passive, dove la presenza di acqua liquida nella neve/firn modifica la brillanza e viene rilevata tramite soglie dinamiche su canali come il 37 GHz H-pol; questo rende possibile una diagnostica quotidiana su lunghi periodi, ma implica anche che “melt extent” sia una metrica di stato (presenza di melt) e non una misura diretta della quantità di fusione. Per questo, affiancare mappe come la Figura 5 a metriche modellate di melt amount o a diagnostiche energetiche è essenziale per non confondere arealità con intensità.
La Figura 6 ripete lo stesso “set diagnostico” della Figura 5, ma per giugno–agosto 2008 (JJA 2008): l’idea è legare dove e per quanto si è verificata la fusione sulla calotta groenlandese (tramite il numero di giorni di fusione, NMD) a quali condizioni atmosferiche medie stagionali possono averla favorita o inibita (campi ERA-Interim e relative anomalie rispetto alla climatologia JJA 2000–2014). È un approccio molto efficace perché mostra, in un’unica tavola, l’interazione tra circolazione sinottica/planetaria (Z500), termodinamica (TCW e T500), flussi vicino-superficie (u10, v10) e dinamica verticale (ω). Questo tipo di lettura integrata è al centro delle analisi sulla variabilità del melt groenlandese nel periodo 2000–2014.
Il primo messaggio che emerge è che il 2008 non appare come un’estate “uniformemente estrema”, ma come una stagione in cui il segnale di fusione è spazialmente selettivo: nel pannello (a) le anomalie di NMD (colori) mostrano aree con incrementi locali dei giorni di fusione e aree con riduzioni, e la sovrapposizione con le isolinee di Z500suggerisce che la distribuzione del melt sia stata modulata più dalla geometria della circolazione (dove si posizionano creste e saccature, quindi quali traiettorie di massa d’aria prevalgono verso la calotta) che da un semplice “riscaldamento diffuso” della colonna atmosferica. Questo è coerente con ciò che si osserva spesso in Groenlandia: gli anni con fusione più intensa e più estesa sono generalmente associati a un blocking groenlandese più marcato (anomalie positive di Z500/GBI), ma la risposta dipende molto anche dall’ampiezza delle anomalie termiche e dal posizionamento del pattern.
Il pannello (b) (NMD totale in JJA 2008) ricostruisce la “climatologia operativa” della fusione: i valori più alti si concentrano lungo le fasce marginali (soprattutto alle basse quote), mentre l’interno e l’altopiano centrale rimangono su valori molto inferiori. Questo non è un dettaglio grafico: significa che la metrica NMD è fortemente controllata dal fatto che le zone costiere e di ablazione possono rimanere in condizioni di melt per molte settimane, mentre l’interno spesso sperimenta melt solo in presenza di intrusioni calde eccezionali. In pratica, un’estate può produrre un segnale significativo in NMD lungo i margini anche senza spingere la fusione in modo massiccio sul plateau.
Il pannello (c) (anomalia di TCW, total column water, kg m⁻²) è particolarmente utile perché il contenuto di umidità della colonna atmosferica è uno dei modi più diretti per “vedere” il potenziale radiativo dell’atmosfera sopra la calotta: più vapore e/o più nubi (specie a bassa quota) tendono ad aumentare la radiazione infrarossa verso il basso, riducendo l’efficienza del raffreddamento radiativo della superficie e rendendo più probabile superare la soglia di fusione, soprattutto vicino a 0 °C. Nei lavori sul melt groenlandese, l’incremento di downwelling longwave e la presenza di masse d’aria più calde/umide sono spesso richiamati come ingredienti chiave degli anni con melt marcato. Nel 2008 la mappa di TCW mostra un pattern non omogeneo: ci sono settori con anomalie debolmente positive e altri con anomalie negative, suggerendo che l’umidità colonnare non abbia sostenuto un forcing “pan-ice-sheet”, ma abbia agito soprattutto su specifiche regioni, coerentemente con un segnale NMD non uniformemente positivo.
Questa lettura viene rafforzata dal pannello (d), che mostra l’anomalia di temperatura a 500 hPa (K). La temperatura in media troposfera è un indicatore utile perché aiuta a capire se le masse d’aria in arrivo sono, in media stagionale, più calde o più fredde della norma. Nel 2008 si nota un contrasto tra aree con anomalia più mite e aree più fredde, e questo si presta a una lettura fisica molto concreta: dove la colonna è più fredda, il sistema deve “lavorare di più” (via radiazione o advezione) per portare la superficie alla fusione, e quindi NMD tende a ridursi o a non crescere in modo significativo; dove la colonna è più mite, anche un forcing radiativo moderato può diventare sufficiente ad aumentare i giorni di melt nelle fasce predisposte (margini, basse quote). Anche nella reportistica climatica su Groenlandia, il 2008 viene descritto come un anno con segnali di melt importanti soprattutto in specifiche porzioni della calotta (in particolare più a nord), dentro un contesto che non implica necessariamente un’anomalia calda uniforme ovunque.
I pannelli (e–h) servono a decodificare come la circolazione vicino-superficie possa aver redistribuito calore e umidità. u10 (m s⁻¹) rappresenta la componente zonale del vento a 10 m (ovest–est), mentre v10 la componente meridiana (sud–nord); le mappe “mean” (e, g) descrivono il campo medio in JJA 2008, e le mappe “anomaly” (f, h) mostrano lo scostamento dalla climatologia JJA 2000–2014. Nel 2008 le anomalie dei venti non suggeriscono un’unica firma dominante su tutta la calotta, ma piuttosto una modulazione regionale: questo è coerente con la “mappa a mosaico” dell’NMD. In particolare, dove le anomalie indicano una riduzione della componente meridionale favorevole (cioè meno trasporto verso nord di aria mite) o un rafforzamento relativo di flussi più settentrionali, ci si aspetta una maggiore probabilità di contenimento del melt, specialmente nelle aree interne/alte quote dove la soglia di fusione è più difficile da raggiungere. Viceversa, dove il pattern di vento favorisce una maggiore frequenza di traiettorie capaci di convogliare aria relativamente mite o umida (anche senza un “super-blocking”), il NMD può aumentare in modo localizzato lungo le fasce marginali. È esattamente il tipo di meccanismo che, nella letteratura, viene usato per spiegare perché alcuni anni mostrino segnali forti in regioni nord-orientali o lungo margini specifici pur senza produrre un’estensione massima eccezionale dell’area in melt.
Il pannello (i) (anomalia di ω a 500 hPa, Pa s⁻¹) chiude il cerchio distinguendo ciò che è “verticale” da ciò che è principalmente “orizzontale”. In meteorologia, ω positivo indica subsidenza (moti discendenti), ω negativo indica sollevamento; su scala sinottica/stagionale la subsidenza può favorire riscaldamento adiabatico e cieli più stabili, mentre il sollevamento è spesso associato a raffreddamento adiabatico e maggiore nuvolosità/precipitazione. Nella Figura 6, come già accade spesso in queste analisi, non emerge un accoppiamento semplice “ω ↔ NMD”: e questo è un risultato importante, perché suggerisce che nel 2008 la fusione sia stata controllata più dal bilancio energetico guidato da temperatura e umidità in colonna + traiettorie di massa d’aria che non da un singolo meccanismo di riscaldamento/raffreddamento adiabatico stagionale. Questa conclusione è coerente con studi che, usando dataset satellitari e rianalisi, identificano nel trasporto di masse d’aria calde (e spesso più umide) e nei regimi di blocking/NAO il driver principale della variabilità della fusione groenlandese, con ω che agisce più come “modulatore” regionale che come controllore unico.
Un’ultima nota (molto pratica) riguarda il fatto che i campi di melt extent/NMD usati in questi lavori derivano frequentemente da microonde passive: sono preziosi perché forniscono una copertura quotidiana e di lungo periodo, ma restano misure di “stato di fusione” (presenza di acqua liquida nella neve/firn) più che misure dirette della quantità di melt. Per questo, leggere la Figura 6 insieme a variabili atmosferiche e, quando disponibile, a stime modellistiche di melt amount/runoff è essenziale per non confondere arealità con intensità. La documentazione MEaSUREs/NSIDC e studi metodologici recenti discutono in modo esplicito punti di forza e limiti dei dataset PMW (risoluzione spaziale, sensibilità alla wet snow, omogeneità sensoriale) e mostrano perché questi prodotti siano comunque centrali per ricostruire la variabilità del melt groenlandese.
In sintesi, il significato della Figura 6 è questo: nel 2008 la fusione estiva groenlandese non è stata un “evento uniforme”, ma il risultato di un forcing atmosferico che ha creato contrasti regionali nella colonna termica e nell’umidità (T500/TCW) e una circolazione vicino-superficie (u10/v10) capace di aumentare o ridurre localmente la frequenza di condizioni pro-melt, soprattutto lungo i margini. Proprio perché il segnale è regionale e non pan-ice-sheet, il 2008 diventa un ottimo caso studio per capire come si possa avere una stagione con anomalie di melt significative in alcune aree senza produrre necessariamente un massimo areale “da record”, e perché le metriche cumulative (giorni di fusione) e quelle di picco (estensione massima) possano divergere.
La Figura 7 riassume, per giugno–agosto 2010, come la fusione (qui descritta dal numero di giorni di fusione, NMD) sia stata il prodotto di un assetto atmosferico favorevole soprattutto lungo i margini della calotta, più che di un singolo forcing capace di “spingere” in modo sistematico la fusione verso l’interno. Tutte le anomalie atmosferiche mostrano scostamenti calcolati rispetto alla media JJA 2000–2014 e derivano da ERA-Interim, una rianalisi ampiamente usata per collegare circolazione e processi di superficie in Groenlandia.
Nel pannello (a) la distribuzione delle anomalie di NMD (colori) evidenzia un segnale nettamente positivo concentrato soprattutto nel settore sud-occidentale e, più in generale, lungo porzioni marginali della calotta, mentre il nord-estmostra segnali più deboli o localmente negativi; questo è coerente con l’idea che nel 2010 la stagione di fusione sia stata “sbilanciata” verso il versante occidentale e meridionale. Le isolinee di Z500 (anomalie di geopotenziale a 500 hPa) suggeriscono un assetto troposferico tendenzialmente favorevole a condizioni miti/stabili sul dominio groenlandese, un aspetto che la letteratura collega frequentemente a estati con melt significativo: non perché il geopotenziale alto “produca” di per sé la fusione, ma perché ridisegna la circolazione e aumenta la probabilità di traiettorie d’aria capaci di alimentare temperature in colonna più elevate e/o trasporti umidi verso la calotta.
Il pannello (b) (NMD totale in JJA 2010) mette in chiaro la geometria del processo: la maggior parte dei giorni di fusione si accumula nelle fasce a bassa quota e costiere, mentre l’interno resta molto meno “attivo” in termini di NMD. Questo spiega perché un anno possa risultare “forte” nelle metriche cumulative (giorni di fusione, estensione cumulativa) senza necessariamente avere un comportamento altrettanto estremo nella massima estensione areale: se la fusione resta prevalentemente confinata ai margini, può durare a lungo e accumularsi nel tempo senza occupare stabilmente l’altopiano interno. Una diagnostica indipendente basata su serie storiche di melt extent conferma che nel 2010 la durata della fusione è stata eccezionalmente lunga nel sud-ovest e più corta nel nord-est, coerentemente con il dipolo che si intravede anche qui.
I pannelli (c) e (d) sono, di fatto, il “motore termodinamico” del 2010: l’anomalia di TCW (contenuto totale d’acqua in colonna) risulta positiva su gran parte del dominio, mentre la temperatura a 500 hPa mostra anomalie calde diffuse, con un’eccezione più fredda verso il nord-est. Questa combinazione è molto efficace nel sostenere il melt marginale perché (i) una colonna più calda aumenta la probabilità che gli strati prossimi alla superficie si avvicinino alla soglia di fusione, e (ii) più vapore/nubi (che TCW spesso sottintende) tende ad aumentare la radiazione infrarossa verso il basso e a ridurre l’efficienza del raffreddamento radiativo, rendendo più “persistente” la fusione nelle regioni predisposte. È anche il contesto che aiuta a capire perché nel 2010 il melt sia stato intenso pur senza mostrare necessariamente i grandi “spike” areali tipici degli anni con forti intrusioni calde fino al plateau.
I campi di vento a 10 m (pannelli e–h) chiariscono come si sia distribuito il trasporto vicino-superficie. Le mappe medie (e, g) mostrano il pattern atteso attorno alla Groenlandia, mentre le anomalie (f, h) indicano che nel 2010 le deviazioni dai regimi medi non sono ovunque fortissime, ma sono coerenti con il segnale NMD lungo i margini: nel pannello (h)l’anomalia di v10 è tendenzialmente positiva sul margine orientale (northerlies più deboli del normale, quindi minore ventilazione fredda lungo costa est) e negativa su porzioni del margine occidentale (southerlies più deboli). Il fatto che, nonostante ciò, l’NMD resti positivo anche a ovest è un indizio che nel 2010 la “leva” principale lungo quel settore sia stata la componente termico-umida (TCW e T500 positive) più che una singola anomalia di vento; in altre parole, quando la massa d’aria è mediamente più calda e più umida, la fusione marginale può rimanere sostenuta anche se il vento non è anomalmente favorevole in modo uniforme.
Il pannello (i) (anomalia di ω a 500 hPa) completa il quadro mostrando che non c’è un accoppiamento semplice e diretto tra moti verticali medi stagionali e pattern di NMD: le anomalie di ω sono relativamente deboli e frammentate, e non “ricalcano” in modo netto le aree a maggiore fusione. Questo è un risultato importante perché suggerisce che, nel 2010, la fusione sia stata controllata soprattutto da advezione e termodinamica in colonna (calore + umidità) e dai feedback di superficie, più che da un forcing adiabatico stagionale facilmente leggibile in ω.
Infine, per capire perché il 2010 sia spesso richiamato come anno di melt notevole, la figura va letta insieme ai feedback di superficie: studi dedicati al 2010 evidenziano che innevamento ridotto e comparsa precoce di ghiaccio nudopossono aver abbassato l’albedo, aumentando l’assorbimento di radiazione solare e rafforzando la fusione in un circuito di feedback positivo, particolarmente efficace nelle aree di ablazione del sud-ovest. In sintesi, la Figura 7 “racconta” un 2010 dominato da stagione lunga e marginale, alimentata da T500 e TCW positive e da un assetto di circolazione coerente con condizioni favorevoli al melt, con un nord-est relativamente meno coinvolto: un caso quasi ideale per capire perché le metriche cumulative possano risultare elevate anche quando l’estensione massima non è necessariamente estrema.
La Figura 8 condensa, per giugno–agosto 2012, il “set minimo” di variabili necessario a spiegare perché quell’estate sia stata eccezionale non solo per l’area massima in fusione, ma soprattutto per la diffusione spaziale e la persistenzadei giorni di melt su gran parte della calotta. Il pannello (a) mostra che le anomalie positive di NMD (giorni di fusione) sono ampie e coerenti con una forte anomalia nella circolazione di media troposfera: le isolinee di Z500 indicano un assetto anticiclonico/bloccante sul dominio groenlandese, cioè una configurazione che tende a deviare il flusso medio e a favorire l’afflusso di masse d’aria più miti verso la calotta, oltre a condizioni favorevoli al riscaldamento della colonna atmosferica. Questo tipo di “Greenland blocking” è un tratto ricorrente delle estati recenti ad alto melt ed è spesso associato a fasi di NAO estiva negativa e a un GBI elevato; nel 2012 tale segnale è stato particolarmente marcato e persistente.
Il pannello (b) (NMD totale in JJA 2012) evidenzia che l’incremento dei giorni di fusione non resta confinato alle sole fasce marginali: pur mantenendo i massimi lungo i bordi (dove la stagione è naturalmente più lunga), il segnale invade in misura anomala anche settori interni e settentrionali. La “spina dorsale” termodinamica di questa risposta è mostrata dai pannelli (c) e (d): l’anomalia positiva di TCW (contenuto totale di acqua nella colonna) è estesa e robusta, e l’anomalia di temperatura a 500 hPa è fortemente positiva su gran parte della Groenlandia. In altre parole, nel 2012 non si tratta soltanto di un riscaldamento superficiale locale: è la colonna troposferica ad essere anormalmente calda e umida, una combinazione che aumenta la probabilità di superare la soglia di fusione anche dove normalmente il melt è raro. Questa cornice è coerente con la letteratura sul 2012, che documenta l’evento estremo di metà luglio (fusione “a” o “vicino” alla superficie su circa 98.6% dell’estensione della calotta il 12 luglio) e collega il riscaldamento diffuso a un trasporto sinottico molto efficiente di calore/umidità, includendo il contributo di un Atmospheric River intenso sul lato occidentale groenlandese.
I campi di vento a 10 m (pannelli e–h) vanno letti come risposta/manifestazione della circolazione anomala più che come unico driver “meccanico” del melt: la struttura di u10 e v10 e delle loro anomalie è compatibile con una circolazione anticiclonica rafforzata attorno alla Groenlandia, che rimodula le traiettorie delle masse d’aria e contribuisce a sostenere anomalie calde e umide sul dominio. In questo contesto, l’effetto del vapore acqueo e delle nubi diventa cruciale per il bilancio energetico superficiale: non è necessario invocare solo cieli sereni, perché nel luglio 2012 è stato mostrato che nubi basse liquide possono amplificare la fusione aumentando la radiazione infrarossa discendente e le temperature prossime alla superficie durante le fasi critiche dell’evento, rafforzando quindi l’estensione e la durata del melt.
Infine, il pannello (i) (anomalia di ω a 500 hPa) è importante proprio perché non “spiega tutto”: la distribuzione delle anomalie di velocità verticale non risulta chiaramente co-localizzata con le massime anomalie di NMD, suggerendo che, su scala stagionale, il contributo di riscaldamento/raffreddamento adiabatico non sia stato il controllo dominante. Questo non esclude che a scala sinottica possano esistere aree di sollevamento forzato e subsidenza associate al blocking, ma indica che nel 2012 la fusione è stata sostenuta soprattutto dalla combinazione di circolazione bloccante, riscaldamento profondo della troposfera, anomalia positiva di umidità in colonna e feedback radiativi (inclusi quelli nuvolosi), che insieme rendono coerente la firma “quasi pan-Groenlandia” dei giorni di fusione mostrata dalla figura.
La Figura 9 (JJA 2013) è, di fatto, la “controprova dinamica” di quanto visto nel 2012: mentre l’anno precedente fu dominato da blocking e riscaldamento troposferico diffuso, qui la Groenlandia entra in una configurazione che sopprime la fusione su scala quasi pan-ice-sheet. Non a caso il 2013 viene spesso descritto come un ritorno a condizioni di melt più vicine alla climatologia di lungo periodo dopo l’estremo 2012, con estensione e frequenza di fusione nettamente ridimensionate.
Il segnale emerge già nei pannelli (a) e (b). In (a) le anomalie di NMD sono prevalentemente negative, coerenti con un campo di Z500 anomalo che suggerisce un regime più ciclonico/di saccatura sull’Artico centrale e sulla Groenlandia rispetto al 2012; in (b) l’NMD integrato su JJA resta sostanzialmente confinato alle fasce marginali e non “penetra” verso l’interno come negli anni ad alto melt. Questo assetto è compatibile con quanto riportato in letteratura: il 2013 è uno dei casi in cui la fusione ridotta si associa al segno “atteso” della teleconnessione, cioè una NAO estiva positiva, che tende a sfavorire il blocking groenlandese e quindi le intrusioni calde verso la calotta.
I pannelli (c) e (d) chiariscono perché il segnale di fusione sia così depresso: il TCW mostra anomalie per lo più negative (colonna più secca) e la temperatura a 500 hPa è diffusamente sotto media (troposfera più fredda). Questa combinazione è particolarmente efficace nel ridurre la probabilità di superamento della soglia di fusione, soprattutto sulle alte quote: aria più fredda in colonna rende più difficile portare la temperatura prossima alla superficie verso 0 °C, e una colonna più secca riduce il contributo “serra” del vapore acqueo alla radiazione infrarossa discendente, facilitando il raffreddamento radiativo della superficie (in media stagionale) e mantenendo più spesso condizioni sub-freezing. In termini di meccanismi, è esattamente il tipo di quadro che studi sinottico-climatologici associano a stagioni di melt contenuto: meno avvezione calda, colonna più fredda e meno supporto radiativo/termico alla fusione persistente.
Nei pannelli (e–h) le anomalie dei venti a 10 m risultano complessivamente deboli, ma il messaggio è coerente con l’idea di una ridotta alimentazione “mite” verso l’interno: dove diminuisce l’efficacia della componente meridionale favorevole (southerlies), l’NMD tende a calare, e questo è particolarmente penalizzante nel settore centrale e sul plateau nord-orientale, dove la fusione è spesso un equilibrio delicato tra advezione e bilancio radiativo. Il pannello (i) (ω a 500 hPa) non mostra una corrispondenza netta con l’NMD: anche quando compaiono zone di subsidenza/sollevamento, il segnale di fusione resta governato soprattutto dalla massa d’aria (fredda e relativamente secca) e dal pattern di circolazione su larga scala, più che da un forcing adiabatico stagionale “leggibile” direttamente in ω. Questo tipo di non-collocazione tra ω e melt è un risultato abbastanza tipico quando la variabilità del melt è dominata da trasporto di calore/umidità e dal bilancio radiativo più che da moti verticali medi stagionali.
Un ulteriore elemento che rende il 2013 “opposto” al 2012 è lo stato della superficie: analisi indipendenti indicano nel 2013 albedo più alta lungo le zone costiere rispetto ad anni recenti, compatibile con meno neve bagnata e minore esposizione di ghiaccio nudo (ablation area ridotta), cioè un indebolimento del feedback positivo dell’albedo che nel 2012 aveva amplificato la fusione. Anche il ruolo delle nubi, quando quantificato in modo esplicito, appare secondario nel senso che può modulare il melt ma non invertire il segno del regime: ad esempio, stime basate su diagnostiche radiative indicano che nel 2013 l’effetto “cloud-driven” sull’estensione di fusione è piccolo rispetto al 2012, coerentemente con un anno dominato soprattutto da colonna fredda/secca e circolazione sfavorevole.
3.3. Statistiche tra le variabili di fusione e le condizioni atmosferiche
In questa analisi viene quantificato quanto e come quattro diversi descrittori della fusione estiva sulla calotta groenlandese (NMD: numero di giorni di fusione; estensione cumulativa della fusione; estensione massima; quantità di fusione modellata) rispondano a un set di forzanti atmosferiche che includono indici di circolazione a grande scala (NAOI e GBI), diagnosi della colonna atmosferica (TCW, cioè il contenuto totale di acqua nella colonna: vapore + idrometeore) e informazioni sulla provenienza delle masse d’aria, oltre alle componenti del vento al suolo. Il quadro che emerge è molto istruttivo perché mostra subito un punto chiave: non esiste una “variabile meteorologica unica” che spieghi allo stesso modo tutte le metriche di melt, e anzi la sensibilità dipende dal fatto che alcune variabili integrano l’intera stagione (NMD, estensione cumulativa, melt modellato) mentre altre catturano soprattutto episodi sinottici rapidi (estensione massima). Questa distinzione di scala temporale è coerente con la letteratura che separa i controlli stagionali legati a persistenza di ridge/blocchi e contenuto di umidità, dai controlli “event-driven” legati a singole intrusioni calde/umide o a configurazioni bariche di breve durata.
Dal punto di vista statistico, il risultato più robusto è il dominio del TCW: la covariazione con le variabili di fusione integrate stagionalmente è molto alta (per NMD estivo medio su GrIS il valore riportato di r·|r| è ~0,83; ~0,84 per l’estensione cumulativa; ~0,82 per la fusione modellata). In termini fisici questo è plausibile perché TCW è un proxy sintetico di un “ambiente favorevole al melt” che combina (i) trasporto di umidità e aria mite, (ii) presenza e persistenza di nubi con contenuto di acqua liquida sufficiente a incrementare la radiazione infrarossa verso la superficie e (iii) condizioni sinottiche che sostengono tali configurazioni per più giorni. La letteratura sulla radiativa delle nubi in Artico/GrIS spiega bene perché, in presenza di albedo elevata (tipica di gran parte della calotta e del plateau interno), il bilancio netto delle nubi tenda spesso a essere di riscaldamento: la componente longwave (effetto serra delle nubi) può dominare sulla riduzione shortwave, specialmente quando le nubi sono otticamente spesse nel longwave (tipicamente legato a water path liquido “sufficiente”).
Questo meccanismo è documentato con misure e sintesi molto solide: Shupe & Intrieri mostrano il ruolo congiunto di proprietà delle nubi, albedo superficiale e geometria solare nel determinare il segno del cloud radiative forcing in ambiente artico; Miller et al. evidenziano che a Summit, proprio per l’albedo persistentemente alta, l’effetto netto delle nubi è dominato dal longwave e risulta spesso positivo (riscaldante) anche in estate; Van Tricht et al. quantificano su scala di calotta che le nubi possono aumentare il runoff da fusione in modo sostanziale (ordine “frazione rilevante” del totale rispetto a cieli sereni), legandolo a un cloud radiative effect medio positivo. Su porzioni di ablation zone occidentale, van den Broeke et al. mostrano come la radiazione e la copertura nuvolosa modulino la disponibilità energetica per la fusione con forte eterogeneità spaziale, rendendo naturale che le correlazioni con TCW (che incorpora anche informazione su nubi/umidità) risultino più “pulite” e coerenti di quelle con indici di circolazione più indiretti.
Gli indici di circolazione (NAOI e GBI) mostrano invece relazioni più deboli e soprattutto più “localizzate” nello spazio. La NAO estiva, quando negativa, è classicamente associata a configurazioni che favoriscono blocking/alte pressioni in area groenlandese e avvezione calda meridionale lungo la Groenlandia occidentale, con impatti documentati su temperature superficiali e fusione. Tuttavia, nello specifico, la correlazione tra NAOI medio JJA e NMD risulta in gran parte non significativa nel periodo analizzato: un’interpretazione fisicamente sensata è la non-linearità della metrica NMD, che “risponde” in modo marcato solo quando la temperatura oscilla attorno allo zero; se una regione è frequentemente sopra soglia, NMD tende a saturare e perde potere discriminante (e viceversa, se resta quasi sempre sotto soglia, NMD resta basso). Questo spiega anche perché ai margini meridionali, dove NMD può risultare vicino al massimo anche nelle estati relativamente “meno fusionali”, le correlazioni con Z500/GBI/NAO si indeboliscano o scompaiano: non perché la circolazione non conti, ma perché la variabile osservata ha poca varianza utile (problema statistico di “range restriction”). L’idea che gli indici di circolazione catturino meglio la predisposizione sinottica (ridging/blocco) che non la risposta locale soglia-dipendente è coerente con molti studi sulla fusione groenlandese e sui ridge a 500 hPa.
Il GBI, che misura la tendenza al blocking sulla Groenlandia tramite la pressione/altezza geopotenziale a media troposfera su un dominio dedicato, presenta una correlazione positiva con NMD ma più debole e confinata (bande strette su Nord e Ovest). È un esito in linea con l’idea che il blocking sia spesso un “abilitatore” (crea persistenza di ridge, cieli più stabili, traiettorie di avvezione) ma non determini da solo l’intensità della fusione se non è accompagnato da adeguato apporto di calore/umidità e da feedback superficiali (albedo, stato del firn, ecc.). Inoltre, il fatto che GBI e NAOI siano fortemente anticorrelati, specialmente su alcune stagioni, è ampiamente documentato e rende naturale che i pattern di correlazione spaziale NMD–GBI e NMD–NAOI risultino quasi “speculari”, pur senza trasferirsi automaticamente alle altre metriche di melt.
Un elemento importante è la coerenza tra Z500 medio su GrIS e NMD: la correlazione positiva elevata, massima verso nord, suggerisce che l’intensificazione/persistenza di un ridge a media troposfera sia un controllo primario della “stagionalità della fusione” in quelle regioni, in accordo con gli studi che collegano direttamente la frequenza di configurazioni di ridge/blocco a un aumento della fusione e del melt extent su scala pluridecennale. In questo contesto rientrano bene sia gli eventi estremi (ad esempio il 2012, in cui configurazioni anomale di circolazione hanno accompagnato un melt eccezionale), sia gli episodi “asymmetric” come il 2015, quando una struttura di alta pressione/cutoff e la modulazione del getto hanno favorito fusione record a nord mentre la limitavano a sud: è un caso emblematico di come la circolazione possa redistribuire spazialmente la fusione, rendendo alcune metriche o aree molto sensibili e altre meno.
Il fatto che la temperatura a 500 hPa (T500) mostri coerenza qualitativa anno-per-anno con anomalie di NMD ma non produca correlazioni locali robuste è un altro segnale di complessità: la fusione superficiale non è funzione lineare della temperatura in libera atmosfera, perché passa attraverso (i) accoppiamento boundary layer–superficie, (ii) nubi e flussi radiativi, (iii) turbolenza e avvezione a bassa quota, e (iv) stato della superficie (albedo, presenza di acqua liquida nel firn, rugosità, ecc.). In altre parole, T500 può essere un indicatore sinottico utile, ma il “collo di bottiglia” energetico della fusione si decide spesso nei termini del bilancio energetico superficiale e della microfisica/struttura delle nubi, rendendo di nuovo TCW (e quindi il contesto termo-igrometrico colonnare) una variabile più efficiente dal punto di vista diagnostico.
Infine, l’estensione massima della fusione risulta la meno spiegabile da singoli predittori stagionali, con un segnale significativo soprattutto legato alla provenienza delle masse d’aria verso Summit da nord-est. Questo è coerente con l’idea che il massimo melt extent sia spesso determinato da finestre sinottiche di pochi giorni (intrusioni calde/umide, configurazioni di flusso particolari, talvolta legate anche a trasporti intensi di vapore come gli atmospheric rivers), mentre le medie JJA di indici e variabili smussano proprio ciò che rende “massimo” l’evento. In anni recenti diversi lavori hanno mostrato che picchi di fusione, specie in Groenlandia nord-orientale, possono seguire episodi di trasporto intenso di vapore e configurazioni dinamiche capaci di innescare warming per compressione/foehn sul versante discendente della calotta, un quadro compatibile con la sensibilità a traiettorie d’aria specifiche più che a un singolo indice stagionale.
Tabella 2 – Relazioni statisticamente significative tra variabili di fusione e forzanti atmosferiche
La Tabella 2 riassume, in forma molto “selettiva”, quali legami tra atmosfera e fusione risultano davvero robusti dal punto di vista statistico. I valori riportati sono correlazioni “al quadrato con segno”, cioè una misura che conserva la direzione della relazione (positiva o negativa) ma ne enfatizza l’intensità, così da mettere in evidenza i controlli più forti. Inoltre, la tabella include solo le associazioni significative: * indica p < 0,01, ** indica p < 0,05. Un dettaglio metodologico importante è che NMD (numero di giorni di fusione) e le variabili atmosferiche sono mediati su JJA (giugno–agosto), mentre le metriche di estensione (cumulativa e massima) e la fusione modellata sono considerate sull’intera stagione di melt: questo rende già plausibile che alcune variabili spieghino meglio i “totali stagionali” e altre gli “estremi” di breve durata.
Il segnale dominante è il TCW (Total Column Water, mediato sulla calotta), che mostra correlazioni molto alte con le metriche più “integrate” della fusione: 0,83* con NMD estivo (JJA), 0,84* con l’estensione cumulativa e 0,82* con la fusione modellata. In termini fisici, un contenuto colonnare d’acqua elevato tende a coincidere con masse d’aria più miti e umide e, soprattutto, con condizioni nuvolose capaci di aumentare la radiazione infrarossa verso la superficie (cloud longwave effect), un meccanismo ben documentato in ambiente artico quando l’albedo è alta e il contributo longwave delle nubi può superare la perdita shortwave. Questa chiave radiativa aiuta a capire perché TCW “batta” gli indici dinamici: TCW è più vicino al processo energetico che porta la superficie a fondere (bilancio radiativo e turbolento) rispetto a NAO/GBI, che descrivono più indirettamente la configurazione sinottica. La letteratura osservativa e modellistica supporta fortemente questa lettura: il forcing radiativo delle nubi sull’Artico e la sua dipendenza da proprietà microfisiche e albedo è un risultato classico.
Gli indici di circolazione entrano, ma con un ruolo più “condizionale” e spazialmente eterogeneo. Il GBI compare in tabella solo per NMD, con 0,38*: segnala che estati con maggiore blocking groenlandese tendono ad avere più giorni di fusione, coerentemente con l’idea che i ridge a 500 hPa e la persistenza anticiclonica favoriscano condizioni favorevoli al melt (stagnazione, subsidenza, advezione calda sul fianco occidentale, ecc.). Tuttavia, il fatto che il GBI non risulti significativo per l’estensione cumulativa (nella tabella non compare) suggerisce che il blocking, da solo, non garantisce l’espansione progressiva della fusione su tutta la stagione: serve che la configurazione dinamica si accompagni a un adeguato apporto termo-igrometrico e a feedback superficiali (albedo/firn) che trasformano “potenziale sinottico” in fusione effettiva. Questo quadro è molto coerente con studi che hanno collegato l’aumento recente della fusione groenlandese alla maggiore frequenza di ridge in media troposfera e a pattern di circolazione favorevoli al melt.
Per la NAO, la tabella mostra una relazione più debole ma significativa con NMD (0,26)** e una relazione più marcata con la fusione modellata (0,42*). Il messaggio implicito è che la NAO agisce come “modulatore di contesto” (ad esempio favorendo avvezioni calde e pattern di pressione utili al melt), ma la risposta osservata nei soli giorni di fusione può essere smorzata da non-linearità e saturazioni della metrica (zone che fondono quasi sempre o quasi mai in JJA). Al contrario, una variabile modellata di melt integra quantità energetiche e idrologiche e può “seguire” più fedelmente la variabilità sinottica-climatologica legata alla NAO. Non a caso, lavori osservativi e diagnostici mostrano che fasi più negative della NAO sono associate a maggiore blocking e a un aumento della frequenza di flussi meridionali caldi lungo l’ovest della calotta, condizioni ricorrenti in estati ad alta fusione.
Infine, la metrica più “estrema”, la massima estensione di fusione, risulta spiegata in modo significativo solo dalla frequenza di masse d’aria di origine nord-orientale a Summit, con una correlazione negativa di −0,30 (p < 0,05). Tradotto: quando a Summit è più frequente l’arrivo di aria da NE, l’estensione massima di fusione tende a ridursi, un segnale compatibile con un contributo più freddo/secco o comunque meno favorevole agli episodi di melt più estesi. Questo è un punto molto importante perché ricorda che i picchi di melt non sono necessariamente “la media stagionale portata all’estremo”: spesso dipendono da specifiche traiettorie d’aria e configurazioni sinottiche di breve durata (intrusioni calde/umide, cutoff highs, ridge posizionati in modo particolare), come mostrato da analisi dedicate ai grandi eventi di fusione e alle loro traiettorie lagrangiane e, per casi recenti, al ruolo di cutoff high artici nel redistribuire la fusione verso latitudini più alte.
Figura 10 – Correlazioni spaziali tra i giorni di fusione estivi (NMD) e la circolazione/termodinamica atmosferica sulla Groenlandia (2000–2014)
La Figura 10 mappa, punto per punto sulla calotta groenlandese, quanto il numero di giorni di fusione in estate (NMD, JJA) covari con cinque descrittori atmosferici: NAO, Greenland Blocking Index (GBI), geopotenziale a 500 hPa (Z500), temperatura a 500 hPa (T500) e Total Column Water (TCW). La scelta di rappresentare la correlazione come “correlazione al quadrato con segno” serve a leggere in un colpo solo sia direzione (rosso = relazione positiva; blu = negativa) sia forza del legame; inoltre, la figura segnala che valori con magnitudine > |0,25| sono significativi al 95%. In altre parole: dove compaiono colori intensi e diffusi, lì c’è un’impronta atmosferica statisticamente robusta sui giorni di fusione; dove il segnale sbiadisce o si annulla, la relazione è debole, non lineare o mascherata da fattori locali (soglia di fusione, albedo, stato del firn, ecc.).
Nel pannello (a) NMD–NAO domina un segnale negativo (tinte azzurre), coerente con l’idea che fasi estive più negative della NAO favoriscano configurazioni sinottiche capaci di aumentare la fusione sulla Groenlandia, soprattutto attraverso una riorganizzazione del flusso atlantico che facilita intrusioni più miti verso il settore occidentale e condizioni di blocco/ridging in area groenlandese. Questa connessione tra NAO estiva, anomalie termiche superficiali e melt extent è ampiamente discussa nella letteratura basata su MODIS e reanalisi, in cui la NAO negativa viene spesso associata a estati “calde” e a feedback di albedo più efficaci.
Il pannello (b) NMD–GBI ribalta il segno: la mappa è prevalentemente positiva (rossi), indicando che un blocking più marcato sulla Groenlandia tende ad accompagnarsi a più giorni di fusione. È un risultato molto fisico: il blocking in media troposfera implica una dorsale quasi-stazionaria, subsidenza, maggiore persistenza sinottica e spesso un’anomalia di spessore/temperatura della colonna che “spinge” il sistema verso condizioni favorevoli al melt. La centralità del blocking nel modulare estati estreme (con il caso 2012 come riferimento classico) è stata documentata con analisi dinamiche dedicate, che collegano ridge persistenti a episodi di fusione eccezionale su larga scala. Un punto importante, che la figura lascia intuire, è che NAO e GBI non sono due “leve” indipendenti: su molte stagioni mostrano una forte anticorrelazione (più netta in inverno, ma presente anche in estate), quindi una parte del segnale NAO in (a) è, di fatto, l’altra faccia del segnale di blocking in (b).
Il pannello (c) NMD–Z500 è, dal punto di vista dinamico, il cuore della figura: le correlazioni positive risultano diffuse e spesso più convincenti che non per la NAO, perché Z500 è una misura più “diretta” della presenza di una dorsale/anticiclone in media troposfera rispetto a un indice teleconnettivo. Qui il messaggio è chiaro: quando l’estate groenlandese è dominata da anomalie positive di geopotenziale a 500 hPa (ridging), i giorni di fusione aumentano su ampie porzioni della calotta. Questo si allinea con studi che ricostruiscono i driver sinottici del melt groenlandese tramite classificazioni dei pattern di Z500 (ad esempio con SOM o analoghi), mostrando che frequenza e intensità dei regimi anticiclonici spiegano molta parte della variabilità interannuale della fusione estiva.
Il pannello (d) NMD–T500 mostra invece un segnale più debole e meno coerente rispetto a Z500/GBI. È un punto molto interessante perché “smonta” l’idea (troppo semplice) che basti una troposfera più calda a 500 hPa per garantire più melt al suolo. La fusione superficiale è una variabile soglia, governata dal bilancio energetico alla superficie: radiazione shortwave/longwave, flussi turbolenti, stato dell’albedo e presenza di acqua liquida nel firn possono amplificare o attenuare l’effetto di una anomalia termica in libera atmosfera. In pratica T500 aiuta a descrivere l’ambiente sinottico, ma non è necessariamente il “collo di bottiglia” energetico che decide quanti giorni superi davvero lo zero alla superficie.
Il pannello (e) NMD–TCW è quello che, in molte aree (in particolare lungo il settore occidentale e su parti dei margini), mostra il legame più netto: più acqua nella colonna atmosferica = più giorni di fusione. Questo è coerente con una linea di evidenza molto robusta: in ambiente artico e su superfici ad albedo elevata, le nubi (e più in generale l’umidità colonnare) possono produrre un effetto radiativo netto di riscaldamento perché l’incremento di longwave verso il suolo può superare la perdita shortwave, soprattutto con determinate proprietà microfisiche e geometrie solari. Sulla Groenlandia, inoltre, osservazioni e simulazioni mostrano che le nubi possono aumentare in modo sostanziale il runoff da fusione rispetto a condizioni di cielo sereno, proprio tramite un cloud radiative effect medio positivo; e nelle zone di ablazione occidentali la modulazione radiativa e nuvolosa è stata misurata direttamente con serie pluriennali di bilancio radiativo. In questo senso TCW funziona spesso da proxy “integrato” (umidità + potenziale nuvoloso + trasporto) più vicino al processo fisico che produce melt rispetto ai soli indici dinamici.
Guardando la figura nel suo insieme, emerge una gerarchia interpretativa piuttosto pulita: Z500/GBI descrivono l’assetto dinamico che rende probabile la fusione (persistenza anticiclonica e ridging), mentre TCW cattura la componente termodinamico-radiativa che spesso decide la “conversione” di quell’assetto in giorni effettivi di melt. Il fatto che alcune aree mostrino correlazioni più deboli o discontinue (specie in porzioni del plateau interno o in margini dove la variabilità dei giorni di fusione è ridotta) è perfettamente compatibile con limitazioni statistiche e fisiche: se NMD varia poco (perché quasi sempre vicino a zero o vicino al massimo stagionale), la correlazione può collassare anche se il controllo atmosferico esiste; inoltre, la risposta locale dipende da feedback di albedo e stato del firn che non sono riassunti da una singola variabile in quota o da un solo indice. In sintesi, la Figura 10 non dice soltanto “quale indice conta di più”, ma mostra che i giorni di fusione estivi sulla Groenlandia sono l’esito di un accoppiamento tra dinamica sinottica (ridging/blocco) e forzanti radiative-umide (nubi/TCW), con un ruolo spesso secondario (o comunque non lineare) della sola temperatura a 500 hPa.
4. Discussione
In questo studio vengono esaminate le relazioni tra quattro descrittori della fusione della calotta groenlandese (GrIS) — numero di giorni di fusione (NMD), quantità di fusione modellata, estensione cumulativa della fusione ed estensione massima della fusione — e le condizioni atmosferiche concomitanti nel corso di 15 stagioni di melt. L’idea di fondo, supportata dai risultati, è che queste variabili non siano intercambiabili: pur descrivendo tutte “la fusione”, lo fanno su scale temporali e spaziali diverse, e quindi rispondono in modo differente alle forzanti atmosferiche. Questo aspetto diventa particolarmente rilevante in un periodo recente contraddistinto da un riscaldamento marcato della GrIS e da un aumento della frequenza di configurazioni sinottiche favorevoli al melt, spesso collegate a fasi più negative della NAO e a maggiore blocking groenlandese.
Una parte importante della letteratura precedente ha affrontato stagioni eccezionali o singoli eventi di fusione usando una sola metrica, oppure tramite confronti fra simulazioni e osservazioni; altri studi hanno invece privilegiato un approccio “da serie temporale” per individuare tendenze nelle variabili di melt. Nel lavoro di Fettweis e colleghi sul periodo 1988–2003, ad esempio, emergono trend positivi nelle aree di fusione stimate da microonde passive e prodotti modellistici, ma con la nota cruciale che anni di massimi e minimi delle diverse metriche non coincidono necessariamente, suggerendo già allora una variabilità interannuale non univoca tra “quantità” e “geometria” della fusione. Su una finestra più lunga, Mernild et al. (1960–2010) mostrano l’aumento di estensione, frequenza e durata della fusione, evidenziando l’estremità del 2010, ma senza entrare nel dettaglio della dinamica atmosferica associata alle stagioni meno estreme.
L’analisi qui discussa chiarisce bene perché le metriche possono divergere: NMD, fusione modellata ed estensione cumulativa sono grandezze integrate che “accumulano” informazione lungo tutta la stagione di melt, mentre l’estensione massima è tipicamente legata a eventi episodici su scala di pochi giorni. Di conseguenza, è frequente osservare stagioni in cui NMD, fusione modellata ed estensione cumulativa risultano elevati, ma l’estensione massima rimane relativamente contenuta: in tali casi la fusione è intensa ma resta confinata soprattutto ai margini (zone di ablazione e basse quote), senza propagarsi in modo esteso verso l’interno e l’altopiano centrale. Questa asimmetria è coerente con la fisica del sistema: per “portare” la fusione sull’alta calotta serve spesso un evento sinottico ben organizzato (intrusione calda/umida, ridge persistente con avvezione meridionale, nubi liquide sottili ma radiativamente molto efficaci), mentre per aumentare NMD ai margini può bastare un’anomalia stagionale moderata ma persistente.
Un punto raffinato, spesso sottovalutato, è che gli eventi che determinano l’estensione massima sono tanto più efficaci quanto più la stagione precedente ha predisposto la superficie: se gli strati superiori del manto nevoso/firn sono già prossimi alla temperatura di fusione, una successiva avvezione calda e umida ha maggiore probabilità di produrre un’espansione rapida del melt verso quote elevate. Il caso del 2012 è emblematico: un episodio di riscaldamento e umidificazione della colonna atmosferica, associato a nubi liquide in quota e condizioni sinottiche particolari, portò temperature superficiali oltre lo zero anche a Summit, con fusione prossima al pan-ice-sheet. Questo tipo di evento è oggi frequentemente interpretato nel quadro dei trasporti intensi di vapore (atmospheric rivers) e delle loro interazioni con la topografia della calotta e i processi di nube, che possono amplificare la radiazione infrarossa discendente e, in alcune regioni, innescare anche effetti favonici sul lato discendente.
Per quanto riguarda NMD, la figura sinottica che emerge è quella di un legame più stretto con le condizioni atmosferiche medie estive. Negli anni con anomalie positive di NMD, le anomalie tendono a co-localizzare con: valori elevati di TCW (contenuto totale d’acqua nella colonna), geopotenziale a 500 hPa più alto (ridging/blocco), temperature troposferiche più miti fino a 500 hPa e spesso venti meridionali a bassa quota. Questo insieme di segnali è coerente con estati dominate da NAO negativa e GBI positivo, cioè con una circolazione più ondulata e propensa al blocking groenlandese, che facilita sia l’avvezione calda sia la persistenza delle configurazioni favorevoli alla fusione.
Dentro questo quadro, TCW risulta un predittore particolarmente “efficiente” perché condensa informazione termodinamica e radiativa: valori elevati di TCW implicano spesso masse d’aria più miti e umide e una maggiore probabilità di nubi con contenuto d’acqua liquida sufficiente a incrementare la radiazione infrarossa verso la superficie. In ambiente artico, il segno del cloud radiative forcing dipende da proprietà delle nubi, albedo e geometria solare; su superfici molto riflettenti, l’effetto longwave delle nubi può dominare e produrre un riscaldamento netto della superficie, soprattutto fuori dal picco estivo o in presenza di specifiche microfisiche. Sulla Groenlandia, osservazioni e approcci combinati (satelliti attivi, modelli climatici e modelli di neve) mostrano che le nubi possono aumentare sensibilmente il runoff da fusione rispetto a cieli sereni, quantificando un effetto radiativo medio positivo non trascurabile.
La distribuzione spaziale del segnale è altrettanto informativa: la correlazione tra TCW e NMD tende a intensificarsi lungo i margini occidentali e settentrionali, aree in cui TCW può diventare episodicamente elevato per incursioni di aria calda e umida legate a traiettorie provenienti dal Nord America e dall’Atlantico nord-occidentale. In questa lettura, TCW funziona come proxy di “aria calda e nuvolosa” e rende coerente il collegamento proposto tra avvezione calda, forcing radiativo positivo delle nubi e riscaldamento recente particolarmente marcato nel settore nord-occidentale della Groenlandia. Un elemento ulteriore, utile per interpretare gli anni 2007–2014, è che l’aumento della fusione può essere associato a un incremento della radiazione infrarossa discendente anche in presenza di una riduzione della frazione di copertura nuvolosa, scenario compatibile con nubi meno frequenti ma mediamente più “efficaci” radiativamente perché più ricche in acqua liquida: questa chiave termodinamica è stata discussa esplicitamente come driver primario dell’aumento del melt in analisi dedicate, insieme all’aumento della frequenza di ridge a 500 hPa e al trasporto meridionale di calore e umidità.
Questi risultati sono coerenti anche con quanto riportato da Cullather e Nowicki (2018), che per la Groenlandia settentrionale associano la fusione superficiale ad anomalie positive della radiazione infrarossa discendente (downward longwave), anomalie negative della radiazione solare discendente (downward shortwave) e anomalie positive della frazione di copertura nuvolosa, in linea con un’ampia letteratura sul fatto che, su superfici ad albedo elevata, l’effetto “serra” delle nubi può prevalere sulla schermatura shortwave. Tuttavia, nel settore nord-orientaledella calotta la fusione risulta correlata nello stesso tempo con longwave discendente positiva, shortwave discendente positiva e copertura nuvolosa in diminuzione, un’apparente combinazione “mista” che mette in evidenza quanto l’impatto delle nubi sulla fusione non sia riducibile a un solo segno o a un solo meccanismo. In pratica, il segnale osservato può riflettere sia cambiamenti nella tipologia di nubi (fase, quota, temperatura di base, contenuto di acqua liquida) sia una diversa importanza relativa dei termini radiativi e turbolenti lungo il ciclo diurno, come mostrato in studi artici e groenlandesi che evidenziano come le nubi a bassa quota, liquide o miste, siano spesso quelle con maggiore capacità di aumentare il longwave al suolo.
A supporto di questa complessità, nello studio di caso sul forte evento di fusione del luglio 2012, Solomon et al. (2017)distinguono tra l’effetto radiativo diretto delle nubi e una serie di effetti indiretti dovuti all’interazione tra nubi, trasferimento radiativo, mescolamento turbolento tra cloud-top e superficie nevosa e flusso conduttivo nel manto nevoso. Sia Solomon et al. sia Van Tricht et al. (2016) insistono sull’importanza dell’effetto integrato sull’intero ciclo diurno, includendo la “precondizione” notturna del manto nevoso tramite il longwave nuvoloso: di notte, l’incremento del longwave può ridurre il raffreddamento radiativo della superficie e mantenere gli strati superiori più vicini a 0 °C, rendendo più facile superare la soglia di fusione durante il giorno successivo. Nel caso 2012, inoltre, lavori osservativi indicano un ruolo chiave di nubi liquide a bassa quota nel rafforzare il longwave discendente durante la fase critica dell’evento, favorendo temperature prossime alla fusione anche ad alta quota.
È utile anche separare concettualmente l’effetto delle nubi (e più in generale del contenuto d’acqua colonnare, TCW) sulla fusione superficiale dall’effetto sul bilancio di massa superficiale (SMB). In condizioni con anomalie negative di TCW la fusione tende a ridursi, ma tali configurazioni possono essere associate a regimi anticiclonici con precipitazioni ridotte e evaporazione/sublimazione più efficiente, elementi che possono comunque favorire un SMB più negativo (ad esempio riducendo l’accumulo o aumentando le perdite non legate alla fusione). Inoltre, mentre l’NMD medio sulla calotta risulta fortemente legato alla forza dell’anomalia di TCW soprattutto in prossimità delle coste (dove le intrusioni caldo-umide sono più frequenti), l’estensione della fusione tende a dipendere di più da quanto l’anomalia di TCW riesce a espandersi sull’intera GrIS: in altre parole, non conta solo “quanto” è anomala la colonna, ma anche “dove” e “quanto lontano” arriva, aspetti controllati dalla situazione sinottica e dall’origine delle masse d’aria, incluse le traiettorie di trasporto di vapore su scala sinottica.
Il ruolo del feedback dell’albedo nella fusione della Groenlandia non viene affrontato direttamente nell’analisi, pur essendo riconosciuto come un amplificatore importante degli estremi recenti. L’albedo diminuisce con l’invecchiamento e la metamorfosi della neve; nelle aree di ablazione lo scioglimento espone ghiaccio più scuro, che assorbe più radiazione solare e accelera ulteriormente la fusione. Per il periodo 2000–2010 sono state documentate tendenze alla diminuzione dell’albedo, in particolare nella zona di ablazione, con interpretazioni che chiamano in causa combinazioni di ridotta neve estiva, maggiore import di calore (anche via avvezione) e crescente esposizione di superfici “scure” e/o ricche di impurità e caratteristiche superficiali (polveri, black carbon, crioconite, alghe, acqua superficiale, crepacci), capaci di ridurre l’albedo e aumentare l’assorbimento shortwave. In questo contesto, lavori più recenti evidenziano anche che, man mano che l’albedo dell’ablation zone diminuisce, cambia il peso relativo dei termini radiativi e diventa ancora più cruciale quantificare correttamente il cloud radiative effect (competizione tra schermatura shortwave e intrappolamento longwave) per stimare la produzione di meltwater runoff. Gli autori motivano l’esclusione dell’albedo dalle analisi anche osservando che, al di sotto della linea di equilibrio, l’albedo è sensibile alla temperatura dell’aria prossima al suolo, mentre è meno sensibile alla temperatura superficiale stessa quando questa è “limitata” vicino a 0 °C durante la fusione.
Infine, la continuazione si inserisce nella cornice più ampia — già nota in letteratura — che associa NAO estiva negativa e blocking groenlandese positivo a stagioni con forte fusione. Studi diagnostici mostrano che le maggiori anomalie positive di fusione tendono a verificarsi insieme ad anomalie positive di geopotenziale in media troposfera in prossimità della Groenlandia e a pattern di vento/temperatura/umidità indicativi di intenso trasporto meridionale di calore e umidità. Fettweis et al. (2013) quantificano inoltre come cambiamenti nelle frequenze dei tipi di circolazione (con maggiore ricorrenza di configurazioni compatibili con NAO estiva negativa) possano spiegare una parte molto ampia del riscaldamento a 700 hPa sopra la Groenlandia nel periodo 1993–2012. Su scala giornaliera, Häkkinen et al. (2014) mostrano che il legame tra blocking e melt vale anche per eventi quotidiani, ma sottolineano che il solo anticiclone di blocco non basta a “spingere” la fusione verso quote elevate senza concomitanti anomalie termiche e avvezioni calde. Per il 2012 e le estati calde dal 2007 in poi, Hanna et al. (2014) descrivono un quadro di alta pressione di blocco persistente e NAO negativa, con avvezione calda sui fianchi occidentali e formazione di una “heat dome” favorevole alla fusione diffusa. In parallelo, analisi su lungo periodo come Rajewicz e Marshall (2014) mettono in relazione i positive degree days e le metriche di melt (NMD ed estensione da microonde passive) con le condizioni atmosferiche estive derivate da rianalisi, fornendo un ponte tra statistiche stagionali e dinamica sinottica dei regimi anticiclonici groenlandesi. Rajewicz e Marshall quantificano in modo molto esplicito il “peso” relativo dei driver atmosferici sulla variabilità estiva della fusione groenlandese: circa il 38–49% della variabilità dell’estensione di fusione in estate viene attribuita a anomalie di circolazione (comparsa di dorsali anticicloniche, configurazioni di blocking e relativa avvezione calda), mentre il riscaldamento di fondo della bassa troposfera spiegherebbe un ulteriore 13–27%. Questo tipo di scomposizione è importante perché separa (i) la componente “dinamica”, legata alla geometria e persistenza dei pattern sinottici che favoriscono il trasporto meridionale di calore/umidità e condizioni radiative favorevoli, da (ii) una componente “termodinamica” più lenta, che alza il livello di base su cui gli eventi sinottici si innestano. Un risultato compatibile è quello di Mioduszewski et al., che collegano le maggiori anomalie di melt a anomalie positive di altezza geopotenziale in prossimità della Groenlandia e a pattern di vento/temperatura/umidità indicativi di intenso trasporto meridionale; inoltre evidenziano un aumento di frequenza delle dorsali a 500 hPa nelle ultime decadi considerate, coerente con un contesto più favorevole al melt.
Un passaggio interessante della discussione riguarda la dinamica verticale: benché ci si aspetterebbe che anomalie di pressione e configurazioni anticicloniche si accompagnino spesso a segnali di subsidenza (e quindi potenziale riscaldamento adiabatico), qui non emerge una collocazione spaziale robusta tra anomalie di velocità verticale a 500 hPa e anomalie di NMD. L’interpretazione proposta – che il moto verticale sopra la calotta sia spesso più “imposto” dalla topografia (sollevamenti forzati, circolazioni orografiche e risposta del boundary layer al rilievo) che non dalle sole anomalie sinottiche – è in linea con lavori che mostrano come l’orografia groenlandese moduli fortemente nubi, circolazione regionale e risposta del melt al forcing su larga scala. In pratica, anche quando un contributo adiabatico esiste, può essere facilmente “coperto” da altri termini del bilancio energetico e termico (ad esempio avvezioni fredde compensanti o variazioni radiative), rendendo debole un legame semplice tra anticiclone e riscaldamento per subsidenza nella statistica spaziale.
Sul confronto tra indici di circolazione e variabili “termodinamico-radiative”, i risultati ribadiscono l’importanza di GBI e NAOI, ma con una differenza cruciale rispetto ad alcune letture precedenti: qui la relazione appare più netta con il numero di giorni di fusione (NMD) che non con l’estensione di fusione. Questo porta a un punto metodologico non banale: trattare l’NMD da rianalisi come semplice proxy dell’estensione di fusione osservata da satellite può essere utile in media, ma non è equivalente quando si vuole attribuire causalmente i driver atmosferici, perché NMD e melt extent possono rispondere a forzanti solo parzialmente sovrapposte (metriche cumulative/soglia-dipendenti contro metriche più “geometriche” e sensibili a eventi sinottici estesi). Nel periodo analizzato, le correlazioni tra NMD e GBI/NAOI e tra melt modellato e NAOI esistono ma sono comunque più deboli di quelle che legano le variabili di fusione al TCW/TCW-like forcing (acqua colonnare come proxy di masse d’aria miti/umide e condizioni nuvolose radiativamente efficaci). Il fatto che il melt modellato non correli con GBI e che l’estensione cumulativa non correli né con NAOI né con GBI rafforza l’idea che gli indici descrivano un “contesto dinamico” necessario ma non sufficiente, mentre l’umidità colonnare e i processi di nube/radiazione risultano più prossimi al meccanismo energetico che determina la fusione. Un riscontro molto solido in letteratura è che le nubi possono aumentare significativamente il runoff da fusione rispetto a cieli sereni tramite un cloud radiative effect medio positivo.
Un dettaglio regionale che emerge come peculiare è l’elevata correlazione locale tra Z500 e NMD nella Groenlandia nord-orientale (valori di r·|r| > ~0,60), suggerendo che in quell’area le variazioni del campo di pressione/geopotenziale su larga scala siano più determinanti per i giorni di fusione rispetto ad altre porzioni della calotta. Questo risultato, presentato come non ampiamente riportato in precedenza, è coerente con l’idea che nel settore settentrionale/NE piccoli spostamenti di dorsali e traiettorie d’aria possano cambiare in modo marcato sia l’avvezione in quota sia le condizioni di nuvolosità e radiazione, in un ambiente dove la frequenza di fusione è normalmente più bassa e quindi più sensibile a “finestre” favorevoli. La stessa letteratura recente sul blocking groenlandese sottolinea che non tutti i blocchi sono equivalenti: posizione e geometria della dorsale determinano impatti diversi su energia superficiale e trasporto meridionale.
Per il vento, la discussione è coerente con un quadro “a due livelli”: statisticamente, su scala interannuale, le correlazioni tra componenti del vento e NMD risultano deboli; però, nei singoli anni, le anomalie positive di NMD tendono comunque a co-occorrere con segnali di venti meridionali più forti (o venti settentrionali più deboli), cioè con configurazioni che implicano trasporto di aria più mite. In alcuni casi specifici – come il 2008 lungo i margini orientali – le anomalie positive di NMD risultano invece associate a un rafforzamento del vento zonale al suolo, interpretato come segnale di downslope winds lungo il versante orientale della calotta: qui il riscaldamento per compressione/adiabatico e, soprattutto, l’aumento del rimescolamento nello strato limite possono innalzare la temperatura vicino alla superficie e favorire la fusione costiera, talvolta dominando sull’apporto di aria più fredda proveniente dall’altopiano interno. Questo ruolo combinato di riscaldamento e mixing durante venti discendenti è coerente con studi sui venti catabatici/downslope e la loro interazione con le temperature near-surface, e con lavori che documentano eventi intensi di downslope in Groenlandia sud-orientale.
L’estensione massima della fusione, rispetto a NMD, estensione cumulativa e melt modellato, resta invece la metrica meno “spiegabile” con un singolo predittore meteorologico, il che è coerente con la sua natura: dipende spesso da episodi di pochi giorni, innescati da combinazioni sinottiche specifiche (advezioni calde/umide verso l’altopiano, nubi radiativamente efficaci, precondizionamento termico del manto nevoso). Rimane comunque un segnale statisticamente significativo: una correlazione negativa tra estensione massima e frequenza di masse d’aria nord-orientali a Summit, suggerendo che determinate traiettorie d’aria siano meno favorevoli agli estremi di estensione. Il confronto 2012–2013 descritto nella continuazione (maggiore peso di masse sud-occidentali nelle aree elevate del NE nel 2012 e maggiore quota di masse NE/SE nel 2013) è qualitativamente in linea con l’interpretazione di Hanna et al. sul 2012, dove blocking persistente e NAO negativa favorirono flussi caldi sui fianchi occidentali e condizioni favorevoli a fusione estesa.
In filigrana resta anche il tema della “predisposizione” della superficie, in particolare tramite feedback di albedo: la diminuzione dell’albedo (invecchiamento/metamorfosi della neve, esposizione di ghiaccio più scuro e impurità) amplifica l’assorbimento shortwave e può rendere più efficace l’impatto di eventi sinottici successivi, soprattutto nelle zone di ablazione. Questo tipo di feedback è stato quantificato e discusso con approcci basati su satelliti e modelli regionali, e viene spesso richiamato per interpretare gli estremi del decennio 2000–2010 e il record 2010.
Infine, la discussione richiama implicitamente una distinzione che in climatologia della Groenlandia è diventata centrale: capire la fusione richiede guardare a più metriche (perché “quanto”, “dove” e “per quanto tempo” non rispondono allo stesso modo), e collegare i driver del melt ai driver del bilancio di massa superficiale richiede includere anche variabili legate a precipitazione (quantità e fase) e teleconnessioni estive che modulano insieme temperatura, precipitazione e SMB.
0 commenti