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

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

Satellite altimeter remote sensing of ice caps

This thesis investigates the use of satellite altimetry techniques for measuring surface elevation changes of ice caps. Two satellite altimeters, Radar Altimeter 2 (RA-2) and Geoscience Laser Altimeter System (GLAS) are used to assess the surface elevation changes of three Arctic ice caps. This is the first time the RA-2 has been used to assess the elevation changes of ice caps - targets much smaller than the ice sheets which are the instrument’s primary land ice targets. Algorithms for the retrieval of elevation change rates over ice caps using data acquired by RA-2 and GLAS are presented. These algorithms form a part of a European Space Agency (ESA) glacier monitoring system GlobGlacier. A comparison of GLAS elevation data to those acquired by the RA-2 shows agreement between the two instruments. Surface elevation change rate estimates based on RA-2 are given for three ice caps: Devon Ice Cap in Arctic Canada (−0.09 ± 0.29 m/a), Flade Isblink in Greenland (0.03 ± 0.03 m/a) and Austfonna on Svalbard (0.33 ± 0.08 m/a). Based on RA-2 and GLAS measurements it is shown that the areas of Flade Isblink below the late summer snow line have been thinning whereas the areas above the late summer snow line have been thickening. Also GLAS observed dynamic thickening rates of more than 3 m/a are presented. On Flade Isblink and Austfonna RA-2 measurements are compared to surface mass balance (SMB) estimates from a regional atmospheric climate model RACMO2. The comparison shows that SMB is the driver of interannual surface elevation changes at Austfonna. In contrast the comparison reveals areas on Flade Isblink where ice dynamics have an important effect on the surface elevation. Furthermore, RACMO2 estimates of surface mass budget at Austfonna before the satellite altimeter era are presented. This thesis shows that both traditional radar and laser satellite altimetry can be used to quantify the response of ice caps to the changing climate. Direct altimeter measurements of surface elevation and, in consequence volume change of ice caps, can be used to improve their mass budget estimates

A Reconciled Estimation of the State of Cryospheric Components in the Southern Andes and California Using Geospatial Techniques

Glaciers are the essential source of fresh water not only to human sustenance, but it is also vital for all lifeforms on earth. Glaciers are also key components in understanding rapid changes in climate. This makes understanding of glacier mass, extent, and overall state essential. In this dissertation, the objective was to analyze the state of snow and ice masses in the mid (California) and low latitude (Chile/Argentina) western American regions using geospatial technology. This study also analyzed the effects of anomalies in snow mass on the regional agricultural practices in California’s Central Valley. In the Southern Andes, the digital elevation models from Shuttle Radar Topographic Mission (SRTM) (the year 2000) were compared with the elevation footprints from the Geoscience Laser Altimeter System (GLAS) campaign for the years 2004 through 2008. Generally, in all sub-regions, the elevation values were lower than the elevation for the year 2000, which demarcates continuous recession of ice mass in the Andean region. Also, this study quantified snow cover extent and mass balance variation in the Sierra Nevada and Mt. Shasta regions in California. To unearth anomalies in snow mass, study used digital elevation models generated from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) between the year 2000 and 2015. A remarkable reduction in snow cover extent of about 80% was observed in the studied watersheds of California. Lastly, the impacts of snow mass anomalies on the total water storage (TWS) and agriculture land cover in the California’s Central Valley were quantified and geo-visualized. The study noticed the change in the land cover area of about 20% (6993 sq.km) due to the alteration of Agriculture land to impervious land covers. Most of the change in the agriculture land cover of about 4402 sq.km occurred in the San Joaquin and Tulare Basins of southern Central Valley region. This dissertation concludes that the increased temperature in the Andes and California has adversely impacted Cryosphere components in the region in the past decade. Besides, it provides valuable insights into the changing state of cryosphere components and highlights impacts of anomalies in TWS on a billion-dollar agricultural industry

Recent glacier and lake changes in High Mountain Asia and their relation to precipitation changes

We present an updated, spatially resolved estimate of 2003–2008 glacier surface elevation changes for the entire region of High Mountain Asia (HMA) from ICESat laser altimetry data. The results reveal a diverse pattern that is caused by spatially greatly varying glacier sensitivity, in particular to precipitation availability and changes. We introduce a spatially resolved zonation where ICESat samples are grouped into units of similar glacier behaviour, glacier type and topographic settings. In several regions, our new zonation reveals local differences and anomalies that have not been described previously. Glaciers in the Eastern Pamirs, Kunlun Shan and central TP were thickening by 0.1–0.7 m a−1, and the thickening anomaly has a crisp boundary in the Eastern Pamirs that continues just north of the central Karakoram. Glaciers in the south and east of the TP were thinning, with increasing rates towards southeast. We attribute the glacier thickening signal to a stepwise increase in precipitation around ∼1997–2000 on the Tibetan Plateau (TP). The precipitation change is reflected by growth of endorheic lakes in particular in the northern and eastern TP. We estimate lake volume changes through a combination of repeat lake extents from Landsat data and shoreline elevations from ICESat and the Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) for over 1300 lakes. The rise in water volume contained in the lakes corresponds to 4–25 mm a−1, when distributed over entire catchments, for the areas where we see glacier thickening. The precipitation increase is also visible in sparse in situ measurements and MERRA-2 climate reanalysis data but less visible in ERA-Interim reanalysis data. Taking into account evaporation loss, the difference between average annual precipitation during the 1990s and 2000s suggested by these datasets is 34–100 mm a−1, depending on region, which can fully explain both lake growth and glacier thickening (Kunlun Shan) or glacier geometry changes such as thinning tongues while upper glacier areas were thickening or stable (eastern TP). The precipitation increase reflected in these glacier changes possibly extended to the northern slopes of the Tarim Basin, where glaciers were nearly in balance in 2003–2008. Along the entire Himalaya, glaciers on the first orographic ridge, which are exposed to abundant precipitation, were thinning less than glaciers in the dryer climate of the inner ranges. Thinning rates in the Tien Shan vary spatially but are rather stronger than in other parts of HMA

High-temporal-resolution water level and storage change data sets for lakes on the Tibetan Plateau during 2000–2017 using multiple altimetric missions and Landsat-derived lake shoreline positions

The Tibetan Plateau (TP), known as Asia's water tower, is quite sensitive to climate change, which is reflected by changes in hydrologic state variables such as lake water storage. Given the extremely limited ground observations on the TP due to the harsh environment and complex terrain, we exploited multiple altimetric missions and Landsat satellite data to create high-temporal-resolution lake water level and storage change time series at weekly to monthly timescales for 52 large lakes (50 lakes larger than 150 km2 and 2 lakes larger than 100 km2) on the TP during 2000-2017. The data sets are available online at https://doi.org/10.1594/PANGAEA.898411 (Li et al., 2019). With Landsat archives and altimetry data, we developed water levels from lake shoreline positions (i.e., Landsat-derived water levels) that cover the study period and serve as an ideal reference for merging multisource lake water levels with systematic biases being removed. To validate the Landsat-derived water levels, field experiments were carried out in two typical lakes, and theoretical uncertainty analysis was performed based on high-resolution optical images (0.8 m) as well. The RMSE of the Landsat-derived water levels is 0.11m compared with the in situ measurements, consistent with the magnitude from theoretical analysis (0.1-0.2 m). The accuracy of the Landsat-derived water levels that can be derived in relatively small lakes is comparable with most altimetry data. The resulting merged Landsat-derived and altimetric lake water levels can provide accurate information on multiyear and short-term monitoring of lake water levels and storage changes on the TP, and critical information on lake overflow flood monitoring and prediction as the expansion of some TP lakes becomes a serious threat to surrounding residents and infrastructure

Use of remote sensing to assess supra-glacial lake depths on the Greenland Ice Sheet.

The influence that supra-glacial lakes have had in the recent mass loss at the margins of the Greenland ice sheet has been widely studied. Lakes can drain to the lase of a glacier, lubricating the bed, and enhancing acceleration of the glacier and hence ice thinning. Recent studies suggested that melt extent is not directly linked to the dynamic loss but it has been broken to be linked to peak summer speed ups of the ice sheet front. Large volumes of water are necessary to propagate cracks to the glacial bed via hydrofractures. Hydrological models showed that lakes above a critical volume can supply the necessary water for this process, so the ability to measure water depth in lakes remotely is important to study these processes. The aim of this thesis was to test the current models used for water depth calculations based on the optical properties of water. An optimisation model to estimate water depths was developed. Atmospherically-corrected data from ASTER and MODIS were used as an input to the water reflectance model. As a reference dataset, ICESat measurements were used to obtain lake geometries over empty lakes. Differences between modelled and reference depths are used in a minimisation model to obtain parameters for the water-reflectance model, yielding optimised lake depth estimates. The key contribution of this research was the development of a Monte Carlo simulation. This method allows the quantification of uncertainties in water depth and hence water volume, for the first time. This robust analysis provided better understanding of the sensitivity of the model to the input parameters. There is scope to improve current models of depth estimations if more extensive held observations are done

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

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings