La position du Burundi, dans le nord-est de la ceinture orogénique Kibarienne, lui confère un fort potentiel en ressources minérales. En effet, une variété de roches intrusives recoupe les roches métasédimentaires du Supergroupe Burundien, l’équivalent local de la chaîne Kibarienne. Ces roches intrusives sont composées de granites accompagnés de pegmatites minéralisées en Sn, Nb-Ta, W, Li, Be, de roches mafiques minéralisées en Fe-Ti-V, de carbonatites et syénites minéralisées en terres rares et en zircon. Outre les gîtes magmatiques, des gîtes hydrothermaux composés de filons de quartz et de brèches aurifères apparaissent dans des zones de cisaillement, à proximité des failles inverses. Des veines métasomatiques minéralisées en terres rares se retrouvent dans une zone de failles normales, situées dans le horst du rift est-africain. Suite aux conditions climatiques tropicales, des latérites nickélifères se sont développées sur des roches ultramafiques et constituent un gisement important de nickel latéritique (± Cr, Cu, EGP).

La formation de carbonates de magnésium et/ou calcium à l’aide du CO₂ atmosphérique (carbonatation minérale) permet la séquestration permanente du CO₂. Puisque la réaction se produit naturellement dans certains résidus miniers, il est possible d’utiliser ces derniers dans un esprit de développement durable. Le but de cette recherche est d’optimiser la carbonatation minérale de résidus miniers par des manipulations physico-chimiques et par l’ajout de sels organiques neutres. Cet ajout devrait permettre de faciliter la mise en solution des ions magnésium contenus dans le résidu minier et ainsi faciliter la précipitation de carbonates de magnésium. Des tests en eudiomètre ont été réalisés en variant différents paramètres tels que la concentration de CO₂ dans la phase gazeuse, la teneur en eau des résidus, le type de résidus utilisé et en ajoutant divers sels organiques tels que citrate, oxalate et EDTA de sodium en différentes concentrations afin de favoriser la carbonatation minérale. Les essais ont été réalisés sur différents résidus miniers riches en silicates/oxydes/hydroxydes de magnésium provenant des mines de Black Lake, Asbestos, Dumont, Raglan et Renard. Les résultats indiquent que le résidu plus réactif est de loin Dumont, suivi par Asbestos et Black Lake. Raglan et Renard sont très peu réactifs. La teneur en eau des résidus pour obtenir une réaction optimale varie de 20% à 60%. L’eau interstitielle est un solvant nécessaire pour la mise en solution des ions Mg, mais trop d’eau nuit à la diffusion du CO₂. Une relation linéaire entre la concentration en CO₂ dans la phase gazeuse et le taux d’absorption a été observée. En remplaçant l’eau par des solutions de chélates, l’absorption de CO₂ a pu être augmentée de 24% avec une solution d’EDTA d’un pH de 8,35 concentrée à 0,19M. Cependant, la meilleure solution varie en nature autant qu'en concentration d'un résidu minier à l'autre. Ces données montrent qu’il est possible de favoriser la séquestration naturelle de CO₂ de manière significative.

La région de l’Outaouais (13760 km²), au Québec, connait une forte croissance démographique et les eaux souterraines sont davantage sollicitées. La perméabilité du socle rocheux, essentiellement cristallin, est contrôlée par des fractures entrecroisées et donc hétérogène. Pour identifier des zones-cibles de forage, nous avons confronté des données de capacité spécifique du Système d’Information Hydrogéologique du Québec, validés à l’aide d’essais de pompage réalisés par l’Université Laval, avec des caractéristiques du roc, parmi lesquelles la lithologie, l’épaisseur de roc exploitée par les puits et la présence de linéaments. Il apparait que la texture et la minéralogie des roches ont une influence sur le développement de fractures perméables. On constate aussi que ces fractures se raréfient en profondeur et qu’il n’est donc pas toujours opportun de forer des puits à des profondeurs importantes (< 80 m). En revanche, aucune corrélation entre la présence de linéaments et le rendement des captage situés à proximité n’a pu être observée. / The Outaouais Region (13760 km²) is one of the most rapidly growing communities in Quebec, and where groundwater is increasingly being relied upon as a water supply. The permeability of the mostly crystalline bedrock is controlled by interconnected fractures and then heterogeneous. In order to identify target areas for installing new pumping wells, specific capacity data from the SIH database, which were validated using pumping tests carried out by Université Laval, were compared with rock characteristics including lithology, thickness of rock exploited by each well and proximity to lineaments. The texture and mineralogy of the rock seemed to have an important influence on the development of permeable fractures. We also found that the fracture density decreases with depth and therefore it is not always profitable to drill to deeper depths (> 80 m) to try to increase productivity. In contrast, proximity to lineaments didn’t seem to have any impact on the productivity of wells.

Une étude hydrogéologique conjointe de l'Université Laval et de la Commission Géologique du Canada (CGC) a été réalisée dans le bassin versant de la rivière Chaudière. Le but de cette étude est d’accroître les connaissances sur le système hydrogéologique du bassin versant de la rivière Chaudière qui s’étend de la frontière américaine jusqu’au fleuve Saint-Laurent. La zone étudiée a une superficie de 6675 km2. Il s’agit d’une région densément peuplée (~ 390,000 hab.) à vocation principalement agricole où 65 % de la population utilise l’eau souterraine comme principale source d’eau potable. L’aquifère régional se compose de deux unités principales : le substratum rocheux de la province géologique des Appalaches et certaines unités quaternaires granulaires. Le confinement de l’unité rocheuse est contrôlé par des sédiments glaciaires, marins ou lacustres fins d’épaisseurs et de distribution variables. L’écoulement dans le roc se fait dans sa partie supérieure où la fracturation est la plus importante. La première étape des travaux a d’abord consisté à définir les propriétés hydrauliques des formations géologiques. Les données utilisées proviennent de divers rapports, de données du système d’information hydrogéologique (SIH), du Ministère du Développement Durable, de l’Environnement, de la Faune et des Parcs du Québec (MDDEFP), de rapports de consultants ainsi que d’une campagne de travaux de terrain effectuée à l’été 2007. Parmi ces données, la conductivité hydraulique (K) constitue la principale propriété à définir. Plusieurs valeurs de K ont été extraites directement ou indirectement de rapports existants. Une partie importante des valeurs de la conductivité hydraulique a été estimée à partir de la base de données SIH. Parmi celles-ci, 1333 valeurs ont été évaluées à partir de la capacité spécifique (Bradbury et al, 1985); 89 de ces mesures ont étés qualifiées de données de qualités car elles respectent toutes les conditions de la méthode, à l’exception du temps de pompage souvent trop court (94% des cas). En se basant sur les mesures recueillies, la moyenne géométrique des données de conductivité hydraulique (K) s’élève à 5 x 10-6 m/s avec un écart-type du logarithme de K égal à 0.60. La lithologie ne semble pas avoir une grande influence sur les valeurs de la conductivité hydraulique car toutes les valeurs moyennes de chacune des lithologies se trouvent à l’intérieur d’un même ordre de grandeur. Un modèle numérique a été construit pour l’ensemble de la zone d’étude à l’aide du logiciel de modélisation par éléments finis Watflow (Molson et al., 2002) afin d’interpréter le système d’écoulement de l’eau souterraine au niveau régional. Le modèle tridimensionnel est composé de 28 couches dont 18 représentent les dépôts meubles et 10 le roc. Le maillage 3D comprend 1 / A joint hydrogeological assessment study between Université Laval and the Geological Survey of Canada (GSC) has been conducted in the Chaudière River watershed. The main objective of the study was to gain further understanding into the groundwater resources within the region. The watershed extends over ~ 6700 km2 from the United States border northwards to the Saint Lawrence River in the province of Québec (Canada). Groundwater is an important source of drinking water as well as for irrigation and industrial use, however the supply is finite and there are potential risks of conflicts arising from different users. It is thus important to understand the aquifer systems throughout the watershed in order to avoid possible conflicts and to help decision makers to better manage the resource. Within the catchment, up to 65% of the population relies on groundwater as the primary drinking water resource, whereby most of the water is used for agricultural purposes. The regional aquifer system is dominated by the underlying rock substratum of the Appalachian province and Quaternary aquifer units. The degree of confinement of the fractured aquifers is controlled by the overlying glacial, marine and fine lacustrine sediments of variable thicknesses. The groundwater circulation within the fractured aquifers predominantly takes place in the uppermost part of the fractured water-bearing units, where fracture density is highest. As a first step, the hydraulic properties of the fractured basement were defined, for which the hydraulic conductivity (K) of the formations can be considered as the most important hydraulic characteristic. A literature review of the available reports as well as the SIH and the MDDEFP data bases was performed in order to assemble the data. Additional field work was completed by the GSC during the spring of 2007 to supplement the data base. Most of the values of the hydraulic conductivity K were taken from the existing reports while others were estimated from the SIH data base. The data base contains 1333 values which were evaluated using Bradbury’s method, of which 83 satisfied the “quality criteria”, since they respected suggested conditions regarding the minimum pumping (which was not met in 94% of the cases). The geometric mean value of the measured hydraulic conductivities was 5x10-6 m/s with a standard deviation of 0.60. The maximum difference between the hydraulic conductivities of the different lithological units was about a factor of 100x. A 3D numerical finite element model was built with the Watflow model (Molson et al., 2002) for simulating the regional groundwater flow system. The model domain extends over the entire watershed and is subdivided into 28 layers, 18 for the unconsolidated deposits and 10 within the basement rock. The 3D triangular prismatic mesh contains 1,694,672 nodes and 3,344,516 elements. The topographic DEM (digital elevation model) forms the uppermost surface of the model. For the model boundary conditions, the Chaudière River and major tributaries and lakes are set as first-type boundaries with constant hydraulic heads. The model base, as well as the lateral limits, are considered as no-flow boundaries. The net recharge over the model surface was estimated using the infiltration model HELP, combined with estimates of the surface-distributed local water withdrawals. The major pumping wells are represented in the model domain as point sinks. The horizontal hydraulic conductivity together with the vertical anisotropy values for each formation were used as the model calibration parameters, using the observation data from 68 observation wells. The results of the calibration of the model showed that the mass balance error is close to zero which implies that the sum of inflows is equal to the sum of outflows. Different scenarios with recharge variations were simulated with the calibrated model. The simulations allowed identifying the principal sensitive zones which need piezometric and water quality control. It was also possible to determine the areas most sensitive to climate variations. Changes in the rate of pumping and recharge (± 20%) did not have a significant impact on the modeled system mainly in relation to groundwater level (observed variations were ≤5m).

Dans le cadre de la présente étude, un modèle géothermique tridimensionnel du pergélisol a été développé pour la toute première fois à l’échelle locale d’un village nordique afin d’évaluer les impacts des changements climatiques sur le pergélisol. Ce village est celui de Salluit au Québec nordique dans la zone de pergélisol continu. La majeure partie des infrastructures de ce village sont fondées sur un pergélisol salin riche en glace vulnérable aux changements climatiques. Le code numérique d’HydroGeoSphere a été modifié pour simuler la transmission de la chaleur par conduction avec changement de phase et tenir compte de la variation des propriétés thermiques du sol avec la température. Différents scénarios de changements climatiques lors du prochain siècle ont été utilisés pour alimenter le modèle géothermique. La variabilité spatiale du régime thermique du pergélisol simulée à l’échelle d’un village nordique est majeure lorsque l’hétérogénéité des dépôts, leur contenu en glace et les différentes conditions de surface observées sur le terrain sont considérés dans une modélisation numérique tridimensionnelle. L’influence des flux de chaleur latéraux dans cette modélisation produit des résultats qui sont dans certains cas significativement différents de ceux simulés au moyen d’un modèle unidimensionnel, particulièrement pour les sols sous d’importantes congères. Selon les scénarios de changements climatiques jugés pessimistes pour le maintien du pergélisol, un talik permanent se formera à partir de la seconde moitié du 21ième siècle sous les congères, en bordure d’un remblai routier où la neige et l’eau de fonte s’accumule à la surface d’un sol mal drainé et dans l’argile saline du secteur nord-est du village de Salluit. L’épaisseur maximale du talik à la fin du 21ième siècle varie entre 1,3 et 19,3 m respectivement dans l’argile et le roc sous les congères, alors qu’elle est de 2,8 m dans l’argile saline du secteur nord-est. Seuls les tills et le roc recouverts par d’importantes congères seront affectés par la présence d’un talik à la toute fin du 21ième siècle advenant un réchauffement climatique optimiste. Une dégradation plus importante du pergélisol sous les remblais routiers que pour un même sol en terrain naturel sous un couvert nival mince est à prévoir selon les scénarios climatiques pessimistes. Les infrastructures contruites sur des remblais seront potentiellement affectées par des tassements différentiels.