# Evidence for an extensive ice shelf in northern Baffin Bay during the Last Glacial Maximum

Pierre-Olivier Couette, Patrick Lajeunesse, Jean-François Ghienne, Boris

Dorschel, Catalina Gebhardt, Dierk Hebbeln, Etienne Brouard

► To cite this version:

Pierre-Olivier Couette, Patrick Lajeunesse, Jean-François Ghienne, Boris Dorschel, Catalina Gebhardt, et al.. Evidence for an extensive ice shelf in northern Baffin Bay during the Last Glacial Maximum. Communications Earth & Environment, 2022, 3 (1), pp.225. 10.1038/s43247-022-00559-7. hal-03806433

HAL Id: hal-03806433

https://hal.science/hal-03806433

Submitted on 8 Oct 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Evidence for an extensive ice shelf in northern Baffin Bay during the Last Glacial Maximum

Pierre-Olivier Couette<fnref n="1" />,²✉, Patrick Lajeunesse<fnref n="1" />, Jean-François Ghienne☐², Boris Dorschel☐³, Catalina Gebhardt³, Dierk Hebbeln☐⁴ & Etienne Brouard⁵,⁶

The glaciological significance of ice shelves is relatively well established for the stability of modern ice sheets of Antarctica. Past ice shelves of the Arctic, however, are poorly documented while their role for the stability of former ice sheets remains mostly unknown. Here we present swath bathymetry data and seismostratigraphic profiles that reveal a large moraine system extending along the continental slope off Baffin Island, demonstrating that a 500-m thick ice shelf covered northern Baffin Bay during the last glacial episode. We suggest that this ice shelf had a profound impact on the stability of a series of major ice streams that drained the interior of the Laurentide, Innuitian and Greenland ice sheets. Climate warming and global sea-level rise in the early stage of deglaciation possibly contributed to a large-scale break-up of the ice shelf, which led to the destabilisation and reorganisation of tributary ice streams from these three ice sheets.

<footnote n="1">¹Département de Géographie, Université Laval, Québec, Canada. ²Institut Terre Environnement de Strasbourg (ITES), UMR 7063, CNRS—Université de Strasbourg, Strasbourg, France. ³Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung (AWI), Bremerhaven, Germany. ⁴MARUM, Center for Marine Environmental Sciences, University of Bremen, Bremen, Germany. ⁵Département des Sciences de la Terre et de l'Atmosphère & GEOTOP, Université du Québec à Montréal, Montréal, Canada. ⁶Geological Survey of Canada, Natural Resources Canada, Ottawa, Canada. ✉email: pierre-olivier.couette.1@ulaval.ca</footnote>
Ice shelves are critical components of marine-terminating ice sheets due to their buttressing effect that controls mass balance by stabilizing the ice margins facing the ocean¹⁻⁷. Retreat and break-up of ice shelves in Antarctica²,³,⁶,⁸⁻¹¹ occur today at a time of important climate warming and may lead to a drawdown of the West Antarctic Ice Sheet, resulting in a rapid transfer of ice to the ocean and thereby contributing to global sea-level rise¹,²,¹¹,¹². Identifying former ice shelves and reconstructing their evolution through time—particularly for past ice sheets of the Northern Hemisphere where they remain poorly documented³,¹²,¹³—will provide an essential analogue for understanding the stability of modern ice sheets and their interplay with climate and global sea-level change.

Baffin Bay forms a 450-km-wide and 2000-m-deep embayment between the Canadian Arctic Archipelago and Greenland (Fig. 1). At the Last Glacial Maximum (LGM; locally ~25–16 ka BP), Baffin Bay was located at the confluence of three ice sheets—the Laurentide (LIS), Innuitian (IIS) and Greenland (GrIS) ice sheets¹⁴⁻¹⁷—forming a continuous belt of ice streams draining the interior of Northern North America and Greenland and extending at or near the continental shelf break in many sectors¹⁸⁻²⁶. However, the extension of these ice sheets beyond their grounding zones to form floating ice shelves remains elusive. The scenario of ice shelves in the northern Hemisphere have been hypothesized repeatedly since it was first introduced by Mercer²⁷ as a comparison between the Arctic Ocean and West Antarctica. Hughes et al.²⁸ later expanded the idea to include floating ice shelves in Baffin Bay and eastern North America. The scenario of an ice shelf extending in Baffin Bay from a southern source of ice grounded in Davis Strait and associated with the Hudson Strait Ice Stream was speculated by Hulbe²⁹ as a starting point to explain the mechanisms behind Heinrich events. Although Hulbe et al.³⁰ modified their original stance to instead support fringing ice shelves along the coasts of Eastern Canada, the idea of an

extensive ice shelf sealing Baffin Bay from the Labrador Sea was adopted in the following decades³¹,³², but this is still debated as a trigger mechanism for Heinrich events³³. However, this scenario was recently refuted by Jennings et al.³⁴ who argued that such an ice shelf would have impeded Labrador Sea water advection, a key component for biological productivity observed at the LGM in cores off western Greenland. They suggested instead that the GrIS margin was protected by either a belt of fringing ice shelves and/or perennial sea-ice—a scenario in line with Hulbe et al.³⁰—that prevented deposition of ice-rafted debris (IRD) and enabled marine advection and food supply during the LGM. Alternatively, a scenario of an ice shelf from a northern source of ice grounded in northern Baffin Bay has been so far unaddressed, in a context where the marine-terminating Lancaster Sound and Smith Sound ice streams captured several major tributaries from the LIS, IIS and GrIS¹⁷,²¹, creating ideal conditions for the formation of an ice shelf that potentially controlled the stability of other ice streams outflowing from the eastern Canadian Arctic Archipelago and western Greenland¹⁷,²⁴,²⁵.

Solving the ice shelf issue in Baffin Bay is therefore necessary for evaluating the effects of buttressing on ice flow, ice streams stability and ice sheet mass balance changes during a key stage of deglaciation of the Arctic and a period of major meltwater flow into the ocean. Here we present geophysical (swath bathymetry, acoustic profiles) datasets and a compilation of previously published sediment core data collected on the outer sector of the cross-shelf trough off Clyde Inlet (Baffin Island) and its adjacent continental slope of western Baffin Bay that together demonstrate the existence of an ice shelf in northern Baffin Bay—the Northern Baffin Bay Ice Shelf (NBBIS)—during the last glacial episode (Fig. 1). This geomorphological evidence and the reinterpretation of sediment cores allow defining the extent of the LIS on the continental shelf at the LGM and examining the potential impact of an ice shelf collapse on the evolution of ice streams of the three marine-terminating ice sheets in Baffin Bay.

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Location map of the investigated area. Location of Clyde Trough (CT), sediment cores and maximum extent of grounded ice on western Baffin Bay (thick line), compiled and modified from previous works²¹,²⁵,²⁶ and this paper. BT Buchan Trough; ST Scott Trough; SFT Sam Ford Trough. Red rectangles represent study areas shown in greater details in Figs. 2 and 5. (Inset) Location of Baffin Bay north of the Labrador Sea and between Greenland and the Canadian Arctic Archipelago (CAA).</figcaption></figure>

## Results and discussion

The geomorphological signature of the NBBIS. Three ridges running in parallel to the continental slope off Clyde Trough between 300 and 420 m below sea-level (mbsl) are observed on the swath bathymetry imagery (Fig. 2a, b). These ridges are independent from grounding zone wedges (GZWs) constructed by outlet glaciers and ice streams emanating from Baffin Island. The central ridge is the largest; it is 30-m high, ~300-m wide and extends almost continuously for >30 km along the slope off Clyde Trough (Figs. 2a and 3). It has an asymmetric profile showing a gentler ocean-facing slope and is observed at progressively decreasing depths from north (330 mbsl) to south (300 mbsl). Although it appears to be streamlined along the continental slope, one of the segments of the central ridge (Fig. 2b and c) shows evidence of sediment accretion by ice push in a NW-SE direction (i.e., towards Baffin Island). In addition, two 5–10 m-high subdue ridges occur on both sides of the main ridge at 420 and 310 mbsl, extending for <5 km in front of Clyde Trough and showing geomorphic similarities to the neighbouring central ridge (Fig. 2b). Similar ridges are also observed at 350 mbsl off Sam Ford Trough and 410 mbsl south of Scott Trough, 40 and 100 km northwest of Clyde Trough, respectively (Fig. 4). These ridges probably correspond to the central ridge observed on the upper continental slope off Clyde Trough as they are also ~30 m high with a gentler ocean-facing slope.

We rule out the possibility of these landforms being terminal moraines of the LIS on the upper continental slope as its geometry does not appear to be influenced by the location of cross-shelf troughs. If they indeed originated from flowing ice
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Geophysical data along the Baffin Shelf margin off Clyde Trough. **a** Bathymetry of the outer Clyde Trough and adjacent upper continental slope (see Fig. 1). Dashed lines represent GZW1 and GZW2 positions. **b** Bathymetry along the continental slope of off Clyde Inlet showing the ice shelf lateral moraines and the deep iceberg ploughmarks. (Inset in **b**) Bathymetric profile across the ice shelf lateral moraine. **c** Zoom on the central ridge of the ice shelf lateral moraine. **d** Parasound profile across the ice shelf lateral moraine.</figcaption></figure>

across the Baffin Island continental shelf, a bulge should be expected in front of the trough. Instead, we observe a curvature in the orientation of these ridges into Clyde Trough (Fig. 2a). These ridges are therefore interpreted as lateral moraines of the NBBIS on the basis of (a) their geomorphology similar to previously reported ice-shelf moraines in the Canadian Arctic Archipelago³⁵⁻⁴⁰, and (b) their location, depth and orientation

along the continental slope off Baffin Island. The progressively decreasing depths of these ice shelf lateral moraines comply with a paleoglaciological reconstruction including a seaward thinning of the ice shelf. The smaller ridges could represent different phases and/or thickness changes of the NBBIS during the LGM, although their relatively subdued appearance could also suggest their formation during prior glaciations (i.e., MIS-6).
<figure><img src="image_5.png" /><figcaption>**Fig. 3** The ice shelf lateral moraines on the continental slope off Clyde Trough. 3D imagery generated using the QPS Fledermaus software showing the ice shelf lateral moraines and iceberg ploughmarks.</figcaption></figure>

In outer Clyde Trough (Figs. 2a and 3), the absence of gullies on the trough-mouth fan and of streamlined glacial bedforms at the continental shelf edge suggest that the Clyde Ice Stream (CIS) did not reach the continental shelf break at the LGM—in contrast with other ice streams on northeastern Baffin Island²⁵—and was probably not connected to the NBBIS. This disconnection between the CIS and the NBBIS in outer Clyde Trough is supported by the presence of two GZWs located ~10 km landward from the continental shelf edge and ~15 km west of the slope ridge system (Figs. 2a and 5), leaving an open area between the CIS margin in the outer trough and the grounded NBBIS on the slope.

More to the north, a GZW located ~20 km off Scott Trough at ~450 mbl and associated with streamlined bedforms oriented parallel to the continental slope indicate grounded ice flowing southeastward along the shelf (Figs. 4a and b). The streamlined bedforms are similar to mega-scale glacial lineations⁴¹⁻⁴³ or sub-ice shelf keel scours⁶,⁴⁴ and were recognized to be the product of ploughing across the seafloor by an ice shelf. This system therefore suggests that ice flowing out of the Scott Trough was either (1) deflected from its along-trough original flow by an obstacle (i.e., a floating ice shelf in Baffin Bay); or (2) a large independent ice body expanding all the way up to Clyde Trough (i.e., a fringing ice shelf). The first hypothesis is more probable as there is no explanation for the divergence of the Scott Ice Stream from an along-trough to along-slope orientation without a buttressing mechanism. These landforms, together with the ice shelf lateral moraines observed at shallower depths than the bottom of Scott Trough (>600 mbl; Supplementary Fig. 1), suggest that, in contrast to Clyde Trough, the Scott Ice Stream— and probably other systems of northeastern Baffin Island— merged with the NBBIS. Mega-scale glacial lineations, ice stream lateral moraine and ‘lift-off’ moraines identified by Brouard and Lajeunesse²⁵ on the Scott trough-mouth fan in turn indicate an

ice streaming phase extending to the shelf edge (Fig. 4a). Initially interpreted as features of LGM ages²⁵, they could instead correspond to a local subsequent readvance prior to 14.2 kyr BP⁴⁵.

In addition, the continental slope along Baffin Island is heavily affected by curvilinear scours between 5 and 10 m deep and down to >750 mbl (Figs. 2a, 3 and 4a). These scours, interpreted as iceberg ploughmarks, indicate extensive calving from a thicker body of ice in northern Baffin Bay probably after the LGM¹⁵.

The available swath bathymetry and seismostratigraphic data do not allow identifying ice shelf lateral moraines south of Clyde Trough with confidence and, therefore, precisely addressing the NBBIS extent toward southern Baffin Bay. Faint bedforms similar to those presented above are observed on the bathymetry some 120 km to the south, but their association with those off Clyde Trough is ambiguous as they are located at greater depths (~400 mbl; Supplementary Fig. 2). The timing of formation of this ice shelf lateral moraine is also unknown as no sediment core, thus no direct dating, were retrieved on this bedform. In this context, the presence of an extensive ice shelf on the Arctic ocean during MIS-6⁴²,⁴⁶ could suggest contemporaneity on Baffin Bay. Such an ice shelf is supported by a raised marine sediment sequence identified on Eastern Baffin Island⁴⁷,⁴⁸ that was dated at this period⁴⁹, suggesting a more extensive ice cover during the penultimate glaciation. In addition, several studies provided evidence for the presence of the LIS along the eastern margin of Baffin Island throughout MIS-4 to MIS-2⁵⁰⁻⁵⁵. However, sediment records in Baffin Bay and along the Baffin Island continental shelf allow discussing the potential timing and extent of the NBBIS.

**The sedimentary signature of the NBBIS.** The previously published sediment records from Baffin Bay (Fig. 6)¹⁴,²¹,⁴⁵,⁵⁶⁻⁶⁰ are compatible with the presence of an ice shelf during the LGM. The
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Geophysical data along the northwestern Baffin Shelf margin. a Bathymetry of the upper continental slope between Scott Trough and Clyde Trough (see Fig. 1). b Bathymetry along the continental slope of Baffin Island showing a GZW south of Scott Trough aligned with the ice shelf lateral moraine and sub-ice shelf keel scours. (Inset in b) bathymetric profile of the ice shelf lateral moraine near Sam Ford Trough. c Airgun profile 78029_AG_274_1516 showing the ice shelf lateral moraine near Scott Trough (location in Fig. 4b).</figcaption></figure>

buildup of ice sheet prior to the LGM is marked by a layer of IRD and carbonate-rich sediments originated from northern Baffin Bay termed Baffin Bay Detrital Carbonate (BBDC) layer, which is identified basin-wide and dated between 26 and 24.7 kyr BP

(BBDC-2<sup>45,56,57,60</sup>). The quasi-absence of IRD throughout the LGM (~24.7–~16.5 kyr BP) in most cores of Baffin Bay suggests an extensive ice cover, either as perennial sea ice<sup>61</sup> or an ice shelf<sup>6,34</sup>. The occasional clasts observed in some cores of northern
<figure><img src="image_7.png" /><figcaption>Fig. 5 Seismic evidence of ice extent in Clyde Trough. Airgun profile 80028_AG_RAYT_257_0200 showing GZW1 and deeply buried glacigenic debris-flows in outer Clyde Trough (location in Fig. 2a).</figcaption></figure>

Baffin Bay could be explained by their relatively short distance from the ice grounding line as it is affected by mass flows and the rain out of ice shelf basal debris⁶,³³,⁶²⁻⁶⁴. Provenance study on core 16PC, located 100 km east of Clyde Trough in central Baffin Bay (Fig. 1), shows uniform mineralogic composition dominated by a Greenland-sourced sediment assemblage and a restricted Baffin Island-sourced input during the LGM interval⁵⁷. Localized and uniform mineralogic assemblages are typical of sub-ice shelf setting reflecting uniform conditions in which transport is dominated by meltwater plumes⁶,⁶⁵. A coarsening-upward trend with frequent IRD in cores of Baffin Bay starts at ~16.5 kyr BP⁵⁷,⁵⁸,⁶⁰ and is followed by the deposition of carbonate-rich BBDC-1 between ~14.5 and ~13.7 kyr BP⁵⁷,⁶⁰. Such a coarsening-upward sequence typifies the progressively more open water conditions that usually follow ice shelf collapse events⁶,⁶⁴. This sedimentary record is therefore in support to the hypothesis of an ice shelf in Baffin Bay during the LGM. In addition, it provides a temporal framework suggesting that break-up started at ~16.5 kyr BP and ended with the onset of BBDC-1 at ~14.5 kyr BP¹⁵,⁵⁷,⁶⁰.

Sediment cores collected on the Baffin Island continental slope (Fig. 6) allow defining the position of the LIS margin on the continental shelf off Baffin Island during the LGM and early deglaciation. The base of most cores is marked by ice-proximal dark grey-brown diamicton interpreted as glacigenic debris-flow deposits and correlated to an ice advance at ~25 kyr BP⁴⁵. In 77PC, a gold brown mud layer located halfway along the core was correlated to BBDC-2²⁶,⁴⁵ and therefore correspond to the ice advance observed in the other cores along the continental slope. Provenance study on this core indicates that a localized source of carbonate sediments, originating from a Paleozoic outcrop east of the Home Bay area⁶⁶,⁶⁷, also contributed to the composition of the gold brown mud layer⁶⁸. Cores from the continental slope offshore Baffin Island (Fig. 1) typically consist of laminated red brown mud with rare or absent IRD during the LGM interval. Jenner et al.⁴⁵ noted, however, that coarse IRD-rich beds from the basal diamicton gradually decrease in abundance upwards into laminated red brown mud and very fine sand with IRD in cores 61PC off Scott Trough. These laminated red brown muds are interpreted as lateral sediment supply from NE Baffin Island, as they thin in deeper water off the trough mouth fans¹⁴,⁴⁵, an interpretation corroborated by their mineral composition⁶⁰. They are, however, absent in core 73PC (Fig. 1), suggesting the absence of grounded ice at the shelf break in Clyde Trough during the LGM, in line with the receded position of GZWs and the presence of predominantly Greenland-sourced sediments in core 16PC⁵⁸. Except for 77PC, located off Home Bay in southern Baffin Bay, the laminated red brown muds are characterized by an absence of bioturbation and microfossils⁴⁵. According to the ice-shelf

sedimentary model of Smith et al.⁶, this sequence might represent an ice advance to the continental shelf edge at the LGM (stratified diamicton), followed by ice shelf buildup with frequent iceberg calving (turbidites interbedded with IRD layers) and, ultimately, a full ice shelf cover (laminated red brown mud). A second diamicton overlying the red brown laminated muds is observed immediately below BBDC-1 in cores north of Clyde Trough⁴⁵,⁶⁹, representing glacigenic debris-flow deposits correlated to a smaller localized advance from northeastern Baffin Island ice streams to the shelf break prior to 14.2 kyr BP⁴⁵.

To our knowledge, no sediment core has been collected or published from the Greenland side of northern Baffin Bay to allow direct comparison with Baffin Island continental slope assemblages. In southern Baffin Bay, sediment core data show similar sequence as those from the Baffin Island continental slope³⁴,⁷⁰,⁷¹. However, widespread foraminifera assemblages and bioturbation during the LGM interval suggest the presence of perennial sea-ice⁶¹ and/or a nearby ice shelf margin⁶²,⁷²,⁷³. A thorough discussion on sediment assemblages and their relation to ice shelf and/or perennial sea-ice on the continental slope of western Greenland was provided by Jennings et al.³⁴.

**Evidence of ice shelf inception in northern Baffin Bay.** The presence of an ice shelf in northern Baffin Bay at the LGM and during the early stage of deglaciation is supported by invoking the Marine Ice Cliff Instability (MICI) theory, which stipulates that fracturing due to longitudinal stresses exceeding the yield strength initiate ice cliff failures when a marine-based ice-margin reaches ~1 km in thickness⁷⁴⁻⁷⁶. Till wedges, GZWs and mega-scale glacial lineations identified on the upper continental slope off Lancaster Sound²¹,²⁵ and Smith Sound (Supplementary Fig. 2) indicate grounded ice as deep as 1350 mbl. Accounting for a coeval higher relative sea-level (RSL) of ~75 m based on the modelled RSL values from a nearby Baffin Island site⁷⁷, ice thickness can be estimated using:

$$H_{\text{LGM}} = (\rho_w/\rho_i) \cdot (D + \Delta D_{\text{LGM}}) \quad (1)$$

where $H_{\text{LGM}}$ is the minimum ice thickness during the LGM, $\rho_w = 1028 \text{ kg m}^{-3}$ and $\rho_i = 910 \text{ kg m}^{-3}$ are the mean densities of sea water and ice, respectively, D is the water depth and $\Delta D_{\text{LGM}}$ is the correction for water depth at the LGM⁷⁴. In the case of the Lancaster Sound and Smith Sound trough-mouth fans, grounded ice had to be >1550-m thick. The Lancaster Sound till wedges and the Smith Sound GZW could not have been deposited without an ice shelf buttressing—and temporarily stabilising—the front of these ice streams⁴,⁷⁶,⁷⁸. In order to maintain a stable position without an ice shelf in northern Baffin Bay, the ice margin would have needed to retreat until it reached the theoretical depth of
<figure><img src="image_8.png" /><figcaption>**Fig. 6 Sediment assemblages in Baffin Bay and along the Baffin Shelf.** Lithology of piston cores from Baffin Bay modified from various published datasets (see 'Methods'). Radiocarbon ages are reinterpreted in calendar years using the online software Calib 8.2 with the Marine20 radiocarbon age calibration curve<sup>87</sup> and a local reservoir correction (ΔR) of 87 ± 20<sup>88</sup>. In addition, relative paleointensity (*) chronology is also presented for core 16PC<sup>59</sup> and 77PC<sup>26</sup>. Correlations are based on Baffin Bay detrital carbonate (BBDC) beds and are compiled from original publications. BBDC events are represented by the brown shaded areas—ages are from Simon et al.<sup>57</sup> and Jackson et al.<sup>60</sup>; BBDC-2: -25–24.7 kyr BP; BBDC-1: -14.2–13.7 kyr BP; BBDC-0: -12.7–11 kyr BP. The grey shaded horizon represents ice shelf setting on Baffin Bay; note the quasi-absence of IRD during that interval.</figcaption></figure>
~800 m, hence continually calving icebergs throughout the LGM⁷⁶. Such a constant calving in northern Baffin Bay is, however, not compatible with the sedimentary record from the slope and deep basin off Baffin Island (Fig. 6) where a quasi-absence of IRD was observed during this interval. Moreover, numerical models indicate that ice grounded at 1500 mbls is expected to form a ~500-m-thick ice shelf⁷⁹; a thickness consistent with observations off Clyde Trough where ice shelf lateral moraines are observed at ~400 mbls (i.e., ~475 m at the LGM).

**Implications of a NBBIS.** Glacial landforms and deposits off the Baffin Island continental shelf combined with the Baffin Bay sedimentary record provide definitive evidence for the presence of a 500-m-thick ice shelf in northern Baffin Bay at and shortly after the LGM. The scenario of an ice shelf fed by ice streams located in northern Baffin Bay does not require grounded ice in Davis Strait sealing inflow of food supply from Labrador Sea waters, therefore reconciling with biological productivity observations on the continental slopes of southern Baffin Bay during the LGM³⁴,⁴⁵. This scenario provides a framework for revising the glaciation development around Baffin Bay and reassessing the potential impacts of the ice shelf collapse on the reorganisation of major ice streams draining the LIS, IIS and GrIS during the last glacial cycle (Fig. 7).

**Ice shelf inception at the LGM.** Based on the presence of BBDC-2 and diamicton along the margin of eastern Baffin Island, the NBBIS probably started building up with an advancing ice shelf margin across Baffin Bay between ~26 and ~24.7 kyr BP. The presence of the short-lived BBDC-3 event⁴⁴ suggests an earlier buildup with frequent break-up of the NBBIS margin. This earlier buildup would be in line with the presence of the LIS on the Baffin Island continental shelf—and possibly Lancaster Sound—during MIS-3⁵⁰⁻⁵⁵. An extensive LIS during MIS-3 would likely have facilitated an earlier advance of the ice margin across the other shelves surrounding Baffin Bay. Regardless of the exact timing of ice shelf inception in northern Baffin Bay, grounded ice flowing from Lancaster Sound and Smith Sound fed the NBBIS between at least ~24.7 and ~16.5 kyr BP, which in turn buttressed peripheral glaciers and ice streams of the eastern Canadian Arctic Archipelago and possibly Western Greenland. The available geological record allows identifying the southward extent of the NBBIS at least down to the latitudes of Clyde Trough. The relatively shallow depths of the ice shelf lateral moraine in front of Clyde Inlet (<350 m) suggests that the ice shelf was thinning and its margin was possibly located nearby, therefore making it difficult to extend farther south to Home Bay and preserve a constant thickness of ~500 m for >100 km. Although geophysical data from the Greenland continental shelf do not allow confirmation, it is probable that the ice shelf had a similar extent across the entire Baffin Bay. The relatively deep troughs (<500 mbls) on the Greenland side of northern Baffin Bay possibly prevented the formation of ice shelf lateral moraines. Similarly to northeastern Baffin Island, ice streams located off the GrIS may have merged with the NBBIS. The exact extent of the NBBIS in southern Baffin Bay remains unknown; sediment assemblages⁴⁵,⁶⁰ show similarities to those of northern Baffin Bay, while the record of marine productivity during the LGM suggests openings in sea-ice as polynias and/or advection under a proximal ice-shelf margin³⁴. It is therefore possible that, while northern Baffin Bay was covered by a basin-wide ice shelf system, southern Baffin Bay was characterized by fringing ice shelf and/or perennial sea-ice along both continental shelves, which could also explain the presence of ice shelf lateral moraines off Home Bay.

**Ice shelf break-up and collapse.** Rising sea level and/or enhanced incursion of relatively warm Atlantic waters via the Western Greenland Current³⁴,⁸⁰ may have provoked a destabilization and collapse of the NBBIS front between ~16.5 and ~14.2 kyr BP (Supplementary Fig. 3). The initiation of this break-up is recorded by an increase in IRD and coarse-grained sediments in cores from Baffin Bay⁵⁷,⁶⁰. Such a break-up would have been accelerated by warming summer temperatures⁸¹ providing the necessary heat to produce sufficient meltwater to fill crevasses and induce hydrofracturing¹⁰,¹¹, in turn favouring the production of small icebergs rather than large tabular icebergs¹,⁴,³⁰. The final collapse of the NBBIS corresponds to the onset of widespread delivery of northern Baffin Bay-sourced detrital carbonates in Baffin Bay—BBDC-1¹⁵,⁵⁸,⁶⁰.

**Ice streams acceleration.** As the NBBIS rapidly collapsed by ~14.2 kyr BP, ice streams of northeastern Baffin Island became unbuttressed, thereby thinning due to their acceleration, similarly to observations made on present-day ice shelves in Antarctica¹,⁵,⁷,⁸. The accelerated discharge from northeastern Baffin Island ice streams probably led to the deposition of the younger diamicton identified on the continental slope north of Clyde Trough⁴⁵. The MSGLs and ice stream lateral moraines at the shelf edge in Buchan and Scott troughs²⁵ further supports this acceleration of ice streams. The presence of 'lift-off' moraines at the shelf break²⁵ suggests tidally related buoyancy at the grounded ice sheet margin⁸², which could potentially be a sign of ice stream thinning. The absence of the second diamicton south of Clyde Trough suggests that no readvance nor glacier acceleration occured in Clyde Trough, probably because the CIS was not initially buttressed by the NBBIS. The occurrence of ridges similar to 'lift-off' moraines on Melville Bay trough-mouth fans²⁴ could hint to a similar history as for the northeastern Baffin Island ice streams. Ice streams acceleration around northern Baffin Bay are coeval to a period marked by exceptionally high rates of global sea-level rise (Meltwater Pulse 1A¹²,⁸³,⁸⁴) suggesting contribution of the North America High Arctic ice sheets to this event.

**Ice retreat in cross-shelf troughs.** Ice stream fronts retreated landward within their troughs at ~14.2 kyr BP⁴⁵, consistent with previous age constraints of ice retreat from the Baffin Island continental shelf¹⁹,²⁰,⁴⁵. The previous event of ice purge may have resulted in their episodic retreat in northeastern Baffin Island due to their destabilisation²⁵, whereas the initially unbuttressed CIS margin underwent a slower retreat. Furthermore, ice retreat in Lancaster Sound and Smith Sound led to the increased delivery of IRD-rich carbonates and deposition of BBDC-1 in Baffin Bay¹⁵,⁴⁵,⁵⁶⁻⁵⁸,⁶⁰.

## Conclusion

Identifying landforms associated with modern ice shelves remains a challenge due to their scarcity, the difficulty to access sub-ice-shelf environments and technological limitations of sonar imaging systems. However, the last decades witnessed improvements in ship-bound technologies (i.e., high-resolution multibeam echosounder, autonomous underwater vehicle and 2D/3D seismic imagery) allowing the identification of submarine bedforms diagnostic of past ice shelves, such as sub-ice-shelf keel scours/ mega-scale glacial lineations⁴¹⁻⁴³, GZWs⁴² and corrugation ridges⁴⁴. A more comprehensive understanding of such landforms is of particular interest as they provide direct and key evidence for past ice shelves in the Arctic, which may have considerable implications on glacial dynamics at the margin of continental ice sheets that was so far little studied.
<figure><img src="image_9.png" /><figcaption>**Fig. 7 Conceptual model for early deglaciation including the NBBIS and the LIS outlets in Baffin Bay.** Reconstruction of ice sheets in Baffin Bay during the LGM and subsequent deglaciation. LIS Laurentide Ice Sheet; IIS Innuitian Ice Sheet; and GrIS Greenland Ice Sheet. CT Clyde Trough; ST Scott Trough. Blue arrows indicate an advancing or stable ice margin and black triangles indicate a retreating ice margin. Brown lines represent ice shelf lateral moraines along the Baffin Island continental shelf.</figcaption></figure>
The ice shelf scenario in northern Baffin Bay provides a satisfactory explanation for the landforms observed along the continental slope between Scott and Clyde troughs, as well as in Lancaster and Smith sounds. It is also in agreement with sediment assemblages identified on the continental slopes off Baffin Island and Greenland, as well as in Baffin Bay itself during the LGM. With a surface of at least 150,000 km²—an area comparable to some of the largest ice shelves of present-day Antarctica—the NBBIS constitutes an analogue from the past for understanding future changes to the cryosphere in a warming world in a context of disintegration of buttressing ice shelves. The demise of the NBBIS would have provided an ideal setting for triggering the acceleration of ice streams and the rapid loss of glacial ice to the ocean, thereby influencing drastically the ice flow regimes of the LIS, IIS and GrIS. These findings could further support the hypothesis that the North America High Arctic ice sheets did not have considerable contribution to rising sea level before Meltwater Pulse 1A⁸³,⁸⁴. Defining the extent of the NBBIS and the chronology of its history is, however, necessary for assessing a scenario that should include the acceleration of major Arctic ice streams following ice shelf collapse when identifying sources for Meltwater Pulse 1a. This paper provides the basis for further investigations aiming at documenting—both offshore and onshore—lateral ice shelf moraines of past ice shelves of the Arctic seas.

## Methods

**Geophysical data along the Baffin Island continental slope.** Swath bathymetric data was collected in 2017 during expedition MSM66 of the RV Maria S. Merian⁸⁵ using a Kingsberg Simrad EM122 (12 kHz) system, which allowed the coverage of Clyde Trough. This dataset is complemented by swath bathymetric data from ArcticNet cruises collected during transits from 2003 to 2016 onboard the CCGS Amundsen, equipped with a Kingsberg Simrad EM302 (30 kHz). These datasets were processed for anomalous data points and artefacts removal using Caris Hips and Sips software and gridded at a 10 m cell-size resolution. They were then imported into ESRI ArcGIS software for geomorphological mapping and landforms identification. The International Bathymetric Chart of the Arctic Ocean (IBCAO) data⁸⁶, gridded at a 500 m cell-size resolution, was also used to show overall bathymetry of Baffin Bay and its adjacent shelves. Parasound profiles were acquired onboard the RV Maria S. Merian using a Teledyne-Reson Parasound DS P-70 (secondary low-frequency ca. 4 kHz). The raw data were recorded into ASD and PS3 format and then converted into SEGY using ps32segy software of Dr. Hanno Keil (University of Bremen). The shallow-acoustic profiles were then imported into the IHS Markit Kingdom software for processing and interpretation. Seismic lines 78029_AG_274_1516 and 80028_AG_RAYT_257_0200 were acquired by the Geological Survey of Canada during airgun surveys in 1978 and 1980, respectively. Extraction and interpretation were done using the LizardTech GeoViewer® software. Both acoustic and seismic data were then transferred into Adobe Illustrator® for figure production. Thicknesses and water depth were calculated using a velocity of 1500 m/s.

**Sediment cores and chronology.** Sediment cores were compiled from previous studies to describe the lithostratigraphy along the Baffin Island continental slope and in Baffin Bay. Jenner et al.⁴⁵ noted that cores 64PC, 61PC, 74PC and 77PC (all included in this study) typify the regional lithostratigraphy along the northeast Baffin margin from north to south. In addition, core 73PC⁴⁵ was included in the compilation as it was located nearby Clyde Trough. Cores located directly in front of the Baffin Island troughs in which the LGM assemblages were mostly dominated by glacial debris-flows and turbidites were not included. Cores from Baffin Bay were selected to allow a north to south representativity and to avoid redundant information. In general, the lithofacies provided a good correlation between them on the basis of radiocarbon dating. For full description of the cores, readers are referred to the original publications—16PC⁵⁷; 34PC and 36PC⁵⁶; 61PC, 64PC, 73PC, 74PC and 77PC⁴⁵; 67PC²¹; SL174⁶⁰.

Radiocarbon ages were calibrated within the age-depth modelling process using the online software Calib 8.2 with the Marine20 radiocarbon age calibration curve⁸⁷ and a local reservoir correction (ΔR) of 87 ± 20⁸⁸ commonly used on northeastern Baffin Island²⁶,⁴⁵,⁵⁷,⁶⁰ (Supplementary Table 1). This reservoir correction is in line with the one currently used for radiocarbon date calibrations of central West Greenland³⁴,⁷¹, that indicate a reservoir correction (ΔR) of 140 ± 35 years⁸⁹. We recognize that an offset caused by a well-mixed ocean was most likely variable through time and the appropriate ΔR reservoir correction for the region might have been larger than the one used and provide younger ages.

However, this would not considerably change the interpretation of the ice shelf scenario.

## Data availability

The swath bathymetric data collected in 2017 during expedition MSM66⁸⁵ have been deposited in the PANGAEA repository, available under https://doi.org/10.1594/PANGAEA.902341. The swath bathymetric data from ArcticNet cruises can be visualized on the Université Laval Géindex + website (http://geoindex-plus.bibl.ulaval.ca). The Parasound profiles collected in 2017 during expedition MSM66⁸⁵ have been deposited in the PANGAEA repository, available under https://doi.pangaea.de/10.1594/PANGAEA.944843. The seismic reflection data along with the acquisition specifics are available on the Marine Data Holding public repository of National Resources Canada (http://geogratis.gc.ca/).

Received: 7 February 2022; Accepted: 22 September 2022; Published online: 03 October 2022

## References

1. Joughin, I. & Alley, R. B. Stability of the West Antarctic ice sheet in a warming world. *Nat. Geosci.* **4**, 506–513 (2011).

2. Pritchard, H. et al. Antarctic ice-sheet loss driven by basal melting of ice shelves. *Nature* **484**, 502–505 (2012).

3. DeConto, R. & Pollard, D. Contribution of Antarctica to past and future sea-level rise. *Nature* **531**, 591–597 (2016).

4. Wise, M. G., Dowdeswell, J. A., Jakobsson, M. & Larter, R. D. Evidence of marine ice-cliff instability in Pine Island Bay from iceberg-keel plough marks. *Nature* **550**, 506–510 (2017).

5. Reese, R., Gudmundsson, G. H., Levermann, A. & Winkelmann, R. The far reach of ice-shelf thinning in Antarctica. *Nat. Clim. Change* **8**, 53–57 (2018).

6. Smith, J. A. et al. The marine geological imprint of Antarctic ice shelves. *Nat. Commun.* **10**, 1–16 (2019).

7. Bart, P. J. & Tulczyk, S. A significant acceleration of ice volume discharge preceded a major retreat of a West Antarctic paleo–ice stream. *Geology* **48**, 313–317 (2020).

8. Scambos, T. A., Bohlander, J. A., Shuman, C. A. & Skvarca, P. Glacier acceleration and thinning after ice shelf collapse in the Larsen B embayment, Antarctica. *Geophys. Res. Lett.* **31**, L18402 (2004).

9. Rignot, E., Mouginot, J., Morlighem, M., Seroussi, H. & Scheuchl, B. Widespread, rapid grounding line retreat of Pine Island, Thwaites, Smith, and Kohler glaciers, West Antarctica, from 1992 to 2011. *Geophys. Res. Lett.* **41**, 3502–3509 (2014).

10. Lai, C. Y. et al. Vulnerability of Antarctica's ice shelves to meltwater-driven fracture. *Nature* **584**, 574–578 (2020).

11. Gilbert, E. & Kittel, C. Surface melt and runoff on Antarctic ice shelves at 1.5 °C, 2 °C, and 4 °C of future warming. *Geophys. Res. Lett.* **48**, e2020GL091733 (2021).

12. Lin, Y. et al. A reconciled solution of Meltwater Pulse 1A sources using sea-level fingerprinting. *Nat. Commun.* **12**, 1–11 (2021).

13. Briner, J. P. et al. Rate of mass loss from the Greenland Ice Sheet will exceed Holocene values this century. *Nature* **586**, 70–74 (2020).

14. Aksu, A. E. & Piper, D. J. Late quaternary sedimentation in Baffin Bay. *Can. J. Earth Sci.* **24**, 1833–1846 (1987).

15. Andrews, J. T., Kirby, M. E., Aksu, A., Barber, D. C. & Meese, D. Late Quaternary detrital carbonate (DC-) layers in Baffin Bay marine sediments (67°–74°N), correlation with Heinrich events in the North Atlantic?. *Quat. Sci. Rev.* **17**, 1125–1137 (1998).

16. Dyke, A. S. et al. The Laurentide and Innuitian ice sheets during the Last Glacial Maximum. *Quat. Sci. Rev.* **21**, 9–31 (2002).

17. Margold, M., Stokes, C. R. & Clark, C. D. Reconciling records of ice streaming and ice margin retreat to produce a palaeogeographic reconstruction of the deglaciation of the Laurentide Ice Sheet. *Quat. Sci. Rev.* **189**, 1–30 (2018).

18. Briner, J. P., Miller, G. H., Davis, P. T. & Finkel, R. C. Cosmogenic radionuclides from fiord landscapes support differential erosion by overriding ice sheets. *Geol. Soc. Am. Bullet.* **118**, 406–420 (2006).

19. Briner, J. P., Overeem, I., Miller, G. H. & Finkel, R. C. The deglaciation of Clyde Inlet, northeastern Baffin Island, Arctic Canada. *J. Quat. Sci.* **22**, 223–232 (2007).

20. Praeg. D. B., MacLean, B., & Sonnichsen, G. Quaternary Geology of the Northeast Baffin Island Continental Shelf, Cape Aston to Buchan Gulf (70° to 72°N) (Geological Survey of Canada, 2007).

21. Li, G., Piper, D. J. & Calvin Campbell, D. The Quaternary Lancaster Sound trough-mouth fan, NW Baffin Bay. *J. Quat. Sci.* **26**, 511–522 (2011).
22. Ó Cofaigh, C. et al. An extensive and dynamic ice sheet on the West Greenland shelf during the last glacial cycle. *Geology* **4**, 219–222 (2013).

23. Dowdeswell, J. A. et al. Late Quaternary ice flow in a West Greenland fjord and cross-shelf trough system: submarine landforms from Rink Isbrae to Uummannaq shelf and slope. *Quat. Sci. Rev.* **92**, 292–309 (2014).

24. Slabon, P. et al. Greenland ice sheet retreat history in the northeast Baffin Bay based on high-resolution bathymetry. *Quat. Sci. Rev.* **154**, 182–198 (2016).

25. Brouard, E. & Lajeunesse, P. Maximum extent and decay of the Laurentide Ice Sheet in Western Baffin Bay during the Last Glacial episode. *Sci. Rep.* **7**, 10711 (2017).

26. Lévesque, Y., St-Onge, G., Lajeunesse, P., Desiage, P. A. & Brouard, E. Defining the maximum extent of the Laurentide Ice Sheet in Home Bay (eastern Arctic Canada) during the Last Glacial episode. *Boreas* **49**, 52–70 (2020).

27. Mercer, J. H. A former ice sheet in the Arctic Ocean? *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **8**, 19–27 (1970).

28. Hughes, T. J., Denton, G. H. & Grosswald, M. G. Was there a late-Würm Arctic ice sheet. *Nature* **266**, 596–602 (1977).

29. Hulbe, C. L. An ice shelf mechanism for Heinrich layer production. *Paleoceanography* **12**, 711–717 (1997).

30. Hulbe, C. L., MacAyeal, D. R., Denton, G. H., Kleman, J. & Lowell, T. V. Catastrophic ice shelf breakup as the source of Heinrich event icebergs. *Paleoceanography* **19**, PA1004 (2004).

31. Álvarez-Solas, J. et al. Heinrich event 1: an example of dynamical ice-sheet reaction to oceanic changes. *Clim. Past* **7**, 1297–1306 (2011).

32. Marcott, S. A. et al. Ice-shelf collapse from subsurface warming as a trigger for Heinrich events. *Proc. Natl Acad. Sci. USA* **108**, 13415–13419 (2011).

33. Alley, R. B., Andrews, J. T., Barber, D. C. & Clark, P. U. Comment on “catastrophic ice shelf breakup as the source of Heinrich event icebergs” by C.L. Hulbe et al. *Palaeoceanography* **20**, PA1009 (2005).

34. Jennings, A. E. et al. Baffin Bay paleoenvironments in the LGM and HS1: resolving the ice-shelf question. *Mar. Geol.* **402**, 5–16 (2018).

35. Hodgson, D. A. & Vincent, J. S. A 10,000 yr BP extensive ice shelf over Viscount Melville Sound, Arctic Canada. *Quat. Res.* **22**, 18–30 (1984).

36. Hodgson, D. A. Episodic ice streams and ice shelves during retreat of the northwesternmost sector of the late Wisconsinan Laurentide Ice Sheet over the central Canadian Arctic Archipelago. *Boreas* **23**, 14–28 (1994).

37. England, J., Bradley, R. S. & Miller, G. H. Former ice shelves in the Canadian High Arctic. *J. Glaciol.* **20**, 393–404 (1978).

38. England, J. H., Furze, M. F. & Doupé, J. P. Revision of the NW Laurentide Ice Sheet: implications for paleoclimate, the northeast extremity of Beringia, and Arctic Ocean sedimentation. *Quat. Sci. Rev.* **28**, 1573–1596 (2009).

39. Furze, M. F., Pienkowski, A. J., McNeely, M. A., Bennett, R. & Cage, A. G. Deglaciation and ice shelf development at the northeast margin of the Laurentide Ice Sheet during the Younger Dryas chronozone. *Boreas* **47**, 271–296 (2018).

40. England, J. H., Coulthard, R. D., Furze, M. F. A. & Dow, C. F. Catastrophic ice shelf collapse along the NW Laurentide Ice Sheet highlights the vulnerability of marine-based ice margins. *Quat. Sci. Rev.* **286**, 107524 (2022).

41. Niessen, F. et al. Repeated Pleistocene glaciation of the East Siberian continental margin. *Nat. Geosci.* **6**, 842–846 (2013).

42. Jakobsson, M. et al. Evidence for an ice shelf covering the central Arctic Ocean during the penultimate glaciation. *Nat. Commun.* **7**, 1–10 (2016).

43. Riedel, M. et al. Mega-scale glacial lineations formed by ice shelf grounding in the Canadian Beaufort Sea during multiple glaciations. *Earth Surf. Process. Landf.* **46**, 1568–1585 (2021).

44. Graham, A. G. et al. Seabed corrugations beneath an Antarctic ice shelf revealed by autonomous underwater vehicle survey: origin and implications for the history of Pine Island Glacier. *J. Geophys. Res.: Earth Surf.* **118**, 1356–1366 (2013).

45. Jenner, K. A., Campbell, D. C. & Piper, D. J. W. Along-slope variations in sediment lithofacies and depositional processes since the Last Glacial Maximum on the northeast Baffin margin, Canada. *Mar. Geol.* **405**, 92–107 (2018).

46. Jakobsson, M. et al. An Arctic Ocean ice shelf during MIS 6 constrained by new geophysical and geological data. *Quat. Sci. Rev.* **29**, 3505–3517 (2010).

47. Løken, O. H. Baffin Island refugia older than 54,000 years. *Science* **153**, 1378–1380 (1966).

48. Davis, P. T., Briner, J. P., Coulthard, R. D., Finkel, R. W. & Miller, G. H. Preservation of Arctic landscapes overridden by cold-based ice sheets. *Quat. Res.* **65**, 156–163 (2006).

49. Refsnider, K. A., Miller, G. H., Fréchette, B. & Rood, D. H. A chronological framework for the Clyde Foreland Formation, Eastern Canadian Arctic, derived from amino acid racemization and cosmogenic radionuclides. *Quat. Geochronol.* **16**, 21–34 (2013).

50. Steig, E. J., Wolfe, A. P. & Miller, G. H. Wisconsinan refugia and the glacial history of eastern Baffin Island, Arctic Canada: coupled evidence from cosmogenic isotopes and lake sediments. *Geology* **26**, 835–838 (1998).

51. Wolfe, A. P., Frechette, B., Richard, P. J., Miller, G. H. & Forman, S. L. Paleocology of a >90,000-year lacustrine sequence from Fog Lake, Baffin Island, Arctic Canada. *Quat. Sci. Rev.* **19**, 1677–1699 (2000).

52. Miller, G. H. et al. The Goldilocks dilemma: big ice, little ice, or “just-right” ice. *Quat. Sci. Rev.* **22**, 33–48 (2002).

53. Briner, J. P., Axford, Y., Forman, S. L., Miller, G. H. & Wolfe, A. P. Multiple generations of interglacial lake sediment preserved beneath the Laurentide Ice Sheet. *Geology* **35**, 887–890 (2007).

54. Axford, Y. et al. Recent changes in a remote Arctic lake are unique within the past 200,000 years. *Proc. Natl Acad. Sci. USA* **106**, 18443–18446 (2009).

55. Miller, G. H. et al. Last interglacial lake sediments preserved beneath Laurentide and Greenland Ice sheets provide insights into Arctic climate amplification and constrain 130 ka of ice-sheet history. *J. Quat. Sci.* **37**, 979–1005 (2022).

56. Aksu, A. E. & Piper, D. J. Late quaternary sedimentation in Baffin Bay. *Can. J. Earth Sci.* **24**, 1833–1846 (1987).

57. Simon, Q., St-Onge, G. & Hillaire-Marcel, C. Late Quaternary chronostratigraphic framework of deep Baffin Bay glaciomarine sediments from high-resolution paleomagnetic data. *Geochem. Geophys. Geosyst.* **13**, Q0AO03 (2012).

58. Simon, Q., Hillaire-Marcel, C., St-Onge, G. & Andrews, J. T. North-eastern Laurentide, western Greenland and southern Innuitian ice stream dynamics during the last glacial cycle. *J. Quat. Sci.* **29**, 14–26 (2014).

59. Simon, Q. et al. Authigenic 10Be/9Be ratios and 10Be-fluxes (230Thxsn-normalized) in central Baffin Bay sediments during the last glacial cycle: Paleoenvironmental implications. *Quat. Sci. Rev.* **140**, 142–162 (2016).

60. Jackson, R. et al. Asynchronous instability of the North American-Arctic and Greenland ice sheets during the last deglaciation. *Quat. Sci. Rev.* **164**, 140–153 (2017).

61. Dowdeswell, J. A. et al. An origin for laminated glacimarine sediments through sea-ice build-up and suppressed iceberg rafting. *Sedimentology* **47**, 557–576 (2000).

62. Domack, E. et al. Stability of the Larsen B ice shelf on the Antarctic Peninsula during the Holocene epoch. *Nature* **436**, 681–685 (2005).

63. Prothro, L. O., Simkins, L. M., Majewski, W. & Anderson, J. B. Glacial retreat patterns and processes determined from integrated sedimentology and geomorphology records. *Mar. Geol.* **395**, 104–119 (2018).

64. Jennings, A. et al. Modern and early Holocene ice shelf sediment facies from Petermann Fjord and northern Nares Strait, northwest Greenland. *Quat. Sci. Rev.* **283**, 107460 (2022).

65. Anderson, J. B., Kurtz, D. D., Domack, E. W. & Balshaw, K. M. Glacial and glacial marine sediments of the Antarctic continental shelf. *J. Geol.* **88**, 399–414 (1980).

66. Tippett, C. R. Glacial dispersal train of Paleozoic erratics, central Baffin Island, NWT, Canada. *Can. J. Earth Sci.* **22**, 1818–1826 (1985).

67. Jackson, G. D. & Berman, R. G. Precambrian metamorphic and tectonic evolution of northern Baffin Island, Nunavut, Canada. *Can. Mineral.* **38**, 399–421 (2000).

68. Andrews, J. T., Klein, A. J., Jenner, K. A., Jennings, A. E. & Campbell, C. The variability of Baffin Bay seafloor sediment mineralogy: the identification of discrete glacial sediment sources and application to Late Quaternary downcore analysis. *Can. J. Earth Sci.* **55**, 620–639 (2018).

69. Andrews, J. T., Jenner, K. A. & Campbell, C. Linking marine core lithofacies and mineral and grain-size compositions on the Baffin Island margin: Changes in provenance and transport. *J. Sediment. Res.* **90**, 763–775 (2021).

70. Ó Cofaigh, C. et al. Glacimarine lithofacies, provenance and depositional processes on a West Greenland trough-mouth fan. *J. Quat. Sci.* **28**, 13–26 (2013).

71. Sheldon, C. et al. Ice stream retreat following the LGM and onset of the west Greenland current in Uummannaq Trough, west Greenland. *Quat. Sci. Rev.* **147**, 27–46 (2016).

72. Riddle, M. J., Craven, M., Goldsworthy, P. M. & Carsey, F. A diverse benthic assemblage 100 km from open water under the Amery Ice Shelf, Antarctica. *Paleoceanography* **22**, PA1204 (2007).

73. Post, A. L. et al. Modern sedimentation, circulation and life beneath the Amery Ice Shelf, East Antarctica. *Cont. Shelf Res.* **74**, 77–87 (2014).

74. Weertman, J. Stability of the junction of an ice sheet and an ice shelf. *J. Glaciol.* **13**, 3–11 (1974).

75. Bassis, J. N. & Walker, C. C. Upper and lower limits on the stability of calving glaciers from the yield strength envelope of ice. *Proc. R. Soc. A. Math. Phys. Eng. Sci.* **468**, 913–931 (2012).

76. Pollard, D., DeConto, R. M. & Alley, R. B. Potential Antarctic Ice Sheet retreat driven by hydrofracturing and ice cliff failure. *Earth Planet. Sci. Lett.* **412**, 112–121 (2015).

77. Tarasov, L., Dyke, A. S., Neal, R. M. & Peltier, W. R. A data-calibrated distribution of deglacial chronologies for the North American ice complex from glaciological modeling. *Earth Planet. Sci. Lett.* **315**, 30–40 (2012).
78. Jakobsson, M. et al. The Holocene retreat dynamics and stability of Petermann Glacier in northwest Greenland. *Nat. Commun.* **9**, 2104 (2018).

79. Gasson, E. G., DeConto, R. M., Pollard, D. & Clark, C. D. Numerical simulations of a kilometre-thick Arctic ice shelf consistent with ice grounding observations. *Nat. Commun.* **9**, 1510 (2018).

80. Lambeck, K., Rouby, H., Purcell, A., Sun, Y. & Cambridge, M. Sea level and global ice volumes from the Last Glacial Maximum to the Holocene. *Proc. Natal Acad. Sci. USA* **111**, 15296–15303 (2014).

81. Kindler, P. et al. Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core. *Clim. Past* **10**, 887–902 (2014).

82. Elvenes, S. & Dowdeswell, J. A. Possible 'lift-off moraines' at Grounded Ice-sheet Margins, *North Norwegian Shelf Edge* 247–248 (Geological Society, 2016).

83. Carlson, A. E. Geochemical constraints on the Laurentide Ice Sheet contribution to meltwater pulse 1A. *Quat. Sci. Rev.* **28**, 1625–1630 (2009).

84. Harrison, S., Smith, D. E. & Glasser, N. F. Late Quaternary meltwater pulses and sea level change. *J. Quat. Sci.* **34**, 1–15 (2018).

85. Dorschel, B. et al. WESTBAFF Reconstruction of the Laurentide Ice Sheet Drainage Into the Northwest Baffin Bay and the Palaoceanography of the West Baffin Bay 52 (MARIA S. MERIAN-Berichte, 2017).

86. Jakobsson, M. et al. The International Bathymetric Chart of the Arctic Ocean (IBCAO) Version 3.0. *Geophys. Res. Lett.* **39**, L12609 (2014).

87. Heaton, T. J. et al. Marine20—the marine radiocarbon age calibration curve (0–55,000 cal BP). *Radiocarbon* **62**, 779–820 (2020).

88. Coulthard, R. D., Furze, M. F., Pieńkowski, A. J., Nixon, F. C. & England, J. H. 2010, New marine ΔR values for Arctic Canada. *Quat. Geochronol.* **5**, 419–434 (2010).

89. Lloyd, J. et al. A 100 yr record of ocean temperature control on the stability of Jakobshavn Isbrae, West Greenland. *Geology* **39**, 867–870 (2011).

## Acknowledgements

We are thankful to the captains, crews and scientific participants on the different cruises for their help and support during the expeditions. This project was funded by the ArcticNet Network of Centers of Excellence, Natural Sciences and Engineering Research Council of Canada (NSERC) and Sentinelle Nord (Apogée Canada) grants to P.L., the Deutsche Forschungsgemeinschaft (DFG) grants to B.D., C.G. and D.H. for ship-time and to D.H. for the International Research Training Group “ArcTrain” (IRTG 1904), as well as Université de Strasbourg grants to P.O.C. We also thank ArcTrain Canada for providing financial support to this project. This paper benefited from fruitful discussions with D. B. Praeg and C. Ó Cofaigh. The final version of this manuscript was greatly improved by thoughtful reviews from M. Margold, J. T. Andrews and two anonymous reviewers.

## Author contributions

P.O.C. conceived the study and initiated this project in cooperation with P.L. and J.F.G. P.O.C. interpreted the geophysical data sets, wrote the paper and prepared the figures. P.O.C., P.L., J.F.G., B.D., C.G., D.H. and E.B. were involved in the discussion of the data and contributed to the final version of this paper.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s43247-022-00559-7.

**Correspondence** and requests for materials should be addressed to Pierre-Olivier Couette.

**Peer review information** *Communications Earth & Environment* thanks John Andrews and Martin Margold for their contribution to the peer review of this work. Primary Handling Editors: Joe Aslin, Heike Langenberg.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022