Search

This study reviews and synthesises existing information generated within the SCOPSCO (Scientific Collaboration on Past Speciation Conditions in Lake Ohrid) deep drilling project. The four main aims of the project are to infer (i) the age and origin of Lake Ohrid (Former Yugoslav Republic of Macedonia/Republic of Albania), (ii) its regional seismotectonic history, (iii) volcanic activity and climate change in the central northern Mediterranean region, and (iv) the influence of major geological events on the evolution of its endemic species. The Ohrid basin formed by transtension during the Miocene, opened during the Pliocene and Pleistocene, and the lake established de novo in the still relatively narrow valley between 1.9 and 1.3 Ma. The lake history is recorded in a 584 m long sediment sequence, which was recovered within the framework of the International Continental Scientific Drilling Program (ICDP) from the central part (DEEP site) of the lake in spring 2013. To date, 54 tephra and cryptotephra horizons have been found in the upper 460 m of this sequence. Tephrochronology and tuning biogeochemical proxy data to orbital parameters revealed that the upper 247.8 m represent the last 637 kyr. The multi-proxy data set covering these 637 kyr indicates long-term variability. Some proxies show a change from generally cooler and wetter to drier and warmer glacial and interglacial periods around 300 ka. Short-term environmental change caused, for example, by tephra deposition or the climatic impact of millennial-scale Dansgaard–Oeschger and Heinrich events are superimposed on the long-term trends. Evolutionary studies on the extant fauna indicate that Lake Ohrid was not a refugial area for regional freshwater animals. This differs from the surrounding catchment, where the mountainous setting with relatively high water availability provided a refuge for temperate and montane trees during the relatively cold and dry glacial periods. Although Lake Ohrid experienced significant environmental change over the last 637 kyr, preliminary molecular data from extant microgastropod species do not indicate significant changes in diversification rate during this period. The reasons for this constant rate remain largely unknown, but a possible lack of environmentally induced extinction events in Lake Ohrid and/or the high resilience of the ecosystems may have played a role.

The DEEP site sediment sequence obtained during the ICDP SCOPSCO project at Lake Ohrid was dated using tephrostratigraphic information, cyclostratigraphy, and orbital tuning through the marine isotope stages (MIS) 15-1. Although this approach is suitable for the generation of a general chronological framework of the long succession, it is insufficient to resolve more detailed palaeoclimatological questions, such as leads and lags of climate events between marine and terrestrial records or between different regions. Here, we demonstrate how the use of different tie points can affect cyclostratigraphy and orbital tuning for the period between ca. 140 and 70 ka and how the results can be correlated with directly/indirectly radiometrically dated Mediterranean marine and continental proxy records. The alternative age model presented here shows consistent differences with that initially proposed by Francke et al. (2015) for the same interval, in particular at the level of the MIS6-5e transition. According to this new age model, different proxies from the DEEP site sediment record support an increase of temperatures between glacial to interglacial conditions, which is almost synchronous with a rapid increase in sea surface temperature observed in the western Mediterranean. The results show how a detailed study of independent chronological tie points is important to align different records and to highlight asynchronisms of climate events. Moreover, Francke et al. (2016) have incorporated the new chronology proposed for tephra OH-DP-0499 in the final DEEP age model. This has reduced substantially the chronological discrepancies between the DEEP site age model and the model proposed here for the last glacial-interglacial transition.

A tephrostratigraphic record covering the Marine Isotope Stages (MIS) 1–15 was established for the DEEP site record of Lake Ohrid (Macedonia and Albania). Major element analyses (energy dispersive spectroscopy (EDS) and wavelength-dispersive spectroscopy (WDS)) were carried out on juvenile fragments extracted from 12 tephra layers (OH-DP-0115 to OH-DP-2060). The geochemical analyses of the glass shards of all of these layers suggest an origin in the Italian volcanic provinces. They include the Y-3 (OH-DP-0115, 26.68–29.42 ka cal BP), the Campanian Ignimbrite–Y-5 (OH-DP-0169, 39.6 ± 0.1 ka), and the X-6 (OH-DP-0404, 109 ± 2 ka) from the Campanian volcanoes, the P-11 of Pantelleria (OH-DP-0499, 133.5 ± 2 ka), the Vico B (OH-DP-0617, 162 ± 6 ka) from the Vico volcano, the Pozzolane Rosse (OH-DP-1817, 457 ± 2 ka) and the Tufo di Bagni Albule (OH-DP-2060, 527 ± 2 ka) from the Colli Albani volcanic district, and the Fall A (OH-DP-2010, 496 ± 3 ka) from the Sabatini volcanic field. Furthermore, a comparison of the Ohrid record with tephrostratigraphic records of mid-distal archives related to the Mediterranean area allowed the recognition of the equivalents of other less known tephra layers, such as the TM24a–POP2 (OH-DP-0404, 102 ± 2 ka) recognized in the Lago Grande di Monticchio and the Sulmona Basin, the CF-V5–PRAD3225 (OH-DP-0624, ca. 163 ± 22 ka) identified in the Campo Felice Basin and the Adriatic Sea, the SC5 (OH-DP-1955, 493.1 ± 10.9 ka) recognized in the Mercure Basin, and the A11/12 (OH-DP-2017, 511 ± 6 ka) sampled at the Acerno Basin, whose specific volcanic sources are still poorly constrained. Additionally, one cryptotephra (OH-DP-0027) was identified by correlation of the potassium X-ray flourescence (XRF) intensities from the DEEP site with those from a short core of a previous study from Lake Ohrid. In these cores, a maximum in potassium is caused by glass shards, which were correlated with the Mercato tephra (8.43–8.63 ka cal BP) from Somma–Vesuvius. The tephrostratigraphic work presented here allows, for the first time, the extension of a consistent part of the Middle Pleistocene tephrostratigraphy of Italian volcanoes as far as the Balkans. The establishment of the tephrostratigraphic framework for the Lake Ohrid record provides important, independent tie points for the age–depth model of the DEEP site sequence, which is a prerequisite for palaeoclimatic and palaeoenvironmental reconstructions. Furthermore, this age–depth model will help to improve and re-evaluate the chronology of other, both undated and dated tephra layers from other records. Thus, the Lake Ohrid record may potentially become the template for the central Mediterranean tephrostratigraphy, especially for the hitherto poorly known and explored lower Middle Pleistocene period.

Ancient lakes, such as lakes Ohrid and Prespa on the Balkan Peninsula, have become model systems for studying the link between geological and biotic evolution. Recently, the scientific deep-drilling project Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) was initiated to better understand the environmental, climatic, and limnological evolution of the lake. It revealed that Lake Ohrid experienced a number of environmental disturbances during its ca. 2.0 million year long history. These are comprised of disturbances that lasted over longer periods of time (“press events”) such as glacial–interglacial cycles and Heinrich events, as well as sudden and short disturbances (“pulse events”) like the deposition of landslides, earthquakes, and volcanic ash depositions. The latter includes one of the most severe volcanic episodes during the Late Pleistocene: the eruption of the Campanian Ignimbrite (known as Y-5 marine tephra layer) from the Campi Flegrei caldera, dated to 39.6 ± 0.1 thousand years ago. The event is recorded by the deposition of a ca. 15 cm thick tephra layer in sediment cores of lakes Ohrid (DEEP-5045-1) and Prespa (Co1204). Coincidently, this pulse event is superimposed by the Heinrich H4 event, 40.4–38.4 thousand years ago. In the current paper, diatoms were used as proxies to compare the responses of these lakes to the Y-5 (pulse) and the H4 (press) disturbances. Based on stratigraphically constrained incremental sum of squares cluster (CONISS) and unconstrained Partitioning Around Medoids (PAM) analyses, we found little evidence that diatom community compositions in either lake responded to the H4 event. However, the Y-5 influx caused clear and rapid diatom community changes. After the initial response, community compositions in Lake Ohrid and, to a lesser extent, in Lake Prespa slowly returned to their quasi pre-disturbance state. Moreover, there is no evidence for disturbance-related extinction events. The combined evidence from these findings suggests that lakes Ohrid and Prespa likely did not experience regime shifts. It is therefore concluded that both lakes show resilience to environmental disturbance. However, it seems that Lake Ohrid is more resilient than Lake Prespa, as the recovery of diatom communities is more pronounced and its estimated recovery time is only ca. 1100 years vs. ca. 4000 years in Lake Prespa. The reasons for the differential responses remain largely unknown, but differences in geology, lake age, limnology, and intrinsic parameters of the diatom proxies may play an important role.

Lake Ohrid (Macedonia, Albania) is thought to be more than 1.2 million years old and host more than 300 endemic species. As a target of the International Continental scientific Drilling Program (ICDP), a successful deep drilling campaign was carried out within the scope of the Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) project in 2013. Here, we present lithological, sedimentological, and (bio-)geochemical data from the upper 247.8 m composite depth of the overall 569 m long DEEP site sediment succession from the central part of the lake. According to an age model, which is based on 11 tephra layers (first-order tie points) and on tuning of bio-geochemical proxy data to orbital parameters (second-order tie points), the analyzed sediment sequence covers the last 637 kyr. The DEEP site sediment succession consists of hemipelagic sediments, which are interspersed by several tephra layers and infrequent, thin (< 5 cm) mass wasting deposits. The hemipelagic sediments can be classified into three different lithotypes. Lithotype 1 and 2 deposits comprise calcareous and slightly calcareous silty clay and are predominantly attributed to interglacial periods with high primary productivity in the lake during summer and reduced mixing during winter. The data suggest that high ion and nutrient concentrations in the lake water promoted calcite precipitation and diatom growth in the epilimnion during MIS15, 13, and 5. Following a strong primary productivity, highest interglacial temperatures can be reported for marine isotope stages (MIS) 11 and 5, whereas MIS15, 13, 9, and 7 were comparably cooler. Lithotype 3 deposits consist of clastic, silty clayey material and predominantly represent glacial periods with low primary productivity during summer and longer and intensified mixing during winter. The data imply that the most severe glacial conditions at Lake Ohrid persisted during MIS16, 12, 10, and 6, whereas somewhat warmer temperatures can be inferred for MIS14, 8, 4, and 2. Interglacial-like conditions occurred during parts of MIS14 and 8.

Gamma ray (GR) fluctuations and potassium (K) values from downhole logging data obtained in the sediments of Lake Ohrid from 0 to 240 m below lake floor (b.l.f). correlate with fluctuations in δ18O values from the global benthic isotope stack LR04 (Lisiecki and Raymo, 2005). GR and K values are considered a reliable proxy to depict glacial–interglacial cycles, with high clastic input during cold and/or drier periods and high carbonate precipitation during warm and/or humid periods at Lake Ohrid. Spectral analysis was applied to investigate the climate signal and evolution over the length of the borehole. Linking downhole logging data with orbital cycles was used to estimate sedimentation rates and the effect of compaction was compensated for. Sedimentation rates increase on average by 14 % after decompaction of the sediment layers and the mean sedimentation rates shift from 45 cm kyr−1 between 0 and 110 m to 30 cm kyr−1 from 110 to 240 m b.l.f. Tuning of minima and maxima of gamma ray and potassium values versus LR04 extrema, in combination with eight independent tephrostratigraphical tie points, allows establishing of a robust age model for the downhole logging data over the past 630 kyr.

Leaching experiments were carried out on fresh ash samples from Popocatépetl 2012, Etna 2011, and Etna 2012 eruptions, in order to investigate the release of compounds in both double-deionized and lake (Lake Ohrid, FYR of Macedonia) waters. The experiments were carried out using different grain sizes and variable stirring times (from 30 min to 7 days). Results were discussed in the light of changing pH and release of compounds for the different leachates. In particular, Etna samples induced alkalinization, and Popocatépetl samples induced acidification of the corresponding leachates. The release of different elements does not show correlation with the stirring time, with the measured maximum concentrations reached in the first hours of washing. General inverse correlation with grain size was observed only for Na+, K+, Cl−, Ca2+, Mg2+, SO42−, and Mn2+, while the other analysed elements show a complex, scattering relationship with grain size. Geochemical modelling highlights leachates' saturation only for F and Si, with Popocatépetl samples sometimes showing saturation in Fe. The analysed leachates are classified as undrinkable for humans on the basis of European laws, due to excess in F−, Mn2+, Fe, and SO42− (the latter only for Popocatépetl samples). Finally, the Etna 2012 and Popocatépetl leachates were used for toxicity experiments on living biota (Xenopus laevis). They are mildly toxic, and no significant differences exist between the toxic profiles of the two leachates. In particular, no significant embryo mortality was observed; while even at high dilutions, the leachates produced more than 20 % of malformed larvae.

Lakes Prespa and Ohrid, in the Balkan region, are considered to be amongst the oldest lakes in Europe. Both lakes are hydraulically connected via karst aquifers. From Lake Ohrid, several sediment cores up to 15 m long have been studied over the last few years. Here, we document the first long sediment record from nearby Lake Prespa to clarify the influence of Lake Prespa on Lake Ohrid and the environmental history of the region. Radiocarbon dating and dated tephra layers provide robust age control and indicate that the 10.5 m long sediment record from Lake Prespa reaches back to 48 ka. Glacial sedimentation is characterized by low organic matter content and absence of carbonates in the sediments, which indicate oligotrophic conditions in both lakes. Holocene sedimentation is characterized by particularly high carbonate content in Lake Ohrid and by particularly high organic matter content in Lake Prespa, which indicates a shift towards more mesotrophic conditions in the latter. Long-term environmental change and short-term events, such as related to the Heinrich events during the Pleistocene or the 8.2 ka cooling event during the Holocene, are well recorded in both lakes, but are only evident in certain proxies. The comparison of the sediment cores from both lakes indicates that environmental change affects particularly the trophic state of Lake Prespa due to its lower volume and water depth.

Four cores from Balkans lakes Ohrid and Prespa were examined for recognition of tephra layers and cryptotephras, and the results presented along with the review of data from other two already published cores from Lake Ohrid. The six cores provide a previously unrealised tephrostratigraphic framework of the two lakes, and provide a new tephrostratigraphic profile (composite) for the Balkans, which spans from the end of the Middle Pleistocene to the AD 472. A total of 12 tephra layers and cryptotephras were recognised in the cores. One is of Middle Pleistocene age (131 ka) and correlated to the marine tephra layer P-11 from Pantelleria Island. Eight volcanic layers are Upper Pleistocene in age, and encompass the period between ca. 107 ka and ca. 31 ka. This interval contains some of the main regional volcanic markers of the central Mediterranean area, including X-6, X-5, Y-5 and Y-3 tephra layers. The other layers of this interval have been related to the marine tephra layers C20, Y-6 and C10, while one was for the first time recognised in distal areas and correlated to the Taurano eruption of probable Vesuvian origin. Three cryptotephras were of Holocene age. Two of which have been correlated to Mercato and AD 472 eruptions of Somma-Vesuvius, while the third has been correlated to the FL eruption from Mount Etna. These recognitions provide a link of the Ohrid and Prespa lacustrine successions to other archives of the central Mediterranean area, like south Adriatic, Ionian, and south Tyrrhenian seas, lakes of southern Italy (Lago Grande di Monticchio, Pantano di San Gregorio Magno and Lago di Pergusa) and Balkans (Lake Shkodra).