Abstract
Introduction/setting
The Last Glacial Maximum (LGM) represents the last major global glaciation, which caused a drop of sea level by ~130 m when global ice volumes reached their maximum (e.g. Lambeck et al., 2014). The LGM lasted from 27.5 to 23.3 ka BP (Hughes and Gibbard, 2015; Rasmussen et al., 2014) or according to Clark et al. (2009) between 26.5 and 19 ka BP. In the area of the modern central and southern North Sea (cf. Figure 1) this lowered sea level resulted in an emerged landscape, “Doggerland” (e.g. Coles, 1998), which connected the actual British Isles with the European mainland.

Map of the three sampling sites A, B, and C in the central North Sea at the Dogger Bank investigated in January 2021 and their position relative to that of the salt dome “Britta.” Data for mapped area with salt-dome overlying shallow gas (gas indications <1000 m sediment depth) were taken from Müller et al. (2018) and for the salt dome from Arfai et al. (2014). The seafloor depressions were previously published in Römer et al. (2021a). Also shown is an abandoned well. The image taken with a camera attached to the BGR’s “Golden Eye” electromagnetic profiler gives an impression on how the patchy peats crop out at the sampling sites. The inlet shows the borders of the Exclusive Economic Zones (EEZ) in the central North Sea including the German EEZ and the area in the Northwest called “Entenschnabel” and position of the area investigated in this study (red star). The green lines denote the approximate coastlines at 10,000 and 9,000 a BP (from Walker et al., 2020).
With the subsequent main phase of deglaciation after the LGM bracketed between 16.5 and 8.2 ka BP according to Lambeck et al. (2014), global sea level rose continuously. The interplay of these processes lead to a global increase in atmospheric methane from the escape from boreal wetlands (Fischer et al., 2008; Loulergue et al., 2008). The landscape of Doggerland that had since developed from a permafrost dominated glacial foreland into a vegetated area structured by older moraine sediments forming hills, scattered lakes, a river drainage system, and wide peat environments especially along the coastlines, slowly submerged again. The sea approached from both the North and the English Channel in the Southwest causing slowly and continuously the development of bays, peninsulas and islands in the topographic highest areas (“Dogger Island”), today’s “Dogger Bank,” in the central North Sea. A more recent reconstruction of the drowning landscape was given by Sturt et al. (2013).
Since the pioneering work of Jelgersma (1961), the reconstruction of sea level during the Holocene in the North Sea area has been further developed, but is still based on precise dating of so-called basal peats. Formation of these peats started with the rising sea level and subsequently with the rising ground water along the slowly submerging coastlines. Coastal peat formed extensively around local mean high-water level - in the landward-shifting zone (e.g. Hijma and Cohen, 2010; Krüger et al., 2017). Precise dating of peats by the radiocarbon method allowed determination of sea level index points (SIP’s) providing a framework of where and when relative sea level was located. Recent reconstructions of relative sea level development are from Vink et al. (2007) for the German Bight and from Hijma and Cohen (2010) for the Netherlands and the southern North Sea. For some sea level reconstructions glacial-isostatic-adjustment models were applied to compensate for effects of isostatic glacial rebound of the foreland from the ice masses during the LGM (e.g. Sturt et al., 2013; Vink et al., 2007).
Apart from their use as SIP’s, basal peat deposits record also the development of the post glacial environment in the form of pollen for example. In addition, compositions of biomarker in peats and compound-specific (stable) carbon isotope signatures are excellently suited to characterize their depositional environment (e.g. Pancost et al., 2000). While in the geological transformation from organic matter toward coals for example, much of the biomarker information is blurred by chemical and physical modification reactions, the relatively immature peats still contain the majority of primary signals. Combined biomarker studies have therefore contributed crucially to a better understanding of peats (and peat-derived lignites) from various geological periods in the Cenozoic, including those covering important environmental perturbations like the Paleocene-Eocene Thermal Maximum (PETM; Elvert et al., 2016; Pancost et al., 2002, 2007; Petersen et al., 2009; Schouten et al., 2007; Venkatesan et al., 1986; Weijers et al., 2006; Xie et al., 2000). Beside the higher plant organic matter that is recorded mostly in mid- and long-chained odd-even predominated
Methanotrophic bacteria appear to play a crucial role for methane turnover in peats and particularly during geological times like the PETM. Intense terrestrial methane cycling is demonstrated by relative 13C-depletions in specific hopanoids (Pancost et al., 2007) and such depletions were highest during the PETM and the Holocene after the LGM (Inglis et al., 2019). Analysis of pollen contained in the so-called basal peats from “Doggerland” demonstrated their origin from variable plant communities typical for boreal settings (Krüger et al., 2017). This confirms previous geochemical and microfacies studies of related coastal lowlands of Northern Germany particularly seeking a better understanding of floral changes (Gerdes et al., 2003; Köller, 2002). Relatively minor studied is, however, the role of peats for methane generation and other microbial processes in the recent North Sea sediments. One recent publication determined the microbial methane turnover and suggested a relatively minor role of peat-associated hydrocarbons for methane seepage in the area (Lippmann et al., 2021).
In this study, five peat samples from three sites were taken from mini-multicorer and gravity core sampling in the German sector of the central North Sea. At two of the sites peat samples from two depths (5 and 10 cm and 10 and 20 cm below seafloor (cmbsf), respectively) were taken to obtain geochemical insights into fluctuations during peat-forming and other processes recorded in these samples. Also sampled and analyzed were the recent gaseous hydrocarbon composition in peat pore waters and the overlying bottom sea waters.
Materials and methods
Peat, sediment, bottom and pore water samples were taken in the so-called “Entenschnabel” region of the German EEZ (Figure 1) close to the border with the Dutch EEZ during cruise MSM98 with RV MARIA S. MERIAN (VARIOSEEP project) in January 2021 (Römer et al., 2021b).
Sampling and sample preparation
Sediment cores containing peat-rich intervals were collected during four stations at three sites (Table 1). Peat, peat-enriched sediment, bottom and pore water samples were taken with a mini-multicorer (MIC) and a gravity corer (GC; all at ~43 m water depth) and originate from the central North Sea close to salt dome “Britta” (Figure 1 and Table 1). Sampling was achieved in so far unexplained seafloor depressions (see for more details Römer et al., 2021a; Figure 1). Sampling was preferentially done to understand pore water gas geochemistry and geomicrobiology in those areas, and peats were incidentally recovered as well. With regard to primary scientific goals of the project, the sediment cores were completely processed immediately upon recovery on deck and not subjected to a detailed sedimentological description. In total, five pure peat samples were prepared from the MIC and GC liners and samples were stored frozen until further processing for geochemical and biomarker analyses in the home laboratory.
Bulk and compound-specific geochemical data of peat samples from the central North Sea taken at 43 m water depth at the three sites A, B, and C.
LIMS: BGR lab code number; cmbsf: cm below sea floor; TOC: Total organic carbon; HI: hydrogen index; HC: hydrocarbon; OI: oxygen index; F/H: Fernenes/Hopenes (
Peat sample from different MUC tube or same depth in GC from same site.
Bottom water samples were taken from the MIC liners with a flexible silicon tube. The water samples were transferred into 122 ml glass bottles allowing only minor contact with air and water turbulence and avoiding gas bubbling. The samples were acidified (poisoned) with 2 ml 37% HCl and stored in the dark at +5°C until processing in the home laboratory. For the compositional analyses of hydrocarbons dissolved in pore waters, sediment samples were taken with a cut-off 3 ml plastic syringe through pre-drilled holes (5 cm distance) in the MIC liners. Sediment samples were transferred into 22 ml glass vials, poisoned with 17 ml of saturated NaCl solution (resulting in a headspace volume of 2 ml) and immediately crimped with a Teflon-coated butyl rubber seal. Analysis of stable carbon isotopic compositions (δ13C) of methane from the peat pore waters could not be realized with samples from the MIC due to insufficient sample volume having remained after compositional analyses. Therefore pore waters from peats sampled in parallel with a GC (station GeoB24443-1) were taken for these investigations. For δ13C-analysis of methane, 3 ml of bulk sediment was collected with cut-off syringes from the GC and transferred into 20 ml glass vials prefilled with 5 ml of 1M NaOH, thereby creating a headspace volume of 12 ml. Subsequently, the glass vials were sealed gas-tight with butyl rubber stoppers and crimp caps.
Analysis of bottom water and pore water gas geochemistry
Concentrations of volatile compounds dissolved in bottom water samples, that were taken with a hose from the MIC liners above the sediment surface, were determined applying a headspace equilibration technique described in detail by Schloemer et al. (2018). Twenty five milliliter of the 122 ml water sample were replaced by laboratory grade Helium and subsequently equilibrated at a temperature of 20°C for at least 2 h on a laboratory shaker. After equilibration, the headspace pressure was measured using a pressure transducer for the calculation of compound-specific concentrations considering Henry’s law (see below).
Sediment samples were not treated other than equilibrating to ambient temperature. The headspace volume of only 2 ml precluded a measurement of the headspace pressure. Therefore, it is reasonable to assume a normal (ambient) pressure of 101 kPa as absolute pressure.
Headspace gas composition was measured on a Trace 1310 gas chromatograph (Thermo Fisher Scientific) with heated valve system and column switching. The system and all included sample loops were evacuated (<6 mbar) with a diaphragm pump via several valves. One milliliter of sample was then injected into the sample loops. The individual components were quantified in parallel on three channels.
On channel 1, pre-separation of hydrocarbons (C1 through C6) was performed on a non-polar polysiloxane polymer column (Restek MX-1, 15 m, 0.28 mm ID, film thickness 3 µm). Molecular weight components >C7 were back-flushed. Full separation was performed on the downstream installed 50 m Al2O3 capillary column (0.32 mm ID, film thickness 5 µm). Both columns were operated non-isothermally. All components were detected on a Flame Ionization Detector (FID). Carrier gas was He.
On channel 2 and 3 other gas species were determined (i.e. CO2, He, Ne, H2, Ar, O2, N2, CH4, and CO). Because data from channels 2 and 3 are not reported, analytical details are not presented.
The concentrations of methane and of higher hydrocarbons dissolved in the bottom water samples were calculated using the partial pressure derived from fractional concentration and total headspace pressure, temperature of the sample, volume of headspace gas (25 ml), and remaining water (98 ml) applying the Henry’s Law constant of the gases under consideration. A correction for the salting-out effect was applied using a total salinity of 0.59 mol/l and the Setchenov constant for the analyzed components (Masterton and Lee, 1970; Morrison and Billett, 1952).
Calculation of hydrocarbons dissolved in the pore water of the sediment samples is similar, but requires additional corrections. Since the headspace gas was air, the methane concentration measured in the headspace was corrected for the methane concentration in air (1.85 ppm). Subsequently, the corrected dissolved hydrocarbon concentrations were normalized to the true pore water volume in the vial. Given that 3 ml of sediment were introduced to the vial and applying measured depth dependent porosities between 80% and 90% (see Table 2) the gas volume originated from 2.4 to 2.7 ml of pore water. The salting-out effect was calculated using a total salinity of 7 mol/l since a saturated NaCl-solution was added to the sample.
Gas geochemical data of pore waters and bottom waters taken from mini-multicores (MIC) from the central North Sea from the three sites A, B, and C in January 2021.
cmasf: cm above sea floor; cmbsf: cm below sea floor; nd: not detected.
The relative error of the gas chromatographical analysis is around ±3%, and for the total analyses of waters typically ~±10%. Due to the additional sources of error for the sediment samples (correction for air-derived methane, pore volume normalization), the overall error of concentrations of hydrocarbons dissolved in pore water is certainly higher, but cannot be accurately stated.
The gas chromatography-isotope ratio mass spectrometry (GC-IRMS) analyses of dissolved methane were performed using about 8 ml of 12 ml headspace of the NaOH treated 3 ml peat samples. Analytical details can be found elsewhere (Stawiarski et al., 2019). Due to the analytical procedure, which was established for (marine) samples with generally high methane concentrations, our data have a significant uncertainty. The peat-associated methane was analyzed from an air headspace, including 1.85 ppm atmospheric methane (δ13C of ~−47‰; Ganesan et al., 2018). As total methane concentrations in the headspace were ~4 ppm (data not shown), a significant fraction of the analytes represented atmospheric methane. Assuming an about 40/60 (air methane/sample methane) relationship, isotopic data were corrected with the atmospheric methane fraction. Uncertainties in the concentration data are, however, high and isotope corrections have thus to be taken with care.
Peat bulk geochemistry (C/S analyses and Rock-Eval pyrolysis)
Aliquots of dried (40°C for 24 h) sediment and peat samples were crushed and ground. Total organic carbon (TOC), carbonate carbon (Ccarb), and sulfur contents were measured separately on a LECO CS-230 (Leco Instrumente, Germany). TOC was determined on decalcified (10% HCl at 80°C; dried at 40°C for 18 h) sample aliquots. Between 10 and 40 mg of each sample was combusted in a high-frequency induction furnace in an oxygen atmosphere and concentrations were determined by use of the absorption signal at the IR detector. The instrument was calibrated using commercially available standards (LECO). Reproducibility of the measurements was ±0.02%. Rock-Eval pyrolysis was performed on a Rock-Eval 6 analyzer using a standard program described earlier (Espitalié et al., 1977; Lafargue et al., 1998): isothermal start at 300°C for 3 min, then applying a heating rate of 25 °C/min up to 650°C. Initial sample weights were ~15 mg to prevent oversaturation of the FID. Hydrocarbons released isothermally at 300°C are presented as “S1” and those released between 300°C and 650°C as “S2” yields. CO2 released up to 400°C is presented as “S3.” Reproducibility of the hydrocarbons and CO2 determination was better than 5%. TOC-δ13C values of the peats were determined after decalcification with HCl with an elemental analyzer coupled with a Thermo Delta V advantage IRMS.
Biomarker compositions and compound-specific stable carbon isotopic compositions (δ13C)
Organic compounds were extracted three times from 3 g of the ground peat samples using a 10 ml mixture of dichloromethane-methanol (DCM:MeOH, 8:2; v:v) for 15 min in an ultrasonic bath. The triplicate extracts were combined and dried under a gentle N2 stream at 40°C. Elemental sulfur was subsequently removed from the extract diluted in DCM by treatment with activated (10% HCl at 60°C for ~1 h) granulated copper, before they were dried and weighed. The extracts were fractionated on chromatographic columns filled with activated silica gel (240°C for 12 h) into aliphatic, aromatic, and polar fractions using isohexane, isohexane/DCM (2:1; v:v), and DCM/methanol (2:1; v:v) as solvents, respectively. All fractions were dried under a N2 stream and weighed. The distribution of compounds contained in the aliphatic fractions was determined by gas chromatograph-FID/mass spectrometer (MS) with an Agilent-7890 gas chromatograph equipped with a Programmable Temperature Vaporizing inlet splitting on two 50 m DB-1 columns (Agilent; 0.2 mm inner diameter; 0.11 μm film thickness), coupled to an FID and to an Agilent 7000 QQQ MS, each. He was used as carrier gas. Measurements were carried out in full scan mode. Compounds were identified by comparison of mass spectra, retention times, and measurement of standard compounds (e.g. hop-17(21)-ene, hop-22(29)-ene, neohop-13(18)-ene). Retention times of hopenes on the DB-1 column was slightly different than on other gas chromatography-columns (e.g. Elvert et al., 2016). Stable carbon isotope analysis (δ13C) of aliphatic compounds was made via gas chromatograph-C-IRMS using a coupled Trace GC Ultra, ConFlow IV, and a Delta V Advantage IRMS (all Thermo Fisher Scientific). The gas chromatograph was equipped with a split/splitless inlet and a DB-1 column (Agilent; 30 m, length; 0.25 mm, inner diameter; 0.25 μm film thickness) coupled to the IRMS via a Combustion III interface (Thermo Fisher Scientific). Measurements were carried out in duplicate. Accuracy of the analysis was checked by replicate measurements of an
Results
Bulk geochemistry
Table 1 shows bulk geochemical and compound-specific data of the five peat samples. With the exception of the peat sample from Site C (41-01; 17.2 wt.% TOC) all samples shared TOC >30 wt.%. Sulfur was abundant in all samples (5–19 wt.%). Rock-Eval pyrolysis confirmed that samples were composed of Type III kerogen and Tmax values demonstrated that they are immature as expected (peat maturity stage; Table 1). Sample 41-01 from Site C showed a relatively slightly elevated hydrogen index (HI) and Tmax (Table 1). TOC δ13C values ranged between −27.4‰ and −29.6‰ PDB.
Dissolved hydrocarbons in peat and sediment pore and bottom waters
Bottom waters overlying peat-rich sediments had concentrations of dissolved methane of ~2 nM. Within the sediment and peat, dissolved methane was found to be elevated compared to bottom waters and, where samples from different depths were retrieved, concentrations increased with sediment depth toward the peat layer at the bottom of the cores (Figure 2 and Table 2). Calculated concentrations of dissolved hydrocarbons were found to be highest in the peat layer. At Site A (Figure 1), changes in the sediment with depth could be studied in detail: In the peat layer, methane concentrations were ~250 nM and ethane, ethene, propane, propene, and

Site A gas geochemistry from a MIC (station 39-02): size of circles for compositions of hydrocarbons dissolved in porewaters correspond with concentrations of dissolved methane (not holding for bottom water). Calculations of concentrations of dissolved CH4 in sedimentary pore water base on porosity data from samples of core 39-02 and for peat from the GC taken for porosity determinations (slightly different to the samples studied for dissolved hydrocarbons). For color codes for individual hydrocarbons refer to the large data circle.
Concentrations of methane in the headspaces of the studied samples were generally too low after compositional analyses for reliable δ13C analyses with the analytical setup available. However, six larger gas samples (with 12 ml air headspace) from peat samples from GC station no. GeoB24443-1 were available for δ13C of methane determinations of peat-related pore waters. δ13C values of the methane were ~−46‰ (five samples) to ~−60‰ PDB (one sample).
Composition of aliphatic biomarkers and carbon isotopic signatures
The aliphatic extract fractions were dominated by long chained
δ13C values of aliphatic biomarkers, such as
Discussion
Rising temperatures during the early Holocene caused the melting of ice sheets and rise of sea level and of groundwater along the coastlines. This also caused flooding of the former “Doggerland” in the today’s central North Sea. During this interval peats were widely formed within coastal fens and the hydrological and climatological change is recorded in peats and surrounding sediments. Peats used for this study were not dated, but they were deposited between close-by 14C-dated shallow basal peat layers (samples “Dogger Bank,” “Elbe palaeovalley (N)”) from Vink et al. (2007) and presumably formed at a similar age interval. Such peats were dated between 9.7 cal. ka BP (“Dogger Bank”) and 10.9–9.2 cal. ka BP (“Elbe palaeovalley (N)”). Peats in the vicinity had similar ages and ages of the basal peats decreased toward the German coast (Vink et al., 2007). A peat sample (“RE29”; GPDN, 2013) that was taken even closer to our samples confirms that our studied peats have approximate ages of 10 cal. ka BP.
Modern biogeochemical processes in peat-rich sediments in the central North Sea
Recent peats host various microbial groups and this holds also for buried submarine peats like those in the North Sea, which also host microoorganisms executing methane formation and oxidation (Lippmann et al., 2021). However, bottom water directly overlying the peat-containing sites, which are located a few hundreds of meters NNE of the salt dome “Britta” that were investigated in this study, was relatively poor in dissolved methane. Small amounts of ethane (and ethene; and no C2+ hydrocarbons) were also observed in one out of the two bottom water samples at station 39-02 taken near “Britta”. The relatively low methane concentrations in the bottom water – during summer 2019 concentrations were found ~10 times higher (Römer et al., 2021a) – can be best explained by a high-turbulence winter effect. During winter season a thermocline, which typically prevents methane discharging from methane-laden sedimentary pore waters into shallow water layers, is not developed. Indeed, pore waters in shallow sediments were enriched in dissolved methane compared to bottom waters (2 nM CH4) overlying the sediments. For instance, at Site A (station 39-02) CH4 concentrations were ~60 nM at 2 cmbsf and increased with depth up to ~250 nM in the peat layer at 12 cmbsf (Figure 2 and Table 2). However, methane concentrations are not elevated in comparison with marine sediments elsewhere and also with other sediments in the central North Sea. Lippmann et al. (2021) reported average CH4 concentrations of sediment pore waters of 2.1 µM at ~2 to 5 mbsf. Irrespective of the different sampling depths applied in the two studies, this difference in methane concentrations suggests that methane in the region NNE of salt dome “Britta” is not much enriched compared to the situation in the broader area. However, CH4 concentrations in peat pore waters investigated in this study are higher than those at similar sediment depths in background cores collected 4 km SSW and 15 km NE to the sites studied herein (e.g. for core MSM98 48-01 (55°06′46.84″N; 04°46′49.836″E)) and core MSM98 49-01 (55°15′53.19″N; 04°56′19.44″E) we observed ~10–30 nM CH4 for both cores), which were virtually devoid of peats and showed much lower TOC contents (0.1–0.2 wt.%, not shown).
In addition to the methane enrichments, higher hydrocarbons (C2–C4), although in very low concentrations, were also observed in our peat samples (Table 2). This holds particularly for ethane (and ethene). Higher hydrocarbons were also found in sediment samples overlying peat layers, but the relative abundance of methane to ethane increased slightly with decreasing depth and also concentrations differed between samples with the lowest in the bottom waters (Table 2). Because of those differences the presence of an actively gas seeping system underlying the study sites is unlikely, as this would likely lead to more uniform hydrocarbon concentrations and compositions. Further, a simple two endmember mixing (bottom waters versus pore waters from peat) cannot explain the hydrocarbon composition in the shallow sediments in between, for example, at 2 and 7 cmbsf (Station A; Figure 2). Classical thermogenic hydrocarbon formation in the peat can be excluded, because of the low thermal maturity (see Tmax in Table 1) and their depositional history including permanent shallow burial depth in relatively shallow waters. On the other hand, the hydrocarbon composition in the peat is different to what is expected in normal marine sediments as light hydrocarbons in shallow marine sediments are commonly nearly entirely composed of biogenic methane (with C1/C2 ratios >1000; e.g. Bernard et al., 1977). Unfortunately, concentrations of methane in the distinct sediment and pore water samples used for this study, which remained after compositional analyses, were too low for analyses of stable C and H isotopes. However, methane in the headspaces overlying peat samples, which were taken for other analyses during MSM98, allowed for an assessment of methane’s gas geochemical characteristics in peats sampled with a GC (station GeoB 24443-1; no further gas data available) at the same position as station A. It is therefore likely that peats in the GC and those at 12 cmbsf at station A (Table 1 and Figure 2) represent the same peat layer. However, δ13C-CH4 values ranging from ~−46‰ to −60‰ PDB, with the most negative value being observed in the shallowest peat sample (data not shown), points at an inconclusive origin. It can either be explained by a mixed thermogenic/biogenic, by a degraded biogenic source (Whiticar, 1999), by a (not-well understood) low-temperature thermal gas generation (Rowe and Muehlenbachs, 1999), or by a biogenic methane specific to (boreal) peats. Indeed, biogenic methane from boreal wetlands is obviously less depleted in 13C (−53‰ to −71‰; Ganesan et al., 2018) than methane from strictly marine sediments (Whiticar et al., 1986). Biogenic higher hydrocarbons, particularly ethane, were occasionally observed in marine sediments (Bernard et al., 2013; Hinrichs et al., 2006; Oremland et al., 1988; Pape et al., 2020; Xie et al., 2013), but the relative abundance of methane to ethane (C1/C2) in samples reported to date (e.g. ~300 in Bernard et al., 2013) was much higher than the ratio which we observed in the peat samples from close to the salt dome “Britta” (C1/C2: ~5–20; Table 2). A normal, widespread marine microbial source is, however, unlikely to be the exclusive origin of the light hydrocarbons, which is also supported by the occurrences of propane for example. As pointed out, it is unlikely that the studied sediment column is currently charged with thermogenic gas even though it cannot be excluded from the geochemical data that part of the gas stems from an underlying oil-mature source rocks (e.g. the Kimmeridge clay; Müller et al., 2020). Another option would be that the gas was desorbed from charcoal, which is known to occur in basal peats at “Doggerland” (Krüger et al., 2017). We consider this process unlikely as well, because adsorption is commonly too strong to release gases by the methane desorption method applied in this study (exclusive addition of NaCl-saturated water; see method section). The occurrence of unsaturated light hydrocarbons (Table 2) additionally conflicts with a thermogenic origin. Ethene and propene are relatively abundant in our samples and are (biogenic) microbial precursors of ethane and propane, which often co-occur in sediments (Bernard et al., 2013; Kvenvolden and Redden, 1980; Kvenvolden et al., 1981; Vogel et al., 1982). The underlying microbial processes are still insufficiently understood, but ethene to ethane conversion has been shown for methanogenic cultures (Xie et al., 2013).
In conclusion, the partly slightly negative δ13C values of methane (as low as −60‰), the co-occurrence of low molecular weight hydrocarbons with their unsaturated analogs and their enrichment in the peat layers make an origin of the hydrocarbons from the peat or peat-influenced microbial turnover most likely. The presence of saturated/unsaturated hydrocarbon doublets in recent peats (from Pennsylvania) was previously shown (Kim, 1978) and support this interpretation.
The peat-forming environment on “Doggerland”
Four of the five peat-rich sediments jointly demonstrated relatively low HI (~150–190 mg HC/g TOC) and high OI (~110–130 mg CO2/g TOC) values (Figure 3 and Table 1). The general narrow ranges of the HI and OI values indicate that the samples were formed in comparable terrestrial settings, which are generally characterised by Type III organic matter. However, the only patchy recovery of generally small peat samples (of few centimeters thickness) does not allow studying chronological changes on “Doggerland” in this study in detail. It, however, presents geochemical snapshots from three spatially related sites.

Rock-Eval data of peat samples from the North Sea. Samples are relatively homogenous except for the sample from Site C, which demonstrated a slightly enhanced Hydrogen Index (HI).
The peat sample from Site C (sample 41-01) had a slightly higher HI value (258 mg HC/g TOC), which corresponds to a higher Tmax value (and the relatively lowest TOC content; Table 1). As maturities of all samples are likely similar, this may be an indication of slight local differences in the paleosetting and, thus, composition of the kerogen.
All peats contain significant amounts of sulfur maximizing at ~19 wt.% (sample 40-01 from Site A) and recent microbial sulfate reduction can be excluded as driver for the sulfur formation. This is evident from the recent relatively low sulfate reduction rates (Martin Krüger; personal communication) and from the fact that peat-surrounding sediments were relatively poor in sulfur (<0.1 wt.% S; not shown). Holocene peats from lake settings in NW Germany and also from coastal environments, all in NW-Germany, were described to be partly highly enriched in (pyrite-bound) sulfur (Dellwig et al., 1998; Gerdes et al., 2003; Regnery et al., 2013). This characteristic is explained by effective sulfate-reduction after flooding of peats in the lagoon-like settings and release of H2S, followed by the precipitation of pyrite from the H2S and Fe, the latter being transported via rivers from the Hinterland (Dellwig et al., 2001). It is likely that the peats from “Doggerland” experienced a similar geochemical history and that the high sulfur contents are an indicator of the flooding with sea water in the aftermath of postglacial peat formation.
Indications for a changeable environment on a small spatial scale
Distributions of hydrocarbon biomarkers in companion with compound-specific stable carbon isotope signatures excellently record major peat-forming plants and related microbial processes. In accordance with other studies, we found evidence for a spatially variable boreal setting, which fostered the growth of peat-forming and other higher plants during the Late Glacial and Early Holocene as has been documented by pollen (Krüger et al., 2017) and

Biomarker composition of peat extracts from the three sites shown from total ion currents during mass spectrometric analyses. Numbers at compounds indicate numbers of carbon atoms.

Detailed view of biomarkers extracted from sample 39-02 MIC from Site A, shown as total ion current (TIC, top) and by selected ions to highlight crucial components. *Identified by the mass spectrum from Howard et al. (1984) and the elution order published in Nytoft and Larsen (2001). Hop-17(21)-ene and diploptene were identified by comparisons with mass spectra and co-elution with authentic standards.
Bacterial hopenes and hopanes are another important class of lipids, which are present in peats from “Doggerland” (Figure 5). While the suite of compounds with the predominating trisnorhop-17(21)-ene, hop-17(21)-ene and diploptene (hop-22(29)-ene) was found in all samples, relative abundances slightly differ between the samples. Fernenes of which fern-7(8)-ene was tentatively identified as major compound of the fernene compound class. Fernenes are structurally related to hopenes, are widespread in certain bacterial groups (e.g. Howard et al., 1984), and were occasionally reported from microbially influenced environmental settings (Brassell et al., 1981; Nytoft and Larsen, 2001). The very low recent methane concentrations, the stable carbon isotopic composition of methane in the recent pore waters, and the high concentrations of hopenes and fernenes in the peats, however, strongly argue against recent methanotrophy as source for these compounds. Bacterial hopanoids present in the peats are, thus, likely and among other bacteria produced by aerobic bacteria involved in the methane turnover during peat formation. Peats are settings in which methanogenesis as well as (aerobic) methanotrophy are important processes. The latter process, and the magnitude of methanotrophy and/or water table fluctuations (Huang et al., 2018) during the time of peat formation, is recorded in the stable carbon isotope signature of hopenes (Pancost et al., 2002). We observed strongly varying δ13C values of hopenes minimizing at ~−63‰ at Site B (Figure 6). Relative depletions in 13C were also found for other hopenes and the tentative fern-7(8)-ene (as negative as −58‰) demonstrating the role of their bacterial source organisms in the methane turnover. As demonstrated by strongest relative depletions in 13C, methanotrophic biomass is particularly abundant at Site B. Nevertheless, bacterial compounds were also slightly depleted in 13C at the other sites.

Compound-specific stable carbon isotope signatures of selected biomarkers from North Sea peats from the three studied sites A, B, and C. Sampling depths in cm at individual sites are given in the legend.
Interestingly, the isotopic composition of individual compounds at the different sites differed, whereas distributions were relatively similar. This indicates variable settings with physiologically flexible bacterial compositions responding on fluctuating methane flux and/or other factors like water table changes (Zhang et al., 2022). The calculated Δδ between hopenes (diploptene) and plant wax alkanes of ~30‰ is at the lower end of other Middle-to-Early Holocene and Late Glacial peats in a global data set (Inglis et al., 2019). It is, however, much higher than for Recent (<2 ka) peats (Inglis et al., 2019), further hinting at the Late Glacial/Early Holocene age of the peats studied herein. Our data demonstrate that 13C-depleted (methanogenic) methane was abundant at certain sites on “Doggerland,” where it fostered activity of aerobic methanotrophic bacteria.
The peat sample from Site C (sample 41-01) documents an additional paleoenvironment. Whereas the majority of biomarkers is also found in the other samples, this sample showed the lowest TOC (17.2 wt.%; Table 1). Moreover, the relatively highest HI values argues for contributions from aquatic algae/bacteria to the bulk organic matter. Indeed, a likely and for peat-settings uncommon diatom-stemming highly branched isoprenoid (C20 HBI) was found in abundances that are similar to plant waxes (Figure 4). This compound was first ascribed to green algae (Rowland et al., 1985) and to cyanobacterial mats (Kenig et al., 1995). The C20 HBI tentatively identified in this study was, however, described from different Holocene lake settings (Aichner et al., 2010; Elvert et al., 2016). Recently it was suggested that it originate from the freshwater diatom
Nevertheless, our biomarker data demonstrate the variability of the deposition setting, which was controlled by peat-forming plants, bacterial methane-turnover, development of freshwater lakes, and a brackish/marine overprint. Further work aiming at extending the knowledge from our snapshot analyses and, thus, a more detailed understanding of methane turnover in peats from “Doggerland” should also involve polar lipid constituents, a robust time resolution and may be accompanied by palynomorph data.
Conclusion
Peats were an integral part at “Doggerland” and other coastal areas in the today’s central North Sea deposited during the Holocene (after the Last Glacial Maximum). For the first time, biomarkers contained in peats cored at three sites in the German Exclusive Economic zone at ~43 m water depth were studied. Relatively high sulfur concentrations (up to 19 wt.%) along with biomarkers specific for sulfur-driven transformation processes typify peat bulk geochemistry and document extensive past marine sulfate reduction after flooding caused by sea-level rise. Biomarker analyses record peat-forming plants, freshwater diatoms, and particularly bacterial methanotrophic turnover (e.g. by 13C-depleted hopenes and fernenes) during deposition at times, when boreal wetlands extended, global atmospheric methane concentrations were relatively high and methane became progressively 13C-depleted (Fischer et al., 2008; Loulergue et al., 2008).
Concentrations of dissolved methane in the sediment pore waters were generally low, but slightly enriched in the peats (~250 nM). Methane was accompanied by ethane, ethene, propane, and propene, which are most likely originate from microbial/diagenetic processes as well. The exact source and the role of peats for hydrocarbon concentrations in the water column overlying peat-rich settings remains a topic for further studies in the North Sea or at other peat-rich flooded settings.
While further investigations should focus on the changes during geological times in more detail, our study presents novel organic geochemical snapshots into the former methane cycle during Holocene peat formation in the today’s North Sea.
