Fluid Origin and Evolution of the Roc Blanc Silver Deposit (Jebilet Massif, Variscan Belt, Morocco): Constraints from Geology and Fluid Inclusions

Pb-Zn-Ag-Au vein deposit is located in the NW of Marrakech, in the Central Jebilet massif. It is spatially related to Bramram-Tabouchennt-Bamega (BTB) granodioritic pluton (ca. 330 Ma) metamorphism aureole. The main veins hosted in black shales are oriented N-S to NNW-SSE. Pb-Zn-Ag-Au ore is associated with quartz, chlorite, sericite, and carbonate gangue minerals. Two major stages of ore deposition were distinguished. The preore stage (stage I) comprises two quartz-mineralised vein generations with Fe, As, Zn, and Cu ores (vg 1 and vg 2 ). The main ore stage (stage II) consists mainly on Ag, Au, Pb, Zn, Cu, and Sb ores, which is hosted by carbonaceous vein (vg 3 ) and by two late quartz generations veins (vg 4 and vg 5 with a geodic quartz). Three types of ﬂ uid inclusions have been recognized in silver mineralisation bearing quartz veins according to petrographic investigations, microthermometry, and Raman spectroscopy studies: (i) liquid-rich H 2 O-N 2 -CH 4 ± CO 2 -(salt) ﬂ uid inclusions (type 1), (ii) vapour-rich H 2 O-CO 2 -CH 4 -N 2 -(salt) ﬂ uid inclusions (type 2), and (iii) aqueous H 2 O-(salt) ﬂ uid inclusions (type 3). The interpretation of ﬂ uid inclusion data shows a mixing of two ﬂ uids that are metamorphic and surface to subsurface origin, trapped at boiling state. The ﬁ rst mineralised stage was deposited at 350 ± 20 ° C (this temperature of ore deposition was supported also by chlorite geothermometry) with salinity of 13.7wt% NaCl equiv., while the deposition of the argentiferous stage, which consists of the main economic mineralisation of the Roc Blanc deposit, occurs during decreasing temperature at 150 ° C with a salinity of 12.1 wt% NaCl equiv. The all-mineralised ore was deposited at relatively low pressure, below ~ 1-1.1kbar. So, ﬂ uid dilution and cooling are probably the main factor for silver deposition in the Roc Blanc polymetallic vein deposit. In addition, ﬂ uid inclusion studies reveal that the mineralising ﬂ uid corresponds to a mixture of metamorphic ﬂ uid (H 2 O-CH 4 -N 2 -CO 2 ) with surface to subsurface aqueous gas-free ﬂ uids (H 2 O-salt, meteoric, or brine).


Introduction
The Roc Blanc Pb-Zn-Ag-Au deposit is located on the central unit of the Jebilet Hercynian massif, which is the most important part of the Jebilet Pb-Zn-Ag-Au metallogenic province. It is hosted by metamorphic Visean-Namurian shales (spotted schists) in the contact metamorphic aureole of the Bramram-Tabouchennt-Bamega (BTB) granodioritic plutons ( Figure 1). The Roc Blanc silver deposit is spatially related to the BTB granodioritic plutons and probably displays a typical concentric base and precious metal zonation of Sn, W, Au, Cu, Pb, Zn, Ag, and Au, covering an area of approximately 10 km [1,2]. The Pb-Zn-Ag-Au mineralisation of the Roc Blanc is similar to some other deposits in the central Jebilet (i.e., Koudiat el Hamra polymetallic mineralisation) which have been mined for silver in the past [2].
The Roc Blanc deposit is located in the central part of the Jebilet massif, 20 km north of Marrakech city. Discovered in 1925, the deposit is a silver-rich ore, with indicated and inferred resources which have been estimated around 200,000 t at 297 g/t and >56,000 t at 234 g/t Ag, respectively. Resources contain about 0.8% Pb and 1.01% Zn and a significant amount of copper and gold [3].
The formation of the Roc Blanc deposit has been subject of many studies that are complementary, but often controversial [4,7], questioning the nature of the mineralising fluids. The ore-forming fluids play a major role in any hydrothermal mineralised system. Therefore, the characteristics, origin, and evolution of these fluids are important for metallogenic studies. Several mechanisms, like boiling, fluid mixing, cooling, decrease of fluid salinity and pressure, and fluid rock-interactions, have been considered for mineral precipitation factors [8,9]. More of them are taken account for the initiation of polymetallic silver precipitation, such as boiling, mixing, cooling, salinity, and pressure decreases (e.g., [4,[10][11][12][13][14][15][16][17][18]). In the present paper, we present paragenetic sequences of the Pb-Zn-Ag-Au veins of Roc Blanc deposit, with new analysis of fluid inclusions (petrography and microthermometry coupled with Raman spectroscopy analysis). The purpose of this study is to discuss the fluid composition and sources since fluid inclusion studies, timing, and metallogenic nature of the deposit through estimating the PT conditions of ore vein formation.  [19]. The Meseta domain is structurally limited by the Haut Atlas Mountain in the south and the Rif domain in the north. It is separated into two sectors (the Western and the Eastern Meseta) by the Mesozoic mountains of the Middle Atlas [20]. The Jebilet massif is located in the western part of the Moroccan Meseta, at 7 km north of the Marrakech city (Figure 1(a)). The massif is outcropping in E-W direction with 170 km in long and 40 km in wide.

Geological Setting
The Jebilet massif consists of three tectonostratigraphic units, subdivided into eastern, central, and western parts (Figure 1(b)), separated by major shear zones [21,22]. The western Meseta Shear Zone (WMSZ) separates the western units from the central one [23], and the Marrakech Shear Zone (MSZ) separates the central from the eastern unit [24]. Their original sedimentary relationships are masked by these complex tectonized contacts. The eastern unit is composed of Visean rocks including allochtonous terranes of Ordovician to Devonian ages. The western units consist of weakly deformed and metamorphosed Cambrian and Ordovician rocks, whereas the central units consist mainly of very low-grade metamorphism of the Visean Sarhlef Formation. Recently, new biostratigraphic results on metamorphosed conodonts studied by Lazreq et al. [25] have been attributed to the Sarhlef Formation. These new biostratigraphic investigations have yielded an age ranging from the Frasnian to late Famennian for the Sarhlef serie [25]. According to these authors, this new age is more compatible with the previous radiometric data of the intrusive magmatic rocks [26][27][28][29][30]. The central units are intruded by three generations of syntectonic intrusions. The earlier igneous activity has been started at c.a. 330 Ma [31] with the emplacement of an early bimodal mafic-ultramafic to felsic suite, comprising olivine to quartz-bearing gabbro and peridotite locally intruded by a ca. 279 ± 7 Ma gabbro dike. The whole was intruded by late Variscan granitoid stocks, leucogranite dike, and microdiorite sheets. Granitoid intrusions correspond to a syntectonic I-type and/or S-type, strongly peraluminous, largely evidenced by the Bamega-Tabouchennt-Bramram stocks and Oulad Ouaslam laccolith [32][33][34]. These plutonic rocks have been dated to 327 ± 4 Ma (whole-rock Rb-Sr ages; [30]). Two highly differentiated micas ± tourmaline-bearing leucogranites have been recognized in the studied area, and they correspond to a NW-trending, subvertical sheet dikes, dated at 296 ± 12 Ma (Rb-Sr isotopes; [30,35]). However, Ar-Ar dating method on white micas from regional shear zones that also affects the leucogranite, gives an oldest age 305:9 ± 0:9 Ma [7]. Microdioritic dikes correspond to the youngest intrusive rocks, dated from ca. 255 -241 ± 5 Ma (K-Ar and 40 K- 40 Ar on kaersutite [36][37][38]) to ca. 235 Ma (SHRIMP U-Th-Pb on zircon [27]). Before these intrusive magmatic rocks, a volcanic rock has been described in Central Jebilet, associated with pyrrhotite-rich massive sulphide deposits, attributed to Upper Visean (e.g., [22,29,[39][40][41]).

Geofluids
Aureoles of hornblende and pyroxene hornfels metamorphic facies are particularly developed in the vicinity of felsic intrusions [33]. For deformation events in the central Jebilet unit, they are distributed unevenly, depending mainly on the relative competence of the rocks [43]. So, many deformation phases have been noticed in this central part of the massif. Four deformation phases named as D 0 , D 1 , D 2 , and D 3 have been recognized ( [22,[46][47][48] and references therein). D 0 was identified as an early synsedimentary event related to extensive-transtension tectonics, allowing the setting up of the Devono-Carboniferous basin [46]. The D 1 and D 2 deformation phases were related to the closure of the basin. Those episodes begin in Upper Visean-Namurian (~325 Ma). The D 3 event was related to the emplacement of late Variscan, syn-to postkinematic granodioritic intrusions, which developed an HT-LP metamorphic aureole.
2.2. The Roc Blanc Deposit Geological Setting. The Roc Blanc deposit is hosted by black shales of the Sarhlef schist serie. This serie covers most of the central unit (Figure 1(b)) and is composed of sedimentary rocks weekly metamorphosed to greenschist facies. The lithology is made up mainly of pelites (black shales), sandstones, and calcareous beds alternating with volcanic and subvolcanic rocks of mafic and felsic composition [22]. The schists are intruded by Bramram-Tabouchennt-Bamega (BTB) granodioritic plutons, leucogranites, and microdioritic dikes [1] at the western boundary of the Marrakech Shear Zone (MSZ), which limits the central part and the eastern part of the Jebilet massif. The deposit is located in the contact metamorphism aureole of the BTB plutons ([1], Figures 1(b) and 1(c)). These plutons contain numerous quartz-carbonate veins, pegmatitic veins, and aplitic rocks.
The Roc Blanc orebodies have been described as quartzcarbonate veins cutting the Sarhlef schists [12,22]. The mineralisation is characterised by [4,5] (i) a complex paragenesis; (ii) modest opening; (iii) direction substantially parallel to the major Hercynian fold axial plan, often slightly inclined on the schistosity; and (iv) its richness in silver and antimony minerals. Mineralogically, the Roc Blanc deposit is particularly distinguished by the presence of owyheeite and andorite, in addition to other Ag-rich bearing minerals (argentiferous tetrahedrite, ullmannite, miargyrite, pyrargyrite, polybasite, argentite, and native silver). Ag-rich minerals commonly are expressed as inclusions in galena, chalcopyrite, and other sulphide and occur sometimes as microfractures filling sulphides and quartz.

Sampling and Analytical Methods
Samples were taken from outcrops and from thirteen drill cores (RB-1, RB-2, RB-3, up to RB-13; Figure 2(a)) of the Roc Blanc deposit for ore descriptions. The mineralogical and petrological studies, already published by Arbaoui et al. [4], are based on results of a sample set, including eighty polished thin sections. Those samples were studied under transmitted-and reflected-light, polarized microscope, and scanning electron microscope (SEM), with an energy dispersive spectrometer (EDS, SE, and BSE) at the REMINEX Research Centre and Laboratory (Morocco). SEM operating conditions were an accelerating voltage of 30 kV, a beam current of 20 nA, and count times of 20 s.
Five samples were selected for fluid inclusion studies, carried out on double polished, 100-300 mm thick sections from quartz bearing mineralisation, including argentiferous phases. Various inclusions in quartz were selected for micro-thermometry and Raman spectroscopy analysis based on their petrographic characteristics. The petrographic studies of inclusions were mainly conducted according to Roedder [51], Van den Kerkhof and Hein [52], and Diamond [53].
Microthermometric studies were performed on thick sections using a Linkam THMS-600 heating-freezing stage mounted on an Olympus BX-50 microscope at the Laboratory of Metallogeny and Magmatism in the Geology Department of STHB University (Algeria). Liquid nitrogen was used for the freezing phase. The machine temperatures range from -196 to 600°C. For all inclusions, the accuracy was 0.2°C during the cooling stage and 0.5°C during the heating stage. The volumetric fraction of the aqueous phase was visually estimated using standard charts [51,54]. The salinities of aqueous gas-free inclusions were calculated using the final ice melting temperatures with the Bodnar equation [55,56]. Molar fractions of CO 2 , CH 4 , and N 2 were determined in individual fluid inclusions by semiquantitative analysis of Raman spectroscopy spectra and performed with a Labram Raman spectrometer at CNRS-CREGU (Nancy, France). The scanning regions were in the range of 1000-1500 cm -1 for CO 2 , 2250-2750 cm-1 for N 2 and H 2 S, and 2750-2950 cm -1 for CH 4 . Finally, isochoric curves of aqueous gas-free inclusions were calculated using the model of Zhang and Frantz [57] to estimate the PT condition of the main mineralising fluid.  [4]. Seven major ore structures (in extension) are identified in the deposit area, occurring as discontinued veins of several hundred meters long, hosted by black shales (Figure 2(a)). They show two major trends that are N135 to N150 (V 3 , V 4 , V 6 , and V 7 ) and N0 to N10 (V 1 , V 2 , and V 5 ) (Figure 2(a)). In most cases, hydrothermal veins in the study area are mainly oriented around N135 and dip from 45°to 65°E (Figure 2(b)). All of the major ore structures are located in the sinistral SSEtrending Marrakech Shear Zone. The mean extension of major ore structures varies from 300 to 600 m long and from 0.3 to 0.5 m of thickness (Figure 2(a)). V 5 is the largest one, being more than 2 000 m long with an average thickness of 0.7 m. The veins cut across the S 1 cleavage, which strikes N10 to N45 and dips 70°-80°E (Figure 2(b)). Stretching lineations (L e ) are also observed; they are oriented toward the NNE direction, along with S 1 cleavage. The stratification in the deposit area is generally transposed into S 1 cleavage plan; however, some ESE strikes are occasionally observable. At map scale, microdioritic dikes, oriented NNE, also cut across all S 1 cleavage, Marrakech Shear Zone, and postdate the major ore structures. These dikes are clearly intersected by sinistral and dextral faults with local breccia mineralisation at their intersection ( Figure 2(b)). The well relationship between these faults and mineralised veins is not observed, but generally, mineralised veins are stopped in front of these faults (e.g., V 2 and V 5 ). 4 Geofluids
All of these mineral assemblages in the Roc Blanc deposit were regrouped into three main ore forming stages ( [4,5]; Figure 4): (1) The first stage is composed of Fe, As, Zn, and Cu mineralisation associated with quartz (Qtz1, 2), chlo-rite, and sericite; it is represented by vg 1 and vg 2 mineral assemblages. Sphalerite and chlorite of this stage are Fe-rich. The mean value of Fe in sphalerite is 5.71 wt%. Chlorite corresponds to a ripidolites with a Fe # (=Fe/(Fe + Mg) value ranging from 0.66 to 0.84 (Table 2 and Figure 5(a)). Kranidiotis and Mac Lean [58] geothermometry of these chlorites reveal temperatures between 330 and 380°C ( Figure 5(b)). The quartz-chlorite-mineralised veins could be formed at 350 ± 10°C, similar to those associated with the Kettara pyrrhotite-rich massive sulphide deposit [50,59] (2) The main ore stage is characterised by Ag, Au, Pb, Zn, Cu, and Sb mineralisation. The precipitation of Ag-bearing minerals is preceded by base metal sulphide minerals. This stage is formed by vg 3 , vg 4 (Qtz3), and vg 5 6 Geofluids and >91 wt% Ag, respectively (Table 1 and Figure 5(c)). Ag-Au alloy is intergrown closely with argentite and typically cuts or replaces galena, argentiferous tetrahedrite, and chalcopyrite. It contains up to 80.21 wt% Au (Table 1) (3) The postore stage is characterised by barren fluid circulation causing some alterations with chlorite, sericite, and quartz minerals. In addition, the primary ore of the Roc Blanc deposit is overprinted by supergene alteration marked by the formation of secondary pyrite, marcasite, magnetite, hematite, goethite, and malachite/azurite

Fluid Inclusions
Fluid inclusion studies were conducted on quartz in the studied five mineralised veins, including three early-stage and two late-stage samples. On the one hand, only one sample was collected at the early stage of the mineralised vg1type which is characterised by quartz 1 + chlorite ± Fe-rich sphalerite ± galena. On the other hand, the other two samples were sampled from the vg2-type consisting of quartz 2 + siderite/ankerite + pyrite + arsenopyrite. The two samples of the late stage are from vg 4 with quartz (Qtz3) + chalcopyrite + ullmannite + (Ag-Sb) sulphosalts and geodic quartz vein (vg 5 ) with quartz (Qtz4) + calcite + Fe-poor sphalerite 7 Geofluids + argentiferous tetrahedrite + Ag-and Au-bearing minerals. Carbonate minerals associated with the mineralisation are poor in fluid inclusions, and the scarce one observed is too small in size (<5 μm) and was not suitable for microthermometric studies.

Fluid Inclusion Petrography.
Fluid inclusion studies were limited to primary fluid inclusions. They form isolated or fluid inclusion assemblage (FIA), [53] in host minerals and do not belong to any fractures. Based on distribution and shape, the vapour-filling ratio (R Flv ), microthermometric, and Raman spectroscopy analyses, three main fluid inclusion types are identified in the quartz-mineralised veins (Qtz1 (vg 1 ), Qtz2 (vg 2 ), Qtz3 (vg 4 ), and Qtz4 (vg 5 )) of the Roc Blanc deposit. Type 1 and type 2 are two-phase (L and V) aqueousgas-bearing fluid inclusions and are observed in Qtz1, Qtz2, and Qtz3. However, type 1 is liquid-rich (H 2 O-N 2 -CH 4 ±CO 2 -(salt)) while type 2 is relatively vapour-rich (H 2 O-CO 2 -CH 4 -N 2 -(salt)). Type 3 is essentially twophase (L and V) aqueous fluid inclusions (H 2 O-(salt)), existing in all quartz generations (Figures 6(a)-6(f)). These inclusions are associated with scarce one-phase fluid inclusions in some quartz. The three main types are distributed randomly and would correspond to primary fluid inclusions compared to the secondary fluid inclusions, which are arranged along transgranular trails [52]. These late inclusions consist of two-phase fluid inclusions, are small in size (<5 μm), and were not suitable for microthermometric studies (Figure 6(d)).

Microthermometry and Raman
Spectroscopy. Application of laser Raman spectroscopy indicates the presence of volatile components in the vapour phase of type 1 and type 2 fluid inclusions. The Raman spectra collected correspond mainly to nitrogen, methane, and carbon dioxide (Figures 7(a)-7(d)). Type 3 fluid inclusions are essentially gas-free, whereas one inclusion analysed from a scarce one-phase fluid inclusions contains methane and nitrogen (Qtz2, 57% CH 4 , and 43% N 2 ). Although the volatile phases (CH 4 , N 2 , and CO 2 ) are detected in many inclusions, twophase transitions have been noticed during freezing and heating runs: ice melting temperature T mðiceÞ and bulk homogenisation temperature (T h ) ( Table 3). Clathrates form generally in the freezing stage of these inclusion types; however, clathrate melting temperatures were difficult to observe.

Type 1 Fluid
Inclusions. These inclusions have been identified in Qtz1, 2, and 3. Their microthermometric data show a variable ice melting temperature. T mðiceÞ are generally between -15 and -1°C in Qtz1, Qtz2, and Qtz3, with an average value of -6.9, -7.6, and -10°C, respectively (Figure 8(a)). In the Qtz1 and Qtz3, scarce fluid inclusions yield a low T mðiceÞ up to -18 and -21.8°C, respectively. During the heating stage, all of these inclusions have been homogenised into the liquid phase. T h are relatively similar in the Qtz1 and Qtz3; they range from 104 to 237°C in the Qtz1 and from 117 to 248°C in Qtz3, with average values of 164 and 190°C, respectively (Figure 8(b)). In the Qtz2, T h are slightly higher, with values between 215 and 288°C, with a mean around 250°C. Raman spectroscopy analysis was carried out in these inclusion types in the Qtz1 and Qtz2. Their vapour phase is mainly composed of variable proportions of nitrogen and methane, whereas carbon dioxide is detected as traces in two inclusions of Qtz1 (Figure 9). According to the semiquantitative composition (X in mole percent) of gases calculated from Raman spectra areas, in the Qtz1, nitrogen is relatively dominant; X N2 ranges from 59 to 77 mol%, and X CH4 ranges from 23 to 39 mol%. Their average composition is approximately 70 mol% N 2 and 30 mol% CH 4 . In the two inclusions where carbon dioxide is detected, X CO2 is 2 mol%. In the Qtz2, methane is prevalent, X CH4 ranges from 52 to 74 mol% with a mean of 60 mol%, and X N2 is between 26 and 48 with an average value of 40 mol%. CO 2 is missing in all inclusions of type 1 in the Qtz2 or its content is lower than the detection limit.    Qtz2 Qtz1 (f) 14 Geofluids According to the abundance of nitrogen and methane, a relative absence of CO 2 , and the ice melting temperature variation, the type 1 fluid inclusions are considered to belong to H 2 O-N 2 -CH 4 ±CO 2 -(salt) system.

Type 2 Fluid
Inclusions. They show a relative lower ice melting and higher bulk homogenisation temperatures than the other fluid inclusion types of the Roc Blanc mineralising systems (Figures 8(c) and 8(d)). Contrary to the type 1, the type 2 fluid inclusions have almost homogenised into vapour phase. In the Qtz1, T mðiceÞ ranges from -15.8 to -8°C and T h from 250 to 310°C, with means of -10.8 and 280°C, respectively. In the Qtz2, these inclusions bear a moderate saline fluid according to their low ice melting temperature; T mðiceÞ are between -24 and -12.8°C with an average value of -18°C. Their T h are also higher than those of the type 1 with values ranging from 300 up to 377°C and a mean around 340°C. This type in Qtz3 shows T mðiceÞ between -23 and -5°C, mean T mðiceÞ = −12°C, and their T h values ranging from 250 to 309°C, with a mean of about 280°C (Figures 8(c) and 8(d)). The type 2 is also aqueous-gas-bearing fluid inclusions, with gas composed of CO 2 , N 2 , and CH 4 ; their average compositions are, respectively, 40, 30, and 30 mol% in the Qtz1; 36, 27, and 37 mol% in the Qtz2; and 78, 13, and 9 mol% in Qtz3 ( Figure 9). Compared to the type 1, it is not argued with the few Raman spectroscopy data that CO 2 is systematically present in all inclusions of the type 2; however, its content here is relatively dominant ( prograde metamorphic reactions. The black shales of the central Jebilet are enriched in organic carbon and phyllosilicates, so the production of CH 4 -N 2 fluids could be possible during prograde metamorphism, and they entrapment during ore forming process which postdate the pic of regional metamorphism. In summary and based on the arguments presented above, the aqueous-gas-bearing fluid inclusions at the Roc Blanc deposit are consistent with the contribution of metamorphic fluids in the mineralising system and were probably mixed with aqueous relatively low temperature fluids. The mineralised fluids were eventually trapped at boiling state. 6.2. Timing of Mineralising Events and P-T Conditions. The mineralised veins postdate S 1 cleavage, which is axial planar to folds generated during the post-Visean tectonometamorphic event ( [45], [4]). So, the mineralisation was formed during or after the retrograde cooling phase of the metamorphism; indeed, the composition of aqueous gas-bearing fluid inclusions yields partly a metamorphic origin as emphasised above. The presence of aqueous fluid inclusions (type 3) in all quartz-mineralised veins argued the retrogression metamorphic cooling, which largely involves rehydration reactions caused by waters from outside [63,68] and that have probably flooded the sink site.

Qtz 4 Vg5
Th mean 17 Geofluids pressure of about 1 and 1.1 kbar, respectively, relative to the deposition of the argentiferous stages ( Figure 11). This must represent the pressure around which the Roc Blanc deposit was formed, potentially near 500 bars. Indeed, the mineralised veins at Roc Blanc deposit cross-cut each other, showing successive emplacement in time, but their inclusions yielded similar fluid compositions, which show that veins were formed relatively in the same context short spaced times.
On the other hand, the highest homogenisation temperatures are generally recorded in the Qtz2 (Figures 8(b), 8(d), 8(f), and 10(d)); the T h of its fluid inclusions would represent an upper limit temperature of the Roc Blanc mineralising fluids (350 ± 20°C). These conditions work with the formation of the first quartz vein generations (vg 1 and vg 2 ). For the argentiferous phases, its deposition occurs at c. 100 to 210°C according to the T h value of aqueous gasfree fluid inclusions of Qtz3 (vg 4 ) and Qtz4 (vg 5 ).
For ore deposition factors, no major characteristic has been observed under fluid inclusion studies to explain the difference in mineralogy of the ore stages. Considering all quartz-mineralised veins, a few significant observations are done in the type 3 fluid inclusions (observed in all quartz) where the mean T h increases from Qtz1 (vg 1 ) to Qzt2 (vg 2 ) and decreases progressively towards Qtz3 (vg 4 ) and Qtz4 (vg 5 ) (Figure 10(d)). The same temperature evolution is shown also by type 1 and 2 fluid inclusions in vg 1 , vg 2 , and vg 4 quartz. These statements can indicate a rapid heating (vg 1 to vg 2 ) followed by progressive cooling (vg 4 to vg 5 ) in the general system. Note that silver-bearing minerals form in the latest veining system (Qtz3 and Qtz4); their genesis can then be linked to the cooling of mineralising fluids.
6.3. Metallogenic Type. The Roc Blanc Pb-Zn-Ag-Au veins are located at 2 km towards the NNE of the BTB Hercynian granodiorites. It is then spatially associated with the contact metamorphic aureole of this intrusion which reaches 5 km in diameter. The granodiorites and their surrounding aureoles (black shales) host several mineralised veins that are Sn-W mineralisation and Cu, Au, Sb, Pb, Ag, and Zn mineralisation, respectively, suggesting a metallic zoning around the BTB granitoids [2,22]. This zoning could be represented by likely: Sn-W (1) → Au-Cu (2) → Pb-Cu-Zn (3) → Sb-Ag-Au (4) zonations according to the model intrusion-related mineralisations established by Blevin [69].
The mineralisation of the Roc Blanc deposit containing As, Zn, Cu, Pb, Ag, Sb, Au, and Ni can then be related between the two later metallic zones around the BTB granodiorites: Pb-Cu-Zn-As (3) → Sb-Ag-Au ± Ni (4) mineralisation.
According to the structural relationship in the Roc Blanc area, the genesis of the deposits took place after regional schistosity elaboration; then, it occurred between the emplacements of late Variscan BTB granodiorite (327 ± 4 Ma; [30]) and post Variscan microdioritic dykes (255 to 241 Ma, [27,[36][37][38]). After Huon [70] and Huon et al. [71], the ore deposition occurs more probably around 290 Ma as the bulk of mineralisation formed syn-to postpeak greenschist-amphibolite facies metamorphism dated at 290 ± 4 Ma, while for Essaifi et al. [7], the genesis of mineralisation is related to the greisenhydrothermal system associated with emplacement of the Stype leucogranites at 305:9 ± 0:9 Ma, which also lead to the metal zonation around BTB intrusions.
According to our results, any direct link has been established between magmatism events and the deposit genesis (leak of clear magmatic fluids in studied inclusions).

Geofluids
However, Pb and Sr isotopes from ore minerals indicate a magmatic origin of these elements [4], while calculated δ 18 O fluid (15 to 18‰) and C stable isotopes (−7:1 ± 1:2‰) obtained on carbonate mineral gangues are consistent with interactions of metamorphic fluids and organic-rich black shale [4]. Therefore, with respect to these data and the results from fluid inclusion studies presented here, two hypotheses could be set for the formation of the Rock Blanc deposit: (i) The first is the mobilisation of metals by metamorphic fluid (i.e., type 2 H 2 O-CO 2 -CH 4 -N 2 -salt fluid inclusions) from magmatic rocks and/or its surrounding host black shales and the deposition of ore into N-S fractures along the Marrakech Shear Zone at low pressure (below 1 to 1.1 kbar). This fluid was probably mixed with surface to subsurface fluids (i.e., type 3 aqueous fluid inclusions), with eventually cooling and boiling in the system, leading to the deposition of ore (ii) The second hypothesis is the circulation of surface to subsurface fluids (meteoric or brine) as convection cells through magmatic and host metamorphic rocks during retrograde cooling, where they have been enriched in metals and volatile components (CO 2 -CH 4 -N 2 ). This can lead to the formation of immiscible fluids (types 1 and 2). The ore deposition will occur by progressive cooling when mineralised fluid is driven toward the surface. After Pohl [63], this could take place during uplift and extensional tectonics where convection cells are open to the surface. So, this hypothesis will bring the formation of the Roc Blanc deposit near to Triassic basin opening In all cases, the precipitation of the first mineralisation (preores Fe-As stage I, vg 1 , and vg 2 ) is linked to reducing conditions characterised by high temperature (c. 350°±20°C), and the argentiferous mineralisation (main vein-filling stage II; vg 3 , vg 4 , and vg 5 ) precipitated during decreasing temperature (below 200°C), where all the mineralised systems were cooling. The salinity of these two main mineralising fluids, stage I and stage II, is similar with values ranging from 3.7 to 19.4 wt% NaCl equiv. (mean = 13:7 wt%NaCl equiv.) and 4.1 to 19.6 wt% NaCl equiv. (mean = 12:1 wt%NaCl equiv.), respectively. The ore-forming fluids were trapped into N-S fractures, near surface to subsurface conditions at pressure below 1.1 kbar, potentially around 500 bars and at boiling state. According to Chouhaïdi [12] and Benchekroun and Essarraj [11], the precipitation of silver-bearing minerals is likely due to the decrease of salinity, fluid temperature, and ƒO 2 and the increase of fluid pH. However, in our case, just a slight decrease of salinity was observed from vg 2 to vg 4 .
The age of ore deposition remains a main question of debate. According to structural relationship, it may occur around the period of the leucogranite dyke emplacement or more after all Variscan-related magmatic events, during the extensional phase. Our data combined with literature review argues more with a link with late Variscan leucogranite intrusions emplacement, e.g., [7]. In fact, an abrupt increasing temperature in the mineralising system is recorded in the Qtz2 vein (vg 2 ) (Figure 10(d)). This would be related to a rapid heating which can be caused only by a shock temperature of magmatic rock emplacements because the hosting rocks were in retrograde metamorphism sequence. On the other hand, any magmatic fluid has been detected through fluid inclusion studies; therefore, ore forming fluid circulated around leucogranite intrusion when it was cooling. Then, the leucogranite dyke contributed by its temperature and may also provide some metallic elements to ore-forming fluid.
To conclude, a metamorphic fluid would have played a major role for metals collecting from intrusions and their surrounding host organic-rich black shales and their deposition at near surface environment.

Conclusion
The Pb-Zn-Ag-Au deposit of Roc Blanc occurs as N-S to NNW-SSE veins hosted in the contact metamorphism area of the BTB intrusions. The mineralisation was deposited in three main stages (through five generations of quartz-and carbonate-mineralised veins): (1) a preore stage represented by vg 1 and vg 2 , with Fe, As, Zn, and Cu mineralisation (stage I); (2) the main stage composed mainly by Ag, Au, Pb, Zn, Cu, and Sb bearing minerals, observed in vg 3 , vg 4 , and vg 5 (stage II); and (3) the postore stage corresponding to an unmineralised phase made up of chlorite, sericite, and quartz minerals (stage III).
Fluid inclusion studies carried out on quartz-bearing ores (Qtz1, Qtz2, Qtz3, and Qtz4, from vg 1 , vg 2 , vg 4 , and vg 5 , respectively) show the presence of three fluid types which consist of H 2 O-N 2 -CH 4 ±CO 2 -(salt) fluid (type 1), H 2 O-CO 2 -CH 4 -N 2 -(salt) fluid (type 2), and H 2 O-(salt) fluids (type 3). Microthermometric data of the type 2 aqueous gas-bearing fluid inclusions (H 2 O-CO 2 -CH 4 -N 2 -salt) gave a homogenisation temperature ranging from 250°C to 380°C and a relative moderate salinity estimated at 14 wt% NaCl equiv.; this fluid is related to metamorphic origin. For the type 3 aqueous fluid (H 2 O-(salt)), the homogenisation temperatures and salinities give an average value of 150-160°C and 13 wt% NaCl equiv., respectively. These characteristics are consistent with epithermal fluids, originating from surface to subsurface. The type 1 fluid inclusion (H 2 O-N 2 -CH 4 ±CO 2 -(salt)) would result from the mixing of these two latest types which are aqueous-gas-bearing and aqueous gas-free fluids. This process tends to equilibrate the two fluids generating an intermediate fluid composition and temperatures as in the type 1 (T h = 104 − 288°C).
Interpretation of fluid inclusion data in the mineralised veins indicates that fluid is resulted from the mixing of at least two fluids trapped at boiling state. The P-T conditions of ore formation are below 1.1 kbar, and 350 ± 20°C for stage I and c. 150°C for stage II. Their mean salinities are 13.7 and 12.1 wt% NaCl equiv., respectively. Dilution and cooling were the main mechanism of argentiferous ore mineral precipitation in the Roc Blanc deposit, caused by the mixing probably of metamorphic fluid (H 2 O-CH 4 -N 2 -CO 2 ) with surface to subsurface aqueous gas-free fluids (H 2 O-salt, meteoric, or brine).