Applicazione di tecniche remote sensing per lo studio dell'evoluzione e della dinamica criosferica in aree remote e di alta quota

Glaciers are sensitive climate indicators because they adjust their size in response to changes in climate (e.g. temperature and precipitation). The attention paid by the scientists to mountain glacier change is increasing as there are robust evidence of a global glacier shrinkage over the past five decades, which in turn is the consequence of global warming. Understanding the glacier response to climate change is of tremendous importance not only for improving scientific knowledge, but also to predict and manage water resources and natural risks for the people living in mountain areas in the short (e.g. glacier lake outburst floods), and long term (e.g. droughts). In this thesis are analysed different cryospheric elements (mainly glaciers and snow coverage) to describe their recent evolution and to look for relations, if any, with climate trends. Firstly, the focus is put on the Karakoram glaciers. Although a general worldwide retreat of mountain glaciers has been acknowledged by the scientific community, the Karakoram region represents an exception in this sense. Indeed, the net mass balance of the glaciers here in the early twenty-first century was slightly positive, and even some are expanding and thickening. This anomalous behaviour is known as Karakoram Anomaly. More precisely the study area is the Central Karakoram National Park (CKNP), a protected national park in Northern Pakistan representative of the glaciation of the whole Karakoram Range. The westerlies represent the dominant wind system and they occur during winter, while the neighbour Himalayan region is mainly influenced by the summer moonson. A comprehensive description of the state of the CKNP glaciers and of their recent evolution is presented. This was made after the compilation of the glacier inventory of the park for the years 2001 and 2010, which is also presented. Moreover, the analysis of the regional climate change in the recent years is also discussed and related to the actual glacier change, in order to understand the causes behind the Karakoram Anomaly. The glacier area change of the 711 glaciers mapped in the study zone during 2001–2010 was only -0.4 ± 202.9 km2 (over 4605.9 ± 86.1 km2 in 2001), evidencing a general stability. The climate analysis supports glacier stability in the area. A slight increase in late summer snow cover area during 2001–2010 was observed from MODIS snow data. At the same time, the available weather stations revealed an increase of snowfall events and a decrease of mean summer air temperatures since 1980, which would translate into more persistent snow cover during the melt season. These results support an enhanced glacier preservation in the ablation areas due to a long-lasting snow cover, and stronger accumulation at higher altitudes, pushing towards positive net balances. The other major aim of the present work is to provide a simple model to evaluate ice melt at the glacier surface. As the supraglacial-debris cover can alter ice ablation close to the glacier surface depending on its thickness, the model was made up of two parts: one which computes the ice melt over the bare ice areas using an enhanced T-index formula; and one for the debris-covered areas using a conductive heat flux module. For the debris-covered parts, the debris thickness map is produced and then provided to the model as input for the computation, other than the distributed shortwave incoming radiation. For the bare ice areas, the modeled air temperature and shortwave incoming radiation are derived from the automatic weather stations present in the CKNP and given to the model. The other model requirement is the digital elevation model. In particular, the meteorological input data were distributed starting from data acquired at Askole automatic weather station, located within the CKNP. The meteorological distribution was validated by comparison with data from other two AWS in the same park limits (Urdukas and Concordia). The modeled ablation data were in strong agreement with measurements collected in the field during 2011 on Baltoro glacier, which is representative of CKNP glaciers. Two sets of the same ablation dataset collected in the field in the CNKP area were used separately for calibration and validation. Snow melt was neglected since snow data in the study area was not systematically available. The model was run against the peak ablation season (23 July–9 August 2011), when meltwater mainly comes from ice melt, with snow thaw playing a minor role in this region. The total freshwater from the ablation areas of CKNP glaciers estimated by the model was 1.963 km3 (0.109 km3 d–1 on average). The meltwater from the debris-covered parts was 0.223 km3 (0.012 km3 d–1 on average; min–max 0.006–0.016 km3 d–1), and 1.740 km3 (0.097 km3 d–1 on average; min–max 0.041–0.139 km3 d–1) from debris-free sectors. The estimated total freshwater corresponds to 14% of the water contained in a large strategic dam along the Indus River, of which all the CKNP glaciers are tributaries. The sensitivity tests suggest that any increase in the extent of debris coverage (which will likely occur due to augmented macrogelivation processes and rockfall events), will affect melt depending on new debris thickness, and melting will increase largely if summer air temperature increases. The second major focus of this research is put on the snow cover variability of the Chilean Andes. A parallel major aim of this research work is to implement a methodology based on remote sensing to study the snow cover variation on an acceptable spatio-temporal resolution. The MODIS sensor was chosen as the most suitable for this purpose and a methodology for deriving snow maps automatically from it is described and applied for analyzing the SCA variation over 18 watersheds of the central Andes in Chile during 2008–2011. The same methodology was then adopted for the climate analysis in the CKNP as mentioned. The study area was divided into three sub-zones (Northern, Central, and Southern), for easing the computation of the snow analysis. Overall, SCA decreased during the four considered years. The maximum SCA was found in the Central Zone, while the topographic and climatic features (i.e. lower altitudes in the South, and a drier climate in the North), limited snow deposition elsewhere. The snow line was found higher in the Northern zone due to the presence of the plateau, while it decreases southwards. In the Northern Zone the minimum SCA was reached sooner than elsewhere, and it lasted for a longer period (November to March), probably because of the drier climate. West aspects showed the maximum of SCA in all zones throughout the study period. Finally, some examples of application of remote sensing to glacier related studies is presented for glaciers of various typology, size, and localization. Six case studies are shown, amongst which there are three alpine glaciers (Miage, Freney, Aletcsh), equatorial glaciers (the Kilimanjaro glaciers), the Harding Icefield in Alaska, and an Antarctic glacier (the Drygalsky Ice Tongue).

I ghiacciai sono efficaci indicatori climatici poichè si modificano in risposta ai cambi del clima (es. temperatura e precipitazioni). L'attenzione sui ghiacciai di montagna sta aumentando tra la comunità scientifica per via del loro sempre più evidente arretramento a scala globale negli ultimi cinquant'anni. Ciò è conseguenza del riscaldamento globale. Comprendere il comportamento dei ghiacciai in risposta al cambio climatico è di enorme importanza non solo per arricchire la conoscenza scientifica, ma anche per poter meglio gestire in futuro le situazioni di rischio naturale che possono colpire le popolazioni che vivono nelle zone montuose, sia nel breve termine (es. GLOF), sia nel lungo (es. Siccità). Questa tesi di dottorato analizza differenti aspetti della criosfera (ghiacciai e neve) per descriverne la variabilità recente e le relazioni con la dinamica climatica. Inizialmente ci si è concentrati sul Karakorum. Questa è un’area particolare per gli studi criosferici, che non segue i trend globali di regresso; infatti, in questa zona il bilancio di massa netto dei ghiacciai nei primi anni del ventunesimo secolo è stato leggermente positivo, con anche taluni casi di espansione. Questa eccezionale situazione è riconosciuta con il nome di Anomalia del Karakorum (Karakoram Anomaly). Più precisamente il presente elaborato si focalizza sulla zona del Central Karakoram National Park (CKNP), un'area protetta nel nord del Pakistan, rappresentativa della glaciazione dell'intera catena del Karakorum. In questa regione, i venti occidentali rappresentano il sistema di venti dominante e sono presenti nella stagione invernale, mentre la confinante regione Himalayana è sotto l'influenza predominante dei monsoni, che sono venti estivi. Il presente lavoro descrive in maniera completa lo stato dei ghiacciai del CKNP e la loro recente evoluzione. Ciò è stato possibile a seguito della compilazione del catasto glaciale del parco per gli anni 2001 e 2010, a sua volta descritto nel dettaglio nel presente elaborato. Inoltre è discussa l'analisi dei cambiamenti climatici poi messa in relazione con quelli glaciali, per poter comprendere le cause dietro l'Anomalia del Karakorum. Il cambiamento areale dei 711 ghiacciai mappati nell'area di studio è stato -0.4 ± 202.9 km2 (su 4605.9 ± 86.1 km2 nel 2001), il che evidenzia una generale situazione di stabilità. Anche l'analisi climatica supporta tale condizione di stabilità. Durante il periodo 2001–2010 si è osservato grazie ai dati del sensore MODIS un leggero aumento delle aree coperte da neve a fine estate. Allo stesso tempo, dati meteo dalle stazioni disponibili hanno rivelato un aumento delle nevicate e una diminuzione della temperatura media dell'aria in estate fin dal 1980, il che si tradurrebbe in coperture nivali più persistenti durante la stagione ablativa. Questi risultati vanno a favore della preservazione glaciale nelle zone di ablazione dovuta a una copertura di neve più duratura, e un maggiore accumulo a quote più alte, presupponendo bilanci di massa netti tendenti al segno positivo. L'altro principale obiettivo del presente elaborato di tesi è quello di fornire un modello di semplice utilizzo per quantificare l'ablazione di ghiaccio alla superficie glaciale. Dal momento che una copertura detritica sopraglaciale è in grado di alterare la fusione del ghiaccio vicino alla superficie in funzione dello spessore, il modello tiene conto di due diversi casi: una parte stima l'ablazione per le aree di ghiaccio scoperto con un metodo definito enhanced T-index; l'altra stima la fusione per le zone coperte da detrito, utilizzando un modello di flusso di calore conduttivo. Per quanto concerne le parti coperte da detrito, è stata prodotta una mappa degli spessori detritici che è poi stata usata come input per il modello, assieme alla radiazione solare entrante distribuita. Per le aree scoperte da detrito, sono state derivate la temperatura dell'aria e la radiazione entrante distribuite attraverso i dati delle stazioni meteo automatiche presenti nell'area, in seguito usate come input. L'altro parametro necessario è un modello di elevazione del terreno. In particolare, la distribuzione degli input meteorologici è stata validata con dati di altre due stazioni presenti all'interno del CKNP (le stazioni di Urdukas e Concordia). L'ablazione modellata è risultata essere fortemente concorde con le misurazioni effettuate sul ghiacciaio del Baltoro nel 2011, ghiacciaio rappresentativo di tutto il CKNP. Due campioni dello stesso set di dati di fusione misurati su terreno sono stati usati ciascuno rispettivamente in sede di calibrazione e validazione. La fusione nivale è stata ignorata dal momento che mancavano dati di neve sistematici nell'area di studio. Il modello è stato fatto girare durante il picco della stagione ablativa (23 luglio–9 agosto 2011), durante il quale l'acqua di fusione deriva primariamente dalla fusione glaciale, mentre quella nivale ha un ruolo decisamente minore in questa regione. Il modello ha calcolato un totale di acqua da fusione glaciale pari a 1.963 km3 (0.109 km3 al giorno in media). Quella derivante dalle parti coperte da detrito ammonta a 0.223 km3 (0.012 km3 al giorno in media; min–max 0.006–0.016 km3 al giorno), mentre per le parti a ghiaccio scoperto è 1.740 km3 (0.097 km3 al giorno in media; min–max 0.041–0.139 km3 al giorno). Tale quantità è paragonabile al 14% di tutta l'acqua contenuta in una grande diga strategica lungo il fiume Indo, di cui i ghiacciai del CKNP sono tributari. I test di sensitività del modello suggeriscono che un aumento delle superfici coperte da detrito sui ghiacciai (probabile per via dell'aumento di eventi di macrogelivazione e di frane) avrà un notevole impatto sulla fusione effettiva in funzione dei nuovi spessori detritici, e l'ablazione aumenterà sensibilmente se la temperatura dell'aria dovesse alzarsi. Successivamente l'attenzione del presente elaborato di tesi è concentrata sulle Ande Cilene e sulla variabilità della copertura nevosa. Un obiettivo principale parallelo della presente ricerca è stato infatti quello di individuare una metodologia basata sul telerilevamento per studiare la variazione della copertura nevosa ad una risoluzione spazio-temporale accettabile. Il sensore MODIS si è rivelato il più idoneo allo scopo ed è stata implementata una metodologia che permettesse di estrarre mappe di copertura di neve in maniera automatica dalle informazioni raccolte dal sensore stesso. In particolare, sono stati studiati diciotto bacini idrografici di montagna delle Ande centrali in Cile durante il periodo 2008–2011. La stessa metodologia è stata esportata e adottata per l'analisi della neve nel CKNP come detto. L'area di studio è stata divisa in tre sotto-zone (Settentrionale, Centrale, Meridionale), per alleggerire il carico di calcolo dell'analisi. In generale, l'area coperta da neve è diminuita nel corso dei quattro anni di riferimento. I valori massimi sono stati ritrovati nella zona centrale, mentre fattori topografici e climatici (i.e. quote basse più a sud e un clima più arido nel nord), hanno limitato la deposizione della neve nelle altre zone. La linea della neve è più alta nella zona settentrionale a causa della presenza dell'altopiano, e si abbassa via via verso la zona merdionale. Nella zona settentrionale i minimi di copertura nivale vengono raggiunti prima che nelle altre zone e durano più a lungo (da novembre a marzo), probabilmente a causa del clima più arido. Durante l'intero periodo i valori massimi di copertura nevosa si ritrovano verso ovest. Al termine dell'elaborato e pertinente al tema principale delle applicazioni del telerilevamento allo studio della criosfera, sono presentati alcuni esempi di analisi di ghiacciai di diversa tipologia, dimensione e area geografica. Si tratta di sei casi, fra cui sono presenti tre ghiacciai alpini (Miage, Freney, Aletcsh), ghiacciai equatoriali (i ghiaccia del Kilimajaro), l'Harding Icefield in Alaska e un esempio di ghiacciaio antartico (la Drygalsky ice Tongue).