3.1 Sediment core OD1507-03TC-41GC-03PC

Sediment core OD1507-03TC-41GC-03PC is stored at the Oregon State University Marine and Geology Repository at 2.8 ^∘C. Samples were taken approximately every 5 cm for the analysis of different organic constituents and approximately every 10 cm for benthic and planktonic foraminiferal analyses. The core represents a spliced record of a trigger core (TC), a gravity core (GC), and a piston core (PC) recovered in outer Petermann Fjord as part of the Petermann 2015 expedition (Fig. 1) (Table 1). The splice has a total length of 556 cm and was recovered 80 km from the 2015 PG grounding zone at an average water depth of 976 m (Reilly et al., 2019). Based on computed tomography (CT) scans, continuous sections with minimal disturbance were chosen for the splice, and correlation was performed using X-ray fluorescence (XRF) Ti $/$ Ca ratios, CT slice images, and CT numbers (Reilly et al., 2019) (Supplement Table S1). Three sedimentary units are distinguished in OD1506-03TC-41GC-03PC, based on physical properties and XRF elemental data (Reilly et al., 2019). Sedimentary unit 3 (ca. 518–555 cm) is a massive diamict composed of a sandy mud with abundant coarse clasts and XRF Ti $/$ Ca ratios around 0.05 (Reilly et al., 2019). The clasts found in this unit likely do not represent ice-rafted debris (IRD) but instead are related to proximity of the grounding zone (Reilly et al., 2019). Unit 2, a clayey laminated mud with no or very low concentrations of coarser material (IRD), is found between ca. 398–518 cm (Reilly et al., 2019). The topmost lithofacies, unit 1 (0–398 cm), is a bioturbated clayey mud with isolated sand and coarser particles (Reilly et al., 2019). Based on the XRF Ti $/$ Ca ratios and the abundance of coarse material, unit 1 has been divided into three subsections (A (0–53 cm), B (53–164 cm), C (164–398 cm)) (Reilly et al., 2019). Subunit 1C is marked by a decreasing trend from high to intermediate IRD occurrence, with intermediate IRD abundances continuing during subunit 1B, followed by low IRD abundance during subunit 1A (Reilly et al., 2019).

The age model for OD1507-03TC-41GC-03PC, established by Reilly et al. (2019), is based on radiocarbon (¹⁴C) dating of benthic and planktonic foraminifera and calibration of radiocarbon ages using the Marine13 curve (Reimer et al., 2013) and MatCal MATLAB tools (Lougheed and Obrochta, 2016). The ΔR value was constrained using paleosecular variation (PSV) stratigraphy, with the best fit determined for a constant ΔR choice of 770 years (Reilly et al., 2019). No ages were determined for the interval between 408 and 556 cm, due to lack of radiocarbon dates. Regional constraints, however, suggest that these sediments are younger than 7600 years, corresponding to the timing of the retreat of PG into the fjord (Jakobsson et al., 2018).

Table 1Metadata of sediment cores used with the OD1507-03TC-41GC-03PC splice (Reilly et al., 2019).

Download Print Version | Download XLSX

3.2 Sea ice biomarker methodology

Reconstructions of past sea ice conditions in Petermann Fjord rely on the identification of source-specific Arctic sea ice biomarkers, together with the identification of common sterol biomarkers, TOC, and planktonic and benthic foraminiferal abundances. This is a qualitative method to determine past sea ice dynamics. For clarity, we will describe past sea ice conditions using the following categories in order of increasing average sea ice concentration: ice free, reduced seasonal sea ice, enhanced seasonal sea ice, and near-perennial.

Source-specific sea ice biomarkers include a mono- and a di-unsaturated highly branched isoprenoid (HBI), termed IP₂₅ (ice proxy with 25 carbon atoms), and HBI II (Belt, 2018; Belt et al., 2007). IP₂₅ is produced by a number of spring sea ice dwelling diatoms, including Haslea spicula, H. kjellmanii, and Pleurosigma stuxbergii var. rhomboides (Brown et al., 2014). Thus, its presence in Arctic marine sediments provides evidence for past seasonal sea ice occurrence, while the absence of IP₂₅ occurs in year-round ice-free environments as well as under perennial sea ice cover (Belt, 2018; Belt et al., 2007; Brown et al., 2014; Navarro-Rodriguez et al., 2013; Xiao et al., 2015). Salinity changes might exert an additional control on IP₂₅ production, with Ribeiro et al. (2017) showing that its production appears to be suppressed by meltwater in fjords of northeastern Greenland. Given its co-production in H. spicula, H. kjellmanii, and Pleurosigma stuxbergii var. rhomboides (Brown et al., 2014), HBI II co-varies with IP₂₅ in the Arctic realm. Thus, the typically higher HBI II concentrations can provide additional information during times of low/absent sedimentary IP₂₅.

In recent years, biomarker-based reconstructions of past sea ice dynamics have also commonly reported a tri-unsaturated HBI (HBI III), produced by diatoms characteristic of the spring sea ice edge bloom in the marginal ice zone (MIZ) (Belt et al., 2015, 2017; Smik et al., 2016). In Petermann Fjord, however, the identification of the HBI III peak in gas chromatography–mass spectrometry (GC–MS) chromatograms was compromised due to interference with a neighboring peak (Fig. S2). This is specific to sediment samples from Petermann Fjord, as HBI III could be identified in reference sediment from Young Sound extracted alongside the samples from Petermann Fjord (see Sect. 3.3.2). Instead, we use a range of sterol biomarkers (campesterol, cholesterol, brassicasterol, dinosterol, and β-sitosterol) and the abundance of benthic and planktonic foraminifera to gain insights into the regional primary productivity regime and the addition of terrestrial organic carbon to Petermann Fjord. Sterols are common compounds in eukaryotic cell membranes, abundant in both marine and terrestrial organic carbon. This complicates their use as unambiguous tracers of a specific organism and/or environmental regime (Belt and Müller, 2013; Rontani et al., 2014; Volkman, 1986). Nonetheless, the relative abundance of certain sterols can vary according to the predominant organic carbon source (Rontani et al., 2014; Volkman, 1986; Yunker et al., 1995) (Table 2).

Table 2Overview of the dominant sources of different sterol biomarkers measured in this study.

Download Print Version | Download XLSX

Depending on the predominant primary producers, high concentrations of brassicasterol and dinosterol are reported from ice-free regions, the MIZ, and under varying concentrations of seasonal sea ice (Xiao et al., 2015). In general sterol concentrations are low under perennial sea ice in the Arctic Ocean, where all primary productivity is impeded due to limited light availability (Xiao et al., 2015).

Here we use brassicasterol, dinosterol, and cholesterol to represent phytoplankton primary productivity in Petermann Fjord, with different relative concentrations, potentially indicating changes in the ecosystem composition. Campesterol and β-sitosterol, on the other hand, are used to gain insight into terrestrial organic carbon input. The multiproxy study of IP₂₅, HBI II, sterol biomarkers, and the benthic and planktonic foraminiferal abundance (see Sect. 3.4) thus allows interpretations of past sea ice dynamics and terrestrial versus marine organic carbon input to Petermann Fjord.

3.3 Analysis of total organic carbon (TOC), carbon isotopic composition of TOC (δ¹³C_org), and sea-ice-related biomarkers in sediments from OD1507-03TC-41GC-03PC

Prior to the analysis of organic constituents, the samples were freeze-dried (−45 ^∘C; 0.2 mbar; 48 h) at Aarhus University using a Christ Alpha 1-4 LSC freeze drier. The dried samples were weighed, and dry bulk densities (DBDs) were determined from the samples' respective volume and dry weight. Subsequently, the samples were homogenized with a dichloromethane (DCM) cleaned pestle and mortar.

3.4 Planktonic and benthic foraminiferal abundances

The OD1507-03TC-41GC-03PC benthic and planktonic foraminiferal abundances were determined on 58 samples. The benthic foraminiferal counts include both calcareous and agglutinated species. Where sufficient core material was available, the sample depths correspond to the TOC and biomarker samples. The foraminiferal samples were weighed and wet sieved at 63 µm. The > 63 µm fraction was counted wet and submerged in a “storage” solution of 70 % distilled water and 30 % ethanol with baking soda to preserve fragile calcareous and agglutinated tests. A wet splitter was used when necessary to achieve a count of at least 200–300 benthic foraminifera. Planktonic foraminifera were counted in the benthic split. Equivalent dry weights of the foraminiferal samples were calculated using the wet weights of the samples and the wet and dry weights of other samples from the same depths. That way, the numbers of benthic and planktonic foraminifera per gram of dry sediment could be calculated without drying the samples. Foraminiferal fluxes were calculated using the samples DBD and LSR and reported in specimens cm⁻² yr⁻¹. Where the depth of the biomarker and foraminiferal samples were not identical, the DBD was linearly interpolated between neighboring biomarker samples.

4.1 HBI biomarker concentrations

IP₂₅ and HBI II were analyzed in 100 samples with an average depth resolution of 5.5 ± 1.8 cm and temporal resolution of 99 ± 45 years. Both biomarkers were present in all samples, with an overall range of 0.3–30.0 ng g⁻¹ sed for IP₂₅ and 0.6–83.0 ng g⁻¹ of sediment for HBI II and fluxes of 0.2–3.0 and 0.3–8.0 ng cm⁻² yr⁻¹, respectively (Fig. 3f, g). HBI concentrations normalized to the amount of sediment follow the biomarker concentrations normalized to the TOC content of the samples (Fig. S3). A significant positive correlation was found for IP₂₅ and HBI II (R²=0.95 [0.84; 0.97], n= 100), in line with co-production of these two biomarkers (Brown et al., 2014). Sedimentary unit 3 is marked by very low IP₂₅ and HBI II concentrations (Fig. 3f, g). Both biomarker concentrations increase stepwise, with a first increase at 515 cm corresponding to the transition from sedimentary units 3 to 2 (Reilly et al., 2019). The second increase in biomarker concentrations at 405 cm precedes the sedimentary unit boundary (1C–2) by ∼ 5 cm and is followed by an interval of highly variable IP₂₅ and HBI II concentrations and fluxes until 320 cm (Fig. 3f, g). From 260 cm, a steep increase in IP₂₅ and HBI II concentrations and fluxes culminates in peak concentrations at 199 cm followed by peak fluxes at 189 cm (Fig. 3f, g). Subsequently a two-step decrease in all HBI concentrations and fluxes is observed between 199 and 154 cm, spanning the transition of sedimentary units 1C and 1B at ∼ 165 cm (Reilly et al., 2019) (Fig. 3f, g). At 154 cm, a sharp increase in IP₂₅ marks the onset of a second interval with sustained high concentrations between 154 and 70 cm. IP₂₅ fluxes increase simultaneously but decrease prior to the concentrations at 89 cm in connection with a decrease in the LSR (Fig. 3a, g). HBI II concentrations remain relatively low throughout lower sedimentary unit 1B, while the fluxes follow the IP₂₅-flux trend (Fig. 3f, g). At 70 cm, however, IP₂₅ and HBI II concentrations and fluxes drop simultaneously, preceding the transition from sedimentary units 1B to 1A by ∼ 7 cm (Fig. 3f, g). HBI biomarker concentrations and fluxes remain low throughout unit 1A, with a minor increase around 25 cm (Fig. 3f, g).

Figure 3HBI biomarker, TOC, and foraminiferal results from OD1507-03TC-41GC-03PC. From the top to the bottom: (a) linear sedimentation rates (LSR; black), the XRF Ti $/$ Ca ratio (grey), and the CT > 2 mm clast index (brown) (Reilly et al., 2019). (b) TOC carbon isotopic values (δ¹³C_org ‰, VPDB; light green). (c) TOC fluxes (green, filled-in area) and TOC concentrations (green line). (d) Benthic foraminiferal fluxes (purple, filled-in area) and benthic foraminiferal abundance in individuals per gram of sediment (purple line with diamonds). (e) Planktonic foraminiferal fluxes (crimson, filled-in area) and planktonic foraminiferal abundance in individuals per gram of sediment (crimson line with diamonds). (f) HBI II fluxes (turquoise, filled-in area) and absolute concentration (turquoise line with dots) normalized to the amount of extracted sediment (ng g⁻¹ sed). (g) IP₂₅ fluxes (blue, filled-in area) and absolute concentration (blue line with dots) normalized to the amount of extracted sediment (ng g⁻¹ sed). Fluxes do not extend beyond ∼ 400 cm, as this depth corresponds to the lowermost available radiocarbon date (Reilly et al., 2019). The vertical background fill indicates the lithological units (1A–1C, 2, 3) at OD1507-03TC-41GC-03PC (Reilly et al., 2019).

Download

The ratio of HBI II and IP₂₅ (DIP₂₅) has previously been used as an indicator for sea surface temperature (SST) and thus as a tracer of warmer water masses (Hörner et al., 2016; Xiao et al., 2013), in line with higher temperatures being more favorable for the synthesis of double bonds. Other studies, however, did not find a relationship between DIP₂₅ and SST and propose instead that a steady DIP₂₅ reflects stable sea ice conditions, while a variable DIP₂₅ indicates more unstable sea ice conditions (Belt and Müller, 2013; Cabedo-Sanz et al., 2013). In outer Petermann Fjord, the DIP₂₅ ratio is low throughout sedimentary unit 3 and rises before the unit 3–unit 2 boundary (Fig. S4). This is followed by the highest recorded DIP₂₅ values between 520 and 430 cm in unit 2. At 430 cm a sharp decline in DIP₂₅ values marks the onset of a 110 cm long section with relatively low but variable DIP₂₅ values (σ²= 0.41). At 320 cm, a small increase in DIP₂₅ is followed by a relatively steady decline throughout upper unit 1C, unit 1B, and unit 1A (Fig. S4). This coincides with reduced variance in the data (σ²= 0.26).

4.2 Sterol biomarker concentrations

Sterol biomarkers were determined on 94 samples with a depth and temporal resolution of 5.7 ± 2.3 cm and 102 ± 51 years, respectively. We measured brassicasterol, dinosterol, and cholesterol (hereinafter grouped as marine sterols) and campesterol and β-sitosterol (hereinafter grouped as terrestrial sterols). All sterols were present consistently throughout the core (Table 3). There are minor differences in the temporal evolution of sterol concentrations normalized to the amount of sediment extracted and to the TOC content of the samples (Fig. S3). These are especially apparent in the dinosterol, campesterol, and β-sitosterol concentrations. The differences are focused on the lower part of the record (500–400 cm), where Petermann Fjord is influenced by enhanced influx of terrestrial organic carbon (see Sect. 5.1). Thus, we use sterol concentrations normalized to the extracted sediment mass and the resulting fluxes to make inferences about past environmental changes. One sample at ∼ 319 cm has β-sitosterol concentrations 3 times higher than the average, while all other sterol concentrations fall within the reported range (Fig. S5). It is not clear why the β-sitosterol concentrations are so high in this specific sample. Thus, we have excluded this sample from all further analyses and interpretations. Cholesterol is significantly correlated with dinosterol (R²= 0.90 [0.64; 0.98], n= 93) and brassicasterol (R²= 0.93 [0.75; 0.98], n= 93), in line with a predominantly marine source of cholesterol in Petermann Fjord. Further, the marine sterols show a good correlation with the IP₂₅ concentrations (IP₂₅-brassicasterol: R²= 0.82 [0.62; 0.92], n= 93, IP₂₅-cholesterol: R²= 0.89 [0.69; 0.96], n= 93, IP₂₅-dinosterol: R²= 0.82 [0.56; 0.93], n= 93), suggesting an important role of sea ice in the regional marine primary productivity throughout the record.

Similar to the HBIs, the concentrations of all sterols increase at the transition from sedimentary units 3 to 2, with the most significant increase observed in β-sitosterol and campesterol (Fig. 4b, c). Both terrestrial sterols are characterized by a double peak throughout sedimentary unit 2, followed by a decline at the 1C–2 boundary, while the concentrations of marine sterols remain relatively stable throughout this interval (Fig. 4). Throughout unit 1C an increase in the cholesterol and dinosterol concentrations is observed from 350 cm, while their fluxes, alongside all other sterol concentrations, remain low until 260 cm (Fig. 4d, e). Peak concentrations are reached between 199 and 189 cm, followed by a slight decrease in all sterol concentrations and a recovery during lower sedimentary unit 1B (Fig. 4). Simultaneously, sterol fluxes are at their highest between 194 and 94 cm, corresponding to the interval of maximum LSR (Fig. 4). While the sterol fluxes decrease from 94 cm, concentrations peak from 84–74 cm, followed by a sharp decrease prior to the transition from sedimentary units 1B to 1C. Unit 1A is characterized by overall low sterol concentrations and fluxes (Fig. 4).

Figure 4Sterol biomarker results at OD1507-03TC-41GC-03PC. From the top to the bottom: (a) marine sterol index (sum of marine sterols divided by the sum of all sterols) (purple line). (b) β-Sitosterol fluxes (brown, filled-in area) and absolute concentration (brown line with dots) normalized to the amount of extracted sediment (ng g⁻¹ sed). (c) Campesterol fluxes (light brown, filled-in area) and absolute concentration (light brown line with dots) normalized to the amount of extracted sediment (ng g⁻¹ sed). (d) Cholesterol fluxes (orange-brown filled-in area) and absolute concentration (orange-brown line with dots) normalized to the amount of extracted sediment (ng g⁻¹ sed). (e) Dinosterol fluxes (dark green, filled-in area) and absolute concentration (dark green line with dots) normalized to the amount of extracted sediment (ng g⁻¹ sed). (f) Brassicasterol fluxes (light green, filled-in area) and absolute concentration (light green line with dots) normalized to the amount of extracted sediment (ng g⁻¹ sed). Fluxes do not extend beyond ∼ 400 cm, as this depth corresponds to the lowermost radiocarbon date (Reilly et al., 2019). The vertical background fill indicates the lithological units (1A–C, 2, 3) at OD1507-03TC-41GC-03PC (Reilly et al., 2019).

Download

For further insight into the environmental factors driving sterol variability, the marine sterol index (sum of marine sterols/sum of all sterols) was determined (Stein et al., 2017). This indicates that sedimentary units 3, 2, and the lowermost 1C have a higher relative concentration of terrestrial sterols, with a slight increase in the relative concentration of marine sterols just prior to the 2–1C boundary (Fig. 4a). The relative concentration of marine sterols increases between 314 and 59 cm, followed by a decrease at the boundary of sedimentary units 1A–1B (Fig. 4a).

Table 3Ranges of sterol biomarker concentrations and fluxes in OD1507-03TC-41GC-03PC.

Download Print Version | Download XLSX

4.3 Total organic carbon (TOC) concentration and TOC carbon isotopes (δ¹³C_org)

The overall amount of TOC in OD1507-03TC-41GC-03PC is very low, varying between 0.1 wt % and 0.3 wt % (Fig. 3c). Between 394–289 cm and at 109 cm, samples recorded high TOC (0.3 wt % –1.8 wt % ). Based on the available evidence, this suggests incomplete inorganic carbonate removal during sample processing. The interval between 394 cm and 289 cm is associated with high Ti $/$ Ca ratios and IRD delivery to the core site, indicating increased glacial sedimentation in outer Petermann Fjord (Reilly et al., 2019) (Fig. 3a), which might have been associated with input of detrital carbonates. Thus, these samples have been removed from the TOC and δ¹³C_org record (Fig. 3b and c).

The lowest TOC is encountered in sedimentary unit 3, followed by a sharp increase prior to the boundary of sedimentary units 3 and 2 at 520 cm depth. The lower unit 2 is characterized by large variability in the TOC content with two peaks at 520 and 490 cm, while the upper sediments of the same unit have a more stable TOC content around 0.15 wt % (Fig. 3c). From 409 cm a steady increase in TOC into sedimentary unit 1C is observed. The TOC content in sedimentary units 1C, 1B, and 1A varies between 0.17 wt % and 0.27 wt %, with local maxima around 250 and 94 cm and minima at 159 and 45 cm (Fig. 3c). TOC fluxes throughout these units vary between 109 and 422 µg cm⁻² yr⁻¹ with maximum fluxes from 194–94 cm, corresponding to the interval of maximum LSR. After 94 cm TOC fluxes decrease towards the boundary of sedimentary units 1B and 1A with minimum fluxes recorded at 54 cm and flux values similar to unit 1C throughout unit 1A (Fig. 3c).

The δ¹³C_org, measured on the same material as TOC, varies between −24.5 ‰ and −28.2 ‰. The lowest values are observed in sedimentary units 3 and 2 with minima at 530 and 490 cm and a relatively steady rise throughout the latter (Fig. 3b). Sedimentary units 1C, 1B, and 1A are marked by δ¹³C_org values between −26.4 ‰ and −24.5 ‰ with local minima at 249, 159, and 10 cm (Fig. 3b).

4.4 Planktonic and benthic foraminiferal abundances

The benthic and planktonic foraminiferal abundances vary between 0 and 429 specimens per gram of sediment and 0 and 92 specimens per gram of sediment, respectively (Fig. 3d, e). This corresponds to benthic and planktonic fluxes of 3–32 and 0–8 specimens cm⁻² yr⁻¹ (Fig. 3d, e). In sedimentary units 2 and 3, benthic and planktonic foraminiferal abundances are low. Both benthic and planktonic abundances increase near the unit 2–unit 1C boundary, at 404 and 384 cm depth, respectively. Unit 1C is characterized by the highest overall foraminiferal abundances and fluxes, with planktonic abundance and fluxes peaking between 312 and 304 cm prior to the benthic abundance, which peaks at 254 cm (Fig. 3d and e). Planktonic foraminifera abundances and fluxes decline from 250 cm in unit 1C and continue at low values through unit 1B and are nearly absent in unit 1A (Fig. 3e). Benthic foraminiferal abundances and fluxes are more variable but reach their lowest abundance in unit 1A (Fig. 3d).

5.1 Holocene variability in organic carbon sources and sea ice dynamics in Petermann Fjord

Organic carbon (C_org) in Petermann Fjord is derived from a variety of sources, including in situ pelagic and sympagic production, advection of marine organic matter, and input of both modern and ancient C_org from the surrounding landmasses. While HBI concentrations are only marginally influenced by the addition of terrestrial C_org, sterols are important constituents in both marine and terrestrial primary producers (Belicka et al., 2004; Rontani et al., 2014; Volkman, 1986, 2003; Volkman et al., 1993, 2000; Yunker et al., 1995) (Table 2). Thus, in combination with δ¹³C_org, they can provide information on the dominant C_org source.

Sedimentary unit 3 represents grounding zone proximal sedimentation, associated with the deglacial retreat of PG into the fjord (Jakobsson et al., 2018; Reilly et al., 2019). This is in line with very low biomarker and TOC concentrations and near absence of benthic and planktonic foraminifera (Fig. 5a–d), suggesting reduced primary productivity under an ice shelf or thick sea ice cover and/or high sedimentation rates close to the grounding zone diluting the concentrations of organic constituents and marine microfossils close to the margin of a retreating glacier. Both sedimentary units 3 and 2 are beyond the lowermost available radiocarbon age, which is why the influence of sedimentation rates on the biomarker and TOC concentrations cannot be assessed. The low marine sterol index throughout sedimentary unit 3 indicates that the C_org input was likely dominated by terrestrial sources, supported by low δ¹³C_org values, characteristic of terrestrial vegetation in high latitudes (−26 ‰ to −28 ‰) (Ruttenberg and Goñi, 1997) (Fig. 3b). Where Washington and Hall Land border Petermann Fjord they are characterized by sedimentary rocks of lower Paleozoic age, primarily composed of shallow marine carbonates and evaporites (Dawes et al., 2000; Harrison et al., 2011). These contain ancient biomass that can be traced in marine sediments. The thermal maturity of these rocks, however, indicates that fossil sterols will have been, to a large part, degraded to steranes (Parnell et al., 2007), not measured as part of this study. Thus, the increased abundance of terrestrial sterols in sedimentary unit 3 is most likely related to input of fresh terrestrial organic material. Campesterol and β-sitosterol are dominant in vascular plants (such as herbs) but are also found in mosses and lichens (Matsuo and Sato, 1991; Safe et al., 1975), common to the high Arctic tundra. The deglacial retreat of PG occurs relatively late during the deglaciation of northern Greenland (after 7600 cal yr BP) (Jakobsson et al., 2018), falling into the regional Holocene Thermal Maximum (HTM; ca. 11 000–5500 cal yr BP) (Kaufman et al., 2004) when adjacent Washington Land was already deglaciated (Ceperley et al., 2020) and pollen records from northern Greenland and Ellesmere Island suggest a more productive terrestrial Arctic ecosystem (Gajewski, 2015; Mode, 1996). This supports enhanced input of fresh terrestrial C_org to Petermann Fjord in response to glacial erosion and enhanced meltwater drainage of surrounding landmasses during the deglacial retreat of PG.

Figure 5Temporal changes in the environmental conditions in Petermann Fjord. From the top to the bottom: (a) benthic foraminiferal fluxes (purple, filled-in area) and benthic foraminiferal abundance in individuals per gram of sediment (purple line with diamonds). (b) Planktonic foraminiferal fluxes (crimson, filled-in area) and planktonic foraminiferal abundance in individuals per gram of sediment (crimson line with diamonds). (c) Brassicasterol fluxes (light green, filled-in area) and absolute concentration (light green line with dots) normalized to the amount of extracted sediment (ng g⁻¹ sed). (d) IP₂₅ fluxes (blue, filled-in area) and absolute concentration (blue line with dots) normalized to the amount of extracted sediment (ng g⁻¹ sed). (e) Marine sterol index (sum of marine sterols divided by the sum of all sterols) (purple line). (f) TOC fluxes (green, filled-in area) and TOC concentrations (wt %, green line). (g) XRF Ti $/$ Ca ratios (grey) and the CT > 2 mm clast index (brown) (Reilly et al., 2019). (h) Reconstruction of the Holocene ice tongue extent (minimum extent in dark blue, maximum extent in light blue) in Petermann Fjord in kilometers from the grounding zone (Reilly et al., 2019). OD1507-03TC-41GC-03PC is located 80 km from the present-day grounding zone (black dashed line). Lithogenic units (1A–C, 2, 3) in sediment core OD1507-03TC-41GC-03PC are indicated with dashed vertical lines; grounding zone proximal sedimentation (unit 3) is indicated with a grey vertical bar. The red triangles on the x axes indicate the timing of the ice tongue breakup (∼ 6900 cal yr BP) and the late Holocene inception of a small ice tongue (∼ 2200 cal yr BP) (Reilly et al., 2019).

Download

Mapped submarine landforms suggest that the retreat of PG into the fjord was rapid, driven by ice cliff instability and promoted by the retrograde slope of the outer fjord sill (Jakobsson et al., 2018). An inner sill (Tinto et al., 2015) ∼ 30 km seaward of the present-day grounding zone (Reilly et al., 2019) likely halted the retreat, allowing for the formation of an extensive ice tongue associated with the laminated, IRD-poor lithofacies of unit 2 (Reilly et al., 2019). The transition from units 3 to 2 is marked by an increase in all biomarkers, TOC, and benthic and planktonic foraminiferal abundances (Fig. 5), suggesting an overall increase in C_org input to outer Petermann Fjord and/or lower sedimentation rates in response to the larger distance of the core site from the grounding zone. The marine sterol index remains low, suggesting continued dominance of terrestrial versus marine sterols, while increasing δ¹³C_org values indicate enhanced importance of marine C_org sources to Petermann Fjord. Thus, compared to sedimentary unit 3, the organic matter in sedimentary unit 2 was likely derived from more varied sources, including in situ/advected marine C_org and terrestrial C_org. This supports a marginal/sub-ice tongue regime in outer Petermann Fjord with no or only intermittent sympagic and pelagic primary productivity in the fjord and low food supply to sustain benthic productivity (Jennings et al., 2020).

The early Holocene ice tongue of PG collapsed around 6900 cal yr BP (unit 2–unit 1C boundary), marked by the abrupt appearance of IRD clasts in sediments across Petermann Fjord (Reilly et al., 2019) (Fig. 5g). This is associated with an increase in the marine sterol index and TOC (Fig. 5), suggesting an increase in the C_org delivery to outer Petermann Fjord and relatively more input of marine versus terrestrial sterols compared to previous sedimentary units. A second, larger increase in the marine sterol index is evident at 5100 cal yr BP (Fig. 5e), contemporaneously with decreasing XRF Ti $/$ Ca ratios and IRD flux (Reilly et al., 2019). This indicates reduced glacial erosion and calving activity or a decrease in the delivery of erosional products to outer Petermann Fjord from 5100 cal yr BP and dominant input of marine over terrestrial C_org, supported by δ¹³C_org (Figs. 5, 3). The marine sterol index remains high until ∼ 600 cal yr BP, while δ¹³C_org is high throughout the rest of the core, suggesting primarily delivery of marine-derived organic matter to outer Petermann Fjord during the middle to late Holocene.

The breakup of the early Holocene ice tongue is also associated with a small increase in the sedimentary IP₂₅ and HBI II concentrations, as well as an increase in the benthic and planktonic foraminiferal abundances. Between 6900 and 5500 cal yr BP, IP₂₅ fluxes are low but variable, indicating an unstable/variable sea ice regime with predominantly reduced seasonal sea ice concentration during spring and low rates of sympagic productivity. An unstable sea ice regime is further supported by variable DIP₂₅ values (Fig. S4). Alternatively, increased meltwater runoff related to the collapse of the ice tongue and associated retreat of PG might have caused low/variable IP₂₅ fluxes by creating environmental conditions unfavorable for IP₂₅-producing diatom species during times of increased freshwater discharge (Ribeiro et al., 2017). The simultaneous increase in TOC, benthic, and especially planktonic foraminiferal fluxes, however, supports a reduced seasonal sea ice cover and a shift from a regime marginal to or below an ice tongue towards ameliorated conditions with a prolonged open water season and enhanced pelagic primary productivity in the fjord, allowing planktonic foraminifera to thrive (Fig. 5b). A reduced seasonal sea ice cover is further supported by evidence of low sea ice concentrations in Hall Basin until at least 6000 cal yr BP (Jennings et al., 2011b). Regional records from northern Greenland and the western Canadian Arctic Archipelago (CAA) demonstrate a HTM between 11 000 and 5500 cal yr BP (Belt et al., 2010; Briner et al., 2016; England et al., 2008; Funder et al., 2011; Jennings et al., 2011b; Kaufman et al., 2004; Knudsen et al., 2008; Lecavalier et al., 2017; Ledu et al., 2010; Vare et al., 2009), associated with regional mean annual surface air temperatures 3 ± 1 ^∘C higher than pre-industrial (1750 CE) (Lecavalier et al., 2017). Thus, reduced seasonal sea ice concentrations in Petermann Fjord between 6900 and 5500 cal yr BP likely represent the late stages of the HTM in the northern Nares Strait region.

From ca. 5800 cal yr BP the increase in planktonic and benthic foraminiferal fluxes steepens, associated with less variable but still low IP₂₅ fluxes from 5500 cal yr BP (Fig. 5a, b, and d). While benthic foraminifera respond to sea ice changes via its influence on marine productivity and the resulting changes in food supply to the seafloor (Seidenkrantz, 2013), planktonic foraminiferal abundances have been shown to be highest in the open water region and along the ice margin, with only few individuals occurring under persistent sea ice (Carstens et al., 1997; Mayot et al., 2020; Pados and Spielhagen, 2014). In the modern environment of the northern Nares Strait and Petermann Fjord, planktonic foraminiferal abundances are very low in the fjord and much higher in the mobile sea ice regime of Nares Strait (Jennings et al., 2020), which is characterized by a shorter seasonal sea ice season than outer Petermann Fjord (Fig. 2). Thus, high planktonic foraminiferal fluxes between 5800 and 3600 cal yr BP suggest sustained periods of seasonally open waters during summer, while the continuously low IP₂₅ fluxes are consistent with reduced seasonal sea ice concentrations and the absence of a sympagic spring bloom in outer Petermann Fjord (Fig. 5d). Alternatively, IP₂₅ production might have been suppressed due to lasting meltwater discharge into the fjord, resulting from the influence of the retreating PG on the outer fjord environment. However, both the > 2 mm clast index and the XRF Ti $/$ Ca ratios decrease throughout this interval (Reilly et al., 2019), suggesting a gradual reduction of the influence of glacial activity on the outer fjord (Fig. 5g).

From 3900 cal yr BP, IP₂₅ fluxes increase steeply, accompanied by an increase in all sterol fluxes (Figs. 4, 5d). This suggests a shift in the ecosystem in outer Petermann Fjord, associated with a transition from a regime dominated by pelagic primary productivity towards a regime characterized by enhanced sympagic productivity. Especially after 3600 cal yr BP, rapidly decreasing planktonic foraminiferal fluxes and steadily increasing IP₂₅ fluxes indicate a progressively longer sea ice season and enhanced sympagic productivity in outer Petermann Fjord (Fig. 5b, d). This falls into a period of long-term declining regional atmospheric temperatures recorded at Agassiz ice cap (Lecavalier et al., 2017) and in lake records in NW Greenland (Axford et al., 2019; Lasher et al., 2017). Neoglacial cooling has been observed in numerous marine and terrestrial archives around Greenland and the wider North Atlantic region (e.g., England et al., 2008; Hansen et al., 2020; Jennings et al., 2011a; Limoges et al., 2020; Vare et al., 2009), as a response to decreasing Northern Hemisphere summer insolation (Marcott et al., 2013). Thus, enhanced seasonal sea ice conditions in Petermann Fjord from 3900 cal yr BP (Fig. 5d) likely record the onset of Neoglacial cooling in the northern Nares Strait region.

Peak IP₂₅ fluxes around 2500 cal yr BP are associated with high fluxes of marine and terrestrial sterols as well as benthic foraminifera, while planktonic foraminiferal fluxes are low (Figs. 4, 5a, b, d), indicating a prolonged seasonal sea ice cover with only short periods of open water during summer. From 2500–2100 cal yr BP a two-stepped decrease in IP₂₅ fluxes is observed, of which the second decline is accompanied by a decrease in the TOC and benthic foraminiferal flux and a small decrease in the (already low) planktonic foraminiferal flux (Fig. 5a, b, f). A contemporaneous decline in all productivity indicators and sea ice biomarkers suggests a restriction in the pelagic and sympagic primary productivity alike, most likely as a response to further lengthening of the sea ice season to near-perennial sea ice cover. This interval precedes the late Holocene inception of a small ice tongue in Petermann Fjord at 2200–2100 cal yr BP, inferred from Ti $/$ Ca ratios and the stacked > 2 mm clast index from four cores in Petermann Fjord (Reilly et al., 2019). After 2100 cal yr BP, sea ice biomarkers, TOC, and benthic foraminiferal fluxes recover, while the planktonic foraminiferal abundance remains low (Fig. 5a, d, f), indicating an ecosystem dominated by sympagic productivity and enhanced seasonal sea ice cover, similar to the interval 3600–2500 cal yr BP. Compared to the deglacial PG ice tongue, the late Holocene ice tongue (after 2100 cal yr BP) (Reilly et al., 2019) does not seem to be associated with increased input of terrestrial organic matter to outer Petermann Fjord (Fig. 5). A possible explanation could be lower atmospheric temperatures compared to the early Holocene (Lecavalier et al., 2017), which were associated with a less diverse and more sparse terrestrial flora in the high Arctic (Gajewski, 2015) and decreased meltwater input. Another reason could be increased distance of OD1507-03TC-41GC-03PC from the PG grounding zone during the late Holocene (Reilly et al., 2019), resulting in reduced delivery of meltwater-derived C_org.

Around 1300 cal yr BP a sharp decline in IP₂₅, sterol, and TOC fluxes is observed, while the concentration of IP₂₅ decreases more gradually, and the concentrations of sterols and TOC increase (Fig. 5). The sharp decrease in fluxes at this time coincides with a large drop in the LSR (Fig. 3a), which might be biasing the flux data. Instead, the decline in biomarker fluxes at 950 cal yr BP, accompanied by a decrease in the biomarker concentrations, seems to be a more reliable feature (Fig. 5). This is associated with declining benthic foraminiferal fluxes and followed by decreases in the planktonic foraminiferal and TOC fluxes at 700 cal yr BP (Fig. 5a, b, f), suggesting a return to near-perennial sea ice conditions with reduced pelagic and sympagic primary productivity, similar to the interval from 2500–2100 cal yr BP. At ca. 600 cal yr BP a rapid extension of the Petermann ice tongue to its modern limits (Reilly et al., 2019) resulted in (at least intermittent) cover of the core site, in line with low biomarker and foraminiferal fluxes between 600 cal yr BP and the top of the core (Fig. 5). The latter indicates low rates of primary productivity in a fjord, which is nearly completely covered by an ice tongue. This is further supported by a decrease in the marine sterol index between 600 cal yr BP and the top of the core, suggesting a relative decrease in the input of marine organic matter (Fig. 5e). However, while there is a small decrease in the TOC at this time, there is no corresponding decrease in the δ¹³C_org (Fig. 3b), suggesting that marine C_org was still the main source of TOC in outer Petermann Fjord.

5.2 Nares Strait sea ice dynamics over the last 7000 cal yr BP

Thus far, only one other biomarker-based sea ice reconstruction, from station Kane2b in northwestern Kane Basin, exists in Nares Strait for comparison with our records (Fig. 1) (Georgiadis et al., 2020). We propose that depending on the ice arch configuration in Nares Strait, Kane Basin and outer Petermann Fjord likely experience opposing conditions related to the proximity of the ice edge during spring–early summer (Fig. S1), which is the dominant productivity season of sea ice biomarkers. Kane2b is located near/under the southern ice arch in Nares Strait. Thus high sympagic productivity and IP₂₅ fluxes occur during times of a stable ice arch in Smith Sound, and sea ice edge conditions exist during spring–summer (Georgiadis et al., 2020) (Fig. 6b, d). Georgiadis et al. (2020) interpret HBI III as generally indicating ice-laden/fresh surface waters, which can occur following the breakup of the southern ice arch and during times of mobile sea ice in Nares Strait. Thus, the southern ice arch scenario is likely associated with variable HBI III fluxes at Kane2b, depending on the seasonal timing of ice arch breakup. During the southern ice arch scenario, sea ice in Petermann Fjord does not break up until late summer–early fall (Fig. 2), likely hindering a pronounced in-ice bloom related to MIZ conditions during spring–summer, resulting in relatively low sea ice and phytoplankton biomarker concentrations (Fig. 6b, d). In years when only the northern ice arch forms, sea ice formed locally in Nares Strait either remains mobile throughout the winter or breaks up during early spring. Under these conditions HBI III fluxes at Kane2b can be relatively high depending on the sea ice flux through Nares Strait during spring–summer, but no pronounced sympagic spring bloom occurs in Kane Basin (Fig. 6c). Instead, outer Petermann Fjord experiences spring MIZ conditions, as the formation of landfast ice in Petermann Fjord is independent of the formation of landfast ice in Nares Strait (Figs. 2, S1). This is likely associated with a significant spring sympagic bloom in outer Petermann Fjord and enhanced primary productivity related to the vicinity of the ice edge, resulting in increased concentrations of sea ice and primary productivity biomarkers in outer Petermann Fjord (Fig. 6c). Periods with contemporaneously low concentrations of sea ice and primary productivity biomarkers in Kane Basin and Petermann Fjord point towards low spring sea ice concentration in the entire Nares Strait, likely associated with a failure of both ice arches (Fig. 6a). In addition to the Kane Basin record, sea ice dynamics around Ellesmere Island and NE Greenland have been inferred from records of driftwood delivery and beach ridge formation (England et al., 2008; Funder et al., 2011), providing context with respect to Holocene sea ice conditions in the Arctic Ocean and Lincoln Sea. Driftwood is transported with Arctic multiyear ice and is deposited along the coastlines of northern Greenland and Ellesmere Island when landfast ice breaks up during summer (Funder et al., 2011). The formation of beach ridges also depends on sufficient wave action and open water along the coast. Thus, abundant driftwood delivery to Ellesmere Island and northeastern Greenland together with abundant formation of beach ridges is indicative of seasonally open waters along the coast (Fig. 6a). Sea ice conditions with year-round landfast ice along the coast, on the other hand, will result in little or no driftwood landings and reduced formation of beach ridges (Fig. 6d). Lastly, variable ice conditions in Lincoln Sea result in occasional/little driftwood landings and reduced/variable formation of beach ridges (Fig. 6b, c).

Figure 6Schematic of different ice arch scenarios in Nares Strait and the resulting spring sea ice conditions and sea ice biomarker (IP₂₅) and primary productivity indicator (Phyto.Marker: HBI III at Kane2b, sterols at OD1507-03TC-41GC-03PC) sedimentation patterns at OD1507-03TC-41GC-03PC (yellow star) and AMD14-Kane2b (red star). Biomarker concentrations are illustrated as relative changes (in order of increasing concentrations: $-+$, $++$, $+++$; $-/+$ indicates variable concentrations). Further, relative changes in the driftwood delivery to CMI and beach ridge formation along the coast of northeastern Greenland are illustrated, based on the different ice arch scenarios. In addition to the access to the coast, driftwood delivery also depends on the multiyear ice conditions in the Arctic Ocean; this is not considered in this simplified schematic. (a) No ice arch formation in Nares Strait resulting in year-round mobile sea ice. (b) Recurrent southern ice arch with landfast sea ice in Nares Strait and an open NOW. (c) Recurrent northern ice arch with locally formed sea ice remaining mobile in Nares Strait year-round. (d) Recurrent northern and southern ice arch with landfast sea ice in Nares Strait and a stable NOW. The time periods in the top left corner of each panel correspond to the proposed Holocene interval that experienced the respective conditions.

While the sediment core from station Kane2b covers the last 9000 cal yr BP (Georgiadis et al., 2020) (Fig. 1), Petermann Fjord did not deglaciate until ∼ 7600 cal yr BP (Jakobsson et al., 2018), and the age model for OD1507-03TC-41GC-03PC is only constrained for the last 7000 cal yr BP. Thus, we focus on the interval from 6900 cal yr BP to the present. Between 6900 and 5500 cal yr BP, corresponding to the late stages of the regional HTM (Kaufman et al., 2004), sea ice biomarker and foraminiferal fluxes in Petermann Fjord and Kane Basin (Georgiadis et al., 2020) indicate reduced seasonal sea ice occurrence (Fig. 7f, g). This corresponds to sparse driftwood delivery but abundant beach ridge formation around NE Greenland from 8500–6000 cal yr BP, suggesting a minimum in Arctic multi-year ice and seasonally open water along the coast (Funder et al., 2011; Möller et al., 2010). Seasonally open water is further supported by maximum driftwood delivery to Disraeli Fiord and Clements Markham Inlet (CMI) on Ellesmere Island (England et al., 2008) (Figs. 1, 7e). This is in line with reduced seasonal sea ice occurrence in Nares Strait, suggesting predominantly mobile sea ice with no or only occasional ice arch formation and export of Arctic sea ice through Nares Strait (Fig. 6a) (Georgiadis et al., 2020).

Figure 7Holocene environmental changes in the wider Nares Strait region. From the top to the bottom: (a) regional mean annual air temperature (MAAT) anomaly based on δ¹⁸O at Agassiz ice cap at 25-year resolution (Lecavalier et al., 2017) (grey line) including a 2σ uncertainty envelope (light grey filled-in area) and a 5 pt running mean (red line), together with the 21 June insolation at 80^∘ N (dark grey). (b) Inferred variability in little auk numbers in seabird colonies at Annikitsoq and Qeqertaq (Davidson et al., 2018). (c) Cumulative probability distribution of calibrated ¹⁴C ages of driftwood in NE Greenland north of 80^∘ N (pink filled-in area) (Funder et al., 2011). The numbers above the shaded area indicate the number of ages contributing to each block. (d) Cumulative probability distribution of calibrated ¹⁴C ages of driftwood at Ward Hunt Ice Shelf and in Disraeli Fiord (yellow filled-in area) (England et al., 2008). The numbers above the shaded area indicate the number of ages contributing to each block. (e) Cumulative probability distribution of calibrated ¹⁴C ages of driftwood in Clements Markham Inlet (yellow filled-in area) (England et al., 2008). The numbers above the shaded area indicate the number of ages contributing to each block. (f) IP₂₅ fluxes (red, filled-in area) and absolute concentrations (red line with diamonds) normalized to the amount of sediment extracted at AMD14-Kane2b (Georgiadis et al., 2020). (g) IP₂₅ fluxes (blue, filled-in area) and absolute concentrations (blue line with dots) normalized to the amount of sediment extracted at OD1507-03TC-41GC-03PC. (h) Reconstruction of the Holocene ice tongue extent (minimum extent in dark blue, maximum extent in light blue) in Petermann Fjord in kilometers from the grounding zone (Reilly et al., 2019). OD1507-03TC-41GC-03PC is located 80 km from the present-day grounding zone (black dashed line). Corresponding to Fig. 6, interpretations of past ice arch dynamics in Nares Strait are marked with dashed vertical lines. The letters at the top correspond to the schematic scenarios in Fig. 6. The red triangles on the x axes indicate the timing of the ice tongue breakup (∼ 6900 cal yr BP) and the late Holocene inception of a small ice tongue (∼ 2200 cal yr BP) (Reilly et al., 2019).

Download

At 5500 cal yr BP increasing sea ice biomarker fluxes in Kane Basin indicate a shift towards later sea ice retreat and ice edge productivity, interpreted to reflect recurrent formation of the southern ice arch in Smith Sound and Kane Basin between 5500 and 3000 cal yr BP (Georgiadis et al., 2020) (Fig. 7f). Studies of seabird colonies in the NOW region suggest the arrival of little auk colonies in NW Greenland at 4400 cal yr BP, supporting the opening of the NOW lee of the Smith Sound ice arch during spring–summer (Davidson et al., 2018; Ribeiro et al., 2021) (Fig. 7b). Little auk are zooplanktivore, feeding on the abundant copepod population of the NOW. Thus, large colonies of little auk in Greenland are only found in the vicinity of productive polynyas, where open water is available for foraging upon their arrival in spring (Davidson et al., 2018). Productive and strong NOW conditions from 4400 cal yr BP are also inferred based on foraminiferal assemblages from the central polynya region (Jackson et al., 2021). Simultaneously, landfast sea ice started to form in Disraeli Fiord on northern Ellesmere Island from ∼ 5500 cal yr BP, covering large parts of the coast by ∼ 3500 cal yr BP (England et al., 2008), and shorter periods of open water and restricted beach ridge formation occur around NE Greenland (Funder et al., 2011). Even though this suggests an intensification of landfast sea ice in the Lincoln Sea and the Arctic Ocean, open water along the NE Greenland coast was more abundant compared to today until at least 4500 cal yr BP (Funder et al., 2011). The period from 5500–3500 cal yr BP thus marks the transition from early Holocene warmth to Neoglacial cooling in the wider Nares Strait region. Compared to Kane2b, biomarker fluxes do not increase significantly at 5500 cal yr BP in outer Petermann Fjord (Fig. 7g). Instead, benthic and planktonic foraminiferal fluxes are at their highest between 5500 and 3500 cal yr BP (Fig. 6a, b), interpreted to reflect reduced seasonal sea ice concentrations, the absence of a sympagic spring bloom, prolonged periods of seasonally open waters during summer, and increased pelagic primary productivity. Alternatively, the continuous influence of the retreating PG during this interval, potentially associated with increased meltwater discharge into the fjord, might have caused low concentrations of sea ice biomarkers in outer Petermann Fjord (Ribeiro et al., 2017). Nonetheless, the high foraminiferal fluxes indicate conditions that differ from what would be expected in outer Petermann Fjord under a stable southern ice arch scenario (Fig. 6b). A possible reason for this could be enhanced seasonality, in particular enhanced winter cooling, due to the increasing sea ice extent in the Arctic Ocean (England et al., 2008; Funder et al., 2011), while summer insolation was still relatively high (though decreasing) compared to the late Holocene (Fig. 7). Modeling studies have shown that the loss of sea ice in the Arctic Ocean during the early Holocene counteracted the increased seasonality prescribed by insolation forcing, due to enhanced ocean–atmosphere heat flux during winter (Fischer and Jungclaus, 2011). The Arctic Ocean thus acted as a heat reservoir with increased latent and sensible heat flux during winter as a result of the reduced sea ice cover (Fischer and Jungclaus, 2011). Conversely, decreasing summer insolation and the associated increase in Arctic Ocean sea ice cover during the mid-Holocene (Fig. 7) strengthened the insulating effect of sea ice on the ocean and led to pronounced cooling during fall–winter (Fischer and Jungclaus, 2011). Enhanced seasonality could explain the observed proxy patterns in Nares Strait between 5500 and 3500 cal yr BP, with an early seasonal breakup of the southern ice arch, leading to increased pelagic primary productivity in northern Nares Strait. This way, the spring ice edge bloom might still have occurred in the southern Nares Strait, followed by a rapid sea ice retreat and open water conditions during summer–fall, as recorded in outer Petermann Fjord. An early seasonal breakup of the southern ice arch is also supported by high HBI III fluxes, especially between 4500 and 3500 cal yr BP (Georgiadis et al., 2020). Thus, we suggest that this interval represents the transition from reduced sea ice conditions and the lack of ice arches in Nares Strait during the late HTM towards more stable sea ice conditions associated with the seasonal formation of a recurrent southern ice arch from at least ∼ 4400 cal yr BP.

From 3900 cal yr BP, increasing IP₂₅ fluxes in Petermann Fjord attest to a shift in the primary productivity regime from predominantly pelagic to MIZ and sympagic primary productivity at 3500 cal yr BP. A similar increase in IP₂₅ concentrations from ∼ 4000 cal yr BP is also observed in the western CAA, in the Barrow, Victoria, and Dease straits (Belt et al., 2010; Vare et al., 2009), in line with enhanced seasonal sea ice formation and the onset of regional Neoglacial cooling. Conversely, IP₂₅ fluxes in Kane Basin decrease from ∼ 4000 cal yr BP until 1400 cal yr BP (Fig. 7f), suggesting no stable southern ice arch in Smith Sound. Instead, the driftwood delivery to CMI ceases at 3500 cal yr BP (England et al., 2008). This likely indicates a regime shift in Nares Strait, with the northern ice arch becoming more prominent between 3500 and 1400 cal yr BP (Fig. 6c). In years where only the northern ice arch forms, along-strait winds in Nares Strait keep locally formed sea ice mobile year-round (Vincent, 2019), creating MIZ conditions in Petermann Fjord during spring–summer (Fig. S1) in line with high IP₂₅ concentrations. Between 2500 and 2100 cal yr BP biomarker fluxes in outer Petermann Fjord suggest a shift towards near-perennial sea ice conditions (Fig. 7g). This coincides with some driftwood landings in CMI (England et al., 2008) and a period of no driftwood or beach ridge occurrences in NE Greenland, suggesting perennial landfast ice (Funder et al., 2011). The numbers of little auks at Qeqertaq (Salve Ø) are variable during this time period, with a local peak in abundance at ∼ 2200 cal yr BP (Davidson et al., 2018). In Kane Basin, IP₂₅ fluxes are relatively low until ca. 1400 cal yr BP; however a minor peak (more pronounced in IP₂₅ concentrations) can be observed between ca. 2300 and 2000 cal yr BP. This coincides with a period of low HBI III fluxes from ca. 2700–2100 cal yr BP at Kane2b (Georgiadis et al., 2020). Georgiadis et al. (2020) interpret HBI III as indicating ice-loaded fresh surface waters due to mobile sea ice/ice arch breakup in Nares Strait. A decrease in HBI III fluxes and concomitant small increase in IP₂₅ could thus result from intermittent formation of the southern ice arch in Kane Basin and opening of the NOW, consistent with the local abundance peak of little auks and sea ice records from outer Petermann Fjord and northeastern Greenland. Foraminiferal assemblages from the central NOW region (91-039-008P, Fig. 1) demonstrate a shift towards increased bottom water ventilation around 2600 cal yr BP (Knudsen et al., 2008), which could point towards strengthened polynya conditions and intensified brine formation in the NOW region. Unfortunately, however, no sediments were recovered at site 91-039-008P between 2450 and 1100 cal yr BP (Knudsen et al., 2008). Increased bottom water ventilation is also inferred from measurements of total sulfur at CASQ1 (Fig. 1) (Jackson et al., 2021). Simultaneously, calcareous foraminiferal abundances are low between 2500–2000 cal yr BP at CASQ1 (Fig. 1) (Jackson et al., 2021). While this hinders unambiguous environmental interpretations, it could also be a sign of strong polynya conditions and formation of CO₂-rich brines, causing dissolution of the calcareous foraminiferal test and preservation of agglutinated species only (Jackson et al., 2021). Georgiadis et al. (2020) suggest that the period of mobile sea ice in Nares Strait between 3500 and 1400 cal yr BP was driven by a change in the dominant phase of the Arctic Oscillation (AO). The shift towards a positive AO between 3000 and 1200 cal yr BP (Darby et al., 2012) might have been associated with stronger along-strait winds in Nares Strait. As ice arch formation in Nares Strait is a function of ice thickness, local wind stress, and atmospheric temperatures (Barber et al., 2001; Samelson et al., 2006), distinct atmospheric cooling spikes between ∼ 2500 and ∼ 1900 cal yr BP, recorded at Agassiz ice cap (Lecavalier et al., 2017), might have been responsible for the enhanced sea ice concentration and the formation of an intermittent southern ice arch in Nares Strait between 2500 and 2100 cal yr BP (Fig. 7a).

After 2100 cal yr BP and until 1400 cal yr BP, the combined biomarker records in Nares Strait support a return to a more northern ice-arch-dominated sea ice regime. This is supported by records from the central NOW region, suggesting weaker polynya conditions from ca. 2000–1200 cal yr BP (Jackson et al., 2021; Ribeiro et al., 2021). From 1400 cal yr BP onwards, IP₂₅ fluxes in Kane Basin increase (Georgiadis et al., 2020), while they decline in Petermann Fjord (Fig. 7f, g) concomitant with high numbers of little auks at Qeqertaq (Salve Ø) (Davidson et al., 2018), suggesting a transition towards a more stable southern ice arch. These conditions are further consolidated from 950 cal yr BP, where the sudden drop in sea ice and primary productivity indicators in Petermann Fjord (Figs. 5, 7g) indicate near-perennial landfast sea ice in Nares Strait promoted by the recurrent formation of a southern ice arch, which was associated with stable NOW conditions from ∼ 800 cal yr BP (Ribeiro et al., 2021). Driftwood is absent in CMI from 2200–400 cal yr BP (England et al., 2008), suggesting that the northern arch might have been a stable feature as well (Figs. 7e, 6d). Periods of double ice arching in Nares Strait are particularly stable (Ryan and Münchow, 2017), in line with the coldest Holocene temperatures recorded at Agassiz ice cap throughout this interval (Lecavalier et al., 2017). At ∼ 600 cal yr BP a large jump in the extent of the late Holocene ice tongue in Petermann Fjord to its modern (pre-2010) limits, leads to a regime marginal to or under an ice tongue in outer Petermann Fjord (Reilly et al., 2019) (Fig. 7h). Thus, further interpretations of the sea ice environment in Nares Strait are hindered due to the influence of the ice tongue on primary productivity in Petermann Fjord.

5.3 Implications of sea ice for the formation and stability of the Petermann ice tongue

A recent study by Reilly et al. (2019) demonstrates that an extensive ice tongue in Petermann Fjord is a relatively recent feature compared to large parts of the Holocene. The final stage of the deglaciation in Petermann Fjord was characterized by the breakup of a deglacial ice tongue around 6900 cal yr BP, followed by a period of enhanced IRD and Ti $/$ Ca ratios in the fjord sediments (6900–3500 cal yr BP; Fig. 5g), indicating increased glacial sedimentation in outer Petermann Fjord. First at ∼ 2200 cal yr BP a small ice tongue re-emerged, followed by gradual growth and a rapid expansion to its modern extent around 600 cal yr BP (Reilly et al., 2019). Thus, the middle to late Holocene was marked by a ∼ 4700-year interval where no stable ice tongue existed in Petermann Fjord (Reilly et al., 2019). The presence of an ice tongue in Petermann Fjord is determined by the interplay of several factors, including the surface mass balance of PG, the loss of glacial ice by calving, and the basal melt rates. Local changes in sea ice dynamics have the potential to influence both the basal melt rates and the stability of the calving front at PG. The latter depends on the formation and seasonal duration of landfast sea ice in Petermann Fjord, creating an ice mélange that stabilizes the calving front and reduces the length of the calving season (Amundson et al., 2010; Carr et al., 2015; Robel, 2017; Todd and Christoffersen, 2014). Along with the increase in IP₂₅ fluxes in outer Petermann Fjord, indicating a transition towards enhanced sea ice conditions, from at least 3500 cal yr BP, a cessation of the IRD flux is observed in OD1507-03TC-41GC-03PC (Fig. 5d, g) (Reilly et al., 2019). This suggests a gradual decrease in the calving of icebergs from PG during the mid-Holocene. We propose that this was partly in response to the formation of seasonal landfast ice in the fjord, shortening the calving season (Figs. 6, 7). Sea ice also has the ability to influence the oceanic heat transport to Petermann Fjord by regulating the wind stress at the atmosphere–ocean interface in Nares Strait, which in turn modifies the circulation-driven inflow of AW to the fjord (Shroyer et al., 2017). Under modern conditions, AW inflow is enhanced at times of mobile sea ice in Nares Strait, while landfast sea ice results in decreased oceanic heat flux to Petermann Fjord (Shroyer et al., 2017). In addition to Nares Strait sea ice conditions, subglacial runoff and glacial isostatic uplift of the outer fjord sill likely also influenced the inflow of AW to Petermann Fjord across the Holocene (Bennike, 2002; Cai et al., 2017; Washam et al., 2019). Recurrent formation of a southern ice arch and seasonal formation of landfast sea ice may have persisted in Nares Strait from 5500–3500 cal yr BP, 2500–2100 cal yr BP, and after 1400 cal yr BP, based on the assessment of regional sea ice records (England et al., 2008; Funder et al., 2011; Georgiadis et al., 2020; this study) (Fig. 6b, d). The interval from 5500–3500 cal yr BP marks the transition from the HTM to Neoglacial cooling in the Nares Strait region. While records from the southern Nares Strait region suggest an opening of the NOW and spring ice edge conditions in Kane Basin from at least 4400 cal yr BP (Davidson et al., 2018; Georgiadis et al., 2020; Jackson et al., 2021; Ribeiro et al., 2021), records from outer Petermann Fjord indicate low sympagic but high pelagic productivity, suggesting prolonged periods of open water during summer (Fig. 5a, b, d). Increased seasonality, with early breakup of the southern ice arch and late seasonal sea ice formation, might be able to explain the observed patterns (Sect. 5.3). Thus, unlike today, the formation of a southern ice arch between 5500 and 3500 cal yr BP is unlikely to have resulted in a prolonged landfast ice season in Nares Strait. Consequently, the longer mobile sea ice season in Nares Strait might have contributed to AW inflow to Petermann Fjord, preventing the formation of an ice tongue throughout this interval.

Between 2500 and 2100 cal yr BP regional records of sea ice (England et al., 2008; Georgiadis et al., 2020; this study), seabird colonies (Davidson et al., 2018) and NOW dynamics (Jackson et al., 2021) suggest intermittent formation of a southern ice arch in Nares Strait, likely in response to regional atmospheric cooling spikes (Lecavalier et al., 2017). This interval spans the inception of a small late Holocene ice tongue in Petermann Fjord (Fig. 7h), suggesting that increased formation of landfast ice in Nares Strait may have reduced the oceanic heat flux to Petermann Fjord and stabilized the glacier front. Together with atmospheric cooling, this likely aided the formation of a small ice tongue at ∼ 2200 cal yr BP (Reilly et al., 2019). Similar mechanisms might have been at play during the interval after 1400 cal yr BP spanning the rapid growth of the Petermann ice tongue to its modern limits at ca. 600 cal yr BP (Reilly et al., 2019). Regional sea ice records indicate the onset of a more recurrent southern arch from ∼ 1400 cal yr BP, with a further decrease in sea ice and primary productivity markers in Petermann Fjord and a concurrent increase in IP₂₅ in Kane Basin at 950 cal yr BP implying the formation of more stable southern ice arch conditions and potentially even double ice arching. Thus, both the inception of the ice tongue in Petermann Fjord at ∼ 2200 cal yr BP and its rapid extension at 600 cal yr BP are preceded by a transition towards stable southern ice arch conditions in the Nares Strait by 300–350 years (Figs. 7, 6b, d). More severe sea ice, with longer periods of landfast ice in Nares Strait, thus likely promoted ice tongue growth in Petermann Fjord, either directly via stabilizing the calving front or indirectly via changes to the inflow in AW. In light of the recent changes in Nares Strait ice arch configuration (Vincent, 2019), our data suggest that the emerging dominance of the northern ice arch associated with year-round mobile sea ice in Nares Strait and a shorter landfast sea ice season in Petermann Fjord will likely contribute to the destabilization of the ice tongue. Thus, a future reduction in landfast sea ice in Nares Strait and adjacent fjords would likely contribute to enhanced mass loss from the GrIS.