Major and rare earth element mineral chemistry of low-grade assemblages inform dynamics of hydrothermal ocean-floor metamorphism in the Dinaridic Neotethys

Abstract This contribution provides insights into ocean-floor hydrothermal metamorphism of the fast-evolving Dinaridic Neotethys. Mineralogical, geochemical and Sr isotope data collected from altered ophiolites and non-ophiolite basalts/andesites and tuffs of the active continental margin are consistent with hydrothermal alteration trajectories that reflect the host-rock composition. This suggests that hydrothermal fluxes were restricted within a simple closed seawater-fed system. Based on the initial isotopic ratios of Sr, two fluid–rock interaction trends are established: (a) low-to-medium degrees of metasomatism in pre-Middle Jurassic anorogenic ophiolites that progressively abated, and (b) increased intensities of metasomatism in post-Middle Jurassic orogenic ophiolites. This agrees with chlorite thermometry and Ca-Al-(Fe)-silicate phase chemistry. The metamorphic assemblages belong to the zeolite, prehnite-pumpellyite, prehnite-actinolite and greenschist facies. The facies is reliant on the temperature of hydrothermal systems and their fluid chemistry. Rare earth element (REE) phase geochemistry shows (a) variable fluid–rock ratios in chlorite and pumpellyite dependent on fluid temperatures, (b) prominent Eu and Ce anomalies that reflect the fluid oxidation state, (c) light REE/heavy REE mobilization attributed to prevalent ligand complexation, and (d) multi-phase fluid percolation across reaction zones of heterogeneous permeability. This study proposes initiation of simple hydrothermal system(s) at or near a spreading centre(s) in the infancy of the Dinaridic Neotethys. Such a system became more complex during Middle Jurassic and Early Cretaceous time with reactive hydrothermal fluids passing the recharge area and reaching the hot reaction zone. An abrupt obliteration of the established high-temperature regime ensued, following the final closure of the Neotethys.


Introduction
Hydrothermal metamorphism that takes place in the oceanic lithosphere facilitates our understanding of crustal evolution and offers insights into processes that characterize hydrothermal systems (Alt, 2009;Alt & Teagle, 2000). The hydrothermal interaction of seawater and oceanic crust controls the heat exchange between lithosphere and hydrosphere, which ultimately affects the chemistry of seawater as well as the upper sequences of the oceanic crust. The convective flow of fluids in the crust is governed by the heat released during formation and subsequent cooling of various mafic lithologies (Fyfe & Lonsdale, 1981;Alt & Teagle, 2000;Gillis et al. 2005;Wilson et al. 2006). Chemical changes taking place in hydrothermal systems have been reported for a broad range of temperatures, pressures and compositions of catalytic fluids at varying seawater/rock ratios (e.g. Mottl, 1983). Mafic extrusive and intrusive rocks of the oceanic crust react with circulating and percolating hydrothermal fluids (i.e. tempered liquid phases), which results in the transformation of primary magmatic minerals into deuteric assemblages. Alternatively, the latter may precipitate directly from the hydrothermal fluid that strongly reflects the geochemistry of affected crystalline rocks. In this contribution, such hydrothermal phases are referred to as 'metamorphic' since hydrothermal metamorphism commonly includes partial recrystallization and the formation of clay minerals; these are essentially out of thermodynamic equilibrium and result in broad compositional ranges (Arbiol et al. 2021). These hydrothermal processes are best understood as being metasomatic and allochemical, and prone to the selective addition and removal of chemical components (Harlov & Marschall, 2009;Paoli et al. 2019). The hydrothermal metamorphic assemblages may therefore provide information on the alteration setting and post-emplacement alteration processes. Furthermore, over the last few decades research on seawater interaction with mafic and ultramafic rocks has sparked academic interest, which has established a hypothesis on the role of hydrothermal alteration of mid-ocean-ridge basalts in the global cycling of elements (e.g. Seyfried & Bischoff, 1981;Mottl & McConachy, 1990;German & Lin, 2004;Lowell et al. 2013).
In the European Western Vardar ophiolites (in the sense of Schmid et al. 2008; Fig. 1), or more precisely in the southwestern segment of the Zagorje-Mid-Transdanubian Zone (ZMTDZ; Pamić & Tomljenović, 1998) or Sava Unit (Haas et al. 2000), Mesozoic mafic intrusive, extrusive and pyroclastic rocks crop out in northwestern Croatia (Slovenec et al. 2011;Fig. 2). These rocks are incorporated as slices and metre-to-kilometre-sized blocks in ophiolite mélange (Slovenec & Šegvić, 2021) known as the Kalnik Unit (KU) (Haas et al. 2000) (Fig. 2). However, the geotectonic affiliation of the KU is still debated; while some argue that it is part of the Western Vardar ophiolites ( Fig. 1; e.g. Schmid et al. 2008) based on the occurrence of non-ophiolite basalts/andesites and tuffs of the active continental margin (Slovenec & Šegvić, 2021), others suggest that it forms the most southern portion of the Southern Alps (Haas et al. 2000;Schmid et al. 2008; Fig. 1). The former proposal considers the KU to be an integral part of the Jurassic oceanic accretionary prism (e.g. Schmid et al. 2008), while the latter suggests its Middle Triassic continental volcanism to be spatially and genetically related to the onset of the formation of the Tethyan Mesozoic Adriatic-Dinaridic carbonate platform(s) (e.g. Vlahović et al. 2005). Such diverse geological histories likely resulted in various hydrothermal assemblages that are dependent on the host lithology, basement architecture, porosity and permeability of the oceanic crust (Alt, 2009).
In this contribution, we provide an overview of the mineralogy as well as the phase and Sr isotope geochemistry of ocean-floor metamorphic parageneses. Compositional variabilities of hydrothermal assemblages were examined in the context of the geological setting of the host rocks (i.e. igneous versus pyroclastic rocks) where particular attention was paid to the rare earth element (REE) mobility during hydrothermal processes, which has long been a matter of concern because of its bearing on the evolution of the oceanic crust (e.g. Michard et al. 1983;Michard & Albarède, 1986;Bau, 1991;Klinkhammer et al. 1994;Bau & Dulski, 1999). Chlorite and pumpellyite are omnipresent in the oceanic crust preserved in the KU regardless of lithology (Slovenec et al. 2012). Acting as either a clay-size replacement of a pyroclastic/effusive matrix and/or a coarse vein infill, chlorite in particular serves as an excellent medium to investigate the mobility of REEs in hydrothermal media controlled by temperature, redox potential, acidity and ionic complexation of pervasive fluids (Brookins, 1989;Tetiker et al. 2015). The purpose of this investigation is to (a) find mineralogical and geochemical clues regarding how and to what extent hydrothermal parageneses and their phase chemistry are controlled by the primary lithology, and (b) document the REE distribution in hydrothermal, texturally diverse, chlorite and pumpellyite, in order to test a possible relationship between the REE geochemistry of these minerals and the prevailing hydrothermal conditions.

Geological setting
Mesozoic mafic and pyroclastic rocks that crop out at the Intra-Pannonian Inselgebirge of the Sava Unit (Figs 1, 2) are classified as ophiolites and non-ophiolite basalts/andesites and tuffs of the active continental margin. Ophiolites are of Middle-Late Triassic, Jurassic or Cretaceous age and make up an integral part of the subduction-related tectonic mélange of the north Croatian mountains of Kalnik, Ivanščica, Medvednica and Samoborska Gora (Slovenec et al. 2011;Lugović et al. 2015). These constitute members of the Kalnik Unit (Fig. 2), which represents a chaotic sedimentary/tectonic complex composed of pebbles and slivers, and up to hectometre-to-kilometre-sized blocks of sedimentary and igneous rocks embedded in a pervasively sheared continentderived Lower-Middle Jurassic pelitic matrix (Fig. 2a-d). According to these authors, the formation of the mélange took place during Middle-Late Jurassic to Early Cretaceous time. This is the period of rapid transition of the Dinaridic branch of Neotethys, which evolved from an active ridge magmatism to an intraoceanic subduction environment and island/back-arc volcanism (Šegvić et al. 2014). The sedimentary component of the mélange is made up of Triassic and Jurassic sandstone and chert as well as Lower Cretaceous limestone and clastic rocks. Metreto-kilometre-sized blocks of magmatic rocks render a fragmented middle-Lower Triassic and Jurassic oceanic crust that constitute upper members of an ophiolite sequence. Those blocks are largely composed of massive and pillow basalts intersected by mafic dikes (Slovenec et al. 2011). Tholeiitic cumulate gabbro is a minor lithology (Lugović et al. 2015), while alkali basalts and a high-grade metamorphic basement are only locally present (Slovenec et al. 2011;Šegvić et al. 2016). Five geochemically distinct types of mafic extrusives occur in the magmatic blocks of the KU (Slovenec et al. 2011;Lugović et al. 2015): (i) within-plate alkaline basalts (WPAB; Middle Triassic, Illirian); (ii) enriched mid-ocean-ridge basalts (E-MORB; Middle Triassic, late Fassanianearly Longobardian); (iii) transitional mid-ocean-ridge basalts (T-MORB; Late Triassic, middle-late Carnian); (iv) very different varieties of essentially normal mid-ocean-ridge basalts (N-MORB; Late Triassic -Midddle Jurassic, latest Carnian-Bajocian); and (v) island-arc basalts (IAT; Middle Jurassic, late Bajocian-Oxfordian). On the other hand, intrusive rocks occur as (i) N-MORB (Early Jurassic, Pliensbachian-Aalenian); (ii) IAT (Middle-Late Jurassic, Callovian-Tithonian); and (iii) back-arc basin basalts (BABB; Early Cretaceous, Albian). Ophiolite rocks are commonly affected by hydrothermal activity, overprinted by metamorphism up to the greenschist facies (Slovenec et al. 2012).
Triassic non-ophiolite basalts/andesites and tuffs of the active continental margin emerge as interstratified, hectometre-to-kilometre-sized blocks in volcanosedimentary successions of the north Croatian mountains of Ivanščica, Strahinjščica, Kuna Gora, Desinić Gora, Ravna Gora and Žumberak (Fig. 2). These rocks constitute part of the Southern Alps (Schmid et al. 2008; Fig. 1). Volcanic and volcanoclastic rocks are intercalated with Middle Triassic siliciclastic and carbonate marine sediments, and consist of fine-grained and strongly altered basic, intermediate to acid volcanic and pyroclastic lithologies (Marci, 1987;Goričan et al. 2005;Slovenec & Šegvić 2021). Alteration processes documented in nonophiolite basalts/andesites and tuffs are attributed to burial metamorphism and hydrothermal activity, and account for the volcanic glass devitrification, plagioclase degradation and occurrence of low-grade assemblages of prehnite-pumpellyite facies (Marci, 1987). The volcanic activity for that period is represented by submarine andesite to basaltic lava flows accompanied by multiple explosive eruptions of volcaniclastic material. The Middle Triassic volcanosedimentary succession of the Sava Unit can be correlated with the Triassic volcanosedimentary series of the Meliata-Maliak arc-back-arc system (Goričan et al. 2005) and analogue successions from the Dinarides (Pamić, 1984;Smirčić et al. 2018).

3.a. Materials
A set of 110 samples of Mesozoic intrusive (gabbro), extrusive (basalt and andesite) and pyroclastic rocks were selected for mineralogical, chemical and Sr isotopic analyses (Table 1). All rock samples were taken from outcrops in the north Croatian mountains of Kalnik, Ivanščica, Strahinjščica, Kuna Gora, Desinić Gora, Ravna Gora, Žumberak, Medvednica and Samoborska Gora (Fig. 2). The ophiolites and non-ophiolite basalts/andesites and tuffs account for 65 and 45 samples, respectively. All samples were first investigated with a petrographic microscope to define the texture of the rocks, as well as primary and secondary mineral assemblages.

3.b. Electron microprobe analyses
Electron microprobe analyses (EMPA) and the elemental X-ray study were carried out at the Institute of Geosciences (Universität Heidelberg, Germany) using a CAMECA SX51 microanalyser equipped with five wavelength-dispersive spectrometers.
In total, 39 representative samples were analysed. Measurements were performed using an accelerating voltage of 15 kV, beam current of 20 nA, beam size of c. 1 μm (for feldspars, 10 μm) and 10 s counting time for all elements. Natural oxides and silicates were used as standards and for calibration. Raw data were corrected for matrix effects with the PAP algorithm (Pouchou & Pichoir, 1984) provided by CAMECA. Mineral formulas were calculated using the software package MINPET (LR Richard, Gatineau, Québec, Canada).
Two different chlorite geothermometers are utilized in this contribution: the empirical one proposed by Kranidiotis & MacLean (1987) and the semi-empirical one of Inoue et al. (2018). The former considers the amount of tetrahedral Al in the structure of Al-rich Mg-Fe chlorite, which is considered to be temperature dependent (details in Arbiol et al. 2021). The absolute accuracy of empirical geothermometers has been widely challenged (Arbiol et al. 2021); the geothermometer of Kranidiotis & MacLean (1987) has therefore only been used tentatively to show potential differences between the empirical and semi-empirical chlorite thermometers applied in hydrothermal systems. The semiempirical chlorite geothermometer of Inoue et al (2018) is used  Schmid et al. 2020). The upper inset map shows the research area in more detail.
here because it was initially developed for application to hydrothermal systems and assumes the equilibrium of quartz-chlorite-water. Ideally, the Fe 3þ /ΣFe should be measured or estimated, as required for the application of most semi-empirical geothermometers; however, the Inoue et al. (2018) geothermometer permits the calculation of the formation temperature of chlorite assuming ΣFe = Fe 2þ , which is the approach followed in this contribution.

3.c. Scanning electron microscopy
Scanning electron microscopy (SEM) analyses were completed at the Microscopy Center of Texas Tech University using a Zeiss Crossbeam 540 apparatus equipped with an energy dispersive spectrometer (EDS). Carbon-coated polished thin-sections were utilized for this investigation. The measurements were performed at high vacuum, 15 kV and c. 1 nA, with two silicon drift energy dispersive X-ray detectors from Oxford instruments. Highresolution backscatter electron (BSE) images were acquired using a four-quadrant-backscatter detector. The spectra acquisition time was 20 s. Zeiss Aztec software was used for the quantification of EDS spectra in a standardless mode. Chemical data were used as atomic percentages and were normalized to 100%.

3.d. Sr isotope analyses
Strontium isotopic compositions of 38 rocks were measured at the Research Centre for Petrography and Geochemistry (CRPG, Université de Lorraine, France) using a Triton Plus mass spectrometer. Normalizing ratios of 86 Sr/ 88 Sr = 0.1194 were assumed. The 87 Sr/ 86 Sr ratio for the NBS 987 Sr standard for the period of measurement was 0.710242 ± 0.000030 (2σ). Total procedural blanks were c. 500 pg.

3.f. Laser ablationinductively coupled plasmamass spectrometry
The REE content of chlorite and pumpellyite was measured in three polished thin-sections of representative basaltic rocks using an Agilent 7500cs quadrupole mass spectrometer coupled with a New Wave UP-213 solid-state laser with dual-volume cell installed at the Geoanalytical Laboratory of Texas Tech University. For each sample, three to five measurements of phases of interest were undertaken. Individual measurements typically comprised 30-slong background measurements followed by 20-to 40-s-long acquisitions with the activated laser. The choice of measurement spot locations was facilitated using a BSE image from the SEM. The laser was operated using a frequency of 15 Hz, a spot size of 40 μm and a measured fluence of 6-7 J cm −2 . For every measurement, corresponding time-resolved charts displaying countsper-second (cps) data were examined for potential mineral contaminations. The GSD-1 Universal Geotechnical Standard and Si abundances from EDS microanalyses were used as external and internal standards, respectively. The United States Geological Survey rock standard BHVO-2 served to monitor instrument performance, precision and accuracy (Jochum et al. 2005).

4.b. XRD mineralogy and EMPA chemistry of hydrothermal phases
The distribution and relative abundances of hydrothermal and weathering secondary phases in analysed ophiolite and nonophiolite basalts/andesites and tuffs are provided in Table 2, while representative phase chemistry and calculated mineral formulas are given in online Supplementary Tables S1a-e (available at http://journals.cambridge.org/geo). Alteration minerals comprise about two-thirds of the studied rocks. The most abundant alteration products are albite, chlorite and quartz. Minor hydrothermal phases are amphibole, pumpellyite, zeolite, prehnite, epidote/ clinozoisite, titanite, rutile, illite/mica, kaolinite, calcite, Fe-oxides, sulphides, carbonates, serpentine-phases, sepiolite/palygorskite, stilpnomelane and various expandable clay minerals. In comparison to ophiolite rocks, the non-ophiolite basalts/andesites and tuffs show notable differences in their hydrothermal mineralogy, which contain more quartz, calcite, illite/mica, illite-smectite and pumpellyite. Conversely, these rocks seem to be comparatively impoverished in titanite, rutile, ilmenite and epidote.

4.b.3. Feldspars
In ophiolites, hydrothermal albite (Table 2) is largely a product of a complete alteration of plagioclase  Ab 82.9-26.1 Or 1.6-0.1 ) (Figs 4, 6), while in non-ophiolite andesite and basalts it grows on a substrate consisting of K-feldspar (Fig. 3h). In albite, the contents of SiO 2 , Al 2 O 3 and FeO vary within a narrow range with the An content spanning 0.69-9.52 (online Supplementary Table S1c; Fig. 7b). Fine-grained muscovite as well as the mineral aggregate saussurite are commonly associated with albite in the form of hydrothermal alteration products of igneous plagioclase. The latter is especially common in Jurassic rocks rich in chlorite.

4.b.4. Chlorite
Chlorite is a common hydrothermal alteration phase in all rock types, and is the most abundant secondary mineral after albite (Table 2; 4,6). Chlorite is characterized by a narrow range of SiO 2 and Al 2 O 3 contents, which are the lowest in the youngest (Lower Cretaceous) BABB gabbroic rocks (online Supplementary Table S1b). As per octahedral layer cation occupancy (5.49-6.33 apfu), chlorite is tri-octahedral (Type 1, Fig. 7c). The nature of dominant divalent cation (i.e. clinochlore versus chamosite) was found to reflect the geotectonic affiliation of host rocks; namely, the ratio of Fe 2þ /Mg steadily increased starting with the WPAB ( (Fig. 7c, d).
Radial to fibrous chlorite commonly occurs in the interstices of pyroxene and plagioclase (Fig. 3b, e, f). In ophiolitic rocks, chlorite is accompanied by epidote, prehnite, calcite and neoalbite (Fig. 3g). Where associated with pumpellyite, it commonly occurs on the peripheral parts of basaltic pillows. In non-ophiolite basalts/andesites and tuffs, however, fine-layered chlorite is a product of devitrification of weathered volcanic glass (Fig. 3h). Chlorite in those rocks commonly forms a substrate of the hypocrystalline matrix or fills monomineralic amygdules (online Supplementary  Fig. S1a, b); alternatively, it may be associated with vein carbonates (online Supplementary Table S1b). The matrix is a mixture of c. 10-μm-sized Mg-chlorite and altered feldspar. Numerous, 0.1-2-mm-long elliptical to circular amygdules filled with Mg-chlorite dominate the matrix (online Supplementary  Fig. S1b). The contact of amygdaloidal chlorite and rock matrix is sharp, but with a visible 1-5-μm-thick reaction front (online

4.b.5. Prehnite
Prehnite occurs in ophiolites as well as in most of the Middle Triassic non-ophiolite basalts (Table 2). However, its abundance is rather modest except in cumulate gabbro and pillow lavas where it is commonly associated with pumpellyite (Figs 3e, f, 6c-e). Prehnite occurs as needle-to plate-like aggregates (Fig. 3e). The  Table S1d) with FeO showing a broad compositional range (0.33-5.30 wt%). Prehnite occurs not only as an aggregate between plagioclase and amphibole, but it also fills millimetre-sized veins (Fig. 3g).

4.b.6. Pumpellyite
Pumpellyite is present in all varieties of extrusive and some intrusive rocks (online Supplementary Table S1); however, like prehnite, the amount of pumpellyite is significant only in pillow lavas that occur in hectometre-to-kilometre-sized fragments of oceanic lithosphere in the ophiolite mélange. Pumpellyite occurs as a bluish to greenish needle-like to fine-grained aggregate (Fig. 3b, e, f) commonly associated with prehnite and chlorite (Figs 4c, 6c, Fig. 8a). The albitization of plagioclase delivers the necessary amount of Al for the crystallization of pumpellyite (Fig. 6d). With a somewhat lower Al/Fe ratio (online Supplementary Table S1d), pumpellyite was likely derived from the alteration of clinopyroxene (Figs 4c, 6f) or alternatively formed through the reaction of Fe(III)-rich prehnite and chlorite, giving rise to pumpellyite and quartz (Deer et al. 2013).

4.b.7. Zeolite
Zeolite is only present in Middle Triassic E-MORB and Middle Jurassic N-MORB ophiolite extrusive rocks (Table 2). Zeolite in basalts is fine-grained and whitish in colour, with a platy to radial habit characteristic of laumontite (Fig. 6b). Zeolite occurs as monomineralic 1.5-mm-wide veins or albitized domains of primary plagioclase. Ca-zeolite is characterized by a low content of total alkalis (Na 2 Oþ K 2 O < 0.42 wt%; online Supplementary Table S1d; Fig. 8c). As per the Si/Al ratio (2.05-2.13), laumontite corresponds to the intermediate-silica type of zeolite (Campbell et al. 2016).

4.b.8. Other accessory phases
Titanite is an accessory phase in all ophiolitic intrusive and extrusive rocks, while in non-ophiolite rocks it is rare (Table 2). It may be found as fine-grained sub-to anhedral aggregates (up to 0.5 mm in size) associated with ilmenite and magnetite (Figs 3f, h, 4, 6). Titanite is commonly altered and transformed into dense,  (Girardeau & Mevel, 1982).
cryptocrystalline aggregates of leucoxene within plagioclase and chlorite (Fig. 3f). The composition of titanite defines it as grothite (Deer et al. 2013), which is a variety rich in Al, Fe and structural vacancies (Fig. 8b). The latter is reflected in low sums of total cations per formula unit (online Supplementary Table S1e).
Since leucoxene commonly accompanies titanite, it is reasonable to hypothesize that leucoxene represents accumulations of cryptocrystalline titanite and Ti-oxides. Anatase has been identified by XRD in several samples of N-MORB, IAT and BABB Jurassic and Cretaceous basaltic rocks (Table 2). Pyrite and pyrrhotite (online Supplementary Table S1e) occur in ophiolite and non-ophiolite lithologies (Table 1) as c. 50-μm-sized particles (Fig. 6b). Haematite and Ti-haematite, comparatively more rare than the sulphides, occur in N-MORB Upper Triassic -Middle Jurassic basalts (online Supplementary Table S1e; Fig. 6f).

4.c. Sr isotopic composition
The Sr isotopic composition of 38 Mesozoic mafic intrusives, extrusives and pyroclastic rocks is shown in Table 3. The initial isotopic ratios of Sr are calculated for the corresponding age of each sample (i.e. the assumed age of crystallization; 245-103 Ma). This age span is based on Ar/Ar and K/Ar isotopic ages of Lugović et al.

4.d. REE geochemistry of chlorite and pumpellyite
Chlorite and, to a certain degree, pumpellyite occur in most of the analysed samples (Table 2). These phases were used to investigate the mobility of REE during alteration processes that postdate the formation of igneous host rocks. Three samples of non-ophiolite basalts were deemed representative because of the various textures involving chlorite and pumpellyite. Textural differences are likely the result of a multitude of nascent hydrothermal events, which show contrasting REE signatures. The REE geochemistry of chlorite and pumpellyite are given in online Supplementary Table S2. REE concentrations of chlorite (online Supplementary Fig. S1a-c) display consistent values where the standard deviation is on average less than 15% for each element (online Supplementary  Table S2). The REE content of chlorite is low (45-50 ppm), ranging over c. 0.45-3.17 relative to chondrite (Boynton, 1984;Fig. 10). The Eu anomaly is moderate (Eu chon /Eu chon * = 0.73-1.11; Eu chon * = (Sm chon þGd chon ) × 0.5) and the REE normalization curves are somewhat similar to N-MORB with a substantial depletion in LREE (La chon /Yb chon = 0.26-0.33).

5.a. Mineralogical changes and their implications
Taking into account all ophiolites and non-ophiolitic volcanic lithologies, the intensity of hydrothermal alteration seems to be strongest in the ophiolite extrusive rocks ( Table 2). The alteration intensity of the rocks, which can be equated to an increase in the seawater-rock interaction rate, may be monitored through the oscillation of the initial Sr isotopic ratios. This feature serves as an alteration index (Bickle & Teagle, 1992) of studied Mesozoic lithologies (Fig. 9). This is recognizable in Mesozoic mafic intrusive and extrusive rocks; namely, the oldest ophiolite rocks are part of anorogenic domains (WPAB, E-MORB, T-MORB, N-MORB) that are related to the extension of the oceanic lithosphere and the formation of the western branch of the Neotethys (Slovenec et al. 2011; Fig. 11a,b). These rocks are weakly to moderately altered with a dwindling Middle Triassic -Middle Jurassic Sr alteration trend (Fig. 9). This seafloor spreading system was succeeded by Neotethyan intraoceanic subduction during the Middle Jurassic Epoch (Bathonian; Fig. 11c). Related suprasubduction processes during the Late Jurassic Epoch resulted in the development of orogenic domains (IAT, BABB; Slovenec et al. 2011). A distinctly increasing Sr alteration trend is documented in these rocks, delineating a growing impact of seawater that lasted through latest Jurassic and Early Cretaceous time (Fig. 11d). Described trends are well correlatable with changes in seawater composition during the Mesozoic Era ( Fig. 9; Kovács et al. 2020). In contrast to ophiolite rocks, Middle Triassic non-ophiolite basalts and tuffs (CAB; Fig. 11a) from the active continental margin indicate a broad range of alteration intensities devoid of any obvious trend (Fig. 9). Another useful parameter that may provide insights into the alteration dynamics is the Sr/Ca ratio (Scambelluri et al. 2001). Its relatively high values (25-74) in Middle Triassic -Middle Jurassic ophiolite mafic intrusives and extrusives rocks that crystallized in anorogenic domains (WPAB, E-MORB, T-MORB, N-MORB), compared with somewhat lower Sr/Ca values  in Upper Jurassic to Lower Cretaceous suprasubduction series rocks, suggest a gradual decrease in the degree of albitization and plagioclase alteration, which is in line with the earlier findings of Slovenec et al. (2011). This trend may further outline a progressive dwindling in depths at which hydrothermal processes occurred, thus suggesting an ocean-floor transition from the greenschist to zeolite facies metamorphism during Middle Triassic -Early Cretaceous time.

Dinaridic Neotethys hydrothermal ocean-floor metamorphism
Bathonian-Valangian highly altered ophiolite rocks (IAT and BABB; Fig. 11c,d) are practically devoid of sulphides, while oxide minerals are commonly present. Sulphide precipitation must have been triggered by a temperature decrease of the mineralizing fluid, coupled with decreased sulphide solubility below 300°C (Henley et al. 1984). Cooling is normally caused by a drop in pressure associated with fluid ascending in the hydrothermal system or cold seawater mixing (Capuano & Cole, 1982). This in turn makes the pH of hydrothermal fluid decrease following sulphide precipitation (Akaku et al. 1991). The moderately to highly reducing character of sulphide-producing fluids and reduced pH causes redox-sensitive phases such as Fe oxides to dissolve (Henley et al. 1984). This likely explains the lack of haematite and magnetite in sulphide-bearing samples ( Table 2). Sulphide-producing hydrothermal fluids stem from deep-seated portions of an active ocean ridge while parageneses free of sulphides are related to the relatively oxidizing environment prevailing at shallow levels of the hydrothermal system (Fouquet et al. 2018). Because the latter pertains to Bathonian-Valangian orogenic ophiolite rocks, it is hypothesized that an intensive suprasubduction tectonism facilitated access to the hot recharge zone of the hydrothermal system without changing the redox character of the mineralizing fluid. This is consistent with Sr/Ca values in Upper Jurassic -Lower Cretaceous rocks, which called for a progressive shrinkage of depths associated with seafloor metamorphism.

5.a.1. Basaltic ophiolitic and non-ophiolitic rocks
The intrinsic textural and structural characteristics of basaltic rocks were largely preserved during hydrothermal metamorphism, despite significant mineralogical changes. Submarine effusive eruptions resulted in an instantaneous quenching of lava to form pillow lavas (e.g. Fujibayashi et al. 2014). As the new seafloor was produced at mid-ocean ridges during Ladinian-Bajocian time (Fig. 11b), the high permeability of the newly formed oceanic crust allowed ample quantities of water to penetrate and acquire heat from the crust (Kuhn et al. 2017). At high temperatures, hydrothermal upwelling fluids eventually reacted with the oceanic crust to form new phases or recrystallized existing phases, which accommodated ascending pathways of venting fluids. Hydrothermal mineral zones contain typical parageneses of very-low-(200-360°C) to low-grade (360-500°C) metamorphic facies (Slovenec et al. 2012).
In the hydrothermal seafloor setting, albitization of plagioclase (online Supplementary Table S1c) commonly takes place in zones of high seawater activity at moderate temperatures and/or upwelling zones where Si-Na-rich and Ca-poor fluids, which stem from high-temperature zones at depth, are cooled down during their ascent (Von Damm et al. 1985;Ray et al. 2009). The albitization-related release of Ca facilitated crystallization of hydrothermal Ca-bearing phases such as calcite, epidote, Al-pumpellyite, Ca-zeolite, zoisite/clinozoisite and titanite. The high amounts of chlorite in primary plagioclase (Fig. 4f) indicate that the albitization of plagioclase and chloritization were penecontemporaneous and mutually dependent processes. Based on chlorite geothermometry (Inoue et al. 2018), temperatures ranged over 92-489°C (interstitial chlorite, online Supplementary Fig. S2; Fig. 7d). The phase chemistry of the hydrothermal chlorite reflects the geotectonic affinity of host rocks that can be traced during a span of c. 150 Ma from the formation of the initial oceanic lithosphere during Middle Triassic time, through the spreading of Neotethys which lasted until Middle Jurassic time and, lastly, until Late Jurassic -Early Cretaceous intraoceanic subduction Table 3.  (Table 2). Reconciling the chlorite phase chemistry and its crystallization temperatures, an inference can be made about high-temperature post-magmatic hydrothermal conditions (prehnite-actinolite facies; Mg no. chl > 57; online Supplementary Table S1b) that existed during Middle-Late Jurassic time during an intensive intraoceanic subduction in the western segment of the Neotethys (Fig. 11c, d). This resulted in a great diversity of hydrothermal parageneses in the rocks of the ophiolite mélange of the Sava Unit (Slovenec et al. 2011) with abundant crystallization of prehnite (online Supplementary Table S1d) in the range 190-340°C (Liou et al. 1985), which implies an origin at an early stage in the alteration history of the Sava Unit oceanic lithosphere. A somewhat lower hydrothermal temperature then prevailed (180-250°C, i.e. prehnite-pumpellyite facies; Mg no. chl < 57; online Supplementary Table S1b) during Late Jurassic -Early Cretaceous time in a newly established back-arc extensional regime (Fig. 11d). Under these conditions in effusive ophiolite rocks (dominantly in pillow lavas) along with chlorite, the dominant alteration phase is pumpellyite (online Supplementary  Table S1d). As per the criteria of Liou et al. (1985) the studied Al-poor pumpellyite (Al ≤ 2.290 apfu) formed at temperatures ranging over 160-320°C, thus corroborating prehnite-pumpellyite facies conditions (Fig. 8a). The pressure needed for the crystallization of prehnite-pumpellyite-bearing assemblages is 0.1-0.2 GPa (Schiffman & Liou, 1983). Conversely, less common Al-pumpellyite (Al > 2.290 apfu) with a high ratio of Al/Fe tot (4.14-6.90) present in zeolite-free parageneses likely formed at a temperature not significantly lower than 250°C and at pressures of 0.3-0.5 GPa.   (Boynton, 1984) for hydrothermal chlorite and pumpellyite as well as the relict clinopyroxene and whole-rock basalt rocks from the study area. Dashed lines represent normalization plots of black-smoker solutions (East Pacific Rise, Michard et al. 1983) and seawater (Høgdahl et al. 1968).
This supports the idea of prehnite-actinolite facies conditions at comparatively higher alteration depths (Liou et al. 1985). The variations in the Al content of pumpellyite are largely the result of the Al-Fe(III) substitution that indicates that, in addition to pressure and temperature, fO 2 exercises a further control on pumpellyite stability (Inoue & Utada, 1991). Indeed, the formation of Alpumpellyite occurs at lower fO 2 (Schiffman & Liou, 1983), which is consistent with a proposed crystallization at high depths. It is therefore reasonable to infer that during the Mesozoic Era (Middle Triassic WPAB and E-MORB through Upper Jurassic suprasubduction extrusive rocks) pumpellyite became progressively impoverished in Al ( Fig. 8a; online Supplementary  Table S1d). This indicates a gradual decrease in crystallization temperatures at lower depths of hydrothermal activity. In a near-surface environment in contact with hydrothermal fluids, during Triassic-Jurassic ocean crust spreading as well as Upper Jurassic subduction, titanitization must have taken place under conditions of prehnite-pumpellyite facies metamorphism (Fig. 8b). This is corroborated by the titanite chemistry (online Supplementary Table S1e) from effusive rocks that contain relatively high abundances of Al (0.12-0.18 apfu) and Fe (0.05-0.08 apfu) as well as low sums of total cations. Al and Fe commonly substitute for Ti, which takes place at temperatures ranging over 180-320°C as suggested by Enami et al. (1993) who revealed that, in metamorphic systems, the average Al þ Fe 3þ content of titanite increases systematically with decreasing temperatures.
In a continuous sequence of ocean-floor metamorphism characterized by conditions of low-pressure (< 1 GPa) and temperature (≤ 240°C) (Liou et al. 1985) paragenesis, ophiolite effusive rocks of the investigated area are sporadically zeolitized to produce laumontite (Table 2; Fig. 8c). This process was detected during the Middle Triassic -Middle Jurassic oceanization (Fig. 11b). The origin of this Ca-zeolite is related to the fluid-induced re-equilibration and the breakdown of Ca-rich plagioclase to secondary phases (Deer et al. 2013). Considering that the lowest temperature zones of thermal metamorphic aureoles in volcanic terranes typically contain laumontite, it may be hypothesized that its formation temperature was c. 125°C. This correlates well with a crystallization temperature of diagenetic laumontite from deep-buried clastic reservoirs that grew on the volcanoclastic substrate (e.g. Ilijima & Utada, 1972).

5.a.2. Gabbroic ophiolitic rocks
In gabbroic ophiolite rocks, however, which are thoroughly affected by the seafloor hydrothermal activity, a multi-phase alteration sequence has been identified with distinct alteration assemblages. The first phase records deuteric changes at high subsolidus temperatures where clinopyroxene (augite) is replaced by hornblende (Fig. 3b). The formation of deuteric magnesiohornblende (Figs 4a,b,d,f,5a) represents the late magmatic crystallizate formed at temperatures not exceeding 500°C. This event marked the onset of the ocean-floor metamorphism (Coleman, 1977), and is related to the period of Triassic-Jurassic oceanization, Upper Jurassic subduction and Lower Cretaceous back-arc extension ( Fig. 11b-d). In some samples of gabbro, a second reaction rim consisting of anthophyllite occurs around hornblende (online Supplementary Table S1a; Figs 4d, 5b). It follows that the altered gabbro made of deuteric hornblende þ tremolite/actinolite þ anthophyllite þ albite þ chlorite ± titanite formed under pumpellyite-actinolite to middle greenschist facies conditions (200-450°C, 0.2-0.4 GPa;Liou et al. 1985). These conditions stem from intensive tectonic activity and crust emplacement in the Neotethyan arc and back-arc settings, which formed during Late Jurassic and Early Cretaceous time (Slovenec et al. 2011;Lugović et al. 2015;Fig. 11d). The next alteration stage is typical of the greenschist facies and includes the development of albite, tremolite, epidote, chlorite and serpentine that grew on a primary substrate of plagioclase, clinopyroxene and/or tschermakitemagnesiohornblende and olivine (Figs 3b, 4a-f; online Supplementary Tables S1a-S1c). Interstitial chlorite formation temperatures largely of c. 250°C (online Supplementary Fig. S2) are in line with such inferences. This stage is related to Upper Jurassic subduction and Lower Cretaceous back-arc extension (Fig. 11c, d). Epidote-bearing parageneses containing albite, chlorite, prehnite ± tremolite (Fig. 3d, e) and their microtexture, as well as substitution values of Fe 3þ in the octahedral site (i.e. pistacite component) of epidote (Fe 3þ /(Fe 3þ þAl) = 0.29-0.31) are all indicative of a relatively high crystallization temperature (320-420°C) irrespective of the host rock and low rock to seawater ratios (Liou et al. 1985). This correlates well with the experimental results of Berndt et al. (1989), and conforms with low-pressure conditions ranging over 0.1-0.2 GPa. Accordingly, the process of epidotization belongs to the transitional zone between the prehnite-actinolite and lower greenschist facies. Epidote and associated chlorite are therefore suggested to be products of reactive components derived from plagioclase and clinopyroxene on the one hand, and high-temperature hydrothermal fluids on the other (Deer et al. 2013). This phase (clinopyroxene) further includes the re-equilibration of Ti-rich spinel through the exsolution of minute pseudobrookite lamellae. The final alteration step produced prehnite, pumpellyite, titanite and illite/mica (Fig. 4c, f) during Upper Jurassic -Lower Cretaceous back-arc extension (Fig. 11d).
Pumpellyite is either Fe-or Al-rich, meaning it was formed after clinopyroxene or plagioclase, respectively (Ishizuka, 1991). The formation of pumpellyite, irrespective of its chemistry, denotes the transition between the zeolite and greenschist facies (e.g. Kamimura et al. 2012).

5.a.3. Controls on the ocean-floor hydrothermal mineralogy
While the fluid temperature, pressure, composition and oxygen fugacity exercise an important role in the evolution of hydrothermal assemblages, it is the chemistry of magmatic host lithologies that sets the trajectories of mineralogical alterations documented here. The difference in chemistry between the two major mafic rock suites (i.e. ophiolite and non-ophiolite basalts/andesites and tuffs of the active continental margin; Table 2) is evident. The SiO 2 content for ophiolites ranges over c. 46-50 wt%, TiO 2 is c. 2 wt%, while CaO and K 2 O contain up to 11 and 1.5 wt%, respectively (e.g. Slovenec et al. 2011;Slovenec & Šegvić, 2018). Conversely, the composition of non-ophiolitic basalts/andesites and tuffs (e.g. Goričan et al. 2005;Slovenec & Šegvić, 2021) is marked by higher levels of SiO 2 (c. 53-57 wt%) and K 2 O (c. 3-6 wt%), whereas CaO and TiO 2 do not exceed c. 5.5 and 1 wt%, respectively. Compositional differences are reflected in the mineralogy of the hydrothermal parageneses of non-ophiolite basalts/andesites and tuffs, which are dominated by quartz and K-bearing phyllosilicates, while ophiolite rocks are virtually free of quartz and mica, and are more enriched in TiO 2 -phases and Ca-silicates such as prehnite, pumpellyite and laumontite ( Table 2). This is consistent with a closed architecture of the hydrothermal systems of the Dinaridic Neotethys where fluid fluxes were likely confined within a single lithological unit, yielding relatively uniform but contrasting hydrothermal assemblages in ophiolite on the one hand, and non-ophiolite basalts/andesites and tuffs on the other. A plausible explanation may be the lack of deep normal faulting that channelizes fluid flows to deeper crustal levels (e.g. Schiffman & Staudigel, 1994), possibly as a result of rapid transitions in the Dinaridic Neotethys geotectonic setting, changing from active ridge magmatism to an intraoceanic subduction environment and island-arc volcanism (Slovenec et al. 2011;Šegvić et al. 2014). Finally, petrographic and EMPA/SEM observations yielded no evidence of alteration of hydrothermal phases formed in ophiolite and non-ophiolite basalts/andesites and tuffs, which is consistent with the closed hydrothermal system hypothesis advocated above (e.g. Morad et al. 2010). Furthermore, our XRD investigation shows that illite-smectite (or smectite) commonly occurs in non-ophiolitic rocks, whereas Mg-phyllosilicates (chlorite-smectite, vermiculite, lizardite, palygorskite/sepiolite) dominate ophiolite rocks ( Table 2). Because of their submicron size, both illite-smectite and Mg-phyllosilicates formed diagenetically through the weathering of mica and chlorite/amphibole, respectively. This took place following the emplacement of these rocks, which is consistent with a recent investigation of the ophiolite mélange of the study area that, based on the Kübler and Árkai indices, indicates sheet silicates of the mélange are solely diagenetic in origin (Judik et al. 2005).

5.b. Chlorite and pumpellyite REE geochemistry
The normalized REE content of seawater (Høgdahl et al. 1968) and that of modern hydrothermal fluids (black smoker) (Michard et al. 1983;Michard & Albarède, 1986) with reference to REE abundances in analysed minerals is shown in Figure 10. Relatively similar REE levels are inferred for the amygdaloidal chlorite mineralization (online Supplementary Fig. S1a, b) and seawater with the exception of MREE where abundances are somewhat higher in the former (Fig. 10). Conversely, the REE patterns of Fe-chlorite and pumpellyite phenocrysts (online Supplementary Fig. S1c, d) are similar to those of black smokers; the REE concentrations are, however, higher in chlorite/pumpellyite by several orders of magnitude (Fig. 10). An inference can be made that in Fe-chlorite and pumpellyite virtually all of the REEs stem from the host rock, while for Mg-rich chlorite amygdules this is not necessarily the case. Research on hydrothermal fluids collected from various geothermal fields shows that the REE balance of percolating fluids increases at higher temperatures (> 230°C) and when pH decreases (Michard et al. 1983), while in low-temperature solutions, the REE content remains low and will not alter regardless of any reasonable fluid/rock ratios (Sanjuan et al. 1988). This means that Mg-chlorite from basaltic amygdules originated from the cold seawater that penetrated down in the recharge zone of the hydrothermal system but did not reach the hot reaction zone (e.g. Valsami & Cann, 1992) or, alternatively, the hydrothermal system was temporarily inactive or cooled down. Such inference is corroborated by the formation temperatures of chlorite (72-100°C), which fills the amygdules and microveins in rock matrix (vein chlorite, online Supplementary Fig. S2). In the areas of active crustal tectonism, repetitive changes from the active to inactive hydrothermal regime are typical (Moeck, 2014), while periods of activity may be separated by thousands of years of quiescence (Heasler et al. 2009). The LREE values shown by amygdaloidal chlorite (online Supplementary Fig. S1a, b) represent a further line of evidence of its low-temperature origin. This illation ensues from a consensus on high-temperature fluids collected from deep-seated hydrothermal areas that are LREE enriched (Klinkhammer et al. 1994;Şener et al. 2017). Also, the analogue LREE enrichment is documented in hydrothermal chlorite formed from such high-temperature solutions (Tetiker et al. 2015 , while the latter is commonly complexed by CO 3 2− (Valsami & Cann, 1992;Li et al. 2019;Gong et al. 2021). The abundance of Upper Palaeozoic and Triassic carbonates interlayered with the studied extrusive rocks and tuffs (Goričan et al. 2005;Slovenec & Šegvić, 2021) must have contributed to the saturation of the ocean-floor seawater by carbonate ions. These solutions are suggested to have entered ancient hydrothermal systems, interacted with host rocks in the recharge zone, and shortly thereafter produced chlorite amygdules marked by elevated HREE/LREE ratios. This likely explains the characteristic N-MORB-like REE patterns of this type of chlorite (Fig. 10).
REE concentrations of Fe-chlorite (online Supplementary  Fig. S1c) and, to a lesser extent, pumpellyite (online Supplementary Fig. S1d) are comparable to those of REE in their extrusive host rocks, which is contrary to the trend exhibited by amygdaloidal chlorite (Fig. 10; D. Smirčić, pers. comm., 2019). SEM images of the former offers an explanation regarding the behaviour of REEs in primary plagioclase (online Supplementary  Fig. S1d) and clinopyroxene (online Supplementary Fig. S1e, f), which commonly occur as inclusions in Fe-chlorite. Eroded grain boundaries of mineral inclusions in chlorite are a sign of fluidfacilitated fast chemical reactions that are not in equilibrium with the fluid (e.g. Jamtveit et al. 2016). Triassic extrusive rocks of the Western Neotethys affected by hydrothermal activity are largely composed of plagioclase, clinopyroxene, K-feldspar and multiple minor phases (Marci, 1987;Slovenec & Šegvić, 2021). With the exception of plagioclase and clinopyroxene, none of the other phases show any significant REE partitioning (Rollinson & Pease, 2021). Moreover, compared with clinopyroxene, plagioclase is two to three times more abundant in the studied rocks (Marci, 1987) and dissolves preferentially in hydrothermal fluids (Slovenec & Šegvić, 2021). It then follows that in these rocks plagioclase is the main carrier of REEs, and therefore effectively controls the REE content of deep-seated, high-temperature hydrothermal fluids (e.g. Klinkhammer et al. 1994;Bau & Dulski, 1999). Several studies have demonstrated remarkable similarities in REE abundances between hydrothermal fluids and plagioclase phenocrysts (Schnetzler & Philpotts, 1970;Papike et al. 1996) characterized by the LREE enrichment and a slightly negative Eu anomaly. Both are marked by REE patterns of host extrusive rocks and Fe-chlorite phenocrysts (Fig. 10). A depletion in Yb and Lu in Fe-chlorite is attributed to the strong complexation of REEs by Cl 2− in deep-seated acidic hydrothermal fluids, whereby HREEs tend to be less fractionated in hydrothermal precipitates (Flynn & Wayne Burnham, 1978;Allen & Seyfried, 2005). Gadolinium, in contrast, prefers the solid phase (Fig. 10), which corroborates the acidic nature of the percolating fluids (Zielinski & Frey, 1974). An extremely prominent negative anomaly of Ce shown by Fe-chlorite (Fig. 10) reflects the composition of its parental fluid that was sufficiently oxidizing for immobile Ce 4þ to form, which in turn led to the impoverishment of Ce in the fluid relative to other REE (Sanematsu et al. 2011). According to the mechanism proposed by Valsami & Cann (1992), ocean water that enters the hydrothermal recharge zone has to move rapidly and still be oxidizing when attaining high-temperature parts of the system. Since a negative Ce anomaly does not occur in high-temperature hydrothermal pumpellyite, a heterogeneous permeability in the recharge and reaction zone is hypothesized, which gives rise to slow-percolating fluids that are readily reduced prior to the reaction with the rock in which a Ce anomaly is absent (Fig. 10). This line of reasoning is corroborated by the random distribution of sulphides in non-ophiolite basalts/andesites and tuffs (Table 2); specifically, their presence calls for the slower movement of mineralizing fluid, which is followed by temperature decrease leading to dissolution of Fe oxides and eventually sulphide precipitation (e.g. Henley et al. 1984).
The REE distribution in phenocrysts of Fe-chlorite is controlled by the temperature, pH, oxidation state, REE complexation by the ligands in the fluid and the fluid/mineral partitioning. Chlorite crystallized from high-temperature solutions will accommodate most of the REE budget of the fluid as a result of the elevated REE distribution coefficients and a relative dominance of chlorite in hydrothermal alteration assemblages (Table 2; Xiao & Chen, 2020;Tan et al. 2021).
The REE patterns of pumpellyite phenocrysts are analogous to those of Fe-chlorite (Fig. 10), and therefore a similar origin through a high-temperature hydrothermal dissolution of plagioclase and, to a lesser degree, clinopyroxene is proposed. The absence of a Ce anomaly and depleted values of HREEs are attributed to the heterogeneity of the deep-seated hydrothermal recharge zone where fluids percolated at a different speed, which likely constrained the redox potential. Indeed, the positive Eu anomaly in pumpellyite is characteristic of acidic high-temperature reducing fluids where trivalent Eu is reduced to Eu 2þ and thus fractionated from the other REEs (e.g. Bau, 1991;Bau & Dulski, 1999).
In summary, the REE distribution in chlorite and pumpellyite of hydrothermal origin can be used as an indicator of the extent and conditions of alteration in fossil hydrothermal systems. We propose a range of hydrothermal solutions that produced amygdaloidal chlorite and phenocrysts of Fe-chlorite and pumpellyite. These phases formed under contrasting temperature conditions and from different solutions, meaning that very distinctive REE patterns resulted (Fig. 10). In the case of chlorite and pumpellyite phenocrysts, minor differences in their REE abundances (Fig. 10)
are attributed to the variations in the redox state and ligand complexation of fluids from which these phenocrysts were derived.

5.c. Geodynamics of the Mesozoic ocean-floor metamorphism in the northwestern segment of Neotethys
The active northwestern continental margin of Palaeotethys, which is composed of non-ophiolitic basaltic (CAB) to acidic rocks and pyroclastites (Fig. 11a), was affected by hydrothermal ocean-floor metamorphism during Middle Triassic time (Anisian-Ladinian; e.g. Lustrino et al. 2019). This led to large-scale albitization and chloritization and, to a smaller extent, prehnitization and titanitization of precursor rocks. Geodynamic processes of initial rifting of the northwestern segment of Neotethys, which commenced during late Anisian and early Ladinian time, resulted in outflows of primitive alkaline WPAB and OIB lavas ( Fig. 11a; Slovenec et al. 2011).
Intensive ocean dynamics continued during Late Triassic and Middle Jurassic (Bajocian) time, having transitioned to rapid ocean spreading and formation of E-, T-and N-MORB oceanic lithosphere ( Fig. 11b; Slovenec et al. 2011). This was subjected to hydrothermal fluxes and seawater ( Fig. 9) and gave rise to abundant crystallization of hydrothermal phases under prehnite-actinolite, prehnite-pumpellyite and zeolite facies conditions. Albite, chlorite, prehnite, pumpellyite and titanite are the most abundant hydrothermal phases with minor magnesiohornblende, tremolite, haematite and rare laumontite. The Middle Jurassic vergence transition from extensional to compressional (Bathonian; Fig. 11c) resulted in an intraoceanic convergence and eventually subduction of the oceanic lithosphere, followed by the arc and back-arc IAT and BABB magma generation during Late Jurassic and Early Cretaceous time ( Fig. 11d; Slovenec et al. 2011). This seems to have created an environment with high rates of water/rock interaction which, in turn, promoted hydrothermal processes in effusive and intrusive basic rocks. Low-to very-low-grade hydrothermal alteration phases resulted in the formation of chlorite, prehnite, pumpellyite, albite and titanite with minor tremolite, epidote, haematite, pyrrhotite, chalcopyrite and pyrite.

Conclusions
• Hydrothermal metamorphic processes in Mesozoic ophiolite and non-ophiolite basalts/andesites and tuffs of the active continental margin of the Western Neotethys took place at the ocean floor and, to a lesser degree, during the emplacement of the rocks. • Major hydrothermal alteration processes that resulted from the fluid-rock interaction include: albitization, chloritization, prehnitization, pumpellyitization, titanitization, epidotization, uralitization and zeolitization. In some locations these processes were limited, leaving relict primary igneous phases. • Hydrothermal metamorphic parageneses reflect a wide range of temperatures and pressures (from c. 450 to < 100°C and from c. 0.4 to < 0.001 GPa) spanning the greenschist facies through to pumpellyite-actinolite and prehnite-pumpellyite to the zeolite facies and diagenesis. • Initial Sr isotope ratios of the studied rocks show two trends in the intensity of the ocean-floor hydrothermal alteration: (a) during Triassic -Middle Jurassic (Bajocian) time, anorogenic ophiolite rocks (WPAB, E-MORB, T-MORB, N-MORB) formed during extension of oceanic lithosphere, leading to low to medium degrees of alteration, with intensity abating during Middle Triassic -Middle Jurassic (Bajocian) time; and (b) during Middle (Bathonian) to Late Jurassic through Lower Cretaceous time, orogenic ophiolites (IAT, BABB) formed following the establishment of the Neotethyan suprasubduction setting, where there was an apparent increase in the intensity of the hydrothermal alteration of the oceanic lithosphere. During Middle Triassic time, non-ophiolite basalts/andesites and tuffs (CAB) of active continental margin outlined a comparatively broad range of alteration intensities, devoid of any apparent decreasing/ increasing alteration trend. • The chemistry of magmatic host lithologies (i.e. non-ophiolite basalts/andesites and tuffs and ophiolite) dictated the hydrothermal alteration trajectories, leading to the prevalence of quartz and K-bearing phyllosilicates in the former, while the latter is enriched in Ca-silicates (prehnite, pumpellyite and laumontite) and rutile. • Chlorite, which is a widespread hydrothermal mineral, shows a broad range of compositions that depend on three variables: (a) protolith composition (clinopyroxene versus plagioclase), (b) octahedral/brucite layer cation exchange, and (c) geotectonic affinity of host rocks and rate of hydrothermal alteration. IAT and BABB Upper Jurassic and Lower Cretaceous extrusive and intrusive rocks show the highest degree of chloritization, which is consistent with the Sr isotopic data and petrographic observations that suggest enhanced alteration of clinopyroxene and plagioclase. • The phase chemistry of pumpellyite is indicative of a continuous decrease in crystallization temperatures at shallower depths of hydrothermal activity. • Epidote and deuteric amphibole occur almost exclusively in Upper Jurassic/Lower Cretaceous IAT and BABB rocks characterized by a high tectonic activity of the oceanic lithosphere and elevated crystallization temperatures of hydrothermal phases (e.g. chlorite) that belong to the prehnite-pumpellyite and pumpellyite-actinolite facies. • REEs in chlorite and pumpellyite in Triassic non-ophiolite basalts/andesites and tuffs of active continental margin call for very different hydrothermal solutions to explain their origin. The REE content of these solutions largely stems from the hydrothermal dissolution of plagioclase (and clinopyroxene); once mobilized by hydrothermal fluids, REEs are readily incorporated in newly formed chlorite and pumpellyite. Their variable REE contents are controlled by fluid temperature, composition, redox state and ligand complexation. • The major and REE geochemistry of critical hydrothermal phases that formed during ocean-floor metamorphism offers important insights into the evolution of the oceanic crust. For the Western Neotethys, where c. 150 Ma of hydrothermal history is elucidated starting with Middle Triassic ocean rifting and spreading until Middle Jurassic time, that evolved into Late Jurassic -Early Cretaceous intraoceanic subduction and back-arc extension, ocean-floor metamorphism peaked during the Middle-Late Jurassic intraoceanic subduction. Concurrent intensive tectonism likely facilitated the penetration of seawater into the recharge zone of the hydrothermal system to reach the hot reaction zone. This produced hightemperature bedrock metasomatism at prehnite-pumpellyite and pumpellyite-actinolite facies conditions. The Mesozoic hydrothermal ocean-floor metamorphism occurred prior to the obduction of the ophiolite mélange onto the passive continental margins of Adria during Early Cretaceous time.