Automatic Detection of Subglacial Water Bodies in the AGAP Region, East Antarctica, Based on Short-Time Fourier Transform

Abstract

Subglacial water bodies are critical components in analyzing the instability of the Antarctic ice sheet. Their detection and identification normally rely on geophysical and remote sensing methods such as airborne radar echo sounding (RES), ground seismic, and satellite/airborne altimetry and gravity surveys. In particular, RES surveys are able to detect basal terrain with a relatively high accuracy that can assist with the mapping of subglacial hydrology systems. Traditional RES processing methods for the identification of subglacial water bodies mostly rely on their brightness in radargrams and hydraulic flatness. In this study, we propose an automatic method with the main objective to differentiate the basal materials by quantitatively evaluating the shape of the A-scope waveform near the basal interface in RES radargrams, which has been seldom investigated. We develop an automatic algorithm mainly based on the traditional short-time Fourier transform (STFT) to quantify the shape of the A-scope waveform in radargrams. Specifically, with an appropriate window width applied on the main peak of each A-scope waveform in the RES radargram, STFT shows distinct and contrasting frequency responses at the ice-water interface and ice-rock interface, which is largely dependent upon their different reflection characteristics at the basal interface. We apply this method on 882 RES radargrams collected in the Antarctic’s Gamburtsev Province (AGAP) in East Antarctica. There are 8822 identified A-scopes with the calculated detection value larger than the set threshold, out of the overall 1,515,065 valid A-scopes in these 882 RES radargrams. Although these identified A-scopes only takes 0.58% of the overall A-scope population, they show exceptionally continuous distribution to represent the subglacial water bodies. Through a comprehensive comparison with existing inventories of subglacial lakes, we successfully verify the validity and advantages of our method in identifying subglacial water bodies using the detection probability for other basal materials of theoretically the highest along-track resolution. The frequency signature obtained by the proposed joint time–frequency analysis provides a new corridor of investigation that can be further expanded to multivariable deep learning approaches for subglacial and englacial material characterization, as well as subglacial hydrology mapping.