Research Article
Print
Research Article
Paleoclimate changes and ecosystem responses of the Bulgarian Black Sea zone during the last 26000 years
expand article infoMariana Filipova-Marinova, Danail Pavlov§, Krasimira Slavova|
‡ Museum of Natural History Varna, Varna, Bulgaria
§ Medical University of Varna, Varna, Bulgaria
| Institute of Oceanology, Bulgarian Academy of Sciences, Varna, Bulgaria
Open Access

Abstract

Multi-proxy analysis (spore-pollen, dinoflagellate cysts, other non-pollen palynomorphs (NPPs), radiocarbon dating and lithology) was performed on marine sediments from three new cores retrieved during the two cruise expeditions on board the Research Vessel “Akademik” in 2009 and 2011. The Varna transect comprises three cores extending from the outer shelf, continental slope and deep-water zone. The record spans the last 26000 years (all ages obtained in this study are given in calendar years BP (cal. yrs BP)). The pollen record reveals the spreading of steppe vegetation dominated by Artemisia and Chenopodiaceae, suggesting cold and dry environments during the Late Pleniglacial – Oldest Dryas (25903–15612 cal. yrs BP). Stands of Pinus and Quercus reflect warming/humidity increase during the melting pulses (19.2–14.5 cal. ka BP) and the Late Glacial interstadials Bølling and Allerød. The Younger Dryas (13257–11788 cal. yrs BP) coldest and driest environments are clearly demonstrated by the maximum relative abundance of Artemisia and Chenopodiaceae. During the Early Holocene (Preboreal and Boreal chronozones, 11788–8004 cal. yrs BP), Quercus appeared as a pioneer species and, along with other temperate deciduous arboreal taxa, formed open deciduous forests as a response to the increased temperature. The rapid expansion of these taxa indicates that they survived in Glacial refugia in the coastal mountains. During the Atlantic chronozone (8004–5483 cal. yrs BP), optimal climate conditions (high humidity and increased mean annual temperatures) stimulated the establishment of species-rich mixed temperate deciduous forests. During the Subboreal chronozone (5483–2837 cal. yrs BP), Carpinus betulus and Fagus expanded simultaneously and became more important components of mixed oak forests and probably also formed separate communities. During the Subatlantic chronozone (2837 cal. yrs BP to pre-industrial time), climate-driven changes (an increase of humidity and a cooling of the climate) appear to be the main drivers of the specific vegetation succession expressed by increased abundance of Alnus, Fraxinus excelsior and Salix along with lianas, suggesting formation of flooded riparian forests (so called ‘Longoz’) lining the river valleys along the Black Sea coast. The first indicators of farming and other human activities have been recorded since 7074 cal. yrs BP. The dinoflagellate cyst (dinocyst) assemblages have been analysed to assess the changes in the Black Sea environment over the last 26000 years in terms of fluctuation in paleoproduction and surface water conditions related to changes in climate, freshwater input and Mediterranean water intrusion. Two major dinocyst assemblages were distinguished: one dominated by stenohaline freshwater/brackish-water species and a successive one dominated by euryhaline marine species. The changes in the composition of the assemblages occurred at 7668 cal. yrs BP. The abrupt decrease of stenohaline freshwater/brackish-water species Pyxidinopsis psilata and Spiniferites cruciformis was followed upwards by a gradual increase in euryhaline marine species, such as Lingulodinium machaerophorum, Spiniferites belerius, S. bentorii and acritarch Cymatiosphaera globulosa. The first occurrence of euryhaline marine species took place synchronously with the onset of sapropel deposition. Modern marine conditions were established after 6417 cal yrs BP when an abundance of Mediterranean-related species, such as Operculodinium centrocarpum and Spiniferites mirabilis, along with other heterotrophic species, occurred. After the stable cold and dry environment during the Last Glacial Maximum, the phytoplankton record of core AKAD 11-17 shows that Pediastrum boryanum var. boryanum has a cyclical abundance associated with the deposition of four red-brown clay layers between 19.2 and 14.5 cal. ka BP. This event is associated with the major melting phase of European Ice drained by the Danube and Dnieper Rivers in response to climate warming observed after the end of the Last Glacial Maximum. During the Early Holocene, P. psilata, characterised by a preference to warmer temperatures, demonstrates its ecological optimum for growth concerning SST reaching maximum relative abundance of 94% between 11072 and 8638 cal. yrs BP. This maximum was interrupted by an abrupt significant short-term decrease in the relative abundance of P. psilata centred between 8500 and 8300 cal. yrs BP reflecting cold conditions similar to those of Younger Dryas. This finding, also confirmed by the rapid significant decrease of arboreal pollen, particularly of Quercus in the same studied core, is considered a regional expression of the well-known ‘8.2 ka cold event’ which is commonly linked to a meltwater-related perturbation of the Atlantic Meridional Overturning Circulation (AMOC) and associated collapse of oceanic northward heat transport. Our fossil pollen and dinocyst data confirm that the high amplitude temperature anomaly associated with ‘the 8.2 ka cold event’ may have also occurred in south-eastern Europe, at lower latitudes of the western Black Sea coastal area, most probably due to atmospheric transition and/or river discharge.

Key words

Dinoflagellate cyst, non-pollen palynomorphs, radiocarbon dating, spore-pollen analysis

Introduction

The Black Sea, as an almost isolated marginal sea, is particularly sensitive to paleoenvironmental changes and, therefore, Black Sea sediments provide an excellent opportunity for high-resolution studies of past climatic, vegetation, human activity and hydrological changes in the catchment (Bahr et al. 2005). In contrast to shelf records which are affected by erosion during lowstands, pollen and dinocyst records from the continental slope and deep water Black Sea cores are of particular interest as they can provide almost uninterrupted sequences covering the Late Pleistocene and Holocene and can be used to obtain an independent record of regional climate change and land-based interpretation of the reconstructed vegetation and paleohydrological regime. These reconstructions are able to describe the interaction between climate and vegetation and also to clarify the role of coastline and other geomorphological changes, salinity and impacts of human activities in the Black Sea region (Cordova et al. 2009).

The Black Sea sediments have been intensively investigated by multi-proxy analysis during the last five decades. The biostratigraphic investigations of Quaternary marine sediments taken by the Scientific-Research vessels ``Atlantis-2” and “Glomar-Challenger’’ established a baseline chronostratigraphy. Palynological investigations of sediments from the deep-water zone allowed Traverse (1974, 1978a, 1978b), Koreneva and Kartashova (1978) and Koreneva (1980) to outline the stratigraphy of sediments, vegetation dynamics and climate changes along the Black Sea coast from the end of the Pliocene through the Pleistocene at a very low resolution of tens of thousands of years. The first detailed marinopalynological investigations of the western Black Sea shelf are those of Roman (1974), Bozilova et al. (1979) and Komarov et al. (1979). Based on marinopalynological data of the western Black Sea sediments, the climatic changes during the Late Glacial and Holocene were estimated and considered to be the main drivers for vegetation changes along the Bulgarian Black Sea coast (Shimkus et al. 1977; Komarov et al. 1978; Bozilova et al. 1979; Chernyishova 1980; Filipova and Dimitrov 1987; Atanassova 1990, 1995, 1999, 2005; Atanassova and Bozilova 1992; Mudie et al. 2007). Palaeoecological changes during the Quaternary and possible sea-level fluctuations during the Holocene were traced by Filipova-Marinova (2003a, 2003b, 2007); Filipova-Marinova and Christova (2004), Filipova-Marinova et al. (2004) and Hiscott et al. (2007). Mudie et al. (2007) showed that organic-walled microfossils including pollen, spores and dinocysts are well-preserved and abundant in deep-water Pleistocene and Holocene sediments. These authors presented pollen assemblages for the last 33000 years and reported the first high resolution Holocene pollen influx data for the SW Black Sea shelf. Filipova-Marinova et al. (2013) reported the first high-resolution palynostratigraphy of Late Quaternary sediments and the chronologically defined vegetation stages with their specific features, particularly a short-term cooling of the Holocene climate, associated with the ‘8.2 ka cold event’ identified for the first time in marine records from the central Bulgarian Black Sea area. Vegetation and environmental dynamics in northern Anatolia during the last 18000 years are reconstructed using multi-proxy records from the southern Black Sea shelf (Shumilovskikh et al. 2012). Mudie et al. (2002) showed the close correlation between marine pollen assemblage zones in the Marmara and southern Black Sea with lakes in northern Turkey and Bulgaria. Quantitative paleovegetation reconstruction and correlation of pollen data of 99 sequences from the Black, Marmara and Azov Seas are presented by Cordova et al. (2009). First palynostratigraphy of western Black Sea sediments was proposed by Filipova-Marinova (2006a), based on 12 representative cores.

The pioneer work of Wall et al. (1973) enabled the use of dinocysts as potential paleoceanographic proxies. Since the 1980s, numerous studies from the Black Sea permitted the definition of distribution patterns in regional levels and have led to determination of the relationship between dinocyst assemblages and sea surface water conditions (De Vernal et al. 1994; Marret 1994; Matthiessen 1995; Mudie and Harland 1996; Mudie et al. 2001; Mudie et al. 2002; Marret and Zonneveld 2003; Marret et al. 2009; Mertens et al. 2012). The studies of van Geel (2001), Brenner (2001), Mudie et al. (2002) and Marret et al. (2009) describe non-pollen palynomorphs (fungal spores and remains of hyphae) as important markers of salinity, nutrient loading and human activity, including ballast discharge, farming and soil erosion. Mudie et al. (2002) and Marret et al. (2009) consider fungal remains as an index of terrigenous sediment influx and transport by large rivers. The Late Quaternary history of connection of the Black Sea to the eastern Mediterranean, especially the timing and conditions during the Holocene transition of the Black Sea from freshwater/brackish-water to marine has been intensively debated: Ryan et al. (1997, 2003), Aksu et al. (2002), Major et al. (2002), Bahr et al. (2005, 2006), Hiscott et al. (2007), Yanko-Hombach (2007); Londeix et al. (2009), Marret et al. (2009), Verleye et al. (2009) and Herrle et al. (2018). In addition to the Holocene connection with the Marmara Sea, the Black Sea experienced a period of connection to the Caspian Sea when meltwater from Scandinavian ice sheets raised the water level of the Caspian Sea over the Manych depression. This period is documented in the Black Sea by the deposition of a series reddish-brown clay layers (Major et al. 2002; Bahr et al. 2005, 2008).

Recent studies in the western Black Sea are focused mainly on palaeoecological changes during the Late Glacial and Holocene. There is a lack of uninterrupted Late Quaternary sediments from the northern Bulgarian Black Sea area adjacent to Varna. Therefore, in order to obtain appropriate records which would allow a better more detailed description of the pollen and dinocyst stratigraphy and more precise palaeoecological reconstructions of the Bulgarian sector of the Black Sea, two expeditions by the Research Vessel “Akademik” in 2009 and 2011 were carried out. A total of 105 new samples were taken for multi-proxy analysis of sediments from three representative cores: AKAD 11-17 (deep-water zone, water depth: 1805 m, core length: 228.5 cm), AKAD 09-10 (zone of the continental slope, water depth: 1000 m, core length: 240 cm) and AKAD 09-15 (outer shelf zone, water depth 164 m, core length 377 cm).

The aim of this study is to establish the palaeoclimatic, palaeohydrological and environmental dynamics of the Bulgarian Black Sea coastal area during the Late Pleistocene and Holocene, as well as to evaluate the timing and extent of the passage between various events previously described in other studies, based on multi-proxy analysis.

Regional setting

Present-day Oceanography of the Black Sea

The Black Sea (Fig. 1) is located between south-eastern Europe and Asia. It has an area of approximately 432,000 km2 and a maximum water depth of 2258 m. It is the largest semi-enclosed inland basin in the world, where the deep-waters and surface waters are not mixed (Nikishin et al. 2003; Murray et al. 2007). The oxygenated surface layer overlies an anoxic deeper layer with elevated H2S concentrations up to 380 µM (Murray et al. 2007). The Black Sea has an extremely large drainage basin of more than two million km2, collecting the water from almost all the European countries, except the westernmost ones. The hydrologic configuration of the Black Sea is controlled by basin bathymetry and also by fluvial inputs from the discharge of the largest European rivers including the Danube, Dnieper, Dniester and the Southern Bug. Ukrainian rivers Dnieper, Southern Bug and Dniester contribute about 65 km3/yr and the Danube River is with a mean water discharge of about 200 km3/yr (Panin 2008). These inputs, plus those from smaller rivers from the Bulgarian and north Turkish margin, create a surface low salinity layer across the whole of the Black Sea. The net freshwater input to the Black Sea is about 300 km3/yr and derives from the river input plus precipitation. The Black Sea is presently connected to the Global Ocean through the narrow and shallow Bosphorus and the Dardanelles Straits that limit the salinity and oxygen provided to the Black Sea (Özsoy et al. 1994; Kerey et al. 2004). In addition to this marine connection, the Black Sea currently also has connection with the freshwater Sea of Azov through the narrow Strait of Kerch that has a constant outflow of waters to the Black Sea (Kosarev et al. 2007). The water exchange occurs between the Black and Marmara Seas via the Bosphorus Strait as a two-layer water flow (Özsoy et al. 1994). About 600 km3/yr cooler (5–15 °C) and less saline (18–20‰) water mass flows out of the Black Sea through the Bosphorus Strait and, from there, through the Sea of Marmara into the Dardanelles Strait as a surface outflow into the Mediterranean (Özsoy et al. 1995). This surface outflow is compensated by about 300 km3/yr of warmer (15–20 °C) and more saline (38–39‰) mixed Marmara and Aegean Sea waters that flow into the Black Sea through the Bosphorus as a deep inflow. This produces a density-stratified water column due to the largely varying salinities entering the region. The surface layer 0–50 m is well oxygenated, while the deeper water layer 100–2243 m is highly anoxic and rich in sulphides. The approximate permanent location of the halocline is 50 to 200 m (Murray et al. 2007). Today, the surface water circulation in the Black Sea consists of two large cyclonic (counterclockwise) central gyres that define the eastern and western basins. The gyres are bounded by the wind-driven ‘Rim current’ that flows along the edge of the continental shelf and above the continental slope around the whole basin (Oguz et al. 1993). Data from autonomous profiling floats have shown currents typically have a velocity of 15 cm/s at 200 m depth along the ‘Rim Current’ jet around the basin. At depths of 750 m and 1500 m, current velocities of 5 cm/s have been recorded with the deeper current closely following the topography along the southern margin of the Black Sea. Outside the Rim Current, numerous quasi-permanent coastal eddies are formed as a result of upwelling around the coastal apron and “wind curl” mechanisms.

Figure 1.

Location of the studied cores (black triangle) AKAD 11-17, AKAD 09-10 and AKAD 09-15 and four reference cores (black square) discussed.

Sedimentology of the Black Sea

The most recent pelagic sediment layers in the Black Sea can be divided into three units:

  • Unit I is represented by micro-laminated sediments, rich in plankton-derived carbonates finely laminated coccolith-bearing ooze. The calcareous material is derived from marine coccolithophorid Emiliania huxleyi, which have formed part of the Black Sea plankton since 2700 years BP (Jones and Gagnon 1994). Blooms of E. huxleyi still occur in the Black Sea each year during late spring and summer. This unit was deposited in oxygen-depleted bottom waters.
  • Unit II sediments contained micro-laminated sapropels, deposited under anoxic marine conditions between 2720 ± 160 14 C yrs BP and 7540 ± 130 14C yrs BP. The base of the sapropel has been dated to 7540 ± 130 14Cyrs BP in cores running in depth from 400 to 2200 m suggesting that anoxia developed at the same time throughout the interior Black Sea (Jones and Gagnon 1994). The onset of Unit II is characterised by the occurrence of finely-laminated layers rich in aragonite crystals (Ross and Degens 1974; Jones and Gagnon 1994; Soulet et al. 2011) and by a sharp increase in Total Organic Carbon (Bahr et al. 2008).
  • Unit III sediments are older than 7540 ± 130 14Cyrs BP and were deposited when the Black Sea was a freshwater to brackish lake and are characterised by a mix of organic-poor clays and silts (Izdar and Ergün 1991; Hay et al. 1991). Unit III sediments have organic contents < 1%. Sediments consist of homogenous mostly centimetre-scale laminated muddy clay deposited under lacustrine conditions.

Geomorphology of the Bulgarian sector of the Black Sea

The large continental shelf in the north-western Black Sea narrows in the southerly direction. On the basis of the relief, shape, time of formation, character and speed of sedimentological processes, three geomorphological zones can be outlined in the western Black Sea shelf: littoral or inner, central and peripheral or outer (Dimitrov 1979). The littoral zone is considered to be of Holocene age. It extends from the coast to a depth of 20 to 50 m of the central part of the Bulgarian Black Sea coastal area. Active wave impact, erosion and accumulation are characteristic processes for this zone (Khrischev 1984). The littoral zone is separated from the central one by a depression 17–20 m deep on the northern Black Sea shelf and 65–70 m deep on the southern Bulgarian Black Sea shelf. The central zone lies between 50 and 72 m in depth. Within it, three subzones run parallel to the coast: an inner depression, an area of depositional bars and a depositional plain. This zone experiences a high sedimentation rate, typically about 2.5 m/kyr (Khrischev 1984). The peripheral zone extends to depths between 90 m and 120 m. It is subdivided into an outer depression and an area of barrier bars and its low sedimentation rate is due to sediment removal by strong bottom currents (Dimitrov 1979). The north-western Black Sea shelf is covered only by a very thin blanket of Holocene sediments and, in some parts, the Holocene layers are completely absent due to the strong (> 50 cm.s-1) cyclonic ‘Rim current’ (Oguz and Besiktepe 1999) which transports most of the suspended sediments alongside the coast to the south (Panin еt al. 1999).

According to the general morphostructural plan of the Black Sea deep-water basin, the continental slope covers 25% of its surface. The transition of the shelf to the continental slope is gentle and has a convex-up profile. The continental slope of the Bulgarian Black Sea zone is characterised by deeply-indented relief including land-sliding complexes, fault slopes, ledges and submerged valleys and canyons. Nine systems of submerged valleys are established in the area (Alexiev 2002). The transition from the continental slope to the deep basin increases in the western Black Sea from 1100 m in the north to 2000 m in the south.

The continental foot is formed by the confluence of the sedimentation materials of the submerged delta valleys. The gentle transition of the steep continental slope to the abyssal plain is accomplished by its slightly undulating plain surface. The formation of the modern shape of the continental slope took place mainly during the Pleistocene. The large input of terrestrial sediments that are the main constructive material for the continental foot was disrupted by the breaking-off access of the coastal rivers to the outer zone of the continental slope during the Holocene. The limit between the continental foot and the abyssal plain is difficult to be located. However, the isobaths 2000–2100 m localise a typical abyssal plain which declines slightly towards the deepest part of the Black Sea Basin. The abyssal plain is the earliest formed morphological element of the Black Sea (Alexiev 2002).

Climate

According to Velev (2002), the Bulgarian Black Sea coast belongs to the Continental-Mediterranean climatic area and is influenced by three different climatic regimes. Climate in the northern part and the Eastern Danube lowland is affected by strong continental influences. Prevailing winds are northeasterly and mean annual precipitation is about 450–500 mm, with a maximum in June and a minimum in February. Mean January temperature is around 0 °C, dropping to -2 °C inland. In the south, the climate is transitional Mediterranean. Mean annual precipitation is estimated at about 500–600 mm, with rainfall mostly in the autumn-winter seasons. Mean January temperature is 2–3 °C and, in July, it is 22 °C. The dry summer period lasts from July to September. Winds blow mostly from the southeast and rarely from the northeast. Mean annual precipitation over the mountain areas of the coast (500 m a.s.l.) is about 600–1000 mm (Velev 2002). The western part of the coast is under the influence of the humid conditions of central Europe with precipitation more than 1000 mm/year.

Recent vegetation

The Bulgarian Black Sea coast covers a narrow strip of land located to the west of the Black Sea coastline. It is 375 km long and 30 to 50 km wide and includes Southern Dobrudzha, the Eastern Stara Planina Mountains (Balkan Range), Burgas Plain and the Eastern Strandzha Mountains. Climate controls and local topography play a dominant role in determining the pattern of highly-varied natural vegetation. This area is considered as a major pollen source area for the investigated core sediments with consideration of wind pattern, river input and gyre systems in the Black Sea. In addition to the vegetation distribution map, simplified characteristics of vegetation types are presented by Bondev (1991). According to this author, the study area falls within the Black Sea region of the Euxinian province of the European deciduous forest. The vegetation is represented by steppe vegetation, different types of temperate deciduous forests, reed, psammophytic and halophytic vegetation. Steppe vegetation of natural origin dominated by Stipa capillata L., Agropyron brandzae Pantu et Solac., Koeleria brevis Stev., Stipa lessingiana Trin. et Rupr., Artemisia lerchiana Weber, Adonis vernalis L., Adonis volgensis Steven ex DC and Paeonia tenuifolia L., amongst others, is preserved only in the northern Bulgarian Black Sea coastal area, in the South Dobrudzha region (Cape Shabla and Cape Kaliakra areas). Xerothermic forest communities dominated by Quercus cerris L. and Q. frainetto Ten. are widespread in the eastern Balkan Range area (Eastern Stara Planina Mts.) (Cape Emine area). Restricted areas on northern slopes and lower moisture ravines are occupied by stands of southeuxinian taxa such as Quercus polycarpa Schur. and Fagus orientalis Lipsky along with Carpinus betulus L., Acer campestre L., Q. cerris and Tilia tomentosa Moench. Xeromesophytic communities comprised of Q. polycarpa and Carpinus betulus are spread in lowland sites and on hilltops, as well as in the Strandzha Mountains. Mediterranean elements, such as Quercus pubescens Willd., Carpinus orientalis Mill., Fraxinus ornus L., Phillyrea latifolia L., Celtis australis L. and Colutea arborescens L. are distributed along the southern Black Sea coast. Relic southeuxinian forests of Fagus orientalis, with an undergrowth of evergreen shrubs (Rhododendron ponticum L., Ilex aquifolium L. and Daphne pontica L.) cluster the more humid ravines of the Strandzha Mts. (Veleka River, Rezovska River and Sozopol area). Riparian forests (so called Longoz) line rivers flowing into the Black Sea and coastal lakes. The main components of these forests are Fraxinus oxycarpa Willd., Ulmus minor Mill., Carpinus betulus, Quercus pedunculiflora C. Koch and Alnus glutinosa (L.) Gaerth. The most characteristic for these periodically flooded forests is the presence of lianas such as: Hedera helix L., Periploca graeca L., Clematis vitalba L., Vitis vinifera L. and Smilax excelsa L. Reed vegetation represented by communities dominated by Phragmites australis (Cav.) Trin. ex Steud., Typha angustifolia L., T. latifolia L. and Schoenoplectus lacustris (L.) Palla are spread along the rivers and the periphery of coastal lakes (Kochev and Jordanov 1981). Psammophytic vegetation includes communities mainly of Leymus racemosus (Lam.) Tzvel. ssp. sabulosus (Bieb.) Tzvel., Ammophilla arenaria (L.) Link, Centaurea arenaria Bieb. ex Willd., Galilea mucronata (L.) Parl and shrub communities with Cionura erecta (L.) Grsb. growing on the sandy beaches and dunes. Halophytic vegetation presented by communities of Salicornia europaea has limited distribution in habitats with high salinity.

Material and methods

Coring and sampling

Sediments pertinent to this study were collected during expeditions in 2009 and 2011 onboard the Research Vessel Akademik owned by the Institute of Oceanology of the Bulgarian Academy of Sciences. The two cruise expeditions recovered a series of sediment cores on a number of shallow-to-deep transects from the Bulgarian Black Sea area. The Varna transect consists of three cores taken from the shelf, continental slope and deep-water zone of the Black Sea (Fig. 1). Sediment core Akad 09-10 was collected with a Vibracore device, because of the feature of the Black Sea consolidated shelf sediments and cores Akad 09-15 and Akad 11-17 were collected with Gravity Corer device. The both devices have the same dimensions of the diameter of a metal tube and a length of metal tube, respectively Ø12 cm and 4 m. The recovered sediment cores were cut in 1-m sections on board and carried to the lab. The sediment was pulled from the steel tube in its plastic steeve, split lengthwise in two half cylinders, photographed, visually lithologically described and sampled. Samples of 1 cm3 were removed from the centre of the split section to avoid contamination with any younger sediment smeared downwards during the insertion of the sediment into the plastic steeve. The cores were lithologically described (Fig. 2) and sampled at every 10 cm on the board of the Research Vessel “Akademik” for further analyses.

Figure 2.

Lithological and geochronological correlations of the studied cores AKAD 11-17, AKAD 09-10 and AKAD 09-15.

Core Akad 11-17 (42°51'13.50"N, 29°01'08.50"E) was recovered from a water depth of 1805 m in the Black Sea deep-water zone (Fig. 1). The investigated length of the core is 229 cm (Fig. 2). From the core base up to the depth of 59 cm, the core sediments are typical lacustrine clay sediments with alternately light and dark layers, deposited under freshwater to brackish conditions (Unit III). From 229 to 133 сm, red-brown clay with grey clay is deposited in several pulses. These red-brown clay layers can be found in the sediments from the whole western Black Sea Basin. In the interval 155–153 cm were described six sand silt lamination layers. Above the red-brown clay was described an alternation of light grey clay with dark grey clay with a band of sulphides. At 69, 72, 72.5, 76, 80, 84, 90, 109–110 сm was described 1 mm thin sand silt lamination. The marine units II and I contain organic-rich microlaminated sapropel sediments and carbonate-rich finely laminated coccolith ooze and are ascertained in the core intervals from 59–27 cm and 27–0 cm to the top, respectively (Fig. 2). The base of the sapropel deposition in the core Akad 11-17 was described at 59 сm. The base of the coccolith ooze was described at 27 cm.

Core Akad 09-10 (42°54.8'N, 28°45.6'E) was recovered from a water depth of 1000 m on the Bulgarian Black Sea continental slope (Fig. 1). The investigated length of the core is 242 cm (Fig. 2). From 242 to 122 сm, sediments are typically lacustrine (Unit III) with homogeneous light to dark grey clay deposited under freshwater to brackish conditions. An important lithological feature of Unit III is the deposition of red-brown clay layers within the range from 236 to 216 сm. The onset of marine sapropel sediment deposition is lithologically determined at 126 cm (Unit II). The sapropel mud is finely laminated, more firm, dark green-grey with three white laminae in the interval from 122 to 123 cm in the lower part of marine Unit II. The boundary between Unit II and Unit I is delimited by the presence of coccolith ooze as well as the stable appearance of Еmiliania huxleyi at 50 cm.

Core Akad 09-15 (42°58.628'N, 28°33.147'E) was recovered from a water depth of 164 m on the Bulgarian Black Sea shelf (Fig. 1). The investigated length of the core is 380 cm. In this core, only lacustrine sediments (Unit III, lutite) were established (Fig. 2). From 380 to 8 cm, sediments are light to dark grey clay with rare whole disarticulated small freshwater Dreissena sp. and more fragmented shells. At 180 cm is described a layer of shell hash. In the interval from 8 to 0 cm, sediments are of grey sand.

Pollen analysis and zonation

A total of 105 samples were selected for palynological, dinoflagellate cyst and other non-pollen palynomorph analyses. Each sample consists of 1 cm3 of wet sediment. The sampling interval was 10 cm. Sediments characterised by much higher sedimentation rates were sampled with varying resolution (Figs 28). All samples were processed according to the standard procedure of Faegri and Iversen (1989). The laboratory technique includes treatment with hot 10% hydrochloric acid (HCl), cold 40% hydrofluoric acid (HF), solution of zinc chloride (ZnCl2), 10% potassium hydroxide (KOH), glacial acetic acid) CH3COOH), 2 min acetolysis and glacial CH3COOH and ethanol (C2H5OH). The removal of mineral components was performed using sodium pyrophosphate and hydrofluoric acid (Birks and Birks 1980). The obtained suspension is stored in glycerine and then used for microscope analysis. A minimum of 500 up to 1000 identifiable pollen grains from terrestrial plants and a minimum of 100 up to 300 dinoflagellate cysts and non-pollen palynomorphs were counted in each sample to ensure statistical significance.

Pollen types were identified using the reference collection of modern pollen types of the Museum of Natural History of Varna, keys in Erdtman et al. (1961), Beug (1961, 2004), Moore and Webb (1978), Faegri and Iversen (1989) and photographs of Reille (1992, 1995). Dinoflagellate cysts taxonomy is based on Wall et al. (1973) and Marret et al. (2004). Other non-pollen palynomorphs (NPPs) were identified by the keys in van Geel et al. (1981), van Geel (2001) and Jankovska and Komarek (2000). The total pollen sum (PS) used for percentage calculations of the individual percentages includes arboreal pollen (AP) and non-arboreal pollen (NAP). Pollen of aquatics, spores, dinoflagellate cysts and other NPPs are excluded from the PS. Their presence was expressed as percentage of the total PS (van Geel et al. 1981). The software TILIA v.1.17.16 (Grimm 2011) was used for pollen and dinocyst percentage calculations and construction of spore-pollen and dinocyst diagrams. A 10× exaggeration of the horizontal scale was used to show changes of low-percentage taxa. The spore-pollen diagrams were subdivided into fifteen local pollen assemblage zones (LPAZ) (Figs 3, 5, 7) and two local dinoflagellate assemblage zones (LDAZ) (Figs 4, 6, 8) using visual inspection of the main changes in taxa composition (described in details in Table 2). The cluster analysis programme CONISS (Grimm 1987) is applied for more precise zonation (Figs 3, 5, 7). These LPAZ are correlated with the established regional pollen assemblage zones (RPAZ) IV to IX after Filipova-Marinova (2006a) (Table 3). In the spore-pollen diagrams, taxa are ordered by their plant functional type (Prentice et al. 1996).

Figure 3.

Percentage spore-pollen diagram of Core AKAD 11-17 (Black Sea deep-water zone).

Figure 4.

Percentage diagram of dinocysts, other algae and non-pollen palynomorphs (NPPs) in Core AKAD 11-17 (Black Sea deep-water zone).

Figure 5.

Percentage spore-pollen diagram of Core AKAD 09-10 (Black Sea continental slope).

Figure 6.

Percentage diagram of dinocysts, other algae and non-pollen palynomorphs (NPPs) in Core AKAD 09-10 (Black Sea continental slope).

Figure 7.

Percentage spore-pollen diagram of Core AKAD 09-15 (Black Sea shelf).

Figure 8.

Percentage diagram of dinocysts, other algae and non-pollen palynomorphs (NPPs) in Core AKAD 09-15 (Black Sea shelf).

Correlations and interpretation of data from spore-pollen analysis

The stratigraphic subdivision of sediments from the western Black Sea area is based on qualitative interpretation of the pollen and spore assemblages and the vertical and spatial distribution of selected indicator taxa. The pollen assemblage zones distinguished are based entirely on the percentage abundances of the predominant and indicator pollen and spores in the assemblages. Pollen spectra delimited for each assemblage zone were obtained from several samples in each sediment core and provide a picture of vegetation changes for the period represented by sediments. According to Birks (1973), the assemblage zones for an individual core have to be considered as local pollen assemblage zones (LPAZ). As these zones are present in two to several sediment cores in adjacent areas, they are further delimited as regional pollen assemblage zones (RPAZ). These RPAZ can be proposed as a regional biostratigraphy and they can be correlated in time and space with concurrent chronostratigraphic scales (Berglund 1983). These LPAZ and RPAZ (Table 3) are tentatively correlated to the Regional Black Sea stratigraphic scale of Shopov (1991), Archaeological Chronology of Todorova (1986) and to the traditional northern European climatostratigraphy of Blytt–Sernander (1876–1908).

Radiocarbon dating and age modelling

Radiocarbon (14C) dating was performed on 15 selected sediment layers (Table 1) spanning all three lithological units and an “Age vs. Depth” model was established for all three cores (Fig. 2). Nine sediment samples from Core Akad 09-10 and Akad 09-15 were selected for 14C AMS dating of bulk organic carbon at the National Ocean Sciences Accelerator Mass Spectrometry (NOSAMS) Facility, WHOI, USA. Six more sediment samples from Core Akad 11-17 were selected for 14C dating at the Gliwice Radiocarbon Laboratory of the Silesian Institute of Technology, Poland. An “Age vs. Depth” model (Fig. 2) was developed for all three cores by 2σ-range calibration of the available radiocarbon dates (Table 1, Fig. 2) with the online calibration programme CALIB 7.1. (Stuiver et al. 2017), using the MARINE13 curve (Reimer et al. 2013) with application of a Reservoir effect (ΔR) of 363 ± 41 yrs (Siani et al. 2000) in order to correlate the obtained results with the available geochronological and archaeological data. The sedimentation rate for this model was calculated by Tilia software (Grimm 2011). This method allows each pollen sample as well as the LPAZ boundaries to be assigned to the relevant calendar age in years BP. All the dates in the pollen diagrams and the text are given as calendar years BP (cal. yrs BP).

Table 1.

Chronology of the studied cores AKAD 11-17, AKAD 09-10 and AKAD 09-15.

Core AKAD Depth (cm) Lab. No Material dated Uncalibrated yrs BP Calibrated yrs BC (2σ range) Calendar yrs BP*
11-17 27 GdA-2598 Bulk 3345±20 909–732 2759
11-17 58.5 GdA-2599 Bulk 7470±25 5722–5537 7584
11-17 134.5 GdA-2601 Bulk 12850±35 12170–11866 13971
11-17 145.5 GdA-2602 Bulk 13660±40 13687–13233 15395
11-17 179.5 GdA-2603 Bulk 15860±45 16587–16175 18346
11-17 228.5 GdA-2604 Bulk 16950±50 17798–17343 19534
09-10 126 OS-79014 Bulk 6920±40 5271–4978 7074
09-10 180 OS-79016 Bulk 19850±170 21500–20592 22985
09-10 212 OS-79017 Bulk 21100±170 23073–22054 24451
09-10 240 OS-79018 Bulk 22400±230 24444–23510 25903
09-15 10 OS-79753 Mollusk 10300±50 9151–8757 10906
09-15 120 OS-79754 Mollusk 10950±45 10181–9582 11877
09-15 170 OS-90808 Mollusk 11500±40 10807–10603 12648
09-15 300 OS-79756 Mollusk 12900±50 12234–11892 14021
09-15 375 OS-74855 Mollusk 13000±60 12562–11954 14144

Results

On the pollen diagrams, six LPAZ from Core AKAD 11-17 (LPAZ AKAD 11-17 1–6) and four LPASZ (LPASZ AKAD 11-17 2 a, b, c, d) (Fig. 3); Six LPAZ from Core AKAD 09-10 (LPAZ AKAD 09-10 1–6) (Fig. 5); Three LPAZ from Core AKAD 09-15 (LPAZ AKAD 09-15 1–3) (Fig. 7) were recognised. They reflect successive changes in vegetation development in the study area. Two local dinoflagellate cyst (dinocyst) assemblage zones (LDAZ 1 and 2) were distinguished (Figs 4, 6, 8). The second dinocyst zone is subdivided into two subzones (LDASZ-2a and LDAZ-2b). All zones are described in detail in Table 2.

Table 2.

Description of the local pollen assemblage zones and subzones (LPA(S)Z) and local dinoflagellate cyst assemblage zones and subzones (LDA(S)Z) from cores AKAD 11-17, AKAD 09-10 and AKAD 09-15.

LPAZ AKAD11-17-1 (228.5–147.5 cm) LPASZ AKAD11-17-2a (147.5–136.5 cm) LPASZ AKAD11-17-2b (136.5–134.5 cm)
19546-15612 cal. yrs BP 15612-14295 cal. yrs BP 14295-14036 cal. yrs BP
Artemisia-Chenopodiaceae-Pinus Pinus-Artemisia-Chenopodiaceae Artemisia-Chenopodiaceae
Dominant non-arboreal pollen (NAP) (up to 72%), mainly Artemisia (around 50%) and Chenopodiaceae (10%). Continuous Poaceae (ca. 2%), Aster-t. (1%) and Brassicaceae (1%). Sporadic Achillea-t., Caryophyllaceae and Scleranthus. Arboreal pollen (AP) dominated by Pinus diploxylon-t. (up to 25%). Regular presence with low values (up to 1%) of Picea, Abies, Juniperus and Ephedra distachya. Gradual decrease of Quercus from 12% up to 3% and of Corylus (3.8-0.5%) at zone top. Continuous low Betula and Ulmus (up to 0.8%). Sporadic Alnus, Carpinus betulus, Fagus and Tilia. Maximum of Pinus diploxylon-t. (up to 42%), coincident with decrease of Artemisia (50 to 35%), constant Chenopodiaceae (to 12.8%). Constant presence of Quercus (ca. 4%), Corylus, Ulmus and Salix (ca. 1%). Sporadic Hippophae. Poaceae, Cichoriaceae, Aster-t., Achillea-t., Brassicaceae and Centaurea jacea-t. are continuously presented. One spectrum: Rapid rise of Artemisia (35-70%) and Chenopodiaceae (up to 16%). Presence of Poaceae (2%) and Achillea-t. (2%). Rapid decrease in Pinus diploxylon-t. (from 42% to 6%). Slight increase in Juniperus and Ephedra distachya. Decline of Quercus (from 4% to 1.5%).
LPASZ AKAD11-17-2c (134.5–125.5 cm) LPASZ AKAD11-17-2d (125.5–99.5 cm) LPASZ AKAD11-17-3 (99.5–63 cm)
14036-13257 cal. yrs BP 13257-11072 cal. yrs BP 11072-8004 cal. yrs BP
Pinus-Artemisia Artemisia-Chenopodiaceae Quercus-Artemisia-Chenopodiaceae
Two spectra: Rapid rise of Pinus diploxylon-t. (from 6% to 25%), decrease of Artemisia (from 70% to 58%) and Chenopodiaceae (from 15% to 10%). Continuous Poaceae (around 2%). Presence of Picea, Quercus, Corylus, Juniperus and Ephedra distachya. Sharp rise of NAP (up to 86%), including highest values of Artemisia (around 67%), max. at the zone bottom. (up to 82%). Gradual increase of Chenopodiaceae (11-16%) and Poaceae (1-3%) at the zone top. Pinus diploxylon-t. after sharp decrease is stabilised below 11%. Constant Quercus (around 2-3%), increased Corylus (up to 1%) at the zone top. Regular low presence of Picea, Juniperus, Ericaceae, Ephedra distachya, Betula, Salix and Hippophae. Presence of numerous heliophytes including Helianthemum, Scleranthus, Caryophyllaceae, Centaurea jacea-t. Gradual rise of Quercus (1-16%), coincidently with decrease of Artemisia (from 64% to 34%). Abrupt decrease in Pinus diploxylon-t. (up to 5%) at the zone top. Trend to rise of Picea, Betula, Corylus and Salix (up to 2%) at the zone top. First significant presence of Carpinus betulus (6%). Presence of Ephedra distachya (0.8%) and Hippophae (0.2%). Typha angustifolia/Sparganium-t., Athyrium and Dryopteris/Thelypteris-t. appear.
LPAZ AKAD11-17-4 (63–44.5 cm) LPAZ AKAD11-17-5 (44.5–27.5 cm) LPAZ AKAD11-17-6 (27.5–3.5 cm)
8004-5483 cal. yrs BP 5483-2837 cal. yrs BP 2837 cal. yrs BP-preindustrial time
Quercus-Corylus-Carpinus betulus-Ulmus-Cerealia-Triticum Quercus-Carpinus betulus-Corylus-Fagus-Carpinus orientalis-Cerealia Quercus-Alnus-Ulmus-Carpinus betulus-Fagus
Dominant Quercus (max. 29%) and Corylus (max. 27%). Consistent increase in Ulmus (0.5-7 %) and Tilia (up to 1.2%). Abundance of Salix and frequent Betula (about 1.1%). Rise of Alnus (up to 3.1%) and Fagus (up to 5.4%) at the zone top. Sporadic Fraxinus excelsior-t., Acer, Pistacia and Hippophae. Constant presence of Hedera. Carpinus orientalis and Fraxinus ornus appear. Decrease in Pinus haploxylon-t. (from 20% to 2.4%). Increase in Picea (up to 1.1%). Abies, Juniperus and Ericaceae are sporadic. Decrease in Artemisia (from 64.4 to 13%) and Poaceae (from 8.5 to 2.6%). Decline of Chenopodiaceae (from max. 24.9% to 2.4%). Appearance of Cerealia-t. (1.3%) and Triticum. Achillea-t., Aster-t., Cichoriaceae, Brassicaceae, Scleranthus, Plantago lanceolata, Polygonum aviculare, Caryophyllaceae and Centaurea jacea-t. present. Polypodiaceaе (around 0.8%), Typha angustifolia/Sparganium-t., Typha latifolia and Cyperaceae are continuously present. Dominant Quercus (about 28.4%). Gradual rise and maximal presence of Carpinus betulus (12.2-18.4%). Decrease in Corylus within the zone (16.8 to 8.1%). Constant presence of Alnus (around 6.1%), Tilia (1.8%), Carpinus orientalis (0.7%), Fagus (6%), Betula (0.1%), Pinus diploxylon-t. (9.5%). Decrease inof Ulmus (6-2.3%) at the zone top. Rise of Ericaceae (up to 0.9%). Sporadic Salix, Acer and Ephedra distachya. Continuous Hedera; Humulus/Cannabis appear. Reduced Artemisia (34% to 8%), Chenopodiaceae (21% to 3.6%) and Poaceae (6.6% to 3.6%). Cerealia-t. and Triticum are continuously presented (around 0.6%). Presence of Aster-t., Achillea-t., Cichoriaceae, Carduus-t., Plantago lanceolata, Polygonum aviculare, Scleranthus, Filipendula and Apiaceae. Polypodiaceae, Cyperaceae, Typha angustifolia/Sparganium-t., Myriophyllum spicatum and Alisma are sporadic. Gradual rise of AP (69 to 84%). Dominant Quercus (around 20%). Rise of Alnus (max. 20%) at the zone top. Continuous Ulmus (around 3.7%) and Fagus (around 5%). Increase in Carpinus betulus (up to 8%) and Corylus (up to 7%). Constant Carpinus orientalis and Tilia (around 1%). Constant Hedera (0.3%). Vitis and Humulus/Cannabis appear. Reduced Artemisia (to 8.4%) at the zone top. Chenopodiaceae and Poaceae are low (2-3%). Continuous Cerealia-t. (1%) and Triticum (0.2%). Plantago lanceolata, Polygonum aviculare, Centaurea cyanus and Rumex are also presented.
LDAZ AKAD11-17-1 (228.5–59 cm) LDASZ AKAD11-17-2a (59–27.5 cm) LDASZ AKAD11-17-2b (27.5–3.5 cm)
19546-7668 cal. yrs BP 7668-2837 cal. yrs. BP 2837 cal. yrs. BP – present
Pyxidinopsis psilata-Spiniferites cruciformis Lingulodinium machaerophorum-Spiniferites belerius-Spiniferites bentorii Lingulodinium machaerophorum-Spiniferites ramosus
Dominant dynoflagelate cysts of Pyxidinopsis psilata (two max. 93.7 and 80.5%), then disappearing. Continuous Spiniferites cruciformis (two max. 20.9% and 23.2% at the zone top), then disappearing. Rise of Pediastrum boryanum var. boryanum (up to 5.1%), Pediastrum simplex var. simplex (around 3.7%), Pediastrum simplex var. sturmii and Pyxidinopsis reticulata are also presented. Other NPPs such as acritarchs Cymatiosphaera globulosa and Hexasteria problematica and fungal spores of Bactrodesmium-type, Sordaria-type, Valsaria valsaroides, Ascospore-type, Biscriate conidium of Alternaria and animal remains of Anthoceras are sporadic. Increase of Spiniferites bentorii (max. 14%) and S. belerius (max. 10.5%) at the zone bottom, followed by rapid increase in Lingulodinium machaerophorum with long processes (max. 107.7%). Substantial presence of S. ramosus, S. mirabilis and S. hyperacanthus. Operculodinium centrocarpum (3.4%), Brigantedinium cariacoense (3.8%), Impagidinium centrocarpum (1.4%) and Echinidinium transparantum (1.5%) appear at the zone bottom and then decrease to 0.1%. Peridinium ponticum, Spiniferites mirabilis, Protoperidinum stellatum and Spiniferites hyperacanthus are sporadic. Increase in acritarchs Cymatiosphaera globulosa (up to 17.2%), then decreasing. First appearance of Micrhystridium cf. ariakense (max. 11.5%), Polykrikos kofoidii (0.6%) and Copepod eggs (0.5%). Sporadic Polykrikos schwartzii, Pleurospora-t.3B, Spirogyra sp. and Selenopemphix quanta. Other NPPs such as Biscriate conidium of Alternaria, Pentaphrasodinium dalei and Glomus-type are rare Decrease of Lingulodinium machaerophorum (48.9% to 1.7%), Spiniferites bentorii (to 0.4%) and S. belerius (to 0.2%). Increase in S. ramosus (up to 4.7%) and Peridium ponticum (to 4.5%). Sporadic presence of Ataxodinium choane, Achomosphaera cf. andalousiense, Fungal spores-type 980, Biscriate conidium of Alternaria, Pediastrum boryanum var. boryanum, Pediastrum simplex var. sturmii. Constant values of acritarchs Cymatiosphaera globulosa (15-5.7%). Presence of Micrhystridium cf. ariakense (max. 5.9%)
LPAZ AKAD09-10-1 (240–160 cm) LPAZ AKAD09-10-2 (160–142 cm) LPAZ AKAD09-10-3 (142–130 cm)
25903-17092 cal. yrs BP 17092-11788 cal. yrs BP 11788-8253 cal. yrs BP
Artemisia-Chenopodiaceae-Pinus Pinus-Artemisia-Chenopodiaceae Quercus-Artemisia-Corylus
Dominant Artemisia (28-49%), constant Chenopodiaceae (6-16%), high Pinus diploxylon-t. (max. 33%). Sporadic Picea and Abies. Maximum of Juniperus (5-10%). Low Quercus (2.7%) and Corylus (1.6-4.3%). Sporadic Ericaceae, Betula, Carpinus betulus, Salix, Ulmus and Alnus. Continuous Polypodiaceae (2%). Many Late-Glacial heliophytes are presented. Dominant Pinus diploxylon-t., decreasing Juniperus (from 7.1% to 1.7%). Appearance of Pinus haploxylon-t. and Ephedra distachya. Decrease in Artemisia (from 42% to 24.2%) and Poaceae (from 7.1% to 3.8%). Constant Chenopodiaceae (around 10%). Rise of Quercus (from 4.3% to 19.8%), constant Corylus (around 2%). Two spectra: Pinus diploxylon-t. strongly reduced to ca. 8% at the zone top. High Quercus (max. 32.5%), increased Corylus (from 5.5% to 15%), Ulmus (up to 4%) and Fagus (to 1.6%). Frequent Carpinus betulus and Betula. Sporadic Humulus/Cannabis and Ephedra distachya. Decreased Artemisia (up to 17.3%), Chenopodiaceae (up to 4.8%) and Poaceae (up to 2.5%).
LPAZ AKAD09-10-4 (130–100 cm) LPAZ AKAD09-10-5 (100–36 cm) LPAZ AKAD09-10-6 (36–10 cm)
8253-ca.5500 cal. yrs BP ca.5500-ca.2800 cal. yrs BP ca.2800 cal. yrs BP – present
Quercus-Corylus-Carpinus betulus-Ulmus-Cerealia-Triticum Quercus-Carpinus betulus-Fagus-Corylus-Carpinus orientalis-Cerealia Quercus-Alnus-Ulmus-Carpinus betulus-Fagus
Dominant Quercus (around 30%), Corylus at max. (28.3%) in the middle of the zone. Increase of Carpinus betulus (from 1.2% to 16.3%), Fagus (from 0.9% to 7.1%) and Alnus (from 1.2% to 4.6%) at the zone top. Continuous Ulmus (around 7%), Fraxinus excelsior-t. (2.8%) ant Tilia (1.9%). Constant Hedera (ca. 0.5%). Appearance of Cornus mas and Acer. Ephedra distachya disappear. Reduced Artemisia (from 24% to 0.6%), Chenopodiaceae (from 3.8% to 0.7%) and Poaceae (from 2.5% to 1.4%). Substantial Cerealia-t. (0.8%) and Triticum (0.8%). Frequent Aster-t. and Achillea-t. (up to 1%). Appearance of Plantago lanceolata and Urtica. Sporadic Filipendula and Scleranthus. Dominant Quercus (from 20% to 28%), Carpinus betulus at max. (17.5%), followed by a decrease to 1.5%. Maximum of Fagus (7.5%), followed by decrease to 2.8%. Decline of Corylus (from 20% to 5.8%) and Ulmus (from 3.7% to 1.6%). First appearance of Carpinus orientalis (1.3%). Continuous presence of Pinus diploxylon-t. (10%). Increase in-Artemisia (from 6.5% to 24.6%). Chenopodiaceae and Poaceae are low (ca. 3-4%). Substantial Cerealia-t. (1%) and Triticum (1%). Presence of Plantago lanceolata (1.2%), Polygonum aviculare, Centaurea jacea and Filipendula. Dominant Quercus (26%); increasing of Alnus (up to 10%) and Ulmus (up to 3%). Constant Carpinus betulus (17%), Fagus (3%), Fraxinus excelsior (1.6%) and Hedera (1.3%). Decrease in Artemisia (from 24.6% to 7.7%), Chenopodiaceae (from 5.5% to 1.1%). Increase in Poaceae (ca. 7.6%). Significant Cerealia-t. (1.6%) and Triticum (0.6%). Appearance of Plantago lanceolata, Scleranthus and Carduus-t.
LDAZ AKAD09-10-1 (240–128 cm) LDASZ AKAD09-10-2a (128–40 cm) LDASZ AKAD09-10-2b (40–10 cm)
25903-7663 cal. yrs BP 7663-ca.2800 cal. yrs. BP ca.2800 cal. yrs. BP – present
Pyxidinopsis psilata-Spiniferites cruciformis Lingulodinium machaerophorum-Spiniferites belerius-Spiniferites bentorii Lingulodinium machaerophorum-Peridinium ponticum
Dominant dinoflagelate cysts of Pyxidinopsis psilata (max. 28.2%) and Spiniferites cruciformis (1.4 to 10.7%). Constant presence of green algal species Pediastrum boryanum var. boryanum (2.2% to 3.4%), Pediastrum simplex var. sturmii (around 3.7%), Botryococcus sporadic. High Glomus-t.207 (1.5% to 6.3%). Multiplicasphaeridium-t., Pleurospora-t.3B, Ascospores-t.20, Achomosphaera cf. andalousiense and Cymatiosphaera globulosa are sporadic. First appearance of Spiniferites ramosus at the zone top. Lingulodinium machaerophorum appear sporadically at the zone top. Rapidly increasing dominant Lingulodinium machaerophorum with long processes (two max. 446% and 166%). L. machaerophorum f. clavate at max 28.7% at the zone bottom, followed by sharp decrease. Spiniferites belerius (max. 26.2%) and S. bentorii (max. 12.9%) form peaks at the zone bottom, then decreasing up to around 6%. High Brigantedinium cariacoense (around 10%). Substantial S. mirabilis, S. membranaceus and S. hyperacanthus. Constant Polykrikos kofoidii (0.7-5%), S. ramosus (2.2%), Operculodinium centrocarpum (0.9-7.2%), Echinidinium transparantum (0.6-5.4%) and Impagidinium aculeatum (0.5-1.8%). Sporadic Tectatodinium pellitum, Polykrikos schwartzii, Protoperidinum stellatum and Bitectatodinium tepikiense. Rising of acritarch Cymatiosphaera globulosa (18.3%), Micrhystridium cf. ariakense (max. 20%). Pleurospora-t.3B and Copepod eggs are also presented. Dominant Peridinium ponticum (10%). Decrease in Lingulodinium machaerophorum (3.7%), Spiniferites bentorii and S. belerius. Polykrikos schwartzii increase (up to 2%). Protoperidinum nudum and Operculodinium centrocarpum are sporadic. Decrease in acritarchs Cymatiosphaera globulosa (1.8%). Pleurospora sp.-t.3B, Fungal spores-t.200, Sordaria and Copepod eggs are also presented.
LPAZ AKAD09-15-1 (377–320 cm) LPAZ AKAD09-15-2 (320–200 cm) LPAZ AKAD09-15-3 (200–10 cm)
14147-14054 cal. yrs BP 14054-12965 cal. yrs BP 12965-10906 cal. yrs BP
Artemisia-Chenopodiaceae-Pinus Pinus-Artemisia-Chenopodiaceae Artemisia-Chenopodiaceae
Dominant Artemisia (max. 44%), constant Chenopodiaceae (11%), high Pinus diploxylon-t. (39-42%). Continuously low Quercus (2.7%), Corylus (1.5%), Salix (1%), Betula (1%), Ulmus (0.6%), Juniperus (2%), Hippophae (0.5%). Presence of Ephedra distachya (up to 0.8%). Continuous Poaceae (around 7%). Numerous Late-Glacial heliophytes are presented. Dominant Pinus diploxylon-t. (max. 61%) decreasing Artemisia (from 31.5% to 12%); Chenopodiaceae peak in the middle of the zone (23.6%). Increase in Poaceae to 10%. Many Late-Glacial heliophytes are presented. Constant Quercus (4%), Salix (2%), Betula (2%), Corylus (1%). Dominant Artemisia (max. 67%) and Chenopodiaceae (24.6%); constant Poaceae (6-8%). Late-Glacial heliophytes are abundant. Typha angustifolia/Sparganium and Cyperaceae appear. Increase of Polypodiaceae (1-2%).
LDAZ AKAD09-15-1 (377–10 cm)
14147-10906 cal. yrs BP
Pyxidinopsis psilata-Spiniferites cruciformis
Dominant dynoflagelate cysts of Pyxidinopsis psilata (10%) and Spiniferites cruciformis (4%). Constant presence of green algal species Pediastrum boryanum var. boryanum (2%). Sporadic acritarchs of Pseudoschizaea circula. Constant presence of Glomus-type 207 (max. 11% at the zone bottom).

Discussion

Pollen assemblages

Regional Pollen Assemblage Zone IV (ArtemisiaChenopodiaceaePinus)

25903–15612 cal. yrs BP

RPAZ IV has a Late Pleniglacial-Oldest Dryas/Upper Neoeuxinian age (25903–15612 cal. yrs BP). It is represented in cores AKAD 09-10 (LPAZ AKAD 09-10 1 and 2) (Fig. 5, Tables 2, 3) and AKAD 11-17 (LPAZ AKAD 11-17 1) (Fig. 3, Tables 2, 3). Sediments with Late Pleniglacial and Oldest Dryas age cannot be separated by pollen analysis in cores studied, because of the similarity in pollen spectra.

Table 3.

Correlation between local and regional pollen assemblage zones and subzones (modified after Filipova-Marinova (2006a)).

cal. kyrs. ВР (Fig. 2) Northerneuropean climatostratigraphy (Blytt 1876Sernander 1908) Regional stages and substages (Shopov 1991) Archaeological Chronology (Todorova 1986) Regional PAZ subzones Pollen assemblages (Filipova-Marinova 2006) 2345 09-15 GGC18 09-10 544 11-17
Local PAZ Pollen assemblages (Filipova-Marinova 2003b) Local PAZ Pollen assemblages (Table 2) Local PAZ Pollen assemblages (Filipova-Marinova et al. 2013) Local PAZ Pollen assemblages (Table 2) Local PAZ Pollen assemblages (Filipova et al. 1989) Local PAZ Pollen assemblages (Table 2)
2.8 H O L O C E N E Subatlantic B L A C K S E A New Black Sea Iron Epoch IX Q-U-Al-Cb-Sa-F 8 Q-U-Al-F 6 Q-Al-Cb-U-F 6 Q-Al-U-Cb-F 6 Q-Al-Cb-U-F 6 Q-Al-U-Cb-F
5.5 Subboreal Old Black Sea Early Bronze Age VIII Q-Cb-Co-F-Cor-Ce 7 Q-Cb-Co-F 5 Q-Cb-F 5 Q-Cb-F-Co-Cor-Ce 5 Q-Cb-Co-F 5 Q-Cb-Co-F
8.2 Atlantic Transi-tional Period VII Q-Co-Cb-U-Tr-Ce 6 Q-Co-U-Cb-Ti 4 Q-Co-Cb-U-F- 4 Q-Co-Cb-U-Ce-Tr 4 Q-Co-U-F-Cb 4 Q-Co-Cb-U-Ce-Tr
Late Eneolithic
Neolithic
11.7 Boreal VI b Q-U-Co-Art stratigraphic hiatus 3 Q-U-Art 3 Q-Art-Co 3 Q-Co-U-Cb-Art 3 Q-Art-Ch
Preboreal a Q-P-Art 2 Q-P-Art
13 P L E I S T O C E N E Late Würm Late Glacial Younger Dryas Upper-New-euxinian V d Art-Ch 5 Art-Ch 3 Art-Ch 1 Art-Ch-P 2 P-Art-Ch 2 Art-Ch 2d Art-Ch
14 Allerød c P-Art-Ch 4 P-Art-Ch 2 P-Art-Ch 1 P-Art-Ch 2c P-Art
14.3 Older Dryas b Art-Ch-Po 3 Art-Ch-Po 1 Art-Ch-P 2b Art-Ch
15.6 Bølling a P-Art 2 P-Art 2a P-Art-Ch
25.9 Oldest Dryas IV Art-Ch-P 1 Art-Ch-P 1 Art-Ch-P
Late Pleniglacial 1 Art-Ch-P

Paleovegetation reconstruction, based on typical high values of non-arboreal taxa and the presence of cold-resistant and heliophilous taxa such as Artemisia and Chenopodiaceae, suggest spreading of cold and dry steppes. Different taxa of Poaceae were also important elements in the steppe communities along the coast together with other heliophilous taxa from Asteraceae, Cichoriaceae, Apiaceae, Brassicaceae, Caryophyllaceae and Helianthemum. According to Prentice et al. (1993), the growing season soil moisture deficit and low winter temperature would have maintained an open vegetation. However, climate aridity along the Bulgarian Black Sea coast was not so extreme because the presence of the desert shrub Ephedra distachya, that is an indicator of dry continental climate and an important constituent of the glacial flora, is constant but low. Patches of eurythermic conifers such as Pinus grow along the coast in valleys where favourable microclimatic conditions may have prevailed. The high percentage values of Pinus diploxylon-type pollen could also be an effect of long distance transport which became significantly more pronounced in an environment with sparse vegetation (Pardoe et al. 2010). Warming/humidity increase during the melting pulses (19.2 to 14.5 cal. ka BP) after the Last Glacial Maximum is expressed by the open coniferous forests development. The constant maximum presence of Glomus spores suggests significant erosion processes.

Temperate deciduous arboreal taxa such as Quercus, Corylus, Ulmus, Betula and Alnus show constant presence in the pollen diagrams. The presence of single pollen grains of Tilia, Carpinus betulus, Abies, Ulmus, Fraxinus excelsior and fern spores of Polypodiaceae suggests an increase of temperature that is seen in the paleotemperature record of Greenland Ice sheet-2 (GISP-2) (Blunier and Brook 2001). Probably, warming/humidity increase during these melting pulses is expressed by the open coniferous forest development. Isolated patches of these taxa must have survived in micro-environmentally favourable pockets. These data suggest that these small pockets in the Black Sea coastal mountains play an important refugial role in the survival of temperate arboreal taxa along south-eastern Europe during the Last Full Glacial (Beug 1975; Bottema 1980; Bennett et al. 1990; Tzedakis 1993; Willis 1994). Our results corroborate with data of core 22-GC3 from the southern Black Sea region (Shumilovskikh et al. 2012), Ioanina (north-western Greece) (Tzedakis 1993) and Thenagi Phillipon (Bottema 1980). Similar pollen spectra were found at ca. 27990 ± 300 14C yrs BP by Frenzel (1964) at Stilfied in Austria and were related to the Paudorf Interstadial. This interstadial is also established, although less clearly and dated at 27295 ± 1120 14C yrs BP for the Black Sea by Komarov et al. (1979). Mudie et al. (2002, 2007) recognised a Pleniglacial pollen assemblage zone with relatively low AP values suggesting a steppe-forest vegetation and moderately high precipitation in winter for Marmara Sea core MAR 95-04, with age of 29540 ± 1540 to 21950 ± 310 14C yrs BP and in the south-western Black Sea Core MAR 98-04 which has a radiocarbon age of 33550 ± 330 14C yrs BP (Aksu et al. 2002).

The Paleoclimate model of Peyron et al. (1998) for the Last Glacial maximum (LGM) shows that the climate in southern Europe during the Pleniglacial was characterised by extremely low temperatures and humidity throughout the year. Reconstructed temperatures for the area north of the Mediterranean Sea were lower than today: -15 ± 5 °C for the coldest month and -10 ± 5 °C for the annual mean temperature. In Greece, the available moisture was 20% lower with a precipitation anomaly of ca. 600 ± 200 mm. Kutzbach and Webb (1993) also suggested the Full Glacial conditions in eastern Europe with extremely cold and arid climate and predicted January temperatures as low as -20 °C.

These marine deposits are found below the 30 m isobath in almost all investigated cores of the western Black Sea shelf. In the peripheral (outer) shelf zone, they form clearly defined depositional bodies of coastal or barrier type sediments at the depth of 100 to 120 m (Khrischev and Shopov 1978). Chepalyga (2002) points out that the eustatic regression of the World oceans had led to loss of the two-way connection with the Mediterranean Sea and a subsequent evaporative drop in the level of the Black Sea to a strand-line of -90 m or lower.

Regional Pollen Assemblage Zone V

15612–11788 cal. yrs BP

This zone is distinguished in all cores studied and could be correlated with the Late Glacial/Neoeuxinian age. All stadials and interstadials were determined palynologically and four subzones were presented.

Regional Pollen Assemblage Subzone Va (PinusArtemisia)

15612–14295 cal. yrs BP

Subzone Va corresponds to the Bølling Interstadial of the European Late Glacial and can be referred to as the Upper Neoeuxinian. It is clearly separated only in core AKAD 11-17 (LPASZ AKAD 11-17 2a) (Fig. 3, Tables 2, 3), but cannot be separated in core AKAD 09-10 (LPAZ AKAD 09-10 2) (Fig. 5, Tables 2, 3) because the pollen analysis in this core section is not of high resolution. The basic characteristic of the vegetation succession during this period of moderate warming is the restriction of areas occupied by steppes dominated by Artemisia and the existence of open forests and scattered patches with Pinus and Quercus along the coast. The contribution of the thermophilous broad-leaf arboreal taxa and the decrease of spread of the desert shrub Ephedra distachya suggest climate amelioration. Although some of the Pinus diploxylon-type pollen could be of long distance origin, these high percentages (more than 50%) suggest a local and regional expansion of Pinus at the expense of shrubs and herbs, as is clearly recorded. This rapid interstadial expansion testifies that Pinus survived along the coast during the Late Glacial period. Detached stands of some temperate arboreal taxa, such as Quercus, Ulmus, Corylus and Salix, may also have been preserved in suitable localities. Increase of Pinus diploxylon-type and some thermophilous pollen taxa was also reported for the Black Sea sediments by other authors (Shimkus et. al. 1977; Traverse 1978b; Komarov et al. 1979) and for the Marmara Sea (Mudie et al. 2002).

Regional Pollen Assemblage Subzone Vb (ArtemisiaChenopodiaceaePoaceae)

14295 – 14036 cal. yrs BP

Subzone Vb suggests stadial environmental conditions. The available Age vs. Depth model (Fig. 2) allows this subzone to be correlated to the Older Dryas Stadial of the Würm Late Glacial and to the Upper Neoeuxinian. It is represented in cores AKAD 09-15 (LPAZ AKAD09-151) (Fig. 7, Tables 2, 3) and AKAD 11-17 (LPASZ AKAD11-17 2b) (Fig. 3, Tables 2, 3). The Older Dryas was very short-lasting and cool. It is always hard to recognise it in sediments on the basis of paleoflora. This is the reason many authors accept Bølling and Allerød as an interstadial complex in the terrestrial pollen records from central and eastern Europe (Lang 2003). The precise distinction and correlation between cores is difficult and obscured due to the differences in both the processes of sediment formation and the succession of local vegetation between particular parts of the coast. Pollen spectra are dominated by Artemisia and Chenopodiaceae, suggesting the wide spread of xerophytic and halophytic herb communities along the coast. The percentage presence of Pinus diploxylon-type is slightly higher in LPAZ AKAD 09-15 1 from the Black Sea shelf for the time span 14174 to 14059 cal. yrs BP compared to other cores studied from the deep-water zone. High values of non-arboreal pollen suggest a cold and dry climate which is also confirmed by the presence of Juniperus and Ephedra distachya pollen.

Regional Pollen Assemblage Subzone Vc (PinusArtemisiaChenopodiaceae)

14036 – 12965 cal. yrs BP

Subzone Vc is considered as an interstadial sequence analogous with the Allerød Interstadial of the Late Glacial and correlated to the Upper Neoeuxinian. It is represented in cores AKAD 09-15 (LPAZ AKAD 09-15 2) (Fig. 7, Tables 2, 3) and AKAD 11-17 (LPASZ AKAD11-17 2c) (Fig. 3, Tables 2, 3). The steppe vegetation is still rich in plant taxa, but more restricted in area, particularly the vegetation of dry habitats. The pollen record shows AP increase, particularly Pinus diploxylon-type pollen and presence of temperate arboreal taxa, such as Quercus, Corylus, Ulmus, Salix and Betula. The maximum of Pinus diploxylon-type pollen occurring in Core AKAD 09-15 might signal climate changes (increase of temperature and humidity). The increase in temperate arboreal pollen indicates the limited migration of some trees from the south-eastern European refugia. van der Hammen et al. (1971) and Beug (1982) also proposed the existence of refugia in southern Europe at elevations from 500 to 800 m a.s.l.

Regional Pollen Assemblage Subzone Vd (ArtemisiaChenopodiaceae)

12965 – 11788 cal. yr BP

Subzone Vd is associated with the last most significant rapid climate deterioration of the last Late Glacial Stage, i.e. the Younger Dryas Stadial and has an Upper Neoeuxinian age. This cold period is of global importance and is recognised everywhere in Europe as an episode of pronounced cooling (Berglund et al. 1994). The subzone Vd is represented in cores AKAD 09-15 (LPAZ AKAD 09-15 3) (Fig. 7, Tables 2, 3) and AKAD 11-17 (LPASZ AKAD 11-17 2d) (Fig. 3, Tables 2, 3). The reversal from interstadial towards stadial conditions occurred after 12965 cal. yrs BP, corresponding to the termination of the Allerød phase in western Europe and to the GISP-2 18O climate record (Stuiver et al. 1995). This phase is also clearly recognised in the Marmara Sea (Mudie et al. 2002, 2007) and in Lake Van (Wick et al. 2003).

This succession reflects the expansion of xerophytic herb (steppe) vegetation. Palynological data show that, in addition to the predominant light-demanding xerophytic and halophytic taxa such as Artemisia and Chenopodiaceae, many other taxa, such as Poaceae, Aster-type, Achillea-type, Centaurea, Thalictrum, Apiaceae and Caryophyllaceae have also participated in these steppe communities. The extremely high percentages of Artemisia suggest that the Younger Dryas climate of the Bulgarian Black Sea coast is analogous to that of the Last Glacial maximum. Chenopodiaceae is always subdominant to Artemisia during the Pleniglacial and the Late Glacial Stadials (Atanassova 2005; Mudie et al. 2007; Shumilovskikh et al. 2012; Filipova-Marinova et al. 2013). This confirms that Younger Dryas is the coldest and driest period along the northern Bulgarian Black Sea coast during the whole Late Glacial. Strong evidence of this climatic deterioration can be seen in the increase and continuous presence of Juniperus and the indicator of the cold and dry climate Ephedra distachya. High percentage values of Artemisia and Chenopodiaceae may have also been a result of the wide extent of halophytic species such as Salsola ruthenica, Suaeda maritima, Salicornia europaea and Artemisia maritima growing on the beach area and on the part of the modern shelf after the withdrawal of the sea to -100 m during the Neoeuxinian regression (Chepalyga 2002). There is a steep decline in AP, reflected by a decrease of pollen of Pinus diploxylon-type from 13257 cal. yrs BP onwards and restricted occurrence of almost all deciduous arboreal taxa. Stands of Pinus, Quercus, Corylus, Carpinus betulus, Ulmus, Tilia and Betula were sparsely distributed in favourable moisture localities in the Eastern Balkan Range (Stara Planina Mts.). This supports the hypotheses of Beug (1975) that refugia of deciduous trees would have been located at mid-altitude sites where the precipitation would have been higher than on the plains during this arid glacial period. This information is of great importance in tracing the main migration routes of these deciduous trees along the Bulgarian Black Sea coast after the last glaciation. Probably, stands of Pinus occupied the higher parts of the coastal plateaux, where conditions were more favorable for the growth of trees because of higher atmospheric humidity.

Multiple factors may have been important in determining the vegetation changes in this region, including climate oscillations (Cordova et al. 2009). Such a vegetation type can be attributed to the reduced pollen production under the influence of cold and dry glacial conditions (Willis 1994). Independent evidence provided by isotope analyses makes the distinction that timing and duration of the Late Glacial and Early Holocene aridity along the Black Sea and eastern Mediterranean is strongly linked with the summer insolation (Wright et al. 2003; Tzedakis 2007). Palynological records from the Black Sea region also confirm the results from climate modelling experiments of Kutzbach and Webb (1993) that summer insolation maxima override the effects of the North Atlantic circulation in the continental lowlands of the Balkans. Similar pollen assemblages appear also in other marine cores from the western Black Sea shelf and deep-water zone and show that this climate oscillation has an ambiguous signal in palynological records and is clearly apparent along the whole Bulgarian Black Sea coast (Atanassova 2005; Shumilovskikh et al. 2012; Filipova-Marinova et al. 2013). The Younger Dryas event in the western Black Sea area is also dated by the Artemisia maximum at ~ 10660 cal. yrs BP (Mudie at al. 2007). A large increase in non-arboreal pollen (NAP) with peaks in Artemisia and Ephedra and increased Chenopodiaceae is found in Core MAR 98-12 of the nearby Marmara Sea during the Younger Dryas (Mudie et al. 2002; Caner and Algan 2002). According to Niklewski and van Zeist (1970) and Connor et al. (2013), aridity, rather than temperature during the glacials, seems to have been the key factor limiting lowland forest development during this steppe phase.

Regional Pollen Assemblage Zone VI (QuercusPinusArtemisia)

11788 – 8004 cal. yrs BP

This zone could be correlated with the Preboreal-Boreal chronozone/Old Black Sea Substage and is represented in cores AKAD 09-10 (LPAZ AKAD 09-10 3) (Fig. 5, Tables 2, 3) and AKAD 11-17 (LPAZ AKAD 11-17 3) (Fig. 3, Tables 2, 3). The studied sediments with Preboreal age cannot be separated from that of the Boreal age in Cores AKAD 09-10 and AKAD 11-17, because the pollen analysis in this section is not of high resolution. This zone is marked by the abrupt change from ArtemisiaChenopodiaceae to a QuercusPinusArtemisia assemblage. The increased values of total AP and deciduous arboreal taxa suggest an establishment of pioneer formations mainly of Quercus and Pinus, most probably in a mosaic structure. This change would be in line with the notion of Giesecke et al. (2011) that the beginning of the Holocene is a period of pronounced and frequent change in the vegetation composition. The pollen records from the Black Sea continental slope and deep-water zone suggest a forest-steppe vegetation succession. The vegetation was relatively open, based on the diversity of xerophytic and mesophytic herbs in the pollen records. Abundance of non-arboreal taxa, such as Artemisia, Chenopodiaceae, Poaceae and Asteraceae was still high and may indicate that humidity was a major limiting factor in the widespread distribution of forests along the coast or that there were still large areas of disturbed open soil. The first step in the afforestation along the western Black Sea coast began after 11788 cal. yrs BP with the rapid increase of some arboreal taxa.

The most characteristic feature for the Early Holocene vegetation palaeosuccession is the early appearance of Quercus as a pioneer element in open pine forests, while, in central and northern Europe, the light-demanding species Corylus avellana started to spread in open forests dominated by Pinus and Betula where interspecies competition was probably of little importance (Birks and Line 1993; Tallantire 2002; Finsinger et al. 2006). The example of Corylus avellana shows that the same species can show different behaviour with respect to its regional expansion. Betula does not play a prominent role along the Bulgarian Black Sea coast in contrast to high mountain areas (Bozilova and Tonkov 2000; Stefanova et al. 2006). The pollen diagrams from the southern Black Sea region (north-western Turkey) suggest greater occurrence of Betula during the Early Holocene than in the modern vegetation cover (Bottema 1990).

The characteristic expansion of Quercus is due to the increase in temperatures and humidity. Probably, different oak species, such as Q. cerris, Q. frainetto, Q. pubescens and Q. polycarpa, took part in the composition of these forests. In addition to Quercus, several temperate taxa, such as Ulmus and Tilia, were also present in these forests. The presence of Ulmus supports the assumption of Stojanov (1950) that, in the past, Ulmus forests were the most important components of the vegetation of the lowlands. Carpinus betulus is still an insignificant component of the oak forests, as well as Corylus, Tilia and Fagus.

The first increase of deciduous arboreal pollen, mainly of Quercus is dated at about 9630 ± 520 14C yrs BP at Core-544 of the deep-water zone (Filipova et al. 1989) and about 9945 ± 160 14C yrs BP at Core-149 of the nearby estuary of the Veleka River (Filipova-Marinova 2003a). Mudie et al. (2002) also reportеd that Quercus, Carpinus betulus and Ulmus were present in the Early Holocene and began to expand after 8500 14C yrs BP. Similar changes in vegetation composition along the coast during that time have also been documented by Atanassova (2005) and Shumilovskikh et al. (2012). Such vegetation change is classically observed for the beginning of the Holocene at the Mediterranean area (Tzedakis 2007). In Romanian Carpathians, the spread of oak forests occurred after that of Ulmus at about 10800 cal. yrs. BP (Feurdean et al. 2010), while, at higher elevations of Bulgarian mountains, it started together with the expansion of Betula forests after 11200 cal. yrs BP (Tonkov et al. 2002). Early pollen increase in Quercus, Ulmus, Tilia and Acer in marine records indicates rapid expansion of these temperate taxa during the Early Holocene, likely from local sources and suggests that they probably had Glacial refugia in the Strandzha Mountains and the Stara Planina Mountains (the Balkan Range). According to Lang (1985), it is not certain if Quercus, Ulmus, Tilia, Fraxinus and Acer occurred together in the mixed oak forests of Europe, because the distribution of these taxa depends on edaphic conditions and topography. Most probably, in the Bulgarian coastal region, Quercus, Ulmus and Corylus were distributed on the richest soils of the southern slopes and the hills of the Strandzha Mountains and the Stara Planina Mountains (the Balkan Range), while Fagus and Carpinus betulus were spread on northern slopes and along the humid ravines. The early expansion of temperate arboreal taxa and the presence of pollen grains of Fagus, Rhododendron, Ericaceae and Juglans support the notion of van der Hammen et al. (1971), Huntley and Birks (1983) and Bennett et al. (1990) that there were refugia of these taxa in the Strandzha Mountains during the glaciations.

The Sofular Cave (Zonguldak Province, north-western Turkey) record also suggests a fast re-vegetation with trees and shrubs at the onset of the Early Holocene (Fleitmann et al. 2009). These results are in accordance with the presumption of Leroy and Arpe (2007), that parts of the Black Sea mountains were Glacial refugia for temperate trees, which facilitated their rapid re-advance at the onset of the Holocene (11600 cal. yrs BP).

Brewer et al. (2002) consider that the climate of the Early Holocene acted as the strongest controlling factor on the spread of the oak. According to Davis et al. (2003), the Early Holocene warming and later equilibrium has been mainly modulated by increased winter temperatures. Shimkus et al. (1977) and Komarov et al. (1979) also suggest climate warming and increase of arboreal taxa along the Black Sea coast from 10737 ± 315 14C yrs BP. Around the Marmara Sea, the amelioration of climate occurred around 10200 14C yrs BP (Mudie et al. 2002).

Warming at the onset of the Holocene also allowed a rapid spread of oak along the Atlantic coast of Europe (Brewer et al. 2002). According to Berglund et al. (1984), such vegetation change is probably a biotic reaction of the climate improvement before about 10000 years including the rise of temperature. The delay of spreading of other arboreal species, such as Carpinus betulus, Corylus, Fagus and Alnus, along the coast was probably due to the low humidity of the Early Holocene climate.

Pinus diploxylon is noted as the major contributor to the Early Holocene pollen assemblages (Mudie et al. 2002). However, along the Bulgarian Black Sea coast, percentage values of Pinus diploxylon-type declined, while the deciduous arboreal taxa expanded their presence in the pollen spectra. Taking into account the over-representation of Pinus pollen in modern surface samples from the western and central Balkan Range (Stara Planina Mts) (Filipovitch et al. 1997) and from the Black Sea coast (Pardoe et al. 2010), the high values of Pinus pollen could be partly due to the long distance transport. The formation of isolated stands of Pinuson higher localities on the Bulgarian Black Sea coast could not be excluded.

A rapid and very short-term steep decline of arboreals is established in marine core AKAD 11-17 from the deep-water zone at 8500 to 8300 cal. yrs BP (70–65 cm) and in core AKAD 09-10 from the continental slope at 8253 cal. yrs BP at the transition of Boreal and Atlantic chronozones (130 cm) (Figs 3, 5). This event could be caused by the changes in climate parameters or short climate extrusion that would trigger synchronous shifts in vegetation composition. For the first time, evidence for such an abrupt cooling event was identified in the pollen diagram from Core GGC-18 from the south-western Black Sea continental slope (Filipova-Marinova et al. 2013). This short-term arid phase in Black Sea coastal palaeoclimate is confirmed by a decrease in arboreal pollen, particularly of Quercus and Ulmus simultaneously with the increase of Artemisia, presence of the indicator of cold climate Ephedra distachya and low sedimentation rate between 8500 and 8253 cal. yrs BP. This event should be considered as a vegetation response to the “8200 yrs BP cold event” (Magny et al. 2003; Alley and Ágústsdóttir 2005; Bahr et al. 2005; Le Grande et al. 2006) as other pollen diagrams from Germany and Finland also show such synchronous pronounced short decline of AP between 8500 and 8300 cal. yrs BP (Giesecke et al. 2011). According to the pollen-based temperature reconstructions, a significant drop in average winter temperatures in the order of 4 °C has been estimated during this rapid climate event, while summer temperatures appear to have been comparable to those currently prevailing in the north-western Carpathians (Feurdean et al. 2008). In Thenagi Philippon (northern Greece) during the “8200 yrs BP event”, a significant reduction in tree pollen is observed, that is representative of a decline in winter temperatures of more than 4 °C caused by perturbation of North Atlantic circulation (Pross et al. 2009). Tinner and Lotter (2001) consider this event as a very rapid synchronous response of southern European vegetation to suborbital climate change. In the Balkan Region, this rapid climate change was also identified in the Holocene sequences from the Carpathians (north-western Romania) (Feurdean et al. 2010), Thenagi Philippon (north Greece) (Müller et al. 2011) and Lake Prespa (Macedonia) (Panagiotopoulos et al. 2012). In the Bulgarian high altitudes in the Rila Mts., this event is manifested by the decline of arboreal pollen accumulation rate to ca. 500 grains cm-2/yr-1 at 8230 cal. yrs BP, particularly by the decline in Pinus diploxylon-type and less of Betula (Tonkov et al. 2016). The decrease in temperate forests most probably reflects a decrease in humidity (Combourieu Nebout et al. 2009; Dormoy et al. 2009). A shift in moisture regime led to a shift in the vegetation composition (Tinner and Lotter 2006). According to Roberts et al. (2011), deciduous oak species require high soil moisture during the summer, probably a decrease in soil moisture level at summer contributing to the Quercus declining trend observed in the cores studied.

Regional Pollen Assemblage Zone VII

(QuercusCorylusCarpinus betulusUlmusTriticumCerealia)

8004 – 5483 cal. yrs BP

This RPAZ can be correlated to the Atlantic chronozone of the Middle Holocene that corresponds to the Old Black Sea Substage. It is represented in cores AKAD 09-10 (LPAZ AKAD 09-10 4) (Fig. 5, Tables 2, 3) and AKAD 11-17 (LPAZ AKAD 11-17 4) (Fig. 3, Tables 2, 3). The temperate deciduous arboreal pollen taxa reach their maximum, indicating rapid spread of stable mixed oak forests. Their maximal distribution is recorded at 6417 cal. yrs BP in LPAZ AKAD 11-17 4 and at 5745 cal. yrs BP in LPAZ GGC18 4. The high values of AP suggest that forests become denser. Quercus appears to be the major arboreal taxon in the mixed oak forests with abundant Corylus, Carpinus betulus, Ulmus, Fraxinus excelsior, Tilia and Acer. The optimal climate conditions as well as moisture balance of the Black Sea Region were determining factors that stimulated the extensive spreading of these forests. This expansion of mixed oak forests was also attributed to the rise of the Black Sea level during this period resulting in an increase in atmospheric moisture content. According to Tonkov and Bozilova (1995), high humidity on the Balkan Peninsula was reached after 8000 14C yrs BP and the climate was determined by the transport of air masses from the Atlantic Ocean. The presence of the indicator taxa, such as Hedera and Humulus/Cannabis, confirms an increase in humidity and rise in mean annual temperatures.

Corylus expanded in the local stand mainly at the expense of Ulmus and became widespread from 7584 to 5483 cal. yrs BP. This taxon has high pollen productivity in open areas. Pollen data suggest the great extent of monodominant communities of Corylus in open areas, but probably also as an undergrowth of the oak forests. The maximum percentage values of Corylus could be associated with a short-term fluctuation of climate parameters, but also with a clearance of mixed oak forests for enlargement of cultivated areas along the coast as is seen by the synchronous maximum values of Cerealia-type pollen (Filipova-Marinova 2006a). Interspecies competition could not be excluded, because the rapid reduction of Corylus values corresponds to the beginning of the continuous curves of Carpinus betulus and Fagus in all pollen diagrams from the Bulgarian Black Sea coast from the end of Atlantic around 5500 cal. yrs BP. The increase of Carpinus betulus started after 6417 cal. yrs BP in LPAZ AKAD 11-17 4 (Fig. 3). The first regular occurrence of Fagus suggests the presence of this taxon at about 9000 cal. yrs BP. The expansion of Fagus started as early as 8283 cal. yrs BP (8355 ± 75 14C yrs BP) (Filipova-Marinova 2003b) at the southern part of the Bulgarian Black Sea coast, while at the northern Bulgarian Black Sea coast, the areas occupied by Fagus expanded after 6500 cal. yrs BP (LPAZ AKAD11-17 4) reflecting favourable conditions for its spreading and it started to be better represented in suitable localities at lower altitudes. The maximum of Fagus occurred around 3260 cal. yrs BP (3070 ± 100 yrs BP) (Filipova 1985). Nowadays, communities of this species could be found in the Eastern Balkan Range at altitudes between 170 and 450 m a.s.l. (Bondev 1991). Another broad-leaf tree pollen type in the diagrams is Acer. A regular occurrence of Acer pollen from 7750 cal. yrs BP is indicative of a local presence of this taxon.

The presence of Alnus, together with several occasional pollen grains of Hedera, confirms the increase in humidity and temperature along the coast. Submediterranean elements such as Carpinus orientalis and Fraxinus ornus also occurred near the coastline. Carpinus orientalis appeared and probably occupied some areas after the degradation of mixed oak forests due to a human impact that influenced the natural vegetation. The first occurrence of Juglans is registered at 7584 cal. yrs BP in LPAZ AKAD 11-17 4. The earliest appearance of several occasional pollen grains of Juglans along the southern Bulgarian Black Sea coast during the Holocene is registered for Preboreal, ca. 10000 cal. yrs BP (Filipova-Marinova 2003a), confirming the possible relic origin of this taxon in the Balkan Peninsula (Bottema 1980).

The first appearance of pollen of anthropophytes, such as Cerealia-type, Triticum, Plantago lanceolata and Polygonum aviculare, coincides with the decline of Corylus and Ulmus, marking human impact during the Late Eneolithic period, 6790–6320 cal. yrs BP (Bozilova and Beug 1992, 1994; Filipova-Marinova and Bozilova 2002, 2003). Human impact was significant along the western Black Sea coast as reflected in pollen diagrams from the Bulgarian Black Sea coastal lakes (Todorova 1986; Bozilova and Filipova 1991; Bozilova and Beug 1994). Terrestrial pollen records from the eastern Black Sea coast reveal anthropogenic impact starting at 4000 14C yrs BP (Connor et al. 2007).

Regional Pollen Assemblage Zone VIII

(QuercusCarpinus betulusCorylusFagusCarpinus orientalisCerealia)

5483 – 2837 cal yrs BP

The characteristic vegetation succession and the available Age vs. Depth model allow correlation of this RPAZ to the Subboreal chronozone of the Middle Holocene that also corresponds to the Old Black Sea Substage. It is represented in cores AKAD 09-10 (LPAZ AKAD 09-10 5) (Fig. 5, Tables 2, 3) and AKAD 11-17 (LPAZ AKAD 11-17 5) (Fig. 3, Tables 2, 3). The pollen record shows that mixed oak forests still dominate, but a change in the forest composition is registered after 5483 to 2857 cal. yrs BP. The most prominent feature for this zone is the increase of Carpinus betulus. The maximum of Carpinus betulus pollen is also dated at 5680 ± 65 14C yrs BP in the Arkutino Lake (Bozilova and Beug 1992) and at 5650 ± 100 14C yrs BP in the Shabla-Ezeretz Lake (Filipova 1985). According to Filipovitch et al. (1998), this species formed a separate belt on higher areas in the Balkan Range at about the same time. The hornbeam belt was also common on eastern Carpathians during the Subboreal and dated to around 4210 ± 35 14C yrs BP (Tantau et al. 2011). The presence of Corylus, Ulmus, Fraxinus excelsior, Tilia and Acer suggest that these taxa also appeared in the mixed oak forests. The decline of Corylus started at the same time as Fagus and Carpinus betulus became established around 5500 cal. yrs BP. The expansion of Fagus during the Subboreal must have been triggered by an external factor. This factor may have been climatic change or human influence which started to increase at least in lowland areas. Forests of Fagus orientalis expanded in the moist ravines of the Strandzha Mountains (Filipova-Marinova 2003a), that is also reflected in marine sediments from cores studied. The local restriction of Ulmus is probably connected to the Fagus expansion which may have out-competed Ulmus within the local forest stand. The decrease of Ulmus and the increase of Carpinus orientalis is due to the destructive human activities during the Early Bronze Age (3200–2600 yr BP) (Todorova 1986) that is confirmed by the increase of anthropogenic indicators, such as Cerealia-type, Triticum, Plantago lanceolata and Polygonum aviculare.

Regional Pollen Assemblage Zone IX

(QuercusUlmusAlnusCarpinus betulusSalixFagus)

2837 cal. yrs BP – pre-industrial time

This RPAZ can be correlated with the Subatlantic chronozone of the Late Holocene and coincides with the New Black Sea Substage. It is represented in cores AKAD 09-10 (LPAZ AKAD 09-10 6) (Fig. 5, Tables 2, 3) and AKAD 11-17 (LPAZ AKAD 11-17 6) (Fig. 3, Tables 2, 3). There is a slight reduction of pollen of deciduous arboreal taxa although they are still dominant. Modern plant communities began to form after 2837 cal. yrs BP. The decrease in mixed oak and hornbeam forests at the beginning of Subatlantic is probably due to the human impact during the Iron Epoch, that is also confirmed by the persistent presence of Cerealia-type and Triticum, as well as anthropophytes Plantago lanceolata and Polygonum aviculare. In areas with erosion, mixed oak forests were replaced by communities of Carpinus orientalis and Quercus pubescens. The increase of Corylus may also reflect deforestation. However, climate-driven changes cannot be excluded. Prominent cooling and increase in humidity are recorded in northern Europe from 2800 to 2500 cal. yrs BP that broadly coincides with the transition of Subboreal to Subatlantic. According to van Geel and Berglund (2000), the cooling during the last 3000 years can be related to the abrupt decrease in solar activity around 2900–2800 cal. yrs BP. Probably the increase of humidity and cooling of the climate along the Black Sea coast were the main reasons for the specific succession and increased abundance of Alnus, Fraxinus excelsior-type and Salix along with lianas and formation of flooded riparian forests (‘Longoz’) lining the river valleys along the Black Sea coast, dominated nowadays by Alnus glutinosa, Fraxinus oxycarpa, Ulmus minor, Carpinus betulus and Quercus pedunculiflora.

Dinocyst assemblages

For a more correct and detailed reconstruction of the natural environment of the Black Sea, the changes in the relative abundance of dinoflagellate species and their assemblages recorded in three cores from the NW Black Sea were studied (Fig. 1) since marine cores are very sensitive and indicative of changes in environmental parameters, such as sea surface temperature (SST), sea surface salinity (SSS) and nutrient variability (Popescu et al. 2009). Two of the cores studied AKAD 09-10 (Fig. 6, Table 2) from the continental slope and AKAD 11-17 (Fig. 4, Table 2) from the deep-water zone, provide long and well-dated continuous Late Quaternary records from the Black Sea spanning the Pleniglacial, Late Glacial and Holocene, while Core AKAD 09-15 (Fig. 8, Table 2) from the peripheral shelf zone represents only a Late Glacial record due to removal of Holocene sediments by strong bottom currents (Dimitrov 1979). Two distinct local dinoflagellate cyst (dinocyst) assemblage zones (LDAZ) were recorded: one dominated by endemic stenohaline freshwater/brackish-water species Pyxidinopsis psilata and Spiniferites cruciformis (LDAZ 1) and a successive one composed of euryhaline brackish to hypersaline species (LDAZ 2). The second zone is subdivided into two subzones (LDASZ 2a and LDASZ 2b). Variations in dinocyst assemblages from the studied cores are very similar and comparable to those studied previously from the western Black Sea for the same period (Wall et al. 1973; Wall and Dale 1974; Mudie et al. 2001, 2002, 2004, 2011; Filipova-Marinova 2006b; Marret et al. 2009; Verleye et al. 2009; Shumilovskikh et al. 2013).

Regional dinocyst assemblage zone 1 (RDAZ 1)

(Pyxidinopsis psilataSpiniferites cruciformis)

25903 – 7668 cal. yrs BP

This zone comprises sediments deposited between 25903 and 7668 cal. yrs BP and could be correlated to the Late Pleniglacial, Late Glacial and Early Holocene (Neweuxinian stage). This assemblage is represented in all three cores studied: AKAD 11-17 from 19546 to 7668 cal. yrs. BP (230–59 cm) (Fig. 4, Table 2); AKAD 09-10 from 25903 to 7663 cal. yrs BP (240–128 cm) (Fig. 6, Table 2); AKAD 09-15 from 14147 to 10906 cal. yrs BP (377–10 cm) (Fig. 8, Table 2). The zone is characterised by a very low number of species and the dominance of two stenohaline freshwater/brackish-water dinoflagellate species Pyxidinopsis psilata and Spiniferites cruciformis accompanied by fresh-water algal taxa, represented by coenobia of the colonial chlorococcalean algae Pediastrum and rarely Botryococcus and Achomosphaera cf. andalusiensis. The spherical Pyxidinopsis psilata is a well-known brackish-water species and an important indicator of cooler SST during glacial periods (Matthiessen and Brenner 1996). However, Brenner (2001) found that occurrence and abundance of P. psilata of annual varves of Baltic Sea sediments is controlled by additional as yet unknown factors other than salinity. Based on modern distribution of P. psilata, Zonneveld and Pospelova (2015) considered this species as a euryhaline temperate one with wide environmental parameter range considering SST and SSS. The ecological affinities and the significance of Spiniferites cruciformis as a crucial indicator of sea surface salinity (SSS) has been an ongoing debate not only recently, but since its discovery. It has been defined classically as a fresh-water/brackish-water species, whose salinity appears not to have exceeded 7‰ (Wall and Dale 1974). Kouli et al (2001) identified S. cruciformis in lacustrine Late Glacial sediments and suggested that S. cruciformis is essentially a freshwater taxon with morphological variations of the cysts that only partly may be linked to salinity variations. Eaton (1996) suggests that the cruciform shape of S. cruciformis must be a stress effect of low salinity waters of less than 12‰. In modern surface sediments, S. cruciformis and P. psilata have been recorded only in the Caspian and Aral Seas at SSS of around 12–13‰ (Marret et al. 2004) and P. psilata in the Baltic Sea at salinity of 12‰ (Yu and Berglund 2007). Based on analogous correlations, most of the authors suggest a maximum SSS of around 12‰ during the lacustrine stage of the Black Sea (Wall et al. 1973; Wall and Dale 1974; Mudie et al. 2001, 2002, 2004; Filipova-Marinova 2006b; Marret et al. 2009; Verleye et al. 2009; Shumilovskikh et al. 2013). The data from our cores revealed that, in addition to the dominant stenohaline brackish-water dinoflagellate species P. psilata and S. cruciformis, the well-known freshwater indicators, such as Pediastrum boryanum var. boryanum, P. simplex var. sturmii, P. kawraskyi, P. integrum and rare Botryococcus, A. cf. andalusiensis and acritarch Multiplicasphaeridium-type, are also present even in modest percentages (note that dinocysts and algae percentages were calculated based on the pollen sum). The investigation of ecological affinities of Pediastrum indicates that this taxon is common to rare at a salinity of about 5‰ (Hiscott et al. 2007) although records from brackish habitats are also documented (Brenner 2001). Detailed studies of surface sediment samples from the Baltic Sea show that the dominant species P. boryanum var. boryanum and P. kawraskyi are stenohaline species associated with a wide range of temperature and trophic conditions, dominant in salinity from 6–8‰ (a full range is 5–9‰), P. simplex occurring in sediments of less than 3–5‰ salinity while Botryococcus is a more south-tolerant taxon sometimes living in brackish waters (Matthiessen and Brenner 1996). Achomosphaera cf. andalusiensis is a species of cold north-temperate climate and low salinity melt-water fluvial conditions (Harland 1983). All these data support our estimation that the limit of SSS for the Neweuxinian stage of the Black Sea falls in the range of 5–12‰.

The dinocyst record shows the low abundance or absence of stenohaline brackish-water species P. psilata from 19 to 15.5 ka yrs BP in core AKAD 11-17 and to 17.1 ka yrs BP in core AKAD 09-10 probably connected with the lack of a favourable environment for the growth of this species due to the strong melt-water and terrigenous input during that time. Pediastrum boryanum var. boryanum (considered a good indicator of freshwater input) has cyclical abundance associated with the deposition of four red-brown clay layers between 19 and 14 ka BP. The red-brown layers have been previously distinguished and dated from the north-western Black Sea shelf from 18.3 to 15.5 ka BP, in several pulses (Major et al. 2002; Bahr et al. 2006). Four red-brown layer intervals were identified and dated by Soulet et al (2016) from 17.2 to 14.8 ka BP. From isotopic and geochemical proxies, Denton et al. (2010) argue that the deposition of the red-brown layers was linked to major melting phases of European ice, drained by the north-western Danube and Dnieper Rivers in response to the climate warming observed after the end of the Last Glacial Maximum (LGM) from 21 to 18 ka BP. Bahr et al. (2005) also indicate four major episodes equivalent to four layers of Caspian water spilling into the Black Sea. According to Major et al. (2002), it is difficult to support a direct northern source for the red-brown layers associated with the European melt-water event. During the time of the red-brown layers’ deposition, the Scandinavian Ice sheet dammed large lakes and diverted rivers and melt-water drainage southwards towards the River Volga and into the Caspian Sea. These melt-waters caused the Caspian Sea to overflow over the Manych Depression (+20 m a.s.l.) into the Black Sea (Mangerud et al. 2004). Bahr et al. (2005) assume that these sediments have been transported to the shelf and the slope by strong surface currents such as the present-day ‘Rim Current’ gyre. The occurrence of acritarch Glomus-type 207 in the sediments provides evidence for terrigenous input and erosional rate increase during the deposition of these red-brown layers. The decrease in the freshwater algae Pediastrum is marked between 14259 and 14036 cal. yrs BP. According to its stratigraphic position at the Older Dryas sediments, it may reflect a drier with reduced river input and/or cooler climate at that time. The increase of P. psilata and S. cruciformis during the interstadials Bølling and Allerød, identified in dinocyst record of core AKAD 11-17, indicates increased phytoplankton activity and rather brackish Black Sea waters during that time, characterised by a marked increase in temperate arboreal pollen suggesting warmer and humid conditions along the Black Sea coast. The upward increase of P. psilata, S. cruciformis and P. boryanum var. boryanum between 13257 and 11072 cal. yrs BP indicates Younger Dryas SST colder than during the Bølling/Allerød in agreement with the pollen record of this core and other cores from the surrounding regions (Bottema et al. 1995; Mudie et al. 2002; Filipova-Marinova 2003b), as well as with the return to almost LGM values of δ18Oprec established by Bahr et al. (2006). The maximum relative abundance of P. psilata up to 94% in the Early Holocene record from core AKAD 11-17 is of considerable interest. Similar abundance of P. psilata was also found in most records from the north-western Black Sea, for example, 93–97% in core MO2-45 (Marret et al. 2009); 91–97% in core 22-GCC3 (Shumilovskikh et al. 2013); 92% in core GeoB7625-2 (Verleye et al. 2009). The explanation of why P. psilata considered a cool brackish-water species shows maximum relative abundance during the Early Holocene is complicated. Preboreal and Boreal chronozones of the Early Holocene are characterised by warming and humidity increase, that is previously stated in pollen records (Komarov et al. 1979; Mudie et al. 2007; Filipova-Marinova et al. 2013) and in the varved δ18O paleo-term records from the Sofular Cave that tracks the isotope signature of Black Sea surface water (Göktürk et al. 2011). During the Early Holocene, characterised by a warmer climate, the stenohaline freshwater/brackish-water species P. psilata that have a preference for warmer temperatures and ice-free conditions (Zonneveld and Pospelova 2015), together with S. cruciformis, demonstrate their ecological optimum in growth concerning SST reaching maximum relative abundance between 11072 and 8638 cal. yrs BP. This maximum was interrupted by an abrupt significant short-term decrease in the relative abundance of P. psilata which reached values close to those characteristic of glacial periods (from 84–20%), as well as S. cruciformis (from 21–3%). This shift, centring between 8500 and 8300 cal. yrs BP in the deep-water core AKAD 11-17, is considerable and reflects a strong cooling of climate. This episode could be attributed to a marked short-term (with a duration of about 200 years) climate deterioration against the background of the climate amelioration during the Early Holocene. This event was also revealed in the same core by pollen data, particularly based on ecological requirements of the vegetation, represented by the ratio of relative abundance of thermophilous arboreal taxa and steppe elements including cold-resistant and heterophyllous taxa. The pollen-driven decrease in temperatures correlates with the decrease in phytoplankton productivity. The same cooling phase was previously assigned in sequence of GGC-18 from the south-western Black Sea continental slope indicating the complexity of environmental changes throughout the Basin (Filipova-Marinova et al. 2013) and confirming that the climate of the Early Holocene acts as the strongest controlling factor (Brewer et al. 2002). Bahr et al. (2006) also established that prevailing high temperatures during the Early Holocene led to authigenic calcite precipitation through increased phytoplankton activity, interrupted by the ‘8.2 ka cold event’. The findings reported above represent the regional expression of the well-known in the Northern Hemisphere ‘8.2 ka cold event’ in the Greenland δ18O ice which is commonly linked to a melt-water-related perturbation of the Atlantic Meridional Overturning Circulation (AMOC) (Magny et al. 2003; Alley and Ágústsdóttir 2005; Bahr et al. 2005; Le Grande et al. 2006) and associated collapse of oceanic northward heat transport (Ellison et al. 2006). This climate anomaly is also confirmed in the Aegean Region by multi-proxy studies. Early to Middle Holocene season-specific SST and δ12O seawater, based on dinocyst and foraminiferal records established in marine sediments revealed a prominent short-term (approximately 150 yrs) cooling event in the central Aegean Sea centred on 8200 cal. yrs BP, coeval to the 8.2 ka BP cold event (Marino et al. 2009). Pollen record from Thenaghi Philippon (NE Aegean Sea) provides evidence for a massive climate-induced turnover in terrestrial ecosystems of the Aegean Region associated with the 8.2 cal. ka BP cold event as a response to the North-Atlantic thermohaline circulation slowdown (Pross et al. 2009). Using a phytoplankton-based SSS record, based on the Emiliania huxleyi transfer function, Herrle et al. (2018) found a reduced outflow of low salinity waters from the Black Sea into the north Aegean Sea corresponding to the 8.2 cal. ka BP cold event that was caused by cooler and drier climatic conditions over Europe. All these data support the statement that high amplitude temperature anomaly linked to the 8.2 cal. ka cold event may have also occurred in lower latitudes much further in the southerly direction than previously assigned. The Black Sea like the Aegean Sea is sufficiently distant and isolated from the North Atlantic not to be directly affected by its oceanic circulation, therefore, signal transmission must have been atmospheric (Marino et al. 2009). The 8.2 cal. ka cold event in many ways mimics the distribution and anomaly type of the Younger Dryas. It is marked, in general, by cooling with summer-time values in the order of 1°C by a strong shift in vegetation (Alley and Ágústsdóttir 2005).

Regional dinocyst assemblage subzone 2a (RDASZ 2a)

(Lingulodinium machaerophorumSpiniferites beleriusSpiniferites bentorii)

7668 – 2837 cal. yrs. BP

A prominent change of the composition of dinocyst assemblages from freshwater/brackish-water is observed at the boundary of LDAZ AKAD 11-17-1 and LDAZ AKAD 11-17-2a. The abrupt decrease of stenohaline freshwater/brackish-water species S. cruciformis and P. psilata at 7668 cal yrs BP indicating, according to Mudie et al. (2001), the inability of these species to survive the abrupt salinity change to values as high as 10–12‰ (Deuser 1972) or even 18‰ (Wall and Dale 1974) was followed upward by a gradual increase of euryhaline marine species. The dinoflagellate cyst assemblage was dominated by euryhaline marine species, such as Lingulodinium machaerophorum, Spiniferites belerius, Spiniferites bentorii and acritarch Cymatiosphaera globulosa and is characterised by high species diversity. The change of assemblage occurs coincidently with the lithological change from banded lutite to sapropel which has been dated also at 7668 cal. yrs BP. Modern marine conditions, influenced by Mediterranean waters, were detected after 6417 cal. yrs BP, when an abundance of Mediterranean-related species Operculodinium centrocarpum, Spiniferites mirabilis and acritarch Mychristidium (occurring in hypersaline waters) along with euryhaline heterotrophic species, such as Brigantedinium cariacoense, Echinidinium transparantum and Polykrikos kofoidii, occurred. The immigration of marine species to the north-western Black Sea zone from 7668 to 6417 cal. yrs BP indicates a progressive increase in salinity due to more open connection with the Mediterranean. All these species are characteristic components of climate amelioration in Quaternary sequences. Such dinoflagellate cyst assemblages tolerate salinity values of 17–9‰ (Deuser 1972). The isotopic salinity estimates are in close agreement with the present day SSS of 18‰ (Mudie et al. 2002). The recorded sporadic occurrence of acritarchs, such as Hexasteria problematica and Fungal spores Sordaria-type, Cercophora-type 207 and Ascospores-type 121 in the studied records, appeared to be the best index of terrigenous input from soil erosion and markers of farming practice (Mudie et al. 2002). The most abundant euryhaline species L. machaerophorum responds to the amelioration of climate and shows three peaks in the record of core AKAD 11-17: at 5016, 3459 and 1436 cal. yrs BP and in the core AKAD 09-10 at 220, 190 and 140 cm. These peaks could be related to the increased SST, SSS, a stratified water column (for instance, calm periods between upwelling events) and nutrient input due to the increased human activity during the Late Eneolithic, Early Bronze Age and Iron epoch along the Black Sea coast, confirmed also by pollen data. According to Mudie et al. (2011), blooms of L. machaerophorum cysts may indicate water-column stratification and nutrient enrichment more strongly than either salinity or sea level. Peaks of another characteristic component of climate amelioration, such as Cymatiosphaera globulosa, are synchronous with those of L. machaerophorum and also reflect mainly nutrient levels and stratification of the water column. The first appearance of Peridinium ponticum, a species related geographically to the Black Sea, is recorded from 6417 cal. yrs BP. This species is considered a good proxy for the reconstruction of Holocene salinity variations since its relative abundance fluctuates asynchronously with L. machaerophorum with long processes (Mudie et al. 2011).

Regional dinocyst assemblage subzone 2b (LDASZ 2b)

(Lingulodinium machaerophorumSpiniferites ramosusPeridinium ponticum)

2837 cal. yrs. BP – pre-industrial time

The composition of the Late Holocene LDASZ AKAD 11-172b is identical to that of LDASZ AKAD 11-17-2a, but the decrease in the relative abundance of almost all marine euryhaline species since 2837 cal. yrs BP is noticeable. The increase in abundance of species that tolerate low-salinity water conditions, such as P. ponticum and S. ramosus, indicate certain freshening. Since P. ponticum is a heterotrophic species (Dale 1996), nutrient availability also influences its abundance (Verleye et al. 2009) as recorded by pollen data, particularly by the increase inhuman activity indicators. The sporadic occurrence of freshwater algal species P. boryanum var. boryanum and P. simplex var. sturmii, as well as Spirogyra and A. cf. andalusiensis, suggests increasing river input associated with the cooling of climate and increase of humidity confirmed by pollen data from the studied cores. The period from 2800 to 2500 cal yrs BP coincides with a prominent cooling and increase of humidity in north Europe (Wanner et al. 2008). According to van Geel and Berglund (2000), the cooling in the Third millennium BP can be related to the abrupt decrease in solar activity around 2900 to 2800 cal. yrs BP. Freshening of the Black Sea surface water during the Late Holocene after 2500 cal yrs BP was previously assigned in the Core GGC-18 of the south-western Black Sea by Filipova-Marinova et al. (2013). Decrease in the abundance of fully marine dinoflagellate species and occurrence of the brackish-water indicator P. ponticum indicates seawater freshening after 2500 cal. yrs BP, also recorded in core 22-GC3 (Shumilovskikh et al. 2013). Substantial freshening of the Black Sea surface waters from 29 to approximately 19‰ after 3000 cal. yrs BP is evidenced by lower δD values of C37 alkenones produced by haplophyte algae Emiliania huxleyi (van der Meer et al. 2008). Additionally, combined lipid biomarkers and fossil DNA analyses of alkenones of E. huxleyi from the GGC-18 sediment core indicate a gradual cooling from 19–15 °C after 2570 cal. yrs BP (Coolen et al. 2009).

Conclusions

Vegetation successions and environmental changes along the north-western Black Sea coastal area during the last 26000 years were reconstructed by multi-proxy analysis including radiocarbon dating of sediments from three new marine cores. The following main conclusions from this study are: (1) The coastal landscape during the Late Pleniglacial (25903–17092 cal. yrs BP) was dominated by steppe vegetation composed of Artemisia, Chenopodiaceae, Poaceae and other cold-resistant and heliophilous herbs suggesting cold and dry environments. Sparse stands of Pinus and Quercus, partly enlarged during the melting pulses (19.2–14.5 cal. ka BP) and during the Late Glacial interstadials Bølling and Allerød reflecting warming and humidity increase. (2) During the Younger Dryas (13257–11788 cal. yrs BP), enlargement of steppe vegetation dominated by Artemisia and the shrubland of Juniperus-Ephedra indicates return to the coldest and driest climate. (3) In the Early Holocene (Preboreal-Boreal) (11788–8004 cal. yrs BP), pioneer forests of Quercus with groups of Ulmus, Tilia, Alnus and Betula spread and clearly confirm the presence of refugia of these taxa in the coastal mountains and their rapid migration due to the climate warming. (4) The short-term decline of arboreal pollen, particularly manifested by Quercus, between 8.5–8.3 ka BP can be explained as a vegetation response to the known in north Atlantic region ‘8.2 ka cold event’. This climatic oscillation is confirmed for the second time in Black Sea sediments. (5) During the Atlantic chronozone (8004–5483 cal. yrs BP), species-rich mixed oak temperate deciduous forests developed in the lowlands following climate optimal conditions (high humidity and increased mean annual temperatures). (6) During the Subboreal chronozone (5483–2837 cal. yrs BP), mixed oak forests dominate alongside a slight enlargement of Carpinus betulus. (7) During the Subatlantic chronozone (2837 cal. yrs BP – pre-industrial time), a specific vegetation succession manifested by the increased abundance of Alnus, Fraxinus excelsior and Salix along with lianas and the formation of flooded riparian forests (e.g. ‘Longoz’) lining the river valleys along the Black Sea suggests a climate shift (an increase of humidity and a cooling of the climate). (8) The first indications of farming and other human activities along the Black Sea coast were recorded during the Late Eneolithic (6790–6320 yrs BP). (9) Two main dinoflagellate cyst assemblages were distinguished: one dominated by stenohaline freshwater/brackish-water species and the successive one dominated by euryhaline marine species. (10) During the Early Holocene, Pyxidinopsis psilata revealed a wide ecological range and demonstrated its ecological optimum of growth concerning the increased sea surface temperature reaching a maximum relative abundance at 9475 cal yrs BP. (11) An abrupt short-term cooling centred between 8.5 and 8.3 ka BP associated with the ‘8.2 ka cold event’ is evidenced by an abrupt decline in the abundance of P. psilata and Spiniferites cruciformis to values close to those characteristic for the Younger Dryas. This climate oscillation is described for the first time in dinocyst records from Black Sea sediments (Core Akad 11-17). This finding confirms that the high amplitude temperature anomaly, associated with the ‘8.2 ka cold event’ may have also occurred in a southern direction possibly through atmospheric transmission of the signals. (12) The change in the composition of dinocyst assemblages occurred at 7668 cal yrs BP. The abrupt disappearance of freshwater/brackish-water species Pyxidinopsis psilata and Spiniferites cruciformis was followed upwards by a gradual increase in euryhaline marine species Lingulodinium machaerophorum, Spiniferites belerius, S. bentorii and acritarch Cymatiosphaera globulosa. (13) A certain freshening of the Black Sea waters after 2837 cal. yrs BP has been established.

Additional information

Conflict of interest

The authors have declared that no competing interests exist.

Ethical statement

No ethical statement was reported.

Funding

The studied gravity cores Akad 09-10, Akad 09-15 and Akad 11-17 were collected by an international scientific team during both expeditions, with financial support by the Bulgarian National Science Fund within the Project DO 02-337, led by Prof. Petko Dimitrov. The study was partly supported by the Project “Upgrading of distributed scientific infrastructure – Bulgarian Network for Long-Term Ecosystem Research” (LTER-BG), (agreement with Ministry of Education and Science, DO1-163/28.07.2022).

Author contributions

All authors have contributed equally.

Author ORCIDs

Mariana Filipova-Marinova https://orcid.org/0000-0002-0786-9476

Danail Pavlov https://orcid.org/0000-0001-7382-2054

Krasimira Slavova https://orcid.org/0000-0002-0622-8490

Data availability

All of the data that support the findings of this study are available in the main text.

References

  • Aksu A, Hiscott R, Kaminski M, Mudie P, Gillespie H, Abrajano T, Yaşar D (2002) Last glacial–Holocene paleoceanography of the Black Sea and Marmara Sea: Stable isotopic, foraminiferal and coccolith evidence. Marine Geology 190(1–2): 119–149. https://doi.org/10.1016/S0025-3227(02)00345-6
  • Alexiev G (2002) The Black Sea. In: Kopralev I (Ed.) Geography of Bulgaria. ForKom Publishers, Sofia, 90–92 . [in Bulgarian]
  • Atanassova J (1990) Vegetation Development during the Late Quaternary based on Spore-Pollen Analyses of Sediments from the Western Section of the Black Sea. PhD Thesis, Sofia University, Bulgaria. [in Bulgarian]
  • Atanassova J (1995) Palynological data of three deep water cores from the western part of the Black Sea. In: Bozilova E, Tonkov S (Eds) Advances in Holocene palaeoecology in Bulgaria, Pensoft Publishers, Sofia-Moscow, 68–83.
  • Atanassova J (1999) Palaeoecological investigations of Late Quaternary sediments from the western part of the Black Sea. Acta Palaeobotanica 2: 545–548.
  • Atanassova J, Bozilova E (1992) Palynological investigations of marine sediments from the western sector of the Black Sea. Proceedings of the Institute of Oceanology (Varna) 1: 97–103.
  • Bahr A, Arz HW, Lamy F, Wefer G (2006) Late glacial to Holocene paleoenvironmental evolution of the Black Sea, reconstructed with stable oxygen isotope records obtained on ostracod shells. Earth and Planetary Science Letters 241(3–4): 863–875. https://doi.org/10.1016/j.epsl.2005.10.036
  • Bahr A, Lamy F, Arz H, Major C, Kwiecien O, Wefer G (2008) Abrupt changes of temperature and water chemistry in the late Pleistocene and early Holocene Black Sea. Geochemistry, Geophysics, Geosystems 9(1): 1–16. https://doi.org/10.1029/2007GC001683
  • Bennett K, Tzedakis P, Willis K (1990) Quaternary refugia of north European trees. Journal of Biogeography 18(1): 103–115. https://doi.org/10.2307/2845248
  • Berglund B (1983) Palaeoecological studies in lakes and mires – a palaeoecological research strategy. In: Gregory K (Ed.) Background to Palaeoecology. John Wiley & Sons, Chichester, 237–249.
  • Berglund B, Lemdahl G, Liedberg-Jönsson B, Persson T (1984) Biotic response to climatic changes during the time span 13,000-10,000 BP – a case study from SW Sweden. In: Mörner NA, Karlén W (Eds) Climatic changes on a Yearly to Millennial basis. D. Reidel Publishing Company, Dordrecht, 25–36. https://doi.org/10.1007/978-94-015-7692-5_3
  • Berglund B, Bjorck S, Lemdahl G, Bergsten H, Nordberg K, Kolstrup E (1994) Late Weichselian environmental change in southern Sweden and Denmark. Journal of Quaternary Science 9(2): 127–132. https://doi.org/10.1002/jqs.3390090206
  • Beug HJ (1961) Leitfagen der Pollenbestimmung. Gustav Fischer Verlag, Stuttgart, 63 pp.
  • Beug HJ (1975) Man as a factor in the vegetational history of the Balkan Peninsula. Proceedings, First International Symposium on Balkan Flora and Vegetation. Varna, Bulgaria, 72–78.
  • Beug HJ (1982) Vegetation history and climatic changes in central and southern Europe. In: Harding A (Ed.) Climatic change in later prehistory. Edinburgh University Press, Edinburgh, 85–102.
  • Beug HJ (2004) Leitfaden der Pollenbestimmung für Mitteleuropa und angrenzende Gebiete. Verlag Dr. Friedrich Pfeil, München, 542 pp.
  • Birks HJ (1973) Past and Present vegetation of the Isle of Skye – a palaeoecological study. Cambridge University Press, Cambridge, 427 pp.
  • Birks HJB, Birks HH (1980) Quaternary Palaeoecology. Edward Arnold, London, 289 pp.
  • Birks H, Line J (1993) Glacial refugia of European trees – a matter of chance? Dissertationes Botanicae 196: 283–291.
  • Blytt A (1876) Essay on the Immigration of the Norwegian Flora During the Alternating Rainy and Dry Period, Alb. Cammermayer, Christiana, Oslo.
  • Bondev I (1991) Vegetation of Bulgaria: 1:600,000 map with explanatory text. “St. Kliment Ohrdski” University Press, Sofia, 184 pp. [in Bulgarian]
  • Bottema S, Woldring H, Aytuğ B (1995) Late Quaternary vegetation history of northern Turkey. Palaeohistoria 35/36: 13–72.
  • Bozilova E, Beug HJ (1992) On the Holocene history of vegetation in SE Bulgaria (Lake Arkutino, Ropotamo region). Vegetation History and Archaeobotany 1(1): 19–32. https://doi.org/10.1007/BF00190698
  • Bozilova E, Beug HJ (1994) Studies on the vegetation history of Lake Varna Region, Northern Black Sea coastal area of Bulgaria. Vegetation History and Archaeobotany 3(3): 143–154. https://doi.org/10.1007/BF00202022
  • Bozilova E, Filipova M (1991) Palynological and palaeoethnobotanical evidence about the human impact on the vegetation along the Black Sea coast from the Neolithic till the Greek Colonisation. Thracia pontica 4: 87–96.
  • Bozilova E, Tonkov S (2000) Pollen from Lake Sedmo Rilsko reveals southeast European postglacial vegetation in the highest mountain area of the Balkans. The New Phytologist 148(2): 315–325. https://doi.org/10.1046/j.1469-8137.2000.00754.x
  • Bozilova E, Filipova M, Dimitrov P (1979) Pollen analysis of Late Quaternary sediments from the peripheral part of the Black Sea. Bulletin du Musée National de Varna 15(30): 157–162. [in Russian]
  • Brenner WW (2001) Organic-walled microfossils from the central Baltic Sea, indicators of environmental change and base for ecostratigraphic correlation. Baltica 14: 40–51.
  • Chepalyga A (2002) Marine basins. The Black Sea. In: Velichko A (Ed.) Dynamics of Terrestrial Landscape Components and Inner Marine Basins of Northern Eurasia during the Last 130 000 Years. GEOS Publishers, Moscow, 170–182. [in Russian]
  • Chernyishova M (1980) Palynological investigations of sediments on the continental terrace. In: Kuprin P (Ed.) Geological and Geophysical researches in the Bulgarian Black Sea zone. Publishing house of the Bulgarian Academy of Sciences, Sofia, 213–222. [in Russian]
  • Combourieu Nebout N, Peyron O, Dormoy I, Desprat S, Beaudouin C, Kotthoff U, Marret F (2009) Rapid climatic variability in the west Mediterranean during the last 25000 years from high resolution pollen data. Climate of the Past 5(3): 503–521. https://doi.org/10.5194/cp-5-503-2009
  • Connor SE, Thomas I, Kvavadze EV (2007) A 5600-year history of changing vegetation, sea levels and human impacts from the Black Sea coast of Georgia. The Holocene 17(1): 25–36. https://doi.org/10.1177/0959683607073270
  • Connor SE, Ross SA, Sobotkova A, Herries AIR, Mooney SD, Longford C, Iliev I (2013) Environmental conditions in the SE Balkans since the Last Glacial Maximum and their influence on the spread of agriculture into Europe. Quaternary Science Reviews 68: 200–215. https://doi.org/10.1016/j.quascirev.2013.02.011
  • Coolen MJL, Saenz JP, Giosan L, Trowbridge NY, Dimitrov P, Dimitrov D, Eglinton TI (2009) DNA and lipid molecular stratigraphic records of haptophyte succession in the Black Sea during the Holocene. Earth and Planetary Science Letters 284(3–4): 610–621. https://doi.org/10.1016/j.epsl.2009.05.029
  • Cordova C, Harrison S, Mudie P, Riehl S, Leroy S, Ortiz N (2009) Pollen, plant macrofossil, and charcoal records for paleovegetation reconstruction in the Mediterranean-Black Sea corridor since the Last Glacial Maximum. Quaternary International 197(1–2): 12–26. https://doi.org/10.1016/j.quaint.2007.06.015
  • Dale B (1996) Dinoflagellate cyst ecology: Modeling and geological applications. In: Jansonius J, McGregor DC (Eds) Palynology: Principles and Applications. AASP Foundation, Salt Lake City, 1249–1275.
  • Davis B, Brewer S, Stevenson A, Guiot J (2003) The temperature of Europe during the Holocene reconstructed from pollen data. Quaternary Science Reviews 22(15–17): 1701–1716. https://doi.org/10.1016/S0277-3791(03)00173-2
  • De Vernal A, Turon JL, Guiot J (1994) Dinoflagellate cyst distribution in high-latitude marine environments and quantitative reconstruction of sea-surface salinity, temperature, and seasonality. Canadian Journal of Earth Sciences 31(1): 48–62. https://doi.org/10.1139/e94-006
  • Deuser W (1972) Late Pleistocene and Holocene history of the Black Sea as indicated by stable isotope studies. Journal of Geophysical Research 77(6): 1071–1077. https://doi.org/10.1029/JC077i006p01071
  • Dimitrov P (1979) Peculiarities of the composition and distribution of bottom sediments in the Black Sea Shelf between Cape Kaliakra and Cape Emine. Oceanologia 3: 23–33. [in Bulgarian]
  • Dormoy I, Peyron O, Combourieu Nebout N, Goring S, Kotthoff U, Magny M, Pross J (2009) Terrestrial climate variability and seasonality changes in the Mediterranean region between 15000 and 4000 years BP deduced from marine pollen records. Climate of the Past 5(4): 615–632. https://doi.org/10.5194/cp-5-615-2009
  • Ellison CRW, Chapman MR, Hall IR (2006) Surface and deep ocean interactions during the cold climate event 8200 years ago. Science 312(5782): 1929–1932. https://doi.org/10.1126/science.1127213
  • Faegri K, Iversen J (1989) Textbook of pollen analysis. John Wiley & Sons, Chichester, 338 pp.
  • Feurdean A, Klotz S, Mosbrugger V, Wohlfarth B (2008) Pollen-based quantitative reconstructions of Holocene climate variability in NW Romania. Palaeogeography, Palaeoclimatology, Palaeoecology 260(3–4): 494–504. https://doi.org/10.1016/j.palaeo.2007.12.014
  • Feurdean A, Willis KJ, Parr C, Tantau I, Farcas S (2010) Postglacial patterns in vegetation dynamics in Romania: Homogenization or differentiation? Journal of Biogeography 37(11): 2197–2208. https://doi.org/10.1111/j.1365-2699.2010.02370.x
  • Filipova M, Dimitrov P (1987) Palynological and stratigraphical data about the Quaternary of the Bay of Balchik. Bulletin du Musée National de Varna 23: 215–229. [in Bulgarian]
  • Filipova M, Bozilova E, Dimitrov P (1989) Palynological investigation of the Late Quaternary deep-water sediments from the Southwestern part of the Black Sea. Bulletin du Musée National de Varna 25: 177–186.
  • Filipova-Marinova M (2003a) Postglacial vegetation dynamics in the coastal part of the Strandza Mountains, Southeastern Bulgaria. In: Tonkov S (Ed.) Aspects of Palynology and Palaeoecology: Festschrift in honour of Elissaveta Bozilova. Pensoft Publishers, Sofia-Moscow, 213–231.
  • Filipova-Marinova M (2003b) Palaeoenvironmental changes in the Southern Bulgarian Black Sea area during the last 29,000 years. Phytologia Balcanica 9: 275–292.
  • Filipova-Marinova M (2006a) Palynostratigraphy of Pleistocene and Holocene sediments from the western Black Sea. Journal of Environmental Micropalaeontology. Microbiology and Meiobenthology 3: 91–112.
  • Filipova-Marinova M (2006b) Late Pleistocene/Holocene dinoflagellate cyst assemblages from the Southwestern Black Sea shelf. In: Ognjanova-Rumenova N, Manoylov K (Eds) Advances in Phycological Studies Festschrift in Honour of Prof. Dobrina Temniskova-Topalova. Pensoft Publishers, Sofia, 267–281.
  • Filipova-Marinova M (2007) Archaeological and paleontological evidence of climate dynamics, sea-level change, and coastline migration of the Bulgarian sector of the Circum-Pontic region. In: Yanko-Hombach V, Gilbert A, Panin N, Dolukhanov P (Eds) The Black Sea Flood Question: Changes in Coastline, Climate, and Human Settlement, Springer, Dordrecht, 453–482. https://doi.org/10.1007/978-1-4020-5302-3_19
  • Filipova-Marinova M, Bozilova E (2002) Palaeoecological conditions in the area of the praehistorical settlement in the Bay of Sozopol during the Eneolithic. Phytologia Balcanica 8(2): 133–143.
  • Filipova-Marinova M, Bozilova E (2003) Palaeoecological evidence of the vegetation history and human occupation in the coastal area of Sozopol (Southeastern Bulgaria). Dobroudzha 21: 279–291.
  • Filipova-Marinova M, Christova R (2004) Sea-level fluctuations in the western part of the Black Sea during the Holocene. Journal of Environmental Micropalaeontology. Microbiology and Meiobenthology 1: 122–135.
  • Filipova-Marinova M, Christova R, Bozilova E (2004) Palaeoecological conditions in the Bulgarian Black Sea area during the Quaternary. Journal of Environmental Micropalaeontology. Microbiology and Meiobenthology 1: 136–154.
  • Filipova-Marinova M, Pavlov D, Coolen M, Giosan L (2013) First high-resolution marinopalynological stratigraphy of Late Quaternary sediments from the central part of the Bulgarian Black Sea area. Quaternary International 293: 170–183. https://doi.org/10.1016/j.quaint.2012.05.002
  • Filipovitch L, Stefanova I, Lazarova M, Petrova M (1997) Holocene vegetation in Stara Planina (the Balkan Range) I. Phytologia Balcanica 3: 15–39.
  • Filipovitch L, Stefanova I, Lazarova M, Petrova M (1998) Holocene vegetation in Stara Planina (the Balkan Range). II. Vegetational development and changes in the vertical distribution of forests. Phytologia Balcanica 4(1–2): 13–25.
  • Finsinger WT, van der Knaap WO, Ammann B (2006) The expansion of hazel (Corylus avellana L. ) in the southern Alps: A key for understanding its Early Holocene history in Europe? Quaternary Science Reviews 25(5–6): 612–631. https://doi.org/10.1016/j.quascirev.2005.05.006
  • Fleitmann D, Cheng H, Badertscher S, Edwards R, Mudelsee M, Göktürk OM, Fankhauser A, Pickering R, Raible CC, Matter A, Kramers J, Tüysüz O (2009) Timing and climatic impact of Greenland interstadials recorded in stalagmites from northern Turkey. Geophysical Research Letters 36(19): L19707. https://doi.org/10.1029/2009GL040050
  • Giesecke T, Bennett KD, Birks HJB, Bjune AE, Bozilova E, Feurdean A, Finsinger W, Froyd C, Pokorný P, Rösch M, Seppä H, Tonkov S, Valsecchi V, Wolters S (2011) The pace of Holocene vegetation change - testing for synchronous developments. Quaternary Science Reviews 30(19–20): 2805–2814. https://doi.org/10.1016/j.quascirev.2011.06.014
  • Göktürk OM, Fleitmann D, Badertscher S, Cheng H, Edwards RL, Leuenberger M, Fankhauser A, Tüysüz O, Kramers J (2011) Climate on the southern Black Sea coast during the Holocene: Implications from the Sofular Cave record. Quaternary Science Reviews 30(19–20): 2433–2445. https://doi.org/10.1016/j.quascirev.2011.05.007
  • Grimm E (2011) TGView. Version 1.5.12. Illinois State Museum.
  • Harland R (1983) Distribution maps of recent dinoflagellate cysts in bottom sediments from the North Atlantic Ocean and adjacent seas. Palaeontology 26: 321–387.
  • Hay BJ, Arthur MA, Dean WE, Neff ED, Honjo S (1991) Sediment deposition in the Late Holocene Abyssal Black-Sea with Climatic and chronological implications. Deep-Sea Research. Part A, Oceanographic Research Papers 38: 1211–1235. https://doi.org/10.1016/S0198-0149(10)80031-7
  • Hiscott RN, Aksu AE, Mudie PJ, Marret F, Abrajano T, Kaminski M, Evans J, Cakiroglu A, Yasar D (2007) A gradual drowning of the southwestern Black Sea shelf: Evidence for a progressive rather than abrupt Holocene reconnection with the eastern Mediterranean Sea through the Marmara Sea Gateway. Quaternary International 167: 19–34. https://doi.org/10.1016/j.quaint.2006.11.007
  • Huntley B, Birks H (1983) An atlas of past and present pollen maps for Europe 0–13000 years ago. Cambridge University Press, Cambridge, 668 pp.
  • Izdar E, Ergün M (1991) Recent geological evolution of the Black Sea: an overview. In: Izdar E, Murray JW (Eds) Black Sea Oceanography. NATO ASI Series, Series C: Mathematical and Physical Sciences. Kluwer Academic Publishers, London, 351: 379–387. https://doi.org/10.1007/978-94-011-2608-3_23
  • Jankovska V, Komarek J (2000) Indicative value of Pediastrum and other coccal green algae in palaeoecology. Folia Geobotanica 35(1): 1–59. https://doi.org/10.1007/BF02803087
  • Kerey IE, Meriç E, Tunoğlu C, Kelling G, Brenner RL, Doğan AU (2004) Black Sea–Marmara Sea Quaternary connections: New data from the Bosphorus, İstanbul, Turkey. Palaeogeography, Palaeoclimatology, Palaeoecology 204(3–4): 277–295. https://doi.org/10.1016/S0031-0182(03)00731-4
  • Khrischev K (1984) Sedimentogenesis of South Bulgarian Black Sea Shelf. Dr. Sciences Thesis, Geological Institute of the Bulgarian Academy of Sciences, Sofia. [in Bulgarian]
  • Khrischev K, Shopov V (1978) Quaternary stratigraphy of the Outer South Bulgarian Black Sea shelf. Geologica Balcanica 8(2): 21–40.
  • Kochev H, Jordanov D (1981) Vegetation of Water Basins in Bulgaria. Publishing house of the Bulgarian Academy of Sciences, Sofia, 183 pp. [in Bulgarian]
  • Komarov A, Filipova M, Oudintzeva O (1978) Palynological evidence of stratigraphy of upper Quaternary shelf sediments of the western Black Sea area. Theoretic and methodologic basis of the complex shelf researches. Geological Institute of the Academy of Sciences of USSR, Moscow, 93–95. [in Russian]
  • Komarov A, Bozilova E, Filipova M, Oudintzeva O (1979) Palynological spectra and their stratigraphic interpretation. In: Malovitskiy Y (Ed.) Geology and Hydrology of the Western Section of the Black Sea. Publishing House of the Bulgarian Academy of Sciences, Sofia, Bulgaria, 85–91. [in Russian]
  • Koreneva E (1980) Palynological investigations of the Black Sea on deep-water coring data. In: Zaklinskaia E, Neisthadt M, Medvedeva A, Rovnina L (Eds) Palynology in USSR. Nauka Publishers, Moscow, USSR, 103–104. [in Russian]
  • Koreneva E, Kartashova G (1978) Palynological study of samples from holes 379A, 380A, Leg 42B. In: Ross D, Neprochnov P (Eds) Initial Reports of Deep-Sea Drilling Project, v. 42 (2). U.S. Government Printing Office, Washington, USA, 951–992. https://doi.org/10.2973/dsdp.proc.42-2.143.1978
  • Kosarev A, Kostianoy AG, Shiganova TA (2007) The Sea of Azov. In: Kostianoy A, Kosarev A (Eds) Handbook of Environmental Chemistry, Part 5Q Water Pollution: The Black Sea Environment. Springer Verlag, Berlin – Heidelberg, 63–89. https://doi.org/10.1007/698_5_091
  • Kutzbach J, Webb T (1993) Conceptual basis for understanding Late Quaternary climates. In: Wright HE, Kutzbach J, Webb T, Ruddiman WF, Street-Perrott FA, Bartlein PJ (Eds) Global Climates since the Last Glacial Maximum. Univ. Minnesota Press, Minneapolis-London, 5–11.
  • Lang G (1985) Palynologic and stratigraphic investigations of Swiss lake and mire deposits. A general view over a research program. Dissertationes Botanicae 87: 107–114.
  • Lang G (2003) Immigration of Tilia in Europe since the last Glacial. In: Tonkov S (Ed.) Aspects of Palynology and Palaeoecology: Festschrift in honour of Elissaveta Bozilova. Pensoft Publishers, Sofia-Moscow, 21–41.
  • Le Grande AN, Schmidt GA, Shindell DT, Field CV, Miller RL, Koch DM, Faluvegi G, Hoffmann G (2006) Consistent simulations of multiple proxy responses to an abrupt climate change event. Proceedings of the National Academy of Sciences of the United States of America 103(4): 837–842. https://doi.org/10.1073/pnas.0510095103
  • Leroy S, Arpe K (2007) Glacial refugia for summer-green trees in Europe and south-west Asia as proposed by ECHAM3 time-slice atmospheric model simulations. Journal of Biogeography 34(12): 2115–2128. https://doi.org/10.1111/j.1365-2699.2007.01754.x
  • Londeix L, Herreyre Y, Turon J-L, Fletcher W (2009) Last Glacial to Holocene hydrology of the Marmara Sea inferred from a dinoflagellate cyst record. Review of Palaeobotany and Palynology 158 (1–2): 52–71.Magny M, Begeot C, Guiot J, Peyron O (2003) Contrasting patterns of hydrological changes in Europe in response to Holocene climate cooling phases. Quaternary Science Reviews 22: 1589–1596. https://doi.org/10.1016/j.revpalbo.2009.07.004
  • Magny M, Bégeot C, Guiot J, Peyron O (2003) Contrasting patterns of hydrological changes in Europe in response to Holocene climate cooling phases. Quaternary Science Reviews 22(15–17): 1589–1596. https://doi.org/10.1016/S0277-3791(03)00131-8.
  • Major C, Ryan W, Lericolais G, Hajdas I (2002) Constraints on the Black Sea outflow to the Sea of Marmara during the last glacial-interglacial transition. Marine Geology 190(1–2): 19–34. https://doi.org/10.1016/S0025-3227(02)00340-7
  • Mangerud J, Jakobsson M, Alexanderson H, Astakhov V, Clarke GKC, Henriksen M, Hjort C, Krinner G, Lunkka JP, Moller P (2004) Ice-dammed lakes and rerouting of the drainage of northern Eurasia during the Last Glaciation. Quaternary Science Reviews 23(11–13): 1313–1332. https://doi.org/10.1016/j.quascirev.2003.12.009
  • Marino G, Rohling EJ, Sangiorgi F, Hayes A, Casford JL, Lotter AF, Kucera M, Brinkhuis H (2009) Early and middle Holocene in the Aegean Sea: Interplay between high and low latitude climate variability. Quaternary Science Reviews 28(27–28): 3246–3262. https://doi.org/10.1016/j.quascirev.2009.08.011
  • Marret F (1994) Distribution of dinoflagellate cysts in recent marine sediments from the East Equatorial Atlantic (Gulf of Guinea). Review of Palaeobotany and Palynology 84(1–2): 1–22. https://doi.org/10.1016/0034-6667(94)90038-8
  • Marret F, Leroy S, Chalie F, Gasse F (2004) New organic-walled dinoflagellate cysts from recent sediments of Central Asian seas. Review of Palaeobotany and Palynology 129(1–2): 1–20. https://doi.org/10.1016/j.revpalbo.2003.10.002
  • Marret F, Mudie P, Aksu A, Hiscott RN (2009) A Holocene dinocyst record of a two-step transformation of the Neoeuxinian brackish water lake into the Black Sea. Quaternary International 197(1–2): 72–86. https://doi.org/10.1016/j.quaint.2007.01.010
  • Matthiessen J (1995) Distribution patterns of dinoflagellate cysts and other organic-walled microfossils in recent Norwegian-Greenland Sea sediments. Marine Micropaleontology 24(3–4): 307–334. https://doi.org/10.1016/0377-8398(94)00016-G
  • Matthiessen J, Brenner W (1996) Chlorococcalalgen und dinoflagellaten Zysten in rezenten sedimenten des greifswalder boddens (südliche Ostesee). Senckenbergiana Maritima 27(1/2): 33–48.
  • Mertens KN, Bradley LR, Takano Y, Mudie PJ, Marret F, Aksu AE, Hiscott RN, Verleye TJ, Mousing EA, Smyrnova LL, Bagheri S, Mansor M, Pospelova V, Matsuoka K (2012) Quantitative estimation of Holocene surface salinity variation in the Black Sea using dinoflagellate cyst process length. Quaternary Science Reviews 39: 45–59. https://doi.org/10.1016/j.quascirev.2012.01.026
  • Moore P, Webb E (1978) An illustrated guide to pollen analysis. Hodder and Stoughton, London-Sydney-Auckland-Toronto, 192 pp.
  • Mudie PJ, Harland R (1996) Aquatic Quaternary. In: Jansonius J, McGregor DC (Eds) Palynology: Principles and Applications. AASP Foundation, Salt Lake City, 843–877.
  • Mudie P, Aksu A, Yaşar D (2001) Late Quaternary dinoflagellate cysts from the Black, Marmara and Aegean seas: Variations in assemblages, morphology and paleosalinity. Marine Micropaleontology 43(1–2): 155–178. https://doi.org/10.1016/S0377-8398(01)00006-8
  • Mudie P, Rochon A, Aksu A (2002) Pollen stratigraphy of Late Quaternary cores from Marmara Sea: Land-sea correlation and paleoclimatic history. Marine Geology 190(1–2): 233–260. https://doi.org/10.1016/S0025-3227(02)00349-3
  • Mudie P, Rochon A, Aksu A, Gillespie H (2004) Late glacial, Holocene and modern dinoflagellate cyst assemblages in the Aegean-Marmara-Black Sea corridor: Statistical analysis and re-interpretation of the early Holocene Noah’s Flood hypothesis. Review of Palaeobotany and Palynology 128(1–2): 143–167. https://doi.org/10.1016/S0034-6667(03)00117-9
  • Mudie P, Marret F, Aksu A, Hiscott R, Gillespie H (2007) Palynological evidence for climatic change, anthropogenic activity and outflow of Black Sea water during the Late Pleistocene and Holocene: Centennial- to decadal-scale records from the Black and Marmara Seas. Quaternary International 167–168: 73–90. https://doi.org/10.1016/j.quaint.2006.11.009
  • Mudie PJ, Leroy SAG, Marret F, Gerasimenko NP, Kholeif SEA, Sapelko T, Filipova-Marinova M (2011) Nonpollen palynomorphs: indicators of salinity and environmental change in the Caspian – Black Sea – Mediterranean corridor. In: Buynevich IV, Yanko-Hombach V, Gilbert AS, Martin RE (Eds) Geology and Geoarchaeology of the Black Sea region: Beyond the flood hypothesis. The Geological Society of America, Colorado, 89–115. https://doi.org/10.1130/2011.2473(07)
  • Müller U, Pross J, Tzedakis C, Gamble C, Kotthoff U, Schmiedl G, Wulf S, Christanis K (2011) The role of climate in the spread of modern humans into Europe. Quaternary Science Reviews 30(3–4): 273–279. https://doi.org/10.1016/j.quascirev.2010.11.016
  • Murray JW, Stewart K, Kassakian S, Krynytzky M, DiJulio D (2007) Oxic, suboxic, and anoxic conditions in the Black Sea. In: Yanko-Hombach V, Gilbert AS, Panin N, Dolukhanov PM (Eds) The Black Sea Flood Question: Changes in Coastline, Climate, and Human Settlement. Springer, Dordrecht, 1–26. https://doi.org/10.1007/978-1-4020-5302-3_1
  • Nikishin A, Korotaev M, Ershov M, Brunet M-F (2003) The Black Sea basin: Tectonic history and Neogene–Quaternary rapid subsidence modelling. Sedimentary Geology 156(1–4): 149–168. https://doi.org/10.1016/S0037-0738(02)00286-5
  • Oguz T, Besiktepe S (1999) Observations on the Rim Current structure, CIW formation and transport in the western Black Sea. Deep-Sea Research I. Deep-sea Research. Part I, Oceanographic Research Papers 46(10): 1733–1753. https://doi.org/10.1016/S0967-0637(99)00028-X
  • Oguz T, Latun VS, Latif MA, Vladimirov VV, Sur HI, Markov AA, Özsoy E, Kotovshchikov BB, Eremeev VV, Ünluata Ü (1993) Circulation in the surface and intermediate layers of the Black Sea. Deep-sea Research. Part I, Oceanographic Research Papers 40(8): 1597–1612. https://doi.org/10.1016/0967-0637(93)90018-X
  • Özsoy E, Latif MA, Beşiktepe Ş, Oğuz T, Güngör H, Ünlüata Ü, Ganines AF, Tuğrul S, Baştürk Ö, Yılmaz A, Yemenicioğlu S (1994) Monitoring Via Direct Measurements of the Model of Mixing and Transport of Wastewater Discharges into the Bosphorus Underflow (Hydrography. Sea-Level. Current and Flux Measurements in the Bosphorus Strait and Acoustical Chemical and Rhodamine-B Dye Tracer Studies of the Ahırkapı Waste Discharge). METU-IMS, Erdemli, Içel, Turkey.
  • Özsoy E, Latif M, Tugrul S, Ünluata U (1995) Exchanges with the Mediterranean, fluxes, and boundary mixing processes in the Black Sea. Bulletin de l'Institut Océanographique 15: 1–25.
  • Panin N (2008) The Black Sea: geology, environment and natural hazards. International Seminar on Natural Hazards in the Marine Area, Bucharest (Romania), 28–29 July 2008.
  • Panin N, Jipa DC, Gomoiu MT, Secrieru D (1999) Importance of Sedimentary Processes in Environmental Changes: Lower River Danube – Danube Delta – Western Black Sea System. In: Beşiktepe ST, Ünlüata Ü, Bologa AS (Eds) Environmental Degradation of the Black Sea: Challenges and Remedies. NATO Science Series, 56. Springer, Dordrecht, 23–41. https://doi.org/10.1007/978-94-011-4568-8_3
  • Pardoe H, Giesecke T, van der Knaap W, Svitavská-Svobodová H, Kvavadze E, Panajiotidis S, Gerasimidis A, Pidek I, Zimny M, Święta-Musznicka J, Latałowa M, Noryśkievicz A, Bozilova E, Tonkov S, Filipova-Marinova M, van der Leeuwen J, Kalniņa L (2010) Comparing pollen spectra from modified Tauber traps and moss samples: Examples from a selection of woodlands across Europe. Vegetation History and Archaeobotany 19(4): 271–283. https://doi.org/10.1007/s00334-010-0258-y
  • Peyron O, Guiot J, Cheddadi R, Tarasov P, Reille M, Beaulieu JL, Bottema S, Andrieu V (1998) Climatic reconstruction in Europe for 18,000 yr BP from pollen data. Quaternary Research 49(2): 183–196. https://doi.org/10.1006/qres.1997.1961
  • Popescu S-M, Dalesme F, Jouannic G, Escarguel G, Head M, Melinte-Dobrinescu M-C, Süto-Szentai M, Bakrac K, Clauzon G, Suc JP (2009) Galeacysta etrusca complex: Dinoflagellate cyst marker of Paratethyan influxes to the Mediterranean Sea before and after the peak of the Messinian Salinity Crisis. Palynology 33(2): 105–134. https://doi.org/10.2113/gspalynol.33.2.105
  • Prentice IC, Guiot J, Harisson S (1993) Mediterranean vegetation, lake levels and paleoclimate at the Last Glacial Maximum. Nature 360(6405): 658–670. https://doi.org/10.1038/360658a0
  • Prentice C, Guiot J, Huntley B, Jolly D, Cheddadi R (1996) Reconstructing biomes from palaeoecological data: A general method and its application to European pollen data at 0 and 6 ka. Climate Dynamics 12(3): 185–194. https://doi.org/10.1007/BF00211617
  • Pross J, Kotthoff U, Muller UC, Peyron O, Dormoy I, Schmiedl G, Kalaitzidis S, Smith AM (2009) Massive perturbation in terrestrial ecosystems of the Eastern Mediterranean region associated with the 8.2 kyr B.P. climatic event. Geology 37(10): 887–890. https://doi.org/10.1130/G25739A.1
  • Reille M (1992) Pollen et spores d’Europe et d’Afrique du nord. Laboratoire de botanique historique et palynologie, Marseille, 520 pp. [in French]
  • Reille M (1995) Pollen et spores d’Europe et d’Afrique du nord. Supplement 1. Laboratoire de botanique historique et palynologie, Marseille, 327 pp. [in French]
  • Reimer PJ, Bard E, Bayliss A, Beck JW, Blackwell PG, Bronk Ramsey C, Buck CE, Cheng H, Edwards RL, Friedrich M, Grootes PM, Guilderson TP, Haflidason H, Hajdas I, Hatt C, Heaton TJ, Hogg AG, Hughen KA, Kaiser KF, Kromer B, Manning SW, Niu M, Reimer RW, Richards DA, Scott EM, Southon JR, Turney CSM, van der Plicht J (2013) IntCal13 and MARINE13 radiocarbon age calibration curves 0-50000 years calBP. Radiocarbon 55(4): 1869–1887. https://doi.org/10.2458/azu_js_rc.55.16947
  • Roberts N, Eastwood WJ, Kuzucuoglu C, Fiorentino G, Caracuta V (2011) Climatic, vegetation and cultural change in the eastern Mediterranean during the mid-Holocene environmental transition. The Holocene 21(1): 147–162. https://doi.org/10.1177/0959683610386819
  • Roman S (1974) Palynoplanktologic analysis of some Black Sea cores. In: Degens E, Ross D (Eds) The Black Sea – Geology, Chemistry and Biology. American Association of Petroleum Geologists, Tulsa, Oklahoma, USA, 396–410.
  • Ross DA, Degens ET (1974) Recent sediments of Black Sea. In: Degens E, Ross D (Eds) The Black Sea – Geology, Chemistry and Biology. American Association of Petroleum Geologists, Tulsa, Oklahoma, USA, 183–199. https://doi.org/10.1306/M20377
  • Ryan W, Pitman W III, Major C, Shimkus K, Moskalenko V, Jones G, Dimitrov P, Görür N, Sakinç M, Yüce H (1997) An abrupt drowning of the Black Sea shelf. Marine Geology 138(1–2): 119–126. https://doi.org/10.1016/S0025-3227(97)00007-8
  • Sernander R (1908) On the evidence of post-glacial changes of climate furnished by the peat mosses of northern Europe. Geologiska Föreningens i Stockholm Förhandlingar 31: 423–448. https://doi.org/10.1080/11035890909443099
  • Shimkus K, Komarov A, Grakova I (1977) On the stratigraphy of deep-water Upper Quaternary sediments of the Black Sea. Oceanologia 7: 675–678. [in Russian]
  • Shopov V (1991) Biostratigraphy (on the basis of the Molluscan Fauna) of Quaternary Sediments from the Bulgarian Black Sea Shelf. Dr. of Sciences Thesis, Geological Institute of the Bulgarian Academy of Sciences, Sofia, Bulgaria. [in Bulgarian]
  • Shumilovskikh L, Tarasov PE, Arz HW, Fleitmann D, Marret F, Nowaczyk NR, Plessen B, Schlütz F, Behling H (2012) Vegetation and environmental dynamics in the southern Black Sea region since 18 kyr BP derived from the marine core 22-GC3. Palaeogeography, Palaeoclimatology, Palaeoecology 337–338: 177–193. https://doi.org/10.1016/j.palaeo.2012.04.015
  • Shumilovskikh L, Marret F, Fleitmann D, Arz HW, Nowaczyk N, Behling H (2013) Eemian and Holocene sea-surface conditions in the southern Black Sea: Organic-walled dinoflagellate cyst record from core 22-GC3. Marine Micropaleontology 101: 146–160. https://doi.org/10.1016/j.marmicro.2013.02.001
  • Siani G, Paterne M, Arnold M, Bard E, Métivier B, Tisnerat N, Bassinot F (2000) Radiocarbon reservoir ages in the Mediterranean Sea and Black Sea. Radiocarbon 42(2): 271–280. https://doi.org/10.1017/S0033822200059075
  • Soulet G, Menot G, Lericolais G, Bard E (2011) A revised calendar age for the last reconnection of the Black Sea to the global ocean. Quaternary Science Reviews 30(9–10): 1019–1026. https://doi.org/10.1016/j.quascirev.2011.03.001
  • Soulet G, Skinner LC, Beaupré SR, Galy V (2016) A Note on reporting of reservoir 14C disequilibria and age offsets. Radiocarbon 58(1): 205–211. https://doi.org/10.1017/RDC.2015.22
  • Stefanova I, Atanassova J, Delcheva M, Wright H (2006) Chronological framework for the Late Glacial pollen and macrofossil sequence in the Pirin Mountains, Bulgaria: Lake Besbog and LakeKremensko-5. The Holocene 16(6): 877–892. https://doi.org/10.1191/0959683606hol979rp
  • Stojanov N (1950) Geography of the vegetation. Nauka i Izkustvo Publishers, Sofia, 530 pp. [in Bulgarian]
  • Stuiver M, Grootes P, Brazunas T (1995) The GISP2 δ18O Climate record of the past 16500 years and the role of the sun, ocean, and volcanoes. Quaternary Research 44(3): 341–354. https://doi.org/10.1006/qres.1995.1079
  • Stuiver M, Reimer PJ, Reimer RW (2017) CALIB 7.1. http://calib.org [accessed 2017-10-29]
  • Tallantire PA (2002) The early-Holocene spread of hazel (Corylus avellana L.) in Europe north and west of the Alps: An ecological hypothesis. The Holocene 12(1): 81–96. https://doi.org/10.1191/0959683602hl523rr
  • Tantau I, Feurdean A, de Beaulieu JL, Reille M, Farcas S (2011) Holocene vegetation history in the upper forest belt on the Eastern Romanian Carpathians. Palaeogeography, Palaeoclimatology, Palaeoecology 309(3–4): 281–290. https://doi.org/10.1016/j.palaeo.2011.06.011
  • Tinner W, Lotter AF (2006) Holocene expansions of Fagus silvatica and Abies alba in Central Europe: Where are we after eight decades of debate? Quaternary Science Reviews 25(5–6): 526–549. https://doi.org/10.1016/j.quascirev.2005.03.017
  • Todorova H (1986) Eneolithic in Bulgaria: Vth millennia BC. Nauka i Izkustvo Publishers, Sofia, 97 pp. [in Bulgarian]
  • Tonkov S, Bozilova E (1995) Palaeoecological data about the end of the Holocene climatic optimum in Bulgaria – IV millennium BC. Annual of Sofia University. Biology (Basel) 87(2): 5–17.
  • Tonkov S, Panovska H, Possnert G, Bozilova E (2002) The Holocene vegetation history in the Northern Pirin Mountain, southwestern Bulgaria: Pollen analysis and radiocarbon dating of core from Lake Ribno Banderishko. The Holocene 12(2): 201–210. https://doi.org/10.1191/0959683602hl535rp
  • Traverse A (1974) Palynological investigation of two Black Sea cores. In: Degens E, Ross D (Eds) The Black Sea – Geology, Chemistry and Biology. American Association of Petroleum Geologists, Tulsa, Oklahoma, USA, 267–275.
  • Traverse A (1978a) Palynological analysis of DSDP LEG 42 B (1975) cores from the Black Sea. In: Ross D, Neprochnov P (Eds) Initial Reports of Deep-Sea Drilling Project, v. 42 (2). U.S. Government Printing Office, Washington, USA, 993–1015. https://doi.org/10.2973/dsdp.proc.42-2.144.1978
  • Traverse A (1978b) Supplementary palynological information for site reports for DSDP Leg 42 B. In: Ross D, Neprochnov P (Eds) Initial Reports of Deep-Sea Drilling Project, v. 42 (2). U.S. Government Printing Office, Washington, USA, 29–305.
  • Tzedakis C (1993) Long-term tree populations in northwest Greece through multiple Quaternary climatic cycles. Nature 364(6436): 437–440. https://doi.org/10.1038/364437a0
  • van der Hammen T, Wijmstra T, Zagwijn W (1971) The floral record of the Late Cenozoic of Europe. In: Turekian K (Ed.) The Late Cenozoic glacial ages. Yale University Press, New Haven, Connecticut, 391–424.
  • van der Meer M, Sangiorgi F, Baas M, Brinkhuis H, Damste J, Schouten S (2008) Molecular isotopic and dinoflagellate evidence for Late Holocene freshening of the Black Sea. Earth and Planetary Science Letters 267(3–4): 426–434. https://doi.org/10.1016/j.epsl.2007.12.001
  • van Geel B (2001) Non-pollen palynomorphs. In: Smol JP, Birks HGB, Last WM (Eds) Tracking environmental change using lake sediments. Volume 3: Terrestrial, algal and silicaceous indicators. Kluwer, Dordrecht, 99–119. https://doi.org/10.1007/0-306-47668-1_6
  • van Geel B, Berglund BE (2000) A causal link between a climatic deterioration around 850 cal BC and a subsequent rise in human population in NW-Europe? Terra Nostra 7: 126–130.
  • van Geel B, Bohncke SJP, Dee H (1981) A palaecological study of an upper Late Glacial and Holocene sequence from ‘De Borchert’, The Netherlands. Review of Palaeobotany and Palynology 31: 367–448. https://doi.org/10.1016/0034-6667(80)90035-4
  • Velev S (2002) Climatic regions. In: Kopralev I (Ed.) Geography of Bulgaria. ForKom Publishers, Sofia, 155–157. [in Bulgarian]
  • Verleye TJ, Mertens KN, Louwye S, Arz HW (2009) Holocene salinity changes in the southwestern Black Sea: A reconstruction based on dinoflagellate cysts. Palynology 33(1): 77–100. https://doi.org/10.2113/gspalynol.33.1.77
  • Wall D, Dale B (1974) Dinoflagellates in Late Quaternary deep-water sediments of Black Sea. In: Degens E, Ross D (Eds) The Black Sea – Geology, Chemistry and Biology. American Association of Petroleum Geologists, Tulsa, Oklahoma, USA, 364–380.
  • Wall D, Dale B, Harada K (1973) Description of new fossil dinoflagellates from the Late Quaternary of the Black Sea. Micropaleontology 19(1): 18–31. https://doi.org/10.2307/1484962
  • Wanner H, Beer J, Bütikofer J, Crowley TJ, Cubasch U, Flückiger J, Goosse H, Grosjean M, Joos S, Kaplan JO, Küttel M, Müller SA, Prentice IC, Solomina O, Stocker TF, Tarasov P, Wagner M, Widmann M (2008) Mid- to Late Holocene climate change: An overview. Quaternary Science Reviews 27(19–20): 1791–1828. https://doi.org/10.1016/j.quascirev.2008.06.013
  • Wick L, Lemcke G, Sturm M (2003) Evidence of Late Glacial and Holocene climatic change and human impact in eastern Anatolia: High-resolution pollen, charcoal, isotopic and geochemical records from the laminated sediments of Lake Van, Turkey. The Holocene 13(5): 665–675. https://doi.org/10.1191/0959683603hl653rp
  • Wright HE, Ammann B, Stefanova I, Atanassova J, Margalitadze N, Wick L, Blyakharchuk T (2003) Late-glacial and Early-holocene dry climates from the Balkan Peninsula to southern Siberia. In: Tonkov S (Ed.) Aspects of Palynology and Palaeoecology: Festschrift in honour of Elissaveta Bozilova. Pensoft Publishers, Sofia-Moscow, 127–136.
  • Yanko-Hombach VV (2007) Controversy over Noah’s Flood in the Black Sea: Geological and foraminiferal evidence from the shelf. In: Yanko-Hombach V, Gilbert A, Panin N, Dolukhanov P (Eds) The Black Sea Flood Question: Changes in Coastline, Climate, and Human Settlement, Springer, Dordrecht, 149–204. https://doi.org/10.1007/978-1-4020-5302-3_7
  • Yu S-Y, Berglund BE (2007) A dinoflagellate cyst record of Holocene climate and hydrological changes along the southeastern Swedish Baltic coast. Quaternary Research 67(2): 215–224. https://doi.org/10.1016/j.yqres.2006.12.004
login to comment