Revealing the former bed of Thwaites Glacier using sea-floor bathymetry

Abstract. The geometry of the sea floor beyond Thwaites Glacier (TG) is a major control on the routing of warm ocean waters towards the ice stream’s grounding zone, which has led to increased mass loss through sub-ice-shelf melting and resulting accelerated ice flow. Nearshore topographic highs act as pinning points for the Thwaites Ice Shelf and potentially provide barriers to warm water incursions. To date, few vessels have been able to access this area due to persistent sea-ice and iceberg cover. This critical data gap was addressed in 2019 during the first cruise of the International Thwaites Glacier Collaboration (ITGC) project, with more than 2000 km2 of new multibeam echo-sounder data (MBES) were acquired offshore TG. Here, these data along with legacy MBES datasets are compiled to produce a set of standalone bathymetric grids for the inner Amundsen Sea shelf beyond both Pine Island and Thwaites glaciers. At TG, the bathymetry is dominated by a > 1200 m deep, structurally-controlled trough and discontinuous ridge, on which the Eastern Ice Shelf is pinned. The geometry and composition of the ridge varies spatially with some parts having distinctive flat-topped morphologies produced as their tops were planed-off by erosion at the base of the seaward-moving Thwaites Ice Shelf, suggesting a positive feedback mechanism for ice-shelf ungrounding. Knowing that this offshore area is a former bed for TG, we applied a novel spectral approach to investigate bed roughness and find that derived power spectra can be approximated using an inverse-square law, a result that is consistent with spectra for bed profiles from the modern TG. Using existing ice-flow theory, we also make a first assessment of the form drag (basal drag contribution) for ice flow over this topography. Ice flowing over the sea-floor troughs and ridges would have been affected by similarly high basal drag to that acting in the grounding zone today. We show that the sea-floor bathymetry is an analogue for extant bed areas of TG and that more can be gleaned from these 3D bathymetric datasets regarding the likely spatial variability of bed roughness and bed composition types underneath TG. Comparisons with existing regional bathymetric compilations for the area show that high-frequency (finer than 5 km) bathymetric variability beyond Antarctic ice shelves can only be resolved by observations such as MBES and that without these data calculations of the capacity of bathymetric troughs, and thus oceanic heat flux, may be significantly underestimated. This work meets the requirements of recent numerical ice-sheet and ocean modelling studies that have recognised the need for accurate and high-resolution bathymetry to determine warm water routing to the grounding zone and, ultimately, for predicting glacier retreat behaviour.