Overview of the Bosphorus Depositional Fan from Data Sets Recovered on the Black Sea Shelf off the Strait of Istanbul

Since the preliminary mapping of the northern Bosphorus outlet of the Black Sea continental shelf by the 90’s NATO project, an important data set bringing morphological, seismic and sedimentological information is now available to describe there a shallow fan system consisting of a main channel and short distributary branches, levees, and sediment waves, correlative to a southern source flowing to the north. This paper brings complementary quantitative data from landscape at the eastern and western shelf area of the shallow fan system. It confirms that the Bosphorus outlet fan preserved on the shelf surface provides an important record of surface water flow. To deduce the history of the Marmara/ Black Sea last connexion, it is fundamental to understand how surface water flow could have generated the observed morphology. To date, morphological studies propose only very short time estimates for the longevity of the necessary hydrological events responsible for this deposit. Here we propose a compilation of multibeam bathymetry data from previous studies completed with multibeam bathymetric data collected from new areas not previously studied, together with high resolution seismic data (Chirp) and gravity cores from the confluence of the Bosphorus and Black Sea. This compilation shows that this shallow fan delta is in all points comparable to a depositional fan with a branching pattern of channels mimicking a mostly exposed fan-shaped body prograding from the foot of a significant break in slope and extending from the outlet of the strait to the shelf edge. Although we do not entirely resolve the ongoing controversies, we provide extra data to obtain a wider view to this key area. Using analogue cases and flume studies, it is acceptable to propose that the distinct morphology of this form of fan delta could have been the result of a single event of basin-filling on a timescale of decades after a relatively strong Mediterranean-originated northerly flow into the partly sub aerially exposed shelf area earlier than 8.5ky BP reshaped by present day saline undercurrent processes.


Introduction
Although the structure of Bosphorus outlet submarine channels had been visually evidenced by analogue echo-sounding [1], the first multibeam images were obtained in 1999 by the NATO RV Alliance and the Turkish Navy survey ship Çubuklu, in the framework of a NATO SACLANT Undersea Research project [2]. These shelf features were initially described as grooves by Özsoy et al. [3]. In 2002, a survey carried out on board the Ifremer RV Le Suroit for BlaSON project [4] increased the multibeam mapping of this underwater channel fan-delta and canyons system. Since then, several other surveys were performed to complete this Black Sea anastomosed channel network mapping [5,6], also described as a subaqueous shelf fan [7] or a channelled shelf fan [8]. The mapped set of shelf-cut channels reaches the heads of recent submarine canyons which was recognized in 2002 ( Figure 1) [4,9].
The Bosphorus (or Istanbul Strait) provides a sole outlet for the semi-enclosed Black Sea basin, connecting it to the Marmara Sea. It is almost 30km long and 3.4km wide at the northern outlet and 2.8km at the southern inlet. In between, the Bosphorus strait displays a minimum width of 700m [10,11]. Its average depth is 65m varying from approximately 15m (depth of the channel rims) to 110m, its deepest part located in midstream of the channel. The inner channel begins off the southern sill (27m deep) and continues to the Black Sea outlet where it reaches 105m ( Figure 2).

Two oppositely flowing currents
a) The upper layer with a salinity ranging from 18 to 20psu coming from the Black Sea and b) The lower layer which salinity is ranging between 36 and 38psu coming from the Marmara Sea, characterize the physical oceanography of this conduit. Flow exchange is affected mainly by the hydraulic conditions generated by the geometry of the strait. One specific water mass through the strait is the cold intermediate water (CIW) observed in the Black Sea during the summer months [12,13]. The warm and more saline lower layer, called Mediterranean Sea Water (MSW), flows into the Black Sea [3,12,14]. For Özsoy et al. [3], the shelf edge structure of the delta and the channel morphology, i.e. curved and overflow channels, evoke river morphology, which corresponds today to the deep salt water flowing out of the Bosphorus channel into the Black Sea. For these authors, it was not clear whether these features were related to an ancient fresh-water stream.
The main sill corresponds to the shallowest part of the Bosphorus strait and is present at its southern entrance as a sedimentary deposit forming an elongated sill reaching a depth of only 27m [10] (Figure 2). This sill had the ultimate control of the Marmara Black Sea connection. The Istanbul Strait multibeam bathymetric data ( Figure 2) indicate that the morphology of the Strait is deepening northwards and its seafloor morphology is mainly characterized by erosion [10].  [10]. The elongated Sill at the southern entrance of the Bosphorus strait corresponds to its shallowest part. These multibeam data clearly show that the morphology of the Strait is deepening northwards.

International Journal of Environmental Sciences & Natural Resources
The Black Sea has been subjected to various depositional and erosional events due to the changing global sea level, during glacial and interglacial stages of Pleistocene. These sea level changes of the Black Sea are complicated and do not fully track the global sea level, due to its regional climatic conditions, confined geographical conditions and limited connection to the world ocean [15,16]. Its brackish-fresh conditions have changed to marine condition, since Holocene with the northward flowingundercurrent of Mediterranean waters from the Bosphorus [17][18][19][20][21][22]. The last re-entrance of the Mediterranean waters into the Black Sea has been controversial and a matter of debate in terms of its mode and timing in the published studies carried out on the North-western and South-western shelf of the Black Sea, and in the Marmara Sea, from 1997 to 2011 [23][24][25][26][27][28][29][30][31][32][33][34]. Since 2011, the date of the last reconnection seems to have been widely agreed and would have occurred with the Initial Marine Inflow (IMI) dated around 8.9 -8.5ky BP [16,[35][36][37] or 9.3-9ky BP [38,39]. The common results of some of these studies indicated a shelf-wide erosional surface with prograding wedge system (lowstand) of lacustrine shelf deposits at the shelf edge. This unconformity is overlain by various bedforms and sedimentary deposition and local erosional unconformities. The sedimentary deposits and bedforms on the shelf area were interpreted as formed in a shallow water-coastal environment [23,[28][29][30]40].
The confluence of the Bosphorus and Black Sea at the inner shelf area is the critical site for understanding the processes that occurred in the recent past. Flood et al. [5] published a detailed morphological map of the channel network from the Bosphorus outlet to the shelf edge and discussed the sedimentological process in its evolution. Hiscott et al. [6] have documented the interaction of the overspill at the outlet with sea bed topography. And recently, Ryan et al. [8] considered that the synchronicity between the arrival of Mediterranean fauna and the onset of the Black Sea outlet fan construction suggests an origin linked to the connection of the Black Sea's lake with the global ocean.
In this study, we propose a compilation of multibeam bathymetry data from afore mentioned studies completed with multibeam bathymetric data collected from the surrounding areas, together with high resolution seismic data (Chirp) and gravity cores from the confluence of the Bosphorus and Black Sea. These data were collected in 2007 and 2008 by RV K. Piri Reis of Dokuz Eylül University, in the frame of a research project supported by TÜBİTAK, under the Programme of PIA-BOSPHORUS. The present detailed morphological and seismic data obtained at this particular confluence (Figure 1), will be compared to analogue features as alluvial fan or eroding rill network from the berm to the beach to supplement the previous published studies on the Bosphorus outlet shelf fan. These features are normally fleeting concentrated flow paths functioning not only as a source of sediment but also as sediment distribution systems creating slope erosion. They are interpreted as being formed by a stream competence loss flowing onto a flood plain or to a high-gradient stream reaching a body of water [41]. Such features form when the flow reaches its basal level and, as the power of the stream slows down, sediments are deposited at the foot of the slope when rapid water flows penetrate a shallow area. The dynamics then slow down and the sediments settle in the bed of the streambed creating patterns similar to the beach patterns visible on the beach when tide recedes ( Figure 3).

The bosphorus channel
The Bosphorus (Istanbul Strait) is a zigzagging channel formed by conjugated pairs of oblique faults inside a meandering strait [42]. The geological evolution of the Bosphorus is formed by combine effects of tectonics of the region and a pre-existing fluvial system, and has an opening history that starts during the Pliocene with a series of tectonic movements of the North Anatolian Fault (NAF) extending to the Marmara Sea. Meandering physiography of the strait was shown to be resulted from the trends of the major faults along the both sides of the Bosphorus by various investigators [11,42,43]. The blocks bounded by conjugated pairs of oblique faults started to move away from each other, initially forming a channel, the Bosphorus began to shape as a zigzagging waterway, when the North Anatolian Fault extended into the Marmara region during the Pliocene-Pleistocene period [42]. The basement in the Bosphorus channel is constituted by a sediment infill locally reaching a thickness of more than 130m in the southern part overlying a Paleozoic and Upper-Cretaceous substratum. Previous drilling borehole data in the Central part of the Bosphorus indicate the intermittent deposition over the basement rock during the last 26ky BP with the appearance of first euryhaline molluscs dated as young as 5.3 14C ky BP [44,45], suggesting an erosional event. The Bosphorus strait multibeam bathymetric data ( Figure 2) revealed that the upper portions of the deeper channel deposits are composed of parallel layers, which were mainly eroded after their deposition [10]. The last erosion surface corresponds to the bottom of the strait, and the volume of eroded material is estimated at 2.108m 3 and km 3 . Such volume of eroded material excavated from the Bosphorus Strait is comparable in magnitude to the volume of material constituting the fan. The present day seafloor morphology indicates that the direction of the flow at the origin of this erosion was from the south to the north. The latest morphology of the strait floor was established at about 5-4ky BP. by the Mediterranean bottom current of the modern two-way flow system [10].

Oceanographical setting
Previous studies, aforementioned, indicate that the prominent cut channel of the shelf fan is directly connected to the Istanbul Strait. This channel turns abruptly northwest when it reaches the mid-shelf of the Black Sea and the stream is deflecting to the left. This is in contradiction with the theory demonstrating that geostrophic currents exiting an outlet in the northern hemisphere and flowing on a sloping sea-floor are deflecting to the right. But no observation was presently able to recognise a Mediterranean plume flowing eastward of the Bosphorus Strait outlet. On the other hand, Stanev (2005) and Stanev et al. (2011) [46] have clearly evidenced a plume persistently flowing northwestwardly within the Black Sea shelf.
Today concerning the surface water circulation, the Black Sea is divided into two basins characterised by two large anticlockwise cyclonic central gyres. These gyres are marked out by a wind-driven current flowing around the whole basin, along the continental shelf edge and above the continental slope [47]. Such anticlockwise "Rim current" flows eastward along the coast of Anatolia on the South western Black Sea shelf in separating anticyclonic coastal eddies from the central basin cyclonic gyres [48][49][50] (Figure 4). The basin-scale model "Black Sea Physics Analysis and Forecast" from MyOcean catalog (http://marine. copernicus.eu/) can be used for continuous analysis and forecast of the Black Sea circulation and stratification. The model output reveals strong mesoscale eddies playing an important role in the sediment transport from coastal zones to deep water (the mean rim current is represented Figure 4). For instance, several anticyclonic eddies of interest for this study were described near the coast. Particularly, the sub-basin scale eddies located near the Bosphorus strait outlet and along the Crimean Peninsula play an important role in the present day sediment dynamic. The ARGO float profile data published by Stanev et al. [51] indicate that the current velocity measured at 200m depth along the "Rim Current" is about 15cm.s -1 and decreases to 5cm.s -1 at depths of 750m and 1500m [52]. This slower deeper current closely follows the topography along the Black Sea southern margin. The numerous quasi-permanent coastal eddies described outside the "Rim Current" (Figure 4) are interpreted to be formed from upwelling and wind stress curl mechanisms occurring along the coastal shelf. These data have also shown strong salinity gradient stem from the contrast between the more saline and denser lower layer and the low salinity surface layer.

International Journal of Environmental Sciences & Natural Resources
The Black Sea has layers of deep and surface waters that do not intermix and therefore is considered as the largest meromictic basin in the World. The freshwater input to the Black Sea is supplied from smaller rivers such as the north Turkish margin rivers and essentially from large rivers: i.e. the Danube River (ca 200km 3 .yr -1 of mean water discharge) and the Ukrainian rivers Dniepr, Southern Bug and Dniestr (about 65km 3 . yr -1 of mean water discharge) [21]. Altogether, the net freshwater input to the Black Sea from the river added to the precipitation is estimated to reach about 300km 3 .yr -1 . In addition to this net freshwater inflow, around 300km 3 .yr -1 of higher salinity water is provided by the deep inflow entering the Black Sea through the Bosphorus [53]. This net input is compensated by about 600km 3 . yr -1 of low salinity surface water outflowing the Bosphorus to the Marmara Sea [53]. At around 50-150m, the upper and lower water of the Black Sea are separated by a steep density gradient (or pycnocline) originated from the gradual increase of salinity with depth below 200m. The saline water inflow entering the Black Sea from the Bosphorus mixes with the Cold Intermediate Layer (CIL) mainly before reaching the shelf break. This CIL located between the pycnocline and the warmer surface water coming from the summer river input, is formed during winters when the upper layer of relatively fresh water of the Black Sea becomes cold. During cold winters more Cold Intermediate Layer water is formed. The outflowing plume out the Bosphorus is mixed downwards and contributes to the deep-water salinity [47].

SW shelf area
Previous published seismic studies carried out in the Black Sea shelf area indicated a widespread shelf-crossing erosional surface [7,23,[28][29][30][31][32][33][34]39,[54][55][56], extending to the shelf edge at about 90-120m water depth. This erosional surface "α" (α and α1 terms are from Aksu et al. [23] is almost flat and corresponds to subaerial exposure of the shelf during Pleistocene lowstand stages of the Black Sea. The sediments unconformably deposited over this "α" surface have been defined as Unit 1. Underneath the "α" surface, four units are recognized: Unit 2, 3, 4 and 5 ( Figure  5) [7,23]. Unit 1 lying above α surface and divided into two subunits, 1b and 1a, has been described by Aksu et al. [23]. Subunits, 1b and 1a are separated by the second discontinuity surface, α1 ( Figure  5). These subunits have been interpreted by Aksu et al. [23] as lowstand, transgressive and highstand systems tracts deposited during the last glacial low stand and Holocene by a gradual sea level rise beginning ~11 14 C ky BP. Such relatively thick deposits with parallel and almost even reflections, including prograding clinoforms are also observed off the Sakarya River in the southern Black Sea shelf [54].

Bathymetrical data
Data completing the Bosphorus outlet survey (A, B, D areas on Figure 1) are derived from Multibeam bathymetric data collected: a) On board of RV Le Suroit, in 2002 using a SIMRAD-EM 300 echo-sounder with a 300kHz frequency and a single Mills Cross array geometry. Swath widths are selectable using either 111 (shallow mode) to 135 beams. The EM 300 system operates from 10 to 5000m of water depth with an International Journal of Environmental Sciences & Natural Resources optimum use between 30m and 500m. This bathymetric data set was processed, using Ifremer "Caraibes" and "Globe" softwares. Contour extraction, "spline" curves and bi-dimensional digital filtering were applied to remove multibeam artifacts (C on Figure 1). b) On board of RV K. Piri Reis with a SeaBeam 1050D echosounder having a frequency of 180kHz and providing 126 beams. Data processing was also performed with Ifremer "Caraibes" and "Globe" softwares and compiled with the other bathymetrical data (E and F on figure 1).

High resolution seismic data (CHIRP)
Locations of the seismic profiles used for this study are displayed Figure 6. Two profiles P36 ( Figure 5) and P37 ( Figure  7), were shot during BlaSON 2 survey [57] using the hull-mounted CHIRP sonar single channel of RV Le Suroit. This system has an operating frequency ranging from 1.8 to 5.3kHz and acquisition was done using a Delph PC based system. Obtained data were processed and presented in envelope traces. Part of profile 36 has been published by Okay et al. [7].  SeaBed 3010 and Bathy 2010 CHIRP sub-bottom profilers were used on board the R/V K. Piri Reis during cruises in 2007 and 2008, respectively. The Seabed Model 3010 is a threechannel acoustic profiling system offering three transmitting frequency ranges: 3 -7Khz, 10 -14kHz and 200kHz. The system was deployed using the first mode. Bathy 2010 system has 9 transducers with a sweep frequency from 2.75 to 6.75kHz and performs real time heave corrections. Collected data were band pass filtered and illustrated as amplitude envelopes.

Sediment cores
We studied thirteen gravity cores collected on board the R/V K. Piri Reis and the R/V Arar (Istanbul University, Institute of Marine Sciences and Management) and labelled Kpx, and two box-cores (labelled B2KIxx) recovered by the R/V Le Suroit. Core coordinates, length and water depth at location are presented in Table 1. The gravity cores and the box corers were photographed and lithological described. Subsequently, seven cores (Kp1, Kp4, Kp5, Kp11, Kp12, B2KI03 and B2KI05) ( Figure 8) were systematically sub-sampled at 2cm intervals. For benthic foraminiferal analysis, samples every 10cm were selected and for each sample about 10g of dried sub-samples were soaked in water for 24 hours, and then wet sieved through a 0.063mm sieve. The residue was dried in air and split with a microsplitter, where foraminifers were abundant, until a subsample containing minimum 300 benthic foraminifer specimens was obtained. In the case of low density (less than 300 specimens), the total residue were examined. All the benthic foraminiferal specimens were identified and counted under a binocular microscope. Age determination of selected mollusc samples was carried out by radiocarbon dating at the Isotope Geochemistry Laboratory of Arizona University and Geochrone Laboratory (USA) ( Table 2). For Table 2, 14 C ages are calibrated using IntCal13 [58] and program Calib Rev 7.0.4 [59]. Reservoir age is assumed to have been constant to 415 ± 90 14 C yrs [60] through the Holocene. Nevertheless, all ages reported in the text are uncalibrated radiocarbon dates with no reservoir correction, in order to compare with previous works in the area.  Table 1. Water depth (m) of cores is given within parentheses below the core labels. Red lines indicate sharp lithological boundaries within the cores. Ages are indicated by "ky" corresponding to 14 C age BP. Correspondence to Calibrated yr BP is given in Table 2.   14 C ages corrected for reservoir age were calibrated using IntCal13 [58] and program Calib Rev 7.0.4 [59].

Bosphorus outlet morphology
The morphological features observed in the study area are briefly defined from the bathymetrical map ( Figure 1) and discussed together with the high resolution CHIRP profiles. From a physiographical point of view; the study area can be defined into two parts: (1) the Bosphorus Channel, and (2) the Shelf area.

Bosphorus channel
The Bosphorus channel at the outlet of the Istanbul Strait is incised to a depth of about -85m, as a narrow ravine and extends into the inner shelf (Figure 1 & 2). A sill with elliptical shape located at a depth of -60m within this ravine has a length of 4.2km. The channel extends 8 km along N-NE direction, turns sharply to NW direction. From the abrupt bend, it is divided into branches with a prominent one flowing northward to the shelf edge. There are short branches to NW and N-NE. The N-NE trending branch is bifurcated, terminating at about -90m. It is blocked by levees and hence, it might have been formed before the main channel (avulsion). The channel becomes shallower and is carved into shallower depths, about 1 to 3m deep from mid-shelf to shelf edge. Flow traces and small bars can be noticed in this shallow NW extension of the channel. Okay et al. [7] have described these branches according to their formation/activity order. Levees occur at both sides of the channel, particularly after its sharp turn to NW direction, from -60 to -70m water depths.

International Journal of Environmental Sciences & Natural Resources
Seismic profiles indicate that the channel axis has irregular topography and an asymmetrical profile, from the outlet of Istanbul Strait to its sharp turn to NW (Figure 1). The western side of the channel bottom becomes deeper. The banks along the channel at -20 to -45m water depths have in general, acoustically weak reflections (Figure 9). From seismic profiles P5 to P8 (Figure 9 & 10), these banks can be seen as rock formations rising to sea floor, particularly at the east side of the channel. Reflections from the channel bottom itself also do not show any sedimentary deposit. The west side levees have some internal reflections and display the same character towards the north. Drapes in both side of the channel characterised by an acoustically transparent thin veneer, are also noticed where the channel shallows ( Figure 10). The basement is exposed on the sea floor, indicating the active erosional process of the channel (Figure 10). A conical mound which might be induced by rising gas or mud [5,23] is also found on the right side of the shallow channel (Chirp profile P17, Figure 10).

Shelf area
The shelf area where the channel incised is a smooth plain includes linear bedforms (ridges), sediment waves and conical mounds ( Figure 11). The linear bedforms at the east of the channel are ranged in the same direction of the channel at a water depth of 75-80m. They have 1 to 5m heights and 5.5-6km lengths. At the west side of the channel, linear bedforms are located at 80-85m water depths and slightly curved. They have up to 5km length. Apart from these morphological features, small hollows cut to a depth of 90-95m are observed at -85m water depth (Figure 11).
The linear ridge-shaped elements at the eastern and western shelf area of the channel appear to be substratum outcropping on the sea floor, and/or accumulating sediments over such

International Journal of Environmental Sciences & Natural Resources
topographical protrusions (Figure 10 & 12). The sediment waves in the vicinity of the levees have mostly asymmetrical shape with their steep sides leaning to N-NE (Figures 10 & 12). They have different heights from the sea floor varying from 3 to 8m and present internal stratifications ( Figures 5, 7 & 12).

Seismic stratigraphy of the SW shelf
The almost flat widespread shelf-crossing erosional surface here-before mentioned ( §2.3, surface "") is clearly recognized on the profiles of this study. Four to five seismic units are truncated by this unconformity (Figure 5). Some specific features are noticeable on the seismic profiles, i.e. Unit 3 is incised by paleochannels. Unit 2 is only observed at the shelf edge. Unit 1 is the youngest unit, deposited above α surface and International Journal of Environmental Sciences & Natural Resources divided into two subunits, 1b and 1a. As described by Aksu et al. [23], Unit 1a is clearly separated from Unit 1b by the second discontinuity surface, α1 ( Figure 5 & 7). This minor erosional surface "α1" is observed in the vicinity of the channel and is less apparent, compared to "α" surface. Unit 1b is represented by a lowstand shelf wedge (thickness; 20ms) at the shelf edge, deposited over previous wedges ( Figure 5B). It has parallel and almost even reflections and is relatively thick, and includes prograding clinoforms. However, it becomes very thin under a chaotic reflection toward the inner and mid-shelf, where it has a limited occurrence ( Figure 5A). Unit 1a is defined between α1 surface and present seafloor ( Figure 5 & 7).
The seismic stratigraphical relation of the channel-levee complex can be better seen after its sharp NW turn. Major erosional surface α is mostly flat, but can be undulating too ( Figures 5, 7 & 12). α1 is a local and weak discontinuity and intersects α surface beneath the levees, where the inclined structure of the basement displays irregularities or undulations. Unit 1b is very thin and located between these two surfaces, and is absent where two erosional surfaces merge (Figures 5 & 10). However, at the SW outer shelf of this study area, this unit is observed as parallel reflectors. The internal structure of the levees is mainly chaotic, but locally displays inclined reflectors at their top. The maximum thickness occurs near the channel, decreasing slowly with increasing distance from the channel (Figures 5B & 9).

Sediment data Additional information from cores
The locations of the sediment cores are shown in Figure 6 and Table 1 gives information on their coordinates, length and water depth at each location. Some of the cores did not recover materials due to the nature (sandy and/or substratum) of the sea-bottom. However, these last cores confirm the seismic interpretations.
On the left side of the channel after its sharp bend at a water depth of 63m, Core Kp12 (Chirp Profile 7, Figure 9) recovered the thin sediment layer covering the substratum, corresponding to the weak reflection interpreted on the seismic profile. Core Kp12 is 220cm long and has three different lithological units (Figure 8). From the bottom, there are 20cm of gray-green mud in which a Dreissena shell (from 214-211cm) yield an age of ~7.8 14 Cky BP (Table 2). In the interval from 200 to 173cm , silt and sand content increases. At 173cm there is a sharp boundary and above the discontinuity a gray-green shelly mud continues until 58cm. A Mytilus shell from this layer (130-128cm) was dated to ~6.2 14 C ky BP ( Table 2). The top 58cm is laminated black-gray mud, Core Kp5 was collected from a water depth of 80m (Chirp Profile P11, Figure 10), next to a sediment wave and above the α1 surface. From the bottom, 75cm to 50cm a gray mud with large shell characterises the lithology. Above 50cm mud with abundant shell fragments occurs up to 20cm. The uppermost 20cm contains dark gray mud with fine sand (Figure 8).
Core Kp4 was recovered from 77m depth over a small sediment wave (Chirp Profile P10, Figure 10). It consists of 14cm of sandy mud at the bottom overlapped after a sharp boundary by 145cm of gray-green mud with occasional shells and shell fragments. Two ages obtained at 141-138 and 145-142cm gave ages of 6.8 and 6.6 14 C ky BP (Table 2).
Cores Kp8 and Kp9, respectively recovered in the Channel and on the lee side of the channel at the outlet of the strait (Figures 9 & 10), recovered only a small amount of sand. Core Kp11 (133cm of recovery) realised in 60m of water depth and located at the right side of the channel (Figure 9, profile P5) consists of gray-green sandy muds (Figure 8). This layer is the thin sediment cover over the substratum, interpreted from Chirp profile P5 (Figure 9). Core Kp6 on the right side of the channel (Chirp Profile P8, Figure 10) has only 5cm-length and consists of sand. These cores could not penetrate further into the substratum.
Core Kp3 (Chirp Profile P12, Figure 10) recovered only a small amount of shell hash, probably reflecting the erosional surface.
Core Kp1 represents the thin veneer (Unit 1b) from 93m water depth above the erosional surface α (Chirp Profile P17, Figure 10). The lowermost 100 cm is gray-green mud gradually changed upward in a laminated black gray mud (Figure 8).
Short box cores collected during BlaSON 2 survey: B2KI03 and B2KI05 (~80cm of length, Figure 8) all consist of gray green mud, including large Mytilus shells and their fragments.

Benthic fauna
Benthic foraminiferal assemblages within core Kp12 follow more or less the lithological boundaries ( Figure 13). Total benthic foraminifera density (TBF) is low from the bottom of the core until cm 125 where it first increases between 120 to 110, followed by a second higher peak between 80 to 60cm ( Figure 13). The uppermost laminated black-gray mud contains the lowest total benthic foraminifera density and poor diversity. Foraminifera diversity decreases from the bottom to the top of the core. This assemblage change is also observed in mollusc (pelecypoda) fauna. A fresh-brackish form, Dreissena, and juvenile euryhaline molluscs are found together below 140cm, whereas only

International Journal of Environmental Sciences & Natural Resources
Mytilus shells are found above this level. The lowermost 80cm of the core represents the transition from fresh-brackish to marine conditions, based on changing mollusc fauna.
In comparison with previous data, benthic foraminiferal abundances of the short box core B2KI03 ( Figure 13) Figure 13). Core B2KI05, which was analysed and the results published by Kirci-Elmas et al. [61], show a higher diversity (number of species: 25-43) compared to other cores and its total benthic foraminifera reduce at the upper parts of both cores.

Discussion
Previous studies, mentioned above, indicate that a prominent cut channel on the seabed has a direct connection with the Bosphorus outlet. This channel turns abruptly northwest when reaching the Black Sea mid-shelf, but in the northern hemisphere, the theory of geostrophic currents demonstrates that the stream exiting the strait and flowing on a sloped bottom should have deflected to the right. Today, all observation carried out in this area have failed to trace the plume of the incoming water from the Marmara Sea flowing eastward. For instance, Stanev et al. [12] & Stanev [46] have clearly evidenced a persistently northwestward plume within the shelf.  [1,2,49,62] have initially proposed that the deflection to the left of the plume could be the result of a guidance from the narrow underwater channel. As already mentioned, this deflection could be related with the anticlockwise rotating 'Rim current' flowing eastward along the Turkish coasts. The anticyclonic coastal eddies are consequently separated from the cyclonic gyres located in the center of the basin [48][49][50] (Figure 4). In the present study, the seafloor morphology confirms that after the plume reaches the continental slope, it turns to the left, and when reaching the shelf break, entrainment dominates the plume dynamics leaving very little propagation along the slope.

A shallow fan
Submarine channels in lakes, seas or open oceans can often be kilometres wide extending for thousands of kilometres across the seafloor and reaching hundreds of metres deep. Their formation results from underwater flows of density currents along the sea floor; like lake water, water denser than sea water or by flows of sand or mud. These channels correspond to the main transport pathway to the deep for moving sediments and are very important as they form large deposits of sediments. Today, laboratory experiments are manly the only way to understand such flow processes and the knowledge of their deposits is obtained from ancient examples now exposed on land or from recent analogue studies like this study.

International Journal of Environmental Sciences & Natural Resources
Holmes & Nemec and Steel [63,64] have defined a "Fan delta" :"a prism of sediment, delivered by an alluvial fan, and deposited mainly or entirely subaqueously at the interface between the active fan and a standing body of water". Such depositional systems are produced in both low and high energy conditions, containing coarse and fine-grained deposits, erosional bedforms (chutes, flutes) and hence are related to episodic events occurring under both riverine and submarine processes [65,66]. Modern fan delta formations are found in general, in front of tectonically active coasts with high relief, and on wave-dominated shelves [63,66] and also can be formed with a sudden catastrophic flood [67,68]. The seismic stratigraphical relations together with the morphology indicate that the present extension of the channel on the shelf area and associated sedimentary depositions are geologically very young and were developed by a relatively fast process, as they were well preserved [39]. This does not imply to the origin of the Istanbul Strait itself, which dates back to Pliocene [11,42] Incised and filled paleochannels observed in older units below α (Figures 7 & 9) surface suggest a slightly different route or routes than that of present, during Pleistocene glacial periods [7]. The rock formations exposing on the sea floor on the east side of the channel (Figure 9) might have caused the channel to redirect its route from NE to NW. Around N41°29' and E 29°02' (Figure 1), mid-channel islands are recognizable, whilst Hiscott et al. [6] have described modern day surficial sediment processes evidence, these islands display classic lemniscate or tear-drop morphology and are generally formed during highmagnitude water flow [69,70].
The presented fan delta defined previously by Giosan et al. [24]; "Bosphorus-shelf fan") and Okay et al [7]; "shelf fan") consists of a channel with asymmetrical section and is relatively wide at the strait outlet, presenting levees located on its both sides where some avulsion channels cut previous levees, downstream. The levees are characterized by mostly chaotic seismic reflectors. These relatively thicker deposits are accumulated along the main channel, decreasing their thickness away from the channel. Asymmetrical sediment waves occur mainly on the right side of the channel. As the channel becomes shallower and branches towards the north, levees appear to be cut by the active channel and its short branches (distributaries) [7].
Core data shows that the proximal area, low amplitude regions of the fan (channel outlet) are eroded to the substratum and/or overlaid by coarse grained (sand-rich) sediments. In this area, deposition and/or preservation of mud must have been prevented by the enhanced erosion. Levees and sediment waves probably consist of relatively coarse-grained sediments accompanied by shells, and are capped by mud drape (cores Kp12 and Kp4). At the distal part of the fan, deposition is characterized by high amplitude and continuous seismic character. There are various mud depositions related to lateral accretions and ongoing activity of the channel at the outer shelf [5].
The morphology of this shallow fan system indicates a southern source flowing to north, in other words from Istanbul Strait to the Black Sea shelf area. Division into smaller branches and decreasing channel depth after its sharp bend might be related with the decreasing energy of its erosion. The distal part (mid-to outer shelf) and lateral portions of the proximal part are subject to minor erosion and lower velocities as indicated by the fine-grained sediment deposition (core KP1; Figure 8 & 10).

Evolution mechanism
The crucial point of this shallow fan system is its formation in respect to sea level; in other words, whether it has formed subaerially by the abrupt flooding (or catastrophic flooding) of Mediterranean waters [34,39,71], or subaqueously below the Black Sea outflow [5,23,27]. If we consider the volume of the fan itself, it can be compared in magnitude to the volume eroded from the entire strait of the Bosporus which is empty of sediments today [8,13]. Ryan et al. [8] proposed that when the Bosphorus sill was breached by the Mediterranean waters, the resulting outburst of saltwater, gaining energy when enlarging the inlet was responsible of the onset of this depositional fan, thus fed from the Strait of Istanbul. Since this system was developed above the both erosional surfaces α and α1, the origin of the younger surface; α1 signifies the effective processes on its formation. The origin of the α1 surface was shown to have formed by subaerial exposure [23], but later revised to a subaqueous erosion in Early Holocene by Hiscott et al. [27]. They indicated a continuous sedimentation during the last 9.3 14C ky BP without a hiatus at α1 from a core collected at a depth of -69m, from ~70km west of the Bosphorus. This erosion was either due to the entrance of the Mediterranean inflow into the Black Sea, or beneath an early and vigorous Rim Current. The evolution of the channel and its associated sedimentary bedforms were shown to have occurred under entirely subaqueous conditions by Flood et al. [5], following the interpretation of α1 surface by Hiscott et al. [27]. However, they showed local erosional truncation at α1 in the outer shelf and a hiatus between 10 and 8.2 14C ky BP (from their Figure 12) at a water depth of 91m. This observation is concordant with the dunes on the desiccated north-western Black Sea outer shelf at about 80-100m water depth which had formed between 10 and 8.5 14C ky BP [15,16]. The occurrence of buried, small incised valleys and preservation of these sand dunes covered by a marine mud drape characterise a rapid transgression leaving no time to ravinement processes related to the rise of the water level to erode significantly the seafloor. This last pre-Holocene lake level fall was estimated at -100m by Lericolais et al. [15,[30][31][32] on the NW shelf, implying a rapid transgression, whereas -40m by Giosan et al. [24] off Danube delta and at -70m by Hiscott et al. [27] on the SW shelf. Parallel bedding above α surface [5,23,26] indicate a more gradual transgression after the described low stand. Numerous conical mounds are located mainly in the west side of the channel at the outer shelf. Most of them were described as mud volcanoes [5,7],

International Journal of Environmental Sciences & Natural Resources
but can also correspond to sediment waves similar to the ones that are located perpendicular to the bifurcated N-NE branch. Apart from these morphological features, small hollows cut to a depth of 90-95m are observed at -85m water depth. These hollows are expressions of scouring over exposed basement ridge as indicated by Flood et al. [5].

Stratigraphy of the south-western shelf
The almost flat widespread shelf-crossing erosional surface here-before mentioned ( §2.3, surface "α"), interpreted as sub aerial exposure of the shelf, is also clearly recognized on the profiles of this study. Four to five seismic units are truncated by this unconformity (Figure 5). These units can be correlated with the geological formations exposed on land. Unit 5 (Figures 5 &  7) is comparable with the Upper Cretaceous-Eocene volcanics exposed onland of both sides of the northern part of the Istanbul Strait [43]. Units 4, 3 and 2 can be correlated with Oligocene-Miocene aged formations on land. Paleochannels incising Unit 3 indicate that the area was subjected to erosion at that time. The presence of a huge Messinian Canyon in the southwest shelf area also supports this correlation [72]. Unit 2 observed only at the shelf edge might represent the deposition of sediments at the sea during this erosional period.
As described here and in accordance with Aksu et al. [23], Unit 1 is separated from Unit 2 by the α surface and is composed of the two subunits, 1b and 1a separated by the second discontinuity surface, α1. However, in this study, our observation from the inner part of the shelf and at the strait outlet is slightly different than observations of Hiscott et al. [27] & Flood et al. [5] and complementary in some way: a) Unit 1b between the two erosional surfaces α, and α1 becomes thin and disappears towards the inner shelf ( Figures 5 & 7), b) α1 intersects α surface beneath levees in the proximal part of the shallow fan system (Figures 7 & 9), and c) The extension of Unit 1b represents a wedge deposit at the shelf edge and is interpreted as a lowstand wedge ( Figure 5).
Unit 1b is the extension of parallel beddings above α surface in Hiscott et al. [27] and Flood et al. [5]. If subaqueous conditions were prevailing in the shelf area before 9 000y BP, Unit 1b should have been also observed in the inner shelf as well. However, this study demonstrates that Unit 1b is absent in the inner shelf. In this case, α1 surface can be interpreted either as a non-depositional discontinuity or as an erosional surface. Non depositional surfaces record episodes of reduced sedimentation, synchronous with a decrease in sand supply and/or storm wave energy [73] and can be formed; a) With no relative change in water depth and can be accounted for by either controls on decreasing sediment supply or storm-wave regime; b) With a minor increase in water depth during which the sand supply and storm-wave energy regime remained uniform, c) With a minor landward shift in facies driven by a sealevel rise [74,75]. The presence and morphology (shape) of the shallow fan demonstrates that there was an important sediment supply, and hence α1 surface must have been formed during a minor transgression (case 3). On the other hand, considering α1 as an erosional surface, absence of Unit 1b may be related to the strength of erosion which was much more effective at the outlet and inner shelf, but leaving the parallel strata in the distal parts. This change could be caused by the changing density contrast between the Black Sea and Mediterranean waters as suggested by Flood et al. [5] under subaqueous conditions. At the beginning of the inflow, the density difference is greater relative to today and hence increases the inflow velocity, creating a strong erosion, but reduces velocities as the differences becomes less. In other words, the flow sourced from the Bosphorus was initially strong and gradually reduced its strength.
There is an alternative evolution mechanism that should be considered to explain the changing erosion strength by time; changing base level or sea level [15]. This can be traced from merging of α and α1 surfaces at shallower depth than -70m, which signifies the subaerial part of the shelf area at the beginning of this erosion. Mediterranean waters entering into the Black Sea shelf along the Bosphorus channel were able to erode sediments along its route in the inner and mid-shelf and down to the substratum locally. Such strong erosion is evident from the floor of the channel at the outlet which is devoid of sediment infill ( Figure 9) and/or overlaid by coarse grained (sand-rich) sediments. The erosion down to basement rock was also shown by Gökaşan et al. [10] from a seismic profile located very close to the strait outlet (in their Figure 11d), although without direct dating. The levees along the channel belonging to Unit 1b were deposited later above α surface in the proximal part and blocked the abandoned N-NE trending branch. At a sea level of -70m [27], most of the levees of this shallow fan in the inner shelf would have been located above sea level (Figure 11), hence could have been deposited during a flood from the Bosphorus. The flooding saline Mediterranean waters penetration had an effect on sediment transport on the shoreface, generating density-driven shoreward currents. These currents may cause a significant amount of sediment transport beyond the surf zone in water depths at least up to -20m, allowing dispersion of sediment to offshore. Their mostly chaotic internal structure ( Figures 5, 7 & 10) might be related to their formation at subaerial conditions under high density driven shoreward currents. Partly inclined structures at their tops must have been formed after their immediate submergence. This is an alternative explanation for the initially strong but gradually reducing strength of flow by submergence of the shelf area.

International Journal of Environmental Sciences & Natural Resources
Lowstand wedge deposits of Unit 1b must have been deposited during the pre-Holocene lowstand and the subsequent initial rise of sea level as observed elsewhere on the shelf edges of the Black Sea [30,31,34], however dating could not be achieved from this deposit. From a long piston core collected from SW shelf edge of the Black Sea, radiocarbon date of sediments deposited after a marine transgression above a lowstand sediment wedge yielded 8.1 14C ky BP [75]. According to previous studies on the Black Sea shelf, transgression started after around ~8.5 14C ky BP [30][31][32]34,36,38,39,71]. The timing of the first entrance of the Mediterranean waters ranges from 8.4 to 8.2 14C ky BP [34,71,72,76]. The formation of channel and levee system were dated ~8-7.5 14C ky BP in Flood et al. [5] at -91m water depth, implying the sea level of the Black Sea was at the same level with Mediterranean by that time, as the inflow of water was nearly instantaneous lasting no longer than 40 years and possibly as little as a decade according to hydraulic calculations [39]. This must be correct for that part of the fan at the outer shelf, but should be slightly earlier in the inner and mid shelf. The mollusc fauna displays a mixed character of brackish and juvenile euryhaline species and was deposited until ~7.8 14C ky BP (considering the core KP12 does not reach the bottom of this sedimentary layer). Brackish mollusc forms characterize the latest lacustrine stage of the Black Sea (Neoeuxinic period), hence mixed fauna represents the influence of Mediterranean waters in this environment. The entrance of the Mediterranean waters must be earlier than 7.8 14C ky BP, due to time necessary for salinity threshold for adaptation of euryhaline species.
Unit 1b, displaying parallel stratification below α1 (between -92 and -98m water depths) at the outer shelf [5] might be an evidence of the deposition of the sediments at the submerged part of the shelf that deeper than -70m. In the NW shelf, Giosan et al. [24] indicated a sea level of -40m during the Danube was building a delta lobe before the marine reconnection of the Black Sea, by examining a core to 42m depth. However, Lericolais et al. [31] pointed out that the bottom of their core corresponds to a proximal part of the forced regression sequence which have started to deposit after 11 14 C ky BP and argued that the location of this core might have been from the upstream in the delta plain. Even at a sea level of -40m, banks at the both sides of the channel would have been subaerially exposed ( Figure 11). Following the submergence of the shelf as the sea level rises, continuous Mediterranean undercurrent processes operate on the final form of the fan and sediment distribution. The absence of Unit 1a in the distal part of the channel fan (Chirp Profile 17 on Figure 10) at the outer shelf and exposure of surface α1 indicate continuing Mediterranean undercurrent and non-deposition along its present route. Asymmetrical sediment waves are possibly formed by overspill of this flow. Coarse to fine grained sediment deposition from proximal to distal areas indicate the reducing flow velocities and also splitting into channel network.
The benthic foraminiferal fauna within the sandy sediments consists of various Mediterranean immigrant species transported with sands from the Bosphorus and only some of them could accommodate to this environment after ~6 14 C ky BP. However increasing population of benthic foraminifers and dominance of euryhaline mollusc fauna occurs after ~6 14 C ky BP, indicating a continuous flow and also relatively increasing salinity conditions at the bottom water of the shelf area. Decreasing diversity of euryhaline benthic foraminifer fauna from base to upper part of the core indicates that most of the Mediterranean immigrant species could not tolerate to adaptation or not reach to high population. Uncalibrated radiocarbon dates from Dreissena and Mytilus shells are 7.8 14 C ky BP and 6.2 14 C ky BP, from intervals 214-211cm and 130 128cm, respectively (Calibrated ages in Table 2). Most of the Mediterranean immigrant benthic foraminifer species diminish and/or are not found at ~6 ky BP. Following this change, euryhaline (‰ 1-26) A. tepida and A. parkinsoniana become more abundant first, and polyhaline (‰18-26; [78]) A. compacta prevails later at ~3.9ky BP (by assuming a constant sedimentation rate and extrapolating from dates in Table 2). These three species are the most common benthic foraminifers occurring in the recent sediments of the Black Sea [77][78][79]. The colonization order of these species indicates a relative increase in salinity after the transition to marine conditions as inferred from mollusc fauna. However, benthic foraminiferal abundance decreases at the uppermost 50cm and is represented by dominantly Ammonia genus.
This changing salinity conditions are in agreement with ecozones based on the calcareous nannoplankton [80] and benthic faunal assemblages [79]. Recently, the sedimentary sequence discovered at archaeological excavations in ancient Byzantium Harbour (Istanbul) at the northern coastal plain of the Marmara Sea indicated that the sea level was much lower than -6m, when Neolithic culture settled in the area at 8.4 14 C ky BP and lived at least a millennia [81], until inundation of this coastal area at about 6 14 C ky BP. This final maximum rise to present-day sea level occurring simultaneously at the Black Sea might have benefited this fauna.
Longer coring, particularly at the proximal part of the fan and at the lowstand wedge is necessary to resolve and to validate the ages and natures of Unit 1 and α surfaces. Spatial distribution of seismic unit thicknesses by mapping would also help for a conclusive interpretation of the process of reconnection between the Marmara Sea and the Black Sea [82].
Additional output of this study is the recognition of the banks and linear bedforms at the eastern and western shelf area of the shallow fan system (Figure 11). Banks along at both side of the channel at 20 to 45m water depth and linear ridge-shaped features are found to be closely related with the substratum formations rising to the sea floor ( Figure 11). It can be suggested that the occurrence of these formations at the eastern side of the channel might have caused the channel to redirect sharply its route from NE to NW. This should also be considered as a research topic for future investigations.

International Journal of Environmental Sciences & Natural Resources
Conclusion A shallow fan system, which consists of a channel and its short distributary branches, levees and sediment waves, is defined at the confluence of the Bosphorus and the Black Sea shelf, based on high resolution bathymetrical and seismic data. The volume of the materials constituting the fan is comparable in magnitude to the volume of the Bosporus Strait empty of sediments if we consider the material of the fan was excavated from the floor of the Strait. This fan is built by an aggrading and avulsing channel which flows from the outlet of the Bosphorus into the shelf area. It has developed above the major erosional surface α and less marked one, α1. Unit 1b is a seismic unit between two erosional surfaces, representing a lowstand wedge at the shelf edge, whereas it is absent after thinning in the inner shelf. Depositional features belonging to this fan system are defined as Unit 1a. The levees are characterized by mostly chaotic seismic reflectors and are thicker deposits accumulated in proximal areas along the main channel. As the channel becomes shallower and branches towards north, these depositional features appear to be cut by the active channel and its short branches (distributaries). The main channel is a very young formation as the recent sediments have not filled it, and it has erosive evidences of Mediterranean undercurrent due to present oceanography. The banks at both sides of the strait outlet and linear bedforms in the western and eastern shelf area of the channel are associated with the outcropping formation from the substratum.
The evolution mechanism for this fan system is shown to have initiated after a relatively strong Mediterranean-originated northerly flow into the partly subaerially exposed shelf area in the Black Sea. Mollusc and benthic foraminiferal assemblages indicated a transition from fresh/brackish to marine conditions earlier than 7.8 14 C ky BP (8 300 -8 000cal. a. BP) and continuous marine conditions since ~6 14 C ky BP (~7 000Cal a. BP).