Holocene sea level and environmental change at the southern Cape – an 8.5 kyr multi-proxy paleoclimate record from Lake Voëlvlei, South Africa

. South Africa is a key region to reconstruct and understand past changes in atmospheric circulation, i.e. temperate westerlies and tropical easterlies. However, due to the scarcity of natural archives, South Africa’s environmental evolution during the late Quaternary remains highly debated. Many available sediment archives are peri-coastal lakes and wetlands; however, the paleoenvironmental signals in these archives are often overprinted by sea-level changes during the Holocene. This study presents a new record from the coastal wetland Voëlvlei, which is situated in the year-round rainfall zone of South Africa on the southern Cape coast. It presents an ideal sedimentary archive to investigate both sea level and environmental changes. A 13 m long sediment core was retrieved and analysed using a multi-proxy approach. The chronology reveals a basal age of 8440 + 200 / − 250 cal BP. Paleoecological and elemental analyses indicate marine in-cursions from ca. 8440 to ca. 7000 cal BP with a salinity opti-mum occurring at 7090 + 170 / − 200 cal BP. At ca. 6000 cal BP, the basin of Voëlvlei was in-ﬁlled with sediment resulting in an intermittent (sporadically desiccated) freshwater lake similar to present.

P. Strobel et al.: Holocene sea level and climate variability at Lake Voëlvlei more reliable, robust paleoclimate reconstructions using direct hydrological proxies are necessary but have been rare in South Africa up to now (Haberzettl et al., 2014). Southern Africa's past and present climate has been driven by complex interactions between two major oceanic and atmospheric circulation systems, i.e. the Benguela and Agulhas currents, and the westerlies and easterlies (Tyson and Preston-Whyte, 2000) (Fig. 1a). Today, three major rainfall zones occur in South Africa. While the eastern and central parts of the country receive most rainfall (> 66 %) from tropical moisture-bearing atmospheric circulation systems during austral summer (summer rainfall zone, SRZ), a narrow belt along the west coast receives most rainfall (> 66 %) from temperate westerlies during the austral winter (winter rainfall zone, WRZ) (Fig. 1a). An intermediary area between the SRZ and WRZ receives rainfall from both systems throughout the year (year-round rainfall zone, YRZ) ( Fig. 1a) (Engelbrecht et al., 2015;Scott and Lee-Thorp, 2004), and this includes the southern Cape coast, which is the focus area for this study.
The YRZ has been the focus of most paleoenvironmental and associated paleoclimatic research in this region. Here, the southern Cape coast, in particular the Wilderness area with its numerous coastal lakes including Bo Langvlei (du Plessis et al., 2020), Eilandvlei (Kirsten et al., 2018a, b;Quick et al., 2018;Reinwarth et al., 2013;Wündsch et al., , 2016b, Groenvlei (Martin, 1959(Martin, , 1968Wündsch et al., 2016a) and Swartvlei (Birch et al., 1978;, has yielded multiple paleoenvironmental records (Fig. 1b). These coastal lakes have formed between large coastal dune cordons that lie parallel to the coast. However, the terrestrial climate signals in these coastal archives are often overprinted by marine water intrusions induced by relative sea level change during the Holocene (Martin, 1959(Martin, , 1968Reinwarth et al., 2013;Wündsch et al., 2018Wündsch et al., , 2016a or anthropogenic impacts, e.g. farming, water abstraction and dredging operations . Further paleoenvironmental information from this area is available from, for example, peatlands (Quick et al., 2016;Strobel et al., 2019), speleothems (Braun et al., 2018(Braun et al., , 2020Talma and Vogel, 1992), rock hyrax middens (Chase et al., 2019(Chase et al., , 2020 and marine sediments (Hahn et al., 2017) (Fig. 1a). However, the climate evolution of South Africa is still debated, reflecting potential spatial variability in climate drivers at the regional scale, and the application of different methodological approaches in various studies Strobel et al., 2019). Consequently, our understanding of environmental dynamics and changing interactions between tropical and temperate climate systems affecting the YRZ is limited.
Compound-specific stable isotope analyses of hydrogen and carbon isotopes of long-chain n-alkanes (≥ C 25 ; δ 2 H n-alkane , δ 13 C n-alkane ) are valuable proxies that complement established methodological approaches related to paleoenvironmental and paleohydrological changes in sediment  Talma and Vogel, 1992;Eilandvlei: EV, Quick et al., 2018;Wündsch et al., 2018;Groenvlei: GV, Wündsch et al., 2016a;Seweweekspoort: SWP, Chase et al., 2017; Verlorenvlei: VER, Kirsten et al., 2020;and Vankervelsvlei: VVV, Strobel et al., 2019). Additionally, the circumpolar westerlies, the tropical easterlies, the Agulhas Current (AC) and the Benguela Current (BC) are depicted. (c) Voëlvlei, its catchment, the coring position as well as a paleo-sea level highstand (+5 m) and the recent shoreline (data sources -rainfall seasonality: WorldClim 2 dataset, Fick and Hijmans, 2017; circulation systems after Chase and Meadows, 2007; DEM: SRTM 1 arcsec; ∼ 30 m, Jarvis et al., 2008). archives, e.g. grain size analyses, geochemistry and palynology. Long-chain n-alkanes are leaf waxes produced by higher terrestrial plants and serve as valuable biomarkers as they remain well preserved in soils and sediments over millennia because of their low water solubility and high resistance to degradation (Eglinton and Eglinton, 2008;Sachse et al., 2012;Sessions, 2016). In South Africa, the δ 2 H n-alkane signal shows potential to reconstruct the isotopic signal of precipitation and thus directly refers to the precipitation source (Herrmann et al., 2017;Strobel et al., 2020). However, δ 2 H n-alkane has rarely been used in terrestrial archives at the southern Cape coast (Strobel et al., 2019).
δ 13 C n-alkane is a suitable proxy to infer past changes in the vegetation composition (e.g. Diefendorf and Freimuth, 2017) as well as variations in plant water use efficiency and thus drought stress (Diefendorf and Freimuth, 2017;Struck et al., 2020). Therefore, using the climatic information from both leaf wax isotopes enables climate reconstructions based on the isotopic signal of precipitation, which can reflect local water availability. Further, leaf wax-derived n-alkanes can be used as a chronological marker since they can be dated using 14 C analyses (Bliedtner et al., , 2018Douglas et al., 2014;Gierga et al., 2016;Haas et al., 2017).
Here we present a sediment record from Voëlvlei, today an intermittent lake (sporadically desiccated), to reconstruct past sea level and environmental changes at the southern Cape coast of South Africa. A multi-proxy approach has been applied to the sediments, comprising compound-specific stable isotope analyses of leaf waxes as well as sedimentological assays of inorganic and organic elemental compositions, and paleoecological analyses on fossil remains. Specifically we aim to i. establish a robust chronology based on diverse dating approaches on different sediment compounds, ii. disentangle marine and climate influences during the development of Voëlvlei, and iii. infer variations in local moisture availability and the source of precipitation.

Site description
Voëlvlei is situated ca. 40 km east of Still Bay and 30 km west of Mossel Bay, at an elevation of 5 m above present sea level (a.s.l.) ca. 10 km inland of the Indian Ocean coast. Today, Voëlvlei has an area of 3.8 km 2 (max. length 4.2 km; max. width 0.7 km) and the catchment has an area of 165 km 2 (database: SRTM 1 arcsec; Jarvis et al., 2008) (Fig. 1c). The catchment altitudes vary between 5 and 333 m a.s.l.
The potential natural vegetation (i.e. that which would be expected to occur in the absence of human impact) consists of variations of fynbos, and only small areas along the main drainage system would be covered by Albany thicket (Mucina and Rutherford, 2006). Today, large areas of the catchment are used for agriculture. The steep slopes of the drainage system are mainly unmanaged and covered by plant communities of the fynbos and Albany thicket.
Mean annual precipitation at the study site is 450 mm yr −1 (Fick and Hijmans, 2017), and rainfall is almost equally distributed throughout the year. Winter precipitation is linked to the temperate westerlies related to the Atlantic ocean as moisture source, and summer precipitation is associated with the tropical easterlies and the Indian Ocean as moisture source (Engelbrecht and Landman, 2016). Moreover, orographic rainfall occurs from local sources due to onshore flows related to ridging anticyclones (Weldon and Reason, 2014). The isotopic composition of precipitation (δ 2 H p ) is 2 H-depleted during winter and 2 H-enriched during summer periods, with a modelled annual mean of −13±1 ‰ (Table 1) (Bowen, 2018;Bowen et al., 2005;Braun et al., 2017;Harris et al., 2010). Mean annual temperature is 17.6 • C and slightly higher temperatures during summer (22 • C) lead to semi-arid climatic conditions at the study site today (Fick and Hijmans, 2017).

Material and methods
For this study, the 13 m long sediment core (VOV16) was retrieved from Voëlvlei (34.259 • S; 21.826 • E; Fig. 1c) in 2016 using a percussion hammer coring system (inner core diameter 5 cm) and transported to the Physical Geography laboratory of the Friedrich Schiller University Jena where it was stored under dark and cool (ca. 4 • C) conditions until processing. Cores were opened, cleaned and photo-documented, and their sedimentological properties and sediment colour were described in detail following the standard protocols of the Physical Geography laboratory at Friedrich Schiller University Jena.

Chronology
3.1.1 Radiocarbon dating of macro particle, bulk TOC and n-alkane samples The chronology of the sediment record is based on 14 C ages from 1 organic plant-macro particle, 3 charcoal samples, 15 bulk organic samples and 7 n-alkane samples (compoundclass). In total, 12 of the bulk organic samples and the organic plant-macro particle were analysed using accelerator mass spectrometry (AMS) at the Poznan Radiocarbon Laboratory, Poland. Three bulk organic, three charcoal and seven n-alkane samples (see Sect. 3.7 for sample extraction prior to measurement) were analysed with the Mini Carbon Dating System (MICADAS) AMS coupled to an element analyser (Ruff et al., 2010;Salazar et al., 2015;Szidat et al., 2014) at the LARA AMS Laboratory, University of Bern, Switzerland. 14 C results from the LARA AMS were reported as F 14 C and corrected for cross-contamination and constant contamination after Salazar et al. (2015).  (Bowen, 2018;Bowen et al., 2005).

OSL dating
A split (half) section of the core was sub-sampled under red light conditions at the University of Leicester luminescence dating laboratory. The upper 5-6 mm of the sediment surface was removed, and the core section was sampled over a depth range of 70 mm. Sediment within 6-7 mm of the core tube inner surface was not sampled but was used for an estimation of sample water content. This (dry) material as well as material from the upper surface was homogenized and used for dose rate analysis. The sediment for equivalent dose analysis was soaked in sodium hexametaphosphate and then wet sieved. The core sediments yielded very limited amounts of sandsized material, most of which was less than 100 µm. This necessitated the use of the fine sand range 55-90 µm for equivalent dose analysis. This material was prepared using standard methods (e.g. Aitken, 1985;Preusser et al., 2008;Wintle, 1997). This involved treatment with dilute (10 %) hydrochloric acid (to remove carbonates) and (32 %) hydrogen peroxide (to remove organics). The sample was then dried and density separated to isolate the < 2.7 g cm −3 and > 2.58 g cm −3 (quartz) fraction, before etching for 45 min in 48 % hydrofluoric acid, washing in dilute hydrochloric acid (HF, to remove fluorides) and dry sieving (55 µm). The HF etch removed any remaining K feldspars and also served to etch the (alpha irradiated) outer surface of the quartz grains. Dose rates were determined using the same core material used to estimate water content via inductively coupled plasma mass spectrometry (ICP-MS) for U and Th and ICP-OES for K analyses at the University of Leicester. The concentrations of U, Th and K were converted to annual dose rates following Guérin et al. (2011) with corrections for grain size (Mejdahl, 1979), water content (Aitken, 1985) and HF etching (Bell, 1979). Cosmic dose rates were determined using the reported sample depth following Prescott and Hutton (1994) with a 5 % relative uncertainty included. Final age uncertainties incorporate 3 % relative uncertainties for the dose rate conversion factors, grain size attenuation factor, water attenuation and HF etching, propagated in quadrature. HF etching is assumed to have entirely removed the α-irradiated outer portion of the quartz grains. It was assumed that the as-measured water content was appropriate with a 3 % (absolute content) uncertainty propagated to the final dose rate uncertainty. In the absence of robust direct evidence for a substantially different water content in the past, we have utilized the modern sample water content. The 3 % absolute uncertainty associated with this will account for some fluctuations through time. For reference, a 10 % (absolute) change in water content results in an age difference of approximately 700 years.
All luminescence measurements were performed on a Risø DA20 TL/OSL reader. Stimulation (40 s at 125 • C) was provided by blue LEDs (stimulation wavelength 470 nm) with OSL signals detected with an EMI 9235QA photomultiplier tube via a U-340 detection filter. Laboratory irradiations were delivered by a 90 Sr beta source with a dose rate (at the time of measurement) of ca. 7.58 Gy min −1 . Sample equivalent doses (D e ) were determined using the singlealiquot regeneration (SAR) protocol (Murray andWintle, 2000, 2003;Wintle and Murray, 2006). All single-aliquot measurements were carried out on small (1 to 2 mm) aliquots, which given the grain size fraction analysed means there are likely > 1000 grains per aliquot (e.g. Duller, 2008).
A dose recovery preheat experiment was used to assess the suitability of the SAR protocol in general, and the most appropriate preheating conditions. The overall dose recovery ratio across all preheating temperatures (160-260 • C) was 0.97 ± 0.03 (n = 22; zero overdispersion), with a ratio of 1.00 ± 0.02 (n = 3) for the chosen preheat temperature combination of 240 • C for 10 s (natural regeneration points) and a 220 • C cut heat for the test dose measurements. All SAR analyses comprised a 7-regeneration point sequence, which included a repeated (recycling) regeneration dose point, an IR depletion regeneration dose point to check for K feldspar contamination (Duller, 2003) and a zero-dose point. The 4 (unique) point dose response curve was generated using the initial 0.64 s of stimulation, with a background signal from the last 8 s. Analyses were carried out in the "Analyst" software. Dose response curves were fitted with saturating exponential fits, with D e uncertainties incorporating counting statistics, curve fitting uncertainties and a 1 % systematic uncertainty (Duller, 2007) (all calculated within the Risø Analyst software). The uncertainty in the final D e estimate also incorporates an additional beta source calibration uncertainty (3 %).

Age modelling
For the final age-depth modelling, all aforementioned dates were used. Bulk organic and n-alkane 14 C ages from the terrestrial part of the sediment core were calibrated with the SHcal20 data calibration curve (Hogg et al., 2020), whereas 14 C ages from the marine part of the record were calibrated with the Marine20 calibration curve  and additionally corrected for a marine reservoir effect using a R of 134 ± 38 14 C years as previously reported by Wünd-sch et al. (2016b) (Table 1). For the compound-class n-alkane samples, the SHCal20 calibration curve (Hogg et al., 2020) was applied due to a predominance of terrestrial synthesized long-chain n-alkanes (C 29 , C 31 , C 33 ) in the samples (see Sect. 4.3; Table 1). The terrestrial organic plant-macro particle and charcoal found in the marine part of the record were calibrated with the SHCal20 curve (Hogg et al., 2020) (Table 1). All calibrations were done with the online version of the Calib 8.2 software (Stuiver et al., 2020). The final agedepth profile was modelled with the R software package Bacon 2.4.3 (Blaauw and Christen, 2011), using the same calibration datasets. In the following, ages are reported as median ages including the upper and lower limit of the 95 % confidence interval.

Grain size analyses
For grain size measurement 50 mg sample aliquots were treated with H 2 O 2 (10 %) to remove organic matter, HCl (10 %) for the destruction of carbonates and Na 4 P 2 O 7 /Na 2 CO 3 for dispersion. The grain size distribution of each sample was determined with a Laser Particle Sizer (FRITSCH ANALYSETTE 22 Microtec, FRITSCH, Germany) and a wet dispersion unit at the Physical Geography laboratory of the University of Greifswald. Before the measurement, samples were treated with ultrasound for 60 s and subsequently measured in duplicate. The grain size distribution is calculated in 99 classes between 0.08-2000 µm. The mean and median grain size as well as fractions of clay, silt and sand of each sample were calculated.

Paleontological analyses
A selection of 23 sediment samples representing all lithological units and focusing on the assumed transition from an estuarine to a freshwater environment were processed and analysed micropaleontologically. About 3-7 mL of sediment from 1 cm thick sediment slices of the core was washed with tap water through stacked sieves of 63 and 200 µm mesh size.
After drying the sieve residues on a heating plate at ca. 50 • C, all microfossils were picked under a low-power stereomicroscope and transferred to microfossil slides for later identification and counting. Microfossils, fragments of macrofossils and charcoal were documented semi-quantitatively as rare (1-2 specimens), common (> 2), abundant (> 10) and very abundant (> 100). Identification relies on Benson and Maddocks (1964), Martens et al. (1996), and Fürstenberg et al. (2017) for Ostracoda and on Schmitt-Sinns (2008) and Fürstenberg et al. (2017) for Foraminifera. Additional macrofossils found during the lithological description, sampling of the core and processing the microfossils were identified, counted and used for paleoecological interpretation. These were all snails with identification relying mainly on Branch et al. (2010) and the mollusc section of the World Register of Marine Species (http://www.marinespecies.org/, last access: 8 July 2020). Bivalves could not be identified because they occurred only as fragments. Paleoecological information was drawn from the papers listed above for identification plus Murray (2006) and Kirsten et al. (2018a). All paleontological material will be stored at the South African Museum of Natural History in Cape Town.

Elemental analyses
XRF data (Al, Br, Zr) were collected every 1 cm downcore using two generator settings (30 kV, 1 mA, 15 s; 15 kV, 0.2 mA, 15 s) for detection of different elemental groups with the XRF Core Scanner II (AVAATECH Serial No. 2) at MARUM at the University of Bremen. The split core surface was covered with a 4 µm thin SPEXCerti Prep Ultralene1 foil to avoid contamination of the XRF measurement unit and desiccation of the sediment. The data reported here were acquired by a Canberra X-PIPS Silicon Drift Detector (SDD; Model SXD 15C-150-500) with 150 eV X-ray resolution, the Canberra Digital Spectrum Analyzer DAS 1000 and an Oxford Instruments 50 W XTF5011 X-Ray tube with rhodium (Rh) target material. Raw data spectra were processed by the analysis of X-ray spectra by iterative least square software (WIN AXIL) package from Canberra Eurisys. The data were normalized by elemental Zr counts, i.e. Al/Zr and Br/Zr, and plotted as log ratios, primarily to eliminate sediment matrix errors (water content, surface roughness and grain size variations) (Weltje and Tjallingii, 2008). Moreover, 150 sample aliquots at 8 cm intervals were freeze dried (−53 • C, for > 48 h), ground and sieved to a size < 40 µm. Aluminium (Al), calcium (Ca) and sodium (Na) concentrations were measured with an ICP-OES 725-ES (VARIAN, USA) at the Physical Geography laboratory of the Friedrich Schiller University Jena. A total of 0.2 g of the samples was processed using a microwave-assisted modified aqua regia digestion of 2 mL HCl (32 %) and 4 mL HNO 3 (65 %). Error estimates were based on triple measurements of three samples (VOV_243, VOV_651, VOV_1186; numbers indicate sample depth) (relative error: Al 3.7 %; Ca: 5.4 %; Na: 5.3 %). Samples of the reference material LGC6 187 (river sediment) were measured as well to calculate the relative analytical error which was 1.1 % for Al and Ca, and 10 % for Na.
Total nitrogen (TN) and total carbon (TC) were analysed with a CNS analyser (EuroVector EA 3000, HEKAtech GmbH, Germany) at the Physical Geography laboratory of the University of Greifswald. Concentrations of total organic carbon (TOC) were determined with the same device after treatment with 3 % and 20 % HCl at 80 • C to remove carbonates. Error estimates were based on triple measurements of 30 samples (mean relative error: N: 8.9 %; TOC: 2.9 %). Total inorganic carbon (TIC) was calculated as difference between TC [%] and TOC [%].

Biogenic silica (BiSi)
Biogenic silica (BiSi) concentrations were determined following Ohlendorf and Sturm (2008). BiSi and Al concentrations were measured using an ICP-OES 725-ES (VAR-IAN, USA) at the Physical Geography laboratory of the Friedrich Schiller University Jena. A correction factor of one (BiSi : Al; 1 : 1) was applied to account for dissolving aluminosilicates.

Powder X-ray diffraction (XRD)
The identification of the mineral composition of 50 powdered samples (e.g. Pecharsky and Zavalij, 2009) from representative core sections was carried out using an X-ray diffractometer (D8-Discover, Bruker AXS) equipped with a CuKα X-ray tube and a gas proportional counter (HI-STAR area detector, Bruker AXS) at Friedrich Schiller University Jena. The qualitative analyses and interpretation of the diffractograms was conducted at Trier University using Bruker DIFFRAC plus 5.0 software. The occurrence of specific mineral phases (low, medium, high) was roughly estimated based on XRD peak intensities in conjunction with the elaborated geochemistry of the investigated sediment sections. For all measurements the quartz peak at 3.342 Å was accepted as internal standard.

Biomarker analyses of n-alkane distributions and stable isotope compositions
Total lipids of the sediment samples (14.5-31.4 g) were extracted with 40 mL dichloromethane (DCM) and methanol (MeOH) (9/1, v/v) using an ultrasonic bath over three 15 min cycles. The total lipid extract was separated by solid phase extraction using aminopropyl silica gel (Supelco, 45 µm) as the stationary phase. The n-alkanes were eluted with 4 mL hexane and further purified over coupled silvernitrate (AgNO 3− ) -zeolite (Geokleen) pipette columns. The n-alkanes trapped in the zeolite were subsequently dissolved in HF and recovered by liquid-liquid extraction using n-hexane. An Agilent 7890 gas chromatograph equipped with an Agilent HP5MS column (30 m, 320 µm, 0.25 µm film thickness) and a flame ionization detector (GC-FID) was used for identification and quantification of the n-alkanes, relative to external n-alkane standards (n-alkane mix n-C 21n-C 40 , Supelco).
n-Alkane concentrations were calculated as the sum of C 25 to C 35 and are given in µg g −1 dry weight. Odd-overeven predominance (OEP) values (Eq. 1) were determined following Hoefs et al. (2002). Low values (< 5) indicate an enhanced state of degradation Zech et al., 2010). The average chain length (ACL) (Eq. 2) was calculated from the odd-numbered n-alkanes (Poynter et al., 1989). OEP = nC 27 + nC 29 + nC 31 + nC 33 nC 26 + nC 28 + nC 30 + nC 32 (1) ACL = 27 · nC 27 + 29 · nC 29 + 31 · nC 31 + 33 · nC 33 nC 27 + nC 29 + nC 31 + nC 33 (2) Compound-specific stable hydrogen isotope analyses of the C 31 and C 33 n-alkanes were carried out on an Iso-Prime vision IRMS, coupled to an Agilent 7890A GC via a GC5 pyrolysis or combustion interface operating in pyrolysis modus with a MaxChrome and silver wool packed reactor at 1050 • C. Samples were injected with a split-splitless injector. The GC was equipped with 30 m fused silica column (HP5-MS, 0.32 mm, 0.25 µm). Each sample was analysed in triplicate, except for single measurements of three samples (VOV_915, VOV_128, VOV_111; numbers indicate sediment depth) due to insufficient compound abundance. δ 2 H n-alkane was measured against calibrated H 2 reference gas and all values are reported in ‰ against VSMOW. 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 in duplicate every nine runs. All measurements were corrected for drift and amount dependency, relative to the standard values in each sequence. Triplicates for the C 31 and C 33 alkanes had a standard deviation of < 4.0 ‰; the analytical error for the standards was < 1.7 ‰ (n = 68). The H + 3 factor was checked every 2 d and stayed stable at 4.40 ± 0.03 during measurements.
Compound-specific stable carbon isotope analyses of C 31 and C 33 n-alkanes were carried out on an IsoPrime vision IRMS, coupled to an Agilent 7890A GC via a GC5 pyrolysis or combustion interface operating in combustion modus with a CuO and silver wool packed reactor at 850 • C. Samples were injected with a split-splitless injector. The GC was equipped with 30 m fused silica column (HP5-MS, 0.32 mm, 0.25 µm). δ 13 C n-alkane values were calibrated against CO 2 reference gas of known isotopic composition and all carbon isotope values are given in ‰ against VPDB. Triplicate injections were conducted for each sample and measurement accuracy was controlled in the same way as for the δ 2 H n-alkane analyses. Triplicates for the C 31 and C 33 alkanes had a standard deviation of < 0.2 ‰; the analytical error for the standards was < 0.2 ‰ (n = 83).

Lithology and chronology
The sediment sequence is 13 m long and consists of three lithological units defined by differences in sediment colour and grain size (Fig. 2). Unit A from 13 to 4.51 m depth consists of dark greyish material and is dominated by silt. Within this unit, several thin greyish clayey and yellowish sandy layers are present. Unit B from 4.51 to 0.78 m sediment depth consists of brownish to reddish silty material. In the lower parts of Unit B thin greyish clayey layers are observed. Unit C from 0.78 m sediment depth to the top is dominated by dark brown silty material. Moreover, two layers (0.53-0.44; 0.31-0 m sediment depth) can be separated through changes in colour and granulometric structure in Unit C (Fig. 2).
The chronology reveals a basal age of 8440 +200 / −250 cal BP (Fig. 2, Table 2) and considering a 95 % confidence interval 76 % of the ages overlap with the age-depth model. In Unit A, 14 C ages from the bulk samples are stratigraphically consistent, except for two samples at 11.19 and 5.66 m sediment depth, which are too young for their stratigraphic position. The 14 C age ranges of a terrestrial organic plant-macro particle and a reservoir-corrected bulk sediment sample at 9.63 and 9.61 m overlap. At 4.57 m sediment depth, two charcoal samples and a reservoir-corrected bulk sediment sample distinctly overlap as well. Above the bulk sediment sample at 4.57 m sediment depth, four bulk sediment samples are too old for their stratigraphic position and are not in stratigraphic order (Units B and C). In contrast, 14 C ages obtained from the n-alkane samples are distinctly younger and are in stratigraphic order in Units B and C (Table 2; Fig. 2).
Consistent with the good dose recovery results the OSL sample produced a bright, rapidly decaying quartz OSL signal (reduced to ca. 2 % of the initial signal within 1 s). None of the analysed aliquots exhibited recycling ratios or IR ratios beyond/below (respectively) 10 % of unity or recuperation (zero dose signal) > 1 % of the natural signal. The equivalent dose distribution however is overdispersed (40 ± 5 %; Table 3), even after removal of one very high (125 Gy) outlier, and despite the considerable signal averaging that is likely given the use of the fine sand fraction for analysis. In the context of the analysed core, its likely antiquity, and the otherwise excellent performance of the sample implied by the dose recovery experiment and internal checks within the SAR protocol, a parsimonious explanation of this broad distribution is the presence of unbleached or incompletely bleached grains. Indeed, the age obtained from the central age model equivalent dose estimate is 10 ± 1 ka, which is implausible given the stratigraphic position of the sample (Fig. 2). Application of a three-component minimum age model (MAM) (Galbraith et al., 1999) administered in the R package "Luminescence" (Burrow, 2019) produces a D e estimate of 20.3 ± 2.4 Gy and an age of 6.8 ± 0.8 ka (Table 3). While caution is required when applying this approach to multi-grain aliquots (especially those with considerable signal averaging as here), the minimum age estimate is much more concordant with the radiocarbon ages from this section of the core and the broader age-depth model for the whole sequence (Fig. 2). In addition, the MAM-derived OSL age is in stratigraphic order with the 14 C ages derived from the n-alkane and charcoal samples (Fig. 2).

Geochemical and paleontological analyses
Significant correlations (Pearson's r; α < 0.05) were obtained for the log-normalized XRF scanning data for Br/Zr and Al/Zr compared to the quantitative elemental contents of Na and Al, respectively (r: Log(Br/Zr)/Na = 0.82; r: Log(Al/Zr)/Al = 0.71; Fig. 3). Na concentrations range from 2090±110 to 10 390±550 ppm over the sequence. Both Na contents and Br/Zr ratios show highest values in Unit A (13-4.51 m sediment depth). Several minima in the Br/Zr ratios derived from XRF-scanning data are not visible in the Na concentrations due to the higher sample resolution of the XRF data (Fig. 3). Both Na concentration and Br/Zr ratios distinctly decrease above 4.51 m sediment depth, but an increase is noticeable in Unit C (< 0.78 m sediment depth) (Fig. 3). Al concentrations range from 39 300 ± 1500 to 82 900 ± 3100 ppm, and both Al contents and Al/Zr ratios show high values over large parts of the sediment record. As mentioned previously, Al/Zr ratios show minima in Unit A, which are not visible in the Al concentrations due to the higher sample resolution of the XRF-scanning data (Fig. 3). Noticeable are distinct low concentrations of Al and Al/Zr ratios from 4.51 to 2.91 m sediment depth and at 1.82 m sediment depth, with an increasing trend to the top of the record (Fig. 3).
TOC and N range from 0.16 % to 3.94 % and 0.03 % to 0.48 % respectively; they show high values in the lower parts of the core (> 4.51 m sediment depth) and distinctly decrease above 4.51 m sediment depth and increase thereafter again (Fig. 3). The C/N molar ratio ranges between 1.9 and 13.6 over the whole sequence with high values between 5.15 and 4.27 m sediment depth (Fig. 3). TIC ranges from 0 % to 0.86 % with noticeable peaks between 6.91 and 6.51, and at 4.83, 2.59 and 1.74 m sediment depth (Fig. 3). Ca concentrations show a similar pattern to TIC and range from 2830±150 to 28 900±1600 ppm (Fig. 3). BiSi concentrations range from 0.8±0.02 % to 2.8±0.08 % with low values in the lower parts of the record (> 6.99 m sediment depth), except for high values between 10.98 and 10.10 m sediment depth (1.6 ± 0.05 % to 2.2 ± 0.07 %) (Fig. 3). They show high values for the upper part of the record < 6.91 m sediment depth (Fig. 3).
Apart from the selected mineral components of the XRD measurements displayed in Fig. 3, all samples have a predominance of quartz, feldspars and micas. Calcite follows the patterns of TIC and Ca and shows highest abundances between 11.46 and 10.18 as well as 3.87 and 2.19 m sediment depth (Fig. 3). The other sections have only medium to low proportions of calcite (Fig. 3). The presence of gypsum is restricted to depths between 9.62 and 9.14 m as well as 5.47 and 5.23 m (Fig. 3). However, an absence of calcite and/or gypsum at the distinct peaks of TIC and Ca (e.g. 2.59 and 1.74 m sediment depth) is likely due to the lower sample resolution of the XRD measurements compared to the elemental data (Fig. 3). High proportions of halite are present in   (Zr), aluminium (Al) and Log 10 ratio of Al and Zr, contents of total organic carbon (TOC), total nitrogen (N), molar C/N ratio, total inorganic carbon (TIC), calcium (Ca), and biogenic silica (BiSi) derived from the sediment core VOV16 from Voëlvlei. The presence of selected mineral components is also depicted. large parts of the sediment record (> 2.19 m sediment depth; Fig. 3).
n-Alkane concentrations range from 0.22 to 5.23 µg g −1 with high values in Unit A and distinctly lower concentrations in Unit B and C (Fig. 4). All samples show a distinct odd-over-even predominance (4.0-15.5), and the ACL ranges from 29.9 to 31.4. Consequently, C 29 , C 31 and C 33 are the predominant n-alkane chain-lengths in all samples (Fig. 4). δ 13 C n-alkane values range from −28.10 ± 0.16 ‰ to −22.72 ± 0.06 ‰ with more negative values at the bottom of Unit A (13.0-12.58 m composite depth), followed by less negative values (12.58-8.85 m composite depth), which show a decreasing trend upwards (Fig. 4). Notably are three higher δ 13 C n-alkane values in the upper parts of Unit A between 6.35 and 4.51 m sediment depth (Fig. 4). δ 2 H n-alkane values range from −154.4 ± 1.3 ‰ to −129.6 ± 1.1 ‰ with Table 2. Conventional radiocarbon ages as well as 2σ calibrated age ranges and median calibrated ages (Calib 8.2) (Stuiver et al., 2020) using the SHCal20 and Marine20 calibration curve Hogg et al., 2020) of dated bulk sediment, n-alkanes and organic macro particle (OMP) samples from the VOV16 record. Samples from the marine and brackish parts of the record were reservoir corrected using a R of 134±38 14 C years after Wündsch et al. (2016b). 1: SHCal20; 2: Marine20 and reservoir corrected using a R of 134±38 14 C years after Wündsch et al. (2016b  less negative values at the bottom followed by a trend to more negative values in Unit A. Unit B shows less negative δ 2 H n-alkane values at the bottom followed by more negative δ 2 H n-alkane values at the top and in Unit C (Fig. 4). Paleontological analyses of the 33 macrofossil samples and 23 microfossil samples revealed 56 snail specimens representing six species. In descending order of abundance these are Turritella capensis (Krauss, 1848), which is clearly dom-inant, the rare Assiminea globulus (Conolly, 1939), Hydrobia sp., Natica tecta (Anton, 1839), Nassarius kraussianus (Dunker, 1846) and a fragment of an unidentified gastropod species (Fig. 5). The microfossil associations show a higher diversity. The dominating foraminifer taxon is Ammonia parkinsoniana (d'Orbigny, 1839), followed by Quinqueloculina sp., and the rare Trochammina inflata (Montagu, 1808) and Haynesina sp. All other foraminifer taxa are represented by only one test each (Bolivina sp., Cribroelphidium articulatum (d'Orbigny, 1839), Spirillina sp. and an unidentified trochamminid). The most abundant ostracod species are Sulcostocythere knysnaensis (Benson and Maddocks, 1964) and juveniles of Loxoconcha parameridionalis? (Benson and Maddocks, 1964); rare ostracods are a myodocopid specimen and juvenile Aglaiella valves as well as the freshwater species Sarscypridopsis aculeata (Costa, 1847), Cyprilla humilis (Sars, 1924) and an unidentifiable fragment of a larger species. Microfossils identified on a group level only are ephippia of cladocerans, gyrogonites of charophytes, Plumatella-like bryozoan statoblasts, mollusc and insect fragments, fish bone remains, fruits, and seeds as well as unidentifiable plant remains. Charcoal was found in large quantities in many samples.
Unit A1 shows a high diversity and contains most of the marine-brackish snails and high numbers of foraminifera. The ostracod fauna is dominated by brackish water taxa. Shell and fish bone fragments are abundant. Unit A2 looks similar to Unit A1 but is generally more variable in abundances and diversity. Macrofossils, i.e. marine-brackish snails, become rare. Salt marsh foraminifera, fruits and seeds occur for the first time. Units B and C lack marine-brackish snails, brackish ostracods and salt marsh foraminifera, and foraminifera in general are documented with a single test at the base of the zone only. In Unit B and C, freshwater ostracods occur for the first time and freshwater taxa in general dominate. Fragments of shells, plants, insects and charcoal disappear in Unit B and C (Fig. 5).

Chronostratigraphy
Micropaleontological, organic and inorganic analyses show that three different depositional settings exist for the sediment sequence from Voëlvlei (Figs. 3, 5). In the lower part of the sequence (Unit A), the high concentrations of Na, TOC and N; high Br/Zr ratios; and a marine/brackish gastropod species assemblage indicate a marine/brackish depositional setting from 13 to 4.51 m sediment depth (Figs. 2, 3, 5). Corresponding reservoir-corrected bulk sediment 14 C ages as well as the terrestrial organic plant-macro particle and charcoal 14 C ages are all stratigraphically consistent and range from 8510 +280 / −200 to 6310 +450 / −210 cal BP suggesting that the sediments were rapidly deposited (Fig. 2). The overlap of the bulk and the plant-macro particle age and the bulk and charcoal ages confirm the suitability of the reservoir correction and application of the Marine20 calibration curve  for bulk samples in Unit A. The only exception are the bulk 14 C ages at 11.19 and 5.66 m sediment depth, which are too young for their stratigraphic position. This is probably due to an increased input of terrestrial organic carbon during this time, indicated by corresponding high inputs of sand, reduced Al contents and lower Al/Zr ratios (Figs. 2, 3). This likely added less 14 C-depleted material, affecting the marine organic carbon stock and diluting the reservoir effect. As we cannot calculate the precise contribution of terrestrial organic carbon, these 14 C ages were excluded from the age model in a second modelling iteration Hogg et al., 2020) Figure 5. Distribution of micropaleontological taxa and charcoal in core VOV16. Units are based on distribution of micropaleontological taxa and lithological characteristics. Additional macrofossils were picked from the lowest part of the core where they are relatively abundant. All abundances are given semi-quantitatively (0 -absent, 1 -rare, 2 -common, 3 -abundant, 4 -very abundant). Single occurrences of taxa are indicated by empty circles and abbreviated name (Nassarius kraussianus, Natica tecta, Bolivina sp., Spirillina sp., Cribroelphidium articulatum, trochamminid foraminifer, Aglaiella, myodocopid ostracod, Cyprilla humilis).
( Fig. 2). At Voëlvlei, n-alkanes show a clear dominance of the C 31 and C 33 homologues and thus are of terrestrial origin Chambers et al., 2014;. We therefore calibrated them with the terrestrial SHCal20 calibration curve. We note that a minor contribution of the shorter chain lengths (< C 25 ) to the dated compoundclass n-alkane samples, which are potentially synthesized by aquatic plants, for example, and may show a marine reservoir effect, leading to overly old 14 C ages when calibrated using a SHCal20 calibration curve. The aforementioned dominance of long-chain n-alkanes, however, makes a terrestrial origin likely and the distinct overlap with marine-calibrated and reservoir-corrected bulk 14 C ages support the application of the SHCal20 calibration curve to the compound-class n-alkane samples (Fig. 2).
At 4.51 m sediment depth, a distinct shift in the depositional setting occurred with decreasing Na concentrations, Br/Zr ratios, and TOC and N concentrations, pointing towards a change from a marine/brackish environment to a fresher lacustrine environmental setting (Figs. 3,5). At these depths bulk 14 C ages are much older than their stratigraphic position, possibly reflecting the deposition of degraded, pre aged sediments via erosion of old organic carbon from deeper soil horizons within the catchment Douglas et al., 2018;Haas et al., 2020). The OSL age (3.50 m sediment depth) is stratigraphically more consistent than the bulk 14 C ages. Since a MAM was applied it is likely that this age more closely reflects the timing of deposition, but the presence of partially or unbleached sand grains also potentially implies an input from older sediments at this time, supporting the interpretation of the bulk 14 C age overestimations. Compound-class n-alkane 14 C ages are stratigraphically consistent in Units B and C. We therefore infer that leaf waxes are derived from topsoils in the catchment and are rapidly deposited in the sediment archive Haas et al., 2017). Therefore, the incorporated climate signal should be close to that of the timing of deposition, and n-alkane-based proxies, i.e. δ 2 H n-alkane and δ 13 C n-alkane , yield paleoenvironmental information that can also be interpreted robustly.

Marine influence and lake development
South Africa, especially the southern Cape coast, is known to have experienced distinct environmental changes related to relative sea level fluctuations during the Holocene (Cooper et al., 2018;Kirsten et al., 2018b;Marker and Miller, 1993;Reddering, 1988;Wündsch et al., 2018Wündsch et al., , 2016a. Voëlvlei can contribute to our understanding of sea level changes during the Holocene due to its present location at an elevation of 5 m a.s.l. and the high temporal resolution of its Early and Middle Holocene depositional record. Inferred from lithological characteristics, paleoecological and elemental analyses, which indicate variable intrusion of marine water, the record suggests three eco-zones, of which one can be subdivided in two subzones (Units A1, A2, and B and C) (Fig. 6).
During the period between 8440 +200 / −250 and 7070 +160 / −200 cal BP (Unit A1), the dominant gastropod is Turritella capensis, which is a common species in sandy coastal lagoons in South Africa today (Branch et (Compton, 2006;Cooper et al., 2018) and (g) the lake floor of Voëlvlei relative to the present sea level.
al., 2010), preferring the mid-intertidal zone (Walters and Griffiths, 1987) (Figs. 5, 6). The rarer gastropod Assiminea globulus is known to be abundant on upper intertidal mudflats of South African estuaries (Barnes, 2018), whereas the small gastropod Hydrobia prefers the upper salt marshes of South Africa (Branch et al., 2010). The single occurrences of Natica tecta and Nassarius kraussianus point to estuarine mudflats (Branch et al., 2010) (Figs. 5, 6). Abundant brackish water ostracods and foraminifera indicate permanent water cover at the coring site during this period. In summary, these fauna reflect an estuarine and shallow subtidal environment close to intertidal mudflats. The abundance of plant remains and fragments of shells and fish bones point to shallow water as well. Therefore, we assume the respective core depths reflect an elevation slightly below the past sea level, which is in line with the paleo-surface of Lake Voëlvlei, i.e. 8.0 to 3.5 m below present sea level (b.s.l.) from 8440 +200 / −250 to 7070 +160 / −200 cal BP (Fig. 6). The low diversity of foraminifera is typical for lower salinity conditions; the dominant small and unornamented Ammonia species, Quinqueloculina and Haynesina, are consistent with this brackish water inference (Murray, 2006) (Figs. 5, 6). The same applies to the brackish water ostracod fauna, which are dominated by the estuarine species Sulcostocythere knysnaensis and Loxoconcha parameridionalis (Fürstenberg et al., 2017;Kirsten et al., 2018a) (Figs. 5, 6). The complete absence of open marine microfossils, especially planktonic foraminifera and echinoderm fragments, indicates an inner estuarine position in the core locality without direct marine inflow, even under high-energy conditions. The few bryozoan statoblasts and gyrogonites of charophytes all derive from freshwater environments (Frenzel, 2019;Kirsten et al., 2018a) and point to an unconfined exchange with river water. The highest salinity for the entire core is probably reached at 7090 +170 / −200 cal BP when the abundance and diversity of the foraminifera reach a maximum and a myodocopid ostracod was found (Figs. 5, 6). The coring position was ca. 3.65 m b.s.l. during this time (Fig. 6). Therefore, sea level was likely at the present height or slightly lower.
Between 7070 +170 / −200 and 6420 +130 / −140 cal BP (Unit A2), continued inner estuarine, brackish water conditions are indicated by a similar assemblage of taxa as observed in Unit A1. Intertidal gastropods typical for Unit A1, however, are now very rare, suggesting decreasing tidal influence and probably decreasing salinity (Figs. 5, 6). Salt marsh foraminifera occur only in Unit A2, albeit in low numbers (Fig. 6). This points to a very shallow water depth under marine/brackish conditions and a close shoreline (Strachan et al., 2017), which is in good agreement with the occurrence of fruits and seeds (only occurring in this unit), and the abundant plant remains. However, the high variability in the abundance of many taxa indicates unstable conditions compared to Unit A1. Freshwater inflow is implied by floating bryozoan statoblasts and charophyte gyrogonites (Frenzel, 2019) (Figs. 5, 6). The respective paleosurface of Lake Voëlvlei was ca. 3.5 m b.s.l. to 0.5 m a.s.l. from 7070 +170 / −200 to 6420 +130 / −140 cal BP (Fig. 6), indicating a sea level higher than the present during this time.
Elemental (Br/Zr ratios, Na contents) and mineralogical data support phases of marine water intrusion at Voëlvlei. Bromide salts are common in sea water but occur in very low concentrations in freshwater systems (Song and Müller, 1993). In aqueous environments the dominant species is Br − , which substitutes the salt constituent chloride (Cl − ) in the sea salt lattice during crystallization (Ullman, 1995). Halite (NaCl) also originates from marine waters, and thus Br, Na and halite can be used as indicators of marine water intrusion (Babel and Schreiber, 2014;Olsen et al., 2012;Wündsch et al., 2018). High values of the marine indicators from the elemental analysis (Na, Br/Zr ratios) and the presence of halite in the sediments are in good agreement with our micropaleontological data indicating marine water intrusions in Units A1 and A2 (8440 +200 / −250 to 6420 +130 / −140 cal BP; Fig. 6). Furthermore, the occurrence of calcite in both units as well as gypsum from 7260 +110 / −210 to 7180 +150 / −200 cal BP and 6590 +150 / −160 to 6550 +150 / −160 cal BP also points towards shallow water conditions as also indicated by the micropaleontology results (Fig. 6). Such shallow water conditions support bioproductivity (BiSi) increases in Unit A1 and A2 which reached a maximum at ca. 6740 +170 / −170 cal BP (Fig. 6). The elevation of the lake floor rapidly increases from ca. 8 m b.s.l. at 8440 +200 / −250 (base of the VOV16 record) to ca. 0.5 m a.s.l. at 6420 +130 / −140 cal BP (Unit A1 and A2; Fig. 6). Therefore, results imply a higher sea level compared to today until 6420 +130 / −140 cal BP and thus lend further support for a proposed local/regional Holocene relative sea level maximum of about +3.8 m a.s.l. (ca. 7600-5800 cal BP) (Cooper et al., 2018) on the southern Cape.
The micropaleontological associations of Units B and C (6420 +130 / −140 cal BP until today) are completely different from those of Unit A1 and A2 (Figs. 5, 6). Many samples are devoid of microfossils and freshwater taxa dominate other samples. One single test of the foraminifer Ammonia sp. found in the lowermost part of unit B is the only brackish water taxon (Murray, 2006). This was likely reworked from older sediments where this species is very abundant (Figs. 5,  6). Therefore, we assume freshwater or athalassic conditions for samples with aquatic taxa. New to Unit B and C are the freshwater ostracod species Sarscypridopsis aculeata and ephippia of cladocerans, which are typical of non-permanent water bodies (Frenzel, 2019;Meisch, 2000). They likely reflect the transformation of Voëlvlei to a non-permanent lake and terrestrial habitat. The lower part of Unit B contains variable amounts of plant remains and charcoal, indicating river transport. Their later disappearance points to an isolation of the basin while silting up.
Overall, the results only provide evidence of marine water intrusion in the Voëlvlei system rather than sea level index points and thus cannot be used to generate or corroborate an exact relative sea-level curve (e.g. Compton, 2006;Cooper et al., 2018). However, in comparison to local/regional records the results of this study are in line with the findings from Eilandvlei (Kirsten et al., 2018b;Wündsch et al., 2018) and Groenvlei (Wündsch et al., 2016a), which indicate a rising sea level during the Early Holocene and a high stand during the Middle Holocene. The proxies used in this study generally support the assumed sea-level evolution at the southern Cape coast of South Africa (Kirsten et al., 2018b;Marker and Miller, 1993;Reddering, 1988;Wündsch et al., 2018Wündsch et al., , 2016a. In (supra-)regional comparison, these results are also in line with studies from the west coast of South Africa (Baxter and Meadows, 1999;Carr et al., 2015;Kirsten et al., 2020) and Namibia (Compton, 2006), as reviewed in Cooper et al. (2018), showing a rapid sea level rise to a maximum of +3.8 m a.s.l. (ca. 7600-5800 cal BP), followed by a decrease to +1 m a.s.l. (5300-4200 cal BP) and a relatively constant sea level, comparable to the present, thereafter (Fig. 6).

Paleoenvironmental and paleoclimate evolution
The main driver of the δ 2 H n-alkane signal at Voëlvlei is the δ 2 H variability of the precipitation source  although it has to be noted that variations in the vegetation composition, e.g. varying biosynthetic fractionation, evapotranspirative enrichment and water use efficiency, may all potentially alter the δ 2 H n-alkane signal (Hou et al., 2007;Sachse et al., 2012). However, modern precipitation suggests that westerly-derived precipitation is 2 Hdepleted and easterly-derived precipitation is 2 H-enriched (Bowen, 2018;Bowen et al., 2005;Braun et al., 2017;Harris et al., 2010; Table 1). This is also imprinted in modern reference material from topsoils Strobel et al., 2020) (mean δ 2 H n-alkane values from topsoils from WRZ and SRZ are −145 ± 10 ‰ and −135 ± 10 ‰, respectively), pointing to the "source effect" as the domi-nant driver of the δ 2 H n-alkane signal at Voëlvlei. δ 2 H n-alkane values from Voëlvlei are moderate from 8440 +200 / −250 to 7070 +160 / −200 cal BP, which imply a year-round precipitation regime, i.e. contributions of both westerly-and easterlyderived precipitation (Fig. 7a).
Year-round precipitation is accompanied by less negative δ 13 C n-alkane values between 8440 +200 / −250 and 7070 +160 / −200 cal BP, likely indicating overall dry conditions for this period at Voëlvlei (Fig. 7e). These dry conditions potentially led to a sparse vegetation cover and runoff (induced by occasional events), which likely carried extremely variable grain sizes and amounts of allochthonous input (Al concentration, Al/Zr ratios) (Fig. 7g). Minimum values in Al concentrations and Al/Zr ratios are assumed to be caused by layers consisting almost exclusively of sand (SiO 2 ) and thus imply high allochthonous inputs (Figs. 2,7h) at this time. Similarly at Vankervelsvlei, input of reworked soil material during this time is interpreted to be the result of low vegetation cover in the catchment (Strobel et al., 2019). This interpretation of dry conditions between 8440 +200 / −250 and 7070 +160 / −200 cal BP is also in line with the findings from Eilandvlei, where low Afrotemperate forest (AFT) pollen percentages are interpreted as overall dry conditions between ca. 8500 and 7000 cal BP   (Fig. 7d), and a study at Seweweekspoort, where less positive δ 15 N values are also interpreted to indicate moderately dry conditions   (Fig. 7c) resulting in a consistent moisture signal in the YRZ during this time.
At 7020 +170 / −200 cal BP, moister conditions are indicated by a distinct shift to more negative δ 13 C n-alkane values and evidence for increased allochthonous inputs (Al content, Al/Zr ratios), which lasts until the present. However, there are decreases in allochthonous inputs from 6420 +130 / −140 to 6080 +320 / −240 cal BP, and at 4480 +580 / −540 and 3530 +500 / −470 cal BP (Fig. 7g). The latter event is accompanied by distinctly negative δ 13 C values (Fig. 7e) whereas the earlier one only shows slightly more negative values. This potentially indicates a denser vegetation cover and reduced runoff-induced sediment input, which are supported by the presence of freshwater taxa during this time (Fig. 5). However, drier conditions inferred from reduced allochthonous input could also be possible due to drought-induced changes of the vegetation composition indicated by a higher proportion of plants using CAM metabolism. These could possibly reduce drought stress and therefore the influence of water use efficiency to 13 C-enrichment in the n-alkanes (Diefendorf and Freimuth, 2017). Overall it is difficult to interpret the variations in the core between 4480 +580 / −540 and 3530 +500 / −470 cal BP. They interrupt a trend towards moister conditions but cannot be ultimately attributed to moister or drier conditions. Noticeable in the record are three data points at 6670 +170 / −170 , 6560 +150 / −160 and 6440 +130 / −140 cal BP showing the least negative δ 13 C n-alkane values, which exceed all other measurements and are accompanied by dis-tinctly negative δ 2 H n-alkane values (grey dots in Fig. 7a, d). For those three samples, the climate signal is possibly overprinted by local effects of plants that grow at shallower water depths and/or at the shoreline of the lake. While less negative δ 13 C n-alkane values point to the presence of plants using C 4 and/or CAM metabolism Carr et al., 2015;Diefendorf and Freimuth, 2017), the δ 2 H n-alkane signal with more negative values can be biased by the photosynthetic mode and also salinity (Aichner et al., 2017;Feakins and Sessions, 2010;Ladd and Sachs, 2012;Sachse et al., 2012). Therefore, disentangling the regional climate and local overprinting is challenging for the three aforementioned data points. However, since allochthonous input is high (Al concentration, Al/Zr ratios; Fig. 7g), we hypothesize overall moist conditions since 7020 +170 / −200 cal BP.
Moister conditions are accompanied by more negative δ 2 H n-alkane values from 7020 +170 / −200 to 6420 +130 / −140 cal BP, which are likely indicative of an increased contribution of westerly-derived precipitation (Fig. 7a). At the same time as the aforementioned reduction of the allochthonous inputs, there is also a marked shift to less negative δ 2 H n-alkane values at 6420 +130 / −140 implying an increased contribution of easterly-derived precipitation (Fig. 7a). Afterwards, δ 2 H n-alkane values are variable but have an overall trend from less negative to more negative δ 2 H n-alkane values lasting until 2060 +140 / −200 cal BP, which is likely suggest a return from an easterly-to a more dominantly westerly-derived precipitation regime (Fig. 7a). However, already at the beginning, this period is interrupted by a short phase of somewhat more negative values from 5560 +430 / −470 to 5150 +340 / −330 , indicating more westerly-derived precipitation (Fig. 7a). From 2060 +140 / −200 cal BP until the present, less negative δ 2 H n-alkane values imply an increased contribution of easterly-derived precipitation (Fig. 7a), which occurs within a trend towards moister conditions (δ 13 C n-alkane ) (Fig. 7e).
Moist conditions between ca. 7000 and ca. 4700 cal BP are in line with studies from the central southern Cape coast Strobel et al., 2019;. At ca. 4700 cal BP, a rapid shift in climate was also detected in numerous studies located in the YRZ. In the Wilderness area, dry conditions and low wind-driven evapotranspiration were reconstructed during this time Strobel et al., 2019;Wündsch et al., 2018) (Fig. 7d). After 4700 cal BP, the overall trend in increasing moisture availability at the central southern Cape coast is again in line with our findings at Voëlvlei (du Plessis et al., 2020;Quick et al., 2018;Strobel et al., 2019;Wündsch et al., 2018). δ 2 H n-alkane and δ 13 C n-alkane from marine sediments recovered off the mouth of the Gouritz River (GeoB18308-2; Fig. 1) show a very similar pattern to Voëlvlei over the past 4000 cal BP (Hahn et al., 2017) (Fig. 7). Those findings thus support our findings from Voëlvlei and suggest a consistent regional picture for the southern Cape coast. Figure 7. (a) δ 2 H n-alkane from Voëlvlei compared to (b) δ 2 H n-alkane from GeoB18308-1 (Hahn et al., 2017), (c) δ 15 N from Seweweekspoort , (d) Afrotemperate forest pollen (AFT) from Eilandvlei , (e) δ 13 C n-alkane from Voëlvlei, (f) δ 13 C n-alkane from GeoB18308-1 (Hahn et al., 2017), and (g) aluminium (Al) content and Log Al/Zr ratio from Voëlvlei. The locations of the studies are depicted in Fig. 1. At Seweweekspoort, ca. 100 km inland, δ 15 N values follow the pattern of δ 2 H n-alkane from Voëlvlei, and climate is thought to have been dry from ca. 7000 to ca. 4700 cal BP and ca. 3000 to ca. 1000 cal BP due to reduced easterly-related precipitation   (Fig. 7c). Intense short-term variations to wetter conditions at Seweweekspoort contradict the aforementioned results from the Wilderness area at ca. 5700 and ca. 4900 cal BP but are well reflected by spikes in the δ 2 H n-alkane record from Voëlvlei which shows shifts towards less negative δ 2 H n-alkane values and increased easterlyrelated precipitation contributions (Fig. 7a, c).
In this context, previous studies hypothesize that the nondominant component of the temperate (westerly)/tropical (easterly) dynamic determines hydro-climatic variability in southern Africa (e.g. Chase et al., 2017Chase et al., , 2015. More precisely, this means that although easterly-derived (westerlyderived) precipitation is relatively low in the western (eastern) parts of the YRZ, variability in the overall moisture signal is thought to strongly depend on moisture contribution from this minor component (e.g. Chase et al., 2017Chase et al., , 2015. The δ 2 H n-alkane values of this study imply increased westerly-derived precipitation during drier phases at Seweweekspoot and increased easterly-derived precipitation during moister phases at Seweweekspoort (Fig. 7c), at a time generally characterized by a trend towards moister conditions at Voëlvlei (Fig. 7e). Moist conditions between 7000 and 5500 cal BP were also found in the north-southtrending Cape Fold Mountains in the Western Cape (Chase et al., 2019, and references therein), indicating an overall moist period for south-western South Africa during this time which is perfectly in line with the increased westerly precipitation contribution at Voëlvlei.
However, an opposite pattern is proposed to exist between sites located inland of the Cape Fold mountains (e.g. Seweweekspoort, Katbakkies Pass) and coastal sites from the southern Cape coast (e.g. Eilandvlei) Quick et al., 2018). There are still discrepancies due to the comparison of different proxies in these studies compared to those from the southern Cape coast. Comparing our δ 2 H n-alkane record from Voëlvlei, which is representative of fluctuations in the major moisture sources, with the δ 15 N from Seweweekspoort , in detail many similarities are apparent, confirming the assumption that hydrological variations at Seweweekspoort are distinctly driven by variations in easterly-derived precipitation. However, on the South Coast of South Africa including Voëlvlei the overall moisture evolution is not restricted to one moisture source and thus driven by a combination of both westerly-and easterly-derived precipitation.
The topmost sediments of Voëlvlei mainly consist of two fluvial flooding events (at 830 +110 / −110 and 10 +10 / −10 cal BP) which deposited large amounts of sediment (ca. 15 and ca. 30 cm, respectively; Figs. 2, 7). Considering the chronological uncertainties (dating and modelling errors), the most recent flooding event can likely be associated with the so called "Laingsburg flood" from January 1981 (Damm and Hagedorn, 2010 and references therein). In marine sediments recovered off the mouth of the Gouritz River (GeoB18308-2; Fig. 1) major flooding events were detected between 650 and 300 cal BP and thus support our findings (Hahn et al., 2017).

Conclusions
Our multi-proxy record from Voëlvlei provides new insights into the sea level and paleoclimate history of the past 8.5 ka at the southern Cape Coast. Various dating approaches including OSL and 14 C of bulk TOC, organic macrofossils, charcoal and n-alkanes were applied to establish a robust chronology. Three temporal eco-zones can be characterized for Voëlvlei. These are related to sea level variations and are in line with results of regional investigations in implying a subtidal marine-brackish followed by an intratidal brackish system from 8440 +200 / −250 to 6420 +130 / −140 cal BP, and a freshwater to terrestrial system due to silting up of Voëlvlei thereafter. In terms of climate, the combined biomarkersedimentological approach allows the additional identification of changing precipitation sources, in combination with an estimation of moisture availability. Dry conditions accompanied by a year-round precipitation contribution prevail from 8440 +200 / −250 to 7070 +160 / −200 cal BP followed by moister conditions and more westerly-derived precipitation contributions from 7070 +160 / −200 to 6420 +130 / −140 cal BP, This was followed by a distinct shift to an easterly dominance at 6420 cal BP, before an overall trend to a westerly dominance, which lasted until 2060 cal BP. This is followed by a trend towards an easterly dominance, but both phases show several intense short-term variations.
These results are in good agreement with previous studies with regard to sea-level change on the western and southern coast of South Africa. The paleohydrological evolution of Voëlvlei is comparable to previous investigations from the Southern Cape coast during the past 8.5 ka, demonstrating an interplay of both westerly-and easterly-derived precipitation contribution. Hydrological variations at Seweweekspoort, located within the interior of central southern South Africa, a region hypothesized to show distinct hydroclimatic trends (relative to the coastal zones) through the Holocene, show many similarities to the Voëlvlei δ 2 H n-alkane record, suggesting a need to further consider in more detail regional-scale complexities and coastal-inland contrasts in paleoclimatic trajectories. Thus, the Voëlvlei δ 2 H n-alkane record provides valuable insights into the source of precipitation at the Southern Cape coast during the past 8.5 ka.
Data availability. All data within the paper will be available at PANGAEA.
Author contributions. PS, RZ and TH designed the study. PS, MB, ASC, PF, BK, GS, JS and SS carried out the laboratory and analytical analyses. All authors contributed to the writing of the paper and data discussion.
Competing interests. The authors declare that they have no conflict of interest.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements. Sediment cores were recovered within the RAiN project funded by the German Federal Ministry of Education and Research (BMBF). Paul Strobel gratefully acknowledges the support by a fellowship from the state of Thuringia (Landesgraduiertenstipendium). Particularly acknowledged are Cristopher Berndt, Nico Blaubach, Tobias Eggert, Christian Gregori, Patricia Rauh, Mike Steinich and Magdalena Wagner for assistance in the lab. Gerhard Daut is thanked for supporting the XRD measurements. Nadia du Plessis, Jussi Baade, Sayed Hess, Christian Gregori, Roland Mäusbacher and Michael E. Meadows are thanked for their fieldwork contributions and constructive discussions. The participants of the 2019 course on Paleoenvironmental Analysis at Friedrich Schiller University Jena supported us in micropaleontological sample preparation and picking of microfossils for an initial batch of subsamples.
Review statement. This paper was edited by Keely Mills and reviewed by Hayley Cawthra and three anonymous referees.