Numerical simulation of groundwater flow in regional rock aquifers, southwestern Quebec, Canada

The St. Lawrence Lowlands platform, Quebec, Canada, is a densely-populated area, heavily dependent on groundwater resources. In 1999, the Geological Survey of Canada initiated a large-scale hydrogeological assessment study over a 1,500 km² region northwest of Montreal. The objectives were to define the regional groundwater flow, and to give quantitative estimates of the groundwater dynamic parameters and of the available groundwater resources. The applied approach consisted of defining the hydrogeologic framework, hydraulic properties of the aquifer units, and groundwater dynamic components. Lower Paleozoic sedimentary rocks represent regional aquifer units. Coarse Quaternary fluvio-glacial sediments locally overlay the rock sequence and constitute an interface aquifer unit. Fine marine sediments confine most of the regional aquifers. Collected GIS based information was synthesized in a finite element numerical model. The regional saturated steady-state flow was calibrated under current stress conditions assuming an equivalent porous medium approach. Water budget calculations show that the total groundwater flow in regional aquifers amounts to 97.7 Mm³/y. Infiltration from precipitation provides 86.6% of the groundwater supply, while 9.6% comes from subsurface inflow and the remaining 3.8% is induced recharge from surface waters. Discharge from regional aquifers occurs through flow to streams (76.9%), groundwater withdrawal (18.4%), and underground outflow (4.7%).

---

Resumen La plataforma de Tierras Bajas San Lorenzo, Quebec, Canadá, es un área densamente poblada que depende fuertemente de recursos de agua subterránea. En 1999 el Servicio Geológico de Canadá inició un estudio de evaluación hidrogeológica a gran escala sobre un área de 1,500 km² en la región noroeste de Montreal. Los objetivos fueron definir el flujo regional de agua subterránea y aportar estimados cuantitativos de los parámetros dinámicos de agua subterráne y de los recursos disponibles de agua subterránea. El enfoque aplicado consistió en definir el marco hidrogeológico, propiedades hidráulicas de las unidades acuíferas, y los componentes dinámicos de agua subterránea. Rocas sedimentarias del Paleozoico Inferior representan unidades regionales de acuíferos. Sedimentos fluvio-glaciares Cuaternarios gruesos sobreyacen localmente la secuencia rocosa y constituyen una unidad acuífera interfacial. Sedimentos marinos finos confinan la mayoría de acuíferos regionales. Información colectada basada en SIG se sintetizó en un modelo numérico de elemento finito. El flujo regional saturado en régimen permanente se calibró bajo condiciones de stress asumiendo un enfoque de medio poroso equivalente. Los cálculos de balance hídrico muestran que el flujo total de agua subterránea en acuíferos regional alcanza 97.7 Mm³/año. Infiltración a partir de lluvia aporta el 86.6% del abastecimiento al agua subterránea, mientras que el 9.6% proviene de entradas subsuperficiales y el restante 3.8% consiste de recarga inducida a partir de aguas superficiales. La descarga proveniente de acuíferos regionales ocurre a través de flujo a ríos (76.9%), utilización de agua subterránea (18.4%), y salida subterránea (4.7%).

---

Résumé La plateforme Lowllands du Saint-Laurent, Québec, Canada, est une aire densément peuplée, dépendant grandement des ressources en eau souterraine. En 1999, le Service Géologique du Canada a initié une étude hydrogéologique à grande-échelle sur 1500 km² au Nord-Ouest de Montréal. Le sobjectifs ont été de définir la dynamique de lécoulement régional, et de donner des estimations quantitatives des paramètres dynamiques des ressources disponibles en eau souterraine. Lapproche utilisée consista à définir le cadre hydrogéologique de travail, les propriétés hydrauliques des unités aquifères, et les composantes dynamiques des eaux souterraines. Les roches sédimentaires du Paléozoïque Inférieur représentent les unités aquifères régionales. Les sédiments marins fins confinent la plus grande partie des aquifères régionaux. Les informations de base, collectées dans un SIG, ont été synthétisées dans un modèle numérique aux éléments finis. Lécoulement permanent régional, en zone saturée, a été calibré en conditions de stress en assumant une approche de milieu poreu équivalent. Les calculs du bilan hydrologique montrent que lécoulement total des eaux souterraines dans les aquifères régionaux atteind 97,7 Mm³/an. Linfiltration à partir des précipitations apporte 86.6% de leau souterraine exploitée; sachant que 9,6% proviennent découlement de subsurfaces, et que les 3,8% restants proviennent de la recharge via les eaux de surface. Le débit pompé des aquifères régionaux apparaît à travers lécoulement des cours deau, le rabattement des eaux souterraines (18.4%), et lécoulement ascendant (4.7%).

     

Effects of fracture lineaments and in-situ rock stresses on groundwater flow in hard rocks: a case study from Sunnfjord, western Norway

Systematic field mapping of fracture lineaments observed on aerial photographs shows that almost all of these structures are positively correlated with zones of high macroscopic and mesoscopic fracture frequencies compared with the surroundings. The lineaments are subdivided into zones with different characteristics: (1) a central zone with fault rocks, high fracture frequency and connectivity but commonly with mineral sealed fractures, and (2) a damage zone divided into a proximal zone with a high fracture frequency of lineament parallel, non-mineralized and interconnected fractures, grading into a distal zone with lower fracture frequencies and which is transitional to the surrounding areas with general background fracturing. To examine the possible relations between lineament architecture and in-situ rock stress on groundwater flow, the geological fieldwork was followed up by in-situ stress measurements and test boreholes at selected sites. Geophysical well logging added valuable information about fracture distribution and fracture flow at depths. Based on the studies of in-situ stresses as well as the lineaments and associated fracture systems presented above, two working hypotheses for groundwater flow were formulated: (i) In areas with a general background fracturing and in the distal zone of lineaments, groundwater flow will mainly occur along fractures parallel with the largest in-situ rock stress, unless fractures are critically loaded or reactivated as shear fractures at angles around 30° to σ_H; (ii) In the influence area of lineaments, the largest potential for groundwater abstraction is in the proximal zone, where there is a high fracture frequency and connectivity with negligible fracture fillings. The testing of the two hypotheses does not give a clear and unequivocal answer in support of the two assumptions about groundwater flow in the study area. But most of the observed data are in agreement with the predictions from the models, and can be explained by the action of the present stress field on pre-existing fractures.     

Hydraulic properties and scale effects investigation in regional rock aquifers,south-western Quebec,Canada

This paper reports on the characterization of hydraulic properties of regional rock aquifers carried out within a groundwater resources assessment project in the St. Lawrence Lowlands of south-western Quebec. To understand the aquifer behavior at both the fracture level and at field scale, hydraulic investigations were carried out using various aquifer tests. The groundwater flow at the local scale is controlled mostly by the fracture system. Results of the constant-head injection tests show a weak decreasing trend of hydraulic conductivity with depth indicating that a major part of the groundwater flow occurs in the first meters of the rock sequence. At the regional scale, the equivalent porous media approach is applicable. The hydraulic conductivity measurements were correlated to the scale of the aquifer tests expressed with the investigated aquifer volume. A simple interpolation procedure for the hydraulic conductivity field was developed based on the distance between field measurements and the tested aquifer volumes. The regional distribution of the hydraulic conductivity for the major fractured aquifer units indicates that dolostone is the most permeable whereas sandstone and crystalline rocks are the least permeable units.

---

Resumen Este artículo trata de la caracterización de las propiedades hidráulicas en acuíferos regionales rocosos, la cual se llevó a cabo dentro del proyecto de evaluación de los recursos de agua subterránea en St. Lawrence Lowlands al suroeste de Quebec. Para entender el comportamiento del acuífero tanto a nivel de fractura como a escala del campo, se ejecutaron investigaciones hidráulicas usando varias pruebas de acuífero. El flujo del agua subterránea a escala local está controlado principalmente por el sistema de fracturas. Los resultados de las pruebas de inyección con cabeza constante muestran una tendencia decreciente débil de la conductividad hidráulica con la profundidad, indicando que la mayor parte del flujo de agua subterránea sucede en los primeros metros de la secuencia rocosa. A escala regional se puede aplicar la aproximación equivalente de medios porosos. Las medidas de la conductividad hidrálica fueron comparadas con la escala de las pruebas del acuífero, expresadas con el volumen investigado del acuífero. Se desarrolló un proceso de interpolación simple para la conductividad hidráulica de campo, con base en la distancia entre las medidas de campo y los volúmenes probados del acuífero. La distribución regional de la conductividad hidráulica, para las unidades acuíferas fracturadas mayores, indica que la dolomita es más permeable mientras que la arenisca y las rocas cristalinas son las unidades menos permeables.

---

Résumé Larticle présente les propriétés hydrauliques dun aquifère régional rocheux qui ont été mesurées dans le cadre dun projet concernant lestimation de la ressource en eau de la plaine de St. Lawrence située dans la partie sud-ouest de Quebec. Affin de comprendre le comportement de l› aquifère tant à l› échelle de fracture qu› à l› échelle régionale on a mené des investigations hydrauliques en utilisant des essais differents. À l› échelle locale lécoulement est déterminé en principal par le système des fractures. Les résultats des essais dinjection au niveau constant ont montré une tendence de décroissance de la conductivité hydraulique avec la profondeur ce qui indique que l› écoulement ait lieu en principal dans les premiers mètres de la structure rocheuse. À l› échelle règionale on peut appliquer lapproche de milieux poreux equivalent. Les mesures de conductivité hydraulique ont été corrélées avec l› échelle des essais de pompage exprimée en termes de volume investigué de l› aquifère. On a mis au point une méthode simple dinterpolation pour le champ de la conductivité hydraulique basée sur la distance entre les essais in situ et les volume investigués de l› aquifère. La distribution règionale de la conductivité hydraulique montre que pour la majorité des unités de roche fracturé les plus perméable sont les roches dolomitiques alors que les grés et les roches cristalines presentent la plus faible perméabilité.

     

水文地质概念模型的水动力场研究方法——"流场分析法"

水文地质概念模型建立的正确与否是地下水资源评价计算的关键,水文地质概念模型的水动力场研究方法——“流场分析法”是以地下水动力学和水均衡原理为基础,通过对水头函数的时空分布特征分析,判断出地下水系统的结构和功能,据此,科学概化出地下水系统的水文地质概念模型,为地下水资源的正确评价计算奠定基础。     

Small scale study of groundwater flow in a fractured carbonate-rock aquifer at the St-Eustache quarry, Québec, Canada

A small-scale hydrogeological study was conducted in a fractured carbonate-rock aquifer located in a quarry to relate groundwater flow to the fracture network. The field study in the St-Eustache quarry, which integrates structural surveys, well logging and hydraulic tests, showed that the most important features that affect groundwater flow in the sedimentary aquifer are high hydraulic conductivity horizontal bedding planes. Vertical fractures are abundant in the quarry and throughout the region, but they have a minor effect on groundwater flow. To have a significant impact on the flow regime and lead to vertical groundwater flow, the permeability of all vertical joints need to be enhanced compared to what was generally observed at the site. Such an increase in permeability could potentially occur where dissolution and fracturing is more intense or at stress release locations such as near the surface in the quarry.     

Contact zone permeability at intrusion boundaries: new results from hydraulic testing and geophysical logging in the Newark Rift Basin, New York, USA

Hydraulic tests and geophysical logging performed in the Palisades sill and the underlying sedimentary rocks in the NE part of the Newark Rift Basin, New York, USA, confirm that the particular transmissive zones are localized within the dolerite-sedimentary rock contact zone and within a narrow interval below this contact zone that is characterized by the occurrence of small layers of chilled dolerite. Transmissivity values determined from fluid injection, aquifer testing, and flowmeter measurements generally fall in the range of 8.1E-08 to 9.95E-06 m²/s and correspond to various scales of investigation. The analysis of acoustic and optical BHTV images reveals two primary fracture sets within the dolerite and the sedimentary rocks—subhorizontal fractures, intersected by subvertical ones. Despite being highly fractured either with subhorizontal, subvertical or both fracture populations, the dolerite above and the sedimentary rocks below the contact zone and the zone with the layers of chilled dolerite are significantly less conductive. The distribution of the particular conductive intervals is not a function of the two dominant fracture populations or their density but rather of the intrusion path of the sill. The intrusion caused thermal fracturing and cracking of both formations, resulting in higher permeability along the contact zone.     

The importance of alternative conceptual models for simulation of concentrations in a multi-aquifer system

Four different conceptual models based on alternative geological interpretations were formulated for a shallow 600 km² aquifer system in Denmark comprising Quaternary deposits. Each of the four models was calibrated against groundwater heads and discharge measurements through inverse modeling. Subsequently, the transport capabilities of the four models were compared to 32 concentration measurements of environmental tracers (tritium ³H, helium-3 ³He, chlorofluorocarbons CFC11, CFC12 and CFC113). The flow simulations showed only minor differences in spatial head distribution associated with alternative conceptualizations despite the complexity of the aquifer system and the significant differences in geological interpretations. The models, however, showed major differences in predictions of the age of the groundwater and environmental tracer concentrations, differences that are seen as an effect of model structure uncertainty, because no additional calibrations to these data were performed. A single conceptualization may be adequate in characterizing the natural behavior of a field system after calibration, because the calibration procedure is able to compensate for errors in the data or in the conceptual model through biased parameter values. However, once extrapolation beyond the calibration base is attempted, different conceptual model formulations result in significantly different results. Consequently, it is crucial to take model conceptual uncertainty into account when making predictions beyond the calibration base.     

A conceptual hydrogeological model of ophiolite hard-rock aquifers in Oman based on a multiscale and a multidisciplinary approach

Ophiolites are found all over the world: from the Alps to the Himalayas, in Cuba, Papua-New Guinea, New Caledonia, Newfoundland, etc. They are composed of hard rocks—basalt, dolerite, gabbro and peridotite, which are formed at the mid-oceanic ridges, with specific ridge-related tectonic fracturing and intense hydrothermal alteration. Their geological and thus their hydrogeological properties differ from those of both granite or classical gabbro and classical basaltic lava. A conceptual hydrogeological model of these hard-rock aquifers was developed based on the convergent results of a multidisciplinary approach at several spatial scales, from rock-sample (centimetre) to catchment (kilometre), on well-preserved ophiolite rocks in Oman. In ophiolite rocks, groundwater circulation takes place mostly in the fissured near-surface horizon (50 m thick), and, to a lesser degree, in the tectonic fractures. Hydrograph analysis (Water Resour Res 34:233–240, 1977), interpretation of numerous pumping tests using both classical Theis and dual porosity models [Water Resour Res 32:2733–2745, 1996; Comput Geosci J (in press)], and mercury porosity and hydraulic conductivity lab-measurements support the aquifer parameter estimates. The hydraulic conductivity K of the fissured horizon is estimated at 10^–5 to 10^–6 m/s for gabbro and dolerite, and 10^–7 m/s for peridotite. The storage coefficient S of the peridotite aquifer is estimated at 10^–3 and appears to be controlled mainly by microcracks (20 to 100 m wide). Tectonic fractures in the ophiolite have similar hydrodynamic properties regardless of lithology (10^–1<T<10^–4 m²/s and 10^–1<S<10^–3) though the probability of obtaining productive wells is two to three times greater in gabbro and dolerite than in peridotite. Some of the tectonic fractures produce small hydrothermal, hyperalkaline springs in the peridotite. The water budget and hydrochemistry of the Oman ophiolite are characterized and support the conceptual hydrogeological model. Despite low annual rainfall, a relatively low hydraulic conductivity and a significant storage coefficient explain why most of the streams in peridotite are perennial.

---

Résumé Les ophiolites se trouvent partout dans le monde : des Alpes à lHymalaya, à Cuba, Papouasie Nouvelle-Guinée, Nouvelle Calédonie, etc. Les ophiolites sont composées de roches dures basaltiques, de dolérites, gabbros et péridotites, qui sont formées aux rides mi-océaniques, avec, en relation avec les rides, une fracturation tectonique et une altération hydrothermale intense. Leurs propriétés géologiques et hydrogéologiques diffèrent de celles des granites ou des gabbros classiques ou des laves basaltiques classiques. Un modèle hydrogéologique conceptuel des aquifères de hard-rock a été développé sur les résultats convergents dune approche multidisciplinaire à différentes échelles spatiales, de léchantillon de roche de quelques centimètres à la taille du bassin versant en kilomètres, des ophiolites bien conservées dOman. Dans les roches ophiolitiques, la circulation des eaux souterraines est surtout localisée dans les fissures dun horizon de surface de plus dune cinquantaine de mètres et à un degré moindre dans les fractures tectoniques. Lanalyse des hydrographes (Water Resour Res 34:233–240, 1977), linterprétation des nombreux essais de pompage utilisant des modèles classiques de Theis et des modèles à double porosité (Water Resour Res 32:2733–2745, 1996; Comput Geosci J (Soumis), la porosité au mercure et la conductivité hydraulique calculée en laboratoire supportent lestimation des paramètres de laquifère. La conductivité hydraulique K de lhorizon fissuré est estimée entre 10^–5 et 10^–6 m/s pour les gabbros et les dolérites, et 10^–7 m/s pour les péridotites. Le coefficient demmagasinement S de laquifère de péridotite est estimé à 10^–3 et semble être contrôlé par les micro-fractures (20 à 100 m de largeur). Les fractures tectoniques dans les ophiolites possèdent des propriétés hydrodynamiques en regard de la lithologie (10-1<T<10-4 m²/s et 10-1/S/10-3) bien que la probabilité dobtenir des puits productifs est trois plus grand dans le grabbro et la dolérite que la péridétite. Quelques unes des fractures tectoniques produisent de petites sources hydrothermales hyperalcalines dans la péridotite. Le bilan hydrologique et hydrochimique de lophiolite dOman caractérise le modèle hydrogéologique conceptuel. Malgré des faibles précipitations annuelles, une relativement basse conductivité hydraulique et un coefficient demmagasinement significatif expliquent pourquoi la plus part des cours deau sur les péridotites ne sont jamais taris.

---

Resumen Las rocas ofiolitas se encuentran en muchas partes del mundo desde los Alpes a los Himilaya, Newfoundland, Nueva Caledonia etc. Son rocas duras compuestas de basalto, dolerita, gabro y peridotita que son formadas por los cordones medio-oceánicos y por fracturación relacionada con los cordones tectónicos específicos y alteración hidrotermal intensa. La geología y los propiedades hidrogeológicas varian entre las de granito y gabro clásico y lava de basalto clásico. Se desarrolló un modelo conceptual de los acuíferos de las rocas duras Basado en los resultados convergentes de un enfoce multidisciplinario a varias escalas a partir de muestras de roca (centímetro) a cuenca kilómetro) enfocado en las rocas ofilíticas bien preservados en Oman. En las rocas ofiolíticas la circulación de agua ocurre principalmente en el horizonte de fisuras cercanas a la superficie (>50 m de espesor) y en un grado menos importante en fracturas tectónicas. Los estimados de los parámetros de estos acuíferos se basan en análisis de hidrógrafos (Water Resour Res 34:233–240, 1977), en la interpretación de pruebas de bombeo usando los modelos clásicos de Theis y de porosidad doble (Hamm y Bidaux 1996; Lods y Gouze, [en prensa]), y porosidad de mecurio y conductividad hidraúlica medida en el laboratorio. La conductividad hidraúlica K del horizonte fisurado se estima en 10E-5 a 10E-6 m-s para gabro y dolerita y 10E-7 para peridotita. El coeficiente de abastecimiento S del acuífero peridotita se estima en 10–3 y parece ser controlado principalmente por microfracturas (20 to 100 m de ancho). Las fracturas tectónicas en la ofiolita tienen propiedades hidrodinámicas similares independientes de la litología (10–1<T<10–4 m²/s and 10–1<S<10–3) pero la probabilidad de obtener pozos productivos es dos o tres veces mayor en gabro que en dolerita o peridotita. Algunas de las fracturas tectónicas producen pequeñas descargas de agua hidrotermal hiperalcalina en la peridotita. El presupuesto e hidroquímica de las ofilitas de Omán se caracteriza y apoya en el modelo conceptual hidrogeológico. La baja lluvia anual, una conductividad hidraúlica relativamente baja y un coeficiente de abastecimiento importante explican porqué la mejoría de los drenajes en peridotita es peremne.

     

3D geologic framework models for regional hydrogeology and land-use management: a case study from a Quaternary basin of southwestern Quebec, Canada

During a regional hydrogeologic survey in the St. Lawrence Lowlands, Canada, a computer-based 3D Geologic Framework Model (GFM) was constructed to obtain a consistent representation of this typical Quaternary glaciated basin over a 1,400 km² area. Such a detailed stratigraphic reconstruction was needed because the Quaternary sediments control the recharge to the underlying regional fractured rock aquifer and also because buried granular aquifers are partly connected to the regional system. The objectives of this geomodeling effort are 1) to improve understanding of subsurface conditions above the regional aquifer and; 2) to provide a common stratigraphic framework for hydrogeologic applications. The method draws on knowledge-driven discrete modeling using gOcad, as well as standardization and quality control procedures to maximize the use of a multisource database. The resulting model represents the bedrock topography and the complex stratigraphic architecture of overlying sediments. The regional till aquitard, the marine clay aquiclude and the buried granular aquifers have been modeled with unprecedented details thus providing a well-constrained 3D hydrostratigraphic framework. The recharge zones of the rock aquifer represent about 35% of the study area. Buried granular aquifers are directly connected to the regional aquifer system over about 10% of the area. The model allows several applications such as assessing aquifer vulnerability and areal groundwater recharge rates; improving the GFM inter-operability with groundwater modeling systems would be the next logical step.

---

Resumen Se construyó un Modelo del Marco Geológico (GFM) 3D basado en computadora durante un levantamiento hidrogeológico regional en las Tierras Bajas de St. Lawrence, Canada para obtener una representación consistente de esta cuenca glaciar Cuaternaria típica en un área de 1,400 km². Se necesitó tal grado de reconstrucción estratigráfica debido a que los sedimentos Cuaternarios controlan la recarga del acuífero rocoso fracturado regional subyacente y también porque los acuíferos granulares enterrados están parcialmente conectados con el sistema regional. Los objetivos de este esfuerzo de modelizado geológico son: 1) mejorar el entendimiento de las condiciones subsuperficiales por encima del acuífero regional y; 2) aportar un marco estratigráfico común para aplicaciones hidrogeológicas. El método se basa en el conocimiento de modelizado discreto utilizando gOcad, así como también en estandarización y procedimientos de control de calidad para maximizar el uso de bancos de datos de fuentes múltiples. El modelo obtenido representa la topografía del macizo rocoso y la arquitectura de la compleja estratigrafía de los sedimentos sobreyacentes. El acuitardo regional de tillita, el acuicludo arcilloso marino y los acuíferos granulares enterrados se han modelizado con detalles sin precedentes aportando de este modo un marco hidroestratigráfico 3D bien definido. Las zonas de recarga del acuífero rocoso representan casi el 35% del área de estudio. Los acuíferos granulares enterrados están conectados directamente al sistema acuífero regional sobre aproximadamente el 10% del área. El modelo permite varias aplicaciones tal como evaluar la vulnerabilidad de acuíferos y ritmos de recarga areales de agua subterránea; el paso lógico siguiente sería mejorar la inter-operabilidad del GFM con los sistemas de modelizado de agua subterránea.

---

Résumé Dans le cadre dune étude dhydrogéologie régionale dans les basses terres du Saint-Laurent (Canada), un Modèle surfacique géologique en 3D (MSG) a été développé pour obtenir une représentation cohérente par ordinateur de ce bassin quaternaire sur une superficie de plus de 1,400 km². Une telle reconstitution stratigraphique était nécessaire en raison du contrôle quexercent les sédiments quaternaires sur la recharge de laquifère régional fracturé sous-jacent et aussi parce que les aquifères granulaires enfouis sont partiellement connectés au système régional. Les objectifs de cet effort de géomodélisation sont: 1) daméliorer la compréhension des conditions de sous-surface au-dessus de laquifère régional et; 2) de fournir un cadre stratigraphique commun à des fins de caractérisation hydrogéologique. La méthode repose sur lintégration de linterprétation géologique dans la modélisation discrète effectuée à laide du logiciel gOcad, ainsi que sur des procédures duniformisation et de contrôle de la qualité des données afin doptimiser lutilisation dune base de données multisources. Le MSG représente la topographie du roc et larchitecture stratigraphique des sédiments sus-jacents. Laquitard régional (till), laquiclude dargile marine et les aquifères granulaires ont été modélisés à un niveau de détail sans précédent, fournissant ainsi un cadre hydrostratigraphique solidement établi. Les zones de recharge de laquifère rocheux représentent environ 35% de la zone détude. Les aquifères granulaires enfouis sont connectés directement au système aquifère régional sur au moins 10% du territoire. Le MSG permet de multiples usages tels lévaluation de la vulnérabilité des aquifères et lestimation de la distribution de la recharge; la prochaine étape logique serait lamélioration de linteropérabilité avec les systèmes de modélisation numérique de lécoulement.

     

Three-dimensional hydrogeological modeling to assess the elevated-water-table technique for controlling acid generation from an abandoned tailings site in Quebec,Canada

The Manitou Mine sulphidic-tailings storage facility No. 2, near Val D’Or, Canada, was reclaimed in 2009 by elevating the water table and applying a monolayer cover made of tailings from nearby Goldex Mine. Previous studies showed that production of acid mine drainage can be controlled by lowering the oxygen flux through Manitou tailings with a water table maintained at the interface between the cover and reactive tailings. Simulations of different scenarios were performed using numerical hydrogeological modeling to evaluate the capacity of the reclamation works to maintain the phreatic surface at this interface. A large-scale numerical model was constructed and calibrated using 3 years of field measurements. This model reproduced the field measurements, including the existence of a western zone on the site where the phreatic level targeted is not always met during the summer. A sensitivity analysis was performed to assess the response of the model to varying saturated hydraulic conductivities, porosities, and grain-size distributions. Higher variations of the hydraulic heads, with respect to the calibrated scenario results, were observed when simulating a looser or coarser cover material. Long-term responses were simulated using: the normal climatic data, data for a normal climate with a 2-month dry spell, and a simplified climate-change case. Environmental quality targets were reached less frequently during summer for the dry spell simulation as well as for the simplified climate-change scenario. This study illustrates how numerical simulations can be used as a key tool to assess the eventual performance of various mine-site reclamation scenarios.

     

Correspondence analysis applied to a comparison of some rhyolitic zones in the Noranda area (Quebec, Canada)

Correspondence analysis was applied to 177 chemical analyses (major elements) of rocks from some rhyolitic zones of the Noranda area (Quebec). The objective was to grade the individual sample sections with respect to their attractiveness considering the chemical alterations related to mineralization. The results obtained show that it is possible to differentiate sections from different rhyolitic horizons.     

DRASTIC-Fm: a modified vulnerability mapping method for structurally controlled aquifers in the southern Gulf Islands,British Columbia,Canada

DRASTIC, the methodology for mapping the intrinsic vulnerability of aquifers, is modified to incorporate the structural characteristics of fractured bedrock aquifers. In these aquifers, groundwater flow is predominantly through fractures, with large-scale fracture zones and faults acting as primary conduits for flow at the regional scale. The methodology is applied to the southern Gulf Islands region of southwestern British Columbia, Canada. Bedrock geology maps, soil maps, structural measurements, mapped lineaments, water-well information and topographic data, assembled within a comprehensive GIS database, form the basis for assigning traditional DRASTIC indices, while adding the structural indices necessary for capturing the importance of regional structural elements in recharge and well capture zone determinations.     

Linking groundwater quality to residence times and regional geology in the St. Lawrence Lowlands,southern Quebec,Canada

The assessment of groundwater quality in shallow aquifers is of high societal relevance given that large populations depend directly on these water resources. The purpose of this study was to establish links between groundwater quality, groundwater residence times, and regional geology in the St. Lawrence Lowlands fractured bedrock aquifer. The study focuses on a 4500 km² watershed located in the St. Lawrence Lowlands of the province of Quebec in eastern Canada. A total of 150 wells were sampled for major, minor, and trace ions. Tritium (³H) and its daughter element, ³He, as well as radiocarbon activity (A¹⁴C) were measured in a subset of wells to estimate groundwater residence times. Results show that groundwater evolves from a Ca–HCO₃ water type in recharge zones (i.e., the Appalachian piedmont) to a Na–HCO₃ water type downgradient, toward the St. Lawrence River. Locally, barium (Ba), fluoride (F), iron (Fe), and manganese (Mn) concentrations reach 90, 2, 18, and 5.9 mg/L respectively, all exceeding their respective Canadian drinking water limits of 1, 1.5, 0.3, and 0.05 mg/L. Release of these elements into groundwater is mainly controlled by the groundwater redox state and pH conditions, as well as by the geology and the duration of rock–water interactions. This evolution is accompanied by increasing ³H/³He ages, from 4.78 ± 0.44 years upgradient to more than 60 years downgradient. Discrepancies between calculated ³H/³He and ¹⁴C water ages (the latter ranging from 280 ± 56 to 17,050 ± 3410 years) suggest mixing between modern water and paleo-groundwater infiltrated through subglacial recharge when the Laurentide Ice Sheet covered the study area, and during the following deglaciation period. A linear relationship between ³H activity and corrected ¹⁴C versus Mg/Ca and Ba support a direct link between water residence time and the chemical evolution of these waters. The Ba, F, Fe, and Mn concentrations in groundwater originate from Paleozoic rocks from both the St. Lawrence Platform and the Appalachian Mountains. These elements have been brought to the surface by rising hydrothermal fluids along regional faults, and trapped in sediment during their deposition and diagenesis due to reactions with highly sulfurous and organic matter-rich water. Large-scale flow of meltwater during subglacial recharge and during the subsequent retreat of the Laurentide Ice Sheet might have contributed to the leaching of these deposits and their enrichment in the present aquifers. This study brings a new and original understanding of the St. Lawrence Lowlands groundwater system within the context of its geological evolution.     

加拿大拉布拉多地槽DSO型铁矿床成矿规律及找矿方法研究

DSO(Direct Shipping Ore)型铁矿床即直运型铁矿床,指加拿大拉布拉多地槽区Fe品位大于50%(w(SiO2)18%)的赤铁矿矿床,因其矿石品位高、经济效益突出而备受关注。对该地区DSO矿床区域成矿背景、典型矿床研究表明:(1)DSO矿体均赋存于苏克曼组含铁建造中;(2)DSO矿体的产出位置主要受构造控制,构造为热液流体的运移提供了通道,也为矿质的重结晶及沉淀提供了空间;(3)软质DSO成因属深部热液-表生改造型,硬质DSO成因则属热液型。根据DSO的成矿特点,可将拉布拉多地槽DSO找矿靶区分为谢弗维尔软质DSO找矿靶区、铁矿湖-阿提卡玛根软质DSO找矿靶区和阿斯特雷—索耶湖硬质DSO找矿靶区。找矿实践中,应通过地质填图、物探等方法对拉布拉多地槽中部地区的岩性和构造进行控制,重点关注NW向大断裂与次级断裂、次级褶皱、玄武岩与含铁建造岩性转换面的叠加处和地球物理低磁高重区。综合地质和地球物理方法优选靶区进行槽探、钻探相结合的快速验证,可有效提高找矿效率,实现拉布拉多地槽区DSO找矿工作的快速突破。     

An empirical model for estimating hydraulic conductivity of highly disturbed clastic sedimentary rocks in Taiwan

This paper presents the establishment of an empirical HC model for estimating rock mass hydraulic conductivity of highly disturbed clastic sedimentary rocks in Taiwan using high-resolution borehole acoustic televiewer and double packer hydraulic tests. Four geological parameters including rock quality designation (RQD), depth index (DI), gouge content designation (GCD), and lithology permeability index (LPI) were adopted for establishing the empirical HC model. To verify rationality of the proposed HC model, 22 in-situ hydraulic tests were carried out to measure the hydraulic conductivity of the highly disturbed clastic sedimentary rocks in three boreholes at two different locations in Taiwan. Besides, the model verification using another borehole data with four additional in-situ hydraulic tests from similar clastic sedimentary rocks was also conducted to further verify the feasibility of the proposed empirical HC model. The field results indicated that the rock mass in the study area has a conductivity between the order of 10^− 10 m/s and 10^− 6 m/s at the depth between 34 m and 275 m below ground surface. Results demonstrate that the empirical HC model may provide a useful tool to predict hydraulic conductivity of the highly disturbed clastic sedimentary rocks in Taiwan based on measured HC-values.     

A conceptual model of hydraulic networks in carbonate rocks, illustrated by examples from Poland

 Hydraulic networks of carbonate rocks usually consist of three types of spaces: pores, fissures, and caverns; sometimes filled forms also are present. At sufficiently large scales, the spaces are assumed to be homogeneous with arbitrary boundaries and to be characterized by mean values of parameters. A conceptual model and its schematic presentation demonstrate the importance of particular parameters for flow and solute transport, and they explain the responses of pumping rates and water levels to changes in boundary conditions. Three examples characterize some typical networks: (1) chalk in eastern Poland, representing fissured aquifers with a high storage volume for solutes in stagnant matrix water; (2) Middle Triassic limestones in southern Poland, representing fissured-cavernous aquifers; however, a large storage volume for solute transport is in mostly stagnant matrix water; and (3) Triassic dolomites in the same area, representing porous-fissured-cavernous aquifers with a considerable storage volume in the matrix, which can partly be released (20%), considerably contributing to the total specific yield. For the last two types of aquifers, the hydraulic conductivity near pumping wells and mine galleries usually quickly increases in time due to suffosion of filled forms, which leads to a temporal large increase of the inflows. Received, April 1997 Revised, January 1998, August 1998 Accepted, May 1998     

A comparison of electrical and electromagnetic methods for the detection of hydraulic pathways in a fractured rock aquifer,Clare Valley,South Australia

Within fractured rock, the irregular and often unpredictable distribution and geometry of hydraulically conductive fractures produces large spatial variations in bore yield and groundwater quality. As fractures act as conduits for flow of both groundwater and electrical charge, methods which can efficiently detect the distribution of electrical pathways can be used to infer characteristics of significant hydrological parameters. This study compares the capabilities and limitations of electrical data obtained from direct current (DC) and electromagnetic (EM) surface azimuthal measurements, and from DC borehole-to-surface and cross-borehole measurements, for the interpretation of major hydrological structures in Clare Valley, South Australia. Electrical and EM surface methods are limited by poor depth sensitivity and the presence of conductive overburden, but provide useful tools for determining directional variations in resistivity at sites lacking bedrock exposure and boreholes. Application of borehole-to-surface methods yielded a better-resolved interpretation of sub-vertical fracture strike and was useful in identifying lateral variations in bedrock heterogeneity. Improved flexibility and sensitivity to measurements at depth permitted cross-borehole electrical tomography data to be used in reconstructing the spatial distribution of sub-horizontal, laterally extensive, electrically conductive zones. While the technique is restricted to small-scale sites with multiple boreholes, inferences can be made on fluid connections over a much larger regional scale. It is important to note, however, that while electrical methods provide valuable information about in-situ hydraulic pathways, they do not provide a complete hydraulic characterisation. Such a task requires integration of surface and borehole geophysics, geologic mapping, sampling and pumping tests of wells with packed-off intervals.

---

Resumen La geometría, y con frecuencia impredecible, distribución irregular de fracturas hidráulicamente conductivas dentro de roca fracturada genera variaciones espaciales grandes en producción de pozos y calidad de agua subterránea. Debido a que las fracturas actúan como conductos de flujo de agua y carga eléctrica, los métodos que detectan eficientemente la distribución de trayectorias eléctricas pueden utilizarse para inferir las características de parámetros hidrológicos significativos. Este estudio compara las capacidades y limitaciones de datos eléctricos obtenidos de mediciones azimutales superficiales electromagnéticas (EM) y de corriente directa (DC), y de mediciones de DC realizadas en la superficie y pozos así como mediciones realizadas entre pozos para la interpretación de estructuras hidrológicas principales en el Valle Clare, sur de Australia. Los métodos superficiales eléctricos y EM están limitados por sensitividad de profundidad pobre y la presencia de cubierta conductiva, pero aportan herramientas útiles para determinar variaciones direccionales de resistividad en sitios que carecen de pozos y afloramientos rocosos. La aplicación de métodos superficiales y de pozos aportan una mejor interpretación del rumbo de fracturas sub-verticales y fue útil en identificar variaciones laterales en la heterogeneidad del macizo rocoso. El mejoramiento de flexibilidad y sensitividad en las mediciones profundas permitió que los datos de tomografía eléctrica de los pozos fuera utilizado en la reconstrucción de la distribución espacial de zonas eléctricas conductivas, sub-horizontales y lateralmente extensas. Aunque la técnica está restringida a sitios de pequeña escala con múltiples pozos, puede realizarse inferencias sobre relaciones entre fluidos en una escala regional mucho más grande. Sin embargo, es importante notar que aunque los métodos eléctricos aportan información valiosa acerca de las trayectorias hidráulicas in-situ, aún no proporcionan una caracterización hidráulica completa. Esta tarea requiere integrar geofísica superficial y de pozos, mapeo geológico, muestreo y pruebas de bombeo en pozos con intervalos sin empaque.

---

Résumé Au sein des roches fracturées, lirrégularité et limprédictible distribution et géométrie des fractures par lesquelles sécoulent leau souterraine produit de larges variations spatiales entre les débits des forages et les paramètres de la qualité de leau. Comme les fractures conduisent et leau et les charges électriques, les méthodes qui peuvent de manière efficiente détecter la distribution des courants électriques peuvent être utilisées pour analyser les caractéristiques des principaux paramètres hydrologiques. Cette étude compare les capacités et les limites des données produites par Courant Direct (DC en Anglais) et par mesure Electromagnétique azimutale de surface (EM en Anglais), et par DC entre surface et forage et entre plusieurs forages, pour linterprétation de la structure hydrologique de la Vallée de Clare, Australie du Sud. Les méthodes électriques et électromagnétiques sont limitées par la faible sensitivité à la profondeur et la présence dune couverture conductrice, mais procure des outils utiles pour déterminer les variations directionnelles en terme de résistivité là où la roche nest pas affleurante et dans les puits. Lapplication de la méthode dinvestigation «forage à surface» apporte une meilleure interprétation des fractures sub-verticales et des hétérogénéités latérales. Les tomographies électriques entre forage bénéficient de la flexibilité et de la sensibilité des mesures en profondeur, et permettent de délimiter lextension latérale des hétérogénéités sub-horizontales de zones conductrices. Alors que la technique est restreinte à de petits sites comprenant de nombreux forages, il est possible de reconstituer les connections hydrauliques à des échelles régionales. Il est important de noter que les méthodes fournissent des informations intéressantes mais pas des caractérisations hydrauliques complètes. Pour cela les données pourraient être complétées par des études plus poussées intégrant les différentes prospections géophysique, les données des cartes géologiques, des échantillonnages et des essais de pompage à différents intervalles de profondeur.

     

Development of a conceptual hydrogeological model for the evaluation of residence times of water in soil and groundwater: the state of Hesse case study, Germany

A regional conceptual hydrogeological model has been developed for evaluating residence times of both, percolate water in the unsaturated zone and groundwater in upper aquifers. The model is based on digitally geo-data bases available at the regional level and has been applied for the entire Federal State of Hesse (Germany) with a spatial resolution of 60?×?60?m. Residence times determined for unconsolidated rock areas typically ranged between 10 and 25?years, whereas residence times of <5?years were assessed for consolidated rock areas. With regard to the implementation of the EU Water Framework Directive, the determined residence times may help to assess the time periods between the introduction of well-targeted groundwater protection measures and their impact on groundwater and surface water quality, respectively.     

Temporal analysis of groundwater nitrate concentrations from wells in Prince Edward Island,Canada: application of a linear mixed effects model

Changes in nitrate concentration in groundwater from wells in Prince Edward Island, Canada were investigated over time using two datasets. Temporal trends in groundwater nitrate concentrations were assessed annually during 1981–1996 (1,299 observations), and both seasonally and monthly during 1988–1991 (1,868 observations). Data were analysed using linear mixed models with random effects and correlation structures. The average nitrate concentration in the monthly dataset was 3.99 mg/L as NO₃–N, with January, May, and November concentrations being higher (p?=?0.018). A seasonal effect was present when season was combined with land use type in an interaction term (p?=?0.004). Wells located in agricultural areas had greater nitrate concentrations than urban areas, which in turn, had greater values than low human-impact areas. Row-cropped areas had higher groundwater nitrate concentrations in the summer, whereas manure storage areas were higher in the spring and autumn. Nitrate in groundwater in areas with low human impact and with centralized sewage disposal infrastructure remained relatively low and stable throughout the seasons. There was no significant annual trend (p?=?0.954), but for individual sites, 9.6% significantly increased in nitrate concentration over time, and 6.6% significantly decreased over time.