Sedimentary Characteristics and Basin Evolution of a Compartmentalized Foreland Basin—Internal Ionian Zone, Western Greece

: This study is based on the detailed sedimentological analysis of eleven sections and one well through the late Eocene–Oligocene flysch formation of the Internal Ionian Zone (IIZ) in Western Greece. The sections are spread from the northern parts of Epirus to the north and Aitolokarnania to the south. Sedimentological data combined with biostratigraphic analyses resulted in a five-stage evolutionary model for the basin. Unit I corresponds to the lower part of the examined sections, indicating the onset of clastic sedimentation. Regarding depositional environments, it is regarded as a basin plain where lobe distal fringe accumulations occur. Unit II consists almost exclusively of heterolithic facies, marking the advance of a lobe complex system. Massive sandstone facies dominate unit III and can be considered a more proximal submarine fan system. Unit IV reflects a calm period of the basin, where mud-dominated heterolithics and hemipelagic mudstones were deposited. Hemipelagic mudstone facies with intervals of heterolithics, conglomerates, and deformed and massive sandstone facies characterize unit V. The architecture resembles a slope system incised by canyons and channels. The sand-rich intervals in Units III and V could act as the most favorable reservoir levels. In contrast, the sand-rich intervals in Unit II are considered less promising due to their higher heterogeneity.


Introduction
In foreland basins, gravity flow deposits are crucial in determining sedimentation and depositional environments.These flows are responsible for depositing large volumes of sediment and shaping the architecture of submarine fans.It is worth noting that such basin fills are of significant economic importance, as they are known to host some of the first submarine fan architectural schemes has resulted in the proposal of new models that are applicable at different scales and geological settings [13][14][15][16][17][18].With the advancement of seismic and well data quality, in combination with high-resolution outcrop data, there has been a significant improvement in our understanding of the architecture of submarine fans as well as of other depositional environments [19][20][21].The analysis of these data has enabled us to unravel the complex processes that control the formation of the genetic elements of submarine fans, leading to a more precise and detailed understanding of their characteristics.Recent studies on submarine channels and lobe complexes have provided new insights into the processes by which sediments are spread via gravity flow.New depositional models have given us further information not only about the dispersal of submarine fans but also about the evolution of basins [22][23][24][25][26][27][28][29][30][31][32][33][34][35][36] In Greece, submarine fan deposits occur mainly in the central and western continental areas.From a sedimentological point of view, flysch deposits were first assessed as submarine fan deposits by Piper [37].Other researchers supported this approach by studying several exposures of flysch in the area [38][39][40][41][42][43][44][45][46].The characterization of the study area as a foreland was introduced by Underhill [47].The evolution of the basin is mainly controlled by crustal thickening of the lithosphere due to the advancement of the External Hellenides [45] in an NW direction [47,48].Two rotational stages, as evidenced by palaeomagnetic data [49], contribute to today's geological and tectonic structure of Greece.

Materials and Methods
The scope of the present study is to assess the paleogeographic evolution of the middle Eocene to the late Oligocene foreland basin fill based on a detailed sedimentological analysis of the exposed submarine fan deposits and to perform a gross evaluation of the potential reservoir units.Detailed sedimentological identification of deep-water facies associations was carried out at 200 locations in 11 cross-sections across the basin (Figure 1).Data acquisition included recognition of basic lithological types (sandstones, mudstones, conglomerates), bed thickness measurements, and identification of thickening and thinning trends.Structural measurements (bedding planes strike and dip, fold axes, fault planes) completed the dataset.The measured thickness of the section presented in the study was corrected in areas of extensive folding (e.g., Metsovo Section).Additionally, Gamma Ray logs and cutting descriptions from Ayios Georgios 3 well [50] were re-examined, and the results are displayed along with the studied cross sections.The well is located near Dafnoti village (Figure 1), and for the convenience of this study, it will be mentioned from now on as the Dafnoti Section.Ayios Georgios well is one of the few available subsurface data sources in the area.
The sedimentological interpretation was carried out by identifying the genetic lithological types on a bed scale.These lithological types were then grouped into discrete facies that were used as a key to recognize facies associations that corresponded to genetic elements of submarine fan deposits.Paleoflow patterns of submarine fans and dating data from previous research in the area [51,52] completed the available data set for the analysis, providing the information to identify facies associations and their distribution and potential reservoir units.

Geological Setting
The study area (Figure 1) is part of the Ionian zone in Western Greece.The Ionian zone has been extensively studied, as it has attracted interest for hydrocarbon exploration for more than half a century.Structural, sedimentological, geochemical and paleontological studies [50,[53][54][55][56][57][58][59][60][61][62][63][64][65] cover almost the whole extent of the exposed stratigraphy to evaluate the hydrocarbon potential of the Ionian zone not only in Western Greece but also in Albania and similar settings in Italy.The Ionian zone consists of sedimentary rocks ranging from Triassic evaporites to Jurassic-Upper Eocene carbonates and minor cherts and shales (Figure 2).[66] regarding the tectonostratigraphic and salt tectonics processes, suggested by various researchers [66][67][68][69][70][71].The red rectangle corresponds to the measured stratigraphy.

Figure 2. Simplified chronostratigraphic correlation of available interpretations for the Ionian zone (modified from Vakalas
Thick middle Eocene [41] to Miocene flysch deposits [72] overlie the aforementioned deposits.In a regional context, the Ionian zone is one member of a series of zones (Paxos, Ionian, and Gavrovo zones) named the External Hellenides belonging to the Hellenic Thrust and Fold Belt, which is part of the Alpine Thrust and Fold Belt.The External Hellenides mainly consist of late Paleozoic-Cenozoic sedimentary rocks deposited in a series of platforms and basins located at the margins of the Tethys Ocean [73].The Africa-Eurasia convergence in late Mesozoic/Cenozoic times, as it was expressed in the collision of the Apulia plate with Europe, resulted in the inversion of Mesozoic basins forming a series of thrust sheets (External Hellenides) [71], resulting in the formation of a foreland (Pindos Foreland).The Pindos foreland is a Paleogene basin fill trending parallel to the External Hellenides and was formed in front of the Pindos Thrust [47,74] ((Figure 1).Submarine fan accumulation resulted from the deformation and uplift of the External Hellenides, which migrated westwards.During this migration, the Gavrovo and Ionian zones acted as a foreland basin [40,47,48,75].Therefore, Gavrovo and Internal Ionian zone deposits are considered a uniform genetic system [76].The age of Pindos' foreland sediments is still a matter of discussion.B.P. [77] proposed an early Miocene to middle Miocene age, explaining the presence of Oligocene fauna as a product of large-scale erosion and reworking of older sediments during the Miocene.The IGSR and IFP [78] suggested a late Eocene to early Miocene age for the basin fill, while Fleury [79], Leigh [80], Wilpshaar [81] and Bellas [82] assigned it an Oligocene age.Vakalas [52] based on calcareous nannofossils originating from 11 lithostratigraphic sections across the Internal Ionian zone, proposes a middle Eocene to late Oligocene age.Avramidis [72] proposed a middle Eocene to early Miocene age, using nannofossil zones from three studied cross sections in the Klematia-Paramythia basin (middle Ionian zone).Triantaphyllou [83] proposed a late Eocene to late Oligocene age for the Internal Ionian zone flysch deposits.

Sedimentology
The measured sections were assembled by categorizing recurring depositional facies to propose a scheme for identifying sedimentary facies correlations.

Lithofacies
The first step was to identify distinct lithofacies to rationalize the lithostratigraphic architecture described in the previous section and construct a depositional and evolutionary model for reservoir identification.For this purpose, three major lithofacies and eight sub-lithofacies have been distinguished and documented in Table 1 and illustrated in Figures 3-6.The flow type interpretation follows Talling's [84] approach.

Mudstones
Mudstones have been distinguished in three categories: massive, laminated, and hemipelagic.Observation scale limitations make difficult further recognition of graded or ungraded mudstones in the manner proposed by Piper [37] and Talling [84].
Massive mudstone (mM): Structureless, bluish to grayish mudstone packets up to 20 m, which in places might be intercalated with intervals of very thin bedded siltstones and fine-grained sandstones (Figure 3b,d).Foraminifera and trace fossils are rare or absent.Interpretation: The structureless nature and large thickness of such lithofacies resemble homogenites, as have been described by Tripsanas [85] in the NW Gulf of Mexico, Campos et al. [86] from the Gulf of Corinth, and Mulder et al. [87] from a Cretaceous deep-sea fan sequence in the Western Pyrenees.Such lithofacies are common in foreland and confined depositional settings.However, it is impossible based on the quality of the outcrops to distinguish whether these deposits are indeed homogenites or a thick succession of distal, muddy turbidites [37,88] Laminated mudstone (lM): Laminated mudstone beds occur either on top of sandstone beds or as intercalated intervals within massive mudstone intervals (mM).The laminations are typically less than 2 mm thick, forming layers less than 10 cm thick (Figure 3c,e).More compacted and better-preserved beds of mudstone reflect higher silt content.Silty beds with planar lamination probably correspond to interval Td, as introduced by Piper [37] (Figure 3e).In these beds, very fine sand may also be present.Abundant organic material of terrestrial origin oriented parallel to the lamination is also common (Figure 3f).Interpretation: According to Stow & Bowen [89,90] and Tripsanas [91], such lithofacies are attributed to deposition by low-density, fine-grained turbidity currents, which represent deposition from the most diluted parts and fine-grained parts of turbidity currents.

Hemipelagic mudstone (hM):
Usually light grayish in color, structureless, rich in microfossils and nannoplankton [52].The light color of hemipelagic mudstone (Figure 3a) reflects the increased content of calcareous material (calcareous organisms such as foraminifera and coccoliths).Common intercalations of sandstone and mudstone beds of gravity flow origin are observed.Interpretation: Such lithofacies are interpreted as the depositional product of hemipelagic/pelagic sedimentation by settling suspended oceanic material.

Sandstones
Sandstones have been categorized into massive, planar laminated, and cross laminated.

Massive sandstones (mS):
This bed type is characterized by beds of >0.6 m thick finegrained to coarse-grained sandstone.Typically, beds are structureless, with thin divisions at bed tops that grade to very fine-grained sandstone.Subtle grain-size breaks or thin discontinuity surfaces (Figure 4f) highlight amalgamation surfaces between beds.Mud clasts are common in mS lithofacies, randomly distributed or along distinct horizons.Flakes of plant material are also commonly observed throughout the mS lithofacies.Interpretation: Such beds are interpreted to be the depositional product in a layer-by-layer fashion by high-density turbidity currents (Ta) (Figure 4b,c).However, in some mud-clast rich intervals (e.g., Figure 4g) an en masse-fashion deposition by (liquefied) debris flows cannot be excluded [84,[92][93][94][95][96][97].

Planar-laminated sandstones/siltstones (plS).
In most cases, planar-laminated sandstones characterize fine-grained beds.Planar-laminated sandstones may occur as individual beds or at the top of massive sandstones or the base of cross-laminated sandstones (Tb2 and Tb1 intervals in Talling [84] (Figures 4b,d,e and 5c,d).The average bed thickness ranges from 0.2 to 1 m.Rarely spaced planar-laminated sand (Figure 4a) occurs with laminar bands ranging in thickness from 0.5 to 10 cm.In this case, the bed thickness may be up to 1.5 m.Each lamina is characterized by an upward trend followed by an abrupt increase in grain size as the next laminar interval is deposited.Interpretation: It is generally difficult to define the type of flow involved in forming planar-lamination exactly.In general, low-density and high-density turbidity currents can produce these sedimentary structures.

Cross-Laminated sandstones/siltstones (clS):
This lithological type is mainly represented by thinner sandstone beds (up to 0.5 m), which are characterized by low-angle lamination and/or climbing current-ripple cross-lamination corresponding to the Tc interval as it was introduced by Bouma [98] (Figure 5a,b,f).These beds occur either at the tops of sandstone (mS or lS) beds or interbedded as distinct beds in mudstone lithofacies.Convolute lamination is also common (Figure 5c,d).Concerning grain size, cross-laminated sandstones are mainly fine-grained.Interpretation: There is a strong consensus based on numerous flume experiments [84] that cross-laminated sandstones have been deposited by dilute and fully turbulent flows, with relatively low rates of sediment fallout.Orthoconglomerates (oCg): Stratified, well-sorted, clast-supported conglomerates of a monomictic or polymictic character.The stratified layers are up to 0.2 m (Figure 6c).Imbrication is also common.Interpretation: Orthoconglomerates are probably the result of deposition in a submarine deltaic environment by debris or grain flow processes [99].

Facies Associations
The lithological types described above are organized in recurring facies associations.These facies associations are recognized across the measured sections and in the available well logs.
CSF: calcareous shaly facies.Fine-grained, light bluish to grayish marly shales (Figure 7-3A) consisting mainly of hemipelagic mudstones (hM) interbedded with thin structureless mudstones (mM) and limestone beds (<0.2m thick).Thin beds with discrete planar-laminated (plS) and cross-laminated sandstones (clS) are common, suggesting deposition by low-density turbidites at the more distal parts of the basin.Hemipelagic mudstones have a characteristic light color, indicative of increased carbonate content.Calcareous mudstone facies occur at the base of clastic succession, formed in a transitional environment where carbonate sedimentation gradually ceases, and an increased influx of clastic material occurs (Figure 7-3B).Diluted turbidity currents transport them and have a total thickness of 30 m.

MHF: Mudstone-Dominated Heterolithic facies.
Alternations of bluish to grayish massive mudstones (mM) with very thin to thin-bedded (up to 2-3 cm) siltstones or very fine-grained sandstones (Figure 7-4A,B).Rarely, hM lithofacies are recognized.Siltstones are characterized by discrete planar lamination (plS), while sandstones are characterized by rippled tops (clS).The bases of siltstone and sandstone beds are sharp and present a good lateral extent.The sandstone/siltstone to mudstone ratio ranges from 1/2 to 2/1 (sandstone bed thickness < 0.3 m) forming packets that may exceed 100 m in thickness.
These deposits usually occur at the distal fringes of turbiditic lobes or on channel overbank areas.
HF: Heterolithic facies.Heterolithic facies consist of intercalations of sandstone beds (lS, and in a lesser degree, mS) and mudstones (mM, lM, and in a lesser degree, hM).These facies form a monotonous trend that can build units several meters thick (in some cases, more than 200 m).Sandstone beds generally present a thickness that ranges from 0.02-0.03m to 0.3 m.Rarely, beds up to 0.5 m may occur.Thicker sandstones are characterized by a massive basinal layer that gradually passes to a laminated or rippled top.The base of the beds is generally sharp.Flute and groove marks and trace fossils of the Nereites type are common.The sandstone-mudstone ratio ranges from 1:1 ( SHF: Sandstone-Dominated Heterolithic facies.Sandstone-dominated heterolithic facies consist predominantly of sandstone beds (<1 m thick) (mS, lS), interbedded with thin (commonly <0.2 m) beds of mM or lM.The sandstone-mudstone ratio ranges from 3:1 (Figure 7-2A) to 9:1 (Figure 7-2B).The thicker sandstone beds consist mainly of mS, which in most cases results from amalgamation processes as indicated by characteristic amalgamation surfaces.Thick mS gradually or abruptly passes upwards to lS beds and/or mudstone lithofacies.Scour marks (flutes and grooves) are common features at the base of the sandstone beds, and they are more prominent in erosive beds.Such sedimentary characteristics are characteristic of (1) the axis and off-axis parts of turbiditic lobes and (2) inner levees and proximal channel overbank areas.Sandstone-dominated heterolithic facies may form successions up to 200 m thick.
MS: Massive sandstone bed facies.Massive sandstone bed facies (Figure 8-3) consist of packets 50 to 100 m thick that mainly comprise wide (commonly > 1 m thick) mS beds.At the top of the beds, planar-laminated sandstones occur (plS and clS).The base of the beds is either erosive, characterized by poor lateral extent (Figure 8-3A), or sharp and can be traced at distances of more than 50 m (Figure 8-3C).Internally, the beds are rich in organic material and mud clasts.Grain size varies from fine-grained to coarse-grained sandstone.Pebbly sandstone lags may be present (oCg lithofacies).The presence of massive sandstones (Figure 8-3B) that gradually pass to a laminated top probably represents top-cut-out Bouma sequences and is perhaps the result of successive events that eroded the upper parts of the beds and transported finer sediment to the more distal parts of the basin.Massive sandstone beds represent deposition in a high-energy environment, such as channel fairways and channel-lobe complexes.Various gravity flow events, like highdensity turbidity currents or grain flow processes, may involve sediment transport.DF: Deformed facies.Deformed facies consist of blocks of contorted, folded, and deformed sandstone and mudstone facies of various lithological types (mS, lS, mM) (Figure 8-4a-c), with evidence of brecciation.DF facies may form 1-2 m up to 30 m thick intervals.In most cases, the blocks float in a muddy or silty matrix.These facies are recognized as separate units in HMF facies and rarely in HF ones.Similar facies have been discussed by Tripsanas [100] and are interpreted to represent landslide and debris-flow deposits.
HMF: Hemipelagic mudstone facies.Structureless, bluish to grayish mudstones consisting mainly of massive mM and hemipelagic (hM) mudstone lithofacies.Intervals of very thin to thin-bedded siltstones (up to 0.05 m) may occur, characterized by rippled top beds (clS lithofacies).HMF facies build very thick successions (more than 300 m) (Figure 8-1a) and are common in the upper stratigraphic levels of the basin.
PMF: Mud-matrix-supported pebbly mudstone facies.PMF facies consist almost exclusively of matrix-supported pebbly to cobble mudstones (pCg lithofacies) forming intervals that range in thickness from 4 to 20 m (Figure 8-2A-C).The matrix is composed mainly of grayish muddy material with a small percentage of very fine-grained sand.The average size of clasts ranges from 0.01m to <0.25 m.Clasts are subrounded to subangular, and their composition varies, ranging from limestone deriving from the calcareous basement to cherts and sandstones originating from the deformed Pindos zone rocks [101].
The poor organization of the deposits reflects deposition by debris flow events.The subrounded subangular shape of the pebbles indicates previous reworking in a terrestrial environment before entering the marine system.
Sand-matrix and clast-supported conglomerate facies (MCF) are organized into three types: MOC: Sand-matrix-supported organized conglomerates facies.Matrix-supported organized conglomerate facies are composed mainly of alternations of pCg (granular to pebbly concerning clast size) and with coarse-grained mS beds that, at the top, gradually pass to plS beds (Figure 9-1A,C).The conglomerates are polymictic, composed mainly of fragments of cherts, sandstones, limestones, and ophiolites.The clasts are supported by a sandy matrix and, in some cases, might be imbricated (Figure 9-1B).Sandstone beds have a thickness that ranges from 0.1 to 0.75 m.MOC facies can be recognized as 1-2 m intervals or form packets up to 10 m.Such deposits are interpreted to originate from high-density turbidity currents and grain flows.They most commonly occur in proximal channel-lobe and channel infill settings.9-2), pointing to the high energy of the environment and its erosive character that caused the scratching of the soft sediments of the seafloor and their enclosure in the load of the downslope-moving flow.Fluctuations in the flow energy are expressed by sandwiched massive coarse-grained sandstone beds (mS) observed in the dominant conglomeratic mass.The matrix is sandy.Erosive bases characterize internally the facies and correspond to successive high-energy flow events.This type of facies has been recognized in the northern and southern margins of the basin and corresponds to submarine canyon/valley settings.Depositional processes probably involve debris flows, traction carpets, and high-energy turbidity currents.
CsWoC: Well-organized clast-supported conglomerates.Well-organized clast-supported conglomerates are characterized by inclined downlapping planes consisting mainly of oCg lithological type (Figure 9-3A).The thickness of these conglomerates is up to 50 m (Figure 9-3B).The clasts are well rounded and platy with well-expressed imbricated texture, showing a size which ranges between 0.02 m and 0.05 m.The matrix is sandy.Clasts consist mainly of limestone fragments and, to a lesser degree, of cherts and sandstones.Small-scale erosive channels of poorly sorted matrix-supported conglomerates have been observed, showing a channelized geometry.Massive sandstone beds (mS) 15 cm thick are also present at the top of the facies.These facies probably correspond to a delta slope environment where delta foresets are developing [102,103].
Depositional processes probably involve debris flows, grain flows, traction carpets, and high-energy turbidity currents.
A summary of the proposed facies scheme is presented in Table 2.
Table 2. Summary of the applied facies scheme.

Mud-Dominated Heterolithic facies (MHF)
Alternations of bluish to grayish massive mudstones (mM) with very thin to thin-bedded (up to 2-3 cm) siltstones or very fine-grained sandstones.Thickness up to 100 m mM, plS, clS, hM

Heterolithic facies (HF)
Heterolithic facies consist of sandstones (mainly laminated sandstones (lS) and, to a lesser degree, massive sandstones (mS) that are interbedded with massive (mM) and laminated mudstones (lM) and to a lesser degree by hemipelagic mudstones (hM).Thickness up to 200 m plS, clS, mS, mM, lM, hM

Massive sandstone beds facies (MS)
Massive sandstone bed facies consist of packets 50 to 100 m thick formed mainly by massive sandstone beds (mS) 10-100 cm thick.Rarely, at the top of the beds, planar-laminated sandstones occur (plS) in stratification.Pebbly sandstone lags may be present (oCg lithofacies).

Deformed facies (DF)
Deformed facies consist of blocks of contorted, folded, deformed sandstone and mudstone facies of various lithological types (mS, lS, mM), with evidence of brecciation.Thickness up to 30 m.

Depositional Settings
The basin in-fill pattern comprises five different depositional settings in the study area.These settings are identified based on the participation percentage and combination of the intercalated facies associations.

Basinal Deposits
Basinal deposits consist exclusively of hemipelagic mudstone (HMF) and mudstonedominated heterolithic facies (MHF), reflecting fluctuations in the influx of clastic material from the emerging foreland and proximity to the sediment entrance points in the basin floor.Such deposits appear to overlie Eocene limestones immediately, are commonly based on calcareous mudstone facies association (CSF), and can be up to 200 m thick.Basinal facies are also observed in multiple other stratigraphic intervals of the studied sections.However, they are characterized by the absence of CSF.

Lobe Complexes
Lobe complexes are characterized by variable facies, forming intervals ranging from 50 to 200 m thick.Vertical trends considering sandstone/mudstone ratio, bed thickness variations, grain size, bed geometry, and lateral extent are used to identify lobe complex components and depositional environment.In this effort, the scheme suggested by Prélat [104] has been applied, dividing lobe complexes into the following components: lobe axis, lobe off-axis, lobe fringe, and distal lobe fringe deposits.This fourfold division has been applied to several outcrop studies [105][106][107][108][109].
Lobe axis: Lobe axis is dominated by massive sandstone facies (MS) with rare intervals of sand-dominated heterolithic facies (SHF).Thickening and coarsening-up sandy cycles with a thickness that ranges from 1.5 to 5 m are typical.There is more than 90% involvement of sandstone beds in this setting.Mudstone of SHF facies is recognized mainly at the lower intervals of each lobe.Individual sandstone bed thickness ranges from 0.1 to 0.5 m, commonly occurring on top of each other through multiple zones of amalgamation occurring parallel to the bedding plane.Sandstone bases of both MS and SHF facies are sharp with the extended presence of flute and groove marks.Trace fossils (Nereites, palaiodictyon, thallasinoeides, Spyroraphe) are also present.Lobe axis packets form up to 10 m thick intervals.
Off-axis lobes: Off-axis lobes consist mainly of SHF facies that gradually pass to MS facies close to the lobe axis or HF facies in more distal areas.Planar-and cross-laminated sandstones of SHF facies dominate with laminated mudstone intervals.Thickening and coarsening upward cycles can still be observed, and they are 2 to 4 m thick.The contribution of the sandstone beds is lower in this setting, ranging from 50 to 90%.Off-axis lobes form packets ranging from 5 to 10 m thick.
Lobe fringe: Lobe fringe deposits comprise a range of facies depending on their location on the lobe axis.HF facies are common in the more proximal locations, consisting mainly of thin-to medium-bedded planar-and cross-laminated sandstones.To a lesser degree, massive sandstones are interbedded with laminated and massive mudstones.Bedding plane bases are sharp and planar.The participation of sandstone beds in such a setting is less than 50%.In more distal locations, the sandstone-mudstone ratio is reducing, with values ranging from 1:2 to 1:3.Lobe fringe deposits are organized into monotonously repeated packets ranging in thickness from 0.5 to 1 m.The packets may be successive or separated by intervals of HMF.
Lobe distal fringe: The lobe distal fringe environment is dominated by MHF facies where thin-bedded, and very fine sandstones and siltstones are intercalated with thinbedded mudstones.Planar-laminated and cross-laminated sandstones are typical.In many cases, hemipelagic mudstones are recognized.Sandstone/siltstone presence is reduced to less than 20%.MHF facies of lobe distal fringe environments can form thick accumulations that exceed 150 m in thickness.Lobe distal fringe deposits gradually pass to hemipelagic facies (HMF) in more distal locations.

Channels and Canyons
Channels and canyons comprise a variety of facies depending on the geometrical properties and the flow types involved in the depositional processes [110][111][112].
Canyon/channel complexes: Submarine fan canyons and channel complexes crosscut slope to base of slope deposits and are characterized by dominant to rare heterolithic facies associations with intercalations of sand-matrix-and mud-matrix-supported conglomerates (PMF, CsPoC), massive sandstone (MS), and sandstone-dominated heterolithic (SHF) facies associations (e.g., upper part of Metsovo BB' section).Strongly erosional basal surfaces characterize the coarser facies associations.Thinning and fining upward cycles of MS facies associations, which might be based on conglomeratic facies, are common within such depositional settings.Amalgamated, massive sandstone (MS) packages might be 100 m thick.Deformed facies (DF) might also be present and are interpreted as collapses of the banks of the canyons/channels.The dominance or absence of HF associations might be viewed as an indication of the bypass nature of the flows [113,114] and thus provide a criterion for the proximity of the fairways to the source area (canyons vs. channel complexes).
Channel margins, levees and overbank deposits: Channel margins, levees and overbanks consist mainly of HF facies and can form packages several meters thick.SHF and MHF facies also occur, reflecting the relevant distance from the channel axis.The sedimentological properties of the aforementioned genetic elements of channel architecture are similar to those described for lobe fringe deposits.

Slope Deposits
Slope deposits form a thick succession (up to 1000 m) at the upper stratigraphic levels of the measured sections.Slope deposits consist predominantly of hemipelagic mudstone facies (HMF) and mudstone-dominated heterolithic (MHF) facies associations, which might be intercalated with HF and SHF intervals and thin to thick PMF and DF packages, suggesting the proximity of steeply inclined slopes prone to sediment failures.

Shelf-Edge Deltas
Shelf-edge deltas were identified at the upper stratigraphic levels of the northern sections of the study area.Deltas consist mainly of well-formed foresets (15-50 cm thick) dominated by well-organized clast-supported conglomerate facies (CsWoC) that pass gradually at each cycle to massive sandstone facies (15-25 cm thick) (MS).The total thickness of shelf-edge deltas has been estimated to 50 m.The clasts are well rounded and platy with well-expressed imbricated texture, showing a size which ranges between 2 and 5 cm.The matrix is sandy.Small-scale erosive channels of poorly sorted matrix-supported conglomerates (CsPoC facies) probably indicate the transition to a canyon-channel setting.

Biostratigraphic Framework
The biostratigraphic correlations presented in the next paragraph were based on the work of Vakalas [41] and Vakalas et al. [52], which cover the measured sections.The nannofossil marker species found in the samples were classified using the biozones proposed by Martini [115].The results of the biostratigraphic analysis are presented in Figure 10 and Table 3.The age determination of the geological cross sections was accomplished using three criteria: (1) Accept the presence of reworked fauna in most of the samples due to fine-to coarse-grained turbidites.Thus, it is accepted that the analyzed samples can be younger than the biostratigraphic analysis indicates but not older.(2) Lithostratigraphic criteria have been applied in samples with large age ranges and poor fauna to narrow the age determination.(3) In some rare occasions, lithostratigraphic criteria override biostratigraphic analysis but only assume that the samples can be younger but not older.
According to that scheme, the studied cross-sections span from the Upper Eocene to the Upper Oligocene.The Metsovo, Amphilochia, and Hellinika sections capture the entire stratigraphic interval.The Petrovouni, Palaiopyrgos, and Kavasila sections only capture the lower part of the stratigraphic interval, ranging from the Upper Eocene to Mid Oligocene.The age determination of the Dafnoti section was entirely based on its lithostratigraphic criteria due to the absence of any biostratigraphic analysis, which also yields an age range from Upper Eocene to Mid Oligocene.The lithostratigraphy in the thick (~4 km) Pramanta section displays various depositional settings, ranging from basinal to lobe to channel to slope settings.However, biostratigraphic analysis from all samples from this section points to an Upper Eocene age for the entire interval.Such great thickness and facies distribution discrepancy compared to the other cross sections suggests a significant uncertainty in the results of the biostratigraphic analysis, which could be attributed to reworked fauna.Therefore, the age determination of the Pramanta section was based almost entirely on lithostratigraphic criteria, yielding ages ranging from Upper Eocene to Upper Oligocene.The rest of the geological cross-section only captures the upper part of the stratigraphic interval of the basin, ranging from Mid to Upper Oligocene).

Facies Association Distribution and Stratigraphic Correlations
Based on the facies associations scheme described above and the available biostratigraphic data, the measured stratigraphy was distinguished into five units (Figures 11 and 12).

Unit I
Unit I represents the lower part of the examined sections located stratigraphically in contact with the underlying Eocene limestones.Unit I comprises the base calcareous shaly facies (CSF), which gradually pass upwards to mudstone-dominated heterolithic facies (MHF) and have been interpreted to represent basin plain deposits with a slightly variable thickness across the studied area.In the Metsovo and Pramanta sections and the Dafnoti well, where calcareous shaly facies are well-exposed, an estimated thickness of the transitional facies (CSF) is 30 to 50 m.In Petrovouni, the thickness is reduced to less than 10 m, while in the southern parts of the basin, the thickness is restricted to 10-20 m.In the Hellinika section, CSF facies are intercalated with the mudstone-dominated heterolithics.The overlying MHF facies with an average thickness of 50-100 m are interpreted to represent the distal part of submarine fans.In the Agnanta, Petas, Kompoti, and Empesos sections, Unit I cannot be recognized since the contact of Eocene limestones with the studied submarine fans succession is covered by recent alluvials.Considering the prograding pattern of the system, Unit I represents basin plain deposits overlapped by distal lobe fringe accumulations.

Unit II
Unit II consists almost exclusively of heterolithic facies presenting the range of the introduced end members of the specific facies (mudstone-dominated heterolithics to sandstone-rich heterolithics).In the Metsovo area, Unit II comprises a total thickness of 810 m and represents a higher energy environment, where coarser sandy material is deposited in the basin.Two sub-units have been recognized in this part of the section.The first unit consists of thin-bedded heterolithic facies with sheet-like geometry of the sandstone beds and good lateral extent.The second sub-unit is characterized by SHF facies organized in repeated thickening-coarsening upwards of packets of 1.5 to 2.5 m.Even though the specific sub-unit is tectonically deformed, the geometrical properties of the beds are identifiable (sharp bases and good lateral extent).In the Petrovouni section, an estimated thickness of 850 m is expected for Unit II.The basal to middle part is not well exposed.SHF facies and HF facies are dominant in two discrete sub-units.The top of the unit is marked by 100 m of mud-dominated heterolithics.The Pramanta and Dafnoti sections depict a similar stratigraphic thickness for Unit II (800 m).In Pramanta, two sub-units of SHF facies and two of HF facies have been distinguished in the stratigraphic record, presenting similar geometrical and sedimentological properties with the SHF and HF sub-units of Petrovouni and Metsovo.Occurrences of MS facies differentiate Pramanta and Dafnoti from the northern sections.At the central part of the basin in the Amphilochia section, UNIT II is very restricted (up to 70 m) and dominated by hemipelagic facies.In the Hellinika section, Unit II comprises a total stratigraphic thickness of 600 m.It consists entirely of HF facies forming cycles of thickening and coarsening upward cycles.The geometrical properties of the sandstone beds (sheet-like geometry) resemble those described in the northern sections.In terms of depositional settings, Unit II comprises elements of a lobe complex system.MHF and HF facies correspond to lobe fringe and lobe distal fringe deposits, whereas SHF and MS facies represent off-axis and axial lobe elements.

Unit III
The main attribute of Unit III is the dominance of MS facies.In the Metsovo section, MS, SHF, HF, and sand-matrix-supported organized conglomerate (MOC) facies resembling geometrical properties of a channel system compose a stratigraphic succession of 800 m.Three distinct sub-units have been identified in this part of the section.The first sub-unit consists of SHF facies, organized in thinning and fining-upward sandstonemudstone cycles.The second sub-unit consists of a matrix-supported fining-upward conglomerate, alternating with coarse-grained sandstone beds.The third sub-unit consists of HF facies.In the Petrovouni section, Unit III covers 700 m of stratigraphic thickness with MS facies characterized by erosive bases and poor lateral extent alternating with heterolithic facies.Deformed facies (DF) and matrix-supported fining-upward conglomerates (MOC) punctuate the increase in system energy.Hemipelagic facies intercalated with deformed facies dominate the upper parts of the unit.On the other hand, the Pramanta section is constructed by similar facies to the Metsovo section, dominated by MS, SHF, and MOC facies beds.Its thickness is slightly thinner (450 m).In the Dafnoti section, the setting is similar to Pramanta (unit thickness up to 500 m).It is uncertain how much to the south this system extends because unit III is not exposed at the central sections (Agnanta, Petas, and Kompoti).In the Amphilochia and Hellinika sections, unit III covers 550 and 450 m of the exposed stratigraphy, respectively.The main difference of Unit III in the Amphilochia and Hellinika areas, compared to the northern area, is the dominance of heterolithics and, to a lesser degree, the participation of MS facies.In the north part of the studied area (Metsovo to Dafnoti sections), unit III represents more proximal parts of submarine fan systems, especially the facies involved in channel fairway architecture.Confined MS facies correspond to channels, and the various heterolithic facies represent channel margin, levee, and overbank deposits.Deformed facies mark the steep geometry of channel flanks.MOC facies probably indicate channel lag deposits in the thalweg of channel fairways.On the other hand, to the south (the Amphilochia and Hellinika sections), Unit III consists of sedimentary facies that are more consistent with lobe complexes.

Unit IV
Unit IV, in the whole extent of the basin, is dominated by mudstone-dominated heterolithics (MHF) and hemipelagic mudstones (HMF).Rarely, HF and thin SHF facies occur in the studied section, and they are interpreted as either channel-overbank deposits in a slope or base of slope environment or distal fan facies in a coarse-grained sedimentstarved depositional setting.Unit IV has a stratigraphic thickness of 600 m in the Metsovo area, while Pramanta has up to 1200 m.The thickness of the Amphilochia and Hellinika sections is restricted to 500 and 400 m, respectively.In the rest of the measured sections (Petrovouni, Agnanta, Petas, Kompoti, Empesos) and Dafnoti well, UNIT IV is only partially exposed, presenting the same characteristics and facies described previously.

Unit V
Unit V is dominated by hemipelagic mudstone facies with intervals of heterolithics, conglomerates, deformed facies, and massive sandstone facies.Unit V exposes 1700 m of stratigraphic thickness in the Metsovo area organized in two discrete sub-units.The lower sub-unit (700 m thick) consists of mudstone-dominated heterolithic facies with intervals of massive sandstones and matrix-supported conglomerates (MOC facies).MS and MOC facies are accompanied by SHF or HF facies.Deformed facies (DF) and pebbly mudstone facies (PMF) are present at various levels of the stratigraphic extent.Hemipelagic mudstone facies with intervals of MOC facies dominate the second sub-unit.The upper stratigraphic levels are dominated by hemipelagic mudstones with significant packets of poorly (CsPoC) and well-organized clast-supported conglomerate facies (CsWoC).The architecture of the stratigraphy in the Metsovo section resembles a prograding slope system where the lower sub-unit facies associations correspond to a lower-slope environment.In contrast, the second sub-unit suggests deposition in an upper-slope to shelf-edge setting according to the presence of the delta-fan CsWoC facies [102,103].In the Pramanta section, Unit V presents a stratigraphic thickness of 1200 m assembled in two sub-units similar to the Metsovo section.The lower sub-unit (400 m) is characterized by mudstonedominated heterolithic facies with rare intervals (a few tens up to 100 m thick) of MS and MOC facies.Deformed facies are also common, resembling a lower-slope environment.Hemipelagic mudstone facies with rare intervals of DF facies dominate the upper subunit.Concerning the central section measured in the basin, Unit V in the Agnanta area presents a total stratigraphic thickness of more than 1800 m.Two sub-units, similar to those in the Metsovo and Pramanta areas, compose the specific stratigraphic interval.The lower sub-unit (1100 m thick) is dominated by the whole spectrum of heterolithic facies (MHF, HF, and SHF) with some intervals of massive sandstone facies.Intercalations of deformed facies are common in this subunit.The second sub-unit (700 m) consists almost exclusively of hemipelagic mudstones.An extremely deformed zone (800 m thick) overlaps the upper sub-unit, corresponding probably to a tectonically superimposed segment of the lower sub-unit.In the Petas section, Unit V comprises HF facies intercalated with MS, PMF, and DF facies.Deformed facies and massive sandstones dominate the lower parts of the unit.In Kompoti, Unit V has similar characteristics to those presented in the Pramanta section.Two sub-units can be proposed: (a) The lower sub-unit consists of sanddominated HF facies alternating with mud-dominated HF facies.MS facies accompanied by MOC facies are present at the upper parts of the sub-unit.The exposed thickness of the lower sub-unit was estimated to be 750 m.(b) The upper sub-unit is dominated by hemipelagic facies (1200 m thick), indicating either a diminishing coarse-grained sediment input in the area or an intense bypass with no deposition.In Figure 12 (the Petas-Agnanta A-A; C-C', D-D,' and E-E'), a series of cross sections showing the lower interval of Unit V is presented moving from the west to the east in the Petas-Agnanta area.The sections are separated by eastward-dipping thrusts within the flysch.Therefore, it is a safe assumption that these sections present a west-east distribution of the depositional system of the foreland basin.Interestingly, sandstone-dominated facies are only observed in the western and eastern parts of the basin, whereas mudstone-dominated facies characterize its central part.This distribution of the sandstone-dominated facies association suggests the concomitant existence of a dual channel-fairway system in the foreland basin.In the Amphilochia section, Unit V (2000 m thick) is also assembled by two sub-units with similar properties to those presented in the Kompoti section.The main difference is the presence of a significant interval (250 m) of deformed facies (DF) at the upper parts of the lower sub-unit.Mudstone-dominated heterolithic facies also dominate the top of the upper subunit.In the Empesos section, the lower parts of Unit V are exposed, characterized by sandrich HF, MS, and PMF facies.In Hellinika, Unit V is characterized by a thick upper subunit (800 m thick) consisting mainly of hemipelagic facies.In comparison, the lower subunit is 200 m thick and is dominated by massive sandstones with MOC and PMF facies intervals.According to the facies association distribution described for Unit V in terms of depositional environments, it can be inferred that the specific unit represents deposition in a slope-to-shelf-edge setting [116].Similar facies distributions have been described to occur at the slopes and basin of the active rift basin of the western Corinth Gulf during the late Quaternary [117][118][119].

Paleocurrent Data
Paleocurrent data (flute and groove marks) [51] was used to further analyze sediment distribution in the basin.(Figure 13).The study of the paleocurrent data provided the following concluding remarks about the sediment distribution in the basin: At the lower stratigraphic positions of the Metsovo A-A', Petrovouni, Pramanta, Amphilochia, and Hellinika sections, the paleoflow direction is mainly towards the west.This fact, in combination with the kind of environment recognized in these stratigraphic levels (lobe complex setting), proves that the basin was not yet confined so that the flows could move unconstrained towards the west.In the Pramanta, Petas, Kompoti, and Hellinika sections, in the middle stratigraphic levels (unit III), the paleoflow regime changes parallel to the direction of the basin axis.Considering that inner fan deposits have been recognized at these levels, it seems that the geometry of the basin must have been changing to a more trough-like character, where the confined architecture of the basin controlled the flows.The eastern margin of the trough is expressed by the Pindos thrust, and the western by the Internal Ionian thrust.In the Amphilochia section, although the flow direction is parallel to the basin axis in the western part of the section the flow shows a northeast trend, while at the eastern end, the flow trend is towards the northwest.This could be explained by different submarine fan systems in this area.At the upper stratigraphic levels of the Metsovo B-B' section, the flow pattern is orthogonal to the Pindos thrust, probably indicating some lateral sediment sources originating from the thrust.Especially in the Metsovo section, where coarsegrained deposits have been observed, such a source point seems to be genetically related to the Kastaniotikos line [120], which, in the late Mesozoic and early Cenozoic, acted as a transform fault.The thrust activity would have increased the instability in these areas, triggering submarine failures.

Basin Evolution
A five-stage basin evolutionary model is proposed based on facies associations, paleoflow patterns, and biostratigraphic correlations.Although the lower two units are captured in the Palaiopyrgos and Kavasila sections to the north, their considerable distance from the Metsovo section does not allow any direct correlation (Figures 11 and 12).

Unit I (NP16-18/Middle-Late Eocene):
The first stage represents the initiation of clastic sedimentation at the northern part of the basin (Palaiopyrgos, Metsovo, and Pramanta areas).Calcareous shaly facies (CSF) gradually pass to mud-dominated heterolithic facies in a basin-plain setting, declaring the arrival of clastic material from the progradation of the Fold and Thrust system in the east.The basin-plain deposits have a thickness of not more than 100 m, which is common in active continental margins [121].Interestingly, the thickness of Unit I and the transitional zone is variable and nearly absent in the cross sections of Petrovouni and Amphilochia (Figure 14).This is a clue that the foreland basin configuration was complex, consisting of several sub-basins.Strike slips and inherited Mesozoic basin architecture could be responsible for the complex configuration of the foreland basin.Strike zones are common figures in fold and thrust settings [122,123].The intercalations of mudstone-dominated heterolithics in the basal CSF transitional zone at the Hellinika cross section (Figures 14 and 15), combined with its much larger thickness, indicate the proximity of this location to lobe complexes.This contradicts the more basinal facies at the central and north part of the study area, suggesting differences in the depositional settings along the foreland basin through the very initial stages of its formation (Figure 15a).
Unit II (NP18-20/late Eocene): At this stage, the advance of submarine fan deposits in the northern part of the basin is evidenced by sand-rich heterolithic facies corresponding to lobe fringe, off-axis lobes, and axial lobes.The paleoflow pattern points towards the west, indicating an unconfined, widespread foreland basin (Figure 13).Lobe axis deposits are most common in the cross sections of Pramanta and Dafnoti.In contrast, off-axis and fringe lobe deposits are more dominant to the north (Figure 15).It is striking that the offaxis and fringe lobe deposits at the Hellinika cross-section are separated from the northern depocenter by the basinal deposits in the Amphilochia cross-section.This is interpreted by the potential existence of a basinal high, which separated two major depocenters in the foreland basin (Figure 14).Gravity flows always follow topographic depressions, which leads to the domination of hemipelagic sedimentation on structural highs.Unit II sandrich heterolithics have also been described by Botziolis [38] southwards of the Amphilochia section.
Unit III (NP21-22/Early Oligocene): Channel sandstones and conglomerates are dominant in the northern of the foreland basin, extending from the Metsovo to Dafnoti sections (Figures 14 and 15).However, more mixed facies are observed in the Petrovouni section, where thinner sandstone and conglomerate packages and deformed facies are intercalated within hemipelagic sedimentary facies.Such facies associations are typical in slope settings.This indicates that at this time interval, the northern depocenter becomes dissected into two depocenters, one to the north expressed by the Metsovo section and one to the south defined by the Pramanta and Dafnoti sections (Figure 14).In a similar manner, the lobe complexes characterize the Amphilochia section to the south.In contrast, off-axis and fringe lobe deposits describe the Hellinika section, suggesting the presence of two additional depocenters to the south.The breaking of the initial unconfined foreland basin into multiple depocenters is attributed to the development of strike slips through the westward advancing fold and thrust system (Figure 15).The most likely candidates are the Kastaniotikos and Alevradas strike slips (Figure 15) Paleoflow trends (Figure 13) during this time interval are oriented parallel to the elongated basin axis (NNW-SSE).The differentiation in paleoflow trends from units I and II to Unit III is attributed to the evolution of the initial unconfined foreland basin into a confined basin, with a well-expressed trough geometry with the major axis trending NW-SE.This confinement is probably due to the initiation of the Gavrovo Thrust at the eastern part of the study area.The tectonic deformation of the upper part of Unit II in the Metsovo section supports such a conclusion.The only exception in the paleoflow trends is observed in the Metsovo section, where a WSW flow direction is observed.Such a flow orientation is consistent with the work of [72] suggesting an eastward sediment source for early Oligocene flysch formation in the Botsara syncline, Internal Ionian zone, to the west of the Metsovo section.

Unit IV (NP22-23/Early-Late Oligocene):
The early to late Oligocene stage marks the prevalence of fine-grained deposits (Figures 14 and 15d) consisting mainly of hemipelagic facies with sand-starving heterolithic facies (MHF to HF) which can be related to two factors: (a) The advancement of Gavrovo Thrust, which might have resulted in the formation of piggy-back basins to the east and, thus, in the entrapment of coarse-grained material in these depositional settings.(b) There was a significant regression event during the early to late Oligocene, expressed by a sea-level rise of 150 m [124].Abundant deformed facies beds in the Petrovouni section indicate that, at least in this area, the gradient of slopes increased, promoting sediment failures in the nearby slopes.This is consistent with the implied Kastaniotikos strike-slip going through this region during the early Oligocene (unit III).

Unit V stage (NP23-25/Late Oligocene):
The common characteristic of this unit throughout its thickness in all cross sections is the presence of rare to typical deformed and conglomeratic facies indicating proximity to steep slopes with ample slope instability events.This is the first clue of a late Oligocene placement of the Gavrovo Thrust within an incipient foreland basin.The tectonically deformed facies within the upper part of unit IV in the Metsovo and Pramanta sections support such a conclusion.Unit V consists of two upper and lower subunits (Figure 14).The large facies variability in the lower subunit among the studied cross-sections indicates the further breakage of initial foreland basins into several depocenters.Intercalations of sand-matrix conglomerates (MOC) with mudstone-dominated heterolytics, grading upwards to delta slope conglomerates (CsWoC) is consistent with a prograding slope and consistent with a feeding system to areas further to the west, as reported by the eastward sediment supply of the flysch formation in the Klimatia-Paramythia basin [72].On the other hand, a dual channel-fairway system is documented in the Petas to Agnanta sections to the western and eastern parts of the basin, which are separated by more hemipelagic facies.Paleoflow trends in both channel fairways are northwards (Figure 15e).The feeder of the western channel fairway appears to be through the Empesos area, where typical channel sandstone and conglomerates have been documented.On the other hand, the channel fairway to the east is attributed to the emplacement of the Gavrovo Thrust within the basin, generating an SSE-NNW corridor (Figure 15e).An analog to such a tectonic setting comes from the Quaternary Ionian Basin between Lefkas and Corfu islands (Figure 16).In this area, the emplacement of a Pliocene thrust in the foreland basin has resulted in its dissection into two, a N-S-oriented basin to the west and a SSE-NNW corridor to the east, which might act as a sediment fairway from the east.Delta deposits on the hanging wall of this thrust support this theory.
The mudstone-dominated facies in the lower subunit in the Pramanta section suggest the elevation of this part of the basin by the westward-advancing fold and thrust system.However, some rare and thin conglomerates and sandstone facies indicate that this area was still partially acting as a southward feeding system.
Similarly, the mudstone-dominated facies with some rare and thin sandstone intervals in the Amphilochia section suggest the elevation of this area by the westward-advancing fold and thrust system.The presence of thick and abundant deformed and mud-matrix conglomerates upwards within unit V indicates proximity to a steepening slope with amplitude of sediment failures.On the other hand, channel complex deposits (conglomerates and massive sandstones) characterize the Hellinika section.
The upper subunit of unit V comprises hemipelagic facies with only some rare and thin intercalations of heterolithics (Figure 14).This sudden change in the depositional system is attributed to the complete advancement of the fold and thrust system and the emplacement of the Internal Ionian Thrust [72].The tectonically deformed facies in Agnanta E-E' section support such a conclusion.As the tectonic activity migrates westwards and the basin is transformed to a piggyback configuration, the monotonous sheet-like heterolithics could be regarded as ponded basins (Brandford, 2003).The northern and southern parts of the basin are characterized by submarine canyons that incise slope and transport sediment westwards through pathways formed by significant strike-slip zones [120,125].Shelf-edge deltaic deposits mark the uplift and gradual transition to a shallow marine environment.As the convergence of the African and Eurasian plates continues, sedimentation migrates westwards, where seismic profiles verify the proposed evolutionary model (Figure 16).

Reservoir Evaluation
The investigation of submarine fan deposits from the middle Eocene to late Oligocene within the Internal Ionian zone has demonstrated the presence of two promising reservoir units, Unit III and Unit V.At the same time, Unit II is deemed less favorable, as explained in the following paragraphs.The diverse trapping mechanisms may include salt tectonics, compressional tectonics that facilitate the development of anticlinal structures, and the primary formation of stratigraphic traps [127].The present study focuses mainly on the primary depositional control and the geometrical relationships of the submarine fan elements.
Unit II: The expansion of sediment load from a confined canyon-channel system results in the formation of a complex lobe system.Sand-rich, sheet-like deposits describe the sedimentological profile of Unit II exposures.Sand-dominated heterolithic facies and massive sandstone beds, organized in repetitive cycles of 2.5 to 8 m, forming stacked lobe intervals of 10 to 100 m thick, could suggest a potential reservoir unit.These deposits are recognized in the northern sections (Metsovo, Petrovouni, Pramanta, and Dafnoti well) and interpreted as lobe axis to off-axis lobe elements.The lateral extent of these deposits is difficult to trace due to the physical constraints of outcrop continuity.The nature of facies associations involved in the specific unit favors the formation of a highly heterogeneous reservoir by various studies in similar settings [128].The heterogeneity is primarily controlled by two factors: in a mesoscopic scale, heterolithic facies contribution at a vertical and lateral scale, where mudstone content ranging from 20% to 80%, reduces the net to gross values.From a microscopic point of view, massive sandstone facies are affected by the involvement of fine-grained components (silts and very fine-grained sands) that increase matrix content and reduce primary permeability.The tectonic setting may also reduce the volume of permeable reservoir units, considering that reverse faults trending NW-SE could compartmentalize the sand bodies and cut off the link with the feeding channel-canyon system.Compartmentalization is also favored by the W-E lobe axis orientation (see paleocurrent analysis), which is vertical to sub-vertical to reverse fault strikes.
Unit III: Unit III has been interpreted as the most sand-dominated interval of the measured stratigraphy, corresponding to a channel and channel-lobe system architecture.Confined sandstone bodies with clast-supported conglomeratic facies form successions up to 500 m thick with high net-to-gross values.Laterally, sand-rich to mud-dominated heterolithics corresponding to channel margin, levee, and overbank deposits stratigraphically confine the sand bodies [36] MS facies of channelized deposits are expected to be "clean" [26] compared to the MS facies of Unit II.Channel architecture, where highly permeable zones (MOC and MS facies) are intercalated with heterolithics, results in a highly heterogeneous reservoir.This architecture can result in the isolation of sandstone bodies and reduce connectivity.Various studies in similar reservoirs [129][130][131] showed that facies architecture in channel systems has a negative impact on recovery.Various techniques have been proposed to enhance production [129,132].Concerning post-deposition effects on reservoir quality, especially the impact of tectonics, NW-SE-trending reverse faults may result in the confinement of channel fairways that, in many cases, have the same orientation in the study area.This, according to some researchers, may inhibit reservoir formation [133,134] while others suggest that at this configuration, sand-and gravel-rich reservoirs are likely to form [135].
Unit V: Unit V has been attributed to a slope-apron to slope setting.According to [22], the highest concentration of reservoirs along muddy slopes occurs across mid-slope and toe-of-slope regions in a variety of accommodation types, including ponded, healedslope, and incised submarine valleys, which seems to be the case in the study area.Unit V comprises two potential reservoir sub-units.The lower sub-unit resembles Unit III properties.Although sand and conglomerate percentages are higher, the steep geometry of the channels favors the development of DF facies with increased mud content that can act as flow barriers [22,24,26].The upper sub-unit consists mainly of various conglomeratic facies confined by extended hemipelagic mudstone facies.Conglomerates are deposited in the axis of canyons, cutting the slope.Facies involved in canyon architecture again introduce primary constraints in reservoir quality due to the type of matrices for both matrixand clast-supported conglomeratic types.Mud clasts of significant size reduce the overall porosity, especially at canyon thalwegs.Tectonically, the orientation of canyons is vertical to sub-vertical to the major reverse fault systems, favoring the segmentation of and the reduction in reservoir volumes.

Correlation with Mesozoic Ionian Zone Petroleum System
The Ionian Zone flysch has been regarded as part of the Ionian Zone petroleum system [67], playing the role of a sealing unit and favoring the formation of anticlinal traps.The main reservoirs in Western Greece consist of the Upper Cretaceous to Eocene carbonates, which have an estimated thickness of 300-700m [67].Eocene carbonates pass gradually to flysch.Fracture and vuggy porosity describe the main attributes of the aforementioned carbonate reservoirs in the area, as also identified by studies in the extent of the Ionian zone in Albania [136,137].The literature has underexplored the sealing capacity and integrity of the middle Eocene to Oligocene flysch of the Ionian zone.In the present study, it can be implied that the extent and thickness of Unit I deposits, where calcareous hemipelagic deposits are overlapped by mud-dominated heterolithic facies (thickness up to 200 m), act as an efficient sealing formation for the Upper Cretaceous to Eocene carbonate reservoirs.The efficiency of middle Eocene to Oligocene flysch has been evaluated in terms of a producing oil field in Albania [136].However, further analysis is required to assess the brittleness of mudstones under geostatic and tectonic stresses [138,139].
Finally, studies from the area [57] imply that middle Eocene to Oligocene flysch accumulations or its younger analogs in the middle Ionian zone [69] can be regarded as autonomous petroleum systems, where specific organic-rich mudstone intervals would play the role of source-rock horizons.The lack of sufficient subsurface data weakens verifying the hypothesis mentioned above.

Conclusions
The study analyzed sediment in 11 sections and a well in Western Greece's Internal Ionian Zone during the late Eocene-Oligocene period.
The measured sections were systematically assembled by categorizing recurring depositional facies into groups.This helped identify sedimentary facies correlations and propose a facies associations scheme for submarine fan environments.Using three basic bed types and eight subcategories, 11 facies were proposed to assemble facies associations.The proposed facies can be correlated into five distinct units.Unit I lies stratigraphically in contact with the underlying Eocene limestones and is regarded as a basin plain where lobe distal fringe accumulations occur.Unit II comprises heterolithic facies and suggests deposition in a lobe complex system.Unit III is dominated by massive sandstone facies, and Unit IV is dominated by mud-dominated heterolithics and hemipelagic mudstones.Unit V has hemipelagic mudstone facies with intervals of heterolithics, conglomerates, deformed facies, and massive sandstone facies, resembling a slope system incised by canyons and channels.
A five-stage evolutionary model is finally proposed based on facies distribution, paleocurrent data, and biostratigraphic correlations.The initial stage highlights the beginning of clastic sedimentation due to the emergence of the thrust and fold belt.The second stage is marked by the advance of lobe complex deposits succeeded by the channel fairways of the third stage.The gradual compartmentalization of the basin is recorded due to (a) internal thrusting and (b) strike-slip faulting.In the fourth stage, muddy deposits accumulate due to either the advance of mud-rich slope deposits or significant sea-level rise.The compartmentalization of the basin characterizes the fifth stage into three zones.The central part of the basin is isolated, and sediment supply is reduced as the major Pindos thrust gradually ceases.In contrast, the northern and southern parts of the basin are incised by submarine canyons transporting sediments westwards through pathways formed by major strike-slip zones.
Regarding reservoir characterization, sand-rich Units III and V are the most favorable.At the same time, Unit II is less promising due to the high heterogeneity resulting in expected lower net-to-gross values.
Funding: This research received no external funding.

Figure 1 .
Figure 1.Internal Ionian Zone (IIZ) submarine fans extent (blue color).The major structural elements of the area are also marked, as well as the location of the measured sections.

Figure 5 .
Figure 5. (a,b) Typical examples of cross-laminated sandstones (yellow arrows, Tc).(c,d) Wellformed convolute lamination on top of laminated sandstones.(e) Successive cross-laminated siltstones (Tc) on top of laminated fine-grained sandstones (Tb1) (the base of the succession is at the top of the photo).f) Cross laminated sandstones (Tc).4.1.3.Conglomerates (Cg) Granural, pebbly, and cobble extraformational conglomerates of a monomictic or polymictic character, depending on the source rocks.Based on the texture, two types of conglomerates have been distinguished: Paraconglomerates (pCg): Matrix-or clast-supported, with a chaotic texture corresponding to Talling [84] Dm2 interval (Figure 6a,b), forming packets up to 15 m thick.Interpretation: Paraconglomerates are considered debris flow deposits based on their structureless and poorly sorted texture.At the top of the aforementioned structureless paraconglomerates, normal grading suggests the transformation of debris flow into turbidity currents (Figure 6a).Orthoconglomerates (oCg): Stratified, well-sorted, clast-supported conglomerates of a monomictic or polymictic character.The stratified layers are up to 0.2 m (Figure6c).Imbrication is also common.Interpretation: Orthoconglomerates are probably the result of deposition in a submarine deltaic environment by debris or grain flow processes[99].

Figure 10 .
Figure 10.Biostratigraphic range of the measured sections of the Internal Ionian zone according to Vakalas [41,52].The light reddish part of the bars indicates the uncertainty of dating due to the presence of reworked nanofossils.In the Palaiopyrgos section, a stratigraphic gap (question mark) due to a barren sample is also marked.Concerning the Pramanta section, the type of uncertainty is described in the text.

Figure 11 .
Figure 11.Stratigraphic correlation of the measured sections of the Internal Ionian Zone (IIZ) where the proposed units are also marked.The red square marks the stratigraphic extend of Agnanta E-E' section presented in Figure 12.

Figure 13 .
Figure 13.Paleocurrent directions map.Vectors are classified based on the facies associations units.Numbers in circles refer to the measured sections.The stratigraphic level of each dataset is marked in the corresponding lithostratigraphic column.

Figure 14 .
Figure 14.A schematic cross-section in an N-S orientation illustrates the five units' proposed distribution and the various facies' lateral correlation.The dark blue pattern indicates the calcareous Eocene basement.The Pindos thrust trace is also shown.

Figure 16 .
Figure 16.A modern analog of the proposed evolutionary model.A transfer zone favoring the transportation of sediments at the deeper parts of the basin.Deltaic facies are dominant at the shallower eastwards-located part of the basin.Notice also the segmentation of the basin due to internal thrusting, the scarp produced by the westernmost thrust, and the produced mass transport complex deposits (MTC) (Seismic profile ref.: Bellas [126].The bathymetric map originates from the Emodnet database: https://portal.emodnet-bathymetry.eu/, accessed on 17 February 2016)

Table 1 .
Summary of the basic bed types and their subcategories.
Well-organized sand-matrix clast-supported conglomerates are characterized by inclined bedding planes consisting mainly of oCg lithological type.Massive sandstone beds (mS) 15 cm thick are also present at the top of the facies.Thickness up to 50 m.
). Intervals of very thin to thinly bedded siltstones (up to 5 cm) may occur, characterized by rippled top beds (clS lithofacies)