From marginal outcrops to basin interior: a new perspective on the sedimentary evolution of the eastern Pannonian Basin

Sedimentary successions exposed at basin margins as a result of late-stage inversion, uplift and erosion usually represent only a limited portion of the entire basin fill; thus, they are highly incomplete records of basin evolution. Small satellite basins, however, might have the potential of recording more complete histories. The late Miocene sedimentary history of the Șimleu Basin, a north-eastern satellite of the vast Pannonian Basin, was investigated through the study of large outcrops and correlative well-logs. A full transgressive–regressive cycle is reconstructed, which formed within a ca. 1 million-year time frame (10.6–9.6 Ma). The transgressive phase is represented by coarse-grained deltas overlain by deep-water lacustrine marls. Onset of the regressive phase is indicated by sandy turbidite lobes and channels, followed by slope shales, and topped by stacked deltaic lobes and fluvial deposits. The deep- to shallow-water sedimentary facies are similar to those deposited in the central, deep part of the Pannonian Basin. The Șimleu Basin is thus a close and almost complete outcrop analogue of the Pannonian Basin’s lacustrine sedimentary record known mainly from subsurface data, such as well-logs, cores and seismic sections from the basin interior. This study demonstrates that deposits of small satellite basins may reflect the whole sequence of processes that shaped the major basin, although at a smaller spatial and temporal scale.


Introduction
In the late Miocene, the paleogeography and climate of Central Europe was shaped by giant Lake Pannon covering the area between the Alps, Carpathians, and Dinarides (e.g. Kázmér 1990;Bruch et al. 2006). This Caspian-type, endorheic, brackish, and very deep lake and its highly endemic biota emerged at the beginning of the late Miocene when the evolution of the orogenic belts isolated the Pannonian Basin System from the Paratethys Sea (ter Borgh et al. 2013). In the nineteenth and early twentieth centuries, efforts to understand the paleogeographic evolution of the lake were focused on surface outcrops located in zones where the lacustrine layers were uplifted and eroded due to basin inversion (Ruszkiczay-Rüdiger et al. 2020 and references therein). Beginning from the 1930s, hydrocarbon exploration boreholes revealed that the thickness of the lacustrine and associated deltaic and fluvial sequence attains 5 km in the deepest subbasins. Correlation between the thick basinal successions and the exposed, thin marginal records, however, remained ambiguous until the basin fill geometry was imaged by seismic stratigraphy in the 1980s (Pogácsás 1984;Trkulja and Kirin 1984;Marton 1985;Pogácsás et al. 1988), providing a basis for the first dynamic depositional models (Révész 1980;Bérczi and Phillips 1985;Pogácsás and Révész 1987;Mattick et al. 1988).
Seismic data revealed that the basin of Lake Pannon was gradually filled with sediments transported by large fluvial systems mostly from the surrounding Alps and Carpathians. The western and central part of the basin was filled by the paleo-Danube system, deriving sediments from the Alps and from the Western Carpathians, whereas sediments into the eastern part of the basin were transported from the Eastern Carpathians and from the Apuseni Mts. mostly by the paleo-Tisza system (Pogácsás 1984;Pogácsás and Révész 1987;Vakarcs et al. 1994;Magyar et al. 2013;Fig. 1A). The spatial and temporal evolution of the paleo-Danube system is relatively well known, and integration of outcrop and subsurface data into uniform sedimentological models was a subject of several studies (e.g. Sacchi et al. 1998;Magyar et al. 2007;Cziczer et al. 2009;Kováč et al. 2011;Sztanó et al. 2013a;Šujan et al. 2016;Sebe et al. 2020). The southwestward sediment transport by the paleo-Tisza and tributaries, however, was studied only in the deep basins of SE Hungary, such as the Derecske Trough (Vakarcs and Várnai 1991;Vakarcs et al. 1994;Lemberkovics et al. 2005;Balázs et al. 2016) and Békés Basin (Teleki et al. 1994;Csató et al. 2015;Fig. 1A). The sedimentary evolution of Lake Pannon in the NE Pannonian Basin in Slovakia, Ukraine, Hungary and Romania remained largely unexplored. Though scattered outcrop data are available from this area, no attempt has been made so far to integrate them into a depositional model, although sedimentary processes in this part of the Pannonian Basin played a key role in the early evolution of the lake.
The objective of this paper is the reconstruction of the sedimentary history of the Șimleu Basin, a northeastern peripheral depression of the Pannonian Basin in Romania, and a better understanding of its contribution to the largescale evolution of Lake Pannon. We conducted field work and collected well data in and around the Șimleu Basin to identify the major elements of the sedimentary system, the direction of sediment transport, and timing of the processes. By establishing correlation between this basin and the central, deep Derecske Trough, we explored the spatial and temporal relationship between sedimentary processes in a marginal and in a central depression of the same source-tosink system.

Geological setting
The development of the Pannonian Basin System is coeval with that of the surrounding Carpathian orogenic belt (Horváth et al. 2006). The opening of deep subbasins, controlled by low-angle normal faults and strike-slip faults, was related to the back-arc extension and rotation of blocks (Horváth and Royden 1981;Fodor et al. 1999;Cloetingh et al. 2005;Horváth et al. 2015;Matenco et al. 2016). Extensional deformation during the late early and early late Miocene caused the migration of depocenters in space and time (Horváth et al. 2006;Balázs et al. 2016), and it was followed by a similar shift in the formation of inversion-related structures (Bada et al. 2007). In contrast, the late Neogene evolution of the Transylvanian Basin was controlled by compressional structures, gravitational tectonics over mid-Miocene salt, and salt diapir development (Ciulavu et al. 2002;Krézsek and Bally 2006;Tiliţă et al. 2013).
The tectonic events triggered substantial paleogeographic re-organization on the surface. During the middle Miocene, the Pannonian Basin System was part of the large, inland Paratethys Sea. The uplift of the surrounding orogenic belt isolated the basin from the sea at the beginning of the late Miocene, and a large lake was formed (Rögl 1999;Magyar et al. 1999a;ter Borgh et al. 2013). The brackish-water Lake Pannon, characterized by an endemic biota, covered the Pannonian Basin System from the Vienna Basin in the west to the Transylvanian Basin in the east, encompassing an area of ca. 230,000 km 2 (Neubauer et al. 2016). During the first few million years of its geological history, the lake had a highly articulated shoreline with many islands and peninsulas. Most of them were later flooded and eventually buried 1 3 under sediments as river deltas advanced into the lake basin from the surrounding orogen.
The Şimleu Basin ( Fig. 1) is located between the central Pannonian Basin and the Transylvanian Basin. Geographically, it is bounded by the Plopiș and Meseș Mts. in the SW and SE, respectively, whereas its rolling hills pass into the flat landscape of the Great Hungarian Plain towards the E and NE (Fig. 1B). The present topography has a relatively high altitude of 120-320 m above the level of the Great Hungarian Plain as a consequence of the lithospheric-scale uplift of the Transylvanian Basin (Krézsek and Bally 2006;Tiliţă et al. 2013).
The evolution of the Şimleu Basin started with middle Miocene (Badenian) extension, resulting in deposition of foraminifera-bearing deep-water grey marls, while the overlying shallow-marine carbonates and siliciclastics represent the onset of post-rift subsidence. The marine sedimentary record is intercalated with felsic tuffs dated to 14.8-15.1 Ma (Szakács et al. 2012). The upper middle Miocene (Sarmatian) brackishwater deposits attain considerable thickness in some parts of the Şimleu Basin (Nicorici 1972;Clichici 1973). A short-term compressional event led to uplift and subaerial exposure, thus the upper Miocene (Pannonian) lacustrine deposits overlie unconformably on Sarmatian or crystalline basement rocks. As neotectonic uplift is still ongoing, the landscape is changing quickly; hence, landslides detaching in Pannonian clays often create new outcrops and cover up older ones.

Materials and methods
Standard field observations, wireline log interpretation, and biostratigraphic correlation were integrated into a basinscale stratigraphic and sedimentary model. More than 20 Traction on angle-of-response slopes outcrops were measured, and the facies were interpreted in terms of depositional processes (Table 1). The facies were grouped into facies associations to unravel architectural elements of the depositional systems (Table 2). Wireline logs (typically old-vintage GR, SP, and resistivity) from eight wells were digitized and sedimentologically interpreted. Characteristic log shapes of the typical Pannonian basin-fill formations were recognized (e.g. Juhász 1994; Pigott and Radivojević 2010;Sztanó et al. 2016). Finally, outcrop and well data from the Şimleu Basin were correlated through seismic profiles to the Der-I well in the central, deep Derecske Trough of the Pannonian Basin. Fifteen samples of fossil molluscs were collected from eight outcrops, and the archive collection of the Paleontology-Stratigraphy Museum of the Babeş-Bolyai University, Cluj-Napoca (BBU) from three localities was also studied. The 43 identified taxa are summarized in Supplement 2. Micropaleontological samples were processed with hydrogen peroxide from about 250 g of air-dried sediments. Twelve samples were barren of ostracods but eleven samples from five outcrops contained ostracod carapaces and single valves, representing 14 taxa (Supplement 3). SEM images were taken at the Botanical Department of the Hungarian Natural History Museum in Budapest.
In Cehei, FA1a unconformably overlies an erosional surface carved into the metamorphic basement with several meter relief, and it is capped by the horizontal beds of the muddy FA2. FA1a is a mixture of various sandy and gravelly facies units (Tables 1, 2). The most common are 5-20-m thick bed-sets of massive, poorly sorted, imbricated, clast-supported conglomerates (F15), coarsegrained, fossiliferous, pebbly sands (F12), planar and trough cross-stratified sand (F10), and several-cm-thick mudstone interbeds (F02). The bed-sets consist of laterally thinning beds with a thickness of a few decimeter, and dip towards the S-SW. Although there are lenticular units separated by erosional surfaces, most bed-sets show a convex geometry and downlap on older, flat bed-set boundaries or onlap on the slightly inclined ones (F02), revealing the delicate architecture of foreset terminations.

Interpretation
Debris flows, grain avalanches and grain flows alternated on the steep, gravelly, sandy foresets. Lobe switching resulted in deposition of mudstones, followed by the onlap/downlap of new lobes. Upwards, with decrease of dip angle and on the flat lying topsets, bedload was deposited as small dunes in shallow channels and as downstream migrating bars. This facies association is interpreted as steep front of Gilbert-type deltas and their delta plain (cf . Postma 1990;Nemec 1990;Gawthorpe and Colella 1990).
The Cehei outcrop reveals the toe and lower portion of foresets, which were covered by offshore mudstones when base-level rise outpaced the locally high sediment input. In the Sâg outcrop, the sediment transport mechanism was dominated by traction of bed-load, evidenced by the dm-scale trough cross-stratification and the low-angle crude-stratification of gravel and sand. The a(p)a(i)-type imbrication in the gravel indicates mass-flow events (cf. Postma 1990). The interfingering of gravel with crossstratified and sheet-like sands demonstrates the lateral variability of the mouth bars near the topset-foreset transitional zone. Upper-flow regime structures and erosional surfaces point to the effect of occasional floods. Unidirectional currents created shallow channel-form architecture. Wave-induced ripples also occur in this shallow-water topset environment.
The foresets of the Sâg Gilbert-type delta and the channel-fill sediments show the same transport direction towards the N-NW. Based on the progradational directions and the mixed composition of the clasts, the sediment source was in the S-SE, in Plopiș Mts.
The Porț and Cehei outcrops show a slightly different depositional setting. The moderately sorted monomictic gravelly bed-sets in Cehei show convex downlapping geometry and can be interpreted as a lower segment of deltaic lobes, below wave-base, near the foreset-bottomset transitional zone. In Porț, the pebbly sand and massive graded pebbly sandstone of the delta slope might have been transported by shortdistance gravity-flows initiated by fluctuating discharge and sediment charge of the feeding local fluvial systems (Gobo et al. 2015).
Fossils of the littoral dwellers Melanopsis fossilis and Congeria hemiptycha (e.g. Pavlović 1928) in the sand of the Cehei outcrop were probably reworked from the eroded topsets. The molluscs and ostracods of the clay-rich intercalations indicate a sublittoral environment with low-energy conditions (e.g. Cziczer et al. 2009). These deposits formed by fallout from suspension and capture periods when active transport occurred on other lobes in laterally offset positions.
A more than 40-m-thick succession of FA2 is exposed in the Varșolț claypit (Fig. 4), where neither the under-nor the overlaying beds are exposed. It consists of brownishgrey, cm-thick, laminated or bluish grey, structureless clay beds (F01), which contain mm-thick silt and up to 20 cmthick very fine sandstone intercalations (F02, F03, F07, F09) with erosional base, normal gradation and planar-to cross-lamination. Sharp-top sandstone intercalations and rare clay-filled vertical burrows of mm size also occur. The frequency and thickness of silty-sandy interbeds increase upwards (Fig. 4A, B). The paleocurrent directions measured from ripples indicate a NW-SW sediment transport.

Interpretation
Based on the sedimentary facies, FA2 was deposited relatively far from the sediment supply system, where deposition was mostly related to suspension fallout. The offshore mudstones in the upper part of the Cehei outcrop were formed when base-level rise outpaced the locally high sediment input, thus transgression of the shores resulted in the full flooding of the local island.
The clayey succession at Varșolţ points to different events. The laminated character of the clay, as well as the low abundance of molluscs and ostracods indicate alternations of dysoxic and aerated bottom conditions. The mollusc fauna consists of species that were adapted to the profundal habitats of Lake Pannon (e.g. Geary et al. 2000). The dominance of early juvenile ostracod specimens over adult ones is commonly observed in profundal zones of recent lakes, due to the high mortality rate of instars, transported from the shallower region, under the cold, temporarily oxygen-depleted bottom water conditions (Zhai et al. 2015).
The thin silty to sandy intercalations and the normally graded sand beds are interpreted as products of low-volume, low-density, silty to sandy turbidity currents. Nongradational sharp-top sandstones may indicate bypassing high-energy turbidity currents. The upward increasing frequency and thickness of the turbidites point to gradual increase of clastic sediment input at the Varșolț section. The continuous sheet-like geometry of the sand beds may reveal either an area very far from distal lobes on the basin floor, or an overbank far from the main slope-channel. The latter option is supported by regional dip observations and the stratigraphic position between outcrops of turbidites and delta-lobes.

Description
FA3 is formed by several-meter-thick successions of sandstones (F09 and F11), conglomerates (F13) and heterolithic alternations of thin sandstones (F03, F07) and mudstones (F01) (Table1). Panic-S and the lower part of Panic-N are sand dominated with a set of related structures, including erosive base, normal gradation, structureless sand containing rip-up mud clasts, horizontal-and cross-lamination, climbing cross-lamination and waterescape structures. Lateral and vertical facies changes are common in the upper part of Panic-N ( Fig. 6A-C), consisting of heterolithics (F01, F03), medium-bedded sandstone laterally transforming into thin-bedded sandstone (F07) and matrix-supported mud-clast conglomerates (F13). The sandstones show channel-form geometry: their width varies between 2 and 10 m, while the depth of scouring is up to 1 m. The large (1-40 cm in diameter), rounded to angular, imbricated rip-up mud clasts, floating in medium-to coarse-grained sandstone matrix, are situated above the erosional scours. The mud clasts consist of the heterolithics (F01, F03) and display soft-sediment folding.
Only a very few, small, indeterminate shell fragments of molluscs and plant remains were found in these sands. Paleocurrent measurements from ripples indicate a westward sediment transport (Fig. 6).

Interpretation
Most facies in the Panic outcrops (F03, F07, F09, F11, and F13) reveal sediment gravity flows. Very thin-, mediumand thick-bedded turbidites, some with Ta, Tcd or Tabcd Bouma-sequences indicate transport by low-and highdensity turbidity currents, respectively. As bed thickness decreases, Tcd becomes more common. Thick-bedded mud-clast conglomerates are interpreted as debrites occurring usually not far from the slope (cf. Lowe 1982;Walker 1992;Talling et al. 2012). Low-relief erosional scours, thick, amalgamated, and laterally continuous beds with dewatering structures and climbing ripples indicate rapid accumulation of the sandy sediments, and point to a proximal lobe axis either between channels or in a very shallow aggrading channel. In contrast, the upper part of Panic-N succession is characterized by small yet distinct erosional scours, with a complex channel-fill lithology.
The scours were formed by bypassing high-energy erosive currents. High-density currents eroded mud clasts from the previous channel margins or inner levees, and suppressed turbulence led to rapid freezing and deposition (cf. Johnson et al. 2001). The overlying normally graded, mediumbedded turbidites are associated with waning currents (cf. Kneller and Branney 1995). Thin-bedded turbidites indicate the termination of the filling phase. These repeated cut-and-fill successions are interpreted as channel-fill storey elements developed by cyclic waxing and waning flows driven by climate forcing in the endorheic lake (cf. Tőkés et al. 2021).
Taking into account the presence of massive lobate sands associated with scour-fills of mixed origin and internal levees, we assume that the deep-water sediments near Panic were deposited in the proximal part of a turbidite system, probably close to the channel-lobe transitional zone (e.g. Postma et al. 2016;Brooks et al. 2018).

Description
Deposits of FA4 form 5 to 15 m thick, sand-dominated, coarsening upward successions. In Nușfalău, Pericei and Zăuan, sandstone bed thickness varies between a few centimeters and a few meters, while mudstones are a few dm thick. The fine-grained intervals consist of horizontally laminated silty mudstone (F01) and interbeddings of very fine-grained sand (F02, F04, F06) with abundant small shell fragments and vertical burrows (Fig. 7F). In the overlying, commonly wedge-shaped sandstones, a wide variety of sedimentary facies occur. The thin beds show planeparallel lamination (F02) and symmetrical and asymmetrical types of cross-lamination (F06) (Fig. 7G), which were formed by rapid flows and wave and current ripples, respectively (Fig. 7E, F). Intensity of bioturbation varies in these beds, from distinct, sand-filled, simple vertical burrows of dm length and cm diameter, through moderate density of burrowing to full homogenization of the sediment. The thick beds show tabular and trough cross-bedding (Fig. 7H) (F08, F10), but structureless, amalgamated, i.e., fully bioturbated beds also occur (Fig. 7F). Small mud clasts, plant debris, and shell fragments are common at the lower part of foresets (Fig. 7D).
The Camăr, Sălăjeni, and Ip outcrops are characterized by thick, tabular and trough cross-bedded, mediumgrained sandstones (F08, F10) containing mud clasts and coquinas at their base. These sand bodies display a lowrelief channel shape, and interfinger with parallel-laminated fine-grained sand.

Interpretation
FA4 is the most commonly exposed facies association in the Şimleu Basin (Table 2.). The sand-dominated, coarseningupward sequences in Nușfalău, Pericei and Zăuan indicate an upward increasing water agitation by currents and waves, reflecting a change from the transitional zone to the foreshore, with the dominance of current-induced structures. The downlapping or wedge-shaped bed-sets with laterally varying thickness indicate lobate geometry, which is typical for the shallow-water mouth-bars of deltas (Martini and Sandrelli 2015). The alternation of thin sandy and muddy beds, associated with different ripples and burrows, indicates the distal/lower delta-front, where the combined effect of waves and river floods prevail (Plink-Björklund and Steel 2004; Jorissen et al. 2018). The thickening of sand beds reflects delta front progradation as the increasing effect of river discharge (Olariu and Bhattacharya 2006). Sandy successions from erosional scours, coquina or mud-clast accumulations, overlain by thick cross-bedded sandstones represent channels on the mouth bars (cf. Jordan and Pryor 1992;Olariu and Bhattacharya 2006). The effect of waves was more pronounced on the lower delta front than on the upper one, therefore these deltas are classified as river-dominated. Yet, the parallel-laminated, sheet-like sandstones at Bilghez deposited on a high-energy wave-modified bar top (cf. Reading and Collinson 1996;Schomacker et al. 2010), pointing to local variations.
The mollusc fauna of FA4 indicates some variety and patchiness in the depositional environment. The Ip fauna suggests a vegetated nearshore environment based on the dominance of small-sized melanopsids (e.g. Müller and Szónoky 1990). In Nușfalău, Sălăjeni, and Zăuan, higherenergy environment can be inferred from the abundance of the large-sized and thick-shelled molluscs (M. vindobonensis, C. hemiptycha, and L. cf. schedelianum). The Nușfalău fauna also includes Melanopsis with irregular, knobby sculpture, described by Pavlović (1928) from Vrčin, Serbia as M. vindobonensis karagacensis (Fig. 5i). We assume that this unusual pattern is a paleopathological feature rather than an ecological response.
The ostracod genus Cyprideis is known to prefer the littoral zone in marine and lacustrine environments (e.g. Morkhoven 1963;van Harten 1990;Boomer et al. 2005;Beker et al. 2008;Stoica et al. 2013). Its massive, less ornate carapace may reflect adaptation to high-energy bottom water conditions and to the sand-dominated substrate.

Facies association 5 (FA5a and FA5b): lower and upper delta plain (Derna, Cuzap, Șărmășag, and Voivozi)
Description FA5 is composed of alternating mud-and sand-prone facies, with very thin lignites in FA5a and mineable ones in FA5b. In the central part of the basin, at the outcrops of Nușfalău, Camăr, Bocșa, Șărmășag and Ip, FA5a overlies FA4 and consists of alternating layers of clay (F01), laminated silty sand (F04), and cross-laminated fine-grained sand (F06) with 2-3 mm wide and a few cm long, sand-filled vertical burrows. It comprises a few meter-thick, coarsening upwards units, topped by organic-rich clays or cm-thick lignite seams (F05) ( Table 1).

Interpretation
The mudstones with thin sandy interbeds filled the interdistributary bays of the lower delta plain during subsequent floods, resulting in the formation of swamps (Table 2). Thick cross-bedded sands formed as active distributary channel fills, while organic-rich deposits represent abandoned channels, ponds or swamps on the upper delta plain. Distribution of sediments and thickness of lignites on the delta plain was strongly influenced by compaction-related subsidence (cf. Phillips and Bustin 1996;Reading and Collinson 1996).
A specimen of the pearl mussel Pseudunio flabellatiformis from Derna (BBU #15.273) indicates a freshwaterfluvial environment. Lubenescu et al. (1967) reported a small fauna with Tinnyea vasarhelyii and small-and large-sized melanopsids from the Șărmășag locality, suggesting a wellvegetated environment with changing salinity conditions.

Facies association 6 (FA6): fluvial channels (Corund, Derșida, and Cohani)
Description FA6 is composed of medium-to coarse-grained sandstone and subordinate mud layers. The Cohani outcrop displays a 20-25 m thick succession of laterally discontinuous, trough cross-bedded, coarse-grained sandstone (F10) and tabular cross-bedded sandstone (F08) packages. Thickness of the bed-sets varies between 0.4 and 1.5 m. Above the basal surfaces, imbricated pebbles occur. Some thin-bedded, silty, very fine sand beds (F06) are associated with completely weathered tuff horizons. Pieces of silicified wood fragments also occur. Sandstone series are bounded by major erosional surfaces with a relief of 2 m. Similar cross-bedded sandstone occurs in the small Corund outcrop.

Interpretation
The cross-bedded sets, separated by erosional surfaces, represent vertically and laterally stacked, aggradational multistorey fluvial channel-fills (cf. Collinson 1996). Each storey contains channel-fill bars separated by minor erosional surfaces. The fine-grained sediments are restricted to the bar tops. Mudstones of the surrounding floodplains are not exposed.

Correlation of facies associations to the subsurface
The depositional systems represented by the exposed facies associations are also recognized in the subsurface. In the Săcueni, Bobota, Mișca, Tăuteu and Chiraleu wells, serrated and blocky log shapes indicate coarse sand and conglomerate of up to 20-50 m thickness, overlying the crystalline basement or middle-Miocene deposits. These correspond to FA1 recognized in outcrops. Elsewhere the lacustrine sedimentary record starts with about 30 m thick offshore marls (FA2).
Above the marls, up to 30 m thick blocky and barrelshaped sandstone sequences separated by 2-10 m thick mudstones occur in the Zalău, Nușfalău and Săcueni wells. These are interpreted as sand-prone lobes of a turbidite system. The corresponding FA3 in outcrops of Panic show the channel and/or channel-lobe transitional zone.
The overlying, up to 300 m thick unit, being present in all wells of the study area, displays smooth or finely serrated log shapes. It is interpreted as mudstones deposited on the shelf-edge slope (FA2). It thickens to 400-700 m towards the basin interior in the west west (Săcueni, Chiraleu).
Funnel log shapes correspond to 20-50 m thick coarsening upward successions and reveal prograding sandy delta lobes (FA4 and FA5) representing the upper part of the basin fill.
The youngest unit of the upper Miocene lacustrine sequence shows a serrated log pattern interrupted by up to 20-m thick intervals of bell and blocky log shapes. These are interpreted as representing a clayey-silty alluvial plain, cross-cutted by sandy channel-belt deposits (e.g. FA6, Cohani). The thickness of this unit increases to 500 m westwards in the Cheţ and Derecske Troughs (e.g. Chiraleu).

Biostratigraphy and age
Earlier mollusc biostratigraphic studies from the Şimleu Basin correlated the lacustrine succession with the lower part of the Pannonian, stressing that the oldest Pannonian biozones are missing in the area (Matyasovszky 1879;Lubenescu et al. 1967;Marinescu 1985; but see Chivu et al. 1966).
The low-diversity mollusc fauna recovered from the coarse-grained deltas in Cehei is indicative of the Lymnocardium conjungens Zone (11.0-9.6 Ma; Fig. 9). The bivalve Congeria hemiptycha (also known by its junior synonym as Congeria pancici Pavlović 1927) first appears in "Zone D" of the Vienna Basin (Papp 1985), which was dated 10.6-10.4 Ma by Harzhauser et al. (2004). The Fig. 9 State-of-the-art biozonation of Lake Pannon deposits with the age of the studied Pannonian localities of the Șimleu Basin. Mollusc zonation follows Magyar and Geary (2012), ostracod zonation is modified after Krstić (1985). Paleomagnetic chrons and European Neogene mammal biozones are used after Hilgen et al. (2012). Blue line indicates profundal, green sublittoral, red littoral, and purple freshwater-fluvial environment. Turquoise bracket shows the 10.6-9.6 Ma time interval supposed for the main transgressive-regressive cycle in the Șimleu Basin age of the Cehei conglomerates and sandstones is thus restricted to the interval of 10.6-9.6 Ma (Fig. 9).
Based on the common occurrence of Lymnocardium conjungens, L. hantkeni, Congeria hemiptycha, Melanopsis vindobonensis etc. (see in Lőrenthey 1893;Strausz 1941;Lubenescu et al. 1967;Nicorici and Karácsonyi 1983), the shallow-water lacustrine deltaic sands (FA4 and FA5) belong to the upper part of the Lymnocardium conjungens Zone (ca. 10.2-9.6 Ma; Fig. 9). The appearance of Dreissena auricularis in these assemblages, however, points to the youngest part of this zone. Papp (1950) postulated that D. auricularis evolved from Congeria gitneri, and the stratigraphic distribution of these two species in the western part of the Pannonian Basin is in accord with this hypothesis (Magyar et al. 1999b). The earliest occurrence of D. auricularis in Pezinok (Horusitzky 1907) was dated by mammal stratigraphy as early MN10, ca. 9.6-9.7 Ma (Joniak 2016).
The ostracod Amplocypris abscissa, occurring in FA1 and FA2 in Cehei and in FA4 in Nușfalău, was considered by Krstić (1985) as a biostratigraphic marker of the Amplocypris abscissa and the overlying Hemicytheria croatica zones. A. abscissa occurs in the Hennersdorf outcrop next to Vienna (Danielopol et al. 2011), which was interpreted by Harzhauser et al. (2004) as being 10.2-10.4 Ma old.
The most common fluvial mussel, Pseudunio flabellatiformis, is well-known from the fluvial deposits of the Pannonian Basin throughout the entire late Miocene (as "Unio wetzleri" or "Margaritifera flabellatiformis" in the older literature; here we follow the revision of Lyubas et al. 2019); thus, it does not constrain the age of FA5 and FA6 occurrences. In Derșida and in the lignite-bearing deltaic layers of Derna, however, Anancus arvernensis remains were found (Codrea and Margin 2009;Gasparik pers. comm. 2020). The oldest known European occurrences of this large-sized mastodon are Messinian (MN12, from 7.2 Ma on; Konidaris and Roussiakis 2018). This dating implies a significant stratigraphic gap between the older-than-9.6-Ma lacustrine and the younger-than-7.2-Ma fluvial deposits, but this gap is not supported by our field observations. This contradiction thus remains unresolved until a more detailed study of the unconformities in the area or a revision of the large mammal stratigraphy brings along new perspectives.

Changes of the depositional environment
The oldest lacustrine sediments exposed in the study area are the coarse-grained conglomerates and sandstones, deposited in locally-sourced progradational Gilbert-type deltas that rimmed the islands during the initial flooding of the Șimleu Basin (Fig. 10A). Their age is known only at Cehei (10.6-9.6 Ma), but they could have developed diachronously along the margins.
As base-level rose and islands submerged, offshore mudstones were deposited by suspension fallout. These preserved the fossils of the profundal Congeria banatica fauna, which was widespread in Lake Pannon between ca. 11 and 9.6 Ma (e.g. Botka et al. 2019). The low abundance and low diversity of the fauna indicate a partially isolated, oxygen-depleted deep-water environment (Fig. 10B).
The first distally sourced clastics, indicating the onset of regression in the Șimleu Basin, were thin-bedded turbidites in the upper part of the profundal mudstones. The coeval slope might have prograded from E or SE (Fig. 10B), implying that the Meseș Mts. was not a confining barrier at that time: the central Pannonian Basin and the Transylvanian Basin were geographically connected through the Șimleu Basin.
A more robust evidence of the ESE to WNW advancing slope progradation is provided by the sandy turbidites at Panic, which can be interpreted as lobe axis and/or lobe-channel fill deposits (Fig. 10C). Few 10 s of metre thick stacked turbidite sequences also appear in wells near Zalău, Nușfalău, Bobota, Chiraleu, and Săcueni. The relatively small thickness, coupled with the probably large extent of these sand lobes is typical for unconfined toeof-slope turbidite systems (cf. Sztanó et al. 2013b;Tőkés and Patacci 2018).
Turbidites are overlain by 200-300-m thick mudstones (Nusfalău and Bobota wells), which can be attributed to the slope. A seismic section from Căuaş, 25 km to the N (Fig. 7 in Ciulavu et al. 2002), displays ca. 150 m high clinoforms. Progradation of similar, relatively small clinoforms might have filled the Şimleu Basin rapidly. The advance of the shelf-edge was supplied by river-dominated prograded from E to W across the Simleu Basin not later than 9.6 Ma ago. E Paleotransport directions based on field data in the Şimleu Basin. F Timing of shelf-edge progradation in Lake Pannon (after Magyar et al. 2013) deltas at around 9.7-9.6 Ma. Measurements indicate paleotransport directions to WNW and NW, i.e. towards the central parts of Lake Pannon (Fig. 10E).
Based on outcrop data, the thickness of individual deltaic parasequences varies between 3 and 10 m, and they build up to 50 m thick deltaic sequences displayed by well-logs. The height of the parasequences and deltaic bodies in the Șimleu Basin is in the same range as in the paleo-Danube or paleo-Tisza systems (Juhász 1994;Sztanó et al. 2013a;Magyar et al. 2019). These packages indicate moderate-amplitude short-term base-level changes and/or autocyclic avulsions of the feeding rivers.

Correlation and comparison with the central Pannonian Basin
Both the eastern part of the deep, central Pannonian Basin and the Şimleu Basin were filled by sediments derived from the Eastern Carpathians. Their evolution was parallel, but infill of the Pannonian Basin lasted for ca. 7 million years (Magyar et al. 2013), while that of the Şimleu Basin took only less than 1 million years. To understand the role that the Şimleu Basin played in the infill of the large system, correlation was established to the Derecske Trough, one of  1 3 the deepest depocenters of the Pannonian basin ca. 70 km to the W of the Şimleu Basin (Fig. 11). The base of the Neogene gradually ascends from several hundred meter depth in the Şimleu Basin to more than 5 km in the Derecske Trough. The largest throw is observed at the Suplacu de Barcău fault zone, which played a crucial role in the accumulation of the Suplacu oil field, one of the largest in the Pannonian Basin (e.g. Panait-Patica et al. 2006).
Exploration wells in the region are generally located above the basement highs, such as the Chiraleu and Săcueni highs, whereas the basin interiors are mostly known from seismic data (Fig. 11). The 5200 m deep Derecske-I hydrocarbon exploration well ("Well-5" in Vakarcs et al. 1994) was drilled close to the center, yet the deepest part of the basin was not reached (Balázs et al. 2016).
The late Miocene lacustrine succession in the Derecske-I well (Table 3) starts with deep-water marls. The overlying basin-center confined turbidite system is made up of lobe complexes with paleotransport directions from NE to SW (cf. Sztanó et al. 2013b). These are overlain by slope deposits, i.e. thin unconfined turbidites and monotonous mudstones. The slope related clinoform height is ca. 650 m in the Derecske Trough (Balázs et al. 2018). It decreases to ca. 500 m in the Cheţ and Abrămuţ basins in the east, whereas in the Şimleu Basin it is not higher than 200 m. Considering compaction, the tallest slopes during deposition could have exceeded 1000 m height in the Derecske Trough (Balázs et al. 2018) and 300 m in the Şimleu Basin. Deltaic successions are built up of 30-50 m thick, coarsening-upward lobe cycles. The overlying late Miocene to Quternary fluvial suit is extremely thick in the Derecske Trough. The elements of the lacustrine sedimentary succession are almost identical in the Derecske Trough and Şimleu Basin, except for their thickness and timing.
The clinoforms provide another aspect for comparison. The direction of progradation was from E-SE to W-NW in the Şimleu Basin, as measured in the field within the turbidites and delta sediments. In contrast, in the eastern depocenters of the central Pannonian Basin, on both sides of the Hungarian-Romanian border, slope progradation was from N-NE to S-SW, as evidenced by seismic clinoforms (Tulucan 2007;Răbăgia 2009;Horváth and Pogácsás 1988;Vakarcs et al. 1994;Lemberkovics et al. 2005;Magyar et al. 2013). This pattern indicates a mergence of different fluvial feeder systems immediately W of the Şimleu Basin. Similar features were observed in other parts of the Pannonian Basin as well (Pogácsás and Révész 1987;Magyar et al. 2013).
In an earlier study (Magyar et al 2013), seismic surfaces were calibrated as biozone boundaries across the Pannonian Basin, and they were dated according to the biochronostratigraphic system of Magyar and Geary (2012) (Fig. 9). The 8.6 Ma surface, corresponding to the base of the Lymnocardium decorum Zone, reaches the shelf edge 4 km to the SW of Derecske-I (Fig. 10F). The correlation of this surface towards the Şimleu Basin is straightforward up to the Suplacu fault zone; coeval deposits in the Şimleu Basin, however, have been eroded (Fig. 11). This correlation sets a lower limit (a minimum) on the age of the lacustrine and deltaic succession of the Şimleu Basin. Thus, there is a 1 million-year time shift in the formation of the shelf edge and related facies associations between the Şimleu Basin and Derecske-I well.
The shelf-edge advanced across the Şimleu Basin at about 9.7-9.6 Ma, probably in a very restricted time interval (a few 100 kys at most). Between 9.7 and 8.7 Ma, the shelf-edge advanced ca. 70 km to the west, across the Abrămuţ Basin and Cheţ Trough to the Derecske Trough, into increasingly deeper water. This rate is similar to the rates observed in other parts of Lake Pannon (cf. Sztanó et al. 2013b) and in supply-dominated, moderately deepwater marine settings worldwide ). The shelf accretion was fed by deltas across the Şimleu Basin (at least up to 9.6 Ma). When delta lobes filled the available accommodation to lake level in the Şimleu area, deposition of turbidites occurred in a water depth somewhat larger than 1000 m in the Derecske Trough. Thus, the late Miocene paleorelief can also be assessed from these data. The roughly estimated sedimentation rates in the studied successions of the Şimleu Basin and in the basin interior are on the same order (Table 3).
The outcrops of the relatively thin Pannonian succession of the Şimleu Basin thus display the whole variety of processes that shaped the Pannonian Basin, including deep-water turbidite and slope formation. This pattern is noteworthy because most Pannonian outcrops expose only a limited segment of the basin fill. The most commonly exposed units are either coarse-grained deltas at the base of the succession (Sztanó et al. 2010;Budai et al. 2019) or deltaic and fluvial deposits (e.g. Kováč et al. 1998Kováč et al. , 2018Sacchi et al. 1998;Uhrin and Sztanó 2007;Sztanó et al. 2013a;Magyar et al. 2017;Pavelić and Kovačić 2018), neither of which represent the full history of basin evolution. Generally, deep-water deposits, such as offshore marls and turbidites, are particularly rare in surface outcrops (Kovačić et al. 2004;Krézsek et al. 2010;Bartha et al. 2015;Tőkés et al. 2021). The Şimleu Basin owes its uniqueness to its paleogeographic position close to the margin of the Pannonian Basin but well within Lake Pannon, to its subsidence and deepening that was mostly coeval with that of the main basin, to its size and topography that allowed the formation of a deep-water environment, and to its early sedimentary infill that took place soon after the onset of normal regression of Lake Pannon.

Conclusion
The upper Miocene facies associations in the Şimleu Basin, a northeastern marginal depression of the Pannonian Basin, represent a major transgressive-regressive cycle deposited in Lake Pannon. The Şimleu Basin is unique, because each important element of a source-to-sink system was present in this relatively small and shallow basin. Fluvial, deltaic, shelf-edge slope, and open-to-deep lacustrine depositional settings were identified in outcrops and shallow boreholes. For instance, coarse-grained deltas, deposited above active fault scarps during the diachronous initial transgression of Lake Pannon, and deep-water turbidites are equally accessible in outcrops. Biochronostratigraphic considerations constrain the age of the lacustrine and deltaic sediments of the Şimleu Basin to 10.6-9.6 Ma. Compared to the upper Miocene succession of the Derecske Trough, an adjacent deep sub-basin of the central Pannonian Basin, the analogous units are an order of magnitude thinner in the Şimleu Basin, yet sedimentation rate was similarly high. Having been located closer to the source, however, the shelf-edge and related facies associations are one million years older in the Şimleu Basin than in the Derecske Trough. Both basins were filled by sediments derived from the Eastern Carpathians, and both funneled them further to the deep basin interior.
Funding Open access funding provided by Eötvös Loránd University.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.