Articles | Volume 16, issue 4
Clim. Past, 16, 1127–1143, 2020
Clim. Past, 16, 1127–1143, 2020

Research article 03 Jul 2020

Research article | 03 Jul 2020

Surface and subsurface Labrador Shelf water mass conditions during the last 6000 years

Surface and subsurface Labrador Shelf water mass conditions during the last 6000 years
Annalena A. Lochte1,2, Ralph Schneider1, Markus Kienast3, Janne Repschläger4, Thomas Blanz1, Dieter Garbe-Schönberg1, and Nils Andersen5 Annalena A. Lochte et al.
  • 1Institute of Geoscience, Kiel University, Ludewig-Meyn Str. 10, 24118 Kiel, Germany
  • 2HOSST Graduate School, GEOMAR Helmholtz-Centre for Ocean Research Kiel, Wischofstraße 1–3, 24148 Kiel, Germany
  • 3Department of Oceanography, Dalhousie University, 1355 Oxford Street, Halifax, Canada
  • 4Department of Climate Geochemistry, Max Planck Institute for Chemistry, Hahn Meitner Weg 1, 55128 Mainz, Germany
  • 5Leibniz Laboratory for Radiometric Dating and Stable Isotope Research, Kiel University, Max-Eyth-Str. 11–13, 24118 Kiel, Germany

Correspondence: Annalena A. Lochte (


The Labrador Sea is important for the modern global thermohaline circulation system through the formation of intermediate Labrador Sea Water (LSW) that has been hypothesized to stabilize the modern mode of North Atlantic deep-water circulation. The rate of LSW formation is controlled by the amount of winter heat loss to the atmosphere, the expanse of freshwater in the convection region and the inflow of saline waters from the Atlantic. The Labrador Sea, today, receives freshwater through the East and West Greenland currents (EGC, WGC) and the Labrador Current (LC). Several studies have suggested the WGC to be the main supplier of freshwater to the Labrador Sea, but the role of the southward flowing LC in Labrador Sea convection is still debated. At the same time, many paleoceanographic reconstructions from the Labrador Shelf focussed on late deglacial to early Holocene meltwater run-off from the Laurentide Ice Sheet (LIS), whereas little information exists about LC variability since the final melting of the LIS about 7000 years ago. In order to enable better assessment of the role of the LC in deep-water formation and its importance for Holocene climate variability in Atlantic Canada, this study presents high-resolution middle to late Holocene records of sea surface and bottom water temperatures, freshening, and sea ice cover on the Labrador Shelf during the last 6000 years. Our records reveal that the LC underwent three major oceanographic phases from the mid- to late Holocene. From 6.2 to 5.6 ka, the LC experienced a cold episode that was followed by warmer conditions between 5.6 and 2.1 ka, possibly associated with the late Holocene thermal maximum. While surface waters on the Labrador Shelf cooled gradually after 3 ka in response to the neoglaciation, Labrador Shelf subsurface or bottom waters show a shift to warmer temperatures after 2.1 ka. Although such an inverse stratification by cooling of surface and warming of subsurface waters on the Labrador Shelf would suggest a diminished convection during the last 2 millennia compared to the mid-Holocene, it remains difficult to assess whether hydrographic conditions in the LC have had a significant impact on Labrador Sea deep-water formation.

1 Introduction

Since the early Holocene, the Labrador Sea has been an important region of the global thermohaline circulation system through the formation of intermediate Labrador Sea Water (LSW) in the central basin (Clarke and Gascard, 1983; Myers, 2005; Rhein et al., 2011; Yashayaev, 2007), which forms the upper part of North Atlantic Deep Water (NADW) and thereby contributes to the strength of the Atlantic Meridional Overturning Circulation (AMOC). Modern-type winter convection in the Labrador Sea (e.g. Hillaire-Marcel et al., 2001; Hoogakker et al., 2015) was likely established through a combination of diminished meltwater run-off following the final deglaciation of the Laurentide Ice Sheet (LIS) at about 7 ka (e.g. Jennings et al., 2015; Ullman et al., 2016) and the early Holocene strengthening of the West Greenland Current (WGC) bringing a larger proportion of warm, more saline Atlantic waters into the Labrador Sea (e.g. Lloyd et al., 2005; Lochte et al., 2019a; Seidenkrantz et al., 2013; Sheldon et al., 2016). In turn, high surface water buoyancy may have been a limiting factor in LSW formation during the early Holocene (e.g. Hillaire-Marcel et al., 2001; Hoogakker et al., 2015). Today, the rate of Labrador Sea convection is controlled by two main processes: heat loss to the atmosphere and the relative supply of buoyant freshwater to the convection region. While changes in the heat loss are predominantly related to atmospheric forcing, Labrador Sea salinity is controlled mainly by two freshwater supply routes (Wang et al., 2018): the eastern route via the East Greenland Current (EGC), which mixes with the Irminger Current (IC) south of Greenland to form the WGC, and the western route via the Labrador Current (LC), the extension of the Baffin Current (BC; Fig. 1a). While early studies suggested the LC as the main source of central Labrador Sea freshwater (Lazier, 1973, 1988; Khatiwala et al., 1999), more recent work advocated the WGC as the main supplier of freshwater (Cuny et al., 2002; Straneo, 2006; Schmidt and Send, 2007).

Figure 1(a) Map of the North Atlantic showing modern surface circulation with cold (blue) and warm (red) currents (adopted from Lazier and Wright, 1993; Rhein et al., 2011). A shaded area indicates the convection region of LSW. Core MSM45-31-1 is marked by a white dot; other cores referred to are marked by black dots. NF: Newfoundland. (b) Depth profiles of temperature (blue) and salinity (red) of the water column near the core site at 5428.53 N, 5604.33 W obtained from CTD deployment during cruise MSM45 in August 2015.

In order to assess the relative roles of the WGC and LC in regulating the freshwater budget and convection of the Labrador Sea, paleoceanographic reconstructions of these two main currents are critical. While several studies have provided middle to late Holocene reconstructions of the northward flowing WGC along the western coast of Greenland (Krawczyk et al., 2010; Lloyd et al., 2007; Moros et al., 2016; Møller et al., 2006; Seidenkrantz et al., 2007; Sha et al., 2017) or the IC (Moffa-Sánchez et al., 2014; Moffa-Sánchez and Hall, 2017), there is less evidence for water mass conditions in the southward flowing LC, as paleoceanographic reconstructions from the Labrador Shelf have mainly focussed on late deglacial to early Holocene meltwater run-off from the LIS (e.g. Jennings et al., 2015; Hoffman et al., 2012; Lewis et al., 2012; Hillaire-Marcel et al., 2007) and its impact on LC strength (Rashid et al., 2017). Conversely, most reconstructions of sea surface temperature and sea ice cover focus on the region around Newfoundland (Keigwin et al., 2005; Solignac et al., 2011; Sicre et al., 2014; Sheldon et al., 2015, 2016) and Orphan Knoll (Hoogakker et al., 2011, 2015), whereas very little information is available about middle to late Holocene water mass conditions in the central LC, as high-resolution late Holocene sediment records from the Labrador Shelf and slope are sparse.

Here, we provide detailed reconstructions of LC temperature and salinity changes since the middle Holocene, in order to discuss the current's potential role in central Labrador Sea convection. Paleoceanographic reconstructions are based on a high-resolution sediment record from the southern Labrador Shelf, providing information about sea surface and bottom water temperatures, freshening and sea ice cover of Labrador Shelf waters during the last 6000 years.

2 Oceanographic setting

The Labrador Sea is the coldest and freshest basin of the North Atlantic (Yashayaev, 2007) and is one of the major areas of open-ocean convection (Marshall and Schott, 1999). It plays a crucial role in influencing the strength and variability of the AMOC through the formation of LSW, which builds the upper part of NADW (Böning et al., 2006; Schmidt and Send, 2007; Schmitz and McCartney, 1993; Yashayaev, 2007; Yashayaev and Loder, 2009, 2016). Intermediate LSW forms through convection following an intense cooling of central Labrador Sea surface waters during particularly cold winters in North America (Clarke and Gascard, 1983). The annual variability of LSW formation was found to correlate with the North Atlantic Oscillation (NAO), the dominant mode of atmospheric variability in the North Atlantic region (Hurrell, 1995). The mode of the NAO wields a major influence on temperature and precipitation patterns across the North Atlantic region and is defined by the difference in atmospheric pressure at sea level between the Icelandic low and the Azores high-pressure systems (Hurrell, 1995). Positive NAO phases (NAO+) are generally associated with strong north-westerly winds over the Labrador Sea. These strong winds bring Arctic air southward, promoting enhanced winter cooling in the Labrador Sea region and thus stimulating deep convection (Dickson et al., 1996) and fostering the southward advection of sea ice (Dickson et al., 2000; Deser et al., 2000; Drinkwater, 1996).

The Labrador Sea receives freshwater through its bordering currents: the WGC and the LC (Wang et al., 2018). As an extension of the BC, the LC flows south-eastwards along the continental margin of Labrador and Newfoundland and exports buoyant freshwater from the Arctic Ocean to the subtropical North Atlantic (Fig. 1a). The LC is divided into two branches: a shallow inner branch that flows over the Labrador Shelf and an outer branch that is centred along the upper continental slope (Lazier and Wright, 1993). The inner LC is the immediate continuation of the BC (Cuny et al., 2002), augmented by Arctic water outflow from Hudson Bay (Straneo and Saucier, 2008), while the outer LC receives a significant proportion of Atlantic waters by the westward retroflection of the WGC, which itself is formed from the cold EGC, and the warm, saline IC (Cuny et al., 2002). Mixing of the inner and outer LC branches occurs in the saddles between the prominent banks on the Labrador Shelf (Vilks, 1980; see Figs. 1b, 2). Just south of Cartwright Saddle, in the vicinity of Flemish Cap and Newfoundland, the LC flows into the subpolar North Atlantic (Fig. 1a).

Figure 2Comparison of different temperature profiles of the water column near the studied core site. The grey line shows the water profile at 5428.53 N, 5604.33 W obtained during cruise MSM45 in August 2015. Coloured lines mark the profiles at the nearest location (27 km from the core site) at 5437.50 N, 5612.5 W, obtained from the World Ocean Atlas 2013 (WOA13; 1955–2012; Locarnini et al., 2013), showing annual mean (black), August (red) and May (green) temperatures. August data obtained during MSM45 and from WOA13 show similar temperature profiles with a warming of up to 6–8 C at the sea surface. The upper 50 m of the water column reflect atmospheric temperatures and therefore differ in the August, May and annual mean profiles from WOA13. All profiles show temperatures of −1 to 1 C between ca. 40 and 200 m, reflecting the cold, inner LC and higher temperatures below 200 m depth, rising up to 4 C at the shelf bottom, which reflects the influence of the warmer, outer LC strongly influenced by the WGC.

3 Material and methods

3.1 Material

Gravity core MSM45-31-1 (1150 cm core length) was recovered from the Labrador Shelf at 5424.74 N, 5600.53 W, at 566 m water depth (Fig. 1a), during the R/V Maria S. Merian cruise MSM45 in August 2015 (Schneider et al., 2016). The core site is located next to the Cartwright Saddle, off Hamilton Inlet, in a depression associated with an over 600 m deep basin that is part of the Labrador Marginal Trough along the inner Labrador Shelf. The sediment column consists of homogenous olive grey, silty, clayey mud and was subsampled at 5 cm intervals on board the research vessel. A CTD (conductivity–temperature–depth) deployment at nearby station 30 at 5428.53 N, 5604.33 W provides depth profiles of temperature and salinity (Fig. 1b) of southern Labrador Shelf waters at the time of collection. A downward increase in temperature and salinity below 200 m depth shows that the inner LC overlies warmer and more saline waters of the outer LC (Figs. 1b, 2).

Table 1Twelve AMS radiocarbon dates, calibrated in Calib7.1 (Stuiver and Reimer, 1993; Stuiver et al., 2017) using the Marine13 dataset (Reimer et al., 2013) with a reservoir correction of ΔR=144±38 years (McNeely et al., 2006).

Relative probability: a 0.192; b 0.808, according to Calib7.1.

Download Print Version | Download XLSX

Figure 3Age–depth model of core MSM45-31-1. Black triangles mark the calibrated ages of the AMS radiocarbon dates. The black dotted line marks the sedimentation rate (cm yr−1).


3.2 AMS radiocarbon dating

The stratigraphy of MSM45-31-1 is based on 12 accelerator mass spectrometry (AMS) 14C measurements on mixed calcareous benthic foraminifera at the Leibniz Laboratory of Kiel University (CAU), Germany. Each sample contained over 1000 specimens (about 5 mg), which were picked from 1 cm thick sediment slices at 1 m intervals down-core. The 14C dates were calibrated using Calib 7.1 (Stuiver and Reimer, 1993; Stuiver et al., 2017) and the Marine13 dataset (Reimer et al., 2013) with a reservoir correction (ΔR) of 144±38 years, based on the present Labrador Shelf marine radiocarbon correction (McNeely et al., 2006). Dates are reported in calibrated years BP (before present) and are presented in Table 1.

3.3 Alkenone biomarkers

In total, 218 bulk sediment samples were analysed at 5 cm intervals at the Biomarker Laboratory, Institute of Geosciences, Kiel University. Long-chained alkenones (C37) were extracted from 2 to 3 g homogenized bulk sediment, using an accelerated solvent extractor (Dionex ASE-200) with a mixture of 9:1 (vv) of dichloromethane : methanol (DCM:MeOH) at 100 C and 100 bar N2 (g) pressure for 20 min. Extracts were cooled at −20C and brought to near dryness by Syncore polyvap at 40 C and 490 mbar. For the identification and quantification of C37:2, C37:3 and C37:4, we used a multi-dimensional, double-column gas chromatography (MD-GC) set-up with two Agilent 6890 gas chromatographs (Etourneau et al., 2010). The addition of an internal standard prior to extraction (cholestane [C27H48] and hexatriacontane [C36H74]) allowed the quantification of the organic compounds that are reported in nanograms per gram dry bulk sediment. The proportion of each unsaturated ketone was obtained by the integration of peak areas of the different compounds in the respective chromatograms. Alkenone concentrations (sum of C37:2, C37:3 and C37:4) are used as an indication of alkenone productivity.

3.3.1 Sea surface temperature estimates

Sea surface temperature (SST) estimates are based on the U37K index (Brassell et al., 1986) and the U37K index (Prahl and Wakeham, 1987). The U37K index was calculated according to Brassell et al. (1986): U37K= (C37:2-C37:4)/(C37:2+C37:3+C37:4). The U37K index was calculated using the calibration of Prahl and Wakeham (1987): U37K= (C37:2)/(C37:2+C37:3). Both result in standard deviations of 0.01. Rosell-Melé (1998) and Bendle and Rosell-Melé (2004) point out that U37K-based SST estimates down to 6 C are more congruent with modern observations than those based on U37K. However, U37K-based SST estimates could be biased with respect to absolute values when C37:4 exceeds 5 % (Rosell-Melé, 1998), which is the case here, and should therefore be interpreted with caution.

3.3.2 %C37:4 – a proxy for meltwater or sea ice cover?

Several studies have suggested an increase in the production of tetra-unsaturated C37 (C37:4) in polar and subpolar surface waters at lower salinities (Bendle et al., 2005; Blanz et al., 2005; Rosell-Melé, 1998; Rosell-Melé et al., 2002; Sicre et al., 2002). Thus, the proportion of C37:4 relative to the sum of alkenones [%C37:4=C37:4/(C37:2+C37:3+C37:4)] in sedimentological records from the North Atlantic is often used to qualitatively indicate changes in salinity. However, several studies have observed no or even an inverse relationship between %C37:4 and salinity (Schwab and Sachs, 2011; Sikes and Sicre, 2002; Filippova et al., 2016; Theroux et al., 2010; Toney et al., 2010). Bendle et al. (2005) found that higher proportions of C37:4 are produced in low-productivity polar waters, probably because C37:4 proportions vary in different alkenone-producing haptophyte species (Conte et al., 1995; Marlowe et al., 1984). A culture experiment by Chivall et al. (2014) showed that certain, non-coccolith-bearing, prymnesiophytes like Chrysotila lamellosa produce between 30 % and 44 % C37:4 in a positive relationship to salinity. Nonetheless, according to Bendle et al. (2005), the highest C37:4 concentrations in surface sediments from the polar Atlantic delineate the mean position of the sea ice margin off Greenland. In addition, preliminary core top studies imply that high (> 15 %) abundances of C37:4 are restricted to surface sediments within the limits of seasonal sea ice extent in Baffin Bay and the Labrador Sea (Thomas Blanz, unpublished data). Although the exact environmental factors that affect the production of C37:4 remain ambiguous, we interpret higher proportions of C37:4 (> 15 %) as an indicator of sea ice margin conditions.

3.4Mg∕Ca measurements

For Mg∕Ca measurements, the benthic foraminifera species Nonionellina labradorica was used. Due to its infaunal lifestyle, the bottom water carbonate ion concentration has been shown to have a limited influence and this species has, in the Arctic, been deemed a suitable recorder of Mg∕Ca in shelf regions (Barrientos et al., 2018). Per sample, 20–70 specimens of N. labradorica were handpicked from the > 315 µm size fraction, weighed and crushed between two glass plates. Two-thirds of the crushed sample were transferred into pre-leached Eppendorf vials and cleaned following the full protocol of Martin and Lea (2002), including a reductive and oxidative cleaning step and a final leaching step with 0.001 N HNO3. After dissolving and diluting the samples in 0.1 N HNO3, they were measured with an inductively coupled plasma–optical emission spectroscopy (ICP-OES) instrument with radial plasma observation at the Institute of Geoscience, Kiel University. The analytical error for Mg∕Ca analyses was 0.1 % relative standard deviation. The JCP-1 coral reference material (Hathorne et al., 2013) was measured each sixth sample for accuracy and drift correction. Additional accuracy control was provided by measurements of standard BAM RS3 (Greaves et al., 2008) each 60th sample. Trace elements (Fe, Al and Mn) were monitored to exclude possibly contaminated or coated samples from the dataset. Mg∕Ca values exceeding 2.4 mmol mol−1 were considered unreliable and are marked as flyers in the dataset (Fig. 4d). Based on seven duplicate down-core sample measurements, we obtained an analytic precision of 0.2 mmol mol−1Mg∕Ca, which translates into 2 C in respective temperature estimates. Bottom water temperature (BWT) reconstructions are based on the calibration of Skirbekk et al. (2016) for the subpolar North Atlantic region. As the calibration is based on a temperature range of 1–4 C, estimates exceeding 4 C may be less accurate, which is the case for about 35 % of the samples reported here. As a result, we refrain from interpreting individual data points exceeding 4 C. However, given the consistency of average Mg∕Ca values for different sections of the core, we are confident in our interpretation of changing average BWT for specific time intervals (see Sect. 5.2.3 below).

Figure 4Down-core record of MSM45-31-1 showing (a) alkenone sum, (b) U37K index displayed with 0.01 uncertainty, and respective SST estimates, (c) U37K index displayed with 0.01 uncertainty and respective SST estimates, (d) %C37:4 displayed with 0.8 % uncertainty, (e) Mg∕Ca (mmol mol−1) of N. labradorica displayed with 0.19 mmol mol−1 uncertainty and respective BWT estimates, (f) δ18Ocalcite (‰-V-PDB) of N. labradorica displayed with 0.11 ‰ uncertainty, and (g) δ18Ow-ivc (‰-V-SMOW) displayed with 0.35 ‰ uncertainty. Grey vertical bars highlight periods of pronounced oceanographic change.


3.5 Stable isotope measurements

For stable isotope measurements, the remaining one-third of the crushed samples of N. labradorica were cleaned with ethanol absolute, decanted and dried at 40 C. Stable isotope analyses were carried out at the Leibniz Laboratory for Radiometric Dating and Stable Isotope Research in Kiel. A Finnigan MAT 253 mass spectrometer coupled with a Kiel IV carbonate preparation device was used. Results were calibrated to the Vienna Pee Dee Belemnite (V-PDB) scales, which is applied to all calcitic δ18O (δ18Oc) analyses. Based on 11 down-core duplicate measurements, the standard deviation is 0.11 ‰. To estimate the δ18O signature of seawater (δ18Ow, ‰ Vienna Standard Mean Ocean Water (V-SMOW)), the value of δ18Oc was first corrected for vital effects of N. labradorica by subtracting 0.15 ‰ (Rasmussen and Thomsen, 2009) and then applied in the Shackleton et al. (1974) equation together with the Mg∕Ca-derived temperature estimates. A sea level correction of 0.0083 ‰ m−1 (1 ‰ for 120 m of sea level), derived from the relative sea level curve by Austermann et al. (2013), was applied, to correct for the effect of global ice volume on the ocean's δ18O signature (δ18Ow-ivc). Based on duplicates, the standard deviation of δ18Ow-ivc is 0.35 ‰.

4 Results

4.1 Chronology

The chronology of core MSM45-31-1 is based on a linear interpolation of the calibrated median probability ages that are presented in Table 1. Based on the lowermost and uppermost age constraints, the core covers at least the period from 5801 (+67/-58) to 201 (+63/-29) years BP. The age–depth model is presented in Fig. 3. The sedimentation rate varies between 0.11 and 0.35 cm yr−1 and is lowest in the deeper part of the core, from 1150 to 900 cm depth (Fig. 3), probably due to compressional effects resulting from gravity coring. The chronological resolution is between 3 and 9 yr cm−1, so the core provides a sub-decadal record of the late Holocene.

4.2 Temporal variability in southern Labrador Shelf waters

Alkenone concentrations range from 20 to 80 ng g−1 and, despite major fluctuations, gently increase on average towards the top of the core (Fig. 4a). The U37K index ranges from 0.30 to 0.43 and respective SST estimates range from 8 to 12 C (Fig. 4b). The U37K index and SST estimates range from 0.05 to 0.25 and from 0 to 6 C, respectively (Fig. 4c). The proportion of C37:4 ranges from 8 % to 22 % and shows the lowest values between 5.6 and 2.1 ka (Fig. 4d). Mg∕Ca ratios range from 1.3 to 3 mmol mol−1 (excluding values exceeding 2.4 mmol mol−1; see Sect. 3.4 above) (Fig. 4e). BWT estimates range from 0.6 to 9.9 C and indicate a shift to warmer temperatures at 2.1 ka (Fig. 4e). The δ18Oc (‰-V-PDB) of N. labradorica ranges from 3.2 ‰ to 4.0 ‰ and does not show any significant minima or maxima (Fig. 4f). The temperature and ice-volume-corrected δ18Ow-ivc (‰-V-SMOW) record ranges from −0.3 ‰ to 2.9 ‰ and shows a small shift to more positive values after 2.1 ka (Fig. 4g). According to major changes in surface and bottom water conditions, the record has been divided into three main environmental intervals: from 6.2 to 5.6 ka, from 5.6 to 2.1 ka and from 2.1 ka to the present (Fig. 4).

4.2.1 From 6.2 to 5.6 ka (1150–1085 cm)

Between 6.2 and 5.6 ka, the alkenone concentration is relatively low with average values of 33 ng g−1. The U37K and U37K indices and relative SST estimates show the lowest values of the record, suggesting a period of cool surface waters on the Labrador Shelf, while the proportion of C37:4 is high at about 18 %. Mg∕Ca levels and BWT estimates are also relatively low at 1.6 mmol mol−1 and about 3 C, respectively. The δ18Oc is at about 3.5 ‰ and the δ18Ow-ivc is at about 0 ‰.

4.2.2 From 5.6 to 2.1 ka (1080–435 cm)

Between 5.6 and 2.1 ka, alkenone concentrations vary between 25 and 80 ng g−1. Two peaks up to 60 and 55 ng g−1 at 5.4 and 3.8–3.0 ka are evident (Fig. 4a). U37K indices are in medium ranges, while U37K indices are in maximum ranges in this interval, leading to SST estimates of 8–10 and 2–6 C, respectively. U37K-based SST estimates show two peaks of 5.8 and 6 C at 5.4 and 3.8–3.0 ka. The proportion of C37:4 shows the lowest values in this interval with an average of 12 % and two minima of 9 % and 8 % at 5.4 and 3.1 ka. Despite the warming of surface waters, BWT reconstructions remain at average values of 3.1 C, which corresponds to modern conditions at the core site (Figs. 1b, 2). However, the BWT record shows two warm peaks of 7.7 and 9.8 C at 5.4 and 5.2 ka, respectively, and a third warm peak to 7 C at 3.8 ka. The δ18Oc record shows a minor long-term increase with average values of 3.6 ‰. The δ18Ow-ivc record shows an increasing trend towards positive values of about 0.5 ‰, except for a maximum of about 2 ‰ at 5.2 ka, suggesting a short increase in bottom water salinity.

4.2.3 From 2.1 ka to the present (435–0 cm)

The interval after 2.1 ka shows an increase in alkenone concentration that peaks at 75 ng g−1 at 0.7 ka, with average values of 48 ng g−1. The U37K-based SST record shows higher fluctuation in this interval with temperatures between 8.5 and 11.5 C, while the U37K-based SST record reveals a general decrease with fluctuations between 0 and 6 C and average values of 2.2 C. The %C37:4 record is characterized by a shift to values of 17 % on average, suggesting an increase in sea ice after 2.1 ka. Especially in this interval, the use of either U37K or U37K, results in opposite SST trends. Nonetheless, the use of the U37K index results in SST estimates closer to the range of modern values in the upper water column (Fig. 2) than those provided by the U37K index. The BWT reconstructions display higher values after 2.1 ka, with average temperatures of 6 C. The δ18Oc record continues with a slight increase to more positive levels, while the δ18Ow-ivc record shows higher fluctuations and generally more positive values in this interval as well as a peak to nearly 3 ‰, in correspondence with the bottom water warm peaks.

5 Discussion

5.1 The reliability of sea surface and bottom water temperature estimates with regard to modern conditions

U37K- and U37K-based SST estimates differ substantially in this region, due to relatively large proportions of C37:4. While U37K-based SST estimates with an average of 9 C appear too warm for this region, U37K-based SST estimates with an average of 3 C are congruent with modern observations of SSTs at the core site. Although the SST measured in August 2015 was 6 C (Fig. 1b), the annual mean SST at the core site is below 2 C (Fig. 2). However, alkenones neither reflect peak summer (i.e. August) nor annual mean temperatures but rather record the average temperatures of the plankton blooming season, which on the Labrador Shelf, today, is between May and November (Frajka-Williams and Rhines, 2010; Harrison et al., 2013). Thus, at our core site, average U37K-based SST estimates of our down-core record (3.3 C) as well as of the most recent core top sample (3.3 C) agree very well with the modern average SST between May and November of 3.5 C (WOA13; 1955–2012; Locarnini et al., 2013). Therefore, we assume that the U37K index provides the most accurate SST estimates for our record, although we are aware that alkenone-based temperature estimates become less robust with increasing values of %C37:4 and need to be interpreted with caution (Rosell-Melé, 1998).

Figure 5Surface (red) and bottom (blue) water temperature reconstructions are based on the alkenone U37K index and Mg∕Ca  ratios of benthic foraminifera N. labradorica, respectively. In contrast to the surface record, the BWT record is displayed with a three-point running mean due to the 2 C uncertainty of individual data points. Both records show a similar temperature range between 6.1 and 2.1 ka. After 2.1 ka, the BWT record displays higher fluctuations and a shift to warmer average temperatures, leading to a temperature reversal in the water column with generally warmer bottom waters underlying colder surface waters in this interval. The warming of bottom waters is associated with a stronger inflow of the westward retroflection of the WGC reaching our core site at 566 m depth underneath the cold inner LC. The grey vertical bar highlights a period of pronounced oceanographic change.


In our record, the down-core mean BWT of 3.8 C is slightly above the mean value of 3.3 C in U37K-based SST estimates, indicating a water column temperature inversion (Fig. 5). Although our BWT estimates may be biased due to the small temperature range of the calibration used, they are consistent with modern annual mean temperatures above 3.5 C at about 500 m depth (Fig. 2). While a water column temperature inversion appears counterintuitive, it is likely caused by inflowing warmer and more saline subsurface waters of the outer branch of the LC that receives a large proportion of the westward retroflection of the WGC, which, in turn, is a mixture of the cold EGC and the warm, saline IC (Fig. 1a). The outer LC, fed by the WGC, reaches the deeper regions on the Labrador Shelf, while the upper 200 m of the water column are dominated by the cold and shallow inner LC, which carries a higher proportion of Arctic waters from Baffin Bay. Thus, we interpret the U37K-based SST estimates as reflecting the inner LC, while Mg∕Ca-based BWT estimates instead reflect conditions of the outer LC and the westward retroflection of the WGC.

5.2 Middle to late Holocene phases of Labrador Shelf oceanography

Based on the surface and bottom water records presented here, we differentiate between three main climatic and oceanographic intervals in Labrador Shelf waters during the middle to late Holocene. From 6.2 to 5.6 ka, the LC experienced a cool period with a strong sea ice cover that is also evident along the western coast of Greenland (Disko Bugt; e.g. Moros et al., 2016). This cold interval was followed by generally warmer conditions in Labrador Shelf surface waters between 5.6 and 2.1 ka associated with the Holocene thermal maximum in the western Arctic (Kaufman et al., 2004). In addition to the gradual cooling trend of the last 3000 years, our %C37:4 record implies a recurrence of stronger sea ice cover after 2.1 ka. Simultaneously, a shift to warmer BWT after 2.1 ka likely corresponded to an increased supply of warmer Atlantic waters with the westward retroflection of the WGC reaching deeper regions on the Labrador Shelf.

5.2.1 Cold Labrador Current (6.2–5.6 ka)

From 6.2 to 5.6 ka, low sea surface and bottom water temperatures suggest that both surface and subsurface waters were dominated by cold, Arctic water masses from Baffin Bay, with a significant sea ice cover indicated by relatively high %C37:4 (Fig. 6e, f, g). We associate this period with a cold inner LC caused by strong advection of Arctic waters and sea ice from Baffin Bay. Colder conditions with a longer annual sea ice cover have also been detected in a northern Labrador fjord between 7 and 5.8 ka (Richerol et al., 2016). A similar scenario is also seen south of our core in Trinity Bay, Newfoundland, between 7.2 and 5.5 ka (Sheldon et al., 2015), with a period of cold surface waters and seasonal sea ice cover, most likely sourced from the Arctic via the LC. Moreover, Solignac et al. (2011) estimated cold SSTs around Newfoundland during the mid-Holocene and ascribed these conditions to a stronger than present LC, which may have been caused by accelerated atmospheric (westerly wind) circulation patterns (Jessen et al., 2011) combined with strong melting in the Arctic. In contrast to this scenario, however, a decline in sortable silt (SS) mean size on the northern Labrador Shelf has been interpreted as reflecting a relatively weak LC during this period (Site Hu2006040-40, Rashid et al., 2017, Fig. 6c).

Figure 6(a) The NAO index (Olsen et al., 2012) with positive (red) and negative (blue) phases. (b) %C37:4 in core 343300 (Moros et al., 2016). (c) Sortable silt mean (µm) in core Hu2006-040-040 (Rashid et al., 2017). (d) SST record in Newfoundland core AI07-12G (Sicre et al., 2014). (e) SST record based on the U37K index (this study). (f) Higher levels of %C37:4 are interpreted as reflecting a longer sea ice season or meltwater. (g) BWT estimates are based on Mg∕Ca ratios in benthic foraminifera N. labradorica and reflect subsurface conditions. (h) Lower abundances of polar water planktic foraminifera species N. pachyderma sinistral in central Labrador Sea core RAPiD-35-COM indicate warming due to stronger IC inflow (Moffa-Sánchez and Hall, 2017). (i) Mg∕Ca temperature estimates representing the Deep Western Boundary Current (Marchitto and deMenocal, 2003). Grey vertical bars highlight periods of pronounced oceanographic change.


On the north-eastern rim of the Labrador Sea, in Disko Bugt (Greenland), high amounts of %C37:4 were also evident around 6 ka (Moros et al., 2016; Fig. 6b), which are in agreement with our findings. The authors interpreted these high proportions of C37:4 as indicating strong meltwater run-off. At the same time, between 6.2 and 5.5 ka, they identified a severe drop in warm water benthic foraminifera species Islandiella norcrossi, suggesting an abrupt subsurface water cooling (Moros et al., 2016). The conditions in Disko Bugt may well have been related to a cold WGC associated with a cooling in the EGC at about this time (Perner et al., 2015; Müller et al., 2012; Ran et al., 2006). A distinct cooling episode between 7 and 5 ka is also seen in Greenland ice core data (O'Brien et al., 1995), which was likely related to changes in atmospheric circulation leading to increased storminess in the north-eastern United States (Noren et al., 2002). A cold spell is evident in the Camp Century ice core record as well, with a pronounced decrease in δ18O values at 5.8–5.6 ka (Vinther et al., 2009), consistent with a temporary cooling–freshening in oceanic conditions affecting north-east Baffin Bay. Additionally, a marked temperature drop in the northern North Atlantic (e.g. Moros et al., 2004; Telesiński et al., 2014) is likely linked to the most pronounced North Atlantic Holocene ice-rafted debris event (Bond et al., 2001). The coincidence of cool periods in records from the western, south-western and eastern Labrador Sea, as well as in many other North Atlantic records (e.g. Bond et al., 2001; Moros et al., 2004; Telesiński et al., 2014), suggests that it was a widespread phenomenon in the western North Atlantic region.

5.2.2 Sudden warming followed by gradual cooling of the Labrador Current (5.6–2.1 ka)

The start of this interval at 5.6 ka is marked by a pronounced increase in SSTs to over 6 C and a decline in %C37:4 to about 10 % (Fig. 6e, f), implying a substantial warming of Labrador Shelf surface waters and a decrease in sea ice cover. Two warm peaks can also be seen in bottom waters at the beginning of this interval, coeval with SST maxima (Figs. 5, 6e, g). Warmer conditions on the Labrador Shelf may be related to the shallowing of the Arctic channels at about 6 ka, which led to a decrease in Arctic water flow to Baffin Bay and instead enhanced the export through the Fram and Denmark straits along the East Greenland coast (Williams et al., 1995). Warmer conditions were also indicated by the benthic foraminiferal assemblage of Trinity Bay, Newfoundland, that showed increases in fresh food supply and primary productivity between 5.7 and 4.85 ka (Sheldon et al., 2015). Further south, off Cape Hatteras, Cléroux et al. (2012) described increased surface water salinities between 5.2 and 3.5 ka, which they attributed to the decreased export of colder and fresh waters from the north, allowing a northward shift of the Gulf Stream. In contrast to this proposed weakening of the LC, Rashid et al. (2017) identified a short-lived increase in SS mean size in northern Labrador Shelf core Hu2006-040-040 (Fig. 6c), suggesting an intensification of the LC. This latter observation is in agreement with dinocyst assemblages of two Newfoundland records from Bonavista and Placentia Bay (Solignac et al., 2011) suggesting fresher and colder sea surface conditions between 5.6 and 4 ka, which the authors link to enhanced meltwater supply during a warm period in the Arctic until 4.5–4 ka.

In Disko Bugt, western Greenland, the interval from 5.5 to 3.5 ka is marked by a return to relatively warm subsurface conditions and a gradual decrease in %C37:4 values following the peak at 6 ka (Moros et al., 2016; Fig. 6b), which is in agreement with our findings. Diatom and dinocyst assemblages both show relatively mild surface water temperatures despite a gradual trend towards cooler conditions, which is also evident in our SST record between 5.5 and 3.5 ka. A reduction in the abundance of N. labradorica indicates lower surface water productivity that is also suggested by dinocyst assemblages (Ouellet-Bernier et al., 2014), which Moros et al. (2016) link to the increased strength and/or warmth of the WGC. The enhanced influence of the WGC may have resulted from a strong and relatively warm IC that is reported from the East Greenland shelf (Jennings et al., 2002, 2011) and south-west and south of Iceland (e.g. Knudsen et al., 2008). The strong influence of the WGC in Disko Bugt may also have been the result of a largely land-based ice sheet and reduced meltwater run-off from the Greenland Ice Sheet (GIS) after 6 ka (Briner et al., 2010; Weidick and Bennike, 2007; Weidick et al., 1990) with a minimum ice extent from ca. 5 to 3 ka (Briner et al., 2014). This period also corresponds to relatively warm air temperatures recorded in Greenland ice cores (e.g. Camp Century, Moros et al., 2016), which were likely a result of the Holocene thermal maximum between 5 and 3 ka over central to southern Greenland (Kaufman et al., 2004). Additionally, several lake records near Jakobshavn Isbræ display high productivity under relatively warm terrestrial conditions, and one of the lakes, North Lake (Axford et al., 2013), indicates relatively high chironomid-based temperatures. High lake levels linked to warmer conditions are also reported in the Kangerlussuaq region, just south of Disko Bugt (Aebly and Fritz, 2009).

In our record, a second warm peak at 3 ka is evidenced by the lowest %C37:4 and the highest SSTs (Fig. 6e, f). The warm temperatures around 3 ka correlate with a period of negative NAO conditions (Olsen et al., 2012; Fig. 6a), which generally result in weaker north-westerly winds and reduced transport of Arctic waters through the Canadian gateways, causing a weaker LC, milder winters and shorter sea ice seasons. Warmer sea surface conditions are also seen in Placentia Bay, Newfoundland, indicated by the lowest numbers of sea ice indicator species Islandinium minutum between 3.2 and 2.2 ka (Solignac et al., 2011). These warmer conditions around Newfoundland correlate with the lowest %C37:4 in our record, suggesting that both regions experienced a reduction in sea ice cover, possibly caused by generally weak north-westerly winds. In agreement with our records, warmer annual mean temperatures were also detected in Greenland ice core GISP2 between 4 and 2 ka (Alley et al., 1999), potentially also related to the predominantly negative NAO conditions at that time. In contrast to this warm episode with likely low sea ice export from the Arctic region, northern Labrador Shelf core Hu2006-040-040 has been interpreted as implying a relatively strong current at this time (Rashid et al., 2017; Fig. 6c). Although the strength of the LC has generally been tied to meltwater drainage and enhanced freshwater fluxes during the early Holocene, no such linkage has been found for the middle to late Holocene (Rashid et al., 2017). Thus, the role of sea ice and freshwater export on the strength of the LC remains uncertain.

After about 3 ka and coeval with the neoglaciation (Kaufman et al., 2004; Vinther et al., 2009), LC surface waters cooled gradually. The onset of the neoglaciation between 3.2 and 2.7 ka was associated with a cooling in the Arctic that reduced meltwater and sea ice export (Scott and Collins, 1996) and/or iceberg production (Andresen et al., 2011). Correspondingly, reduced meltwater and cold bottom water conditions were detected in Disko Bugt (Greenland) between 3.5 and 2 ka (Lloyd et al., 2007). A similar cooling trend was implied by a general decrease in Atlantic water species in Disko Bugt after 3.5 ka (Perner et al., 2011). Lake deposits from southern Greenland indicate dry and cold conditions from 3.7 ka (Andresen et al., 2004) and shortly after 3.0 ka (Kaplan et al., 2002), respectively. On the north Iceland shelf, temperatures dropped significantly at ca. 3.2 ka (Eiríksson et al., 2000; Jiang et al., 2002), and in the North Atlantic, Bond et al. (1997) identified a strong ice-rafted debris event at 2.8 ka. Hence, the gradual cooling seen in our surface temperature record after 3 ka is in agreement with many North Atlantic records that show colder conditions associated with the neoglaciation.

5.2.3 Subsurface warm water intrusions on the Labrador Shelf despite further surface cooling (2.1 ka to the present)

After 2.1 ka, the SST record shows a continuation of the gradual cooling trend that started around 3 ka with the onset of the neoglaciation. This continued cooling corresponds to the shift from predominantly negative to predominantly positive NAO conditions (Olsen et al., 2012; Fig. 6a), which would have resulted in generally stronger north-westerly winds and colder winters in the Labrador Sea region. The colder surface waters with stronger sea ice cover evident in our record during the last 2100 years suggest that the inner LC received a strong supply of Arctic waters from Baffin Bay in response to the neoglacial cooling and generally positive NAO conditions. In support of our findings, colder conditions and an increase in the sea ice duration were also recorded in Trinity Bay, Newfoundland, after 2.1 ka (Sheldon et al., 2015). Additionally, a gradual surface water cooling of about 2 C during the last 2000 years was also evident in another core from Newfoundland (AI07-12G; Sicre et al., 2014; Fig. 6d), which further supports the surface water cooling trend observed in our record. Coeval with this neoglacial cooling trend, northern Labrador Shelf core Hu2006-040-040 shows a gradual decrease in SS mean size after about 3 ka (Fig. 6c), implying an LC weakening. The general positive correlation of our SST record with the SS mean size record in Hu2006-040-040 implies that colder LC surface temperatures correspond to a weaker current. Although LC strength was tied to early Holocene freshwater discharge, no such linkage is apparent during the middle to late Holocene (Rashid et al., 2017). The authors, however, suggested that the LC weakening of the last 3000 years was related to the neoglacial cooling and enhanced sea ice production, which would have locked up larger amounts of freshwater, resulting in a diminished supply of freshwater and a weakening of the LC. The negative correlation of our %C37:4 record and the SS record of core Hu2006-040-040 does support this scenario of enhanced sea ice cover reducing LC strength. Another explanation for the LC weakening and cooling during the last 3000 years may be the shift to predominantly positive NAO-like conditions (Olsen et al., 2012; Fig. 6a) caused by stronger north-westerly winds that promoted the southward advection of sea ice (Drinkwater, 1996) while also causing a cooling of Labrador Sea surface waters. This would, in turn, lead to a reduced density of offshore waters relative to inshore waters, which would decrease the baroclinic pressure gradient and weaken LC transport (Dickson et al., 1996).

Despite the gradual cooling trend of the LC during the last 3000 years (Fig. 6e), bottom waters show a shift to warmer conditions after 2.1 ka (Fig. 6g). The shift to higher BWTs occurred simultaneously with the increase in %C37:4 (Fig. 6f), implying a strongly stratified water column that was influenced by two different water masses. This scenario can be compared to modern conditions of warmer waters underlying cold waters at our core site (Figs. 1b, 2). Today, the cold inner LC only dominates Labrador Shelf waters down to about 200 m depth (Lazier and Wright, 1993), while the deeper part of the shelf, such as our core site at 566 m depth, is influenced by the warmer outer LC that is largely supplied by the westward retroflection of the WGC. Hence, we ascribe the warmer conditions seen in Labrador Shelf bottom waters after 2.1 ka to a strengthening or warming of the WGC carrying warmer Atlantic waters to the deeper regions on the Labrador Shelf. Episodes of an increased influx of Atlantic sourced water from the IC were also seen in the central Labrador Sea, indicated by decreases in the abundances of polar water planktic foraminifera species Neogloboquadrina pachyderma sinistral (%Nps) at about 2.3–2.1 and 1.6 ka (RAPiD-35-COM; Moffa-Sánchez and Hall, 2017; Fig. 6h). These correspond quite well to the subsurface warming peaks in our record. Despite the generally warmer conditions in our subsurface water record during the last 2100 years, the record displays a trend returning to colder temperatures in this interval. A general cooling trend after 2.2 ka is also evident in core RAPiD-35-COM and has been linked to a weakening in LSW formation (Moffa-Sánchez and Hall, 2017).

While our surface records imply the formation of a cold inner LC with a larger seasonal sea ice cover during the last 2100 years, central Labrador Sea records document enhanced inflow of the warm and saline IC (Moffa-Sánchez and Hall, 2017) that is reflected by subsurface warming and higher salinity on the Labrador Shelf (Fig. 4). According to Gelderloos et al. (2012), cold and fresh surface waters are thought to inhibit convection in the Labrador Sea, while modern observations of warming and increased salinity in central Labrador Sea subsurface waters correlate with episodes of weak Labrador Sea convection (Yashayaev, 2007). Furthermore, extending observational data by proxy records for the last 1200 years, Thornalley et al. (2018) report a strong reduction in LSW formation and Deep Western Boundary Current (DWBC) flow during the last 150 years that seems to correspond to a trend to warmer and saltier subsurface waters in the Labrador Sea (Moffa-Sánchez et al., 2014). According to these studies, the reconstructed surface cooling and subsurface warming and salinity increase on the Labrador Shelf starting at 2.1 ka would strongly suggest that LSW formation in the western Labrador Sea was reduced during the last 2 millennia compared to the mid-Holocene. In contrast, Moffa-Sánchez and Hall (2017) infer the opposite scenario for the last 2000 years with a stronger subpolar gyre (SPG) and intensified LSW formation during warmer eastern Labrador Sea conditions (see Fig. 6h). The latter would be congruent with the observations during more positive NAO conditions resulting in stronger convection (e.g. Zantopp et al., 2017).

Despite the ambiguous evidence of the coupling between LC temperature and LSW formation, it appears that a colder inner LC was coupled with a warmer and saltier IC and WGC during the last 2100 years (see also discussion in Moffa-Sánchez et al., 2019). This link can be explained by a reduced LC strength – as more freshwater would have been locked up in sea ice – resulting in a diminished supply of cold, Arctic waters into the North Atlantic Current (NAC), which, in turn, would have led to a warmer IC. Nonetheless, it is currently still difficult to assess whether hydrographic conditions in the LC have had a significant impact on deep-water formation in the central basin or just responded to general atmospheric conditions across the Labrador Sea region, which may have been the most important driver controlling LSW production during the late Holocene.

Interestingly, Marchitto and deMenocal (2003) reported a significant Little Ice Age (LIA) cooling at intermediate water depths along the path of the DWBC south of Newfoundland, corresponding to the cooling seen in our subsurface water record from the Labrador Shelf (Fig. 6i). This correspondence might indicate a coupling between outer LC and DWBC temperature changes at a centennial timescale with the LSW formation being the link between both water masses. However, overall short-term variability in DWBC temperatures during the last 4000 years does not exactly match our record (Fig. 6i). Therefore, to fully assess the interaction between LC variability and LSW formation, more high-resolution records of changes in LSW production from the western Labrador Sea are required.

6 Conclusions

Overall, middle to late Holocene conditions in Labrador Shelf waters display variability that can partly be related to atmospheric forcing, such as NAO-like conditions. Between 6.2 and 5.6 ka, our records imply a cold episode with strong sea ice cover on the Labrador Shelf that has also been evident in other cores from the North Atlantic region. This cold episode was followed by a generally warmer interval from about 5.6 to 2.1 ka corresponding to a late Holocene thermal maximum, interrupted by colder surface temperatures between about 4.5 and 3.5 ka. After about 3 ka, our SST record shows a gradual cooling trend in association with the neoglaciation, which has also been observed in other surface records from the Labrador Sea region. However, at 2.1 ka our record shows an abrupt shift to enhanced sea ice cover and warmer bottom waters on the Labrador Shelf that corresponds to a shift from predominantly negative to predominantly positive NAO-like conditions. We associate the cooling seen in surface waters with stronger north-westerly winds and harsher winters in the region during positive NAO, while the warming in bottom waters was possibly related to a stronger inflow of the westward retroflection of the WGC in response to a stronger supply of the IC that was seen in the central Labrador Sea. Our record implies that phases of enhanced sea ice cover on the Labrador Shelf corresponded to reduced LC strength during the last 6000 years, which may be have been related to atmospheric conditions controlling the length of the sea ice season as well as regulating the baroclinic pressure gradient that drives LC velocities. During the last 2100 years, a reduced current and diminished freshwater supply via the LC – as more freshwater would have been locked up in sea ice – would have led to a diminished supply of freshwater through the western route to the SPG, which, in turn, could have resulted in saltier and denser NAC and IC (Fig. 1). At the same time, intense cooling of dense SPG surface waters during generally positive NAO conditions would have promoted winter convection and strengthened LSW formation in the central basin. Therefore, it appears that positive NAO conditions may have been responsible for both increased Labrador Shelf sea ice cover during the last 2100 years and reduced freshwater supply and LC strength, while enhancing winter cooling and convection in the central Labrador Sea. The potential indirect effect of enhanced sea ice cover and limited freshwater supply via the LC may thus have played an additional role in enabling deeper convection, forming colder and saltier LSW, a positive feedback mechanism during positive NAO conditions.

Data availability

The data reported in this paper are archived in PANGAEA and available at (Lochte et al., 2019b).

Author contributions

RS, JR and MK guided sediment and water column sampling at sea on the RV Maria S. Merian. AAL performed the Mg∕Ca and stable isotope sample preparation and established the age model. TB performed the biomarker sample preparations and measurements. DGS led the ICP-OES measurements. NA performed the stable isotope measurements. AAL led the data compilation, interpretation and writing of the paper in discussions with RS, JR and MK. All co-authors contributed to improving the paper.

Competing interests

The authors declare that they have no conflict of interest.


We wish to thank the captain and crew of the RV Maria S. Merian as well as the scientific team for their great help during the MSM45 cruise. We thank Karen Bremer and Silvia Koch for technical assistance with the ICP-OES and biomarker analyses, respectively.

Financial support

This study was supported by a PhD fellowship to Annalena A. Lochte through the Helmholtz Research School on Ocean System Science and Technology at the GEOMAR Helmholtz Centre for Ocean Research Kiel (grant no. VH-KO-601) and Kiel University.

Review statement

This paper was edited by Bjørg Risebrobakken and reviewed by three anonymous referees.


Aebly, F. A. and Fritz, S. C.: Palaeohydrology of Kangerlussuaq (Søndre Strømfjord), west Greenland during the last ∼8000 years, Holocene, 19, 91–104,, 2009. 

Alley, R. B., Clark, P. U., Keigwin, L. D., and Webb, R. S.: Making sense of millennial-scale climate change, Geophysical Monograph-American Geophysical Union, 112, 385–394, 1999. 

Andresen, C. S., Björck, S., Bennike, O., and Bond, G.: Holocene climate changes in southern Greenland: evidence from lake sediments, J. Quaternary Sci., 19, 783–795,, 2004. 

Andresen, C. S., McCarthy, D. J., Valdemar Dylmer, C., Seidenkrantz, M. S., Kuijpers, A., and Lloyd, J. M.: Interaction between subsurface ocean waters and calving of the Jakobshavn Isbræ during the late Holocene, Holocene, 21, 211–224,, 2011. 

Austermann, J., Mitrovica, J. X., Latychev, K., and Milne, G. A.: Barbados-based estimate of ice volume at Last Glacial Maximum affected by subducted plate, Nat. Geosci., 6, 553–557,, 2013. 

Axford, Y., Losee, S., Briner, J. P., Francis, D. R., Langdon, P. G., and Walker, I. R.: Holocene temperature history at the western Greenland Ice Sheet margin reconstructed from lake sediments, Quaternary Sci. Rev., 59, 87–100,, 2013. 

Barrientos, N., Lear, C. H., Jakobsson, M., Stranne, C., O'Regan, M., Cronin, T. M., Gukov, A. Y., and Coxall, H. K.: Arctic Ocean benthic foraminifera Mg∕Ca ratios and global Mg∕Ca temperature calibrations: New constraints at low temperatures, Geochim. Cosmochim. Ac., 236, 240–259,, 2018. 

Bendle, J. and Rosell-Melé, A.: Distributions of UK37 and UK37 in the surface waters and sediments of the Nordic Seas: Implications for paleoceanography, Geochem. Geophy. Geosy., 5, 11,, 2004. 

Bendle, J., Rosell-Mele, A., and Ziveri, P.: Variability of unusual distribution of alkenones in the surface waters of the Nordic seas, Paleoceanography, 20, PA2001,, 2005. 

Blanz, T., Emeis, K. C., and Siegel, H.: Controls on alkenone unsaturation ratios along the salinity gradient between the open ocean and the Baltic Sea, Geochim. Cosmochim. Ac., 69, 3589–3600,, 2005. 

Böning, C. W., Scheinert, M., Dengg, J., Biastoch, A., and Funk, A., Decadal variability of subpolar gyre transport and its reverberation in the North Atlantic overturning, Geophys. Res. Lett., 33, 21,, 2006. 

Bond, G., Showers, W., Cheseby, M., Lotti, R., Almasi, P., deMenocal, P., Priore, P., Cullen, H., Hajdas, I., and Bonani, G.: A pervasive millennial-scale cycle in North Atlantic Holocene and glacial climates, Science, 278, 1257–1266,, 1997. 

Bond, G., Kromer, B., Beer, J., Muscheler, R., Evans, M. N., Showers, W., Hoffmann, S., Lotti-Bond, R., Hajdas, I., and Bonani, G.: Persistent solar influence on North Atlantic climate during the Holocene, Science, 294, 2130–2136, ttps://, 2001. 

Brassell, S. C., Eglinton, G., Marlowe, I. T., Pflaumann, U., and Sarnthein, M.: Molecular stratigraphy: a new tool for climatic assessment, Nature, 320, 129–133,, 1986. 

Briner, J. P., Stewart, H. A. M., Young, N. E., Philipps, W., and Losee, S.: Using proglacial-threshold lakes to constrain fluctuations of the Jakobshavn Isbræ ice margin, western Greenland, during the Holocene, Quaternary Sci. Rev., 29, 3861–3874,, 2010. 

Briner, J. P., Kaufman, D. S., Bennike, O., and Kosnik, M. A.: Amino acid ratios in reworked marine bivalve shells constrain Greenland Ice Sheet history during the Holocene, Geology, 42, 75–78,, 2014. 

Chivall, D., M'Boule, D., Sinke-Schoen, D., Sinninghe Damsté, J. S., Schouten, S., and van der Meer, M. T. J.: The effects of growth phase and salinity on the hydrogen isotopic composition of alkenones produced by coastal haptophyte algae, Geochim. Cosmochim. Ac., 140, 381–390,, 2014. 

Clarke, R. A. and Gascard, J.-C.: The formation of Labrador Sea water. Part I: Large-scale processes, J. Phys. Oceanogr., 13, 1764–1778,<1764:TFOLSW>2.0.CO;2, 1983. 

Cléroux, C., Debret, M., Cortijo, E., Duplessy, J. C., Dewilde, F., Reijmer, J., and Massei, N.: High-resolution sea surface reconstructions off Cape Hatteras over the last 10 ka, Paleoceanogr. Paleocl., 27,, 2012. 

Conte, M. H., Thompson, A., Eglinton, G., and Green, J. C.: Lipid biomarker diversity in the coccolithophorid Emiliania huxleyi (Prymnesiophyceae) and the related species Gephyrocapsa oceanica, J. Phycol., 31, 272–282,, 1995. 

Cuny, J., Rhines, P. B., Niiler, P. P., and Bacon, S.: Labrador Sea boundary currents and the fate of the Irminger Sea Water, J. Phys. Oceanogr., 32, 627–647,<0627:LSBCAT>2.0.CO;2, 2002. 

Deser, C., Walsh, J. E., and Timlin, M. S.: Arctic sea ice variability in the context of recent atmospheric circulation trends, J. Climate, 13, 617–633,<0617:ASIVIT>2.0.CO;2, 2000. 

Dickson, R., Lazier, J., Meincke, J., Rhines, P., and Swift, J.: Long-term coordinated changes in the convective activity of the North Atlantic, Prog. Oceanogr., 38, 241–295, 1996. 

Dickson, R. R., Osborn, T. J., Hurrell, J. W., Meincke, J., Blindheim, J., Adlandsvik, B., Vinje, T., Alekseev, G., and Maslowski, W.: The Arctic ocean response to the North Atlantic oscillation, J. Climate, 13, 2671–2696,<2671:TAORTT>2.0.CO;2, 2000. 

Drinkwater, K.: Atmospheric and oceanic variability in the Northwest Atlantic during the 1980s and early 1990s, J. NW Atlantic fishery Sci., Dartmouth NS, 18, 77–97, 1996. 

Eiríksson, J., Knudsen, K. L., Haflidason, H., and Henriksen, P.: Late-glacial and Holocene palaeoceanography of the North Icelandic shelf, J. Quaternary Sci., 15, 23–42,<23::AID-JQS476>3.0.CO;2-8, 2000. 

Etourneau, J., Schneider, R., Blanz, T., and Martinez, P.: Intensification of the Walker and Hadley atmospheric circulations during the Pliocene–Pleistocene climate transition, Earth Planet. Sc. Lett., 297, 103–110,, 2010. 

Filippova, A., Kienast, M., Frank, M., and Schneider, R. R.: Alkenone paleothermometry in the North Atlantic: A review and synthesis of surface sediment data and calibrations, Geochem. Geophy. Geosy., 17, 1370–1382,, 2016. 

Frajka-Williams, E. and Rhines, P. B.: Physical controls and interannual variability of the Labrador Sea spring phytoplankton bloom in distinct regions, Deep-Sea Res. Pt. I, 57, 541–552,, 2010. 

Gelderloos, R., Straneo, F., and Katsman, C. A.: Mechanisms behind the temporary shutdown of deep convection in the Labrador Sea: lessons from the great-salinity anomaly years 1968–71, J. Climate, 25, 6743–6755,, 2012. 

Greaves, M., Caillon N., Rebaubier, H., Bartoli, G., Bohaty, S., Cacho, I., Clarke, L., Cooper, M., Daunt, C., Delaney, M., deMenocal, P., Dutton, A., Eggins, S., Elderfield, H., GarbeSchoenberg, D., Goddard, E., Green, D., Groeneveld, J., Hastings, D., Hathorne, E., Kimoto, K., Klinkhammer, G., Labeyrie, L., Lea, D. W., Marchitto, T., Martínez‐Botí, M. A., Mortyn, P. G., Ni, Y., Nuernberg, D., Paradis, G., Pena, L., Quinn, T., Rosenthal, Y., Russell, A., Sagawa, T., Sosdian, S., Stott, L., Tachikawa, K., Tappa, E., Thunell, R., and Wilson, P. A.: Interlaboratory comparison study of calibration standards for foraminiferal Mg∕Ca thermometry, Geochem. Geophy. Geosy., 9, 8,, 2008. 

Harrison, W. G., Børsheim, K. Y., Li, W. K., Maillet, G. L., Pepin, P., Sakshaug, E., Skogen, M. D., and Yeats, P. A.: Phytoplankton production and growth regulation in the Subarctic North Atlantic: A comparative study of the Labrador Sea-Labrador/Newfoundland shelves and Barents/Norwegian/Greenland seas and shelves, Prog. Oceanogr., 114, 26–45,, 2013. 

Hathorne, E., Gagnon, A., Felis, T., Adkins, J., Asami, R., Boer, W., Caillon, N., Case, D., Cobb, K. M., Douville, E., deMenocal, P., Eisenhauer, A., Garbe‐Schönberg, D., Geibert, W., Goldstein, S., Hughen, K., Inoue, M., Kawahata, H., Kölling, M., Cornec, F. L., Linsley, B. K., McGregor, H. V., Montagna, P., Nurhati, I. S., Quinn, T. M., Raddatz, J., Rebaubier, H., Robinso, L., Sadekov, A., Sherrell, R., Sinclair, D., Tudhope, A. W., Wei, G., Wong, H., Wu, H. C., and You, C.-F.: Interlaboratory study for coral Sr/Ca and other element/Ca ratio measurements, Geochem. Geophy. Geosy., 14, 3730–3750,, 2013. 

Hillaire-Marcel, C., De Vernal, A., Bilodeau, G., and Weaver, A. J.: Absence of deep-water formation in the Labrador Sea during the last interglacial period, Nature, 410, 1073,, 2001. 

Hillaire-Marcel, C., De Vernal, A., and Piper, D. J.: Lake Agassiz final drainage event in the northwest North Atlantic, Geophys. Res. Lett., 34, 15,, 2007. 

Hoffman, J. S., Carlson, A. E., Winsor, K., Klinkhammer, G. P., LeGrande, A. N., Andrews, J. T., and Strasser, J. C.: Linking the 8.2 ka event and its freshwater forcing in the Labrador Sea, Geophys. Res. Lett., 39, 18,, 2012. 

Hoogakker, B. A., McCave, I. N., Elderfield, H., Hillaire-Marcel, C., and Simstich, J.: Holocene climate variability in the Labrador Sea, J. Geol. Soc., 172, 272–277,, 2015. 

Hoogakker, B. A. A., Chapman, M. R., McCave, I. N., Hillaire-Marcel, C., Ellison, C. R. W., Hall, I. R., and Telford, R. J.: Dynamics of North Atlantic Deep Water masses during the Holocene, Paleoceanography, 26, PA4214,, 2011. 

Hurrell, J. W.: Decadal trends in the north atlantic oscillation: regional temperatures and precipitation, Science, 269, 676–679,, 1995. 

Jennings, A., Andrews, J., and Wilson, L.: Holocene environmental evolution of the SE Greenland Shelf North and South of the Denmark Strait: Irminger and East Greenland current interactions, Quaternary Sci. Rev., 30, 980–998,, 2011. 

Jennings, A., Andrews, J., Pearce, C., Wilson, L., and Ólfasdótttir, S.: Detrital carbonate peaks on the Labrador shelf, a 13–7 ka template for freshwater forcing from the Hudson Strait outlet of the Laurentide Ice Sheet into the subpolar gyre, Quaternary Sci. Rev., 107, 62–80,, 2015. 

Jennings, A. E., Knudsen, K. L., Hald, M., Hansen, C. V., and Andrews, J. T.: A mid-Holocene shift in Arctic sea-ice variability on the East Greenland Shelf, Holocene, 12, 49–58,, 2002. 

Jessen, C. A., Solignac, S., Nørgaard-Pedersen, N., Mikkelsen, N., Kuijpers, A., and Seidenkrantz, M. S.: Exotic pollen as an indicator of variable atmospheric circulation over the Labrador Sea region during the mid to late Holocene, J. Quaternary Sci., 26, 286–296,, 2011. 

Jiang, H., Seidenkrantz, M. S., Knudsen, K. L., and Eríksson, J.: Late-Holocene summer sea-surface temperatures based on a diatom record from the north Icelandic shelf, Holocene, 12, 137–147,, 2002. 

Kaplan, M. R., Wolfe, A. P., and Miller, G. H.: Holocene environmental variability in southern Greenland inferred from lake sediments, Quaternary Res., 58, 149–159,, 2002. 

Kaufman, D. S., Ager, T. A., Anderson, N. J., Anderson, P. M., Andrews, J. T., Bartlein, P. J., Brubaker, L. B., Coats, L. L., Cwynar, L. C., Duvall, M. L., Dyke, A. S., Edwards, M. E., Eisner, W. R., Gajewski, K., Geirsdóttir, A., Hu, F. S., Jennings, A. E., Kaplan, M. R., Kerwin, M. W., Lozhkin, A. V., MacDonald, G. M., Miller, G. H., Mock, C. J., Oswald, W. W., Otto-Bliesner, B. L., Porinchu, D. F., Rühland, K., Smol, J. P., Steig, E. J., and Wolfe, B. B.: Holocene thermal maximum in the western Arctic (0–180 W), Quaternary Sci. Rev., 23, 529–560,, 2004. 

Keigwin, L. D., Sachs, J. P., Rosenthal, Y., and Boyle, E. A.: The 8200 year BP event in the slope water system, western subpolar North Atlantic, Paleoceanography, 20, 2,, 2005. 

Khatiwala, S. P., Fairbanks, R. G., and Houghton, R. W.: Freshwater sources to the coastal ocean off northeastern North America: Evidence from H18O∕H16O, J. Geophys. Res., 104, 18241–18255,, 1999. 

Knudsen, K. L., Stabell, B., Seidenkrantz, M. S., Eiriksson, J., and Blake, W.: Deglacial and Holocene conditions in northernmost Baffin Bay: sediments, foraminifera, diatoms and stable isotopes, Boreas, 37, 346–376,, 2008. 

Krawczyk, D., Witkowski, A., Moros, M., Lloyd, J., Kuijpers, A., and Kierzek, A.: Late-Holocene diatom-inferred reconstruction of temperature variations of the West Greenland Current from Disko Bugt, central West Greenland, Holocene, 20, 659–666,, 2010. 

Lazier, J. and Wright, D.: Annual velocity variations in the Labrador Current, J. Phys. Oceanogr., 23, 659–678,<0659:AVVITL>2.0.CO;2, 1993. 

Lazier, J. R. N.: The renewal of Labrador Sea Water, Deep-Sea Res., 20, 341–353,, 1973. 

Lazier, J. R. N.: Temperature and salinity changes in the deep Labrador Sea, 1962–1986, Deep-Sea Res., 35, 1247–1253,, 1988. 

Lewis, C. F. M., Miller, A. A. L., Levac, E., Piper, D. J. W., and Sonnichsen, G. V.: Lake Agassiz outburst age and routing by Labrador Current and the 8.2 cal ka cold event, Quatern. Int., 260, 83–97,, 2012. 

Lloyd, J. M., Park, L. A., Kuijpers, A., and Moros, M.: Early Holocene palaeoceanography and deglacial chronology of Disko Bugt, west Greenland, Quaternary Sci. Rev., 24, 1741–1755,, 2005. 

Lloyd, J. M., Kuijpers, A., Long, A., Moros, M., and Park, L. A.: Foraminiferal reconstruction of mid-to late-Holocene ocean circulation and climate variability in Disko Bugt, West Greenland, Holocene, 17, 1079–1091,, 2007. 

Locarnini, R. A., Mishonov, A. V., Antonov, J. I., Boyer, T. P., Garcia, H. E., Baranova, O. K., Zweng, M. M., Paver, C. R., Reagan, J. R., Johnson, D. R., Hamilton, M., and Seidov, D.: World Ocean Atlas 2013, Volume 1: Temperature, edited by: Levitus, S. and Mishonov, A., NOAA Atlas NESDIS 73, 40 pp.,, 2013. 

Lochte, A. A., Repschläger, J., Seidenkrantz, M. S., Kienast, M., Blanz, T., and Schneider, R. R.: Holocene water mass changes in the Labrador Current, Holocene, 4, 676–690,, 2019a. 

Lochte, A. A., Schneider, R. R., Repschläger, J., Kienast, M., Blanz, T., Garbe-Schönberg, D., and Andersen, N.: Age model, alkenones and stable isotopes of sediment core MSM45-31-1 (MSM45_431-1), PANGAEA,, 2019b. 

Marchitto, T. M. and deMenocal, P. B.: Late Holocene variability of upper North Atlantic Deep Water temperature and salinity, Geochem. Geophy. Geosy., 4, 12,, 2003. 

Marlowe, I. T., Green, J. C., Neal, A. C., Brassell, S. C., Eglinton, G., and Course, P. A.: Long chain (n-C37–C39) alkenones in the Prymnesiophyceae. Distribution of alkenones and other lipids and their taxonomic significance, Brit. Phycolog. J., 19, 203–216,, 1984. 

Marshall, J. and Schott, F.: Open-ocean convection: Observations, theory, and models, Rev. Geophys., 37, 1–64,, 1999. 

Martin, P. A. and Lea, D. W.: A simple evaluation of cleaning procedures on fossil benthic foraminiferal Mg∕Ca, Geochem. Geophy. Geosy., 3, 8401,, 2002. 

McNeely, R., Dyke, A. S., and Southon, J. R.: Canadian marine reservoir ages: preliminary data assessment, Geol. Survey of Canada, 2006. 

Moffa-Sánchez, P. and Hall, I. R.: North Atlantic variability and its links to European climate over the last 3000 years, Nat. Commun., 8, 1726,, 2017. 

Moffa-Sánchez, P., Hall, I. R., Barker, S., Thornalley, D. J., and Yashayaev, I.: Surface changes in the eastern Labrador Sea around the onset of the Little Ice Age, Paleoceanography, 29, 160–175,, 2014. 

Moffa-Sánchez, P., Moreno-Chamarro, E., Reynolds, D. J., Ortega, P., Cunningham, L., Swingedouw, D., Amrhein, D. E., Halfar, J., Jonkers, L., Jungclaus, J. H., Perner, K., Wanamaker, A., and Yeager, S.: Variability in the northern North Atlantic and Arctic oceans across the last two millennia: A review, Paleoceanogr. Paleocl., 34, 1399–1436,, 2019. 

Moros, M., Emeis, K., Risebrobakken, B., Snowball, I., Kuijpers, A., McManus, J., and Jansen, E.: Sea surface temperatures and ice rafting in the Holocene North Atlantic: climate influences on northern Europe and Greenland, Quaternary Sci. Rev., 23, 2113–2126,, 2004. 

Moros, M., Lloyd, J. M., Perner, K., Krawczyk, D., Blanz, T., Kuijpers, A., Jennings, A. E., Witkowski, A., Schneider, R., and Jansen, E.: Surface and sub-surface multi-proxy reconstruction of middle to late Holocene palaeoceanographic changes in Disko Bugt, West Greenland, Quaternary Sci. Rev., 132, 146–160,, 2016. 

Møller, H. S., Jensen, K. G., Kuijpers, A., Aagaard-Sørensen, S., Seidenkrantz, M.-S., Prins, M., Endler, R., and Mikkelsen, N.: Late-Holocene environment and climatic changes in Ameralik Fjord, southwest Greenland: evidence from the sedimentary record, Holocene, 16, 685–695,, 2006. 

Müller, J., Werner, K., Stein, R., Fahl, K., Moros, M., and Jansen, E.: Holocene cooling culminates in sea ice oscillations in Fram Strait, Quaternary Sci. Rev., 47, 1–14,, 2012. 

Myers, P. G.: Impact of freshwater from the Canadian Arctic Archipelago on Labrador Sea water formation, Geophys. Res. Lett., 32, 6,, 2005. 

Noren, A. J., Bierman, P. R., Steig, E. J., Lini, A., and Southon, J.: Millennial-scale storminess variability in the northeastern United States during the Holocene epoch, Nature, 419, 821–824,, 2002. 

O'Brien, S. R., Mayewski, P. A., Meeker, L. D., Meese, D. A., Twickler, M. S., and Whitlow, S. I.: Complexity of Holocene climate as reconstructed from a Greenland ice core, Science, 270, 1962–1964,, 1995. 

Olsen, J., Anderson, N. J., and Knudsen, M. F.: Variability of the North Atlantic Oscillation over the past 5,200 years, Nat. Geosci., 5, 808–812,, 2012. 

Ouellet-Bernier, M. M., de Vernal, A., Hillaire-Marcel, C., and Moros, M.: Paleoceanographic changes in the Disko Bugt area, West Greenland, during the Holocene, Holocene, 24, 1573–1583,, 2014. 

Perner, K., Moros, M., Lloyd, J. M., Kuijpers, A., Telford, R. J., and Harff, J.: Centennial scale benthic foraminiferal record of late Holocene oceanographic variability in Disko Bugt, West Greenland, Quaternary Sci. Rev., 30, 2815–2826,, 2011. 

Perner, K., Moros, M., Lloyd, J. M., Jansen, E., and Stein, R.: Mid to late Holocene strengthening of the East Greenland Current linked to warm subsurface Atlantic water, Quaternary Sci. Rev., 129, 296–307,, 2015. 

Prahl, F. G. and Wakeham, S. G.: Calibration of unsaturation patterns in long-chain ketone compositions for palaeotemperature assessment, Nature, 330, 367–369,, 1987. 

Ran, L., Jiang, H., Knudsen, K. L., Eiríksson, J., and Gu, Z.: Diatom response to the Holocene climatic optimum on the North Icelandic shelf, Mar. Micropaleontol., 60, 226–241,, 2006. 

Rashid, H., Piper, D. J., Lazar, K. B., McDonald, K., and Saint-Ange, F.: The Holocene Labrador Current: Changing linkages to atmospheric and oceanographic forcing factors, Paleoceanography, 32, 498–510,, 2017. 

Rasmussen, T. L. and Thomsen, E.: Stable isotope signals from brines in the Barents Sea: Implications for brine formation during the last glaciation, Geology, 37, 903–906,, 2009. 

Reimer, P. J., Bard, E., Bayliss, A., Beck, J. W., Blackwell, P. G., Ramsey, C. B., Buck, C. E., Cheng, H., Edwards, R. L., Friedrich, M., Grootes, P. M., Guilderson, T. P., Haflidason, H., Hajdas, I., Hatté, C., Heaton, T. J., Hoffmann, D. L., Hogg, A. G., Hughen, K. A., Kaiser, K. F., Kromer, B., Manning, S. W., Niu, M., Reimer, R. W., Richards, D. A., Scott, E. M., Southon, J. R., Staff, R. A., Turney, C. S. M., and van der Plicht, J.: IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP, Radiocarbon, 55, 1869–1887,, 2013. 

Rhein, M., Kieke, D., Hüttl-Kabus, S., Roessler, A., Mertens, C., Meissner, R., Klein, B., Böning, C. W., and Yashayaev, I.: Deep water formation, the subpolar gyre, and the meridional overturning circulation in the subpolar North Atlantic, Deep-Sea Res. Pt. II, 58, 1819–1832,, 2011. 

Richerol, T., Fréchette, B., Rochon, A., and Pienitz, R.: Holocene climate history of the Nunatsiavut (northern Labrador, Canada) established from pollen and dinoflagellate cyst assemblages covering the past 7000 years, Holocene, 26, 44–60,, 2016. 

Rosell-Melé, A.: Interhemispheric appraisal of the value of alkenone indices as temperature and salinity proxies in high latitude locations, Paleoceanography, 13, 694–703,, 1998. 

Rosell-Melé, A., Jansen, E., and Weinelt, M.: Appraisal of a molecular approach to infer variations in surface ocean freshwater inputs into the North Atlantic during the last glacial, Global Planet. Change, 34, 143–152,, 2002. 

Schmidt, S. and Send, U.: Origin and Composition of Seasonal Labrador Sea Freshwater, J. Phys. Oceanogr., 37, 1445–1454,, 2007. 

Schmitz, W. J. and McCartney, M. S.: On the north Atlantic circulation, Rev. Geophys., 31, 29–49,, 1993. 

Schneider, R., Blanz, T., Evers, F., Gasparatto, M.-C., Gross, F., Hüls, M., Keul, N., Kienast, M., Lehner, K., Lüders, S., Mellon, S., Merl, M., Reissig, S., Repschläger, J., Salvatteci, R., Schulten, I., Schwarz, J. -P., Schönke, M., Steen, E., Tietjens, A., and van Nieuwenhove, N.: MERIAN-Reports, Paleoclimate Understanding Labrador Sea (PULSE), Cruise No. MSM45, available at: (last access: 14 June 2020), 2016. 

Schwab, V. F. and Sachs, J. P.: Hydrogen isotopes in individual alkenones from the Chesapeake Bay estuary, Geochim. Cosmochim. Ac., 75, 7552–7565,, 2011. 

Scott, D. B. and Collins, E. S.: Late mid-Holocene sea-level oscillation: a possible cause, Quaternary Sci. Rev., 15, 851–856,, 1996. 

Seidenkrantz, M. S.: Benthic foraminifera as palaeo sea-ice indicators in the subarctic realm–examples from the Labrador Sea–Baffin Bay region, Quaternary Sci. Rev., 79, 135–144,, 2013. 

Seidenkrantz, M. S., Aagaard-Sørensen, S., Sulsbrück, H., Kuijpers, A., Jensen, K. G., and Kunzendorf, H.: Hydrography and climate of the last 4400 years in a SW Greenland fjord: implications for Labrador Sea palaeoceanography, Holocene, 17, 387–401,, 2007. 

Sha, L., Jiang, H., Seidenkrantz, M.-S., Li, D., Andresen, C. S., Knudsen, K. L., Liu, Y., and Zhao, M.: A record of Holocene sea-ice variability off West Greenland and its potential forcing factors, Palaeogeogr. Palaeocl., 475, 115–124,, 2017. 

Shackleton, N. J.: Attainment of isotopic equilibrium ocean water and the benthonic foraminifera genus Uvigerina: Isotopic changes in the ocean during the last glacial, Cent. Natl. Rech. Sci. Colloq. Int., 219, 203–209, 1974. 

Sheldon, C., Seidenkrantz, M.-S., Frandsen, P., Vinther Jacobsen, H., Van Nieuwenhove, N., Solignac, S., Pearce, C., Graae Palitzsch, M., and Kuijpers, A.: Variable Influx of West Greenland Current Water into the Labrador Current through the Last 8000 Years, Based on a Multiproxy Study from Trinity Bay, NE Newfoundland, Arktos, 1, 1,, 201. 

Sheldon, C. M., Seidenkrantz, M. S., Pearce, C., Kuijpers, A., Hansen, M. J., and Christensen, E. Z.: Holocene oceanographic changes in SW Labrador Sea, off Newfoundland, Holocene, 26, 274–289,, 2016. 

Sicre, M. A., Bard, E., Ezat, U., and Rostek, F.: Alkenone distributions in the North Atlantic and Nordic sea surface waters, Geochem. Geophy. Geosy., 3, 2,, 2002. 

Sicre, M.-A., Weckström, K., Seidenkrantz, M.-S., Kuijpers, A., Benetti, M., Masse, G., Ezat, U., Schmidt, S., Bouloubassi, I., Olsen, J., Khodri, M., and Mignot, J.: Labrador current variability over the last 2000 years, Earth Planet. Sc. Lett., 400, 26–32,, 2014. 

Sikes, E. L. and Sicre, M. A.: Relationship of the tetra-unsaturated C37 alkenone to salinity and temperature: Implications for paleoproxy applications, Geochem. Geophy. Geosy., 3, 1–11,, 2002. 

Skirbekk, K., Hald, M., Marchitto, T. M., Junttila, J., Klitgaard Kristensen, D., and Aagaard Sørensen, S.: Benthic foraminiferal growth seasons implied from Mg∕Ca-temperature correlations for three Arctic species, Geochem. Geophy. Geosy., 17, 4684–4704,, 2016. 

Solignac, S., Seidenkrantz, M. S., Jessen, C., Kuijpers, A., Gunvald, A. K., and Olsen, J.: Late-Holocene sea-surface conditions offshore Newfoundland based on dinoflagellate cysts, Holocene, 21, 539–552,, 2011. 

Straneo, F.: Heat and Freshwater Transport through the central Labrador Sea, J. Phys. Oceanogr., 36, 606–628,, 2006. 

Straneo, F. and Saucier, F.: The outflow from Hudson Strait and its contribution to the Labrador Current, Deep-Sea Res. Pt. I, 55, 926–946,, 2008. 

Stuiver, M. and Reimer, P. J.: Extended 14C data base and revised CALIB 3.0 14C age calibration program, Radiocarbon, 35, 215–230,, 1993. 

Stuiver, M., Reimer, P., and Reimer, R.: CALIB 7.1, available at:, last access: 27 August 2017. 

Telesiński, M. M., Spielhagen, R. F., and Lind, E. M.: A high-resolution Lateglacial and Holocene palaeoceanographic record from the Greenland Sea, Boreas, 43, 273–285,, 2014. 

Theroux, S., D'Andrea, W. J., and Toney, J.: Phylogenetic diversity and evolutionary relatedness of alkenone-producing haptophyte algae in lakes: implications for continental paleotemperature reconstructions, Earth Planet. Sc. Lett., 300, 311–320,, 2010. 

Thornalley, D. J. R., Oppo, D. W., Ortega, P., Robson, J. I., Brierley, C. M., Davis, R., Hall, I. R., Moffa-Sanchez, P., Rose, N. L., Spooner, P. T., Yashayaev, I., and Keigwin, L. D.: Anomalously weak Labrador Sea convection and Atlantic overturning during the past 150 years, Nature, 556,, 227–230,, 2018.  

Toney, J. L., Huang, Y., and Fritz, S. C.: Climatic and environmental controls on the occurrence and distributions of long chain alkenones in lakes of the interior United States, Geochim. Cosmochim. Ac., 74, 1563–1578,, 2010. 

Ullman, D. J., Carlson, A. E., Hostetler, S. W., Clark, P. U., Cuzzone, J., Milne, G. A., Winsor, K., and Caffee, M.: Final Laurentide ice-sheet deglaciation and Holocene climate-sea level change, Quaternary Sci. Rev., 152, 49–59,, 2016. 

Vilks, G.: Postglacial basin sedimentation on Labrador Shelf, Geological Survey of Canada, Department of Energy, Mines and Resources, 78–28, 17 pp., 1980. 

Vinther, B. M., Buchardt, S. L., Clausen, H. B., Dahl-Jensen, D., Johnsen, S. J., Fisher, D. A., Koerner, R. M., Raynaud, D., Lipenkov, V., Andersen, K. K., Blunier, T., Rasmussen, S. O., Steffensen, J. P., and Svensson, A. M.: Holocene thinning of the Greenland ice sheet, Nature, 461, 385–388,, 2009. 

Wang, H., Legg, S., and Hallberg, R.: The Effect of Arctic Freshwater Pathways on North Atlantic Convection and the Atlantic Meridional Overturning Circulation, J. Climate, 31, 5165–5188,, 2018. 

Weidick, A. and Bennike, O.: Quaternary glaciation history and glaciology of Jakobshavn Isbræ and the Disko Bugt region, West Greenland: a review, Vol. 14, Copenhagen: Geol. Survey of Denmark and Greenland, 2007. 

Weidick, A., Oerter, H., Reeh, N., Thomsen, H. H., and Thorning, L.: The recession of the Inland Ice margin during the Holocene climatic optimum in the Jakobshavn Isfjord area of West Greenland, Palaeogeogr. Palaeocl., 82, 389–399,, 1990. 

Williams, K., Short, S., Andrews, J., Jennings, A., Mode, W., and Syvitski, J.: The eastern Canadian Arctic at ca. 6 ka BP: a time of transition, Géogr. Phys. Quat., 49, 13–27,, 1995. 

Yashayaev, I.: Hydrographic changes in the Labrador Sea, 1960–2005, Prog. Oceanogr., 73, 242–276,, 2007. 

Yashayaev, I. and Loder, J. W.: Enhanced production of Labrador Sea water in 2008, Geophys. Res. Lett., 36, 1,, 2009. 

Yashayaev, I. and Loder, J. W.: Recurrent replenishment of Labrador Sea Water and associated decadal-scale variability, J. Geophys. Res.-Oceans, 121, 8095–8114,, 2016. 

Zantopp, R., Fischer, J., Visbeck, M., and Karstensen, J.: From interannual to decadal: 17 years of boundary current transports at the exit of the Labrador Sea, J. Geophys. Res.-Oceans, 122, 1724–1748,, 2017. 

Short summary
The Labrador Sea is important for the modern global thermohaline circulation system through the formation of Labrador Sea Water. However, the role of the southward flowing Labrador Current in Labrador Sea convection is still debated. In order to better assess its role in deep-water formation and climate variability, we present high-resolution mid- to late Holocene records of sea surface and bottom water temperatures, freshening, and sea ice cover on the Labrador Shelf during the last 6000 years.