Long-term Satellite-derived Bathymetry of Arctic Supraglacial Lake from ICESat-2 and Sentinel-2

Abstract. The polar ice sheets serve as natural thermostats, regulating Earth’s temperatures. The Greenland Ice Sheet (GrIS), the second-largest ice sheet, is a critical indicator of climate change and global warming. Estimating the volume of supraglacial lakes on the GrIS, which is directly linked to the extent of melting in the Arctic ice sheet, requires information on both lake area and water depth. Conventional bathymetric methods (i.e., airborne bathymetric LiDAR, shipborne echo-sounder) are commonly used for accurate water depth measurement. However, polar supraglacial lakes face challenging conditions, leading to uncertainties in their spatial and temporal distribution. To overcome the limitations, this study combines Sentinel-2 and ICESat-2 (Ice, Cloud, and Land Elevation Satellite-2) to estimate bathymetry and detect changes in lake volume on the GrIS from 2019 to 2023. Firstly, Sentinel-2 images were pre-processed, and ICESat-2 single-photon LiDAR points were extracted using the DBSCAN (Density-Based Spatial Clustering of Applications with Noise) method, followed by the bathymetric corrections as the training data. Subsequently, three bathymetry models (i.e., log-linear, log-ratio, and BP (Back Propagation) neural network) were constructed using Sentinel-2 images and ICESat-2 data. Lastly, the highresolution ArcticDEM (Arctic Digital Elevation Model) was used as the validation data to assess the satellite-derived bathymetry accuracy. In this study, the log-ratio model yielded the best results with the R2, RMSE, and MAE of 0.92, 0.79 m (lower than 10% of the maximum depth), and 0.62 m. The results demonstrate the feasibility of the integrated active and passive remote sensing approach for bathymetry in Arctic supraglacial lakes.