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

. 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 20 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 km 2 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 25 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 30 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

composition) of topographic highs as former grounding zones and ice-shelf pinning points, and to better resolve sea-floor troughs as potential modified CDW (CDW) pathways to the modern Thwaites grounding zone. In the second part of the paper, we use a spectral approach to relate flow-line roughness to the drag contribution for an overriding ice mass using bathymetric profiles from the inner ASE shelf, and bed profiles from upstream areas of Pine Island and Thwaites glaciers. To reflect these 105 two components, (1) we describe the observational (geophysical) datasets used and interpret the new bathymetric dataset, which is also provided as a publically-available standalone grid and (2) we describe the spectral approach and use it to quantitatively examine former and modern TG beds. We demonstrate that the offshore area just seaward of Thwaites Ice Shelf is an appropriate analogue for the modern grounding zone of TG, both in terms of its bed characteristics and in the effect of its rugged bed topography on ice flow by calculating drag contribution for different scales of roughness. We highlight the 110 importance of high-resolution MBES observational data for constraining gravity inversions and regional bathymetry compilations, which are essential boundary conditions for predictive numerical modelling experiments, and for accurately calculating the flux of warm water to the grounding zone.

Multibeam echosounder datasets (MBES) 115
During ITGC cruise NBP19-02, the marine areas in front of TG were unusually clear of sea-ice and icebergs providing a unique opportunity for bathymetric data acquisition at the margin of Thwaites Ice Shelf . MBES data was acquired using a hull-mounted 1° x 1° Kongsberg EM122 echo sounder with 288 across-track beams and an operational frequency in the range 11.25-12.75 kHz. The MBES was configured with "High Density Equidistant" beam spacing meaning that more than one sounding can be produced per beam (up to 432) effectively increasing across-track resolution, and in "Dual Ping" 120 mode which ensures equal across-and along-track sounding spacing. As an example, in 600 m water depth with a 60° beam angle (typical conditions and settings on NBP19-02), this results in a sounding spacing on the sea floor of ~5-7 m. Here, in addition to this new dataset, we have compiled all available MBES data in the area from UK, US, German, Swedish and Korean expeditions (Table 1; Fig. 2) to produce gridded bathymetric data products.
For NBP19-02, data processing was performed on board using MB-System (Caress and Chayes, 1996;Caress et al., 2019) in 125 order to apply optimal sound velocity profiles (SVPs) to each data file and to remove erroneous soundings. Ray tracing and sea-floor depths were calculated using SVPs generated from conductivity-temperature-depth and expendable bathythermograph casts made during NBP19-02 . Most of the other bathymetry datasets were also initially ping-edited during each respective cruise; however, minor additional cleaning was performed in MB-System and QPS Fledermaus (Mayer et al., 2000) after the datasets were collated when clear outliers could be easily identified. Ultimately, the 130 data were gridded in MB-System using a Gaussian weighted mean filter algorithm to produce an isometric 50-m digital elevation model (DEM) for the sea floor on the southern ASE shelf. A degree of interpolation was applied to the final grids in data gaps only, filling areas six cell widths away from cells with real soundings, i.e. for a 50-m grid interpolation will fill cells https://doi.org/10.5194/tc-2020-25 Preprint. Discussion started: 3 February 2020 c Author(s) 2020. CC BY 4.0 License. up to 300 m away from a cell with real soundings. A lower resolution DEM (500-m grid cells) was produced for studies that typically use coarser bathymetric information. Together, the 50-m and 500-m DEMs are presented as a standalone regional 135 mid-resolution bathymetric dataset. In addition, in order to examine the nature of specific sea-floor features (e.g. ice-shelf pinning points), higher-resolution grids (30-m grid cells) were produced where sounding densities allowed (e.g. Fig. 4). Final grids were visualised and analysed in QPS Fledermaus and ArcGIS 10.6. Sea-floor trough and channel metrics (including widths, depths, symmetry, form ratio, u-/v-shape characterisation) were derived using the methods described in Kirkham et al. (2019). 140

Bed profiles for investigating bed roughness and basal drag
The gridded bathymetric datasets that we have produced represent the former bed of an expanded TG. Palaeo-flow lines across this former bed were picked manually by tracing lines parallel to subglacial landforms taken to indicate palaeo-ice flow directions (e.g. Fig. 3b). We also perform this exercise for bed profiles from three additional areas for comparison with the Thwaites offshore area (Figs. S1, S2): (i) modern along-flow bed profiles for Thwaites Glacier and (ii) Pine Island Glacier,145 and (iii) an area of "smooth" bed topography on the middle continental shelf in the Dotson-Getz palaeo-ice stream trough, offshore from the Getz A Ice Shelf, representing a sedimentary palaeo-ice stream bed characterised by mega-scale glacial lineations (MSGL) Spagnolo et al., 2014). Onshore bed profiles were extracted from the AGASEA (Holt et al., 2006) and Operation IceBridge (Cochran & Bell, 2010, updated 2014 datasets for Thwaites and Pine Island glacier beds, respectively; the profile from the Dotson-Getz Trough was extracted from an MBES dataset fully described by Larter et 150 al. (2009) and Graham et al. (2009).

New bathymetric compilation for the inner Amundsen Sea Embayment shelf
Our bathymetric compilation includes more than 2000 km 2 of new MBES data between the EIS and the TGT and west of the TGT that provides near-continuous bathymetric coverage for ~40 km north of the present day ice-shelf margin (Figs. 2, 3).
Data gaps remain in the western part of the area in front of the ice mélange, east of Crosson Ice Shelf, as icebergs and bergy 155 bits released from the ice mélange persistently covered this region during NBP19-02. In addition, perennial fast ice and huge icebergs, such as the 80-km long and 45-km wide iceberg B-22A, that calved periodically from the TGT and then moved slowly north-and northwest, thereby episodically running aground, have generally prevented survey between TGT and Crosson Ice Shelf (see area with B-22A outlines on Fig. 1).
The sea-floor bathymetry offshore Thwaites Ice Shelf is dominated by an elongate depression oriented NNE-SSW (T1) and a 160 series of topographic highs (H1-H3) along its southern margin (Fig. 3a). The depression is characterized by water depths of 1100-1250 m, which is 200 to 400 m deeper than the sea floor on its flanks; it typically has a relatively flat or gently-inclined floor (Figs. 3a,c). Although the depression appears to be continuous for at least 75 km, and connects with areas of deep (>1300 m) sea floor directly north of the EIS and in Pine Island Trough further east (Fig. 2a), its width varies significantly along its https://doi.org/10.5194/tc-2020-25 Preprint. Discussion started: 3 February 2020 c Author(s) 2020. CC BY 4.0 License. length and the flanks are discontinuous in form. As a result, we do not define this feature as a "channel", which implies incision 165 by water flow, but rather as a small trough. The topographic highs that make up the southern flank of the trough (H1-H3) decrease in height (above the surrounding sea floor) from east to west and appear to form a broad (>15 km) but discontinuous ridge, also oriented NNE-SSW. The large, discontinuous ridge is the extension of a bedrock ridge in Pine Island Trough to the east that has the same orientation (Figs. 2, 3a;Nitsche et al., 2013). The most prominent high (H1) occurs immediately north of the EIS with a shallowest recorded water depth of just 82 m and extends at least 40 km in a north-south direction. The 170 northern part of the EIS is pinned, about 45 km downstream of the grounding zone (Rignot et al., 2001;Tinto and Bell, 2011;MacGregor et al., 2012), on the southern part of this high, forming an ice rumple in the EIS (Matsuoka et al., 2015). The shallowest water depths on the two highs H2 and H3 (from NNE to SSW) are 362 m and 611 m, respectively. Indeed, remotesensing data had identified a pinning point for TGT on H2  but the TGT must have fully retreated from that point prior to NBP19-02. North of the trough and ridge (north of 74°30' S) is an area of rugged morphology characterised 175 by shallow sea floor ridges and deep basins (>1400 m); this area merges with similar terrain in Pine Island Trough described by Nitsche et al. (2013) as their "area 2" (Fig. 2). East of the EIS, the eastern flank of H1 is just exposed but the bathymetry is dominated by a deep (1000-1200 m), rugged area of sea floor, bounded on its eastern edge by a bedrock high in Pine Island Bay that continues southward towards the grounding line (Fig. 2a). The deepest part of this area appears to form a poorlydefined bathymetric trough oriented NNE-SSW (T4; Fig. 3a); oceanographic measurements and models confirm that this deep 180 acts as a pathway for CDW towards grounding lines in western Pine Island Bay (Fig. 2a) (Dutrieux et al., 2014;Nakayama et al., 2019).

Glacial landforms
The large-scale morphology of the sea floor offshore Thwaites Ice Shelf is overprinted by linear features oriented sub-parallel to the trough and ridge. Streamlined subglacial landforms occur in areas where bedrock crops out at the sea floor, on either 185 topographic highs or interrupting the smooth trough floor or rugged basin floors (Figs. 3a,b). These features are identified as crag-and-tails by their tapering form and are 750-4000 m long, 200-2000 m wide and 25-250 m high (Figs. 3,4), and their morphological similarity to crag-and-tails from other deglaciated terrains (e.g. Dowdeswell et al., 2016;Maclean et al., 2016;Nitsche et al., 2016). The wider southern ends of the crag-and-tails are rugged and have significant relief, suggesting bedrock composition, whereas the northern ends are smooth and elongate, suggesting sedimentary composition. These landforms are 190 produced subglacially as glacier ice flows across bedrock obstacles producing the characteristic morphology through erosion and deposition (Benn & Evans, 2010;Nitsche et al., 2016). The tapering of these features north-northeastward indicates palaeoice flow of TG in this direction towards Pine Island Trough (Fig. 3b); east of the EIS the orientation of these landforms varies from N-S to NW-SE depicting flow around the H1 high (Fig. 3b). Curved or semi-circular moats -crescentic scours (Lowe and Anderson, 2003;Graham et al., 2009;Graham and Hogan, 2016) -occur around the southern (upstream) sides of some 195 crag-and-tails (Figs. 3,4) and are suggestive of erosion by meltwater (or a till slurry) upstream of bedrock obstacles (cf. Graham & Hogan, 2016). On the H2 and H3 highs, subtle elongate ridges separated by linear grooves also occur. These glacial https://doi.org/10.5194/tc-2020-25 Preprint. Discussion started: 3 February 2020 c Author(s) 2020. CC BY 4.0 License.
lineations have parallel sides, lack a pronounced and wider end, and are 1000-2500 m long, 100-200 m wide and 2-10 m high.
These features, which were produced subglacially, have relatively short elongation ratios (<10:1) and are oriented parallel to crag-and-tails also on the tops of the highs but are slightly oblique to crag-and-tails in the trough (Fig. 3). A distinct but 200 discontinuous scarp is mapped on the top of the H2 high, upstream from its frontal edge (Figs. 3,4b,4d). The scarp has a curved planform and a steep (3-4°) northern and gentle (0.5°) southern slope; glacial lineations occur on the gentle back-slope of this feature which extends for about 15 km in a SW-NE direction. The asymmetric geometry and lineated back-slope of this landform identify it as a grounding-zone wedge (GZW) (Alley et al., 1989;Larter and Vanneste, 1995;Dowdeswell & Fugelli 2012;Jakobsson et al., 2012), i.e. a wedge of sedimentary material that built up at the grounding line when it was stable for a 205 time on the H2 high. This feature is interrupted, however, by a ~40 m deep groove or small channel with bedrock ridges visible on either side suggesting that the sedimentary wedge is not thick enough to fully bury the underlying topography, i.e. it is only tens of metres thick. Discrete 6-10 m deep linear to curvilinear furrow with small berms (4-6 m) were mapped and on H3 and north of the GZW front scarp (Figs. 4a,b); these are interpreted as iceberg ploughmarks.
Sea-floor highs in the area are also variously gouged and streamlined, resulting in a pattern of grooves and bedrock ridges. 210 Grooves on the H1-H3 ridges occur on their southern (upstream) parts and exhibit a range of orientations that probably relate to the structure of the underlying bedrock exploited by glacial erosion. The grooves are typically <20 m deep, <200 m wide and <6000 m long. The surfaces of sea-floor highs north of the H1-H3 ridge have a more streamlined appearance resulting from shallow, semi-parallel shallow grooves that are preferentially aligned with the crag-and-tails (Figs. 2,3). This type of streamlined bedrock highs is typical of inner-shelf morphologies around Antarctica (e.g. Wellner et al., 2006;Livingstone et 215 al., 2013) including in the adjacent eastern part of Pine Island Bay . Taken together, the orientations of the streamlined glacial landforms (crag-and-tails, glacial lineations, bedrock grooves/ridges) define the former ice-flow directions of an expanded TG (depicted by red arrows in Fig. 3b).

Trough and channel metrics
Multiple bathymetric troughs and bedrock channels were mapped and analysed in the area beyond Thwaites Ice Shelf as part 220 of this study (Fig. 3a). Troughs and channels in the adjacent eastern part of Pine Island Bay have been described comprehensively by Nitsche et al. (2013) and Kirkham et al. (2019). The larger troughs in our study area, which have also been identified on gravity-derived regional bathymetry maps (Millan et al., 2017;Jordan et al., 2020.), are considered important as potential pathways for the transport of CDW towards the grounding zone of TG and warrant full description here. We distinguish these comparatively larger troughs, based on their size, connectivity and variable flank form (as described above), 225 from the channels which are notably smaller in scale ( Fig. S3) with continuous, parallel sided-flanks, undulating thalwegs, and incise into rugged sea-floor areas interpreted as bedrock. It is widely accepted that channels of the latter type were eroded by pressurised subglacial water flow (Lowe and Anderson, 2003;Nitsche et al., 2013;Kirkham et al., 2019), whereas the troughs likely relate, at least in part, to underlying geological structures such as mafic dykes and tectonic deformations (cf. ) that have been variously modified by ice.
Cross-sectional analyses of the troughs reveal their large scale, with average widths, depths and cross-sectional areas being 2090 m, 90 m and 144 000 m 2 , respectively (Fig. S3). The troughs are typically 10-30 times as wide as they are deep, although we note that there is a significant size difference between the main NE-SW trough (T1) and the remaining troughs ( Fig. 3a;   Fig. S3). By comparison, the bedrock channels are on average 520 m wide, 50 m deep with a cross-sectional area of 18 000 m 2 . The channels are generally 5-10 times wider than they are deep. The derived b values, which characterise cross-sectional 235 shape (Pattyn and Van Huele, 1998), suggest that the bedrock channels are between v-and u-shaped, whereas the larger troughs have no dominant cross-section shape (Fig. S3). In general, the trough floors are flat or inclined in cross-profiles (Fig. 3c), and are gently undulating in along-trough profiles.
The main T1 trough is oriented NNE-SSW, which is oblique to the northerly palaeo ice-flow directions immediately in front of Thwaites Ice Shelf (Fig. 3b). This indicates that, at the time that the subglacial landforms were produced, ice was thick 240 enough not to be steered by even major elements of the bed topography. The two southernmost troughs (T2 and T3 on Fig. 3a) are oriented perpendicular to TI (i.e. NNW-SSE), and the troughs north and east of the EIS are generally aligned with palaeoice flow directions (Figs. 2, 3a). The T2 and T3 troughs, whose floors have water depths of 800-900 m, separate the H1-H3 bathymetric highs and are of interest as potential pathways for CDW to the TG grounding line.

Bathymetric highs and ridges 245
Owing to the importance of sea-floor highs in front of the Thwaites Ice Shelf as barriers to CDW inflow, and as former ice shelf/sheet pinning points, we examine the morphology of the discontinuous NNE-SSW trending ridge in detail (Fig. 4). The ridge comprises the H1-H3 highs separated by the two troughs described above (T2 and T3; Fig. 3a). The width of the ridge varies significantly, from 6 km in the SW of the study area to at least 40 km over H1, although we acknowledge that data coverage is limited. In places, the bathymetric highs are strikingly flat-topped. These planar features are accentuated in maps 250 of the first derivative of bathymetry, slope, which reveals both low slopes (<2°) (Figs. 4c, d), and low roughness over H2, H3 and the western part of H1. The continuation of the ridge further north into Pine Island Trough has a similar surface expression but is generally narrower (Fig. 3a). These areas with low surface slopes are atypical when compared with other bathymetric highs in the area, which have rugged surface morphologies characterised by bedrock grooves and channels (Fig. 2, 4d) Arndt et al., 2018;Kirkham et al., 2019). Instead, the low slope values are similar to those derived for the base of 255 the troughs in front of Thwaites Ice Shelf (Fig. 4a) and the sediment-filled basins just seaward of the Pine Island Ice Shelf front Kuhn et al., 2017). At least three distinct levels of flat-topped surfaces occur at 400 m, 570 m and 620 m water depth (Fig. 5). We interpret this morphology as being generated during a prolonged period of overriding by a formerly expanded Thwaites Ice Shelf (with the necessary ice thickness). Support for this comes from the orientation of glacial lineations on the tops of the highs (Fig. 3b), which are in line with modern ice-velocity vectors  but 260 oblique to the orientation of crag-and-tails in the troughs. A similar interpretation was made for the lineated surface of a former pinning point of the Pine Island Ice Shelf (PIIS) that has been recently exposed by ice-shelf calving events (Arndt et al., 2018), although that feature was not planed off to form a flat-topped high but rather has a stepped and rugged surface morphology https://doi.org/10.5194/tc-2020-25 Preprint. Discussion started: 3 February 2020 c Author(s) 2020. CC BY 4.0 License. albeit with some gently-sloping parts (see their Fig. S3). The morphology of the highs in our study area is striking and notably rare for sea-floor highs around Antarctica. However, we do identify a similar but less pronounced terrain over some highs 265 along the same structural ridge in eastern Pine Island Bay, i.e. east of the EIS (Fig. 5b), and a solitary planed-off high, with comparable dimensions to those offshore TG, is visible <5 km north of the Getz-A Ice Shelf (Fig. 5c). Outwith these rare examples, the best analogy for this morphology probably comes from a set of "iceberg terraces" on terminal moraines at the mouth of a Svalbard fjord, which display remarkably flat-topped surfaces at several bathymetric levels. These are interpreted to have formed as tabular, flat-based icebergs over-topped and eroded morainal sediments (Noormets et al., 2016). It should 270 be noted, however, that sediments of this morainal bank complex probably consist of unconsolidated material that has not been overridden (or compacted) by grounded ice, meaning that they are likely more erodible than the basement highs in front of Thwaites Ice Shelf.
New geomorphic information is also revealed by the flanks of the highs (Figs. 4a, b). The northern (ice distal) flank of the H2 and H3 highs are characterised by subtle down-slope trending gullies that transition into a smooth but inclined sea floor in the 275 troughs at a distinct break of slope. In addition, there are a few, discrete semi-circular indentations in the scarp between the surface of the highs and their flanks (Fig. 4a). The gullies have simple non-branching geometries, small dimensions (widths 150-700 m, depths 5-50 m, lengths <2 km), and typically define broad u-shapes in cross-section, although some v-shaped forms are present. In addition, their long-profiles have negative slopes (i.e., deepen consistently down the flank of the highs) and their form is consistent with other Antarctic submarine gully systems (e.g. Gales et al., 2013;Post et al, 2019). Thus, we 280 interpret the gullies as the result of the down-slope mass movement of material from the tops and sides of the H2 and H3 highs via gravitational processes into the small sediment fans at the base of the slope (Fig. 4a). The semi-circular indentations may be the headwalls of small slide scars (cf. Noormets et al., 2009;Gales et al., 2013). One 6 km by 2.6 km segment of the H3 high is somewhat detached from other parts of the ridge and appears particularly fragmented on its flanks (Fig. 4a). About 1500 m south of this, on the main H3 high, is a distinct break in slope with the same planform shape as the back of the detached 285 segment that we refer to as a "block" (bl in Figs. 4a,c). Similarly, 2 km west of this block is another somewhat isolated 3 km by 1.8 km block of H3 that is incised by gullies on its northern front and lacks lineations on its surface (Fig. 4a). We consider several interpretations for these features: First, it is possible that they are detached slide blocks of the H3 high that had been displaced down-slope over a short distance (black arrows in Figs. 4a, c), remaining largely intact, but subsequently affected by some gravitational collapse of their flanks. Slide "megablocks" with similar dimensions (or larger), non-crystalline 290 compositions, and degraded flanks are known from, for example, the Hinlopen Slide scar on the northern Barents Sea margin (Vanneste et al., 2006;Hogan et al., 2013). The second possibility is that the blocks are small bedrock highs that have been variously mantled by and surrounded by glacigenic sediment as it was transported down-slope (towards the sediment fans) from the H3 high (see blue arrows on Fig. 4a). If the latter case is true, then the pronounced semi-circular indentation on the eastern block (labelled in Fig. 4a) may be an erosional scour mark formed by persistent bottom water motion, as is evidenced 295 by the small channel on the eastern side of the block (white arrows in Fig. 4c).

Spectral analysis of bed roughness
The high-resolution bathymetric grid that we have produced (Fig. 2) represents the former beds of expanded Thwaites and Pine Island glaciers. As such, it contains information about bed geometry and roughness that can be used to investigate the 300 behaviour of glacier ice over this terrain, specifically how roughness at the bed (e.g. bedrock topography and subglacial landforms) might generate "form drag" (Schoof, 2002) within the ice. The bathymetry contains roughness at various wavelengths, λ, each equivalent to a spatial frequency f=1/λ. We use power spectra of the high-resolution bathymetry (Figs. 6, 7), together with a theoretical expression for form drag (Schoof, 2002, equation 67) to recover the drag caused by undulations at each spatial frequency. The theoretical expression applies to flow lines and was applied to palaeo-and modern flow lines 305 for TG and the Dotson-Getz palaeo-ice stream as described in Section 2.2 (locations in Fig. S1).
To analyse the variance of roughness at each wavelength scale we computed the power spectrum of the bathymetric topography (Figs. 6, 7). As an example, Figure 6a shows a bed elevation profile (location in Fig. S1) sampled at 25-m intervals along a flow line. Figure 6b shows the derived one-sided periodogram ( ), evaluated at a set of equally-spaced spatial frequencies = −1 = / , where is the length of a moving window. The periodogram was obtained using Welch's method (Welch, 310 1967), with a Hamming window of length =6.4 km and 50% overlap between consecutive windows along the bed profile.
The periodogram shows no strong peaks at any particular preferred scales of roughness (Fig. 6b). Instead, the variance of roughness decreases continuously as spatial frequency increases and the reduction of variance with increasing spatial frequency approximately follows an inverse-square power law so that the periodogram can be approximated as = ( ) = −2 .
This is the power spectrum of a fractal 'Brown-noise' or random-walk elevation profile, as would be produced by adding an 315 uncorrelated random offset in the vertical for each unit-length step taken in the horizontal. The inverse-square law provides a good approximation to all power spectra observed (Figs. 7, SF4), so the variance of roughness at any particular frequency scales inversely to the square of that frequency and proportionately to the coefficient, . The red line in Figure 6b corresponds to a value = 0.1 m and the spectra are remarkably consistent across many of the profiles considered (Figs. 7, SF4) with the exception of the smoother MSGL area that is discussed in section 4.2. 320 The power-law approximation to the power spectrum provides reasonable agreement for multiple flow lines from the offshore bathymetry (Figs. 7, SF4). It also agrees closely to power spectra of profiles of bed elevation from airborne radar flown over TG (Fig. 7c), so the MBES data provide an analogue to the subglacial undulations that control the sliding of Thwaites Glacier today. Despite improvements in the methodology of high-resolution radar surveys of the active subglacial bed (King et al., 2016;Bingham et al., 2017), the MBES provides a more detailed view of the shorter spatial scales than airborne or ground-325 based radar.

Relating bed topography to basal drag
In this section, we will use a theory of ice flow to assess how the roughness at these different scales will affect the drag that opposes glacier sliding. As demonstrated above, the power-law approximation is appropriate for bed profiles from both the bathymetric data as well as for bed profiles from airborne radar data. This allows us to consider the contribution that the 330 observed sea-floor bathymetric roughness at short wavelengths would make to form drag when covered by flowing ice. The theory of ice flow over an undulating bed (Schoof, 2002) provides an approximate expression for the form drag  , expressed as a basal shear stress that acts to resist sliding: In this expression, is the speed of ice flow averaged over some horizontal length scale (here, chosen as , the window length, i.e. the longest wavelength considered in the spectral analysis) and β is the drag coefficient. According to Schoof (2002), the is the sum of contributions , each of which is caused by roughness at a different wavelength scale . As per the spectral analysis above, these wavelengths correspond to equally-spaced spatial frequencies = −1 = 340 / , and spatial wavenumbers = 2 . Schoof (2002) provides an expression for the drag contributed by each spatial frequency: In this expression, ĥ is the Fourier component of the bed roughness at spatial frequency . The scaling constants are * = 345 / and * = 1/ , where and are representative viscosity and thickness, respectively, averaged over the length scale (Schoof, 2002). We estimate the squared magnitude of the Fourier component of the roughness using the one-sided periodogram,|ĥ | 2 = 1 2 −2 −1 . This corresponds to the definition of the roughness Fourier series used by Schoof (2002).
We consider two limiting cases for the function : These are derived by Schoof (2002) (equations 65 and 66 for small and large bed roughnesses, respectively). Note that we 355 have used dimensional quantities throughout, so / * , ĥ / and ( / * ) in our notation are equivalent to nondimensionalised quantities , ĥ and ( ) respectively in Schoof's notation. For the same example bed profile in Figure 6a, values of  / * are plotted against spatial frequency in Figure 6c for a representative ice thickness = 1 km. For sufficiently small features, so that ≫ * , the functions 1 ( / * ) 2 ( / * ) both tend to unity. Then the form drag becomes insensitive to ice thickness and to the choice of function used (Schoof, 360 2002). In those circumstances, we can use = 2 = 2 / , with the inverse-square law = −2 = 2 −2 , to show that the drag contribution will grow approximately linearly with frequency according to: Even though the amplitude of roughness decreases at shorter wavelengths, our results, combined with the theory, show that those short scales will be more effective at causing form drag. This is because the factor ( / * ) 3 increases faster than the inverse square law decreases. The resulting linear dependence is shown as a red line on Figure 6c. As a consequence of the linear increase in the drag contribution with frequency, the total drag would become unbounded, and the sliding would stop, unless the bathymetry becomes smooth at scales smaller than some wavelength . This wavelength provides an upper 370 bound to the spatial frequency = −1 = / . Under this assumption, the total drag can be approximated by: ≪ , so that ≫ 1, this gives the approximation: 375 Therefore, if the observed inverse-square law power spectra applies down to the finest length scale , the drag will be determined by that scale, along with the viscosity , and the coefficient of the power law recovered from the periodograms. 380 It is common in sliding theories to define the slip length, = / . This lets us make a dynamical distinction between regions of fast sliding with little internal deformation > and regions of slow sliding and shearing flow < . For the situation described above, the ratio of slip length to ice thickness is / = 2 /(8 3 ). For a value of = 0.1 m and a typical ice thickness scale of = 1 km, this suggests that features on scales smaller than = 150 m would provide sufficient drag to induce significant vertical shearing within the ice. Since features on this scale are well resolved by the bathymetric profiles 385 (e.g. Figs. 3, 4) and fall within the range of frequencies where the inverse square power law applies, we conclude that the form drag produced by the observed subglacial roughness would have produced significant shearing within the flow of the grounded ice as it retreated over the highs and ridges surveyed by the MBES. This suggests that it is important to include the effects of form drag caused by basal roughness over such terrain, and by extension over the extant parts of TG today. profile is exceptionally smooth. The spectral analysis confirms this (Figs. S4p, q), and the coefficient that best fits the observations is some two orders of magnitude below the more generally applicable value of A = 0.1 m. Repeating the above analysis with = 0.001 m shows that the power law would have to apply down to horizontal scales smaller than = 15 m.
Features on this scale are not well resolved in our bed profiles (e.g. the MBES grids have cell sizes of 50 m). This means that in contrast to the more general case, it remains possible that the MSGL terrain is so smooth that the resulting form drag 395 produces little vertical shearing within the ice that flows over it, making the ice dynamics of this area more akin to the flow described for slippery-based ice streams by MacAyeal (1989). This result is consistent with our understanding of how MSGLs form, i.e. via the self-organisation of deforming sediment at the bed under fast-flowing ice (e.g. Spagnolo et al., 2014). We also repeated the analysis in the direction perpendicular to elongated features (Figs. 7d, S4k-o, q). There is no evidence that ice flowed in this direction, but the theory can nevertheless compute the contributions to form drag that would arise in that 400 hypothetical situation. For most of the across-flow lines (Fig. S4), the power spectra are similar to the along-flow direction, and the drag contributions at each frequency are similar. This suggests that the drag coefficient over much of the surveyed terrain is not especially sensitive to the flow direction.

Implications from the new bathymetric and geomorphological data 405
Our results provide the first observation-based, high-resolution geomorphic characterisation of the coastal bathymetry at TG, a former bed for the glacier. These data allow us to investigate bathymetric controls on ocean circulation towards the modern grounding zone, as well as to identify the locations, water depths and substrate compositions of ice-shelf pinning points and former grounding zones. The dominant bathymetric features, a NNE-SSW trending trough and landward flanking discontinuous ridge (Fig. 2a), represent a subtly different morphologic terrain from highly rugged, basin-dominated areas north 410 and east of the EIS (Fig. 3a) or the moderate-relief areas of lineated terrain with fewer bathymetric highs in eastern Pine Island Bay (Fig. 2a) Arndt et al., 2018;Kirkham et al., 2019). The continuity and orientation of the trough and ridge relates to the structure of basement rocks on the inner shelf. NNE-SSW to ENE-WSW structural lineaments have been identified by previous aeromagnetic surveys (Gohl, 2012;, and gravity-derived bathymetries all resolve a broad NNE-SSW ridge coincident with H1-H3, as well as deeper troughs on either side of the ridge (Tinto & Bell, 2011;415 Millan et al., 2017;Jordan et al., 2020.). Thermochronological analyses of onshore rock samples also infer a NNE-SSW trending tectonic rift structure (Spiegel et al., 2016).
Our results highlight several key differences when compared with the most recent regional bathymetric datasets (Fig. 8). For the area of new MBES data in front of TG, gravity-derived bathymetry generally underestimates sea-floor depths (Figs. 8a, b), whereas the IBCSO bathymetry, which is based on real sea-floor soundings but relies on gravity-inversion elevations and 420 interpolation in this area (Arndt et al., 2013), generally overestimates sea-floor depths (Fig. 8c). All of these datasets fail to https://doi.org/10.5194/tc-2020-25 Preprint. Discussion started: 3 February 2020 c Author(s) 2020. CC BY 4.0 License. capture the higher-frequency topographic variability revealed by the new MBES data. Although sea-floor highs are sometimes >100 m shallower than the regional products predict, this effect is most notable for the troughs, which are in reality 100 to 550 m deeper than gravity-derived bathymetries, and 50 to 250 m deeper than in the IBCSO dataset (Fig. 8d). When we consider cross-sections of three troughs that are potential pathways for CDW to the grounding zone (T2, T3, T4; locations marked by 425 asterisks in Fig. 8a), the depth errors are up to 500 m, 400 m and 250 m from east to west. Using a conservative top-CDW depth of 500 m for the TG area (based on hydrographic data acquired during NBP19-02; B. Queste, pers. comm.; cf. Nakayama et al., 2013), we calculate the cross-sectional area that CDW occupies in these troughs from our MBES data and from the grid of Millan et al. (2017). We find that the gravity-derived bathymetry underestimates the cross-sectional areas by 9%-38% for two of the three troughs, and that trough T2 between H1 and H2 (Fig. 3a) is not resolved at all on the Millan et al. (2017) grid 430 (Table S1; Fig. S5). This discrepancy would lead to a significant underestimation of the volume of CDW that can be transported through the troughs to the ice-shelf cavity and grounding zone of TG (cf. Walker et al., 2007;Arneborg et al., 2012; Table S1), which has implications for ocean circulation modelling of oceanic heat flux, ice-shelf melting, warm-water access to the grounding zone and, ultimately, projected mass loss from TG and the WAIS. The discrepancy underlines the importance of high-resolution observational datasets like MBES for capturing high-amplitude bathymetric variations at short to medium 435 wavelengths (i.e. λ <10 3 km), particularly in areas close to ice-shelf cavities and the grounding zone.
Having produced a new gravity-derived bathymetry model for the Thwaites, Crosson and Dotson-Getz area, constrained by the NBP19-02 MBES data, Jordan et al. (subm.) similarly determined that high-resolution observations are necessary in areas where knowledge of the bed at scales of less than a few kilometres is required. The need for high-resolution bathymetry was underscored by recent predictive modelling studies of Antarctic outlet glaciers, which conclude that the shape of the ice-shelf 440 cavity and knowledge of small, kilometre-scale pinning points are both key to improving predictions of ice-sheet retreat and sea-level change Favier et al., 2016). Similarly, the latest high-resolution ocean models demonstrate that warm deep water reaches the grounding zone of TG through topographically-constrained pathways, again highlighting the critical need for high-resolution bathymetry in making accurate predictions (Nakayama et al., 2019).
The geometry and detailed morphology of the H1-H3 ridge also provide insight on ice-shelf pinning points. Historical 445 grounding-line positions, as mapped from remotely-sensed ice-shelf tidal response, confirm that the Thwaites Ice Shelf is still pinned on high H1 and was pinned on high H2 as recently as 1992 and 2011 . By 2011, the area of grounding had reduced to <0.5 km 2 (Fig. 3a) but the recent configuration of the TGT suggests that at least some parts of the TGT remain in (ephemeral) contact with the sea floor. Thus, the exposed H1 and H2 sea-floor highs, and by analogy H3, can be studied as current or recent pinning points for the Thwaites Ice Shelf. Glacial lineations and (or) grounding-zone wedges 450 (GZWs) on the surface of H2 and H3, as well as rare iceberg ploughmarks (Fig. 4b), confirm that these pinning points are mantled by some amount of unconsolidated sediment that can be ploughed or moulded by ice. Sub-bottom profiles over the H2 and H3 highs support this as they show either an incredibly smooth sea-floor response, strongly indicative of unconsolidated sediment cover, or up to 10-15 m of unconsolidated sedimentary units (Fig. S6); furthermore, coring of the top of H3 recovered after grounded ice had retreated or lifted-off from the high, the presence of a GZW and on H2 (Figs. 4b, d) suggests that at least the uppermost part of this high may have been constructed via sedimentation at the grounding zone (cf. Alley et al., 1989;2007). The potential effects of this grounding-zone sedimentation are two-fold: when the TG grounding zone had retreated onto these highs during the Holocene (i.e., at some time after 10.3 cal. ka BP; Hillenbrand et al., 2013), GZW formation could have temporarily slowed further retreat (Alley et al., 2007). Second, continued pinning of an ice shelf on the high and GZW, 460 when most of the grounding line had eventually retreated further landward, would have buttressed the grounded upstream section of TG. The new MBES dataset we present here and the sea floor landforms it reveals, supported by core recovery and sub-bottom profiles, indicate that more sediment is present in this area than is typical of other Amundsen Sea inner shelf environments that experienced rapid ice-sheet retreat, including the adjacent Pine Island Bay and the Dotson-Getz palaeo-ice stream trough (e.g. Larter et al., 2007;Graham et al., 2009;Nitsche et al., 2013Nitsche et al., , 2016. This is likely because this area 465 immediately offshore TG represents a location where the grounding zone was positioned for a relatively long period of time. Our only constraints on grounding-zone retreat through this area (during the Holocene) are from the core on the H1 high which shows that grounded ice withdrawal from that high by 10.3 cal. ka BP  and grounding zones mapped from satellite-era datasets . Thus, the TG grounding zone must have been located between H1 and the current grounding zone, potentially on the sea-floor ridges identified here, for thousands of years delivering a significant 470 volume of sediment to the area.
The submarine landforms observed on and around the sea-floor highs raises the question of the composition of these features.
Landforms on the flanks of the pinning points (gullies, slide scars, isolated blocks; Fig. 4) may indicate that these highs consist, at least in part, of an erodible material with a probable (hard) bedrock core. In marine settings, slide scars and gullies incise large, pronounced sedimentary scarps like the shelf edge (e.g. Noormets et al., 2009;Gales et al., 2013) or the headwalls of 475 major submarine slides (e.g. Laberg and Vorren, 2000;Vanneste et al., 2006) but do not characterise hard bedrock (crystalline) settings. Further evidence comes from the new observation of planed-off (flat-topped) morphology of the H2 and H3 highs confirming that the upper part of these (down to the level of planing off) consists of a lithology that is erodible by the motion of an ice shelf. Other examples of flat glacial erosion surfaces planed off by ice shelves or flat-based tabular icebergs from the Arctic all document erosion into sedimentary substrates (e.g. Vogt et al., 1994;Jakobsson et al., 2010;Noormets et al., 2016). 480 Small (<100 m high, <5 km wide), flat-topped mounds in the Ross Sea are also thought to consist of unconsolidated volcanogenic deposits rather than crystalline bedrock and, interestingly, GZWs have built up on them, indicating that these features slowed grounding-zone retreat in that area (Lawver et al., 2012;Greenwood et al., 2018). Our interpretation of a proportion of unconsolidated sedimentary substrate, and thus low density material, on the H2 and H3 highs may explain why bathymetries derived from gravity over-estimate the height of some of these features (Figs. 8a, b). In contrast, the morphology 485 of H1 is consistent with it having a crystalline composition, which is supported by the small differences between the water depths of our MBES dataset and the gravity-derived datasets at the H1 location (see proximal 25 km in profile x-x' of Fig. 8d).
This very shallow, rugged feature is cross-cut by bedrock grooves and channels typical of hard rock exposures on the inner Antarctic shelf (Lowe and Anderson, 2002;Livingstone et al., 2013) and has no planed-off sections or glacial lineations on its https://doi.org/10.5194/tc-2020-25 Preprint. Discussion started: 3 February 2020 c Author(s) 2020. CC BY 4.0 License. surface (Fig. 3b). We also note the relative scarcity of bedrock channels or other landforms related to subglacial meltwater 490 flow in the TG MBES dataset with the crescentic scours being the exception. As an example, Kirkham et al. (2019) mapped more than 1000 subglacial channels in Pine Island Bay whereas we map only 175 forms here, albeit over a smaller area. It is not clear whether evidence of previous meltwater routing is buried by sediment in the deep troughs, or has been destroyed by ice flow over the highs. Physical-property and geochemical analyses on cores from the area, acquired as part of ITGC, should shed light on the frequency and magnitude of meltwater during the retreat of grounded ice over the sea-floor highs. Schroeder 495 et al. (2013) identified a transition from a distributed channel network with ponded water behind ridges at the modern grounding zone, to a system of concentrated channels downstream. It is possible that a similar configuration for the basal hydrological system occurred as ice retreated over the offshore highs and that evidence is preserved in the marine sedimentary record.
Given the apparent shaping and fragmentation of the H2 and H3 highs (Figs. 4a, b), this observation highlights another potential 500 mechanism influencing glacier dynamics, because if the substrate of a pinning point is soft enough to be moulded by an ice shelf, it might be eroded over time. Erosion of material from the surface of a pinning point, as it is planed off, in conjunction with retrogressive failures at its seaward flank and possibly larger slide or slump events may reduce its height until it cannot serve as a pinning point for the ice shelf (and glacier ice upstream) any longer. This process of unpinning would be exacerbated by increased flow velocities (leading to increased erosion) and/or by ice-shelf thinning, due to either flow acceleration or sub-505 ice shelf melting, (leading to ungrounding). The result, in a setting with erodible pinning points, is the potential for rapid grounding-zone destabilisation and resulting mass loss, whereby increased ice-flow velocities accelerate pinning point destruction which, in conjunction with simultaneous ice-shelf thinning in response to sub-ice-shelf melting, then promotes ungrounding earlier than would occur on a corresponding "hard", less erodible pinning point.

Implications from the new bed roughness data 510
One major objective of our research is to assess the deglaciated subaqueous terrain offshore from TG as an analogue for the modern bed to gain new insights on TG bed characteristics. The consistency of derived power spectra and drag contributions for bed profiles from the inner ASE shelf and for upstream areas of Pine Island and Thwaites glaciers (Fig. 7) indicates that the roughness properties of the offshore and onshore areas are comparable across all resolvable frequencies. Furthermore, observations confirm that recent grounding-zone retreat affecting TG has occurred over a series of bedrock ridges with the 515 loss of pinning points and formation of new cavities (Tinto & Bell, 2011;Milillo et al., 2019). Further upstream, about 100 km from the recent grounding zone, analyses of radar specularity suggest that the modern TG bed is characterised by high roughness attributed to bedrock cropping out at the glacier base (Schroeder et al., 2014). Our data reveal that the bedrock ridges and intervening troughs underlying the modern grounding zone (Holt et al., 2006;, with length scales of tens of kilometres and amplitudes of several hundreds of metres, constitute a morphological terrain similar to the 520 coastal bathymetry (Figs. 2, 3). Further, we demonstrate that this rugged terrain would exert the same strong influence on basal drag for an overriding ice mass (assuming no cavitation) (Figs. 7,SF4). This is consistent with results from inverse methods https://doi.org/10.5194/tc-2020-25 Preprint. Discussion started: 3 February 2020 c Author(s) 2020. CC BY 4.0 License.
to determine basal drag for the modern TG bed, which derive high basal shear stresses near the grounding zone that may explain the steep surface profile of the glacier (Joughin et al., 2009;Arthern et al., 2015). We also note that crag-and-tail landforms (which form subglacially) extend down to the floors of the deep troughs (e.g. Fig. 3, 4b). This confirms that, at least 525 at the time when these features formed, ice of an expanded TG was grounded in the troughs as well as on the highs, and probably experienced high basal shear similar to ice at the present-day grounding zone. The orientation of the crag-and-tails also confirms that ice flow was not directed along troughs but rather overrode the existing topography; this finding is consistent with cosmogenic exposure data from Bear Peninsula (Fig. 1) showing that the ice-sheet surface rose above the top of this terrain during the last glacial period (Johnson et al., 2017). 530 For shorter wavelengths of bed topography, we can consider the form of the individual sea-floor highs over length scales of several kilometres. We note that these features exhibit the same correlation of morphology with bed type that has been described from on-ice seismic reflection profiles for TG (Muto et al., 2019a, b). Specifically, the correlation is between "hard" beds on the stoss sides of topographic highs, associated with crag-and-tail landforms, and "soft" or sedimentary beds on the lee sides of these features. Although this correlation was based on data from a single along-flow and two across-flow seismic-535 reflection lines in the upstream area of the glacier, the inference is that this correlation can be applied to the bed in downstream and grounding-zone areas (Muto et al., 2019a, b). This pattern is clearly replicated over the H2 and H3 highs, which have rugged upstream ends with crag-and-tail landforms, glacial lineations over their tops, and sedimentary "tails" on their downstream ends (Fig. 4). This landform pattern is also similar to the 3D morphology of a topographic high near the grounding zone of the Rutford Ice Stream that was imaged using a grid of densely-spaced radar lines (Fig. 4 in King et al., 2016). The 540 correlation of bed types with sea-floor highs (and ridges) holds true for several other areas of the inner shelf around West Antarctica, where streamlining of bedrock highs has often produced landforms with sedimentary tails on the lee sides of bedrock obstacles (e.g. Larter et al., 2009Livingstone et al., 2013;Nitsche et al., 2013;2016), although this is not always the case. Thus, the variability in bed types on topographic highs in offshore regions may provide useful constraints on bed type variability in onshore areas. 545 Regarding the spectral analysis of roughness and basal drag contributions presented here (Figs. 6, 7, S4), we acknowledge that these only provide an order of magnitude assessment of the contribution to basal drag from the different wavelength scales resolved by the bathymetric DEM (Fig. 2). Analysis beyond the simple 2D-flow line theory used here (see Section 4.2) would be needed to account for 3D-flow effects, and for the nonlinear dependence of ice viscosity on stress (Glen, 1955). Here, we have not been specific about the physical mechanism that controls the smallest roughness wavelength ( ). Candidate 550 mechanisms that might limit the influence of roughness at small spatial scales include cavitation (Fowler, 1986), fracture and plucking of crystalline or sedimentary rocks, bulldozing of unconsolidated sediment, or regelation flow around small obstacles (Weertman, 1957). More sophisticated theories accounting for the potential of ice to form cavities in the lee of obstacles could be deployed similarly, but the drag contribution would then depend also on water pressure (Fowler, 1986;Schoof, 2005).
Process models of subglacial hydrology, phase change, fracture and sediment transport could all be incorporated in to a more 555 elaborate analysis using MBES datasets as input.
It is clear from our results that the increased spatial resolution of the MBES data is critical for capturing the high-frequency bathymetric variability on the inner continental shelf seaward of TG, which is necessary to understand warm water incursions into sub-ice shelf cavities (Figs. 8, S5;Nakayama et al., 2019). The strong correlation of our observations with interpretations of the present bed conditions of TG and, therefore, the robustness of this deglaciated terrain as an analogue for the modern 560 bed, further demonstrates that more information can be gleaned from this type of marine dataset (i.e. near-continuous bathymetry with spatial resolution better than 0.5 km) For example, the 3D-nature of MBES (with approximately equal resolution in all directions) means that bathymetric variability could be examined in any direction, not only along survey lines, as is the case with onshore radar or seismic-reflection profiles of extant bed topography, and over a variety of spatial scales.
These analyses contribute to our understanding of across-flow contributions to basal drag or hydraulic potential (e.g. Muto et 565 al., 2019a), and allow us to constrain the spatial variability of bed types (sediment vs. hard beds, thickness of deformable till, grain size distribution of subglacial sediments), particularly if the sea-floor sediments are cored for ground-truthing. Similarly, the application of theories of subglacial processes as discussed above to high-resolution bathymetric datasets will increase our understanding of ice flow over high-frequency bed roughness, particularly if combined with ultra-high resolution (sub-metre resolution) bathymetries from AUV surveys (e.g. Davies et al., 2017). Indeed, AUV surveys and (or) a dense grid of seismic 570 soundings (only obtainable from non-crevassed ice shelves) are the only way to determine bed geometry in ice-shelf cavities.
New techniques such as swath-radar that can image the present glacier bed in 3D, albeit in narrow swaths, have already been trialled on TG (Holschuh et al., 2019) and could be used in conjunction with offshore bathymetric data to build a betterinformed, more complete and more uniform resolution picture of basal conditions under TG and at its grounding zone.

Conclusion 575
New 3D bathymetric data from just offshore Thwaites Glacier reveals that the coastal bathymetry is dominated by a ~65 kmlong, ~1200 m-deep trough and discontinuous ridge with water depths of 600 m to <100 m. Spatial variations in the morphology of the ridge segments/highs suggest differences in substrate composition along the ridge, with the two southernmost highs having a significant "erodible" component, which is presumably sedimentary in composition. The geometry (flat tops) and landform evidence (glacial lineations, gullies, sediment fans) indicate that the bathymetric highs were 580 planed off and variously eroded by the action of Thwaites Ice Shelf as it flowed over them, presumably reducing the height of these former pinning points over time. A feedback mechanism during unpinning may have occurred, whereby as the ice shelf started to lose contact with the high and frontal buttressing weakened, the resultant increase in flow velocities exacerbated erosion of the high and facilitated further unpinning of Thwaites Glacier.
We present three lines of evidence that this coastal bathymetry provides a good analogue for the modern grounding zone of 585 Thwaites Glacier. First, on length scales of several tens of km the ridge and trough morphology is consistent with the bed topography of the grounding-zone area based on available DEMs and over-ice geophysical data. Second, our spectral decomposition of roughness and basal drag over this rugged, deglaciated terrain is consistent with similar spectral https://doi.org/10.5194/tc-2020-25 Preprint. Discussion started: 3 February 2020 c Author(s) 2020. CC BY 4.0 License. decompositions, and inversions of basal drag, for profiles from the modern grounding zone area and for areas of the Thwaites bed, where bedrock crops out subglacially (e.g. Schroeder et al., 2014). In contrast, smooth beds, characterised by thick 590 sedimentary substrates and linear glacial landforms, produce distinctly different power spectra and drag contributions. Third, the distribution of landforms and substrate types (unconsolidated sediment vs. bedrock) over the ridge indicates that it displays the same correlation of bed type with topography that has been described for upstream bed areas and inferred for the grounding zone (Muto et al., 2019a, b). As such, further analyses of this deglaciated terrain may provide realistic constraints on acrossflow roughness and bed type distribution, and should inform geophysical observations of the modern TG bed that will be 595 acquired as part of ITGC.
As discussed above, observational datasets like MBES are required seaward of Antarctic ice shelves in order to capture the high-frequency variability that characterises the bathymetry of nearshore areas. Modelled bathymetry cannot adequately reproduce the km-to sub-km-scale features that are important for accurately calculating inflows of warm ocean water in troughs, and for defining the topographic highs that may act as pinning points for ice shelves and as barriers to warm water 600 incursions.
Author contributions. KAH, RDL, RA, TAJ, AGCG and FON developed the concept of the paper. KAH, RDL, AGCG, RTM, JDK, RC and VF acquired and performed initial processing on the bathymetry data during NBP19-02; JDK performed all channel metric analyses. KG and JEA, and JH provided bathymetric data from German and Korean cruises, respectively. KAH 605 compiled, part-processed and gridded all bathymetric datasets, and wrote the first draft of the paper with substantial contributions from RA and RDL. RA developed the methodology for, performed spectral analyses and basal drag estimations, and wrote the text for these sections; KAH and TAJ provided the profile data. All authors contributed to the development of the final paper and data visualisation.

610
Competing interests. The authors declare no competing interests.    (Table 1). The dark blue coastline illustrates the ice-shelf and ice-mélange extent during survey on NBP19-02 and was digitised from Landsat 8 imagery.