1. 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, Reference Alt, Frey and Robinson2009; Alt & Teagle, Reference Alt, Teagle, Dilek, Moores, Elthon and Nicolas2000). 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, Reference Fyfe and Lonsdale1981; Alt & Teagle, Reference Alt, Teagle, Dilek, Moores, Elthon and Nicolas2000; Gillis et al. Reference Gillis, Coogan and Pedersen2005; Wilson et al. Reference Wilson, Teagle, Alt, Banerjee, Umino, Miyashita, Acton, Anma, Barr, Belghoul, Carlut, Christie, Coggon, Cooper, Cordier, Crispini, Durand, Einaudi, Galli, Gao, Geldmacher, Gilbert, Hayman, Herrero-Bervera, Hirano, Holter, Ingle, Jiang, Kalberkamp, Kerneklian, Koepke, Laverne, Vasquez, Maclennan, Morgan, Neo, Nichols, Park, Reichow, Sakuyama, Sano, Sandwell, Scheibner, Smith-Duque, Swift, Tartarotti, Tikku, Tominaga, Veloso, Yamasaki, Yamazaki and Ziegler2006). 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, Reference Mottl1983). 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. Reference Arbiol, Layne, Zanoni and Šegvić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, Reference Harlov and Marschall2009; Paoli et al. Reference Paoli, Dini, Petrelli and Rocchi2019). 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, Reference Seyfried and Bischoff1981; Mottl & McConachy, Reference Mottl and McConachy1990; German & Lin, Reference German and Lin2004; Lowell et al. Reference Lowell, Farough, Hoover and Cummings2013).
In the European Western Vardar ophiolites (in the sense of Schmid et al. Reference Schmid, Bernoulli, Fügenschuh, Matenco, Schefer, Schuster, Tischler and Ustaszewski2008; Fig. 1), or more precisely in the southwestern segment of the Zagorje-Mid-Transdanubian Zone (ZMTDZ; Pamić & Tomljenović, Reference Pamić and Tomljenović1998) or Sava Unit (Haas et al. Reference Haas, Mioč, Pamić, Tomljenović, Árkai, Bérczi-Makk, Koroknai, Kovács and Rálisch-Felgenhauer2000), Mesozoic mafic intrusive, extrusive and pyroclastic rocks crop out in northwestern Croatia (Slovenec et al. Reference Slovenec, Lugović, Meyer and Šiftar2011; Fig. 2). These rocks are incorporated as slices and metre-to-kilometre-sized blocks in ophiolite mélange (Slovenec & Šegvić, Reference Slovenec and Šegvić2021) known as the Kalnik Unit (KU) (Haas et al. Reference Haas, Mioč, Pamić, Tomljenović, Árkai, Bérczi-Makk, Koroknai, Kovács and Rálisch-Felgenhauer2000) (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. Reference Schmid, Bernoulli, Fügenschuh, Matenco, Schefer, Schuster, Tischler and Ustaszewski2008) based on the occurrence of non-ophiolite basalts/andesites and tuffs of the active continental margin (Slovenec & Šegvić, Reference Slovenec and Šegvić2021), others suggest that it forms the most southern portion of the Southern Alps (Haas et al. Reference Haas, Mioč, Pamić, Tomljenović, Árkai, Bérczi-Makk, Koroknai, Kovács and Rálisch-Felgenhauer2000; Schmid et al. Reference Schmid, Bernoulli, Fügenschuh, Matenco, Schefer, Schuster, Tischler and Ustaszewski2008; Fig. 1). The former proposal considers the KU to be an integral part of the Jurassic oceanic accretionary prism (e.g. Schmid et al. Reference Schmid, Bernoulli, Fügenschuh, Matenco, Schefer, Schuster, Tischler and Ustaszewski2008), 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. Reference Vlahović, Tišljar, Velić and Matičec2005). 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, Reference Alt, Frey and Robinson2009).
Simplified geotectonic map showing major units of the Alps, Carpathians and Dinarides (simplified after Schmid et al. Reference Schmid, Fügenschuh, Kounov, Maţenco, Nievergelt, Oberhänsli, Pleuger, Schefer, Schuster, Tomljenović, Ustaszewski and van Hinsbergen2020). The upper inset map shows the research area in more detail.
Simplified geological map of (a) Mounts Ivanščica and Strahinjščica, (b) Mounts Kuna Gora and Desinić Gora, (c) Mount Kalnik and (d) Mount Žumberak (modified after Pletničar et al. Reference Pletničar, Premru and Herak1975; Šikić et al. Reference Šikić, Basch and Šimunić1977; Šimunić et al. Reference Šimunić, Pikija and Hećimović1982; Aničić & Jureša, Reference Aničić and Jureša1984). Legend: 1 – Neogene and Pleistocene sedimentary rocks; 2 – Middle Cretaceous breccias, conglomerates, shales, marls, limestones and cherts; 3 – Upper Jurassic – Lower Cretaceous limestones; 4 – Jurassic limestones, cherts and dolomites; 5 – intercalation of Cretaceous calcarenite, marl and silicified marly limestone, shale and siltite, then ophiolite mélange with blocks of Triassic/Jurassic/Cretaceous basalt, gabbro and tectonite peridotite (green field) and Triassic–Jurassic radiolarites, sandstones and shales (not separated on the maps); 6 – Triassic limestones, dolomites, radiolarites, clastic rocks intersected by a series of andesite-basalts and/or tuffs (green field); 7 – reverse or thrust faults; 8 – normal faults; 9 – discordance line, tectonic-erosion discordance; 10 – sample location (a, ophiolitic rocks; b, non-ophiolitic rocks); 11 – picture break; 12 – quarry.
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. Reference Michard, Albarède, Michard, Minster and Charlou1983; Michard & Albarède, Reference Michard and Albarède1986; Bau, Reference Bau1991; Klinkhammer et al. Reference Klinkhammer, Elderfield, Edmond and Mitra1994; Bau & Dulski, Reference Bau and Dulski1999). Chlorite and pumpellyite are omnipresent in the oceanic crust preserved in the KU regardless of lithology (Slovenec et al. Reference Slovenec, Lugović and Slovenec2012). 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, Reference Brookins1989; Tetiker et al. Reference Tetiker, Yalçın and Bozkaya2015). 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.
2. 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. Reference Slovenec, Lugović, Meyer and Šiftar2011; Lugović et al. Reference Lugović, Slovenec, Schuster, Schwarz and Horvat2015). 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 continent-derived 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. Reference Šegvić, Kukoč, Dragičević, Vranjković, Brčić, Goričan, Babajić and Hrvatović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. Metre-to-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. Reference Slovenec, Lugović, Meyer and Šiftar2011). Tholeiitic cumulate gabbro is a minor lithology (Lugović et al. Reference Lugović, Slovenec, Schuster, Schwarz and Horvat2015), while alkali basalts and a high-grade metamorphic basement are only locally present (Slovenec et al. Reference Slovenec, Lugović, Meyer and Šiftar2011; Šegvić et al. Reference Šegvić, Lugović, Slovenec and Meyer2016). Five geochemically distinct types of mafic extrusives occur in the magmatic blocks of the KU (Slovenec et al. Reference Slovenec, Lugović, Meyer and Šiftar2011; Lugović et al. Reference Lugović, Slovenec, Schuster, Schwarz and Horvat2015): (i) within-plate alkaline basalts (WPAB; Middle Triassic, Illirian); (ii) enriched mid-ocean-ridge basalts (E-MORB; Middle Triassic, late Fassanian – early 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. Reference Slovenec, Lugović and Slovenec2012).
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. Reference Schmid, Bernoulli, Fügenschuh, Matenco, Schefer, Schuster, Tischler and Ustaszewski2008; 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, Reference Marci1987; Goričan et al. Reference Goričan, Halamić, Grgasović and Kolar-Jurkovšek2005; Slovenec & Šegvić Reference Slovenec and Šegvić2021). Alteration processes documented in non-ophiolite 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, Reference Marci1987). 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. Reference Goričan, Halamić, Grgasović and Kolar-Jurkovšek2005) and analogue successions from the Dinarides (Pamić, Reference Pamić1984; Smirčić et al. Reference Smirčić, Kolar-Jurkovšek, Aljinović, Barudžija, Jurkovšek and Hrvatović2018).
3. Materials and analytic methods
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.
Overview of the studied rock samples from NW Croatian mountains.
AN – andesite; AB – andesite-basalt; BA – basalt; BABB – back-arc basin basalts; BO – boninite; CA – calc-alkaline; DG – Desinić Gora Mount; DI – diabase; E-MORB – enriched mid-ocean-ridge basalts; EMPA – electron microprobe analyses; GB – gabbro; IAT – island-arc tholeiites; IV – Ivanščica Mount; KA – Kalnik Mount; KG – Kuna Gora Mount; LA-ICP-MS – laser ablation inductively coupled plasma mass spectrometry; ME –Medvednica Mount; N-MORB – normal mid-ocean-ridge basalts; OM – optical microscopy; RG – Ravna Gora Mount; SG – Samoborska Gora Mount; ST – Strahinjščica; T-MORB – transitional mid-ocean ridge basalts; TU – tuff (andesite/dacite/rhyolite); WPAB – within-plate alkaline basalts; XRD – X-ray diffraction; ŽU – Žumberak Mount.
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, Reference Pouchou and Pichoir1984) 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 (Reference Kranidiotis and MacLean1987) and the semi-empirical one of Inoue et al. (Reference Inoue, Inoué and Utada2018). 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. Reference Arbiol, Layne, Zanoni and Šegvić2021). The absolute accuracy of empirical geothermometers has been widely challenged (Arbiol et al. Reference Arbiol, Layne, Zanoni and Šegvić2021); the geothermometer of Kranidiotis & MacLean (Reference Kranidiotis and MacLean1987) has therefore only been used tentatively to show potential differences between the empirical and semi-empirical chlorite thermometers applied in hydrothermal systems. The semi-empirical chlorite geothermometer of Inoue et al (Reference Inoue, Inoué and Utada2018) is used here because it was initially developed for application to hydrothermal systems and assumes the equilibrium of quartz-chlorite-water. Ideally, the Fe3+/ΣFe should be measured or estimated, as required for the application of most semi-empirical geothermometers; however, the Inoue et al. (Reference Inoue, Inoué and Utada2018) geothermometer permits the calculation of the formation temperature of chlorite assuming ΣFe = Fe2+, 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. High-resolution 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 86Sr/88Sr = 0.1194 were assumed. The 87Sr/86Sr 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.e. X-ray diffractometry
X-ray powder diffraction (XRD) was carried out on 36 representative samples using a Bruker D8-Advanced diffractometer in the Department of Geosciences at Texas Tech University. The measurements were run using a step scan in the Bragg-Brentano geometry with CuKα radiation. The settings included 45 kV and 40 mA with sample mounts scanned from 3 to 70 °2θ. Measurements were completed under air-dried conditions at a counting time of 2.5 s per 0.02 °2θ. XRD data were processed using a Bruker EVA diffraction suite following recommendations on phyllosilicate interpretation detailed in Arbiol et al. (Reference Arbiol, Layne, Zanoni and Šegvić2021).
3.f. Laser ablation – inductively coupled plasma – mass 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-s-long 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 counts-per-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. Reference Jochum, Willbold, Raczek, Stoll and Herwig2005).
4. Results
4.a. Light microscope petrography
Petrographic studies of intrusive, effusive and pyroclastic rocks have previously been conducted (e.g. Slovenec et al. Reference Slovenec, Lugović, Meyer and Šiftar2011; Slovenec & Šegvić, Reference Slovenec and Šegvić2021) and, as a result, only the most pertinent information is provided here. Intrusive rocks from the KU ophiolite mélange are represented by cumulate and non-cumulate medium- to coarse-grained gabbro that shows a preserved igneous fabric and heteroadcumulate poikilitic, and/or anhedral granular texture and homogenous structure (Fig. 3a–c). The crystallization sequence of cumulate gabbro commenced with olivine, clinopyroxene (augite and diopside) and plagioclase (bytownite and anorthite), whereas amphibole (tschermakite, edenite and magnesiohornblende) and Fe-Ti oxide (magnetite-ulvöspinel) are late-stage crystallization products. The gabbro consists of plagioclase, augite and Fe-Ti oxide (Fig. 3b, c). Gabbroic rocks of the KU mélange bear an omnipresent prehnite-pumpellyite to greenschist facies metamorphic overprint that produced hydrothermal assemblages consisting of albite, tremolite/actinolite, chlorite, serpentine, prehnite, pumpellyite, titanite and muscovite (Fig. 3b, c).
Microphotographs of gabbros: (a) sample mc-2, N+; (b) sample vsk-228, N+; (c) sample be-3, N+ and basalts; (d) sample vsk-215/4, N+; (e) sample vk-372, N+; (f) sample got-6, N+; (g) sample vs-220, N+ from Mounts Kalnik, Ivanščica and Medvednica ophiolite mélange and andesite-basalt; (h) sample tskg-1/2b N+ from Mount Kuna Gora. Legend: Ab – albite; Afs – alkali feldspar; Am – amphibole; Chl – chlorite; Cpx – clinopyroxene (augite); Ep – epidote; Ill – illite; Ilm – ilmenite; Leu – leucoxene; NAb – neoalbite; Pl – plagioclase; Pmp – pumpellyite; Prh – prehnite; Py – pirite; Qtz – quartz; Ttn – titanite; Vg – volcanic glass.
Extrusive rocks of the KU ophiolite mélange consist of basaltic pillow lavas and massive lavas (Fig. 3d–g). Petrographic evidence suggests the following crystallization order: ± spinel (ulvöspinel-magnetite, chromite) → plagioclase (andesine to labradorite) → clinopyroxene (augite) + plagioclase (andesine to labradorite) ± Fe-Ti oxides (ilmenite). However, in the IAT-type lavas, plagioclase phenocrysts may enclose clinopyroxene indicating a somewhat different crystallization order: ± spinel (ulvöspinel-magnetite, chromite) → clinopyroxene (augite) → plagioclase (andesine to labradorite) ± Fe-Ti oxides (ilmenite). The primary assemblages are commonly altered to albite or peristerite, illite/mica, saussurite, calcite, chlorite, zoisite, epidote, leucoxene, analcime, prehnite, pumpellyite, quartz and zeolite (Fig. 3d–f), indicating the ocean-floor hydrothermal overprint in the zeolite – prehnite-pumpellyite – lower greenschist facies (Mével, Reference Mével2003).
Non-ophiolite Middle Triassic basalts/andesites and tuffs derived from an active continental margin are interstratified in volcanosedimentary successions of the north Croatian mountains (Figs 1, 2). These rocks include massive basic, intermediate to acid lavas and spatially associated fine-grained, pyroclastic rocks (i.e. tuffs). Volcanic rocks are characterized by ophitic through porphyritic and glomeroporphyritic seriate textures, which may gradually change to a fluidal texture (Fig. 3h). On the other hand, pyroclastic rocks show vitroclastic, vitro-crystalloclastic, litho-crystalloclastic and vitro-litho-crystalloclastic ash flow tuffs (i.e. ignimbrite) of crystallo-vitrophyre structure and fluidal texture. The volcanic rock microcrystalline matrix is holocrystalline to hypocrystalline and consists of devitrified volcanic glass and microlites of plagioclase/albite and K-feldspar with some minor Fe-Mg phases (Fig. 3h). The primary crystallization sequence is K-feldspar (sanidine) ± plagioclase (bytownite to anorthite) ± clinopyroxene (augite) ± amphibole ± Fe-Ti oxides (ilmenite, magnetite).
Accessory phases are zircon, apatite, rutile, haematite, pyrite, chalcopyrite and illite. The primary assemblage is commonly altered to illite/mica, albite, chlorite, calcite, leucoxene, quartz, and various clay minerals (Fig. 3h), while the glassy matrix altered into a mixture of fine-grained chlorite and mica whose occurrence is indicative of hydrothermal metamorphism (e.g. Tillick et al. Reference Tillick, Peacor and Mauk2001; Wang et al. Reference Wang, Lin, Xi, Cao, Wang, Yuan, Kashif and Song2018; Arbiol et al. Reference Arbiol, Layne, Zanoni and Šegvić2021). Conversely, the presence of clay-mineral-rich assemblages calls for somewhat lower temperatures that correspond to diagenetic and very low-grade metamorphic conditions (Alt, Reference Alt, Frey and Robinson2009; Šegvić et al. Reference Šegvić, Zanoni and Moscariello2020).
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 non-ophiolite 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.
Distribution and XRD/EMPA relative abundances of hydrothermal alteration phases in the studied Mesozoic ophiolite and non-ophiolite basalts/andesites and tuffs.
Ab – albite; AB – andesite-basalt; Am – amphibole; AN – andesite; Ant – anatase; Ap – apatite; BA – basalt; BABB – back-arc basin basalts; CAB – calc-alkaline basalts; Cal – calcite; Chl – chlorite; C-S – chlorite-smectite; Czo – clinozoisite; DI – diabase; Dol – dolomite; E-MORB – enriched mid-ocean-ridge basalts; Ep – epidote; GB – gabbro; Hal – halloysite; Hm – haematite; Hpy – chalcopyrite; IAT – island-arc tholeiites, Ilm – ilmenite; I/M – illite/mica; I-S – illite-smectite; Kln – kaolinite; Liz - lizardite; Mt – magnetite; N-MORB – normal mid-ocean-ridge basalts; Pal – palygorskite; Pmp – pumpellyite; Po – pyrrhotite; Prh – prehnite; Py – pyrite; Qtz – quartz; Rt – rutile; Sep – sepiolite; Sme – smectite; Stp – stilpnomelane; T-MORB – transitional mid-ocean-ridge basalts; Ttn – titanite; TU – tuff; Ver – vermiculite; WPAB – within-plate alkaline basalts; Zeo – zeolite; ++ – major phases; + – minor phases; * – phases present, but not unequivocally confirmed by XRD; x – phases detected by EMPA in samples not measured by XRD. Mineral abbreviations after Kretz (Reference Kretz1983).
4.b.1. Amphibole
Amphibole occurs in N-MORB (Early Jurassic, Pliensbachian–Aalenian), IAT (Middle–Late Jurassic, Callovian–Tithonian) and BABB (Early Cretaceous, Albian) isotropic and cumulate gabbro and gabbrodiorite from the KU mélange (Table 2; online Supplementary Table S1a). It largely consists of brown subhedral pleochroic magnesiohornblende (up to 2 mm) and anhedral pale green tremolite (Figs 3a, b, 4a–f, 5a). Rare fibres of Fe-anthophyllite marked by a faint pleochroism and minute size (< 0.02 mm) also occur (Figs 4d, 5b).
BSE microphotographs of gabbros: (a) sample vsk-228; (b) sample be-3; (c) sample vsi-8/1; (d) sample gtc-3; (e) sample gtc-5; (f) sample vsk-229/1 from Mounts Kalnik and Ivanščica ophiolite mélange. Numbered acquisition points correspond to the data in online Supplementary Tables S1a–S1e. Legend: Ab – albite; Am – amphibole; Ap – apatite; Chl – chlorite; Cpx – clinopyroxene (augite); Ep – epidote; Ilm – ilmenite; Pmp – pumpellyite; Prh – prehnite; Qtz – quartz; Rt – rutile; Stp – stilpnomelane; Ttn – titanite; Zrn – zircon.
Classification and discrimination diagrams for Ca monoclinic and orthorhombic Mg-Fe-Mn amphibole from the gabbroic and basaltic rocks from Mounts Kalnik and Ivanščica ophiolite mélange. (a, b) A(Na+K+2Ca – C(Al+Fe3++2Ti) plot (Hawthorne et al. Reference Hawthorne, Oberti, Harlow, Maresch, Martin, Schumacher and Welch2012) and (c) Na+K – Ti plot (Girardeau & Mevel, Reference Girardeau and Mevel1982).
In general, hydrothermal amphibole is characterized by a narrow range of Al2O3 (< 3.54 wt%; AlIV = 0.122–0.537 apfu), high FeO (17.00–34.82 wt%) and relatively low Mg no. (27.7–58.0; the lowest being in anthophyllite), as well as low total alkalis (Na+K < 0.5 apfu) and Ti (≤ 0.136 apfu) (Fig. 5c). Compared with tremolite, magnesiohornblende has a higher content of Cr2O3 and TiO2, while TiO2 and K2O abundances are the lowest in anthophyllite (Table 2; online Supplementary Table S1a).
4.b.2. Epidote
Epidote occurs only in IAT (Lower Jurassic) isotropic gabbro and BABB (Upper Jurassic – Lower Cretaceous) gabbrodiorite (Fig. 4a). Rarely, it occurs in N-MORB (Late Triassic – Middle Jurassic) and IAT (Late Jurassic) basalts from the KU ophiolite mélange (Table 2; online Supplementary Table S1a; Figs 3d, e, 6d). Epidote has constant concentrations of Fe3+ (0.880–0.944 apfu), AlVI (2.050–2.198 apfu) and Ca (1.882–1.952 apfu). The Mn-content of epidote is low (0.012–0.014 apfu) and is coupled with high concentrations of Fe3+ (Fig. 7a). In gabbro, xenomorphic epidote is accompanied by prehnite and chlorite to fill rare millimetre-sized cross-cutting veins.
BSE microphotographs of basalts: (a) sample vhk-11; (b) sample h-23; (c) sample vhk-47; (d) sample vsk-215/4; (e) sample vsk-221; (f) sample vsi-9 from Mounts Kalnik and Ivanščica ophiolite mélange. Numbered acquisition points correspond to the data in online Supplementary Tables S1a–1e. Legend: Ab – albite; Cal – calcite; Ccp – chalcopyrite; Chl – chlorite; Cpx – clinopyroxene (augite); Ep – epidote; Hm – haematite; Mag – magnetite; Pl – plagioclase; Pmp – pumpellyite; Prh – prehnite; Qtz – quartz; Rt – rutile; Stp – stilpnomelane; Ttn – titanite; Zeo – zeolite.
Classification diagram for (a) epidote group minerals (Armbruster et al. Reference Armbruster, Bonazzi, Akasaka, Bermanec, Chopin, Gieré, Heuss-Assbichler, Liebscher, Menchetti, Pan and Pasero2006); (b) feldspar (Ab–An–Or plot; Deer et al. Reference Deer, Howie and Zussman2013) from gabbroic and basaltic rocks from Mounts Kalnik and Ivanščica ophiolite mélange; and chlorite classification diagrams (c) Mg–Fe–VI(Al+Cr+Fe3+) (Zane & Weiss, Reference Zane and Weiss1998) and (d) Fe/(Fe+Mg)–Si (after Sun et al. Reference Sun, Zhang, Li, Hao, Wang, Li, Zhang, Zhang and Guo2019).
4.b.3. Feldspars
In ophiolites, hydrothermal albite (Table 2) is largely a product of a complete alteration of plagioclase (An73.9-15.5Ab82.9-26.1Or1.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 SiO2, Al2O3 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; Figs 3b–h, 4, 6). Chlorite is characterized by a narrow range of SiO2 and Al2O3 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 Fe2+/Mg steadily increased starting with the WPAB (Middle Triassic, Illirian), through E-MORB (Middle Triassic, upper Fassanian – lower Longobardian), T-MORB (Upper Triassic, middle–upper Carnian), N-MORB (Upper Triassic, upper Carnian – middle Norian; Middle Jurassic, Bajocian), IAT (Middle Jurassic, upper Bajocian – Oxfordian; Upper Jurassic, Callovian–Tithonian) to BABB (Lower Cretaceous, Albian) gabbro and/or basalt (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 Supplementary Fig. S1b). Amygdaloidal as well as matrix chlorite likely formed concomitantly with intensive albitization (Fig. 6e) and sericitization. Carbonate veins commonly cross-cut chlorite amygdules as well as the matrix (online Supplementary Fig. S1a), thus marking the youngest alteration event. These veins are made of Fe-dolomite, graphite and kaolin minerals, as well as inclusions of the matrix of the host rocks.
The empirical and semi-empirical geothermometers of Kranidiotis & MacLean (Reference Kranidiotis and MacLean1987) and Inoue et al. (Reference Inoue, Inoué and Utada2018) were applied to the compositions of chlorite. For that purpose, two broadly defined textural types of chlorite were identified (interstitial and vein). The former occurs in the interstices of pyroxene and plagioclase/albite or as their alteration product (Figs. 3b, e, f), while the latter fills circular amygdules and is present in microveins in the hypocrystalline matrix of various rock types (Fig. 6e). The geothermometer of Inoue et al. (Reference Inoue, Inoué and Utada2018) yields two distinct populations that are related to chlorite microtexture; interstitial chlorite shows temperatures of 92–489 °C (mostly 112 and 251 °C) and vein chlorite yields temperatures of 72–100 °C (online Supplementary Fig. S2). Geothermometric calculations using the Kranidiotis & MacLean (Reference Kranidiotis and MacLean1987) geothermometer returned comparable formation temperatures for vein chlorite, whereas calculated temperatures of interstitial chlorite are c. 50–100 °C lower than those yielded by the geothermometer of Inoue et al. (Reference Inoue, Inoué and Utada2018). Accordingly, we rely on the temperatures of Inoue et al. (Reference Inoue, Inoué and Utada2018) because the empirical geothermometer of Kranidiotis & MacLean (Reference Kranidiotis and MacLean1987) gives unrealistically low temperatures (online Supplementary Fig. S2).
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 stoichiometry of prehnite is practically ideal (online Supplementary 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, d, f). Fe-pumpellyite in BABB gabbro is characterized by a moderately high content of SiO2 (≤ 43.80 wt%) and variable amounts of FeO and MgO (2.78–13.91 wt% and 0.32–4.95 wt%, respectively), while the content of CaO is relatively constant (17.40–22.74 wt%) (online Supplementary Table S1d; 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. Reference Deer, Howie and Zussman2013).
Classification diagram for (a) pumpellyite, (b) titanite and (c) zeolite from gabbroic and basaltic rocks from Mounts Kalnik and Ivanščica ophiolite mélange. Pumpellyite Fetot–Mg–Al plot after Coombs et al. Reference Coombs, Nakamura and Vaugnat1976). Compositional fields of pumpellyite from the East Taiwan Ophiolite (zeolite facies), the Olympic Peninsula (prehnite–pumpellyite facies) and the Taveyannaz Formation (upper prehnite–pumpellyite and pumpellyite– actinolite facies) are taken from Coombs et al. (Reference Coombs, Nakamura and Vaugnat1976) and Rahn et al. (Reference Rahn, Mullis, Erdelbrock and Frey1994), respectively. The discrimination diagram (Ti/10)–Fe–Al for titanite with compositional fields of metamorphic facies/environments is taken from Calderón et al. (Reference Calderón, Prades, Hervé, Avendaño, Fanning, Massonne, Theye and Simonetti2013). Zeolite CaO+MgO – Na2O+K2O – 10×(SiO2/Al2O3) plot with compositional fields after Bastias et al. (Reference Bastias, Fuentes, Aguirre, Hervé, Demant, Deckart and Torres2016). Anl – analcime; Barr – barrerite; Cbz – chabazite; Cl – clinoptilolite; Dach – dachiardite; Epist – epistilbite; F – faujasite; Gi – gismondine; Hul – heulandite; Lmt – laumontite; Mes – mesolite; Mor – mordenite; Ntr – natrolite; O – offretite; Ph – phillipsite; Sc – scolecite; St – stellerite; Stb – stilbite; Tmp – thomsonite; Yu – yugawaralite.
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 (Na2O+ K2O < 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. Reference Campbell, Charnock, Dyer, Hillier, Chenery, Stoppa, Henderson, Walcott and Rumsey2016).
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, cryptocrystalline aggregates of leucoxene within plagioclase and chlorite (Fig. 3f). The composition of titanite defines it as grothite (Deer et al. Reference Deer, Howie and Zussman2013), 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. (Reference Lugović, Slovenec, Schuster, Schwarz and Horvat2015), Slovenec et al. (Reference Slovenec, Lugović, Meyer and Šiftar2011) and Slovenec & Šegvić (Reference Slovenec and Šegvić2021). The 87Sr/86Sr ratios show large variations (0.703208–0.718603) based on (a) the location of samples, (b) the type of geotectonic environment, and (c) the age of crystallization of the rocks and the degree of alteration.
Sr isotopic data of Mesozoic ophiolitic and non-ophiolitic basalts/andesites and tuffs from the mountains of NW Croatia.
a 87Rb/86Sr ratios calculated from Rb and Sr ICP-MS concentrations (displayed in ppm) following: 87Rb/86Sr = (Rb/Sr) × [2.6939 + 0.2832 × (87Sr/86Sr)].
b Initial 87Sr/86Sr(t) calculated assuming λ Rb = 1.42 × 10−11 a−1.
c Part of the data from Slovenec et al. (Reference Slovenec, Lugović, Meyer and Šiftar2011, Reference Slovenec, Lugović and Slovenec2012).
AB – andesite-basalt; BA – basalt; BABB – back-arc basin basalts; BN – boninite; CAB – calc-alkaline basalts; DI – diabase; E-MORB – enriched mid-ocean-ridge basalts; GB – gabbro; IAT – island-arc tholeiites; KA – Kalnik Mount; KG – Kuna Gora Mount; IV – Ivanščica Mount; ME – Medvednica Mount; N-MORB – normal mid-ocean-ridge basalts; SG – Samoborska Gora Mount; ST – Strahinjščica; T-MORB – transitional mid-ocean-ridge basalts; TU – tuff; WPAB – within-plate alkaline basalts; ŽU – Žumberak Mount.
The most prominent variations in Sr initial isotopic ratios ((87Sr/86Sr)t = 0.703195–0.710252; Fig. 9) are of Middle Triassic non-ophiolite basalts/andesites and tuffs, and pyroclastic rocks (CAB). A somewhat narrower variation range of Sr initial isotopic ratios (0.703725–0.708406) occurs in Middle–Upper Jurassic and Lower Cretaceous intrusive and extrusive orogenetic rocks (IAT, BABB; Fig. 9). On the other hand, Middle Triassic – Middle Jurassic mafic ophiolite intrusive and extrusive rocks that crystallized in anorogenic domains (WPAB, E-MORB, T-MORB, N-MORB) are characterized by a more narrow range of Sr initial isotopic ratios ((87Sr/86Sr)t = 0.703122–0.705632; Fig. 9). In these rocks, the (87Sr/86Sr)t values are generally the lowest, reflecting their weakest alteration intensities.
Initial 87Sr/86Sr isotope ratios diagram for Mesozoic ophiolitic and non-ophiolitic basalts/andesites and tuffs from the NW Croatia mountains. The record of Mesozoic seawater 87Sr/86Sr taken from Peucker-Ehrenbrink & Fiske (Reference Peucker-Ehrenbrink and Fiske2019). Part of the data originate from), and Slovenec et al. (Reference Slovenec, Lugović, Meyer and Šiftar2011, Reference Slovenec, Lugović and Slovenec2012).
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, Reference Boynton and Henderson1984; Fig. 10). The Eu anomaly is moderate (Euchon/Euchon* = 0.73–1.11; Euchon* = (Smchon+Gdchon) × 0.5) and the REE normalization curves are somewhat similar to N-MORB with a substantial depletion in LREE (Lachon/Ybchon = 0.26–0.33).
The chondrite-normalized plots (Boynton, Reference Boynton and Henderson1984) 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. Reference Michard, Albarède, Michard, Minster and Charlou1983) and seawater (Høgdahl et al. Reference Høgdahl, Melsom and Bowen1968).
Chondrite-normalized REE curves of euhedral chamosite phenocrysts show a REE enrichment (1–100 × chondrite), pronounced negative Eu anomalies (Eun/Eun* = 0.35–0.46) and negative Ce anomalies (Cechon/Cechon* = 0.001–0.010; Cechon* = (Lachon+Prchon) × 0.5) (Fig. 10). Chlorite is of hydrothermal origin in the non-ophiolite basalts (Table 2), as euhedral pumpellyite (online Supplementary Fig. S1d) in the hypocrystalline matrix consisting of chlorite and albite is also present. Chondrite-normalized REE patterns of pumpellyite are similar to those of boninite from ophiolite complexes (Saccani & Tassinari, Reference Saccani and Tassinari2015; Slovenec & Šegvić, Reference Slovenec and Šegvić2018) outlining concave-upwards (U-shape) profiles with a slight light REE (LREE; Lachon/Smchon = 3.14–3.71) and heavy REE (HREE; Tbchon/Luchon = 0.38–0.89) enrichment at levels of 9.01–36.06 times relative to chondrite. Abundances of medium REE (MREE) in this case are depleted with regards to HREE (Fig. 10). A prominent positive Eu anomaly (Euchon/Euchon* = 1.17–1.74) is also present.
5. Discussion
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, Reference Bickle and Teagle1992) 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. Reference Slovenec, Lugović, Meyer and Šiftar2011; 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. Reference Slovenec, Lugović, Meyer and Šiftar2011). 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. Reference Kovács, Demangel, Richoz, Hippler, Baldermann and Krystyn2020). 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. Reference Scambelluri, Rampone and Piccardo2001). 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 (16–49) 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. (Reference Slovenec, Lugović, Meyer and Šiftar2011). 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.
Schematic geodynamic model for interaction of (a) active continental margin magmatic activity (ensialic volcanic arc magmatism) and roughly contemporaneous proto back-arc rifting, (b) active spreading, (c) subduction of an active ocean ridge and (d) evolved subduction stage with initiation of an island arc and formation of the subsequent ensialic marginal (back-arc) basin in Meliata-Maliak-Dinaric-Vardar ocean system (modified from Slovenec et al. Reference Slovenec, Lugović, Meyer and Šiftar2011, Lugović et al. Reference Lugović, Slovenec, Schuster, Schwarz and Horvat2015; Slovenec & Šegvić, Reference Slovenec and Šegvić2021). Scale is approximate. IA – island arc (infant proto-arc/island-arc system); BAB – back-arc basin; AP – accretionary prism; WPAB – within-plate alkaline basalts; BABB – back-arc basin basalts; CAB – calc-alkaline basalts; E-MORB – enriched mid-ocean-ridge basalts; IAT – island-arc tholeiites; N-MORB – normal mid-ocean-ridge basalts; OIB – ocean-island basalts; T -MORB – transitional mid-ocean ridge basalts.
Hydrothermal sulphide-oxide systems are essentially controlled by mineral-fluid equilibria, which may provide valuable information on the prevailing chemical and physical characteristics of the mineralizing fluid such as elemental composition, temperature, pH, fO2, dilution and cooling rates (e.g. Rottier et al. Reference Rottier, Kouzmanov, Wälle, Bendezú and Fontboté2016). Middle Triassic – Lower Cretaceous ophiolite rocks show variable amounts of sulphides and oxides (Table 2; online Supplementary Table S1e). Weakly to moderately altered non-orogenic Ladinian–Bajocian ophiolite rocks (WPAB, E-MORB, T-MORB and N-MORB; Fig. 11a,b) contain either sulphides (pyrite or pyrrhotite) or oxides (haematite or magnetite). Conversely, orogenic 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. Reference Henley, Truesdell, Barton and Whitney1984). Cooling is normally caused by a drop in pressure associated with fluid ascending in the hydrothermal system or cold seawater mixing (Capuano & Cole, Reference Capuano and Cole1982). This in turn makes the pH of hydrothermal fluid decrease following sulphide precipitation (Akaku et al. Reference Akaku, Reed, Yagi, Kai and Yasuda1991). 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. Reference Henley, Truesdell, Barton and Whitney1984). 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. Reference Fouquet, Pelleter, Konn, Chazot, Dupré, Alix, Chéron, Donval, Guyader, Etoubleau, Charlou, Labanieh and Scalabrin2018). 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. Reference Fujibayashi, Asakura, Hattori and Allen2014). 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. Reference Kuhn, Versteegh, Villinger, Dohrmann, Heller, Koschinsky, Kaul, Ritter, Wegorzewski and Kasten2017). 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. Reference Slovenec, Lugović and Slovenec2012).
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. Reference Von Damm, Edmond, Grant, Measures, Walden and Weiss1985; Ray et al. Reference Ray, Mevel and Banerjee2009). 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. Reference Inoue, Inoué and Utada2018), 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 (Fig. 11a–d). Accordingly, a clear downwards trend in the Mg content (clinochlore) and an increase in Fe2+ (chamosite) is obtained starting with WPAB (Middle Triassic, Illirian) through E-MORB (Middle Triassic, upper Fassanian – lower Longobardian), T-MORB (Upper Triassic, middle–upper Carnian), N-MORB (Upper Triassic, upper Carnian – middle Norian to Middle Jurassic, Bajocian), IAT (Middle Jurassic, upper Bajocian – Oxfordian; Upper Jurassic, Callovian–Tithonian) and BABB (Lower Cretaceous, Albian; Fig. 7c, d) gabbro and/or basalt. This correlates well with the petrographic observations indicating a more intensive alteration of clinopyroxene and albitization of plagioclase coupled with increased abundances of chlorite in Upper Jurassic and Lower Cretaceous extrusive and intrusive rocks compared with those of Middle Triassic and Middle Jurassic age (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. Reference Slovenec, Lugović, Meyer and Šiftar2011) with abundant crystallization of prehnite (online Supplementary Table S1d) in the range 190–340 °C (Liou et al. Reference Liou, Maruyama and Cho1985), 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. (Reference Liou, Maruyama and Cho1985) 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, Reference Schiffman and Liou1983). Conversely, less common Al-pumpellyite (Al > 2.290 apfu) with a high ratio of Al/Fetot (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. This supports the idea of prehnite-actinolite facies conditions at comparatively higher alteration depths (Liou et al. Reference Liou, Maruyama and Cho1985). 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, fO2 exercises a further control on pumpellyite stability (Inoue & Utada, Reference Inoue and Utada1991). Indeed, the formation of Al-pumpellyite occurs at lower fO2 (Schiffman & Liou, Reference Schiffman and Liou1983), 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. (Reference Enami, Suzuki, Liou and Bird1993) who revealed that, in metamorphic systems, the average Al + Fe3+ 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. Reference Liou, Maruyama and Cho1985) 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. Reference Deer, Howie and Zussman2013). 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, Reference Ilijima and Utada1972).
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, Reference Coleman1977), 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. Reference Liou, Maruyama and Cho1985). 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. Reference Slovenec, Lugović, Meyer and Šiftar2011; Lugović et al. Reference Lugović, Slovenec, Schuster, Schwarz and Horvat2015; 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 tschermakite-magnesiohornblende 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 Fe3+ in the octahedral site (i.e. pistacite component) of epidote (Fe3+/(Fe3++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. Reference Liou, Maruyama and Cho1985). This correlates well with the experimental results of Berndt et al. (Reference Berndt, Seyfried and Janecky1989), 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. Reference Deer, Howie and Zussman2013). 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, Reference Ishizuka1991). The formation of pumpellyite, irrespective of its chemistry, denotes the transition between the zeolite and greenschist facies (e.g. Kamimura et al. Reference Kamimura, Hirajima and Fujimoto2012).
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 SiO2 content for ophiolites ranges over c. 46–50 wt%, TiO2 is c. 2 wt%, while CaO and K2O contain up to 11 and 1.5 wt%, respectively (e.g. Slovenec et al. Reference Slovenec, Lugović, Meyer and Šiftar2011; Slovenec & Šegvić, Reference Slovenec and Šegvić2018). Conversely, the composition of non-ophiolitic basalts/andesites and tuffs (e.g. Goričan et al. Reference Goričan, Halamić, Grgasović and Kolar-Jurkovšek2005; Slovenec & Šegvić, Reference Slovenec and Šegvić2021) is marked by higher levels of SiO2 (c. 53–57 wt%) and K2O (c. 3–6 wt%), whereas CaO and TiO2 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 TiO2-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, Reference Schiffman and Staudigel1994), 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. Reference Slovenec, Lugović, Meyer and Šiftar2011; Šegvić et al. Reference Šegvić, Kukoč, Dragičević, Vranjković, Brčić, Goričan, Babajić and Hrvatović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. Reference Morad, El-Ghali, Caja, Sirat, Al-Ramadan and Mansurbeg2010). 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. Reference Judik, Árkai, Horváth, Dobosi and Al2005).
5.b. Chlorite and pumpellyite REE geochemistry
The normalized REE content of seawater (Høgdahl et al. Reference Høgdahl, Melsom and Bowen1968) and that of modern hydrothermal fluids (black smoker) (Michard et al. Reference Michard, Albarède, Michard, Minster and Charlou1983; Michard & Albarède, Reference Michard and Albarède1986) 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. Reference Michard, Albarède, Michard, Minster and Charlou1983), while in low-temperature solutions, the REE content remains low and will not alter regardless of any reasonable fluid/rock ratios (Sanjuan et al. Reference Sanjuan, Michard and Michard1988). 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, Reference Valsami and Cann1992) 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, Reference Moeck2014), while periods of activity may be separated by thousands of years of quiescence (Heasler et al. Reference Heasler, Jaworowski, Foley, Young and Norby2009). 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. Reference Klinkhammer, Elderfield, Edmond and Mitra1994; Şener et al. Reference Şener, Şener and Uysal2017). Also, the analogue LREE enrichment is documented in hydrothermal chlorite formed from such high-temperature solutions (Tetiker et al. Reference Tetiker, Yalçın and Bozkaya2015). Our line of reasoning therefore advocates a deficiency of LREEs in low-temperature solutions from which chlorite amygdules formed. Should LREEs have been abundant they would, depending upon favourable water/mineral partitioning, be accommodated by nascent chlorite either in its structure or at reactive surfaces (Xiao & Chen, Reference Xiao and Chen2020; Tan et al. Reference Tan, Mao, Yu, Sun and Lv2021). Since this is not the case, it is reasonable to hypothesize that parental fluid for the formation of chlorite was low-temperature and mildly acidic. In such solutions, REEs are mostly complexed; LREEs are however less complexed than HREEs with the former displaying considerable complexing by Cl−, SO4 2− and CO3 2−, while the latter is commonly complexed by CO3 2− (Valsami & Cann, Reference Valsami and Cann1992; Li et al. Reference Li, Webb, Algeo, Kershaw, Lu, Oehlert, Gong, Pourmand and Tan2019; Gong et al. Reference Gong, Li, Lu, Wang and Tang2021). The abundance of Upper Palaeozoic and Triassic carbonates interlayered with the studied extrusive rocks and tuffs (Goričan et al. Reference Goričan, Halamić, Grgasović and Kolar-Jurkovšek2005; Slovenec & Šegvić, Reference Slovenec and Š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 fluid-facilitated fast chemical reactions that are not in equilibrium with the fluid (e.g. Jamtveit et al. Reference Jamtveit, Austrheim and Putnis2016). Triassic extrusive rocks of the Western Neotethys affected by hydrothermal activity are largely composed of plagioclase, clinopyroxene, K-feldspar and multiple minor phases (Marci, Reference Marci1987; Slovenec & Šegvić, Reference Slovenec and Šegvić2021). With the exception of plagioclase and clinopyroxene, none of the other phases show any significant REE partitioning (Rollinson & Pease, Reference Rollinson and Pease2021). Moreover, compared with clinopyroxene, plagioclase is two to three times more abundant in the studied rocks (Marci, Reference Marci1987) and dissolves preferentially in hydrothermal fluids (Slovenec & Šegvić, Reference Slovenec and Š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. Reference Klinkhammer, Elderfield, Edmond and Mitra1994; Bau & Dulski, Reference Bau and Dulski1999). Several studies have demonstrated remarkable similarities in REE abundances between hydrothermal fluids and plagioclase phenocrysts (Schnetzler & Philpotts, Reference Schnetzler and Philpotts1970; Papike et al. Reference Papike, Fowler, Shearer and Layne1996) 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 Cl2− in deep-seated acidic hydrothermal fluids, whereby HREEs tend to be less fractionated in hydrothermal precipitates (Flynn & Wayne Burnham, Reference Flynn and Wayne Burnham1978; Allen & Seyfried, Reference Allen and Seyfried2005). Gadolinium, in contrast, prefers the solid phase (Fig. 10), which corroborates the acidic nature of the percolating fluids (Zielinski & Frey, Reference Zielinski and Frey1974). 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 Ce4+ to form, which in turn led to the impoverishment of Ce in the fluid relative to other REE (Sanematsu et al. Reference Sanematsu, Moriyama, Sotouky and Watanabe2011). According to the mechanism proposed by Valsami & Cann (Reference Valsami and Cann1992), 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. Reference Henley, Truesdell, Barton and Whitney1984).
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, Reference Xiao and Chen2020; Tan et al. Reference Tan, Mao, Yu, Sun and Lv2021).
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 Eu2+ and thus fractionated from the other REEs (e.g. Bau, Reference Bau1991; Bau & Dulski, Reference Bau and Dulski1999).
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. Reference Lustrino, Abbas, Agostini, Gaggiati, Carminati and Gianolla2019). 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. Reference Slovenec, Lugović, Meyer and Šiftar2011). 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. Reference Slovenec, Lugović, Meyer and Šiftar2011). 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. Reference Slovenec, Lugović, Meyer and Šiftar2011). 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.
6. 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 high-temperature 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.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S0016756822001030
Acknowledgments
This work was supported by the Croatian Science Foundation (project IP-2019-04-3824) and the Croatian Ministry of Science, Education and Sport (grant no. 181-1951126-1141). Further support was obtained from the Geosciences Clay Laboratory at Texas Tech University. We are especially grateful to Boško Lugović (University of Zagreb) and Hans-Peter Meyer (Universität Heidelberg) for providing excellent microprobe data, as well as Illona Fin (Universität Heidelberg) for preparing high-quality polished thin-sections. We extend our appreciation to Kevin Werts for his guidance in the GeoAnalytical Laboratory at Texas Tech University. We are also indebted to James Browning (Texas Tech University) for his help with sample preparation for LA-ICP-MS analyses. Bruno Tomljenović (University of Zagreb) provided valuable materials needed for the preparation of Figure 1. We also acknowledge Carlos Arbiol (COREM, Québec) for fruitful discussions on chlorite chemistry. Aleksandar Ristić and Kevin Byerly are thanked for assistance with the English language. Critical comments and constructive reviews by Ben Tutolo, Dionysis Foustoukos and Ondrey Nemec as well as editorial assistance by Paul Spry and Tim Johnson have contributed significantly to the manuscript quality.
Open access
