Sequence Stratigraphy of the Late Cretaceous–Paleocene Gurpi Formation in Southwest Iran

GeoArabia, 2014, v. 19, no. 2, p. 89-102 Gulf PetroLink, Bahrain

Sequence stratigraphy of the Late Cretaceous–Paleocene Gurpi Formation in southwest Iran

Bijan Beiranvand, Ebrahim Ghasemi-Nejad, Mohammad Reza Kamali and Akram Ahmadi

ABSTRACT

Facies associations, microplanktonic diversity, palynofacies variations, geochemical data, and natural gamma-ray logs were analyzed from the Danial and Gurpi sections of the Campanian–Selandian Gurpi Formation in the Zagros Mountains, southwest Iran. The biostratigraphic data indicate that deposition across the Cretaceous/Paleogene boundary was continuous in the Danial Section. In contrast, a minor stratigraphic break seems to be present in the Gurpi Section, where several planktonic foraminiferal subzones are not identified. Nine depositional sequences were interpreted and correlated between the two sections. They are apparently of great lateral extent because they closely correlate to the global sea-level cycles. The Maastrichtian maximum flooding surface MFS K180 (68 Ma) of the Arabian Plate was also identified. Detailed palynofacies analysis, integrated with standard tropical/subtropical planktonic foraminifera, indicate warm Neo-Tethyan upper-bathyal to middle-shelf depositional environments for the Gurpi Formation.

INTRODUCTION

The Upper Cretaceous–Paleocene Gurpi Formation consists of deep-marine shales, marls and argillaceous lime mudstones that crop out in the Zagros Mountains, southwest Iran (Figure 1). The petroleum geology, stratigraphy and sedimentology of the formation have been presented in several papers (e.g. James and Wynd, 1965; Wynd, 1965; Sampo, 1969; Setudehnia, 1972, 1978; Stoneley, 1974, 1990; Motiei, 2003; Ghasemi-Nejad et al., 2006; Darvishzad et al., 2007). In contrast to previous studies, this paper seeks to integrate the biostratigraphy and facies analysis of the formation in a chrono- and sequence-stratigraphic framework. It is the third of a series that include Beiranvand and Ghasemi-Nejad (2013) and Beiranvand et al. (2013).

In the first paper, Beiranvand and Ghasemi-Nejad (2013) documented the high-resolution planktonic foraminiferal biostratigraphy of the Gurpi Formation in the Danial Section, located on the Payun Anticline, northeast Izeh in Iran (Figure 1). The study was focused on a 15 m stratigraphic interval that spans the Cretaceous/Paleogene boundary (K/Pg; see Table 1 for abbreviations), traditionally referred to as the Cretaceous/Tertiary boundary (KTB). This boundary is recognized as a major regional unconformity, referred to as the “pre-Cenozoic” or “base Tertiary” unconformity in the Arabian Plate (e.g. Sharland et al., 2001). In the Danial Section, however, a complete succession was shown to cross the boundary as evident by the occurrence of the uppermost Cretaceous Plummerita hantkeninoides Taxon-range Zone (TRZ), and lowermost Paleocene Guembelitria cretacea (P0) and Parvularugoglobigerina eugubina zones (see Enclosures 1 to 4).

The second paper of the series by Beiranvand et al. (2013) interpreted genetic stratigraphic sequences and sea-level cycles based on the response of palynomorphs in the complete section of the Gurpi Formation in the Danial Section. The present paper interprets the depositional sequences and sea-level cycles in the Gurpi Anticline, where the type section of the Gurpi Formation is defined (James and Wynd, 1965). It then shows that the nine Gurpi Formation sequences interpreted in the Danial Section can be correlated to the Gurpi Section over a distance of about 70 km (Figure 1).

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 Beiranvand et al.

48ºE 49º 50º51º 52º Mountain Bala Rud Gurpi Section Front Fault Fault Zone Danial Section Lali Dezful N Kal-Chenar H Imbricated Belt i 0 100 Shush g h Z km 32ºN Shushtar a 32º Pyun gr Dezful Embayment Masjed os t Izeh Faul Soleyman Zagros Main Thrust Fault Ahvaz (Ahwaz) Izeh IRAQ Fault Zone

Khuzestan Ramhormoz 31º Region 31º (Al Ahwaz) Imbricated IRAN Mountai Bandar Front Fault High ZagrosBel Faultt Mahshahr n Abadan

48º 49º 50º51º

Figure 1: Geographic location of the central Zagros Mountains and the two City studied sections. The Danial Section is located about 145 km northeast of Road Ahvaz city in southwest Iran. It was sampled at Pyun Anticline, about 25 km Thrust/reverse fault northeast of Izeh, and about 3 km north of the Kal-Chenar village (base Strike-slip fault located at 32°25'10"N and 49°58'04"E). The Gurpi Section, the type section of the Gurpi Formation, is located at 32°26'42"N and 49°15'54"E. Oil ﬁeld Gas ﬁeld

GEOLOGICAL SETTING

During the late Cretaceous and Paleocene, the northeastern region of the Arabian Plate was covered by the Neo-Tethys Ocean (Ziegler, 2001; Sharland et al., 2001; Motiei, 2003; Alavi, 2004). The region of the present-day Zagros Mountains was a NW-trending foredeep basin situated southeast of the Zagros Suture Zone (Figure 1; Ziegler, 2001; Alavi, 2004). The deepest part of the basin was located next to the suture and a broad shallow-marine shelf extended to the west over most of the Arabian Plate (Koop and Stoneley, 1982). During that time, an upwelling system of currents developed in the southeastern Neo-Tethys Ocean resulting in the deposition of thick organic-rich, relatively deep- marine shales and argillaceous lime mudstones in the Zagros Foredeep (Kolodny, 1980; Parrish and Curtish, 1982; Almogi-Labin et al., 1990, 1993; Eshet et al., 1994; Ziegler, 2001). The Gurpi Formation was deposited during the late Campanian to late Danian in this asymmetric foredeep basin southeast of the Zagros Suture (Bahroudi and Talbot, 2003).

The Gurpi Formation is about 340 m thick in the Gurpi and Danial sections (Figure 1, Enclosures 1 and 2), and consists of dark-gray to gray marly shales, marls and marly limestones. In these sections it overlies the Ilam Formation and the contact is a sharp disconformity surface with no evidence of subaerial erosion or weathering. It is conformably overlain by purple clayey shale and marl layers of the basal part of the Pabdeh Formation. The carbonates of the Emam Hassan Member occur in the middle part of the Gurpi Formation, and extend throughout the Zagros Basin with some distinctive lateral facies variations. In the study area, the member consists of thin-bedded pelagic and hemipelagic marls and argillaceous lime mudstones with planktonic foraminiferal assemblages (Wynd, 1965). This member has a gradational contact with the underlying strata making them practically indifferentiable except for their lithologic color.

STUDY PROGRAM

The study program started in the field by measuring the thickness of the two sections and recording their lithology, sedimentary structures and geometry. Natural gamma-ray (NGR) was

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 Stratigraphy of Gurpi Formation, Late Cretaceous–Paleocene, southwest Iran

Table 1 Abbreviations

AOM = amorphous organic matter OI = oxygen index Al = aluminium GR = gamma-ray PMI = palynological marine index Ba = barium HI = hydrocarbon index PRSZ = partial-range subzone Ca = calcium HST = highstand systems tract Rm = marine palynomorphs richness Cs = caesium IZ = interval zone Rt = terrestrial palynomorphs richness Fe = iron km = kilometer SB = sequence boundary Ir = iridium K/Pg = Cretaceous/Paleogene boundary S/C = spiniferites/cyclonephelium ratio K = potassium KTB = Cretaceous/Tertiary boundary SEM = scanning electron microscope Mg = magnesium Kyr = thousand years T-R = transgressive-regressive Mo = molybdenum LOZ = lowest-occurrence zone Tmax = maximum temperature Rb = rubidium LOSZ = lowest-occurrence subzone TOC = total organic carbon Si = silica LST = lowstand systems tract TRZ = taxon-range zone Sn = stannum m = meter TS = transgressive surface Th = thorium Ma = million years before present TST = transgressive systems tract Ti = titanium MFS = maximum ﬂooding surface XRD = X-ray diffraction U = uranium Myr = million year XRF = X-ray ﬂuorescence Zn = zinc NGR = natural gamma-ray Zr = zirconium

measured using a SURVEY mode in GR-130 portable gamma in the two sections (677 m thick) with a 0.25 m sampling interval (Enclosures 1 and 2). A total of 299 rock samples were collected from the sections, most of which consist of dark-gray to gray, laminated and slightly bioturbated marly shales, marls, and shaley marls with subordinate marly limestones.

A total of 180 samples were selected for thin-section preparation and 106 samples were processed for planktonic foraminiferal biostratigraphy, following the standard method of Keller et al. (1995). In addition, a total of 210 samples were processed for stratigraphic and organic-matter investigations, following the palynological technique described in Wood et al. (1996). A total of 15 samples were prepared for X-ray diffraction (XRD) analysis and 25 powdered bulk samples underwent standard analytical methods of Espitalie et al. (1986) using Rock-Eval II plus TOC (total organic carbon) module to determine TOC, Tmax (maximum temperature), OI (oxygen index), and HI (hydrocarbon index). The distribution of major and trace elements was determined on 52 samples by XRF (x-ray fluorescence) analysis.

Thin sections were studied by optical microscopy, and the data were processed using the routine grain-size distribution analysis (Flügel, 2004) for petrography and microfacies analysis so as to correlate the facies between the sections. Dating of time lines was established using a detailed planktonic foraminiferal biozonation and dinocyst events in the sections (Enclosures 3 and 4). X-ray techniques (XRD, XRF and SEM) were used as the main methods to identify mineral phases, clay content, chemical composition (major and trace element analysis), and microplanktonic species diversity. Palynofacies analysis, geochemical data (TOC, Tmax, OI, and HI), and natural gamma- ray spectroscopy, which are strongly dependent on variation in depositional environments (Tyson, 1995), were used as multi-desciplinary approaches to interpret sea-level cycles. The interpretation of this data is described in more detail below.

Facies Analysis

Field descriptions of lithology, sedimentary structures, trace fossils and strata surfaces, together with microfacies analysis of the formation, led to the identification of four main sedimentary facies related to sea-level variations.

(1) Gray marls and argillaceous limestones to carbonate mudstone/wackestone containing minor debris of planktonic and rare benthic foraminifera; this facies occurs in shallowing upward parasequences of the lowstand systems tract (LST) of each depositional sequence. It is consolidated following a period of erosion and/or non-deposition after the falling stage of the sea-level cycle.

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 Beiranvand et al.

(2) Dark-gray to gray marls and marly shales to carbonate wackestone, wackestone-packstone and packstone, with abundant planktonic and rare benthic foraminifers; this facies occurs in deepening upward parasequences of the transgressive systems tract (TST) during the rising stage of the sea-level cycle. (3) Gray fossiliferous marls and shaley marls containing abundant planktonic and rare benthic foraminifera occur in the parasequences of the early highstand systems tract (HST). (4) Light-gray to cream marls interbedded with marly limestones, characterized by planktonic foraminifera wackestone, containing common trace fossils with abundant oxidized minerals suggest deposition during period of low sea level within outer-neritic to upper-slope settings.

A hardground surface at the top of each depositional sequence marks a major sedimentological change from the gray marly limestones and marls below, to the marls and clayey marls above. The latter contains several resistant layers of cream marly limestones. This pattern suggests a rising sea level, sometimes with fluctuations.

Planktonic Foraminifera

Nutrient supply, oxygen, temperature and salinity are the main environmental parameters that control variations in the relative abundances of planktonic foraminifers’ species (Huber et al., 1999; Gebhardt et al., 2004). It has also been demonstrated that the eustatic sea level can fluctuate with temperature. Therefore, shifts in the percentage and morphotype groups of planktonic foraminifera indicate changing environmental conditions, which can be linked to water depth (e.g. Boersma and Shackleton 1981; van der Zwaan et al., 1990; Huber et al., 1995, 1999; Nederbragt et al., 1998; West et al., 1998; Premoli Silva and Sliter, 1999; Gebhardt et al., 2004). Accordingly, four main planktonic foraminiferal morphogroups are recognized:

(1) Heterohelicids are interpreted as opportunists indicating high productivity or unstable conditions (Nederbragt et al., 1998; West et al., 1998; Premoli Silva and Sliter, 1999), which may point to the deepest habitats (Huber et al., 1995, 1999). (2) Hedbergellids (non-keeled group) are interpreted as surface dwellers (Huber et al., 1995). These are considered to be open-marine species, which dominate when favorable conditions at greater depths do not exist (shallow-water depths or oxygen minimum zones; Leary et al., 1989; Koutsoukos and Hart, 1990; West et al., 1998; Premoli Silva and Sliter, 1999). (3) The Whiteinella and Praeglobotruncana group indicate intermediate water-depth habitats (Huber et al., 1995, 1999). (4) The Rotalipora and Dicarinella (keeled genera) group are interpreted as deep open-marine dwelling species (Huber et al., 1995, 1999; Premoli Silva and Sliter, 1999), and can represent different ecological niches. These are incorporated here for environmental interpretation.

Thus, high percentages of planktonic foraminifera, with frequent keeled specimens, are characteristic of deep water. The percentages within the four main environmental index groups indicate values vary widely between 0% and 61%, with a generally increasing trend according to the environmental habitat of each group during sea-level cycles.

Palynomorphs

The main palynological parameters that can be used to interpret relative sea-level changes include:

(1) amount of terrestrially derived organic matter (Posamentier and Vail, 1988; Roberts and Coleman, 1988; Leithold and Bourgeois, 1989); (2) marine to terrestrial palynomorph ratio (Gregory and Hart, 1992); (3) the ratio of recycled palynomorphs (Eshet et al., 1988a, b); (4) depositional organic facies (Habib and Miller, 1989; Habib et al., 1992; Edet and Nyong, 1993); (5) dinoflagellate species diversity (Habib and Miller, 1989; Habib et al., 1992; Moshkovitz and Habib, 1993).

In this study, the abundance and diversity Phytoclasts of marine palynomorphs, represented by 100 0

the Palynological Marine Index (PMI), and Redox plus proximity to ﬂuvial source s 90 10 Spiniferites/Cyclonephelium Ratio (S/C) are I used as the main palynological parameters 80 20 to determine changes in depositional 70 II 30

environments and transgressive-regressive 60 40 sequences (Enclosure 1). III 50 IVa 50 The percentages of five main palynological 40 VI IVb 60 elements (1) amorphous organic matter (AOM), 30 70 (2) dinocyst, and (3) foraminiferal test linings 20 IX 80 s (as marine palynomorphs), (4) phytoclast and Redox plus proximity to ﬂuvial source VII V 10 90 (5) sporomorphs (as terrestrial components) VIII vary from 0–96%. They show a generally 0 100 increasing trend for marine elements from the 100 90 80 70 60 50 40 30 20 10 0 transgressive surface (TS) to the maximum AOM Redox plus masking effect Palynomorphs flooding surface (MFS), and a decreasing trend for phytoclast and sporomorphs in Lower Gurpi Fm Upper Gurpi Fm (Cretaceous) each sequence. The calculated percentage of Middle Gurpi Fm Upper Gurpi Fm (Paleocene) palynomorphs is largly controlled by AOM and Figure 2: Tyson (1995) “AAP” ternary diagram phytoclast dilution. Large percentages of AOM representing concentrations at IX and VIII in the succession resulted from a combination zones that show deep basin and stratiﬁed shelf of good preservation (directly related to sea deposits of the Gurpi Formation. dysoxic–anoxic conditions) and low-energy environments.

Low percentages of terrestrial elements (phytoclasts and sporomorphs) are mostly related to deep- marine depositional conditions and long-distance transportation of the particles. A very small ratio of sporomorphs to phytoplankton reflects the proximal–distal trend and indicates deep-marine condition for the succession. Plotting these data on ‘AAP’ ternary diagram of Tyson (1995; Figure 2) shows a concentration in the IX and VIII zones implying deep-water basinal setting and stratified marine shelf deposits. On the other hand, determining species diversity of dinocysts indicated that values vary widely between 1 and 42 with a normal increasing trend during relative sea-level rises. Generally, the abundance and diversity of marine palynomorphs increase during deposition of TSTs and the highest diversity is recorded at the MFSs (Enclosure 1).

The PMI that was proposed by Helenes et al. (1998) to support the interpretation of depositional environments is calculated using the formula:

Palynological Marine Index (PMI in %) = 100 x [(Rm/Rt) + 1]

Rm is richness of marine palynomorphs (dinoflagellates cysts, acritarchs and foraminiferal test linings) and Rt is the richness of terrestrial palynomorphs (pollens and spores) counted per sample. In the present study, Rm and Rt were expressed as number of species per sample. High values of PMI are interpreted as indicative of normal marine depositional conditions. When the samples have no marine palynomorphs (Rm = 0), the PMI value is 100. The results of determining PMI from counting 210 samples are shown in Enclosure 1. PMI values range widely between 167 and 1,600 with a generally increasing trend from the transgressive surface to the MFSs (Enclosure 1).

Generally, cyst types with different lengths and complexity of processes can be used in environmental analysis (Williams, 1977). The association with dominant long and complex process chorate cysts (such as Spiniferites, Achomosphaera, Oligosphaeridium, Hystrichosphaeridium, and Hystrichodinium) is considered to indicate open-marine shelf environments (Davey and Rogers, 1975; Wall et al., 1977; Brinkhuis and Zachariasse, 1988; Carvalho et al., 2006), and the association of proximate cysts with short, stout and berbed-process (such as Cyclonephelium, Exochosphaeridium, Cleistosphaeridium, and Micrhystridium) reflects a more coastal to nearshore relatively restricted

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 Beiranvand et al.

environment (Brinkhuis and Zachariasse, 1988; Eshet et al., 1992). Optima of the Spiniferites group are in most cases associated with a sea-level high and/or an energy low condition. However, the S/C ratio, as additional evidence of the palynologic response to sea-level change, can be defined and taken to reflect a more detailed, smaller-scale change of depositional environments compared to species diversity alone. In general, the S/C ratio shows an increase in the seaward direction and is a good indicator of environmental change in the offshore direction (Harland, 1973). The calculated ratio for 180 samples show a wide range between 0–10, with a typical increasing trend from the TS to the MFS (Enclosure 1). The subsequent change (increase, and then decrease) of the S/C ratio, corresponding with a similar change of species richness, highlights the T-R depositional sequences.

The results of palynofacies analysis indicate that dinoflagellate abundance and diversity, PMI, S/C ratio, and other sedimentary organic-matter parameters are related to the migration of the various depositional environments in response to sea-level changes. They contribute to the recognition of systems tracts and key horizons such as flooding surfaces and sequence boundaries (SB). Accordingly, each palynological cycle develops during a T-S sequence; it starts with a dramatic increase in PMI index, species diversity and S/C ratio in the TST, and the decrease of these criteria during the HST, thus reflecting a complete cycle of relative sea level.

Geochemical Analysis, Natural Gamma-ray Log and TOC Content

Inorganic whole rock geochemical data can be used to define depositional facies, stratigraphic surfaces, cycles and sequences reflecting changes in relative sea level. Nine major elements (Ca, Si, Al, Fe, Mg, K, Ba, Rb, and Ti) and six trace elements (Ir, Mo, Zn, Sn, Cs, and Zr) were determined

in this study. Data are reported as weight percent oxides (CaO, SiO2, Fe2O3, Al2O3, K2O, MgO, BaO,

RbO, and TiO2) of major and trace elements. The CaO compound shows the most abundance in the

successions (46.8 wt% on average of all samples) and SiO2 and Al2O3 are the other main components

with 23.6 wt% and 11.2 wt%, respectively. Others, Fe2O3, K2O, MgO, BaO, RbO, and TiO2 have

positive correlations with Al2O3 content. SiO2 contents in hemipelagic marine sediments consist of

both biogenic and detrital silica. The SiO2/Al 2O3 ratio has a clear negative correlation with Al2O3 content (Enclosure 1).

The samples with high Si/Al ratios, representing high biogenic silica, correspond to those with high TOC content, thus suggesting that TOC content should mainly be derived from sea- surface production of siliceous phytoplankton dominated by diatoms. In this study multi- element geochemistry has been used as a supplement for stratigraphic studies. The part of chemostratigraphy in Enclosure 1 shows the variations of the principal components versus depth to help discriminate between the different chemical signatures of the nine sequences.

The total of natural gamma-ray (NGR) intensity is a function of the combined contributions of K, U, and Th in sediments, matrix density, and matrix lithology. Clay mineral content is often diluted by other components such as biogenic silica. Because of this dilution, if the NGR signal in the sedimentary sequence is to be used to reconstruct the environmental record, it is important to know which chemical components relate to the NGR signal. Thus, NGR measurements can help determine the mineralogy and abundance of clay and other radioactive minerals including micas and feldspars (Serra, 1985). They can also help characterize depositional environments and diagenetic processes in sediments. The NGR intensity of the succession has been plotted with analyzed chemical results in Enclosure 1. The changes of U/Th ratio within the depositional sequences are indicative of the T-R sequences. Changes with depth occur where the parent rocks and organic carbon affect the radioactive element concentrations in specific layers.

Total organic carbon (TOC) in the sediments may be interpreted as a proxy for organic influx. In general, increase in TOC correlates with increasing water depth, but here the organic carbon and Rock-Eval pyrolysis data (Enclosure 1) indicate that TOC values are generally low in the TST and higher values (generally associated with hardgrounds) occur in the HST. The coincident increase in carbonate and proportion of planktonic foraminifera is related to the planktonic foraminifera (mainly heterohelicids) being the main sources of carbonate in the organic-rich shales and marls (Leine, 1986). The pyrolysis data further suggest that all of the kerogens are thermally immature,

Table 2 Correlation of Gurpi and global sequences

Gurpi Sequence (this study) Sneddon and Liu (2011)

Boundary and Type Section Danial Section Rate Surface Age Sequence (meter) (meter) meter/kyr

Top Gurpi SB Th1 59.3 Gu9 27 20 1.5 MFS Sel1 60.6 SB Gu9 SB Sel1 61.1 Gu8 20 14 1.5 MFS Da4 61.7 SB Gu8 SB Da4 62.1 Gu7 22 18 1.0–2.0 MFS Da2 64.5 SB Gu7 SB Da1 65.1 Gu6 50 33 1.0–8.0 MFS Ma5 66.3 SB Gu6 SB Ma5 67.7 Gu5 58 24 1.0–2.0 MFS Ma4 68.0 SB Gu5 SB Ma4 69.0 Gu4 68 94 8.0–9.5 MFS Ma2 69.7 SB Gu4 SB Ma2 70.1 Gu3 28 28 2.5–3.5 MFS Ma1 70.5 SB Gu3 SB Ma1 70.8 Gu2 33 56 3.0–3.5 MFS Cam9 72.1 SB Gu2 SB Cam9 73.1 Gu1 32 40 1.5–2.0 MFS Cam8 74.5 SB Gu1 SB Cam8 75.6

as indicated by the low to slightly elevated temperature maximum (Tmax) values of 422–457˚C and a terrigenous source (Type III). However, these data indicate that total organic carbon (TOC) values are generally very low and rarely exceed 0.5 wt%, with average values about 0.23 wt% (Enclosure 1).

SEQUENCE STRATIGRAPHY

Haq et al. (1988), in their “Mesozoic–Cenozoic Cycle Chart”, interpreted Campanian–Maastrichtian UZA-4 and Maastrichtian–Paleocene TA-1 supercycles, and dated them between 80.0 and 58.5 Ma in their time scale. They divided these supercycles into cycles or third-order sequences. On the basis of biostratigraphy the Gurpi Formation correlates to third-order sequences UZA-4.4 (75–71 Ma), UZA-4.5 (71–68 Ma) and TA-1.1 to TA-1.3 (68–60 Ma). The cycles of Haq et al. (1988) were revised by Hardenbol et al. (1998) and compiled and recalibrated in the revised geological time scale GTS 2004 (Gradstein et al., 2004) by Snedden and Liu (2011). In the latter updated chart, the Gurpi Formation is calibrated between 75.6 and 59.3 Ma, and correlates to sequences Campanian Ca8 and Ca9, Maastrichtian Ma1 to Ma5, Danian Da1 to Da4, and Selandian Sel1. The sea-level cycles and coastal onlap curves for these sequences are shown in Enclosure 1.

Beiranvand et al. (2013) interpreted the Gurpi Formation in the Danial Section in terms of nine depositional sequences, denoted Gu1 to Gu9 in ascending order (Enclosure 1). The nine Gurpi sequences are interpreted based on facies associations, palynofacies variations, microplanktonic diversity, and geochemical analysis. In the Danial Section, no evidences of subaerial erosion or unconformities is recognized across the Cretaceous/Paleogene boundary, and in general the outer- shelf and deeper-marine marls contain typical correlative conformities (Enclosure 2). Eustasy is believed to have controlled the sea-level cycles, and the sequences are characterized by sequence boundaries (SB), lowstand, transgressive and highstand systems tracts (LST, TST and HST), and maximum flooding surfaces (MFS). The sequence boundaries are given the same name as the overlying sequence (Table 2).

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 Beiranvand et al.

The ages of the sequences are determined from high-resolution planktonic foraminiferal biostratigraphy (Beiranvand and Ghasemi-Nejad, 2013; Enclosures 3 and 4), and correlated to the global sequences by stratigraphic position (Snedden and Liu, 2011). In Table 2, the correlations are listed together with the ages for SBs and MFSs of Snedden and Liu (2011). The sediment accumulation rate is calculated using the thickness of the sequence and depositional age of the global sequence. Note the K/Pg boundary is calibrated at 65.5 Ma in Snedden and Liu (2011), and 66.0 Ma in GTS 2012. In our study the K/Pg is interpreted at the top of Sequence Gu6.

Gurpi Sequence Gu1

SB Gu1 is taken at the base of the studied sections at the boundary between the Gurpi and Ilam formations. It occurs in the upper part of the Radotruncana calcarata TRZ and is correlated to SB UZA-4.4 and correlative SB Cam8, which is a major eustatic sea-level drop in the late Campanian (Haq et al., 1988; Hardenbol et al., 1998). Sequence Gu1 is correlated to Campanian Sequence Cam8 between 75.6–73.1 Ma (Snedden and Liu, 2011).

The TST Gu1 is characterized by decreasing in relative abundance of phytoclast particles and increased abundance of AOM directly related to dysoxic–anoxic conditions and low-energy environments (Enclosure 1). The increased abundance of heterohelicids and keeled genera planktonic foraminifera is interpreted as an open-marine deep-water habitat. Rock-Eval pyrolysis data indicate that TOC values are generally low and kerogen type is dominantly Types II and III (Enclosure 1). Geochemical analysis indicate that the amount of terrigenous minerals (including K, Al, and related elements) increase progressively during the TST. In addition, the U/Th ratio changes are directly related to the sea-level change. These criteria apparently reflect a rapid marine TST during a sea-level rise in outer-neritic conditions, or hydrodynamically lower-energy settings.

The MFS of Sequence Gu1 is identified by the maximum abundance of heterohelicids and keeled planktonic foraminifera, optimum species diversity, PMI, AOM, CaO components and highest ratios of S/C and U/Th, and scarcity of phytoclast particles. Additionally, GR and NGR logs show the greatest deflections at the MFS. This late Campanian MFS is positioned near the top of the Radotruncana calcarata TRZ, immediately below the base of the Globotruncanella havanensis PRZ. The dominance of carbonate in these beds reflects maximum accommodation and highest relative sea level, with minimum siliciclastic sediment influx in an overall shelf depositional setting. The HST exhibits a significant and progressive decrease in planktonic foraminiferal percentages and increase in coarser siliciclastic sediments, indicating an overall decrease in water depth and progradation.

Major element compounds, except for SiO2, show decreasing trends, and the GR and NGR logs deflect to the left. The subsequent decreases of U/Th ratio indicate the regression.

Campanian Gurpi Sequences Gu2 and Gu3

Dinocyst-based events and associated species, in correlation with the planktonic foraminiferal biozones, imply a Late Campanian age for the sequences Gu2 and Gu3, which are correlated to sequences Cam9 and possibly Ma1 (Snedden and Liu, 2011). Sequence Gu3 is apparently coincident with the brief Ma1 sea-level cycle (ca. 1 Myr) and may have been deposited during a lowstand near the Campanian/Maastrichtian boundary (e.g. Haq et al., 1988; Barrera et al., 1997; Li et al., 1998).

The TSTs of both sequences are characterized by low TOC and increasing PMI and microplanktonic diversity. Their MFSs are characterized by an abrupt change in trend of S/C and U/Th ratios, an increasing AOM, and relative abundance and diversity of planktonic foraminifera and dinocysts (100 dinocyst species belonging to 66 genera). The HSTs are reflected by a marked progressive decrease in planktonic foraminifera and dinocycst diversity and an increase in coarser phytoclasts, which indicate decreasing of water depths during regression and progradation (Enclosure 1). All the evidence indicates outer-neritic, deep-marine conditions for these two sequences.

Maastrichtian Gurpi Sequence Gu4

Sequence Gu4 consists of gray marl and marly limestone of the Emam Hassan Member (Enclosures 1 and 2). A total of 61 dinoflagellate cyst species belonging to 53 genera reveals a moderate diversity in the sequence. Based on the bioevents, it has a Maastrichtian age and is correlated to sequences Ma2 and Ma3 (Snedden and Liu (2011). It has the highest rate of sedimentation, 8–9.5 cm/Kyr, of all the Gurpi sequences (Table 2). SB Gu4 corresponds to a minimum of the U/Th ratio and occurs in the lower part of the Gansserina gansseri IZ. The TST occurs in the same interval zone and the increase in water depth is recognized by the increased clay content in the marls, and the abundance and diversity of the microfossil assemblages (Enclosure 1). The MFS occurs in the middle part of the Emam Hassan Member. In the Gurpi Section the MFS is often only recognized by a change in the relative abundance and diversity of planktonic foraminifera and dinocysts. The early HST is characterized by a sharp decrease in PMI values. Light gray to cream marls with carbonate planktonic foraminifera wackestone contain abundant oxidized minerals, suggesting deposition during low sea level within outer neritic to upper slope. The late HST reflects the shallowest environments and the greatest regression.

Late Maastrichtian Gurpi Sequence Gu5

Sequence Gu5 consists of gray marls and marly shales (Enclosures 1 and 2). Based on the planktonic foraminiferal biozones (CF5 and CF4) and correlatable dinoflagellate cysts bioevents, it is mid-Late Maastrichtian and correlated to Sequence Ma4 (Snedden and Liu, 2011). SB Gu5 is characterized by a minimum of the U/Th ratio, a sharp decrease in palynological content and a hardground. The LST is ca. 3 m thick and likely to be equivalent to the uppermost Contusotruncana contusa IZ The TST is recognized as an increase in water depth characterized by the vertical lithologic changes and the components of the microfossil assemblages; it occurs in the Pseudotextularia intermedia PRZ. The MFS corresponds to the base of the Abathomphalus mayaroensis TRZ, and is correlated to MFS K180 at 68 Ma (Sharland et al., 2001). The overlying sediments are interpreted as a relatively long-lasting HST with drops in dinoflagellate species diversity and an increase in planktonic foraminifera assemblage diversity. The paleoenvironmental indicators represent open-marine, outer-neritic conditions.

Late Maastrichtian Gurpi Sequence Gu6

Sequence Gu6 is the youngest Cretaceous sequence with the K/Pg boundary defining its upper boundary (Enclosures 1 and 2). It comprises dark gray marly shales, and planktonic foraminiferal biozones and dinoflagellate cysts bioevents indicate a latest Maastrichtian age. Sequence Gu6 is correlated to the late Maastrichtian Sequence Ma5 (Snedden and Liu, 2011). Note that the latter authors calibrated the age of 65.1 Ma for the K/Pg boundary as in GTS 2008; however in GTS 2012 it has an estimated age of 66.0 Ma.

SB Gu6 is characterized by a minimum of the U/Th ratio, and a sharp decrease in palynological content. It occurs at the base of the Pseudoguembelina hariaensis PRSZ. The LST is ca. 2 m thick, and probably occurs within the basal Pseudoguembelina hariaensis PRSZ. The TST occurs within this PRSZ and is characterized by a slight increase in the values of the PMI reaching the highest PMI and diversities in its upper part. The MFS is recognized by increased planktonic foraminifera and dinocysts of moderate diversity (79 dinocyst species belonging to 71 genera) and abrupt increasing of S/C, U/Th ratios, and amount of AOM. The HST corresponds to the latest Maastrichtian Pseudoguembelina palpebra PRSZ and Plummerita hantkeninoides TRSZ. It may contain a minor regressive phase (sea-level fall) in its middle part, indicated by the low percentage of planktonic foraminifera and slight increase in phytoclasts. At that time, the water depth reached a level where the oxygen minimum zone almost covered the area and permitted the preservation of organic matter. Other minor sea-level falls thought to have occurred at this time are, however, not evidenced by the presence of corresponding ‘‘shallow-water’’ benthic foraminifera. The paleoenvironmental evidence indicates open-marine outer-neritic conditions.

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 Beiranvand et al.

Unlike the Danial locality, where the section is apparently biostratigraphically continuous across the K-Pg boundary, the very latest Maastrichtian and very earliest Paleocene planktonic foraminiferal biozones (the P. hantkeninoides, G. cretacea, and P.eugubina zones) are apparently absent in the Gurpi (northeast Lali) locality. The time period represented by all these zones is only approximately 750,000 years in total and therefore the absence of any or all of them may be due to insufficiently closely spaced sampling.

Paleocene Gurpi Sequence Gu7

Sequence Gu7 is the oldest Paleocene sequence, and is characterized by the highest PMI values, particularly in the lower part, and lowest carbonate content. SB Gu7 is correlated to the K/Pg boundary and SB Da1 (Li et al., 1999; Snedden and Liu, 2011). The LST is ca. 1.5 m thick, and corresponds to the basalmost Paleocene planktonic foraminifera Parvularugoglobigerina eugubina TRZ (Pa) and Guembelitria cretacea PRZ (P0). These two foraminiferal biozones coincide with the lowest occurences of Damassadinium californicum and Senoniasphaera inornata dinocyst events that are indices for the K/Pg boundary. In addition to the bioevents, there is a correlation between the species richness, iridium anomaly, TOC, quartz, and calcite components (Enclosure 1) that also represent the K/Pg boundary.

The early Danian TST ends at the MFS, which is marked by an increase in frequency and diversity of planktonic foraminifera and dinocyst species (119 dinoflagellate cyst species belonging to 83 genera), and abrupt increase of PMI and AOM amounts and S/C and U/Th ratios. The HST corresponds to the Parasubbotina pseudobulloides PRSZ (P1a), and Subbotina triloculinoides LOSZ (P1b). It contains abrupt decreases in the abundance and diversity of biota and periodic increases in phytoclasts, which reflect minor sea-level falls. The highest TOC concentration occurs in this systems tract and is probably due to minor sea-level falls because of the occurrence of type III kerogen in these sediments. These minor falls are, however, not evidenced by the presence of corresponding ‘‘shallow-water’’ benthic foraminifera. The bathymetric investigations indicate a marine, mainly outer-neritic, with occasional transgressions, to upper-bathyal conditions for the sequence.

Paleocene Gurpi Sequence Gu8 and Gu9

These two sequences comprise dark gray to reddish gray marly shales. Based on the bioevents they are dated as late Danian and Selandian and correlated to sequences Da4 and Sel1 (Snedden and Liu, 2011). The two sequences are characterized by the highest PMI values, particularly in their lower parts. Their MFSs are marked by an increase in frequency and diversity of planktonic foraminifera and dinocyst species and abrupt increases of the PMI and AOM, and S/C and U/Th ratios. The HSTs show probable minor sea-level falls indicated by the abruptly decreasing abundance and diversity of biota and periodic increases of phytoclasts. The most abundant TOC occurs in Sequence Gu9 and is concentrated in the HST. Bathymetric investigations indicate a marine, mainly outer- neritic, with occasional transgressions to upper-bathyal conditions.

CONCLUSIONS

Four main lithofacies, integrated with biofacies response to paleoenvironmental changes, indicate that the Upper Campanian–Paleocene Gurpi Formation consists of at least nine depositional sequences. These were deposited in the Zagros Foredeep in settings that varied between middle- to outer-neritic and particularly upper-bathyal depths and hydrodynamically lower-energy conditions. The Gurpi sequences can be correlated to the global sequences compiled by Snedden and Liu (2011).

The Cretaceous/Paleogene (K/Pg) boundary, traditionally referred to as the Cretaceous/Tertiary boundary (KTB), is recognized as at the boundary between Gurpi sequences Gu6 and Gu7. The lowstand systems tract of basalmost Paleocene Gurpi Sequence Gu7 is characterized by planktonic foraminifera Parvularugoglobigerina eugubina TRZ (Pa) and Guembelitria cretacea PRZ (P0), the

lowest occurence of Damassadinium californicum and Senoniasphaera inornata dinocyst events, and a correlation between the species’ richness, iridium anomaly, TOC, quartz, and calcite components – all representataive indices for the K/Pg boundary (Enclosure 1).

Arabian Plate maximum flooding surface Maastrictian MFS K180 (Sharland et al., 2001) is a distinctive mappable surface with an estimated age of at 68 Ma. It coincides to the base of the Abathomphalus mayaroensis TRZ in Gurpi Sequence Gu6, which is correlated to global Maastrichtian Sequence Ma5 (Snedden and Liu, 2011).

The low TOC content of the Gurpi depositional sequences in the TSTs and higher values in HSTs can be interpreted as due to the floodings of planktonic foraminifera (mainly heterohelicids) as the main sources of carbonate in the organic-rich shales and marls. Furthermore, the pyrolysis-data suggest that all of the kerogens present here are thermally immature, as indicated by the low to slightly elevated temperature maximum (Tmax) values of 422–457°C.

ACKNOWLEDGEMENTS

This study was supported by the Research Institute of Petroleum Industry (RIPI) of Iran. We thank Mr. S.S. Hendi, Head of Exploration and Production Research Center for supporting the work and Mrs. E.H. Tavakoli and H. Alinaghian for processing the samples. The authors thank Roger Davies, Mike Bidgood and Moujahed Al-Husseini for their careful review and many constructive remarks that greatly improved the quality of this paper. Kathy Breining is thanked for proofreading the text and GeoArabia’s Production Co-manager, Arnold Egdane, for designing the paper for press.

REFERENCES

Alavi, M. 2004. Regional stratigraphy of the Zagros fold-thrust belt of Iran and its proforeland evolution. American Journal of Science, v. 304, p. 1-20. Almogi-Labin, A., A. Bein and E. Sass 1990. Agglutinated foraminifera in organic-rich neritic carbonates (Upper Cretaceous, Israel) and their use in identifying oxygen levels in oxygen-poor environments. In C. Hemleben, M.A. Kaminski, W. Kuhnt and D.B. Scott (Eds.), Palaeoecology, Biostratigraphy, Palaeoceanography, and Taxonomy of Agglutinated Foraminifera. Kluwer, Dordrecht, p. 565-585. Almogi-Labin, A., A. Bein, and E. Sass 1993. Late Cretaceous upwelling along the Southern Tethys margin (Israel): Interrelationship between productivity, bottom water environments and organic matter preservation. Paleoceanography, v. 8, p. 671-690. Bahroudi, A. and C.J. Talbot 2003. The configuration of basement beneath the Zagros Basin. Journal of Petroleum Geology, v. 26, no. 3, p. 257-282. Barrera, E., S.M. Savin, E. Thomas and C.E. Jones 1997. Evidence for thermohaline-circulation reversals controlled by sea- level change in the latest Cretaceous. Geology, v. 25, p. 715-718. Beiranvand, B. and E. Ghasemi-Nejad 2013. High resolution planktonic foraminiferal biostratigraphy of the Gurpi Formation, K/Pg boundary of the Izeh Zone, SW Iran. Revista Brasileira de Paleontologia, v. 16, no. 1, p. 5-26. Beiranvand, B., E. Ghasemi-Nejad and M. Reza Kamali 2013. Palynomorphs’ response to sea-level fluctuations: A case study from Late Cretaceous–Paleocene, Gurpi Formation, SW Iran. Journal Geopersia, v. 3, no. 1, p. 11-24. Boersma, A. and N.J. Shackleton 1981. Oxygen- and carbon-isotope variations and planktonic foraminifer depth habitat, Late Cretaceous to Paleocene, Central Pacific. In J. Thiede, T.L. Valuer, et al., Initial Reports of the DSDP 62, Washington (U.S. Govt. Printing Office), p. 513-526. Brinkhuis, H. and W.J. Zachariasse 1988. Dinoflagellate cysts, sea level changes and planktonic foraminifers across the Cretaceous-Tertiary boundary at El Haria, northwest Tunisia. Marine Micropaleontology, v. 13, p. 153-191. Carvalho, M.A., J.G.M. Filho and T.R. Menezes 2006. Palynofacies and sequence stratigraphy of the Aptian–Albian of the Sergipe Basin, Brazil. Sedimentary Geology, v. 192, p. 57-74. Darvishzad, B., E. Ghasemi-Nejad, S. Ghourchaei and G. Keller 2007. Planktonic foraminiferal biostratigraphy and faunal turnover across the Cretaceous-Tertiary boundary in southwestern Iran. Journal of Sciences, Islamic Republic of Iran, v. 18, no. 2, p. 139-149. Davey, R.J. and J. Rogers 1975. Palynomorph distribution in recent offshore sediments along two traverses off south west Africa. Marine Geology, v. 18, p. 213-225. Edet, J.J. and E.E. Nyong 1993. Depositional environments, sea-level history and palaeobiogeography of the late Campanian-Maastrichtian on the Calabar Flank, SE Nigeria. Palaeogeography, Palaeoclimatology, Palaeoecology, v. 102, p. 161-175. Eshet, Y., H.L. Cousminer and D. Habib 1988a. A model for using reworked palynomorphs as sedimentological and environmental indicators. Abstract, Palynology, v. 12, p. 236-237. Eshet, Y., H.L. Druckman, Y. Cousminer, D. Habib and W.S. Drugg 1988b. Reworked palynomorphs and their use in the determination of sedimentary cycles. Geology, v. 16, p. 662-665. Eshet, Y., S. Moshkovitz, D. Habib, C. Benjamin and M. Mogaritz 1992. Calcareous nannofossil and dinoflagellate stratigraphy across the Cretaceous/Tertiary boundary at Hor Hahar, Israel. Marine Micropaleontology, v. 18, p. 199-228.

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 Beiranvand et al.

Eshet, Y., A. Almogi-Labin and A. Bein 1994. Dinoflagellate cysts, paleoproductivity and upwelling systems: A Late Cretaceous example from Israel. Marine Micropaleontology, v. 23, p. 231-240. Espitalie, J., G. Deroo and F. Marquis 1986. La pyrolyse Rock-Eval et ses applications. Partie 3. Revue de l’Institut Français du Pétrole, Rueil-Malmaison, v. 41, no. 1, p. 73-89. Flügel, E. 2004. Microfacies of Carbonate Rocks: Analysis, Interpretation and Application. Springer-Verlag, Berlin, Heidelberg, New York, 976 p. Gebhardt, A.C., B. Gaye-Haake, D. Unger, N. Lahajnar and V. Ittekkot 2004. Recent particulate organic carbon and total suspended matter fluxes from the Ob and Yenisei Rivers into the Kara Sea (Siberia). Marine Geology, v. 207, p. 225-245. Ghasemi-Nejad, E., M.H. Hobbi and P. Schiøler 2006. Dinoflagellate and foraminiferal biostratigraphy of the Gurpi Formation (upper Santonian–upper Maastrichtian), Zagros Mountains, Iran. Cretaceous Research, v. 27, no. 6, p. 828- 835. Gradstein, F.M., J.G. Ogg, A.G. Smith, F.P. Agterberg, W. Bleeker, R.A. Cooper, V. Davydov, P. Gibbard, L. Hinnov, M.R. House (†), L. Lourens, H-P. Luterbacher, J. McArthur, M.J. Melchin, L.J. Robb, J. Shergold, M. Villeneuve, B.R. Wardlaw, J. Ali, H. Brinkhuis, F.J. Hilgen, J. Hooker, R.J. Howarth, A.H. Knoll, J. Laskar, S. Monechi, J. Powell, K.A. Plumb, I. Raffi, U. Röhl, P. Sadler, A. Sanfilippo, B. Schmitz, N.J. Shackleton, G.A. Shields, H. Strauss, J. Van Dam, J. Veizer, Th. van Kolfschoten and D. Wilson 2004. A Geologic Time Scale 2004. Cambridge University Press, 589 p. Gregory, W. and G.F. Hart 1992. Towards a predictive model for the palynologic response to sea-level changes. Palaios, v. 7, no. 1, p. 3-33. Habib, D and J.A. Miller 1989. Dinoflagellate species and organic facies evidence of marine transgression and regression in the Atlantic coastal plain. Palaeogeography, Palaeoclimatology, Palaeoecology, v. 74, p. 23-47. Habib, D., S. Moshkovitz and C. Kramer 1992. Dinoflagellate and calcareous nannofossil response to sea level change in Cretaceous-Tertiary boundary sections. Geology, v. 20, p. 165-168. Haq, B.U., J. Hardenbol and P.R. Vail 1988. Mesozoic and Cenozoic chronostratigraphy and cycles of sea-level change. In C.K. Wilgus, B.S. Hastings, C.G.St.C. Kendall, H.W. Posamentier, C.A. Ross and J.C. Van Wagoner (Eds.), Sea-level Changes: An Integrated Approach. Society of Economic Paleontologists and Mineralogists Special Publication no. 42, p. 71-108. Hardenbol, J., J. Thierry, M.B. Farley, T. Jacquin, P.-C. de Graciansky and P.R. Vail 1998. Mesozoic and Cenozoic sequence chronostratigraphic framework of European basins. In P.C. Graciansky, J. Hardenbol, T. Jacquin and P.R. Vail (Eds.), Mesozoic and Cenozoic Sequence Stratigraphy of European Basins. Society of Economic Paleontologists and Mineralogists Special Publication no. 60, p. 3-13, charts 1-8. Harland, R. 1973. Dinoflagellate cysts and acritarchs from the Bearpaw Formation (upper Campanian) of southern Alberta, Canada. Palaeontology, v. 16, p. 665-706. Helenes, J., C. de-Guerra and J. Vásquez 1998. Palynology and chronostratigraphy of the upper Cretaceous in the subsurface of the Barinas area, western Venezuela. The American Association of Petroleum Geologists Bulletin, v. 82, p. 1308-1328. Huber, B.T., D.A. Hodell and C.P. Hamilton 1995. Middle– late Cretaceous climate of the southern high latitudes: Stable isotopic evidence for minimal equator-to-pole thermal gradients. Geological Society of America Bulletin, v. 107, p. 1164- 1191. Huber, B.T., R.M. Leckie, R.D. Norris, T.J. Bralower and E. CoBabe 1999. Foraminiferal assemblage and stable isotopic change across the Cenomanian–Turonian boundary in the subtropical North Atlantic. Journal of Foraminiferal Research, v. 29, p. 392-417. James, G.A. and J.G. Wynd 1965. Stratigraphic nomenclature of the Iranian oil consortium agreement area. American Association of Petroleum Geologists Bulletin, v. 49, p. 2182-2245. Keller, G., L. Li and N. MacLeod 1995. The Cretaceous/Tertiary boundary stratotype section at El Kef, Tunisia: How catastrophic was the mass extinction? Palaeogeography, Palaeoclimatology, Palaeoecology, v. 119, p. 221-254. Kolodny, Y., 1980. Carbon isotopes and depositional environment of a high productivity sedimentary sequence - the case of the Mishash–Ghareb formations, Israel. Israel Journal of Earth-Sciences, v. 29, p. 147-156. Koop, W.J. and R. Stoneley 1982. Subsidence history of the Middle East Zagros Basin, Permian to Recent. In P. Kent, M.H.P. Bott, D.P. McKenzie and C.A. Williams (Eds.), The Evolution of Sedimentary Basins. Philosophical Transactions of the Royal Society of London, Part A, v. 305, p. 149-168. Koutsoukos, E.A.M. and M.B. Hart 1990. Cretaceous foraminiferal morphogroup distribution patterns, palaeocommunities and trophic structures: A case study from the Sergipe Basin, Brazil. Transactions of the Royal Society of Edinburgh, Earth Sciences, v. 81, p. 221-246. Leary, P.N., G.A. Carson, M.K.E. Cooper, M.B. Hart, D. Horne, I. Jarvis, A. Rosenfeld and B.A. Tocher 1989. The biotic response to the late Cenomanian oceanic anoxic event; integrated evidence from Dover, SE England. Journal of the Geological Society of London, v. 146, p. 311-317. Leine, L. 1986. Geology of the Tarfaya oil shale deposits, Morocco. Geologie en Mijnbouw, v. 65, p. 57-74. Leithold, E.L. and J. Bourgeois 1989. Sedimentation, sea-level change, and tectonics on an early Pleistocene continental shelf, northern California. Geological Society of America Bulletin, v. 101, p. 1209-1224. Li, L. and G. Keller 1998. Maastrichtian climate, productivity and faunal turnovers in planktic foraminifera in South Atlantic DSDP Sites 525 and 21. Marine Micropaleontology, v. 33, p. 55-86. Li, L., G. Keller and W. Stinnesbeck 1999. The Late Campanian and Maastrichtian in northwestern Tunisia: Paleoenvironmental inferences from lithology, macrofauna and benthic foraminifera. Cretaceous Research, v. 20, p. 231- 252. Moshkovitz, S. and D. Habib 1993. Calcareous nannofossil and dinoflagellate stratigraphy of the Cretaceous–Tertiary boundary, Alabama and Georgia. Micropaleontology, v. 39, p. 167-191. Motiei, H. 2003. Stratigraphy of Zagros, Treatise on the geology of Iran. Tehran, Iran, Geologicial Survey Press, 583 p. Nederbragt, A.J., R.N. Erlich, B.W. Fouke and G.M. Ganssen 1998. Palaeoecology of the biserial planktonic foraminifer Heterohelix moremani (Cushman) in the late Albian to middle Turonian Circum-North Atlantic. Palaeogeography, Palaeoclimatology, Palaeoecology, v. 144, p. 115-133. Parrish, J.T. and R.L. Curtish 1982. Atmospheric circulation, upwelling and organic-rich rocks in the Mesozoic and Cenozoic eras. Palaeogeography, Palaeoclimatology Palaeoecology, v. 40, p. 31-66.

100

Posamentier, H.W and P.R. Vail 1988. Eustatic controls on clastic deposition II - Sequence and systems tract models. In C.K. Wilgus, B.S. Hastings, C.G.St.C. Kendall, H.W. Posamentier, C.A. Ross and J.C. Van Wagoner (Eds.), Sea-level Changes: An Integrated Approach. Society of Economic Paleontologists and Mineralogists, Special Publication no. 42, p. 125-154. Premoli Silva, I. and W.V. Sliter 1999. Cretaceous paleoceanography: Evidence from planktonic foraminiferal evolution. In E. Barrera and C. Johnson (Eds.), Evolution of the Cretaceous Ocean-Climate System: Boulder, Colorado. Geological Society of America Special Publication 332, p. 301-328. Roberts, H.H. and J.M. Coleman 1988. Lithofacies characteristics of shallow expanded and condensed sections of the Louisiana distal shelf and upper slope. Gulf Coast Association of Geological Societies Transactions, v. 38, p. 291-301. Sampo, M. 1969. Microfacies and microfossils of the Zagros Area, southwestern Iran (from pre-Permian to Miocene). International Sedimentary Petrographical Series, v. 12, E. J. Brill, Leiden, 74 p., 105 pls. Serra, O. 1985. Sedimentary Environments from Wireline Logs, Schlumberger, Houston, 211 p. Setudehnia, A. 1972. Stratigraphic Lexicon of Iran. Union International des Sciences Geologiques, vol. 3., ASIE, southwest Iran. Setudehnia, A. 1978. The Mesozoic sequence in southwest Iran and adjacent areas. Journal of Petroleum Geology, v. 1, no. 1, p. 3-42. Sharland, P.R., R. Archer, D.M. Casey, R.B. Davies, S.H. Hall, A.P. Heward, A.D. Horbury and M.D. Simmons 2001. Arabian plate sequence stratigraphy. GeoArabia Special Publication 2, Gulf PetroLink, Bahrain, 371p., with 3 charts. Snedden, J. W. and C. Liu 2011. Recommendations for a uniform chronostratigraphic designation system for Phanerozoic depositional sequences. American Asssociation of Petroleum Geologists Bulletin, v. 95, no. 7, p. 1095-1122. Stoneley, R. 1974. Evolution of the continental margins bounding a former southern Tethys. In C.A. Burk and C.L. Darke (Eds.), The Geology of Continental Margins. Springer-Verlag, New York, p. 889-903. Stoneley, R. 1990. The Middle East Basin: A Summary Overview. In J. Brooks (Ed.), Classic Petroleum Provinces. Geological Society of London, Special Publication no. 50, p. 293-298. Tyson, R.V. 1995. Sedimentary Organic Matter, Organic Facies and Palynofacies. Chapman and Hall, London, 615 p. van der Zwaan, G.J., F.J. Jorissen and H.C. de Stigter 1990. The depth dependency of planktonic/benthic foraminiferal ratios: Constraints and applications. Marine Geology, v. 95, p. 1-16. Wall, D., B. Dale, G.P. Lohmann and W.K. Smith 1977. The environmental and climatic distribution of dinoflagellate cysts in modem marine sediments from regions in the North and South Atlantic oceans and adjacent seas. Marine Micropaleontology, v. 2, p. 121-200. West, O.L.O., R.M. Leckie and M. Schmidt 1998. Foraminiferal paleoecology and paleoceanography of the Greenhorn Cycle along the southwestern margin of the Western Interior Sea. Society of Economic Paleontologists and Mineralogists Concepts Sedimentology and Paleontology, no. 6, p. 79-99. Williams, G.L. 1977. Dinoflagellate cysts, their classification, biostratigraphy and palaeoecology. In A.T.S. Ramsay (Ed.), Oceanic Micropalaeontology. Academic Press, London, p. 1231-1325. Wood, G.D., A.M. Gabriel and J.C. Lawson 1996. Palynological techniques-processing and microscopy. In J. Jansonius and D.C. McGregor (Eds.), Palynology: Principles and Applications. American Association of Stratigraphic Palynologists Foundation, Dallas, p. 29-50. Wynd, J.G. 1965. Biofacies of the Iranian Oil Consortium Agreement area. Iranian Oil Operating Companies, Geological and Exploration Division, Report 1082, 89 p. Ziegler, M.A. 2001. Late Permian to Holocene paleofacies evolution of the Arabian Plate and its hydrocarbon occurrences. GeoArabia, v. 6, no. 3, p. 445-504.

ABOUT THE AUTHORS

Bijan Beiranvand is a Basin Analyst with the Research Institute of Petroleum Industry (RIPI) at the National Iranian Oil Company (NIOC). He received a BSc (1992) in Geology (Hydrogeology) from Shiraz University, Iran and an MSc (1995) and PhD in Geology (Stratigraphy) from the University of Tehran, Iran. Bijan joined the Bureau of Engineering of the Iranian Air force as a Supervisor of the Soil Mechanic Laboratory from 1995 to 1997, then he worked as a Supervisor of the Exploration Department of the Bureau of Lorestan and Chahar Mahal-o Bakhtiari Provinces (1997–2000). Bijan joined RIPI in 2000 and now he is working for the Petroleum Geology Department of this company. His interests are stratigraphy and basin analysis, basin modeling and reservoir geology. [email protected]

101

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 Beiranvand et al.

Ebrahim Ghasemi-Nejad is currently a professor of Stratigraphy and Marine Palynology with the Department of Geology of the University of Tehran, Iran. He received a BSc degree in Geology from Ferdowsi University of Mashhad and an MSc degree in the field of Stratigraphy and Micropaleontology from the University of Tehran. He then continued his studies for a PhD degree at the Department of Earth Sciences of the University of Saskatchewan, Canada from 1990 to 1995. His PhD thesis is on the Palynology and Paleoenvironment of the Jurassic strata of the Jura Mountain in Switzerland. After finishing his PhD (1996) he returned to Iran to work as an academic with the University of Tehran. Since 1996 he has been working there, teaching graduate and undergraduate students of Geology and handling research projects in the fields of marine palynology (dinoflagellates), palynostratigraphy, foraminiferal biostratigraphy and evaluation of petroleum potential of different rock units via using palynology, palynofacies and organic geochemistry analyses. Ebrahim has presented and published many papers and is currently active as the chief editor of the Geopersia journal and at the same as editorial board member for several national geological journals. His publications include papers on Triassic, Jurassic and Cretaceous marine palynology, palynostratigraphy, paleoenvironment, and petroleum potential evaluation of different rock units. He was promoted to the level of professorship in 2009. [email protected]

Mohammad Reza Kamali is Associate Professor and Deputy, International Business Development, Center for Exploration and Production–Research Institute of Petroleum Industry, Iran. He received his BSc and MSc in Geology from the University of Mysore (India) in 1982 and 1984 respectively. After joining to National Iranian Oil Company in 1985 as a Senior Geologist, Mohammad continued his PhD in Petroleum Geology in Adelaide University, Australia and graduated in 1996. Mohammad has published six books on Petroleum Geochemistry and Geology is author and co-author for more than 95 research papers published in national and international journals. His research interests are petroleum geochemistry, basin modelling and reservoir geology. He is an active member of AAPG, SPE, EAGE, Geological Society of Iran and Petroleum Engineering Society of Iran. [email protected]

Akram Ahmadi worked as a Petroleum Geologist with the Research Institute of Petroleum Industry (RIPI) at the National Iranian Oil Company from 2004 to 2009. She received a BSc (2002) in Geology from Kharazmi University, Iran and an MSc (2005) in Petroleum Geology from the University of Tehran, Iran. Akram’s interests are reservoir evaluation and petrophysical data interpretation. [email protected]

Manuscript received October 1, 2012

Revised July 25, 2013 Accepted October 29, 2013

102

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 GeoArabia Beiranvand, B., E. Ghasemi-Nejad, M.R. Kamali and A. Ahmadi 2014. Sequence stratigraphy of the Late Cretaceous–Paleocene Gurpi Formation in southwest Iran. GeoArabia, v. 19, no. 2, p. 89-102, 4 enclosures Enclosure 1: Danial Section, Kuh-e Pyun, Iran

Biostratigraphy Lithostratigraphy Gamma-Ray Log Chemostratigraphy Microplanktonic Characteristics Sea-Level Cycles and Sequences

Minor Palynological Palynological Species Global Coastal Onlaps Sea-Level Lithology Major Element Planktonic Element Elements Indices Diversity This Study (Snedden and Liu, 2011) History Marly Limestone Laboratory U/Th Foraminifera Organic NGR 0.15 0.55 Enviromental Geo- Shaly Marl Rb P Sporomorph Index Planktonic 20 (cps) 25 K/Al Group chemistry Marl Ti Cs GTS Foraminifera Biozonation Foraminiferal Sedimentation Rate 0.2 1.0 Hedbergelids 2012 (Beiranvand and Ghasemi-Nejad, 2013) Uranium Ba Sn Lining Si/Al 0 (cm/ky) 10 Age (Ma) 5 (ppm) 12 Dinocyst Thickness (m) Texture K Mo 1 9 Heterohelicids Group System This Study Member Hydrocarbon Coastal Onlaps Formation Thorium SiO2 Mg/Al Tracts Mg Zn Foraminifera Index +200 (m) 0 Polarity Chron Phytoclast Sample Number Non-Keeld Paleomagnetic Lithology Field 20 (ppm) 40 0 0.8 PMI 0 300 Group Group System Tracts Record Gamma-Ray Al Zr Global Mean (API) Mg/Ca 0 2,000 Cr30n Normal Potassium CaO Amorph Dinocyst Keeld Total Organic Sea-Level Changes Fe Ir Sequence Terminology Sequence Terminology 0 0.2 Group S/C Group Carbon Landward Basinward Snedden and Liu, 2011 Cr30r Reverse 0 100 0.5 (%) 1.2 0 (%) 100

Mudstone Wackestone Packstone Grainstone 0 (%) 30 0 (%) 1 0 (%) 100 0 8 0 (%) 100 0 (%) 100 0 (%) 1 0.5 0.4 0.3 0.2 0.1 0 +200 (m) 0 59.2 TST 59.2 Ma 340 P4 P4a G. pseudobulloides–P. var. interval subzone SB 59.3 P3b I. albeari lowest-occurrence subzone 145 LHST HST C26r P3 I. pusilla 140 EHST Gu9 Sel1 Selandian P3a MFS 320 partial-range subzone MFS 60.6 61.6 TST TST C27n P2 P. uncinata lowest-occurrence zone 135 LST SB LHST 61.1 300 EHST MFS PALEOCENE Gu8 Da4 C27r G. compress–P. uncinata 130 MFS P1c TST LST 61.7 lowest-occurrence subzone TST SB Danian C28n P1 62.1 125 LHST MFS Da3 C28r Subbotina triloculinoides P1b Gu7 63.4 lowest-occurrence subzone 120 EHST 64.5 Da2 280 C29n MFS Pα P1a P. pseudobulloides–S. triloculinoides interval subzone TST LST Cretaceous/ TST Da1 66.0 CF1 Plummerita hantkeninoides 115 SB K/Pg P0 Paleogene LHST 65.1 66.0 Ma Pseudoguembelina palpebra Boundary HST 110 EHST C29r CF2 partial-range subzone MFS MFS 66.3 260 105 Gu6 Ma5 TST 100 TST CF3 Pseudoguembelina hariaensis partial-range subzone 95 LST C30n SB 240 LHST 67.7

R. fructicosa 90 HST C30r CF4 partial-range subzone EHST Gu5 Abathomphalus mayaroensis Ma4 MFS 85 MFS TST 68.0 220 CF5 Pseudotextularia intermedia partial-range subzone TST LST C31n SB 69.0 80 HST MFS Contusotruncana contusa CF6 LHST 69.2 Ma3 interval zone 200 75 TST

69.5 70

180 Maastrichtian

EHST HST Gu4 Gurpi 65

160 MFS Ma2 MFS 69.7

Emam Hassan 60 C31r

TST 140 55 TST Gansserina gansseri CF7 interval zone

50 LST SB LATE CRETACEOUS LATE 70.1 120 LHST 45 HST EHST Gu3 MFS MFS 70.5 Ma1 72.1 40 TST TST 100 LST SB 70.8 35 lowstand LST systems tract LHST Transgressive TST systems tract 80 30 LHST Late highstand systems tract HST EHST Early highstand C32n systems tract Cam9 Globotruncana aegyptiaca 25 EHST Gu2 CF8 interval zone 60 MFS MFS 20 72.1 TST TST Campanian LST 40 15 SB 73.1 LHST Globotruncanella havanensis CF9 partial-range zone 10 HST 20 C33n Gu1 Cam8 EHST Maximum MFS Flooding 05 Surface MFS Radotruncana MFS CF10 SB Sequence 74.5 calcarata taxon-range zone TST Boundary 0 SB Enclosure 1: Correlation chart of different data to recognize sea-level changes and identify systems tracts and correspondence to the global sea-level cycles and systems tracts. See table 1 for abbreviations.

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021 GeoArabia Beiranvand, B., E. Ghasemi-Nejad, M.R. Kamali and A. Ahmadi 2014. Enclosure 2a: Danial Reference Section, Enclosure 2b: Sequence stratigraphy of the Late Cretaceous–Paleocene Gurpi Formation in southwest Iran. GeoArabia, v. 19, no. 2, p. 89-102, 4 enclosures Pyun Anticline, Northeast Izeh (Iran) Gurpi Type Section, Gurpi Anticline, Northeast Lali (Iran)

Sequence Biostratigraphy Lithostratigraphy Sea-Level Cycles and Sequences Stratigraphy Lithostratigraphy Biostratigraphy

Lithology Global Coastal Onlaps Lithology This Study (Snedden and Liu, 2011) This Study Marly Limestone Marly Limestone

Shaly Marl Shaly Marl Marl Marl GTS GTS Planktonic Planktonic 2012 Dinocyst Events 2012 Foraminifera Biozonation Foraminifera Biozonation Age (Ma) (This Study) Age (Ma) (Beiranvand and Ghasemi-Nejad, 2013) (This Study) Thickness (m) Texture System System Texture Thickness (m)

Member Coastal Onlaps Member

Formation Tracts Tracts Formation Polarity Chron Paleomagnetic Sample Number Lithology Field Field Lithology Sample Number Record Gamma-Ray System Tracts Gamma-Ray Cr30n Normal (API) (API) Sequence Terminology Landward Basinward Sequence Terminology Sequence Terminology Cr30r Reverse 0 100 0 100

Mudstone Wackestone Packstone Grainstone 0.5 0.4 0.3 0.2 0.1 0 Mudstone Wackestone Packstone Grainstone 59.2 TST 59.2 340 340 P4 P4a G. pseudobulloides–P. var. ISZ SB 59.3 SB G. pseudobulloides–P. var. P4a P4 I. albeari C. cornuta LHST I. albeari P3b 145 LHST P3b lowest-occurrence subzone HST lowest-occurrence subzone C26r P3 D. colligrum/ 130 P3 EHST Gu9 I. pusilla C. gracile 140 EHST Gu9 I. pusilla

Selandian MFS Sel1 Selandian 320 P3a partial-range subzone MFS MFS partial-range subzone P3a 320 O. pulcherrimumum 60.6 TST 125 61.6 TST 61.6 TST LST SB P. uncinata C27n P2 P. uncinata C. diebelli 135 LST LHST P2 lowest-occurrence zone SB 120 lowest-occurrence zone LHST 61.1 Gu8 300 EHST MFS 300 PALEOCENE C27r Gu8 Da4 G. compress–P. uncinata H. bulbosum 130 MFS MFS G. compress–P. uncinata P1c TST LST 61.7 P1c PALEOCENE lowest-occurrence subzone TST SB lowest-occurrence subzone SB LHST 115 Danian C28n P1 62.1

C. inodes D. acuminatum 125 LHST MFS Da3 Danian C28r Subbotina triloculinoides P1b D. califonicum Gu7 63.4 EHST P1 lowest-occurrence subzone S. inornata 120 EHST 64.5 Da2 280 C29n MFS Gu7 MFS S. triloculinoides 280 Pα P1a P. pseudobulloides–S. triloculinoides ISZ C. cornuta TST LST P1b Cretaceous/ TST Da1 110 lowest-occurrence subzone 66.0 S. inornata TST CF1 Plummerita hantkeninoides 115 SB P0 A. reticulata A. euclaense Paleogene LHST 65.1 LST SB K/Pg 66.0 Pseudoguembelina palpebra C.utinensis Boundary HST D. spp 110 EHST LHST 66.0 Ma C29r CF2 partial-range subzone MFS 105 MFS 66.3 Pseudoguembelina palpebra 260 105 CF2 260 H. bulbosum Gu6 Ma5 EHST partial-range subzone M. seelandica TST 100 TST MFS 100 CF3 Pseudoguembelina hariaensis partial-range subzone 95 LST Gu6 95 C30n SB 240 LHST 67.7 240 TST Pseudoguembelina hariaensis R. fructicosa 90 HST 90 partial-range subzone CF3 C30r CF4 partial-range subzone EHST Gu5 Abathomphalus mayaroensis Ma4 MFS 85 MFS 68.0 85 CF5 Pseudotextularia intermedia PRSZ TST 220 TST LST SB 220 C31n P. magniﬁcum LST SB 69.0 LHST M. druggii 80 80 HST Abathomphalus mayaroensis MFS R. fructicosa Contusotruncana contusa CF4 CF6 LHST 69.2 Ma3 partial-range subzone interval zone EHST 75 200 75 200 TST Gu5 MFS 70 Pseudotextularia intermedia 69.5 CF5 70 partial-range subzone

180 Maastrichtian 180 TST 65 EHST HST Gu4 Gurpi Gurpi

65 60 Maastrichtian LST Contusotruncana contusa CF6 160 SB interval zone 160 MFS Ma2 MFS 69.7

Emam Hassan 60 LHST 50 C31r P. infusoroides P. cingulatum A. senonensis TST 140 55 TST 140 Gansserina gansseri 45 CF7 interval zone

50 LST Gu4 Emam Hassan SB EHST LATE CRETACEOUS LATE 70.1

120 LHST CRETACEOUS LATE 120 HST 45 40 EHST Gu3 MFS MFS 70.5 Ma1 72.1 40 TST MFS TST 35 100 TST 100 LST Gansserina gansseri CF7 SB interval zone 70.8 LST SB 35 lowstand O. operculata LST systems tract LHST LHST Transgressive 30 TST systems tract EHST T. suspectum and O. porifera 80 30 Late highstand Gu3 MFS 72.1 80 LHST HST systems tract TST X. ceratioides EHST Early highstand 25 C32n systems tract Cam9 Globotruncana aegyptiaca EHST Gu2 LST CF8 25 SB interva zone LHST 60 20 60 MFS C. verrucosa MFS 20 72.1 EHST Gu2 TST TST MFS Globotruncana aegyptiaca

Campanian 15 LST CF8 A. coronata TST interval zone 40 15 SB 40 73.1

LHST LST Campanian Globotruncanella havanensis SB 10 CF9 P. cingulatum CF9 partial-range zone 10 EHST Globotruncanella havanensis partial-range zone M. Druggii HST 20 C33n Gu1 Cam8 EHST 20 P. pyrophorum EHST Maximum Gu1 MFS Flooding MFS 05 Surface MFS 5 Radotruncana calcarata Radotruncana TST CF10 CF10 MFS SB Sequence 74.5 taxon-range zone calcarata taxon-range zone TST Boundary 0 SB 0 Enclosure 2: Sequence stratigraphic correlation of the Gurpi depositional sequences between the Danial and Gurpi sections. See table 1 for abbreviations.

Lithology Nannofossil Event Marly Limestone

Shaly Marl Marl

GTS Planktonic 2012 Lithology Texture Foraminifera Biozonation Age (Ma) (Beiranvand and Ghasemi-Nejad, 2013) Member Formation Thickness (m) Polarity Chron Sample Number

sp.

Paleomagnetic Record Cr30n Normal

Cr30r Reverse Mudstone Wackestone Packstone Grainstone

Archaeoglobigerina cretacea

Archaeoglobigerina blowi

Abathomphalus mayaroensis

Acarinina strabocella

Acarinina subsphaerica Acarinina mckannai

Radotruncana calcarata Globotruncana lapparenti Globotruncana bulloides Contusotruncana fornicata Contusotruncana plummerae Globotruncana ventricosa Globotruncana linneiana Globigerinelloides alvarezi Radotruncana subspinosa Globotruncanita stuartiformis Globotruncana arca Globotruncana mariei Globotruncana falsostuarti Globotruncana Hedbergella monmouthensis Hedbergella holmdelensis Heterohelix globulosa Heterohelix labellosa Heterohelix striata Heterohelix punctulata Laeviheterohelix dentata Laeviheterohelix pulchra Rugoglobigerina rugosa

Globotruncanella havanensis Globotruncanita stuarti Planoglobulina carseyae Planoglobulina brazoensis

Pseudoguembelina costulata Globotruncana aegyptiaca Globotruncanella petaloidea Rugoglobigerina macrocephala Gansserina gansseri Globotruncana dupeublei Rugoglobigerina milamensis Plummerita reicheli Pseudoguembelina excolata Rugoglobigerina hexacamerata Pseudotextularia elegans Racemiguembelina powelli

Contusotruncana contusa Pseudotextularia intermedia

Racemiguembelina fructicosa Trinitella scotti Pseudoguembelina hariaensis Pseudoguembelina palpebra Guembelitria irregularis Guembelitria cretacea Plummerita hantkeninoides Woodringina hornerstownensis Eoglobigerina fringa Chiloguembelina midwayensis Parvularugoglobigerina eugubina Guembelitria trifolia Subbotina trivialis Globotruncanita conica Eoglobigerina eobulloides Praemurica pseudoinconstans Praemurica taurica Globoconusa daubjergensis Parasubbotina pseudobulloides Globanomalina planocompressa Praemurica inconstans Subbotina triloculinoides Globanomalina compressa Parasubbotina varianta Globanomalina imitata Eoglobigerin spiralis Praemurica uncinata Praemurica praecursoria Morozovella praeangulata Subbotina cancellata

Morozovella angulata Igorina pusilla Parasubbotina variospira

Subbotina triangularis 59.2 340 P4 P4a G. pseudobulloides–P. var. ISZ P3b I. albeari 145 lowest-occurrence subzone C26r P3 140 I. pusilla Selandian 320 P3a partial-range subzone 61.6 C27n 135 P2 P. uncinata lowest-occurrence zone P α TRZ 300 PALEOCENE C27r 130 G. compress–P. uncinata P1c lowest-occurrence subzone

Danian C28n P1 125 C28r Subbotina triloculinoides

P. eugubina P. P1b 120 lowest-occurrence subzone 280 C29n P1a P. pseudobulloides–S. triloculinoides ISZ Cretaceous/Paleogene 66.0 K/Pg 115 66.0 Ma CF1 Plummerita hantkeninoides Boundary

110 Pseudoguembelina palpebra C29r CF2 partial-range subzone 260 105

100 CF3 Pseudoguembelina hariaensis partial-range subzone C30n 95

240 G. cretacea P0 PRZ

90 R. fructicosa C30r CF4 partial-range subzone Abathomphalus mayaroensis

85 Pseudotextularia intermedia PRSZ 220 CF5 C31n 80

CF6 Contusotruncana contusa interval zone 200 75

180 Maastrichtian Gurpi 65

160

Emam Hassan 60 C31r

140 55 Gansserina gansseri CF7 interval zone

50 LATE CRETACEOUS LATE 120 45

72.1 40 100

80 30

C32n 25 Globotruncana aegyptiaca CF8 interval zone 60

20 Campanian

40 15

Globotruncanella havanensis CF9 10 partial-range zone

20 C33n

05 Radotruncana CF10 calcarata taxon-range zone 0 Enclosure 3: Full range chart of the Gurpi Formation in the Danial Section. See table 1 for abbreviations.

Lithology Nannofossil Event Marly Limestone

Shaly Marl Marl

GTS Planktonic 2012 Lithology Texture Foraminifera Biozonation Age (Ma) (This Study) Member Formation Thickness (m) Polarity Chron Sample Number

Paleomagnetic Record Cr30n Normal

Cr30r Reverse Mudstone Wackestone Packstone Grainstone

Archaeoglobigerina cretacea

Archaeoglobigerina blowi

Abathomphalus mayaroensis

Acarinina strabocella

Acarinina subsphaerica Acarinina mckannai

Radotruncana calcarata Globotruncana lapparenti Globotruncana bulloides Contusotruncana fornicata Globotruncana ventricosa Globotruncana linneiana Globigerinelloides alvarezi Radotruncana subspinosa Globotruncanita stuartiformis Globotruncana arca Globotruncana mariei Globotruncana falsostuarti Hedbergella monmouthensis Hedbergella holmdelensis Heterohelix globulosa Heterohelix labellosa Heterohelix striata Heterohelix punctulata Laeviheterohelix dentata Laeviheterohelix pulchra Rugoglobigerina rugosa

Globotruncanella havanensis Globotruncanita stuarti

Globotruncana aegyptiaca Globotruncanella petaloidea Rugoglobigerina macrocephala Gansserina gansseri Globotruncana dupeublei Rugoglobigerina milamensis Plummerita reicheli Rugoglobigerina hexacamerata Pseudotextularia elegans Racemiguembelina powelli Globotruncanita conica Contusotruncana contusa Pseudotextularia intermedia

Racemiguembelina fructicosa Trinitella scotti Pseudoguembelina hariaensis Pseudoguembelina palpebra Eoglobigerina eobulloides Praemurica pseudoinconstans Praemurica taurica Globoconusa daubjergensis Parasubbotina pseudobulloides Globanomalina planocompressa Praemurica inconstans Subbotina triloculinoides Globanomalina compressa Parasubbotina varianta Globanomalina imitata Eoglobigerina spiralis Praemurica uncinata Praemurica praecursoria Morozovella praeangulata Subbotina cancellata Subbotina triangularis Morozovella angulata Igorina pusilla Parasubbotina variospira

Morozovella conicotruncana Igorina albeari Morozovella apanthesma Subbotina velascoensis Globanomalina chapmani Morozovella velascoensis Morozovella acutispira Globanomalina pseudomenardii Morozovella occlusa 59.2 340 P4 P4a G. pseudobulloides–P. var. P3b I. albeari lowest-occurrence subzone C26r 130 P3 I. pusilla Selandian 320 P3a partial-range subzone 125 61.6 C27n P. uncinata 120 P2 lowest-occurrence zone

300 C27r G. compress–P. uncinata PALEOCENE P1c lowest-occurrence subzone 115 C28n Danian C28r P1 C29n S. triloculinoides 280 P1b 110 lowest-occurrence subzone Cretaceous/Paleogene 66.0 K/Pg 66.0 Ma Boundary C29r 105 Pseudoguembelina palpebra 260 CF2 partial-range subzone

100

C30n 95 240 Pseudoguembelina hariaensis 90 CF3 partial-range subzone C30r

85 220 C31n 80 Abathomphalus mayaroensis R. fructicosa CF4 75 partial-range subzone 200

70 CF5 Pseudotextularia intermedia partial-range subzone 180 65 Gurpi

Maastrichtian 60 CF6 Contusotruncana contusa 160 interval zone

50 C31r

140 45 Emam Hassan

120 CRETACEOUS LATE 40

35 100 CF7 Gansserina gansseri interval zone

80 72.1

C32n 25

60 20

15 Globotruncana aegyptiaca CF8 interval zone 40

Campanian 10 CF9 Globotruncanella havanensis partial-range zone 20 C33n

5 Radotruncana calcarata CF10 taxon-range zone 0 Enclosure 4: Full range chart of the Gurpi Formation in the Gurpi Section. See table 1 for abbreviations.

Downloaded from http://pubs.geoscienceworld.org/geoarabia/article-pdf/19/2/89/4567952/beiranvand.pdf by guest on 02 October 2021