Mineralogy and Geochemistry of Ultramafic Rocks from Rachoni Magnesite Mine, Gerakini (Chalkidiki, Northern Greece)

: The importance of magnesite for the EU economy and industry is very high, making the understanding of their genesis for the exploration for new deposits a priority for the raw materials scientific community. In this direction, the study of the magnesite-hosting ultramafic rocks can be proved very useful. For the present study, ultramafic rock samples were collected from the magnesite ore-hosting ophiolite of the Gerakini mining area (Chalkidiki, Greece) to investigate the consecutive alteration events of the rocks which led to the metallogenesis of the significant magnesite ores of the area. All samples were subjected to a series of analytical methods for the determination of their mineralogical and geochemical characteristics: optical microscopy, XRD, SEM, EMPA, ICP–MS/OES and CIPW normalization. The results of these analyses revealed that the ultramafic rocks of the area have not only all been subjected to serpentinization, but these rocks have also undergone carbonation, silification and clay alteration. The latter events are attributed to the circulation of CO 2 -rich fluids responsible for the formation of the magnesite ores and locally, the further alteration of the serpentinites into listvenites. The current mineralogy of these rocks was found to be linked to one or more alteration event that took place, thus a significant contribution to the metallo- and petrogenetic history of the Gerakini ophiolite has been made. Furthermore, for the first time in literature, Fe inclusions in olivines from Greece were reported.


Introduction
Magnesite ore deposits are closely associated with ophiolitic ultramafic rocks-more specifically dunites and hartzburgites. Magnesite is the main source of magnesium for industry and the raw material needed for the production of various kinds of magnesia (MgO). The large number of applications of magnesia make magnesite one of the most important minerals for industry and the economy. In the European Union, magnesite is considered a critical material under the Raw Materials Initiative (RMI). Apart from magnesite and its industrial products, recently many research efforts have been in progress for the valorization of the by-products produced during magnesite mining and processing [1][2][3][4][5][6]. These mine wastes, namely serpentinized peridotites, are not considered hazardous for the environment or public health; nevertheless, their large volumes make their valorization a critical issue for the achievement of zero waste production and circular economy policies of the European Union. The importance of magnesite ores, their products and mining by-products, make the study of their hosting of ultramafic rocks extremely important for the understanding of their genesis and the exploration of new deposits.
Ultramafic rocks in ophiolite complexes are often subjected to alteration by the infiltration of aqueous fluids [7], resulting in the formation of serpentinites. The infiltration of CO2-rich fluids can produce magnesite through the interaction with a serpentinite or a peridotite. Until today, the source of CO2 remains debatable. Abu-Jaber and Kimberley [8] suggest that Mg is supplied by the ultramafic host, while the source of carbon remains uncertain. A possible source could be an influx of metamorphic-degassing CO2 and CH4 into groundwater which flows into the ultramafic rocks. CO2rich water would produce magnesite after interacting with a serpentinite or a peridotite. This model requires a source of carbon in depths greater than 10 km (T > 300 °C) for the liberation of volatiles and their mixing with a carbon-poor solution. This mechanism is also supported by Pohl [9], but without excluding the possibility of meteoric water and hydrothermal fluids mixing. Other researchers [10][11][12][13] considered the oxidation of organic matter in sedimentary formations, which is then infiltrated by meteoric water into the ultramafic rocks in contact as the main source of CO2-rich fluids. Finally, more recently, Kelemen and Matter [14] and Kelemen et al. [15] supported the exclusively meteoric origin of CO2 and excluded magmatic/metamorphic procedures and the organic matter oxidation as possible sources, with their claim mainly based on the absence of 14 C from magnesite ores. Magnesium is provided by the ultramafic host rock by the hydrolysis of Mg-rich minerals. Serpentine, in particular, was demonstrated to be soluble at low T and P in pure water [16], and to release Mg 2+ ions during weathering [17]. Moreover, an extensive alteration of ultramafic rocks by CO2-rich fluids, can result in "listvenitization", leading to the replacement of the ultramafic protolith by a quartz-carbonate rock denominated listvenite [18]. According to Tsirampides and Filippidis [19], the indicative magnesite ore deposits of Greece are almost 300,000 tons with their value being estimated at ~EUR 10 billion. Exploitable magnesite deposits in Greece are found mainly in two areas: -Chalkidiki (e.g., Gerakini, Ormilia, Vavdos, Vasilika); and -Evia island (e.g., Mantoudi, Limni, Troupi, Petisounas, Afrati, Papades) Other magnesite occurrences in Greece include the areas of Gomati and Nea Roda (Chalkidiki), Nigrita, Kozani, Grevena, Atalanti, Ermioni and Lesvos. In this work, new insights into the mineralogy and geochemistry of the magnesite-hosting ultramafic rocks of the Gerakini area are presented.

Geologic Setting
The magnesite deposits of the Chalkidiki (Northern Greece) area are hosted within ophiolitic massifs. According to Kaufmann et al. [20], the Chalkidiki ophiolites are placed into the Circum Rhodope Belt (CRB). The Circum Rhodope Belt sensu stricto comprises low-grade metamorphosed Triassic and Jurassic sedimentary rocks fringing the high-grade metamorphic rocks of the Serbo-Macedonian and Rhodope massifs in northern Greece. Main outcrops occur in the easternmost part of the Vardar suture zone in the Chalkidiki peninsula (Melissochori Formation; formerly Svoula flysch) and in Thrace (Makri unit and Melia Formation) [21].
The tectonostratigraphic relationship between the CRB and the high-grade metamorphics has not yet been clarified: initial research works claimed that the CRB represents the original Mesozoic stratigraphic cover of the Serbo-Macedonian crystalline basement, whereas more recently, the existence of two distinct greenschist facies metamorphic events is proposed [21].
The Circum Rhodope Belt is further divided into three geologic formations: Nteve Koran-Doubia Formation, Melissochori-Cholomontas Formation and the Aspri Vrisi-Chortiatis Formation. Ophiolitic massifs hosting magnesite deposits are found excessively in the Aspri Vrisi-Chortiatis Formation, i.e., the ophiolites of Vavdos and Gerakini-Ormilia areas. This ophiolitic complex is a part of the Eastern Ophiolite Belt of Greece.
The Gerakini-Ormilia ophiolite mainly consists of ultramafic and mafic rocks. The lower part of the ophiolite is dominated by serpentinized harzburgites and dunites which are locally intruded by pyroxenitic and gabbroic veins. These peridotites host the stockwork-type magnesite ores which are currently exploited. Additionally, locally some chromite bodies are found which were exploited in the past. Moreover, near faults, massifs of listvenites are preset, hosting vein-type magnesite ores. The upper part of the ophiolite is dominated by massive pyroxenites which are also intruded by magmatic veins.
For the present work, several rock samples of serpentinized ultramafic rocks (dunites and harzburgites) were taken from the Rachoni open pit site in Gerakini ( Figure 1).

Analytical Methods
Samples were collected for the present study, including sixteen samples (W1-W16) from the Rachoni magnesite open pit mine in Gerakini. Additionally, one listvenite sample (W17) was collected from the nearby Ugo magnesite open pit mine for comparison. Polished sections and powdered material were prepared for each sample. The mounted polished sections were examined under optical microscope.
The mineralogical characterization of these samples was performed by X-Ray Powder Diffraction (XRPD) using a Bruker D8 Advance diffractometer at the Chemical Engineering Department of the Aristotle University of Thessaloniki. Data were collected at 2θ from 5° to 70°, with a scan time of 0.2 s and increment of 0.02. The filament X-ray tube of Cu with a wavelength 1.5418 Å and the detector LYNXEYE (1D mode) were used. The identification of the existing mineral phases in XRPD patterns was performed using the PCPDF database while for the quantification, the Rietveld method was applied through the Fullprof software.
Scanning electron microscopy (SEM) images and microprobe analyses (EMPA) were conducted with a JEOL 8200 (JEOL Ltd., Akishima, Japan) equipped with a wavelength dispersive system (WDS) at the Earth Sciences Department of the University of Milan. The microprobe system operated using an accelerating voltage of 15 kV, a sample current on brass of 15 nA, a counting time of 20 s on the peaks, and 10 s on the background. The approximate detection limit was 0.01 wt. % for each element. A series of natural and synthetic standards were used for analysis: Si, Ti, Al, Cr, Fe, Ni, Mn, Mg, Ca, Na and K. EMPA data were corrected for matrix effects by applying the PRZ algorithm included in JEOL software.
The major and trace elements of the whole rock samples were analyzed using a Perkin Elmer ICP-OES (Medtech, Waltham, MA, USA) and a Perkin Elmer Sciex Elan 9000 ICP-MS (Medtech, Waltham, MA, USA), following the LiBO2/LiB4O7 fusion and HNO3 digestion of the fused solid sample at the external collaborating laboratories (Bureau Veritas Commodities Canada Ltd., Vancouver, BC, Canada). Based on the chemical analyses of the rock samples, their CIPW normative mineralogy was subsequently calculated.

Petrography and Mineralogy
Petrographic and mineralogical features are based on XRPD and microscopical observations. Ultramafic rocks comprised of 13 samples of serpentinized peridotites, and one serpentinite sample with plagioclase and a late clay alteration (W3), with three listvenite samples, two from Rachoni mine (W1 and W2) and one from Ugo mine (W17). The three lithologies, which are well distinguished based on their mineralogical assemblage, are reported in Table 1. Microscopic observations on the major mineral phases identified were verified by the XRPD results (Figures 2 and 3).   Peridotites are characterized by the abundance of olivine and orthopyroxene relicts in a serpentine matrix. Minor and accessory phases are clinopyroxene, chlorite, carbonates, amphiboles and spinels. The mineralogical content of the peridotites, as determined by the Rietveld method, is presented in Figure 4. Olivines occur in roundish and intensely fractured relics within the serpentine matrix with typical mesh texture ( Figure 5A-D,F,H). Irregular patches of a Fe-rich olivine are also present in some of the samples ( Figure 5B). Orthopyroxenes occur as fractured grains with equigranular shape ( Figure  5C,F,J). Serpentine is present partially replacing olivine along rims and as an interstitial phase ( Figure  5A-I). Fe-rich serpentine, when present, forms thin irregular rims bordering serpentine grains ( Figure  5D,G). The most abundant carbonate phase is magnesite, forming colorless anhedral grains ( Figure  5A,E). Three varieties of spinels can be observed: large unaltered chromite grains ( Figure 5F,I), altered ferrian chromite grains with a characteristic spongy texture ( Figure 5G) and small magnetite grains mainly bordering olivine crystals ( Figure 5D,H). Accessory phases are clinopyroxenes, occurring as colorless intensely fractured grains ( Figure 5I,J), rare amphiboles ( Figure 5G), calcite and dolomite ( Figure 5I). Listvenites, found in association with peridotites, show intense alteration and are crosscut by magnesite veins ( Figure 6A). Listvenitized samples show variable degrees of carbonation and silicification and different mineralogical associations. Their recurrent mineralogical association is quartz + carbonates. The most abundant carbonate is dolomite, followed by calcite and minor magnesite. Biotite mica occurs in some samples as a minor listvenite phase. Relict phases, witnesses of the ultramafic protolith are partially spinel (chromite and Fe-chromite), serpentine, chlorite, and amphibole. These minerals are also indicative of two processes predating or co-occurring with listvenitization: serpentinization and chloritization. Quartz is the most abundant phase. It forms aggregates of roundish crystals characterized by wavy extinction ( Figure 6B). Plagioclase, when present, forms euhedral crystals, crossed by carbonate-filled fractures ( Figure 6C). Dolomite forms euhedral crystals ( Figure 6D), from colorless to rosy and slightly pleocroic, with typical polysynthetic geminations that appear as lamellae oriented according to the short diagonals of the rhombohedral cleavage. Serpentine is present either in small veins, as aggregates of pseudomorphic crystals, or bordering quartz crystals in Fe-rich variety greenalite ( Figure 6D). Opaque phases are mainly Crspinels with little to absent Fe-chromitization ( Figure 6D). Sample W3 is constituted by a serpentine matrix with abundant plagioclase, carbonates and amphiboles. Plagioclase forms euhedral crystals often presenting the characteristic albite polysynthetic twinning ( Figure 7A,B,D). Magnesite is the prevalent carbonate and forms colorless anhedral grains ( Figure 7A,B). Serpentine is present in the aggregates of pseudomorphic crystals with a lamellar colorless to pale green structure, with low relief and low birefringence ( Figure 7A,B). The few primary silicate relicts are mainly represented by rounded and intensely fractured olivine relics ( Figure 7C).

Bulk-Rock Geochemistry
Bulk-rock compositions of the samples are reported in Table 2.

Major Elements
Listvenites are characterized by high variability in major element composition, with SiO2 ranging between 64.9 and 72.0 wt.%, Al2O3 from 0.1 to 10.7 wt.%, MgO from 5.8 to 8.4 wt.%, CaO from 3.9 to 7.4 wt.%, and loss of ignition (LOI) from 6.9 to 12.5 wt.%, depending on the relative abundance of mineral phases. Sample W17, from Ugo mine, has the highest SiO2 and the lowest Al2O3 contents. Peridotites present lower variability in their major element composition, with SiO2 ranging from 35 to 43.6 wt.%, Al2O3 from 0.2 to 0.7 wt.%, and MgO from 40.7 to 45.2 wt.%. Only LOI shows a much wider range of values that appear to be correlated to the degree of serpentinization of the samples. Lowest LOI values for least serpentinized samples is 3 wt.%, whereas the highest value in the strongly serpentinized sample is 14.4 wt.%. Sample W3 has a bulk rock geochemistry comparable to peridotites, at least concerning major elements. SiO2 content is 40.6 wt.%, Al2O3 is higher than in peridotites, with a value of 2.7 wt.%, and an MgO concentration of 33.3 wt.%, intermediate between peridotites and listvenites.

Trace Elements
All samples have generally low to very low trace elements contents (mostly below detection limits). All samples have relatively high contents of V, Cr, Co, and Ni, due to their ultramafic origin. Listvenites present high levels of Ga, Rb, Sr, Y, Zr, Nb, Mo, Ba, U, Th and rare earth elements (REE) with respect to peridotites. Rare earth elements (REE) contents are up to 1000 times higher in listvenites (up to 157.8 ppm) than in peridotites (<0.14 ppm). REE patterns normalized to primitive mantle are shown in Figure 8. Listvenites W1 and W2 and altered serpentinite W3 present an enrichment in LREE (Light REE) with respect to HREE (Heavy REE), also showing a slight negative Eu anomaly. W17 listvenite has a much lower REE content and less pronounced LREE enrichment, limited to La only. REE profiles could not be plotted for the other samples which were generally below detection limits, something expected for ultramafic ophiolitic rocks.

CIPW Normative Mineralogy
Normative mineralogy calculations were performed on partially serpentinized harzburgite rock samples in order to produce an idealized mineralogical association starting from whole rock analyses; the results are shown in Table 3. The application of these calculations was made to classify the studied samples; these samples have undergone various metasomatic events and so the current mineral assemblage cannot help under these circumstances. FeO was calculated from the total Fe2O3 using the relation: FeO = Fe2O3 × 0.8998. In all the samples, olivine is the most abundant mineral, ranging between 62.9 and 86.1 wt.%. Orthopyroxene is the second mineral in order of abundance, ranging between 0.6 and 26.6 wt.%. Other minor components are clinopyroxene (0.7-3.9 wt.%) and chromite (0.6-1.0 wt.%). The amount of plagioclase up to 2.2 wt.% is probably an artifact due to the limitations of the calculating formulae. From the normative mineralogic calculations, four samples (W5, W6, W7 and W11) were classified as dunites, while the remaining samples were harzburgites. However, due to the partial serpentinization of the samples, the differences in olivine and pyroxene abundances could partially be an artifact. As serpentine has roughly the same SiO2 amount as olivine, but less SiO2 amount than pyroxenes, the replacement of pyroxenes by serpentine results in the CIPW calculation in a systematic overestimation of olivine at the highest LOI%.

Serpentinized Peridotites
Peridotites are characterized by abundant olivine and orthopyroxene relicts in a serpentine matrix. THE mineral chemistry of major and minor phases is reported in Table 4

Serpentinite with Clay Minerals
Sample W3 has a peculiar mineralogical composition (Table 6), being a serpentinite with abundant plagioclase and clay minerals. The dominant mineral phases are serpentine, plagioclase, chlorite and amphibole, whereas clay minerals are mainly vermiculites.
Plagioclase (An43.5) has a composition comparable to the one in listvenites. Amphibole is characterized by an homogeneous composition, comparable to those of actinolites in listvenites. MgO contents range between 20.56 and 21.30 wt.%, FeO between 4.59 and 5.24 wt.%, CaO between 11.36 and 12.03 and Al2O3 between 3.02 and 3.31 wt.%.

Discussion
Based on the mineralogical assemblage, bulk rock geochemical data and EMPA analyses, we can trace the succession of events that affected the Rachoni massif peridotites. The peculiar mineral assemblage is indicative of several metasomatic events affecting the area. After a first stage of serpentinization, CO2-rich fluids percolated through the rocks replacing the mineralogical assemblage with a suite of new formed minerals. Locally, the carbonation of ultramafic rocks led to the formation of listvenites ( Figure 6A), with the precipitation of quartz and K-mica along with carbonates.

Genesis of the Protolith
The composition of unaltered chromites in the Rachoni peridotites is variable, with Cr2O3 contents between 34.53 and 53.90 wt.%, Al2O3 between 14.68 and 34.63 wt.% and TiO2 up to 0.33 wt.%. Their composition is comparable to the one of spinels plotting within the ophiolite compositional field ( Figure 9A,B). Rachoni peridotite spinels are characterized by relatively high XCr and low XMg, and plot in the forearc peridotite genetic field ( Figure 9C). Such character is in agreement with previous works on chromitite spinels from Gerakini [23] and Vavdos [24] ophiolites, that are within the podiform chromitites field plot in the area of low XMg. Compositional fields for abyssal peridotites [25][26][27], forearc peridotites [26] and podiform chromitites [23,24,28] are shown for comparison.

Fe-Chromitization
Alteration into ferrian chromite in the Rachoni peridotites is unevenly distributed and limited to a few samples. The general ferrian chromitization trend shows a depletion in Al, with a relative enrichment in Cr and Fe 3+ , coupled with an important Mg depletion ( Figure 10). Ferrian chromitization is generally confined to chromite rims, and only in rare cases it affects the whole chromite grain. Cr-chlorite development, which is generally associated to ferrian chromite alteration [29], is absent in our samples. The uneven distribution of ferrian chromitization, not correlated to that of listvenites and not as widespread as serpentinization, suggests that it can be related to the limited circulation of low-T fluids.

Serpentinization
Textural and petrographic observations indicate that serpentinization is the first alteration stage of the peridotites. The starting hydration reaction occurring at the first stages of serpentinization, at temperatures of ~400 °C, is the conversion of olivine to serpentine following reaction (1) [30]: According to the XRD data and petrographic observations, we can infer that peridotite samples from the Rachoni mine underwent serpentinization at high silica activity, which prevented the formation of brucite, which was not detected.
Orthopyroxene, as well as olivine, is also involved in hydration reactions. The reaction of orthopyroxene with silica-rich and water-rich fluids results in the precipitation of serpentine and talc, following reaction ( The development of Fe-rich hydrated phases during serpentinization, observed in some samples, suggests that during initial serpentinization, iron partitions from olivine and pyroxene into lizardite and brucite, as observed also in other harzburgites [31], according to reaction (3) With proceeding serpentinization, partial iron release from early formed Fe-hydrated phases led to the formation of magnetite ( Figure 5D) [30,32,33].

Fe Inclusions in Olivine
During the examination with SEM of sample W5, a very rare phenomenon was observed: olivine grains had inclusions (Figure 11), which the EDS and WDS analyses confirmed consisted of native Fe. The study of the mineral chemistry of olivines from sample W5 (Table 7) did not reveal any notable differences between the olivine grains hosting Fe inclusions and those not. All olivines were found with a high Fo content; Fo shows a constant value of 0.90 (Fo0.90Fa0.10). The NiO content of olivines is rather low (<0.44 wt.%) and no Ni alloys were found during the microscopic study of the sample. Native Fe has been found in ophiolitic rocks such as serpentinized dunites, peridotites, serpentinites and chromitites [34]. In most cases, native Fe presence is associated with serpentinization, e.g., [35][36][37][38][39]. Sakai and Kuroda [39] suggested that during serpentinization, Fe is released from olivine forming magnetite; hydrogen is also produced during serpentinization as indicated by the reaction [35]: 6 Mg1.5Fe0.5SiO4 + 7H2O = 3Mg3Si2O5(OH)4 + Fe3O4 + H2 (4) Hydrogen could theoretically reduce magnetite to Fe, but it is uncertain if the low temperatures during serpentinization would allow magnetite to be affected. Furthermore, native Fe in the aforementioned cases [35][36][37][38][39] is found in serpentine veins; in our sample, Fe was found as inclusions in olivine.
In order to explain the Fe inclusions in olivine, we propose that they could have segregated from magmas prior to or along with the crystallization of olivine. The segregation of these inclusions, although a minor amount, could result in Fe depletion in the evolved magmas [40]. Nevertheless, subsolidus equilibration between olivine and chromite (Mg 2+ diffusion from chromite to olivine and Fe 2+ from olivine to chromite at a subsolidus temperature) is the most probable factor controlling the olivine chemical composition and its high Fo content.

Carbonation
Carbonation of Rachoni harzburgites is characterized by the precipitation of carbonate phases, mainly magnesite, calcite and dolomite. The paragenetic sequence of carbonated ultramafic rocks has been studied by several authors, who suggest that a serpentinite can be transformed into a talc-rich rock and/or into a quartz-carbonate-bearing rock by interaction with CO2-rich hydrothermal fluids [41,42], derived from seawater and the atmosphere [42]. When serpentinized peridotites come into contact with such fluids, they are partially or completely replaced by talc-magnesite and/or quartzmagnesite rocks [43][44][45][46]. Textural observations reveal that magnesite in the Rachoni harzburgites preferentially replaces serpentine and pyroxenes, as observed in other carbonated serpentinites [43]. The reaction suggested for this process involves the interaction of lizardite/chrysotile with CO2-rich fluids to form antigorite and magnesite [31,47].
Minor carbonate phases found in the Rachoni peridotites are Fe-rich magnesite, dolomite and calcite. Fe-rich magnesite, found only locally, results from the replacement of olivine instead of serpentine [42], whereas dolomite and calcite are the result of the alteration of clinopyroxene crystals [31] by the same CO2-rich fluids.

Plagioclase Formation
Sample W3 shows a peculiar mineralogical association that could be related to a further alteration stage before listvenitization. The XRPD and mineral chemistry data reveal the presence of serpentine, magnesite, amphibole, chlorite and talc, very similar to the association in partially serpentinized harzburgites, but also of plagioclase. While the major element content of sample W3 is similar to the one in harzburgites, the REE content is higher, more similar to the content found in the listvenite of the Rachoni complex. This enrichment is due to Na-rich percolating fluids, also responsible for the precipitation of plagioclase.

Listvenitization
Areas of strong listvenitization can be found within the ultramafic rocks of Rachoni mine, crosscut by later magnesite veins ( Figure 6A). The main effect is carbonation in the presence of Kbearing fluids, with the replacement of primary minerals by an assemblage of Mg-Fe(-Ca) carbonates, quartz and chromian muscovite or sericite [48]. There is a mineralogical distinction between talc-carbonate rocks and listvenites, which is the presence of potassium in listvenites [49]. However, listvenitization is also commonly used to describe carbonated ultramafic rocks, and the term has already been used to describe similar lithologies in the ophiolite mélange of Iti mountain, in central continental Greece [49].
Quartz is present in all listvenite samples but is the most abundant in sample W17. Carbonation reactions and the influx of CO2 can lead to quartz oversaturation and its co-precipitation with magnesite [31]. Quartz and magnesite can also result from the reaction of talc and CO2-rich fluids. The absence of talc in the Rachoni listvenites could indicate its complete consumption in this second reaction. In contrast to "classic" listvenite, K-bearing minerals in Rachoni are not Cr-muscovite or sericite, but biotite and K-feldspar. In Ugo mine (W17), no K-bearing phase was found.
REE patterns (Figure 8) in listvenitic samples are generally flat but show an enrichment in LREE with respect to HREE. The two listvenites from Rachoni mine are highly enriched in LREE, with a content similar to plagioclase-enriched sample W3, whereas sample W17 (Ugo mine), although containing the highest SiO2 content, is less enriched in rare earth elements. Rachoni listvenites and plagioclase-enriched serpentinite, moreover, present a small to accentuated Eu-negative anomaly, which is absent in listvenite W17.
The REE enrichment in samples W1 to W3 could have been inherited from the same circulating fluids, percolating through the rock along fractures and pores and enriching the host lithology in incompatible elements. The enrichment in LREE also suggests a strong contribution of carbonates, since CO3 −2 ions fractionate LREE relative to HREE [50].

Conclusions
According to the results of our study, the Gerakini ophiolite shows signs of having undergone various stages of metamorphism/matasomatism. The Fe inclusions observed in some olivine grains must have segregated from magmas prior to or along with the crystallization of olivine and the segregation of these inclusions, which could result in Fe depletion in the evolved magmas.
After the initial serpentinization event (i.e., ocean floor metamorphism), CO2-rich fluids circulated into these rocks resulting in the further alteration of these rocks and the formation of new mineral phases with magnesite being the most important one due its economic interest. Furthermore, locally, the rocks were affected by silicification, carbonation and clay alterations, resulting in the formation of listvenites and clay mineral phases in the serpentinites.
The serpentinization of the rocks resulted in the formation of serpentine, whereas brucite was not found in these samples; the latter along with the presence of olivine relics, implies that serpentinization stopped before the final stage of brucite formation.
The post-serpentinization circulation of CO2-rich fluids in these rocks resulted in the formation of antigorite and magnesite from the carbonation of lizardite and/or chrysotile and the formation of Fe-magnesite from the carbonation of olivine, whereas dolomite and calcite were formed by the effect of the same fluids into clinopyroxenes.
During carbonation, the local development of the listvenites, resulting in the depletion of Mg, Mn, HREE, and Pd, and to the enrichment in Si, Ca, Ba, Zr, LREE, in the listvenites; Ti, Cr, Ni and Fe have been conserved. Two different fluids were involved in the evolution of these rocks as inferred by the differences in the REE contents of listvenites from two different localities of the Gerakini ophiolites.