Abstract

The Szeged Basin (S Hungary) occupies a relatively central position within the European Alpine–Carpathian–Dinaride orogenic belt. An ongoing controversy about the tectonic position of the study area indicates that its evolution is still not fully understood; however, several important hydrocarbon occurrences are known in the fractured basement reservoirs. The main aim of this contribution is to investigate the petrographic features and possible Alpine metamorphic conditions of volcanic/volcanoclastic and siliciclastic rocks from the Kelebia basement area. Due to the outcrop conditions and poor exposure, study samples are obtained from cores and core chips resulting from oil exploration. Based on an evaluation of petrographic (including also cathodoluminescence analysis) and microstructural features, joined with mineralogical and metamorphic data such as “illite crystallinity” and K-white mica crystallite size obtained by X-ray powder diffractometry (XRPD), a very low- to low-grade (ca. 300°C) Alpine metamorphic imprint of this portion of the basement can be proposed. Several deformation characteristics (deformation lamellae in quartz, deformation twins in dolomite, fragmented porphyroclasts, and strain shadows) were recognized in the studied samples, showing a weakly to moderately developed disjunctive foliation in the Permian rocks, as well as quartz veinlets, microcracks, and fluid inclusion planes in the Lower Triassic sandstones. Most likely, one of the Cretaceous orogenic events, namely, the “Turonian” phase (Early–Late Cretaceous nappe stacking), resulted in the prograde greenschist facies metamorphism in the study area, instead of the burial depth. We propose that the Permo-Triassic cover succession was also affected by shearing episodes accompanied by fluid migrations along the contact zone between the tectonic units. The scientific approach and dataset provided here are examples of how the application of XRPD parameters of phyllosilicates and micropetrographic observations can help to understand the evolution of an orogen and improve knowledge about the basement structure.

1. Introduction

The European Alpine–Carpathian–Dinaride (ACD) orogenic belt was the focus of several relevant projects that resulted in advanced concepts for the tectonometamorphic and thermal evolution of the region (e.g., [110]). Furthermore, correlational studies led to improved versions of the tectonic map of the ACD orogenic system [1117]. In this overall compressional setting, a large rift-related extensional basin, namely, the Pannonian Basin in eastern Central Europe, has developed from the Oligocene [14]. As an integral part of the ACD orogenic system, the back-arc Pannonian Basin is superimposed on two distinct tectonic megaunits or composite terranes with complex internal structures (e.g., [3, 7, 11, 1417]): (i) the so-called ALCAPA (Alpine–West Carpathian–Pannonian) and (ii) the Tisza–Dacia (Figure 1(a)).

The study area, corresponding to the Codru Nappe System (Figure 1(b)), occupies a relatively central position within the ACD orogenic system; however, its crystalline basement and Permo-Mesozoic cover rocks are largely covered by thick Cenozoic sedimentary sequences (up to ~7000 m, average thickness: 2000–3000 m) (Figure 1(a)). Isolated basement outcrops are known in the surrounding small hills and mountains (Mecsek and Villány Hills in Hungary, the Slavonian Mts in Croatia, and the northern Apuseni Mts in Romania) (e.g., [5, 6, 16, 17]). Therefore, local to subregional correlations between the main tectonic units are difficult and are generally based on distinct subsurface maps [12, 14, 17], necessarily having uncertainties, assumptions, and tentative interpretations. An ongoing controversy about the tectonic position of the study area (Figure 1(b)) indicates that the Alpine evolution of the central part of the ACD orogenic system is still not fully understood. One reason for this uncertainty is the fact that the pre-Cenozoic basement areas, especially the siliciclastic Permian and Triassic cover formations, are still poorly studied in terms of their timing and petrological characterization. The lack of comprehensive lithostratigraphic and diagenetic/metamorphic control, consequently the scarcity of basement rock samples from the exploration wells, is also one of the important limiting factors.

During the periodic hydrocarbon (HC) exploration projects, it was beyond the scope of the reports to discuss the internal structure of the basement from a pure geological point of view. Standard thin sections with cover glass were prepared from selected basement core samples; however, detailed mineralogical and petrographic data are not available on these rocks. The principal aim of this contribution is to provide a mineralogical and petrological characterization of the Permo-Triassic noncarbonate cover succession in the Kelebia basement area, Szeged Basin. Regarding this study area, the key results together with a compilation of literature data can lead to a better understanding of the tectonic evolution and Alpine structure of the pre-Cenozoic basement and allow for further constraints on the ACD orogenic system.

2. Materials and Methods

2.1. Geological Setting

The Szeged Basin is a subbasin on the southern margin of the Pannonian Basin in Hungary (Figure 1(a)). Within it, the pre-Cenozoic basement belongs structurally to the Tisza megaunit (Tisia Composite Terrane, as referred elsewhere [11, 18]), which traditionally comprises four main nappe systems of Alpine origin formed in the Albian to Turonian (from bottom to top): Mecsek, Bihor, Codru, and Biharia [3, 1416, 19]. In the Hungarian part of the Tisza megaunit, the Mecsek, Bihor (Villány–Bihor), and Codru (Békés–Codru) Nappe Systems are only found [1416, 1921]. The Mecsek and Bihor tectonic units generally consist of a Variscan metamorphic basement intruded by Carboniferous granitoids, as well as Permo-Mesozoic volcanic and sedimentary cover, whereas the Codru Nappe System predominantly consists of nonmetamorphic Permo-Mesozoic rocks of variable thickness which are locally preserved. In the Biharia Nappe System, the pre-Alpine terranes, composed of metamorphic sequences, are occasionally covered by Alpine metamorphosed Permo-Triassic sequences (Figure S1) [8, 16, 1921]. It is important to note, however, that the highest Biharia Nappe System is taken as a constituent of the Dacia megaunit according to a different interpretation as discussed elsewhere [8, 12, 13].

On the other hand, Hungary has a long history of exploration and production of HCs [2224]. Among the principal producing areas, the southern part of the Pannonian Basin has an excellent HC potential. Several commercially important hydrocarbon occurrences are also known in the Szeged Basin, such as the Ásotthalom North (Ás–É), Dorozsma (Do), Kelebia North (Kel), Mórahalom (Móra), and the Szeged petroleum and natural gas fields (Figure 2). The most productive conventional HC reservoirs reside in the Neogene sediments (e.g., Miocene conglomerate and sandstone beds); nevertheless, important ones exist in the fractured crystalline (Kel and Do fields) and pre-Cenozoic volcanic (Kel field, fractured Permian rhyolite) and sedimentary (Móra and Szeged fields, brecciated Triassic dolomite) basement areas, about 1000–3000 m below the present surface [2125]. Fortunately, during the crude oil exploration projects (from 1968 to 1985), several wells were drilled with intermittent usage of coring equipment in the Kelebia basement area (wells Kel and Ás–É). Some of them penetrated the basement top, which essentially consists of Variscan metamorphic rocks (e.g., gneiss and mica schist). Permian cover succession is also discovered locally [24].

The southern part of the basement of the Szeged Basin (Figures 1 and 2) contains blocks made up of medium-grade Variscan metamorphic rocks, Permian to Lower Triassic siliciclastic and volcanosedimentary formations (Korpád Sandstone, Jakabhegy Sandstone, and Gyűrűfű Rhyolite Formations, respectively), and shallow marine mudrocks and carbonates from the Middle Triassic (Szeged Dolomite Formation). Because of the tectonically induced uplift and erosion from the Late Cretaceous, the post-Triassic Mesozoic deposits are totally absent in the study area [16, 21, 24, 25]. The topographic undulations of the basement were covered by marine sediments during the Miocene (i.e., 20–10 Ma). Subsequently, the paleo-Danube and paleo-Tisza Rivers progressively filled the Pannonian Basin with clastic material, creating a thick Neogene nonmarine succession [4, 7, 2123].

The alluvial Korpád Sandstone Formation (Lower Permian) consists of red siliciclastic sedimentary rocks (dominated by sandstones and mudrocks) [16, 24]. Permian (, using the zircon U–Pb LA–SF–ICP–MS method) Gyűrűfű Rhyolite samples from the Kelebia area contain oriented, sericitized pumices, and broken subhedral (micro)phenocrysts (15–20 vol% quartz, 10–15 vol% altered plagioclase, and ~5 vol% altered biotite). Based on the primary composition and microtextural characteristics, these felsic volcanic rocks are rhyodacitic ignimbrites [26]. The Lower Triassic Jakabhegy Sandstone Formation, represented by gray, red, and lilac fluvial to shallow marine sandstones, covers the crystalline basement and older Permian rocks with unconformity [16, 21, 24]. The Late Anisian to Early Ladinian brecciated dark gray dolomite sequence, corresponding to the Szeged Dolomite Formation, shows evidence for multistage dolomitization [21].

2.2. Sampling

For mineralogical and petrographic analysis of the Permian terrigenous sedimentary and volcanogenic/volcaniclastic rocks, the collections of archive thin sections at the host institution were studied from HC exploration wells near the village of Kelebia (Kel field; Permian red beds underlain by the ignimbrites, Permian ignimbrites, and volcaniclastic sedimentary rocks; Gyűrűfű Rhyolite Formation: wells Kelebia–7, 11, 12, 14, 22, and 23; Korpád Sandstone Formation: wells Kelebia–12 and 14) at depths of ~800–900 m below the surface (Figure 2). Additionally, corresponding to the scarcity of the basement cores, a total of 3 small core chips from the wells Kelebia–7 (877.5–879.5 m and 900.0–903.0 m; Gyűrűfű Rhyolite) and Kelebia–14 (905.0–910.0 m; Permian red beds underlain by the ignimbrites) were newly sampled for this study at the repository of the host institution (Table 1).

For Lower Triassic siliciclastic sedimentary rock characterization, core samples from the exploration well Mórahalom–1 (well Móra–1, Móra field) were examined (Figure 2). The available core section #17 (Jakabhegy Sandstone Formation) from the depth interval between 2306 and 2312 m below the surface was taken for detailed petrographic studies, including diagenetic and microtectonic observations. Additionally, the MOL Group (an integrated oil and gas company, Hungary) provided a covered thin-section collection from wells near the village of Ásotthalom (Ás–É field) at depths of ~2000–2100 m below the surface. The available standard thin sections belong to the Lower Miocene marine conglomerate sequence, including reworked pebbles and cobbles of Jakabhegy Sandstones (Table 1).

2.3. Methodology of Mineralogical and Petrographic Analyses

Determination of the bulk mineralogical composition and characterization of the separated <2 μm grain size fraction of the selected Permian drill core chip samples () was made by X-ray powder diffractometry (XRPD). The representative samples were measured by a Rigaku Ultima IV X-ray diffractometer using the Bragg-Brentano geometry (instrumental parameters: CuKα radiation, graphite monochromator, proportional counter, divergence, and detector slits of 2/3°). Macroscopically homogeneous rock chips were grounded and pulverized in an agate mortar. Random powder mounts were made using ~0.04 g rock powder on a Si single crystal sample holder to determine the bulk mineralogical composition. The qualitative and semiquantitative evaluation of the XRPD spectra was performed by the Rigaku PDXL 1.8 software (ICDD database, PDF2010) and estimated based on the reference intensity ratio method, respectively.

Grain size separation for clay fraction (<2 μm) analysis was achieved by repeated ultrasonic deflocculation and gravitational settling using Stokes’ law. Highly oriented XRPD slides with 3 mg/cm2 density were prepared by repeated sedimentation of the separated fraction on a standard glass sample holder. Both air-dried and ethylene glycol-solvated preparations were scanned at 45 kV/35 mA, from 2 to 50°2θ with goniometer step rate 1°/min and step width 0.1°. Subsequently, the mounts were treated at 350 and 550°C and measured immediately after storage in a desiccator. Semiquantitative composition of the clay fraction was estimated following Moore and Reynolds’ recommendation [27].

For the determination of illite “Crystallinity Index” and calculation of crystallite size, a triplicate scan of the same slide was run at 40 kV/30 mA, from 3 to 14°2θ with goniometer step rate 0.6°/min and step width 0.01°. The determination of “crystallinity” was made using the Crystallinity Index Standards (CIS) of Warr and Rice [28] and the Kübler Index () calculated after Warr and Mählmann [29] as discussed elsewhere [30]. Additionally, the illite crystallite size was determined by a routine single peak method (on the 001 reflection of the K-white mica) using the Scherrer equation embedded in the PDXL 1.8 software. For differentiation between K-, Na-, and Ca-rich white micas (i.e., illite, paragonite, and margarite), the XRPD (00,10) basal reflections around 2 Å were checked using highly oriented slides as discussed elsewhere [31, 32]. For a better resolution and precise determination of the position of the 00,10 K-white mica peaks, triplicate analyses were performed at 50 kV/40 mA, from 40 to 50°2θ with a goniometer scan speed of 0.5°/min and a step width of 0.02°.

Petrographic studies were conducted on small hand specimens () and standard thin sections () that were studied by polarized light microscopy using an Olympus BX41 microscope. Framework grains of sandstones were classified by common micropetrographic methods [33, 34]. The microstructural analysis of the samples is based on Passchier and Trouw [35]. However, notice that there is no information on the definite orientation of the core samples studied. Cathodoluminescence (CL) microscopy of selected Gyűrűfű samples was carried out using a Reliotron VII type cold CL device operating at a 7 kV accelerating voltage and 0.7 mA beam current, mounted on an Olympus BX43 microscope. Photomicrographs were made using an Olympus DP73 digital microscope camera. Conventional mineral abbreviations have been used during the whole paper [36].

3. Results and Discussion

3.1. XRPD Data: Mineralogical Composition and Characteristics of the <2 μm Fraction

Based on semiquantitative XRPD analyses, the bulk volcaniclastic samples (Gyűrűfű Rhyolite) have a rather uniform composition (expressed in mass%). Quartz (40–50%) and white mica (illite±muscovite, 30–50%) are the major components with some K-feldspar (~5%) and trace amount of plagioclase (albite), chlorite, dolomite, and hematite. As far as the red bed sample studied is concerned, it has a bulk composition of high amounts of white mica (70–80%) and quartz (20–30%). Some chlorite, plagioclase, and hematite were also detected, but K-feldspar and dolomite phases are missing relative to the bulk volcaniclastic samples (Figure 3(a)).

The mineralogical composition of the <2 μm fraction of both studied lithologies is rather similar: white mica is the predominant phase (>90%) with minor chlorite and a trace amount of quartz±hematite (Figure 3(b)). Regarding the clay fraction of the red bed sample, the position of the peak 00,10 shows a spacing (Figure 4). This value together with the lack of shoulders on any side of this reflection suggests a near-theoretical “muscovite” (K-white mica) mineralogical composition without any Ca or Na substitutions [37]. After the saturation of the samples with ethylene glycol, the XRPD patterns prove that the K-white mica does not contain a detectable swelling phase (i.e., smectitic interstratification). The Esquevin index value [38] is ~0.26 which indicates a significant Fe substitution and, consequently, phengitic composition for the K-white mica of the studied red bed sample. The determined value refers (Figure 5(a)), close to the anchizone–epizone boundary () [3942]. The main apparent crystallite size of K-white mica corresponds to which is close to the anchizone–epizone transition (~520 Å) [43]. Consequently, XRPD data suggest a very low- to low-grade metamorphism of the Permian sequence at ~300°C peak metamorphic temperature (Figure 5(b)).

3.2. Petrographic Fingerprints of Low-Grade Metamorphic Processes
3.2.1. Permian Siliciclastic Rocks

The fine-grained matrix of the Korpád samples overlain by ignimbrites has a weakly to moderately developed foliation with oriented sericite-rich (white mica) bands. Broken crystals of volcanic quartz, plagioclase, and muscovite occur as silt- to sand-sized detrital grains (Figures 6(a)6(d)). Occasionally, tiny rhombohedral carbonate crystals with a somewhat saddle-like shape (likely dolomite, first generation) appear in the matrix. In the overlying Gyűrűfű ignimbrites, the presence of dolomite was confirmed using the XRPD method (Figure 3(a)). The larger grains of dolomite (~200–400 μm) are strongly undulatory, show curved crystal boundaries, and have fibrous strain shadows subparallel to the foliation (Figures 6(c) and 6(d)). Around the core objects, displacement-controlled carbonate (dolomite, second generation)+hematite fringes were first developed, following quartz fringes. In a ductile deforming rock, such strain shadows were commonly formed by precipitation from the solution in the case of low-T deformation and high fluid pressure [35].

The studied red bed samples, overlying the Gyűrűfű ignimbrites, are composed of relatively well-foliated pale green to red fine-grained rocks, showing a slaty appearance. On the thin-section scale, spaced cleavage is defined by fine-grained micaceous material with detrital quartz and muscovite stacks in microlithons (Figures 6(e) and 6(f)).

3.2.2. Permian Ignimbrites

In the Kelebia area, the samples of the Gyűrűfű Rhyolite Formation contain ca. 30–50 vol% crystal fragments with an assemblage of quartz+K-feldspar+plagioclase+biotite and ca. 5–10% altered pumice clasts (sericitized) that float in a devitrified and strongly recrystallized groundmass. Moreover, there are some accessory minerals such as apatite, zircon, and opaque grains. Feldspars are partially replaced by dolomite and/or reflect features of sericitic alteration. Mafic minerals (biotite, rare pyroxene) are also moderately to strongly altered (Figure 7).

On the thin-section scale, several deformation features were recognized in the studied ignimbrite samples, showing a weakly to moderately developed disjunctive foliation (in the sense of Passchier and Trouw [35]). Quartz phenocrysts are frequently fractured, fragmented, and variably deformed, so evidence of the original crystal shape is generally absent. Crystal fragments are typically angular in shape; however, the resorbed crystal outlines are rarely recognizable. Broken quartz crystal fragments also show some evidence of rotation and displacement. Quartz grains mostly show undulatory extinction and deformation lamellae in cross-polarized light and exhibit a low-intensity blue (core region) to relatively high-intensity violet-red (rim) reverse zoning luminescence pattern in CL (Figure 7 and Figure S2). Furthermore, domino-type and mosaic fragmented quartz clasts (porphyroclasts in the sense of Passchier and Trouw [35] as discussed below) are common (Figure 8 and Figure S2).

Understanding the origin of broken and fragmented crystals in volcanic deposits can often be challenging, as many crystal-rich ignimbrites contain extensively fragmented crystals [35, 44]. According to van Zalinge et al. [44], primary fragmentation caused by rapid decompression is associated with the eruption, whereas secondary fragmentation is related to compaction and welding. Generally, the primary process is characterized by an extensive crystal breakage and expansion of individual crystals (~5 vol%), but their original shapes are preserved. The secondary one is generally associated with enhanced fragmentation localized at points of contact between adjacent fragments, resulting in loss of volume between them [44]. Secondary fragmentation increases the crystal aspect ratio, therefore forming textures that resemble mantled porphyroclasts (i.e., crystals with a rim, which is of the same composition as the clast) in shear zones [35, 4446]. Based on our micropetrographic observations (e.g., presence of deformation lamellae, domino-type and mosaic fragmented crystals with rotated chips, absence of welded pumice clasts, and insignificant amounts of mantled clasts), we suggest that the abovementioned features of brittle fracturing together with the ductile (intracrystalline) deformation of quartz (Figure 8 and Figure S2) reflect signs of tectonic processes (e.g., pure to simple shear) in a variously deformed ignimbrite, instead of compaction and welding.

In the moderately well-foliated samples (e.g., well Kel–7 cores #7 and #8, well Kel–12 core #1, and well Kel–23 core #3), the cleavage domains are rich in white mica (i.e., sericite and/or muscovite) subparallel to the trend of the domains. In microlithons, fibrous or elongated strain shadows are common on both sides of the core object (e.g., white mica strain fringe around apatite and quartz strain shadow around dolomite, respectively). Such strain shadows are composed of white mica around the quartz porphyroclasts, showing a gradual transition between the fringes and the recrystallized matrix (Figure 8 and Figure S2). Mantled quartz clasts also occur locally (sample Kel–23). It is important to note that mechanical twins can be distinguished in replacement dolomite (first generation). Additionally, between the quartz porphyroclast fragments, some undeformed dolomite rhombs (second generation) appear in the recrystallized matrix (Figure 8), reflecting postkinematic mineralizing fluids.

3.2.3. Triassic Siliciclastic Rocks

The Triassic sandstone samples from the Móra field are moderately to well-sorted subarkose to quartz arenite (>95% quartz). They are pink or pale red through hematite pigmentation; additionally, fossils are absent. The samples are composed of rounded mono- and polycrystalline quartz and lesser amounts of K-feldspar (microcline) and volcanic grains. Accessories are zircon and tourmaline. Quartz cement as syntaxial overgrowths is common. Deformation lamellae occur in both the original quartz framework grain and its overgrowth. Microcrystalline silica and carbonate replacements are also typical. Locally, wide (>1 μm) deformation twins were developed in the replacement carbonate crystals. Additionally, the occurrence of tiny quartz veinlets and transgranular microcracks marked by fluid inclusion planes (FIPs) is a characteristic feature of the samples (Figure 9 and Figure S3). These signs reflect that fluids have to be involved in the cementation of brittle dilatational microfaults.

In the Ásotthalom area, belonging to the Ás–É field, the Lower Miocene conglomerate samples are composed of acidic volcanic, metamorphic (quartz-sericite/muscovite schist), and sedimentary rock clasts derived from the pre-Neogene basement (Figure S4). Quartzite and quartz arenite clasts are the main components, showing strong similarities with the Triassic sandstones in the study area (Figure 10). The appearance of deformation lamellae in quartz overgrowth cement together with quartz veinlets and FIPs in quartz arenite clasts clearly indicates that the time of ductile and brittle deformation and contemporaneous fluid migration preceded the Neogene uplift and erosion (i.e., partial subaqueous exposure).

3.3. Diagenesis versus Metamorphism: Evaluation of the Temperature Conditions

To establish the T (and P) conditions reached by the studied rocks of the Permo-Triassic noncarbonate cover succession in the Kelebia basement area, mainly temperature-dependent characteristics are evaluated and interpreted (Table 2). The slaty appearance, together with the disjunctive foliation, is one of the most striking features of the Permian samples subjected to increased T and P, corresponding to very low- to low-grade metamorphic conditions at temperatures greater than 200–250°C [35, 4648]. Based on the XRPD data ( value and the main apparent crystallite size of the K-white mica), corresponding to a progressive trend of phyllosilicate crystal chemical changes [3943, 47, 48], a peak metamorphic temperature of ~300°C (anchizone–epizone boundary) can be proposed.

Additionally, deformation lamellae in tectonically deformed quartz and twin morphology in carbonate minerals, especially in calcite, can be used as a crude geothermometer [35, 45, 46, 49]. The plasticity of quartz is generally restricted to temperatures greater than 200°C, and the presence of deformation lamellae could be taken as a qualitative indicator of relatively high stress levels. At low-grade conditions (300–400°C) undulose extinction of quartz and deformation lamellae in a subparallel, nonplanar arrangements are characteristic structures [35, 46, 49, 50]. On the other hand, plasticity in the form of twinning can appear even at room temperature in calcite (type I). Wider twins that can be optically resolved () dominate above 200°C up to 300°C; however, twinning apparently does not develop below 300°C in dolomite [35, 45, 49]. Mechanical twinning in dolomite requires higher critical resolved shear stress, and so it is activated at higher T compared to calcite [45]. All of these observations coincide with the high anchizone–epizone values mentioned above.

Another important property of quartz is the variability of its CL characteristics in dependence on the specific origin (e.g., different grades of metamorphism) [51, 52]. A typical brown CL of quartz is a characteristic feature of low-T metamorphic rocks, whereas quartz derived from metamorphic conditions with temperatures above 500°C exclusively exhibits visible blue CL that is possibly related to recrystallization processes [52]. The CL study of the Permian ignimbrites from the Kelebia area indicated that homogenous blue CL was not observed, so such high-T conditions can be excluded. Quartz crystals with low-intensity blue CL cores overgrown by a relatively high-intensity violet-red CL rim (reverse zoning; Figure 7) can exhibit the original magmatic features. According to Vernon [46], blue to violet CL colors are associated with enrichment in Ti. Based on evidence from quartz compositional zoning and titanium-in-quartz geothermometry, bright-CL rims are generally Ti-richer than dark-CL cores in ignimbrite units, reflecting a magmatic recharge event as a potential eruption trigger [53, 54]. This phenomenon is consistent with the genetic interpretation of the Gyűrűfű Rhyolite Formation in the study area [26].

It is noteworthy that the determination of metamorphic conditions in very low- to low-grade metamorphic successions of (silici)clastic rocks is not an easy task [42, 48]. However, specific criteria, such as crystallite size values and crystal chemical parameters of the sheet silicates or other temperature-dependent petrographic characteristics summarized in Table 2, can be considered as evidence for the progressive metamorphism of the studied units. Taking notice of the different positions in the geological literature on the existence or absence of prograde metamorphic post-Variscan cover formations within the Codru Nappe System [8, 16, 20, 26, 55], the mineralogical and petrological data presented here clearly reflect that the Permo-Triassic noncarbonate cover succession in the Kelebia basement area suffered a high anchizonal to epizonal (~greenschist facies) metamorphic alteration after their deposition.

The abovementioned results are consistent with the lithological interpretation of wells in northern Vojvodina and in the southern part of the Neogene Pannonian Basin, Serbia, close to the study area. The wells at Palić (~27 km SE of Kelebia) penetrated a low-grade Permian metamorphic sequence composed of schistose, reddish brown to grayish green metasiltstones and metasandstones, alternating with rhyolitic to rhyodacitic volcanic and volcaniclastic rocks [56].

3.4. Time–Temperature Path of the Kelebia Basement Area: Tectonothermal Evolution

Based on detailed petrographic investigations and interpretation, a tentative time–temperature (depth, qualitative) path is summarized in Figure 11. Clearly, the observed Permian and Triassic rocks were deposited in continental environments. Based on well-rounded grains, good sorting, and high textural maturity together with a compositional maturity (i.e., quartz arenites), the studied Triassic sandstone samples, at least partially, are interpreted to be of aeolian origin [34]. A published burial model for Korpád and Jakabhegy Sandstones is not available; however, we suggest that the initial subsidence during the Middle to Late Triassic was relatively rapid (mesogenesis I). The formation of replacive saddle dolomite (Figure 6) is generally a late diagenetic event and requires minimum temperatures of 60 to 80°C (typically 90 to 160°C) [21, 58]. Additionally, the precipitation of syntaxial quartz overgrowths in most sedimentary basins occurs during deep burial (>2.5 km) diagenesis at elevated temperatures (typically 90 to 130°C) [59]. In the study area, Lower Triassic sandstones are overlain by marine sequences, including Middle Triassic nonmetamorphic carbonates (Szeged Dolomite Formation) with a maximum burial depth of approximately 5–6 km (<200°C, corresponding to a normal geothermal gradient of 35°C/km) before the Late Cretaceous to Paleogene inversion and exhumation [21]. Therefore, the maximum burial depth for the oldest studied units (Permian siliciclastics and ignimbrites) is proposed to be at 6–7 km (Figure 11). The burial model proposed for the Upper Permian deposits (Mecsek Hills, southern Hungary, ~160 km W of Kelebia; Figure 1(a)) agrees with a burial that reaches 200–250°C in the Early Jurassic [60].

On the other hand, the appearance of the deformational microstructures (e.g., disjunctive foliation, tectonically deformed quartz, and twin morphology of carbonate minerals; Figures 8 and 10) unequivocally indicates that the Permo-Triassic sequence was subject to tectonic stress before the Early Miocene uplift and erosion which later stage corresponds to an exhumed position close to the coeval surface (Figure 11). Within the ACD system, undoubtedly, the most significant orogenic events are related to the thrusting and nappe stacking during the Cretaceous [412]. Very likely, these tectonic processes favored the prograde reaction series of phyllosilicates (~300°C temperature for the peak metamorphism), instead of burial depth.

Unfortunately, the Alpine orogenic phases that shaped the pre-Neogene basement in the Kelebia basement area are poorly studied in terms of kinematics and timing of deformation. Nevertheless, five tectonic events were proved in the nearby Apuseni Mountains (Figure 1(a)): (i) Late Jurassic obduction of the Transylvanian ophiolites (D0), (ii) late Early Cretaceous (~130 Ma) nappe stacking (D1, “Austrian” phase), (iii) Early–Late Cretaceous (110–90 Ma) nappe stacking (D2, “Intra-Turonian” or “Turonian” phase), (iv) Late Cretaceous extension and differential exhumation (D3), and (v) latest Late Cretaceous–Early Paleogene compression with brittle conditions (D4, “Laramian” phase) [6, 8]. It is important to note that no evidence of the D0 and D1 phases can be found in the Tisza megaunit [6, 8, 17]. During the D1 phase, the Biharia Nappe System, being a part of the Dacia megaunit [12, 13], was buried to a minimum of 8 km [6]. Therefore, despite the metamorphosed Permo-Triassic cover, the Kelebia basement area could not belong to this structural unit.

An Albian to Turonian age is generally proposed for the main deformation phase (D2) that forms the current geometry of the Alpine nappe stack [412]. The D2 phase led to a pervasive retrograde greenschist facies overprint in the structurally lower parts of the nappe pile. Additionally, the Codru Nappe System was also affected by the Turonian metamorphic overprint [6, 8]. Most likely, corresponding to the first significant Alpine orogenic event in the Tisza megaunit, the D2 phase resulted in the metamorphism of the prograde greenschist facies in the study area (Figure 11). This agrees with published mineralogical and geochronological data from wells in SE Hungary (basement high near Szeged and its surroundings), where the Alpine prograde metamorphism accompanied by mylonitization is attributed to Late Cretaceous Ar–Ar ages (mainly 96–82 Ma and somewhat younger) [55, 61, 62].

Furthermore, in the Apuseni Mountains, an approximately 30 km wide and more than 200 km long anastomosing Alpine crustal shear zone (Highiș–Biharia shear zone, HBSZ) was recognized along the boundary between the Tisza and Dacia megaunits that separates the Codru and Biharia Nappe Systems [8, 63]. Within the HBSZ, a minimum of 20 km of sinistral displacement was proposed, accompanied by phyllonite formation and an influx of fluid during greenschist to subgreenschist conditions. The lower part of the Biharia Nappe System was affected by repeated shearing events, and fluids were channeled by Alpine tectonic features [8, 6365]. Additionally, it was suggested that the interval of the fractured crystalline basement rocks near Dorozsma (Do field, Figure 2) could be the western continuation of the HBSZ [25]. Based on the observed microtectonic indicators (e.g., fluid-related fringe structures, veinlets, FIPs, deformation lamellae in quartz, and fragmented porphyroclasts), we propose that the study area could belong to the abovementioned major shear zone, having both brittle and ductile segments. Consequently, the upper part of the Codru Nappe System, including the Permo-Triassic cover succession, was also affected by shearing episodes along the contact zone between the megaunits. It seems to be reasonable that the HBSZ could be dissected during the Miocene extension of the Pannonian Basin, so different deeply buried segments of the shear zone could be expanded and moved away from each other at some distance.

Cretaceous tectonic elements could be easily reactivated during the Neogene rift-related extension, establishing structural control on fluid migration into fractured basement reservoirs. After the Middle Miocene, an intensive subsidence of the basement (mesogenesis II) resulted in the deposition of a thick (up to 4 km) sedimentary succession in the study area. The generation of HCs and the charging of the reservoirs occurred from the Late Miocene (or Pliocene) to the present-day [21, 25].

4. Conclusions

This case study highlights the major effects of diagenetic to metamorphic overprint on the mineralogy and micropetrographic features of three Permo-Triassic lithologies (Permian siliciclastics and ignimbrites, as well as Triassic sandstones) located in the pre-Cenozoic basement, Codru Nappe System (Szeged Basin, Hungary). The following conclusions can be drawn about the studied basement rocks: (1)Mineralogical and petrological data suggest that the Permo-Triassic noncarbonate cover succession in the Kelebia basement area suffered a high anchizonal to epizonal metamorphic alteration after their deposition. Additionally, several deformation features were recognized in the studied samples, showing a weakly to moderately developed disjunctive foliation in the Permian rocks(2)The occurrence of deformation lamellae in quartz overgrowth cement together with quartz veinlets and FIPs in Lower Triassic quartz arenite clasts indicates that the time of ductile and brittle deformation and contemporaneous fluid migration was before the Early Miocene partial surface exposure(3)Most likely, Cretaceous orogenic events favored the prograde reaction series of phyllosilicates (~300°C temperature for the peak metamorphism), instead of the depth of burial. We propose that the “Turonian” phase (Early–Late Cretaceous nappe stacking) resulted in the prograde greenschist facies metamorphism in the study area. Furthermore, the Permo-Triassic cover succession of the Codru Nappe System was also affected by shearing episodes along the contact zone between the Tisza and Dacia megaunits

The results provided here are an example of how knowledge of clay mineralogical and complex petrographic investigations can improve the interpretation of the tectonometamorphic evolution of a deeply buried orogen.

Data Availability

The mineralogical and petrographic data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This work was financially supported by the University of Szeged Open Access Fund (4635), the Bolyai Research Scholarship of the Hungarian Academy of Sciences (BO/266/18), and the National Research, Development, and Innovation Fund (projects no. K 108375 and K 131690). The authors would like to thank Tivadar M. Tóth (head of the Department of Mineralogy, Geochemistry and Petrology, University of Szeged) and Balázs Kiss (MOL Plc, Algyő) for providing the core material and thin-section collections that were essential for the research.

Supplementary Materials

Figure S1: generalized tectonic stratigraphy of the northern Apuseni Mountains [20]. Abbreviations: P: Permian; Tr: Triassic; J: Jurassic; K: Cretaceous. Figure S2: thin-section photomicrographs of the Permian ignimbrites, Kelebia area: alteration and deformation textures (a–e: well Kel–7 core #8, f: well Kel–23 core #3). (a) Deformation lamellae in a resorbed and fractured quartz phenocryst. (b, c) Domino-type and mosaic fragmented quartz crystals. Notice rotated quartz chips that flake off from the original fractured grains. (d, e) White mica strain shadow around an apatite grain. (f) Disjunctive foliation in a deformed ignimbrite sample (the trend of the cleavage domains is indicated by the yellow line). Notice mantled quartz clasts with strain shadows associated with strain caps. Abbreviations: PPL: plane-polarized light; XPL: cross-polarized light; Ap: apatite; Bt: biotite; Dol: dolomite; Ms: muscovite; Qz: quartz. Figure S3: thin-section photomicrographs of the Lower Triassic sandstones (well Móra–1 core #17). (a–d) Framework grains with quartz overgrowths. Quartz overgrowth cement postdated by silica (b) and carbonate (c) replacements. Note: thick (>1 μm) deformation twins in carbonate. The appearance of tiny quartz veinlets (d) and transgranular microcracks is a characteristic feature. Abbreviations: XPL: cross-polarized light; Cb: carbonate; Lv: volcanic lithic fragment; Qm: monocrystalline quartz; Qp: polycrystalline quartz; Qz: quartz. Figure S4: thin-section photomicrographs of the representative clast types of the Lower Miocene conglomerate (well Ás–É–4 core #3). Abbreviations: PPL: plane-polarized light; XPL: cross-polarized light; FIP: fluid inclusion plain; Lm: metamorphic lithic fragment; Lm(s): metasedimentary lithic fragment (quartzite); Ls: sedimentary lithic fragment (dolomite); Lv: volcanic lithic fragment; Qm: monocrystalline quartz; Qp: polycrystalline quartz; Bt: biotite; Ms: muscovite; Ser: sericite; Qz: quartz. (Supplementary Materials)