Numerical analysis of the hydrogeologic controls in a layered coastal aquifer system, Oahu, Hawaii, USA

 The coastal aquifer system of southern Oahu, Hawaii, USA, consists of highly permeable volcanic aquifers overlain by weathered volcanic rocks and interbedded marine and terrestrial sediments of both high and low permeability. The weathered volcanic rocks and sediments are collectively known as caprock, because they impede the free discharge of groundwater from the underlying volcanic aquifers. A cross-sectional groundwater flow and transport model was used to evaluate the hydrogeologic controls on the regional flow system in southwestern Oahu. Controls considered were: (a) overall caprock hydraulic conductivity; and (b) stratigraphic variations of hydraulic conductivity in the caprock. Within the caprock, variations in hydraulic conductivity, caused by stratigraphy or discontinuities of the stratigraphic units, are a major control on the direction of groundwater flow and the distribution of water levels and salinity. Results of cross-sectional modeling confirm the general groundwater flow pattern that would be expected in a layered coastal system. Groundwater flow is: (a) predominantly upward in the low-permeability sedimentary units; and (b) predominantly horizontal in the high-permeability sedimentary units. Received, October 1996 Revised, August 1997 Accepted, September 1997     

Free convective controls on sequestration of salts into low-permeability strata: insights from sand tank laboratory experiments and numerical modelling

Using sand tank experiments and numerical models, local-scale solute-transport processes associated with free convection in both the region surrounding as well as within discrete low-permeability strata are explored. Different permeability geometries and contrasts between high- and low-permeability regions are tested. Results show that two free convective processes occur at different spatial and temporal scales. In the high-permeability region, salinisation was rapid and occurred predominantly by free convective flow around the low-permeability blocks (interlayer convection). A free convection flow field also became concurrently established within the low-permeability lenses (intralayer convection). It was found that upward vertical flow created by the large-scale interlayer free-convective flow field in the high-permeability region retards salinisation of the lenses as buoyant freshwater displacements oppose the downward penetration of dissolved salts. Salinisation of the low-permeability structures eventually takes place from below as saltwater is dragged upwards. This bottom up convective salinisation process of low permeability strata has not been reported in previous literature. These results demonstrate that variable-density sequestration of solutes driven by a source resident above the low-permeability layer is a complicated function of the geometry of the permeability distribution and the permeability contrast between low- and high-permeability regions.     

Estimating aquifer thickness using multiple pumping tests

A method to estimate aquifer thickness and hydraulic conductivity has been developed, consisting of multiple pumping tests. The method requires short-duration pumping cycles on an unconfined aquifer with significant seasonal water-table fluctuations. The interpretation of several pumping tests at a site in India under various initial conditions provides information on the change in hydrodynamic parameters in relation to the initial water-table level. The transmissivity linearly decreases compared with the initial water level, suggesting a homogeneous distribution of hydraulic conductivity with depth. The hydraulic conductivity is estimated from the slope of this linear relationship. The extrapolation of the relationship between transmissivity and water level provides an estimate of the aquifer thickness that is in good agreement with geophysical investigations. The hydraulically active part of the aquifer is located in both the shallow weathered and the underlying densely fractured zones of the crystalline basement. However, no significant relationship is found between the aquifer storage coefficient and initial water level. This new method contributes to filling the methodological gap between single pumping tests and hydraulic tomography, in providing information on the variation of the global transmissivity according to depth. It can be applied to any unconfined aquifer experiencing large seasonal water-table fluctuations and short pumping cycles.     

The effects of geological heterogeneities and piezometric fluctuations on groundwater flow and chemistry in a hard-rock aquifer, southern India

Crystalline aquifers of semi-arid southern India represent a vital water resource for farming communities. A field study is described that characterizes the hydrodynamic functioning of intensively exploited crystalline aquifers at local scale based on detailed well monitoring during one hydrological year. The main results show large water-table fluctuations caused by monsoon recharge and pumping, high spatial variability in well discharges, and a decrease of well yields as the water table decreases. Groundwater chemistry is also spatially variable with the existence of aquifer compartments within which mixing occurs. The observed variability and compartmentalization is explained by geological heterogeneities which play a major role in controlling groundwater flow and connectivity in the aquifer. The position of the water table within the fracture network will determine the degree of connectivity between aquifer compartments and well discharge. The presented aquifer conceptual model suggests several consequences: (1) over-exploitation leads to a drop in well discharge, (2) intensive pumping may contribute to the hydraulic containment of contaminants, (3) groundwater quality is highly variable even at local scale, (4) geological discontinuities may be used to assist in the location of drinking-supply wells, (5) modeling should integrate threshold effects due to water-table fluctuations.     

On the advection of turdidity in the saint lawrence middle estuary

The suspension transport away from the extensive turbidity zone of the St. Lawrence estuary is largely determined by the channel topography. The suspended sediments are advected downstream by a 40 km long turbid plume which takes its source in a turbidity maximum at the head of the estuary and flows downstream partly confined by the South Channel. During the ebbing phase of tidal cycles, the turbid waters of the plume are forced downstream through narrow converging sections of the channel, and slowed down through more opened diverging regions, particularly down-stream of the St. Roch Traverse. These, large fluctuations in stratification modulate the vertical transport of suspended material from the bottom to the surface layer. Midway down the estuary, dispersion of the plume occurs along a frontal zone which seasonally migrates 30 km or more in response to changes in fresh water discharge. The plume is reinforced and the turbidity gradient is intensified by local injections of inshore waters from Ste. Anne Bay, a subtidal platform highly enriched in suspended material by intertidal exchanges with large mudllats. Lateral erosion of the plume and cross-channel transport of suspended matter from the South into the North Channel is made possible by large horizontal shears developing in the central part of the middle estuary during the early flood. These are created by a one-hour tidal phase difference between the North and the South Channel.     

Hydrogeologic framework of Pico Island, Azores, Portugal

Pico, the youngest island of the Azores archipelago, is composed of basaltic volcanic deposits less than 300,000 years old. The principal aquifer system consists mainly of recent lava flows that are very permeable and whose head is influenced by tidal fluctuations. Groundwater abstraction is almost entirely by drilled wells. The hydraulic gradient is very low, about 10^–4, which agrees with observations made on similar volcanic islands. Groundwater also occurs in perched-water bodies, but the spring discharge from them is very low, about 10^–3 L/s. The transmissivity of the volcanic rocks ranges from 9.44×10^–3 to 3.05×10^–1 m²/s, indicating the heterogeneity of the aquifers. The hydraulic diffusivity, estimated from observations of the effects of tidal fluctuations, also confirms the high permeability of the aquifer system; the average value is higher than published values for other volcanic islands. A mixing process for fresh water and seawater, often coupled with ion-exchange mechanisms, explains the groundwater composition, which is mainly of the sodium-chloride type. The water salinity influences the groundwater quality, resulting in a chloride content that exceeds the recommended chloride limit in 91% of the wells . Water–rock interactions are dominant in the chemical evolution of the perched water, which is characterized by bicarbonate-anion type water. Electronic Publication     

Alpine tectonic evolution and thermal water circulations of the Argentera Massif (South-Western Alps)

Three groups of thermal springs with temperatures close to 70 °C discharge both in the core (at Bagni di Vinadio and Terme di Valdieri) and on the external margin (at Berthemont-Les-Bains) of the Argentera Massif. Detailed structural field analysis carried out on the hydrothermal sites allows us to delineate both a model of Alpine tectonic evolution of the Argentera Massif and the patterns of hydrothermal circulation that were active during its final exhumation. The observed fault rock assemblages provide information relative to deformation that occurred in viscous, frictional-to-viscous and frictional crustal regimes. During the Early Miocene, the Bersezio Fault Zone and the Fremamorta Shear Zone, two main mylonitic shear zones, mainly accommodated regional transpression and provided pathways for fluid flow promoting mineral reactions in greenschist facies. During the Late Miocene–Early Pliocene, frictional-to-viscous deformation affected the massif, which underwent predominant transpression in the internal sectors and extension on the external margin. During the Plio-Pleistocene, deformation in frictional condition accompanied the final exhumation of the massif in a transpressive regime and resulted in the development of the NW–SE striking cataclastic zones. The hydraulic properties of these structures mainly influence the patterns of the active thermal circulations and the localization of the recharge and discharge zones. At Berthemont these faults represent conduits, whereas at Vinadio and Valdieri they form complex systems of conduits and barriers. In these two latter sites, the cataclastic faults compose flower structures that constrain laterally the thermal fluid flows while intensely fractured granites sited at depth constitute a highly-transmissive geothermal reservoir. Less permeable migmatitic gneisses overlaying the granites prevent a massive infiltration of the cold fluids at depth. This context favours within the high-permeability fractures granites the development of buoyancy-driven flows which combined with topographically-driven flows, provided the conditions for the upflow of the high-temperature waters.     

深厚覆盖层坝基水文地质结构控渗机理研究

厚度大且为强-弱-强渗透性组合的覆盖层在我国西南河流中分布较广。为研究此类夹层型水文地质结构覆盖层中坝基渗流场的特征及其影响因素,利用Visual MODFLOW模拟了给定覆盖层厚度下强弱渗透性层组的不同渗透性比值、不同厚度比值,和变覆盖层厚度下不同弱透水层厚度这两种类型16种工况下的渗流场,通过监测坝基渗流场中弱透水层顶底板的水头差和水力坡度,研究弱透水层渗透性和厚度对渗流场的影响。结果表明：强-弱-强夹层型水文地质结构的坝基平面渗流场受上层强透水覆盖层性质控制;给定覆盖层厚度工况中,弱透水层顶底板的水头差与强弱透水层渗透性比值相关,两者先表现出正相关的关系,当强弱透水层渗透性比值大于500时,弱透水层顶底板的水头差变化不再明显;变覆盖层厚度工况中,弱透水层顶底板的水头差也呈现出先增大后趋于稳定的趋势,转折点发生在弱透水层厚度20 m的工况条件下;变覆盖层厚度工况中水力坡度则有着随着弱透水层厚度增大而变小的规律。综上所述,在夹层型水文地质结构中,弱透水层厚度与渗透性对渗流场起控制作用。研究结果可为深厚覆盖层河谷区的水电水利工程坝基水害防治提供依据。     

Interactions of diffuse and focused allogenic recharge in an eogenetic karst aquifer (Florida, USA)

The karstic upper Floridan aquifer in north-central Florida (USA) is recharged by both diffuse and allogenic recharge. To understand how recharged water moves within the aquifer, water levels and specific conductivities were monitored and slug tests were conducted in wells installed in the aquifer surrounding the Santa Fe River Sink and Rise. Results indicate that diffuse recharge does not mix rapidly within the aquifer but instead flows horizontally. Stratification may be aided by the high matrix porosity of the eogenetic karst aquifer. Purging wells for sample collection perturbed conductivity for several days, reflecting mixing of the stratified water and rendering collection of representative samples difficult. Interpretive numerical simulations suggest that diffuse recharge impacts the intrusion of allogenic water from the conduit by increasing hydraulic head in the surrounding aquifer and thereby reducing influx to the aquifer from the conduit. In turn, the increase of head within the conduits affects flow paths of diffuse recharge by moving newly recharged water vertically as the water table rises and falls. This movement may result in a broad vertical zone of dissolution at the water table above the conduit system, with thinner and more focused water-table dissolution at greater distance from the conduit.     

Prediction of seasonal water-table fluctuations in La Pampa and Buenos Aires, Argentina

The fluctuation of the water table east of La Pampa province and northwest of Buenos Aires province, Argentina, influences agricultural production in the region because it is closely related to the alternation of dry and wet periods. Sea-surface temperature (SST) anomalies have been used as predictors to forecast atmospheric variables in different regions of the world. The objective of this work is to present a simple model to forecast seasonal rainfall using SST distribution in the Pacific Ocean as a predictor. Once the relationship between precipitation and water-table fluctuations was established, a methodology for the prediction of water-table fluctuations was developed. A good agreement between observed and predicted water-table fluctuations was found when estimating water-table fluctuations in the summer and autumn seasons. Electronic Publication     

Modeling water-table fluctuations in a sloping aquifer with random hydraulic conductivity

To prevent environmental problems like water logging and increase in soil salinity which are responsible for the degradation of the top productive soils, an optimum ditch drainage design is required. For this purpose a knowledge of the spatio-temporal variation of the water table is essential. In this study the spatio-temporal variation of the water table in a sloping ditch drainage system has been modeled from a stochastic point of view, incorporating randomness in hydraulic conductivity to get the expression for the mean and the standard deviation of the water-table height. The hydraulic conductivity has been considered to be a realization of a log-normal distribution. Application of these expressions in the prediction of mean water-table variation with the associated error bounds has been demonstrated with the help of a ditch drainage problem of a sloping aquifer. The sensitivity analysis has also been carried out to see the effect of variability in the hydraulic conductivity on the water-table fluctuations. The error bounds quantified on the water-table height will thus help in the decision-making process for proper drainage design.     

Hydrogeochemical exploration and sources of anomalous waters

Fifty-four surface and twenty drill-hole waters were sampled in the vicinity of the old Eustis mine in the Eastern Townships of Quebec. Some of these waters were tested periodically over a three-month period. The purpose of the investigation was to develop a method for locating sources of anomalous waters using chemical analyses together with geomorphic and geologic evidence. Techniques have been developed which use chemical and hydrologic theory to estimate travel distance and location of the approximate source area of a water. The travel distance can be approximated as a function of hydraulic head, sodium concentration in the water, and an approximate constant which must be derived for the area.Water samples having anomalous concentrations of copper or zinc can be ascribed, in the specific region of this investigation, partly to contamination from mine workings (35%), partly to known mineralization (25%), and partly to unknown causes (40%).     

Multiple-methods investigation of recharge at a humid-region fractured rock site,Pennsylvania, USA

Lysimeter-percolate and well-hydrograph analyses were combined to evaluate recharge for the Masser Recharge Site (central Pennsylvania, USA). In humid regions, aquifer recharge through an unconfined low-porosity fractured-rock aquifer can cause large magnitude water-table fluctuations over short time scales. The unsaturated hydraulic characteristics of the subsurface porous media control the magnitude and timing of these fluctuations. Data from multiple sets of lysimeters at the site show a highly seasonal pattern of percolate and exhibit variability due to both installation factors and hydraulic property heterogeneity. Individual event analysis of well hydrograph data reveals the primary influences on water-table response, namely rainfall depth, rainfall intensity, and initial water-table depth. Spatial and seasonal variability in well response is also evident. A new approach for calculating recharge from continuous water-table elevation records using a master recession curve (MRC) is demonstrated. The recharge estimated by the MRC approach when assuming a constant specific yield is seasonal to a lesser degree than the recharge estimate resulting from the lysimeter analysis. Partial reconciliation of the two recharge estimates is achieved by considering a conceptual model of flow processes in the highly-heterogeneous underlying fractured porous medium.     

Heat tracing to determine spatial patterns of hyporheic exchange across a river transect

Significant spatial variability of water fluxes may exist at the water-sediment interface in river channels and has great influence on a variety of water issues. Understanding the complicated flow systems controlling the flux exchanges along an entire river is often limited due to averaging of parameters or the small number of discrete point measurements usually used. This study investigated the spatial pattern of the hyporheic flux exchange across a river transect in China, using the heat tracing approach. This was done with measurements of temperature at high spatial resolution during a 64-h monitoring period and using the data to identify the spatial pattern of the hyporheic exchange flux with the aid of a one-dimensional conduction-advection-dispersion model (VFLUX). The threshold of neutral exchange was considered as 126 L m^?2 d^?1 in this study and the heat tracing results showed that the change patterns of vertical hyporheic flux varied with buried depth along the river transect; however, the hyporheic flux was not simply controlled by the streambed hydraulic conductivity and water depth in the river transect. Also, lateral flow dominated the hyporheic process within the shallow high-permeability streambed, while the vertical flow was dominant in the deep low-permeability streambed. The spatial pattern of hyporheic exchange across the river transect was naturally controlled by the heterogeneity of the streambed and the bedform of the stream cross-section. Consequently, a two-dimensional conceptual illustration of the hyporheic process across the river transect is proposed, which could be applicable to river transects of similar conditions.     

The effects of Holocene barrier-island evolution on water-table elevations, Hatteras Island, North Carolina, USA

 On Hatteras Island, North Carolina, USA, complex stratigraphy results from surficial sediments being superimposed on an inherited Pleistocene platform. Numerical simulations used to simulate water-table profiles on the island illustrate the influence that the complex stratigraphy has on water-table elevations. Field data consist of water-table profiles collected from a cross-island transect of wells. These profiles show unusually high water-table elevations in the south-central portion of Hatteras Island. Geophysical data are used to constrain the subsurface stratigraphy. Simple analytical solutions to the groundwater flow equation using laterally varying recharge rates cannot match mean water-table profiles. More complex numerical simulations, which do match the profiles, indicate that the elevated water levels result from a low-permeability unit that extends several kilometers along the longitudinal axis of the island. Island geomorphology and ground-penetrating radar (GPR) data suggest that this unit is a former interdunal wetland that has been buried by a series of parabolic dunes. Received, September 1999 / Revised, March 2000 / Accepted, March 2000     

含底坡项浅水方程的一种数值求解方法

激波捕捉格式能处理含强间断的浅水方程,因为必须采用关于守恒变量的守恒形式方程,作为水流实际驱动力的水位梯度项被分解为压力梯度项和底坡项,如果离散格式不能很好地平衡这两项的作用,将导致不准确的结果。基于理论推导,提出了一种新方法,在单元内部采用非守恒形式,保留水位梯度项,在界面处采用守恒格式捕捉间断。利用4个具有代表性的算例对提出的方法进行了验证,结果表明,该方法能够处理变化地形上的强间断和流态变化,并能够准确捕捉扰动的传播过程。     

Helium in soil and overburden gas as an exploration pathfinder — an assessment

Orientation surveys were conducted over five deposits to test the potential of determining helium in overburden gas as a pathfinder for uranium mineralization and other deposits containing uranium or thorium. Samples were collected via fixed tubes emplaced at depths of 6 m in backfilled holes drilled for this purpose. Compared to the atmospheric background value of 5.24 ppm v/v He, a variable weak anomaly (maximum 5.45–5.65 ppm He) was found over part of the Angela uranium deposit, N.T., in an arid area where mineralization is mostly at a depth of 60–90 m, at or below the water-table. Helium contents were mostly at background levels over a uranium deposit in the Officer Basin, W.A., where mineralization is at the water-table at 30–35 m, although radon gave a marked anomaly. Neither helium nor radon indicated the Manyingee deposit, W.A., which has uranium mineralization in a confined aquifer at 60–110 m. Similarly, no helium anomalies were found over the uranium- and thorium-rich Mt. Weld carbonatite or mineral sands at Eneabba.There appeared to be no correlation between helium distributions shown by groundwater and overburden-gas sampling at Manyingee or Mt. Weld. At Mt. Weld, groundwaters contained 0.06 to 13.60 μ/l He and overburden gases 5.24–5.47 ppm He, with the higher gas concentrations over country rock, where waters had background helium contents. It is presumed that equilibration between overburden gas and the atmosphere is far more rapid than that between overburden gas and groundwater, so that any helium released from the water is quickly dispersed.Overburden-gas helium concentrations were found to vary according to overburden type, being 5.24–5.32 ppm where sandy and porous and 5.30–5.50 ppm where clay-rich and less permeable. These background variations, which are greater than the total background-anomaly contrasts reported in the literature for shallow soil gases, have not been accounted for in most trial surveys, nor has the possibility of similar variations being due to analytical error. From the data obtained, there is little evidence that helium can be considered an effective pathfinder for blind or concealed deposits using soil gas or overburden gas as sample media. Previous work on the use of soil and soil-gas helium determinations in uranium exploration is reviewed in the light of these findings and the concept and techniques assessed.     

Hydrogeological properties of fault zones in a karstified carbonate aquifer (Northern Calcareous Alps,Austria)

This study presents a comparative, field-based hydrogeological characterization of exhumed, inactive fault zones in low-porosity Triassic dolostones and limestones of the Hochschwab massif, a carbonate unit of high economic importance supplying 60 % of the drinking water of Austria’s capital, Vienna. Cataclastic rocks and sheared, strongly cemented breccias form low-permeability (<1 mD) domains along faults. Fractured rocks with fracture densities varying by a factor of 10 and fracture porosities varying by a factor of 3, and dilation breccias with average porosities >3 % and permeabilities >1,000 mD form high-permeability domains. With respect to fault-zone architecture and rock content, which is demonstrated to be different for dolostone and limestone, four types of faults are presented. Faults with single-stranded minor fault cores, faults with single-stranded permeable fault cores, and faults with multiple-stranded fault cores are seen as conduits. Faults with single-stranded impermeable fault cores are seen as conduit-barrier systems. Karstic carbonate dissolution occurs along fault cores in limestones and, to a lesser degree, dolostones and creates superposed high-permeability conduits. On a regional scale, faults of a particular deformation event have to be viewed as forming a network of flow conduits directing recharge more or less rapidly towards the water table and the springs. Sections of impermeable fault cores only very locally have the potential to create barriers.     

Cycling of trace metals (Mn, Fe, Mo, U, V, Cr) in deep pore waters of intertidal flat sediments

Trace metals (Mn, Fe, Mo, U, Cr, V) were studied in pore waters of an intertidal flat located in the German Wadden Sea. The study system is an example of a permeable tidal flat system where pore water exchange is affected by tidal driven pressure gradients besides diffusion. Permanently installed in situ samplers were used to extract pore waters down to 5 m depth throughout one year. The samplers were either located close to the tidal flat margin or in central parts of the tidal flat. Despite dynamic sedimentological and hydrological conditions, the general trends with depth in deep tidal flat pore waters are remarkably similar to those observed in deep sea environments. Rates of trace metal cycling must be comparably large in order to maintain the observed pore water profiles. Trace metals further show similar general trends with depth close to the margin and in central parts of the tidal flat. Seasonal sampling revealed that V and Cr vary concurrent with seasonal changes in dissolved organic carbon (DOC) concentration. This effect is most notable close to the tidal flat margin where sulphate, DOC, and nutrients vary with season down to some metres depth. Seasonal variations of Mn, Fe, Mo, and U are by contrast limited to the upper decimetres of the sediment. Their seasonal patterns depend on organic matter supply, redox stratification, and particulate matter deposited on sediment surfaces. Pore water sampling within one tidal cycle provides evidence for pore water advection in margin sediments. During low tide pore water flow towards the creekbank is generated by a hydraulic gradient suggesting that deep pore waters may be seeping out of creekbank sediments. Owing to the enrichment of specific elements like Mn in pore water compared to sea water, seeping pore waters may have an impact on the chemistry of the open water column. Mass balance calculations reveal that the impact of deep pore waters on the Mn budget in the open water column is below 4%. Mn deep pore water discharge of the whole Wadden Sea is estimated to be about 9% of the total dissolved riverine Mn input into the Southern North Sea.     

Relationship between precipitation and water-table fluctuation in a coastal dune aquifer: northeastern coast of the Buenos Aires province, Argentina

The water-table fluctuation (WTF) method is one of the most widely used means to estimate aquifer recharge. In the northeastern coast of the Buenos Aires province, Argentina, the geomorphological and climatic characteristics, as well as the presence of a shallow, homogeneous unconfined aquifer, make it possible to apply this methodology. The relationship between water-table fluctuations and precipitation in a humid climate, considering its seasonal variations, is assessed. Water tables were measured monthly between February 2008 and September 2010 in a monitoring network; rainfall data were analysed. The water table rises when the accumulated precipitation between measurements is more than 53?mm/month in the dry season and more than 97?mm/month in the rainy season. The index, relating water-table fluctuations and precipitation occurring between measurements, shows that values below 0 suggest no increase in the water reserves, whereas higher values entail an increase. In the study area, where there is a lack of historical data, finding out the relationship between water-table fluctuations and precipitation will constitute a tool for groundwater use and management, and set up an early warning system for dry periods. It could also be extrapolated to other regions with similar hydrological conditions lacking in data.