Differentiating between hydrothermal and diagenetic carbonate using rare earth element and yttrium (REE+Y) geochemistry: a case study from the Paleoproterozoic George Fisher massive sulfide Zn deposit, Mount Isa, Australia

Carbonate minerals are ubiquitous in most sediment-hosted mineral deposits. These deposits can contain a variety of carbonate types with complex paragenetic relationships. When normalized to chondritic values (CN), rare-earth elements and yttrium (REE+YCN) can be used to constrain fluid chemistry and fluid-rock interaction processes in both low- and high-temperature settings. Unlike other phases (e.g., pyrite), the application of in situ laser ablation-inductively coupled plasma-mass spectroscopy (LA-ICP-MS) data to the differentiation of pre-ore and hydrothermal carbonates remains relatively untested. To assess the potential applicability of carbonate in situ REE+Y data, we combined transmitted light and cathodoluminescence (CL) petrography with LA-ICP-MS analysis of carbonate mineral phases from (1) the Proterozoic George Fisher clastic dominated (CD-type) massive sulfide deposit and from (2) correlative, barren host rock lithologies (Urquhart Shale Formation). The REE+YCN composition of pre-ore calcite suggests it formed during diagenesis from diagenetic pore fluids derived from ferruginous, anoxic seawater. Hydrothermal and hydrothermally altered calcite and dolomite from George Fisher is generally more LREE depleted than the pre-ore calcite, whole-rock REE concentrations, and shale reference values. We suggest this is the result of hydrothermal alteration by saline Cl--rich mineralizing fluids. Furthermore, the presence of both positive and negative Eu/Eu* values in calcite and dolomite indicates that the mineralizing fluids were relatively hot (>250°C) and cooled below 200–250°C during ore formation. This study confirms the hypothesis that in situ REE+Y data can be used to differentiate between pre-ore and hydrothermal carbonate and provide important constraints on the conditions of ore formation.


Introduction
In many mineral systems, carbonate minerals constitute a significant proportion of the host rock and gangue mineralogy. Complex carbonate parageneses are often found in ore deposits hosted in sedimentary rocks (Cline et al. 2005;Hitzman et al. 2005;Leach et al. 2005), and can derive from detrital, diagenetic, hydrothermal, and tectonic or metamorphic processes. Such carbonate-rich ore systems also typically lack an extensive alteration footprint due to the effective buffering of acidic hydrothermal fluids by carbonate-rich lithologies. The complex carbonate paragenesis and the limited alteration footprint are major challenges when investigating, and exploring for, ore systems in carbonate-rich rocks.
Recent analytical developments are providing exploration geoscientists with new techniques for confronting these challenges. For example, new techniques that enable the generation of large C and O isotope datasets have shown how alteration halos in carbonate-hosted mineral systems may extend beyond bulk-rock lithogeochemical anomalies (e.g., off-axis integrated cavity output spectroscopy; Barker et al. 2013). As a result of being generated using bulk-rock techniques, such isotopic datasets may need to be supplemented by in situ Editorial handling: B. Lehmann * Philip Rieger phil.rieger.geo@gmail.com techniques in order to constrain hydrothermal and background processes in samples with a complex paragenetic history. In low-temperature (e.g., marine or diagenetic) environments, in situ trace element analysis of carbonate minerals (e.g., laser ablation ICP-MS) has attracted increasing usage (e.g., Webb et al. 2009;Himmler et al. 2010). As the partitioning of rare-earth elements and yttrium (REE+Y) between fluids and carbonate minerals is well constrained (e.g., Morgan and Wandless 1980;Tanaka and Kawabe 2006; Gysi 2020) the fractionation of REE+Y can provide important constraints on physicochemical conditions (pH, temperature, ligand type, or fluid-rock interaction; e.g., Bau 1991;Bau and Möller 1992;Barker et al. 2006). The REE+Y composition of carbonate can also vary significantly due to intrinsic processes such as closed-system effects and the carbonate growth rate (e.g., Kontak and Jackson 1995;Barker and Cox 2011). Altogether, there is potential for the REE+Y systematics of diagenetic and hydrothermal carbonates in constraining key aspects of fluid chemistry. Nevertheless, the direct application of in situ REE+Y analysis of carbonate minerals has remained relatively untested for hydrothermal mineral systems (e.g., Magnall et al. 2016;Vaughan et al. 2016), despite the wealth of information that is available from these types of analyses (Debruyne et al. 2016;Smrzka et al. 2019).
The need to differentiate between pre-ore and hydrothermal carbonate signatures is particularly relevant in the Carpentaria province (Fig. 1), where several of the world's largest clastic-dominated (CD-type) Zn-Pb massive sulfide deposits are hosted in sedimentary rocks with multiple types of carbonate (e.g., Perkins and Bell 1998;Large et al. 2005;McGoldrick et al. 2010). Previous bulk-rock lithogeochemical studies have found large-scale (several km) halos of carbonate-associated elements (Mn and Fe) related to these hydrothermal systems (Large and McGoldrick 1998;Large et al. 2000). However, both Fe and Mn can also be enriched in carbonate minerals by background processes in marine or diagenetic environments (e.g., Brand and Veizer 1980;Kah 2000;Wittkop et al. 2020). As a result, bulk-rock and conventional in situ major-element analyses (i.e., electron probe microanalysis) alone may not be sufficient to differentiate between hydrothermal and background carbonate signatures.
In this study, we test the hypothesis that in situ REE+Y data can be used to differentiate diagenetic and hydrothermal carbonate mineral phases. We combine transmitted-light and cathodoluminescence (CL) petrography with LA-ICP-MS trace-element analyses of carbonates from drill-core samples through the main ore bodies at the southern Carpentaria George Fisher deposit (165 Mt at 9.1% Zn, 3.4% Pb, and 55 g/t Ag; Glencore 2019) and from a drill hole that intersected correlative, unmineralized host-rock lithologies (Urquhart Shale Formation). The in situ data are then compared to REE+Y compositions from whole-rock data from the same samples.

Rare-earth elements and yttrium (REE+Y)
The REE+Y are hard cations (high charge/radius ratio) and all have broadly similar chemical behavior, although subtle differences in solubility occur due to decreasing ionic radii from the light to the heavy REE+Y (e.g., Bau 1991;Migdisov et al. Gibson et al. 2017). The locations of the George Fisher, Hilton, and Mount Isa deposits and the unmineralized drill-hole (Shovel Flats) are denoted by stars and a circle respectively. d Stratigraphy of the Mount Isa Group (approximate formation thicknesses adapted from Neudert 1983) 2009; Williams-Jones et al. 2012). As a consequence of the Oddo-Harkins effect, elements with even atomic numbers are more abundant than those with uneven atomic numbers; therefore, the REE+Y are typically normalized to reference values (e.g., post-Archean Australian shale - PAAS, McLennan 1989;chondrite, McDonough and Sun 1995), which allows for better evaluation of the relative fractionation from different reservoirs and among individual REE+Y subgroups (light, middle, and heavy rare earth elements; LREE, MREE, and HREE). This fractionation can then be used to trace biogeochemical processes and to constrain fluid chemistry in lowand high-temperature settings (for reviews see Debruyne et al. 2016 andSmrzka et al. 2019). The REE+Y are mostly trivalent (3+), apart from Ce (3+, 4+), Eu (2+, 3+), and Yb (2+, 3+), which are redox sensitive. The redox sensitivity of Yb is only important at high temperatures (> 420°C), whereas Ce and Eu can undergo valence changes at ambient and lower hydrothermal temperatures under geochemically reasonable fO 2 conditions (for example at pH = 3; Fig. 2b). The fractionation of Ce and Eu relative to their neighbouring REE can, therefore, be used to constrain fluid redox conditions (e.g., De Baar et al. 1983;Sverjensky 1984;Bau and Möller 1992). This contrasting behavior is expressed as empirical anomalies relative to reference values, which can either be determined linearly or geometrically (e.g., Lawrence et al. 2006). The Ce anomalies (Ce/Ce * N ) are typically shale-normalized (e.g., to PAAS; Ce/Ce * SN ) in order to evaluate seawater or pore water redox, whereas Eu anomalies are commonly chondrite-normalized (Eu/Eu * CN ) to avoid overestimation of Eu/Eu * values due to the negative Eu anomaly of PAAS (Fig. 2a).
Yttrium and Ho exhibit similar geochemical behavior, although Ho is preferentially complexed and removed from fluids by Fe-oxides and organic particles; consequently, seawater and seawater-derived fluids preserve super-chondritic (>28) Y/Ho values, whereas crustal fluids have chondritic Y/ Ho values (Zhang et al. 1994;Nozaki et al. 1997;Bau and Dulski 1999). Similarly, there are subtle differences in the adsorption of LREE, MREE, HREE, and Ce onto organic particles, oxide phases, or mineral phases (e.g., German and Elderfield 1989;Sholkovitz et al. 1994;Alibo and Nozaki 1999). These differences affect the absolute and relative abundances of REE+Y in seawater as a function of redox state and water depth; further modification can then occur due to mineral formation or redox reactions in early diagenetic pore waters (e.g., Haley et al. 2004;Abbott et al. 2019; for REE+Y reference profiles see Fig. 2a and for a recent review on REE+ Y fractionation see Smrzka et al. 2019).
Ligand concentration, pH, and temperature can have strong effects on REE+Y solubility in aqueous solutions (e.g., Michard 1989;Craddock et al. 2010;Williams-Jones et al. 2012) and the REE+Y composition of hydrothermal fluids can be extremely variable (Fig. 2a). Fluid-mineral interactions under hydrothermal conditions can also modify the REE+Y signatures of carbonate mineral phases (Vaughan et al. 2016). Commonly, hydrothermal fluids have positive Eu/Eu * and REE+Y CN profiles that are LREE enriched, which has been interpreted to directly result from dissolution of plagioclase during fluid-rock interaction (e.g., Klinkhammer et al. 1994;Douville et al. 1999). Empirical and experimental studies suggest that both the LREE and Eu 2+ are more soluble in high temperature Cl --rich fluids (Migdisov et al. 2009;Craddock et al. 2010;Williams-Jones et al. 2012). Because of the higher solubility of Eu 2+ relative to Eu 3+ , the temperature-dependent redox sensitivity of Eu has a strong effect on Eu/Eu * ; for example, under geochemically reasonable fO 2 and acidic conditions the more soluble Eu 2+ is only dominant over Eu 3+ at temperatures above 200-250°C (Fig. 2b;Sverjensky 1984;Bau 1991;Bilal 1991). As a consequence, positive Eu/Eu * can only develop in hot (>200-250°C) hydrothermal fluids or in extremely reducing environments. Divalent Eu, however, is not readily incorporated into the carbonate crystal lattice, which is why positive Eu/Eu * can only develop in carbonate mineral phases that formed from hot fluids that cooled below 200-250°C in order to stabilize Eu 3+ (Bau and Möller 1992). Many hydrothermal carbonates also preserve LREE depletion relative to PAAS (e.g., Roberts et al. 2009;Debruyne et al. 2013;Magnall et al. 2016). Such LREE-depleted REE+Y signatures in hydrothermal precipitates have been interpreted to result from co-precipitation with LREE-enriched phases (Roberts et al. 2009), LREE scavenging by fluid-mineral interaction along the fluid pathway (e.g., by monazite; Debruyne et al. 2016), inherited REE+Y signatures from fluid-rock interaction (Lüders et al. 1993;Hecht et al. 1999), or LREE retention in the fluid by Cl --complexes (Craddock et al. 2010;Magnall et al. 2016;Perry and Gysi 2018).

Geological background
The Carpentaria Province comprises the Mount Isa Inlier and the McArthur Basin, which formed in an intracontinental setting during the Paleo-to Mesoproterozoic (Betts et al. , 2016Giles et al. 2002). During this time, three unconformitybound superbasins formed due to episodic rifting, sag phase, and basin inversion; those superbasins can be further subdivided into twelve supersequences ( Fig. 1; Jackson et al. 2000;Southgate et al. 2000Southgate et al. , 2013. The George Fisher, Hilton, and Mount Isa massive sulfide deposits are hosted by the Urquhart Shale Formation, which is part of the Mount Isa Group of the Isa Superbasin (Southgate et al. , 2013. Closure of the superbasins was initiated by the onset of the Isan orogeny (ca. 1600 Ma), during which the rocks of the Mount Isa Inlier underwent multiple phases of deformation in at least four separate north-south or east-west directed deformation events (Page and Bell 1986;Connors and Page 1995;Bell and Hickey 1998).
The Urquhart Shale Formation was deposited at ca. 1652 ± 7 Ma and comprises mainly laminated to bedded mudstones and siltstones, which are variably calcareous, dolomitic, carbonaceous, and pyritic (Page and Sweet 1998;Painter et al. 1999;Page et al. 2000). The depositional environment of the Urquhart Shale Formation was interpreted as either a carbonate slope (Neudert 1983;Painter, 2003) or deeper-water environment . Moderate concentrations of organic carbon (< 2 wt.%) and Mo (< 30 ppm) are consistent with anoxic and ferruginous seawater, which is further supported by sulfur isotope values of diagenetic pyrite that do not preserve evidence of extreme sulfate limitation under euxinic conditions (Rieger et al. 2020a(Rieger et al. , 2021a. Carbonates are major rock-forming mineral phases in the Urquhart Shale. There are several types of detrital, diagenetic, and hydrothermal calcite and dolomite. Zoned, porous, siltsized detrital clasts of calcite and dolomite were deposited together with other clastic mineral grains (Painter et al. 1999) and may have been entirely replaced by diagenetic carbonate (Chapman 1999). Diagenetic carbonates comprise calcareous and dolomitic cements, which have been interpreted as supratidal crusts (e.g., Neudert 1983;Painter et al. 1999), seafloor cement , or as carbonate alteration (Chapman 1999). Furthermore, nodular carbonates have been interpreted as diagenetic precipitates (e.g., Chapman 1999Chapman , 2004Domagala et al. 2000;Painter, 2003), as diagenetic pseudomorphs after sulfate minerals (e.g., van den Heuvel 1969;McClay and Carlile 1978;Painter et al. 1999), or as hydrothermal replacement textures (Perkins 1997;Perkins and Bell 1998).
In and around the George Fisher, Hilton, and Mount Isa deposits there are several types of calcite, dolomite, and siderite veins, infill, and breccias associated with the Zn-Pb and Cu ore systems (e.g., Waring 1990;Valenta 1994;Chapman 1999). Generally, carbonate mineral phases are more dolomitic and Fe-and Mn-bearing at the Mount Isa and George Fisher deposits relative to the un-mineralized Urquhart Shale (Painter, 2003;Rieger et al. 2021a). Furthermore, the ore forming systems have likely produced 18 O-depleted carbonate haloes around the mineral deposits in the Urquhart Shale Formation beyond visible alteration (Waring 1990;Chapman 1999). Besides differences in carbonate mineralogy and isotopic composition, the footprints of these CD-type massive sulfide deposits are evident from element enrichment (e.g., Tl, Mn, Ag) and depletion (e.g., Na, Sr) in whole-rock data (Painter, 2003;Rieger et al. 2021a). The main sulfide minerals at George Fisher are formed in multiple paragenetic stages and comprise pyrite, sphalerite, galena, and pyrrhotite, with minor quantities of chalcopyrite (Chapman 1999(Chapman , 2004Murphy 2004;Rieger et al. 2020a). Recently, Rieger et al. (2020a) described pre-ore diagenetic fine-grained pyrite followed by three ore stages: (1) stratabound sphalerite + pyrite ± galena, (2) breccia-hosted galena + sphalerite + pyrite + pyrrhotite, and (3) vein and breccia-hosted pyrite + pyrrhotite + chalcopyrite ± galena and sphalerite.
A whole-rock lithogeochemistry study on Urquhart Shale samples from an unmineralized drill-core and from the George Fisher deposit has shown that REE+Y, Si, Al, Ti, Sc, Nb, and Th were immobile during hydrothermal alteration, and that mineralization at George Fisher resulted in dilution of these immobile elements (Rieger et al. 2021a). The nature of any REE+Y fractionation was not discussed in Rieger et al. (2021a) and the data are presented here in Fig. 3 and are also available in Rieger et al. (2020b), and fully discussed below. The REE+Y CN profiles for Urquhart Shale and massive sulfide samples are LREE enriched relative to chondrite, and have relatively flat MREE and HREE profiles, as well as negative Eu/Eu * and chondritic Y/Ho values (Fig. 3). These REE+Y CN profiles are generally very similar to PAAS and to other formations from the Mount Isa Group, yet, are offset to lower concentrations ( Fig. 3; Nance and Taylor 1976;McLennan 1989).
Drill-core samples were examined using a binocular microscope and key samples were selected for polished thin-section preparation (n = 95). Thin sections were further examined using transmitted-and reflected-light microscopy and representative carbonate samples were selected for cathodoluminescence (CL) petrography (n = 30). The hotcathode optical CL system was operated at 14 keV and 0.15-0.20 mA. To capture zonation and differences in luminescence, exposure times were varied and quartz grains were used as reference for luminescence intensity, because quartz is characterized by very low CL relative to carbonates or feldspars (Marshall 1988).

Electron probe micro-analyzer (EPMA)
Quantitative wavelength-dispersive spectrometry (WDS) was performed on a JEOL Superprobe JXA 8230, Hyperprobe JXA 8500F and Hyperprobe JXA-8530Fplus. The samples were coated with a 20-nm-thick carbon film and quantified for Mg, Ca, Fe, Mn, and Sr (and additionally for Ba and Zn in some samples) in carbonates using an acceleration voltage of 15 kV, a beam current of 10 nA, a probe size of 5-40 μm depending on the grain size and relatively short measurement times between 10 to 50 s (for both backgrounds (±) and peak). Those analytical parameters are reported as optimal for precise electron probe micro analyses of carbonates (Zhang et al. 2019). The φ(ρZ) matrix correction scheme (CITZAF; Armstrong 1995) was applied and natural standards were analysed during the course of the analytical sessions to ensure the quality of measured data. Standards were calcite (Ca), dolomite (Ca, Mg), siderite (Fe), strontianite (Sr), rhodonite (Mn), BaSi2O5 for Ba and ZnS for Zn. Under these conditions, analytical errors (2δ) and detection limits were 1.5% and 300 ppm for Ca, 0.4% and 300 ppm for Mg, Mn, and Fe, 0.1% and 200 ppm for Zn and Sr, and 0.1% and 350 ppm for Ba.  Rieger et al. 2020bRieger et al. , 2021a, whole-rock data from mineralized Urquhart Shale (orange; George Fisher; Rieger et al. 2020bRieger et al. , 2021a, and whole-rock data from massive sulfides (dark red; George Fisher; Rieger et al. 2020bRieger et al. , 2021a Mass spectroscopy (LA-ICP-MS) Laser ablation ICP-MS of 19 polished thin sections was carried out using the Analyte Excite 193 nm ArF* excimer-based laser ablation (LA) system (Teledyne Photon Machines, Bozeman, MT, USA), coupled to the quadrupole-ICP-MS iCAP from Thermo Scientific. The LA-system is equipped with a HelEx II 2-volume ablation cell. Helium was used as a carrier gas for aerosol transport from the sample surface to the ICP and was mixed downstream with Ar as a make-up gas before entering the plasma. Operational parameters of the ICP-MS instrument and LA-unit were tuned for maximum sensitivity, low oxide formation based on the 232 Th 16 O/ 232 Th ratio and low laser-induced elemental fractionation based on the 232 U/ 232 Th ratio using NIST SRM 610. We used 43 Ca as internal standard and the certified reference material NIST610 for calibration for all elements. Samples were ablated with a spot size of 50 μm, for 30 s with a repetition rate of 10Hz and an energy density of 2-3 J/cm 2 . Time intervals for data reduction were selected by visual inspection of each spectrum using Iolite™ (Paton et al. 2011) and the data reduction scheme X_trace_elemets_IS (Woodhead et al. 2007). Uncertainty estimates for the elements measured are based on repeated measurement of the reference material MACS-3 and ECRM-752 and are in general below 10 % (relative percent). Data were then screened for contamination by other mineral phases than carbonates (e.g., silicates, phosphates, sulfides) by using the elements Al, K, P, S, Si, and Ti. Data with high concentrations of these elements (Al + K + P + S + Si + Ti > 1000 ppm) were then rejected.

La Ce Pr Nd Sm Eu Gd Tb Dy Y Ho Er
The Ce/Ce * and Eu/Eu * values were calculated geometrically (cf., Lawrence et al. 2006): Lanthanum was not used to calculate Ce * in order to avoid the influence of La anomalies on Ce/Ce * values (cf., Bau and Dulski 1996).

Carbonate in the Urquhart Shale Formation
Detrital carbonate grains are most abundant in siltstone lithologies and preserve similar grain size (≤ 20 μm) and morphology (sub-rounded to angular) to detrital feldspar and quartz grains. Micritic calcite and dolomite is common throughout all lithologies and is very fine-grained (< 5 μm). These detrital and micritic carbonates were too fine grained for LA-ICP-MS analyses and will not be discussed further. Carbonate nodules are interbedded with laminated, pyritic, carbonaceous siltstones. Nodules are, typically, several cm long and wide in the lateral and vertical direction (Fig. 4a, b). Calcite is the dominant carbonate phase in the carbonate nodules and calcite grains are anhedral to subhedral and >10 μm. Coarse-grained (>100 μm) calcite veins cross-cut all lithologies and often contain euhedral pyrite. Where the veins cross-cut carbonate nodules, euhedral pyrite can occur in the nodules.

Carbonate at the George Fisher deposit
Both calcite and dolomite have been observed in the paragenesis at the George Fisher deposit. When calcite and dolomite occur together, calcite can be differentiated by a brighter luminescence, whereas dolomite commonly has a red, lowintensity luminescence signal (Fig. 4i). Nodular carbonates at George Fisher have similar morphologies to those from Shovel Flats samples, but comprise mostly anhedral to subhedral finer-grained dolomite and minor amounts of coarser-grained calcite. The edges of carbonate nodules are also commonly replaced by sphalerite (Fig. 4h).
In the sulfide paragenesis described by Rieger et al. (2020a, b); i.e., ore stages 1, 2, and 3), there are multiple phases of hydrothermal carbonate. The most abundant type of this carbonate is anhedral, coarse-grained (>100 μm), and Fe-rich dolomite; it is generally observed with anhedral quartz, and is in textural disequilibrium with sphalerite and galena, which is indicated by irregular grain boundaries and partial replacement of the dolomite by the sulfides (e.g., Fig. 5h, i). Dolomite, therefore, likely formed before sulfide precipitation of the individual ore stages.
In ore stage 1, this Fe-rich dolomite occurs in stratabound veins together with massive sulfide (e.g., Fig. 5c, d). In the later ore stages 2 and 3, this dolomite occurs as irregular clasts in massive sulfide breccias (e.g., Fig. 5h, i). The other, less abundant, types of carbonate that are associated with massive sulfides are calcites and dolomites that have planar crystal faces and are in textural equilibrium with quartz and sphalerite (Fig. 6d, e). In the samples studied here, these calcite and dolomite variants have only been found associated with ore stage 2. They are coarsegrained (> 100 μm) and occur in cross-cutting, mineralized carbonate-quartz veins or in carbonate-quartz infill adjacent to massive sulfide (Fig. 6). To simplify the later discussion, the carbonate variants that are associated with the massive sulfides will be grouped into two types: hydrothermal dolomite type A, which comprises dolomite in textural disequilibrium with sulfide minerals; and hydrothermal calcite and dolomite type B, which comprises carbonate in textural equilibrium with sulfide minerals. In the ore stage samples no phosphate-bearing mineral phases were observed.
There are also barren veins with coarse-grained calcite (> 100 μm) at George Fisher that cross-cut massive sphalerite of ore stage 1. In the samples studied here, the paragenetic relationship of these barren calcite veins with ore stage 2 and 3 could not be observed.

Rare-earth elements and yttrium (REE+Y)-LA-ICP-MS
The LA-ICP-MS carbonate data of all samples is presented in Rieger et al. (2021b). Total rare-earth element concentrations ( REE) of carbonates are generally within one order of magnitude of the post-Archean Australian shale (PAAS; Fig. 8, Fig. 9c), and several times higher than REE concentrations in  Fig. 2). Generally, the nodular calcite samples from Shovel Flats have REE+Y CN profiles similar to PAAS (Fig. 8a), whereas nodular calcite and dolomite at George Fisher have more variable REE+Y CN profiles (Fig. 8c). Hydrothermal dolomite type A from ore stages 1, 2, and 3 has REE+Y CN profiles that are LREE-depleted and preserves negative Eu/Eu * (Fig. 8d). Hydrothermal calcite and dolomite type B also have REE+ Y CN profiles that are depleted in LREE, but preserve positive Eu/Eu * (Fig. 8e, f). The highest Eu/Eu * values are preserved in hydrothermal carbonate type B (calcite and dolomite; Fig.  8e, f, Fig. 9e) and in calcite veins from the Shovel Flats drill hole (Fig. 8b, Fig. 9e). The Ce/Ce * values are distributed around unity and are slightly skewed to negative anomalies (Fig. 9b). The Y/Ho values are variable, but are generally distributed between 15 and 50, Pr/Yb SN values are typically below 1 in George Fisher dolomite and calcite, and the greatest variability in both Y/Ho and Pr/Yb SN values is preserved in nodular dolomite and calcite from George Fisher samples (Fig. 9d, e).
The in situ composition of nodular calcite from Shovel Flats samples preserves overall higher REE+Y concentrations than whole-rock samples, whereas carbonates at George Fisher are generally LREE-depleted (Fig. 10a). When the data are normalized to the whole-rock REE+Y value of the respective sample (Fig. 10b), the in situ carbonate data from George Fisher is LREE-depleted, whereas Shovel Flats nodular calcite REE+Y profiles are relatively flat. Based on whole-rock REE+Y concentrations (REE+Y wr in ppm) and total carbonate contents of respective samples (XRD carb-wr in %; Rieger et al. 2021a), the mass balance of the REE+Y in carbonate (REE+Y carb-wr in %) has been calculated from the in situ REE+Y concentrations (REE+Y LA-ICP-MS in ppm): Assuming that the analyzed carbonate is representative for the entire sample, carbonate contains 61 % (average value, n = 2 samples) and 49 % (average value, n = 8 samples) of total REE+Y concentrations in the samples from Shovel Flats and George Fisher respectively.

Pre-ore carbonate
The unmineralized Paleoproterozoic Urquhart Shale Formation consists of carbonate-rich siltstones and mudstones that comprise detrital, micritic, and nodular carbonates (e.g., Neudert 1983;Chapman 1999;Painter et al. 1999). The detrital and micritic carbonate is fine grained (≤ 20 μm) and not amenable to laser ablation studies, but the nodular carbonates are coarser grained and occur both in the barren Urquhart Shale at Shovel Flats and in the George Fisher deposit. In the mineralized units the nodular carbonates pre-date the first phase of mineralization (Fig. 4;Chapman 2004;Rieger et al. 2021a), and therefore, provide constraints on carbonate formation under pre-ore conditions.
The nodular carbonates in the Urquhart Shale have been interpreted as pseudomorphs after sulfate evaporites (e.g.,  (2004), Chapman (1999), and Valenta (1988). Endmember compositions of calcite, dolomite, ankerite, and siderite are shown in grey Fig. 8 a- (McLennan 1989), and grey lines are median values of REE+Y whole-rock data of the respective samples that in situ data was generated from (whole-rock data from Rieger et al. 2020bRieger et al. , 2021a Painter et al. 1999) or as diagenetic precipitates (e.g., Domagala et al. 2000). Indicators for shallow-water, evaporitic conditions, such as halite casts and stromatolites, are preserved in other formations of the Mount Isa Group (McClay and Carlile 1978;Neudert and Russell 1981;Neudert 1983). However, sulfur isotope values of early diagenetic pyrite (cf., Painter et al. 1999;Rieger et al. 2020a) indicate that open-system conditions with respect to sulfate availability were dominant during deposition and early diagenesis of the Urquhart Shale. Such conditions are inconsistent with the closed-system conditions typical of evaporite formation and more typical of diagenetic carbonate precipitation (cf., Domagala et al. 2000;Chapman 2004). The composition of the nodular carbonates from the Urquhart Shale Formation is, therefore, considered representative of diagenetic pore fluids, which predated the hydrothermal system at George Fisher.

La Ce Pr Nd Sm Eu Gd Tb Dy Y Ho Er
The REE+Y systematics in seawater and diagenetic pore fluids are typically controlled by redox processes (for a review see Smrzka et al. 2019). Unlike the modern oxygenated oceans, mid-and Paleoproterozoic oceans were mostly ferruginous (anoxic, non-sulfidic), with euxinic conditions spatially and temporally restricted to highly productive margins and with oxygenated conditions only developed in surface waters (e.g., Canfield 1998;Planavsky et al. 2011;Lyons et al. 2014). Indeed, the REE+Y compositions of Proterozoic carbonates are consistent with REE+Y systematics observed in modern anoxic basins (cf., Bau et al. 1997;Tang et al. 2016;Bellefroid et al. 2019). Such REE+Y systematics are characterized by both variable Ce/Ce * SN and variable Eu/Eu * CN values around unity, variable chondritic to super-chondritic Y/Ho values, and REE+Y CN profiles that are similar to PAAS (e.g., Planavsky et al. 2010;Tang et al. 2016;Bellefroid et al. 2019). Compared to oxygenated seawater (Fig. 2a), these contrasting REE+Y signatures are likely the result of (1) the dissolution of LREE-, Ce-, and Ho-enriched phases, such as Mnoxides, at the chemocline between oxidizing and reducing seawater, and (2) enhanced solubility of Eu 2+ under reducing conditions (Sverjensky 1984;Planavsky et al. 2010).
Altogether, the nodular calcite preserves Ce/Ce * SN , Eu/ Eu * CN , and Y/Ho values (Fig. 9b, e) that are consistent with pore waters derived from anoxic bottom waters during the deposition of the Urquhart Shale Formation. The existence of ferruginous conditions is supported by the sulfur isotope values of early diagenetic pyrite and low Mo concentrations from the Urquhart Shale Formation (Rieger et al. 2020a(Rieger et al. , 2021a and in agreement with growing evidence for widespread ferruginous conditions in Proterozoic oceans (e.g., Poulton et al. 2010;Planavsky et al. 2011;Song et al. 2017). As such, the nodular carbonate chemistry from the unmineralized Urquhart Shale Formation preserves normal background marine and early-diagenetic signatures for the Proterozoic.

Carbonate at George Fisher
The George Fisher deposit was formed from a multi-stage mineralizing system (Chapman 2004;Rieger et al. 2020a), and it is possible that multiple processes contributed to the trace element composition of George Fisher carbonates. The deposit comprises a number of different carbonate generations (this study ;Chapman 1999;Murphy 2004), including evidence of hydrothermal alteration that has resulted in dolomitization of pre-ore calcite (Rieger et al. 2021a). Hydrothermal dolomitization is clearly observed in the comparison between nodular carbonates from unmineralized (Fig. 4a-e; Shovel Flats) and mineralized (Fig. 4f-j; George Fisher) Urquhart Shale samples, whereby samples from the latter are partially replaced by sphalerite and have higher dolomite/calcite ratios (Fig. 4g, h;Rieger et al. 2021a).
Nodular carbonates from the unmineralized rocks also have homogenous CL signals (Fig. 4d), whereas those at George Fisher have irregular signals with relatively bright CL of calcite compared to dull and quenched CL of dolomite (Fig. 4i). Such dull and quenched luminescence is characteristic of carbonate minerals with high Fe and variable Mn concentrations (Pierson 1981;Machel and Burton 1991;Baele et al. 2019), which is also consistent with the major-element chemistry of George Fisher and Shovel Flats carbonates (Fig. 7, Fig. 9a). These irregular CL signals (Fig. 4i) and variable Fe and Mn concentrations (Fig. 9a), together with higher whole-rock dolomite/calcite ratios (Rieger et al. 2021a), are consistent with dolomitization of pre-ore calcite that was associated with the mineralization at George Fisher. Similar to calcite and dolomite in nodular carbonate samples at George Fisher, type A dolomite has dull to quenched luminescence (Fig. 5c, h). It is, therefore, likely that dolomitization was associated with all stages of mineralization at George Fisher.
This dolomitization and hydrothermal alteration has not resulted in fractionation of the whole-rock REE+Y CN composition between Shovel Flats and George Fisher samples, apart from dilution of the REE+Y in the samples containing massive sulfide mineralization ( Fig. 3; Rieger et al. 2021a). In contrast, the REE+Y CN composition of hydrothermal carbonate from George Fisher is fractionated (LREE-depleted) relative to respective whole-rock and in situ pre-ore carbonate REE+Y CN profiles (Fig. 10). One possibility is that REE+Y were fractionated during precipitation of another phase that formed in an assemblage with the hydrothermal carbonate. For example, phosphate minerals (e.g., monazite or apatite) have a high partition coefficient for the REE+Y and could have resulted in REE+Yfractionated carbonate without changing the overall wholerock composition (cf., Debruyne et al. 2016).
A chemical mass balance analysis of the whole-rock composition of samples from Shovel Flats and George Fisher indicates that P was relatively immobile and in fact was slightly depleted during hydrothermal alteration at George Fisher (Rieger et al. 2021a). Co-precipitation of a P-bearing phase would have only locally affected the REE+Y composition of carbonate at George Fisher (e.g., due to local dissolution and precipitation reactions of P-bearing phases). Nevertheless, there is no evidence of hydrothermal phosphate mineral phases in an bright CL (Machel and Burton, 1991) dull CL (Machel and Burton, 1991) quenched CL (Machel and Burton, 1991) no CL extinction (Pierson, 1981) luminescence (Pierson, 1981) Bau et al. 1997), the median value of the Mount Isa Group (5; Nance and Taylor 1976), and the whole-rock median values of Urquhart Shale from the unmineralized Shovel Flats drill core (6), from the George Fisher deposit (7), and from massive sulfide samples from George Fisher (8; Rieger et al. 2020bRieger et al. , 2021a. Note that values of #3, and #3 and #4 are lower than the plotted x-axis for b and c assemblage with the ore stage sulfides at George Fisher, making this scenario unlikely. Rather, it is more likely that hydrothermal carbonate only accounts for a small fraction of the whole-rock REE+Y budget at George Fisher. For example, it is possible that the whole-rock REE+Y composition is dominated by pre-ore carbonate, which has a similar REE+Y CN profile to the whole-rock (Figs. 4, 8, 10), and by phyllosilicate mineral phases, which are commonly important hosts for REE+ Y in fine-grained sedimentary rocks (e.g., Condie 1991;Abbott et al. 2019). Based on these considerations, the following discussion will focus on hydrothermal carbonate at George Fisher and the implications for the hydrothermal system that can be derived from the carbonate REE+Y composition.
The REE+Y composition of hydrothermal dolomite should provide a record of fluid chemistry that includes both pre-and syn-sulfide hydrothermal events. For example, the sulfide replacement textures and irregular grain boundaries with sphalerite ( Fig. 4g-i, Fig. 5c, d, h, i) are indicative of a pre-sphalerite timing for some dolomitization (i.e., type A) at the onset of mineralization (e.g., Chapman 1999). In contrast, the dolomite and calcite from type B formed in textural equilibrium with sphalerite (Fig. 6d, e). On the basis of these paragenetic relationships, the REE+Y chemistry of the pre-sphalerite dolomite (type A) should provide constraints on the hydrothermal fluids that caused the alteration of pre-ore calcite whereas type B carbonate should provide constraints on the hydrothermal fluids during sphalerite formation. In order to constrain the syn-sphalerite end-member of the hydrothermal fluid, the REE+Y composition of this type B carbonate is discussed first.
Type B calcite and dolomite are characterized by (1) LREE depletion relative to PAAS, respective whole-rock, and preore calcite, (2) chondritic and sub-chondritic Y/Ho values, and (3) Eu/Eu * values >1 (Fig. 6f, i, Fig. 8e, f, Fig. 9d, e). Lightrare-earth-element fractionation is a common feature of hydrothermal fluids and there are a number of processes that could have produced the REE+Y CN profiles in the type B carbonates (Fig. 6f, i, Fig. 8e, f). A mineralogical control on this LREE fractionation is unlikely because both calcite and dolomite are LREE-depleted relative to pre-ore calcite (Fig. 8a, e, f). This observation is consistent with other studies on hydrothermal gangue carbonates that indicate no significant REE+Y fractionation between calcite and dolomite (e.g., Hecht et al. 1999;Roberts et al. 2009). The depletion of LREEs in carbonate mineral phases can also be inherited from complexationcontrolled REE+Y-leaching from crustal source rocks (e.g., Lüders et al. 1993) or from LREE-scavenging by mineral phases with high REE-partition coefficients (e.g., monazite) along the fluid pathways (e.g., Debruyne et al. 2013Debruyne et al. , 2016. The source rocks for the mineralization are unknown at George Fisher and whether the footwall units could contain significant amounts of phosphatic minerals is also unknown. However, if the REE+Y signatures were inherited from the metal source or fluid-rock interaction along the flow pathway, they would likely be more homogenous than the variability preserved in the REE+Y profiles of hydrothermal calcite and dolomite at George Fisher (Fig. 6f, i, Fig. 8e, f) Fig. 10 a Chondrite normalized REE+Y profiles for samples from the Shovel Flats drill core and the George Fisher deposit. The lines denote median values for in situ carbonate data from the Shovel Flats drill core (dark blue) and from the George Fisher deposit (red), and for whole-rock REE+Y profiles from the Shovel Flats drill core (light blue; n = 41) and from the George Fisher deposit (orange; Urquhart Shale and massive sulfide; n = 69); the light purple (Shovel Flats) and light red (George Fisher) fields denote the 25-75 percentiles of the in situ carbonate REE+Y profiles. For more detailed whole-rock data see Fig. 3. b REE+ Y profiles (medians and 25-75 percentiles) for individual carbonate subgroups from the Shovel Flats drill core and the George Fisher deposit normalized to whole-rock REE+Y data of the respective sample that in situ data was generated from Y CN profiles (Fig. 3;Nance and Taylor 1976;Hannan et al. 1993); as a consequence, a hydrothermal fluid buffered by these rocks would most likely not be LREE-depleted. Instead, the observed variability may result from closedsystem Rayleigh fractionation during carbonate precipitation (e.g., Kontak and Jackson 1995), which is consistent with variable LREE concentrations in stage B calcites (e.g., Fig.  6i). Similar observations have been made from experimental data, where carbonate grown from static solutions preserved large variability in REE composition between individual growth zones (Barker and Cox 2011). Fluid composition and temperature also provide a major control on REE+Y fractionation (e.g., Michard 1989;Craddock et al. 2010;Williams-Jones et al. 2012). Unfortunately, there are no direct constraints on the composition of the mineralizing fluids at George Fisher (e.g., fluid inclusions). Based on the basin fill in the stratigraphy below the Carpentaria CD-type deposits, Cooke et al. (2000) suggested that the mineralizing fluids were likely oxygenated, neutral to mildly acidic, and saline brines (ca. 25 wt.% NaCl equiv ), which would have been highly effective in transporting Zn and Pb at moderate temperatures (150°C). Fluid inclusion data from the relatively undeformed Century deposit do indicate relatively low-temperature (<125°C) and high-salinity fluids (ca. 23 wt.% NaCl equiv ; Polito et al. 2006). Notably, salinity can have a strong effect on REE+Y solubility and fractionation. For example, experimental studies and numerical modelling of REE+Y partitioning and solubility in hydrothermal fluids suggest that chlorides are the dominant ligands for REE+Y-complexation at hydrothermal temperatures and that Cl --ligands form the strongest complexes with the REE+Y with the lowest charge/radius ratios (Eu 2+ and LREE; e.g., Migdisov et al. 2009;Williams-Jones et al. 2012;Perry and Gysi 2018). If the mineralizing fluids for the George Fisher deposit were saline (cf., Cooke et al. 2000;Polito et al. 2006), Cl --complexes would have been the dominant ligands for REE+Ycomplexation and could have resulted in stronger LREEretention in the fluid during carbonate precipitation. Consequently, REE+Y CN profiles of syn-ore carbonates should generally be LREE-depleted despite the effects of intrinsic processes.
A second common signature of type B carbonates at George Fisher are chondritic and sub-chondritic Y/ Ho values ratios (Fig. 9d, e). Chondritic Y/Ho values are typical of hydrothermal fluids (Bau and Dulski 1999;Douville et al. 1999), whereas sub-chondritic Y/Ho values may indicate fractionation of Y relative to Ho. Experiments at ambient temperatures suggest that Ho partitions more strongly into calcite than Y (e.g., Tanaka and Kawabe 2006). Such differential partitioning is consistent with sub-chondritic Y/Ho values in type B calcite at George Fisher. There is, however, no experimental data available for Y/Ho fractionation during carbonate precipitation at hydrothermal temperatures. If future experiments show differential fractionation of Y and Ho into carbonate minerals at hydrothermal temperatures, subchondritic Y/Ho values may be a distinguishing criterium for hydrothermal carbonates.
A third dominant feature of type B carbonates at George Fisher are Eu/Eu * values >1 (Fig. 6f, i, Fig. 8e, f, Fig. 9b). The major control on Eu/Eu * is temperature ( Fig. 2b; Sverjensky 1984;Bau 1991;Bilal 1991). In oxidized hydrothermal fluids, such as those suggested for the Carpentaria CD-type deposits (SO 4 2-> H 2 S; Cooke et al. 2000), Eu/Eu * values >1 can only develop at high temperatures, because the more soluble Eu 2+ is only stable at temperatures >200-250°C ( Fig. 2b; Sverjensky 1984;Bau 1991;Bilal 1991). It is, therefore, unlikely that the George Fisher deposit formed from fluids that were initially cooler than 200-250°C. Divalent Eu is, however, not readily incorporated into the crystal lattice of carbonate minerals (Bau and Möller 1992); therefore, the hot, Eu 2+dominated hydrothermal fluids (>200-250°C) must have cooled in order to stabilize sufficient Eu 3+ to produce Eu/ Eu * values >1 in the type B carbonate at George Fisher. Such fluid cooling at George Fisher is consistent with thermodynamic models for sulfide formation (Rieger et al. 2020a), calcite twinning temperatures (Murphy 2004), and also with overall temperature constraints for the Urquhart Shale Formation (see summary for temperature constraints in Rieger et al. 2021a). Furthermore, this is consistent with Eu/ Eu * values <1 preserved by type A dolomite from ore stage 2 (Fig. 5j). These Eu/Eu * values <1 in type A dolomites support the hypothesis that the hydrothermal fluids were hot at the onset of mineralization, and then cooled below 200-250°C before type B calcite and dolomite formation.
Dolomite that formed from dolomitization of pre-ore calcite at the onset of mineralization is preserved both in samples of mineralized nodular carbonate and in massive sulfide samples at George Fisher ( Fig. 4f-i, Fig. 5a-i). Partially mineralized nodular carbonates preserve the highest textural variability, together with the highest variability in the REE+Y composition (Fig. 4e, j, Fig. 8c, Fig. 9b-d). Therefore, they likely comprise mixed REE+Y compositions from pre-ore diagenetic and hydrothermal end members, which is consistent with the large ranges of LREE-depletion relative to pre-ore calcite and the large variability of Y/Ho values (Fig. 8c, Fig. 9d). The REE+Y composition of type A dolomite at George Fisher is slightly more homogenous than that of mineralized nodular carbonate (Fig. 8d, Fig. 9d), and this may indicate that the REE+Y composition is more dominated by hydrothermal signatures as a result of progressive dolomitization. The REE+ Y CN profiles of type A dolomite are LREE-depleted relative to PAAS, pre-ore calcite, and to whole-rock REE+Y profiles (Fig. 4j, Fig. 5e, j, Fig. 8d, Fig. 10b). This LREE-depletion is similar to REE+Y CN profiles of type B calcite and dolomite and consistent with Cl --complexation and LREE-retention in the fluid.
Chemical alteration and dissolution-reprecipitation of carbonate is not limited to the high-grade domains of ore deposits, and chemical and isotopic changes in carbonate can be traced beyond zones of visible alteration (e.g., Barker et al. 2013;Vaughan et al. 2016). In fact, isotopic studies indicate that large-scale 18 O-depletion in carbonate is associated with the Mount Isa Cu-system (Waring 1990;Waring et al. 1998a), and similarly, 18 O-depletion has also been reported for carbonate from the George Fisher Zn-Pb system (Chapman 1999). There is, therefore, the potential that the mineralizing fluids at George Fisher have also produced alteration of the REE+Y signatures beyond the highly mineralized zones. For example, sulfide formation would have resulted in the destabilization of chloride complexes and carbonate dissolution via the generation of acid (e.g., ZnCl 2 + HS -→ ZnS + H + + 2Cl -). Such a fluid could have caused even higher solubilities of LREE as Cl --complexes, which could have resulted in LREE-depletion along the fluid migration pathway beyond visible alteration.
Late calcite veins postdate the stratabound Zn mineralization at George Fisher, and late calcite veins cross-cut all lithologies in the unmineralized Shovel Flats drill core. Cross-cutting calcite veins have also been described before at Mount Isa and George Fisher (e.g., Waring et al. 1998b;Murphy 2004). The calcite veins at George Fisher preserve REE+Y signatures similar to the pre-ore calcites and to PAAS (Fig. 8g, Fig. 9), which may indicate that they are rock-buffered. The late calcite veins from Shovel Flats preserve REE+Y signatures that are more similar to type B calcites (Fig. 8b, Fig. 9), which could represent a common hydrothermal affinity. The differences in REE+ Y composition suggest that there are, at least, two generations of late calcite veins. These veins were, however, not the focus of this study and future work is needed to understand their significance and relationship to mineralization.
Overall, the Y/Ho values, Eu/Eu * values, and LREEdepleted REE+Y CN profiles that are characteristic of George Fisher calcite and dolomite are consistent with interaction of a hot, saline (Cl --rich) hydrothermal fluid that subsequently cooled to temperatures below 200-250°C during fluid-rock interaction. Notably, such conditions correspond well with other CD-type systems, for which the thermal evolution of the hydrothermal fluids is well constrained (e.g., MacMillan Pass; Magnall et al. 2016).

Implications for the application of in situ REE+Y carbonate data
There are some important implications of the in situ REE+Y carbonate data from George Fisher, both for ore formation models and future exploration programs. The George Fisher deposit formed from a multi-stage system and all mineralization post-dated deposition of the host rock (Chapman 2004;Rieger et al. 2020a). The mineralizing fluids must have reacted with the host rock, therefore, in order to create sufficient porosity for fluid flow and base metal sulfide precipitation. The fractionation between whole-rock REE+Y profiles and carbonate in situ REE+Y profiles (Fig. 10a, b) indicates that (1) the whole-rock REE+Y composition is still dominated by diagenetic carbonate and other REE+Y-bearing phases after hydrothermal alteration, that (2) carbonate mineral phases were more sensitive to REE+Y alteration relative to the whole-rock, and that (3) fluid-rock ratios were not high enough to alter whole-rock REE+Y compositions (cf., Michard and Albarède 1986;Bau and Möller 1992). Altogether, this suggests that the hydrothermal fluids at George Fisher were highly selective in replacing and altering pre-ore carbonate mineral phases. Moreover, this means that REE+Y carbonate data could provide a more sensitive record of hydrothermal alteration relative to bulk-rock techniques. At this stage, however, the small number of samples analyzed at George Fisher is not sufficient to evaluate the 3D REE+Y footprint of the mineralizing system.
The paragenetic and geochemical data provide evidence of fluid-rock interaction that involved hot, acidic, saline, and base metal-bearing hydrothermal fluids. The resulting replacement and dolomitization of pre-ore calcite led to the development of secondary porosity, fluid cooling, and pH increase, which lowered base metal solubilities and caused sulfide precipitation. Such a diagenetic carbonate replacement model is consistent with (1) carbonate replacement textures (e.g., Fig. 4 ;Chapman 2004;Rieger et al. 2021a), with (2) fluid cooling indicated by Eu/Eu * <1 in type A dolomite opposed to Eu/Eu * >1 in type B dolomite and calcite ( Fig. 8d-f, Fig. 9e), and with (3) Cl --rich hydrothermal fluids indicated by LREE-depleted REE+Y CN profiles in George Fisher carbonates (Fig. 4j, Fig. 5e, j, Fig.  6f, i, Fig. 8c-f, Fig. 10). Notably, carbonate replacement was reported for a number of other CD-type massive sulfide systems in the Carpentaria province (Eldridge et al. 1993;Perkins and Bell 1998;Chapman 2004;Magnall et al. 2021;Spinks et al. 2021), which highlights the importance of the reactive carbonate-rich host rocks throughout the province.
More broadly, it is generally accepted that base-metal deposits in sedimentary basins form from saline hydrothermal brines (Bodnar et al. 2014;Heinrich and Candela 2014); and indeed, previous studies on Zn-Pb-and Cu-systems in the Carpentaria province propose that saline fluids were responsible for base-metal transport (Heinrich et al. 1989;Cooke et al. 2000;Polito et al. 2006). At George Fisher, the LREEdepleted REE+Y CN profiles in calcite and dolomite (Fig. 4j,  Fig. 5e, j, Fig. 6f, i, Fig. 8c-f, Fig. 10a, b) are consistent with LREE-retention in such Cl --rich saline hydrothermal fluids. Furthermore, there is evidence that LREE-depletion is a common feature of carbonate mineral phases from a number of mineral deposits that formed from saline fluids (Roberts et al. 2009;Debruyne et al. 2013;Genna et al. 2014;Magnall et al. 2016), whereas LREE depletion is not evident from carbonates that precipitated from low-salinity fluids that were dominated by HS --or CO 3 2--complexes (Maskenskaya et al. 2015;Vaughan et al. 2016). As a consequence, LREE depletion in calcites or dolomites relative to background REE+Y profiles may be an effective tracer for fluid-rock interaction involving saline hydrothermal fluids. So, if this hypothesis is supported by future studies, the analysis of REE+Y in carbonate mineral phases could be useful for future exploration programs in sedimentary basins in order to identify hydrothermally altered carbonate-rich lithologies.
In summary we show that LA-ICP-MS analyses of REE+Y in carbonate minerals can distinguish hydrothermal carbonate from pre-existing diagenetic carbonate, and this may be a useful exploration tool in combination with whole-rock and isotope geochemistry for mineral deposits hosted in sedimentary rocks. Future studies should test this hypothesis by systematically sampling the 3D grid around an (undeformed) deposit to assess the nature of the distal footprint on carbonate REE+Y systematics.