Elevation change and mass balance of Svalbard glaciers from geodetic data

This thesis uses ground-based, airborne and spaceborne elevation measurements to estimate elevation change and mass balance of glaciers and ice caps on the Svalbard archipelago in the Norwegian Arctic. Remote sensing data are validated against field measurements from annual campaigns at the Austfonna ice cap. A new and more accurate DEM of the ice cap is constucted by combining SAR interferometry with ICESat laser altimetry. The precision of the DEM is sufficient to correct ICESat near repeat-tracks for the cross-track topography such that multitemporal elevation profiles can be compared along each reference track. The calculated elevation changes agree well with more accurate elevation change data from airborne laser scanning and GNSS surface profiling. The average mass balance of Austfonna between 2002 and 2008 is estimated to -1.3 ± 0.5 Gt y-1, corresponding to an area-averaged water equivalent elevation change of -0.16 ± 0.06 m w.e. y-1. The entire net loss is due to a retreat of the tidewater fronts. Earlier time periods are difficult to assess due to limitations in the amount and quality of previous elevation data sets. Other Svalbard regions have been precisely mapped by aerial photogrammetry, so the ICESat profiles from 2003-2008 can be compared with existing topographic maps and DEMs from 1965-1990. The mass balance for this period is estimated to -9.7 ± 0.6 Gt y-1 (or -0.36 ± 0.02 m w.e. y-1), excluding Austfonna. Repeat-track ICESat data are also analysed for the entire Svalbard yielding an average 2003-2008 mass balance of -4.3 ± 1.4 Gt y-1 (or -0.12 ± 0.04 m w.e. y-1) when tidewater front retreat is not accounted for. The most accurate elevation change estimates are obtained using all available ICESat data in a joint regression where surface slope and elevation change are estimated for rectangular planes that are fitted to the data along each track. The good performance of the plane method implies that it can also be used in other Arctic regions where accurate DEMs typically are not available

Quantification and interpretation of glacier elevation changes

Glaciers, ice caps and ice sheets constitute a large reservoir in the global hydrological cycle and provide a coupling between climate and sea-level. Observations of glacial change is important for constraining their contribution to sea-level fluctuations and to better understand the interactions between glaciers and climate. This thesis focuses on glacier observations through measurements of elevation change. The research in this thesis is oriented towards the methodological detection of elevation changes using remote sensing techniques. The quality of glacier elevation change measurements is dependent on controlling the potential errors and biases within the data. Therefore, one aspect is focused on a universal co-registration method for elevation products and further identification and correction of biases that remain, specifically in satellite stereo products. For glaciological studies, elevation changes require conversion into volume and mass changes. This is sometimes complicated when the data available is not spatially continuous and/or temporally consistent. Therefore, another aspect of this thesis explores methods for estimating regional glacier volume change. Specifically, Svalbard glacial contribution to sea-level has been estimated using regionalization techniques from scattered elevation measurements over roughly two time epochs. We observed that Svalbard glaciers over the past few decades have had a negative mass balance, contributing approximately 0.026 mm per year to the oceans. During the past few years, the sea-level contribution from Svalbard glaciers decreased slightly to 0.013 mm per year. Interpretations of elevation changes are convoluted by their dependence on climatic and dynamic forces operating on glacier systems. The last aspect of this thesis experiments with surface mass balance modelling for quantifying the climatic component of an elevation change. Combining this with observed elevation changes using theory of mass continuity can yield estimates of the calving flux of icebergs into the ocean. We observed on one particular fast flowing glacier in Svalbard that the average calving flux in the 1966-1990 epoch increased in the 1990-2007 epoch

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

The glaciers climate change initiative: Methods for creating glacier area, elevation change and velocity products

Glaciers and their changes through time are increasingly obtained from a wide range of satellite sensors. Due to the often remote location of glaciers in inaccessible and high-mountain terrain, satellite observations frequently provide the only available measurements. Furthermore, satellite data provide observations of glacier character- istics that are difficult to monitor using ground-based measurements, thus complementing the latter. In the Glaciers_cci project of the European Space Agency (ESA), three of these characteristics are investigated in detail: glacier area, elevation change and surface velocity. We use (a) data from optical sensors to derive glacier outlines, (b) digital elevation models from at least two points in time, (c) repeat altimetry for determining elevation changes, and (d) data from repeat optical and microwave sensors for calculating surface velocity. For the latter, the two sensor types provide complementary information in terms of spatio-temporal coverage. While (c) and (d) can be generated mostly automatically, (a) and (b) require the intervention of an analyst. Largely based on the results of various round robin experiments (multi-analyst benchmark studies) for each of the products, we suggest and describe the most suitable algorithms for product creation and provide recommendations concerning their practical implementation and the required post-processing. For some of the products (area, velocity) post-processing can influence product quality more than the main-processing algorithm

Error sources and guidelines for quality assessment of glacier area, elevation change, and velocity products derived from satellite data in the Glaciers_cci project

Satellite data provide a large range of information on glacier dynamics and changes. Results are often reported, provided and used without consideration of measurement accuracy (difference to a true value) and precision (variability of independent assessments). Whereas accuracy might be difficult to determine due to the limited availability of appropriate reference data and the complimentary nature of satellite measurements, precision can be obtained from a large range of measures with a variable effort for determination. This study provides a systematic overview on the factors influencing accuracy and precision of glacier area, elevation change (from altimetry and DEM differencing), and velocity products derived from satellite data, along with measures for calculating them. A tiered list of recommendations is provided (sorted for effort from Level 0 to 3) as a guide for analysts to apply what is possible given the datasets used and available to them. The more simple measures to describe product quality (Levels 0 and 1) can often easily be applied and should thus always be reported. Medium efforts (Level 2) require additional work but provide a more realistic assessment of product precision. Real accuracy assessment (Level 3) requires independent and coincidently acquired reference data with high accuracy. However, these are rarely available and their transformation into an unbiased source of information is challenging. This overview is based on the experiences and lessons learned in the ESA project Glaciers_cci rather than a review of the literature

Recent mass balance of the Purogangri Ice Cap, central Tibetan Plateau, by means of differential X-band SAR interferometry

Due to their remoteness, altitude and harsh climatic conditions, little is known about the glaciological parameters of ice caps on the Tibetan Plateau. This study presents a geodetic mass balance estimate of the Purogangri Ice Cap, Tibet's largest ice field between 2000 and 2012. We utilized data from the actual TerraSAR-X mission and its add-on for digital elevation measurements and compared it with elevation data from the Shuttle Radar Topography Mission. The employed data sets are ideal for this approach as both data sets were acquired at X-band at nearly the same time of the year and are available at a fine grid spacing. In order to derive surface elevation changes we employed two different methods. The first method is based on differential synthetic radar interferometry while the second method uses common DEM differencing. Both approaches revealed a slightly negative mass budget of −44 ± 15 and −38 ± 23 mm w.eq. a⁻¹ (millimeter water equivalent) respectively. A slightly negative trend of −0.15 ± 0.01 km² a⁻¹ in glacier extent was found for the same time period employing a time series of Landsat data. Overall, our results show an almost balanced mass budget for the studied time period. Additionally, we detected one continuously advancing glacier tongue in the eastern part of the ice cap

Sensitivity of glacier volume change estimation to DEM void interpolation

Glacier mass balance has been estimated on individual glacier and regional scales using repeat digital elevation models (DEMs). DEMs often have gaps in coverage (“voids”), the properties of which depend on the nature of the sensor used and the surface being measured. The way that these voids are accounted for has a direct impact on the estimate of geodetic glacier mass balance, though a systematic comparison of different proposed methods has been heretofore lacking. In this study, we determine the impact and sensitivity of void interpolation methods on estimates of volume change. Using two spatially complete, high-resolution DEMs over southeast Alaska, USA, we artificially generate voids in one of the DEMs using correlation values derived from photogrammetric processing of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) scenes. We then compare 11 different void interpolation methods on a glacier-by-glacier and regional basis. We find that a few methods introduce biases of up to 20&thinsp;% in the regional results, while other methods give results very close (&lt;1&thinsp;% difference) to the true, non-voided volume change estimates. By comparing results from a few of the best-performing methods, an estimate of the uncertainty introduced by interpolating voids can be obtained. Finally, by increasing the number of voids, we show that with these best-performing methods, reliable estimates of glacier-wide volume change can be obtained, even with sparse DEM coverage.</p

Historical analysis and visualization of the retreat of Findelengletscher, Switzerland, 1859–2010

Since the end of the Little Ice Age around 1850, glaciers in Europe have strongly retreated. Thanks to early topographic surveys in Switzerland, accurate maps are available, which enable us to trace glacier changes back in time. The earliest map for all of Switzerland that is usable for a detailed analysis is the Dufour map from around 1850 with subsequent topographic maps on a ~ 20 year interval. Despite the large public and scientific interest in glacier changes through time, this historic dataset has not yet been fully utilized for topographic change assessment or visualization of historic glacier extents. In this study, we use eleven historical topographic maps and more recent digital datasets for the region of Zermatt to analyze geometric changes (length, area and volume) of Findelengletscher as well as for creating animations of glacier evolution through time for use in public communication. All maps were georeferenced, the contour lines digitized, and digital elevation models (DEMs) created and co-registered. Additional digital data like the SRTM X-band DEM and high resolution laser scanning data were used to extend the analysis until 2010. Moreover, one independent DEM from aerial photogrammetry was used for comparison. During the period 1859–2010, Findelengletscher lost 3.5 km of its length (6.9 km in 2010), 4.42 ± 0.13 km² of its area (15.05 ± 0.45 km² in 2010) and 1.32 ± 0.52 km³ of its volume. The average rate of thickness loss is 0.45 ± 0.042 m yr− 1 for the 151 years period. Four periods with high thickness change from − 0.56 m ± 0.28 yr− 1 (1859–1881), − 0.40 ± 0.08 m yr− 1 (1937–1965), − 0.90 ± 0.31 m yr− 1 (1995–2000) and − 1.18 ± 0.02 m yr− 1 (2000–2005) have been identified. Small positive thickness changes were found for the periods 1890–1909 (+ 0.09 ± 0.46 m yr − 1) and 1988–1995 (+ 0.05 ± 0.24 m yr− 1). During its retreat with intermittent periods of advance, the glacier separated into three parts. The above changes are demonstrated through an animation (available from the supplementary material), which has been created to inform the general public