Eocene to Oligocene nannofossils stratigraphy and environmental conditions in Izeh Province,Zagros Basin,East Tethys

Saeedeh Senemari and Farah Jalili

Abstract Data obtained from the calcareous nannofossils, distributed in the upper part of the Pabdeh Formation(Priabonian-Rupelian) and the lower part of the Asmari Formation (Chattian) in the Bid-Zard section, were used to investigate the Eocene to Oligocene palaeoenvironmental conditions in the southwest of Izeh, southwestern Iran(eastern Tethys). The upper part of the Pabdeh Formation was composed of shale, thin-bedded pelagic limestone and dolostone, which is disconformably overlain by the Asmari Formation. For the first time, 29 species of calcareous nannofossils belonging to 13 genera were identified in the studied section. The calcareous nannofossils in the upper part of the Pabdeh Formation indicate the Isthmolithus recurvus Zone/Sphenolithus pseudoradians Zone(combined zone), Ericsonia subdisticha Zone, Helicosphaera reticulata Zone and Sphenolithus praedistentus Zone,from the Priabonian to the Rupelian. The Sphenolithus ciperoensis Zone of the Chattian was identified in the lower part of the Asmari Formation. Calcareous nannofossil stratigraphy across the upper Eocene-Oligocene interval also reveals a disconformity at the Rupelian/Chattian transition due to a bio-event. Shallowing of the basin and environmental changes in this part of the Tethyan domain could have led to the lithostratigraphic and biostratigraphic changes. In fact, during the late Eocene to late Oligocene, marine phytoplankton was sensitive to climate changes such as decreasing temperature, as well as possibly to a nutrient increase and changes in basin depth.

Keywords: Biostratigraphy,Calcareous nannofossils,Ecological indicators, Tethys domain, Zagros, Paleogene

1 Introduction

Marine sediments of the Pabdeh and Asmari formations at the Eocene-Oligocene transition are one of the most important source rocks and hydrocarbon reservoirs at the margin of the Zagros Basin in southwestern Iran(Kamali et al. 2006; Bordenave and Hegre 2010). The Zagros Basin is located at the collisional zone of the Arabian Plate and Iranian Block, and includes a very thick marine sedimentary succession (Al-Husseini 2000;Lacombe and Mouthereau 2006) (Fig. 1). In fact, the Zagros Mountains of SW Iran are a part of the Alpine-Himalayan orogenic system, which is one of the most important sedimentary areas in the world with a NWSE trend (Fig. 1b). The basin is composed of a gently folded rock succession with several stratigraphic provinces such as Lurestan, Izeh, Fars, and Khuzestan (including the Dezful Embayment) (Motiei 1994;Darvishzadeh 2003; Alavi 2004). The Pabdeh Formation was mostly deposited during the Paleocene-Oligocene and is present in all these areas in both outcrops and subsurface. After the Late Cretaceous compressive phase(Laramid phase), the sea level rose and the Zagros Basin deepened progressively during the Paleogene, which led to the deposition of the Pabdeh Formation in the deeper parts (Aghanabati 2004). Then, as the global sea level fell, the basin became shallower toward the top of the Pabdeh Formation.This was followed at the top of the Pabdeh Formation by a transitional shale facies and by deposition of anhydrite of the Asmari Formation (Amin-Rasouli et al. 2012). Thus, the Pabdeh Formation is overlain by the Asmari Formation.However,disagreement still exists about whether the boundary between the Pabdeh and Asmari formations is continuous or discontinuous, as well as about the exact age of the basal anhydrite and its overlying transitional sediments in the Asmari Formation (e.g., James and Wynd 1965; Motiei 1995; Seyrafian 2000; Vaziri-Moghaddam et al.2006;Amirshahkarami et al.2007).

Fig. 1 a Palaeogeographic map during the late Eocene (37-34 Ma) modified from Popov et al. (2004); b Approximate position of the Alpine-Himalayan orogenic system, which is shown with a gray line from Europe to China, and of the Zagros Basin, which is shown with a red line in southwestern Iran (adapted and modified from World Map-A non-Frame.png on https://en.wikipedia.org/wiki/Alpide_belt)

Fig. 2 a The Zagros study site locality on the Iran map; b Tectonic outline map and sub-divisions of the Zagros Basin in southwestern Iran; c Approximate geographical location of the Bid-Zard section marked with an asterisk

In recent decades, the Pabdeh Formation has been considered by most geologists as one of the important source rocks in the Zagros Basin (e.g., Amirkhani et al.2015). Biostratigraphic studies on the Zagros Basin, especially on the Pabdeh Formation, were mostly based on foraminifera (e.g., Behbahani et al. 2010). These studies recognized a conformable boundary between the Pabdeh and Asmari formations. In addition, the Mesozoic and Cenozoic deposits in the Zagros Basin have been extensively studied using calcareous nannofossils (e.g.,Gholami Fard 2007; Hadavi et al. 2007; Senemari 2007;Sadeghi and Hadavand Khani 2010; Khavari Khorasani et al. 2014; Senemari 2014; Senemari and Bakhshandeh 2014; Ahifar et al. 2015).This paper presents the first results of the nannofossil stratigraphy of the Pabdeh and Asmari formations.

In this study, calcareous nannofossils from the Bid-Zard section in the Dezful Embayment Zone (Fig. 2), located about 50 km southwest of Izeh, are sampled and analyzed. The study region covers the area between longitudes 49°15′E-49°45′E and latitudes 31°30′N-31°40′N. The Bid-Zard section was sampled every 1-2 m from the upper part of the Pabdeh Formation and the lower part of the Asmari Formation (Fig. 3). The main objectives are to document and identify calcareous nannofossil assemblages, establish a biostratigraphic framework,and evaluate depositional/environmental conditions across the Pabdeh/Asmari boundary, in order to investigate source and reservoir rocks in southwestern Izeh.This evaluation, based on the calcareous nannofossils,highlights the significance of biological events anddemonstrates that the boundary between the Pabdeh Formation and the Asmari Formation was discontinuous and that the water depth of the basin decreased at the end of the Rupelian.

Fig. 3 Lithology and distribution of calcareous nannofossils across the upper Eocene to Oligocene from the Bid-Zard section southwest of Izeh,Zagros Basin. Red dots mark the first and last occurrence of the species; black dots are the spread of species. Note that dolomite and anhydrite facies without calcareous nannofossils are shown larger for greater clarity and do not correspond to the scale

2 Material and methods

The studied samples belong to the Bid-Zard section in southwestern Izeh, a 129-m-thick sedimentary succession containing the upper part of the Pabdeh Formation and the lower part of the Asmari Formation.Seventy-eight fresh rock samples from shales and thin-bedded limestones were collected for calcareous nannofossils based on lithological changes.The section was sampled at about 1-1.5 m intervals;combined results of every three samples are shown in Figs.3 and 4. The upper part of the Pabdeh Formation in the section is composed of shale, thin-bedded pelagic limestone,and dolostone. The shale facies is laminated, dark-colored,rich in organic matter,and contains thin-bedded limestones.The abundance of organic matter indicates that the facies formed in the absence of oxygen. The thin-bedded pelagic limestones contain phosphate grains, glauconite, and planktonic foraminifers.Due to these features,the shale facies was inferred to have formed above the carbonate compensation depth(CCD).At the top of the Pabdeh Formation,deep pelagic dolomitic rocks(samples 73-78)were deposited including a 4.2-m-thick microbial facies and stromatolites. The microbial facies consist of dark and light as well as wavy laminae.Stromatolites also consist of dark and light layers,flat to wavy, and contain a fenestral fabric. The boundary between the Pabdeh and Asmari formations is marked by a red paleosol horizon.Subsequently,the lowest part of the Asmari Formation, mainly composed of a basal blocky anhydrite,was deposited. An about-2.45-m-thick layer of the basal anhydrite (samples 79-84) probably formed as a result of a temporary disconnection with the open sea (regression).Eventually, the deep deposits of the transitional zone of the Pabdeh and Asmari formations were formed with shale facies. In fact, during the depositional process in transitional zone following transgression, shale facies was deposited(samples 85-90) (Figs. 3 and 4). Twelve samples were also taken from both dolomite and anhydrite facies in the transition zone of the two formations.

Fig. 4 Diversity and abundance of calcareous nannofossils across the upper Eocene to Oligocene from the Bid-Zard section southwest of Izeh,Zagros Basin. Note that dolomite and anhydrite facies without calcareous nannofossils are shown larger for greater clarity and do not correspond to the scale

Selected samples were prepared following the standard smear slide method (Bown and Young 1998). All the slides were studied under the crossed-and plane-polarized light of a BH2 Olympus microscope with 1000-fold magnification and photographed by a Canon digital camera.Palaeoenvironmental reconstruction of the section was performed based on semi-quantitative analyses of selected calcareous nannofossils. For each slide, five to seven transverse rows were counted to identify the calcareous nannofossils in order to observe changes in the assemblage of each taxon.The abundance of calcareous nannofossils (per fields of view/per FOVs) was assessed as: R=Rare (1 specimen per 11-100 FOVs);F=Few (1 specimen per 2-10 FOVs);C=Common(1-10 specimens per 1 FOVs); A=Abundant (＞10 specimens per 1 FOVs).The most common species are illustrated in Fig.5.For identification of the calcareous nannofossils,the taxonomic concepts of Perch-Nielsen(1985)and the Nannotax website(http://www.mikrotax.org/Nannotax3)were used.For biostratigraphic analysis, zonations of Martini (1971),Okada and Bukry (1980), and Agnini et al. (2014) were applied. For the Paleogene nannoplankton (NP) abbreviations of the Martini(1971)zonation were used,for Paleogene coccoliths (CP) those of the Okada and Bukry (1980) zonation,and for Eocene/Oligocene calcareous nannofossil (CNE/CNO) those of the Agnini et al. (2014) zonation. The bioevents were delineated as the first occurrence (FO) and the last occurrence(LO)of index species.

Calcareous nannofossil taxa identified in this study(Fig. 5) are given based on Perch-Nielsen (1985) in alphabetical order as follows:

- Blackites spinosus (Deflandre and Fert 1954) Hay and Towe 1962

- Braarudosphaera bigelowii (Gran and Braarud 1935)Deflandre 1947

- Coccolithus miopelagicus Bukry 1971

- Coccolithus pelagicus (Wallich 1877)Schiller 1930

- Cribrocentrum isabellae Catanzariti, Rio and Fornaciari in Fornaciari et al. 2010

- Cribrocentrum reticulatum Perch-Nielsen 1971

- Cyclicargolithus floridanus (Roth and Hay, in Hay et al. 1967) Bukry 1971

- Dictyococcites bisectus (Hay, Mohler and Wade 1966) Bukry and Percival 1971

- Discoaster barbadiensis Tan Sin Hok 1927

- Discoaster deflandrei Bramlette and Riedel 1954

- Discoaster saipanensis Bramlette and Riedel(1954)

- Discoaster tanii Bramlette and Riedel 1954

- Ericsonia formosa (Kamptner 1963) Haq 1971=Coccolithus formosus(Kamptner1963)Wise 1973

- Helicosphaera bramlettei (Müller 1970) Jafar and Martini 1975

- Helicosphaera compacta Bramlette and Wilcoxon 1967

- Helicosphaera euphratis Haq 1966

- Pontosphaera enormis (Locker 1967) Perch-Nielsen 1984

- Pontosphaera multipora (Kamptner 1948 ex Deflandre in Deflandre and Fert 1954)Roth et al.1970

- Reticulofenestra bisecta (Hay, Mohler and Wade 1966) Roth et al.1970

- Reticulofenestra minuta Roth et al. 1970

- Reticulofenestra reticulata (Gartner and Smith 1967)Roth and Thierstein 1972

- Reticulofenestra stavensis (Levin and Joerger 1967)Varol 1989

- Reticulofenestra umbilicus (Levin 1965) Martini and Ritzkowski 1968

- Sphenolithus dissimilis Bukry and Percival 1971

- Sphenolithus distentus (Martini 1965) Bramlette and Wilcoxon 1967

- Sphenolithus moriformis (Br?nnimann and Stradner 1960) Bramlette and Wilcoxon 1967

- Sphenolithus praedistentus Bramlette and Wilcoxon 1967

- Sphenolithus pseudoradians Bramlette and Wilcoxon 1967

- Zygrhablithus bijugatus (Deflandre in Deflandre and

Fert 1954)Deflandre 1959

3 Zonal biostratigraphy

The calcareous nannofossil group is a powerful tool for studying stratigraphic correlations on a global and regional scale, and has been discussed in detail by various researchers (e.g., Martini 1971; Okada and Bukry 1980; Perch-Nielsen 1985; Gradstein et al.2012). In the present study, lithologies of the upper part of the Pabdeh Formation and the lower part of the Asmari Formation in the Bid-Zard section were examined. A total of 13 genera and 29 species of calcareous nannofossils were identified, which are well preserved in most slides and document general changes through the section.

The zonal boundaries were determined by the comparison with the standard biozonations of Martini(1971) and Agnini et al. (2014). Accordingly, our data based on calcareous nannofossil events allow the studied section to be assigned to the NP19/NP20 combined zone to the NP25 zone, indicating a Priabonian to Chattian interval. This study also focuses on the environmental conditions of the calcareous nannofossils, which are index fossils at low latitudes during the late Eocene to late Oligocene.

Five biozones in the studied section are introduced in the following. They are equivalent to the NP19/NP20 combined zone to NP23 zone in the upper part of the Pabdeh Formation and to the NP25 zone in the lower part of the Asmari Formation according to the Martini (1971) zonation:

Fig.5 Light microscope images of nannofossils from the Bid-Zard section southwest of Izeh,Zagros Basin.Taxonomic definition follows Perch-Nielsen(1985).a Pontosphaera multipora(Kamptner 1948 ex Deflandre in Deflandre and Fert 1954)Roth et al.1970;b Cribrocentrum isabellae Catanzariti,Rio and Fornaciari in Fornaciari et al. 2010; c Cyclicargolithus floridanus(Roth and Hay, in Hay et al. 1967)Bukry 1971; d Helicosphaera euphratis Haq 1966; e Sphenolithus dissimilis Bukry and Percival 1971; f Discoaster deflandrei Bramlette and Riedel 1954; g Coccolithus pelagicus(Wallich 1877) Schiller 1930; h Dictyococcites bisectus (Hay,Mohler and Wade 1966)Bukry and Percival 1971;i Braarudosphaera bigelowii(Gran and Braarud 1935)Deflandre 1947;j Coccolithus formosus(Kamptner 1963)Wise 1973;k Zygrhablithus bijugatus(Deflandre in Deflandre and Fert 1954)Deflandre 1959;l Helicosphaera bramlettei(Müller 1970)Jafar and Martini 1975;mSphenolithus pseudoradians Bramlette and Wilcoxon 1967; n Sphenolithus praedistentus Bramlette and Wilcoxon 1967; o-p Cribrocentrum reticulatum Perch-Nielsen 1971;q Discoaster barbadiensis Tan Sin Hok 1927;r Discoaster saipanensis Bramlette and Riedel(1954);s Sphenolithus distentus(Martini 1965)Bramlette and Wilcoxon 1967;t Reticulofenestra umbilicus(Levin 1965)Martini and Ritzkowski 1968.Scale bar for(a)to(t):5 μm

3.1 Biozones in the upper part of the Pabdeh Formation

3.1.1 Isthmolithus recurvus Zone/Sphenolithus pseudoradians Zone (NP19/NP20 combined zone) — Top of CNE18, CNE19, CNE20/CP15b

The lowest nannofossil biozone recorded in this study is the Priabonian Isthmolithus recurvus Zone/Sphenolithus pseudoradians Zone (NP19/NP20 combined zone). The NP19 zone is identified from the first occurrence of Isthmolithus recurvus to the first occurrence of Sphenolithus pseudoradians in the Martini (1971) zonation. Isthmolithus recurvus was not found in the studied section, so the lower boundary of the NP19 zone has not been detected. Subsequently, the NP20 zone is defined from the first occurrence of Sphenolithus pseudoradians to the last occurrence of Discoaster saipanensis in the Martini(1971) zonation. According to Martini (1971), since the first occurrence of S. pseudoradians was reported differently and might not be an accurate indicator for determining the lower boundary of the NP20 zone, the NP19 zone cannot be separated clearly from the NP20 zone.Hence, the interval below the last occurrence of Discoaster saipanensis (the marker of the base of both NP21 and CP16 zones), can be defined as the NP19/NP20 combined zone or is equivalent to the CP15b subzone.In this zone, the calcareous nannofossil species are abundant and they are well-preserved. The nannofossil assemblages include Blackites spinosus, Cribrocentrum isabellae, C. reticulatum, Coccolithus pelagicus, Cyclicargolithus floridanus, Discoaster deflandrei, D. barbadiensis, D. saipanensis, Ericsonia formosa (=Coccolithus formosus), Helicosphaera bramlettei, H. compacta, H.euphratis, Pontosphaera multipora, Reticulofenestra reticulata, R. stavensis, R. umbilicus, R. minuta, Sphenolithus praedistentus, S. pseudoradians, S. moriformis, and Zigrhablithus bijugatus. According to the Agnini et al.(2014) zonation, the CNE18 zone is defined as the interval from the top of common Cribrocentrum erbae to the first occurrence of Cribrocentrum isabellae. In the Bid-Zard section, the base of the CNE18 zone cannot be recognized due to the absence of Cribrocentrum erbae, but the top of the zone was determined by the first occurrence of C. isabellae (at 11.15 m). The CNE19 zone is defined from the first occurrence of C. isabellae (at 11.15 m) to the last occurrence of C. reticulatum (at 33.42 m), which is restricted to the Priabonian. In addition, the CNE20 zone is defined by the last occurrence of C. reticulatum, a suitable index used to determine the boundary of the CNE19/CNE20 combined zone (Fornaciari et al. 2010; Agnini et al. 2011). The last occurrence of C. reticulatum is found at the 33.42 m level and the last occurrence of D. saipanensis at 44.56 m, which respectively mark the base of the CNE21 zone and the top of the CP15b subzone in the Okada and Bukry (1980) zonation. Therefore, the NP19/NP20 combined zone corresponds to the upper part of CNE18,CNE19, CNE20,and CP15b with a thickness of 44.56 m.

3.1.2 Ericsonia subdisticha Zone (NP21) —CNE21-CNO1/CP16a-CP16b

The next nannofossil biozone recorded in the upper part of the Pabdeh Formation is the NP21 zone dated as Priabonian/Rupelian. The Ericsonia subdisticha Zone is defined as an interval from the last occurrence of

Discoaster saipanensis to the last occurrence of Ericsonia formosa (=Coccolithus formosus) (Martini 1971). The last occurrence of Discoaster saipanensis and the last occurrence of Ericsonia formosa were found at 44.56 m and 66.84 m, respectively. The most important nannofossil taxa of this zone belong to Coccolithus pelagicus,C. miopelagicus, C. formosus, Cyclicargolithus floridanus,Dictyococcites bisectus, Discoaster deflandrei, Helicosphaera bramlettei, H. compacta, H. euphratis, Pontosphaera multipora, Reticulofenestra minuta,Sphenolithus moriformis and Zygrhablithus bijugatus. It should be mentioned that across this interval, the total abundance of nannofossil increased and they are wellpreserved. In the Agnini et al. (2014) zonation, the CNE21 zone is introduced as an interval from the last occurrence of Discoaster saipanensis to the first occurrence of acme Clausicoccus subdistichus. Subsequently,CNO1 is defined as an interval from the first occurrence of acme Clausicoccus subdistichus to the last occurrence of Ericsonia formosa (=Coccolithus formosus). In the present study, the base of the CNO1 zone cannot be defined due to the absence of C. subdistichus. Therefore,the CNE21/CNO1 zone boundary and the CP16a/CP16b subzone boundary are not specified. The last occurrence of E. formosa (=C. formosus), marks the top of the CNO1 zone and the CP16b subzone (66.84 m). Hence,the NP21 zone corresponds to the CNE21-CNO1 combined zones and also to the CP16a-CP16b combined subzones. The thickness of this zone is 22.28 m.

3.1.3 Helicosphaera reticulata Zone (NP22) — CNO2/CP16c

The Helicosphaera reticulata Zone (NP22) is defined as the interval from the last occurrence of Ericsonia formosa(=Coccolithus formosus)to the last occurrence of Reticulofenestra umbilicus in the Martini (1971) zonation. In this study, the NP22 zone spans an interval from the last occurrence of C.formosus(at 66.84 m)to the last occurrence of Reticulofenestra umbilicus (at 89.12 m) and falls in the Rupelian(32.95-32 Ma)(Gradstein et al.2012).According to Martini(1971),the high occurrence of R.umbilicus defines the base of the NP23 zone. Species identified in this biozone include Braarudosphaera bigelowii, Coccolithus pelagicus, C. miopelagicus, Cyclicargolithus floridanus,Discoaster deflandrei, Dictyococcites bisectus, Helicosphaera euphratis, H. compacta, H. bramlettei, Reticulofenestra umbilicus,R.minuta,Sphenolithus praedistentus,S.moriformis,S.pseudoradians,and Zygrhablithus bijugatus.

The preservation of calcareous nannofossils in this zone is moderate to good.In the Agnini et al.(2014)zonation,the CNO2 and CNO3 zones are defined by the last occurrence of E.formosa(=C.formosus)and the last occurrence of R.umbilicus,which were found in the Bid-Zard section at 66.84 m and 89.12 m, respectively. Therefore, the NP22 zone corresponds to the CNO2 zone and CP16c subzone,already shown by Agnini et al. (2014). The NP22 zone reaches a thickness of 22.2 m.

3.1.4 Sphenolithus praedistentus Zone (NP23) —CNO3 and the lower part of CNO4/CP17 and the lower part of CP18

The uppermost nannofossil biozone recorded in the Pabdeh Formation is the Sphenolithus praedistentus Zone (NP23). According to Martini (1971), the NP23 zone spans an interval from the last occurrence of Reticulofenestra umbilicus to the first occurrence of Sphenolithus ciperoensis. As mentioned above, the last occurrence of Reticulofenestra umbilicus was recorded at 89.12 m, but the first occurrence of Sphenolithus ciperoensis was not confirmed due to the absence of all calcareous nannofossil species in the interval from 111.35 m to 118 m. Nannofossils were absent from the upper part of the Pabdeh Formation composed by the tufa microbial facies, which is a bio-event that coincided with the Rupelian/Chattian interval in the studied section. Therefore, the top of NP23 zone or the transition between the Pabdeh and Asmari formations is simply lacking in nannofossils. Here, the boundary between the formations as well as the upper part of the NP23 zone and the entire NP24 zone were not determined due to the lack of index species(Fig. 3). Species such as Braarudosphaera bigelowii,Coccolithus pelagicus, C. miopelagicus, Cyclicargolithus floridanus, Discoaster deflandrei, Dictyococcites bisectus, Helicosphaera bramlettei, H. euphratis, Reticulofenestra minuta, R. reticulata, Sphenolithus distentus, S.moriformis, S. praedistentus, S. pseudoradians and Zygrhablithus bijugatus were identified in the lower part of the NP23 zone, and were moderately preserved. In the Agnini et al. (2014) zonation, the CNO3 zone is defined from the last occurrence of Reticulofenestra umbilicus to the first occurrence of Sphenolithus distentus. In this study, the first occurrence of S. distentus was recorded at 105.85 m, so the boundaries of CNO3/CNO4 and CP17/CP18 zones were specified. According to Agnini et al. (2014), the CNO4 zone spans the interval from the first occurrence of Sphenolithus distentus to the last occurrence of S. praedistentus. As mentioned above, the last occurrence of S. praedistentus and the first occurrence of S. ciperoensis were not recorded. Therefore, the NP23 zone corresponds to the CNO3 zone and the lower part of CNO4 or the CP17 zone and the lower part of CP18. In this study, the abundance of calcareous nannofossils is also suddenly reduced in this biozone, and no nannofossil was recorded at levels between 111.35 m and 118 m. The lower part of the NP23 zone is defined as Rupelian (32-29.6 Ma)(Gradstein et al. 2012).

3.2 Biozone in the lower part of the Asmari Formation

Although calcareous nannofossils were absent from the transition zone of the Pabdeh and Asmari formations,gradually their abundance increased again in the shale of the Asmari Formation during the Chattian. The biozone identified close to the base of the Asmari Formation is the Sphenolithus ciperoensis Zone.

3.2.1 Sphenolithus ciperoensis Zone (NP25) — The lower part of CNO5/the lower part of CP19b

The last biozone in the section is the Sphenolithus ciperoensis Zone. It was proposed by Bramlette and Wilcoxon(1967) and emended by Martini (1976). According to Martini (1971), this zone is defined between the last occurrence of Sphenolithus distentus and the last occurrence of Sphenolithus ciperoensis. The nannofossil assemblage of the topmost part of the section includes

Braarudosphaera bigelowii, Coccolithus pelagicus, Cyclicargolithus floridanus, Pontosphaera enormis, Sphenolithus dissimilis, S. moriformis, and Zygrhablithus bijugatus. The last occurrence of S. distentus and the first occurrence of S. dissimilis were recorded at 123.5 m, but S. ciperoensis was not identified. Therefore, only the lower part of the NP25 zone was determined. The preservation of calcareous nannofossils in this part is moderate to good. Moreover, the lower part of this zone corresponds to the lower part of the CP19b subzone defined by the last occurrence of Sphenolithus distentus (at

123.5 m) according to the Okada and Bukry (1980) zonation. Based on Agnini et al. (2014), the CNO5 zone spans the interval from the last occurrence of Sphenolithus praedistentus to the last occurrence of Sphenolithus ciperoensis. At the base of the Asmari Formation,the last occurrence of S. predistentus was not recorded.However, based on the last occurrence of S. distentus,the lower part of the CNO5 zone is suggested. The age of this zone is Chattian (late Oligocene).

4 Discussion

The present study focuses on two key topics: (1) identification and abundance of calcareous nannofossils, and relevant bio-events; (2) environmental conditions and nutrient controls in the Eocene-Oligocene transition(EOT).

4.1 Evaluation of calcareous nannofossils and bio-events

Coccolithophores are the group of marine algae that secretes calcium carbonate crystals inside their cells to build individual coccoliths (Thierstein and Young 2004).This group is sensitive to changes in temperature, light,and nutrient concentration and is therefore widely used in paleoecological and environmental studies (Wade and Bown 2006). Also, they play an important role in the global carbon cycle (Rost and Riebesell 2004). Due to their great diversity and extensive geographical distribution in oceanic environments, coccoliths provide one of the most important fossil components of Cenozoic marine sediments (Bopp et al. 2001; Bown et al. 2004). According to Roth and Thierstein (1972), the preservation of species is highly sensitive to alterations in sediments such as recrystallization and dissolution. In other words,the preservation depends on the degree of etching and overgrowth, which provides accurate preservation information (Roth et al. 1970). In the present study, a semiquantitative analysis was carried out throughout the sequence from the late Eocene to the late Oligocene. Here,the calcareous nannofossils are generally common, as well as show good preservation so that many stratigraphic marker species have been recognized, allowing an accurate zonation. However, the diversity of calcareous nannofossils at the transition between the Pabdeh and Asmari formations sharply decreased. In this regard,the present study presents a biostratigraphic framework based on the occurrence of the following taxa: the first occurrences of Cribrocentrum isabellae, Sphenolithus distentus, and the last occurrences of Coccolithus formosus, Cribrocentrum reticulatum, Discoaster saipanensis,Reticulofenestra umbilicus, and Sphenolithus distentus.Species such as Coccolithus pelagicus, Cribrocentrum reticulatum, Cyclicargolithus floridanus, Dictyococcites bisectus, Discoaster deflandrei, Helicosphaera euphratis,Pontosphaera multipora, Reticulofenestra minuta, Reticulofenestra umbilicus, Sphenolithus moriformis, Sphenolithus praedistentus, and Zygrhablithus bijugatus are common in some slides. Discoaster barbadiensis, D. saipanensis, Helicosphaera compacta, Reticulofenestra stavensis, and Sphenolithus pseudoradians are relatively few in numbers. Moreover, Discoaster tanii, Sphenolithus distentus, and Sphenolithus dissimilis are often rare taxa in this study (Fig. 4). Most of the species listed above reflect open marine to shelf conditions (Okada and Bukry 1980; Perch-Nielsen 1985; Bohaty and Zachos 2003;Thierstein and Young 2004; Dunkley Jones et al. 2008;Agnini et al. 2014). In contrast, there are other indicators that indicate a decrease in water depth and stress conditions, including nannofossils such as Braarudosphaera(e.g.,Konno et al. 2007).

4.2 Environmental conditions and nutrient controls

About 34 million years ago, Earth’s climate shifted from relatively greenhouse conditions to icehouse conditions(e.g., Bohaty and Zachos 2003; Dunkley Jones et al.2008). A diverse range of fossil and sedimentary data change in a characteristic way during the Eocene-Oligocene transition, indicating a period of reorganization of the global climate (Blaj et al. 2009). For this period, calcareous nannofossils are precise ecological indicators that reflect environmental variables, as their distribution was affected by ecological parameters such as temperature and nutrient supply (Bown et al. 2004; Blaj et al. 2009). In fact, nannoplankton distribution pattern is strongly related to surface water-temperature and nutrient availability, which in turn are linked to water circulation. Tropical oceanic environments, which are seasonally stable and support the highest diversity but have few calcareous nannofossil assemblages, indicate oligotrophic (low nutrient) conditions (Kanungo et al.2017). In contrast, eutrophic (nutrient-rich) conditions that are associated with the continental shelf, support high abundances dominated by only a few species(Andruleit et al. 2003). Moreover, the recovery of hydrocarbons in Paleogene sediments and identification of source rocks, reservoir rocks, and stratigraphic traps have led to a demand for a detailed biostratigraphy of the Eocene-Oligocene interval (Kanungo et al. 2017).During this period, important Eocene-Oligocene genera include Chiasmolithus, Discoaster, Helicosphaera, Sphenolithus, and Reticulofenestra, which provide evidence for climate changes (Coxall and Pearson 2007; Bordiga et al. 2015). Significant changes in the facies near the boundary of the Pabdeh and Asmari formations in the Bid-Zard section are also due to environmental changes.In the present study, dolomitic facies with microbialites formed in the upper part of the Pabdeh Formation in very shallow water at the boundary of the Pabdeh and Asmari formations in the Oligocene (Rupelian). The ideal condition for tufa formation is water temperature,but some studies have shown that tufa also formed under lowered temperature conditions, which are influenced by seasonal temperature fluctuations (Dipova and Doyuran 2006; Kosun 2012). After the formation of microbial tufa in the upper part of the Pabdeh Formation,anhydrite sediments formed at the base of the Asmari Formation. Anhydrite deposits form in warm climates(Dipova and Doyuran 2006). Subsequently, during marine transgression, shales were deposited on the Asmari anhydrite. Therefore, the presence of evaporite deposits between the Pabdeh and Asmari formations indicates a change in climatic conditions with a decrease in the depth of the sedimentary basin. Due to the depth decrease in the late Rupelian, a disconformity was formed between the Pabdeh and Asmari formations. In addition,the existence of a disconformity was confirmed by a red paleosol horizon. The changing conditions also control the nannofossil assemblages at the boundary. In fact,there is a direct relationship between temperature fluctuations and calcareous nannofossil assemblages. According to Bown et al. (2004), cooling reduced the diversity of calcareous nannofossils during the Oligocene. In contrast, increasing productivity is suggested by the occurrence of calcareous nannofossils. For example,

Coccolithus pelagicus and Cyclicargolithus floridanus

Fig. 6 Comparison of the standard zonation schemes for calcareous nannofossils of the Tethyan realm with the biozonation established in the Bid-Zard section, and the correlation of bio-events in the transition between the Pabdeh and Asmari formations in the Zagros Basin (the red line is the boundary between the Pabdeh and Asmari formations)

occur in high-productivity eutrophic environments and indicate a marked increase of nutrient concentrations in surface waters during the Oligocene (Dunkley Jones et al. 2008; Villa et al. 2008; Blaj et al. 2009; Shcherbinina 2010). In this study, Cyclicargolithus floridanus increases slightly in abundance up-section and just above the Eocene-Oligocene Transition (EOT). Coccolithus pelagicus, another important component of the nannofossil assemblage (Fig. 4), increases in abundance between 22.29 m and 60.5 m, above the EOT, and then decreases from 60.5 m onwards. The presence of Cyclicargolithus floridanus, with an increase in the Rupelian interval, indicates eutrophic conditions and shows that the Oligocene had a cooler climate than the Eocene (Winter et al. 1994; Villa et al. 2008; Oszczypko-Clowes and ?ydek 2012). Moreover, temperate water taxa were also identified in the studied section based on the varying abundance of

Reticulofenestra umbilica and Dictyococcites bisectus (Villa et al. 2008; Bordiga et al. 2015), among which D. bisectus constitutes a significant part of the assemblages throughout the late Eocene (NP19/NP20 combined zone) to the early Oligocene (the lower part of NP23 zone). In contrast, the abundance of Discoaster species(D.barbadiensis,D.saipanensis, D. tanii) decreases below the EOT and never becomes as abundant as in the Priabonian. In addition, taxa such as Braarudosphaera and small Reticulofenestra spp.indicate shallow shelf environments (Kameo 2002; Villa et al. 2008). For example, Braarudosphaera bigelowii, observed in the upper part of the studied section, prefers coastal waters and nutrient-rich environments (Konno et al. 2007). The small Reticulofenestra spp. are eutrophic species, as well as other reticulofenestrids such as Reticulofenestra minuta, and would suggest warm marine conditions (Kameo 2002). Therefore, the abundance of Reticulofenestra minuta is connected with marine nutrient availability at the top of the Pabdeh Formation. Here, the mechanisms for increasing nutrient availability along with decreasing water depth in the Rupelian should be considered. Due to the diversity of species such as Cyclicargolithus, Discoaster, Helicosphaera, Pontosphaera, and Sphenolithus,which are known as very useful ecological indicators, it can be established that the study area was located at low latitudes. Also, helicosphaerids such as Helicosphaera euphratis, common in the Paleogene assemblages, are generally considered to prefer warm to temperate waters and are an indicator for increased nutrient availability (Auer et al. 2014). Thus, this also supports that the basin has been affected by environmental changes.

In the present study, the zonal biostratigraphy based on calcareous nannofossils was also compared with the existing zonations of the Zagros Basin such as the Fars Province (Khavari Khorasani et al. 2014) and the

Khozestan Province (Ahifar et al. 2015), with the standard Eocene-Oligocene zonation published by Perch-Nielsen (1985), and with the Paleogene zonations for the Tethyan realm (Agnini et al. 2014) (Fig. 6). The transition between the Pabdeh and Asmari formations indicates different ages in different parts of the Zagros Basin, allowing us to recognize several biozones within the studied time interval. Thus, the upper part of NP23 and NP24 was not recorded in the uppermost part of the Pabdeh Formation during the Oligocene and a disconformity formed at the boundary between the Pabdeh and Asmari formations. In contrast, in some of the previous studies of the Fars Province, a continuous boundary between the Pabdeh and Asmari formations was recorded from the Chattian (NP24) (e.g., Khavari Khorasani et al.2014), and in some previous studies of the Tang-e-Abolhayat section, a continuous trend was recorded from the Priabonian (NP20) to the Chattian (NP24). In the Khozestan Province (Hatti section), Ahifar et al. (2015)also delineated a continuous boundary between the Pabdeh and Asmari formations in the Rupelian(NP23/CP17).No biozone was recorded in the lower part of the Asmari Formation in the previous studies. In the present section,a biostratigraphic gap could be observed at the Rupelian-Chattian boundary, while the lower part of NP25/CP19b was recorded at the base of the Asmari Formation.

5 Conclusions

This study focused on the distribution of calcareous nannofossils and environmental variable conditions. As a result,29 species belonging to 13 genera of calcareous nannofossils were recorded from the Pabdeh and Asmari formations southwest of Izeh in the Zagros Basin of the eastern Tethys domain. Based on evolutionary changes among the calcareous nannofossils, five biozones were recognized according to standard zonation schemes, including the NP19/NP20 combined zone (CP15b/top of CNE18-CNE20) to NP23(CP18/CNO4) in the upper part of the Pabdeh Formation and the NP25 (CP19b/CNO5) biozone at the base of the Asmari Formation. The different zonations indicate a biostratigraphic gap and no continuous sedimentary record at the Pabdeh/Asmari boundary during the Oligocene, and a lack of nannofossils between the Rupelian and Chattian intervals.However,nannofossils reappeared in the lower part of the Asmari Formation. Biostratigraphic zonations confirm the dating of the studied interval as late Eocene (Priabonian) to late Oligocene (Chattian). Moreover, the presence of a diverse low-latitude nannofossil assemblage in the studied section indicates a tropical climatic.The analysis of calcareous nannofossil assemblages indicates that their abundance and diversity decreased towards the boundary between the Pabdeh and Asmari formations.This decrease is due to environmental changes.Thus,the calcareous nannofossils document relatively warm marine condition, associated with a decreased water depth, which prevailed during the late Eocene-Oligocene of the southwestern Izeh,Zagros Basin,of the eastern Tethys domain.

Abbreviations

A: Abundant; C: Common; CCD: Calcite compensation depth;CNE: Calcareous nannofossils Eocene; CNO: Calcareous nannofossils Oligocene; CP: Coccoliths Paleogene; F: Few; FO: First occurrence; FOV: Field of view; LO: Last occurrence; NP: Nanoplankton Paleogene; R: Rare

Acknowledgements

Both authors would like to thank the editors of the Journal of Palaeogeography for helpful comments on versions of the manuscript, and we thank anonymous reviewers for their help to improve the manuscript.

Authors’ contributions

SS analyzed and interpreted the data regarding the calcareous nannofossils and their biostratigraphy, and was a major contributor in writing the manuscript. FJ performed the sampling and preparation of the slides. Both authors read and approved the final manuscript.

Funding

This study has been supported by Grant 751541 of the Imam Khomeini International University.

Availability of data and materials

All datasets cited and analyzed in this study are available in this published paper.

Declarations

Competing interests

The authors declare that they have no competing interests.

Author details

1Department of Mining, Faculty of Engineering, Imam Khomeini International University, Qazvin, Iran.2Department of Geology, Faculty of Science,University of Birjand, Birjand, Iran.

Received: 10 November 2020 Accepted: 22 April 2021