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

: 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 associatedwiththeenigmaticCimmerianOrogenicepisode,butlimitedgeochronologyhasbeenundertaken 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.


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  [10,11]); (C) Kanarata quarry-sample locality.

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).

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.

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.

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.

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 cm 3 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

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.

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.

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 N 2 (~7 mL/min) in a~1 cm 3 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): 31 P (10), 35 Cl (20) 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 207 Pb-based correction method.

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 (Fe 2+ 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 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 (Mg 5 Al)(Si 3 Al)O 10 (OH) end-member, 12%-49% of the chamosite (Fe 2+ 5 Al)(Si 3 Al)O 10 (OH) end-member, and <3% of the pennantite (Mn 5 Al)(Si 3 Al)O 10 (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.

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 207 Pb/ 206 Pb 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 207 Pb/ 206 Pb initial ratio of 0.772 ± 0.021 ( Figure 4B). All LA-ICP-MS data are listed in Table S5.

Age Significance: Apatite vs. Titanite
Titanite (CaTiSiO 5 ) and apatite (Ca 5 (PO 4 ) 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 207 Pb/ 206 Pb initial ratio of the titanite (suggesting derivation from a radiogenic fluid), also implies separate crystallization episodes.

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]). Multiequilibrium 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 207 Pb/ 206 Pb 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). 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). Temperatures is based on three geothermometers: [42] (red), [43] (green), and [44] (orange).

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. Figure 5. Al IV vs. temperature plot for the chlorites investigated in this study. Temperatures is based on three geothermometers: [42] (red), [43] (green), and [44] (orange).

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].

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 http://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 O 2− 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 O 2− 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 O 2− 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.
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.