Atlantic Water advection vs glacier dynamics in northern Spitsbergen since early deglaciation

Martin Bartels, Jürgen Titschack, Kirsten Fahl, Rüdiger Stein, Marit-Solveig Seidenkrantz, Claude Hillaire-Marcel, Dierk Hebbeln MARUM – Center for Marine Environmental Sciences, University of Bremen, 28359 Bremen, Germany 5 SaM – Senckenberg am Meer, Marine Research Department, 26382 Wilhelmshaven, Germany Alfred Wegener Institute for Polar and Marine Research, 27568 Bremerhaven, Germany Centre for Past Climate Studies and Arctic Research Centre, Department of Geoscience, Aarhus University, 8000 Aarhus C, Denmark GEOTOP – Université du Québec à Montréal, Montreal, H3C 3P8, Canada 10


Introduction
In the context of the ongoing global warming, the Arctic is gaining increasing attention from the scientific community and the general public (e.g.IPCC, 2014;Masson-Delmotte et al., 2013).This is primarily motivated by extreme events like the 2012 sea ice minimum, the declining summer sea ice cover of recent decades (Stroeve et al., 2012) and the significant contribution of melting Arctic glaciers to the rising global sea level (e.g.Alley et al., 2005;Gregory and Huybrechts, 2006;IPCC, 2014;Overpeck et al., 2006).The Arctic is extremely sensitive to climatic changes and due to various feedback mechanisms, it is warming twice as fast as the global mean -a phenomenon commonly called the Arctic amplification (Serreze and Barry, 2011).These drastic changes in the Arctic also likely have far-reaching impacts.Several studies draw a connection to climate phenomena (especially severe winters) at lower latitudes (e.g.Cohen, 2016; Published by Copernicus Publications on behalf of the European Geosciences Union. M Bartels et al.: Atlantic Water advection vs. glacier dynamics Cohen et al., 2014;Kug et al., 2015;Mori et al., 2014;Overland et al., 2015;Overland, 2015), although, causal linkages between these concurrent events are open to discussion (e.g.McCusker et al., 2016;Meleshko et al., 2016;Overland et al., 2016).
While various studies suggest that these recent climate changes are induced by human activities (e.g.Abram et al., 2016;IPCC, 2014), it still remains unclear to which extent natural mechanisms may be involved.Instrumental climate data rarely cover more than the last century; sea ice observations using satellites did not even start before late 1978 (Stroeve et al., 2012).Hence, data reaching far beyond these instrumental time series are needed to assess the full range of natural environmental variability, especially the linkage between oceanic current changes and glacier instability.Sedimentary records covering centuries and millennia offer such archives of natural climate variations during preindustrial times, i.e. without any significant anthropogenic impact.Also in the geological past, the Arctic has experienced periods of warming when ice shelves broke up and glaciers melted dramatically.The transition from the Pleistocene to the early Holocene, and especially the Holocene Thermal Maximum (HTM), when enhanced northward advection of Atlantic Water (AW; e.g.Sarnthein et al., 2003) and maximum insolation were recorded at high northern latitudes (Berger and Loutre, 1991), may help demonstrate ocean-glacier interactions under a fast warming trend, possibly comparable to the ongoing global warming.
Today, the advection of warm ocean currents is a main contributor to basal glacier melting in both hemispheres, which eventually leads to global sea level rise (e.g.Cook et al., 2016;Luckman et al., 2015;Straneo et al., 2010).Being situated at the only deep-sea gateway that connects the Arctic Ocean with the North Atlantic, the Svalbard Archipelago is an ideal location to investigate changes in the inflow of relatively warm AW into the Arctic Ocean and its influence on the sea ice regime, regional glacier activity and climatic conditions.This passage -the Fram Strait -is a major gateway for heat advection to high northern latitudes as well as for sea ice export from the Arctic towards lower latitudes.The majority of Arctic sea ice (∼ 90 %; Rudels, 2009) and about 50 % of its total freshwater (during the time period 1979-2001;Serreze et al., 2006) are exported through the Fram Strait.More recent studies report an increasing liquid freshwater export of almost 4 % but decreasing sea ice export of ∼ 17 % for the time period 2000-2010 compared to the previous 2 decades (e.g.Haine et al., 2015).Due to circulation of the northernmost branch of the North Atlantic Current along its western slope, Svalbard is bathed by warmer waters than any other Arctic location (Fig. 1a, b), providing the area with a relatively mild climate.In the present study, we chose a site at the mouth of Woodfjorden (Fig. 1c), located at the northern margin of Svalbard where the main branch of AW enters the Arctic Ocean, for the documentation of ocean-ice interactions from deglacial times to the present.In this fjord system, which is fed today by tidewater glaciers, enhanced sedimentation rates offer a high temporal resolution.
During the Late Glacial Maximum (∼ 21-18 ka; Svendsen et al., 2004), Svalbard's glaciers and ice caps were part of the Svalbard-Barents Sea Ice Sheet.This grounded ice sheet covered the entire Barents Sea and expanded as far as the shelf break of Svalbard (Landvik et al., 1998).It must have been connected to the Fennoscandian Ice Sheet during the Last Glacial Maximum (Hughes et al., 2016;Landvik et al., 1998;Svendsen et al., 2004).Woodfjorden and its trough constituted one of the main northern gateways for fast flowing ice streams that drained the northwestern sector of these ice sheets during the Last Glacial Maximum (Hormes et al., 2013;Ottesen et al., 2007).Whereas the Holocene and in part the late Weichselian have been relatively well documented off (south)western Svalbard (Hald et al., 2004;Jernas et al., 2013;Majewski et al., 2009;Rasmussen et al., 2007Rasmussen et al., , 2012;;T. L. Rasmussen et al., 2014a;Sarnthein et al., 2003;Skirbekk et al., 2010;Ślubowska-Woldengen et al., 2007;Werner et al., 2011), palaeoceanographic records from northern Svalbard (Chauhan et al., 2014(Chauhan et al., , 2015;;Koç et al., 2002;Kubischta et al., 2011;Ślubowska et al., 2005) are still scarce (Fig. 1b).The present study aims to fill this gap by compiling multi-proxy data, in particular, physical properties of the sediment, the faunal distribution and geochemical composition of benthic foraminifera, and biomarkers.Additionally, we generated a high-resolution record of ice-rafted debris (IRD) based on computer tomography.This set of proxies offers the opportunity to reconstruct the influence of AW in surface, intermediate and bottom waters, environmental conditions in surface waters and at the bottom water-sediment interface (including sea ice coverage and primary and export productivity), and glacier activity.Additionally, the study site is located at the transition zone from a glacier-influenced fjord milieu to an open ocean setting; therefore, very sensitive glacier responses to minor changes in oceanic forcing are expected.
The Woodfjorden system in northern Spitsbergen, the main island of the Svalbard Archipelago, includes the eponymous main fjord and the two tributary fjords Bockfjorden and Liefdefjorden (Fig. 1c).At present, only the latter is fed ) at station MSM02/03-666-3 (same location as the studied sediment core GeoB10817-4; see Sect. 3) derived during cruise MSM02/03 applying a Sea-Bird CTD sensor as described in the cruise report (Lherminier et al., 2009).The profile illustrates subsurface Atlantic Water inflow at the core site (grey shading), framed by less saline surface water on top and colder winter water at the bottom.Winter water forms locally by convection processes during autumn-winter (Cottier et al., 2010).Temperature (> 3 • C) and salinity boundaries (> 34.65) following Cottier et al. (2005).

Material and methods
The 275 cm long sediment core GeoB10817-4 used in this study was retrieved at the mouth of Woodfjorden (N Spitsbergen; 79.    ing halves were sampled with 10 mL syringes for microfossil analyses and grain size measurements at 5 cm intervals as well as for biomarker analyses at a lower resolution at selected depths.

Computer tomography
Archive halves of the sediment core GeoB10817-4 were scanned using a Toshiba Aquilion 64 ™ computer tomograph (CT) at the hospital Klinikum Bremen-Mitte, with an Xray source voltage of 120 kV and a current of 600 mA.
The CT image stacks have a resolution of 0.35 mm in the x and y directions and 0.5 mm resolution in the z direction (0.3 mm reconstruction unit).Images were reconstructed using Toshiba's patented helical cone beam reconstruction technique.The CT data obtained were processed using the ZIB edition of the Amira software (version 2016.25;Stalling et al., 2005).Within Amira, the CT scans of the core sections were merged and core liners, including about 2 mm of the core rims, were removed from the data set.All clasts > ∼ 1 mm and bioturbation traces were quantified in each reconstruction slice with the "Segmentation Editor" (threshold segmentation) and the "Material Statistics" module.The threshold values used were > 1500 for lithic clasts, 601-1499 for matrix sediment, 1-600 for bioturbation traces (open voids within the sediment that were filled by air and/or water were included), and < 1 for the surrounding air and water.With the "Connected Components" module, the individual clasts were separated and subsequently analysed with the "Shape Analysis" module for their parameterization.The determined clast length was further used for a clast size analysis.Therefore, every clast within an interval of 167 CT slices (corresponding to a ∼ 5 cm core interval) was considered and the obtained result was written to the depth of the central slice position.The analysing interval was moved slice by slice.The final results (unit: vol.% of all segmented clasts) were exported in a spreadsheet.Additionally, the X-ray density of the matrix sediment was evaluated by calculating the mean value and its standard deviation of the matrix sediment per slice (measured in Hounsfield units (HU); the matrix sediment segmentation was reduced by two voxels to exclude potential marginal artefacts).

AMS radiocarbon measurements
The chronostratigraphy of core GeoB10817-4 is based on 11 accelerator mass spectrometry (AMS) radiocarbon measurements from either mixed benthic foraminifera or N. labradorica tests, achieved at the CologneAMS facility of the University of Cologne or at the ETH of Zürich, respectively (Table 1).Radiocarbon ages were converted into calibrated years before present (before 1950 CE) using the radiocarbon calibration software CALIB 7.1 (calib.qub.ac.uk/calib/) and the Marine13 calibration curve (Reimer et al., 2013).Reservoir ages may have changed in the past, especially at high northern latitudes as well as during the Younger Dryas (YD) and during Heinrich Stadial 1 (e.g.Austin et al., 1995Austin et al., , 2011;;Bard et al., 1994;Hanslik et al., 2010;Waelbroeck et al., 2001).However, as no information about palaeo-reservoir ages is available for the Svalbard area (see also Ślubowska-Woldengen et al., 2007), we applied an averaged regional reservoir age correction ( R = 98 ± 37 years), based on six R values from the Svalbard region (Mangerud, 1972;Mangerud and Gulliksen, 1975;Olsson, 1980)  Sediment grains larger than 310 µm have not been detected by the laser diffraction particle size analyser, although coarser sand particles and even gravels are detected within the sieved samples of the microfossil analysis as well as within the CT clast size analyses (see Sects. 3.4 and 4.1;Fig. 3a).This bias results from the very small sample volume (∼ 0.5 mL) for the particle size analysis.However, taking into account that > 98 wt.% of the dried bulk sediment is < 150 µm (with the exception of four samples with 3, 4, 6 and 12 wt.% > 150 µm fraction; Bartels, unpublished data), the observed grain size patterns of the siliciclastic fraction are considered as reliable.

Microfossil analyses
Bulk sediment samples were washed through a 63 µm sieve.The coarse fraction was dried and sieved through 100 µm and 150 µm meshes.Faunal analyses were carried out on two fractions: 100-150 µm and > 150 µm (Fig. A1 in the Appendix).Counts from both fractions were summed up (i.e.> 100 µm) for the discussion below.Samples with high numbers of foraminifera were split with a microsplitter.Wherever possible, 200-300 benthic foraminifera were counted in each fraction of every sample to obtain statistically ro-bust results (Murray, 2006;Patterson and Fishbein, 1989).In addition, three replicates of an un-split sample as well as of a split sample were counted (each split was counted four times, including the initially counted sample).Statistical calculations reveal a standard error of 3.52 and 2.95 %, respectively.Most benthic foraminifera were identified to species level, except for agglutinants, which were mainly identified to genera level due to their poorer preservation.The majority of calcareous benthic foraminifera are well preserved, but several specimens show indications of dissolution or transport.Small fragments were not counted since these specimens might be allochthonous.Distinct Stainforthia and Buccella morphotypes were lumped together as Stainforthia loeblichi s.l. and Buccella frigida s.l., respectively, because in each case possibly two species (S. loeblichi s.s. and S. concava as well as B. frigida s.s. and B. tenerrima, respectively) could not be distinguished definitely.Hyaline and agglutinating species that could not be identified are summarized as Rotaliina and Textulariina, respectively.Only a few planktonic foraminifera were present in the material.As postmortem degradation might be the main reason for the distribution of arenaceous taxa (Murray, 2006), percentages of calcareous benthic species were calculated excluding agglutinants to avoid these taphonomic effects.
Accumulation rates of benthic foraminifera were calculated following the equation of Ehrmann and Thiede (1985): where ARBF is the accumulation rate of benthic foraminifera [individuals cm −2 ka −1 ], LSR the linear sedimentation rate [cm ka −1 ], ρ the density of dry sediment [g cm −3 ] and BFN the benthic foraminifera number (per gram dry sediment) [individuals g −1 ].

Stable isotope measurements
Stable carbon and oxygen isotopes were measured at ZMT (Leibniz Centre for Tropical Marine Ecology) in Bremen, using a Finnigan MAT 253 ™ gas isotope ratio mass spectrometer with a Kiel IV automated carbonate preparation device (standard deviations of house standard: 0.02 ‰ for δ 13 C and 0.06 ‰ for δ 18 O), and at MARUM (Bremen), using a Finnigan MAT 251 ™ mass spectrometer with a Kiel I carbonate preparation device (standard deviations of house standard: 0.02 ‰ for δ 13 C and 0.03 ‰ for δ 18 O).The measurements were performed on two to eight tests (test size: 150-500 µm) of the epibenthic foraminifera Cibicides lobatulus (Ahrens et al., 1997;Dubicka et al., 2015) and the endobenthic foraminifera Nonionellina labradorica (Ahrens et al., 1997;Alve and Bernhard, 1995;Loubere and Rayray, 2016) for bottom water and pore water signals, respectively.Stable oxygen isotope values were corrected for a global ice volume effect (Waelbroeck et al., 2002) and adjusted for vital effects: +0.64 ‰ for C. lobatulus (Shackleton, 1974) and −0.2 ‰ for N. labradorica (Duplessy et al., 2005).
M Bartels et al.: Atlantic Water advection vs. glacier dynamics

Biomarker analyses
Selected dried sediment samples were homogenized and used for geochemical analyses.For the determination of the total organic carbon (TOC) content using an ELTRA ™ analyser, 100 µg of sediment was used.Alkenone analyses, carried out to reconstruct sea surface temperatures (SSTs), were extracted using 2 g of sediment with an accelerated solvent extractor (DIONEX ™ , ASE 200; 100 • C, 1000 psi, 15 min, dichloromethane and methanol (99 : 1, v/v) as solvent).The separation of compounds was carried out using open column chromatography (SiO 2 ) using n-hexane, dichloromethane (1 : 1, v/v) and dichloromethane.The composition of alkenones was analysed using gas chromatography (Agilent ™ 7890).Individual alkenone (C 37:4 , C 37:3 , C 37:2 ) identification is based on retention time and the comparison with an external standard.
The alkenone unsaturation index (U K 37 ) as proxy for SST [ • C] was calculated following Brassell et al. (1986) due to the presence of the C 37:4 alkenone.For calculating SST, the equation by Müller et al. (1998) -based on a global core-top calibration (60 • N-60 • S) -was used (replacing U K 37 using U K 37 due to the presence of the C 37:4 alkenones).This equation (Eq.2) yielded more realistic results compared to other equations whose temperature reconstructions were much too high (e.g.Prahl and Wakeham, 1987).The standard error of the calibration is reported as ±1.5 • C. The instrument stability was continuously controlled using re-runs of an external alkenone standard (extracted from Emiliania huxleyi cultures with a known growth temperature) during the analytical sequences.The range of the total analytical error calculated using replicate analyses is less than 0.4 • C.
For the quantification of IP 25 , its molecular ion (m/z 350) in relation to the abundant fragment ion m/z 266 of the internal standard (7-HND) was used (SIM mode).The different responses of these ions were balanced using an external calibration (Fahl and Stein, 2012).Brassicasterol and dinosterol were quantified as trimethylsilyl ethers using the molecular ions m/z 470 and m/z 500, respectively, in relation to the molecular ion m/z 464 of cholesterol-D 6 .For more details about the quantification of IP 25 as well as the sterols, see Fahl and Stein (2012).Accumulation rates of sea ice and phytoplankton biomarkers (IP 25 and brassicasterol, respectively) were calculated following Eq.( 1) by replacing the BFN term accordingly.
For more semi-quantitative estimates of present and past sea ice coverage, Müller et al. (2011) combined the sea ice proxy IP 25 and phytoplankton biomarkers in a phytoplankton-IP 25 index, the so-called PIP 25 index: with c equal to the mean IP 25 concentration divided by the mean phytoplankton biomarker concentration.
As phytoplankton biomarker, brassicasterol was used (for further discussion of advantages and limitations of the PIP 25 approach see Belt and Müller, 2013;Smik et al., 2016;Stein et al., 2012;Xiao et al., 2015).Following the classification scheme of Müller et al. (2011), PIP 25 values between 0.3 and 0.5, between 0.5 and 0.7, and 0.7 point to a reduced sea ice cover, a seasonal sea ice cover including an ice edge situation and an extended to perennial sea ice cover, respectively.

Core description
The sediment of core GeoB10817-4 consists of some coarse sand and gravel, with scattered pebbles (Fig. 3a) embedded in a homogeneous silty clay matrix sediment.The upper ∼ 30 cm of the sediment core are heavily bioturbated, as illustrated by the CT scan.The interval between 31 and 54 cm is characterized by a basal erosional unconformity, weakly visible on the CT picture (Fig. 3b), and high amounts of (extra-and intra-)clasts (Fig. 3b), sand (including high quantities of foraminiferal tests) and bivalve shell fragments.Furthermore, the interval exhibits an increasing X-ray density with an enhanced standard deviation accompanied by a coarsening upward of the siliciclastic fine fraction (< 63 µm; observed in the grain size analyses; Fig. A2 in the Appendix).These sedimentological features point to an allochthonous origin and a deposition by mass wasting.Consequently, this interval is excluded from the subsequent palaeoceanographic reconstructions.Below this disturbance and as deep as 140 cm, a small number of clasts are found, whereas the lower part of the core contains high amounts of coarse debris (Fig. 3a).

Chronology
Assuming linear sedimentation rates between dated depths (ranging between 8 and 121 cm ka −1 ), calibrated ages have been interpolated (Fig. 3c; Table 1).The segment below the 31-54 cm interval (Sect.4.1) covers a time span from ∼ 15.7 to 7.8 cal.ka, whereas the one above ranges from ∼ 1.8 to 0.3 cal.ka (from here on, all ages are given as calibrated ages and data are presented exclusively against age).This indicates that the deposition of the 31-54 cm layer led to a ∼ 6000-year gap in the sequence, spanning a large part of the Holocene.The radiocarbon age at 257 cm depth has been excluded from the age model because it was obtained from a much smaller sample (∼ 0.4 mg carbonate) than the overlying 237 cm deep older sample (∼ 4.3 mg) that has been measured four times with consistent results.

Ice-rafted debris
Clasts > ∼ 1 mm detected using CT analyses (Sect.3.1) were considered as IRD.The interpretation of the IRD abundances is based on the volume percentages and number of clasts observed on CT scans, as the combination of these data depicts the best available IRD distribution (Figs.3a, 4c): only a consideration of volume percentages would overestimate large clasts while the number of clasts alone would overestimate the significance of small debris.For comparison, accumulation rates of counted clasts (> 500 µm; based on Ehrmann and Thiede, 1985) from sieved syringe samples are shown in the Appendix (Fig. A3) illustrating that the CTbased data offer a much more precise IRD record.CT data show that IRD is more abundant at the core bottom, containing high quantities prior to ∼ 14.5 ka (up to 10 vol.%; up to 2.3 clasts cm −3 ), but decrease afterwards (down to 0 % and 0.4 clasts cm −3 ).From ∼ 13.1 ka on, IRD amounts rise again until ∼ 11.5 ka (max.7 %; max.2.8 clasts cm −3 ).From then on, the volume percentages and number of clasts illustrate a low but continuous flux of ice-rafted material (mean: 0.2 %; 0.3 clasts cm −3 ), and some increase particularly from ∼ 8.5 ka on (up to 1.9 % and 0.8 clasts cm −3 ).Above the reworked layer (from ∼ 1.8 ka on), a further increase in the IRD content is identified, especially for larger clasts (up to 7 % and 0.8 clasts cm −3 ; Fig. 4c).Some peaks in volume percent reflect single large clasts (marked with a star in Fig. 4c).

Faunal composition
Planktic foraminifera were rare or absent.Their percentages range between 0 and 2 % of the total fauna apart from one sample dated ∼ 13 ka (∼ 5 %).Thus, only benthic species were taken into account for calculations and interpretation.In total, 83 taxa (61 calcareous species) were identified (7-56 taxa per sample, 7-39 calcareous taxa per sample; Appendix A).Calcareous species represent the major part of the benthic fauna (∼ 80-100 %); agglutinating species occur mainly in the upper part above the reworked layer (up to ∼ 20 %) with decreasing percentages down-core (mean: ∼ 2 %; Fig. A4a in the Appendix).

Biomarkers
Alkenone unsaturation index (U K 37 ) values range from 0.1 to 0.3.They yield SST values ranging from ∼ 1 to ∼ 7 • C. The highest temperatures were attained between ∼ 11 and ∼ 10 ka, whereas the lowest values are calculated prior to ∼ 11.9 ka (Fig. 7a).
The accumulation rates of brassicasterol range between ∼ 1.5 µg cm −2 ka −1 (∼ 14.2 ka) and ∼ 167.5 µg cm −2 ka −1 (∼ 10.9 ka) (Fig. 7d).The accumulation rate of dinosterol shows an almost parallel trend to that of brassicasterol but values are 1 magnitude lower (Fig. A5  Acc.rate benthic foram.In the Woodfjorden area, the Heinrich Stadial 1 (until ∼ 14.6 ka) was characterized by a very high IRD content (Fig. 8b) signifying enhanced calving of the glacier front.The increased percentages of fine siliciclastic sediments within this interval (mean: ∼ 4-5 µm; Fig. 8a) point to some contemporaneous intensified meltwater outflows (in this study, sources of siliciclastic sediments were interpreted as terrigenous).Thus, it is assumed that parts of the Svalbard-Barents Sea Ice Sheet disintegrated during Heinrich Stadial 1 while stable oxygen isotope records of the NGRIP ice core show that atmospheric temperatures remained cold (Fig. 8g; S. O. Rasmussen et al., 2014).Contrastingly, bottom waters might have been relatively warm and/or saline as indicated by the appearance of the benthic foraminiferal species C. neoteretis (Fig. 8d) as well as the high percentages (up to ∼ 80 %) of C. reniforme (Fig. 5h).The occurrence of C. neoteretis has been connected with a strong influence of AW in bottom waters because this species is assumed to respond to a limited temperature and salinity range as well as to an enhanced food supply (Jennings et al., 2004;Jennings and Helgadottir, 1994;Mackensen et al., 1985;Rytter et al., 2002;Seidenkrantz, 1995;Seidenkrantz et al., 2013;Steinsund, 1994).C. reniforme also prefers saline bottom waters (Hald and Korsun, 1997;Jennings et al., 2004) (Brassell et al., 1986;Müller et al., 1998).(b) Phytoplankton (brassicasterol)-IP 25 index (PIP 25 ); 0.3 to 0.5: reduced sea ice cover ( * ), 0.5 to 0.7: seasonal sea ice cover ( * * ), 0.7: extended to perennial sea ice cover ( * * * ) (Müller et al., 2011).(c) Accumulation rate of the sea ice biomarker IP 25 .(d) Accumulation rates of phytoplankton biomarker brassicasterol.Grey vertical shadings; see Fig. 4. Dark grey diamonds: calibrated radiocarbon ages with error ranges.therefore be linked to the intrusion of AW into bottom waters.The advection of AW as a subsurface or even bottom water mass during Heinrich Stadial 1 as well as during earlier Heinrich stadials has already been proposed for the Svalbard area by Rasmussen et al. (2007) and T. L. Rasmussen (2014b).Today, basal melting of tidewater glaciers in Spitsbergen from high advection rates of the relatively warm AW is evident (Luckman et al., 2015).Thus, during Heinrich Stadial 1, the influx of AW in the Woodfjorden area possibly contributed to the destabilization of glacier fronts.Apart from oceanic heat advection, topographic features of the relatively wide Woodfjorden system may have promoted the glacier retreat.Gump et al. (2017) reported a deglaciation of narrower fjords with higher surrounding mountains in southwestern Norway delayed more than 1000 years (see also Stokes et al., 2014).
to sediment plumes loaded by glacial flour from the erosion of Devonian Old Red sandstones.Source rocks, i.e. the Wood Bay Formation, crop out at the northwestern coast as well as at the southeastern coast of Woodfjorden (Reinsdyrflya and Andrée Land, respectively; Fig. 1c; Harland, 1998).Melting glaciers from the former lowland may have been a potential source of these reddish sediments.This interpretation is supported by 10 Be ages of erratic boulders, which suggest that Reinsdyrflya started to deglaciate during Heinrich Stadial 1 and became ice-free during the Bølling-Allerød (Gjermundsen et al., 2013;Hormes et al., 2013).Rasmussen and Thomsen (2013) reported on pink sediment layers derived from Devonian red beds in sediment cores from the northwestern slope of Svalbard.They interpreted those layers as deposited by meltwater plumes during interstadials.Forwick et al. (2010) also suggested that meltwater run-off connected to retreating glaciers was one possible source for fine reddish sediments deposited close to glacier fronts.Correspondingly, a glacier-proximal position is likely at our core site where tidewater glaciers may have successively waned.Percentages of C. neoteretis reached their maximal values during the Bølling (Fig. 8d).Almost simultaneously, this taxon peaked at the southwestern margin (Rasmussen et al., 2007) and at the northern margin of Svalbard ( Ślubowska-Woldengen et al., 2007;Ślubowska et al., 2005) as well as in the Laptev Sea (off northern Siberia; Taldenkova et al., 2010) -although with considerably higher values probably deriving from a more direct influence of AW and greater water depths compared to our core location.Massive meltwater outbursts connected to the retreating glaciers possibly caused a deepening of the halocline (cf.Rasmussen and Thomsen, 2004) as the extended cold and relatively fresh surface water mass forced the more saline and thus denser AW to submerge.Accordingly, cold surface waters have been insulated from those relatively warm bottom waters, enabling an extended sea ice coverage (see SST and sea ice cover; Fig. 8).
In the course of the Allerød (∼ 14-12.7 ka), the grain size distribution shows a general trend towards slightly coarser sediment (mean: ∼ 5.5-6 µm; Fig. 8a).The coarsening may reflect the retreat of some glaciers to more distal positions.Consequently, less fine-grained material from meltwater plumes would have reached the core location.Szczucin and Zaj ączkowski (2012) showed that meltwater-derived suspended particulate matter is mostly deposited within a few kilometres from sources in glacier-influenced fjords of Svalbard, as mixing of meltwater with sea water induces flocculation and a subsequent suspension settling.Ł ącka et al. ( 2015) observed similar phenomena during the deglaciation of the Storfjord Trough (southern Spitsbergen).A further retreat of the glaciers is also suggested by slightly increasing numbers of (most likely ice-rafted) clasts in the core studied (Fig. 8b), indicating enhanced calving.Radiocarbon-dated beach ridges at the eastern coast of Woodfjorden (Andrée Land) point to a glacio-isostatic emergence of ∼ 20 m (Forman et al., 2004) in the course of the Allerød (Fig. 8f).This uplift illustrates a massive retreat of the Svalbard-Barents Sea Ice Sheet, at least following Heinrich Stadial 1, as the decreasing ice sheet load caused a successive continental rebound.Despite this Allerød sea level drop, AW advection may have only been slightly reduced (if at all).The inflow of oceanic water into Woodfjorden is not restricted by a sill in a narrower sense, only on the seaward side; the peninsula Reinsdyrflya and the island Moffen are connected by a shallower ridge (∼ 80-100 m depth; see Fig. 1c), possibly slightly decelerating the AW influx.
Towards the transition to the YD, percentages of C. neoteretis increased again (Fig. 8d).This trend is also in agreement with studies from the western, southwestern and northern Svalbard margin (Rasmussen et al., 2007;Ślubowska-Woldengen et al., 2007;Ślubowska et al., 2005) as well as from East Greenland (Jennings et al., 2006), although percentages of C. neoteretis were significantly higher at these locations during the Allerød-YD transition.The concomitant appearance of C. neoteretis may indicate some AW influx reaching the sea floor (e.g.Lubinski et al., 2001;Steinsund, 1994).

Younger Dryas (∼ 12.7-11.7 ka) -oceanographic transitions and glacier responses
During the early YD, the benthic foraminiferal fauna was still dominated by C. reniforme (Fig. 5h) but percentages of N. labradorica rapidly increased (Fig. 8c).The latter species prefers fresh organic matter and is even able to sequester chloroplasts (Bernhard and Bowser, 1999;Cedhagen, 1991).N. labradorica shows maximum abundances between 1 and 2 cm sediment depth (e.g.Ahrens et al., 1997;Alve and Bernhard, 1995;Loubere and Rayray, 2016;Newton and Rowe, 1995).In general, this species is associated with glacier-distal habitats (Korsun andHald, 1998, 2000).It commonly blooms in connection with enhanced productivity found near oceanic fronts (e.g.Rytter et al., 2002;Sheldon et al., 2016b;Steinsund, 1994).Ślubowska-Woldengen et al. (2007) concluded that rising percentages of N. labradorica reflect the approaching Arctic Coastal Front, which divided Atlantic (i.e. the Svalbard Branch) and Arctic waters.Studies from the northern Svalbard continental slope (Koç et al., 2002;Ślubowska et al., 2005) and from the Bellsund trough (western Spitsbergen; Ślubowska-Woldengen et al., 2007) show similar distributions of N. labradorica during the YD.According to reconstructions of the Svalbard-Barents Sea Ice Sheet, the Barents Sea itself was no longer covered by this ice sheet at the beginning of the YD (Hormes et al., 2013), enabling the East Spitsbergen Current to surround the Svalbard Archipelago, transporting Arctic Water to the western and northern Svalbard margins and thus causing the development of the Arctic Coastal Front.At the study site in Woodfjorden, AW may have risen in the water column, possibly detaching from the sea floor and enabling the formation of denser winter water at the bottom (see modern situation illustrated in Fig. 2).The rise of AW may have been caused by a slightly declined meltwater outflow (Figs. 4a,8a), which no longer forced AW to submerge (see Sect. 5.1), and/or by upwelling events.Modern AW intrusion into the fjords of Spitsbergen occurs intermittently (e.g.Carmack and Wassmann, 2006;Wassmann, 2015) and is amplified by windinduced upwelling (Cottier et al., 2005(Cottier et al., , 2007)).Today, easterlies are the prevailing winds at the northern margin of Svalbard (Lind and Ingvaldsen, 2012).Then, northward Ekman transport of surface waters generates (southward) upwelling of subsurface (or intermediate) AW, which subsequently is able to flood the northern shelf of Svalbard (Falk-Petersen et al., 2015;Lind and Ingvaldsen, 2012).A northward retreat of the sea ice edge even amplifies this mechanism (Falk-Petersen et al., 2015, and references therein).Comparable conditions might have occurred very frequently at the mouth of Woodfjorden during the early YD, especially because the sea ice coverage was temporarily reduced during this time interval (Fig. 7b).Due to the resultant shallower halocline, an oceanic front (i.e. the Arctic Coastal Front) formed, separating AW from Arctic Water.In the Arctic, AW is commonly associated with enhanced nutrition: nutrients are transported from deeper water masses to the surface in the North Atlantic by winter mixing.The subsequent northward advection of these nutrients promotes primary production at higher northern latitudes while Arctic-derived waters contain only low nutrient amounts (Hunt et al., 2016).Thus, the advection of AW distributes nutrients onto the Arctic shelves, e.g. the shelves and fjords of Svalbard (Carmack and Wassmann, 2006).Therefore, the inflow of (sub)surface AW in the photic zone is interpreted as leading directly to a higher productivity (Carmack and Wassmann, 2006;Falk-Petersen et al., 2015;Sakshaug, 1997).Accordingly, nutrient advection and consequentially increasing phytoplankton blooms in spring and early summer may have resulted from such upwelling events during the early YD.High δ 13 C offsets ( δ 13 C; Fig. 8e) between the epi-(C.lobatulus) and the endobenthic species (N.labradorica) during this time interval point to enhanced remineralization rates and consequentially enhanced export production (see, e.g.Mackensen, 2008;Mackensen and Schmiedl, 2016).Therefore, δ 13 C values are considered to be valuable palaeoproductivity proxies (Mackensen et al., 2017).Apart from a generally amplified productivity resulting from the stronger summer insolation of the YD (Fig. 8g; Berger and Loutre, 1991), as suggested by Müller and Stein (2014), the increased (export) production possibly resulted from AW incursion into the photic zone.The subsurface advection of the relatively warm AW may also have contributed to further retreat of local glaciers: maximum amounts of IRD (Fig. 8b) and increased percentages of fine-grained sediments (mean: ∼ 6 µm; Fig. 8a) point to increased calving rates and meltwater outflow.
Following the peak appearance of N. labradorica, a shift in the faunal distribution is recognized around 12.4 ka: percentages of C. reniforme declined while those of E. clavatum rose (Fig. 5e, g, h).The dominance of E. clavatum possibly responded to comparable unstable conditions linked to sediment-laden waters and/or high-frequency IRD deposition connected to glacier melting events.E. clavatum is found in extreme modern environments, e.g.affected by very turbid meltwater plumes and/or relatively low salinities, where C. reniforme is not able to survive (Korsun andHald, 1998, 2000;Steinsund, 1994).Also at the coast of Newfoundland, eastern Canada, a comparable faunal change occurred during a similar time period (Pearce et al., 2014).A shift from relatively stable to unstable conditions during the YD has been inferred from various studies, although this shift may have started at the earliest at ∼ 12.3-12.2ka (Bakke et al., 2009;Lane et al., 2013;Pearce et al., 2013).The slightly earlier inception at Woodfjorden (and off Newfoundland) may result from dating uncertainties and/or regional differences.Lane et al. (2013) showed, based on lacustrine records from Germany and Norway, that timing differences are likely connected to diverging regional conditions.Interestingly, Ł ącka et al. ( 2015) describe a contrary trend to our record with warmer intervals during the late YD in the Storfjord Trough (southern Spitsbergen).Mild summers and a strong seasonality during the YD are also reported from southern and eastern Greenland (Björck et al., 2002;Hall et al., 2008).Cabedo-Sanz et al. (2013) report about ameliorating, sea-ice-free conditions off northern Norway as early as ∼ 11.9 ka.In our record, lower IRD deposition (Fig. 8b) along with gradually declining percentages of E. clavatum (Fig. 5g) are observed during the late YD (from ∼ 12 ka on), which probably indicates a further retreat of glacier fronts into inner fjord positions with a consequential reduction in icebergs melting at the mouth of Woodfjorden.Salvigsen and Høgvard (2006) noted that several glaciers at the head of Bockfjorden, a tributary fjord of the inner Wood-fjorden (Fig. 1c), were disconnected from the main fjord at the end of the YD.Prevailing fine-grained sediments (mean: ∼ 6 µm; Fig. 8a) -presumably connected to meltwater runoff -also point to some disintegration of glaciers.Birgel and Stein (2004) assumed that a "last deglaciation event" on Svalbard occurred around 12 ka, although they proposed a concurrent stable AW advection warming surface waters.In contrast, our data indicate that the spreading of cold and low-salinity surface waters -resulting from enhanced inputs of meltwater -likely deepened the Atlantic-derived waters.Accordingly, AW-influenced bottom waters may have resulted in the observed slightly increasing percentages of C. neoteretis (Fig. 8d) during the late YD (e.g.Jennings et al., 2004;Rytter et al., 2002;Steinsund, 1994).As surface waters were consequently better insulated from relatively warm AW, SSTs probably fell again while sea ice expanded (see Fig. 7a, b).However, a perennial closed sea ice cover in Woodfjorden, as suggested by Brückner and Schellmann (2003), seems unlikely as the data above show continuous deposition of IRD during the entire YD (Fig. 8b), even when maximum PIP 25 values are observed (Fig. 7b).This leads to the inference that broken-up sea ice at least during summer enabled high productivity of sea ice diatoms and iceberg calving and dispersal.Hence, it remains unclear if C. neoteretis responded rather to those food pulses during summer, as proposed by Wollenburg et al. (2004), than to the presence of AW.
The retreat of glaciers during the YD as suggested by our data is further supported by studies based on radiocarbondated beach ridges or moraines (Mangerud and Landvik, 2007;Salvigsen and Høgvard, 2006;Salvigsen and Österholm, 1982) as well as on 10 Be-dated boulders (Reusche et al., 2014).They point to the absence of any YD glacier advance along the western and northern coasts of Spitsbergen.Possible contemporaneous glacier advances in several other regions of Svalbard -especially in Isfjorden (central Spitsbergen) -still remain under debate (Birgel and Stein, 2004;Boulton, 1979;Forwick, 2005;Svendsen et al., 1996).Nonetheless, YD glacier retreats have also been observed in the Disko Bay, western Greenland, where the Jakobshavn Isbrae retreated from ∼ 12.2 ± 0.6 ka onwards (Rinterknecht et al., 2014).Rinterknecht et al. (2014) connect this retreat to warm ocean currents, which destabilized the marine terminating glacier by basal melting.A possibly ongoing deglaciation already at the onset of the YD and a further glacier retreat in the course of the YD is controversially discussed for the Scottish Highlands (Bromley et al., 2014;Small and Fabel, 2016).(∼ 11.7-7.8

ka) -amelioration of environmental conditions
Climatic conditions were ameliorated during the YD-early Holocene transition (starting ∼ 11.7 ka) as summer insolation approached its maximum values (Fig. 8g; Berger and Loutre, 1991).Also, the high offset between δ 18 O values of C. lobatulus and N. labradorica (Fig. 6a) -which already started during the second half of the YD and persists during the entire early Holocene -presumably indicates increasing seasonal differences due to the contrast between high summer and low winter insolation (Fig. 8g), as suggested by Rasmussen et al. (2012) for the Isfjord.Indeed, C. lobatulus calcifies in summer, whereas N. labradorica calcifies in spring or early summer (Zaj ączkowski et al., 2010).Rising δ 13 C and lower δ 13 C values in N. labradorica (Figs.6b, 8e) may indicate higher organic matter inputs deriving from enhanced primary production.High percentages of N. labradorica (Fig. 8c) point to a growing influence of the Arctic Coastal Front with intensified AW circulation in intermediate waters (see Sect. 5.2).Simultaneously, declining IRD amounts suggest a significant reduction in iceberg rafting (Fig. 8b).Grain size measurements illustrate some coarsening (Fig. 8a), likely resulting from the retreat of glaciers to innermost positions in the fjords as smaller amounts of fine-grained suspended material associated with meltwater plumes reached the core location (see Szczucin and Zaj ączkowski, 2012).Between ∼ 11.1 and ∼ 10.8 ka, a short-term shift in the faunal assemblage is recognized: percentages of N. labradorica dropped while those of E. clavatum and C. reniforme increased (Fig. 5e, g, h).This likely illustrates a northward shift of the Arctic Coastal Front (away from the coast) and a consequentially stronger influence of cold and fresh Arctic waters.Also, some surface cooling is indicated by rapidly decreasing SSTs (Fig. 7a).E. clavatum -an opportunistic species (e.g.Korsun and Hald, 2000;Newton and Rowe, 1995) -took advantage of these deteriorated conditions.Increasing percentages of C. reniforme may reflect the presence of some saline AW at the fjord floor as this species is reported to respond to Atlantic (Intermediate) Water, e.g. at the eastern Greenland shelf (Jennings and Helgadottir, 1994).Several studies describe simultaneous and similar faunal changes around Svalbard (Rasmussen et al., 2012;T. L. Rasmussen et al., 2014a;Skirbekk et al., 2010;Ślubowska-Woldengen et al., 2007;Ślubowska et al., 2005), generally associated with the Preboreal Oscillation (PBO) cooling interval (Björck et al., 1996).Some divergence in the timing of the PBO -especially at the Hinlopen trough ( Ślubowska et al., 2005) may be due to chronological uncertainties and possibly distinct reservoir ages, as suggested by T. L. Rasmussen et al. (2014a).Different trigger mechanisms have been suggested for the PBO inception.Meltwater discharge from the disintegrating Fennoscandian Ice Sheet into the Baltic Ice Lake and subsequent drainage into the Nordic Seas have been proposed as one possible cause of the PBO cooling (Björck et al., 1996(Björck et al., , 1997;;Hald and Hagen, 1998).In contrast, Fisher et al. (2002) assumed that a single meltwater flood event from Lake Agassiz into the Arctic Ocean resulted in an expansion and thickening of sea ice as well as in fresher North Atlantic surface waters.Despite this cooling, there is no evidence for any significant glacier re-advance in the Woodfjorden area as IRD deposition did not change much (Fig. 8b).Grain size measurements still illustrate high but steadily decreasing fractions of fine sediments (∼ 7-8.5 µm; Fig. 4a, d), which might reflect continuous meltwater discharge and thus further deglaciation.Nevertheless, increasing percentages of coarser sediments (∼ 20 µm) are also observed during the PBO (Fig. 4a,  d).These coarser sediments probably resulted from sea-icerafted deposition as sediments > 10 µm are interpreted as transported by sea ice (Hebbeln, 2000).During previous time intervals, enhanced percentages of this coarser fraction (∼ 20 µm) were masked because they were diluted by high amounts of meltwater-transported finer sediment (e.g. grain size distribution at ∼ 13.9 ka; Fig. 4d).In contrast, presumed sea-ice-rafted debris is better recognizable during the PBO when the grain size distribution is less influenced by suspended sediments derived from meltwater plumes (see Fig. 4d).Sea ice biomarkers exhibit a decreasing trend indicating a moderate sea ice coverage (Fig. 7b), but peak values of IP 25 (Fig. 7c) illustrate a high number of sea ice diatoms, which possibly characterize a sea ice edge position at the mouth of Woodfjorden.Thus, the deposition of sea-ice-rafted material at the core site seems realistic.Subsequently, maximum accumulation rates of the phytoplankton marker brassicasterol (Fig. 7d) point to algal blooms in spring-summer under open water conditions that persisted throughout the early Holocene (see PIP 25 values; Figs.7b and 8).
At the end of the PBO, N. labradorica rapidly recovered (reaching up to ∼ 50 %; Fig. 8c), while E. clavatum and C. reniforme declined in percentage (Fig. 5g, h).This faunal transition probably reflects the return of the Arctic Coastal Front in the vicinity of the core site, accompanied by rising AW intruding into intermediate water masses in the fjord (see Sect. 5.2).This hydrographic change was likely a response to the onset of the Holocene Thermal Maximum (HTM) in Svalbard, at ∼ 10.8 ka (Miller et al., 2010), when summer temperatures significantly exceeded modern temperatures due to the higher summer insolation (Fig. 8g; Berger and Loutre, 1991).Peaks in the ARBF, with more than 50 000 ind. cm −2 ka −1 , as well as enhanced accumulation rates of brassicasterol (Fig. d), indicate high food supplies.The observed increase in δ 13 C values (Fig. 8e) implies a high export productivity lasting until ∼ 10.4 ka.This may be linked to the presence of AW as a subsurface or even surface water mass, hence amplifying productivity in the photic zone.The appearance of the thermophilous mussel Mytilus edulis at the outer part of Woodfjorden (Reinsdyrflya) supports the assumption of a surface injection of AW at least from ∼ 10.6 ka on (9815 ± 80 14 C a; Salvigsen, 2002).Several studies around Spitsbergen report about a benthic faunal pattern comparable to the record above during the HTM (Groot et al., 2014;Rasmussen et al., 2012;Skirbekk et al., 2010;Ślubowska-Woldengen et al., 2007).Ślubowska-Woldengen et al. (2008) showed that N. labradorica reached maximum percentages at sites influenced by the WSC or the Svalbard Branch.At the mouth of Woodfjorden, percentages of I. helenae (Fig. 5d) rose at the onset of the HTM.This species has been connected to sea ice algae blooms in various studies (mostly lumped with I. norcrossi; e.g.Steinsund, 1994).Because sea ice coverage shows a decreasing trend (see PIP 25 ; Fig. 7b), whereas alkenone-based SSTs reached peak values of ∼ 7 • C (Fig. 7a), it seems more likely that I. helenae did not respond specifically to sea ice algae blooms but rather to an enhanced productivity (probably due to open waters), as already proposed by Polyak et al. (2002), Seidenkrantz (2013) and Wollenburg et al. (2004).
During the HTM, increasing SSTs (based on planktic foraminiferal assemblages) were reconstructed at the western Svalbard continental slope (Hald et al., 2004), south of Spitsbergen at the Storfjorden slope (T.L. Rasmussen et al., 2014a), and at the western Barents slope (Sarnthein et al., 2003).Lacustrine records from NW Spitsbergen indicate maximum temperatures based on alkenones until ∼ 10.5 ka (van der Bilt et al., 2016).Carbonara et al. (2016) and Rigual-Hernández et al. (2017), report about maximum concentrations of diatoms and coccoliths (mainly warm water species) and inferred a causally strengthened advection of AW during the HTM south of Spitsbergen.Thus, a significant warming in the Svalbard area occurred during this time interval, with temperatures exceeding modern temperatures (compare alkenone-based SST data in Fig. 7a with hydrographic data in Fig. 2 and in Rasmussen and Thomsen, 2014;T. L. Rasmussen et al., 2014a).
Iceberg production remained low during the HTM.From ∼ 10.1 ka on, IRD percentages were reduced to a minimum (Fig. 8b) -although slightly increasing at around 9 ka -indicating a retreat of most glaciers to innermost fjord positions while some large tidewater glaciers (e.g.Monacobreen; Fig. 1c) stayed connected to the fjord system during the early Holocene.Radiocarbon-dated beach ridges show a contemporary drop in relative sea level almost down to the present level (Fig. 8f; Forman et al., 2004) in response to the deglaciation of Spitsbergen with the consequential rebound persisting into the HTM.
Grain size distribution during the HTM exhibits maxima around 20 µm (Fig. 4a) that may indicate relatively high fluxes of sea-ice-rafted debris.In this scenario, some winter sea ice cover would have been present, although sea ice biomarker concentration is low during the interval (see PIP 25 minima in Fig. 7b).The PIP 25 index is based on sea ice diatoms and is thought to respond to spring-summer conditions (Müller et al., 2011).However, a coarsening of the sediment due to stronger bottom currents cannot be excluded.
After ∼ 10.1 ka another slight shift in the faunal distribution is recognized (Fig. 5): N. labradorica declined whereas C. reniforme rose again (Fig. 5e, h).This faunal change may be linked to a repeated northward movement of the Arctic Coastal Front away from the core site in concert with submerging AW.Following a significant peak in the ARBF at ∼ 10.1 ka, all proxies for productivity reveal a decreasing trend (i.e.δ 13 C, accumulation rates of brassicasterol and benthic foraminifera; Figs.5a, 7d, 8e) indicating a decreasing AW influence in the photic zone.Slightly cooler conditions have been consistently reported from the Storfjord Trough (southern Spitsbergen) at the end of the early Holocene, probably due to a stronger influence of Arctic Water (Ł ącka et al., 2015).
Around 8.8 ka, a further coarsening of the sediment is observed (maximum percentages ∼ 45-50 µm, mean: ∼ 15 µm; Figs.4a, 8a).This shift in the grain size distribution may indicate higher current velocities.The pronounced skewness towards coarse sediments seen in the grain size distribution points to winnowing and/or bypass of the fine fraction (see grain size distribution at ∼ 8.8 ka in Fig. 4d).The epibenthic species C. lobatulus -commonly connected to strong bottom currents (Hald and Korsun, 1997;Hansen and Knudsen, 1995;Steinsund, 1994) -displays slightly increased proportions during this time period (Fig. 5b).Higher percentages of B. frigida s.l.(Fig. 5c) support these findings because the species develops preferentially in coarse sediments (Steinsund, 1994).This taxon has also been connected to ice edge algal blooms (Hald and Steinsund, 1996;Steinsund, 1994).
Here, however, this seems unlikely; analogous to I. helenae, its presence rather corresponds to reduced sea ice cover conditions during spring-summer, as inferred from declining PIP 25 values (Figs. 5c,7b).Seidenkrantz (2013) suggested that B. frigida thrives right after the retreat of sea ice.Apart from the species' preference for coarse sediments, the occurrence of B. frigida s.l.seems to be mainly controlled by the presence of open waters enabling primary production and high organic matter fluxes to the sea floor.Consistently, B. frigida is reported to feed on fresh particulate organic matter in the Northeast Water Polynya (northeastern Greenland; Newton and Rowe, 1995).An enhanced productivity at the mouth of Woodfjorden is indeed indicated by increasing δ 13 C values (Fig. 8e) simultaneously matching peak percentages of B. frigida s.l. at around 8.8 ka (Fig. 5c).

Late Holocene (∼ 1.8-0.4 ka) -aftermath of the neoglaciation
In comparison with the early Holocene, distinct differences in the benthic foraminiferal fauna mark the last ∼ 1,400 years of our record.The faunal composition suggests frequent AW incursions, temporarily as bottom water mass: the N. labradorica abundance decreased to ∼ 13-22 % (Fig. 8c), whereas C. reniforme exhibits percentages up to ∼ 31 % (Fig. 5h).C. reniforme may indicate more saline conditions in bottom waters linked to the presence of AW at the sea floor (Hald and Korsun, 1997;Jennings et al., 2004;Jennings and Helgadottir, 1994).Enhanced percentages of the opportunistic species E. clavatum (up to ∼ 28 % Fig. 5g) point to unstable conditions, e.g.high inter-annual variability in environmental conditions.Nonetheless, I. helenae reached pro-portions comparable to those of the HTM (> 5 %; Fig. 5d), reflecting still sufficient food availability.Increasing δ 13 C values (Fig. 8e) characterize this interval, likely responding to the temporary influence of AW in subsurface waters transporting nutrients into the photic zone.Late Holocene seasonally ice-free waters (see Fig. 7b) enabled primary productivity despite lower SSTs than those inferred above for the HTM interval (∼ 4 • C; Fig. 7a), very likely due to the decreasing summer insolation and reduced seasonal temperature contrast (Berger and Loutre, 1991;Fig. 8g).Cooling in the Fram Strait area during the late Holocene is reported at various locations and assigned reduced insolation (e.g.Hald et al., 2004;Müller et al., 2012;Rasmussen et al., 2012;Rigual-Hernández et al., 2017;Ślubowska-Woldengen et al., 2007;Werner et al., 2013).Although sea ice biomarkers suggest open waters during summer (Figs. 7b,8), increased proportions of ∼ 20 and ∼ 50 µm grain size fractions (Fig. 4a) may imply the deposition of (winter) sea-ice-rafted debris (see Sect. 5.3).Parallel to this cooling, a slight increase in IRD percentages indicates enhanced glacier activity with intensified iceberg production (Fig. 8b) in the Woodfjorden area.Furrer et al. (1991) report about advances of Monacobreen into Liefdefjorden (Fig. 1c) during the Older Subatlantic (∼ 1.2 ka; 1315 ± 100 14 C a), whose glacier ice front reached Lernerøyane, a small archipelago ∼ 6 km from its modern position.Advancing glaciers have also been reported in central Spitsbergen during a similar time period (Baeten et al., 2010;Humlum et al., 2005).The growth of Monacobreen possibly responded to the neoglacial cooling trend (as described by Porter and Denton, 1967), which led to waxing of Arctic glaciers (e.g.Miller et al., 2010;Solomina et al., 2015Solomina et al., , 2016)), and is reflected in increased IRD deposition in the present study.
In contrast to conditions inferred above, some other studies state a general warming trend (Bonnet et al., 2010;D'Andrea et al., 2012;Jernas et al., 2013;Majewski et al., 2009;Spielhagen et al., 2011;Werner et al., 2011) as well as a strengthened influx of AW during the last 2 millennia in the Svalbard area (Groot et al., 2014;Rasmussen et al., 2012;Sarnthein et al., 2003;Ślubowska et al., 2005).However, as the mid-Holocene is missing in the present record (see Sect. 4.1), we can only compare the early vs. late Holocene.As the late Holocene is marked by cooler conditions than those prevailing during the HTM, a link to the neoglaciation is assumed, but the actual trend from the mid-to the late Holocene is unknown.
Interpreting the appearance of C. neoteretis, numerous studies assume that the intensification of AW inflow along the Svalbard margin occurred relatively simultaneously during the Bølling-Allerød (Klitgaard Kristensen et al., 2013;Rasmussen et al., 2007;Ślubowska et al., 2005).In the Laptev Sea, this species has been connected to Atlanticderived water from ∼ 15.4 ka onwards (Taldenkova et al., 2010).AW presence is also inferred prior to ∼ 14 ka off West Greenland (Sheldon et al., 2016a).At the margin of the Arctic Ocean, north of Svalbard, C. neoteretis showed high abundances from at least 24 ka onwards (Wollenburg et al., 2004).Despite the assumed intrusion of relatively warm water masses, planktic foraminifera assemblages indicate low subsurface temperatures without considerable changes far into the YD along the W Barents slope (Sarnthein et al., 2003) and along the W Svalbard margin (Hald et al., 2004).Thus, a subsurface or bottom water inflow of relatively warm AW is assumed during the Bølling-Allerød, while surface waters remained cold.
The YD cold event has been studied extensively.This cold reversal was most likely initiated by a massive glacial lake drainage into the Arctic Ocean that reduced the Atlantic Meridional Overturning Circulation (Broecker, 2006;Broecker et al., 1989;Fahl and Stein, 2012;Murton et al., 2010;Stein et al., 2012;Tarasov and Peltier, 2005).Nonetheless, the influence of AW is seen at numerous locations around Svalbard, at least during the early YD and possibly temporarily (e.g.Rasmussen et al., 2007;Ślubowska et al., 2005; this study, see Fig. 8c, d).A shift of the Arctic Coastal Front towards the study areas at the northern Svalbard margin characterizes the situation during the YD.The proximal position of this oceanic front is linked to an initial AW intrusion into intermediate waters of Woodfjorden providing higher food supplies (Fig. 8e) while contributing to the retreat of tidewater glaciers (Fig. 8a, b).
During the early Holocene, the Nordic Seas, the Fram Strait and Barents Sea were characterized by a high AW inflow (e.g.Bauch et al., 2001;Carbonara et al., 2016;Groot et al., 2014;Müller and Stein, 2014;Rigual-Hernández et al., 2017;Risebrobakken et al., 2011;Telesiński et al., 2015;Werner et al., 2013Werner et al., , 2016)), coinciding with a general warming trend in the entire Svalbard region -especially after the PBO (e.g.Hald et al., 2004;Rasmussen et al., 2012;Sarnthein et al., 2003;Skirbekk et al., 2010;Ślubowska et al., 2005; this study, see Fig. 7a).This trend is in agreement with the early Holocene warming observed in the entire North Atlantic (Marchal et al., 2002).The interruption of the warming trend during the PBO occurred almost simultaneously around Svalbard but is more pronounced in northern locations (T.L. Rasmussen et al., 2014a;Skirbekk et al., 2010;Ślubowska et al., 2005; this study, see Figs. 5, 7) than along the (south)western margin of Svalbard (Rasmussen et al., 2012;T. L. Rasmussen et al., 2014a;Ślubowska-Woldengen et al., 2007).Accordingly, it seems likely that the PBO was triggered by a vast meltwater flood from the Mackenzie River area into the Arctic Ocean (Fisher et al., 2002) than by a meltwater discharge from the Fennoscandian Ice Sheet into the Nordic Seas as assumed earlier (Björck et al., 1996).Thus, northern sites may have been more strongly affected by the resulting colder and fresher surface waters and the expansion of a thicker pack ice, which concurred to force AW to sink deeper in the water column.
Nevertheless, as discussed above, AW was present at the northern margin of Spitsbergen during the last 2 millennia, although its influx was weaker and its position possibly deeper in the water column when compared to HTM conditions (Fig. 8c).Various studies report similar oceanographic conditions during this time interval (e.g.Hald et al., 2004;T. L. Rasmussen et al., 2014a;Skirbekk et al., 2010;Ślubowska-Woldengen et al., 2007;Ślubowska et al., 2005).Contrary to the HTM and the YD, there are indications of glacier advances regardless of the AW advec-tion during the late Holocene (Fig. 8a, b), possibly resulting from the lower summer insolation of the interval (Berger and Loutre, 1991; see Fig. 8g).Additionally, the AW incursion into Woodfjorden as a relatively deep water mass may have been of minor influence on glacier fronts as tidewater glaciers may only have been present in the shallower inner parts of the fjord system during the late Holocene (Furrer et al., 1991;Salvigsen and Høgvard, 2006).Hence, the glacier dynamics were possibly decoupled from AW advection during the late Holocene.

Conclusions
Multi-proxy analyses of sediment core GeoB10817-4 from the mouth of Woodfjorden (northern Spitsbergen) enabled us to document ocean-ice interactions since the deglaciation (since ∼ 15 500 years) at the northern margin of Svalbard, providing insights into glacier activity as well as into bottom, intermediate and surface water conditions, including sea ice coverage.
-The deglaciation period (until ∼ 12.7 ka; excluding the Younger Dryas) was marked by a disintegrating Svalbard-Barents Sea Ice Sheet, as illustrated by high calving rates and enhanced meltwater inputs.Cold surface waters were covered by extensive sea ice, while early intrusions of Atlantic Water as a bottom water mass occurred (Fig. 9a).
-During the early Younger Dryas (∼ 12.7-12.4ka), Atlantic Water was still advected as an intermediate water mass, probably connected to a coastward shift of the Arctic Coastal Front.Its intrusion was accompanied by warmer surface waters and high iceberg melting rates, which may have been linked to the retreat of some tidewater glaciers towards the fjord heads (Fig. 9b).During the late Younger Dryas (∼ 12.4-11.7 ka), surface waters cooled while sea ice extended and Atlantic Water temporarily penetrated into bottom waters.
-Conditions improved during the early Holocene (∼ 11.7-7.8ka), with an increasing subsurface Atlantic Water inflow, and were briefly interrupted during the Preboreal Oscillation.This deterioration occurred almost simultaneously around Svalbard, but was more pronounced at the northern margin of Svalbard than in the south.During the succeeding Holocene Thermal Maximum, sea surface temperatures reached their maximum, surpassing modern temperatures.Productivity flourished and sea ice coverage declined.Most glaciers retreated to innermost fjord positions during this interval.Only single tidewater glaciers survived (especially Monacobreen) (Fig. 9c).
-The early vs. late Holocene depict distinct environmental conditions.Productivity gradually increased in the course of the late Holocene (∼ 1.8-0.4ka) while Atlantic Water temporarily sank from intermediate to bottom waters.Even though sea surface temperatures were lower than during the early Holocene, sea ice coverage was reduced.High iceberg production possibly responded to some late Holocene glacier advance (neoglaciation) in response to low summer insolation (Fig. 9d).
-From the deglaciation on, the presence of Atlantic Water -in bottom as well as in intermediate waters -is connected to retreating glaciers at the northern Svalbard margin, while during the late Holocene, glaciers might have grown despite relatively constant inflow of Atlantic Water (possibly below surface).
-The Holocene Thermal Maximum includes -apart from a few differences (e.g.considerable higher summer insolation, dropping relative sea level due to the glacioisostatic rebound) -numerous similarities to a future warming: sea surface temperatures were most likely warmer than today, resulting in a reduced sea ice coverage and consequential open waters that enabled high primary production.Already today, reduced sea ice promotes increasing phytoplankton blooms (Arrigo andvan Dijken, 2011, 2015).During the Holocene Thermal Maximum, subsurface intrusion of warm Atlantic Water contributed to the melting and retreat of glaciers to inner fjord positions, a phenomenon that is also observed today, e.g. at Svalbard's tidewater glaciers (Luckman et al., 2015).Thus, environmental conditions during the Holocene Thermal Maximum are comparable with the already ongoing global warming, at least for the Svalbard region.

Figure 2 .
Figure2.Temperature and salinity against depth (in metres below sea level, m b.s.l.) at station MSM02/03-666-3 (same location as the studied sediment core GeoB10817-4; see Sect. 3) derived during cruise MSM02/03 applying a Sea-Bird CTD sensor as described in the cruise report(Lherminier et al., 2009).The profile illustrates subsurface Atlantic Water inflow at the core site (grey shading), framed by less saline surface water on top and colder winter water at the bottom.Winter water forms locally by convection processes during autumn-winter(Cottier et al., 2010).Temperature (> 3 • C) and salinity boundaries (> 34.65) followingCottier et al. (2005).

Figure 3 .
Figure 3. (a) CT scan of the core GeoB10817-4.Left to right: orthogonal profile; interpreted image (clasts and bioturbation; see legend); clast size distribution (0-20 vol.% of clasts: blue to red, respectively).The yellow rectangle marks the disturbed section also shown in (b).(b) CT images (orthogonal and interpreted) of the disturbed section (yellow rectangle; definition of boundaries; see Sect.4.1 and Fig. A2 in the Appendix).(c) Age-depth plot for sediment core GeoB10817-4: the diamonds show calibrated (median) ages with error bars (Table1); the solid black line shows linear interpolation excluding an outlier at 257 cm.The blue solid line shows corresponding sedimentation rates.

Figure 8 .
Figure 8.The blueish horizontal bar on top symbolizes sea ice coverage based on biomarker data (PIP 25 index).Red-blue horizontal bar symbolizes (qualitative) sea surface temperatures (red: warm, blue: cold) based on biomarker data (alkenones).(a) Mean grain size of siliciclastic (terrigenous) sediments: low values illustrate deposition of suspended sediment from meltwater plumes (note descending y axis).(b) Volume percentage of ice-rafted debris (solid red line) and number of clasts per cubic centimetre (reddish shading) signifying glacier activity (iceberg calving).(c) Relative percentages of N. labradorica.High percentages: Atlantic Water intrusion into intermediate waters delivering nutrients into the photic zone.Note the almost parallel increase in export production as illustrated in (e).(d) Relative percentages of C. neoteretis.High percentages: Atlantic Water advection into bottom waters.(e) Offset of epi-and endobenthic δ 13 C values indicating export productivity.(f) Radiocarbon-dated beach ridges in the Woodfjorden area(Forman et al., 2004) illustrating glacio-isostatic emergence (in metres above high tide, m a.h.t.; note descending y axis).Colours mark different locations as indicated in the legend (locations shown in Fig.1c).(g) Stable oxygen isotopes from NGRIP ice core (black line; S. O.Rasmussen et al., 2014)  and June-July insolation at 80 • N (orange line;Berger and Loutre, 1991).Grey vertical shadings; see Fig.4.Dark grey diamonds: calibrated radiocarbon ages with error ranges.

Figure 9 .
Figure 9. Environmental conditions in Woodfjorden during (a) the deglaciation, (b) the Younger Dryas, (c) the early Holocene and (d) the late Holocene.Black bars mark the location of sediment core GeoB10817-4.Glacier retreats or advances are marked by turquoise arrows.Yellow circles symbolize insolation changes(Berger and Loutre, 1991).The crossed circle in (b) signifies easterlies (blowing perpendicular to the profile) triggering a northward Ekman transport of the surface waters, which entailed an upwelling of intermediate waters (i.e.AW; seeLind and Ingvaldsen, 2012).

Figure
Figure A2.(a) Grain size distribution of fine siliciclastic sediments (0-63 µm).(b) Interpreted CT analysis showing bioturbation (green) and clasts (orange).(c) X-ray density (black line) in Hounsfield units (HU) with standard deviation (red line) derived from CT analysis.Yellow rectangle marks disturbed section.Note coarsening upwards of fine sediments (a), high accumulation of clasts (b) and increased X-ray density, as well as respective standard deviation (c) in the disturbed section.

Figure A3 .Figure A4 .Figure A5 .
Figure A3.Comparison of methods to determine the input of ice-rafted debris (IRD): accumulation rates(Ehrmann and Thiede, 1985) of clasts > 500 µm based on count data as well as volume percentages and numbers of clasts (> 1 mm) cm −3 based on CT analyses.

Table 1 )
; the solid black line shows linear interpolation excluding an outlier at 257 cm.The blue solid line shows corresponding sedimentation rates.

Table 1 .
(Reimer et al., 2013)rements and calibrated ages applying the Marine13 calibration curve(Reimer et al., 2013)and an averaged regional R = 98 ± 37 years (calib.qub.ac.uk).All measurements, apart from the measurements at 2.5 and 270.5 cm (CologneAMS, University of Cologne), were carried out at ETH Zürich.Median ages (in bold) were used for the age-depth model.