Abstract

The Roc Blanc 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 (vg1 and vg2). The main ore stage (stage II) consists mainly on Ag, Au, Pb, Zn, Cu, and Sb ores, which is hosted by carbonaceous vein (vg3) and by two late quartz generations veins (vg4 and vg5 with a geodic quartz). Three types of fluid inclusions have been recognized in silver mineralisation bearing quartz veins according to petrographic investigations, microthermometry, and Raman spectroscopy studies: (i) liquid-rich H2O-N2-CH4±CO2-(salt) fluid inclusions (type 1), (ii) vapour-rich H2O-CO2-CH4-N2-(salt) fluid inclusions (type 2), and (iii) aqueous H2O-(salt) fluid inclusions (type 3). The interpretation of fluid inclusion data shows a mixing of two fluids that are metamorphic and surface to subsurface origin, trapped at boiling state. The first mineralised stage was deposited at (this temperature of ore deposition was supported also by chlorite geothermometry) with salinity of 13.7 wt% 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.1 kbar. So, fluid dilution and cooling are probably the main factor for silver deposition in the Roc Blanc polymetallic vein deposit. In addition, fluid inclusion studies reveal that the mineralising fluid corresponds to a mixture of metamorphic fluid (H2O-CH4-N2-CO2) with surface to subsurface aqueous gas-free fluids (H2O-salt, meteoric, or brine).

1. 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].

Mineralisation in the Roc Blanc deposit occurs as veins displaying banded or brecciated structures. It consists mainly of pyrrhotite, pyrite, arsenopyrite, chalcopyrite, sphalerite, galena, Ag-tetrahedrite, Pb-Ag-Sb sulphosalts, argentite, native silver, Ag-Au amalgam, and Ni-Co minerals [4, 5]. These minerals form different paragenesis, and they are precipitated during successive mineralising fluid events [46]. The mineralogy of the Roc Blanc ore corresponds to an intermediate- to low-sulphidation deposit [4].

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, 1018]). 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.

2. Geological Setting

2.1. Central Jebilet Geological Setting

The central Jebilet massif is part of the Moroccan Variscan belt, represented by the Meseta domain and the Palaeozoic blocks of the High Atlas [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.

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 [2630]. 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. 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 [3234]. These plutonic rocks have been dated to (whole-rock Rb-Sr ages; [30]). Two highly differentiated have been recognized in the studied area, and they correspond to a NW-trending, subvertical sheet dikes, dated at (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 [7]. Microdioritic dikes correspond to the youngest intrusive rocks, dated from ca. (K-Ar and 40K-40Ar on kaersutite [3638]) 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, 3941]).

The central Jebilet formations were metamorphosed at low-grade regional metamorphism achieving epizonal to anchizonal conditions [22, 4244]. This metamorphism in the central Jebilet was contemporaneous with post-Visean crustal shortening ductile deformation [45]. Subsequently to the regional metamorphism, a contact metamorphism accompanied the magmatic intrusions emplacement. 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 , , , and have been recognized ([22, 4648] and references therein). was identified as an early synsedimentary event related to extensive-transtension tectonics, allowing the setting up of the Devono-Carboniferous basin [46]. The and deformation phases were related to the closure of the basin. Those episodes begin in Upper Visean-Namurian (~325 Ma). The event was related to the emplacement of late Variscan, syn- to postkinematic granodioritic intrusions, which developed an HT-LP metamorphic aureole.

The Jebilet massif hosts have many metalliferous deposits, mainly represented by massive sulphides, followed by pyrometasomatic deposits and mineralised veins [22]: (i) massive sulphide deposits are the most studied ones. They are base metal pyrrhotite rich (Kettara, Koudiat Aicha, Draa Sfar, etc.) and are differently considered by authors as syngenetic or epigenetic mineralisation (e.g., [29, 39, 40, 44, 49, 50]). Their relative age is Upper Visean-Namurian, which would be reconsidered to Devonian for syngenetic affiliation mineralisation, according to the recent biostratigraphic reconsiderations [25]. (ii) The pyrometasomatic mineralisation (e.g., Frag El Maa) is linked to syntectonic felsic intrusion emplacement [22]. (iii) The mineralised veins are polymetallic and rich in precious metals (e.g., [4, 6, 7]); their emplacement generally postdates the major regional S1 fabrics and before the microodioritic magmatism event which is dated at [36]. The Roc Blanc Pb-Zn-Ag-Au deposit is part of this event.

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 quartz-carbonate 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.

3. 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 microthermometry 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 CO2, CH4, and N2 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 CO2, 2250-2750 cm-1 for N2 and H2S, and 2750-2950 cm-1 for CH4. 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. Ore Descriptions and Typology

4.1. Morphology of Mineralised Veins

The Roc Blanc deposit is a Pb-Zn-Ag-Au mineralisation generally bearing vein structures. It consists of massive open-space fillings of a network subvertical veins, veinlets, and an echelon tension gashes [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 (, , , and ) and N0 to N10 (, , and ) (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 SSE-trending 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)). is the largest one, being more than 2 000 m long with an average thickness of 0.7 m. The veins cut across the cleavage, which strikes N10 to N45 and dips 70°-80°E (Figure 2(b)). Stretching lineations () are also observed; they are oriented toward the NNE direction, along with cleavage. The stratification in the deposit area is generally transposed into cleavage plan; however, some ESE strikes are occasionally observable. At map scale, microdioritic dikes, oriented NNE, also cut across all 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., and ).

4.2. Mineralogy and Paragenesis of the Pb-Zn-Ag-Au Veins

Mineralogical study revealed five vein generations (vg). The first mineralised vein (vg1) is composed of quartz 1 (Qtz1) + chlorite + sericite-muscovite + pyrrhotite ± chalcopyrite ± Fe‐rich sphalerite ± galena (Figure 3(a)). The second vein (vg2) generation is fracture filling quartz 2 (Qtz2) + siderite/ankerite + pyrite + arsenopyrite ore (Figures 3(b) and 3(c)). According to microscopic observations, this vein cut minerals of the first one. Both of these two veins types have been described in central Jebilet, where they are associated with pyrrhotite-rich massive sulphide deposits [50]. The third vein (vg3) is made up of Fe-dolomite + Fe-rich sphalerite + galena + boulangerite + owyheeite + diaphorite + Ag-rich minerals (pyrargyrite, miargyrite, argentopyrite, polybasite, stephanite, Ag-tetrahedrite, argentite, and native Ag; Table 1). The fourth vein (vg4) generates quartz 3 (Qtz3) + chalcopyrite + ullmannite + (Ag-Sb) sulphosalts paragenesis and clearly cross cut the second vein type bearing Qtz2 (Figures 3(c) and 3(d)). The fifth vein (vg5) type contains quartz 4 (Qtz4) + calcite + Fe-poor sphalerite + Ag-tetrahedrite + Ag-Au alloy (Figure 3(d)). Here, Qtz4 occurs as geodic quartz overlapped by secondary pyrite.

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), chlorite, and sericite; it is represented by vg1 and vg2 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 , 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 vg3, vg4 (Qtz3), and vg5 (Qtz4) vein generation minerals. Late sphalerite of this stage is Fe-poor with 3.42 wt% Fe. Textural relationships among silver-rich minerals indicate the following chronological order of precipitation: Ag-tetrahedrite → miargyrite → pyrargyrite → polybasite → stephanite → argentite → native silver. These minerals contain 32 to 35, 36 to 37, 59 to 61, 62 to 63, 65 to 66, 85 to 90, 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

5. 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 vg1-type 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 vg4 with quartz (Qtz3) + chalcopyrite + ullmannite + (Ag-Sb) sulphosalts and geodic quartz vein (vg5) with quartz (Qtz4) + calcite + Fe-poor sphalerite + 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.

5.1. 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 (), microthermometric, and Raman spectroscopy analyses, three main fluid inclusion types are identified in the quartz-mineralised veins (Qtz1 (vg1), Qtz2 (vg2), Qtz3 (vg4), and Qtz4 (vg5)) of the Roc Blanc deposit.

Type 1 and type 2 are two-phase (L and V) aqueous-gas-bearing fluid inclusions and are observed in Qtz1, Qtz2, and Qtz3. However, type 1 is liquid-rich (H2O-N2-CH4±CO2-(salt)) while type 2 is relatively vapour-rich (H2O-CO2-CH4-N2-(salt)). Type 3 is essentially two-phase (L and V) aqueous fluid inclusions (H2O-(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)). (i)Type 1 consists of a liquid-dominated two-phase fluid inclusions with ranging from 10 to 30%. The inclusions are oval or elongated and have regular to subregular shapes with 5 to 15 μm in size. These inclusions are abundant and can represent approximately 40% in single quartz. They form fluid inclusion assemblages (FIA) with other types or can be found isolated (Figures 6(a), 6(b), and 6(e))(ii)Type 2 is composed of two-phase vapour-rich fluid inclusions with an varying from 10 to 50% (av. ). They appear dark, and their size is between 10 and 15 μm. The inclusions have regular shapes and represent approximately 20% of inclusion populations. These inclusions are isolated or occur as clusters throughout quartz crystals. They are closely associated with type 1 and are not all time distinguishable from them by only petrographic observations (Figures 6(a)6(e))(iii)Type 3 consists of two-phase aqueous (liquid + vapour) fluid inclusions showing a bright aspect. These inclusions are generally elongated, regular to subregular with 5-15 μm in size. They have a very consistent low vapour-filling ratio, generally around 5 to 10%. Contrary to types 1 and 2, type 3 inclusions are present in all samples studied. They form generally FIA with type 1 and type 2 in Qtz1, Qtz2, and Qtz3 (Figures 6(a), 6(c), and 6(e)), and they can be dominant with proportion up to 50% or exist alone as in the Qtz4 (Figure 6(f))

5.2. 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% CH4, and 43% N2). Although the volatile phases (CH4, N2, and CO2) are detected in many inclusions, two-phase transitions have been noticed during freezing and heating runs: ice melting temperature and bulk homogenisation temperature () (Table 3). Clathrates form generally in the freezing stage of these inclusion types; however, clathrate melting temperatures were difficult to observe.

5.2.1. Type 1 Fluid Inclusions

These inclusions have been identified in Qtz1, 2, and 3. Their microthermometric data show a variable ice melting temperature. 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 up to -18 and -21.8°C, respectively. During the heating stage, all of these inclusions have been homogenised into the liquid phase. 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, 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 (in mole percent) of gases calculated from Raman spectra areas, in the Qtz1, nitrogen is relatively dominant; ranges from 59 to 77 mol%, and ranges from 23 to 39 mol%. Their average composition is approximately 70 mol% N2 and 30 mol% CH4. In the two inclusions where carbon dioxide is detected, is 2 mol%. In the Qtz2, methane is prevalent, ranges from 52 to 74 mol% with a mean of 60 mol%, and is between 26 and 48 with an average value of 40 mol%. CO2 is missing in all inclusions of type 1 in the Qtz2 or its content is lower than the detection limit. According to the abundance of nitrogen and methane, a relative absence of CO2, and the ice melting temperature variation, the type 1 fluid inclusions are considered to belong to H2O-N2-CH4±CO2-(salt) system.

5.2.2. 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, ranges from -15.8 to -8°C and 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; are between -24 and -12.8°C with an average value of -18°C. Their 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 between -23 and -5°C, , and their 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 CO2, N2, and CH4; 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 CO2 is systematically present in all inclusions of the type 2; however, its content here is relatively dominant (Table 3). According to the microthermometric and Raman spectroscopy data, type 2 fluid inclusions belong to the H2O-CO2-CH4-N2-(salt) system.

5.2.3. Type 3 Fluid Inclusions

They are two-phase (L and V) aqueous bearing fluid inclusions, with very low vapour filling ratio (≈5%) indicating low trapping temperature. In the Qtz1, the ranges from -19.2 to -4.8°C and from 85 to 169°C with mean values of -11.6 and 125°C, respectively (Figures 8(e) and 8(f)). In Qtz2 and Qtz3, have similar values that are between -9.5 and -2.2°C and average values of -6.5 and -6.7°C, respectively. However, their are different and occur at about 160°C in the Qtz2 and 110°C in the Qtz3. In the Qtz4, where only type 3 fluid inclusions have been observed, are variable, ranging from -17.4 to -2.4°C () and are between 74 and 202°C ().

These inclusions are being gas-free; their salinities have then been calculated from ice melting temperatures using the Bodnar equation of state [55, 56]. In the Qtz1, salinities are between 7.9 and 19.4 wt% NaCl equiv. with an average salinity of 16 wt% NaCl equiv. Qtz2 and Qtz3 yield a narrow salinity range (3.7 and 13.4 wt% NaCl equiv.) with a mean around 11 wt% NaCl equiv. For the Qtz4, salinities are between 4.1 and 19.6 wt% NaCl equiv. (). Type 3 fluid inclusions belong therefore to the H2O-(salt) system.

6. Discussion

6.1. Fluid Compositions and Sources

Microthermometric study and Raman spectroscopy analysis show similar fluid inclusion types in three quartz-mineralised veins (Qtz1, Qtz2, and Qtz3) that are H2O-N2-CH4±CO2-(salt) (type 1), H2O-CO2-CH4-N2-(salt) (type 2), and H2O-(salt) (type 3). These similarities indicate the formation of these veins in the same contexts and relatively short interval times. The main fluid inclusion types, forming FIA in Qtz1 (vg1), Qtz2 (vg2), and Qtz3 (vg4) mineralised veins, have not systematically the same compositions. The coexistence of these fluid inclusions with different contents could be explained here by two main hypotheses, which are mixing or boiling. A bivariate plot of (depending on salinities) versus of representative fluid inclusions shows linear distribution of inclusions along axis (Figures 10(a)10(c)). This distribution mode suggests a variation of salinities at nearly constant homogenisation temperature that are characteristic of a boiling process (e.g., [10]). If this process is true for fluid inclusions of the same type in the same mineral, it does not explain the coexistence of all fluid inclusion types due to the different recorder between types 1, 2, and 3 in Qtz1, Qtz2, and Qtz3. However, type 1 fluid inclusions are always plotted between type 2 and type 3 (Figure 10), and results probably from the mixing of these two latest which are aqueous-gas-rich and aqueous gas-free fluids, respectively. This process tends to equilibrate the two fluids, generating an intermediate fluid composition and temperatures as in the type 1 (which is gas-poor, especially in CO2). To resume, the fluid inclusions in the mineralised veins (vg1, vg2, and vg4) would result from the mixing of at less two fluids trapped at boiling state as

Note that CO2 exists only in two inclusions of type 1 in minor proportion (2 mol %).

According with the mixing process and a relative structural time of emplacement of the veins (post cleavages), the source of the aqueous gas-bearing fluid (type 2) can be linked either to retrograde pseudo metamorphic fluids (H2O, CO2, CH4, and N2 mixtures) or to the mesothermal fluid, whereas the aqueous fluid (type 3) is typical for epithermal quartz veins [61, 62]. The type 3 fluid inclusions yield similar salinities and range especially in the Qtz2 and Qtz3 (3.7-13.4 wt% NaCl equiv., 110-160°C, respectively) that are characteristic of epithermal fluids ( and ; [62]). The fluid inclusions of geodic quartz (Qtz4) are typical gas-free and yield also an epithermal characteristic ( and ). For aqueous gas-bearing fluid inclusions which can have metamorphic or mesothermal origin, if the values (250-380°C) argued with these sources, the volatile components (mixtures of CO2, CH4, and N2) and low (-10.8°C) of type 2 fluid inclusions relative to moderated salinities (c. >14 wt% NaCl equiv. estimated from [55]) are not consistent with mesothermal quartz veins. Indeed, after Wilkinson [62], fluids of mesothermal quartz veins are predominantly CO2-N2 or CO2-CH4 mixtures, where intermediate compositions are rare and salinities are low. On the other hand, the contribution of magmatic fluids has not been ruled out for type 2 fluid inclusions. It is possible that these fluids represent magmatic derived or surface to subsurface fluids (meteoric or brine) driven through metamorphic rocks during retrograde cooling, where they have been enriched in volatile components. This could take place during uplift and extensional tectonics, and the resulting convection cells are open to the surface [63]. This can also explain the low recorded, relative to moderate salinities, in addition to the boiling hypothesis for salinity variation. Otherwise, although the type 2 and type 3 coexist in the same FIA (Figures 6(a)6(e)), they display different homogenisation temperatures (e.g., in Qtz2 , 160°C, respectively). This difference can result for fluids trapping at an immiscible state [64] in the aqueous-gas-bearing fluid inclusions (type 2), where reaches 377°C. But these high are consistent with temperature from chlorite geothermometry (). Therefore, the mixing of two fluids with different temperatures is more convincing.

In the type 1 fluid inclusions, the volatile content is dominated by a mixture of N2 and CH4, with minor proportions or absence of CO2, whereas in the type 2, CO2 is abundant and can be dominant, forming a mixture with CH4 and N2 (Figure 10). N2-CH4 mixture indicates reducing conditions which could characterise the mineralisation genesis conditions in vg1 (Qtz1 + chlorite + sericite-muscovite + pyrrhotite ± chalcopyrite ± Fe-rich sphalerite ± galena). This hypothesis is argued by the formation of pyrrhotite which is deemed to form under reducing conditions [65]. A similar fluid composition has been observed in quartz-chlorite-mineralised veins (mineralisation dominated by pyrrhotite) associated to massive sulphide deposit in the central part of Jebilet massif (e.g., Kettara pyrrhotite-rich massive sulphide deposit) where reducing conditions are also advanced [50, 66].

The presence of CH4 in the studied inclusions would derive from the organic matter in the host black shales, explaining the higher proportions of CH4 in the trapped fluids of the area. In the case of N2, according to Bastoul [67], a release of nitrogen occurs by the breakdown of NH4-bearing micas during the prograde metamorphic reactions. The black shales of the central Jebilet are enriched in organic carbon and phyllosilicates, so the production of CH4-N2 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 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.

To constrain the P-T conditions of mineralised vein formations, the batch of isochoric curves was calculated from representative aqueous gas-free fluid inclusions (type 3) using the model of Zhang and Frantz [57] (Figure 11). The Qzt1 (vg1: quartz + chlorite + sericite-muscovite + pyrrhotite ± chalcopyrite ± Fe-rich sphalerite ± galena) and Qtz4 (vg5: geodic quartz + calcite + Fe-poor sphalerite + Ag-tetrahedrite + Ag-Au-alloy) show a large field, giving unrealistic pressure trapping conditions. Contrariwise, the Qtz2 (vg2: quartz + siderite/ankerite + pyrite + arsenopyrite) and Qtz3 (vg4: quartz + chalcopyrite + ullmannite + (Ag-Sb) sulphosalts) yield a narrow field. The intersection of an approximate maximal of type 3 fluid inclusions of Qtz2 (200°C) and Qtz3 (150°C) gives a maximum trapping 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 of its fluid inclusions would represent an upper limit temperature of the Roc Blanc mineralising fluids (). These conditions work with the formation of the first quartz vein generations (vg1 and vg2). For the argentiferous phases, its deposition occurs at c. 100 to 210°C according to the value of aqueous gas-free fluid inclusions of Qtz3 (vg4) and Qtz4 (vg5).

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 increases from Qtz1 (vg1) to Qzt2 (vg2) and decreases progressively towards Qtz3 (vg4) and Qtz4 (vg5) (Figure 10(d)). The same temperature evolution is shown also by type 1 and 2 fluid inclusions in vg1, vg2, and vg4 quartz. These statements can indicate a rapid heating (vg1 to vg2) followed by progressive cooling (vg4 to vg5) 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 (; [30]) and post Variscan microdioritic dykes (255 to 241 Ma, [27, 3638]). 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 , while for Essaifi et al. [7], the genesis of mineralisation is related to the greisen-hydrothermal system associated with emplacement of the S-type leucogranites at , 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). However, Pb and Sr isotopes from ore minerals indicate a magmatic origin of these elements [4], while calculated δ18Ofluid (15 to 18‰) and C stable isotopes () 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 H2O-CO2-CH4-N2-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 (CO2-CH4-N2). 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, vg1, and vg2) is linked to reducing conditions characterised by high temperature (c. ), and the argentiferous mineralisation (main vein-filling stage II; vg3, vg4, and vg5) 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. (.) and 4.1 to 19.6 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 ƒO2 and the increase of fluid pH. However, in our case, just a slight decrease of salinity was observed from vg2 to vg4.

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 (vg2) (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.

7. 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 vg1 and vg2, 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 vg3, vg4, and vg5 (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 vg1, vg2, vg4, and vg5, respectively) show the presence of three fluid types which consist of H2O-N2-CH4±CO2-(salt) fluid (type 1), H2O-CO2-CH4-N2-(salt) fluid (type 2), and H2O-(salt) fluids (type 3). Microthermometric data of the type 2 aqueous gas-bearing fluid inclusions (H2O-CO2-CH4-N2-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 (H2O-(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 (H2O-N2-CH4±CO2-(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 ().

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 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 (H2O-CH4-N2-CO2) with surface to subsurface aqueous gas-free fluids (H2O-salt, meteoric, or brine).

Data Availability

This study constitutes a continuum of the first author’s PhD thesis of Ag-Pb-Zn ores in the Roc Blanc deposit. The studies were carried out at the DLGR Laboratory (Cadi Ayyad University, Marrakech) and at the REMINEX Research Center and Laboratory (Morocco).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

We sincerely thank Research Engineer Marie-Camille Caumon for her assistance in Raman spectroscopy at the University of Lorraine-CNRS-CREGU (Nancy, France). Funding for this project was provided by the Guemassa Mining Company (Managem Group) and the DLGR Laboratory.