Significance of enhanced infiltration due to groundwater extraction on the disappearance of a headwater lagoon system: Toluca Basin, Mexico

Mexico City relies significantly on groundwater resources drawn from the Sistema Lerma well field located in the Toluca Basin, Mexico. Enhanced infiltration caused by groundwater extraction is suspected to be both a prime factor in the disappearance of a lagoon system at the Toluca Basin and a potential risk to long-term groundwater quality. A combined approach of field investigation and numerical modeling was adopted to assess the groundwater-surface water interactions within the lagoon system. Potentiometric data indicate that current downward vertical hydraulic gradients below the lagoon and surrounding wetland area are extremely low suggesting very slow infiltration rates. Geochemical and isotope data from surface water and groundwater sampling also indicate that very little surface water infiltration has occurred. Numerical simulations demonstrate that enhanced surface water infiltration is unlikely to be the primary cause in the significant reduction in size of the lagoon system. Other factors such as modifications to the surface water drainage system and capture of spring flow from the surrounding mountainous regions are likely more significant. Simulations also suggest that contaminants originating in the lagoon system are currently entering nearby production wells although the total contaminant mass flux to the wells is still very low and significantly diluted.

---

Resumen La Ciudad de México depende significativamente de los recursos de aguas subterráneas obtenidos del campo de pozos Sistema Lerma ubicado en la Cuenca Toluca, oeste de la Cuenca de México. Se sospecha que la infiltración estimulada causada por extracción intensiva de aguas subterráneas es un factor principal en la desaparición de un sistema de lagunas en la cabeceras del Río Lerma en el centro de México y un riesgo potencial para la calidad de aguas subterráneas a largo plazo en la región. Se adoptó un enfoque combinado de investigación de campo y modelizado numérico para evaluar las interacciones de agua superficial-agua subterránea dentro del sistema de lagunas. Datos potenciométricos indican que los gradientes hidráulicos verticales descendentes actuales, ubicados por debajo de la laguna y vecina área de humedales, son extremadamente bajos sugiriendo ritmos de infiltración muy lentos. Datos geoquímicos e isotópicos provenientes del muestreo de agua superficial y subterránea también indican que ha ocurrido muy poca infiltración de agua superficial. Las simulaciones numéricas demuestran cuantitativamente que no es probable que la infiltración estimulada de agua superficial sea la causa principal en la reducción significativa del tamaño del sistema de lagunas en la cabecera. Otros factores tal como modificaciones al sistema superficial de drenaje de agua y captura de flujo de manantial en las regiones montañosas vecinas son probablemente más significativos. Las simulaciones también sugieren que los contaminantes que se originan en el sistema de lagunas están actualmente ingresando a pozos de producción cercanos aunque el flujo de masa contaminante total a los pozos es aún muy bajo y diluido significativamente en los volúmenes grandes de extracción. La magnitud de la amenaza a largo plazo a la calidad del agua subterránea en los pozos de producción Sistema Lerma proveniente de contaminantes infiltrados es una preocupación y amerita estudios futuros.

---

Résumé La Cité de Mexico dépend significativement des ressources en eau souterraine pompées au site de captages de Sistema Lerma, dans le Bassin de Toluca, à lOuest du Bassin de Mexico. Laugmentation de linfiltration, causée par lextraction extensive des eaux souterraines, est suspectée dêtre à la fois un facteur primaire de la disparition dun système de lagons en amont de la rivière Lerma au centre de Mexico, et un risque potentiel à long-terme pour la qualité de leau souterraine dans la région. Uen approche combinée détudes de terrain et de modélisation numérique a été adoptée pour estimer les interactions entre les eaux de surface et les eaux souterraines, au sein du système de lagons. Les données potentiométriques indiquent que les gradients hydrauliques verticaux descendants sont très bas au niveau des lagons, ce qui suppose un taux dinfiltration faible. Les données géochimiques et isotopiques provenant des eaux de surface et souterraines, indiquent également que peu deau de surface sinfiltrent. La simulation numérique démontre quantitativement que laugmentation de linfiltration des eaux de surface nest probablement pas la première cause de réduction des lagons. Dautres facteurs, tels que les modifications du drainage des eaux de surface, laménagement en captage des sources situées aux pieds des montagnes environnantes, sont probablement plus importants. Les simulations suggèrent également que les contaminants provenant des lagons se retrouvent également dans les captages les plus proches, bien que le flux total de contaminants entrant dans les captages reste faible. Les volumes dextraction étant très important ils diluent effectivement les contaminants. Les effets à long-terme de lextraction critique au niveau des forages de Sistema Lerma et de la dégradation de la qualité des eaux souterraines par les contaminants, nécessite une grande attention pour le futur.

     

A reactive flow model of the geothermal reservoir Waiwera, New Zealand

A numerical model of the complex interaction between fluid flow, heat transfer, and chemical reactions of the shallow, coastal, hydrothermal system of Waiwera (New Zealand) is used to test the proposed conceptual model of the field.Due to declining water levels, resulting from over-exploitation during the 1970s, the objective here was to set up a coupled fluid flow and heat transfer model to help in enabling a predictive and sustainable use of the resource. The presented conceptual model of the area is based on hydraulic, thermal, and chemical field observations, which date back as far as 1863. The numerical simulations were carried out with the reactive transport code SHEMAT.The inflow of geothermal water at the bottom of the reservoir prevents seawater from entering the aquifer. If seawater intrusion occurs, it is due to over-exploitation and happens in the upper parts of the aquifer. This is in contrast to common seawater-freshwater interfaces where seawater intrudes at the bottom. The numerical investigations emphasize that, after modifications of the production regime in the late 1980s, the system is recovering again.Additionally, geochemical calculations have been conducted to answer the question if dissolution or precipitation reactions might change the hydraulic properties of the geothermal aquifer. Mixing of fresh, geothermal, and seawater, although each of them in thermodynamic equilibrium with calcite, lead to calcite precipitation or dissolution. Nevertheless, the simulations show that the hydraulic properties of the aquifer are not significantly affected by dissolution or precipitation.

---

Resumen Se ha utilizado un modelo numérico de la interacción compleja entre el flujo de fluidos, la transferencia de calor y las reacciones químicas del sistema hidrotermal costero poco profundo de Waiwera (Nueva Zelanda) con el objeto de probar el modelo conceptual del área propuesto. Debido a los niveles de agua decrecientes, el objetivo de este trabajo es establecer un modelo pareado del flujo de aguas y de la transferencia de calor para colaborar en el uso predecible y sostenible del recurso. El modelo conceptual del área presentado está basado en observaciones de terreno hidráulicas, termales y químicas que datan de 1863. Las simulaciones se llevaron a cabo con el código de transporte reactivo SHEMAT. El influjo de agua geotermal al fondo del reservorio impide que el agua de mar entre en el acuífero. Si la intrusión de agua de mar se presenta, esto sucede a causa de la sobreexplotación y ocurre en las partes superiores del acuífero. Esto contrasta con las interfaces comunes de agua de mar-agua dulce donde la intrusión del agua de mar se da al fondo. Las investigaciones numéricas enfatizan que el sistema se está recuperando de nuevo luego de las modificaciones al regimen de producción a fines de los 80. Adicionalmente, los cálculos geoquímicos están dirigidos a contestar si las reacciones de disolución o precipitación pueden cambiar las propiedades hidráulicas del acuífero geotérmico. La mezcla de agua dulce, geotermal y de mar causa precipitación o disolución de calcita a pesar de que cada una de ellas se encuentra en equilibrio termodinámicao con la calcita. Sin embargo, las simulaciones muestran que las propiedades hidráulicas del acuífero no son afectadas siginificativamente por la disolución o precipitación.

---

Résumé Un modèle numérique de linteraction complexe entre lécoulement de leau, le transfert de chaleur et les réactions chimiques du système hydrothermal côtier peu profond de Waiwera (Nouvelle-Zélande) est utilisé afin de tester le modèle conceptuel. Lobjectif est de bâtir un modèle couplé découlement et de transfert de chaleur afin dencourager une approche prévisioniste et une utilisation durable de la ressource en réponse à des niveaux deau en déclin qui résultent dune surexploitation depuis les années 1970. Le modèle conceptuel présenté est basé sur des observations de terrain concernant lhydraulique, la température et la composition chimique de laquifère, des observations qui peuvent dater daussi loin que 1863. Les simulations numériques ont été réalisées avec le modèle de transport SHEMAT. Lentrée deau géothermale à la base du réservoir préviens lintrusion deau de mer dans laquifère. Si lintrusion deau de mer a lieu, elle est due à la surexploitation qui a lieu dans la partie supérieure de laquifère. Cette situation est lopposée des interfaces eau de mer - eau douce traditionnelles où leau de mer pénètre par la base. Les simulations numériques suggèrent que suite aux modifications du régime de production à la fin des années 1980, le système récupère. De plus, des calculs géochimiques ont été effectués afin de répondre à la question à savoir si les réactions de dissolution et de précipitation changent les propriétés hydrauliques de laquifère géothermal. Le mélange de leau douce, leau géothermale et leau de mer, même si elle sont tous en équilibre thermodynamique avec la calcite, provoque la dissolution ou précipitation de la calcite. Néanmoins, les simulations montrent que les propriétés hydrauliques de laquifère ne sont pas significativement affectées par la dissolution ou la précipitation.

     

Groundwater demand management at local scale in rural areas of India: a strategy to ensure water well sustainability based on aquifer diffusivity and community participation

Watershed development programmes provide an opportunity for sustainable management strategies, although currently, they remain largely supply-side mechanisms of water resources development. Hydrogeological conditions, community participation and status of groundwater usage are important in evolving strategies on demand-side groundwater management.Neemkheda aquifer is a typical low-storage, low-hydraulic conductivity aquifer from a watershed in the dryland regions of Madhya Pradesh State of central India. A shallow unconfined aquifer, it consists of an upper coarse, calcareous sandstone unit underlain by a fine-grained sandstone unit. A well commune of seven wells is poised to test the concept of joint groundwater management, wherein wells are mechanisms of tapping a common water source, the Neemkheda aquifer.The strategy for systematic groundwater management in the Neemkheda well commune is based upon the relationship between Transmissivity (T) and Storage coefficient (S), i.e. aquifer diffusivity, and its variation within the aquifer. Wells within a high diffusivity domain tend to dewater more quickly than wells within a low diffusivity domain. A well-use schedule during the dry season, based upon aquifer diffusivity forms the basis of the groundwater management concept. The distribution of local aquifer diffusivities governs the relationship between local and regional aquifer depletion times and forms the basis of the groundwater management exercise being proposed for the Neemkheda aquifer.

---

Resumen Los programas de desarrollo de una cuenca hídrica son una oportunidad para el uso de estrategias de gestión sostenible, aunque hoy en día estas siguen siendo principalmente mecanismos para el desarrollo de recursos hídricos con énfasis en la oferta. Las condiciones hidrogeológicas, la participación comunitaria y la condición de utilización del agua subterránea, son importantes en el desarrollo de estrategias para la gestión del agua subterránea, desde el punto de vista de la demanda.El acuífero de Neemkheda es un acuífero típico de almacenamiento bajo y conductividad hidráulica baja, perteneciente a una cuenca hídrica ubicada en las regiones secas del Estado de Madhya Pradesh, en la parte central de India. Un acuífero de tipo libre, somero, formado por una unidad superior de arenisca calcárea de grano grueso, subyacida por una unidad de arenisca de grano fino.Un campo de pozos comunitario compuesto por siete pozos, pone a prueba el concepto de gestión conjunta del agua subterránea, dentro del cual los pozos son mecanismos para usar una fuente de agua común: El acuífero de Neemkheda.La estrategia usada para la gestión sistemática del agua subterránea, en el campo de pozos comunitario de Neemkheda, se basa en la relación existente entre Transmisividad (T) y Coeficiente de Almacenamiento (S), es decir en la Difusividad del Acuífero y en su variación observada dentro del mismo acuífero. Los pozos cuya difusividad esta dentro de un rango alto, tienden a experimentar un descenso en su nivel mas rápidamente que aquellos pozos con una difusividad baja. El concepto de gestión de agua subterránea se basa en un programa especial de operación para cada pozo durante la estación seca, el cual a su vez se basa en la difusividad del acuífero. La distribución de difusividades locales del acuífero rige las relaciones entre las épocas de descensos de nivel en el acuífero a escala regional y local, y además constituyen la base del ejercicio de la gestión del agua subterránea que se esta proponiendo para el acuífero de Neemkheda.

---

Résumé Les programmes de développement des bassins versants offrent une occasion de développer des stratégies pour une gestion durable, bien qu à présent elles restent en grande mesure le terme source dans le mécanisme de développement de la ressource en eau. Les conditions hydrogéologiques, la participation de la communauté et lusage des eaux souterraines sont importants pour le terme demande dans les stratégies de gestion des eaux souterraines. Laquifère de Neemkheda presents des caractéristiques typiques pour un bassin versant dans les régions arides de létat Madhya Pradesh de lInde centrale en ce qui concerne les faible valeurs de la conductivité hydraulique et du coéfficient demmagasinement. Il sagit dun aquifère phrèatique dont la partie supérieure est constitué par des grès calcaire qui restent sur des grés plus fins. Afin dessayer le concept de gestion integrée, on a réalisé dans laquifère de Neemkheda un captage ayant sept forages qui forment la ressource en eau communale. La stratégie de la gestion systématique du captage est basée sur la relation entre la transmissivité (T), le coéfficient demmagasinement (S), donc la diffusivité hydraulique ainsi que sa variation spatiale. Les puits forés dans les zones à grande diffusivité tendent à sassècher plus vite que ceux creusés dans des zones à faible diffusivité. Le concept de la gestion des eaux souterraines est donc basé sur la distribution spatiale des diffusivités. À partir de ce concept on a développé un programme dexploitation des forages pendant les saisons sèches. La relation entre les temps d épuissment locaux et régionaux est determiné par la distribution spatiale de la diffusivitée et forme la base de lexercice sur la gestion des eaux proposée pour laquifère de Neemkheda.

     

Determination of aquifer roof extending under the sea from variable-density flow modelling of groundwater response to tidal loading: case study of the Jahe River Basin,Shandong Province,China

The main task of studies on salt-water intrusion into coastal confined aquifers is to predict the position of the fresh- salt-water interface, which can be determined from the length of the aquifer roof extending under the sea. Records of groundwater level affected by tides can be used to infer hydrological conditions and determine hydraulic parameters of an aquifer extending under the sea. In this paper, a three-dimensional, variable-density groundwater flow model has been developed to determine the equivalent roof length of an aquifer extending under the sea from the tidal-effected data of groundwater level in the Jahe River Basin, Shandong Province, China. The seaward boundary is obtained by converging hydraulic head fluctuations observed in drill holes with calculated values, and the aquifer parameters in the extending zone are estimated. The impacts of aquifer roof length and aquifer parameters on the fluctuation of tidal groundwater are studied. It is concluded that the length of the aquifer roof extending under the sea should correspond with certain aquifer parameters in the extrapolation zone. Therefore, the seaward boundary determined from tidal-effect information is the equivalent boundary in hydrodynamic characteristics rather than the true boundary of the confined aquifer

---

Resumen El principal objetivo de los estudios sobre intrusiones de agua salada en acuíferos costeros confinados es predecir la posición de la interfase agua dulce-agua salada, la cual puede determinarse a partir de la longitud del techo del acuífero que se extiende por debajo del mar. Los registros de niveles de agua subterránea afectados por las mareas puede utilizarse para inferir las condiciones hidrológicas y determinar los parámetros hidráulicos de un acuífero que se extiende por debajo del mar. En este artículo se ha desarrollado un modelo de flujo tri-dimensional de agua subterránea de densidad variable para determinar la longitud del techo equivalente de un acuífero que se extiende por debajo del mar a partir de datos, afectados por la marea, de niveles de agua subterránea en la Cuenca del Río Jahe, Provincia Shandong, China. El límite hacia el océano se obtiene por convergencia de fluctuaciones de presiones hidráulicas observadas en pozos con valores calculados, y se estiman los parámetros del acuífero en la zona extendida. Se estudian los impactos de la longitud del techo del acuífero y los parámetros del acuífero en la fluctuación del agua subterránea afectada por las mareas. Se concluye que la longitud del techo del acuífero que se extiende por debajo del mar debería corresponder con ciertos parámetros del acuífero en la zona de extrapolación. Por lo tanto, el límite hacia el océano determinado a partir de información de efectos de marea es el límite equivalente en características hidrodinámicas más que el límite real del acuífero confinado.

---

Résumé Les sujets principaux des études dinstrusion saline dans les aquifères confinés en zone côtière sont la prédiction de la position de linterface entre leau salée et leau fraîche, qui peut être déterminée à partir de lextention du toit de laquifère sous la mer. Les enregistrements des niveaux des eaux souterraines influencés par les marées peuvent être utilisés pour préciser les conditions hydrologiques et déterminer les paramètres hydrauliques dun aquifère possédant une extension sous la mer. Dans cet article, un modèle tridimensionnel comprenant des eaux souterraines de densité variable a été développé pour déterminer la longueur équivalente du toit dun aquifère qui sétend sous la mer à partir des données concernant les effets de marée sur les eaux souterraines dans le bassin de la rivière Jahe, dans la province de Shandong, Chine. La limite de salinité est déterminée en faisant converger les fluctuations des hauteurs piézométriques avec les valeurs calculées, et les paramètres de laquifère sont estimés dans la zone sétendant sous la mer. Lincidence de la longueur de laquifère sous la mer sur les fluctuations des niveaux est étudiée. On en conclut que la longueur du toit de laquifère sous la mer peut correspondre à certains aquifères paramètres dans la zone dextrapolation. Par conséquent, la limite de salinité déterminée à partir des effets de marée est léquivalent dune limite hydrodynamique plutôt que la véritable limite de laquifère.

     

Estimation of transport hydraulic parameters in loessic sediment, Argentina: Application of column tests

Dispersivity and porosity parameters needed to simulate flow and non-reactive transport were calculated and validated for loess-like sediments of the Pampean aquifer located in southeastern Buenos Aires province, Argentina, using infiltration experiments in laboratory columns filled with selected undisturbed sediment. To obtain values for these parameters, the experimental curves of arrival for chloride were visually adjusted to the curves obtained from numerical simulation using PHREEQC 2.0. Visual adjustment consisted of trying different parameter values until the most representative one was obtained. The best adjustment was obtained for dispersivity and porosity values of 0.01 m and 0.35, respectively. Since analytical solutions for the flow and transport equations are applicable to simplified cases such as this, an attempt was made to adjust the parameters automatically, using a model based on this type of solution. This automated analytical solution was developed using the Fortran programming language, and the automatic parameter calibration was achieved by minimizing the distance between the observed and calculated values of two parameters. The automated technique gave values of 0.0084 m and 0.39 for dispersivity and porosity, respectively. The sediments studied here constitute the uppermost unit of the Pampean aquifer. Given the relevance of this layer and the growing anthropic influence on it, parameter quantification is important for both laboratory tests and for hydrogeochemical simulations that investigate contaminant transport.

---

Resumen Se han obtenido los parámetros hidráulicos dispersividad y porosidad efectiva para sedimentos loessoides del acuífero denominado Pampeano cuya extensión abarca casi todo el sector central de la República Argentina. Los sedimentos estudiados en este trabajo se localizan en el sudeste de la provincia de Buenos Aires, Argentina. Los parámetros fueron calculados y validados mediante la obtención de curvas de arribo de cloruro en ensayos en columnas de laboratorio con el sedimento no perturbado, y su ajuste visual a curvas obtenidas en modelos numéricos con el programa PHREEQC 2.0. Los mejores ajustes se alcanzaron con valores de porosidad y dispersividad de 0.01 y 0.35 respectivamente. Asimismo, se implementó una solución analítica del problema en código FORTRAM para lograr el ajuste automático de estos dos parámetros mediante la minimización de las distancias en función de los dos parámetros. Los valores obtenidos fueron: dispersividad 0.0084 m y porosidad efectiva 0.39. Los resultados fueron validados mediante la realización de nuevos ensayos y ajustes numéricos. Los sedimentos aquí estudiados son de importancia por constituir la unidad mas superficial del sistema acuífero Pampeano, la cual es usualmente explotada por perforaciones domiciliarias. La cuantificación de sus parámetros hidráulicos es necesaria para ensayos de laboratorio y simulaciones hidroquímicas orientadas al estudio de la contaminación.

---

Résumé Les paramètres de dispersivité et de porosité nécessaires à la simulation des écoulements et des transports non-réactifs, ont été calculés et validés pour des sédiments loessiques de laquifère Pampéien dans le Sud-Est de la province de Buenos Aires, Argentine, au moyen dexpériences dinfiltration sur colonnes expérimentales remplies du sédiment non remanié. Pour obtenir les valeurs des paramètres, les courbes expérimentales de larrivée du chlore (élément conservatif) ont été visuellement calibrées sur les courbes obtenues par la simulation numérique sur Phreeqc 2.0.3 Lajustement visuel consistait à tester différents paramètres jusquà trouver le plus représentatif. Le meilleur ajustement était atteint avec des dispersivités et des porosités respectivement de lordre de 0.01 m et 0.35 m. Dés lors que la solution analytique de lécoulement et du transport est applicable dans des cas simplifiés comme celui-ci, il a également été opté dajuster automatiquement les paramètres, en utilisant un programme Fortran ; le calibrage automatique a été obtenu en minimisant la distance fonction des deux paramètres. Les valeurs obtenues ont été de 0.0084 m et 0.39 pour la dispersivité et la porosité. Le sédiment étudié ici constitue lunité supérieure de laquifère. Cette couche est primordiale notamment du point de vue de lactivité anthropique qui est menée dessus ; la quantification des paramètres de transport est importante pour les tests en laboratoire et les simulations hydrogéochimiques orientées vers les études de transport de soluté.

     

Evaluation of an intrinsic and a specific vulnerability assessment method in comparison with groundwater salinisation and nitrate contamination levels in two agricultural regions in the south of Portugal

The applicability of two vulnerability assessment methods in evaluating the impact of agricultural activities on groundwater quality, is tested in two areas in the south of Portugal with modest results. Intensive citri- and horticulture require large amounts of fertiliser and water supplied by irrigation, which induces groundwater salinisation and contamination by nitrates. The degree of contamination varies highly within and between the study areas and is related to hydrogeological factors as well as intensity of agricultural practices. Vulnerability mapping is performed with the intrinsic DRASTIC method and the specific Susceptibility Index (SI), which is an adaptation of DRASTIC. These methods can constitute useful groundwater management tools, for instance when designating new Nitrate Vulnerable Zones as defined in the European Directive 91/676/EEC. However, in the case of DRASTIC, little correspondence exists between the most vulnerable and the most contaminated areas. This is mainly a result of underestimating the dilution capacity and overemphasising the attenuating potential of the unsaturated zone and aquifer, as both chloride and nitrate prove to be very stable contaminants. By including a parameter for land use, SI manages to produce more reliable results, although in many areas the vulnerability is overestimated.

---

Resumen Se evalúa la aplicabilidad de dos métodos de estimación de vulnerabilidad en evaluar el impacto de actividades agrícolas en la calidad del agua subterránea para dos áreas en el sur de Portugal obteniendo resultados modestos. La horticultura y citricultura intensiva requiere grandes cantidades de fertilizantes y agua abastecida por riego, lo cual induce salinización de agua subterránea y contaminación por nitratos. El grado de contaminación varía fuertemente dentro y entre las áreas de estudio y se relaciona con factores hidrogeológicos así como con la intensidad de las prácticas agrícolas. El mapeo de vulnerabilidad se lleva a cabo con el método intrínsico DRASTIC y el Índice de Susceptibilidad específica (SI), el cual es una adaptación de DRASTIC. Estos métodos pueden constituir herramientas de manejo de aguas subterráneas útiles, por ejemplo al designar nuevas Zonas Vulnerables por Nitratos del modo que se definen en la Directiva Europea 91/676/EEC. Sin embargo, en el caso de DRASTIC, existen poca correspondencia entre las zonas más vulnerables y las áreas más contaminadas. Esto se debe principalmente a la subestimación de la capacidad de dilución y a al sobre énfasis del potencial de atenuación de la zona no saturada y el acuífero, ya que tanto cloruro como nitrato han probado ser contaminantes muy estables. Al incluir un parámetro del uso de la tierra, SI genera resultados más confiables, aunque en muchas áreas se sobrestima la vulnerabilidad.

---

Résumé Lapplication de deux méthodes de calcul de la vulnérabilité permettant dévaluer limpact des activités agricoles sur la qualité des eaux souterraines, est testée dans deux zones du Sud du Portugal, avec des résultats modestes. La citriculture et lhorticulture intensives nécessitent de grandes quantités e fertilisants et deau souterraine pour lirrigation, ce qui induit la salinisation et la contamination des eaux souterraines par les nitrates. Le degré de contamination varie grandement à lintérieur et entre les zones détudes, en fonction des facteurs hydrogéologiques et de lintensité des pratiques agricoles. La cartographie de la vulnérabilité est mise en oeuvre via la méthodologie DRASTIC et lIndex de Susceptibilité (SI) spécifique, qui est une adaptation de la méthode DRASTIC. Ces méthodes êuvent constituer des outils de management des eaux souterraines, par exemple lors de la désignation de nouvelles zones de vulnérabilité aux Nitrates selon la n Directive Européenne 91/676/EEC. Par ailleurs dans le cas de DRASTIC, de petites correspondances existent entre les zones les plus vulnérables et les plus contaminées. Ceci est principalement le résultat dune sous-estimation de la capacité de dilution et de la sur-accentuation du potentiel datténuation de la zone non-saturée de laquifère, car et le chlore et les nitrates sont des contaminants très stables. En incluant un paramètre dutilisation des sols, SI produit des résultats plus réalistes, bien que dans de nombreuses zones la vulnérabilité soit surestimée.

     

Combined interpretation of radar, hydraulic, and tracer data from a fractured-rock aquifer near Mirror Lake, New Hampshire, USA

An integrated interpretation of field experimental cross-hole radar, tracer, and hydraulic data demonstrates the value of combining time-lapse geophysical monitoring with conventional hydrologic measurements for improved characterization of a fractured-rock aquifer. Time-lapse difference-attenuation radar tomography was conducted during saline tracer experiments at the US Geological Survey Fractured Rock Hydrology Research Site near Mirror Lake, Grafton County, New Hampshire, USA. The presence of electrically conductive saline tracer effectively illuminates permeable fractures or pathways for geophysical imaging. The geophysical results guide the construction of three-dimensional numerical models of ground-water flow and solute transport. In an effort to explore alternative explanations for the tracer and tomographic data, a suite of conceptual models involving heterogeneous hydraulic conductivity fields and rate-limited mass transfer are considered. Calibration data include tracer concentrations, the arrival time of peak concentration at the outlet, and steady-state hydraulic head. Results from the coupled inversion procedure suggest that much of the tracer mass migrated outside the three tomographic image planes, and that solute is likely transported by two pathways through the system. This work provides basic and site-specific insights into the control of permeability heterogeneity on ground-water flow and solute transport in fractured rock.

---

Resumen Una interpretación integrada de radar experimental de campo transversal a pozos, trazadores, y datos hidráulicos demuestra el valor de combinar el monitoreo geofísico realizado en periodos de tiempo con mediciones hidrológicas convencionales en la caracterización mejorada de un acuífero rocoso fracturado. Se llevó a cabo tomografía de radar por periodos de tiempo y diferencia de atenuación durante un experimento con trazadores salinos en el sitio de investigación hidrológica de roca fracturada del Servicio Geológico de Estados Unidos cerca del Lago Espejo, Condado Grafton, New Hampshire, USA. La presencia del trazador salino eléctricamente conductivo refleja efectivamente fracturas permeables o trayectorias para imágenes geofísicas. Los resultados geofísicos orientan la construcción de modelos numéricos tri-dimensionales de flujo de agua subterránea y transporte de solutos. En un esfuerzo por explorar explicaciones alternativas para los datos tomográficos y trazadores se considera un conjunto de modelos conceptuales que involucran campos de conductividad hidráulica heterogéneos y transferencias de masa de ritmo limitado. La calibración de datos incluye concentraciones de trazadores, el tiempo de llegada de la concentración pico en la salida, y presión hidráulica en régimen permanente. Los resultados del procedimiento de acoplamiento invertido sugieren que mucho de la masa del trazador migró fuera de los tres planos de imagen tomográfica, y que el soluto es probablemente transportado por dos trayectorias a través del sistema. Este trabajo aporta ideas básicas y específicas del sitio en relación con el control de la heterogeneidad de permeabilidades en el flujo de agua subterránea y transporte de solutos en rocas fracturadas.

---

Résumé Une interprétation intégrée détudes de terrain (radar entre puits, traçages, données hydrauliques) démontre la valeur de la combinaison entre la géophysique des temps finis et les mesures hydrologiques conventionnelles pour une interprétation améliorée dun aquifère de roche fracturée. La tomographie au radar a été mise en uvre durant un traçage artificiel au sel au site de recherche sur lhydrologie des roches fracturées du Service Géologique des US, à proximité du Lac Mirror, Conté de Grafton, Nouvel Hampshire, USA. La présence du traceur électriquement conducteur met en relief, grâce à la géophysique, la présence de fractures ou découlements préférentiels. Les résultats de la géophysique ont permis la construction de modèle hydrogéologique tri-dimensionnel des écoulements et du transport de soluté. Dans loptique dexplorer des interprétations alternatives des données de traçage et de tomographie, différents modèles conceptuels sont utilisés concernant lhétérogénéité des conductivités hydrauliques et des taux limités de transferts de solutés. Les données du calibrage incluent les données de concentration du traceur, le temps darrivée du pic de restitution et les données piézométriques en régime permanent. Les résultats de la procédure dinversion couplée suggèrent quune quantité très importante du traceur migre au delà de la fenêtre de visualisation des tomographies, et que le soluté est transporté via deux voies découlement préférentiel. Ce travail apporte des connaissances de base et spécifiques au site concernant la distribution de la perméabilité dans laquifère et le transport de soluté dans les roches fracturées.

     

The origin of increased salinity in the Cambrian-Vendian aquifer system on the Kopli Peninsula,northern Estonia

Monitoring of the confined Cambrian-Vendian aquifer system utilised for industrial water supply at Kopli Peninsula in Tallinn over 24 years reveals remarkable changes in chemical composition of groundwater. A relatively fast 1.5 to 3.0-fold increase in TDS and in concentrations of major ions in abstracted groundwater is the consequence of heavy pumping. The main sources of dissolved load in Cambrian-Vendian groundwater are the leaching of host rock and the other geochemical processes that occur in the saturated zone. Underlying crystalline basement, which comprises saline groundwater in its upper weathered and fissured portion, and which is hydraulically connected with the overlying Cambrian-Vendian aquifer system, is the second important source of ions. The fractured basement and its clayey weathering crust host the Ca-Cl type groundwater, which is characterised by high TDS values (2–20 g/L). Intensive water abstraction accelerates the exchange of groundwaters and increases the area of influence of pumping. Chemical and isotopic studies of groundwater indicate an increasing contribution of old brackish water from the crystalline basement and rule out the potential implication of an intrusion of seawater into aquifer.

---

Resumen El origen del incremento en salinidad en un sistema de acuíferos Cámbrico-Vendiano en la Península Kopli, norte de EstonieMonitoreo a largo plazo de un sistema de acuíferos confinados, de edad Cámbrico-Vendiano, que se utiliza como fuente de abastecimiento industrial en la Península Kopli, al norte de Estonie, revela cambios notables en la composición química del agua subterránea. Un incremento de 1.5 a 3 veces en TDS y en concentraciones de iones mayores en agua subterránea explotada ha sido ocasionado por bombeo fuerte. Las fuentes principales de carga disuelta en el agua subterránea Cámbrico-Vendiano son la lixiviación de la roca encajonante y los procesos geoquímicos que ocurren en la zona saturada. Basamento cristalino subyacente, que aloja agua subterránea salada en la parte superior intemperizada y fisurada, y está conectado hidráulicamente con el sistema acuífero Cámbrico-Vendiano sobreyacente, es la segunda fuente importante de iones. El basamento fracturado y su corteza de intemperismo arcillosa alojan agua subterránea de tipo Ca-Cl la cual se caracteriza por valores altos de TDS (2–20 g/l). Debido a extracción intensiva se ha acelerado el intercambio de agua subterránea y se ha incrementado el área de influencia del bombeo. Los estudios químicos e isotópicos de agua subterránea indican una contribución creciente por filtración derivada del basamento cristalino. Es evidente una intrusión de agua salada hacia el sistema de acuíferos con implicaciones subsecuentes para la calidad del agua.

---

Résumé Lorigine de la salinité croissante dans le système aquifère du Cambrien-Vendien dans la péninsule de Kopli, nord de lEstonieLe suivi à long terme du système aquifère captif du Cambrien-Vendien utilisé pour lapprovisionnement deaux industrielles dans la Péninsule de Kopli, nord de lEstonie, révèle de remarquables changements dans la composition chimique des eaux souterraines. Une augmentation de facteur 1.5 à 3 de la TDS et des concentrations en ions majeurs dans leau souterraine est la conséquence de pompages intensifs. Les sources principales des charges dissoutes dans les eaux de laquifère du Cambrien-Vendien sont le lessivage des roches et dautres phénomènes géochimiques ayant lieu dans la zone saturée. Le soubassement rocheux cristallin, qui renferme des eaux souterraines salines dans sa partie supérieure altérée et fissurée, et est hydrauliquement connecté avec laquifère supérieur du Cambrien-Vendien, est la deuxième importante source dions. Le soubassement fracturé et le matériel argileux de laltération, renferme leau souterraine de type Ca-Cl, caractérisée par un haut TDS (2–20 g/l). A cause de la mobilisation intensive de leau les échanges deau souterraine est sont accélérés et la zone dinfluence des pompages augmentent. Les études chimiques et isotopiques indiquent une contribution croissante du drainage des eaux du soubassement cristallin. Lintrusion deaux salées de la mer dans le système aquifère nest pas un phénomène évident.

     

Characterising the main karst aquifers of the Alvand basin, northwest of Zagros, Iran, by a hydrogeochemical approach

The Alvand River basin, situated in the northwest of the Zagros mountain range, Iran, drains carbonate aquifers through some important karst springs. The physical, chemical and isotopic characteristics of spring water were studied for two years in order to assess the origin of groundwater and determine the factors driving the geochemical composition. Principal components analysis was used to identify the main factors controlling the water chemistry. Two groups of springs were identified: (1) low mineralisation, ion concentration, especially sulphate, low temperature, light isotope composition and high elevation of the recharge area, and (2) moderate to high mineralisation, especially sulphate, higher temperature, heavy isotope composition and low altitude of the recharge area. The main factors controlling the groundwater composition and its seasonal variations are the geology, because of the presence of evaporite formations, the elevation and the rate of karst development. In both groups, the carbonate chemistry is diagnostic of the effect of karst development. The supersaturation with respect to calcite indicates CO₂ degassing, occurring either inside the aquifer in open conduits, or at the outlet in reservoirs. The undersaturation with respect to calcite shows the existence of fast flow and short residence-time conditions inside the aquifer. A PCA analysis showed that, contrary to most developed karst systems, where dilution occurs during the wet season, leaching of the gypsum-bearing Gachsaran Formation by rainwater produced higher mineralisation.

---

Resumen La cuenca del río Alvand, situada en el noreste de la cadena montañosa Zagros drena acuíferos compuestos por carbonatos a través de manantiales kársticos. Durante dos años se estudiaron las características físicas, químicas e isotópicas del agua de manantial con el objetivo de determinar el origen del agua subterránea y de los factores que controlan su composición geoquímica. Se utililizó el análisis de componentes principales para identificar los factores principales que controlan la química del agua. Se identificaron dos grupos de mantaniales: (1) de baja mineralización, concentración de iones especialmente sulfatos, y baja temperatura, composición de isotopos livianos y elevación alta del área de recarga, y (2) mineralización moderada y alta, especialmente, sulfatos, y temperatura más alta, composición de isotopos pesados y altitud más baja del agua de recarga. Los factores principales que controlan la composición del agua subterránea y sus variaciones estacionales son la geología debido a la presencia de formaciones evaporitas, la elevación y la tasa a la que se desarrolla el karst. En ambos grupos la química de los carbonatos es un diagnóstico del efecto del desarollo de la estructura kárstica. La supersaturación con respecto a la calcita indica presencia de CO₂ que está en proceso de degasificación que ocurre dentro de acuífero en conductos abiertos o a la salida en reservorios. La subsaturación en lo referente a la calcita muestra la existencia de flujo rápido y de tiempos cortos de residencia al interior del acuífero. El análisis de componentes principales muestra que a diferencia de la mayor parte de sistemas kársticos en donde la dilución ocurre durante la temporada de lluvias, la lixiviación de la formación Gachsaran que contiene yeso a causa de lluvia produce mineralizaciones más altas.

---

Résumé Le bassin de la rivière Alvand, situé dans le nord-ouest de la chaîne de montagne Zagros, draine des aquifères carbonatés par le biais dun important réseau de sources karstiques. Les caractéristiques physiques, chimiques et isotopiques de leau provenant des sources ont été étudiées pendant deux ans afin de déterminer lorigine de leau souterraine et les facteurs responsables de sa composition chimique actuelle. Lanalyse des éléments chimiques majeurs a été utilisée pour identifier les processus qui contrôlent la composition chimique de leau. Deux groupes différents de sources ont été identifiés: (1) faible minéralisation et concentration en ions, spécialement en sulfates, faible température, composition isotopique légère et haute élévation de la zone de recharge, et (2) minéralisation moyenne à élevée, particulièrement en sulfates, température élevée, composition isotopique lourde et plus basse altitude de la zone de recharge. Les facteurs principaux qui contrôlent la composition et la variation saisonnière de leau souterraine sont la géologie, en raison de la présence dévaporites, lélévation et le taux de développement des karsts. Dans les deux groupes, la teneur en carbonates est un diagnostic de leffet du développement de la structure des karsts. La sursaturation par rapport à la calcite indique un dégazage en CO₂ qui a lieu soit à lintérieur de laquifère, dans les conduits ouverts, ou à sa sortie, dans les réservoirs. La sous-saturation par rapport à la calcite indique lexistence dun écoulement rapide et un faible temps de résidence dans laquifère. Une analyse des éléments chimiques majeurs a démontré que contrairement aux systèmes de karst les mieux développés où la dilution a lieu lors de la saison humide, le lessivage du gypse de la formation de Gachsaran par leau de pluie produit une minéralisation plus élevée.

     

A local-scale groundwater flow model for groundwater resources management in Dakhla Oasis,SW Egypt

Dakhla Oasis is located in the Western Desert of Egypt. Groundwater exploited from the Nubian Sandstone aquifer is the only available water resource in this area. This resource has been heavily exploited since 1960, which has led to a substantial decline in the potentiometric surface of the aquifer. A regional numerical groundwater flow model, calibrated under unsteady-state conditions, has been developed and used to investigate the hydrodynamic impacts of different groundwater management options on the potentiometry of the aquifer. To account for local details and to allow a precise analysis of pumping and the resulting drawdown in Dakhla Oasis, a local-scale model was developed by refining the grid cells in the calibrated regional model. The local-scale model gave a detailed picture about the expected drawdown due to the different groundwater management options in the next 100 years. The simulated results indicated that the planned increase in groundwater extraction will have a major impact on groundwater flow patterns in the whole area located southwest of Dakhla Oasis.

---

Resumen El oasis Dakhla está situado en el gran desierto occidental de Egipto. El agua subterránea que se extrae del acuífero de roca arenosa nubia es la única fuente de agua disponible en esta área. Este recurso ha sido explotado en medida considerable desde 1960. Este hecho he causado una disminución sustancial en la superficie potenciométrica del acuífero. Se ha desarollado un modelo de flujo de agua subterránea numérico regional calibrado bajo condiciones de estado variable. Este modelo se ha utilizado para investigar el impacto hidrodinámico que diferentes opciones de manejo de aguas subterráneas tienen sobre la potenciometría del acuífero. Se desarolló un modelo de escala local mediante el refinamiento de las células de malla en el modelo regional calibrado con el objecto de involucrar los detalles locales y de permitir un análisis preciso del bombeo y disminución del nivel de aqua resultante en el oasis Dakhla. El modelo de escala local proporcionó un imagen detallada de la disminución de nivel esperada según las diferentes opciones de manejo de agua subterranea en los próximos 100 años. Los resultados simulados indican que el incremento de extración planificado en la agua subterránea tendrá un gran impacto en el patrón de flujo de agua subterranea en toda el area ubicada al suroeste del oasis de Dakhla.

---

Résumé Loasis de Dakhla est située dans la partie occidentale du désert dÉgypte. La seule source deau disponible dans la région provient de lexploitation de laquifère gréseux de Nubian. Cette ressource a largement été exploitée depuis 1960, résultant en la diminution significative de la surface piézométrique de laquifère. Un modèle numérique découlement régional de leau souterraine, calibré en régime transitoire, a été développé et utilisé afin détudier les impacts hydrodynamiques de différents scénarios de gestion de laquifère sur la surface piézométrique. Afin dinclure les particularités locales et de réaliser une analyse précise du pompage et du rabattement qui y est associé dans lOasis de Dakhla, un modèle local découlement a été développé en raffinant le maillage dans le modèle régional préalablement calibré. Le modèle local découlement fournit une image détaillée du rabattement prévu pour les cent prochaines années selon les différents scénarios de gestion de laquifère. Les résultats obtenus par simulation numérique indiquent que laugmentation prévue dans lexploitation de leau souterraine aura un impact majeur sur les patrons découlement de leau souterraine dans toute la région située au sud-ouest de loasis de Dakhla.

     

Determination of hydraulic parameters of an unconfined alluvial aquifer by the floodwave-response technique

Groundwater constitutes a vital component of the water-resource system. Adequate knowledge of aquifer parameters is of utmost importance for proper groundwater management. In the present study, an attempt has been made to critically analyze the river–aquifer interaction and to explore the practicability and reliability of the floodwave-response technique in estimating the hydraulic parameters of the Konan aquifer, Japan. The analysis of the pertinent hydrologic and hydrogeologic data revealed a strong relationship between the Monobe River stage and the groundwater levels over a major portion of the Konan basin. Inverse modeling, using the floodwave-response model at four sites, yielded high values of hydraulic diffusivity that ranged from 16 to 194 m²/s. However, the overall average aquifer diffusivity was found to be reasonable (0.7–3.5 m²/s). All the selected flood events were not found suitable for determining reasonable values of hydraulic diffusivity at a given site. Using the optimal diffusivities at two sites and the aquifer transmissivity from the pumping tests at these sites, the storage coefficient values were estimated as those of truly confined aquifers. In contrast, reasonable values of storage coefficient were obtained based on the average aquifer diffusivity estimates. The sensitivity analysis of the floodwave-response model indicated that the optimal hydraulic diffusivity is very sensitive to the distance parameters of the model, which necessitate precise measurement of these parameters. It is concluded that a judicious application of the floodwave-response technique is essential, especially in the case of unconfined aquifers.

---

Resumen El agua subterránea constituye un componente vital de los sistemas de recursos hídricos. El conocimiento adecuado de los parámetros del acuífero es de suma importancia para el manejo apropiado del agua subterránea. En el presente estudio se ha realizado un intento para analizar críticamente la interacción río-acuífero y explorar la confiabilidad y utilidad de la técnica de respuesta de onda de crecida para estimar los parámetros hidráulicos del acuífero Konan, Japón. El análisis de datos hidrológicos e hidrogeológicos pertinentes muestra una relación fuerte entre el nivel del Río Monobe y los niveles de agua subterránea en gran parte de la cuenca Konan. El modelizado de inversión utilizando el modelo de respuesta de onda de crecida en cuatro sitios generó valores altos de difusividad hidráulica que varían de 16 a 194 m²/s. Sin embargo, se encontró que la difusividad promedio global del acuífero tenía valores razonables (0.7 a 3.5 m²/s). Todos los eventos de crecida seleccionados no fueron encontrados satisfactorios para determinar valores razonables de difusividad hidráulica en un sitio determinado. Utilizando las difusividades óptimas en dos sitios y la transmisividad del acuífero proveniente de pruebas de bombeo en esos sitios se estimaron valores del coeficiente de almacenamiento como los que se obtendrían para acuíferos completamente confinados. En contraste, se obtuvieron valores razonables del coeficiente de almacenamiento basados en los cálculos de difusividades promedio del acuífero. Los análisis de sensitividad del modelo de respuesta de onda de crecida indicaron que la difusividad óptima hidráulica es muy sensitiva a los parámetros de distancia del modelo por lo que se necesitan mediciones precisas de esos parámetros. Se concluye que es esencial una aplicación sensata de la técnica de respuesta de onda de crecida, especialmente en el caso de acuíferos no confinados.

---

Resumé Leau souterraine constitue une composante vitale des systèmes de ressource en eau. Une connaissance adéquate des paramètres des aquifères est importante pour une gestion correcte de ces systèmes. Dans cette étude, la démarche consiste à analyser linteraction rivière–aquifère de manière critique et dexplorer la technique permettant destimer les paramètres hydrauliques de laquifère Konan, Japon. Lanalyse pertinente des données hydrologiques et hydrogéologiques révèle une relation forte entre le niveau de la rivière Monobe et les niveaux piézométriques sur une grande partie du bassin de Konan. La modélisation inverse utilisant le modèle de réponse aux vagues de crue sur quatre sites fournit des valeurs de diffusivité comprises entre 16 et 194 m²/s. Néanmoins la valeur moyenne de diffusivité est apparue raisonnable (entre 0.7 et 3.5 m²/s). Tous les évènements de crue nont pas été retenus comme permettant de calculer avec une bonne fiabilité la diffusivité sur un site donné. En utilisant les données optimales de diffusivité provenant de deux sites, et la transmissivité fournie par les essais de pompage, le coefficient demmagasinement a été déterminé et correspond bien à un aquifère confiné. Par contraste, des valeurs raisonnables de coefficient demmagasinement ont été déterminées par estimation de la diffusivité moyenne de laquifère. Lanalyse de sensibilité du modèle de réponse aux vagues de crue indique que la diffusivité hydraulique optimale est très sensible à la distance utilisée dans le modèle, qui nécessite dés lors une mesure précise. En conclusion il est essentiel que ce type dapplication doit tre judicieusement mis en oeuvre, et spécialement dans le cadre daquifères captifs.

     

Arsenic in groundwater of the Bengal Basin, Bangladesh: Distribution, field relations, and hydrogeological setting

Arsenic contaminates groundwater across much of southern, central and eastern Bangladesh. Groundwater from the Holocene alluvium of the Ganges, Brahmaputra and Meghna Rivers locally exceeds 200 times the World Health Organisation (WHO) guideline value for drinking water of 10 µg/l of arsenic. Approximately 25% of wells in Bangladesh exceed the national standard of 50 µg/l, affecting at least 25 million people. Arsenic has entered the groundwater by reductive dissolution of ferric oxyhydroxides, to which arsenic was adsorbed during fluvial transport. Depth profiles of arsenic in pumped groundwater, porewater, and aquifer sediments show consistent trends. Elevated concentrations are associated with fine-sands and organic-rich sediments. Concentrations are low near the water table, rise to a maximum typically 20–40 m below ground, and fall to very low levels between about 100 and 200 m. Arsenic occurs mainly in groundwater of the valley-fill sequence deposited during the Holocene marine transgression. Groundwater from Pleistocene and older aquifers is largely free of arsenic. Arsenic concentrations in many shallow hand-tube wells are likely to increase over a period of years, and regular monitoring will be essential. Aquifers at more than 200 m below the floodplains offer good prospects for long-term arsenic-free water supplies, but may be limited by the threats of saline intrusion and downward leakage of arsenic.

---

Resumen El arsénico ha contaminado gran parte de las aguas subterráneas en el Sur, centro y Este de Bangla Desh. Su concentración en las aguas subterráneas del aluvial Holoceno de los ríos Ganges, Brahmaputra y Meghna supera localmente en un factor 200 el valor guía del arsénico en el agua potable, establecido por la Organización Mundial de la Salud (OMS) en 10 µg/L. Aproximadamente, el 25% de los pozos de Bangla Desh superan el estándar nacional de 50 µg/L, afectando al menos a 25 millones de personas. El arsénico ha llegado a las aguas subterráneas por la disolución reductora de hidróxidos férricos a los que se adsorbe durante el transporte fluvial. Los perfiles del arsénico en las aguas subterráneas bombeadas, agua de poro y sedimentos del acuífero muestran tendencias coherentes. Las concentraciones elevadas están asociadas a arenas finas y sedimentos ricos en materia orgánica. Las concentraciones de arsénico son bajas cerca del nivel freático, se incrementan hasta un máximo que se localiza generalmente a entre 20 y 40 m bajo la cota del terreno, y disminuyen a valores muy pequeños entre alrededor de 100 y 200 m. El arsénico se encuentra sobretodo en las aguas subterráneas existentes en la secuencia de sedimentación que tuvo lugar en el valle durante la transgresión marina del Holoceno. Las aguas subterráneas del Pleistoceno y acuíferos más antiguos están mayoritariamente libres de arsénico. Es probable que las concentraciones de arsénico aumenten en los próximos años en muchos pozos de tipo tubo perforados manualmente, por lo que será esencial efectuar un muestreo regular. Los acuíferos ubicados a más de 200 m bajo las llanuras de inundación ofrecen buenas perspectivas de abastecimiento a largo plazo sin problemas de arsénico, pero pueden estar limitados por las amenazas de la intrusión salina y de la precolación de arsénico desde niveles superiores.

---

Résumé Larsenic contamine les eaux souterraines dans la plus grande partie du sud, du centre et de lest du Bangladesh. Les eaux des nappes alluviales holocènes du Gange, du Brahmapoutre et de la Meghna dépassent localement 200 fois la valeur guide donnée par lOMS pour leau de boisson, fixée à 10 µg/l darsenic. Environ 25% des puits du Bangladesh dépassent la valeur standard nationale de 50 µg/l, affectant au moins 25 millions de personnes. Larsenic a été introduit dans les nappes par la dissolution par réduction doxy-hydroxydes ferriques sur lesquels larsenic était adsorbé au cours du transport fluvial. Des profils verticaux darsenic dans leau souterraine pompée, dans leau porale et dans les sédiments des aquifères montrent des tendances convergentes. Les concentrations élevées sont associées à des sédiments à sable fin et riches en matières organiques. Les concentrations sont faibles au voisinage de la surface de la nappe, atteignent un maximum typiquement entre 20 et 40 m sous le sol, puis tombent à des niveaux très bas entre 100 et 200 m. Larsenic est surtout présent dans les eaux souterraines de la séquence de remplissage de vallée déposée au cours de la transgression marine holocène. Les eaux souterraines des aquifères pléistocènes et plus anciens sont très largement dépourvus darsenic. Les concentrations en arsenic dans de nombreux puits creusés à la main doivent probablement augmenter au cours des prochaines années ; aussi un suivi régulier est essentiel. Les aquifères à plus de 200 m sous les plaines alluviales offrent de bonnes perspectives pour des alimentations en eau sans arsenic à long terme, mais ils peuvent être limités par les risques dintrusion saline et la drainance descendante de larsenic.

     

Recharge areas and geochemical evolution of groundwater in an alluvial aquifer system in the Sultanate of Oman

A regional hydrogeochemical model was developed to evaluate the geochemical evolution of different groundwaters in an alluvial aquifer system in the Interior of Oman. In combination with environmental isotopes the model is able to extract qualitative and quantitative information about recharge, groundwater flow paths and hydraulic connections between different aquifers. The main source of water to the alluvial aquifer along the flow paths of Wadi Abyadh and Wadi Muaydin in the piedmont is groundwater from the high-altitude areas of the Jabal Akhdar and local infiltration along the wadi channels. In contrast, the piedmont alluvial aquifer along Wadi Halfayn is primarily replenished by lateral recharge from the ophiolite foothills to the east besides smaller contributions from the Jabal Akhdar and local infiltration. Further down gradient in the Southern Alluvial Plain aquifer a significant source of recharge is direct infiltration of rain and surface runoff, originating from a moisture source that approaches Oman from the south. The model shows that the main geochemical evolution of the alluvial groundwaters occurs along the flow path from the piedmont to the Southern Alluvial Plain, where dedolomitization is responsible for the observed changes in the chemical and carbon isotope composition in these waters.

---

Resumen Un modelo hidrogeoquímico regional fue desarrollado, para evaluar la evolución geoquímica de diferentes aguas subterráneas en un sistema acuífero aluvial ubicado en el interior de Omán. El modelo, en combinación con isótopos ambientales, es capaz de obtener información cualitativa y cuantitativa sobre recarga, direcciones de flujo de agua subterránea y sobre las conexiones hidráulicas entre diferentes acuíferos. La fuente principal de agua para el acuífero aluvial en el piedemonte son las aguas subterráneas, a lo largo de las direcciones de flujo del Wadi Abyadh y del Wadi Muaydin, desde las áreas de gran altura de Jabal Akhdar y como infiltración local a lo largo de los canales del wadi. En contraste, el acuífero aluvial del piedemonte a lo largo del Wadi Halfayn, esta alimentado principalmente por recarga lateral desde las partes bajas de las montañas de ofiolitas hacia el Este, además por pequeñas contribuciones desde el Jabal Akhdar e infiltración local. Siguiendo gradiente abajo en el acuífero de la Llanura Aluvial del Sur, allí tanto la infiltración directa de lluvia, como la escorrentía superficial constituyen fuentes significativas de recarga, las cuales se originan en una fuente de humedad que llega a Omán proveniente del Sur. El modelo muestra, que la evolución geoquímica principal de las aguas subterráneas en el aluvión, sucede a lo largo de la dirección de flujo, que va desde el piedemonte hacia la Llanura Aluvial Sur, en donde el proceso de des-dolomitización es el responsable de los cambios observados en la composición química y también del contenido de isótopos de carbono en estas aguas.

---

Résumé On a développé un modèle régional hydrogéochimique pour évaluer lévolution géochimique des différentes eaux souterraines dans le système aquifère de lintérieur dOman. Avec la contribution des isotopes denvironnement, le modèle est capable de fournir des informations qualitatives et quantitatives sur la recharge, les directions découlement et sur les connections hydrauliques entre les différents aquifères du système. La source principale pour laquifère situé au long de la direction découlement des oueds de Abdyah et de Muaydin dans le piémont est leau souterraine provenant des zones de haute altitude de Jabal Akhbar ainsi que linfiltration locale au long des canaux de loued. Par contre, à part dune petite contribution provenant des infiltrations locaux et de Jabal Akhbar, l aquifère alluvial de Wadi Halfayan est en principal alimenté par la recharge latérale provenant des ophiolites situées dans la partie est. Les plus importantes sources de recharge dans laval de laquifère de la Plaine Alluviale sont linfiltration directe et le ruissellement provenant dune source dhumidité qui sapproche dOman par le sud. Le modèle montre que la principale évolution géochimique se produise au long de la direction découlement qui part de piémont de la Plaine Alluviale du Sud, où le procès de dédolomisation est responsable pour les changements observés dans la composition géochimique et en carbone des eaux.

     

Understanding and managing the stressed Mexico-USA transboundary Hueco bolson aquifer in the El Paso del Norte region as a complex system

As the latest drought has swept over the southwestern United States, most local and regional water providers have had to rely on ground water to meet their water demands as surface water supplies have diminished. In areas near major population centers this situation has continued to put strain and stress on already over-drafted regional aquifers. Metropolitan centers throughout the southwestern United States, as exemplified by Las Vegas, Nevada; Phoenix and Tucson, Arizona; and El Paso, Texas, have a history of over drafting their available ground water resources. How these stressed aquifers should be managed is an urgent issue. In El Paso, Texas long-term groundwater mining in the transboundary Hueco bolson aquifer has resulted in large water level drawdowns within the basin, particularly in several well fields. This situation has promoted the intrusion of brackish water into fresh water zones, and consequently reduced the volume of fresh water stored in the aquifer. Management of this aquifer is made more complex as the water resources of the region are shared across the United States/Mexico boundary and the transboundary region continues to grow rapidly. This paper describes the current situation within the Hueco bolson due to over drafting and continued reliance on the fresh water in storage and describes how understanding the dynamics of the aquifer as a complex system can offer managerial perspectives that can be used to establish operational procedures and programs for bringing drawdowns into equilibrium, increasing the longevity of the aquifer, gradually restoring water quality, and offering the hope of utilization of the resource in a sustainable manner.

---

Resumen A medida que la última sequía ha pasado ligeramente sobre el suroeste de Estados Unidos la mayoría de abastecedores locales y regionales de agua han tenido que depender del agua subterránea para satisfacer sus demandas de agua debido a que las fuentes de agua superficial han disminuido. En áreas cerca de centros poblados principales esta situación ha continuado ejerciendo presión en los ya sobre explotados acuíferos regionales. Los centros metropolitanos en todo el suroeste de Estados Unidos, por ejemplo Las Vegas, Nevada; Fénix y Tucson, Arizona; y El Paso, Texas tienen una historia de sobre explotar los recursos de agua subterránea disponibles.La gestión de estos acuíferos estresados constituye un problema urgente. En El Paso, Texas el minado de agua subterránea a largo plazo en el acuífero transnacional Hueco bolsón ha resultado en varios campos de pozos. Esta situación ha promovido la intrusión de agua salobre hacia zonas con agua fresca y reducido consecuentemente el volumen de agua fresca almacenado en el acuífero. El manejo de este acuífero se hace más complejo debido a que los recursos hídricos de la región son compartidos a través del límite Estados Unidos/México y la región transnacional continúa creciendo rápidamente. Este artículo describe la situación actual dentro del Hueco bolsón debido a la sobre explotación y continua dependencia en el agua fresca almacenada y describe como el entendimiento de la dinámica del acuífero como un sistema complejo puede ofrecer perspectivas de gestión que pueden utilizarse para establecer procedimientos operacionales y programas que conduzcan al equilibrio de los descensos, aumentando la longevidad del acuífero, gradualmente restaurando la calidad del agua, y ofreciendo la esperanza de utilizar el recurso de manera sostenible.

---

Résumé Alors que la dernière sécheresse a atteint toute la partie Sud-Ouest des Etats-Unis, la plus part des distributeurs deau, locaux et régionaux, doivent sapprovisionner en eau souterraine dés lors que les niveaux des eaux de surface ont diminués. Dans les zones proches des grands centres densément peuplés, cette situation a continuée a altérée des aquifères déjà surexploités. Les centres métropolitains du Sud-Ouest des Etats-Unis (Las Vegas dans le Nevada, Phoenix et Tucson dans lArizona, El Paso au Texas) sont connus historiquement pour surexploités leurs ressources en eaux souterraines. Il devient urgent de définir un mode de gestion de ses aquifères surexploités. A El Paso dans le Texas, le minage de laquifère transfrontalier engendre un rabattement important dans le bassin, et plus particulièrement autour des zones de captage. Cette situation a induit lintrusion deaux saumâtres dans des zones deaux douces, et a par conséquent réduit le volume deau douce stockée et exploitable dans laquifère. La gestion de cet aquifère est par ailleurs complexifiée par la situation transfrontalière, dautant que la zone frontalière connaît une croissance rapide. Cet article décrit la situation actuelle dans la vallée fermée du Hueco, due au sur-rabattement et à lexploitation continue des eaux douces, et décrit comment la compréhension de laquifère, vu comme un système complexe, peut offrir des perspectives de gestion amenant à des procédures opérationnelles et à des programmes à long terme permettant de ramener à léquilibre la baisse des niveaux deau. Ainsi la longévité de laquifère sera mieux garantie, la qualité de leau sera graduellement restaurée et lutilisation de leau respectera le souci de développement durable.

     

Assessing the extent of induced leakage to an urban aquifer using environmental tracers: an example from Bishkek, capital of Kyrgyzstan, Central Asia

A groundwater residence time study of the deep fluvioglacial aquifer supplying Bishkek, capital of Kyrgyzstan, has found evidence of deep infiltration of recent recharge both in the main periurban wellfield and below the city. Commonly-employed hydrochemical markers detected urban influence in the city-centre to depths of 65–100 m, but gave no indication of the important role of induced river/canal bed leakage, either upgradient in the periurban wellfield or within the city. This was revealed by O and H stable isotope measurements, which showed that local rainfall/snowfall play little part in the aquifer water balance. More remarkably, the universal detection of CFCs and SF₆, including in boreholes with 140–220 m deep upper screens, demonstrated that induced leakage of water just a few decades old had penetrated much deeper into the aquifer system than other hydrochemical markers indicated. A two-dimensional flow model set up to test whether such deep pumping-induced leakage could occur below the periurban wellfield confirmed its feasibility. The results imply vertical infiltration rates of 5–10 m/year and demonstrate that in this not-uncommon intergranular aquifer setting, deep boreholes with deep screen settings do not necessarily abstract old water. Hence, there are major implications for urban groundwater management and protection in such settings.

---

Resumen Un estudio sobre el tiempo de residencia del agua subterránea, en un acuífero fluvioglacial profundo que abastece a Bishkek, capital de Kyrgyzstan, ha encontrado evidencias de infiltración profunda a partir de recarga reciente, tanto en el campo de pozos de las afueras de la ciudad, como también por debajo de la ciudad. Los indicadores hidroquímicos comúnmente usados, detectaron influencia urbana en el centro de la ciudad hasta profundidades de 65–100 m, pero no dieron indicación del papel importante del goteo inducido en el lecho del canal/río, tanto aguas arriba en el campo de pozos de los suburbios, como también dentro de la ciudad. Esto se reveló mediante mediciones de isótopos estables de O y H, las cuales mostraron que las precipitaciones/nevadas de tipo local, tienen una pequeña contribución en el balance de agua del acuífero. De manera más notoria, la detección universal de CFCs y SF6, incluidos en las rejillas superiores de pozos profundos con 140–220 m de profundidad, demostraron que un goteo inducido de agua, de solo unas pocas décadas de antigüedad, ha penetrado mucho más profundo dentro del sistema acuífero, de lo que los otros marcadores hidroquímicos indicaban. Un modelo bidimensional de flujo, realizado para probar si el goteo inducido—bombeo profundo podía suceder debajo del campo de pozos suburbano, confirmó su factibilidad. Los resultados implican tasas de infiltración vertical de 5 a 10 m/año, y demuestran que los pozos profundos con tramos de rejilla también profundos, no necesariamente extraen aguas antiguas, en estos ambientes acuíferos con porosidad intergranular. Por lo tanto, hay implicaciones importantes para la gestión y protección en zonas urbanas, y para esos ambientes de las aguas subterráneas.

---

Résumé Lanalyse du temps de résidence des eaux souterraines dans un aquifère fluvial profond, dorigine glaciaire, qui alimente en eau Bishek, la capitale de Kyrgyzstan a mis en évidence une recharge récente tant dans le captage en eau que sous la ville même. Les traceurs chimiques ont indiqué linfluence urbaine dans le centre de la ville, à une profondeur de 65–100 m, mais ils nont pas mis en évidence le rôle important de lalimentation induite par la couche du fleuve, tant en amont, dans le captage que en aval, dans la ville. Ces effets ont été relevés par des isotopes stables, comme O et h qui ont montré linfluence récente des précipitations dans le bilan hydrique de laquifère. De plus, la détection universelle de CFC et de SF6, y compris dans les forages ayant une profondeur de 140–220 m, montre que lalimentation en eau induite de plusieurs décennies a pénétré à des profondeurs plus grandes que celles indiqués par des traceurs chimiques. Les possibilités dune drainance verticale induite par les pompages dans le puits profonds du captage ont été vérifiées par un modèle bidimensionnel découlement. Les résultats de la modélisation ont conduit à une infiltration verticale de 5–10 m/an et ont démontré aussi que leau extraite par des puits profonds dans un aquifère granulaire nest pas nécessairement dage ancienne. Il sagit donc des implications majeures concernant la gestion et la protection des captages urbaines en eaux souterraines

     

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.

     

A very large scale GIS-based groundwater flow model for the Nubian sandstone aquifer in Eastern Sahara (Egypt,northern Sudan and eastern Libya)

A three-dimensional GIS-based groundwater flow model for the Nubian Sandstone Aquifer in the eastern Sahara was developed and calibrated under steady-state and transient conditions. The model was used to simulate the response of the aquifer to climatic changes that occurred during the last 25,000 years. The simulation results indicated that the groundwater in this aquifer was formed by infiltration during the wet periods 20,000 and 5,000 years b.p. The recharge of groundwater due to regional groundwater flow from more humid areas in the south was excluded. It also indicates that the Nubian Aquifer System is a fossil aquifer, which had been in an unsteady state condition for the last 3,000 years.

---

Resumen Fue desarrollado un modelo de flujo de agua subterránea en tres dimensiones, basado en un SIG, para el Acuífero Arenisca Nubian en el Sahara Oriental, el cual fue calibrado para condiciones de estado estacionario y transitorio. El modelo se usó para simular la respuesta del acuífero a los cambios climáticos que ocurrieron durante los últimos 25000 años. Los resultados de esta simulación indicaron que el agua subterránea en este acuífero, se formó por infiltración, durante los períodos húmedos que hubo hace 20000 y 5000 años, antes del presente. Fue excluida la recarga del acuífero debida a un flujo regional de agua subterránea proveniente de áreas con un clima más húmedo en el sur. El modelo también muestra, que el Sistema Acuífero Nubian es un acuífero fósil, el cual ha permanecido en una condición de estado no estacionario, durante los últimos tres mil años.

---

Résumé Pour laquifère gréseux Nubien de Sahara -Est on a mis au points un modèle tridimensionnel, basé sur GIS. Le modèle a été calibré tant pour lécoulement stationnaire que pour lécoulement transitoire. On a simulé après la réponse de laquifère aux changements climatiques des derniers 25000 ans. Les résultats des simulations indiquent que la nappe a été rechargée par des infiltrations pendant une période humide qui sétend 5000 et 20000 ans, dès temps actuel. On na pas pris en compte la recharge de laquifère par la zone plus humide située dans sa partie sud. Le modèle indique aussi que leau de laquifère Nubien est une eau fossile qui a eu un écoulement transitoire pendant les derniers 3000 ans.

     

Numerical investigation of seawater intrusion at Gooburrum,Bundaberg, Queensland,Australia

Seawater intrusion in coastal agricultural areas due to groundwater abstraction is a major environmental problem along the northeastern coast of Australia. Management options are being explored using numerical modelling, however, questions remain concerning the appropriate level of sophistication in models, choice of seaward boundary conditions, and how to accommodate heterogeneity and data uncertainty. The choice of seaward boundary condition is important since it affects the amount of salt transported into the aquifers and forms the focus of the present study. The impact of this boundary condition is illustrated for the seawater-intrusion problem in the Gooburrum aquifers, which occur within Tertiary sedimentary strata. A two-dimensional variable-density groundwater and solute-transport model was constructed using the computer code 2DFEMFAT (Cheng et al. 1998). The code was tested against an experiment for a steady-state freshwater-saltwater interface and against the Elder (Elder 1967) free-convection problem. Numerical simulations show that the imposition of the commonly-used equivalent hydrostatic freshwater heads, combined with a constant salt concentration at the seaward boundary, results in overestimated seawater intrusion in the lower Gooburrum aquifer. Since the imposition of this boundary condition allows water flow across the boundary, which subsequently takes salt into the aquifer, a careful check is essential to estimate whether too much mass of salt is introduced.

---

Resumen La intrusión de agua salada en áreas agrícolas costeras a causa de la abstracción de aguas subterráneas es un problema medioambiental mayor a lo largo de la costa noreste de Australia. Se están explorando opciones de manejo mediante el uso de modelos numéricos, sin embargo, quedan pendientes cuestionamientos con respecto al nivel apropiado de sofisticación de los modelos, la elección de condiciones de límites hacia el mar y cómo acomodar la heterogeneidad y la incertidumbre de los datos. La elección de límites hacia el mar es importante ya que afecta la cantidad de sal transportada a los acuíferos y forma el centro de interés del presente estudio. Se ilustra el impacto de esta condición de límites para el problema de la intrusión de agua de mar en los acuíferos de Gooburrum, Queensland, Australia. Se construyó un modelo de transporte de solutos y de densidad variable de aguas subterráneas de dos dimensiones mediante el uso del código de computación 2DFEMFAT. El código se probó contra un experimento para una interfase de agua salada-agua dulce de estado constante y contra el problema de Elder de convexión libre. Las simulaciónes numéricas muestran que las cimas de agua dulce hidrostáticas equivalentes utlizadas comúnmente combinadas con una concentración de sal constante en el límite hacia el mar resultan en una sobreestimación de la intrusión de agua de mar en el acuífero bajo de Gooburrum. Puesto que la imposición de esta condición de límites permite flujo de agua a través de los límites, el cual subsecuentemente transporta sal al acuífero, una verificación cuidadosa es esencial para estimar si se introduce demasiada masa de sal.

---

Résumé Investigations numériques de lintrusion deau de mer à Gooburrum, Bundaberg, Queensland. Lintrusion deau salée due au captage deau souterraine dans les régions agricoles situées le long de la côte Nord-Est de lAustralie est un grave problème environnemental. Des solutions à ce problème sont présentement étudiées à laide de la modélisation numérique. Cependant des questions subsistent concernant le niveau de sophistication du modèle, le choix des limites le long de lintrusion salée et la manière daccommoder lhétérogénéité et lincertitude des données. Le choix de la limite le long de lintrusion deau salée est important puisque cette limite détermine la quantité de sel transporté dans laquifère et constitue le sujet de la présente étude. Leffet du choix de cette limite est bien illustré par le problème dintrusion deau salée dans les aquifères de Gooburrum, Queensland, Australie. Un modèle bidimensionnel découlement à densité variable et de transport de solutés a été conçu en utilisant le code 2DFEMFAT. Ce code a été testé avec une expérimentation ayant pour objet lécoulement en régime permanent dune interface eau salée/eau douce. Il a aussi été mis à lépreuve par le problème de convection libre de Elder. Les simulations numériques révèlent que limposition de la limite communément utilisée, quon appelle la «charge hydrostatique équivalente deau douce», conduit, dans le cas o lon a une concentration constante en sels le long de lintrusion deau salée, à une surestimation de lintrusion deau salée dans laquifère inférieur de Gooburrum. Puisque que cest limposition de cette limite qui permet lécoulement de leau au-delà de la frontière de laquifère, ce qui introduit subséquemment du sel dans laquifère, une vérification minutieuse doit tre effectuée afin de contrôler la quantité de sel introduit.