Metamorphic and Metasomatic Kyanite-Bearing Mineral Assemblages of Thassos Island (Rhodope, Greece)

The Trikorfo area (Thassos Island, Rhodope massif, Northern Greece) represents a unique mineralogical locality with Mn-rich minerals including kyanite, andalusite, garnet and epidote. Their vivid colors and large crystal size make them good indicators of gem-quality materials, although crystals found up to now are too fractured to be considered as marketable gems. The dominant lithology is represented by a garnet–kyanite–biotite–hematite–plagioclase ± staurolite ± sillimanite paragneiss. Thermodynamic Perple_X modeling indicates conditions of ca. 630–710 ◦C and 7.8–10.4 kbars. Post-metamorphic metasomatic silicate and calc-silicate (Mn-rich)-minerals are found within (i) green-red horizons with a mineralogical zonation from diopside, hornblende, epidote and grossular, (ii) mica schists containing spessartine, kyanite, andalusite and piemontite, and (iii) weakly deformed quartz-feldspar coarse-grained veins with kyanite at the interface with the metamorphic gneiss. The transition towards brittle conditions is shown by Alpine-type tension gashes, including spessartine–epidote–clinochlore–hornblende-quartz veins, cross-cutting the metamorphic foliation. Kyanite is of particular interest because it is present in the metamorphic paragenesis and locally in metasomatic assemblages with a large variety of colors (zoned blue to green/yellow-transparent and orange). Element analyses and UV-near infrared spectroscopy analyses indicate that the variation in color is due to a combination of Ti4+–Fe2+, Fe3+ and Mn3+ substitutions with Al3+. Structural and mineralogical observations point to a two-stage evolution of the Trikorfo area, where post-metamorphic hydrothermal fluid circulation lead locally to metasomatic reactions from ductile to brittle conditions during Miocene exhumation of the high-grade host-rocks. The large variety of mineral compositions and assemblages points to a local control of the mineralogy and fO2 conditions during metasomatic reactions and interactions between hydrothermal active fluids and surrounding rocks.


Introduction
Dispersed metal ions in substitution for Al and Si may lead to a large variety of colors in calc-and alumino-silicate minerals, e.g., [1,2]. For instance, the most common natural color of kyanite is blue but it can also be green, yellow, orange, white, black, grey or colorless as a function of the nature of elemental substitutions [3,4]. Kyanite is composed of usually >98 wt. % Al 2 O 3 and SiO 2 . The~2 wt. % left is generally dominated by Fe, Ti, Mg, Mn, Cr and V. The variation in blue colors is attributed to Ti and/or Fe substitution and intervalence charge transfers within the crystal lattice [5], while the orange color of kyanite is generally attributed to the presence of Mn 3+ [1,[6][7][8][9]. The incorporation of trace elements is interpreted as the result of (i) oxygen fugacity during crystal growth (Fe), (ii) metamorphic grade, temperature of formation (Ti), or the nature of the protolith (Cr and V) [4]. In the same way, Mn 3+ may substitute for Al 3+ in andalusite [10,11] resulting in the dark green Mn-rich variety of andalusite, formerly known as viridine [2,12]. This may also be the case for epidotes e.g., [13], garnets e.g., [14] and other metamorphic and metasomatic Ca-Al silicates. Among the interesting trace elements to produce gem quality calc-and alumino-silicate minerals, Mn seems to play an important role. High-PT manganian silicate and calc-silicate metamorphic assemblages are very rare on Earth and require very specific conditions such as the preservation of high f O 2 during the entire metamorphic cycle, e.g., [15,16]. The presence in the protolith of Mn-oxides capable of buffering f O 2 to high levels during metamorphism as well as low fluid/rock ratio during the whole geodynamic evolution are necessary [17,18]. The best way to produce gem-quality Mn-rich silicates seems thus hydrothermal fluid circulation through an already Mn-rich protolith.
The geology of Greece is marked by a recent alpine belt, which results in outstanding localities for minerals and gemstones owing to high PT metamorphic conditions and subsequent magmatism and hydrothermalism, e.g., [19][20][21][22][23]. Among these localities, the island of Thassos, in the Rhodope Belt, represents a unique mineralogical locality for kyanite and uncommon varieties of manganiferous/manganian and magnesian silicate and calc-silicate minerals, including garnet, andalusite and epidote among others [19,[24][25][26][27][28][29][30][31]. These crystals are often large (several centimeters long) and vividly colored, translucent but commonly fractured, which makes them usually inappropriate for faceting gemstones, but more proper for cabochon-shaped material. However, considering the extent of the geological units, their occurrence is promising for gemstones exploration in the area. In this article, we use the term "near gem-quality" as a reference to their aesthetic aspect (large size and vivid color) rather than to their potential as facetable material.
Thassos Island is part of the Southern Rhodope Metamorphic Core Complex (SRCC), with PT conditions recorded by the garnet-kyanite-biotite-hematite-plagioclase ± staurolite ± sillimanite assemblage of the paragneisses of the Trikorfo area (intermediate unit) in the range 600-650 • C, 4-7 kbars [24]. At the Trikorfo area, four main distinct lithologies, often bearing Ca-, Mg-and Mn-rich silicates, are distinguished with (i) the dominant metamorphic lithology composed of metasedimentary rocks, mostly mica schists and paragneisses, with intercalation of carbonate layers especially towards the base and top of the formation [24], (ii) green to red horizons parallel to the regional foliation enriched in silicate and calc-silicate minerals always showing a typical mineralogical zonation from the contact with the paragneisses/mica schists to the center of the layers, (iii) kyanite-andalusite bearing mica schists, and (iv) kyanite-quartz ± feldspar ± andalusite weakly deformed coarse-grained veins. The transition to brittle deformation conditions is shown by quartz-clinochlore-ilmenite-adularia-albite alpine-type tension gashes cross-cutting the metamorphic foliation [25]. Locally, when cross-cutting Mn-rich layers, these brittle veins can host spessartine and Mn-bearing clinochlore, epidote and hornblende.
On the basis of the description made by Voudouris et al. [20], the aim of the present study is (i) to obtain a complete inventory of the metamorphic and metasomatic minerals and their assemblages found at the Trikorfo site, (ii) to evaluate the cation substitution responsible for the color variations, (iii) to reevaluate the PT conditions attained by the gneiss host-rock on the basis of a comprehensive petrological phase equilibria using the Perple_X program, and (iv) to discuss the conditions of formation, metamorphic versus metasomatic, of these exceptional parageneses and assemblages of (Mn-rich)-silicate and calc-silicate minerals. Kyanite, as a ubiquitous silicate found in both metamorphic and localized metasomatic assemblages, is the referent mineral described in the manuscript.

The Southern Rhodope Core Complex and Thassos Island
The Aegean domain, which extends from the Rhodope Massif in the north to the island of Crete in the south of the Hellenides, is a broad metamorphic domain whose formation started with the subduction of the Adriatic microplate (the northern part of the Africa Plate) below the Eurasia Plate, in late Cretaceous to Eocene times, e.g., [32][33][34][35]. The progressive southward retreat of the downgoing slab since the Oligocene is responsible for the exhumation of metamorphic core complexes (MCC) along low-angle detachment faults [36]. The Southern Rhodope Core Complex (SRCC) [37][38][39][40][41] and the Attic-Cycladic Metamorphic Complex (ACMC) [42][43][44] are the two main metamorphic domains resulting from Aegean extension during the Oligocene to Miocene. The Aegean domain is now in a back-arc position and the subduction front is currently located to the south of the island of Crete [35]. The SRCC, located in the northern part of Greece, displays an almost triangular shape and extends to south to southern Bulgaria [40]. It is mainly constituted by the Vertiskos Gneiss Complex from the Serbo-Macedonian Massif and to the west by the Chalkidiki Peninsula dominated by gneisses from the Sideroneron Massif ( Figure 1). The area is marked by Jurassic and Cretaceous age ultra-high-pressure metamorphic events, including kyanite eclogites, followed by Barrovian type metamorphism, e.g., [37][38][39][40][41]. Geological map and regional tectonic framework of northern Greece, modified after Melfos and Voudouris [45].
Thassos Island, with an area of 400 km 2 , exposes the south-east end of the SRCC. The geology of the island is dominated by the juxtaposition of shallow-dipping units mainly made of marbles, orthogneisses and paragneisses/mica schists ( Figure 2). It is assumed that orthogneisses are the result of metamorphism of Hercynian plutonic units of the Greek Rhodope Massif [36,46] and that the protoliths of marbles and mica schists correspond to the Mesozoic sedimentary cover overlying the Hercynian plutonic basement.  [36] and Brun and Sokoutis [40]. The location of the Trikorfo area is indicated with an orange star. Thassos Island is represented by three main metamorphic units with, from base to top, (i) the lower metamorphic unit dominated by mylonitic marbles with intercalations of mica schists and paragneisses forming metamorphic domes to the east of the island, (ii) an intermediate metamorphic unit composed at the base of gneisses overlapped by marbles with paragneiss, mica schist and amphibolite intercalations, (iii) the upper unit is considered as the hanging wall, consisting of non-mylonitic metamorphic rocks dominated by gneisses, migmatites, pegmatites and locally ultramafic rocks and marbles [36,40]. Metamorphic rocks from the lower and intermediate units reached amphibolite facies conditions [24,40,47,48], achieved during crustal thickening [49] at 21-18 Ma and 26-23 Ma (Rb-Sr mica ages from the main metamorphic fabric) for the lower and intermediate units, respectively [36].
Retrograde metamorphism is marked by the development of mylonite marbles, common boudinage along a shallow-dipping foliation with NE-SW stretching lineation with a top to SW sense of shear, and widespread brittle tension gashes [36,40,48,50]. Rb-Sr (on white mica and biotite) and U-Pb (on monazite and xenotime) dating by the above authors revealed a continuous cooling history under ductile to brittle conditions from 700 • C to 300 • C from 26 Ma to 12 Ma. The exhumation PT path is interpreted as the result of the progressive uplift of the metamorphic rocks during the formation of the metamorphic core complex [36,51]. The emplacement of the low-angle detachment fault in the Oligocene-Miocene is also associated with the intrusion of syntectonic plutons at the SRCC scale, where the Symvolon pluton crops out 30 km to the NW of Thassos Island [52][53][54] (Figure 1).

The Trikorfo Area
The study area at Trikorfo is located along the major SW-dipping low-angle normal detachment separating the lower from the intermediate unit of Thassos Island. Mica schists and paragneisses predominate and tectonically overlay the marbles of Theologos (i.e. lower unit). Muscovite, phlogopite, or two mica quartzo-feldspathic schists (mica schists or paragneisses) with plagioclase (oligoclase) porphyroblasts are the most common lithologies, locally interbedded with carbonate/dolomitic rocks towards the base and top of the formation [24]. In many occurrences, these rocks are rich in minerals from the epidote and garnet groups together with common tourmaline (schorl mainly) and opaque minerals, mostly hematite and ilmenite. Metamorphic Al-silicate minerals (kyanite, fibrolitic sillimanite, andalusite, staurolite) are also documented. These mica schists are thus likely derived from calcareous pelitic sediments [24]. Interestingly, this unit also locally presents uncommon Mn-rich mineral assemblages attributed to post-metamorphic hydrothermal recrystallization [20,24,26]. The garnet-kyanite-biotite-hematite-plagioclase ± staurolite ± sillimanite paragenesis of the mica schists, contemporaneous with the ductile fabric related to extensional deformation, has been specifically interpreted as having attained PT conditions of 5.5 ± 1.5 kbars and 600 ± 50 • C, without any evidence of partial melting [24]. Voudouris et al. [20] documented at the top of the formation the presence of an orthogneiss at the contact with the paragneisses/mica schists. Although not exposed, a granitoid is considered to be genetically related to the widespread Miocene gold mineralization in the area [27].

Materials and Methods
Representative samples of host-rocks, mineral assemblages and veins were selected for bulk rock analyses and conventional 30 µm polished sections for petrographic observations and electron probe microanalyses (EPMA). Individual kyanite crystals were prepared for bulk rock analyses (powder) and for EPMA and Laser Ablation-Inductively Coupled Plasma Mass Spectroscopy (LA-ICPMS) on 200-300 µm thick sections, and UV-visible-near infrared absorption spectroscopy on 1 mm thick sections. The metamorphic paragenesis of the dominant gneiss of the area was modeled by Perple_X for PT conditions estimation. This work represents a compilation of samples collected by different authors during successive field campaigns, resulting in a large variety of sample labelling. All samples come from the Trikorfo area (intermediate unit of Thassos Island) ( Figure 2).

Bulk-Rock Analyses
Host-rock samples CS16_292b and CS16_297 and individual kyanite crystals were powderized so as to obtain a grain-size of about 80 µm. Kyanite crystals were cut as to not contain any apparent alteration minerals, mainly micas, and cleaned with deionized water. Bulk compositions (oxides and trace elements) were obtained by ICP-Optical Emission Spectrometry and ICP-Mass Spectrometry (LiBO 2 fusion) at the CRPG-CNRS laboratory (Nancy, France). Sample preparation, analytical conditions and limits of detection are given in Carignan et al. [55].

Electron Probe MicroAnalyses (EPMA)
The chemical composition of minerals (samples CS16_292b, CS297, KT01-03, KT07 and MB13) was determined using a Cameca SX100 electron microprobe analyzer (EPMA) equipped only with wavelength dispersive spectrometers at the GeoRessources laboratory (Nancy, France) operating with an emission current of 20 nA, an acceleration voltage of 15 kV, and a beam diameter of 1 µm.

Laser Ablation-Inductively Coupled Plasma Mass Spectroscopy (LA-ICPMS)
Kyanite in-situ quantitative LA-ICPMS elemental analyses were performed at GeoRessources laboratory (Nancy, France) using an Agilent 7500c quadrupole ICPMS coupled with a 193 nm GeoLas ArF Excimer laser (MicroLas, Göttingen, Germany). Laser ablation was performed continuously with a speed of 2 µm/s with a constant 5 Hz pulse frequency and a constant fluence of 4.20 J/cm 2 by focussing the beam at the sample surface, from edge to edge of the crystals, along the same profiles analyzed by EPMA. Diameter of the ablation spot was 32 µm. Helium was used as carrier gas to transport the laser-generated particles from the ablation cell to the ICPMS and argon was added as an auxiliary gas via a flow adapter before the ICP torch. Typical flow rates of 0.5 L/min for He and 1 L/min for Ar were used. EPMA analyses of Si and Al served as internal standards. The certified reference materials NIST SRM 610 and 612 (concentrations from Jochum et al. [57]) were used as external standards for calibration of all analyses; they were analyzed twice at the beginning and at the end for each kyanite sample, following a bracketing standardization procedure. 27 Al and 29 Si were measured with a dwell-time of 10 ms; the following 17 isotopes were measured with a dwell-time of 20 ms for each: 7 Li (limit of detection of 20 ppm), 23 Na (70 ppm), 24 Mg (1 ppm), 39 K (50 ppm), 43 Ca (6000 ppm), 47 Ti (20 ppm), 51 V (1 ppm), 53 Cr (15 ppm), 55 Mn (15 ppm), 57 Fe (150 ppm), 63 Cu (10 ppm), 66 Zn (10 ppm), 88 Sr (0.5 ppm), 107 Ag (2 ppm), 137 Ba (3 ppm), 197 Au (1 ppm), and 208 Pb (0.5 ppm). All data were reduced off-line with the limits of detection calculated using the commercial version of Iolite (Version 3.71, https://iolite-software.com) data reduction software [58] running with Igor Pro. The uncertainty is calculated for each element for each analysis as a function of the uncertainty of the reference standards, the absolute element content and the time of integration [58]. The resulting value of uncertainty is then varying for a single element between each analysis and is in the order of 10% for high concentration to 80% when the element concentration is measured close to the limit of detection.

UV-Visible-Near Infrared Absorption Spectra
UV-visible-near infrared absorption spectra were measured at room temperature using a Perkin-Elmer 1050 dual beam spectrophotometer at the Institut des Matériaux Jean Rouxel (University of Nantes, France) on blue-green (sample THA16), yellow (THA08), and orange (THA13) kyanite crystals. Spectra were acquired at a resolution of 1 nm at a rate of 120 nm/min. As our samples were significantly fractured and hence translucent but not transparent, we prepared them as 1 mm thick slices. At this thickness, they became sufficiently transparent so that the light beam could go through, and they also remained colored enough to characterize spectral features related to color.

Phase Diagram Calculation
PT conditions of the main paragneiss host-rock (sample CS16_292b) were modeled on the basis of the paragenesis described in Section 4.2. The phase diagram was calculated using Perple_X Version 6.8.5 (http://www.perplex.ethz.ch) [59] and the internally consistent end-member data set of Holland and Powell [60]. Calculations were undertaken in a set of chemical systems ranging from NCKFMASH to MnNCKFMASHT(O) (MnO-Na 2 O-CaO-K 2 O-FeO-MgO-Al 2 O 3 -SiO 2 -H 2 O-TiO 2 -O 2 ) in order to constrain the effect of TiO 2 , Fe 2 O 3 and MnO-bearing minerals.

Tectonic and Structural Setting of Al-Silicates
The main lithology of the Trikorfo area is dominated by garnet-kyanite-biotite-hematiteplagioclase ± staurolite ± sillimanite bearing paragneisses/mica schists. At the scale of Thassos Island, the orientation of the main foliation varies widely due to post-foliation tectonic events, including large-scale folding. In the Trikorfo area, rare relics of an early foliation S n−1 are observed. The paragneisses/mica schists mainly show a sub-horizontal (azimuth N140-300, dip up to 20 • W) S n foliation, bearing kyanite mineral lineation, oriented N255 ± 15 • , and gently dipping to the WSW. Late brittle deformation is evidenced by tension gashes that cross-cut S n . These veins are oriented N330 ± 15 • and are generally close to vertical. The variation in the orientation of the tension gashes is the same as the regional lineation underlined by synmetamorphic kyanite crystals within 30 • that may reflect their continuous formation during counter-clockwise rotation of Thassos Island [40].

Lithological Units and Rock Sampling
Samples CS16_292b (Figure 3a) and CS16_292a represent the main paragneiss/schist host-rock, apparently unaltered at the outcrop scale. An increasing number of green layers parallel to the gneiss foliation are observed towards the top of the formation (samples CS16_297 (Figure 3b) and Th06). Although not analyzed in our samples, the presence of braunite in the mica schists of the Trikorfo area was shown by previous studies [26,28]. Green to red calc-silicate rocks (up to 80 cm wide) are observed, at the top of the formation, interbedded or at the contacts between quartzo-feldspathic mica schists and metacarbonate rocks ( Figure 3c). The rocks typically show a mineral zonation with tremolitic hornblende, diopside, quartz, anorthite/bytownite, epidote, grossular, titanite and Mn-bearing clinozoisite from an outer band adjacent to host mica schists towards the center of the layer (samples Th01, Th02, Th03, KT01 and KT02 (Figure 3c-e)). Sample Th04 is a spessartine-epidote-braunite-hornblende mica schist. Kyanite is found as a metamorphic mineral, underlining the mineral lineation of the main paragneisses/mica schists rocks as in sample CS16_292b. Locally at the top of the formation, kyanite-piemontite ± andalusite ± muscovite ± spessartine ± braunite mica schists (where kyanite is distributed without any specific orientation) are intercalated within the spessartine-kyanite-bearing mica schists (samples Th05, KT03 (Figure 4a), KT05, and KT07 (Figure 3f)). Towards the top of the formation, pluricentimetric kyanite crystals are observed at the contact between metamorphic quartzo-feldspathic kyanite-bearing rocks and coarse-grained quartz-feldspar veins, which are transposed within or secant to the main foliation. In this case, kyanite is generally weakly oriented parallel to the regional lineation. In general the crystals are zoned with dark blue cores and green-transparent (Figure 4b,c) to yellow rims ( Figure 4d). Less common are orange kyanite-bearing samples showing association with green andalusite, spessartine and piemontite (Figure 4e,f). An orthogneiss exposed in the area is composed of orthoclase, oligoclase/andesine along with phlogopite, quartz, hematite, epidote and allanite and is surrounded by a transitional zone with the same mineralogy as the above mentioned calc-silicate layers (i.e. amphibole, plagioclase, and epidote mineral groups from granitoid outwards) [20]. The presence of pluricentimetric crystals of tourmaline (schorl mainly) was also shown in quartz veins and late brittle fissures [20].

Mineralogy and Mineral Chemistry of Paragneisses/Mica Schists of the Metamorphic Unit at Trikorfo
The composition of two host-rock samples (CS16_292a and _292b) is shown in Table S1a. The structure of the rock is granolepidoblastic with quartz-plagioclase-dominated layers and a schistosity underlined by biotite locally retrogressed to chlorite. C'/S structures indicate top N228 sense of shear. The mineralogy is made of biotite, garnet, quartz, kyanite, staurolite, plagioclase, apatite and Ti-hematite ( Figure 6a). Rare late muscovite and epidote are observed without any preferential orientation ( Figure 6a). Biotite has the composition of phlogopite with X Mg in the range 0.73-0.75 (Table S2a). Plagioclase is oligoclase/andesine with a composition between An 23 and An 34 (Table S3a). Garnet and staurolite crystals are isolated in the mica-rich layers. Garnet grains are always corroded remnants of preexisting porphyroblasts (Figure 6a). Garnet chemistry (Table S4a) is dominated by X Fe (up to 39.25) and X Mn (up to 44.06); X Mg is significant, up to 18.82, and X Ca is never higher than 7.72. Garnet composition thus lies close to the border between almandine-(pyrope) and spessartine fields (Figure 7a). No significant and systematic chemical variation was noticed from border to border among all analyzed crystals (Table S4a). Staurolite, also corroded, contains significant amounts of MnO, up to 1.47 wt. % (Table S5). Although Dimitriadis [24] indicated that kyanite seems to postdate garnet growth (no kyanite inclusions in garnet and garnet early porphyroblast remnants might be enveloped by kyanite crystals), the relations are not so clear in sample CS16_292b because most of the crystals are isolated from each other. Our observations did not provide convincing evidence that kyanite postdates garnet growth and both minerals are thus considered as cogenetic. Fe-rich (up to 1.23 wt. % FeO) syn-metamorphic kyanite crystals are oriented parallel to the main foliation S n . Chlorite analyses revealed a composition between clinochlore and amesite ( Figure A1; Table S11). All opaque minerals analyzed by EPMA were Ti-rich hematite (Table S6a). Rare ilmenite was identified by optical microscopy.

Mineralogy and Mineral Chemistry of Green-Colored Horizons
Sample CS16_297 is representative of the many green-colored horizons parallel to the foliation of the main schist/paragneiss unit of the Trikorfo area (Figure 3b). Element analyses are reported in Table S1a. The structure is mainly granoblastic with the Sn foliation locally underlined by phlogopite alignment. Identified minerals include biotite, quartz, plagioclase, garnet, epidote, Ti-hematite, titanite, rutile and apatite (Figure 6b). Magnetite with hematite exsolution was noticed. The mineralogy is dominated by quartz with layers rich in secondary spessartine and/or epidote/clinozoisite. Biotite has the composition of phlogopite with X Mg in the range 0.72-0.74 (Table S2b). Plagioclase is oligoclase with a composition of An 15 to An 21 , slightly depleted in CaO compared to CS16_292b (Table S3b). The content of TiO 2 in Ti-hematite does not exceed 8.15 wt. % ( Table S6b). The composition of titanite and rutile is given in Table S9a. Chlorite, staurolite, kyanite and sillimanite are not found in this rock.
The shape of garnet crystals is usually irregular. X Fe and X Mg are much lower (up to 21.69 and 3.24, respectively), and X Mn and X Ca (up to 58.00 and 37.29, respectively) significantly higher than garnet from sample CS16_292b. No significant and systematic chemical variation was noticed from border to border among all analyzed crystals (Table S4b) and garnet composition falls in the spessartine field ( Figure 7a). Epidote is abundant and shows a composition depleted in Mn 3+ ranging from epidote end-member to Fe 3+ -rich clinozoisite ( Figure 7b; Table S7a,j) and is responsible for the green coloration within these rocks.

Mineralogy and Mineral Chemistry of Calc-Silicate-Dominated Horizons
Moving towards the top of the formation, the thickness of the green layers increases and green-red horizons up to 80 cm thick are visible (Figure 3c,d). These beds are intercalated within the main foliation of the kyanite-bearing paragneisses/mica schists at proximity of carbonate/dolomitic rocks. They typically show calc-silicate assemblages more or less enriched in Mg-and Mn-bearing minerals with typical mineralogical zonation showing tremolitic hornblende, diopside, quartz, anorthite/bytownite, clinozoisite/epidote, grossular, titanite and pink clinozoisite from an outer band adjacent to host mica schists towards the center of the layer (Figure 3c,d) [20,24].
The dark green external part of the calc-silicate is dominated by Mn-bearing varieties of diopside (up to 1.44 wt. % MnO) and hornblende (up to 2.75 wt. % MnO) ( Table S8a,c,d). The progressive change in color from pale green to red (samples KT01-02, Th01-03) is mainly due to variation in epidote composition, whose EPMA analyses are presented in Table S7b,e-g, and reported in the Mn-Al-Fe ternary plot of Figure 7b.
Green epidote is reported in green layers parallel to rock foliation as in sample CS16_297 previously described and at the outer margins of the calc-silicate layers along with hornblende, anorthite/bytownite (Table S3d) and titanite (Table S9b) (Figure 6c,d). The composition evolves from epidote to clinozoisite end-members (Figure 7b). It may contain up to 1.04 wt. % MnO. LREE were not analyzed for epidote mineral group; their presence may mark the low total of some of the analyses.
Mn-poor clinozoisite exhibits subhedral to euhedral crystals. The pale pink crystals are usually found in the intermediate zone of the calc-silicate layers mostly associated with quartz, anorthite, titanite, diopside and hornblende. Available microprobe data revealed Fe 2 O 3 content between 1.36 and 6.11 wt. % (sample KT02) and between 1.54 and 11.06 wt. % (samples Th01 and Th03), and minor Mn 2 O 3 content below 0.67 wt. % (Table S7b,e,g).
Near gem-quality Mn-poor clinozoisite, with the same optical and chemical properties as clinothulite variety described by Bocchio et al. [64], is locally observed in the inner part of the calc-silicate layers, intergrown with Mn-grossular, titanite and quartz, both in the matrix of the rock and in fissures cross-cutting the foliation. The crystals occur in deep pink to red colors and form translucent subhedral to euhedral crystals up to 10 cm (Figures 3e and 6e; Table S7b,f). EPMA data revealed low contents of Fe 2 O 3 (<3.15 wt. % in sample KT02 and <2.90 wt. %, in sample Th02) and Mn 2 O 3 (<0.72 wt. % in sample KT02 and <0.21 wt. % in sample Th02), however high enough to be responsible for the pink-red coloration. (a) Garnet composition considering spessartine, grossular and almandine + pyrope end-members. Garnet composition is X 3 Al 2 Si 3 O 12 with the X site occupied by divalent cations, Mn 2+ for spessartine, Ca 2+ for grossular, Fe 2+ for almandine and Mg 2+ for pyrope, e.g., [65]. (b) Epidote group minerals, Ca 2 M 3+ Al 2 Si 3 O 12 (OH), classified according to one M 3+ cation substitution with Al in octahedral coordination with M 3+ = Al 3+ , Mn 3+ or Fe 3+ , corresponding to the clinozoisite, piemontite and epidote end-members, respectively [66,67]. Composition analyses for deep pink to red clinozoisite (clinothulite) are shown. (c) Andalusite-like minerals, Al 3+ (Al 3+ Mn 3+ Fe 3+ )SiO 5 , classified according to octahedral Al substitution with Mn 3+ or Fe 3+ [2,68].
In the calc-silicate lenses, pure grossular to Mn-bearing grossular form euhedral yellowish crystals of near gem-quality in close association with pink to red colored Mn-bearing clinozoisite and quartz (Figures 3e and 6e). MnO content ranges from 0.08 to 3.60 wt. % (Table S4c,e). All data are plotted in the ternary Mn-Ca-Fe + Mg diagram (Figure 7a).

Mineralogy and Mineral Chemistry of Mn-rich Schist Layers and Quartz Lenses
Near the top of the Trikorfo area, the formation is dominated by schist layers (samples Th04, Th05 and KT07) and intercalations of quartz lenses (sample KT08) particularly enriched in Mn-bearing minerals as shown by the assemblage of piemontite and Mn-rich epidote together with Mn-rich andalusite and spessartine (Figures 3f, 4d and 5f). Although EPMA analyses were not entirely satisfying (MnO between 27.99 and 45.31 wt. % and maximum total of all analyzed elements of 51.80 wt. %), probably reflecting the heterogeneity and variable Mn oxidation state of the minerals, Mn-oxides could be identified (Table S10). In both lithologies, piemontite and Mn-rich epidote form dark pink/purple-colored prismatic crystals with typical piemontite pleochroism. Microprobe analyses revealed piemontite in the andalusite-bearing mica schists with Mn 2 O 3 content up to 12.69 wt. % and Mn-rich epidote in association within the spessartine-braunite-quartz aggregates with Mn 2 O 3 in the range of 0.13 to 6.59 wt. % (Figure 7b, Table S7c,h,i). Phlogopite has X Mg of 0.96-0.97 (Table S2d). Plagioclase is oligoclase/andesine with a composition between An 22 and An 34 (Table S3c).
Accessory phases include monazite (relatively abundant in discrete layers), titanite and zircon. Tourmaline is very abundant in quartz lenses within mica schists, where it occurs as black crystals (schorl) up to 10 cm and associated with muscovite and chlorite (clinochlore). However not analyzed, olive green-colored tourmaline (crystals up to 3 cm) is observed. Tourmaline is also present in fissures intergrown with clinochlore, titanite and quartz. Finally, up to 10 cm tremolite white crystals were documented in sheared marbles [20].

Different Colors of Kyanite Crystals
Among the many spectacular minerals observed at the Trikorfo area, kyanite, found at the contact between quartzo-feldspathic coarse-grained veins and the paragneiss/mica schist rocks, towards the top of the series, is of particular interest. The crystals show a general orientation parallel to the regional foliation but may also be unoriented (Figure 4a). The up to 20 cm crystals are slightly bent as a result of weak plastic deformation. Over a small area of a few km 2 , besides the Fe-rich syn-metamorphic kyanite already described in the metamorphic paragenesis (Table S12a), we report at least four dominant types of kyanite based on their color: (i) dark blue-black (Figure 4a, Table S12b), (ii) zoned crystals with dark blue cores and transparent-green rims (Figure 4b,c), (iii) zoned crystals with dark blue cores and yellow rims (Figure 4d) and (iv) orange crystals generally associated with green Mn-rich andalusite (Figure 4e; Table S12c). Most of the crystals are inappropriate for faceting as marketable gemstones due to common fractures, but are at least proper for cabochon-shaped material. In addition, with regard to their overall aesthetic aspect, mainly large size and vivid color, these kyanite crystals are referred to as being of near gem-quality.

Bulk Element Concentration of Individual Kyanite Crystals
When observed down their <100> faces, most kyanite crystals show blue-green-colorless bands parallel to the c-axis. Representative kyanite samples have been selected for bulk element concentration analyses (zoned kyanite THA16_1, _2, _3, yellow THA08, pale orange THA12 and deep orange THA13; Figure 8). The elements are classified as oxides and trace elements according to their content in the host rock sample CS16_292b (* symbol in Figures 9a, 10 and 11).     (Table S1b), which appear to be correlated but randomly distributed between the samples. This is likely due to alteration phases in microfissures. The total amount of Fe 2 O 3 , MgO, MnO and TiO 2 increases with coloration from blue-green (1.10-1.21 wt. %), to yellow (1.51 wt. %) and orange (1.42-1.61 wt. %). Yellow and orange samples show markedly higher amounts of MnO (up to one order of magnitude) and MgO (Figure 9b).
Minor elements from the bulk rock analyses are found with amounts lower than 100 ppm in kyanite samples. Ba, Zn and Ga are the most abundant elements. Blue-green kyanite is significantly enriched in V, Nb, Sc and Ta. Yellow and orange samples of kyanite contain predominant amounts of As, Cr, Ni, Cu, W, Co, Sn, Th, Cs and Y. For all other trace element contents, there is no significant difference among the analyzed samples ( Figure 10; Table S1c).
The rare earth element concentration measured in each sample was normalized to chondritic abundance ( Figure 11). All kyanite samples have comparable chondrite-normalized REE patterns as CS16_292b host-rock with a slight decrease until Gd and then relatively flat towards heavy REE. As for all other elements, REE are less abundant in blue-green samples than in yellow and orange kyanite. The yellow (THA08) sample is characterized by a remarkable enrichment in light REE compared to other samples. This is to be correlated with the higher content of this sample with the elements attributed to alteration phases in microfissures such as K, Na, Li, Sr and Ba.

EPMA X-ray Maps and Element Distribution
For all the samples, Al and Si variations are difficult to interpret. Both elements appear relatively homogeneous throughout the crystals and their maximum/minimum concentration seems to be slightly anti-correlated, likely due to a problem of surface homogeneity. A decrease of Al is noticed on the edges of THA16 marking alteration or re-equilibration with minerals in the surrounding rock. The element Mn also does not show any obvious variation in the different analyzed areas. The only point of interest regarding Mn is that, accordingly with bulk rock analyses, the concentration is lower in blue THA16 (<0.10 wt. % MnO) than in the yellow THA08 (<0.18 wt. % MnO) and orange THA13 (<0.16 wt. % MnO) kyanite samples.
Fe concentration however shows very contrasting behavior. In kyanite samples THA16 (<1.27 wt. % FeO) and THA08 (<1.59 wt. % FeO), the variations are straight and parallel to the c-axis. At this scale, the higher concentration seems to be correlated with the clearer parts of the kyanite crystals. In sample THA13 (<1.56 wt. % FeO), Fe distribution seems to be oblique regarding the orientation of the c-axis and could be related to plastic deformation of the analyzed crystal ( Figure 12). Ti distribution is discontinuous and does not seem to be related with crystallographic orientation. Three marked areas of higher concentration are highlighted in THA16 (<0.05 wt. % TiO 2 ) and seem to correspond to the blue area of the crystal, at least for the left one. It is also the case to the right of the crystal, but the blue coloration is barely visible. The direct correspondence between blue coloration and Ti concentration is however difficult to establish firmly as most of the blue area is below the sample surface. In sample THA08 (<0.02 wt. % TiO 2 ) also, the irregular distribution of Ti reveals enrichment at the central right part of the crystal where a thin blue zonation is observed. On THA13 (<0.06 wt. % TiO 2 ), Ti distribution follows that of Fe, though anti-correlated.

Orange Kyanite and Green Andalusite Assemblages
This atypical polymorphic association is observed towards the top of the formation (i) in kyanite-phlogopite-piemontite-quartz-plagioclase ± muscovite ± spessartine ± braunite mica schists, and (ii) in quartz-feldspar-muscovite lenses. Mn-rich andalusite forms dark green-colored subhedral to euhedral crystals up to 7 cm. In both cases, kyanite and andalusite seem to have grown at equilibrium (Figure 6f). Green andalusite is commonly zoned likely due to twinning, and exhibits deformation lamellae [29]). EPMA data reveal elevated Mn 2 O 3 contents, up to 15.98 wt. % (Figure 7c; Table S13a,b).

LA-ICPMS Profiles
The limits of detection of LA-ICMPS are significantly lower than those of EPMA and many elements can be analyzed sequentially during the same analytical run. Edge to edge sections were conducted at the same position, or approximately, as that for the EPMA profiles on the three kyanite samples (Figure 13). Al and Si variations serve as internal standards to provide quantitative elemental evolution through kyanite and are not reported on the figures. All analyses of Ca, Cu, Ag and Au were below the limits of detection. For each sample, the elements are reported according to three different groups with (i) the dominating elements (Fe, Mn, Ti, Mg and Zn) on top, (ii) the elements generally enriched in alteration minerals (K, Na, Li, Sr and Ba) in the middle, and (iii) the other elements, mainly transition metals, only found as traces in the host-rock (Cr, Pb and V) on the lower part.
Kyanite sample THA08 shows large variations in the concentration of the elements K, Na, Li, Sr and Ba. Moreover, it can be noticed that these elements are mainly concentrated either on the sides of the crystals or along fractures inside the crystals, which are common in THA08. Accordingly, these elements are interpreted as the result of alteration or late equilibration with host-rock minerals (feldspars and micas like muscovite and phlogopite).
In all samples, Fe is by far the most common element besides Al and Si. Local concentration decreases are observed where the amount of alteration elements increases, as is clearly visible on sample THA08. Mn and Ti are in the same order of magnitude in THA16 and THA08 but with Ti slightly higher in the blue-green kyanite. Mn concentration does not correlate with any color change, whereas that of Ti shows a relatively good correspondence with the blue color bands, even if they are below the sample surface. This is especially well expressed in THA08 where Ti > Mn at the intersection with the blue band. It is also visible on THA16, where the three Ti highs seen in the EPMA mapping are distinguishable, and correspond to a blue color on the left. However, it must be kept in mind that blue color is mostly below the crystal surfaces. Also, it must be noted that the laser ablation itself down-cuts through the sample (ablation spot of 32 µm) and hence what was actually analyzed may largely be below the surface. In sample THA13, Mn dominates largely over Ti and is one order of magnitude higher as in samples THA16 and THA08. The Mg concentration variation seems to partly correlate with edges and fractures and can be interpreted at least in part due to equilibration with host-rock phlogopite. Along the profiles, Zn is equally present and its concentration evolution is marked by rapid variations up to one order of magnitude. In some parts of the kyanite crystals, Zn is almost as abundant as Mn, Ti and Mg. These variations are not related to any characteristic feature of the crystals. The other elements do not show any obvious correlations. V is found above the limit of detection (1 ppm) only in sample THA16. No systematic evolution of Cr could be observed, because of concentration too close to the limit of detection (15 ppm). Pb shows the same evolution as the elements attributed to alteration and post-crystallization re-equilibration, with higher concentrations at the edges of the crystals and within the fractures.

UV-Near-Infrared Spectroscopy
All spectra (blue, yellow and orange kyanite samples) show an absorption continuum rising from infrared or mid-visible to UV, making absorption range from 315 nm to below (Figure 14). They also all show sharp and weak absorption peaks at 370, 382, 432, 446 nm. Spectra of blue samples are dominated by large absorption bands at 520 and 620 nm that, combined with the previous spectral features, generate a transmission window around 480 nm, in the blue sample. Spectra of orange and yellow kyanite samples are similar. In addition to the common spectral features described initially, both show a band at about 470 nm and a weaker band at 890 nm. These two bands are weaker in the yellow than in the orange kyanite sample. As a consequence, transmission increases from the green (550 nm) toward the infrared, hence the orange and yellow colors of the samples. The allocation of these spectral features to specific ions or defects is proposed in the Discussion section.

Phase Diagram Calculation of the Paragneiss
The PT pseudo-section is presented in Figure 15 for the PT range of ca. 550-750 • C and ca.  Figure A2). In consequence the phase diagram section was computed with this XFe 2 O 3 molar ratio. Considering the P-T range used for this pseudosection (550-750 • C; 5.5-10.5 kbars), hematite-rich ilmenite with ca. 50-70% of hematite end-member is stabilized in each stability field of ilmenite and, thus, labelled as "Ti-Hem" in Figure 15. For the specific PT field of the described paragenesis (ca. 670 • C and 9.1 kbars), the calculated ilmenite in the ilmenite-hematite solid solution is around 33%, consistent with the 10-35 % measured by EPMA (Table S6a).  Chlorite thermometry [70] made on retrogressed phlogopite of the paragneiss CS16_292b and schist sample KT03 indicates temperatures in the range 398-334 • C (Table S2c for phlogopite and  Table S11a,b for chlorite composition). The values given by the clinochlore analyzed within sample piemontite-andalusite-muscovite-kyanite schist KT07 are in the range 280-326 • C (Table S11c), while the values from Mn-bearing clinochlore from the late brittle vein MB13 are from 233 to 315 • C (Table S11d).

Discussion
The Trikorfo area on Thassos Island represents a unique mineralogical locality where near gem-quality Mn-rich minerals are described. A detailed structural/mineralogical study permits to describe metamorphic parageneses and local metasomatic assemblages. The dominant formation is made of paragneisses/mica schists containing a metamorphic paragenesis, including kyanite. Towards the top of the formation, different units are described parallel to the metamorphic foliation, with some of them bearing large and colored minerals: (i) calc-silicate dominated horizons, (ii) Mn-rich schist layers and quartz lenses, (iii) kyanite-bearing quartz-feldspar coarse-grained veins, and (iv) late brittle veins cross-cutting the metamorphic foliation. The discussion is divided into three parts where (i) the near gem-quality minerals are described and the elemental substitutions leading to their color variation are discussed, (ii) the distinction between metamorphic parageneses and local metasomatic assemblages is made, (iii) a PT-deformation exhumation path for the conditions of formation of the different units is proposed. Kyanite is found within the metamorphic paragenesis and as near gem-quality in quartz-feldspar coarse-grained veins; it is therefore the guidance mineral of the study and the discussion.

Colored Minerals and Cation Substitution in the Trikorfo Area, Thassos Island
Dispersed metal ions in silicate and calc-silicate crystal structures lead to color variation of particular interest for mineral collectors. Silicate and calc-silicate minerals of the Trikorfo area are large and vividly colored, close to gem-quality. Moreover, thanks to trivalent cation substitutions (Mn 3+ and Fe 3+ ) with Al, the area hosts an outstanding variety of colors within the minerals. This chapter aims to depict the main substitutions observed within kyanite, andalusite, epidote and garnet minerals.
The common natural color of kyanite (Al 2 SiO 5 ) is blue but it can also be green, yellow, orange, white, black, grey or colorless as a function of the nature of elemental substitutions [4]. Moreover, several alternating colors can be found along bands parallel to the c-axis in a single kyanite crystal, marking the presence of trace elements in the crystal structure. Kyanite has two Al sites in six-coordination and one Si site in four-coordination. The occupying site and bonding state of transition elements determine the color and physical properties of the crystals [71]. In the Trikorfo area, besides syn-metamorphic crystals, kyanite (zoned blue to green/yellow-transparent, yellow and orange) is also observed in quartz-feldspar coarse-grained veins or in quartz lenses. Bulk chemical composition analyses reveal that each type of color kyanite has its own signature with a characteristic distribution of measured oxides and trace elements in the minerals (Figures 9 and 10). V and Cr distribution is markedly different between zoned blue (up to~10 ppm and~2 ppm, respectively) and orange kyanite (below detection limit and up tõ 15 ppm, respectively), suggesting a local control from the host-rock composition. However, the low concentrations point to a similar felsic origin [4] that should be attributed to interaction with the main paragneiss unit of the area. Orange kyanite localities are very rare on Earth. The Loliondo deposit (Tanzania) is characterized by elevated MnO content up to 0.23 wt. %. Polarized optical absorption spectra show that the orange coloration is governed by crystal field d-d transitions of Mn 3+ replacing Al in six-fold coordinated triclinic sites of the kyanite structure [9]. Even if orange kyanite from the Trikorfo area is richer in MnO than the yellow and blue zoned crystals, their MnO content (0.07 wt. % for the bulk sample analysis and locally up to 0.16 wt. % for EPMA analyses) remain far below that of orange kyanite from Tanzania. This is consistent with an orange color less intense in kyanite from Trikorfo than from Loliondo. All kyanite samples show Zn in a significant amount but without any evident correlation with crystal structure or color. The total Cr content is higher in the yellow/orange samples than in the blue kyanite ( Figure 10). However, LA-ICPMS sections do not show conclusive information regarding the distribution of Cr due to concentration close to the limit of detection of 15 ppm ( Figure 13). Moreover, kyanite colored by Cr 3+ shows absorption spectra dominated by two broad absorption bands centered at 405-420 nm and 595-625 nm [72,73]. No such absorptions are observed in any of our samples, making Cr 3+ an unrealistic chromophore in this case. This is consistent with Cr content very low in our samples (typically 1 to 20 ppm; Figure 10; Table S1c) compared to blue kyanite colored by dispersed Cr 3+ (typically Cr 2 O 3 > 1 wt. %; [73]). Hereafter, only Mn 3+ , Fe 3+ and Fe 2+ -Ti 4+ substitutions with Al 3+ are thus discussed. The kyanite crystals analyzed by UV-near infrared spectroscopy show weak absorption bands at 370, 381, 432 and 446 nm due to isolated Fe 3+ [7]. The large absorption bands at 520 and 620 nm, responsible for the blue color of THA16, are due to Fe 2+ -Ti 4+ inter-valence charge transfer (IVCT; [74][75][76]). No spectral features due to the Fe 2+ -Fe 3+ IVCT were found. In the spectra of yellow (THA08) and orange (THA13) kyanite samples, the band at 470 nm, that importantly contributes to color, is consistent with the main band attributed to isolated Mn 3+ [9]. Hence, the origin of the yellow to orange color in our samples is attributed to isolated Mn 3+ , with some minor contribution of isolated Fe 3+ . The contribution of the weak 890 nm band remains unattributed.
EPMA X-ray maps and LA-ICPMS profiles show some relation between the blue zones and the Ti content of the zoned kyanite crystals THA16 and THA08 (Figure 12a,b and Figure 13a,b). Interestingly, orange kyanite (THA13) contains more TiO 2 than blue, zoned crystals. This shows that the presence of Ti alone is not sufficient to generate a blue color: Ti 4+ atoms need to be crystallographically linked to Fe 2+ atoms for the charge-transfer phenomenon to occur, and hence the blue color to appear. If for any reason, the charge transfer does not occur (titanium atoms may be isolated, or iron ions may be trivalent), the blue color due to IVCT does not appear. A similar behavior is documented in other colored minerals such as corundum: some yellow sapphires sometimes contain more titanium than blue sapphires colored by Fe 2+ -Ti 4+ IVCT, e.g., [77]. Again, in this case, titanium concentration alone is not correlated with the blue color.
Orange kyanite is associated with green andalusite in quartz lenses (Figure 3f). Mn 3+ might substitute for Al 3+ in andalusite mineral formula, resulting in deep green, Mn-rich andalusite (Al, Mn) 2 SiO 5 [10,78]. This andalusite variety is named viridine in reason of its brilliant green color [12]. Viridine occurrences (Germany, Belgium, Japan, Australia, Brazil and Tanzania) usually provide small or altered crystals [16,78]. Large viridine crystals of the Trikorfo area reach up to ten centimeters large, are vivid-to-deep green and unaltered. Their occurrence is promising for exploration of gem-quality, marketable viridine in the area. Compared to kyanite and sillimanite, andalusite incorporates maximum contents of Mn 3+ due to its crystal structure [2]. Experimental works quantified that the Mn 2 SiO 5 end-member cannot be higher than 6 mol% for kyanite and 1.7 mol % for sillimanite, but can reach 22 mol % for andalusite [2,79]. In the present study, Mn 2 O 3 can reach up to 15.98 wt. %, yielding 17.5 mol % Mn 2 Si 2 O 5 , close to the maximum possibly accepted by andalusite.
In the Trikorfo area, epidote group minerals are represented by green epidote, pink to red Mn-rich epidote and Mn-rich clinozoisite (variety clinothulite), purple piemontite, and allanite, which is only reported at the thin-section scale by Voudouris et al. [20]. Epidote, clinozoisite and piemontite represent the Fe-, Al-and Mn-rich end-members respectively of the epidote mineral family. Solid solution involves substitution of trivalent cations (Al 3+ , Fe 3+ or Mn 3+ ) in three non-equivalent octahedral sites in the basic formula for epidote Ca 2 M 3+ Al 2 Si 3 O 12 (OH) [66,67]. True piemontite (up to 12.69 wt. % Mn 2 O 3 ) is restricted to Mn-rich andalusite bearing samples, indicating that these mica schist layers are locally particularly enriched in manganese. Mn-bearing epidote (up to 7.76 wt. % Mn 2 O 3 ) is usually associated with spessartine-Mn oxides-quartz assemblages. The Mn-poor clinozoisite/epidote and Mn-poor clinozoisite (clinothulite) from the Trikorfo area contain up to 0.72 wt. % Mn 2 O 3 which is responsible for their pink to red color. In accordance to Bonazzi and Menchetti [80] and Franz and Liebscher [67] the coloration in pink epidote/clinozoisite is due to Mn 3+ and not Mn 2+ in the mineral structure.
Four types of garnet are recognizable on the basis of their colors: (i) syn-metamorphic millimetric red almandine (± pyrope)-spessartine, (ii) plurimillimetric Ca-Fe-rich yellow to orange spessartine within the green horizons, (iii) centimetric almost pure pale yellowish brown grossular associated with epidote and clinozoisite in calc-silicate aggregates, and (iv) centimetric deep amber orange almost pure spessartine associated with Mn-rich andalusite, muscovite and orange kyanite. The colors of these garnet crystals directly reflect their main dominant cation (Ca 2+ , Fe 2+ or Mn 2+ ). Mn-rich garnet is found within the metamorphic assemblage and the composition to almost pure spessartine is only very localized to certain mica schist layers in association with orange kyanite and green andalusite, and to brittle veins cross-cutting Mn-rich horizons close to the top of the Trikorfo area.

Metamorphic Parageneses versus Metasomatic Assemblages
The Trikorfo area is dominated by garnet-kyanite-biotite-hematite-plagioclase ± staurolite ± sillimanite paragneisses/mica schists (sample CS16_292b), with intercalations of metacarbonate units near the base and top of the formation and local calc-silicate layers and Mn-rich mica schists, also close to the top of the formation. Field and mineralogical observations suggest that the epidote-grossular-bearing calc-silicate layers represent local inhomogeneities in the bedding of the sedimentary protoliths and/or compositional changes that have been produced by localized metasomatic processes during regional metamorphism, as suggested for similar layers at Therapio, Evros area [81]. In this sense, the mineralogy of the calc-silicate layers could have developed during prograde metamorphism of a Mn-rich, calcareous pelitic protolith, followed by vein formation and local metasomatic reactions during retrograde metamorphism accompanying the exhumation of Thassos Island during the Oligocene-Miocene.
Alternatively, the skarn-similar mineralogy of the calc-silicate layers (e.g., grossular, diopside, hornblende and epidote) could have been formed by fluids released by granitoids during contact metamorphism with the carbonate (calcite and dolomite) rocks found at the top of the Trikorfo area. With the lack of radiometric data, the exact age (i.e. Hercynian?) of the orthogneiss at Trikorfo is not known but is likely older than alpine deformation. Although not exposed on Thassos Island, Oligo-Miocene magmatism is described nearby in the Rhodope area with the intrusion of syntectonic plutons 30 km to the north of Thassos Island, e.g., [52][53][54]. Melfos and Voudouris [27] considered that widespread gold mineralization in the area is related to this Miocene magmatic activity. Tourmaline at Trikorfo also occurs as a retrograde mineral, resulting in tourmalinite and quartz-tourmaline veins cross-cutting the metamorphic foliation. According to van Hinsberg et al. [82] metasomatic introduction of boron on the retrograde path is most commonly associated with the intrusion of late granites in orogenic belts.
The Trikorfo area shows the superposition of late secondary minerals over the metamorphic paragenesis. These silicates and calc-silcates, including a large variety of garnet and epidote composition, kyanite, andalusite, diopside, hornblende, Ti-hematite, magnetite and tourmaline indicate post-metamorphic metasomatic reactions where hot, chemically active fluids that altered the metamorphic assemblage. Furthermore Mn-rich mineralogical assemblages are found. This type of assemblage is relatively rare on Earth because it requires highly oxidized environments [15][16][17][18].  [18], garnet which coexists with piemontite is uniformly more spessartine-rich.
Garnet and epidote are minerals showing large Fe-Al-Mn solid solutions sensitive to PT and f O 2 environmental conditions. In the case of epidote composition, oxygen fugacity is more important than temperature [18]. Due to extensive solid solution in terms of Mn/Fe and (Mn + Fe)/Al, Mn-rich garnet and epidote may coexist with element ratio of the different phases sensitive to temperature and degree of oxidation. With decreasing f O 2 , piemontite becomes poorer in Mn and garnet and epidote minerals might tend to become more aluminous with increasing temperature [18].
Spessartine can be stable over a wide range of f O 2 and PT conditions from greenschist to amphibolite-granulite facies and is evidence of elevated amounts of Mn in the rock [15,83,84]. Natural occurrences of piemontite have shown evidence for metamorphism at high oxygen fugacity [18]. A highly oxidizing environment is shown by the presence of Ti-hematite, Fe 3+ -and Mn 3+ -rich epidote, Mn 3+ -rich andalusite, Fe 3+ -Mn 3+ -rich kyanite, and so on. The intercalation of decimeter-thick layers of epidote-spessartine-Mn-rich andalusite in less oxidized levels with epidote-spessartine(almandine) assemblages implies (i) significant f O 2 gradients at restricted scale, and (ii) internally controlled f O 2 resulting from oxidized protoliths containing minerals such as Mn-oxides capable of buffering f O 2 to high pressure levels during metamorphism [17,18]. Preservation of the f O 2 gradients supports evidence for a reduced mobility of oxygen throughout high-pressure metamorphism [17] or the localized circulation of highly oxidized fluids during metasomatism. The strong variability of f O 2 dependent mineral assemblages, at all scales, within the Trikorfo area thus indicates a strong local buffering of f O 2 during post-metamorphic metasomatic reactions. As for f O 2 , the mineralogical assemblages also suggest a local control in Mn concentration. Indeed, Mn-rich layers, which are mostly found at the top of the series, are scarce and never exceed a few centimeters thick.
Orange kyanite was found coexisting together with Mn-rich andalusite in quartz-feldspar coarse-grained veins and as post-deformation assemblages in mica schists towards the top of the unit in the Trikorfo area. Mn-rich andalusite may also have formed during decompression in veins within host rocks containing kyanite, e.g., [85]. However, petrographic observations suggest growth at equilibrium between kyanite and andalusite. Co-existing Al 2 SiO 5 polymorphs occur in various metamorphic rocks and may form either during regional metamorphism or a combination of regional and contact metamorphism [86,87]. When kyanite, sillimanite and andalusite are free of chemical impurities, their equilibrium coexistence is invariant. However, these three polymorphs show extensive solid solution and can incorporate significant amounts of Fe 3+ and Mn 3+ , substituting for Al, in highly oxidized rocks [16,88,89]. This results in an increase of the variance of the system and the possible coexistence of kyanite, andalusite and sillimanite in mutual equilibrium over a measurable range of PT conditions [90]. The presence of Mn-bearing andalusite, and especially the association with braunite, quartz and spessartine, indicates high f O 2 [15].
Structural and mineralogical observations thus point to a two-stage evolution of the Trikorfo area, where metamorphic units were affected locally by metasomatic events during exhumation. Fluid circulation leads to fluid-rock interactions and crystallization of large, vividly-colored minerals. The chemistry of the metasomatic mineralogical assemblages is controlled by the chemistry of the protolith more or less enriched in manganese and by local f O 2 buffer.

PT-Deformation Conditions of Metamorphic Equilibrium and Metasomatic Reactions
The Rhodope domain exposes a large variety of HP-rocks of the Variscan continental crust and Mesozoic sediments that were subjected to Alpine subduction (Early Cretaceous-Eocene) and subsequent exhumation (Early Oligocene-Miocene) [91]. The recorded conditions can be as high as >19 kbars and 750-800 • C with local partially amphibolitized kyanite eclogites, e.g., [37,91].
Thassos Island, belonging to the South Rhodope metamorphic core complex, displays Oligocene to Miocene metamorphic and exhumation history [36,40]. In response to crustal thickening during Alpine tectonics, the geothermobarometric analyses of amphibole gneisses suggest metamorphic conditions of approximately 620 • C/4.7 kbars for the lower unit and 580 • C/2.4 kbars for the intermediate unit [47,48]. On the basis of the observation of a garnet-kyanite-sillimanite assemblage, Dimitriadis [24] estimated the conditions achieved by the mica schists of the Trikorfo area around 5.5 ± 1.5 kbars and 600 ± 50 • C. Brun and Sokoutis [40] reported an unreferenced personal communication of Kostopoulos (2004) stating that the garnet amphibolites yielded temperatures ranging from 545 to 660 • C, as described in Dimitriadis [24], but at significantly higher pressures from 6.5 to 9.5 kbars. On the basis of the garnet-staurolite-kyanite-plagioclase-quartz-phlogopite-Ti-hematite assemblage, our thermodynamic modeling indicates metamorphic conditions at Trikorfo of ca. 9.1 ± 1.3 kbars and 670 ± 40 • C, at significantly higher P and T conditions than previously estimated. The field given in Figure 15 is the only one where kyanite and staurolite are stable together. Instead of staurolite, Dimitriadis [24] reported the presence of sillimanite growing at the expense of kyanite within the metamorphic gneiss. The conditions of the Trikorfo formation should rather stand in the high-temperature/low-pressure part of the domain indicated by our stable metamorphic assemblage, with sillimanite post-dating kyanite. The MnO content of the calculated equilibrated garnet composition is a bit higher than that measured in sample CS16_292b. It is probably due to the fact that the solid solution model used in our thermodynamic calculation is not fully appropriate. The addition of MnO in the system has a small effect on the garnet mode but a great effect on the garnet stability field. Indeed, in agreement with White et al. [92], the incorporation of MnO in our phase diagram calculations (in the MnNCKFMASHT(O) system) leads to a shift to lower PT conditions of the garnet-in line, compared to a system without MnO. The PT conditions obtained for the peak metamorphic assemblage recorded by the kyanite-bearing paragneisses are in good agreement with the geological context and therefore we are confident with the calculated phase diagram, even if the modeled garnet composition is not strictly the same than the measured one. Indeed, this MP-HT metamorphism recorded by paragneisses from the intermediate unit is comparable to those estimated at the Thassos scale. This amphibolite metamorphism, which is associated with the continuous development of the extensive low-angle S n foliation since ductile to brittle deformation, could mark the beginning of the nappe stack exhumation.
Kyanite and kyanite-andalusite assemblages are found at the interface between metamorphic host rock and quartzo-feldspathic veins. The crystals are generally weakly deformed and aligned parallel to regional lineation. This indicates a local reaction where a SiO 2 -rich fluid interacts with the Al 2 O 3 -rich host-rock [93], as indicated by their Cr/V ratio [4], to form either kyanite or kyanite-andalusite as a function of the PT conditions and the variance of the system linked to available MnO [90]. This reaction occurs at ductile conditions while the rock still records stretching lineation. The temperatures derived from the Ti-in kyanite thermometer from Müller et al. [4] are in the range 670-790 • C coherent with the Perple_X modeling. This would then indicate that metasomatic kyanite form at or close to the metamorphic peak conditions. However, EPMA and LA-ICPMS data (Figures 12 and 13) show that Ti distribution is variable within the crystals and these temperatures should be considered with caution. Except for some of the kyanite and andalusite crystals, the other metasomatic silicate and calc-silicate minerals of Thassos are usually not oriented. Dimitriadis [24] suggested retrogression under static conditions as evidenced by undeformed matted fibrolite and unoriented growth of chlorite. However, continuous deformation during the main kinematic event is evidenced by coherent stretching and mineral lineation and brittle vein opening to temperatures as low as 230 • C, as indicated by clinochlore growth in late brittle veins. The presence of non-oriented minerals should be the result of crystal growth from hydrothermal fluids during metasomatic reactions. The presence of calc-silicate minerals (grossular, diopside, hornblende, epidote), especially located to the top of the series where the paragneisses/mica schists are in contact with marble units, implies decarbonation during metasomatic reactions.

Conclusions
The garnet-kyanite-biotite-hematite-plagioclase ± staurolite ± sillimanite paragneiss of the Trikorfo area (Thassos Island, Rhodope Massif, northern Greece) recorded metamorphic conditions of ca. 670 ± 40 • C and ca. 9.1 ± 1.3 kbars. Structural, petrographic and mineralogical observations indicate that localized metasomatic reactions occurred during the exhumation of the HP unit locally containing layers enriched in Mn with a strong local f O 2 buffering. Metasomatism is likely the result of Miocene intrusion of granitoids, which are related to ore deposits on Thassos Island. Metasomatic reactions first occurred under ductile conditions in an extensive context close to the kyanite-andalusite stability curve and continued until purely brittle conditions as indicated by the presence of late veins cross-cutting the metamorphic foliation. The result is a mineralogical site that could be regarded as a unique geotope with uncommon (Mn-rich) silicate and calc-silicate mineral assemblages similar to what is observed in skarn deposits where hot active hydrothermal fluids interact with carbonate rocks.
Crystals of kyanite, green andalusite, garnet (grossular and spessartine) and red epidote-clinozoisite are large (up to several centimeters), show vivid colors, and are suitable for cabochon-shaped gemstones. As such, the Trikorfo locality can be regarded as a promising area for the exploration of true, facetable gemstones. Their genesis due to metasomatic reactions also underlines the important role of metasomatism for gemstone formation in general, as previously noted in the literature [94][95][96].
Our combined petrographic, geochemical and spectroscopic observations permit us to demonstrate that the orange color of kyanite in the Trikorfo area is mainly due to Mn 3+ substitution. Although the Mn content is significantly lower than that in the Loliondo (Tanzania) deposit, Thassos Island (Rhodope, Greece) can now be classified as the second locality worldwide where Mn-rich orange kyanite is reported.