Tracing the cryptic Sardic (Ordovician) metamorphism across Alpine Europe: the Krndija region in the Slavonian Mountains, Croatia

Results of a combined petrological, geochemical and geochronological study suggest that metasedimentary rock units in the Krndija region of the Slavonian Mountains, Croatia, were affected by at least three major tectonometamorphic imprints: during the Middle Ordovician (Sardic event), the early Carboniferous (Variscan event), and the Cretaceous (Alpine event). All three metamorphic phases are established by electron microprobe-based in-situ U–Th–Pb dating of monazite grains. The Sardic metamorphic event is additionally confirmed by a precise Lu–Hf garnet-whole-rock isochron age of 466.0 ± 2.3 Ma. Taken together, the data unveil a relatively large and well-preserved piece of the cryptic Sardic orogen in central Krndija, that we name the Kutjevo Zone. A Sardic subduction-related metamorphic event (ca. 540-580 ℃, 8–11 kbar) at ca. 466 Ma is manifested in the mineral paragenesis Ca-rich garnet plus rutile. A low degree of retrograde reequilibration suggests a subsequent fast exhumation. Low-Ca cores in some garnets and staurolite relics record a pre-HP metamorphic event that involves isobaric heating from 570 to 610 ℃ at ~ 7 kbar. We attribute this (so far undated) event to mid-crustal contact metamorphism caused by early Sardic magmatism. Southern parts of Krndija (the Gradište Zone) experienced an (additional?) clockwise PT evolution in Variscan times at ca. 350 Ma. Garnet formed with ilmenite during a PT increase from 580 ℃/5 kbar to 600 ℃/6 kbar and underwent later strong retrograde resorption. Slow Variscan exhumation resulted in andalusite formation at < 550 ℃/ < 3.8 kbar. Penetrative Alpine metamorphism was observed in low-grade phyllites in the north. The lithology and metamorphic history of the Kutjevo Zone is similar to what has been reported from the Sardic Strona-Ceneri Zone in the western Alps. Both areas expose metapelitic (metagreywacke) rocks with a pre-middle Ordovician formation age. These metasedimentary rocks are inter-layered with numerous small amphibolitic units as well as metagranitoids and were likely deposited along the active Gondwana margin, perhaps in a fore-arc position, prior to their subduction during the middle Ordovician. According to recent palaeogeographic reconstructions, both the Kutjevo Zone and the Strona-Ceneri Zone have once resided in an eastern sector of the northern Gondwana margin (i.e., in E-Armorica). We conclude that in the Middle Ordovician, important subduction activities took place in this E-Armorican segment of north Gondwana, which is today exposed in the Alps. The W-Armorican segment of north Gondwana (now exposed in the French, German, and Czech Variscides) had probably already mutated from a (Cadomian) subduction setting to an extensional (transtensional–transpressional) setting by the late Cambrian.


Introduction
Significant relics of pre-Permian basement can be found in many places throughout the Alpine orogen (e.g., Schaltegger and Gebauer 1999;Thöni 1999;Eichhorn et al. 2001;von Raumer et al. 2013;Neubauer et al. 2022). These rocks constitute important mosaic pieces in the reconstruction of the Palaeozoic evolution of central Europe. In simple terms, the pre-Permian basement rocks of the Alps are often referred to as "Variscan basement", assuming that they formed (or at least received a major tectonometamorphic imprint) during the Variscan orogeny. However, geochronological research increasingly shows that many of these rocks have a prominent Early Palaeozoic metamorphic and magmatic history (e.g., von Quadt 1992; Bussy and von Raumer 1994;Poller et al. 1997;Zurbriggen et al. 1997;Schulz et al. 2004;Franz and Romer 2007;von Raumer et al. 2013;Neubauer et al. 2022).
The Early Palaeozoic rocks represent important relics of the northern Gondwana margin, from where the central European crust is mainly derived (Linnemann et al. 2004;Nance et al. 2008;von Raumer et al. 2013). A widely accepted assumption is that the northern Gondwana margin was in a state of extension in the Early Palaeozoic (Franke et al. 2000;Franke et al. 2017) and was affected by high heat flow and rift-related magmatism (Pin and Marini 1993). However, the existence of Ordovician Barrow-type metamorphism in Alpine basement units (e.g., Biino 1995;Zurbriggen et al. 1997;Faryad et al. 2002;Franz and Romer 2007) implies that subduction processes must have played a role as well, at least locally. These Ordovician subduction activities, summarized as the Sardic tectonometamorphic phase by Handy et al. (1999), constitute a cryptic feature in the evolution of the central European crust. The debate on the significance of an orogenic phase of "Caledonian" age in Variscan Europe was first introduced by Stille (1939) to account for a Middle-Ordovician gap in the stratigraphy of SW Sardinia (Teichmüller 1931).
In this paper, we introduce another key area for the study of this "pseudo-Caledonian" orogenic phase, namely the Krndija region of the Slavonian Mountains, Croatia. Our multidisciplinary geochronological-petrological-geochemical study clearly reveals that this area contains well-preserved Sardic basement, which shares many similarities with the Strona-Ceneri Zone, the type locality for Sardic metamorphism in the Alps. Parts of Krndija were also affected by strong Variscan and by Alpine metamorphism. The systematic combination of in-situ monazite dating with thermobarometry (Schulz 2021) allows us to resolve the polymetamorphic history of this part of the Slavonian Mountains.

Geological background
The Slavonian Mountains in Croatia, subdivided geographically into the Psunj, Papuk, and Krndija regions (Fig. 1, inset), form one of the largest crystalline basement complexes within the Tertiary Pannonian Basin (Fodor et al. 1999;Pamić and Jurković 2002). In terms of the Alpine plate tectonic framework, they belong to the Tisia-Unit. This tectonic unit lies adjacent to the Southern Alps and the Dinarides (Fig. 1), and it is commonly regarded as a lithospheric fragment broken off from the European plate during the Middle Jurassic opening of the eastern branch of the Alpine Tethys (Géczy 1973;Csontos 1995;Schmid et al. 2008). The Tisia-Unit came into its present-day position after complex translational and rotational movements during the Mesozoic and Cenozoic (Csontos 1995;Fodor et al. 1999;Csontos and Vörös 2004). These motions were controlled by the major tectonic contact zones of the Alpine-Carpathian-Dinaridic orogenic system, most of them representing oceanic sutures (Schmid et al. 2008). The internal structure of the Tisia-Unit encompasses three internal, southward dipping, Alpine nappe systems, called Mecsek, Villány-Bihor, and Békés-Codru, each comprising igneous and metamorphic basement rocks and post-Variscan overstep sequences (Haas and Péró 2004;Csontos and Vörös 2004;Schmid et al. 2008). The Slavonian Mountains belong to the Villány-Bihor nappe system (Fig. 1, inset).
The main crystalline lithologies of the Slavonian Mountains are mica schists, paragneisses, amphibolites, migmatites, and granites (Jamičić and Brkić 1987;Jamičić 1989). The abundant granitic rocks ( Fig. 1) are generally considered to be Variscan intrusive rocks, based on field relationships and geochronological data (Pamić et al. 1988;Pamić and Lanphere 1991;Pamić and Jurković 2002). The formation ages of the metamorphic rocks ( Fig. 1) have been controversially assessed in the past. Historically, an Early Palaeozoic/ or late Neoproterozoic (Baikalian) age and a pre-Variscan metamorphism were envisaged for, at least, some of the rocks based mainly on stratigraphical constraints (Jamičić 1983(Jamičić , 1989Jamičić and Brkić 1987). In more recent work, it was mostly suggested that the metamorphic evolution of the Slavonian Mountains is essentially Variscan (Pamić et al. 1988;Pamić and Lanphere 1991;Pamić and Jurković 2002). The significance of a very low-to low-grade metamorphic overprint during the Cretaceous has also been recently highlighted (Biševac et al. 2009(Biševac et al. , 2011Balen et al. 2013;Starijaš Mayer et al. 2013).

Results of previous studies
Our current knowledge about the metamorphic rocks of the Slavonian Mountains relies mainly on the detailed mapping work of Jamičić (1983Jamičić ( , 1989 and Jamičić and Brkić (1987), as well as several petrographic studies by Pamić and coauthors (e.g., Pamić et al. 1988Pamić and Lanphere 1991). Jamičić (1983Jamičić ( , 1988 distinguished three tectonometamorphic units in the Slavonian Mountains, termed the Radlovac, the Psunj, and the Papuk metamorphic complexes (Fig. 1).
The Radlovac metamorphic complex comprises a very low-grade metamorphic sequence of slates, metagreywackes, metaconglomerates, and subordinate phyllites, with some intrusive bodies of metadiabase and metagabbro (Pamić and Jamičić 1986). Recent studies have shown that the regional metamorphic overprint is most probably of Cretaceous (Alpine) age (Biševac et al. 2009(Biševac et al. , 2011. Chloritoid schists occurring at the base of the Radlovac complex formed at conditions of 340-380 ℃/3.5-4.0 kbar . The rocks of the Radlovac complex overly the higher grade Psunj metamorphic complex and are themselves unconformably overlain by a clastic-carbonate succession of Late Permian and Triassic age (Jamičić 1983;Jamičić and Brkić 1987).
The Papuk metamorphic complex consists of migmatites and migmatic gneisses arranged around a core of S-type granite (Pamić et al. 1988(Pamić et al. , 1996Pamić and Lanphere 1991). The migmatitic rocks grade into an amphibolite facies metamorphic sequence mapped as biotite muscovite gneisses, which also contain garnetiferous amphibolites, paragneisses, and mica schists. According to Pamić and Lanphere (1991), the metamorphic evolution in Papuk is contemporaneous with that in Psunj, representing deeper Variscan crust affected by migmatization and granite formation.
The concept of a monocyclic Variscan metamorphic evolution of the entire Slavonian Mountains was challenged in a study of Balen et al. (2006), who investigated amphibolite, paragneiss, and mica schist samples from a locality in the Psunj metamorphic complex, ~ 3 km north of Kutjevo (Fig. 1). Only pre-Variscan monazite (420-460 Ma) was identified using electron microprobe-based monazite dating methodology. Based 1 Location studied by Balen et al. (2006) 2 Location studied by Horváth et al. (2010) Sample locations this study

Psunj metamorphic complex
Garnet-staurolite gneisses with amphibolite bodies Metagabbro 1 Location studied by Balen et al. (2006) 2 Location studied by Horváth et al. (2010) Jamičić et al. 1986;Jamičić and Brkić 1987;Jamičić 1988) with inset-map showing the position of the Slavonian Mountains (black) within the Tisia Unit and the Alpine-Carpathian-Dinaric tectonic framework, respectively. Inset map slightly modified after Schmid et al. (2008) on their data, Balen et al. (2006) suggested that the amphibolite facies metamorphism seen at this locality (~ 600 ℃, 10 kbar) has a late Ordovician or early Silurian age. Mica schists from another locality in the Psunj metamorphic complex near Gradište (Fig. 1) were studied by Horváth et al. (2010). In this study, however, only metamorphic monazite of Variscan (Early Carboniferous) age was found. Horváth et al. (2010) stated that the Variscan metamorphic evolution of this micaschist involved an early stage at ~ 630 ℃ and 7-8 kbar (recorded by garnet cores) and evolved clockwise into a penetrative crystallization stage at ~ 550 ℃ and 3-5 kbar recorded in the growth of staurolite and andalusite.

D I N A R I D E S S O U T H E R N A L P S
Considered together, the two local studies of Balen et al. (2006) and Horváth et al. (2010) lead to two important conclusions: 1) The Psunj metamorphic complex is obviously polymetamorphic with a Variscan and a pre-Variscan history.
2) The Variscan PT evolution appears to be different from the pre-Variscan one, in that it involves a higher thermal gradient.

Objectives of the present study
The first aim of our study is to better define the regional extent of pre-Variscan and Variscan metamorphism in the metapelite series of Krndija, using the method of in-situ electron microprobe (EMP) dating of monazite grains in thin section. For this purpose, about 20 samples of metapelitic rocks were collected (Fig. 1). A second objective is to confirm the age of the pre-Variscan metamorphism by an independent dating method to obtain more precise geochronometric information than is possible with the EMP monazite dating method. To this end, we carried out Lu-Hf isochron dating of garnet fractions separated from a chlorite-muscovite-garnet schist from north of Kutjevo. Third, geochemical data are presented, which provide new constraints on the protoliths of the investigated metasedimentary rocks, as well as on the palaeotectonic environment during deposition. Based on our new geochronological data and by carefully considering the petrography of the metamorphic rocks, we can finally provide new constraints on the (poly)metamorphic history of the Slavonian Mountains. Constraints on the palaeoplate-tectonic history are also discussed.

Th-U-Pb monazite dating with the electron microprobe
Monazite analyses were carried out with a wavelength dispersive system on a JEOL-JX8600 electron microprobe at the Paris Lodron University of Salzburg, following routines described in detail in Krenn et al. (2008) and Krenn and Finger (2004). The analytical procedure considers suggestions of Pyle et al. (2002) and Jercinovic and Williams (2005) of how to optimize the precision and accuracy of the method. Single point ages and errors were calculated after Montel et al. (1996). Weighted mean ages were calculated using the Isoplot program 3.41b (Ludwig 2007) and the software IsoplotR of Vermeesch (2018).

Lu-Hf dating of garnet
Ca. 2-3 kg of chlorite-muscovite-garnet schist with euhedral garnet were processed in a jaw crusher and sieved. Three garnet fractions, weighing between 70 and 400 mg, from the 125-250 µm sieved material were handpicked under a binocular microscope, and grains with visible inclusions were excluded. These handpicked fractions were then digested in closed Teflon ® vials on a 120 °C hotplate to selectively dissolve the garnet and avoid a contamination by potential zircon and rutile inclusions. The method used is described in detail by Lagos et al. (2007) and Zeh and Gerdes (2014). Note that high-pressure Parr ® bombs were not used. After rinsing with Milli-Q H 2 O, the garnet solutions were spiked with a mixed 176 Lu/ 180 Hf tracer for Lu and Hf concentration determinations and then digested as follows: the garnet was decomposed in HF-HNO 3 -HClO 4 and then 10 M HCl, and the samples were dried at high temperature (fuming HClO 4 ) between each step (Lagos et al. 2007). Separation of Lu and Hf was achieved using an ion-exchange column containing Eichrom Ln-Spec resin . Lu and Hf isotope ratios were measured using a Finnigan Neptune MC-ICP-MS equipped with a Cetac ARIDUS™ sample introduction system at Goethe University, Frankfurt. This instrumental setup ensured high sensitivity, enabling the precise measurement of Hf isotope compositions of samples having as little as 10 ng of Hf at a signal intensity of ~ 350 mV of 176 Hf for a 10 ppb Hf solution. Because only 50-80% of the Yb were separated from Lu during the purification technique employed here, a correction is necessary for the interference of 176 Yb on 176 Lu (Blichert-Toft et al. 1997). The Yb interference was monitored by measuring two interference-free Yb isotopes ( 173 Yb and 171 Yb). Their ratio was used to apply an instrumental mass bias correction to measured 176 Lu/ 175 Lu values, assuming a 173 Yb/ 171 Yb of 1.129197 and the exponential law ). This mass bias correction was also applied to the 176 Yb/ 171 Yb used for correcting the 176 Yb interference on 176 Lu. In-run statistics for the measured and corrected 176 Lu/ 175 Lu range from 0.006 to 0.01% 2 S.E. Mass bias on the Hf isotope ratios was corrected using 179 Hf/ 177 Hf = 0.7325 and the exponential 1 3 law. The Hf solutions were virtually free of any Yb and Lu but contained various amounts of Ta and W. All isobaric interferences on Hf isotopes were monitored and corrected using the mass bias corrected 173 Yb/ 176 Yb, 175 Lu/ 176 Yb, 180 Ta/ 181 Ta and 183 W/ 180 W values. During the course of this study, the Lu standard yielded 176 Lu/ 175 Lu = 0.026550 ± 17 (2σ ≈ 0.07%, n = 6) and the Hf Standard JMC-475 yielded 176 Hf/ 177 Hf = 0.282155 ± 12 (2σ ≈ 0.6ε, n = 6). In-run analytical uncertainty for 176 Hf/ 177 Hf was typically around ± 0.3ε or better. For the calculation of the mineral isochrons, the ISOPLOT program (Ludwig 2007) was used with a 176 Lu = 1.867 × 10 −11 a −1 Söderlund et al. 2004). For the 176 Hf/ 177 Hf uncertainties, we used quadratic additions of the reproducibility of the JMC-475 (2σ = 0.6ε) and the in-run precisions (2σ errors). Uncertainties on 176 Lu/ 177 Hf were propagated from the reproducibility of the Lu standard and the uncertainty in the spike calibration (0.15%) and multiplied by an error magnification factor that depends on the measured 176

Whole-rock geochemistry
Twenty X-ray fluorescence analyses were performed on glass beads (major elements) and pressed pellets (trace elements) using a Bruker S4 pioneer crystal spectrometer device equipped with a 4 kW rhodium tube at the Paris Lodron University of Salzburg. The analyses refer to ca. 1 kg sample material. Major elements were generally determined at reduced tube energies. Counting times were chosen, such that the relative 2σ uncertainties were better than 1% for SiO 2 and Al 2 O 3 , and better than 5% for elements occurring at the 1-10 wt. % concentration level. For the determination of trace elements, tube conditions and counting times were optimized automatically up to 4 kW and 400 s per element to obtain a detection limit of at least 3 ppm (3σ). Typical errors (2σ) from the counting statistics are 1-2 ppm at low concentrations (< 10 ppm), ~ 5 ppm at the 100 ppm concentration level, and better than 50 ppm at the 1000 ppm level. Calibration is based on a set of 30 international rock standards (mainly USGS standards).

Thin-section studies, mineral analyses, and pressure-temperature (PT) calculations
Thin sections were first studied carefully under the optical microscope. Mineralogical and microstructural details were further investigated using a scanning electron microscope (SEM: Zeiss Ultraplus) at the Paris Lodron University of Salzburg. Mineral analyses were partly performed with an Oxford X-max 50 detector coupled to the Zeiss Ultraplus, partly with the electron microprobe Jeol JX 8600 in wavelength dispersive mode. Multiple control analyses carried out repeatedly on the same minerals of the same samples showed that both methods yield consistent results.
PT calculations were performed using various methods. For most samples, peak PT conditions were estimated using conventional geothermobarometers for the mineral assemblage biotite-garnet-plagioclase-muscovite (calibrations of Battacharya et al. 1992;Holland and Powell 1990). For these calculations, we used garnet rim compositions (unmodified by retrograde diffusion) and nearby grains of plagioclase (rims), biotite, and muscovite. Reported PT variations result from the analyses of at least five different microdomains per sample (see Table 1). In addition to conventional geothermobarometry, average PT conditions were estimated with the THERMOCALC software 3.33 (Powell and Holland 1998), which is based on the internally consistent thermodynamic data set file "tc-ds55.txt" and on a-X relationships described in Holland and Powell (2003) and Coggon and Holland (2002).
To obtain PT path information during garnet core and rim formation, PT pseudosections were calculated in the system MnO-Na 2 O-CaO-K 2 O-FeO-MgO-A l 2 O 3 -SiO 2 -H 2 O-TiO 2 ( MnN CFMASHT) using the THERIAK-DOMINO software and the thermodynamic data set of Holland and Powell (1998) transcribed for, and implemented into, THERIAK-DOMINO by Gaidies et al. (2008). Further detailed information on calculation parameters and the used solid-solution models are given in Gaidies et al. (2008). PT data for different garnet core and rim domains were estimated by comparison between measured/observed and calculated mineral assemblages, modal and individual garnet compositions including mainly the grossular and pyrope components (see Figs. 4 and 5). The geochemical input data for THERIAK-DOMINO result from analyses of the ca. 2 × 3 × 1 cm rock off-cuts from the thin-section preparation (XRF analyzed) and are thus representative of the probed thin sections. H 2 O was assumed to be in excess, because the samples contain chlorite and are free of carbonate. Ferric iron has not been considered, because no free Table 1 Petrographic, thermobarometric, and geochronological data for the investigated samples 5.6 ± 1.0 351 ± 14 (n = 10/16) ferric iron phases like hematite or magnetite were found, and the hematite component in ilmenite is negligible. In addition, the Y-in-monazite thermometer of Heinrich et al. (1997) was used.

Petrographic data
The metapelitic rocks from the Krndija area can be divided into four petrographic groups. Each group is defined by a distinct mineral assemblage and crops out in a distinct regional zone (Fig. 2).
Group I metapelites comprise fine-grained phyllites from the northern part of the Psunj metamorphic complex. The regional distribution of these phyllites is well defined in the existing geological map (phyllite zone in Fig. 1). Parts of the phyllite zone are presently unaccessible due to landmines from the Yugoslavian war. We collected two samples (BS 502 and BS 606) of finely folded greenish-gray phyllites (Fig. 2). Both consist of ~ 40% quartz, ~ 20% chlorite, ~ 20% muscovite, and 10-20% albite. Accessory minerals are apatite, rutile, tourmaline, zircon, xenotime, and monazite. PT estimates carried out with THERMOCALC and the yttriumin-monazite thermometer (Heinrich et al. 1997) indicate lower-to-middle greenschist facies metamorphic conditions of about 400-450 °C (Table 1), consistent with the phyllitic rock texture (Fig. 3a) and the general absence of biotite and garnet.
Group II metapelites: This rock type outcrops along the road from Kutjevo to Orahovica (Figs. 1, 2). The finegrained, greenish-gray rocks contain small garnet porphyroblasts (3-4 vol. %) embedded in a chlorite-muscovite-quartz matrix (Fig. 3b). A few albite crystals and accessory epidote grains are present in the matrix as well, but there is no biotite. Further accessory minerals are rutile, ilmenite, apatite, and zircon. Rutile grains are commonly mantled by ilmenite.
From our sampling (Fig. 2), it would appear that group II metapelites form a ~ 300 m-thick sheet in the footwall of the phyllite zone. However, these rocks are not delineated as a discrete lithological unit in official maps and it is, therefore, unclear whether they extend along the entire southern margin of the phyllite zone or constitute only a local lens. Pamić and Lanphere (1991) briefly referred to this rock type as a lower-T equivalent of the Grt-Bt-Ms schists, which are exposed further to the south (our group III metapelites), and transitional in metamorphic grade to the phyllites to the north.
The garnet in group II metapelites is almandine-rich (Alm = 55-63 mol%). It shows a weak zoning with an inner domain that is slightly lower in calcium (Grs = 13 mol%) than the rim domain (Grs = 17 mol.%; Fig. 4a; Table 2). The spessartine component decreases relatively steadily from core to rim (25-13 mol%), whereas the almandine and pyrope components increase, as does the Fe/Mg. The observed zonation patterns indicate prograde garnet growth. Results of geothermobarometry and PT pseudosection modeling define a nearly isothermal pressure increase during garnet growth from 530 ℃/7.5 bar to 540-570 ℃ at 9.0 kbar (cf. Figure 4a, 5a, Table 1).
Group III metapelites are fine-to medium-grained mica schists to gneisses, which typically contain the assemblage Grt-Bt-Ms-Chl-Pl-Qtz-Rt, with quartz and plagioclase being volumetrically dominant (Table 1). In rare places, minor amounts of staurolite are also present. Group III mica schists occur along a ca. 25 km-long and 2 km-wide zone north of Kutjevo (Fig. 2). They grade into the biotite-free group II metapelites toward the north. The sampling site of Balen et al. (2006) refers to a locality in the middle part of this group III metapelite zone.
The abundance of biotite (3-19 vol.%) and the ubiquitous (2-10 vol.%), strongly zoned garnet prophyroblasts with Fig. 2 Geological map of the Krndija region (eastern part of the Psunj metamorphic complex) showing the regional distribution of group I, II, III, and IV metapelites as defined in this study. Rock legend after Jamičić and Brkić (1987). Note that group II, III, and IV metapelites were previously mapped as one unit. The general dip of the rocks is toward the north. The phyllite zone occupies the hanging wall position of the complex 1 3 euhedral grossular-poor cores (Grs = 7-9 mol.%) surrounded by grossular-rich rims (Grs = 20-25 mol.%) distinguish group III metapelites from group II metapelites. Muscovite contents vary mostly between 3 and 20 vol.% and chlorite contents between 1 and 10 vol.%. Euhedral rutile crystals are a characteristic accessory mineral in group III metapelites and are up to 500 µm in size. Four ca. 100 µm-sized staurolite grains are identified in two samples (BS 419b, 420). They have ilmenite inclusions(!) and are marginally resorbed and surrounded by newly formed muscovite and chlorite. In addition, a tiny 20 µm staurolite inclusion was identified in a Grs-poor garnet core.
The garnet crystals of group III metapelites are commonly euhedral (Fig. 3) with no substantial replacement by mica or chlorite obvious. Microtextures suggest that the garnet rims grew in equilibrium with most of the matrix minerals (except staurolite). The typically euhedral, angular, and unresorbed, Grs-poor garnet cores are of variable size and, in places, occupy more than two-thirds of the garnet porphyroblast (Fig. 3). Some of these cores contain small ilmenite inclusions. No rutile inclusions were observed. In many cases, the cores are internally weakly zoned with grossular and spessartine contents decreasing from the center outwards (Fig. 4b), while the pyrope component increases. A severe chemical discontinuity is everywhere observed at the core-rim contact; grossular contents suddenly increase four-or fivefold, whereas both pyrope and almandine components significantly decrease. Interestingly, there is no great change in the spessartine content between core and rim, implying that the older garnet generation was not subjected to significant resorption before the grossularrich rim started to grow.
The observed garnet zonation patterns, together with the results of geothermobarometry and PT pseudosection calculations, point to a fairly complex PT history of the host rock. Good agreement between measured and modeled garnet compositions is obtained for the low-Ca garnet cores by assuming a nearly isobaric temperature increase from 575 to 600 °C at 7 kbar (event M1). The zonation pattern of the Ca-rich garnet rim (event M2) is consistent with a clockwise PT evolution from 570 °C/10 kbar over 575 °C/11 kbar to 585 °C/10 kbar (thermal peak). The peak PT conditions at the end of the clockwise PT loop overlap, within error, with those obtained by the conventional geothermobarometry for unretrogressed group III metapelites (570-620 ℃/9-11 kbar, see Table 1). These PT data are in good agreement with those published by Balen et al. (2006). It could be argued that the samples from the north of the "Kutjevo" zone ( Fig. 2), collected close to the group II metapelites (BS 412-416), record slightly lower temperatures than samples Pseudosections for group II, group III and group IV metapelites (samples BS 418, BS 203a, BS 411c). Black arrows refer to P-T data derived from garnet zonation and observed mineral reactions (c: garnet cores, cf. Figure 4). Two possible P-T evolutions 1 and 2 are shown in b and discussed in the text. c P-T evolution of group IV metapelites during the Variscan. The garnet cores in these rocks may be of Ordovician age and likely traveled an independent P-T loop. Thin dashed arrows: P-T paths suggested by Horváth et al. (2010) for group III and IV metapelites  from the south; however, the difference is not considered enough to be conclusive. The derived PT path sections for group III metapelites imply that staurolite + ilmenite (observed in 2 samples) were stable phases during formation of the low-Ca garnet cores but became unstable during formation of the Ca-rich garnet rims (see Fig. 5b). This interpretation is consistent with the above-mentioned observation that staurolite has only an isolated occurrence in the rocks, is everywhere partially replaced by muscovite and chlorite, and, in one case, is preserved as an inclusion in a garnet core.
It is important to note that there is not much retrograde alteration of most of the group III samples. Garnet porphyroblasts are typically perfectly euhedral and biotite is only rarely chloritized. This indicates that rock exhumation was under (most likely dry) conditions that did not favor retrograde mineral reactions and fast enough to prevent significant volume diffusion (see also the discussion in Horváth et al. 2010). The only evidence for a retrograde reaction in group III metapelites is the partial replacement of rutile by ilmenite, which must have happened during decompression.
Significant retrogression of group III metapelites is only observed along discrete shear zones/bands, which cut the main foliation at an oblique angle (samples BS 423 and 505 in Table 1). These shear bands are commonly filled with synkinematic chlorite, muscovite, quartz, and albite, and biotite is absent. Garnet porphyroblasts are commonly replaced by chlorite. Retrograde deformation also produced chains of small, aligned garnet clasts within the shear bands.
Group IV metapelites: Mica schist samples collected in the southern part of the Krndija area are macroscopically not much different from group III metapelites (note the similar modal contents in Table 1). Nevertheless, these rocks carry a completely different type of garnet with low grossular content. Importantly, rutile is widely absent and ilmenite is abundant in these rocks, which sets them apart from group III metapelites. A further petrographic difference to group III metapelites is the widespread presence of small amounts (< 1 vol.% in most samples) of late-stage andalusite porphyroblasts intergrown with large chlorite crystals. In general, group IV metapelites have a finer grain size than group III metapelites (Fig. 3), and the peak mineral assemblage exhibits a comparably much stronger degree of retrogression. Garnet is commonly surrounded by a biotite-chlorite reaction seam. Staurolite, where present, occurs exclusively as small grains in the matrix and nowhere is it included in garnet. Hence, it is interpreted to be younger than the garnet.
The more common, core-free garnet is typically weakly zoned with a slight increase of the pyrope component, and  Wimmenauer (1984), and in CaO vs SiO 2 and Rb/Sr vs SiO 2 diagrams. Major elements are normalized to 100% before being plotted. Metapelites of group I to IV have different protolith compositions, implying that they belong to different stratigraphic series. The Radlovac Formation, which is shown for comparison, comprises sedimentary rocks with significantly higher Rb/Sr ratios decrease in the spessartine and grossular components from the center to the rim of the grains. A slight increase of the spessartine component at the edges of these garnet grains is considered a chemical modification due to back-diffusion. The rims around the high Ca cores show comparable chemical zonation patterns as in the core-free garnet grains (Fig. 4).
Results of geothermobarometry and PT pseudosection modeling suggest that the rare Ca-rich cores formed at PT conditions of 580 °C/11 kbar, and the dominant low-Ca garnet porphyroblasts (or overgrowths) formed during a clockwise PT evolution from 530 °C/ 5 kbar to 570 °C/6 kbar (Fig. 5). The widespread growth of tiny staurolite and andalusite grains in the matrix is evidence that the metamorphic crystallization continued along a retrograde path from the staurolite into the andalusite field, whereby both garnet and staurolite became resorbed and replaced by chlorite and mica. This alteration very likely caused the observed backdiffusion of Mn into the outermost garnet rims. Due to continuous recrystallization of the rocks during decompression, results of conventional thermobarometry must be viewed with caution. Nevertheless, the obtained PT data (Table 1) are consistent with those obtained from the garnet zonation patterns (Fig. 5).

Geochemical composition of the metasedimentary rocks of the Krndija region
The SiO 2 /Al 2 O 3 vs K 2 O/Na 2 O discrimination diagram of Wimmenauer (1984) suggests mainly pelitic protoliths for group II metapelites and greywacke protoliths for groups I, III and IV metapelites (Fig. 6). However, notable chemical differences are discernible between group III and group IV metapelites in terms of CaO contents and Rb/Sr ratios (Fig. 6, Table 3). CaO and Sr contents are significantly higher in group III metapelites and are taken to reflect a higher proportion of detrital plagioclase in the sedimentary protoliths.
For comparison, we have also plotted data from metasedimentary rocks in the Radlovac formation (Biševac et al. 2009(Biševac et al. , 2013Balen et al. 2013). These are generally clay-rich, feldspar-poor pelites (Fig. 6). Group I metapelites, unfortunately represented by only two samples, are Ca poor similar to the Radlovac metasedimentary rocks but lack the very high Rb/Sr ratios that characterize the latter.
The diagram of Bhatia and Crook (1986) indicates that most of the precursors to the Krndija metapelites were deposited in an active plate margin environment or, more specifically, in a continental island arc environment (Fig. 7).

Th-U-Pb electron microprobe dating of monazite
Group I metapelites: The two samples taken from the phyllite zone contain several, although generally small, monazite grains. Measured Th-U-Pb dates (Table 5) are combined to yield weighted mean ages of 91 ± 31 Ma (BS 502) and 97 ± 26 Ma (BS 606), constraining monazite crystallization to the Cretaceous. Yttrium contents in the monazite grains are generally low (< 0.5 wt.% Y 2 O 3 ). Application of Y-inmonazite thermometry (Heinrich et al 1997) indicates that the monazite formed at temperatures of around 400-450 °C, consistent with the other thermobarometric estimates for the rocks (Table 1). The Cretaceous monazite ages are interpreted to date the timing of low-grade metamorphism of the phyllite zone.
Group II metapelites: Although five samples have been carefully screened, monazite was not detected in group II metapelites.
Group III metapelites: Accessory monazite is common in these samples and grains from seven samples were analyzed. The grains are typically 5-20 µm in size, round-tosubhedral, and are evenly disseminated throughout the rock matrix. Some of the monazite grains are enclosed in garnet rim zones. Calculated single-point dates vary from 263 (± 123) Ma to 631 (± 129) Ma, but by far most dates are between 430 and 540 Ma with associated individual errors (2σ) typically of ± 70-90 Ma (Table 5). An exception is the mylonitic sample BS 423, where monazite ages below 400 Ma are frequently found (see below). Because the single-point dates generally overlap within error, it is interpreted that the analyses define one age population in each sample.
Calculated weighted mean ages for single samples also overlap within error and are between 465 ± 34 Ma (sample BS 203b) and 481 ± 26 Ma (sample BS 413a). The overall mean of all samples is 471 ± 15 Ma and is interpreted to constrain monazite growth in group III metapelites to the middle Ordovician. Notably, these monazite ages are significantly older than those published by Balen et al. (2006) for their locality from the same zone (444 ± 19 and 428 ± 25 Ma).
Xenotime is generally absent in group III metapelites (again with the exception of sample BS 423). The strongly varying and commonly moderately low Y contents in the monazite analyses indicate a general yttrium undersaturation. Applying Y-in-monazite thermometry, monazite grains with the highest measured Y contents (Y 2 O 3 = 2 wt.%) yield crystallization temperatures of around 600 °C, similar to those obtained using other geothermobarometric methods for group III metapelitic rocks (see above).
The monazite grains in sample BS 423 (mylonitic shear band) mostly yield Variscan ages (weighted mean of 372 ± 26 Ma). In some monazite grains from this sample, remnant Ordovician cores were detected (Fig. 8d, Table 5). The Variscan monazite grains in sample BS 423 all have low Y 2 O 3 contents (< 0.8 wt.%). Given the fact that xenotime is also present in the sample (Fig. 8d), these low Y contents in the Variscan monazite are consistent with growth at lowtemperature conditions of around 400-450 °C. Interestingly, the monazite grains in sample BS 423 show apatite-epidote coronas indicating a late-stage, post-crystallization monazite consumption, similar to that described by Finger et al. (1998).
Group IV metapelitic rocks: Accessory monazite was found in two samples (BS 408a, BS 411c) and is of similar size and shape to that in group III metapelites. However, in group IV samples, many monazite grains are overgrown by an apatite-allanite corona. Most monazite crystals occur in the matrix, and some are enclosed in low-Ca garnet. There is a large variation in Th contents with values ranging from only ~ 1 wt.% (and extremely large age errors) up to about 17 wt.% ThO 2 ( Table 5) . The measured Th-U-Pb ages vary between 274 and 447 Ma, but by far most dates are between 300 and 400 Ma with individual errors of ~ 30-80 Ma (2σ). Almost identical weighted mean ages of 354 ± 16 Ma and 351 ± 14 Ma are calculated for the two investigated samples (Table 1). These age data are interpreted to record a single Variscan event, although it cannot be excluded that some older monazite relics remain unrecognized due to large analytical errors. The ages presented here overlap the monazite age obtained by Horváth et al. (2010) from the Gradište locality (356 ± 23 Ma). We interpret that the Variscan monazite crystallized (or recrystallized from a pre-existing monazite) close to the thermal peak of the Variscan PT evolution, and slow cooling in the andalusite stability field led to the formation of the secondary apatite/allanite coronas around the monazites (Finger et al. 1998).

Lu-Hf isotope analyses and garnet dating
Three garnet fractions and one whole-rock sample from the group II metapelitic rocks (sample BS 417) were analyzed for Lu-Hf isotopes ( Table 4). The garnets are slightly zoned. However, the zonation is continuous and is interpreted to record one prograde growth phase between 7 kbar/830 ℃ and 9 kbar/850 ℃ (see Figs. 4a, 5a). All four analyses plot on an isochron ( Fig. 9) with an age of 466.0 ± 2.3 Ma (MSWD = 1.8) and the latter is accepted as geologically meaningful (for more details, see discussion below). The whole-rock analysis yields εHf 466Ma = − 4.7 ± 0.7 (Table 4), indicating that the protolith of the metapelitic rock must have contained relatively old crustal components. This is reflected in a two-stage hafnium model age of 1.62 Ga. Fig. 7 Th-Sc-Zr/10 ternary diagram after Bhatia and Crook (1986) showing data from metapelites in the Kutjevo Zone. A-Oceanic island arc; B-Continental island arc; C-Active continental margin; D-Passive margin

The polymetamorphic architecture of the Krndija region
Contrary to the frequently cited model of a monocyclic Variscan metamorphic evolution (Pamić and Lanphere 1991), it now becomes increasingly obvious that the Krndija area in the Slavonian Mountains records a complex, polymetamorphic history involving Sardic (Ordovician), Variscan, and Alpine metamorphic imprints. The EMP monazite age dates from this study not only highlight but place some regional constraints on this polymetamorphic evolution. We find the following monazite age distribution in the Krndija region (Figs. 2, 10):  1 3 1) Alpine monazite ages occur in the northern part of Krndija, in the so-called phyllite zone (i.e., group I metapelites). 2) Variscan monazite ages are recorded in group IV metapelites on the southern flank of the Krndija metamorphic complex that define a zone (herein termed the Gradište Zone) of penetrative, amphibolite facies grade, Variscan metamorphism. 3) Ordovician monazite ages are prominent in the intervening crustal strip (i.e., group III metapelites).
Notably, despite its only moderate precision, EMP-based monazite dating has proved to be a key method in revealing the complex geological history of the study area. Admittedly, it would appear at first sight that up to half of the measured single-point monazite ages are insignificant due to very high individual age uncertainties. However, even without any data screening, the overall statistics of the measured monazite dates fundamentally records a polymetamorphic evolution with Alpine, Variscan, and Ordovician age peaks (Fig. 10a). The geochronological constraints become tighter when veryhigh-error dates are omitted (Fig. 10b) and weighted average ages for single samples are considered (Fig. 10c, Table 5). However, as a precaution, the geological significance of such statistically optimized EMP monazite dates should always be confirmed by another geochronological method. In this study, a precise Lu-Hf garnet age of 466 ± 2.3 Ma was obtained for the (monazite-free) metapelites of group II, which reside between group III metapelites and the phyllite zone. This age is consistent with the average monazite age of 471 ± 15 Ma for the adjacent group III metapelites.
Based on all the geochronological information, we propose that a large geologically discrete composite zone, with well-preserved remnant Sardic metamorphic basement, is present in the central part of Krndija and has a minimum length of 23 km and is approximately 2 km wide (Fig. 2). Small amphibolite and metagranitoid lenses intercalated in group III metapelites (Balen et al. 2015) are most probably  Table 5). See text for further explanation also pieces of Sardic metamorphic basement. This composite Sardic complex in the Krndija zone, which we name the Kutjevo Zone, complements the Strona-Ceneri Zone (Zurbriggen et al. 1997), and should be recognized as a key region in any study of the "Sardic orogeny".

The Sardic record of the Kutjevo Zone
Our data support the pioneering conclusions of Balen et al. (2006) that the Slavonian Mountains host a pre-Variscan metamorphic component. However, the monazite ages published in Balen et al. (2006) were rather in favor of a late Ordovician or Silurian metamorphic event at about 440 Ma (although with very large error), whereas our data clearly support a middle Ordovician metamorphic event, which coincides with the so-called Sardic phase described in the other parts of the Alps (see Introduction section). The peak T stage of the Sardic metamorphism in the Kutjevo Zone, which is represented by growth of garnet in group II and garnet rim zones in group III metapelites, took place at 540-620 °C and at relatively high-pressures of 9-12 kbar typical of a subduction zone setting. Based on the Lu-Hf garnet isochron age for sample BS 417, we propose that this Barrovian metamorphic event occurred around 466 Ma. Note that the dated garnets from this sample can be expected to provide a geologically meaningful Lu-Hf isochron age, because they are continuously zoned and likely formed one-phase prograde between 7 kbar/830 °C and 9 kbar/850 °C (see Fig. 4a, 5a).
A previous (slightly or perhaps even significantly older) metamorphic event (M1) with a different PT evolution is Weighted average ages (bold) are at the 95% confidence level 1 3 recorded in garnet cores in the group III metapelites. The composition of these cores reveals a pre-HP stage of prograde isobaric heating at 7 kbar (Figs. 4b, 5b). The typically discontinuously zoned garnets experienced a growth break when passing through the 7-10 kbar pressure interval. Their rim zones record the Sardic high-pressure event in form of a clockwise PT evolution to peak T conditions of 580-620 °C and 9-10 kbar (Figs. 4b, 5b). Metastable staurolite, which is locally found in group III metapelites, is considered an M1 relict. Horváth et al. (2010), who studied a sample of group III metapelites (also with low-Ca garnet cores), proposed a simple prograde PT evolution from 550 ℃ at 6 kbar to 640 ℃/12 kbar for this rock (thin dashed arrow in Fig. 5). However, such a PT trajectory is clearly not consistent with the observed M1 garnet core zoning (which shows no compression) and, is indeed inconsistent with the presence of early staurolite in the rock. Furthermore, a continuous zoning in Xgrs, Xsps, and Xpy would be expected in garnet growth during a simple prograde PT path but, instead, the garnet porphyroblasts in group III metapelites show an abrupt compositional change at the core-rim transition.
The age of the M1 metamorphism is so far unconstrained. It could be pre-Sardic (for instance of Cadomian age), in which case an exhumation stage between the M1 and the M2 event may be a realistic assumption (path 2 in Fig. 5b). On the other hand, the M1 event could also be explained in terms of a mid-crustal contact metamorphism caused by early Sardic magmatism (perhaps represented by amphibolite and orthogneiss lenses in the metapelites). In the latter case, an anticlockwise compressional PT evolution (path 1 in Fig. 5b) is very feasible. The observation that the M1 garnet cores are widely unresorbed (euhedral) may lend some support to the latter model, but any further discussion on this matter must await robust and precise geochronological data on the formation ages of the surrounding amphibolites and orthogneisses, which are as yet not available.
The accurate time difference between M1 and M2 could perhaps be resolved in the future by attempting to directly date the M1 staurolite using U-Pb methods or by attempting to carefully separate the Ca-poor garnet cores for Lu-Hf dating. Also, the monazite in group III metapelites could record both events. Considering the work of Spear (2010), we assume that metasedimentary rocks like group III metapelites, with low-to-moderate CaO contents < 1.8 wt. % (Fig. 6), have a wide monazite stability field, which spans the PT conditions of both M1 and M2. Thus, a first monazite growth event could theoretically have occurred during M1 and a second growth or recrystallization event during M2. However, based on the observation that only a few monazite grains are enclosed in the garnet (and only in garnet rim zones), we conclude that the main phase of monazite growth in group III metapelites took place close to the T climax of the M2 event, i.e., after the M2 peak-P stage. The calculated weighted average ages (Fig. 10c) should, therefore, be tentatively interpreted as such, i.e., until further data are available, for example from more precise (e.g., SIMS-based) in-situ dating of monazite.
Another question that requires clarification in a future study concerns the apparent lithological division of the Kutjevo Zone into chlorite-rich (biotite-free) group II metapelites and biotite-bearing group III metapelites. This is unlikely to be just the result of variations in metamorphic grade, as stated by Pamić and Lanphere (1991). Our study shows that there is no pronounced hiatus in the metamorphic grade: the peak metamorphic conditions recorded in the biotite-free, chlorite-bearing group II metapelites are only slightly below those recorded in the biotite-bearing group III metapelites (540-570 °C/8-9 kbar vs. 570-620 °C/8-12 kbar). It is proposed here that the given petrographic contrasts (chlorite vs biotite assemblages) result mainly from bulk rock chemical differences (Fig. 6, Table 3), i.e., higher Al-contents and lower X Fe in group II metapelites have stabilized chlorite instead of biotite. Obviously, we are dealing with two different sedimentary formations on the active Gondwana margin, as such, determination of the deposition ages via detrital zircon dating would be worthwhile.

The Variscan event
Traces of the Variscan orogeny are mainly recorded in the Gradište Zone. Our petrological data suggest that this zone underwent a clockwise PT evolution at medium-to-lower grade amphibolite facies conditions during the early Carboniferous (at ca. 355 Ma). The Variscan PT loop is characterized by a prograde PT increase from 530 ℃/ 5 kbar to 570 ℃/ 6 kbar (metamorphic peak), followed by a PT decrease to 550 ℃ at 3 kbar (Table 1).
Rare Ca-rich cores in some garnet porphyroblasts of group IV metapelitic rocks (Figs. 4c, 5c) represent probable relics of a pre-Variscan higher-P metamorphism. Considering the PT conditions, we assume that these garnet cores are Sardic in age. It is thus proposed that the Gradište Zone originally belonged, at least in part, to the Sardic metamorphic complex.
An interesting question arises as to how the Variscan orogeny affected the Kutjevo Zone. As stated earlier, monazite grains in lower-to-middle greenschist facies mylonitic rocks of the Kutjevo Zone record Variscan ages (sample BS423). This implies that at least parts of the Kutjevo Zone experienced low-grade metamorphic conditions during the Variscan orogeny (Table 1). However, the Variscan reheating obviously did not result in an appreciable recrystallization of the Kutjevo Zone rocks, as can be inferred from the low degree of retrogression of the Sardic high-P assemblage and the preservation of Ordovician monazite ages.
Another interesting question is whether the phyllite zone experienced Variscan metamorphism. High-grade metamorphic mineral relics are not observed in these rocks, so any metamorphic overprint must have been low grade. If the protoliths of the phyllites were indeed Upper Devonian sedimentary rocks, as proposed by Jamičić (1983) and Jamičić and Brkić (1987), then a continuous Variscan metamorphic zonation could exist in Krndija from amphibolite facies grade in the south (Gradište Zone) over low grade (greenschist facies) in the Kutjevo Zone to very low grade in the northern Phyllite Zone. On the other hand, it is also possible that the protoliths of the phyllites were deposited in post-Variscan times; for example, during the Upper Carboniferous, Permian or Mesozoic. Dating of detrital zircons would be necessary to solve this question.

The Alpine event
Results of monazite dating and petrographic observations indicate that the phyllites in northern Krndija crystallized in the Cretaceous under low-to-middle greenschist facies conditions (Table 1).
We speculate that the older, higher grade basement units of Krndija (the Kutjevo and Gradište Zones) also experienced a low-grade heating event during the Cretaceous, although this is difficult to confirm using the current geochronometric methods. Local chlorite-bearing shear zones with an S2 foliation in the basement rocks (Balen et al. 2006) could well have formed during the Alpine event (but alternatively also during the Variscan). Nevertheless, the generally unaltered nature of the Sardic mineral parageneses in the Kutjevo zone (group II and III metapelites) allows for the conclusion that Alpine metamorphism only barely affected the crystalline basement of Krndija.

Comparison of the Kutjevo Zone with the Strona-Ceneri Zone and palaeogeographic implications
It is undeniable that there are striking petrographical parallels between the Strona-Ceneri Zone and the Kutjevo Zone. First, there is the lithological association, which is in both cases clearly dominated by metapelitic (metagreywacke) rocks with a pre-middle Ordovician formation age. These metasediments are inter-layered with numerous small amphibolitic units as well as metagranitoids. According to Zurbriggen et al. (2017), the latter may represent syntectonic I-type granitoids of Ordovician age.
In addition, both areas received their main deformation during an Ordovician Barrow-type metamorphic event and lack signs of significant post-orogenic (Ordovician) reequilibration at P < 7 kbar. Ordovician eclogite facies metamorphism as recorded in the Strona-Ceneri Zone (Franz and Romer 2007) is so far not found in Krndija. However, it should be mentioned that the amphibolitic rocks of Krndija, which are considered the most likely to contain eclogite relics, are still not well investigated. Garnet core compositions in the Kutjevo Zone metapelites may be interpreted in terms of syn-compressional contact metamorphism and syntectonic magmatism, respectively (see above). This would match the tectonothermal situation that Zurbriggen (2017) has proposed for the Strona-Ceneri Zone. Zurbriggen et al. (1997Zurbriggen et al. ( , 2017 interpreted the Strona-Ceneri Zone as a Cambro-Ordovician fore-arc complex that was subducted in the middle Ordovician. A similar scenario is also feasible for the Kutjevo Zone, because (i) the metapelitic rocks in the Kutjevo Zone are also associated with amphibolites with MORB signatures, and (ii) they were, according to geochemical data, most likely deposited along an active continental margin (Fig. 7).
There is wide agreement today that the pre-Variscan intra-Alpine basement units are mainly derived from an eastern segment of the northern Gondwana margin. They may have resided either adjacent to South China (von Raumer et al. 2013) or in front of the Arabian-Nubian Shield and the Sahara craton (Siegesmund et al. 2021;Finger and Riegler 2022). Also, there is wide agreement that large parts of the northern Gondwana margin had an Andean-type character during the Cadomian orogeny (e.g., Nance et al. 2008;von Raumer et al. 2013). However, it remains controversial whether and to what extent this Cadomian Andean-type setting persisted into the Ordovician. Many geologists believe that subduction had widely ceased by the late Cambrian and that most of the northern Gondwana margin became a passive margin in the Ordovician (see, e.g., Žák and Sláma 2018 and further references therein). Indeed, evidence for Ordovician subduction-related metamorphism is fairly rare in central Europe. Until now, such evidence was reported only from a few intra-Alpine regions with the Strona-Ceneri Zone as the most prominent example (Faryad et al. 2002;Biino 1995;Zurbriggen et al. 1997;Handy et al. 1999).
We can now present further data in favor of a Sardic subduction scenario from the Slavonian Mountains in Croatia. This region is almost 1000 km away from the Strona-Ceneri Zone. However, considering the palaeogeographic reconstructions of Stampfli et al. (2013), both units were situated relatively closely together in Ordovician times (Fig. 11). In the palaeogeographic map for the middle Ordovician (461 Ma) by Stampfli et al. (2013), the positions of both units are located somewhat remote from the Gondwana plate edge, to the east behind a large Ordovician back-arc rift and adjacent to South China (Fig. 11). However, prior to the back-arc rifting, they could have been positioned much closer to the Gondwana plate margin (note arrows in Fig. 11) and an Alaskan-type subduction-accretion scenario as proposed by Zurbriggen (2015) is absolutely feasible.
It is important to note that Stampfli et al. (2013) support the notion that an extended Ordovician Cordillera did exist along the northern Gondwana margin (Fig. 11). They suggest that much of this Cordillera was later detached from Gondwana along with the Hunic ribbon terrane, and finally became part of Asia (Fig. 11). Based on this model, the intra-Alpine basement terranes with Sardic metamorphism could be interpreted as Cordillera relics that remained connected to Gondwana throughout the Early Palaeozoic and later became part of the Variscan collision zone.
However, the Hunia terrane model as such can be fundamentally challenged, since several distal (i.e., ocean near) segments of the Cadomian orogen are distributed over both the Alps (Siegesmund et al. 2021;Neubauer et al. 2022) and the Bohemian Massif (Žák and Sláma 2018), a fact which would argue against a Cordillera loss. Apart from that, the concept of a persistent Cadomian-Ordovician subduction zone along the entire north Gondwana margin ) is disclaimed by numerous geologists, who argue instead in favor of a predominantly passive margin setting in the Ordovician (e.g., Franke et al. 2017). Indeed, regions with Sardic orogenic imprints outside of the Alps (e.g., Alvaro et al. 2016;Casas 2010;Cocco and Funedda 2017) do not show evidence for subduction and could be explained in terms of transpressional-transtensional tectonics along the north Gondwana margin as well.
A viable compromise that we suggest here is to invoke an (Alaskan-type?) Ordovician subduction setting only for a distinct eastern segment of the north Gondwana margin, which we may term E-Armorica following Finger and Riegler (2022) and which is now mainly exposed in the Alps (Fig. 11). In this scenario, the more western parts of frontal north Gondwana (W-Armorica according to Finger and Riegler 2022, see Fig. 11) had already undergone a change from an active to a transtensional-transpressional passive margin setting by the end of the Cambrian.

Conclusions
1) The Krndija region of the Slavonian Mountains, Croatia, was affected by at least three major tectonometamorphic imprints: during the Middle Ordovician (Sardic event), the Early Carboniferous (Variscan event), and the Cretaceous (Alpine event).

3
2) Based on monazite dating and thermobarometric studies on metapelitic rocks, three tectonometamorphic zones can be delineated in Krndija: the Kutjevo Zone with well-preserved medium-grade Ordovician basement, the Gradište Zone with penetrative medium-grade Variscan metamorphism, and the Phyllite Zone with low-grade Alpine metamorphism.
3) The Sardic metamorphic event in the Kutjevo Zone is additionally confirmed by a precise Lu-Hf garnetwhole-rock isochron age of 466.0 ± 2.3 Ma. 4) The Kutjevo Zone is recognized as a key region for the study of the Sardic (Ordovician) orogeny. It shows striking similarities to the Strona-Ceneri Zone in the western Alps in terms of lithology and metamorphic history. 5) We suggest a tectonic scenario where the E-Armorican (now intra-Alpine) sector of the northern Gondwana margin experienced a subduction setting in the middle Ordovician, while the W-Armorican segment (now exposed in the French, German and Czech Variscides) had already mutated from a (Cadomian) subduction setting to an extensional (transtensional-transpressional) setting by the late Cambrian.