the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.

Impact of deoxygenation and hydrological changes on the Black Sea nitrogen cycle during the Last Deglaciation and Holocene
Nicole Bale
Rick Hennekam
Bingjie Yang
Darci Rush
Gert-Jan Reichart
Ellen C. Hopmans
Stefan Schouten
The marine nitrogen (N) cycle profoundly impacts global ocean productivity. Amid rising deoxygenation in marine environments due to anthropogenic pressures, understanding the impact of this on the marine N cycle is vital. The Black Sea's evolution from an oxygenated lacustrine basin to an anoxic marine environment over the Last Deglaciation and Holocene offers insight into these dynamics. Here, we generated records of the organic biomarkers heterocyte glycolipids (HGs), crenarchaeol, and bacteriohopanetetrol, associated with various water column microbial N-cycle processes, which indicate a profound change in Black Sea N-cycle dynamics at ∼7.2 ka when waters became severely deoxygenated. This transition substantially reduced Thaumarchaeota-driven nitrification and enhanced loss of bioavailable nitrogen through anaerobic ammonium oxidation (anammox). In contrast, other climatic changes over the Last Deglaciation and Holocene, such as freshwater input, water-level variations, and temperature changes, did not impact these processes. Cyanobacterial nitrogen fixation in surface waters proved more responsive to changes in salinity, which affected species composition, and associated water column stratification, which reduced the vertical transport of nutrients. Our results indicate that future deoxygenation in certain marine environments may enhance bioavailable nitrogen loss by anammox and reduce nitrification by Thaumarchaeota, while enhanced stratification may increase cyanobacterial nitrogen fixation in the surface waters.
- Article
(4648 KB) - Full-text XML
-
Supplement
(1492 KB) - BibTeX
- EndNote
The marine nitrogen (N) cycle is a significant control of biological productivity in the global oceans. It is directly connected to the fixation of atmospheric carbon dioxide (CO2) and carbon export from the ocean's surface, influencing atmospheric CO2 levels over geological timescales (Falkowski et al., 1998). As the marine N cycle is strongly regulated by biology, the (de)oxygenation of the ocean determines the microorganisms involved in these biogeochemical cycles and the aerobic/anaerobic pathways that occur. Under anoxic conditions, the loss of bioavailable nitrogen is substantial, attributed to anaerobic ammonium oxidation (anammox) and denitrification (Kuypers et al., 2003; Dalsgaard et al., 2012). With deoxygenation in marine environments increasing due to anthropogenic climate and environmental changes (i.e. Keeling et al., 2010; Bopp et al., 2013) and with research linking deoxygenation to changes in the marine N cycle (Kalvelage et al., 2013; Naafs et al., 2019), it is important to enhance our understanding of how the marine N cycle may respond to future deoxygenation and what the associated feedbacks on carbon fixation might be.
Marine basins that have experienced changes in oxygenation in the past can provide perspective on the current deoxygenation of modern global oceans and the associated feedbacks in the marine N cycle, particularly on timescales beyond the observational record. Today, the Black Sea is the world's largest permanently stratified anoxic basin, with limited connection to the global ocean through the Bosporus Strait, and its redox gradient is a hotspot of diverse microbial populations and metabolisms, experiencing many of the same crucial microbial biogeochemical cycle processes as the global ocean (Kusch et al., 2022). However, over the Last Deglaciation and Holocene (approximately the last 20 ka), the Black Sea has experienced large hydrological changes. The basin was an oxygenated freshwater lacustrine environment during the Last Glacial Maximum (LGM) (Schrader, 1979) and experienced many environmental changes during the subsequent deglaciation, including temperature changes (Bahr et al., 2005, 2008; Ion et al., 2022), water-level variations (Ivanova et al., 2007; Nicholas et al., 2011; Piper and Calvert, 2011), and changes in freshwater input into the basin, through the melting of both Eurasian ice sheets and alpine glaciers after the LGM and changes in regional precipitation (Bahr et al., 2005, 2006, 2008; Badertscher et al., 2011; Shumilovskikh et al., 2012). It became reconnected to the global ocean at ∼9.6 ka when post-glacial sea-level rise caused an initial marine inflow (IMI) over the Bosporus sill (Aksu et al., 2002; Major et al., 2006; Bahr et al., 2008; Ankindinova et al., 2019), leading to enhanced salinity of the water column (Marret et al., 2009; Verleye et al., 2009; Filipova-Marinova et al., 2013) and euxinic deep waters developing in the basin after ∼7.2 ka (Arthur and Dean, 1998; Eckert et al., 2013). Thus, sedimentary records of the Black Sea may provide a unique perspective of the impacts of deoxygenation, changing temperature, and salinity on the marine N cycle.
Diagnostic lipid biomarkers of microbes preserved in the geological record can offer a unique insight into past changes in the N cycle (Rush and Sinninghe Damsté, 2017, and references cited therein; Elling et al., 2021; van Kemenade et al., 2023), although the limitations of extrapolating modern findings to ancient climates must be acknowledged, as past ecosystems may have operated under different dynamics that are not fully captured by contemporary analogues. Nitrogen-fixing heterocytous cyanobacteria play a crucial role in transforming nitrogen gas (N2) into bioavailable nitrogen (NH3) and sustaining primary productivity in both marine and freshwater environments (Villareal, 1992; Ploug, 2008). Identification of their diagnostic biomarkers, heterocyte glycolipids (HGs), in the geological record are a widely used proxy for exploring past changes in nitrogen fixation by these microbes (Bauersachs et al., 2009, 2010; Sollai et al., 2017; Bale et al., 2019; Elling et al., 2021; Pérez Gallego et al., 2025). The structure of the sugar moiety of HGs is a useful indicator of palaeoenvironmental conditions, as HGs with a hexose (C6) headgroup are typically found in free-living heterocystous cyanobacteria (Bauersachs et al., 2009; Wörmer et al., 2012), while HGs with a pentose (C5) headgroup are associated with endosymbiotic heterocystous cyanobacteria in marine diatoms (diatom–diazotroph associations, DDAs) (Schouten et al., 2013; Bale et al., 2015). Consequently, both hexose and pentose HGs have been applied as specific palaeo-biomarkers for the presence of N2-fixing cyanobacteria in marine and lacustrine geological records (Bauersachs et al., 2010; Sollai et al., 2017; Bale et al., 2019; Elling et al., 2021). Nitrification, the microbial two-step conversion of ammonia (NH3) and/or ammonium () to nitrate (), is a central part of the marine N cycle. Archaea of the phylum Thaumarchaeota (also known as Nitrososphaerota) are among the most abundant and widespread marine prokaryotes (Karner et al., 2001; Francis et al., 2005), playing a crucial role in nitrification in the Black Sea (Lam et al., 2007) by aerobically oxidizing ammonia to nitrite (Könneke et al., 2005; Wuchter et al., 2006). As Thaumarchaeota are the exclusive producers of the membrane-spanning lipid, crenarchaeol (Sinninghe Damste et al., 2002), this biomarker can be used to identify Thaumarchaeota in the geological record and explore the palaeo-marine N cycle. Another critical part of the N cycle is the loss of bioavailable nitrogen to N2. Under anoxic conditions, bioavailable nitrogen (, , NH3, and ) can be lost through two processes in subsurface waters: anammox (van de Graaf et al., 1997; Kuypers et al., 2003) and denitrification (Kuenen and Robertson, 1988). It is possible to explore past changes in anammox activity in the sedimentary record using the unique ladderane fatty acids (Sinninghe Damste et al., 2002), but these are relatively poorly preserved in sediments (Jaeschke et al., 2007). Alternatively, the ratio of bacteriohopanetetrol (BHT)-34S (which is ubiquitously synthesized by aerobic bacteria) to the later-eluting stereoisomer BHT-x (which is predominately synthesized by marine anammox bacteria, i.e. Ca. Scalindua spp.) (Rush et al., 2014; Schwartz-Narbonne et al., 2020; van Kemenade et al., 2023) can be used to trace past anammox activity. Denitrification is performed by a large range of organisms (Knowles, 1982), but, at present, there are no associated diagnostic lipid biomarkers (Rush and Sinninghe Damsté, 2017).
In this study, we used lipid biomarkers of microbes widely involved in the N cycle across various marine and freshwater environments in combination with other geochemical records from a sediment core located in the western Black Sea spanning the Last Deglaciation and Holocene (∼20 ka–present) to better constrain and assess the sensitivity of the marine N cycle under changing hydrological and oxygenation conditions and explore its potential links to broader global climate dynamics.
The Black Sea is a large meromictic marginal basin connected to the Mediterranean Sea via the Turkish Straits (the Bosporus Strait, the Sea of Marmara, and the Dardanelles Strait) (Fig. 1). The Black Sea has a net outflow into the Aegean Sea via the Turkish Straits and is primarily supplied by three major rivers: the Danube, the Dnieper, and the Don. With freshwater flowing out of the basin and dense, highly saline waters flowing in, the water column is highly stratified with respect to salinity (density). An oxygenated colder surface layer (0–50 m) overlies warmer, anoxic, sulfidic, hypersaline deep waters (100–2300 m), separated by a suboxic layer (50–100 m) (Murray et al., 1989, 1995). The general circulation of Black Sea surface waters is a basin-scale cyclonic boundary current encompassing large eastern and western cyclonic gyres, with several smaller, anticyclonic coastal eddies (Fig. 1) (Özsoy and Ünlüata, 1997). The modern Black Sea water column is characterized by the presence of strong redox gradients. In the basin's western gyre, the oxic zone (0–75 m depth range) has an oxygen concentration of ∼121 µmol kg−1 and salinity of 19.4 psu at 50 m depth, with no detectable sulfide concentrations (Sollai et al., 2018; Bale et al., 2021). The suboxic zone lies below (75–115 m depth range), where salinity increases with depth, and traces of sulfide are detected at the bottom of this layer at ∼110 m (Sollai et al., 2018; Bale et al., 2021). Beneath this lies the euxinic zone (115–2000 m depth range), where salinity and sulfide concentrations increase with depth, reaching 22.3 psu and ∼400 µmol L−1 at a depth of 2000 m (Sollai et al., 2018; Bale et al., 2021). Ammonium concentrations are low in the oxic and upper suboxic zones (<0.1 µmol L−1), increasing below 90 m, with highest concentrations (∼100 µmol L−1) at 2000 m water depth (Sollai et al., 2018). Nitrite concentrations are highest in the oxic zone, peaking at 50 m (∼0.08 µmol L−1), decreasing through the suboxic zone, with the exception of a peak at 85 m, while, in the euxinic zone, nitrite is below the limit of detection (Sollai et al., 2018). Nitrate concentrations are highest in the lower oxic and upper suboxic zones, peaking between 70–80 m (∼2.5 µmol L−1), decreasing with depth, and reaching the limit of detection below 105 m (Sollai et al., 2018; Bale et al., 2021).
During the cruise with the RV Pelagia in April 2017, piston core 64PE418 (235 cm length) was recovered from 1970 m below sea level (m b.s.l.) depth in the Black Sea (42°56 N, 30°02 E) (Fig. 1). A total of 44 sediment samples were taken at 5 cm intervals along the depth of the core.
3.1 Biomarker extraction and analysis
Lipids were extracted from these samples using a modified Bligh and Dyer extraction method as described previously (Bale et al., 2021). Using a mixture of methanol (MeOH), dichloromethane (DCM), and phosphate buffer (, v:v), the sediment was twice extracted ultrasonically (10 min). The combined supernatants were phase-separated by adding DCM and phosphate buffer to create a solvent ratio of (v:v). The organic phase was collected, and the aqueous phase was re-extracted three times using DCM. All extraction steps were then repeated on the residue but with a mixture of MeOH, DCM, and aqueous trichloroacetic acid solution (TCA) pH 3 (, v:v). Finally, the organic extracts were combined and dried under an N2 gas stream. A deuterated betaine lipid {1,2-dipalmitoyl-sn-glycero-3-O-4'-[N,N,N-trimethyl(d9)]-homoserine; Avanti Lipids} internal standard was added to each sample before filtering the extract through 0.45 µm cellulose syringe filters (4 mm diameter; BGB, USA). Extraction blanks were performed alongside the sediment extractions, using the same glassware, solvents, and extraction methodology, but without sediment. Analysis of the extracts was performed using the following UHPLC-HRMS reversed-phase method. An Agilent 1290 Infinity I UHPLC instrument was used, equipped with a thermostatted auto-injector and column oven, and coupled to a Q Exactive Orbitrap MS instrument with an Ion Max source using a heated electrospray ionization (HESI) probe (Thermo Fisher Scientific, Waltham, MA). Separation was achieved using an Acquity BEH C18 column (Waters, 2.1 mm×150 mm, 1.7 µm) maintained at 30 °C. The eluent composition was (A) MeOH/H2O/formic acid/14.8 M NH3aq [ (v:v)] and (B) IPA/MeOH/formic acid/14.8 M NH3aq [ (v:v)]. The elution programme was 95 % A (for 3 min) followed by a linear gradient to 40 % A (at 12 min) and then to 0 % A (at 50 min), which was maintained until 80 min. The flow rate was 0.2 mL min−1. Positive-ion HESI settings were as follows: capillary temperature, 300 °C; sheath gas (N2) pressure, 40 arbitrary units (AU); auxiliary gas (N2) pressure, 10 AU; spray voltage, 4.5 kV; probe heater temperature, 50 °C; S lens, 70 V. Lipids were analysed with a mass range of 350–2000 (resolving power 70 000 ppm at 200), followed by data-dependent tandem MS/MS (resolving power 17 500 ppm), in which the 10 most abundant masses in the mass spectrum were fragmented successively. Optimal fragmentation was achieved with a stepped normalized collision energy of 15, 22.5, and 30 (isolation width, 1.0 ) for IPL analysis (Bale et al., 2021) and 22.5 and 40 (isolation width 1.0 ) for BHP analysis (Hopmans et al., 2021). The Q Exactive was calibrated within a mass accuracy range of 1 ppm using the Thermo Scientific Pierce LTQ Velos ESI Positive Ion Calibration Solution. During analysis, dynamic exclusion was used to temporarily exclude masses (for 6 s) to allow the selection of less abundant ions for MS/MS.
Biomarkers were identified based on their retention time, exact mass, and fragmentation spectra. Integrations were performed on (summed) mass chromatograms of relevant molecular ions (, , and ) and, in the case of crenarchaeol, also the second isotope peak for each of the three adducts. Due to the coelution of BHT-34S, BHT-x isomer, and an unknown nitrogen containing a compound with the same mass, identification and integration of BHT-34S and BHT-x were conducted using the 529.462 dehydrated in-source product (-H2O). Isoprenoidal glycerol dialkyl glycerol tetraether (isoGDGT) crenarchaeol, monohexose crenarchaeol, and a crenarchaeol isomer were all integrated and combined as “crenarchaeol”. The lipid biomarker records are all presented as peak area per gram of total organic carbon (TOC).
3.2 Total organic carbon, total nitrogen, and δ15Nbulk measurements
Freeze-dried sediments were analysed for TOC, total nitrogen (TN), and bulk δ15N (δ15Nbulk) using a Thermo Scientific Flash EA DELTA V Plus IRMS instrument. Flow was 100 mL min−1, and the temperatures for oxidation, reduction, and the oven were 900, 680, and 40 °C, respectively. Nitrogen isotopic measurements were calibrated to atmospheric air (AIR), and values are expressed in per mille (‰) units. Inorganic carbon was removed from the sediment prior to TOC analysis using HCl (2 mol), cleaned with bi-distilled water, then freeze-dried.
3.3 Age model
Accelerator mass spectrometry (AMS) 14C ages of bulk organic matter were measured from core 64PE418 (n=7) to create a chronology on the 64PE418 depth scale. Samples were weighed and freeze-dried at NIOZ. The AMS 14C measurements () were determined using a Compact Carbon AMS system at the Poznań Radiocarbon Laboratory, Poland. The sediment samples were pre-treated with 0.25 M HCl (room temperature overnight, then 80 °C, 1+ h) and rinsed with deionized water until reaching a pH of 7. Samples were then combusted in closed (sealed under vacuum) quartz tubes, together with CuO and Ag wool (900 °C, 10 h). The CO2 released was then dried in a vacuum line and reduced with H2 using 2 mg of iron (Fe) powder as a catalyst. The obtained carbon and Fe mixture was then pressed into an aluminium holder (Czernik and Goslar, 2001). The measurement was performed by comparing intensities of ionic beams of 14C, 13C, and 12C measured for each sample and for standard samples (with “Oxalic Acid II” used as modern standard and “coal” used as background standard of 14C-free carbon). In each AMS run, 30–33 samples of unknown age were measured, alternated with measurements of 3–4 samples of modern standard and 1–2 samples of background standard. The measured ratios were corrected for isotopic fractionation and reported as conventional radiocarbon age according to Stuiver and Polach (1977).
Seven bulk organic matter 14C dates were used in the production of the age model for core 64PE418 (Table S1 and Fig. S1 in the Supplement). Six of these were from this core, with an additional bulk organic carbon 14C date from the widely acknowledged Unit I/II boundary of core KNR 134-08 BC17, which was used to further refine the age model for the upper part of the core (Jones and Gagnon, 1994). Core KNR 134-08 BC17 was sourced from the same location and water depth as 64PE418, and this boundary was identified in our core using the same significant colour and elemental changes described in previous studies (Figs. 2 and 3) (i.e. Arthur and Dean, 1998; Bahr et al., 2005). While seven 14C measurements were conducted on core 64PE418, one was excluded from the age model due to an age reversal (142.5 cm), likely due to the presence of reworked material, which would have made the age model out of line with other studies (Fig. S2 in the Supplement). Variable reservoir ages were added to our calibration (Table S1), using those calculated by Kwiecien et al. (2008) for intermediate water depths in the Black Sea over the Last Deglaciation and Holocene which were deemed the most suitable for our site location. The 14C dates were calibrated using the Marine20 calibration curve (Heaton et al., 2020) for the upper three samples (24.5, 39, 76.5 cm), which reflect the period after the infiltration of marine water; this is based on the colour and elemental changes in the core which indicate that these samples fall within Units I and II (Arthur and Dean, 1998; Bahr et al., 2005). The lower four samples (118.5, 158.5, 183.5, and 217.5 cm) were calibrated using the IntCal20 calibration curve (Reimer et al., 2020), as they reflect the period prior to the marine infiltration when the Black Sea was a lacustrine environment, as indicated by colour and elemental signatures in the core (Arthur and Dean, 1998; Bahr et al., 2005). Using the R code CLAM (Blaauw, 2010), the age–depth model was created based on the seven 14C dates. Our age model shows that the 64PE418 biomarker records span the last 19.5 ka, with an average resolution of ∼450 years. The following transitions are identified in our core by colour (Fig. 2) and elemental changes (Fig. 3) and are dated by our age model as follows: the onset of the IMI (138 cm) is at 9.6 ka ± 237 yr, the boundary of Unit II/III (96 cm) is dated at 7.2 ka ± 202 yr, and the Unit I/II boundary (39 cm) is dated at 2.6 ka ± 402 yr. Despite the complexity of dating Black Sea cores, due to lack of suitable material for dating and dispute over reservoir age corrections (Kwiecien et al., 2008; Soulet et al., 2011; Yanchilina et al., 2017), the dates of these boundaries align well with previously published calibrated ages for these transitions (i.e. Jones and Gagnon, 1994; Ankindinova et al., 2019; Huang et al., 2021), as shown in Fig. S2.

Figure 2Scan of core 64PE418 showing colour changes and the depth of unit boundaries. Unit boundaries are defined according to Arthur and Dean (1998) and are identified by colour changes and XRF core scan changes in Ti and Ca (Fig. S1).
4.1 TOC, TN, and colour changes
Sedimentary bulk TOC (%), bulk TN (%), and δ15Nbulk (‰) range between 0.3 %–22.8 % for TOC, 0.05 %–1.9 % for TN, and 5.2 ‰–0.0 ‰ for δ15Nbulk (Fig. 4). There are significant colour changes in the core, as shown in Fig. 2, which correspond to changes in TOC, TN, and the elemental composition (Fig. 3). In the lower part of the core (19.5–9.6 ka), values are relatively low for TOC and TN, at ∼0.84 % and ∼0.10 %, respectively. At 9.6 ka, there is an appreciable change in the elemental composition of the core, with increases in , K, and V and a decrease in , which corresponds with a transition to darker sediments and an increase in TOC and TN to ∼2.41 % and ∼0.26 %, respectively. At 7.2 ka, there is another major change in the colour and bulk elemental composition of the core, with an increase in redox-sensitive elements U, V, and Mo and a decrease in and K (Fig. 3), which corresponds with darker sediments and increasing TOC values. TOC peaks between 6.6–4.6 ka (∼21 % for TOC and ∼1.7 % for TN), declining towards the top of the core. δ15Nbulk shows a general decline in values from the upper to the lower part of the core. This decline is small between 19.5–7.7 ka (4.9 ‰–3.3 ‰), before a more significant decrease to 1.2 ‰ at 6.6 ka (3.3 ‰–1.2 ‰). Values increase to 3.7 ‰ at 6.1 ka before declining to 0.0 ‰ at 3.9 ka, increasing slightly towards the top of the core to values of 1.3 ‰.
4.2 Biomarkers
We examined a number of lipid biomarkers related to the N cycle in Black Sea core 64PE418 (Fig. 4). HGs were identified in all samples (with the exception of 215 cm (16.4 ka)). These include HGs with a hexose (C6) headgroup, i.e. hexose C26 diol, hexose C28 diol, hexose C28 triol, and hexose C30 triol, which are specific to free-living cyanobacteria found in predominately freshwater and brackish environments (Bauersachs et al., 2009; Wörmer et al., 2012). In addition, we detected those with a pentose (C5) headgroup, i.e. pentose C30 diol, pentose C30 triol, and pentose C32 triol, which are specific to cyanobacteria symbiotic with marine diatoms (Schouten et al., 2013; Bale et al., 2015). Hexose HGs are present throughout the core, increasing substantially in abundance between 9.6–6.6 ka, reaching maximum values at 9.6 ka. Pentose HGs are detected from 4.3 ka onwards, increasing in abundance at the top of the record and coinciding with a low abundance of hexose HGs. Crenarchaeol, a marker for Thaumarchaeota, was identified throughout our record, showing high values in the early part of the record ( peak area per g TOC) until 6.9 ka, abruptly shifting to lower values of peak area per g TOC thereafter. The BHT-x ratio, a biomarker for anammox bacteria, is low in the early part of our record (<0.3) due to a low abundance of BHT-x. The BHT-x ratio increases after 6.9 ka to values around 0.3 due to a higher abundance of BHT-x and a lower abundance of BHT with a 34S stereoconfiguration.
Finally, to reconstruct levels of oxygen in the subsurface waters of the Black Sea, isorenieratene was identified (as described in Bale et al., 2021). Isorenieratene is a marker of the brown strains of the photosynthetic green sulfur bacteria, Chlorobiaceae, which are anoxygenic photoautotrophs that require light and hydrogen sulfide (H2S); their presence indicates photic zone euxinia, whereby anoxic sulfidic waters reach the photic zone (Sinninghe Damste et al., 1993; Koopmans et al., 1996). Isorenieratene was identified in many of our samples after 9.5 ka, peaking between 5.6–4.3 ka (reaching 3.39×1012 per g TOC at 5.6 ka), but was not detected between 3.9–2.7 ka.
Based on clear changes in TOC (Fig. 4), colour, and elemental signatures (Figs. 2 and 3), we divided core 64PE418 into three widely acknowledged units, in line with previous studies (Jones and Gagnon, 1994; Arthur and Dean, 1998; Bahr et al., 2005). Unit III spans ∼20–7.2 ka, covering the period where the Black Sea was a lacustrine environment, disconnected from the global ocean, and also the transition interval, where the basin moved towards a marine environment after the IMI over the Bosporus sill at ∼9.6 ka (Aksu et al., 2002; Major et al., 2006; Bahr et al., 2008; Ankindinova et al., 2019). Unit II (∼7.2–2.6 ka) and Unit I (∼2.6 ka–present) span the period where the Black Sea had become an anoxic brackish-to-marine environment. Diagenesis and preservation biases can affect lipid biomarker records, but, since the records do not all mirror the oxygenation of the water column and since organic carbon contents remained relatively high, it is unlikely that diagenesis and preservation significantly influenced these records.
5.1 Oxic lacustrine phase (19.5–9.6 ka)
Throughout the Last Deglaciation and early Holocene (19.5–9.6 ka), TOC and TN levels are low, likely due to poor preservation of organic material, caused by the well-ventilated, oxygenated freshwater environment that existed in the basin at this time (Schrader, 1979). Isorenieratene is not detected during this period, while elements that accumulate in sediment under anoxic conditions (i.e. Algeo and Li, 2020) also remained low (i.e. U, V, Mo; Fig. 2), and this points to a well-oxygenated environment. Freshwater/brackish conditions prevailed throughout this time, as shown by previous studies (Fig. S3 in the Supplement; Filipova-Marinova et al., 2013; Ion et al., 2022; Huang et al., 2021). Throughout this period, the abundance of Thaumarchaeota, indicated by crenarchaeol abundance, and anammox, indicated by the BHT-x ratio, remained relatively steady (Fig. 4). This stability is remarkable, since the region experienced significant climatic changes which led to large variations in the surface water temperatures of the Black Sea, varying from ∼10 °C during the Bølling–Allerød, ∼7 °C during the Younger Dryas and ∼14 °C by the early Holocene (Ménot and Bard, 2012), along with changes in the input of freshwater into the basin due to regional precipitation variability and the melting of Eurasian ice sheets and alpine glaciers (Bahr et al., 2005, 2006, 2008; Badertscher et al., 2011; Shumilovskikh et al., 2012; Filipova-Marinova et al., 2013; Ion et al., 2022). The subsurface stability is likely due to the stratification of the basin, where significant climatic shifts primarily impacted the surface waters, and the limited vertical mixing minimized the influence on subsurface waters, creating a stable environment for the subsurface nitrogen cycle. In contrast, changes in HG abundance and distribution suggest that surface-dwelling nitrogen-fixing cyanobacteria were sensitive to surface water hydrological changes in the Black Sea over this period (Fig. 5). The dominant HG structure varies between hexose C26 diol, hexose C28 diol, and hexose C30 triol, and, after 11 ka, hexose C28 triol becomes present, which has been shown to be the major HG in members of the Rivulariaceae family (i.e. Calothrix sp.) (Bauersachs et al., 2009). The warmer, wetter conditions of the early Holocene may have provided a trigger for this change in HG abundance and composition. Indeed, an increase in the abundance of the genus Rivularia was also noted in coastal regions of SW India during this period, coinciding with an increasingly warm and wet climate (Limaye et al., 2017). Another cause for this shift may have been related to changes in nutrient availability, with members of the Rivulariaceae family typically occurring in environments with highly variable phosphorus availability (Whitton and Mateo, 2012).
5.2 Transition phase (9.6–7.2 ka)
In line with existing research (Arthur and Dean, 1998; Bahr et al., 2006, 2008), the IMI occurred at ∼9.6 ka, leading to a significant change in colour (Fig. 2) and elemental composition of the sedimentary record (Fig. 3), along with a substantial increase in the abundance of HGs (Fig. 4). This increase does not coincide with higher TOC content, suggesting that enhanced preservation of HGs was not the cause. It is possible that these lipid biomarkers were transported fluvially to this site from lakes within the catchment basin of the Black Sea due to the warm/wet conditions at this time (Göktürk et al., 2011; Shumilovskikh et al., 2012; Filipova-Marinova et al., 2013). This, however, appears unlikely, as our site is located a substantial distance from the mouths of major rivers (>230 km) and the BIT index remains low during this period (∼0.08; Bingjie Yang, personal communications, 2024), indicating only a minor contribution of terrestrial organic matter at our site (Hopmans et al., 2004). Furthermore, as the subsequent period (7–5.6 ka) was also warm and wet (Göktürk et al., 2011; Shumilovskikh et al., 2012; Filipova-Marinova et al., 2013), we would expect the continuation of this peak if the HGs were being sourced from surrounding lacustrine environments. Instead, these high values decline abruptly after 6.6 ka.
It is therefore likely that the peak abundance in nitrogen-fixing cyanobacteria is related to warmer Black Sea surface temperatures during the early to mid-Holocene (Bahr et al., 2008) in combination with surface water stratification (Bahr et al., 2006). This stratification will have slowed the upward supply of fixed nitrogen, reducing nutrient availability in the surface waters, thereby promoting nitrogen fixation by cyanobacteria (Hutchins and Fu, 2017). This stratification may have been driven in part by enhanced freshwater influx due to wetter conditions but may also have been triggered by the IMI through the Bosporus Strait at ∼9.6 ka (Major et al., 2006; Bahr et al., 2008; Ankindinova et al., 2019). This IMI likely led to the gradual salinization of the water column over this transition interval and intermittent build-up of anoxia in the water column. This, in turn, led to periods of higher preservation of organic matter compared to the preceding period, as indicated by the slight increase in TOC after 9.6 ka. The presence of isorenieratene after 9.4 ka indicates that anoxia reached the photic zone at intermittent periods during this transition interval, thereby providing sufficient conditions for the presence of the anoxygenic photoautotrophs, Chlorobiaceae. While the peak in nitrogen-fixing cyanobacteria occurs ∼2 ka before anoxia intermittently entered the photic zone, the initial influx of dense saline water may have led to some reduction in vertical circulation, which reduced the amount of fixed nitrogen upwelled to the upper water column, leading to the presence of nitrogen-fixing cyanobacteria at 9.6 ka. This also coincides with a change in the distribution of HGs in our record between 9.7–6.9 ka, where hexose C28 diol and hexose C30 triol increase in abundance and hexose C28 triol declines in relative abundance and is no longer present after 9.1 ka, coinciding with the presence of isorenieratene. These changes may reflect a shift in species composition, linked to the gradual salinization and periodic anoxification of the water column after the IMI. The IMI at ∼9.6 ka appears, however, to have had little impact on the abundances of anammox and Thaumarchaeota. This is possibly because basin-wide water column stratification and the permanent build-up of anoxia did not occur until later in the record, meaning that neither process instantaneously reacted to the IMI at ∼9.6 ka.
5.3 Shift to anoxic brackish-to-marine mode of operation: a critical N-cycle threshold (∼7.2 ka to present)
After 7.2 ka, there was a substantial increase in TOC and TN and an abrupt shift in parts of the subsurface N cycle. The latter is shown by an increase in the BHT-x ratio, indicating an intensification of anammox, which is coeval with a decrease in crenarchaeol, indicating that there was a decline in Thaumarchaeota-driven nitrification. Studies have shown that, by ∼7.2 ka, anoxia had built up in the water column, as indicated by changes in redox elements (Fig. S2; Eckert et al., 2013; Wegwerth et al., 2018), and water column salinity had significantly increased (Fig. S3; Hiscott et al., 2007; Marret et al., 2009; Soulet et al., 2011; Filipova-Marinova et al., 2013), following the IMI from the Sea of Marmara at ∼9.6 ka (Major et al., 2002, 2006; Bahr et al., 2005, 2008; Ankindinova et al., 2019). This is supported by the presence of isorenieratene in our record during this time, which indicates that anoxia penetrated the photic zone. This water column anoxia likely led to the enhanced preservation of TOC and TN and triggered a shift in the subsurface N cycle, which crossed a threshold from an oxygenated lacustrine mode of operation to an anoxic brackish-to-marine mode of operation. The anoxic water column enabled anammox bacteria to expand their habitat from the anoxic sediments, where they were likely confined when the basin was an oxygenated freshwater environment, up into the suboxic/anoxic water column. This may therefore have commenced part of the modern-day N cycle in the Black Sea, where anammox activity occurs in the lower suboxic zone (∼100 m b.s.l.), where O2 is (nearly) depleted and H2S is absent (Jensen et al., 2008), with anammox bacteria consuming ammonium diffusing from the deep sea and utilizing the nitrite produced by both Thaumarchaeota and ammonia-oxidizing bacteria (AOB) (Kuypers et al., 2003; Lam et al., 2007). Consequently, it may be that the abundance of anammox bacteria increased as a result of the coupling to nitrite production by other microbes in the suboxic zone, whilst benefitting from ammonium diffusing upwards from the deep sea. The increased anammox after 7.1 ka likely indicates that more bioavailable nitrogen was lost from the system after the switch to the anoxic brackish-to-marine mode of operation. At the same time, Thaumarchaeota abundance declined, which may be in part due to the build-up of anoxia in the water column that reduced the niche of these aerobic microbes and the nitrification performed by them. Once these processes crossed a threshold from an oxygenated lacustrine mode of operation to an anoxic brackish-to-marine mode of operation, they appear to have remained steady for the remainder of the Holocene despite changes in the salinity of the basin (van der Meer et al., 2008; Mertens et al., 2012; Coolen et al., 2013) and significant changes in regional temperature and precipitation (Göktürk et al., 2011; Shumilovskikh et al., 2012; Filipova-Marinova et al., 2013). This shows that deoxygenation was the main driver of the observed change in anammox and in archaeal nitrification and that they were not affected by hydrological changes mainly occurring at the surface. While our record shows centennial-scale changes in the N cycle, we acknowledge that there may have also been rapid or short-term variations in N-cycle dynamics over this period that may not have been captured by the resolution of this record.
At 6.1 ka, the abundance of the HGs substantially declined, coinciding with an increase in δ15Nbulk, indicating a reduction in nitrogen fixation. This corresponds with the high-resolution record of Fulton et al. (2012), which also shows a decline in δ15Nbulk values at this time. As this decline in HG abundance and increase in δ15Nbulk does not coincide with a reduction in TOC, it is unlikely that reduced preservation of HGs played a role here. As nitrogen-fixing cyanobacteria inhabit the upper surface layer, it is likely that this change is linked to the salinization of the surface waters, with many studies demonstrating the disappearance of many freshwater mollusc, ostracod, and dinoflagellate cyst species at this time, which were replaced by an increased abundance of euryhaline Mediterranean species (Hiscott et al., 2007; Marret et al., 2009; Filipova-Marinova et al., 2013; Ivanova et al., 2015). At 6.1 ka, hexose C26 diol and hexose C28 diol are the only HGs present in the record, which may reflect the dominance of genera in the Nostocaceae family (i.e. Anabaena sp., Aphanizomenon sp., Nodularia sp., Nostoc sp.), as these members demonstrate a dominance of the hexose C26 diol and also contain varying amounts of hexose C28 diol (Gambacorta et al., 1999; Bauersachs et al., 2009). This distribution is similar to that of the Baltic Sea after ∼7.2 ka, when a series of weak intrusions of saline water led to the basin becoming fully brackish (Sollai et al., 2017). It is therefore possible that the peak in HGs in our Black Sea record between 9.6–6.9 ka represents a transition from the dominance of freshwater-tolerant nitrogen-fixing cyanobacteria to more brackish species, with brackish species dominating the surface waters after 6.6 ka. After 6.1 ka, δ15Nbulk gradually decreases, indicating a rise in nitrogen fixation, with this pattern also shown in previous records (Blumenberg et al., 2009; Fulton et al., 2012). While previous studies have shown riverine nitrogen to be a major source of fixed nitrogen in the modern Black Sea (McCarthy et al., 2007), due to the remoteness of our core site from the coast, our records are unlikely to have been significantly influenced by riverine input. It should be noted that a previous study (Fulton et al., 2012) suggested, based on compound specific measurements of pyropheophytin, that sedimentary δ15N in the Black Sea is primarily derived from eukaryotic algae rather than cyanobacteria, which exhibits a different fractionation of nitrogen isotopes. This means the use of δ15Nbulk as a nitrogen-fixation signal must be used with caution. HGs, however, are only derived from N-fixing cyanobacteria and are therefore an unambiguous biomarker of nitrogen fixation. Interestingly, at 4.3 ka, pentose HGs are detected, coinciding with lowest δ15Nbulk, indicating the presence of marine nitrogen-fixing cyanobacteria found in symbiosis with marine diatoms. This indicates that the surface water salinity had reached a threshold which enabled these marine microbes to survive, with research indicating salinity reached ∼17 ‰ during the deposition of Unit I (Ion et al., 2022), and freshwater/brackish species had disappeared by this time (Fig. S3; Filipova-Marinova et al., 2013). Indeed, reported increases in the number of euryhaline species at this time also point to the increasing salinity of the surface waters (Marret et al., 2009; Bradley et al., 2012), which may be linked to warmer/drier conditions that reduced freshwater influx and/or enhanced evaporation (Göktürk et al., 2011). Between 3.9–2.7 ka, isorenieratene is not detected in the samples, reflecting the findings of previous studies (Sinninghe Damsté et al., 1993). It has been suggested that this resulted from the erosion of the chemocline (Sinninghe Damsté et al., 1993), while other research shows a short reoccurrence of freshwater/brackish species (Fig. S3; Filipova-Marinova et al., 2013), which may indicate that enhanced freshwater input was responsible for lowering the chemocline below the photic zone. The disappearance of hexose HGs after 0.6 ka indicates that surface water salinities may more recently have become too high for the proliferation of brackish nitrogen-fixing cyanobacteria.
This study shows a relatively stable subsurface N cycle in the Black Sea over the Last Deglaciation and Holocene, with the exception of a critical threshold observed at 7.2 ka when the basin shifted from an oxygenated lacustrine environment to an anoxic brackish-to-marine basin. At this time, the loss of bioavailable nitrogen through anammox activity was enhanced and Thaumarchaeota-driven nitrification was reduced. Prior to and after this transition, the subsurface N cycle was remarkably stable despite various climatic and hydrological changes that impacted the basin during the deglaciation and the Holocene. Both the amount of nitrogen fixation by cyanobacteria and the composition of these microbes in the surface waters, however, appear to be much more dynamic and sensitive to hydrological changes over this period, responding in particular to salinity and temperature changes and stratification of the water column. Consequently, these records provide important insight into how future deoxygenation in marine environments may affect the microorganisms involved in the N cycle, possibly leading to enhanced loss of bioavailable nitrogen by anammox and reduced nitrification by Thaumarchaeota. Furthermore, in areas where enhanced water column stratification limits the supply of fixed nitrogen in the surface waters, localized enhanced cyanobacterial nitrogen fixation in the surface waters may occur. These changes may have associated feedbacks on nutrient cycling and carbon fixation, with implications for the future global carbon budget.
All data generated for this study are archived and publicly available in the Mendeley Data repository online at https://doi.org/10.17632/4c9fg7jf5d.1 (Cutmore et al., 2024).
The supplement related to this article is available online at https://doi.org/10.5194/cp-21-957-2025-supplement.
AC: conceptualization, formal analysis, investigation, data curation, visualization, writing (original draft and review and editing); NB: conceptualization, methodology, investigation, supervision, writing (review and editing); RH: resources, formal analysis, investigation, writing (review and editing); DR: formal analysis, writing (review and editing); BY: formal analysis, investigation, writing (review and editing); GJR: resources, supervision, writing (review and editing); ECH: supervision; SS: conceptualization, supervision, funding acquisition, writing (review and editing).
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
We thank the Chief Scientist Laura Villanueva and the captain and crew of the RV Pelagia for the collection of core 64PE418. We would like to thank Jaap Sinninghe Damsté for useful discussions. For laboratory support, we thank Anchelique Mets, Denise Dorhout, and Monique Verweij. Research cruise 64PE418 was funded by the SIAM Gravitation Grant (grant no. 024.002.002) from the Dutch Ministry of Education, Culture and Science (OCW). This study was funded by the Netherlands Earth System Science Centre (grant no. 024.002.001) from the Dutch Ministry of Education, Culture and Science (OCW).
This research has been supported by the Netherlands Earth System Science Centre (grant no. 024.002.001).
This paper was edited by Erin McClymont and reviewed by two anonymous referees.
Aksu, A., Hiscott, R. N., Kaminski, M. A., Mudie, P. J., Gillespie, H., Abrajano, T., and Yasar, D.: Last glacial–Holocene paleoceanography of the Black Sea and Marmara Sea: stable isotopic, foraminiferal and coccolith evidence, Mar. Geol., 190, 119–149, https://doi.org/10.1016/S0025-3227(02)00345-6, 2002.
Algeo, T. J. and Li, C.: Redox classification and calibration of redox thresholds in sedimentary systems, Geochim. Cosmochim. Ac., 287, 8–26, https://doi.org/10.1016/j.gca.2020.01.055, 2020.
Ankindinova, O., Hiscott, R. N., Aksu, A. E., and Grimes, V.: High-resolution Sr-isotopic evolution of Black Sea water during the Holocene: Implications for reconnection with the global ocean, Mar. Geol., 407, 213–228, https://doi.org/10.1016/j.margeo.2018.11.004, 2019.
Arthur, M. A. and Dean, W. E.: Organic-matter production and preservation and evolution of anoxia in the Holocene Black Sea, Paleoceanogr. Paleoclimatol., 13, 395–411, https://doi.org/10.1029/98PA01161, 1998.
Badertscher, S., Fleitmann, D., Cheng, H., Edwards, R. L., Göktürk, O. M., Zumbühl, A., Leuenberger, M., and Tüysüz, O.: Pleistocene water intrusions from the Mediterranean and Caspian seas into the Black Sea, Nat. Geosci., 4, 236–239, https://doi.org/10.1038/ngeo1106, 2011.
Bale, N., Hopmans, E. C., Zell, C., Sobrinho, R. L., Kim, J.-H., Sinninghe Damsté, J. S., Villareal, T. A., and Schouten, S.: Long chain glycolipids with pentose head groups as biomarkers for marine endosymbiotic heterocystous cyanobacteria, Org. Geochem., 81, 1–7, https://doi.org/10.1016/j.orggeochem.2015.01.004, 2015.
Bale, N., Ding, S., Hopmans, E. C., Arts, M. G. I., Villanueva, L., Boschman, C., Haas, A. F., Schouten, S., and Sinninghe Damsté, J. S.: Lipidomics of Environmental Microbial Communities. I: Visualization of Component Distributions Using Untargeted Analysis of High-Resolution Mass Spectrometry Data, Front. Microbiol., 12, 1–15, https://doi.org/10.3389/fmicb.2021.659302, 2021.
Bale, N. J., Hennekam, R., Hopmans, E. C., Dorhout, D., Reichart, G.-J., van der Meer, M. T. J., Villareal, T. A., Sinninghe Damsté, J. S., and Schouten, S.: Biomarker evidence for nitrogen-fixing cyanobacterial blooms in a brackish surface layer in the Nile River plume during sapropel deposition, Geology, 47, 1088–1092, https://doi.org/10.1130/G46682.1, 2019.
Bahr, A., Lamy, F., Arz, H., Kuhlmann, H., and Wefer, G.: Late glacial to Holocene climate and sedimentation history in the NW Black Sea, Mar. Geol., 214, 309–322, https://doi.org/10.1016/j.margeo.2004.11.013, 2005.
Bahr, A., Arz, H., Lamy, F., and Wefer, G.: Late glacial to Holocene paleoenvironmental evolution of the Black Sea, reconstructed with stable oxygen isotope records obtained on ostracod shells, Earth Planet. Sc. Lett., 241, 863–875, https://doi.org/10.1016/j.epsl.2005.10.036, 2006.
Bahr, A., Lamy, F., Arz, H., Major, C., Kwiecien, O., and Wefer, G.: Abrupt changes of temperature and water chemistry in the late Pleistocene and early Holocene Black Sea, Geochem. Geophy. Geosy., 9, 1–16, https://doi.org/10.1029/2007GC001683, 2008.
Bauersachs, T., Compaoré, J., Hopmans, E. C., Stal, L. J., Schouten, S., and Sinninghe Damsté, J.: Distribution of heterocyst glycolipids in cyanobacteria, Phytochemistry, 70, 2034–2039, https://doi.org/10.1016/j.phytochem.2009.08.014, 2009.
Bauersachs, T., Speelman, E. N., Hopmans, E. C., Reichart, G.-J., Schouten, S., and Sinninghe Damsté, J.: Fossilized glycolipids reveal past oceanic N2 fixation by heterocystous cyanobacteria, Earth Planet. Sc. Lett., 107, 19190–19194, https://doi.org/10.1073/pnas.1007526107, 2010.
Blaauw, M.: Methods and code for “classical” age-modelling of radiocarbon sequences, Quat. Geochronol., 5, 512–518, https://doi.org/10.1016/j.quageo.2010.01.002, 2010.
Blumenberg, M., Seifert, R., Kasten, S., Bahlmann, E., and Michaelis, W.: Euphotic zone bacterioplankton sources major sedimentary bacteriohopanepolyols in the Holocene Black Sea, Geochim. Cosmochim. Ac., 73, 750–766, https://doi.org/10.1016/j.gca.2008.11.005, 2009.
Bopp, L., Resplandy, L., Orr, J. C., Doney, S. C., Dunne, J. P., Gehlen, M., Halloran, P., Heinze, C., Ilyina, T., Séférian, R., Tjiputra, J., and Vichi, M.: Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models, Biogeosciences, 10, 6225–6245, https://doi.org/10.5194/bg-10-6225-2013, 2013.
Bradley, L. R., Marret, F., Mudie, P. J., Aksu, A. E., and Hiscott, R. N.: Constraining Holocene sea-surface conditions in the south-western Black Sea using dinoflagellate cysts, J. Quaternary Sci., 27, 835–843, https://doi.org/10.1002/jqs.2580, 2012.
Coolen, M. J. L., Orsib, W. D., Balkema, C., Quince, C., Harris, K., Sylva, S. P., Filipova-Marinova, M., and Giosan, L.: Evolution of the plankton paleome in the BlackSea from the Deglacial to Anthropocene, P. Natl. Acad. Sci. USA, 110, 8609–8614, https://doi.org/10.1073/pnas.1219283110, 2013.
Cutmore, A., Bale, N., Hennekam, R., Yang, B., Rush, D., Reichart, G.-J., Hopmans, E. C., and Schouten, S.: Lipid biomarker records from Black Sea core 64PE418 spanning the Last Deglaciation and Holocene, Mendeley Data, V1 [data set], https://doi.org/10.17632/4c9fg7jf5d.1, 2024.
Czernik, J. and Goslar, T.: Preparation of Graphite Targets in the Gliwice Radiocarbon Laboratory for AMS 14C Dating, Radiocarbon, 43, 283–291, https://doi.org/10.1017/S0033822200038121, 2001.
Dalsgaard, T., Thamdrup, B., Farías, L., and Revsbech, N. P.: Anammox and denitrification in the oxygen minimum zone of the eastern South Pacific, Limnol. Oceanogr., 57, 1331–1346, https://doi.org/10.4319/lo.2012.57.5.1331, 2012.
Eckert, S., Brumsack, H.-J., Severmann, S., Schnetger, B., März, C., and Fröllje, H.: Establishment of euxinic conditions in the Holocene Black, Geology, 41, 431–434, https://doi.org/10.1130/G33826.1, 2013.
Elling, F. J., Hemingway, J. D., Kharbush, J. J., Becker, K. W., Polik, C. A., and Pearson, A.: Linking diatom-diazotroph symbioses to nitrogen cycle perturbations and deep-water anoxia: Insights from Mediterranean sapropel events, Earth Planet. Sc. Lett., 571, 1–11, https://doi.org/10.1016/j.epsl.2021.117110, 2021.
Falkowski, P. G., Barber, R. T., and Smetacek, V.: Biogeochemical controls and feedbacks on ocean primary production, Science, 281, 200–206, https://doi.org/10.1126/science.281.5374.200, 1998.
Filipova-Marinova, M., Pavlov, D., Coolen, M., and Giosan, L.: First high-resolution marinopalynological stratigraphy of Late Quaternary sediments from the central part of the Bulgarian Black Sea area, Quatern. Int., 293, 170–183, https://doi.org/10.1016/j.quaint.2012.05.002, 2013.
Francis, C. A., Roberts, K. J., Beman, J. M., Santoro, A. E., and Oakley, B. B.: Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean, P. Natl. Acad. Sci. USA, 102, 14683–14688, https://doi.org/10.1073/pnas.0506625102, 2005.
Fulton, J. M., Arthur, M. A., and Freeman, K. H.: Black Sea nitrogen cycling and the preservation of phytoplankton δ15N signals during the Holocene, Global Biogeochem. Cy., 26, 1–15, https://doi.org/10.1029/2011GB004196, 2012.
Gambacorta, A., Trincone, A., Soriente, A., and Sodano, G.: Chemistry of glycolipids from the heterocysts of nitrogen-fixing cyanobacteria, Phytochemistry, 2, 145–150, 1999.
Göktürk, O. M., Fleitmann, D., Badertscher, S., Cheng, H., Edwards, R. L., Leuenberger, M., Fankhauser, A., Tüysüz, O., and Kramers, J.: Climate on the southern Black Sea coast during the Holocene: implications from the Sofular Cave record, Quaternary Sci. Rev., 30, 2433–2445, https://doi.org/10.1016/j.quascirev.2011.05.007, 2011.
Heaton, T. J., Köhler, P., Butzin, M., Bard, E., Reimer, R. W., Austin, W. E. N., Bronk Ramsey, C., Grootes, P. M., Hughen, K. A., Kromer, B., Reimer, P. J., Adkins, J., Burke, A., Cook, M. S., Olsen, J., and Skinner, L. C.: Marine20 – The Marine Radiocarbon Age Calibration Curve (0–55 000 cal BP), Radiocarbon, 62, 779–820, https://doi.org/10.1017/RDC.2020.68, 2020.
Hennekam, R., van der Bolt, B., van Nes, E. H., de Lange, G.-J., Scheffer, M., and Reichart, G.-J.: Early-warning signals for marine anoxic events, Geophys. Res. Lett., 47, 1–9, https://doi.org/10.1029/2020GL089183, 2020.
Hiscott, R. N., Aksu, A. E., Mudie, P. J., Marret, F., Abrajano, T., Kaminski, M. A., Evans, J., Çakiroğlu, A.İ., and Yaşar, D.: A gradual drowning of the southwestern Black Sea shelf: Evidence for a progressive rather than abrupt Holocene reconnection with the eastern Mediterranean Sea through the Marmara Sea Gateway, Quatern. Int., 167–168, 19–34, https://doi.org/10.1016/j.quaint.2006.11.007, 2007.
Hopmans, E. C., Weijers, J. W. H., Schefuß, E., Herfort, L., Sinninghe Damsté, J. S., and Schouten, S.: A novel proxy for terrestrial organic matter in sediments based on branched and isoprenoid tetraether lipids, Earth Planet. Sc. Lett., 224, 107–116, https://doi.org/10.1016/j.epsl.2004.05.012, 2004.
Hopmans, E. C., Smit, N. T., Schwartz-Narbonne, R., Sinninghe Damste, J. S., and Rush, D.: Analysis of non-derivatized bacteriohopanepolyols using UHPLC-HRMS reveals great structural diversity in environmental lipid assemblages, Org. Geochem., 160, 1–17, https://doi.org/10.1016/j.orggeochem.2021.104285, 2021.
Huang, Y., Zheng, Y., Heng, P., Giosan, L., and Coolen, M. J. L.: Black Sea paleosalinity evolution since the last deglaciation reconstructed from alkenone-inferred Isochrysidales diversity, Earth Planet. Sc. Lett., 564, 1–9, https://doi.org/10.1016/j.epsl.2021.116881, 2021.
Hutchins, D. A. and Fu, F.: Microorganisms and ocean global change, Nat. Microbiol., 2, 17058, https://doi.org/10.1038/nmicrobiol.2017.58, 2017.
Ion, G., Briceag, A., Vasiliu, D., Lupaşcu, N., and Melinte-Dobrinescu, M.: A multiproxy reconstruction of the Late Pleistocene-Holocene paleoenvironment: New insights from the NW Black Sea, Mar. Geol., 443, 1–19, https://doi.org/10.1016/j.margeo.2021.106648, 2022.
Ivanova, E. V., Murdmaa, I. O., Chepalyga, A. L., Cronin, T. M., Pasechnik, I. V., Levchenko, O. V., Howe, S. S., Manushkina, A. V., and Platonova, E. A.: Holocene sea-level oscillations and environmental changes on the Eastern Black Sea shelf, Palaeogeogr. Palaeocl., 246, 228–259, https://doi.org/10.1016/j.palaeo.2006.09.014, 2007.
Ivanova, E. V., Marret, F., Zenina, M. A., Murdmaa, I. O., Chepalyga, A. L., Bradley, L. R., Schornikov, E. I., Levchenko, O. V., and Zyryanova, M. I.: The Holocene Black Sea reconnection to the Mediterranean Sea: New insights from the northeastern Caucasian shelf, Palaeogeogr. Palaeocl., 427, 41–61, https://doi.org/10.1016/j.palaeo.2015.03.027, 2015.
Jaeschke, A., Hopmans, E. C., Wakeham, S. G., Schouten, S., and Sinninghe Damsté, J. S.: The presence of ladderane lipids in the oxygen minimum zone of the Arabian Sea indicates nitrogen loss through anammox, Limnol. Oceanogr., 52, 780–786, https://doi.org/10.4319/lo.2007.52.2.0780, 2007.
Jensen, M. M., Kuypers, M. M. M., Lavik, G., and Thamdrup, B.: Rates and regulation of anaerobic ammonium oxidation and denitrification in the Black Sea, Limnol. Oceanogr., 53, 23–36, https://doi.org/10.4319/lo.2008.53.1.0023, 2008.
Jones, G. A. and Gagnon, A. R.: Radiocarbon chronology of Black Sea sediments, Deep-Sea Res. Pt. I, 41, 531–557, https://doi.org/10.1016/0967-0637(94)90094-9, 1994.
Kalvelage, T., Lavik, G., Lam, P., Contreras, S., Arteaga, L., Löscher, C. R., Oschlies, A., Paulmier, A., Stramma, L., and Kuypers, M. M. M.: Nitrogen cycling driven by organic matter export in the South Pacific oxygen minimum zone, Nat. Geosci., 6, 228–234, https://doi.org/10.1038/ngeo1739, 2013.
Karner, M. B., DeLong, E. F., and Karl, D. M.: Archaeal dominance in the mesopelagic zone of the Pacific Ocean, Nature, 409, 507–510, https://doi.org/10.1038/35054051, 2001.
Keeling, R. F., Körtzinger, A., and Gruber, N.: Ocean Deoxygenation in a Warming World, Annu. Rev. Mar. Sci., 2, 199–229, https://doi.org/10.1146/annurev.marine.010908.163855, 2010.
Knowles, R.: Denitrification, Microbiol. Rev., 46, 43–70, https://doi.org/10.1128/mr.46.1.43-70.1982, 1982.
Könneke, M., Bernhard, A. E., de la Torre, J. R., Walker, C. B., Waterbury, J. B., and Stahl, D. A.: Isolation of an autotrophic ammonia-oxidizing marine archaeon, Nature, 437, 543–546, https://doi.org/10.1038/nature03911, 2005.
Koopmans, M. P., Köster, J., Van Kaam-Peters, H. M. E., Kenig, F., Schouten, S., Hartgers, W. A., de Leeuw, J. W., and Sinninghe Damsté, J. S.: Diagenetic and catagenetic products of isorenieratene: Molecular indicators for photic zone anoxia, Geochim. Cosmochim. Ac., 60, 4467–4496, https://doi.org/10.1016/S0146-6380(97)00025-9, 1996.
Kuenen, J. G. and Robertson, L. A.: Ecology of nitrification and denitrification, in: The Nitrogen and Sulphur Cycles, edited by: Cole, J. A. and Ferguson, S. J., Cambridge University Press, Cambridge, 161–218, ISBN 0521351995, 1988.
Kusch, S., Wakeham, S. G., and Sepúlveda, J.: Bacteriohopanepolyols across the Black Sea redoxcline trace diverse bacterial metabolisms, Org. Geochem., 172, 1–18, https://doi.org/10.1016/j.orggeochem.2022.104462, 2022.
Kuypers, M. M. M., Sliekers, A. O., Lavik, G., Schmid, M., Barker Jørgensen, B., Kuenen, J. G., Sinninghe Damsté, J. S., Strous, M., and Jetten, M. S. M.: Anaerobic ammonium oxidation by anammox bacteria in the Black Sea, Nature, 422, 608–611, https://doi.org/10.1038/nature01472, 2003.
Kwiecien, O., Arz, H. W., Lamy, F., Wulf, S., Bahr, A., Röhl, U., and Haug, G. H.: Estimated Reservoir Ages of the Black Sea Since the Last Glacial, Radiocarbon, 50, 99–118, https://doi.org/10.1017/S0033822200043393, 2008.
Lam, P., Jensen, M. M., Lavik, G., McGinnis, D. F., Müller, B., Schubert, C. J., Amann, R., Thamdrup, B., and Kuypers, M. M. M.: Linking crenarchaeal and bacterial nitrification to anammox in the Black Sea, P. Natl. Acad. Sci. USA, 104, 7104–7109, https://doi.org/10.1073/pnas.0611081104, 2007.
Limaye, R. B., Padmalal, D., and Kumaran, K. P. N.: Cyanobacteria and testate amoeba as potential proxies for Holocene hydrological changes and climate variability: Evidence from tropical coastal lowlands of SW India, Quatern. Int., 443, 99–114, https://doi.org/10.1016/j.quaint.2016.09.044, 2017.
Major, C., Ryan, W., Lericolais, G., and Hajdas, I.: Constraints on Black Sea outflow to the Sea of Marmara during the last glacial–interglacial transition, Mar. Geol., 190, 19–34, https://doi.org/10.1016/S0025-3227(02)00340-7, 2002.
Major, C., Goldstein, S., Ryan, W., Lericolais, G., Piotrowski, A. M., and Hajdas, I.: The co-evolution of Black Sea level and composition through the last deglaciation and its paleoclimatic significance, Quaternary Sci. Rev., 25, 2031–2047, https://doi.org/10.1016/j.quascirev.2006.01.032, 2006.
Marret, F., Mudie, P., Aksu, A., and Hiscott, R. N.: A Holocene dinocyst record of a two-step transformation of the Neoeuxinian brackish water lake into the Black Sea, Quatern. Int., 197, 72–86, https://doi.org/10.1016/j.quaint.2007.01.010, 2009.
McCarthy, J. J., Yilmaz, A., Coban-Yildiz, Y., and Nevins, J. L.: Nitrogen cycling in the offshore waters of the Black Sea, Estuar. Coast. Shelf S., 74, 493–514, https://doi.org/10.1016/j.ecss.2007.05.005, 2007.
Ménot, G. and Bard, E.: A precise search for drastic temperature shifts of the past 40 000 years in southeastern Europe, Palaeoceanography, 27, 1–13, https://doi.org/10.1029/2012PA002291, 2012.
Mertens, K. N., Bradley, L. R., Takano, Y., Mudie, P. R., Marret, F., Aksu, A. E., Hiscott, R. N., Verleye, T. J., Mousing, E. A., Smyrnova, L. L., Bagheri, S., Mansor, M., Pospelova, V., and Matsuoka, K.: Quantitative estimation of Holocene surface salinity variation in the Black Sea using dinoflagellate cyst process length, Quaternary Sci. Rev., 39, 45–59, https://doi.org/10.1016/j.quascirev.2012.01.026, 2012.
Murray, J. W., Jannasch, H. W., Honjo, S., Anderson, R. F., Reeburgh, W. S., Top, Z., Friederich, G. E., Codispoti, L. A., and Izdar, E.: Unexpected changes in the oxic/anoxic interface in the Black Sea, Nature, 337, 411–413, https://doi.org/10.1038/338411a0, 1989.
Murray, J. W., Codispoti, L. A., and Friederich, G. E.: Oxidation-reduction environments: the suboxic zone in the Black Sea, in: Aquatic Chemistry: Interfacial and Interspecies Processes, edited by: Huang, C. P., O'Melia, C. R., and Morgan, J. J., ACS Advances in Chemistry, vol. 224, ACS Publications, Washington, DC, USA, 157–176, https://doi.org/10.1021/ba-1995-0244.ch007, 1995.
Naafs, B. D. A., Monteiro, F. M., Pearson, A., Higgins, M. B., Pancost, R. D., and Ridgwell, A.: Fundamentally different global marine nitrogen cycling in response to severe ocean deoxygenation, P. Natl. Acad. Sci. USA, 116, 24979–24984, https://doi.org/10.1073/pnas.1905553116, 2019.
Nicholas, W. A., Chivas, A. R., Murray-Wallace, C. V., and Fink, D.: Prompt transgression and gradual salinisation of the Black Sea during the early Holocene constrained by amino acid racemization and radiocarbon dating, Quaternary Sci. Rev., 30, 3769–3790, https://doi.org/10.1016/j.quascirev.2011.09.018, 2011.
Özsoy, E. and Ünlüata, Ü.: Oceanography of the Black Sea: A review of some recent results, Earth-Sci. Rev., 42, 231–272, https://doi.org/10.1016/S0012-8252(97)81859-4, 1997.
Pérez Gallego, R., von Meijenfeldt, B., Bale, N. J., Sinninghe Damsté, J. S., and Villanueva, L.: Emergence and evolution of heterocyte glycolipid biosynthesis enabled specialized nitrogen fixation in cyanobacteria, Microbiology, 122, e2413972122, https://doi.org/10.1073/pnas.2413972122, 2025.
Piper, D. Z. and Calvert, S. E.: Holocene and late glacial palaeoceanography and palaeolimnology of the Black Sea: Changing sediment provenance and basin hydrography over the past 20 000 years, Geochim. Cosmochim. Ac., 75, 5597–5624, https://doi.org/10.1016/j.gca.2011.07.016, 2011.
Ploug, H.: Cyanobacterial surface blooms formed by Aphanizomenon sp. and Nodularia spumigena in the Baltic Sea: Small-scale fluxes, pH, and oxygen microenvironments, Limnol. Oceanogr., 53, 914–921, https://doi.org/10.4319/lo.2008.53.3.0914, 2008.
Reimer, P. J., Austin, W. E. N., Bard, E., Bayliss, A., Blackwell, P. G., Bronk Ramsey, C., Butzin, M., Cheng, H., Edwards, R. L., Friedrich, N., Grootes, P. M., Guilderson, T. P., Hajdas, I., Heaton, T. J., Hogg, A. G., Hughen, K. A., Kromer, B., Manning, S. W., Muscheler, R., Palmer, J. G., Pearson, C., van der Plicht, J., Reimer, R. W., Richards, D. A., Scott, E. M., Southon, J. R., Turney, C. S. M., Wacker, L., Adolphi, F., Büntgen, U., Capano, M., Fahrni, S. M., Fogtmann-Schulz, A., Friedrich, R., Köhler, P., Kudsk, S., Miyake, F., Olsen, J., Reinig, F., Sakamoto, M., Sookdeo, A., and Talamo, S.: The IntCal20 Northern Hemisphere Radiocarbon Age Calibration Curve (0–55 cal kBP), Radiocarbon, 62, 725–757, https://doi.org/10.1017/RDC.2020.41, 2020.
Rush, D. and Sinninghe Damsté, J. S.: Lipids as paleomarkers to constrain the marine nitrogen cycle, Environ. Microbiol., 19, 2119–2132, https://doi.org/10.1111/1462-2920.13682, 2017.
Rush, D., Sinninghe Damsté, J. S., Poulton, S. W., Thamdrup, B., Garside, A. L., Acuña González, J., Schouten, S., Jetten, M. S. M., and Talbot, H. M.: Anaerobic ammonium-oxidising bacteria: A biological source of the bacteriohopanetetrol stereoisomer in marine sediments, Geochim. Cosmochim. Ac., 140, 50–64, https://doi.org/10.1016/j.gca.2014.05.014, 2014.
Schouten, S., Hopmans, E. C., and Sinninghe Damsté, J. S.: The organic geochemistry of glycerol dialkyl glycerol tetraether lipids: A review, Org. Geochem., 54, 19–61, https://doi.org/10.1016/j.orggeochem.2012.09.006, 2013.
Schrader, H.-J.: Quaternary Paleoclimatology of the Black Sea basin, Sediment. Geol., 23, 165–180, https://doi.org/10.1016/0037-0738(79)90013-7, 1979.
Schwartz-Narbonne, N., Schaeffer, P., Hopmans, E. C., Schenesse, M., Charlton, E. A., Jones, D. M., Sinninghe Damsté, J. S., Farhan, M., Haque, U., Jetten, M. S. M., Lengger, S. K., Murrell, J. C., Normand, P., Nuijten, G. H. L., Talbot, H. M., and Rush, D.: A unique bacteriohopanetetrol stereoisomer of marine anammox, Org. Geochem., 143, 1–10, https://doi.org/10.1016/j.orggeochem.2020.103994, 2020.
Shumilovskikh, L. S., Tarasov, P., Arz, H. W., Fleitmann, D., Marret, F., Nowaczyk, N., Plessen, B., Schlütz, F., and Behling, H.: Vegetation and environmental dynamics in the southern Black Sea region since 18 kyr BP derived from the marine core 22-GC3, Palaeogeogr. Palaeocl., 337–338, 177–193, https://doi.org/10.1016/j.palaeo.2012.04.015, 2012.
Sinninghe Damste, J. S., Wakeham, S. G., Kohnen, M. E. L., Hayes, J. M., and de Leeuw, J. W.: A 6000 year sedimentary molecular record of chemocline excursions in the Black Sea, Nature, 362, 827–829, https://doi.org/10.1038/362827a0, 1993.
Sinninghe Damsté, J. S., Schouten, S., Hopmans, E. C., van Duin, A. C. T., and Geenevasen, A. J. A.: Crenarchaeol, J. Lipid Res., 43, 1641–1651, https://doi.org/10.1194/jlr.m200148-jlr200, 2002.
Sollai, M., Hopmans, E. C., Bale, N. J., Mets, A., Warden, L., Moros, M., and Sinninghe Damsté, J. S.: The Holocene sedimentary record of cyanobacterial glycolipids in the Baltic Sea: an evaluation of their application as tracers of past nitrogen fixation, Biogeosciences, 14, 5789–5804, https://doi.org/10.5194/bg-14-5789-2017, 2017.
Sollai, M., Villanueva, L., Hopmans, E. C., Reichart, G.-J., and Sinninghe Damsté, J. S.: A combined lipidomic and 16S rRNA gene amplicon sequencing approach reveals archaeal sources of intact polar lipids in the stratified Black Sea water column, Geobiology, 17, 91–109, https://doi.org/10.1111/gbi.12316, 2018.
Soulet, G., Ménot, G., Lericolais, G., and Bard, E.: A revised calendar age for the last reconnection of the Black Sea to the global ocean, Quaternary Sci. Rev., 30, 1019–1026, https://doi.org/10.1016/j.quascirev.2011.03.001, 2011.
Stuiver, M. and Polach, H. A.: Discussion Reporting of 14C Data, Radiocarbon, 19, 355–363, https://doi.org/10.1017/S0033822200003672, 1977.
van de Graaf, A. A., de Bruijn, P., Robertson, L. A., Jetten, M. S. M., and Kuenen, J. G.: Metabolic pathway of anaerobic ammonium oxidation on the basis of 15N studies in a fluidized bed reactor, Microbiology, 143, 2415–2421, https://doi.org/10.1099/00221287-143-7-2415, 1997.
van der Meer, M. J. M., Sangiorgi, F., Baas, M., Brinkhuis, H., Sinninghe Damsté, J. S., and Schouten, S.: Molecular isotopic and dinoflagellate evidence for Late Holocene freshening of the Black Sea, Earth Planet. Sc. Lett., 267, 426–434, https://doi.org/10.1016/j.epsl.2007.12.001, 2008.
van Kemenade, Z. R., Cutmore, A., Hennekam, R., Hopmans, E. C., van der Meer, M. T. J., Mojtahid, M., Jorissen, F. J., Bale, N. J., Reichart, G.-J., Sinninghe Damsté, J. S., and Rush, D.: Marine nitrogen cycling dynamics under altering redox conditions: insights from deposition of sapropels S1 and the ambiguous S2 in the Eastern Mediterranean Sea, Geochim. Cosmochim. Ac., 354, 197–210, https://doi.org/10.1016/j.gca.2023.06.018, 2023.
Verleye, T. J., Mertens, K. N., Louwye, S., and Arz, H. W.: Holocene salinity changes in the southwestern black sea: A reconstruction based on dinoflagellate cysts, Palynology, 33, 77–100, https://doi.org/10.1080/01916122.2009.9989666, 2009.
Villareal, T. A.: Marine Nitrogen-Fixing Diatom-Cyanobacteria Symbioses, in: Marine Pelagic Cyanobacteria: Trichodesmium and other Diazotrophs, edited by: Carpenter, E. J., Capone, D. G., and Rueter, J. G., Springer, Dordrecht, 163–175, https://doi.org/10.1007/978-94-015-7977-3_10, 1992.
Wegwerth, A., Eckert, S., Dellwig, O., Schnetger, B., Severmann, S., Weyer, S., Brüske, A., Kaiser, J., Köster, J., Arz, H. W., and Brumsack, H.-J.: Redox evolution during Eemian and Holocene sapropel formation in the Black Sea, Palaeogeogr. Palaeocl., 489, 249–260, https://doi.org/10.1016/j.palaeo.2017.10.014, 2018.
Whitton, B. A. and Mateo, P.: Rivulariaceae, in: Ecology of Cyanobacteria II, edited by: Whitton, B. A., Springer, Dordrecht, 561–591, https://doi.org/10.1007/978-94-007-3855-3_22, 2012.
Wörmer, L., Cires, S., Velazquez, D., Quesada, A., and Hinrichs, K.-U.: Cyanobacterial heterocyst glycolipids in cultures and environmental samples: Diversity and biomarker potential, Limnol. Oceanogr., 57, 1775–1788, https://doi.org/10.4319/lo.2012.57.6.1775, 2012.
Wuchter, C., Abbas, B., Coolen, M. J. L., Herfort, L., van Bleijswijk, J., Timmers, P., Strous, M., Teira, E., Herndl, G. J., Middelburg, J. J., Schouten, S., and Sinninghe Damsté, J. S.: Archaeal nitrification in the ocean, P. Natl. Acad. Sci. USA, 103, 12317–12322, https://doi.org/10.1073/pnas.0600756103, 2006.
Yanchilina, A. G., Ryan, W. B. F., McManus, J., Dimitrov, P., Dimitrov, D., Slavova, K., and Filipova-Marinova, M.: Compilation of geophysical, geochronological, and geochemical evidence indicates a rapid Mediterranean-derived submergence of the Black Sea's shelf and subsequent substantial salinification in the early Holocene, Mar. Geol., 383, 14–34, https://doi.org/10.1016/j.margeo.2016.11.001, 2017.