How dry was the Younger Dryas ? Evidence from a coupled δ 2 H – δ 18 O biomarker paleohygrometer applied to the Gemündener Maar sediments , Western Eifel , Germany

Causes of the Late Glacial to Early Holocene transition phase and particularly the Younger Dryas period, i.e. the major last cold spell in central Europe during the Late Glacial, are considered to be keys for understanding rapid natural climate change in the past. The sediments from maar lakes in the Eifel, Germany, have turned out to be valuable archives for recording such paleoenvironmental changes. For this study, we investigated a Late Glacial to Early Holocene sediment core that was retrieved from the Gemündener Maar in the Western Eifel, Germany. We analysed the hydrogen (δ2H) and oxygen (δ18O) stable isotope composition of leaf-wax-derived lipid biomarkers (n-alkanes C27 and C29) and a hemicellulose-derived sugar biomarker (arabinose), respectively. Both δHn-alkane and δOsugar are suggested to reflect mainly leaf water of vegetation growing in the catchment of the Gemündener Maar. Leaf water reflects δ2H and δ18O of precipitation (primarily temperaturedependent) modified by evapotranspirative enrichment of leaf water due to transpiration. Based on the notion that the evapotranspirative enrichment depends primarily on relative humidity (RH), we apply a previously introduced “coupled δHn-alkane–δOsugar paleohygrometer approach” to reconstruct the deuterium excess of leaf water and in turn Late Glacial–Early Holocene RH changes from our Gemündener Maar record. Our results do not provide evidence for overall markedly dry climatic conditions having prevailed during the Younger Dryas. Rather, a two-phasing of the Younger Dryas is supported, with moderate wet conditions at the Allerød level during the first half and drier conditions during the second half of the Younger Dryas. Moreover, our results suggest that the amplitude of RH changes during the Early Holocene was Published by Copernicus Publications on behalf of the European Geosciences Union. 714 J. Hepp et al.: How dry was the Younger Dryas? more pronounced than during the Younger Dryas. This included the occurrence of a “Preboreal Humid Phase”. One possible explanation for this unexpected finding could be that solar activity is a hitherto underestimated driver of central European RH changes in the past.

Abstract.Causes of the Late Glacial to Early Holocene transition phase and particularly the Younger Dryas period, i.e. the major last cold spell in central Europe during the Late Glacial, are considered to be keys for understanding rapid natural climate change in the past.The sediments from maar lakes in the Eifel, Germany, have turned out to be valuable archives for recording such paleoenvironmental changes.
For this study, we investigated a Late Glacial to Early Holocene sediment core that was retrieved from the Gemündener Maar in the Western Eifel, Germany.We analysed the hydrogen (δ 2 H) and oxygen (δ 18 O) stable isotope composition of leaf-wax-derived lipid biomarkers (n-alkanes C 27 and C 29 ) and a hemicellulose-derived sugar biomarker (arabinose), respectively.Both δ 2 H n-alkane and δ 18 O sugar are suggested to reflect mainly leaf water of vegetation growing in the catchment of the Gemündener Maar.Leaf water re-flects δ 2 H and δ 18 O of precipitation (primarily temperaturedependent) modified by evapotranspirative enrichment of leaf water due to transpiration.Based on the notion that the evapotranspirative enrichment depends primarily on relative humidity (RH), we apply a previously introduced "coupled δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach" to reconstruct the deuterium excess of leaf water and in turn Late Glacial-Early Holocene RH changes from our Gemündener Maar record.
Our results do not provide evidence for overall markedly dry climatic conditions having prevailed during the Younger Dryas.Rather, a two-phasing of the Younger Dryas is supported, with moderate wet conditions at the Allerød level during the first half and drier conditions during the second half of the Younger Dryas.Moreover, our results suggest that the amplitude of RH changes during the Early Holocene was Published by Copernicus Publications on behalf of the European Geosciences Union.

Introduction
In order to evaluate the relevance of man-made climate change in the future, it is of great importance to study and understand large and rapid climate fluctuations in the past.Many studies have focused on the Late Glacial to Early Holocene transition phase, a period with various expressions in temperature, atmospheric circulation and hydrology worldwide (Alley, 2000;Brauer et al., 2008;Denton et al., 2010;Partin et al., 2015;Wagner et al., 1999).Explanation for the Younger Dryas (YD) period, i.e. the major last cold spell in central Europe during the Late Glacial just before the onset of the Holocene warm period (Denton et al., 2010;Heiri et al., 2014;Isarin and Bohncke, 1999), has long been considered crucial for understanding rapid natural climate change in the past (Alley, 2000).The sediments from maar lakes in the Eifel, Germany, have turned out to be valuable archives for paleoenvironmental reconstruction by providing high-resolution palynological, sedimentological and geochemical records for climate, vegetation and landscape history (Brauer et al., 2008;Brunck et al., 2015;Litt et al., 2003;Litt and Stebich, 1999;Sirocko et al., 2013;Zolitschka, 1998).
Within this study, we aimed at applying the coupled δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach to the Late Glacial-Early Holocene sediment cores of the Gemündener Maar.More specifically, we addressed the following objectives: (i) source identification of the sedimentary organic matter and the investigated n-alkanes and sugars (aquatic vs. terrestrial), (ii) reconstructing leaf water isotope composition based on compound-specific δ 2 H and δ 18 O values of the nalkane and sugar biomarkers, (iii) reconstructing RH changes using the coupled δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach, and (iv) inferring implications for central European paleoclimate history from the established Gemündener Maar RH record.

The Gemündener Maar and sampling
The Gemündener Maar is located in the Eifel volcanic fields in western Germany at an altitude of 407 m a.s.l.(50 • 10 39.853 N, 6 • 50 12.912E; Fig. 1a and b; Sirocko et al., 2013).The maar was formed during a phreatomagmatic explosion within the local Devonian siltstone (greywacke) around 20-25 ka (Büchel, 1993).The lake is 39 m deep at its maximum and has a diameter of roughly 300 m.Due to its formation conditions the lake is almost circular with a lake surface area of 75 000 m 2 and is surrounded by a small catchment (Fig. 1b), with an area of 430 000 m 2 (Scharf and Menn, 1992).The lake is fed by precipitation and groundwater (no surface inflow and outflow present).The sediments are, accordingly, not affected by fluvial sediment input.The catchment area is furthermore steep and densely vegetated with broadleaved trees (Fig. 1c).The investigated samples were taken from the 8 m Gemündener Maar core (GM1), which was taken at approximately 20 m water depth near the centre of the maar (Fig. 1b) with a Livingston piston corer (UWITEC, Mondsee, Austria).The GM1 core was retrieved from a terrace on the steep slope of the maar exactly in a fan of an underwater erosion gully structure.The core is part of the Eifel Laminated Sediment Archive Project of the Institute for Geoscience at Johannes Gutenberg University Mainz (Sirocko et al., 2013(Sirocko et al., , 2016)).

Bulk analysis and pollen analysis
Bulk analyses were carried out on 112 samples, covering a section of 606 to 727 cm depth of the Gemündener Maar GM1 core.Total carbon (TC) and nitrogen (N), bulk carbon isotope composition (δ 13 C TC ), and nitrogen isotope composition (δ 15 N) were determined at the Institute of Agronomy and Nutritional Sciences, Soil Biogeochemistry, Martin Luther University Halle-Wittenberg, using EuroVector EA 3000 elemental analyser (Hekatech, Wegberg, Germany) coupled via a Conflo III Interface to a Delta V Advantage isotope ratio mass spectrometer (IRMS; both from Thermo Fisher Scientific, Bremen, Germany).Additionally, total organic carbon (TOC) and bulk δ 13 C of the total organic carbon (δ 13 C TOC ) were assessed after removal of carbonate with 32 % hydrochloric acid (HCl) fumigation followed by a neutralisation step with moist sodium hydroxide, both for 24 h under 60 • C and vacuum conditions.This allows calculating TOC/N atomic ratios.Laboratory standards from the International Atomic Energy Agency (IAEA) as well as from the United States Geological Survey (USGS) with known total carbon, nitrogen, 13 C and 15 N contents (IAEA N2, IAEA CH6, IAEA NO3, IAEA CH7, IAEA 305A, USGS 41) were used for calibration.The 13 C and 15 N contents are expressed in the common δ notation as relative to an international standard (δ 13 C: Vienna Pee Dee Belemnite, VPDB; δ 15 N: atmospheric N 2 , air).
For pollen analysis, 16 samples were investigated covering the relevant depth section.Each sample covered a depth range of 1 cm.Preparation was conducted by Frank Dreher according to standard procedures at the laboratory of the Group of Climate and Sediments, Institute of Geosciences, Johannes Gutenberg University of Mainz, using potassium hydroxide, HCl and hydrofluoric acid (Sirocko et al., 2016).Afterwards, acetic acid and a mixture of acetic anhydride and sulfuric acid (9 : 1) were used for acetolysis.The samples were then centrifuged at 3000 to 3500 rpm for 5 min and then sieved over a 200 and a 10 µm sieve.Afterwards, the samples were fixed with anhydrous glycerol for reliable identification, and a maximum magnification of 600 was used for counting the remains.Pollen results are reported in relative percentages (%).The investigated sediments are partially laminated.The first tie point to establish a chronology for the Gemündener Maar core is a radiocarbon-dated piece of charcoal in 727 cm core depth, which dates to 13 800 ± 110 a cal BP (Fig. 2d).This age is derived from a 14 C age of 11 950 ± 65 a BP as part of the Supplement of Sirocko et al. (2013), calibrated using Cal-Pal software (Weninger and Jöris, 2008) calculated with the IntCal13 calibration curve (Reimer et al., 2013).The uncertainty of the calibrated 14 C age represents the 68 % probability range.The second tie point is the clearly visible Laacher See Tephra between 673 and 680 cm core depth (Fig. 2a and  d).The latter can be used as chronological marker due to the varve-counted age of 12 880 a BP in the adjacent Meerfelder Maar (Brauer et al., 1999).The onset of the Younger Dryas period was set to 12 680 a BP (varve-counted in Meerfelder Maar sediments; Brauer et al., 1999;Litt et al., 2009) identified at a depth of 670 cm in the GM1 core due to a clear colour change (Fig. 2a and d).The onset of the Preboreal (Holocene) was found to date to 11 590 a BP in Meerfelder Maar by varve counting (Litt et al., 2009).This was used to wiggle-match the distinct changes in the pollen spectra (decreasing Poaceae, peaking Artemisia, increasing Pinus and Betula; Fig. 2b and c), the clear rise in TOC (Fig. 3a) and the colour change (Fig. 2a), which were identified at 643 cm depth (Fig. 2d).The Late Glacial to Preboreal (Holocene) transition is commonly well recorded in maar sediments from the Eifel region, i.e. clear changes in deposition as well as pollen pattern (Brauer et al., 1999;Litt et al., 2001Litt et al., , 2003;;Litt and Stebich, 1999), dated to 11 600 a BP, e.g. in Holzmaar, by a combination of varve counting and 14 C dating (Zolitschka, 1998).The last time marker used to constrain the age model is the middle of the sharp increase in Corylus (hazel) pollen at 622 cm depth (Fig. 2b, c and d).We used this sharp increase as a marker for the Preboreal to Boreal transition, which is varve-counted by Litt et al. (2009) to 10 740 a BP in the Meerfelder Maar sediments.The offset of 60 years to the varve-counted Holzmaar record of Zolitschka (1998), as it is presented by Litt et al. (2009), is within the uncertainty of placing the onset of the Preboreal in the Gemündener Maar Corylus curve.
The investigated core section from 607 to 694 cm depth therefore covers the time between ∼ 13 150 and 10 140 a BP, i.e. the Allerød, the Younger Dryas, the Preboreal, and the beginning of the Boreal, with regard to the biomarkers (Fig. 2a  and d).Assuming constant sedimentation rates between the markers, an average resolution of 51 a cm −1 can be calculated; the minimum and maximum resolution are 19 and 124 a cm −1 , respectively.The part above the Laacher See Tephra reveals a lower mean resolution (55 a cm −1 ) than the section below (30 a cm −1 ).

Biomarker and compound-specific isotope analysis
For δ 2 H analyses of n-alkanes as well as δ 18 O analyses of sugars, 59 samples were prepared from 607 to 694 cm depth of the Gemündener Maar GM1 core, in order to cover the core section with already high TOC content and the Late Glacial to Holocene transition (Figs. 2 and 3a).n-Alkanes were extracted from 1 to 6 g freeze-dried and ground samples by microwave extraction at 100 • C for 1 h, using 15 mL of solvent (dichloromethane and methanol, at a ratio of 9 : 1).The resultant total lipid extracts were separated over aminopropyl silica gel (Supelco 45 µm) filled pipette columns.Non-polar compounds (including n-alkanes) were eluted with n-hexane.The fraction was spiked with a known amount of 5α-androstanone, used as an internal standard.Identification and quantification was carried out on an Agilent MS 5975 (EI) interfaced with an Agilent 7890 GC equipped with a 30 m fused silica capillary column (HP5-MS 0.25 mm i.d., 0.25 µm film thickness) and a split-splitless injector operating in splitless mode at 320 • C. Carrier gas was helium and the temperature program was 1 min at 50 • C, from 50 to 200 • C at 30 • C min −1 , from 200 to 320 • C at 7 • C min −1 , and 5 min at 320 • C. Data recording comprised the total ion count (scan mode from m/z 40 to 600) and single-ion monitoring (m/z 57, 71, 85 and 99).Concentrations were calculated relative to the internal standard and to an external standard (n-C 21 to n-C 40 alkane mixture, Supelco), injected in different concentrations (40, 4, 1, 0.4 ng µL −1 ).
Prior to compound-specific isotope analyses, the n-alkanes were further purified.The non-polar fractions were passed over a pipette column filled with activated AgNO 3 impregnated silica gel and a pipette column filled with zeolite (Geokleen).After drying, the zeolite was removed using hydrofluoric acid and the n-alkanes were recovered by liquidliquid extraction with hexane.The purified n-alkane fractions were measured for their compound-specific stable hydrogen isotope composition (δ 2 H).The measurements were performed at the Institute of Geography, University of Bern on an IsoPrime 100 IRMS, coupled to an Agilent 7890A GC via a GC5 pyrolysis or combustion interface operating in pyrolysis modus with a Cr (ChromeHD) reactor at 1000 • C. Samples were injected with a split-splitless injector.The GC was equipped with 30 m fused silica column (HP5-MS, 0.32 mm inner diameter, 0.25 µm film thickness).The precision was checked by co-analysing a standard alkane mixture (n-C 27 , n-C 29 , n-C 33 ) with known isotope composition (Arndt Schimmelmann, University of Indiana), injected twice every six runs.The samples were analysed in three repetitions (except from the samples in 622 and 672 cm depth), and the analytical precision was generally better than 5 ‰.The stable hydrogen isotope compositions are given in the δ notation (δ 2 H n-alkane ) versus Vienna Standard Mean Ocean Water (VSMOW).The H + 3 -correction factor was checked every 2 days and stayed stable over the course of measurements at 3.14.The δ 2 H n-alkane values refer to the area-weighted  mean of the δ 2 H values of n-alkanes with 27 and 29 carbon atoms (n-C 27 , n-C 29 ) because of their relatively high abundance in the samples (Fig. 4a).
The sample preparation for δ 18 O analyses of hemicellulose-or polysaccharide-derived sugars followed standard procedures at the Institute of Agronomy and Nutritional Sciences, Soil Biogeochemistry, Martin Luther University Halle-Wittenberg, according to the method of Zech and Glaser (2009).The monosaccharide sugars were hydrolytically extracted from samples containing approximately 10 mg total organic carbon with 10 mL of 4 M trifluoroacetic acid at 105 • C for 4 h, applying the method described by Amelung et al. (1996).After filtration over glass fibre filters, the extracted sugars were cleaned using XAD-7 (to remove humic-like substances) and Dowex 50WX8 columns (to remove interfering cations).Afterwards, the purified samples were freeze-dried and derivatised by adding methylboronic acid (4 mg in 400 µL pyridine) for 1 h at 60 • C.
The compound-specific δ 18 O measurements were performed using a Trace GC 2000 coupled to a Delta V Advantage IRMS via an 18 O-pyrolysis reactor (GC IsoLink) and a ConFlo IV interface (all devices from Thermo Fisher Scientific, Bremen, Germany).Each sample was measured in threefold repetition, embedded in-between co-derivatised sugar standards at various concentrations and known δ 18 O values.The δ 18 O values of the samples are expressed in δ notation (δ 18 O sugar ) versus VSMOW.The measured δ 18 O sugar values were corrected for drift, amount and area dependency and also for the hydrolytically introduced oxygen atoms that form carbonyl groups with the C1 atoms of the sugar molecules (Zech and Glaser, 2009).Mean standard errors for the triplicate measurements of all 59 samples are 0.6 ‰, 0.7 ‰ and 0.7 ‰ for arabinose, fucose and xylose, respectively.The δ 18 O sugar values refer to the δ 18 O values of the monosaccharides arabinose, fucose and xylose (Fig. 4b).Rhamnose areas, or concentrations, were too low for reliable isotope measurements in most samples.

Conceptual framework of the coupled
δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach The coupled δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach was described in detail by Tuthorn et al. (2015) and Zech et al. (2013a).The most fundamental assumption of the approach is that the isotope composition of leaf water can be reconstructed by applying biosynthetic fractionation factors on the measured δ 2 H n-alkane and δ 18 O sugar values (Fig. 5).
The concept is furthermore based on the observation that the isotope composition of global precipitation plots typically close to the global meteoric water line (GMWL; δ 2 H prec = 8 • δ 18 O prec + 10; Dansgaard, 1964).In Germany, a local meteoric water line (LMWL Germany ) slightly deviating from GMWL was described by Stumpp et al. (2014) (δ 2 H prec = 7.72±0.13•δ 18O prec +4.90±0.01;Fig. 5), which we used as the baseline for our calculations.The quite similar LMWLs for Trier (δ 2 H prec = 7.81 ± 0.08 • δ 18 O prec + 5.06 ± 0.60) and Koblenz (δ 2 H prec = 7.80 ± 0.07 • δ 18 O prec + 2.68 ± 0.53) as well as the GMWL are additionally displayed in Fig. 5 for comparison.The local precipitation is the source for soil water and shallow groundwater, which in turn acts as source water for plants.During daytime, however, leaf water is typically 2 H-and 18 O-enriched compared to the source water due to evapotranspiration through the stomata (Fig. 5; Allison et al., 1985;Bariac et al., 1994;Walker and Brunel, 1990).The leaf water reservoir at the evaporative sites quickly achieves steady-state conditions (Allison et al., 1985;Bariac et al., 1994;Gat et al., 2007;Walker and Brunel, 1990).Thus, the isotope composition of the transpired water vapour is equal to the isotope composition of the source water utilised by the plants during the transpiration process.The evaporative enrichment of leaf water under steady-state conditions can be described via a Craig-Gordon model (e.g.Flanagan et al., 1991;Roden and Ehleringer, 1999) by the following expression (e.g.Barbour et al., 2004): where δ e , δ s and δ a are the hydrogen and oxygen isotope compositions of leaf water at the evaporative sites, in source water and in atmospheric water vapour, respectively; ε * is the equilibrium enrichment expressed as (1 − 1/α L/V ) × 10 3 where α L/V is the equilibrium fractionation between liquid and vapour in ‰; and ε k is the kinetic fractionation parameters for water vapour diffusion from intracellular air space through the stomata and the boundary layer, both for 2 H and 18 O, respectively; and e a /e i is the ratio of atmospheric vapour pressure to intracellular vapour pressure.When leaf temperature is equal to air temperature, the e a /e i ratio represents the RH of the local atmosphere.If the plant source water and the local atmospheric water vapour are in isotopic equilibrium, the term δ a −δ s can approximated by −ε * .Thus, Eq. ( 1) can be reduced to The kinetic fractionation parameters (ε k ) are typically related to stomatal and boundary layer resistances with respect to water flux (Farquhar et al., 1989).Since direct measurements of those plant physiological parameters can be hardly assessed in a paleo application, we used the kinetic enrichment parameters C k instead, derived from a more generalised form of the Craig-Gordon model, for describing the kinetic isotope enrichment for 2 H and 18 O, which leads to Eq. ( 3) (Craig and Gordon, 1965;Gat and Bowser, 1991): In a δ 2 H-δ 18 O diagram, the hydrogen and oxygen isotope composition of leaf and source water can be described as a local deuterium (d) excess = δ 2 H − 7.72 • δ 18 O (Stumpp where d e and d s are the d excess values of the leaf water at the evaporative sites and the source water, respectively, and the equilibrium (ε * 2 and ε * 18 ) and kinetic enrichment parameters (C 2 k and C 18 k ) are expressed for both isotopes.From Eqs. (1) to (4) the primary control of RH on the isotope composition of the leaf water is demonstrated when stomata are open through transpiration.If d e can be derived from compound-specific δ 2 H and δ 18 O measurements of the nalkane and sugar biomarkers, which derive δ 2 H e and δ 18 O e values for the purpose of calculating d e values via the equa-tion d e = δ 2 H e − 7.72 • δ 18 O e , the d s can also be approximated from the d excess of the LMWL Germany (= 4.9).Accordingly, Eq. ( 4) can be rearranged in order to calculate the RH of the local atmosphere normalised to leaf temperature as given by Eq. ( 5) (Hepp et al., 2017;Tuthorn et al., 2015;Zech et al., 2013a):  Merlivat (1978), who reported maximum values during the molecular diffusion process of water through a stagnant boundary layer.The assumption that maximum kinetic fractionation occurs seems to be most suitable for sedimentological application where a signal averaging over decades can be assumed (see above and discussion in Zech et al., 2013a).It should also be noted that ε k values of broadleaf trees and shrubs over broad climatic conditions are well within the range of the C 2 k and C 18 k values used, revealing 23.9 ‰ (±0.9) and 26.7 ‰ (±1.0) for δ 2 H and δ 18 O, respectively (derived from the Supplement of Cernusak et al., 2016).
The numerator of Eq. ( 5) describes the parallel distance between the d excesses of LMWL and leaf water at the evaporative sites, which is converted into RH values, while the denominator is a combination of the slopes of LMWL and the local evaporation line (LEL).This means in turn that the quantification with Eq. ( 5) is done by obtaining the distance between the source water points, calculated via the intersects between the individual LELs and the LMWL Germany , and the leaf water points.The underlying slope of those LELs can be derived from Eq. ( 6) via the Craig-Gordon model using the same assumptions as outlined above in a rearranged form (Eq. 6; Zech et al., 2013a).When using the fractionation parameters from above, the slope of the LEL is constant over time, independent of RH and equal to ∼ 2.8 (Eq.6).This agrees well with field and laboratory studies (Allison et al., 1985;Bariac et al., 1994;Gat et al., 2007;Tipple et al., 2013;Walker and Brunel, 1990).
In order to provide an uncertainty interval in terms of measurement precision covering the Gemündener Maar RH record, we calculated an error propagation for d e values according to Eq. ( 7), by using the analytical standard errors (SEs).Maximum and minimum values were then applied to Eq. ( 5) resulting in a lower and upper RH limit (blue-shaded area in Fig. 7a).For basic sedimentological characterisation, TOC, N, δ 15 N, δ 13 C TC and δ 13 C TOC values as well as the TOC/N atomic ratios (Fig. 3a to e) are displayed from 605 to 727 cm depth.TOC values range from 0.6 % to 19.7 %.N ranges from 0.1 % to 1.4 % and correlates highly significantly with TOC (r = 0.99, p < 0.001, n = 110).Higher TOC contents during the Allerød, Preboreal and Boreal likely reflect warmer conditions being favourable for terrestrial and aquatic biomass production, whereas lower TOC values during the Younger Dryas likely reflect less favourable conditions for biomass production and possibly increasing minerogenic sedimentation.Interestingly, the Late Glacial-Early Holocene TOC patterns seem not to be the same for all maar lakes, because the Meerfelder Maar shows a clear TOC two-phasing during the Younger Dryas (Brauer et al., 1999) and the Holzmaar is lacking an Allerød TOC maximum (Lücke et al., 2003).The δ 15 N values of the Gemündener Maar record range from 0 ‰ to 5 ‰, showing the maximum and minimum within the Allerød period.δ 13 C TC and δ 13 C TOC reveal values between −31 ‰ and −17 ‰ and −36 ‰ and −24 ‰, respectively.While δ 13 C TC shows maximum values at 703 cm depth, δ 13 C TOC is decreasing continuously from the beginning to the end of the Allerød, followed by increasing values during the Younger Dryas and the Preboreal and Boreal, interrupted by a short decrease around the beginning of the Holocene.δ 13 C TC clearly shows the presence of carbonate between 690 and 727 cm depth with less negative δ 13 C TC values compared to δ 13 C TOC values.TOC/N atomic ratios range from 6 to 19 with the end of the Allerød revealing increasing ratios, while the late Younger Dryas shows slightly decreasing ratios and the Preboreal is marked by the highest ratios.
The source of organic matter in lacustrine sediments of small lakes, as one of the most crucial questions and challenges when interpreting organic proxies from lacustrine sedimentary records (Meyers and Ishiwatari, 1993), can either be autochthonous (aquatic origin) or allochthonous (terrestrial origin).The TOC/N ratio and δ 13 C values are the most common proxies for sedimentary source determination.While non-vascular aquatic organisms often reveal C/N atomic ratios between 4 and 10 (due to low amounts of cellulose and lignin), vascular plants commonly show C/N atomic ratios of 20 and higher (Meyers and Ishiwatari, 1993).According to Meyers (2003), a TOC/N atomic ratio of 10 is often used as threshold for identifying aquatic versus terrestrial input (Fig. 3e).Accordingly, the input from terrestrial organic matter increased during the Allerød, decreased slightly during the Younger Dryas and was highest during the Holocene.The Gemündener Maar δ 13 C TOC values (Fig. 3d) are well within the range of C 3 land plants and lacustrine algae (Meyers and Lallier-Vergés, 1999); evidence for the occurrence of C 4 land plants is missing.Overall, no clear additional information about the sedimentary organic matter origin of the Gemündener Maar sediments can be inferred neither from δ 13 C TOC alone (cf.Lücke et al., 2003) nor by combining δ 13 C TOC with TOC/N ratios (cf.Meyers and Lallier-Vergés, 1999).When considering that both δ 13 C TOC and TOC/N values of terrestrial organic matter are additionally affected by mineralisation and degradation, resulting in more positive δ 13 C TOC values and lower TOC/N ratios (e.g.Zech et al., 2007), a straightforward interpretation of those proxies seems to be challenging.Similarly, δ 15 N has been investigated as proxy for sedimentary organic matter origin (Meyers and Ishiwatari, 1993;Meyers and Lallier-Vergés, 1999;Wolfe et al., 1999).However, numerous processes that may have an influence, like nitrogen uptake by plants, various nitrogen sources, discrimination during denitrification and diagenesis, complicate the use of δ 15 N as a direct source determination proxy.
Despite the uncertainties presented above, concerning the origin of bulk sedimentary organic matter in the Gemündener Maar, the origin of the sedimentary biomarkers, namely nalkanes and sugars, needs to be addressed.This is crucial because aquatic biomarkers incorporate the isotope composition of lake water, whereas terrestrial biomarkers incorporate the isotope composition of leaf water (Huang et al., 2004;Kahmen et al., 2013;Mügler et al., 2008;Sachse et al., 2004Sachse et al., , 2012;;Sauer et al., 2001;Tuthorn et al., 2014;Zech et al., 2013bZech et al., , 2014b)).With regard to the n-alkane biomarkers, high amounts of the chain lengths n-C 27 and n-C 29 are characteristic of the Gemündener Maar sediments.Such patterns are typical of epicuticular leaf wax layers of higher terrestrial plants (e.g.Eglinton and Hamilton, 1967).With regard to the sugar biomarkers, they were previously studied in detail by Hepp et al. (2016).According to the authors' own results and a compilation from the literature (including, e.g., Jia et al., 2008;Prietzel et al., 2013;Zech et al., 2012Zech et al., , 2014b)), relatively high amounts of arabinose are a good indicator of a primarily terrestrial origin (higher vascular plants) of the sugars.This interpretation is in agreement with the Gemündener Maar being a small lake with densely forested steep crater walls (Fig. 1c).We therefore conclude and suggest that arabinose as well as n-C 27 and n-C 29 in our Gemündener Maar record are primarily of terrestrial rather than  4b).They reveal similar trends overall (arabinose vs. xylose: r = 0.7, p < 0.001, n = 59; arabinose vs. fucose: r = 0.8, p < 0.001, n = 59; xylose vs. fucose: r = 0.8, p < 0.001, n = 59).All sugar records show a clear shift to more positive values at the Younger Dryas-Holocene transition.While xylose and fucose exhibit a change of ∼ 8 ‰ and 7 ‰, arabinose δ 18 O values show a less pronounced shift of ∼ 3 ‰ (changes are based on the mean δ 18 O values for the Younger Dryas compared to the Preboreal/Boreal period).Xylose is however slightly more negative throughout the Allerød and Younger Dryas compared to arabinose and fucose.Consistently less pronounced changes can be observed for the Allerød-Younger Dryas transition of 1.9 ‰, 1.7 ‰ and 0.9 ‰ for xylose, fucose and arabinose, respectively (based on the mean δ 18 O values for the Allerød compared to the Younger Dryas).A distinct minimum during the early Preboreal (633 cm depth) characterises all three δ 18 O sugar records.
From the study of tree rings, it is known that stem cellulose does not show the full leaf water 18 O enrichment signal.Barbour and Farquhar (2000) related this signal dampening to the proportion of unenriched source water contributing to the local synthesis water (p x ) and to the proportion of exchangeable oxygen during cellulose synthesis (p ex ).The latter is often assumed to be rather constant around 0.40, as estimated from leaf and wood cellulose of Eucalyptus globulus and values compiled from the literature (Cernusak et al., 2005), meaning that around 40 % of the oxygens in the stem cellulose exchanged.Already Helliker and Ehleringer (2002) compared the signal transfer from leaf water to the cellulose of tree stems with the signal transfer occurring in grasses.And Liu et al. (2016) reported signal dampening in the range between 34 % and 53 % for the C 4 grass Cleistogenes squarrosa.
Figure 4c illustrates that Poaceae pollen concentrations ranged between 11 % and 33 % during the Allerød and the Younger Dryas in the Gemündener Maar record.Hence, a correction for the 18 O signal dampening may be required in order to take these vegetation changes into consideration.A respective correction procedure based on mass balance considerations is given in Eq. ( 10) in order to adjust δ 18 O l to δ 18 O l # : The correction presented in Eq. ( 10) is based on assumptions that 40 % (0.4) of the leaf water enrichment is lost during hemicellulose biosynthesis of grass leaves, which is well within the range of values presented in the literature for cellulose synthesis in tree rings and grasses (Cernusak et al., 2005;Liu et al., 2016).Furthermore, the Poaceae pollen concentration in percentage is used to calculate the fraction of non-grassy pollen (f non-grasses = (100 − Poaceae)/100) corresponding to the non-grassy biomarker contribution, which may serve as a rough approximation.For a paleo application, δ 18 O s remains a priori unknown.Therefore, the intercept between the individual LELs (Eq. 6) and the LMWL of Germany were used to generate δ 18 O s values.Note that the signal dampening effect described here for cellulose synthesis is likely not fully applicable to our approach using the sugar biomarker arabinose.In fact, pentoses like arabinose are biosynthesised via decarboxylation of the carbon at position six (C6) from glucose (Altermatt and Neish, 1956;Burget et al., 2003;Harper and Bar-Peled, 2002).Waterhouse et al. (2013) showed that the oxygens at C6 position in glucose moieties are most strongly affected by the exchange with a local water medium of 80 %, as indicated by heterotrophic cellulose synthesis.Thus, arabinose has lost a strongly exchanged (dampened) oxygen and the remaining pentose shows less 18 O signal dampening.With regard to the ε 2 bio value of −160 ‰, this biosynthetic fractionation factor is confirmed by climate chamber studies of dicotyledonous plants (Kahmen et al., 2011(Kahmen et al., , 2013;;Tipple et al., 2015).However, the latter studies also reveal a range of ∼ 35 ‰, interpreted as species-specific effects during nalkane biosynthesis.The difference between dicotyledonous and monocotyledonous C 3 plants is much more pronounced the regarding the degree to which the leaf water isotope enrichment is transferred into leaf n-alkanes (Gamarra et al., 2016;Kahmen et al., 2013).While dicotyledonous plants show signal transfer rates of 96 % on average (Kahmen et al., 2013), a larger range of between 38 % and 61 % is found for monocotyledonous plants (Gamarra et al., 2016).The latter implies that 39 % to 62 % of the 2 H leaf water enrichment is not recorded by the n-alkanes of grasses.Hence, like for δ 18 O, a correction may be requested to account for grassderived n-alkanes: where δ 2 H l * are the grass-corrected δ 2 H l values.The δ 2 H s values and the non-grassy pollen fraction are defined as in Eq. ( 10).The mass balance correction presented in Eq. ( 11) is based on assumptions that only 50 % of the leaf water enrichment is incorporated by the n-alkanes during biosynthesis in grass leaves.
In summary, the discussion outlined above allows reconstructing δ 2 H/δ 18 O leaf-water (and thus RH results with Eq. 5) for four scenarios (see also Using the Gemündener Maar d l values as input for Eq. ( 5), RH values during daytime and vegetation period (RH dv ) can be calculated (scenario 1 in Table 1).Reconstructed RH dv values range from 32 % to 82 % (Fig. 6b).The error bars covering the Gemündener Maar RH dv record, calculated using pooled d e standard errors ranging from 3.2 ‰ to 44.4 ‰ according Eq. ( 7), result in an RH uncertainty range of 1.7 % to 23.4 %.The RH dv record shows quite large variability with no clear trend during the Allerød and the first half of the Younger Dryas.The late Younger Dryas and the early and the middle Preboreal are characterised by lower RH values.By contrast, the middle Preboreal reveals the most pronounced RH maximum.The mean reconstructed RH dv value is 53 % (mean RH dv upper limit = 45 %; mean RH dv lower limit = 62 %; see Sect.2.5).For comparison, the modern RH dv value (06:00 to 19:00 CET from April to October) from the adjacent meteorological station Nürburg-Barweiler (approx.25 km northeast of Gemündener Maar (GM); hourly data from 1995 to 2015 from Deutscher Wetterdienst, 2016) is 67 % (Fig. 6b).In addition, the range of the reconstructed RH dv values of 50 % agrees well with the modern RH dv variability of 45 %, within a range of 48 % to 93 % (definition and meteorological station details as above).As proposed, in the previous chapter, three correction scenarios can be applied when reconstructing d l and RH dv values in order to account for 2 H and 18 O signal dampening occurring in grasses.
Accordingly, the full correction for grass-derived alkane and sugar biomarkers (scenario 4 in Table 1) results in 0.0 % to 6.3 % (mean 1.8 %) lower RH dv values (RH dv # * in Fig. 6b).This corresponds to d l decreases of 0.0 ‰ to −12.0 ‰ (d l # * in Fig. 6a).Such small changes are still far below the pooled analytical standard errors.When only correcting for the 18 O signal dampening (scenario 3 in Table 1), d l values decrease by 0.0 ‰ to −22.7 ‰, corresponding to RH decreases of 0.0 % to −12.0 % (d l # and RH dv # in Fig. 6a and b, respectively).By contrast, when only correcting for the 2 H signal dampening (scenario 2 in Table 1), this leads to 0.0 ‰ to 10.6 ‰ more positive and 0.0 % to 5.6 % higher RH dv values (d l * and RH dv * in Fig. 6a and b).Overall, these results suggest that the reconstructed RH dv values are not strongly affected by 2 H and 18 O signal dampening of grasses.
We are aware that microclimatic conditions with higher RH values often develop in lower canopy levels of forests (Graham et al., 2014;Parker, 1995).This may result in RH overestimations when applying the coupled δ 2 H n-alkaneδ 18 O sugar paleohygrometer approach.However, most leaf biomass is produced at higher canopy levels, which are exposed to sunlight and free-air RH values.This is in agreement with a study of Zech et al. (2015), who investigated n-alkanes in soils of the tropical montane rainforest of Mt Kilimanjaro.There, n-alkanes reflect δ 2 H leaf-water as calculated from free-air RH rather than as calculated from nearly saturated ground-level RH.
A basic assumption of our coupled δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach is isotopic equilibrium between www.clim-past.net/15/713/2019/Clim.Past, 15, 713-733, 2019 Table 1.Scenarios 1-4 used for reconstructing deuterium (d) excess of leaf water and corresponding RH dv values in order to assess/estimate the effect of variable grass contributions on the reconstructed Gemündener Maar RH record (see also Fig. 6).
Scenario Leaf water reconstructed Equations used for leaf Resulting d excess Relative air humidity from n-alkane or sugar water reconstruction of leaf water as during daytime and biomarkers input for Eq. ( 5) vegetation period according Eq. ( 5) 8) and ( 9) δ 2 H l * /δ 18 O l (8) and ( 9 plant source water and water vapour.In order to test the robustness of this assumption and respective effects on reconstructed RH values, we used data of Jacob and Sonntag (1991), who measured the isotope composition of precipitation and of atmospheric water vapour in Heidelberg, Germany, during the period 1981 to 1989.The mean difference between the annual weighted means of precipitation (≈ plant source water) and the water vapour averaged over the vegetation period (April-October) was therefore calculated.Such derived apparent fractionation (ε ap ) amounts to 18.3 ‰ and 1.57 ‰ on average for 2 H and 18 O, respectively.We used this ε ap in Eq. ( 1) instead of the difference δ a − δ s and recalculated the RH values.This recalculation leads to an average RH change of only −1.7 % (±0.9), which is far below the analytical errors of the d excess of leaf water.
Finally, the stability of the d excess and slope of the LMWL Germany through the past needs to be discussed.According to Stumpp et al. (2014), the long-term d excess of precipitation from 28 sites in Germany does not show pronounced relationships to local climate conditions of the site.All reported values are close to 10 ‰, which indicates that Atlantic air masses are the main moisture source for Germany (e.g.Rozanski et al., 1993).In addition, the d excess of precipitation from the stations Trier and Koblenz,  2010), respectively.RAPiD-15-4P N. pachyderma data from Thornalley et al. (2011).Note that each record is plotted on its own timescale (for planktonic Mg/Ca data, see Thornalley et al., 2009Thornalley et al., , 2010; for 10 Be data on GICC05, see Rasmussen et al., 2006; for 14 C data on IntCal13 calibration curve, see Reimer et al., 2013; for RH dv data on GM age-depth model, see Fig. 2d).AL: Allerød, LST: Laacher See Tephra, YD: Younger Dryas, PB: Preboreal and BO: Boreal.which are close to the Gemündener Maar, reveal rather small variability on a monthly, annual and long-term basis.For Trier monthly averaged d excess values (March to October) range from 5.3 ‰ to 8.7 ‰.Annually weighted mean d excess values range from 1.9 ‰ to 10.6 ‰, and the long-term weighted mean is 6.7 ‰ (±2.2); for Koblenz the d excess values range between 2.1 ‰ and 6.4 ‰ and 1.4 ‰ and 8.7 ‰, and the long-term weighted mean is 4.1 ‰ (±1.8) (derived from IAEA/WMO, 2018).Finally, d excess variability in Greenland and Antarctic ice cores does not exceed 4 ‰ over the timescale relevant here (Masson-Delmotte et al., 2005;Stenni et al., 2010).In addition, paleowater samples from Europe suggest that the d excess of precipitation was rather constant throughout the past 35 000 years, which implies that the principle atmospheric circulation patterns over the European continent did not change substantially (Rozanski, 1985).In summary, the variations in the slope of the LMWL of Germany are assumed to be rather small over longer timescales.
The detailed discussions in the above three sections address numerous uncertainties when using the coupled δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach.Conclusively, the reconstructed RH dv history of the Gemündener Maar seems, however, robust enough to infer reliable paleoclimatic or hydrologic conclusions.

How dry was the Younger Dryas in western Europe?
While it is well known that the Younger Dryas was a cold spell occurring in the Northern Hemisphere during the Late Glacial (Denton et al., 2010;Heiri et al., 2014;Isarin and Bohncke, 1999), there is much less clear evidence concerning moisture supply or availability and RH changes during the Younger Dryas.The Gemündener Maar RH dv record suggests quite some variability but on average moderate RH dv conditions of ∼ 56 % during the end of the Allerød and the first half of the Younger Dryas.This is within the range of modern RH dv values (Fig. 6b).In the second half of the Younger Dryas, a clear RH dv decrease of ∼ 11 % occurred (Fig. 7a).Such a two phasing of the Younger Dryas has been suggested previously based on multiproxy climate data for western Europe (Isarin et al., 1998).In more detail, Isarin et al. (1998) reported a cold and humid first phase being followed by drier and warmer conditions.It is moreover speculated that a shift in the mean sea-ice margin during winter in the North Atlantic Ocean slightly to the north could have caused this two phasing.Reduced cyclonic activity and precipitation thereby primarily affected western Europe because this region was situated at the southern margin of the main storm tracks during the first Younger Dryas period (Isarin et al., 1998).The authors also presented evi-  (Brauer et al., 2008).
In line with this, Brauer et al. (1999) interpreted high biogenic opal contents and Pediastrum remains concentrations during the early Younger Dryas as enhanced aquatic productivity due to an increased nutrient supply caused by soil erosion and the reworking of littoral sediments.The varve formation throughout the second Younger Dryas period (between 12 240 and 11 590 varve a BP) is interpreted to be mainly controlled by snowmelt-driven surface runoff (Brauer et al., 1999).Moreover, the authors speculate if during that time the Meerbach began to drain into the Meerfelder Maar, which could be possibly linked to enhanced precipitation amounts.In summary, the interpretations derived from the Younger Dryas sediments of the Meerfelder Maar by Brauer et al. (2008Brauer et al. ( , 1999) seem neither to be in accordance with the results of Isarin et al. (1998) nor with the established RH dv record of the Gemündener Maar (Fig. 7a).
Recently, Rach et al. (2017) reconstructed RH changes and generally dry Younger Dryas climatic conditions by investigating δ 2 H of terrestrially versus aquatically derived nalkanes (published in Rach et al., 2014) from the Meerfelder Maar archive.At the current state of research, it can only be speculated about the reasons for this discrepancy, with our Gemündener Maar RH record not corroborating an overall dry Younger Dryas.While the uncertainties of the coupled δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach were discussed in detail in the previous sections, in our opinion the most important uncertainties affecting the dual-biomarker approach of Rach et al. (2014Rach et al. ( , 2017) ) are the following.First, lake water is assumed to reflect δ 2 H of precipitation.Indeed, Holzmaar, which seems to be comparable to the Meerfelder Maar at least for the drainage conditions via one creek, shows a difference of 7.4 ‰ in δ 2 H between inflow and lake water (Sachse et al., 2004).This lake water enrichment is likely to have been variable in the past, especially when including the speculation concerning the drainage of the Meerbach during the Younger Dryas (Brauer et al., 1999).Second, n-C 23 is interpreted to be of aquatic origin (from Potamogeton) and used for reconstructing δ 2 H lake-water .However, there is increasing evidence that n-C 23 is also of terrestrial origin (Rao et al., 2014).For instance, Aichner et al. (2018) have recently shown for a lake in Poland that n-C 23 shows a variable mixture of aquatic and terrestrial origin in those Late Glacial and Early Holocene sediments.And birch as a pioneering and one of the dominant tree species during Late Glacial reforestation of central Europe is known to produce considerable amounts of mid-chain n-alkanes (Tarasov et al., 2013).Although they are not included in the latter publica-tion, n-C 23 concentrations of Betula exilis and Betula pendula reached 653 and even 2323 µg g −1 in that study.This is highly relevant, because the biosynthetic fractionation factor of aquatic n-alkanes is much smaller than the one of terrestrial n-alkanes.Minor changes in the contribution of terrestrial vs. aquatic n-alkanes will thus have a considerable impact on the reconstructed δ 2 H n-C 23 record and in turn on reconstructed RH values when applying the dual-biomarker approach.Finally, it may worth acknowledging that Sachse et al. (2004) found no significant correlation for δ 2 H of n-C 23 and lake water and precipitation along a European lake surface transect.
Also recently and also applying the dual-biomarker approach, Muschitiello et al. (2015) studied Younger Dryas lake sediments from Hässeldala Port in southern Sweden.
Here, the authors used δ 2 H of n-C 21 as a proxy for lake water and summer precipitation.The calculated difference between terrestrial and aquatic n-alkane δ 2 H values suggests more humid conditions at the beginning of the Younger Dryas followed by a more or less steady trend towards drier conditions, peaking around 11 700 a BP (Muschitiello et al., 2015).Within age uncertainties, this would be in line with the Gemündener Maar RH dv minimum between ∼ 11 700 and 11 900 a BP.Last but not least, Gázquez et al. (2018) analysed triple oxygen and hydrogen isotopes of gypsum in the southern Pyrenees and thus reconstructed RH changes.Again, more humid conditions are reported for the beginning of the Younger Dryas.
In search of possible drivers or mechanisms for the observed Gemündener Maar RH dv record, we came across the 14 C production and 10 Be flux rates (Fig. 7b), derived from IntCal13 and the Greenland ice cores (GRIP, GISP2), respectively (Muscheler et al., 2014).These records are commonly interpreted in terms of solar activity (and thus insolation) changes (Stuiver and Braziunas, 1988;Vonmoos et al., 2006) and reveal striking similarities with our Gemündener Maar RH dv record.For instance, all three records reveal quite high centennial-scale variability during the Allerød and the first half of the Younger Dryas.Generally low RH dv values during the second half of the Younger Dryas and the Early Preboreal coincide with high solar activity, whereas the pronounced RH dv maximum from 11 260 to 11 050 a BP coincides within age uncertainties with a pronounced solar activity minimum (Fig. 7).We dub this wet period the "Preboreal Humid Phase", which should not be confused with the Preboreal Oscillation (Björck et al., 1997).The Preboreal Oscillation is a short cold event recorded in Greenland ice cores ∼ 11 400 ka (Rasmussen et al., 2007) and led to more arid conditions at least in the Netherlands according to palynological results (Bos et al., 2007;van der Plicht et al., 2004).These pollen records also show the existence of a pronounced humid phase thereafter, thus corroborating the Preboreal Humid Phase.Widespread glacial advances in the Alps are also attributed to the Preboreal Oscillation (Moran et al., 2017).However, given the dating uncertainties they may actually rather reflect increased precipitation during the Preboreal Humid Phase.
It should be emphasised, that the described similarities between the Gemündener Maar RH dv record and the solar activity records do not allow an a priori causality interpretation.It is widely accepted that the Younger Dryas and the Preboreal Oscillation are related to freshwater forcing in the North Atlantic (e.g.Fisher et al., 2002;Murton et al., 2010;Muschitiello et al., 2015).However, the causes and mechanisms responsible for climate and environmental changes during the rest of the Holocene remain vague, and more research including paleoclimate modelling is clearly needed and encouraged to investigate the possible influence of solar activity (Renssen et al., 2007;Rind, 2002).We propose that both the North Atlantic Ocean temperature and solar activity (the latter triggering solar insolation) were the two main drivers for the RH dv variability in central Europe.A key example might be the Preboreal Humid Phase.It can be expected that the North Atlantic Ocean, the main moisture source for central Europe, already had considerably higher temperatures during the Preboreal Humid Phase compared to the Younger Dryas, as indicated by a consistent ∼ 2 • C increase in Mg/Ca temperatures derived from planktonic foraminifera (Globorotalia inflata, Globorotalia bulloides and Neogloboquadrina pachyderma) in a marine sediment core south of Iceland (Fig. 7b, Thornalley et al., 2009, 2010, 2011).This led to an enhanced moisture content of the atmosphere.When these wet air masses were transported onto continental Europe, where low solar insolation inhibited warming up and drying of these air masses, pronounced humid climate conditions were established.

Conclusions
Referring to the underlying research questions and based on the presented results and the outlined discussion (including the cited literature), the following conclusions have to be drawn.
The terrestrial vs. aquatic origin of bulk sedimentary organic matter cannot be determined unambiguously for the Gemündener Maar.This is caused by the bulk proxies (TOC/N, δ 13 C and δ 15 N) not being straightforwardly interpretable.By contrast, the alkane biomarkers with the chainlength n-C 27 and n-C 29 and the sugar biomarker arabinose can be most likely associated with the epicuticular leaf wax layers and the hemicellulose structures of higher terrestrial plants, respectively.Therefore, they are interpreted as originating primarily from leaf material of the Gemündener Maar catchment.
δ 2 H/δ 18 O leaf-water could be reconstructed from δ 2 H n-alkane (n-C 27 and n-C 29 ) and δ 18 O sugar (arabinose) by applying biosynthetic fractionation factors.We acknowledge that the assumption of constant fractionation factors introduces uncertainty as highlighted by the broad literature discussion.A correction for the signal dampening of leaf water 2 H/ 18 O enrichment occurring in grasses is possible but seems negligible in the case of the Gemündener Maar record.
The detailed discussion of possible uncertainties of the applied coupled δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach suggests that robust RH reconstructions are possible for the Gemündener Maar record.The reconstructed RH values refer to daytime and vegetation period (RH dv ).
The established Gemündener Maar RH dv record supports a two phasing of the Younger Dryas with moderate wet conditions at Allerød level during the first half and drier conditions during the second half of the Younger Dryas.Overall, dry climatic conditions characterising the Younger Dryas could not be corroborated.Unexpectedly, the amplitude of RH dv changes during the Early Holocene was more pronounced than during the Younger Dryas and includes a pronounced Preboreal Humid Phase occurring from ∼ 11 260 to 11 050 a BP.We propose North Atlantic Ocean temperature and solar activity (and thus insolation) as the main drivers for Late Glacial-Early Holocene RH changes in central Europe and encourage respective paleoclimate modelling studies in order to validate or falsify our proposition.
Data availability.The data are available in the Supplement, including sampling depth, calculated ages, measured δ 2 H n-alkane and δ 18 O sugar values, reconstructed δ 2 H/δ 18 O leaf-water , d excess of leaf water, RH dv , and δ 2 H/δ 18 O source-water values.
Author contributions.JH, LW and MZ wrote the paper; MZ and RZ acquired financial support.FS was responsible for lake coring and provided the chronology and stratigraphy.LW, JH, TB, MB and IKS carried out laboratory work and did data evaluation.MB, IKS, BG, KR, FS and RZ contributed to the discussion of the data and commented on the paper.
Competing interests.The authors declare that they have no conflict of interest.

Figure 1 .
Figure 1.(a) Location of the Gemündener Maar in the Eifel region in Germany (generated using OpenStreetMap homepage, ©OpenStreetMap contributors, https://www.openstreetmap.org, last access: 15 August 2016).(b) Digital terrain model and drainage system of the immediate surroundings of the Gemündener Maar, with maar borders according to Büchel (1994) representing the size of the crater.In addition, the core position is displayed (GM1; 50 • 10 39.853 N, 6 • 50 12.912E) along with the short core named GMf (not part of this study) marked as ELSA drill sites.Both cores are part of the Eifel Laminated Sediment Archive Project (ELSA project).(c) Photo of Gemündener Maar showing the steep and densely forested catchment (by Michael Zech, reproduced with his permission).

Figure 2 .
Figure 2. (a)Photo of the investigated GM1 core section, with regard to the biomarkers (607 to 694 cm depth), displaying the position of the Laacher See Tephra (LST), varve-counted to 12 880 a BP in the adjacent Meerfelder Maar (cf.Brauer et al., 1999).(b) Defined pollen zones according toBrauer et al. (1999) andLitt et al. (2009).(c) Pollen profiles of pollen groups, which were used for defining the pollen zones.Pollen analysis was carried out by Frank Dreher (Johannes Gutenberg University of Mainz).(d) Age-depth model of the full investigated GM1 section (606 to 727 cm depth) consisting of a 14 C-dated piece of charcoal, the LST and the onsets of the Younger Dryas, Preboreal and Boreal (Holocene).Additionally, the biomarker sampling points are displayed (black points).The error bars of the 14 C age and the LST represent the uncertainty of the calibration (68 % probability range) and the error during of the varve counting (±40 a;Brauer et al., 1999), respectively.

Figure 4 .
Figure 4. (a) Depth profiles of compound-specific stable hydrogen isotope composition of the individual alkanes n-C 27 and n-C 29 and the weighted mean (δ 2 H n-alkane ).(b) Compound-specific stable oxygen composition of the individual sugars arabinose, xylose and fucose (δ 18 O sugar ).Error bars show analytical standard errors; bold lines show three-point moving averages.(c) Depth profile of Poaceae pollen.Additionally, the resampled data points (black points) used for the grass correction procedures (Eqs.10 and 11) are displayed.In addition, the GM1 core picture with the age markers used is displayed.AL: Allerød; LST: Laacher See Tephra; YD: Younger Dryas; PB: Preboreal; BO: Boreal.
identification of bulk organic matter and of the investigated n-alkane and sugar biomarkers

Figure 5 .
Figure 5. Conceptual framework of the coupled δ 2 H n-alkane -δ 18 O sugar paleohygrometer approach displayed as δ 18 O-δ 2 H diagram showing the measured n-alkanes (weighted mean of n-C 27 and n-C 29 ) and sugar (arabinose) biomarkers (black crosses), the reconstructed leaf water (open circles), the global meteoric water line (GMWL, green line), and the local meteoric water lines of Germany (LWML Germany , black line), Trier (LWML Trier , yellow line) and Koblenz (LWML Koblenz , blue line).The black arrows indicate natural processes of evapotranspirative enrichment of leaf water along local evaporation lines (LELs), biosynthetic fractionation during biomarker synthesis and the temperature effect on the source water isotope composition (∼ precipitation).Grey lines indicate the parallel distance between the individual reconstructed evaporative site leaf water points and the LMWL Germany , expressed as d = δ 2 H − 7.72 • δ 18 O.The difference between the d excesses of the leaf water and source water can serve as proxy for mean daytime vegetation period relative humidity (RH dv ; red double arrow).

Figure 6 .
Figure 6.(a) Deuterium (d) excess depth profiles of reconstructed leaf water: d l (black line): no correction for grasses; d l * (light blue line): δ 2 H corrected for grasses; d l # (light red line): δ 18 O corrected for grasses; d l # * (light green line): δ 2 H and δ 18 O corrected for grasses.The error bars of d l values are calculated according to Eq. (7).(b) Reconstructed RH dv records.Modern RH variability during daytime and vegetation period (RH dv ) is displayed as a box plot derived from the adjacent meteorological station Nürburg-Barweiler, using monthly means from April to October between 06:00 and 19:00 CET (based on hourly data from 1995 to 2015; Deutscher Wetterdienst, 2016).The numbers within the box plot represent the maximum, median and minimum values.(c) Depth profile of Poaceae pollen.Additionally, the resampled data points (black points) are displayed.The GM1 core picture with the used age markers are displayed on the left.AL: Allerød; LST: Laacher See Tephra; YD: Younger Dryas; PB: Preboreal; BO: Boreal.
www.clim-past.net/15/713/2019/Clim.Past, 15, 713-733, 2019 dence for the strengthening of the westerly winds in western Europe as consequence of northward-shifted North Atlantic Ocean sea-ice margin during the late Younger Dryas period.This contradicts, however, with the interpretation of the Meerfelder Maar sedimentary record.Here, the thicker varves during the early Younger Dryas (between 12 680 and 12 240 varve a BP) are used along with geochemical results as indicator of stronger winter winds