Saftic.indd Tertiary Subsurface Facies, Source Rocks and Hydrocarbon Reservoirs in the SW Part of the Pannonian Basin (Northern Croatia and South-Western Hungary) Bruno SAFTIĆ1, Josipa VELIĆ1, Orsolya SZTANÓ2, Györgyi JUHÁSZ3 and Željko IVKOVIĆ4 1. INTRODUCTION Hydrocarbon exploration resulted in an enormous amount of data about the structure, sedimentary fill and evolution of the Pannonian Basin (HORVÁTH & ROYDEN, 1981; ROYDEN & HORVÁTH, 1988; BÉRCZI et al., 1988; SZENTGYÖRGYI & JUHÁSZ, 1988; JUHÁSZ, 1991, 1994; HORVÁTH & TARI, 1999; PRELOGOVIĆ et al., 1998; LUČIĆ et al., 2001; VELIĆ et al., 2002). Some of the most important results from the Croatian part of the study area regretfully remained unpublished (ŠIMON, 1980); they have been used here together with other data from the archives of INA, Zagreb (Croatian national oil company). Above the historically complex basement the struc- ture of the study area is complicated (CSONTOS et al., 1992; CSONTOS 1995; CSONTOS & NAGYMARO- SY, 1998): elongated basement highs and narrow depressions developed during Mid-Miocene rifting (HORVÁTH & ROYDEN, 1981; HORVÁTH, 1993), refigured by several phases of basin inversion prior to, during and after the subsequent thermal subsidence (HORVÁTH, 1995; HORVÁTH & CLOETINGH, 1996; PRELOGOVIĆ et al., 1998; TARI & PAMIĆ, 1998; HORVÁTH & TARI, 1999; FODOR et al., 1999; PAVELIĆ, 2001; TOMLJENOVIĆ & CSONTOS, 2001; CSONTOS et al., 2002). A Neogene sedimentary succession up to 5–7 km in thickness was accumulated in these depressions which are separated by uplifted and partially eroded units. The aforementioned depressions (sub basins) – the most northern Mura, the western and the eastern part of the Drava, the Sava and the Slavon- ija–Srijem depressions (Fig. 1), are hosts of important source rocks and petroleum reservoir rocks. Some of the basement highs are exposed as Palaeozoic–Mesozo- ic Mountains (e.g. Mecsek, Villány, Papuk, Psunj, etc.), others are still covered by about a 1–2 km thick pile of Neogene sediments. Conditions for the accumulation and preservation of large quantities of marine and terrigeneous organic matter were most favourable during the Badenian, Sarmatian and Early Pannonian in deep basin settings (HORVÁTH et al., 1988; DANK, 1988; ALAJBEG et al., 1990; BARIĆ et al., 1989, 1992; HERNITZ et al., 1995a, b). The generation of hydrocarbons was Geologia Croatica 56/1 101–122 13 Figs. ZAGREB 2003 Key words: Paratethys, Lake Pannon, Miocene, Pliocene, Quaternary, Subsurface Facies, Source Rocks, Hydrocarbon Reservoirs, Sedimentation, Tectonics, Pannonian Basin, Croatia, Hungary. 1 Faculty of Mining, Geology and Petroleum Engineering, University of Zagreb, Pierottijeva 6, 10000 Zagreb, Croatia; e-mail: bsaft@rgn.hr 2 Department of Geology, Eötvös University, Pázmány Péter st. 1/c, 1117 Budapest, Hungary; e-mail: sztano@ludens.elte.hu 3 Hungarian Geological Institute, Stefánia út 14, 1143 Budapest, Hungary; e-mail: juhaszgy@mafi.hu 4 INA industrija nafte dd., Naftaplin, Šubićeva 29, 10000 Zagreb, Croatia; e-mail: zeljko.ivkovic@ina.hr Abstract The Neogene sedimentary successions of the Drava, Sava and Slavon- ija–Srijem depressions in the SW part of the Pannonian Basin System are built up of three 2nd order megacycles separated by four major erosional unconformities. The first megacycle contains terrestrial to marine syn-rift and early post-rift sediments of Early to Mid-Miocene age. The second is built up of Late Miocene Lake Pannon deposits, while the third contains those sediments which were deposited in the remnants of Lake Pannon and in the subsequent fluvial systems, in areas of continuous subsidence associated with basin inversion from the Pliocene onwards. Most of the petroleum source rocks and reservoir rocks are of Miocene age and were formed during the first and second depositional megacycle. Conditions for the accumula- tion and preservation of large quantities of marine and terrigenous organic matter were most favourable during the Badenian, Sarmatian and Early Pannonian, in deep basin settings, partly associated with rifting. The generation of hydrocarbons was promoted by relatively high geothermal gradients during the initial and subsequent thermal subsidence. Various sedimentary environments produced deposits with good reservoir characteristics: e.g. fault-related talus breccia (mainly Lower Miocene), reefs (mainly Badenian), coastal, shallow marine (Karpatian, Badenian) and deltaic (Pannonian–Pontian) sand bodies or turbiditic sand lobes (mainly Pannonian). The hydrocarbon (HC) migration paths were often provided by the major unconformi- ties bounding the three megacycles, as well as by faults, particularly around the basement highs. 102 Geologia Croatica 56/1 Pannonian Basin is a historically complex basin devel- oped parallel with the Alpine–Carpathian orogen. Its northern part, the Alcapa unit, escaped from the Alpine zone and was juxtaposed to the southern Tisza unit along the Mid-Hungarian Line during the Late Eocene– Early Miocene (CSONTOS & NAGYMAROSY, 1998, and references therein). Thereafter the Neogene Panno- nian Basin proper, opened as a back arc basin in which extension was manifested through a set of low angle normal (TARI et al., 1999) and strike-slip faults. Exten- sion was coeval with the compression within the Car- pathian thrust belt (HORVÁTH & ROYDEN, 1981). Paratethys was an enclosed sea existing from Oli- gocene to Middle Miocene times consisting of a chain of basins of various tectonic origin (BÁLDI, 1980). These basins were covered most of the time by the same mass of water, sharing a similar aquatic biota, communicating between the Mediterranean and the Indopacific realm. When Paratethys became separated from other seas in the beginning of the Late Miocene, the Lake Pannon individualized with decreasing salinity (cf. MAGYAR & GEARY, 1999; MÜLLER et al., 1999). As both Paratethys and Lake Pannon were characterized by an increasing rate of endemism, a local chronostratigraphic scheme was developed as used in Figs. 2–4 (for detailed discussion the reader is referred to VASS et al., 1987; RÖGL & DAXNER-HÖCK, 1996; VASS, 1999; MAGYAR et al., 1999a, b). promoted by relatively high geothermal gradients (between 3–4 and 7ºC/100 m – HORVÁTH et al., 1988) during the subsequent thermal subsidence. Various sedimentary environments produced deposits with good reservoir characteristics (Fig. 2), e.g. fault-related talus breccia (mainly Lower Miocene), reefs (mainly Badenian), coastal, shallow and open marine (Karpatian, Badenian, Sarmatian) and deltaic (Pannonian–Pontian) sand bodies or turbiditic sand lobes (mainly Pannonian). These are only the most frequent occurrences; almost all of the mentioned environments existed contemporaneously. The HC migration paths were often provided by the major unconformities as well as by faults, particularly around the basement highs. The late stage inversion of the basin promoted accumulation of the hydrocarbons in areas of active uplift (HORVÁTH, 1995). The aim of this paper is to review the sedimen- tary facies and depositional environments with special regard to the potential source and reservoir rocks in the SE part of the Pannonian Basin. 2. GEOLOGICAL SETTING The studied sediments from SW Hungary and N Croatia were deposited in the Pannonian Basin, the largest intramontaine basin of Central Paratethys (Fig. 1). The Fig. 1 Location map. 103 S aftić, Velić, S ztanó, Juhász & Ivković: Tertiary S ubsurface Facies, S ource R ocks... Fig. 2 S chem atic section – lithostratigraphic units and the three N eogene depositional m egacycles. SAVA DEPRESSION S N SIROKO POlJE Fm. nDSTAR IVANIC Pm. IVANIC GRAD Fm. Legend: ~ ~-... -_ --.i L ~ ~ OSFm. ,, ~--- ~ vokanic~ ~ ~ ~ v v pyroclastic rocks ~ , I~~h~~ 1J a s e Neog D sandstone - arkoIe ~ Oo~ gravel .,oa[ D sandstone - graywacke II ... d ~-.- D=l I::~=I clay Effi3 calcareous ~ m.d Neogene depositional megacycles s '.6 lULoGo.RA, FlIl. c:: W Post-rift 0 ? 0 subsidence ::; Z KLOiiTARIVANICFm. ~ a: z w 0 Z a. z a. IVANIe GRAD Fm. 8: ::;) ~ ,~ ~ IBADENIAN KARPATIAN -J OTTNA. NGw<1 I _ _ 18.3 EGGI:.N- I BURGIAN EGERIAN - DRAVA DEPRESSION N VERAFm UJ?ALU'{InI.· • o.:ijc:k ...... u > o ~ Jr' A~" f,lPl6 '£I1Jl. ~. ,:;."~ " ~ > Correlation horizons: ~ major and minor ~ unconformities EJ stI"',.,r-"" conformable ""'",,, 104 Geologia Croatica 56/1 1993; VAKARCS et al., 1994; VAKARCS, 1997; SAC- CHI et al., 1999), whereas it happened in Croatia during the Pontian (7–5 Ma). The third phase of basin evolution is represented by repeated compressional events in relation to basin inversion from the Pliocene onwards (HORVÁTH, 1995; CSONTOS et al., 2002). All of the compres- sional events resulted in the reactivation of former faults mainly as reverse types, and gentle folding of the overlying strata (see Fig. 5 and Figs. 6–9). Between the areas of active uplift, subsidence continued mainly in the deepest depressions resulting in the formation of an occasionally thick Pliocene–Quaternary succession (HORVÁTH & CLOETINGH, 1996; PRELOGOVIĆ et al., 1998; HORVÁTH & TARI, 1999; LUČIĆ et al., 2001; VELIĆ et al., 2002). 3. SEDIMENTARY FACIES AND ENVIRONMENTS An overview of the major lithostratigraphic units in the study area, their composition and contacts is given in Fig. 2, with more details and description of depositional characteristics in Figs. 3 and 4. The Neogene to Quater- nary basin fill in Croatia starts with Egerian–Eggenbur- gian sediments documented in the Mura depression, in the Hrvatsko zagorje region and in the nearby western- most part of the Drava depression. This was referred to as the “Western marginal part of the Pannonian basin” (LUČIĆ et al., 2001) interpreted to be the place of formation of the first depocentres, concurrent with the initial tectonic transport in the formation of the Pannon- ian Basin System. In the basal sections coarse clastics, sandstones and clays with coal were observed, depos- ited in paralic environments. The marine tuffaceous sandstones (“Macelj sandstones”) are interpreted to be of Eggenburgian age (LUČIĆ et al., 2001). These sedi- ments were not drawn in Figs. 2–4 because they belong to the end of the Palaeogene basin cycle and not to the rifting of PBS. Anyway, they are present only in the NW corner of the study area and we wanted to empha- size the most important similarities between the Hun- garian and Croatian stratigraphic data. The most impor- tant of these is that the overlying Ottnangian sediments are found in a much wider area – in all of the structural depressions in the southern part of the Pannonian basin (Figs. 3 and 4). They are usually coarse-grained clastics deposited in river-to-lake environments. Lacustrine sedimentation was influenced with pyroclastics. Lakes prevailed in the upper part of these successions but they still probably remained isolated until the main marine transgression occurred in Karpatian times. Based on the data from outcrops (e.g. between Mt. Psunj and Mt. Papuk after PAVELIĆ et al., 1998) the Ottnangian lake was gradually transformed into a marine environ- ment during the Karpatian. Karpatian marine sediments – composed of sandy and clayey marl with tuff inter- calations – were found in the Hrvatsko zagorje region In the Pannonian Basin, Early Miocene subsidence and sedimentation began as a result of the first syn-rift extensional tectonic phase (HORVÁTH & ROYDEN, 1981; TARI, 1994; TARI & PAMIĆ, 1998; LUČIĆ et al., 2001; PAVELIĆ, 2001). This extension was also synchronous with the major extension of two microcon- tinents: Alcapa and Tisza. Subsidence from the end of the Karpatian corresponds to the eastward motion and possible rotation of the Tisza unit under a transtensive stress field (CSONTOS et al., 1992, 2002; CSONTOS, 1995). In SW Hungary in the vicinity of the Mecsek Mts. the presence of syn-sedimentary half grabens of Karpatian–Badenian age was proved by several stud- ies (POGÁCSÁS, 1984; BERGERAT & CSONTOS, 1988; TARI, 1992; BENKOVICS, 1997; CSONTOS & NAGYMAROSY, 1998; FODOR et al., 1999). Fol- lowing a short terrestrial period, marine deposition occurred in these depressions connected by the huge water mass of Paratethys. The sea-floor had significant relief: in small, but deep depressions the palaeo-water depth often exceeded 500 m (BÁLDI et al., 2002); these depressions were surrounded by shallow submarine platforms and an archipelago (HÁMOR, 2001). Gener- ally in this time period, bathymetric changes were con- trolled principally by movements of the basin floor. The first compressional event indicates the end of rifting at about the Late Badenian (cf. POGÁCSÁS, 1984; TARI, 1994), followed by only moderate rates of subsidence. The subsequent palaeobathymetric changes, still within the Badenian, were mainly controlled by global sea level variations (BÁLDI et al., 2002). The first early post-rift event of basin inversion took place during the Sarmatian (HORVÁTH & TARI, 1999) resulting in a widespread pre-Pannonian unconformity. Afterwards quiet and relatively slow thermal subsidence occurred (HORVÁTH & ROYDEN, 1981), coupled with uplift and erosion of the surrounding Alpine– Carpathian–Dinaric mountain belt. This supplied huge amounts of sediment via large fluvial to deltaic systems into Lake Pannon (BÉRCZI et al., 1988; SZENT- GYÖRGYI & JUHÁSZ, 1988; JUHÁSZ, 1991, 1994), resulting in the accumulation of post-rift sedimentary successions several thousand metres thick. This unit was subdivided into several 3rd order depositional sequences (POGÁCSÁS et al, 1988; UJSZÁSZI & VAKARCS, 1993; VAKARCS et al., 1994; VAKARCS, 1997; TÓTH-MAKK, 2003), which seemingly are in phase with global sea level changes. Others argue for simple climatic correlation (SACCHI et al., 1999), while a tectonically controlled stratigraphic subdivision is suggested by HORVÁTH & TARI (1999). Because the facies boundaries within the lacustrine succession are time-transgressive (POGÁCSÁS et al., 1988, 1993), a combination of litho-, bio-, magneto- and (seismo-) sequence stratigraphy is used to establish reliable corre- lation between sub basins. Bypassing of the delta front through the northern part of the study area (southern Transdanubia) is estimated to have occurred between 9–7 Ma (Late Pannonian – UJSZÁSZI & VAKARCS, 105Saftić, Velić, Sztanó, Juhász & Ivković: Tertiary Subsurface Facies, Source Rocks... brownish marls and platy marls (LUČIĆ et al., 2001). In the Mura and Drava depressions occasionally very thick sequences of lava and pyroclastics of andesitic to basaltic composition were documented, sometimes in two intervals (Fig. 3). In SW Hungary, the fish-scale bearing lacustrine clayey marl and the Congeria-bear- ing sands represent the transition from Karpatian brack- ish water to Badenian fully marine conditions. These are overlain by coarse clastics (sands, sandstones, and conglomerates) that were deposited from foreshore to nearshore areas (the Budafa Formation and the overly- ing “regressive complex”) and by dark grey, sandy silt- stones and marls (Tekeres Schlier) rich in foraminifera, which formed in an open marine environment (Fig. 2). The thickness of the schlier rarely exceeds 200–400 m, the coeval nearshore clastics may reach 600–700 m (HÁMOR, 2001). The Tekeres Schlier was generally believed to be deposited in a rather shallow environ- ment, referred to mostly as being upper neritic or sub- littoral (HÁMOR, 1970, 2001; SZENTGYÖRGYI & JUHÁSZ, 1988). Recent quantitative palaeoecological analysis of foraminifera, however, proved that it was and in the central and eastern part of the Drava depres- sion (both in wells and at outcrops – PAVELIĆ, 2001). Karpatian marine sedimentation was also documented in the NW marginal part of the Sava depression, on the margins of Mt. Psunj and also in the Požega valley. The first Neogene sediments overlying the Pal- aeozoic–Mesozoic basement in SW Hungary (in the Mecsek Mts. and the surrounding basins both to north, west and south) are Ottnangian to Karpatian terrestrial deposits – conglomerates, sandstones, siltstones, even coal-bearing clay – of large alluvial fans and lakes (Fig. 3). These are interfingered with the “Lower” Rhyolite Tuff which largely contributed to the infilling of the continental basin (HÁMOR, 1970; TARI, 1992). Thick intercalations of Dacite Tuffs (15.3–17.4 Ma – ÁRVA- SÓS & MÁTHÉ, 1992) mark the Karpatian/Badenian boundary. The thickest epiclastic–volcanoclastic suc- cessions can reach 500–1500 m. The Early Badenian sea-level rise covered the Karpatian deposits in Croatia. A pronounced facies differentiation is recorded by interfingering of calcare- nitic biogenic sandstones, turbidites with blue-grey and Fig. 3 Stratigraphic column of the Drava depression. 106 Geologia Croatica 56/1 matically decreased and is reflected in the local deposi- tion of pyroclastics (Fig. 3). Following the Sarmatian, the contemporaneous structural changes in the surrounding Alpine–Carpathi- an mountain belt led to the final isolation of the sea and formation of Lake Pannon (for the best recent summa- ries see MAGYAR & GEARY, 1999; MAGYAR et al., 1999a, b; MÜLLER et al., 1999). Depending on the palaeogeographic position, namely the tectonically determined coeval relief of the basin floor, Lake Pannon sediments either conform- ably overlie Sarmatian layers (Figs. 3 and 4), or they unconformably cover older sediments, or occasionally even the crystalline basement on uplifted blocks. Lake Pannon deposits are described from north to south because this approximates the general direction of sedi- ment transport and the progradation of the facies belts. The southernmost depressions (Sava and partly the Slavonija–Srijem) were partly separated by inselbergs situated perpendicular to the main sediment transport direction (PAVELIĆ, 2001), thus the palaeogeographic situation (in sense of the basin fill directions) becomes more complex to the south. deposited in bathyal depths (BÁLDI et al., 2002). The Late Badenian is represented by Turritella–Cor- bula bearing clay-marls (Szilágy Fm.) – up to 100 m in thickness – indicating a shallow open marine realm. The clay–marl interfingers with the Lithothamnion-bearing lower and upper Leitha Limestones which were formed in the nearshore to littoral zone. Both limestones are made up of carbonate cemented conglomerates, mol- lusc-bearing to oolitic calcarenites and reefs up to about 50 m thick. The Badenian was the last fully marine period in the life of Paratethys. The Late Badenian sea-level rise, as a consequence of re-opening of the Indo-pacific seaway (RÖGL, 1996), caused flooding of uplifted blocks on many places in northern Croatia. The boundary in the base of Upper Badenian deposits may be compared with the acknowledged Middle Badenian unconformity, and is considered as a syn-rift/post-rift boundary (HORVÁTH, 1995; PAVELIĆ, 2001). In northern Croatia, thin Sarmatian deposits are found in conformable successions. The thin-bedded calcareous marls (occasionally with some sandstones) contain a characteristic small-sized brackish faunal assemblage. In addition, varve-like laminites (LUČIĆ et al., 2001; PAVELIĆ, 2001) are also known. Volcanic activity dra- Fig. 4 Stratigraphic col- umn of the Sava depression. 107Saftić, Velić, Sztanó, Juhász & Ivković: Tertiary Subsurface Facies, Source Rocks... Silty clayey marls and siltstones are widespread all over the basin above the turbiditic succession, in con- nection with the approaching delta slope (Algyő Fm. – Fig. 2). In the Drava depression deltaic progradation occurred from the N–NW. Fine grained slope sediments contain up to 2–40 m thick lentoid sandy intercalations of different mass flow origin. The thickness of the strongly progradational delta slope reflects palaeo- water depth (POGÁCSÁS, 1984) and varies from 50 m to several hundred metres. According to reflection seis- mic sections the slope angle also changed from rather gentle and ramp-like, to steep in the deepest depressions (UJSZÁSZI & VAKARCS, 1993). The overlying well-developed, thick sandstones were accumulated mostly on the delta front and delta plain as distributary mouth bar, channel fill or wave- reworked shore-face deposits. Due to the high rate of subsidence a very thick (700–1000 m) aggradational unit (Újfalu Fm.) was formed reflecting several epi- sodes of minor relative lake-level oscillations. The flooded delta plain gradually gave way to an allu- vial plain, on which thin-bedded, often variegated silty clays, occasionally lignites alternate with sand beds a few metres thick (Zagyva Fm.). The fine grained depos- its are of floodplain, wet plain and shallow lacustrine (pond) origin, while sands mainly represent fluvial channel fills. The thickness of this unit is about 300– 400 m in the Drava depression, but elsewhere it was considerably eroded during the Pliocene (Fig. 2). As illustrated in Fig. 4, Late Miocene sediments in the southern part of the study area are generally the same as in the north, but with local differences. A pro- nounced unconformity was identified at the base of the Pannonian sediments, as described from exposures in the region between Mt. Moslavačka gora and Mt. Psunj (BLAŠKOVIĆ, 1982; BLAŠKOVIĆ et al., 1984). In the northern part of the Drava depression, the transgressive deposits usually begin with conglomerates of destructive coasts, fining up into nearshore sand- stones, mainly above the pre-Neogene basement highs (BÉRCZI et al., 1988; SZENTGYÖRGYI & JUHÁSZ, 1988). Above this thin coarse-grained veneer, dark grey calcareous marls that gradually turn into clayey marls (Endrőd Fm. – Fig. 2) of hemipelagic origin are found basin-wide, indicating a longer period of basin starvation. The clayey marl is already an indicator of the slowly approaching siliciclastic input. The colour of the marls varies from black to brown indicating an alternation of anoxic to dysoxic bottom conditions. It was found that the colour lightens towards the base- ment highs (SZENTGYÖRGYI & JUHÁSZ, 1988). It is difficult if not impossible to estimate water depth by palaeoecological tools because a completely endemic mollusc and ostracod fauna developed in Lake Pannon (JUHÁSZ & MAGYAR, 1992; MÜLLER et al., 1999). Despite these difficulties, the depositional environ- ment is thought to be generally deep during this period, particularly because several series of fine-grained sand- stones of turbiditic origin can be found in the imme- diate overlying sequence (Szolnok Fm.). These are supposed to be formed mainly on turbidite fan lobes and sheets (BÉRCZI et al., 1988; POGÁCSÁS et al., 1988; SZENTGYÖRGYI & JUHÁSZ, 1988; JUHÁSZ, 1994). The thickness of individual sandstone bodies usually varies in the range of 3–10 m, and they are separated by 2–20 m thick intercalations of marls. The thickness of the whole turbiditic succession, however, shows great variation depending on the topography of the basin floor: it pinches out towards the basin margins and is significantly thinner above basement highs, while it can reach 1500 m in the deepest parts of the depres- sions (Fig. 5). Fig. 5 Geological cross-section based on correlation of well-log data (location in Fig. 1). 108 Geologia Croatica 56/1 SW Pannonian basin where sedimentation continued in the same environments (Figs. 3 and 4). Infilling of the remnants of Lake Pannon occurred with a variety of coarse clastics mixed with clay, and lignite seams formed in delta plain marshes. These were followed by meandering river systems. As climatic change coupled with basin inversion and uplift of mountainous areas took place in the Quaternary (FRISCH et al., 1998), not only sandstones but also gravel beds appeared as channel fills. The oldest gravels in the Mura and Drava river valley are the Belvedere gravels (oldest Qua- ternary or Pliocene–Quaternary) while the Sava river brought mainly carbonate-clasts as sandy gravels only in the Holocene. During the Pleistocene glacial periods loess mainly accumulated in the area and there are also aeolian quartz sands deposited in the western part of the Drava depression. Lacustrine–marsh sediments are characteristic for the warmer periods – more detailed recent descriptions of these sediments can be found in the work of VELIĆ & SAFTIĆ (1991), VELIĆ & DURN (1993), VELIĆ et al. (1999) and BAČANI et al. (1999), where results of mainly hydrogeological explo- rations are used for stratigraphic interpretation. 4. NEOGENE SEDIMENTARY MEGACYCLES The Neogene to Quaternary basin fill in the SW mar- ginal parts of the Pannonian basin usually has an aver- age thickness in the range of 500–1500 m, except in the central parts of depressions where much bigger values – around 4000 m in the Slavonija–Srijem depression, 5500 m in the Sava and Mura depressions, and close to 7000 m in the Drava depression are found (Fig. 10). This huge sedimentary pile can be subdivided into three main sequences. The cyclic character of the Neogene sedimen- tary successions of the region was recognized long ago (FILJAK et al., 1969; HÁMOR, 1970) and the three units were later defined by ŠIMON (1980) as the “macro-rhythms of sedimentation”. The comparable units were later named “megacycles” (VELIĆ et al., 2002), as their inferred time-span is of 6.8, 5.9 and 5.6 Ma respectively (measured on the geochronologi- cal scale of RÖGL, 1996; cf. Figs. 2, 3 and 4). Each of these megacycles encompasses rocks deposited during 2nd order transgressive–regressive cycles, influenced mainly by tectonics (VAIL et al., 1991). The megacycles are separated by major unconformi- ties, as follows: – Base Neogene unconformity, which usually sepa- rates the historically complex basement from either the Early, Mid- or even Late Miocene sediments. On seismic sections it is an onlap surface (Figs. 6, 8 and 9); in outcrops or borehole successions it is usu- ally covered by a transgressive lag of coarse clastics of various ages. This is usually regarded as the syn- rift unconformity (cf. HORVÁTH & TARI, 1999), Lower Pannonian sediments – mainly calcareous marls and limestones (Prkos or Valpovo Fm. – Fig. 2) – were deposited in a shallow lacustrine environment (Figs. 3 and 4), but there are also minor deltaic bodies along the faulted margins of depressions. The composition of these sandstones is closely associated with the underly- ing rocks so they are interpreted to be locally derived. In the western part of the Drava depression in Croatia (Fig. 3) this is already the time when the first turbidite sandstones were transported from the NW (Koprivnica Ssts. – Fig. 2), found to be lithologically identical with the overlying units. The following, mainly Upper Pan- nonian to Pontian sandstone–shale sequence (Ivanić Grad, Kloštar Ivanić and Široko Polje Formations in the Sava depression – see ŠIMON, 1973, 1980) was depos- ited in locally very deep troughs, especially in the west- ern part of the Drava and Sava depressions. The troughs were largely isolated from the rest of the Pannon- ian basin by islands or by subaqueous basement highs. These sediments are of fairly uniform lithological com- position: a sequence of sand/sandstone (subgreywacke and calcarenite–subgreywacke) bodies interlayered with silts and marls (Fig. 2). Sands and sandstones are grey-coloured, with a predominance of quartz grains, and some lithoclasts (usually limestone), platy miner- als and feldspar grains. The mineralogical composition of the sandstones (subgreywackes) is almost the same and is interpreted as being a consequence of an external and distant sedimentary source (in contrast to the older Miocene sediments – ŠĆAVNIČAR, 1979). The mor- phology of these multiple sandstone bodies indicates a turbidite fan to delta system. Sedimentary transport from the N–NW resulted in progradation of firstly the turbidites and then the prodelta, delta slope to delta front clastics, with the aerial distribution of lithofacies largely influenced by the elongated shape of the depo- sitional basins. A depositional model was constructed with the source of sediments in the Alps, and transport by turbidity currents (ŠIMON, 1980). In this process, coarser grains were deposited in the deepest parts of the depressions and fine-grained material at the base- ment highs. In the Upper Pannonian sediments of the western Sava depression this was later confirmed by the investigations of VRBANAC (2002a, b), substantiated by a large data set of core investigations. Analysing the aerial distribution of the Lower Pontian turbidites in the western part of the Sava depression SAFTIĆ (1998) confirmed that regionally recognisable E-log markers delimit numerous sandstone members as genetic units. The markedly uniform petrographic composition of both sandstones and marls, together with their sedimen- tological characteristics, are congruent with the external sedimentary source. More detailed dating of high- resolution units still remains a problem nevertheless, it appears that the deltaic environment was formed almost ubiquitously in the Late Pontian, except in the deepest parts of the Sava depression. During the Pliocene and partly in the Pleistocene, there remained several deep or shallow lakes in the 109Saftić, Velić, Sztanó, Juhász & Ivković: Tertiary Subsurface Facies, Source Rocks... however, the onlapping younger deposits (e.g. Late Miocene) indicate that this surface is often a compos- ite unconformity. The following Basal Late Badenian unconformity is usually considered to be the syn-rift/ post-rift boundary (ROYDEN & HORVÁTH, 1988; TARI, 1994; HORVÁTH, 1995; VAKARCS, 1997; PAVELIĆ, 2001). It is not as straightforward to iden- tify seismically as the younger one. – The Base Pannonian is the next major unconform- ity. Lake Pannon deposits are onlapping either the basement (Figs. 5 and 9) or older Miocene deposits (Figs. 6–8) with a transgressive lag (POGÁCSÁS et al., 1988; SZENTGYÖRGYI & JUHÁSZ, 1988; HÁMOR, 2001). Seismic sections clearly reveal the erosional character of this surface, as it cuts the previ- ously tilted Badenian–Sarmatian successions (Figs. 5 and 8). The Base Pannonian unconformity is a result of uplift and erosion caused by the first compression- al (inversion) event, proved both by seismic and out- crop-scale structural data (TARI, 1994; HORVÁTH & TARI, 1999; CSONTOS et al., 2002). – The Base Pliocene unconformity seems to be the result of another compressional event, generating both uplift of the pre-existing basement highs with the overlying Lake Pannon deposits, formation of new highs and subsidence between these areas. Thus sediments above this unconformity are onlapping in character (Figs. 5, 6, 8 and 9). In areas of continu- ous subsidence an uninterrupted Miocene–Pliocene sedimentation is recorded. Even a fifth, Quaternary unconformity exists, as the reflection of the youngest basin-wide and local Quaternary uplift/compressional event (DUNKL et al., 1994; TOMLJENOVIĆ & CSONTOS, 2001; CSONTOS et al., 2002). In some places it can be well observed on seismic sections (Fig. 6). A description of the main characteristics of distribution/composition of rock units grouped in the three megacycles follows, illustrated by thickness maps (Figs. 10–13). 4.1. Terrestrial to marine deposits (Base Neogene – top Sarmatian) The succession from the Base Neogene unconformity to the erosional top of the Sarmatian sediments is referred to as the 1st megacycle, which is the sedimentary response to syn-rift and early post-rift structural evolu- tion of the Pannonian Basin. On the thickness map of this unit (Fig. 11) the darkest spots mark regions where all the underlying older formations are exposed. These range from the oldest, Palaeozoic crystalline rocks, through a variety of Mesozoic formations (mostly in the inselbergs) to the youngest pre-basinal unit – Creta- ceous–Eocene flysch in the contact zone with the Inner Dinarides (W of Varaždin – “Vz” and SE of Karlovac – “Ka”). Areas where the thickness exceeds 1 km are Fig. 6 Section A – western Sava depression (see location in Figs. 10–13). 110 Geologia Croatica 56/1 5 AVANIĆ, R., PAVELIĆ, D. & KOVAČIĆ, M. (1995): Građa površin- skih ekvivalenata formacija Ivanić-Grad, Kloštar Ivanić i Široko Polje.– Unpubl. report, INA–Naftaplin Archive, Zagreb. 6 PIKIJA, M. & ŠIKIĆ, K. (1996): Stratigrafsko-facijesna interpreta- cija naslaga miocena (zaključno sa sarmatom) Dravske potoline.– Unpubl. report, INA–Naftaplin Archive, Zagreb. 7 PIKIJA, M. & ŠIKIĆ, K. (1999): Stratigrafsko-facijesne odlike sta- rijeg i srednjeg miocena Murske depresije i obodnih područja.– Unpubl. report, INA–Naftaplin Archive, Zagreb. The deepest troughs, thus the thickest successions, are found in the northern part of the Drava depression (Fig. 11, also sections in Figs. 5 and 9). These half grabens are thought to be associated with stretching of large lithospheric units during eastward motion and rotation (CSONTOS et al., 1992; CSONTOS, 1995; CSONTOS & NAGYMAROSY, 1998). Mainly neu- tral–acidic volcanism was related to the stretching, filling the half grabens with lavas and volcanoclastics up to a thickness of 1500 m. A series of partly isolated depocentres scattered along the southern margin of the present Drava depression and along the northern margin of the Sava depression are also remarkable. These occur along the strike of the major longitudinal fault systems along which rifting started in Ottnangian times (TARI, 1994; HORVÁTH, 1995; TARI & PAMIĆ, 1998; HOR- VÁTH & TARI, 1999; FODOR et al., 1999; PAVELIĆ, 2001). The strong influence of tectonics on sedimenta- tion in the syn-rift phase resulted in a variety of depo- sitional settings and lithofacies. Obviously, 3rd order depositional sequences – not to be discussed here – are tectonically overprinted (cf. BÁLDI et al., 2002). In the Croatian part of the study area, lithological composition and depositional characteristics of the 1st megacycle are described based on the maps and reports from the INA archives (AVANIĆ et al., 19955 and PIKIJA & ŠIKIĆ, 19966, 19997). Ottnangian freshwater coarse clastics were found throughout the Drava depres- sion, all around the Slavonian Mts. (Papuk, Krndija and Psunj), and also in the western Sava depression. Karpa- tian marine sediments were documented in the eastern Drava depression, at outcrops in practically all moun- tains of northern Croatia and NW of the Mecsek Mts. in Hungary. At about the Karpatian/Badenian boundary, huge areas became inundated by normal salinity sea water and the former terrestrial depocentres became the loci of the deepest (bathyal) basins. Due to this Early Badenian sea-level rise the majority of the 1st megacycle consists of Badenian rocks. Upper Badenian coastal to shallow marine deposits are found around the coeval basement highs: Lithothamnion and bioclastic limestones prevail in the western part of the study area; conglomerates in the centre and extending to the Slavonija–Srijem depression; while marls are found on the slopes and in the central parts of all the depressions. The youngest part of the 1st megacycle is the Sarmatian, mainly con- formably overlying Badenian sediments in many places throughout the study area. The basin was then starved loci of the main depocentres and are shaded by grades of light grey. Fig. 7 Section B – central Sava depression (location in Figs. 10–13). 111Saftić, Velić, Sztanó, Juhász & Ivković: Tertiary Subsurface Facies, Source Rocks... Fig. 8 Section C – Bjelovar sag and Drava depression (location in Figs. 10–13). Fig. 9 Section D – northern part of the Drava depression (location in Figs. 10–13). u t:: o .... o CD rJ) (s) .lM.l z o c: o .... o Q) en (S)lMl 112 Geologia Croatica 56/1 graphic units which start with turbidites, followed by progradation of deep-water deltas gradually shallowing up to the alluvial sediments. 4.3. Inversion-related lacustrine to fluviatile deposits (Base Pliocene – surface) Sediments of the 3rd megacycle are found only in those depressions which are situated between the most elevated basement units. They begin with lacustrine deposits and marshes in the remnants of Lake Pan- non in Pliocene times. Fluvial beds become more frequent in the Pleistocene (in the Drava river val- ley) and Holocene (Sava river valley). The end of the cycle is marked by gravels, loess and unconsolidated surface layers (Figs. 2–4). As illustrated in Fig. 13, the thickness of this unit is negligible in Hungary (below 200–300 m), but reaches well over 1000 m in several places to the south. Local thickness variations can be quite abrupt, producing differences in relatively closely spaced seismic sections (compare the NE end of Section C in Fig. 8 with the southern end of Section D in Fig. 9). In the central part of the Sava depression, over 1500 m were drilled and the sequence reaches 1800 m in the central Drava depression between the Slavonian Mts. and the Mecsek–Villány Mts. Contacts with underly- ing rocks are frequently marked with onlaps (Section A, Fig. 6; Section D, Fig. 9), and there are regions where these sediments cover all the older Neogene rocks and directly overlie Palaeozoic rocks (N and W slopes of Mt. Moslavačka gora). As a rule, in comparison to the previous megacycle the depocentres are shifted in a southerly direction (western Sava) or to the SE (Bjelovar sag and central Drava depression). Different relationships were determined in the eastern part of the study area where areas of the maximal thickness of the 2nd and 3rd megacycles either overlap each other (east- ern Drava and Slavonija–Srijem depression), or there is a local shift of the depocentre in a westerly direction (eastern Sava – S of Mt. Psunj). 5. HYDROCARBON POTENTIAL 5.1. Source rocks within the Neogene basin fill Studies of hydrocarbon source rocks in the Pannon- ian basin were intensified in the last two decades of the 20th century. A number of papers were published both in Croatia and in Hungary by the joint efforts of organic geochemists and petroleum geologists (BARIĆ & RADIĆ, 1988; DANK, 1988; HORVÁTH et al., 1988; KŐRÖSSY, 1988, 1989; BARIĆ et al., 1989; PUTNIKOVIĆ et al., 1989; RADIĆ et al., 1989; ALAJBEG et al., 1990; SZALAY & KONCZ, 1991; CLAYTON & KONCZ, 1994; HERNITZ et al., 1995a; HORVÁTH & TARI, 1999). Neogene source rocks with type II and III kerogen content were deposited as marine or lacustrine clays from the sediment supply and was later strongly eroded on some elevated blocks, so the Sarmatian unit is rarely thicker than 50 m. These layers were mostly deposited in a shelf environment with occasional coastal influence along the southern margin of the Sava depression and NE of Mt. Krndija. The offshore Sarmatian sediments were documented in the southern and western mar- ginal parts of the Pannonian Basin, i.e. NE of Karlovac (“Ka”), W of Krapina (“Kr”), as well as in the eastern part of the Sava depression. During the Late Badenian very rapid subsidence ceased, only slow post-rift subsidence occurred. At the end of the Sarmatian the first very pronounced compressive event resulted in rejuvenation of former extensional structural elements in a compressive regime (cf. TARI & PAMIĆ, 1998; HORVÁTH & TARI, 1999; CSONTOS et al., 2002). Some basement blocks were uplifted and to a various extent erosion of the syn-rift and early post-rift strata took place (Figs. 6 and 8). 4.2. Lacustrine, deltaic to alluvial deposits (Base Pannonian – top Pontian) Sediments of the 2nd megacycle were all deposited in Lake Pannon during an interval of thermal subsidence that followed rifting (HORVÁTH & ROYDEN, 1981; ROYDEN & HORVÁTH, 1988; HORVÁTH & TARI, 1999). As subsidence commenced after the basin inver- sion, the “lacustrine” transgression quickly inundated both the wide depressions and the basement highs. As a result, a large underfed lake with an inherantly uneven basin floor topography was born (ŠIMON, 1973, 1980; BÉRCZI et al., 1988; KŐRÖSSY, 1989). The thickness of the 2nd megacycle (Fig. 12) varies from 2000 m (western Sava depression) through 3000 m north of the Iharosberény high, to 4–5000 m along the SW margin of the Drava depression. As a rule, the maximum thick- ness is close to basin margin faults. These are also the places where sandstones of turbiditic origin show their maximum thickness (cf. Fig. 5). Turbidites, accumulat- ing mainly in the deepest depressions, particularly in front of, and on the slopes of basement highs, gradually levelled the basin floor topography. Above the turbidites deltaic progradation fed from the N, NW occurred. Due to persisting high rates of subsidence coupled with enormous rates of sediment supply, delta front-, delta plain-, coastal plain- to alluvi- al plain-series accumulated mainly in an aggradational style (Figs. 2–4), resulting in a relatively large uniform thickness. In the northern part of the Drava depression, the top of the successions infilling Lake Pannon were eroded mainly above the previously existing base- ment highs (Figs. 5 and 9), suggesting rejuvenation of former structural elements again in a reverse sense (cf. CSONTOS et al., 2002). According to the present regional seismic interpretation in the Sava depression (IVKOVIĆ et al., 2000), the entire Pannonian to Qua- ternary sequence can be subdivided into 6 seismostrati- 113 S aftić, Velić, S ztanó, Juhász & Ivković: Tertiary S ubsurface Facies, S ource R ocks... Fig. 10 Isopach m ap of the entire N eogene to Q uaternary sequence w ith H C accum ulations. on 171 ~ 6300000 , Thickness of Ng+Q basin fill ~2;1'\ 5~~o/ 0==~ (depth of pre-rift Palaeozoic-Mesozoic basement) ~-~~'-1.S ........ - ............. + \ ~-) o Scale 50km Co~ined seismic sactiDn~"-"'-, I I borders I m I Oil field I c;::,1 Oil and gas field I I Gas field 6400000 6450000 6500 000 114 Geologia Croatica 56/1 115Saftić, Velić, Sztanó, Juhász & Ivković: Tertiary Subsurface Facies, Source Rocks... and marls. The source rocks are dark grey to black marls, calcareous marls, clayey limestones, siltstones and shales. In these sediments the organic matter usually appears uniformly dispersed in the matrix, but it can also be found in the form of lenses, bands or as centimetre to millimetre thick laminae. These source rocks were deposited during the 1st and in the beginning of the 2nd megacycle – mostly between the Base Neogene and the Lake Pannon turbiditic successions. They were formed mainly during the Karpatian and Badenian, and a minor portion during the Sarmatian and Early Pannonian (Figs. 3 and 4). According to SZALAY & KONCZ (1991), in the Hungarian part of the basin, the oil generation window is situated between 2.4–4.3 km, while TOC averages 1–2 wt.% (up to 5 wt.%). There are type II and III kero- gens, thus both oil and gas generation have occurred. In the southern part of the study area (in Croatia) there is a local variation in composition of organic matter. Two main types appear – either the hydrogen-rich kerogen type II (liptinite) mixed with some type I (alginite), or kerogen type III (vitrinite). Both types originate from the bacterially altered marine algal–lipid matter more or less mixed with terrigenous matter and depos- ited in the dysoxic–anoxic environment of the bottom waters. The source rocks are mostly mature, because they reached the conditions for expulsion at depths exceeding 2000 m and temperatures of over 100ºC. Oil–Source correlation determined that the C13/C12 isotope ratio of the above source rocks has the same value in oils of the area, which means that the analysed hydrocarbons really originate from the described source rocks (PUTNIKOVIĆ et al., 1989; RADIĆ et al., 1989; HERNITZ et al., 1995a). When compared to the pro- duced quantities and the total HC reserves, the Mid- dle Miocene source rocks are evaluated as being thick enough and having a very good generation potential. The situation in the western part of the Drava depression in Croatia is different because of the pres- ence of large gas and gas condensate fields (Fig. 10). These are the deepest reservoirs of the Pannonian basin in Croatia. According to BARIĆ et al. (1992) the sediments rich in organic matter are found there in two intervals. The older one is of Early Miocene age (Fig. 3) and contains siltstones and mudstones at a depths of 3200–4500 m. Hydrogen indices are significantly reduced – this is the kerogen type III, and rocks are presently in the phase of increased thermal maturity (mature to over-mature) and capable of generation of dry and wet gas and gas condensate. It was found that part of the source rocks from this interval has generated only dry gas, which is thermally highly altered. The younger interval of source rocks in this area is formed of Badenian sediments and Pannonian fossiliferous calcareous marls at a depth of 2000–3000 m. Opti- cal analyses of the maceral composition of kerogen, revealed predominance of vitrinite and debris of wood structure, as well as sporadic concentrations of exinite which, of course, originate from the land, but which were deposited in lacustrine anoxic conditions. These rocks are presently situated within the oil window. There are also older, pre-Neogene source rocks in the nearby parts of the Pannonian basin. In the Hungar- ian Zala basin, in the NW part of the study area, two– three types of oils were identified by isotopic, chemical and biological signals, which indicate that the Upper Triassic black shales (Kössen Marl) were the main source rocks, and that Miocene source rocks produced only minor amounts. Oil generation and expulsion both from Triassic and Miocene source rocks began during the Late Miocene extension (CLAYTON & KONCZ, 1994). In the Drava basin this Kössen Marl cannot be found in the pre-Neogene basement, therefore the main source rocks must be of Neogene age, but that does not mean that similar future findings can be ruled out, espe- cially in the contact zone of the Dinarides and Eastern Alps. In the same area, there are the remnants of the Slovenian Trough with lots of flysch sediments and also Cretaceous–Palaeogene flysch sediments under the most western part of the Sava depression – but this is beyond the scope of this paper. 5.2. Main petroleum reservoirs Matured hydrocarbons migrated along the major uncon- formities and fault zones from the deepest parts of the depressions towards the structural traps situated usually above the basement highs or in en-echelon folds in the vicinity of the most faulted margins of the depressions. Neogene sand beds in folded and compacted anticlines over basement highs are the most obvious hydrocarbon traps (KRANJEC, 1972; DANK, 1988). Occasion- ally the fractured basement together with the gener- ally overlying basal conglomerates bears hydrocarbons. Stratigraphic traps are difficult to find and only minor occurrences of this type have been discovered. Most of the reservoirs are within the 1st and 2nd megacycles – Ottnangian breccias, Karpatian, Badenian or Pannonian basal conglomerates and Lake Pannon prodelta turbiditic lobes, sand sheets, slumps, delta front mouth bars and delta plain channel fills have the most favourable reservoir characteristics (Figs. 3 and 4). The seals are usually provided by the associated intralobe mudstones or interdistributary clays. The most important lithostratigraphic units are illustrated in Fig. 2, details can be found in the works of ŠIMON (1973) and KŐRÖSSY (1989). The oldest coarse-grained clastic member of the Ivanić Grad Formation – IVA Sandstones, contains 30% of the total discovered HC reserves in Croatia. Most thoroughly investigated in the Sava depression (VRBANAC, 2002b), it was depo- sited in a deep basinal, prodelta environment. This finding can be extrapolated in other sub-basins. More detailed descriptions can be found in the works of ĐUREKOVIĆ (1995), KOVAČEVIĆ (1996), LARVA (1996) and MATASOVIĆ (1996) and also the first part of this paper. The most important characteristic for HC accumulation in this unit is that porosity varies from 114 G eologia C roatica 56/1 115 S aftić, Velić, S ztanó, Juhász & Ivković: Tertiary S ubsurface Facies, S ource R ocks... Fig. 11 Isopach m ap of the 1st m egacycle. ~ § § ~ .... - ~ J I~' ( ;iI''' ZU.,. , 0 !.. Kr .... :, ... ,. ! ~I ~, .... C<:\\, Con1tJined seismic sections~ 6350000 '. ". r--~--.: : .f' J. .(" ,.', .' ,i ;. ...... .,.... ...... -:~/. (~:/.);. ,. 6400 000 6450000 10t megacycle Terrestrial to marine deposits (Base Neogene - Base Pannonian) .. I' 6500000 • rC;{--" .,.,-t;~ (,. ..... ..... _! ......... /'. , {;~;! '. C', ,-. \ i ....... :. \. f' . . -.- 6600 000 116 G eologia C roatica 56/1 117 S aftić, Velić, S ztanó, Juhász & Ivković: Tertiary S ubsurface Facies, S ource R ocks... Fig. 12 Isopach m ap of the 2nd m egacycle. 8 ~ ~ o Scale SOlan ! ! 6400000 6450000 2nd megacycle Lacustrine - deltaic - alluvial deposits (Base Pannonian - Base Pliocene) .. .. "" 6550000 116 G eologia C roatica 56/1 117 S aftić, Velić, S ztanó, Juhász & Ivković: Tertiary S ubsurface Facies, S ource R ocks... Fig. 13 Isopach m ap of the 3rd m egacycle. !5 ~ Bj 0 8 ? Scai,e ~Okm borders 6 350 000 6400 000 6450000 3'd megacycle Inversion-related lacustrine to terrestrial deposits (Base Pliocene - surface) .. • '" ~ ~~ o. o ------------o.~ c::-'-:;---, . , ~ o· ; - \.,1 "'; ..... _! ......... /'. I~?~J '-, \., .... ~.,. ( ... ~,._.ri._ ..... , dJ{"---' ". , (~ 0'-'_. ' . 2P 9 6500000 ._-.... .;;~: ... ;: ,~j 0' .... _ . ..,. " ._.l'.~.~. + "'o.s" 6 550 000 8 ..... _ ..... ) -\ > .. --- .... 000 118 Geologia Croatica 56/1 119Saftić, Velić, Sztanó, Juhász & Ivković: Tertiary Subsurface Facies, Source Rocks... 11% to 25% and even 30%. It was also determined that porosity declines with depth, mainly because of compaction and because the ratio of the silt-sized grains increases in this direction. The horizontal perme- ability was also measured on cores and it is between 0.1x10-3m2 and 380x10-3m2. Reservoir characteristics usually vary over relatively short distances, which must be precisely defined in order to obtain a better recovery ratio, especially in the phase of secondary production (injection of water or gas) or in the planned tertiary pro- duction methods. The Drava depression is both the largest and the deepest sub-basin in the study area. The deepest part lies between Virovitica and Slatina (“Vi” and “Sl” in Fig. 10), where the thickness of Neogene and Quater- nary sediments can reach 7000 m. In the Croatian part of the Drava depression, 30 gas and oil fields were discovered, of which 23 are presently producing oil, gas condensate or gas. In the western part of the Drava depression, the HC traps are within conglomerates and sandstones of the Lower and Middle Miocene (1st megacycle) and Pannonian sandstones (basal part of the 2nd megacycle). Above, the Lower Pontian sandstones also have good reservoir properties. Nevertheless, the most valuable hydrocarbon accumulations are in the fractured Palaeozoic and Mesozoic rocks with second- ary porosity; the reservoirs at Molve and Kalinovac (“M” and “K” in Fig. 10) are the largest gas and gas condensate fields in Croatia. In the eastern part of the Drava depression, the largest HC accumulation is the Beničanci oil field (marked with “B”, where more than 18x106m3 of oil has been produced since 1972). Here the reservoir rocks are Middle Miocene basal carbonate breccias with excellent petrophysical properties (HER- NITZ et al., 1995b). According to TIŠLJAR (1993), the Beničanci reservoir is composed of wedge-shaped bodies of rockfall and debris breccias, formed by the accumulation of large masses of dolomite detritus from the shore zone. The cliffs of the steep shoreline were situated along the margins of the faulted dolomite mas- sif that was gradually uplifted. In this unit there are also thinner bodies of breccias, conglomerates and sand- stones built up by detritus of coral reefs (VELIĆ et al., 2000). They were also deposited in the shoreline zone during the Early Badenian sea-level rise and subsequent Late Badenian transgression. On the Hungarian side of the Drava depression two basement highs are of primary interest: the Inke–Iha- rosberény high at the north-western part of the basin and the Görgeteg–Babócsa high to the south (Fig. 10). Görgeteg is also a historical place, because the first well log measurements in Hungary were carried out here in 1935. In this area both oil and gas bearing rocks are confined to Lake Pannon turbiditic sand sheets and in one example to a delta plain channel-fill sandstone. The excellent reservoir quality sandstones of the delta front and delta plain hardly ever contain traps because there are no good seals in the extremely thick wave- reworked, delta-front successions. Not far from Gör- geteg a minor oil pool was also found in the Badenian Leitha (Rákos) Limestone, and north of the high, minor indications of oil reservoirs were detected in Badenian (tuffaceous) sandstones, and conglomerates of shal- low marine origin. South of the main Görgeteg high, a minor anticline was found, in which gas reservoirs are situated at the top of the Palaeozoic basement and in the overlying Badenian transgressive coarse clastics. At the Inke high small isolated sand bodies of Badenian to Lower Pannonian are producing minor amounts of gas. In addition, occasionally the top of the fractured Trias- sic carbonates is a CO2–CH4 reservoir. Between the two major highs, the Somogyudvarhely–Gyékényes depres- sion is found, which acted rather as a source area; here both the fractured crystalline basement and the overly- ing Karpatian to Badenian sandstones contain indica- tions of oil and CO2 (KŐRÖSSY, 1989). The Mura depression is in the NW part of the study area. The most valuable hydrocarbon reservoirs in the southern part of it occur in breccias, conglomerates and sandstones of the Early and Middle Miocene, and Lake Pannon sandstones of different depositional facies. Reservoir rocks were deposited in a variety of depo- sitional settings, ranging from the subaqeous slumps and turbidites (Ottnangian–Badenian) to flood plain and slope lithoarenites and biocalcarenites, mixed with siliciclastic detritus that were formed in the vicinity of patch reefs (mostly Badenian). The described sediments form HC reservoirs of several smaller and two larger fields – Veliki Otok and Legrad (“VO” and “L” in Fig. 10) which, according to the cumulative production, are amongst the 10 largest gas fields in Croatia. Similarly to the Hungarian part of the Drava depression, quite large, but CO2-rich gas accumulations were recently discovered in the Triassic basement highs in the very NW corner of the study area, close to the border with Slovenia. The Slavonija–Srijem depression in the SE corner of the study region is relatively shallow and covers the smallest area. Only the western part lies in Croatia, where the entire Neogene–Quaternary sedimentary and volcanic succession has a maximum thickness of 4000 m (Fig. 10). Three HC fields were discovered in this area, with relatively shallow reservoirs (approx. 1000 m) formed of the fractured and weathered Palaeozoic and Mesozoic rocks together with overlying Lower and Middle Miocene clastics (1st megacycle). The Sava depression lies along the SW margin of the Pannonian basin and the thickest sequence of the Neogene–Quaternary sediments (over 5000 m) can be found in its central part, S of Mt. Moslavačka gora (Fig. 10). The most important reservoir rocks are Pannonian and Pontian sandstones which start with the transgres- sive locally sourced sandstones, soon pass into tur- bidites and end with large masses of deltaic sediments (mostly sandstones). There are HC reservoirs also in the weathered and fractured Palaeozoic granite and gneiss of the covered western slopes of Mt. Moslavačka gora, 118 Geologia Croatica 56/1 119Saftić, Velić, Sztanó, Juhász & Ivković: Tertiary Subsurface Facies, Source Rocks... but the most important fields are in Neogene reservoirs – Stružec, Žutica and Okoli fields (marked with “S”, “Ž” and “O” in Fig. 10). 6. CONCLUSIONS 1) In the Neogene sedimentary succession three depo- sitional megacycles can be identified in the SW part of the Pannonian Basin. Each encompasses rocks deposited during 2nd order cycles influenced mainly by the structural evolution of the basin. The megacycles are separated by major unconformities – the Base Neogene, the Base Pannonian and the Base Pliocene unconformities, which usually separate the onlapping deposits of the younger cycle from uplifted, tilted and eroded older rocks. These uplifts interrupting basin subsidence reveal compressional events in association with inversion tectonics. Such a subdivision is suitable for the subsurface geological data, in the sedimentary record the actual syn-rift/post-rift boundary is identified as the prominent unconformity in the base of the Late Badenian sediments. Even a fifth, Quaternary unconformity exists, as the reflection of the youngest uplift of the mountainous areas. 2) The 1st megacycle is composed of Early to Middle Miocene syn-rift and early post-rift sediments. Terrestrial sandstones, subordinate coal seams, coastal talus breccias, conglomerates and sandstones are overlain by shallow to deep marine marls and shales of varying carbonate content and with thin sandy intercalations, while biogenic and bioclastic limestones were formed in the nearshore zone. The end of the megacycle is characterised mostly by fine-grained deposition in the starved brackish- water basin. 3) The 2nd megacycle is connected to the thermal subsidence of the Late Miocene Pannonian basin, thus is built up by the post-rift sediments of Lake Pannon. Apart from local variations due to pre- Pannonian topography the sedimentary succession begins with littoral limestones and a nearshore trans- gressive lag overlain by hemipelagic calcareous and clayey marls basin-wide. The deepest depressions were filled by lacustrine turbidite lobes and channel fills of considerable thickness, thus the initial basin floor topography gradually became levelled. Turbiditic successions are overlain by shale-prone delta slope and sandy delta front to coastal plain sediments. The northern part of the Drava depression was filled to lake level and the youngest deposits of this megacycle were formed in fluvial environments. To the south, deposition continued on the delta front and in the prodelta region respectively. 4) The 3rd megacycle was formed during Pliocene to Quaternary basin inversion resulting in subsidence in the deepest zones, associated with uplift and erosion of the most elevated blocks. At that time Lake Pannon occupied only the southern part of the study area where basin fill up to fluvial deposits continued in a similar manner as previously occurred. 5) The main source rocks with type II and III kerogens are Middle Miocene marine clays and marls and Upper Miocene hemipelagic marls of Lake Pannon. They generally occur in the oil generation window or in the wet and dry gas zone depending on their present-day position in the basin. 6) The main hydrocarbon reservoirs are above struc- tural highs in the studied basins, either in the frac- tured Neogene basement, associated with basal conglomerates and sandstones, or in the Neogene sedimentary succession. Lower to Middle Miocene coarse clastics – rockfall breccias composed mostly of Mesozoic dolomite fragments, conglomerates and sandstones accumulated in the coastal to shallow marine zones during the Badenian, form reservoirs with good characteristics. The most important res- ervoirs, however, are sandstone bodies of differ- ent (turbiditic or deltaic) origin within the Upper Miocene 2nd megacycle. Their reservoir character- istics are highly variable. Acknowledgements We are indebted to László CSONTOS for discussions on structural events and also for providing some unpub- lished data. Imre MAGYAR is thanked for helping us with the Lake Pannon stratigraphy and for having read the manuscript with the utmost care. Davor PAVELIĆ and Boris VRBANAC kindly helped us in correlation/ interpretation of the data acquired from the surface with the subsurface data sources. For discussions regarding the Pannonian basin development in Croatia we are thankful to Jožef ŠIMON and Eduard PRELOGOVIĆ. Financial support was provided by the Croatian Min- istry of Science and Technology (Scientific Project no. 0195035) and by Hungarian Scientific Research Grants (OTKA) nos. T 037724 and T 035168. 7. REFERENCES ALAJBEG, A., BRITVIĆ, V., ŠVEL-CEROVEČKI, S., CORNFORD, Ch., TODORIĆ, A., RAJKOVIĆ, J., BARIĆ, G. & PUTNIKOVIĆ, A. (1990): Geochemi- cal study of the oils and source rocks in the Pannonian Basin.– Organic Geochemistry, 16, 339–352. ÁRVA-SÓS, E. & MÁTHÉ, Z. (1992): Mineralogical and petrographic study of some Neogene layers of the Mecsek Mts. (south Hungary) and their K–Ar dating.– Földtani Közlöny, 35/2, 177–192. BAČANI, A., ŠPARICA, M. & VELIĆ, J. (1999): Quaternary deposits as the hydrogeological system of Eastern Slavo- nia.– Geol. Croatica, 52/2, 141–152. 120 Geologia Croatica 56/1 121Saftić, Velić, Sztanó, Juhász & Ivković: Tertiary Subsurface Facies, Source Rocks... BÁLDI, T. (1980): The early history of Paratethys (in Hungar- ian with English abstract).– Földtani Közlöny, 110/3–4, 456–472. BÁLDI, K., BENKOVICS, L. & SZTANÓ, O. (2002): Bade- nian (Middle Miocene) basin development in SW Hunga- ry: subsidence history based on quantitative paleobathym- etry of foraminifera.– Int. J. Earth. Sci., 91, 490–504. BARIĆ, G. & RADIĆ, J. (1988): Identifikacija materinskih porod jugozapadnoj časti Pannonskogo baseina, Dinari- dov i Adriatičeskogo podmorja.– Petrolgeohim, Dokladi po geohimičeskim vaprosam, 1, 155–165, Szolnok. BARIĆ, G., MARIČIĆ, M. & RADIĆ, J. (1989): Geochemical correlation of source rocks and oils in southern margins of Pannonian Basin.– 14th International meeting on organic geochemistry, Abstracts, 218, Paris. BARIĆ, G., MESIĆ, I. & JUNGWIRTH, M. (1992): Sedi- mentary and geochemical characteristics of rocks and fluids of the western part of the Drava depression.– Nafta, 43/4–5, 225–238, Zagreb. BENKOVICS, L. (1997): Etude structurale et géodynamique des Monts Buda, Mecsek et Villány (Hongrie). Thèse de doctorat (Structural and geodynamic study of the Buda Hills, Mecsek and Villány Mountains, in French with English and Hungarian summary).– Unpubl. PhD Thesis, University of Lille, 230 p. BÉRCZI, I., HÁMOR, G., JÁMBOR, Á. & SZENTGYÖR- GYI, K. (1988): Neogene sedimentation in Hungary.– In: ROYDEN, L.H. & HORVÁTH, F. (eds.): The Pannonian Basin – A Study in Basin Evolution. AAPG Memoir, 45, 57–67. BERGERAT, F. & CSONTOS, L. (1988): Brittle tectonics and paleo-stress field in the Mecsek and Villány Mts. (Hungary): correlation with the opening mechanism of the Pannonian Basin.– Acta Geol. Hung., 31/1–2, 81–100. BLAŠKOVIĆ, I. (1982): The Neogene of the Ilova river depresssion (northern Croatia).– Acta Geologica, 12/2, 23–67, Zagreb. BLAŠKOVIĆ, I., TIŠLJAR, J., DRAGIČEVIĆ, I. & VELIĆ, J. (1984): Razvoj sedimentacijskih okoliša miocenskih naslaga na zapadnim obroncima Psunja (sjeverna Hrvat- ska).– Geološki vjesnik, 37, 11–32. CLAYTON, J.L. & KONCZ, I. (1994): Petroleum geochemis- try of the Zala Basin, Hungary.– AAPG Bull., 78/1, 1–22. CSONTOS, L. (1995): Tertiary tectonic evolution of the Intra- Carpathian area.– Acta Vulcanologica, 7, 1–13. CSONTOS, L. & NAGYMAROSY, A. (1998): The Mid-Hun- garian line: a zone of repeated tectonic inversions.– Tec- tonophysics, 297, 51–71. CSONTOS, L., NAGYMAROSY, A., HORVÁTH, F. & KOVAČ, M. (1992): Tertiary evolution of the Intra-Car- pathian area: a model.– Tectonophysics, 208, 221–241. CSONTOS, L., BENKOVICS, L., BERGERAT, F., MANSY, J-L. & WÓRUM, G. (2002): Tertiary deformation history from seismic section study and fault analysis in a former European Tethyan margin (the Mecsek–Villány area, SW Hungary).– Tectonophysics, 297, 51–71. DANK, V. (1988): Petroleum geology of the Pannonian Pasin, Hungary: an overview.– In: ROYDEN, L. & HORVÁTH F. (eds.): The Pannonian Basin – A Study in Basin Evolu- tion. AAPG Memoir, 45, 319–331. DUNKL, I., ÁRKAI, P., BALOGH, K., CSONTOS, L. & NAGY, G. (1994): A hőtörténet modellezése fis- sion track adatok felhasználásával – a Bükk hegység kiemelkedéstörténete.– Földtani Közlöny, 124/1, 1–24. ĐUREKOVIĆ, M. (1995): Taloženje i rezervoarska obilježja IVA-pješčenjaka na naftnom polju Ivanić.– Unpubl. MSc Thesis, University of Zagreb, 59 p. FILJAK, R., PLETIKAPIĆ, Ž., NIKOLIĆ, D. & AKSIN, V. (1969): Geology of petroleum and natural gas from the Neogene complex and its basement in the southern part of the Pannonian basin, Yugoslavia.– Conference of the Insti- tute of Petroleum and AAPG, Brighton, 1969, 113–130. FODOR, L., CSONTOS, L., BADA, G., GYÖRFI, I. & BENKOVICS, L. (1999): Tertiary tectonic evolution of the Pannonian Basin system and neighbouring orogens: a new synthesis of palaeostress data.– In: DURAND, B., JOLIVET, L., HORVÁTH, F. & SÉRANNE, M. (eds.): The Mediterranean Basins: Tertiary extention within the Alpine Orogene. Geol. Soc. London, Spec. Publ., 156, 295–334. FRISCH, W., KUHLEMANN, J., DUNKL, I. & BRÜGEL, A. (1998): A palinspastic reconstruction and topographic evolution of the Eastern Alps during the Tertiary tectoni- cal evolution.– Tectonophysics, 297, 1–15. HÁMOR, G. (1970): The Miocene of the East Mecsek Mts (in Hungarian).– Ann. Inst. Geol. Hung., 53/1, 1–483. HÁMOR, G. (2001): Miocene palaeogeography of the Car- pathian basin: Explanatory notes to the Miocene palaeo- geographic maps of the Carpathian Basin 1:3.000,000.– Publication of the Geological Institute of Hungary, Buda- pest, 71+VII p.. HERNITZ, Z., VELIĆ, J. & BARIĆ, G. (1995a): Origin of hydrocarbons in the eastern part of the Drava Depression (eastern Croatia).– Geol. Croatica, 48/1, 87–95. HERNITZ, Z., VELIĆ, J., KALAC, K. & SAFTIĆ, B. (1995b): Miocenske taložine u širem području Beničanaca i njihovo naftnogeološko značenje (istočna Hrvatska).– In: VLAHOVIĆ, I., VELIĆ, I. & ŠPARICA, M. (eds.): First Croatian Geological Congress, Opatija, 1995, Proceed- ings, 1, 195–200. HORVÁTH, F. (1993): Towards a mechanical model for the formation of the Pannonian Basin.– Tectonophysics, 225, 333–358. HORVÁTH, F. (1995): Phases of compression during the evo- lution of the Pannonian basin and its bearing on hydro- carbon exploration.– Marine and Petroleum Geology, 12, 837–844. HORVÁTH, F. & CLOETINGH, S. (1996): Stress-induced late-stage subsidence anomalies in the Pannonian basin.– Tectonophysics, 266, 287–300. HORVÁTH, F. & ROYDEN, L.H. (1981): Mechanism for formation of the intra-Carpathian basins: A review.– Earth Evol. Sci., 1/3–4, 307–316. HORVÁTH, F. & TARI, G. (1999): IBS Pannonian Basin Project: a review of the main results and their bearings on hydrocarbon exploration.– In: DURAND, B., JOLIVET, L., HORVÁTH, F. & SÉRANNE, M. (eds.): The Mediter- 120 Geologia Croatica 56/1 121Saftić, Velić, Sztanó, Juhász & Ivković: Tertiary Subsurface Facies, Source Rocks... ranean Basins: Tertiary extension within the Alpine Oro- gen. Geol. Soc. London, Spec. Publ., 156, 195–213. HORVÁTH, F., DÖVÉNYI, P., SZALAY, Á. & ROYDEN L.H. (1988): Subsidence, thermal and maturation history of the Great Hungarian Plain.– In: ROYDEN, L. & HOR- VÁTH, F. (eds.): The Pannonian Basin – A Study in Basin Evolution. AAPG Memoir, 45, 355–375. IVKOVIĆ, Ž., MATEJ, S. & ŠKOKO, M. (2000): Seizmos- tratigrafska interpretacija gornjomiocenskih i pliocenskih sedimenata Savske depresije (Seismostratigraphic inter- pretation of Upper Miocene and Pliocene sediments of the Sava depression).– In: VLAHOVIĆ, I. & BIONDIĆ, R. (eds.): Second Croatian Geological Congress, Cavtat – Dubrovnik, 2000, Proceedings, 219–222 (in Croatian with English abstract). JUHÁSZ, Gy. (1991): Lithostratigraphical and sedimento- logical framework of the Pannonian (s.l.) sedimentary sequence in the Hungarian Plain (Alföld), Eastern Hun- gary.– Acta Geologica Hungarica, 34/1–2, 53–72. JUHÁSZ, Gy. (1994): Comparison of the sedimentary sequen- ces in Late Neogene subbasins in the Pannonian Basin, Hungary (Magyarországi neogén medencerészek pannóni- ai s.l. üledéksorának összehasonlító elemzése).– Földtani Közlöny, 124/4, 341–365. JUHÁSZ, Gy. & MAGYAR, I. (1992): A pannóniai (s.l.) litofáciesek és molluszka-biofáciesek jellemzése és kor- relációja az Alföldön (Review and correlation of the Late Neogene (Pannonian s.l.) lithofacies and mollusc biofacies in the Great Plain, eastern Hungary).– Földtani Közlöny, 122/2–4, 167–194. KŐRÖSSY, L. (1988): Hydrocarbon geology of the Zala Basin in Hungary.– General Geological Review, 23, 3–162. KŐRÖSSY, L. (1989): Hydrocarbon geology of the Dráva Basin in Hungary.– General Geological Review, 24, 3–122. KRANJEC, V. (1972): Subsurface structures in the Sava and Drava river depressions and their classification.– Bull. sci. Cons. Acad. Yougosl., (A), 17/11–12, 362–363, Zagreb. KOVAČEVIĆ, A. (1996): IVA-pješčenjaci na naftnom polju Ježevo.– Unpubl. BSc Thesis, University of Zagreb, 35 p. LARVA, O. (1996): IVA-pješčenjaci na naftno-plinskom polju Voloder.– Unpubl. BSc Thesis, Univ. of Zagreb, 39 p. LUČIĆ, D., SAFTIĆ, B., KRIZMANIĆ, K., PRELOGOVIĆ, E., BRITVIĆ, V., MESIĆ, I. & TADEJ, J. (2001): The Neogene evolution and hydrocarbon potential of the Pannonian Basin in Croatia.– Marine and Petroleum Geology, 18, 133–147. MAGYAR, I. & GEARY, D.H. (eds.) (1999): Fossils and strata of Lake Pannon, a long-lived lake from the Upper Miocene of Hungary.– Acta Geologica Hungarica, 42/1, 108 p. MAGYAR, I., GEARY, D.H. & MÜLLER, P. (1999a): Paleo- geographic evolution of the Late Miocene Lake Pannon in Central Europe.– Palaeogeography, Palaeoclimatology, Palaeoecology, 147, 151–167. MAGYAR, I., GEARY, D.H., SÜTŐ-SZENTAI, M., LAN- TOS, M. & MÜLLER, P. (1999b): Integrated biostratigra- phic, magnetostratigraphic and chronostratigraphic correla- tions of the Late Miocene Lake Pannon deposits.– Acta Geologica Hungarica, 42, 5–31. MATASOVIĆ, M. (1996): Petrološke i sedimentološke značajke pješčanih tijela – ležišta “C” naftno-plinskog polja Mramor Brdo.– Unpubl. BSc Thesis, University of Zagreb, 33 p. MÜLLER, P., GEARY, D.H. & MAGYAR, I. (1999): The endemic molluscs of the Late Miocene Lake Pannon: their origin, evolution, and family-level taxonomy.– Lethaia, 32, 47–60. PAVELIĆ, D. (2001): Tectonostratigraphic model for the North Croatian and North Bosnian sector of the Miocene Pannonian Basin System.– Basin Research, 13, 359–376. PAVELIĆ, D., MIKNIĆ, M. & SARKOTIĆ ŠLAT, M. (1998): Early to Middle Miocene facies succession in lacustrine and marine environments on the southwestern margin of the Pannonian basin system (Croatia).– Geologia Car- pathica, 49/6, 433–443. POGÁCSÁS, Gy. (1984): Seismic stratigraphic features of the Neogene Sediments in the Pannonian Basin.– Geophysi- cal transactions, 30/4, 373–410. POGÁCSÁS, Gy., LAKATOS, L., UJSZÁSZI, K., VAKA- RCS, G., VÁRKONYI, P., VÁRNAI, P. & RÉVÉSZ, I. (1988): Seismic facies, electro facies and Neogene Sequ- ences, Chronology of the Pannonian basin.– Acta Geol. Hung., 31/3–4, 175–207. POGÁCSÁS, Gy., MÜLLER, P. & MAGYAR, I. (1993): The role of seismic stratigraphy in understanding biolo- gical evolution in the Pannonian lake (SE Europe, Late Miocene).– Geologia Croatica, 46, 63–69. PRELOGOVIĆ, E., SAFTIĆ, B., KUK, V., VELIĆ, J., DRA- GAŠ, M. & LUČIĆ, D. (1998): Tectonic activity in the Croatian part of the Pannonian basin.– Tectonophysics, 297, 283–293. PUTNIKOVIĆ, A., BRITVIĆ, V., DRAGAŠ, M. & MAJSTO- ROVIĆ, J. (1989): Geološka interpretacija i usporedba geokemijskih podataka u svrhu sagledavanja međusobne povezanosti matičnih stijena i ugljikovodika u Panonskom bazenu.– DIT INA–Naftaplin, 26–27, 87–94, Zagreb. RADIĆ, J., BRITVIĆ, V. & DRAGAŠ, M. (1989): Matične stijene Panonskog bazena (Istražno područje INA–Nafta- plin).– DIT INA–Naftaplin, 26–27, 77–85, Zagreb. RÖGL, F. (1996): Stratigraphic correlation of the Paratethys Oligocene and Miocene.– Mitt. Ges. Geol. Bergbaustud. Österr., 41, 65–73. RÖGL, F. & DAXNER-HÖCK, G. (1996): Late Miocene Paratethys Correlations.– In: BERNOR, R.L., FAHL- BUSCH, V. & MITTMANN, H.-W. (eds): The Evolution of Western Eurasian Neogene Mammal Faunas. Columbia University Press, New York, 47–55. ROYDEN, L.H. & HORVÁTH, F. (eds.) (1988): The Panno- nian Basin – A Study in Basin Evolution.– AAPG Memoir 45, 394 p. SACCHI, M., HORVÁTH, F. & MAGYARI, O. (1999): Role of unconformity-bounded units in the stratigraphy of the continental record: a case study from the Late Miocene of the western Pannonian Basin, Hungary.– In: DURAND, B., JOLIVET, L., HORVÁTH F. & SÉRANNE, M. (eds.): The Mediterranean Basins: Tertiary Extension within the 122 Geologia Croatica 56/1 Alpine Orogene. Geol. Soc. London. Spec. Publ., 156, 357–390. SAFTIĆ, B. (1998): Genetska stratigrafska sekvencijska analiza u pontskim naslagama zapadnoga dijela Savske depresije.– Unpubl. PhD Thesis, University of Zagreb, 136 p. SZALAY, Á. & KONCZ, I. (1991): Genetic relation of hydro- carbons in the Hungarian part of the Pannonian Basin.– In: SPENCER, A.M. (ed.): Generation, Accummulation and Production of Europe’s Hydrocarbons. EAPG Special Publication, 1, Oxford, 317–322. SZENTGYÖRGYI, K. & JUHÁSZ, Gy. (1988): Sedimen- tological characteristics of the Neogene sequences in SW Transdanubia, Hungary.– Acta Geol. Hung., 31/3–4, 209–225. ŠĆAVNIČAR, B. (1979): Pješčenjaci miocena i pliocena sav- ske potoline.– 3. god. naučni skup Sekcije za primjenu geologije, geofizike i geokemije Znanstvenog savjeta za naftu Jugosl. akademije znanosti i umjetnosti, Novi Sad 22.–24.3.1977., Zbornik radova 2, 351–382, JAZU, Radovi ZSN, Serija A, knjiga 6/2, Zagreb (in Croatian with an English abstract). ŠIMON, J. (1973): O nekim rezultatima regionalne korelacije litostratigrafskih jedinica u jugozapadnom području Panon- skog bazena.– Nafta, 24/12, 623–630. ŠIMON, J. (1980): Prilog stratigrafiji u taložnom sustavu pješčanih rezervoara Sava-grupe naslaga mlađeg tercijara u Panonskom bazenu sjeverne Hrvatske.– Unpubl. PhD Thesis, University of Zagreb, 66 p. TARI, G. (1992): Late Neogene transpression in the North- ern Thrust zone Mecsek Mts. Hungary.– Annales of the Eötvös University Budapest, Sect. Geol., 29, 165–187. TARI, G. (1994): Alpine tectonics of the Pannonian Basin. Unpubl. PhD Thesis, Rice University, Houston, 501 p. TARI, G., DÖVÉNYI, P., DUNKL, I., HORVÁTH, F., LENKEY, L., STEFANESCU, M., SZAFIÁN, P. & TÓTH, T. (1999): Lithospheric structure of the Pannon- ian basin derived from seismic, gravity and geothermal data.– In: DURAND, B., JOLIVET, L., HORVÁTH, F. & SÉRANNE, M. (eds.): The Mediterranean Basins: Tertiary extension within the Alpine Orogen. Geol. Soc. London, Spec. Publ., 156, 215–250. TARI, V. & PAMIĆ, J. (1998): Geodynamic evolution of the northern Dinarides and the southern part of the Pannonian Basin.– Tectonophysics, 297, 269–281. TIŠLJAR, J. (1993): Sedimentary bodies and depositional models for the Miocene oil-producing areas of Ladislavci, Beničanci and Obod.– Nafta, 44/10, 531–542, Zagreb. TOMLJENOVIĆ, B. & CSONTOS, L. (2001): Neogene– Quaternary structures in the border zone between Alps, Dinarides and Pannonian Basin (Hrvatsko zagorje and Karlovac Basins, Croatia).– Int. J. Earth. Sci., 90, 560– 578. TÓTH-MAKK, Á. (2003): Szekvencia sztratigráfiai vizs- gálatok a Duna–Tisza köze középső és déli részének késő miocén rétegsorában (Sequence stratigraphy of Late Miocene deposits of the Danube–Tisza interfluve).– Unpubl. PhD thesis, ELTE Budapest, 156 p. UJSZÁSZI, K. & VAKARCS, G. (1993): Sequence strati- graphic analysis in the south Transdanubian region, Hun- gary.– Geophysical transactions, 38, 69–87. VAIL, P.R., AUDEMARD, F., BOWMAN, S.A., EISNER, P.N. & PEREZ-CRUZ, C. (1991): The stratigraphic sig- natures of tectonics, eustasy and sedimentology – an over- view.– In: EINSELE et al. (eds.): Cycles and Events in Stratigraphy, 617–659, Springer. VAKARCS, G. (1997): Sequence stratigraphy of the Cenozoic Pannonian Basins, Hungary.– Unpubl. PhD thesis, Rice University, Houston, Texas, 514 p. VAKARCS, G., VAIL, P.R., TARI, G., POGÁCSÁS, Gy., MATTICK, R.E. & SZABÓ, A. (1994): Third-order Mid- dle Miocene–Early Pliocene depositional sequences in the prograding delta complex of the Pannonian Basin.– Tec- tonophysics, 240, 81–106. VASS, D. (1999): Numeric age of the Sarmatian boundaries (Suess 1866).– Slovak Geol. Mag., 5, 227–232. VASS, D., REPČOK, I., BALOGH, K. & HALMAI, J. (1987): Revised radiometric time-scale for the Central Paratethyan Neogene.– Annals of the Hungarian Geolo- gical Institute, 70, 423–434, Budapest. VELIĆ, J. & DURN, G. (1993): Alternating lacustrine–marsh sedimentation and subareal exposure phases during Qua- ternary: Prečko, Zagreb, Croatia.– Geol. Croatica, 46/1, 71–90. VELIĆ, J. & SAFTIĆ, B. (1991): Subsurface spreading and facies characteristics of Middle Pleistocene deposits between Zaprešić and Samobor.– Geol. vjesnik, 44, 69–82. VELIĆ, J., SAFTIĆ, B. & MALVIĆ, T. (1999): Lithologic composition and stratigraphy of Quaternary sediments in the area of the “Jakuševec” waste depository (Zagreb, Northern Croatia).– Geol. Croatica, 52/2, 119–130. VELIĆ, J., TIŠLJAR, J., DRAGIČEVIĆ, I. & BLAŠKOVIĆ, I. (2000): Shoreline cross-bedded biocalcarenites (Middle Miocene) in the Podvrško–Šnjegavić area, Mt. Psunj, and their petroleum significance (Požega subdepression – Eastern Croatia).– Geol. Croatica, 53/2, 281–293. VELIĆ, J., WEISSER, M., SAFTIĆ, B., VRBANAC, B. & IVKOVIĆ, Ž. (2002): Petroleum-geological charac- teristics and exploration level of the three Neogene depo- sitional megacycles in the Croatian part of the Pannonian basin.– Nafta, 53/6–7, 239–249, Zagreb. VRBANAC, B. (2002a): Chronohorizons based on resistivity curve variations – Upper Miocene sediments of the Ivanić Grad Formation in the Sava depression, NW Croatia.– Geol. Croatica, 55/1, 11–23. VRBANAC, B. (2002b): Facies and facies architecture of the Ivanić Grad Formation (Upper Pannonian) – Sava depres- sion, NW Croatia.– Geol. Croatica, 55/1, 57–77. Manuscript received April 16, 2003. Revised manuscript accepted June 20, 2003.