Abstract

During the Albian, the hyperextension of the Pyrenean passive margin led to a hyperthinning of the continental crust and the subsequent subcontinental mantle exhumation. The giant Trimouns talc-chlorite deposit represents the most prominent occurrence of Albian metasomatism in the Pyrenees, with the occurrence of the largest talc deposit worldwide. Consequently, this deposit, which is located on a fault zone and a lithological contact, represents one of the major drains at the scale of the Pyrenees and one of the best geological targets in order to determine the origin(s) of the fluid(s) that circulated during this period. Talc-chlorite ore is characterized by the presence of brines trapped in dolomite, quartz, and calcite fluid inclusions in the vicinity of the talc-rich zone. Considered as being responsible for the formation of talc, these fluids may be interpreted in several ways: (i) primary brines expelled from Triassic evaporites, (ii) secondary brines produced through halite leaching by diagenetic/metamorphic fluids, and (iii) brines derived from seawater serpentinization of mantle rocks. Stable isotope analyses (δ13C, δ18O, δD, and δ37Cl) and Cl/Br ratio measurements in fluid inclusions and their host minerals were carried out in order to determine the origin of the fluid(s) involved in the formation of the ore deposit. The data are consistent with a primary brine origin for the mineralizing fluid, which could have been expelled from the Triassic levels. Other hypotheses have been tested, for example, the production of brines via the seawater concentration during serpentinization. The geochemical proxies used in this study provide equivocal results. The first hypothesis is by far the most realistic one considering the numerous occurrences of Trias formations nearby, their deformation during the extension, and the drainage of the expulsed brines as evidenced by the high-salinity fluid inclusions found all around the deposit. Alternatively, the exhumation of the mantle is considered as a major source of heat and stress that favored brine migration along the major shear zones. Our results fit well with brine circulation in a hyperextensional geodynamic context, which is related to the formation of the talc-chlorite ore, the thinning of the continental crust, and the exhumation of the subcontinental mantle, in accordance with recent works.

1. Introduction

The circulation of fluids is known to play a crucial role both in the deformation of rocks and the transfer of heat and elements throughout the lithosphere. In the Pyrenees, there is a lot of evidence of widespread metasomatism attesting to large-scale fluid circulation in the crust. Three main types of metasomatism are recognized: (i) the albitization of the plutonic basement [1], (ii) the serpentinization of peridotite bodies ([2, 3]; Saint Blanquat et al., 2016), and (iii) the magnesian metasomatism expressed by the dolomitization, talcification, and chloritization of the basement ([48]; Core et al., 2017; [9]). All of these fluid/rock interaction systems were active between at least 122 Ma and 96 Ma [2, 1014] and reached a climax during the Albo-Cenomanian regional hyperextension of the Pyrenean passive margin which led to mantle exhumation [15] and the high-temperature-low-pressure (HT/LP) metamorphism of the North Pyrenean Zone (NPZ) [11, 1619]. This synchronism suggests strong links between fluid circulation, deformation, and heat transfers at the crustal scale in the Pyrenees as already proposed by some authors [11].

Beyond the effects of fluids on the thermomechanical behavior of the crust, the nature and origin of the different fluids that have circulated during the Albo-Cenomanian hyperextension remain poorly known. The occurrences of serpentinites, ophicalcites, albitites, talcitites, and chloritites all reflect the location of the main fluid pathways that were active during the Albo-Cenomanian period. Consequently, these evidences of metasomatism are excellent geological targets to characterize the fluid circulation systems active during this period.

This study is focused on the Trimouns talc-chlorite deposit (Figure 1) which represents the most prominent occurrence of Albo-Cenomanian magnesian metasomatism in the Pyrenees and one of the biggest talc deposits worldwide. Roughly 400000 tons of talc are produced annually, and close to around 20 million tons have already been extracted from this deposit. The current resources are estimated to be approximately 20 million tons [20]. Therefore, the large amount of fluid needed to form this talc deposit suggests that the Trimouns deposit represents one of the major drains at the scale of the Pyrenees.

In this paper, we characterize the fluid circulation system at the origin of the magnesian metasomatic system expressed at Trimouns by the widespread formation of dolostone, talc, and chlorite. Stable isotope analyses (δ13C, δ18O, δD, and δ37Cl), the Cl/Br ratio, and chlorinity measurements taken on fluid inclusions (FI) and their host minerals were used to document the paleofluid features and their possible origin. The objective is to relate the fluid geochemical signatures to the major fluid reservoirs and to constrain the potential migration pathways at the crustal scale that favored subsequent fluid-rock interaction at the origin of the deposit.

2. Geological Setting

The Pyrenees fold and thrust belt (Figure 1) result from two successive orogenies, the Variscan orogeny at the end of the Paleozoic, characterized by a polyphase compression history, and the Mesozoic-Cenozoic Alpine orogeny. The latter results from the inversion of a transcurrent hyperextended rift linked to the opening of the Bay of Biscay during the Cretaceous period [2124]. This extreme thinning of the continental crust led to the local exhumation of the subcontinental mantle at the base of the Albo-Cenomanian basins [19, 25], occurring today as relicts of serpentinized peridotites in the North Pyrenean Zone [15]. The Albo-Cenomanian basins were formed during the opening of the Bay of Biscay which induced transtensive sinistral movements in the NPZ [22]. The basin formation was associated with a HT-LP metamorphism between ca. 113 and 85 Ma [1618, 2628]. To date, the consensus is to attribute this thermal anomaly to very high geothermal gradients related to the intense crustal and lithospheric thinning [11, 1619].

In the NPZ, the circulation of large amounts of fluid under hydrothermal conditions is evidenced by the presence of numerous occurrences of talc-chlorite [68, 10], albitite [1, 12, 13], and serpentinized peridotite ([2]; Saint Blanquat et al., 2016). The Na-Ca metasomatism (albitization) is well dated between 110 and 92 Ma [12, 13] whereas the Mg metasomatism depicts widespread ages throughout the Pyrenees. The main talc formation episode is Cretaceous in age, dated between 122 and 96 Ma [10, 14] based on the U-Pb dating of REE-bearing minerals, particularly allanite formed in dolostone geodes [29, 30]. Lastly, the serpentinization event is supposed to be Albo-Cenomanian in age, even if it has still not been clearly dated yet, given the consistency of this age with the hyperthinning model for the crust and associated mantle exhumation [2].

The Trimouns deposit is localized at the eastern edge of the Saint Barthelemy Massif (SBM) (Figure 1). From the base to the top, the SBM is composed of (i) amphibolitic to granulitic gneiss, (ii) migmatites and mica schists, and (iii) middle to upper Paleozoic series including Ordovician to Devonian schists, quartzite, and Silurian black schists and marbles [3133]. Dolostones occur in the Paleozoic succession, but their origin remains a matter of debate. Some authors suggest that they are Ordovician [31] whereas others suggest that they were formed during Albo-Cenomanian times, via the metasomatism of the Silurian marbles [6]. The contact between the mica schists and carbonate rocks and the Silurian black schists is tectonic and corresponds to a fault considered to be a decollement level located within the Silurian black schists [10]. Talc and chlorite occur as lenses that are discontinuously distributed along this tectonic interface where mineralizing fluid(s) has/have been channelized. In further detail (Figure 2), chlorite ore is localized above the footwall of the fault where mica schists and, locally, pegmatites occur. Talc ore is mainly localized at the hanging wall of the fault composed of dolostones, marbles, and schists, from the Siluro-Devonian series (Figure 1). This spatial distribution of the two ores results from different interlinked chemical reactions driven by the nature of the metasomatized parent rock.

3. Constraints on the Fluid Circulation System

3.1. The Metasomatic Alteration

At Trimouns, talc can be formed following two main reactions: (1)The interaction between the dolostones and a silica-rich fluid [8] (2)A direct fluid oversaturation with respect to the talc, e.g., a fluid enriched in Mg and Si [8]

The Mg(Fe) chlorite formation is dependent from both Mg supply and Al availability, which is known to be poorly mobile. This explains why chlorite ore (mainly clinochlore, [29]) forms mostly at the expense of mica schists and pegmatites, which are rich in aluminous minerals (mainly muscovite and plagioclase) and quartz. According to Moine et al. [8], the aluminum needed to form chlorite primarily comes from the dissolution of muscovite in mica schists following the reaction: whereas in pegmatites, the main source of aluminum is plagioclase [29]. The formation of chlorite from plagioclase can be expressed as follows:

As shown above, a silica supply is needed to form both chlorite and talc. Fortuné (1980) and De Parseval [29] suggested that silica was leached from quartz present in the hosting rocks of the deposit (pegmatites and mica schists). The dissolution of most silicates, except tourmaline, and their replacement by talc confirm that silicates are all unstable in the presence of mineralizing fluids. The excess silica precipitated around the talc body as discrete euhedral quartz in the dolostone geodes, together with euhedral calcite, dolomite, and talc (Figure 3).

Ca2+ and CO2 leached out during talc formation from dolostones induce the formation of calcite ([7, 29], the so-called hydrothermal calcite from [6]), which occurs in geodes or crosscutting dolostones as veins (Figure 3). Based on the similar REE compositions, δ18O and δ13C values, and 87Sr/86Sr ratios, Boulvais et al. [6] showed that hydrothermal calcite formed from a single fluid event under a variable fluid-rock ratio.

Based on reactions (2), (3), (4), and (5), an additional supply of magnesium is needed to precipitate the talc and to form the chlorite via the alteration of the pegmatites and mica schists. The source of magnesium is still a matter of debate: either the magnesium is provided by the fluid itself (Fortuné, 1980; [10, 29]) or it is acquired by the mineralizing fluid interacting with encountered percolated rock from the source to the deposit.

3.2. The Mineralizing Fluid

Several studies have attempted to characterize the nature of the mineralizing fluid(s) and its physical circulation conditions. Based on a microthermometric study of the fluid inclusions hosted by the euhedral quartz localized in pods and inferred to be synchronous with Mg metasomatism, Moine et al. [8] proposed that the main mineralizing fluid was a brine (~25 wt% eq. NaCl) that circulated at about 300°C and 1.5 kbar. De Parseval et al. [31] came to the same conclusions based on a study of the fluid inclusions hosted by (i) the quartz pods in chloritite ore, (ii) the fluorapatite in the chloritite ore, and (iii) the allanite observed in the dolostone geodes (as shown in Figure 3). In these two studies, the P-T conditions rely solely on the consideration of fluid inclusion isochors and no independent pressure or temperature estimates are available which means that the P-T path is a hypothetical estimation. The main mineralizing fluids were evaporated seawater saturated with respect to halite on the basis of the chlorinities and the Cl/Br ratio (values from 140 to 565) that circulated between 250-280°C and 0.4 to 2 kbar [4, 5, 34]. These brines contain traces of N2- and a low magnesium content [5].

All of these studies conclude that primary brines are probably responsible for the talc ore formation. These brines are supposed to be synchronous with the deposition of the nearby Triassic evaporites. A stable isotope study (δ13C, δ18O, and δD) of the whole metasomatic sequence, from unaltered marbles to dolostones to hydrothermal minerals (talc, chlorite, calcite, and quartz), suggested that both the dolomitization and talc-chlorite ore formation stages resulted from the same continuous hydrothermal event [6]. This implies that the supply of an allochthonous magnesium source was transported by the mineralizing fluid. Even though the authors suggest that the main fluid reservoir was Albian seawater, an (minor) involvement of secondary brines is also suggested. Lastly, Boutin [7] provided an extended thermometric study of the different chlorites occurring at Trimouns in order to discuss the range of temperatures for the formation of the deposit. He concluded that the chlorite in the deposit formed between 275°C and 350°C with a first mode at around 275°C and a second one at around 325°C.

To summarize, from previous studies, it appears that the main fluid(s) at the origin of the talc-chlorite ore deposit probably was/were brine(s) that circulated at a temperature of [4, 5, 7, 8, 31, 34]. However, the origin(s) of the brine(s) has/have not been clearly addressed. It has also been suggested that the whole metasomatic sequence, from dolomitization to hydrothermal mineral formation (talc, chlorite, calcite, and quartz), resulted from the same continuous hydrothermal event [6].

The main objectives of this work are to characterize the type of brines (primary vs. secondary) involved in the formation of the talc-chlorite deposit and to determine if these brines have been mixed with other types of fluids. The spatial extent of the fluid circulation system is specified by analyzing the regional sample. Following the hypothesis of Boulvais et al. [6], the fluids at the origin of dolomitization and talc-chlorite formation are compared in order to determine if the whole metasomatic sequence resulted from the same hydrothermal event. Lastly, the possible sources of fluids involved in the magnesian metasomatic system recorded at Trimouns are discussed in order to constrain their potential migration pathways at the crustal scale.

4. Methods

4.1. Sampling Strategy

Several minerals considered to have been synchronous with the main talc formation stage were sampled: (i) euhedral quartz and hydrothermal calcite, both occurring in geodes hosted by dolostones and (ii) calcite veins crosscutting dolostones (Figure 3). Euhedral quartz and calcite are considered synchronous with the main talc formation since, in geodes, they are spatially and texturally associated with euhedral talc and with REE minerals (mainly allanite). The calcite veins crosscutting the dolostones are either (i) thin fractures (≤2 cm) filled by massive calcite or (ii) thicker veins (≥10 cm) filled by euhedral calcite and considered to be tension gashes. As explained previously, these calcite occurrences are supposed to be by-products of the transformation of dolostone into talc [6, 7, 29] and are consequently considered as hydrothermal calcite [6]. In the following sections, the term “hydrothermal calcite” will refer to the thin and thick calcite veins crosscutting the dolostone sampled for this study. In order to demonstrate the genetic link between the dolomitization stage and the talc-chlorite ore formation stage suggested by Boulvais et al. [6], a series of dolostones was sampled. These samples are massive and undeformed dolostones localized in the hanging wall of the deposit. The intergranular porosity is filled by talc, chlorite, and organic matter. As described above, these dolostones host geodes and calcite veins.

Lastly, the quartz-calcite veins occurring in the Devonian schist located at the southern part of the deposit (Caussou, Figure 1) were sampled in order to characterize the possible spatial extension of the mineralizing fluid circulation. They represent further evidence of fluid circulation in the close vicinity of the deposit; however, no geochronological constraints are available. This sample comes from a geometrically complex network of quartz-calcite veins localized along a main fault in a deformation zone affecting the Devonian schist (Figure 4). This network attests to the circulation of a significant amount of fluid in this deformation zone located just to the south of the North Pyrenean Fault.

4.2. Fluid Inclusions
4.2.1. Salinity Measurements and Absolute Cl Content

Fluid inclusions were studied by microthermometry using a Linkam heating-freezing stage. The fluid inclusions used for the calibration were a pure CO2 natural standard fluid inclusion from Camperio (triple point at -56.6°C) and H2O-NaOH synthetic fluid inclusions (ice melting temperature at -0.4°C). The accuracy at high temperature is ±2°C. The salinities were determined via either the ice melting temperature or the melting temperature of the salt cubes when present.

In dolostones, fluid inclusions are very small and the microthermometric measurements for Tm ice are rather difficult to obtain. Thus, the salinity was obtained by Raman spectroscopy at room temperature using the methodology published by Dubessy et al. [35] and Caumon et al. [36]. The typical model uncertainty is ±0.3 mass% NaCl (after the revision of the calibration curve in [36]).

The average Cl content of the fluid inclusions for the various samples was estimated based on the average salinity value calculated from the whole salinity dataset for each type of fluid inclusion.

4.2.2. Cl/Br Ratio

The bulk Cl and Br contents of the water trapped in the fluid inclusions were determined at the Université de Lorraine, GeoRessources, Vandoeuvre-lès-Nancy, France, using the bulk crush-leach technique described by Banks and Yardley [37], Banks et al. [38], and Gleeson et al. [39]. The quartz and carbonate veins were crushed to a grain size fraction of at least 1–2 mm. The pure mineral separates were then picked and placed in an ultrasonic bath for a maximum of 1 hour. Then, the samples were cleaned by boiling twice in Milli-Q water. The cleaned quartz and carbonate grains were then crushed to a fine powder in an agate pestle and mortar and leached by Milli-Q water.

Leachates from the fluid inclusions were analyzed for Cl and Br with a Dionex DX-500 ion chromatograph. The detection limit for the anions in the fluid inclusion leachates was 3 ppb for Cl and 1 ppb for Br. The results are expressed as a Cl/Br molar ratio.

4.2.3. δ37Cl

Stable isotope analyses were performed on the same bulk leachate as those used to measure the Cl/Br molar ratio as explained in the previous section. The sample preparation technique for the δ37Cl analysis using CH3Cl as a working gas has been described in detail elsewhere [4043]. It consists of (1) the precipitation of chloride as AgCl, (2) the methylation of AgCl to gaseous CH3Cl, and (3) the purification of CH3Cl via separation from the excess CH3I in the two GC columns filled with Porapak Q. The δ37Cl values are then determined on purified CH3Cl gas with a DI-IRMS (Thermo Finnigan, Delta plus XP) at the Laboratoire de Géochimie des Isotopes Stables (Institut de Physique du Globe de Paris and Université Paris 7). δ37Cl is expressed as the permil deviation relative to the standard SMOC (Standard Mean Ocean Chloride). Every seventh sample, a series of three seawater standard isotopic compositions was measured in order to estimate the reproducibility of the analyses. For each session (10 measurements), the δ37Cl values of the samples were corrected by the mean seawater isotopic composition of the session. For this study, at least two measurements were taken from the same leachates during the various analytical sessions. Lastly, the isotopic compositions of the samples are the mean value of the duplicate. The uncertainty associated with each sample corresponds to the standard deviation and is lower than 0.1‰. During the course of this study, the external reproducibility on the seawater standard was (, ), as routinely obtained at IPGP for over two decades now [40, 44, 45].

4.2.4. δD

Fluids trapped in the quartz, hydrothermal calcite, and dolostones were extracted by crushing under vacuum. Grains weighing from 1.0 to 6.7 g and measuring several millimeters to centimeter were loaded into steel tubes and degassed overnight at 120°C under vacuum in order to release any water adsorbed at the mineral surface. H2O was separated cryogenically from other gases and reduced into H2 through a uranium reactor at 800°C. The D/H ratios of the H2 gases were measured using a GV IsoPrime mass spectrometer at CRPG (Vandoeuvre-lès-Nancy, France). The results are expressed in δ notation in permil relative to the V-SMOW-SLAP scale. Two CRPG internal water standards were analyzed along with each weekly batch of water samples obtained from the fluid inclusions. To detect any potential effect of memory on the D/H ratios within the uranium reactor, the samples were duplicated whenever possible. This effect was not detected, but slight variations in δD for some samples were recorded. Given that these variations probably reflect variations in the δD composition within the fluid inclusion population, all of our duplicate measurements are presented here. They are denoted by adding the letter a, b, or c to the sample name in Table 1. The reproducibility for the standards was estimated to be lower than 2‰.

4.3. Host Minerals
4.3.1. δ18O and δ13C in Dolostones and Hydrothermal Calcite

Isotopic measurements were taken at the CRPG laboratory (Vandoeuvre-lès-Nancy, France). The carbon and oxygen isotopic compositions of the dolostones and hydrothermal calcite were measured using a GasBench autosampler coupled to a Thermo Scientific MAT253 isotope ratio mass spectrometer (IRMS). The residual dolostone and hydrothermal calcite powders were collected after crushing in order to obtain the δD measurements for the fluid inclusions. An aliquot between 200 and 400 μg of fine powder was reacted with 2 mL of supersaturated orthophosphoric acid at 70°C for at least 10 hours in a He atmosphere. Ten measurement cycles of the produced CO2 isotopic composition were performed with the Thermo Scientific MAT 253 continuous flow isotope ratio mass spectrometer. The isotopic compositions are given in the delta notation in ‰ relative to V-PDB for carbon and converted to V-SMOW for oxygen. The sample measurements were adjusted to the internal reference material (Merck synthetic calcite) calibrated using international standards IAEA CO-1, IAEA CO-8, and NBS 19.

In order to remain consistent with the D/H measurements, the samples were duplicated and the values are presented (denoted as “a,” “b,” or “c”). The reproducibility for the standard was less than 0.05‰ for δ13C and 0.1‰ for δ18O.

4.3.2. δ18O of Quartz

The δ18O isotopic analysis of the quartz samples was performed at the Laboratoire de Géochimie des Isotopes Stables (Institut de Physique du Globe de Paris and Université Paris 7). The isotopic compositions are given in the standard delta notation in ‰ relative to V-SMOW for oxygen. The samples were finely crushed and reacted with BrF5 as an oxidizing agent overnight in Ni tubes at 550°C following the method described by Clayton and Mayeda [46]. The O2 was directly analyzed on a Thermo Fisher DELTA V mass spectrometer. During the course of the analyses, the NBS28 quartz standard measurements were used to confirm the general reliability of the protocol (mean , ). In order to stay consistent with the δD measurements, the samples were duplicated and both values are presented (denoted as “a” or “b”). Each sample corresponds to the residual quartz powder after crushing.

5. Results

5.1. Description and Salinity of the Fluid Inclusions

Fluid inclusions (FI) from the euhedral quartz and hydrothermal calcite display similar features (Figure 3): they are aqueous, most of them are saturated at an ambient temperature with respect to halite (i.e., they contain a halite cube), and their liquid/vapor ratio ranges between 10% and 30%. The salinity measurements display a narrow and similar range for both the euhedral quartz and hydrothermal calcite with an average salinity of 28.6 wt% eq. NaCl and 30.5 wt% eq. NaCl, respectively (Figure 5 and Table 2). The homogenization temperature (Th) of the fluid inclusion in the quartz and calcite ranges from 130° to 230°C. Traces of gases, in particular N2 (>60 mol% in the volatile phase and minor amounts of CO2 and CH4), were found by Raman spectroscopy.

There are numerous but small fluid inclusions in dolostones () making it difficult to carry out petrographic and microthermometric studies. The studied fluid inclusions are aqueous (L-V inclusions), with a high salinity ranging between 20.8 and 29.6 wt% eq. NaCl with a mean value at 25.0 wt% eq. NaCl (Figure 5 and Table 2). Due to the small size of the fluid inclusions in the dolostones, it has been impossible to determine the traces of gas.

For the regional samples (Caussou 1), the aqueous fluid inclusions in the quartz do not have halite cube (Figure 4(b)). The mean salinity value is at 25.6 wt% eq. NaCl (Table 2). The Th values range between 190 and 200°C. Gases are only present as traces of predominant CO2 (>95 mol% in the volatile phase and several mol% of N2).

For the different minerals, it appears that the features of the fluid inclusions are similar with high salinity, a homogenization temperature close to 200°C, and, when it has been measured, traces of gas.

5.2. Cl/Br, δ37Cl, and δD of the Fluid Inclusions and δ18O and δ13C of the Host Minerals

All of the results are presented in Table 1. The fluid inclusions hosted in the euhedral quartz, hydrothermal calcite, and dolostone display very similar isotopic and Cl/Br ratios. Figure 6 shows the Cl/Br ratios as a function of the chlorinity. Based only on the Cl/Br ratio, all of the samples fall in the field of evaporated seawater saturated with respect to halite. The fluid inclusions hosted in the dolostones and hydrothermal calcites have similar Cl/Br ratios centered at around 350 whereas the fluid inclusions hosted in the quartz have a higher Cl/Br ratio of up to 609. The fluid from the Caussou calcite sample is the only one to fall above the Cl/Br value of seawater reaching 682 (Table 1). Based on Figure 6, no correlation is identified between the Cl/Br and chlorinity variations in the fluid inclusions hosted in the different samples.

For the δ37Cl values, most of the samples fall in the field of evaporated seawater between 0‰ and -1‰ (Figure 7). The fluids from one quartz and two hydrothermal calcites shifted to positive δ37Cl values with only one sample significantly shifting up to 0.87‰ (Table 1). This sample is a hydrothermal calcite sampled in the La Porteille area, the northern extension of the deposit (Figure 1). As represented in Figure 7, when combined with the Cl/Br ratio, numerous samples plot out of the evaporated seawater field. This could suggest the mixing of evaporated seawater saturated with respect to halite with seawater and/or with brine formed by the dissolution of halite (secondary brine).

For the δD and δ18O values, most of the samples (except two from Caussou) have similar δ18O values ranging between 9‰ and 12.2‰ whereas the δD values are more variable and range between -25.6‰ and -66.6‰. Variations in the δD and δ18O values are correlated to the mineral in which the fluid inclusions are hosted. The euhedral quartz has higher δ18O values, the dolostones have intermediate values, and the hydrothermal calcites have the lowest δ18O values. For the δD values, the euhedral quartz and dolostone have similar values whereas the hydrothermal calcites display the lowest ones (around -65‰). These relationships are well illustrated in Figure 8 which displays the δD values directly measured in the fluid inclusions and the δ18O values of the fluid recalculated using the oxygen isotopic composition of the hosted minerals and the adapted isotopic fractionation coefficients [47, 48] at 300°C. As mentioned previously, only the quartz and calcite samples from Caussou have significantly shifted from the main field with δ18O values of 23.5‰ and 20.7‰, respectively, and δD values of -45‰ for the quartz and between -36.5‰ and -38.3‰ for the calcite.

For δ13C values, both hydrothermal calcite and dolostone have similar values (between -3.2‰ and -2.2‰ and between -0.8‰ and -0.3‰, respectively) than those presented and discussed by Boulvais et al. [6]. Calcite from Caussou has a carbon isotopic composition of -0.8‰.

6. Discussion

The first part of this section deals with the characterization of the nature of the mineralizing fluid(s) at the origin of the magnesian metasomatism expressed at Trimouns through widespread dolomitization and talc-chlorite formation. Although previous studies have argued that the main fluid(s) at the origin of the talc-chlorite ore deposit was/were probably brine(s) that circulated at a temperature of [4, 5, 7, 8, 31, 34], the nature of the brine(s) (primary vs. secondary) and the possible involvement of other types of fluid(s) have not been clearly addressed. The interaction of the mineralizing brine(s) with the hosting rock of the deposits is also discussed in order to determine if the dolomitization and talc-chlorite formation result from a single fluid circulation event as previously suggested by Boulvais et al. [6]. Finally, the signification of the fluid event documented in Caussou (Figure 1) is discussed in order to determine the potential spatial extent of the mineralizing brine circulation system at a multikilometer scale.

The second part of this discussion deals with the determination of the source of the mineralizing brine(s) based on the particular geological context of the area during the Albo-Cenomanian period. The goal of this work is to propose a fluid circulation model at the crustal scale that is focused on the Trimouns deposit.

6.1. Nature, Interaction, and Spatial Extent of the Mineralizing Brine(s) in the Trimouns Area
6.1.1. Primary Signature of the Mineralizing Brine(s)

(1) Chlorine Proxies (Salinity, Cl/Br, and δ37Cl). The fluid inclusions hosted by the minerals that are synchronous with the talc primarily have high salinities between 27.5-32.0 wt% eq. NaCl and 29.2-33.5 wt% eq. NaCl for the euhedral quartz and hydrothermal calcite, respectively (Figure 5 and Table 2). However, few moderate salinities are also measured in the fluid inclusions (quartz: 10-15 wt% eq. NaCl; calcite: 12-18 wt% eq. NaCl). There are two ranges for salinity in the dolostones: a major one corresponding to highly saline fluids at approximately 25 wt% eq. NaCl (20.8 and 29.6 wt% eq. NaCl, max at 29.6 wt% eq. NaCl) and a second minor one corresponding to a more diluted fluid at about 12-15 wt% eq. NaCl. The general high salinity of the fluids is consistent with the fact that the main mineralizing fluid has passed the halite saturation point. However, moderate salinities should reflect minor mixing with another more diluted fluid.

The δ37Cl signature and the Cl/Br ratio of a fluid are powerful proxies to determine the nature of the brines [4951]. Fontes and Matray [51] have demonstrated that during the evaporation of seawater, the Cl/Br ratio has a typical evolution correlated with the increasing chlorinity of the residual solution. Here, our data, gathered as boxes, plot close to the evaporation line (Figure 6). The box for the dolostones is consistent with the area of the primary brines saturated with respect to halite. For the hydrothermal calcite and euhedral quartz, the Cl/Br ratios are also consistent with those expected for a primary brine having passed the halite saturation point; however, the chlorinities are slightly to significantly higher (up to ). Two main processes may increase the chlorinity. The first one is a reequilibration of the primary brines with halite at higher temperatures—the solubility of halite increases with temperature. In this case, the Cl content is increased but not the Br content and the data are shifted to higher Cl/Br ratios. The second process corresponds to an uptake of water. The replacement of anhydrous minerals (dolomite, feldspars, and plagioclases) by hydrated minerals (talc and chlorite) increases the chlorinity during the alteration of the crystalline basement by fluids. This process has already been suggested in unconformity U-type deposits hosted in the Athabasca Basin in Canada [52, 53], in the Modum Complex in Norway [39], and in the Schwarzwald district in Germany [54, 55].

If we consider the δ37Cl values, all of the fluid inclusions hosted in the dolostones, euhedral quartz, and hydrothermal calcite have isotopic compositions (-0.58‰ to 0.14‰, Table 1) that are consistent with those expected for evaporated seawater (, [49, 50]). Only one sample (PORT-2B, ) displays a significantly higher value, which may reflect the influence of a more intense fluid/basement interaction since biotite and amphibolite are known to be the main Cl-bearing minerals in igneous and metamorphic rocks and to have positive δ37Cl values ([53] and references herein). With respect to the combination of the Cl/Br ratio and the chlorine isotopic compositions (Figure 7), the whole magnesian metasomatic sequence is either in or close to the field of evaporated seawater.

(2) Hydrogen and Oxygen Isotopic Systems. Figure 8 represents the hydrogen and oxygen isotopic composition of the water trapped in fluid inclusions hosted in the studied samples calculated at a temperature of 300°C. This figure also shows the position of (i) the seawater, (ii) the present-day Global Meteoric Water Line, (iii) the seawater evaporation trend [5658], (iv) the isotopic composition of the talc and chlorite from Trimouns [6], and (v) the water in isotopic equilibrium with the talc (oxygen isotope data from [59]; hydrogen isotope data from [60]).

All of the waters from the dolostones, euhedral quartz, and hydrothermal calcite are broadly either in or close to the ideal range of isotopic compositions of the evaporated seawater saturated with respect to halite. The range of isotopic compositions of the water in equilibrium with the talc from Trimouns [6] also plots in this field. In more detail, all of our data are slightly depleted in deuterium and euhedral quartz and are also slightly shifted to more positive δ18O values. This general trend toward more negative δD values can be explained by mixing evaporated seawater with a low δD fluid. The possible source of this type of fluid is discussed in the following section. For δ18O and more specifically for the euhedral quartz, the shift to more positive values can be explained by a buffering of the mineralizing fluid by the talc and chlorite due to a low fluid/rock ratio. The oxygen (and hydrogen) isotopic composition of the water in isotopic equilibrium with the euhedral quartz falls in, or is really close to, the chlorite and talc field (Figure 8) from the deposit. Thus, the H2O uptake during the talc and chlorite formation is possibly at the origin of both the increase in chlorinity recorded in the quartz and calcite and the isotopic buffering of the fluid by the ore. Surprisingly, calcite, which is considered to be one of the last minerals to have precipitated, does not seem to be particularly affected by this process, suggesting that it formed under a higher fluid/rock ratio than quartz.

6.1.2. Evidence of Mixing between Primary Brines and Dilute Formation Waters

Chlorine proxies and stable isotope proxies provide consistent signatures indicating that the primary brines saturated with respect to halite are the brines that form the mineralizing fluid. However, the chlorinity ranges that include the lowest values close to 10% eq. NaCl, certain Cl/Br ratio values close to 650, and the general trend toward low δD values are all indications of partial mixing with other fluids.

Therefore, a mixing between more dilute waters may explain the slight scattering of certain values. Among the dilute waters, several endmembers may be considered: (i) recharge waters fed by meteoric inputs: such recharge waters may have leached the Triassic evaporites in some instances, which in turn could have resulted in an increase in Cl and are then the so-called secondary brines, (ii) (Albian) seawater, and (iii) formation waters which are known to have low δD values (up to -70‰) and oxygen isotopic compositions ranging between -10‰ and +10‰ [61]. These waters underwent various processes (e.g., a fluid-rock interaction implying chemical and isotopic modification, salt dissolution, and ion ultrafiltration) which altered the initial seawater composition [62]. In the geological context of this study, such waters are expected to have been hosted in the porosity of the overlying sediments. Based on the general trend of low δD values, the hypothesis of the mixing of primary brines with either secondary brines that have a meteoric origin or seawater is unlikely. Conversely, a mixing of evaporated seawater with formation waters should easily explain the low δD values and the slight salinity and Cl/Br data discrepancies.

6.1.3. From Dolomitization to Talc Formation

The dolomitization of the Devonian calcitic marbles is considered by Boulvais et al. [6] as belonging to a same and continuous water/rock interaction process that finally led to the formation of the talc. Our data support this interpretation and suggest that the whole magnesian metasomatic sequence (dolomitization-talcification-chloritization) results from the interaction between primary brines and rocks (calcitic marbles, pegmatites, and mica schists). Accordingly, the dolomitization of the calcitic marbles precedes the formation of the talc and chlorite. A significant supply of magnesium is therefore needed. As shown by Fontes and Matray [51], primary brines saturated with respect to halite are significantly enriched in magnesium, up to 50500 mg·L-1. We therefore propose that at least part of the magnesium needed for Mg-rich mineral formation may have been supplied by the evaporitic brine itself. This hypothesis could be consistent with the fact that, to date, the primary brines trapped in the deposit are depleted in magnesium [5].

As shown previously, a few fluid inclusions with lower salinities are also identified in the dolostones and show a large range of salinities from 18 to 28 wt% eq. NaCl and a few moderate saline fluid inclusions (Figure 5) reflecting mixing with formation waters. The changes in the degree of mixing between the dolomitization stage and the talc ore formation stage, in particular the increase in salinity and subsequently the cation concentration of the fluid (including the Mg2+ content), may possibly explain the changes in the main precipitated assemblages. It is also interesting to consider the role of CO2 during the magnesian metasomatism. Boiron et al. [4, 5] identified brines with traces of N2- in the euhedral quartz from the geodes. The presence of a low-density volatile phase dominated by N2, e.g., reflecting a rather low PCO2 value, could be consistent with the results of experimental studies dealing with talc formation from dolostone. Indeed, Wan et al. [63] experimentally determined that the talc formation rate increases with increasing temperature and decreases with increasing PCO2.

6.1.4. What about the Regional Circulation of Primary Brines?

One regional sample, coming from a quartz-calcite vein system localized in the Devonian schists from Caussou, has recorded the circulation of similar brines as those studied in the Trimouns deposit, suggesting a regional extent of this fluid circulation. The average fluid salinities measured in in quartz and calcite from the vein are similar to those from Trimouns (25.6 wt% eq. NaCl) (Table 2). The plot of the Cl/Br ratio as a function of the Cl content (Figure 6) confirms that they are also primary brines saturated with respect to halite. In the Cl/Br vs. δ37Cl plot (Table 1, Figure 7), the brines plot really close to the field of evaporated seawater. The δD values are also close to those corresponding to the Trimouns brines. However, the calculated δ18O value pairs are significantly shifted with respect to the Trimouns data (Table 1, Figure 8). This difference could be due to (i) the involvement of another fluid, (ii) the use of an inappropriate fractionation coefficient between minerals and fluids, or (iii) the use of an inappropriate temperature (300°C) to calculate the δ18Ofluid composition. Based on the similarity of the Cl/Br ratios, the salinities, and the δ37Cl and δD values of the mineralizing fluid documented in Trimouns with the primary brine signature, it is unlikely that another source of fluid is involved. Different isotopic fractionation coefficients between quartz, calcite, and H2O are available in the literature and are relatively consistent between them. The isotopic fractionation coefficient is between 8.2‰ and 6.9‰ for quartz and H2O at 300°C [47, 6468] and between 3.8‰ and 5.8‰ for calcite and H2O at 300°C [48, 6972]; an exception for this is provided by Kim and O’Neil [73] who predict a fractionation coefficient of -1 between calcite and H2O at the same temperature. However, their slight differences cannot explain the large variation in the δ18Ofluid values compared to the primary brine composition field (~10‰). If two minerals occur together, it is possible to postulate that they were in isotopic equilibrium when they formed from the mineralizing fluid and to calculate their temperature of formation using the value. In this way, temperatures of ~120°C are calculated using the fractionation coefficient between quartz and calcite given by Zheng [48], ~150°C using that provided by Clayton and Keiffer [74], and 280°C using the one cited by Sharp and Kirschner [66]. The calculated oxygen isotopic composition of the fluid (using [47, 48]) falls in (at 120°C) or close to (at 150°C) the range of the primary brines saturated with respect to halite (: and calcite: and : and calcite: ). Even if this range of temperatures is rather somewhat low, the fluid documented in Caussou has very similar features to those of the mineralizing fluid and may reflect the cooling of the system at the end of the process.

6.2. Origin(s) of the Primary Brines
6.2.1. Primary Brines Expelled from the Triassic

The evaporitic Triassic was both widespread and distributed all along the Hercynian basement. Therefore, the involvement of Triassic brines in various geological contexts postdating the Trias evaporite formation is often inferred. Two main processes are often proposed. The first one consists of the expelling of the interstitial fluids (the salt porosity is up to several volume percentages, [7578]) hosted in the evaporites, which are the residual fluids after seawater evaporation and salt precipitation. This type of brine is called a primary brine. The second one consists of the dissolution of salt by a less saline fluid generating the so-called secondary brines. Based on our results, we have demonstrated that the fluids at the origin of the deposit were primary brines saturated with respect to halite mixed with a minor amount of formation waters. We propose that such primary brines could have been expelled from the Triassic (Figure 9) during the Cretaceous by (i) Albian sediment overloading on the underlying Mesozoic rocks and/or (ii) the shearing of the Triassic levels induced by the gravity gliding of the Mesozoic series [79].

6.2.2. What about the Influence of Mantle Metasomatism?

The recent literature underlines the fact that salt enrichment in fluids could originate from other processes, such as H2O removal when forming hydrated minerals from anhydrous ones ([39, 52, 54, 55, 58, 80]; Hovland et al., 2018). Carrying this idea a step further, Scribano et al. [81] suggested that the formation of a certain subseafloor giant salt deposit in deep marine basins may be genetically related to the serpentinization of ultramafic rocks. In the present context of passive margin hyperextension, mantle exhumation, and subsequent fluid circulation, several authors have already proposed that seawater-driven serpentinization may contribute to the transport of dissolved elements (Mg, Cr, and V; [82]). Recently, Corre et al. [83] suggested that such fluids circulated in the Saraillé Massif based on the Cr content of the chlorite localized in the detachment fault (average Cr content of 1.44 wt%). It is also interesting to note that in the latter case, these fluids are extremely saline with up to 27 wt% eq. NaCl.

Given the current proximity of the serpentinized peridotite bodies from the deposit (Figure 1), we suggest estimating the relevance of the oxygen and deuterium isotopic systems and chlorine proxies in order to test the hypothesis of mantle rock serpentinization and the subsequent generation of brines as a potential origin of the Trimouns deposit (Figure 9).

From an isotopic point of view, the generation of brine by seawater evaporation or serpentinization refers to two distinct processes leading to different fractionation effects. The first step consists of comparing the fluid isotopic composition documented here with the one expected for a fluid in isotopic equilibrium with serpentine. Unfortunately, the oxygen and deuterium isotopic compositions of the Pyrenean serpentines are not available. Despite this lack of data in the literature, we propose to consider the large range of modern oceanic serpentinites compiled by Pope et al. [84]. Considering a temperature of 300°C, it is possible to calculate the expected δD and δ18O values of water in isotopic equilibrium with these serpentines [60]. As displayed in Figure 8, the large field of composition of such water overlaps the field of evaporated seawater and therefore the Trimouns data. This may indicate that stable isotopes cannot discriminate between the two processes.

Bonifacie et al. [85] and Barnes and Cisneros [86] measured the chlorine isotopic compositions of the oceanic serpentinites, which have δ37Cl values ranging between -1.4‰ and +1.8‰. As explained by Bonifacie et al. [85], most serpentinites overlain by sediments show negative δ37Cl values interpreted as resulting from the interaction with sediment pore waters [87] that have negative δ37Cl values [45, 8890]. This implies that serpentinization does not disturb the chlorine isotopic composition of the fluid, which keeps its negative signature. Consequently, in our case, the δ37Cl value does not really help us to determine the origin of the brines.

As for the Cl/Br ratios, the chlorine behavior during the serpentinization of the peridotites is the only well constrained in the literature. The fresh oceanic lithosphere has low chlorine content (e.g., [91, 92]), whereas secondary minerals formed after a water-rock interaction can have high chlorine content, e.g., as serpentine ([93] and references herein). Consequently, serpentine can be considered as a “sink” for chlorine even if only a small amount of Cl can be incorporated at . With regard to the serpentinization-derived brine hypothesis, the interaction between the formation water (derived from Albian seawater) and the exhumed mantle should have led to a decreased chlorine content in the residual fluid. Thus, the expected Cl/Br ratio of the residual fluid should be lower than the initial one. Considering that the formation water probably was initially Albian seawater, the expected Cl/Br ratio of the formation water after the serpentinization (the residual fluid) should be lower than 650, a value that is consistent with our data. However, given that the behavior of the bromine during the serpentinization is poorly known, the expected Cl/Br composition of the residual fluid remains highly hypothetical.

7. New Insights into the Trimouns Magnesian Metasomatose Model

In this study, using complementary geochemical proxies (δD, δ18O, δ37Cl, and Cl/Br), we demonstrate that the formation of the talc-chlorite Trimouns deposit results from the interaction of primary brines saturated with respect to halite with the rocks lying on both sides of the fault located on the eastern border of the Saint Barthelemy Massif. Based on the chemical similarities of the fluid contained in the syn-deposit minerals (quartz and calcite) and antemineralization rocks (dolostones), we suggest that the whole magnesian metasomatic sequence (; pegmatite or mica ) results from successive pulses of primary brine saturated with respect to halite during a relatively long period of time (at least between 122 Ma and 96 Ma). We suggest that the evaporitic Triassic succession, which is both widespread and distributed on top of the Hercynian basement in the Pyrenees, is the main source of this fluid. In some cases, the chlorinity ranges that include the lowest values close to 10 wt% eq. NaCl, certain Cl/Br ratio values close to 650, and the general trend toward low δD values are all indications of partial mixing with lower salinity fluids such as formation waters hosted within the overlying sediments.

Given the particular geological context (hyperextension, mantle exhumation) in which the formation of the Trimouns deposit took place, we have also examined the possible influence of mantle metasomatism on the formation of primary brines. We conclude that, although our results are not entirely inconsistent with this hypothesis, the link with serpentinization is not demonstrated. Moreover, if brines derived from the serpentinization have been involved in the formation of the deposit, their circulation at the deposit level should be temporally disconnected from those of the brines derived from seawater evaporation. As suggested in the literature [54], it is difficult to invoke simultaneous downward and upward flows of the various fluids in order to explain the evidences for mixing in the deposit. From a geochemical point of view, additional studies are needed to clearly address the question of whether serpentinization can influence the formation of brines in the Pyrenees. First, a deuterium and oxygen isotopic survey of the Pyrenean serpentinites is required, and second, an equivalent survey focusing on the chlorine isotopic composition of the Pyrenean serpentinites would be also useful to test this hypothesis. Lastly, the characterization of the magnesium isotopic composition of the dolostones, talc, chlorite, and serpentinites would provide some interesting insights. As has already been suggested by several authors, serpentinization is not isochemical [9496] and magnesium should be transferred into the reacting fluid from peridotites during their alteration. This idea remains a matter of debate, and the comparing of the magnesium isotopic composition of the Pyrenean talc, chlorite, dolostones, and serpentinites could help to determine if the mantle provided part of the magnesium needed to form the deposit. However, to date, only evaporated seawater saturated with respect to halite is identified as the most likely source of magnesium.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was financially and technically supported by Total S.A. We thank Thibaut Klein (Imerys Talc) for giving us access to the Luzenac mining site. Serge Fourcade, Philippe de Parseval, and Alexandre Boutin are warmly thanked for stimulating discussions on the field. Antonin Richard and Pierre Martz are warmly thanked for the fruitful discussions about brines in general. We also thank Yves Lagabrielle and Philippe Boulvais for the general discussions about the geological history of the Pyrenees. Sylvain Calassou is thanked for his scientific advices about the geometry of the Pyrenean hyperextended passive margin. Marina Kouadio and Frederic Diot are thanked for their technical support on the fluid inclusion study and ion chromatography, respectively. Enrique Gomez-Rivas is thanked for his constructive review that helped the authors improve the clarity of the manuscript.