Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

WEKU – a GIS-Supported stochastic model of groundwater residence times in upper aquifers for the supraregional groundwater management

 The supraregional GIS-supported stochastical model, WEKU, for the determination of groundwater residence times in the upper aquifers of large groundwater provinces is presented. Using a two-dimensional analytical model of groundwater flow, groundwater residence times are determined within two extreme cases. In the first case, maximal groundwater residence times are calculated, representing the part of groundwater, that is drained by the main surface water of a groundwater catchment area. In the second case, minimal groundwater residence times for drainage into the nearest surface water are determined. Using explicit distribution functions of the input parameters, mean values as well as potential ranges of variations of the groundwater residence times are derived. The WEKU model has been used for the determination of groundwater residence times throughout Germany. The model results – mean values and deviations of the groundwater velocity and the maximal and minimal groundwater residence times in the upper aquifers – are presented by general maps and discussed in detail. It is shown that the groundwater residence times in the upper aquifer vary regionally, differentiated between less than 1 year and more than 2000 years. Using this information, the time scales can be specified, until measures to remediate polluted groundwater resources may lead to a substantial groundwater quality improvement in the different groundwater provinces of Germany. With respect to its supraregional scale of application, the WEKU model may serve as a useful tool for the supraregional groundwater management on a state, federal or international level. Received: 15 August 1995 · Accepted: 15 October 1995     

天津市地下水开采对地面沉降影响的多元回归分析

地下水的过量开采是天津市引起地面沉降的主要原因。因此天津市提出了“压缩地下水开采量”、“地下水人工回灌”、“调整地下水开采层次”等控制地面沉降的3大技术措施。经过多年的努力,控制地面沉降效果明显。如何解决地下水开发与控制地面沉降的关系,更好的贯彻这3大技术措施,是该文编写的初衷。即在开采同样地下水量的情况下,如何使地面沉降量最小;或在地面沉降量容许的情况下,如何开采最大量的地下水。压缩地下水开采量是治理地面沉降的根本措施,亦即如何压缩采水量或调整开采层次会达到最好效果。论文对天津市某区各个地下水开采层的多年累计开采量、累计沉降量进行数据统计分析,建立了该地区累计沉降量及各个地下水开采层的多元相关方程。在此基础上,分析了各个地下水开采层对地面沉降影响的相关程度。以此为该区控制地面沉降的提供依据。     

Use of unit response functions for management of regional multilayered aquifers: application to the North Aquitaine Tertiary system (France)

Although there are many numerical applications used in aquifer management, few methodologies seem to be sufficiently adapted to provide solutions for regional-scale problems. Their applications for multilayered aquifers are also too limited. A new groundwater management tool adapted to such large physical systems has been developed, using the unit response function (URF) approach. This tool is based on the link between a numerical groundwater flow model and a mathematical optimisation tool. Pre- and post processing are based on a geographical information system (GIS). The developed tool has been applied to a real management case—the regional Aquitaine multilayered aquifer (France). A representative hydrodynamic model was built using MODFLOW 2000 code. First tests show that the linearity condition needed for the application of URF methodology is respected despite the high complexity of the system. Two management scenarios were tested. Optimal solutions thus obtained allow for viable exploitation alternatives to be proposed, which integrate complex environmental constraints.

---

Résumé Malgré un nombre important d’applications existant dans le domaine de la gestion des aquifères, peu de méthodologies semblent adaptées pour proposer des solutions à des problèmes réels de dimension régionale. Leurs application pour des systèmes multicouches sont souvent trop limitées. Un outil de gestion adapté à ces systèmes de large taille a été développé, basé sur la théorie des Fonctions de Réponses Unitaires (FRU). Cet outil est basé sur le couplage d’un modèle numérique d’écoulement et d’un code d’optimisation mathématique. Les entrées et sorties du modèle sont réalisées à l’aide d’un Système d’Information Géographique (SIG). L’outil développé est appliqué à un cas de gestion réel—le système aquifère régional multicouche aquitain (France). Un modèle hydrodynamique représentatif a été construit, basé sur le code MODFLOW 2000. Les premiers tests montrent que les conditions de linéarité nécessaires à l’application de la méthode FRU sont respectées, malgré l’importante complexité du système. Deux scénarii de gestion sont testés. Les solutions optimales ainsi obtenues permettent de proposer des solutions d’exploitation alternatives viables, intégrant des contraintes environnementales complexes.

---

Resumen Aunque hay muchas aplicaciones numéricas usadas en la gestión de acuíferos, unas pocas metodologías parecen estar adaptadas suficientemente para dar las soluciones a los problemas de escala regional. También son limitadas sus aplicaciones para los acuíferos multicapa. Se ha desarrollado una nueva herramienta para el manejo del agua subterránea, adaptada a tales sistemas físicos grandes, usando el avance de la Función de Respuesta Unitaria (URF). Esta herramienta se basa en la combinación entre un modelo numérico de flujo de agua subterránea y una herramienta de optimización matemática. El pre-proceso y el post-proceso se basan en un Sistema de Información Geográfico (GIS). La herramienta desarrollada se ha aplicado a un caso real de manejo en el acuífero regional multicapa de Aquitania (Francia). Se construyó un modelo hidrodinámico representativo, que usó el código MODFLOW 2000. Las primeras pruebas muestran que la condición de linearidad, necesaria para la aplicación de la metodología URF, se respeta a pesar de la complejidad alta del sistema. Se probaron dos escenarios de manejo. Las soluciones óptimas así obtenidas permiten la proyección de alternativas viables de explotación, las cuales integran restricciones medioambientales complejas.

     

Contributions of boron isotopes to understanding the hydrogeochemistry of the coastal detritic aquifer of Castellón Plain,Spain

The Castellón Plain alluvial aquifer, Spain, is intensively exploited to meet the demand for agricultural irrigation and industrial water supply. The geochemistry of its groundwater shows complex salinization in the northern and southern parts of the aquifer, with significant pollution from human origin in the central portion. Boron content and B isotope geochemistry are useful for distinguishing between various sources of pollution and their relative importance in different parts of this aquifer. Boron concentrations in the groundwater vary between 0.01 and 0.85 mg/L. In the more saline groundwaters, found at the northern and southern ends of the study area, the presence of B is linked to inputs from seawater and water with a calcium-magnesium sulphate facies, which feed the aquifer and clearly influence the chemistry of its waters. Evidence of B adsorption processes in some samples is shown by the low B/Cl ratios and the high values of δ¹¹B. In the central portion of the aquifer, the high B/Cl ratios and the strongly negative δ¹¹B are related to pollution of human origin.     

Hydrogeology and hydrochemistry of a shallow alluvial aquifer,western Saudi Arabia

A hydrogeological and hydrochemical study was conducted on a shallow alluvial aquifer, Wadi Wajj, in western Saudi Arabia to assess the influence of protection measures on groundwater quality. The hydrochemistry was assessed up-gradient and down-gradient from potential contamination sources in the main city in dry and wet seasons prior to and after the installation of major drainage and wastewater facilities. Wadi Wajj is an unconfined aquifer where water is stored and transmitted through fractured and weathered bedrock and the overlying alluvial sediments. Natural recharge to the aquifer is about 5% of rainfall-runoff. Hydrochemistry of the aquifer shows temporal and seasonal changes as influenced by protection measures and rainfall runoff. Both groundwater and runoff showed similar chemical signature, which is mostly of chloride-sulfate-bicarbonate and sodium-calcium type. Groundwater downstream of the city, though of poorer quality than upstream, showed significant improvement after the installation of a concrete runoff tunnel and a wastewater treatment plant. Concentrations of many of the groundwater quality indicators (e.g., TDS, coliform bacteria, and nitrate) exceed US Environmental Protection Agency drinking-water standards. Heavy metal content is, however, within allowable limits by local and international standards. The chemical analyses also suggest the strong influence of stream runoff and sewage water on the groundwater quality.

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.

---

Résumé Une étude hydrogéologique et géochimique a été menée sur l’aquifère phréatique alluviale Wadi Wajj dans l’Ouest de l’Arabie Saoudite afin d’évaluer l’influence de mesures de protection sur la qualité de l’eau souterraine. L’hydrogéochimie de l’eau a été étudiée en amont et en aval de sources potentielles de contamination dans la ville principale, pendant les saisons sèches et humides, avant et après l’installation de réseaux majeurs de drainage et d’eaux usées. L’aquifère Wadi Wajj est libre, l’eau est stockée et s’écoule dans les sédiments alluviaux et dans le socle fracturé et altéré sous-jacent. La recharge naturelle de l’aquifère représente 5% des eaux de pluie et de ruissellement. L’hydrogéochimie de l’eau de l’aquifère montre que les changements saisonniers et temporaires sont influencés par les mesures de protection et par le ruissellement des eaux pluviales. L’eau souterraine et l’eau de ruissellement ont présenté le même faciès chimique, de type bicarbonatée-sulfatée-chlorée et calco-sodique. En aval de la ville, l’eau souterraine, bien que de moins bonne qualité qu’en amont, a présenté une nette amélioration après l’installation d’un système de récupération et d’une station de traitement des eaux usées. Les concentrations de plusieurs paramètres indicateurs de la qualité de l’eau (tels que la charge totale dissoute, les coliformes, et les nitrates) dépassent les normes de potabilité de l’eau de consommation de l’agence américaine de la protection de l’environnement. Les teneurs en métaux lourds en revanche n’excèdent pas les normes locales et internationales. Les analyses chimiques indiquent aussi l’influence importante de l’écoulement par ruissellement et des eaux usées sur la qualité de l’eau souterraine.

---

Resumen Un estudio hidrogeológico e hidroquímico fue hecho en un acuífero somero de Wadi Wajj, en Arabia Saudi oeste para evaluar la influecia de medidas de protección en la calidad del agua subterránea. La hidroquímica fue evaluada gradiente-arriba y gradiente-abajo de las fuentes potenciales de contaminación de la ciudad principal durante las estaciones seca y lluviosa, antes y después de la instalación de sistemas principales de drenaje y aguas servidas. Wadi Wajj es un acuífero no-confinado donde el agua es almacenada y transmitida a través de roca fracturada y meteorizada, y los sedimentos aluviales que le sobreyacen. La recarga natural del acuífero es de cerca del 5% de la precipitación-escorrentía. La hidroquímica del acuífero muestra cambios temporales y estacionales influenciados por las medidas de protección y la escorrentía de precipitación. Ambas, agua subterránea y escorrentía mostraron composición química similar, siendo mayoritariamente de tipos cloruro-sulfato-bicarbonato y sodio-calcio. El agua subterránea aguas arriba de la ciudad, aunque de calidad más pobre que aguas abajo, mostró significante mejoría después de la instalación de un tunel de concreto para escorrentía y una planta de tratamiento de aguas servidas. Las concentraciones de muchos de los indicadores de calidad de agua subterránea (e.g., STD, coliformes, y nitrato) exceden los estándares de la Agencia de Protección Ambiental USA para agua potable. El contenido de metales pesados está, sin embargo, dentro de los límites permisibles de los estándares locales e internacionales. Los análisis químicos también sugieren la fuerte influencia de la escorrentía y aguas residuales en la calidad del agua subterránea.

     

Occurrence of salt water above fresh water in dynamic equilibrium in a coastal groundwater flow system near De Panne, Belgium

A salt water lens is found above fresh water under the shore between Dunkerque (France) and Nieuwpoort (Belgium). This inverse density distribution is in a dynamic equilibrium. It develops due to the infiltration of salt water on the back shore during high tide. Under this salt water lens, water infiltrated in the adjacent dune area flows towards the sea and discharges at the seabed. This water quality distribution differs from the classic salt water wedge under fresh water described in the literature. Here, the evolution to this water quality distribution is simulated with a density dependent numerical model. A large tidal range, shore morphology and a permeable groundwater reservoir are the main conditions for the observed water quality distribution. By altering these conditions, intermediate water quality distributions between the classic salt water wedge and the one discussed here develop. Based on these simulations, it is expected that similar kinds of inverse density distribution could be present in a number of coastal areas, which have tides, a gently sloping shore and a permeable substratum.     

Hydrochemistry and isotope geochemistry as tools for groundwater hydrodynamic investigation in multilayer aquifers: a case study from Lomellina, Po plain, South-Western Lombardy, Italy

A multicriteria approach in studying hydrodynamics of a multilayer aquifer system has been used in the Lomellina region (Northern Italy). It involves the reconstruction of the hydrogeological framework coupled to the definition of the hydrochemical and isotopic features of the aquifers. A shallow phreatic aquifer, reaching depths of about 60–80 m from the surface, and deeper aquifers containing confined groundwater, were distinguished. Groundwater generally shows mineralisation decreasing with depth; dissolved ions depict calcium-bicarbonate hydrochemical facies and stable isotopes define the recharge mechanisms, the origin of groundwater, and the hydraulic confinement of deep aquifers. The phreatic aquifer is fed by local infiltration and by streams and irrigation channels. Tritium and Carbon-14 groundwater dating indicate long residence times (on the order of thousands of years) for confined aquifers. The confined aquifers show essentially passive hydrodynamic conditions and maintain a higher piezometric level than the phreatic aquifer. This inhibits the possibility of recent water penetrating far below the surface. The hydrogeological setting of the Lomellina region displays features which are common to other sectors of the Po plain. As a consequence, the results of this study, although conducted on a restricted area, are highly illustrative of groundwater hydrodynamics in large sedimentary aquifers.     

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     

Steady, annual, and monthly recharge implied by deep unconfined aquifer flow

We consider the response of a deep unconfined horizontal aquifer to steady, annual, and monthly recharge. A groundwater divide and a zero head reservoir constrain the aquifer, so that sinusoidal monthly and aperiodic annual recharge fluctuations create transient specific discharge near the reservoir and an unsteady water table elevation inland. One existing and two new long-term data sets from the Plymouth-Carver Aquifer in southeastern Massachusetts calibrate and confirm hydraulic properties in a set of analytical models. [Geohydrology and simulated groundwater flow, 1992] data and a new power law for tritiugenic helium to tritium ratios calibrate the steady recharge that drives the classical parabolic model of steady hydraulics [Applied Hydrogeology, 2001]. Observed water table and gradient fluctuations calibrate the transient recharge models. In the latter regard, monitoring wells within 1 km of Buttermilk Bay exhibit appreciable specific discharge and reduced water table fluctuations. We apply [Trans Am Geophys Union 32(1951)238] periodic model to the monthly hydraulics and a recharge convolution integral [J Hydrol 126(1991)315] to annual flow. An infiltration fraction of 0.79 and a consumptive use coefficient of 1.08×10⁻⁸ m/s °C relate recharge to precipitation and daylight weighted temperature across all three time scales. Errors associated with this recharge relation decrease with increasing time scale.