Géologie de la région du lac de la Fourche, Province de Grenville, régions de la Mauricie et du Saguenay–Lac-Saint-Jean, Québec, Canada
Projet visant le feuillet 32A03
Francis Talla Takam, Abdelali Moukhsil
BG 2023-02
Publié le
Au cours de l’été 2022, un levé géologique à l’échelle 1/50 000 a été effectué dans la région du lac de la Fourche (feuillet SNRC 32A03), au NW de La Tuque. Ce levé couvre partiellement les régions du Saguenay–Lac-Saint-Jean et de la Haute-Mauricie.
Les affleurements de roche cartographiés appartiennent à l’Allochtone de la Province de Grenville. La Suite plutonique de Bardeau (1364 ±9 Ma), qui affleure à l’extrême est, constitue le socle cristallin de la région cartographiée. Elle est en contact cisaillé avec la Suite plutonique de Pope (1350 ±4 Ma), qui couvre la plus grande superficie, et sur laquelle repose l’unité supracrustale du Complexe de Wabash (<1204 ±4 Ma). Ces suites et ce complexe sont coupés par d’autres intrusions nouvellement définies : les Intrusions de Windigo, le Pluton de Wenache et le Batholite de Cristal. Ce dernier est polyphasé et à texture porphyroïde. D’autres intrusions, dont la Suite plutonique de Bonhomme (1030 ±19 Ma), la Suite intrusive de Rhéaume (1019 ±2 Ma) et la Suite de Roc (>1012 ±1 Ma) affleurent dans la région. Toutes les roches sont coupées par des dykes de pegmatite granitique à syénitique.
Le grain structural régional est principalement orienté NE-SW et a accommodé la mise en place des intrusions. À l’ouest, cependant, un couloir de déformation N-S de ~25 km de longueur sur 8 km de largeur est présent. Les roches de la région sont d’âge mésoprotérozoïque et ont connu plusieurs phases de déformation. Elles sont coupées par des failles et des zones de cisaillement. Le métamorphisme régional observé varie du faciès des granulites au faciès supérieur des amphibolites. Par endroits, un métamorphisme rétrograde est observé avec une ouralitisation de l’orthopyroxène, localement coronitique à bordure de clinopyroxène ± amphibole ± grenat. Ce phénomène se rencontre à l’ouest de la région à l’étude, plus spécifiquement dans le couloir de déformation. Un phénomène d’altération hydrothermale se superpose à ce dernier.
Plusieurs zones minéralisées favorables à l’exploration ont été mises au jour dans ce couloir d’altération hydrothermale et dans les roches du Batholite de Cristal :
- la zone du Lac du Relais, avec une minéralisation en ETR-Nb-Ta dans des dykes de pegmatite granitique;
- la zone du Lac Fogh, qui se caractérise par la présence de dykes décimétriques granitiques ou syénitiques à ETR ± Th;
- la zone du Lac Benès, avec une minéralisation en graphite ± zinc;
- la minéralisation de Fe-Ti ± P ± V dans les gabbronorites du Batholite de Cristal.
Méthode de travail
La région a été cartographiée en utilisant la méthode établie pour les levés effectués dans les zones forestières desservies par un réseau de chemins secondaires. Les travaux de cartographie géologique ont été réalisés par une équipe composée de deux géologues, d’un stagiaire en géologie et de deux étudiantes, du 3 juin au 19 août 2022. La cartographie et la synthèse du projet Lac de la Fourche ont permis de produire et de mettre à jour les éléments d’information présentés dans le tableau ci-contre.
Laboratoire de terrain
Le projet était muni d’un laboratoire de terrain simplifié. La susceptibilité magnétique est la seule propriété physique à avoir été mesurée selon le protocole établi par Christian Dupuis (Université Laval). Sous la supervision du géologue responsable, les étudiantes préalablement formées ont réalisé la prise de mesures sur les échantillons de roche en cassure fraîche et en face sciée en continu durant le programme de cartographie. Les mesures ont été réalisées sur la majorité des lithologies principales observées en affleurement et sur certaines lithologies secondaires jugées significatives, comme celles d’origine volcanique ou minéralisées.
| Élément | Nombre |
|---|---|
| Affleurement décrit (géofiche) | 546 affleurements |
| Analyse lithogéochimique totale | 96 échantillons |
| Analyse lithogéochimique des métaux d’intérêt économique | 27 échantillons |
| Analyse géochronologique | 5 échantillons |
| Lame mince standard | 95 |
| Lame mince polie | 26 |
| Fiche stratigraphique | 9 |
| Fiche structurale | 0 |
| Fiche de substances minérales métalliques | 2 |
| Fiche de substances minérales non métalliques | 1 |
| Mesure de susceptibilité magnétique | 322 |
Travaux antérieurs
Le tableau ci-dessous présente une liste des travaux réalisés dans le secteur à l’étude depuis 1965. Il inclut aussi les références citées dans le rapport. Une liste plus exhaustive peut être trouvée dans la base de données documentaire EXAMINE.
| Auteur(s) | Type de travaux | Contribution |
|---|---|---|
| Laurin et Sharma, 1975 | Cartographie et reconnaissance géologique à l’échelle 1/250 000 | Premiers travaux d’inventaire géologique couvrant la région à l’étude; projet de cartographie échelonné sur plusieurs années (1965 à 1967) |
| Benahmed et Intissar, 2015 | Levé magnétique aéroporté du secteur est du réservoir Gouin | Couverture géophysique à haute résolution |
| Cartographie géologique à l’échelle 1/50 000 | Géologie des feuillets 31P09 à 31P16, 32A01, 32A02, 32A07 et 32A08 |
Stratigraphie
Cette partie présente succinctement les différentes unités cartographiées de la région du lac de la Fourche dans un cadre stratigraphique et temporel. Le schéma stratigraphique et la frise chronologique ci-contre accompagnent cette discussion. La description lithologique de toutes les unités peut être consultée à partir de leur fiche stratigraphique respective via l’hyperlien associé à leur nom.
La chronologie relative de mise en place des unités est illustrée dans le schéma stratigraphique qui tient compte des relations de recoupement décrites sur le terrain. Ces observations sont supportées par quelques données géochronologiques U-Pb sur zircon. Ce schéma constitue une section très idéalisée de la géologie de la région étudiée montrant les relations de recoupement des unités ainsi que leurs positions relatives. La dimension des unités représentées sur le schéma respecte plus ou moins les superficies cartographiées en surface.
Comme mentionné précédemment, le secteur cartographié du lac de la Fourche se situe dans les régions de la Mauricie et du Saguenay–Lac-Saint-Jean, dans la partie allochtone de la Province de Grenville. Cette région regroupe des roches dont l’âge est mésoprotérozoïque (1364 à 1012 Ma).
Socle mésoprotérozoïque
Les suites plutoniques de Bardeau (mPbad, Moukhsil et Daoudene, 2019) et de Pope (mPpop, Moukhsil et Côté, 2017), lesquelles sont en contact cisaillé, constituent les plus anciennes unités stratigraphiques de la région cartographiée. La Suite plutonique de Bardeau, dont l’âge de mise en place est de 1364 ±9 Ma (Papapavlou, 2020), affleure essentiellement à l’extrême est du feuillet 32A03. La Suite plutonique de Pope, quant à elle, est constituée de deux sous-unités, mPpop1 et mPpop2. Elle est datée à 1350 ±4 Ma (David, 2018) et couvre la plus grande superficie. La sous-unité mPpop2 est prédominante. Ces sous-unités sont composées de roches intrusives intermédiaires à felsiques communément d’aspect gneissique. Elles sont constituées de charnockite, de granite, de mangérite, d’enderbite, de gneiss granitique et de monzodiorite (mPpop1). La sous-unité mPpop2 contient en plus de la monzonite. La Suite plutonique de Bardeau est composée de syénogranite, de granite à feldspath alcalin, de monzodiorite avec ou sans quartz, de jotunite, de syénite quartzifère, ainsi que de proportions mineures de granodiorite et de mangérite.
Unités Supracrustales
Le Complexe de Wabash (mPwab, âge maximum de dépôt de 1204 ±2 Ma, Moukhsil et al., 2015) est l’unité supracrustale qui repose sur le socle. Des quatre unités reconnues (mPwab1, mPwab2, mPwab3 et mPwab4) de ce complexe, seule la sous-unité mPwab2 est la plus abondante et affleure principalement à l’ouest du feuillet 32A03, dans le couloir de déformation de Windigo.
Intrusions mésoprotérozoïques
Les ensembles de roches intrusives et sédimentaires énumérés ci-dessus (mPbad, mPpop et mPwab) encaissent divers plutons, batholites et suites plutoniques d’envergure kilométrique constitués de roches felsiques, intermédiaires et mafiques.
Les nouvelles intrusions cartographiées sont les Intrusions de Windigo (mPwid), le Pluton de Wenache (mPWen) et le Batholite de Cristal (mPctl). Les Intrusions de Windigo sont un ensemble de charnockite, de granite, de granite d’anatexie et de granite à feldspath alcalin à hypersthène-grenat de teinte rouge ou rose lilas, ainsi que de syénite quartzifère à feldspath alcalin localisé dans le couloir de déformation de Windigo, au contact avec les roches de la Suite plutonique de Pope et du Complexe de Wabash (partie ouest du feuillet 32A03). Le Pluton de Wenache est constitué de granite à feldspath alcalin, de mangérite, de syénogranite et de granite avec une proportion mineure de charnockite. De forme sigmoïdale, le pluton affleure au centre sud de la région entre les roches de la Suite plutonique de Pope, de la Suite intrusive de Rhéaume et du Batholite de Cristal. Ce dernier est polyphasé; il affleure au centre nord de la carte et prend la forme d’une ellipse étirée NNE-SSW. Il est subdivisé en trois sous-unités, soit deux de composition intermédiaire à felsique (mPctl1, mPctl2) et une de composition est mafique à intermédiaire (mPctl3). La sous-unité mPctl1 est constituée de syénite quartzifère, de mangérite et de granite à feldspath alcalin porphyroïdes à phénocristaux de feldspath avec une proportion moindre de syénite et de syénite à hypersthène. La sous-unité mPctl2 se compose de mangérite, de syénite quartzifère et de syénite porphyroïdes à phénocristaux de feldspath avec une proportion moindre de granite à feldspath alcalin. La sous-unité mPctl3 regroupe plutôt une leuconorite, un gabbronorite à oxydes de Fe-Ti ± P et une proportion moindre de mangérite.
D’autres intrusions affleurent également dans la région cartographiée : la Suite Plutonique de Bonhomme (mPboh, Moukhsil et Daoudene, 2019; 1030 ±19 Ma, Papapavlou, 2019), au NE, la Suite intrusive de Rhéaume (mPrhe, Moukhsil et Côté, 2017; 1019 ±2 Ma, David, 2018), au SW, et les masses lenticulaires de la Suite de Roc (mPsro, >1012 ±1 Ma, Moukhsil et al., 2015) dans les roches intrusives tout comme dans les roches métasédimentaires.
Lithogéochimie
La lithogéochimie des unités géologiques de la région du lac de la Fourche est présentée séparément sous forme de tableaux.
Géologie structurale
Les roches de la région sont généralement déformées, surtout dans les zones de contact entre les intrusions et dans les couloirs et zones de cisaillement. Les plissements sont principalement observés dans les roches métasédimentaires du Complexe de Wabash, qui apparaissent sous forme de lentilles et de radeaux transposés à l’intérieur des roches plutoniques. Les plis y sont isoclinaux avec des plans axiaux subparallèles aux fabriques planaires des roches plutoniques. Les structures primaires de ces roches ont été transposées lors d’une première phase de déformation D1 et ont été oblitérées par les évènements de déformation des orogenèses subséquentes, dont l’Orogenèse grenvillienne. En effet, les roches de la région ont subi plusieurs phases de déformation et sont coupées par des failles majeures et des zones de décrochement.
Le grain structural régional est principalement orienté NE-SW et a accommodé la mise en place des intrusions. Cependant, à l’ouest de la région, il y a la présence d’un couloir de déformation d’orientation N-S de ~25 km de longueur sur 8 km de largeur. Sur la base de l’attitude des linéaments interprétés des cartes aéromagnétiques ainsi que de la distribution des structures planaires (D2/G2) et linéaires mesurées dans les roches, la région du lac de la Fourche (feuillet 32A03) a été subdivisée en quatre domaines structuraux, dont trois nouveaux : les domaines structuraux de Cristal (DSctl), de Wenache (DSwen) et de Windigo (DSwid); et le Domaine structural de Ministic (DSmin) défini par Moukhsil et Daoudene (2019) et Moukhsil et El Bourki (2019). Les limites entre ces domaines correspondent généralement au changement d’orientation de la fabrique planaire principale (Sn) jumelé à l’orientation des plans axiaux régionaux.
Le Domaine structural de Ministic (DSmin) a été défini dans les feuillets adjacents 32A02 et 32A07 (Moukhsil et Daoudene, 2019) et 32A10 (Moukhsil et El Bourki, 2019). Dans la région à l’édude, celui-ci occupe la moitié est du feuillet 32A03. Il montre une fabrique planaire variable avec une direction moyenne ENE-WSW et un pendage modéré à fort vers l’ESE. Les linéations minérales portées par les structures planaires se concentrent essentiellement dans le quadrant NE et sont à composante horizontale. Ce domaine se caractérise également par des plis déversés vers le NW et trois familles de failles; la première est d’orientation N-S et la seconde est orientée NW-SW avec un mouvement apparent senestre. La troisième famille est d’orientation ENE-WSW et est symétrique à une zone de cisaillement senestre majeure qui coupe les roches de la Suite plutonique de Pope.
Le Domaine structural de Cristal (DSctl) occupe le centre nord du feuillet 32A03. Il se caractérise principalement par la mise en place du Batholite de Cristal. La projection stéréographique des pôles des plans de la foliation (fabrique planaire) montre une disposition radiaire qui définit un pli conique ayant un plan axial subvertical orienté NNE-SSW avec un pendage vers l’ESE. Les linéations minérales mesurées ont une attitude très variable.
Le Domaine structural de Wenache (DSwen) se limite au Pluton de Wenache qui occupe le centre sud du feuillet 32A03. La projection stéréographique des pôles des plans de la foliation définit un pli vertical à axe subvertical dont le plan est orienté E-W. Ces structures planaires portent une linéation minérale qui est de plan axial.
Le Domaine structural de Windigo (DSwid) occupe la portion ouest du feuillet 32A03. Il comprend principalement la fabrique planaire des unités de la Suite plutonique de Pope, du Complexe de Wabash, des Intrusions de Windigo et de la Suite intrusive de Rhéaume. Ce domaine présente une fabrique planaire moyenne orientée N-S qui est généralement caractérisée par une foliation minérale secondaire (tectonométamorphique), une foliation ou un rubanement mylonitique et une gneissosité avec un fort pendage vers l’est. La projection stéréographique des pôles des plans de foliation montre une distribution radiale avec une forte concentration dans les côtés ouest et est, ce qui définit un pli droit d’axe N-S. Cette distribution fait aussi ressortir une interférence de pli E-W. Les linéations minérales portées par les structures planaires se disposent le long du plan axial N-S du plissement régional. Elles sont donc de plan axial. Ce domaine est caractérisé également par des failles à mouvement indéterminé de direction N-S à NW-SE. La Zone de cisaillement de Windigo à mouvement senestre qui coupe ce domaine du nord au sud le délimite en deux sous-domaines caractérisés à l’est par des plis déversés vers le NW, et à l’ouest essentiellement par des plis droits d’axe N-S.
Afin de montrer les relations entre les différentes unités lithostratigraphiques de la région du lac de la Fourche, une coupe structurale simplifiée dénotée AB dans le secteur du Batholite de Cristal a été réalisée selon une orientation E-W.
Métamorphisme
L’étude de l’évolution métamorphique de la région du lac de la Fourche a été faite à partir des observations minéralogiques et structurales à l’échelle macroscopique et microscopique (121 lames minces d’échantillons représentatifs des unités). La majorité des roches de la région montre des assemblages minéralogiques qui indiquent un métamorphisme régional prograde aux faciès des amphibolites et des granulites (carte métamorphique simplifiée).
Certains paragneiss du Complexe de Wabash contiennent une paragenèse à biotite (BO), sillimanite (SM) et grenat (GR) indiquant des conditions au faciès supérieur des amphibolites. Les roches de ce complexe sont migmatitisées à différents degrés de fusion partielle.
La présence d’orthopyroxène dans des roches intrusives felsiques à intermédiaires atteste de l’atteinte du faciès des granulites. La présence de structures perthitiques et mésoperthitiques (inclusions d’albite) dans le feldspath potassique contenu dans ces roches témoigne également des conditions de haute température.
Les roches de la région présentent également un métamorphisme rétrograde au faciès inférieur des amphibolites. Ceci est indiqué par un remplacement partiel ou total de l’orthopyroxène par l’amphibole (ouralitisation). Localement, celui-ci est coronitique à bordure de clinopyroxène ± amphibole ± grenat. Ce phénomène est observé à l’ouest de la région, dans le couloir de déformation de Windigo, où l’on observe un enrichissement potassique et ferrifère dans la roche et, localement, de la sillimanite dans le granite d’anatexie (lessivage de la roche). Ceci est probablement le résultat d’une altération hydrothermale. Au moins trois phases de métamorphisme ont été datées à 1090 ±41 Ma (Papapavlou, 2020), 1088 ±10 Ma (David, 2018) et 1012 ±1 Ma (Moukhsil et al., 2015) dans les unités intrusives cartographiées.
Géologie économique
La région du lac de la Fourche contient plusieurs zones minéralisées et zones favorables pour trois types de minéralisation :
- minéralisation en éléments de terres rares ± Th ± Nb ± Ta associée à des dykes de pegmatite granitique et à des roches intrusives syénitiques;
- minéralisation de graphite ± Zn associée aux roches métasédimentaires;
- minéralisation magmatique de Fe-Ti (± P ± V) associée aux roches intrusives mafiques.
Le tableau des zones minéralisées ci-dessous présente les résultats d’analyses pour deux nouvelles zones minéralisées découvertes dans le cadre de nos travaux.
| Nom | Teneurs | Commentaires |
|---|---|---|
| Minéralisation associée aux pegmatites granitiques | ||
| FT-1057 | 13 990 ppm ETR (G) | |
| FT-1065 | 3489 ppm ETR (G) | |
| Substance non métallique | ||
| René | Teneurs : L’analyse de l’échantillon 2022-FT-1081A a donné des valeurs de 2,9 % C graphitique, 8,84 % Al2O3 et 0,0205 % Zn. La zone minéralisée contient du grenat millimétrique de couleur rose lilas (jusqu’à 5 %) en association avec de la sillimanite (~1 %). | |
(G) : Échantillon choisi.
Le tableau des analyses lithogéochimiques des métaux d’intérêt économique donne la localisation, la description et les résultats d’analyse pour 27 échantillons choisis dans le but d’évaluer le potentiel économique de la région.
Minéralisations connues de la région d’étude
La région cartographiée ne présente aucune minéralisation connue.
Minéralisations méconnues et découvertes lors des présents travaux
Nouvelle zone minéralisée et nouvelle zone favorable en éléments de terres rares ± Nb ± Ta associées à des dykes de pegmatite granitique
À la suite des travaux de cartographie de l’été 2022, l’affleurement 22-FT-1057 a révélé des lectures de 1822 cps et de 264 ppm Th par spectrométrie de rayonnement gamma (TP_RS-230) sur un dyke métrique de pegmatite granitique coupant un granite à feldspath alcalin. La zone minéralisée FT-1057 est localisée à ~2,15 km au NE du lac du Relais. L’observation d’une lame mince en lumière polarisée et analysée montre que la roche est constituée de : quartz, feldspath potassique, allanite, titanite, zircon et minéraux opaques (pyrite). L’allanite, qui prend la forme de cristaux automorphes qui atteignent 2 mm de longueur, est associée à la titanite. Une cartographie au µXRF M4-Tornado d’une lame mince montre que les terres rares sont contenues dans l’allanite. L’analyse de l’échantillon 22-FT-1057-B1 a retourné des valeurs de 13990 ppm ETR (total) (7270 ppm Ce, 3810 ppm La, 1940 ppm Nd et 651 ppm Pr) dont 13849 ppm ETR légères.
En effet, une zone minéralisée en terres rares (ETR) a été mise en évidence sur la propriété du Lac du Relais (feuillet 32A04; Ouellet, 2015). Cette zone minéralisée en ETR lourdes, Nb et Ta est localisée dans des pegmatites granitiques en contact avec une faille majeure. L’échantillon composite 267263 (Ouellet, 2015) prélevé sur ce réseau de pegmatites a donné des valeurs de 0,33 % ETR lourdes, de 0,49 % Nb et de 229 ppm Ta. Nous proposons d’étendre cette zone vers l’est dans le feuillet 32A03 sous le nom de zone favorable du Lac du Relais.
Nouvelles zones minéralisée et favorable en éléments de terres rares ± Th associées à des dykes de syénite quartzifère à feldspath alcalin et de pegmatite granitique
La zone FT-1065, minéralisée en terres rares légères (ETR), est associée à un dyke décimétrique de syénite quartzifère à feldspath alcalin. Les phases minérales majeures sont le feldspath potassique, le plagioclase et le quartz. La biotite est mineure, chlorite et apatite sont des minéraux accessoires. Une cartographie au µXRF M4-Tornado d’une lame mince montre que les terres rares seraient probablement associées à l’apatite. L’analyse de l’échantillon 22-FT-1065-B1 a donné des valeurs de 795 ppm La, de 1650 ppm Ce, de 657 ppm Nd et de 185 ppm Pr pour un total de 3489 ppm ETR et de 402 ppm Th. Cet affleurement se situe dans la région de l’indice du Lac Fogh (>5 % ETR, Ouellet, 2015) dans le feuillet 33A04, à ~240 m à l’est de la zone du lac Fogh (quart SE du feuillet 32A03). Cette nouvelle analyse justifie de prolonger la zone favorable du Lac Fogh vers l’est (feuillet 32A03).
Nouvelles zones minéralisée et favorable pour le graphite et le zinc
La zone favorable du Lac Benès est associée aux roches métasédimentaires du Complexe de Wabash (mPwab), en contact avec les roches des Intrusions de Windigo. La roche hôte de la minéralisation est localement migmatitisée et peut contenir jusqu’à 10 % de graphite dans des niveaux rouillés. Le graphite, qui prend la forme de paillettes gris métallique de taille millimétrique, est disséminé ou en amas de grains submillimétriques formant des rubans millimétriques au sein de la roche. Ces zones minéralisées contiennent généralement du grenat millimétrique rose lilas (jusqu’à 5 %) en association avec de la sillimanite (~1 %). L’analyse de l’échantillon 22-FT-1081-A1 (Indice René) a donné des valeurs de 2,9 % C graphitique et de 205 ppm Zn.
Nouvelle zone favorable pour une minéralisation en Fe-Ti ± P ± V associée aux roches intrusives mafiques
Les unités mafiques du Batholite de Cristal constituent une zone favorable pour les minéralisations magmatiques de Fe-Ti-P-V. La zone favorable du Lac Cristal correspond à un ensemble de roches leuconoritiques et gabbronoritiques fortement magnétiques. Ces roches forment une bande circulaire associée à une anomalie magnétique positive sur la carte aéromagnétique (Benahmed et Intissar, 2015). L’analyse de l’échantillon 22-SC-2073-A1 a donné des valeurs de 20,5 % Fe2O3(T), de 4,18 % TiO2 et de 2,81 % P2O5.
La nouvelle carte géologique à l’échelle 1/50 000 améliore significativement les descriptions stratigraphiques, structurales et métallogéniques de la région du lac de la Fourche. Une étude approfondie du couloir de déformation de Windigo nouvellement identifié sera nécessaire afin de déterminer son lien structural, métamorphique et hydrothermal avec les roches environnantes.
| Auteurs | Francis Talla Takam, géo., Ph. D. francis.tallatakam@mrnf.gouv.qc.ca Abdelali Moukhsil, géo., Ph. D. abdelali.moukhsil@mrnf.gouv.qc.ca |
| Géochimie | Olivier Lamarche, géo., M. Sc. |
| Géophysique | Rachid Intissar, géo., M. Sc. |
| Évaluation de potentiel | Fabien Solgadi, géo., Ph. D. |
| Logistique | Marie Dussault, coordonnatrice |
| Géomatique |
Kathleen O’Brien Dominique Plante |
| Conformité du gabarit et du contenu | François Leclerc, géo., Ph. D. |
| Accompagnement /mentorat et lecture critique |
Abdelali Moukhsil, géo., Ph. D. Mhamed El Bourki géo. M. Sc. |
| Organisme | Direction générale de Géologie Québec, Ministère des Ressources naturelles et des Forêts, Gouvernement du Québec |
Remerciements :
Ce Bulletin GéologiQUE est le fruit de la collaboration de nombreuses personnes qui ont activement pris part aux différentes étapes de la réalisation du projet. Nous tenons à remercier le stagiaire en géologie Samuel Coulombe ainsi que les étudiantes Catherine Mai-Ly Daigneault et Alix Pailhon. Nous aimerions souligner l’excellent accueil de l’équipe du Relais 44, Louise Gosselin et Yvan Gauthier. Les échanges et discussions géologiques avec Dr. Louise Corriveau et Dr. Bertrand Rottier ont été très profitables.
Références
Publications du gouvernement du Québec
Autres publications
17 novembre 2023