The Karavansalija Mineralized Center at the Rogozna Mountains in SW Serbia: Magma Evolution and Time Relationship of Intrusive Events and Skarn Au ± Cu–Pb–Zn Mineralization

Abstract

The Karavansalija Mineralized Center (KMC) with its Au–Cu skarn mineralization associated with the Rogozna Mountains magmatic suite in southwestern Serbia belongs to the Oligocene Serbo-Macedonian magmatic and metallogenic belt (SMM-MB). Samples from intrusive and volcanic rocks at the KMC show typical arc signatures of subduction-derived magmas through enrichment in large-ion lithophile elements (LILE) and depletion of high–field strength elements (HFSE). The magmas developed a high-K (calc-alkaline) fractionation trend and evolved toward shoshonitic compositions. Whole-rock trace element data suggest plagioclase-absent, high-pressure amphibole ± garnet fractionation that generates adakite-like hydrous magmas during evolution in lower crustal magma chambers. Zircon LA–ICP–MS and high-precision CA–ID–TIMS dating together with zircon trace elements and Hf isotope measurements were carried out in order to couple the geochronologic and geochemical evolution of the KMC. The results suggest that magmatism starts around 29.34 Ma with granitic to rhyodacitic subvolcanic intrusions followed by a more evolved magmatic intrusion that was emplaced into Cretaceous limestone, generating a widespread skarn alteration at ca. 28.96 Ma. After a period of quiescence of about 1.2 My, either another magma body evolved or the same upper crustal magma chamber was recharged and also likely partly reactivated older plutonic rocks as indicated by xenocrysts. The REE ratios shift from apatite, titanite ± amphibole-dominated fractionation of the older magmatic event to crystallization of allanite, efficiently depleting the LREE and Th/U in the younger upper crustal magma. After a lamproite-like melt was injected, the increased heat and fluid pressure led to the expulsion of a quartz-monzonite porphyritic stock at ca. 27.72 Ma, strongly interacting with the skarns and established a fertile hydrothermal system. Soon after a non-mineralized second pulse of some porphyry dykes cut the previous phenocryst-rich “crowded” porphyries and skarns at ca. 27.60 Ma, thus bracketing the maximum timespan of ore mineralization to about 112 ± 45 Ka. Increased contribution of a lamproite-like melt is inferred from the presence of phlogopite micro-phenocrysts, phlogopitization of biotite, and diopside clusters in the latest porphyry dykes. There is a trend of increased crustal assimilation from the oldest volcanic phase to the emplacement of the youngest porphyry dykes recorded by ɛ-Hf of the zircons. Oligocene occurrences of significant base metal mineralization within Serbia, northern Macedonia, and Greece, e.g., Crnac, Rudnik, Veliki Majdan, Stratoniu, or the Cu–Au porphyry at Buchim (northern Macedonia), are all associated with trachy-andesitic (quartz latitic) porphyry dykes, which originated through post-collisional tectonic settings or upper plate extension involving reworking of crustal arc-derived rocks and partial melting of the mantle wedge. This study demonstrates that on the basis of field relationships and the application of high-precision CA-ID-TIMS zircon age data, pulses of porphyry dykes of a 10ka age range can be distinguished, and the timing of mineralization can be parenthized.

Introduction

Zircon geochemistry and geochronology are widely applied to improve our understanding of crustal magmatic systems and constrain timescales of geological processes, such as those forming events of economically important magmatic–hydrothermal Cu–Au–Mo deposits.

Metal deposits (Pb, Zn, Cu, and Au, among others) provide successive samples of magma extracted from crustal-scale hydrous magmatic systems, skarn reaction between granitic magmas and marbles including fluid-enriched solutions, or/and remobilization of the base metals out of the skarn into new-growing granodioritic dikes. They are, therefore, a critical source of information about the processes and rates of magma ascent, magma storage, and fluid generation, bridging those of volcanism and pluton formation. The identification of the processes that lead to skarn deposit formation (Meinert et al., 2005) and later to enrichment of “porphyry-like” intrusion dikes and the timescales on which they operate can provide us with valuable information about magmatic processes but could also potentially help in discriminating fertile magmatic systems from ubiquitous infertile systems, resulting in barren intrusions or volcanic eruptions.

Recent advances in high-precision zircon geochronology by chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA–ID–TIMS (Mattinson, 2005; Bowring et al., 2011; McLean et al., 2011; Condon et al., 2015; McLean et al., 2015) now allow dating of multiple intrusions associated with ore formation with unprecedented precision. The improved precision permits the constraint of rapid events, such as individual skarn and hydrothermal mineralization phases (<100 kyr; von Quadt et al., 2011; Buret et al., 2016; Tapster et al., 2016), that typically occur during a long term period of volcanism and intrusive magma emplacement extending over several million years (Halter et al., 2004; Rohrlach and Loucks, 2005; Maksaev et al., 2006; Lee et al., 2017). The integration of the chemical information gained from zircon can yield time-calibrated information about magma chemistry, thermal evolution, and mixing processes during zircon crystallization in magmatic systems (Chelle-Michou et al., 2014; Samperton et al., 2015; Buret et al., 2016; Szymanowski et al., 2017).

Geological timescales for magmatic and hydrothermal processes involved in ore formation have been suggested based on in situ LA–ICP–MS U–Pb data (Banik et al., 2017; Lee et al., 2017; Garwin, 2000) and CA–ID–TIMS geochronology, which became increasingly precise (von Quadt et al., 2011; Lehmann et al., 2013; Chelle-Michou et al., 2014; Tapster et al., 2016; Large et al., 2018). However, several studies applying multiple techniques to the same sample sets have resulted in variable dates (von Quadt et al., 2011; Chiaradia et al., 2013; Chiaradia et al., 2014; Chelle-Michou et al., 2014). The discrepancy demands a more detailed understanding of the precision and accuracy of the techniques and statistical data treatment that are applied to derive a geological age.

For the present study, we obtained a large dataset of zircon geochemistry and geochronology by laser ablation inductively coupled plasma mass spectrometry (LA–ICP–MS) followed by high-precision geochronology of the same zircon crystals, utilizing the CA–ID–TIMS technique. This study focuses on a regional suite of predominantly calc-alkaline granitoid rocks, which crop out along the border between the easternmost Vardar zone and Serbo-Macedonian Massif (Figure 1). We focus on the Karavansalija Mineralized Center (KMC) that lies within the base metal–dominated area of the Serbo-Macedonian magmatic and metallogenic belt (SMMMB) surrounded by the Cenozoic magmatic intrusives and extrusives of the Rogozna mountains (Serbia; Figures 1–3). The KMC was mentioned as an occurrence of Au–Cu mineralization with secondary base-metal enrichments, which is distal to known deposits in the Rogozna Mountains (Urosević, M. et al., 1973; Kandić, M. et al., 1987). The ore mineralization is described by Budinov et al. (2015) as hosted mainly by garnet–pyroxene–epidote skarns and the nearby quartz–epidote rocks. However, the KMC stands out due to its Au-bearing skarn mineralization. In particular, we shall discuss the much debated question as to whether the Oligocene–Miocene intrusions of the SMMMB can be followed into the Rhodope Massif and to the Cassandra ore district and whether these intrusions represent the back-arc magmatism with respect to the occurring Pb–Zn–Au deposits, or alternatively, whether they are related to the post-orogenic collapse of an overthickened crust in the internal Dinarides (Koroneos et al., 2010), mantle delamination, and/or to the extension in the Pannonian basin (e.g., Cvetković et al., 2007).

FIGURE 1

FIGURE 2

FIGURE 3

These coupled data from the KMC base-metal deposit allow us to resolve the chemical evolution and the changing physical state of the magma reservoir over time and the timescales of hydrothermal processes.

Regional Geology

The Western Tethyan magmatic belt, whose remnants occur today in SE Europe, can be traced through the countries of Slovakia, Hungary, Romania, Serbia, Kosovo, northern Macedonia, Greece, and Bulgaria (Figure 1). The magmatic belt is known as one of the Europe’s most significant ore-bearing belt including Cu–Au and Pb–Zn–Ag metal deposits. The tectonic setting and geodynamic evolution along the 3,500-km strike of the Western Tethyan belt are complex and broadly evolved from the Upper Jurassic to the Late Cretaceous (Gallhofer et al., 2015). Along the Sava suture zone, subduction-related arc magmatism including locally extensional-related magmatism was fairly minor until the Eocene (Ustaszewski et al., 2010; Stojadinovic et al., 2013; Gelder et al., 2015; Stojadinovic et al., 2017; Toljić et al., 2018); the reactivating of magmatism during the early Eocene (<42 Ma) probably started with the transitioning from collision to post-orogenic extension (late Eocene to early Oligocene) and a further widespread extension-related magmatism during the Miocene (Seghedi et al., 2004; Richards, 2015). This resulted in the formation of the late Eocene to Oligocene sedimentary basins, volcanic centers of calc-alkaline, high-K calc-alkaline, and shoshonitic composition, and extension-related epithermal deposits that formed between 30 and 35 Ma (Moritz et al., 2014, Figure 2). The Late Cretaceous Apuseni–Banat–Timok–Srednogorie belt (ABTS belt) (von Quadt et al., 2005; Gallhofer et al., 2015) is dominated of Cu-rich deposits, with a total endowment of ∼20 million tons (Mt) of Cu compared to ∼15 Mt of Cu deposits in the Cenozoic, whereas the Cenozoic porphyry deposits contain the largest concentration of the Au content (63.6 Moz) (Baker, 2019). Furthermore, there are a few major economic high-sulfidation epithermal deposits directly associated with porphyry deposits in the Cenozoic. The Cenozoic Au ± Cu deposits span greater ranges in age, distribution, and style than the Cretaceous deposits (Baker, 2019; Marchev and Singer 2002; Marchev et al., 2010). In a very simplistic view, the age variation of the Cenozoic deposits show several clusters at 35–31 Ma, 28–24 Ma, 15–12 Ma, and <9 Ma (Baker, 2019) in comparison to the Cretaceous age cluster with its Cu–Au deposits between 92 and 80 Ma (von Quadt et al., 2005; Jelenković et al., 2016; Banješević et al., 2019).

The central Balkan Peninsula is characterized by a complex geotectonic framework which mainly developed during Mesozoic convergence and collision processes. There is a general agreement that the Tethys finally closed during the late Mesozoic and that the West Dinarides and the East Vardar units represent relicts of former oceans. The Serbo-Macedonian magmatic and metallogenic belt (SMMMB, Figure 1) represents one major part, traditionally been considered as the Pb–Zn–Bi–Ag–Au dominant ore district in Cenozoic time (Janković, 1997; Borojević Šoštarić et al., 2012). The occurring Cenozoic granitoids associated with ophiolites provide important constrains on the geodynamic setting. The Late Eocene to earliest Miocene magmatic activity is widespread in the Dinarides and is thought to be related to rollback and possibly slab break-off of the Adriatic/African plate (Figure 1) (Andrić et al., 2018). The Eocene to Oligocene core complex formation as observed in the Kopaonik area is abundant in the Dinarides and well documented (Table 1) for the northern Bosnia, Cer, Bukulja, Fruška Gora, and Jastrebac areas (Ustaszewski et al., 2010; Matenco and Radivojević, 2012; Stojadinovic et al., 2013; Toljić et al., 2013; Mladenović et al., 2015; Erak et al., 2017) and the Central Rhodopes (Marchev et al., 2005; Rohrmeier et al., 2013). Published geochronological data (Table 1) of the Rogozna Mt. magmatic rocks allow for distinguishing two main magmatic events: 1) the andesitic to dacitic volcanics occur in the eastern part of the Rogozna Mountain magmatic complex with Ar/Ar ages of ∼30 Ma and 2) the quartz-latites lasted from 29.5 ± 0.1 to 27.3 ± 0.1 Ma (Borojević Šoštarić et al., 2012). The process of roll-back of the subducted oceanic crust in the Carpathians and the Aegean area (Rhodope Massif) induced local clockwise rotation of the Dinarides and caused the sharp bend of the Sava suture zone south of Belgrade in map view (Figure 2) as documented by paleomagnetic studies (Lesić et al., 2019).

This region has been traditionally considered as a Pb–Zn–Ag dominant district (Janković, 1997; Borojević Šoštarić et al., 2012), but in fact, it has a large Au endowment (∼25 Moz; (Baker, 2019), as well. Gold is associated with a number of Oligocene–Miocene porphyries (Skouries, Illovitza, Buchim, and Tulare) and carbonate replacement deposits (Olympias and Piavitsa), as well as with the Plavica high-sulfidation epithermal deposit. Notable are the Trepča Pb–Zn deposit in Kosovo with still 60 million tons of ore reserves left; the Lece district in southeastern Serbia with its epithermal gold mineralization (Pamic and Balen, 2001); the Perama Hill high sulfidation Au–Ag–Te–Se deposit in northeastern Greece hosted in altered volcanic rocks; and Skouries, the high grade Au–Cu porphyry deposit in northern Greece (Melfos and Voudouris, 2017; Voudouris et al., 2019).

Geology of the Karavansalija Mineralized Center (KMC)

The KMC is subdivided into three mineralized areas called Shanac, Gradina, and Copper Canyon (Figures 3, 4). The most dominant lithology on the surface is represented by volcaniclastic rocks that correspond to the earlier stage of the Cenozoic magmatism. Shallow-intrusive bodies and dykes of probably the Cenozoic age intrude the volcanic and sedimentary rocks. The skarns are mainly exposed in the southern part (Figures 3, 4) of the KMC, and they are described as products of metasomatic alteration (Urosević et al., 1973). The Upper Cretaceous limestones reacted with the intruding of shoshonitic extrusive rocks during the last stage of the Oligocene–Miocene volcanism, and they are composed of epidote, grossular, zoisite, calcite, pyroxene (hedenbergite), and rarely quartz. Proximal garnet, diopside, epidote ± wollastonite, followed by a pyroxene–epidote zone with subordinate garnet skarn mineral composition and distal skarnized limestone with dominant epidote and occasional pyroxene are described by Dusanić (2006; citation in Budinov et al., 2015).

FIGURE 4

The KMC was mentioned as an occurrence of gold–copper mineralization with secondary base-metal enrichments, which is distal to known deposits in the Rogozna Mountains (Urosević, et al., 1973; Kandic et al., 1987). The ore mineralization is described by Budinov et al. (2015). It is hosted mainly by garnet–pyroxene–epidote skarns and the nearby quartz–epidote rocks and the overlying volcanic and volcaniclastic rocks. Disseminated cobalt–nickel sulfides from the gersdorffite–krutovite mineral series, cobaltite, and pyrite–marcasite–chalcopyrite–base metal aggregates are defined in the epidosites. The skarn sulfide mineralization is characterized by chalcopyrite, pyrite, pyrrhotite, bismuth-phases (bismuthinite and cosalite), arsenopyrite, gersdorffite, and sphalerite (Budinov et al., 2015). Arsenopyrite aggregates with subordinate amounts of pyrite and chalcopyrite have been described in the overlying volcaniclastic rocks. Gold is rarely developed in visible aggregates of native gold but occurs mainly “invisible” in arsenopyrite. Fluid inclusion micro-thermometry with homogenization temperature in the range of roughly 150–400°C and varying salinities (0.5–8.5 wt% NaCl_eq for two-phase low-density fluid inclusions and 15–41 wt% NaCl _eq for two- and three-phase high-salinity fluid inclusion) are interpreted as evidence for two fluid sources of the Cu–Au mineralization—magmatic and meteoric. The primary magmatic source of mineralizing fluids is confirmed by sulfur isotope data ranging from 0.4‰ to +3.9‰ ^δ34S_CDT in arsenopyrite, chalcopyrite, pyrite, and marcasite from the different sulfide assemblages (Budinov et al., 2015).

Relationships of Skarn, Hydrothermal Alteration, and Ore Mineralization

The KMC Cu–Au skarn deposit is a result of the multiphase hydrothermal activity caused by discrete magmatic pulses. The reaction products of skarnitization mainly occur in the Copper Canyon (Figures 3, 4); they are observed within the contact zone between the granitic porphyry (e.g., EOKSC 1677) and the Cretaceous limestone (Figures 5A, D, E, F). Endoskarns and exoskarns are distinguished at the KMC (in the porphyries and the limestones, respectively). Garnet is the dominant rock-forming constituent of the extensive exoskarns, appearing as relatively small (≤1 mm) grossular crystals. Insights from SEM-EDS (Burkhard, 2017) suggest that the primary garnet belongs to the grossular–andradite solid solution. It may be observed as the “atoll garnet” (Galuskina et al., 2007; Taghipour et al., 2013)—turbid to almost opaque cores of garnet grains, due to fracturing and replacement, with rather clear unaffected rims (Figures 5A, C). Common associates of the grossular garnet in exoskarns are grains of hedenbergite, clinopyroxene, and epidote. Corrosion of the garnets and formation of chlorite and/or epidote highlight the retrograde overprint of the garnet skarn. Moreover, calcite is abundant as an interstitial phase between garnet grains but also in fractures through the garnets. In some samples, e.g., EOKSC 1361A, 182 m (Figure 5E , Figure 6C), in a strongly base-metal–mineralized dark-green to black rock mass, ilvaite was found as a dominating alteration mineral. The latter is considered a retrograde skarn mineral replacing garnet and clinopyroxene (Burt, 1989). Fe-rich axinite occurred as a vein constituent along with quartz and calcite (EOKSC 1361A, 51 and 164 m; Figure 5B). The first intrusion of the granitic porphyry shows itself a phyllic overprinting, i.e., feldspars, amphiboles, and biotite have all been altered to sericite and minor quartz (observed at EOKSC 1677, 104 m). This hydrothermal alteration along with the extensive skarnitization field, prograde, and retrograde mineral associations provides evidence for high syn-to post-intrusion fluid activity which would be required for extensive metasomatism.

FIGURE 5

FIGURE 6

Intrusions of the “crowded” quartz-monzodioritic porphyry (CP; description in next chapter 3) discretely cut through the exoskarn, proving that CP is younger than the latter. However, the porphyry is itself prone to strong skarn reactions along fractures/veins. Although the contacts of this porphyry to the exoskarns are typically clear, intrusive juxtaposition of those chemically very different lithologies did not leave the intruding porphyry unaffected (Figures 5C, D, F). Bleaching (Figure 5B) associated with weak (i.e., incipient) skarn reactions is observed. Endoskarn veins contain garnet, but the abundance of epidote and, especially, chlorite and calcite points to a generally retrograde skarn characteristic.

In endoskarns, incipient reactions in porphyries could be recognized as the matrix taking a turbid greenish hue in haloes around thin (≤1 cm) skarn veins, or along contacts to the exoskarn. The coloring effect is accompanied by mineral reactions by which plagioclase is replaced by alkali feldspar and amphiboles and biotite are transformed into calcite, minor clinopyroxene, and alkali feldspar, whereas epidote and calcite precipitate in small amounts along fractures.

The overall volume of CP intrusions remains unspecified; however, they clearly occur as dykes. Thus, they are expected to lack the necessary heat for significantly extending the existing exoskarns field (related to the earlier emplacement of the granitic porphyry). On the other hand, the frequent occurrence of endoskarn sections does prove a hydrothermal event induced by the CP intrusion. All lines of evidence point to the CP intrusion being the pulse that delivers and precipitates ore minerals as the later intrusion of porphyry dykes may reveal weak endoskarn veining but does not cause significant skarnitization (Figures 5A, C–F) neither host economic mineralization.

The main ore minerals at the KMC are pyrite, arsenopyrite, pyrrhotite, chalcopyrite, sphalerite, and galena. They are commonly found in exoskarn and endoskarn samples. Magnetite is notably present in EOKSC 1361A (182 m) and 0,831 (340 and 508 m). Hematite was found only in accessory amounts, and there were occurrences of marcasite together with pyrite. Pyrite and arsenopyrite seem to be early phases; they occur along the boundaries of skarn rocks to veins, or alternatively, they have sub-to euhedral shapes in skarn masses. Sphalerite and galena likely appear the latest in the mineralization sequence; pyrrhotite and chalcopyrite are, thus, thought to form before them. However, these sulfide ore phases are often intimately associated with each other, which preclude the identification of a clear precipitation sequence.

The skarn formation with its propylitic and phyllic alteration and the contact relationships between skarns and the different porphyries, skarn veins and mineralization indicate three main phases of hydrothermal activity (Hoerler, J. 2017): 1) formation of extensive prograde exoskarns (garnet skarn) in the sedimentary host rocks, emplacement of a comparably larger subvolcanic granitic intrusion; 2) cooling of the skarn system that leads to retrograde overprinting of the garnet skarn; and 3) intrusion of the quartz-monzonitic porphyry dykes that partly reactivates the skarn system and results in vein-bound mineralization of base metals and, especially Au in exoskarn and endoskarn veins (Burkhard, 2017). Subsequent quartz-monzonitic porphyry intrusion remains unmineralized, but it does also show less pronounced endoskarn veining.

Analytical Methods, Sample Location, and Petrographic Description

The details of the methods and sample location are documented in the Supplementary Figure S1. All selected samples are listed in Table 2, including the number of the borehole, the depth, the QAP classification, and the rock type; the selected samples are shown on a geological map, and the relationship to host rocks is documented in two geological profiles of the KMC target area (Figure 4B).

The porphyritic rocks have been split in types according to the visible texture and also their relative timing (Figure 5): a more phenocryst-rich (“crowded”) quartz-monzonitic porphyry (CP) variety (Figures 5B, D) and a quartz-monzonitic porphyry with less phenocrysts (“proper”) in a fine-grained matrix (PP; Figures 5C, E, F). The phenocrystal content of the quartz-monzonitic porphyries ranges between 45–60% and 20–33%, respectively. Almost all quartz-monzonitic rocks are strongly skarn-altered. The volcanic samples (EOKSC 1687, 17103) consist mostly of feldspar (sanidine), amphibole, biotite, and quartz as finer grains instead of larger and resorbed phenocrysts. The texture is fine-grained, and some flow bands are recognized in the hand specimen.

The porphyritic rocks (CP) EOKSC 1361A and 0,723 and the porphyritic rocks (PP) 1361B, 0,831, and 1,681 (Table 2) are classified as quartz-monzonitic, quartz-monzodioritic, and granitic (QAP classification) and appear as dykes of several meters to 100 m in thickness. The main phenocrysts of the PP rocks are quartz, plagioclase, K-feldspar, amphibole, biotite, and sanidine. The PP-type intrusions are less encountered in the drill-cores and show a smaller extent than the CP (Figure 5). They are as well strongly altered and bleached (Figure 5B). The PP dykes contain phenocrysts of quartz, plagioclase, K-feldspar (sanidine), amphibole, diopside, biotite, and phlogopite; they also contain small diopside clusters and evidence of phlogopitization of biotite and diopside selvages around quartz. The granitic porphyry EOKSC 0727 (LP) is the least altered sample compared to all andesitic and granitic/rhyodacitic porphyries marked consequently as late porphyries-LP, Figures 9, 10; cross-cutting dykes are often strongly bleached in the typical white appearance. The skarnitization occurrence is documented in Figures 4–6 between the intrusion contact zone and the marble units; the volcanic rocks on top of the marble and the skarn area in the Shanac area (EOKSC 1677, 1,687, 17103, and 1,259) have a strong sericite–pyrite alteration and replacement of the feldspar into a fine-grained mass of sericite, therefore interpreted as predating the hydrothermal alteration and mineralization.

The knowledge of the geological and temporal relationship of the Eldorado Gold Corporation from field and drill cores and our observations on the relationships of the magmatic rocks to the skarns, hydrothermal alteration, and ore vein were used to select the samples for precise CA–ID–TIMS dating: a phenocryst-rich quartz-monzonitic porphyry (EOKSC 1361A; CP), two apparently younger quartz-monzonitic/granitic porphyries (EOKSC 1361B and 0,831; PPs), a granite porphyry dyke (EOKSC 0727), and a subvolcanic (shallow intrusive body) granite (EOKSC 1677).

Results

Geochemistry of Magmatic Rocks

On the total alkalis vs. silica diagram, the quartz-monzonitic porphyries, quartzz-monzodioritic, and granite/rhyodacite plot in the trachy-andesite field, whereas some granitic rocks are shifted to the andesitic area (Figure 7A, Supplementary Table S3). All extrusive samples (volcanics) show a trachy-andesitic composition (Figure 7A). Available geochemical data of magmatic rocks which are in close relationships to ore deposits, e.g., Buchim (North Macedonia), Rudnik, and Veliki Majdan (Serbia), are of broader compositions, ranging from basaltic-trachy-andesitic to trachy-andesitic and further to andesitic composition (Figure 7A); the K₂O vs. SiO₂ plot shows that most samples are of shoshonitic to high-K composition (Figure 7B).

FIGURE 7

To avoid the strong influence of hydrothermal alteration on the rock classification, the discrimination diagram of Winchester and Floyd (1977; Figure 7C) was used as it takes the immobile trace elements into account. The majority of intrusive and volcanic rocks of the KMC area plots in the trachy-andesitic, trachytic, and rhyodacitic field, whereas some samples plot into the phonolitic field. However, the latter samples still contain quartz and, therefore, should not be defined as phonolites; they are similar to rocks from Veliki Majdan that represent hybrids of mixing between acid and lamproitic magmas (Prelevic´ et al., 2004). Additional data from other locations with magmatic–hydrothermal ore deposits (Buchim, Rudnik, and Veliki Majdan) from the Cenozoic SMM-MB are plotted on Figure 7 for comparison and listed in Table 1. The results of our study mostly coincide with the ones from the Buchim porphyry (Cu ± Au; (Lehmann et al., 2013), with a slight tendency to a more alkaline, phonolitic composition. Rocks from the KMC area are more silicic and less potassium-rich in contrast to the rocks from Rudnik and Veliki Majdan which mainly plot in the shoshonite series field.

The rare Earth element (REE) patterns show enrichment of light REE (LREE) and depletion of the heavy REE (HREE) with a weak Eu anomaly (Figure 8A); the multi-element diagram of the KMC rocks reveals high–field-strength element (HFSE) depletion for Ta, Nb, and Ti, whereas Ba, Sr, and Pb and other large-ion lithophile elements (LILEs), especially U and Th are enriched (Figure 8B). These patterns are generally indicative for subduction-related magmatism or for such that is generated in the subduction-modified lithosphere. The plots of La/Sm vs. Yb and Sr/Y vs. Y aim to distinguish between adakite-like and normal arc magmas (Figures 8C,D); adakite-like signatures for the ratios La/Sm and Sr/Y higher than 20 and Yb and Y concentration less than 1.9 and 18 ppm are used according to Defant and Drummond (1993) and Richards and Kerrich (2007). The porphyritic dykes of the KMC area (Figure 7C) plot on the border between normal arc- and adakite-like affinities. The analyses of the rocks from the Rudnik and Veliki Majdan area mostly plot in the adakite-like field. The intrusive rocks of Buchim and Rogozna show a trend between adakitic to normal arc on higher La/Sm and Sr/Y values. The (La/Sm)_n and (Dy/Yb)_n ratios (Figures 8E, F) are used to identify mineral differentiation paths (Davidson et al., 2007; Richards and Kerrich, 2007). The (La/Sm)_N ratio of all samples generally increases with the increasing SiO₂ content (amphibole–garnet fractionation) but shows a slight trend of possible monazite, allanite, or apatite fractionation at SiO₂ of 60% and higher. The (Dy/Yb)_N ratio decreases with increasing SiO₂, which implies amphibole fractionation in all KMC samples. Both diagrams with KMC data show a similar geochemical evolution to those from the Rudnik and Veliki Majdan areas, except the (Dy/Yb)_N rise above 62 wt% SiO₂ in the Rudnik samples, indicating garnet fractionation.

FIGURE 8

LA–ICP–MS and CA–ID–TIMS Zircon U/Pb Geochronology

The U/Pb LA–ICP–MS data, including ²⁰⁶Pb/²³⁸U, ²⁰⁷Pb/²³⁵U, and ²⁰⁷Pb/²⁰⁶Pb U/Pb ratios, are listed in Supplementary Table S4, and the obtained data sets are show in Figure 9.

FIGURE 9

The analyzed samples of the KMC show heterogenous zircon population, indicating the existence of inherited zircon population with Cretaceous, Triassic, and Silurian ages (Supplementary Table S4); the obtained pre-Cenozoic ages are not plotted but listed in Supplementary Table S4. The weighted mean dates of concordant zircons of the LA–ICP–MS study that could be geologically meaningful are plotted in Figure 9. All obtained and calculated ²⁰⁶Pb/²³⁸U ages of the LA–ICP–MS measurements are overlapping within their errors. The calculated 2σ errors of the weighted ²⁰⁶Pb/²³⁸U ratios are partly less than 0.8%; taking into account the uncertainty of the standard calibration, instability of the laser, and the estimated external error, approximal uncertainties show an approximation uncertainty of 1.5%, which is calculated for each data set (Supplementary Appendix Figure S2). For the full data tables, the readers are referred to the Supplementary Appendix Figure S2. The presumably older volcanics (17,103, 1,687, 1,677, and 1,259) cannot be distinguished simply based on the LA–ICP–MS data set (Figure 11).

The sub-volcanic rocks (Figure 9) have an age of 28.56 ± 0.56 Ma and 28.45 ± 0.54 Ma based on two samples (EOKSC 17103 and EOKSC 1687). Two granitic/rhyodacitic porphyries (EOKSC 1677 and EOKSC 1259) which belong to the older geological unit (Phase I, Figures 9, 10) have weighted mean ages of 28.16 ± 0.53 Ma and 28.49 ± 0.53 Ma. Two pre-ore quartz-monzonitic/granitic dykes (Figures 9, 10) yield ages of 27.35 ± 0.59 Ma and 27.38 ± 0.41 Ma [(Phase II), EOKSC 1361A and EOKSC 0723]. The syn-ore granitic, quartz-monzodioritic, and quartz-monzonitic dykes have ages of 27.47 ± 0.60 Ma, 27.37 ± 0.41 Ma, and 27.42 ± 0.59 Ma based on samples (EOKSC 0831, EOKSC 1681, and EOKSC 1361B, Figures 9, 10), respectively.

One post-ore granitic dyke (PDMC 0727-LP) has a weighted mean age of 26.95 ± 0.40 Ma (MSWD = 1.13). The obtained ²⁰⁶Pb/²³⁸U ages of the older magmatic group (Phase I, Figures 9, 10) show an overlapping within their uncertainty with the data of the syn- and post-ore formation (Figure 10). The obtained dates of Phase I (Figure 11), Phase II, ore-formation, and post-ore formation are overlapping although the Phase I with their older mean values indicate a small offset to a previous event.

FIGURE 10

FIGURE 11

High-precision CA–ID–TIMS data (Supplementary Table S5,) help to resolve the time relationship of rocks and distinguish the presumably older volcanics rhyodacitic/granitic intrusions (Phase I) from the younger quartz-monzodioritic and granitic dykes (phase II) (Figures 9, 10). Zircons of all intrusive and extrusive rocks likely record continuous zircon growth in an evolving magma chamber (Bachmann et al., 2007; Barboni et al., 2013; Buret et al., 2016; Keller et al., 2018). The use of the youngest zircon to determine the age of eruption or magma emplacement has been tested by several studies and accepted in the research approach (von Quadt et al., 2011; Schoene et al., 2013; Wotzlaw et al., 2013; Condon et al., 2015; Buret et al., 2016; Large et al., 2018). This approach of course could potentially introduce uncorrected errors (e.g., through zones of Pb-loss on a single zircon grain); this effect, however, should be minimal as the studied zircons are young, and the chemical abrasion technique should be highly efficient.

Several zircon grains of the rhyodacitic-granitic rock (EOKSC 17103, and 1,677) have been selected for TIMS analysis. The youngest zircon age of every sample is regarded as the porphyry emplacement age (Wotzlaw et al., 2013; Szymanowski et al., 2017; Keller et al., 2018). The weighted mean LA–ICP–MS ages and the CA–ID–TIMS dates are listed in Table 3; further details are available in Figures 9, 10 and Supplementary Table S4, and Supplementary Table S5. The CA–ID–TIMS dates of the granite and volcanics are notably older than the calculated LA–ICP–MS average mean dates (Table 3). The discrepancy of U/Pb zircon ages obtained by the CA–ID–TIMS and LA–ICP–MS techniques has been observed at several geological locations (von Quadt et al., 2014; Large et al., 2020; Rottier et al., 2020).

Phase I: The youngest zircon grain (EOKSC 17103) records a date of 29.343 ± 0.03 Ma, and the youngest zircon grain (EOKSC 1677) records a date of 28.964 ± 0.01 Ma (CA-ID-TIMS, Figure 10). Phase II: The youngest statistically equivalent zircon grains (n = 2) yield a weighted mean age of 27.759 ± 0.013 Ma (MSWD = 3.4). Ore formation (bracket between the altered and mineralized CP quartz-monzonitic dyke EOKSC 1361A and the PP quartz-monzonitic dykes EOKSC 0831 and 1361B): eight, six, and seven zircon grains of the quartz-monzonitic dykes (EOKSC 1361A), (EOKSC 1361B), and (EOKSC 0831) have been selected for TIMS analysis (Figure 10); the youngest zircon records a date of 27.607 ± 0.02 Ma (EOKSC 0831), and the youngest zircons record dates of 27.759 ± 0.013 Ma (EOKSC 1361B), respectively. Post-ore formation: The post-ore granitic porphyry (EOKSC 0727-LP) shows based on the youngest grain an age of 27.583 ± 0.031 Ma (Figure 10).

Zircon Trace Element Budget and Hf Isotopes

The analyses of zircon trace element data can be used for distinguishing magmatic pulses on their geochemical characteristics, determine the crystallization temperature, and identify the role of assimilation and fractionation components (Ferry and Watson, 2007; Wotzlaw et al., 2013; Szymanowski et al., 2017).

The REE patterns in zircon are characterized with LREE depletion with respect to HREE together with a weak negative Eu and a strong positive Ce anomaly (Figure 11; calculation, Supplementary Figure S1). The REE content on zircon rims yields a higher positive Ce anomaly and a slightly weaker Eu anomaly (Eu* = 0.38–0.47) than that on the zircon cores (Eu* = 0.41–0.51; Figure 11B). The REE content variations between the samples are not visible or at most seen as slightly higher LREE values (excluding Ce) for the porphyry stock and dykes than all other samples (Hörler, 2017).

Zircon trace element analyses of the granitic, quartz-monzonitic, quartz-monzo-dioritic, and rhyodacitic intrusive/extrusive rocks (EOKSC 727, 1,259, 1,677, 1,687, 17103, 0,831, 1,681, 723, 1361A, and 1361B) have been plotted on Figures 12, 13; the data set is listed in Supplementary Table S7. We have used the Th/U, Ce/Sm, Yb/Gd, Ce, Yb, Ti ppm, and Eu-anomaly to demonstrate possible fractionation trends (Halter et al., 2004; Reid et al., 2011; Wotzlaw et al., 2013), mixing mineral fractionation scenario, and crystallization temperature. Several minerals, e.g., apatite, titanite, and allanite have a strong influence by the REE residual melt concentration (Klimm et al., 2008) on the fractionation process in magmatic systems (Szymanowski et al., 2017).

FIGURE 12

FIGURE 13

Ti-in-zircon temperatures were calculated using the calibration of Ferry and Watson (2007). For αSiO₂, a value of 1 was used as in all rocks, the quartz occurs as phenocrysts. Ti-bearing minerals are rarely seen in the fresh rock samples, but the trace element patterns propose the presence of titanite as it is also seen in other porphyritic rocks (e.g., in Argentina, Peru, South America) (Chelle-Michou et al., 2014; Buret et al., 2016); consequently, an αTiO₂ of 0.7 was used. Estimated temperatures (zircon rims) for quartz-monzonite, quartz-monzodiorite, and for the granite are low, ranging from 640°C up to about 700°C (Figures 12G,H); the estimated zircon core temperature, ranging from 650°C up to 750°C (Figures 12G,H). The granite 0,727 (Figure 12E) shows a narrow T field between 650 and 670°C, and all rim analyses (EOKSC 0727-LP) show T data below 650°C. A change in the TiO₂ or SiO₂ activity of 0.1 will result in minor differences of ±25°C (including the models uncertainty). These values plot near and probably also below the zircon solidus, similar to other studies on porphyritic rocks (Buret et al., 2016; Chelle-Michou et al., 2014). The role of plagioclase fractionation is visible between rocks of older and younger magmatic evolution; the Eu anomaly (0.4–0.6) of the quartz-monzonite and other sub-volcanic rocks drops down to values <0.4 of the post-ore–related rocks like quartz-monzonite and granite suites caused by a strong plagioclase fractionation (Table 3). The Th/U ratios show a decreasing trend from 0.6 to 0.8 of the volcanic/sub-volcanic rocks to values <0.4 of the post-ore granite, quartz-monzonite, and porphyries of the Buchim porphyry, as well (Lehmann et al., 2013).

The average Hf isotopic compositions of all magmatic rocks are plotted together with those from the Osogovo intrusion (Grozdev et al., 2017); the time-corrected data range from ɛ-Hf of 2.43 to −3.35. The older magmatic suite (Figure 13, Phase I) shows positive ɛ-Hf between 2.43 and 1.44, whereas granitic, quartz-monzonitic, and quartz-monzo-dioritic rocks have ɛ-Hf between −1.29 and −3.35, which demonstrates the descending trend from older to younger magmatic rocks (Figure 13; Supplementary Table S8). The data of sub-/volcanic rocks are positive in ɛ-Hf; all intrusive rocks are similar in the ɛ-Hf range, and the older intrusion (29.5 Ma) of Osogovo (Grozdev et al., 2017) has a stronger crustal input with an ɛ-Hf values at −7.

Discussion

Different Timescales of Processes Related to Magmatic Evolution of the KMC, the Mineralization, and Skarnitization Events

The CA–ID–TIMS data of the older magmatic group (Figure 10), the granitic porphyry (EOKSC 1677; age of 28.964 ± 0.01 Ma), and the granite (EOKSC 17103; age of 29.343 ± 0.03 Ma) yield ²⁰⁶Pb/²³⁸U ages that do not overlap within error uncertainties with the data of the porphyries of the syn- and post-ore formation (Figure 10). The next magmatic event (Phase II, Figure 10) starts with the CP quartz-monzonitic porphyry intrusion (EOKSC 1361A) with an emplacement age of 27.733 ± 0.03 Ma; the individual magmatic lifetime of this porphyry (EOKSC 1361A) with a range of 414 ± 45 ka is similar to that of Phase I (Figure 10). The time gap (1.231 Ma, Figure 10) between the quartz-monzonite EOKSC 1361A (Phase II) and the older magmatic group (Phase I, Figure 10) is remarkable and probably points to two independent geological events. The intrusion age of the quartz-monzonite porphyry (EOKSC 1361A, CP) defines the start of the mineralization at the KMC (Figure 10, 15). The quartz-monzonitic porphyries (EOKSC 1361B, and 0,831, PP) crosscut the porphyry EOKSC 1361A and mark the end of the mineralization event at 27.607 ± 0.02 Ma (Figure 10). The difference in age (112 ± 45 ka) between both, the youngest zircon grains of the EOKSC 1361A and the EOKSC 0831 rock, defines the lifetime of the ore formation process. The last and the youngest magmatic phase of the granitic rock EOKSC 0727 (LP) shows an emplacement age of 27.583 ± 0.031 Ma, which overlaps within the uncertainty with two members of the ore formation phase (EKOSC 1361B and EKOSC 0831, Figure 10).

The obtained lifetime of the ore formation at the KMC with 112 ± 45 ka (Figure 10) is only a maximum duration because ore mineralization could be on much shorter timescales (Buret et al., 2016; Large et al., 2021). These numbers are consistent with models discussed by Weis et al. (2012) where timescales of porphyry ore formation within a porphyry stock are calculated between 50 and 100 Ka. To date, no clear relationship between the duration of magmatic–hydrothermal activity and the size of epithermal/hydrothermal deposits was identified from studies applying high-precision CA–ID–TIMS geochronology. Comparison of published datasets (Buret et al., 2016; Tapster et al., 2016; Large et al., 2018; Large et al., 2021) reveals maximum durations of metal-forming events of between tens of thousands of years and 10⁵ yr. Although these studies are so far constrained to deposits of <10 Mt of contained Cu, these cover deposits that range over at least one order of magnitude in size (Koloula vs. Batu Hijau) (Tapster et al., 2016; Large et al., 2020). A correlation between the duration of the mineralizing event and the total mass of the deposited metal had been previously suggested based on compilations of different geochronological data sets (Lehmann et al., 2013; Chelle-Michou et al., 2017; Chiaradia, 2020). Longer durations of ore formation (> 1 Myr) were suggested based on Re–Os geochronology on molybdenite at the giant porphyry deposits and deposit clusters in Chile (> 50 Gt Cu; El Teniente), (Cannell et al., 2005 and Maksaev et al., 2004; Rio Blanco, Deckart et al., 2012) and Chuquicamata, (Barra et al., 2013). Mineralizing timescales of copper (and gold) were calculated by subtracting the youngest from the oldest Re–Os date. Thus, the large tonnage of these deposits could be the result of the superimposition of several ore-forming mid-to upper crustal magmatic systems. As the correlation of the deposit size and timescales of shallow magmatic–hydrothermal systems are currently ambiguous, we would argue that other variables could be the dominant factors controlling the deposit size, such as magma reservoir size, magma or fluid chemistry, fluid release, geological-tectonic situation and focusing mechanisms, or the metal precipitation efficiency. In that region (SW Serbia, North Macedonia), there are several magmatic occurrences with U–Pb ages between 23 and 35 Ma; the active Cu-porphyry mine at Buchim–Damjan–Borov Dol (North Macedonia) shows an Oligocene age of 24.5 Ma (Lehmann et al., 2013a), and the Cu-porphyry at Ilovitsa in SE (North Macedonia) has a magmatic Eocene/Oligocene age of 35 Ma.

Zircon crystallization over ∼ 200 kyr before the onset of ore formation recorded at the KMC is consistent with other high-precision geochronological studies on ore deposits (Figure 14: Buret et al., 2016; Tapster et al., 2016; Large et al., 2018, Large et al., 2021). The lack of variation observed in these deposits suggests the necessity of a long-lived and continuously crystallizing magma reservoir preceding the formation of economically significant ore deposits. The recorded ∼ 200 kyr of protracted zircon crystallization could indicate a period of volatile enrichment as a result of fractional crystallization and cooling of the magma reservoir before porphyry emplacement. The geochronological data from the KMC region are further evidence that magmatic emplacement (Phase I, Figure 10), reactivation of magma emplacement (Phase II), and ore formation (< 100 kyr) are the products of a longer-term evolution (a few 100 kyr, Figures 10, 14) of a large magma reservoir underlying the ore deposit that is the main driver of ore formation (von Quadt et al., 2011; Chelle-Michou et al., 2014; Buret et al., 2016; 2017; Tapster et al., 2016; Large et al., 2018). Magma reservoirs capable of forming porphyries and epithermal deposits are in turn part of a longer-term (several Myr) evolution of lithosphere-scale magma systems (Sasso, 1998; Rohrlach et al., 2005; Longo et al., 2010; Rezeau et al., 2016), which is consistent with the 2 Myr record of intrusive rocks preceding porphyry emplacement and ore formation recorded in e.g., the Batu Hijau district (Garwin, 2000; Wawryk and Foden, 2017).

FIGURE 14

Reconstructing the Chemical and Physical Evolution of the Magma Reservoir

Th/U and Yb/Gd and Ce/Sm ratios and Ti concentrations of zircon grains are commonly used as proxies for the degree of crystal fractionation within a magma reservoir (Samperton et al., 2015; Large et al., 2021). The decreasing Th/U ratios and increasing Yb/Gd ratios between samples and from cores to rims are indicative of progressive melt differentiation during zircon crystallization.

Trace element populations of zircons from the two rhyodacitic/granitic rocks demonstrate that the crystallizing magma at the time of emplacement of the older magmatics (EOKSC 1677, 1,687, and 17103) was hot and less evolved magma (Figure 12G) than at the time of emplacement of the younger quartz-monzonite dyke suite (i.e., 1.231 Myr after emplacement of the old porphyries, respectively). The good correlation of proxies indicating progressive differentiation (Th/U and Yb/Gd) with decreasing Ti concentrations (Figures 12E,F) indicates that the magma reservoir cooled during concurrent crystallization and melt evolution. In situ analyses of cores and rims are evidence of an evolving magma reservoir over the course of individual zircon crystallization (decreasing Ti; Figures 12E,F). Core–rim systematics of zircons from the old granite–rhyodacite further demonstrate cooling during protracted zircon growth (Figure 12). The coherent zircon trace element systematics recording melt differentiation over time are commonly inferred to result from zircon crystallization within a homogenous magma that evolved continuously (Figures 12, 14). The numerical model (Karakas et al., 2017) describes a system made up of a magma chamber in the upper crust that is connected with further deeper chambers in the lower crust. In principle, the model might indicate that all zircons analyzed in this study have crystallized from the same reservoir, where the zircons from the granite and rhyodacite reflect the earliest crystallized zircons from the least evolved melt. Hot source magma forms in the mantle or lower crust before rising through cracks and chimneys into the upper magma chamber (upper crust). It forms a reservoir, which cools down and partially crystallizes but can survive as a crystal mush for hundreds of thousands of years (Karakas et al., 2017).

However, the > 1.231 Myr time gap (Figure 10) is longer than the thermal lifetime of any recognized upper-crustal magmatic body (200–400 ky) (Wotzlaw et al., 2013; Large et al., 2020). We, therefore, consider two possible models: either the zircons within the rhyodacite/granite porphyry represent a separate crustal magmatic system not directly related to the younger ore-forming quartz-monzonitic intrusion system or the latter resulted from the same chamber that evolved through magma replenishment (model shown on Figure 15). Our results also require that the majority of zircon crystallization occurred in zircon-saturated reservoirs at deeper crustal levels prior to final magma emplacement and cooling, which has implications for using zircon U–Pb geochronology to infer upper crustal magma residence times (Barboni et al., 2015).

FIGURE 15

Irregular zircon trace element systematics in other intrusive magmatic settings have been explained with crystallization in non-homogenized and small melt batches sometimes, with contemporaneous incremental magma in addition to the mushy magma reservoir (Schoene et al., 2012; Tapster et al., 2016; Buret et al., 2016). Geochemically, similar zircon chemistries of granite, quartz-monzodiorite, or quartz–monzonites could also result from chemical stability as the magma reservoir reached the “petrological trap” at a crystallinity of ∼ 55–65% (Caricchi and Blundy, 2015), where the crystal fraction does not change over a broad temperature interval. Rim analyses that plot more outside the mineral co-crystallization trends than the respective core analyses (Figures 12B,C) could suggest latest crystallization within a nearly solidified magma that can be characterized by unsystematically variable trace element signatures (Buret et al., 2016; Lee et al., 2021). Alternatively, they could indicate thermal and possibly chemical rejuvenation of the magma. The latter would help for explaining the recorded thermal stability over tens of thousands of years. It is not possible to unambiguously identify one of the two mechanisms as dominant, and a concurrence of both is feasible. We, therefore, propose that in-between the emplacement of the old and young granite–quartz–monzonite, the underlying magma reservoir was in a thermally and chemically stable and crystal-rich state and was most likely affected by incremental magma underplating. This suggests that Au deposits form after a few hundred thousand years of cooling and crystallization, potentially within an originally melt-rich magma reservoir.

The ratios of Eu/Eu* and Ce/Sm in zircon are suggested in various studies (Ballard et al., 2002; Chelle-Michou et al., 2014) to record the oxidation state of the magma. Ce⁴⁺ is much more compatible in zircon than Ce³⁺, which is why zircon commonly presents a positive Ce anomaly (Hoskin and Schaltegger, 2003). The Cu–(Au) systems are often associated with oxidized magmas (Richards, 2011), which document the high Ce anomaly accompanied by decreasing Ti concentrations in case of the most evolved magmas. Additionally, due to the the higher uptake of MREEs than LREEs, in some minerals, the residual melt will enrich in Ce relative to Sm. However, an opposite trend at the KMC of decreasing Ce/Sm with decreasing age is recognized over all samples. The same trend is best visualized in the Th/U vs. Ce/Yb plot, where the rim values decrease from the volcanics to older sub-volcanics and finally to the quartz-monzonitic/granitic dykes. This strong depletion of Ce in combination of Th is most probably due to allanite fractionation (Figure 12). The observed weak Eu anomaly is more consistent due to commonly strong plagioclase fractionation, which nevertheless, is reduced to some extent by coeval titanite/apatite crystallization (Deering and Bachmann, 2010). The Hf isotopes support the presence of a cooling magma chamber with increased crustal assimilation with time because the ɛ-Hf values evolve from more “primitive” +3.5 ± 1.3 for the older volcanics and granite stocks/dykes to −0.8 ± 0.9 for the younger more evolved quartz-monzonitic-granitic rocks (Figure 13); the Late Oligocene porphyries at the Cu–Au Buchim deposit (North Macedonia) or Oligocene rocks in the Osogovo base metal region (Bulgaria) show similar evolution with crustal-dominated source material (Figure 13).

Model for the Formation of the Karavansalija Mineralized Center

Conclusion

The Pb–Zn–Ag–Au mineralization of the KMC at Rogozna Mt. in southwestern Serbia is linked to a shallow intrusive porphyritic intrusion, which is typical for the abundant metal ore deposits and minor precious metal deposits in the SMM-MB of the Oligocene age. The geochemical evolution of the ore-forming magma chamber could be tracked by analyzing zircon from initial zircon saturation through the emplacement of the various porphyries. In situ LA–ICP–MS and high-precision CA–ID–TIMS analyses lead to the discovery of the presence of two evolved magmas with distinct trace element patterns. A first stage of violent eruptive and shallow intrusive volcanic activity of granitic/rhyodacitic composition at ca. 29.34 Ma is followed by the intrusion of another porphyritic body at ca. 28.96 Ma trending toward shoshonitic composition into Cretaceous limestone, producing a skarn reaction zone. A separate magma chamber evolved with increased crustal assimilation recorded by the Hf isotopes and with a strong allanite fractionation trend seen on zircon trace element ratios. The input of a new melt into the base of the magma chamber increased the fluid pressure and likely provided additional fluids, metals, and sulfur before the emplacement of the quartz-monzonitic dykes into the skarn zone, precipitating the gold in a retrograde hydrothermal system starting at ca. 27.76 Ma. About 112 ka later (27.62 Ma), the last pulse of the quartz-monzonitic-porphyritic magma cuts the mineralization and records the maximum duration of mineralization at the KMC. These melts show signs of increased temperature and disequilibrium through the injection of a potassium-rich melt recorded from the abundant phlogopite micro-phenocrysts, small diopside clusters, and diopside reaction-rims on quartz recognized under the microscope. Nearby occurrences of significant base-metal mineralization, e.g., Crnac, Rudnik, or Veliki Majdan, are all thought to be associated with trachy-andesitic (quartz latitic) porphyry dykes but are not very precisely dated. The quartz-monzonitic dykes responsible for the mineralization at the KMC display many similarities with the described porphyries at Crnac (e.g., crowded phenocrysts, shoshonitic character, disequilibrium textures, and trace element pattern), and they likely share the same magmatic source responsible for their mineralization.

The geochronological data set shows that without CA–ID–TIMS dating the individual magmatic pulse responsible for the mineralization could not be identified; however, the LA–ICP–MS data help to see variation within a single zircon population and could reveal trends between core and rim measurements. Additionally, the laser technique records the homogeneity of the ablated material and therefore helps to interpret the ID–TIMS data where for example the input of xenocrystic or inherited cores cannot be determined. So far not many Cenozoic ore deposits in the Balkan area have been dated with the CA–TIMS U/Pb dating methods. This study has proven the benefits of a high-precision dating resolution in resolving different magmatic pulses within a single magmatic system. Consequently, the mineralization, which is typically associated with a magmatic intrusive event, can be traced and timed more precisely. Further work in the Balkan region involving U/Pb TIMS dating will help to delineate the various mineralizations with respect to their time of formation, link to magmatism, and the tectonic frame in which they were built.

References

• 1

  AndrićN.VogtK.MatencoL.CvetkovićV.CloetinghS.GeryaT. (2018). Variability of Orogenic Magmatism during Mediterranean-Style continental Collisions: A Numerical Modelling Approach. Gondwana Res.56, 119–134.

  • Google Scholar

• 2

  BachmannO.CharlierB. L. A.LowensternJ. B. (2007). Zircon Crystallization and Recycling in the Magma Chamber of the Rhyolitic Kos Plateau Tuff (Aegean Arc). Geol.35, 73–76. 10.1130/g23151a.1

  • CrossRef
  • Google Scholar

• 3

  BakerT. (2019). Gold ± Copper Endowment and Deposit Diversity in the Western Tethyan Magmatic Belt, Southeast Europe: Implications for Exploration. Econ. Geol.114, 1237–1250. 10.5382/econgeo.4643

  • CrossRef
  • Google Scholar

• 4

  BallardJ. R.PalinM. J.CampbellI. H. (2002). Relative Oxidation States of Magmas Inferred from Ce(IV)/Ce(III) in Zircon: Application to Porphyry Copper Deposits of Northern Chile. Contrib. Mineral. Petrol.144, 347–364. 10.1007/s00410-002-0402-5

  • CrossRef
  • Google Scholar

• 5

  BanikT. J.CobleM. A.MillerC. F. (2017). Porphyry Cu Formation in the Middle Jurassic Yerington Batholith, Nevada, USA: Constraints from Laser Raman, Trace Element, U-Pb Age, and Oxygen Isotope Analyses of Zircon. Geosphere13, 1113–1132. 10.1130/ges01351.1

  • CrossRef
  • Google Scholar

• 6

  BanješevićM.CvetkovićV.von QuadtA.ObradovićD. L.VasićN.PačevskiA.et al (2019). New Constraints on the Main Mineralization Event Inferred from the Latest Discoveries in the Bor Metallogenetic Zone (BMZ, East Serbia). Minerals9, 672. 10.3390/min9110672

  • CrossRef
  • Google Scholar

• 7

  BarboniM.SchoeneB.OvtcharovaM.BussyF.SchalteggerU.GerdesA. (2013). Timing of Incremental Pluton Construction and Magmatic Activity in a Back-Arc Setting Revealed by ID-TIMS U/Pb and Hf Isotopes on Complex Zircon Grains. Chem. Geology.342, 76–93. 10.1016/j.chemgeo.2012.12.011

  • CrossRef
  • Google Scholar

• 8

  BarboniM.AnnenC.SchoeneB. (2015). Evaluating the Construction and Evolution of Upper Crustal Magma Reservoirs with Coupled U/Pb Zircon Geochronology and thermal Modeling: A Case Study from the Mt. Capanne Pluton (Elba, Italy). Earth Planet. Sci. Lett.432, 436–448. 10.1016/j.epsl.2015.09.043

  • CrossRef
  • Google Scholar

• 9

  Borojević ŠoštarićS.CvetkovićV.NeubauerF.PalinkašL. A.BernroiderM.GenserJ. (2012). Oligocene Shoshonitic Rocks of the Rogozna Mts. (Central Balkan Peninsula): Evidence of Petrogenetic Links to the Formation of Pb–Zn–Ag Ore Deposits. Lithos148, 176–195.

  • Google Scholar

• 10

  BarraF.AlcotaH.RiveraS.ValenciaV.MunizagaF.MaksaevV. (2013). Timing and Formation of Porphyry Cu–Mo Mineralization in the Chuquicamata District, Northern Chile: New Constraints from the Toki Cluster. Miner. Deposita48, 629–651. 10.1007/s00126-012-0452-1

  • CrossRef
  • Google Scholar

• 11

  BowringJ. F.McLeanN. M.BowringS. A. (2011). Engineering Cyber Infrastructure for U-Pb Geochronology: Tripoli and U-Pb_Redux. Geochem. Geophys. Geosystems12. 10.1029/2010gc003479

  • CrossRef
  • Google Scholar

• 12

  BudinovZ. D.YonezuK.TindellT.Gabo-RatioJ. A.MilutinovicS.BoyceA. J.et al (2015). Copper–Gold Skarn Mineralization at the Karavansalija Ore Zone, Rogozna Mountain, Southwestern Serbia. Resour. Geol.65, 328–344.

  • Google Scholar

• 13

  BuretY.von QuadtA.HeinrichC.SelbyD.WälleM.PeytchevaI. (2016). From a long-lived upper-crustal magma chamber to rapid porphyry copper emplacement: Reading the geochemistry of zircon crystals at Bajo de la Alumbrera (NW Argentina). Earth Planet. Sci. Lett.450, 120–131. 10.1016/j.epsl.2016.06.017

  • CrossRef
  • Google Scholar

• 14

  BurtD. M. (1989). Compositional and phase relations among rare earth element minerals. Rev. Miner.21, 259–307.

  • Google Scholar

• 15

  BurkhardR. (2017). Style, Relative Timing and Precipitation Conditions of the Au Mineralisation at Karavansalija Ore Zone, Rogozna Mts. ETH Zurich, 54. SerbiaMaster.

  • Google Scholar

• 16

  CannellJ.CookeD. R.WalsheJ. L.SteinH. (2005). Geology, mineralization, alteration, and structural evolution of the El Teniente porphyry Cu-Mo deposit. Econ. Geol.100, 979–1003.

  • Google Scholar

• 17

  CaricchiL.BlundyJ. (2015). The Temporal Evolution of Chemical and Physical Properties of Magmatic Systems. Geol. Soc. Lond. Spec. Publ.422, 1–15. 10.1144/sp422.11

  • CrossRef
  • Google Scholar

• 18

  Chelle-MichouC.ChiaradiaM.OvtcharovaM.UlianovA.WotzlawJ.-F. (2014). Zircon Petrochronology Reveals the Temporal Link between Porphyry Systems and the Magmatic Evolution of Their Hidden Plutonic Roots (The Eocene Coroccohuayco deposit, Peru). Lithos198-199, 129–140. 10.1016/j.lithos.2014.03.017

  • CrossRef
  • Google Scholar

• 19

  ĆebićV. (1990). Konačni Izveštaj o geološko-petrološkom i Geohemijskom Izučavanju Tercijarnog Magmatskog Kompleksa Surduličke Oblasti. Fond Geoinstituta: Beograd.

  • Google Scholar

• 20

  Chelle-MichouC.RottierB.CaricchiL.SimpsonG. (2017). Tempo of Magma Degassing and the Genesis of Porphyry Copper Deposits. Sci. Rep.7, 40566. 10.1038/srep40566

  • CrossRef
  • Google Scholar

• 21

  ChiaradiaM.SchalteggerU.SpikingsR.WotzlawJ.-F.OvtcharovaM. (2013). How Accurately Can We Date the Duration of Magmatic-Hydrothermal Events in Porphyry Systems?--An Invited Paper. Econ. Geology.108, 565–584. 10.2113/econgeo.108.4.565

  • CrossRef
  • Google Scholar

• 22

  ChiaradiaM.SchalteggerU.SpikingsR. (2014). Time Scales of Mineral Systems—Advances in Understanding over the Past Decade.

  • Google Scholar

• 23

  ChiaradiaM. (2020). Stochastic Modelling of Deep Magmatic Controls on Porphyry Copper deposit Endowment: Scientific Reports, 11.

  • Google Scholar

• 24

  CondonD. J.SchoeneB.McLeanN. M.BowringS. A.ParrishR. R. (2015). Metrology and Traceability of U-Pb Isotope Dilution Geochronology (EARTHTIME Tracer Calibration Part I). Geochim. Cosmochim. Acta164, 464–480. 10.1016/j.gca.2015.05.026

  • CrossRef
  • Google Scholar

• 25

  CvetkovicA.IvanovicM.SaricK. (2006). The Role of Streptococcus Mutans Group and Salivary Immunoglobulins in Etiology of Early Childhood Caries. Serbian Dent J.53, 113–123. 10.2298/sgs0602113c

  • CrossRef
  • Google Scholar

• 26

  CvetkovićV.PoliG.ChristofidesG.KoroneosA.PecskayZ.Resimić-SarićK.et al (2007). The Miocene granitoid rocks of Mt. Bukulja (central Serbia): evidence for pannonian extension-related granitoid magmatism in the northern Dinarides. Eur. J. Mineral.19, 513–532.

  • Google Scholar

• 27

  CvetkovićV.PecskayZ.SaricK. (2016). Cenozoic igneous tectonomagmatic events in the Serbian part of the Balkan Peninsula: Inferences from K/Ar geochronology. Acta Vulcanologica25, 113–125.

  • Google Scholar

• 28

  DavidsonJ.TurnerS.HandleyH.MacphersonC.DossetoA. (2007). Amphibole “sponge” in arc crust?. Geology35, 787–790.

  • Google Scholar

• 29

  DeckartK.ClarkA. H.CuadraP.FanningM. (2012). Refinement of the time-space evolution of the giant Mio-Pliocene Rio Blanco-Los Bronces porphyry Cu–Mo cluster, Central Chile: new U–Pb (SHRIMP II) and Re–Os geochronology and 40Ar/39Ar thermochronology data. Miner. Deposita48, 57–79. 10.1007/s00126-012-0412-9

  • CrossRef
  • Google Scholar

• 30

  DeeringC. D.BachmannO. (2010). Trace Element Indicators of crystal Accumulation in Silicic Igneous Rocks. Earth Planet. Sci. Lett.297, 324–331. 10.1016/j.epsl.2010.06.034

  • CrossRef
  • Google Scholar

• 31

  DefantM. J.DrummondM. S. (1993). Mount St. Helens: Potential Example of the Partial Melting of the Subducted Lithosphere in a Volcanic Arc: Geology21, 547–550.

  • Google Scholar

• 32

  ErakD.MatencoL.ToljićM.StojadinovićU.AndriessenP. A. M.WillingshoferE.et al (2017). From Nappe Stacking to Extensional Detachments at the Contact between the Carpathians and Dinarides - the Jastrebac Mountains of Central Serbia. Tectonophysics710-711, 162–183. 10.1016/j.tecto.2016.12.022

  • CrossRef
  • Google Scholar

• 33

  FerryJ. M.WatsonE. B. (2007). New Thermodynamic Models and Revised Calibrations for the Ti-In-Zircon and Zr-In-Rutile Thermometers. Contrib. Mineral. Petrol.154, 429–437. 10.1007/s00410-007-0201-0

  • CrossRef
  • Google Scholar

• 34

  FreiR. (1992). Isotope (Pb, Rb-Sr, S, O, C, U-Pb) geochemical investigations on Tertiary 830 intrusives and related mineralizations in the Serbomacedonian Pb-Zn, Sb+ Cu-Mo metallogenetic province in northern Greece. PhD thesis, ETH Zurich.

  • Google Scholar

• 35

  GallhoferD.QuadtA. V.PeytchevaI.SchmidS. M.HeinrichC. A. (2015). Tectonic, Magmatic, and Metallogenic Evolution of the Late Cretaceous Arc in the Carpathian‐Balkan Orogen. Tectonics34, 1813–1836. 10.1002/2015tc003834

  • CrossRef
  • Google Scholar

• 36

  GarwinS. L. (2000). The Setting, Geometry and Timing of Intrusionrelated Hydrothermal System in the Vicinity of the Batu Hijau Porphyry Copper-Gold deposit. Sumbawa, Indonesia: University of Western Australia.

  • Google Scholar

• 37

  GehrelsG. E.ValenciaV. A.RuizJ. (2008). Enhanced Precision, Accuracy, Efficiency, and Spatial Resolution of U-Pb Ages by Laser Ablation-Multicollector-Inductively Coupled Plasma-Mass Spectrometry. Geochem. Geophys. Geosyst.9. 10.1029/2007gc001805

  • CrossRef
  • Google Scholar

• 38

  GelderI. E.MatencoL.WillingshoferE.TomljenovićB.AndriessenP. a. M.DuceaM. N.et al (2015). The Tectonic Evolution of a Critical Segment of the Dinarides‐Alps Connection: Kinematic and Geochronological Inferences from the Medvednica Mountains, NE Croatia. Tectonics34, 1952–1978. 10.1002/2015tc003937

  • CrossRef
  • Google Scholar

• 39

  GeorgievS.von QuadtA.HeinrichC. A.PeytchevaI.MarchevP. (2012). Time evolution of a rifted continental arc: Integrated ID-TIMS and LA-ICPMS study of magmatic zircons from the Eastern Srednogorie, Bulgaria: Lithos. Bulgaria: Lithos53, 57–67.

  • Google Scholar

• 40

  GrozdevV.PeytchevaI.von QuadtA.VassilevaR.GeorgievS. (2017). Analysis of the Zircon Population and Sr-Isotope Data of the Paleogene Igneous Rocks from Kyustendil Area, SW Bulgaria. Geol. Balc.46, 143–151.

  • Google Scholar

• 41

  HalterW. E.BainN.BeckerK.HeinrichC. A.LandtwingM.VonQuadtA.et al (2004). From Andesitic Volcanism to the Formation of a Porphyry Cu-Au Mineralizing Magma Chamber: the Farallón Negro Volcanic Complex, Northwestern Argentina. J. Volcanol. Geotherm. Res.136, 1–30. 10.1016/j.jvolgeores.2004.03.007

  • CrossRef
  • Google Scholar

• 42

  HörlerJ. (2017). The Karavansalija Ore Zone at Rogozna Mt. In Southwestern Serbia: Magma Evolution and Time Relationship of Intrusive Events, Skarn Mineralization and Overlying Volcanics:

  • Google Scholar

• 43

  HoskinP. W. O.SchalteggerU. (2003). 2. The Composition of Zircon and Igneous and Metamorphic Petrogenesis. Rev. Mineral. Geochem.53, 27–62. 10.1515/9781501509322-005

  • CrossRef
  • Google Scholar

• 44

  JankovićS. (1997). The Carpatho-Balkanides and Adjacent Area: a Sector of the Tethyan Eurasian Metallogenic belt. Miner. Deposita32, 426–433.

  • Google Scholar

• 45

  JelenkovićR.MilovanovićD.KoželjD.BanješevićM. (2016). The Mineral Resources of the Bor Metallogenic Zone: A Review. Geol. Croat.69, 143–155. 10.4154/gc.2016.11

  • CrossRef
  • Google Scholar

• 46

  KandicM.MicicI.MileticG.MarkovicD.MilutinovicS. (1987). Metallogenic Characteristics of the Rogozna Field. Geozavod-IMS, 1–126.

  • Google Scholar

• 47

  KarakasO.DegruyterW.BachmannO.DufekJ. (2017). Lifetime and Size of Shallow Magma Bodies Controlled by Crustal-Scale Magmatism. Nat. Geosci.10, 446–450. 10.1038/ngeo2959

  • CrossRef
  • Google Scholar

• 48

  KellerC. B.SchoeneB.SampertonK. M. (2018). A Stochastic Sampling Approach to Zircon Eruption Age Interpretation. Geochem. Persp. Let., 31–35. 10.7185/geochemlet.1826

  • CrossRef
  • Google Scholar

• 49

  KlimmK.BlundyJ. D.GreenT. H. (2008). Trace Element Partitioning and Accessory Phase Saturation during H2O-Saturated Melting of Basalt with Implications for Subduction Zone Chemical Fluxes. J. Petrol.49, 523–553. 10.1093/petrology/egn001

  • CrossRef
  • Google Scholar

• 50

  KoroneosA.PoliG.CvetkovićV.KrstićD.PECskayZ. (2010). Petrogenetic and tectonic inferences from the study of the Mt Cer pluton (West Serbia). Geol. Mag.148, 89–111.

  • Google Scholar

• 51

  LargeS. J. E.QuadtA. v.WotzlawJ.-F.GuillongM.HeinrichC. A. (2018). Magma Evolution Leading to Porphyry Au-Cu Mineralization at the Ok Tedi Deposit, Papua New Guinea: Trace Element Geochemistry and High-Precision Geochronology of Igneous Zircon. Econ. Geology.113, 39–61. 10.5382/econgeo.2018.4543

  • CrossRef
  • Google Scholar

• 52

  LargeS. J. E.WotzlawJ. F.GuillongM.QuadtA. v.HeinrichC. A. (2020). Resolving the Timescales of Magmatic and Hydrothermal Processes Associated with Porphyry deposit Formation Using Zircon U-Pb Petrochronology: Geochronology Discussions, 1–39.

  • Google Scholar

• 53

  LargeS. J. E.BuretY.WotzlawJ. F.KarakasO.GuillongM.von QuadtA.et al (2021). Copper-mineralised Porphyries Sample the Evolution of a Large-Volume Silicic Magma Reservoir from Rapid Assembly to Solidification. Earth Planet. Sci. Lett.563, 116877. 10.1016/j.epsl.2021.116877

  • CrossRef
  • Google Scholar

• 54

  Le MaitreR. W. B.DudekP.KellerA.LameyreJ.Le BasJ.SabineM. J.et al (1989). A Classification of Igneous Rocks and Glossary of Terms: Recommendations of the International Union of Geological Sciences, Subcommission on the Systematics of Igneous Rocks. International Union of Geological Sciences.

  • Google Scholar

• 55

  LeeR. G.DillesJ. H.TosdalR. M.WoodenJ. L.MazdabF. K. (2017). Magmatic Evolution of Granodiorite Intrusions at the El Salvador Porphyry Copper Deposit, Chile, Based on Trace Element Composition and U/Pb Age of Zircons. Econ. Geology.112, 245–273. 10.2113/econgeo.112.2.245

  • CrossRef
  • Google Scholar

• 56

  LeeR. G.ByrneK.D’AngeloM.HartC. J. R.HollingsP.GleesonS. A.et al (2021). Using Zircon Trace Element Composition to Assess Porphyry Copper Potential of the Guichon Creek Batholith and Highland Valley Copper deposit, South-central British Columbia. Miner Deposita56, 215–238. 10.1007/s00126-020-00961-1

  • CrossRef
  • Google Scholar

• 57

  LehmannS.BarcikowskiJ.von QuadtA.GallhoferD.PeytchevaI.HeinrichC. A.et al (2013). Geochronology, Geochemistry and Isotope Tracing of the Oligocene Magmatism of the Buchim-Damjan-Borov Dol Ore District: Implications for Timing, Duration and Source of the Magmatism. Lithos180-181, 216–233. 10.1016/j.lithos.2013.09.002

  • CrossRef
  • Google Scholar

• 58

  LesićV.MártonE.GajićV.JovanovićD.CvetkovV. (2019). Clockwise Vertical-axis Rotation in the West Vardar Zone of Serbia: Tectonic Implications. Swiss J. Geosci.112, 199–215. 10.1007/s00015-018-0321-8

  • CrossRef
  • Google Scholar

• 59

  LilovP.YanevY.MarchevP. (1987). K/Ar dating of the Eastern Rhodopes Paleogene magmatism. Geologica Balcanica17, 49–58.

  • Google Scholar

• 60

  LongoA. A.DillesJ. H.GrunderA. L.DuncanR. (2010). Evolution of calc-alkaline volcanism and associated hydrothermal gold deposits at Yanacocha, Peru. Econ. Geol.105, 1191–1241.

  • Google Scholar

• 61

  MaksaevV.MunizagaF.McWilliamsM. (2004). New Chronology for El Teniente, Chilean 916 Andes, from U-Pb, ⁴⁰Ar/³⁹Ar, Re-Os, and Fission-Track Dating: Implications for the Evolution of a Supergiant Porphyry Cu-Mo Deposit', in Sillitoe, R.H., Perello, J. and Vidal, C.E. (ed.). Andean Metallogeny: New Discoveries, Concepts, and Updates. Littleton, CO: Society of Economic Geologists, 15–54.

  • Google Scholar

• 62

  MarchevP.Raicheva R.DownesH.VaselliO.ChiaradiaM.MoritzR. (2015). Tectonic, Magmatic, and Metallogenic Evolution of the Late Cretaceous Arc in the Carpathian‐Balkan Orogen. Tectonics34, 1813–1836. 10.1002/2015tc003834

  • CrossRef
  • Google Scholar

• 63

  MaksaevV.MunizagaF.FanningM.PalaciosC.TapiaJ. (2006). SHRIMP U-Pb Dating of the Antucoya Porphyry Copper deposit: New Evidence for an Early Cretaceous Porphyry-Related Metallogenic Epoch in the Coastal Cordillera of Northern Chile. Miner Deposita41, 637–644. 10.1007/s00126-006-0091-5

  • CrossRef
  • Google Scholar

• 64

  MarchevP.Kaiser-RohrmeierM.HeinrichC.OvtcharovaM.von QuadtA.RaichevaR. (2005). 2: Hydrothermal Ore Deposits Related to post-orogenic Extensional Magmatism and Core Complex Formation: The Rhodope Massif of Bulgaria and Greece. Ore Geol. Rev.27, 53–89. 10.1016/j.oregeorev.2005.07.027

  • CrossRef
  • Google Scholar

• 65

  MarchevP.GeorgievS.RaichevaR.PeytchevaI.von QuadtA.OvtcharovaM.et al (2013). Adakitic magmatism in post-collisional setting: An example from the Early-Middle Eocene Magmatic Belt in Southern Bulgaria and Northern Greece. Lithos180, 159–180.

  • Google Scholar

• 66

  MarchevP.SingerB. (2002). 40Ar/39Ar Geochronology of Magmatism and Hydrothermal Activity of the Madjarovo Base-Precious Metal Ore District. Eastern Rhodopes, Bulgaria: Geological Society, London, Special Publications, 204, 137–150.

  • Google Scholar

• 67

  MarchevP.KibarovP.SpikingsR. A.OvtcharovaM.MartonI.MoritzR. (2010). 40Ar/39Ar and U-Pb Geochronology of the Iran Tepe Volcanic Complex. Eastern Rhodopes: Geologica Balcanica, 39, 3–12.

  • Google Scholar

• 68

  MatencoL.RadivojevićD. (2012). On the Formation and Evolution of the Pannonian Basin: Constraints Derived from the Structure of the junction Area between the Carpathians and Dinarides. Tectonics31. 10.1029/2012tc003206

  • CrossRef
  • Google Scholar

• 69

  MattinsonJ. M. (2005). Zircon U-Pb Chemical Abrasion ("CA-TIMS") Method: Combined Annealing and Multi-step Partial Dissolution Analysis for Improved Precision and Accuracy of Zircon Ages. Chem. Geol.220, 47–66. 10.1016/j.chemgeo.2005.03.011

  • CrossRef
  • Google Scholar

• 70

  McLeanN. M.BowringJ. F.BowringS. A. (2011). An Algorithm for U-Pb Isotope Dilution Data Reduction and Uncertainty Propagation. Geochem. Geophys. Geosyst.12. 10.1029/2010gc003478

  • CrossRef
  • Google Scholar

• 71

  McLeanN. M.CondonD. J.SchoeneB.BowringS. A. (2015). Evaluating Uncertainties in the Calibration of Isotopic Reference Materials and Multi-Element Isotopic Tracers (EARTHTIME Tracer Calibration Part II). Geochim. Cosmochim. Acta164, 481–501. 10.1016/j.gca.2015.02.040

  • CrossRef
  • Google Scholar

• 72

  MeinertL. D.DippleG. M.NicolescuS. (2005). World Skarn Deposits:10.5382/av100.11

  • CrossRef
  • Google Scholar

• 73

  MelfosV.VoudourisP. (2017). Cenozoic metallogeny of Greece and potential for precious, critical and rare metals exploration. Ore Geol. Rev.89, 1030–1057.

  • Google Scholar

• 74

  MladenovićA.TrivićB.CvetkovićV. (2015). How Tectonics Controlled post-collisional Magmatism within the Dinarides: Inferences Based on Study of Tectono-Magmatic Events in the Kopaonik Mts. (Southern Serbia). Tectonophysics646, 36–49. 10.1016/j.tecto.2015.02.001

  • CrossRef
  • Google Scholar

• 75

  MoritzR.NoverrazC.MártonI.MarchevP.SpikingsR.FontignieD.et al (2014). Sedimentary-rock-hosted Epithermal Systems of the Tertiary Eastern Rhodopes, Bulgaria: New Constraints from the Stremtsi Gold prospect. Geol. Soc. Lond. Spec. Publ.402, 207–230. 10.1144/sp402.7

  • CrossRef
  • Google Scholar

• 76

  PamicJ.BalenD. (2001). Tertiary magmatism of the Dinarides and the adjoining South Pannonian Basin: an overview. Acta Vulcan13, 9–24.

  • Google Scholar

• 77

  PatonC.WoodheadJ. D.HellstromJ. C.HergtJ. M.GreigA.MaasR. (2010). Improved Laser Ablation U-Pb Zircon Geochronology through Robust Downhole Fractionation Correction. Geochem. Geophys. Geosystems11. 10.1029/2009gc002618

  • CrossRef
  • Google Scholar

• 78

  PecskayZ.EleftheriadisG.KoroneosA.SoldatosT.ChristofidesG. (2003). K/Ar dating, geochemistry and evolution of the Tertiary volcanic rocks (Thrace, northeastern Greece): Mineral exploration and sustainable development. Millpress: Rotterdam. 1229–1232.

  • Google Scholar

• 79

  PrelevicD. (2005). Tertiary Ultrapotassic Volcanism in Serbia: Constraints on Petrogenesis and Mantle Source Characteristics. J. Petrol.46, 1443–1487.

  • Google Scholar

• 80

  PrelevicD.FoleyS. F.CvetkovicV.RomerR. L. (2004). Origin of Minette by Mixing of Lamproite and Dacite Magmas in Veliki Majdan, Serbia. J. Petrol.45, 759–792. 10.1093/petrology/egg109

  • CrossRef
  • Google Scholar

• 81

  ReidM. R.VazquezJ. A.SchmittA. K. (2011). Zircon-scale Insights into the History of a Supervolcano, Bishop Tuff, Long Valley, California, with Implications for the Ti-In-Zircon Geothermometer. Contrib. Mineral. Petrol.161, 293–311. 10.1007/s00410-010-0532-0

  • CrossRef
  • Google Scholar

• 82

  RezeauY.MoritzR.WotzlawJ.-F.TayanR.MelkonyanR.UlianovA.et al (2016). Temporal and genetic link between incremental pluton assembly and pulsed porphyry Cu-Mo formation in accretionary orogens. Geology44, 627–630.

  • Google Scholar

• 83

  RichardsJ. P.KerrichR. (2007). Special Paper: Adakite-Like Rocks: Their Diverse Origins and Questionable Role in Metallogenesis. Econom. Geol.102, 537–576.

  • Google Scholar

• 84

  RichardsJ. P. (2011). Magmatic to Hydrothermal Metal Fluxes in Convergent and Collided Margins. Ore Geol. Rev.40, 1–26. 10.1016/j.oregeorev.2011.05.006

  • CrossRef
  • Google Scholar

• 85

  RichardsJ. P. (2015). Tectonic, Magmatic, and Metallogenic Evolution of the Tethyan Orogen: From Subduction to Collision. Ore Geol. Rev.70, 323–345. 10.1016/j.oregeorev.2014.11.009

  • CrossRef
  • Google Scholar

• 86

  RickwoodP. C. (1989). Boundary Lines within Petrologic Diagrams Which Use Oxides of Major and Minor Elements. Lithos22, 247–263. 10.1016/0024-4937(89)90028-5

  • CrossRef
  • Google Scholar

• 87

  RohrlachB. D.LoucksR. R. (2005). “Multi-million-year Cyclic Ramp-Up of Volatiles in a Lower Crustal Magma Reservoir Trapped below the Tampakan Copper-Gold deposit by Mio-Pliocene Crustal Compression in the Southern Philippines: Super Porphyry Copper and Gold Deposits,” in A Global Perspective: Adelaide. PGC Publishing, 369–407.

  • Google Scholar

• 88

  RohrlachB. D.LoucksR. R. (2005). Multi-million-year cyclic ramp-up of volatiles in a lower crustal magma reservoir trapped below the Tampakan copper-gold deposit by Mio-Pliocene crustal compression in the southern Philippines: Super porphyry copper and gold deposits: A global perspective. Adelaide: PGC Publishing2, 369–407.

  • Google Scholar

• 89

  RohrmeierM. K.von QuadtA.DriesnerT.HeinrichC. A.HandlerR.OvtcharovaM.et al (2013). Post-orogenic Extension and Hydrothermal Ore Formation: High-Precision Geochronology of the Central Rhodopian Metamorphic Core Complex (Bulgaria-Greece). Econ. Geol.108, 691–718. 10.2113/econgeo.108.4.691

  • CrossRef
  • Google Scholar

• 90

  RottierB.KouzmanovK.OvtcharovaM.UlianovA.WälleM.SelbyD.et al (2020). Multiple rejuvenation episodes of a silicic magma reservoir at the origin of the large diatreme-dome complex and porphyry-type mineralization events at Cerro de Pasco (Peru). Lithos376-377, 105766. 10.1016/j.lithos.2020.105766

  • CrossRef
  • Google Scholar

• 91

  SassoA. (1998). The Farallón Negro Group, Northwest Argentina: Magmatic, Hydrothermal and Tectonic Evolution and Implications for Cu-Au Metallogeny in the Andean Back-arc. Soc. Econ. Geol. Newsletter34, 8–18.

  • Google Scholar

• 92

  SampertonK. M.SchoeneB.CottleJ. M.Brenhin KellerC.CrowleyJ. L.SchmitzM. D. (2015). Magma Emplacement, Differentiation and Cooling in the Middle Crust: Integrated Zircon Geochronological-Geochemical Constraints from the Bergell Intrusion, Central Alps. Chem. Geol.417, 322–340. 10.1016/j.chemgeo.2015.10.024

  • CrossRef
  • Google Scholar

• 93

  ScheferS.CvetkovićV.FügenschuhB.KounovA.OvtcharovaM.SchalteggerU.et al (2011). Cenozoic Granitoids in the Dinarides of Southern Serbia: Age of Intrusion, Isotope Geochemistry, Exhumation History and Significance for the Geodynamic Evolution of the Balkan Peninsula. Int. J. Earth Sci. (Geol Rundsch)100, 1181–1206. 10.1007/s00531-010-0599-x

  • CrossRef
  • Google Scholar

• 94

  SchmidS. M.FügenschuhB.KounovA.MaţencoL.NievergeltP.OberhänsliR.et al (2020). Tectonic Units of the Alpine Collision Zone between Eastern Alps and Western Turkey. Gondwana Res.78, 308–374. 10.1016/j.gr.2019.07.005

  • CrossRef
  • Google Scholar

• 95

  SchoeneB.SchalteggerU.BrackP.LatkoczyC.StrackeA.GüntherD. (2012). Rates of Magma Differentiation and Emplacement in a Ballooning Pluton Recorded by U-Pb TIMS-TEA, Adamello Batholith, Italy. Earth Planet. Sci. Lett.355-356, 162–173. 10.1016/j.epsl.2012.08.019

  • CrossRef
  • Google Scholar

• 96

  SchoeneB.CondonD. J.MorganL.McLeanN. (2013). Precision and Accuracy in Geochronology. Elements9, 19–24. 10.2113/gselements.9.1.19

  • CrossRef
  • Google Scholar

• 97

  SeghediI.DownesH.VaselliO.SzakácsA.BaloghK.PécskayZ. (2004). Post-collisional Tertiary-Quaternary Mafic Alkalic Magmatism in the Carpathian-Pannonian Region: a Review. Tectonophysics393, 43–62. 10.1016/j.tecto.2004.07.051

  • CrossRef
  • Google Scholar

• 98

  SerafimovskiT.StefanovaV.VolkovA. V. (2010). Dwarf copper-gold porphyry deposits of the Buchim-Damjan-Borov Dol ore district, Republic of Macedonia (FYROM). Geol. Ore Deposits52, 179–195.

  • Google Scholar

• 99

  SerafimovskiT.BoevB. (1996). Metallogeny of the Kratovo–Zletovo volcano-intrusive complex. Terranes Serbia356, 347–352.

  • Google Scholar

• 100

  SingerB.MarchevP. (2000). Temporal evolution of arc magmatism and hydrothermal activity, including epithermal gold veins, Borovitsa caldera, southern Bulgaria. Econom. Geol.95, 1155–1164.

  • Google Scholar

• 101

  SironC. R.RhysD.ThompsonJ. F. H.BakerT.VeligrakisT.CamachoA. (2018). Structural Controls on Porphyry Au-Cu and Au-Rich Polymetallic Carbonate-Hosted Replacement Deposits of the Kassandra Mining District, Northern Greece. Econom. Geol.113, 309–345.

  • Google Scholar

• 102

  StrmićPalinkašS.A PalinkašL.MandićM.Roller-LutzZ.PecskayZ. (2016). 1039 Origin and K-Ar age of the phreatomagmatic breccia at the Trepča Pb-Zn-Ag skarn deposit, Kosovo: Implications for ore-forming processes. Geologia Croatica69, 121–142.

  • Google Scholar

• 103

  StojadinovicU.MatencoL.AndriessenP. A. M.ToljićM.FoekenJ. P. T. (2013). The Balance between Orogenic Building and Subsequent Extension during the Tertiary Evolution of the NE Dinarides: Constraints from Low-Temperature Thermochronology. Glob. Planet. Change103, 19–38. 10.1016/j.gloplacha.2012.08.004

  • CrossRef
  • Google Scholar

• 104

  StojadinovicU.MatencoL.AndriessenP.ToljićM.RundićL.DuceaM. N. (2017). Structure and Provenance of Late Cretaceous-Miocene Sediments Located Near the NE Dinarides Margin: Inferences from Kinematics of Orogenic Building and Subsequent Extensional Collapse. Tectonophysics710-711, 184–204. 10.1016/j.tecto.2016.12.021

  • CrossRef
  • Google Scholar

• 105

  SunS.-s.McDonoughW. F. (1989). Chemical and Isotopic Systematics of Oceanic Basalts: Implications for Mantle Composition and Processes. Geol. Soc. Lond. Spec. Publ.42, 313–345. 10.1144/gsl.sp.1989.042.01.19

  • CrossRef
  • Google Scholar

• 106

  SzymanowskiD.WotzlawJ.-F.EllisB. S.BachmannO.GuillongM.von QuadtA. (2017). Protracted Near-Solidus Storage and Pre-eruptive Rejuvenation of Large Magma Reservoirs. Nat. Geosci.10, 777–782. 10.1038/ngeo3020

  • CrossRef
  • Google Scholar

• 107

  TapsterS.CondonD. J.NadenJ.NobleS. R.PettersonM. G.RobertsN. M. W.et al (2016). Rapid thermal Rejuvenation of High-Crystallinity Magma Linked to Porphyry Copper deposit Formation; Evidence from the Koloula Porphyry Prospect, Solomon Islands. Earth Planet. Sci. Lett.442, 206–217. 10.1016/j.epsl.2016.02.046

  • CrossRef
  • Google Scholar

• 108

  TaghipourB.MooreF.MackizadehM. A.TaghipourS. (2013). Hydrothermal garnet in porphyry copper related skarn deposits, Ali-Abad, Yazd Province, Iran. Iran. J. Sci. Technol.A1, 11–22.

  • Google Scholar

• 109

  ToljićM.MatencoL.DuceaM. N.StojadinovićU.MilivojevićJ.ĐerićN. (2013). The Evolution of a Key Segment in the Europe–Adria Collision: The Fruška Gora of Northern Serbia. Glob. Planet. Change103, 39–62. 10.1016/j.gloplacha.2012.10.009

  • CrossRef
  • Google Scholar

• 110

  ToljićM.MatencoL.StojadinovićU.WillingshoferE.Ljubović-ObradovićD. (2018). Understanding Fossil Fore-Arc Basins: Inferences from the Cretaceous Adria-Europe Convergence in the NE Dinarides. Glob. Planet. Change171, 167–184.

  • Google Scholar

• 111

  UrosevićM.PavlovićZ.KlisićM.KaramataS.MalesevicM.StefanovićM.et al (1973). Basic Geological Map 1:100,000, Sheet Novi Pazar. Explanatory Note. Belgrade: Federal Geological Institute, 1–77.

  • Google Scholar

• 112

  UstaszewskiK.KounovA.SchmidS. M.SchalteggerU.KrennE.FrankW.et al (2010). Evolution of the Adria-Europe Plate Boundary in the Northern Dinarides: From Continent-Continent Collision to Back-Arc Extension. Tectonics29. 10.1029/2010tc002668

  • CrossRef
  • Google Scholar

• 113

  von QuadtA.MoritzR.PeytchevaI.HeinrichC. A. (2005). 3: Geochronology and Geodynamics of Late Cretaceous Magmatism and Cu-Au Mineralization in the Panagyurishte Region of the Apuseni-Banat-Timok-Srednogorie belt, Bulgaria. Ore Geology. Rev.27, 95–126. 10.1016/j.oregeorev.2005.07.024

  • CrossRef
  • Google Scholar

• 114

  von QuadtA.ErniM.MartinekK.MollM.PeytchevaI.HeinrichC. A. (2011). Zircon Crystallization and the Lifetimes of Ore-Forming Magmatic-Hydrothermal Systems. Geology39, 731–734. 10.1130/g31966.1

  • CrossRef
  • Google Scholar

• 115

  von QuadtA.GallhoferD.GuillongM.PeytchevaI.WaelleM.SakataS. (2014). U-pb Dating of CA/non-CA Treated Zircons Obtained by LA-ICP-MS and CA-TIMS Techniques: Impact for Their Geological Interpretation. J. Anal. Spectrom.29, 1618–1629. 10.1039/c4ja00102h

  • CrossRef
  • Google Scholar

• 116

  VoudourisP.MavrogonatosC.SpryP. G.BakerT.MelfosV.KlemdR.et al (2019). Porphyry and epithermal deposits in Greece: An overview, new discoveries, and mineralogical constraints on their genesis: Ore Geology Reviews, v. 107, p. 654–691. Ore Geol. Rev.107, 654–691.

  • Google Scholar

• 117

  WinchesterJ. A.FloydP. A. (1976). Geochemical Magma Type Discrimination: Application to Altered and Metamorphosed Basic Igneous Rocks. Earth Planet. Sci. Lett.28, 459–469. 10.1016/0012-821x(76)90207-7

  • CrossRef
  • Google Scholar

• 118

  WeisP.DriesnerT.HeinrichC. A. (2012). Porphyry-Copper Ore Shells Form at Stable Pressure-Temperature Fronts Within Dynamic Fluid Plumes. Science338, 1613–1616.

  • Google Scholar

• 119

  WotzlawJ.-F.SchalteggerU.FrickD. A.DunganM. A.GerdesA.GüntherD. (2013). Tracking the Evolution of Large-Volume Silicic Magma Reservoirs from Assembly to Supereruption. Geology41, 867–870. 10.1130/g34366.1

  • CrossRef
  • Google Scholar