Ophiolitic Remnants from the Upper and Intermediate Structural Unit of the Attic-cycladic Crystalline Belt (aegean, Greece): Fingerprinting Geochemical Affinities of Magmatic Precursors

The ophiolitic rocks of the Attic-Cycladic crystalline belt are considered of key importance for understanding the Mesozoic tectonic evolution of the Aegean region. Unresolved questions concern their tectono-stratigraphic relationships across the region. The mode of occurrence of the Cycladic ophiolites varies, as they appear as: (a) dismembered blocks (olistoliths) within the supra-detachment units of Paros and Naxos; (b) mélange formations in the upper structural unit of western Samos and Skyros; and (c) meta-ophiolitic mélange in the Cycladic Blueschist Unit (CBU) from central Samos. The trace element geochemistry and Pb-Sr-Nd isotopes of the mafic ophiolitic rocks indicate four geochemical groups: (a) the upper unit metabasites from Paros and western Samos (Kallithea) display an evolved basaltic composition (Mg# 40.2–59.6), with low Zr/Nb values (5–16) and high Ce/Y values (1.3 to 2.6) compared to MORB, indicating island-arc tholeiite affinities; (b) Naxos upper unit metabasalts show spider diagrams patterns indicating ocean island basalt (OIB-type) affinities; (c) Central Samos metagabbros (CBU) are primitive rocks with Back-Arc Basin basalt affinities; (d) the Skyros metadolerites and Tinos (Mt Tsiknias) and S. Evia (CBU) metagabbros, cluster as a separate geochemical group; they exhibit high MgO values (>10 wt %), very low TiO 2 values (0.1–0.2 wt %), Y and Yb, and depleted trace element N-MORB normalized patterns, similar to volcanic rocks formed in modern oceanic fore-arc settings, such as boninites. A combination of the Pb-and Sr-isotopic compositions of Cycladic metabasites indicate that the Pb and Sr incorporated in the Cycladic ophiolites correspond to mixtures of magmatic fluids Naxos—irrespective of the structural unit—display chemical affinities of ocean floor peridotites formed in a supra-subduction zone. The characteristics of harzburgite relicts in Cycladic serpentinites and Skyros indicate a highly residual nature of the mantle source. Geochemical data from this study shed further light on the complex structure of the oceanic lithosphere from which the Cycladic ophiolites originated.


Introduction
Ophiolitic remnants are widespread in the Cyclades, southern Aegean, occurring in the upper portion of the tectono-stratigraphic stack (Figure 1) [1][2][3][4][5][6].The mode of occurrence of ophiolitic associations varies in different parts of the Cyclades (Figure 1b,c) [7][8][9].The typical mode of occurrence is within low-pressure metamorphosed mélanges in the Upper Cycladic Unit, as on the islands of Anafi [5], Ikaria [10], and Samos (Kallithea nappe) [11], interpreted as subduction-accretion units of Cretaceous age [12] (Figure 1c).In contrast to that is the centrl Cyclades ophiolitic upper unit, which is limited to the occurrence of serpentinized peridotites and minor mafic lithotypes [4].In the central Cyclades islands of Paros and Naxos, supra-detachment basins associated with metamorphic core complexes enclose olistoliths of the latter lithologies overlain by Miocene and Pliocene sedimentary rocks [12,13].Additionally, highly sheared serpentinites at the base of these basins mark the low angle detachment faults of Naxos-Paros, and their tectonic position is interpreted as representing a very weak horizon during progressive exhumation [14] (Figure 1c,d).
The initial age of formation of ophiolitic rocks of the upper unit of the central Cyclades is unknown, although in other Cycladic islands and Crete both Jurassic and Cretaceous age ophiolitic rocks have been reported [5,10,15].Additionally, dismembered metaophiolites occur in the intermediate structural unit enclosed as blocks within mafic mélange units of the Cycladic Blueschists Unit (CBU) [4,6,[15][16][17] (Figure 1b).Based on the geochemical characteristics and structural appearance, these metaophiolitic rocks in the Cycladic islands (e.g., Samos, Syros, Tinos: [15]) have been correlated to the Jurassic ophiolites of the Pindos basin (western ophiolite belts) [5].An important question arising is whether the late Cretaceous ophiolite remnants from some Cycladic Upper Unit occurrences were high-level parts of a long (~25 Ma) accretion-subduction process of which the Eocene blueschist ophiolite remnants are deep-level representatives [4].
In the Hellenides orogenic belt a major Tethyan-type ophiolite emplacement occurred in mid-to late Jurassic times [18].This ophiolite belt ('Eohellenic' ophiolites) is interpreted as originating from a Mesozoic Neo Tethyan oceanic basin, the Pindos Ocean, and subsequently thrust northeastwards onto the Pelagonian passive continental margin (Figure 1; [19,20]).Ophiolite remnants in the upper structural units of the Cycladic islands and Crete have been considered as the southernmost continuation of this belt, based mainly on their tectonic position and geochemical characteristics [5] (Figure 1b).However, it is not yet well understood whether all the Cycladic ophiolites correlate with the Middle Jurassic ophiolites of the Pindos ocean ("Eohellenic nappe") in western Greece [21], or if some occurrences correlate with the Upper Cretaceous (70-95 Ma) ophiolites of western Turkey, and Cyprus [18].Ophiolite rocks from a mélange formation of Skyros and other northern Aegean islands are interpreted, according to their tectonic position, as possible Cretaceous age obducted ophiolite from the internal Vardar ocean [5,22].
Based on their tectonic setting, an ophiolite is emplaced either from down-going oceanic lithosphere via subduction-accretion, or from the upper plate in a subduction zone through trench-continent collision [23,24].In the case of the Aegean ophiolites, both modes of emplacement occur, indicating a change of the tectonic setting in the Aegean region through the Mesozoic time [25][26][27][28][29].
In this paper, we examine in detail the chemical compositions of mafic and associated ultramafic rocks from the structurally upper unit of Paros-Naxos, Samos, and Skyros, as well as minor occurrences of metamafic rocks from the migmatitic gneiss dome of Naxos (Figure 2).Moreover, we make a comprehensive overview of published geochemical data of ophiolitic lithologies from typical ophiolite exposures of the upper structural unit (UU) of Tinos [15] and Ikaria [10], and compare them with the metamorphosed ophiolitic mélange units (CBU) from Tinos, southern Evia, Andros, Samos [4,15,28], and Naxos [29], in order to place the new data into a regional context.
With these new whole-rock geochemical data that include Sr-Nd-Pb isotopic data for metabasites and ultramafics, we aim to: (a) present a basic geochemical characterization from the ophiolite occurrences of different Aegean islands (Cycladic islands and Skyros); and (b) recognize distinguishing chemical characteristics in order to make regional correlations.This research paper is based on sampling and whole-rock (major and trace element) analytical data collected by [30] and augmented later on by isotopic data on selected rock samples.Blueschist Unit (CBU), the Cycladic Basement rocks, and the Pelagonian unit within the ACC is also shown (modified after [12] and references therein).(c) Schematic tectono-stratigraphic sequence of the main Cycladic units after [1,[12][13][14] and Samos nappe piles after [19,26].The stratigraphic location of the ophiolitic rocks is marked with numbers on the columns: 1-Naxos and Paros dismembered supra-detachment ophiolites; 2-Skyros melange; 3-Naxos amphibolites from the CBU; 4-western Samos (Kallithea) mélange; 5-central Samos meta-ophiolitic mélange (CBU).Age abbreviations: Mz-Mesozoic, Ks-Cretaceous; Tr-J-Triassic-Jurassic.(d) Structural cross-sections across the Cycladic metamorphic belt showing the main structural units and the location of ophiolitic rocks, modified after [9].Blueschist Unit (CBU), the Cycladic Basement rocks, and the Pelagonian unit within the ACC is also shown (modified after [12] and references therein).(c) Schematic tectono-stratigraphic sequence of the main Cycladic units after [1,[12][13][14] and Samos nappe piles after [19,26].The stratigraphic location of the ophiolitic rocks is marked with numbers on the columns: 1-Naxos and Paros dismembered supra-detachment ophiolites; 2-Skyros melange; 3-Naxos amphibolites from the CBU; 4-western Samos (Kallithea) mélange; 5-central Samos meta-ophiolitic mélange (CBU).Age abbreviations: Mz-Mesozoic, Ks-Cretaceous; Tr-J-Triassic-Jurassic.(d) Structural cross-sections across the Cycladic metamorphic belt showing the main structural units and the location of ophiolitic rocks, modified after [9].
In the CBU there are widespread tectonic slices of ophiolites, intercalated with marbles, metapelites, mafic volcanites and flysch [3].These exposures include the metamorphic rocks in Andros, Tinos, Sifnos, Syros, Naxos and, further north, the Lavrion Blueschist Unit [38] and the Styra-Ochi unit of southern Evia [34].In most studies these units are considered equivalents of the Pindos ocean and some are correlative with sub-ophiolitic mélanges (Figure 1b).
The age of meta-ophiolitic fragments (olistoliths) from the Upper Unit of the Attic-Cycladic crystalline belt and the Late Miocene supra-detachment sedimentary units is considered to be pre-Late Cretaceous according to the Cretaceous metamorphic ages of those units [32,33].However, an overview of published radiometric data by [39] indicates that both Jurassic and Cretaceous meta-ophiolites occur in the HP mélanges of the CBU, but it is still questionable whether those age groups occur within the same mélange formations or at different lithostratigraphic and/or tectonic levels (Figure 1c).

Local Geology
On Paros, ophiolitic rocks occur in the Marmara nappe, which in places are transgressively covered by Barremian age limestones (Figure 1c).This transgression is considered to represent the timing of emplacement of ophiolites from the Pindos ocean over the Pelagonian continental margins (Eohellenic event) [40].The lower part of the Marmara Nappe (Supra-detachment units of map, Eastern Paros in Figure 2a), consists mainly of serpentinites, which preserve relicts of peridotite and only minor blocks of metabasites (metagabbros, Figure 2a).
On Naxos, mafic-ultramafic rock associations from the Upper unit occur as: (a) blocks and slabs of variable size (hundreds of meters) of serpentinized peridotites containing lenses of metagabbros (Galanados, Figure 2b); (b) small occurrences of metagabbros found along the detachment plane; and (c) slices of brecciated metadolerites enclosed in a matrix of sandstone, conglomerate and chert (Agios Thaleleos, Figure 2b).Additionally, metamorphosed and foliated ultramafic rocks mark the contact between the upper and lower structural units [16].These metamorphosed ultramafic rocks occur as elongated bodies, and often they are marked by the formation of monomineralic zones at their contact with the surrounding gneisses and schists.Serpentinites from two stratigraphic horizons are sampled in this study: (a) the main horizon of meta-ultrabasic bodies of Naxos dome at the boundary of leucogneiss-core with the Lower Series [29], from the localities of Agia and Kourounohori; and (b) the Upper (supra-detachment) Unit serpentinites from the localities of Galanados and Agios Thaleleos (Figure 2b).The aim was to evaluate the similarities/differences between samples from different structural units of the same islands.[8,25]; detailed geological maps of the (c) Samos (modified after [26,40]); (d) Skyros by [27]).Numbers indicate sample locations.
On Samos island six subunits have been distinguished by previous authors (Figures 1c and 2c).These include the metamorphosed Selçuk, Ampelos, Agios Nikolaos, and Kerketas nappes, and the uppermost Kallithea nappe in western Samos [26,40], Figure 2c.An ophiolitic mélange has been described at the uppermost unit of Kallithea nappe (western Samos, Figure 2c).This unit is described as chaotic mixture of pillow basalts, relict peridotites, serpentinites, conglomerates, radiolarites, clastic sediments, and limestones containing ammonitico rosso, having a thickness of 200 to 300 m [11].The size of mafic and ultramafic blocks within this chaotic formation ranges from few centimeters to tens of meters.Moreover, a metamorphosed mélange is described in the Selçuk nappe, at central Samos (Figure 2c).The Selçuk nappe has been correlated to the upper most part of the Cycladic Blueschist Unit [40].For the purpose of this study we have sampled ophiolitic occurrences from the Upper Unit (Kallithea nappe) as well as central Samos metaophiolitic mélange.
On Skyros island, central Aegean (Figure 2d), an ophiolitic mélange overthrusts the Palaeozoic-Middle Jurassic Pelagonian basement rocks and is considered as part of the Eohellenic nappe [5].K-Ar dating on a gabbroic rock (diorite) from Skyros mélange determined a Lower Cretaceous age (125 ± 8 Ma), which has been assigned to the formation age of the magmatic protolith [5].Ophiolite lithologies outcrop in the central-northern part of the island (Figure 2d) and consist of highly deformed serpentinite, massive gabbroic rocks, dolerite dykes, and massive basaltic/andesitic lavas and ophicalcite [41].In the framework of this study, selected mafic and ultramafic rocks from the Skyros mélange are compared with the structurally correlative Upper Cycladic Unit ophiolitic rocks.

Analytical Techniques
Thirty-three rock samples (13 metabasites and 20 serpentinites) were analyzed for major and trace elements.Sample locations are provided in Supplementary Table S1.Of these, 11 samples-selected as representative for each lithological type-were analyzed for Sr and Nd isotope compositions.Lead isotope ratios were measured on representative samples from Naxos, Paros, and Samos metabasic rocks.Lead isotope analyses of the metaperidotites samples were not possible to be determined due to the very low U and Pb concentrations in these samples, except for two samples from Naxos (NA-6, NA-14).
Major and trace element concentrations of the rock samples were determined by the XRF technique at Leicester University (UK).Major elements were determined on glass beads fused in a gold-platinum crucible at 1100 • C from ignited powders, and trace elements were determined on pressed powder pellets (see Supplementary Materials).Major elements were determined using Rh anode excitation, and the trace elements Sc, V, Cr, Ni, Zn, Ga, Rb, Sr, Y, Zr, Nb, La, Ce, Nd, and Ba, on pressed powder briquettes, using either Rh or W excitation, both on the Phillips PW1600 X-ray fluorescence spectrometer (Phillips electronics, Eindhoven, The Netherlands).Details of operating conditions and techniques for this instrument are given in [42].Although the concentrations of some of the trace elements (particularly Nb, La, Ce, and Th) were near or below the limit of detection by XRF techniques in many samples, a number of key, relatively immobile, elements were able to be determined with sufficiently high precision by XRF methods to enable clear deductions concerning the petrogenesis of the Cycladic ophiolitic rocks.
Radiogenic isotopic compositions of Sr, Nd and Pb were measured on 12 whole rock (six metaperidotite and six metabasites) samples from Naxos, Paros, Samos, and Skyros at the University of Geneva, Switzerland.Between 100 and 120 mg of whole rock powder were dissolved during seven days in Savillex ® Teflon (Savillex, Minnetonka, MN, USA) vials using 4 mL of concentrated HF and 1 mL of HNO 3 14 M, at a temperature of 140 • C and with the help of ultrasonication for 30 min twice a day.Subsequently, samples were dried and re-dissolved for three days (also with 30 min ultrasonication twice a day) in 3 mL of HNO 3 14 M, and dried again.Sr, Nd, and Pb were then separated using cascade columns with Sr-Spec, TRU-Spec, and Ln-Spec resins, respectively, according to a protocol modified from [43].Finally, the material was redissolved in 2% HNO 3 solutions and ratios were measured using a Thermo Neptune PLUS Multi-Collector ICP-MS (Thermo Fisher Scientific, Waltham, MA, USA) in static mode.Ratios used to monitor internal fractionation were: 88 Sr/ 86 Sr = 8.375209 for the 87 Sr/ 86 Sr ratio, 146 Nd/ 144 Nd = 0.7219 for the 143 Nd/ 144 Nd ratio, and 203 Tl/ 205 Tl = 0.418922 for the three Pb ratios (a Tl standard was added to the solution).External standards used were SRM987 ( 87 Sr/ 86 Sr = 0.710248, long-term external reproducibility: 10 ppm), JNdi-1 ( 143 Nd/ 144 Nd = 0.512115; [44] long-term external reproducibility: 10 ppm), and SRM981 [45] for Pb (long-term external reproducibility of 0.0048% for 206 Pb/ 204 Pb, 0.0049% for 207 Pb/ 204 Pb and 0.0062% for 208 Pb/ 204 Pb). 87Sr/ 86 Sr, 143 Nd/ 144 Nd and Pb isotope ratios were further corrected for external fractionation (due to systematic differences between measured isotope ratios and those proposed for SRM987, JNdi-1 and SRM981, by [44][45][46] respectively]) by −0.039 , +0.047 and +0.5 amu, respectively.Interferences at masses 84 ( 84 Kr), 86 ( 86 Kr) and 87 ( 87 Rb) were corrected by monitoring 83 Kr and 85 Rb, 144 Sm interference on 144 Nd was monitored on the mass 147 Sm and corrected by using a 144 Sm/ 147 Sm value of 0.206700, and 204 Hg interference on 204 Pb was corrected by monitoring 202 Hg.
All the analyzed metabasites have low U/Pb (ranging from 0.1 to 0.5) as well as Th/Pb ratios, and therefore no corrections due to time-integrated decay of U and Th were made because the addition of radiogenic lead is considered to be insignificant (Table 1).An exception to this is Naxos metabasite NA-3, which displays an elevated Th/Pb value of 1.3, and the serpentinite NA-14, which is characterized by slightly more radiogenic compositions, suggesting a higher U/Pb ratio.

Naxos
The metaophiolitic outcrops of Naxos comprise mafic and ultramafic rocks.Mafic rocks are classified as (a) metagabbros and (b) metavolcanics.Ultramafic rocks comprise dunites and serpentinites.

Dunites and Serpentinites
Dunites show a rather constant mineralogical composition, comprising mainly of olivine with subordinate orthopyroxene and spinel.These fine-grained rocks display mosaic and fine equigranular textures composed of small, polygonal and occasionally unstrained olivine neoblasts (Figure 3c).Spinel is dark red and display anhedral grains heterogeneously distributed in the olivine matrix.Locally, spinel grains exhibit narrow rims of magnetite.The degree of serpentinization varies between 10% and 30%.
Serpentinites display in general mesh texture and consist essentially of serpentine minerals, mainly antigorite with minor chlorite, calcite and talc.The only preserved primary phase is Cr-spinel, which occasionally is totally replaced by magnetite.Pseudomorphic prismatic plates of bastite (lizardite pseudomorphs after pyroxenes), which sometimes mimic the cleavage of earlier pyroxene, are common.

Paros
The collected metaophiolitic mafic and ultramafic rock samples from Paros are classified as: (a) metagabbros and (b) serpentinites.

Metagabbros
The metagabbros display relict granular to granoblastic textures and consist of 30-40 vol % plagioclase, 20-30 vol % amphibole, 5-10 vol % chlorite, ~5 vol % epidote, minor titanite, and accessory magnetite, quartz, and carbonates.Plagioclase laths are replaced by a microcrystalline aggregate of sericite, albite, calcite, epidote and clinozoisite, forming pseudomorphs due to sericitization and saussuritization.Amphibole typically occurs as thin prisms and fibres of actinolite and tremolite, replacing clinopyroxene.Relicts of clinopyroxene occur as inclusions within some amphibole prisms, or more rarely as large xenomorphic crystals up to 5 mm in diameter, surrounded by continuous thin rings of actinolite and chlorite.These metagabbros show better preserved igneous textures relative to the amphibolites of Naxos, and have not acquired a visible metamorphic fabric.

Serpentinites
Serpentinites exhibit mesh and locally hourglass textures.They consist essentially of serpentine minerals, mainly antigorite with minor chlorite, calcite and talc.The only preserved primary phases are relicts of olivine and Cr-spinel, which occasionally is totally replaced by magnetite.Pseudomorphic prismatic plates of bastite (lizardite pseudomorphs after pyroxenes), which sometimes mimic the cleavage of earlier pyroxene, are common (Figure 3d).Sparsely, veinlets of calcite, magnesite and talk crosscut the serpentinites.
Serpentinites exhibit mesh textures.They consist essentially of serpentine minerals, mainly antigorite and chrysotile with minor chlorite, calcite, and talc.Occasionally relicts of olivine and Cr-spinel are presented.Pseudomorphic prismatic plates of bastite (lizardite pseudomorphs after pyroxenes), which sometimes mimic the cleavage of earlier pyroxene, are common.Sparsely, veinlets of calcite, serpentine, and talc crosscut the serpentinites.

Skyros
The collected metaophiolitic mafic and ultramafic rock samples from Skyros are classified as: (a) rodingites and (b) serpentinites.

Rodingites
Rodingites mainly display a blasto-subophitic texture with local relicts of primary subophitic texture, which is strongly indicative of a doleritic protolith (Figure 3f).They contain 30-50 vol % clinopyroxene, 20-30 vol % chlorite, 10-20 vol % garnet with minor calcite, titanite, and ilmenite.Clinopyroxene and chlorite formed probably after actinolite but they retain a crystal habit of an older clinopyroxene.Garnet (grossular) occurs in between the clinopyroxene grains and clearly pseudomorphs original plagioclase (Figure 3f).It comprises weakly birefringent, frequently tiny, micron-size crystals.Its small size and anisotropic character make it difficult to distinguish from vesuvianite, which also commonly appears in rodingites.Fe-Ti oxides with a dusty appearance are dispersed throughout the samples.

Serpentinites
The serpentinites are composed essentially of antigorite, chrysotile, and lizardite, with subordinate amounts of carbonates, chromite, magnetite, magnesite, talc, tremolite, and chlorite.Antigorite is the dominant mineral and occurs as interpenetrating and interlocking elongate plates.Sometimes, the antigorite flakes and blades are arranged in roughly parallel aggregates.Serpentine, which developed as an alteration after olivine, usually has a mesh texture while that formed after orthopyroxene has a bastite texture.These textures indicate dunite and harzburgite parent rocks.
The opaque minerals in the studied serpentinites are represented mainly by dark-brown spinel and magnetite.Spinel occurs as disseminated subhedral crystals and/or irregular grains.Alteration of spinel along grain borders and cracks is very common.Magnetite occurs as anhedral grains in the interstitial spaces between the serpentine minerals or as veinlets.

Whole-Rock Geochemistry
The chemical compositions of analyzed samples are presented in Table 1.Sample locations are indicated in Figure 2a-d.Standard geochemical diagrams have been utilized in order to classify the investigated former magmatic rocks (Figure 4).
Geosciences 2017, 7, 14 15 of 32 opaque minerals in the studied serpentinites are represented mainly by dark-brown spinel and magnetite.Spinel occurs as disseminated subhedral crystals and/or irregular grains.Alteration of spinel along grain borders and cracks is very common.Magnetite occurs as anhedral grains in the interstitial spaces between the serpentine minerals or as veinlets.

Whole-Rock Geochemistry
The chemical compositions of analyzed samples are presented in Table 1.Sample locations are indicated in Figure 2a-d.Standard geochemical diagrams have been utilized in order to classify the investigated former magmatic rocks (Figure 4).

Upper Unit Metabasites
Studied ophiolitic metabasites from the Cycladic islands and Skyros span basalts, gabbros, and minor basaltic andesites and dolerites, ranging from 39 to 54.5 wt % SiO2.Samples from the Cycladic islands Paros, Naxos, and western Samos (Kallithea) include moderately to highly evolved basaltic

Upper Unit Metabasites
Studied ophiolitic metabasites from the Cycladic islands and Skyros span basalts, gabbros, and minor basaltic andesites and dolerites, ranging from 39 to 54.5 wt % SiO 2 .Samples from the Cycladic islands Paros, Naxos, and western Samos (Kallithea) include moderately to highly evolved basaltic rocks with 2.3 to 6.4 wt % MgO and Mg-number (Mg#) values from 40.2 to 59.6 (Figure 4) (Table 1).Only Skyros metabasites include primitive basaltic composition rocks with 10 wt % MgO at 50 wt % SiO 2 , Ni > 200 ppm but Co < 50 ppm, thus indicating early olivine fractionation (Table 1) (Figure 4).TiO 2 contents of the relatively more primitive basaltic rocks (Paros: PA-9; Samos: SA-8, SA-2) are low, i.e., 0.6 to 0.7 wt % when compared to N-and E-type MORB compositions, and resemble typical island arc tholeiites or some Back-Arc Basin basalts, e.g., East Scotia (Table 1) [47,50].In fact, there is a large range of TiO 2 content in the studied rocks that can be attributed to various processes, including fractionation and post-magmatic alteration.
A common feature in all the mafic rocks is the high potassium concentration.Even the less evolved rocks of gabbroic composition from Paros, SE Naoussa, contain up to 2 wt % K 2 O (Table 1).Moreover, the Pb concentrations of all the Cycladic metabasites exceed 1 ppm (i.e., 1-17 ppm), which is high for typical MORB [51].Notwithstanding, values up to 3 ppm may be normal for Back-Arc Basin basalts, e.g., in lavas from Southern Mariana Trough [52], yet most Cycladic rocks display extremely high values, up to 18 ppm, which clearly indicate secondary alteration (Table 1).The petrographic characteristics of the Cycladic metabasites (UU) document ocean-floor hydrothermal alteration (hydration), especially associated with extensive sericitization of plagioclase, which justifies the mobility of alkalis and lead.This is further attested in the spread of 207 Pb/ 206 Pb diagram (see Section 5.4).
Few detailed studies have been made of the Cycladic upper unit metabasites and the discussion in this section is based on the dataset of [15] on Tinos metagabbros (UU) and Ikaria ophiolitic rocks (Faros mélange) [10].In order to elucidate the petrogenesis of the Cycladic rocks, we rely on the relatively immobile key elements.The Total Akali-Silica (TAS) proxy diagram of [51] is considered to provide an effective discrimination of altered ophiolitic rocks, avoiding the wide dispersion of the IUGS TAS [53,54].The Nb/Y ratio is used as a proxy for total alkalis (Na 2 O + K 2 O) and Zr/Ti as a proxy for SiO 2 .In Figure 5 most Cycladic metabasites span the field of tholeiiteic basalts and basaltic andesites.A group of samples from the Naxos Upper Unit (metabasalts, basaltic andesites) and two amphibolites from Ikaria plot in the field of alkali basalts (Nb/Y > 1).Sample SA-8 (meta-andesite) plot as an outlier due to extremely high Nb/Y value (7.85), indicating extreme alkalinity of the magma and possibly a different magmatic type.
High precision Zr-Y determinations are a particularly sensitive method of distinguishing different magma batches in a magmatic suite.Both elements remain essentially incompatible within basaltic magma during fractional crystallization.Hence, a single magma batch will retain an almost constant Y/Zr ratio [42].In Figure 6a there is a reasonable linear correlation between Zr and Y but the Zr/Y ratio is variable.The majority of upper unit metabasites exhibit mainly a tholeiitic affinity (one linear trend), whereas a small group of samples from Naxos shows a different linear trend with higher Zr/Y (e.g., 5.2-5.4).Such relationship cannot be attributed to crystallization of a single magma batch.Similar geochemical systematics have been reported from the sheeted dyke suite and pillow lavas of Troodos complex [42] and basalts from modern Back-Arc Basins [55].

Naxos Amphibolites and Central Samos Metagabbro (CBU)
Two amphibolite samples from the Naxos gneiss dome (NA-7: Kourounohori, NA-9: Agia) show very low TiO 2 contents (TiO 2 : 0.1-0.2wt %) and extremely high MgO (19-27 wt %) for mafic rocks, and very low Al 2 O 3 at SiO 2 ≥ 52% (Table 1).According to the extreme composition of this group of rocks and considering the microscopic observations of thin sections, NA-9 ampibolite is interpreted as a metamorphosed tholeiitic basalt, whereas NA-7 is probably a metasomatised cummulate rock, in particular a former pyroxenite (NA-7).The bulk-rock geochemical characteristics of these amphibolites are comparable to amphibolites from metamorphic soles beneath the ultramafic ophiolites in northern Evia [8].A metagabbro (SA-11) from a metamorphosed ophiolitic mélange in central Samos-Myli, displays primitive basaltic composition with 9 wt % MgO at 50.3 wt % SiO 2 .Ni and Co contents of the latter sample is lower than 250 ppm and 50 ppm, respectively, thus indicating early olivine fractionation [49] (Figure 4a,b).In Figure 6, SA-11 metagabbro plots in the field of MORB tholeiites whereas Naxos amphibolites plot in the field of island arc tholeiites.

Multi-Element Geochemical Patterns
For representation of multi-element spider diagram patterns (Figure 7), only immobile elements are included and ordered in increasing compatibility, in the standard form after [51].Typical compositions of basaltic rocks from modern well defined settings are shown for direct comparisons.In general, relatively immobile elements under seafloor hydrothermal alteration of oceanic crust are considered the High-Field-Strength Element (HFSE) group (e.g., Zr, Nb, Ta, Ti), the REEs, Mg, Y, and especially the transition metals (Sc, V, Cr, Fe, Co, Ni, and Zn) [60,61].On the contrary, the Large Ion Lithophile Elements (LILE) K, Rb, Ba, Sr, and U are considered readily mobile under low-temperature alteration of ocean floor basalts and basaltic glass, and they are commonly added through the process of the formation of clay [62,63].
In these diagrams, N-type MORB [48] is selected as a normalized factor because the resulting flat N-MORB pattern is a useful reference regarding tectonic setting.
The Cycladic samples can be classified into three geochemical groups (Figure 7a-g): (a) Paros metagabbros (Figure 7a) (SE Naoussa-Upper Unit) and Samos (Kallithea nappe) metavolcanites (Figure 7b) display a Zr-Y flat pattern and higher Th/Y ratio relative to MORB or OIB; these characteristics are typical of modern arc lavas from complex tectonic settings (C-MORB, Taitao lavas, Chile; [51] and references therein).The gradient of the patterns for the Paros-W.Samos samples is even higher than that for back-arc basin basalts (BABB) tholeiites (Figure 7a,b).Moreover, the (La/Nb) N ratios, which range from 1.2 to 2.35 (Figure 6b), emphasize the HFSE depletion of this type of magmas relative to LREE, which is another characteristic feature of arc magmas [64,65].(b) Sample SA-11 (metagabbro from CBU), from central Samos, displays a flat pattern typical for MORB (Figure 7b).Following the definition of N-type MORB of [66], this sample exhibits Ce/Nb > 2.5 but Ce N /Yb N = 1 (Table 1) and, combined with initial Sr > 0.7035, it is assigned to an enriched MORB category (E-MORB) rather than N-MORB.However, comparing sample SA-11 with a complete dataset from [15] (Figure 7f) it becomes clear that there is a whole range of compositions typical of BABB, similar to Troodos ophiolitic lavas [51].(c) Trace element patterns of Naxos metabasaltic rocks (NA-3, NA-15) share many characteristics with oceanic island basalts (OIB) (Figure 7c): low Zr/Nb (i.e., ~5) relative to N-MORB [67], high Nb/Y, and the negative gradient in the Ti-Yb(Y) part of the diagram, which are characteristics of enriched ocean-island basalts (OIB).This is distinct from the flat patterns for Ti-Yb(Y) that characterize both N-and EMORB compositions.(d) Naxos amphibolites (CBU) also display a high Th/Y ratio similar to Paros and western Samos metabasalts; however, they also have a strong depletion for most of the trace elements LREE and HFSE, compared to typical MORB patterns (Figure 7c).
(e) Skyros metagabbro shows a typical pattern of boninitic rocks (Figure 7d) from modern oceanic arc environments, such as the fore-arc volcanics of the Mariana-Izu-Bonin [68].One metabasalts and a metagabbro from the same area also show similar boninitic affinities (data from [41], Table 1).

Upper Unit Serpentinites and Serpentinized Peridotites
Major and trace element contents of 20 representative samples from Cyclades and Skyros have been analyzed and the results are listed in Table 1.

a. Cycladic serpentinites
Based on the olivine-orthopyroxene-clinopyroxene CIPW normative classification, most of the Cycladic serpentinite protoliths are assigned to harzburgites and minor olivine orthopyroxenite whereas only one sample from Skyros (SK-9) is assigned to the lherzolite category.However, published data for metaperidotites from the main horizon of Naxos gneiss dome [29] (Katzir et al., 1999) display a broader compostional range from lherzolites to harzburgites and rare dunites, whose variability is not shown in the upper unit peridotites of Naxos possibly due to the scattered occurrence of the upper unit serpentinites (Figure 8).In major element variation diagrams the Mg# is used as a fractionation index (Figure 9a-c).The Mg# values are high in all the Cycladic rocks, ranging from 84 to 91.5.SiO 2 , Al 2 O 3 , and CaO display negative correlation trends defined by the Paros-Naxos peridotites, which are similar to those found for the main horizon of metaperidotites from Naxos.These trends have been attributed to a successive depletion process caused by melt extraction from an original lherzolitic mantle source as indicated in Figure 9d [29].In the binary diagram of MgO/SiO 2 vs. Al 2 O 3 /SiO 2 (Figure 9d), the serpentinized peridotites and serpentinites of our study plot below the "successive magmatic depletion" line, as expected for moderately serpentinized peridotitic rocks [69][70][71], which is attributed to restricted loss of Mg during seafloor weathering after serpentinization [71,72].The Aegean serpentinites are comparable to peridotites of the western ophiolite belt (e.g., east Othris [73].
Geosciences 2017, 7, 14 20 of 32 Major and trace element contents of 20 representative samples from Cyclades and Skyros have been analyzed and the results are listed in Table 1.

a. Cycladic serpentinites
Based on the olivine-orthopyroxene-clinopyroxene CIPW normative classification, most of the Cycladic serpentinite protoliths are assigned to harzburgites and minor olivine orthopyroxenite whereas only one sample from Skyros (SK-9) is assigned to the lherzolite category.However, published data for metaperidotites from the main horizon of Naxos gneiss dome [29] (Katzir et al., 1999) display a broader compostional range from lherzolites to harzburgites and rare dunites, whose variability is not shown in the upper unit peridotites of Naxos possibly due to the scattered occurrence of the upper unit serpentinites (Figure 8).In major element variation diagrams the Mg# is used as a fractionation index (Figure 9a-c).The Mg# values are high in all the Cycladic rocks, ranging from 84 to 91.5.SiO2, Al2O3, and CaO display negative correlation trends defined by the Paros-Naxos peridotites, which are similar to those found for the main horizon of metaperidotites from Naxos.These trends have been attributed to a successive depletion process caused by melt extraction from an original lherzolitic mantle source as indicated in Figure 9d [29].In the binary diagram of MgO/SiO2 vs. Al2O3/SiO2 (Figure 9d), the serpentinized peridotites and serpentinites of our study plot below the "successive magmatic depletion" line, as expected for moderately serpentinized peridotitic rocks [69][70][71], which is attributed to restricted loss of Mg during seafloor weathering after serpentinization [71,72].The Aegean serpentinites are comparable to peridotites of the western ophiolite belt (e.g., east Othris [73].Compositional variation of the Cycladic serpentinites show a dominant chemical group with average low Al2O3 (0.3-1.5 wt %), TiO2 (0.01-0.04 wt %) and CaO (0.03-1.52 wt %), and low V (5-49 ppm), coupled with high Ni (2318 ppm), Co (81-123 ppm), and Cr (1935-3932 ppm), indicating a highly residual nature of the mantle source; they are comparable to the composition of harzburgites from Pindos ophiolite complex [74] and mantle wedge serpentinites from oceanic fore-arc settings [75] (Figure 9).A subordinate group from Paros-Naxos displays a relatively less refractory character testified by the higher Al2O3 (1.79-2.86wt %) and V (50-60 ppm), and lower Co (96-103 ppm) and Cr (2398-2677 ppm) concentrations.Two serpentinite samples, one from Samos (SA-14) and one from Skyros (SK-4), display unusually high Al2O3 values of 4.82 and 8.61, respectively (Figure 9b).Based on thin section observations and mineral analysis data, these values can be attributed to extensive chlorite intergrowth in serpentine.

b. Skyros serpentinites
Skyros serpentinites are classified mostly as hartzburgites, except from one sample as lherzolite and another as olivine websterite (Figure 8).Skyros ultramafic rocks display the highest Mg# (up to 91.5).Ni and Cr abundances are as high as 2720 ppm and 3932 ppm, respectively.

Sr-Nd Results
Sr and Nd isotope ratios were determined for 12 samples.The three peridotite samples failed to produce Nd isotopic data due the very low content of REEs (<1 ppm).The results are summarized in Table 2 and presented in Figure 10a,b.Initial Sr and εNd values have been recalculated at 160 Ma, which corresponds to the available mean ages of the Pindos Jurassic ophiolite, and at 80 Ma, for the case of metaigneous rocks from the CBU [15,[76][77][78].The Cycladic ophiolites exhibit initial 143 Nd/ 144 Nd values lower than 0.5130, which is the lower limit value of most N-type middle ocean ridge basalts (N-type MORB), and higher initial 87 Sr/ 86 Sr ratios than 0.7035 [66].These isotopic values, if primary, suggest that the source of these rocks was isotopically modified from an enriched mantle component (OIB).Generally, the initial Sr-Nd isotopic ratios of Cycladic metabasites display a large scatter considering the number of samples analyzed.In the Sr-Nd isotopic diagram of Figure 10a, a main, b.Skyros serpentinites Skyros serpentinites are classified mostly as hartzburgites, except from one sample as lherzolite and another as olivine websterite (Figure 8).Skyros ultramafic rocks display the highest Mg# (up to 91.5).Ni and Cr abundances are as high as 2720 ppm and 3932 ppm, respectively.

Sr-Nd Results
Sr and Nd isotope ratios were determined for 12 samples.The three peridotite samples failed to produce Nd isotopic data due the very low content of REEs (<1 ppm).The results are summarized in Table 2 and presented in Figure 10a,b.Initial Sr and εNd values have been recalculated at 160 Ma, which corresponds to the available mean ages of the Pindos Jurassic ophiolite, and at 80 Ma, for the case of metaigneous rocks from the CBU [15,[76][77][78].The Cycladic ophiolites exhibit initial 143 Nd/ 144 Nd values lower than 0.5130, which is the lower limit value of most N-type middle ocean ridge basalts (N-type MORB), and higher initial 87 Sr/ 86 Sr ratios than 0.7035 [66].These isotopic values, if primary, suggest that the source of these rocks was isotopically modified from an enriched mantle component (OIB).Generally, the initial Sr-Nd isotopic ratios of Cycladic metabasites display a large scatter considering the number of samples analyzed.In the Sr-Nd isotopic diagram of Figure 10a, a main, horizontal trend is recognized.This trend is defined by the metamorphosed gabbros of central Samos (CBU), and sample SA-11 (Myli location, this study) also plot along this trend together with a serpentinite from Samos and the Tinos (UU) metagabbros (Mt.Tsiknias) [15].This trend is characterized by a shift towards higher 87 Sr/ 86 Sr (ISr80 = 0.7062 to 0.708) at almost constant initial 143 Nd/ 144 Nd (εNd80 = +6.54 to +7.48), calculated for a presumed protolith age of 80 Ma.The Skyros metadolerite (SK-1) also plots together with the Samos and Tinos metabasites.In Figure 10b initial isotopic ratios are recalculated at 160 Ma which is the approximate age of the Pindos ophiolite complex ( [77]; and references therein]).In the latter diagram, Paros (PA-9) and Naxos (NA-3) upper unit metabasites plot close to Avdella mélange gabbros, and define a separate group showing elevated initial Sr (0.70413-0.70563)and distinctly lower εNd (+4.13-+0.32)isotopic values compared to the Samos field.Sr and Nd isotope compositions of metagabbros from Andros CBU (calculated for 160 Ma) can be clearly distinguished from the other Cycladic metabasites (Paros, Naxos, Tinos) samples by even higher 87 Sr/ 86 Sr = 0.70613 to 0.70694 ratios and lower εNd = −3.77 to −1.54 values (Figure 10b), strongly indicating significant crustal involvement from terrigenous sediment or arc crust [15,28].
Upper Unit ophiolites of Paros and Naxos show a strong coupling of the two isotopic systems towards more radiogenic Sr and less radiogenic Nd isotopic ratios.This second trend could possibly be attributed to the combined effect of contamination of the source from a subduction-derived component (hydrous fluids or melt) and seafloor hydrothermal alteration.

Lead Isotope Geochemistry
Lead isotopes, in combination with other isotopic systems (Sr-Nd), have been proven as an important source of information on the source reservoir of magmatic rocks [79].In this study, lead isotope compositions were measured on samples from metavolcanites and metagabbros from the Upper and intermediate Units of the Cyclades.Lead isotope compositions of the metabasites are given in Table 3.
The analyzed metabasites from the Cyclades (metagabbros, metabasalts, and metabasaltic andesites) display a wide range of lead isotope compositions for a relatively small number of samples, especially regarding 207   11a,b).Similarly, high whole-rock 207 Pb/ 204 Pb ratios (Figure 11a) have been reported for the hydrothermally altered parts of the Troodos ophiolite complex (greenschists and epidosites) [80].In the 207 Pb-206 Pb and 208 Pb-206 Pb diagrams, the Cycladic metabasites display a spread beyond the normal MORB field and Troodos glasses [80] and plot within the field of seawater, attesting the effect of seawater alteration on the Cycladic metabasites compositions.Sample NA-14 (serpentinized harzburgite) has a high 206 Pb/ 204 Pb ratio, probably due to time-integrated development of radiogenic lead caused by a higher U/Pb ratio compared to the other samples (not plotted in Figure 11).[15]. 87Rb/ 86 Sr and 147 Sm/ 144 Nd isotopic ratios were calculated on the basis of whole rock Rb/Sr and Sm/Nd ratios (concentrations from XRF analysis), using the applications of GCDKit 4.0 software [81] and equations of [82].Initial isotopic ratios and epsilon Nd values were calculated using the plug-in modules of GCDKit software.n.d: not detected, bdl: below detection limit.* 80 Ma is the reported age of Samos metagabbros from [15]. 87Rb/ 86 Sr and 147 Sm/ 144 Nd isotopic ratios were calculated on the basis of whole rock Rb/Sr and Sm/Nd ratios (concentrations from XRF analysis), using the applications of GCDKit 4.0 software [81] and equations of [82].Initial isotopic ratios and epsilon Nd values were calculated using the plug-in modules of GCDKit software.

Determining the Effect of Alteration on the Cycladic Metabasites
In general, the Cycladic metabasites exhibit K enrichment when compared to gabbroic rocks from recent oceanic crust.For instance, normal MORB compositions do not exceed 0.17 wt % K2O (e.g., [79].This enrichment in K2O is attributed to an addition of potassium (plus Rb and Pb) due to low-T seafloor alteration reactions in basalts, which is commonly observed within the oceanic crust ( [83]; and references therein]).Sericitization of plagioclase is a common feature in all the mafic lithotypes (metagbbros and metavolcanites).Similar K-enrichments are reported in gabbroic rocks and dolerite dykes from southern Aegean ophiolites (Crete, Rhodes, and Karpathos [21]).
Moreover, the 87 Sr/ 86 Sr ratio is considered a very powerful indicator of hydrothermal alteration in the oceanic crust, whereby fresh modern crust is unradiogenic ( 87 Sr/ 86 Sr = 0.7025-0.703)and modern seawater is radiogenic (0.7092) [79].Based on microscopic observations (e.g., Figure 3a) and the formation of secondary minerals, all the Cycladic upper unit samples show evidence of lowtemperature seafloor hydrothermal alteration at greenschist facies.Sr-Nd isotopic variations of Samos (and Tinos, data from [28]) samples display a characteristic shift towards higher 87 Sr/ 86 Sr (up to 0.708) at almost constant initial 143 Nd/ 144 Nd (i.e., εNd80 = +6.54 to +7.48) (Figure 10).This shift indicates a seawater contribution at hydrothermal conditions of upper extrusive oceanic crust, and is a common feature in Mediterranean ophiolite complexes (e.g., Troodos: [83]).Previous authors [83] suggested that interchange of Jurassic and Cretaceous seawater Sr with ocean crust Sr occurs during hydrothermal circulation and may buffer the Sr isotope composition of seawater.Nevertheless, seawater 87 Sr/ 86 Sr has changed with time; thus, before 80 My, seawater was less radiogenic (0.7076) than today and has increased since then [84,85].The value of 0.708 is the maximum value of 87 Sr/ 86 Sr found in the Cycladic metabasites.In this context, the initial 87 Sr/ 86 Sr of the metabasites reflects the mixing ratio of seawater-derived Sr and basalt Sr present in a sample.Thus, the significance of Sr isotope values for determining the characteristics of the magmatic protolith and tectonic setting is limited.From Figure 11a,b it is clear that the Pb-isotopic ratios for Samos and Naxos metabasites show wide variation, from close to MORB into the seawater field.The elevated 207 Pb/ 206 Pb and the

Determining the Effect of Alteration on the Cycladic Metabasites
In general, the Cycladic metabasites exhibit K enrichment when compared to gabbroic rocks from recent oceanic crust.For instance, normal MORB compositions do not exceed 0.17 wt % K 2 O (e.g., [79].This enrichment in K 2 O is attributed to an addition of potassium (plus Rb and Pb) due to low-T seafloor alteration reactions in basalts, which is commonly observed within the oceanic crust ( [83]; and references therein]).Sericitization of plagioclase is a common feature in all the mafic lithotypes (metagbbros and metavolcanites).Similar K-enrichments are reported in gabbroic rocks and dolerite dykes from southern Aegean ophiolites (Crete, Rhodes, and Karpathos [21]).
Moreover, the 87 Sr/ 86 Sr ratio is considered a very powerful indicator of hydrothermal alteration in the oceanic crust, whereby fresh modern crust is unradiogenic ( 87 Sr/ 86 Sr = 0.7025-0.703)and modern seawater is radiogenic (0.7092) [79].Based on microscopic observations (e.g., Figure 3a) and the formation of secondary minerals, all the Cycladic upper unit samples show evidence of low-temperature seafloor hydrothermal alteration at greenschist facies.Sr-Nd isotopic variations of Samos (and Tinos, data from [28]) samples display a characteristic shift towards higher 87 Sr/ 86 Sr (up to 0.708) at almost constant initial 143 Nd/ 144 Nd (i.e., εNd 80 = +6.54 to +7.48) (Figure 10).This shift indicates a seawater contribution at hydrothermal conditions of upper extrusive oceanic crust, and is a common feature in Mediterranean ophiolite complexes (e.g., Troodos: [83]).Previous authors [83] suggested that interchange of Jurassic and Cretaceous seawater Sr with ocean crust Sr occurs during hydrothermal circulation and may buffer the Sr isotope composition of seawater.Nevertheless, seawater 87 Sr/ 86 Sr has changed with time; thus, before 80 My, seawater was less radiogenic (0.7076) than today and has increased since then [84,85].The value of 0.708 is the maximum value of 87 Sr/ 86 Sr found in the Cycladic metabasites.In this context, the initial 87 Sr/ 86 Sr of the metabasites reflects the mixing ratio of seawater-derived Sr and basalt Sr present in a sample.Thus, the significance of Sr isotope values for determining the characteristics of the magmatic protolith and tectonic setting is limited.From Figure 11a,b it is clear that the Pb-isotopic ratios for Samos and Naxos metabasites show wide variation, from close to MORB into the seawater field.The elevated 207 Pb/ 206 Pb and the large scatter of values, which overlaps with the seawater, also suggest that the Pb incorporated in the Cycladic metabasites is derived from the seawater.The positive correlation between 206 Pb/ 204 Pb and 208 Pb/ 204 Pb for all the Cycladic rocks (serpentinite-metagabbros-metabasalts) is attributed to the mixing relationships between the magmatic and seawater end-member components (hydrothermal system).

Comparisons of the Upper Unit and CBU Ophiolites of the Broader Aegean Region
The geochemical characteristics of the Upper Unit ophiolites reveal significant differences between the Cycladic rocks (Paros-Naxos-Western Samos) and Skyros: (a) Random occurrence of ophiolitic lithologies characterizes the Cycladic mélange units, whereas the Skyros ophiolite consists of more coherent slices of metabasites with essentially boninitic affinities.(b) Basalts with Ocean Island type (OIB) affinities and alkaline characteristics are recognized in the Cycladic islands (Naxos-Ikaria), but are also common elsewhere in the Cycladic belt and the mainland ophiolite complexes: the amphibolite soles of central-northern Evia ophiolite [8] and Ikaria [10], as well as the main Western Ophiolite Belt (Avdella mélange-Pindos and Koziakas: [86]).These rocks display extreme HFSE enrichments relative to N-type MORB (high Ti, P and Nb, Zr/Y ~5, Nb/Y > 1).Alkaline basalts (AB) found in ophiolites mainly include rocks formed at seamounts (OIB) [87].They mainly occur as tectonic slices or blocks incorporated into mélanges, that is, as accreted fragments of oceanic terranes.Contrary to the Pindos peridotites, the Cycladic rocks exhibit variable enrichments in alkali earth elements (Sr, Rb, and Ba).These enrichments are a common feature of the Cycladic serpentinites from mélange units, which have been interpreted in different ways: oceanic hydrothermal alteration, synmetamorphic fluid-rock interaction, igneous differentiation, or some combination of these processes (e.g., S. Evia, Syros: [17,34]).Due to the random variation in these mobile elements (Sr, Rb, and Ba) we consider these enrichments a result of ocean floor hydrothermal alteration.(c) The peridotite relicts in Cycladic ophiolites display considerable chemical and petrographic variation and a general refractory nature.Overall their composition correlates with the two groups distinguished by [74] for the ultramafic rocks of the Pindos complex.The observed elemental variations are typical of serpentinites corresponding to harzburgites from mantle wedge settings (Supra Subduction Zones-SSZ [75]).
There is a scarcity of ages for the formation of metaophiolites from the upper unit of the Cycladic islands (Paros, Naxos, and western Samos) and the broader Aegean region.In most cases, available ages are based on isotopic data for metamorphic soles and on paleontological data of closely associated sediments (e.g., [18] and references therein).Due to that fact, regional correlations between the remnant ophiolitic occurrences exposed in the mélange units of the CBU and the upper unit are difficult to establish.Only for a few ophiolite occurrences ion probe U-Pb zircon ages have been reported [15,76,77].In terms of geochemical characteristics, previous work revealed mid-ocean ridge to island arc affinities (e.g., metagabbros and glaucophanites) from Syros, Tinos, and adjacent islands, and this observation has been interpreted to indicate a back-arc setting [88,89].

Regional Correlation and Significance
A regional correlation of Cycladic metabasites from upper unit ophiolite mélanges and CBU mélanges is attempted (a total of 60 analyses) in order to examine possible clustering of compositions.Published data include the supra-detachment mélange units of Ikaria (Faros and Kefala: [10]), Tinos Upper Unit (Mt Tsiknias: [15]), which can be directly correlated to the low-pressure ophiolitic mélange units of the present study.The metaophiolitic mélanges (within CBU) include S. Andros (Late Jurassic), S. Evia (undated), central Samos-Selcuk nappe, and Tinos CBU (Late Createous) [15,28]).
A combination of discriminant plots have been used to examine the diversity of magma types and refine their tectonic setting.In the Zr vs. Zr/Y discriminant plot (Figure 12a) [60], the Zr/Y ratio is used against Zr as a fractionation index, to investigate long term source heterogeneities and distinguish between basalts erupted in within-plate settings, mid-ocean ridges or volcanic arcs.In Figure 12a the Zr vs. Zr/Y systematics (Figure 12a) suggests that the distribution of geochemical groups is random between the upper unit and CBU.Three magmatic types are recognized: normal MORB, Volcanic Arc Basalts (VAB) and Within Plate Basalts (WPB).Skyros and Tinos (Upper Unit) metagabbros (Figure 12a) define a separate field and plot together with Tethyan ophiolite boninitic lavas, late dykes and the Mariana Fore Arc Basalts fall [57].
A combination of discriminant plots have been used to examine the diversity of magma types and refine their tectonic setting.In the Zr vs. Zr/Y discriminant plot (Figure 12a) [60], the Zr/Y ratio is used against Zr as a fractionation index, to investigate long term source heterogeneities and distinguish between basalts erupted in within-plate settings, mid-ocean ridges or volcanic arcs.In Figure 12a the Zr vs. Zr/Y systematics (Figure 12a) suggests that the distribution of geochemical groups is random between the upper unit and CBU.Three magmatic types are recognized: normal MORB, Volcanic Arc Basalts (VAB) and Within Plate Basalts (WPB).Skyros and Tinos (Upper Unit) metagabbros (Figure 12a) define a separate field and plot together with Tethyan ophiolite boninitic lavas, late dykes and the Mariana Fore Arc Basalts fall [57].
True MORB samples are only a few, the majority of Cycladic metabasites plot in the transitional VAB/MORB field.From the WPB group, only Naxos metabasalts, and two amphibolite samples from Ikaria (Faros) [10], can be geochemically classified as alkali basalt: Nb/Y > 1 and high P, Ba, Sr and Ce N /Yb N , higher than E-MORB (i.e., >2) these geochemical features rather resemble OIB-type (see Table 1).Alkali basalts are very commonly associated with tholeiitic rocks in many Mesozoic Tethyan ophiolite mélanges in the Hellenides, and were mostly interpreted as remnants of intra-oceanic islands (OIB) or sea-mounts in the Pindos ocean (Pindos-Avdella: [74], Koziakas: [86]).Alkali basalts in modern oceanic setting are associated with fractures zones, seamounts, and aseismic ridges [50].
The V-Ti discriminant plot of [56] is considered to provide further information for distinguishing between ophiolites already recognized as having a suprasubduction zone origin [51].Based on Ti vs. V relationships (Figure 12b), the Skyros and Tinos (Upper Unit) metagabbros have very low Ti and exhibit arc-like Ti/V <10.These rocks plot in the same area of late stage boninitic basalts and basaltic andesite dykes of the Pindos-Mirdita ophiolite complexes [57,74].In this diagram (Figure 12b) Paros and Andros and part of Samos metagabbros exhibit arc-like Ti/V = 10-20; however, most of Cycladic rocks plot within the broad field of mid-ocean ridge basalts (MORB) and back-arc basin basalts (BABB): Ti/V = 20-50.

Tectonic Setting
Extensive geochronological dating (U-Pb on zircons) and geochemical studies have focused on metagabbros and felsic gneisses within mélange units [15,17,28,90,91] and the marble-schist sequences of the CBU [76,91].Previous studies have shown, on the basis of δO 18 and εHf values of zircons, that Late Cretaceous igneous zircons (Syros and Tinos: [76,87,90,91]) are characterized by MORB-like affinities, suggesting that the Late Cretaceous magmas were produced by partial melting of depleted MORB mantle.Fu et al. [91] also pointed out the significant crustal component in 160 Ma metagabbros from Andros, implying that some Jurassic gabbroic rocks of the Hellenides are not part of SSZ-type (supra-subduction zone) ophiolites that are common elsewhere along the margin of the Pelagonian zone (Western Ophiolite belt).
So far, both harzburgite compositions and those of metabasites (gabbros and dolerites) are linked by the fact that their chemical affinities share many characteristics with supra-subduction zone setting magmatism [66,92].

Conclusions
The results of this study can be summarized as follows:

Figure 8 .
Figure 8. Spinel peridotite norm compositions of the ultramafic rocks from the Cyclades and Skyros plotted on the ultramafic classification diagram of IUGS.

Figure 8 . 32 Figure 9 .
Figure 8. Spinel peridotite norm compositions of the ultramafic rocks from the Cyclades and Skyros plotted on the ultramafic classification diagram of IUGS.

Figure 9 .
Figure 9. Variation diagrams of (a) SiO 2 ; (b) Al 2 O 3 and (c) CaO vs. Mg# as fractionation index for the Aegean ultramafic rocks; (d) Binary diagram of Mg/SiO 2 vs. Al 2 O 3 /SiO 2 with serpentinized peridotites and serpentinites.Successive magmatic depletion line according to [69-71].The shaded field shows the main magmatic trend defined by Paros samples.Skyros serpentinites display the largest scatter, due to alteration effects.
(a) Metavolcanics (dolerites, basaltic andesites, minor basalts) and metagabbros from all the upper unit ophiolitic exposures exhibit wide ranges of major element composition (SiO 2 : 39-56.6 wt %; Mg-number: 37.5-69.2and high K 2 O i.e., >0.2 wt %), indicating derivation from evolved magmas.(b) The petrographic, major and trace element and Sr-isotope data show that all the studied rocks (metagabbros, metavolcanics, and peridotites) underwent ocean floor alteration.Microscopic observations corroborate that low-temperature alteration during the generation of the ocean crust affected the upper unit metavolcanites and metagabbros.(c) Based on the immobile-element TAS proxy classification diagram, Ti/V, Zr/Y vs. Zr discriminant plots and spidergram patterns of N-MORB-normalized compositions, three magmatic groups are recognized among the upper unit (UU) ophiolitic metabasites: (i) metabasites from Paros and western Samos (Kallithea melange) display affinities of island arc tholeiites (IAT); (ii) Naxos basaltic rocks are highly enriched in incompatible elements and display an OIB type affinity, representing an oceanic island setting; (iii) Metagabbros from the northern Skyros ophiolitic-mélange (Fere-Kambos) and Tinos (Upper Unit-Mt.Tsiknias) exhibit boninitic affinities indicating a fore-arc setting.(d) Central Samos metagabbros of the Cycladic Blueschists Unit (CBU) are primitive rocks with Back-Arc Basin affinities.(e) Most peridotite relicts from Samos, Paros, and Naxos-irrespective of the structural unit-display common chemical affinities similar to ocean floor peridotites formed in a supra-subduction zone environment.Major and trace element characteristics of harzburgite relicts in Cycladic serpentinites and Skyros indicate a highly residual nature of the mantle source, comparable to the mantle wedge serpentinites from oceanic fore-arc settings, in accordance with previous studies.(f) A combination of Pb and Sr isotopic characteristics of the Cycladic metabasites shows the importance of seawater hydrothermal alteration on the composition of the studied metaophiolitic rocks.Trace element patterns of Central Samos metabasites (CBU) and Nd isotopes suggest affinities of Back-Arc Basin Basalts (BABB).

Table 1 .
Whole-rock major and trace element compositions of metabasites and ultramafic rocks from the Aegean ophiolites.

Table 2 .
Rb-Sr and Sm-Nd isotope data of representative ophiolitic metabasites and serpentinites from the Aegean islands of Paros, Naxos, Samos, and Skyros.

Table 2 .
Cont. .d:not detected, bdl: below detection limit.* 80 Ma is the reported age of Samos metagabbros from n