Jean-Pierre Suc

Speranta-Maria Popescu

Damien Do Couto

Georges Clauzon

Jean-Loup Rubino

Carmen Mihaela

Frédéric Melinte-Dobrinescu

Nikola Quillévéré

Ivan Dumurdžanov

Vesna Zagorchev

Lesić

Melinte-Dobrinescu

Frédéric Quillévéré

Jean-Pierre Brun

Nikola Dumurdžanov

Ivan Zagorchev

Vesna Lesić

Dragana Tomić

Dimitrios Sokoutis

Bertrand Meyer

Rodica Macaleţ

Helena Rifelj

Marine gateway vs. fluvial stream within the Balkans from 6 to 5 Ma

Keywords: Mediterranean -Paratethys connections, Messinian Salinity Crisis, Marine microfossils, Fluvial erosion, Balkans corridor

come

Introduction

The matter of Neogene relationships between the Mediterranean Sea and the former Paratethys, an adjacent brackish basin (Fig. 1A) has been early tackled by [START_REF] Rögl | Vom Zerfall der Tethys zu Mediterran und Paratethys. Die neogene Paläogeographie und Palinspastik des zirkum-mediterranen Raumes[END_REF] at a time where stratigraphic correlations were insufficiently ensured. They considered a gateway in the area of the İ stanbul area before and after the Messinian Salinity Crisis (MSC), a hypothesis then adopted by [START_REF] Popov | Late Miocene to Pliocene palaeogeography of the Paratethys and its relation to the Mediterranean[END_REF]. However, the assumption of a gateway or a drainage network through the Balkans was suggested by [START_REF] Hsü | Stratigraphy of the lacustrine sedimentarion in the Black Sea[END_REF] (Fig. 2A) and [START_REF] Hsü | History of the Mediterranean salinity crisis[END_REF]Hsü et al. ( , 1978b) ) (Fig. 2B), respectively. A possible corridor was sought along the Struma/Strymon palaeo-river, but [START_REF] Kojumdgieva | Evolution géodynamique du bassin Egéen pendant le Miocène supérieur et ses relations à la Paratéthys Orientale[END_REF] was inconclusive because of the obvious lack of marine or even lacustrine sediments in the Sandanski graben [START_REF] Ognjanova-Rumenova | Neogene Lake Basins in Southwestern Bulgaria -biostratigraphy, palaeoecology and diatom floristic interrelations[END_REF]. Such a hypothesis arose from the old concept of a North Aegean Lake [START_REF] Cvijić | L'ancien lac égéen[END_REF], revisited by [START_REF] Krstić | Pliocene ostracodes of the Paludinian Beds in Pannonian Plain, Serbian Part. The Natural History Museum and Society "Nauka I društvo[END_REF] who studied ostracod faunas. A recent review by [START_REF] Krstić | The Neogene Lakes on the Balkan Land[END_REF] adopted the idea about a gateway between the Pannonian ("Paludinian" Lake) and the Aegean basins via Morava and Vardar but has overestimated to some extent the territories of the "Aegean Lake" and postulated a not-proven Pliocene connection along the Palaeo-Struma. An important turnaround in the concept occurred when marine calcareous nannoplankton was discovered in the Dacic Basin [START_REF] Mǎrunţeanu | The connection between the Dacic and Mediterranean Basins based on calcareous nannoplankton assemblages[END_REF] within Lower and Upper Pontian deposits ascribed to the reverse C3r palaeomagnetic Chron (Snel et al., 2006a; see Figure 4 for correspondence between global stages and Paratethyan Stages). Calcareous nannoflora of the Lower Pontian sediments belongs to the NN11b Subzone (characterized by the occurrence of Discoaster quinqueramus), that of the Upper Pontian sediments belongs to NN12 Zone (characterized by the occurrence of Ceratolithus acutus) (Snel et al., 2006a). After almost similar results, [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF] concluded that high sea-level episodes connected the Mediterranean and Dacic Basin just prior the MSC and just after it, also documented by planktonic foraminifers [START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF]. [START_REF] Popescu | Pollen record and integrated high-resolution chronology of the Early Pliocene Dacic Basin (Southwestern Romania)[END_REF] showed that Mediterranean micro-organisms (calcareous nannofossils, planktonic foraminifers, dinoflagellate cysts) invaded Dacic Basin before 5.4 Ma, an age that has been further specified at 5.45 Ma by [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF]. This marine water input reached the Dacic Basin significantly before the Euxinian Basin (i.e., the Black Sea; Fig. 1A) where it has been recorded at 5.31 Ma [START_REF] Popescu | Pliocene and Lower Pleistocene vegetation and climate changes at the European scale: Long pollen records and climatostratigraphy[END_REF]. In addition to this time-lag, the absence of Mediterranean -Black Sea connection in the area of İ stanbul since ca. 8 Ma (Suc et al., submitted) also supports a marine gateway through the Balkans, as documented by [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF] (Fig. 2C). As a consequence, a new pattern of crossed exchanges of organisms between the Mediterranean and the Dacic Basin during high sea-level episodes can be tested to explain the episodes of Lago Mare 1 (before the peak of the MSC) and 3 (after the peak of the MSC) in the Mediterranean, as argued by [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF], [START_REF] Popescu | Galeacysta etrusca complex, dinoflagellate cyst marker of Paratethyan influxes into the Mediterranean Sea before and after the peak of the Messinian Salinity Crisis[END_REF], [START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF], [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF], [START_REF] Bakrač | Middle and Late Miocene palynological biozonation of the south-western part of Central Paratethys (Croatia)[END_REF] and [START_REF] Couto | Lago Mare and the Messinian Salinity Crisis: Evidences from the Alboran Sea (S. Spain)[END_REF].

This paper aims [START_REF] Gillet | La faune de mollusques du Pontien de trilophos (SW Thessaloniki)[END_REF] to complete the field data sets published in [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF], [START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF] and [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF] that supports the Balkans marine gateway, and (2) to discuss how this passage was affected by the peak of the MSC.

Methodological approach and geographic area of investigation

As soon as marine calcareous nannofossils were recorded in the Dacic Basin for a timeinterval encompassing the MSC [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF]Snel et al., 2006a), an approach has been developed for re-visiting the Balkans region that was a strong candidate to have hosted the Mediterranean -Dacic Basin exchange gateway (see above).

This investigation was done following the present-day back to back fluvial networks (Fig. 3): Struma/Strymon and Vardar/Axios towards the Mediterranean Basin, Morava, Timok and Iskar towards the Danube River; the two obvious options for a corridor being via Iskar -Sofia Lake -Struma/Strymon or Morava -Vardar/Axios; in view of the lack of marine or lacustrine Neogene sediments along most of Palaeo-Struma to the north of the Belasitsa horst (Zagorčev, 1992a;[START_REF] Zagorchev | Late Cenozoic development of the Strouma and Mesta fluviolacustrine systems, SW Bulgaria and northern Greece[END_REF][START_REF] Ognjanova-Rumenova | Neogene Lake Basins in Southwestern Bulgaria -biostratigraphy, palaeoecology and diatom floristic interrelations[END_REF], we concentrated [START_REF] Clauzon | The Pliocene Gilbert-type fan delta of Dračevo (Skopje area, Republic of Macedonia): Paleogeographic inference[END_REF][START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF] on the second option (Morava-Vardar) where a narrow divide must be noticed in the area of Skopje. In this frame, we searched for:

-calcareous nannofossils because of (1) their worldwide stratigraphic significance regularly reasserted [START_REF] Martini | Standard Tertiary and Quaternary calcareous nannoplankton zonation[END_REF][START_REF] Perch-Nielsen | Cenozoic calcareous nannofossils[END_REF][START_REF] Berggren | Geochronology, Time Scales and Global Stratigraphic Correlation: A Unified Temporal Framework for an Historical Geology[END_REF][START_REF] Raffi | A review of calcareous nannofossil astrobiochronology encompassing the past 25 million years[END_REF], and (2) their ability to be the first invasive indicators of previously isolated basins by marine surface waters where they are living [START_REF] Winter | Biogeography of living coccolithophores in ocean waters[END_REF], and -pollen grains and dinoflagellate cysts from some localities in order [START_REF] Gillet | La faune de mollusques du Pontien de trilophos (SW Thessaloniki)[END_REF] to complete the available dataset [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF][START_REF] Popescu | Pollen record and integrated high-resolution chronology of the Early Pliocene Dacic Basin (Southwestern Romania)[END_REF][START_REF] Jiménez-Moreno | Neogene flora, vegetation and climate dynamics in southeastern Europe and the northeastern Mediterranean. Micropalaeontol[END_REF] and (2) to contribute to reconstruct palaeoenvironments.

Data from the Aegean to the Dacic Basin

Both basins and their drainage areas have been subject of detailed and important studies and numerous publications that will be only partially referred to in the present paper. We restrain mostly to new evidence of utmost importance for the problem of the Balkans corridor. Available plus unpublished data are examined from the South to the North, with a focus on the most significant precisely dated localities, and following the different branches of the socalled North Aegean Lake (Fig. 2).

3.1.North Aegean area

The remarkable locality of Akropotamos is located near the Aegean shoreline (locality 4 in Fig. 1B). It is known for a long time because of the occurrence of Messinian gypsums (with NN11 calcareous nannofossils from the intercalated clays; Figs. 5B-C) in a quarry topped by (1) limestones (40° 47' 30.62" N, 24° 01' 56.46" E;[START_REF] Steffens | Field guide to the neogene of Northern Greece (Thessaloniki area and Strimon Basin)[END_REF][START_REF] Dertmizakis | Ecostratigraphic observations on the Messinian deposits of Akropotamos area (Kavala, N. Greece)[END_REF][START_REF] Dertmizakis | Ecostratigraphic observations on the Messinian deposits of Akropotamos area (Kavala, N. Greece)[END_REF]) and (2) a clastic succession, the silty beds of which being ascribed to the calcareous nannoplankton NN12 Zone [START_REF] Steffens | Field guide to the neogene of Northern Greece (Thessaloniki area and Strimon Basin)[END_REF]. The nannoplankton ages have been specified by Snel et al. (2006b) who recorded Ceratolithus acutus and Triquetrorhabdulus rugosus in the silts of the clastic member hence ascribed to the earliest NN12b Subzone (Fig. 4). The issue raised by this section, which does not concern the age as it is confirmed by A. Di Stefano from our sampling (40° 46' 01.63" N, 24° 01' 39.99" E; Fig. 5H), concerns the interpretation of the stratigraphic relationships between the sedimentary units. In Snel et al. (2006b: fig. 3), the section is presented in a vertical log as usually done for boreholes. Indeed, such a drawing covers up the lateral relationships between the units, aiming to instill a continuous transition from Messinian to Zanclean and, as a consequence, to disregard the impact of the MSC. Actually, we observed that the clastic unit is nested of about 70 m within the evaporitic unit (including the topmost limestones) (Fig. 5B). These units are separated by an erosional surface which can be followed from place to place (e.g.: 40° 46' 36.20" N, 24° 01' 50.19" E; 40° 45' 54.57" N, 24° 00' 34.42" E; Figs. 5D-E) and, according to the above-mentioned ages, is the Messinian Erosional Surface (MES), contemporaneous of the massive drawdown of the Mediterranean Sea level between 5.60 and 5.46 Ma according to [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF]. The overlying clastic sediments are those of a Gilbert-type fan delta with all its characteristic constituents (submarine foreset beds and bottomset beds, and subaerial topset beds: Figs. 5B, D-G), attesting a fast marine reflooding (for a detailed description, see: [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF].

In the same area, the MES has also been observed near the statue of the Amphipolis Lion (40° 48' 08.89" N; 23° 50' 30.89" E; locality 3 in Fig. 1B) where gravely to sandy foreset beds truncate Miocene continental deposits (Figs. 6A-B). Because of their dipping orientation, we consider that the foreset beds belong to a tributary of a large Gilbert-type fan delta, the bottomset beds of which are known for their marine Zanclean malacofauna (Syrides, 1998a). The topset beds, exposed along the road going up the Strymon River, provided a pollen flora typical of Zanclean coastal environments in the region (Table 1; [START_REF] Jiménez-Moreno | Neogene flora, vegetation and climate dynamics in southeastern Europe and the northeastern Mediterranean. Micropalaeontol[END_REF]. This Gilbert-type fan delta is clearly nested within the Miocene deposits that, in several places, provided Paratethyan molluscs (dreisseinids and lymnocardiids) (Syrides, 1998b).

The Prinos hydrocarbon Field (Fig. 1B) lies within a pull-apart basin beneath a thick evaporitic unit at the top of which an erosional unconformity has been revealed (Fig. 7; [START_REF] Proedrou | Prinos Field-Greece, Aegean Basin[END_REF][START_REF] Proedrou | Prinos Basin -A model for oil exploration[END_REF]. Thick clastic deposits were deposited below and above the evaporitic unit that also contains several thin clastic intercalations. Biostratigraphic data available from Nestos 1 and 2 wells (location 5 in Fig. 1B), located on its northeast edge, indicate (1) the occurrence of the calcareous nannofossil Discoaster quinqueramus between 1030 and 1050 m bsf (below sea floor) in Nestos 2, and (2) the occurrence of the planktic foraminifer Globorotalia puncticulata between 585 and 600 m bsf in well Nestos 1 and between 800 and 870 m bsf in well Nestos 2 (Fig. 7; [START_REF] Sidiropoulos | Geologie der Griechischen Tertiärbecken Nestos-Prinos und Xanti-Komotini[END_REF]. In well Nestos 2, the benthic foraminifer Uvigerina rutila, characteristic of the Lower Pliocene [START_REF] Iaccarino | Mediterranean Miocene and Pliocene planktic foraminifera[END_REF], has been continuously recorded from 480 to 830 m bsf, supporting that the sediments immediately overlying the evaporitic unit belong to Zanclean. The underlying evaporitic unit is Messinian in age, as attested by the occurrence in-between of Discoaster quinqueramus (Fig. 7), a calcareous nannofossil species that disappeared at 5.54 Ma (Fig. 4; [START_REF] Raffi | A review of calcareous nannofossil astrobiochronology encompassing the past 25 million years[END_REF]. In the Aegean Sea, many seismic profiles document its status of peripheral shallow water basin during the 1 st step of the MSC [START_REF] Clauzon | Alternate interpretation of the Messinian salinity crisis: Controversy resolved?[END_REF][START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF] because of the identification of the MES [START_REF] Schuster | CIDOG Mer Egée/Mer de Crète[END_REF]Hsü et al., 1978a;[START_REF] Anastasakis | Upper Cenozoic stratigraphy and paleogeographic evolution of Myrtoon and adjacent basins, Aegean Sea, Greece[END_REF][START_REF] Crombez | Evolution cénozoïque des basins sédimentaires de la mer Egée (Grèce)[END_REF], well-dated at Site DSDP 378 northward of Crete where it cuts Messinian peripheral evaporites (Hsü et al., 1978a). Accordingly, the erosional unconformity revealed by [START_REF] Proedrou | Prinos Field-Greece, Aegean Basin[END_REF] at the top of the evaporitic unit can be interpreted as the MES (Fig. 7), which was placed at 965 m in the Nestos 2 section in agreement with biostratigraphic data and pollen content [START_REF] Jiménez-Moreno | Neogene flora, vegetation and climate dynamics in southeastern Europe and the northeastern Mediterranean. Micropalaeontol[END_REF].

Near Thessaloniki, at Trilophos (locality 2 in Fig. 1B), we have observed a gently erosional contact separating the grey sands, limestones and clays with Paratethyan molluscs of the Trilophos Formation (Syrides, 1998b) from the underlying reddish continental clays and sands, Late Miocene in age according to mammals (Steffens et al., 1979) (Fig. 6C: 40° 27' 41.73" N, 22° 59' 34.94" E). In its uppermost part, the Trilophos Formation (ca. 20 m in thickness) is marine according to abundant Mactra shells among other marine species [START_REF] Gillet | La faune de mollusques du Pontien de trilophos (SW Thessaloniki)[END_REF]. Two samples from these beds (Fig. 6D; 40° 27' 33.65" N, 22° 59' 10.96" E) provided a diversified calcareous nannoflora with, among other species, Triquetrorhabdulus rugosus and Ceratolithus acutus, which together characterize the early NN12b Subzone (Fig. 4). The gently underlying erosional surface is thus the MES (Fig. 6C). The thin Trilophos Formaton is interpreted as indicating a lagoonal environment deposition affected by a marine incursion. Accordingly, the Trilophos section is a noticeable example of co-occurrence of Paratethyan and marine Mediterranean organisms, as in the Intepe section on the southern shoreline of the Dardanelles Strait [START_REF] Melinte-Dobrinescu | The Messinian Salinity Crisis in the Dardanelles region: Chronostratigraphic constraints[END_REF]. The coastal environment is also supported by a pollen record, with abundant halophytes (Amaranthaceae-Chenopodiaceae) and freshwater plants (Sparganium, Typha), which also points out a dense forest cover back to the delta up to high altitude of the Olympus Mount (Table 1; [START_REF] Biltekin | Vegetation and climate of North Anatolian and North Aegean region since 7 Ma according to pollen analysis[END_REF].

In the Prosilio section (locality 1 in Fig. 1B; 40° 07' 57.28" N, 21° 56' 01.30" E), an erosional surface separates the Lava Member from the overlying Prosilio Member made of sands and conglomerates (Fig. 6E-G; [START_REF] Steenbrink | Late Miocene to Early Pliocene depositional history of the intramontane Florina-Ptolemais-Servia Basin, NW Greece: Interplay between orbital forcing and tectonics[END_REF]. Regional stratigraphic correlations and datations allowed [START_REF] Steenbrink | Late Miocene to Early Pliocene depositional history of the intramontane Florina-Ptolemais-Servia Basin, NW Greece: Interplay between orbital forcing and tectonics[END_REF] to propose a magnetostratigraphy of the section (reminded in Figure 6E). The erosional truncation impacted sediments of the normal Chron C3An.1n. This leads us to believe that this erosional event, even of a relatively mild intensity, is not local but probably corresponds to the erosional episode caused by the peak of the MSC (Figs. 4). The overlying Tomea Eksi Member shows near its top a palaeomagnetic reversal from Chron C3r to the normal Chron C3n.4n that corresponds to the earliest Zanclean (Streenbrink et al., 2006). Marine fish remains from the Mugilidae family have been recorded in the middle part of the Prosilio Member that indicates a marine incursion in the area (M. Böhme and A. Ilg, in litteris). Karistineos andGeorgiades-Dikeoulia (1985-1986) established the presence of marine Pliocene in the Serres Basin (location 6 in Fig. 1B) on the basis of some molluscs such as Pecten benedictus and Ostrea lamellose. They described a prograding deltaic context including some brackish layers with Paratethyan molluscs (i.e., the so-called "Choumnikon beds"; see also : Syrides, 1998b). Unfortunately, relationships between the constituents of the Pliocene deltaic system remained unclear. Road works allowed us to observe in 2009 newly exposed sections showing the obvious nesting of the Pliocene deltaic coarse clastic deposits within the Miocene piedmont of the Vrondou Mount (Figs. 8A-B). There (41° 08' 58.94" N, 23° 32' 13.48" E), near Kato Metochi (Fig. 8A), foreset beds of a sigmoid sedimentary construction, although broken by small transverse faults indicated by [START_REF] Psilovikos | A depositional sedimentary model for the Neogene uraniferous lignites of the Serres graben, Greece. Palaeogeogr[END_REF], show that they belong to a Gilbert-type fan delta (Fig. 8C; [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF]. These foreset beds were displayed by Karistineos andGeorgiades-Dikeoulia (1985-1986) in photographs of their figures 5 and 6 despite a description as cross bedding sediments.

The erosional contact of the foreset beds on the Miocene piedmont deposits was exposed at the same place (Fig. 8D) and refers to the MES (Fig. 8D). Bottom set beds of the Kato Metochi Gilbert-type fan delta have been sampled in a nearby place (41° 08' 52.17" N, 23° 32' 14.46" E; Fig. 8E) where lignites interfinger to clays and silts. They did not provide any calcareous nannofossil, probably because of their upstream location within the Gilbert-type fan delta. Karistineos andGeorgiades-Dikeoulia (1985-1986) pointed out the occurrence of a debris-flow (with granite blocks) at the base of the Pliocene marine sediments. They correspond to the commonly observed coarse deposit marking around the Mediterranean the post-MSC marine reflooding (see: [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF], as we noticed also in the nearby Sidirókastro locality where such debris-flows are covered by large olistoliths made of marble, that recalls the situation in the Roussillon Basin (Clauzon et al., submitted). A pollen analysis was realized on samples from the area of Kato Metochi and ascribed to the Late Miocene according to a poor and inconclusive micromammal fauna [START_REF] Karistineos | Palaeoenvironmental and palaeoclimatic evolution of the Serres Basin (N. Greece) during the Miocene. Palaeogeogr., Palaeoclimatol[END_REF]. As the stratigraphic location of these pollen records is unclear, we performed the palynological analysis of a clayey sample from the above mentioned bottomset beds of the Kato Metochi Gilbert-type fan delta: it indicates a coastal swamp with Taxodiaceae in a regional forested context up to elevated altitudes (Table 1).

According to these new results and the synthesis of available data, we can conclude that the MES is present in the North Aegean area, less marked in the Vardar/Axios River drainage basin than in the Struma/Strymon River one, where thick and large Gilbert-type fan deltas developed as a result of the sudden marine reflooding which ended the MSC. However, several localities contain in their post-MSC deposits both Paratethyan molluscs and marine Mediterranean organisms (molluscs and/or calcareous nannofossils).

The intra-Balkans area

We re-visited the Pliocene of the Skopje region (locality 8 in Fig. 1B; Dumurdanov et al., 2005) where we evidenced and mapped near Dračevo a Gilbert-type fan delta nested within the basement to the South and Upper Miocene sediments to the North (Fig. 9A; [START_REF] Clauzon | The Pliocene Gilbert-type fan delta of Dračevo (Skopje area, Republic of Macedonia): Paleogeographic inference[END_REF]. Since this finding, after the record of calcareous nannofossils with a large biostratigraphic range in the bottomset beds of the Dračevo Gilbert-type fan delta, new researches for microfossils resulted in the record of Ceratolithus acutus among other species near Batinci (41° 55' 09.40" N, 21° 27' 39.10" E), which confirms the post-MSC age of the sediments (Fig. 9D). The three constituents of the sedimentary system (for more details, see: [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF], i.e. the foreset beds (e.g., Fig. 9C), bottomset beds (e.g., Fig. 9D) and topset beds (e.g., Fig. 9E) are exposed in several places. There, we observed two successive marine-continental transitions (vertical passage from foreset beds to topset beds) (Fig. 9B). Such a peculiarity is known in the Roussillon Basin (S. France) and results from a sea-level rise after the marine reflooding of the Mediterranean Basin (Clauzon et al. submitted). In the northern part of the Gilbert-type fan delta, the topset beds include reworked travertine blocks (Fig. 9E). In several localities, the travertine tops the Upper Miocene sediments and is very close to the erosional contact with the Pliocene deposits (i.e. the MES) [START_REF] Clauzon | The Pliocene Gilbert-type fan delta of Dračevo (Skopje area, Republic of Macedonia): Paleogeographic inference[END_REF]. Accordingly, this travertine may be considered as the spring deposit at the time of Messinian erosion. The Dračevo Gilbert-type fan delta was developed in the Messinian valley of the Markova River, a southern tributary of the Vardar River. Because of the weak Quaternary fluvial erosion in this short lateral valley, the Dračevo Gilbert-type fan delta constitutes the larger remaining evidence of the MSC in the middle part of the Balkans.

Between Serres and Skopje, only one area is designated with almost continuous deposits ascribed to Pliocene between Štip and Strumica (locality 7 in Figure 1B; [START_REF] Dumurdzanov | Evolution of the Neogene-Pleistocene Basins of Macedonia[END_REF][START_REF] Dumurdzanov | Cenozoic tectonics of Macedonia and its relation to the South Balkan extensional regime[END_REF]. There, near the Lakavica hamlet (Fig. 10A), one may observe a succession of SE prograding (probably subaqueous) clayey, sandy and conglomeratic deposits (41° 38' 53.16" N, 22° 13' 09.62" E; Figs. 10B-C), then of aggrading (probably continental) sandy, loamy and conglomeratic deposits (41° 38' 40.64" N, 22° 14' 41.55" E; Fig. 10E). The transition between the prograding and aggrading deposits is well exposed (41° 38' 38.82" N, 22° 13' 29.21" E; Fig. 10D). Despite several attempts, we did not find calcareous nannofossils either palynomorphs (pollen grains or dinoflagellate cysts) in the clays of the prograding body. [START_REF] Dumurdzanov | Cenozoic tectonics of Macedonia and its relation to the South Balkan extensional regime[END_REF] showed that these deposits unconformably overlie the Late Miocene ones dated by diatom floras. As a consequence, in spite of no direct dating, we consider the Lakavica deposits as a reliable candidate of Pliocene marine connection between the Skopje and Serres regions.

In the surrounding area of Niš (locality 9 in Fig. 1B), the Gabrovačka locality (Fig. 11A) shows exposed clays and conglomerates. This series, which unconformably lies over schistose basement and Lower-Middle Miocene deposits, is considered Lower Pliocene in age according to floral elements [START_REF] Pavlovic | Moravian freshwater basins and Intracarpathian basins (Niš Basin)[END_REF][START_REF] Kräutner | Geological map of the Carpatho -Balkanides between Mehadia, Oraviţa, Niš and Sofia[END_REF]. We sampled these deposits in two localities: the Gabrovačka Reka clays (43° 18' 09.10" N, 21° 55' 25.80" E; Fig. 11B) provided calcareous nannoplankton with, among several species, Triquetrorhabdulus rugosus and Ceratolithus acutus, attesting ascription to the NN12b Subzone [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF]Fig. 4); clays from the Gabrovačka stadium were barren (43° 18' 00.50" N, 21° 55' 12.60" E; Fig. 11C). It has been concluded that these clays and the overlying fluvial conglomerate (Fig. 11D) belong to a small deltaic system adjacent to the Nišava River (Fig. 11A), itself tributary of the Morava River [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF].

These data show that deltaic deposits of post-MSC age with marine incursions (NN12b calcareous nannoplankton Subzone; Fig. 4) are nested within older rocks (basement or Miocene) in the central Balkans along the Vardar and Morava drainage networks. [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF] and [START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF] evidenced post-MSC Gilbert-type fan deltas in the area of Turnu Severin (localities 11 and 12 in Fig. 1B). Such thick deposits are separated from older sediments or the basement by a strong erosional surface identified as the MES [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF][START_REF] Leever | Messinian sea level fall in the Dacic Basin (Eastern Paratethys): palaeogeographical implications from seismic sequence stratigraphy[END_REF][START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF]. This convergent system of two Gilbert-type fan deltas is a perfect example of a sedimentary succession in subaqueous conditions of selected terrigenous materials, from coarse conglomerates to sands then clays, deposited with a sedimentary dip within a relatively deep available space resulting from the Messinian fluvial erosion (Fig. 12A-D). Its ascription to the calcareous nannoplankton NN12b Subzone is unequivocally established in several localities where bottomset beds are narrowly interlocked with sandy to gravelly foreset beds (Fig. 12E; [START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF]. The "traditional upper Sarmatian age" of the Gura Vaii conglomerates (i.e., our foreset beds) cannot be supported [START_REF] Jipa | Sediment routing in a semi-enclosed epicontinental sea: Dacian Basin, Paratethys domain, late Neogene, Romania[END_REF][START_REF] Stoica | Paleoenvironmental evolution of the East Carpathian foredeep during the late Mioceneearly Pliocene (Dacian Basin; Romania)[END_REF] without careful fieldwork along the MES because there is some risk to inadvertently cross it. An upper Badenian -early Sarmatian (i.e. Serravallian) pollen flora is available in the close area of Gura Vaii, at Valea Morilor [START_REF] Jiménez-Moreno | Utilización del análisis polínico para la reconstrucción de la vegetación, clima y estimación de paleoaltitudes a lo[END_REF][START_REF] Jiménez-Moreno | Neogene flora, vegetation and climate dynamics in southeastern Europe and the northeastern Mediterranean. Micropalaeontol[END_REF] and can be compared to the Hinova one (Figs. 1B,12A), ascribed to the latest Messinianearliest Zanclean according to calcareous nannoplankton, magnetostratigraphy and cyclostratigraphy [START_REF] Popescu | Pollen record and integrated high-resolution chronology of the Early Pliocene Dacic Basin (Southwestern Romania)[END_REF]. The Valea Morilor pollen flora includes abundant pollen grains of Engelhardia with many other tropical and subtropical plants (Euphorbiaceae, Mussaenda type, Rubiaceae, Theaceae, Arecaceae, Sapotaceae, Distylium, Rhodoleia, Platycarya, Rhoiptelea). These thermophilous elements are not present at Hinova where Engelhardia is poorly represented. There are the usual differences between upper Badenianearly Sarmatian (i.e., Serravallian) and early Zanclean pollen floras in the region [START_REF] Jiménez-Moreno | Neogene flora, vegetation and climate dynamics in southeastern Europe and the northeastern Mediterranean. Micropalaeontol[END_REF].

The Dacic Basin

At last, the clarification done by [START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF] concerning the misapprehension of the biostratigraphic and stratigraphic data from Hinova by [START_REF] Krijgsman | Rise and fall of the Paratethys Sea during the Messinian Salinity Crisis[END_REF], was ignored by [START_REF] Stoica | Paleoenvironmental evolution of the East Carpathian foredeep during the late Mioceneearly Pliocene (Dacian Basin; Romania)[END_REF]. The pre-and post-MSC marine incursions by Mediterranean waters in the Dacic Basin are validated by the data published by [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF], Snel et al. (2006a), and[START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF] including planktonic foraminifers.

DISCUSSION

We successively discuss in which context the proposed gateway through the Balkans may have occurred, which way did it follow and how did it probably evolve during the peak of the MSC. We also search for the possible relationships with the primary proposed hypotheses in the area (Figs. 2A-B;[START_REF] Hsü | Stratigraphy of the lacustrine sedimentarion in the Black Sea[END_REF][START_REF] Hsü | History of the Mediterranean salinity crisis[END_REF]Hsü et al., , 1978b)). However, we start with a clarification concerning reliability of calcareous nannofossil biostratigraphy.

4.1.Reliability of Ceratolithus acutus as marker of the Mediterranean reflooding

Some post-MSC sediments have been dated by the occurrence, among others, of the calcareous nannofossil Ceratolithus acutus. This species, easily recognizable, has a very short duration of life (5.345 -5. 04 Ma; Fig. 4; [START_REF] Raffi | A review of calcareous nannofossil astrobiochronology encompassing the past 25 million years[END_REF], which makes it an excellent biostratigraphic marker of the post-MSC marine reflooding of the Mediterranean Basin. However, it is indicated as rare in the Mediterranean [START_REF] Castradori | Calcareous nannofossils in the basal Zanclean of the Eastern Mediterranean Sea: remarks on paleoceanography and sapropel formation[END_REF][START_REF] Di Stefano | Refinements of calcareous nannofossil biostratigraphy at the Miocene/Pliocene Boundary in the Mediterranean region[END_REF]. Its scarcity encouraged [START_REF] Grothe | Black Sea desiccation during the Messinian Salinity Crisis: Fact or fiction?[END_REF] to discredit its significance as a global biostratigraphic fine marker, although it has been recorded by several specialists in many localities around the Mediterranean and peripheral Paratethyan basins (e.g.: [START_REF] Braga | Testing models for the Messinian salinity crisis: The Messinian record in Almería[END_REF]Snel et al., 2006aSnel et al., , 2006b;;[START_REF] Soria | The Messinian-early Pliocene stratigraphic record in the southern Bajo Sedura Basin (Betic Cordillera, Spain): Implications for the Mediterranean salinity crisis[END_REF][START_REF] Melinte-Dobrinescu | The Messinian Salinity Crisis in the Dardanelles region: Chronostratigraphic constraints[END_REF][START_REF] Popescu | Pliocene and Lower Pleistocene vegetation and climate changes at the European scale: Long pollen records and climatostratigraphy[END_REF][START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF][START_REF] Cornée | The Early Pliocene reflooding in the Western Mediterranean: New insights from the rias of the Internal Rif, Morocco[END_REF]. According to their stratigraphic position, all these records immediately followed the post-MSC first inflow of surface marine waters in a Mediterranean peripheral basin or in a Paratethyan sub-basin, which occurred at 5.46 Ma, i.e. prior to the Zanclean Global boundary Stratotype-Section and Point (for detail, see: Bache et al., 2012, and references therein). This neo-species may have benefited from (1) the rush of surface waters, and (2) the missing biological competition in a newly invaded space of life. In this way, the contradiction could be solved between relatively frequent occurrences of Ceratolithus acutus in shallow peripheral deposits and its scarcity in deep-water sediments.

Context of a gateway through the Balkans

Occurrence of Paratethyan organisms in the northern Aegean and evidence of marine Mediterranean species in the Dacic Basin before and after the MSC (prior to their arrival in the Black Sea) led us to explore first the most direct way via the Sofia Basin [START_REF] Kamenov | Stratigraphy of the Neogene in the Sofia basin[END_REF]) and the Blagoevgrad, Simitli and Sandanski grabens along the Struma River (Fig. 1B). This prospect was encouraged by the occurrence of few marine dinoflagellate cysts in the earliest Dacian of the Ravno Pole borehole drilled in the eastern part of the Sofia Basin [START_REF] Drivaliari | Aspects climatiques et paléogéographiques d'un transect latitudinal (de la Roumanie au Delta du Nil)[END_REF][START_REF] Jiménez-Moreno | Neogene flora, vegetation and climate dynamics in southeastern Europe and the northeastern Mediterranean. Micropalaeontol[END_REF]. However, diatoms from the C14 Katina borehole drilled near Sofia do not document any marine influence during the Pontian and early Dacian [START_REF] Ognjanova-Rumenova | Paleoecological development of the Sofia Neogene Basin (Southwestern Bulgaria) based on sedimentological and paleontological evidence[END_REF]. The record of scarce dinoflagellate cysts at Ravno Pole must thus be considered as a reworking. Field observations in the grabens along the Struma River established that no marine influx reached the area [START_REF] Kojumdgieva | Evolution géodynamique du bassin Egéen pendant le Miocène supérieur et ses relations à la Paratéthys Orientale[END_REF]Zagorčev, 1992a;[START_REF] Spassov | Late Neogene stratigraphy, biochronology, faunal diversity and environments of South-west Bulgaria (Struma river valley)[END_REF][START_REF] Zagorchev | Late Cenozoic development of the Strouma and Mesta fluviolacustrine systems, SW Bulgaria and northern Greece[END_REF], as also stated by diatoms [START_REF] Ognjanova-Rumenova | Some aspects and problems concerning diatom biochronology for the Neogene in the region of the Balkan Peninsula[END_REF][START_REF] Ognjanova-Rumenova | Neogene Lake Basins in Southwestern Bulgaria -biostratigraphy, palaeoecology and diatom floristic interrelations[END_REF]. Thus, no passage may be supported along the Struma River system (Fig. 3). Similar field investigations were realized from Prilep to Ptolemais (Pelagonia Basin) and in the Ohrid Lake that came to the same conclusion, also supported by the diatom records [START_REF] Ognjanova-Rumenova | Upper Neogene siliceous microfossils from Pelagonia Basin (Balkan Peninsula)[END_REF], 2006).

Considering the above-described localities with post-MSC marine sediments, we propose a course for the gateway through the Balkans following as much as possible the present-day major fluvial systems shown in Figure 3 (Strymon, Vardar, Morava, Timok; Fig. 13). Three areas are nowadays relatively elevated on the corresponding altimetric profile (Fig. 14): they are the divides delimiting their respective fluvial networks, respectively at 443 m (A), 437 m (B) and 602 m (C) in altitude. So few moderately high points along a profile uplifted of 200 -300 m on the whole support the old age of the fluvial networks as inheritance of the Aegean Lake system, as suggested by [START_REF] Zagorchev | Late Cenozoic development of the Strouma and Mesta fluviolacustrine systems, SW Bulgaria and northern Greece[END_REF] for the southern part of the Balkans.

The present-day altitude envelope of the marine-continental transition (MCT) (i.e., the Zanclean seashore) in the six deltaic systems observed along the profile shows an arching that is homothetic to the altimetric profile of the major fluvial systems (Fig. 14). This indicates that the Zanclean seashore has been deformed since the deposition of the deltaic systems. Along most of the profiles, from Turnu Severin and Lakavica, the altitude difference between the MCTs and the closest river profiles is lower than 130 m. Conversely, close to the Strymon Basin the altitude difference reaches values as high as 240 m (Serres) and 250 m (Akropotamos). As the whole region underwent extension since the Pliocene (Zagorčev, 1992a(Zagorčev, , 1992b;;[START_REF] Burchfiel | Cenozoic Extension in Bulgaria and Northern Greece: the Northern Part of the Aegean Extensional Regime[END_REF][START_REF] Westaway | Late Cenozoic extension in SW Bulgaria: a synthesis[END_REF][START_REF] Brun | Kinematics of the Southern Rhodope Core Complex (North Greece)[END_REF], it cannot be considered that the highest MCT altitudes (Lakavica and Skopje) remained unchanged since the Zanclean. However, taking these highest MCT altitudes as a reference, it is especially relevant to discuss the tectonic meaning of the lowest MCT altitudes toward the Pannonian Basin and the Aegean. Toward the North, the 130 m altitude difference between the Turnu Severin MCT and the Danube River could be entirely due to fluvial incision due to the ongoing subsidence of the Pannonian Basin [START_REF] Csato | A three-dimensional stratigraphic model for the Messinian crisis in the Pannonian Basin, eastern Hungary)[END_REF] whereas the 100-160 m altitude difference between Lakavica-Skopje and Nis-Turnu Severin is more likely due to moderate extension [START_REF] Dumurdzanov | Cenozoic tectonics of Macedonia and its relation to the South Balkan extensional regime[END_REF]. Toward the Southeast, the altitude difference between the Serres-Akropotamos MCT and Strymon River that reaches up to 280 m is related to post-Zanclean extension in the onshore and offshore basins of Strymon and Orfanos, respectively. These basins located along the boundary between the Southern Rhodope Core Complex and the Chalkidiki Block occurred in clockwise rotation [START_REF] Brun | Kinematics of the Southern Rhodope Core Complex (North Greece)[END_REF], as documented by palaeomagnetism [START_REF] Dimitriadis | Dextral rotation and tectonomagmatic evolution of the southern Rhodope and adjacent regions (Greece)[END_REF]. The Chalkidiki Block, which underwent a bulk rotation of 30° since the Eocene, still rotated by 11°-13° since the Pliocene coeval with a 10 to 15 Km of extension in the Strymon Basin. Whereas it is difficult to quantify the amount of extension in the Orfanos Basin, it must be noted that this basin displays an extremely thick post-Messinian depocenter up to 3,000 m (see Fig. 15 in [START_REF] Brun | Kinematics of the Southern Rhodope Core Complex (North Greece)[END_REF].

In summary, the gateway through the Balkans has been controlled by the extensional tectonics that affected the region from the early Late Miocene to the Late Pliocene, i.e. from ca. 10 to 4 Ma, and created new basins on top of the recently exhumed metamorphic and plutonic rocks. The evolution of this Balkans corridor and, in particular, its forthcoming closure, which created obstacles for the free exchange of marine waters between the Paratethys realm and the Mediterranean, apparently resulted since the Zanclean from tectonically controlled subsidence to the North, toward the Pannonian Basin, and to the Southeast, toward the North Aegean.

Is it possible that two opposed marine branches have independently existed, a northern one linked to the Dacic Basin, a southern one linked to the Aegean Sea? In such a palaeogeographic context, the marine waters should have penetrated into the Dacic Basin via the Euxinian Basin (i.e., the Black Sea; Fig. 1A) that is contradicted by [START_REF] Gillet | La faune de mollusques du Pontien de trilophos (SW Thessaloniki)[END_REF] the older age of post-MSC marine incursion in the Dacic Basin than in the Black Sea [START_REF] Popescu | Pollen record and integrated high-resolution chronology of the Early Pliocene Dacic Basin (Southwestern Romania)[END_REF][START_REF] Popescu | Pliocene and Lower Pleistocene vegetation and climate changes at the European scale: Long pollen records and climatostratigraphy[END_REF], (2) the absence of an outlet into the Black Sea for the Palaeo-Tisza River during the peak of the MSC that did not morphologically favoured a marine invasion [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF][START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF], and (3) the barrier of the İ stanbul area which probably isolated the Aegean Sea from the Black Sea [START_REF] Melinte-Dobrinescu | The Messinian Salinity Crisis in the Dardanelles region: Chronostratigraphic constraints[END_REF]Suc et al., submitted). Accordingly, we take into account a marine gateway that crossed the Balkans region, of which few remains are today perceptible (Fig. 13).

For most of them, the deltaic deposits that constitute these remains belong to adjacent branches of the post-MSC marine gateway, as shown by the direction of their prograding sedimentation. They have been protected from the Quaternary erosion. Only the foreset beds of Turnu Severin and the deposits of Lakavica display an axial progradation. They thus benefited from an exceptional preservation.

4.3.Palaeogeographic evolution of the Balkans corridor between 6 and 5 Ma

Marine species of calcareous nannoplankton and/or dinoflagellate cysts have been recorded in the late Pannonian -early Pontian of three localities of the Pannonian Basin (Malunje, Krašić, Krajačići; Fig. 1A), i.e. between 8 and 7 Ma [START_REF] Kovačić | Lacustrine basin to delta evolution in the Zagorje Basin, a Pannonian sub-basin (Late Miocene: Pontian, NW Croatia)[END_REF]. The inferred connection with the Mediterranean Sea was defined by [START_REF] Popescu | Galeacysta etrusca complex, dinoflagellate cyst marker of Paratethyan influxes into the Mediterranean Sea before and after the peak of the Messinian Salinity Crisis[END_REF] as a western branch of the Balkans corridor between Niš and Belgrade, supported by regional geological maps (see: [START_REF] Kräutner | Geological map of the Carpatho -Balkanides between Mehadia, Oraviţa, Niš and Sofia[END_REF]. The marine influence was thus felt up to the western extremity of the Pannonian Basin (in the Mura-Zala Basin; Fig. 1A) where brackish ostracods became momentarily predominant in the late Pontian deposits before being replaced by freshwater species [START_REF] Stevanović | The Pontian sediments in Slovenia[END_REF], indicating a short lasting brackish water impulse. At that time, i.e. prior to the MSC, the Balkans corridor branched northward of Niš in direction of the Pannonian and Dacic basins, respectively (for detail, see: [START_REF] Popescu | Galeacysta etrusca complex, dinoflagellate cyst marker of Paratethyan influxes into the Mediterranean Sea before and after the peak of the Messinian Salinity Crisis[END_REF].

During the sea-level drop of the Mediterranean (peak of the MSC: 5.60 -5 46 Ma; Fig. 4; [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF], fluvial erosion prevailed up to the northern part of the Pannonian Basin (Csato et al., 2007[START_REF] Csato | A three-dimensional stratigraphic model for the Messinian crisis in the Pannonian Basin, eastern Hungary)[END_REF]. The peri-Carpathian Palaeo-Tisza River deeply cut the Iron Gates (Fig. 1B; [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF][START_REF] Leever | Messinian sea level fall in the Dacic Basin (Eastern Paratethys): palaeogeographical implications from seismic sequence stratigraphy[END_REF][START_REF] Suc | The Messinian Salinity Crisis in the Dacic Basin (SW Romania) and early Zanclean Mediterranean -Paratethys high sea-level connection[END_REF] and possibly continued within the former corridor in direction of the Aegean Basin, collecting a lot of tributaries and draining the residual Pannonian and Dacic lakes. Abundant clastics have been transported by the Palaeo-Struma and Palaeo-Mesta rivers into the Serres graben, and farther south, towards the North Aegean trough (Zagorčev, 1992a(Zagorčev, , 1992b;;[START_REF] Zagorchev | Late Cenozoic development of the Strouma and Mesta fluviolacustrine systems, SW Bulgaria and northern Greece[END_REF]. Such a phenomenon could have been at the origin of clastics that were deposited just above the MES in the Prinos Field (Fig. 7; [START_REF] Proedrou | Prinos Field-Greece, Aegean Basin[END_REF]. Prinos evaporites may also have experienced a similar event, but of less amplitude, during the first lowering in sea level of the MSC (5.971 -5.60 Ma), that would explain the clastics intercalations (Fig. 7; [START_REF] Proedrou | Prinos Field-Greece, Aegean Basin[END_REF].

As estimated by [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF], the Dacic Basin reconnected with the Aegean Sea at 5.45 Ma, just after the marine reflooding of the Mediterranean Basin. The high sea-level reconnection allowed crossed exchanges of organisms: dinoflagellates, ostracods and molluscs invaded the Aegean Sea and then the Mediterranean transported by surface water current (3 rd Lago Mare event; see: [START_REF] Couto | Lago Mare and the Messinian Salinity Crisis: Evidences from the Alboran Sea (S. Spain)[END_REF], while foraminifers, calcareous nannoplankton and dinoflagellates entered the Dacic Basin using a subsurface reverse current [START_REF] Clauzon | Influence of the Mediterranean sea-level changes over the Dacic Basin (Eastern Paratethys) in the Late Neogene. The Mediterranean Lago Mare facies deciphered[END_REF]. Influence of the invading marine waters into the Dacic Basin probably reached the westernmost part of the Pannonian Basin where brackish ostracods were again recorded in the Mura-Zala Basin (Fig. 1A; [START_REF] Stevanović | The Pontian sediments in Slovenia[END_REF]. The marine influxes reached the Black Sea at 5.31 Ma [START_REF] Popescu | Pliocene and Lower Pleistocene vegetation and climate changes at the European scale: Long pollen records and climatostratigraphy[END_REF][START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF].

It is possible to hypothesize that lagoons and lakes developed around the Balkans corridor during episodes of connection at high sea level, giving a sense to the so-called North Aegean Lake. Some of these scattered lakes may have been filled by brackish waters, providing another example of coeval nearby brackish (with dreisseinids, lymnocardiids and ostracods) and marine (with calcareous nannofossils) environments as evidenced by [START_REF] Couto | Lago Mare and the Messinian Salinity Crisis: Evidences from the Alboran Sea (S. Spain)[END_REF] in southern Spain.

According to [START_REF] Mǎrunţeanu | Mediterranean calcareous nannoplankton in the Dacic Basin[END_REF] and their calcareous nannoplankton investigations in the Dacic Basin, intermittent connections between the Mediterranean Sea and Dacic Basin may have existed up to 4 -3 Ma (Romanian Stage; Fig. 4). Using geochemistry, [START_REF] Sprovieri | Ba/Ca evolution in water masses of the Mediterranean late Neogene[END_REF] observed the inflow of Paratethyan low salinity waters into the Mediterranean from the Tortonian up to the Late Pliocene, except for the MSC duration. They hypothesized a Mediterranean -Paratethys connection during this time interval. The occurrence of Spiniferites cruciformis, a Paratethyan dinoflagellate cyst, at 4.2 Ma at Ptolemais (Fig. 1B; [START_REF] Kloosterboer-Van Hoeve | A short-term cooling event, 4.205 million years ago, in the Ptolemais Basin, northern Greece. Palaeogeogr., Palaeoecol[END_REF]) may be explained by water exchanges through the Balkans corridor. The search for calcareous nannofossils at Site 380 (SW Black Sea) evidenced two sets of discontinuous invasions [START_REF] Popescu | Pliocene and Lower Pleistocene vegetation and climate changes at the European scale: Long pollen records and climatostratigraphy[END_REF]: (1) from 5.31 to 4 Ma, expected via the Balkans corridor and the Dacic Basin [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF]; (2) from 2 to 0.9 Ma, possibly via the Marmara Sea formed since the propagation of the North Anatolian Fault [START_REF] Armijo | Westward propagation of the North Anatolian fault into the northern Aegean: Timing and kinematics[END_REF][START_REF] Armijo | Asymmetric slip partitioning in the Sea of Marmara pull-apart: a clue to propagation processes of the North Anatolian Fault?[END_REF][START_REF] Karakaş | Propagation et évolution structurale de la faille nord-anatolienne dans la région de Marmara[END_REF] and occasional flows over a proto-Bosphorus sill.

At last, propagation of the North Anatolian Fault into the Aegean Basin induced a new tectonic style resulting in the depression of the North Aegean Trough [START_REF] Armijo | Quaternary evolution of the Corinth Rift and its implications for the Late Cenozoic evolution of the Aegean. geophys[END_REF][START_REF] Westaway | Late Cenozoic extension in SW Bulgaria: a synthesis[END_REF] that probably deepened the Prinos peripheral evaporitic series. This crucial geodynamic change probably closed the Balkans corridor.

CONCLUSION

This synthesis of the available data on the Late Neogene deposits from the North Aegean Sea to the Dacic Basin emphasizes the impact of the MSC and of the subsequent marine reflooding. Before field and micropalaeontological investigation of the Balkans, Mediterranean marine and Paratethyan brackish organisms were both identified on the North Aegean coastline and in the Dacic Basin. Intermediate significant localities were found in Skopje and Niš while a passage over the İstanbul area was discarded.

It is concluded that a marine gateway crossed the Balkans region before the MSC thanks to a propitious extensional tectonic context. Its closure resulted from tectonically controlled subsidence to the north and south of the Skopje region, respectively. Such a gateway gives sense to the so-called North Aegean Lake: a mass of small lagoons and lakes probably surrounding the corridor during the high sea-level phases. During the MSC and particularly during its paroxysm, it was cut by fluvial erosion which probably originated from the Iron Gates and resulted in (1) transport of clastics in the North Aegean, and (2) making it easier to be flooded by marine waters at 5.46 Ma. The present-day dense and close fluvial networks of the region remember this corridor. At last, it appears essential to recall the importance of palaeogeographic reconstructions to decipher the Messinian Salinity Crisis.

ACKNOWLEDGEMENTS

Georges Clauzon, who passed away in March, 2013, actively participated to the field trips in the region with most of us. M. Böhme and A. Ilg are acknowledged for information on fossil fishes from the Prosilio section and for arousing our interest in erosional surfaces in this sedimentary succession. A. Di Stefano analysed our samples from the Akropotamos section. A, According to [START_REF] Hsü | Stratigraphy of the lacustrine sedimentarion in the Black Sea[END_REF]. B, According to [START_REF] Hsü | History of the Mediterranean salinity crisis[END_REF]Hsü et al. ( , 1978b)). C, According to [START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF].

Relief is from GeoMapApp© System http://www.geomapapp.org/ [START_REF] Ryan | Global Multi-Resolution Topography synthesis[END_REF]. Table 1. Results of some pollen analyses.

Fig. 1 .

 1 Fig. 1. Location maps. A, Location of the study area. Palaeogeography of the Central and Eastern Paratethys before and just after the Messinian Salinity Crisis is from Krijgsman et al. (2010), modified. Proposed gateway(s), subject of this paper, are not shown. Pannonian -Early Pontian quoted localities: a, Krašić; b, Malunje; c, Krajačići. Pannonian -Late Pontian quoted locality: d, Mura-Zala Basin. B, Studied localities. 1, Prosilio; 2, Trilophos; 3, Amphipolis Lion; 4, Akropotamos; 5, Nestos 1 and 2 wells; 6, Serres; 7, Lakavica; 8, Dračevo; 9, Gabrovačka; 10, Gura Vaii -Turnu Severin ; 11, Hinova.Relief is from GeoMapApp© System http://www.geomapapp.org/[START_REF] Ryan | Global Multi-Resolution Topography synthesis[END_REF].

Fig. 2 .

 2 Fig. 2. Previously mapped gateways within the Balkans.A, According to[START_REF] Hsü | Stratigraphy of the lacustrine sedimentarion in the Black Sea[END_REF]. B, According to[START_REF] Hsü | History of the Mediterranean salinity crisis[END_REF]Hsü et al. ( , 1978b)). C, According to[START_REF] Bache | A two-step process for the reflooding of the Mediterranean after the Messinian Salinity Crisis[END_REF].Relief is from GeoMapApp© System http://www.geomapapp.org/[START_REF] Ryan | Global Multi-Resolution Topography synthesis[END_REF].

Fig. 3 .

 3 Fig. 3. Present-day fluvial network in the Balkans region, with opposed drainage in direction of the Aegean Sea on the one hand, of the Danube River on the other hand.Relief is from GeoMapApp© System http://www.geomapapp.org/[START_REF] Ryan | Global Multi-Resolution Topography synthesis[END_REF].

Fig. 4 .

 4 Fig. 4. Standard global nannoplankton bio-chronostratigraphy (Raffi et al., 2006) with respect to the Global Polarity Time Scale (Lourens et al., 2004) from 8.5 to 3 Ma. NN Zones are from Martini (1971) amended by Perch-Nielsen (1985) and Berggren et al. (1995). See also Melinte-Dobrinescu et al. (2009) for more details concerning this nannoplankton chart. Correlations with Paratethys Stages are from Vasiliev et al. (2004) and Krijgsman et al. (2010).

Fig. 5 .

 5 Fig. 5. Record of the MSC at Akropotamos (Northern Aegean). A, GoogleEarth three-dimensional view showing the location of the studied section (yellow line) and localities near Akropotamos (relief is exaggerated three times). B, Interpretative sketch of the Akropotamos section. 1, Miocene continental deposits (sands, conglomerates); Messinian: 2, Gypsums; 3, Clays; 4, Limestones; 5, Messinian Erosional Surface; Post-MSC Gilbert-type fan delta: 6, Clayey bottomset beds; 7, Sandy and conglomeratic foreset beds; 8, Topset beds (limestone). C, Gypsum quarry northward of Akropotamos. D, The Messinian Erosional Surface separating the Gilbert-type fan delta foreset beds from the Messinian limestone along the road to Kariani. E, Detail of this contact. F, The Messinian Erosional Surface separating the Gilbert-type fan delta foreset beds from the Messinian limestone southward of Akropotamos. G, Proximal part of the Akropotamos Gilbert-type fan delta. H, Distal part of the Akropotamos Gilbert-type fan delta.

Fig. 6 .

 6 Fig. 6. Examples of the Messinian Erosional Surface. A, Near the statue of the Amphipolis Lion. B, Detail of the foreset beds at the same place. C, Near Trilophos. D, Marine clays topping the Trilophos Formation at the same place. E, View of the middle part of the Prosilio section. F, Detail of the same section. G, Detail of photograph F.

Fig. 7 .

 7 Fig. 7. Re-interpreted synthetic section through the Prinos hydrocarbon field (from Proedrou and Sidiropoulos, 1993).

Fig. 8 .

 8 Fig. 8. The Serres area. A, GoogleEarth three-dimensional view showing the location of the Serres Zanclean Gilbert-type fan delta nested within the Miocene piedmont of the Vrondou Mount (relief is exaggerated three times). B, View of the foreset beds of the Serres Zanclean Gilbert-type fan delta nested within the Miocene piedmont. C, Zoom on the foreset beds of the Serres Zanclean Gilbert-type fan delta nested within the Miocene piedmont. D, The Messinian Erosional Surface. E, Bottomset beds with lignite beds of the Serres Zanclean Gilbert-type fan delta.

Fig. 9 .

 9 Fig. 9. The Dračevo Zanclean Gilbert-type fan delta. A, Skopje area: GoogleEarth three-dimensional view showing the location of the Dračevo Zanclean Gilbert-type fan delta with exposed sections near Batinci and Varvara (relief is exaggerated three times). B, Double marine-continental transition. C, Sandy to conglomeratic foreset beds. D, Clayey bottomset beds. E, Conglomeratic topset beds with travertine blocks.

Fig. 10 .

 10 Fig. 10. The Lakavica series. A, GoogleEarth three-dimensional view showing the location of the sections (relief is exaggerated three times). B, Prograding deposits. C, Detail of the conglomeratic prograding deposits. D, Transition from prograding to aggrading deposits. E, Aggrading deposits.

Fig. 11 .

 11 Fig. 11. The Niš area. A, GoogleEarth three-dimensional view showing the location of the Gabrovačka sections (relief is exaggerated three times). B, Clays with marine microfossils at Gabrovačka Reka. C, Lacustrine clays back the Gabrovačka stadium. D, Overlying fluvial conglomerates at Gabrovačka Reka.

Fig. 12 .

 12 Fig. 12. Major features of the Turnu Severin Gilbert-type fan delta. A, GoogleEarth three-dimensional view showing the location of the illustrated sections (relief is exaggerated three times). B, Detail of the conglomeratic foreset beds. C, Gravely and sandy foreset beds. D, Clayey bottomset beds. E, Interlocked bottomset and foreset beds at Pd. Trestenic. F, Detail of the clayey bottomset beds.

Fig. 13 .

 13 Fig. 13. Proposed course of the Balkans corridor with the present-day watersheds A, B, and Cwith respect to the Strymon, Vardar, Morava and Timok drainage basins, respectively (see Figure3).

Fig. 14 .

 14 Fig. 14. Present-day altitudinal profile of the proposed Balkans corridor with highest places corresponding to the drainage divides (A, between the Strymon and Vardar systems; B, between the Vardar and Morava systems; C, between the Morava and Timok systems) on which is reported the present-day altitude envelope (in blue) of the marine-continental transition (i.e., the Zanclean seashore) in the six evidenced Zanclean deltaic systems.

The main part of fieldwork and micropalaeontological analyses was realized during the ANR-EGEO Project that provided a significant financial support. Two anonymous referees significantly contributed to improve this manuscript.