One Decade of Glacier Mass Changes on the Tibetan Plateau Derived from Multisensoral Remote Sensing Data

The Tibetan Plateau (TP) with an average altitude of 4,500 meters above sea level is characterized by many glaciers and ice caps. Glaciers are a natural indicator for climate variability in this high mountain environment where meteorological stations are rare or non-existent. In addition, the melt water released from the Tibetan glaciers is feeding the headwaters of the major Asian river systems and contributes to the rising levels of endorheic lakes on the plateau. As many people directly rely on the glacier melt water a continuous glacier monitoring program is necessary in this region. In situ measurements of glaciers are important, but are spatial limited due to large logistical efforts, physical constrains and high costs. Remote sensing techniques can overcome this gap and are suitable to complement in situ measurements on a larger scale. In the last decade several remote sensing studies dealt with areal changes of glaciers on the TP. However, glacier area changes only provide a delayed signal to a changing climate and the amount of melt water released from the glaciers cannot be quantified. Therefore it is important to measure the glacier mass balance. In order to estimate glacier mass balances and their spatial differences on the TP, several remote sensing techniques and sensors were synthesized in this thesis. In a first study data from the Ice Cloud and Elevation Satellite (ICESat) mission were employed. ICESat was in orbit between 2003 and 2009 and carried a laser altimeter which recorded highly accurate surface elevation measurements. As in mid-latitudes these measurements are rather sparse glaciers on the TP were grouped into eight climatological homogeneous sub-regions in order to perform a statistical sound analysis of glacier elevation changes. To assess surface elevation changes of a single mountain glacier from ICESat data, an adequate spatial sampling of ICESat measurements need to be present. This is the case for the Grosser Aletschgletscher, located in the Swiss Alps which served as a test site in this thesis. In another study data from the current TanDEM-X satellite mission and from the Shuttle Radar Topography Mission (SRTM) conducted in February 2000 were employed to calculate glacier elevation changes. In a co-authored study, these estimates could be compared with glacier elevation changes obtained from the current French Pléiades satellite mission. In order to calculate glacier mass balances, the derived elevation changes were combined with assumptions about glacier area and ice density in all studies. In this thesis contrasting patterns of glacier mass changes were found on the TP. With an ICESat derived estimate of -15.6±10.1 Gt/a between 2003 and 2009 the average glacier mass balance on the TP was clearly negative. However, some glaciers in the central and north-western part of the TP showed a neutral mass balance or a slightly positive anomaly which was also confirmed by data from the current TanDEM-X satellite mission. A possible explanation of this anomaly in mass balance could be a compensation of the temperature driven glacier melt due to an increase in precipitation

Perennial Supraglacial Lakes in Northeast Greenland Observed by Polarimetric SAR

Supraglacial liquid water at the margins of ice sheets has an important impact on the surface energy balance and can also influence the ice flow when supraglacial lakes drain to the bed. Optical imagery is able to monitor supraglacial lakes during the summer season. Here we developed an alternative method using polarimetric SAR from Sentinel-1 during 2017-2020 to distinguish between liquid water and other surface types at the margin of the Northeast Greenland Ice Stream. This allows the supraglacial hydrology to be monitored during the winter months too. We found that the majority of supraglacial lakes persist over winter. When comparing our results to optical data, we found significantly more water. Even during summer, many lakes are partly or fully covered by a lid of ice and snow. We used our classification results to automatically map the outlines of supraglacial lakes, create time series of water area for each lake, and hence detect drainage events. We even found several winter time drainages, which might have an important effect on ice flow. Our method has problems during the peak of the melt season, but for the rest of the year it provides crucial information for better understanding the component of supraglacial hydrology in the glaciological system

Seasonal Observations at 79°N Glacier (Greenland) From Remote Sensing and in situ Measurements

This study investigates seasonal ice dynamics of Nioghalvfjerdsfjorden or 79°N Glacier, one of the major outlet glaciers of the North East Greenland Ice Stream. Based on remote sensing data and in-situ GPS measurements we show that surface melt water is quickly routed to the ice-bed interface with a direct response on ice velocities measured at the surface. From the temporally highly resolved GPS time series we found summer peak velocities of up to 22 faster than their winter baseline. These average out to 9 above winter velocities when relying on temporally lower resolved velocity estimates from TerraSAR-X intensity offset tracking. From our GPS time series we also found short term ice acceleration after the melt season. By utilizing optical satellite imagery and interferometrically derived digital elevation models we were able to link the post melt season speed-up to a rapid lake drainage event (<24 h) with an estimated drainage volume of 28x10⁶ m³. We further highlight that GPS measurements are needed to resolve short term velocity fluctuations with low amplitudes, whereas remote sensing estimates are rather useful for the calculation of general trends in velocity behavior

Arctic Meltponds: Automated Detection Algorithm Using Enhanced Machine Learning

Melt ponds are pools of water that form during summer on the surface of the arctic ice. Due to the lower albedo, melt ponds absorb more solar radiation than surrounding ice and hence have higher temperature. This causes more water to melt, creating a feedback loop. This means that melt pond fraction in ice sheets is an important factor to consider in global climate and sea ice models. In situ measurements are difficult and expensive in terms of time and labor. Furthermore, these measurements can only cover limited areas. This makes using Earth Observation methods for this task particularly attractive. Until today, there is no sophisticated global melt pond data set available: Accurate methods may exist for determining melt ponds from Sentinel-2 data. The downside of using Sentinel-2 is that parts of the High Arctic are not covered by this mission. MODIS data covers the whole globe at least once every three days, but the downside of it is that MODIS resolution is much coarser (250m vs. 10m). Since melt ponds are in general much smaller than 250m, it means that accurately capturing melt pond fraction from these data is difficult. We propose to address these issues by employing Deep Learning techniques. Namely, we use Sentinel-2 data to train a model to super-resolve MODIS images to higher resolution and to use all available MODIS bands and their surrounding pixels for information context when predicting melt pond and open water fractions. In addition, a thorough uncertainty quantification (UQ) will be applied by using the UQ Toolbox

Calving Induced Speedup of Petermann Glacier

This study assesses the response on ice dynamics of Petermann Glacier, a major outlet glacier in northern Greenland, to the 2012 and a possible future calving event. So far Petermann Glacier has been believed to be dynamically stable as another large calving event in 2010 had no significant impact on flow velocity or grounding line retreat. By analyzing a time series of remotely sensed surface velocities, we find an average acceleration of 10% between winter 2011/2012 and winter 2016/2017. This increase in surface velocity is not linear but can be separated into two parts, starting in 2012 and 2016 respectively. By conducting modeling experiments, we show that the first speedup can be directly connected to the 2012 calving event, while the second speedup is not captured. However, on recent remote sensing imagery newly developing fractures are clearly visible ∼12 km upstream from the terminus, propagating from the eastern fjord wall to the center of the ice tongue, indicating a possible future calving event. By including these fracture zones as a new terminus position in the modeling domain, we are able to reproduce the second speedup, suggesting that surface velocities remain on the 2016/2017 level after the anticipated calving event. This indicates that, from a dynamical point of view, the terminus region has already detached from the main ice tongue

High Mountain Asian glacier response to climate revealed by multi-temporal satellite observations since the 1960s

Funding: This study was supported by the Strategic Priority Research Program of Chinese Academy of Sciences (XDA20100300) and the Swiss National Science Foundation (200021E_177652/1). NN received funding from the European Union’s Horizon 2020 programme (No. 689443).Knowledge about the long-term response of High Mountain Asian glaciers to climatic variations is paramount because of their important role in sustaining Asian river flow. Here, a satellite-based time series of glacier mass balance for seven climatically different regions across High Mountain Asia since the 1960s shows that glacier mass loss rates have persistently increased at most sites. Regional glacier mass budgets ranged from −0.40 ± 0.07 m w.e.a−1 in Central and Northern Tien Shan to −0.06 ± 0.07 m w.e.a−1 in Eastern Pamir, with considerable temporal and spatial variability. Highest rates of mass loss occurred in Central Himalaya and Northern Tien Shan after 2015 and even in regions where glaciers were previously in balance with climate, such as Eastern Pamir, mass losses prevailed in recent years. An increase in summer temperature explains the long-term trend in mass loss and now appears to drive mass loss even in regions formerly sensitive to both temperature and precipitation.Publisher PDFPeer reviewe

Ice flow and the conditions of the ice-bed interface at the onset of the Northeast Greenland Ice Stream

The Northeast Greenland Ice Stream (NEGIS) is an important dynamic component contributing to the total mass balance of the Greenland ice sheet, as it reaches up to the central divide and drains 12% of the ice sheet. The ice stream geometry and surface velocities in the onset region of the NEGIS are not yet sufficiently well reproduced by ice sheet models. We present an assessment of the basal conditions of the onset region in a systematic analysis of airborne ultra-wideband radar data. Our data yield a new detailed model of ice-thickness distribution and basal topography in the upstream part of the ice stream. We observe a change from a smooth to a rougher bed where the ice stream widens from 10 to 60 km, and a distinct roughness anisotropy, indicating a preferred orientation of subglacial structures. The observation of off-nadir reflections that are symmetrical to the bed reflection in the radargrams suggests that these structures are elongated subglacial landforms, which in turn indicate potential streamlining of the bed. Together with basal water routing pathways, our observations hint to two different zones in this part of the NEGIS: an accelerating and smooth upstream region, which is collecting water, with reduced basal traction, and in the further downstream part, where the ice stream is slowing down and is widening, with a distribution of basal water towards the shear margins. Our findings support the hypothesis that the NEGIS is strongly interconnected to the subglacial water system in its onset region, but also to the subglacial substrate and morphology

Different mechanisms of Arctic first-year sea-ice ridge consolidation observed during the MOSAiC expedition

Sea-ice ridges constitute a large fraction of the ice volume in the Arctic Ocean, yet we know little about the evolution of these ice masses. Here we examine the thermal and morphological evolution of an Arctic firstyear sea-ice ridge, from its formation to advanced melt. Initially the mean keel depth was 5.6 m and mean sail height was 0.7 m. The initial rubble macroporosity (fraction of seawater filled voids) was estimated at 29% from ice drilling and 43%–46% from buoy temperature. From January until mid-April, the ridge consolidated slowly by heat loss to the atmosphere and the total consolidated layer growth during this phase was 0.7 m. From mid-April to mid-June, there was a sudden increase of ridge consolidation rate despite no increase in conductive heat flux. We surmise this change was related to decreased macroporosity due to transport of snow-slush to the ridge keel rubble via adjacent open leads. In this period, the mean thickness of the consolidated layer increased by 2.1 m. At the peak of melt in June–July we suggest that the consolidation was related to the refreezing of surface snow and ice meltwater and of ridge keel meltwater (the latter only about 15% of total consolidation). We used the morphology parameters of the ridge to calculate its hydrostatic equilibrium and obtained a more accurate estimate of the actual consolidation of the keel, correcting from 2.2 m to 2.8 m for average keel consolidation. This approach also allowed us to estimate that the average keel melt of 0.3 m, in June–July, was accompanied by a decrease in ridge draft of 0.9 m. An ice mass balance buoy in the ridge indicated total consolidation of 2.8 m, of which 2.1 m was related to the rapid mode of consolidation from April to June. By mid-June, consolidation resulted in a drastic decrease of the macroporosity of the interior of keel while the flanks had little or no change in macroporosity. These results are important to understanding the role of ridge keels as meltwater sources and sinks and as sanctuary for ice-associated organisms in Arctic pack ice