A 300 000-year record of cold-water coral mound build-up at the East Melilla Coral Province (SE Alboran Sea, western Mediterranean)

. This study provides a detailed reconstruction of cold-water coral mound build-up within the East Melilla Coral Province (southeastern Alboran Sea), more precisely at the northern part of Brittlestar Ridge I, over the last 300 kyr. The multiproxy investigation of core MD13-3462G reveals that mound build-up took place during both interglacial and glacial periods at average aggradation rates


Introduction
Cold-water coral (CWC) reefs are diverse marine ecosystems that are widespread in the world's oceans (Freiwald et al., 2004;Roberts et al., 2009). The most important reefbuilding CWC species in the Atlantic Ocean and Mediterranean Sea are the scleractinian species Desmophyllum pertusum (formerly known as Lophelia pertusa, see Addamo et al., 2016) and Madrepora oculata ). These predominantly suspension-feeding organisms depend on enhanced hydrodynamic regimes that provide food to their polyps (White et al., 2005;Mienis et al., 2007;Carlier et al., 2009;Davies et al., 2009;Roberts et al., 2009;Hanz et al., 2019). The role played by internal waves (i.e. waves that occur at the interface between two water masses of different densities) in the proliferation of CWCs is vital, since these oscillations increase turbulence, and hence nutrient supply, and accumulate particulate organic matter due to their sharp density gradient (White et al., 2005;Davies et al., 2009;Pomar et al., 2012;Wang et al., 2019). Physicochemical properties of the ambient water (e.g. salinity, temperature, dissolved oxygen concentrations, pH, density) also affect CWC growth (Freiwald et al., 2004;Dullo et al., 2008;Davies and Guinotte, 2011;Hanz et al., 2019). If favourable conditions are maintained over longer periods, successive reef generations may build CWC mounds through the interaction between coral growth and sediment accumulation (Wilson, 1979;Roberts et al., 2006;Foubert and Henriet, 2009;Roberts et al., 2009;Hebbeln et al., 2016). Consequently, CWC mounds can reach considerable heights of over 300 m and spread for kilometres in width and length at their base (De Mol et al., 2002;Kenyon et al., 2003;Huvenne et al., 2005). Mound development may span from thousands to millions of years and attain important mound aggradation rates, e.g. ± 400 cm kyr −1 in the East Melilla Coral Province (EMCP; Frank et al., 2009;López Correa et al., 2012;Fink et al., 2013;Stalder et al., 2015;Wienberg et al., 2018). As such, and in spite of mound formation being generally discontinuous, CWC mounds are valuable environmental and climate archives (Rüggeberg et al., 2007;Roberts et al., 2009). Moreover, the sensitivity of CWCs to climate change is useful to monitor variations in environmental conditions (e.g. water mass variability, surface productivity, bottom current velocity; Rüggeberg et al., 2007;Huvenne et al., 2009;Hebbeln et al., 2016;Wienberg et al., 2018Wienberg et al., , 2020.
The predominant occurrence of CWCs during interglacial periods noticed by Krengel (2020) in the EMCP is also observed during the last four interglacial-glacial cycles in the South Cabliers Mound Province, although mound build-up in the latter is the most important during deglacials and tem-perate interstadials (Corbera et al., 2021). In contrast, over the last 400 kyr, the Tunisian Coral Mound Province situated in the central Mediterranean experienced its most marked phase of coral development during the last glacial period, whereas interglacial periods were particularly scarce in coral occurrences (Corbera et al., 2022). These different temporal distributions suggest that mound build-up in the western and central Mediterranean Sea does not follow a uniform https://doi.org/10.5194/cp-18-1915-2022 Clim. Past, 18, 1915-1945 pattern. Although the development of coral communities at the EMCP during the last 30 kyr is well-documented and novel long-term records are emerging from the western and central Mediterranean Sea (Krengel, 2020;Corbera et al., 2021Corbera et al., , 2022, the long-term environmental forcing affecting the EMCP still remains little documented. The aims of this study are hence (1) to constrain the environmental parameters driving CWC mound formation in the EMCP over the last 300 kyr and (2) to assess the heterogeneities in long-term CWC mound formation within the western Mediterranean Sea.

Geological setting
The Alboran Sea is the westernmost basin of the Mediterranean Sea and is closely connected to the Atlantic Ocean by the Strait of Gibraltar. It is approximately 400 km long, with a width of 200 km, an average depth of 1300 m, and a maximum depth of 1800 m (Olivet et al., 1973;Comas et al., 1999). The Alboran Sea's metamorphic basement is intruded by a number of volcanic plateaus and seamounts formed through the extensional processes that took place between 17 and 8 million years ago (Comas et al., 1999;Duggen et al., 2008). One of these shallow volcanic plateaus, the Provençaux Bank (ca. 200 m depth), extends in a series of three ridges colonized by CWCs named the "Brittlestar ridges" (BRI, BRII, BRIII) (Comas et al., 2009;Fink et al., 2013). They are part of the larger EMCP nestled at depths between 250 and 450 m. The ridges are 3 to 20 km in length and vary in height from 50 to 150 m . These mounds are characterized by dead coral framework with some living corals at their summits and erosional moats at their base, supporting the presence of dynamic currents in the area   (Fig. 1c).

Oceanography
Low-salinity (ca. 36.5 PSU), low-density Atlantic Water enters the Mediterranean through the Strait of Gibraltar. This inflowing water mass mixes with Mediterranean water while crossing the Strait of Gibraltar to form the Modified Atlantic Water (MAW), the dominant surface water mass in the Alboran Sea (La Violette, 1983;Millot, 2009). In addition, evaporation also exceeds river runoff and precipitation; hence, MAW becomes saltier and denser journeying east and finally sinks in the Levantine, Aegean, Adriatic, and Liguro-Provençal sub-basins (Millot et al., 2006). Intermediate waters consist of the highly saline (ca. 38.5 PSU) and warm (ca. 13.5 • C) Levantine Intermediate Water (LIW) that forms in the Levantine basin and flows from east to west, entering the western Mediterranean through the Strait of Sicily to finally exit through the Strait of Gibraltar (Millot, 2013). Levantine Intermediate Water contributes to ca. 70 % of the total out-flow of Mediterranean Outflow Water (MOW; Millot, 2013) and flows between 200 and 600 m water depth in the Alboran Sea, whilst the core of the LIW is situated at approximately 400 m depth (Millot, 2009). It is important to note that, as it moves towards the west, the LIW receives contributions from other water masses, and hence its characteristics gradually change as it gets closer to the Strait of Gibraltar (Millot, 2013). Moreover, intermediate waters appear to differ between the northern and southern Alboran Sea (Fig. 2). Brittlestar Ridge I lays in the depth range of LIW (Fig. 2). Western Mediterranean Deep Water makes up the deepest water mass, flowing under the LIW (Millot and Taupier-Letage, 2005). It forms in the Gulf of Lion and flows westward to finally exit through the Strait of Gibraltar and contributes to the deeper MOW (Millot et al., 2006). In the Alboran Sea, WMDW principally circulates along the Moroccan margin (Millot and Taupier-Letage, 2005).
The surface MAW extends down to approximately 200 m depth (Katz, 1972) and enters the northwestern Alboran Sea as a jet (1.6 Sv; 1 Sv = 10 6 m 3 s −1 ; Lanoix, 1974). This jet triggers the formation of the quasi-permanent anti-cyclonic Western Alboran Gyre that contributes to mixing between surface MAW and underlying LIW (Heburn and La Violette, 1990;Lafuente et al., 1998). When the waters of the Western Alboran Gyre reach the African coast, they separate into two branches: one flows back westward along the coast towards the Strait of Gibraltar, while the other flows towards the eastern part of the basin to form the Eastern Alboran Gyre (La Violette, 1983;Viúdez and Tintoré, 1995). This second non-permanent gyre also contributes to the mixing process between surface and intermediate water masses. The Provençaux Bank and Brittlestar Ridge I are situated in the path of the westward-circulating branch of the Eastern Alboran Gyre (Lanoix, 1974;Viúdez and Tintoré, 1995;Fig. 1). The mixing between surface and intermediate water masses occurs down to ca. 30 m water depth (Heburn and La Violette, 1990). The Strait of Gibraltar is a shallow (ca. 300 m depth) and narrow (ca. 20 km wide) crossing point for entering lower-salinity MAW and exiting higher-salinity MOW (Heburn and La Violette, 1990;Millot, 2009). Thus, the Strait of Gibraltar plays a key role in controlling water mass exchanges between the semi-enclosed Mediterranean Sea and the Atlantic Ocean. The importance of the water exchange varies between glacial and interglacial periods as a function of sea level change. Moreover, the narrow width and depth of the Strait of Gibraltar, together with the geometry of the Alboran basin and the Coriolis force, affects the formation, mean position, and shape of the Alboran gyres (Heburn and La Violette, 1990). Thus, this will in turn affect mixing between surface and intermediate water masses in the Alboran Sea.

Sample collection
This study is based on the multiproxy analysis of gravity core MD13-3462G (35 • 26.531 N, 2 • 31.073 W; 327 m depth; 926 cm long) recovered during the EUROFLEETS cruise MD194 Gateway on board the R/V Marion-Dufresne II (Van Rooij et al., 2017). Cores were split frozen and sedimentary facies descriptions were made at the University of Fribourg prior to sampling. These descriptions include a detailed investigation of the texture, grain size, and colour of the matrix sediment, together with the identification and assessment of the preservation state of major macrofaunal components (Fig. 3). All data were plotted using the ggplot2 package for R (Wickham, 2016;R Core Team, 2018).

Macrofaunal quantification
X-ray computed tomography (CT) imaging was carried out on whole-round sections using a Siemens Somatom Definition AS64 at the Institute of Forensic Sciences at the University of Bern (Switzerland). Core sections were scanned using an X-ray source operating at 120 kV. The images were reconstructed with a slice thickness of 0.6 mm taking into account an increment of 0.3 mm, whilst the pixel resolution of the slices is 0.3 mm. The Avizo 9.4 software was used to visualize, segment, and quantify the volumes of the main macrofaunal components (coral, bryozoan, and bivalve-brachiopod fragments). Prior to segmentation, images were filtered to remove noise in the matrix using a non-local means filter. Bra-chiopods and bivalves were segmented manually. Corals, matrix, pores, and bryozoans were segmented through the combination of dual thresholding and watershed segmentation. Labelled fragments smaller than 5 voxels were filtered prior to quantification. The MaterialStatistics module was used to quantify the vol % of faunal fragments per slice, and the same volume of interest was selected for each core section.

Geochemical logging
Geochemical logging was performed using the Itrax highresolution X-ray fluorescence (XRF) core scanner on split cores at the Institute of Geological Sciences, University of Bern (Switzerland). Measurements were taken at 5 mm intervals using an integration time of 20 s at 30 kV and 45 mA. To counter potentially biased measurements linked to the uneven surface of CWC cores, such as the direct measurement of air or of CWC skeletons, a three-step post-treatment of the dataset was carried out. First, X-ray fluorescence values with argon counts higher than 6000, representing the measurement of air and thereby more porous and/or cracked media not representative of changes in matrix sediment composition, were removed from the final dataset. Secondly, each individual measurement point was compared to high-resolution core images to assess if the measurement was taken on the matrix sediment or not. Finally, elemental counts were normalized by a conservative (minor) element of the background sediment (i.e. aluminium). Aluminium can be used effectively to counter variations in coral content (Löwemark et al., 2011). Normalization of the minor elements with Al is effec- Figure 3. Core description, stratigraphy, and macrofaunal composition of core MD13-3462G. Stratigraphy is based on the planktonic (G. bulloides) and benthic (L. lobatula) δ 18 O records (‰ VPDB), the uranium-series ages of coral fragments, and the epibenthic foraminiferal radiocarbon ages for the first metre of the core (see Fig. 4). Clim. Past, 18, 1915-1945 https://doi.org/10.5194/cp-18-1915-2022 tive when the detrital-terrestrial contribution to the sediment is high. Indeed, aluminium generally behaves conservatively and can hence be used to assess the relative variations of specific elements in sedimentary records (Calvert and Pedersen, 2007, and references therein;Löwemark et al., 2011;Rodrigo-Gámiz et al., 2011;Martinez-Ruiz et al., 2015) In this study, we use the log 10 normalized (Gregory et al., 2019) Si/Al ratio as a proxy for terrestrial (fluvial and aeolian) input. Indeed, the Saharan region is the dominant source of aeolian dust in the Mediterranean Sea and is essentially composed of silicates with high quartz content (Guieu and Thomas, 1996;Caquineau et al., 1998Caquineau et al., , 2002. Moreover, these are rare in Alboran Sea sediments (Masqué et al., 2003), hence the applicability of Si/Al ratio to track variations in terrestrial input. In combination with information provided by benthic foraminiferal assemblages, the log 10 normalized the Si/Al ratio provides a robust and valuable record of terrestrial input.

Grain size analysis and organic geochemistry
The grain size of the siliciclastic fraction was analysed using the Malvern Mastersizer 3000 at the Department of Geology, Ghent University (Belgium). The core was sampled with a small spoon (1 cm 3 ) every 5 cm. Large clasts (> 1 cm), such as coral or bryozoan fragments, were sieved out prior to analysis. Samples were placed in 35 % H 2 O 2 to remove organic matter and boiled until the reaction ended. Following this first step, samples were boiled in 10 % HCl for 2 min to dissolve CaCO 3 . Prior to measurement, samples were placed in 2 % sodium polymetaphosphate and boiled to assure complete disaggregation. Any remaining particle larger than 2 mm was sieved out before measurement. A total of 87 size classes were measured (from 0.01 to 2000 µm). Each sample was measured three times and results were then averaged. The mean grain size of the siliciclastic fraction GS (Folk and Ward, 1957) was calculated for the entire dataset with the Rysgran package for R (Gilbert et al., 2015;R Core Team, 2018). The sortable silt mean size SS, as defined by McCave et al. (1995, i.e. the mean of the 10-63 µm grain size range), was also calculated following the same procedure. Furthermore, following McCave and Hall (2006), the percentage of sortable silt (SS %) in the total <63 µm fraction was calculated. This percentage, together with the sortable silt mean size, was used as an indication of bottom current velocity (McCave and Hall, 2006;Toucanne et al., 2012). It has to be mentioned that the use of SS as a proxy for bottom current velocity on cores recovered from CWC mounds may be biased (e.g. Eisele et al., 2011). Indeed, the baffling effect of the coral framework can locally reduce bottom current velocity and favour the deposition of fine sediments (Huvenne et al., 2009;Titschack et al., 2009;Fentimen et al., 2020b), thus leading to an underestimation of SS during periods with high CWC content. Because of this, only relative increases in SS are considered in combination with results obtained from other proxies.
Total organic carbon (TOC, wt %) and mineral carbon (MinC, wt %) contents were determined on matrix sediments every 10 cm using the Rock-Eval6 technique at the Laboratory of Sediment Geochemistry at the University of Lausanne (Fantasia et al., 2019). The Rock-Eval6 technique produces an oxygen and hydrogen index, respectively corresponding to the quantity of CO 2 relative to TOC and the quantity of pyrolyzable organic compounds relative to TOC (Fantasia et al., 2019). These two indices give an indication of the origin of the organic matter present in the samples (Van Krevelen, 1993).

Microfaunal and macrofaunal investigations
The core was sampled (sliced) every 10 cm for micropaleontological analysis. Samples were weighed dry, washed through a 63 µm mesh sieve, and dried at 30 • C. Each fraction was then dry-sieved through a series of 63, 125, and 2000 µm mesh sieves and weighed. A target number of 300 benthic foraminifera were identified from the fraction larger than 125 µm for each sample. If the residue contained more than 600 specimens, it was split using a dry microsplitter. Relative abundances (%) of benthic species were calculated from the total benthic foraminiferal assemblage. The benthic foraminiferal density was calculated by dividing the total number of foraminifera of a given sample by the sample fraction's weight. The diversity Shannon index (H ) was computed using the PRIMER6 software (Clarke and Gorley, 2006).
Samples prepared for micropaleontological analysis were further used to identify bryozoan species and genera at the Department of Biological, Geological and Environmental Sciences, University of Catania (Italy), on the 12 µm to 2 mm and > 2 mm sized fractions. Key intervals with high bryozoan content, previously identified by CT imagery, were selected. Dominant scleractinian corals and main brachiopod and bivalve species were identified at the lowest taxonomic level possible on the > 2 mm sized fraction at the Department of Geosciences, University of Fribourg (Switzerland).

Radiometric dating
Radiocarbon dating was performed on benthic foraminifera from three samples from the upper first metre of core MD13-3462G at the Laboratory of Ion Beam Physics, ETH Zürich, Switzerland (Table 1). The epibenthic foraminifera species Discanomalina coronata, Lobatula lobatula, and Cibicides refulgens were picked in order to obtain between 4 and 10 mg of pure carbonate. The samples were first dissolved in phosphoric acid. The resulting extracted CO 2 was then converted to graphite and measured by the accelerator mass spectrometry (AMS) technique using the MICADAS dedicated instrument (Synal et al., 2007). Results were corrected for 13 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).
Uranium-series dating was carried out on 24 CWC fragments (D. pertusum and M. oculata) using a multicollector inductively coupled plasma source mass spectrometer (MC-ICPMS; Thermo Fisher Scientific Neptune plus ) coupled with a dissolver (Aridus I) at the Institute of Environmental Physics, Heidelberg University (Table 2). In order to constrain the chronostratigraphy of the core, well-preserved coral fragments were selected at the upper and lower boundaries of coral-rich units. These were identified based on visual core descriptions and CT analysis (macrofaunal quantification; Fig. 3). Coral fragments were physically cleaned with a Dremel ® drill tool as well as by sand-blasting and further chemically cleaned using a weak acid leaching prior to measurements. The detailed sample protocol is described by Frank et al. (2004), while spectrometry and chemical U and Th extraction as well as purification followed Wefing et al. (2017). Uranium-series coral ages were used to calculate mound aggradation rates.

Oxygen and carbon stable isotope analysis
Stable oxygen and carbon isotope compositions were measured on 5 to 12 specimens of the planktonic foraminifera Globigerina bulloides and the benthic foraminifera L. lobatula from the size fraction 212-250 µm in order to prevent any ontogenic effect on the measurements (Schiebel and Hemleben, 2017). The specimens were first cleaned three times with distilled water in an ultrasonic bath for 2 s. 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 (Switzerland) according to the 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 δ 13 C and 0.06 ‰ for δ 18 O values based on eight replicate analyses of standards in each sequence of 40 samples.

Chronostratigraphy
The chronostratigraphy of core MD13-3462G is based on the combination of the coral ages (U-series dating), the planktonic and benthic stable oxygen isotope records, and the foraminiferal radiocarbon ages for the top first metre of the core (Figs. 3 and 4). The U-series coral ages indicate that core MD13-3462G extends approximately from 300 ka BP (Marine Isotope Stage 9) to the Holocene (Figs. 3 and 4, Table 2). Coral ages have been widely used to define the chronology of cores recovered from coral mounds. This approach provides satisfying results, although age reversals down-core have to be taken into account (e.g. Rüggeberg et al., 2007;Frank et al., 2009;Matos et al., 2017). Indeed, reefs are fragile structures and can collapse, topple, and fragment through the action of bioerosion, strong bottom currents, and gravity-driven processes, resulting in transport and redeposition of coral fragments (Beuck et al., 2005;White, 2007). In contrast, constructing a continuous age model based on stable isotope records is generally considered untrustworthy for cores collected from coral mounds since sedimentation is intermittent (e.g. Dorschel et al., 2005). However, coral ages at the upper and lower boundaries of coral build-up phases in core MD13-3462G (e.g. at 390 and 507 cm depth) correspond to changes in the stable oxygen isotope records (Fig. 3), which in turn match the changes between Marine Isotope Stages (MISs; Lisiecki and Raymo, 2005). As such, the stable oxygen isotope records can, in the case of core MD13-3462G and in conjunction with coral ages, indicate important stratigraphic boundaries (Fig. 4). This is particularly relevant during times when CWCs did not grow and hence cannot serve to construct a timeframe. Stable oxygen isotope records were hence correlated with the reference LR04 stack (Lisiecki and Raymo, 2005) for the δ 18 O stratigraphy (see Fig. 4a). Tie points were visually identified and a best correlation coefficient determined using Lineage in the software package AnalySeries v. 2.0.8 (Paillard et al., 1996; Table 3). A clear subdivision into glacial-interglacial stages and substages was possible for Marine Isotope Stages (MISs) 1 to 6 with a Pearson's correlation coefficient of r = 0.79 (p<0.001), in good agreement with the CWC U-series age dates (see Fig. 4b).
However, the correlation became difficult below 650 cm core depth (>150 ka) due to the resolution of sampling (10 cm), the lower sedimentation rate, and possible hiatuses and reworking units flattening the curve. Therefore, U-series ages were used to correlate the lower part of the core, resulting in a Pearson's correlation coefficient of r = 0.51 (p<0.009). The foraminiferal δ 18 O records still follow the LR04 stack until late MIS 7, but the signal remains at relatively light δ 18 O values for the bottom ∼ 100 cm of MD13-3462G covering a time span of around 100 kyr (Fig. 4b).
The stratigraphic boundaries from the base of the core to ca. 650 cm depth were defined based on the U-series coral ages, as planktonic stable oxygen isotope compositions show little variation. The boundaries of MIS 8 are the most poorly defined. Due to difficulties precisely defining the stratigraphy of this section of the core, it will not be considered in detail during this study. In contrast, the planktonic and benthic δ 18 O values and the coral ages do constrain the stratigraphic boundaries from MIS 6 to MIS 1 (Fig. 4). Contrary to sediment records from CWC mounds of the North Atlantic, where no clear glacial or interglacial δ 18 O values are reported (e.g. Dorschel et al., 2005;Rüggeberg et al., 2007;Eisele et al., 2008;Mienis et al., 2009), core MD13-3462G  Table 3) between the LR04 benthic δ 18 O stack of Lisiecki and Raymo (2005) and the benthic (L. lobatula) and planktonic (G. bulloides) δ 18 O record of MD13-3462G. (b) The δ 18 O stratigraphy has good correlation for the younger part (0-150 ka, r = 0.79, p<0.001), but the lower part (>150 ka), which is based on the U-series CWC dates, has a weak correlation (r = 0.51, p<0.009). The comparison between the U-series-and δ 18 O-stratigraphy-based age-depth correlations indicates good coherence. The resulting δ 18 O-stratigraphy-based linear sedimentation rate (LSR) may serve as an indication for changes in the sedimentary regime but shows similar values and/or trends as the CWC-age-based aggradation rate (AR), with higher rates during MIS 3, MIS 6, and late MIS 7. Marine Isotope Stages (MISs) follow boundaries defined by the LR04 stack (Lisiecki and Raymo, 2005). Table 1. Radiocarbon ages of epibenthic foraminifera (species selected: Lobatula lobatula, Cibicides refulgens, and Discanomalina coronata). Ages are corrected for a reservoir age of 390 ± 80 years (Siani et al., 2000).  Table 2. Uranium-series isotope measurements (U / Th) carried out on 24 coral fragments. All errors are 2σ of the mean analytical uncertainty. Ratios were determined using a Th-U spike calibrated to a secular equilibrium reference material (HU-1 at the IUP). Uncorrected, closed-system age was calculated using the decay constants of Jaffey et al. (1971)   presents typical interglacial and glacial δ 18 O values of the Alboran Sea for both planktonic (<1 ‰ and ∼ 3 ‰) and benthic (∼ 1.5 ‰ and ∼ 4‰) foraminifera (e.g. Cacho et al., 1999Cacho et al., , 2006Stalder et al., 2015). Therefore, low planktonic and benthic δ 18 O values correspond to interglacial periods, whilst high planktonic and benthic δ 18 O values correspond to the two last glacial periods (Figs. 3 and 4).

Sediment characterization
The sediment in core MD13-3462G consists mostly of macrofaunal remains (essentially corals and bryozoans) embedded in a clay-to silt-sized carbonate-siliciclastic matrix.
No important variation in the matrix sediment is observed throughout the core. Total organic carbon content in the sediment varies between 0.16 wt % and 1.13 wt % (Fig. 5).
The highest TOC value is measured during late MIS 3 (1.13 wt %), whilst the lowest is recorded during MIS 8 (0.16 wt %; Fig. 5). The most important shifts to higher TOC values are observed during MIS 5, MIS 3, and at the transition between MIS 2 and MIS 1 (Fig. 5). The sediment samples are further characterized by low hydrogen index values (<300 mg HC g −1 TOC; Fig. 6), indicating that the organic matter is oxidized and essentially of terrestrial origin (Espitalié et al., 1985). The mean sortable silt grain size of the siliciclastic fraction (SS) varies between ca. 19 and ca. 26 µm (Fig. 5). Overall, a decrease in SS marks the passage from interglacial to glacial periods. This is particularly noticeable at the transition from MIS 7 to MIS 6, when SS decreases abruptly from approximately 25 to 19 µm (Fig. 5). Conversely, an increasing trend is observed from ca. 550 to ca. 375 cm depth, corresponding to the passage from the later phases of MIS 6 to the end of MIS 5 (Fig. 5).

Stable carbon isotopes and elemental geochemistry
The range of δ 13 C values of the planktonic G. bulloides goes from −2.2 ‰ at 12 cm to −0.5 ‰ at 292 cm, whereas that of the benthic L. lobatula goes from 0.9 ‰ at 872 cm to 1.8 ‰ at 362 (Fig. 5). The planktonic δ 13 C record has higher variability compared to the benthic δ 13 C record (Fig. 5). During MIS 6, the benthic δ 13 C is relatively high (ca. 1.5 ‰), whilst the planktonic δ 13 C record fluctuates between −0.6 ‰ and −1.5 ‰. A decrease in the planktonic δ 13 C record (from −0.7 ‰ to −1.5 ‰) marks the middle of MIS 5. In contrast, the benthic δ 13 C remains stable and low (ca. 1.2 ‰) throughout MIS 5 (Fig. 5). The passage from MIS 4 to MIS 3 is characterized by a shift from the low planktonic δ 13 C recorded during MIS 4 (−1.5 ‰) to higher planktonic δ 13 C (−0.5 ‰). Conversely, benthic δ 13 C values shift from high (1.8 ‰) to lower values (1.3 ‰). The passage from MIS 2 to MIS 1 is marked by a sharp decrease in planktonic and benthic δ 13 C (from −1.2 ‰ to −2.2 ‰ and from 1.8 ‰ to 1.0 ‰, respectively). The last two glacial intervals, in particular MIS 4, are overall marked by more negative 13 C values than during interglacials (Fig. 4). Variations in Si/Al are more marked during MIS 7 and the last glacial period in comparison with the more stable values recorded during MIS 6 and MIS 5. The transitions from MIS 7 to MIS 6 and from MIS 5 to MIS 4 are characterized by fluctuating Si/Al values (Fig. 5).

Macrofauna
The major macrofaunal fragments present in core MD13-3462G are scleractinian corals, bryozoans, brachiopods, and bivalves ( Figs. 3 and 7). Sea urchins, gastropods, serpulids, and gorgonian fragments are more sporadically distributed. Although the dominant coral species in the core is D. pertusum, it is replaced in the upper 20 cm by M. oculata (Figs. 3 and 7). A third and solitary species, Desmophyllum dianthus, is scarcely distributed (Fig. 3). Higher CWC content is observed during interglacial periods (22.2 vol % average), whilst lower content characterizes glacial periods (14.5 vol % average) (Fig. 3). However, coral content shows an uneven distribution during MIS 3, with a range of values from less than 10 vol % to ca. 27 vol % (Fig. 3). Mound aggradation rates ranging between 5 and 18 cm kyr −1 during MIS 3 and early MIS 6 are well in coherence with the linear sedimentation rate based on the foraminifera δ 18 O stratigraphy of the background sediment (Fig. 4b). In contrast, lower mound aggradation rates characterize MIS 5 (ca. 2 cm kyr −1 ) together with MIS 1, 2, and 4 (ca. 4 cm kyr −1 ) (Fig. 4b).  In total 23 genera of bryozoans were identified. Buskea dichotoma is by far the dominant bryozoan species (Fig. 8). Accessory species and genera are mainly represented by Reteporella sparteli, Tubuliporina sp., and Palmiskenea sp. Bryozoan content varies in general between 10 vol % and 20 vol % (Fig. 3). Very high content is, however, observed during MIS 2, reaching nearly 70 vol %. The fragments, although delicate and fragile, are well-preserved, large-sized, and unworn (Fig. 7). Bryozoans are absent during most of MIS 5. This absence corresponds to the time interval when coral content is the most important (Fig. 3). Conversely, the maximum abundance of bryozoans during MIS 2 correlates with a minimum in coral content (Fig. 3).
Brachiopods are mainly represented by the co-occurrence of the species Gryphus vitreus and Terebratulina retusa (Fig. 7). These two brachiopods are regularly associated with the bivalve Bathyarca pectunculoides (Fig. 7). These three invertebrates have been formerly reported from Mediterranean CWC environments. G. vitreus and T. retusa are also recorded from Pleistocene CWC deposits from Rhodes, Greece (Bromley, 2005), whilst B. pectunculoides was found at the Santa Maria di Leuca CWC province (Mastrototaro et al., 2010;Negri and Corselli, 2016). G. vitreus was also found associated with "white corals" between 235 and 255 m depth off the coast of the Hyères Islands, France (Emig and Arnaud, 1988). Although fragile, the shells are wellpreserved (Fig. 7). The brachiopods and bivalves concentrate as layers and demonstrate a non-continuous distribution (Figs. 3 and 7). They reach their highest abundance during glacial periods, in particular at the end of MIS 3 (30 vol % at 80 cm). Brachiopods and bivalves are completely absent during the last two interglacial periods (Fig. 3).

Benthic foraminiferal assemblages
Shannon diversity ranges between ca. 2.8 at 652 cm and 3.6 at 782 cm (Fig. 8). High Shannon diversity values between 3.4 and 3.6 are recorded during interglacial periods (Fig. 8). The lowest Shannon diversity values (between 2.8 and 3.0) are associated with glacial periods (Fig. 8).

Variability of coral mound build-up within the western Mediterranean
Long-term coral mound formation at the location of core MD13-3462G took place during both interglacial and glacial periods (Fig. 4b).
The highest mound aggradation rates of ca. 6 and 10 cm kyr −1 are respectively reached during the middle of MIS 6 and MIS 3, with a short peak of 18 cm kyr −1 during MIS 3. Mound aggradation rates do not exceed ca. 4 cm kyr −1 during interglacial periods and generally range between 1 and 2 cm kyr −1 (Fig. 4b). These rates are comparable to inactive reefs in the Porcupine Seabight (<5 cm kyr −1 ; Frank et al., 2011) and are below the 15 cm kyr −1 threshold set by Frank et al. (2009) for active CWC reef and mound formation, thus suggesting that CWCs did not thrive at the site of core MD13-3462G but rather developed under stressful environmental conditions. The planktonic and benthic δ 18 O values recorded for the last two interglacial and glacial periods, which demonstrate typical interglacial-glacial variations ( Fig. 4; Cacho et al., 1999;Lisiecki and Raymo, 2005;Cacho et al., 2006), are a clear indication that the studied mound location demonstrates a slow albeit continuous build-up history across this time period. Mound aggradation rates for core MD13-3462G are lower than the rates of 17 and 25 cm kyr −1 respectively calculated for MIS 9 and 7 in the neighbouring BRI core GeoB18118-2 (Krengel, 2020) and those of 49 and 83 cm kyr −1 respectively calculated for MIS 7 and 5 in core GeoB18116-2 on Dragon Mound ( Fig. 1c; Krengel, 2020). Similar to the very low mound aggradation rates observed during MIS 5 in core MD13-3462G (ca. 2 cm kyr −1 , Fig. 4b), Krengel (2020) noticed an absence of coral mound build-up during MIS 5 on BRI. With only 9.2 m accumulated over the last ca. 300 kyr against ca. 32 m for the same time period at the site of core GeoB18118-2 (Krengel, 2020), the overall mound aggradation at the site of core MD13-3462G is particularly low. This threefold difference in mound aggradation may suggest that the northern part of BRI (core MD13-3462G) was submitted to more mass-wasting events and/or erosional processes than the southern area (core GeoB18118-2; Krengel, 2020), resulting as such in an overall reduced mound buildup. However, as discussed previously, the benthic and plank-tonic δ 18 O values recorded for the last two interglacial and glacial periods, together with the absence of erosional features downcore, allow us to dismiss this hypothesis. Krengel (2020) also observed that core GeoB18118-2 was in stratigraphic order and only showed minor signs of erosional processes. MD13-3462G (this study) and GeoB18118-1 (Krengel, 2020) are both situated on the crest of BRI at 327 and 332 m depth, respectively, at a distance of ca. 1.3 km (Fig. 1c). Therefore, the observed differences in the timing and rates of long-term mound build-up at BRI are likely driven by local rather than regional and/or basin-wide environmental variability. We propose that disparities in hydrodynamic regimes could lead to reduced food supply and/or low sediment input at the northern part of the ridge where core MD13-3462G was recovered. Supplementary investigations for both cores using bottom current proxies are, however, needed to test this hypothesis. A number of studies have demonstrated that the EMCP and West Melilla Coral Province experienced a rapid phase of mound build-up during the Bølling-Allerød interstadial and the early Holocene (Fink et al., 2013;Stalder et al., 2015;Wang et al., 2019;Fentimen et al., 2020a;Krengel, 2020) with mound aggradation rates varying between 75 and 420 cm kyr −1 . Similar mound aggradation rates between 44 and 203 cm kyr − 1 were calculated by Corbera et al. (2021) at the Cabliers Mound Province. However, in contrast with these observations, mound deposits recovered within core MD13-3462G do not demonstrate such a rapid build-up phase (Fig. 4b). We hypothesize that this deviation from the Alboran Bølling-Allerød to early Holocene mound build-up trend is further evidence that the coral communities situated at the northern part of BRI developed under unfavourable environmental conditions.
In the South Cabliers Mound Province, Corbera et al. (2021) identified four mound formation phases covering the last 400 kyr spread between MIS 9, 7, 6, and 5 at rates of 4, 5, 3.5, and 20 cm kyr −1 , respectively, whereas long-term mound build-up at the Tunisian Coral Mound Province in the central Mediterranean essentially took place during MIS 2 at rates of ca. 20 cm kyr −1 (Corbera et al., 2022). These different mound build-up phases were separated by periods of mound stagnation (Corbera et al., 2021(Corbera et al., , 2022. The contrasting observations made at the Cabliers Mound Province, Tunisian Coral Mound Province, Dragon Mound, and BRI (Krengel, 2020;Corbera et al., 2021Corbera et al., , 2022 suggest that, together with local discrepancies at BRI, the timing of long-term mound build-up in the western and central Mediterranean is not concurrent and does not follow a clear interglacial-glacial pattern as in the North Atlantic (Dorschel et al., 2005;Rüggeberg et al., 2007;Frank et al., 2009Frank et al., , 2011Matos et al., 2015Matos et al., , 2017. The temporal distribution of western and central Mediterranean CWC mounds is rather comparable to mounds situated off the coasts of Angola and Mauritania, where mound build-up took place during both interglacial and glacial periods (Wienberg and Titschack, 2016;Wefing et al., 2017;Wienberg et al., 2018). Overall, the striking disparity in the timing of mound build-up across the western and central Mediterranean hints at strong differences in regional and local environmental forcing.

Glacial mound build-up: a recurrent
Mediterranean trend?
Core MD13-3462G provides the first record of consistent coral growth during the last glacial period in the EMCP and more generally in the Alboran Sea (Fig. 3). Previous observations made by Krengel (2020) at BRI and Dragon Mound evidence very scarce occurrences during MIS 6, with two corals dated at 145.7 and 142.5 ka at BRI and another individual dated at 171.9 ka at Dragon Mound. The Cabliers Mound Province is also characterized by an absence of CWCs during the last glacial, despite demonstrating a phase of mound build-up during MIS 6 (Corbera et al., 2021). Thus, the last glacial occurrence of CWCs in core MD13-3462G stands out, all the more so given that mound aggradation rates reach their highest values during this time (Fig. 4b). This observation contrasts with the complete absence of last glacial coral occurrences in the neighbouring core GeoB18118-1 investigated by Krengel (2020) and demonstrates once again the important heterogeneity in the timing of long-term coral mound build-up along BRI and more generally in the Alboran Sea. Last glacial occurrences of CWCs and mound aggradation rates of 18.3-21.6 cm kyr −1 have recently been reported at the Tunisian Coral Mound Province in the central Mediterranean (Corbera et al., 2022). Coral mound formation is essentially concentrated during MIS 2 (Corbera et al., 2022), unlike at the northern part of BRI (core MD13-3462G) where the most important mound build-up phase occurred during MIS 3 (Fig. 4b). Corbera et al. (2022) argue that increased productivity was a main driver behind this MIS 2 mound formation phase. Likewise, coral growth during the last glacial has been reported from the Gulf of Cádiz and is also suggested to be promoted by increased paleo-productivity linked to strengthened aeolian dust import (Wienberg et al., 2009). Thus, coral growth during the last glacial period spans from the Gulf of Cádiz to the western (EMCP) and central (Tunisian Coral Mound Province) Mediterranean and appears to be a recurrent pattern. The benthic and planktonic foraminiferal δ 18 O and δ 13 C values from core MD13-3462G suggest that environmental conditions were particularly unstable during the last glacial period, as suggested by previous studies Martrat et al., 2004;Pérez-Folgado et al., 2004;Cacho et al., 2006;Bout-Roumazeilles et al., 2007). Moreover, high numbers of the infaunal benthic foraminifera G. subglobosa and C. laevigata (Fig. 8) would indicate that MIS 3 was marked by phases of increased productivity (Schmiedl and Mackensen, 1997;Martins et al., 2006), hence similar to the environmental conditions during the last glacial at the Tunisian Coral Mound Province and Gulf of Cádiz coral mounds.
In contrast with other long-term Mediterranean coral mound records (Krengel, 2020;Corbera et al., 2021Corbera et al., , 2022, the mound deposits situated at the northern part of BRI (site MD13-3462G) show a high contribution of the erect cheleistome bryozoan B. dichotoma (Fig. 3). High abundance of this species during the Bølling-Allerød has previously been reported from the EMCP, where it reached approximately 20 % of the total macrofaunal assemblage (Stalder et al., 2015). Fentimen et al. (2020a) also documented B. dichotoma abundances of up to 30 vol % at the end of the last glacial period at BRI (in core MD13-3455G, see Fig. 1c). With the exception of MIS 5, the mound deposits recovered in core MD13-3462G demonstrate that B. dichotoma was present in numbers throughout the last 300 kyr of mound development and was particularly abundant during the last glacial (ca. 70 vol %; Fig. 3). Based on mound aggradation rates and macrofaunal content, we propose that B. dichotoma communities favoured mound formation at the site of core MD13-3462G, noticeably during the last glacial, by capturing finegrained sediments in a similar way as CWCs do. As such, the investigated mound deposits stand out and may be considered a mixed B. dichotoma-CWC framework rather than a CWC mound per se.

Environmental controls on coral proliferation during the last two interglacial periods
During interglacial periods, benthic foraminiferal assemblages at BRI are marked by high abundances of the infaunal Bulimina spp., U. mediterranea, and to a lesser extent B. spathulata. Several authors describe Bulimina spp. as characteristic for eutrophic and dysoxic environments (Phleger and Soutar, 1973;Lutze and Coulbourn, 1984;Jorissen, 1987;Schmiedl et al., 2000). In the Mediterranean Sea, they are dominant in the vicinity of the Po River delta in the northern Adriatic Sea and close to the Rhône River delta (Jorissen, 1987; Mojtahid et al., 2009). The shallow infaunal U. mediterranea and the opportunistic B. spathulata are known to demonstrate a positive correlation with organic matter flux (De Rijk et al., 2000;Schmiedl et al., 2000;Fontanier et al., 2002Fontanier et al., , 2003Drinia and Dermitzakis, 2010). Moreover, Bulimina spp. and U. mediterranea are reported to be able to feed on fresh but also more refractory organic matter (De Rijk et al., 2000;Koho et al., 2008;Dessandier et al., 2016). Based on these observations, the benthic foraminiferal assemblage during interglacials would support high organic matter export to the seafloor. The overall higher TOC levels during interglacials confirm that the sediment during these periods was relatively enriched in organic matter in comparison to glacial periods (Fig. 5). High abundance of the shallow infaunal G. subglobosa has been linked to the deposition of fresh phytodetritus on the seafloor after bloom events (Gooday, 1993;Fariduddin and Loubere, 1997;Suhr et al., 2003;Sun et al., 2006). It is typically found in highenergy (e.g. steep flanks, ridges) and well-oxygenated envi-ronments (Mackensen et al., 1995;Milker et al., 2009), and it is a common taxon of the Alboran Platform and of CWC environments (Margreth et al., 2009;Milker et al., 2009;Spezzaferri et al., 2014). Mackensen et al. (1995) noted that G. subglobosa dominated in areas of the South Atlantic Ocean where the organic carbon flux did not exceed 1 g m −2 yr −1 . In contrast, in the Mediterranean Sea, B. marginata is restricted to sites with an organic carbon flux >2.5 g m −2 yr −1 , whilst B. aculeata is associated with a flux of 3 g m −2 yr −1 (De Rijk et al., 2000). The last two interglacials (MIS 7 and MIS 5) are marked by an increased abundance of G. subglobosa at early stages followed by a general decline. Buliminids follow a converse trend, particularly during MIS 5, with lower abundances at early stages (Fig. 8). This suggests that conditions during the later stages of interglacials became increasingly eutrophic and in turn less oxygenated at the sediment-water interface, as the consumption of organic matter led to oxygen depletion. These more environmentally stressful conditions resulted in decreased foraminiferal diversity and a proliferation of opportunistic taxa (Fig. 8). Indeed, the lower abundances of Miliolids, which are typically found in well-oxygenated environments (Murray, 2006), further confirm eutrophication coupled to lower oxygenation at the seafloor during interglacials, specifically towards the end of such periods (Fig. 8). Yet, the absence of deep infaunal benthic foraminifera (e.g. Chilostomella spp. or Globobulimina spp.) implies that seafloor oxygenation was never at a minimum, such as during the restricted intervals prior to and after sapropel events in the eastern Mediterranean (Jorissen, 1999;Schmiedl et al., 2003). Schmiedl et al. (2010) link the high abundance of U. mediterranea in the Aegean Sea to humid climatic conditions and increased river runoff. Increased fluvial input has been widely linked in the eastern Mediterranean to more humid continental conditions during interglacial times in response to a northern shift of the African monsoon (e.g. Gasse, 2000;Gasse and Roberts, 2005;Osborne et al., 2008;Coulthard et al., 2013). In contrast, the Alboran Sea lies below the maximum Intertropical Convergence Zone northward position and is sheltered by the Atlas Mountains (Rohling et al., 2002;Tuenter et al., 2003;Lavaysse et al., 2009). Modernday observations show that rainfall over the northwestern Atlas Mountains is generally associated with baroclinic activity over the North Atlantic (Knippertz et al., 2003;Braun et al., 2019). The south of the Atlas Mountains has one of the highest cyclonic activities in the Mediterranean borderlands, whilst the largest fraction of cyclones entering the Mediterranean Sea arrives from the Atlantic (Lionello et al., 2016). Pasquier et al. (2018) noticed that periods of increased input of organic matter from sediment-laden rivers occur during warm substages of the last 200 kyr. These authors relate these pluvial events to negative North Atlantic Oscillationlike conditions (Pasquier et al., 2018). The EMCP is located 50 km away from the mouth of the Moulouya River, which takes its source in the High Atlas Mountains (Snousi, 2004;Emelyanov and Shimkus, 2012;Tekken and Kropp, 2012). The basin of the Moulouya River covers approximately 54 000 km 2 , hence representing the largest river basin in northwestern Africa (Emelyanov and Shimkus, 2012;Tekken and Kropp, 2012). We propose that the influence of warm and moist Atlantic air masses during interglacial periods promoted warmer and more humid conditions over northwestern Africa and torrential rainfall. This would have led to a strengthening of the Moulouya River's flow rate, hence triggering episodes of important terrestrial organic matter input at BRI. These events may have in turn caused eutrophication and oxygen depletion at the seafloor, compatible with the observed benthic foraminiferal assemblages (Fig. 8).
In addition, the variations in the ventilation of the eastern Mediterranean Sea and the formation of LIW may also drive the observed benthic foraminiferal assemblages and reduced mound aggradation rates during interglacial periods. Indeed, at the EMCP Stalder et al. (2015) reported a benthic foraminiferal assemblage of similar composition as the interglacial benthic foraminiferal community within core MD13-3462G. The assemblage described by Stalder et al. (2015) also demonstrates a high abundance of B. marginata, B. aculeata, and C. laevigata. It coincides with periods of D. pertusum absence and the deposition of sapropel S1 in the eastern Mediterranean, hence suggesting that sapropel-related perturbations to the thermohaline circulation triggered oxygen depletion at the seafloor. Sapropel events have also been shown to be concurrent with coral demise in the Cabliers Mound Province, at least since sapropel 7 at the end of MIS 7 (Corbera et al., 2021), and in the southern Adriatic Sea (Taviani et al., 2019). During sapropel events, water column stratification in the eastern Mediterranean led to a dwindling of LIW formation and consequently its reduced circulation in the western Mediterranean Basin (Toucanne et al., 2012;Bahr et al., 2015;Filippidi et al., 2016). Similar to the Cabliers Mound Province (Corbera et al., 2021) and as suggested by benthic foraminiferal assemblages, oxygen depletion linked to reduced LIW circulation and important organic matter input may have resulted in unfavourable and stressful environment conditions for coral growth during interglacial periods, hence explaining the low mound aggradation rates recorded within core MD13-3462G (Fig. 4).
During interglacial periods, the high sea level and the increased evaporation in the Mediterranean led to a more important inflow of low-salinity MAW through the Strait of Gibraltar (Sierro et al., 2005). Thus, surface waters in the Alboran Sea are, in comparison to glacial periods, warmer and less dense. This is also noticed in the planktonic δ 18 O record (Fig. 3). The enhanced MAW flow during interglacials triggers stronger Western and Eastern Alboran Gyres, resulting in better mixing and downwelling. Knowing that the Provençaux Bank and BRI are situated at relatively shallow water depths and in the path of the westward-circulating branch of the Eastern Alboran Gyre (Lanoix, 1974;Viúdez and Tintoré, 1995;Fig. 10), as well as that mixing between surface and intermediate water masses is documented to occur down to ca. 300 m water depth (Heburn and La Violette, 1990), it is conceivable that the corals currently living at 327 m depth were bathed by or situated at the limit of mixing between surface and intermediate water masses during interglacial periods. Wang et al. (2019) suggest that the same phenomenon occurred during the Bølling-Allerød interstadial and the early Holocene. Higher input of MAW into the Alboran Sea would lead to an increased contribution of surface waters to LIW (Fig. 10) and a deepening of the pycnocline. This would promote the formation of internal waves and increase turbulence at the seafloor of BRI, as suggested by the slightly higher SS values at the end of interglacials (Fig. 5), and would have favoured coral proliferation by increasing lateral nutrient supply (Fig. 10). A better mixing of surface and intermediate water masses is suggested by the decreased 13 C and 18 O during the last two interglacial periods (Fig. 5). We suggest that increased food availability during interglacial periods may have enabled coral communities to develop despite oxygen-depleted seafloor conditions in the same way as in the oxygen minimum zones on the Angolan and Namibian margins (Hanz et al., 2019). These conditions would, however, have been detrimental for bryozoan communities.

Environmental conditions during the last two glacial periods
With the exception of MIS 8, for which the boundaries are poorly defined, glacial periods are marked by a change in macrofaunal composition with lower coral and higher bryozoan content in comparison to interglacial periods (Fig. 3). This higher bryozoan content at BRI is in tune with previous observations made at the Great Australian Bight, where bryozoan proliferation during glacial periods promoted the formation of mounds (James et al., 2000;Holbourn et al., 2002). Conversely, higher temperatures and downwelling during interglacials halted bryozoan extension at the Great Australian Bight (James et al., 2000;Holbourn et al., 2002). Rigid erect branching bryozoans such as B. dichotoma are known to be fragile and hence to prefer low-energy environments, being unable to withstand strong bottom currents and turbulence (Scholz and Hillmer, 1995;Bjerager and Surlyk, 2007). Eutrophic environments dominated by infaunal benthic foraminifera (e.g. Bulimina spp.) are unfavourable for erect bryozoans due to the high concentration of suspended food particles clogging up their feeding apparatus (Holbourn et al., 2002). Low SS values and reduced TOC content in the sediment confirm that glacial periods were marked by weaker bottom current velocities and organic matter flux (Fig. 5). The presence of brachiopod and bivalve layers dominated by the brachiopod G. vitreus also characterizes the glacial macrofauna (Fig. 3). This species is found between 160 and 250 m depth along the Mediterranean continental margin and thrives in areas dominated by moderate bottom currents (Emig and Arnaud, 1988). Thus, the co-occurrence of this species with bryozoans confirms that variations in sea level stand, hydrodynamics, and trophic conditions govern the change in macrofaunal dominance at BRI. Low organic matter flux during glacial periods has been related to predominantly arid conditions over North Africa, in association with a weak North African monsoon (Gasse, 2000;Sierro et al., 2005). The reduced precipitation and retreat of vegetation would have triggered a dwindling of terrestrial input during the last glacial period at BRI, as evidenced by a generally lower Si/Al elemental ratio (Fig. 5).
Glacial benthic foraminiferal assemblages are characterized by the dominance of large epibenthic suspensionfeeding foraminifera, such as L. lobatula and D. coronata, together with the infaunal C. laevigata (Fig. 8). This follows observations made by Stalder et al. (2018), who noticed increased abundances of Cibicides spp., D. coronata, and C. laevigata during glacial periods at BRI. These species share a preference for high-quality fresh marine organic matter (De Rijk et al., 2000;Milker et al., 2009;Stalder et al., 2018). Lobatula lobatula and D. coronata have been described to prefer oxygen-rich bottom waters (Linke and Lutze;1993;Margreth et al., 2009). In the Arctic basins and Norwegian-Greenland Sea, the dominance of the epibenthic Cibicides wuellerstorfi (a relative of L. lobatula) reflects a relatively low flux of organic matter (Linke and Lutze;1993) as this species tolerates vertical flux rates <2 g cm −2 yr −1 (Altenbach, 1989). The dominance of L. lobatula, D. coronata, C. laevigata, and Miliolids would thus indicate that the seafloor during glacial periods received less but higherquality organic matter and became more oxygenated in response to the stronger influence of intermediate and deepwater masses (Fig. 10). These observations suggest that more arid conditions during glacial periods led to a reduced influence of terrestrial input on benthic communities (Fig. 10). We propose that weaker but comparatively fresher organic matter input allowed the development of CWC communities, particularly during the last glacial period, and the bryozoan B. dichotoma. This assumption is supported by experimental observations demonstrating how erect bryozoans essentially feed on diatoms and that suspension-feeding foraminifera use the same food sources (Winston, 1977(Winston, , 1981Best and Thorpe, 1994;Goldstein, 1999). It can be hypothesized that there may be a threshold in the quality and quantity of organic matter determining whether D. pertusum or B. dichotoma dominates the benthic environment at BRI. Wang et al. (2019) relate low off-mound GS and high benthic foraminiferal δ 13 C values at BRI during glacials to a dominant influence of MAW coinciding with a low sea level stand. However, whilst the benthic foraminiferal δ 13 C values from core MD13-3462G are indeed relatively high during glacial periods, the planktonic foraminiferal δ 13 C values do not follow the same trend (Fig. 5). The overall low 13 C values during the two last glacial periods, noticeably during  Lanoix (1974), La Violette (1983), and Viúdez and Tintoré (1995). Sea level of interglacial periods corresponds to the current sea level, whilst a 100 m lower sea level stand, following observations made by Rabineau et al. (2006), illustrates glacial periods. The LIW and WMDW flows have been simplified and thus do not represent their exact dynamics. The schematic models are not to scale, although relative depth limits between MAW and LIW have been respected. GEBCO_2019 gridded bathymetric data were used to construct the diagrams. MIS 4, suggest that water mass stratification was greater than during interglacial periods and that the seafloor was not under the direct influence of surface MAW. During glacial periods, the flow of MAW was reduced due to lower sea level and the reduced evaporation over the Mediterranean (Sierro et al., 2005). This would have reduced the contribution of MAW to LIW and weakened the Western and Eastern Alboran Gyres, which would have in turn led to less mixing between surface and intermediate water masses, whilst conversely increasing stratification (Fig. 10). Modern observations show that recently formed dense waters do not necessarily reach the deep western Mediterranean but may, in contrast, be located at intermediate water depths above 1500 m depth (Sparnocchia et al., 1995;Millot, 1999;Ercilla et al., 2016). Ercilla et al. (2016) further revealed that WMDW can be identified at depths shallower than 500 m depth along the Moroccan margin and that it contributes to the overlying LIW, whilst deepwater overturning and ventilation peaked during MIS 2 (Cacho et al., 2006;Toucanne et al., 2012). Increased oxygenation of the seafloor, as evidenced by the benthic foraminiferal assemblage (Fig. 8), may suggest that the contribution of well-ventilated deep and intermediate water masses at BRI was more important during glacials than during interglacials (Fig. 10). The physical shape and structure of BRI possibly plays a role in the shoaling of deep waters during glacial periods. In addition, the overall higher benthic 13 C values and the abundance of fresh-organic-matter-feeding foraminifera (L. lobatula and D. coronata) during glacial periods could indicate that these waters were also nutrient-rich. Although stratification between surface and intermediate water masses was greater during glacials, the stronger flow of well-ventilated WMDW at BRI would explain the higher oxygen availability at the seafloor. Overall during glacial periods, and in particular during the Last Glacial Maximum (LGM), enhanced contribution of nutrition-rich and well-ventilated WMDW to overlying LIW would have promoted mound aggradation.

Conclusions
The multiproxy study of core MD13-3462G reveals that mound build-up at the northern part of Brittlestar Ridge I (East Melilla Coral Province, SE Alboran Sea) took place during both interglacial and glacial periods. A number of key observations can be underlined.
1. Average coral mound aggradation rates are particularly low, varying between 1 and 10 cm kyr −1 , whilst maximum aggradation rates are recorded during MIS 3 (18 cm kyr −1 ). These rates suggest that corals never thrived in this sector of Brittlestar Ridge I but rather developed under stressful environmental conditions. We propose that weak bottom-water oxygenation linked to sapropel-related events and/or increased precipitation over North Africa led to the slow development of coral communities during interglacial periods. Intensified cir-culation of Levantine Intermediate Water and the import of fresh organic matter would have provided suitable conditions for bryozoan and coral communities during glacial periods.
2. Core MD13-3462G provides the first record of consistent coral growth during the last glacial period in the East Melilla Coral Province and more generally in the Alboran Sea. This conspicuous observation, in conjunction with other records of Mediterranean longterm coral mound build-up, suggests that coral mound development does not follow a clear-cut interglacialglacial pattern in the western Mediterranean. Furthermore, regional-and local-scale environmental variability appears to play a decisive role in mound build-up in the eastern Alboran Sea.
3. The planktonic and benthic δ 18 O records of cold-water coral mound sediments at Brittlestar Ridge I show typical interglacial-glacial variations since early MIS 6. This is in contrast with δ 18 O records generally recovered from coral mound deposits and highlights that the northern part of Brittlestar Ridge I experienced reduced albeit continuous build-up.
From a wider perspective, the build-up of cold-water coral mounds situated at Brittlestar Ridge I during both interglacial and glacial periods stresses how cold-water coral communities are capable of withstanding important environmental changes as well as surviving and adapting to different climatic conditions. This study further shows that the role of associated species, such as rigid erect bryozoans, may be linked to the resilience of coral ecosystems.  (Fentimen et al., 2020c).

Sample availability.
Archive halves of all core sections investigated for this study are available at the Department of Geosciences, University of Fribourg (Switzerland). The sediment residues and the splits of each sample analysed for benthic foraminiferal assemblages are stored at the Department of Geosciences, University of Fribourg (Switzerland). Bryozoans identified in this study are available at the Palaeontological Museum of the University of Catania (Italy).
Author contributions. RF was responsible for writing (original draft), visualization, conceptualization, core sampling, and investigation (benthic foraminiferal assemblages, main macrofaunal fragments, particle size analysis, stable isotope measurements assisted by TV, and radiocarbon dating assisted by IH). EF was responsible for conceptualization, writing (review and editing), XRF inveshttps://doi.org/10.5194/cp-18-1915-2022 Clim. Past, 18, 1915-1945 tigation (assisted by HV), and preparation of samples for uraniumseries dating and Rock-Eval6 pyrolysis. ARü was responsible for conceptualization, writing (review and editing), and supervision. EH was responsible for investigation (CT analysis, macrofaunal quantification). VR was responsible for writing (review and editing) and visualization. TV was responsible for writing (review and editing), investigation (stable isotope measurements), and resources. IH was responsible for writing (review and editing), investigation (radiocarbon dating), and resources. ARo was responsible for writing (review and editing) and investigation (bryozoan taxonomy). DVR was responsible for writing (review and editing) and resources. TA was responsible for writing (review and editing), investigation (Rock-Eval6 pyrolysis), and resources. HV was responsible for writing (review and editing), investigation (XRF), and resources. NF was responsible for writing (review and editing) and investigation (uranium-series dating). AF was responsible for investigation (core description, CT data analysis, XRF data analysis), conceptualization, writing (review and editing), project administration, funding acquisition, and supervision.