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

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

---

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

---

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

     

Deuterium composition and flow path analysis as additional calibration targets to calibrate groundwater flow simulation in a coastal wetlands system

A conceptual model of the Lake Warden coastal wetlands system, Western Australia, was developed using hydraulic, chemical and stable isotopic data, and formed the basis for a groundwater flow model using the finite element numerical code (FEFLOW). The system to be modeled is complex. The surface water and groundwater within the wetlands system show varying salinity and isotopic composition over short distances and time frames. As a first step, the flow model was calibrated to observed groundwater levels measured since 2001 for both steady state and transient stresses. Particle tracking analysis was conducted to test the source areas of water discharging to the lakes within the wetlands system. The analysis was able to delineate the connectivity between the lakes in the wetland and the flow path. Enrichment of isotopic concentration is evident along a NE–SW transect and the data set provides a means for calibrating a detailed transport model. The study incorporates the varying deuterium composition of the water bodies directly into a transport model and a good match between observed and simulated temporal variations along the transect indicates that the model closely simulated the dynamics of water exchange between the lakes and groundwater within the system.     

A semi-distributed groundwater recharge model for estimating water-table and water-balance variables

A semi-distributed groundwater recharge model is presented, which estimates water-table fluctuation and water-balance variables. The model is expressed by the water-balance concept linking atmospheric and hydrogeological parameters to different water uses (industrial, agricultural, domestic, etc). It was calibrated and validated using 5 years of data collected in the Dogo Plain in Japan. A 3-year dataset, from 2000 to 2002, was used in the calibration, while a 2-year dataset, from 2003 to 2004, was used for the validation. Calibration of the model was achieved by the shuffled complex evolution automatic optimization of model parameters to match simulated results with measured water-table depth. Square roots of relative error (R²) are 0.88 and 0.90 for calibration and validation processes, respectively. Monthly evolution of water storage change was then estimated and the water-table drawdown in different pumping scenarios was simulated. Finally, the groundwater-pumping amount planned by the government for future sustainable groundwater utilization was evaluated. The government-planned groundwater-pumping amount is feasible in most regions while the midstream region should be paid more attention. This study offers a scientific basis to control and prevent depletion of groundwater resources.     

Numerical studies on saltwater intrusion in a coastal aquifer in northwestern Germany

Saltwater intrusion in coastal aquifers depends on the distribution of hydraulic properties, on the climate, and on human interference such as land reclamation. In order to analyze the key processes that control saltwater intrusion, a hypothetical steady-state salt distribution in a representative cross-section perpendicular to the coastline was calculated using a two-dimensional density-dependent solute transport model. The effects of changes in groundwater recharge, lowering of drainage levels, and a rising sea level on the shape and position of the freshwater/saltwater interface were modeled in separate simulations. The results show that the exchange of groundwater and surface water in the marsh areas is one of the key processes influencing saltwater intrusion. A rising sea level causes rapid progression of saltwater intrusion, whereas the drainage network compensates changes in groundwater recharge. The time scale of changes resulting from altered boundary conditions is on the order of decades and centuries, suggesting that the present-day salt distribution does not reflect a steady-state of equilibrium.     

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

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

Development of a sharp interface model that simulates coastal aquifer flow with the coupled use of GIS

A simulation program, which works seamlessly with GIS and simulates flows in coastal aquifers, is presented in the present paper. The model is based on the Galerkin finite element discretization scheme and it simulates both steady and transient freshwater and saltwater flow, assuming that the two fluids are separated by a sharp interface. The model has been verified in simple cases where analytical solutions exist. The simulation program works as a tool of the GIS program, which is the main database that stores and manages all the necessary data. The combined use of the simulation and the GIS program forms an integrated management tool offering a simpler way of simulating and studying saline intrusion in coastal aquifers. Application of the model to the Yermasogia aquifer illustrates the coupled use of modeling and GIS techniques for the examination of regional coastal aquifer systems.

---

Resumen Se ha desarrollado un modelo casi tridimensional de elementos finitos para simular el flujo de agua dulce y salada, tanto en régimen estacionario como en transitorio, en sistemas acuíferos costeros, bajo la hipótesis de separación por medio de una interfaz abrupta. Las ecuaciones del modelo han sido discretizadas mediante un esquema de Galerkin de discretización en elementos finitos. El modelo ha sido verificado en casos sencillos para los que existe solución analítica. Todos los datos necesarios se introducen y gestionan con un Sistema de Información Geográfica [SIG] por ordenador. El programa de simulación forma parte del programa de SIG, constituyendo una herramienta integrada de gestión para estudiar la intrusión salina en acuíferos costeros. La aplicación del modelo al acuífero de Yermasogia ilustra el uso acoplado de las técnicas de modelación y de SIG con el fin de examinar sistemas acuíferos costeros a escala regional.

---

Résumé Pour étudier un système aquifère côtier, nous avons développé un modèle aux éléments finis en quasi 3-D qui simule les écoulements deau douce et deau salée en régime aussi bien permanent que transitoire. Les équations qui les régissent sont discrétisées par un schéma de discrétisation de Garlekin aux éléments finis. Le modèle a été vérifié dans des cas simples où il existe des solutions analytiques. Toutes les données nécessaires sont introduites et gérées grâce à un logiciel de gestion de SIG. Le programme de simulation est utilisé comme un outil du logiciel de SIG, constituant ainsi un outil de gestion intégrée dont le but est de simuler et détudier lintrusion saline dans les aquifères côtiers. Lapplication du modèle à laquifère de Yermasogia illustre lutilisation couplée de la modélisation et des techniques de SIG pour létude des systèmes aquifères côtiers régionaux.

     

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.     

Effects of intraborehole flow on groundwater age distribution

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

Numerical simulation of sea water intrusion near Beihai, China

 A leaky aquifer system occurs in the coastal plain near Beihai, China. Seawater intrusion into the confined aquifer took place along the northern coast. Chloride concentrations at some observation wells increased steadily from 1988 and were at their peak in 1993. A quasi-three-dimensional element model has been developed to simulate the spatial and temporal evolution of hydraulic heads and chloride concentrations of the groundwater near the northern coast. The simulation model was based on the transition zone approach, which requires simultaneous solution of the governing water flow and solute transport equations. An irregular grid of a quadrangle was used to discretize the flow domain. Various aquifer parameters were verified with the numerical model in order to obtain satisfactory matches between computed values and observed data from an investigation. Three pumpage schemes were designed to use the calibrated model for prediction of future changes in water levels and chloride concentrations in groundwater in the study area. Results show that seawater intrusion would worsen in the confined aquifer if the current rates of groundwater pumpage continue. The alternative, to eliminate pumpage in the intruded area and to moderate pumpage rates from water supply wells far from the seashore, may limit seawater intrusion significantly and is considered attractive in the area. Received: 27 September 1999 · Accepted: 27 December 1999     

Numerical simulation of groundwater flowing to horizontal seepage wells under a river

A horizontal seepage well, consisting of an interconnected vertical well, galleries, chambers and small-diameter radiating bores, is used to acquire relatively clean water that has been filtered through natural alluvial deposits in a riverbed. It has wide application, especially in arid and semi-arid areas. The lack of calculation formulae or models for horizontal seepage wells, up until now, has resulted in several false applications. Based on the analysis of groundwater flow characteristics, it has been concluded that several flow regimes coexist and hydraulic head loss exists in the horizontal seepage well. To avoid the difficulty of confirming the flux or head distribution in such a complex system, the model boundary of the whole horizontal seepage well has been moved to that of just the vertical well, and the well-aquifer system was treated as a heterogeneous medium, where the horizontal seepage well itself is a highly permeability medium. A mathematical model has been developed, based on the coupled seepage-pipe flow, by the introduction of equivalent hydraulic conductivity according to different flow regimes. Then a three-dimensional finite difference numerical model, based on the mathematical model, was developed and applied to a horizontal seepage well in China. The numerical model verified the groundwater flow characteristics of the horizontal seepage well. An erratum to this article can be found at     

Groundwater flow modelling within a coastal alluvial plain setting using a high-resolution hydrofacies approach; Bells Creek plain, Australia

Ground penetrating radar (GPR) has proved to be an extremely useful geophysical tool, in conjunction with direct geological data, to develop a realistic, macroscopic, subjective-based conceptual model of aquifer architecture within a shallow coastal alluvial plain. Subsequent finite-difference groundwater modelling has not only enabled determination of the dominant groundwater flow paths for the plain, but has also quantified the effects of within-facies and between-facies sedimentary heterogeneity on those flow paths. The interconnection of narrow, unconfined alluvial channels and a broad, semi-confined alluvial delta is ensuring that most fresh groundwater that enters the plain in the form of precipitation or recharge from lateral bedrock hills, is discharged into the eastern coastal wetlands via that alluvial delta aquifer.     

Evaluation of groundwater quality in coastal areas: implications for sustainable agriculture

Seawater intrusion is a problem in the coastal areas of Korea. Most productive agricultural fields are in the western and southern coastal areas of the country where irrigation predominantly relies on groundwater. Seawater intrusion has affected agricultural productivity. To evaluate progressive encroachment of saline water, the Korean government established a seawater intrusion monitoring well network, especially in the western and southern part of the peninsula. Automatic water levels and EC monitoring and periodic chemical analysis of groundwater help track salinization. Salinization of fresh groundwater is highly associated with groundwater withdrawal. A large proportion of the groundwaters are classified as Na–Cl and Ca–Cl types. The Na–Cl types represent effects of seawater intrusion. The highest EC level was over 1.6 km inland and high Cl values were observed up to 1.2 km inland. Lower ratios of Na/Cl and SO₄/Cl than seawater values indicate the seawater encroachment. A linear relation between Na and Cl represents simple mixing of the fresh groundwater with the seawater. The saline Na–Cl typed groundwaters showed Br/Cl ratios similar to or less than seawater values. The Ca–HCO₃ type groundwaters had the highest Br/Cl ratios. Substantial proportions of the groundwaters showed potential for salinity and should be better managed for sustainable agriculture.     

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

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

---

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

---

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

---

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

     

Numerical simulation of groundwater and early warnings from the simulated dynamic evolution trend in the plain area of Shenyang,Liaoning Province (P.R. China)

Groundwater level is the most direct factor reflecting whether groundwater is in a virtuous cycle. It is the most important benchmark for deciding whether a balance can be struck between groundwater discharge and recharge and whether groundwater exploitation will trigger problems pertinent to environment, ecology and environmental geology. According to the borehole and long-term monitoring wells data in the plain area of Shenyang, a numerical groundwater model is established and used to identify and verify the hydrogeological parameters and balanced items of groundwater. Then the concept of red line levels, the control levels of groundwater is proposed, the dynamic evolution trend of groundwater under different scenarios is analyzed and predicted and groundwater alerts are given when groundwater tables are not between the lower limit and the upper limit. Results indicated: (1) The results of identification and verification period fitted well, and the calculation accuracy of balanced items was high; (2) with the implementation of shutting wells, groundwater levels in urban areas of Shenyang would exceed the upper limit water level after 2020 and incur some secondary disasters; (3) under the recommended scenario of water resources allocation, early-warnings for groundwater tables outside the range would occur in the year of 2020, 2023, 2025 respectively for successive wet, normal and dry years. It was imperative to reopen some groundwater sources and enhance real-time supervision and early-warning to prevent the occurrence of potential problems.     

Multi-aquifer pumping test to determine cutoff wall length for groundwater flow control during site excavation in Tokyo, Japan

Multi-aquifer pumping tests, using a multi-screen pumping well and multi-level piezometers, were carried out for groundwater flow control in a large-scale excavation site in Tokyo, Japan. The site was underlain by multi-layered confined aquifers. In the tests, pumping was carried out using a multi-aquifer pumping well in which a screen depth was chosen arbitrarily. Changes in groundwater pressure heads in each aquifer were measured at each screen position of the multi-aquifer pumping well. Hydraulic conductivity (K) and specific storage (S _s) of not only aquifers, but also for low permeability layers between the aquifers, were estimated using the Cooper-Jacob method, and calibrated by a finite element method (FEM) groundwater model. Four different cutoff wall lengths were assumed for final excavation depth, and correlations among wall length, pumping discharge and drawdown at the back of the cutoff wall were obtained from simulations using the K and S _s parameters in the FEM model. Then, the most suitable wall length was selected based on the simulated correlations considering environmental condition, construction period and cost of the cutoff wall.     

Sustainability analysis of groundwater resources in a coastal aquifer, Alabama

Fort Morgan Peninsula is an attached portion of a dynamic barrier complex in the northern Gulf of Mexico and is a large tourist area that brings in a significant amount of revenue for Alabama. Many of the hotels and tourist attractions depend on the groundwater as their water supply. The over-withdrawal of groundwater and saltwater intrustion will have a negative impact on the ecology, tourism and economy if groundwater resources are not properly monitored and managed. In this study a calibrated groundwater flow model was used to analyze the sustainability of groundwater resources at Fort Morgan Peninsula. Detailed flow budgets were prepared to check the various components of inflow and outflow under different water use and climatic conditions. The results indicated the locations where groundwater was over-pumped and subjected to saltwater intrusion, or will be subjected to saltwater intrusion under a range of projected water use and climatic conditions.     

Rates of salinization by free convection in high-permeability sediments: insights from numerical modeling and application to the Dutch coastal area

Numerical modeling and dimensional analysis is used to study the salinization of thick, high-permeability aquifers by free convection from a salt source at the surface. Current understanding of this process mainly concerns the initial stages of salinization only (boundary-layer development, break-up into fingers and initial phase of finger descent). In the modeling, special attention is paid to the role of two processes in the long-term salinization rate: (1) the progressive loss of salt from fingers by lateral diffusion, and (2) the coalescence of fingers during their descent. From the numerical simulations a relationship is derived that describes the development of the horizontally averaged salinity with depth and time as a function of permeability and initial-density contrast for aquifer Rayleigh numbers up to Ra =6,000. This relationship is consistent with and provides an extension to previous generalized relationships of the rate of finger descent. Its applicability to real-world aquifers (Ra >10⁵) that include complexities due to anisotropy, heterogeneity, and mechanical dispersion is discussed. Application to the Pleistocene coastal aquifer of the Netherlands (thickness 200 m, permeability 10^-11 m²) suggests that salinization of the aquifer during historic episodes of inundation by seawater occurred within decades.

---

Resumen Se utiliza modelos numéricos y un análisis dimensional para estudiar la salinización de acuíferos potentes de alta permeabilidad por convección libre a partir de una fuente salina superficial. El conocimiento acutal de este proceso se limita a las fases iniciales de la salinización (desarrollo de la capa de contorno, creación de digitaciones y fase inicial de la progresión de éstas). En la modelación, se presta atención especial al papel desempeñado por dos procesos de salinización a largo plazo: (1) la pérdida progresiva de sal por difusión lateral desde las digitaciones, y (2) la coalescencia de las digitaciones durante su avance. A partir de las simulaciones numéricas, se obtiene una relación que describe el desarrollo de la salinidad con la profundidad y el tiempo, promediada horizontalmente, el cual depende de la permeabilidad y del contraste inicial de densidad para números de Rayleigh inferiores a 6.000. Esta relación es coherente con índices previos generalizados del avance de las digitaciones, y representa una extensión a estos. Se discute su aplicabilidad a acuíferos reales (con números de Rayleigh superiores a 105), que tienen complejidades asociadas a la anisotropía, la heterogeneidad y la dispersión mecánica. La aplicación al acuífero costero Pleistoceno de los Países Bajos (20 m de potencia y 10-11 m2 de permeabilidad) sugiere que la salinización tuvo lugar en décadas, debido a episodios históricos de inundación por aguas marinas.

---

Résumé Une modélisation numérique et une analyse dimensionnelle ont été mises en oeuvre pour étudier la salinisation d'aquifères épais et à forte perméabilité par convection libre d'une source de sel en surface. La compréhension habituelle de ce processus concerne principalement les étapes initiales de la salinisation seule (développement d'une couche limite, partition en digitations et phase initiale de développement des digitation). Dans la modélisation, une attention particulière a été portée au rôle de deux processus du taux de salinisation à long terme: (1) la perte progressive de sel dans les digitations par diffusion latérale et (2) la coalescence des digitations au cours de leur développement. À partir de simulations numériques, une relation a été obtenue qui permet de décrire l'extension de la salinité horizontalement en profondeur et au cours du temps en fonction de la perméabilité et du contraste initial de densité pour des nombres de Rayleigh de l'aquifère jusqu'à Ra =6,000. Cette relation est compatible avec ces résultats et fournit une extension des relations précédemment généralisées du taux de développement des digitations. On discute son applicabilité à des aquifères réels (Ra >105) incluant des complexités liées à l'anisotropie, l'hétérogénéité et la dispersion mécanique. L'application à l'aquifère côtier du Pléistocène des Pays-Bas (épaisseur environ 200 m, perméabilité environ10–11 m²) laisse penser que la salinisation de cet aquifère au cours d'épisodes historiques d'inondation par la mer s'est produit durant des décennies.

---

     

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

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