Arsenic Geochemistry and Hydrostratigraphy in Midwestern U.S. Glacial Deposits

Arsenic concentrations exceeding the U.S. EPA's 10 μg/L standard are common in glacial aquifers in the midwestern United States. Previous studies have indicated that arsenic occurs naturally in these aquifers in association with metal-(hydr)oxides and is released to groundwater under reducing conditions generated by microbial oxidation of organic matter. Despite this delineation of the arsenic source and mechanism of arsenic mobilization, identification of arsenic-impacted aquifers is hindered by the heterogeneous and discontinuous nature of glacial sediments. In much of the Midwest, the hydrostratigraphy of glacial deposits is not sufficiently characterized to predict where elevated arsenic concentrations are likely to occur. This case study from southeast Wisconsin presents a detailed characterization of local stratigraphy, hydrostratigraphy, and geochemistry of the Pleistocene glacial deposits and underlying Silurian dolomite. Analyses of a single core, water chemistry data, and well construction reports enabled identification of two aquifers separated by an organic-rich aquitard. The upper, unconfined aquifer provides potable water, whereas arsenic generally exceeds 10 μg/L in the deeper aquifer. Although coring and detailed hydrostratigraphic characterization are often considered impractical, our results demonstrate that a single core improved interpretation of the complex lithology and hydrostratigraphy. This detailed characterization of hydrostratigraphy facilitated development of well construction guidelines and lays the ground work for further studies of the complex interactions among aquifer sediments, hydrogeology, water chemistry, and microbiology that lead to elevated arsenic in groundwater.     

Correlating Bedrock Folds to Higher Rates of Arsenic Detection in Groundwater,Southeast Wisconsin,USA

Arsenic in private drinking water wells is a significant problem across much of eastern Wisconsin, USA. The release mechanism and stratigraphic distribution of sulfide and iron (hydr)oxide sources of arsenic in bedrock aquifers are well understood for northeastern Wisconsin. However, recent geologic mapping has identified numerous small bedrock folds to the south, and the impact of these geologic structures on local groundwater flow and well contamination has been little studied. This paper examines the hydrologic and structural effects of the Beaver Dam anticline, southeast Wisconsin, on arsenic in groundwater in the region. Multivariate logistic regression shows wells near the Beaver Dam anticline are statistically more likely to detect arsenic in groundwater compared to wells farther away. Structural and hydrologic changes related to folding are interpreted to be the cause. Core drilled near the fold axis is heavily fractured, and many fractures are filled with sulfides. Elevated hydraulic conductivity estimates are also recorded near the fold axis, which may reflect a higher concentration of vertical fractures. These structural and hydrologic changes may have led to systematic changes in the distribution and concentration of arsenic-bearing mineral hosts, resulting in the observed detection pattern. For areas with similar underlying geology, this approach may improve prediction of arsenic risk down to the local level.     

Assessing hydrogeologic controls on dynamic groundwater storage using long‐term instrumental records of water table levels

This study analyzes a long‐term regional compilation of water table response to climate variability based on 124 long‐term groundwater wells distributed across New England, USA, screened in a variety of geologic materials. The New England region of the USA is located in a humid‐temperature climate underlain by low‐storage‐fractured metamorphic and crystalline bedrock dissected by north–south trending valleys filled with glacial and post‐glacial valley fill sediments. Uplands are covered by thin glacial till that comprises more than 60% of the total area. Annual and multi‐annual responses of the water table to climate variability are assessed to understand how local hydraulic properties and hydrogeologic setting (located in recharge/discharge region) of the aquifer influence the hydrologic sensitivity of the aquifer system to climate variability. This study documents that upland aquifer systems dominated by thin deposits of surface till comprise ~70% of the active and dynamic storage of the region. Total aquifer storage changes of +5 to ?7 km³ occur over the region during the study interval. The storage response is dominated by thin and low permeability surficial till aquifer that fills and drains on a multi‐annual basis and serves as the main mechanism to deliver water to valley fill aquifers and underlying bedrock aquifers. Whereas the till aquifer system is traditionally neglected as an important storage reservoir, this study highlights the importance of a process‐based understanding of how different landscape hydrogeologic units contribute to the overall hydrologic response of a region.     

Arsenic in glacial drift aquifers and the implication for drinking water--lower Illinois River Basin

The lower Illinois River Basin (LIRB) covers 47,000 km2 of central and western Illinois. In the LIRB, 90% of the ground water supplies are from the deep and shallow glacial drift aquifers. The deep glacial drift aquifer (DGDA) is below 152 m altitude, a sand and gravel deposit that fills the Mahomet Buried Bedrock Valley, and overlain by more than 30.5 m of clayey till. The LIRB is part of the USGS National Water Quality Assessment program, which has an objective to describe the status and trends of surface and ground water quality. In the DGDA, 55% of the wells used for public drinking-water supply and 43% of the wells used for domestic drinking water supply have arsenic concentrations above 10 micrograms/L (a new U.S. EPA drinking water standard). Arsenic concentrations greater than 25 micrograms/L in ground water are mostly in the form of arsenite (AsIII). The proportion of arsenate (AsV) to arsenite does not change along the flowpath of the DGDA. Because of the limited number of arsenic species analyses, no clear relations between species and other trace elements, major ions, or physical parameters could be established. Arsenic and barium concentrations increase from east to west in the DGDA and are positively correlated. Chloride and arsenic are positively correlated and provide evidence that arsenic may be derived locally from underlying bedrock. Solid phase geochemical analysis of the till, sand and gravel, and bedrock show the highest presence of arsenic in the underlying organic-rich carbonate bedrock. The black shale or coal within the organic-rich carbonate bedrock is a potential source of arsenic. Most high arsenic concentrations found in the DGDA are west and downgradient of the bedrock structural features. Geologic structures in the bedrock are potential pathways for recharge to the DGDA from surrounding bedrock.     

Aquifer Vulnerability Assessment Based on Sequence Stratigraphic and 39Ar Transport Modeling

A large‐scale groundwater flow and transport model is developed for a deep‐seated (100 to 300 m below ground surface) sedimentary aquifer system. The model is based on a three‐dimensional (3D) hydrostratigraphic model, building on a sequence stratigraphic approach. The flow model is calibrated against observations of hydraulic head and stream discharge while the credibility of the transport model is evaluated against measurements of ³⁹Ar from deep wells using alternative parameterizations of dispersivity and effective porosity. The directly simulated 3D mean age distributions and vertical fluxes are used to visualize the two‐dimensional (2D)/3D age and flux distribution along transects and at the top plane of individual aquifers. The simulation results are used to assess the vulnerability of the aquifer system that generally has been assumed to be protected by thick overlaying clayey units and therefore proposed as future reservoirs for drinking water supply. The results indicate that on a regional scale these deep‐seated aquifers are not as protected from modern surface water contamination as expected because significant leakage to the deeper aquifers occurs. The complex distribution of local and intermediate groundwater flow systems controlled by the distribution of the river network as well as the topographical variation (Tóth 1963) provides the possibility for modern water to be found in even the deepest aquifers.     

Estimation of groundwater vulnerability to pollution based on DRASTIC in the Niipele sub-basin of the Cuvelai Etosha Basin,Namibia

Surface water is a scarce resource in Namibia with about sixty percent of Namibia's population dependent on groundwater for drinking purposes. With increasing population, the country faces water challenges and thus groundwater resources need to be managed properly. One important aspect of Integrated Water Resources Management is the protection of water resources, including protection of groundwater from contamination and over-exploitation. This study explores vulnerability mapping as a basic tool for protecting groundwater resources from pollution. It estimates groundwater vulnerability to pollution in the upper Niipele sub-basin of the Cuvelai-Etosha in Northern Namibia using the DRASTIC index. The DRASTIC index uses GIS to estimate groundwater vulnerability by overlaying different spatially referenced hydrogeological parameters that affect groundwater contamination. The study assesses the discontinuous perched aquifer (KDP) and the Ohangwena multi-layered aquifer 1 (KOH-1). For perched aquifers, point data was regionalized by a hydrotope approach whereas for KOH-1 aquifer, inverse distance weighting was used. The hydrotope approach categorized different parts of the hydrogeological system with similar properties into five hydrotopes. The result suggests that the discontinuous perched aquifers are more vulnerable than Ohangwena multi-layered aquifer 1. This implies that vulnerability increases with decreasing depth to water table because contaminants have short travel time to reach the aquifer when they are introduced on land surface. The nitrate concentration ranges between 2 and 288 mg/l in perched aquifers while in Ohangwena multi-layered aquifer 1, it ranges between 1 and 133 mg/l. It was observed that perched aquifers have high nitrate concentrations than Ohangwena 1 aquifer, which correlates well with the vulnerability results.     

Surface water-groundwater exchange dynamics in buried-valley aquifer systems

Groundwater is a primary source of drinking water worldwide, but excess nutrients and emerging contaminants could compromise groundwater quality and limit its usage as a drinking water source. As such contaminants become increasingly prevalent in the biosphere, a fundamental understanding of their fate and transport in groundwater systems is necessary to implement successful remediation strategies. The dynamics of surface water-groundwater (hyporheic) exchange within a glacial, buried-valley aquifer system are examined in the context of their implications for the transport of nutrients and contaminants in riparian sediments. High conductivity facies act as preferential flow pathways which enhance nutrient and contaminant delivery, especially during storm events, but transport throughout the aquifer also depends on subsurface sedimentary architecture (e.g. interbedded high and low conductivity facies). Temperature and specific conductance measurements indicate extensive hyporheic mixing close to the river channel, but surface water influence was also observed far from the stream-aquifer interface. Measurements of river stage and hydraulic head indicate that significant flows during storms (i.e., hot moments) alter groundwater flow patterns, even between consecutive storm events, as riverbed conductivity and, more importantly, the hydraulic connectivity between the river and aquifer change. Given the similar mass transport characteristics among buried-valley aquifers, these findings are likely representative of glacial aquifer systems worldwide. Our results suggest that water resources management decisions based on average (base) flow conditions may inaccurately represent the system being evaluated, and could reduce the effectiveness of remediation strategies for nutrients and emerging contaminants.     

Detecting a Defective Casing Seal at the Top of a Bedrock Aquifer

An improperly sealed casing can produce a direct hydraulic connection between two or more originally isolated aquifers with important consequences regarding groundwater quantity and quality. A recent study by Richard et al. (2014) investigated a monitoring well installed in a fractured rock aquifer with a defective casing seal at the soil–bedrock interface. A hydraulic short circuit was detected that produced some leakage between the rock and the overlying deposits. A falling‐head permeability test performed in this well showed that the usual method of data interpretation is not valid in this particular case due to the presence of a piezometric error. This error is the direct result of the preferential flow originating from the hydraulic short circuit and the subsequent re‐equilibration of the piezometric levels of both aquifers in the vicinity of the inlet and the outlet of the defective seal. Numerical simulations of groundwater circulation around the well support the observed impact of the hydraulic short circuit on the results of the falling‐head permeability test. These observations demonstrate that a properly designed falling‐head permeability test may be useful in the detection of defective casing seals.     

Groundwater flow and implications for groundwater contamination north of Prewitt,New Mexico,U.S.A.

In the southern San Juan Basin, New Mexico, strata of Permian and younger age dip gently toward the center of the basin. Most previous investigators believed that recharge to these strata occurred by precipitation on the outcrops and groundwater flowed downdip to the north and northeast. Recent water-level measurements in an undeveloped part of the basin near Prewitt, New Mexico, show that groundwater at shallow depths in alluvium and bedrock flows southward, opposite to the dip direction, and toward a major ephemeral drainage in a strike valley. North of this area, groundwater in deep bedrock aquifers does appear to flow northward. This information suggests that there are two groundwater circulation patterns; a shallow one controlled by topography and a deeper one controlled by geologic structure.Significant amounts of recharge to sandstone aquifers by infiltration through outcrops is unlikely due to the near-vertical exposures on cliffs, the gentle dip of the strata, and small annual precipitation. Numerical model results suggest that recharge to bedrock aquifers may be from downward leakage via aquitards over large areas and leakage from narrow alluvial aquifers in the subcrop area. The recharge mechanism is controlled by the hydraulic conductivity of the strata.As the flow path is controlled by hydraulic conductivity contrasts, geologic structure, and topography, contamination movement from surface impoundments is likely to be difficult to predict without a thorough hydrogeological site investigation.     

Patterns and Rates of Ground-Water Flow on Long Island, New York

Increased ground-water contamination from human activities on Long Island has prompted studies to define the pattern and rate of ground-water movement. A two-dimensional, fine-mesh, finite-element model consisting of 11,969 nodes and 22,880 elements was constructed to represent ground-water flow along a north-south section through central Long Island. The model represents average hydrologic conditions within a corridor approximately 15 miles wide. The model solves discrete approximations of both the potential and stream functions. The resulting flownet depicts flow paths and defines the vertical distribution of flow within the section. Ground-water flow rates decrease with depth. Sixty-two percent of the water flows no deeper than the upper glacial (water-table) aquifer, 38 percent enters the underlying Magothy aquifer, and only 3.1 percent enters the Lloyd aquifer. The limiting streamlines for flow to the Magothy and Lloyd aquifers indicate that aquifer recharge areas are narrow east-west bands through the center of the island. The recharge area of the Magothy aquifer is only 5.4 miles wide; that of the Lloyd aquifer is less than 0.5 miles. The distribution of ground-water traveltime and a flownet are calculated from model results; both are useful in the investigation of contaminant transport or the chemical evolution of ground water within the flow system. A major discontinuity in traveltime occurs across the streamline which separates the flow subsystems of the two confined aquifers. Water that reaches the Lloyd aquifer attains traveltimes as high as 10,000 years, whereas water that has not penetrated deeper than the Magothy aquifer attains traveltimes of only 2,000 years. The finite-element approach used in this study is particularly suited to ground-water systems that have complex hydrostratigraphy and cross-sectional symmetry.     

Field verification of stable perched groundwater in layered bedrock uplands

Data substantiating perched conditions in layered bedrock uplands are rare and have not been widely reported. Field observations in layered sedimentary bedrock in southwestern Wisconsin, USA, provide evidence of a stable, laterally extensive perched aquifer. Data from a densely instrumented field site show a perched aquifer in shallow dolomite, underlain by a shale-and-dolomite aquitard approximately 25 m thick, which is in turn underlain by sandstone containing a 30-m-thick unsaturated zone above a regional aquifer. Heads in water supply wells indicate that perched conditions extend at least several kilometers into hillsides, which is consistent with published modeling studies. Observations of unsaturated conditions in the sandstone over a 4-year period, historical development of the perched aquifer, and perennial flow from upland springs emanating from the shallow dolomite suggest that perched groundwater is a stable hydrogeologic feature under current climate conditions. Water-table hydrographs exhibit apparent differences in the amount and timing of recharge to the perched and regional flow systems; steep hydraulic gradients and tritium and chloride concentrations suggest there is limited hydraulic connection between the two. Recognition and characterization of perched flow systems have practical importance because their groundwater flow and transport pathways may differ significantly from those in underlying flow systems. Construction of multi-aquifer wells and groundwater withdrawal in perched systems can further alter such pathways.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

Boundary depletion rate and drawdown in leaky wedge‐shaped aquifers

This study presents analytical solutions of the three‐dimensional groundwater flow to a well in leaky confined and leaky water table wedge‐shaped aquifers. Leaky wedge‐shaped aquifers with and without storage in the aquitard are considered, and both transient and steady‐state drawdown solutions are derived. Unlike the previous solutions of the wedge‐shaped aquifers, the leakages from aquitard are considered in these solutions and unlike similar previous work for leaky aquifers, leakage from aquitards and from the water table are treated as the lower and upper boundary conditions. A special form of finite Fourier transforms is used to transform the z‐coordinate in deriving the solutions. The leakage induced by a partially penetrating pumping well in a wedge‐shaped aquifer depends on aquitard hydraulic parameters, the wedge‐shaped aquifer parameters, as well as the pumping well parameters. We calculate lateral boundary dimensionless flux at a representative line and investigate its sensitivity to the aquitard hydraulic parameters. We also investigate the effects of wedge angle, partial penetration, screen location and piezometer location on the steady‐state dimensionless drawdown for different leakage parameters. Results of our study are presented in the form of dimensionless flux‐dimensionless time and dimensionless drawdown‐leakage parameter type curves. The results are useful for evaluating the relative role of lateral wedge boundaries and leakage source on flow in wedge‐shaped aquifers. This is very useful for water management problems and for assessing groundwater pollution. The presented analytical solutions can also be used in parameter identification and in calculating stream depletion rate and volume. Copyright © 2011 John Wiley & Sons, Ltd.     

Bioremediation in Fractured Rock: 1. Modeling to Inform Design,Monitoring, and Expectations

Field characterization of a trichloroethene (TCE) source area in fractured mudstones produced a detailed understanding of the geology, contaminant distribution in fractures and the rock matrix, and hydraulic and transport properties. Groundwater flow and chemical transport modeling that synthesized the field characterization information proved critical for designing bioremediation of the source area. The planned bioremediation involved injecting emulsified vegetable oil and bacteria to enhance the naturally occurring biodegradation of TCE. The flow and transport modeling showed that injection will spread amendments widely over a zone of lower‐permeability fractures, with long residence times expected because of small velocities after injection and sorption of emulsified vegetable oil onto solids. Amendments transported out of this zone will be diluted by groundwater flux from other areas, limiting bioremediation effectiveness downgradient. At nearby pumping wells, further dilution is expected to make bioremediation effects undetectable in the pumped water. The results emphasize that in fracture‐dominated flow regimes, the extent of injected amendments cannot be conceptualized using simple homogeneous models of groundwater flow commonly adopted to design injections in unconsolidated porous media (e.g., radial diverging or dipole flow regimes). Instead, it is important to synthesize site characterization information using a groundwater flow model that includes discrete features representing high‐ and low‐permeability fractures. This type of model accounts for the highly heterogeneous hydraulic conductivity and groundwater fluxes in fractured‐rock aquifers, and facilitates designing injection strategies that target specific volumes of the aquifer and maximize the distribution of amendments over these volumes.     

Estimating Preferential Flow in Karstic Aquifers Using Statistical Mixed Models

Karst aquifers are highly productive groundwater systems often associated with conduit flow. These systems can be highly vulnerable to contamination, resulting in a high potential for contaminant exposure to humans and ecosystems. This work develops statistical models to spatially characterize flow and transport patterns in karstified limestone and determines the effect of aquifer flow rates on these patterns. A laboratory‐scale Geo‐HydroBed model is used to simulate flow and transport processes in a karstic limestone unit. The model consists of stainless steel tanks containing a karstified limestone block collected from a karst aquifer formation in northern Puerto Rico. Experimental work involves making a series of flow and tracer injections, while monitoring hydraulic and tracer response spatially and temporally. Statistical mixed models (SMMs) are applied to hydraulic data to determine likely pathways of preferential flow in the limestone units. The models indicate a highly heterogeneous system with dominant, flow‐dependent preferential flow regions. Results indicate that regions of preferential flow tend to expand at higher groundwater flow rates, suggesting a greater volume of the system being flushed by flowing water at higher rates. Spatial and temporal distribution of tracer concentrations indicates the presence of conduit‐like and diffuse flow transport in the system, supporting the notion of both combined transport mechanisms in the limestone unit. The temporal response of tracer concentrations at different locations in the model coincide with, and confirms the preferential flow distribution generated with the SMMs used in the study.     

Application of Direct Push Methods to Investigate Uranium Distribution in an Alluvial Aquifer

The U.S. EPA 2000 Radionuclide Rule established a maximum contaminant level (MCL) for uranium of 30 µg/L. Many small community water supplies are struggling to comply with this new regulation. At one such community, direct push (DP) methods were applied to obtain hydraulic profiling tool (HPT) logs and install small diameter wells in a section of alluvial deposits located along the Platte River. This work was conducted to evaluate potential sources of elevated uranium in the Clarks, Nebraska drinking water supply. HPT logs were used to understand the hydrostratigraphy of a portion of the aquifer and guide placement of small diameter wells at selected depth intervals. Low-flow sampling of the wells provided water quality parameters and samples for analysis to study the distribution of uranium and variations in aquifer chemistry. Contrary to expectations, the aquifer chemistry revealed that uranium was being mobilized under anoxic and reducing conditions. Review of the test well and new public water supply well construction details revealed that filter packs extended significantly above the screened intervals of the wells. These filter packs were providing a conduit for the movement of groundwater with elevated concentrations of uranium into the supply wells and the community drinking water supply. The methods applied and lessons learned here may be useful for the assessment of unconsolidated aquifers for uranium, arsenic, and many other drinking water supply contaminants.     

Steady state flow with hydraulic conductivity change around large diameter wells

Darcian flow law in aquifers assumes that the aquifer hydraulic conductivity is constant and the groundwater movement is due only to the piezometric level changes through hydraulic gradient. In practice, after the well development the aquifer just around the well has comparatively larger hydraulic conductivity and gradient. Patchy aquifer solutions in the literature consider sudden hydraulic conductivity changes with distance for the steady state flow. The change of transmissivity is demonstrated by the application of slope‐matching procedure to actual field data. It is the main purpose of this paper to derive simple analytical expressions for aquifer parameter evaluations with steadily decreasing hydraulic conductivity around the well. Spatial nonlinear hydraulic conductivity changes around a large‐diameter well within the depression cone of a confined aquifer are considered as exponentially decreasing functions of the radial distance. Copyright © 2010 John Wiley & Sons, Ltd.     

Semi-analytical solutions of groundwater flow in multi-zone (patchy) wedge-shaped aquifers

Alluvial fans are potential sites of potable groundwater in many parts of the world. Characteristics of alluvial fans sediments are changed radially from high energy coarse-grained deposition near the apex to low energy fine-grained deposition downstream so that patchy wedge-shaped aquifers with radial heterogeneity are formed. The hydraulic parameters of the aquifers (e.g. hydraulic conductivity and specific storage) change in the same fashion. Analytical or semi-analytical solutions of the flow in wedge-shaped aquifers are available for homogeneous cases. In this paper we derive semi-analytical solutions of groundwater flow to a well in multi-zone wedge-shaped aquifers. Solutions are provided for three wedge boundary configurations namely: constant head–constant head wedge, constant head–barrier wedge and barrier–barrier wedge. Derivation involves the use of integral transforms methods. The effect of heterogeneity ratios of zones on the response of the aquifer is examined. The results are presented in form of drawdown and drawdown derivative type curves. Heterogeneity has a significant effect on over all response of the pumped aquifer. Solutions help understanding the behavior of heterogeneous multi-zone aquifers for sustainable development of the groundwater resources in alluvial fans.     

Hydraulic head applications of flow logs in the study of heterogeneous aquifers

Permeability profiles derived from high-resolution flow logs in heterogeneous aquifers provide a limited sample of the most permeable beds or fractures determining the hydraulic properties of those aquifers. This paper demonstrates that flow logs can also be used to infer the large-scale properties of aquifers surrounding boreholes. The analysis is based on the interpretation of the hydraulic head values estimated from the flow log analysis. Pairs of quasi-steady flow profiles obtained under ambient conditions and while either pumping or injecting are used to estimate the hydraulic head in each water-producing zone. Although the analysis yields localized estimates of transmissivity for a few water-producing zones, the hydraulic head estimates apply to the far-field aquifers to which these zones are connected. The hydraulic head data are combined with information from other sources to identify the large-scale structure of heterogeneous aquifers. More complicated cross-borehole flow experiments are used to characterize the pattern of connection between large-scale aquifer units inferred from the hydraulic head estimates. The interpretation of hydraulic heads in situ under steady and transient conditions is illustrated by several case studies, including an example with heterogeneous permeable beds in an unconsolidated aquifer, and four examples with heterogeneous distributions of bedding planes and/or fractures in bedrock aquifers.     

Numerical modeling of geothermal groundwater flow in karst aquifer system in eastern Weibei,Shaanxi Province,China

The quantitative assessment of geothermal water resources is important to the exploitation and utilization of geothermal resources. In the geothermal water systems the density of groundwater changes with the temperature, therefore the variations in hydraulic heads and temperatures are very complicated. A three-dimensional density-dependent model coupling the groundwater flow and heat transport is established and used to simulate the geothermal water flow in the karst aquifers in eastern Weibei, Shaanxi Province, China. The multilayered karst aquifer system in the study area is cut by some major faults which control the regional groundwater flow. In order to calibrate and simulate the effect of the major faults, each fault is discretized as a belt of elements with special hydrological parameters in the numerical model. The groundwater dating data are used to be integrated with the groundwater flow pattern and calibrate the model. Simulation results show that the calculated hydraulic heads and temperature fit with the observed data well.