Bio-sequence stratigraphy of the Neogene: an example from El-Wastani gas ﬁeld, onshore Nile Delta, Egypt

. Due to modern hydrocarbon development and exploration activities throughout the onshore Nile Delta of Egypt, a high-resolution biochronologic sequence stratigraphy of the Neogene sequence was conducted to illustrate the gas-bearing reservoirs’ depositional sequences. Our study used a multidisciplinary approach comprising biostratigraphy, facies analysis, geophysical logs, and seismic data to shed light on the Neogene stratigraphic framework. The biostratigraphic analysis of planktonic foraminifera and calcareous nannofossils allowed the recognition of six zones and/or subzones and nine zones, respectively. An open-shelf environment was suggested for the Middle–Upper Miocene Sidi Salem Formation, while the Upper Miocene Qawasim and Abu Madi formations were deposited under stressed environmental conditions interpreted as estuary facies. The Lower Pliocene (Zanclean) succession deposited in the middle to outer shelf domains, including the upper-bathyal environments. Seven depositional sequences bounded by six major sequence boundaries were recognized from Serravallian to Zanclean times. These boundaries signiﬁcantly inﬂuenced changes in reservoir properties and architecture of the incised valley ﬁlls. The sea-level oscillations are interpreted by correlating the sequence boundaries and ﬂooding surfaces with global eustatic charts


Introduction
The architecture of the Nile Delta system is critical to understand the eastern Mediterranean plate tectonics.It is located along the northern margin of the African Plate, which extends from the Tethys subduction zones to the rift of the Red Sea, which separates it from the Arabian Plate.The geological study of the Nile Delta was initiated in the early 1960s, when oil corporations began drilling in the area.Due to modern exploration and use of geological and geophysical modelling approaches, this province provided some hidden hydrocarbon resources (EGPC, 1994).The Nile Delta is characterized by a thick sequence of Neogene clastic deposits, which represent one of the most significant gas provinces in Egypt (El-Heiny et al., 1990;Abdel al et al., 2000;El-Barkooky and Helal, 2002;Ismail et al., 2010;El-Kahawy et al., 2022;Selim et al., 2022).The Messinian onshore Nile Delta reservoirs were sedimentologically assessed, and their significant potential as reservoirs was confirmed (Selim et al., 2021).
The delta architecture is characterized by several palaeohighs and palaeo-lows, thereby influencing facies distribution of the Nile Delta.Correlation of the different sub-basins using biostratigraphy is a powerful tool in establishing the sequence stratigraphy framework of the Nile Delta.Our aim is to clarify the Nile Delta stratigraphic framework using wellconstrained biostratigraphic data integrated with 2D seismic profiles, wireline logs, and sequence stratigraphic analysis.We also introduced a quantitative analysis of foraminifera and calcareous nannofossils to investigate their relationships to system tracts compared to gamma rays and CycloLog patterns.

Geologic background and lithostratigraphy
The Nile Delta (Fig. 1a) is considered a significant hydrocarbon reservoir in the south-eastern Mediterranean region.The Nile Delta is composed of two significant provinces (Sestini, 1995), a deep offshore and an onshore sector.The deep Published by Copernicus Publications on behalf of The Micropalaeontological Society.offshore domain received substantial sediment input (3500 m in 4-5 Myr) during the Plio-Pleistocene.It was affected by extensive Pliocene listric faults; prominent rollovers; rotated blocks; and slumping, particularly oriented toward the northnorth-east and north-east.The onshore sector was dissected into northern and southern sub-basins to the south and to the north by a flexure zone termed the hinge line.The geometry of the pre-Miocene strata is influenced by the east-westtrending hinge line, through the mid-delta area (Said, 1981).
The southern-delta province (10 600 km 2 ) displays northward-dipping middle Eocene carbonates, indicated by asymmetric folds associated with the Syrian arc system that extends to the southern part of the Nile Delta and the Western Desert (Kamel et al., 1998).Deep pre-Tortonian faults with an east-west trend are most probably related to Eocene or Late Cretaceous tectonic events.Another fault system corresponds to post-Messinian deformations, mostly due to sediment loading along an unstable Nile Delta margin.These are normal faults also related to slumping and diapir rises in uncompacted Pliocene sediments and Messinian evaporites (Kamel et al., 1998).
From the Neogene to the Quaternary, the geological history of the Nile Delta has been subdivided into three cycles: Miocene, Plio-Pleistocene, and Holocene in age (Rizzini et al., 1978).The base of the Miocene cycle is represented by the Sidi Salem Formation, which is overlain by the Messinian Qawasim, Abu Madi, and Rosetta formations with fluviomarine facies (Fig. 1d).The second cycle is characterized by clay deposits of Plio-Pleistocene age, with the Kafr El-Sheikh, El-Wastani, and Mit Ghamr formations (Fig. 1d).The third cycle consists of the Holocene Bilqas Formation composed of clays and subordinate sandstone beds.The Sidi Salem, Qawasim, Abu Madi, and Kafr El-Sheikh rock units (from older to younger) have been investigated (Fig. 1d).

Materials and methods
A collection of 151 subsurface ditch cutting samples was obtained from five wells located in the El-Wastani field (Fig. 1a,  b, and c), spanning the Neogene successions of the Nile Delta.
The foraminiferal content was studied in 20 g of drilled cutting samples.A hydrogen peroxide solution (15 %) was used to soak the crushed samples in water.The solution was then sieved with a 63 µm mesh to filter the finest fractions.The residue was dried at low temperature (60 • C) and screened through a series of sieves (63, 125, 250, and 400 µm).The foraminifera species have been analysed using a Leica stereo-binocular microscope of the sieve-filtered fractions > 63 µm.The identification process is based on the schemes of Bolli and Saunders (1985) and Loeblich and Tappan (1988).The full taxonomic list of the planktonic and benthic foraminiferal taxa is provided in Appendix A.
To study calcareous nannofossils, samples have been prepared according to the method of Fornaciari et al. (1996).To capture the scarcer calcareous nannofossil taxa, each sediment sample was subjected to three replicates, from the bottom, middle, and top of the water column in the beaker.During microscopy, six traverses were examined; each traverse consisting of over 500 fields of view (FOVs) has been screened, and all the encountered specimens were counted after a quick slide analysis to detect the rare taxa.Relative abundance is estimated as follows: > 40 specimens in 500 FOVs is considered abundant, 30-40 is considered common, 20-30 is considered frequent, 10-20 is considered rare, and < 10 specimens is considered very rare (see File S2 in the Supplement).The taxonomic list of calcareous nannofossil taxa is presented in Appendix A.
The zonation schemes for planktonic foraminifera and calcareous nannofossils used in this study are based on Lirer et al. (2019), Agnini et al. (2017), Backman et al. (2012), Fornaciari et al. (1996), and Raffi et al. (2003), as established for the Mediterranean region.These schemes are correlated with the standard-zonation nannofossil schemes of Martini (1971) and Okada and Bukry (1980).
The scanning electron microscope (SEM) was used to capture the foraminiferal species, while we used a phasecontrast, light-polarized Nikon microscope at ×1000 magnification using a Tucsen digital camera to image calcareous nannofossil species.
The foraminiferal and calcareous nannofossil samples analysed in this study are curated at the Paleontological Museum of the Geology Department at Cairo University.
In this study, the gamma-ray log was analysed using the CycloLog software, as outlined by Nio et al. (2005).This software examines the frequency patterns of geophysical logs that are associated with facies, providing cyclic and stratigraphic approaches.By applying the CycloLog software to the gamma-ray log, two curves are generated: PEFA (prediction error filter analysis) and INPEFA (integrated prediction error filter analysis).

Planktonic foraminiferal biostratigraphy
Six foraminiferal zones and/or subzones spanning the Middle Miocene to Lower Pliocene sequence have been identified based upon the recognition of 54 species and 34 genera of both planktonic (Figs.2-3) and benthic  foraminifera.We followed the biozonal scheme for planktonic foraminifera proposed by Lourens et al. (2004), Iaccarino (2007), and Lirer et al. (2019).The Mediterranean biozones were compared with other similar zones in geographically adjacent areas (Fig. 6).The similarities and dissimilarities between different authors are also discussed.The stratigraphic ranges of the recognized species (planktonic and benthic) in the investigated wells are quantitatively presented in File 1 in the Supplement.
From bottom to top, the planktonic foraminiferal zones recorded in the studied wells are as follows.
Occurrence.This subzone is 60 m thick and is only encountered in well EWE-2 at depths ranging from 3050 to 2950 m (Fig. 7e; see File 3 in the Supplement).
Assemblage.This planktonic foraminifer subzone comprises a rich and diversified benthic assemblage (File 1 in the Supplement).
Remarks.This zone corresponds to the Serravallian.In the Atlantic, it correlates with zone N12 of Blow (1969) and with the upper parts of N12 of Srinivasan and Kennett (1983) and Berggren et al. (1995).It is equivalent to M9b of Wade et al. (2011), based on the occurrence of Paragloborotalia siakensis.In the Nile Delta Basin, it corresponds to the Globoro-talia mayeri zone of Abdou et al. (1984) and to the Globorotalia fohsi zone of Rizzini et al. (1976).

Neogloboquadrina acostaensis interval zone (MMi11)
Definition.It is defined as the interval between the LO of Neogloboquadrina acostaensis and the LO of the Globigerinoides extremus.
Occurrence.The biozone has only been encountered in wells EWW1-St2 and EW-5.It occurs at depths from 3300 to 3200 m in well EWW1-St2 and from 3300 to 3100 m in well EW-5, reaching a thickness of 100 and 200 m, respectively (Fig. 7e, d).
Assemblage.See File 1 in the Supplement.

Non-distinctive interval zone (NDZ)
Definition.According to Iaccarino and Salvatorini (1982), this zone was originally defined as an interval determined by the transition in the coiling directions of Neogloboquadrina acostaensis from sinistral to dextral at its base and the first occurrence of Globorotalia suterae at the top.However, planktonic foraminifera lack this interval in our study.Subsequently, the zone is identified here based on a few sporadic occurrences of benthic foraminifera (Ammonia beccarii and Elphidium sp.) and ostracods (Cyprideis spp.).
Occurrence.The zone is recorded in wells EWE-2, EW-8, EW-12, EW-5, and EWW1-St2 at depths of 2950 to 2630 m, 3120 to 2670 m, 3105 to 2725 m, 3100 to 2650 m, and 3200 to 2805 m, respectively (Fig. 7a, b, c, d, and e).The respective thicknesses in the five wells are thus 320, 450, 380, 450, and 395 m.Remarks.This zone corresponds to zone N17 of Blow (1969).In the Mediterranean region, the NDZ is characteristic of the Messinian stage, as described by Iaccarino and Salvatorini (1982).Many works dealing with the eastern Mediterranean (i.e.Iaccarino et al., 2007;Lirer et al., 2019) have shown the absence of planktonic foraminifers in the NDZs.In the Atlantic domain, Berggren et al. (1995) recognized the Late Messinian age by the occurrence of Globigerinoides extremus (M14) and the lowermost part of the Gt.tumida/Gt.nepenthes (PL1), according to Wade et al. (2011).The non-distinctive interval zone characterizes the Late Miocene Messinian stage in the Nile Delta Basin (Rizzini et al., 1976;Farouk et al., 2014).

Sphaeroidinellopsis seminulina interval zone (MPl1)
Definition.This interval zone comprises the LO of Sphaeroidinellopsis spp., as a common occurrence, and the LO of Globorotalia margaritae.
Occurrence.The zone has been encountered in the basal part of the Kafr El-Sheikh Formation from wells EWE-2, EW-8, EW-12, and EWW1-St2.The monitored depths in the four wells are from 2805 to 2788 m, 2670 to 2660 m, 2725 to 2717 m, and 2630 to 2615 m, respectively (Fig. 7a-c and  e).Therefore, the thickness of this zone varies as follows: 17, 15, 8, and 15 m.
Assemblage.See File 1 in the Supplement.
Remarks.Following Cita (1975), Iaccarino and Salvatorini (1982), and Iaccarino (1985), the occurrence of Sphaeroidinellopsis seminulina in the Mediterranean corresponds to zone N18 of Blow (1969) of Early Pliocene age.This zone also corresponds to the Mediterranean zone MPl1, where the top horizon has been placed according to the first frequent occurrence of Globorotalia margaritae (Lourens et al., 2004;Iaccarino et al., 2007;Lirer et al., 2019).This biozone also corresponds to the lower part of the Globorotalia tumida/Globoturborotalita nepenthes (PL1) in the Atlantic and Pacific oceans (Wade et al., 2011).The same zone was recorded in the Early Pliocene by Rizzini (1978) and Abdou et al. (1984) in the Nile Delta.Conversely, Ismail et al. (2010) determined the Early Pliocene with the Globorotalia margaritae zone.
Assemblage.A quite rich and distinctive planktonic and benthic foraminiferal assemblage has been recorded.See File 1 in the Supplement.
Remarks.According to Blow (1969), this zone coincides with the lower part of his zone N19.It is comparable to the Globorotalia margaritae zone, which is found in the Mediterranean area (Borsetti et al., 1979;Iaccarino and Salvatorini, 1982).According to Lirer et al. (2019) and Lourens et al. (2004), it corresponds to zone MPl2 of Zanclean age (Early Pliocene).In the Nile Delta Basin, the recognized zone belongs to the Globorotalia margaritae zone (Rizzini et al., 1976;Abdel-Kireem et al., 1984), whereas it corresponds to the lowermost part of the Globorotalia puncticulata zone of Abdou et al. (1984).
4.1.6Globorotalia puncticulata/Globorotalia margaritae concurrent range zone (MPl3) Definition.It is defined as the interval characterized by the concomitant presence of Globorotalia puncticulata (after its LO) and Globorotalia margaritae (before its HO).
Occurrence.This zone has been recorded only in the Kafr El-Sheikh Formation of well EW-5 and spans the interval from 2555 to 2460 m (95 m thick; Fig. 7d).Due to the poor preservation state of the foraminifera specimens, this zone is not recorded in the other wells.
Assemblage.The benthic foraminifers recovered from this interval are highly diversified and abundant.See File 1 in the Supplement.
Remarks.It is suggested that the proposed biozone spans the upper zone N19, in accordance with Blow (1969).In the Mediterranean realm, the recognized zone is consistent with zone MPl3 of Lirer et al. (2019) and others.It also matches with the lower part of the Globorotalia puncticulata zone in the eastern Mediterranean region (Borsetti et al., 1979).The identified zone has been correlated in part with the Globorotalia puncticulata zone of Rizzini (1978), Abdel-Kireem et al. (1984), andEl-Heiny et al. (1990) of the Nile Delta Basin and of the Red Sea (Selima, 1998).

Calcareous nannofossil biostratigraphy
The five studied Middle Miocene-Lower Pliocene successions show a diversified calcareous nannofossil assemblage (Figs. 8 and 9).The identification of 54 species, along with their stratigraphic distribution, allowed the definition of 9 calcareous nannofossil zones.Such zones are compared to those published for the Mediterranean realm (Fig. 6).The distribution of various calcareous nannofossil species, their equivalent zones/subzones, and proposed ages in the five wells are illustrated in File 2 in the Supplement.The recorded calcareous nannofossil zones are from bottom to top.

Calcidiscus premacintyrei interval zone (CNM8)
Definition.This zone was defined as the interval from the HO of Sphenolithus heteromorphus to the HO of Calcidiscus premacintyrei.
Occurrence.This zone has been recorded only in well EWE-2, spanning depths from 3050 to 2990 m, where the Sidi Salem Formation has a thickness of ∼ 60 m (Fig. 7a).
Assemblage.See File 2 in the Supplement.NN6 of Martini (1971) and zone CN5a of Okada and Bukry (1980), respectively.Also, the zone corresponds to MNN6a of Fornaciari et al. (1996) and Raffi et al. (2003) for the Mediterranean region.
Age. Middle Miocene (Serravallian).Occurrence.This zone was recorded only in the Sidi Salem Formation of well EWE-2, encompassing depths from 2990 to 2950 m, with a total thickness of ∼ 40 m (Fig. 7a).
Assemblage.See File 2 in the Supplement.
Remarks.Zone CNM9 (Discoaster variabilis) of Backman et al. (2012) corresponds to the upper parts of zones NN6 and CN5a of Martini (1971) and Okada and Bukry (1980), respectively.In the Mediterranean region, this zone is the equivalent of MNN6b of Fornaciari et al. (1996) and Raffi et al. (2003).
Assemblage.See File 2 in the Supplement.
Remarks.Zone CNM14 (Reticulofenestra pseudoumbilicus) of Backman et al. (2012) is correlated with the lower parts of zones NN10 and CN8 (CN8a) of the standard schemes of Martini (1971) and Okada and Bukry (1980) and with MNN10a of Raffi et al. (2003) in the Mediterranean.This zone corresponds to the lowermost part of zone CN8.

Discoaster bellus interval zone (CNM15)
Definition.This zone was defined as the interval from the LO of R. pseudoumbilicus to the LO of Discoaster berggrenii.
Occurrence.This interval zone has only been recorded in well EW-5 and encompasses a depth interval from 3300 to 3200 m (100 m thick) within the Sidi Salem Formation (Fig. 7d).
Assemblage.See File 2 in the Supplement.
4.1.12Amaurolithus primus interval zone (CNM17) Definition.This zone was drawn from the LO of Amaurolithus primus to the LO of Nicklithus amplificus.
Assemblage.See File 2 in the Supplement.Remarks.Zone CNM17 (Amaurolithus primus) of Backman et al. ( 2012) is comparable to the middle part of zone NN11 of Martini (1971) and to the lower part of zone CN9b of Okada and Bukry (1980).In the Mediterranean region, this zone corresponds to zone MNN11b of Raffi et al. (2003).
4.1.13Nicklithus amplificus total range zone (CNM18) Definition.This zone corresponds to the total range of N. amplificus.
Assemblage.See File 2 in the Supplement.
Remarks.Zone CNM18 (Amaurolithus primus) of Backman et al. ( 2012) corresponds to the upper parts of zone NN11 of Martini (1971) and subzone CN9b (Amaurolithus primus) of Okada and Bukry (1980).It is the equivalent of MNN11c defined by Raffi et al. (2003) in the Mediterranean region.
Assemblage.See File 2 in the Supplement.

Sphenolithus neoabies partial range zone (CNPL2)
Definition.It is a partial range between the HO of C. armatus and the lowest common occurrence (LCO) of Discoaster asymmetricus.
Assemblage.See File 2 in the Supplement.
4.1.16Discoaster asymmetricus/Reticulofenestra pseudoumbilicus concurrent range zone (CNPL3) Definition.The zone is defined based on the consecutive occurrence of the lowest common occurrence of D. asymmetricus and the HO of R. pseudoumbilicus.Occurrence.Zone CNPL3 was recorded in wells EW-8 (from 2630 to 2585 m), EW-12 (from 2570 to 2550 m), and EW-5 (from 2470 to 2460 m), corresponding to the Kafr El-Sheikh Formation.The corresponding thickness is 45, 20, and 10 m, respectively (Fig. 7b-c and e).
Assemblage.See File 2 in the Supplement.

Palaeoenvironmental reconstruction
The integration of the foraminifera and calcareous nannofossil assemblages allowed reconstruction of the palaeoenvironmental conditions in the study area based on proxies such as benthic foraminiferal depth indicators, P /B ratio, species richness and abundance, and calcareous nannofossil taxa characteristics of peculiar environmental conditions.In the Serravallian, the occurrence of the calcareous nannofossil C. floridanus is reported by some authors to indicate open-marine conditions (Wei and Wise, 1990;Aubry, 1992;Monechi et al., 2000).The benthic foraminiferal association which occurred in this time interval (see File 1 in the Supplement) and is illustrated in Fig. 10 suggests a middleto outer-neritic environment (Murray, 1991;Holbourn et al., 2013).This result suggests that the depositional environment for the Sidi Salem Formation during the Serravallian was an open shelf (middle-to outer-neritic).
In the Messinian, the calcareous nannofossil assemblages were dominated by Sphenolithus spp.and by small reticulofenestrids, along with a reduced Discoaster percentage and the sporadic occurrence of A. beccarii.The occurrence of rare ostracod shells might indicate marginal estuarine conditions, as suggested by Wade and Bown (2006) and Chakraborty et al. (2021).
In the Early Pliocene (Zanclean), a peak in abundance of calcareous nannofossils and foraminifera (both benthic and planktonic) is recorded in the studied cores.A marked increase in the Discoaster percentage might indicate an openmarine environment (Chakraborty et al., 2021).This interpretation is further supported by the occurrence of planktonic foraminifera such as Sphaeroidinellopsis spp. that are also considered to be open-marine taxa (Wei et al., 1992;Iaccarino et al., 2007).Such conditions are consistent with a global sea-level rise already suggested by Snedden and Liu (2010).Warm-water indicators, such as C. armatus, C. cristatus, O. universa, and Globigerinoides spp.(Sinaci, 2013;Verducci et al., 2009), are also recorded in the earliest Zanclean.Consequently, calcareous nannofossil and foraminifer records might confirm that the Kafr El-Sheikh Formation was deposited under open-marine conditions.

Sequence stratigraphy
Our high-resolution biostratigraphic investigation allows determination of the depositional sequences and system tracts of the Neogene sediments of the Nile Delta system.Facies types are also characterized using a detailed lithologic description incorporated with wireline logs analysis.The use of CycloLog software to generate PEFA and INPEFA patterns has enabled sequence identification (Nio et al., 2005).Many criteria, such as biozonal breaks and subaerial exposures, have been used to inspect sequence boundaries, whereas palaeobathymetry was used to clarify system tracts based on benthic biofacies, planktonic-benthic ratio, abundance, and diversity.Several studies have shown the usefulness of these parameters (e.g.Armentrout, 1996;Leckie and Olson, 2003;Miller et al., 2008).Biostratigraphy allowed the dating of flooding and maximum flooding events.The ages of these events were calibrated with the Neogene timescale of Backman et al. (2012) and Lourens et al. (2004) by using the marker species of both calcareous nannofossils and planktonic foraminifera.The flooding events and depositional sequences have been recognized and correlated with the global cycle chart adopted by Snedden and Liu (2010).Also, the sequence boundaries of this study were correlated with those published for the Nile Delta and the Mediterranean by Dolson et al. (2005), Kellner et al. (2009), andFarouk et al. (2014); other regional sequence boundaries are shown by Haq et al. (1987), Hardenbol et al. (1998), andWornardt (1999) (Fig. 11).The investigated succession has been subdivided into seven major third-order depositional sequences, bounded by six sequence boundaries related to euhttps://doi.org/10.5194/jm-42-147-2023 Figure 10.Upper depth limit of benthic foraminifera with respect to P /B ratio after Van Morkhoven et al. (1986), Van der Zwann et al. (1990), Murray (1991), andHolbourn et al. (2013).static sea-level changes and/or local tectonics.This is based on the results of the integration between wireline logging, seismic data, and high-resolution biostratigraphic data.

Sequence boundaries of El-Wastani field
Figure 7 shows six sequence boundaries that were correlated with the global sequences proposed by Snedden and Liu (2010) (see File 3 in the Supplement).The unconformity surfaces could be regarded as sequence boundaries (SBs) (Kassab et al., 2015), which define seven depositional sequences in our study area.These sequence boundaries are unconformity surfaces related to subaerial exposure named type 1, such as SB1, SB2, SB3, SB4, and SB5, or submarine erosional unconformity named type 2, such as SB6, as shown in Fig. 7 and File 3 in the Supplement.
The upper SB of the Serravallian sequence (Serseq) is an erosive unconformity surface of type 1 (SB1), which could be dated to 12.18 Ma through correlation with the paper of Kellner et al. (2009).The lower SB could be correlated with the previously defined SB at 12.70 Ma as published by Hardenbol et al. (1998) and Kellner et al. (2009) (Fig. 11).
The upper SB (SB2) of the first Tortonian sequence (Torseq1) is represented by an erosional surface and angular unconformity due to uplift, forming a biostratigraphic break in well EWW1-St2, and could be correlated with the SB at 8.8 Ma identified by Wornardt (1999); the top of this sequence is a maximum flooding surface (MFS).The lower SB of the Torseq1 was not reached in our study, but it might be attributed to the SB described by Hardenbol et al. (1998) and Kellner et al. (2009) at 9.26 Ma (Fig. 11).
The upper SB of Torseq2 (SB3) is characterized by an angular unconformity between the tilted blocks of the Sidi Salem Formation and the overlying Qawasim Formation, as shown in the seismic section of well EW-5 (Fig. 11).The lower SB of Torseq2 was not reached; however, it may be correlated with the sequence boundary defined by Wornardt (1999) at 8.60 Ma (Fig. 11).
The upper SB (SB4) of the first Messinian sequence (Meseq1) is characterized by local unconformity inferred from the incised valley truncation by fluvial input in the lower Abu Madi Formation, as shown by the NNE-SSWtrending seismic line that crosses the valley direction (NNW-SSE) (Fig. 12a-b).The Meseq1 upper SB (local unconformity) is observed in all the studied wells and could be correlated with the SB identified by Hardenbol et al. (1998) and Kellner et al. (2009) at 6.98 Ma.The lower SB was not reached in wells EW-8 and EW-12, but it could be correlated with the SB at 7.24 Ma defined by Farouk et al. (2014) in the Nile Delta (Fig. 11).
The upper SB of sub-cycle Meseq2a represents another phase of valley incision, which truncated older beds, as shown by seismic sections (Fig. 12a-b).The upper SB (SB5) of Meseq2 corresponds to a type 1 erosive surface due to subaerial erosion (zones CNM19 and CNM20), dated to the latest Messinian.
The upper SB (SB6) is a submarine erosive surface (type 2), as indicated by lithofacies and gamma-ray signals.This SB is correlated with the ones recognized at 4.36 Ma by Hardenbol et al. (1998) and Kellner et al. (2009) and at 4.37 Ma by Snedden and Liu (2010) (Fig. 11).

Depositional sequences of El-Wastani field
The five studied wells contain seven depositional sequences encompassing the Serravallian to Zanclean periods (Fig. 7).These are Serseq1, Torseq1, Torseq2, Meseq1 (a + b), Meseq2 (a +b), Zaseq1, and Zaseq2.The description of each depositional sequence and sequence boundary in the five wells is as follows.

I. Serravallian sequence (Serseq)
The lowermost Miocene depositional sequence (Serseq) occurs in the Sidi Salem Formation and corresponds to foraminiferal subzone MMi7b and calcareous nannofossil zones CNM8 and CNM9.Serseq is recorded in well EWE-2 (Fig. 7a).Serseq is composed of shale beds that reflect an open-shelf (middle-to outer-neritic) environment as inferred from benthic foraminifer and calcareous nannofossil contents (Fig. 7).It corresponds to retrograding parasequences defining the TST (transgressive system tract).A gap occurred corresponding to the uppermost Serravallian-lower Tortonian sequences, as also shown by the absence in the studied well of some foraminiferal zones (MMi8, MMi9, MMi10), as well as of some calcareous nannofossil zones (CNM8 to CNM13).
SB1 is marked by a long hiatus occurring between the Serravallian Sidi Salem and the overlying Messinian Qawasim formations.The regional tectonic uplift of the Serravallian tilted blocks in the eastern part of Africa was the main driver for this hiatus, easily seen in the seismic profile (Harms and Wary, 1990).A shift occurs in the depositional environment, from middle-to outer-neritic deposits in the Serravallian to fluvial-marine and lagoon facies of the Messinian age, as inferred from the rare and poorly diversified benthic foraminifera (see File 1 in the Supplement).Subzone MMi7b and nannofossil zones CNM8 and CNM9 are recorded below SB1.The latter surface is overlain directly by the Messinian NDZ and by zones CNM17 and CNM18.The hiatus and the observed change in environmental conditions were likely driven by local tectonics in the Nile Delta area due to the rifting in the Gulf of Suez (Farouk et al., 2014).
The TST is characterized by the deposition of middle-to outer-neritic shales, rich in foraminifera and calcareous nannofossils (subzone MMi7b, zones CNM17 and CNM18), because of a relative sea-level rise.The foraminiferal assemblages show a P /(P + B) ratio of 39 % to 67 %, high species richness (25-36), and a moderate to good preservation state (see File 1 in the Supplement).

II. First Tortonian sequence (Torseq1)
This sequence is recorded in well EWW1-St2 (Fig. 7e).It corresponds to the upper Sidi Salem Formation (zone MMi11 and the lower part of zone CNM14).The Torseq1 is characterized by a regional angular unconformity between the tilted blocks of the Sidi Salem Formation and the Qawasim Formation, as shown by the NW-SE seismic section (Fig. 12a-b) and previously inferred for the Nile Delta by El-Heiny and Morsi (1992).Torseq1 corresponds to a TST made of retrograding parasequences (Fig. 7e).
The TST is composed of middle-to outer-neritic shales yielding high planktonic foraminiferal and calcareous nannofossil contents.This interval shows a P /(P + B) ratio of ∼ 76 % and moderate to high foraminiferal diversity (benthic and planktonic, ∼ 17 species), with moderate to good preservation.The benthic foraminiferal assemblage comprises deep-water taxa (see File 1 in the Supplement).High abundance and diversity of calcareous nannofossils are also observed in the TST interval.

III. Second Tortonian sequence (Torseq2)
It is recorded in well EW-5 (Fig. 7d) and encompasses the uppermost Sidi Salem Formation (zone MMi11, zone CNM15).The top of Torseq2 occurs near the Tortonian-Messinian boundary and is accompanied by significant change in foraminiferal assemblages, which shift from moderate to high abundances to a barren interval.Some plankhttps://doi.org/10.5194/jm-42-147-2023tonic foraminiferal and calcareous nannofossil zones are missing, namely MMi11 to MMi13 and CNM17 to CNM18.Torseq2 contains only the highstand system tract (HST) prograding parasequences.
The HST is composed of inner-to middle-neritic shale yielding moderately abundant and diversified planktonic foraminifera and calcareous nannofossils.It is characterized by a P /(P + B) ratio of ∼ 45 %, with a moderate foraminiferal diversity (benthic and planktonic, 8-17 species) and a moderate preservation state (see File 1 in the Supplement).

IV. First Messinian sequence (Meseq1)
It corresponds to the lowermost Messinian Qawasim Formation in the five studied wells (Fig. 7) and to zones NDZ and CNM17.The Meseq1 sequence comprises the lowstand system tract (LST)/TST and HST as retrograding/prograding and aggrading parasequences.It is well represented in well EW-12 (Fig. 7c).
The LST is a prograding sedimentary body, as shown by the deposits of the lower Messinian Qawasim succession.Lithofacies are mainly composed of sandy deposits with grains of variable size, corresponding to a submarine canyon formed in the context of sea-level drop.These LST deposits are well recognized in well EW-8 (Fig. 7c).An interval barren in planktonic and benthic foraminifera and with low to moderate calcareous nannofossil abundance was observed in some wells (EW-8, EW-12, EWW1-St2) (see Files 1 and 2 in the Supplement).
The TST is dominated by sands with shale interbeds.The TST deposits occur in the five wells (EWE-2, EWW1-St2, EW-5, EW-8, EW-12).Calcareous nannofossils show a trend of increasing abundance upward, whereas only shallow benthic foraminifera are recorded (see Files 1 and 2 in the Supplement).
The MFS could be placed in the first four wells and represents the maximum abundance and diversity of the calcareous nannofossils.The overlying HST deposits exhibit low to moderate nannofossil abundances.

V. Second Messinian sequence (Meseq2)
It is recorded in wells EWW1-St2, EWE-2, EW-5, EW-8, and EW-12 and corresponds to the calcareous nannofossil zones CNM17 and CNM18 and to the uppermost part of the foraminiferal NDZ of the Abu Madi Formation.The Meseq2 depositional sequence is composed of two subcycles, Meseq2a and Meseq2b (Fig. 7).The lower and upper SBs of the Meseq2a and Meseq2b sub-cycles present multiple phases of valley incision and truncation of the beds (Fig. 12a-b).Due to such features, reworked calcareous nannofossils commonly occur in the sediments at the base of each sub-cycle.
According to Catuneanu et al. (2011), the LST follows a regressive sequence, which deposited at the beginning of sea-level drop and the subsequent relative sea-level rise.The fluvial-deltaic system throughout the delta region is represented by a prograding, lowstand clastic unit (Abu Madi Fm.), which is interpreted as having deposited during sealevel fall (Abdel Hady et al., 2013).High percentages of reworked microfossils (from Cretaceous, Palaeogene, and Neogene) are recorded in the Messinian Abu Madi Fm.Catuneanu et al. (2011) showed that the infilling of incised valleys began in the final stage of the LST.In the five studied wells, the lowstand facies in the two sub-cycles of Meseq2 represent the thickest part of the entire sequence, being composed of a sandstone unit (Fig. 7).Because of its porosity, the LST of the Abu Madi Fm. represents the most important Miocene gas reservoir in the Nile Delta system.
The TST is dominated by shales and sand interbeds, as shown by the gamma-ray signal and INPEFA log pattern (Fig. 7).Also, a slight increase in the calcareous nannofossil abundance occurs during the retrograding parasequences.Conversely, planktonic foraminifera are absent, and only A. beccarii characterizes the benthic forms.The TST facies occur in the five investigated wells, in both the Meseq1a and Meseq2b (Fig. 7).
The Meseq2b MFS has been tentatively placed at the top of the transgressive sequence based on the record of high calcareous nannofossil abundances and diversity.The overlying regressive sequence (HST) contains slightly fewer calcareous nannofossils and some shallow-water benthic foraminifers such as Elphidium sp. and A. beccarii.

VI. First Zanclean sequence (Zaseq1)
Zaseq1 is the oldest Pliocene depositional sequence that comprises the Kafr El-Sheikh Formation and spans zones MPl1 and MPl2 and the lower part of zone MPl3, correlated with calcareous nannofossil zones CNPL1 and CNPL2.
The TST is placed in the lower part of the Kafr El-Sheikh Formation (zones MPl1, MPl2, and CNPL1).The TST retrograding facies is recorded from four wells (EWW1-St2, EWE-2, EW-8, and EW-12), as shown by gamma-ray signals and INPEFA log patterns (Fig. 7).The highest abundance https://doi.org/10.5194/jm-42-147-2023J. Micropalaeontology, 42, 147-169, 2023 and species diversity of both planktonic foraminifera and calcareous nannofossils are observed in Zaseq1.The high occurrence of deep-water benthic foraminifer indicators (see File 1 in the Supplement) suggests an outer-neritic-upper-bathyal environment, also compatible with a P /(P + B) ratio up to 95 %.The MFS is identified within the Kafr El-Sheikh Formation in four wells (Fig. 7).This surface corresponds to the highest abundance of planktonic and benthic foraminifera and calcareous nannofossils.Foraminifera specimens are coated by pyrite, whose formation is compatible with dysoxic conditions within the sediment.The MFS is correlated with the analogous event described by Haq et al. (1987) and Dolson et al. (2005) at 5.0 Ma and by Wornardt (1999) at 5.05 Ma (Fig. 11).The MFS of Zaseq1 is thus one of the most relevant surfaces near the top of zone MPl1 and the calcareous nannofossil zone CNPL1 at 5.0 Ma.
The HST facies is recorded from all five wells (Fig. 7).The HST is composed of neritic shale, with few sands and silt interbeds.The interpretation of electric logs and CycloLog data suggests the occurrence of a continuous sea-level fall and a change from deep, open-marine facies to middle-and outer-neritic environments.The P /(P + B) ratio shows values of 53 %-68 %, with a moderate to high species diversity of foraminifera and calcareous nannofossils, as well as a good preservation of benthic forms.

VII. Second Zanclean sequence (Zaseq2)
This depositional sequence occurs in wells EW-5 and EW-12 (Fig. 7) in the lowermost part of the Kafr El-Sheikh Formation.In well EW-5, the depositional sequence is dated to the upper part of MPl3 and the calcareous nannofossil zone CNPL3.Well EW-12 corresponds to the upper part of zones MPl2 and CNPL3.
The LST/TST interval in the Kafr El-Sheikh Formation is composed of light to brownish-grey shales with traces of organic matter.The P /(P + B) ratio is 45 %-67 %, with an upwards increase in relative abundance and species diversity of both foraminifera and calcareous nannofossils.Benthic foraminifera are dominant in this interval (see File 1 in the Supplement) and suggest deposition in an inner-to middleneritic environment.

Correlation of the depositional sequences with previous works
The Middle-Late Miocene (Serravallian/Tortonian) boundary is characterized by a small-amplitude sea-level fall in the Tethyan realm (Hardenbol et al., 1998;Snedden and Liu, 2010).In the Serravallian, the eastern part of Egypt was affected by tectonic activity (Dolson et al., 2001).This caused uplift and tilting to the north of the older rock units.Accordingly, the oldest Serravallian sequence (Serseq) is only recorded in well EWE-2, corresponding to the south-eastern part of the El-Wastani field (Fig. 7).This tectonic activity is responsible for a long hiatus between the Serravallian sequence, which is lacking in the northern wells, and the overlying Messinian sequence.The Serseq sequence is interpreted here as being made of TST deposits.This interpretation is in contradiction to Farouk et al. (2014), who interpreted HST facies for the same interval.On the north-western side of the study area, wells EWW1-St2 and EW-5 are composed of Tortonian sediments.Both wells were affected by tectonic activity that occurred in the Serravallian, hence creating the upper SBs of Torseq1 and Torseq2.The two Tortonian sequences are correlated with the ones recorded by Farouk et al. (2014), even though these authors interpreted the lower part of the sequence as an LST.
During the Messinian time, the drawdown of the Mediterranean and, in turn, of the Nile Delta incised channels cut the upper parts of the Serravallian-Tortonian Sidi Salem Formation; such channels were infilled by the Qawasim unit (Meseq1).The sequence thickness increases in the central part of the field (well EW-8), which may represent the deepest part of the channel, in agreement with Farouk et al. (2014).The Abu Madi Formation was formed under the influence of multiphases of fluvial channels cut, which form the upper SB of Meseq2.The thickness of Meseq2 is important along the margin of the field; thus the main channel pass was likely located between wells EWW1-St2 and EWE-2, whereas minor channels occurred in the central field.Meseq2 is also recorded by Farouk et al. (2014) in the same time interval.
The lowermost Zanclean (Zaseq1, Kafr El-Sheikh Formation) likely corresponds to the first arrival of marine water from the Atlantic Ocean through the Gibraltar seaway (Roveri et al., 2014) (Fig. 7).This sequence corresponds to TST/HST in four of the investigated wells, excluding EW-5, which only bears a HST due to the lack of zone MPl1 and zone CNPL1.This lack may be due to the rarity of the zonal markers in well EW-5.The Zaseq1 thickness is important in wells EW-5 and EW-12, probably due to fault initiation in the Pliocene (Kamel et al., 1998) or because the deepest part of the basin was located around these wells.Well Zaseq1 matches with the sequence reported by Farouk et al. (2014).

Sequence stratigraphy and hydrocarbon potential
In the present work, the analysis of the lower and upper Messinian rock units (Qawasim and Abu Madi formations) in the framework of a sequence stratigraphy study allows an effective correlation for exploration with other onshore concessions of the Nile Delta.The LSTs correspond to the most important gas reservoirs due to the occurrence of sandstone deposits (Catuneanu et al., 2011).These thick sandstone units have been observed in the LST of the Meseq1 and Meseq2 of the Messinian Qawasim and Abu Madi formations of the present study.Hydrocarbon occurrence in a lowstand sequence is influenced by charge, seals, and source rocks.According to Catuneanu et al. (2011), fluvial, coastal, and shallow-water lowstand reservoirs are sealed by overlying evaporite layers and transgressive shales.In the El-Wastani field, three prograding lowstand wedges were identified, two belonging to the Abu Madi Formation and the third one belonging to the Qawasim Formation, based on our high-resolution biostratigraphy and geophysical log patterns.The LST deposits are sealed by marine shales.The Messinian channel infilling (Qawasim and Abu Madi Fm.) in a bounded-fault context, locally bearing unconformities due to occurrence of fluvial channels, acted as stratigraphic traps (Fig. 12a, b).The MFS deposits may have acted as source rocks for hydrocarbon generation (Dolson et al., 2002).

Conclusions
A detailed biostratigraphic study based upon planktonic and benthic foraminifera and calcareous nannofossils was conducted in five wells located within the El-Wastani field, onshore Nile Delta, Egypt, and allowed identification of Middle Miocene-Early Pliocene sequences.Six biozones based on planktonic foraminifera markers and nine biozones based upon calcareous nannofossils were established.These allowed correlation with the most relevant Miocene-Pliocene biozonal schemes from the Mediterranean region and middle latitudes.Benthic foraminifera and calcareous nannofossil assemblages were used to reconstruct depositional environments, which range from lagoon to upper-bathyal environments.The integration of sequence stratigraphy, wireline logs (gamma rays), and biostratigraphy allowed the identification of seven depositional sequences.These sequences highlighted the reservoir-bearing intervals of the Qawasim and Abu Madi formations that proved to be incised valleys over older rock units, as shown by the NE-SW seismic lines.Our study will enable the future exploration in the Nile Delta provinces.Paragloborotalia mayeri Cushman and Ellisor, 1939, Figs. 3.3-3.4The specimens illustrated possess five inflated chambers.The first chamber is small and clearly apparent in ventral view.The spiral suture is slightly curved and depressed.These specimens display an increasing chamber size, which is characteristic of the Middle Miocene P. mayeri.These criteria are compatible with the holotype description of P. mayeri.

Figure 6 .
Figure 6.Biostratigraphic correlation chart of the proposed Middle Miocene-Early Pliocene biozones with other biozones in the Mediterranean realm.

4. 1 . 9
Discoaster variabilis partial zone (CNM9)Definition.According toWade et al. (2011) andBackman et al. (2012), this partial range zone corresponds to a body of rock located between the HO of a given taxon and the LO of another taxon.Consequently, we defined the Discoaster variabilis zone as a partial range zone comprising the HO of Calcidiscus premacintyrei and the LO of Cyclicargolithus floridanus.
It comprises the HO of Discoaster hamatus and the LO of Discoaster loeblichii.

Figure 11 .
Figure 11.Sequence stratigraphic boundaries of the onshore Nile Delta in comparison with global and other related geographical regions.Chronostratigraphic scale after Ogg et al. (2016).

Figure 12 .
Figure 12.NE-SW seismic-to-well-tie profile illustrating the rock units, stratum terminations, and stratigraphic boundary of each rock unit of wells EWW1-St2 and EW-5 in the El-Wastani gas field.