Journal cover Journal topic
Biogeosciences An interactive open-access journal of the European Geosciences Union
Journal topic
Biogeosciences, 15, 5909-5928, 2018
https://doi.org/10.5194/bg-15-5909-2018
Biogeosciences, 15, 5909-5928, 2018
https://doi.org/10.5194/bg-15-5909-2018

Research article 08 Oct 2018

Research article | 08 Oct 2018

# Coastal primary productivity changes over the last millennium: a case study from the Skagerrak (North Sea)

Coastal primary productivity changes over the last millennium
Anna Binczewska1, Bjørg Risebrobakken2, Irina Polovodova Asteman2,a, Matthias Moros3, Amandine Tisserand2, Eystein Jansen2,4, and Andrzej Witkowski1 Anna Binczewska et al.
• 1Faculty of Geosciences, University of Szczecin, Szczecin, Poland
• 2Uni Research Climate, Bjerknes Centre for Climate Research, Bergen, Norway
• 3Leibniz Institute for Baltic Sea Research (IOW), Warnemünde, Germany
• 4Department of Earth Science, University of Bergen, Norway
• acurrently at: Marin Mätteknik (MMT) Sweden AB, Gothenburg, Sweden
Abstract

A comprehensive multi-proxy study on two sediment cores from the western and central Skagerrak was performed in order to detect the variability and causes of marine primary productivity changes in the investigated region over the last 1100 years. The cores were dated by Hg pollution records and AMS 14C dating and analysed for palaeoproductivity proxies such as total organic carbon, δ13C, total planktonic foraminifera, benthic foraminifera (total assemblages as well as abundance of Brizalina skagerrakensis and other palaeoproductivity taxa) and palaeothermometers such as Mg∕Ca and δ18O. Our results reveal two periods with changes in productivity in the Skagerrak region: (i) a moderate productivity at  CE 900–1700 and (ii) a high productivity at  CE 1700–present. During  CE 900–1700, moderate productivity was likely driven by the nutrients transported with the warm Atlantic water inflow associated with a tendency for a persistent positive NAO phase during the warm climate of the Medieval Climate Anomaly, which continues into the LIA until  CE 1450. The following lower and more variable temperature period at  CE 1450–1700 was likely caused by a reduced contribution of warm Atlantic water, but stronger deep-water renewal, due to a generally more negative NAO phase and a shift to the more variable and generally cooler climate conditions of the Little Ice Age. The productivity and fluxes of organic matter to the seafloor did not correspond to the temperature and salinity changes recorded in the benthic Melonis barleeanus shells. For the period from  CE 1700 to the present day, our data point to an increased nutrient content in the Skagerrak waters. This increased nutrient content was likely caused by enhanced inflow of warm Atlantic water, increased Baltic outflow, intensified river runoff, and enhanced human impact through agricultural expansion and industrial development. Intensified human impact likely increased nutrient transport to the Skagerrak and caused changes in the oceanic carbon isotope budget, known as the Suess effect, which is clearly visible in our records as a negative shift in δ13C values from  CE 1800. In addition, a high appearance of S. fusiformis during the last 70 years at both studied locations suggests increased decaying organic matter at the sea floor after episodes of enhanced primary production.

1 Introduction

Growth of marine microalgae is stimulated by enhanced concentrations of major nutrients like nitrogen and phosphorus in the photic zone (e.g. Sigman and Hain, 2012). Microalgae in the oceans are primary producers, which provide food for consumers at higher trophic levels and oxygen for respiration (e.g. Micheli, 1999). Through photosynthesis and the biological pump, marine primary producers also extract CO2 from the atmosphere. Carbon, nutrients and trace elements, which are fixed by primary producers, are further ingested by higher organisms or sink to the ocean floor where they are stored in the form of organic matter. The organic matter will eventually be remineralised, releasing the carbon, nutrients and trace elements to the bottom water (e.g. Sigman and Hain, 2012). A fraction is also buried in sediments. Supply of carbon to the oceanic ecosystems is an important part of the biogeochemical cycles, which are presently being disturbed by the human impact. Via photosynthesis, primary producers can help to remove CO2 from the atmosphere through the biological pump. However, increasing levels of dissolved CO2 in the oceans associated with so-called “ocean acidification” may negatively impact ocean's carbonate producers (e.g. Doney et al., 2009; Haynert et al., 2012) by influencing their survival and fitness (e.g. Thomsen et al., 2017). Excessive export of organic matter furthermore changes the oxygen condition at the sea floor due to decay of organic matter, which lowers the dissolved oxygen content, in turn negatively influencing the benthic life regime (e.g. Kristiansen and Aas, 2015, and references therein).

Coastal zones are among the most productive marine regions, characterised by high atmospheric CO2 uptake, organic matter accumulation and decomposition (e.g. Hjalmarsson et al., 2010). The Skagerrak, located between the North Sea and the Baltic Sea and in the close proximity to land, has many potential nutrient sources, such as the North Atlantic, Baltic Sea, and North Sea, as well as continental discharge and river runoff (Aure and Dahl, 1994; Andersson, 1996; Gustafsson and Stigebrandt, 1996). Upwelling and precipitation further increase the nutrient supply to the surface waters, additionally stimulating productivity in this region (Pingree et al., 1982; Aure and Dahl, 1994; Fonselius, 1996). The North Sea and the Skagerrak absorb significant quantities of atmospheric CO2 via the biological pump (1.38 and 1.2 mol C m−2 yr−1, respectively) and thus play an important role in the carbon cycle (Thomas et al., 2005; Hjalmarsson et al., 2010). The Skagerrak acts as a main depositional basin for about half of the refractory carbon produced in the North Sea and for a high amount of labile organic matter either imported with waters of the near-bottom current from the Danish coast or produced by intense algal blooms (Boon et al., 1999). The input of nutrients largely regulates food webs, which makes nutrients of great economic importance for the coastal areas worldwide (Micheli, 1999; FAO, 2016). In the Skagerrak nutrients are particularly important for the Nordic fisheries (Hop et al., 1992; Iversen et al., 2002; Olsen et al., 2004; Skogen et al., 2007). Fisheries and aquaculture sectors of the Skagerrak, commercially valuable for the Scandinavian nations, also make a relevant contribution to a growing global demand for food, which will largely rely on coastal regions to host a major part of food production in the future (e.g. FAO, 2016). Effects of increased primary production range from positive impacts on growth rate, size and reproduction of fish and shellfish populations to disruptive alterations in the food webs, thus yielding or reducing the profit rates of fisheries (Hop et al., 1992; Micheli, 1999; Iversen et al., 2002; Olsen et al., 2004, 2005; Breitburg et al., 2009; FAO, 2016). Disruptive changes in trophic levels of the Skagerrak ecosystem have been attributed to overfishing (Cardinale and Svedäng, 2004) but ongoing studies alert about adverse impact of increased nutrient inputs driving heavy phytoplankton blooms and eutrophication in the region (e.g. Baden et al., 1990; Aure et al., 1996; Breitburg et al., 2009). Eutrophication causes high demand and depletion of oxygen in the bottom waters, which affect species diversity, morphology and population growth and forces organisms to migrate (Rosenberg et al., 1990; Conley, 2009). Thus, to better understand the ongoing and possible future productivity changes and associated environmental effects, more historical studies are needed.

Although the Skagerrak is a well-investigated area, retrospective studies from this region have focused on climate change instead of productivity and, hence, past primary productivity changes are not well known (Polovodova Asteman et al., 2018). Previous studies suggested that among the processes driving primary productivity changes in the Skagerrak are (1) the North Atlantic Oscillation (NAO), (2) an input of anthropogenic nutrients from the Skagerrak catchment area and (3) nutrients transported from the Baltic Sea. Thus, while the NAO influences the inflow of the nutrient-rich Atlantic water mass into the Skagerrak (Gil et al., 2005; Brückner and Mackensen 2008), anthropogenic activities in the Skagerrak catchment area, such as land-use changes, also cause an increased nutrient input and high organic matter flux to the basin (e.g. Filipsson and Nordberg, 2010). In addition, eutrophication of the Baltic Sea and the resulting transport of nutrient-rich water with the Baltic Current through surface water exchange processes potentially represents an additional nutrient source (Andersson, 1996; Hjalmarsson et al., 2010; Krossa et al., 2015; Filipsson et al., 2017; Polovodova Asteman et al., 2018). In its turn, a nutrient-overloaded ecosystem in the Skagerrak may further supply dissolved inorganic nitrogen (DIN) to the Kattegat bottom waters (Carstensen et al., 2006) and increase the nutrient level in the coastal waters of Norway (Rosenberg et al., 1987; Rydberg et al., 2006) (Fig. 1). The last two processes involve human-induced influences on the ecosystem, while the NAO relates to effects following variability in the climate system. Over the last 1100 years the interaction between these processes has changed, as human influence has increased. In perspective of the last 4500 years, it is seen that productivity and its variability has increased from 50 years before the Common Era (BCE) until the present day (Polovodova Asteman et al., 2018). The increased productivity occurred as a response to enhanced local runoff coinciding with high winter rainfall and general cooling in Scandinavia, as well as intensified Baltic outflow, which may have significantly contributed to the nutrient supply in the Skagerrak in the past (Polovodova Asteman et al., 2018).

Figure 1Map of core locations (stars) in the Skagerrak in the north-eastern North Sea (modified from Polovodova Asteman et al., 2018). The Norwegian Trench is outlined by a thin grey dotted curve along the coast of Norway. Major current systems and water masses in the Skagerrak are indicated by arrows: the Baltic Current (BC), Northern Jutland Current (NJC), Southern Jutland Current (SJC), Norwegian Coastal Current (NCC), North Sea Water (NSW) and Atlantic Water (AW). The light grey arrows show surface water and the black arrows show subsurface water.

Our study aims to detect the variability and causes of marine primary productivity changes in the Skagerrak over the last 1100 years. We address this aim through a comprehensive multi-proxy study of two sediment cores from the central and the western Skagerrak, integrating records of the total organic carbon (TOC), foraminiferal assemblage data, stable carbon and oxygen isotopes (δ13C, δ18O) and trace element ratio (Mg∕Ca). In addition, we evaluate our results in the context of the long-term productivity changes described previously in Polovodova Asteman et al. (2018).

2 Study area

The Skagerrak basin is located in the north-eastern part of the North Sea, connected to the Baltic Sea through the Kattegat (Fig. 1). The basin has a mean water depth of 210 m and a sill depth of 270 m. With a maximum depth of 700 m, the Skagerrak represents the deepest part of the Norwegian Trench (Rodhe, 1996). The area is characterised by an anticlockwise circulation and complex hydrography. The surface circulation (< 30 m) is to a large extent dominated by a surface current consisting of inflowing saline water from the southern North Sea and the North Atlantic and outflowing less saline water from the Baltic Sea (Danielssen et al., 1997). The inflowing nutrient-rich surface water flows along the Danish coast driven by the Southern Jutland Current (SJC) and the Northern Jutland Current (NJC) while the outflowing Baltic Sea water (BW ∼20–30 psu) flows as the Baltic Current (BC) along the Swedish west coast towards the north-eastern Skagerrak where it merges with the NJC and turns to the north-west as the low-salinity Norwegian Coastal Current (NCC) (Rodhe, 1996; Rydberg et al., 1996). The water flowing from the Skagerrak towards the Norwegian Sea (NCC) partly recirculates to the western Skagerrak (Rodhe, 1996). The surface water has a high nutrient concentration mostly due to the freshwater input via rivers draining from the Norwegian south coast, German and Danish east coasts, and the Baltic catchment area, but the upwelling of the underlying nutrient-rich Atlantic water is also considered to be an additional nutrient supply (Gustafsson and Stigebrandt, 1996; Rodhe, 1996). As a consequence of the mixing of different water types and the high freshwater input enhanced by precipitation, the upper layer of the surface water has low salinity (25–32 psu) and is determined as the Skagerrak Coastal Water (4.5–10 C). The intermediate water layer (30–270 m) is referred to as Skagerrak Water (32–35 psu, 4.5–10 C) and is driven by the subsurface circulation (Andersson, 1996). The deep-water layer below sill depth (> 270 m) is dominated by Atlantic Water and is recognised as the Skagerrak Basin Water (> 35 psu, 5.5–6.5 C) (Aure and Dahl, 1994).

The subsurface circulation (below 30 m water depth) consists of nutrient-rich Atlantic deep water (AW > 35 psu, 5.5–8.5 C) flowing through the northern North Sea, and the water from the central and southern North Sea (NSW ∼31–35 psu) (Rodhe, 1996). The inflowing water follows the southern side of the Norwegian Trench, entering the Skagerrak in its central part where this water is mixed with fresh surface water and flows out as the NCC (Winther and Johannessen, 2006).

Large-scale atmospheric systems and regional meteorological factors (e.g. precipitation and storms) influence the flow regime, creating a highly dynamic system in the upper layer of the water column where water mixing is largely caused by the south-westerly winds (Gustafsson and Stigebrandt, 1996). At the same time, calmer hydrographic conditions are typical for the intermediate layer and the deep water down to ∼400 m, with a maximum water residence time of 3 months (Andersson, 1996). This is in contrast to the renewal of the deepest water mass below ∼400 m, which is replenished every 1 to 3 years depending on the strength of the Atlantic water inflows (Aure and Dahl, 1994; Rodhe 1996), closely correlating with the NAO index (Brückner and Mackensen, 2006).

3 Material and methods

Two gravity cores (GCs) were retrieved from the Skagerrak during a R/V Elisabeth Mann-Borgese cruise in May 2013. Core EMB046/20-3GC (4.8 m long) was taken from the central Skagerrak (south-east Norwegian Trench; 5831.75 N, 0929.13 E; 533 m water depth), while core EMB046/10-4GC (4.62 m) comes from the western Skagerrak (south-west Norwegian Trench; 5749.73 N, 0717.62 E; 457 m water depth) (Fig. 1). Both cores were taken below sill depth within the deep waters of the Norwegian Trench. The cores were cut in 1 m sections on board before being split and subsampled every 1 cm ashore. This study is based on results from the upper 170.5 and 164.5 cm of cores EMB046/20-3GC and EMB046/10-4GC respectively, which corresponds to the last 1100 years. CTD measurements are available in the Supplement (Fig. S1). A part of the foraminiferal dataset and the TOC data were previously published in Polovodova Asteman et al. (2018). Here we present new stable isotopes (δ13C, δ18O), and trace element ratio (Mg∕Ca) data covering the last 1100 years, in combination with foraminiferal assemblage data and multivariate statistics. Both cores consist of mostly homogeneous soft organic-rich clay, have olive-grey colour and show no significant changes in grain size and lithology throughout the studied intervals. The TOC content was determined using “Rapid CS cube –Elementar” analyser (Department of Geology and Paleogeography, University of Szczecin, Poland) with a measurement accuracy of 0.01 % (95 % confidence level (CL) at 99.5 % detection limit (DL)). For detailed methodology of geochemistry measurements (TOC) see Polovodova Asteman et al. (2018).

For the stable carbon and oxygen isotopes (δ13C, δ18O) and trace element analyses, well-preserved shells of benthic foraminifera Melonis barleeanus were picked from the dried sediment (fraction > 150 µm). Melonis barleeanus was selected for the analyses due to its relatively high abundance throughout the investigated intervals at both sites and its well-known potential for geochemical palaeoreconstructions (Mackensen et al., 2000; Kristjánsdóttir et al., 2007; Brückner and Mackensen, 2008; Butruille et al., 2017). Stable isotope measurements were performed at 1 cm intervals and the trace element analyses were done at 1–2 cm intervals from both gravity cores down to 170 cm (EMB046/20-3GC) and 164.5 cm (EMB046/10-4GC). No stable isotope or trace element analyses were done between 4.5 and 7.5 cm in EMB046/10-4GC, due to lack of material, because most of the foraminifera from that interval where used for 14C AMS dates.

Stable isotope analyses were run on a Finnigan MAT 253 mass spectrometer equipped with an automatic Kiel device at FARLAB of the University of Bergen. Prior to measurement the tests of M. barleeanus were lightly crushed, cleaned with methanol (≥98.8 %) using an ultrasonic bath and dried at 60 C. For each measurement 2–4 specimens were used. All results are reported in ‰ versus Vienna Pee Dee Belemnite standard (V-PDB), using the National Bureau of Standards (NBS) 19 and 18, in combination with the internal lab standard CM12. The long-term analytical uncertainty is ±0.08 ‰ and ±0.03 ‰ (95 % CL) for oxygen and carbon isotopes, respectively.

Shells of M. barleeanus were also cleaned and analysed for trace elements at the Trace Element Lab (TELab) at Uni Research Climate, Bergen (Norway). For each analysis, approximately 15–20 specimens were gently crushed between two glass plates under a microscope to allow the contaminants to be removed. The samples were cleaned following the procedure described by Barker et al. (2003). The cleaning method includes clay removal steps, oxidation of the organic matter and surface leaching. Samples containing enough material were mixed and split into two subsamples to allow duplicate analysis. All samples were dissolved in trace metal pure 0.1 M HNO3 (prepared from HNO3 TraceSELECT®) and diluted to a final concentration of 40 ppm of calcium (Ca).

The trace elements were measured on an Agilent 720 inductively coupled plasma optical emission spectrometer (ICP-OES). Six standards have been prepared at the TELab and they have a composition similar to foraminiferal carbonate (0.50–7.66 mmol mol−1). Fe∕Ca and Al∕Ca values have been checked and showed no significant correlation with the measured Mg∕Ca values. The correlation coefficients (R2) between the Mg∕Ca ratio and the Fe∕Ca and Al∕Ca ratios are 0.012 and 0.095 for the core EMB046/20-3GC and 0.020 and 0.067 for the core EMB046/10-4GC, indicating no systematic contamination due to insufficient cleaning. The Mn∕Ca values in our samples are higher than the recommended maximum (< 105 µmol mol−1) (Boyle, 1983), indicating that diagenetic coatings might also affect our results. The Mn∕Ca values, however, show no significant correlation with the measured Mg∕Ca values (R2=0 for EMB046/20-3GC and R2=0.021 for EMB046/10-4GC). A standard solution with Mg∕Ca of 5.076 mmol mol−1 is analysed after every eight samples to correct for instrumental biases and analytical drift of the instrument. The long-term Mg∕Ca analytical precision, based on the standard solution, is ±0.016 mmol mol−1 (1σ standard deviation, SD) or 3.11 % (relative SD). The average reproducibility of duplicate measurements (pooled SD, dof = 41) is equivalent to an overall average precision of 4.09 %. The average Mg∕Ca of long-term international limestone standard (ECRM752-1) measurements was 3.76 mmol mol−1 (1σ= 0.07 mmol mol−1) with the average published value of 3.75 mmol mol−1 (Greaves et al., 2008).

The Melonis spp. Mg∕Ca bottom water temperature (BWT) is calculated from the measured Mg∕Ca using the new Melonis spp. calibration (Mg∕Ca= 0.113 (±0.005)*BWT + 0.792 (±0.036)), based on core top data covering a Mg∕Ca range of 0.68–3.66 mmol mol−1 for a temperature range of 0.89–15.58C (Hasenfratz et al., 2017). According to the calibration uncertainty, a 1σ temperature error (95 % CL) of ±0.9 to ±1.7C has to be taken into consideration for the temperature range (4.1–9 C) covered by EMB046/20-3GC. For EMB046/10-4GC the 1σ temperature error (95 % CL) is similar and ranges between ±1 and ±1.6C for the temperatures 5–8.5 C. Further discussion on Mg∕Ca-derived BWT within this article will utilise the Hasenfratz's calibration (Hasenfratz et al., 2017).

The results from Mg∕Ca-derived BWT and δ18O were used to calculate the seawater isotopic composition (δw) using the following equation: δw= ((δ18Otopδ18Odown)–(TtopTdown)*0.23) + 0.3. The δ18Otop is the δ18O value measured on the foraminiferal shells from the uppermost sample (0–1 cm core depth) and the Ttop is the temperature taken from CTD measured at the time of coring. The δ18Otop equals 2.04 and 2.06 (‰) while Ttop was 5.72 and 5.34 C for the EMB046/10-4GC and EMB046/20-3GC, respectively. The δ18Odown and Tdown is the down core δ18O and temperature as measured on the foraminifera samples. Temperature estimates based on δ18O follow Shackleton (1974), where 0.23 ‰ equals 1 C in the temperature interval estimated for these sites. A constant of 0.3 ‰ is used to correct for the difference between V-PDB and δw. No correction for an ice volume effect was applied as this is considered negligible over the last 1000 years.

The foraminiferal analysis was carried out on 5–10 g wet sediment, gently sieved over a 63 µm sieve and wet-counted for foraminifera immediately afterwards. The two Skagerrak records were counted, over the targeted time interval covering the last 1100 years, at 1–3 cm resolution, with the exception of a 7 cm interval between 89 cm and 96 cm ( CE 1030–1055) for EMB046/20-3GC, and at 1–2 cm resolution for the EMB046/10-4GC record. In the > 63 µm fraction, at least 300 benthic and 300 planktonic (where possible) specimens were counted under a stereomicroscope and identified to a species level. Both relative (%) and absolute (individuals per gram in wet sediments, denoted as “ind. g−1 wet sed.”) abundances were calculated. The benthic foraminiferal species were categorised depending on their relative abundance in the assemblage as dominant (> 10 %), accessory (5 %–10 %) or rare (< 5 %). Only dominant and accessory species are discussed in this study (Fig. 3, Table 2). Benthic species with a relative abundance of > 5 % in at least two samples were subject to multivariate statistics using simple CABFAC factor analysis with varimax rotation (Table 2), performed by the PAST software (Hammer et al., 2001). This statistical tool provides a reliable method to distinguish the statistically significant foraminiferal units dominated by different species (e.g. Polovodova Asteman et al., 2013). In addition, benthic foraminiferal species indicative of increased organic matter fluxes to the sea floor, and hence, algal blooms, were grouped as “palaeoproductivity fauna” and included Alabaminella weddelensis, Brizalina skagerrakensis, Bulimina marginata, Epistominella spp., Nonionella iridea and Uvigerina spp. (Polovodova Asteman et al., 2018, and references therein). The planktonic foraminifera are presented as total planktonic individuals.

4 Chronology

A common age model has previously been established for the two cores (Polovodova Asteman et al., 2018). The two age models were set at a common depth scale based on 30 available AMS 14C dates, as well as a correlation among total inorganic carbon, relative abundance of B. skagerrakensis and mercury (Hg) records from both cores. All dates were calibrated using Calib 7.10 (Stuiver et al., 2017), Marine13 (Reimer et al., 2013) and $\mathrm{\Delta }R=\mathrm{0}±\mathrm{50}$. The age model of Polovodova Asteman et al. (2018) is reasonable when investigating the longer-term trends over the last 4.5 ka. However, when focusing on the last 1100 years we found that there was a need for an improvement of the age model over this time period. This was achieved by increasing the number of dates and by a fine-tuning of the reservoir age. The new age model of EMB046/20-3GC, reaching back to CE 295, is based on eight 14C AMS dates in addition to the initial Hg increase at CE 1900 (Moros et al., 2017). The core top was set to 2013, the year of coring. A modern core top age is confirmed by a post-bomb age at 5.5 cm, as well as a recording of the Cs-137 signal associated with the nuclear weapons testing period (not shown). Considering all of the information from the upper part of the core in detail makes it clear that the use of a $\mathrm{\Delta }R=\mathrm{200}±\mathrm{50}$ provides a better transfer, avoiding an unlikely jump in sedimentation rate, between the modern ages and the 14C ages than when using a reservoir age of $\mathrm{\Delta }R=\mathrm{0}±\mathrm{50}$. Hence, when establishing the new age model, the 14C AMS dates were calibrated using Calib 7.10 (Stuiver et al., 2017), Marine13 (Reimer et al., 2013) and $\mathrm{\Delta }R=\mathrm{200}±\mathrm{50}$, as well as linear interpolation between the established tie points (Table 1; Fig. 2). The new ages calculated for EMB046/10-4GC are based on the new age model of EMB046/20-3GC and the previously established common depth scale for EMB046/10-4GC and EMB046/20-3GC (Polovodova Asteman et al., 2018) (Fig. 2). Due to the established relationship between the depth scales of the two cores (Polovodova Asteman et al., 2018), the age model for EMB046/20-3GC can be, and are, used to create the new age model for EMB046/10-4GC.

Table 1Information about the new chronology of core EMB046-20-3GC over the last ca. 1600 years

5 Results

## 5.1 Organic geochemistry of bulk sediment

Both records show low TOC values until  CE 1700, around 1.7 %–2.1 % and 1.5 %–1.8 % in EMB046/20-3GC and EMB046/10-4GC, respectively (Fig. 5). From  CE 1700, the TOC content strongly increases towards the core tops, where the ranges 1.85 %–2.5 % (EMB046/20-3GC) and 1.75–2.3 (EMB046/10-4GC) are recorded. When comparing the two cores, the TOC values are higher for the EMB046/20-3GC record than for EMB046/10-3GC (Fig. 5).

Figure 2(a) Linear interpolation between the established tie points for the EMB046/20-3GC age model. (b) The relationship between B. skagerrakensis of EMB046/10-4GC (red curve) and EMB046/20-3GC (blue curve).

## 5.2 Carbon isotopes

Both δ13C records show similar long-term variations through the study interval. Between CE 1500 and 1700 there is, however, a distinct increase in δ13C values for EMB046/20-3GC and a decrease for EMB046/10-4GC. Mean δ13C values of the time interval between  CE 900 and 1700 are generally higher (0.53 ‰ on average in EMB046/10-4GC and 0.44 ‰ on average in EMB046/20-3GC) than during  CE 1700–2000 when mean δ13C values of 0.73 ‰ (EMB046/10-4GC) and 0.58 ‰ (EMB046/20-3GC) are observed (Fig. 5). From  CE 1700 towards the present, both records show a strong decreasing δ13C trend from ca. 0.2 ‰ to 1.6 ‰, where generally lower absolute δ13C values are recorded in EMB046/10-4GC than in EMB046/20-3GC (Fig. 5).

## 5.3 Oxygen isotopes

Both Skagerrak records display similar δ18O values, ranging from around 1.7 ‰ to 2.7 ‰ (Fig. 4). In general, the δ18O in EMB046/20-3GC shows lower values at  CE 1050–1350 and somewhat higher or more variable values between  CE 1400 and 1550, followed by a decrease until  CE 1700. The overall lower δ18O values in EMB046/10-4GC between  CE 900 and 1550, characterised by the lowest recorded δ18O at  CE 1350–1550, is interrupted by an increase at  CE 1200–1350. Consequently, the δ18O long-term trend is not common for both records; however, there is one distinct period in both records of relatively high δ18O values between  CE 1700 and 1800, after which the δ18O gradually decrease until  CE 1850, followed by a steady increased δ18O in EMB046/20-3CG and more variable but overall lower δ18O values in EMB046/10-4GC towards the core top (Figs. 4–5).

Figure 3Foraminiferal assemblages including dominant and accessory benthic species for both cores (EMB046/10-4GC and EMB046/20-3GC) and CABFAC results. Absolute abundance is shown as grey shading, while relative abundance as black curve with symbols. The absolute abundance of total benthic foraminifera is a sum of all species: agglutinated (Agglut.) and calcareous (Calc.). The dashed line divides record into two periods of the most pronounced palaeoproductivity changes, which are discussed in the text.

## 5.4 Mg∕Ca analyses and BWT

The Mg∕Ca values vary from 1.33 to 1.87 mmol mol−1 in the EMB046/20-3GC record and from 1.37 to 1.97 in EMB046/10-4GC, in general giving an estimated BWT range between 4.7 and 8.1 C, which is within the range of instrumentally recorded temperatures (ICES, 2010; Fig. 4b). This comparison and further interpretation refer to the smoothed data, while the raw data are mostly within the range of the instrumental data, but not completely. Through the records, there is in general a good correlation between the BWT changes and the variability of the δ18O values, showing similar patterns with periods of higher BWT corresponding to those with decreased δ18O values and vice versa, except between  CE 1000 and 1150 in EMB046/20-3GC where both proxies show relatively high values (Fig. 4). Hence, there is consistency between the proxies within each core but not between the cores (Figs. 4–5). Overall, BWT in EMB046/10-4GC is characterised by relatively little variability between  CE 900 and 1550 followed by higher variability in the records. In contrast, BWT in EMB046/20-3GC first shows a decreasing trend until  CE 1500, somewhat higher but variable values between  CE 1550 and  CE 1700, and a drop in values at  CE 1700–1800, and again warmer BWT is shown for the youngest part of the record ( CE 1800–2000) (Figs. 4–5).

## 5.5 Water isotopic composition

The changes in δw records follow the pattern in the BWTs curves in both records. The decrease in δw values corresponds to periods of low temperature and vice versa. Lower correlation between δw and BWT was observed in the early part of the EMB046/20-3GC record ( CE 1000–1150) (Figs. 4–5).

## 5.6 Foraminiferal assemblages

Eight planktonic foraminiferal species are identified in both records. Globigerinita uvula and G. glutinata are the most abundant species, while Globorotalia inflata, Globigerina bulloides, Neogloboquadrina pachyderma, Neogloboquadrina incompta, Turborotalia quinqueloba and Orbulina universa are less abundant. However, in this study all planktonic species are presented together as total individuals per gram of wet sediments (Fig. 3) due to their overall low absolute abundance, varying between 0 and 49 ind. g−1 wet sed. (EMB046/20-3GC) and 1.6–110 ind. g−1 wet sed. (EMB046/10-4GC). The planktonic foraminifera are most abundant in the interval between  CE 900 and 1700 in EMB046/10-4GC and at  CE 900–1550 in EMB046/20-3GC, after which they decrease towards the top of the cores and almost disappear after  CE 1850 in both records (Figs. 3, 5).

Table 2List of dominant (bold) and accessory benthic foraminiferal species. The species names marked by “*” represent foraminiferal species with relative abundance > 5 % in only one sample and thus were excluded from statistic classification.

Table 3The factor results from a CAB-FAC factor analyses.

Table 4The varimax scores for factors 1–3. The bold numbers indicate foraminiferal species with absolute value of factor scores > 1.

Figure 4(a) Stable oxygen isotopes (δ18O), (b) Mg∕Ca-derived bottom water temperature (Mg∕Ca-derived BWT), and (c) stable oxygen isotope composition of seawater (δw) from both studied cores against age. The errors bands represent 1 SD uncertainties of the records. (d) The same parameters as in (a–c) (δ18O, Mg∕Ca-derived BWT, δw); however, only measurements from core depths where all proxies are available are shown (symbols). The curves correspond to a five-point running average. The light blue box indicates the range of instrumentally recorded temperatures of the time period between 2009 and 1924 years from the area between 57 17–58 N and 8–9 79 E at 300–340 dbar (ICES 2010).

The benthic foraminiferal record from core EMB046/10-4GC is characterised by consistently high absolute abundances (123–455 ind. g−1 wet sed.) in contrast to overall lower values in core EMB046/20-3GC (43–527 ind. g−1 wet sed., where the highest value represents an individual peak above 361 ind. g−1 wet sed. significantly standing out from the rest of the record). In EMB046/20-3GC, the absolute abundance of benthic foraminifera is high until  CE 1400. At CE 1400, then values drop below 150 ind. g−1 wet sed. for the next ∼300 years. After  CE 1700 the absolute abundances gradually increase to reach the highest recorded values between CE 1800 and 1900. A similar absolute abundance trend is shown for agglutinated foraminifera; however, those appear in higher numbers in core EMB046/10-4GC than in EMB046/20-3GC (Fig. 3).

The benthic foraminiferal assemblages consist of up to 61 and 57 species in the cores EMB046/20-3GC and EMB046/10-4GC, respectively. Among those, eight species are dominant (> 10 %) and nine are accessory (5 %–10 %) in core EMB046/20-3GC while the EMB046/10-4GC record has six dominant and five accessory taxa (Table 2). The common dominant species for both cores include Brizalina skagerrakensis, Cassidulina laevigata, Eggereloides medius, Nonionella iridea, Pullenia osloensis and Stainforthia fusiformis (for a full list of dominant and accessory species see Table 2). Among benthic foraminiferal species Brizalina skagerrakensis shows the most prominent and consistent changes when comparing both records (Fig. 3).

The CABFAC factor analysis distinguished three factors for each of the cores (Fig. 3), which together explain 95 % (EMB046/20-3GC) and 97 % (EMB046/10-4GC) of the total variance (Table 3). The foraminiferal species with absolute value of factor scores > 1 are considered to contribute significantly to the defined foraminiferal assemblages (Table 4) and are used to name the distinguished factors (assemblages). Thus, “Pullenia osloensis assemblage” associated with Factor 1 explains 81 % (EMB046/20-3GC) and 86 % (EMB046/10-4GC) of the variance, and includes species P. osloensis and Nonionella iridea, defined for both records, with an addition of Cassidulina laevigata in the EMB046/10-4GC dataset. The “Brizalina skagerrakensis assemblage” associated with Factor 2 is dominated by species B. skagerrakensis and explains ∼11 % (EMB046/20-3GC) and ∼9 % (EMB046/10-4GC) of the variance. Finally, factor 3 explains ∼3 % (EMB046/20-3GC) and ∼1.8 % (EMB046/10-4GC) of the variance and includes N. iridea as a common species for both records, with P. osloensis as the second dominant species for the EMB046/20-3GC record and Stainforthia fusiformis for the EMB046/10-4GC record, consequently resulting in “N. iridea–P. osloensis” (EMB046/20-3GC) and “N. iridea–S. fusiformis” (EMB046/10-4GC) assemblages (Table 3). The individual factor weights (importance) for each counted sample are expressed by factor loading (Fig. 3). Factors with loadings above 0.5 are considered to be the most significant. The factor analysis shows that both records are defined by a clear dominance of the P. osloensis factor alternating with the B. skagerrakensis factor between CE 900 and 1700. The most pronounced changes in the foraminiferal assemblages occur between CE 1700 and the present day when the P. osloensis factor is to a large extent replaced by the B. skagerrakensis factor. Similar long-term variability is seen in the `palaeoproductivity fauna' group due to a strong dominance of B. skagerrakensis in this group (Fig. 3). In addition, palaeoproductivity fauna appears in higher abundance between  CE 900 and 1200. In the uppermost part ( CE 1950) of the EMB046/10-4GC record the N. iridea- S. fusiformis factor distinctly increased, while the N. iridea- P. osloensis factor of the EMB046/20-3GC record shows less variability (Fig. 3).

Figure 5Comparison of absolute abundances of planktonic foraminifera, total organic carbon (TOC), stable carbon isotope (δ13C), absolute abundance of palaeoproductivity fauna, the CABFAC results, stable oxygen isotope composition of seawater (δw), Mg∕Ca-derived bottom water temperature (Mg∕Ca-derived BWT), and stable oxygen isotopes (δ18O) of the two studied sediment cores EMB046/10-4GC (green curves) and EMB046/20-3GC (black curves). Reconstructions of the NAO index (yellow curve – Trouet et al., 2009, orange curve – Olsen et al. (2012), red curve – Faust et al., 2016) and winter precipitation (blue curve – Bakke et al. (2008)). The thicker curves correspond to a five-point running average. The errors bands represent uncertainties of the records. Dashed line divide record into two periods of the most pronounced palaeoproductivity changes, which are discussed in the text.

6 Discussion

## 6.1 Productivity changes in the last millennium

All dominant species in our benthic foraminiferal assemblages, grouped into factors, have documented association with quality (e.g. fresh or decaying) and availability of organic matter at the sea floor (e.g. Conradsen et al., 1994; Alve and Murray, 1995, 1997; Alve, 2003; Gustafsson and Nordberg, 2001; Duffield et al., 2015). The absolute abundance of planktonic foraminifera, stable carbon isotopes and total organic carbon also provide information on past variability of productivity. We combine these proxies to assess productivity changes in the Skagerrak. Two periods with different productivity in the Skagerrak region are identified: (i) moderate productivity between  CE 900 and 1700 and (ii) high productivity from  CE 1700 until the present (Figs. 3, 5). For each defined period, we discuss the level of and changes in primary productivity and potential causes behind this productivity variability. Throughout the discussion, we also refer to the smoothed data of stable isotopes and Mg∕Ca-derived BWT records as palaeothermometry proxies.

### 6.1.1 Moderate primary productivity (∼ CE 900–1700)

The highest absolute abundance of planktonic foraminifera and a clear dominance of the P. osloensis factor are recorded between  CE 900 and 1700, indicating a period of moderate primary productivity (Fig. 3). The TOC values and δ13C values do not show any major changes within this interval and until  CE 1500, respectively (Fig. 5). A high abundance of planktonic foraminifera is strongly correlated with nutrient-rich water, making them a good proxy for productivity changes (Boltovskoy and Correa, 2016). The P. osloensis factor includes the species P. osloensis and N. iridea, with an addition of C. laevigata in EMB046/10-4GC, in line with the previously identified C. laevigata–P. osloensis cluster (Erbs-Hansen et al., 2011) and N. iridea–C. laevigata category (Alve, 2010). These three species have an ecological preference for nutrient-rich environments, preferably with oxic bottom water conditions (Alve, 2003, 2010; Duffield et al., 2015). Hence, overall nutrient-rich conditions likely prevailed in the Skagerrak during this period.

Furthermore, the B. skagerrakensis factor and palaeoproductivity fauna also peak occasionally in the early and the late part of this interval ( CE 900–1200,  CE 1600–1700) (Figs. 3, 5). The B. skagerrakensis factor relates to ecological preferences of the epifaunal to shallow infaunal benthic species B. skagerrakensis, a species associated with high fresh phytodetritus fluxes to the sea floor accompanied by a continuously high oxygen content in the sediments (Duffield et al., 2015). The abundant occurrence of B. skagerrakensis in the Skagerrak and Oslofjord area is restricted to water masses with temperatures between 5 and 7 C and salinity around 35 PSU (Qvale and Nigam, 1985, and references therein; Alve and Murray, 1995, 1997; Duffield et al., 2015). In contrast to P. osloensis, N. iridea and C. laevigata, the taxon B. skagerrakensis does not feed on decaying organic matter but prefers freshly settled algal material (Duffield et al., 2015). Hence, the appearance of B. skagerrakensis between  CE 900 and 1200 and  CE 1600–1700 in both records likely indicates a period of well-oxygenated bottom water conditions with high fresh phytodetritus fluxes, while the period from  CE 1200 to 1600, characterised by dominance of the P. osloensis factor (Fig. 3), suggests that bottom water oxygen conditions in the Skagerrak were somewhat less favourable for B. skagerrakensis and/or a change to a more food-competitive environment where the herbivorous B. skagerrakensis was likely out-competed by the more omnivorous to detritivorous species C. laevigata, P. osloensis and N. iridea, which are all able to feed on both fresh and decaying organic matter (Alve, 2010; Duffield et al., 2015).

Because the stable carbon isotope composition recorded in calcareous benthic foraminiferal shells can be used to reconstruct past bottom water environments modified by fluxes of organic matter (Rohling and Cooke, 1999; Ravelo and Hillaire-Marcel, 2007), corresponding changes of the benthic δ13C similar to those in foraminiferal assemblages would be expected. Marine organisms preferentially take up more 12C than 13C in their biomass. When this organic matter disintegrates after it is deposited at the ocean floor, more 12C is released to the surrounding water. Hence, enhanced degradation at the bottom, e.g. that related to enhanced primary productivity in the surface waters, will increase the 12C concentrations in the bottom or pore water, in addition to increasing the nutrient content. Foraminifera that calcify in such a 12C-enriched water mass will record lower δ13C values than if less degradation of organic matter took place (Ravelo and Hillaire-Marcel, 2007; Filipsson and Nordberg, 2010). Hence, while the changes in the benthic and planktonic foraminiferal assemblages at  CE 900–1700 suggest overall constant nutrient-rich conditions characteristic of the Skagerrak region, the δ13C display light values and little variability until  CE 1500, supporting our interpretation of moderate primary production within the interval. This corresponds to studies by Hebbeln et al. (2006), who recorded increasing δ13C values after a period of high productivity in the southern Skagerrak and relatively little variability in δ13C in the northern Skagerrak between  CE 700 and 1500.

### 6.1.2 Causes of moderate primary productivity (∼ CE 900–1700)

Within the period of moderate primary production common for both sites, the δw, BWT and δ18O records reveal low correlation between both cores, pointing to different water conditions in the central and western Skagerrak (Fig. 5). At the time of higher BWT in the earliest ∼250 years of the EMB046/20-3GC record followed by a decreasing temperature trend until  CE 1450, the temperature in EMB046/10-4GC is lower and less variable between  CE 900 and 1350, after which it increases. The two events with the most contrasting temperatures between the cores were found at  CE 900–1100 and  CE 1350–1500. In the first instance, warming in the central Skagerrak seen from higher BWT is not indicated by the δ18O, which instead shows higher values and corresponds to decreases in δw. Since the δ18O can be induced by salinity and temperature, periods of good correlation between δ18O and δw but contradictory to the BWT pattern may suggest the influence of salinity on the temperature signal (Brückner and Mackensen, 2006). In contrast, good correlation between all three proxies is believed to give a fair estimation of temperature and salinity changes (Fig. 4). During the second instance ( CE 1350–1500), there is a general good correlation between the BWT changes and variability of δ18O values, showing similar patterns of higher BWT and δw, corresponding to a drop in δ18O values in EMB046/10-4CG records, and low BWT and δw when δ18O values increase in EMB046/20-3GC (Fig. 4). These changes a reflect colder bottom water temperature and lower salinity from  CE 1350 to 1500 in the central Skagerrak than in the western Skagerrak, which interestingly coincide with an event of minimum surface salinity in the north-eastern Skagerrak interpreted as enhanced outflow of low saline Baltic Sea water (Hebbeln et al., 2006). Since our cores were retrieved bellow 400 m within the deep waters of Norwegian Trench we would expect to obtain similar temperature and salinity signals for both sites. Instead, the cooling observed in the central Skagerrak most likely resulted from a renewal of the deep water by inflowing colder and denser North Sea waters, which apparently did not reach the shallower-located EMB046/10-4GC (Ljøen and Svansson, 1972). Thus, it is possible that the higher BWT of EMB046/10-4GC rather reflected the temperature of warm Atlantic water occupying the western Skagerrak basin (Brückner and Mackensen, 2006; Butruille et al., 2017).

The above-described changes in oxygen isotopes and BWT do not appear to correspond to variability seen in foraminiferal assemblages and the δ13C records. While it is difficult to find a good match between changes in palaeoproductivity and palaeotemperature proxies, the high absolute abundance of planktonic foraminifera recorded at both sites at  CE 900–1700 (Fig. 5) suggests that the primary productivity was driven by nutrient-rich Atlantic water and abundant phytodetritus fluxes rather than enhanced nutrients entering the area through the Baltic outflow feeding the NCC. This interpretation is supported by a study from the northern North Sea, where Klitgaard-Kristensen and Sejrup (1996) argued that the Atlantic water is favourable for planktonic foraminifera, while the low salinity of the NCC reduces their abundance. Our interpretation is further supported by previous studies based on foraminiferal (Erbs-Hansen et al., 2011) and diatom (Gil et al., 2006) assemblages, as well as a multi-proxy study by Hebbeln et al. (2006), which all report on an onset of enhanced Atlantic water advection to the Skagerrak at  CE 900. Gil et al. (2006) documented an increase in diatom species associated with high-salinity water in the Skagerrak and argued for enhanced inflow of nutrient-rich water via the NJC. Moreover, water with suspended sediments and low salinity are not favourable for planktonic foraminifera (Murray, 1976); therefore the higher absolute abundance of planktonic foraminifera in EMB046/10-4GC than in EMB046/20-3GC can be explained by an advantageous exposure to Atlantic water and smaller contribution of the low-salinity Baltic Sea water within the NCC in the western Skagerrak than in the central Skagerrak. It has to be noted, however, that planktonic foraminifera, due to their ability to inhabit the water column down to ∼200 m (for species found in our records), can reflect the character of both surface and upper intermediate water layers (Jonkers et al., 2010; Schiebel et al., 2017).

Differences in water conditions between both sites are further supported by the differences in foraminiferal assemblages at CE ∼900–1700 between the cores. It is likely that the higher abundance of C. laevigata in core EMB046/10-4GC than in the EMB046/20-3GC reflects a higher contribution of well-oxygenated Atlantic waters to the western Skagerrak compared to its central part (Fig. 3). C. laevigata and P. osloensis are mostly recorded in the Skagerrak and Norwegian Trench area and are associated with Atlantic water influence (Van Weering and Qvale, 1983; Conradsen et al., 1994; Alve and Murray, 1995; Klitgaard-Kristensen et al., 2002; Wollenburg et al., 2004). In EMB046/20-3GC, C. laevigata is largely replaced by N. iridea, which is commonly present in the Skagerrak and the Scandinavian fjord waters with a salinity > 35 PSU and a temperature range of 6–6.5 C (Polovodova Asteman et al., 2013, and references therein). In addition, N. iridea is capable of growth under hypoxic–suboxic conditions (Duffield et al., 2015). Hence, the dominance of N. iridea over C. laevigata in EMB046/20-3GC may be related to less favourable bottom water oxygen conditions in the central Skagerrak, which is more exposed to the brackish and nutrient-rich water of the BC as well as to enhanced river runoff.

Based on the combined high absolute abundance of planktonic foraminifera, the intermediate abundance of benthic palaeoproductivity species, occasional peaks in the B. skagerrakensis factor, no major changes in TOC and relatively little variability in δ13C, we conclude that the time interval CE 900–1700 was characterised by moderate palaeoproductivity in the Skagerrak. Furthermore, we argue that palaeoproductivity does not show coherence to changes in palaeothermometry proxies (Mg∕Ca, δ18O) and that the moderate productivity at that time was largely sustained by nutrient-rich Atlantic water bathing the sites, as deduced from the appearance of planktonic foraminifera.

From the discussion above it is seen that several processes took place between  CE 900 and 1700, including changes in bottom water circulation, oxygen and salinity fluctuations, or carbon fluxes. Each of these may in turn have been influenced by anthropogenic, climatic and/or oceanic factors (Brückner and Mackensen, 2008; Filipsson and Nordberg, 2010). Interestingly, the beginning of the moderate productivity period characterised by peaks in the B. skagerrakensis factor and palaeoproductivity species ( CE 900–1200) corresponds well to the overall stable and relatively warm temperatures at the Northern Hemisphere (CE 830–1100) (PAGES 2k Consortium, 2013) associated with the early stage of the Medieval Climate Anomaly (MCA) (e.g. Hass, 1996). Hass (1996) argued that the MCA lasted until CE 1300 in the Skagerrak area. Based on granulometric analyses he suggested that the MCA was associated with a decreased strength of south-westerly winds as a result of a more northerly located cyclonic track causing weaker bottom currents. In contrast, the North Atlantic Oscillation reconstructions by Trouet et al. (2009), Olsen et al. (2012), Faust et al. (2016), among others, all suggest a tendency for prevailing positive NAO conditions during the MCA and, hence, south-westerlies dominating the meteorological regimes during winter (Fig. 5). During a positive NAO phase, strong south-westerlies result in warm and wet winters over northern Europe (Hurrell, 1995; Hurrell et al., 2001; Trouet et al., 2009) and coincide with intensification of water mass exchange (inflows and outflows) in the Skagerrak (Winther and Johannessen, 2006). Predominant positive NAO conditions would, however, also intensify the river runoff and the outflow of brackish water from the Baltic Sea to the Skagerrak due to increased precipitation over the catchment area. Increased riverine input and Baltic Sea outflow would in turn enhance the nutrient supply to the surface waters of the Skagerrak and, hence, increase the primary production. However, our palaeoproductivity proxies do not provide information on increased productivity during  CE 900–1500. Neither does the high absolute abundance of planktonic foraminifera recorded at the same time in both Skagerrak sites (Figs. 3, 5) support enhanced riverine input and increased Baltic Sea outflow, since their abundance usually decreases in areas with increased brackish water conditions and decreased water transparency due to, for example, runoff (Murray, 1976; Klitgaard-Kristensen and Sejrup, 1996). Therefore, the higher winter precipitation during the MCA reported from south-western Norway was likely blocked by the mountain ranges in southern Norway resulting in less runoff reaching our study sites (Fig. 5, Bakke et al., 2008).

The long period of high planktonic foraminiferal abundance, depleted δ18O and overall warm BWT in the central Skagerrak corresponds well with a pronounced positive NAO phase reconstructed by Trouet et al. (2009) and Olsen et al. (2012) lasting until  CE 1450, which suggests a strong advection of warm Atlantic water (Fig. 5). Around that time, a period characterised by a deep water warming and weaker deep and cold North Sea water inflows (deep-water renewal) was suggested to take over in the Skagerrak (Butruille et al., 2017). The following drop in temperature and thus cooling in the central Skagerrak is consistent with a temperature decline in the North Atlantic around  CE 1400, marking the onset of the Little Ice Age (LIA) in northern Europe (Berstad et al., 2003; Brückner and Mackensen, 2006; Büntgen et al., 2011; Erbs-Hansen et al., 2011). Among other reconstructions, Trouet et al. (2009) proposed that during the LIA the NAO index became more negative, and that the associated weaker westerly airflow resulted in a reduction of the Atlantic water inflow. Reduced advection of warm Atlantic water coinciding with a deep-water renewal (discussed above) to the central Skagerrak may explain the colder bottom water conditions seen around  CE 1500 in the EMB046/20-3GC record. The changing climate conditions during the transition between MCA and the LIA and later during the LIA, likely accompanied by increased storminess (Gil et al., 2006), were also an important additional factor behind the contradicting bottom water conditions between the central and western Skagerrak.

To conclude, it is likely that the moderate productivity during  CE 900–1450 was primarily driven by the enhanced influence of nutrient-rich Atlantic water, likely related to the predominant positive NAO associated with the warm Medieval Climate Anomaly. The Atlantic inflow was stronger than the Baltic Sea outflow, creating favourable conditions for planktonic foraminifera in the region. The second part of the moderate productivity period ( CE 1450–1700) coincides with variable climate conditions, which are characteristic of the Little Ice Age, where a predominantly more negative NAO likely reduced the warm Atlantic water inflow and trigged deep-water renewal.

### 6.1.3 High primary productivity (∼ CE 1700–present)

A prominent change in benthic and planktonic foraminiferal assemblages, the δ13C and TOC records took place around CE 1700 in both cores, suggesting a shift in environmental conditions in the Skagerrak. The palaeoproductivity-related foraminiferal fauna increased and the values of the B. skagerrakensis factor largely replaced the P. osloensis factor (Figs. 3, 5). The planktonic foraminifera content decreases towards the top of the cores and almost disappear after  CE 1850 in both records. These changes in foraminiferal assemblages were accompanied by a gradual reduction of δ13C and a continuous increase in TOC in both cores (Fig. 5). All proxies suggest enhanced primary productivity from CE 1700 until the present day.

### 6.1.4 Causes of high primary productivity (∼ CE 1700–present)

Consistent with the previous period, comparison of the two cores shows that the periods of high primary productivity are common for both study sites while the BWT and oxygen isotopes generally provide negative correlation in temperature and salinity between the central and western Skagerrak. However, from  CE 1700 onward the δ18O and δw patterns show some similarities as comparably low values in both records and reflect increased values around  CE 1800, coinciding with generally low BWTs, which are again warmer at the western Skagerrak side (Figs. 4–5). These changes suggest a colder climate followed by warmer and more saline bottom water conditions. The cooling appears during a distinctly negative NAO period, showing a similar relation between changes in the water masses and the NAO as the one seen at  CE 1350–1500. However, we suggest that at this time deep-water renewal reached to the western Skagerrak, consequently lowering temperature and salinity at both sites.

The following warming of the bottom water after  CE 1800 reflects naturally induced environmental changes in the Skagerrak region accompanied by gradually increasing human activity. The warming seen in our records is consistent with intensified heat transport to the Northern Hemisphere (Brückner and Mackensen, 2006), an overall warm climate in Fennoscandia (Briffa et al., 1992), and warm spring conditions recorded between CE 1750 and 1920 off the Norwegian continental margin (Berstad et al., 2003). Furthermore, our results are supported by similar δ18O changes recorded in the same area and at the same time by Brückner and Mackensen (2006), and Hass (1996). The bottom water warming was associated with the termination of the coldest LIA phase in the Skagerrak region (Berstad et al., 2003; Brückner and Mackensen, 2006), a transition to a more positive NAO mode that would entail wetter and warmer winters over Scandinavia (Hurrell et al., 2001), or both. Either way, our data demonstrate a long-term intensification of nutrient supply, likely due to increased inflows of Atlantic water to the Skagerrak.

Associated with a generally more positive NAO phase, changes in the atmospheric and oceanic circulation systems are expected to result in enhanced surface water outflow from the Baltic Sea (e.g. Gustafsson and Stigebrandt, 1996; Zorita and Laine, 2000). As a consequence of increased precipitation and thus enhanced river runoff over the large Baltic Sea catchment area, the outflowing low saline Baltic Sea water would supply the Skagerrak surface water with nutrients (Svansson, 1975; Aure et al., 1998; Krossa et al., 2015). The decreased, almost disappearing, abundance of planktonic foraminifera in our records is in line with lower salinity and transparency in the upper water layers. The total freshwater riverine discharge from the Baltic Sea together with the contributions from the major Norwegian rivers to the Skagerrak contribute much fewer nutrients to the Skagerrak waters than nutrient transport via the inflows from the North Sea (Danielssen et al., 1997). However, Krossa et al. (2015) still showed a good correlation between the increased Baltic outflow and enhanced productivity in the Skagerrak on timescales longer than 1100 years, based on an increased alkenone C37:4 concentration, a proxy for the influence of brackish water. Furthermore, Polovodova Asteman et al. (2018) documented increased palaeoproductivity over the last 1700 years that corresponded in time with the increased alkenone C37:4 concentration (Krossa et al., 2015). Both studies argued that the nutrient-rich water causing enhanced productivity in the central Skagerrak was to a large extent of Baltic origin. Thus, the Baltic outflow can play an important role in the Skagerrak nutrient budget.

Zillén et al. (2008) showed that a widespread oxygen deficiency in the Baltic Sea was stimulated by an increased nutrient loading associated with a growing population and intensification of land use changes, which all began around CE 1600 and were followed by industrial development at around CE 1800. In more recent times (after  CE 1900), land use changes in Scandinavia caused an increased terrestrial runoff through either sparsely cultivated lands after massive deforestation or due to extensive farming in southern Sweden (Zillén et al., 2008, and references therein; Kaplan et al., 2009). Hence it is likely that, when superimposed on the natural variability in volume of the outflowing Baltic water, the concentration of nutrients in the outflowing water has changed after  CE 1700 until the present due to a gradually increased human impact.

At the same time as increased primary production caused eutrophication in the Baltic Sea, our data show a clear decrease in the δ13C values at both Skagerrak sites (Fig. 5). This distinct reduction in δ13C provides evidence for a change from the oceanic–atmospheric relationship established for the preceding periods towards an additional contribution of the lighter carbon isotope (12C) to the seawater from the atmosphere due to the increase in atmospheric CO2 concentration caused by anthropogenic emissions, known as the Suess effect (e.g. Cage and Austin, 2010; Filipsson and Nordberg, 2010; Eide et al., 2017). Both the increased primary productivity and the Suess effect may cause a reduction of benthic δ13C. However, the decrease of ca 0.9 ‰ seen in our δ13C records from  CE 1800 until the present is likely to a large extent explained by the Suess effect, which is in line with the ca 0.8 ‰ δ13C decrease between the pre-industrial and modern periods observed in the North Atlantic Ocean (Eide et al., 2017).

To summarise, from  CE 1700 until the present day, changes in the palaeoproductivity proxies indicate increased primary productivity likely caused by an enhanced nutrient content driven by a combination of the influence of nutrient-rich warm Atlantic water, enhanced Baltic outflow, intensified river runoff, and enhanced human impact through agriculture expansion and industrial development. This increase in primary production occurred during a high-variability period in temperature and NAO index.

### 6.1.5 Changes in the last ∼70 years

During the last 70 years B. skagerrakensis started to decrease in favour of higher S. fusiformis abundances (Fig. 3). S. fusiformis is considered as an indicator of trophic changes in Scandinavian waters due to its high tolerance of oxygen-depleted and organic-rich conditions (e.g. Alve, 2003, and references therein). This opportunistic species has the highest reproduction (up to 7 times per month) and growth rates after the phytoplankton blooms, followed by an enhanced food supply to the sea floor and decay of organic matter (Gustafsson and Nordberg, 2001). This explains the taxon's food preferences recognised as both fresh phytodetritus and microbes associated with the degradation of organic matter (Duffield et al., 2015). Therefore, the increased abundance of S. fusiformis during the last 70 years and simultaneous drop in B. skagerrakensis suggest changes in the quality of organic matter at the sea floor after high-productivity episodes, causing increased enrichment of organic matter and enhanced degradation in the sediments. At the same time, continuously high nutrients content coinciding with a gradual decline in oxygen concentration has been shown for the Skagerrak fjords with sluggish bottom water circulation (e.g. Rosenberg, 1990; Johannessen and Dahl, 1996; Alve, 2003; Filipsson and Nordberg, 2004). There are no hydrographic studies reporting on low-oxygen conditions in the deep Skagerrak basin during the last 70 years; thus, high appearance of S. fusiformis in our records will indicate a rather higher amount of degraded food than depleted oxygen.

7 Conclusions

This study provides evidence for changes in primary productivity in the Skagerrak during the last millennium. Our multi-proxy records show that the time interval  CE 900–1700 was characterised by moderate primary production with nutrients largely sustained by warm Atlantic water, as revealed by a high abundance of planktonic foraminifera. The first part of this interval ( CE 900–1450) was likely associated with the warm Medieval Climate Anomaly, during which a persistent positive NAO strengthened the westerlies, resulting in more frequent warm Atlantic water inflows. After  CE 1450, continuously moderate productivity at both sites are indicated by a dominance of the P. osloensis factor, high abundance of planktonic foraminifera, relatively stable TOC and overall stable δ13C. This continuously moderate productivity at both sites coincided with the variable climate conditions characteristic of the Little Ice Age. Episodes of negative NAO-triggered deep-water renewal resulted in colder bottom water temperature in the central Skagerrak, while the western Skagerrak seems to be more resistant to this cooling and instead reflects the temperature of warm Atlantic water that occupied this site.

Finally, the high primary productivity period between  CE 1700 and 2000 is documented by an increase in both TOC and the B. skagerrakensis factor, high absolute abundance of palaeoproductivity fauna and decreased δ13C values. Enhanced nutrient availability was likely caused by a stronger Baltic Sea outflow, increased river runoff, intensified inflows of the nutrient-rich Atlantic water, together with agricultural and industrial expansion. Simultaneously, an increase in human-induced CO2 emissions caused a great change in the oceanic carbon isotope budget, indicated by the Suess effect, shown in our records by strongly negative δ13C values since  CE 1800. The most pronounced increase in primary production at  CE 1800–2000 occurred during a warm period with a more positive NAO, wetter and warmer winters in Scandinavia as shown by an increase in BWT, and a decrease in δ18O in our records.

A comparison between the two records show an slight difference in species composition overall. This difference in species composition is likely due to a more favourable habitat in the western Skagerrak, with less exposure to low-salinity nutrient-rich NCC water, while the higher TOC in the central Skagerrak mostly results from an exposure to the Baltic outflow, nutrient-rich water reaching the site via NJC and terrestrial runoff due to a more inland location. The productivity and the fluxes of organic matter to the seafloor appear to not correspond to the temperature and salinity changes recorded in the benthic Melonis barleeanus shells.

Data availability
Data availability.

The presented data are available at https://www.pangaea.de/, last access: 2 October 2018, https://doi.org/10.1594/PANGAEA.894130 (Binczewska et al., 2018).

Supplement
Supplement.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

This research is a part of the ClimLink project, which was funded by Norway Grants: POL-NOR/199763/92/2014 in the Polish-Norwegian Programme operated by the National Centre of Research and Development of Poland. We thank the captains, chief scientists and crews of R/V Elisabeth Mann-Borgese for logistical and technical assistance. We also thank Małgorzata Bąk (University of Szczecin) for coordinating the project. Joanna Sławińska, Ryszard Borówka and staff of the Laboratory of the Department of Geology and Palaeogeography (University of Szczecin, Poland) performed geochemistry analyses, while Are Olsen (University of Bergen) and Jeroen Groeneveld (University of Bremen) contributed with valuable comments regarding the δ13C – the Suess effect relationship and the Mn∕Ca ratio, respectively. Rocio Castano Primo helped with archiving data at Pangaea. Finally, we thank Marit-Solveig Seidenkrantz and one anonymous reviewer for constructive feedback.

Edited by: Markus Kienast
Reviewed by: Marit-Solveig Seidenkrantz and one anonymous referee

References

Alve, E.: A common opportunistic foraminiferal species as an indicator of rapidly changing conditions in a range of environments, Estuar. Coast. Shelf S., 57, 501–514, 2003.

Alve, E.: Benthic foraminiferal responses to absence of fresh phytodetritus: a two-year Experiment, Mar. Micropaleontol., 76, 67–75, 2010.

Alve, E. and Murray, J. W.: Benthic foraminiferal distribution and abundance changes in Skagerrak surface sediments: 1937 (Höglund) and 1992/1993 data compared, Mar. Micropaleontol., 25, 269–288, 1995.

Alve, E. and Murray, J. W.: High benthic fertility and taphonomy of foraminifera: a case study of the Skagerrak, North Sea, Mar. Micropaleontol., 31, 157–175, 1997.

Andersson, L.: Trends in nutrient and oxygen concentrations in the Skagerrak-Kattegat, J. Sea Res., 35, 63–71, 1996.

Aure, J. and Dahl, E.: Oxygen, nutrients, carbon and water exchange in the Skagerrak Basin, Cont. Shelf Res., 14, 965–977, 1994.

Aure, J., Danielssen, D., and Sætre, R.: Assessment of eutrophication in Skagerrak coastal waters using oxygen consumption in fjordic basins, ICES J. Mar. Sci., 53, 589–595, 1996.

Aure, J., Danielssen, D., and Svendsen, E.: The origin of Skagerrak coastal water off Arendal in relation to variations in nutrient concentrations, ICES J. Mar. Sci., 55, 610–619, 1998.

Baden, S. P., Loo, L. O., Pihl, L., and Rosenberg, R.: Effects of eutrophication on benthic communities including fish: Swedish west coast, AMBIO, 19, 113–122, 1990.

Bakke, J., Lie, O., Dahl, S. O., Nesje, A., and Bjune, A. E.: Strength and spatial patterns of the Holocene wintertime Westerlies in the NE Atlantic region, Global Planet Change, 60, 28–41, 2008.

Barker, S., Greaves, M., and Elderfield, H.: A study of cleaning procedures used for foraminiferal Mg∕Ca paleothermometry, Geochem. Geophy. Geosy., 4, 8407, https://doi.org/10.1029/2003GC000559, 2003.

Berstad, I. M., Sejrup, H. P., Klitgaard-Kristensen, D., and Haflidason, H.: Variability in temperature and geometry of the Norwegian Current over the past 600 yr; stable isotope and grain size evidence from the Norwegian margin, J. Quaternary Sci., 18, 591–602, 2003.

Binczewska, A., Risebrobakken, B., Polovodova Asteman, I., Moros, M., Tisserand, A., Jansen, E., and Witkowski, A.: Multi-proxy study on two sediment cores in western and central Skagerrak, https://doi.org/10.1594/PANGAEA.894130, 2018.

Boltovskoy, D. and Correa, N.: Foraminifers (Planktonic), in: Encyclopedia of Marine Geosciences, edited by: Harff, J., Meschede, M., Petersen, S., and Thiede, J., Springer, the Netherlands, 1, 255–261, 2016.

Boon, A. R., Duineveld, G. C. A., and Kok, A.: Benthic Organic Matter Supply and Metabolism at Depositional and Non-depositional Areas in the North Sea, Estuar. Coast. Shelf S., 49, 747–761, 1999.

Boyle, E. A.: Manganese carbonate overgrowths on foraminifera, Geochim. Cosmochim. Ac., 47, 1815–1819, 1983.

Breitburg, D. L., Hondorp, D. W., Davias, L. A., and Diaz R. J.: Hypoxia, Nitrogen, and Fisheries: Integrating Effects Across Local and Global Landscapes, Annu. Rev. Mar. Sci., 1, 329–349, 2009.

Briffa, K. R., Jones, T. S., Bartholin, D., Eckstein, F. H., Schweingruber, W. K., Zetterber, M., Karlén, W., Zetterberg, P., and Eronen, M.: Fennoscandian summers from ad 500: temperature changes on short and long timescales*, Clim. Dynam., 7, 111–119, 1992.

Brückner, S. and Mackensen, A.: Deep-water renewal in the Skagerrak during the last 1200 years triggered by the North Atlantic Oscillation: evidence from benthic foraminiferal δ18O, The Holocene, 16, 331–340, 2006.

Brückner, S. and Mackensen, A.: Organic matter rain rates, oxygen availability, and vital effects from benthic foraminiferal δ13C in the historic Skagerrak, North Sea, Mar. Micropaleontol., 66, 192–207, 2008.

Büntgen, U., Tegel, W., Nicolussi, K., McCormick, M., Frank, D., Trouet, V., Kaplan, J. O., Herzig, F., Heussner, K. U., Wanner, H., Luterbacher, J., and Esper, J.: 2500 Years of European Climate Variability and Human Susceptibility, Science, 3, 578–582, https://doi.org/10.1126/science.1197175, 2011.

Butruille, C., Krossa, V. R., Schwab, C., and Weinelt, M.: Reconstruction of mid- to late-Holocene winter temperatures in the Skagerrak region using benthic foraminiferal Mg∕Ca and δ18O, The Holocene, 27, 63–72, 2017.

Cage, A. G. and Austin, W. E. N.: Marine climate variability during the last millennium: The Loch Sunart record, Scotland, UK, Quaternary Sci. Rev., 29, 1633–1647, 2010.

Cardinale, M. and Svedäng, H.: Modelling recruitment and abundance of Atlantic cod, Gadus morhua, in the eastern Skagerrak–Kattegat (North Sea): evidence of severe depletion due to a prolonged period of high fishing pressure, Fish. Res., 69, 263–282, 2004.

Carstensen, J., Conley, D. J., Andersen, J. H., and Ærtebjerg, G.: Coastal eutrophication and trend reversal: A Danish case study, Limnol. Oceanogr., 51, 398–408, 2006.

Conley, D. J., Carstensen, J., Vaquer-Sunyer, R., and Duarte, C. M.: Ecosystem thresholds with hypoxia, Hydrobiologia, 629, 21–29, doi10.1007/s10750-009-9764-2, 2009.

Conradsen, K. H., Bergsten, K. L., Knudsen, K. N., Nordberg, K., and Seidenkrantz, M. S.: Recent benthic foraminiferal distribution in the Kattegat and the Skagerrak, Scandinavia, Cushman Foundation Special Publication, 32, 53–68, 1994.

Danielssen, D. S., Edler, L., Fonsehus, S., Hernroth, L., Ostrowski, M., Svendsen, E., and Talpsepp, L.: Oceanographic variability in the Skagerrak and northern Kattegat, May–June, 1990, ICES J. Mar. Sci., 54, 753–773, 1997.

Doney, S. C., Fabry, V. F., Feely, R. A., and Kleypas, J. A.: Ocean Acidification The Other CO2 Problem, Annu. Rev. Mar. Sci., 1, 169–192, 2009.

Duffield, C. J., Hess, S., Norling, K., and Alve, E.: The response of Nonionella iridea and other benthic foraminifera to “fresh” organic matter enrichment and physical disturbance, Mar. Miropaleontol., 120, 20–30, 2015.

Eide, M., Olsen, A., Ninnemann, U. S., and Eldevik, T.: A global estimate of the full oceanic 13C Suess effect since the preindustrial, Global Biogeochem. Cy., 31, 492–514, https://doi.org/10.1002/2016GB005472, 2017.

Erbs-Hansen, D. R., Knudsen, K. L., Gary, A. C., Gyllencreutz, R., and Jansen, E.: Holocene climatic development in the Skagerrak, eastern North Atlantic: foraminiferal and stable isotope evidence, The Holocene, 22, 301–312, 2011.

FAO (Food and Agriculture Organization of the United Nations): The State of World Fisheries and Aquaculture 2016, Contributing to food security and nutrition for all, Rome, Italy, 107–137, 2016.

Faust, J. C., Fabian, K., Milzer, G., Giraudeau, J., and Knies, J.: Norwegian fjord sediments reveal NAO related winter temperature and precipitation changes of the past 2800 years, Earth Planet. Sc. Lett., 435, 84–93, 2016.

Filipsson, H. L. and Nordberg, K.: Climate Variations, an Overlooked Factor Influencing the Recent Marine Environment. An Example from Gullmar Fjord, Sweden, Illustrated by Benthic Foraminifera and Hydrographic Data, Estuaries, 27, 867–881, 2004.

Filipsson, H. L. and Nordberg, K.: Variations in organic carbon flux and stagnation periods during the last 2400 years in a Skagerrak fjord basin, inferred from benthic foraminiferal δ13C, in: Fjord Systems and Archives, edited by: Howe, J. A., Austin, W. E. N., Forwick, M., and Paetzel, M., Geol. Soc. Sp., London, 344, 261–270, https://doi.org/10.1144/SP344.18, 2010.

Filipsson, H. L., McCorkle, D. C., Mackensen, A., Berhard, J. M., Andersson, L. S., Naustvoll, L. J., Cabellero-Alfonso, A. M., Nordberg, K., and Danielssen, D. S.: Seasonal variability of stable carbon isotopes (δ13CDIC) in the Skagerrak and the Baltic Sea: Distinguishing between mixing and biological productivity, Palaeogeogr. Palaeocl., 483, 15–30, https://doi.org/10.1016/j.palaeo.2016.11.031, 2017.

Fonselius, S.: The upwelling of nutrients in the central Skagerrak, Deep-Sea Res. II, 43, 57–71, 1996.

Gil, I. M., Abrantes, F., and Hebbeln, D.: The North Atlantic Oscillation forcing through the last 2000 years: Spatial variability as revealed by high-resolution marine diatom records from N and SW Europe, Mar. Micropaleontol., 60, 113–129, 2006.

Greaves, M., Caillon, N., Rebaubier, H., Bartoli, G., Boharty, S., Cacho, I., Clarke, L., Cooper, M., Daunt, C., Delaney, M., DeMenocal, P., Dutton, A., Eggins, S., Elderfield, H., Garbe-Schoenberg, D., Goddard, E., Green, D., Groeneveld, J., Hastings, D., Hathorne, E., Kimoto, K., Klinkhammer, G. P., Labeyrie, L., Lea, D. W., Marchitto, T., Martinez-Boti, M. A., Mortyn, P. G., Ni, Y., Nürnberg, D., Paradis, G., Pena, L., Quinn, T., Rosenthal, Y., Russell, A., Sagawa, T., Sosdian, S., 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, Q08010, https://doi.org/10.1029/2008GC001974, 2008.

Gustafsson, B. and Stigebrandt, A.: Dynamics of the freshwater-influenced surface layers in the Skagerrak, J. Sea Res., 35, 39–53, 1996.

Gustafsson, M. and Nordberg, K.: Living (stained) benthic foraminiferal response to primary production and hydrography in the deepest part of the Gullmar Fjord, Swedish west coast, with comparisons to Höglund's 1927 material, J. Foramin. Res., 31, 2–11, 2001.

Hammer, Ø., Harper, D. A. T., and Ryan, P. D.: Past: Paleontological Statistics Software Package for Education and Data Analysis, Palaeontol. Electron., 4, 1–9, 2001.

Hasenfratz, A. P., Schiebel, R., Thornalley, D. J. R., Schönfeld, J., Jaccard, S. L., Martínez-García, A., Holbourn, A., Jennings, A. E., Kuhnt, W., Lear, C. H., Marchitto, T. M., Quillmann, U., Rosenthal, Y., Yu, J., and Haug, G. H.: Mg∕Ca-temperature calibration for the benthic foraminifera Melonis barleeanum and Melonis pompilioides, Geochim. Cosmochim. Ac., 217, 365–383, https://doi.org/10.1016/j.gca.2017.08.038, 2017.

Hass, H. C.: Northern Europe climate variations during late Holocene: Evidence from marine Skagerrak, Palaeogeogr. Palaeocl., 123, 121–145, 1996.

Haynert, K., Schönfeld, J., Polovodova-Asteman, I., and Thomsen, J.: The benthic foraminiferal community in a naturally CO2-rich coastal habitat of the southwestern Baltic Sea, Biogeosciences, 9, 4421–4440, https://doi.org/10.5194/bg-9-4421-2012, 2012.

Hebbeln, D., Knudsen, K. L., Gyllencreutz, R., Kristensen, P., Klitgaard-Kristensen, D., Backman, J., Scheurle, C., Jiang, H., Gil, I., Smelror, M., Jones, P. D., and Sejrup, H. P.: Late Holocene coastal hydrographic and climate changes in the eastern North Sea, The Holocene, 16, 987–1001, 2006.

Hjalmarsson, S., Anderson, L. G., and She, J.: The exchange of dissolved inorganic carbon between the Baltic Sea and the North Sea in 2006 based on measured data and water transport estimates from a 3D model, Mar. Chem., 121, 200–205, 2010.

Hop, H., Gjøsæter, J., and Danielssen, D. S.: Seasonal feeding ecology of cod (Gadus morhua L.) on the Norwegian Skagerrak coast, ICES J. Mar. Sci., 49, 453–461, https://doi.org/10.1093/icesjms/49.4.453, 1992.

Hurrell, J. W.: Decadal trends in the North Atlantic Oscillation: regional temperatures and precipitation, Science, 269, 676–679, 1995.

Hurrell, J. W., Kushnir, Y., and Visbeck, M.: The North Atlantic Oscillation, Science 291, 603–605, https://doi.org/10.1126/science.1058761, 2001.

ICES: Oceanographic DataSet – temperature (57.17–58 N, 8–9.79 E), Copenhagen, Denmark, available at: http://www.ices.dk (9 September 2011), 2010.

Iversen, S. A., Skogen, M., and Svendsed, E.: Availability of horse mackerel (Trachurus trachurus) in the north-eastern North Sea, predicted by the transport of Atlantic water, Fish. Oceanogr., 11, 245–250, 2002.

Johannessen, T. and Dahl, E.: Declines in oxygen concentrations along the Norwegian Skagerrak coast, 1927-1993: A signal of ecosystem changes due to eutrophication?, Limnol. Oceanogr., 41, 766–778, 1996.

Jonkers, L., Brummer, G.-J. A., Peeters, F. J. C., van Aken, H. M., and Femke De Jong, M.: Seasonal stratification, shell flux, and oxygen isotope dynamics of left-coiling N. pachyderma and T. quinqueloba in the western subpolar North Atlantic, Paleoceanography, 25, PA2204, https://doi.org/10.1002/palo.20018, 2010.

Kaplan, J. O., Krumhardt, K. M., and Zimmermann, N.: The prehistoric and preindustrial deforestation of Europe, Quaternary Sci. Rev., 28, 3016–3034, 2009.

Klitgaard-Kristensen, D. and Sejrup, H. P.: Modern benthic foraminiferal biofacies across the Northern North Sea, Sarsia, 81, 97–106, 1996.

Klitgaard-Kristensen, D., Sejrup, H. P., and Haflidason H.: Distribution of recent calcareous benthic foraminifera in the northern North Sea and relation to the environment, Polar. Res., 21, 275–282, 2002.

Kristiansen, T. and Aas, E.: Water type quantification in the Skagerrak, the Kattegat and off the Jutland west coast, Oceanologia, 57, 177–195, 2015.

Kristjánsdóttir, G. B., Lea, D. W., Jennings, A. E., Pak, D. K., and Belanger, C.: New spatial Mg∕Ca-temperature calibrations for three Arctic, benthic foraminifera and reconstruction of north Iceland shelf temperature for the past 4000 years, Geochem. Geophy. Geosy., 8, Q03P21, https://doi.org/10.1029/2006GC001425, 2007.

Krossa, V. R., Moros, M., Blanz, T., Jansen, E., and Schneider, R.: Late Holocene Baltic Sea outflow changes reconstructed using C37:4 content from marine cores, Boreas, 44, 81–93, https://doi.org/10.1111/bor.12093, 2015.

Ljøen, R. and Svansson, A.: Long-term variations of subsurface temperatures in the Skagerrak, Deep-Sea Res., 19, 277–288, 1972.

Mackensen, A., Schumacher, S., Radke, J., and Schmidt, D. N.: Microhabitat preferences and stable carbon isotopes of endobenthic foraminifera: Clue to quantitative reconstruction of oceanic new production?, Mar. Micropaleontol., 40, 233–258, 2000.

Micheli, F.: Eutrophication, Fisheries, and Consumer-Resource Dynamics in Marine Pelagic Ecosystems, Science, 285, 1396–1398, 1999.

Moros, M., Andersen, T. J., Schulz-Bull, D., Häusler, K., Bunke, D., Snowball, I., Kotilainen, A., Zillén, L., Jensen, J. B., Kabel, K., Hand, I., Leipe, T., Lougheed, B. C., Wagner, B., and Arz, H. W.: Towards an event stratigraphy for Baltic Sea sediments deposited since AD 1900: approaches and challenges, Boreas, 46, 129–142, 2017.

Murray, J. W.: A method of determining proximity of marginal seas to an ocean, Mar. Geol., 22, 103–119, 1976.

Olsen, E. M., Knutsen, H., Gjøsæter, J., Jorde, P. E., Knutsen, J. A., and Stenseth, N. C.: Life-history variation among local populations of Atlantic cod from the Norwegian Skagerrak coast, J. Fish Biol., 64, 1725–1730, 2004.

Olsen, E. M., Lilly, G. R., Heino, M., Morgan, M. J., Brattey, J., and Dieckmann, U.: Assessing changes in age and size at maturation in collapsing populations of Atlantic cod (Gadus morhua), Can. J. Fish Aquat. Sci., 62, 811–823, 2005.

Olsen, J., Anderson, N. J,, and Knudsen, M. F.: Variability of the North Atlantic Oscillation over the past 5200 years. Nat. Geosci., 5, 808–812, https://doi.org/10.1038/ngeo1589, 2012.

PAGES 2k Consortium: Continental-scale temperature variability during the past two millennia, Nat. Geosci., 6, 339–346, https://doi.org/10.1038/ngeo1797, 2013.

Pingree, R. D., Holligan, P. M., Mardell, G. T., and Harris, R. P.: Vertical distribution of plankton in the Skagerrak in relation to doming of the seasonal thermocline, Cont. Shelf Res., 1, 209–219, 1982.

Polovodova Asteman, I., Nordberg, K., and Filipsson, H. L.: The Little Ice Age: evidence from a sediment record in Gullmar Fjord, Swedish west coast, Biogeosciences, 10, 1275–1290, https://doi.org/10.5194/bg-10-1275-2013, 2013.

Polovodova Asteman, I., Risebrobakken, B., Moros, M., Binczewska, A., Dobosz, S., Jansen, E., Sławinska, J., and Bąk, M.: Late Holocene palaeoproductivity changes: a multi-proxy study in the Norwegian Trench and the Skagerrak, North Sea, Boreas, 47, 238–255, https://doi.org/10.1111/bor.12264, 2018.

Qvale, G. and Nigam, R.: Bolivina skagerrakensis, a new name for Bolivina cf. robusta, with notes on ecology and distribution, J. Foramin. Res., 15, 6–12, 1985.

Ravelo, A. C. and Hillaire-Marcel, C.: The Use of Oxygen and Carbon Isotopes of Foraminifera in Paleoceanography, in: Developments in Marine Geology, edited by: Hillaire-Marcel, C. and de Vernal, A., Elsevier, Amsterdam, 1, 735–764, 2007.

Reimer, P. J., Bard, E., Bayliss, A., Beck, J. W., Blackwell, P. G., Bronk Ramsey, C., 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., Hoffman, 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.

Rohde, J.: On the dynamics of the large-scale circulation of the Skagerrak, J. Sea Res., 35, 9–21, 1996.

Rohling, E. J. and Cooke, S.: Stable oxygen and carbon isotopes in foraminiferal carbonate shells, in: Modern Foraminifera, edited by: Sen Gupta, B. K., Kluwer Academic Publishers, Dordrecht/Boston/London, 297, 239–258, 1999.

Rosenberg, R., Gray, J. S., Josefson, A. B., and Pearson, J. H.: Petersen's benthic stations revisited. II. Is the Oslofjord and eastern Skagerrak enriched?, J. Exp. Mar. Biol. Ecol., 105, 219–251, 1987.

Rosenberg, R., Elmgre, E., Fleischer, S., Jonsson, P., Persson, G., and Dahlin, H.: Marine eutrophic case studies in Sweden, AMBIO, 19, 102–108, 1990.

Rydberg, L., Haamer, J., and Liungman, O.: Fluxes of water and nutrients within and into the Skagerrak, J. Sea Res., 35, 23–38, 1996.

Rydberg, L., Ærtebjerg, G., and Edler, L.: Fifty years of primary production measurements in the Baltic entrance region, trends and variability in relation to land-based input of nutrients, J. Sea Res., 56, 1–16, 2006.

Schiebel, R., Spielhagen, R. F., Garnier, J., Hagemann, J., Howa, H., Jentzen, A., Martínez-Garcia, A., Meilland, J., Michel, e., Repschläger, J., Salter, I., Yamasaki, M., and Haug, G.: Modern planktic foraminifers in the high-latitude ocean, Mar. Micropaleontol., 136, 1–13, 2017.

Shackleton, N. J.: Attainment of isotopic equilibrium between ocean water and the benthonic foraminifera genus Uvigerina: isotopic changes in the ocean during the last glacial, Colloques Internationaux du C.N.R.S., 203–209, 1974.

Sigman, D. M. and Hain, M. P.: The Biological Productivity of the Ocean, Nature Education, 3, 1–16, 2012.

Skogen, M. D., Budgell, W. P., and Rey, F.: Interannual variability in Nordic seas primary production, ICES J. Mar. Sci., 64, 889–898, doi.org/10.1093/icesjms/fsm063, 2007.

Stuiver, M., Reimer, P. J., and Reimer, R. W.: CALIB 7.1 (wwwprogram), available at: http://calib.org, last access: 16 February 2017.

Svansson, A.: Physical and chemical oceanography of the Skagerrak and the Kattegat, Fishery Board of Sweden, Institute of Marine Research, Report, 1, 1–88, 1975.

Thomas, H., Bozec, Y., de Baar, H. J. W., Elkalay, K., Frankignoulle, M., Schiettecatte, L.-S., Kattner, G., and Borges, A. V.: The carbon budget of the North Sea, Biogeosciences, 2, 87–96, https://doi.org/10.5194/bg-2-87-2005, 2005.

Thomsen, J., Stapp, L. S., Haynert, K., Schade, H., Danelli, M., Lannig, G. K., Wegner, M., and Melzner, F.: Naturally acidified habitat selects for ocean acidification–tolerant mussels, Sci. Adv., 3, e1602411, https://doi.org/10.1126/sciadv.1602411, 2017.

Trouet, V., Esper, J., Graham, N. E., Baker, A., Scourse, J. D., and Frank, D. C.: Persistent positive North Atlantic oscillation mode dominated the Medieval Climate Anomaly, Science, 324, 78–80, 2009.

Van Weering, T. C. E. and Qvale, G.: Recent sediments and foraminiferal distribution in the Skagerrak, northeastern North Sea, Mar. Geol., 52, 75–99, 1983.

Winther, N. G. and Johannessen, J. A.: North Sea circulation: Atlantic inflow and its destination, J. Geophys. Res.-Oceans, 111, C12018, https://doi.org/10.1029/2005JC003310, 2006.

Wollenburg, J. E., Knies, J., and Mackensen, A.: High-resolution paleoproductivity fluctuations during the past 24 kyr as indicated by benthic foraminifera in the marginal Arctic Ocean, Palaeogeogr. Palaeocl., 204, 209–238, 2004.

Zillén, L., Conley, D. J., Andrén, T., Andrén, E., and Björck, S.: Past occurrences of hypoxia in the Baltic Sea and the role of climate variability, environmental change and human impact, Earth-Sci. Rev., 91, 77–92, 2008.

Zorita, E. and Laine, A.: Dependence of salinity and oxygen concentrations in the Baltic Sea on large-scale atmospheric circulation, Clim. Res., 14, 25–41, 2000.