# Cold-Water Coral Mound Archive Provides Unique Insights Into Intermediate Water Mass Dynamics in the Alboran Sea During the Last Deglaciation

Robin Fentimen, Eline Feenstra, Andres Rüggeberg, Torsten Vennemann, Irka Hajdas, Thierry Adatte, David van Rooij, Anneleen Foubert

<figure><img src="image_1.png" /></figure>

## To cite this version:

Robin Fentimen, Eline Feenstra, Andres Rüggeberg, Torsten Vennemann, Irka Hajdas, et al. Cold-Water Coral Mound Archive Provides Unique Insights Into Intermediate Water Mass Dynamics in the Alboran Sea During the Last Deglaciation. Frontiers in Marine Science, 2020, 7, 10.3389/fmars.2020.00354. hal-04278641

HAL Id: hal-04278641

https://univ-angers.hal.science/hal-04278641v1

Submitted on 10 Nov 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Cold-Water Coral Mound Archive Provides Unique Insights Into Intermediate Water Mass Dynamics in the Alboran Sea During the Last Deglaciation

Robin Fentimen¹*, Eline Feenstra¹, Andres Rüggeberg¹, Torsten Vennemann², Irka Hajdas³, Thierry Adatte², David Van Rooij⁴ and Anneleen Foubert¹

¹ Department of Geosciences, University of Fribourg, Fribourg, Switzerland, ² Faculty of Geosciences and Environment, University of Lausanne, Lausanne, Switzerland, ³ Laboratory of Ion Beam Physics, ETH Zürich, Zurich, Switzerland, ⁴ Department of Geology, Ghent University, Ghent, Belgium

## OPEN ACCESS

**Edited by:**
Ana Colaço,
Marine Research Institute (IMAR),
Portugal

**Reviewed by:**
Furu Mienis,
Royal Netherlands Institute for Sea
Research (NIOZ), Netherlands
Lorenzo Angeletti,
Institute of Marine Science (CNR), Italy

*Correspondence:
Robin Fentimen
robin.fentimen@unifr.ch;
robin.fentimen@gmail.com

**Specialty section:**
This article was submitted to Deep-Sea Environments and Ecology, a section of the journal Frontiers in Marine Science

Received: 09 January 2020

Accepted: 27 April 2020

Published: 16 June 2020

**Citation:**
Fentimen R, Feenstra E, Rüggeberg A, Vennemann T, Hajdas I, Adatte T, Van Rooij D and Foubert A (2020) Cold-Water Coral Mound Archive Provides Unique Insights Into Intermediate Water Mass Dynamics in the Alboran Sea During the Last Deglaciation. Front. Mar. Sci. 7:354. doi: 10.3389/fmars.2020.00354

The Alboran Sea is widely recognized to host numerous cold-water coral ecosystems, including the East Melilla Coral Province. Yet, their development through time and response to climatic variability has still to be fully understood. Based on a combined investigation of benthic foraminiferal assemblages, foraminiferal stable isotope compositions, grain size analysis, sediment geochemistry, and macrofaunal quantification, this study identifies key events and processes having governed cold-water coral development at the East Melilla Coral Province between Greenland Stadial 2.1 and the Early Holocene. The transition from Greenland Stadial 2.1 to Greenland Interstadial 1 is associated to a decline of bryozoan communities and their replacement by cold-water corals, together with changes in benthic foraminiferal assemblages and a decrease in the sediment mean grain size. These results suggest that a rapid decrease in bottom currents and the establishment of dysoxic and mesotrophic conditions at the seafloor, possibly associated to enhanced fluvial input, resulted in the decline of bryozoans as the dominant suspension feeding organisms and their replacement by a thriving cold-water coral community. This transition from a bryozoan to a coral dominated environment is concomitant with the beginning of the African Humid Period, confirming that increasing fluvial input could have been a main factor triggering the establishment of cold-water corals in the East Melilla Coral Province during Greenland Interstadial 1. A change in benthic foraminiferal communities and an increase in the sediment mean grain size mark the passage from the Early to Late Greenland Interstadial 1. The current velocity of intermediate water masses is suggested to have increased during the Early to Late Greenland Interstadial 1, whilst simultaneously fluvial input would have reduced. Such changes suggest that the climate became more arid during the second phase of Greenland Interstadial 1.

**Keywords:** benthic foraminifera, stable isotopes, grain size, Mediterranean, paleoenvironment, Bølling-Allerød, Greenland Interstadial 1
# INTRODUCTION

Scleractinian cold-water corals (CWCs) shape worldwide, unique and highly diverse environments (Freiwald et al., 2004; Roberts et al., 2009). They settle in areas favorable for their suspension feeding life style, such as steep cliffs, rock overhangs or pre-existing seamounts, forming with time large CWC reefs (Freiwald, 2002; Roberts et al., 2009). The most prevalent scleractinian CWC species in the North Atlantic Ocean and Mediterranean Sea are *Desmophyllum pertusum* (formerly known as *Lophelia pertusa* (Linnaeus, 1758), see Addamo et al., 2016) and *Madrepora oculata* (Freiwald et al., 2004; Wheeler et al., 2007). Their development is promoted by high food availability coupled to enhanced hydrodynamic regimes, whilst the physico-chemical properties of the ambient seawater (e.g., density, salinity, temperature) also play a role (White et al., 2005; Mienis et al., 2007; Dullo et al., 2008; Davies et al., 2009; Davies and Guinotte, 2011; Hanz et al., 2019). A key characteristic of CWCs are their capacity to baffle suspended sediments entering their framework (Wilson, 1979; Mienis et al., 2009). The interplay between this sediment baffling and coral growth leads, over time, to the formation of CWC mounds (Roberts et al., 2009). These mound structures are essentially confined to continental margins between 200 and 1000 m depth, where they can form build-ups of over 300 m in height (Mienis et al., 2007; Roberts et al., 2009; Wienberg and Titschack, 2015; Hebbeln et al., 2016, 2019). The development of coral mounds depends primarily on the presence of dynamic bottom current regimes. Indeed, vigorous bottom currents increase lateral food supply, enabling corals to thrive, and support mound aggradation through sediment supply (White et al., 2005; Mienis et al., 2007; Hebbeln et al., 2016). The equilibrium between sediment and nutrient supply is fundamental for a coral mound to maintain sustained development through time. Although coral mounds demonstrate a discontinuous record, often punctuated by a number of hiatuses, they present much higher sedimentation rates than surrounding off-mound settings where sedimentation is limited (Rüggeberg et al., 2007; Eisele et al., 2008; Margreth et al., 2011; Raddatz et al., 2011; Stalder et al., 2018; Wienberg et al., 2018). As such, they possess a strong potential to recognize short term paleoenvironmental changes and may serve as unique high-resolution sedimentary records, preserving sediments that are otherwise eroded or non-deposited in neighboring open environments (Thierens et al., 2013).

Living scleractinian CWCs, mostly represented by *M. oculata*, are presently found in the Mediterranean Sea. They have, for example, been observed in the Northern Ionian Sea (Savini et al., 2014), in the Sicily channel (Freiwald et al., 2009), in the Tyrrhenian Sea (Taviani et al., 2017, 2019), in the Adriatic and Aegean Seas (Vafidis et al., 1997; Freiwald et al., 2009), in the Gulf of Lions (Zibrowius, 2003; Orejas et al., 2007; Fabri et al., 2014), and at the summit of ridges in the Alboran Sea (Hebbeln et al., 2009; Fink et al., 2013; Lo Iacono et al., 2014; Stalder et al., 2015, 2018; Wang et al., 2019). The oldest occurrences of fossil CWC mounds in the Mediterranean Sea date back to the Late Pliocene/Early Pleistocene and are related to the

recolonization of the Mediterranean Sea by CWCs after the Messinian Salinity crisis (Taviani et al., 2005; Vertino et al., 2014). Mediterranean coral mounds are concentrated in the Alboran Sea, along the Moroccan margin (Margreth et al., 2011; van Rooij et al., 2013; Lo Iacono et al., 2014). The largest CWC mound field, the Melilla mound field, covers an area in excess of 500 km² parallel to the continental margin (van Rooij et al., 2013; **Figure 1**). It can be subdivided into two coral provinces, the West and East Melilla Coral Provinces, respectively 35 km northwest and 15 km northeast of the Cape Tres Forcas (Lo Iacono et al., 2014; **Figure 1**). Coral mounds situated in the West Melilla Coral Province dwell at depths ranging from 300 to 430 m and reach heights between 8 and 21 m above the seafloor (Lo Iacono et al., 2014). Coral mounds that are part of the East Melilla Coral Province demonstrate higher morphological variability. The mounds that are located to the North of the East Melilla Coral Province possess the particularity of being situated along three steep volcanic ridges (Brittlestar Ridges I, II and III). The summits of these Brittlestar Ridges are at ca. 250 m depth, whilst their base is found at ca. 450 m depth (Hebbeln et al., 2009; Lo Iacono et al., 2014; Wang et al., 2019; **Figure 1**). They are between 100 and 250 m wide and extend up to 3 km in length (Hebbeln et al., 2009; Fink et al., 2013).

The recent (last 20 ka) history of coral proliferation of the East Melilla Coral Province has been documented by a number of studies in the past years (Fink et al., 2013, 2015; Stalder et al., 2015, 2018; Wang et al., 2019; Feenstra, 2020). The conclusions drawn from these different studies indicate that corals experienced their most intense recent proliferation phase during Greenland Interstadial 1 (GI-1) and during the Early Holocene, whilst coral growth was reduced or halted during the Mid to Late Holocene, Younger Dryas and Last Glacial Maximum. Fink et al. (2013, 2015) and Stalder et al. (2015, 2018) suggested that increased productivity was the main factor driving coral growth during GI-1 and the Early Holocene. Wang et al. (2019), based on high benthic foraminiferal accumulation rates, further suggested that high productivity is a decisive parameter sustaining mound development during GI-1 and the Early Holocene. Moreover, these authors suggest that a rearrangement of surface and intermediate water masses would also have triggered sustained coral growth at these times, whilst Feenstra (2020) further stress the importance of fluvial input and changes in water mass configuration for the initiation of coral growth during GI-1 and the Early Holocene. Based on these combined observations, it appears clear that GI-1 and its upper and lower boundaries are key periods for recent coral growth history in the Southeastern Alboran Sea. However, previous studies lack the sufficient resolution to enable a detailed investigation of the interaction between environmental changes and coral growth during GI-1. The aims of this study are thus: (1) to provide a high-resolution paleoenvironmental reconstruction of mound formation during the last deglaciation at the East Melilla Coral Province, with a particular focus on GI-1; (2) to assess in detail the variability of intermediate water masses during GI-1, and at the lower and upper GI-1 boundaries, in the Southeastern Alboran Sea.
<figure><img src="image_4.png" /></figure>

**FIGURE 1** | Location and oceanographic circulation of the study area. (A) General map of the Mediterranean Sea and location of the western Mediterranean (red rectangle). (B) Bathymetric map of the western Mediterranean. The red box east of the Cape Tres Forcas indicates the location of the East Melilla Coral Province (EMCP). The bold brown arrows illustrate the general circulation of the inflowing surface Modified Atlantic Water (MAW), which describes two anticyclonic gyres: the Western Alboran Gyre (WAG) and the Eastern Alboran Gyre (EAG). The Atlantic Anticyclonic Gyre (AAG) is illustrated by the thin brown arrows. The surface water circulation is based on a compilation of results from Lanoix (1974), Tintore et al. (1988), Viúdez and Tintoré (1995) and L'Helguen et al. (2002), whilst the bathymetry uses the GEBCO_2019 gridded bathymetric data. (C) Precise location of the Banc des Provençaux and Brittlestar Ridge I and II (BR I and BR II). The location of the studied core MD13-3455G is illustrated by the red dot, whereas the two yellow dots indicate the location of the neighboring cores MD13-3459G and MD13-3462G sampled during the same cruise. This map uses MultiBeam Echosounder data retrieved during the "GATEWAY" cruise (N°194) onboard the *Marion Dufresne II* Research Vessel (van Rooij et al., 2013).
# Regional Setting

## Geographical and Geological Setting

The Alboran Sea is the westernmost basin of the Mediterranean Sea and is approximately 400 km long and 200 km wide, with an average depth of 1300 m and a maximum depth of ca. 1800 m (Olivet et al., 1973; Comas et al., 1999). It is limited to the North by the Iberian Peninsula and to the South by North Africa. The Alboran Sea demonstrates an intricate seafloor morphology, with the presence of ridges, seamounts, plateaus, troughs and a number of sub-basins (Figure 1; Comas et al., 1999). This complex morphology is a result of the extension and subsidence that occurred in the region during the Early to Middle Miocene (Comas et al., 1999; Faccenna et al., 2004; Do Couto et al., 2016). The extension processes and associated crustal thinning led to the intrusion of the metamorphic basement by volcanic material (López Casado et al., 2001; Duggen et al., 2008). The shallow Banc des Provençaux, situated at ca. 200 m depth, is one of these volcanic edifices (Figure 1). Its limits extend out to form the three Brittlestar ridges (I, II and II) that are partly built up by CWCs and belong to the East Melilla Coral Province (Figure 1). Brittlestar Ridge I is situated at ca. 50 km from the mouth of the Moulouya river, which is one of North Africa's largest rivers (Emelyanov and Shimkus, 2012). The bases of the Brittlestar Ridges are surrounded by erosional moats, testimony to the strong currents affecting the area (Comas et al., 2009). Video surveys along the ridges revealed a diverse fauna consisting of numerous organisms, such as ophiuroids, soft corals and sponges. However, only scarce unregularly distributed living scleractinian coral colonies (D. pertusum and M. oculata) were observed, concentrated on the crest of the ridges (Hebbeln et al., 2009).

## Oceanographic Circulation

The circulation patterns in the semi-enclosed Mediterranean Sea are governed by the exchange between inflowing low salinity Atlantic Water and outflowing highly saline Mediterranean Outflow Water at the Strait of Gibraltar. The Alboran Sea is dominated by three water masses. The surface water consists of low salinity (< 36.5 psu) Modified Atlantic Water (MAW) that extends down to approximately 200 m depth and becomes increasingly saltier as it flows eastwards (La Violette, 1983; Millot, 2009). High saline (ca. 38.5 psu) and warm (ca. 13.3°C) Levantine Intermediate Water flows beneath the MAW between 200 and 600 m depth. It enters the Western Mediterranean through the Straits of Sicily to then exit through the Strait of Gibraltar, contributing to approximately 70% of the total Mediterranean Outflow Water (Millot, 2013). The Levantine Intermediate Water circulates along the Spanish margin, whilst the Southern Alboran Sea would be rather dominated by Shelf Water, i.e., a mixture of MAW and Western Mediterranean Deep Water (WMDW; Ercilla et al., 2016). Western Mediterranean Deep Water, colder (ca. 12.8°C) and slightly fresher (ca. 38.4 psu) than the above lying Levantine Intermediate Water and Shelf Water, further supplies the Mediterranean Outflow Water and makes up the deepest water mass in the region (Millot et al., 2006; Millot, 2013). The WMDW circulates essentially along the Moroccan margin in the Alboran Sea (Hernández-Molina et al., 2002; Ercilla et al., 2016). The surface circulation in the Alboran Sea varies

seasonally. The MAW enters the Northeastern Alboran Sea as a jet (1.6 Sv; 1 Sv = $10^6 \text{ m}^3 \cdot \text{s}^{-1}$) (Katz, 1972; Lanoix, 1974; La Violette, 1983). As a result of the entry of this jet, the surface circulation in the Alboran Sea describes two anticyclonic gyres, namely the Western Alboran Gyre and the Eastern Alboran Gyre (Heburn and La Violette, 1990; Lafuente et al., 1998; **Figure 1**). The strength of these two Alboran gyres fluctuates through time as a function of atmospheric and oceanographic parameters, such as the strength of westerly winds and the exchange rates between Modified Atlantic Water inflow and Mediterranean Outflow Water outflow (Heburn and La Violette, 1990). During winter, the Eastern Alboran Gyre tends to weaken and is replaced by a jet circulating eastward along the African coast (Tintore et al., 1988; Heburn and La Violette, 1990). A number of eddies develop along the gyres, such as the Almeria Cyclonic Gyre or the Atlantic Anticyclonic Gyre (L'Helguen et al., 2002; **Figure 1**). The East Melilla Coral Province is directly influenced by the Eastern Alboran Gyre and the Atlantic Anticyclonic Gyre (Lanoix, 1974; Viúdez and Tintoré, 1995; L'Helguen et al., 2002; **Figure 1**). Both gyres contribute to mixing between surface and intermediate water masses down to ca. 300 m depth (Heburn and La Violette, 1990; Lafuente et al., 1998).

# MATERIALS AND METHODS

## Sample Collection

This study is based on the detailed investigation of the 491 cm long MD13-3455G gravity core (35°26.337′N, 2°30.926′W, 319 m depth) recovered from Brittlestar Ridge I (East Melilla Coral Province) during the EUROFLEETS cruise MD194 Gateway "The Mediterranean-Atlantic Gateway Code: The Late Pleistocene Carbonate Mound Record" onboard the R/V *Marion-Dufresne II* (van Rooij et al., 2013; **Figure 1**). A 1 cm slice of each half core section was taken every 5 cm for the investigation of benthic foraminiferal assemblages. The same 5 cm sampling interval was chosen to pick planktonic and benthic foraminifera for the analysis of their stable oxygen and carbon isotopic composition and for the particle-size analysis of the sediment matrix. The core was sampled every 10 cm to measure the Total Organic Carbon of the sediment matrix. All results were plotted with the ggplot2 package for R (Wickham, 2016; R Core Team, 2018).

## Benthic Foraminiferal Assemblages

After drying and weighing, samples were wet sieved through a series of 200, 125, and 63 µm mesh sieves and dried at room temperature. Samples did not possess the same volume. The resulting fractions were then weighed separately. A minimum target number of 300 individuals was picked and identified per sample on the fraction > 125 µm. The sample was split with a microsplitter when necessary. Benthic foraminiferal densities were calculated by dividing the total number of foraminifera (number of individuals per aliquot multiplied by the split) by the dry fraction's weight. The benthic foraminiferal relative abundances were calculated for each sample and converted to percentages. The Shannon diversity index was calculated for
each sample using the PRIMER6 software (Clarke and Gorley, 2006). Multivariate statistical analysis was performed using the software PRIMER6 (Clarke and Gorley, 2006). Foraminiferal composition data was root transformed in order to down-weigh the contribution of the most abundant species, thus increasing the contribution of rarer species (Field et al., 1982). Bray-Curtis similarity coefficients were then calculated (Bray and Curtis, 1957) and the resulting similarity matrix was used to perform a hierarchical cluster analysis. The hierarchical cluster analysis groups samples that share high similarity levels and is used to define benthic foraminiferal assemblages that are specific to a given depth range/epoch. A Similarity Percentage (SIMPER) analysis was computed at two different similarity levels in order to precisely identify the species responsible for the similarity among and the dissimilarity between clusters identified in the hierarchical cluster analysis (Clarke and Gorley, 2006).

## Oxygen and Carbon Stable Isotope Analysis

A portion of the material prepared for the assessment of benthic foraminiferal assemblages was dry sieved through 250 and 212 µm mesh sieves. Between 5 and 12 specimens of the planktonic foraminifera *Globigerina bulloides* and the benthic foraminifera *Cibicides lobatulus* were picked from the 212–250 µm size fraction in order to measure their stable oxygen and carbon isotope compositions. This fraction was chosen to avoid any ontogenic effect on the measurements (Schiebel and Hemleben, 2017). Specimens were first cleaned three times with distilled water in an ultrasonic bath for 2 s to remove any attached particles or infilling. The measurements were then made using a Thermo Fisher Scientific GasBench II connected to a Thermo Finnigan Delta Plus XL isotope ratio mass spectrometer at the Stable Isotope Laboratory of the University of Lausanne according to a method adapted from Spötl and Vennemann (2003). Results are reported in the conventional δ-values in permil (‰) relative to the Vienna Pee Dee Belemnite (VPDB) standard. Analytical standard deviations (1σ) average 0.04‰ for δ¹³C and 0.06‰ for δ¹⁸O values based on 8 replicate analyses of standards in each sequence of 40 samples.

## Grain Size Analyses

Grain size was determined on the bulk material and on the siliciclastic fraction using the Malvern Mastersizer 3000 at the department of Geology of the University of Ghent. Large clasts (e.g., coral and bryozoan fragments) were first hand picked out of each sample. Samples were then placed in 35% H₂O₂ and boiled until the reaction stopped to remove organic matter. One half of each sample was boiled for 2 min in 10% HCl. This allowed dissolution of the carbonate material, hence only the siliciclastic fraction remained. Before measurement, all samples were placed in 2% sodium polymetaphosphate and boiled to assure total particle separation. A 2 mm mesh sieve placed above the receiver of the Malvern Mastersizer 3000 retained any remaining large particles. Each sample was measured three times and an average value was taken as the final result. The mean particle-size of the bulk and siliciclastic fractions, together with the mean sorted silt

size (i.e., the mean of the 10–63 µm siliciclastic fraction; McCave et al., 1995), were calculated using the Rysgran Package for R (Gilbert et al., 2015; R Core Team, 2018). The mean sorted silt size served as a proxy for bottom water current velocity (McCave and Hall, 2006), in combination with changes in benthic foraminiferal assemblages. Indeed, using solely mean sorted silt size as a bottom current velocity indicator is possibly biased in CWC environments, since coral framework can locally reduce bottom current speed and lead to the deposition of fine material (Huvenne et al., 2009; Eisele et al., 2011; Fentimen et al., 2020), thus leading to an underestimation of bottom current velocities.

## Organic Geochemistry

Total Organic Carbon (TOC) and Mineral Carbon (MinC) content were determined at the Laboratory of Sediment Geochemistry at the University of Lausanne using the Rock-Eval6 pyrolysis technique. The RockEval6 pyrolysis technique determines a Hydrogen Index (HI) and an Oxygen Index (OI), respectively corresponding to the quantity of pyrolyzable organic compounds relative to TOC and the quantity of CO₂ relative to TOC (Baudin et al., 2015). The graphical representation of HI vs. OI (i.e., the modified Van Krevelen diagram) has regularly been used to determine the type of kerogen and the origin of the organic matter present in samples, although the exactness of this method has been subject to discussion (Katz, 1983; Delvaux et al., 1989; van Krevelen, 1993; Tyson, 1995). The measured OI values in the samples presented in this study are extremely high and are probably the result of cracking at low temperature of poorly crystallized CaCO₃ (Baudin et al., 2015). Thus, the OI was not used to determine the type of kerogen. Instead, HI vs. Tmax was plotted in order to determine the kerogen type and the origin of the organic matter (Delvaux et al., 1989). Tmax being the temperature at which the maximum hydrocarbon yield resulting from kerogen cracking occurs (Baudin et al., 2015). Moreover, the high OI implies that TOC values are possibly overestimated (Behar et al., 2001).

## Radiometric Dating

Uranium-series isotope measurements (²³⁰Th/U) were carried out on 6 cold-water coral fragments (*D. pertusum* and *M. oculata*) using a multicollector inductively coupled plasma source mass spectrometer MC-ICPMS (*Thermo Fisher Scientific Neptuneplus*) coupled with a dissolver (*Aridus I*) at the Institute for Environmental Physics at Heidelberg University (Table 1). Prior to measurements, pristine-looking coral fragments were selected, physically cleaned by sand blasting and further chemically cleaned using a weak acid leaching. More details about the protocol can be found in Wefing et al. (2019) and in Frank et al. (2004).

Radiocarbon dating was performed at the ETH in Zürich on 5 benthic foraminiferal samples (Table 2). The epibenthic foraminiferal species *Cibicides lobatulus*, *Cibicides refulgens* and *Discanomalina coronata* were picked in order to obtain between 4 and 10 mg of pure carbonate. The samples were first dissolved in phosphoric acid. The extracted CO₂ was then converted to graphite and measured by Accelerator Mass Spectrometry (AMS) technique using the *MICADAS* dedicated instrument
<table><caption>**TABLE 1** | Uranium-series isotope (U/Th) ages of cold-water coral fragments.</caption><thead><tr><th>LAB ID</th><th>Depth (cm)</th><th>*S</th><th>**Age (years BP)</th><th>±</th><th><sup>238</sup>U (µ g/g)</th><th>±</th><th><sup>232</sup>Th (ng/g)</th><th>±</th><th>δ <sup>234</sup>U (‰)</th><th>±</th><th>δ <sup>234</sup>Ui (‰)</th><th>±</th></tr></thead><tbody><tr><td>IUP-8487</td><td>1</td><td><i>M</i></td><td>2328</td><td>0.013</td><td>3.655</td><td>0.00009</td><td>0.215</td><td>0.00054</td><td>147.75</td><td>0.43</td><td>148.75</td><td>0.43</td></tr><tr><td>IUP-8488</td><td>12</td><td><i>D</i></td><td>6759</td><td>0.024</td><td>3.210</td><td>0.00011</td><td>0.349</td><td>0.00071</td><td>146.88</td><td>0.58</td><td>149.74</td><td>0.59</td></tr><tr><td>IUP-8489</td><td>35</td><td><i>D</i></td><td>13725</td><td>0.049</td><td>3.404</td><td>0.00010</td><td>0.245</td><td>0.00073</td><td>144.02</td><td>0.56</td><td>149.75</td><td>0.59</td></tr><tr><td>IUP-8490</td><td>280</td><td><i>D</i></td><td>14501</td><td>0.054</td><td>3.214</td><td>0.00010</td><td>0.188</td><td>0.00050</td><td>143.12</td><td>0.54</td><td>149.13</td><td>0.56</td></tr><tr><td>IUP-8491</td><td>328</td><td><i>D</i></td><td>14519</td><td>0.052</td><td>3.294</td><td>0.00009</td><td>0.230</td><td>0.00072</td><td>144.44</td><td>0.58</td><td>150.52</td><td>0.60</td></tr><tr><td>IUP-8492</td><td>445</td><td><i>D</i></td><td>14747</td><td>0.052</td><td>3.278</td><td>0.00013</td><td>0.205</td><td>0.00061</td><td>145.34</td><td>0.75</td><td>151.55</td><td>0.78</td></tr></tbody></table>

<table><caption>**TABLE 2** | Radiocarbon ages of epibenthic foraminifera (species selected: *Cibicides lobatulus*, *Cibicides refulgens* and *Discanomalina coronata*).</caption><thead><tr><th>LAB ID</th><th>Depth (cm)</th><th><sup>14</sup>C age (years BP)</th><th>± 1σ</th><th>2σ lower (years BP)</th><th>2σ upper (years BP)</th><th>2σ median (years BP)</th></tr></thead><tbody><tr><td>ETH-87756</td><td>37</td><td>12273</td><td>30</td><td>13170</td><td>13510</td><td>13350</td></tr><tr><td>ETH-87755</td><td>157</td><td>13033</td><td>35</td><td>13900</td><td>14650</td><td>14170</td></tr><tr><td>ETH-87753</td><td>382</td><td>13538</td><td>38</td><td>14800</td><td>15510</td><td>15180</td></tr><tr><td>ETH-87752</td><td>447</td><td>13758</td><td>32</td><td>15220</td><td>15790</td><td>15500</td></tr><tr><td>ETH-87751</td><td>472</td><td>14566</td><td>38</td><td>16330</td><td>16970</td><td>16650</td></tr></tbody></table>

(Synal et al., 2007). Results were corrected for <sup>13</sup>C and calibrated using the Marine13 calibration curve (Reimer et al., 2013) and the software OxCal v4.2.4 (Ramsey, 2017). A reservoir age of 390 ± 80 years was applied to all ages (Siani et al., 2000).

# RESULTS

## Chronology

The chronostratigraphic framework of core MD13-3455G was constructed based on absolute radiocarbon and Uranium-series isotope (U/Th) ages (Tables 1, 2) combined with the planktonic (G. bulloides) and benthic (C. lobatulus) oxygen isotope records. High planktonic and benthic δ<sup>18</sup>O values, respectively between ca. 2 to 3‰ and 3 to 4‰, characterize the base of the core from 490 to 440 cm. Considering the U/Th and radiocarbon ages, this interval is interpreted as Greenland Stadial 2.1 as defined by Rasmussen et al., 2014 (GS-2.1; Figure 2). Greenland Interstadial 1 is defined as the period with lower planktonic and benthic δ<sup>18</sup>O values between 440 and 15 cm (Rasmussen et al., 2014). The following decrease in planktonic and benthic δ<sup>18</sup>O values (respectively from about 2 to 0.5 and from 3 to 1.5‰) from 15 cm to the top of the core represents the Holocene (Figure 2). The Younger Dryas would have possibly been eroded or not deposited, although a short excursion toward higher benthic δ<sup>18</sup>O values, respectively from about 2 to 3‰, is situated between 25 and 15 cm. Low planktonic δ<sup>18</sup>O values (about 2‰) would however suggest that this interval does not correspond to the Younger Dryas and that these deposits are indeed absent in this core. Such discontinuity and the presence of hiatuses is a regularly observed feature in cores recovered from CWC mounds (Dorschel et al., 2005; Rüggeberg et al., 2007; Frank et al., 2009; Raddatz et al., 2011; Stalder et al., 2018). The resulting chronology hence indicates that a great part of the core, from ca. 440 to

15 cm, corresponds to GI-1 (Figure 2). Greenland Stadial 2.1 and Greenland Interstadial 1 can be related respectively to the time periods formerly known as the Oldest Dryas and the Bolling-Allerød. The latest chronology defined by Rasmussen et al. (2014) will be used throughout this study.

## Core Description and Macrofauna

Scleractinian CWCs are absent during GS-2.1 (Figure 2). This interval is characterized by a dense accumulation of the bryozoan *Buskea dichotoma* associated to the brachiopod species *Gryphus vitreus* (Figure 2). The upper part of the core, corresponding to GI-1 and the Holocene, consists of framework building scleractinian CWCs (*D. pertusum* and *M. oculata*) embedded in a matrix of mixed, clay to silt sized carbonate/siliciclastic material (Figure 2; Feenstra, 2020). Sedimentary and macrofaunal variations observed in core MD13-3455G are reported in more detail by Feenstra (2020).

## Micropaleontology

### Relative Abundances and Diversity

A total number of 167 benthic foraminifera species were recognized (Supplementary Material). Benthic foraminiferal densities vary between 653 ind.g<sup>-1</sup> during the Late GI-1 (32 cm) and 10560 ind.g<sup>-1</sup> during the Early GI-1 (302 cm; Supplementary Material). The Shannon diversity index (Figure 3A) reaches minimum values during the GS-2.1 (1.831) and maximum values during the Holocene (3.589). Higher benthic foraminiferal diversity also characterizes the Early GI-1, whilst the Late GI-1 (from 230 to 25 cm) is marked by lower values (Figures 3A,B). The distribution of the main benthic foraminifera species is illustrated in Figures 3A,B. As all Bolivinid species (*Bolivina alata*, *B. difformis*, *B. dilatata*, *B. pseudoplicata*, *B. pseudopunctata*, *B. spathulata*, *B. striatula* and *B. variabilis*) follow the same distribution pattern throughout

<footnote>*S, coral species (D, Desmophyllum pertusum; M, Madrepora oculata). **Ages are corrected for the influence of seawater and/or detrital derived <sup>230</sup>Th using a <sup>230</sup>Th/<sup>232</sup>Th ratio of 8 ± 4. Errors are 2σ.</footnote>

<footnote>Ages are corrected for a reservoir age of 390 ± 80 years (Siani et al., 2000).</footnote>
<figure><img src="image_5.png" /></figure>

**FIGURE 2** | Stratigraphy, core description and macrofaunal content of core MD13-3455G. The core descriptions and macrofaunal content were previously reported by Feenstra (2020). The stratigraphy is constructed by combining the *G. bulloides* planktonic (light blue curve) and *C. lobatulus* benthic (dark blue curve) δ¹⁸O records (‰ VPDB), U/Th coral and calibrated median ¹⁴C foraminifera ages (Tables 1, 2). The resulting stratigraphic column is given to the far right (according to Rasmussen et al., 2014). The benthic foraminiferal assemblages (BFA) identified at 47% similarity and at 62% similarity thanks to the hierarchical cluster analysis (Figure 4) are indicated.

the core, they are grouped together as *Bolivina* spp. For the same reasons, the same grouping is applied to Buliminid species (*Bulimina aculeata*, *B gibba*, *B. marginata* and *B. striata*) and to Miliolids (*Figures 3A,B*). *Miliolinella subrotunda* is the most abundant species of the Miliolid group, making up on average ca. 64% of the total abundance. The GS-2.1 period is marked by high abundances of *C. lobatulus* (ca. 55%), *Astrononion stelligerum* (ca. 18%), *Trifarina angulosa* (ca. 12%), *D. coronata* (ca. 10%) and *Hyalinea balthica* (ca. 9%). Greenland Stadial 2.1 and GI-1 are set apart by an important change in the benthic foraminiferal fauna (**Figure 3**). Species typical of GS-2.1 disappear or strongly retreat and are replaced by newcomers. Infaunal species such as *Bulimina* spp. (ca. 50%), *Bolivina* spp. (ca. 25%), *Cassidulina laevigata* (ca. 24%), *Chilostomella oolina* (ca. 20%), *Spiroplectammina wrightii* (ca. 10%) and *Cassidulinoides bradyi* (ca. 7%), reach their highest abundance during GI-1 (**Figure 3**). Other abundant species are the epibenthic *C. lobatulus* (ca. 35%), *Rosalina bradyi* (ca. 13%) and

*Gavelinopsis praegeri* (ca. 10%). The Late GI-1, in comparison to the Early GI-1, is marked by higher *Bolivina* spp., *C. bradyi*, *C. laevigata* and *C. lobatulus* abundances, whilst *G. praegeri*, *R. bradyi* and *S. wrightii* demonstrate lower abundances. The Holocene is characterized by an increase in the abundance of the infaunal *Globocassidulina subglobosa* (from ca. 4 to 13%) and *Uvigerina mediterranea* (from ca. 1 to 10%), whilst the infaunal *C. laevigata* (from ca. 15 to 7%) and *Bolivina* spp. (from ca. 10 to 4%) abundances decrease (**Figure 3**). However, the Holocene only represents a thin horizon in core MD13-3455G, thus not allowing a precise appreciation of the benthic foraminiferal community present at the time. Miliolids are present throughout the core, reaching minimum and maximum abundances during GI-1 (ca. 3 and ca. 17%; **Figure 3B**).

## Benthic Foraminiferal Assemblages

The hierarchical cluster analysis allowed to distinguish 2 clusters at 47% similarity and 5 clusters and an outlier sample at 62%
similarity (Figure 4). Each cluster relates to a specific benthic foraminiferal assemblage (BFA). At 47% similarity, clusters separate between a GS-2.1 assemblage (BFA.A) and a combined GI-1 and Holocene assemblage (BFA.B). The most contributing species to BFA.A are *C. lobatulus* (ca. 30%), *A. stelligerum* (ca. 12%) and *T. angulosa* (ca. 9%), while *C. lobatulus* (ca. 14%), *C. laevigata* (ca. 10%) and *M. subrotunda* (ca. 6%) contribute the most to BFA.B. More detailed results of the SIMPER analysis are given in Tables 3, 4. At 62% similarity, clusters separate between a GS-2.1 assemblage (BFA.b), a transitional GS-2.1/GI-1 assemblage from 437 to 447 cm (BFA.d), an Early GI-1 assemblage from 432 to 232 cm (BFA.f), a Late GI-1 assemblage from 232 to 17 cm (BFA.e) and a Holocene assemblage from 12 cm to the core top (BFA.c) (Figure 4). The 5 species contributing the most to each cluster are reported in Table 5, whereas the species contributing the most to dissimilarity between clusters are given in Table 6. Overall, the clustering confirms and provides more details about benthic foraminiferal distribution patterns described in section “Relative Abundances and Diversity.”

## Stable Carbon Isotopic Records

The δ¹³C values of the benthic *C. lobatulus* have a range between 0.7‰ during the Late GI-1 (157 cm) and 1.4‰ during the Late GI-1 (337 cm), while the δ¹³C values of the planktonic *G. bulloides* range between -2.1‰ during the Late GI-1 (367 cm) and -0.8‰ during the Holocene (7 cm; **Figure 5**). The planktonic δ¹³C values have a higher variability compared to the benthic δ¹³C values (**Figure 5**). Average planktonic and benthic δ¹³C values are comparatively high during GS-2.1 (ca. -1.2 and 1.3‰ respectively) (**Figure 5**). The planktonic δ¹³C values have more variability during GI-1. Greenland Interstadial 1 can be divided into two phases; a first phase between ca. 440 and 210 cm and a second phase between 210 and 25 cm, respectively corresponding to average lower (ca. -1.6‰) and higher (-1.2‰) δ¹³C values. The benthic δ¹³C record can also be divided into the same two phases during GI-1, although in contrast to the planktonic δ¹³C values, the early phase corresponds to average higher benthic δ¹³C values and the later phase to average lower values. These two phases correspond to two distinct benthic foraminiferal assemblages:

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /></figure>

**FIGURE 3** | (A) Benthic foraminiferal Shannon diversity and distribution of the main benthic foraminifera species (expressed as percentages). The benthic foraminiferal assemblages (BFA) identified at 47% similarity and at 62% similarity thanks to the hierarchical cluster analysis (Figure 4) are illustrated. The stratigraphic column defined in Figure 2, the *G. bulloides* planktonic and *C. lobatulus* benthic $\delta^{18}\text{O}$ records and the U/Th coral and calibrated median $^{14}\text{C}$ benthic foraminifera ages are further indicated. (B) Distribution of the main benthic foraminifera species (expressed as percentages). The benthic foraminiferal assemblages (BFA) identified at 47% similarity and at 62% similarity thanks to the hierarchical cluster analysis (Figure 4) are illustrated. The stratigraphic column defined in Figure 2, the *G. bulloides* planktonic and *C. lobatulus* benthic $\delta^{18}\text{O}$ records and the U/Th coral and calibrated median $^{14}\text{C}$ benthic foraminifera ages are further indicated.

f and e (Figure 5). Planktonic and benthic $\delta^{13}\text{C}$ values show an overall decoupled trend, noticeably during the Late GI-1 (between ca. 160 and 180 cm; Figure 5). An increase in planktonic $\delta^{13}\text{C}$ values (from ca. -1.6 to -0.8‰) marks the transition toward the Holocene, also corresponding to a change in benthic foraminiferal assemblage (BFA.e to BFA.c). Benthic $\delta^{13}\text{C}$ values remain relatively high and invariable during the Holocene (ca. 1.2‰; Figure 5).

## Sediment Characterization

### Grain Size

The mean grain size of the bulk material varies between ca. 5 and 60 µm (Figure 5). The coarsest intervals are limited to GS-2.1 (ca. 60 µm), Late GI-1 (ca. 36 µm) and Holocene (ca. 40 µm). A number of coarser intervals are observed during GI-1 (Figure 5). The mean grain size of the siliciclastic fraction ranges

from ca. 5.3 to 9.4 µm, i.e., in the fine silt size range (Figure 5). The transition from GS-2.1 to GI-1 is characterized by a decrease in the mean grain size of the siliciclastic fraction from ca. 7.0 to 5.3 µm. Greenland Interstadial 1 is divided into different phases based on the mean grain size of the siliciclastic fraction. Between 440 and 280 cm, mean grain size is low (ca. 6 µm). This first phase is followed by a slight increase from ca. 6 to 7 µm between 280 and 180 cm (Figure 5). An important coarsening takes place at 180 cm (from ca. 7 up to 9.4 µm). This coarser interval (average mean grain size of 9 µm) spans from 180 to 120 cm (Figure 5). At 120 cm, the mean grain size decreases to ca. 8 µm and remains high until the end of GI-1, although the latest part of GI-1 (from 50 to 15 cm) shows a slight increase in mean grain size (ca. 8.5 µm). The passage from GI-1 to the Holocene is marked by a decrease in mean grain size of the siliciclastic fraction, from ca. 8.6 to 6.2 µm (Figure 5).
<figure><img src="image_8.png" /><figcaption>**FIGURE 4** | Hierarchical cluster analysis obtained from the benthic foraminifera Bray-Curtis similarity matrix. Benthic foraminifera abundances were root transformed. Two clusters are distinguished at 47% similarity (A,B). At a higher similarity level (62%), the cluster analysis allows to isolate 5 clusters and one outlier sample. Each cluster corresponds to a specific benthic foraminiferal assemblage. The species contributing the most to the similarity among and the dissimilarity between clusters are given in **Tables 3–6**.</figcaption></figure>

<table><caption>**TABLE 3** | List of the 5 most contributing species to the similarity within clusters A and B, which separate at 47% similarity on the hierarchical cluster analysis (Figure 5).</caption><thead><tr><th colspan="5">Average similarity: 63.08</th></tr><tr><th>Species</th><th>Av.Abund</th><th>Av.Sim</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>29.62</td><td>18.12</td><td>28.72</td><td>28.72</td></tr><tr><td>Astrononion stelligerum</td><td>11.57</td><td>8.34</td><td>13.23</td><td>41.95</td></tr><tr><td>Trifarina angulosa</td><td>8.66</td><td>7.14</td><td>11.32</td><td>53.27</td></tr><tr><td>Cassidulina laevigata</td><td>6.56</td><td>5.2</td><td>8.24</td><td>61.51</td></tr><tr><td>Discanomalina coronata</td><td>6.18</td><td>4.18</td><td>6.63</td><td>68.14</td></tr></tbody><thead><tr><th colspan="5">Average similarity: 58.22</th></tr><tr><th>Species</th><th>Av.Abund</th><th>Av.Sim</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>14.09</td><td>10.51</td><td>18.04</td><td>18.04</td></tr><tr><td>Cassidulina laevigata</td><td>10.05</td><td>7.28</td><td>12.51</td><td>30.55</td></tr><tr><td>Miliolinella subrotunda</td><td>6.21</td><td>4.59</td><td>7.89</td><td>38.45</td></tr><tr><td>Rosalina brady</td><td>5.82</td><td>3.92</td><td>6.73</td><td>45.17</td></tr><tr><td>Chilostomella oolina</td><td>5.63</td><td>3.47</td><td>5.97</td><td>51.14</td></tr></tbody></table>

The average abundance, similarity, contribution (%) and cumulative contribution (%) are indicated for each species. The average similarity of the cluster is also given.

The sortable silt mean grain size ranges from ca. 16.9 to 21.6 µm (Figure 5). The trend in mean@silt silt matches the one observed for the siliciclastic fraction (Figure 5). Highest values are recorded during GS-2.1 and the second part of GI-1 (between 180 to 120 cm and 30 to 15 cm), whilst lowest values relate to the first part

of GI-1, corresponding to BFA.f. Decreases in the mean@silt silt characterize the transitions between GS-2.1/GI-1 and GI-1/HC, respectively corresponding to the transitions BFA.A/BFA.B and BFA.e/BFA.c (Figure 5). An increase that takes place at 180 cm, mirroring the trend noticed in the mean@siliciclastic fraction, can be related to the change in
<caption>**TABLE 4** | List of the species contributing to the dissimilarity between clusters A and B, which separate at 47% similarity on the hierarchical cluster analysis (Figure 5).</caption>

Clusters A and B

<table><caption>Average dissimilarity = 63.87</caption><thead><tr><th>Species</th><th>Cluster B<br>Av. Abund</th><th>Cluster A<br>Av. Abund</th><th>Av. Diss</th><th>Diss/SD</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>14.09</td><td>29.62</td><td>9.33</td><td>1.16</td><td>14.61</td><td>14.61</td></tr><tr><td>Astrononion stelligerum</td><td>1.19</td><td>11.57</td><td>5.22</td><td>2.09</td><td>8.17</td><td>22.79</td></tr><tr><td>Trifarina angulosa</td><td>1.51</td><td>8.66</td><td>3.58</td><td>2.79</td><td>5.60</td><td>28.39</td></tr><tr><td>Discanomalina coronata</td><td>0.08</td><td>6.18</td><td>3.05</td><td>1.98</td><td>4.77</td><td>33.16</td></tr><tr><td>Rosalina brady</td><td>5.82</td><td>0.18</td><td>2.82</td><td>1.68</td><td>4.42</td><td>37.58</td></tr><tr><td>Chilostomella oolina</td><td>5.63</td><td>0.03</td><td>2.80</td><td>1.41</td><td>4.39</td><td>41.97</td></tr><tr><td>Bulimina aculeata</td><td>5.70</td><td>0.72</td><td>2.56</td><td>0.89</td><td>4.00</td><td>45.97</td></tr><tr><td>Cassidulina laevigata</td><td>10.05</td><td>6.56</td><td>2.48</td><td>1.27</td><td>3.89</td><td>49.86</td></tr><tr><td>Cibicides refulgens</td><td>0.21</td><td>3.43</td><td>1.62</td><td>2.23</td><td>2.54</td><td>52.40</td></tr><tr><td>Bolivina spathulata</td><td>3.41</td><td>0.21</td><td>1.61</td><td>1.13</td><td>2.52</td><td>54.92</td></tr><tr><td>Gavelinopsis praegeri</td><td>4.11</td><td>1.04</td><td>1.57</td><td>1.45</td><td>2.45</td><td>57.37</td></tr><tr><td>Bulimina striata</td><td>3.01</td><td>0</td><td>1.50</td><td>1.38</td><td>2.35</td><td>59.73</td></tr><tr><td>Miliolinella subrotunda</td><td>6.21</td><td>4.00</td><td>1.49</td><td>1.33</td><td>2.34</td><td>62.07</td></tr><tr><td>Hyalinea balthica</td><td>1.08</td><td>3.66</td><td>1.45</td><td>1.15</td><td>2.27</td><td>64.33</td></tr><tr><td>Spliroplectammina wrightii</td><td>2.65</td><td>0.60</td><td>1.30</td><td>1.04</td><td>2.04</td><td>66.37</td></tr><tr><td>Bolivina striatula</td><td>2.67</td><td>0.26</td><td>1.26</td><td>0.90</td><td>1.98</td><td>68.35</td></tr><tr><td>Bulimina marginata</td><td>2.25</td><td>0</td><td>1.12</td><td>0.90</td><td>1.76</td><td>70.11</td></tr></tbody></table>

The average abundance in each cluster, average dissimilarity, average dissimilarity/standard deviation ratio (Diss/SD), contribution (%) and cumulative contribution to the total dissimilarity are indicated for each species. A high Diss/SD ratio indicates that the considered species is a good discriminating species. The average dissimilarity between the two clusters is also given.

benthic foraminiferal assemblage (BFA e to f) that separates the Early and Late GI-1.

## Organic Geochemistry

Total organic carbon (TOC) varies between 0.18 wt% during GS-2.1 (460 cm) and 0.84 wt% at the transition between GI-1 and the Holocene (20 cm; **Figure 5**). Lowest values are recorded during GS-2.1 (0.18 wt%) and during the Holocene (0.34 and 0.37 wt%). The transition from GS-2.1 to GI-1 is characterized by an increase followed by a gradual decrease in TOC from 0.77 to 0.44 wt% (**Figure 5**). Low TOC values (ca. 0.6 wt%) are recorded for the Early phase of GI-1. It is followed by a gradual increase to higher TOC values, from ca. 0.6 to 0.8 wt%, toward the end of GI-1 (between 280 and 210 cm; **Figure 5**). This slight increase in TOC matches the first coarsening phase noticed in the sorted silt and siliciclastic fraction mean grain size, together with the change in benthic foraminiferal assemblage (**Figure 5**). Total organic carbon content remains unchanged and is relatively high (ca. 0.8%) until the end of GI-1.

The HI vs. Tmax diagram provides information about the origin of the organic matter and its relation to changes in benthic foraminiferal assemblages (**Figure 6**). All samples, except for two outliers (15 and 460 cm), are restricted to HI and Tmax values situated respectively between 123 and 286 mgHC/g TOC), and 409 and 428°C (**Figure 6B**). All samples, except sample 15, plot in the mixed marine-terrestrial/terrestrial deltaic domain (**Figure 6A**). The Greenland Stadial 2.1 samples, corresponding to BFA.A, are associated to the highest HI values (248 to 286 mgHC/g TOC) and are associated to the most marine

related organic matter. The Late GI-1 (corresponding to BFA.e) and the Early GI-1 (corresponding to BFA.f) samples relate to respectively intermediate (188 to 270 mgHC/g TOC) and low (123 to 229 mgHC/g TOC) HI values (**Figure 6A**). The Early GI-1 is associated with more terrestrial deltaic organic matter, whereas the Late GI-1 is rather affiliated to more mixed terrestrial/marine organic matter (**Figure 6A**). Higher Tmax (426 and 428°C) and intermediate (173 and 251 mgHC/g TOC) HI values characterize the Holocene (corresponding to BFA.c). Thus the organic matter during the Holocene appears to have a rather mixed terrestrial/marine origin, although the number of samples analyzed does not allow any precise interpretation.

## DISCUSSION

### Greenland Stadial 2.1

The Greenland Stadial 2.1 benthic foraminiferal assemblage (BFA.A) resembles the glacial assemblage described by Stalder et al. (2018) at BRI, with high abundances of the epibenthic *C. lobatulus* and *D. coronata* and infaunal *C. laevigata* (**Figures 3A,B** and **Table 3**). The filter feeding *D. coronata* and *C. lobatulus* preferentially cling to elevated substrata (e.g., coral branches, dropstones, hydroids) in well-oxygenated waters dominated by vigorous bottom currents (Linke and Lutze, 1993; Schönfeld, 1997; Margreth et al., 2009; Schönfeld et al., 2011). Together with the infaunal *C. laevigata* and *H. balthica*, they are typically associated with the availability of fresh, high quality marine organic matter (De Rijk et al., 2000; Schmiedl et al., 2000;
<table><caption>**TABLE 5** | List of the 5 most contributing species to the similarity within clusters b, c, d, e and f, which separate at 62% similarity on the hierarchical cluster analysis (Figure 5).</caption><thead><tr><th colspan="5">Average similarity: 65.11</th></tr><tr><th>Species</th><th>Av.Abund</th><th>Av.Sim</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>31.06</td><td>18.63</td><td>28.61</td><td>28.61</td></tr><tr><td>Astrononion stelligerum</td><td>13.22</td><td>11.02</td><td>16.92</td><td>45.53</td></tr><tr><td>Trifarina angulosa</td><td>8.68</td><td>6.97</td><td>10.7</td><td>56.24</td></tr><tr><td>Cassidulina laevigata</td><td>6.14</td><td>4.82</td><td>7.4</td><td>63.64</td></tr><tr><td>Discanomalina coronata</td><td>5.51</td><td>3.65</td><td>5.61</td><td>69.25</td></tr></tbody><thead><tr><th colspan="5">Average similarity: 69.51</th></tr><tr><th>Species</th><th>Av.Abund</th><th>Av.Sim</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Globocassidulina subglobosa</td><td>11.16</td><td>9.19</td><td>13.22</td><td>13.22</td></tr><tr><td>Cibicides lobatulus</td><td>10.6</td><td>8.99</td><td>12.93</td><td>26.15</td></tr><tr><td>Cassidulina laevigata</td><td>7.18</td><td>6.46</td><td>9.29</td><td>35.44</td></tr><tr><td>Uvigerina mediterranea</td><td>5.93</td><td>3.85</td><td>5.54</td><td>40.97</td></tr><tr><td>Miliolinella subrotunda</td><td>4.79</td><td>3.83</td><td>5.51</td><td>46.48</td></tr></tbody><thead><tr><th colspan="5">Average similarity: 70.25</th></tr><tr><th>Species</th><th>Av.Abund</th><th>Av.Sim</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>20.19</td><td>18.22</td><td>25.93</td><td>25.93</td></tr><tr><td>Gavelinopsis praegei</td><td>8.77</td><td>8.31</td><td>11.83</td><td>37.76</td></tr><tr><td>Miliolinella subrotunda</td><td>9.19</td><td>7.74</td><td>11.01</td><td>48.77</td></tr><tr><td>Spliroplectammina wrightii</td><td>8.46</td><td>6.86</td><td>9.76</td><td>58.53</td></tr><tr><td>Cassidulina laevigata</td><td>7.64</td><td>6.37</td><td>9.06</td><td>67.59</td></tr></tbody><thead><tr><th colspan="5">Average similarity: 66.73</th></tr><tr><th>Species</th><th>Av.Abund</th><th>Av.Sim</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>17.09</td><td>13.04</td><td>19.55</td><td>19.55</td></tr><tr><td>Cassidulina laevigata</td><td>13.82</td><td>11.94</td><td>17.9</td><td>37.44</td></tr><tr><td>Bolivina spathulata</td><td>5.32</td><td>3.8</td><td>5.69</td><td>43.13</td></tr><tr><td>Chilostomella oolina</td><td>5.85</td><td>3.65</td><td>5.47</td><td>48.6</td></tr><tr><td>Bulimina aculeata</td><td>5.6</td><td>3.46</td><td>5.18</td><td>53.78</td></tr></tbody><thead><tr><th colspan="5">Average similarity: 65.27</th></tr><tr><th>Species</th><th>Av.Abund</th><th>Av.Sim</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>10.74</td><td>8.56</td><td>13.12</td><td>13.12</td></tr><tr><td>Rosalina brady</td><td>8.38</td><td>6.94</td><td>10.64</td><td>23.76</td></tr><tr><td>Miliolinella subrotunda</td><td>7.74</td><td>6.35</td><td>9.73</td><td>33.49</td></tr><tr><td>Cassidulina laevigata</td><td>6.47</td><td>4.83</td><td>7.4</td><td>40.89</td></tr><tr><td>Chilostomella oolina</td><td>6.17</td><td>4.2</td><td>6.44</td><td>47.33</td></tr></tbody></table>

The average abundance, similarity, contribution (%) and cumulative contribution (%) are indicated for each species. The average similarity of the cluster is also given.
<table><caption>**TABLE 6** | List of the species contributing to the dissimilarity between clusters b, c, d, e and f, which separate at 47% similarity on the hierarchical cluster analysis (Figure 5).</caption><caption>Clusters d and b</caption><thead><tr><th rowspan="2">Species</th><th colspan="6">Average dissimilarity = 54.35</th></tr><tr><th>Cluster d<br>Av.Abund</th><th>Cluster b<br>Av.Abund</th><th>Av.Diss</th><th>Diss/SD</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>20.19</td><td>31.06</td><td>9.19</td><td>1.46</td><td>16.92</td><td>16.92</td></tr><tr><td>Astrononion stelligerum</td><td>3.25</td><td>13.22</td><td>4.99</td><td>2.74</td><td>9.17</td><td>26.09</td></tr><tr><td>Spliroplectammina wrightii</td><td>8.46</td><td>0.06</td><td>4.20</td><td>3.94</td><td>7.73</td><td>33.82</td></tr><tr><td>Gavelinopsis praegeri</td><td>8.77</td><td>0.79</td><td>3.99</td><td>10.06</td><td>7.34</td><td>41.16</td></tr><tr><td>Trifarina angulosa</td><td>3.28</td><td>8.68</td><td>2.84</td><td>1.76</td><td>5.22</td><td>46.38</td></tr><tr><td>Rosalina brady</td><td>5.61</td><td>0.08</td><td>2.77</td><td>2.04</td><td>5.09</td><td>51.46</td></tr><tr><td>Miliolinella subrotunda</td><td>9.19</td><td>4.16</td><td>2.52</td><td>2.31</td><td>4.63</td><td>56.10</td></tr><tr><td>Discanomalina coronata</td><td>0.72</td><td>5.51</td><td>2.39</td><td>1.65</td><td>4.40</td><td>60.50</td></tr><tr><td>Globocassidulina subglobosa</td><td>4.94</td><td>1.38</td><td>1.78</td><td>1.65</td><td>3.28</td><td>63.78</td></tr><tr><td>Hyalinea balthica</td><td>0.60</td><td>3.99</td><td>1.73</td><td>1.26</td><td>3.19</td><td>66.97</td></tr><tr><td>Cibicides refulgens</td><td>0.98</td><td>3.45</td><td>1.34</td><td>1.72</td><td>2.46</td><td>69.43</td></tr><tr><td>Cassidulina laevigata</td><td>7.64</td><td>6.14</td><td>1.21</td><td>1.22</td><td>2.22</td><td>71.65</td></tr></tbody></table>

<table><caption>Clusters f and d</caption><thead><tr><th rowspan="2">Species</th><th colspan="6">Average dissimilarity = 45.93</th></tr><tr><th>Cluster f<br>Av.Abund</th><th>Cluster d<br>Av.Abund</th><th>Av.Diss</th><th>Diss/SD</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>10.74</td><td>20.19</td><td>4.79</td><td>2.17</td><td>10.44</td><td>10.44</td></tr><tr><td>Chilostomella oolina</td><td>6.17</td><td>0.47</td><td>2.85</td><td>1.59</td><td>6.20</td><td>16.64</td></tr><tr><td>Bulimina aculeata</td><td>6.39</td><td>1.56</td><td>2.45</td><td>0.66</td><td>5.34</td><td>21.98</td></tr><tr><td>Spliroplectammina wrightii</td><td>4.35</td><td>8.46</td><td>2.23</td><td>1.76</td><td>4.86</td><td>26.84</td></tr><tr><td>Bulimina marginata</td><td>4.22</td><td>0.00</td><td>2.11</td><td>1.74</td><td>4.59</td><td>31.43</td></tr><tr><td>Bulimina striata</td><td>4.15</td><td>0.00</td><td>2.07</td><td>2.36</td><td>4.51</td><td>35.94</td></tr><tr><td>Gavelinopsis praegeri</td><td>5.10</td><td>8.77</td><td>1.85</td><td>1.85</td><td>4.02</td><td>39.96</td></tr><tr><td>Rosalina brady</td><td>8.38</td><td>5.61</td><td>1.83</td><td>1.28</td><td>3.98</td><td>43.94</td></tr><tr><td>Globocassidulina subglobosa</td><td>2.36</td><td>4.94</td><td>1.54</td><td>1.48</td><td>3.34</td><td>47.28</td></tr><tr><td>Cassidulina laevigata</td><td>6.47</td><td>7.64</td><td>1.50</td><td>1.27</td><td>3.26</td><td>50.54</td></tr><tr><td>Miliolinella subrotunda</td><td>7.74</td><td>9.19</td><td>1.34</td><td>1.38</td><td>2.92</td><td>53.46</td></tr><tr><td>Trifarina angulosa</td><td>1.49</td><td>3.28</td><td>1.25</td><td>1.15</td><td>2.71</td><td>56.17</td></tr></tbody></table>

<table><caption>Clusters e and f</caption><thead><tr><th rowspan="2">Species</th><th colspan="6">Average dissimilarity = 47.06</th></tr><tr><th>Cluster e<br>Av.Abund</th><th>Cluster f<br>Av.Abund</th><th>Av.Diss</th><th>Diss/SD</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Cibicides lobatulus</td><td>17.09</td><td>10.74</td><td>4.02</td><td>1.26</td><td>8.54</td><td>8.54</td></tr><tr><td>Cassidulina laevigata</td><td>13.82</td><td>6.47</td><td>3.87</td><td>1.99</td><td>8.22</td><td>16.76</td></tr><tr><td>Rosalina brady</td><td>3.59</td><td>8.38</td><td>2.51</td><td>1.65</td><td>5.33</td><td>22.09</td></tr><tr><td>Bulimina aculeata</td><td>5.60</td><td>6.39</td><td>2.49</td><td>0.75</td><td>5.29</td><td>27.38</td></tr><tr><td>Bolivina striatula</td><td>4.79</td><td>0.81</td><td>2.13</td><td>1.65</td><td>4.52</td><td>31.89</td></tr><tr><td>Chilostomella oolina</td><td>5.85</td><td>6.17</td><td>2.05</td><td>1.16</td><td>4.37</td><td>36.26</td></tr><tr><td>Bolivina spathulata</td><td>5.32</td><td>1.70</td><td>2.00</td><td>1.49</td><td>4.24</td><td>40.50</td></tr><tr><td>Spliroplectammina wrightii</td><td>0.62</td><td>4.35</td><td>1.88</td><td>1.64</td><td>4.00</td><td>44.50</td></tr><tr><td>Miliolinella subrotunda</td><td>4.64</td><td>7.74</td><td>1.87</td><td>1.48</td><td>3.97</td><td>48.48</td></tr><tr><td>Bulimina marginata</td><td>0.54</td><td>4.22</td><td>1.86</td><td>1.52</td><td>3.96</td><td>52.44</td></tr><tr><td>Bulimina striata</td><td>2.18</td><td>4.15</td><td>1.41</td><td>1.56</td><td>3.00</td><td>55.44</td></tr><tr><td>Gavelinopsis praegeri</td><td>2.83</td><td>5.10</td><td>1.33</td><td>1.46</td><td>2.82</td><td>58.26</td></tr><tr><td>Cassidulinoides bradyi</td><td>2.29</td><td>0.02</td><td>1.14</td><td>1.23</td><td>2.42</td><td>60.68</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 6 | Continued</caption><caption>Clusters c and e<br>Average dissimilarity = 51.77</caption><thead><tr><th>Species</th><th>Cluster c<br>Av. Abund</th><th>Cluster e<br>Av. Abund</th><th>Av. Diss</th><th>Diss/SD</th><th>Contrib%</th><th>Cum.%</th></tr></thead><tbody><tr><td>Globocassidulina subglobosa</td><td>11.16</td><td>3.53</td><td>3.81</td><td>2.48</td><td>7.37</td><td>7.37</td></tr><tr><td>Cibicides lobatulus</td><td>10.6</td><td>17.09</td><td>3.79</td><td>1.21</td><td>7.33</td><td>14.69</td></tr><tr><td>Cassidulina laevigata</td><td>7.18</td><td>13.82</td><td>3.31</td><td>1.97</td><td>6.4</td><td>21.09</td></tr><tr><td>Chilostomella oolina</td><td>0.34</td><td>5.85</td><td>2.75</td><td>1.34</td><td>5.31</td><td>26.41</td></tr><tr><td>Uvigerina mediterranea</td><td>5.93</td><td>1.21</td><td>2.36</td><td>1.58</td><td>4.55</td><td>30.96</td></tr><tr><td>Bolivina striatula</td><td>0.10</td><td>4.79</td><td>2.34</td><td>1.83</td><td>4.52</td><td>35.48</td></tr><tr><td>Bulimina aculeata</td><td>1.87</td><td>5.6</td><td>2.02</td><td>1.16</td><td>3.90</td><td>39.38</td></tr><tr><td>Bolivina spathulata</td><td>1.98</td><td>5.32</td><td>1.84</td><td>1.40</td><td>3.56</td><td>42.94</td></tr><tr><td>Nutallides umbonifera</td><td>3.01</td><td>0.00</td><td>1.5</td><td>34.45</td><td>2.90</td><td>45.84</td></tr><tr><td>Non-ionella turgida</td><td>0.23</td><td>2.71</td><td>1.24</td><td>2.11</td><td>2.39</td><td>48.23</td></tr><tr><td>Miliolinella subrotunda</td><td>4.79</td><td>4.64</td><td>1.04</td><td>1.40</td><td>2.01</td><td>50.24</td></tr><tr><td>Spliroplectammina wrightii</td><td>2.66</td><td>0.62</td><td>1.03</td><td>2.62</td><td>1.99</td><td>52.23</td></tr><tr><td>Gavelinopsis praegeri</td><td>4.29</td><td>2.83</td><td>0.94</td><td>1.54</td><td>1.82</td><td>54.05</td></tr><tr><td>Cassidulinoides bradyi</td><td>0.51</td><td>2.29</td><td>0.94</td><td>1.07</td><td>1.82</td><td>55.87</td></tr><tr><td>Rosalina brady</td><td>3.03</td><td>3.59</td><td>0.91</td><td>1.35</td><td>1.76</td><td>57.63</td></tr><tr><td>Bulimina marginata</td><td>2.13</td><td>0.54</td><td>0.85</td><td>2.39</td><td>1.63</td><td>59.26</td></tr></tbody></table>

The average abundance in each cluster, average dissimilarity, average dissimilarity/standard deviation ratio (Diss/SD), contribution (%) and cumulative contribution to the total dissimilarity are indicated for each species. A high Diss/SD ratio indicates that the considered species is a good discriminating species. The average dissimilarity between each pair of clusters is also given.

Milker et al., 2009; Stalder et al., 2018). This assemblage, together with the HI vs. Tmax diagram, confirms that the seafloor at BRI received the input of high quality marine derived organic matter during GS-2.1 (Figure 6). Surface productivity was important during GS-2.1, as suggested by the high planktonic δ¹³C values (Figure 5). Based on Al-normalized elemental ratios, Feenstra (2020) suggest that arid continental conditions prevailed at the end of the last glacial period and that eolian dust was an important fertilizer for marine communities during GS-2.1.

The presence of vigorous bottom currents is directly recorded by the high sortable silt mean grain size (Figure 5). These increased bottom currents would have been favorable for the large epibenthic foraminifera *C. lobatulus* and *D. coronata*, and the infaunal *T. angulosa* (Figures 3A,B and Table 3). *Trifarina angulosa* is reported from current-swept areas and is resistant to permanent winnowing (Mackensen et al., 1995; Schönfeld, 2002). The association of *C. lobatulus*, *D. coronata*, *T. angulosa*, *Cassidulina* spp. and CWCs is reported from the Porcupine Seabight (Margreth et al., 2009; Schönfeld et al., 2011; Fentimen et al., 2018), Norwegian margin (Spezzaferri et al., 2013; Stalder et al., 2014) and the Alboran Sea (Margreth et al., 2011; Stalder et al., 2015, 2018). However, this assemblage is here observed in relation to a thick accumulation of the bryozoan *Buskea dichotoma* associated to the brachiopod *G. vitreus*, whilst CWCs are absent (Figure 2). Thus, despite seemingly favorable conditions such as high surface productivity and vigorous bottom currents, CWCs did not thrive during GS-2.1 (Figure 2). Bryozoans, CWCs and *G. vitreus* share the same need for sustained bottom currents, although bryozoans and *G. vitreus* prefer moderate velocities and do not cope with excessively

strong bottom currents (Emig and Arnaud, 1988; Scholz and Hillmer, 1995; Bjerager and Surlyk, 2007; Roberts et al., 2009). Likewise, bryozoans and CWCs thrive when fresh organic matter is available, although CWCs are more tolerant to varying food sources than bryozoans (Best and Thorpe, 1994; Holbourn et al., 2002; Duineveld et al., 2004; van Oevelen et al., 2016). Lower temperatures coupled to lower sea level stand and increased upwelling promoted the development of bryozoan reef mounds at the Great Australian Bight (James et al., 2000; Holbourn et al., 2002). Similar cold conditions and low sea level stand during GS-2.1 at BRI are suggested by the high planktonic δ¹⁸O values and the presence of the brachiopod *G. vitreus* (Figure 2). Indeed, *G. vitreus* is found today along the Mediterranean continental margin, at depths between ca. 100 and 250 m (Emig, 1988). It is associated to “white coral” debris between 235 and 255 m depth near the Hyères Islands (SE France), where it reaches densities between 5 and 10 ind.m⁻² (Emig and Arnaud, 1988). Below this depth range, *G. vitreus* becomes rarer and quickly disappears from the faunal community (Emig and Arnaud, 1988). Thus, these observations would suggest that the sea level during GS-2.1 was 75 ± 10 m lower than today. Additionally, *A. stelligerum* and *H. balthica*, two foraminifera species affiliated to the GS-2.1 assemblage (Figures 3A,B), are reported to prefer cold conditions (Husum and Hald, 2004; Murray, 2006). Hence, the presence of bryozoans and the absence of corals during GS-2.1 may be partly induced by lower temperatures at BRI. Overall, the environment during GS-2.1 is characterized by cold temperatures, high surface productivity and relatively vigorous but not excessively strong bottom currents. Increased bottom current velocity during glacial periods in the Western Mediterranean has been linked
<figure><img src="image_9.png" /><figcaption>**FIGURE 5** | Planktonic (G. *bulloides*) and benthic (C. *lobatulus*) δ¹³C values, Total Organic Carbon (%), mean grain size of bulk material (µm), siliciclastic (µm) and sorted silt fraction (corresponding to the 10 to 63 µm siliciclastic grain size range, expressed in µm). The benthic foraminiferal assemblages (BFA) identified at 47% similarity and at 62% similarity thanks to the hierarchical cluster analysis (Figure 4) are illustrated. The stratigraphic column defined in Figure 2, the G. *bulloides* planktonic and C. *lobatulus* benthic δ¹⁸O values and the U/Th coral and calibrated ¹⁴C benthic foraminifera ages are provided as supporting information.</figcaption></figure>

to a stronger intensity of the Levantine Intermediate Water (Rogerson et al., 2008; Toucanne et al., 2012). We suggest that such increased intermediate water flow occurred during the cold GS-2.1 at BRI (Figure 7). Moreover, Cacho et al. (2006) reported that deep water density and ventilation reached a maximum during MIS 2. Recently formed dense waters in the Mediterranean may lie at intermediate depths, above 1500 m, hence never reaching the deepest parts of the western Mediterranean (Sparnocchia et al., 1995; Millot, 1999; Ercilla et al., 2016). In addition, WMDW has been traced at depths shallower than 500 m along the Moroccan margin (Ercilla et al., 2016) and contributes significantly to the intermediate waters in the region. Thus, it is conceivable that shoaling of WMDW along the Moroccan margin during the last glacial period, especially during GS-2.1, led to a higher admixture of this water mass to overlying intermediate waters. The contribution of these nutrient-rich and well-ventilated waters at intermediate depths, together with an important input of eolian dust, would

have favored the proliferation of bryozoans and associated epibenthic foraminifera.

## The Early Greenland Interstadial 1
### Transition Between Greenland Stadial 2.1 and Greenland Interstadial 1

The transition from GS-2.1 to GI-1 is marked by a significant change in the benthic macro and microfauna (Figures 2, 3). The larger epibenthic *C. lobatulus* and *D. coronata* retreat, together with the infaunal *T. angulosa* and the cold water related *H. balthica* and *A. stelligerum*. Simultaneously, CWCs become the dominant macro-organism, to the detriment of bryozoans (Figure 2). The benthic foraminiferal assemblage at the transition between GS-2.1 and GI-1 (BFA.d) is characterized by high abundances of the epibenthic agglutinated *S. wrightii*, together with the small epibenthic hyaline *G. praegeri* and *R. bradyi* (Tables 3, 4). *Spiroplectammina wrightii* has been associated to modern CWC habitats dominated by dead corals at the
<figure><img src="image_10.png" /></figure>

**FIGURE 6** | (A) Hydrogen Index (mgHC/g TOC) vs. Tmax (°C) diagram. (B) Close-up of the dashed rectangle illustrated in part A. The distribution of samples allows to identify the origin of the organic matter present in the sediment. The origin of organic matter becomes more terrestrial with decreasing HI values. Benthic foraminiferal assemblages defined by the hierarchical cluster analysis are overlaid. The color code used for samples and overlays corresponds to the one defined for benthic foraminiferal assemblages in **Figure 4**.
<figure><img src="image_11.png" /></figure>

**FIGURE 7** | Three-step diagram illustrating the environmental changes affecting Brittlestar Ridge I during Greenland Stadial 2.1 (GS-2.1) and Greenland Interstadial 1 (GI-1). Parts (A–C) are close-ups of the seafloor at Brittlestar Ridge I for each period. The illustrated foraminifera are selected based on the dissimilarity values provided by the Similarity Percentage analysis (Table 4). The selected foraminifera do not encompass the complete assemblage of each time period and are a simplified interpretation. The illustrated macrofauna and microfauna are not to scale.

Porcupine Seabight, NE Atlantic (Margreth et al., 2009). Its ecological requirements are not well constrained, however a related species, *S. biformis*, is known to tolerate severe dysoxia

and to respond positively to increased nutrient availability in the Norwegian Drammensfjord (Alve, 1991; Gross, 2000). Szymańska et al. (2017) found it to be the most abundant genus
in the Adventfjorden and Hornsund fjords of Svalbard, whilst it is also a common occupant of seep environments (Bernhard et al., 2001; McGann and Conrad, 2018). The epiphytic *R. bradyi* and *G. praegeri* have previously been observed together with CWCs at BRI (Stalder et al., 2015, 2018). *Rosalina bradyi* is frequently found along shallow shelf sites of the Mediterranean and dwells attached to elevated substrates, such as vegetation (Jorissen, 1987; Cimerman and Langer, 1991; Barras et al., 2014). It also shows a preference for warmer waters (Murray, 2006). Laboratory experiments have demonstrated that it responds positively to sporadic pulses of high amounts of organic matter and can tolerate extreme salinity changes and harsh oxygen depletion (Heinz et al., 2001; Fontanier et al., 2008). Thus, BFA.d indicates that the transition from GS-2.1 to GI-1 was marked by decreasing seafloor oxygenation and eutrophication. Rapidly decreasing bottom current velocities, inferred from the decrease in sortable silt mean grain size, would have further characterized the passage from GS-2.1 to GI-1. We link this decrease in bottom current velocity to a slowdown of intermediate water mass circulation (Figure 7). A substantial reduction in WMDW flow is observed in the deep western Mediterranean at 14.5 ka (Rogerson et al., 2008). We propose that the change in faunal composition and the abrupt decrease in sortable silt mean grain size at the transition between GS-2.1 and GI-1 is linked to this reduction of WMDW flow. A reduced WMDW flow would have weakened the ventilation of intermediate water masses. This reduced ventilation and current velocity at BRI may have in turn resulted in an important faunal turnover. These observations suggest that the decrease in the intensity of WMDW flow in the Southeastern Alboran Sea started as early as ca. 14.8 ka ago. Furthermore, the planktonic and benthic δ¹⁸O values support that surface and bottom waters also became warmer (Figure 2). This observation is supported by the disappearance of cold water related benthic foraminifera species (*A. stelligerum* and *H. balthica*) and the appearance of the warmer water dwelling *R. bradyi* (Figure 3).

## Increased Terrestrial Input During the Early Greenland Interstadial 1

The environmental changes that happened at the GS-2.1/GI-1 transition intensified during the Early GI-1. High abundances of the infaunal *Bulimina* spp. and *C. oolina*, together with increased *R. bradyi* abundances, characterize BFA.f (Figure 3 and Tables 5, 6). *Bulimina* spp. is a typical genus for dysoxic and eutrophic environments and is dominant in the Mediterranean Sea near the Po and Rhône river deltas (Phleger and Soutar, 1973; Lutze and Coulbourn, 1984; Jorissen, 1987; Schmiedl et al., 2000; Mojtahid et al., 2009). *Bulimina* spp. is adapted to feed on degraded organic matter (De Rijk et al., 2000; Schmiedl et al., 2000). The intermediate or deep infaunal *C. oolina* has been described as a species adapted to suboxic conditions and feeding preferentially on degraded organic matter (Corliss, 1985; Sen Gupta and Machain-Castillo, 1993; Bernhard and Sen Gupta, 1999; Schmiedl et al., 2000; Fontanier et al., 2002). Elevated organic matter concentrations in the deeper parts of the sediment offer favorable conditions for the proliferation of this species (Rathburn and Corliss, 1994). Thus, the association of *Bulimina*

spp., *C. oolina* and *R. bradyi* suggests that the eutrophication trend that began at the transition between GS-2.1 and GI-1 was reinforced throughout the Early GI-1. Moreover, BFA.f relates to a more terrestrial organic matter origin (Figure 6). Organic-rich layers (ORLs), presenting a similar assemblage composition as BFA.f, have previously been described in the deep Alboran basins (Comas et al., 1996; Cacho et al., 2002; Rogerson et al., 2008). Organic-rich layers boast high TOC levels (up to 0.9 wt%), although noticeably lower than sapropel deposits in the eastern Mediterranean (> 2 wt% TOC) (Emeis et al., 2000; Cacho et al., 2002; Martínez-Ruiz et al., 2003; Rogerson et al., 2008). The sediment during the Early GI-1 at BRI does not have higher TOC content compared to GS-2.1. However, based on the benthic foraminiferal assemblage, we propose that it can still be related to ORL1, although not representing an ORL deposit *per se*.

The warming that started at the transition between GS-2.1 and GI-1 persisted throughout the Early GI-1, as indicated by the planktonic and benthic δ¹⁸O values (Figure 2). Higher abundances of *R. bradyi* during the Early GI-1 confirm this trend (Figure 3). Permanent low bottom current velocity during the Early GI-1, implied by the stable low sortable silt mean grain size, suggests that the weakening of intermediate water flow persisted throughout the Early GI-1 (Figure 7). Furthermore, the lower abundance of the epibenthic *C. lobatulus*, *S. wrightii*, *G. praegeri* and the infaunal *T. angulosa* in BFA.f in comparison to BFA.d may even indicate that the intermediate water mass flow decreased during the Early GI-1 (Table 4). This would in turn suggest a persistent low or even dwindling WMDW flow.

## Environmental Parameters Triggering Coral Proliferation

The African continent experienced important climatic changes during the Late Pleistocene and Holocene (Jolly et al., 1998; Gasse, 2000). A warm and humid period, known as the African Humid Period, generally reported to have lasted from 11 to 5.5 ka, had a major impact on marine and terrestrial systems, such as faunal and floral communities, as well as human migration pathways (Barth, 1857; Kuper and Kröpelin, 2006; deMenocal and Tierney, 2012 and references therein). East African lakes were tens or hundreds of meters higher than today, whilst large river systems flowed through the Western Sahara (Talbot and Livingstone, 1989; Skonieczny et al., 2015). The abrupt change from cold hyper-arid to warm humid conditions began as soon as 14.8 ka BP off Cap Blanc, Mauritania (deMenocal et al., 2000). The composition of BFA.f (Tables 3, 4) and the type of organic matter present in the sediment (Figure 6) could suggest that the East Melilla Coral Province experienced an increase in fluvial input during the GS-2.1/GI-1 transition and the Early GI-1. Moreover, the coral fragment dated at the transition GS-2.1/GI-1 has an age of 14.814 ± 0.052 ka (Table 1 and Figure 2). This age can be considered near to concomitant with the beginning of the African Humid Period as defined by deMenocal et al. (2000) at IODP Site 658C and the beginning of ORL1 (ca. 14.800 ka; Cacho et al., 2002; Jimenez-Espejo et al., 2008; Rodrigo-Gámiz et al., 2011; Martinez-Ruiz et al., 2015). Thus, this would confirm that the switch from cold arid to warm humid conditions over the North African continent, related to the beginning of the
African Humid Period, led to a sharp increase in fluvial input at BRI. At the same time hydrodynamic changes were taking place in the Mediterranean Sea. The contribution of MAW in the Alboran Sea increases with rising sea level and deepening of the Gibraltar sill (Sierro et al., 2005). Simultaneously, changes in deep-water mass dynamics occurred. Rogerson et al. (2019) proposed that during the last deglacial, dense waters formed in the Gulf of Lions were injected at intermediate depths rather than at deeper depths. This would have isolated water masses above 2000 m from water bodies situated beneath this limit (Rogerson et al., 2019). Following the observations made by Ercilla et al. (2016), we suggest that this boundary may have been situated at shallower depths along the Moroccan margin, separating BRI from the influence of deeper well-ventilated waters. The important freshwater inflow from the North African continent, combined to the injection of more MAW in the Alboran Sea, further resulted in weak intermediate and deep-water ventilation and increased stratification (Sierro et al., 2005; Rogerson et al., 2008; Toucanne et al., 2012). Thus, the equilibrium at BRI appears to shift between a system primarily influenced by exchanges between intermediate and deep water masses during GS-2.1 to a system dominated by surface water dynamics and separated from the influence of deep water masses during the Early GI-1. The combination of rapidly increasing fluvial input, temperature and sea level rise, and the slowdown of intermediate and deep water mass flow, triggered a change in benthic habitats and the advent of a thriving CWC ecosystem (Figure 7). Weaker bottom water currents, combined to the high fluvial input and the presence of a dense coral framework, led to dysoxic conditions at the sediment/water interface at BRI (Figure 7). These observations confirm that CWCs can withstand and even thrive in weakly ventilated bottom waters if sufficient food is available.

## The Late Glacial Interstadial 1

The Late GI-1 assemblage (BFA.e), in comparison to the Early GI-1 assemblage (BFA.f), is characterized by higher *C. bradyi*, *C. laevigata*, *C. lobatulus*, *G. subglobosa* and Bolivinids. abundances and lower *Bulimina* spp., *G. praegegi*, *R. bradyi* and *S. wrightii* abundances (Figures 3, 4 and Tables 5, 6). The opportunistic *Bolivina* spp. has an infaunal living strategy and dwells in highly productive areas associated to low seafloor oxygenation (Lutze and Coulbourn, 1984; Mackensen et al., 1995; Schmiedl et al., 2000; Murray, 2006). The association of *Bolivina* spp. and *C. oolina* suggests that the seafloor during the Late GI-1 experienced periods of oxygen depletion. However, in contrast with the Early GI-1, higher abundances of the epibenthic *C. lobatulus* and the infaunal *C. laevigata*, *C. bradyi* and *G. subglobosa* indicate that the seafloor was overall better oxygenated during the Late GI-1. The species *G. subglobosa* is regularly found along the Alboran platform and prefers high energy and well-oxygenated environments (Mackensen et al., 1995; Milker et al., 2009). It responds favorably and quickly to the input of fresh phytodetritus (Gooday, 1993; Suhr et al., 2003; Sun et al., 2006). It can be inferred from the higher *C. lobatulus* and *G. subglobosa* abundances and increased sortable silt mean grain size that bottom currents picked up during the Late GI-1 (Figures 3, 6). We propose that these observations can be related

to an increase in velocity of intermediate water masses during the later phase of the GI-1, together with a re-establishment of deep water circulation in the Southeast Alboran Sea (Figure 7). In parallel to these changes taking place at intermediate and possibly deep water depths, surface waters would also have been affected by varying conditions. The origin of the organic matter changes from more terrestrial to more marine derived (Figure 6). The higher contribution of phytodetritus feeding benthic foraminifera, such as *C. lobatulus* and *G. subglobosa*, confirms that BRI received less terrestrial material and was experiencing an increased input of fresher marine derived organic matter. This change in the origin of the organic matter reaching the seafloor and in the benthic foraminiferal community indicates that fluvial input was overall reduced during the Late GI-1 (Figure 7). Such a reduction of fluvial input may be linked to a transition toward less humid and more arid conditions above the continent. We propose that the Late GI-1, as observed in core MD13-3455G, could represent a precursor phase to the arid conditions that prevailed during the Younger-Dryas in the region (Fink et al., 2015; Stalder et al., 2015). Furthermore, the slightly higher planktonic and benthic $\delta^{18}$O values during the Late GI-1 may suggest a minor cooling of surface and intermediate water masses. Coral proliferation phases are less sustained during the Late GI-1 in comparison to the Early GI-1 (Figure 2). Thus, the weaker fluvial input at BRI during the Late GI-1 had a negative effect on the coral community, although it did not lead to a full retreat (Figure 7). Moreover, reduced fluvial input in the Southeast Alboran Sea may have favored the re-establishment of intermediate and deep water mass circulation. It can be inferred from these observations that intermediate and deep water mass flows reached a minimum from ca. 14.8 to 14.3 ka, whilst it picked up after 14.3 ka (Figure 7).

Environmental conditions during the Late GI-1 were not stable. A short interval during the middle of the Late GI-1 (ca. 120 to 60 cm) corresponds to a decrease in sortable silt mean grain size and higher abundance of *C. oolina* (Figures 3, 6), thus implying a slowdown in bottom current velocity and an increased input of degraded organic matter. We propose that this time interval is associated to a short period of increased humidity over the continent and to reduced intermediate and deep water mass flow, similar to the Early GI-1, although of a shorter duration.

## The Holocene

The Holocene is known to show significant short-term climatic variability in the Alboran Sea (Bond et al., 2001; Frigola et al., 2008; Ausín et al., 2015; Fink et al., 2015). However, the reduced size of the Holocene horizon in core MD13-3455G does not allow a precise interpretation of changes having occurred during this period. Thus, we can only track the general environmental trends that characterized the Holocene at BRI. These would correspond to major rearrangements in water masses, which led to important changes in benthic communities. Rapidly increasing surface and bottom water temperatures during the Holocene are inferred from the decrease in planktonic and benthic $\delta^{18}$O values and lower *H. balthica* abundance compared to the Late GI-1 (Figures 2, 3). This matches previous observations made in the Alboran
Sea (Cacho et al., 2001; Ausín et al., 2015). This general warming at shallow and intermediate water depths is associated to a shift in the coral community, *M. oculata* becoming the dominant species instead of *D. pertusum*, and in the benthic foraminiferal assemblage (Figure 2). This change from *D. pertusum* to *M. oculata* dominated environments is documented at a number of different areas in the East and West Melilla Coral Provinces (Fink et al., 2013, 2015; Stalder et al., 2015, 2018; Wang et al., 2019; Feenstra, 2020), and is hence a widespread phenomenon. The decrease in portable silt mean grain size together with the retreat of the epibenthic foraminifera species *C. lobatulus* indicates that the seafloor experienced a decrease in bottom current strength during the Holocene. We relate this decrease in bottom water currents, comparable to the passage from GS-2.1 to the Early GI-1, to a reduced flow of intermediate water masses. Reduced rates of Levantine Intermediate Water flow into the Alboran Sea during interglacials have previously been documented in the Corsica Through and in the Sicily Channel (Incarbona et al., 2011; Toucanne et al., 2012). In contrast, the influence of MAW during interglacials increases as a result of sea level rise (Sierro et al., 2005). The increase in MAW flow could explain higher abundances of the phytodetritus feeding infaunal *G. subglobosa* and *U. mediterranea*. The species *U. mediterranea* is considered opportunistic in the Mediterranean Sea, responding positively to the important supply of fresh organic matter (De Rijk et al., 2000; Schmiedl et al., 2000; Fontanier et al., 2006). Thus, increasing MAW inflow would favor surface productivity during the Holocene, in turn promoting the development of phytodetritus feeding foraminifera at the seafloor.

## CONCLUSION

A range of measurements made on the MD13-3455G gravity core demonstrate that a combination of short-term fluctuations in continental and marine systems control coral development and demise at the East Melilla Coral Province. Furthermore, it is likely that fluvial input and nutrient supply are the primary environmental factor driving cold-water coral proliferation in the Southeast Alboran Sea, whilst they may influence the dynamics of intermediate and deep water masses in the region. An important change in benthic communities took place at the transition from Greenland Stadial 2.1 to Greenland Interstadial 1. Cold conditions, lower sea-level and stronger bottom currents during Greenland Stadial 2.1 would have favored the development of an ecosystem dominated by the bryozoan *Buskea dichotoma* associated to the brachiopod *Gryphus vitreus*. These environmental conditions were tightly linked to enhanced intermediate and deep water mass flow. Rapidly increasing fluvial input and surface and bottom water temperatures, coupled to a weakening of bottom currents, led to the decline of bryozoans and their replacement by cold-water corals during the Early Greenland Interstadial 1. Changes in the benthic foraminiferal fauna, such as the lower abundances of large epibenthic (e.g., *C. lobatulus* and

*D. coronata*) and increased numbers of infaunal species (e.g., Buliminids and *C. oolina*), also reflect the environmental changes from Greenland Stadial 2.1 to Greenland Interstadial 1. We propose that this switch is related to the beginning of the African Humid Period, establishing more humid and warmer conditions over the continent and initiating continental runoff. Simultaneously, the flow of intermediate and deep water masses would have rapidly decreased in the Southeast Alboran Sea. Warmer bottom waters and fluvial related eutrophication were detrimental for bryozoans but proved optimal for cold-water coral growth. Benthic foraminiferal assemblages suggest that the sediment/water interface was poorly oxygenated during the Early Greenland Interstadial 1.

During the mid-Greenland Interstadial 1, bottom currents strengthened, clearly separating this period between an early and a later phase. At the same time, the origin of the organic matter becomes more marine, whilst the benthic foraminiferal fauna is subject to compositional changes. We suggest that the impact of fluvial input on benthic communities was reduced during the later phase of Greenland Interstadial 1. More arid continental conditions together with a re-establishment of dynamic intermediate and deep water mass circulation would have characterized this period. However, our results stress that conditions were unstable at the seafloor during the Late Greenland Interstadial 1, with periods of increased fluvial input resulting in eutrophication and low oxygen availability. In contrast, bottom currents were weaker during the Holocene whilst water temperatures rapidly increased. This led to the development of renewed microfaunal and coral communities, *M. oculata* becoming the dominant species instead of *D. pertusum*.

## DATA AVAILABILITY STATEMENT

All the datasets used in this study are available at the open-access repository PANGEA using the following link: doi.pangaea.de/10.1594/PANGAEA.914702.

## AUTHOR CONTRIBUTIONS

RF acquired the data, interpreted the data, wrote and coordinated the manuscript. EF and TA participated in the data acquisition. AR participated in data interpretation and in writing the manuscript. TV and IH participated in data acquisition, interpretation and in writing the manuscript. DV head organizing member of the research cruise, participated in writing the manuscript. AF participated in data interpretation and in writing the manuscript, provided financial support.

## ACKNOWLEDGMENTS

We thank the SNF (Swiss National Science Foundation) funded projects FN-20020_153125 "Faunal assemblages from active, declining and buried cold-water coral ecosystems", FN-200021_175587 "Unconventional carbonate factories in
the Eastern Mediterranean: cold-water coral ecosystems and seeps" and FN-200021_149247 "4D-Diagenesis@Mound: Understanding the temporal and spatial variability of early diagenesis in carbonate mounds" for funding this research. Mr. Efraim Hall is thanked for his help with the CT quantification of the bioclasts and sedimentary core descriptions. We also thank Prof. Norbert Frank (Institute of Environmental Physics, University of Heidelberg, Germany) for performing and providing the Uranium-series isotopemeasurements. We are grateful for the ship time provided by IPEV on the R/V Marion Dufresne II within the framework of the EuroFLEETS GATEWAYS

project (grant agreement 228344). We further thank Dr. Tim Collart for the help provided with the Rysgran Package for R and Prof. Sébastien Bertrand (Department of Geology, Ghent University, Belgium) for the assistance with the grain-size measurements.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.00354/full#supplementary-material

## REFERENCES

Addamo, A. M., Vertino, A., Stolarski, J., García-Jiménez, R., Taviani, M., and Machordom, A. (2016). Merging scleractinian genera: the overwhelming genetic similarity between solitary Desmophyllum and colonial Lophelia. *BMC Evol. Biol.* 16:108. doi: 10.1186/s12862-016-0654-8

Alve, E. (1991). Foraminifera, climatic change, and pollution: a study of the late Holocene sediments in Drammensfjord, southeast Norway. *Holocene* 1, 243-261. doi: 10.1177/095968369100100306

Ausín, B., Flores, J. A., Sierro, F. J., Cacho, I., Hernández-Almeida, I., Martrat, B., et al. (2015). Atmospheric patterns driving Holocene productivity in the Alboran Sea (Western Mediterranean): a multiproxy approach. *Holocene* 25, 583-595. doi: 10.1177/0959683614565952

Barras, C., Jorissen, F. J., Labrune, C., Andral, B., and Boissery, P. (2014). Live benthic foraminiferal faunas from the French Mediterranean Coast: towards a new biotic index of environmental quality. *Ecol. Indic.* 36, 719-743. doi: 10.1016/j.ecolind.2013.09.028

Barth, H. (1857). *Travels and Discoveries in North and Central Africa*. New York, NY: Harper and Brothers.

Baudin, F., Disnar, J.-R., Aboussou, A., and Savignac, F. (2015). Guidelines for Rock-Eval analysis of recent marine sediments. *Organ. Geochem.* 86, 71-80. doi: 10.1016/j.orggeochem.2015.06.009

Behar, F., Beaumont, V., De, B., and Penteado, H. L. (2001). Rock-Eval 6 technology: performances and developments. *Oil Gas Sci. Technol.* 56, 111-134. doi: 10.2516/ogst:2001013

Bernhard, J. M., Buck, K. R., and Barry, J. P. (2001). Monterey Bay cold-seep biota: assemblages, abundance, and ultrastructure of living foraminifera. *Deep Sea Res.* 1 48, 2233-2249. doi: 10.1016/S0967-0637(01)00017-6

Bernhard, J. M., and Sen Gupta, B. K. (1999). "Foraminifera of oxygen-depleted environments," in *Modern Foraminifera*, ed. B. K. Sen Gupta (Philip Drive Norwell, MA: Kluwer Publishers), 201-216. doi: 10.1007/0-306-48104-9_12

Best, M. A., and Thorpe, J. P. (1994). "Particle size, clearance rate and feeding efficiency in marine Bryozoa," in *Biology and Palaeobiology of Bryozoans*, eds P. J. Hayward, J. S. Ryland, and P. D. Taylor (Fredensborg: Olsen and Olsen), 9-14.

Bjerager, M., and Surlyk, F. (2007). Benthic palaeoecology of Danian deep-shelf bryozoan mounds in the Danish Basin. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 250, 184-215. doi: 10.1016/j.palaeo.2007.03.008

Bond, G., Kromer, B., Beer, J., Muscheler, R., Evans, M. N., Showers, W., et al. (2001). Persistant solar influence on north atlantic climate during the holocene. *Science* 294, 2130-2136. doi: 10.1126/science.1065680

Bray, J. R., and Curtis, J. T. (1957). An ordination of the upland forest communities of Southern Wisconsin. *Ecol. Monogr.* 27, 325-349. doi: 10.2307/1942268

Cacho, I., Grimalt, J. O., Canals, M., Sbaffi, L., Shackleton, N. J., Schönfeld, J., et al. (2001). Variability of the western Mediterranean Sea surface temperature during the last 25,000 years and its connection with the Northern Hemisphere climatic changes. *Paleoceanography* 16, 40-52. doi: 10.1029/2000PA000502

Cacho, I., Shackleton, N., Elderfield, H., Sierro, F. J., and Grimalt, J. O. (2006). Glacial rapid variability in deep-water temperature and 818O from the Western Mediterranean Sea. *Q. Sci. Rev.* 25, 3294-3311. doi: 10.1016/j.quascirev.2006.10.004

Cacho, I., Sierro, F., Shackleton, N. J., Elderfield, H., and Grimalt, J. (2002). Reconstructing Alboran Sea hydrography during the last organic rich layer formation, *Geochim. Cosmochim. Acta* 66:A115.

Cimerman, F., and Langer, M. R. (1991). *Mediterranean Foraminifera*. Ljubljana: Historia Naturalis.

Clarke, K. R., and Gorley, R. N. (2006). PRIMER v6: User Manual/Tutorial (Plymouth Routines in Multivariate Ecological Research). Plymouth: PRIMER-E.

Comas, M., Pinheiro, L. M., Ivanov, M., and TTR-17 Leg 1 Scientific Party (2009). "Deep-water coral mounds in the Alboran Sea: the Melilla mound field revisited," in *Geo-Marine Research on the Mediterranean and European-Atlantic Margins. International Conference and TTR-17 Post-Cruise Meeting of the Training-through-Research Programme*, eds M. Comas and A. Y. Suzyumov (Granada, Spain: UNESCO), 8-9.

Comas, M., Zahn, R., Klaus, A., and Party, S. S. (1996). *Proceedings of the Ocean Drilling Program, Initial Reports*. College Station, TX: Ocean Drilling Program. doi: 10.2973/odp.proc.ir.161.1996

Comas, M. C., Platt, J. P., Soto, J. I., and Watts, A. B. (1999). "The origin and tectonic history of the Alboran Basin: insights from Leg 161 results," in *Proceedings of the Ocean Drilling Program 161 Scientific Results*, College Station, TX. doi: 10.2973/odp.proc.sr.161.262.1999

Corliss, B. H. (1985). Microhabitats of benthic foraminifera within deep-sea sediments. *Nature* 314, 435-438. doi: 10.1038/314435a0

Davies, A. J., Duineveld, G., Lavaleye, M., Bergman, M., van Haren, H., and Roberts, J. (2009). Downwelling and deep-water bottom currents as food supply mechanisms to the cold-water coral *Lophelia pertusa* (Scleractinia) at the Mingulay Reef Complex. *Limnol. Oceanogr.* 54, 620-629. doi: 10.4319/lo.2009.54.2.0620

Davies, A. J., and Guinotte, J. M. (2011). Global habitat suitability for framework-forming cold-water corals. *PLoS One* 6:e18483. doi: 10.1371/journal.pone.0018483

De Rijk, S., Jorissen, F., Rohling, E. J., and Troelstra, S. R. (2000). Organic flux control on bathymetric zonation of Mediterranean benthic foraminifera. *Mar. Micropaleontol.* 40, 151-166. doi: 10.1016/S0377-8398(00)00037-2

Delvaux, D., Martin, H., Leplat, P., and Paulet, J. (1989). Geochemical characterization of sedimentary organic matter by means of pyrolysis kinetic parameters. *Adv. Organ. Geochem.* 16, 175-187. doi: 10.1016/0146-6380(90)90038-2

deMenocal, P., Ortiz, J., Guilderson, T., Adkins, J., Sarnthein, M., Baker, L., et al. (2000). Abrupt onset and termination of the African Humid Period: rapid climate responses to gradual insolation forcing. *Q. Sci. Rev.* 19, 347-361. doi: 10.1016/S0277-3791(99)00081-5

deMenocal, P. B., and Tierney, J. E. (2012). Green sahara: african humid periods paced by Earth's orbital changes. *Nat. Educ. Knowl.* 3:12.

Do Couto, D., Gorini, C., Jolivet, L., Debret, N., Augier, R., Gumiaux, C., et al. (2016). Tectonic and stratigraphic evolution of the Western Alboran Sea Basin in the last 25Myrs. *Tectonophysics* 677-678, 280-311. doi: 10.1016/j.tecto.2016.03.020

Dorschel, B., Hebbeln, D., Rüggeberg, A., and Dullo, C. (2005). Carbonate budget of a cold-water coral carbonate mound: propeller mound, porcupine Seabight. *Int. J. Earth Sci.* 96, 73-83. doi: 10.1007/s00531-005-0493-0
Duggen, S., Hoernle, K., Klügel, A., Geldmacher, J., Thirlwall, M., Hauff, F., et al. (2008). Geochemical zonation of the Miocene Alborán Basin volcanism (westernmost Mediterranean): geodynamic implications. *Contrib. Mineral. Petrol.* 156, 577-593. doi: 10.1007/s00410-008-0302-4

Duineveld, G. C. A., Lavaleye, M. S. S., and Berghuis, E. M. (2004). Particle flux and food supply to a seamount cold-water coral community (Galicia Bank, NW Spain). *Mar. Ecol. Prog. Ser.* 277, 13-23. doi: 10.3354/meps277013

Dullo, W.-C., Flögel, S., and Rüggeberg, A. (2008). Cold-water coral growth in relation to the hydrography of the Celtic and Nordic European continental margin. *Mar. Ecol. Prog. Ser.* 371, 165-176. doi: 10.3354/meps07623

Eisele, M., Frank, N., Wienberg, C., Hebbeln, D., López Correa, M., Douville, E., et al. (2011). Productivity controlled cold-water coral growth periods during the last glacial off Mauritania. *Mar. Geol.* 280, 143-149. doi: 10.1016/j.margeo.2010.12.007

Eisele, M., Hebbeln, D., and Wienberg, C. (2008). Growth history of a cold-water coral covered carbonate mound — Galway Mound, Porcupine Seabight, NE-Atlantic. *Mar. Geol.* 253, 160-169. doi: 10.1016/j.margeo.2008.05.006

Emeis, K. C., Sakamoto, T., Wehausen, R., and Brumsack, H.-J. (2000). The sapropel record of the eastern Mediterranean Sea - results of Ocean Drilling Program Leg 160. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 158, 371-395. doi: 10.1016/S0031-0182(00)000539-0

Emelyanov, E. M., and Shimkus, K. M. (2012). *Geochemistry and Sedimentology of the Mediterranean Sea*. Cham: Springer Science and Business Media.

Emig, C. C. (1988). Les brachiopodes actuels, indicateurs (paléo) bathymétriques? Paléobathymétrie, éustasisme et séquences de dépôts, Marseille. *Géol. Médicerranéenne* 1, 65-71. doi: 10.3406/geolm.1988.1395

Emig, C. C., and Arnaud, P. M. (1988). Observations en submersible sur la densité des populations de *Gryphus vitreus* (Brachiopode) le long de la marge continentale de Provence (Méditerranée nord-occidentale). *C. R. Acad. Sci. Paris* 306, 501-505.

Ercilla, G., Juan, C., Hernández-Molina, J., Bruno, M., Estrada, F., Alonso, B., et al. (2016). Significance of bottom currents in deep-sea morphodynamics: an example from the Alboran Sea. *Mar. Geol.* 378, 157-170. doi: 10.1016/j.margeo.2015.09.007

Fabri, M. C., Pedel, L., Beuck, L., Galgani, F., Hebbeln, D., and Freiwald, A. (2014). Megafauna of vulnerable marine ecosystems in French mediterranean submarine canyons: spatial distribution and anthropogenic impacts. *Deep Sea Res. I* 104, 184-207. doi: 10.1016/j.dsr2.2013.06.016

Faccenna, C., Piromallo, C., Crespo-Blanc, A., Jolivet, L., and Rossetti, F. (2004). Lateral slab deformation and the origin of the western Mediterranean arcs. *Tectonics* 23:TC1012. doi: 10.1029/2002TC001488

Feenstra, E. J. (2020). *Carbonate Preservation Pathways in Cold-Water Coral Mounds*. PhD thesis, University of Fribourg, Switzerland.

Fentimen, R., Lim, A., Rüggeberg, A., Wheeler, A. J., Van Rooij, D., and Foubert, A. (2020). Impact of bottom water currents on benthic foraminiferal assemblages in a cold-water coral environment: the Moira Mounds (NE Atlantic). *Mar. Micropaleontol.* 154:101799. doi: 10.1016/j.marmicro.2019.10.1799

Fentimen, R., Rüggeberg, A., Lim, A., El Kateb, A., Foubert, A., Wheeler, A. J., et al. (2018). Benthic foraminifera in a deep-sea high-energy environment: the Moira Mounds (Porcupine Seabight, SW of Ireland). *Swiss J. Geosci.* 111, 561-572. doi: 10.1007/s00015-018-0317-4

Field, J. G., Clarke, K. R., and Warwick, R. M. (1982). A practical strategy for analysing multispecies distribution patterns. *Mar. Ecol. Prog. Ser.* 8, 37-52. doi: 10.3354/meps008037

Fink, H. G., Wienberg, C., De Pol-Holz, R., and Hebbeln, D. (2015). Spatiotemporal distribution patterns of Mediterranean cold-water corals (*Lophelia pertusa* and *Madrepora oculata*) during the past 14,000 years. *Deep Sea Res. I* 103, 37-48. doi: 10.1016/j.dsr.2015.05.006

Fink, H. G., Wienberg, C., De Pol-Holz, R., Wintersteller, P., and Hebbeln, D. (2013). Cold-water coral growth in the Alboran Sea related to high productivity during the Late Pleistocene and Holocene. *Mar. Geol.* 339, 71-82. doi: 10.1016/j.margeo.2013.04.009

Fontanier, C., Jorissen, F. J., Lansard, B., Mouret, A., Buscail, R., Schmidt, S., et al. (2008). Live foraminifera from the open slope between Grand Rhône and Petit Rhône Canyons (Gulf of Lions, NW Mediterranean). *Deep Sea Res. I* 55, 1532-1553. doi: 10.1016/j.dsr.2008.07.003

Fontanier, C., Jorissen, F. J., Licari, L., Alexandre, A., Anschutz, P., and Carbonel, P. (2002). Live benthic foraminiferal faunas from the Bay of Biscay: faunal density, composition, and microhabitats. *Deep Sea Res. I* 49, 751-785. doi: 10.1016/S0967-0637(01)00078-4

Fontanier, C., Mackensen, A., Jorissen, F. J., Anschutz, P., Licari, L., and Griveaud, C. (2006). Stable oxygen and carbon isotopes of live benthic foraminifera from the Bay of Biscay: microhabitat impact and seasonal variability. *Mar. Micropaleontol.* 58, 159-183. doi: 10.1016/j.marmicro.2005.09.004

Frank, N., Paterne, M., Ayliffe, L., van Weering, T., Henriet, J.-P., and Blamart, D. (2004). Eastern North Atlantic deep-sea corals: tracing upper intermediate water $\Delta^{14}C$ during the Holocene. *Earth Planet. Sci. Lett.* 219, 297-309. doi: 10.1016/S0012-821X(03)00721-0

Frank, N., Ricard, E., Lutringer-Paquet, A., van der Land, C., Colin, C., Blamart, D., et al. (2009). The Holocene occurrence of cold water corals in the NE Atlantic: implications for coral carbonate mound evolution. *Mar. Geol.* 266, 129-142. doi: 10.1016/j.margeo.2009.08.007

Freiwald, A. (2002). "Reef-forming cold-water corals," in *Ocean Margin Systems*, eds G. Wefer, D. Billet, D. Hebbeln, B. B. Jørgensen, M. Schlüter, and T. C. E. van Weering (Berlin: Springer), 365-385. doi: 10.1007/978-3-662-05127-6_23

Freiwald, A., Beuck, L., Rüggeberg, A., Taviani, M., and Hebbeln, D. (2009). The white coral community in the centra Mediterranean Sea revealed by ROV surveys. *Oceanography* 22, 58-74. doi: 10.5670/oceanog.2009.06

Freiwald, A., Fosså, J. H., Grehan, A., Koslow, T., and Roberts, J. M. (2004). *Cold-Water Coral Reefs*. Cambridge, MA: UNEP_WCMC.

Frigola, J., Moreno, A., Cacho, I., Canals, M., Sierro, F. J., Flores, J. A., et al. (2008). Evidence of abrupt changes in Western Mediterranean Deep Water circulation during the last 50kyr: a high-resolution marine record from the Balearic Sea. *Q. Int.* 181, 88-104. doi: 10.1016/j.quaint.2007.06.016

Gasse, F. (2000). Hydrological changes in the African tropics since the Last Glacial Maximum. *Q. Sci. Rev.* 19, 189-211. doi: 10.1016/S0277-3791(99)00061-X

Gilbert, E. R., Camargo, M. G., and Sandrini-Neto, L. (2015). *Rysgran: Grain Size Analysis, Textural Classifications and Distribution of Unconsolidated Sediments, R Package Version 2.1.0*. Available online at: https://CRAN.R-project.org/package=rysgran

Gooday, A. J. (1993). The biology of deep-sea foraminifera: a review of some advances and their applications in paleoceanography. *Palaios* 9, 14-31. doi: 10.2307/3515075

Gross, O. (2000). Influence of temperature, oxygen and food availability on the migrational activity of bathyal benthic foraminifera: evidence by microcosm experiments. *Hydrobiologia* 426, 123-137. doi: 10.1023/A:10039308 31220

Hanz, U., Wienberg, C., Hebbeln, D., Duineveld, G., Lavaleye, M., Juva, K., et al. (2019). Environmental factors influencing cold-water coral ecosystems in the oxygen minimum zones on the Angolan and *Namibian margins*. *Biogeosci. Discuss.* 16, 4337-4356. doi: 10.5194/bg-16-4337-2019

Hebbeln, D., Bender, M., Gaide, S., Titschack, J., Vandorpe, T., Van Rooij, D., et al. (2019). Thousands of cold-water coral mounds along the Moroccan Atlantic continental margin: distribution and morphometry. *Mar. Geol.* 411, 51-61. doi: 10.1016/j.margeo.2019.02.001

Hebbeln, D., Van Rooij, D., and Wienberg, C. (2016). Good neighbours shaped by vigorous currents: cold-water coral mounds and contourites in the North Atlantic. *Mar. Geol.* 378, 171-185. doi: 10.1016/j.margeo.2016.01.014

Hebbeln, D., Wienberg, C., Beuck, L., Freiwald, A., and Wintersteller, P. (2009). *Report and Preliminary Results of RV POSEIDON Cruise POS 385 "Cold-water Corals of the Alboran Sea (Western Mediterranean Sea)", Faro-Toulon, May 29-June 16, 2009*. Bremen: University of Bremen.

Heburn, G. W., and La Violette, P. E. (1990). Variations in the structure of the anticyclonic gyres found in the Alboran Sea. *J. Geophys. Res.* 95, 1599-1613. doi: 10.1029/JC095iC02p01599

Heinz, P., Kitazato, H., Schmiedl, G., and Hemleben, C. (2001). Response of deep-sea benthic foraminifera from the Mediterranean Sea to simulated phytoplankton pulses under laboratory conditions. *J. Foraminiferal Res.* 31, 210-227. doi: 10.2113/31.3.210

Hernández-Molina, F. J., Somoza, L., Vazquez, J. T., Lobo, F., Fernández-Puga, M. C., Llave, E., et al. (2002). Quaternary stratigraphic stacking patterns on the continental shelve of the southern Iberian Peninsula: their relationship with global climate and palaeoceanographic changes. *Q. Int.* 92, 5-23. doi: 10.1016/S1040-6182(01)00111-2
Holbourn, A., Kuhnt, W., and James, N. (2002). Late Pleistocene bryozoan reef mounds of the Great Australian Bight: isotope stratigraphy and benthic foraminiferal record. *Paleoceanography* 17, 11–14. doi: 10.1029/2001PA000643

Husum, K., and Hald, M. (2004). Modern foraminifera distribution in the subarctic Malangen fjord and adjoining shelf, northern Norway. *J. Foraminiferal Res.* 34, 34–48. doi: 10.2113/0340034

Huvenne, V. A. I., Masson, D. G., and Wheeler, A. J. (2009). Sediment dynamics of a sandy contourite: the sedimentary context of the Darwin cold-water coral mounds, Northern Rockall Trough. *Int. J. Earth Sci.* 98, 865–884. doi: 10.1007/s00531-008-0312-5

Incarbona, A., Sprovieri, M., Lirer, F., and Sprovieri, R. (2011). Surface and deep water conditions in the Sicily channel (central Mediterranean) at the time of sapropel S5 deposition. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 306, 243–248. doi: 10.1016/j.palaeo.2011.04.030

James, N. P., Feary, D. A., Surlyk, F., Toni Simo, J. A., Betzler, C., Holbourn, A. E., et al. (2000). Quaternary bryozoan reef mounds in cool-water, upper slope environments: great Australian Bight. *Geology* 28, 647–650. doi: 10.1130/0091-7613(2000)28<647:QBRMIC>2.0.CO;2

Jimenez-Espejo, F. J., Martinez-Ruiz, F., Rogerson, M., González-Donoso, J. M., Romero, O. E., Linares, D., et al. (2008). Detrital input, productivity fluctuations, and water mass circulation in the westernmost Mediterranean Sea since the Last Glacial Maximum. *Geochem. Geophys. Geosyst.* 9, 1–19. doi: 10.1029/2008GC002096

Jolly, D., Prentice, I. C., Bonnefille, R., Ballouche, A., Bengo, M., Brenac, P., et al. (1998). Biome reconstruction from pollen and plant macrofossil data for Africa and the Arabian peninsula at 0 and 6000 years. *J. Biogeogr.* 25, 1007–1027. doi: 10.1046/j.1365-2699.1998.00238.x

Jorissen, F. J. (1987). The distribution of benthic foraminifera in the Adriatic Sea. *Mar. Micropaleontol.* 12, 21–48. doi: 10.1016/0377-8398(87)90012-0

Katz, B. J. (1983). Limitations of “Rock-Eval” pyrolysis for typing organic matter. *Organ. Geochem.* 4, 195–199. doi: 10.1016/0146-6380(83)90041-4

Katz, E. J. (1972). The levantine intermediate water between the strait of sicily and the strait of gibraltar. *Deep Sea Res.* 19, 507–520. doi: 10.1016/0011-7471(72)90018-6

Kuper, R., and Kröpelin, S. (2006). Climate-controlled holocene occupation in the sahara: motor of Africa's evolution. *Science* 313, 803–807. doi: 10.1126/science.1130989

La Violette, P. E. (1983). The advection of submesoscale thermal features in the Alboran Sea Gyres. *J. Phys. Oceanogr.* 14, 550–565. doi: 10.1175/1520-0485(1984)014<0550:TAOSTF>2.0.CO;2

Lafuente, J. G., Camno, N., Vargas, M., Rubin, J. P., and Hernández-Guerra, A. (1998). Evolution of the Alboran Sea hydrographic structures during July 1993. *Deep Sea Res. I* 45, 39–65. doi: 10.1016/S0967-0637(97)00216-1

Lanoix, R. (1974). *Projet Alboran: Etude Hydrologique et Dynamique de la Mer d'Alboran*. Technical Report 66. Brussels: NATO.

L'Helguen, S., Le Corre, P., Madec, C., and Morin, P. (2002). New and regenerated production in the Almeria-Oran front area, eastern Alboran Sea. *Deep Sea Res. I* 49, 83–99. doi: 10.1016/S0967-0637(01)00044-9

Linke, P., and Lutze, G. F. (1993). Microhabitat preferences of benthic foraminifera—a static concept or a dynamic adaptation to optimize food acquisition? *Mar. Micropaleontol.* 20, 215–234. doi: 10.1016/0377-8398(93)90034-U

Lo Iacono, C., Gràcia, E., Ranero, C. R., Emelianov, M., Huvenne, V. A. I., Bartolomé, R., et al. (2014). The West Melilla cold water coral mounds, Eastern Alboran Sea: morphological characterization and environmental context. *Deep Sea Res. II* 99, 316–326. doi: 10.1016/j.dsr2.2013.07.006

López Casado, C., Sanz de Galdeano, C., Molina Palacios, S., and Henares Romero, J. (2001). The structure of the Alboran Sea: an interpretation from seismological and geological data. *Tectonophysics* 338, 79–95. doi: 10.1016/S0040-1951(01)00059-2

Lutze, G. F., and Coulbourn, W. T. (1984). Recent benthic foraminifera from the continental margin of Northwest Africa: community structure and distribution. *Mar. Micropaleontol.* 8, 361–401. doi: 10.1016/0377-8398(84)90002-1

Mackensen, A., Schmiedl, G., Harloff, J., and Giese, M. (1995). Deep-sea foraminifera in the South Atlantic Ocean: ecology and assemblage generation. *Micropaleontology* 41, 342–358. doi: 10.2307/1485808

Margreth, S., Gennari, G., Rüggeberg, A., Comas, M. C., Pinheiro, L. M., and Spezzaferri, S. (2011). Growth and demise of cold-water coral ecosystems on

mud volcanoes in the West Alboran Sea: the messages from the planktonic and benthic foraminifera. *Mar. Geol.* 282, 26–39. doi: 10.1016/j.margeo.2011.02.006

Margreth, S., Rüggeberg, A., and Spezzaferri, S. (2009). Benthic foraminifera as bioindicator for cold-water coral reef ecosystems along the Irish margin. *Deep Sea Res. I* 56, 2216–2234. doi: 10.1016/j.dsr.2009.07.009

Martinez-Ruiz, F., Kastner, M., Gallego-Torres, D., Rodrigo-Gámiz, M., Nieto-Moreno, V., and Ortega-Huertas, M. (2015). Paleoclimate and paleoceanography over the past 20,000 yr in the Mediterranean Sea Basins as indicated by sediment elemental proxies. *Q. Sci. Rev.* 107, 25–46. doi: 10.1016/j.quascirev.2014.09.018

Martínez-Ruiz, F., Paytan, A., Kastner, M., González-Donoso, J. M., Linares, D., Bernasconi, S. M., et al. (2003). A comparative study of the geochemical and mineralogical characteristics of the S1 sapropel in the western and eastern Mediterranean. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 190, 23–37. doi: 10.1016/S0031-0182(02)00597-7

McCave, I. N., and Hall, I. R. (2006). Size sorting in marine muds: processes, pitfalls, and prospects for paleoflow-speed proxies. *Geochem. Geophys. Geosyst.* 7, 1–37. doi: 10.1029/2006GC001284

McCave, I. N., Manighetti, B., and Robinson, S. G. (1995). Sortable silt and fine sediment size/composition slicing: parameters for palaeocurrent speed and palaeoceanography. *Paleoceanography* 10, 593–610. doi: 10.1029/94PA03039

McGann, M., and Conrad, J. E. (2018). Faunal and stable isotopic analyses of benthic foraminifera from the Southeast Seep on Kimki Ridge offshore southern California, USA. *Deep Sea Res. II* 150, 92–117. doi: 10.1016/j.dsr2.2018.01.011

Mienis, F., de Stigter, H. C., de Haas, H., and van Weering, T. C. E. (2009). Near-bed particle deposition and resuspension in a cold-water coral mound area at the Southwest Rockall Trough margin, NE Atlantic. *Deep Sea Res. I* 56, 1026–1038. doi: 10.1016/j.dsr.2009.01.006

Mienis, F., de Stigter, H. C., White, M., Duineveld, G., de Haas, H., and van Weering, T. C. E. (2007). Hydrodynamic controls on cold-water coral growth and carbonate-mound development at the SW and SE Rockall Trough Margin, NE Atlantic Ocean. *Deep Sea Res. I* 54, 1655–1674. doi: 10.1016/j.dsr.2007.05.013

Milker, Y., Schmiedl, G., Betzler, C., Römer, M., Jaramillo-Vogel, D., and Siccha, M. (2009). Distribution of recent benthic foraminifera in shelf carbonate environments of the Western Mediterranean Sea. *Mar. Micropaleontol.* 73, 207–225. doi: 10.1016/j.marmicro.2009.10.003

Millot, C. (1999). Circulation in the western Mediterranean Sea. *J. Mar. Syst.* 20, 423–442. doi: 10.1016/S0924-7963(98)00078-5

Millot, C. (2009). Another description of the Mediterranean Sea outflow. *Prog. Oceanogr.* 82, 101–124. doi: 10.1016/j.poccean.2009.04.016

Millot, C. (2013). Levantine intermediate water characteristics: an astounding general misunderstanding! *Sci. Mar.* 77, 217–232. doi: 10.3989/scimar.03518.13A

Millot, C., Candela, J., Fuda, J.-L., and Tber, Y. (2006). Large warming and salinification of the Mediterranean outflow due to changes in its composition. *Deep Sea Res. I* 53, 656–666. doi: 10.1016/j.dsr.2005.12.017

Mojtahid, M., Jorissen, F., Lansard, B., Fontanier, C., Bombled, B., and Rabouille, C. (2009). Spatial distribution of live benthic foraminifera in the Rhône prodelta: faunal response to a continental–marine organic matter gradient. *Mar. Micropaleontol.* 70, 177–200. doi: 10.1016/j.marmicro.2008.12.006

Murray, J. W. (2006). *Ecology and Applications of Benthic Foraminifera*. Cambridge, MA: Cambridge University Press. doi: 10.1017/CBO9780511535529

Olivet, J. L., Auzende, J. M., and Bonnin, J. (1973). Structure et évolution tectonique du bassin d'Alboran. *Bull. Soc. Géol. France* 7, 491–495. doi: 10.2113/gsgfbull.S7-XV.2.108

Orejas, C., Gori, A., and Gili, J. M. (2007). Growth rates of live Lophelia pertusa and Madrepora oculata from the Mediterranean Sea maintained in aquaria. *Coral Reefs* 27, 255–255. doi: 10.1007/s00338-007-0350-7

Phleger, F. B., and Soutar, A. (1973). Production of benthic foraminifera in three east pacific oxygen minima. *Micropaleontology* 19, 110–115. doi: 10.2307/1484973

R Core Team (2018). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing.

Raddatz, J., Rüggeberg, A., Margreth, S., and Dullo, W.-C. (2011). Paleoenvironmental reconstruction of Challenger Mound initiation in the Porcupine Seabight, NE Atlantic. *Mar. Geol.* 282, 79–90. doi: 10.1016/j.margeo.2010.10.019
Ramsey, C. (2017). *OxCal 4.2.4. Electronic Document*. Available online at: https://c14.arch.ox.ac.uk/oxcal.html (accessed February 4, 2020).

Rasmussen, S. O., Bigler, M., Blockley, S. P., Blunier, T., Buchardt, S. L., Clausen, H. B., et al. (2014). A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy. *Quat. Sci. Rev.* 106, 14–28. doi: 10.1016/j.quascirev.2014.09.007

Rathburn, A. E., and Corliss, B. H. (1994). The ecology of living (stained) deep-sea from the Sulu Sea. *Paleoceanography* 9, 87–150. doi: 10.1029/93PA02327

Reimer, P. J., Bard, E., Bayliss, A., Beck, J. W., Blackwell, P. G., Ramsey, C. B., et al. (2013). IntCal13 and marine13 radiocarbon age calibration curves 0–50,000 years cal BP. *Radiocarbon* 55, 1869–1887. doi: 10.2458/azu_js_rc.55.16947

Roberts, J. M., Wheeler, A. J., Freiwald, A., and Cairns, S. (2009). *Cold-Water Corals*. Cambridge, MA: Cambridge University Press. doi: 10.1017/CBO9780511581588

Rodrigo-Gámiz, M., Martínez-Ruiz, F., Jiménez-Espejo, F. J., Gallego-Torres, D., Nieto-Moreno, V., Romero, O., et al. (2011). Impact of climate variability in the western Mediterranean during the last 20,000 years: oceanic and atmospheric responses. *Q. Sci. Rev.* 30, 2018–2034. doi: 10.1016/j.quascirev.2011.05.011

Rogerson, M., Cacho, I., Jimenez-Espejo, F., Reguera, M. I., Sierro, F. J., Martinez-Ruiz, F., et al. (2008). A dynamic explanation for the origin of the western Mediterranean organic-rich layers. *Geochem. Geophys. Geosyst.* 9, 1–20. doi: 10.1029/2007GC001936

Rogerson, M., Dublyansky, Y., Hoffmann, D. L., Luetscher, M., Töchterle, P., and Spötl, C. (2019). Enhanced Mediterranean water cycle explains increased humidity during MIS 3 in North Africa. *Clim. Past* 15, 1757–1769. doi: 10.5194/cp-15-1757-2019

Rüggeberg, A., Dullo, C., Dorschel, B., and Hebbeln, D. (2007). Environmental changes and growth history of a cold-water carbonate mound (Propeller Mound, Porcupine Seabight). *Int. J. Earth Sci.* 96, 57–72. doi: 10.1007/s00531-005-0504-1

Savini, A., Vertino, A., Marchese, F., Beuck, L., and Freiwald, A. (2014). Mapping cold-water coral habitats at different scales within the Northern Ionian Sea (Central Mediterranean): an assessment of coral coverage and associated vulnerability. *PLoS One* 9:e87108. doi: 10.1371/journal.pone.0087108

Schiebel, R., and Hemleben, C. (2017). *Planktic Foraminifers in the Modern Ocean*. Berlin: Springer-Verlag. doi: 10.1007/978-3-662-50297-6

Schmiedl, G., De Bovée, F., Buscail, R., Charriere, B., Hemleben, C., Medernach, L., et al. (2000). Trophic control of benthic foraminiferal abundance and microhabitat in the bathyal Gulf of Lions, western Mediterranean Sea. *Mar. Micropaleontol.* 40, 167–188. doi: 10.1016/S0377-8398(00)00038-4

Scholz, J., and Hillmer, G. (1995). Reef-bryozoans and bryozoan-microreefs: control factor evidence from the Philippines and other regions. *Facies* 32, 109–144. doi: 10.1007/BF02536866

Schönfeld, J. (1997). The impact of the Mediterranean Outflow Water (MOW) on benthic foraminiferal assemblages and surface sediments at the southern Portuguese continental margin. *Mar. Micropaleontol.* 29, 211–236. doi: 10.1016/S0377-8398(96)00050-3

Schönfeld, J. (2002). Recent benthic foraminiferal assemblages in deep high-energy environments from the Gulf of Cadiz (Spain). *Mar. Micropaleontol.* 44, 141–162. doi: 10.1016/S0377-8398(01)00039-1

Schönfeld, J., Dullo, W.-C., Pfannkuche, O., Freiwald, A., Rüggeberg, A., Schmidt, S., et al. (2011). Recent benthic foraminiferal assemblages from cold-water coral mounds in the Porcupine Seabight. *Facies* 57, 187–213. doi: 10.1007/s10347-010-0234-0

Sen Gupta, B. K., and Machain-Castillo, M. L. (1993). Benthic foraminifera in oxygen-poor habitats. *Mar. Micropaleontol.* 20, 183–201. doi: 10.1016/0377-8398(93)90032-S

Siani, G., Paterne, M., Arnold, M., Bard, E., Métivier, B., Tisnerat, N., et al. (2000). Radiocarbon reservoir ages in the Mediterranean Sea and Black Sea. *Radiocarbon* 42, 272–280. doi: 10.1017/S0033822200059075

Sierro, F. J., Hodell, D. A., Curtis, J. H., Flores, J. A., Reguera, I., Colmenero-Hidalgo, E., et al. (2005). Impact of iceberg melting on Mediterranean thermohaline circulation during Heinrich events. *Paleoceanography* 20:A2019. doi: 10.1029/2004PA001051

Skonieczny, C., Paillou, P., Bory, A., Bayon, G., Biscara, L., Crosta, X., et al. (2015). African humid periods triggered the reactivation of a large river system in Western Sahara. *Nat. Commun.* 6:8751. doi: 10.1038/ncomms9751

Sparnocchia, S., Picco, P., Manzella, G., Ribotti, A., Copello, S., and Brasey, P. (1995). Intermediate water formation in the Ligurian Sea. *Oceanol. Acta* 18, 151–162.

Spezzaferri, S., Rüggeberg, A., Stalder, C., and Margreth, S. (2013). Benthic foraminifer assemblages from norwegian cold-water coral reefs. *J. Foraminiferal Res.* 43, 21–39. doi: 10.2113/gsjfr.43.1.21

Spötl, C., and Vennemann, T. W. (2003). Continuous-flow IRMS analysis of carbonate minerals. *Rapid Commun. Mass Spect.* 17, 1004–1006. doi: 10.1002/ rcm.1010

Stalder, C., El Kateb, A., Vertino, A., Rüggeberg, A., Camozzi, O., Pirkenseer, C. M., et al. (2018). Large-scale paleoceanographic variations in the western Mediterranean Sea during the last 34,000 years: from enhanced cold-water coral growth to declining mounds. *Mar. Micropaleontol.* 143, 46–62. doi: 10.1016/j. marmicro.2018.07.007

Stalder, C., Spezzaferri, S., Rüggeberg, A., Pirkenseer, C., and Gennari, G. (2014). Late Weichselian deglaciation and early Holocene development of a cold-water coral reef along the Lopphavet shelf (Northern Norway) recorded by benthic foraminifera and ostracoda. *Deep Sea Res. II* 99, 249–269. doi: 10.1016/j.dsr2. 2013.07.009

Stalder, C., Vertino, A., Rosso, A., Ruggeberg, A., Pirkenseer, C., Spangenberg, J. E., et al. (2015). Microfossils, a key to unravel cold-water carbonate mound evolution through time: evidence from the eastern Alboran Sea. *PLoS One* 10:e0140223. doi: 10.1371/journal.pone.0140223

Suhr, S. B., Pond, D. W., Gooday, A. J., and Smith, C. R. (2003). Selective feeding by benthic foraminifera on phytodetritus on the western Antarctic Peninsula shelf: evidence from fatty acid biomarker analysis. *Mar. Ecol. Prog. Ser.* 262, 153–162. doi: 10.3354/meps262153

Sun, X., Corliss, B. H., Brown, C. W., and Showers, W. J. (2006). The effect of primary productivity and seasonality on the distribution of deep-sea benthic foraminifera in the North Atlantic. *Deep Sea Res.* *I* 53, 28–47. doi: 10.1016/j. dsr.2005.07.003

Synal, H. A., Stocker, M., and Suter, M. (2007). MICADAS: a new compact radiocarbon AMS system. *Nucl. Instr. Methods Phys. Res.* *B* 259, 7–13. doi: 10.1016/j.nimb.2007.01.138

Szymanśka, N., Pawłowska, J., Kucharska, M., Kujawa, A., Łącka, M., and Zajączkowski, M. (2017). Impact of shelf-transformed waters (STW) on foraminiferal assemblages in the outwash and glacial fjords of Adventfjorden and Hornsund, Svalbard. *Oceanologia* 59, 525–540. doi: 10.1016/j.oceano.2017. 04.006

Talbot, M. R., and Livingstone, D. A. (1989). Hydrogen index and carbon isotopes of lacustrine organic matter as lake level indicators. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* 70, 121–137. doi: 10.1016/0031-0182(89)90084-9

Taviani, M., Angeletti, L., Canese, S., Cannas, R., Cardone, F., Cau, A., et al. (2017). The "Sardinian cold-water coral province" in the context of the Mediterranean coral ecosystems. *Deep Sea Res. II* 145, 61–78. doi: 10.1016/j.dsr2.2015.12.008

Taviani, M., Angeletti, L., Cardone, F., Montagna, P., and Danovaro, R. (2019). A unique and threatened deep water coral-bivalve biotope new to the Mediterranean Sea offshore the Naples megalopolis. *Sci. Rep.* 9:3411. doi: 10.1038/s41598-019-39655-8

Taviani, M., Freiwald, A., and Zibrowius, H. (2005). Deep coral growth in the Mediterranean Sea: an overview. *Cold Water Corals Ecosyst.* 3, 137–156. doi: 10.1007/3-540-27673-4_7

Thierens, M., Browning, E., Pirlet, H., Loutre, M. F., Dorschel, B., Huvenne, V. A. I., et al. (2013). Cold-water coral carbonate mounds as unique palaeo-archives: the Plio-Pleistocene Challenger Mound record (NE Atlantic). *Q. Sci. Rev.* 73, 14–30. doi: 10.1016/j.quascirev.2013.05.006

Tintore, J., Violette, P. E., Blade, I., and Cruzado, A. (1988). A study of an intense front in the eastern alboran sea: the almeira-oran front. *J. Phys. Oceanogr.* 18, 1384–1397. doi: 10.1175/1520-0485(1988)018<1384:ASOAID>2.0.CO;2

Toucanne, S., Jouet, G., Ducassou, E., Bassetti, M.-A., Dennielou, B., Angue Minto'o, C. M., et al. (2012). A 130,000-year record of Levantine Intermediate Water flow variability in the Corsica Trough, western Mediterranean Sea. *Q. Sci. Rev.* 33, 55–73. doi: 10.1016/j.quascirev.2011.11.020

Tyson, R. V. (1995). *Sedimentary Organic Matter - Organic Facies and Palynofacies*. London: Chapman and Hall. doi: 10.1007/978-94-011-0739-6
Vafidis, D., Koukoras, A., and Voultsiadou-Koukoura, E. (1997). Actinaria, corallimorpharia, and scleractinia (Hexacorallia, Anthozoa) of the Aegean Sea, with a checklist of the eastern Mediterranean and Black Sea species. *Isr. J. Zool.* 43, 55–70.

van Krevelen, D. W. (1993). *Coal: Typology-Physics-Chemistry-Constitution*, 3rd Edn. Amsterdam: Elsevier Science Publishers.

van Oevelen, D., Mueller, C. E., Lundålv, T., and Middelburg, J. J. (2016). Food selectivity and processing by the cold-water coral Lophelia pertusa. *Biogeosciences* 13, 5789–5798. doi: 10.5194/bg-13-5789-2016

van Rooij, D., Hebbeln, D., Comas, M., Vandorpe, T., and Delivet, S. (2013). *EUROFLEETS Cruise Summary Report. The Mediterranean-Atlantic Gateway Code: The Late Pleistocene Carbonate Mound Record*, 214.

Vertino, A., Stolarski, J., Taviani, M., and Bosellini, F. (2014). “Mediterranean corals: from Miocene to recent,” in *The Mediterranean Sea: Its History and Present Challenges*, eds S. Goffredo and Z. Dubinsky (Cham: Springer), 257–274. doi: 10.1007/978-94-007-6704-1_14

Viúdez, Á, and Tintoré, J. (1995). Time and space variability in the eastern Alboran Sea from March to May 1990. *J. Geophys. Res.* 100, 8571–8586. doi: 10.1029/94JC03129

Wang, H., Lo Iacono, C., Wienberg, C., Mitschack, J., and Hebbeln, D. (2019). Cold-water coral mounds in the southern Alboran Sea (western Mediterranean Sea): internal waves as an important driver for mound formation since the last deglaciation. *Mar. Geol.* 412, 1–18. doi: 10.1016/j.margeo.2019.02.007

Wefing, A.-M., Arps, J., Blaser, P., Wienberg, C., Hebbeln, D., and Frank, N. (2017). High precision U-series dating of scleractinian cold-water corals using an automated chromatographic U and Th extraction. *Chem. Geol.* 475, 140–148. doi: 10.1016/j.chemgeo.2017.10.036

Wheeler, A. J., Beyer, A., Freiwald, A., de Haas, H., Huvenne, V. A. I., Kozachenko, M., et al. (2007). Morphology and environment of cold-water coral carbonate mounds on the NW European margin. *Int. J. Earth Sci.* 96, 37–56. doi: 10.1007/s00531-006-0130-6

White, M., Mohn, C., De Stigter, H. C., and Mottram, G. (2005). “Deep-water coral development as a function of hydrodynamics and surface productivity around the submarine banks of the Rockall Trough, NE Atlantic,” in *Cold-water Corals and Ecosystems*, eds A. Freiwald and J. M. Roberts (Berlin: Springer-Verlag), 503–514. doi: 10.1007/3-540-27673-4_25

Wickham, H. (2016). *ggplot2: Elegant Graphics for Data Analysis*. New York, NY: Springer-Verlag. doi: 10.1007/978-3-319-24277-4

Wienberg, C., and Mitschack, J. (2015). Framework-forming scleractinian cold-water corals through space and time: a late quaternary north atlantic perspective, in *Marine Animal Forests*, eds S. Rossi, L. Bramanti, A. Gori, and C. Orejas Saco del Valle (Cham: Springer), 1–34. doi: 10.1007/978-3-319-17001-5_16-1

Wienberg, C., Mitschack, J., Freiwald, A., Frank, N., Lundålv, T., Taviani, M., et al. (2018). The giant Mauritanian cold-water coral mound province: oxygen control on coral mound formation. *Q. Sci. Rev.* 185, 135–152. doi: 10.1016/j.quascirev.2018.02.012

Wilson, J. B. (1979). 'Patch' development of the deep-water coral *Lophelia pertusa* (L.) on Rockall Bank. *J. Mar. Biol. Assoc. U. K.* 59, 165–177. doi: 10.1017/S0025315400046257

Zibrowius, H. (2003). *The "White Coral Community", Canyon and Seamount Faunas of the Deep Mediterranean Sea*. London: RAC/SPA.

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Fentimen, Feenstra, Rüggeberg, Vennemann, Hajdas, Adatte, Van Rooij and Foubert. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.