Marine-terminating glaciers dominate the evolution of the Greenland Ice Sheet (GrIS) and its contribution to sea-level rise. Widespread glacier acceleration has been linked to the warming of ocean waters around the periphery of Greenland but a lack of information on the bathymetry of the continental shelf and glacial fjords has limited our ability to understand how subsurface, warm, salty ocean waters of Atlantic origin (AW) reach the glaciers and melt them from below. Here, we employ high-resolution, airborne gravity data (AIRGrav) in combination with multibeam echo sounding (MBES) data, to infer the bathymetry of the coastal areas of Northwest Greenland for NASA’s Ocean Melting Greenland (OMG) mission. High-resolution, AIRGrav data acquired on a 2 km spacing, 150 m ground clearance, with 1.5 mGal crossover error, is inverted in three dimensions to map the bathymetry. To constrain the inversion away from MBES data, we compare two methods: one based on the Direct Current (DC) shift of the gravity field (absolute minus observed gravity) and another based on the density of the bedrock. We evaluate and compare the two methods in areas with complete MBES coverage. We find the lowest standard error in bed elevation (±60 m) using the DC shift method. When applied to the entire coast of Northwest Greenland, the three-dimensional inversion reveals a complex network of connected sea bed channels, not known previously, that provide natural and varied pathways for AW to reach the glaciers across the continental shelf. The study demonstrates that the gravity approach offers an efficient and practical alternative to extensive ship mapping in ice-filled waters to obtain information critical to understanding and modeling ice-ocean interaction along ice sheet margins.
Pre-Wisconsin stratigraphy and paleoclimates off Atlantic Canada, and its bearing on glaciation in Québec