skip to main content


Title: Evaluation of high‐resolution DEMs from satellite imagery for geomorphic applications: A case study using the SETSM algorithm
Abstract

High‐resolution digital elevation models (DEMs) have revolutionized research in geomorphology by allowing for detailed quantitative analysis of Earth's surface. Satellite stereo images offer the promise of expanding the availability of high‐resolution DEMs over broad areas, but rigorous evaluation of the scientific application of these datasets remains limited. In this study, we consider DEMs built using stereo pairs of high‐resolution (0.5 m) satellite imagery and the open‐source DEM extraction algorithm, Surface Extraction from TIN Space‐search Minimization (SETSM). We selected locations across a range of landscapes to evaluate the application of these DEMs to geomorphic problems, with particular attention to hillslope analyses where high spatial resolution has been shown to be important for revealing topographic signatures of tectonic and environmental processes. We compared the quality of SETSM 2 m DEMs to LiDAR‐derived DEMs and the widely available SRTM‐30 m and ALOS‐30 m DEMs by comparing the elevation data and derivative products (e.g., slope, aspect, and curvature). We found that SETSM DEMs performed noticeably better than SRTM and ALOS DEMs, but with systematic biases relative to LiDAR DEMs in regions with vegetation. Moreover, noise in the initial SETSM elevation data is amplified with every subsequent derivative, significantly decreasing quality. Finally, we evaluated the potential use of SETSM products for change detection. Applying DEM differencing to a major landslide, we found volume and sediment thickness from SETSM DEMs were similar to volumes and thicknesses from other studies. This example illustrates the capabilities of SETSM and other satellite‐based stereo‐photogrammetry for contributing to rapid response after natural disasters. Overall, we conclude that DEMs derived from satellite image stereo‐photogrammetry can markedly improve on lower resolution global elevation data for terrain analysis and can open possibilities for change detection, but that care needs to be taken in their application especially in regions with significant vegetation.

 
more » « less
NSF-PAR ID:
10363955
Author(s) / Creator(s):
 ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Earth Surface Processes and Landforms
Volume:
47
Issue:
3
ISSN:
0197-9337
Page Range / eLocation ID:
p. 706-722
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Bedrock river width is an essential geometric parameter relevant to understanding flood hazards and gauging station rating curves, and is critical to stream power incision models and many other landscape evolution models. Obtaining bedrock river width measurements, however, typically requires extensive field campaigns that take place in rugged and steep topography where river access is often physically challenging. Although prior work has turned to measuring channel width from satellite imagery, these data present a snapshot in time, are typically limited to rivers ≥ 10–30 m wide due to the image resolution, and are physically restricted to areas devoid of vegetation. For these reasons, we are generally data limited, and the factors impacting bedrock channel width remain poorly understood. Due to these limitations, researchers often turn to assumptions of width‐scaling relationships with drainage area or discharge to estimate bedrock channel width. Here we present a new method of obtaining bedrock channel width at a desired river discharge through the incorporation of a high‐resolution bare‐earth digital elevation model (DEM) using MATLAB Topotoolbox and the HEC‐RAS river analysis system. We validate this method by comparing modeled results to US Geological Survey (USGS) field measurements at existing gauging stations, as well as field channel measurements. We show that this method can capture general characteristics of discharge rating curves and predict field‐measured channel widths within uncertainty. As high‐resolution DEMs become more available across the United States through the USGS three‐dimensional elevation program (3DEP), the future utility of this method is notable. Through developing and validating a streamlined, open‐source, and freely available workflow of channel width extraction, we hope this method can be applied to future research to improve the quantity of channel width measurements and improve our understanding of bedrock channels.

     
    more » « less
  2. Abstract

    The availability of high‐resolution digital elevation models (DEMs) derived from airborne light detection and ranging (LiDAR) surveys has spurred the development of several methods to identify and map fluvial terraces. The post‐glacial landscape of the Sheepscot River watershed, Maine, where land‐use change has produced fill terraces upstream of historic dam sites, was selected to implement a comparison between terrace mapping methodologies. At four study sites within the watershed, terraces were manually mapped on LiDAR‐DEM‐derived hillshade images to facilitate the comparison among fully and semi‐automated DEM‐based procedures, including: (1) spatial relationships between interpreted terraces and surrounding natural topography, (2) feature classification algorithms, and (3) the TerEx terrace mapping toolbox. Each method was evaluated based on its accuracy and ease of implementation. The four study sites have varying longitudinal slope (0.0008–0.006 m/m), channel width (< 5–30 m), surrounding landscape relief (20–80 m), type and density of surrounding land use, and mapped surficial geologic units. All methods generally overestimate terrace areas (average predicted area 210% of manually defined area) with the most accurate results achieved within confined river valleys surrounded by the steep hillslopes. Accuracy generally decreases for study sites surrounded by low‐relief landscapes (predicted areas ranged 4–953% of manual delineations). We conclude with the advantages and drawbacks of each method tested and make recommendations for the scenarios where the use of each method is most appropriate. Copyright © 2016 John Wiley & Sons, Ltd.

     
    more » « less
  3. Abstract. Here we describe the curriculum and outcomes from a data-intensivegeomorphic analysis course, “Geoscience Field Issues Using High-ResolutionTopography to Understand Earth Surface Processes”, which pivoted to virtualin 2020 due to the COVID-19 pandemic. The curriculum covers technologies formanual and remotely sensed topographic data methods, including (1) GlobalPositioning Systems and Global Navigation Satellite System (GPS/GNSS)surveys, (2) Structure from Motion (SfM) photogrammetry, and (3) ground-based(terrestrial laser scanning, TLS) and airborne lidar. Course content focuseson Earth-surface process applications but could be adapted for othergeoscience disciplines. Many other field courses were canceled in summer2020, so this course served a broad range of undergraduate and graduatestudents in need of a field course as part of degree or researchrequirements. Resulting curricular materials are available freely within theNational Association of Geoscience Teachers' (NAGT's) “Teaching with Online Field Experiences” collection. Theauthors pre-collected GNSS data, uncrewed-aerial-system-derived (UAS-derived) photographs, and ground-based lidar, which students then used in courseassignments. The course was run over a 2-week period and had synchronousand asynchronous components. Students created SfM models that incorporatedpost-processed GNSS ground control points and created derivative SfM and TLSproducts, including classified point clouds and digital elevation models(DEMs). Students were successfully able to (1) evaluate the appropriatenessof a given survey/data approach given site conditions, (2) assess pros andcons of different data collection and post-processing methods in light offield and time constraints and limitations of each, (3) conduct error andgeomorphic change analysis, and (4) propose or implement a protocol to answera geomorphic question. Overall, our analysis indicates the course had asuccessful implementation that met student needs as well as course-specificand NAGT learning outcomes, with 91 % of students receiving an A, B, or Cgrade. Unexpected outcomes of the course included student self-reflectionand redirection and classmate support through a daily reflection anddiscussion post. Challenges included long hours in front of a computer,computing limitations, and burnout because of the condensed nature of thecourse. Recommended implementation improvements include spreading the courseout over a longer period of time or adopting only part of the course andproviding appropriate computers and technical assistance. This paperand published curricular materials should serve as an implementation andassessment guide for the geoscience community to use in virtual or in-personhigh-resolution topographic data courses that can be adapted for individuallabs or for an entire field or data course. 
    more » « less
  4. Abstract. Ocean-induced basal melting is directly and indirectly responsible for much of the Amundsen Sea Embayment ice loss in recent decades, but the total magnitude and spatiotemporal evolution of this melt is poorly constrained. To address this problem, we generated a record of high-resolution Digital Elevation Models (DEMs) for Pine Island Glacier (PIG) using commercial sub-meter satellite stereo imagery and integrated additional 2002&ndash;2015 DEM/altimetry data. We implemented a Lagrangian elevation change (Dh/Dt) framework to estimate ice shelf basal melt rates at 32&ndash;256-m resolution. We describe this methodology and consider basal melt rates and elevation change over the PIG shelf and lower catchment from 2008&ndash;2015. We document the evolution of Eulerian elevation change (dh/dt) and upstream propagation of thinning signals following the end of rapid grounding line retreat around 2010. Mean full-shelf basal melt rates for the 2008&ndash;2015 period were ~82&ndash;93&thinsp;Gt/yr, with ~&thinsp;200&ndash;250&thinsp;m/yr basal melt rates within large channels near the grounding line, ~&thinsp;10&ndash;30&thinsp;m/yr over the main shelf, and ~&thinsp;0&ndash;10&thinsp;m/yr over the North and South shelves, with the notable exception of a small area with rates of ~&thinsp;50&ndash;100&thinsp;m/yr near the grounding line of a fast-flowing tributary on the South shelf. The observed basal melt rates show excellent agreement with, and provide context for, in situ basal melt rate observations. We also document the relative melt rates for km-scale basal channels and keels at different locations on the shelf and consider implications for ocean circulation and heat content. These methods and results offer new indirect observations of ice-ocean interaction and constraints on the processes driving sub-shelf melting beneath vulnerable ice shelves in West Antarctica.

     
    more » « less
  5. ArcticDEM provides the public with an unprecedented opportunity to access very high-spatial resolution digital elevation models (DEMs) covering the pan-Arctic surfaces. As it is generated from stereo-pairs of optical satellite imagery, ArcticDEM represents a mixture of a digital surface model (DSM) over a non-ground areas and digital terrain model (DTM) at bare grounds. Reconstructing DTM from ArcticDEM is thus needed in studies requiring bare ground elevation, such as modeling hydrological processes, tracking surface change dynamics, and estimating vegetation canopy height and associated forest attributes. Here we proposed an automated approach for estimating DTM from ArcticDEM in two steps: (1) identifying ground pixels from WorldView-2 imagery using a Gaussian mixture model (GMM) with local refinement by morphological operation, and (2) generating a continuous DTM surface using ArcticDEMs at ground locations and spatial interpolation methods (ordinary kriging (OK) and natural neighbor (NN)). We evaluated our method at three forested study sites characterized by different canopy cover and topographic conditions in Livengood, Alaska, where airborne lidar data is available for validation. Our results demonstrate that (1) the proposed ground identification method can effectively identify ground pixels with much lower root mean square errors (RMSEs) (<0.35 m) to the reference data than the comparative state-of-the-art approaches; (2) NN performs more robustly in DTM interpolation than OK; (3) the DTMs generated from NN interpolation with GMM-based ground masks decrease the RMSEs of ArcticDEM to 0.648 m, 1.677 m, and 0.521 m for Site-1, Site-2, and Site-3, respectively. This study provides a viable means of deriving high-resolution DTM from ArcticDEM that will be of great value to studies focusing on the Arctic ecosystems, forest change dynamics, and earth surface processes. 
    more » « less