TanDEM-X Digital Elevation Model Research Articles

Generation of Dense and High-Precision Digital Elevation Model Using Low-Cost Unmanned Aerial Vehicle and Space-Borne TanDEM-X to Measure Exposed Area Change Due to Tidal Invasion

Tidal flats are internationally protected areas because of their environmental, geological, and economic value. Nevertheless, these are areas where casualties are common among visitors. To ensure the safety of visitors, efforts must be made to understand the geomorphologic characteristics of tidal flats. It is necessary to use three-dimensional topographic data to clearly determine when an exposed tidal flat rapidly disappears as a result of periodic tidal fluctuations. Digital elevation models (DEMs) created using unmanned aerial vehicle (UAV) imaging and commercial photogrammetric software to measure tidal flats are produced using ground control points (GCPs). It is difficult to conduct field surveys and readings of image coordinates that correspond to these because tidal flat areas are difficult to access. The distribution of GCPs affects the accuracy of a DEM because the entire test region cannot be covered during a survey. Without GCPs, the DEM generated with the UAV images and photogrammetric software using the structure from motion (SfM) technique probably has nonlinear distortion such as a bowl shape. This paper proposes a practical method of least-squares three-dimensional (3D) surface matching with a polynomial model (LS3D-PM) to correct a distorted DEM. For this LS3D-PM matching, a global TanDEM-X DEM was used instead of a field survey to collect the GCPs in a tidal flat. Practical testing was conducted using the proposed method. The root mean square error with respect to the height decreased from approximately 1.0 m to 0.1 m, and the bowl-effect error was eliminated. With the DEM corrected using the proposed method and the classified orthoimages, the exposed tidal flat changes due to seawater invasion were estimated at different times during flood tide. This method has non-trivial technical values for quick and relatively precise DEM generations at a low cost compared to that of other techniques, particularly in emergency situations such as tsunamis, landslides, and earthquakes.

Spatiotemporal Dynamics and Geodetic Mass Changes of Glaciers With Varying Debris Cover in the Pangong Region of Trans-Himalayan Ladakh, India Between 1990 and 2019

Glaciers across the Himalayan arc are showing varying signs of recession. Glaciers in the eastern and western parts of the Himalayan arc are retreating more rapidly as compared to other regions. This differential retreat is often attributed to climatic, topographic, and geologic influences. The glaciers in the Trans-Himalayan region of Ladakh are believed to be relatively stable as compared to other parts of the western Himalaya. The present study ascertained the area changes and frontal retreat of 87 glaciers in the Pangong Region between 1990 and 2019 using satellite data. The geodetic mass changes were also assessed using SRTM and TanDEM-X digital elevation models of 2000 and 2012 respectively. Besides, the glacier outlines were delineated manually and compared with existing regional and global glacier inventories that are available over the region. The GlabTop model was used to simulate the glacier-bed overdeepenings of four glaciers that are associated with a proglacial lake. The study also analyzed the impact of topographic influences and varying debris cover on glacier recession. This analysis indicated deglaciation of 6.7 ± 0.1% (0.23% a−1) from 1990 to 2019 over the Pangong Region with clean-ice glaciers showing a higher retreat (8.4 ± 0.28%) compared to the debris-covered glaciers (5.7 ± 0.14%). However, the overall recession is lower compared to other parts of northwestern Himalayas. The glacier recession showed a positive correlation with mean glacier slope (r = 0.3) and debris cover (r = 0.1) with bigger size glaciers having retreated at a lesser pace compared to smaller ones. This underpins the need for in-situ data about debris thickness to precisely ascertain the role of debris on glacier recession in the Trans-Himalayan Ladakh where debris thickness data is absent. The mean glacier elevation did not indicate any influence on glacier recession. From 2000 to 12, the glaciers lost an ice mass amounting to 0.33 ± 0.05 m we. per year. The formation of four new proglacial lakes, although small (<6 ha), need to be monitored using remote sensing data while the infrastructure development activities should not be permitted given glacial lake outburst flood risk.

Parameterized Modeling and Calibration for Orbital Error in TanDEM-X Bistatic SAR Interferometry over Complex Terrain Areas

The TerraSAR-X add-on for Digital Elevation Measurements (TanDEM-X) bistatic system provides high-resolution and high-quality interferometric data for global topographic measurement. Since the twin TanDEM-X satellites fly in a close helix formation, they can acquire approximately simultaneous synthetic aperture radar (SAR) images, so that temporal decorrelation and atmospheric delay can be ignored. Consequently, the orbital error becomes the most significant error limiting high-resolution SAR interferometry (InSAR) applications, such as the high-precision digital elevation model (DEM) reconstruction, subway and highway deformation monitoring, landslide monitoring and sub-canopy topography inversion. For rugged mountainous areas, in particular, it is difficult to estimate and correct the orbital phase error in TanDEM-X bistatic InSAR. Based on the rigorous InSAR geometric relationship, the orbital phase error can be attributed to the baseline errors (BEs) after fixing the positions of the master SAR sensor and the targets on the ground surface. For the constraint of the targets at a study scene, the freely released TanDEM-X DEM can be used, due to its consistency with the TanDEM-X bistatic InSAR-measured height. As a result, a parameterized model for the orbital phase error estimation is proposed in this paper. In high-resolution and high-precision TanDEM-X bistatic InSAR processing, due to the limited precision of the navigation systems and the uneven baseline changes caused by the helix formation, the BEs are time-varying in most cases. The parameterized model is thus built and estimated along each range line. To validate the proposed method, two mountainous test sites located in China (i.e., Fuping in Shanxi province and Hetang in Hunan province) were selected. The obtained results show that the orbital phase errors of the bistatic interferograms over the two test sites are well estimated. Compared with the widely applied polynomial model, the residual phase corrected by the proposed method contains little undesirable topography-dependent phase error, and avoids unexpected height errors ranging about from −6 m to 3 m for the Fuping test site and from −10 m to 8 m for the Hetang test site. Furthermore, some fine details, such as ridges and valleys, can be clearly identified after the correction. In addition, the two components of the orbital phase error, i.e., the residual flat-earth phase error and the topographic phase error caused by orbital error, are separated and quantified based on the parameterized expression. These demonstrate that the proposed method can be used to accurately estimate and mitigate the orbital phase error in TanDEM-X bistatic InSAR data, which increases the feasibility of reconstructing high-resolution and high-precision DEM. The rigorous geometric constraint, the refinement of the initial baseline parameters, and the assessment for height errors based on the estimated BEs are investigated in the discussion section of this paper.

The Global Water Body Layer from TanDEM-X Interferometric SAR Data

The interferometric synthetic aperture radar (InSAR) data set, acquired by the TanDEM-X (TerraSAR-X add-on for Digital Elevation Measurement) mission (TDM), represents a unique data source to derive geo-information products at a global scale. The complete Earth’s landmasses have been surveyed at least twice during the mission bistatic operation, which started at the end of 2010. Examples of the delivered global products are the TanDEM-X digital elevation model (DEM) (at a final independent posting of 12 m × 12 m) or the TanDEM-X global Forest/Non-Forest (FNF) map. The need for a reliable water product from TanDEM-X data was dictated by the limited accuracy and difficulty of use of the TDX Water Indication Mask (WAM), delivered as by-product of the global DEM, which jeopardizes its use for scientific applications, as well. Similarly as it has been done for the generation of the FNF map; in this work, we utilize the global data set of TanDEM-X quicklook images at 50 m × 50 m resolution, acquired between 2011 and 2016, to derive a new global water body layer (WBL), covering a range from −60∘ to +90∘ latitudes. The bistatic interferometric coherence is used as the primary input feature for performing water detection. We classify water surfaces in single TanDEM-X images, by considering the system’s geometric configuration and exploiting a watershed-based segmentation algorithm. Subsequently, single overlapping acquisitions are mosaicked together in a two-step logically weighting process to derive the global TDM WBL product, which comprises a binary averaged water/non-water layer as well as a permanent/temporary water indication layer. The accuracy of the new TDM WBL has been assessed over Europe, through a comparison with the Copernicus water and wetness layer, provided by the European Space Agency (ESA), at a 20 m × 20 m resolution. The F-score ranges from 83%, when considering all geocells (of 1∘ latitudes × 1∘ longitudes) over Europe, up to 93%, when considering only the geocells with a water content higher than 1%. At global scale, the quality of the product has been evaluated, by intercomparison, with other existing global water maps, resulting in an overall agreement that often exceeds 85% (F-score) when the content in the geocell is higher than 1%. The global TDM WBL presented in this study will be made available to the scientific community for free download and usage.

Effusion rate estimation based on solidified lava flows: Implications for volcanic hazard assessment in the Negros de Aras monogenetic volcanic field, northern Chile

Effusion rate is the instantaneous lava flow output by a vent. It is one of the most important factors that govern the emplacement and dynamics of lava flows and can be determined by direct measurements or estimated through modeling. The Negros de Aras monogenetic volcanic field, located in the Central Volcanic Zone of the Andes, is the largest volcanic field in northern Chile. It is situated north of the active Socompa volcano and south of the Salar de Atacama basin, consisting of scoria cones and lava flows. Here, we estimate the effusion rates of three overlapping and representative solidified lava flows from this volcanic field applying an iterative procedure using the Q-LavHA GIS plugin, where the effusion rate is treated as an unknown parameter and the solution is reached by comparing the simulated results with the real extension of the flows. The pre-eruptive surface was reconstructed using the 12 m resolution TanDEM-X digital elevation model. Other required input parameters for the modeling (e.g., channel ratio, viscosity, temperature) were estimated through geomorphometric, petrographic, and geochemical analysis. The estimated effusion rates vary between 14 and 113 m3/s, comparable with rates measured elsewhere. The method was validated by applying it to two directly measured (Nyamuragira 2006; Mount Etna 2001) and one calculated effusion rate (Lentiscal 2450 BP, Canary Islands). The estimated effusion rates gave an average underestimation of 10%. Considering this percentage of adjustment, the effusion rates for the three studied lava flows from the Negros de Aras monogenetic volcanic field vary from 15 m3/s to 124 m3/s. Our results provide valuable information for mitigating volcanic risk in the event of future lava flow eruptions, which could affect nearby villages and mining operations.

TanDEM-X PolarDEM 90 m of Antarctica: generation and error characterization

Abstract. We present the generation and validation of an updated version of the TanDEM-X digital elevation model (DEM) of Antarctica: the TanDEM-X PolarDEM 90 m of Antarctica. Improvements compared to the global TanDEM-X DEM version comprise filling gaps with newer bistatic synthetic aperture radar (SAR) acquisitions of the TerraSAR-X and TanDEM-X satellites, interpolation of smaller voids, smoothing of noisy areas, and replacement of frozen or open sea areas with geoid undulations. For the latter, a new semi-automatic editing approach allowed for the delineation of the coastline from DEM and amplitude data. Finally, the DEM was transformed into the cartographic Antarctic Polar Stereographic projection with a homogeneous metric spacing in northing and easting of 90 m. As X-band SAR penetrates the snow and ice pack by several meters, a new concept for absolute height adjustment was set up that relies on areas with stable penetration conditions and on ICESat (Ice, Cloud, and land Elevation Satellite) elevations. After DEM generation and editing, a sophisticated height error characterization of the whole Antarctic continent with ICESat data was carried out, and a validation over blue ice achieved a mean vertical height error of just −0.3 m ± 2.5 m standard deviation. The filled and edited Antarctic TanDEM-X PolarDEM 90 m is outstanding due to its accuracy, homogeneity, and coverage completeness. It is freely available for scientific purposes and provides a high-resolution data set as basis for polar research, such as ice velocity, mass balance estimation, or orthorectification.

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.

Comprehensive vertical accuracy analysis of freely available DEMs for different landscape types of the Rur catchment, Germany

This study evaluates the vertical accuracy of nine freely available digital elevation models (DEMs) with a medium spatial resolution for a typical mid-latitude flat to hilly landscape in the Rur catchment in Germany. All datasets were evaluated with two reference datasets, a 1 m DEM and a highly precise set of elevation control points. The accuracy was evaluated with regard to different slopes, landforms and land use. The results reveal in flat areas an average error of 2–4 m which increases about two to three times in undulated terrain. Areas with dense tree cover show a mean error of 6–10 m. The highest accuracies were achieved by the SRTM DEM, followed by other DEMs that mainly implemented SRTM data. TanDEM-X and Copernicus DEM showed ambivalent results showing a higher ratio of existing errors in the elevation product for vegetation and water areas that prevent a higher overall accuracy of these DEMs.

High-resolution topography of the Antarctic Peninsula combining the TanDEM-X DEM and Reference Elevation Model of Antarctica (REMA) mosaic

Abstract. The Antarctic Peninsula (AP) is one of the widely studied polar regions because of its sensitivity to climate change and potential contribution of its glaciers to global sea level rise. Precise digital elevation models (DEMs) at a high spatial resolution are much demanded for investigating the complex glacier system of the AP at fine scales. However, the two most recent circum-Antarctic DEMs, the 12 m TanDEM-X DEM (TDM DEM) from bistatic interferometric synthetic aperture radar (InSAR) data acquired between 2013 and 2014 and the Reference Elevation Model of Antarctica mosaic (REMA mosaic) at an 8 m spatial resolution derived from optical data acquired between 2011 and 2017 have specific individual limitations in this area. The TDM DEM has the advantage of good data consistency and few data voids (approx. 0.85 %), but there exist residual systematic elevation errors such as phase-unwrapping errors in the non-edited DEM version. The REMA mosaic has high absolute vertical accuracy, but on the AP it suffers from extended areas with data voids (approx. 8 %). To generate a consistent, gapless and high-resolution topography product of the AP, we fill the data voids in the TDM DEM with newly processed TDM raw DEM data acquired in austral winters of 2013 and 2014 and detect and correct the residual systematic elevation errors (i.e., elevation biases) in the TDM DEM with the support of the accurately calibrated REMA mosaic. Instead of a pixelwise replacement with REMA mosaic elevations, these provide reference values to correct the TDM elevation biases over entire regions detected through a path propagation algorithm. The procedure is applied iteratively to gradually correct the errors in the TDM DEM from a large to small scale. The proposed method maintains the characteristics of an InSAR-generated DEM and is minimally influenced by temporal or penetration differences between the TDM DEM and REMA mosaic. The performance of the correction is evaluated with laser altimetry data from Operation IceBridge and ICESat-2 missions. The overall root mean square error (RMSE) of the corrected TDM DEM has been reduced from more than 30 m to about 10 m which together with the improved absolute elevation accuracy indicates comparable values to the REMA mosaic. The generated high-resolution DEM depicts the up-to-date topography of the AP in detail and can be widely used for interferometric applications as well as for glaciological studies on individual glaciers or at regional scales.

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.

Tectonic geomorphology and active faults in the Bolivian Amazon

The Bolivian Amazon foreland basin is controlled by i) the sediment input from the Andes and the topographic load; ii) the basinal uplift in the forebulge; and iii) the uplifting Fitzcarrald Arch in the north. These boundary conditions should lead to bending and active faulting of the crust. However, although the presence of active faults has been inferred from geomorphological studies, neither seismicity studies nor geodetic methods and geomorphological constraints have been successful in identifying active faults so far. Reasons for this are mainly the dense vegetation and the extremely dynamic fluvial landscape of the Bolivian Amazon. In this study we use the TanDEM-X and SRTM digital elevation models to show that there are at least five faults, several tens of kilometres long, that are active. We present offset measurements from fault scarps and combine these data with river system analyses. Our results suggest that tilting of tectonic blocks occurred in the Late Pleistocene or Holocene. The reversal of fluvial longitudinal profiles and incipient creek erosion inside paleochannels, markers of young tectonic activity, allow us to reconstruct the tilting directions. We find two orthogonal main faulting directions: WNW-ESE trending faults parallel to the Andean uplift and NNE-SSW trending faults tangent to the bulge of the uplifting Fitzcarrald Arch. Thus, our study area records the interplay between these two main drivers of crustal deformation. We show that some fault scarps were caused by large single earthquakes with magnitudes of around M W 7, while others have recorded multiple strong seismic events. These magnitudes exceed by far the magnitudes known from instrumental seismicity in the last 40 years. Our results shed light on the seismic hazard of the region and highlight that catastrophic river course changes due to sudden uplift in this low-relief region are an underestimated hazard. • We have identified four previously unknown active faults in the Bolivian Amazon. • These faults suggest that here neotectonics are driven by the interplay between the forebulge and the Fitzcarrald Arch. • Faults and tilting blocks control the evolution of the fluvial landscape of the Bolivian Amazon. • Seismic risk in the Bolivian Amazon is higher than previously thought.

Remote Sensing Investigation of the Offset Effect between Reservoir Impoundment and Glacier Meltwater Supply in Tibetan Highland Catchment

This article presents multi-source remote sensing measurements to quantify the water impoundment and regulation of the Zhikong Reservoir (ZKR) and Pangduo Reservoir (PDR), together with the estimation of the glacier mass balance to explore whether the increased glacier meltwater supply can buffer the influences of the reservoir impoundment to some degree in the Tibetan highland catchment. The ZKR and PDR are two reservoirs constructed on the upper Lhasa River that originate from the Nyainqentanglha glaciers in the remote headwater in the Tibetan Plateau (TP) and lacks historical in situ hydrological observations in the long term. Therefore, the Joint Research Center (JRC) Global Surface Water dataset (GSW), and the Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) data were used for estimating the total amount of water storage of the two reservoirs, and the SRTM and TanDEM-X DEMs were used for estimating the glacier mass balance. The result shows that the total amount of water impounded by reservoirs is 0.76 Gt, roughly 54% of their design capacities. The mass balance of the glaciers is estimated by comparing the elevation changes between the SRTM and TanDEM-X DEMs. The glaciers in this region melt at an average rate of 0.09 ± 0.02 Gt·year−1 from 2000 to circa 2013, and the impounded water of these reservoirs is comparable to the amount of glacier-fed meltwater in eight years.

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.

Detection of Low Elevation Outliers in TanDEM-X DEMs With Histogram and Adaptive TIN

Low elevation outliers (LEOs) in TanDEM-X (TDX) digital elevation models (DEMs) have a great effect on the accuracy of digital terrain models (DTMs) generated by filtering DEM data. A novel method, combining local elevation histograms and adaptive triangulated irregular network (HATIN), was developed to detect LEOs in TDX DEMs. The performance of HATIN was evaluated by comparing LEOs identified by median, morphological, HATIN, and iterative inverse distance weighted (IIDW) interpolation methods with LEOs defined by differences between TDX DEMs and light detection and ranging (LiDAR) DTMs. The HATIN method detected almost all large clustered LEOs in the downtown area of Miami, USA with many high-rise buildings, and reduced the root-mean-square error (RMSE) of the resulting DTM from 10.37 to 1.11 m. In the same area, the IIDW method missed some large LEOs, producing a DTM with an RMSE slightly larger than that of the HATIN method. In contrast, the median and morphological methods made more omission and commission errors, failing particularly in the detection of many large LEOs, and producing larger RMSEs (2.43–2.91 m). In low-rise residential areas with isolated LEOs, the HATIN method also identified all large LEOs, producing the fewest commission errors of all methods. However, differences among DTM RMSEs for all four methods were minor. Therefore, it is recommended that, before TDX DEMs are filtered to generate DTMs for urban areas, isolated LEOs should be removed using the computationally simple median or morphological methods, whereas clustered LEOs should be removed by the complex and computationally intensive HATIN or IIDW methods.

Investigating the bias of TanDEM-X digital elevation models of glaciers on the Tibetan Plateau: impacting factors and potential effects on geodetic mass-balance measurements

Abstract The TanDEM-X DEM is a valuable data source for estimating glacier mass balance. However, the accuracy of TanDEM-X elevation over glaciers can be affected by microwave penetration and phase decorrelation. To investigate the bias of TanDEM-X DEMs of glaciers on the Tibetan Plateau, these DEMs were subtracted from SPOT-6 DEMs obtained around the same time at two study sites. The average bias over the studied glacier areas in West Kunlun (175.0 km2) was 2.106 ± 0.012 m in April 2014, and it was 1.523 ± 0.011 m in Geladandong (228.8 km2) in October 2013. By combining backscatter coefficients and interferometric coherence maps, we found surface decorrelation and baseline decorrelation can cause obvious bias in addition to microwave penetration. If the optical/laser data and winter TanDEM-X data were used as new and historic elevation sources for mass-balance measurements over an arbitrary observation period of 10 years, the glacier mass loss rates in West Kunlun and Geladandong would be potentially underestimated by 0.218 ± 0.016 and 0.158 ± 0.011 m w.e. a−1, respectively. The impact is therefore significant, and users should carefully treat the bias of TanDEM-X DEMs when retrieving a geodetic glacier mass balance.

GIS-based rockfall hazard zones modeling along the coastal Gulf of Aqaba Region, Egypt

Rockfall is a natural hazard in mountainous areas not to be underestimated. Mass activities differing in rock volume may cause considerable economic damage. Accomplishing qualitative appraisal of high-potential zones for rockfall is a first step towards implementing mitigation strategies. Nowadays, Geographical Information Systems (GIS) are the state-of-the-art tool for a fast and economic approach of identifying potential hazard zones rather than using conventional mapping with in-situ field data. Primarily, current research focuses on designing and implementing user-friendly tools delineating potential rockfall hazard zonation (RHZ). The constructed model examines triggering factors like slope, aspect, elevation, lithology, structural lineament, rainfall intensity, and seismic activity focal depth of a mountainous coastal region (Gulf of Aqaba, Egypt). The extracted geomorphological parameters were based on a high-resolution TanDEM-X Digital Elevation Model. The enhanced Landsat ETM + 7 was used to generate the lithological and structural lineament parameters, while the rainfall data were collected from NASA project tool. The zonation model was implemented by means of ESRI’s ArcGIS Pro ModelBuilder. Google Earth Pro orthophotos compared with the generated rockfall hazard zonation map indicate the potential RHZ with high reliability. The achieved results show that 15 % of the study area qualifies as a high rockfall hazard zone. As the RHZs generated by the model depend on the input data and the selected rating scores and weights, obtaining ground truth is essential to get a trustworthy result. Finally, this study recommends employing the built RHZ model on similar terrains worldwide to support decision-makers involving any sustainable development projects.

A new look at the Butic Canal, Egypt

Abstract. The Butic Canal – a Roman period transversal route across the northern Nile Delta – was the longest artificial watercourse in the Nile Delta, yet it remains very poorly understood. To date, the canal has not yet been verified by archeological excavations. The route of the eastern section of the canal has been indirectly identified based on a linear elevated feature most likely representing earth from the excavation of the canal. This study combines the analysis of historical sources and remote sensing data, such as satellite imagery and the TanDEM-X digital elevation model, in order to discuss its date of construction, route, and functions. Based on the data of the digital elevation model, new constructional features are visible in the eastern delta providing the first detailed route of a Roman-era artificial watercourse in Egypt. It is suggested that the canal's construction is placed in the context of imperial investments in the infrastructure of the eastern part of the Roman empire.

Natural Pans as an Important Surface Water Resource in the Cuvelai Basin—Metrics for Storage Volume Calculations and Identification of Potential Augmentation Sites

Numerous ephemeral rivers and thousands of natural pans characterize the transboundary Iishana-System of the Cuvelai Basin between Namibia and Angola. After the rainy season, surface water stored in pans is often the only affordable water source for many people in rural areas. High inter- and intra-annual rainfall variations in this semiarid environment provoke years of extreme flood events and long periods of droughts. Thus, the issue of water availability is playing an increasingly important role in one of the most densely populated and fastest growing regions in southwestern Africa. Currently, there is no transnational approach to quantifying the potential storage and supply functions of the Iishana-System. To bridge these knowledge gaps and to increase the resilience of the local people’s livelihood, suitable pans for expansion as intermediate storage were identified and their metrics determined. Therefore, a modified Blue Spot Analysis was performed, based on the high-resolution TanDEM-X digital elevation model. Further, surface area–volume ratio calculations were accomplished for finding suitable augmentation sites in a first step. The potential water storage volume of more than 190,000 pans was calculated at 1.9 km3. Over 2200 pans were identified for potential expansion to facilitate increased water supply and flood protection in the future.

Geomorphological and Climatic Drivers of Thermokarst Lake Area Increase Trend (1999–2018) in the Kolyma Lowland Yedoma Region, North-Eastern Siberia

Thermokarst lakes are widespread in Arctic lowlands. Under a warming climate, landscapes with highly ice-rich Yedoma Ice Complex (IC) deposits are particularly vulnerable, and thermokarst lake area dynamics serve as an indicator for their response to climate change. We conducted lake change trend analysis for a 44,500 km2 region of the Kolyma Lowland using Landsat imagery in conjunction with TanDEM-X digital elevation model and Quaternary Geology map data. We delineated yedoma–alas relief types with different yedoma fractions, serving as a base for geospatial analysis of lake area dynamics. We quantified lake changes over the 1999–2018 period using machine-learning-based classification of robust trends of multi-spectral indices of Landsat data and object-based long-term lake detection. We analyzed the lake area dynamics separately for 1999–2013 and 1999–2018 periods, including the most recent five years that were characterized by very high precipitation. Comparison of drained lake basin area with thermokarst lake extents reveal the overall limnicity decrease by 80% during the Holocene. Current climate warming and wetting in the region led to a lake area increase by 0.89% for the 1999–2013 period and an increase by 4.15% for the 1999–2018 period. We analyzed geomorphological factors impacting modern lake area changes for both periods such as lake size, elevation, and yedoma–alas relief type. We detected a lake area expansion trend in high yedoma fraction areas indicating ongoing Yedoma IC degradation by lake thermokarst. Our concept of differentiating yedoma–alas relief types helps to characterize landscape-scale lake area changes and could potentially be applied for refined assessments of greenhouse gas emissions in Yedoma regions. Comprehensive geomorphological inventories of Yedoma regions using geospatial data provide a better understanding of the extent of thermokarst processes during the Holocene and the pre-conditioning of modern thermokarst lake area dynamics.

Reconstruction and Evaluation of DEMs From Bistatic Tandem-X SAR in Mountainous and Coastal Areas of China

TerraSAR-X add-on for Digital Elevation Measurements (TanDEM-X) mission is designed to generate 3-D images of the Earth as the first bistatic synthetic aperture radar (SAR). However, few quantitative studies of TanDEM-X digital elevation model (DEM) quality validation have been conducted specifically in China. This article presents an iterative method to generate high-resolution TanDEM-X DEMs and assesses the vertical accuracy with high-accuracy GPS observations, 1 arc second global DEMs available [Advanced Land Observing Satellite World 3-D-30 m v2.2, Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM) v3, Shuttle Radar Topography Mission (SRTM) v3.1, NASADEM, and X-band SRTM], and TanDEM-X 90 m DEM. The results demonstrate remarkable elevation quality and consistency in coastal areas with root mean square error of 1.7 m and 90% linear error (LE90) of 0.4 m, whereas 3–4 times weaker accuracies in steep mountainous areas. A positive bias of 1–2 m for an overall LE90 measure exists in the dense vegetation and steep-slope mountainous areas. TanDEM-X DEM-based SAR interferometry deformation uncertainty simulation indicates a low or even negligible topographic error contribution of 2–4 mm in mountainous areas and less than 1 mm in coastal areas. It indicates that the TanDEM-X DEM performs better than other global DEMs overall and shows a better elevation consistence with SRTM C-band DEM in the coastal area. As an excellent source of up-to-date information, the TanDEM-X DEM are expected be an advantage for understanding dynamic land use changes and improving identification and delineation of coastal lands, mountainous landslides, and earthquakes disasters.