Abstract
Benthic foraminiferal assemblages from a core southwest of Disko Bugt provide a Holocene perspective (last ~7 ka BP) on ice-sheet/ocean interactions between the West Greenland Current (WGC) and the West Greenland ice sheet. Changes in the fauna reveal significant variations in the water mass properties (temperature and salinity) of the WGC through time. From 7.3 to 6.2 ka BP, a relatively warm/strong WGC influences ice-sheet melt in Disko Bugt and causes enhanced meltwater production, resulting in low surface-water productivity. The most favourable oceanographic conditions occur from 5.5 to 3.5 ka BP, associated with ‘thermal optimum-like’ conditions, encompassing minimum ice sheet extent in the Disko Bugt area. These conditions are attributed to: (1) reduced meltwater influence as the ice sheet is land based and (2) enhanced contribution of warm/saline water masses from the Irminger Current to the WGC. The transition into the late Holocene (last ~3.5 ka BP) is characterized by a cooling of oceanographic conditions, caused by increased advection of cold/low-salinity water masses from the East Greenland Current. A longer-term late-Holocene cooling trend within the WGC is attributed to the onset of Neoglacial cooling within the North Atlantic region. Superimposed on this cooling trend, multicentennial-scale variability within the WGC matches reconstructions from a nearby coring site in Disko Bugt as follows: (1) cooling at ~2.5 ka BP, linked to the 2.7 ka BP ‘cooling event’; (2) a warm phase centred at 1.8 ka BP, associated with the ‘Roman Warm Period’; (3) slight warming between 1.4 and 0.9 ka BP, linked to the ‘Medieval Climate Anomaly’; (4) severe cooling of the WGC after 0.9 ka BP, culminating at 0.3 ka BP during the ‘Little Ice Age’. We show that multicentennial-scale palaeoceanography variability along the West Greenland margin is driven by ocean forcing, i.e. variations in the relative contribution of Atlantic (Irminger Current) and Polar (East Greenland Current) water masses to the WGC during the last ~7 ka BP, influencing ice sheet dynamics.
Introduction
Over the past two decades a negative mass balance and increased surface melting of the Greenland Ice Sheet (GIS), accompanied by enhanced acceleration of many of Greenland’s marine-terminating outlet glaciers (e.g. Jakobshavn Isbræ, Helheim, Kangerdlugssuaq) has been identified (Howat et al., 2007, 2008, 2011; Joughin et al., 2008; Moon and Joughin, 2008; Rignot and Kanagaratnam, 2006; Straneo et al., 2010; Zwally et al., 2002). Recent studies have suggested that ocean forcing may exert an important control on modern ice sheet dynamics by the influence of warmer ocean conditions (e.g. Bindschadler, 2006; Holland et al., 2008; Joughin et al., 2004; Thomas, 2004).
The Disko Bugt region in West Greenland is a key area to investigate the influence of ocean forcing on GIS behaviour. This region has a relatively wide shelf area and contains high-resolution sedimentary archives with the potential for records to obtain the interaction between oceanographic variability and ice sheet behaviour. The modern hydrographic conditions of the region are dominated by the West Greenland Current (WGC, Figure 1). The water mass composition of the WGC is linked to the large-scale North Atlantic climate system. Recent studies show that on a multidecadal timescale, temperature changes within the WGC have a profound impact on subsurface melting of Disko Bugt outlet glaciers (e.g. Jakobshavn Isbræ), at least during the last 60 years (Holland et al., 2008; Lloyd et al., 2011; Rignot et al., 2010). Ice sheet limits in central West Greenland are uncertain during the last glacial maximum (LGM), but it has been suggested that ice streams may have extended to the shelf edge along deep cross-shelf troughs in the Disko Bugt region (see Funder et al., 2011 and references therein). Limited evidence reported re-appearance of the WGC in Disko Bugt from c. 9 to 10 ka BP (Funder and Weidick, 1991; Lloyd et al., 2005). Marine (Hogan et al., 2011; Lloyd et al., 2005) and terrestrial (Briner et al., 2010; Weidick and Bennike, 2007; Young et al., 2011) studies suggest that Jakobshavn Isbræ had retreated into the Isfjord by c. 7.8 ka BP. There are, however, no high-resolution records available spanning the period of ice retreat (after 8 ka BP) and minimum ice sheet extent (6–4 ka BP; e.g. Briner et al., 2010; Weidick and Bennike, 2007) documenting the palaeoceanography evolution of Disko Bugt (West Greenland).

Schematic bathymetric map of Disko Bugt, adapted from Jakobsson et al. (2008), showing the location of core 343300 (circle) and 343310 (star) in southwestern Egedesminde Dyb and present-day oceanographic setting of the study area. The insert shows the oceanographic setting around Greenland. EGC: East Greenland Current; IC: Irminger Current; WGC: West Greenland Current; LC: Labrador Current.
To assess and understand the link between oceanic forcing and ice sheet dynamics, a longer-term palaeoceanographic perspective is needed. In this study we aim to investigate the oceanographic conditions off West Greenland during the important period when the ice sheet retreated from western to eastern Disko Bugt, and subsequently into the Isfjord. Specifically we aim to assess the potential link between ice margin retreat and ocean temperatures during this period. We focus on a new marine sediment core, from the shelf southwest of Disko Bugt, spanning the last c. 7 ka BP. The coring site is located directly below the flow path of the WGC and is expected to record: (1) meltwater influence from the GIS, received along the West Greenland margin and from Disko Bugt; (2) shifts in the relative contribution of the relatively warm/saline Atlantic (Irminger Current) and colder/fresher Polar (East Greenland Current) water masses to the WGC. Palaeoenvironmental reconstructions are inferred from benthic foraminiferal assemblage data. We use groupings of Atlantic and Arctic water species as a proxy to identify qualitative changes in bottom water mass properties of the WGC (e.g. Lloyd et al., 2011; Perner et al., 2011). Our marine-based reconstructions will provide a high-resolution longer-term Holocene perspective on the palaeoceanographic development of Disko Bugt.
Modern environmental setting of Disko Bugt
Disko Bugt is a large marine embayment (40,000 km2) in central West Greenland (Figure 1). Shallow water depths, varying between 200 and 400 m are typically found, with maximum water depths up to 900 m in Egedesminde Dyb, a deep water trough of glacial origin (Long and Roberts, 2003; Roberts and Long, 2005; Zarudski, 1980). Jakobshavn Isbræ, one of Greenland’s largest outlet glaciers, flows into Disko Bugt and currently drains about 7% of the GIS (Bindschadler, 1984). The present-day oceanographic setting of West Greenland is dominated by the WGC, which is formed by a combination of: (1) relatively warm and saline Atlantic-sourced water from the Irminger Current (IC), a side branch of the North Atlantic Current (NAC); (2) polar-sourced cold, low-salinity water from the East Greenland Current (EGC; Buch, 1981); and (3) local meltwater discharge along the SW Greenland coast (Figure 1). The WGC enters Disko Bugt from the southwest and flows northwards exiting primarily through the Vaigat into Baffin Bay. A branch of the WGC is deflected into Baffin Bay west of Disko Island, while the main current continues to flow along West Greenland northward into northern Baffin Bay (Andersen, 1981; Bâcle et al., 2002; Ribergaard et al., 2006). Temperature and salinity data from Disko Bugt (Andersen, 1981; Buch, 1981; Buch et al., 2004; Harff et al., 2007; Lloyd, 2006) show that the WGC (3.5–4°C, 34.2–34.4 PSU) forms the bottom waters within the bay. Andersen (1981) found no indications of admixture of deep Baffin Bay waters below 300 m water depth, penetrating into Disko Bugt. Surface waters, however, are influenced by meltwater flux from land, icebergs and the previous season’s pack ice, as well as relatively low-salinity polar surface water advected from Baffin Bay. At present, the Arctic sea-ice edge, formed annually in Baffin Bay between September and March, is found just north of Disko Bugt in spring (Tang et al., 2004) and influences the surface water productivity in the area (e.g. Hansen et al., 1999; Levinsen et al., 2000).
Material and methods
This study focuses on the composite record of a multi and gravity core obtained from site MSM 343300 (68°28,311′N/ 54°00,119′W, Figure 1) southwest of Disko Bugt (cruise MSM05/03 of the R/V Maria S. Merian; Harff et al., 2007). The multi core (length: 28.5 cm) and gravity core (total length: 1132 cm, this study focuses only on the depths 0–400 cm) were retrieved from 519 m water depth in Egedesminde Dyb.
Age control is provided by accelerator mass spectrometry AMS14C dates on mollusc shells and benthic foraminifera (Table 1, Figure 2). The chronology of the composite record (last 7 ka BP) is based on 17 AMS 14C dates. The AMS radiocarbon dates were calibrated using the Marine09 (Reimer et al., 2009) calibration curve in CALIB 6.0.2 (Stuiver and Reimer, 1993). Following the results from Lloyd et al. (2011), we applied a marine reservoir age correction ΔR of 140±35 years, which represents the modern ΔR value for the Disko Bugt area.
Radiocarbon dates for gravity core 343300. Uncertainties include 68% of the probability distribution.

Lithological characterization and age–depth model of gravity core 343300. (a) Sediments deposited in the >200 μm fraction (%), dark line, number of counted ice-rafted detritus (IRD) >2 mm, grey line; (b) the 63–200 μm fraction (%); (c) total organic carbon content (TOC %); (d) age model for gravity core 343300, dark line, with depths of the AMS14C radiocarbon dates, dashed lines indicate the lower and upper limit of the age model.
Foraminiferal analysis was carried out on a standard volume of 5 ml fresh sediment, soaked in deionized water overnight and gently sieved at 63 µm just before counting. Multi core samples were counted at 1–2 cm intervals and gravity core samples at 4 cm intervals. Calcareous and agglutinated foraminifera were counted on a squared picking tray and identified to species level under a stereomicroscope from the wet residue >63 µm to reduce the loss of the more fragile arenaceous species caused by drying out of sediment. The total number of specimens counted per sample ranges from 300 to about 750. The Shannon-Wiener Index (S(H)), a measure of faunal diversity, was calculated separately for the total (S(H)total; Figure 3), the agglutinated (S(H)agglutinated; Figure 4), and the calcareous (S(H)calcareous; Figure 5) assemblage.

Total foraminiferal assemblage (calcareous and agglutinated) from site 343300 versus age. Foraminiferal frequencies are expressed as a percentage of the total specimens counted. Only species with abundance greater than 10% are included. Additionally, the total number of benthic foraminifera counted, number of benthic foraminifera counted per millilitre wet sediment, the ratio of calcareous versus agglutinated specimens, number of test linings and grouping of AtlW (red) and AW (blue) indicator species, are presented. The light blue (AW) and light red (AtlW) coloured plots show the foraminiferal concentration (no. of specimens per ml wet sediment) of the respective groups. Additionally, we present the Shannon-Wiener-Index (S(H)) calculated for the total benthic foraminiferal assemblage. Colour figure available online.

Agglutinated foraminiferal assemblage from site 343300 versus age. Foraminiferal frequencies are expressed as a percentage of total agglutinated specimens counted. Only species with abundance greater than 2% are included. Red (blue) colored species are included in the AtlW (AW) group. In addition, the Shannon-Wiener-Index (S(H)) was calculated based on the agglutinated fauna. Colour figures available online.

Calcareous foraminiferal assemblage of site 343300. Foraminiferal frequencies are expressed as a percentage of total calcareous specimens counted. Only species with abundance greater than 2% are included. Red (blue) colored species are included in the AtlW (AW) group. In addition, the Shannon-Wiener-Index (S(H)) was calculated based on the calcareous fauna. Colour figures available online.
Sediment samples (1 cm interval) were wet sieved at the 63 and 200 µm grain size fraction to determine the sand content. The fraction 63–200 µm is used as an approximate measure of the WGC’s current strength and sediments deposited in the fraction >200 µm provide indications of ice-rafted debris (IRD) at the coring site. Additionally, counts of IRD (>2 mm) on x-ray images were carried out. The total organic carbon content (TOC) of the bulk sediment (2 cm interval) is determined indirectly by subtracting the total inorganic carbon form the total carbon content.
Results
Age model and lithology
Sediments are composed of mottled olive/greenish grey to moderate olive brown organic-rich clay with occasional shell fragments and drop stones. The depth–age model was fitted to the calibrated 14C dates using mixed effect modelling (Heegaard et al., 2005). The age model of the multi core is based on linear interpolation between the AMS14C date at 26.5 cm depth (Table 1) and the modern age of 2007, the sampling year, for the core top (see Lloyd et al., 2011). Figure 2 presents the age model for the last 7.3 ka BP of the gravity core. According to our age model, a gap of c. 500 years exists between the multi and gravity core. Loss of the upper sediments is presumably due to the gravity coring technique. The sedimentation rate in the gravity core averages 0.44 cm/yr between 7.3 and 3.5 ka BP, increasing to 0.85 cm/yr between 3.5 and 1 ka BP.
The content of IRD (>2 mm fraction) is low within the core, supported by low values of the >200 µm % fraction throughout the last c. 7 ka BP (Figure 2a). The sand content (63–200 µm % fraction) averages about 8% between 7.3 and 7 ka BP, then increases to maximum values between 7 and 6.3 ka BP (averages 25%; Figure 2b). From about 6 ka BP, there is a gradual decrease in the sand content, reaching an average of about 5% after 3.5 ka BP. The TOC content is initially low with an average of ~1% between 7.3 and 5.5 ka BP, then increases gradually, reaching an average of about 2.5% after 2.5 ka BP (Figure 2c). X-ray radiographs (A Jennings, unpublished data, 2011) reveal two turbidites in the record, one at ~3 ka BP (199.5–201.5 cm core depth; peak in sand content; 0.5% drop in TOC content) and one at ~5 ka BP (291–296.5 cm core depth). Data from these depths have been excluded from the record and the following discussion.
The benthic foraminiferal assemblage and ecology
Benthic foraminifera were counted from 120 samples from the combined cores. A total of 52 benthic foraminiferal species were identified: 17 agglutinated and 35 calcareous species (see Appendix 1 for the complete faunal list). Both calcareous and agglutinated specimens were well preserved and showed minimal evidence of postmortem (dissolution) changes throughout the core. This is supported by relatively low counts of test linings per sample (Figure 3). Following previous studies on high-resolution sites from the Disko Bugt area (e.g. Lloyd et al., 2011; Perner et al., 2011), we present changes in the total benthic foraminiferal assemblage (agglutinated and calcareous assemblage) along with summary curves of a chilled Atlantic water group (AtlW) and an Arctic water group (AW). These groupings are based on environmental preferences of the species (associated directly or indirectly with salinity and temperature; e.g. Murray, 1991; Rytter et al., 2002; Sejrup et al., 2004) and are used to identify changes in the relative temperature and salinity of the WGC associated with changes in the respective water mass composition (IC/EGC influence) during the Holocene.
In Table 2, we present a list of species included in the AtlW and AW group along with references supporting species allocations. The AtlW group includes species such as Islandiella norcrossi and Cassidulina reniforme (calcareous) and Adercotryma glomerata, Reophax pilulifer, and Ammoscalaria pseudospiralis (agglutinated). In the presented study, I. norcrossi is the most abundant species of the AtlW group, indicating a relatively strong IC component of the WGC. The highest abundance of this species is associated with a stable salinity during the most ameliorated oceanic conditions of the Holocene. The AW group includes species such as Cuneata arctica, Spiroplectammina biformis, and Textularia torquata (agglutinated) and Elphidium excavatum f. clavata and Islandiella helenae (calcareous). These species are indicative of relatively fresh and cold water mass characteristics (strong EGC component and/or local meltwater influence from the GIS) within the WGC. The total benthic foraminiferal assemblage also contains various productivity indicator species such as Nonionellina labradorica, Globobulimina auriculata arctica, Buccella frigida, Melonis barleeanus, Epistominella vitrea, Stainforthia loeblichi, Trifarina fluens and Pullenia sp. The occurrence of these species is often related to high productivity at the sea surface, which enhances food supply to the sea floor and the availability of more or less degraded/altered organic matter within the sediments (e.g. Caralp, 1989; Jennings et al., 2004; Mudie et al., 1984; Polyak and Solheim, 1994).
Benthic foraminifera included in the chilled Atlantic water species (AtlW) and Arctic water species (AW).
The relative abundance (%) of the dominant agglutinated and calcareous species based on the total assemblage over the last 7.3 ka BP is presented in Figure 3. Based on changes within the total assemblage, along with changes in AtlW and AW groupings, an informal subdivision into three zones (A–C) has been made. The relative abundance (%) of the agglutinated and calcareous fauna is also plotted separately in Figures 4 and 5. Additionally, we provide absolute abundance of the agglutinated and calcareous fauna (plotted as total specimens counted per millilitre wet sediment) separately in the supplementary section (Figure S1-agglutinated species; Figure S2-calcareous species, available online).
The agglutinated assemblage
Deuterammina ochracea is the most abundant species of the agglutinated assemblage during the last 7.3 ka BP ranging from 40% to 60%. This is a cosmopolitan species that is found widely in Arctic and sub-Arctic environments and provides limited information on palaeoenvironmental change within this study. In zone A, from 7.3 to 6.2 ka BP, AW species (e.g. C. arctica, S. biformis) are relatively common averaging 10% of the fauna. This suggests relatively cold and fresh bottom water mass conditions. At about 6.5 ka BP, there is a slight increase in AtlW taxa (R. pilulifer, S. diflugiformis and A. pseudospiralis, A. glomerata), which is possibly linked to occurrence of relatively warmer water masses at that time.
In zone B, E. advena increases to 15–20% between 6.2 and 5.6 ka BP, and Cribrostomoides sp. averages 10% during this time interval. A minor rise noted in AW species at c. 6.3 ka BP suggests a short-term cooling and freshening of bottom water conditions. From c. 5.5 to 3.5 ka BP, we find the lowest occurrence of the AW group over the period studied, which indicates minimum influence from either local meltwater sources in Disko Bugt or from the EGC component to the WGC during this time interval.
In zone C, from c. 1.7 ka BP onwards, we note a gradual rise in overall abundances and diversity of agglutinated specimens (S(H) Index exceeds 1.5). AtlW species (A. glomerata, A. pseudospiralis, and S. diflugiformis Reophax sp.) increase rapidly from c. 1.6 ka BP, reaching peak abundance of about 30% (Figure 4). Subsequently, from c. 1.6 to 1.3 ka BP, a slight reduction in abundance of this group occurs to about 20%. AW species show a gradual increase from c. 1.5 ka BP, reaching peak values of approximately 30% from c. 0.9 ka BP onwards indicating subsequent freshening and cooling.
The calcareous assemblage
Between 7.3 and 6.3 ka BP (zone A), AtlW species dominate (40–80%), but co-occurrence with AW species (10–40%) presumably reflects mixed (warm and fresh) bottom water conditions. The overall abundance of productivity indicator species (e.g. N. labradorica, B. frigida) is relatively low through this interval. Relatively high abundance of M. barleeanus, compared with following intervals, suggests that relatively old/degraded organic material is present at the site, with limited replenishment from surface productivity at this time. At c. 7 ka BP, peak abundance of the AtlW group (Figure 5; S2) reflects relatively warmer bottom water conditions.
In zone B, peak abundance of N. labradorica, centred at 6 ka BP, indicates a prolonged period of high surface water productivity, causing enhanced supply of fresh phytodetritus to the sea floor. A pronounced rise in AW species (e.g. Elphidium excavatum f. clavata, I. helenae) is also noted with peak abundance at c. 5.5 ka BP, suggesting possible bottom water cooling. Subsequently, a pronounced decline in the abundance of E. excavatum f. clavata is noted, coinciding with an increase in AtlW species. From about c. 5.5 to 3.5 ka BP, the AtlW group dominates the assemblage, documenting relatively warm and ameliorated bottom water conditions in the area (Figure 5; S2).
In zone C, we observe a rise in AW species, in particular E. excavatum f. clavata, but also a minor increase in I. helenae, reaching a peak at c. 2.7 ka BP indicating cooling of bottom water conditions. The abundance of the AtlW group declines through this interval, although a sudden peak in abundance is seen at c. 1.8 ka BP, comparable with abundances found between c. 5.5 and 3.5 ka BP. This indication of a relatively warm WGC is accompanied by low abundance of AW species (e.g. E. excavatum f. clavata) at c. 1.8 ka BP.
The most prominent feature in the uppermost part of the record is the abrupt decrease of E. excavatum f. clavata, I. norcrossi and C. reniforme from about 0.9 ka BP onwards, implying a change to harsher environmental conditions than previously.
Discussion
Long-term Holocene changes in oceanographic variability of the WGC
Our new benthic foraminiferal record, from the shelf southwest of Disko Bugt, illustrates the permanent influence of the WGC, influencing the oceanographic conditions within the area over the last c. 7 ka BP. From our data, the grouping of AtlWcalc and AWagg indicator species, we can extract two main factors that influence water mass properties of the WGC over time: (1) meltwater influence from the GIS, received along the West Greenland margin and from local sources in the Disko Bugt region; (2) shifts in the relative contribution of warm/saline Atlantic (IC) and of colder/fresher Polar (EGC) water masses to the WGC. The variations in local influence from the GIS and regional ocean forcing will be discussed below.
Meltwater influence from the GIS (7.3 to c. 6.2 ka BP) in the Disko Bugt area
Postglacial reappearance of the WGC is reported from the West Greenland area and the Canadian Arctic already after c. 9 ka BP (e.g. Hillaire-Marcel et al., 2001; Knudsen et al., 2008a; Lloyd et al., 2005; Ren et al., 2009). Evidence from marine and terrestrial records suggest that the GIS retreated from the shelf west of Disko Bugt to the eastern part of the embayment by c. 10.2 ka BP (e.g. Lloyd et al., 2005; Long and Roberts, 2002; Long and Roberts, 2003; Young et al., 2011).
Strong melting from the GIS, translated to runoff into the ocean, is reported between 8 and 6 ka BP (Alley and Anandakrishnan, 1995), and is thought to be in response to atmospheric forcing, i.e. pronounced temperature rise over the GIS (Dahl-Jensen et al., 1998; Vinther et al., 2009). In addition to this proposed atmospheric forcing, our reconstructions (zone A) provide evidence of an ocean forcing, with a relatively warm WGC entering Disko Bugt from c. 7.3 to 6.2 ka BP, which presumably supports ice sheet melt in the area. This lower part of the record is also characterized by strong variability in the flow strength of the WGC, displayed by the sand content (Figure 6c). Initially, between 7.3 and 7 ka BP, decreasing sand content suggests a weaker flow of the WGC (Figure 6c). The fauna is dominated by AtlW species, but with variable amounts of AW species, which presumably reflects a mixed (relatively warm and fresh) WGC (Figure 3). A relatively weaker flow of the WGC might be related to meltwater influence from the GIS, received during the travel of the current along the West Greenland margin. This interval coincides also with reduced deposition of biogenic carbon (Figure 6e; MD99-2322) and low abundance of AtlW species on the East Greenland shelf south of Denmark Strait, attributed to continuing and enhanced melting from the GIS (Jennings et al., 2011). Following this initial period, from c. 7 ka BP, significant increase in sand content reflects a strengthening of the WGC flow, and a maximum is seen between 6.5 and 6.2 ka BP (Figure 6c). During this period, postglacial initiation of deep convection is reported from the Labrador Sea (Hillaire-Marcel et al., 2001). The relatively low TOC content and relatively high abundance of M. barleeanus, feeding on altered organic matter (e.g. Caralp, 1989), indicates low in situ marine productivity at this time. We postulate that relatively warmer bottom water flow into Disko Bugt supports/enhances ice sheet retreat and caused enhanced meltwater production that circulated within the embayment at this time. In turn, this meltwater discharge influences surface water productivity, which is translated into reduced concentration of foraminiferal tests (Figure 3) and lower abundance of productivity indicator species feeding on fresh phytodetritus (e.g. N. labradorica, Figures 5 and S2, available online). This influence of surface water productivity is prominently seen, when comparing the relative abundance and concentration (number of specimens per millilitre) of the AtlW group (see Figure 3, red and light red graphs – colour figures available online).

Summary of results compared with other regional data sets. (a) TOC (%) content (343300); (b) age–depth model (dark line, dashed lines show lower and upper limit of the model) and AMS14C dates (343300); (c) sand content (% fraction >63–200 μm; 343300); (d) number of calcareous Atlantic water specimens (AtlWcalc) per millilitre wet sediment, bold line displays data from site 343300 and light line data from nearby site 343310; (e) biogenic carbon (%) content of sediments from site MD99-2322, Denmark Strait (Jennings et al., 2011); (f) number of agglutinated Arctic water species (AWagg) per millilitre wet sediment, bold line displays data from site 343300 and light line data from nearby site 343310; (g) drift ice proxy data (Quartz%) from core site MD99-2269, NW Iceland (Moros et al., 2006a). Known historical climatic events such as the ‘Roman Warm Period’ (RWP), the ‘Medieval Climate Anomaly’ (MCA) and the ‘Little Ice Age’ (LIA) are marked. The black arrows indicate the position of the two turbidites found in the sediments based on x-ray-radiographs (A Jennings, unpublished data, 2011).
The increased flow of the WGC, accompanied by a rise in the concentration of the AtlWcalc (Figure 6d) from c. 7 ka BP, correlates with enhanced deposition of biogenic carbon and an abrupt rise in AtlW species in a core from the southern Denmark Strait after c. 7 ka BP, which Jennings et al. (2011) relate to a warmer and stronger IC (Figure 6e). Our results suggests that meltwater runoff from the GIS started to decline after 7 ka BP and had reduced influence on the water mass composition of the WGC along the West Greenland margin. This is supported by terrestrial reconstructions, which report a largely land-based ice sheet by c. 7 ka BP, which had retreated behind its present margin in eastern Disko Bugt (Weidick and Bennike, 2007).
Thermal optimum oceanographic conditions (6.2 to ~3.5 ka BP) in the Disko Bugt area
A significant shift in oceanographic conditions is observed from 6.2 ka BP onwards (zone B). This period starts with a prominent peak in productivity, shown by the dominance of N. labradorica, between 6.2 and 5.8 ka BP (Figure 5). This most likely relates to enhanced surface-water productivity leading to increased supply of fresh phytodetritus (food supply) to the sea floor. Such a strong productivity event indicates a shift in the time period, when Arctic sea-ice breakup occurs in spring months at Disko Bugt. We assume that prior to c. 6 ka BP the Arctic sea-ice edge was positioned further south of Disko Bugt and breakup of the sea-ice edge occurred during summer months at Disko Bugt, causing a greater annual bloom over the core site. This is further supported by the occurrence of I. helenae and S. loeblichi between 6.2 and 5.5 ka BP (Figure 5; Polyak and Solheim, 1994). From c. 5.5 ka BP onwards, the distribution of N. labradorica remains at a continuously relatively lower level, suggesting reduced sea-ice edge influence southwest of Disko Bugt.
Geomorpholocial studies in the eastern Disko Bugt area report a largely land-based ice sheet and reduced meltwater runoff from the GIS after c. 6 ka BP (Briner et al., 2010; Weidick and Bennicke, 2007; Weidick et al., 1990). As the ice sheet in the Disko Bugt area was now largely land-based, entrainment of the relatively warm WGC into the embayment could have no direct impact on the ice sheet and force enhanced melting. Consequently, local meltwater discharge from the GIS in the Disko Bugt area will have a more limited influence on the oceanographic record at our core site southwest of Disko Bugt. From this time onwards, the WGC presumably displayed the dominant water mass in Disko Bugt, and it is likely to have influenced surface water properties as well. Henceforth, shifts in abundance of the AtlWcalc and AWagg groups would now indicate changes in the qualitative water mass contributions from the WGCs source currents (EGC and IC, respectively). The overall dominance of the AtlWcalc fauna, between c. 5.5 and 3.5 ka BP (Figure 6d), highlights strong contribution from the IC, and hence a relatively warm and saline WGC in Disko Bugt. A strong and relatively warm IC is also reported from the East Greenland shelf (Figure 6e; Jennings et al., 2002, 2011) and to the south of Iceland during this time interval (e.g. Knudsen et al., 2008b).
Already after 6 ka BP, we observe decreasing current strength of the WGC at the coring site by decreasing sand content to about 15% at 4.0 ka BP (Figure 6c). The reduced flow speed of the WGC registered at the core site does not affect the relative warmth of the WGC, as shown by the dominance and high concentration of the AtlWcalc group at that time (Figure 6d). We assume that the main core of the WGC moved away from the sea floor at the core site. This is accompanied by a slight increase in sedimentation rate, coincident with a increasing TOC (%) content, which is most likely related to enhanced in situ marine productivity at the sea surface (Figures 6a, b).
Between c. 5.5 and 3.5 ka BP, we find low abundance of total AW species (e.g. E. excavatum f. clavata; C. arctica, and S. biformis; Figures 4, 5, 6f), reflecting minimum influence of cool freshwater either from the GIS meltwater and/or the EGC during this period. Accordingly, regional oceanographic conditions had stabilized between c. 5.5 and 3.5 ka BP and the warmest, most ameliorated oceanographic bottom water conditions, reflecting ‘thermal optimum-like’ conditions prevailed. As noted above, a relatively warm WGC was already present in the area from c. 7.3 ka BP onwards, but until c. 5.5 ka BP, the regional oceanographic signal was overprinted by the influence of ice sheet meltwater on the WGC at the site.
The faunal data, presented here, tend to support an interpretation of extended ‘thermal optimum-like’ conditions from c. 7 to 3.5 ka BP. The assumption of a longer optimum would be in accordance with previous reconstructions of thermal optimum conditions from terrestrial and marine studies in the Disko Bugt area (Fredskild, 2000; Funder and Weidick, 1991; Lloyd et al., 2007; Young et al., 2011). Relatively late thermal optimum conditions in the bottom waters southwest of Disko Bugt also support the spatially variable nature of Holocene Thermal Maximum (HTM) conditions, caused by local ice sheet melt influence (see discussion by Kaplan and Wolfe, 2006; Kaufman et al., 2004). Nonetheless, increased occurrence of drift ice on the North Iceland Shelf from 5.5 ka BP onwards (Figure 6g), attributed to an expansion of the EGC, indicates a climatic shift, involving widespread circulation changes in the high latitude North Atlantic (Figure 6g; Moros et al., 2006a). This might consequently also affect the flow strength of the WGC southwest of Disko Bugt.
Onset of Neoglacial cooling (3.5 ka BP to present)
From about 3.5 ka BP onwards (zone C) the benthic foraminiferal assemblage documents cooling/deterioration of the oceanographic/environmental conditions at the core site. We observe cooling/freshening (increase in total abundance of agglutinated specimens and AWagg fauna; Figure 6f) and weakening (decrease in sand content; Figure 6c) of the WGC. This cooling trend is consistent with results from a previous study of nearby coring site 343310 (Figure 1; Perner et al., 2011) and is attributed to an enhanced influence of the EGC. This increased contribution of relatively cold and fresh water masses to the WGC correlates with the onset of Neoglacial cooling in West Greenland.
After 3.5 ka BP, the sedimentation rate increases from an average rate of 0.44 cm/yr to about 0.85 cm/yr, accompanied by a 1% rise in TOC content (Figure 6a, b). This is linked to enhanced in situ productivity at the site, which may also be a consequence of the notably weaker WGC, compared with the preceding intervals (Figure 6c). The observed deterioration of environmental conditions, initiated around 3.5 ka BP, corresponds well with the reported termination of relatively warm conditions on the NW Iceland Shelf (Jiang et al., 2002), in the Denmark Strait (decrease in biogenic carbon content, Figure 6e; Jennings et al., 2011) and in western/south Greenland (e.g. Bennike, 2004; Kaplan et al., 2002; Kuijpers et al., 2003; Lassen et al., 2004; Lloyd et al., 2007; Møller et al., 2006; Moros et al., 2006b; Seidenkrantz et al., 2007, 2008). In addition to this oceanic cooling, terrestrial reconstructions suggest an advance of the GIS after 4 ka BP within the Disko Bugt area (Briner et al., 2010; Weidick and Bennike, 2007; Young et al., 2011), as well as cooling over the GIS (Dahl-Jensen et al., 1998), coinciding with decreasing summer solar insolation (Berger and Loutre, 1991).
Superimposed on this late-Holocene cooling trend we find multicentennial-scale variability in the WGC, correlating with the variability reported by Perner et al. (2011). A cooling around 2.5 ka BP (increased abundance of AW species I. helenae; Figures 5 and S2, available online), is associated with the 2.7 ka BP ‘cooling event’, which has been recorded in various marine and terrestrial records in the North Atlantic region (e.g. Hall et al., 2004; Moros et al., 2004; Oppo et al., 2003; Risebrobakken et al., 2003). A pronounced relatively warm phase is seen at 1.8 ka BP (marked increase in the AtlWcalc fauna; Figure 6d), can be linked to the ‘Roman Warm Period’ (RWP). This period records oceanographic conditions comparable with, or perhaps slightly warmer than, the ‘thermal optimum-like’ conditions seen between c. 5.5 and 3.5 ka BP of West Greenland. Significant warming during the RWP also correlates with findings from the Reykjanes Ridge in the central North Atlantic (Moros et al., 2012).
The gradual cooling of bottom water conditions, which becomes more pronounced from c.1.8 ka BP onwards (gradual rise in AWagg species; Figure 6f), indicates enhanced freshwater forcing from the EGC (cf. Perner et al., 2011) and reduction in IC contribution to the WGC. Nonetheless, a minor warming of bottom water conditions is seen between 1.4 and 0.9 ka BP (Figure 6d), which corresponds to the MCA, suggesting a continuous significant IC contribution to the WGC during this time period. The foraminiferal fauna shows that the MCA warming is less pronounced than that of the RWP, and it is also evident that the RWP records the warmest phase during the late Holocene (see Figure 6d; Perner et al., 2011). Relative cooling of oceanographic conditions from the RWP to the MCA, between 1.8 and 0.7 ka BP, agrees well with previous findings from the Disko Bugt area (Perner et al., 2011) and with other marine records from the northern North Atlantic (e.g. Andrews and Giraudeau, 2003; Moros et al., 2004, 2006a, 2012; Richter et al., 2009).
No sediments are recovered in the present record between 0.7 and 0.3 ka BP (gap in composite record of core 343300). However, severe cooling is seen after 0.9 ka BP culminating in the ‘Little Ice Age’ at 0.3 ka BP related to a continuous increase in the EGC component of the WGC (rise in overall abundance/diversity of agglutinated species and the AWagg fauna, Figure 6f, and supplemented with data from core 343310 (Perner et al., 2011); Figure 6d – light line, 6f – light line). Enhanced EGC contribution to the WGC during the LIA is supported by studies from the East Greenland Shelf, North Icelandic Shelf, Denmark Strait and Southeast and West Greenland, reporting expansion and intensification of the EGC by enhanced contribution of relatively fresher/colder Polar water masses and increased drift ice (e.g. Figure 6g) within the EGC (Andrews et al., 1997; Eirίksson et al., 2004; Jennings et al., 2011; Kuijpers et al., 2003; Moros et al., 2006a; Sha et al., 2011). This oceanic cooling encompasses the reported re-advance of Jakobshavn Isbræ to its LIA maximum position, 20 km west of its current position (Weidick et al., 1990).
Reconstructions by Kaufman et al. (2009) from terrestrial archives document a strong increase in arctic summer air temperature during the last c. 100 years. Contrary to this our data show that oceanographic conditions (WGC) are relatively cool and remained cooler during the last 100 years than during the MCA. Similar results were obtained by Sha et al. (2011) from a site further south on the West Greenland Shelf. These environmental conditions are presumably determined by a relatively strong contribution from the cold EGC during the last 100 years, and are in agreement with the freshening of Baffin Bay between 1916 and 2003, reported by Zweng and Münchow (2006).
Summary and conclusions
A new Holocene benthic foraminiferal record from southwest Disko Bugt provides detailed information of multicentennial-scale bottom water (WGC) variability off West Greenland and interaction of the WGC with the West Greenland ice sheet. Our reconstructions provide a long-term Holocene perspective on the influence of ocean forcing on the Greenland Ice Sheet in the Disko Bugt area.
Between c. 7 and 6.2 ka, we observe a relatively warm and strong WGC, which is likely to support ice sheet melt in Disko Bugt, leading to increased meltwater production, and consequently results in low surface water productivity. A subsequent strong productivity event at c. 6 ka BP suggests that the Baffin Bay Arctic sea-ice edge migrated from its location to the south of Disko Bugt northwards across the site, presumably linked to the influence of the relatively warm and strong WGC. A prolonged relatively warm/stable phase of the WGC is indicated by persisting dominance of the AtlWcalc fauna from c. 5.5 to 3.5 ka BP, reflecting ‘thermal optimum-like’ conditions off West Greenland.
Most likely a relatively warm and strong WGC was continuously present on the shelf of Disko Bugt during the entire period between c. 7 and 3.5 ka BP, albeit its signal is diluted/deflected (overprinted) in our data by the melting ice sheet, which in turn resulted in relatively cooler and variable environmental conditions.
From about 3.5 ka BP, benthic foraminifera identify a long-term late-Holocene cooling of the bottom waters, associated with the onset of Neoglacial cooling. This indication of a long-term cooling agrees well with studies from West/East Greenland fjord and shelf areas, which report gradual cooling of the WGC, along with glacial re-advances. It is likely that cooling of oceanographic conditions favoured the observed re-advance of the ice sheet in the Disko Bugt area after c. 3.5 ka BP. Superimposed on this cooling trend, we reconstruct marked multicentennial-scale variability within the WGC: (1) a cooling at c. 2.5 ka BP, related to the 2.7 ka BP ‘cooling event’; (2) a relatively warm phase at c. 1.8 ka BP, corresponding to the ‘Roman Warm Period’; (3) only a slight warming in bottom waters at the transition into the ‘Medieval Climate Anomaly’ and; (4) strong cooling from c. 1.8 ka BP culminating in the ‘Little Ice Age’ cold period.
Cooling of bottom waters, confirmed by a gradual rise in AWagg species, is linked to an enhanced influence of fresher/cooler water mass contribution from the EGC to the WGC. Agglutinated species dominate the benthic foraminiferal assemblage also during the last 100 years, corroborating a persistent strong influence of the EGC on the WGC, consistent with previous studies from Disko Bugt and the West Greenland shelf.
Footnotes
Appendix
List of identified foraminiferal taxa.
| Foraminiferal list |
|---|
| Adercotryma glomerata (Brady, 1878) |
| Ammoscalaria pseudospiralis (Williamson, 1858) |
| Astrononion gallowayi Loeblich & Tappan, 1953 |
| Bolivina pseudopuncata Höglund, 1947 |
| Buccella frigida (Cushman, 1922) |
| Buccella frigida calida Cushman & Cole, 1930 |
| Cassidulina neoteretis Seidenkrantz, 1995 |
| Cassidulina reniforme Nørvang, 1945 |
| Cibicides lobatulus (Walker & Jacob, 1798) |
| Cribrostomoides crassimargo (Norman, 1858) |
| Cribrostomoides jeffreysi (Williamson, 1958) |
| Cribrostomoides spp. |
| Cuneata arctica (Brady, 1881) |
| Dentalina spp. |
| Deuterammina ochracea (Williamson, 1858) |
| Discorbina spp. |
| Eggerella advena (Cushman, 1922) |
| Elphidium albiumbilicatum (Weiss, 1954) |
| Elphidium bartletti Cushman, 1933 |
| Elphidium hallandese Brotzen, 1943 |
| Elphidium excavatum (Terquem) forma clavata Cushman, 1930 |
| Elphidium spp. |
| Epistominella vitrea Parker, 1952 |
| Fissurina spp. |
| Glandulina spp. |
| Globobulimina auriculata arctica (Höglund, 1954) |
| Islandiella helenae Feyling-Hanssen & Buzas, 1976 |
| Islandiella norcrossi (Nørvang, 1945) |
| Lenticulina spp. |
| Melonis barleeanus (Williamson, 1858) |
| Nonionella turgida var. digitata Nørvang, 1945 |
| Nonionella auricula Heron-Allen & Earland, 1930 |
| Nonionellina labradorica (Dawson, 1860) |
| Parafissurina spp. |
| Portatrochammina spp. |
| Pullenia bulloides (d'Orbigny, 1826) |
| Pullenia osloensis Feyling-Hanssen, 1954 |
| Procerolagena spp. |
| Quinqueloculina spp. |
| Recurvoides turbinatus (Brady, 1881) |
| Reophax fusiformis (Williamson, 1858) |
| Reophax gracilis (Kiear, 1900) |
| Reophax pilulifer Bandy, 1884 |
| Robertina arctica d'Orbigny, 1846 |
| Saccammina difflugiformis (Brady, 1879) |
| Spiroplectammina biformis (Parker & Jones, 1865) |
| Stainforthia feylingi Knudsen & Seidenkrantz, 1994 |
| Stainforthia loeblichi (Feyling-Hanssen, 1964) |
| Stainforthia spp. |
| Textularia earlandi Phleger, 1952 |
| Textularia torquata Parker, 1952 |
| Trifarina fluens (Todd, 1947) |
Acknowledgements
We thank the Captain and Crew of the R/V Maria S. Merian for their excellent work during cruise MSM05/03. We also thank Marit-Solveig Seidenkrantz for fruitful discussion of the benthic foraminiferal assemblage and their interpretation, Tomasz Goslar from Poznań Radiocarbon Laboratory and also Richard Telford from the University of Bergen, for performing the age–depth modelling of the gravity core.
Funding
The authors thank the Deutsche Forschungsgemeinschaft (DFG) for funding the project ‘Disko Climate’ (MO1422/2-1).
