Next Article in Journal
The Viscosity and Atomic Structure of Volatile-Bearing Melilititic Melts at High Pressure and Temperature and the Transport of Deep Carbon
Previous Article in Journal
Provenance and Implication of Carboniferous–Permian Detrital Zircons from the Upper Paleozoic, Southern Ordos Basin, China: Evidence from U-Pb Geochronology and Hf Isotopes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Stage Late Jurassic to Early Cretaceous Hydrothermal Activity in the Sakar Unit of Southeastern Bulgaria

1
Institute of Earth Sciences, Faculty of Natural Sciences, University of Silesia in Katowice, Będzińska 60, 41-200 Sosnowiec, Poland
2
Institute of Geological Sciences PAS, 00-818 Warsaw, Poland
3
Department of Geology, School of Natural Sciences, Trinity College Dublin, 2 Dublin, Ireland
4
The Earth and Man National Museum in Sofia, 1421 Sofia, Bulgaria
5
Hristo Smirnenski 37-41, G, 51 Topolovgrad, Bulgaria
*
Author to whom correspondence should be addressed.
Minerals 2020, 10(3), 266; https://doi.org/10.3390/min10030266
Submission received: 24 December 2019 / Revised: 23 February 2020 / Accepted: 4 March 2020 / Published: 15 March 2020
(This article belongs to the Section Mineral Geochemistry and Geochronology)

Abstract

:
Southeastern Bulgaria is composed of a variety of rocks from pre-Variscan (ca. 0.3 Ga) to pre-Alpine sensu lato (ca. 0.15 Ga) time. The Sakar Unit in this region comprises a series of granitoids and gneisses formed or metamorphosed during these events. It is cut by a series of post-Variscan hydrothermal veins, yet lacks pervasive Alpine deformation. It thus represents a key unit for detecting potential tectonism associated with the enigmatic Cimmerian Orogenic episode, but limited geochronology has been undertaken on this unit. Here we report age constraints on hydrothermal activity in the Sakar Pluton. The investigated veins contain mainly albite–actinolite–chlorite–apatite–titanite–quartz–tourmaline–epidote and accessory minerals. The most common accessory minerals are rutile and molybdenite. Apatite and titanite from the same vein were dated by U–Pb LA–ICP-MS geochronology. These dates are interpreted as crystallization ages and are 149 ± 7 Ma on apatite and 114 ± 1 Ma on titanite, respectively. These crystallization ages are the first to document two stages of hydrothermal activity during the late Jurassic to early Cretaceous, using U–Pb geochronology, and its association with the Cimmerian orogenesis. The Cimmerian tectono-thermal episode is well-documented further to the east in the Eastern Strandja Massif granitoids. However, these are the first documented ages from the western parts of the Strandja Massif, in the Sakar Unit. These ages also temporally overlap with previously published Ar–Ar and K–Ar cooling ages, and firmly establish that the Cimmerian orogeny in the studied area included both tectonic and hydrothermal activity. Such hydrothermal activity likely accounted for the intense albitization found in the Sakar Unit.

1. Introduction

The Balkan Peninsula consists of a series of variably metamorphosed terranes that predominantly formed in the Phanerozoic due to the collision of Gondwana against Baltica, then Laurussia, and later proto-Africa against proto-Europe during the Caledonian, Variscan, and Alpine orogeny [1]. Southern Bulgaria is divided into the Rhodope and Strandja massifs, with the relationship between them remaining uncertain [2]. These massifs are composed of highly deformed orthogneisses and paragneisses. The highly metamorphosed and deformed pre-Variscan basement is intruded by Variscan granites, which formed as a result of the collision of Gondwana and Laurussia during Pangea assembly [3]. Later, Cimmerian-aged events preceded the main collision of Africa and Europe during the Alpine Orogeny [4]. The Cimmerian Orogeny (~200–150 Ma; [5] and references therein) occurred to the south and east of Bulgaria, with most Cimmerian tectonism documented between Turkey and Central Asia [4,5,6,7]. The enigmatic Cimmerian Orogeny in Southeastern Europe, however, is difficult to identify mostly due to the pervasive later Alpine orogenic overprint [8]. In Southeastern Bulgaria, the late Variscan Sakar pluton intrudes the pre-Variscan basement [9] and is in turn affected by later events of deformation and metamorphism, which in the past have been attributed to the Alpine Orogeny [10]. Many of the Sakar granitoids are affected by a series of veins and pervasive albitization, connected with Na-metasomatism of undetermined age [11,12]. Herein, we provide the first U–Pb LA-ICP-MS dates from apatite and titanite in one vein type, documenting two episodes of hydrothermal activity within the Sakar Unit of the Strandja Massif, which critically escaped significant Alpine overprinting.

2. Geological Setting

The Sakar Unit of the Strandja Zone [13], also known as the Strandja (Strandzha) Massif (SM; [14,15] and references therein) or Istranca Massif (IM; [16]), is located in the southeastern part of Bulgaria, close to the border with Turkey and Greece (Figure 1A). It belongs to the East Srednogorie Zone of the Balkan Orogenic System [17,18]. The Sakar Unit is exposed in the Sakar Mountains, which neighbor the Rhodope and Strandja Mountains, and is bounded by the Maritsa and Tundzha rivers (Figure 1A). The Sakar Unit consists of the Sakar Pluton, which intruded into high-grade meta-igneous and metasedimentary rocks, and in turn are unconformably overlain by Triassic to middle Jurassic metasedimentary rocks (Figure 1B). The Sakar Pluton is an East–West elongated magmatic body, ca. 20 km long and up to ca. 15 km wide, forming the major component of the Sakar Unit. In the west and north of the unit, the contact between the pluton and the host rocks is covered by Neogene sedimentary rocks (Figure 1), whereas in the northeast and south, the intensely sheared contacts of the pluton are exposed [13,19,20]. The Sakar Pluton consists of various granitoids divided into three groups: equigranular granitoids, porphyritic granitoids also with feldspar megacrysts, and aplitic leucogranites. The Sakar granitoids contain xenoliths of the country rock and are cut by numerous aplite veins and some mafic dykes [11]. Locally, the granitoids are affected by extensive albitization (e.g., documented from the northeastern part of pluton, south of the Orlov dol village and Kanarata quarry near Hlyabovo [11,12]). The timing of albitization is unknown; however, this process also affected the low-grade Triassic metasedimentary rocks, and is therefore post-Triassic [11,12].
The crystallization age of the Sakar Pluton, as well as the age of the meta-igneous and metasedimentary rocks has been the subject of long debate. Isotope dating by the Rb–Sr method from the Sakar Pluton has yielded contradictory results of 320 ± 18 Ma [21], 499 ± 70 Ma [21], and 250 ± 35 Ma [22], which is discussed by Ivanov et al. [23], who proposed a post-middle-Triassic emplacement age for the Sakar Pluton. K–Ar ages range between ca. 144 Ma and ca. 111 Ma [21,24,25,26,27] and likely reflect a regional thermal overprint. However, preliminary U–Pb zircon dating results presented by Peytcheva et al. [28], Bonev et al. [9], and Pristavova et al. [20] yield late-Carboniferous to early Permian (ca. 305 Ma to ca. 295 Ma) crystallization ages for the Sakar Pluton. These new results contradict the timing of supposed Jurassic to early Cretaceous emplacement of the granitoid presented by Ivanov [23] which was supported by field observations of Gerdjikov [13]. These new late-Variscan (ca. 300 Ma) ages also provide a minimum age for the host rocks of the Sakar Pluton (the Strandja Massif). The Strandja Massif has widely been accepted as Precambrian in age [29], although there are no modern geochronological studies, apart from the recent study by Bonev et al., [9], which indicates the presence of an Ordovician magmatic event in the massif. Thus the Paleozoic, pre-Carboniferous evolution of the Strandja Massif thus remains uncertain.

3. Sample Locality and Methods

Rock samples for this study were collected from the Kanarata quarry (Figure 1C), located ca. 3 km north of Hlyabovo village and ca. 5 km east of Topolovgrad (Haskovo District), within the northeastern part of the Sakar Pluton. The quarry is located near the contact with the country rock (<500 m to the west of the quarry), and contains granitoids with abundant amphibolite xenoliths, which are a common feature of the margin of the pluton. The quarry is currently active, with the granites being exploited as dimension stone. The locality is also known for the occurrence of gemstone-quality apatite ([30]; Figure 2A).
Pale-blue apatite and olive-green titanite crystals (Figure 2B) from the Hlyabovo (Alpine-type) fissures were collected for U–Pb laser ablation–inductively coupled plasma mass spectrometry (LA-ICP-MS) dating.

3.1. Microscopy

Petrographic analyses of thin sections, using light microscopy, were undertaken at the Faculty of Natural Sciences in the University of Silesia in Katowice, Sosnowiec, Poland, using an Olympus BX−51 microscope to constrain mineralogical, textural, and microstructural relationships of the samples and to determine the presence of dateable accessory minerals. The petrographic observations were used to select representative minerals for subsequent electron probe microanalysis (EMPA) and whole-rock major and trace-element geochemical analyses.

3.2. Electron Probe Micro-Analyses (EMPA)

Microprobe analyses of the main rock-forming and accessory minerals were carried out at the Inter-Institutional Laboratory of Microanalyses of Minerals and Synthetic Substances, University of Warsaw, Warsaw, using a CAMECA SX−100 electron microprobe. The analytical conditions employed an accelerating voltage of 15 kV, a beam current of 20 nA, counting times of 4 s for peak and background, and a beam diameter of 1–5 µm. The full analytical procedures are listed in Table S1.

3.3. LA–ICP-MS U–Pb Apatite and Titanite Dating

U–Pb apatite and titanite ages were determined by using a Photon Machines Analyte Exite 193 nm ArF Excimer laser-ablation system with a Helex 2-volume ablation cell coupled to an Agilent 7900 ICP-MS, using an Aerosol Rapid Introduction System (ARIS; Teledyne Photon Machines), at the Department of Geology, Trinity College, Dublin, Ireland. The ICP-MS was tuned, using NIST 612 standard glass to yield Th/U ratios of unity and low oxide production rates (ThO+/Th+ typically <0.15%). A 0.15 L/min He carrier gas was fed into the cell body, and 0.25 L/min He was fed into the cup, the aerosol was subsequently mixed with 0.6 L/min Ar make-up gas and a small volume of N2 (~7 mL/min) in a ~1 cm3 signal-smoothing device coupled to the ARIS, to enhance signal sensitivity and reduce oxide formation. For apatite, the following 30 isotopes were measured (with their respective dwell times in milliseconds listed in parentheses): 31P (10), 35Cl (20), 43Ca (10), 51V(7.5), 55Mn (75), 88Sr (7.5), 89Y (5), 90Zr (10), 137Ba (10), 139La (5), 140Ce (5), 141Pr (5), 146Nd (5), 147Sm (7.5), 153Eu (10), 157Gd (10), 159Tb (10), 163Dy (10), 165Ho (10), 166Er (10), 169Tm (10), 172Yb (10), 175Lu (15), 202Hg (5), 204Pb (5), 206Pb (75), 207Pb (75), 208Pb (5), 232Th (5), and 238U (75), corresponding to a total dwell time of 456 ms and a total cycle time of 510 ms. For titanite, the following 25 isotopes were measured (with their respective dwell times in milliseconds listed in parentheses): 43Ca (10), 49Ti(10), 90Zr (5), 139La (5), 140Ce (5), 141Pr (5), 146Nd (5), 147Sm (5), 153Eu (5), 157Gd (5), 159Tb (5), 163Dy (5), 165Ho (7), 166Er (7), 169Tm (7), 172Yb (7), 175Lu (10), 202Hg (25), 204Pb (25), 206Pb (50), 207Pb (70), 208Pb (20), 232Th (10), and 238U (30), corresponding to a total dwell time of 343 ms and a total cycle time of 379.5 ms. For both titanite and apatite analyses, the laser fluence was 2.5 J/cm2, with a repetition rate of 15 Hz, a 47 μm spot size, and an analysis time of 20 s, followed by an 8 s background measurement.
The raw isotope data were reduced, using the “VizualAge UcomPbine” data-reduction scheme (DRS) of Chew et al. [31], a modification of the U–Pb geochronology “VizualAge” DRS of Petrus and Kamber [32] that can account for the presence of variable common Pb in the primary standard materials. The DRS runs within the freeware IOLITE package of Paton et al. [33]. In IOLITE, user-defined time intervals are established for the baseline correction procedure to calculate session-wide baseline-corrected values for each isotope. The time-resolved fractionation response of individual standard analyses is then characterized, using a user-specified down-hole correction model (such as an exponential curve, a linear fit, or a smoothed cubic spline). The VizualAge data-reduction scheme then fits this appropriate session-wide “model” U–Th–Pb fractionation curve to the time-resolved standard data and the unknowns. Sample-standard bracketing is applied after the correction of down-hole fractionation to account for long-term drift in isotopic or elemental ratios by normalizing all ratios to those of the U–Pb reference materials. Common Pb in the apatite and titanite standards was corrected by using the 207Pb-based correction method.
Blocks of nine standards and two NIST612 standard glass analyses were followed by 20 unknown samples. For apatite analyses, the Madagascar natural mineral standard [34] was used (weighted mean ID-TIMS concordia age of this crystal is 473.5 ± 0.7 Ma). McClure Mountain syenite apatite (weighted mean 207Pb/235U age of 523.51 ± 2.09 Ma [35], and Durango apatite (31.44 ± 0.18 Ma [36], were used as secondary standards. As the in-house standards, McClure Mountain apatite yielded a U–Pb Tera–Wasserburg concordia lower intercept age of 523.5 ± 3.0 Ma (MSWD = 2.4), anchored using a 207Pb/206Pb value of 0.88198 derived from an apatite ID-TIMS total U–Pb isochron [35]. The Durango apatite standard yielded a U–Pb Tera–Wasserburg concordia lower intercept age of 33.43 ± 0.62 Ma (MSWD = 1.8), anchored using a 207Pb/206Pb value of value of 0.838 derived from the Stacey and Kramers [37] terrestrial Pb evolution model at an age of 31.44 Ma. For titanite analyses, MKED1 titanite (207Pb/206Pb TIMS age of 1521.02 ± 0.55 Ma; [38]) was used as the primary standard. OLT−1 titanite (U–Pb TIMS concordia age of 1014.8 ± 2.0 Ma; [39]) and BLR titanite (U–Pb TIMS age of 1047.1 ± 0.4 Ma; [40] were used as the secondary U–Pb calibration standards and yielded weighted average 207Pb-corrected ages of 1017.4 ± 6.7 Ma (n = 17, MSWD = 3.4) and 1067.9 ± 6.2 Ma (n = 14, MSWD = 2.4), respectively.

4. Results

4.1. Petrography

The studied rocks are hosted by the so-called Aleksandrovo-type granites (see [12] and references therein), which in Kanarata quarry are represented by medium-grained albitized granitoids (Figure 3A). These rocks contain amphibolite xenoliths from tens of centimeters up to 3 m in length (Figure 3B). All of the xenoliths have sharp contacts with the country rock and exhibit a distinctive mineral lineation (Figure 3B,C).
All investigated veins are Alpine-type fissures, varying from a few centimeters to up to 1 m in length (Figure 3A,C). These veins are rich in albite, actinolite, chlorite, apatite, and titanite. The albite is composed of milky to translucent crystals, with an anorthite content up to 3.5 wt.% (Table S2). The albite crystals are euhedral, and up to 1 cm in length (Figure 2D,E). Actinolite (Mg-rich hornblende; Table S3), forms different varieties of pale- to dark-green crystals. They typically have a needle-like habit, with long prismatic crystals up to 2–3 cm in length and an aspect ratio as high as 100:1. The actinolite commonly occurs as inclusions in albite. The chlorite in the vein contains 50%–82% of the clinochlore (Mg5Al)(Si3Al)O10(OH) end-member, 12%–49% of the chamosite (Fe2+5Al)(Si3Al)O10(OH) end-member, and <3% of the pennantite (Mn5Al)(Si3Al)O10(OH) end-member (classification after Bailey [41]; Table S4). The chlorite flakes partially overgrow albite and actinolite, and mostly form aggregates up to 300 µm in length. Chlorite thermometry, using the Cathelineau and Nieva’s [42], Jowett’s [43], and Karaniditis et al.’s [44] methodologies, yielded temperature ranges of 298–340 °C, 286–329 °C, and 272–300 °C, respectively (Table S4; Figure 5). Apatite generally forms well-shaped hexagonal milky to bluish crystals, which commonly are prismatic, dipyramidal, and stubby. Sporadically, flat or tabular plates were also noted (Figure 2A). The apatite crystals are up to 5–7 cm in length, with an aspect ratio of 4:1. Twinned crystals are lacking. All investigated crystals are fluorapatite (Ca-apatite; Table 1), and imaging by EMPA reveals the chemical homogeneity of the investigated apatite crystals. Titanite from the same locality usually forms euhedral crystals, twinned on {221}. Some of the crystals may be in flattened into tabular or even prismatic forms, with pointed terminations and complex dipyramidal crystals. The titanite crystals are pale to dark green, and are up to 2 cm in length (Figure 2B). The crystals do not reveal any chemical zonation (Table S2). The titanite in most cases is not overgrown by any associated minerals. Titanite crystals contain more chlorite inclusions than the apatite, regardless of the crystal size. Moreover, both phases do not form intergrowths between each other. The presence of rutile is limited to prismatic metallic crystals, with lengths up to 200 µm. The rutile crystals are often strongly weathered. This type of vein (the Alpine-type fissures) is often located close to amphibolite xenoliths.
Sporadically, veins composed mainly of quartz, epidote–clinozoisite, tourmaline (schörl), albite, rutile, and molybdenite are also noted, although these were not dated in this study. The groundmass of these veins is mainly quartz. Barren quartz sporadically forms euhedral crystals (up to 4–5 cm in length) that are smoky to in color transparent. The more common massive quartz is mostly gray–milky in color (Figure 2F,G). Quartz is associated only with this type of vein. Albite forms white crystals, which are orientated perpendicular to the contact with the host rock. A similar situation is observed if rutile crystals are present, with the rutile crystals reaching up to 4 cm in length (Figure 2E). Prismatic, elongated crystals of rutile, with a prominent trigonal prism and pyramid, were observed. The surface of the rutile crystals appears to be oxidized and coated by a thin film of Fe-hydroxides (goethite/limonite?). Epidote group minerals are present in this vein-type as radial aggregates of acicular to columnar crystals, having pale-green to green color (Figure 2E). Molybdenite in this vein type is up to 2.5 cm across, and the crystals are commonly tabular- and barrel-shaped (Figure 2F). Tourmaline group minerals are represented by schörl, which forms prismatic to acicular crystals, and sporadically may be modified along the [0001] plane, with a prominent trigonal prism and pyramid. The crystals are up to 8 cm in length, with an aspect ratio of 10:1. Tourmaline mostly occurs only with quartz.

4.2. U–Pb Geochronology

The homogeneity and size of the titanite (Table 1) and apatite (Table 2) grains in this study make them useful for U–Pb dating. Sixty titanite grains were analyzed, and all the data points are highly discordant due to high proportions of common Pb relative to radiogenic Pb. All data lie on an unanchored Tera–Wasserburg discordia line with a lower intercept age of 114 ± 1 Ma (MSWD = 1.3) and a 207Pb/206Pb initial ratio of 0.6799 ± 0.011 (Figure 4A).
Apatite grains were also analyzed, and the 49 data points were also highly discordant due to the high proportion of common Pb compared to radiogenic Pb. The data points all fall on an unanchored Tera–Wasserburg discordia, with a lower intercept age of 149 ± 7 Ma (MSWD = 1.7) and a 207Pb/206Pb initial ratio of 0.772 ± 0.021 (Figure 4B). All LA–ICP-MS data are listed in Table S5.

5. Discussion

5.1. Age Significance: Apatite vs. Titanite

Titanite (CaTiSiO5) and apatite (Ca5(PO4)3(Cl/F/OH)) are common accessory minerals in a compositionally wide range of igneous, metamorphic (especially in amphibolite facies), and hydrothermal rocks. Both minerals are also important carriers for numerous trace elements, including U, Th, and Pb. Titanite commonly contains trace amounts of U and largely excludes initial Pb, thus making it useful for radiometric dating using the U–Pb–Th isotopic system (e.g., [45,46,47]). Multi-equilibrium thermobarometry suggests that the titanite closure temperature falls within a restricted range of 660–700 °C (e.g., [48]). For apatite, experimental determination of the diffusion parameters of Pb by high-temperature annealing experiments [49] and ion implantation with Rutherford backscattering techniques [50] imply closure temperatures of between 450 and 550 °C. Apatite is chemically stable at middle-amphibolite facies conditions (i.e., temperatures above its closure temperature). Under such conditions, it is thought the U–Pb systematics of apatite are controlled predominantly by volume diffusion rather than by the new growth or recrystallization [51]. Taking the above factors into consideration, the titanite and apatite ages represent separate crystallization episodes in the same vein type, as the calculated ages lie significantly outside of their mutual analytical uncertainties. In addition, the significantly lower 207Pb/206Pb initial ratio of the titanite (suggesting derivation from a radiogenic fluid), also implies separate crystallization episodes.
The Alpine-type fissure veins contain abundant apatite, titanite, albite, actinolite, and chlorite, and they exhibit similarities with vein systems across most orogenic zones connected with the Alpine collision in most of Southern Europe, from Italy to Greece, to Turkey, and onward to Iran and China (e.g., [52,53,54,55,56]). This veining is associated with hydrothermal activity, which is interpreted in this study as being the cause of Na-albitization in the country rocks. The occurrence of two mineral phases of different ages in one druse indicates two discrete hydrothermal crystallization episodes. Chlorite temperature calculations show that the most probable temperature for hydrothermal activity occurred in the range of ca. 290–330 °C (Figure 5).

5.2. Geotectonic Implications of the Apatite and Titanite Ages

The margins of the late-Carboniferous to early Permian Sakar Pluton (305–295 Ma [9,12,28]; Figure 6A) are affected by extensive Na-albitization, as documented by Kamenov et al., [11] and Pristvova et al., [12], along with the vein systems such as those are observed at Kanarata quarry (Figure 6B,C). The obtained U–Pb LA-ICP-MS ages of 149 ± 7 Ma on apatite and 114 ± 1 Ma on titanite, respectively, suggest two discrete crystallization ages within the Sakar Pluton. The hydrothermal veins, therefore, represent either a long-lived (i.e., ca. 35 Myr) hydrothermal episode within the Sakar Unit in the late-Jurassic to early Cretaceous, or two different episodes of fracturing associated with hydrothermal activity. Albitization from Kanarata quarry was documented by Pristavova et al., [12]; however, the timing of this alteration process remained unknown. As field relationships suggest the spatial and temporal link of albitization and veining here, we suggest the hydrothermal activity within the Sakar Unit is coeval with K–Ar ages between ca. 144 and ca. 111 Ma [21,22,24,25,26,27]. Similar late-Jurassic to early Cretaceous ages are also documented from the eastern part of the Strandja Zone, where a two-stage process of ductile and brittle deformation is documented by Elmas et al., [57], between ca. 156 and ca. 143 Ma and between ca. 136 and ca. 119 Ma, respectively. This deformation was interpreted by Elmas et al., [57] as being coincident with core complex formation during extension. However, this contradicts previous interpretations by Dabovski et al., [18] and Okay et al., [14], who connected ca. 155 Ma Rb–Sr ages as being related to regional metamorphism and thrusting during orogenesis.
The ages presented here are coeval with other regional events, such as those documented further west in the Rhodope Massif [58,59], to the north in the Carpathians [60,61,62], and to the east in Turkey and Iran [63], and they are likely connected to the Cimmerian Orogeny [4]. The Cimmerian Orogeny is mostly documented from Central Asia, where oceanic island arcs were subducted beneath Asia prior to the arrival of India and the of the Neo-Tethys basin [4]. However, recent work shows that the orogen extended much further to the west, including into Iran, Turkey and into Bulgaria [8]. In Şengör [4,64], the Cimmerian Orogeny was defined by the continent–continent collision of a Gondwana-derived continental ribbon on to Laurasia. However, more recent work indicates a series of Cimmerian tectonic events [62,65,66,67,68,69]. The Cimmerian Orogeny is succeeded by the Alpine Orogeny, which stretches further to the west, and includes the Pyrenees, Alps, Carpathians, Balkanides, Pontides, Zagros, and ultimately, the Himalayas, as the bulk of Africa and India collided with Europe and Asia, respectively. However, distinguishing between Cimmerian- and Alpine-related structures is difficult. The Sakar Unit, along with the majority of the Strandja Massif further to the east, appears to have escaped any significant Alpine orogenesis, sensu stricto [8]. Although the final stages of the Cimmerian Orogeny likely overlapped in age with the Alpine Orogeny sensu lato, the ages from this study predate the Alpine Orogeny and are not affected by it. The Alpine Orogeny is documented further to the north in the Carpathians and Srednogorie mountains [70,71] and to the west, in the Rhodope Massif [72].

6. Conclusions

1. The U–Pb ages in this study are interpreted as distinct crystallization episodes, at 149 ± 7 Ma on apatite and 114 ± 1 Ma on titanite, respectively. They are the first evidence for two stages of hydrothermal activity during the late-Jurassic to early Cretaceous in Southeastern Bulgaria and are broadly coeval with the Cimmerian orogenesis. The ages obtained in this study are broadly coeval with other units from the eastern part of the Strandja Massif, where two stages of deformation are documented.
2. The U–Pb apatite concordia systematics are robust, with both the U–Pb apatite and U–Pb titanite data implying they were unaffected by open system processes above their respective closure temperatures.
3. The study also confirms the lack of any significant Alpine orogenic footprint on the western part of the Strandja Massif, making it distinct from the neighboring Rhodope and Balkanides mountains.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-163X/10/3/266/s1. Table S1: EMPA patterns, crystals, and reference materials. Table S2: Representative electron microprobe analyses (EMPA) of albite and the number of a.p.f.u. calculated on the basis of 8 O2− of the investigated sample from the Kanarata quarry, Bulgaria. Table S3: Representative electron microprobe analyses (EMPA) of actinolite and the number of a.p.f.u. calculated on the basis of 12 O2− of the investigated sample from the Kanarata quarry, Bulgaria. Table S4: Representative electron microprobe analyses (EMPA) of chlorites and number of a.p.f.u. calculated on the basis of 12 O2− of the investigated sample from the Kanarata quarry, Bulgaria. Table S5: Representative LA-ICP-MS U–Pb apatite and titanite data from the Kanarata quarry, Bulgaria.

Author Contributions

Conceptualization, K.S., A.S., and A.P.G.; methodology, K.S., D.C., A.P.G., A.Z., and E.D.; formal analysis, K.S., D.C., and K.B.; investigation, all authors.; writing—original draft preparation, K.S., A.P.G., N.G. and A.S.; writing—review and editing, K.S., A.S., P.P. and A.G.; visualization, K.S. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

Anna Sałacińska acknowledges the financial support from an Early Career Grant (grant number EC−45065R−18) from the National Geographic Society, entitled “On the edge of two paleo-continents: in search of the eastern end of the European Variscan orogeny. David Chew acknowledges support from research grant 13/RC/2092 from Science Foundation Ireland, which is co-funded under the European Regional Development Fund and by PIPCO RSG and its member companies.

Acknowledgments

Open access for this publication was paid through funds from internal founding, at INoZ, WNP, UŚ, in Poland.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dewey, J.F.; Burke, K. Tibetan, Variscan and Precambrian basement reactivation: Products of continental collision. J. Geol. 1973, 81, 683–692. [Google Scholar] [CrossRef]
  2. Hsü, K.J.; Nachev, I.K.; Vuchev, V.T. Geologic evolution of Bulgaria in light of plate tectonics. Tectonophysics 1977, 40, 245–256. [Google Scholar] [CrossRef]
  3. Carrigan, C.W.; Mukasa, S.B.; Haydoutov, I.; Kolcheva, K. Age of Variscan magmatism from the Balkan sector of the orogen, central Bulgaria. Lithos 2005, 82, 125–147. [Google Scholar] [CrossRef]
  4. Şengör, A.M.C.; Altıner, D.; Cin, A.; Ustaömer, T.; Hsü, K.J. Origin and Assembly of the Tethyside Orogenic Collage at the Expense of Gondwana Land; Special Publications; Geological Society: London, UK, 1988; Volume 27, pp. 119–181. [Google Scholar] [CrossRef]
  5. Zanchi, A.; Berra, F.; Mattei, M.; Ghassemi, M.R.; Sabouri, J. Inversion tectonics in central Alborz, Iran. J. Struct. Geol. 2006, 28, 2023–2037. [Google Scholar] [CrossRef]
  6. Wilmsen, M.; Fürsich, F.T.; Seyed-Emami, K.; Majidifard, M.R.; Taheri, J. The Cimmerian Orogeny in northern Iran: Tectono-stratigraphic evidence from the foreland. Terra Nova 2009, 21, 211–218. [Google Scholar] [CrossRef]
  7. Angiolini, L.; Jadoul, F.; Leng, M.J.; Stephenson, M.H.; Rushton, J.; Chenery, S.; Crippa, G. How cold were the Early Permian glacial tropics? Testing sea-surface temperature using the oxygen isotope composition of rigorously screened brachiopod shells. J. Geol. Soc. 2009, 166, 933–945. [Google Scholar] [CrossRef]
  8. Cattò, S.; Cavazza, W.; Zattin, M.; Okay, A.I. No significant Alpine tectonic overprint on the Cimmerian Strandja Massif (SE Bulgaria and NW Turkey). Int. Geol. Rev. 2018, 60, 513–529. [Google Scholar] [CrossRef]
  9. Bonev, N.; Filipov, P.; Raicheva, R.; Moritz, R. Timing and tectonic significance of Paleozoic magmatism in the Sakar unit of the SakarStrandzha Zone, SE Bulgaria. Int. Geol. Rev. 2019, 61, 1957–1979. [Google Scholar] [CrossRef]
  10. Gerdjikov, I. Alpine Metamorphism and Granitoid Magmatism in the Strandja Zone: New Data from the Sakar Unit, SE Bulgaria. Turkish J. Earth Sci. 2005, 14, 167–183. [Google Scholar]
  11. Kamenov, B.K.; Vergilov, V.; Dabovski, C.; Vergiloc, I.; Ivchinova, L. The Sakar Batolithpetrology, geochemistry and magmatic evolution. Geochemistry. Mineral. Petrol. 2010, 48, 1–37. [Google Scholar]
  12. Pristavova, S.; Tzankova, N.; Gospodinov, N.; Filipov, P. Petrological study of metasomatic altered granitoids from Kanarata deposit, Sakar Mountain, southeastern Bulgaria. J. Min. Geol. Sci. 2019, 62, 53–61. [Google Scholar]
  13. Gerdjikov, I. Transformation of porphyritic granite into banded eye gneiss—An example from the NW margin of Sakar Pluton. C. R. Acad. Bulg. Sci. 2003, 56, 51–56. [Google Scholar]
  14. Okay, A.I.; Satur, M.; Tüysüz, O.; Akyüz, S.; Chen, F. The tectonics of Strandja Massif: Late-Variscan and mid-Mesozoic deformation and metamorphism in the northern Aegean. Int. J. Earth Sci. 2001, 90, 217–233. [Google Scholar] [CrossRef]
  15. Natal’in, B.; Sunal, G.; Gün, E.; Wang, B.; Zhiqing, Y. Precambrian to Early Cretaceous rocks of the Strandja Massif (northwestern Turkey): Evolution of a long lasting magmatic arc. Can. J. Earth Sci. 2016, 53, 1312–1335. [Google Scholar] [CrossRef]
  16. Bedi, Y.; Vasilev, E.; Dabovski, C.; Ergen, A.; Okuyucu, C.; Dogan, A.; Kagan Tekin, U.; Ivanova, D.; Boncheva, I.; Lakova, I.; et al. New age data from the tectonostratigraphic units of the Istranca “Massif” in NW Turkey: A correlation with SE Bulgaria. Geol. Carpath. 2013, 64, 255–277. [Google Scholar] [CrossRef] [Green Version]
  17. Dabovski, C.; Zagorchev, I. Bulgarian lands in the Alpine tectonic models of the Balkan Peninsula and Eastern Mediterranean region. In Geology of Bulgaria; Mesozoic, G., Zagorchev, I., Dabovski, C., Nikolov, T., Eds.; Prof. Marin Drinov Academic Publishing House: Sofia, Bulgaria, 2009; Volume 2, pp. 15–20. (In Bulgarian) [Google Scholar]
  18. Dabovski, C.; Boyanov, I.; Khrischev, K.; Nikolov, T.; Sapounov, L.; Yanev, Y.; Zagorchev, I. Structure and Alpine evolution of Bulgaria. Geol. Balc. 2002, 32, 9–15. [Google Scholar]
  19. Gerdjikov, I.; Gautier, P.; Cherneva, Z.; Kostopoulos, D. Tectonic setting of ultrahigh-pressure metamorphic rocks form the Chepelare area, Central Rhodope. In Proceedings of the Bulgarian Geological Society Annual Conference, Sofia, Bulgaria, 4–5 November 2003; pp. 44–45. [Google Scholar]
  20. Zagorchev, I. Alpine evolution of the pre-Alpine amphibolites facies basement in south Bulgaria. Mitt. Österr. Geol. Ges. 1993, 86, 9–23. [Google Scholar]
  21. Lilov, P. Rb-Sr and K-Ar dating of the Sakar granitoid pluton. Geol. Balc. 1990, 20, 53–60. (In Russian) [Google Scholar]
  22. Skenderov, G.; Skenderova, T. Subduction of the Vardar oceanic crust at the end of Jurassic and its role for the Alpine tectonic-magmatic development of parts of the Balkan Peninsula. Rev. Bulg. Geol. Soc. 1995, 56, 45–63. (In Russian) [Google Scholar]
  23. Ivanov, Z.; Gerdjikov, I.; Kounov, A. New data and considerations about structure and tectonic evolution of Sakar region, SE Bulgaria. Annu. Univ. Sofia Geol. Geogr. 2001, 91, 35–80. (In Bulgarian) [Google Scholar]
  24. Palshin, I.; Skenderov, G.; Bojkov, I.; Michailov, Y.; Kotov, E.; Bedrinov, I.; Ivanov, I. New geochronological data for the Cimmeridian and Alpine magmatic and hydrothermal products in the Srednogorie and Stara Planina Zones in Bulgaria. Rev. Bulg. Geol. Soc. 1989, 50, 75–92. (In Bulgarian) [Google Scholar]
  25. Skenderov, G.; Palshin, I.; Michailov, Y.; Bojkov, I.; Savova, L. On the age of the Sakar granite pluton (South-Eastern Srednogorie). Geochem. Mineral. Petrol. 1986, 22, 69–81. (In Russian) [Google Scholar]
  26. Boyadjiev, S.; Lilov, P. On the K-Ar dating of the South Bulgarian granitoids from Srednogorie and Sakar-Strandja Zones. Proceedings of the Geological Institute, ser. Geochem. Mineral. Petrogr. 1972, 26, 121–220. (In Bulgarian) [Google Scholar]
  27. Firsov, L. On the age of South-Bulgarian granites in the Rhodopes, Srednogorie and Sakar-Strandja areas. Geol. Geophys. 1975, 1, 27–34. (In Russian) [Google Scholar]
  28. Peytcheva, I.; Georgiev, S.; von Quadt, A. U/Pb ID-TIMS dating of zircons from the Sakar-Strandzha Zone: New data and old questions about the Variscan orogeny in SE Europe. In Proceedings of the Annual Conference of the Bulgaria Geological Society “Geosciences 2016”, Sofia, Bulgaria, 7–8 December 2016; Bulgarian Geological Society: Sofia, Bulgaria, 2016; pp. 71–72. [Google Scholar]
  29. Zagorchev, I.; Lilov, P.; Moorbath, S. Results of the Rubidium-Strontium and Potassium-Argon radiogeochronological studies on the metamorphic and magmatic rocks in South Bulgaria. Geol. Balc. 1989, 19, 41–54. [Google Scholar]
  30. Dencheva, S. Apatite from Sakar Mountain, Bulgaria—Morphology and physical properties. In Proceedings of the National Conference with International Participation “Geosciences 2017”, Sofia, Bulgaria, 7–8 December 2017; Bulgarian Geological Society: Sofia, Bulgaria, 2017; pp. 21–22. [Google Scholar]
  31. Chew, D.M.; Petrus, J.A.; Kamber, B.S. U-Pb LA-ICPMS dating using accessory mineral standards with variable common Pb. Chem. Geol. 2014, 363, 185–199. [Google Scholar] [CrossRef]
  32. Petrus, J.A.; Kamber, B.S. VizualAge: A Novel Approach to Laser Ablation ICP-MS U-Pb Geochronology Data Reduction. Geostand. Geoanalytical Res. 2012, 36, 247–270. [Google Scholar] [CrossRef]
  33. Paton, C.; Helistrom, J.; Paul, B.; Woodhead, J.; Herqt, J. Iolite: Freeware for the visualisation and processing of mass spectrometric data. J. Anal. At. Spectrom. 2011, 26, 2508–2518. [Google Scholar] [CrossRef]
  34. Thomson, S.N.; Gehrels, G.E.; Ruiz, J.; Buchwaldt, R. Routine low-damage apatite U–Pb dating using laser ablation-multicollector-ICPMS. Geochem. Geophys. Geosyst. 2012, 13, Q0AA21. [Google Scholar] [CrossRef]
  35. Schoene, B.; Bowring, S.A. U-Pb systematics of the McClure Mountain syenite: Thermochronological constraints on the age of the Ar-40/Ar-39 standard MMhb. Contrib. Mineral. Petrol. 2006, 151, 615–630. [Google Scholar] [CrossRef]
  36. McDowell, F.W.; McIntosh, W.C.; Farley, K.A. A precise 40Ar-39Ar reference age for the Durango apatite (U-Th)/He and fission-track dating standard. Chem. Geol. 2005, 214, 249–263. [Google Scholar] [CrossRef]
  37. Stacey, J.S.; Kramers, J.D. Approximation of terrestrial lead isotope evolution by a two-stage model. Earth Planet. Sci. Lett. 1975, 26, 207–221. [Google Scholar] [CrossRef]
  38. Spandler, C.; Hammerli, J.; Sha, P.; Hilbert-Wolf, H.; Hu, Y.; Roberts, E.; Schmitz, M. MKED1: A new titanite standard for in situ analysis of Sm-Nd isotopes and U–Pb geochronology. Chem. Geol. 2016, 425, 110–126. [Google Scholar] [CrossRef]
  39. Kennedy, A.K.; Kamo, S.L.; Nasdala, L.; Timms, N.E. Grenville skarn titanite: Potential reference material for SIMS U–Th–Pb analysis. Can. Mineral. 2010, 48, 1423–1443. [Google Scholar] [CrossRef]
  40. Aleinikoff, J.N.; Wintsch, R.; Tollo, R.P.; Unruh, D.M.; Fanning, C.M.; Schmitz, M.D. Ages and origins of rocks of the Killingworth dome, south-central Connecticut: Implications for the tectonic evolution of southern New England. Am. J. Sci. 2007, 307, 63–118. [Google Scholar] [CrossRef]
  41. Bailey, S.W. Structures of Layer Silicates: In Crystal Structures of Clay Minerals and Their X-Ray Identification; Brindley, G.W., Brown, G., Eds.; Mineralogical Society: London, UK, 1980; pp. 1–123. [Google Scholar]
  42. Cathelineau, M.; Nieva, D. A chlorite solid solution geothermometer. The Los Azufres (Mexico) geothermal system. Contrib. Mineral. Petrol. 1985, 91, 235–244. [Google Scholar] [CrossRef]
  43. Jowett, E. Fitting iron and magnesium into the hydrothermal chlorite geothermometer. In Abstract Book. In Proceedings of the GAC/MAC/SEG Joint Annual Meeting, Toronto, ON, Canada, 27–29 May 1991; Volume 16, p. A62. [Google Scholar]
  44. Kranidiotis, P.; MacLean, W.H. Systematics of chlorite alteration at the Phelps Dodge massive sulfide deposit, Matagami, Quebec. Econ. Geol. 1987, 82, 1898–1911. [Google Scholar] [CrossRef]
  45. Tilton, G.R.; Grünenfelder, M. Sphene: Uranium-lead ages. Science 1968, 159, 1458–1461. [Google Scholar] [CrossRef]
  46. Tucker, R.D.; Raheim, A.; Krogh, T.E.; Corfu, F. Uranium-lead zircon and titanite ages from the northern portion of the Western Gneiss region, south-central Norway. EPSL 1987, 81, 203–211. [Google Scholar] [CrossRef]
  47. Heaman, L.; Parrish, R. U-Pb geochronology of accessory minerals. In Applications of Radiometric Isotope Systems to Problems in Geology; Heaman, L., Ludden, J.N., Eds.; Mineralogical Association of Canada: Quebec, QC, Canada, 1991; Volume 19, pp. 59–102. [Google Scholar]
  48. Scott, D.J.; St-Onge, M.R. Constraints on Pb closure temperature in titanite based on rocks from the Ungava orogen, Canada: Implications for U-Pb geochronology and P-T-t path determinations. Geology 1995, 23, 1123–1126. [Google Scholar] [CrossRef]
  49. Watson, E.B.; Harrison, T.M.; Ryerson, F.J. Diffusion of Sm, Sr and Pb in fluoroapatite. Geochim. Cosmochim. Acta 1985, 49, 1813–1823. [Google Scholar] [CrossRef]
  50. Chemiak, D.J.; Lanford, W.A.; Ryerson, F.J. Lead diffusion in apatite and zircon using ion implantation and Rutherford backscattering techniques. Geochim. Cosmochim. Acta 1991, 55, 1663–1674. [Google Scholar]
  51. Chamberlain, K.R.; Bowring, S.A. Apatite–feldspar U–Pb thermochronometer: A reliable, mid-range (∼450 °C), diffusion-controlled system. Chem. Geol. 2000, 172, 173–200. [Google Scholar] [CrossRef]
  52. Stalder, H.A. Petrographische und mineralogische Untersu-chungen im Grimselgebiet (Mittleres Aarmassiv). Schweiz. Min. Alogischen Petrogr. Mitt. 1964, 44, 187–398. [Google Scholar]
  53. Schmid, S.M.; Fügenschuh, B.; Kissling, E.; Schuster, R. Tectonic map and overall architecture of the Alpine orogen. Eclogae Geol. Helv. 2004, 97, 93–117. [Google Scholar] [CrossRef]
  54. Mullis, J. Fluid inclusion studies during very low-grade metamorphism. In Low Temperature Metamorphism; Frey, M., Ed.; Blackie: London, UK, 1988; pp. 162–199. [Google Scholar]
  55. Mullis, J.; Dubessy, J.; Poty, B.; O’Neil, J. Fluid regimes during late stages of a continental collision: Physical, chemical, and stable isotope measurements of fluid inclusions in fissure quartz from a geotraverse through the Central Alps, Switzerland. Geochim. Cosmochim. Acta. 1994, 58, 2239–2267. [Google Scholar] [CrossRef]
  56. Taylor, R.P.; Fryer, B.J. Multiple-stage hydrothermal alteration in porphyry copper systems in northern Turkey: The temporal interplay of potassic, propylitic, and phyllic fluids. Can. J. Earth. Sci. 1980, 17, 901–926. [Google Scholar] [CrossRef]
  57. Elmas, A.; Yilmaz, I.; Yigitbas, E.; Ullrich, T. A Late Jurassic–Early Cretaceous metamorphic core complex, Strandja Massif, NW Turkey. Int. J. Earth Sci. 2011, 100, 1251–1263. [Google Scholar] [CrossRef]
  58. Stampfli, G.; Kozur, H. Europe from Variscan to the Alpine cycles. In European Lithosphere Dynamics, 32th ed.; Geological Society of London: London, UK, 2006; Volume 32, pp. 57–82. [Google Scholar]
  59. Stampfli, G.; Mosar, J.; Favre, P.; Pillevuit, A.; Vannay, J.-C. Permo-Mesozoic evolution of the western Tethyan realm: The Neotethys/EastMediterranean connection. In PeriTethys Memoir 6: Peritethyan Rift/Wrench Basins and Passive Margins; Me ´Moires du Museum Nationale d’Histoire; Ziegler, P., Cavazza, W., Robertson, A., Crasquinsoleau, S., Eds.; IGCP 369: Lisbon, Portugal, 2001; Volume 186, pp. 51–108. [Google Scholar]
  60. Golonka, J.; Krobicki, M.; Oszczypko, N.; Ślaczka, A.; Słomka, T. Geodynamic evolution and palaeogeoraphy of the Polish Carpathians and adjacent areas during Neo-Cimmerian and preceding events (latest Triassic-earliest cretaceous). Geol. Soc. Spec. Publ. 2003, 208, 138–158. [Google Scholar] [CrossRef]
  61. Szopa, K.; Włodyka, R.; Chew, D. LA-ICP-MS U-Pb apatite dating of Lower Cretaceous rocks from teschenite-picrite association in the Silesian Unit (southern Poland). Geol. Carpathica 2014, 65, 273–284. [Google Scholar] [CrossRef] [Green Version]
  62. Gawęda, A.; Szopa, K.; Chew, D.; O’Sullivan, G.J.; Burda, J.; Klötzli, U.; Golonka, J. Variscan post-collisional cooling and uplift of the Tatra Mountains crystalline block constrained by integrated zircon, apatite and titanite LA-(MC)-ICP-MS U-Pb dating and rare earth element analyses. Chem. Geol. 2018, 484, 191–209. [Google Scholar] [CrossRef] [Green Version]
  63. Ballato, P.; Nowaczyk, N.R.; Landgraf, A.; Strecker, M.R.; Friedrich, A.; Tabatabaei, S.H. Tectonic control on sedimentary facies pattern and sediment accumulation rates in the Miocene foreland basin of the southern Alborz mountains, northern Iran. Tectonics 2008, 27, TC6001. [Google Scholar] [CrossRef] [Green Version]
  64. Şengör, A.M.C. Mid-Mesozoic closure of Permo-Triassic Tethys and its implications. Nature 1979, 279, 590–593. [Google Scholar] [CrossRef]
  65. Stampfli, G.; Borel, G. The TRANSMED transects in space and time: Constraints on the Paleotectonic evolution of the Mediterranean Domain. In The TRANSMED Atlas: The Mediterranean Region from Crust to Mantle; Cavazza, W., Roure, F., Spakman, W., Ziegler, P., Eds.; Springer: Berlin, Germany, 2004; pp. 53–80. [Google Scholar]
  66. Stampfli, G.; Hochard, C. Plate tectonics of the Alpine realm. Geol. Soc. London. In Ancient Orogens and Modern Analogues; Special Publications; The Geological Society of London: London, UK, 2009; Volume 327, pp. 89–111. [Google Scholar]
  67. Okay, A.I.; Sunal, G.; Sherlock, S.; Altiner, D.; Tüysüz, O.; Kylander-Clark, A.R.C.; Aygül, M. Early Cretaceous sedimentation and orogeny on the southern active margin of Eurasia: Central Pontides, Turkey. Tectonics 2013, 32, 1247–1271. [Google Scholar] [CrossRef]
  68. Okay, A.I.; Altiner, D.; Kiliç, A.M. Triassic limestone, turbidites and serpentinite–the Cimmeride orogeny in the Central Pontides. Geol. Mag. 2015, 152, 460–479. [Google Scholar] [CrossRef] [Green Version]
  69. Topuz, G.; Göçmengil, G.; Rolland, Y.; Çelic, Ö.F.; Zack, T.; Schmitt, A.K. Jurassic accretionary complex and ophiolite from northeast Turkey: No evidence for the Cimmerian continental ribbon. Geology 2013, 41, 255–258. [Google Scholar] [CrossRef] [Green Version]
  70. Hrubcová, P.; Środa, P.; Grad, M.; Geissler, W.H.; Guterch, A.; Vozár, J.; Hegedűs, E. From the Variscan to the Alpine Orogeny: Crustal structure of the Bohemian Massif and the Western Carpathians in the light of the SUDETES 2003 seismic data 2010. Geophys. J. Int. 2010, 183, 611–633. [Google Scholar] [CrossRef] [Green Version]
  71. Georgiev, S.; von Quadt, A.; Heinrich, C.A.; Peytcheva, I.; Marchev, P. Time evolution of a rifted continental arc: Integrated ID-TIMS and LA-ICPMS study of magmatic zircons from the Eastern Srednogorie, Bulgaria. Lithos 2012, 154, 53–67. [Google Scholar] [CrossRef] [Green Version]
  72. Lips, A.; White, S.; Wijbrans, J. Middle-Late Alpine thermotectonic evolution of the southern Rhodope Massif, Greece. Geodin. Acta 2000, 13, 281–292. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The Sakar Mountains in Bulgaria—study area: (A) geological map of the Strandja Massif (Southeast Bulgaria) and neighboring regions; (B) geological sketch map of the Sakar Unit (modified after [10,11]); (C) Kanarata quarry—sample locality.
Figure 1. The Sakar Mountains in Bulgaria—study area: (A) geological map of the Strandja Massif (Southeast Bulgaria) and neighboring regions; (B) geological sketch map of the Sakar Unit (modified after [10,11]); (C) Kanarata quarry—sample locality.
Minerals 10 00266 g001
Figure 2. Example of minerals within the investigated veins. (A) Blue apatite crystal on albite and chlorite, (B) green titanite, (C) albite crystals, (D) actinolite overgrown by albite, (E) epidote aggregate, (F) molybdenite flake on milky quartz, and (G) rutile crystals in quartz. Scale bars are 1 cm.
Figure 2. Example of minerals within the investigated veins. (A) Blue apatite crystal on albite and chlorite, (B) green titanite, (C) albite crystals, (D) actinolite overgrown by albite, (E) epidote aggregate, (F) molybdenite flake on milky quartz, and (G) rutile crystals in quartz. Scale bars are 1 cm.
Minerals 10 00266 g002
Figure 3. Photographs of the most common rock types in the investigated Kanarata quarry. (A) monzonite with an albite nest, (B) amphibolite xenolith within the contact zone, (C) an example of Alpine-type fissures, and (D) the second vein type with massive quartz and epidote. Ap = apatite, Chl = chlorite, Ab = albite, Ep = epidote, and Qtz = quartz.
Figure 3. Photographs of the most common rock types in the investigated Kanarata quarry. (A) monzonite with an albite nest, (B) amphibolite xenolith within the contact zone, (C) an example of Alpine-type fissures, and (D) the second vein type with massive quartz and epidote. Ap = apatite, Chl = chlorite, Ab = albite, Ep = epidote, and Qtz = quartz.
Minerals 10 00266 g003
Figure 4. Tera–Wasserburg concordia plots for LA–ICP-MS U–Pb titanite (A) and apatite (B) analyses from the Kanarata quarry. Data-point error ellipses are 2σ.
Figure 4. Tera–Wasserburg concordia plots for LA–ICP-MS U–Pb titanite (A) and apatite (B) analyses from the Kanarata quarry. Data-point error ellipses are 2σ.
Minerals 10 00266 g004
Figure 5. AlIV vs. temperature plot for the chlorites investigated in this study. Temperatures is based on three geothermometers: [42] (red), [43] (green), and [44] (orange).
Figure 5. AlIV vs. temperature plot for the chlorites investigated in this study. Temperatures is based on three geothermometers: [42] (red), [43] (green), and [44] (orange).
Minerals 10 00266 g005
Figure 6. Simplified model of the hydrothermal activity in the Sakar Unit of Southeastern Bulgaria, based on the geochronological data presented in this study and field observations. (A) Formation of the Variscan host rock (~305–295 Ma). (B) The first, older hydrothermal event forms the primary mineralization (ca. 149 Ma. (C) The second, younger hydrothermal episode forms titanite (ca. 114 Ma age).
Figure 6. Simplified model of the hydrothermal activity in the Sakar Unit of Southeastern Bulgaria, based on the geochronological data presented in this study and field observations. (A) Formation of the Variscan host rock (~305–295 Ma). (B) The first, older hydrothermal event forms the primary mineralization (ca. 149 Ma. (C) The second, younger hydrothermal episode forms titanite (ca. 114 Ma age).
Minerals 10 00266 g006
Table 1. Representative electron microprobe analyses (EMPA) of apatite and number of a.p.f.u. are calculated on the basis of 25 O2− of the investigated sample from the Kanarata quarry, Bulgaria.
Table 1. Representative electron microprobe analyses (EMPA) of apatite and number of a.p.f.u. are calculated on the basis of 25 O2− of the investigated sample from the Kanarata quarry, Bulgaria.
Compound12345678910
CaO (wt.%)55.6955.8255.5755.5755.8355.7755.7556.1356.1156.26
SiO20.230.10b.d.l.b.d.l.b.d.l.b.d.l.b.d.l.b.d.l.b.d.l.b.d.l.
P2O541.5142.3241.8342.5642.5842.3141.8742.6142.4242.07
SO30.120.13b.d.l.b.d.l.b.d.l.b.d.l.b.d.l.b.d.l.b.d.l.b.d.l.
F2.832.823.052.552.272.881.992.792.642.40
H2O0.390.410.270.540.680.370.800.430.500.61
O=F,Cl1.191.191.291.070.961.220.841.181.111.01
Total99.58100.4299.44100.14100.41100.1199.57100.79100.56100.33
Ca (a.p.f.u.)10.3710.2710.3610.2010.2110.3110.2610.2910.3010.34
Si0.040.02--------
P6.116.156.166.186.156.186.096.176.156.11
S0.020.02--------
XF1.551.531.681.381.221.571.081.511.431.30
XOH0.450.470.320.620.780.430.920.490.570.70
Note: b.d.l.—below detection limit.
Table 2. Representative electron microprobe analyses (EMPA) of titanite and number of ion a.p.f.u. calculated on the basis of 18 O2− of the investigated sample from the Kanarata quarry, Bulgaria.
Table 2. Representative electron microprobe analyses (EMPA) of titanite and number of ion a.p.f.u. calculated on the basis of 18 O2− of the investigated sample from the Kanarata quarry, Bulgaria.
Compound12345678910
SiO2 (wt.%)30.2030.1929.9330.1830.0430.1230.3130.1930.2329.91
TiO238.6538.4738.4736.8838.3838.8138.2138.5438.5438.59
CaO28.3328.2628.0728.2428.2528.4128.2428.2428.2628.21
Al2O31.111.231.021.750.950.871.181.030.840.86
MnOb.d.l.b.d.l.b.d.l.b.d.l.b.d.l.b.d.l.0.120.12b.d.l.b.d.l.
Fe2O30.680.600.760.640.610.740.710.340.760.43
Y2O30.130.160.110.160.17b.d.l.0.140.15 0.14
V2O50.190.250.270.260.300.240.290.190.310.24
F0.240.180.100.410.120.020.220.160.260.10
O=F0.100.080.040.170.050.010.090.070.110.04
Total99.4299.2598.6998.3598.7899.2199.3298.8899.0998.44
Si (a.p.f.u.)3.583.583.573.623.583.573.603.603.603.58
Ti3.453.433.453.333.443.463.413.453.453.47
Ca3.603.593.583.633.613.613.593.603.603.61
Al0.160.170.140.250.130.120.170.140.120.12
Mn------0.010.01--
Fe0.060.050.070.060.060.070.060.030.070.04
Y0.010.010.010.010.01 0.010.01 0.01
V0.020.020.020.020.020.020.020.020.020.02
F0.090.070.040.160.040.010.080.060.100.04
Total10.9610.9210.8711.0810.9010.8510.9610.9210.9510.88
Note: b.d.l.—below detection limit.

Share and Cite

MDPI and ACS Style

Szopa, K.; Sałacińska, A.; Gumsley, A.P.; Chew, D.; Petrov, P.; Gawȩda, A.; Zagórska, A.; Deput, E.; Gospodinov, N.; Banasik, K. Two-Stage Late Jurassic to Early Cretaceous Hydrothermal Activity in the Sakar Unit of Southeastern Bulgaria. Minerals 2020, 10, 266. https://doi.org/10.3390/min10030266

AMA Style

Szopa K, Sałacińska A, Gumsley AP, Chew D, Petrov P, Gawȩda A, Zagórska A, Deput E, Gospodinov N, Banasik K. Two-Stage Late Jurassic to Early Cretaceous Hydrothermal Activity in the Sakar Unit of Southeastern Bulgaria. Minerals. 2020; 10(3):266. https://doi.org/10.3390/min10030266

Chicago/Turabian Style

Szopa, Krzysztof, Anna Sałacińska, Ashley P. Gumsley, David Chew, Petko Petrov, Aleksandra Gawȩda, Anna Zagórska, Ewa Deput, Nikolay Gospodinov, and Kamila Banasik. 2020. "Two-Stage Late Jurassic to Early Cretaceous Hydrothermal Activity in the Sakar Unit of Southeastern Bulgaria" Minerals 10, no. 3: 266. https://doi.org/10.3390/min10030266

APA Style

Szopa, K., Sałacińska, A., Gumsley, A. P., Chew, D., Petrov, P., Gawȩda, A., Zagórska, A., Deput, E., Gospodinov, N., & Banasik, K. (2020). Two-Stage Late Jurassic to Early Cretaceous Hydrothermal Activity in the Sakar Unit of Southeastern Bulgaria. Minerals, 10(3), 266. https://doi.org/10.3390/min10030266

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop