Using temperature to test models of flow near Yucca Mountain,Nevada

Ground water temperatures in the fractured volcanic aquifer near Yucca Mountain, Nevada, have previously been shown to have significant spatial variability with regions of elevated temperatures coinciding roughly with near-vertical north-south trending faults. Using insights gained from one-dimensional models, previous investigators have suggested upwelling along faults from an underlying aquifer as a likely explanation for this ground water temperature pattern. Using a three-dimensional coupled flow and heat-transport model, we show that the thermal high coinciding with the Paintbrush fault zone can be explained without significant upwelling from the underlying aquifer. Instead, the thermal anomaly is consistent with thermal conduction enhanced slightly by vertical ground water movement within the volcanic aquifer sequence. If more than approximately 400 m3/day of water enters the volcanic aquifer from below along a 10 km fault zone, the calculated temperatures at the water table are significantly greater than the measured temperatures. These results illustrate the potential limitations in using one-dimensional models to interpret ground water temperature data, and underscore the value in combining temperature data with fully coupled three-dimensional simulations.     

Influence of ancient thrust faults on the hydrogeology of the Blue Ridge Province

The Blue Ridge Province contains ubiquitous northeast-southwest-trending thrust faults or smaller thrust "slivers" that greatly impact the nature and character of ground water flow in this region. Detailed investigations at a field site in Floyd County, Virginia, indicate that high-permeability zones occur in the brittle crystalline rocks above these thrust faults. Surface and borehole geophysics, aquifer tests, and chlorofluorocarbon and geochemical data reveal that the shallow saprolite aquifer is separated from the deeper fault-zone aquifer by a low-fracture permeability bedrock confining unit, the hydraulic conductivity of which has been estimated to be six orders of magnitude less than the conductivity of the fault zones at the test site. Within the Blue Ridge Province, these fault zones can occur at depths of 300 m or more, can contain a significant amount of storage, and yield significant quantities of water to wells. Furthermore, it is expected that these faults may compartmentalize the deep aquifer system. Recharge to and discharge from the deep aquifer occurs by slow leakage through the confining unit or through localized breach zones that occur where quartz accumulated in high concentrations during metamorphism and later became extensively fractured during episodes of deformation. The results of this investigation stress the importance of thrust faults in this region and suggest that hydrogeologic models for the Blue Ridge Province include these ancient structural features. Faults in crystalline-rock environments may not only influence the hydrology, they may dominate the flow characteristics of a region.     

Delineation of saltwater intrusion zones using the time domain electromagnetic method: the Nabeul–Hammamet coastal aquifer case study (NE Tunisia)

The time domain electromagnetic method (TDEM) is applied to monitor, to delineate and to map the saltwater intrusion zones in the Mediterranean Plio‐Quaternary aquifer. Forty‐two TDEM soundings were carried out in the coastal plain of Nabeul–Hammamet region (NE Tunisia). TDEM resistivity data were correlated with the existing borehole logging data to assign them to a particular lithology and to provide information about the position of the freshwater–seawater transition zone. The geoelectric sections showing the vertical configuration of seawater intrusion, with the brackish‐salty‐saturated zones, have a resistivity ranging from ~0.1 to 5 Ω?m and are detected at a depth lower than 1.5 m. The salinized zones are located at Nabeul (Sidi Moussa, Sidi El Mahrsi, Al Gasba and Mrazgua) and at Hammamet (Touristic zone of Hammamet north and south, Baraket Essahel) and reached a distance of 4 km from the coastline, indicating a severe state for the aquifer in these zones. These TDEM results are confirmed by the increase of chloride concentration content in the analysed water samples of monitoring wells. Moreover, in the northeastern part, the presence of a saltwater front located far from the coast and along the NW–SE major surface fault can be explained by two hypothesis: (i) this fault seems to provide a conduit for seawater to move readily towards the water wells and (ii) the clay and gypsum infiltration of marine Messinian deposits through the fault plane leads to low resistivities. Finally, it comes out from this study that TDEM survey has successfully depicted salinized zones of this coastal aquifer. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Fracture control of ground water flow and water chemistry in a rock aquitard

There are few studies on the hydrogeology of sedimentary rock aquitards although they are important controls in regional ground water flow systems. We formulate and test a three-dimensional (3D) conceptual model of ground water flow and hydrochemistry in a fractured sedimentary rock aquitard to show that flow dynamics within the aquitard are more complex than previously believed. Similar conceptual models, based on regional observations and recently emerging principles of mechanical stratigraphy in heterogeneous sedimentary rocks, have previously been applied only to aquifers, but we show that they are potentially applicable to aquitards. The major elements of this conceptual model, which is based on detailed information from two sites in the Maquoketa Formation in southeastern Wisconsin, include orders of magnitude contrast between hydraulic diffusivity (K/S(s)) of fractured zones and relatively intact aquitard rock matrix, laterally extensive bedding-plane fracture zones extending over distances of over 10 km, very low vertical hydraulic conductivity of thick shale-rich intervals of the aquitard, and a vertical hydraulic head profile controlled by a lateral boundary at the aquitard subcrop, where numerous surface water bodies dominate the shallow aquifer system. Results from a 3D numerical flow model based on this conceptual model are consistent with field observations, which did not fit the typical conceptual model of strictly vertical flow through an aquitard. The 3D flow through an aquitard has implications for predicting ground water flow and for planning and protecting water supplies.     

Hengill,a high temperature thermal area in Iceland

The results of geological, geochemical and isotopic investigations in the Hengill thermal area in Iceland are presented. This area may be regarded as typical of the Icelandic high temperature areas. The Hengill area is mainly built up of pillow lavas and hyaloclastites, which were piled up in subglacial eruptions. Its western half is traversed by a very active NE-SW trending fault zone about 5 km broad, within which are several eruptive fissures of postglacial age. The volcanic rocks are basalts of various kinds, but minor occurrences of intermediate and rhyolitic rocks are also found. Solfataras and fumaroles occur mainly at higher altitude. Within the fault zone they have a linear distribution along faults. Hot water springs are found in the lowlands on the southern border of the thermal area. The concentrations of chlorine, tritium and deuterium in the water indicate that the hot springs on the lowland are fed by a deep lying aquifer which receives water from a recharge area situated in the western branch of the active volcanic zone upto 70 km to the northeast. The hot springs at higher altitude are the outlets of a secondary circulation system fed by local rain water. This water derives its heat partly from the heat content of the rocks penetrated and partly by steam and gas boiled off from the deeper lying aquifer. Preliminary results from other high temperature areas in Iceland indicate the presence of similar hydrological systems.     

Hydrogeological interactions between fault zones and alluvial aquifers in regional flow systems

The hydrological influence of fault zones in tectonic areas is usually difficult to depict from field data. Numerical simulation allows representation of such flow systems and an estimation of flow lines and rates. This paper reports on simulations of the groundwater flow in a range‐and‐basin area affected by a regional fault zone, which may drain or recharge an overlaying alluvial aquifer. Different hydraulic conductivity values for the range rocks, the fault‐zone, and the sedimentary infill of the basin are considered, as well as different fault‐zone widths and boundary conditions. Results show that upward and downward fluxes develop in the upper part of the fault zone, controlled by the action of the alluvial aquifer, influencing the recharge of the sedimentary basin. This paper shows the hydrological efficiency of fault zones as preferential flow; it also analyses the constraints that determine groundwater recharge to the surrounding basins. These results contribute to the understanding of hydrogeological dynamics in tectonic areas. Copyright © 2008 John Wiley & Sons, Ltd.     

Stochastic ground water flow simulation with a fracture zone continuum model

A method is presented for incorporating the hydraulic effects of vertical fracture zones into two-dimensional cell-based continuum models of ground water flow and particle tracking. High hydraulic conductivity features are used in the model to represent fracture zones. For fracture zones that are not coincident with model rows or columns, an adjustment is required for the hydraulic conductivity value entered into the model cells to compensate for the longer flowpath through the model grid. A similar adjustment is also required for simulated travel times through model cells. A travel time error of less than 8% can occur for particles moving through fractures with certain orientations. The fracture zone continuum model uses stochastically generated fracture zone networks and Monte Carlo analysis to quantify uncertainties with simulated advective travel times. An approach is also presented for converting an equivalent continuum model into a fracture zone continuum model by establishing the contribution of matrix block transmissivity to the bulk transmissivity of the aquifer. The methods are used for a case study in west-central Florida to quantify advective travel times from a potential wetland rehydration site to a municipal supply wellfield. Uncertainties in advective travel times are assumed to result from the presence of vertical fracture zones, commonly observed on aerial photographs as photolineaments.     

High-resolution slug testing

The hydraulic conductivity (K) variation has important ramifications for ground water flow and the transport of contaminants in ground water. The delineation of the nature of that variation can be critical to complete characterization of a site and the planning of effective and efficient remedial measures. Site-specific features (such as high-conductivity zones) need to be quantified. Our alluvial field site in the Kansas River valley exhibits spatial variability, very high conductivities, and nonlinear behavior for slug tests in the sand and gravel aquifer. High-resolution, multilevel slug tests have been performed in a number of wells that are fully screened. A general nonlinear model based on the Navier-Stokes equation, nonlinear frictional loss, non-Darcian flow, acceleration effects, radius changes in the wellbore, and a Hvorslev model for the aquifer has been used to analyze the data, employing an automated processing system that runs within the Excel spreadsheet program. It is concluded that slug tests can provide the necessary data to identify the nature of both horizontal and vertical K variation in an aquifer and that improved delineation or higher resolution of K structure is possible with shorter test intervals. The gradation into zones of higher conductivity is sharper than seen previously, and the maximum conductivity observed is greater than previously measured. However, data from this project indicate that well development, the presence of fines, and the antecedent history of the well are important interrelated factors in regard to slug-test response and can prevent obtaining consistent results in some cases.     

Stream water infiltration, bank storage, and storage zone changes due to stream-stage fluctuations

During a flood period, stream-stage increases induce infiltration of stream water into an aquifer; subsequent declines in stream stage cause a reverse motion of the infiltrated water. This paper presents the results of the water exchange rate between a stream and aquifer, the storage volume of the infiltrated stream water in the surrounding aquifer (bank storage), and the storage zone. The storage zone is the part of aquifer where groundwater is replaced by stream water during the flood. MODFLOW was used to simulate stream–aquifer interactions and to quantify rates of stream infiltration and return flow. MODPATH was used to trace the pathlines of the infiltrated stream water and to determine the size of the storage zone. Simulations were focused on the analyses of the effects of the stream-stage fluctuation, aquifer properties, the hydraulic conductivity of streambed sediments, regional hydraulic gradients, and recharge and evapotranspiration (ET) rates on stream–aquifer interactions. Generally, for a given stream–aquifer system, larger flow rates result from larger stream-stage fluctuations; larger storage volumes and storage zones are produced by larger and longer-lasting fluctuations. For a given stream-stage hydrograph, a lower-permeable streambed, an aquitard, or an anisotropic aquifer of low vertical hydraulic conductivity can significantly reduce the rate of infiltration and limit the size of the storage zone. The bank storage solely caused by the stage fluctuation differs slightly between gaining and losing streams. Short-term rainfall recharge and ET loss in the shallow groundwater slightly influence on the flow rate, but their effects on bank storage in a larger area for a longer period can be considerable.     

Efficiency verification of a horizontal flow barrier via flowmeter tests and multilevel sampling

The remediation strategy for an industrial site located in a coastal area involves a pump and treat system and a horizontal flow barrier (HFB) penetrating the main aquifer. To validate the groundwater flow conceptual model and to verify the efficiency of the remediation systems, we carried out piezometric measurements, slug tests, pumping tests, flowmeter tests and multilevel sampling. Flowmeter tests are used to infer vertical groundwater flow directions, and base exchange index is used to infer horizontal flow directions at a metric scale. The selected wells are located both upstream and downstream of the HFB. The installation of the HFB produced constraints to the groundwater flow. A stagnant zone of contaminated freshwater floating over the salt wedge in the upper portion of the aquifer is detected downstream of the HFB. This study confirms that the adopted remediation system is efficiently working in the area upstream of the HFB and even downstream in the bottom part of the aquifer. At the same time, it has also confirmed that hot spots are still present in stagnant zones located downstream of the HFB in the upper part of the aquifer, requiring a different approach to accomplish remediation targets. The integrated approach for flow quantification used in this study allows to discriminate the direction and the magnitude of groundwater fluxes near an HFB in a coastal aquifer. Copyright © 2012 John Wiley & Sons, Ltd.     

Geological structure of hydrothermal fields in the taupo volcanic zone,New Zealand

Wairakei hydrothermal field is underlain by an acid volcanic sequence consisting of the following units: Recent pumice cover, Wairakei Breccia, Huka Falls Formation, Haparangi Rhyolite, Waiora Formation, Waiora Valley Andesite, Wairakei Ignimbrites, Ohakuri Group. The stratified volcanic sequence is draped over a basement horst and thickens eastward and westward into adjoining volcanotectonic depressions. These major depressions have grown progressively during the Quaternary by differential subsidence along active faults and were not the direct result of collapse following major ignimbrite and rhyolite eruptions. The bulk of the steam production is obtained from a thick aquifer of pumice breccias (Waiora aquifer), which is capped by lacustrine shales of the Huka Falls Formation. The best production is obtained from fault zones intersecting the lower part of this aquifer. The Ohakuri Group, underlying the Wairakei Ignimbrites, constitutes a lower, pumice breccia aquifer, not yet directly exploited by drill-holes. Hydrothermal water at 265°C is fed from the lower to the upper aquifer through linear fissures in the Wairakei Ignimbrites, principally at the cress of a small structural dome. These fissures are related to active, north-east striking, normal faults, having a small dextral transcurrent component. Major zones of heat liberation have been localised at the intersection of secondary, north-west cross faults. Fossil, hydrothermal, mud-flow, conglomerates intercalated in the mid-Pleistocene Huka Falls Formation, suggest that hydrothermal activity at Wairakei is at least 500,000 years old. Isopachs on the early Pleistocene Waiora Valley Andesite indicate that the andesite plug and associated flows were crupted at the same fault intersections, now controlling heat flow into the Waiora aquifer. The inferred great age of the hydrothermal system poses problems in maintaining heat flow, and a model of a large, semi-permanent granitic batholith leaking super-critical hydrothermal fluids up active fault zones to heat near-surface aquifers of meteoric water is proposed. Temperature and pressure equilibria of such a magma body may be maintained by gaseous diffusion. Structural and heat flow evidence suggest that the most probable location of the magma body is to the south of the hydrothermal field below a large, late Pleistocene, rhyolitic eruptive centre. Other hydrothermal fields in the Taupo Volcanic Zone — Orakeikorako, Waiotapu and Kawerau — have similar geological sequences and structures to Wairakei. Important features in common include:

1. 1)
   Proximity to a late Pleistocene, rhyolitic eruptive centre typified by ignimbrites, non-welded ash flows and rhyolite domes and sills, which is considered to directly overlie a large, semipermanent magma chamber.     
   
   

Patterns and dynamics of river–aquifer exchange with variably-saturated flow using a fully-coupled model

The parallel physically-based surface–subsurface model PARFLOW was used to investigate the spatial patterns and temporal dynamics of river–aquifer exchange in a heterogeneous alluvial river–aquifer system with deep water table. Aquifer heterogeneity at two scales was incorporated into the model. The architecture of the alluvial hydrofacies was represented based on conditioned geostatistical indicator simulations. Subscale variability of hydraulic conductivities (K) within hydrofacies bodies was created with a parallel Gaussian simulation. The effects of subscale heterogeneity were investigated in a Monte Carlo framework. Dynamics and patterns of river–aquifer exchange were simulated for a 30-day flow event. Simulation results show the rapid formation of saturated connections between the river channel and the deep water table at preferential flow zones that are characterized by high conductivity hydrofacies. Where the river intersects low conductivity hydrofacies shallow perched saturated zones immediately below the river form, but seepage to the deep water table remains unsaturated and seepage rates are low. Preferential flow zones, although only taking up around 50% of the river channel, account for more than 98% of total seepage. Groundwater recharge is most efficiently realized through these zones. Subscale variability of K_sat slightly increased seepage volumes, but did not change the general seepage patterns (preferential flow zones versus perched zones). Overall it is concluded that typical alluvial heterogeneity (hydrofacies architecture) is an important control of river–aquifer exchange in rivers overlying deep water tables. Simulated patterns and dynamics are in line with field observations and results from previous modeling studies using simpler models. Alluvial heterogeneity results in distinct patterns and dynamics of river–aquifer exchange with implications for groundwater recharge and the management of riparian zones (e.g. river channel-floodplain connectivity via saturated zones).     

Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model

In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.     

Hydraulics of a partially penetrating well with skin zone in a confined aquifer

A steady/quasi-steady model is developed for predicting flow into a partially penetrating well with skin zone in a confined aquifer overlying an impervious layer. The model takes into account flow through the bottom of the wellbore, finite skin thickness and finite horizontal and vertical extent of the aquifer. Moreover, the solution can be easily extended to include the mixed-type boundary condition at the well face, where a Dirichlet in the form of a specified hydraulic head and a Neumann in the form of zero flux coexist at the same time at different portions of the well face. The validity of the proposed solution is tested by comparing a few results obtained from the developed model with corresponding results obtained by analytical and numerical means. The study shows that, among other factors remaining constant, both the horizontal and vertical extent of an artesian aquifer, thickness of the skin zone, bottom flow and conductivity contrast of the skin and formation zones, play an important part in deciding flow to a well dug in the aquifer, and hence these factors must be considered while analyzing the problem. The model proposed here can be used to estimate skin thickness as well as hydraulic conductivities of the skin and formation zones of a well with skin zone in an artesian aquifer underlain by an impervious layer by utilizing pumping test data falling in the steady or quasi-steady state of a typical pumping test. As the proposed solution is of a general nature in the sense that it can handle, apart from partial penetration and bottom flow, the finite size skin zone and finite horizontal and vertical extent of an artesian aquifer together with the mixed-type boundary condition at the well face, it is hoped that the predictions coming out of the model will be more realistic than those obtained using solutions developed with more stringent assumptions.     

红河断裂带北段温泉水循环深度与地震活动性的关系探讨

利用SiO_2地热温标计算了红河断裂带北段14个温泉水的热储温度,并结合当地地温梯度计算了温泉水的循环深度,同时统计分析了886～1987年间的地震频度和震源深度。作者从水对断裂及其围岩的弱化的角度讨论了温泉水深循环对断裂强度、断裂摩擦力和地震活动性的影响。结果表明,各亚段地震活动性差异与水对各亚段的弱化程度有关。认为温泉水循环深度最深的洱源亚段相当于“滑移区”,而温泉水循环深度较浅的剑川亚段和弥渡亚段及温泉稀少的大理亚段则相当于“闭锁区”,强震将孕育于温泉水循环深度较深的地区的两端     

The influence of geologic structures on deformation due to ground water withdrawal

A 62 day controlled aquifer test was conducted in thick alluvial deposits at Mesquite, Nevada, for the purpose of monitoring horizontal and vertical surface deformations using a high-precision global positioning system (GPS) network. Initial analysis of the data indicated an anisotropic aquifer system on the basis of the observed radial and tangential deformations. However, new InSAR data seem to indicate that the site may be bounded by an oblique normal fault as the subsidence bowl is both truncated to the northwest and offset from the pumping well to the south. A finite-element numerical model was developed using ABAQUS to evaluate the potential location and hydromechanical properties of the fault based on the observed horizontal deformations. Simulation results indicate that for the magnitude and direction of motion at the pumping well and at other GPS stations, which is toward the southeast (away from the inferred fault), the fault zone (5 m wide) must possess a very high permeability and storage coefficient and cross the study area in a northeast-southwest direction. Simulated horizontal and vertical displacements that include the fault zone closely match observed displacements and indicate the likelihood of the presence of the inferred fault. This analysis shows how monitoring horizontal displacements can provide valuable information about faults, and boundary conditions in general, in evaluating aquifer systems during an aquifer test.     

地震前后井-含水层系统潮汐参数变化特征分析——以云南弥勒井为例

本文对弥勒井2004～2016年井水位同震响应形态及地震前后井水位潮汐参数变化关系特征进行分析。结果表明,在弥勒井水流运动方向总体是以垂向流为主、径向流为辅的前提下,天然因素(地震)和人工因素(洗井、装置改造)均可改变含水层水流运动方向。当震前潮汐因子和相位差反向变化时,能量较大的地震可以改变其水流运动方向,使震后潮汐因子和相位差呈同向变化,水流运动方向由以垂向为主变为以径向为主;能量较小的地震只是改变其含水层渗透系数,没有改变其水流运动方向,震后仍以垂向流为主。当震前潮汐因子和相位差呈同向变化时,能量较大的地震发生可以使径向含水层渗透系数增大,水流运动方向不变,震后仍以径向为主。通过分析潮汐参数变化,可以进一步得知井-含水层径、垂向流渗透系数的变化,为进一步精确计算含水层参数提供了参考。     

The effect of low‐permeability fault zones on groundwater flow in a compartmentalized system. Experimental evidence from a carbonate aquifer (Southern Italy)

The hydrogeological behaviour of fault zones in carbonate aquifers is often neglected in conceptual and numerical models. Furthermore, no information is available regarding the relationships between piezometric levels when significant compartmentalization occurs due to the occurrence of low‐flow fault zones. The aim of this study was to refine the conceptualization of subsurface flow in faulted carbonate aquifers and to analyse relationships between sub‐basins within a compartmentalized aquifer system in Southern Italy. The interactions between compartments that straddle low‐flow faults were investigated over four hydrologic years using a statistical approach to compare (i) the hydraulic heads within two wells located up‐ and down‐gradient of tectonic discontinuities as well as (ii) the rainfall and piezometric levels. The results of this study suggest that a set of barriers exists between the wells, and, therefore, the total head loss observed between the wells (approximately 80 m) should be distributed across several aquitards, with one aquitard exhibiting a relatively high permeability or low degree of integrity. Due to slight differences in permeability, transient conditions in aquitards can occur over relatively short periods, which is in agreement with the results of the statistical data analysis. Consequently, rather than being caused by pure aquitards, aquifer system compartmentalization likely results from slight differences in the permeability between lower‐permeability fault zones and adjacent higher‐permeability protoliths. Copyright © 2014 John Wiley & Sons, Ltd.     

The permeability of the Elkhorn fault zone,South Park,Colorado

The purposes of this study are to use both field and modeling approaches to characterize the permeability of a fault and to assess the role of the fault on regional ground water flow. The study subject is the Elkhorn fault, a low-angle reverse fault that brings Precambrian crystalline rocks over the sediments of Colorado's South Park Basin. The fault is hypothesized to act as a low-permeability barrier to flow, restricting interaction between the crystalline aquifer and the basin sediments. To test this hypothesis and to better predict the permeability structure of the fault, we synthesized geologic data to create a geologic model of the fault, conducted aquifer tests to estimate the hydrogeologic properties of the fault zone, and used ground water modeling to test the influence of a range of hydraulic properties for the fault zone on ground water flow in the region. Our study suggests that the fault is a low-permeability feature. Estimated heads are best matched to observations by modeling the fault as a 10-foot-thick interval of low-permeability fault gouge. Steady-state flow models show that much of the flow in the study area is topographically driven near land surface. Flow rates decrease with depth in the aquifers. In the footwall, ground water moves updip in the Michigan-San Isabel syncline to discharge in the South Park Basin. In the hanging wall, ground water moves east to a regional ground water divide. Sensitivity analyses indicate that hydraulic heads are most sensitive to changes in hydraulic conductivity and recharge.