Metallogenie, Imiter, Bou Azzer, Bou Madine, Atlas, Maroc
Metallogenie, Imiter, Bou Azzer, Bou Madine, Atlas, Maroc
Metallogenie, Imiter, Bou Azzer, Bou Madine, Atlas, Maroc
Je voudrais enfin remercier mes parents, qui m'ont laissé choisir les
chemins que je voulais suivre fier et fort de leur soutien. Je leur dois
aujourd’hui d’être là et qui je suis.
Merci à mon fafa de fréro, d’être lui même,…. mais non je ne t’en veux
pas pour les playmobiles ! ! Un grand merci à ma soeurette et à François de
leur soutien de tous les jours et d’avoir fait entrer une tornade de bonheur
dans nos vie. A Hugo, mon éternel soleil…
Structural analysis of the Ag-Hg Imiter deposit indicates two tectonic regimes that are
contemporaneous with the silver epithermal episode. Firstly a N-S extensional event associated
with quartz veining phase and rhyolitic intrusions, followed by sinestral movement associated
with a dolomitic phase. U/Pb dating shows two distinct magmatic events at 572 Ma (Taouzzakt
granodiorite) and 550 Ma (Takhatert rhyolitic intrusion). 40Ar/39Ar dating confirms the existence
of two magmatic and additional hydrothermal events related to base metal mineralization and to a
silver epithermal mineralization. Sulphur isotope compositions indicate two sources for this
element, a magmatic one and a metasedimentary one. Re/Os isotopic analyses reveal a mantle
origin of the silver mineralization. The Ag-Hg Imiter deposit is interpreted as a neutral epithermal
deposit, defined by two important parameters, the rhyolitic intrusions and the N-S extension
which controlled the mineralization emplacement.
Mineralization dating of the Bou Madine (553 Ma) and Bou Azzer deposits (533 Ma and 218
Ma) confirms the importance of a metallogenic event during the Precambrian/Cambrian transition
and the existence of potential post-Cambrian mineralizations.
The integration of these results in a regional geodynamic model illustrates the Anti-Atlas
evolution from an active continental margin to a passive continental margin during the
Precambrian/Cambrian transition.
Sommaire
PARTIE I
Chapitre 1 : Introduction générale.....................................................................................1
PARTIE II
3-1-3 Conclusion................................................................................................................85
3-2 Les principaux métallotectes des minéralisations .........................................................87
3-2-1 Contrôle lithologique ...............................................................................................87
3-2-2 Contrôle structural de la minéralisation..................................................................87
3-3 Conclusions sur l'étude structurale ................................................................................88
minéralisation................................................................................................................107
5-1 Problématique et méthode d'analyse............................................................................107
5-2 Discussion et signification des résultats des analyses isotopiques du soufre.............107
5-2-1 Les intrusions rhyolitiques .....................................................................................107
5-2-2 L'épisode tectono-métamorphique .........................................................................109
5-2-3 L'épisode hydrothermale à métaux de base (phase Qtz 1).....................................109
5-2-4 L'épisode épithermale argentifère..........................................................................109
5-3 Discussion et conclusion sur la géochimie du soufre ..................................................111
PARTIE III
géochronologique..........................................................................................................137
1-1 Le cadre géologique régional.........................................................................................137
1-1-1 Le socle métasédimentaire néoprotérozoïque ........................................................137
1-1-2 La couverture volcano-sédimentaire néoprotérozoïque terminale ........................138
1-1-3 La couverture sédimentaire paléozoïque ...............................................................140
1-2 Le gisement de Bou Madine ..........................................................................................141
1-3 Les données isotopiques du Plomb et du soufre ..........................................................141
1-3-1 Les analyses isotopiques du plomb ........................................................................141
1-3-2 Les analyses isotopiques du soufre ........................................................................142
1-4 Les modèles génétiques existants ..................................................................................142
1-5 Les intrusions rhyolitiques ............................................................................................142
1-6 Datation U/Pb de la rhyolite (BM99-1).........................................................................144
Sommaire
géochronologique..........................................................................................................149
2-1 Le cadre géologique régional.........................................................................................149
2-2 La minéralisation............................................................................................................151
2-3 Les paragénèses d’altération .........................................................................................153
2-4 Datation du gisement .....................................................................................................153
2-5 Les modèles génétiques existants ..................................................................................154
2-6 Le filon 7..........................................................................................................................155
2-7 Nouvelles données géochimiques et géochronologiques..............................................157
2-7-1 Analyses isotopiques du soufre ..............................................................................157
2-7-2 Analyses cristallographiques des minéraux d’altération ......................................158
2-7-3 Les datations 40Ar/39Ar de la minéralisation du filon 7 .........................................161
2-7-4 Les datations 40Ar/39Ar du filon de quartz de Tamdrost ........................................161
2-7-5 Datation U/Pb du Trachyte d'Aghbar ....................................................................161
2-8 Conclusion : contribution au modèle génétique de Bou Azzer ..................................163
Références bibliographiques............................................................................................171
Liste des figures
PARTIE I
Chapitre 1
Figure 1.1 Carte de distribution des minéralisations dans l'Anti-Atlas . ................................................................ 1
Chapitre 2
Figure 2.1 Le diagramme concordia (Wetherill, 1956), communément utilisé pour interpréter les analyses
isotopiques du système isotopique U/Pb. ........................................................................................................ 6
Figure 2.2 Diagramme concordia de Tera-Wasserburg, (1972). Les rapports isotopiques ne sont pas corrigés
par rapport au plomb commun....................................................................................................................... 7
Figure 2.3 Procédures de préparation des zircons pour l’analyse isotopique U/Pb à la microsonde ionique..... 10
Figure 2.4 Précautions de montage des grains de zircons sur la bague de résine................................................ 11
Figure 2.5 Schéma simplifié de la microsonde ionique CAMECA IMS 1270 du CRPG-CNRS, Nancy. ............... 13
Figure 2.6 Diagramme de calibration pour la rhyolite de Takhatert (IM99-72) sur une session d'une durée de
vingt-deux heures. Douze analyses de standards ont été effectuées. Une droite de corrélation linéaire est
calculée pour chaque échantillon avec les paramètres de la droite et le coéfficient de corrélation. ............ 15
Chapitre 3
Figure 3.1 Fractionnement isotopique à l’équilibre des sulfures communs en fonction de H2S. ......................... 19
Figure 3.2 Courbes de fractionnement des composés isotopiques du soufre ........................................................21
Figure 3.3 Evolution de la composition isotopique des sulfures et des sulfates au cours des temps géologiques.22
Figure 3.4 Composition isotopique des différents réservoirs terrestres et extra-terrestres .................................. 24
Figure 3.5 Stabilité des mesures de composition isotopique du soufre sur le standard de pyrite (CAR 123) au
cours de période d’analyse de 3 jours........................................................................................................... 27
Chapitre 4
Figure 4.1 Calcul d'un âge modèle d'une roche crustale. La ligne représente l'évolution chondritique du
187
manteau avec un rapport initial Os/188Os de 0.096 et un rapport 187
Re/188Os de 0.402. La ligne en
pointillée représente l'évolution d'une roche crustale ayant un rapport 187Re/188Os de 19.3 et un âge modèle
d'extraction à partir du manteau de 3.5 Ga. ................................................................................................. 31
Figure 4.2 Schémas d'un tube de Carius ----------------------------------------------------------------------------------- 35
Figure 4.3 Descriptif simplifié de la procédure chimique utilisée pour la séparation et l'analyse de l'Os/Re---- 36
Figure 4.4 Principe de la microdistillation -------------------------------------------------------------------------------- 37
Figure 4.5 Schémas du spectromètre de masse Finnigan MAT 262 (TIMS). ........................................................ 38
PARTIE II
Chapitre 1
Figure 1.1 Schéma géologique simplifié de la chaîne de l'Anti-Atlas et localisation des gisements étudiés ------ 41
Figure 1. 2 Interprétations successives des subdivisions du Précambrien dans l'Anti-Atlas. ----------------------- 43
Figure 1.3 Données stratigraphiques et géochronologiques dans l'Anti-Atlas Occidental et central-------------- 45
Liste des figures
Figure 1.4 Données stratigraphiques et géochronologiques dans l'Anti-Atlas Central (Bou Azzer, Ouarzazate) 46
Figure 1.5 Composition de la série volcano sédimentaire du Néoprotérozoïque terminal de l'Anti-Atlas, dans le
diagramme SiO2 vs Na2O+K2O --------------------------------------------------------------------------------------- 47
Figure 1.6 Corrélation stratigraphique et géochronologique (U/Pb sur zircon) de la limite
Précambrien/Cambrien dans l'Anti-Atlas, la Sibérie et Avalone .......................................................................... 48
Figure 1.7 Comparaison des âges des roches volcaniques mésozoïques au Maroc par rapport aux âges
déterminés sur les marges atlantiques de l'Afrique de l'Ouest et des USA.-------------------------------------- 53
Figure 1. 8 Cadre géologique de la boutonnière d'Imiter, modifié. ---------------------------------------------------- 55
Figure 1.9 Composition de la série néoprotérozoïque terminale dans le diagramme TAS -------------------------- 58
Figure 1.10 Successions paragénétiques----------------------------------------------------------------------------------- 62
Figure 1.11 Successions paragénétiques de la minéralisation épithermale argentifère ----------------------------- 62
Figure 1.12 distribution des analyses Pb/Pb dans le diagramme de Zartman et Doe (1981)----------------------- 67
Figure 1. 13 Modèle génétique issu de la synthèse de Baroudi et al. (1999). ----------------------------------------- 68
Chapitre 2
Figure 2.1 Carte géologique et coupe N-S de la partie Ouest de la boutonnière d'Imiter (protrusion rhyolitique
de Takhatert). ................................................................................................................................................ 70
Figure 2.2 A Diagramme TAS de Le Maître (1984) des intrusions rhyolitiques et de la série volcanique
néoprotérozoïque terminal d'Imiter, données compilées Ouguir (1997) et Bajja (1987 ; 1998). ------------- 72
B Diagramme AFM des intrusions rhyolitiques et de la série volcanique néoprotérozoïque terminal d'Imiter,
données compilées Ouguir (1997) et Bajja (1987 ; 1998). -------------------------------------------------------- 72
C Diagramme SiO2 vs Log (CaO / (Na2O-K2O)) des intrusions rhyolitiques et de la série volcanique
néoprotérozoïque terminal d'Imiter, données compilées Ouguir (1997) et Bajja (1987 ; 1998)................... 72
Figure 2.3 A Diagramme Zr/TiO2 vs SiO2 d'après Winchester et Floyd (1977) des intrusions rhyolitiques et de la
série volcanique néoprotérozoïque terminal d'Imiter -------------------------------------------------------------- 74
B Diagramme (Nb/Y )/SiO2 vs SiO2 d'après Winchester et Floyd (1977) des intrusions rhyolitiques et de la série
volcanique néoprotérozoïque terminal d'Imiter --------------------------------------------------------------------- 74
C Diagramme (Nb/Zr)N vs Zr d'après Thieblement et Tergrey (1994) des intrusions rhyolitiques et de la série
volcanique néoprotérozoïque terminal d'Imiter ............................................................................................ 74
Figure 2.4 A Spectre étendu de variations des terres rares dans les intrusions rhyolitiques à Imiter. Normé au
manteau primitif, McDonough et Sun (1995) ----------------------------------------------------------------------- 75
B Spectre de variations des terres rares dans les intrusions rhyolitiques à Imiter. Normé aux chondrites, Evensen
(1978) ............................................................................................................................................................ 75
Figure 2.5 Comparaison des spectres moyens de variations des terres rares dans les intrusions rhyolitiques par
rapport à la série volcanique néoprotérozoïque terminal d'Imiter. Normé aux chondrite, Evensen (1978). 76
Chapitre 3
Figure 3.1 A Carte géologique de surface du gisement d'Imiter I (d'après les documents SMI, modifiés) et régime
de déformation----------------------------------------------------------------------------------------------------------- 78
B Coupe N-S du gisement d'Imiter I déformation (d'après les documents SMI, modifiés) .................................... 78
Liste des figures
Figure 3.2 A Carte géologique de surface du gisement de la colline B3 (d'après les documents SMI, modifiés) et
régime de déformation-------------------------------------------------------------------------------------------------- 80
B Carte de la galerie d'exploitation niveau 1350 m. --------------------------------------------------------------------- 80
C Coupe N-S du gisement de la colline B3 (d'après les documents SMI, modifiés ............................................... 80
Figure 3.3 A Carte géologique de surface du gisement d'Imiter IV (d'après les documents SMI, modifiés) et
régime de déformation-------------------------------------------------------------------------------------------------- 82
B Carte de la galerie d'exploitation niveau 1320 m ---------------------------------------------------------------------- 82
C Coupe N-S du gisement d'Imiter IV (d'après les documents SMI, modifiés)...................................................... 82
Figure 3.4 A Carte géologique de surface du gisement d'Imiter V (d'après les documents SMI, modifiés) et
régime de déformation-------------------------------------------------------------------------------------------------- 84
B Coupe N-S du gisement d'Imiter V (d'après les documents SMI)....................................................................... 84
Figure 3.5 Synthèse des régimes de déformations A et B définis à l'échelle du faisceau de la faille d'Imiter....... 86
Chapitre 4
Figure 4.1 Succession paragénétique de la minéralisation épithermale argentifère (cette étude)........................ 93
Figure 4.2 Diagramme ternaire de la sphalérite-------------------------------------------------------------------------- 96
Figure 4.3 Diagramme de différenciation des sulfoarseniures à cobalt nickel --------------------------------------- 97
Figure 4.4 Diagramme d'évolution Ag vs Hg dans la minéralisation épithermale argentifère............................. 98
Figure 4.5 A Diagramme de caractérisation cristallochimique des phengites ................................................... 100
Figure 4.6 A Diagramme de caractérisation cristallochimique des chlorites----------------------------------------102
B Courbe de densité (histogramme de densité) des températures -----------------------------------------------------102
Figure 4.7 Courbes d'équilibre de différentes réactions entre les sulfures dans un espace as2-T ----------------103
Figure 4.8 Diagramme de stabilité de l'argent à 250°C en fonction du pH et de la fugacité en oxygène........... 105
Chapitre 5
Figure 5.1 Distribution des valeurs de δ34SCDT par phase de minéralisation et par type de sulfures ................. 110
Chapitre 6
Figure 6.1 Diagramme d'évolution du rapport 187Os/188Os en fonction du temps............................................... 115
Chapitre 7
206
Figure 7.1 Diagramme Pb/238U vs 207
Pb/235U des analyses à la sonde ionique CAMECA IMS 1270 sur
monograin de zircon de la granodiorite de Taouzzakt ................................................................................ 121
206
Figure 7.2 Diagramme Pb/238U vs 207
Pb/235U des analyses à la sonde ionique CAMECA IMS 1270 sur
monograin de zircon de la protrusion rhyolitique de Takhatert.................................................................. 122
206
Figure 7.3 Diagramme Pb/238U vs 207
Pb/235U des analyses à la sonde ionique CAMECA IMS 1270 sur
monograin de zircon de l'intrusion rhyolitique de Tachkakacht ................................................................. 123
Figure 7.4 Diagramme de caractérisation géochimique HfO2 vs Y2O3 d'après Pupin (1992) des grains de zircon
de la protrusion de Takhatert et de la granodiorite de Taouzzakt .............................................................. 124
Liste des figures
Figure 7.5 Spectres d'âge des échantillons dans l'auréole de minéralisation et âges en spot analyses dans un
feuillet de l'échantillon IM00-53 ................................................................................................................. 127
Chapitre 8
Figure 8.1 Colonne stratigraphique du Néoprotérozoïque et du cambrien dans le gisement d'Imiter. Les âges de
la colonne de gauche proviennent de la bibliographie (voir Fig. 1.3 et 1.4 pour les références)... ............ 133
Figure 8.2 Nouveau modèle génétique, épithermal neutre, pour le gisement Ag-Hg d'Imiter ............................ 135
PARTIE III
Chapitre 1
Figure 1.1 Carte géologique de la boutonnière de Bou Madine, et localisation des échantillons analysés------138
Figure 1.2 Colonne stratigraphique de la région de Bou Madine. ----------------------------------------------------139
Figure 1.3 Schéma synthétique montrant les relations entre les différents faciès mis en place avant les
minéralisations --------------------------------------------------------------------------------------------------------143
Figure 1.4 Composition chimique d'une rhyolite stérile, d'une rhyolite minéralisée et de la série volcanique
néoprotérozoïque terminale dans le diagramme SiO2 vs Na2O + K2O------------------------------------------144
206
Figure 1.5 Diagramme Pb/238U vs 207
Pb/235U des analyses à la sonde ionique sur monograin de zircon de
l'intrusion rhyolitique de Bou Madine-------------------------------------------------------------------------------145
Figure 1.6 Distribution des valeurs de δ34SCDT de la littérature et les données nouvelles de cette étude à Bou
Madine................................................................................................................................................................. 147
Chapitre 2
Figure 2.1 Carte géologique de la boutonnière de Bou Azzer - El Graara, et localisation des échantillons
analysés. ...................................................................................................................................................... 150
Figure 2.2 Succession paragénétique des minéralisations de Bou Azzer............................................................ 152
Figure 2.3 Modèle génétique de la minéralisation. Succession chronologique des altérations qui ont affecté la
diorite quartzite de Bou Azzer. ---------------------------------------------------------------------------------------154
Figure 2.4 Carte structurale simplifiée et phases de déformation observées au filon 7 ----------------------------155
Figure 2.5 Succession minéralogique du filon 7. Ce filon renferme apparemment des paragenèses correspondant
au stade à arséniures --------------------------------------------------------------------------------------------------156
Figure 2.6 Distribution des valeurs de δ34SCDT de la littérature et les données nouvelles de cette étude à Bou
Azzer ---------------------------------------------------------------------------------------------------------------------158
Figure 2.7 Diffractogramme de Rayons-X de l’échantillon BA99-8 --------------------------------------------------159
Figure 2.8 A Diagramme de détermination cristallochimique des chlorites -----------------------------------------160
B Formule structurale moyenne calculée sur 11 atomes d’oxygène---------------------------------------------------160
C Courbe de densité (histogramme de distribution) des températures------------------------------------------------160
40
Figure 2.9 Datation Ar/39Ar des adulaires de la chloritite associée à la minéralisation à (As Fe Co Ni) du
filon7 et des muscovites du filon de quartz qui recoupe la rhyolite au niveau de la mine de Tamdrost.----161
Liste des figures
Figure 2.10 Diagramme 206Pb/238U vs 207Pb/235U des analyses à la sonde ionique sur monograin de zircon du sill
trachyte d'Aghbar. ....................................................................................................................................... 163
Chapitre 3
Figure 3.1 Nouvelles colonnes stratigraphiques du Néoprotérozoïque et du Cambrien dans l'Anti-Atlas central et
oriental.. ..................................................................................................................................................... 166
Figure 3.2 Scénario schématique de l'évolution géodynamique de l'Anti-Atlas central et oriental du
Néoprotérozoïque terminal au Cambrien inférieur ..................................................................................... 169
Liste des tableaux
PARTIE I
Chapitre 2
Tableau 2.1 Ions analysés à la microsonde ionique pour la datation U/Pb, avec les temps de comptage typiques
et les ordres de grandeur des signaux enregistrés................................................................................ 14
Chapitre 3
Tableau 3.1 Equations de fractionnement des composés isotopiques du soufre. .................................................. 21
Tableau 3.2 Résolution de masse nécessaire pour s'affranchir des interférences isobariques sur les isotopes 32S
et 34S du soufre. Les interférences principales sont soulignées. ........................................................... 26
Tableau 3.3 Les standards internes utilisés.......................................................................................................... 27
Chapitre 4
Tableau 4.1 Caractéristiques de l’Osmium et du Rhénium. .................................................................................. 29
Tableau 4.2 Teneur et rapport isotopiques de différents réservoirs...................................................................... 31
PARTIE II
Chapitre 5
Tableau 5.1 Tableau synthétique des résultats de δ34S........................................................................................ 108
Chapitre 6
Tableau 6.1 Composition isotopique en Os et concentration en Os et Re des échantillons de la minéralisation
argentifère et des intrusions rhyolitiques analysés. (2) : duplicata.................................................... 114
Tableau 6.2 Contribution de l'Os dérivé des black shales dans la minéralisation argentifère en fonction des
rapports isotopiques initiaux calculés à 550 Ma. Le rapport 187Os/188Os initial de l'échantillon IM00-
37, est irréaliste, aussi nous avons utilisé le rapport mesuré pour obtenir une estimation maximale
pour cet échantillon............................................................................................................................ 117
Chapitre 7
Tableau 7.1 Résultats analytiques de la datation U/Pb par sonde ionique des zircons de la granodiorite de
Taouzzakt. Les rapports isotopiques du plomb sont corrigés du plomb commun............................... 120
Tableau 7.2 Résultats analytiques de la datation U/Pb par sonde ionique des zircons de l'intrusion rhyolitique de
Takhatert. Les rapports isotopiques du plomb sont corrigés du plomb commun. .............................. 122
Tableau 7.3 Résultats analytiques de la datation U/Pb par sonde ionique des zircons du dyke rhyolitique de
Tachkakacht. Les rapports isotopiques du plomb sont corrigés du plomb commun .......................... 123
PARTIE III
Chapitre 1
Tableau 1.1 Résultats analytiques de la datation U/Pb par sonde ionique des zircons de l'intrusion rhyolitique de
Bou Maadine. Les rapports isotopiques du plomb sont corrigés du plomb commun. Les rapports U/Pb
sont déterminés en utilisant la corrélation linéaire déterminée sur les standards. ............................ 145
Tableau 1.2 Valeurs de δ34SCDT calculé de la littérature (pyrite, sphalérite, blende) (spectromètre à source
gazeuse, Dagallier, 1989) et les données nouvelles (pyrite) de cette étude (Sonde ionique CAMECA
IMS 3F) à Bou Madine. ...................................................................................................................... 146
Chapitre 2
Tableau 2.1 Valeurs de δ34SCDT de la littérature (pyrite, chalcopyrite, barytine) (Maacha et al., 1998) et les
données nouvelles (pyrite, galène) de cette étude............................................................................... 157
Tableau 2.2 Résultats analytiques de la datation U/Pb par sonde ionique des zircons du Trachyte d'Aghbar. Les
rapports isotopiques du plomb sont corrigés du plomb commun. Les rapports U/Pb sont déterminés en
utilisant la corrélation linéaire déterminée sur les standards............................................................ 162
Partie I Introduction générale et méthodes analytiques
Partie I
Introduction générale
et méthodes analytiques
Partie I Introduction générale
Chapitre 1
Introduction générale
L’Anti-Atlas marocain est une vaste province métallogénique, qui renferme un grand
nombre de gisements métallifères du groupe Ag-Au-Ni-Co-Cu-Pb-Zn dont l’étude de la
distribution dans l’espace et dans le temps demeure la base nécessaire à l’élaboration de
modèles conceptuels indispensables aux programmes d’exploration. Les gisements d’Imiter,
de Bou Madine et de Bou Azzer en particulier ont fait l’objet de nombreuses études depuis
plusieurs décennies visant à la mise au point d’un modèle génétique satisfaisant au mieux les
exigences de la recherche de réserves exploitables. Cependant une certaine confusion existe
toujours aujourd’hui quant à l’interprétation génétique de ces gisements. Dans le but de lever
cette incertitude, deux paramètres géologiques importants doivent être mieux pris en compte,
(i) l’âge absolu de dépôt des minéralisations (déterminations radiométriques), (ii) la source
des métaux (traçages isotopiques), visant à intégrer la distribution géographique des
gisements dans leur contexte géologique et tectonique.
1
Partie I Introduction générale
Dans cet esprit nous avons entrepris de réaliser dans un premier temps, sur le gisement
d'Imiter, une analyse structurale détaillée et le calage de la minéralisation par rapport aux
événements magmatiques néoprotérozoïques. Dans un second temps nous avons réalisé la
40
datation de la minéralisation et des roches associées par les méthodes Ar/39Ar et U/Pb, et
d'identifier les sources des métaux (sulfures et métaux natifs) avec les isotopes du soufre et du
système Rhénium-Osmium. Enfin nous nous sommes intéressés aux gisements de Bou-Azzer
et de Bou-Madine afin de contribuer à la compréhension de la genèse de ces gisements et de
proposer une vision originale de la métallogénèse régionale de l'Anti-Atlas. Un des objectifs
visé est de mieux situer la distribution dans le temps et dans l'espace des métaux précieux au
niveau de la transition Néoprotérozoïque / Paléozoïque.
L'Anti-Atlas est composé de terrains datés du Birimien à l'actuel, témoins d'une
histoire géologique complexe. Dans un tel contexte polyphasé, seul un mode de datation in
situ sur des minéraux particulièrement résistant à l'érosion physique et chimique peut
permettre de déterminer des âges de mise en place. La méthode que nous avons choisi
d'utiliser est la datation U-Pb à la sonde ionique de monograins de zircons.
40
La datation de la minéralisation a été réalisée in situ en Ar/39Ar indirectement sur
monograins de muscovites.
La minéralisation à Imiter est complexe, ce qui rend la détermination de sa source
difficile. C'est pourquoi nous avons décidé de coupler deux types d'analyses isotopiques
distinctes : (i) une étude isotopique in situ du soufre à la sonde ionique afin de pouvoir
analyser les différentes phases de minéralisations et (ii) une étude du couple Re-Os sur
fraction minérale pure
Le plan retenu pour ce rapport final tient compte de ces préoccupations ; il comprend
trois grandes parties :
La première partie consiste en une description des méthodes analytiques utilisées lors
de ce travail.
La seconde partie concerne les résultats obtenus sur le gisement (Ag-Hg) d'Imiter et la
proposition d'un nouveau modèle génétique pour ce gisement.
La troisième partie aborde les aspects géochronologiques et métallogénique des
gisements Co-Ni-As-Au-Ag de Bou Azzer (Anti-Atlas central) et BPGC de l'Ougnat (Anti-
Atlas oriental). A partir de ces études locales nous proposons un modèle génétique original et
de nouveaux guides de prospections à l'échelle de l'Anti-Atlas.
2
Partie I Introduction générale
3
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
Chapitre 2
206
Pb* e λ238t − 1 238
U
Age 207 = x
Pb* e λ235t − 1 235
U
5
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
Les âges obtenus avec les trois géochronomètres peuvent être soit identiques soit
différents. Le premier à avoir observé cette différence, dite aussi discordance, a été Ahrens
(1955), mais ce fut Wetherill (1956) qui l’expliqua par des épisodes successifs de
fractionnement entre U et Pb. Wetherill proposa également le concept de concordia (Fig 2.1).
La concordia est une courbe ou les âges obtenus par les trois méthodes sont identiques, ils
sont dits concordants. Les trois âges ne sont concordants que si le minéral analysé a
fonctionné comme un système chimiquement clos, sans perte ou gain de U ou de Pb.
Figure 2.1 Le diagramme concordia (Wetherill, 1956), communément utilisé pour interpréter les analyses
isotopiques du système isotopique U/Pb.
Lorsque les âges sont discordants, c'est à dire qu'ils ne se positionnent pas sur la concordia, ils
peuvent former un nuage de points ou s’aligner selon une droite qui intercepte la concordia en
deux points. Ces deux points sont appelés intercept supérieur et intercept inférieur. Les âges
6
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
Figure 2.2 Diagramme concordia de Terra-Wasserburg, (1972). Les rapports isotopiques ne sont pas corrigés
par rapport au plomb commun.
commun tend à déplacer les résultats sur une ligne presque verticale, tandis que les pertes de
plomb radiogénique induisent un décalage horizontal. En outre, ce diagramme est bien adapté
238
pour les âges phanérozoïques, puisque le rapport U/206Pb varie beaucoup dans cette la
7
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
gamme d’âge. De ce fait, les deux diagrammes sont complémentaires. De fait, pour les âges
207
récents, le rapport Pb/235U ne varie pas beaucoup, donc le diagramme concordia n'est pas
adapté.
Les minéraux utilisables pour les datations U/Pb, sont le zircon (ZrSiO4), l’apatite
(CaPO4), la baddeleyite (ZrO2), la monazite (CePO4), la titanite (CaTiSiO3), le rutile (TiO2),
la perovskite (CaTiO3), l’allanite (CaFeAl2Si3O12), le xenotime (YPO4). Ces minéraux sont
favorables car ils ne renferment que peu ou pas de très gros cations pouvant être substitués par
Pb2+ et sont riches en terres rares auxquelles se substitue facilement l’Uranium. Le minéral le
plus utilisé pour la méthode U/Pb est le zircon.
8
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
9
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
Figure 2.3 Procédures de préparation des zircons pour l’analyse isotopique U/Pb à la microsonde ionique
L'échantillon est broyé à l’aide d’un concasseur à mâchoires puis d’un broyeur à
cylindres, de façon à obtenir des particules de taille inférieure à 1 mm. Ces particules sont
tamisées en milieu humide avec des tamis métalliques de 50, 200, et 500 µm. La fraction 50-
200 µm subit une séparation par liqueurs denses au bromoforme, puis à l’iodure de
méthylène. La fraction lourde récupérée à l’issue du passage à l’iodure de méthylène est triée
par séparation magnétique (séparateur magnétique Cook). Les échantillons sont triés plusieurs
fois en modifiant les réglages du séparateur magnétique. La fraction finale est non
magnétique. Un tri manuel sous la loupe binoculaire est effectué sur la fraction non
magnétique. Les zircons sont séparés selon des critères physiques et optiques (morphologie,
état de surface, héritages et/ou de zonations, limpidité, couleur, structure, inclusions,
10
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
fractures, etc.). A l’issue de ce dernier tri les grains destinés à être datés, sont montés sur un
bague de résine epoxy®, avec des fragments de zircons standards d 'âge connu. Le standard
utilisé au CRPG de Nancy est le standard 91500, d’âge 1062,4 ± 0,4 Ma (Wiedenbeck et al.,
1995). Pour faciliter le polissage chaque pastille ne contient que des grains de taille
homogène. Les grains sont positionnés le long d'une ligne parallèlement à leur axe C (Fig.
2.4). La pastille obtenue est abrasée de façon à exposer les centres des grains puis polie.
L’étude en microscopie optique des grains de zircons est
Figure 2.4 Précautions de montage des grains de zircons sur la bague de résine.
complétée par des observations au MEB (en électrons retrodiffusés). Cette étude permet
d'identifier précisément les grains ne pouvant pas etre analysés à la microsonde ionique. Pour
l'analyse, la pastille est positionnée dans le porte-échantillon de telle façon que le plus grand
axe du faisceau d'analyse soit parallèle à l'alignement des grains. Cette préparation
particulière permet d'optimiser le positionnement du spot d'analyse sur les grains de très petite
taille (Fig. 2.4). Lors du polissage de grains de petite taille la résine epoxy est abrasée plus
rapidement que le grain. Ce dernier forme un léger relief qui sera amplifié lors des
11
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
12
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
Figure 2.5 Schéma simplifié de la microsonde ionique CAMECA IMS 1270 du CRPG-CNRS, Nancy.
13
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
calibration sont données dans Deloule et al. (2001). Nous y reviendrons plus loin dans la
présentation du traitement des résultats.
Tableau 2.1 Ions analysés à la microsonde ionique pour la datation U/Pb, avec les temps de comptage typiques
et les ordres de grandeur des signaux enregistrés.
206 Pb 206
Pb 206 Pb 206
Pb
=
238
238 238 238
U éch. U std . mesuré U éch. U std . vrai
206
Le rapport Pb/238U mesuré du standard est déterminé pour chaque analyse avec le
rapport 238U16O/238U de l’échantillon, en utilisant les paramètres de la droite de calibration des
standards selon la formule :
14
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
206
Pb
206
Pb 206
Pb 238
U s tan dard
= *
238
U vrai 238
U mesuré UO
+ B
A*
U mesuré
Figure 2.6 Diagramme de calibration pour la rhyolite de Takhatert (IM99-72) sur une session d'une durée de
vingt-deux heures. Douze analyses de standards ont été effectuées. Une droite de corrélation linéaire est
calculée pour chaque échantillon avec les paramètres de la droite et le coefficient de corrélation.
15
Partie I Chapitre 2 La méthode de datation U/Pb à la sonde ionique
indiqué de réaliser des analyses de standard le plus régulièrement possible (au moins tous les
trois échantillons) afin de bien contraindre la droite de calibration et d'identifier le cas échéant
de possibles variations de celle-ci.
16
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
Chapitre 3
3-1. Généralités
Le soufre possède quatre isotopes stables dont les abondances sont :32S = 95.0% ; 33
S
= 0.76% ; 34S = 4.22% ; 36S = 0.014% (Dean, 1979). Les erreurs analytiques sont petites par
rapport aux variations dans les milieux naturelles (< ‰ ; > %) . Les analyses sont effectuées
soit par spectromètre de masse (MS à source gazeuse) soit par sonde ionique. Les
34
compositions isotopiques sont toujours exprimées sous forme d’un quotient S/32S. La
différence entre le quotient de l’échantillon (Qx) et celui du standard (Qs) est exprimée sous
forme de δ en ‰, suivant la formule :
Q
δ 34S ‰ = x − 1 *1000
Qs
La valeur de δ peut être soit positive soit négative, ce qui signifie respectivement un
enrichissement ou un appauvrissement de l’isotope lourd 34S dans l’échantillon par rapport au
standard. Le standard international est la troilite (FeS) de Canyon Diablo (CDT), dont le
quotient 34S/32S est de 0.044994 (Thode et al., 1961).
La composition isotopique du soufre varie dans le milieu naturel. Cette variation est
due en partie au fractionnement isotopique lors des réactions chimiques. Nous allons présenter
successivement les grandes lois qui régissent le comportement isotopique du soufre et les
applications qui en découlent.
17
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
composant dont les énergies de liaisons atomiques sont les plus élevées. Le facteur de
fractionnement isotopique (α) entre deux composés, A et B, est égal au quotient des isotopes
lourds sur les isotopes légers (R) des composés A et B. Dans les conditions d’équilibre, le
facteur de fractionnement et le quotient (Q) sont liés par la constante d’équilibre (K)
1+ δ A
α A− B =
RA
= 1000 = 1000 + δ A
RB 1 + δ B 1000 + δ B
1000
(3)
Les mesures d’abondance des isotopes stables sont données sous forme de quotient (Q) ; il
faut donc tenir compte des déviations du facteur de fractionnement (α) par rapport au standard
suivant l’équation, ci-dessous, qui dérive des équations (1) et (2).
La valeur α est toujours proche de 1 (1 ± 0.0x à 1 ± 0.00x). La différence de entre
deux composés soufrés est notée ∆. La valeur ∆ est similaire à la valeur de la transformation
logarithmique du facteur de fractionnement associé.
∆ A− B = δ A ‰ - δ B ‰ ≈ 1000lnα A -B
(4)
Afin d’illustrer les équations théoriques précédentes, prenons comme exemple réel la réaction
d’équilibre isotopique entre le dioxyde de soufre (SO2) et le sulfure d’hydrogène gazeux
(H2S):
32
SO2 + H234S = 34SO2 + H232S (5)
18
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
Figure 3.1 Fractionnement isotopique à l’équilibre des sulfures communs en fonction de H2S.
D’après Ohmoto et Rye, (1979).
19
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
3-3-1. La géothermométrie
Les équations de fractionnements isotopiques en fonction de la température des
minéraux les plus communs dans les gisements hydrothermaux sont présentées Tableau 3.1 et
Figure 3.2. En-dessous de 350°C les réactions entre sulfates et sulfures sont particulièrement
lentes et l’équilibre isotopique n’est jamais atteint. Certains paramètres comme un temps de
résidence long, un refroidissement lent, un pH bas et un fluide riche en soufre favorisent les
réactions d’équilibres isotopiques (Ohmoto et Rye, 1979 ; Ohmoto et Lasagna, 1982). Les
différents couples isotopiques de sulfures ne réagissent pas de la même façon selon les
conditions du milieu (Ohmoto et Rye, 1979). La forme générale des équations de
géothermométrie est:
Cette équation indique que les effets du fractionnement donnés par 1000lnα A-B, varient avec
l’inverse du carré de la température (en degrés Kelvin). Les valeurs positives ou négatives de
1000lnα montrent respectivement un enrichissement ou un appauvrissement des isotopes
lourds dans le premier membre du couple. Un fractionnement isotopique important caractérise
la réaction entre SO42- et H2S. Ce fractionnement est du à la différence d’état d’oxydation
entre SO42- et H2S. L’isotope lourd, 34
S, enrichit préférentiellement le membre le plus oxydé
du couple (Fig. 3.2, courbe 1 et 2).
20
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
Figure 3.2 Courbes de fractionnement des composés isotopiques du soufre (Ohmoto et Rye,
1979 ; Ohmoto et Lasaga, 1982)
21
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
Les météorites sont les objets qui présentent les compositions isotopiques de soufre les
plus homogènes, avec une moyenne de δ34SCDT = 0.2 ± 0.2 ‰ (Thode et al., 1961). Cette
Figure 3.3 Evolution de la composition isotopique des sulfures et des sulfates au cours des temps géologiques.
(Schidlowski et al., 1983)
22
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
homogénéité conduit de nombreux auteurs à supposer que le δ34SCDT moyen terrestre est
probablement très proche de (δ34SCDT ≈ 0 ‰) (Thode et al., 1961). Cette hypothèse a été
confirmée par l’analyse de roches d’origine mantellique (MORB). Schidlowski et al. (1983)
montrent que les différents réservoirs crustaux ont commencé à s’éloigner de δ34SCDT ≈ 0 ‰,
il y a 2,7 Ga d’années (Fig. 3.3). Cet âge coïncide approximativement avec le développement
de la vie sur terre et donc l’apparition probable des fractionnements biogéniques associés à la
réduction bactérienne des sulfates en sulfures. Cette hypothèse est bien illustrée par les études
sur le Précambrien de Gorjan et al. (2000). En effet, ces auteurs mettent en évidence les
relations entre le taux de réduction des sulfates et les cycles glaciaires (diminution de la vie).
Claypool et al. (1980) montrent l'importante variabilité des valeurs de δ34SCDT des sulfates
marins (δ34SCDT ≈ 10 à 30 ‰) depuis le Phanérozoïque. Leur valeur actuelle est de
δ34SCDT ≈ 20 ±1‰.
Les sulfures produits par les émanations volcaniques/géothermales sont en général
isotopiquement homogènes. Mais ils peuvent présenter des variations remarquables lors de
succession de séquence redox (SO2 ≈ -8 à +18 ‰ ; soufre natif ≈ -15 à +16 ‰ ; H2S ≈ -9
à 6 ‰). Les valeurs de δ34SCDT des intrusions basiques sont en général proches de ≈ 0 ‰
(gabbros de Sudbury et de Stillwater ; Thode et al., 1962). La dispersion des valeurs de
δ34SCDT des sulfures dans certaines intrusions basiques (δ34SCDT ≈ -6 à +10 ‰) est expliquée à
la fois par de la contamination par des sources sédimentaires et par les réactions redox dans la
chambre magmatique. Les sulfures sédimentaires d’origine diagénétique - syngénétique ou
épigénétique ont une distribution des valeurs de δ34SCDT extraordinairement large (δ34SCDT ≈ -
50 à 50 ‰). Ces variations sont causées par le fractionnement chimique entre sulfates et
sulfures (δ34SCDT ≈ - 20 à -50 ‰). Cette dispersion est accrue lorsque se surimpose le
fractionnement lié à la réduction biogénique de SO4-2 en H2S (Ohmoto and Rye, 1979). Si
34
l’appauvrissement en S des sulfures est considéré classiquement comme un effet des
processus biogéniques, l’enrichissement en 34S est associé à des facteurs tel que la source, le
réservoir de sulfates, et les conditions physico-chimiques dans le système hydrothermal.
En conclusion, les grandes différences de δ34SCDT des principaux réservoirs terrestres
(Fig. 3.4) font du soufre un élément potentiellement très puissant pour l'identification de la
source de cet élément.
23
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
Figure 3.4 Composition isotopique des différents réservoirs terrestres et extra-terrestres (Field et Lombardi,
1972 ; Field et Gustafson., 1976 ; 1983 ; 1984 ; Ohmoto et Rye, 1979 ; Claypool et al., 1980 ; Hoefs, 1980 ;
Sakai et al., 1984).
la réduction de sulfates marins. Dans les gisements épithermaux les sulfures précipitent dans
des conditions de pH et de fO2 où H2S est l'espèce dominante. Quand les fluides se sont
24
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
équilibrés avec les roches encaissantes, le pH et/ou la fO2 augmentent, ce qui engendre la
précipitation de sulfure avec des valeurs plus basses (négatives). L’évidence d’une source
magmatique pour les sulfures à ≈ 0‰ est une généralisation dangereuse dans les
environnements épithermaux. Cette généralisation n’est vraisemblable que si la quantité de
sulfures précipités permet de supposer que le système n’est composé que de soufre réduit
(H2S).
Pour cette étude, nous avons utilisé différents moyens d’analyse. La majeure partie des
résultats présentés a été obtenue sur la sonde ionique CAMECA IMS 3F. Pour des problèmes
d’interférence de masse, les analyses sur la blende ont été réalisées avec la sonde ionique
CAMECA IMS 1270. Enfin des analyses par spectrométrie de masse à ablation laser ont été
réalisées à l’Université de Queen’s par Kurt Kyser. Dans cette partie nous ne détaillerons que
les procédures que nous avons personellement réalisées, c’est à dire les analyses par sonde
ionique CAMECA IMS 3F et CAMECA IMS 1270.
25
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
Tableau 3.2 Résolution de masse nécessaire pour s'affranchir des interférences isobariques sur les isotopes 32 et
34 du soufre. Les interférences principales sont soulignées.
Pour cette étude les isotopes du soufre sont mesurés en mono-collection par multiplicateur
d'électron. Le faisceau primaire d’ions d’oxygène a un diamètre de 50 µm. Le faisceau
secondaire n’est pas filtré en énergie. Le fractionnement instrumental est de -63.6 ± 0.9 ‰, -
66.6 ± 0.5 ‰, -70.7 ± 1.1 ‰ respectivement pour les standard internes de pyrite,
d'arsenopyrite et de galène. Tous les résultats sont présentés en δ34SCDT ‰. La reproductibilité
26
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
des analyses sur le standard de pyrite, d'arsenopyrite et de galène est inférieure à 1.2 ‰. La
tension du faisceau primaire est de -10 kV. L’intensité sur l’échantillon est de 40 nA. Les
temps de comptage de 32S et le 34S sont respectivement de 3 et 10 secondes. La séquence de
mesure est répétée 100 fois par bloc de dix cycles chacun. Le bruit de fond mesuré à la masse
31.5 est compté pendant 50 sec au début de chaque analyse. La valeur du bruit de fond est
toujours inférieur à 0,1 coups par seconde, et sa valeur moyenne lors d’une session de 5 jours
(soit 200 analyses de standard et d’échantillon de nature diverse) est de 0.015, avec un écart
type de 0.025.
La liste des standards utilisés ainsi que la valeur de leur δ34SCDT sont présentées dans
le tableau 3.3. Pour calibrer les mesures, une série de standards est analysée au début et à la
fin de chaque séance d’analyses ainsi qu’à chaque changement d’échantillon. Le quotient
entre le nombre de standards et le nombre d’échantillons analysés est proche de 3. Dans des
conditions normales, la reproductibilité sur les standards est de l’ordre de 1‰ (Fig. 3.5). Sur
une période de plusieurs jours la reproductibilité varie un peu plus (1.5‰) du fait des
instabilités instrumentales.
Figure 3.5 Stabilité des mesures de composition isotopique du soufre sur le standard de pyrite (CAR 123) au
cours d'une période d’analyse de 3 jours.
27
Partie I Chapitre 3 La méthode d'analyse des isotopes du soufre
Les analyses isotopiques de blende nécessitent une résolution de masse de 6300 pour
s'affranchir des interférences isobariques sur les isotopes 68Zn2 / 34S. La sonde ionique IMS 3F
ne permettant pas d'obtenir une telle résolution, aussi les analyses ont été réalisées à la sonde
ionique CAMECA IMS 1270. Les conditions d’analyses sur la sonde ionique CAMECA IMS
1270 sont différentes de la sonde IMS 3F. Le courant primaire, associé au canon à électron,
est de +10kV et l’intensité sur l’échantillon est de 20 nA. Le diamètre du spot d’analyse est de
10 µm. Les ions sont analysés avec la cage de Faraday. Le faisceau secondaire n’est pas filtré
en énergie. Les isotopes analysés sont le soufre 32S et 34
S simultanément en 10 cycles de 10
secondes. Un nettoyage de la surface d'une durée de 5 minutes est effectué avant chaque
analyse.
28
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
Chapitre 4
29
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
La propagation des erreurs se fait à partir des erreurs analytiques des rapports mesurés et de la
constante de désintégration (Sambridge et al., 1997)
187
(ii) la construction d’une isochrone avec les rapports Re/187Os vs 187
Os/188Os, issus de
roches ou de minéraux cogénétiques. L'âge est déterminé à partir de la pente de la droite
passant par l'ensemble des points (régression des moindres carrés), en utilisant la relation : m
= eλt-1. L'intersection de la droite et de l'axe des ordonnées donne directement le rapport
isotopique initial 187Os/188Os. Si le MSWD calculé est faible (<2), les incertitudes sur l'âge et
la composition isotopique initiale sont uniquement contrôlées par les erreurs analytiques. Un
MSWD fort suggère que le système isotopique est perturbé (métamorphisme,
30
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
Figure 4.1 Calcul d'un âge modèle d'une roche crustale. La ligne représente l'évolution chondritique du
187
manteau avec un rapport initial Os/188Os de 0.096 et un rapport 187
Re/188Os de 0.402. La ligne en pointillée
représente l'évolution d'une roche crustale ayant un rapport 187Re/188Os de 19.3 et un âge modèle d'extraction à
partir du manteau de 3.5 Ga.
Tableau 4.2 Teneur et rapport isotopiques de différents réservoirs d'après Shirey et Walker (1998).
31
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
4-2-1. Le noyau
Le noyau métallique terrestre est considéré comme le réservoir majeur d'Os et de Re
sur Terre (Allègre et Luck, 1980 ; Brandon et al., 1998). Les analyses effectuées sur des
météorites métalliques permettent d'estimer la concentration de Re et d'Os dans le noyau à une
valeur de l'ordre du ppm (10-6 g/g) (Luck et al., 1980 ; Horan et al., 1992 ; Morgan et al.,
1995; Shen et al., 1996). Cette valeur se différencie nettement de celle du manteau ou de la
croûte, dont la concentration est de l'ordre respectivement du ppb (10-9 g/g) ou du ppt (10-12
g/g).
32
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
Turekian, 1992 ; Horan et al., 1994 ; Pierson-Wickmann et al., 2000). La forte teneur en Re
187
de ces formations engendre de l'Os radiogénique, ce qui augmente le rapport Os/188Os.
Nous verrons une application de ce principe dans notre étude.
33
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
La méthode utilisée dans cette étude est la mesure faite sur un spectromètre de masse à
thermo-ionisation en mode négatif (N-TIMS). C'est actuellement la méthode la plus fiable et
la plus précise pour déterminer des rapports isotopiques de l'osmium (Creaser et al., 1991) et
les teneurs en rhénium et en osmium . L'attaque acide de l'échantillon se fait par la méthode
du tube de Carius et l'extraction de l'osmium et du rhénium suivant la méthode de Birck et al.
(1997).
34
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
35
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
Figure 4.3 Descriptif simplifié de la procédure chimique utilisée pour la séparation et l'analyse de l'Os/Re selon
de Birck et al (1997)
brome, alors que le Re reste dans la phase aqueuse. Une fois refroidi, l'Os contenu dans le
brome est extrait de la solution et mis dans un bêcher contenant de l'eau distillée pour éviter
qu'il ne s'évapore. 1cc de Br2 est ajouté à la solution initiale et cette phase de
chauffage/extraction est répétée une seconde fois. Une fois la seconde extraction réalisée, l'eau
distillée recouvrant le brome extrait est enlevée et 0.5 cc de HBr sont ajoutés. La solution
obtenue est évaporée pour subir la microdistillation (Roy Barman, 1993). Le principe de la
36
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
microdistillation est schématisé sur la figure 4.4. Le résidu d'évaporation placé dans le
couvercle d'un bécher conique, est recouvert d'une goutte de solution de CrO3. Une demi-
goutte de HBr est déposée dans le fond du bêcher conique. La présence de CrO3 permet
l'oxydation de OsBr62- en OsO4 qui est volatil à basse température. Le bêcher est entouré de
papier aluminium pour favoriser la condensation de l'Os au niveau de la goutte de HBr. La
distillation se fait à 80°C pendant 3 heures. La goutte de HBr est alors séchée jusqu'à un
volume d'environ 2 ml. L'Os est alors prêt pour l’analyse au spectromètre de masse.
37
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
38
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
contribution du rhénium. La solution la plus simple et la plus fiable est d'avoir un rapport
185
Os/188Os mesuré le plus faible possible.
Les valeurs obtenues au spectromètre de masse sont corrigées (1) des interférences
isobariques des molécules contenantes isotopes lourds de l'oxygène (17O, 18
O) sur les
molécules contenants les isotopes légers de l'osmium ; (2) du fractionnement de masse ; (3) de
la contribution du traceur.
(1) L'abondance de 186Os16O3- est estimée à partir de la mesure sur les masses 235 (186Os16O2
I7
O- + 187
0s16O3-) et 236 (186Os16O218O- + 187
Os160217O- + 188
Os16O3-). La valeur du rapport
naturel 188Os/186Os est estimée à 8.30910 (Luck et Allègre, 1982). L'équation suivante permet
d'estimer l'intensité de 186Os16O3- :
( ) −( )
17
188 187 O
Os16O3 M
Os16O3 M
x3 x 18
186
Os16O3− = O
2
188Os 18
O 17O
186 + 3 x 16 − 9 x 16
Os R O O
où M représente les intensités mesurées aux masses 236 et 235 et R le rapport 188Os/186Os réel
17
Les valeurs des rapports O/16O et 18
O/16O sont respectivement 0.0003708 et 0.002045
(Creaser et al. (1991), d'après Nier (1950)
Les corrections les différents pics se font selon les équations ci-dessous :
( ) ( )
17
O
187
Os16O3− = (masse 235) M − 186
Os16O3− R x3 x 16
O
( ) ( ) ( )
17 18
O O
188
Os16O3− = (masse 236) M − 186
Os16O3− x 3 x 16
− 186
Os16O3− R x3 x 16
R
O O
( ) ( ) ( )
17 18
O O
189
Os16O3− = (masse 237) M − 188
Os16O3− x 3 x 16
− 187
Os16O3− R x3 x 16
R
O O
( ) ( ) ( )
17 18
O O
190
Os16O3− = (masse 238) M − 189
Os16O3− x 3 x 16
− 188
Os16O3− R x3 x 16
R
O O
( ) ( )
17
O
192
Os16O3− = (masse 240) M − 190
Os16O3− x 3 x 16
R
O
192
(2) Le fractionnement de masse est corrigé avec le rapport Os/188Os (rapports de masses
240/236) qui est égal à 3.08271 (Luck et Allègre, 1982). Les rapports individuels de chaque
39
Partie I Chapitre 4 La méthode d'analyse du système Rhénium / Osmium
(3) Le traceur isotopique d'Os n'est pas une solution pure d'190Os. Les résultats déterminés sur
les échantillons sont corrigés de la participation du traceur. Les concentrations d'osmium sont
calculées à partir du rapport 192Os/190Os.
Les standards sont analysés en début de session de mesures sur 5 blocs de 10 mesures.
Les dépôts se font comme pour les échantillons. Le rapport isotopique d'Os a une valeur
moyenne de 0.17393 ± 0.00038 (2σ) sur la période de l’étude.
1.667
FDM =
187 Re
185
Re std
La précision de ces rapports, basée sur la statistique de comptage, varie entre 0.15% et 1.31%
pour les rapports mesurés (2σm). Le rapport 187Re/185Re corrigé de l'échantillon est obtenu en
multipliant le rapport 187Re/185Re brut par le FDM. Le rapport 187Re/185Re corrigé est introduit
dans l'équation de dilution isotopique pour calculer les concentrations de rhénium.
40
Partie II Etude géologique et isotopique du gisement Ag-Hg d'Imiter
Partie II
Chapitre 1
Avant de présenter la géologie des différents gisements étudiés dans le détail, ce qui
permettra de mettre en évidence les variations locales, il convient de rappeler les grands traits
de la géologie de l'Anti-Atlas.
Le Maroc est subdivisé en cinq grands domaines géologiques du Nord au Sud, le Rif,
la Meseta, l’Atlas, l’Anti-Atlas et le Sahara. La chaîne de l’Anti-Atlas fait partie de la ceinture
orogénique panafricaine. Elle s’étend sur 750 km selon une direction WSW-ENE, de
l’Atlantique, où elle prolonge la chaîne de Zemmour, jusqu’au Tafilalet. L’Anti-Atlas
constitue le domaine structural majeur du Sud marocain. Il est limité au Sud par le bassin
carbonifère de Tindouf. D'un point de vue géographique cette limite est matérialisée par la
Figure 1.1 Schéma géologique simplifié de la chaîne de l'Anti-Atlas et localisation des gisements étudiés
41
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
42
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Figure 1.2 Interprétations successives des subdivisions du Précambrien dans l'Anti-Atlas d'après Hassenforder
(1987), complété.
43
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
44
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Figure 1.3 Données stratigraphiques et géochronologiques dans l'Anti-Atlas occidental et central (d'apres
Hassenforder, 1987 ; modifié)
45
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
d'aulacogènes fonctionnant indépendamment les uns des autres et dont le degré d'évolution
diffère (Fekkak et al., 1999).
Figure 1.4 Données stratigraphiques et géochronologiques dans l'Anti-Atlas Central (Bou Azzer, Ouarzazate)
46
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Figure 1.5 Composition de la série volcano sédimentaire du Néoprotérozoïque terminal de l'Anti-Atlas, dans le
diagramme SiO2 vs Na2O+K2O de Le Maître et al., 1989. Données compilées de Garcia (1979) ; Chikaoui
(1981) ; Gyssel, (1982) ; Bodinier et al. (1984) ; Baouch (1984) ; Lasri (1984) ; Bajja (1987, 1998) Freton,
(1988) ; Abia (1991) ; Ait Saadi (1992) ; Ouguir (1997).
continentale, à la fin de la collision panafricaine (Bajja, 1998 ; Ouguir, 1997 ; Aït Saadi,
1992) (Fig. 1.4). Des essais de modèles pétrogénétiques réalisés par Aït Saadi, (1992) dans
47
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
48
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
49
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
50
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
de la base au sommet, les âges déterminés sont de 534 ± 10 Ma pour les trachytes du Jbel
Boho interstratifiés dans les dolomies inférieures à la base de la série Lie de vin (Ducrot et al.,
1977 ; U/Pb par dilution isotopique sur population de zircons), 522 ± 2 Ma et 521 ± 7 Ma
pour la seconde partie de la série Lie de Vin (Landing et al., 1998 ; U/Pb sur monograins de
zircon par dilution isotopique ; Compston et al., 1992 ; U/Pb sur monograins de zircon à la
microsonde ionique SHRIMP) et 517 ± 2Ma pour le sommet des calcaires supérieurs
(Landing et al., 1998 ; U/Pb sur monograins de zircon par dilution isotopique), qui représente
la limite entre le Cambrien inférieur et le Cambrien moyen. Ces âges sont en accord avec les
observations sur les microfossiles et les analyses d'isotopes stables réalisées respectivement
par Lathan et al., 1990 et Tucker, (1986), qui placent la limite Précambrien/Cambrien sous le
contact entre dolomie inférieure et la série Lie de Vin.
51
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
52
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Figure 1.7 Comparaison des âges des roches volcaniques mésozoïques au Maroc par rapport aux âges
déterminés sur les marges atlantiques de l'Afrique de l'Ouest et des USA.
Ages K/Ar sur roche totale : -Leblanc, 1973 : Foum Zguid ; -Hailwood et al., 1971 : Foum Zguid, Bas
Draa ; -Schermerhorn et al., 1978 : Foum Zguid, Messajana dyke ; -Westphal et al., 1979 : Meseta, Atlas Nord,
Atlas Sud ; -Conrad et al., 1975 : Zemour, Reggane ; -Manspeizer et al., 1978 : Atlas, New Jerzey, Connecticut,
Nova Scotia ; -Dupuy, 1988 : Libéria, Nova Scotia
Ages Ar/Ar sur plagioclase : -Sebai et al., 1991: Messajana, Taoudenni, Haut Atlas, Foum Zguid. ;
-Fiechtner et al., 1992 : Haute Mouloudia, Haut Atlas.
Ages U/Pb sur zircon par dilution isotopique : -Dunning et al., 1990 : Palissade
Est et Ouest de l'Atlantique (Fig. 1.7). Sebai et al., (1991) et Fiechtner et al., (1992) ont daté
40
la base et le sommet des volcanites du Haut-Atlas (Maroc) par la méthode Ar/39Ar sur
plagioclase. Les résultats obtenus sont respectivement de 201 Ma et 210 Ma pour la base et de
197 Ma et 196 Ma pour le sommet. Le dyke de Foum Zguid est daté par la même méthode à
196 Ma (Sebai et al., 1991).
53
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
sur la route reliant Errachidia à Ouarzazate (Fig. 1.1). La mine se trouve à une altitude de
1500 m ou règne un climat semi-désertique.
54
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
où la production d'argent n'est pas associée à celle d'autres métaux de base tels que le cuivre,
le zinc ou le plomb.
Figure 1.8 Cadre géologique de la boutonnière d'Imiter, modifié. (Leistel et Qadrouci, 1991)
55
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
56
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
A partir des analyses réalisées par Camara (1993) sur les amphiboles des différents
massifs, nous avons déterminé la pression de cristallisation de ces minéraux à l'aide des
formules définies par Schmid (1992). Les pressions minimales déterminées sur les massifs de
Bou teglimt et d'Igoudrane sont respectivement de 3 kb et 5 kb.
57
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Figure 1.9 Composition de la série néoprotérozoïque terminale dans le diagramme de Le Maître et al. (1989).
Données compilées de Bajja (1988) et Ouguir (1997).
58
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
1997). La série néoprotérozoïque terminale est interprétée d'un point de vue géodynamique,
dans contexte complexe associant marge continentale active, et province volcanique, ce qui
expliquerait la coexistence d’un magma crustal et d’un magma mantellique peu évolué.
Une partie des dykes basiques à intermédiaires recoupe les formations du Paléozoïque.
59
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Imiter I-Grande Carrière est le gisement le plus ancien. Il regroupe la carrière des
anciens et la grande carrière. Les corps minéralisés sont tous encaissés dans la formation
métasédimentaire du Néoprotérozoïque moyen. Leur extension horizontale est de 800 m pour
une extension verticale de 250 m. La teneur moyenne oscille entre 800 et 1200 g/t Ag.
Colline B3 est le gisement qui présente la teneur moyenne la plus élevée, 2500 g/t Ag.
La minéralisation est encaissée dans la formation métasédimentaire du Néoprotérozoïque
moyen ainsi que dans le conglomérat de base et les tufs de la formation volcano-sédimentaire
du Néoprotérozoïque terminal. L'extension latérale de la minéralisation est de l'ordre de 350
m pour une extension verticale reconnue de 300 m.
Plusieurs études paragénétiques détaillées ont été réalisées par Vargas (1983), Guillou
60
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
et al. (1985, 1988), Popov et al. (1986), Popov (1995), Baroudi et al. (1999). Les résultats
sont assez discordants comme le montrent les deux tableaux suivants (Fig. 1.10-11).
Néanmoins nous effectuerons dans ce paragraphe introductif une tentative de synthèse des
deux tableaux de Guillou et al. (1985) et Vargas (1983) et de Popov (1995).
61
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Figure 1.10 Successions paragénétiques selon Vargas (1983) et Guillou et al. (1985)
Figure 1.11 Successions paragénétiques de la minéralisation épithermale argentifère selon Popov (1995).
62
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
1-9-4. L'épisode sulfuré hydrothermal précoce (Guillou et al.,, 1988 ; Vargas, 1983 ;
Popov, 1995)
Il est caractérisé par des fentes de quartz en mosaïque et chlorite. Ces fentes recoupent
les minéralisations syngénétiques et métamorphiques et sont elles-mêmes recoupées par les
stockwerks à minéralisation argentifère épithermale. La paragenèse est composée de pyrite en
plages xénomorphes corrodées par les autres sulfures (Ag 0,18-0,34%, Hg 0,48-1,22%, As
0,86%), de sphalérite à inclusion de chalcopyrite (Fe 3,32-7,8%, Hg 0,2-0,6%, As 0,2%, Ni
0,25%, Bi traces), de galène (Ag 600 ppm à 0,25%, Bi 0,6-1,7% Hg 0,46%) et d'arsénopyrite
en cristaux automorphes.
63
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
64
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Ce paragraphe résume les résultats de l'étude de Baroudi et al. (1999). Ces auteurs ont
travaillé sur deux types de quartz, qu'il associe respectivement aux phases de minéralisation
structuro-métamorphique ou tectono-métamorphique et à gangue dolomitique ou épithermale
argentifère. Le fluide associé à la première phase (structuro-métamorphique ou tectono-
métamorphique) est caractérisé par une composition carbonique à aquo-carbonique dans le
système H2O-CO2-CH4 ± NaCl ± H2S. Les conditions de pressions et de températures de
piégeage sont estimées par l'auteur de 0.8 kb à 5 kb et de 300°C à 400°C. Les quartz
argentifères de la deuxième phase (à gangue dolomitique ou épithermale argentifère) ont
piégé simultanément deux types de fluide. Le premier fluide est chaud (T°homogénéisation :
290°C à 180°C) et de salinité faible (T°fusion de la glace : 0°C à -15°C). Ce second fluide
présente des températures d'homogénéisation entre 190°C et 110°C, pour des températures de
fusion de la glace réparties entre -15°C et -30°C. Le second fluide, hypersalin, est formé dans
un système complexe H2O - NaCl - CaCL2 - MgCl2 - FeCl2 - LiCl. Baroudi et al. (1999)
proposent un modèle de mélange de fluides pour expliquer la formation du gisement d'Imiter.
Dans ce modèle, la précipitation de la paragenèse argentifère liée au stade à quartz est due à
un mélange de deux fluides, l'un associé au volcanisme (290-180°C et faible salinité) et l'autre
d'origine superficielle et hypersalin (190-110°C et saumures complexes).
65
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Une première étude isotopique (δ13C, δ18O et Pb/Pb) du gisement a été menée par
Pasava (1994) sur des fentes à carbonates et des galènes.
Les compositions isotopiques du carbone et de l'oxygène des fentes à carbonates ne
varient pas entre les échantillons stériles et les échantillons argentifères. Les valeurs de δ13C
et de δ18O sont comprises respectivement entre -10.5‰ et -8.1‰ et entre -13.4‰ et -0.1‰, ce
qui illustre, selon l'auteur, le rôle important de la matière organique dans l'origine des
dolomies et la présence de différents types de fluides.
Dans le diagramme de Zartman et Doe (1981) (Fig. 1.12), les résultats obtenus sur les
isotopes du plomb dans les galènes se répartissent en deux groupes compris dans les champs
des rides océaniques et des orogènes. Le premier groupe est interprété par l'auteur comme issu
d'un mélange entre une source mantellique océanique et les black shales du
Néoprotérozoïques moyen. Le second groupe, plus radiogénique, montre des similitudes avec
les données obtenues sur les gisements de Jbel Aouam (Maroc) et de El Hammam (Maroc),
66
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
Pasava (1994) suggère ainsi une origine en relation avec le volcanisme felsique
néoprotérozoïque terminal, ou avec des fluides associés à un événement plutonique ultérieur
(Varisque).
Figure 1.12 Distribution des analyses Pb/Pb (Pasava, 1994) dans le diagramme de Zartman et Doe (1981)
67
Partie II Chapitre 1 Généralités sur l'Anti-Atlas
68
Partie II Chapitre 2 Le volcanisme rhyolitique intrusif tardi néoprotérozoïque terminal
Chapitre 2
69
Partie II Chapitre 2 Le volcanisme rhyolitique intrusif tardi néoprotérozoïque terminal
Figure 2.1 Carte géologique et coupe N-S de la partie Ouest de la boutonnière d'Imiter (protrusion rhyolitique
de Takhatert)
70
Planche 2.1
71
Partie II Chapitre 2 Le volcanisme rhyolitique intrusif tardi néoprotérozoïque terminal
Figure 2.2 A Diagramme TAS de Le Maître (1989) des intrusions rhyolitiques et de la série volcanique
néoprotérozoïque terminal d'Imiter, données compilées Ouguir (1997) et Bajja (1987 ; 1998).
B Diagramme AFM des intrusions rhyolitiques et de la série volcanique néoprotérozoïque terminal d'Imiter,
données compilées Ouguir (1997) et Bajja (1987 ; 1998).
72
Partie II Chapitre 2 Le volcanisme rhyolitique intrusif tardi néoprotérozoïque terminal
C Diagramme SiO2 vs Log (CaO / (Na2O-K2O)) des intrusions rhyolitiques et de la série volcanique
néoprotérozoïque terminal d'Imiter, données compilées Ouguir (1997) et Bajja (1987 ; 1998)
terminal publiées antérieurement (Fig. 2.2 A). Les teneurs de K2O comprises entre 4.7 % et 8
% souligne le caractère fortement potassique de ce faciès. La teneur en alcalins (Na2O + K2O)
comprise entre 8.4% et 7.9% indique une affinité calco-alcaline à sub-alcaline. Les variations
du rapport Fe2O3/MgO entre 2,75 et 5,81 traduisent la présence de sulfures et/ou d'oxydes
dans la roche. Ces valeurs sont comparables à de nombreuses autres roches calco-alcalines
données en référence par Dostal et al. (1977). La mobilité attribuée aux alcalins dans de telles
paléo-volcanites limite de fait l'interprétation des résultats obtenus. Les teneurs en oxydes
(MgO, CaO, Fe2O3, MnO, TiO2 et P2O5) diminuent proportionnellement à l'augmentation de
la teneur en SiO2. Cette variation reflète la proportion et probablement la composition des
minéraux ferromagnésiens, ferrotitanés, des apatites et des plagioclases.
Le report des compositions de ces roches dans diverses représentations confirme leur
désignation pétrographique, ainsi que la nature calco-alcaline des faciès acides de la série
volcanique néoprotérozoïque terminale et tardive (Fig. 2.2 A-B-C). Les variations des teneurs
en éléments majeurs montrent que les différents groupes de faciès de la série
néoprotérozoïque terminal et les intrusions rhyolitiques tardives pourraient correspondre à une
lignée magmatique unique, dont l'évolution géochimique est interprétée en terme de
différenciation cristal/liquide.
73
Partie II Chapitre 2 Le volcanisme rhyolitique intrusif tardi néoprotérozoïque terminal
Figure 2.3 A Diagramme Zr/TiO2 vs SiO2 d'après Winchester et Floyd (1977) des intrusions rhyolitiques et de la
série volcanique néoprotérozoïque terminal d'Imiter
B Diagramme (Nb/Y )/SiO2 vs SiO2 d'après Winchester et Floyd (1977) des intrusions rhyolitiques et de la série
volcanique néoprotérozoïque terminal d'Imiter
C Diagramme (Nb/Zr)N vs Zr d'après Thieblement et Tergrey (1994) des intrusions rhyolitiques et de la série
volcanique néoprotérozoïque terminal d'Imiter
(3,38 <La/Sm)N< 5,74 ; 1,61 <Gd/Yb< 1,18) et une anomalie négative en europium
accentuée avec l'augmentation de la teneur en SiO2. Cette anomalie marque le fractionnement
précoce et important des feldspaths qui, d'après Steinberg et al. (1979) serait dû à une forte
variation de la fugacité en oxygène du magma.
74
Partie II Chapitre 2 Le volcanisme rhyolitique intrusif tardi néoprotérozoïque terminal
Figure 2.4 A Spectre étendu de variations des terres rares dans les intrusions rhyolitiques à Imiter. Normé au
manteau primitif, McDonough et Sun (1995)
B Spectre de variations des terres rares dans les intrusions rhyolitiques à Imiter. Normé aux chondrites, Evensen
et al. (1978)
75
Partie II Chapitre 2 Le volcanisme rhyolitique intrusif tardi néoprotérozoïque terminal
Figure 2.5 Comparaison des spectres moyens de variations des terres rares dans les intrusions rhyolitiques par
rapport à la série volcanique néoprotérozoïque terminal d'Imiter. Normé aux chondrite, Evensen et al.
(1978).
variations des éléments majeurs et traces mesurées dans les intrusions rhyolitiques indiquent
une affinité calco-alcaline. Comparées aux mesures réalisées sur les roches de la série
volcanique néoprotérozoïque terminal, les valeurs en terres rares légères déterminées sur les
intrusions rhyolitiques présentent un enrichissement relatif (Fig. 2.5) Cette observation tend à
distinguer la nature des intrusions rhyolitiques de celle de la série volcanique
néoprotérozoïque terminale, notamment les ignimbrites.
76
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
Chapitre 3
Au niveau du gisement d'Imiter I (Fig. 3.1 ; planche 3.1), les mesures ont été
effectuées dans la grande carrière et dans la galerie E-W du niveau 1370 m. La zone d'Imiter I
correspond à un nœud tectonique et se présente comme un complexe faillé d'allongement
77
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
Figure 3.1 A Carte géologique de surface du gisement d'Imiter I (d'après les documents SMI, modifiés) et régime
de déformation
B Coupe N-S du gisement d'Imiter I déformation (d'après les documents SMI, modifiés)
78
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
général E-W (Leistel et Qadrouci, 1991). Dans la carrière, il est difficile de distinguer une
faille principale, c'est un véritable réseau de failles anastomosées, dont l'expression
cartographique indique un mouvement normal décrochant sénestre. La faille nord matérialise
le contact entre les métasédiments du Néoprotérozoïque moyen et le bloc le plus effondré du
Néoprotérozoïque terminal. Ce contact est marqué par une caisse de failles mylonitisées de
puissance métrique. La direction de ce plan de faille évolue de N62°E à N110°E, ce qui lui
donne un aspect gauchi. Le bloc sud est haché par des failles de directions moyenne N60°E,
qui décalent en décrochement sénestre les structures satellites de la "faille d'Imiter" et
viennent se connecter sur cette dernière. L'ensemble a un aspect en tresses sigmoïdes.
45 plans de faille et leurs stries ont été mesurés. De nombreux plans de direction E-W
témoignent des jeux polyphasés successivement normaux, puis décrochants sénestres. Le
pendage principal des différentes structures est voisin de 90° vers le Nord (75-80°N).
Localement, ces structures sont associées à des failles antithétiques à pendage sud.
Une première famille correspondant au régime A, regroupe 27 failles de direction
moyenne N80° à fort pendage principalement Nord. Le mouvement déterminé sur les plans de
faille est un jeu normal avec un pitch de stries voisin de 90°. Cette famille de failles définit un
régime purement normal (σ1 = N38°E- 71°S ; σ2 = N67°E- 17°E; σ3 = N154°E - 9°S ; R =
(σ2−σ3)/(σ2−σ1) = 0,24).
Une seconde famille de failles correspondant au régime B, est composée de 18 plans
de faille. Elle est constituée de failles de direction N60° à N105°E avec un pendage fort (80° à
90°) vers le Nord ou le Sud. Les mouvements observés sur les plans de failles sont
décrochants sénestres (pitch de 15° W) sur les failles de direction E-W, et décrochants
normaux (pitch 66° W) sur les failles de direction N050°. Cette famille est composée d'une
partie des failles du régime A, qui ont rejoué en décrochement et de nouvelles failles N60°
qui accommodent le décrochement et jouent en faille décrochante normale (Fig. 3.2). Le
régime de contraintes déterminé pour cette population de failles (σ1 = N19°E - 54°S ; σ2 =
N29°E- 36°N; σ3 = N119°E- 5°N ; R = 0,79) est normal décrochant sénestre, avec une
direction d'extension WNW. Cette composante transtensive est bien illustrée par la valeur de
R > 0,5.
79
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
Figure 3.2 A Carte géologique de surface du gisement de la colline B3 (d'après les documents SMI, modifiés) et
régime de déformation
B Carte de la galerie d'exploitation niveau 1350 m.
C Coupe N-S du gisement de la colline B3 (d'après les documents SMI, modifiés)
80
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
Il s'agit du plus important gisement actuellement reconnu à Imiter. Le corps principal, d'une
puissance pluridécamétrique, présente des teneurs moyennes en Ag° de 2500 g.t-1. Un dyke
rhyolitique a été cartographié au niveau 1325 m du Puits III de la colline B3 dans une caisse
de faille normale (régime A) de direction N70° et de pendage 75°N reprise en jeu sénestre
(régime B) dans la série du Néoprotérozoïque moyen (planche 3.1). Deux générations de
stries ont été mesurées sur le plan de faille. La génération la plus ancienne présente un pitch
de 75°W, et la plus récente, un pitch de 15°W. La dernière génération de stries a également
été mesurées sur les épontes du dyke rhyolitique. Le dyke est tronqué par une faille relais
N40°E décrochante sénestre. L'ensemble de ces observations suggère que l'injection du dyke
est synchrone de l'extension N-S qui définit le régime A, repris ensuite dans le régime B.
Cette observation démontre que la mise en place des filons rhyolitiques est antérieure au
régime B. Ceci constitue un argument essentiel dans la détermination de la chronologie
relative des événements volcaniques et hydrothermaux à Imiter et sera repris lors de la
synthèse finale au chapitre 8.
Une première famille correspondant au régime A, est composée de 22 plans de failles
normales de direction moyenne N80° et de pendage proche de 90° Nord et Sud. Le régime de
contraintes déterminé est purement distensif (σ1 = N170°E - 84°N ; σ2 = N63°E - 2 E; σ3 =
N153°E - 6°S ; R = 0,27).
La seconde famille de failles correspondant au régime B, est composée de deux types
de failles. Le premier type correspond à un rejeu en décrochement sénestre des failles
normales du régime précédent. Le second type de failles correspond à des failles de direction
N40°E à pendage fort voisin de 90° vers le NW, qui se connectent aux failles E-W, ou les
recoupent en jeu décrochant sénestre. Le régime de déformation déterminé à partir de cette
famille de faille est décrochant sénestre transtensif (σ1 = N22°E - 17°S ; σ2 = N33°E- 73°N;
σ3 = N123°E - 4°E ; R = 0,81). Deux familles de failles supplémentaires de directions
respectives N100°E à N140°E et N30°E à N50°E ont été mesurées dans le Puits III. Les plans
de failles de direction moyenne N120°E présentent un jeu polyphasé décrochant sénestre et
normal. Les plans de failles de direction moyenne N40°E présentent des jeux inverses. Dans
les deux cas, le faible nombre de mesures ne permet pas de définir un régime de déformation
statistiquement représentatif. Les failles de direction N120°E sont postérieures au jeu distensif
du régime A. Ouguir (1994) a défini un régime "C" sur des failles de même direction.
L'existence de failles de direction N40°E inverses semble incongrue dans un ensemble de
régimes distensifs. En fait, elles s’intègrent parfaitement à la structure particulière d’une zone
81
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
Figure 3.3 A Carte géologique de surface du gisement d'Imiter IV (d'après les documents SMI, modifiés) et
régime de déformation
B Carte de la galerie d'exploitation niveau 1320 m
C Coupe N-S du gisement d'Imiter IV (d'après les documents SMI, modifiés)
82
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
Dans le gisement d'Imiter IV (planche 3.1), les observations et les mesures micro
structurales ont été réalisées en surface et au niveau 1320 m. Sur les quatre galeries de
reconnaissance ouvertes au Puits IV (Fig. 3.3), seul le niveau 1320 m recoupe la faille
d'Imiter. Cette galerie est une coupe N-S, longue de 240 m ; elle recoupe 20 m de tuf
néoprotérozoïque terminal et 40 m de conglomérat de base. Les zones minéralisées sont
réparties dans le Néoprotérozoïque moyen et le conglomérat de base. Les corps sont de
puissance métrique et sont associés à des zones de brèches. Les contacts entre les formations
sont soulignés par des caisses bréchiques de failles normales ouvertes, à pendage nord
(planche 3.1). Un dyke rhyolitique est injecté dans le contact conglomérat/néoprotérozoïque
moyen. Une quarantaine de couples faille/strie dont les caractéristiques de mouvement sont
bien déterminées, a été mesuré. Suivant la chronologie relative établie lors des observations
de terrains (striation double), les couples sont répartis en deux familles distinctes. La première
famille correspondant au régime A, est composée de 17 failles de direction moyenne N100°E
à pendage fort principalement Nord. Le mouvement déterminé sur les plans de failles est un
jeu normal avec un pitch voisin de 90°. Cette famille de failles définit un régime purement
normal (σ1 = N150°E- 60°S ; σ2 = N86°E- 18°E ; σ3 = N177°E- 13°N ; R = 0,49). La
seconde famille de failles correspondant au régime B, présente des directions de plans de
failles moins homogène. Elle est constituée de failles de direction N50°E à N110°E avec un
pendage variable de 90° à 50°N. Les jeux observés sur les plans de failles sont décrochants
sénestres (pitch de 15°W) pour des failles de direction E-W, et décrochants normaux (pitch
66° W) pour les failles de direction N50°E. Cette famille est composée d'une partie des failles
antithétiques du régime A qui ont rejoué en décrochement et de nouvelles failles N50°E, qui
accommodent le décrochement et jouent en faille normale (Fig. 3.3). Le régime de contraintes
déterminé sur cette population de failles est normal décrochant sénestre (σ1 = N36°E- 69°S ;
σ2 = N30°E- 19°N ; σ3 = N121°E- 5°E ; R = 0,48).
83
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
Figure 3.4 A Carte géologique de surface du gisement d'Imiter V (d'après les documents SMI, modifiés) et
régime de déformation
B Coupe N-S du gisement d'Imiter V (d'après les documents SMI)
84
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
Le Puits V (Igoudrane ; planche 3.1) n'est pas situé sur la faille d'Imiter mais sur une
structure satellite au Sud de celle-ci. En effet, au niveau du Puits IV la structure principale
amorce une virgation de direction N60°, pour recouper finalement la granodiorite d'Igoudrane
au nord du Puits V. En 2000, le Puits V est réouvert et donc de nouveau accessible, après
avoir été longtemps ennoyé. L'état d'oxydation des galeries a limité, de fait, les observations
réalisables. Le seul niveau accessible était le niveau 1270 m (-240 m). Il s'agit d'une galerie de
reconnaissance, qui recoupe le contact entre les métasédiments du Néoprotérozoïque moyen
et la diorite d'Igoudrane (Fig. 3.4). Le Puits V est caractérisé par des couloirs de failles
bréchiques de puissance métrique à pluri-métrique dans les métasédiments. Ces couloirs
contrôlent la mise en place des corps minéralisés et des dykes andésitiques (Fig. 3.4 ; planche
3.1). Trente couples faille-strie ont été mesurés. Les directions de plans de failles et de
pendages, toutes deux trés homogènes, sont respectivement E-W et 80°N. La majorité des
failles montre un jeu biphasé, normal puis décrochant. Le jeu normal (régime A) est mesuré
sur 18 plans de failles, il permet de différencier un régime de déformation purement normal
(σ1 = N169°E - 69°S ; σ2 = N75°E - 10°E; σ3 = N161°E - 19°N ; R = 0,17) résultant d'une
distension NNW-SSE. Le jeu décrochant correspondant (régime B ) est mesuré sur 12 plans
de failles. Ce sont des plans de failles polyphasés qui appartiennent aux deux régimes. Le
régime de contrainte déterminé sur cette population est décrochant sénestre transpressif (σ1 =
N42°E - 17°S ; σ2 = N173°E - 65°S ; σ3 = N127°E - 18°E; R = 0,17). Il se manifeste à
l'affleurement par des couloirs de brèches métriques, localisés au niveau des failles E-W. Des
failles de direction N50°E à N00°E ont été observées dans les travaux miniers et en surface,
mais il a été impossible d'identifier et de mesurer une strie ou tout autre tectoglyphe témoin du
jeu de ces structures.
3-1-3. Conclusion
L’étude des plans de failles striés dans les gisements, le long du faisceau de failles
d'Imiter, a permis de mettre en évidence deux régimes de déformation majeurs. Le premier
(régime A) est un régime distensif de direction N-S à NNW-SSE. Il a été défini sur
l'ensemble de la zone de failles d'Imiter et de la zone minéralisée. Le second régime (régime
B) est de type décrochant sénestre avec des axes σ1 et σ3 de directions générales respectives
NNE-SSW et SE. Nous confirmons ici les résultats d'Ouguir et al. (1994). L'évolution du
régime A vers le régime B, se traduit par une rotation horaire des axes de contrainte. Cette
évolution est caractérisée par des variations locales de l'intensité de la composante
85
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
Figure 3.5 Synthèse des régimes de déformations A et B définis à l'échelle du faisceau de la faille d'Imiter.
86
Planche 3.1
III et qui ont été associées au jeu décrochant sénestre du faisceau de faille d'Imiter. Le régime
C correspond selon l'auteur à un régime distensif tardif. Le nombre limité de mesures
réalisées sur ce type de failles ne nous permet pas de déterminer un régime de déformation
statistiquement réaliste. Toutefois, nos observations permettent de confirmer le caractère
tardif de ces failles par rapport aux régimes A et B. L'aspect basculé des corps minéralisés en
direction du Sud Ouest témoigne de l'influence de ce régime sur la minéralisation. Des failles
décrochantes dextres (pitch 15°E) ont été identifiées de façon éparse sur l'ensemble du
faisceau de failles. Ces stries correspondent systématiquement au jeu le plus tardif. Le faible
nombre de mesures réalisées sur ce type de plan de faille ne permet pas d'identifier le régime
de déformation auquel elles appartiennent.
87
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
normal plus marqué du régime B dans la partie Ouest de la boutonnière. En outre, le régime
décrochant sénestre conduit à une forme sygmoïde des corps minéralisés.
Les deux régimes de déformations A et B que nous venons d'identifier et de d'écrire
sont responsables de la mise en place de la minéralisation argentifère sous la forme de veines
à gangue de quartz et/ou de dolomite. Ces veines correspondent
1- au plan principal d'ouverture N70-90°E du régime A.
2- aux plans de cisaillement (N70-90°E) du régime A réactivés (brèches) et à des
fentes de tension N40-50°E décrochantes sénestres du régime B.
Ces veines minéralisées constituent les structures principales à Imiter. Elles sont
accompagnées par des structures minéralisées de moindre importance comme les veines
horizontales observées, par exemple, au Puits III. Ces veines horizontales correspondent aux
zones de push-up associées aux failles décrochantes du régime B (Ouguir et al. 1994).
Au total, les failles principales et les satellites conduisent, à toutes les échelles, à une
intense fracturation des roches encaissantes et par conséquent à la dissémination des fluides
minéralisateurs dans l'enveloppe de structures principales.
Cette analyse des fentes de tensions nous a permis de mettre en évidence une première
relation spatiale et temporelle entre la minéralisation et le réseau de failles d'Imiter. Dans un
premier temps la minéralisation est supportée par une gangue à quartz, puis à dolomite,
suivant une extension N-S évoluant de façon horaire en N130°E (passage du régime A au
régime B). Ces régimes de déformation conditionnent la mise en place de la minéralisation.
88
Partie II Chapitre 3 Etude structurale de la bordure Nord de la boutonnière d'Imiter
89
Partie II Chapitre 4 Etude paragénétique
CHAPITRE 4
Etude paragénétique
4-1. Introduction
Les échantillons analysés ont été prélevés dans les amas minéralisés des différents
gisements de la boutonnière d'Imiter, dans le socle métasédimentaire néoprotérozoïque moyen
et la couverture volcano-sédimentaire néoprotérozoïque terminale. Nous avons complété nos
observations en microscopie optique, par des analyses fines au microscope électronique à
balayage et par des analyses ponctuelles à la microsonde électronique, qui ont permis de
préciser la nature exacte des différentes phases. Cette étude a été menée en collaboration avec
Ph. Marion (ENSG-LEM). Une attention plus particulière est portée sur les minéralisations
hydrothermales.
Cette étude vise à la caractérisation et à l'établissement d'un schéma global des
diverses associations et successions de dépôts, ce qui permettra de discuter les conditions de
formation du minerai et servira de base à la réalisation des études isotopiques.
Au cours de cette étude nous avons distingué trois épisodes successifs dans le district d'Imiter.
91
Partie II Chapitre 4 Etude paragénétique
dehors du halo de minéralisation, ce qui était le cas également pour les auteurs précédents
(planche 4.1)
Par exemple, l'échantillon IM99-03 a été prélevé dans une veine de quartz-dolomite à pyrite,
chalcopyrite, galène, freibergite, muscovite (Ms1-2) et chlorite (chl1) à remplissage centripète, encaissée dans
les tufs rhyolitiques lités de la base de la série néoprotérozoïque terminale (Puits III, 1450 m). L'étude
microscopique permet de déterminer la chronologie relative suivante, quartz dendritiques à extinction onduleuse,
pyrite, galène, chalcopyrite et muscovite (Ms1), chlorites (chl1), puis freibergite, muscovite "en flamme" (Ms2)
et dolomite. L'association à freibergite, muscovite (Ms2) et dolomite correspond à l'épisode épithermal
argentifère ultérieur (planche 4.3)
L'échantillon IM99-37 a été prélevé dans une fente à quartz (Qtz1), dolomite, chlorite
(chl1) et pyrite dans des black shales (Puits IV, niveau 1320 m). La teneur en Ag de la zone de prélèvement est
supérieure à 1000 g/t. Les observations pétrographiques mettent en évidence deux générations de muscovite
(Ms1-2), mais pas d'argent métal. Les muscovites (Ms1) sont en position inter-granulaire entre les quartz (Qtz1),
en inclusion dans les sulfures. Les muscovites (Ms2) tapissent le long des épontes des fractures à dolomite.
(planche 4.3).
L'échantillon IM99-58 a été prélevé dans une fente à quartz (Qtz1), chlorite (chl1), pyrite
et sphalérite, dans les black shales, (Puits III, niveau 1325 m). La zone de prélèvement présente des teneurs en
Ag de 1500 g/t. La minéralisation argentifère n'est pas visible dans la lame. Les muscovites (Ms1) sont de très
petites tailles. Elles se présentent soit dans les joints inter-granulaires, soit en inclusions dans les quartz (Qtz1).
(planche 4.3).
92
Partie II Chapitre 4 Etude paragénétique
L'échantillon IM00-53 a été prélevé dans une cornéenne de contact entre la lame de
granodiorite de la colline B3 et les black shales du Néoprotérozoïque moyen au Puits III, niveau 1335 m.
L'échantillon se compose de deux zones pétrographiques différentes. La première est composée uniquement de
quartz (Qtz1) et de muscovite (Ms1). Le quartz présente des textures graphiques typiques de filons
pegmatitiques. La seconde partie correspond à la zone black shales. Elle est caractérisée par une "imprégnation"
de quartz et de muscovite dans la mésostase des black shales.
Les observations en microscopie à balayage (MEB), avant irradiations des muscovites Ms1 (IM00-53 ;
IM99-37) analysées, montrent des zones sombres le long des fractures, des bordures et des espaces interfoliaires.
Nous interprétons ces zones comme des zones de déstabilisation ou de recristallisation.
93
Partie II Chapitre 4 Etude paragénétique
Par exemple, l'échantillon IM00-18 a été prélevé dans le Puits III au niveau 1454 m. Il s'agit d'une veine à
quartz dans les black shales du Néoprotérozoïque moyen avec des injections locales de dolomites. La paragénèse
est composée de chalcopyrite, galène, arsénopryrite, sulfoarséniure de Co-Ni, pyrite et argentite. Les
sulfoarseniures de Co Ni sont associés aux injections de dolomites (planche 4.1)
Par exemple, l'échantillon IM00-22 a été prélevé dans le Puits III au niveau 1454 m. Il s'agit d'une veine
de dolomite dans le micro conglomérat de base de la série néoprotérozoïque terminale. Un liseré fin de quartz
dendritique non minéralisé marque le contact entre le micro-conglomérat et la veine de dolomite. La
minéralisation est composée de galène, d'arsenopyrite, de sphalérite de chalcopyrite et d'argentite. La galène
présente des inclusions d'argentite. La galène et la sphalérite sont tardives par rapport à l'arsenopyrite. La plage
d'argentite en baguette cristallise en dernier (planche 4.1)
94
Planche 4.1 Photographies au MEB en électrons rétrodiffusés
Photographie 2 IM99-35 Puits III, niveau 1374m. Veine à quartz argentifère (Qtz2) et
dolomite dans les tufs lités du Néoprotérozoïque terminal. La dolomite est tardive par
rapport au quartz. La minéralisation est essentiellement composée d'amalgame
d'argent associé au quartz en mosaïque.
Photographie 4 IM99-54 Puits III, niveau 1335 m. Cet échantillon est typomorphe de la
phase à dolomite. La minéralisation présente des textures en cupules composées de
sulfoarséniures de Co-Ni à coeur de pyrite et amalgame d'argent. L'amalgame d'argent
remplace la pyrite dans les cupules. Les relations entre l'amalgame d'argent et les
sulfoarséniures de Co-Ni ne sont pas claires.
Photographie 5-6 IM00-22 Puits III, niveau 1320 m. Veine à dolomite, argentite, galène,
arsenopyrite et sphalérite. L 'argentite cristallise en fibres avec des intercalations de
galène. La galène forme des plages polycristallines avec la sphalérite et l'arsénopyrite
qu'elle remplace.
Planche 4.3
Photographie 2 -IM99-03, Puits III niveau 1450m. : Veine à quartz (Qtz1), à remplissage
centripète de pyrite, chalcopyrite, chlorite (Chl1), sur lesquels croissent des
muscovites (Ms2) puis de la dolomite.
Photographie 4 -IM99-37, Puits IV niveau 1325m. Veine à quartz (Qtz1) avec veine de
chlorite (Chl1), sur lesquelles se développent des muscovites (Ms2) dans des zones de
réouvertures à dolomite.
Photographie 5 et 6 -IM99-58, Puits III, niveau 1320m. Veine à quartz (Qtz1) dans des
métasédiments néoprotérozoïque moyen, à plages polyminérales de pyrite et sphalérite
embalées dans les chlorites (Chl1).
NB. Les cercles noirs (photographies 3 et 4) correspondent aux zones repérées avant l'analyse
microsonde électronique et MEB. Les ellipses noirs (photographie 5) correspondent au
spot d'analyse de la microsonde ionique (IMS 3F ou IMS 1270)
Planche 4.4 : Photographies au MEB en électrons rétrodiffusés
95
Partie II Chapitre 4 Etude paragénétique
96
Partie II Chapitre 4 Etude paragénétique
La loellingite : ((Fe, Co, Ni) AsS) n'est repérée que dans la paragénèse carbonatée,
elle se présente en petites plages autour des galènes, ou associée à de l'argentite.
La cobaltite : ((Fe, Co, Ni) AsS) n'a été mise en évidence qu'en microscopie
électronique à balayage. Elle se présente sous divers aspects, en dépôts eutectiques avec de
l'amalgame ou de l'argentite, ou en bordure de structures mirmékitiques en "rosette" avec au
cœur de la pyrite et/ou de l'argentite. La composition de la cobaltite présente de faibles
teneurs de Ni (4 à 15%). La composition de la cobaltite d'Imiter est comparable à celle définie
pour le gisement de Cobalt au Canada, par Petruk et al. (1971).
La gersdorphite : ((Ni, Co, Fe) As S) n'a été identifiée qu'en microscopie
électronique. Elle se présente en très petites plages en association avec l'argentite et
l'amalgame.
L'amalgame d'argent (AgHg) : c'est la phase la plus abondante. Son expression la
plus spectaculaire est sous forme de plaquettes dans les veines carbonatées et/ou squelettique
dans les veines de quartz (quartz gris). A son contact la pyrite, la blende et l’arsénopyrite sont
toujours corrodées. Il se présente en liseré de déstabilisation de certains sulfures argentifères,
en inclusion dans l’arsénopyrite, l’argentite, et les sulfosels d’argent. L’amalgame présente
une fluctuation importante des teneurs en Ag (77-95%). La teneur en Hg est toujours
inférieure à celle de la Kongsbergite (40% Hg). Cette composition permet de considérer
97
Partie II Chapitre 4 Etude paragénétique
l’amalgame d’Imiter comme de l’argent mercureux α-phase, ou le mercure entre en tant que
solution solide dans la structure d’argent. Les zonations observées en microscopie
électronique illustrent la variation du rapport Ag/Hg. Les analyses à la microsonde révèlent la
présence épisodique d’éléments traces (Bi, Te) dans l'ensemble des occurrences et d'impuretés
de Fe (<1%), Co (<0.7%) et Ni (<0.6%) dans la phase carbonatée. L’amalgame présente un
accroissement progressif de la teneur en Hg entre la phase à quartz et la phase carbonatée
(Fig. 4.4). L'augmentation de la teneur en mercure est corrélée avec l'apparition du cobalt et
du nickel dans les solutions au début de la phase carbonatée. De l’argent pur (1%Hg) a été
observé principalement dans les dépôts à gangue silicatée et en lisière des sulfo-arséniures.
Dans le second cas, il pourrait s’agir de déstabilisation des sulfosels. Des fils d'argent
probablement d'origine supergène, ont également été observés dans les géodes de quartz
tardif.
L'argentite (Ag2S) : elle apparaît sous sa forme stable l'acanthite, mais nous
utiliserons le nom courant d'argentite dans ce qui suit. Les argentites observées en
microscopie et au MEB présentent toutes un système cristallin de type dodécaèdrique. En
lame mince, elle apparaît sous forme interstitielle entre les grains de quartz ou de plages
irrégulières polymétalliques. Associée à l'amalgame, l’argentite se présente en grandes plages
polymétalliques, remplaçant les sulfures précoces (pyrite, blende, arsénopyrite) et le
polybasite. Elle est pure dans la phase à quartz, et présente des traces de Cu et d'As dans la
phase à dolomite. Dans cette dernière phase, elle se dépose également au sein de nids de Co-
Ni sulfoarséniures (ou cupule) en association avec Ag-Hg présente en association avec la
pyrite, ou les (Co, Ni) sulfoarséniures. Ces textures sont interprétées comme le résultat de
dépôts simultanés, eutectiques, traduisant des modifications de chimisme des solutions.
98
Partie II Chapitre 4 Etude paragénétique
99
Partie II Chapitre 4 Etude paragénétique
reconnus sont représentés par la malachite, l'azurite, l'érythrite, les oxydes de Mn, la cérusite,
la limonite et la kaolinite ainsi que par des fils d'argent dans des géodes de quartz laiteux
tardifs. Ce dernier dépôt d'argent est limité en volume et n'a rien de commun avec le dépôt des
minéralisations associées à l'épisode épithermal argentifère. D'une manière générale
l'existence d'un épisode d'enrichissement supergène (Vargas, 1983 ; Guillou et al., 1988 ;
Popov 1995) ne resiste pas à l'examen des associations paragénétiques que nous venons de
décrire. En effet, le synchronisme du dépôt de l'Argent natif et des autres sulfures et sulfosels
(notamment à Co-Ni) correspond à un dépôt primaire de l'argent dans des structures
minéralisées (cf. photos). En outre, l'augmentation progressive de la teneur en mercure des
minéralisations es veines de quartz aux veines de dolomites implique nécessairement le dépôt
de la phase argentifère lors de le formation des structures minéralisées et non pas lors d'une
phase ultérieure supergène de remobilisation (cémentation).
4-3. Thermométrie
Figure 4.5 A Diagramme de caractérisation cristallochimique des phengites selon Velde, (1965)
B Diagramme de détermination de température des phengites d'après Monier et Robert, (1986)
100
Partie II Chapitre 4 Etude paragénétique
0.71 et 0.97 a./u.f, un contenu en AlIV compris entre 0.60 et 0.91, et un contenu en Mg et Fe
respectivement de 0.03 à 0.26 et de 0.07 à 0.21. Dans le diagramme de Velde (1965), ces
micas blancs sont tous répartis sur la droite reliant le pôle phengite au pôle muscovite avec
une bonne distinction entre Ms1 proche du pôle muscovite et Ms2 plus proche du pôle
phengite (Fig. 4.5A). Les micas contiennent peu de Na (0,01<Na<0,09) et pas de Ca. La
contribution du pôle paragonitique est donc très faible.
Dans le diagramme de Monier et Robert (1986) nos résultats sont répartis en deux
groupes (Fig. 4.5B). Le premier groupe est constitué par des muscovites (Ms1) des
cornéennes de contact de la lame de la granodiorite de la colline B3 (IM00-53) et des veines
quartz (Qtz1) (IM99-03-37-58). On observe une gamme de variation de température entre
650°C et 400°C pour ce premier groupe. L'échantillon (IM00-53), prélevé au contact de la
lame granodioritique de la colline B3, présente les températures les plus élevées, de l'ordre de
~600°C. Il semble exister une correlation entre la température de cristallisation des
muscovites Ms1 et leur distance par rapport à la lame de granodiorite de la colline B3. Le
second groupe est composé de muscovites Ms2 (IM99-03-37). Les températures déterminées
sur ces muscovites (Ms2) sont de ~350°C à 275°C respectivement pour les échantillons
IM99-37 et IM99-03. Les analyses effectuées sur les muscovites de l'échantillon IM99-03
sont en dessous du domaine de validité défini par les auteurs (<300°C).
101
Partie II Chapitre 4 Etude paragénétique
(T°moy = 350 ± 15°C) et non au niveau 1320 m (T°moy = 320 ± 15°C). Le moteur thermique
responsable de la circulation des fluides ne semble pas se situer en profondeur.
Figure 4.6 A Diagramme de caractérisation cristallochimique des chlorites selon Bayliss, (1985)
B Courbe de densité (histogramme de densité) des températures calculées selon la formule de
Cathelineau, (1988)
Sur les inclusions fluides, Baroudi et al. (1999) déterminent une gamme de
température comprise entre 290°C et 180°C donc inférieure à la gamme de température
déterminée sur les chlorites et les muscovites de l'épisode hydrothermal à métaux de base.
La limite basse de température à 180°C déterminée par Baroudi et al. (1999) est
confirmée par la forme dodécaèdrique de l'acanthite (observation microscopique et MEB)
dont la température minimale de cristallisation est estimée à 177°C (Kracek, 1946).
102
Partie II Chapitre 4 Etude paragénétique
Figure 4.7 Courbes d'équilibre de différentes réactions entre les sulfures dans un espace as2-T°, d'après Barton
et Skinner, (1967)
103
Partie II Chapitre 4 Etude paragénétique
Les erreurs sur la température et sur la valeur du logarithme sont respectivement de l'ordre de
±10°C et de ±0.2. En utilisant, la gamme de températures déterminées précédement, la
variation de fO2 des gaz à l'équilibre avec le système minéralisateur peut être estimée de 10-30
atm (à 350°C) à 10-39 atm (à 230°C).
104
Partie II Chapitre 4 Etude paragénétique
Figure 4.8 Diagramme de stabilité de l'argent à 250°C en fonction du pH et de la fugacité en oxygène selon
Gammons et Barnes (1989).
4-7. Conclusion
105
Partie II Chapitre 4 Etude paragénétique
106
Partie II Chapitre 5 Analyses isotopiques du soufre dans les intrusions rhyolitiques et la minéralisation
Chapitre 5
et la minéralisation
107
Partie II Chapitre 5 Analyses isotopiques du soufre dans les intrusions rhyolitiques et la minéralisation
1980). Sakai et al. (1984) et Mandeville et al. (1998) ont étudié les variations des isotopes du
soufre, respectivement dans les basaltes du Kilauea et dans la série éruptive du Kratatau. Ils
mettent en évidence dans le cas du volcanisme subaérien l'importance du dégazage de SO2 sur
le fractionnement isotopique. Le dégazage rapide de SO2 dans des magmas ascendants induit
un déséquilibre isotopique entre les sulfates et les sulfures dans les liquides résiduels et les
Tableau 5.1 Tableau synthétique des résultats de δ34S. Abréviations : Py, pyrite ; Gal, galène ; Sph, sphalérite ;
Cpy, chalcopyrite ; Pyrar, pyrargyrite ; Apy, arsenopyrite ; Cin, cinabre ; IMS, sonde ionique CAMECA
IMS3F ; SG, spectromètre de masse à source gazeuse VG6002; IRMS, "Isotope Ratio Mass Spectrometry".
108
Partie II Chapitre 5 Analyses isotopiques du soufre dans les intrusions rhyolitiques et la minéralisation
vapeurs. Ce déséquilibre induit à son tour un fractionnement important des sulfures dans les
liquides résiduels. Lors de l'ascension brutale des magmas dans des conditions de fO2 faibles
le dégazage de SO2 enrichit le liquide résiduel en isotope léger. Les résultats présentés par les
auteurs suggèrent qu'il n'y a pas de ré-équilibration isotopique dans les magmas qui ont
rapidement dégazés juste avant leur solidification. La protrusion de Takhatert montre des
évidences de dégazage (flammes importantes), le dyke de Tachkakacht non. Les valeurs de
δ34SCDT obtenues sur ces deux échantillons sont respectivement de ca -7 ‰ et ca -2 ‰. Dans
le cadre du modèle proposé par Sakai et al. (1984), ces valeurs sont en accord avec les
observations géologiques. La variation des valeurs de δ34SCDT mesurées est également
comparable aux variations déterminées par Sakai et al. (1984) de l'ordre de 7‰.
109
Partie II Chapitre 5 Analyses isotopiques du soufre dans les intrusions rhyolitiques et la minéralisation
Figure 5.1 Distribution des valeurs de δ34SCDT par phase de minéralisation et par type de sulfures
110
Partie II Chapitre 5 Analyses isotopiques du soufre dans les intrusions rhyolitiques et la minéralisation
111
Partie II Chapitre 5 Analyses isotopiques du soufre dans les intrusions rhyolitiques et la minéralisation
(négatives ou trop élevées) avec les températures estimées pour cette phase épithermale (350-
180°C). Ceci implique que l'équilibre isotopique n'est jamais atteint dans le système. Il
convient donc de trouver une autre explication à cette dérive du δ34SCDT à Imiter. Une partie
des variations isotopiques observées dans le dyke rhyolitique de Tachkakacht et la protrusion
de Takhatert (-2 à -7 ‰) peut s'expliquer par un dégazage de SO2 des fluides ascendants
entraînant un fractionnement préférentiel de l'isotope 34S et donc un enrichissement relatif en
32
S dans les sulfures (Sakai et al. 1984). L'ordre de grandeur de cette dérive mise en évidence
à Imiter est comparable à celle observée dans la série éruptive du Krakatoa (Mandeville et al.
1998). Aux variations près dont nous venons de discuter, les valeurs de δ34SCDT déterminées
sur les pyrites des intrusions rhyolitiques, représenteraient les valeurs des sulfures à l'équilibre
avec le réservoir magmatique. Ce sont les valeurs les plus lourdes mesurées à Imiter et leur
limite supérieure est -2 ‰.
Toutes les valeurs mesurées dans les deux épisodes minéralisés (épisode hydrothermal
à métaux de base et épisode épithermal argentifère Qtz2 et dolomite) varient entre cette valeur
de -2 ‰ et les valeurs les plus légères déterminées dans les pyrites de la phase tectono-
métamorphique qui s'étalent entre -30 ‰ et -38 ‰. Il est donc tout à fait réaliste d'interpréter
la dérive des valeurs dans les phases minéralisées comme représentative d'un mélange de
deux réservoirs équilibrés, l'un avec le soufre magmatique et l'autre avec le soufre des black
shales du Néoprotérozoïque moyen. A Imiter, l'étude des isotopes du soufre met donc
clairement en évidence deux sources possibles pour cet élément : l'une magmatique, liée soit
au volcanisme rhyolitique (épisode épithermal argentifère) soit aux granodiorites (épisode
hydrothermal à métaux de base), l'autre métasédimentaire liée aux black shales du
Néoprotérozoïque moyen et surtout exprimée dans la phase hydrothermale argentifère.
En conclusion, cette étude des isotopes du soufre ne permet pas de trancher clairement
sur l'origine de la minéralisation d'Imiter et laisse la porte ouverte aux deux sources possibles.
Il nous est donc apparu indispensable d'utiliser d'autres traceurs isotopiques pour essayer de
résoudre cette question. Dans ce but, nous avons choisi le couple Re-Os.
112
Partie II Chapitre 6 Analyse isotopique Re/Os
Chapitre 6
113
Partie II Chapitre 6 Analyse isotopique Re/Os
et IM00-37(2) contiennent des quantités d'Os inférieures à celle du blanc et donc ne seront pas
discutés. Les blancs de Re varient de façon aléatoire entre 2 et 22 pg. Les incertitudes données
pour [Re], et les rapports isotopiques 187Re/188Os et 187Os/188Osinitiale prennent en compte cette
variation des valeurs de blanc de Re.
114
Planche 6.1
Figure 6.1 Diagramme d'évolution du rapport 187Os/188Os en fonction du temps (Les échantillons prélevés dans
la phase argentifère sont représentés en gras)
être expliquée soit par un problème d'ordre analytique lié à des rendements d'attaque et de
187
dissolution variables et qui faussent les rapports Re/188Os, soit par une contribution
115
Partie II Chapitre 6 Analyse isotopique Re/Os
importante du blanc pour les échantillons pauvres en Os, soit par le fonctionnement du
système isotopique en système ouvert (Foster et al., 1996). Cette dernière hypothèse est
favorisée dans le cas des échantillons de black shales. Dans les black shales l'Os et le Re ne
sont pas concentrés dans les sulfures mais pourraient être associés à la matière organique et
donc plus solubles dans les fluides météoriques et/ou hydrothermaux (Pierson - Wickmann et
al., 1999). Ces auteurs et Jaffe et al. (2000) proposent une perte de Re et d'Os par
remobilisation et altération (Xiong et Wood, 1999). Peuker-Ehrenbrink et Hannigan, (2000)
montrent que les black shales altérés perdent de 45% à 90% de leur teneur en Os initial. Dans
notre contexte géologique cette hypothèse semble la plus réaliste. Dans le cas d'un système
isotopique clos depuis la cristallisation (par exemple : IM99-74) le manque de reproductibilité
des analyses dupliquées est probablement expliqué par l'addition d'une quantité inappropriée
du traceur à l'échantillon lors de la première analyse du fait de l'absence de données préalables
sur les teneurs en Os et Re de l'échantillon. La nature hétérogène des phases analysées et la
répartition en "pépite" de l'Os et du Re dans les phases sulfurées contribuent également à
compliquer l'équilibre échantillon-traceur.
Pour tous les échantillons de l'épisode épithermal argentifère, quelle que soit la phase
analysée (sulfure ou amalgame d'argent, IM99-51B ; IM00-37-49-50-52), les rapports
187
Os/188Os sont localisés sur la droite d'évolution du manteau, ce qui indiquent sans
ambiguïté la dominance d'une source mantellique des métaux (Fig. 6.1). Ceci est un résultat
majeur obtenu grâce à l'utilisation du traceur Re/Os. Ces analyses complètent et précisent les
résultats obtenus par les analyses isotopiques du soufre.
Il est possible d'estimer la contribution éventuelle du réservoir isotopique
187
métasedimentaire du Néoprotérozoïque moyen en supposant que la valeur Os/188Os de la
minéralisation argentifère est le résultat de la contribution de deux pôles purs, représentés par
le manteau d'un côté (187Os/188Os = 0.12) et par les black shales de l'autre (187Os/188Os 550Ma =
inconnu).
La signature isotopique 187Os/188Os des mers précambriennes est inconnue. Aussi nous
supposons que la valeur 187Os/188Os initiale (à ca. 800 Ma) minimale est identique à celle du
manteau (187Os/188Os = 0.12). La contribution radiogénique en 187Os, entre 800 Ma et 550 Ma,
187
est déterminée en supposant que la valeur Os/188Os actuelle correspond à 800 Ma de
187
contribution radiogénique dans un système non perturbé. La valeur du rapport Os/188Os à
116
Partie II Chapitre 6 Analyse isotopique Re/Os
550 Ma est estimée à ~0.822, ce qui représente une valeur minimale. En fonction de ces
hypothèses nous avons déterminé pour chaque échantillon caractérise de l'épisode de la
minéralisation argentifère la contribution maximale d'Os dérivée du réservoir black shale à
partir des rapports isotopiques 187Os/188Os initiaux. Les résultats sont présentés dans le tableau
6.2. La contribution maximale estimée est de l'ordre de 10 %, elle est probablement trés sur-
évaluée.
Tableau 6.2 Contribution de l'Os dérivé des black shales dans la minéralisation argentifère en fonction des
187
rapports isotopiques initiaux calculés à 550 Ma. Le rapport Os/188Os initial de l'écahantillon IM00-37, est
irréaliste, aussi nous avons utilisé le rapport mesuré pour obtenir une estimation maximale pour cet échantillon
117
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
Chapitre 7
et de la minéralisation
7-1. Problématique
Trois échantillons ont été analysés; ils proviennent des intrusions rhyolitiques de
Tachkakacht et de Takhatert et de la granodiorite de Taouzzakt. La position stratigraphique
relative et les caractéristiques pétrographiques des ces objets sont décrites dans le chapitre 3
de la partie II de ce manuscrit.
119
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
La granodiorite de Taouzzakt
207
Seize grains de zircons ont été analysés. Les âges Pb/235U et 206
Pb/238U sont
discordants pour la plupart des analyses réalisées. La distribution des âges 206Pb/207Pb dans le
diagramme concordia révèle trois groupes différents (Fig. 7.1 ; tableau 7.1). Les deux
premiers groupes sont constitués d'âges concordants répartis respectivement entre 569 ± 4 Ma
et 579 ± 4 Ma (âge moyen pondéré : 573 ± 4 Ma ; MSWD : 2,6 ; n = 4), et entre 653 ± 3 Ma
et 677 ± 6 Ma. Le troisième groupe est constitué de sept âges discordants répartis le long
Tableau 7.1 Résultats analytiques de la datation U/Pb par sonde ionique des zircons de la granodiorite de
Taouzzakt. Les rapports isotopiques du plomb sont corrigés du plomb commun. Les rapports U/Pb sont
déterminés en utilisant la corrélation linéaire déterminée sur les standards.
120
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
206
Figure 7.1 Diagramme Pb/238U vs 207
Pb/235U des analyses à la sonde ionique CAMECA IMS 1270 sur
monograin de zircon de la granodiorite de Taouzzakt.
d'une discordia dont l'intercept supérieur à 572 ± 5 Ma (MSWD = 0,6 ; n = 10) est localisé au
niveau du premier groupe à (573 ± 4 Ma). Cet âge est considéré comme l'âge de mise en place
de la granodiorite de Taouzzakt. L'intercept inférieur de la discordia tend vers l'origine de la
courbe concordia. La droite discordia représente une perte en plomb continue. Les deux âges
concordants à ca. 650 Ma sont réalisés sur des zircons monophasés. Ces grains de zircons sont
des xénocristaux qui témoignent de l'existence d'un événement magmatique contemporain de
la déformation B1 (Leblanc et Lancelot, 1980). Ces sont les âges les plus anciens mis en
évidence dans l'Anti-Atlas oriental.
121
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
Ma et 650 Ma, Leblanc et Lancelot, 1980). Les cristaux présentant des cœurs hérités n'ont pas
été analysés.
Tableau 7.2 Résultats analytiques de la datation U/Pb par sonde ionique des zircons de l'intrusion rhyolitique de
Takhatert. Les rapports isotopiques du plomb sont corrigés du plomb commun. Les rapports U/Pb sont
déterminés en utilisant la corrélation linéaire déterminée sur les standards.
206
Figure. 7.2 Diagramme Pb/238U vs 207
Pb/235U des analyses à la sonde ionique CAMECA IMS 1270 sur
monograin de zircon de la protrusion rhyolitique de Takhatert
122
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
Tableau 7.3 Résultats analytiques de la datation U/Pb par sonde ionique des zircons du dyke rhyolitique de
Tachkakacht. Les rapports isotopiques du plomb sont corrigés du plomb commun. Les rapports U/Pb
sont déterminés en utilisant la corrélation linéaire déterminée sur les standards.
206
Figure. 7.3 Diagramme Pb/238U vs 207
Pb/235U des analyses à la sonde ionique CAMECA IMS 1270 sur
monograin de zircon de l'intrusion rhyolitique de Tachkakacht
123
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
La granodiorite de Taouzzakt
Comparées aux données réalisées par Pupin (1992), les valeurs de Zr/Hf sont
similaires aux valeurs déterminées pour les zircons de granites calco-alcalins. Dans le
diagramme HfO2 vs Y2O3 de Pupin (1992) nos analyses sont réparties de façon homogène
entre les domaines calco-alcalins (CA) et calco-alcalins potassiques (CCA) (Fig. 7.4).
Figure 7.4 Diagramme de caractérisation géochimique HfO2 vs. Y2O3 d'après Pupin (1992) des grains de zircon
de la protrusion de Takhatert et de la granodiorite de Taouzzakt
124
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
125
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
générations Ms1. Les échantillons ont été prélevés dans des zones minéralisées
particulièrement riches.
Echantillon IM99-37 (Fig. 7.5a) : C'est la seule analyse qui permet de déterminer un
âge plateau à 560 ± 5 Ma (100% d'39Ar libéré, en 4 étapes). Cet âge est identique à l'âge
intégré calculé à 560 ± 9 Ma. L'âge déterminé sur l'étape de plus haute température est de
565 ± 8 Ma. L'âge minimum déterminé en base de spectre est de 540 ± 17 Ma.
126
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
Figure 7.5 Spectres d'âge des échantillons de la phase Qtz1 et âges in-situ en spot analyses dans un feuillet de
l'échantillon IM00-53
Les âges minima qui forment la base du spectre sont localisés entre 60 % et 80 %
d'39Ar libéré. Les âges déterminés sur ces étapes sont de 516 ± 26 Ma et de 528 ± 17 Ma.
127
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
Echantillon IM00-53 (Fig. 7.5c) : Un spectre en onze étapes a été réalisé sur un
monograin de muscovite. Les âges paliers se répartissent entre 482 ± 4 Ma et 573 ± 3 Ma. Le
spectre est perturbé et présente une forme en selle. Les âges les plus bas sont localisés entre
50 % et 80 % d'39Ar libéré. Ils sont compris entre ca. 555 Ma et ca. 559 Ma. L'étape de plus
haute température révèle un âge de 563 ± 5 Ma. L'âge intégré calculé sur l'ensemble du
spectre donne un âge de 560 ± 17 Ma (MSWD = 49).
Echantillon IM98-03 (Fig. 7.5d) : L'âge intégré défini sur cette muscovite est de 564 ±
06 Ma (MSWD = 82). Les quatre premières étapes de basse température correspondent à un
phénomène de dégazage superficiel. Le premier âge de basse température est de 520 ± 2 Ma.
Les étapes suivantes montrent un accroissement constant des âges jusqu'aux étapes de haute
température. L'âge déterminé sur la dernière étape est de 573 ± 1 Ma. Le spectre en escalier de
l'échantillon IM98-03 peut être interprété suivant un modèle de diffusion de Turner (1968),
caractéristique d'une perte d'40Ar lors d'un événement thermique plus récent.
Echantillon IM00-53 (âge spot) (Fig. 7.5e) : Une série de dix spot-fusion a été réalisée
sur un monograin de muscovite Ms1. Cette muscovite est composée de deux zones distinctes,
une première claire et translucide et une seconde, au centre du grain, grisâtre et recoupée par
des micro fractures. Les âges réalisés se répartissent en deux groupes. Le premier groupe est
constitué d'âges inférieurs à 551 Ma avec des marges d'erreur entre 17 et 19 Ma. Ils
correspondent aux analyses réalisées dans la zone fracturée de la muscovite. Le second
groupe est composé d'âge à ca. 561 Ma avec des marges d'erreur plus faibles (5 Ma). Les
spots sont localisés dans des zones "saines" du feuillet de muscovite. Les marges d'erreur
deviennent plus importantes lorsque le spot est localisé en bordure du feuillet.
128
Partie II Chapitre 7 Datations du plutonisme, du volcanisme et de la minéralisation
L'interprétation des perturbations observées sur les spectres 40Ar/39Ar peut également
apporter des informations importantes sur l'âge de l'épisode épithermal argentifère.
En effet, les âges les plus bas obtenus sur les paliers de moyenne température des
échantillons IM99-37 -58 et IM00-53 donnent des âges compris entre 516 ± 26 Ma et 555 ± 6
Ma. Ces observations peuvent être interprétées soit comme le résultat de dégazages successifs
dans différents domaines du cristal analysé, soit comme le résultat d'une perturbation
thermique et d'une recristallisation (Cheilletz et al., 1991). Ces effets ont pu être amplifiés par
39
le recul de Ar lors de l'irradiation, ce qui augmente l'aspect en scelle du spectre au niveau
des paliers de moyennes températures (Ruffet et al., 1991). Compte tenu de l'existence de
deux générations de muscovite ainsi que de la forme générale des spectres, une réinitialisation
partielle du système isotopique accompagnée d'une recristallisation des muscovites primaires
Ms1 semble une hypothèse acceptable. Cette hypothèse est validée par la répartition des âges
spots dans l'échantillon IM00-53. Les âges spots réalisés dans la zone claire et dans la zone
sombre de l'échantillon IM00-53 sont, respectivement, comparables aux âges intégrés et aux
âges de base de scelle déterminés pour les différents échantillons. Cette expérience permet de
confirmer l'existence de deux groupes d'âges.
Ainsi les analyses 40Ar/39Ar effectuées sur les muscovites Ms1 mettent en évidence un
événement thermique ultérieur qui perturbe le système isotopique K-Ar des minéraux. Les
températures de 350-180°C estimées pour l'épisode épithermal argentifère font de cet
événement un bon candidat à l'origine de cette perturbation. En outre, les âges de base de selle
compris entre 555 et 540 Ma, les analyses spots dans la zone fracturée comprises entre 534 et
551 Ma sont similaires à l'âge U/Pb des intrusions rhyolitiques à Imiter.
40
Les données Ar/39Ar et U/Pb établissent donc définitivement le synchronisme de
l'épisode rhyolitique et de l'épisode épithermal argentifère, déjà mis en évidence lors de
l'étude structurale (Partie II, chapitre 3).
129
Partie II Chapitre 8 Conclusions sur l'étude du gisement d'Imiter
Chapitre 8
Les études paragénétiques sont importantes car elles ont permis de distinguer deux
épisodes minéralisateurs, (i) un hydrothermalisme précoce développé dans l'auréole de
métamorphisme de contact des granodiorites et (ii) un épisode épithermal argentifère induit
par le volcanisme rhyolitique néoprotérozoïque terminal et qui se superpose spatialement au
premier épisode.
L'estimation, sur un certain nombre d'échantillons, des conditions de température, de
fugacité en soufre et en oxygène ainsi que du pH, a permis de fixer les caractéristiques
physico-chimiques correspondants au domaine de stabilité des assemblages minéralogiques
observés.
131
Partie II Chapitre 8 Conclusions sur l'étude du gisement d'Imiter
40
Deux méthodes géochronologiques différentes (U/Pb sur zircon, et Ar/39Ar sur
muscovite) ont été utilisées pour dater les intrusions plutoniques et volcaniques et la
minéralisation.
Les analyses U/Pb ont permis de dater précisément deux événements magmatiques
distincts à 572 ± 5 Ma (granodiorite de Taouzzakt) et 550 ± 3 Ma (rhyolite de Takhatert). Les
40
analyses Ar/39Ar ont fourni des résultats qui confirment l'existence des deux événements
magmatiques et hydrothermaux à Imiter. Il faut noter cependant que nous ne disposons pas
d'âge direct de la minéralisation, faute de minéraux analysables.
132
Partie II Chapitre 8 Conclusions sur l'étude du gisement d'Imiter
Figure 8.1 Colonne stratigraphique du Néoprotérozoïque et du cambrien dans le gisement d'Imiter. Les âges de
la colonne de gauche proviennent de la bibliographie (voir Fig. 1.3 et 1.4 pour les références).
Il ressort de cette synthèse lithostratigraphique les deux points suivants :
133
Partie II Chapitre 8 Conclusions sur l'étude du gisement d'Imiter
La figure 8.2 présente la synthèse des résultats obtenus sur le gisement d'Imiter. Nous
proposons d'expliquer la concentration exceptionnelle d'argent et de métaux associés, à l'aide
d'un modèle épithermal relié à la mise en place du magmatisme rhyolitique à la fin du
Néoprotérozoïque terminal. Dans ce modèle il n'est pas nécessaire de faire appel à un
phénomène de reconcentration supergène ultérieur (cémentation). En outre ce modèle ouvre
de nouvelles possibilités d'extension en profondeur du gisement comme c'est le cas à Acupan
(Cooke et al., 1996) par exemple, un autre gisement épithermal ou la colonne minéralisée est
continue sur 900 m de hauteur.
Les données des traceurs isotopiques en particulier le Re/Os suffisent pour envisager
une source profonde probablement mantellique pour l'essentiel du stock métal. Les porteurs
des minéralisations à Imiter sont des phases sulfurées en plus de l'argent natif. Dans l'état
actuel des données disponibles sur les isotopes du soufre, nous mettons en évidence deux
réservoirs possibles pour cet élément sans que l'on puisse, à l'aide de ce seul traceur, trancher
quant à l'importance relative des deux sources.
Néanmoins, l'importance des cristallisations de pyrite accompagnant le volcanisme
rhyolitique (protrusion de Takhatert et dyke de Tachkakacht) ainsi que la nature du
complexant probable de l'argent en solution (SH)-2 militent en faveur d'une origine en
majorité magmatique hydrothermale pour le soufre.
Nous proposons de classer le gisement d'Imiter dans la catégorie des gisements
épithermaux neutres (White et Hedenquist, 1995).
134
Partie II Chapitre 8 Conclusions sur l'étude du gisement d'Imiter
Figure 8.2 Nouveau modèle génétique, épithermal neutre, pour le gisement Ag-Hg d'Imiter.
Une retombée indirecte de ce travail est que, dans les gisements épithermaux les
études isotopiques du soufre ne sont pas déterminantes pour l'identification des sources. Ceci
est probablement dû au trop grand nombre de paramètres intervenants sur le fractionnement
de cet élément, en particulier : le dégazage des fluides, leurs températures et l'effet tampon
induit par les encaissants notamment, à Imiter, les black shales.
Les fluides de bassin chlorurés et sursalés mis en évidence par Baroudi et al. (1999) à
l'aide des inclusions fluides ont uniquement un rôle de refroidissant du système épithermal
permettant la précipitation des minéralisations. Aucune contribution d'éventuels sulfates
transportés par ces fluides n'est révélée par les analyses d'isotopes du soufre (tous les δ34SCDT
mesurés sont négatifs).
135
Partie II Chapitre 8 Conclusions sur l'étude du gisement d'Imiter
136
Partie III Etude géochronologique et géochimique des gisements (Pb-Zn-Cu-Ag-Au) de Bou Madine et
(Co-Ni-As-Au-Ag) de Bou Azzer : Implications géodynamiques
Partie III
Etude géochronologique
Implications géodynamiques
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
Chapitre 1
apport de la géochronologie
137
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
Figure 1.1 Carte géologique de la boutonnière de Bou Madine, et localisation des échantillons analysés (d'après
Paile, 1983).
138
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
lignée magmatique unique d'affinité calco-alcaline évoluant vers une composition tholéitique
pour les basaltes tardifs. L'extension latérale limitée des coulées volcaniques rend difficile
toute tentative de corrélation avec les formations volcaniques néoprotérozoïques tardives dans
Figure 1.2 Colonne stratigraphique de la région de Bou Madine (d'après Paile, 1983).
139
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
l'Anti-Atlas en général et avec la série de Ourazazate dans le Saghro en particulier. Cette série
volcanique est composée de quatre formations superposées (Fig. 1.2 ; Paile, 1983 ; Freton,
1988 ; Aït Saadi, 1992):
(1) Le conglomérat de base, à éléments anguleux et hétéro-métriques, affleure en
niveaux discontinus. Il est constitué de produits de démantèlement de la série
métasédimentaire et des granitoïdes. Il s'agit d'accumulations désordonnées le long des
grandes fractures N30°E (Aït Saadi, 1992).
(2) La formation de Tamerzaga est constituée d'un empilement de nappes
ignimbritiques et d'une coulée andésitique intercalaire. Elle se termine par la mise en place
d'intrusions rhyolitiques (chonolite) considérées par l'ensemble des auteurs comme
synchrones des minéralisations. La mise en place de la formation s'est effectuée en contexte
cisaillant sénestre selon la direction majeure N30°E (Freton, 1988)
(3) La formation d'Oufroukh regroupe des faciès détritiques (conglomérats, arkose,
grés), des tufs, des niveaux de rhyolites, d'andésites, de calcaires, de dolomies et de cherts. Ils
se mettent en place dans de petits bassins jalonnant les structures N30°E (Aït Saadi, 1992).
(4) La formation d'Ajouane-Aissa comporte des unités à dominante ignimbritique,
dont la mise en place est contrôlée par les structures N30°E. Elle se termine par la mise en
place des seuls basaltes en coulée de la région. La fin de ce cycle est marquée par des
injections de dykes et de stocks basiques et acides le long de fractures N-S à N160°E. La
transition avec la formation d'Oufroukh n'est pas nette. La série détritique pourrait se
poursuivre et constituer la base de la formation d'Ajouane-Aissa.
L'ensemble des faciès de la couverture néoprotérozoïque terminale est affecté par une
propylitisation régionale à quartz, albite, chlorite, calcite, rutile et épidote. Cette altération
régionale est oblitérée au niveau de la zone minière par une association à quartz et séricite
dominants.
140
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
sont structurés par une déformation hercynienne, qui réactive les accidents néoprotérozoïques
N30°E et développe des plis à grand rayon de courbure de même direction.
141
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
Les intrusions rhyolitiques (Fig. 1.4) montrent des formes ovoïdes ordonnées suivant
des fractures de direction N160°E à subméridienne. Ces roches sont composées de
phénocristaux de quartz et de feldspaths (0.5 à 1 mm) supportés par une mésostase
microlitique et microcristalline formée de lattes d'albite, d'orthose, et d'une pâte granuleuse à
quartz et chlorite comportant de rares cristaux de zircon et de minéraux opaques. Les
phénocristaux de quartz sont subautomorphes ou automorphes et présentent des golfes de
croissance typiques des faciès rhyolitiques. Ces cristaux sont souvent entourés d'auréoles de
142
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
Figure 1.3 Schéma synthétique montrant les relations entre les différents faciès mis en place avant les
minéralisations (d'après Dagallier et al., 1988)
Figure 1.4 Composition chimique d'une rhyolite stérile (chonolite), (BM99-1) d'une rhyolite minéralisée (Fe-Zn)
(BM99-4) et de la série volcanique néoprotérozoïque terminale dans le diagramme SiO2 vs Na2O + K2O de Le
Maître et al. (1989). Données compilées de Paile (1983), Freton (1988) et Ait Saadi (1992).
143
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
Les grains de zircons analysés sont de petite taille (50-200 µm). Tous les grains sont
automorphes, prismatiques ou aciculaires (ou des fragments de cristaux automorphes). Ils sont
homogènes ou présentent une fine zonation magmatique, certains grains sont perlés de fines
inclusions d'apatite et de thorite. Un grain présente des évidences de structures héritées
(Fig. 1.5). Les rapports isotopiques et les âges calculés sont présentés dans le tableau 1.1
Quatorze zircons ont été analysés. Pour les mêmes raisons que celles évoquées dans le
chapitre 7, nous privilégions ici les âges 206Pb/238U, qui présentent une erreur plus faible que
les âges 207Pb/235U. En utilisant cette grille de lecture, la distribution des âges 206Pb/238U dans
le diagramme concordia révèle deux groupes principaux (Fig. 1.5) constitués d'âges
concordants répartis respectivement entre 532 ± 2 Ma et 567 ± 5 Ma (âge moyen pondéré :
553 ± 16 Ma ; MSWD = 69 ; n = 6) et entre 581 ± 1 Ma et 582 ± 2 Ma (âge moyen pondéré :
582 ± 2 Ma ; MSWD = 0.2 ; n = 2). L'âge moyen pondéré 206Pb/238U, de 553 ± 16 Ma, calculé
sur les 6 âges concordants du premier groupe peut être considéré comme l'âge de mise en
place de la rhyolite. Les âges plus anciens représentent probablement des xénocrystaux hérités
des formations encaissantes (les granitoïdes et la série volcanique néoprotérozoïque
terminale). Il est possible d'identifier les âges déterminés sur la série volcanique
néoprotérozoïque terminale (ca. 580Ma), ainsi que les âges des déformations Pan-africaines
B1 (ca. 650 Ma, Leblanc et Lancelot, 1980). L'âge à 2223 ± 13 Ma, comme les âges
discordants, sont interprétés comme des âges de mélange entre un cœur hérité et une zone de
cristallisation postérieure. Ces âges ne peuvent donc être interprétés directement en terme
d'âge géologique.
144
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
Tableau 1.1 Résultats analytiques de la datation U/Pb par sonde ionique des zircons de l'intrusion rhyolitique de
Bou Maadine. Les rapports isotopiques du plomb sont corrigés du plomb commun. Les rapports U/Pb sont
déterminés en utilisant la corrélation linéaire déterminée sur les standards
206
Figure 1.5 Diagramme Pb/238U vs 207
Pb/235U des analyses à la sonde ionique sur monograin de zircon de
l'intrusion rhyolitique de Bou Madine
145
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
La méthode (sonde ionique CAMECA IMS3F) et les conditions d'analyse sont décrites
en détail dans la partie I chapitre 3 de ce manuscrit. Les analyses sont réalisées sur une
minéralisation massive de pyrite et de sphalérite dans une rhyolite (BM99-4). Les analyses
effectuées sur la pyrite donnent un δ34SCDT moyen de -0.6 ± 2.4 ‰. Ces résultats sont en
accord avec ceux de Dagallier (1989). L'étude in situ permet de montrer que la variation de
δ34SCDT des pyrites de l'ordre de 4 ‰ à l'échelle du gisement se retrouve à l'échelle de
l'échantillon (Tableau 1.2).
Les valeurs de δ34S calculées pour les fluides H2S à partir de l'ensemble des données
(Dagallier, 1989 ; cette étude) dans la gamme de températures (150°C à 360°C) déterminées
sur les inclusions fluides et sur la minéralisation (Ait Saadi, 1992 ; Abia et al., 1999)
présentent des valeurs moyennes comprises respectivement entre -1.9 ‰ et -1.3 ‰.
Ces valeurs se situent dans la gamme de variation des valeurs connues pour le soufre
magmatique (Ohmoto, 1974). Les valeurs de δ34S déterminées à Bou Madine sont interprétées
comme les témoins d'une source magmatique pour le soufre. Les variations de δ34S observées
sont interprétées comme l'expression de conditions locales particulières (mélange de fluides,
dégazage de SO2, variation locale de fO2) caractéristiques de ce type de système
hydrothermal, lié à un phénomène pulsatif de fracturation hydraulique/cicatrisation.
Tableau 1.2 Valeurs de δ34SCDT calculées de la littérature (pyrite, sphalérite, blende) (spectromètre à source
gazeuse, Dagallier, 1989) et les données nouvelles (pyrite) de cette étude (Sonde ionique CAMECA IMS 3F) à
Bou Madine.
146
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
Figure 1.6 Distribution des valeurs de δ34SCDT de la littérature et les données nouvelles de cette étude à Bou
Madine
147
Partie III Chapitre 1 Le gisement (Pb-Zn-Cu-Ag-Au) de Bou Madine : apport de la géochronologie
type porphyre cuprifère. A ce titre, une comparaison du gisement de Bou Madine avec le
gisement de Bingham (Cox 1986 ; Cunningham et al., 2001) ou mieux encore celui de Butte
(Marcoux et Jebrack, 1999) pourrait être très fructueuse. C'est peut-être dans un cadre
conceptuel différent qu'il convient de poser maintenant le problème de la récupération de l'or
à Bou Madine.
148
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
Chapitre 2
apport de la géochronologie
149
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
Figure 2.1 Carte géologique de la boutonnière de Bou Azzer - El Graara, (d’après Leblanc, 1975) et
localisation des échantillons analysés.
150
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
(5) La dernière expression magmatique reconnue dans la région est représentée par un
dyke doléritique. Les datations 40Ar/39Ar sur monograin de plagioclase réalisées sur ce dyke
donnent un âge de 196 ± 1.8 Ma (Sebei et al., 1991).
2-2. La minéralisation
Des études paragénétiques détaillées ont été réalisées par Maacha et al. (1998) et
Ennaciri (1995) (Fig. 2.2). Ces auteurs différencient deux types de minéralisation, des filons
parallèles au contact entre les serpentinites et la granodiorite et des filons sécants sur ce
contact. La succession minéralogique dans les deux cas est similaire, elle débute par une
séquence anté-arséniure composée de silicates (quartz, amphibole, chlorite et séricite)
associés parfois avec des sulfures. Le stade à arséniures et sulfoarséniures est composé
successivement de tri et diarséniures de Ni et de Co, de diarséniures de Fe, de sulfoarséniures
de Ni et de Co, de sulfoarséniures de Co-Fe et de sulfoarséniures de Fe. Cette séquence est
supportée par une gangue de quartz. La fin de ce stade à arséniure est marquée par la présence
de sulfures (sphalérite, galène, pyrite, chalcopyrite, bismuthinite, digénite et bornite), de
sulfosels (tétraédrite et tennantite) et localement de l’argent et du bismuth natifs. La mise en
place de fissures à barytine, la talcification des gangues carbonatées et l’altération des
arséniures en hétérogénite et érythrine constituent un événement tardi à post-minéralisation.
Les analyses multi-échelles, de la lame mince au district, révèlent une zonation de la
minéralisation. A l’échelle du gisement, cette zonation est illustrée par la succession du cœur
vers la périphérie des tri-diarséniures, sulfoarséniures et sulfures. Au niveau du district une
zonation Est-Ouest est exprimée par des concentrations prédominantes de safflorite et de
skuttérudite, à l’Ouest et de rammelsbergite au centre et à l’Est. Cette zonation régionale est
également exprimée dans l’étude des fluides, avec des températures d’homogénéisation plus
chaudes à l’Ouest. Ennaciri (1995) montre que la nature des fluides dans les inclusions
primaires et secondaires de la minéralisation de contact, est identique dans l’ensemble du
district (50 km). Il s’agit de saumures complexes (NaCl + CaCl2), de salinité entre 34,5% et
40,5% pds eq. La température d’homogénéisation moyenne est de 225°C à l’Ouest et de
190°C au centre et à l’Est. Ennaciri (1995) ; Ennaciri et al. (1997) expliquent cette zonation
régionale en terme de zonation verticale, les blocs Ouest et Centre-Est sont décalés de façon
sénestre par un accident post-précambrien, la variation des températures entre les blocs
suggère un niveau de mise en place plus profond à l’Ouest. Essaraj (1999) propose l’existence
de deux autres types de fluides. Un fluide tardif lié au dépôt de l’argent, dont la salinité est
151
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
plus variable, de 19% à 40.5% pds eq. NaCl + CaCl2. Les températures de piégeage sont
estimées entre 140°C et 250°C à des pressions de 20 à 100 MPa. Chronologiquement l’auteur
associe ce type de fluides aux fractures NE-SW et NW-SE respectivement précambriennes à
cénozoïques et hercyniennes. Le second type de fluides présente des salinités faibles, de 0.3 à
17% pds. Eq. NaCl + CaCl2 et des températures fortes, autour de 280°C. Ce fluide est supposé
antérieur au dépôt de la minéralisation. Une analyse par écrasement-lessivage semble indiquer
une origine profonde des saumures similaires aux saumures de bassins pétroliers (Essaraj,
1999).
Figure 2.2 Succession paragénétique des minéralisations de Bou Azzer, (Maacha et al., 1998). L'épaisseur du
trait est fonction de l'abondance relative du minéral.
Une étude isotopique du soufre et de l’oxygène a été réalisée sur les sulfures et les
sulfates et sur les gangues de dolomite et de calcite (Maacha et al., 1998). Les valeurs de δ34S
obtenues sont de -6.2 ‰ à + 4.9 ‰ pour les sulfures et de +25.7 ‰ pour les sulfates. Les
152
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
valeurs de δ18O sont de -2.9 ‰ à + 1.4 ‰ dans les carbonates. Ces résultats sont interprétées
par les auteurs comme la contribution isotopique de sulfates marins et d’eaux exogènes dans
le système minéralisateur. Ils proposent un modèle de mélange de fluides exogènes avec des
fluides endogènes chargés en arsenic. Dans ce modèle le cobalt et le nickel sont lessivés à
partir des péridotites serpentinisées du complexe ophiolitique de Bou Azzer.
Les profils d’altération au contact de la diorite et de la serpentinite ont été étudiés par
Leblanc et Lbouabi (1988). Les termes successifs de la zone d’altération sont : une diorite
saine, une diorite cataclasée, une cataclasite metasomatisée en K-Na, une rodingite
(métasomatisme Mg-Ca), une chloritite et une serpentinite. Les minéraux d’altération
prépondérants sont les phyllosilicates (chlorites) et les feldspaths (adulaire et albite). Les
chlorites sont bien développées dans la zone à arséniures. Essaraj (1999) identifie deux
familles de chlorites : des brunsgivites et des brunsgivites-clinochlores. Cet auteur décrit
également plusieurs générations d’albite ainsi que des filonnets d’adulaire et d’albite tardifs
recoupant la granodiorite.
153
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
Les modèles génétiques proposés pour le district de Bou Azzer sont nombreux. Les
premiers modèles de Jouravsky (1952) et Goloubinov (1956) proposent un lessivage du Co et
Ni à partir des serpentinites par des fluides hydrothermaux lors de l’orogenèse hercynienne.
Figure 2.3 Modèle génétique de la minéralisation selon Leblanc et Lbouabi, (1988). Succession chronologique
des altérations qui ont affecté la diorite quartzite de Bou Azzer.
Kroutov et al. (1989) propose un modèle similaire, mais il associe les minéralisations au cycle
de distension adoudounien. Leblanc (1975) identifie trois stades de minéralisations (Fig. 2.3):
154
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
2-6. Le filon 7
Figure 2.4 Carte structurale simplifiée et phases de déformation observées dans le filon 7. (d’après Maacha et
al., 1998)
155
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
présente une puissance variant entre 1 et 20 m et il est reconnu jusqu’à une profondeur de
-340 m. Le levé de surface permet à Ennaciri (1995) d’identifier la chronologie de dépôts
suivante : (1) ankérite issue en partie de la carbonatation des serpentines ; (2) quartz laiteux ;
(3) calcite spathique blanche ; (4) barytine ; (5) quartz translucide. Selon les travaux de
Goloubinov (1956), la minéralisation est soit parallèle au contact, soit associée à des accidents
de direction NE-SW. La minéralisation s’exprime en surface sous forme d’érythrine dans les
fentes à quartz orientées N115°E et d’hétérogénéite dans les fentes de quartz laiteux
orientéesN40°E. Ces deux familles sont exprimées au niveau -145 m respectivement par des
arséniures associés à des quartz et carbonates orientés N75°E et par des fentes plus tardives de
direction N45°E (Fig. 2.4).
L’inventaire minéralogique et la succession paragénétique selon Ennaciri (1995) sont
présentés sur la figure 2.5. Les minéraux prépondérants sont la skuttérudite et la
(clino)safflorite. La succession chronologique est la suivante (du plus ancien au plus récent) :
- le stade A est constitué de triarséniures (rare rammelsbergite), diarséniures [(clino)safflorite,
skuttérudite], puis sulfoarséniures (arsenic et cobaltite).
- Le stade B est constitué de sulfures. Ils remplissent les fissures dans les arséniures et les
sulfoarséniures.
Figure 2.5 Succession minéralogique du filon 7. (Ennaciri, 1995). Ce filon renferme apparemment des
paragenèses correspondant au stade à arséniures de Maacha et al., 1998
156
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
Tableau 2.1 Valeurs de δ34SCDT de la littérature (pyrite, chalcopyrite, barytine) (Maacha et al., 1998) et les
données nouvelles (pyrite, galène) de cette étude
La température calculée à partir des valeurs de δ34SCDT du couple barite-pyrite est de 265 ±
10°C. Cette valeur est équivalente aux températures de cristallisation déterminées sur des
157
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
chlorites de Bou Azzer. Elle n’est toutefois acceptable que si ces deux minéraux sont
cogénétiques. La succession paragénétique à Bou Azzer apparaît très complexe et limite
l'utilisation des traceurs isotopiques. Par exemple, le calcul de la température de formation à
partir des valeurs de δ34SCDT du couple pyrite-chalcopyrite de l’échantillon 3043 (Maacha et
al., 1998) est impossible parce que le résultat est négatif. On peut donc en conclure que soit
les deux minéraux ne sont pas cogénétiques, soit que les conditions de cristallisation n'ont pas
permis d’atteindre l’équilibre isotopique.
Figure 2.6 Distribution des valeurs de δ34SCDT de la littérature et les données nouvelles de cette étude à Bou
Azzer
158
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
(1) Les chlorites : Les analyses à la microsonde électronique des chlorites montrent que ce
sont principalement des Mg-chamosites selon la classification de Bayliss (1985) (Fig. 2.8A)
La formule structurale moyenne calculée sur 11 atomes d’oxygène est donnée dans la figure
2.8B. Les températures de cristallisation de ces chlorites calculées avec la formule de
Cathelineau (1988) montre une gamme de température de 250°C à 274°C. Ces températures
sont similaires aux températures basses que nous avons calculées à partir des analyses
réalisées par Essaraj (1999) sur les mêmes minéraux (Fig. 2.8C). Les deux groupes
cristallographiques de chlorites identifiés par Essaraj (1999) (Mg-chamosite et Mg-
chamosite/clinochlore) forment respectivement deux modes de températures distincts à
270°Cet 320°C. La mise en évidence de ces deux groupes de chlorites et des températures
associées suggèrent l’existence de deux événements hydrothermaux distincts. La chlorite
étudiée dans cette étude appartient à l’événement dont la gamme des températures présente un
mode à 270°C. Les chlorites constituant la roche présentent une cristallochimie homogène et
une température de cristallisation correspondant à l’événement hydrothermal tardif proposé
par Essaraj (1999).
(2) Les feldspaths potassiques : Les analyses MEB réalisées sur les feldspaths
potassiques présentent des compositions de type adulaire. Ces feldspaths adulaires ont été
séparés pour datation 40Ar/39Ar.
159
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
160
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
40
Figure 2.9 Datation Ar/39Ar des adulaires de la chloritite associée à la minéralisation à (As Fe Co Ni) du
filon7 et des muscovites du filon de quartz qui recoupe la rhyolite au niveau de la mine de Tamdrost.
161
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
mise en place des sills permet donc de déterminer l'âge maximum de la minéralisation. La
roche est composée de phénocristaux d'orthose et d'albite dans une mésostase trachytique
formée de lattes d'albite, d'orthose, et d'une pâte granuleuse à quartz comportant de rares
cristaux de zircon et de nombreux minéraux opaques. Les phénocristaux d'orthoses sont
séricitisés. Le trachyte présente des teneurs en SiO2 et Na2O + K2, respectivement de 64,33%
et 8,07% ainsi qu'un enrichissement en éléments incompatibles caractéristiques des
associations calco-alcalines à alcalines.
Les résultats analytiques et les âges déterminés sur les 12 grains de zircons analysés
sont présentés dans le tableau 3.1 et dans le diagramme concordia (Fig. 2.10). Les âges
206
individuels Pb/238U obtenus varient entre 533 ± 2 Ma et 1631 ± 16 Ma avec trois modes
principaux à 580 Ma, 640 Ma et 730 Ma. Le mode à 580 Ma est constitué de 7 âges
concordants entre 533 ± 2 Ma et 592 ± 11 Ma. Le mode à 640 Ma est constitué de 3 âges
concordants entre 637 ± 8 Ma et 660 ± 6 Ma. Le troisième mode est constitué d'un seul âge à
731 ± 8 Ma. Enfin une analyse a donné un âge à 1631 ± 16 Ma. Les âges obtenus présentent
une grande dispersion mais sont tous concordants. Les zircons analysés sont monophasés et
par conséquent chaque grain de zircon est représentatif d'un événement géologique. Ce
résultat est cohérent avec la nature volcanique de la roche, qui a incorporé une partie de la
croûte sous-jacente lors de sa mise en place. L'âge le plus jeune représente donc l'âge de mise
en place du trachyte d'Aghbar à 533 ± 2 Ma. Cet âge cambrien inférieur est similaire à l'âge
Tableau 2.2 Résultats analytiques de la datation U/Pb par sonde ionique des zircons du Trachyte d'Aghbar. Les
rapports isotopiques du plomb sont corrigés du plomb commun. Les rapports U/Pb sont déterminés en utilisant
la corrélation linéaire déterminée sur les standards
162
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
Figure 2.10 Diagramme 206Pb/238U vs 207Pb/235U des analyses à la sonde ionique sur monograin de zircon du sill
trachyte d'Aghbar.
déterminé par Ducrot et Lancelot (1977) pour le volcanisme alcalin du Jbel Boho à 534 ± 10
Ma (U/Pb par dilution isotopique de population de zircons). Il est également cohérent avec les
datations réalisées par Compston et al. (1992) et par Landing et al. (1998) sur la "Série Lie de
Vin" stratigraphiquement plus jeune et datée à 522 ± 2 Ma (U/Pb sur monograin de zircon par
sonde ionique SHRIMP).
163
Partie III Chapitre 2 Le gisement (Co-Ni-As-Au-Ag) de Bou Azzer : apport de la géochronologique
la datation du trachyte d'Aghbar donne un âge de 533 ± 2 Ma, qui est conforme avec la
position stratigraphique du sill intercalé dans la formation du Cambrien inférieur. Cet âge est
similaire à l'âge mesuré sur les trachytes du Jbel Boho (Ducrot et Lancelot, 1977).
Du point de vue des minéralisations Co-Ni de Bou Azzer, celles-ci recoupent le
trachyte d'Aghbar, ce qui implique donc que ces minéralisations ont un âge maximum de
533 ± 2 Ma
Les datations réalisées sur les adulaires de la chloritite BA99-8 du filon 7 donnent un
âge plateau 40Ar/39Ar à 218 ± 8 Ma (Norien). Ce résultat est important car il est voisin de la
datation Pb/Pb sur brannérite obtenue par Ledent (1960) à 240 ± 10 Ma. Il y a donc à Bou
Azzer un événement hydrothermal minéralisateur plus récent, dont l'âge serait triasique.
Il est important également de rappeler l'existence au niveau de l'Anti-Atlas d'autres
occurrences de cobalt/nickel, dont la chronologie relative suggère un âge fini-triasique.
- L'indice de Tizi Asdrem situé sur le bord oriental du Siroua, le long d'une faille NE-SW,
jalonné par un dyke de dolérite. La faille est interprétée comme un décrochement
précambrien, ayant rejoué en faille normale (Choubert et al., 1974).
- L'indice à BPGC avec arséniures et sulfoarséniures de Ni-Co à gangue de sidérite de Mefis
(Jbel Debouah) dans le Tafilalet. Les filons de minéralisations remplissent les fractures
hercyniennes, qui sont également empreintées par des dykes de dolérites (Leblanc, 1975).
- Les occurrences de Ni-Co du Méfis sont situées sur le prolongement oriental de l'Accident
Majeur de l'Anti-Atlas (comme Bou Azzer), et présente une liaison géométrique étroite avec
les dolérites supposées jurassiques (Agard 1958 ; Jouravsky 1952)
En outre, des données récentes mettent en évidence un événement hydrothermal
hercynien daté à 300 Ma (40Ar/39Ar sur muscovite) dans la zone minéralisées aurifère
d'Iourirne, (Anti-Atlas central) (Gasquet et al. 2001b).
Ainsi le gisement de Bou Azzer apparaît clairement comme le résultat d'une histoire
polyphasée se déroulant de la base du Cambrien inférieur jusqu'au Trias. Cette histoire
polyphasée explique la complexité et l'enchaînement des successions paragénétiques dont le
décryptage constitue le préalable indispensable à toute étude géochimique fine. Du point de
vue géodynamique un des éléments caractéristique de l'Anti-Atlas central est la présence d'un
accident lithosphérique majeur (Accident majeur de l'Anti-Atlas) dont les rejeux répétés
(Guiraud et al., 2000) depuis le Protérozoïque pourraient expliquer la permanence des
phénomènes métallogéniques dans cette région.
164
Partie III Chapitre 3 Conclusions : implications géodynamiques
Chapitre 3
165
Partie III Chapitre 3 Conclusions : implications géodynamiques
Figure 3.1 Nouvelles colonnes stratigraphiques du Néoprotérozoïque et du Cambrien dans l'Anti-Atlas central et
oriental.
166
Partie III Chapitre 3 Conclusions : implications géodynamiques
l'épisode volcanique représentée par les rhyolites est datée à 550 Ma. Ce résultat souligne la
situation émergée, jusqu'au Cambrien moyen, de l'Anti-Atlas oriental tandis que l'Anti-Atlas
central est progressivement envahi par la mer infracambrienne.
2- les datations sur les trachytes d'Aghbar (cette étude) et celles du Jbel Boho (Ducrot
et Lancelot, 1977) indiquent clairement qu'une partie de la formation d'Adoudou est
précambrienne.
3- les âges hérités sont identiques à ceux obtenus sur les granitoïdes de l'Anti-Atlas
central (Ducrot et Lancelot 1977 ; Leblanc et Lancelot 1980 ; Mifdal et Peucat 1985) et
permettent de mieux contraindre l'âge des deux phases de déformation panafricaines.
Il apparaît à la lecture de cette colonne stratigraphique synthétique que la transition
Précambrien/Cambrien est caractérisée par une interaction d'événements magmatiques
(volcanisme effusif calco-alcalin et alcalin), tectoniques (distension généralisée),
sédimentaires (transgression infracambrienne puis cambrienne) et métallogéniques (gisements
de Bou Madine, d'Imiter et pro parte Bou Azzer en particulier). Il est donc important à ce
stade de proposer un scénario de l'évolution géodynamique globale de l'Anti-Atlas, qui prenne
en compte l'ensemble de ces observations nouvelles.
Après avoir examiné en détail les résultats de notre étude sur la limite
Précambrien/Cambrien dans l'Anti-Atlas, nous allons à présent en examiner les conséquences
géodynamiques sur une plus grande échelle de temps cette fois-ci, allant du Néoprotérozoïque
moyen jusqu'au méso-Cénozoïque. Le but de ce paragraphe n'est bien évidemment pas de
refaire la géologie du Maroc mais d'apporter quelques éléments nouveaux permettant
d'enrichir le débat actuel.
Le scénario d'une évolution géodynamique probable pour la fin du Protérozoïque est
présenté sur la figure 3.2. Deux étapes seulement ont été schématisées dans un but de
simplification. Ce modèle d'évolution concerne essentiellement une section NNW-SSE
passant dans l'Anti-Atlas oriental. Il est évident que, dans le détail ainsi que sur l'ensemble de
l'Anti-Atlas, ces schémas mériteraient d'être multipliés et complétés.
167
Partie III Chapitre 3 Conclusions : implications géodynamiques
3-2-2. La deuxième étape de cette évolution n'est pas développée sur la figure 3.2. Elle
concerne le serrage de l'ensemble marge continentale - bassin d'avant arc, au cours du
Panafricain dans l'intervalle de 685 - 610 Ma (phases B1 et B2 des auteurs). Certes, il est à
noter que dans l'Anti-Atlas central, ce serrage aboutit à l'obduction de la croûte océanique et à
la formation de l'écaille ophiolitique de Bou Azzer (Saquaque et al. 1989a ; Ennih et Liegeois
2001). Le résultat d'ensemble de cette phase de convergence est un important
surépaississement crustal de la bordure NW du craton Ouest africain.
168
Partie III Chapitre 3 Conclusions : implications géodynamiques
al. 1999). Ces sédiments sont, en outre, associés à des intercalations de volcanisme alcalin
d'origine profonde (Aghbar, Jbel Boho).
Ces schémas géodynamiques (Fig. 3.2) permettent de coupler, dans l'espace et dans le
temps, l'évolution magmatique de cette marge continentale avec les phénomènes
métallogéniques.
169
Partie III Chapitre 3 Conclusions : implications géodynamiques
170
Références bibliographiques
Références
bibliographiques
Références bibliographiques
Références bibliographiques
171
Références bibliographiques
172
Références bibliographiques
Becker H, 2000 Re-Os fractionation in eclogites and blueshists and the implications for
recycling of oceanic crust into the mantle. Earth and Planetary Science Letters 177:
287-300.
Benisek A et Finger F 1993. Factors controlling the development of prism faces in granite
zircons : a microprobe study. Contrib. Mineral. Petrol. 114: 441-451.
Benssaou M et Hamoumi N 1999. Paléoenvironnements et minéralisations de l'Anti-Atlas
occidental marocain au Cambrien précoce. Chronique de la Recheche Minière, 536-
537, 113-121.
Benziane F, Prost AE et Yazidi A 1983. Le passage du Précambrien au Cambrien précoce
volcanique et sédimentaire de l'Anti-Atlas oriental ; comparaisons avec l'Anti-Atlas
occidental. Bull. Soc. Géol. France, 7, XXV, 4, 549-556.
Bertrand H et Prioton JM 1975. Les dolérites marocaines et l'ouverture de l'Atlantique :
étude pétrologique et géochimique. Thèse, Univ. Lyon. p.321.
Bertrand H, Dostal J et Dupuy C 1985. Geochemistry of Early Mesozoic tholeiites from
Morocco. Earth and Planetary Sciance Letter. 58, p. 228-239.
Birck JL et Allègre CJ 1994. Contrasting Re/Os magmatic fractionation in planetary basalts.
Earth and Planetary Science Lettres, 124, 139-148.
Birck JL, Roy-Barman M, et Campas F, 1997. Re-Os isotopic measurments at the
femtomole level in natural samples. Geostandards Newsletter, 20, p.139-148.
Bodinier JL, Dupuy C et Dostal J 1984. Geochemistry of Precambrian ophiolites from Bou
Azzer, Morocco. Contributions to Mineralogy and Petrology, 78, 43-50.
Brandon AD, Walker RJ, Morgan JW, Norman MD and Prichard HM, 1998 Coupled
186
Os and 187Os evidence for Core-Mantle Interaction. Science, 280, 1570-1573.
Cahen L, Choubert G, Delhal J et Ledent D 1965. Premières données sur l'âge du
Précambrien ancien de l'Anti-Atlas (Maroc) par la méthode Strontium-Rubidium.
Notes Serv. Géol. Maroc, 24, 183, p. 63-74.
Caironi V, Colombo A, Tunesi A et Gritti C 2000. Chemical variations of zircon compared
with morphological evolution during magmatic crystallization : an example from the
Valle del Cervo Pluton (Western Alps). Eur. J. Mineral. 12, 779-794.
Camara LS 1993. Etude cartographique, pétrographique, minéralogique et géochimique des
granitoïdes du district minier d’Imiter. Thèse de troisième cycle. Univ. Marrakech, p.
143.
Carlson RW et Irving AJ 1994. Depletion and enrichment history of subcontinental
lithospheric mantle: An Os, Sr, Nd, and Pb isotopic study of ultramafic xenoliths from
173
Références bibliographiques
the northwestern Wyoming Craton. Earth and Planetary Science Letters 126: 457-
472.
Cathelineau M 1988. Cation site occupancy in chlorites and illites as a function of
temperature. Clay Minerals, 23, 471-485.
Cathelineau M, Nieva D, 1985. A chlorite solid solution geothermometer. The Los Azufres
geothermal system (Mexico). Contrib. Mineral. Petrol., 91, p. 235-244.
Chaabane A 1991. Les roches vertes du Protérozoïque supérieur de Khzama (Siroua Anti-
Atlas, Maroc). Un exemple précambrien d'ophiolite d'avant-arc formée en contexte
transformant. Thèse d'état ès-Sciences, Univ. Marrakèch; 570p.
Charlot R 1978. Caractérisation des événement éburnéens et pan-africains dans l’Anti-Atlas
marocain. Apport de la méthode géochronologique Rb-Sr. Thèse Université Rennes I,
p. 220.
Chaussidon M, Demange JC, 1988. Instrumental mass fractionation in ion microprobe
studies of sulphur isotopic ratios. In : Secondary Ion Mass Spectrometry VI (eds. A.
Benninghoven, A.M. Huber and H.W. Werner), John Wiley & Sons, 937-940.
40
Cheilletz A, Ruffet G, Marignac C, Kolli O, Gasquet D et Feraud G 1991. Ar/39Ar
dating of shear zones in the Variscan basement of Greater Kabylia (Algeria). Evidence
of an Eo-Alpine event at 128 Ma (Hauterivian-Barremian boundary) : geodynamic
consequences. Tectonophysics, 306, 97-116.
Cheilletz A, Ruffet G, Marignac C, Kolli O, Gasquet D, Feraud G (1999) 40Ar/39Ar dating
of shear zones in the Variscan basement of Greater Kabylia (Algeria). Evidence of an
Eo-Alpine event at 128 Ma (Hauterivian-Barremian boundary) : geodynamic
consequences. Tectonophysics, 306, 97-116.
Chen JH, Papanastassiou DA et Wasserburg GJ 1998. Re/Os systematics in chondrites
and the fractionation of the platinum group elements in the early solar system.
Geochimica et Cosmochimica Acta. 62, 3379
Cherniak DJ et Watson EB 2001. Pb diffusion in zircon. Chemical Geology, 172, 5-20.
Cherniak DJ, Hanchar JM et Watson EB 1997. Rare-earth diffusion in zircon. Chemical
Geology. 134, 289-301.
Chesley JT, Integrative geochronology of Ore Deposits : New insights into the duration and
timing of hydrothermal circulation. In Application of radiogenic isotopes to ore
deposit research and exploration. Reviews in Economic Geology, Vol 12, 5, 115-141.
174
Références bibliographiques
175
Références bibliographiques
176
Références bibliographiques
Dunning GR et Hodyck JP 1990. U/Pb zircon and baddeleyite ages for the Palissades and
Gettysburg sills of the northeastern United States: Implications for the age of the
Triasic/Jurassic boundary. Geology 18, p. 795-798.
Dupuy C, Marsh J, Dostal J, Michard A et Testa S 1988. Astenospheric and lithospheric
source for mesozoic dolerite from Liberia (Africa) : trace element and isotopic
evidence. Earth and Planetary Science Letters, 87, 100-110.
Duthou JL, Amberger A et Lassere M 1976. Résultats graphiques et interprétations de
mesures isotopiques di plomb de galène et de minéraux oxydés du Maroc. Mémoires
hors série Soc. Géol. Fr., 7, 221-226.
El Aouli H, Gasquet D et IkenneM 2001. Le magmatisme basique de la boutonnière
d'Igherm Anti-Atlas occidental, Maroc) : un jalon des distensions néoprotérozoïques
sur la bordure nord du craton ouest-africain. Bull. Soc. géol. France, 172, 3, pp. 309-
317.
El Boukhari A 1991. Magmatisme et metasédiments associés au Protérozoïque supérieur de
la zone N'Kob (Siroua SE, Anti-Atlas Central, Maroc). Une ophiolite formée et mise
en place sur la marge du craton Ouest-Africain. Thèse Doct. ès-Science, Univ. Cadi
Ayyad, Marrakèch, Maroc, 511p.
Eldridge CS, Compston W, Williams JL, Walshe JL, Both RA, 1987. In situ
34
microanalysis for S/32S ratios using the ion microprobe SHRIMP. Internat. Jour.
Mass. Spectrometry Ion Processes. 76, 65-83.
Ennaciri A 1995. Contribution à l’étude du district à Co, As, (Ni, Au, Ag) de Bou Azzer,
Anti-Atlas (Maroc). Données minéralogiques et géochimiques; études des inclusions
fluides. Thèse Université d’Orléans, 238 p.
Ennaciri A, Barbanson L et Touray JC 1997. Brine inclusions from the Co-As-(Au) Bou
Azzer district, Anti-Atlas Mountains, Morocco. Econ. Geol., 92, 360-367.
Ennaciri A, Barbanson L, Lancelot J et Touray JC 1996. Distribution et âge de mise en
place des minéralisations aurifères du district de Bou Azzer (Anti-Atlas, Maroc) 16ème
RST, Orléans, France.
Ennih N et Liègeois JP sous presse. The Morocco Anti-Atlas: the West African craton
passive margin with limited Pan-African activity. Implication for the northern limit of
the craton. Precambrian Research
Essaraj S 1999. Circulations fluides associées aux minéralisations argentifères de l'Anti-Atlas
central exemple des gisements de Zgounder (Ag-Hg) et Bou Azzer (Co-Ni-As-Au-
Ag). Thèse d'état, Université Cadi Ayyad, Marrakech, 358p.
177
Références bibliographiques
Esser B et Turkian K 1993. The Osmium isotopic composition of the continental crust.
Geochimica et Cosmochimica Acta 57: 3093-3104.
Evensen NM, Hamilton PJ et O'Nions RK 1978. Rare earth abundances in chondritic
meteorites. Geochim. Cosmochim. Acta., 42, 1199-1212.
Fekkak A, Pouclet A, Ouguir H, Badra L et Gasquet D 1999. Le groupe du
Néoprotérozoïque inférieur de Kelaat Mgouna (Saghro, Anti-Atlas, Maroc) : témoin
d'un stade précoce de líextension pré-panafricaine. Bull. Soc. géol. France, 170, p.
789-797.
Fiechtner L, Friedrichsen H et Hammerschmidt K 1992. Geochemistry and
geochronology of early Mesozoic tholeiites from Central Morocco. Geol Rundsch.,
81p. 45-62.
Field CW et Fifarek R 1985. Light stable isotope systematics in the epithermal environment.
In Geology and geochemistry epithermal systems. Berger BR et Bethke PM. Review in
Economic Geology, 2, 99-128.
Field CW, Gustafson LB, 1976. Sulfur isotopes in the porphyry copper deposit at El
Salvador, Chile : Economic Geology, v.71, p. 1533-1548.
Field CW, Lombardi G, 1972. Sulfur isotopic evidence for the supergene origin of alunite
deposits, Tolfa district, Italy : Mineralium Deposita, v.7, p. 113-125.
Field CW, Rye RO, Dymond JR, Whelan F, Senechal RG, 1983. Metalliferous sediments
of the East Pacific ; in Shanks, W.C., III (eds), Cameron Volume on Unconventional
Mineral Deposits : Society of Mining Engineer, New York, p. 133-156.
Field CW, Sakai H, Ueda A, 1984. Isotopic contraints on the origin of sulfur in oceanic
igneous rocks ; in Syngenesis and epigenesis in the formation of mineral Deposits :
Springer – Verlag, Berlin-Heidelberg, p.573-589.
Foster JG, Lambert DD, Frick LR et Maas R 1996. Re-Os isotopic evidence for genesis of
Archean nickel ores from uncontaminated komatiites. Nature, 382, 703-706.
Francin B et Marion P 1984. Rapport d'étude minéralogiques; gisement de Bou Madine
(Maroc). Rapport BRPM, p.53.
Freton R 1988. Contribution à l'étude métallogénique du district de Bou Madine (Anti-Atlas,
Maroc) : environnement tectonique et concentrations epithermales BPGC à Ag et Au.
Thèse INPL, Nancy, p.246.
Freydier C, Ruiz J, Chessley J, McCandless T et Munizaga F 1997. Re/Os isoptope
systematics of sulfides from felsic igneous rocks: Application to base metal porphyry
mineralization in Chile. Geology 25: 775-778.
178
Références bibliographiques
Gammons CH et Barnes HL 1989. The solubility of Ag2S in near neutral aquous sulfide
solutions at 25 to 300°C. Geochimica et Cosmochimica Acta, 53, 279-290
Garcia D 1979. Altérations de roches encaissantes de gisements cobaltiféres du district de
Bou Azzer (Anti-Atlas, Maroc). Thèse Ecole des Mines de St Etienne, 127 p.
Gasquet D 1991. Genèse d'un pluton composite tardi-Hercynien. Le massif du Tichka, Haut
Atlas occidental (Maroc). These d'état, Nancy 1, 413p.
Gasquet D, Roger J, Chalot-Prat F, Hassenforder B, Baudin T, Chèvremont P,
Benlakhdim A et Razin P 2000b. Notice explicative, carte géologique Maroc (1/50
000), feuille de Tamazrar. MEM/BRGM, p78.
Gasquet D, Roger J, Chalot-Prat F,, Hassenforder B, Baudin T, Chèvremont P,
Cocherie A, Benlakhdim A, Zouhair M, Cheilletz A et Archibald D 2001a. Le
magmatisme polycyclique des boutonnières de la Tagragra d'Akka et du Kerdous-
Tafeltast (Anti-Atlas occidental, Maroc). Abstract Magmatisme Metamorphisme et
minéralisations associées, Marrakèch, Maroc, p 30.
Goloubinov R 1956. Prospection géochimique du cobalt à Bou Azzer (Maroc). Sci. Terre, 4,
Nancy, 211-226.
Gorjan P, Veevers JJ et Walter MR 2000. Neoproterozoic sulfur-isotope variation in
Australia and global implications. Precambrian Reserch, 100, 151-179.
Gradstein F et Ogg J 1996. A Phanerozoic time scale. Episodes 19: 3
Grotzinger JP, Bowring SA, Saylor BZ et Kaufman AJ 1995. Biostratigraphic and
geochronologic constraints on early animal evolution. Science, 270, 598-604.
Guillou JJ, Monthel J, Picot P, Pillard F, Protas J et Samama JC 1985. L'imitérite,
Ag2HgS2, nouvelle espèce minérale; propriétés et structure cristalline. Bull. Minéral.,
108,pp. 457-464.
Guillou JJ, Monthel J, Samama JC et Tijani A 1988. Morphologie et chronologie relative
des associations minérales du gisement mercuro-argentifères d'Imiter (Anti-Atlas,
Maroc). Notes et Mém. Serv. Géol. Maroc, Rabat 44, 334,pp. 215-228.
Guiraud R, Doumnang Mbaigane JC, Carretier S et Dominguez S 2000. Evidence for a
6000 Km lengh NW-SE striking lineament in northern Africa: the Tibesti Lineament.
J. of Geol. Soc. London, Vol 157, p. 897-900.
Gyssel JL 1982. Volcanisme et minéralisations cuprifère à Issougrid (région de Ouarzazate,
Maroc) Thèse de 3ème cycle, Paris XI Orsay.
179
Références bibliographiques
Hafid A, Sagon JP, Fonteilles M et Moutte J 1998. Existance de deux séries tholéiitiques
distinctes parmi les dolérites du Protérozoïque supérieur de la boutonnière d'Irhem
(Anti-Atlas occidental, Maroc). C.R.Acad. Sci. Paris, 236, p. 21-26.
Haillwood EA et Mitchel JG 1971. Paleomagnetic and radiometric dating results from
jurassic intrusions in South Morocco. Geophys. J. R. Astron. Soc., 24, 351-364.
Halden NM et Hawthorne FC 1993. The fractal geometry of oscillatory in crystals :
Application to zircon. Am.Mineral. 78, 1113-1116.
Hall HT 1967. The pearceite and polybasite series. American Mineralogist, 52, 1311-1321.
Hanchar JM et Miller CF 1993. Zircon zonation patterns as revealed by
cathodoluminescence and backscattered electron images: Implications for
interpretation of complex crustal histories. Chemical Geology 110: 1-13.
Hartman LA, Takehara L, Leite JAD, McNaughton NJ et Vasconcellos MAZ 1997.
Fracture sealing in zircon as evaluated by electron microprobe analyses and back-
scattered electron imaging. Chemical Geology 141: 67-72.
Hassenforder B 1987. La tectonique panafricaine et varisque de l'Anti-Atlas dans le massif
du Kerdous (Maroc). thèse Sci. Univ. L. Pasteur, Strasbourg., p. 249.
Hauri EK et Hart SR 1993. Re/Os isotope systematics of HIMU and EMII oceanic island
basalts from the south Pacific Ocean. Earth and Planetary Science Letters 114, 353-
371.
Hauri EK et Hart SR 1997. Rhenium abondances and systematics in oceanic basalts.
Chemical geology, 139, 185-205.
Hefferan K, Karson J et Saquaque A 1991. Proterozoic collisional basins in a Pan-African
suture zone, Anti-Atlas Mountains, Morocco, Precambrian Research, 54, p. 295-319.
Hoefs J, 1980. Stable isotope Geochemistry, Second Edition : Springer-Verlag, Berlin-
Heidelberg, New York, 208p.
Horan MF, MorganW, Grauch RI, Coveney RM, Murowchick JB et Hulbert LJ 1994.
Rhenium and Osmium isotopes in Black shales and Ni-Mo-PGE-rich sulfide layers,
Yukon Territory, Canada, and Hunan and Guizhou provinces, China. Geochimica et
Cosmochimica Acta 58(257-265).
Horan MF, MorganW, Walker RJ, and Grossman JN, 1992 Rhenium-Osmium isotope
constraints on the age of iron meteorites. Science, 255, 1118-1121.
Hoskin P 2000. Patterns of chaos :Fractal statistics and oscillatory chemistry of zircon.
Geochimica et Cosmochimica Acta 64, (11) 1905-1923..
180
Références bibliographiques
181
Références bibliographiques
182
Références bibliographiques
Lee JKW, Williams IS et Ellis DJ 1997. Pb, U and Th diffusion in natural zircon. Nature,
390, 159-161.
Leistel JM et Qadrouci A 1991. Le gisement argentifère d'Imiter (Protérozoïque supérieur
de l'Anti-Atlas, Maroc). Contrôles des minéralisations, hypothèse génétique et
perspectives pour l'exploration. Chron Rech Min 502: 5-22.
Leroux H, Reimold WU, Koeberl C, Hornemann U et Doukhan JC 1999. Experimental
shock deformation in zircon : a transmission electron microscopic study. Earth and
Planetary Science Letters. 169, 291-301.
Levchenkov OA, Rizvanoa NG, Maslenikov AV, Makeev AF, Bezmen N et Levsky LK
1996. Kinetics of Pb and Loss from Metamict Zircon under Different P-T-X
Conditions. Geochemistry International, 36, 11, 1006-1013.
Luck JM et Allegre CJ, 1982 The study of molybdenites through the 187Re-187Os
chronometer : Earth and Planetary Science Letters 31: 291-296.
187
Luck JM, Birck JL et Allegre CJ 1980. Re-187Os systematics in meteorites: Early
chronology of the solar system and age of the galaxy. Nature, 283, 256-259.
Ludwig KR 1999. Isoplot/Ex 2.0X. A Geochronological Toolkit for Microsoft Excel. Berkley
Geochronology Center Special Publication N°1a.
Maacha L, Azizi-Samir MR et Bouchta R 1998. Gisements cobaltifères du district de Bou
Azzer (Anti-Atlas) : structure, minéralogie et conditions de genèse. Chron. Rech. Min.,
531-532, 65-75.
Mandeville CW, Sasaki A, Saito G, Faure K, King R et Hauri E 1998. Open system
degassing of sulfur from Krakatau 1883 magma. Earth and Planetary Science Letters,
160, 709-722.
Manspeier W, Puffer JH et Cousminer HL 1978. Separation of Morocco and eastern North
America : A Triasic-Liasic Stratigraphic record. Geol. Soc. Am. Bull., 89, p. 901-920.
Marcantonio F, Reiseberg L, Zindler A, Wyman D et Hulbert L 1994. An isotopic study
of the Ni-Cu-PGE-rich Wellgreen intrusion of the Wrangellia Terrane: Evidence for
hydrothermal mobilization of rhenium and Osmium. Geochimica et Cosmochimica
Acta 58(2): 1007-1017.
Marcoux E et Jébrak M 1999. Butte, Montana, Etats-Unis: un exemple de téléscopage et de
district métallogénique stationnaire. Chronique de la Recherche Minière, 535: 45-58.
Marini F et Ouguir H 1990. Un nouveau jalon dans l'histoire de la distension pre-
panafricaine au Maroc: le Precambrien II des boutonnières du Jbel Sghro nord-orient
(Anti-Atlas, Maroc). Comptes rendus de l'Academie des siences Série 2, 310, 577-582.
183
Références bibliographiques
184
Références bibliographiques
Ohmoto H, Lasagna AC, 1982. Kinetics of reactions between aquous sulfates and sulfides in
hydrothermal systems : Geochimica et Cosmochimica Acta, v.46, p. 1727-1745.
Ouguir H 1991. Analyse sédimentologique, structurale et lithogéochimique de la série du
Protérozoïque supérieur (PII) de l’environnement de la mine d'argent d'Imiter (Anti-
Atlas Oriental, Maroc) : implications géodynamiques. Thèse I.N.P.L, Nancy 197p.
Ouguir H 1997. Contexte géologique du gisement argentifère d'Imiter (Anti-Atlas oriental,
Maroc). Contrôle volcanique et structural de la mise en place des concentrations
métalliques à Ag-Hg., Doctorat d’état es-Science, Faculté des Sciences de Meknès,
Maroc: 215p.
Ouguir H, Macaudiere J et Dagallier G 1996. Le Protérozoïque supérieur d’Imiter, Saghro
oriental, Maroc:un contexte géodynamic d’arrière-arc. J. of African Earth Sci. 22
(2):173-189.
Ouguir H, Macaudière J, Dagallier G, Qadrouci A et Leistel JM 1994. Cadre structural du
gîte Ag-Hg d'Imiter (Anti-Atlas, Maroc); Implication métallogénique. Bull. Soc. géol.
France., 165, pp. 233-248.
Paile Y 1983. Etudes des séries volcaniques de Précambrien III de l'Ougnat (Anti-Atlas
oriental, Maroc) et des minérailisations plombo-zincifères complexes associées (gïte
de Bou Madine). Thèse 3ème cycle, Paris XI Orsay, p.290.
Pasava J 1994. Geochemistry and rôle of anoxic sediments in the origin of the Imiter silver
deposit in Morocco. Czech Geol. Surv. Bull., 69, 1-11.
Petruck W, Harris DC et Stewart JM 1971.The silver arsenides deposits of the Cobalt-
Gondwana region, Ontario : Geochemistry of the Ores. Canadian Mineralogist, 11,
150-186.
Peucker-Ehrenbrink B et Hannigan G 2000. The marine osmium isotope record. Terra
Nova, 12, 205-219.
Pierson-Wickmann AC, Reisberg L et France-Lanord C, 1999. The Os isotopic
composition of Himalayan river bedloads aand bedrocks : Importance of Black shales.
Earth and Planetary Science Letters 176, 201-216.
Pimminger M, Grassenbauer M, Schroll M, Cerny I, 1984. Microanalysis in galena by
secondary mass spectrometry for determination of sulfur isotopes. Anal. Chem. 56,
407-411.
Piqué A, Bouabdelli M, Soulimaïni A, Youbi N et Iliani M, 1999. Les conglomérats du PIII
(Néoprotérozoïque supérieur) de l'Anti Atlas (Sud du Maroc) : molasses panafricaines,
ou marqueurs d'un rifting fini-protérozoïque. C. R. Acad. Sci. Paris, 328 (6), 409-414.
185
Références bibliographiques
Poller U 1997. U/Pb single zircon study of gabbroic and granitic rocks of Val Barlas-ch
(Silvretta nappe, Switzerland). Schweiz. Mineral. Petrogr. Mitt. 77: 351-359.
Popov AG 1995. Gisement argentifère d'Imiter : Etude minéralogique, paragenese et zonalité
du gisement. Rapport Interne SMI, p. 62.
Popov AG, Belkasmi A et Qadrouci A 1986. Le gisement argentifère d'Imiter synthèse
géologique, résultats des recherches de 1985 et bilan des réserves au 31/12/1985.
Rapport Interne SMI, p.168.
Pupin JP 1980. Zircon and Granite Petrology. Contributions to Mineralogy and Petrology
73: 207-220.
Pupin JP 1992. Les zircons des granites océaniques et continentaux : couplage typologie-
géochimie des éléments en traces. Bull. Soc. géol. France 163(4): 495-507.
Ramdohr P 1969. The Ore Minerals and their Intergrowths. Pergamon, New York
Ravizza G et Turekian KK 1992. The Osmium isotopic composition of organic-rich marine
sediments. Earth and Planetary Science Letters 110: 1-6.
Reisberg L, Zindler A, Marconio F, White W et Wyman D 1993. Os isotope systematics
in ocean island basalts. Earth and Planetary Science Letters 120: 149-167.
Rosman KJR et Taylor PDP 1998. Isotopic Compositions of the Elements 1997, Pure and
Applied Chemistry, 70, 217.
Roy Barman M 1993. Mesure du rapport 187Os/186Os dans les basaltes et les pérodotites :
Contribution à la systématique 187Re/187Os dans le manteau. PhD Thesis, Paris VII.
Roy-Barman M, Wasserburg G, Papanastassiou D et Chaussidon M 1998. Osmium
isotopic compositions and Re-Os concentrations in sulfide globules from basaltic
glasses.” Earth and Planetary Science Letters 154: 331-347.
40
Ruffet G, Féraud G, Amouric M 1991. Comparison of Ar/39Ar conventional and laser
dating of biotites from the North Trégor Batholith. Geochim. Cosmochim. Acta 55,
1675-1680.
Rye RO, Ohmoto H, 1974. Sulfur and Carbon isotopes and ore genesis – in review :
Economic Geology, v.69, p. 826-842.
Sakai H, 1968. Isotopic properties of sulfur compounds in hydrothermal processes :
Geochimical Journal (Japan), v.2, p. 29-49.
Sakai H, Des Marais DJ, Ueda A, Moore JG, 1984. Concentrations and isotope ratios of
carbon, nitrogen, and sulfur in ocean-floor basalts : Geochimica et Cosmochimica
Acta, v.48, p. 2433-2441
186
Références bibliographiques
187
Références bibliographiques
Shirey SB et Walker RJ, 1995 The Re-Os Isotope system in cosmochemistry and high-
temperature geochemistry. Rev. Earth Planet. Sci. 20: 423-497.
Shirey SB et Walker RJ, 1998. The Re-Os Isotope system in cosmochemistry and high-
temperature geochemistry. Rev. Earth Planet. Sci. 20: 423-497.
Singh S, Trivedi JR et Krishnaswami S 1999. Re/Os isotope systematics in balck shales
187
from the Lesser Himalaya: Their chronology and role in the Os/188Os evolution
seawater. Geochimica et Cosmochimica Acta 63(16): 2381-2392.
Sinha AK, Wayne DM et Hewitt DA 1992. The hydrothermal stability of
zircon :Preliminary experimental and isotopic studes. Geochim. Cosmochim. Acta. 56,
3551-3560.
Slodzian 1964. Etude d’une méthode d’analyse locale chimique et isotopique utilisant
l’émission ionique secondaire. Thèse d’Université, Orsay, p58.
Smoliar MI, Walker RJ et Morgan JW 1996 Re-Os ages of group IIA, IIIA, IVA and IVB
iron meteorites. Science 271, 1099-1102.
Speer JA 1980. Zircon In :Ortho-silicates, Ribbe P. (Ed) Mineralogical Society of America,
Washington.
St Gal de Pons A 1975. Le gisement filonien polymetallique de Bou Madine de l'Ougnat.
Rapport interne BRPM 431-15-GA-1-42.
Stacey JS et Kramers JD 1975. Approximation of terrestrial lead isotope evolution by a two-
stage model. Earth Planet. Sci. Lett., 26, pp. 207-221.
Steiger RH et Jäger E 1977. Subcommission on geochronology :Convention on the use of
decay constants in geo- and cosmochronology. Earth and Planetary Science Letters,
36, 359-362.
Stein HJ, John W, et Scherstin A 2000. Re/Os Dating of Low-Level Highly Radiogenic
(LLHR) Sulfides: The Harnas Gold Deposit, Southwest Sweden, Records Continental
Scale Tectonic Events. Economic geology, 95, 1657.
Stein HJ, Sunblad K, Morgan JW, Markey RJ et Motuza G 1998. Re-Os ages for
Archean molybdenite and pyrite, kuitilla, Finland and Proterozoic molybdenite,
Kabeliai, Lithania: a metamorphic and metasomatic test for the geochronometer :
Mineralium deposita, 33: 329-345.
Steinberg M, Treuil M, Touray JC, 1979. Geochimie principes et méthodes. Tome II :
Cristallochimie et éléments en traces. Doin (Eds.) Paris. 599p.
188
Références bibliographiques
Terra F et Wasserburg GF 1972. U-Th-Pb systematic in three Apollo 14 basalts and the
problem of the initial Pb in lunar rocks. Earth and Planetary Science Letters, 14, 281-
304.
Thieblement D et Tegyey M 1994. Une discrimination géochimique des roches différenciées
témoins de la diversité d'origine et de la situation tectonique des magmas calco-
alcalins. C.R. Acad. Paris, 319, série II, 87-94.
Thode HG, Dunford HB, Shima M, 1962. Sulfur isotope abondances in rocks of the
sudbury district and their geological signifiance. Economic Geology. 57, 565-578.
Thode HG, Monster J, Dunford HB, 1961. Sulfur isotope geochemistry. Geochimica et
Cosmochimica Acta. 25, 159-174.
Touil A, El-Boukhari A, Bilal E et Moutte J 1999. Les tholéiites a affinité alcaline du
secteur ouest du Siroua (Anti-Atlas central, Maroc): témoins d'une distension an
néoprotérozoïque. Journal of African Earth Sciences, vol. 29, 4, p. 699-713.
Tucker ME 1986. Carbon isotope excursions in Precambrian/Cambrian boundary beds,
Morocco. Nature 319, 48-50.
Turner G (1968) The distribution of potassium and argon in chondrites. In Origin and
distribution of the elements (ed. L.H. Ahrens), 387-398.Pergamon, London.
Vargas JM 1983. Etude métallographique des minéralisations mercuro-argentifère d'Imiter,
Fondation Scientifique de la géologie et de ses applications, Nancy, 60p.
Vavra G 1990. On the kinematics of zircon growth and its petrogenetic significance: a
cathodoluminescence study. Contrib. Mineral.Petrol. 106: 90-99.
Vavra G 1994. Systematics of internal zircon morphology in major Variscan granitoide
types. Contrib. Mineral. Petrol. 117: 331-344.
Velde B 1965. Phengite micas: Synthesis, stability and natural occurrences. Amer. J. Sci.,
n°263, 886-913.
Vialon P, Ruhland M et Grolier L 1976. Elements de tectonique analytique. Masson, 118 p.
Villemin T et Charlesworth H 1992. An interactive computer program to determine
paleostress axes using data from striated faults. Cordillera transect and cordilleran
tectonic workshop, Alberta.
Villeneuve M et Cornée JJ 1994. Evolution paleogeographique de la marge nord-ouest de
l'Afrique du Cambrien à la fin du Carbonifère (du Maroc au Liberia). Can. J. Earth
Sci., 28, 1121-1130.
189
Références bibliographiques
Walker RJ, Naldrett AJ, Li C et Fassett JD 1991. Re/Os isotope systematics of Ni-Cu
sulfide ores, Sudbury Igneous Complex, Ontario : evidence for a major crustal
component. Earth and Planetry Science Letters 105: 416-429.
Walker RJ, Shirey S, Hanson GN, Rajamani V et Horan MF 1989. Re/Os, Rb/Sr, and O
isotopic systematics of the Archean Kolar Schist belt, Karnataka, India. Geochimica et
Cosmochimica Acta 53: 3005-3013.
Walker RJ, Shirey SB, Horan MF et Echearria LM 1991. Re/Os isotopic contraints on the
origin of volcanic rocks, Gorgona Island, Colombia : Os isotopic evidence for ancient
heterogeneities in the mantle. Contribution to Mineralagy and Petrology 107: 150-
162.
Watson EB 1980. Some experimentally determined zircon/liquid partition coefficients for the
rare earth elements. Geochim. Cosmochim.Acta. 44, 895-897.
Watson EB et Cherniak DJ 1995. Oxygen diffusion in zircon. Earth and Planetary Science
Letters. 148, 527-544.
Watson EB et Cherniak DJ 1997. The incorporation of Pb into zircon. Chemical Geology
141: 19-31.
Watson EB et Harrison TM 1983. Zircon saturation revisited : temperature and composition
effects in a variety of crustal magma types. Earth and Planetary Science Letters. 64,
295-304.
Wayne M et Sinha K 1988. Physical and chemical response of zircons to deformation.
Contrib. Mineral. Petrol. 98: 109-121.
Westphal M, Montigny R, Thuizar R, Bardon C, Bossert A, Hamzeh R et Rolley JP
1979. Paléomagnétisme et datation du volcanisme permien, triasique et crétacé du
Maroc. Can. J. Earth Sci. 16, p. 2150-2164.
Wetherill GH 1956. An interpretation of the rhodesia and Witwatersand age patterns.
Geochim. Cosmochim.Acta, 9, 290-292.
White NC et Hedenquist JW 1995. Epithermal gold deposits: styles, characteristics and
exploration. SEC Newletters, Octobre 1995, 23, 1-13.
Wiedenbeck M, Allé P, Corfu F, Griffin WL, Meier M, Oberli F, Von Quadt A, Roddick
JC et Spiegel W 1995. Three natural zircon standards for U-Th-Pb, Lu-Hf, trace
elements and REE analysis. Geostandard Newsletters 19, 1-23.
Winchester JA et Floyd PA 1977. Geochemical discrimination of different magma series
and their differenciation products using immobile elements. Chem. Geol. 20, 325-343.
190
Références bibliographiques
Xiong Y et Wood SA 1999. Experimental determination of ReO2 and the dominant oxidation
state of rhenium in hydrothermal solutions. Chemical Geology 158: 245-256.
Zartman RE, Doe BR, 1981. Plumbotectonics : The model. Tectonophysics, 138, p. 135-
162.
Zotov AV, Kudrin Av, Levin KA, Shikina ND et Var'yash LN 1995. Experimental studies
of the solubility and complexing of the slected ore elements (Au, Ag, Cu, Mo, As, Sb,
Hg) in aqous solutions. In : "Fluids in the crust", Shmulovitch KI, et Yardley BWD
eds, p. 95-137.
Zouita F 1986. Etude de la distribution des terres rares et d’autres éléments en traces dans les
roches carbonatées du district minier de Bou Azzer-Bleida (Anti-Atlas, Maroc):
conséquences génétiques et variations à l’approche des minéralisations. Thèse
Univ.GrenobleI. p. 170.
191