Modeling stream-aquifer interactions in a shallow aquifer, Choele Choel Island, Patagonia, Argentina

A groundwater/surface-water interaction model was developed for the shallow alluvial aquifer of the Choele Choel Island in Patagonia, Argentina. In this semiarid climate, agriculture is sustained by an irrigation/drainage system. During the irrigation season, seepage losses through unlined distribution canals in irrigated fields contribute to elevated groundwater levels, jeopardizing fruit productivity in some areas. Moreover, high stream stages during the irrigation season interfere with groundwater drainage. The model utilized MODFLOW and its stream package, and was successfully calibrated for a historical irrigation season. Modeling results indicate that drainage through streams is significantly higher than drainage through artificial drains. The stream/aquifer relationship proved very responsive to water table rises caused by irrigation. This response manifested as changes in the gaining/losing character of stream reaches. A synthetic run aimed at isolating the effect of streamflow changes on groundwater levels showed that the effect of higher streamflows dissipates toward the interior of the island, disappearing completely at the island center. Even though some results were qualitative, the model helped to provide a better understanding of the coupled system to elucidate some of the causes of a rising water table on the island.

A regional-scale groundwater flow model for the Leon-Chinandega aquifer,Nicaragua

Groundwater flow in the Leon-Chinandega aquifer was simulated using transient and steady-state numerical models. This unconfined aquifer is located in an agricultural plain in northwest Nicaragua. Previous studies were restricted to determining groundwater availability for irrigation, overlooking the impacts of groundwater development. A sub-basin was selected to study the groundwater flow system and the effects of groundwater development using a numerical groundwater flow model (Visual MODFLOW). Hydrological parameters obtained from pumping tests were related to each hydrostratigraphic unit to assign the distribution of parameter values within each model layer. River discharge measurements were crucial for constraining recharge estimates and reducing the non-uniqueness of the model calibration. Steady-state models have limited usefulness because of the major variation of recharge and agricultural pumping during the wet and dry seasons. Model results indicate that pumping induces a decrease in base flow, depleting river discharge. This becomes critical during dry periods, when irrigation is highest. Transient modeling indicates that the response time of the aquifer is about one hydrologic year, which allows the development of management strategies within short time horizons. Considering further development of irrigated agriculture in the area, the numerical model can be a powerful tool for water resources management.     

East African topography and volcanism explained by a single,migrating plume

Anomalous topographic swells and Cenozoic volcanism in east Africa have been associated with mantle plumes.Several models involving one or more fixed plumes beneath the northeastward migrating African plate have been suggested to explain the space-time distribution of magmatism in east Africa.We devise paleogeographically constrained global models of mantle convection and,based on the evolution of flow in the deepest lower mantle,show that the Afar plume migrated southward throughout its lifetime.The models suggest that the mobile Afar plume provides a dynamically consistent explanation for the spatial extent of the southward propagation of the east African rift system(EARS),which is difficult to explain by the northeastward migration of Africa over one or more fixed plumes alone,over the last≈45 Myr.We further show that the age-progression of volcanism associated with the southward propagation of EARS is consistent with the apparent surface hotspot motion that results from southward motion of the modelled Afar plume beneath the northeastward migrating African plate.The models suggest that the Afar plume became weaker as it migrated southwards,consistent with trends observed in the geochemical record.     

Effects of intraborehole flow on groundwater age distribution

Environmental tracers are used to estimate groundwater ages and travel times, but the strongly heterogeneous nature of many subsurface environments can cause mixing between waters of highly disparate ages, adding additional complexity to the age-estimation process. Mixing may be exacerbated by the presence of wells because long open intervals or long screens with openings at multiple depths can transport water and solutes rapidly over a large vertical distance. The effect of intraborehole flow on groundwater age was examined numerically using direct age transport simulation coupled with the Multi-Node Well Package of MODFLOW. Ages in a homogeneous, anisotropic aquifer reached a predevelopment steady state possessing strong depth dependence. A nonpumping multi-node well was then introduced in one of three locations within the system. In all three cases, vertical transport along the well resulted in substantial changes in age distributions within the system. After a pumping well was added near the nonpumping multi-node well, ages were further perturbed by a flow reversal in the nonpumping multi-node well. Results indicated that intraborehole flow can substantially alter groundwater ages, but the effects are highly dependent on local or regional flow conditions and may change with time.     

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%).

     

Effect of advective mass transfer on field scale fluid and solute movement: Field and modeling studies at a waste disposal site in fractured rock at Oak Ridge National Laboratory, Tennessee, USA

Advective mass transfer is a pore scale mass-transfer process that affects fluid and solute movement between pore domains such as fracture and matrix in a structured porous medium. Mechanistically similar to advection in the advection-dispersion of solutes in non-structured porous medium, it redistributes solutes by moving solute and solvent simultaneously between pore domains. While there is much research on diffusive mass transfer that is often referred to as matrix diffusion, there is a lack of information and study for advective mass transfer in the literature. The objective of this research is to study the effects of advective mass transfer on fluid and solute movement between pore domains. First, field hydraulic measurements at a waste disposal site in fractured rock at Oak Ridge National Laboratory (ORNL), Tennessee, USA, are used to calibrate a fracture-matrix, two-pore-domain groundwater flow model. Latin-hypercube sensitivity analysis suggests that the uncertainty of the calibrated model parameters is small and the calibrated flow model is nearly the optimal. Fracture spacing thus obtained is used to calculate diffusive mass transfer coefficients. The individual effects of advective and diffusive mass transfer on solute movement are then quantitatively evaluated. The calculations indicate that pore structure conceptual models may significantly affect the role of advective mass transfer in field and pore-scale mass transfer. In the particular ORNL field site and with a fracture-matrix pore structure model, contribution of advective mass transfer to solute mass movement is about three to eight orders of magnitude smaller than that of diffusive mass transfer.

---

Resumen La transferencia de masa advectiva es un proceso de transferencia de masa en escala intersticial que afecta el movimiento de fluido y soluto entre ámbitos porosos tal como fractura y matriz en un medio estructurado poroso. Este proceso, mecánicamente similar a la advección en la dispersión-advección de solutos de medios porosos no estructurados, redistribuye los solutos mediante el movimiento simultáneo de soluto y solvente entre ámbitos porosos. Mientras que existe bastante investigación en transferencia difusiva de masa que frecuentemente se conoce como difusión en matriz, existe falta de información y estudio de transferencia advectiva de masa en la literatura. El objetivo de esta investigación es estudiar los efectos de la transferencia advectiva de masa en el movimiento de fluido y soluto entre ámbitos porosos. Primero se utilizaron mediciones hidráulicas de campo en un sitio de depósito de residuos en roca fracturada en el Laboratorio Nacional Oak Ridge (ORNL), Tennessee, USA, para calibrar un modelo de flujo de agua subterránea de ámbito de dos poros fractura-matriz. Análisis de sensitividad hipercúbico-latino sugieren que la incertidumbre de los parámetros del modelo calibrado es pequeña y que el modelo de flujo calibrado es aproximadamente el óptimo. El espaciamiento de fracturas así obtenido se utiliza para calcular los coeficientes de transferencia de masa difusiva. Luego se evalúa cuantitativamente los efectos individuales de transferencia de masa advectiva y difusiva en el movimiento de soluto. Los cálculos indican que los modelos conceptuales de estructura porosa pueden afectar significativamente el papel de transferencia de masa advectiva en escalas de campo e intersticial de transferencia de masa. En el sitio de campo específico ORNL y con un modelo de estructura porosa de matriz-fractura, la contribución de transferencia de masa advectiva al movimiento de masa soluto es aproximadamente tres a ocho órdenes de magnitud más pequeño que la contribución por transferencia de masa difusiva.

---

Résumé Résumé Le transfert advectif de masse est un processus de transfert de masse à léchelle du pore qui affecte les mouvements du fluide et du soluté entre les différents domaines de pores, tel les fractures et la matrice dans un milieu poreux structuré. Mécaniquement similaire à ladvection dans le concept advection-dispersion de solutés dans les milieux non structurés, ce transfert redistribue les solutés simultanément avec le fluide entre les différents domaines poreux. Alors que de nombreuses recherches portent sur les transferts de masse par diffusion, se référant généralement à une diffusion par la matrice, il y a un grand manquement dinformations et détudes sur les transferts de masse par advection dans la littérature. Lobjectif de cette recherche est détudier leffet du transfert advectif de masse de fluide et de soluté entre les différents domaines poreux. Premièrement, les mesures hydrauliques de terrain sur la décharge en milieu fracturé du laboratoire national dOak Ridge ORNL, Tennessee, USA, sont utilisées pour calibrer un modèle hydrogéologique à double porosité fracture-matrice. Lanalyse de sensibilité latin-hypercube suggère que lincertitude sur les paramètres du modèle est faible et que le calibrage est pratiquement optimal. Lespace de fracture résultant permet de calculer les coefficients de transfert de masse par diffusion. Les effets individuels de ladvection et de la diffusion sur les mouvements de solutés sont dés lors évalués. Les calculs indiquent que le modèle conceptuel de la structure des pores peuvent significativement affectés le rôle du transfert advectif de masse à léchelle du pore et du terrain. Dans le cas du site de lORNL et avec un modèle structuré fracture-matrice, la contribution de ladvection au transfert de masse est de lordre de trois-huitième du transfert de masse par diffusion.

     

Sustainable development and management of an aquifer system in western Turkey

This study presents the establishment of sustainable development and management policies for the Küçük Menderes River Basin aquifer system in western Turkey. Geological, hydrogeological, and geophysical data are used conjunctively to define various hydrogeological units and their geometry. Distributions of hydraulic-parameter values and recharge are estimated by geostatistical methods and hydrologic simulations, respectively. A finite-difference groundwater flow model is used to represent the unconfined flow in the aquifer system. The model has been calibrated under steady state and transient conditions. The resulting model was used to test seven management scenarios for a planning period of 21 years to determine the so-called safe yield and sustainable yield of the aquifer system and to investigate the potential impacts of four planned surface water reservoirs on groundwater resources in the basin. The results demonstrate that the continuation of the present pumping rates exceeds both the safe yield and the sustainable yield of the aquifer system. Consequently, the growing need for irrigation water should be met by the construction of the planned surface water reservoirs and the implementation of efficient water management policies and plans.

---

Resumen Este estudio presenta el establecimiento de políticas sustentables de desarrollo y gestión en el sistema acuífero de la cuenca del río Küçük Menderes, al Oeste de Turquía, para lo que se ha utilizado datos geológicos, hidrogeológicos y geofísicos de forma conjunta de cara a definir diversas unidades hidrogeológicas y su geometría. La distribución de los parámetros hidráulicos y de la recarga ha sido estimada mediante métodos geoestadísticos y simulaciones hidrológicas, respectivamente. Se ha empleado un modelo de las aguas subterráneas en diferencias finitas para representar el flujo no confinado en el sistema acuífero, el cual se ha calibrado bajo condiciones estacionarias y transitorias. El modelo resultante ha sido usado para contrastar siete escenarios de gestión durante un período de planificación de 21 años con el fin de determinar el punto de explotación segura y sustentable del sistema acuífero, así como para investigar los impactos potenciales sobre los recursos subterráneos de los cuatro embalses superficiales que se hallan en proyecto. Los resultados demuestran que el mantenimiento de las tasas actuales de extracción del acuífero supera tanto el régimen de bombeo seguro como el sustentable del sistema. Por consiguiente, se debería compensar la necesidad creciente de agua para riego con la construcción de los reservorios superficiales previstos y con la implementación de políticas y planes de gestión eficiente del agua.

---

Résumé Cette étude présente la proposition dune politique de développement et de gestion durables du système aquifère du bassin du Petit Mendérès dans louest de la Turquie. Des données géologiques, hydrogéologiques et géophysiques ont été utilisées conjointement pour définir les différentes unités hydrogéologiques et leur géométrie. Les distributions des paramètres hydrauliques et de la recharge ont été estimées respectivement par des méthodes géostatistiques et des simulations hydrologiques. Un modèle découlement souterrain aux éléments finis a été utilisé pour représenter lécoulement non captif dans le système aquifère. Le modèle a été calibré dans des conditions de régimes permanent et transitoire. Le modèle résultant a servi à tester sept scénarios de gestion pour une période de programmation de 21 ans, afin de déterminer les débits de prélèvement sûr et durable dans le système aquifère et détudier les impacts potentiels de quatre réservoirs deau de surface en projet sur les eaux souterraines du bassin. Les résultats montrent que la poursuite des prélèvements au débit actuel excède aussi bien le débit dexploitation de sécurité que celui durable pour le système aquifère. Par conséquent, les besoins croissants deau pour lirrigation doivent être satisfaits par la construction des réservoirs projetés et par la mise en place de politiques et de plans de gestion de leau efficace.

     

Hydrochemical tracers in the middle Rio Grande Basin,USA: 2. Calibration of a groundwater-flow model

The calibration of a groundwater model with the aid of hydrochemical data has demonstrated that low recharge rates in the Middle Rio Grande Basin may be responsible for a groundwater trough in the center of the basin and for a substantial amount of Rio Grande water in the regional flow system. Earlier models of the basin had difficulty reproducing these features without any hydrochemical data to constrain the rates and distribution of recharge. The objective of this study was to use the large quantity of available hydrochemical data to help calibrate the model parameters, including the recharge rates. The model was constructed using the US Geological Surveys software MODFLOW, MODPATH, and UCODE, and calibrated using ¹⁴C activities and the positions of certain flow zones defined by the hydrochemical data. Parameter estimation was performed using a combination of nonlinear regression techniques and a manual search for the minimum difference between field and simulated observations. The calibrated recharge values were substantially smaller than those used in previous models. Results from a 30,000-year transient simulation suggest that recharge was at a maximum about 20,000 years ago and at a minimum about 10,000 years ago.

---

Resumen La calibración de un modelo de aguas subterráneas con el apoyo de datos hidroquímicos ha demostrado que la recarga relativamente baja en la cuenca media del Río Grande es probablemente responsable de una depresión de aguas subterráneas en el centro de la cuenca y de la presencia de una cantidad considerable de agua del Río Grande en el acuífero del Grupo Santa Fe. Los modelos propuestos con anterioridad para la cuenca tenían dificultades para reproducir estas características ya que no tenían datos hidroquímicos que permitieran delimitar los ritmos y distribución de recarga. El objetivo del presente estudio consistió en utilizar una gran cantidad de datos hidroquímicos disponibles para ayudar a calibrar los parámetros del modelo, incluyendo los ritmos de recarga. El modelo se construyó utilizando los modelos MODFLOW, MODPATH, y UCODE del USGS, mientras que la calibración se realizó en base a concentraciones de ¹⁴C y a la posición de ciertas zonas definidas con los datos hidroquímicos. La estimación de parámetros se realizó en base a una combinación de técnicas de regresiones no lineares y a una búsqueda a viva fuerza del error mínimo entre los datos observados y los simulados. Los valores de recarga calibrados fueron significativamente más bajos que los estimados en los modelos anteriores. Los resultados de una simulación transitoria de 30,000 años sugieren que la recarga durante la última glacial máxima (LGM) fue diez veces el ritmo moderno, pero que la recarga que ocurrió inmediatamente después de la LGM fue más baja que el ritmo moderno.

---

Résumé Le calibrage dun modèle hydrogéologique avec laide de données hydrochimiques a démontré que la recharge relativement faible dans le Grand Bassin du Middle Rio est vraisemblablement responsable dune dépression des eaux souterraines dans le centre du bassin et de la présence dune quantité substantielle deau du Rio Grande dans laquifère du Groupe de Santa Fe. Les modèles antérieurs avaient des difficultés à reproduire ses conclusions sans laide de données hydrochimiques pour contraindre les taux et la distribution de la recharge. Lobjectif de cette étude était dutiliser une grande quantité de données hydrochimiques permettant de calibrer les paramètres du modèle, et notamment les taux de recharge. Le modèle a été construit avec les logiciels MODFLOW, MODPATH et UCODE, et calibré en utilisant les concentrations en ¹⁴C et la position de certaines zones définies par les données hydrochimiques. Lestimation de certains paramètres a été réalisée en utilisant une combinaison de techniques de régression non linéaire et une méthode de recherche exhaustive (Brute Force Search) de lerreur minimum entre les résultats des observations et les simulations. Les valeurs de la recharge calibrée sont substantiellement plus basses que celles estimées dans les modèles antérieurs. Les résultats dune simulation en régime transitoire sur 30.000 ans suggèrent que la recharge au maximum de la dernière glaciation (last glacial maximum, LGM) était 10 fois supérieure au taux actuel, mais que la recharge qui a suivit la LGM était plus bas que la recharge actuelle.

     

Numerical modeling of contaminant transport resulting from dissolution of a coal-tar pool in an experimental aquifer

A previously developed two-dimensional numerical model is further developed for simulating the transport of dissolved contaminants originating from dissolution of a coal tar pool in a stratified, saturated porous medium. The model is used to simulate contaminant transport resulting from a rectangular-shaped coal-tar-pool dissolution experiment conducted in a large-scale experimental aquifer. The experimental porous medium consists of two sand strata, a high-hydraulic-conductivity upper stratum and a low-hydraulic-conductivity bottom stratum. The experiment was conducted to a time of 354 days and the groundwater velocity was changed several times within this duration. Model simulations show good agreement against observed contaminant concentrations, and simulations show that dissolved solute below the pool migrated deeper into the bottom stratum as compared to the upper stratum. Furthermore, simulations also suggest that contaminant concentrations in the lower stratum never reached quasi steady-state during the experimental time frame.     

Groundwater flow through Pleistocene glacial deposits in the rapidly urbanizing Rouge River–Highland Creek watershed,City of Scarborough,southern Ontario,Canada

The city of Scarborough lies on the eastern margin of the Greater Toronto Area of southern Ontario, Canada, along the northern coastline of Lake Ontario. The City has a population of 500,000 and is presently one of the fastest growing communities in Canada. The City is expanding northwards onto rural land on the south slope of the large Pleistocene glacial Oak Ridges Moraine system. The moraine system is underlain by a thick (150 m) succession of tills, sands and gravels and is a regionally-significant recharge area for three principle aquifer systems that discharge to numerous watercourses that flow to Lake Ontario. Protection of deeper aquifers from surface-generated urban contaminants is a particular concern. A groundwater flow model using Visual MODFLOW was developed for the 350-km² Rouge River–Highland Creek (RRHC) drainage basin using an extensive GIS-based collection of subsurface geological, geophysical and hydrogeological data, maps of land use and surficial geology. The RRHC model was calibrated against point water level data, known potentiometric surfaces of the principal aquifers and baseflow measurements from streamflow gauging stations and determined to be within acceptable limits. Water balance calculations indicate that 70% of the basin recharge (106,000 m³/day) enters the Upper Aquifer along the crest and immediate flanks of the Oak Ridges Moraine. To the south, Upper Aquifer water moving through fractured till aquitards accounts for more than 75% of recharge to deeper aquifers. Water quality data confirm previous observations that urban- and rural-sourced contaminants (chlorides and nitrates) present in Upper Aquifer waters are moving rapidly into deeper aquifers. Some 83% of total RRHC recharge water is ultimately discharged as baseflow to creeks draining to Lake Ontario; the remainder discharges to springs and along eroding lakeshore bluffs. Model results demonstrate that deeper aquifers are poorly protected from urban contaminants and that long-term protection of ground and surface water quality has to be a priority of municipal planners if the resource is not to be severely degraded. Electronic Publication     

Three-dimensional variable-density flow simulation of a coastal aquifer in southern Oahu, Hawaii, USA

Three-dimensional modeling of groundwater flow and solute transport in the Pearl Harbor aquifer, southern Oahu, Hawaii, shows that the readjustment of the freshwater–saltwater transition zone takes a long time following changes in pumping, irrigation, or recharge in the aquifer system. It takes about 50 years for the transition zone to move 90% of the distance to its new steady position. Further, the Ghyben–Herzberg estimate of the freshwater/saltwater interface depth occurred between the 10 and 50% simulated seawater concentration contours in a complex manner during 100 years of the pumping history of the aquifer. Thus, it is not a good predictor of the depth of potable water. Pre-development recharge was used to simulate the 1880 freshwater-lens configuration. Historical pumpage and recharge distributions were used and the resulting freshwater-lens size and position were simulated through 1980. Simulations show that the transition zone moved upward and landward during the period simulated.Previous groundwater flow models for Oahu have been limited to areal models that simulate a sharp interface between freshwater and saltwater or solute-transport models that simulate a vertical aquifer section. The present model is based on the US Geological Surveys three-dimensional solute transport (3D SUTRA) computer code. Using several new tools for pre- and post-processing of model input and results have allowed easy model construction and unprecedented visualization of the freshwater lens and underlying transition zone in Hawaiis most developed aquifer.

---

Electronic Supplementary Material Supplementary material is available in the online version of this article at .

---

Resumen La modelación tridimensional del flujo de agua subterránea y del transporte de solutos en el acuífero de Pearl Harbor, en la parte sur de Oahu, Hawaii, muestra que el reajuste de la zona de transición agua dulce–agua salada, toma un largo tiempo a partir de cambios en el bombeo, irrigación o recarga en el sistema acuífero. Le toma alrededor de 50 años, a la zona de transición, moverse el 90% de la distancia hacia su nueva posición estacionaria. Además, el estimativo de Ghyben–Herzberg, sobre la profundidad de la interfase agua dulce–agua salada, se encuentra entre el 10 y el 50% en los contornos simulados de concentración de agua salada, de una manera compleja, durante 100 años de la historia de bombeo del acuífero. Por tanto, no es este un buen predictor de la profundidad del agua potable. Se utilizó una recarga pre – desarrollo, para simular la configuración del lente de agua dulce en 1880. Fueron utilizadas las distribuciones históricas del bombeo y de la recarga y se simularon el tamaño y posición resultantes del lente de agua dulce hasta 1980. Esas simulaciones muestran que la zona de transición se movió tierra adentro y hacia arriba, durante el periodo que se simuló.Los anteriores modelos de flujo para agua subterránea en Oahu, han sido limitados a modelos areales, que simulan una interfase abrupta entre agua dulce y agua salada, o bien han sido modelos de transporte de solutos que simulan una sección vertical del acuífero. El modelo presente está basado en el programa de computador del US Geological Survey (3D SUTRA), para transporte de solutos en tres dimensiones. Mediante el uso de varias herramientas nuevas para pre – procesamiento y post – procesamiento de las entradas y resultados del modelo, se ha permitido una construcción fácil del mismo y una visualización sin precedentes del lente de agua dulce y de la zona de transición subyacente en el acuífero más desarrollado de Hawaii.

---

Résumé La modélisation tridimensionnelle de lécoulement et du transport dans la partie sud de laquifère Oahu-Hawai montre que le temps de réajustement de la zone de transition entre leau douce et leau salée est assez long et dépend de la variation des pompages et des irrigations, ainsi que de la recharge du système aquifère. Il sont nécessaires 50 ans pour que la zone de transition parcoure 90% de la distance qui la sépare de sa nouvelle position. La profondeur du biseau estimée par le schéma Ghyben–Herzberg se trouve entre les contours de 10 et 50% de la concentration de leau salée. Ce résultat a été obtenu après la simulation de lhistoire du pompage de laquifère pendant une période de 100 ans. Donc le schéma Ghyben–Herzberg conduit aux valeurs erronées de la profondeur de leau potable. La valeur de la recharge davant lexploitation de laquifère a été utilisée pour simuler la configuration des lentilles deau douce en 1880. En utilisant lhistoire du pompage et la distribution de la recharge ont on a simulé les dimensions et le positions des lentilles deau douce jusqu› en 1980. Les simulation montrent que le mouvement de la zone de transition est ascendant et vers le continent.Les modèles antérieurs de la zone dOahu ont été des modèles locaux qui ont simulé une interface nette eau douce-eau salée ou des modèles de transport bidimensionnels, dans une coupe verticale. Le modèle actuel est basé sur le code 3D-SUTRA, réalisé par le Service Géologique des États-Unis. L› utilisation des différents techniques de traitement des données a permis une construction facile du modèle, ainsi qu› une visualisation sans précédent des lentilles deau douces et de la zone de transition sous-jacente dans le plus grand aquifère du Hawai.

     

Groundwater flow modeling in the Zhangye Basin,Northwestern China

The Zhangye basin is in the middle reaches of the Heihe River, northwestern China. Heavy abstraction of groundwater since the 1970s in the area is for agricultural, industrial and drinking water supplies and has led to a substantial decline in the potentiometric surface. A three-dimensional regional numerical groundwater flow model, calibrated under transient conditions, has been developed and used to predict the drawdown for the period from 2000 to 2030 under two different groundwater management scenarios.     

Relationships between hydraulic conductivity distribution and synsedimentary faults,Rharb-Mamora basin,Morocco; Hydrogeological,geostatistical and modeling approaches

This paper deals with the relationship between the physical characteristics of the Rharb-Mamora basin (Northern Morocco) and the distribution of hydraulic conductivity values in the Plio-Quaternary aquifer. A steady state groundwater model has been calibrated using these values. The distribution of hydraulic conductivity corresponds to a hydrodynamic partition with permeabilities of 2×10^–3 to 5×10^–2 ms^–1 (in the coastal part). A numerical modeling process shows a relative error map between calculated and measured piezometric surfaces, in order to facilitate calibration in steady state.The distribution of hydraulic conductivity values, and the lateral variation of sequences (deduced from engineering geology methods e.g. hydrogeological drilling data, seismic reflection interpretation, and geostatistical analysis), are controlled by synsedimentary fault activity, which deepened the Mio-Pliocene basement resulting in a thickening of overlying permeable rocks.

---

Resumen Esta nota trata la relación existente entre las características físicas y la distribución de la conductividad hidráulica de los acuíferos subyacentes de la cuenca del Rharb-Mamora (norte de Marruecos). Para la simulación del flujo del agua subterránea se ha establecido un modelo en régimen permanente. La distribución de la conductividad hidráulica corresponde a un reparto hidrodinámico cuyo grado de permeabilidad es generalmente alto: desde 2×10^–3 hasta 5×10^–2 ms^–1 (zona costera). Con el objetivo de facilitar la calibración del modelo en régimen permanente, el proceso numérico realizado mostró errores relativos entre la piezometría calculada y medida. La distribución de los valores de permeabilidad y la variación lateral de las unidades permeables ha sido deducida a partir de métodos de geología aplicada (análisis de los sondeos hidrogeológicos, interpretación de los perfiles de sísmica de reflexión y aplicación de métodos geoestadísticos), y están controlados por la actividad de las fallas sinsedimentarias que influyen sobre la topografía del basamento mio-plioceno y el grosor de las unidades acuíferas.

---

Résumé Cette note traite la relation entre les caractéristiques physique du bassin Rharb-Mamora (Nord du Maroc) et la répartition de la conductivité hydraulique du réseau aquifère. Le modèle de simulation des écoulements souterrains a été établi en régime permanent. La distribution de la conductivité hydraulique correspond à un découpage hydrodynamique dont le degré de perméabilité est généralement bon entre 2×10^–3 to 5×10^–2 ms^–1 (zone côtière). Afin de faciliter le calage en régime permanent, le code numérique établi permet de mettre en évidence les erreurs relatives entre la piézométrie calculée et mesurée.La distribution des valeurs de perméabilité et la variation latérale des unités perméables, déduites des méthodes de la géologie appliquée (analyse des forages hydrogéologiques, interprétation des profils sismique réflexion et analyse géostatistique), sont contrôlées par lactivité des failles synsédimentaires qui influencent lapprofondissement du substratum hydrogéologique mio-pliocène et lépaississement des corps perméables.

     

Numerical simulation of groundwater flow for a coastal plain in Japan: data collection and model calibration

Using a three-dimensional finite element model, this study characterizes groundwater flow in a costal plain of the Seto Inland Sea, Japan. The model characterization involved taking field data describing the aquifer system and translating this information into input variables that the model code uses to solve governing equations of flow. Geological geometry and the number of aquifers have been analyzed based on a large amount of geological, hydrogeological and topographical data. Results of study demonstrate a high correlation between the ground surface elevation and the groundwater level in the shallow coastal aquifer. For calibrating the numerical groundwater model, the groundwater flow was simulated in steady state. In addition, the groundwater level and trend in the transient state has also been elucidated. The numerical result provides excellent visual representations of groundwater flow, presenting resource managers and decision makers with a clear understanding of the nature of the types of groundwater flow pathways. Results build a base for further analysis under different future scenarios.     

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Identification and mapping of chromium (VI) plume in groundwater for remediation: A case study at Kanpur,Uttar Pradesh

With only twenty five percent population living in urban areas, India has cities amongst the biggest in the world. Urban growth in most of Indian cities is concomitant with rise in water demand for community, as well as, for industrial purposes. The complex situation resulting from indiscriminate disposal of waste and its severe impact on groundwater quality is set for continuous worsening mainly for want of sustained effort aimed at site-specific remediation. The study, a prerequisite for actual remediation in an industrial city of Kanpur, India, envisages detailed investigation about pollutant transport, evaluation of concept of Bio-remediation and a range of other options and finally full scale implementation of the best suited. Drilling of piezometers and resistivity survey indicates that the area is constituted of alluvial sands, gravels and their various admixtures. Chemical analysis of water samples collected from piezometers and hand pumps shows the presence of hexavalent chromium rich horizons at various depths. The alarming concentration of this carcinogenic heavy metal of the order of 16.3 mg/l against the permissible concentration (of 0.05 mg/l) for drinking water and high concentration within sediments of the area poses a major threat to the entire ecosystem. The projection of migration contaminant plume of hexavalent chromium as depicted in the paper is indicative of a concentrated extent of core zone existing in shallow alluvial aquifer, which may be targeted for interception by remedial measures. The present work, elaborating on the source, potential and monitoring the migration of the pollutant plume is the first field scale study of its kind in the country. The findings of these studies are of strong relevance in addressing the ground water pollution due to indiscriminate disposal practices of hazardous waste in areas located within the alluvial zones.     

Characterization of the shallow groundwater system in an alpine watershed: Handcart Gulch,Colorado, USA

Water-table elevation measurements and aquifer parameter estimates are rare in alpine settings because few wells exist in these environments. Alpine groundwater systems may be a primary source of recharge to regional groundwater flow systems. Handcart Gulch is an alpine watershed in Colorado, USA comprised of highly fractured Proterozoic metamorphic and igneous rocks with wells completed to various depths. Primary study objectives include determining hydrologic properties of shallow bedrock and surficial materials, developing a watershed water budget, and testing the consistency of measured hydrologic properties and water budget by constructing a simple model incorporating groundwater and surface water for water year 2005. Water enters the study area as precipitation and exits as discharge in the trunk stream or potential recharge for the deeper aquifer. Surficial infiltration rates ranged from 0.1–6.2×10^?5 m/s. Discharge was estimated at 1.28×10^?3 km³. Numerical modeling analysis of single-well aquifer tests predicted lower specific storage in crystalline bedrock than in ferricrete and colluvial material (6.7×10^?5–2.0×10^?3 l/m). Hydraulic conductivity in crystalline bedrock was significantly lower than in colluvial and alluvial material (4.3×10^?9–2.0×10^?4 m/s). Water budget results suggest that during normal precipitation and temperatures water is available to recharge the deeper groundwater flow system.     

Predicting long-term contamination potential of perched groundwater in a mine-waste heap using a random-walk method

Mine-waste heaps are potential long-term sources of contamination for surface-water courses and groundwater systems. Application of a novel physically based particle-tracking model to a mine-waste heap in northern England, UK, has enabled predictions to be made of the lifetime of contaminants leaching, revealing a pattern of source-mineral depletion. A mine-waste heap is conceptualised by a series of one-dimensional unsaturated “columns” in which active weathering of source minerals takes place. These columns drain into a saturated zone, through which the contaminants are transported to the heap discharge. Solute transport is simulated within the model by the random-walk method while reaction kinetics are incorporated to account for the timescales of source mineral depletion. Results reveal that the mine-waste heap is likely to remain polluting for several centuries, with the governing factor in the magnitude of pollution being the transport of the reactant, oxygen, to the source-mineral surfaces.     

Comparative analysis between analytical approximations and numerical solutions describing recession flow in unconfined hillslope aquifers

Recession flow of aquifers from a hillslope can be described by the non-linear Boussinesq equation. Under strong assumptions and for specific conceptual formulations, different authors derived analytical approximations or linearized versions to this partial differential equation. A comparative analysis between some analytical approximations of the Boussinesq equation and the numerical solution of the recession flow of an unconfined homogeneous aquifer (horizontal, inclined and concave aquifer floor) was carried out. The objective was to define the range where the analytical solutions approximate the numerical solution. The latter was considered in this study as the reference method, because it requires fewer assumptions. From the considered analytical approximations, exponential decay relationships were found to be mainly valid for fine domain materials when horizontal, mild slopes (less than 2%) and concave aquifer floors were considered, but failed to reproduce coarse aquifer numerical model outflows, in contrast to the quadratic decay relationship, which better reproduce outflows in such domains. On the basis of the comparative analysis, it has been found that recession flows obtained with the considered analytical approximations yield similar values only for certain ranges of aquifer properties and geometries.

---

Résumé L’écoulement de décrue dans les aquifères de versant peut être décrit par l’équation non-linéaire de Boussinesq. En considérant d’importantes hypothèses et des formulations conceptuelles spécifiques, différents auteurs ont dérivé des approximations analytiques ou des versions linéarisées de cette équation partielle différentielle. Une analyse comparative entre des approximations analytiques de l’équation de Boussinesq et une solution numérique de la décrue d’un aquifère homogène libre (limite inférieure de l’aquifère horizontale, inclinée et concave) a été effectuée. L’objectif était de définir la plage pour laquelle les solutions analytiques approchent la solution numérique. Cette dernière a été considérée dans cette étude comme la méthode de référence car elle demande moins d’hypothèses. Parmi les approximations analytiques considérées, les relations de décroissance exponentielle se sont avérées être principalement valides pour des matériaux fins lorsque des couches aquifères horizontales, concaves et de pentes douces (moins de 2%) étaient considérées; en revanche, celles-ci n’ont pas réussi à reproduire les écoulements sortants d’un modèle numérique pour un aquifère à matériaux grossiers, contrairement à la relation de décroissance quadratique qui reproduit mieux les écoulements sortants dans de tels milieux. Sur la base de l’analyse comparative, il s’est avéré que les écoulements en décrue obtenus avec les approximations analytiques considérées présentent des valeurs similaires seulement pour certains critères géométriques et propriétés de l’aquifère.

---

Resumen Puede describirse el flujo de recesión de acuíferos de una ladera mediante la ecuación no-lineal de Boussinesq. Bajo fuertes supuestos y para formulaciones conceptuales específicas, distintos autores han derivado aproximaciones analíticas o versiones lineales para esta ecuación parcial diferencial. Se realizó un análisis comparativo entre algunas aproximaciones analíticas de la ecuación de Boussinesq y la solución numérica del flujo de recesión de un acuífero homogéneo no confinado (piso de acuífero cóncavo, inclinado, y horizontal). El objetivo consistió en definir el rango donde las soluciones analíticas aproximan la solución numérica. En este estudio se consideró la solución numérica como el método de referencia debido a que requiere menores supuestos. De las aproximaciones analíticas consideradas se encontró que las relaciones de desintegración exponencial eran principalmente válidas para materiales de ámbitos finos donde se consideraron pisos de acuíferos horizontal, pendientes medianas (menores a 2%) y cóncavos pero fallaron en reproducir flujos de salida de modelos numéricos de ámbito grueso, en contraste con la relación de desintegración cuadrática la cual reproduce mejor los flujos de salida en tales ámbitos. En base al análisis comparativo se encontró que los flujos de recesión obtenidos con las aproximaciones analíticas consideradas aportaron valores similares solo para ciertos rangos de propiedades y geometrías de acuíferos.

     

Spatial characteristics of groundwater temperature in the Ishikari Lowland, Hokkaido, northern Japan: analytical and numerical applications

In porous sediments of the Ishikari Lowland, there is a gradual increase in the background geothermal gradient from the Ishikari River (3–4 °C 100 m^–1) to the southwest highland area (10 °C 100 m^–1). However, the geothermal gradient at shallow depths differs in detail from the background distribution. In spite of convective heat-flow loss generally associated with groundwater flow, heat flow remains high (100 mW m^–2) in the recharge area in the southwestern part of the Ishikari basin, which is part of an active geothermal field. In the northeastern part of the lowland, heat flow locally reaches 140 mW m^–2, probably due to upward water flow from the deep geothermal field. Between the two areas the heat flow is much lower. To examine the role of hydraulic flow in the distortion of the isotherms in this area, thermal gradient vs. temperature analyses were made, and they helped to define the major components of the groundwater-flow system of the region. Two-dimensional simulation modeling aided in understanding not only the cause of horizontal heat-flow variations in this field but also the contrast between thermal properties of shallow and deep groundwater reservoirs. Electronic Publication