Mineral and isotopic (Nd, Sr) signature of fine-grained deglacial and Holocene sediments from the Mackenzie Trough, Arctic Canada

ABSTRACT Dolomites in Arctic Ocean sediments are widely attributed to erosion and transport of sediments from northern Canada and Greenland. Coarse-grained dolomite-rich ice-rafted debris is often linked to iceberg transport, but the origin of fine-grained dolomite is less well constrained. A presumed source is the Mackenzie River. In this article, we fingerprint the minerogenic and isotopic (Nd and Sr) composition of deglacial and Holocene fine-grained sediments (<38 μm) from the shallow Mackenzie Trough. Sediments from an 81.5-m borehole (MTW01) were analyzed. The borehole is composed of a progradational (deltaic), transitional (transgressive), and marine unit (<9.4 cal. ka BP). The average dolomite content (~7 percent) and ɛNd signals (−13.3) are surprisingly constant in the progradational and marine units. The isotopic signature is inherited from the Interior Platform, the major underlying bedrock region of the Mackenzie River mainstream. The transitional unit contains fluctuations in ɛNd (−11.0 and −14.6), reflecting enhanced input from the North American Cordillera and Canadian Shield that are not associated with elevated amounts of dolomite. Additional studies combining minerogenic and detrital ɛNd analyses from sites proximal to the paleo-icestreams draining the Canadian Arctic are required to ascertain the origin of dolomite enrichment in central Arctic Ocean sediments.


Introduction
The composition of terrigenous marine sediments in the central Arctic Ocean reflect interactions between riverine inputs and terrestrial runoff, ocean circulation, sea ice drift, and the dynamics of circum-Arctic ice sheets, which extended across large parts of the continental shelves during past glacial periods (Stein 2008;Jakobsson et al. 2014;Wegner et al. 2015). On many topographic highs within the Arctic Ocean, terrigenous sediments are largely composed of sea ice and icebergrafted material (Stein 2008), although the contribution from suspended sediments carried in different water masses remains uncertain. The mineral and isotopic (i.e., Nd, Sr, Pb) composition of these sediments can, in some instances, be used to identify their provenance, or source region (e.g., Fagel et al. 2014;Maccali, Hillaire-Marcel, and Not 2018). This provides a powerful tool for paleoceanographic reconstructions. For example, provenance studies of arctic sediments have been used to reconstruct ice drift patterns during glacial periods (Bischof, Clark, and Vincent 1996;Phillips and Grantz 2001), interpret the configuration and calving dynamics of continental ice sheets (Darby et al. 1989(Darby et al. , 2002Spielhagen et al. 2004;Stokes et al. 2005;Polyak et al. 2009;Kaparulina, Strand, and Lunkha 2016;Bazhenova, Fagel, and Stein 2017;Dong et al. 2017), identify the inception of perennial sea ice in the Arctic (Krylov et al. 2008;Darby 2014), and investigate millennial-scale dynamics in oceanic gateway exchanges (Darby et al. 2012).
The clay and heavy mineral composition of surface sediments from some arctic continental shelves (Barents, Kara, Laptev, East Siberian and Alaskan) have been fairly well documented and tied to the geology of the surrounding landmasses (Stein [2008] and references therein). Although no single mineral is uniquely associated to a local region (Stein 2008;Darby et al. 2011), some minerals and mineral assemblages can be convincingly tied to different areas of the Arctic. For example, the clay mineral smectite is most abundant in the Kara Sea region of the Arctic, due to weathering of the Putorana flood basalts (Vogt and Knies 2008). Dolomite in arctic sediments is widely associated with weathering of outcropping Paleozoic dolomites and limestones that extend across the Canadian Arctic Archipelago and exist as far east as Northern Greenland (Krawitz 1982;Harrison et al. 2011; Figure 1), although some ocean bottom outcrops in the Amerasian Basin have also been described (Skolotnev, Fedonkin, and Korniychuk 2017). The input of fine and coarse fraction dolomites from this broad area has resulted in Arctic Ocean sediments being considered some of the most dolomite-rich marine sediments worldwide (Darby et al. 1989;Stein 2008). This is exemplified by the often described and recurrent "pink layers" found in Pleistocene sediment cores from the Amerasian Basin (Clark et al. 1980;Stein et al. 2010;Bazhenova, Fagel, and Stein 2017), which can contain more than 50 wt% of dolomite in the <63 μm size fraction (Dong et al. 2017).
Historically, the dolomite content of central Arctic Ocean sediments was more commonly investigated in the coarser fractions (>63 µm) and used as a proxy for iceberg delivered material calved from the Laurentide and Innuitian ice sheets (Clark et al. 1980;Bischof, Clark, and Vincent 1996;Phillips and Grantz 2001;Spielhagen et al. 2004). Although elevated dolomite abundances in Pleistocene sediments from the Amerasian Basin are generally accompanied by an overall increase in the >63 µm fraction, the dominant proportion of dolomite is still found in the <63 μm fraction (Carrol 1970;Dalrymple and Maas 1987;Darby et al. 1989;Vogt 1997).
To better discriminate between the sources of dolomite in Arctic Ocean sediments, measurements of radiogenic isotopes (Sr, Nd, and Pb) have been investigated in a number of recent studies (Not and Hillaire-Marcel 2012;Hillaire-Marcel et al. 2013;Fagel et al. 2014;Bazhenova, Fagel, and Stein 2017;Maccali, Hillaire-Marcel, and Not 2018;Dong et al. 2020;Wang et al. 2021;Xiao et al. 2021;Schwab et al. 2021). One focus of these studies is on the mineral and isotopic composition of sediments from the Canadian Beaufort Sea shelf and the wider Canadian Arctic Archipelago. The ability to identify a unique Mackenzie River provenance is of particular importance, because the 4,200-km-long Mackenzie River drains an area of 1.8 × 10 6 km 2 ( Figure 1) and delivers the largest annual load of suspended material to the Arctic Ocean of all Arctic rivers (124 × 10 12 g/yr; Carson, Jasper, and Conly 1998;  (Harrison et al. 2011). Yellow area is the Mackenzie River Delta. SL = Great Slave Lake. Base map from International Bathymetric Chart of the Arctic Ocean (IBCAO) (Jakobsson et al. 2008). Holmes et al. 2002). The estimated annual contribution of sediment from the Mackenzie River to the Beaufort Sea (after subtracting estimated delta sedimentation) is 65 × 10 12 g/yr (Macdonald et al. 1998;Rachold et al. 2000). Furthermore, the Mackenzie River has long been implicated as an important route for deglacial meltwater outbursts from the retreating Laurentide Ice Sheet (Fisher, Smith, and Andrews 2002;Tarasov and Peltier 2005;Peltier, Vettoretti, and Stastna 2006;Murton et al. 2010;Condron and Windsor 2012;Keigwin et al. 2018;Klotsko, Driscoll, and Keigwin 2019;Süfke et al. 2022). In this study, mineral and isotopic (Sr and Nd) compositions of fine fraction (<38 μm) sediments are reported from an 81.5-m deglacial to Holocene sediment sequence (MTW01) recovered in the Mackenzie Trough. The aim is to use a source to sink perspective and better constrain the mineralogical and isotopic signature of Mackenzie River sediments and investigate to what degree these changed during deglaciation and the Holocene.

Geological context of the Mackenzie River Basin
The underlying bedrock of the Mackenzie River basin can be divided into three major domains: (1) the folded Proterozoic to Mesozoic sedimentary strata of the North American Cordilleras (NAC), (2) the undeformed Cambrian to Cretaceous sedimentary strata of the Interior Platform (IP), and (3) the Precambrian igneous/metamorphic rocks of the Canadian Shield (CS; Figure 1). Prominent carbonate bedrock areas, including dolostones and limestones, are widespread in the Mackenzie River basin and are more pronounced in the NAC and IP (Figure 1). The Mackenzie River basin was covered by the Laurentide ice sheet during the last glacial period but became gradually ice-free between 18 and 10 ka BP as the Laurentide ice sheet retreated eastwards (Dyke 2004;Stokes and Tarasov 2010;Dalton et al. 2020). Substantial glacial deposits are common all over the Mackenzie basin and constitute an additional source of prospective Arctic Ocean sediments (Millot et al. 2003).
Though the major tributaries of the Mackenzie River cross large areas of outcropping dolomitic bedrock (Figure 1), there exists little data on the actual dolomite content of fine fraction material carried/deposited by the Mackenzie River and how this may have changed through time. The majority of fine fraction mineral studies incorporating sediments from the Mackenzie River or Beaufort Sea either removed the carbonate bearing minerals during pretreatments or did not analyze for or report dolomite contents (Table 1). Measurements of modern suspended sediments from the Mackenzie River indicate that detrital carbonates (<500 µm) constitute 5 wt% of the total sediment delivery during fall and winter (1987/1988) and 2.5 wt% during spring and summer (1988) (O'Brien et al. 2006). Darby et al. (2011) analyzed the dolomite content in the <45 μm fraction of two samples from the Beaufort Sea using X-ray diffraction (XRD), one sample came from the Amundsen Gulf (8.1 percent dolomite), and one sample from the Mackenzie River delta (14.6 percent dolomite; Figure 2 and Table 2). In a more recent study, the bulk mineralogy of thirty-four surface sediments from across the Canadian Beaufort Shelf reveal a large range in dolomite content from 2.7 to 34.8 wt%, with a clear gradient toward higher concentrations moving east from the Mackenzie Trough toward the Amundsen Gulf ( Figure 2; Gamboa et al. 2017).
Radiogenic Sr in the dissolved load of Mackenzie River waters and its tributaries is directly related to the age of the bedrock source. Accordingly, low 87 Sr/ 86 Sr (0.7074-0.7140) in the river waters draining the IP contrasts with high 87 Sr/ 86 Sr (0.7346-0.7517) in the tributaries draining the CS (Figure 3; Millot et al. 2002;Millot et al. 2003). The corresponding 87 Sr/ 86 Sr for the NAC fraction of Mackenzie River tributaries is 0.7066-0.7580, reflecting the heterogeneous age composition of this mountain range (Figure 3; Millot et al. 2003). Surprisingly stable 87 Sr/ 86 Sr ratios (0.7114-0.7115) are recorded in the Mackenzie mainstream along a 2,000km transect extending from the Slave River south of Great Slave Lake to Tsiigehtchic located slightly upstream from the confluence of the Peel and Mackenzie Rivers (Figure 1b; Millot et al. 2003). This stable and low radiogenic signal is explained by low contributions from the CS to the Mackenzie River mainstream due to less intense chemical weathering of the Precambrian bedrock (Millot et al. 2003).
The suspended bulk load of the IP tributaries shows 87 Sr/ 86 Sr that are consistently higher than the dissolved load (between 0.7151 and 0.7204), and the NAC tributaries 87 Sr/ 86 Sr values are between 0.7129 and 0.7368 ( Figure 3; Millot et al. 2003). The 87 Sr/ 86 Sr of the suspended load changes from 0.7199 at Liard River confluence with the Mackenzie River mainstream to 0.7189 at Tsiigehtchic (Millot et al. 2003). This observation is explained by estimates that only 5 percent of the <0.125 mm suspended material derives from the more radiogenic CS tributaries (Carson, Jasper, and Conly 1998). Carson, Jasper, and Conly (1998) pointed out that the majority (65 percent) of the <0.125 mm suspended load is eroded from the mainstream riverbanks of the Mackenzie and Liard rivers; hence, the Interior Platform is the major provenance source of the suspended load.
Nd isotope measurements from the Mackenzie Basin (bedrock, suspended load, or dissolved load) are less abundant than Sr data. However, ɛ Nd of lake sediment from the Mackenzie River delta exhibit a west-to-east gradient ( Figure 3; Vonk et al. 2015). More radiogenic signals (ε Nd from −11.2 to −12.0) in the western parts of the delta are attributed to inputs from the NAC via the Peel River while the eastern parts are fed by the Mackenzie River mainstream (IP) and exhibit lower Nd isotope values (ε Nd from −12.5 to −14.5; Figure 1).

Site description and borehole stratigraphy
The 150-km-long Mackenzie Trough crosses the Canadian Beaufort shelf in a south-southeast direction from the modern Mackenzie River delta (Figure 1b). It is believed to be a glacially excavated trough (Shearer 1971;Blasco et al. 1990) occupied by an ice stream draining the Laurentian ice sheet during the Illinoian/Early Wisconsinan and Late Wisconsinan glaciations (Batchelor, Dowdeswell, and Pietras 2013). Consequently, seismic reflection data reveal two major subglacial till units in the Mackenzie Trough constituting up to 500 m of unconsolidated material upon the bedrock surface (Blasco et al. 1990;Batchelor, Dowdeswell, and Pietras 2013). The topmost seismic unit (MT1: Blasco et al. 1990; Sequence 5: Batchelor, Dowdeswell, and Pietras 2013) consists of 30 to 200 m of deglacial and Holocene deposits, with the thickest deposits on the landward side of the trough (Blasco et al. 1990;Hill et al. 1991;Batchelor, Dowdeswell, and Pietras 2013).
The 81.5-m borehole (MTW01, 69°20′53″N, 137°59′ 13″W) was drilled at 45-m water depth in the Mackenzie Trough in 1984 by the Geological Survey of Canada using an anchored drilling barge ( Figure 1b). Sections of the borehole (0.3-0.8 m) were recovered every 1 to 3 m and subsequently analyzed with respect to sedimentology and physical properties (Moran, Hill, and Blasco 1989). Together with data from a 71-m piezocone penetrometer profile from a nearby (50-m offset) site and a seismic reflection profile crossing the borehole, five lithological units were defined (Figure 4). Units C-E (29-81 mbsf) consist of sandy silts and are interpreted as a progradational sequence deposited in a deltaic setting when sea level was lowered during the last glacial cycle. Unit B (21-29 mbsf) is deposited unconformably onto Unit C during a transgressive phase and consists of laminated silt and clay, which show signs of deformation near the top of the unit. Near the top of Unit B (21-23 mbsf), planktonic foraminifera and ostracods first appear and indicate the onset of a marine environment (O'Regan  Cook et al. (1975) Average <1% Higher concentration in sandy layers Dalrymple and Maas (1987) Canada Basin off Chukchi Plateau and Alpha Ridge Bulk sediment 5 sediment cores (15-202 mbsf; Plio-Pleistocene to modern sediment) Gasometric analysis + XRD Quantification according to Griffin (1971) Average 12.3% and standard deviation 9.6 Darby et al. (1989) Canada . Unit A lays unconformably on Unit B and is a homogeneous foraminifer-bearing fine-grained clay deposited after marine inundation at the site (Moran, Hill, and Blasco 1989). Radiocarbon measurements of shell fragments from mollusks and mixed benthic foraminifer from Units A and B provide an estimate on the age of deposition of 9180 � 310 350 cal. yr BP for the base of Unit A (20.73-21.09 mbsf) and give a linear sedimentation rate of 2.63 ± 0.05 m ka −1 for Unit A ( Figure 4; Hilton et al. 2015;O'Regan et al. 2018). Radiocarbon age calibrations on foraminifera and mollusk fragments were performed with Marine13 (Reimer et al. 2013) and a ΔR of 335 ± 85 years based on twenty-four living mollusk specimens collected before 1956 from the northwestern Canadian Arctic Archipelago (Coulthard et al. 2010). Measurements on particulate organic matter were calibrated using IntCal13 data set (Reimer et al. 2013). The topmost sample in Unit B (21.64-21.99 mbsf) provides an age of 9920 9180 � 370 290 cal yr BP and constitutes the oldest absolute date in the MWT01 sequence. Although no absolute dates constrain the age of the progradational sequence, it was initially reported as being younger than ~14 ka yr BP based upon unpublished pollen assemblage data (Blasco et al. 1990). Radiocarbon measurements of particulate organic carbon from the >63 μm, which are the youngest form of carbon delivered to the site by the Mackenzie River (Hilton et al. 2015), supports this general assertion indicating that it is younger than 24820 þ390 À 380 cal. yr BP and was most likely deposited after glacial ice retreated from the coast at the end of the Toker Point glaciation (sometime between 22 and 16 cal. ka BP; Murton et al. 2007;O'Regan et al. 2018). This inferred age is largely consistent with the most recent ice margin chronologies, which show pronounced inland retreat of the Mackenzie ice stream between 14.9 and 14.2 cal. ka BP (Dalton et al. 2020). The implication is that the Younger Dryas would have occurred during deposition of sediments in MTW01. However, the lack of absolute dates below 21.64-21.99 mbsf precludes identification of its exact location. Murton et al. (2010) suggested that it might have occurred across the Unit C/ B boundary, just prior to marine inundation of the site.

Mineral assemblage investigations
The twenty-three sieved sediment samples (<38 μm) from MTW01 were homogenized using a synthetic sapphire mortar before being measured in an X'PERT X-ray diffractometer at the Department of Geological Sciences, Stockholm University. The diffractometer was equipped with Cu anode (Kα = 1.541) and a fixed divergence slit (1°/1.52 mm) and data were collected at 3° to 85° 2θ with a step size of 0.016° measuring 0.5 seconds per step. Postprocessing in HighScore 3.0x (PANalytical B.V. Almelo, The Netherlands) included baseline determination and peak search using default settings and subsequent peak identification using the ICSD 3.2 mineral database (Supplementary material 1). The mineral assemblage for each sample was quantified by Rietveld refinement using the free Material Analysis Using Diffraction (MAUD) software (http://maud. In order to produce comparable data sets, the Rietveld refinement process was fully automatized including five iterations for each sample. The statistical significance of each refinement is given by the expected R value (R exp = theoretically best possible fit) and weighted profile R value (R wp = the real fit after refinement). Relatively high R exp (22.3-23.1 percent) values were returned, reflecting the heterogeneous character of the samples (Toby 2006). However, the actual R wp values are 37.3 to 41.1 percent, constituting less than twice (R wp � 1:75 � R exp ) the corresponding R exp . This indicates a decent fit of the refinement curve (McCusker et al. 1999).

Sr and Nd isotopic analysis
Ten of the fine-grained (<38 µm) samples were prepared for 87 Sr/ 86 Sr and 143 Nd/ 144 Nd analysis. Measuring the isotope composition of marine sediments provides a combined isotope signal partly derived from the detrital component and partly from a secondary component. The latter typically has an isotope composition derived from the water mass from which they were precipitated and usually appear as Fe-Mn hydroxide coatings (e.g.,   1988;Frank 2002). These can be extracted by leaching procedures to obtain a pure detrital and/or seawater isotopic signal (e.g., Bayon et al. 2002;Gutjahr et al. 2007). However, detrital carbonates (e.g., dolomite and calcite) will dissolve using the traditional methods of leaching and can potentially remove an important component of the isotopic signals detrital signature (Hillaire-Marcel et al. 2013). In order to gain a representative detrital signal that incorporates the carbonate phases, we analyzed Nd and Sr isotopic compositions of the bulk (<38 μm) sediment samples, bearing in mind that part of the isotopic signal likely has a water mass origin.
For the ten selected samples, 0.1 g of dried sediment, a BCR-2 rock standard, and a procedural blank were pretreated with 0.5 mL 8 N HNO 3 and approximately 5 drops of 30 percent H 2 O 2 to oxidize organic material. The samples were then completely digested in a successive treatment using 2 mL concentrated hydrofuran and 20 drops of concentrated HNO 3 at 205°C for 72 hours, 5 × 10 drops concentrated HNO 3 with intermediate evaporation, and finally 5 mL 6 N hydrochloric acid (HCl) at 205°C for 24 hours. Element separation of Nd and Sr was performed according to Aciego et al. (2009). Iron was eliminated using AG1-X8 200 to 400 mesh (Bio-Rad) and TRU resin B 50 to 155 µm (TrisKem) in an HCl environment. Rare Earth Elements (REEs) were separated from Sr and other cations using an AG1-X8 TRU resin B 50 to 150 µm in a HNO 3 environment. Subsequently, Sr was isolated from the matrix using a Sr resin B 50 to 100 µm (TrisKem). The Nd fraction of the REE was isolated using a LN resin C 100 to 150 µm (TrisKem) in an HCl environment. All chemicals were of high-purity quality (Seastar Baseline).
Strontium and Nd isotopes were analyzed at the Swedish Museum of Natural History on a Thermo Scientific TRITON thermal ionization mass spectrometer using total spiking with a mixed 149 Sm/ 150 Nd spike. Analytical details for the Sr and Nd isotope determinations are reported in Table 3. Total procedural blank for Sr was recorded to 1.9 ng. The 87 Sr/ 86 Sr for the BCR-2 standard had a value of 0.705027 ± 0.000015, within the range of the reported normalized value (0.705010 ± 0.000010; Weis et al. 2006).
Total procedural blank for Nd was recorded to be 0.23 ng. The 143 Nd/ 144 Nd for the BCR-2 standard had a value of 0.5126511 ± 0.0000014, close to the reported normalized value (0.512627 ± 0.000012; Weis et al. 2006).

Grain size measurements
All twenty-three samples (<38 µm) used for mineral analysis were measured with respect to grain size distribution using a Malvern 3000 laser diffractometer. Grain size analysis was performed to reveal the distribution between clay-and silt-sized particles in the samples in order to support interpretation of changes in 87 Sr/ 86 Sr ratios. Calculation of mean grain size for the samples was done using the geometric method of moments (Blott and Pye 2001).
Representative minerals were chosen based on the most high-ranked suggestions from HighScore processing. However, it is most likely that various minerals contribute to these mineral groups and therefore the output data are presented as feldspars, illite, and chlorite/kaolinite. By default, the identified mineral assemblages sum to 100 percent. The downhole quantitative mineralogical variation exhibits a significant drop in quartz and calcite content at the Unit B/Unit A boundary along with a corresponding increase in illite ( Figure 6). Feldspar, chlorite/kaolinite, and dolomite fractions are relatively constant throughout the core (Figure 6). Dolomite content varies within 4 to 11 wt% (Figure 7) with an average of 7.6 ± 0.6 wt% throughout MTW01. The average dolomite content is 6.5 ± 0.6 wt% and 8.4 ± 0.6 wt% for Units A-B and Units C-E, respectively.

Sr and Nd isotopes and concentrations
The 87 Sr/ 86 Sr values vary from 0.716068 to 0.719041 in Units C-E but become significantly more radiogenic (0.7244435 to 0.726687) in the top sample from Unit B and throughout Unit A (Table 3 and Figure 7). Concentrations of Sr vary from 180 to 210 ppm in the progradational unit and from 150 to 170 ppm in the transitional and marine units. The Nd concentration varies from 29 to 39 ppm with a higher concentration (~36-38 ppm) in Unit A (Table 3 and Figure 7). The ɛ Nd exhibit a stable downhole pattern with values between −13.6 to −12.8 in both Units C-E and Unit A, whereas notable ɛ Nd excursions are found in Unit B (−11.0 and −14.6; Table 3 and Figure 7).

Grain size distribution
The >38 µm content of samples from MTW01 was first presented by O' Regan et al. (2018) and was based on sample weights after sieving. This data illustrated the expected coarser nature of the progradational units, followed by a fining in the transitional and Holocene marine sediments (Figure 7). Analysis of the clay and silt composition of the <38 µm fraction reveals that the clay fraction is 10 to 20 percent in Units B-E and does not increase significantly until the top of Unit B where it then accounts for 30 to 40 percent (Figure 7 and Table 4). This implies that fine fraction sediments in the transitional unit (Unit B) are much more silt rich than the overlying Holocene marine sediments in Unit A (Figure 7).

Mineral assemblage and dolomite content
The XRD-derived mineral assemblages from MTW01 are consistent with analyses of surface sediments from the modern Mackenzie River delta where quartz, calcite, dolomite, feldspars, chlorite, and illite (order of abundance) are reported (Wagemann, Brunskill, and Graham 1977). Suspended material sampled from the river delta had a similar composition but with chlorite and illite being more abundant than feldspars (Wagemann, Brunskill, and Graham 1977). A quantitative study of <45 µm sediment from surface sediments of the Beaufort Shelf also indicate 46.2 percent quartz, 14.6 percent dolomite, 13.1 percent illite, 11.2 percent feldspars, 5.5 percent chlorite, 4.7 percent calcite, 2.7 percent smectite, and 2 percent minor accessory minerals (Darby et al. 2011). In the Mackenzie Trough, bulk mineral assemblages from a surficial sediment sample included (in descending abundance) quartz, illite, chlorite + kaolinite, dolomite, feldspar, and calcite (Kobayashi et al. 2016). Gamboa et al. (2017) reported similar abundance patterns in three bulk surface sediments from the Mackenzie Trough dominated by (in descending order) phyllosilicates, quartz, plagioclase, K-feldspar, and dolomite, with other minerals contributing <1 wt%. General consistency between these different studies suggests that this assemblage is characteristic for the Mackenzie River, which is estimated to supply >95 percent of the sediments found on the Beaufort Shelf (Hill et al. 1991).
One important finding from this study is that there is a remarkably consistent dolomite content in the sediments from MTW01. During the Holocene (Unit A), the dolomite content has remained stable at 6.5 ± 0.6 wt% of the <38 µm fraction, with slightly higher abundances during the period of progradation (8.4 ± 0.6 wt%). Compared to other proximal records, our average Holocene dolomite content is slightly higher than the three bulk surface sediments from the Mackenzie Trough reported by Gamboa et al. (2017), which ranged between 3.6 and 4.4 wt%, but this difference could be ascribed to different size classes analyzed. However, our data and those from Gamboa et al. (2017) are much lower than the 14.6 wt% dolomite reported from one station close to the Mackenzie Trough by Darby et al. (2011), which appears as an outlier when compared to regional measurements (Figure 2).
Two downcore records of XRD-based mineral assemblages exist from the slope of the Canadian Beaufort Sea, a late Holocene record (AMD0214-03PC/TWC; Kutos, Rochon, and Montero-Serrano 2021), and a late deglacial/Holocene record (<13 ka; AMD0214-02PC; Figure 2; ). In the Late Holocene, the dolomite content in 3PC is <10 wt% (Kutos, Rochon, and Montero-Serrano 2021). Kutos, Rochon, and Montero-Serrano (2021) argued that this is a signature of fine-grained material from the Mackenzie River, an argument that agrees well with the dolomite content of sediments from MTW01. In AMD0214-02PC, located further to the east, sediments younger than 10 ka have an average dolomite content of 1.4 ± 2.4 wt% (n = 18), whereas an underlying laminated and ice-rafted debris (IRD) enriched unit deposited between 13 and 10 ka has a mean of 6.8 ± 7.2 wt% (n = 17). Two discrete intervals with peak dolomite contents of 28 and 15 wt% occur in this laminated unit and are dated to ~10.8 and 12.9 cal. ka BP, respectively . This Ca-enriched laminated unit is mapped across the Beaufort slope and appears thickest in front of the Mackenzie Trough (Klotsko, Driscoll, and Keigwin 2019), implying that the Mackenzie River is the primary source of this material. Klotsko, Driscoll, and Keigwin (2019) also suggested that discrete Ca-enriched IRD intervals in this regionally identified laminated unit are likely derived from other sources. Specifically, these include inputs from ice stream activity and meltwater emanating from the large Amundsen and possibly Mc'Clure straits to the east. This is consistent with the considerably higher dolomite contents in surface sediments from the Amundsen Gulf and presumably also further east (Figure 2). Pronounced glacial retreat occurred in Amundsen and Mc'Clure straits during the Bølling-Allerød and early Younger Dryas (Lakeman et al. 2018), with the Amundsen Gulf being implicated  as key sources of meltwater discharge during the Bølling-Allerød (~14.7-13 ka; Süfke et al. 2022).
In bulk sediment (i.e., nonsieved sediment samples), dolomite quantities from the Amerasian Basin have been reported by Darby et al. (1989), Vogt (1997), Stein et al. (2010), Bazhenova (2012), and Zao (2016). Darby et al. (1989) documented dolomite contents of up to 46 percent in five Plio-Pleistocene sediment cores from ridges and abyssal plains in northwestern Amerasian Basin (Figure 1). Stein et al. (2010) revisited one of these sites and found a similar downcore pattern of dolomite abundance up to 50 percent in the so-called pink layers. Vogt (1997) reported a range of 2 to 20 percent for dolomite abundances in five surficial sediment samples from the Canada Basin with maximum abundances on the Chukchi Plateau (Figure 1). At the Alaska/Chukchi Sea border, high levels (<35 percent) of dolomite associated with pulses of IRD were found in pre-Last Glacial Maximum (LGM) sediment (<250 µm) and were followed by low and stable (<2 percent) dolomite levels during the LGM ; Figure 1). The same core exhibited variations between 0 and 33 percent of dolomite during mid-to late Holocene with an average dolomite content of 4 ± 7 percent ; Figure 1).
The use of clay mineralogy in ocean sediment provenance research is widespread since the late 1960s (e.g., Biscaye 1965;Carroll 1970;Naidu and Mowatt 1983;Viscosy-Shirley et al. 2003), and a number of studies have included samples from the Beaufort Shelf (Naidu, Burrell, and Hood 1971;Clark et al. 1980;Naidu and Mowatt 1983), as well as in adjacent areas in the Amerasian Basin (Carrol 1970;Darby 1975;Naidu, Mowatt, and Creager 1982;Moser and Hein 1984). However, the dolomite content in the clay-sized fraction is seldom reported, sometimes due to dissolution of carbonates in the pretreatment procedures (e.g., Naidu, Burrell, and Hood 1971;Moser and Hein 1984; Table 1). Nevertheless, Carrol (1970) reported an average dolomite content of 5 percent in the clay fraction (<2 µm) from seventeen sediment cores recovered in central Canada Basin and on the Alpha Ridge (Table 2). Dalrymple and Maas (1987) reported an average of 1 percent dolomite in the clay fraction of two sediment cores from the Alpha Ridge covering the Plio-Pleistocene to Holocene period (Table 2). Vogt (1997) reported less than 1 wt% of dolomite content for Chukchi and East Siberian Sea shelf sediments. The <63 μm fraction of sediments from a Mendeleev Ridge sediment core (HLY0503-11MC8) contained a wide range (1-28 percent) of dolomite deposited during the last 250 ka (Table 2; Not and Hillaire-Marcel 2010). East of the Mendeleev Ridge the same fraction is reported to contain 0 to 53 percent during the Pleistocene with peaks in dolomite found in discrete pink layers (Table 2; Dong et al. 2017). Fagel et al. (2014) reported between 0 and 10 percent dolomite in the 20 to 63 μm fraction of a sediment cores from the Mendeleyev Ridge. Analysis of separate grain size fractions (<20 µm, 20-63 µm, and >63 µm) from the HLY0503-18MC sediment core recovered in the intra-basin of the Lomonosov Ridge revealed dolomite contents of 0.5 to 6 percent, 2.5 to 15 percent, and 0 to 25 percent, respectively (Table 2; Not and Hillaire-Marcel 2012).
The relatively constant dolomite content in the <38 μm fraction (6-8 wt%) of the MTW01 borehole agrees relatively well with the low-end estimates of dolomites in the fine fraction studies (<63 μm) from the Amerasian Basin (Table 2). In contrast, the elevated (>15-20 wt%) dolomite levels reported in many of these studies, and especially those from the dolomite-rich pink layers (Table 2), far exceed the levels in deglacial and Holocene sediments from MTW01 and other proximal records influenced by Mackenzie River runoff.
Though clearly the Mackenzie River is a source of fine fraction dolomite to the Arctic, at the resolution of this study, there is little evidence for intervals of pronounced dolomite delivery during the Holocene and the period of deglaciation captured in the progradational sequence. It remains possible that more dolomite-enriched fine fraction material delivered by the Mackenzie River bypassed MTW01 during the progradational phase. However, the depositional age of this period in the borehole remains poorly constrained. It is likely younger than 22 to 16 ka, being deposited after glacial ice retreated from the region (Murton et al. 2007;O'Regan et al. 2018), and was originally suggested to be younger than 14 ka based on unpublished pollen data (Blasco et al. 1990). Future work is needed to quantify the dolomite content in older and more slowly deposited laminated sediments found toward the seaward end of the Mackenzie Trough, as well as the amount of glacially eroded detrital carbonate delivered by icestream and meltwater activity in the larger glacial troughs to the east (Marigold et al. 2015). Historically these have been inferred to be the primary source of coarse-grained (IRD) dolomite found in central Arctic Ocean sediments (Stokes et al. 2005), and given the close proximity to large areas of outcropping dolomite (Figure 1), they are a likely path for the delivery of fine-grained dolomites entrained in meltwater (i.e., Süfke et al. 2022).

Nd isotope signature
The Nd isotope signature of MTW01 is relatively constant, with the exception of two measurements from the transitional Unit B (Figures 7 and 8). If the two extreme data points from Unit B (ɛ Nd = −11.0 and −14.6) are excluded, ɛ Nd values in MTW01 span between −12.8 and −13.6 with an average of −13.3 ± 0.1 for Units C-E and Unit A (Figure 8). Neodymium isotope composition in MTW01 (Table 3) is in close agreement with Mackenzie River suspended sediments (bulk: ɛ Nd = −14.3; Goldstein et al. 1984; and detrital: ɛ Nd = -11.6 to −13.9; Vonk et al. 2015). Modern detrital Nd surface samples of the Beaufort Shelf are also broadly consistent with the MTW01 isotope composition (Figure 8; Asahara et al. 2012;Vonk et al. 2015;Maccelli et al. 2018;Schwab et al. 2021). However, our ɛ Nd results are less negative than much of the published surface sediment data from the region (Figure 9). The compiled surface sediment analyses were performed after removing the carbonate and authigenic (Fe-Mn) fractions, which both have a less negative ɛ Nd compared to the detrital silicate phase (Asahara et al. 2012;Schwab et al. 2021). Authigenic Nd phases would primarily be derived from the waters of the Mackenzie River (ɛ Nd = −12.9 ± 0.3; Zimmerman et al. 2009), with a possible contribution from surface waters in the Canada Basin (ɛ Nd = −6.4 to −8.7; Porcelli et al. 2009). Thus, the reported values from MTW01 overestimate the detrital signal and thus the actual detrital ɛ Nd should be lower.
The relatively invariant ɛ Nd values in Units A and C-E implies a joint and fairly stable sediment source. A similarity to the Nd isotope composition of the eastern delta lake sediment (−12.5 to −14.5) implies a substantial IP imprint on these units ( Figure 8; Vonk et al. 2015). The relatively low ɛ Nd signals (−13.7 and −14.7) in two detrital surface sediment samples at the river mouth (Asahara et al. 2012) are likely also affected by the rivers of Mackenzie mainstream and its IP imprint. This is consistent with a minor contribution to the Mackenzie River mainstream from the east bank tributaries (CS), which are estimated to only account for ~5 percent of the total suspended load today (Carson, Jasper, and Conly 1998). The results from this study suggest that the IP remained the dominant source of fine-grained suspended sediments in the Mackenzie River during deglacial progradation and the Holocene. Reported Nd concentration of both modern river sediments (16.09 ppm; Goldstein et al. 1984) and Beaufort Shelf surficial detrital sediments (27.2-27.6 ppm; Asahara et al. 2012) are lower than the concentrations of Nd in MTW01. However, an equivalent increase in Nd concentration from the deltaic sediments in Units C-E to the marine setting in Unit A is seen in MTW01 (Figure 7). This trend might reflect authigenic uptake of Nd from seawater.

Sr isotope signature
The average 87 Sr/ 86 Sr ratio from Units C-E in this study (0.717206; Table 3) is slightly lower than the modern riverine bulk sediment load at Tsiigehtchic (0.7190; Millot et al. 2003). Modern Beaufort Shelf detrital sediment reveal a 87 Sr/ 86 Sr signature considerably higher (Figures 8 and 9) than Units C-E but is closer to the average 87 Sr/ 86 Sr signature (0.724612) in Units A and B (Figures 7 and 8). Our bulk sediment data will be affected by both carbonate and authigenic phases that have a lower 87 Sr/ 86 Sr ratio than detrital silicate fractions (Asahara et al. 2012;Schwab et al. 2021). For example, sediments from MTW01 in Unit A would likely be affected by an authigenic Sr signal that usually resembles the common 87 Sr/ 86 Sr ocean water signal of 0.70918 (Henderson et al. 1994;Gutjahr et al. 2007). Similar to ɛ Nd , the authigenic and carbonate phases in our samples introduce a substantial offset compared to the pure detrital silicate signal reported for surface sediments (Figures 8 and 9). Overall, the progradational units (C-E) tend to have a signature similar to modern Mackenzie River bulk sediment, whereas the radiogenic Sr signal in the marine unit (A) shifts toward a shelf sediment signature. Furthermore, the strontium concentrations in Units C-E (189-211 ppm) (Figure 7) are in good agreement with reported concentrations for suspended bulk sediment in the main tributaries of the Mackenzie River (169-204 ppm; Millot et al. 2003) but are considerably higher than Beaufort Shelf surficial sediment Sr concentrations (113-118 ppm).
Changes in 87 Sr/ 86 Sr can arise from variations in source area, grain size, or physical erosion (Tütken et al. 2002;Meyer, Davies, and Stuut 2011;Ali et al. 2015). The ɛ Nd data suggest a stable sediment provenance through most of the MTW01 sequence, implying that the change in 87 Sr/ 86 Sr between the progradational and marine unit may not reflect a provenance change. The more radiogenic 87 Sr/ 86 Sr ratios in Unit A are accompanied by a reduction in the Sr concentration, an increase in the fine fraction content, and a clay mineral shift toward increased illite and decreased chlorite concentration (Figures 6 and 7). High relative abundance of 87 Sr in clays (i.e., illite) is common as the high Rb concentration produces 87 Sr through radioactive decay (Meyer, Davies, and Stuut 2011;Ali et al. 2015). The 87 Sr/ 86 Sr ratio in MTW01 is strongly correlated to the overall increase in clay content of the samples ( Figure 10). Therefore, the higher 87 Sr/ 86 Sr ratio in the marine unit and the lower overall Sr concentration likely reflect a relative increase in the clay mineral concentration and composition, compared to the progradational sediments ( Figure 10). Hence, 87 Sr/ 86 Sr is not interpreted to reflect a change in sediment provenance, such as a larger contribution from the CS or NAC tributaries in the Holocene. Instead, we suggest that the increase in 87 Sr/ 86 Sr in Unit A is a consequence of the fining of sediment, suggesting that coarser material has been trapped in the modern river delta and nearshore area following transgression. Notably, calcite tends to be the major Sr-bearing mineral in Mackenzie River sediments (Millot et al. 2003). The general decrease in Sr concentration in Unit A might be explained by the observed decrease in calcite.

Mineral and isotopic signatures of meltwater events
Unlike 87 Sr/ 86 Sr ratios, ɛ Nd is not sensitive to grain size variations (Eisenhauer et al. 1999;Tütken et al. 2002) and has therefore been used extensively as a provenance proxy in the Arctic (e.g., Haley et al. 2008;Fagel et al. 2014;Meinhardt et al. 2016;Bazhenova, Fagel, and Stein 2017). The relatively constant ɛ Nd values throughout Units C-E and further upwards in Unit A are interrupted by two extreme ɛ Nd values in Unit B (Figures 7 and 8). At the base of Unit B (27.6 mblf), an ɛ Nd value of −11.0 suggests an NAC/IP dominated sediment pulse. The topmost sample of Unit B (21.8 mbsf) is characterized by a low radiogenic Nd signal (ε Nd = −14.6) and exhibits the lowest concentration of dolomite (and calcite) throughout the MTW01 sequence ( Figures 5  and 7). Altogether, this suggests a pulse of sediment with a more CS-dominated imprint (Figure 2).  (Asahara et al. 2012;Vonk et al. 2015;Maccali, Hillaire-Marcel, and Not 2018;Schwab et al. 2021). The box plots show the upper quartile (UQ), lower quartile (LQ), median, and range of data points not deemed as outliers. Outliers were identified as data points whose value was >UQ + 1.5 * IQD or <IQ − 1.5 * IQD, where IQD is the interquartile distance (IQD = UQ − LQ). The raw data are overlain to show the distribution. For MTW01, symbols indicate the depositional setting (progradational, transitional, marine). Annotations indicate the ɛ Nd values for sediments from lakes in the east and west of the Mackenzie delta and 87 Sr/ 86 Sr of suspended sediments in major tributaries of the Mackenzie River (from Figure 3).
These excursions occur at interesting stratigraphic points in the borehole. Although the Unit B/C transition has not been dated (O'Regan et al. 2018), Murton et al. (2010) have suggested that this boundary correlates to an erosional surface on the Beaufort Shelf formed during a deglacial meltwater outburst at the onset of the Younger Dryas. The second Nd isotopic excursion occurs at the top of Unit B in a sequence of deformed silty-clay lenses dated to 9180 þ310 À 350 cal. yr BP (Figure 7; O'Regan et al. 2018). It also coincides with a secondary, less extensive glacial meltwater outburst documented and dated by Murton et al. (2010) in the terrestrial stratigraphy of Richards Island and the adjacent coastal plains. Klotsko, Driscoll, and Keigwin (2019) recently identified three meltwater pulses from the Mackenzie River corridor during delgaciation by combining the regional acoustic stratigraphy and X-ray fluorescence scanning elemental abundances of sediment cores from the Beaufort slope. These events were dated at 14.5 to 13.0 cal. ka BP (during the Bølling-Allerød), 12.94 cal. ka BP. and 11.3 cal. ka BP. The oldest event had two phases and is argued to originate from drainage of regional proglacial lakes. The latter two events were attributed to outbursts from Glacial Lake Aggasiz at the onset of the Younger Dryas and Pre-Boreal Oscillation. Süfke et al. (2022), have further refined the chronology and origin of Figure 10. Spatial variations in surface sediment εNd and 87 Sr/ 86 Sr in the Canadian Beaufort Sea (Asahara et al. 2012;Vonk et al. 2015;Maccali, Hillaire-Marcel, and Not 2018;Schwab et al. 2021). The range of MTW01 sediments is shown on the right, alongside relative changes expected between the authigenic and carbonate phases and the pure detrital silicate signal (taken from tests individual extracted phases reported by Asahara et al. (2012) and Schwab et al. (2021). The lower 87 Sr/ 86 Sr ratios and less negative εNd signature of MTW01 sediments can both be partially explained by the influence of authigenic and carbonate phases, which were not removed in our study. meltwater events using authigenic and detrital Pb isotopes ( 208 Pb/ 204 Pb, 206 Pb/ 204 Pb) in three cores from the Amundsen Gulf and Beaufort Sea slope. They documented a near-continuous input of meltwater-derived sediments from the retreating Laurentide ice sheet starting at the onset of the Bølling-Allerød. This began with an "Amundsen Gulf Event," followed by two pronounced meltwater pulses during the Younger Dryas (12.0 and 11.2 cal. ka BP) and a minor event around 9.7 cal. ka BP. Presumably, the Amundsen Gulf Event would be associated with relatively large amounts of fine-grained dolomite, however this has not been quantitatively defined but only qualitatively indicated by elevated Ca/ Sr ratios in X-ray fluorescence scanning data (Süfke et al. 2022). Due to the poor age control, it is not possible to establish whether the Younger Dryas and/or Pre-Boreal Oscillation events are captured and sampled in MTW01. The changes in ɛ Nd seen in the transitional unit are the only notable changes in the isotopic signature of sediments that could be indicative of meltwater pulses. These events were not associated with elevated inputs of dolomite.
Existing studies on ɛ Nd in deglacial and Holocene sediments from the Beaufort Sea are limited to the analysis of bulk sediment leachates used to derive the paleo-seawater signal (Deschamps et al. 2019; AMD0214-02PC; Figure 2). They report highly unradiogenic ɛ Nd (−11.9 to −16.8) in the deglacial and Early to Middle Holocene that they broadly attribute to sediment provenance and inland weathering changes in the Mackenzie drainage basin. Given uncertainties in age models and lack of XRD-based mineral quantification, it remains unclear whether the discrete pulses of dolomite-enriched sediments recorded in cores from the Beaufort Sea slope  are solely attributed to meltwater emanating from the Amundsen Gulf (i.e., the Amundsen Gulf Event) or whether they reflect outflow via the Mackenzie River. Generally, ɛ Nd values become less radiogenic (more negative) as one moves east from the Mackenzie Trough across the Amundsen Gulf and Mc'Clure Strait and through the Canadian Arctic Archipelago (CAA) (Dong et al. 2020). In studying late Quaternary sediments from the Northwind Ridge, Wang et al. (2021) used the less radiogenic ɛ Nd and more dolomite-enriched sediments from around Banks Island as a key feature for discriminating between sediment provenance from the Mackenzie River and the western CAA. Additional work on more proximal deglacial and Holocene records is required to test the robustness of this interpretation, which may have great potential for unraveling the provenance of dolomite rich sediments in the central Arctic.
In studying far-field records from the Arctic Ocean, most authors have looked at the modern signature of Mackenzie River sediments or surface sediments on the shelf to fingerprint possible meltwater events (Not and Hillaire-Marcel 2012;Hillaire-Marcel et al. 2013). A growing body of evidence, including our results from MTW01, highlights the possibility that the mineral and isotopic signature of meltwater pulses will likely differ from a relatively stable and well-defined Mackenzie River baseline signal. This again emphasizes the need to groundtruth mineral and isotopic proxies using records proximal to the major glacial troughs draining the Laurentide ice sheet in the Canadian Arctic.

Conclusions
• Results from this study indicate a relatively stable dolomite, ɛ Nd , and 87 Sr/ 86 Sr signature of finegrained Mackenzie River sediments through the studied deglacial and Holocene sequence. • Dolomite in Holocene marine sediments (Unit A) remained stable at 6.5 ± 0.6 wt% of the <38 µm fraction, with slightly higher abundances during the period of progradation (8.4 ± 0.6 wt %). These are much lower than the dolomiterich layers found in central Arctic cores and discrete intervals in deglacial records from the Beaufort Sea slope. Consequently, it remains plausible that the origin of these dolomite-rich sediments may be from meltwater or ice streaming activity in the glacial troughs found east of the Mackenzie delta. • The average neodymium isotope composition in MTW01 (ɛ Nd = −11.0 to −14.6, average ɛ Nd = −13.2 ± 0.1) is in good agreement with Mackenzie River suspended sediments and modern surface samples of the Beaufort shelf. Two exceptions occur in the transitional unit, deposited immediately after the initial marine transgression at this site. These are not accompanied by increased dolomite content. Although dating of this unit is too uncertain to accurately integrate into regional chronologies of meltwater delivery (Süfke et al. 2022), these results do point towards potentially unique detrital ɛ Nd signatures during meltwater events that could be exploited in future studies. • Distinct increases in 87 Sr/ 86 Sr in the marine Holocene sediments are interpreted as a response to decreased grain size resulting from the more distal location of the Mackenzie River mouth after marine transgression. They do not appear to be driven by a notable change in sediment provenance.

Disclosure statement
No potential conflict of interest was reported by the authors.

Funding
This research was supported by the Swedish Research Council grants DNR-2016-05092 and DNR-2020-04379.