Towards Delineating Channel Meltwater Extent Using Sentinel-1 on Greenland: An Assessment of Challenges and Limitations

Global warming is accelerating the melting of glaciers and ice sheets, leading to the formation of supraglacial meltwater, which accumulates in supraglacial lakes and significantly impacts glacier stability. Supraglacial channels can redistribute meltwater and improve glacier stability and are curre...

Full description

Bibliographic Details
Main Author: van Boldrik, Niels (author)
Other Authors: Wouters, B. (mentor), Lhermitte, S.L.M. (graduation committee), Delft University of Technology (degree granting institution)
Format: Master Thesis
Language:English
Published: 2024
Subjects:
Online Access:http://resolver.tudelft.nl/uuid:cee56088-e8f7-4693-8009-2f59af530962
Description
Summary:Global warming is accelerating the melting of glaciers and ice sheets, leading to the formation of supraglacial meltwater, which accumulates in supraglacial lakes and significantly impacts glacier stability. Supraglacial channels can redistribute meltwater and improve glacier stability and are currently monitored using optical data from Sentinel-2. Using Sentinel-1 synthetic aperture radar data omits issues related to cloud cover and darkness, which is the primary focus of this study. The study focuses on three glacier site subsections: Nioghalvfjerdsbrae Glacier, Humboldt Glacier, and Russell Glacier. Sentinel-1 data is compared to the well-established optical data from Sentinel-2, which serves as the reference for channel delineation. The methodology encompasses both thresholding techniques and a GLCM-assisted random forest regression for co- and cross-polarized Sentinel-1 data, while Sentinel-2 data utilizes a thresholding technique based on Glen et al. (2024) and a path-opening algorithm developed by Yang et al. (2015) to delineate narrow channels from enhanced NDWI images. Google Earth Engine is employed for data access and preprocessing, whereas Python is leveraged to conduct the subsequent analysis. Preliminary studies have shown that individual Sentinel-1 scenes may only partially capture supraglacial channels. Therefore, individual scenes are aggregated by taking the backscatter minimum (assumed to be indicative of water due to specular reflection) and standard deviation (representative of change) over a melt season. Composite images offer a more comprehensive view of the supraglacial drainage network during the melt season and form the basis for GLCM texture feature calculations. A 3-by-3 kernel is used around pixels to generate an additional 64 input features by calculating the contrast, correlation, heterogeneity, and dissimilarity over 4 discrete angles for the HH and HV minima and variance. The backscatter features are fed through a random forest regressor with the per-pixel water fraction (i.e. ...