Lorsqu'ils sont exposés à l'eau de pluie et au dioxyde de carbone, les résidus miniers de l'exploitation du chrysotile capturent le CO2 grâce à une réaction de minéralisation du carbone. Ce phénomène a largement été étudié en laboratoire mais demeure mal compris à l'échelle d'une halde de stériles. Deux forages ont été réalisés au centre d'une halde à résidus miniers de Thetford Mines (Qc, Canada) dans le but de documenter la distribution spatiale de la réaction de carbonatation et d'identier sa signature pétrophysique. Les propriétés physiques des résidus miniers ont été évaluées à partir d'échantillons de la carotte de forage. An d'identier la signature pétrophysique des zones favorables à la carbonatation minérale, des diagraphies de géoradar, résistivité électrique et susceptibilité magnétique ont été réalisées. L'altération liée à la minéralisation du carbone n'est pas uniforme dans la halde, étant donné les variations de perméabilité à l'air et à l'eau. Des chemins d'écoulement préférentiels à l'air facilitent la carbonatation en profondeur en augmentant la quantité de CO2 qui entre dans la halde. Une hausse de résistivité électrique accompagnée d'une baisse de susceptibilité magnétique peuvent servir d'indicateur des zones aux conditions favorables pour la réaction de minéralisation du carbone. En effet, les carbonates qui sont résistifs précipitent entre les résidus de manière à les cimenter. La croissance des carbonates limite également la présence d'eau dans l'espace poreux ce qui contribue à l'augmentation de la résistivité électrique. La baisse de la susceptibilité magnétique, quant à elle, est attribuable à la destruction partielle de la magnétite pendant le processus de carbonatation et de météorisation. Le fer libéré lors de l'altération de la magnétite et de la serpentine est accumulé à la surface des grains et alimente la formation de goethite. Le degré d'altération de la magnétite est un indicateur du niveau de météorisation in-situ des résidus. Ces résultats sont la première étape pour suivre le progrès de la réaction de carbonatation dans le temps et pour quantier le volume de CO2 séquestré dans une pile de résidus miniers. / When exposed to rainwater and carbon dioxide, magnesium rich mining wastes spontaneuously capture CO2 through carbon mineralization. Extensively studied in the laboratory, this reaction is still ill-understood at the scale of industrial mining waste pile. In order to document the spatial distribution of carbon mineralization and identify its petrophysical signature, two holes were drilled at the center of a tailings pile in Thetford Mines, Qc, Canada. The physical properties of the mining residues were evaluated from drill core samples. The petrophysical signature of zones favorable for carbon mineralization was investigated using electrical resistivity, magnetic susceptibility and borehole radar data. These data showed that the variations of air and water permeability cause alterations that are related to carbon mineralization not being uniform throughout mining wastes. Preferential air-ow can enhance carbonation at depth by increasing the amount of CO2 that reaches the interior of the waste pile. The growth of carbonates which are resistive isolate the grains and restrict water ingress and contributes to the observed increases in electrical resistivity. The decrease in magnetic susceptibility during carbonation is attributed to the partial destruction of magnetite. The iron released during the alteration of magnetite and serpentine is accumulated on the grain surface and feeds the formation of goethite. The level of destruction of the magnetite provides insights into the level of alteration caused by in-situ weathering of the mining wastes. Therefore, an increase in electrical resistivity accompanied by a decrease in magnetic susceptibility can serve as an indicator for zones with favorable conditions for the carbon mineralization reactions. These results represent the rst step towards quantifying the volume of CO2 captured in these mining residues and track the progress of the reaction with time.

Au Nunavik, Québec, suite à la déglaciation et aux transgressions marines le long des côtes des Baies d’Hudson et d’Ungava, des sédiments fins gélifs se sont déposés dans ces secteurs côtiers. Lors de l’émersion postglaciaire suivant le relèvement isostatique, ces dépôts marins ont été mis en contact avec l’air froid et un pergélisol riche en glace s’est formé. Les infrastructures des communautés Inuites qui se trouvent aussi le long des côtes du Nunavik peuvent donc être construites sur un pergélisol riche en glace et leurs infrastructures sont peut-être vulnérables à la dégradation du pergélisol. Cette vulnérabilité peut affecter le développement durable de ces communautés au niveau de la performance, des coûts d’entretien et de la pérennité de leurs infrastructures nordiques. Afin de répondre aux besoins en ingénierie nordique d’investigation efficace des environnements pergélisolés pour soit éviter les zones vulnérables à la dégradation du pergélisol ou soit concevoir des méthodes de mitigation de cette dégradation, des objectifs ont été fixés dans le cadre de ce projet de recherche au doctorat: 1) de déterminer un modèle géocryologique conceptuel d’une butte de pergélisol riche en glace en zone de pergélisol discontinu au Nunavik à partir de l’interprétation d’une série de cinq essais de pénétration au cône réalisés dans cette butte, 2) en tenant compte de ce modèle conceptuel et à l’aide de la modélisation directe en géophysique, de développer une approche géophysique originale basée sur différentes méthodes géo-électriques pour l’auscultation du pergélisol, 3) de tester cette approche géophysique sur la butte de pergélisol riche en glace précédente, et 4) de développer des connaissances sur la géocryologie des buttes de pergélisol au Nunavik. Le site d’études est localisé dans une vallée profonde près de la communauté Inuite d’Umiujaq sur la côte est de la Baie d’Hudson. Des environnements pergélisolés caractéristiques du Nunavik dont notamment des buttes de pergélisol riche en glace sont facilement accessibles à partir de cette communauté. La méthodologie de recherche comprend l’utilisation de la modélisation directe à partir d’un modèle géocryologique de la butte de pergélisol étudiée pour concevoir l’investigation géophysique par tomographies de polarisation provoquée et diagraphies de résistivité électrique. Des tomographies bidimensionnelles et tridimensionnelles de polarisation provoquée ainsi que des diagraphies de résistivité électrique ont ensuite été réalisées pour caractériser le pergélisol du site d’étude. L’inversion et l’interprétation de ces levés géo-électriques ont permis de définir des modèles unidimensionnels, bidimensionnels et voire même tridimensionnels de résistivité et de chargeabilité électrique du site d’étude. Le modèle géocryologique initial de la butte de pergélisol étudiée a été amélioré en tenant compte des résultats de l’investigation géo-électrique. Le coeur riche en glace de la butte de pergélisol a été délimité avec précision à l’aide de cette approche d’investigation géo-électrique. En conclusion, l’approche d’investigation géo-électrique développée dans le cadre de ce projet de recherche au doctorat pour l’auscultation bidimensionnelle et tridimensionnelle du pergélisol peut être utilisée pour identifier des zones vulnérables à la dégradation du pergélisol avant la construction d’infrastructures nordiques ou au droit de telles infrastructures déjà construites pour évaluer leur vulnérabilité. / In Nunavik (Québec), frost-susceptible marine sediments were deposited along the coasts of Hudson Bay and Ungava Bay after their deglaciation and during the subsequent marine transgression. Following the postglacial emersion due to the isostatic uplift, the marine deposits came in contact with the cold air and ice-rich permafrost formed in these deposits. Therefore, the Inuit communities located along the coasts of Nunavik may be built on ice-rich permafrost and their infrastructures may be vulnerable to permafrost degradation. This vulnerability may impact the sustainable development of these communities affecting the performance, maintenance costs and service life of their infrastructures. In answer to the needs of efficient permafrost investigation in cold regions engineering to avoid zones vulnerable or use mitigation method to permafrost degradation, the main goals of this doctoral research project are: 1) to assess a conceptual geocryological model of an ice-rich permafrost mound in the discontinuous permafrost zone in Nunavik from the interpretation of five cone penetration tests carried out in the mound, 2) taking into account the previous model and using the forward modelling approach in geophysics, to develop an innovative geophysical approach based on different geo-electrical methods for permafrost investigation, 3) to test this geophysical approach for the characterization of the previous ice-rich permafrost mound, and 4) to contribute to knowledge development on the geocryology of permafrost mounds in Nunavik. The study site is located in a deep valley near the Inuit community of Umiujaq along the east coast of Hudson Bay. Permafrost environments characteristic of Nunavik such as ice-rich permafrost mounds are found in this valley which is accessible from this community. The research methodology consists in using the forward modelling approach in geophysics taking into account the conceptual geocryological model of the studied ice-rich permafrost mound to design the geophysical investigation based of induced polarization tomography and electrical resistivity logging. Two-dimensional and three-dimensional induced polarization tomography and electrical resistivity logging were then performed of the study site. Unidimensional, two-dimensional and three-dimensional models of electrical resistivity and chargeability of the study site were found from the inversion of these geo-electrical surveys. The conceptual geocryological model was then improved by interpretating these geo-electrical models. The ice-rich core of the permafrost mound was accurately delineated using this geophysical approach. In conclusion, the geo-electrical investigation approach developed herein for the two-dimensional and three-dimensional characterization of ice-rich permafrost environments can be used to delineate zones vulnerable to permafrost degradation before the construction of northern infrastructures or along existing infrastructures to assess their vulnerability.

De nombreux échantillons sont prélevés quotidiennement dans les mines d'or et envoyés aux laboratoires d’analyses chimiques pour en déterminer la teneur en or. Habituellement, un délai de 48 heures est nécessaire à l’obtention des résultats obligeant les opérateurs miniers à prendre des décisions opérationnelles avant de connaître les résultats. Pour répondre à cette problématique, la technologie émergente de spectroscopie d’émission de plasma induit par laser (S.L.P.) (dont l’acronyme anglais est Laser Induced Breakdown Spectroscopy LIBS) fournit une analyse quantitative de la concentration en or in situ et en temps réel. Un premier travail de cette technique a été effectué amenant la limite de détection de l’or à environ 0,7 ppm sur tous les types d’échantillons rencontrés en domaine minier (poudres, roches à tout venant, carottes…). Cette valeur correspond à la limite de détection de l’appareil souhaitée par les compagnies minières. En effet, les concentrations trouvées dans les gisements aurifères sont généralement supérieures à 0,7 ppm. Toutefois, certaines mines ont une teneur de coupure à cette valeur, ou dans cet ordre de grandeur, d’où la nécessité d’avoir une limite de détection plus basse. Ce projet de maîtrise a été conçu pour aider à l’étalonnage et à la validation de cette technique. Pour cela, des études pétrographiques et minéralogiques permettant de déterminer des lithologies vont être combinées à une analyse quantitative par fluorescence des rayons X (XRF). L’ensemble de ces résultats permettra de mieux comprendre le signal LIBS. Par la suite, ces résultats LIBS vont être validés en étant comparés avec ceux fournis par les laboratoires externes au projet. Ce mémoire de maîtrise a permis de démontrer que le LIBS est une méthode très prometteuse, avec un écart type relatif (<100%) en adéquation avec celui des laboratoires d’analyse chimique (<100%) pour une méthode plus rapide et non destructive. En effet, en comparant les résultats de teneur en or obtenus par analyse chimique conventionnelle provenant de deux laboratoires différents, les écarts types relatifs sur les échantillons ayant utilisé la courbe d’étalonnage riche en silice du LIBS, et cela quelle que soit la mine, sont inférieurs à 100%, mais deviennent supérieurs à 200% pour les échantillons ayant utilisé la courbe d’étalonnage riche en fer et soufre (mines Casa Berardi, Westwood et LaRonde). Les teneurs en or obtenues avec le LIBS, comparées avec les laboratoires externes, montrent que la technique LIBS est très prometteuse avec un écart type relatif également inférieure à 100 %, indépendamment de la courbe d’étalonnage utilisée et cela, quelle que soit la forme du matériel analysé. L’écart type relatif entre les laboratoires et le LIBS est bien souvent inférieur à celui trouvé entre les deux laboratoires utilisant le même protocole expérimental démontrant le caractère prometteur de la technique LIBS. Il est possible de noter une sous-évaluation des teneurs en or par rapport aux résultats obtenus avec les laboratoires sur les échantillons rocheux et les carottes de forage. Cette dernière s’explique par le caractère hétérogène de la distribution de l’or sur ces surfaces. Cependant, les analyses LIBS d’échantillons en poudre, ayant une meilleure homogénéité, ont un écart type relatif plus faible (<50%) et donc une meilleure précision sur ce type de matériel. Afin de pallier au problème de sous-évaluation sur les roches et carottes, il serait intéressant d’utiliser la méthode des k plus proches voisins permettant de tenir compte et d’étudier statistiquement les données censurées, c’est-à-dire les teneurs en or situées en dessous de la limite de détection de l’appareil LIBS. Actuellement, une généralisation de ces résultats n’est pas encore possible en raison du nombre d’échantillons analysés, mais ils restent très prometteurs pour la suite de ce projet. / Numerous rock samples are taken daily from gold mines and sent to laboratories to measure their gold content. A 48-hour delay in the delivery of results may force mining operators making operational decisions prior to knowing the results. To address this issue, the emerging Laser-Induced Breakdown Spectroscopy (LIBS) technology provides in-situ real-time quantitative gold analysis. First, this technique was performed by Rifai K. (NRCBoucherville postdoctoral fellow), bringing the gold detection limit to about 0.7 ppm on all types of samples encountered in mining operations (powders, rocks to all comers, cores ...). This value corresponds to the detection limit of the device desired by the mining companies. In fact, the concentrations found in the gold deposits are generally greater than 0.7 ppm. However, some mines have a cut-off grade at this value, or in this order of magnitude. This master’s project was set up to calibrate and validate this technique. For this, petrographic and mineralogical studies of representative lithologies will be combined with quantitative X-ray fluorescence (XRF) analysis. All of these results will provide a better understanding of the LIBS signal. Subsequently, these LIBS results will be validated by being compared with fire assay gold measurements provided by laboratories external to the project. This master’s thesis demonstrate that LIBS is a very promising method, with a percentage error (<100%) in line with that of laboratories for a faster and non-destructive method. In fact, by comparing the gold content results obtained by conventional chemical analysis from two different laboratories, the relative standard deviations on the samples that used the silica-rich calibration curve of LIBS, irrespective of the mine, are less than 100%, but becomes greater than 200% for samples using the calibration curve rich in iron and sulfur (Casa Berardi, Westwood and LaRonde mines). The gold contents obtained with LIBS, compared with the external laboratories show that the LIBS technique is very promising with relative standard deviation lower than 100%, regardless of the calibration curve used and the form of the material analyzed. The relative standard deviation between laboratories and LIBS is often lower than that found between the two laboratories. So, the LIBS technique is a very promising method. Similarly, it is possible to note an undervaluation of gold grades on rock samples and drill cores, due to the heterogeneous nature of the distribution of gold on these surfaces. However, the powder samples having a better homogeneity have a lower error percentage (<50%) and therefore a better accuracy on this type of material. In order to solve the problem of undervaluation on rocks and cores, it would be interesting to use the "k-nearest neighbors" method to take into account and statistically study censored data, i.e. those gold grades below the detection limit of the LIBS device. Currently, a generalization of these results is not yet possible due to the number of samples analyzed, but they remain very promising for the continuation of this project.

Lorsqu'ils sont exposés à l'eau de pluie et au dioxyde de carbone, les résidus miniers de l'exploitation du chrysotile capturent le CO2 grâce à une réaction de minéralisation du carbone. Ce phénomène a largement été étudié en laboratoire mais demeure mal compris à l'échelle d'une halde de stériles. Deux forages ont été réalisés au centre d'une halde à résidus miniers de Thetford Mines (Qc, Canada) dans le but de documenter la distribution spatiale de la réaction de carbonatation et d'identier sa signature pétrophysique. Les propriétés physiques des résidus miniers ont été évaluées à partir d'échantillons de la carotte de forage. An d'identier la signature pétrophysique des zones favorables à la carbonatation minérale, des diagraphies de géoradar, résistivité électrique et susceptibilité magnétique ont été réalisées. L'altération liée à la minéralisation du carbone n'est pas uniforme dans la halde, étant donné les variations de perméabilité à l'air et à l'eau. Des chemins d'écoulement préférentiels à l'air facilitent la carbonatation en profondeur en augmentant la quantité de CO2 qui entre dans la halde. Une hausse de résistivité électrique accompagnée d'une baisse de susceptibilité magnétique peuvent servir d'indicateur des zones aux conditions favorables pour la réaction de minéralisation du carbone. En effet, les carbonates qui sont résistifs précipitent entre les résidus de manière à les cimenter. La croissance des carbonates limite également la présence d'eau dans l'espace poreux ce qui contribue à l'augmentation de la résistivité électrique. La baisse de la susceptibilité magnétique, quant à elle, est attribuable à la destruction partielle de la magnétite pendant le processus de carbonatation et de météorisation. Le fer libéré lors de l'altération de la magnétite et de la serpentine est accumulé à la surface des grains et alimente la formation de goethite. Le degré d'altération de la magnétite est un indicateur du niveau de météorisation in-situ des résidus. Ces résultats sont la première étape pour suivre le progrès de la réaction de carbonatation dans le temps et pour quantier le volume de CO2 séquestré dans une pile de résidus miniers. / When exposed to rainwater and carbon dioxide, magnesium rich mining wastes spontaneuously capture CO2 through carbon mineralization. Extensively studied in the laboratory, this reaction is still ill-understood at the scale of industrial mining waste pile. In order to document the spatial distribution of carbon mineralization and identify its petrophysical signature, two holes were drilled at the center of a tailings pile in Thetford Mines, Qc, Canada. The physical properties of the mining residues were evaluated from drill core samples. The petrophysical signature of zones favorable for carbon mineralization was investigated using electrical resistivity, magnetic susceptibility and borehole radar data. These data showed that the variations of air and water permeability cause alterations that are related to carbon mineralization not being uniform throughout mining wastes. Preferential air-ow can enhance carbonation at depth by increasing the amount of CO2 that reaches the interior of the waste pile. The growth of carbonates which are resistive isolate the grains and restrict water ingress and contributes to the observed increases in electrical resistivity. The decrease in magnetic susceptibility during carbonation is attributed to the partial destruction of magnetite. The iron released during the alteration of magnetite and serpentine is accumulated on the grain surface and feeds the formation of goethite. The level of destruction of the magnetite provides insights into the level of alteration caused by in-situ weathering of the mining wastes. Therefore, an increase in electrical resistivity accompanied by a decrease in magnetic susceptibility can serve as an indicator for zones with favorable conditions for the carbon mineralization reactions. These results represent the rst step towards quantifying the volume of CO2 captured in these mining residues and track the progress of the reaction with time.

Un modèle numérique bidimensionnel a été développé afin d’évaluer l’impact de l’écoulement d’eau souterraine sur la dynamique du pergélisol dans un contexte de réchauffement climatique au Québec nordique. Le modèle conceptuel développé concerne une butte de pergélisol située dans la zone de pergélisol discontinu à proximité de la communauté Inuite d’Umiujaq, Nunavik, Québec. Le pergélisol s’est mis en place dans une unité gélive de silts marins qui se trouve audessus de deux unités de sédiments grossiers de sable et de gravier d’origine fluvio-glaciaire et glaciaire qui forment un aquifère confiné par l’unité de silts et le pergélisol où il y a un écoulement d’eau souterraine. Le code numérique HEATFLOW a été utilisé pour simuler l’écoulement d’eau souterraine couplé à la transmission de chaleur par conduction et advection le long d’une coupe 2D orientée dans la direction de l’écoulement de l’eau souterraine au droit de la butte de pergélisol étudiée. En premier lieu, le modèle a été étalonné manuellement à partir de profils de température mesurés dans la butte au cours des 10 dernières années à l’aide de câbles à thermistances et en tenant compte des flux de chaleur mesurés près de la surface du sol. En second lieu, une deuxième simulation a été réalisée en ne considérant que la transmission de chaleur par conduction et en négligeant ainsi l’écoulement d’eau souterraine. La comparaison entre ces deux simulations révèle le rôle important de l’écoulement d’eau souterraine sur la dynamique du pergélisol à Umiujaq. En effet, cet écoulement transporte l’eau plus chaude des zones de recharge vers l’aquifère confiné, ce qui contribue à réchauffer significativement le système en comparaison avec le cas sans écoulement. Une couche de pergélisol beaucoup plus mince est simulée lorsque l’écoulement d’eau souterraine est considéré dans la modélisation numérique. En outre, selon les résultats des simulations, l’énergie se dissipe le long de la ligne d’écoulement d’eau souterraine sous la base du pergélisol ce qui réduit sensiblement les températures du sol et de surface à proximité des zones de résurgence de l’eau souterraine le long d’un ruisseau en comparaison avec les zones de recharge. Finalement, en troisième lieu, le comportement futur du système simulé sous l’effet des changements climatiques est ensuite prédit en générant un scénario de réchauffement climatique selon une augmentation constante de la température de l’air et des précipitations. Les résultats des simulations suggèrent une dégradation du pergélisol par la base à un taux de 80 cm par année, et par le toit à un taux de 12 cm par année, jusqu’à la disparition complète du pergélisol dans le site d’étude d’ici 2040. / A 2D numerical model has been developed to assess the impacts of groundwater flow on permafrost dynamics under a warming climate in northern Québec. The conceptual model developed herein focuses on a small permafrost mound located in the discontinuous permafrost zone near the Inuit community of Umiujaq, Nunavik, Québec. At the study site, permafrost is found in marine silt overlying a deep confined sand and gravel aquifer with active groundwater flow. To better understand the cryo-hydrogeological system, the HEATFLOW numerical code was used to simulate coupled groundwater flow and heat transport by conduction and advection along a 2D cross-section through the permafrost mound and oriented along the assumed direction of groundwater flow. The model was first calibrated manually using temperature profiles in the permafrost mound measured along thermistor cables over the past 10 years and using observed heat fluxes near the ground surface. A second simulation was then performed assuming only conductive heat transfer and neglecting groundwater flow. A comparison between both simulations reveals the important role of groundwater flow on permafrost dynamics at the Umiujaq site. As groundwater flow brings warmer water from recharge areas into the deep confined aquifer, it contributes significantly to warming of the system relative to conduction alone, and significantly decreases permafrost thickness. However, the simulation showed that thermal energy is also lost along the flowpath below the permafrost base which induces a cooling in the discharge areas in comparison to the recharge areas. The future system behavior is then predicted by taking into account a climate change scenario based on increases in temperature and precipitation. The predictive simulation suggests that permafrost will thaw from the base at a rate of about 80 cm per year, and from the permafrost table at a rate of 12 cm per year, until completely thawed by about 2040.

La Sous-province d’Opinaca, située dans le craton du Supérieur, est un important bassin néoarchéen de sédiments métamorphisés. Les modèles tectoniques disponibles pour la formation puis l’évolution de l’Opinaca sont variés, allant du prisme accrétionnel au bassin d’arrière-arc. L’étude tectonométamorphique proposée permet de définir le parcours prograde de la région et ainsi discriminer les modèles probables, chacun impliquant une progression du métamorphisme distincte. Dans le nord-ouest de la Sous-province d’Opinaca, des isogrades tracés à partir de minéraux indicateurs du métamorphisme soulignent une progression du faciès des schistes verts dans le nord-ouest au faciès supérieur des amphibolites au sud-est. L’étude de la chimie minérale des minéraux métamorphiques révèle une homogénéisation généralisée des zonations chimiques. Les relations texturales observées dans le secteur au nord de l’Opinaca et une modélisation de l’équilibre des phases indiquent un chemin P-T-t en sens horaire avec un segment en décompresson de 6 à 4 kbar, à ~ 600 °C. Dans la partie sud du secteur d’étude, le chemin PTt est caractérisé par une décompression isothermique suparsolidus de 9 à 5 kbar, à ~ 800 °C, L’âge de la sédimentation est définie par les populations de zircons détritiques et le recoupement d’intrusifs felsiques entre 2712 Ma et >2690 Ma selon les secteurs. La géochronologie U-Pb sur des monazites indiquent deux pics de métamorphisme, à ~ 2670 et ~ 2645 Ma. La datation des grenats avec le système Lu-Hf indique que la croissance du grenat correspond à la plus jeune génération de monazite, circa 2645 Ma. La première génération de monazite est interprétée comme liée à un événement métamorphique de basse pression n’impliquant pas la croissance de grenat. Ces résultats permettent de proposer pour l’évolution du bassin d’Opinaca un polymétamorphisme dans un environnement géodynamique d’inversion tectonique d’un basin. Après la sédimentation, un premier épisode métamorphisme de haute température affecte le bassin en extension vers 2670 Ma. Vers 2645 Ma, la fermeture du bassin et l’épaississement crustal provoque un deuxième épisode de métamorphisme, de type barrovien, caractérisé par une décompression isothermale possiblement aidée par un certain degré d’extrusion ductile / The Opinaca Subprovince, in the Superior craton, is an important Neoarchean metamorphosed sediments basin. Various models are proposed to explain the formation and evolution of the Opinaca, ranging from accretional prism to back-arc basin. The proposed tectonometamorphic study allows the determination of the prograde evolution of the basin and hence the selection of the most probable regional model, each being characterized by contrasting styles of metamorphic progression. In the north-western Opinaca Subprovince, isograds traced from index minerals highlight a progression from greenschist facies in the north-west to upper amphibolites facies in the south-east. The chemistry of metamorphic minerals shows a global homogenisation of growth zonation. In the northernmost Opinaca Subprovince, textural relations and phase equilibrium modelling indicate a clockwise PTt path with an isothermal decompression segment from 6 to 4 kbar at ~ 600 °C. In the southern part of the study region, the clockwise PTt path is characterized by stronger, suprasolidus isothermal decompression from 9 to 5 kbar at ~ 800 °C. We constrain deposition of the Opinaca greywacke from 2712 to 2690 Ma with the youngest detrital zircon population and crosscutting felsic intrusions. U-Pb monazite geochronology indicates two pulses of metamorphism, at ~2670 and ~ 2645 Ma. Lu-Hf dating of garnet supported by textural analysis and trace element mineral chemistry indicates that garnet growth is coeval with the younger population of monazite, circa 2645 Ma. The first generation of monazite around 2670 Ma is thus interpreted as a low-pressure metamorphic event that did not involve garnet growth. These results point towards a polymetamorphic evolution for the Opinaca Subprovince consistent with the tectonic inversion of a rift-like basin. Clastic sedimentation between 2712 to 2690 Ma in a magmatically active, rift-like basin was followed by regional low-P metamorphism at 2670 Ma. The onset of crustal shortening and thickening in the basin by 2645 Ma resulted in Barrovian-type metamorphism, and involved isothermal decompression that could have been accommodated by some degree of ductile extrusion.

L'implication des émissions de dioxyde de carbone (CO2) anthropiques dans les changements climatiques est aujourd'hui admise et des solutions émergent pour lutter contre l'accumulation de CO2. La minéralisation du carbone, qui permet de séquestrer le CO2 sous forme de carbonates, stables à l'échelle géologique, est une des options envisagées. Parmi les voies de minéralisation du carbone envisagées, la minéralisation passive des résidus miniers ultramafiques permettrait de compenser les émissions en CO2 d'une exploitation minière. Toutefois, les impacts sur la qualité des eaux de lixiviation et l'évolution de la capacité de séquestration en conditions naturelles, à moyenne et grande échelle, sont peu documentés. La compagnie RNC Minerals a pour objectif d'exploiter un gisement de nickel situé dans le Nord-Ouest de la province du Québec. L'exploitation du Projet Dumont Nickel (PDN) aboutirait à la production d'environ 1,7 Gt de résidus miniers ultramafiques. Les différents facteurs qui influencent la capacité de séquestration des résidus du PDN ont été étudiés en laboratoire, à des teneurs en CO2 variables. Dans cette étude, les processus de la minéralisation passive dans les résidus du PDN, sont décrits et la capacité de séquestration en CO2 atmosphérique est estimée à moyenne échelle, en conditions naturelles. Pour étudier les impacts de l'altération météorique des résidus miniers du PDN, deux cellules expérimentales ont été construites et instrumentées. La première EC-1, contient les résidus ultramafiques, qualifiés de stériles (Waste-rock) et la seconde EC-2 a été remplie avec les résidus d'usinage (Tailings). Les propriétés hydrogéologiques et la surface spécifique des résidus des deux cellules sont différentes alors que la minéralogie est similaire. Les résidus sont composés principalement d'antigorite, de lizardite, de chrysotile, de brucite, de magnetite et de chlorite. Entre 2011 et 2015, l'évolution de la concentration en CO2, de la minéralogie, et de la composition chimique des lixiviats ont été enregistrées. Le suivi des concentrations en CO2 permet d'observer une diminution de la concentration en CO2 de la surface (~390 ppmv) vers le fond des parcelles (~100 ppmv). Dans le même temps, la teneur en carbone dans les résidus altérés a augmenté et les analyses minéralogiques révèlent la présence de plusieurs carbonates de magnésium comme l'hydromagnésite. Ces données suggèrent que les résidus séquestrent du CO2 passivement. Dans les cellules expérimentales le CO2 peut provenir de 3 sources : (1) l'atmosphère, (2) la dégradation de la matière organique, et (3) la dissolution des carbonates. Les compositions isotopiques du CO2(g), et des carbonates néoformés ont été mesurées. Ces analyses ont permis de mettre en évidence que la dissolution du CO2(g) dans l'eau interstitielle limite la capacité de séquestration et que le CO2 atmosphérique est la source du CO2 séquestré. Malgré les différences entre les deux cellules expérimentales les même processus contrôlent la séquestration du CO2. Un modèle conceptuel de la réaction de minéralisation du carbone, comprenant l'évolution de la composition isotopique, est proposé. Les lixiviats, récoltés aux bas des cellules expérimentales entre mai et novembre depuis 2011 sont caractérisés par un pH alcalin (~9,5), une alcalinité élevée (~90 à ~750 mg/L) et une forte concentration en magnésium (~50 à ~750 mg/L). Cette composition est en accord avec l'altération des résidus ultramafiques en milieu ouvert au CO2. Depuis 2012, la composition chimique des lixiviats évolue en fonction des saisons. Ces variations saisonnières sont expliquées par : (1) les variations climatiques au cours d'une année et (2) l'augmentation de la précipitation de carbonate entre mai et juillet. La diminution saisonnière de l'alcalinité et de la concentration en magnésium, provoqué par l'augmentation de la précipitation de carbonates, induit une sous-saturation des minéraux carbonatés ce qui limite la capacité de séquestration en CO2. Un taux de séquestration en CO2 atmosphérique de 1,4 (+/- 0.3) kg CO2/tonne/an a été mesuré dans les résidus de concentrateur (EC-2). À l'échelle de l'exploitation minière, les résidus de concentrateur permettraient la séquestration de 21 kt de CO2 atmosphérique par an ce qui correspond à un quart des émissions annuelles de la future mine. Le modèle MIN3P, qui permet de simuler le transport réactif multi composants et multiphasiques dans un milieu poreux insaturé, a été utilisé pour simuler en 1D la réaction de minéralisation au centre de la cellule EC-2. L'ensemble des données récoltées a été utilisé pour calibrer le modèle. Toutefois, aucune des simulations n’a permis de reproduire l'évolution de la géochimie des lixiviats et la concentration en CO2 observés. Plusieurs simplifications du modèle conceptuel pourraient expliquer les différences avec les données observées. / The implication of anthropogenic carbon dioxide (CO2) emissions in climate change is now widely accepted and solutions are emerging in order to limit the accumulation of CO2. Carbon mineralization, which allows the sequestration of CO2 through carbonate precipitation, stable minerals over geological time scales, is one of the options considered. Among the proposed carbon mineralization pathways, passive carbon mineralization in ultramafic mining residues can potentially lead to developing carbon-neutral mines. However, the impacts on leachate water quality and evolution of sequestration capacity in natural conditions, on medium and large scales, are still poorly documented. RNC Minerals plans to mine a nickel deposit located in the northwestern part of Quebec. The operation at the Dumont Nickel Project (DNP) would produce approximately 1.7 Gt of ultramafic mining residues. Several factors which influence the carbon sequestration capacity of the DNP residues have been studied in the laboratory, at variable CO2 concentrations. In this study, the processes of passive carbon mineralization in the DNP mining residues are described and the atmospheric CO2 sequestration capacity is estimated, at the experimental cell scale, under natural conditions. In order to study the impacts of meteoric weathering of the DNP residues, two experimental cells were built and instrumented. The first cell EC-1, contains the ultramafic waste rock, and the second EC-2, was filled with milling residues (Tailings). The hydrogeological properties and surface area of the residues contained in the two cells are different whereas the mineralogy is similar. The main minerals in the residues are chrysotile, lizardite, brucite, chlorite and magnetite. Between 2011 and 2015, changes in CO2 concentrations, mineralogy, and chemical composition of leachate waters were recorded. Monitoring of CO2 concentrations showed a decrease in CO2 concentration from the surface (~ 390 ppmv) to the bottom of the cells (~ 100 ppmv). At the same time, the carbon content in the weathered residues increased and the mineralogical analyses revealed precipitation of several magnesium carbonates such as hydromagnesite. These observations indicate that passive mineral carbonation of the mining residues is occurring within the experimental cells, for which three potential sources of CO2 can be identified : (1) the atmosphere, (2) the CO2(g) produced from organic matter oxydation, and (3) CO2(g) produced from carbonate dissolution. The isotopic compositions of CO2(g) and newly formed carbonates were measured. Using these isotopic compositions it was possible to demonstrate that dissolution of CO2(g) in interstitial water limits the sequestration capacity and that atmospheric CO2 is the main source of the CO2 sequestered. Despite the differences between the two experimental cells the same processes control CO2 sequestration. A conceptual model of the carbon mineralization reactions, including evolution of the isotopic compositions, is proposed. The leachate water sampled at the bottom of the experimental cells, between May and November since 2011, is characterized by an alkaline pH (~9.5), a high alkalinity (~90 to ~750 mg/L CaCO3) and a high concentration of magnesium (~50 at ~750 mg/L). This composition is consistent with weathering of ultramafic rocks in a system open to CO2. Since 2012, the chemical composition of the leachate water was evolved seasonnaly. These seasonal variations are explained by: (1) recharge and temeprature variations over the year and (2) increased carbonate precipitation between May and July. The seasonal decrease of alkalinity and magnesium concentrations, caused by increased carbonate precipitation, induces undersaturation of carbonate minerals. Therefore carbonate precipitation self-limits carbon sequestration through a negative feed-back loop. Since 2011, an estimated 13 kg of atmospheric CO2 was sequestered in the milling residues from EC-2, which corresponds to a mean rate of 1,4 (+/- 0.3) kgCO2/tonne/year. Using this mean rate, during the mining operation the milling residues will sequester about 21 kt of atmospheric CO2 each year, which will represents one quarter of the 127,700 tonnes of CO2 emitted. Using MIN3P, a numerical model which allow to simulate multi-component and multiphase reactive transport in unsaturated porous media, the carbon mineralization reactions were simulated in 1D at the center of cell EC-2. The data collected during the 4 years of monitoring were used to calibrate the numerical model. However, none of the simulations allowed to reproduce the evolution of the leachate water geochemistry and the CO2 concentrations observed in the experimental cell. Several simplifications of the conceptual model could explain the differences with the observed data.