Late Quaternary palaeoenvironmental reconstruction from Lake Ohrid

Lake Ohrid is a large, deep, ancient lake located on the Balkan Peninsula and has an extraordinary degree of endemic biodiversity. A deep drilling campaign was carried out in 2013 as part of the Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) project to better understand the influence of climate and environmental change on evolutionary patterns and endemism. Initial results from a 569 m sediment succession drilled in the centre of the basin indicate continuous lacustrine conditions over the past ca. 1.2 Ma. This thesis presents new stable isotope data from carbonate within the upper 248 m of the composite profile, covering the last ca. 640 ka at a millennial-scale resolution (≈ 0.5 ka). Isotope data are used to provide a long-term palaeoenvironmental reconstruction, which is supported by a quantitative model of modern lake water isotope composition and a detailed multi-proxy investigation of climate evolution since the last glacial period. Water balance calculations confirm Ohrid to be an evaporative system with a complex hydrology. Variations in the isotope composition of modern lake water are suggested to represent long-term fluctuations in water balance. A trend from wetter to drier conditions through the Holocene is consistent with regional and hemispheric processes related to changes in insolation. Over the last ca. 640 ka, endogenic calcite is precipitated in abundance during warm stages, however carbonate is negligible during glacial periods with the exception of discrete horizons comprising early diagenetic authigenic siderite. The oxygen isotope composition of lake water, calculated using calcite and siderite, indicates more evaporated conditions during warm stages and fresher lake water in glacial phases. Low frequency variability shown by calcite isotope data suggests relatively stable conditions before ca. 450 ka, a transition to a wetter climate between ca. 400-250 ka, and a trend to drier climate conditions after ca. 250 ka. Higher frequency millennial-scale oscillations and abrupt climate events observed during warm stages are likely associated with regional climate change as a function of orbital forcing. This study emphasises the potential of Lake Ohrid as a valuable archive of climate change in the central Mediterranean region and demonstrates the efficacy of isotope data for palaeoenvironmental reconstruction at Lake Ohrid

Complexity of diatom response to Lateglacial and Holocene climate and environmental change in ancient, deep and oligotrophic Lake Ohrid (Macedonia and Albania)

© Author(s) 2016. Lake Ohrid (Macedonia and Albania) is a rare example of a deep, ancient Mediterranean lake and is a key site for palaeoclimate research in the northeastern Mediterranean region. This study conducts the analysis of diatoms as a proxy for Lateglacial and Holocene climate and environmental change in Lake Ohrid at a higher resolution than in previous studies. While Lake Ohrid has the potential to be sensitive to water temperature change, the data demonstrate a highly complex diatom response, probably comprising a direct response to temperature-induced lake productivity in some phases and an indirect response to temperaturerelated lake stratification or mixing and epilimnetic nutrient availability in others. The data also demonstrate the possible influence of physical limnological (e.g. the influence of wind stress on stratification or mixing) and chemical processes (e.g. the influence of catchment dynamics on nutrient input) in mediating the complex response of diatoms. During the Lateglacial (ca. 12 300-11 800 cal yr BP), the low-diversity dominance of hypolimnetic Cyclotella fottii indicates low lake productivity, linked to low water temperature. Although the subsequent slight increase in small, epilimnetic C. minuscula during the earliest Holocene (ca. 11 800-10 600 cal yr BP) suggests climate warming and enhanced stratification, diatom concentration remains as low as during the Lateglacial, suggesting that water temperature increase was muted across this major transition. The early Holocene (ca. 10 600-8200 cal yr BP) is characterised by a sustained increase in epilimnetic taxa, with mesotrophic C. ocellata indicating high water-temperature-induced productivity between ca. 10 600-10 200 cal yr BP and between ca. 9500-8200 cal yr BP and with C. minuscula in response to low nutrient availability in the epilimnion between ca. 10 200-9500 cal yr BP. During the middle Holocene (ca. 8200-2600 cal yr BP), when sedimentological and geochemical proxies provide evidence for maximum Holocene water temperature, anomalously low C. ocellata abundance is probably a response to epilimnetic nutrient limitation, almost mimicking the Lateglacial flora apart from the occurrence of mesotrophic Stephanodiscus transylvanicus in the hypolimnion. During the late Holocene (ca. 2600 cal yr BP-present), high abundance and fluctuating composition of epilimnetic taxa are probably a response more to enhanced anthropogenic nutrient input, particularly nitrogen enrichment, than to climate. Overall, the data indicate that previous assumptions concerning the linearity of diatom response in this deep, ancient lake are invalid, and multi-proxy analysis is essential to improve understanding of palaeolimnological dynamics in future research on the long, Quaternary sequence

Quantitative reconstruction of early Holocene and last glacial climate on the Balkan Peninsula using coupled hydrological and isotope mass balance modelling

We investigate the modern hydrology of Lake Ohrid (Macedonia/Albania) using a combined hydrological and isotope-based modelling approach and present a new evaluation of contemporary water balance and palaeoclimate estimates. The combined model is able to estimate hydrological components that cannot be directly measured, and indicates that sublacustrine spring inflow is in the order of 50% higher than previous estimates and groundwater outflow comprises approximately a third of overall water outflow. In combination with sediment core oxygen isotope data, we used the combined model to quantitatively reconstruct past climate, in particular precipitation, during the early Holocene and last glacial period. Calculated precipitation in the early Holocene was higher than the value for present day and was approximately 44% lower than present during the last glacial, assuming the majority of precipitation fell as snow. The estimated amount of precipitation in the last glacial would have been high enough to provide refugial conditions at Lake Ohrid and to support the continuous existence of arboreal vegetation in the catchment. The improved understanding of the modern isotope hydrology of Lake Ohrid is fundamental for explaining the systematics of past isotope variation and providing context for extended sediment records from the lake, which will provide longer-term palaeoclimate reconstructions covering multiple glacial-interglacial cycles

Western Mediterranean climate and environment since Marine Isotope Stage 3: a 50,000-year record from Lake Banyoles, Spain

We present new stable isotope (δ¹⁸Ocalcite and δ¹³ Ccalcite) and diatom data from a 67-m sediment core (BAN II) from Lake Banyoles, northeastern Spain. We reassessed the chronology of the sequence by correlating stable isotope data with a shorter U-series-dated record from the lake, confirming a sedimentological offset between the two cores and demonstrating that BAN II spans Marine Isotope Stages (MIS) 3–1. Through comparison with previous records, the multi-proxy data are used to improve understanding of palaeolimnological dynamics and, by inference, western Mediterranean climate and environmental change during the past ca. 50,000 years. Three main zones, defined by isotope and diatom data, correspond to the MIS. The basal zone (MIS 3) is characterised by fluctuating δ¹⁸Ocalcite and benthic diatom abundance, indicating a high degree of environmental and climate variability, concomitant with large lake-level changes. During the full glacial (MIS 2), relatively constant δ¹⁸Ocalcite and a poorly preserved planktonic-dominated diatom assemblage suggest stability, and intermittently, unusually high lake level. In MIS 1, δ¹⁸Ocalcite and δ¹³Ccalcite initially transition to lower values, recording a pattern of Late Glacial to Holocene change that is similar to other Mediterranean records. This study suggests that Lake Banyoles responds limnologically to changes in the North Atlantic ocean–atmosphere system and provides an important dataset from the Iberian Peninsula, a region in need of longer-term records that can be used to correlate between marine and terrestrial archives, and between the western and eastern Mediterranean

Environmental control on the occurrence of high-coercivity magnetic minerals and formation of iron sulfides in a 640 ka sediment sequence from Lake Ohrid (Balkans)

The bulk magnetic mineral record from Lake Ohrid, spanning the past 637 kyr, reflects large-scale shifts in hydrological conditions, and, superimposed, a strong signal of environmental conditions on glacial–interglacial and millennial timescales. A shift in the formation of early diagenetic ferrimagnetic iron sulfides to siderites is observed around 320 ka. This change is probably associated with variable availability of sulfide in the pore water. We propose that sulfate concentrations were significantly higher before  ∼  320 ka, due to either a higher sulfate flux or lower dilution of lake sulfate due to a smaller water volume. Diagenetic iron minerals appear more abundant during glacials, which are generally characterized by higher Fe / Ca ratios in the sediments. While in the lower part of the core the ferrimagnetic sulfide signal overprints the primary detrital magnetic signal, the upper part of the core is dominated by variable proportions of high- to low-coercivity iron oxides. Glacial sediments are characterized by high concentration of high-coercivity magnetic minerals (hematite, goethite), which relate to enhanced erosion of soils that had formed during preceding interglacials. Superimposed on the glacial–interglacial behavior are millennial-scale oscillations in the magnetic mineral composition that parallel variations in summer insolation. Like the processes on glacial–interglacial timescales, low summer insolation and a retreat in vegetation resulted in enhanced erosion of soil material. Our study highlights that rock-magnetic studies, in concert with geochemical and sedimentological investigations, provide a multi-level contribution to environmental reconstructions, since the magnetic properties can mirror both environmental conditions on land and intra-lake processes

Western Mediterranean climate and environment since Marine Isotope Stage 3: a 50,000-year record from Lake Banyoles, Spain

We present new stable isotope (δ18Ocalcite and δ13Ccalcite) and diatom data from a 67-m sediment core (BAN II) from Lake Banyoles, northeastern Spain. We reassessed the chronology of the sequence by correlating stable isotope data with a shorter U-series-dated record from the lake, confirming a sedimentological offset between the two cores and demonstrating that BAN II spans Marine Isotope Stages (MIS) 3–1. Through comparison with previous records, the multi-proxy data are used to improve understanding of palaeolimnological dynamics and, by inference, western Mediterranean climate and environmental change during the past ca. 50,000 years. Three main zones, defined by isotope and diatom data, correspond to the MIS. The basal zone (MIS 3) is characterised by fluctuating δ18Ocalcite and benthic diatom abundance, indicating a high degree of environmental and climate variability, concomitant with large lake-level changes. During the full glacial (MIS 2), relatively constant δ18Ocalcite and a poorly preserved planktonic-dominated diatom assemblage suggest stability, and intermittently, unusually high lake level. In MIS 1, δ18Ocalcite and δ13Ccalcite initially transition to lower values, recording a pattern of Late Glacial to Holocene change that is similar to other Mediterranean records. This study suggests that Lake Banyoles responds limnologically to changes in the North Atlantic ocean–atmosphere system and provides an important dataset from the Iberian Peninsula, a region in need of longer-term records that can be used to correlate between marine and terrestrial archives, and between the western and eastern Mediterranean

Complexity of diatom response to Lateglacial and Holocene climate and environmental change in ancient, deep and oligotrophic Lake Ohrid (Macedonia and Albania)

Lake Ohrid (Macedonia and Albania) is a rare example of a deep, ancient Mediterranean lake and is a key site for palaeoclimate research in the northeastern Mediterranean region. This study conducts the analysis of diatoms as a proxy for Lateglacial and Holocene climate and environmental change in Lake Ohrid at a higher resolution than in previous studies. While Lake Ohrid has the potential to be sensitive to water temperature change, the data demonstrate a highly complex diatom response, probably comprising a direct response to temperature-induced lake productivity in some phases and an indirect response to temperature-related lake stratification or mixing and epilimnetic nutrient availability in others. The data also demonstrate the possible influence of physical limnological (e.g. the influence of wind stress on stratification or mixing) and chemical processes (e.g. the influence of catchment dynamics on nutrient input) in mediating the complex response of diatoms. During the Lateglacial (ca. 12 300-11 800 cal yr BP), the low-diversity dominance of hypolimnetic Cyclotella fottii indicates low lake productivity, linked to low water temperature. Although the subsequent slight increase in small, epilimnetic C. minuscula during the earliest Holocene (ca. 11 800-10 600 cal yr BP) suggests climate warming and enhanced stratification, diatom concentration remains as low as during the Lateglacial, suggesting that water temperature increase was muted across this major transition. The early Holocene (ca. 10 600-8200 cal yr BP) is characterised by a sustained increase in epilimnetic taxa, with mesotrophic C. ocellata indicating high water-temperature-induced productivity between ca. 10 600-10 200 cal yr BP and between ca. 9500-8200 cal yr BP and with C. minuscula in response to low nutrient availability in the epilimnion between ca. 10 200-9500 cal yr BP. During the middle Holocene (ca. 8200-2600 cal yr BP), when sedimentological and geochemical proxies provide evidence for maximum Holocene water temperature, anomalously low C. ocellata abundance is probably a response to epilimnetic nutrient limitation, almost mimicking the Lateglacial flora apart from the occurrence of mesotrophic Stephanodiscus transylvanicus in the hypolimnion. During the late Holocene (ca. 2600 cal yr BP-present), high abundance and fluctuating composition of epilimnetic taxa are probably a response more to enhanced anthropogenic nutrient input, particularly nitrogen enrichment, than to climate. Overall, the data indicate that previous assumptions concerning the linearity of diatom response in this deep, ancient lake are invalid, and multi-proxy analysis is essential to improve understanding of palaeolimnological dynamics in future research on the long, Quaternary sequence

Northern Mediterranean climate since the Middle Pleistocene: a 637 ka stable isotope record from Lake Ohrid (Albania/Macedonia)

Lake Ohrid (Macedonia/Albania) is an ancient lake with unique biodiversity and a site of global significance for investigating the influence of climate, geological, and tectonic events on the generation of endemic populations. Here, we present oxygen (δ18O) and carbon (δ13C) isotope data from carbonate over the upper 243 m of a composite core profile recovered as part of the Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) project. The investigated sediment succession covers the past ca. 637 ka. Previous studies on short cores from the lake (up to 15 m, < 140 ka) have indicated the total inorganic carbon (TIC) content of sediments to be highly sensitive to climate change over the last glacial–interglacial cycle. Sediments corresponding to warmer periods contain abundant endogenic calcite; however, an overall low TIC content in glacial sediments is punctuated by discrete bands of early diagenetic authigenic siderite. Isotope measurements on endogenic calcite (δ18Oc and δ13Cc) reveal variations both between and within interglacials that suggest the lake has been subject to palaeoenvironmental change on orbital and millennial timescales. We also measured isotope ratios from authigenic siderite (δ18Os and δ13Cs) and, with the oxygen isotope composition of calcite and siderite, reconstruct δ18O of lake water (δ18Olw) over the last 637 ka. Interglacials have higher δ18Olw values when compared to glacial periods most likely due to changes in evaporation, summer temperature, the proportion of winter precipitation (snowfall), and inflow from adjacent Lake Prespa. The isotope stratigraphy suggests Lake Ohrid experienced a period of general stability from marine isotope stage (MIS) 15 to MIS 13, highlighting MIS 14 as a particularly warm glacial. Climate conditions became progressively wetter during MIS 11 and MIS 9. Interglacial periods after MIS 9 are characterised by increasingly evaporated and drier conditions through MIS 7, MIS 5, and the Holocene. Our results provide new evidence for long-term climate change in the northern Mediterranean region, which will form the basis to better understand the influence of major environmental events on biological evolution within Lake Ohrid

Sediment residence time reveals Holocene shift from climatic to vegetation control on catchment erosion in the Balkans

Understanding the evolution of soil systems on geological time scales has become fundamentally important to predict future landscape development in light of rapid global warming and intensifying anthropogenic impact. Here, we use an innovative uranium isotope-based technique combined with organic carbon isotopes and elemental ratios of sediments from Lake Ohrid (North Macedonia/Albania) to reconstruct soil system evolution in the lake's catchment during the last ~16,000 cal yr BP. Uranium isotopes are used to estimated the paleo-sediment residence time, defined as the time elapsed between formation of silt and clay sized detrital matter and final deposition. The chronology is based on new cryptotephra layers identified in the sediment sequence. The isotope and elemental data are compared to sedimentary properties and pollen from the same sample material to provide a better understanding of past catchment erosion and landscape evolution in the light of climate forcing, vegetation development, and anthropogenic land use. During the Late Glacial and the Early Holocene, when wide parts of the catchment were covered by open vegetation, wetter climates promoted the mobilisation of detrital matter with a short paleo-sediment residence time. This is explained by erosion of deeper parts of the weathering horizon from thin soils. Detrital matter with a longer paleo-sediment residence time, illustrating shallow erosion of thicker soils is deposited in drier climates. The coupling between climatic variations and soil erosion terminates at the Early to Mid-Holocene transition as evidenced by a pronounced shift in uranium isotope ratios indicating that catchment erosion is dominated by shallow erosion of thick soils only. This shift suggests a threshold is crossed in hillslope erosion, possibly as a result of a major change in vegetation cover preventing deep erosion of thin soils at higher elevation. The threshold in catchment erosion is not mirrored by soil development over time, which gradually increases in response to Late Glacial to Holocene warming until human land use during the Late Holocene promotes reduced soil development and soil degradation. Overall, we observe that soil system evolution is progressively controlled by climatic, vegetation, and eventually by human land use over the last ~16,000 years

Silicic Acid Cycling in the Bering Sea During the Mid‐Pleistocene Transition

The rate of deep-ocean carbon burial is considered important for modulating glacial-interglacial atmospheric CO2 concentrations and global climate during the Quaternary. It has been suggested that glacial iron fertilization and increased efficiency of the biological pump in the Southern Ocean since the Mid-Pleistocene Transition (MPT) was key in lowering atmospheric pCO2 and facilitating rapid land ice accumulation. There is growing evidence that a similar mechanism may have existed in the subarctic Pacific Ocean, although this has not yet been assessed. Here, the silicon isotope composition of diatoms (δ30Sidiatom) from the Bering Sea upwelling region is used to assess the role of nutrient cycling on the subarctic Pacific biological pump during the MPT. Results show that during and after the “900 kyr event,” the high productivity green belt zone was characterized by low silicic acid utilization but high supply, coincident with the dominance of diatom resting spores. We posit that as nutrient upwelling was suppressed following pack ice growth and expansion of glacial North Pacific Intermediate Water (GNPIW), primary productivity became nitrate-limited and enhanced opal remineralization caused a relative increase in silicic acid supply. However, preferential preservation and higher cellular carbon content of diatom resting spores, as well as increased supply of iron from expanded sea ice, likely sustained the net efficiency of the Bering Sea biological pump through the MPT. Remnant iron and silicic acid may also have propagated into the lower subarctic Pacific Ocean through GNPIW, aiding a regionally efficient biological pump at 900 kyr and during post-MPT glacials