Performance Assessment of TanDEM-X DEM for Mountain Glacier Elevation Change Detection

Evaluating the impact of the Central Chile Mega Drought on debris cover, broadband albedo, and surface drainage system of a Dry Andes glacier

In recent years, Chile has experienced an extraordinary drought that has had significant impacts on both the livelihoods of people and the environment, including the Andean glaciers. This study focuses on analyzing the surface processes of Universidad Glacier, a benchmark glacier for the Dry Andes. Multiple remote sensing datasets are used alongside a novel spectral index designed for mapping of rock material located on the glacier's surface. Our findings highlight the precarious state of the glacier, which serves as a crucial water source for the region. The glacier exhibits locally varied debris accumulation and margin retreat. The most significant impacts are observed on the tongue and secondary accumulation cirques, with the latter at risk of disappearing. The debris cover on the tongue is expanding, reaching higher elevations, and is accompanied by glacier retreat, especially at higher altitudes. The equilibrium line is rapidly shifting upglacier, although the mid-season snow cover still frequently reaches the 2013 equilibrium line, even in 2020. Changes in stream density on the glacier tongue indicate an increased water supply in this area, likely due to enhanced melting of glacial ice. These observed processes align well with meteorological data obtained from reanalysis products. The behavior of dust and debris is influenced by precipitation amount, while the rate of retreat is linked to air temperature.

Remote Sensing Application in Mountainous Environments: A Bibliographic Analysis

Advancement in remote sensing platforms, sensors, and technology has significantly improved the assessment of hard-to-access areas, such as mountains. Despite these improvements, Africa lags in terms of research work published. This is of great concern as the continent needs more research to achieve sustainable development. Therefore, this study applied a bibliometric analysis of the annual production of publications on the application of remote sensing methods in mountainous environments. In total, 3849 original articles between 1973 and 2021 were used, and the results indicate a steady growth in publications from 2004 (n = 26) to 2021 (n = 504). Considering the source journals, Remote Sensing was the top-ranked, with 453 total publications. The University of the Chinese Academy of Sciences was the highest-ranking affiliation, with 217 articles, and China produced the highest number of publications (n = 217). Keywords used between 1973 and 1997, such as "Canada", "alps", and "GIS", metamorphosed into "remote sensing" between 1998 and 2021. This metamorphosis indicates a change in the areas of interest and an increase in the application of remote sensing methods. Most studies were conducted in the Global North countries, and a few were published in low-impact journals within the African continent. This study can help researchers and scholars better understand the progress and intellectual structure of the field and future research directions in the application of remote sensing methods in mountainous environments.

Co-Registration Methods and Error Analysis for Four Decades (1979–2018) of Glacier Elevation Changes in the Southern Patagonian Icefield

The main goal of this paper is to compare two co-registration methods for geodetic mass balance (GMB) calculation in 28 glaciers making up the Upper Santa Cruz River basin, Southern Patagonian Icefield (SPI), from 1979 to 2018. For this purpose, geospatial data have been used as primary sources: Hexagon KH-9, ASTER, and LANDSAT optical images; SRTM digital radar elevation model; and ICESat elevation profiles. After the analyses, the two co-registration methods, namely M1, based on horizontal displacements and 3D shift vectors, and M2, based on three-dimensional transformations, turned out to be similar. The errors in the GMB were analyzed through a k index that considers, among other variables, the error in elevation change by testing four interpolation methods for filling gaps. We found that, in 63% of the cases, the relative error in elevation change contributes 90% or more to k index. The GMB throughout our study area reported that a loss value of −1.44 ± 0.15 m w. e. a−1 (−3.0 Gt a−1) and an ice thinning median of −1.38 ± 0.11 m a−1 occurred within the study period. The glaciers that showed the most negative GMB values were Upsala, with an annual elevation change median of −2.07 ± 0.18 m w. e. a−1, and Ameghino, with −2.31 ± 0.22 m w. e. a−1.

Geodetic Mass Balance of Haxilegen Glacier No. 51, Eastern Tien Shan, from 1964 to 2018

The eastern Tien Shan hosts substantial mid-latitude glaciers, but in situ glacier mass balance records are extremely sparse. Haxilegen Glacier No. 51 (eastern Tien Shan, China) is one of the very few well-measured glaciers, and comprehensive glaciological measurements were implemented from 1999 to 2011 and re-established in 2017. Mass balance of Haxilegen Glacier No. 51 (1999–2015) has recently been reported, but the mass balance record has not extended to the period before 1999. Here, we used a 1:50,000-scale topographic map and long-range terrestrial laser scanning (TLS) data to calculate the area, volume, and mass changes for Haxilegen Glacier No. 51 from 1964 to 2018. Haxilegen Glacier No. 51 lost 0.34 km2 (at a rate of 0.006 km2 a−1 or 0.42% a−1) of its area during the period 1964–2018. The glacier experienced clearly negative surface elevation changes and geodetic mass balance. Thinning occurred almost across the entire glacier surface, with a mean value of −0.43 ± 0.12 m a−1. The calculated average geodetic mass balance was −0.36 ± 0.12 m w.e. a−1. Without considering the error bounds of mass balance estimates, glacier mass loss over the past 50 years was in line with the observed and modeled mass balance (−0.37 ± 0.22 m w.e. a−1) that was published for short time intervals since 1999 but was slightly less negative than glacier mass loss in the entire eastern Tien Shan. Our results indicate that Riegl VZ®-6000 TLS can be widely used for mass balance measurements of unmonitored individual glaciers.

Applicability of Data Acquisition Characteristics to the Identification of Local Artefacts in Global Digital Elevation Models: Comparison of the Copernicus and TanDEM-X DEMs

Several global digital elevation models (DEMs) have been developed in the last two decades. The most recent addition to the family of global DEMs is the TanDEM-X DEM. The original version of the TanDEM-X DEM is, however, a nonedited product (i.e., it contains local artefacts such as voids, spikes, and holes). Therefore, subsequent identification of local artefacts and their editing is necessary. In this study, we evaluated the accuracy of the original TanDEM-X DEM and its improved edited version, the Copernicus DEM, in three major European mountain ranges (the Alps, the Carpathians, and the Pyrenees) using a digital surface model derived from airborne laser scanning data as a reference. In addition, to evaluate the applicability of data acquisition characteristics (coverage map, consistency mask, and height error map) and terrain characteristics (slope, aspect, altitude) to the localization of problematic sites, we modeled their associations with the TanDEM-X DEM error. We revealed local occurrences of large errors in the TanDEM-X DEM that were typically found on steep ridges or in canyons, which were largely corrected in the Copernicus DEM. The editing procedure used for the Copernicus DEM construction was evidently successful as the RMSE for the TanDEM-X and Copernicus DEMs at the 90 m resolution improved from 45 m to 12 m, from 16 m to 6 m, and from 24 m to 9 m for the Alps, the Pyrenees, and the Carpathians, respectively. The Copernicus DEM at the 30 m resolution performed similarly well. The boosted regression trees showed that acquisition characteristics provided as auxiliary data are useful for locating problematic sites and explained 28–50% of deviance of the absolute vertical error. The absolute vertical error was strongly related to the height error map. Finally, up to 26% of cells in the Copernicus DEM were filled using DEMs from different time periods and, hence, users performing multitemporal analysis or requiring data from a specific time period in the mountain environment should be wary when using TanDEM-X and its derivations. We suggest that when filling problematic sites using alternative DEMs, more attention should be paid to the period of their collection to minimize the temporal displacement in the final products.

Use of TanDEM-X and SRTM-C Data for Detection of Deforestation Caused by Bark Beetle in Central European Mountains

The availability of global digital elevation models (DEMs) from multiple time points allows their combination for analysing vegetation changes. The combination of models (e.g., SRTM and TanDEM-X) can contain errors, which can, due to their synergistic effects, yield incorrect results. We used a high-resolution LiDAR-derived digital surface model (DSM) to evaluate the accuracy of canopy height estimates of the aforementioned global DEMs. In addition, we subtracted SRTM and TanDEM-X data at 90 and 30 m resolutions, respectively, to detect deforestation caused by bark beetle disturbance and evaluated the associations of their difference with terrain characteristics. The study areas covered three Central European mountain ranges and their surrounding areas: Bohemian Forest, Erzgebirge, and Giant Mountains. We found that vertical bias of SRTM and TanDEM-X, relative to the canopy height, is similar with negative values of up to −2.5 m and LE90s below 7.8 m in non-forest areas. In forests, the vertical bias of SRTM and TanDEM-X ranged from −0.5 to 4.1 m and LE90s from 7.2 to 11.0 m, respectively. The height differences between SRTM and TanDEM-X show moderate dependence on the slope and its orientation. LE90s for TDX-SRTM differences tended to be smaller for east-facing than for west-facing slopes, and varied, with aspect, by up to 1.5 m in non-forest areas and 3 m in forests, respectively. Finally, subtracting SRTM and NASA DEMs from TanDEM-X and Copernicus DEMs, respectively, successfully identified large areas of deforestation caused by hurricane Kyril in 2007 and a subsequent bark beetle disturbance in the Bohemian Forest. However, local errors in TanDEM-X, associated mainly with forest-covered west-facing slopes, resulted in erroneous identification of deforestation. Therefore, caution is needed when combining SRTM and TanDEM-X data in multitemporal studies in a mountain environment. Still, we can conclude that SRTM and TanDEM-X data represent suitable near global sources for the identification of deforestation in the period between the time points of their acquisition.

Quality Assessment of TanDEM-X DEMs, SRTM and ASTER GDEM on Selected Chinese Sites

Digital elevation models (DEMs) are the basic data of science and engineering technology research. SRTM and ASTER GDEM are currently widely used global DEMs, and TanDEM-X DEM, released in 2016, has attracted users’ attention due to its unprecedented accuracy. These global datasets are often used for local applications and the quality of DEMs affects the results of applications. Many researchers have assessed and compared the quality of global DEMs on a local scale. To provide some additional insights on quality assessment of 12- and 30-m resolution TanDEM-X DEMs, 30-m resolution ASTER GDEM and 30-m resolution SRTM, this study assessed differences’ performance in relation to not only geographical features but also the ways in which DEMs have been created on selected Chinese sites, taking ICESat/GLAS points with 14-cm absolute vertical accuracy but size of 70-m diameter and 12-m resolution TanDEM-X DEM with less than 10-m absolute vertical accuracy as the reference data for comprehensive quality evaluation. When comparing the three 30-m DEMs with the reference DEM, an improved Least Z-Difference (LZD) method was applied for co-registration between models, and Quantile–Quantile (Q-Q) plot was used to identify if the DEM errors follow a normal distribution to help choose proper statistical indicators accordingly. The results show that: (1) TanDEM-X DEMs have the best overall quality, followed by SRTM. ASTER GDEM has the worst quality. The 12-m TanDEM-X DEM has significant advantages in describing terrain details. (2) The quality of DEM has a strong relationship with slope, aspect and land cover. However, the relationship between aspect and vertical quality weakens after data co-registration. The quality of DEMs gets higher with the increasing number of images used in the fusion process. The quality in where slopes opposite to the radar beam is the worst for SRTM, which could provide a new perspective for quality assessment of SRTM and other DEMs whose incidence angle files are available. (3) Systematic deviations can reduce the vertical quality of DEM. The differences have non-normal distribution even after co-registration. For researchers who want to know the quality of a DEM in order to use it in further applications, they should pay more attention to the terrain factors and land cover in their study areas and the ways in which the DEM has been created.

Evaluation of the Impact of Climate Change on Runoff Generation in an Andean Glacier Watershed

Excluding Antarctica and Greenland, 3.8% of the world’s glacier area is concentrated in Chile. The country has been strongly affected by the mega drought, which affects the south-central area and has produced an increase in dependence on water resources from snow and glacier melting in dry periods. Recent climate change has led to an elevation of the zero-degree isotherm, a decrease in solid-state precipitation amounts and an accelerated loss of glacier and snow storage in the Chilean Andes. This situation calls for a better understanding of future water discharge in Andean headwater catchments in order to improve water resources management in glacier-fed populated areas. The present study uses hydrological modeling to characterize the hydrological processes occurring in a glacio-nival watershed of the central Andes and to examine the impact of different climate change scenarios on discharge. The study site is the upper sub-watershed of the Tinguiririca River (area: 141 km2), of which nearly 20% is covered by Universidad Glacier. The semi-distributed Snowmelt Runoff Model + Glacier (SRM+G) was forced with local meteorological data to simulate catchment runoff. The model was calibrated on even years and validated on odd years during the 2008–2014 period and found to correctly reproduce daily runoff. The model was then forced with downscaled ensemble projected precipitation and temperature series under the RCP 4.5 and RCP 8.5 scenarios, and the glacier adjusted using a volume-area scaling relationship. The results obtained for 2050 indicate a decrease in mean annual discharge (MAD) of 18.1% for the lowest emission scenario and 43.3% for the most pessimistic emission scenario, while for 2100 the MAD decreases by 31.4 and 54.2%, respectively, for each emission scenario. Results show that decreasing precipitation lead to reduced rainfall and snowmelt contributions to discharge. Glacier melt thus partly buffers the drying climate trend, but our results show that the peak water occurs near 2040, after which glacier depletion leads to reducing discharge, threatening the long-term water resource availability in this region.

Impact of UAV Surveying Parameters on Mixed Urban Landuse Surface Modelling

Unmanned Aerial Vehicles (UAVs) as a surveying tool are mainly characterized by a large amount of data and high computational cost. This research investigates the use of a small amount of data with less computational cost for more accurate three-dimensional (3D) photogrammetric products by manipulating UAV surveying parameters such as flight lines pattern and image overlap percentages. Sixteen photogrammetric projects with perpendicular flight plans and a variation of 55% to 85% side and forward overlap were processed in Pix4DMapper. For UAV data georeferencing and accuracy assessment, 10 Ground Control Points (GCPs) and 18 Check Points (CPs) were used. Comparative analysis was done by incorporating the median of tie points, the number of 3D point cloud, horizontal/vertical Root Mean Square Error (RMSE), and large-scale topographic variations. The results show that an increased forward overlap also increases the median of the tie points, and an increase in both side and forward overlap results in the increased number of point clouds. The horizontal accuracy of 16 projects varies from ±0.13m to ±0.17m whereas the vertical accuracy varies from ± 0.09 m to ± 0.32 m. However, the lowest vertical RMSE value was not for highest overlap percentage. The tradeoff among UAV surveying parameters can result in high accuracy products with less computational cost.

Detailed Lacustrine Calving Iceberg Inventory from Very High Resolution Optical Imagery and Object-Based Image Analysis

In the field of iceberg and glacier calving studies, it is important to collect comprehensive datasets of populations of icebergs. Particularly, calving of lake-terminating glaciers has been understudied. The aim of this work is to present an object-based method of iceberg detection and to create an inventory of icebergs located in a proglacial lagoon of San Quintín glacier, Northern Patagonia Icefield, Chile. This dataset is created using high-resolution WorldView-2 imagery and a derived DEM. We use it to briefly discuss the iceberg size distribution and area–volume scaling. Segmentation of the multispectral imagery produced a map of objects, which were classified with use of Random Forest supervised classification algorithm. An intermediate classification product was corrected with a ruleset exploiting contextual information and a watershed algorithm that was used to divide multiple touching icebergs into separate objects. Common theoretical heavy-tail statistical distributions were tested as descriptors of iceberg area and volume distributions. Power law models were proposed for the area–volume relationship. The proposed method performed well over the open lake detecting correctly icebergs in all size bands except 5–15 m2 where many icebergs were missed. A section of the lagoon with ice melange was not reliably mapped due to uniformity of the area in the imagery and DEM. The precision of the DEM limited the scaling effort to icebergs taller than 1.7 m and larger than 99 m2, despite the inventory containing icebergs as small as 4 m2. The work demonstrates viability of object-based analysis for lacustrine iceberg detection and shows that the statistical properties of iceberg population at San Quintín glacier match those of populations produced by tidewater glaciers.

A Relief Dependent Evaluation of Digital Elevation Models on Different Scales for Northern Chile

Many geoscientific computations are directly influenced by the resolution and accuracy of digital elevation models (DEMs). Therefore, knowledge about the accuracy of DEMs is essential to avoid misleading results. In this study, a comprehensive evaluation of the vertical accuracy of globally available DEMs from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Shuttle Radar Topography Mission (SRTM), Advanced Land Observing Satellite (ALOS) World 3D and TanDEM-X WorldDEM™ was conducted for a large region in Northern Chile. Additionally, several very high-resolution DEM datasets were derived from Satellite Pour l’Observation de la Terre (SPOT) 6/7 and Pléiades stereo satellite imagery for smaller areas. All datasets were evaluated with three reference datasets, namely elevation points from both Ice, Cloud, and land Elevation (ICESat) satellites, as well as very accurate high-resolution elevation data derived by unmanned aerial vehicle (UAV)-based photogrammetry and terrestrial laser scanning (TLS). The accuracy was also evaluated with regard to the existing relief by relating the accuracy results to slope, terrain ruggedness index (TRI) and topographic position index (TPI). For all datasets with global availability, the highest overall accuracies are reached by TanDEM-X WorldDEM™ and the lowest by ASTER Global DEM (GDEM). On the local scale, Pléiades DEMs showed a slightly higher accuracy as SPOT imagery. Generally, accuracy highly depends on topography and the error is rising up to four times for high resolution DEMs and up to eight times for low-resolution DEMs in steeply sloped terrain compared to flat landscapes.

Combined Use of Terrestrial Laser Scanning and UAV Photogrammetry in Mapping Alpine Terrain

Airborne and terrestrial laser scanning and close-range photogrammetry are frequently used for very high-resolution mapping of land surface. These techniques require a good strategy of mapping to provide full visibility of all areas otherwise the resulting data will contain areas with no data (data shadows). Especially, deglaciated rugged alpine terrain with abundant large boulders, vertical rock faces and polished roche-moutones surfaces complicated by poor accessibility for terrestrial mapping are still a challenge. In this paper, we present a novel methodological approach based on a combined use of terrestrial laser scanning (TLS) and close-range photogrammetry from an unmanned aerial vehicle (UAV) for generating a high-resolution point cloud and digital elevation model (DEM) of a complex alpine terrain. The approach is demonstrated using a small study area in the upper part of a deglaciated valley in the Tatry Mountains, Slovakia. The more accurate TLS point cloud was supplemented by the UAV point cloud in areas with insufficient TLS data coverage. The accuracy of the iterative closest point adjustment of the UAV and TLS point clouds was in the order of several centimeters but standard deviation of the mutual orientation of TLS scans was in the order of millimeters. The generated high-resolution DEM was compared to SRTM DEM, TanDEM-X and national DMR3 DEM products confirming an excellent applicability in a wide range of geomorphologic applications.

Glacier Facies Mapping Using a Machine-Learning Algorithm: The Parlung Zangbo Basin Case Study

Glaciers in the Tibetan Plateau are an important indicator of climate change. Automatic glacier facies mapping utilizing remote sensing data is challenging due to the spectral similarity of supraglacial debris and the adjacent bedrock. Most of the available glacier datasets do not provide the boundary of clean ice and debris-covered glacier facies, while debris-covered glacier facies play a key role in mass balance research. The aim of this study was to develop an automatic algorithm to distinguish ice cover types based on multi-temporal satellite data, and the algorithm was implemented in a subregion of the Parlung Zangbo basin in the southeastern Tibetan Plateau. The classification method was built upon an automated machine learning approach: Random Forest in combination with the analysis of topographic and textural features based on Landsat-8 imagery and multiple digital elevation model (DEM) data. Very high spatial resolution Gao Fen-1 (GF-1) Panchromatic and Multi-Spectral (PMS) imagery was used to select training samples and validate the classification results. In this study, all of the land cover types were classified with overall good performance using the proposed method. The results indicated that fully debris-covered glaciers accounted for approximately 20.7% of the total glacier area in this region and were mainly distributed at elevations between 4600 m and 4800 m above sea level (a.s.l.). Additionally, an analysis of the results clearly revealed that the proportion of small size glaciers (<1 km2) were 88.3% distributed at lower elevations compared to larger size glaciers (≥1 km2). In addition, the majority of glaciers (both in terms of glacier number and area) were characterized by a mean slope ranging between 20° and 30°, and 42.1% of glaciers had a northeast and north orientation in the Parlung Zangbo basin.