Sedimentary Facies Analysis, Reservoir Characteristics and Paleogeography Significance of the Early Jurassic to Eocene Carbonates in Epirus (Ionian Zone, Western Greece)

Sedimentological, micropalaeontological, and marine geological results from the Early Jurassic to Eocene carbonate formations of the Ionian zone, from six localities of Epirus, provide new insights into the basin palaeogeographic evolution and better correlation with coeval analogous tectono-stratigraphic successions along the southern margin of the Neo-Tethys Ocean. Facies analysis allowed the recognition of several microfacies types and their depositional characteristics. During the Early Jurassic, autochthonous carbonates (Pantokrator Limestones) were deposited in shallow-water environment. The overlying (hemi)pelagic Siniais or their lateral equivalent Louros Limestones were deposited to the basin borders and mark the general deepening of the Ionian domain. During Toarcian to Tithonian, the Ionian Basin was characterized by an internal differentiation in small sub-basins with half-graben geometry presenting abrupt thickness and facies changes. The deeper parts were characterized by continuous sedimentation, while the elevated parts were marked by unconformities. The Early Cretaceous marks the homogenization of sedimentation by the deposition of the pelagic Vigla Limestones all over the Ionian zone. The transition from the Early to Late Cretaceous records a significant carbonate diversification in terms of biota assemblages, and related mineralogy due to intense tectonic activity in the region. From Late Cretaceous to Paleogene, allochthonous carbonates were transported to the outer shelf by turbidity currents (calciturbidites) and/or debris flows (limestones with breccia) formed by the gravitational collapse of the platform margin. Additional porosity and bulk density measurements showed that petrophysical behavior of these carbonates are controlled by the depositional environment and further influenced by diagenetic processes. The partly dolomitized neritic Jurassic carbonates, but mainly the Senonian calciturbidites and the microbrecciated Paleocene/Eocene limestones display the higher average porosity values, and therefore present enhanced carbonate reservoir quality.


Introduction
Marine biogenic carbonates are among the most important archives of Earth's history reflecting past changes in ocean chemistry, water mass circulation and the evolution of life [1][2][3][4][5][6][7][8][9][10]. Particularly during the Cretaceous, changes in climate, oceanographic circulation, sea level, tectonic and volcanic activity contributed to triggering widespread crisis events, which resulted in world-wide anoxic  [34,35]) illustrating the principal tectonostratigraphic zones: Pre-Apulian, Ionian, Gavrovo, Pindos. The red box shows the study area, the northwestern part of the Epirus region, where the separation of Ionian zone to external, middle, and internal sub-basins along with the regional locations of the study sections to be also indicated. Legend interpretations are presented in the inset.  [34,35]) illustrating the principal tectonostratigraphic zones: Pre-Apulian, Ionian, Gavrovo, Pindos. The red box shows the study area, the northwestern part of the Epirus region, where the separation of Ionian zone to external, middle, and internal sub-basins along with the regional locations of the study sections to be also indicated. Legend interpretations are presented in the inset.

Lithostratigraphy of Jurassic-Eocene Formations in the Ionian Basin
All the study sections are situated within the middle-internal Ionian Zone, in Epirus region near Ioannina (Figure 2A-E). They represent to the Mesozoic-Paleogene pre-to post-rift sequences of Ionian zone, and particularly the upper part of Early Jurassic Pantokrator Limestones, the pelagic Early Cretaceous Vigla Limestones, the Late Cretaceous Senonian Limestones, and the microbreccious limestones of Paleocene/Eocene age. These formations are described in stratigraphic order as follows.

Pantokrator Limestones and Lateral Equivalents
Early Jurassic (Hettangian to Sinemurian) Pantokrator Limestones represent the upper part of the pre-rift sequence of the Ionian zone [64,72]. This formation can be described as a neritic formation, consisted of limestones and dolomites with calcareous algae and benthic foraminifera of more than 1000 m of total thickness. These shallow-water limestones overlie Early to Middle Triassic evaporites (>2000 m thick) and the Foustapidima Limestones of the Ladinian-Rhaetian. Due to their remarkable facies' homogeneity, such Jurassic marine deposits represent an extensive carbonate platform at that time covering the entire region from the Apulian to the Gavrovo zone [58]. The overlying syn-rift sequence in the deeper part of the half-graben sub-basins begins with the Pliensbachian pelagic Siniais Limestones and its lateral equivalent Louros Limestones (both indicative of the first deepening event of the Ionian basin), overlain in some sub-basins by the succession "Ammonitico Rosso", "Limestones with filaments", and "Upper Posidonia Beds", while in others, the "Ammonitico Rosso" is replaced by the coeval "Lower Posidonia Beds" in the previous sequence, or even the whole sequence is represented by the undifferentiated Posidonia Beds [34,58]. The boundary between the Pantokrator Limestones and Louros or Siniais Limestones is gradational. J. Mar. Sci. Eng. 2020, 8, x FOR PEER REVIEW 5 of 26 sequence in the deeper part of the half-graben sub-basins begins with the Pliensbachian pelagic Siniais Limestones and its lateral equivalent Louros Limestones (both indicative of the first deepening event of the Ionian basin), overlain in some sub-basins by the succession "Ammonitico Rosso", "Limestones with filaments", and "Upper Posidonia Beds", while in others, the "Ammonitico Rosso" is replaced by the coeval "Lower Posidonia Beds" in the previous sequence, or even the whole sequence is represented by the undifferentiated Posidonia Beds [34,58]. The boundary between the Pantokrator Limestones and Louros or Siniais Limestones is gradational.

Vigla Formation
This pelagic formation can be described as sub-lithographic thin-bedded to platy deposits, rich in plaktonic organisms (Calpionelles, Radiolaires, and Globotruncanes) of Cretaceous age. Their continuity is usually interrupted by thin intercalations, chert lenses, or even thin chert layers. Its dominant lithology consists of light grey to yellowish micrites and radiolarian biomicrites (wackestones to packstones). More explicitly, the calcareous beds (Vigla Limestones) of the lower part consist of mudstones-wackestones, biomicrites with foraminifera and radiolaria, and siliceous wackestones and packstones [42]. In the upper part, this formation consists of the Vigla Shale member (also known as "Upper siliceous zone"; [66]) composed by yellow marly or shaly limestones and shales with chert intercalations, and red to green or locally black clay layers. The latter contains the equivalent of the anoxic events of Selli (OAE1a) during the Aptian-Albian, Paquier Evet (OAE1b) of Early Albian age and Bonarelli (OAE2) at the Cenomanian-Turonian boundary in the Ionian zone [12,14,73,74], and extends to Italy and Albania. Vigla Limestones generally feature considerable lateral variations in thickness, mainly in the basin borders where its thickness may reach the double of its average, which is indicative of persistent differential subsidence [34,58,59]. In the internal subbasin, the Vigla Limestones consist of compact, thick-bedded, often bituminous, dolomitic limestones, with lenses of slightly dolomitized microbreccia and thin cherty intercalated layers or nodules. It is also characterized by the presence of several aptychus in its lower part [59,64]. In Epirus, along the eastern border of the central Ionian zone, phosphatic horizons are intercalated with the uppermost horizons of the Vigla Limestone formation, stratigraphically above the Vigla Shale member. Overall, this formation corresponds to the first post-rift sediments of the Ionian zone [58,59,64].

Vigla Formation
This pelagic formation can be described as sub-lithographic thin-bedded to platy deposits, rich in plaktonic organisms (Calpionelles, Radiolaires, and Globotruncanes) of Cretaceous age. Their continuity is usually interrupted by thin intercalations, chert lenses, or even thin chert layers. Its dominant lithology consists of light grey to yellowish micrites and radiolarian biomicrites (wackestones to packstones). More explicitly, the calcareous beds (Vigla Limestones) of the lower part consist of mudstones-wackestones, biomicrites with foraminifera and radiolaria, and siliceous wackestones and packstones [42]. In the upper part, this formation consists of the Vigla Shale member (also known as "Upper siliceous zone"; [66]) composed by yellow marly or shaly limestones and shales with chert intercalations, and red to green or locally black clay layers. The latter contains the equivalent of the anoxic events of Selli (OAE1a) during the Aptian-Albian, Paquier Evet (OAE1b) of Early Albian age and Bonarelli (OAE2) at the Cenomanian-Turonian boundary in the Ionian zone [12,14,73,74], and extends to Italy and Albania. Vigla Limestones generally feature considerable lateral variations in thickness, mainly in the basin borders where its thickness may reach the double of its average, which is indicative of persistent differential subsidence [34,58,59]. In the internal sub-basin, the Vigla Limestones consist of compact, thick-bedded, often bituminous, dolomitic limestones, with lenses of slightly dolomitized microbreccia and thin cherty intercalated layers or nodules. It is also characterized by the presence of several aptychus in its lower part [59,64]. In Epirus, along the eastern border of the central Ionian zone, phosphatic horizons are intercalated with the uppermost horizons of the Vigla Limestone formation, stratigraphically above the Vigla Shale member. Overall, this formation corresponds to the first post-rift sediments of the Ionian zone [58,59,64].

Senonian Limestones
This formation corresponds to thick pelagic strata (1 to 3 m) intercalated with calciturbidites, comprising limestones with fragments of globotruncanids and rudists and fine-grained breccia intervals with limestones and rudist fragments within calcareous cement containing pelagic fauna. Calciturbidites are usually alternated with sub-lithographic pelagic limestones and cherts [34,35]. Overall, this formation corresponds to a period of basinal sedimentation, with significant variations of its lithological characteristics and sedimentary facies from the external (western) to the internal (eastern) parts of the Ionian basin. In the external sub-basin, they have been described as "clastic limestones" [75] containing floatstones, rudstones and rare grainstones and packstones with a micritic matrix, biomicritic intercalations, and rare cherty nodules. In the middle sub-basin, these carbonates are characterized as microclastic, bioclastic or polygenic microbreccias with a micritic matrix (wackestone, packstones), intercalated with (bio)micrites. In the internal sub-basin, the Senonian Limestones are massive, thick bedded microbrecias to breccias containing rudists and coral fragments, with only local appearance of thin-bedded or nodular chert [65,66,76].

Paleocene-Eocene Limestones
This formation also known as "Limestones with microbreccia" [75] is characterized by the same sedimentary facies with the underlying formation, with prominent microbreccia derived from the erosion of Cretaceous carbonates from both the Gavrovo (to the east) and Apulian (to the west) platforms. On Late Paleocene and progressively through the Early Eocene, the supply of calciclastic material diminished significantly, especially in the central Ionian Basin. The main depositional facies during the Eocene were platy wackestone and mudstone with Globigerinidae and chert nodules, analogous to those of the Vigla Limestones, but lack continuous siliceous intervals [34].

Materials and Methods
Six sections (A-F; Figure     The study of the Jurassic to Eocene sediments' distribution was based on tectonic, stratigraphic, paleontologic, and sedimentologic observations. Depending mostly on outcropping conditions, an average sampling interval of approximately 1 m has been used. An overall number of 150 samples were collected, among which 70 used for sedimentary facies analysis and 50 for measuring porosity and bulk density. Microfacies definition and textural characters analysis of the carbonate rocks, including both biogenic and inorganic dominant components, were done according to carbonate classification schemes of [77], which later modified by [78], based on the Standard Microfacies Types (SMF) in the facies zones (FZ) of the rimmed carbonate platform model. Depositional environments were reconstructed based on the derived sedimentological features and through comparison with additional standard facies reconstructions [79][80][81][82][83][84][85]. Therefore, assemblages consisting of several SMF types characterize all depositional environments. Thin sections were prepared in the Historical Geology and Paleontology Laboratory (National and Kapodistrian University of Athens; NKUA), biostratigraphically and sedimentologically studied under a polarized LEICA DM LP microscope, and photos have been performed with OLYMPUS UC30 Microscope Digital Camera. The reservoir potential of these carbonates was also obtained through the examination of porosity and bulk density parameters on representative of each stratigraphic formation samples, with the use of GeoPyc 1360 Envelope Density Analyzer and helium AccuPyc 1330 Pycnometer, respectively, also at NKUA. All measurements were made on dry samples that were previously dried in the oven at 40 °C for 24 h. The petroleum potential information was further compared with existing regional data covering the studied stratigraphic interval. Finally, ArcGIS software was used to visualize the results of the study area, along with additional data known from the literature, outcropped and well data [33,35,42,59,86] expanding the entire Ionian zone within paleogeographic/paleoenvironmental maps for the relevant time interval. The study of the Jurassic to Eocene sediments' distribution was based on tectonic, stratigraphic, paleontologic, and sedimentologic observations. Depending mostly on outcropping conditions, an average sampling interval of approximately 1 m has been used. An overall number of 150 samples were collected, among which 70 used for sedimentary facies analysis and 50 for measuring porosity and bulk density. Microfacies definition and textural characters analysis of the carbonate rocks, including both biogenic and inorganic dominant components, were done according to carbonate classification schemes of [77], which later modified by [78], based on the Standard Microfacies Types (SMF) in the facies zones (FZ) of the rimmed carbonate platform model. Depositional environments were reconstructed based on the derived sedimentological features and through comparison with additional standard facies reconstructions [79][80][81][82][83][84][85]. Therefore, assemblages consisting of several SMF types characterize all depositional environments. Thin sections were prepared in the Historical Geology and Paleontology Laboratory (National and Kapodistrian University of Athens; NKUA), biostratigraphically and sedimentologically studied under a polarized LEICA DM LP microscope, and photos have been performed with OLYMPUS UC30 Microscope Digital Camera. The reservoir potential of these carbonates was also obtained through the examination of porosity and bulk density parameters on representative of each stratigraphic formation samples, with the use of GeoPyc 1360 Envelope Density Analyzer and helium AccuPyc 1330 Pycnometer, respectively, also at NKUA. All measurements were made on dry samples that were previously dried in the oven at 40 • C for 24 h. The petroleum potential information was further compared with existing regional data covering the studied stratigraphic interval. Finally, ArcGIS software was used to visualize the results of the study area, along with additional data known from the literature, outcropped and well data [33,35,42,59,86] expanding the entire Ionian zone within paleogeographic/paleoenvironmental maps for the relevant time interval. Figure 5 illustrates a general view of the studied sections, regarding the Jurassic-Eocene carbonates of the Ionian zone in the study area. Agios Georgios section consists of 50 m of a uniform, largely condensed shallow-water limestone unit (Figure 5a,b) with calcareous algae, benthic foraminifera, and, rarely, brachiopods. The middle part of this succession (samples L15-L20) is recrystallized and slightly dolomitized, a phenomenum that usually occurred during the Jurassic in the Ionian zone. Such Early Jurassic sediments represent the lower part of the carbonate unit in western Greece. The basal 10 m of Perivleptos section contains neritic Pantokrator Limestones, overlain unconformably by reddish and green organic matter-rich shales with thin-to medium-bedded marly limestone interbeds and siliceous lenses, the Vigla Shales (Figure 5c- Figure 5 illustrates a general view of the studied sections, regarding the Jurassic-Eocene carbonates of the Ionian zone in the study area. Agios Georgios section consists of 50 m of a uniform, largely condensed shallow-water limestone unit (Figure 5a,b) with calcareous algae, benthic foraminifera, and, rarely, brachiopods. The middle part of this succession (samples L15-L20) is recrystallized and slightly dolomitized, a phenomenum that usually occurred during the Jurassic in the Ionian zone. Such Early Jurassic sediments represent the lower part of the carbonate unit in western Greece. The basal 10 m of Perivleptos section contains neritic Pantokrator Limestones, overlain unconformably by reddish and green organic matter-rich shales with thin-to mediumbedded marly limestone interbeds and siliceous lenses, the Vigla Shales (Figure 5c-e).  The Vigla section, which has a thickness of 20 m, is characterized by typical Vigla formation deposits, such as pelagic limestones interbedded with centimeter-to decimeter-thick radiolarian chert beds (Figure 5f,g). Koloniati section (60 m thickness) consist of Cretaceous deposits of both Vigla and Senonian Limestones formations separated by an unconformity (Figure 5i). The lower half consists of pelagic massive limestones with rare intercalations of chert beds, and spherical, but also flattened, cherty nodules (Figure 5i,k). The siliceous nodules are beige, or light grey to dark grey color. The lithology of the upper half is described as solid, thick-bedded limestones that can be easily separated from the underlain thin-bedded ones (Figure 5j). In the uppermost part of the section brecciated horizons have been observed. Similar microbreccia bioclastic limestones with roudist fragments were observed in the 10 m thick Asprageli-2 section (Figure 5l,m). Finally, at the top of the studied carbonate succession, the 25 m thick Asprogeli-1 section consists of grey-light micro-to mesoporous, turbiditic limestones with abundant radiolaria and planktonic foraminifera (Figure 5n,o).

Sedimentary Facies Analysis
Section A (Agios Georgios Section): The principal lithofacies of this section is boundstone, composed of red algal communities, and less frequent benthic foraminifera within a micritic clotted matrix (Figure 6a,b). Some lime horizons are dolomitized.
Section B (Perivleptos section): The lower part of this section is characterized by the lithofacies recrystallized grainstone of peloids with signals of dolomitization and fracturing (Figure 6c,d), while the upper part consists of mudstone-wackstone with radiolarians and/or planktonic foraminifera (Figure 6e).
Section F (Asprageli-1 Section): In Section C a packstone with planktonic foraminifera and some scattered mollusc and benthic foraminifera (Figure 6n,o) has been identified.
Generally, high energy environments such as platform, fore-shoal, and intertidal channel display grain supported textures. On the contrary, in low energy environments such as deep shelf and open marine, mud supported textures developed. In between, medium energy settings include platform slope environments, deposits of which are characterized by significant debris of rudists, algae, porcelaneous benthic foraminifera, and rarely variable shell fragments (bivalves, bryozoans). Grainy or micritic texture also strongly affected porosity, fluid flow, and diagenetic processes.
Section C (Vigla Section): Assemblages of small-to-medium-sized morphotypes of Calpionella alpina followed by the acme of Calpionella elliptica at the basal part of the section suggest an Early Berriasian age, while progressively through the top of the section radiolarians and other planktonics (Favusella hauterivica, Hedbergella sigali, Hedbergella dendroensis) have been further identified. In the upper half of the section the Globigerinelloides "acme" and "eclipse" intervals, characterized by the abundant presence and the lack of representatives of this genus, respectively, show the evolution through the Barremian up to the Albian in the topmost meters of the Vigla Limestones.
Section F (Asprageli-1 Section): The genera Subbotina, Acarinina, and Morozovella are the dominant groups of Paleocene-Early Eocene assemblages. More explicitly, the lower part of the section contains carbonates enriched on radiolarians and planktonic foraminifera (Parasubbotina pseudobulloides, Acarinina sp., Subbotina sp., Igorina pussila, Chiloguembelina sp.), which deposited during the Early Paleocene (Selandian). Progressively through the top of the section, the increase of highly and full body ornamented Morozovellid species (e.g., Morozovella aequa and Morozovella subbotinae, Morozovella velascoensis) and quadrate Acarininids (e.g., Acarinina wilcoxensis and Acarinina pseudotopilensis) are indicative to the end of the Paleocene and the start of the Eocene (Ypresian).

Porosity and Bulk Density Measurements
The results of porosity and bulk density measurements of the studied carbonates are listed in Table 1, while the average values per formation in Table 2. With the exception of Vigla Shales, which present the lowest measured bulk density values (2.60-2.63; average 2.62 gr/cm 3 ), the bulk density is quite homogeneous for all carbonate samples, ranging from 2.64 to 2.74 gr/cm 3 , with an average of 2.68-2.70 gr/cm 3 , respectively, per study formation (Table 2). On the contrary the total porosity values evaluated with the pycnometer present significant variability between the different sections and stratigraphic formations. In Agios Georgios section porosity values range between 2.31% and 9.71%, with an average value of 4.97% (Table 1). Significant porosity values above the average value for the Jurassic deposits were recorded in some horizons of this unit (samples AG10-20), where the limestones seem to be dolomitized. These values approximate the maximum porosities (~10%) that have been reported for this formation for the entire Ionian zone [33]. The total average porosity for the Pantokrator Limestones, including both Agios Georgios and the lower part of Perivleptos sections, is 5.26% (Table 2).  Table 1). On the contrary, the typical Vigla Limestones in both Vigla section and the lower part of Koloniati outcrop display the lowest values ranging between 1.81% and 3.86% (average 2.86%) in the former and 1.28% and 5.65% (average 3.28%) in the latter (Table 1). Overall, the average porosity for Vigla Limestones, taking into account the data from both aforementioned areas, equals to 3.02% ( Table 2). The analyzed samples from the upper part of Koloniati section characterized as Senonian Limestones showed steadily high porosity in the range of 3.86% to 9.26% (average 6.58%), especially at the top of this unit, where brecciated limestones with the maximum porosities of about 10% were identified ( Table 1). The same range showing the tendency for increased porosity values (2.60-9.63%, average 6.44%; Table 1) was observed in Asprageli-2 section, with the hemipelagic calciturbidites and fractured resedimented microbreccia. The total average porosity measured for this Late Cretaceous carbonate unit is 6.50% ( Table 2). The porosity of Asprageli-1 brecciated Paleogene limestones is in the range of 2.50-9.35% (Table 1), with the average of 5.63% (Table 2).

Microfacies Types and Depositional Environments
The main textural and compositional characteristics, as well as the sedimentary features of the distinguished microfacies, are summarized in Table 3, corresponding to different depositional environments or facies zones (FZ) defined by [84,85] (Figure 7). More specifically, the Pantokrator Limestones were classified as boundstone of algae (SMF 7) and bioclastic grainstone (SMF 11) in the Agios Georgios and Perivleptos sections respectively, giving evidence of a depositional environment characterized by a platform, with both intertidal and subtidal environments (FZ5-6). Recrystallization-dolomitation and fracturing of some redeposited carbonate clasts points a subaerial exposure of parts of the platform. Vigla Limestones consist of mudstone-wackestone with radiolarians and planktonic foraminifera (SMF 2-3), characterizing a low energy, relatively deep environment, such as the toe of slope (FZ3) and/or deep shelf (FZ2). Senonian Limestones present a variety of lithofacies, which correspond to environments ranging from the slope to deep shelf. They mostly include in-situ wackestone-packstone with planktonic foraminifera along with microbreccia bioclastic packstone with fragments of shallow water fauna (rudists and benthic foraminifera) (SMF 4), and micrite with transported ooids (SMF 13), which in total represent a medium-to-high-energy environment (e.g., slope), possibly due to the transportation of the sediments within the basin from the platform.
The above facies distribution reflects the separation of the deep Ionian Basin into a central topographically-higher area characterized by reduced sedimentation, and two surrounding talus slopes with increased sedimentation [66]. Locally, micritic peloids and dark-gray intraclasts, floating in an overall micritic matrix were also observed within these carbonates. The range of depositional interpretations of these formations includes supratidal settings, vadose-marine inorganic precipitation in inter-and subtidal environments of formation in marine seepage or groundwater springs [85]. Overall, Late Cretaceous calciturbidites suggest relatively deep-slope depositional conditions (FZ3-4). Relatively similar and possibly deeper depositional conditions apply for the Paleogene biomicritic packstones with radiolaria (mostly at the Early Paleocene) and planktonic foraminifera (through the Eocene) (SMF3), suggesting a basinal environment (FZ2).  environment, such as the toe of slope (FZ3) and/or deep shelf (FZ2). Senonian Limestones present a variety of lithofacies, which correspond to environments ranging from the slope to deep shelf. They mostly include in-situ wackestone-packstone with planktonic foraminifera along with microbreccia bioclastic packstone with fragments of shallow water fauna (rudists and benthic foraminifera) (SMF 4), and micrite with transported ooids (SMF 13), which in total represent a medium-to-high-energy environment (e.g., slope), possibly due to the transportation of the sediments within the basin from the platform. The above facies distribution reflects the separation of the deep Ionian Basin into a central topographically-higher area characterized by reduced sedimentation, and two surrounding talus slopes with increased sedimentation [66]. Locally, micritic peloids and dark-gray intraclasts, floating in an overall micritic matrix were also observed within these carbonates. The range of depositional interpretations of these formations includes supratidal settings, vadose-marine inorganic precipitation in inter-and subtidal environments of formation in marine seepage or groundwater springs [85]. Overall, Late Cretaceous calciturbidites suggest relatively deep-slope depositional conditions (FZ3-4). Relatively similar and possibly deeper depositional conditions apply for the Paleogene biomicritic packstones with radiolaria (mostly at the Early Paleocene) and planktonic foraminifera (through the Eocene) (SMF3), suggesting a basinal environment (FZ2).

Reservoir Potential of the Early Jurassic to Eocene Carbonate Rocks of the Ionian Zone
Reservoir characterization deals with physical characteristics of the reservoir, including petrophysics (porosity-bulk density and grain density measurements, capillary pressure measurements), fluid properties (e.g., reservoir fluid saturations) and reservoir drive mechanisms [51,91]. Moreover, reservoir quality is defined as the amount of porosity and bulk density in a reservoir and can be a function of many control factors for both carbonates and sandstones [92][93][94]. In the present study, the reservoir properties, porosity and bulk density, are comparatively examined in the studied sections on a large intra-basin scale in order to assess the quality of carbonate reservoirs of the Ionian zone. Measuring porosity and bulk density of a given reservoir is a direct measure for the storage and flow capacity. Though porosity seems to be a main contributor to the flow capacity, bulk density is mostly controlled by the pore throat distribution [95]. However, they are difficult to predict, since they depend on both initial depositional processes and diagenetic overprinting. Particularly, their complexity in carbonate reservoirs should be attributed to the different interplays, among other factors, of hydrodynamic conditions, carbonate cementation or dissolution, and tectonic setting that form the architecture of the marine setting [96].
The surveyed carbonates showed remarkable average porosities (3.02-6.50%; Table 2) accompanied by even lower bulk densities (2.62-2.70%), and therefore, the quality of the reservoir has been described as poor to fair. Low and homogeneous bulk density values are in good agreement with the literature data for pure calcite rocks [97][98][99]. Average porosity values are quite variable in the study sections, with a tendency to slightly higher values in section E and upper part of section D both of Senonian age, as well as in the upper part of section B corresponding to the Vigla Shales. On the contrary, the lowest porosity (and bulk density) values were reported for the biomicritic mudstone-wackestone depositional facies with pelagic fauna of Vigla formation. The Pantokrator Limestones present significant porosity and bulk density values only on dolomitized horizons, explaining the significant variability observed into this formation. Generally, this intra-zone original porosity evolution characterized by an increased tendency from the Early Cretaceous pelagic limestones to the Late Cretaceous calciturbidites, which are also slightly more permeable, is consistent with previous findings in the Ionian zone of NW Greece [33,100,101]. Particularly, the overall evolution in the study area display a quite variable average porosity pattern, characterized mostly by low to moderate values (~3-7%), while some high individual values around 10% also recorded. The levels with such a sudden porosity increase may be related to burial diagenesis and dolomitization, which increase the reservoir quality to good.
Factors related to sequence architecture, including the occurrence of intervals with clay laminations (e.g., Vigla Shales; promoting chemical compaction and associated cementation) and the distribution of early dolomitization (promoting porosity preservation during burial), generally increases secondary fracture porosity. A higher increase of porous within the marly limestones of Vigla Shales formation than in Senonian calciturbidites, could be caused due to fracturing due to the nodule's development, in the manner described by [43]. Intracrystalline porosity could be also developed within dolomitized intervals and therefore some breccia limestones may form potential reservoirs. Especially for the Ionian zone, it is documented that the dolomitization front has changed through time [33]. In the internal and external parts of the Ionian zone, dolomitization continued well into the Cretaceous, whilst in the central part, it did not continue after the Middle Jurassic. However, it is worth noting that the low porosity and bulk density identified within the study formations could either imply high fluid pressures or fluid migration through permeable "fracture conduits" in the vicinity of fault zones. The observed porosity values correspond to microporosity that does not take into consideration the fracture porosity related to thrusts which considerably increases the carbonate reservoir quality. Some of the porosity variations observed in this study may due to the vicinity/distance of the samples with the tectonic zones. In any case, distribution of Ionian tectonic zones is related to the prevailing tectonic style of the Ionia zone, which is a combination of thick-and thin-skinned deformation [33]. The elucidation of the predominance of tectonic style has not been achieved to date, due to the fact that deep seismic surveillance is hindered by the subsurface Ionian evaporites [34,58,70,71].
The incorporation of the depositional with the carbonate reservoir quality data indicate that there is not a clear view for the reservoir quality that can be associated with a specific depositional environment. Generally, in western Greece, most of the fair to good potential reservoirs are deposited in shallow to restricted platforms (e.g., Gavrovo platform carbonates close to its transition to the Ionian zone, in the thrust sheets of Ionian and pre-Apulian zones) [33,34,40]. In the Ionian zone, similar medium to high energy environments, such as tidal domains, reef barriers, and slopes were recorded during Jurassic (e.g., Pantokrator Limestones, Posidonia beds) and Late Cretaceous to Paleocene/Eocene (e.g., Limestones with microbreccia), respectively [35,40]. However, Jurassic studied sediments do not contain any proven reservoirs with the most significant porosities to be associated with the development of fracture and/or diagenetic zones. This study shows that deeper depositional (deep marine basinal and/or slope) environments can also be associated with the deposition of good potential reservoirs. On this regard, potential reservoir rocks within the Ionian zone further include the upper part of pelagic Vigla Limestones, Senonian Limestones, and the microbrecciated intervals of the Paleocene/Eocene limestones, all presented very good porosity values up to 10%.

Paleogeographic Analysis of the Ionian Zone
In Epirus area only the Jurassic to Eocene carbonate succession occurs. The carbonate platform sediments begin at the base with thick-bedded neritic Jurassic Pantokrator Limestones, which feature remarkable facies homogeneity, indicating that an extensive shallow sea was spread all over the study area during that time. In Perivleptos and Agios Georgios sections, this facies association is mostly represented by biolithitic boundstones and biosparite grainstones with calcareous algae and benthonic foraminifera, implying a carbonate margin platform with both intertidal and subtidal environments ( Figure 8A). These extensive platforms are developed until the Hettangian-Pliensbachian age, when the overlying synrift sequence begins. Pliensbachian Siniais Limestones correspond to the general deepening of the Ionian Basin. The structural differentiation that followed caused the fragmentation of the initial basin into smaller paleogeographic units with half-graben geometry. This is recorded in the abruptly changing thickness of the synrift formations that take the form of syn-sedimentary wedges [34]. In the deeper parts of the half grabens, these wedges include complete Toarcian-Tithonian successions, whereas in the shallower parts of the half grabens, the successions are interrupted by unconformities. However, this topmost part of the Jurassic is not recorded in the studied sections. In Koloniati section, there is an unconformity between the base of the Vigla Shales and the topmost horizons of the underlying Pantokrator formation, which marks a period of uplift and erosion at the beginning of the Toarcian. This led to occasional karstification of Pantokrator Limestones.  The post-rift sequence begins with the pelagic Vigla Limestones, whose deposition was synchronous throughout the Ionian Basin, beginning in the Early Berriasian [73,102]. The basal sequence of the Vigla limestones, consisted of thin layered, sub-lithographic, pelagic limestones, with abundant radiolarian and frequent cherty beds enriched with radiolarian, is related to the Early Cretaceous subsidence caused deepening in the entire basin ( Figure 8B). Towards the upper part of this formation, chert layers become more abundant, containing intercalations of green, red, and locally black shales, named as the "Vigla Shales" member, indicative of basinal sedimentation. The microfacies analysis of the Koloniati (lower part) and Vigla sections suggest that these carbonate sediments represent a low energy, relatively deep environment, like the toe of slope and the deep shelf of the basin. Apart from the halokinetic movements, which probably caused the variation in thickness of Vigla Limestones [34] from the western (external) to the eastern (internal) parts of the basin, the pelagic depositional conditions persisted until the Late Eocene, when flysch sedimentation began.
During the Late Cretaceous, the Senonian Limestones formation consisted of hemipelagic calciturbidites and resedimented microbreccia, reflect similar deep marine slope environments. In particular, sedimentary facies analysis of the Koloniati (upper part) and Asprageli-2 samples suggests that they were deposited in a deep-water toe of slope and platform margin environment, respectively, where miccrobrecciated carbonates were transported and accumulated ( Figure 8C). The allochthonous bioclastic material identified in the analysed samples consists mainly of rudist (typical reef builder) fragments and benthic foraminifera (e.g., forereef dweller Orbitoides; [103] and/or inner platform taxa Cuneolina, Textulariids and Miliolids; [85,103,104]), originated and transported from a nearby shallow shelf environment (e.g., platform or reef). Such shelf margins of the nearby pre-Apulian and Gavrovo platforms [75] and/or internally to the Ionian basin [101] were characterized by high productivity of such skeletal material, transported and redeposited in the deeper parts of the Ionian basin [34]. These bioclasts are also accompanied by pelagic Globotruncanids and radiolarian specimens observed in the in situ micritic matrix. Moreover, ooid lithofacies and some reworked lithoclasts observed in the upper part of the Koloniati Senonian Limestones indicate shallow-water conditions that were exposed during the uppermost interval of the Cretaceous (Maastrichtian). This is further reinforced by extensional tectonics, while possible sea level effects cannot be ruled-out and could also be related to the eustatic sea level low-stand (~150 m sea-level drop) that took place between the Late Cretaceous and the Paleocene [105,106]. Our paleoenvironmental observations from the Late Cretaceous interval in Epirus fully agree with recent sedimentary findings of [42] for the Araxos area (internal Ionian zone), as well as with the previous literature [34,35,75] for the entire Ionian zone. Overall, the facies distribution of the Senonian reflects the separation of the Ionian Basin into a central area (middle and outer part of the Ionian Zone) characterized by deeper water sedimentation and two surrounding talus slopes, issued from western Gavrovo platform and western Apulian platforms. Both platforms provided the clastic carbonate material that was transported by turbidity currents into the Ionian Basin.
The study of the Asprageli-1 samples provides evidence that the supply of clastic material due to tectonics diminished significantly during the Paleocene/Eocene. However, despite the reduced tectonic activity during that time, the slumping of platform edge sediments produced turbidity currents resulting to the deposition of Limestones with microbreccia and calciturbidites. The main depositional facies of platy mudstone-wackestone with Globigerinidae, Globorotaliidae and rare siliceous nodules, analogous to those of the Vigla Limestones, imply that the depositional environments during that period did not change significantly from the Late Cretaceous ( Figure 8D). The greatest thicknesses of the Eocene units can be found in the marginal parts of the Ionian Zone, where the microbreccias are more frequent.

Conclusions
The Ionian zone consists of a heterogeneous multi-layered calciclastic reservoir in Epirus region (western Greece). The identified carbonate formations display various facies ranging in a full spectrum of depositional conditions, from shallow platforms (reefs) to slope (platform margin) environments, even to the open marine settings, with different lithologies, sedimentary features, energy conditions, and diagenetic overprints. This heterogeneity also explains the rock petrophysical/geomechanical variation of these carbonate rocks. The Early Jurassic limestones (biolithites boundstone) do not contain any proven reservoirs due to relatively low porosities, with the exception of the microcrystallized or dolomited horizons, which increase the reservoir quality in a local scale. The Early Cretaceous limestones and cherts (biomicrites mudstone-wackestone) of Vigla formation has been described as the poorest of the studied carbonates, in terms of their reservoir potential. On the contrary, the Late Cretaceous (Senonian Limestones) and the Paleocene/Eocene carbonate units can be considered the primary target for oil/gas exploration in the study area, since they contain calciturbidites deposited mainly in the slope (bioclastic packstone-rudstone with rudist fragments and benthic foraminifera) and the deep shelf (planktonic foraminiferal biomicrites mudstone-wackstone). The highest porosity values recorded in those carbonates may be further associated with the development of fracture networks and/or diagenetic zones. Overall, the results of this study may have implication for reservoir-and/or source-rock-geologists and diagenetic modelling approaches in the presented area and elsewhere within the eastern Mediterranean Sea [42,86,100,101,107,108], implying that sample specific analyses or a very well understood regional diagenetic framework are required for accurate prediction of reservoir quality.