Importance of unsaturated zone flow for simulating recharge in a humid climate

Transient recharge to the water table is often not well understood or quantified. Two approaches for simulating transient recharge in a ground water flow model were investigated using the Trout Lake watershed in north-central Wisconsin: (1) a traditional approach of adding recharge directly to the water table and (2) routing the same volume of water through an unsaturated zone column to the water table. Areas with thin (less than 1 m) unsaturated zones showed little difference in timing of recharge between the two approaches; when water was routed through the unsaturated zone, however, less recharge was delivered to the water table and more discharge occurred to the surface because recharge direction and magnitude changed when the water table rose to the land surface. Areas with a thick (15 to 26 m) unsaturated zone were characterized by multimonth lags between infiltration and recharge, and, in some cases, wetting fronts from precipitation events during the fall overtook and mixed with infiltration from the previous spring snowmelt. Thus, in thicker unsaturated zones, the volume of water infiltrated was properly simulated using the traditional approach, but the timing was different from simulations that included unsaturated zone flow. Routing of rejected recharge and ground water discharge at land surface to surface water features also provided a better simulation of the observed flow regime in a stream at the basin outlet. These results demonstrate that consideration of flow through the unsaturated zone may be important when simulating transient ground water flow in humid climates with shallow water tables.     

Use of temperature profiles and stable isotopes to trace flow lines: Nagaoka area, Japan

In this study, we use borehole temperature data and stable isotopes to delineate the flow system and estimate the effect of urbanization in the Nagaoka area of Japan. Temperature profiles were measured four times in observation wells during the period 2000-2001 and compared with those measured in the same wells during the period 1977-1983 (Taniguchi 1986). Water was sampled in both observation and pumping wells during the same period. The temporal and spatial variability in temperature indicate clearly the effect of urban warming and heavy pumping on the ground water system. Urban warming caused higher temperatures recently as compared to the older values, and pumping caused induced recharge from the river to the ground water. The stable isotope data show the ground water flow system is divided into shallow, intermediate, and deep systems, and that land use and infiltration rate are affecting the shallow flow system.     

Using MODFLOW 2000 to model ET and recharge for shallow ground water problems

In environments with shallow ground water elevation, small changes in the water table can cause significant variations in recharge and evapotranspiration fluxes. Particularly, where ground water is close to the soil surface, both recharge and evapotranspiration are regulated by a thin unsaturated zone and, for accuracy, must be represented using nonconstant and often nonlinear relationships. The most commonly used ground water flow model today, MODFLOW, was originally designed with a modular structure with independent packages representing recharge and evaporation processes. Systems with shallow ground water, however, may be better represented using either a recharge function that varies with ground water depth or a continuous recharge and evapotranspiration function that is dependent on depth to water table. In situations where the boundaries between recharging and nonrecharging cells change with time, such as near a seepage zone, a continuous ground water flux relationship allows recharge rates to change with depth rather than having to calculate them at each stress period. This research article describes the modification of the MODFLOW 2000 recharge and segmented evapotranspiration packages into a continuous recharge-discharge function that allows ground water flux to be represented as a continuous process, dependent on head. The modifications were then used to model long-term recharge and evapotranspiration processes on a saline, semiarid floodplain in order to understand spatial patterns of salinization, and an overview of this process is given.     

Heat flow and ground water movement in the Bohemian cretaceous basin (Czechoslovakia)

Sub-surface temperature fields may be considerably affected by active ground water systems, thereby seriously hampering the interpretation of heat flow data. Quantitative evaluation of the convective component of heat transfer is thus very important in cases such as large sedimentary basins with vast underground water circulation. We propose in this study a simple model of horizontal aquifer. This model was used to examine the effect of the lateral convection on the surface heat flow near the recharge zone of basinal margins. The perturbation of the heat flow field above the aquifer was calculated for various aquifer geometry and various flow velocities and the regional scale dependence of the perturbation on the hydraulic properties of the aquifer was demonstrated. The model was applied to the Bohemian Cretaceous Basin and it was shown that within a few kilometres from the recharge zones the observed surface heat flow may be underestimated by up to several tens of percent. The procedure was tested in two locations in this area, in an attempt to make hydrogeological corrections to the measured heat flow values in several boreholes.     

断裂带地下水位的降雨动态模型及异常识别

穿越断裂带的观测井能够提供特殊的构造运动信息,其水位动态是构造活动、降雨等各种因素综合作用的结果.对由降雨引起的水位动态进行模型解译,有利于更加准确地提取断裂带的构造活动异常.根据断裂带地下水循环的基本特征进行简化,利用响应函数法描述入渗补给的滞后特征,建立了降雨-水位动态的组合水箱模型.该模型被用于分析北京八宝山断裂带大灰厂观测孔2008年上半年的水位动态,发现2008年5月存在显著的异常.这种异常可能与局部或区域的构造活动有关,也可能是断裂带结构参数变化的结果.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

Impact of Recharge Through Residual Oil Upon Sampling of Underlying Ground Water

At an aviation gasoline spill site in Traverse City, Michigan, historical records indicate a positive correlation between significant rainfall events and increased concentrations of slightly soluble organic compounds in the monitoring wells of the site. To investigate the recharge effect on ground water quality due to infiltrating, water percolating past residual oil and into the saturated zone, an in situ infiltration experiment was performed at the site. Sampling cones were set at various depths below a circular test area, 13 feet (4 meters) in diameter. Rainfall was simulated by sprinkling the test area at a rate sufficiently low to prevent runoff. The sampling cones for soil-gas and ground water quality were installed in the unsaturated and saturated zones to observe the effects of the recharge process. At the time of the test, the water table was below the residual oil layer. The responses of the soil-gas and ground water quality were monitored during the recharge and drainage periods, which resulted from the sprinkling.
Infiltrated water was determined to have transported organic constituents of the residual oil, specifically benzene, toluene, ethylbenzene, and ortho-xylene (BTEX), into the ground water beneath the water table, elevating the aqueous concentrations of these constituents in the saturated zone. Soil-gas concentrations of the organic compounds in the unsaturated zone increased with depth and time after the commencement of infiltration. Reaeration of the unconfined aquifer via the infiltrated water was observed. It is concluded that water quality measurements are directly coupled to recharge events for the sandy type of aquifer with an overlying oil phase, which was studied in this work. Ground water sampling strategies and data analysis need to reflect the effect of recharge from precipitation on shallow, unconfined aquifers where an oil phase may be present.     

Modeling decadal timescale interactions between surface water and ground water in the central Everglades, Florida, USA

Surface-water and ground-water flow are coupled in the central Everglades, although the remoteness of this system has hindered many previous attempts to quantify interactions between surface water and ground water. We modeled flow through a 43,000 ha basin in the central Everglades called Water Conservation Area 2A. The purpose of the model was to quantify recharge and discharge in the basin's vast interior areas. The presence and distribution of tritium in ground water was the principal constraint on the modeling, based on measurements in 25 research wells ranging in depth from 2 to 37 m. In addition to average characteristics of surface-water flow, the model parameters included depth of the layer of ‘interactive’ ground water that is actively exchanged with surface water, average residence time of interactive ground water, and the associated recharge and discharge fluxes across the wetland ground surface. Results indicated that only a relatively thin (8 m) layer of the 60 m deep surfical aquifer actively exchanges surface water and ground water on a decadal timescale. The calculated storage depth of interactive ground water was 3.1 m after adjustment for the porosity of peat and sandy limestone. Modeling of the tritium data yielded an average residence time of 90 years in interactive ground water, with associated recharge and discharge fluxes equal to 0.01 cm d⁻¹. ³H/³He isotopic ratio measurements (which correct for effects of vertical mixing in the aquifer with deeper, tritium-dead water) were available from several wells, and these indicated an average residence time of 25 years, suggesting that residence time was overestimated using tritium measurements alone. Indeed, both residence time and storage depth would be expected to be overestimated due to vertical mixing. The estimate of recharge and discharge (0.01 cm d⁻¹) that resulted from tritium modeling therefore is still considered reliable, because the ratio of residence time and storage depth (used to calculated recharge and discharge) is much less sensitive to vertical mixing compared with residence time alone. We conclude that a small but potentially significant component of flow through the Everglades is recharged to the aquifer and stored there for years to decades before discharged back to surface water. Long-term storage of water and solutes in the ground-water system beneath the wetlands has implications for restoration of Everglades water quality.     

Vadose zone-attenuated artificial recharge for input to a ground water model

Accurate representation of artificial recharge is requisite to calibration of a ground water model of an unconfined aquifer for a semiarid or arid site with a vadose zone that imparts significant attenuation of liquid transmission and substantial anthropogenic liquid discharges. Under such circumstances, artificial recharge occurs in response to liquid disposal to the vadose zone in areas that are small relative to the ground water model domain. Natural recharge, in contrast, is spatially variable and occurs over the entire upper boundary of a typical unconfined ground water model. An improved technique for partitioning artificial recharge from simulated total recharge for inclusion in a ground water model is presented. The improved technique is applied using data from the semiarid Hanford Site. From 1944 until the late 1980s, when Hanford's mission was the production of nuclear materials, the quantities of liquid discharged from production facilities to the ground vastly exceeded natural recharge. Nearly all hydraulic head data available for use in calibrating a ground water model at this site were collected during this period or later, when the aquifer was under the diminishing influence of the massive water disposals. The vadose zone is typically 80 to 90 m thick at the Central Plateau where most production facilities were located at this semiarid site, and its attenuation of liquid transmission to the aquifer can be significant. The new technique is shown to improve the representation of artificial recharge and thereby contribute to improvement in the calibration of a site-wide ground water model.     

Ground Water Sampling Bias Observed in Shallow, Conventional Wells

A previous field demonstration project on nitrate-based bioremediation of a fuel-contaminated aquifer used short-screened clustered well points in addition to shallow (10 foot), conventional monitoring wells to monitor the progress of remediation during surface application of recharge. These well systems were placed in the center and at one edge of each of two treatment cells. One cell received recharge amended with nitrate (nitrate cell), and the other received unamended recharge (control cell). Data from the clustered well points were averaged to provide a mean estimate for comparison with the associated conventional monitoring well.
Conservative tracer profiles were similar for each of the four systems, with better fits obtained for well systems located at the edge of the treatment cells. However, aromatic hydrocarbon and electron acceptor profiles varied greatly for the two center well systems, with the conventional monitoring well data suggesting that remediation was proceeding at a much more rapid rate than indicated by the cluster well points. Later tests with an electromagnetic borehole flowmeter demonstrated a significant vertical flow through the well-bore of the conventional monitoring well under simulated operating conditions. This created an artifact during sampling, thought to arise from preferential flow of recharge water from the water table to deeper portions of the contaminated zone resulting in several effects, including an actual decreased residence time of water sampled by the conventional well. These data provide additional evidence that conventional monitoring wells may be inadequate for monitoring remediation in the presence of significant vertical hydraulic gradients, even for fairly shallow homogeneous aquifers.     

On the interpretation of recharge estimates from steady-state model calibrations

Ground water recharge is often estimated through the calibration of ground water flow models. We examine the nature of calibration errors by considering some simple mathematical and numerical calculations. From these calculations, we conclude that calibrating a steady-state ground water flow model to water level extremes yields estimates of recharge that have the same value as the time-varying recharge at the time the water levels are measured. These recharge values, however, are a subdued version of the actual transient recharge signal. In addition, calibrating a steady-state ground water flow model to data collected during periods of rising water levels will produce recharge values that underestimate the actual transient recharge. Similarly, calibrating during periods of falling water levels will overestimate the actual transient recharge. We also demonstrate that average water levels can be used to estimate the actual average recharge rate provided that water level data have been collected for a sufficient amount of time.     

Simulation of submarine ground water discharge to a marine estuary: Biscayne Bay, Florida

Variable density ground water flow models are rarely used to estimate submarine ground water discharge because of limitations in computer speed, data availability, and availability of a simulation tool that can minimize numerical dispersion. This paper presents an application of the SEAWAT code, which is a combined version of MODFLOW and MT3D, to estimate rates of submarine ground water discharge to a coastal marine estuary. Discharge rates were estimated for Biscayne Bay, Florida, for the period from January 1989 to September 1998 using a three-dimensional, variable density ground water flow and transport model. Hydrologic stresses in the 10-layer model include recharge, evapotranspiration, ground water withdrawals from municipal wellfields, interactions with surface water (canals in urban areas and wetlands in the Everglades), boundary fluxes, and submarine ground water discharge to Biscayne Bay. The model was calibrated by matching ground water levels in monitoring wells, baseflow to canals, and the position of the 1995 salt water intrusion line. Results suggest that fresh submarine ground water discharge to Biscayne Bay may have exceeded surface water discharge during the 1989, 1990, and 1991 dry seasons, but the average discharge for the entire simulation period was only approximately 10% of the surface water discharge to the bay. Results from the model also suggest that tidal canals intercept fresh ground water that might otherwise have discharged directly to Biscayne Bay. This application demonstrates that regional scale variable density models are potentially useful tools for estimating rates of submarine ground water discharge.     

Hydrogeochemical processes involved in salt-dissolution zones,texas panhandle,U.S.A.

Permian evaporite deposits have been extensively dissolved beneath the perimeter of the Southern High Plains in the Texas Panhandle. Hydrologic and geochemical data were collected from six test wells to determine hydrogeochemical processes involved and the source and flow paths of ground water moving in salt-dissolution zones. Geochemical similarities and hydraulic-head relationships indicate that ground water dissolving halite and anhydrite moves downward from aquifers in post-Permian formations and follows flow paths influenced by topography. Holocene salt-dissolution rates probably are lower than Tertiary and Pleistocene rates owing to regional changes in physiography and climate that probably decreased the amount of recharge to salt-dissolution zones. Present as well as palaeohydrologic ground-water velocities and salt-dissolution rates are probably less beneath the Southern High Plains than in adjacent, peripheral salt-dissolution zones because of lower hydraulic conductivities and lower hydraulic-head gradients. Salinities in peripheral salt-dissolution zones are low (67 000 to 95 000 mg L^?1) despite high solubility of halite, reflecting relatively open circulation of ground water. In interior salt-dissolution zones beneath the Southern High Plains, ground-water circulation is low and water composition tends to reach halite saturation.     

Distinguishing sources of ground water recharge by using delta2H and delta18O

Stable isotope values of hydrogen and oxygen from precipitation and ground water samples were compared by using a volumetrically based mixing equation and stable isotope gradient to estimate the season and location of recharge in four basins. Stable isotopes were sampled at 11 precipitation sites of differing elevation during a 2-year period to quantify seasonal stable isotope contributions as a function of elevation. Supplemental stable isotope data collected by the International Atomic Energy Association during a 14-year period were used to reduce annual variability of the mean seasonal stable isotope data. The stable isotope elevation relationships and local precipitation elevation relationships were combined by using a digital elevation model to calculate the total volumetric contribution of water and stable isotope values as a function of elevation within the basins. The results of these precipitation calculations were compared to measured ground water stable isotope values at the major discharge points near the terminus of the basins. Volumetric precipitation contributions to recharge were adjusted to isolate contributing elevations. This procedure provides an improved representation of recharge contributions within the basins over conventional stable isotope methods. Stable isotope values from wells and springs at the terminus of each basin were used to infer the elevations of precipitation important for recharge of the regional ground water flow system. Ancillary climatic, geologic, and stable isotope values were used to further constrain the location where precipitation is entering the ground water flow system.     

Salinization of a fresh palaeo-ground water resource by enhanced recharge

Deterioration of fresh ground water resources caused by salinization is a growing issue in many arid and semi-arid parts of the world. We discuss here the incipient salinization of a 10(4) km2 area of fresh ground water (<3,000 mg/L) in the semiarid Murray Basin of Australia caused by widespread changes in land use. Ground water 14C concentrations and unsaturated zone Cl soil water inventories indicate that the low salinity ground water originated mainly from palaeo-recharge during wet climatic periods more than 20,000 years ago. However, much of the soil water in the 20 to 60 m thick unsaturated zone throughout the area is generally saline (>15,000 mg/L) because of relatively high evapotranspiration during the predominantly semiarid climate of the last 20,000 years. Widespread clearing of native vegetation over the last 100 years and replacement with crops and pastures leads to enhancement of recharge rates that progressively displace the saline soil-water from the unsaturated zone into the ground water. To quantify the impact of this new hydrologic regime, a one-dimensional model that simulates projected ground water salinities as a function of depth to ground water, recharge rates, and soil water salt inventory was developed. Results from the model suggest that, in some areas, the ground water salinity within the top 10 m of the water table is likely to increase by a factor of 2 to 6 during the next 100 years. Ground water quality will therefore potentially degrade beyond the point of usefulness well before extraction of the ground water exhausts the resource.     

Dynamics of flood water infiltration and ground water recharge in hyperarid desert

A study on flood water infiltration and ground water recharge of a shallow alluvial aquifer was conducted in the hyperarid section of the Kuiseb River, Namibia. The study site was selected to represent a typical desert ephemeral river. An instrumental setup allowed, for the first time, continuous monitoring of infiltration during a flood event through the channel bed and the entire vadose zone. The monitoring system included flexible time domain reflectometry probes that were designed to measure the temporal variation in vadose zone water content and instruments to concurrently measure the levels of flood and ground water. A sequence of five individual floods was monitored during the rainy season in early summer 2006. These newly generated data served to elucidate the dynamics of flood water infiltration. Each flood initiated an infiltration event which was expressed in wetting of the vadose zone followed by a measurable rise in the water table. The data enabled a direct calculation of the infiltration fluxes by various independent methods. The floods varied in their stages, peaks, and initial water contents. However, all floods produced very similar flux rates, suggesting that the recharge rates are less affected by the flood stages but rather controlled by flow duration and available aquifer storage under it. Large floods flood the stream channel terraces and promote the larger transmission losses. These, however, make only a negligible contribution to the recharge of the ground water. It is the flood duration within the active streambed, which may increase with flood magnitude that is important to the recharge process.     

Controls on Ground Water Chemistry in the Central Couloir Sud Rifain,Morocco

Irrigation, urbanization, and drought pose challenges for the sustainable use of ground water in the central Couloir sud rifain, a major agricultural region in north-central Morocco, which includes the cities of Fès and Meknès. The central Couloir is underlain by unconfined and confined carbonate aquifers that have suffered declines in hydraulic head and reductions in spring flow in recent decades. Previous studies have surveyed ground water flow and water quality in wells and springs but have not comprehensively addressed the chemistry of the regional aquifer system. Using graphical techniques and saturation index calculations, we infer that major ion chemistry is controlled (1) in the surficial aquifer by cation exchange, calcite dissolution, mixing with deep ground water, and possibly calcite precipitation and (2) in the confined aquifer and warm springs by calcite dissolution, dolomite dissolution, mixing with water that has dissolved gypsum and halite, and calcite precipitation. Analyses of ²H and ¹⁸O indicate that shallow ground water is affected by evaporation during recharge (either of infiltrating precipitation or return flow), whereas deep ground water is sustained by meteoric recharge with little evaporation. Mechanisms of recharge and hydrochemical evolution are broadly consistent with those delineated for similar regional aquifer systems elsewhere in Morocco and in southern Spain.     

The Application of In Situ Air Sparging as an Innovative Soils and Ground Water Remediation Technology

Vapor extraction (soil venting) has been demonstrated to be a successful and cost-effective remediation technology for removing VOCs from the vadose (unsaturated) zone. However, in many cases, seasonal water table fluctuations, drawdown associated with pump-and-treat remediation techniques, and spills involving dense, non-aqueous phase liquids (DNAPLS) create contaminated soil below the water table. Vapor extraction alone is not considered to be an optimal remediation technology to address this type of contamination.
An innovative approach to saturated zone remediation is the use of sparging (injection) wells to inject a hydrocarbon-free gaseous medium (typically air) into the saturated zone below the areas of contamination. The contaminants dissolved in the ground water and sorbed onto soil particles partition into the advective air phase, effectively simulating an in situ air-stripping system. The stripped contaminants are transported in the gas phase to the vadose zone, within the radius of influence of a vapor extraction and vapor treatment system.
In situ air sparging is a complex multifluid phase process, which has been applied successfully in Europe since the mid-1980s. To date, site-specific pilot tests have been used to design air-sparging systems. Research is currently underway to develop better engineering design methodologies for the process. Major design parameters to be considered include contaminant type, gas injection pressures and flow rates, site geology, bubble size, injection interval (areal and vertical) and the equipment specifications. Correct design and operation of this technology has been demonstrated to achieve ground water cleanup of VOC contamination to low part-per-billion levels.     

Demonstration of a Microbiologically Enhanced Vertical Ground Water Circulation Well Technology at a Superfund Site

A full-scale ground water circulation well (GCW) system was installed and operated to demonstrate in situ remediation of soil and ground water impacted with a mixture of chlorinated and nonchlorinated organic compounds at a Superfund site in upstate New York. System performance and applicability under site-specific conditions were evaluated based on the system's ability to meet the New York State Department of Environmental Conservation (NYSDEC) cleanup goals for target compounds in ground water and soil. Contaminants from the unsaturated zone were mobilized (volatilized) by one-way vacuum extraction, and treated via enhanced biodegradation (bioventing). In the saturated zone, contaminants were mobilized by soil flushing (solubilized) and treated by a combination of air stripping and biodegradation. An in situ aqueous phase bioreactor, and an ex situ gas phase bioreactor, were integrated into the system to enhance treatment via bioremediation. After 15 months of operation, the mass of target contaminants in soil and ground water combined had been reduced by 75%. Removal by biological mechanisms ranged from 35% to 56% of the total observed mass reduction. The in situ and the ex situ bioreactors mineralized 79% and 76%, respectively, of their target biodegradable contaminant loads. Results indicate that some mass reduction in target contaminants may have been from aerobic and aerobic processes within the circulation cell. Nonchlorinated compounds were relatively easy to mobilize (volatilize, solubilize, and/or transport) and treat when compared to chlorinated compounds. The data collected during the 15-month study indicate that remediation could be accomplished at the Sweden-3 Chapman site using the technology tested.     

Delineation of regional arid karstic aquifers: an integrative data approach

This research integrates data procedures for the delineation of regional ground water flow systems in arid karstic basins with sparse hydrogeologic data using surface topography data, geologic mapping, permeability data, chloride concentrations of ground water and precipitation, and measured discharge data. This integrative data analysis framework can be applied to evaluate arid karstic aquifer systems globally. The accurate delineation of ground water recharge areas in developing aquifer systems with sparse hydrogeologic data is essential for their effective long-term development and management. We illustrate the use of this approach in the Cuatrociénegas Basin (CCB) of Mexico. Aquifers are characterized using geographic information systems for ground water catchment delineation, an analytical model for interbasin flow evaluation, a chloride balance approach for recharge estimation, and a water budget for mapping contributing catchments over a large region. The test study area includes the CCB of Coahuila, Mexico, a UNESCO World Biosphere Reserve containing more than 500 springs that support ground water-dependent ecosystems with more than 70 endemic organisms and irrigated agriculture. We define recharge areas that contribute local and regional ground water discharge to springs and the regional flow system. Results show that the regional aquifer system follows a topographic gradient that during past pluvial periods may have linked the Río Nazas and the Río Aguanaval of the Sierra Madre Occidental to the Río Grande via the CCB and other large, currently dry, upgradient lakes.