Science.gov

Sample records for estimate groundwater flow

  1. 3PE: A Tool for Estimating Groundwater Flow Vectors

    EPA Science Inventory

    Evaluation of hydraulic gradients and the associated groundwater flow rates and directions is a fundamental aspect of hydrogeologic characterization. Many methods, ranging in complexity from simple three-point solution techniques to complex numerical models of groundwater flow, ...

  2. Estimation of the velocities and discharge rates of groundwater flows by helium and tritium distributions

    SciTech Connect

    Ivanov, V.V.; Antonenko, E.F.; Obukhova, S.N.

    1995-11-01

    Modeling is used to show that bends in the curves of groundwater flows and confluence of different-age flows demonstrate themselves in spatial variability in concentration fields of {sup 4}He, {sup 3}H, and {sup 14}C isotopes. The results are used to reveal typical flows in a section of Switzerland and to estimate their velocities and discharge rates.

  3. Radiocarbon determinations for estimating groundwater flow velocities in central Florida

    USGS Publications Warehouse

    Hanshaw, B.B.; Back, W.; Rubin, M.

    1965-01-01

    Carbon-14 activity was determined from HCO3- in samples of groundwater obtained from the principal artesian aquifer in Florida. From these data the "age" of water obtained from a series of wells, each progressively farther down gradient on the piezometric surface, was established. Relative carbon-14 ages indicated a velocity of groundwater movement of 23 feet (7 meters) per year for about 85 miles (137 kilometers) of travel. A velocity of 23 feet per year was calculated independently from Darcy's law.

  4. Simultaneous parameter estimation and contaminant source characterization for coupled groundwater flow and contaminant transport modelling

    USGS Publications Warehouse

    Wagner, B.J.

    1992-01-01

    Parameter estimation and contaminant source characterization are key steps in the development of a coupled groundwater flow and contaminant transport simulation model. Here a methodologyfor simultaneous model parameter estimation and source characterization is presented. The parameter estimation/source characterization inverse model combines groundwater flow and contaminant transport simulation with non-linear maximum likelihood estimation to determine optimal estimates of the unknown model parameters and source characteristics based on measurements of hydraulic head and contaminant concentration. First-order uncertainty analysis provides a means for assessing the reliability of the maximum likelihood estimates and evaluating the accuracy and reliability of the flow and transport model predictions. A series of hypothetical examples is presented to demonstrate the ability of the inverse model to solve the combined parameter estimation/source characterization inverse problem. Hydraulic conductivities, effective porosity, longitudinal and transverse dispersivities, boundary flux, and contaminant flux at the source are estimated for a two-dimensional groundwater system. In addition, characterization of the history of contaminant disposal or location of the contaminant source is demonstrated. Finally, the problem of estimating the statistical parameters that describe the errors associated with the head and concentration data is addressed. A stage-wise estimation procedure is used to jointly estimate these statistical parameters along with the unknown model parameters and source characteristics. ?? 1992.

  5. Estimating evapotranspiration and groundwater flow from water-table fluctuations for a general wetland scenario

    USGS Publications Warehouse

    Carlson Mazur, Martha L.; Michael J. Wiley; Douglas A. Wilcox

    2015-01-01

    The use of diurnal water-table fluctuation methods to calculate evapotranspiration (ET) and groundwater flow is of increasing interest in ecohydrological studies. Most studies of this type, however, have been located in riparian wetlands of semi-arid regions where groundwater levels are consistently below topographic surface elevations and precipitation events are infrequent. Current methodologies preclude application to a wider variety of wetland systems. In this study, we extended a method for estimating sub-daily ET and groundwater flow rates from water-level fluctuations to fit highly dynamic, non-riparian wetland scenarios. Modifications included (1) varying the specific yield to account for periodic flooded conditions and (2) relating empirically derived ET to estimated potential ET for days when precipitation events masked the diurnal signal. To demonstrate the utility of this method, we estimated ET and groundwater fluxes over two growing seasons (2006–2007) in 15 wetlands within a ridge-and-swale wetland complex of the Laurentian Great Lakes under flooded and non-flooded conditions. Mean daily ET rates for the sites ranged from 4.0 mm d−1 to 6.6 mm d−1. Shallow groundwater discharge rates resulting from evaporative demand ranged from 2.5 mm d−1 to 4.3 mm d−1. This study helps to expand our understanding of the evapotranspirative demand of plants under various hydrologic and climate conditions.

  6. Water-balance and groundwater-flow estimation for an arid environment: San Diego region, California

    NASA Astrophysics Data System (ADS)

    Flint, L. E.; Flint, A. L.; Stolp, B. J.; Danskin, W. R.

    2012-03-01

    The coastal-plain aquifer that underlies the San Diego City metropolitan area in southern California is a groundwater resource. The understanding of the region-wide water balance and the recharge of water from the high elevation mountains to the east needs to be improved to quantify the subsurface inflows to the coastal plain in order to develop the groundwater as a long term resource. This study is intended to enhance the conceptual understanding of the water balance and related recharge processes in this arid environment by developing a regional model of the San Diego region and all watersheds adjacent or draining to the coastal plain, including the Tijuana River basin. This model was used to quantify the various components of the water balance, including semi-quantitative estimates of subsurface groundwater flow to the coastal plain. Other approaches relying on independent data were used to test or constrain the scoping estimates of recharge and runoff, including a reconnaissance-level groundwater model of the San Diego River basin, one of three main rivers draining to the coastal plain. Estimates of subsurface flow delivered to the coastal plain from the river basins ranged from 12.3 to 28.8 million m3 yr-1 from the San Diego River basin for the calibration period (1982-2009) to 48.8 million m3 yr-1 from all major river basins for the entire coastal plain for the long-term period 1940-2009. This range of scoping estimates represents the impact of climatic variability and realistically bounds the likely groundwater availability, while falling well within the variable estimates of regional recharge. However, the scarcity of physical and hydrologic data in this region hinders the exercise to narrow the range and reduce the uncertainty.

  7. Estimating groundwater recharge through unsaturated flow modelling: Sensitivity to boundary conditions and vertical discretization

    NASA Astrophysics Data System (ADS)

    Carrera-Hernández, J. J.; Smerdon, B. D.; Mendoza, C. A.

    2012-07-01

    SummaryUnsaturated flow modelling is increasingly being used to estimate potential groundwater recharge. A review of previous studies found that unit-gradient and fixed water table lower boundary conditions have been applied to models of both constant and variable vertical grid spacing (discretization). In order to provide a general guidance, this work studies the effect of both discretization and boundary conditions on simulation times and estimated fluxes at the water table, using one-dimensional models of 2, 4, 6, and 12 m comprised of sand, sandy loam, loamy sand, and loam. The study uses climatological data from the Boreal Plain of northern Alberta, Canada. Because of the long-term average water deficit and the thick unconsolidated glacial deposits, unsaturated flow is expected to be vertical, both downward and upward, and inter-annual changes in water storage will be important. Long-term simulations (1919-2007) that comprised both wet and dry cycles, reveal that when a variable vertical discretization at both the top and bottom of the columns (varying from 0.1 to 10 cm) is utilized, a balance between simulation accuracy and running time can be achieved. It is also found that whenever the unsaturated flow modelling approach is used to estimate potential groundwater recharge, a fixed-head lower boundary condition should be selected because it also allows upward flux from the water table during dry periods, a situation that prevails on both sub-humid and semi-arid areas, where accurate groundwater recharge estimates are needed the most. However, it should be kept in mind that the use of a fixed water table is a simple representation of the regional water table, which in reality interacts with the regional groundwater flow and surface water bodies (e.g., lakes and wetlands).

  8. Estimation of regional-scale groundwater flow properties in the Bengal Basin of India and Bangladesh

    USGS Publications Warehouse

    Michael, H.A.; Voss, C.I.

    2009-01-01

    Quantitative evaluation of management strategies for long-term supply of safe groundwater for drinking from the Bengal Basin aquifer (India and Bangladesh) requires estimation of the large-scale hydrogeologic properties that control flow. The Basin consists of a stratified, heterogeneous sequence of sediments with aquitards that may separate aquifers locally, but evidence does not support existence of regional confining units. Considered at a large scale, the Basin may be aptly described as a single aquifer with higher horizontal than vertical hydraulic conductivity. Though data are sparse, estimation of regional-scale aquifer properties is possible from three existing data types: hydraulic heads, 14C concentrations, and driller logs. Estimation is carried out with inverse groundwater modeling using measured heads, by model calibration using estimated water ages based on 14C, and by statistical analysis of driller logs. Similar estimates of hydraulic conductivities result from all three data types; a resulting typical value of vertical anisotropy (ratio of horizontal to vertical conductivity) is 104. The vertical anisotropy estimate is supported by simulation of flow through geostatistical fields consistent with driller log data. The high estimated value of vertical anisotropy in hydraulic conductivity indicates that even disconnected aquitards, if numerous, can strongly control the equivalent hydraulic parameters of an aquifer system. ?? US Government 2009.

  9. Estimation of the pore pressure distribution from three dimensional groundwater flow model at mine sites in Korea

    NASA Astrophysics Data System (ADS)

    Kang, Sangsoo; Jang, Myounghwan; Kim, Gyoungman; Kim, Donghui; Kim, Daehoon; Baek, Hwanjo

    2016-04-01

    Mining activities continually change the groundwater flow and associated pore pressure distributions within the rockmass around the mine openings or the open-pit bench during the operational periods. As the pore pressure distributions may substantially affect the mechanical behaviour or stability of the rockmass, it is important to monitor the variation of pore pressure incurred by mining operation. The pore pressure distributions within the rockmass can be derived using a two- or three-dimensional finite element groundwater flow model, adopted to simulate the groundwater flow. While the groundwater inflow at mines has generally been dealt with respect to the working environment, detailed case studies on the distribution of pore water pressure related to the stability analysis of mine openings have been relatively rare in Korea. Recently, however, as the health and safety problems are emerged for sustainable mining practice, these issues are of the major concerns for the mining industries. This study aims to establish a three dimensional groundwater flow model to estimate the pore pressure distributions in order to employ as an input parameter for numerical codes such as the FLAC 3D. Also, the groundwater flow simulated can be used for de-watering design at a mine site. The MINEDW code, a groundwater flow model code specifically developed to simulate the complicated hydro-geologic conditions related to mining, has mainly been used in this study. Based on the data collected from field surveys and literature reviews, a conceptual model was established and sensitivity analysis was performed.

  10. Relative efficiency of four parameter-estimation methods in steady-state and transient ground-water flow models

    USGS Publications Warehouse

    Hill, M.C.

    1990-01-01

    Parameters in numerical ground-water flow models have been successfully estimated using nonlinear-optimization methods such as the modified Gauss-Newton (GN) method and conjugate-direction methods. This paper investigates the relative efficiency of GN and three conjugate-direction parameter-estimation methods on two-dimensional, steady-state and transient ground-water flow test cases. The steady-state test cases are included to compare the performance of the algorithm with published examples. The three conjugate-direction methods are the Fletcher-Reeves (FR) and quasi-Newton (QN) regression methods, and combination Fletcher-Reeves quasi-Newton (FR-QN). All three are combined with Newton's method of calculating step size. The numerical ground-water flow model is described by McDonald and Harbaugh.

  11. Estimated Ground-Water Withdrawals from the Death Valley Regional Flow System, Nevada and California, 1913-98

    USGS Publications Warehouse

    Moreo, Michael T.; Halford, Keith J.; La Camera, Richard J.; Laczniak, Randell J.

    2003-01-01

    Ground-water withdrawals from 1913 through 1998 from the Death Valley regional flow system have been compiled to support a regional, three-dimensional, transient ground-water flow model. Withdrawal locations and depths of production intervals were estimated and associated errors were reported for 9,300 wells. Withdrawals were grouped into three categories: mining, public-supply, and commercial water use; domestic water use; and irrigation water use. In this report, groupings were based on the method used to estimate pumpage. Cumulative ground-water withdrawals from 1913 through 1998 totaled 3 million acre-feet, most of which was used to irrigate alfalfa. Annual withdrawal for irrigation ranged from 80 to almost 100 percent of the total pumpage. About 75,000 acre-feet was withdrawn for irrigation in 1998. Annual irrigation withdrawals generally were estimated as the product of irrigated acreage and application rate. About 320 fields totaling 11,000 acres were identified in six hydrographic areas. Annual application rates for high water-use crops ranged from 5 feet in Penoyer Valley to 9 feet in Pahrump Valley. The uncertainty in the estimates of ground-water withdrawals was attributed primarily to the uncertainty of application rate estimates. Annual ground-water withdrawal was estimated at about 90,000 acre-feet in 1998 with an assigned uncertainty bounded by 60,000 to 130,000 acre-feet.

  12. Estimated Ground-water Withdrawals From the Death Valley Regional Flow System, Nevada and California, 1913-98

    SciTech Connect

    M.T. Moreo; K.J. Halford; R.J. LaCamera; and R.J. Laczniak

    2003-09-30

    Ground-water withdrawals from 1913 through 1998 from the Death Valley regional flow system have been compiled to support a regional,three-dimensional, transient ground-water flow model. Withdrawal locations and depths of production intervals were estimated and associated errors were reported for 9,300 wells. Withdrawals were grouped into three categories: mining, public-supply, and commercial water use; domestic water use; and irrigation water use. In this report, groupings were based on the method used to estimate pumpage. Cumulative ground-water withdrawals from 1913 through 1998 totaled 3 million acre-feet, most of which was used to irrigate alfalfa. Annual withdrawal for irrigation ranged from 80 to almost 100 percent of the total pumpage. About 75,000 acre-feet was withdrawn for irrigation in 1998. Annual irrigation withdrawals generally were estimated as the product of irrigated acreage and application rate. About 320 fields totaling 11,000 acres were identified in six hydrographic areas. Annual application rates for high water-use crops ranged from 5 feet in Penoyer Valley to 9 feet in Pahrump Valley. The uncertainty in the estimates of ground-water withdrawals was attributed primarily to the uncertainty of application rate estimates. Annual ground-water withdrawal was estimated at about 90,000 acre-feet in 1998 with an assigned uncertainty bounded by 60,000 to 130,000 acre-feet.

  13. Estimates of ground-water recharge, base flow, and stream reach gains and losses in the Willamette River basin, Oregon

    USGS Publications Warehouse

    Lee, Karl K.; Risley, John C.

    2002-01-01

    Precipitation-runoff models, base-flow-separation techniques, and stream gain-loss measurements were used to study recharge and ground-water surface-water interaction as part of a study of the ground-water resources of the Willamette River Basin. The study was a cooperative effort between the U.S. Geological Survey and the State of Oregon Water Resources Department. Precipitation-runoff models were used to estimate the water budget of 216 subbasins in the Willamette River Basin. The models were also used to compute long-term average recharge and base flow. Recharge and base-flow estimates will be used as input to a regional ground-water flow model, within the same study. Recharge and base-flow estimates were made using daily streamflow records. Recharge estimates were made at 16 streamflow-gaging-station locations and were compared to recharge estimates from the precipitation-runoff models. Base-flow separation methods were used to identify the base-flow component of streamflow at 52 currently operated and discontinued streamflow-gaging-station locations. Stream gain-loss measurements were made on the Middle Fork Willamette, Willamette, South Yamhill, Pudding, and South Santiam Rivers, and were used to identify and quantify gaining and losing stream reaches both spatially and temporally. These measurements provide further understanding of ground-water/surface-water interactions.

  14. Hydraulic-property estimates for use with a transient ground-water flow model of the Death Valley regional ground-water flow system, Nevada and California

    SciTech Connect

    W.R. Belcher; P.E. Elliott; A.L. Geldon

    2001-12-31

    The Death Valley regional ground-water flow system encompasses an area of about 43,500 square kilometers in southeastern California and southern Nevada. The study area is underlain by Quaternary to Tertiary basin-fill sediments and mafic-lava flows; Tertiary volcanic, volcaniclastic, and sedimentary rocks; Tertiary to Jurassic granitic rocks; Triassic to Middle Proterozoic carbonate and clastic sedimentary rocks; and Early Proterozoic igneous and metamorphic rocks. The rock assemblage in the Death Valley region is extensively faulted as a result of several episodes of tectonic activity. This study is comprised of published and unpublished estimates of transmissivity, hydraulic conductivity, storage coefficient, and anisotropy ratios for hydrogeologic units within the Death Valley region study area. Hydrogeologic units previously proposed for the Death Valley regional transient ground-water flow model, were recognized for the purpose of studying the distribution of hydraulic properties. Analyses of regression and covariance were used to assess if a relation existed between hydraulic conductivity and depth for most hydrogeologic units. Those analyses showed a weak, quantitatively indeterminate, relation between hydraulic conductivity and depth.

  15. Hydraulic-property estimates for use with a transient ground-water flow model of the Death Valley regional ground-water flow system, Nevada and California

    USGS Publications Warehouse

    Belcher, Wayne R.; Elliott, Peggy E.; Geldon, Arthur L.

    2001-01-01

    The Death Valley regional ground-water flow system encompasses an area of about 43,500 square kilometers in southeastern California and southern Nevada, between latitudes 35? and 38?15' north and longitudes 115? and 117?45' west. The study area is underlain by Quaternary to Tertiary basin-fill sediments and mafic-lava flows; Tertiary volcanic, volcaniclastic, and sedimentary rocks; Tertiary to Jurassic granitic rocks; Triassic to Middle Proterozoic carbonate and clastic sedimentary rocks; and Early Proterozoic igneous and metamorphic rocks. The rock assemblage in the Death Valley region is extensively faulted as a result of several episodes of tectonic activity. This study is comprised of published and unpublished estimates of transmissivity, hydraulic conductivity, storage coefficient, and anisotropy ratios for hydrogeologic units within the Death Valley region study area. Hydrogeologic units previously proposed for the Death Valley regional transient ground-water flow model were recognized for the purpose of studying the distribution of hydraulic properties. Analyses of regression and covariance were used to assess if a relation existed between hydraulic conductivity and depth for most hydrogeologic units. Those analyses showed a weak, quantitatively indeterminate, relation between hydraulic conductivity and depth.

  16. Death Valley regional groundwater flow model calibration using optimal parameter estimation methods and geoscientific information systems

    USGS Publications Warehouse

    D'Agnese, F. A.; Faunt, C.C.; Hill, M.C.; Turner, A.K.

    1996-01-01

    A three-layer Death Valley regional groundwater flow model was constructed to evaluate potential regional groundwater flow paths in the vicinity of Yucca Mountain, Nevada. Geoscientific information systems were used to characterize the complex surface and subsurface hydrogeological conditions of the area, and this characterization was used to construct likely conceptual models of the flow system. The high contrasts and abrupt contacts of the different hydrogeological units in the subsurface make zonation the logical choice for representing the hydraulic conductivity distribution. Hydraulic head and spring flow data were used to test different conceptual models by using nonlinear regression to determine parameter values that currently provide the best match between the measured and simulated heads and flows.

  17. Scarce water resources and scarce data: Estimating recharge for a complex 3D groundwater flow model in arid regions

    NASA Astrophysics Data System (ADS)

    Gräbe, A. C.; Guttman, J.; Rödiger, T.; Siebert, C.; Merz, R.; Kolditz, O.

    2012-12-01

    Semi-arid to arid regions are usually characterized by a scarcity of precipitation and a lack of stream flow. Especially in desert environments, groundwater is one of the most important fresh water sources and its recharge is basically controlled by two main mechanisms: the direct regional infiltration of precipitation in the mountains and interdrainage areas in the first place and secondly the flood water infiltration through ephemeral channel beds (transmission loss). Due to extensive spatio-temporal data scarcity, direct quantitative estimations of groundwater recharge are often difficult to perform, and numerical models simulating the water fluxes, have to be applied to enable a quantitative approximation of the groundwater recharge. We made an assumption about the quantity of recharge for the subsurface catchment of the western Dead Sea escarpment, which is at the same time the input for the complex groundwater flow model of the Judea Group Aquifer. This can only be suggested if the hydrogeological situation in the tectonically complex region is fully understood. A number of simplified models of the Judea Group aquifer have been formulated and employed using a two-dimensional (one horizontal layered) numerical simulation of groundwater flow (Baida et al. 1978; Goldschtoff & Shachnai, 1980; Guttman, 2000; Laronne Ben-Itzhak & Gvirtzmann, 2005). However, all previous approaches focused only on a limited area of the Judea Group aquifer. We developed a high resolution regional groundwater flow model for the entire western basin of the Dead Sea. Whereas the structural model could be defined using a large geological dataset, the challenge was to generate the groundwater flow model with only limited well data. With the help of the scientific software OpenGeoSys (OGS) the challenge was reliably solved resulting in a simulation of the hydraulic characteristics (hydraulic conductivity and hydraulic head) of the cretaceous aquifer system, which was calibrated using PEST.

  18. Estimation of groundwater consumption by phreatophytes using diurnal water table fluctuations: A saturated-unsaturated flow assessment

    USGS Publications Warehouse

    Loheide, S.P., II; Butler, J.J., Jr.; Gorelick, S.M.

    2005-01-01

    Groundwater consumption by phreatophytes is a difficult-to-measure but important component of the water budget in many arid and semiarid environments. Over the past 70 years the consumptive use of groundwater by phreatophytes has been estimated using a method that analyzes diurnal trends in hydrographs from wells that are screened across the water table (White, 1932). The reliability of estimates obtained with this approach has never been rigorously evaluated using saturated-unsaturated flow simulation. We present such an evaluation for common flow geometries and a range of hydraulic properties. Results indicate that the major source of error in the White method is the uncertainty in the estimate of specific yield. Evapotranspirative consumption of groundwater will often be significantly overpredicted with the White method if the effects of drainage time and the depth to the water table on specific yield are ignored. We utilize the concept of readily available specific yield as the basis for estimation of the specific yield value appropriate for use with the White method. Guidelines are defined for estimating readily available specific yield based on sediment texture. Use of these guidelines with the White method should enable the evapotranspirative consumption of groundwater to be more accurately quantified. Copyright 2005 by the American Geophysical Union.

  19. Ground-water discharge determined from estimates of evapotranspiration, Death Valley regional flow system, Nevada and California

    USGS Publications Warehouse

    Laczniak, Randell J.; Smith, J. LaRue; Elliott, Peggy E.; DeMeo, Guy A.; Chatigny, Melissa A.; Roemer, Gaius J.

    2001-01-01

    The Death Valley regional flow system (DVRFS) is one of the larger ground-water flow systems in the southwestern United States and includes much of southern Nevada and the Death Valley region of eastern California. Centrally located within the ground-water flow system is the Nevada Test Site (NTS). The NTS, a large tract covering about 1,375 square miles, historically has been used for testing nuclear devices and currently is being studied as a potential repository for the long-term storage of high-level nuclear waste generated in the United States. The U.S. Department of Energy, as mandated by Federal and State regulators, is evaluating the risk associated with contaminants that have been or may be introduced into the subsurface as a consequence of any past or future activities at the NTS. Because subsurface contaminants can be transported away from the NTS by ground water, components of the ground-water budget are of great interest. One such component is regional ground-water discharge. Most of the ground water leaving the DVRFS is limited to local areas where geologic and hydrologic conditions force ground water upward toward the surface to discharge at springs and seeps. Available estimates of ground-water discharge are based primarily on early work done as part of regional reconnaissance studies. These early efforts covered large, geologically complex areas and often applied substantially different techniques to estimate ground-water discharge. This report describes the results of a study that provides more consistent, accurate, and scientifically defensible measures of regional ground-water losses from each of the major discharge areas of the DVRFS. Estimates of ground-water discharge presented in this report are based on a rigorous quantification of local evapotranspiration (ET). The study identifies areas of ongoing ground-water ET, delineates different ET areas based on similarities in vegetation and soil-moisture conditions, and determines an ET rate for

  20. Scarce data in hydrology and hydrogeology: Estimation and modelling of groundwater recharge for a numerical groundwater flow model in a semi-arid to arid catchment

    NASA Astrophysics Data System (ADS)

    Gräbe, Agnes; Schulz, Stephan; Rödiger, Tino; Kolditz, Olaf

    2013-04-01

    Water resources are strongly limited in semi-arid to arid regions and groundwater constitutes often the only possibility for fresh water for the population and industry. An understanding of the hydrological processes and the estimation of magnitude of water balance parameters also includes the knowledge of processes of groundwater recharge. For the sustainable management of water resources, it is essential to estimate the potential groundwater recharge under the given climatic conditions. We would like to present the results of a hydrological model, which is based on the HRU- concept and intersected the parameters of climatic conditions, topography, geology, soil, vegetation and land use to calculate the groundwater recharge. This model was primarily developed for humid area applications and has now been adapted to the regional conditions in the semi-arid to arid region. It was quite a challenge to understand the hydrological processes in the semi-arid to arid study area and to implement those findings (e.g. routing [Schulz (in prep.)]) into the model structure. Thus we compared the existing approaches for groundwater recharge estimations (chloride mass balance [Marei et. al 2000], empirical relations such as rainfall and base flow-relation [Goldschmidt 1960; Guttman 2000; Hughes 2008; Issar 1993; Lerner 1990; De Vries et. al 2002]) with the results of our numerical model. References: De Vries, J. J., I. Simmers (2002): Groundwater recharge: an overview of processes and challenges. Hydrogeology Journal (2002) 10: 5-17. DOI 10.1007/s10040-001-0171-7. Guttman, J., 2000. Multi-Lateral Project B: Hydrogeology of the Eastern Aquifer in the Judea Hills and Jordan Valley. Mekorot Water Company, Report 468, p. 36. Hughes, A. G., M. M. Mansour, N. S. Robins (2008): Evaluation of distributed recharge in an upland semi-arid karst system: the West Bank Mountain Aquifer, Middle East. Hydrogeology Journal (2008) 16: 845-854. DOI 10.1007/s10040-008-0273-6 Issar, A. S. (1993

  1. U-234/U-238 ratio: Qualitative estimate of groundwater flow in Rocky Flats monitoring wells

    SciTech Connect

    Laul, J.C.

    1994-02-01

    Groundwater movement through various pathways is the primary mechanism for the transport of radionuclides and trace elements in a water/rock interaction. About three dozen wells, installed in the Rocky Flats Plant (RFP) Solar Evaporation Ponds (SEP) area, are monitored quarterly to evaluate the extent of any lateral and downgradient migration of contaminants from the Solar Evaporation Ponds: 207-A; 207-B North, 207-B Center, and 207-B South; and 207-C. The Solar Ponds are the main source for the various contaminants: radionuclides (U-238, U-234, Pu-239, 240 and Am-241); anions; and trace metals to groundwaters. The U-238 concentrations in Rocky Flats groundwaters vary from <0.2 to 69 pCi/I (IpCi = 3 ug). However, the activity U-234/U-238 ratios are low and range mostly 1.2 to 2.7. The low activity ratios can be interpreted to suggest that the groundwaters are moving slow (estimate of about a few ft/year based on hydrologic parameters [1]. Uranium is probably present in the +6 state, predominantly as a uranyl carbonate complexes UO{sub 2}(CO{sub 3}){sub 2}{sup 2{minus}}, because of the predominant bicarbonate medium.

  2. Estimation of groundwater use for a groundwater-flow model of the Lake Michigan Basin and adjacent areas, 1864-2005

    USGS Publications Warehouse

    Buchwald, Cheryl A.; Luukkonen, Carol L.; Rachol, Cynthia M.

    2010-01-01

    The U.S. Geological Survey, at the request of Congress, is assessing the availability and use of the Nation's water resources to help characterize how much water is available now, how water availability is changing, and how much water can be expected to be available in the future. The Great Lakes Basin Pilot project of the U.S. Geological Survey national assessment of water availability and use focused on the Great Lakes Basin and included detailed studies of the processes governing water availability in the Great Lakes Basin. One of these studies included the development of a groundwater-flow model of the Lake Michigan Basin. This report describes the compilation and estimation of the groundwater withdrawals in those areas in Wisconsin, Michigan, Indiana, and Illinois that were needed for the Lake Michigan Basin study groundwater-flow model. These data were aggregated for 12 model time intervals spanning 1864 to 2005 and were summarized by model area, model subregion, category of water use, aquifer system, aquifer type, and hydrogeologic unit model layer. The types and availability of information on groundwater withdrawals vary considerably among states because water-use programs often differ in the types of data collected and in the methods and frequency of data collection. As a consequence, the methods used to estimate and verify the data also vary. Additionally, because of the different sources of data and different terminologies applied for the purposes of this report, the water-use data published in this report may differ from water-use data presented in other reports. These data represent only a partial estimate of groundwater use in each state because estimates were compiled only for areas in Wisconsin, Michigan, Indiana, and Illinois within the Lake Michigan Basin model area. Groundwater-withdrawal data were compiled for both nearfield and farfield model areas in Wisconsin and Illinois, whereas these data were compiled primarily for the nearfield model

  3. Global scale groundwater flow model

    NASA Astrophysics Data System (ADS)

    Sutanudjaja, Edwin; de Graaf, Inge; van Beek, Ludovicus; Bierkens, Marc

    2013-04-01

    As the world's largest accessible source of freshwater, groundwater plays vital role in satisfying the basic needs of human society. It serves as a primary source of drinking water and supplies water for agricultural and industrial activities. During times of drought, groundwater sustains water flows in streams, rivers, lakes and wetlands, and thus supports ecosystem habitat and biodiversity, while its large natural storage provides a buffer against water shortages. Yet, the current generation of global scale hydrological models does not include a groundwater flow component that is a crucial part of the hydrological cycle and allows the simulation of groundwater head dynamics. In this study we present a steady-state MODFLOW (McDonald and Harbaugh, 1988) groundwater model on the global scale at 5 arc-minutes resolution. Aquifer schematization and properties of this groundwater model were developed from available global lithological model (e.g. Dürr et al., 2005; Gleeson et al., 2010; Hartmann and Moorsdorff, in press). We force the groundwtaer model with the output from the large-scale hydrological model PCR-GLOBWB (van Beek et al., 2011), specifically the long term net groundwater recharge and average surface water levels derived from routed channel discharge. We validated calculated groundwater heads and depths with available head observations, from different regions, including the North and South America and Western Europe. Our results show that it is feasible to build a relatively simple global scale groundwater model using existing information, and estimate water table depths within acceptable accuracy in many parts of the world.

  4. Ground-water system, estimation of aquifer hydraulic properties, and effects of pumping on ground-water flow in Triassic sedimentary rocks in and near Lansdale, Pennsylvania

    USGS Publications Warehouse

    Senior, Lisa A.; Goode, Daniel J.

    1999-01-01

    areas in three drainages, the Wissahickon, Towamencin, and Neshaminy Creeks.Ground-water flow was simulated for different pumping patterns representing past and current conditions. The three-dimensional numerical flow model (MODFLOW) was automatically calibrated by use of a parameter estimation program (MODFLOWP). Steady-state conditions were assumed for the calibration period of 1996. Model calibration indicates that estimated recharge is 8.2 inches (208 millimeters) and the regional anisotropy ratio for the sedimentary-rock aquifer is about 11 to 1, with permeability greatest along strike. The regional anisotropy is caused by up- and down-dip termination of high-permeability bed-oriented features, which were not explicitly simulated in the regional-scale model. The calibrated flow model was used to compare flow directions and capture zones in Lansdale for conditions corresponding to relatively high pumping rates in 1994 and to lower pumping rates in 1997. Comparison of the 1994 and 1997 simulations indicates that wells pumped at the lower 1997 rates captured less ground water from known sites of contamination than wells pumped at the 1994 rates. Ground-water flow rates away from Lansdale increased as pumpage decreased in 1997.A preliminary evaluation of the relation between ground-water chemistry and conditions favorable for the degradation of chlorinated solvents was based on measurements of dissolved-oxygen concentration and other chemical constituents in water samples from 92 wells. About 18 percent of the samples contained less than or equal to 5 milligrams per liter dissolved oxygen, a concentration that indicates reducing conditions favorable for degradation of chlorinated solvents.

  5. U.S. Geological Survey groundwater toolbox, a graphical and mapping interface for analysis of hydrologic data (version 1.0): user guide for estimation of base flow, runoff, and groundwater recharge from streamflow data

    USGS Publications Warehouse

    Barlow, Paul M.; Cunningham, William L.; Zhai, Tong; Gray, Mark

    2015-01-01

    This report is a user guide for the streamflow-hydrograph analysis methods provided with version 1.0 of the U.S. Geological Survey (USGS) Groundwater Toolbox computer program. These include six hydrograph-separation methods to determine the groundwater-discharge (base-flow) and surface-runoff components of streamflow—the Base-Flow Index (BFI; Standard and Modified), HYSEP (Fixed Interval, Sliding Interval, and Local Minimum), and PART methods—and the RORA recession-curve displacement method and associated RECESS program to estimate groundwater recharge from streamflow data. The Groundwater Toolbox is a customized interface built on the nonproprietary, open source MapWindow geographic information system software. The program provides graphing, mapping, and analysis capabilities in a Microsoft Windows computing environment. In addition to the four hydrograph-analysis methods, the Groundwater Toolbox allows for the retrieval of hydrologic time-series data (streamflow, groundwater levels, and precipitation) from the USGS National Water Information System, downloading of a suite of preprocessed geographic information system coverages and meteorological data from the National Oceanic and Atmospheric Administration National Climatic Data Center, and analysis of data with several preprocessing and postprocessing utilities. With its data retrieval and analysis tools, the Groundwater Toolbox provides methods to estimate many of the components of the water budget for a hydrologic basin, including precipitation; streamflow; base flow; runoff; groundwater recharge; and total, groundwater, and near-surface evapotranspiration.

  6. Groundwater flow and transport modeling

    USGS Publications Warehouse

    Konikow, L.F.; Mercer, J.W.

    1988-01-01

    Deterministic, distributed-parameter, numerical simulation models for analyzing groundwater flow and transport problems have come to be used almost routinely during the past decade. A review of the theoretical basis and practical use of groundwater flow and solute transport models is used to illustrate the state-of-the-art. Because of errors and uncertainty in defining model parameters, models must be calibrated to obtain a best estimate of the parameters. For flow modeling, data generally are sufficient to allow calibration. For solute-transport modeling, lack of data not only limits calibration, but also causes uncertainty in process description. Where data are available, model reliability should be assessed on the basis of sensitivity tests and measures of goodness-of-fit. Some of these concepts are demonstrated by using two case histories. ?? 1988.

  7. Groundwater-flow parameter estimation and quality modeling of the Equus Beds aquifer in Kansas, U.S.A.

    USGS Publications Warehouse

    Sophocleous, M.A.

    1984-01-01

    The salinity problems created in the Burrton area as a result of poor oil-field brine disposal practices of the past continue to be a major concern to the area depending on the Equus Beds aquifer for water, including the City of Wichita, Kansas. In this paper, an attempt is made to predict where and how fast the brine plume will move in this area, and what the average chloride concentrations in different parts of the aquifer are. In order to make such predictions, it was necessary to get a calibrated model of the groundwater-flow velocity field. Multiple regression analysis is used for parameter estimation of the steady-state groundwater-flow equation applied in the most critical area of the Equus Beds aquifer. Results of such an analysis produced a correlation coefficient of 0.992 between calculated and observed values of hydraulic head. A chloride transport modeling effort is then carried out despite some serious data deficiencies, the significance of which are evaluated through sensitivity analysis. Thus, starting with the quasi steady-state conditions of the early 1940's, it was possible to match the present chloride distribution satisfactorily. Chloride concentration predictions made for the year 2000 indicate that the quality of the Wichita well-field waters will not generally deteriorate from their present condition by that time. ?? 1984.

  8. Oahu Groundwater Flow Model

    DOE Data Explorer

    Nicole Lautze

    2015-01-01

    Groundwater flow model for the island of Oahu. Data is from the following sources: Rotzoll, K., A.I. El-Kadi. 2007. Numerical Ground-Water Flow Simulation for Red Hill Fuel Storage Facilities, NAVFAC Pacific, Oahu, Hawaii - Prepared TEC, Inc. Water Resources Research Center, University of Hawaii, Honolulu.; Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume VII – Island of Oahu Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008.; and Whittier, R. and A.I. El-Kadi. 2009. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. December 2009.

  9. Global estimates of fresh submarine groundwater discharge

    NASA Astrophysics Data System (ADS)

    Luijendijk, Elco; Gleeson, Tom; Moosdorf, Nils

    2016-04-01

    Fresh submarine groundwater discharge, the flow of fresh groundwater to oceans, may be a significant contributor to the water and chemical budgets of the world's oceans. We present new estimates of the flux of fresh groundwater to the world's oceans. We couple density-dependent numerical simulations of generic models of coastal basins with geospatial databases of hydrogeological parameters and topography to resolve the rate of terrestrially-derived submarine groundwater discharge globally. We compare the model results to a new global compilation of submarine groundwater discharge observations. The results show that terrestrially-derived SGD is highly sensitive to permeability. In most watersheds only a small fraction of groundwater recharge contributes to submarine groundwater discharge, with most recharge instead contributing to terrestrial discharge in the form of baseflow or evapotranspiration. Fresh submarine groundwater discharge is only significant in watersheds that contain highly permeable sediments, such as coarse-grained siliciclastic sediments, karstic carbonates or volcanic deposits. Our estimates of global submarine groundwater discharge are much lower than most previous estimates. However, many tropical and volcanic islands are hotspots of submarine groundwater discharge and solute fluxes towards the oceans. The comparison of model results and data highlights the spatial variability of SGD and the difficulty of scaling up observations.

  10. Ground-water pumpage and artificial recharge estimates for calendar year 2000 and average annual natural recharge and interbasin flow by hydrographic area, Nevada

    USGS Publications Warehouse

    Lopes, Thomas J.; Evetts, David M.

    2004-01-01

    Nevada's reliance on ground-water resources has increased because of increased development and surface-water resources being fully appropriated. The need to accurately quantify Nevada's water resources and water use is more critical than ever to meet future demands. Estimated ground-water pumpage, artificial and natural recharge, and interbasin flow can be used to help evaluate stresses on aquifer systems. In this report, estimates of ground-water pumpage and artificial recharge during calendar year 2000 were made using data from a variety of sources, such as reported estimates and estimates made using Landsat satellite imagery. Average annual natural recharge and interbasin flow were compiled from published reports. An estimated 1,427,100 acre-feet of ground water was pumped in Nevada during calendar year 2000. This total was calculated by summing six categories of ground-water pumpage, based on water use. Total artificial recharge during 2000 was about 145,970 acre-feet. At least one estimate of natural recharge was available for 209 of the 232 hydrographic areas (HAs). Natural recharge for the 209 HAs ranges from 1,793,420 to 2,583,150 acre-feet. Estimates of interbasin flow were available for 151 HAs. The categories and their percentage of the total ground-water pumpage are irrigation and stock watering (47 percent), mining (26 percent), water systems (14 percent), geothermal production (8 percent), self-supplied domestic (4 percent), and miscellaneous (less than 1 percent). Pumpage in the top 10 HAs accounted for about 49 percent of the total ground-water pumpage. The most ground-water pumpage in an HA was due to mining in Pumpernickel Valley (HA 65), Boulder Flat (HA 61), and Lower Reese River Valley (HA 59). Pumpage by water systems in Las Vegas Valley (HA 212) and Truckee Meadows (HA 87) were the fourth and fifth highest pumpage in 2000, respectively. Irrigation and stock watering pumpage accounted for most ground-water withdrawals in the HAs with the sixth

  11. Estimating groundwater recharge

    USGS Publications Warehouse

    Stonestrom, David A.

    2011-01-01

    Groundwater recharge is the entry of fresh water into the saturated portion of the subsurface part of the hydrologic cycle, the modifier "saturated" indicating that the pressure of the pore water is greater than atmospheric.

  12. Groundwater hydrology--coastal flow

    USGS Publications Warehouse

    Sanford, Ward E.

    2010-01-01

    How groundwater flow varies when long-term external conditions change is little documented. Geochemical evidence shows that sea-level rise at the end of the last glacial period led to a shift in the flow patterns of coastal groundwater beneath Florida.

  13. Kauai Groundwater Flow Model

    DOE Data Explorer

    Nicole Lautze

    2015-01-01

    Groundwater flow model for Kauai. Data is from the following sources: Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014.; and Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume IV – Island of Kauai Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2015.

  14. ESTIMATING FLOW AND FLUX OF GROUND-WATER DISCHARGE USING WATER TEMPERATURE AND VELOCITY. (R827961)

    EPA Science Inventory

    The nature of ground water discharge to a stream has important implications for nearby ground water flow, especially with respect to contaminant transport and well-head protection. Measurements of ground water discharge were accomplished in this study using (1) differences bet...

  15. Evaluation of geohydrologic framework, recharge estimates and ground-water flow of the Joshua Tree area, San Bernardino County, California

    USGS Publications Warehouse

    Nishikawa, Tracy; Izbicki, John A.; Hevesi, Joseph A.; Stamos, Christina L.; Martin, Peter

    2005-01-01

    Ground water historically has been the sole source of water supply for the community of Joshua Tree in the Joshua Tree ground-water subbasin of the Morongo ground-water basin in the southern Mojave Desert. The Joshua Basin Water District (JBWD) supplies water to the community from the underlying Joshua Tree ground-water subbasin. The JBWD is concerned with the long-term sustainability of the underlying aquifer. To help meet future demands, the JBWD plans to construct production wells in the adjacent Copper Mountain ground-water subbasin. As growth continues in the desert, there may be a need to import water to supplement the available ground-water resources. In order to manage the ground-water resources and to identify future mitigating measures, a thorough understanding of the ground-water system is needed. The purpose of this study was threefold: (1) improve the understanding of the geohydrologic framework of the Joshua Tree and Copper Mountain ground-water subbasins, (2) determine the distribution and quantity of recharge using field and numerical techniques, and (3) develop a ground-water flow model that can be used to help manage the water resources of the region. The geohydrologic framework was refined by collecting and interpreting water-level and water-quality data, geologic and electric logs, and gravity data. The water-bearing deposits in the Joshua Tree and Copper Mountain ground-water subbasins are Quarternary alluvial deposits and Tertiary sedimentary and volcanic deposits. The Quarternary alluvial deposits were divided into two aquifers (referred to as the 'upper' and the 'middle' alluvial aquifers), which are about 600 feet (ft) thick, and the Tertiary sedimentary and volcanic deposits were assigned to a single aquifer (referred to as the 'lower' aquifer), which is as thick as 1,500 ft. The ground-water quality of the Joshua Tree and Copper Mountain ground-water subbasins was defined by collecting 53 ground-water samples from 15 wells (10 in the

  16. Using remote sensing and GIS techniques to estimate discharge and recharge. fluxes for the Death Valley regional groundwater flow system, USA

    USGS Publications Warehouse

    D'Agnese, F. A.; Faunt, C.C.; Keith, Turner A.

    1996-01-01

    The recharge and discharge components of the Death Valley regional groundwater flow system were defined by remote sensing and GIS techniques that integrated disparate data types to develop a spatially complex representation of near-surface hydrological processes. Image classification methods were applied to multispectral satellite data to produce a vegetation map. This map provided a basis for subsequent evapotranspiration and infiltration estimations. The vegetation map was combined with ancillary data in a GIS to delineate different types of wetlands, phreatophytes and wet playa areas. Existing evapotranspiration-rate estimates were then used to calculate discharge volumes for these areas. A previously used empirical method of groundwater recharge estimation was modified by GIS methods to incorporate data describing soil-moisture conditions, and a recharge potential map was produced. These discharge and recharge maps were readily converted to data arrays for numerical modelling codes. Inverse parameter estimation techniques also used these data to evaluate the reliability and sensitivity of estimated values.

  17. Death valley regional ground-water flow model calibration using optimal parameter estimation methods and geoscientific information systems

    USGS Publications Warehouse

    D'Agnese, F. A.; Faunt, C.C.; Hill, M.C.; Turner, A.K.

    1999-01-01

    A regional-scale, steady-state, saturated-zone ground-water flow model was constructed to evaluate potential regional ground-water flow in the vicinity of Yucca Mountain, Nevada. The model was limited to three layers in an effort to evaluate the characteristics governing large-scale subsurface flow. Geoscientific information systems (GSIS) were used to characterize the complex surface and subsurface hydrogeologic conditions of the area, and this characterization was used to construct likely conceptual models of the flow system. Subsurface properties in this system vary dramatically, producing high contrasts and abrupt contacts. This characteristic, combined with the large scale of the model, make zonation the logical choice for representing the hydraulic-conductivity distribution. Different conceptual models were evaluated using sensitivity analysis and were tested by using nonlinear regression to determine parameter values that are optimal, in that they provide the best match between the measured and simulated heads and flows. The different conceptual models were judged based both on the fit achieved to measured heads and spring flows, and the plausibility of the optimal parameter values. One of the conceptual models considered appears to represent the system most realistically. Any apparent model error is probably caused by the coarse vertical and horizontal discretization.A regional-scale, steady-state, saturated-zone ground-water flow model was constructed to evaluate potential regional ground-water flow in the vicinity of Yucca Mountain, Nevada. The model was limited to three layers in an effort to evaluate the characteristics governing large-scale subsurface flow. Geoscientific information systems (GSIS) were used to characterize the complex surface and subsurface hydrogeologic conditions of the area, and this characterization was used to construct likely conceptual models of the flow system. Subsurface properties in this system vary dramatically, producing

  18. A hybrid framework for improving recharge and discharge estimation for a three-dimensional groundwater flow model.

    PubMed

    Meyer, Scott C; Lin, Yu-Feng; Roadcap, George S

    2012-01-01

    We employed the ArcGIS plug-in package PRO-GRADE (Lin et al. 2009), developed for zonation of recharge/discharge (R/D) for modeling two-dimensional aquifer systems, to develop alternative R/D zonations for an existing three-dimensional groundwater flow model of a complex hydrogeologic setting. Our process began by intersecting PRO-GRADE output with the existing model's 4-zone R/D representation to develop a model having 12 R/D zones (R12) and then calibrating the resulting model using PEST. We then revised the R12 zonation using supplementary GIS data to develop a 51-zone R/D zonation (R51). From R51, we developed a series of daughter models having 40, 30, 28, and 18 R/D zones by removing zones from R51 if calibration resulted in little change in the zone's starting R/D rate and/or if the model was insensitive to the zone's R/D rate. For these models (R40N, R30N, R28N, and R18N), we used the ArcGIS Nibble tool to rapidly and consistently reassign model cells within eliminated zones of R51 to the zone of the nearest model cell in a retained zone having the same starting value. R12, R51, R40N, R30N, R28N, and R18N are all more accurate than the original model (R4), although improvements relative to stream discharge targets exceeded improvements relative to head targets. The models also executed with better numerical stability and less mass balance discrepancy than R4. These improvements demonstrate that R/D estimation in a complex shallow three-dimensional steady-state model can be improved with PRO-GRADE estimates of R/D when guided by calibration statistics and supplemental geographic data. PMID:21797853

  19. Using remote sensing and GIS techniques to estimate discharge and recharge fluxes for the Death Valley regional groundwater flow system, USA

    USGS Publications Warehouse

    D'Agnese, F. A.; Faunt, C.C.; Turner, A.K.

    1996-01-01

    The recharge and discharge components of the Death Valley regional groundwater flow system were defined by techniques that integrated disparate data types to develop a spatially complex representation of near-surface hydrological processes. Image classification methods were applied to multispectral satellite data to produce a vegetation map. The vegetation map was combined with ancillary data in a GIS to delineate different types of wetlands, phreatophytes and wet playa areas. Existing evapotranspiration-rate estimates were used to calculate discharge volumes for these area. An empirical method of groundwater recharge estimation was modified to incorporate data describing soil-moisture conditions, and a recharge potential map was produced. These discharge and recharge maps were readily converted to data arrays for numerical modelling codes. Inverse parameter estimation techniques also used these data to evaluate the reliability and sensitivity of estimated values.The recharge and discharge components of the Death Valley regional groundwater flow system were defined by remote sensing and GIS techniques that integrated disparate data types to develop a spatially complex representation of near-surface hydrological processes. Image classification methods were applied to multispectral satellite data to produce a vegetation map. This map provided a basis for subsequent evapotranspiration and infiltration estimations. The vegetation map was combined with ancillary data in a GIS to delineate different types of wetlands, phreatophytes and wet playa areas. Existing evapotranspiration-rate estimates were then used to calculate discharge volumes for these areas. A previously used empirical method of groundwater recharge estimation was modified by GIS methods to incorporate data describing soil-moisture conditions, and a recharge potential map was produced. These discharge and recharge maps were readily converted to data arrays for numerical modelling codes. Inverse parameter

  20. Monitoring probe for groundwater flow

    DOEpatents

    Looney, B.B.; Ballard, S.

    1994-08-23

    A monitoring probe for detecting groundwater migration is disclosed. The monitor features a cylinder made of a permeable membrane carrying an array of electrical conductivity sensors on its outer surface. The cylinder is filled with a fluid that has a conductivity different than the groundwater. The probe is placed in the ground at an area of interest to be monitored. The fluid, typically saltwater, diffuses through the permeable membrane into the groundwater. The flow of groundwater passing around the permeable membrane walls of the cylinder carries the conductive fluid in the same general direction and distorts the conductivity field measured by the sensors. The degree of distortion from top to bottom and around the probe is precisely related to the vertical and horizontal flow rates, respectively. The electrical conductivities measured by the sensors about the outer surface of the probe are analyzed to determine the rate and direction of the groundwater flow. 4 figs.

  1. Monitoring probe for groundwater flow

    DOEpatents

    Looney, Brian B.; Ballard, Sanford

    1994-01-01

    A monitoring probe for detecting groundwater migration. The monitor features a cylinder made of a permeable membrane carrying an array of electrical conductivity sensors on its outer surface. The cylinder is filled with a fluid that has a conductivity different than the groundwater. The probe is placed in the ground at an area of interest to be monitored. The fluid, typically saltwater, diffuses through the permeable membrane into the groundwater. The flow of groundwater passing around the permeable membrane walls of the cylinder carries the conductive fluid in the same general direction and distorts the conductivity field measured by the sensors. The degree of distortion from top to bottom and around the probe is precisely related to the vertical and horizontal flow rates, respectively. The electrical conductivities measured by the sensors about the outer surface of the probe are analyzed to determine the rate and direction of the groundwater flow.

  2. Estimates of tracer-based piston-flow ages of groundwater from selected sites-National Water-Quality Assessment Program, 1992-2005

    USGS Publications Warehouse

    Hinkle, Stephen R.; Shapiro, Stephanie D.; Plummer, L. Niel; Busenberg, Eurybiades; Widman, Peggy K.; Casile, Gerolamo C.; Wayland, Julian E.

    2011-01-01

    This report documents selected age data interpreted from measured concentrations of environmental tracers in groundwater from 1,399 National Water-Quality Assessment (NAWQA) Program groundwater sites across the United States. The tracers of interest were chlorofluorocarbons (CFCs), sulfur hexafluoride (SF6), and tritium/helium-3 (3H/3He). Tracer data compiled for this analysis primarily were from wells representing two types of NAWQA groundwater studies - Land-Use Studies (shallow wells, usually monitoring wells, in recharge areas under dominant land-use settings) and Major-Aquifer Studies (wells, usually domestic supply wells, in principal aquifers and representing the shallow, used resource). Reference wells (wells representing groundwater minimally impacted by anthropogenic activities) associated with Land-Use Studies also were included. Tracer samples were collected between 1992 and 2005, although two networks sampled from 2006 to 2007 were included because of network-specific needs. Tracer data from other NAWQA Program components (Flow System Studies, which are assessments of processes and trends along groundwater flow paths, and various topical studies) were not compiled herein. Tracer data from NAWQA Land-Use Studies and Major-Aquifer Studies that previously had been interpreted and published are compiled herein (as piston-flow ages), but have not been reinterpreted. Tracer data that previously had not been interpreted and published are evaluated using documented methods and compiled with aqueous concentrations, equivalent atmospheric concentrations (for CFCs and SF6), estimates of tracer-based piston-flow ages, and selected ancillary data, such as redox indicators, well construction, and major dissolved gases (N2, O2, Ar, CH4, and CO2). Tracer-based piston-flow ages documented in this report are simplistic representations of the tracer data. Tracer-based piston-flow ages are a convenient means of conceptualizing groundwater age. However, the piston-flow

  3. Natural recharge estimation and uncertainty analysis of an adjudicated groundwater basin using a regional-scale flow and subsidence model (Antelope Valley, California, USA)

    USGS Publications Warehouse

    Siade, Adam J.; Nishikawa, Tracy; Martin, Peter

    2015-01-01

    Groundwater has provided 50–90 % of the total water supply in Antelope Valley, California (USA). The associated groundwater-level declines have led the Los Angeles County Superior Court of California to recently rule that the Antelope Valley groundwater basin is in overdraft, i.e., annual pumpage exceeds annual recharge. Natural recharge consists primarily of mountain-front recharge and is an important component of the total groundwater budget in Antelope Valley. Therefore, natural recharge plays a major role in the Court’s decision. The exact quantity and distribution of natural recharge is uncertain, with total estimates from previous studies ranging from 37 to 200 gigaliters per year (GL/year). In order to better understand the uncertainty associated with natural recharge and to provide a tool for groundwater management, a numerical model of groundwater flow and land subsidence was developed. The transient model was calibrated using PEST with water-level and subsidence data; prior information was incorporated through the use of Tikhonov regularization. The calibrated estimate of natural recharge was 36 GL/year, which is appreciably less than the value used by the court (74 GL/year). The effect of parameter uncertainty on the estimation of natural recharge was addressed using the Null-Space Monte Carlo method. A Pareto trade-off method was also used to portray the reasonableness of larger natural recharge rates. The reasonableness of the 74 GL/year value and the effect of uncertain pumpage rates were also evaluated. The uncertainty analyses indicate that the total natural recharge likely ranges between 34.5 and 54.3 GL/year.

  4. Natural recharge estimation and uncertainty analysis of an adjudicated groundwater basin using a regional-scale flow and subsidence model (Antelope Valley, California, USA)

    NASA Astrophysics Data System (ADS)

    Siade, Adam; Nishikawa, Tracy; Martin, Peter

    2015-09-01

    Groundwater has provided 50-90 % of the total water supply in Antelope Valley, California (USA). The associated groundwater-level declines have led the Los Angeles County Superior Court of California to recently rule that the Antelope Valley groundwater basin is in overdraft, i.e., annual pumpage exceeds annual recharge. Natural recharge consists primarily of mountain-front recharge and is an important component of the total groundwater budget in Antelope Valley. Therefore, natural recharge plays a major role in the Court's decision. The exact quantity and distribution of natural recharge is uncertain, with total estimates from previous studies ranging from 37 to 200 gigaliters per year (GL/year). In order to better understand the uncertainty associated with natural recharge and to provide a tool for groundwater management, a numerical model of groundwater flow and land subsidence was developed. The transient model was calibrated using PEST with water-level and subsidence data; prior information was incorporated through the use of Tikhonov regularization. The calibrated estimate of natural recharge was 36 GL/year, which is appreciably less than the value used by the court (74 GL/year). The effect of parameter uncertainty on the estimation of natural recharge was addressed using the Null-Space Monte Carlo method. A Pareto trade-off method was also used to portray the reasonableness of larger natural recharge rates. The reasonableness of the 74 GL/year value and the effect of uncertain pumpage rates were also evaluated. The uncertainty analyses indicate that the total natural recharge likely ranges between 34.5 and 54.3 GL/year.

  5. A tidal creek water budget: Estimation of groundwater discharge and overland flow using hydrologic modeling in the Southern Everglades

    NASA Astrophysics Data System (ADS)

    Michot, Béatrice; Meselhe, Ehab A.; Rivera-Monroy, Victor H.; Coronado-Molina, Carlos; Twilley, Robert R.

    2011-07-01

    Taylor Slough is one of the natural freshwater contributors to Florida Bay through a network of microtidal creeks crossing the Everglades Mangrove Ecotone Region (EMER). The EMER ecological function is critical since it mediates freshwater and nutrient inputs and controls the water quality in Eastern Florida Bay. Furthermore, this region is vulnerable to changing hydrodynamics and nutrient loadings as a result of upstream freshwater management practices proposed by the Comprehensive Everglades Restoration Program (CERP), currently the largest wetland restoration project in the USA. Despite the hydrological importance of Taylor Slough in the water budget of Florida Bay, there are no fine scale (˜1 km 2) hydrodynamic models of this system that can be utilized as a tool to evaluate potential changes in water flow, salinity, and water quality. Taylor River is one of the major creeks draining Taylor Slough freshwater into Florida Bay. We performed a water budget analysis for the Taylor River area, based on long-term hydrologic data (1999-2007) and supplemented by hydrodynamic modeling using a MIKE FLOOD (DHI, http://dhigroup.com/) model to evaluate groundwater and overland water discharges. The seasonal hydrologic characteristics are very distinctive (average Taylor River wet vs. dry season outflow was 6 to 1 during 1999-2006) with a pronounced interannual variability of flow. The water budget shows a net dominance of through flow in the tidal mixing zone, while local precipitation and evapotranspiration play only a secondary role, at least in the wet season. During the dry season, the tidal flood reaches the upstream boundary of the study area during approximately 80 days per year on average. The groundwater field measurements indicate a mostly upwards-oriented leakage, which possibly equals the evapotranspiration term. The model results suggest a high importance of groundwater contribution to the water salinity in the EMER. The model performance is satisfactory

  6. Effects of wetlands creation on groundwater flow

    USGS Publications Warehouse

    Hensel, B.R.; Miller, M.V.

    1991-01-01

    Changes in groundwater flow were observed near four Experimental Wetland Areas (EWAs) constructed along a reach of the Des Plaines River in northeastern Illinois. These changes were observed during monthly monitoring of groundwater elevation in nested piezometers and shallow observation wells before and after the wetlands were filled with water. A numerical model was calibrated with observed data and used to estimate seepage from the wetlands into the Des Plaines River. After the wetlands became operational, groundwater levels in adjacent wells increased by about 0.5m, while water levels in wells distant from the wetlands decreased. The increase in groundwater levels near the wetlands is a result of seepage from the wetlands. Numerical predictions of seepage from the wetlands are 60-150 m3 day-1 for two wetlands situated over sand and gravel and less than 1 m3 day-1 for two wetlands situated over clayey till. The difference in seepage rates is attributed to two factors. First, the hydraulic conductivity of the sand and gravel unit is greater than that of the till, and thus there is less mounding and a greater capacity for transmitting water beneath the wetlands overlying this deposit. Secondly, the wetlands located over till are groundwater flow-through ponds, whereas the wetlands over the sand and gravel are primarily groundwater recharge areas. The model was used to estimate that seepage from the wetlands will double groundwater discharge into the Des Plaines River and a tributary relative to pre-operational discharge from the study area. Overall, the wetlands have acted as a constant head boundary, stabilizing groundwater flow patterns. ?? 1991.

  7. Geothermal properties and groundwater flow estimated with a three-dimensional geological model in a late Pleistocene terrace area, central Japan

    NASA Astrophysics Data System (ADS)

    Funabiki, A.; Takemura, T.; Hamamoto, S.; Komatsu, T.

    2012-12-01

    1. Introduction The ground source heat pump (GSHP) is a highly efficient and renewable energy technology for space heating and cooling, with benefits that include energy conservation and reductions in greenhouse gas emissions. One result of the huge Tohoku-oki earthquake and tsunami and the subsequent nuclear disasters is that GSHPs are receiving more attention from the media and they are being introduced by some local governments. Heat generated by underground GSHP installation, however, can pollute the geothermal environment or change groundwater flow patterns . In this study, we estimated possible effects from the use of GSHPs in the Tokyo area with a three-dimensional (3D) geological model. 2. Geological model The Tokyo Metropolitan Area is surrounded by the Late Pleistocene terraces called the Musashino uplands. The terrace surfaces are densely populated residential areas. One of these surfaces, the Shimosueyohi surface, formed along the Tama River during the last deglacial period. The CRE-NUCHS-1 core (Funabiki et al., 2011) was obtained from this surface, and the lithology, heat transfer coefficients, and chemical characteristics of the sediments were analyzed. In this study, we used borehole log data from a 5 km2 area surrounding the CRE-NUCHS-1 core site to create a 3D geological model. In this area, the Pleistocene Kazusa Group is overlain by terrace gravels and a volcanic ash layer called the Kanto Loam. The terrace gravels occur mainly beneath the Kanda, Kitazawa, and Karasuyama rivers , which flow parallel to the Tama River, whereas away from the rivers , the Kanto Loam directly overlies the Kazusa Group sediments. 3. Geothermal disturbance and groundwater flow Using the geological model, we calculated the heat transfer coefficients and groundwater flow velocities in the sediments. Within the thick terrace gravels, which are at relatively shallow depth (8-20 m), heat transfer coefficients were high and groundwater flow was relatively fast. The amount

  8. Modeling groundwater flow on MPPs

    SciTech Connect

    Ashby, S.F.; Falgout, R.D.; Smith, S.G.; Tompson, A.F.B.

    1993-10-01

    The numerical simulation of groundwater flow in three-dimensional heterogeneous porous media is examined. To enable detailed modeling of large contaminated sites, preconditioned iterative methods and massively parallel computing power are combined in a simulator called PARFLOW. After describing this portable and modular code, some numerical results are given, including one that demonstrates the code`s scalability.

  9. Real-time groundwater flow modeling with the Ensemble Kalman Filter: Joint estimation of states and parameters and the filter inbreeding problem

    NASA Astrophysics Data System (ADS)

    Hendricks Franssen, H. J.; Kinzelbach, W.

    2008-09-01

    Real-time groundwater flow modeling with filter methods is interesting for dynamical groundwater flow systems, for which measurement data in real-time are available. The Ensemble Kalman Filter (EnKF) approach is used here to update states together with parameters by adopting an augmented state vector approach. The performance of EnKF is investigated in a synthetic study with a two-dimensional transient groundwater flow model where (1) only the recharge rate is spatiotemporally variable, (2) only transmissivity is spatially variable with σlnT2 = 1.0 or (3) with σlnT2 = 2.7, and (4) both recharge rate and transmissivity are uncertain (a combination of (1) and (3)). The performance of EnKF for simultaneous state and parameter estimation in saturated groundwater flow problems is investigated in dependence of the number of stochastic realizations, the updating frequency and updating intensity of log-transmissivity, the amount of measurements in space and time, and the method (iterative versus noniterative EnKF), among others. Satisfactory results were also obtained if both transmissivity and recharge rate were uncertain. However, it was found that filter inbreeding is much more severe if hydraulic heads and transmissivities are jointly updated than if only hydraulic heads are updated. The filter inbreeding problem was investigated in more detail and could be strongly reduced with help of a damping parameter, which limits the intensity of the perturbation of the log-transmissivity field. An additional reduction of filter inbreeding could be achieved by combining two measures: (1) inflating the elements of the predicted state covariance matrix on the basis of a comparison between the model uncertainty and the observed errors at the measurement points and (2) starting the flow simulations with a very large number of realizations and then sampling the desired number of realizations after one simulation time step by minimizing the differences between the local cpdfs (and

  10. Comparison of methods for estimating ground-water recharge and base flow at a small watershed underlain by fractured bedrock in the Eastern United States

    USGS Publications Warehouse

    Risser, Dennis W.; Gburek, William J.; Folmar, Gordon J.

    2005-01-01

    This study by the U.S. Geological Survey (USGS), in cooperation with the Agricultural Research Service (ARS), U.S. Department of Agriculture, compared multiple methods for estimating ground-water recharge and base flow (as a proxy for recharge) at sites in east-central Pennsylvania underlain by fractured bedrock and representative of a humid-continental climate. This study was one of several within the USGS Ground-Water Resources Program designed to provide an improved understanding of methods for estimating recharge in the eastern United States. Recharge was estimated on a monthly and annual basis using four methods?(1) unsaturated-zone drainage collected in gravity lysimeters, (2) daily water balance, (3) water-table fluctuations in wells, and (4) equations of Rorabaugh. Base flow was estimated by streamflow-hydrograph separation using the computer programs PART and HYSEP. Estimates of recharge and base flow were compared for an 8-year period (1994-2001) coinciding with operation of the gravity lysimeters at an experimental recharge site (Masser Recharge Site) and a longer 34-year period (1968-2001), for which climate and streamflow data were available on a 2.8-square-mile watershed (WE-38 watershed). Estimates of mean-annual recharge at the Masser Recharge Site and WE-38 watershed for 1994-2001 ranged from 9.9 to 14.0 inches (24 to 33 percent of precipitation). Recharge, in inches, from the various methods was: unsaturated-zone drainage, 12.2; daily water balance, 12.3; Rorabaugh equations with PULSE, 10.2, or RORA, 14.0; and water-table fluctuations, 9.9. Mean-annual base flow from streamflow-hydrograph separation ranged from 9.0 to 11.6 inches (21-28 percent of precipitation). Base flow, in inches, from the various methods was: PART, 10.7; HYSEP Local Minimum, 9.0; HYSEP Sliding Interval, 11.5; and HYSEP Fixed Interval, 11.6. Estimating recharge from multiple methods is useful, but the inherent differences of the methods must be considered when comparing

  11. Estimates of consumptive use and ground-water return flow and the effect of rising and sustained high river stage on the method of estimation in Cibola Valley, Arizona and California, 1983 and 1984

    USGS Publications Warehouse

    Owen-Joyce, Sandra J.

    1990-01-01

    In Cibola Valley, Arizona, water is pumped from the Colorado River to irrigate crops and to maintain wildlife habitat. Unused water percolates to the water table and, as groundwater, moves downgradient into areas of phreatophytes, into a drainage ditch, out of the flood plain, and back to the river. In 1983 and 1984, groundwater return flow was negligible because in most of Cibola Valley the river lost water to the aquifer. Evapotranspiration was used as an approximation for consumptive use by vegetation. Evapotranspiration was calculated as the sum of the products of the area of vegetation types and water-use rate by vegetation type. Evapotranspiration was estimated to be 70,100 acre-ft in 1983 and 62,600 acre-ft in 1984. These estimates may be in error because of the effect of sustained inundation on the rate of water use by phreatophytes. The effects cannot be quantified and therefore adjustments to rates calculated for dry-surface conditions could not be made. The method of estimating consumptive use of water by vegetation and groundwater return flow is affected by changing conditions during years of rising and sustained high river stage caused by flood-control releases at Parker Dam. Most of the bank storage that will return to the river when the high river stage subsides did not originate as irrigation water. High river stage caused some areas to be flooded directly or raised groundwater levels above the land surface. No crops could be grown in flooded fields. The decreased depth to water and inundation with fresh water resulted in new phreatophyte growth in some areas. In some areas that were flooded, many phreatophytes died. Changes in the inundated and flooded areas throughout the years made it difficult to estimate the evaporation losses from the increased water surface. (USGS)

  12. Estimation of groundwater flow directions and the tensor of hydraulic conductivity in crystalline massif rocks using information from surface structural geology and mining exploration boreholes

    NASA Astrophysics Data System (ADS)

    Florez, C.; Romero, M. A.; Ramirez, M. I.; Monsalve, G.

    2013-05-01

    In the elaboration of a hydrogeological conceptual model in regions of mining exploration where there is significant presence of crystalline massif rocks., the influence of physical and geometrical properties of rock discontinuities must be evaluated. We present the results of a structural analysis of rock discontinuities in a region of the Central Cordillera of Colombia (The upper and middle Bermellon Basin) in order to establish its hydrogeological characteristics for the improvement of the conceptual hydrogeological model for the region. The geology of the study area consists of schists with quartz and mica and porphyritic rocks, in a region of high slopes with a nearly 10 m thick weathered layer. The main objective of this research is to infer the preferential flow directions of groundwater and to estimate the tensor of potential hydraulic conductivity by using surface information and avoiding the use of wells and packer tests. The first step of our methodology is an analysis of drainage directions to detect patterns of structural controls in the run-off; after a field campaign of structural data recollection, where we compile information of strike, dip, continuity, spacing, roughness, aperture and frequency, we built equal area hydro-structural polar diagrams that indicate the potential directions for groundwater flow. These results are confronted with records of Rock Quality Designation (RQD) that have been systematically taken from several mining exploration boreholes in the area of study. By using all this information we estimate the potential tensor of hydraulic conductivity from a cubic law, obtaining the three principal directions with conductivities of the order of 10-5 and 10-6 m/s; the more conductive joint family has a NE strike with a nearly vertical dip.

  13. Estimating Age Distributions of Base Flow in Watersheds Underlain by Single and Dual Porosity Formations Using Groundwater Transport Simulation and Weighted Weibull Functions

    NASA Astrophysics Data System (ADS)

    Sanford, W. E.

    2015-12-01

    Age distributions of base flow to streams are important to estimate for predicting the timing of water-quality responses to changes in distributed inputs of nutrients or pollutants at the land surface. Simple models of shallow aquifers will predict exponential age distributions, but more realistic 3-D stream-aquifer geometries will cause deviations from an exponential curve. In addition, in fractured rock terrains the dual nature of the effective and total porosity of the system complicates the age distribution further. In this study shallow groundwater flow and advective transport were simulated in two regions in the Eastern United States—the Delmarva Peninsula and the upper Potomac River basin. The former is underlain by layers of unconsolidated sediment, while the latter consists of folded and fractured sedimentary rocks. Transport of groundwater to streams was simulated using the USGS code MODPATH within 175 and 275 watersheds, respectively. For the fractured rock terrain, calculations were also performed along flow pathlines to account for exchange between mobile and immobile flow zones. Porosities at both sites were calibrated using environmental tracer data (3H, 3He, CFCs and SF6) in wells and springs, and with a 30-year tritium record from the Potomac River. Carbonate and siliciclastic rocks were calibrated to have mobile porosity values of one and six percent, and immobile porosity values of 18 and 12 percent, respectively. The age distributions were fitted to Weibull functions. Whereas an exponential function has one parameter that controls the median age of the distribution, a Weibull function has an extra parameter that controls the slope of the curve. A weighted Weibull function was also developed that potentially allows for four parameters, two that control the median age and two that control the slope, one of each weighted toward early or late arrival times. For both systems the two-parameter Weibull function nearly always produced a substantially

  14. West Maui Groundwater Flow Model

    DOE Data Explorer

    Nicole Lautze

    2015-01-01

    Groundwater flow model for West Maui. Data is from the following sources: Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014; and Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume V – Island of Maui Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008.

  15. Hawaii Island Groundwater Flow Model

    DOE Data Explorer

    Nicole Lautze

    2015-01-01

    Groundwater flow model for Hawaii Island. Data is from the following sources: Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume II – Island of Hawaii Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008; and Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014.

  16. East Maui Groundwater Flow Model

    DOE Data Explorer

    Nicole Lautze

    2015-01-01

    Groundwater flow model for East Maui. Data is from the following sources: Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014; and Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume V – Island of Maui Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008.

  17. Death Valley regional groundwater flow system, Nevada and California-Hydrogeologic framework and transient groundwater flow model

    USGS Publications Warehouse

    : Belcher, Wayne R., (Edited By); Sweetkind, Donald S.

    2010-01-01

    by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional groundwater flow system. Groundwater flow in the Death Valley region is composed of several interconnected, complex groundwater flow systems. Groundwater flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional groundwater flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and ET are the dominant natural groundwater discharge processes. Groundwater also is withdrawn for agricultural, commercial, and domestic uses. Groundwater flow in the DVRFS was simulated using MODFLOW-2000, the U.S. Geological Survey 3D finitedifference modular groundwater flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 meters (m) on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins simulated by constant-head boundaries. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient mode

  18. Assimilating ambiguous observations to jointly estimate groundwater recharge and conductivity

    NASA Astrophysics Data System (ADS)

    Erdal, Daniel; Cirpka, Olaf A.

    2016-04-01

    In coupled modelling of catchments, the groundwater compartment can be an important water storage as well as having influence on both rivers and evapotranspirational fluxes. It is therefore important to parameterize the groundwater model as correctly as possible. Primarily important to regional groundwater flow is the spatially variable hydraulic conductivity. However, also the groundwater recharge, in a coupled system coming from the unsaturated zone but in a stand-alone groundwater model a boundary condition, is also of high importance. As with all subsurface systems, groundwater properties are difficult to observe in reality and their estimation is an ongoing topic in groundwater research and practice. Commonly, we have to rely on time series of groundwater head observations as base for any parameter estimation. Heads, however, have the drawback that they can be ambiguous and may not uniquely define the inverse problem, especially if both recharge and conductivity are seen as unknown. In the presented work we use a 2D virtual groundwater test case to investigate how the prior knowledge of recharge and conductivity influence their respective and joint estimation as spatially variable fields using head data. Using the Ensemble Kalman filter, it is shown that the joint estimation is possible if the prior knowledge is good enough. If the prior is erroneous the a-priori sampled fields cannot be corrected by the data. However, it is also shown that if the prior knowledge is directly wrong the estimated recharge field can resemble the true conductivity field, resulting in a model that meets the observations but has very poor predictive power. The study exemplifies the importance of prior knowledge in the joint estimation of parameters from ambiguous measurements.

  19. Geomorphic aspects of groundwater flow

    NASA Astrophysics Data System (ADS)

    LaFleur, Robert G.

    The many roles that groundwater plays in landscape evolution are becoming more widely appreciated. In this overview, three major categories of groundwater processes and resulting landforms are considered: (1) Dissolution creates various karst geometries, mainly in carbonate rocks, in response to conditions of recharge, geologic setting, lithology, and groundwater circulation. Denudation and cave formation rates can be estimated from kinetic and hydraulic parameters. (2) Groundwater weathering generates regoliths of residual alteration products at weathering fronts, and subsequent exhumation exposes corestones, flared slopes, balanced rocks, domed inselbergs, and etchplains of regional importance. Groundwater relocation of dissolved salts creates duricrusts of various compositions, which become landforms. (3) Soil and rock erosion by groundwater processes include piping, seepage erosion, and sapping, important agents in slope retreat and headward gully migration. Thresholds and limits are important in many chemical and mechanical groundwater actions. A quantitative, morphometric approach to groundwater landforms and processes is exemplified by selected studies in carbonate and clastic terrains of ancient and recent origins. Résumé Les rôles variés joués par les eaux souterraines dans l'évolution des paysages deviennent nettement mieux connus. La revue faite ici prend en considération trois grandes catégories de processus liés aux eaux souterraines et les formes associées: (1) La dissolution crée des formes karstiques variées, surtout dans les roches carbonatées, en fonction des conditions d'alimentation, du cadre géologique, de la lithologie et de la circulation des eaux souterraines. Les taux d'érosion et de formation des grottes peuvent être estimés à partir de paramètres cinétiques et hydrauliques. (2) L'érosion par les eaux souterraines donne naissance à des régolites, résidus d'altération sur des fronts d'altération, et l'exhumation r

  20. Geomorphic aspects of groundwater flow

    NASA Astrophysics Data System (ADS)

    LaFleur, Robert G.

    The many roles that groundwater plays in landscape evolution are becoming more widely appreciated. In this overview, three major categories of groundwater processes and resulting landforms are considered: (1) Dissolution creates various karst geometries, mainly in carbonate rocks, in response to conditions of recharge, geologic setting, lithology, and groundwater circulation. Denudation and cave formation rates can be estimated from kinetic and hydraulic parameters. (2) Groundwater weathering generates regoliths of residual alteration products at weathering fronts, and subsequent exhumation exposes corestones, flared slopes, balanced rocks, domed inselbergs, and etchplains of regional importance. Groundwater relocation of dissolved salts creates duricrusts of various compositions, which become landforms. (3) Soil and rock erosion by groundwater processes include piping, seepage erosion, and sapping, important agents in slope retreat and headward gully migration. Thresholds and limits are important in many chemical and mechanical groundwater actions. A quantitative, morphometric approach to groundwater landforms and processes is exemplified by selected studies in carbonate and clastic terrains of ancient and recent origins. Résumé Les rôles variés joués par les eaux souterraines dans l'évolution des paysages deviennent nettement mieux connus. La revue faite ici prend en considération trois grandes catégories de processus liés aux eaux souterraines et les formes associées: (1) La dissolution crée des formes karstiques variées, surtout dans les roches carbonatées, en fonction des conditions d'alimentation, du cadre géologique, de la lithologie et de la circulation des eaux souterraines. Les taux d'érosion et de formation des grottes peuvent être estimés à partir de paramètres cinétiques et hydrauliques. (2) L'érosion par les eaux souterraines donne naissance à des régolites, résidus d'altération sur des fronts d'altération, et l'exhumation r

  1. Death Valley regional ground-water flow system, Nevada and California -- hydrogeologic framework and transient ground-water flow model

    USGS Publications Warehouse

    : Belcher, Wayne R., (Edited By)

    2004-01-01

    provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional ground-water flow system. Ground-water flow in the Death Valley region is composed of several interconnected, complex ground-water flow systems. Ground-water flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional ground-water flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and evapotranspiration (ET) are the dominant natural ground-water discharge processes. Ground water also is withdrawn for agricultural, commercial, and domestic uses. Ground-water flow in the DVRFS was simulated using MODFLOW-2000, a 3D finite-difference modular ground-water flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 m on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins, also simulated by drains. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model gener

  2. Using Groundwater Temperatures and Heat Flow Patterns to Assess Groundwater Flow in Snake Valley, Nevada and Utah, USA

    NASA Astrophysics Data System (ADS)

    Masbruch, M. D.; Chapman, D. S.

    2009-12-01

    The Southern Nevada Water Authority’s (SNWA) proposal to develop groundwater resources in Snake Valley and adjacent basins in eastern Nevada has focused attention on understanding the links between basin-fill and carbonate aquifer systems, groundwater flow paths, and the movement of groundwater between basins. The SNWA development plans are contentious in part because (1) there are few perennial streams that flow into the basins and these surface-water resources are fully appropriated; (2) groundwater resources that sustain streams, springs, wetlands, and the local agricultural economy are also limited; and (3) because Snake Valley straddles the Utah-Nevada state line. We report groundwater temperatures and estimates of heat flow used to constrain estimates of groundwater flow into and through Snake Valley. Thermal logs have been collected from 24 monitoring wells in the Utah part of the valley. Natural, undisturbed geothermal gradients within the Basin and Range are generally 30 °C/km, which correspond to heat flow values of approximately 90 mW/m2. Geothermal gradients in the southern portion of Snake Valley are lower than typical Basin and Range geothermal gradients, with the majority ranging between 10 and 20 °C/km, corresponding to heat flow values of 30 to 60 mW/m2. In the northern portion of the basin, however, geothermal gradients are generally higher than typical Basin and Range geothermal gradients, with thermal logs of two wells indicating gradients of 39 °C/km and 51 °C/km, which correspond to heat flow values of approximately 117 and 153 mW/m2, respectively. These observations suggest heat is being redistributed by groundwater flow to discharge points in northern Snake Valley. This interpretation is also supported by spring temperatures in northern Snake Valley and at Fish Springs National Wildlife Refuge to the northeast that are higher than ambient (12 °C) surface temperature. These thermal data are being used together with water levels and

  3. Using groundwater levels to estimate recharge

    USGS Publications Warehouse

    Healy, R.W.; Cook, P.G.

    2002-01-01

    Accurate estimation of groundwater recharge is extremely important for proper management of groundwater systems. Many different approaches exist for estimating recharge. This paper presents a review of methods that are based on groundwater-level data. The water-table fluctuation method may be the most widely used technique for estimating recharge; it requires knowledge of specific yield and changes in water levels over time. Advantages of this approach include its simplicity and an insensitivity to the mechanism by which water moves through the unsaturated zone. Uncertainty in estimates generated by this method relate to the limited accuracy with which specific yield can be determined and to the extent to which assumptions inherent in the method are valid. Other methods that use water levels (mostly based on the Darcy equation) are also described. The theory underlying the methods is explained. Examples from the literature are used to illustrate applications of the different methods.

  4. Influence of perched groundwater on base flow

    USGS Publications Warehouse

    Niswonger, R.G.; Fogg, G.E.

    2008-01-01

    Analysis with a three-dimensional variably saturated groundwater flow model provides a basic understanding of the interplay between streams and perched groundwater. A simplified, layered model of heterogeneity was used to explore these relationships. Base flow contribution from perched groundwater was evaluated with regard to varying hydrogeologic conditions, including the size and location of the fine-sediment unit and the hydraulic conductivity of the fine-sediment unit and surrounding coarser sediment. Simulated base flow was sustained by perched groundwater with a maximum monthly discharge in excess of 15 L/s (0.6 feet3/s) over the length of the 2000-m stream reach. Generally, the rate of perched-groundwater discharge to the stream was proportional to the hydraulic conductivity of sediment surrounding the stream, whereas the duration of discharge was proportional to the hydraulic conductivity of the fine-sediment unit. Other aspects of the perched aquifer affected base flow, such as the depth of stream penetration and the size of the fine-sediment unit. Greater stream penetration decreased the maximum base flow contribution but increased the duration of contribution. Perched groundwater provided water for riparian vegetation at the demand rate but reduced the duration of perched-groundwater discharge nearly 75%. Copyright 2008 by the American Geophysical Union.

  5. Revised conceptualization of the North China Basin groundwater flow system: Groundwater age, heat and flow simulations

    NASA Astrophysics Data System (ADS)

    Cao, Guoliang; Han, Dongmei; Currell, Matthew J.; Zheng, Chunmiao

    2016-09-01

    Groundwater flow in deep sedimentary basins results from complex evolution processes on geological timescales. Groundwater flow systems conceptualized according to topography and/or groundwater table configuration generally assume a near-equilibrium state with the modern landscape. However, the time to reach such a steady state, and more generally the timescales of groundwater flow system evolution are key considerations for large sedimentary basins. This is true in the North China Basin (NCB), which has been studied for many years due to its importance as a groundwater supply. Despite many years of study, there remain contradictions between the generally accepted conceptual model of regional flow, and environmental tracer data. We seek to reconcile these contractions by conducting simulations of groundwater flow, age and heat transport in a three dimensional model, using an alternative conceptual model, based on geological, thermal, isotope and historical data. We infer flow patterns under modern hydraulic conditions using this new model and present the theoretical maximum groundwater ages under such a flow regime. The model results show that in contrast to previously accepted conceptualizations, most groundwater is discharged in the vicinity of the break-in-slope of topography at the boundary between the piedmont and central plain. Groundwater discharge to the ocean is in contrast small, and in general there are low rates of active flow in the eastern parts of the basin below the central and coastal plain. This conceptualization is more compatible with geochemical and geothermal data than the previous model. Simulated maximum groundwater ages of ∼1 Myrs below the central and coastal plain indicate that residual groundwater may be retained in the deep parts of the basin since being recharged during the last glacial period or earlier. The groundwater flow system has therefore probably not reached a new equilibrium state with modern-day hydraulic conditions. The

  6. Traffic Flow Estimates.

    ERIC Educational Resources Information Center

    Hart, Vincent G.

    1981-01-01

    Two examples are given of ways traffic engineers estimate traffic flow. The first, Floating Car Method, involves some basic ideas and the notion of relative velocity. The second, Maximum Traffic Flow, is viewed to involve simple applications of calculus. The material provides insight into specialized applications of mathematics. (MP)

  7. Groundwater recharge rate and zone structure estimation using PSOLVER algorithm.

    PubMed

    Ayvaz, M Tamer; Elçi, Alper

    2014-01-01

    The quantification of groundwater recharge is an important but challenging task in groundwater flow modeling because recharge varies spatially and temporally. The goal of this study is to present an innovative methodology to estimate groundwater recharge rates and zone structures for regional groundwater flow models. Here, the unknown recharge field is partitioned into a number of zones using Voronoi Tessellation (VT). The identified zone structure with the recharge rates is associated through a simulation-optimization model that couples MODFLOW-2000 and the hybrid PSOLVER optimization algorithm. Applicability of this procedure is tested on a previously developed groundwater flow model of the Tahtalı Watershed. Successive zone structure solutions are obtained in an additive manner and penalty functions are used in the procedure to obtain realistic and plausible solutions. One of these functions constrains the optimization by forcing the sum of recharge rates for the grid cells that coincide with the Tahtalı Watershed area to be equal to the areal recharge rate determined in the previous modeling by a separate precipitation-runoff model. As a result, a six-zone structure is selected as the best zone structure that represents the areal recharge distribution. Comparison to results of a previous model for the same study area reveals that the proposed procedure significantly improves model performance with respect to calibration statistics. The proposed identification procedure can be thought of as an effective way to determine the recharge zone structure for groundwater flow models, in particular for situations where tangible information about groundwater recharge distribution does not exist. PMID:23746002

  8. Estimates of tracer-based piston-flow ages of groundwater from selected sites: National Water-Quality Assessment Program, 2006-2010

    USGS Publications Warehouse

    Shapiro, Stephanie D.; Plummer, L. Niel; Busenberg, Eurybiades; Widman, Peggy K.; Casile, Gerolamo C.; Wayland, Julian E.; Runkle, Donna L.

    2012-01-01

    Piston-flow age dates were interpreted from measured concentrations of environmental tracers from 812 National Water-Quality Assessment (NAWQA) Program groundwater sites from 27 Study Units across the United States. The tracers of interest include chlorofluorocarbons (CFCs), sulfur hexafluoride (SF6), and tritium/helium-3 (3H/3He). Tracer data compiled for this analysis were collected from 2006 to 2010 from groundwater wells in NAWQA studies, including: * Land-Use Studies (LUS, shallow wells, usually monitoring wells, located in recharge areas under dominant land-use settings), * Major-Aquifer Studies (MAS, wells, usually domestic supply wells, located in principal aquifers and representing the shallow drinking water supply), * Flow System Studies (FSS, networks of clustered wells located along a flowpath extending from a recharge zone to a discharge zone, preferably a shallow stream) associated with Land-Use Studies, and * Reference wells (wells representing groundwater minimally impacted by anthropogenic activities) also associated with Land-Use Studies. Tracer data were evaluated using documented methods and are presented as aqueous concentrations, equivalent atmospheric concentrations (for CFCs and SF6), and tracer-based piston-flow ages. Selected ancillary data, such as redox data, well-construction data, and major dissolved-gas (N2, O2, Ar, CH4, and CO2) data, also are presented. Recharge temperature was inferred using climate data (approximated by mean annual air temperature plus 1°C [MAAT +1°C]) as well as major dissolved-gas data (N2-Ar-based) where available. The N2-Ar-based temperatures showed significantly more variation than the climate-based data, as well as the effects of denitrification and degassing resulting from reducing conditions. The N2-Ar-based temperatures were colder than the climate-based temperatures in networks where recharge was limited to the winter months when evapotranspiration was reduced. The tracer-based piston-flow ages

  9. Considering Barometric Pressure in Groundwater Flow Investigations

    SciTech Connect

    Spane, Frank A. )

    2002-06-18

    Well water-level elevation measurements are commonly used as the basis to delineate groundwater-flow patterns (i.e., flow direction and hydraulic gradient). Barometric pressure fluctuations, however, can have a discernable impact on well water-levels. These barometric effects may lead to erroneous indications of hydraulic head within the aquifer. Total hydraulic head within the aquifer, not well water-level elevation, is the hydrologic parameter for determining groundwater-flow direction and hydraulic gradient conditions. For low-gradient, unconfined aquifer sites exhibiting variable vadose zone characteristics (e.g., thickness, pneumatic diffusivity), barometric pressure fluctuations can also cause actual, temporal changes in lateral flow direction and flow velocity. Discrete water-level measurements used to determine the average or long-term groundwater-flow conditions, therefore, may provide non-representative results. Calculation of the barometric response characteristics for individual wells provides the basis to account for the temporal effects of barometric pressure fluctuations from monitor well measurements, so that average, long-term groundwater-flow pattern behavior can be determined.

  10. Comparison of Estimated Areas Contributing Recharge to Selected Springs in North-Central Florida by Using Multiple Ground-Water Flow Models

    USGS Publications Warehouse

    Shoemaker, W. Barclay; O'Reilly, Andrew M.; Sepulveda, Nicasio; Williams, Stanley A.; Motz, Louis H.; Sun, Qing

    2004-01-01

    Areas contributing recharge to springs are defined in this report as the land-surface area wherein water entering the ground-water system at the water table eventually discharges to a spring. These areas were delineated for Blue Spring, Silver Springs, Alexander Springs, and Silver Glen Springs in north-central Florida using four regional ground-water flow models and particle tracking. As expected, different models predicted different areas contributing recharge. In general, the differences were due to different hydrologic stresses, subsurface permeability properties, and boundary conditions that were used to calibrate each model, all of which are considered to be equally feasible because each model matched its respective calibration data reasonably well. To evaluate the agreement of the models and to summarize results, areas contributing recharge to springs from each model were combined into composite areas. During 1993-98, the composite areas contributing recharge to Blue Spring, Silver Springs, Alexander Springs, and Silver Glen Springs were about 130, 730, 110, and 120 square miles, respectively. The composite areas for all springs remained about the same when using projected 2020 ground-water withdrawals.

  11. Regression modeling of ground-water flow

    USGS Publications Warehouse

    Cooley, R.L.; Naff, R.L.

    1985-01-01

    Nonlinear multiple regression methods are developed to model and analyze groundwater flow systems. Complete descriptions of regression methodology as applied to groundwater flow models allow scientists and engineers engaged in flow modeling to apply the methods to a wide range of problems. Organization of the text proceeds from an introduction that discusses the general topic of groundwater flow modeling, to a review of basic statistics necessary to properly apply regression techniques, and then to the main topic: exposition and use of linear and nonlinear regression to model groundwater flow. Statistical procedures are given to analyze and use the regression models. A number of exercises and answers are included to exercise the student on nearly all the methods that are presented for modeling and statistical analysis. Three computer programs implement the more complex methods. These three are a general two-dimensional, steady-state regression model for flow in an anisotropic, heterogeneous porous medium, a program to calculate a measure of model nonlinearity with respect to the regression parameters, and a program to analyze model errors in computed dependent variables such as hydraulic head. (USGS)

  12. Estimated ground-water recharge from streamflow in Fortymile Wash near Yucca Mountain, Nevada

    SciTech Connect

    Savard, C.S.

    1998-10-01

    The two purposes of this report are to qualitatively document ground-water recharge from stream-flow in Fortymile Wash during the period 1969--95 from previously unpublished ground-water levels in boreholes in Fortymile Canyon during 1982--91 and 1995, and to quantitatively estimate the long-term ground-water recharge rate from streamflow in Fortymile Wash for four reaches of Fortymile Wash (Fortymile Canyon, upper Jackass Flats, lower Jackass Flats, and Amargosa Desert). The long-term groundwater recharge rate was estimated from estimates of the volume of water available for infiltration, the volume of infiltration losses from streamflow, the ground-water recharge volume from infiltration losses, and an analysis of the different periods of data availability. The volume of water available for infiltration and ground-water recharge in the four reaches was estimated from known streamflow in ephemeral Fortymile Wash, which was measured at several gaging station locations. The volume of infiltration losses from streamflow for the four reaches was estimated from a streamflow volume loss factor applied to the estimated streamflows. the ground-water recharge volume was estimated from a linear relation between infiltration loss volume and ground-water recharge volume for each of the four reaches. Ground-water recharge rates were estimated for three different periods of data availability (1969--95, 1983--95, and 1992--95) and a long-term ground-water recharge rate estimated for each of the four reaches.

  13. GROUNDWATER FLOW IN LOW-PERMEABILITY ENVIRONMENTS.

    USGS Publications Warehouse

    Neuzil, C.E.

    1986-01-01

    Certain geologic media are known to have small permeability; subsurface environments composed of these media and lacking well developed secondary permeability have groundwater flow systems with many distinctive characteristics. Moreover, groundwater flow in these environments appears to influence the evolution of certain hydrologic, geologic, and geochemical systems, may affect the accumulation of petroleum and ores, and probably has a role in the structural evolution of parts of the crust. Such environments are also important in the context of waste disposal. This review attempts to synthesize the diverse contributions of various disciplines to the problem of flow in low-permeability environments. Problems hindering analysis are enumerated together with suggested approaches to overcoming them. A common thread running through the discussion is the significance of size- and time-scale limitations of the ability to directly observe flow behavior and make significance of size- and time-scale limitations of the ability to directly observe flow behavior and make measurements of parameters.

  14. Groundwater Flow and Arsenic Biogeochemistry in Bangladesh

    NASA Astrophysics Data System (ADS)

    Harvey, C. F.

    2004-12-01

    Although groundwater in Bangladesh is severely contaminated by arsenic, little is known about the complex transient patterns of groundwater flow that flush solutes from aquifers and carry solutes into the subsurface. Hydrologic modeling results for our field site in the Munshiganj district indicate that groundwater flow is vigorous, flushing the aquifer over time-scales of decades and also introducing solute loads into the aquifer with recharge from rice fields, ponds and rivers. The combined hydrologic and biogeochemical results from our field site imply that the biogeochemistry of the aquifer system may not be in steady-state, and that the net effect of competing processes could either increase or decrease arsenic concentrations over the next decades. Modeling results suggest that irrigation has greatly changed the location, timing and chemical content of recharge to the aquifer, drawing large fluxes of anoxic water into the aquifer during the dry season that may mobilize arsenic from oxides in near-surface sediments.

  15. IN SITU FLOW METERS AROUND A GROUNDWATER CIRCULATION WELL (ABSTRACT)

    EPA Science Inventory

    The primary benefit of groundwater circulation well (GCW) technology is the development of strong vertical flows surrounding the treatment well. The extent of significant vertical flow surrounding a circulation well is difficult to establish from traditional groundwater elevation...

  16. IN SITU FLOW METERS AROUND A GROUNDWATER CIRCULATION WELL (PAPER)

    EPA Science Inventory

    The primary benefit of groundwater circulation well (GCW) technology is the development of strong vertical flows surrounding the treatment well. The extent of significant vertical flow surrounding a circulation well is difficult to establish from traditional groundwater elevation...

  17. Geospatial database of estimates of groundwater discharge to streams in the Upper Colorado River Basin

    USGS Publications Warehouse

    Garcia, Adriana; Masbruch, Melissa D.; Susong, David D.

    2014-01-01

    The U.S. Geological Survey, as part of the Department of the Interior’s WaterSMART (Sustain and Manage America’s Resources for Tomorrow) initiative, compiled published estimates of groundwater discharge to streams in the Upper Colorado River Basin as a geospatial database. For the purpose of this report, groundwater discharge to streams is the baseflow portion of streamflow that includes contributions of groundwater from various flow paths. Reported estimates of groundwater discharge were assigned as attributes to stream reaches derived from the high-resolution National Hydrography Dataset. A total of 235 estimates of groundwater discharge to streams were compiled and included in the dataset. Feature class attributes of the geospatial database include groundwater discharge (acre-feet per year), method of estimation, citation abbreviation, defined reach, and 8-digit hydrologic unit code(s). Baseflow index (BFI) estimates of groundwater discharge were calculated using an existing streamflow characteristics dataset and were included as an attribute in the geospatial database. A comparison of the BFI estimates to the compiled estimates of groundwater discharge found that the BFI estimates were greater than the reported groundwater discharge estimates.

  18. Detecting influential observations in nonlinear regression modeling of groundwater flow

    USGS Publications Warehouse

    Yager, R.M.

    1998-01-01

    Nonlinear regression is used to estimate optimal parameter values in models of groundwater flow to ensure that differences between predicted and observed heads and flows do not result from nonoptimal parameter values. Parameter estimates can be affected, however, by observations that disproportionately influence the regression, such as outliers that exert undue leverage on the objective function. Certain statistics developed for linear regression can be used to detect influential observations in nonlinear regression if the models are approximately linear. This paper discusses the application of Cook's D, which measures the effect of omitting a single observation on a set of estimated parameter values, and the statistical parameter DFBETAS, which quantifies the influence of an observation on each parameter. The influence statistics were used to (1) identify the influential observations in the calibration of a three-dimensional, groundwater flow model of a fractured-rock aquifer through nonlinear regression, and (2) quantify the effect of omitting influential observations on the set of estimated parameter values. Comparison of the spatial distribution of Cook's D with plots of model sensitivity shows that influential observations correspond to areas where the model heads are most sensitive to certain parameters, and where predicted groundwater flow rates are largest. Five of the six discharge observations were identified as influential, indicating that reliable measurements of groundwater flow rates are valuable data in model calibration. DFBETAS are computed and examined for an alternative model of the aquifer system to identify a parameterization error in the model design that resulted in overestimation of the effect of anisotropy on horizontal hydraulic conductivity.

  19. Sublacustrine groundwater discharge in esker aquifers; fully integrated groundwater flow modeling compared with novel field techniques

    NASA Astrophysics Data System (ADS)

    Ala-aho, Pertti; Rossi, Pekka M.; Isokangas, Elina; Kløve, Bjørn

    2015-04-01

    Groundwater (GW) discharge to surface water bodies such as streams, lakes and wetlands can greatly affect their water quantity, quality and related aquatic ecology. Therefore better understanding of GW - surface water interaction is needed in integrated management of water resources. Sublacustrine groundwater discharge (SGD) to lakes was studied in a complex unconfined Rokua esker aquifer system. SGD was studied for 12 lakes in the area to better understand water and solute inputs through lake beds and thereby the role of GW on lake water budget and solute concentrations. The locations and fluxes of SGD were simulated using a fully integrated groundwater flow model HydroGeoSphere. The used hydrological simulator allows water to flow and partition into overland and stream flow, evaporation, infiltration, and subsurface discharge into surface water features in a physically-based way, which was needed in simulating SGD of the complex aquifer system. The model was first calibrated for subsurface hydraulic conductivity in steady state using data of measured long-term average groundwater and lake levels and stream baseflow. The model performance in transient simulations was then examined against recorded hydrographs for lake and groundwater levels and stream flow. After model performance was verified, the simulated locations and fluxes of SGD were extracted from the model and compared with results from three independent field methods: airborne thermal imaging, stable isotope water balance and seepage meter measurements. Airborne thermal imaging was used to infer locations of SGD into lakes based on temperature anomalies at lakes shorelines due to discharging cold groundwater. Isotopic composition (H2 and O18) was analysed for lake water, groundwater and the data was used to estimate SGD flux into lakes. Finally, seepage meter measurements were conducted for one of the lakes to establish both locations and fluxes of SGD in detail. The simulated and field-based estimated

  20. Modeling groundwater flow on massively parallel computers

    SciTech Connect

    Ashby, S.F.; Falgout, R.D.; Fogwell, T.W.; Tompson, A.F.B.

    1994-12-31

    The authors will explore the numerical simulation of groundwater flow in three-dimensional heterogeneous porous media. An interdisciplinary team of mathematicians, computer scientists, hydrologists, and environmental engineers is developing a sophisticated simulation code for use on workstation clusters and MPPs. To date, they have concentrated on modeling flow in the saturated zone (single phase), which requires the solution of a large linear system. they will discuss their implementation of preconditioned conjugate gradient solvers. The preconditioners under consideration include simple diagonal scaling, s-step Jacobi, adaptive Chebyshev polynomial preconditioning, and multigrid. They will present some preliminary numerical results, including simulations of groundwater flow at the LLNL site. They also will demonstrate the code`s scalability.

  1. Multiphase groundwater flow near cooling plutons

    USGS Publications Warehouse

    Hayba, D.O.; Ingebritsen, S.E.

    1997-01-01

    We investigate groundwater flow near cooling plutons with a computer program that can model multiphase flow, temperatures up to 1200??C, thermal pressurization, and temperature-dependent rock properties. A series of experiments examines the effects of host-rock permeability, size and depth of pluton emplacement, single versus multiple intrusions, the influence of a caprock, and the impact of topographically driven groundwater flow. We also reproduce and evaluate some of the pioneering numerical experiments on flow around plutons. Host-rock permeability is the principal factor influencing fluid circulation and heat transfer in hydrothermal systems. The hottest and most steam-rich systems develop where permeability is of the order of 10-15 m2. Temperatures and life spans of systems decrease with increasing permeability. Conduction-dominated systems, in which permeabilities are ???10-16m2, persist longer but exhibit relatively modest increases in near-surface temperatures relative to ambient conditions. Pluton size, emplacement depth, and initial thermal conditions have less influence on hydrothermal circulation patterns but affect the extent of boiling and duration of hydrothermal systems. Topographically driven groundwater flow can significantly alter hydrothermal circulation; however, a low-permeability caprock effectively decouples the topographically and density-driven systems and stabilizes the mixing interface between them thereby defining a likely ore-forming environment.

  2. Use of Chemical and Isotopic Tracers for Estimating Ground-Water Recharge, Flow Paths, and Residence Times in the Middle San Pedro Basin, Southeast Arizona

    NASA Astrophysics Data System (ADS)

    Adkins, C. B.; McIntosh, J.; Eastoe, C.; Dickinson, J.

    2008-12-01

    Ground water is often the primary source of water for rapidly growing populations in the semi-arid southwestern United States. In addition, ground-water discharge to streams sustains wildlife in riparian areas. Improved understanding of the sources of ground water, recharge areas, flow paths, and water quality of basin aquifer systems is needed to assess water availability and develop effective water management policies. This study analyzes variations of major ion (Ca, Na, K, Mg, Sr, Fe, Si, Zn, F, Cl, Br, NO3, SO4) and isotope (18O, 2H, 3H, 34S, 13C, 14C) chemistry of ground water, surface waters and precipitation with in conjunction with hydrogeologic data (e.g. hydraulic head, subsurface structure, and stratigraphy) to infer recharge areas, mixing of water sources, and residence times of ground water within the middle San Pedro watershed in southeastern Arizona. The San Pedro basin is bound by crystalline and carbonate rocks of the Whetstone and Rincon Mountains on the west and by crystalline rocks of the Dragoon Mountains to the east. Differences in mineral assemblages of these mountain blocks impart distinct chemical signatures in ground waters through mineral weathering. Potentially, these differences in water chemistry can serve as chemical tracers for identifying ground-water flow paths and mixing relations. Ground-water chemistry variations suggest compartmentalization of waters into an upper and lower alluvial aquifer system comprised of permeable sands and gravels ranging in depth from ten to over one thousand feet in the basin center; the units are separated by confining units of silt and clay in the basin center. Variations include higher fluoride (up to 8 ppm) near the Dragoon Mountains, higher chloride (up to 54 ppm) near the Whetstone Mountains, and higher sulfate (up to 750 ppm) concentrations in both upper and lower sands and gravels owing to interaction with thick Permian or Neocene evaporites. Chloride is generally lower (less than 8 ppm) in the

  3. Effects of intraborehole flow on groundwater age distribution

    USGS Publications Warehouse

    Zinn, B.A.; Konikow, L.F.

    2007-01-01

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

  4. Connections between groundwater flow and transpiration partitioning.

    PubMed

    Maxwell, Reed M; Condon, Laura E

    2016-07-22

    Understanding freshwater fluxes at continental scales will help us better predict hydrologic response and manage our terrestrial water resources. The partitioning of evapotranspiration into bare soil evaporation and plant transpiration remains a key uncertainty in the terrestrial water balance. We used integrated hydrologic simulations that couple vegetation and land-energy processes with surface and subsurface hydrology to study transpiration partitioning at the continental scale. Both latent heat flux and partitioning are connected to water table depth, and including lateral groundwater flow in the model increases transpiration partitioning from 47 ± 13 to 62 ± 12%. This suggests that lateral groundwater flow, which is generally simplified or excluded in Earth system models, may provide a missing link for reconciling observations and global models of terrestrial water fluxes. PMID:27463671

  5. Connections between groundwater flow and transpiration partitioning

    NASA Astrophysics Data System (ADS)

    Maxwell, Reed M.; Condon, Laura E.

    2016-07-01

    Understanding freshwater fluxes at continental scales will help us better predict hydrologic response and manage our terrestrial water resources. The partitioning of evapotranspiration into bare soil evaporation and plant transpiration remains a key uncertainty in the terrestrial water balance. We used integrated hydrologic simulations that couple vegetation and land-energy processes with surface and subsurface hydrology to study transpiration partitioning at the continental scale. Both latent heat flux and partitioning are connected to water table depth, and including lateral groundwater flow in the model increases transpiration partitioning from 47 ± 13 to 62 ± 12%. This suggests that lateral groundwater flow, which is generally simplified or excluded in Earth system models, may provide a missing link for reconciling observations and global models of terrestrial water fluxes.

  6. Groundwater flow systems in mountainous terrain, 1. Numerical modeling technique

    NASA Astrophysics Data System (ADS)

    Forster, Craig; Smith, Leslie

    1988-07-01

    A coupled model of fluid flow and heat transfer is developed to characterize steady groundwater flow within a mountain massif. A coupled model is necessary because high-relief terrain can enhance groundwater flow to depths where elevated temperatures are encountered. A wide range in water table form and elevation expected in high-relief terrain is accommodated using a free-surface method. This approach allows us to examine the influence of thermal conditions on the patterns and rates of groundwater flow and the position of the water table. Vertical fluid flow is assumed to occur within the unsaturated zone to provide a simple basis for modeling advective heat transfer above the water table. This approach ensures that temperatures at the water table, and throughout the domain, are consistent with temperature conditions specified at the bedrock surface. Conventional free-surface methods provide poor estimates of the water table configuration in high-relief terrain. A modified free-surface approach is introduced to accommodate recharge at upper elevations on the seepage face, in addition to recharge at the free surface.

  7. Groundwater flow as a cooling agent of the continental lithosphere

    NASA Astrophysics Data System (ADS)

    Kooi, Henk

    2016-03-01

    Groundwater that flows through the outer shell of the Earth as part of the hydrologic cycle influences the distribution of heat and, thereby, the temperature field in the Earth’s crust. Downward groundwater flow in recharge areas lowers crustal temperatures, whereas upward flow in discharge areas tends to raise temperatures relative to a purely conductive geothermal regime. Here I present numerical simulations of generalized topography-driven groundwater flow. The simulations suggest that groundwater-driven convective cooling exceeds groundwater-driven warming of the Earth’s crust, and hence that groundwater flow systems cause net temperature reductions of groundwater basins. Moreover, the simulations demonstrate that this cooling extends into the underlying crust and lithosphere. I find that horizontal components of groundwater flow play a central role in this net subsurface cooling by conveying relatively cold water to zones of upward groundwater flow. The model calculations suggest that the crust and lithosphere beneath groundwater basins can cool by several tens of degrees Celsius where groundwater flows over large distances in basins that consist of crustal rock. In contrast, groundwater-induced cooling is small in unconsolidated sedimentary settings, such as deltas.

  8. Estimated Water Flows in 2005: United States

    SciTech Connect

    Smith, C A; Belles, R D; Simon, A J

    2011-03-16

    Flow charts depicting water use in the United States have been constructed from publicly available data and estimates of water use patterns. Approximately 410,500 million gallons per day of water are managed throughout the United States for use in farming, power production, residential, commercial, and industrial applications. Water is obtained from four major resource classes: fresh surface-water, saline (ocean) surface-water, fresh groundwater and saline (brackish) groundwater. Water that is not consumed or evaporated during its use is returned to surface bodies of water. The flow patterns are represented in a compact 'visual atlas' of 52 state-level (all 50 states in addition to Puerto Rico and the Virgin Islands) and one national water flow chart representing a comprehensive systems view of national water resources, use, and disposition.

  9. Groundwater velocities at the Nevada Test Site: {sup 14}Carbon-based estimates

    SciTech Connect

    Chapman, J.B.; Hershey, R.L.; Lyles, B.F.

    1995-07-01

    Chemical and isotopic data can be used to constrain and validate groundwater flow models. This study examines probable groundwater flowpaths at the Nevada Test Site (NTS) and estimates groundwater velocities for these flowpaths using water chemistry and carbon isotopes. These velocities are provided for comparison to velocities calculated by a numerical flow model developed by GeoTrans, Inc. Similar to numerical flow models, models of chemical and isotopic evolution are not unique; any number of combinations of reactions can simulate evolution from one water to another, but are no guarantee that the simulation is correct. Knowledge of the hydrology, mineralogy, and chemistry must be combined to produce feasible evolutionary paths.

  10. Coupled Groundwater and Heat Flow in the Tahoe Basin Region

    NASA Astrophysics Data System (ADS)

    Trask, J. C.

    2002-12-01

    We propose that recent developments in available temperature probe technology and improvements in appropriate modeling software, together with the advent of desktop high-speed computing, have enabled the thermal analysis approach to be an inexpensive, robust, and practical way to constrain groundwater flow fields in a wide variety of environments. We present field measurements and numerical models of coupled heat and groundwater flow in the Tahoe Basin region. In montane regions such as the Tahoe Basin, steep topography provides a driving force for deep groundwater flow. Deep groundwater flow re-routes subsurface heat flow, impacting temperature gradients to depth, including the shallow subsurface (<100m depth). In the Tahoe Basin region, the magnitude of deep groundwater flow on the areal or regional scale has been largely unknown. We present examples of borehole temperature profiles that constrain possible areal groundwater flow patterns, including the magnitude of flow beneath the bottom of boreholes probed.

  11. Groundwater Recharge and Flow Processes in Taihang Mountains, a Semi-humid Region, North China

    NASA Astrophysics Data System (ADS)

    Sakakibara, Koichi; Tsujimura, Maki; Song, Xianfang; Zhang, Jie

    2015-04-01

    Groundwater flow/recharge variations in time and space are crucial for effective water management especially in semi-arid and semi-humid regions. In order to reveal comprehensive groundwater flow/recharge processes in a catchment with a large topographical relief and seasonal hydrological variations, intensive field surveys were undertaken at 4 times in different seasons (June 2011, August 2012, November 2012, February 2014) in the Wangkuai watershed, Taihang mountains, which is a main groundwater recharge area of the North China Plain. The groundwater, spring, stream water and reservoir water were taken, and inorganic solute constituents and stable isotopes of oxygen-18 and deuterium were determined on all water samples. Also, the stream flow rate and the depth of groundwater table were observed. The stable isotopic compositions and inorganic solute constituents in the groundwater are depleted and shown similar values as those of the surface water at the mountain-plain transitional area. Additionally, the groundwater in the vicinity of the Wangkuai Reservoir presents clearly higher stable isotopic compositions and lower d-excess than those of the stream water, indicating the groundwater around the reservoir is affected by evaporation same as the Wangkuai Reservoir itself. Hence, the surface water in the mountain-plain transitional area and Wangkuai Reservoir are principal groundwater recharge sources. An inversion analysis and simple mixing model were applied in the Wangkuai watershed using stable isotopes of oxygen-18 and deuterium to construct a groundwater flow model. The model shows that multi-originated groundwater flows from upstream to downstream along topography with certain mixing. In addition, the groundwater recharge occurs dominantly at the altitude from 421 m to 953 m, and the groundwater recharge rate by the Wangkuai Reservoir is estimated to be 2.4 % of the total groundwater recharge in the Wangkuai watershed. Therefore, the stream water and

  12. Patterns in groundwater chemistry resulting from groundwater flow

    NASA Astrophysics Data System (ADS)

    Stuyfzand, Pieter J.

    Groundwater flow influences hydrochemical patterns because flow reduces mixing by diffusion, carries the chemical imprints of biological and anthropogenic changes in the recharge area, and leaches the aquifer system. Global patterns are mainly dictated by differences in the flux of meteoric water passing through the subsoil. Within individual hydrosomes (water bodies with a specific origin), the following prograde evolution lines (facies sequence) normally develop in the direction of groundwater flow: from strong to no fluctuations in water quality, from polluted to unpolluted, from acidic to basic, from oxic to anoxic-methanogenic, from no to significant base exchange, and from fresh to brackish. This is demonstrated for fresh coastal-dune groundwater in the Netherlands. In this hydrosome, the leaching of calcium carbonate as much as 15m and of adsorbed marine cations (Na+, K+, and Mg2+) as much as 2500m in the flow direction is shown to correspond with about 5000yr of flushing since the beach barrier with dunes developed. Recharge focus areas in the dunes are evidenced by groundwater displaying a lower prograde quality evolution than the surrounding dune groundwater. Artificially recharged Rhine River water in the dunes provides distinct hydrochemical patterns, which display groundwater flow, mixing, and groundwater ages. Résumé Les écoulements souterrains influencent les différents types hydrochimiques, parce que l'écoulement réduit le mélange par diffusion, porte les marques chimiques de changements biologiques et anthropiques dans la zone d'alimentation et lessive le système aquifère. Ces types dans leur ensemble sont surtout déterminés par des différences dans le flux d'eau météorique traversant le sous-sol. Dans les "hydrosomes" (masses d'eau d'origine déterminée), les lignes marquant une évolution prograde (séquence de faciès) se développent normalement dans la direction de l'écoulement souterrain : depuis des fluctuations fortes de la

  13. Groundwater flow and solute movement to drain laterals, western San Joaquin Valley, California: 2. Quantitative hydrologic assessment

    USGS Publications Warehouse

    Fio, John L.; Deverel, S.J.

    1991-01-01

    Groundwater flow modeling was used to quantitatively assess the hydrologic processes affecting ground water and solute movement to drain laterals. Modeling results were used to calculate the depth distribution of groundwater flowing into drain laterals at 1.8 m (drain lateral 1) and 2.7 m (drain lateral 2) below land surface. The simulations indicated that under nonirrigated conditions about 89% of the flow in drain lateral 2 was from groundwater originating from depths greater than 6 m below land surface. The deep groundwater has higher selenium concentrations than shallow groundwater. Simulation of irrigated conditions indicates that as recharge (deep percolation) increases, the proportional contribution of deep groundwater to drain lateral flow decreases. Groundwater flow paths and travel times estimated from the simulation results indicate that groundwater containing high concentrations of selenium (greater than 780 μg L−1) probably will continue to enter drain lateral 2 for decades.

  14. A Comparison of Three Stochastic Approaches for Parameter Estimation and Prediction of Steady-State Groundwater Flow: Nonlocal Moment Equations and Monte Carlo Method Coupled with Ensemble Kalman Filter and Geostatistical Stochastic Inversion.

    NASA Astrophysics Data System (ADS)

    Morales-Casique, E.; Briseño-Ruiz, J. V.; Hernández, A. F.; Herrera, G. S.; Escolero-Fuentes, O.

    2014-12-01

    We present a comparison of three stochastic approaches for estimating log hydraulic conductivity (Y) and predicting steady-state groundwater flow. Two of the approaches are based on the data assimilation technique known as ensemble Kalman filter (EnKF) and differ in the way prior statistical moment estimates (PSME) (required to build the Kalman gain matrix) are obtained. In the first approach, the Monte Carlo method is employed to compute PSME of the variables and parameters; we denote this approach by EnKFMC. In the second approach PSME are computed through the direct solution of approximate nonlocal (integrodifferential) equations that govern the spatial conditional ensemble means (statistical expectations) and covariances of hydraulic head (h) and fluxes; we denote this approach by EnKFME. The third approach consists of geostatistical stochastic inversion of the same nonlocal moment equations; we denote this approach by IME. In addition to testing the EnKFMC and EnKFME methods in the traditional manner that estimate Y over the entire grid, we propose novel corresponding algorithms that estimate Y at a few selected locations and then interpolate over all grid elements via kriging as done in the IME method. We tested these methods to estimate Y and h in steady-state groundwater flow in a synthetic two-dimensional domain with a well pumping at a constant rate, located at the center of the domain. In addition, to evaluate the performance of the estimation methods, we generated four unconditional different realizations that served as "true" fields. The results of our numerical experiments indicate that the three methods were effective in estimating h, reaching at least 80% of predictive coverage, although both EnKF were superior to the IME method. With respect to estimating Y, the three methods reached similar accuracy in terms of the mean absolute value error. Coupling the EnKF methods with kriging to estimate Y reduces to one fourth the CPU time required for data

  15. Groundwater Flow Model for Taos, New Mexico

    NASA Astrophysics Data System (ADS)

    Burck, P. W.; Barroll, P. W.; Core, A. B.; Rappuhn, D.

    2003-12-01

    The New Mexico Office of the State Engineer - Hydrology Bureau (OSE) has developed a regional groundwater flow model for Taos, New Mexico. The MODFLOW 2000 model will serve as a tool to evaluate alternatives in settlement negotiations in an on-going water rights adjudication. If current settlement negotiations fail, it is conceivable that the model might be used in support of litigation. OSE produced the model in cooperation with technical representatives of the various parties to the adjudication. Regional hydrogeologic data including well records, aquifer test results, stream flow measurements and seepage studies have been shared relatively freely among the parties. A recent deep drilling program conducted in conjunction with the negotiation effort has added substantially to the hydrogeologic data set. Among the hydrologic processes simulated by the model are mountain front recharge; areal recharge from precipitation; evapotranspiration; discharge from springs; river and stream flow; accretions to groundwater from irrigation return flow, seepage from acequias, canals, and ditches, and deep percolation; and pumping by municipal entities and mutual domestic water users associations. The resulting model files are available for all parties to review and evaluate. Comments are assessed and many have resulted in significant improvements to the model. At this stage, however, it is unclear whether adopting this cooperative approach will increase the likelihood of model acceptance by the parties.

  16. Guidelines for Evaluating Ground-Water Flow Models

    USGS Publications Warehouse

    Reilly, Thomas E.; Harbaugh, Arlen W.

    2004-01-01

    Ground-water flow modeling is an important tool frequently used in studies of ground-water systems. Reviewers and users of these studies have a need to evaluate the accuracy or reasonableness of the ground-water flow model. This report provides some guidelines and discussion on how to evaluate complex ground-water flow models used in the investigation of ground-water systems. A consistent thread throughout these guidelines is that the objectives of the study must be specified to allow the adequacy of the model to be evaluated.

  17. Using Groundwater Age and Other Isotopic Signatures to Delineate Groundwater Flow and Stratification

    SciTech Connect

    Moran, J E; Hudson, G B

    2005-08-31

    Isotopic tracers, such as stable isotopes of the water molecule and tritium, have been used in investigations of groundwater flow and transport and recharge water source for several decades. While these data can place hard constraints on groundwater flow rates, the degree of vertical flow between aquifers and across aquitards, and recharge source area(s), they are rarely used, even for validation, in conceptual or numerical models of groundwater flow. The Groundwater Ambient Monitoring and Assessment Program, sponsored by the California State Water Resources Control Board, and carried out in collaboration with the U.S. Geological Survey, has provided the means to gather an unprecedented number of tritium-helium groundwater ages in the basins of California. As the examples below illustrate, a collection of groundwater ages in a basin allows delineation of recharge areas (youngest ages), bulk flow rates and flowpaths, as well as a means of assessing susceptibility to anthropogenic contaminants.

  18. Assessment of factors influencing groundwater-level change using groundwater flow simulation, considering vertical infiltration from rice-planted and crop-rotated paddy fields in Japan

    NASA Astrophysics Data System (ADS)

    Iwasaki, Yumi; Nakamura, Kimihito; Horino, Haruhiko; Kawashima, Shigeto

    2014-12-01

    Assessing factors that influence groundwater levels such as land use and pumping strategy, is essential to adequately manage groundwater resources. A transient numerical model for groundwater flow with infiltration was developed for the Tedori River alluvial fan (140 km2), Japan. The main water input into the groundwater body in this area is irrigation water, which is significantly influenced by land use, namely paddy and upland fields. The proposed model consists of two models, a one-dimensional (1-D) unsaturated-zone water flow model (HYDRUS-1D) for estimating groundwater recharge and a 3-D groundwater flow model (MODFLOW). Numerical simulation of groundwater flow from October 1975 to November 2009 was performed to validate the model. Simulation revealed seasonal groundwater level fluctuations, affected by paddy irrigation management. However, computational accuracy was limited by the spatiotemporal data resolution of the groundwater use. Both annual groundwater levels and recharge during the irrigation periods from 1975 to 2009 showed long-term decreasing trends. With the decline in rice-planted paddy field area, groundwater recharge cumulatively decreased to 61 % of the peak in 1977. A paddy-upland crop-rotation system could decrease groundwater recharge to 73-98 % relative to no crop rotation.

  19. Analysis of groundwater flow in mountainous, headwater catchments with permafrost

    NASA Astrophysics Data System (ADS)

    Evans, Sarah G.; Ge, Shemin; Liang, Sihai

    2015-12-01

    Headwater catchments have a direct impact on the water resources of downstream lowland regions as they supply freshwater in the form of surface runoff and discharging groundwater. Often, these mountainous catchments contain expansive permafrost that may alter the natural topographically controlled groundwater flow system. As permafrost could degrade with climate change, it is imperative to understand the effect of permafrost on groundwater flow in headwater catchments. This study characterizes groundwater flow in mountainous headwater catchments and evaluates the effect of permafrost in the context of climate change on groundwater movement using a three-dimensional, finite element, hydrogeologic model. The model is applied to a representative headwater catchment on the Qinghai-Tibet Plateau, China. Results from the model simulations indicate that groundwater contributes significantly to streams in the form of baseflow and the majority of groundwater flow is from the shallow aquifer above the permafrost, disrupting the typical topographically controlled flow pattern observed in most permafrost-free headwater catchments. Under a warming scenario where mean annual surface temperature is increased by 2°C, reducing the areal extent of permafrost in the catchment, groundwater contribution to streamflow may increase three-fold. These findings suggest that, in headwater catchments, permafrost has a large influence on groundwater flow and stream discharge. Increased annual air temperatures may increase groundwater discharge to streams, which has implications for ecosystem health and the long-term availability of water resources to downstream regions.

  20. Ground-water flow and solute transport at a municipal landfill site on Long Island, New York; Part 2, Simulation of ground-water flow

    USGS Publications Warehouse

    Wexler, E.J.; Maus, P.E.

    1988-01-01

    Data on the hydrogeology of a 26-sq-mi area surrounding the Brookhaven landfill site in central Suffolk County were collected as part of a hydrologic investigation of solute transport from the site. These data were used to develop a steady-state groundwater flow model of the upper glacial (water table) aquifer in the area. The model accounts for the leakage through confining units underlying the aquifer, seepage to streams, recharge from precipitation, and pumpage and redistribution of water. Refined estimates of aquifer and confining-unit properties were obtained through model calibrations. Water table altitudes generated by the calibrated model were used to determine groundwater velocities and probable flow paths in the vicinity of the site under long-term average hydrologic conditions. Groundwater velocities and probable flow paths in the study area were calculated from simulated water table altitudes generated by the calibrated flow model. Groundwater at the center of the site flows southeastward at a velocity of 1.1 ft/d. The report is the second in a three part series describing the hydrologic conditions and groundwater quality, groundwater flow, and solute transport in the vicinity of the Brookhaven landfill. (USGS)

  1. Uncertainty in global groundwater storage estimates in a Total Groundwater Stress framework

    NASA Astrophysics Data System (ADS)

    Richey, Alexandra S.; Thomas, Brian F.; Lo, Min-Hui; Famiglietti, James S.; Swenson, Sean; Rodell, Matthew

    2015-07-01

    Groundwater is a finite resource under continuous external pressures. Current unsustainable groundwater use threatens the resilience of aquifer systems and their ability to provide a long-term water source. Groundwater storage is considered to be a factor of groundwater resilience, although the extent to which resilience can be maintained has yet to be explored in depth. In this study, we assess the limit of groundwater resilience in the world's largest groundwater systems with remote sensing observations. The Total Groundwater Stress (TGS) ratio, defined as the ratio of total storage to the groundwater depletion rate, is used to explore the timescales to depletion in the world's largest aquifer systems and associated groundwater buffer capacity. We find that the current state of knowledge of large-scale groundwater storage has uncertainty ranges across orders of magnitude that severely limit the characterization of resilience in the study aquifers. Additionally, we show that groundwater availability, traditionally defined as recharge and redefined in this study as total storage, can alter the systems that are considered to be stressed versus unstressed. We find that remote sensing observations from NASA's Gravity Recovery and Climate Experiment can assist in providing such information at the scale of a whole aquifer. For example, we demonstrate that a groundwater depletion rate in the Northwest Sahara Aquifer System of 2.69 ± 0.8 km3/yr would result in the aquifer being depleted to 90% of its total storage in as few as 50 years given an initial storage estimate of 70 km3.

  2. Empirical estimation of groundwater quality changes using remote sensing

    NASA Astrophysics Data System (ADS)

    Gibbons, A.; Thomas, B. F.; Famiglietti, J. S.

    2015-12-01

    Recent groundwater availability studies estimate large-scale aquifer depletion rates and aquifer stress using monthly water storage variations from NASA's Gravity Recovery and Climate Experiment (GRACE) mission. To further evaluate available groundwater resources, assessing potability of groundwater is necessary. Statistical relationships are initially developed at individual well locations to discern our ability to predict groundwater geochemistry as a function of groundwater levels. Next, up-scaled multivariate relationships to estimate total dissolved solid (TDS) concentrations as a function of GRACE-derived subsurface storage anomalies, dominant land use, and other physical parameters are developed in two important aquifer systems in the United States: the High Plains aquifer and the Central Valley aquifer. A goodness of fit test was performed to evaluate model strength. Results demonstrate the potential to characterize global groundwater potability variations using remote sensing.

  3. Tracing Groundwater Flow on a Barrier Island in the North-east Gulf of Mexico

    NASA Astrophysics Data System (ADS)

    Reide Corbett, D.; Dillon, K.; Burnett, W.

    2000-08-01

    Groundwater flow on St. George Island, a barrier island in the north-east Gulf of Mexico, was monitored downfield from wastewater systems using artificial tracer techniques. Sulphur hexafluoride and fluorescein dye were used to determine groundwater flow velocity, hydraulic conductivity, and dispersivity at selected sites on the island. Monthly hydraulic head measurements illustrate the aquifer's dependence on recharge associated with rainfall. However, during periods of little rain, tidal stage also influences the direction and magnitude of groundwater flow within approximately 30 m of the waters edge. Estimated hydraulic conductivities ranged from 3 to 180 m day -1, with an overall estimate of 36 m day -1. Groundwater tracers showed very little dilution and calculated longitudinal dispersivities were approximately 0·1 to 0·5 m, which is in the same range as previous studies of sandy aquifers. The total groundwater flux into adjacent Apalachicola Bay was also evaluated using two independent techniques. Darcy's law was applied using an estimated cross sectional area and the experimental horizontal transport rates to estimate the volumetric flow. In addition, a simple water balance calculation was used, which accounted for all known sources and sinks of water to the aquifer. The two independent approaches agreed very well, with an estimated groundwater flux from the surficial aquifer to Apalachicola Bay between 1-9×10 6 m 3 year -1.

  4. Application of chlorofluorocarbons (CFCs) to estimate the groundwater age at a headwater wetland in Ichikawa City, Chiba Prefecture, Japan

    NASA Astrophysics Data System (ADS)

    Han, Zhiwei; Tang, Changyuan; Piao, Jingqiu; Li, Xing; Cao, Yingjie; Matsumaru, Touma; Zhang, Chipeng

    2014-09-01

    To delineate the groundwater flow system in a basin, the groundwater age was estimated by analyzing chlorofluorocarbons (CFC-11, CFC-12 and CFC-113) in a typical headwater wetland in Ichikawa, Japan. Feasibility of groundwater dating by CFCs was assessed comprehensively based on the concentrations of NO3 -, SO4 2-, Fe2+ and dissolved CH4 in the groundwater, because the CFCs would be degraded under the reduction condition available in a wetland. It was found that the CFC-11 apparent age was much older than that estimated by other CFC species. It showed that CFC-12 and CFC-113 were suitable tracers for groundwater dating because of their stability in the wetland environment. Furthermore, the mixture of groundwater with different age was discussed by CFC-12 and CFC-113 based on the binary mixing model and piston-flow model. As a result, the apparent age of groundwater in the study area is in the range of 38-48 years.

  5. CONCEPTUAL FRAMEWORK FOR REGRESSION MODELING OF GROUND-WATER FLOW.

    USGS Publications Warehouse

    Cooley, Richard L.

    1985-01-01

    The author examines the uses of ground-water flow models and which classes of use require treatment of stochastic components. He then compares traditional and stochastic procedures for modeling actual (as distinguished from hypothetical) systems. Finally, he examines the conceptual basis and characteristics of the regression approach to modeling ground-water flow.

  6. Effects of linking a soil-water-balance model with a groundwater-flow model

    USGS Publications Warehouse

    Stanton, Jennifer S.; Ryter, Derek W.; Peterson, Steven M.

    2013-01-01

    A previously published regional groundwater-flow model in north-central Nebraska was sequentially linked with the recently developed soil-water-balance (SWB) model to analyze effects to groundwater-flow model parameters and calibration results. The linked models provided a more detailed spatial and temporal distribution of simulated recharge based on hydrologic processes, improvement of simulated groundwater-level changes and base flows at specific sites in agricultural areas, and a physically based assessment of the relative magnitude of recharge for grassland, nonirrigated cropland, and irrigated cropland areas. Root-mean-squared (RMS) differences between the simulated and estimated or measured target values for the previously published model and linked models were relatively similar and did not improve for all types of calibration targets. However, without any adjustment to the SWB-generated recharge, the RMS difference between simulated and estimated base-flow target values for the groundwater-flow model was slightly smaller than for the previously published model, possibly indicating that the volume of recharge simulated by the SWB code was closer to actual hydrogeologic conditions than the previously published model provided. Groundwater-level and base-flow hydrographs showed that temporal patterns of simulated groundwater levels and base flows were more accurate for the linked models than for the previously published model at several sites, particularly in agricultural areas.

  7. Simulation of groundwater flow and effects of groundwater irrigation on stream base flow in the Elkhorn and Loup River Basins, Nebraska, 1895-2055-Phase Two

    USGS Publications Warehouse

    Stanton, Jennifer S.; Peterson, Steven M.; Fienen, Michael N.

    2010-01-01

    Regional groundwater-flow simulations for a 30,000-square-mile area of the High Plains aquifer, referred to collectively as the Elkhorn-Loup Model, were developed to predict the effects of groundwater irrigation on stream base flow in the Elkhorn and Loup River Basins, Nebraska. Simulations described the stream-aquifer system from predevelopment through 2005 [including predevelopment (pre-1895), early development (1895-1940), and historical development (1940 through 2005) conditions] and future hypothetical development conditions (2006 through 2033 or 2055). Predicted changes to stream base flow that resulted from simulated changes to groundwater irrigation will aid development of long-term strategies for management of hydrologically connected water supplies. The predevelopment through 2005 simulation was calibrated using an automated parameter-estimation method to optimize the fit to pre-1940 groundwater levels and base flows, 1945 through 2005 decadal groundwater-level changes, and 1940 through 2005 base flows. The calibration results of the pre-1940 period indicated that 81 percent of the simulated groundwater levels were within 30 feet of the measured water levels. The results did not indicate large areas of simulated groundwater levels that were biased too high or too low, indicating that the simulation generally captures the regional trends. Calibration results using 1945 through 2005 decadal groundwater-level changes indicated that a majority of the simulated groundwater-level changes were within 5 feet of the changes calculated from measured groundwater levels. Simulated groundwater-level rises generally were smaller than measured rises near surface-water irrigation districts. Simulated groundwater-level declines were larger than measured declines in several parts of the study area having large amounts of irrigated crops. Base-flow trends and volumes generally were reproduced by the simulation at most sites. Exceptions include downward trends of simulated

  8. Estimating contributions of nitrate and herbicides from groundwater to headwater streams, northern Atlantic Coastal Plain, USA

    USGS Publications Warehouse

    Ator, Scott; Denver, Judith

    2012-01-01

    Groundwater transport often complicates understanding of surface-water contamination. We estimated the regional flux of nitrate and selected herbicides from groundwater to nontidal headwater streams of the Atlantic Coastal Plain (New Jersey through North Carolina) based on late-winter or spring base-flow samples from 174 streams. Sampled streams were selected randomly, and flux estimates are based on resulting population estimates rather than on empirical models, which have been used previously for similar estimates. Base-flow flux in the estimated 8,834 headwater streams of the study area are an estimated 21,200 kg/day of nitrate (as N) and 5.83, 0.565, and 20.7 kg/day of alachlor, atrazine, and metolachlor (and selected degradates), respectively. Base-flow flux of alachlor and metolachlor is <3% of the total base-flow flux of those compounds plus degradates. Base-flow flux of nitrate and herbicides as a percentage of applications is typically highest in well-drained areas and lowest in areas with abundant poor drainage and anoxic conditions. In Coastal Plain watersheds of Albemarle and Pamlico Sounds, <2% of applied nitrogen reaches headwater streams as base flow. On the Delmarva Peninsula part of the Chesapeake Bay watershed, however, more than 10% of such applications are transported through groundwater to streams, and base-flow nitrate flux represents 70% of total nitrogen flux in headwater streams.

  9. Groundwater flow system under a rapidly urbanizing coastal city as determined by hydrogeochemistry

    NASA Astrophysics Data System (ADS)

    Kagabu, Makoto; Shimada, Jun; Delinom, Robert; Tsujimura, Maki; Taniguchi, Makoto

    2011-01-01

    In the Jakarta area (Indonesia), excessive groundwater pumping due to the rapidly increasing population has caused groundwater-related problems such as brackish water contamination in coastal areas and land subsidence. In this study, we adopted multiple hydrogeochemical techniques to demonstrate the groundwater flow system in the Jakarta area. Although almost all groundwater existing in the Jakarta basin is recharged at similar elevations, the water quality and residence time demonstrates a clear difference between the shallow and deep aquifers. Due to the rapid decrease in the groundwater potential in urban areas, we found that the seawater intrusion and the shallow and deep groundwaters are mixing, a conclusion confirmed by major ions, Br -:Cl - ratios, and chlorofluorocarbon (CFC)-12 analysis. Spring water and groundwater samples collected from the southern mountainside area show younger age characteristics with high concentrations of 14C and Ca-HCO 3 type water chemistry. We estimated the residence times of these groundwaters within 45 years under piston flow conditions by tritium analysis. Also, these groundwater ages can be limited to 20-30 years with piston flow evaluated by CFCs. Moreover, due to the magnitude of the CFC-12 concentration, we can use a pseudo age indicator in this field study, because we found a positive correlation between the major type of water chemistry and the CFC-12 concentration.

  10. Groundwater recharge in Wisconsin--Annual estimates for 1970-99 using streamflow data

    USGS Publications Warehouse

    Gebert, Warren A.; Walker, John F.; Hunt, Randall J.

    2011-01-01

    The groundwater component of streamflow is important because it is indicative of the sustained flow of a stream during dry periods, is often of better quality, and has a smaller range of temperatures, than surface contributions to streamflow. All three of these characteristics are important to the health of aquatic life in a stream. If recharge to the aquifers is to be preserved or enhanced, it is important to understand the present partitioning of total streamflow into base flow and stormflow. Additionally, an estimate of groundwater recharge is important for understanding the flows within a groundwater system-information important for water availability/sustainability or other assessments. The U.S. Geological Survey operates numerous continuous-record streamflow-gaging stations (Hirsch and Norris, 2001), which can be used to provide estimates of average annual base flow. In addition to these continuous record sites, Gebert and others (2007) showed that having a few streamflow measurements in a basin can appreciably reduce the error in a base-flow estimate for that basin. Therefore, in addition to the continuous-record gaging stations, a substantial number of low-flow partial-record sites (6 to 15 discharge measurements) and miscellaneous-measurement sites (1 to 3 discharge measurements) that were operated during 1964-90 throughout the State were included in this work to provide additional insight into spatial distribution of annual base flow and, in turn, groundwater recharge.

  11. Combining groundwater quality analysis and a numerical flow simulation for spatially establishing utilization strategies for groundwater and surface water in the Pingtung Plain

    NASA Astrophysics Data System (ADS)

    Jang, Cheng-Shin; Chen, Ching-Fang; Liang, Ching-Ping; Chen, Jui-Sheng

    2016-02-01

    Overexploitation of groundwater is a common problem in the Pingtung Plain area of Taiwan, resulting in substantial drawdown of groundwater levels as well as the occurrence of severe seawater intrusion and land subsidence. Measures need to be taken to preserve these valuable groundwater resources. This study seeks to spatially determine the most suitable locations for the use of surface water on this plain instead of extracting groundwater for drinking, irrigation, and aquaculture purposes based on information obtained by combining groundwater quality analysis and a numerical flow simulation assuming the planning of manmade lakes and reservoirs to the increase of water supply. The multivariate indicator kriging method is first used to estimate occurrence probabilities, and to rank townships as suitable or unsuitable for groundwater utilization according to water quality standards for drinking, irrigation, and aquaculture. A numerical model of groundwater flow (MODFLOW) is adopted to quantify the recovery of groundwater levels in townships after model calibration when groundwater for drinking and agricultural demands has been replaced by surface water. Finally, townships with poor groundwater quality and significant increases in groundwater levels in the Pingtung Plain are prioritized for the groundwater conservation planning based on the combined assessment of groundwater quality and quantity. The results of this study indicate that the integration of groundwater quality analysis and the numerical flow simulation is capable of establishing sound strategies for joint groundwater and surface water use. Six southeastern townships are found to be suitable locations for replacing groundwater with surface water from manmade lakes or reservoirs to meet drinking, irrigation, and aquaculture demands.

  12. Use of a ground-water flow model with particle tracking to evaluate ground-water vulnerability, Clark County, Washington

    USGS Publications Warehouse

    Snyder, Daniel T.; Wilkinson, James M.; Orzol, Leonard L.

    1998-01-01

    A ground-water flow model was used in conjunction with a particle-tracking program to demonstrate a method of evaluating ground-water vulnerability. The study area encompassed the part of the Portland Basin located in Clark County, Washington. A new computer program was developed that interfaces the particle-tracking program with a geographic information system (GIS). The GIS was used to display and analyze the particle-tracking results and to evaluate ground-water vulnerability by identifying recharge areas and their characteristics, determining the downgradient impact of land use at recharge areas, and estimating the age of ground water. The report presents a description of the methods used and the results of the evaluation of ground-water vulnerability.

  13. Wave-Induced Groundwater Flows in a Freshwater Beach Aquifer

    NASA Astrophysics Data System (ADS)

    Malott, S. S.; Robinson, C. E.; O'Carroll, D. M.

    2014-12-01

    Wave-induced recirculation across the sediment-water interface can impact the transport of pollutants through a beach aquifer and their ultimate flux into coastal waters. The fate of nutrients (e.g. from septic and agricultural sources) and fecal indicator bacteria (e.g. E. coil) near the sediment-water interface are of particular concern as these pollutants often lead to degradation of recreational water quality and nearshore ecosystems. This paper presents detailed field measurements of groundwater flows in a freshwater beach aquifer on Lake Huron over periods of intensified wave conditions. Quantifying wave-driven processes in a freshwater beach aquifer enables wave effects to be studied in isolation from density and tidal effects that complicate groundwater flows in marine beaches. Water exchange across the sediment-water interface and groundwater flow patterns were measured using groundwater wells, arrays of vertically nested pressure transducers and manometers. Results show that wave action induces rapid infiltration/exfiltration across the sediment-water interface and a larger recirculation cell through the beach aquifer. Field data is used to validate a numerical groundwater model of wave-induced groundwater flows. While prior studies have simulated the effects of waves on beach groundwater flows, this study is the first attempt to validate these sophisticated modeling approaches. Finally, field data illustrating the impact of wave-induced groundwater flows on nutrient and bacteria fate and transport in beach aquifers will also be presented.

  14. Using 14C and 3H to understand groundwater flow and recharge in an aquifer window

    NASA Astrophysics Data System (ADS)

    Atkinson, A. P.; Cartwright, I.; Gilfedder, B. S.; Cendón, D. I.; Unland, N. P.; Hofmann, H.

    2014-12-01

    Knowledge of groundwater residence times and recharge locations is vital to the sustainable management of groundwater resources. Here we investigate groundwater residence times and patterns of recharge in the Gellibrand Valley, southeast Australia, where outcropping aquifer sediments of the Eastern View Formation form an "aquifer window" that may receive diffuse recharge from rainfall and recharge from the Gellibrand River. To determine recharge patterns and groundwater flow paths, environmental isotopes (3H, 14C, δ13C, δ18O, δ2H) are used in conjunction with groundwater geochemistry and continuous monitoring of groundwater elevation and electrical conductivity. The water table fluctuates by 0.9 to 3.7 m annually, implying recharge rates of 90 and 372 mm yr-1. However, residence times of shallow (11 to 29 m) groundwater determined by 14C are between 100 and 10 000 years, 3H activities are negligible in most of the groundwater, and groundwater electrical conductivity remains constant over the period of study. Deeper groundwater with older 14C ages has lower δ18O values than younger, shallower groundwater, which is consistent with it being derived from greater altitudes. The combined geochemistry data indicate that local recharge from precipitation within the valley occurs through the aquifer window, however much of the groundwater in the Gellibrand Valley predominantly originates from the regional recharge zone, the Barongarook High. The Gellibrand Valley is a regional discharge zone with upward head gradients that limits local recharge to the upper 10 m of the aquifer. Additionally, the groundwater head gradients adjacent to the Gellibrand River are generally upwards, implying that it does not recharge the surrounding groundwater and has limited bank storage. 14C ages and Cl concentrations are well correlated and Cl concentrations may be used to provide a first-order estimate of groundwater residence times. Progressively lower chloride concentrations from 10

  15. Regional groundwater flow and geochemical evolution in the Amacuzac River Basin, Mexico

    NASA Astrophysics Data System (ADS)

    Morales-Casique, Eric; Guinzberg-Belmont, Jacobo; Ortega-Guerrero, Adrián

    2016-05-01

    An approach is presented to investigate the regional evolution of groundwater in the basin of the Amacuzac River in Central Mexico. The approach is based on groundwater flow cross-sectional modeling in combination with major ion chemistry and geochemical modeling, complemented with principal component and cluster analyses. The hydrogeologic units composing the basin, which combine aquifers and aquitards both in granular, fractured and karstic rocks, were represented in sections parallel to the regional groundwater flow. Steady-state cross-section numerical simulations aided in the conceptualization of the groundwater flow system through the basin and permitted estimation of bulk hydraulic conductivity values, recharge rates and residence times. Forty-five water locations (springs, groundwater wells and rivers) were sampled throughout the basin for chemical analysis of major ions. The modeled gravity-driven groundwater flow system satisfactorily reproduced field observations, whereas the main geochemical processes of groundwater in the basin are associated to the order and reactions in which the igneous and sedimentary rocks are encountered along the groundwater flow. Recharge water in the volcanic and volcano-sedimentary aquifers increases the concentration of HCO3 -, Mg2+ and Ca2+ from dissolution of plagioclase and olivine. Deeper groundwater flow encounters carbonate rocks, under closed CO2 conditions, and dissolves calcite and dolomite. When groundwater encounters gypsum lenses in the shallow Balsas Group or the deeper Huitzuco anhydrite, gypsum dissolution produces proportional increased concentration of Ca2+ and SO4 2-; two samples reflected the influence of hydrothermal fluids and probably halite dissolution. These geochemical trends are consistent with the principal component and cluster analyses.

  16. Groundwater dating for understanding nitrogen in groundwater systems - Time lag, fate, and detailed flow path ways

    NASA Astrophysics Data System (ADS)

    Morgenstern, Uwe; Hadfield, John; Stenger, Roland

    2014-05-01

    Nitrate contamination of groundwater is a problem world-wide. Nitrate from land use activities can leach out of the root zone of the crop into the deeper part of the unsaturated zone and ultimately contaminate the underlying groundwater resources. Nitrate travels with the groundwater and then discharges into surface water causing eutrophication of surface water bodies. To understand the source, fate, and future nitrogen loads to ground and surface water bodies, detailed knowledge of the groundwater flow dynamics is essential. Groundwater sampled at monitoring wells or discharges may not yet be in equilibrium with current land use intensity due to the time lag between leaching out of the root zone and arrival at the sampling location. Anoxic groundwater zones can act as nitrate sinks through microbial denitrification. However, the effect of denitrification on overall nitrate fluxes depends on the fraction of the groundwater flowing through such zones. We will show results from volcanic aquifers in the central North Island of New Zealand where age tracers clearly indicate that the groundwater discharges into large sensitive lakes like Lake Taupo and Lake Rotorua are not yet fully realising current land use intensity. The majority of the water discharging into these lakes is decades and up to over hundred years old. Therefore, increases in dairy farming over the last decades are not yet reflected in these old water discharges, but over time these increased nitrate inputs will eventually work their way through the large groundwater systems and increasing N loads to the lakes are to be expected. Anoxic zones are present in some of these aquifers, indicating some denitrification potential, however, age tracer results from nested piezo wells show young groundwater in oxic zones indicating active flow in these zones, while anoxic zones tend to have older water indicating poorer hydraulic conductivity in these zones. Consequently, to evaluate the effect of denitrification

  17. Groundwater balance estimation in karst by using simple conceptual rainfall-runoff model

    NASA Astrophysics Data System (ADS)

    Željković, Ivana; Kadić, Ana; Denić-Jukić, Vesna

    2014-05-01

    The objective of this work is the study of Opačac karst spring which geographically lies in Dalmatia (Croatia). Numerous studies have been carried out in karst aiming the investigation of groundwater regime. The karst spring hydrograph can reflect the groundwater regime and consequently the analysis is based on them. A simple conceptual rainfall-runoff model is proposed for the estimation of groundwater balance components including the influences of time invariant catchment boundaries and intercatchment flows. The proposed parameter estimation procedure merges the soil-moisture balance and the groundwater balance approaches to obtain the complete groundwater budget. The effective rainfall is calculated by using mathematical model based on soil-moisture balance equations i.e. Palmer's fluid mass balance method. The parameters of model of effective rainfall are determined by using simple conceptual rainfall-runoff model consisting of two linear reservoirs representing the fast and slow flow component of the recession. The weight coefficient between the fast and slow component is determined by using BFI (Base Flow Index) analysis of hydrograph. Recession coefficient of the slow flow component and the weight coefficient are determined from hydrograph analysis. Available data from nearby meteorological station includes on daily basis daily average discharge, the amount of precipitation, the average temperature and the humidity from 1995-2010. The average catchment area is also estimated with the average yearly runoff deficit using Turc's method and compared with the values obtained from the application of the rainfall-runoff model. Nash-Sutcliffe model efficiency coefficient for simulated hydrograph is applied to assess the predictive power of model. Calculated groundwater balance shows that the Opačac Spring aquifer contains a significant storage capacity. The application of series of linear reservoirs is a classical and common technique, but the proposed simple

  18. Comparison of different estimation techniques to quantify groundwater recharge in Pirna, Germany

    NASA Astrophysics Data System (ADS)

    Ringleb, Jana; Sallwey, Jana; Stefan, Catalin

    2015-04-01

    Water scarcity in combination with groundwater exploitation is a major concern worldwide because of climate change, population growth and rising water demand. To be able to sustainably manage and protect groundwater resources, it is necessary to quantify the amount of water which leaks through the unsaturated zone and recharges the aquifer naturally. However, quantifying the spatial and temporal distribution of recharge is difficult because of soil heterogeneity and the influence of vegetation. For that reason and because field measurements of recharge are difficult to obtain, models are valuable tools to quantify recharge. Numerical models need a lot of parameters which are hard to measure and hence can only be estimated. Therefore analytical models or empirical equations which use less and / or easier obtainable parameters could estimate groundwater recharge as well as numerical models because of the underlying uncertainty in parameter estimation. Recharge estimation methods which use different model approaches and have varying complexity were compared at Pirna test field site, Germany to select suitable methods which will later be integrated into a web-based Decision Support System (DSS) developed for the sustainable management of groundwater. The complexity of the used methods covers numerical models, analytical models as well as empirical equations. Different model approaches were used to estimate groundwater recharge including amongst others a groundwater flow model, an unsaturated zone model and a watershed model. The resulting groundwater recharge estimates received from the numerical and analytical models and from empirical equations were compared to evaluate whether the methods are suitable to estimate groundwater recharge considering the complexity, data requirements and time-consumption of each method.

  19. Estimation of ecological high flow

    NASA Astrophysics Data System (ADS)

    Lin, Jen-Yang; Chen, Yen-Chang; Hsienshao Tsao, Eric

    2006-02-01

    Floods can destroy fish habitat. During a flood a fish has to seek shelters (refuges) to survive. It is necessary to know the maximum discharge that the fish can sustain against the strong current. Ecological and hydraulic engineers can simulate the flow condition of high flow for designing the refuge when restoring and enhancing the rivers are needed. Based on the average ratio of the mean and maximum velocities invariant with time, discharge and water level, this paper tries to introduce the concept of ecological high flow. The mean-maximum velocity ratio can be used to estimate the mean velocity of the river. If the maximum velocity of the cross section is replaced by the maximum sustained swimming speeds of fish, the mean velocity of ecological high flow can be calculated with the constant ratio. The cross-sectional area can be estimated by the gage height. Then the ecological high flow can be estimated as the product of mean velocity of ecological high flow multiplied by the cross-sectional area. The available data of the upstream of the Dacha River where is the habitat of the Formosan landlocked salmon were used to illustrate the estimation of the ecological high flow. Any restoration project at Sonmou that try to improve the stream habitat can use the ecological high flow to design the hydraulic structure at suitable location to offer refuges for the Formosan landlocked salmon that is an endangered species in Taiwan

  20. Intelligent Flow Friction Estimation

    PubMed Central

    Brkić, Dejan; Ćojbašić, Žarko

    2016-01-01

    Nowadays, the Colebrook equation is used as a mostly accepted relation for the calculation of fluid flow friction factor. However, the Colebrook equation is implicit with respect to the friction factor (λ). In the present study, a noniterative approach using Artificial Neural Network (ANN) was developed to calculate the friction factor. To configure the ANN model, the input parameters of the Reynolds Number (Re) and the relative roughness of pipe (ε/D) were transformed to logarithmic scales. The 90,000 sets of data were fed to the ANN model involving three layers: input, hidden, and output layers with, 2, 50, and 1 neurons, respectively. This configuration was capable of predicting the values of friction factor in the Colebrook equation for any given values of the Reynolds number (Re) and the relative roughness (ε/D) ranging between 5000 and 108 and between 10−7 and 0.1, respectively. The proposed ANN demonstrates the relative error up to 0.07% which had the high accuracy compared with the vast majority of the precise explicit approximations of the Colebrook equation. PMID:27127498

  1. Intelligent Flow Friction Estimation.

    PubMed

    Brkić, Dejan; Ćojbašić, Žarko

    2016-01-01

    Nowadays, the Colebrook equation is used as a mostly accepted relation for the calculation of fluid flow friction factor. However, the Colebrook equation is implicit with respect to the friction factor (λ). In the present study, a noniterative approach using Artificial Neural Network (ANN) was developed to calculate the friction factor. To configure the ANN model, the input parameters of the Reynolds Number (Re) and the relative roughness of pipe (ε/D) were transformed to logarithmic scales. The 90,000 sets of data were fed to the ANN model involving three layers: input, hidden, and output layers with, 2, 50, and 1 neurons, respectively. This configuration was capable of predicting the values of friction factor in the Colebrook equation for any given values of the Reynolds number (Re) and the relative roughness (ε/D) ranging between 5000 and 10(8) and between 10(-7) and 0.1, respectively. The proposed ANN demonstrates the relative error up to 0.07% which had the high accuracy compared with the vast majority of the precise explicit approximations of the Colebrook equation. PMID:27127498

  2. Uncertainty in global groundwater storage estimates in a Total Groundwater Stress framework

    PubMed Central

    Richey, Alexandra S.; Thomas, Brian F.; Lo, Min‐Hui; Swenson, Sean; Rodell, Matthew

    2015-01-01

    Abstract Groundwater is a finite resource under continuous external pressures. Current unsustainable groundwater use threatens the resilience of aquifer systems and their ability to provide a long‐term water source. Groundwater storage is considered to be a factor of groundwater resilience, although the extent to which resilience can be maintained has yet to be explored in depth. In this study, we assess the limit of groundwater resilience in the world's largest groundwater systems with remote sensing observations. The Total Groundwater Stress (TGS) ratio, defined as the ratio of total storage to the groundwater depletion rate, is used to explore the timescales to depletion in the world's largest aquifer systems and associated groundwater buffer capacity. We find that the current state of knowledge of large‐scale groundwater storage has uncertainty ranges across orders of magnitude that severely limit the characterization of resilience in the study aquifers. Additionally, we show that groundwater availability, traditionally defined as recharge and redefined in this study as total storage, can alter the systems that are considered to be stressed versus unstressed. We find that remote sensing observations from NASA's Gravity Recovery and Climate Experiment can assist in providing such information at the scale of a whole aquifer. For example, we demonstrate that a groundwater depletion rate in the Northwest Sahara Aquifer System of 2.69 ± 0.8 km3/yr would result in the aquifer being depleted to 90% of its total storage in as few as 50 years given an initial storage estimate of 70 km3. PMID:26900184

  3. New method to estimate variability in groundwater volumes

    NASA Astrophysics Data System (ADS)

    Wendel, JoAnna

    2014-04-01

    As the warming climate and increasing population put stress on the world's water supply, it has become increasingly important to have a global understanding of how groundwater volumes vary from season to season and from year to year. Current global hydrological models do not include lateral groundwater flow, which plays a significant role in providing water to plants and in recharging lakes, rivers, and streams.

  4. Comparison of wetland evapotranspiration estimates using diurnal groundwater fluctuations and measurements of a groundwater lysimeter

    NASA Astrophysics Data System (ADS)

    Fahle, Marcus; Dietrich, Ottfried; Lischeid, Gunnar

    2013-04-01

    Sound water management in wetlands requires knowledge of on-going processes and estimates of the water balance components. Specifically evapotranspiration is of crucial importance, as it is often the main water extracting quantity. To avoid elaborate and expensive equipment, which is often required for estimating actual values, potential evapotranspiration is frequently used, which can be easily derived from standard meteorological measurements. However, the potential values may under- or overestimate actual evapotranspiration significantly. A cheap and easy-to-use method for estimating actual values in shallow groundwater environments relies on diurnal groundwater fluctuation. Basically the 24 hours groundwater level decline, considering in some way the prevalent groundwater recovery, is multiplied by the readily available specific yield. Various varieties of this approach have been employed for that purpose, above all differing in their assumptions on groundwater recovery, i.e. lateral or vertical in- or outflow. The objective of our study is therefore to compare these different methods. For this purpose we use data of a weighable groundwater lysimeter situated at a ditch drained grassland site in the Spreewald wetland in Northeastern Germany. The groundwater level in the lysimeter was adjusted to a reference gauge and simulated the conditions of the surrounding area. Hence the lysimeter reflected near natural conditions and provided measurements of all water balance components with high temporal resolution (up to 10 minute intervals). Suitable days, i.e. with a pronounced diurnal fluctuation, of the vegetation periods 2011 and 2012 are chosen and used to prove common assumptions about groundwater recharge, e.g. if the values remain constant during the day or if diurnal variations resulting from gradient changes exist. Finally, based on the lysimeter measurements, the evapotranspiration estimates gained from different approaches that employ diurnal groundwater

  5. Wellbore and groundwater temperature distribution eastern Snake River Plain, Idaho: Implications for groundwater flow and geothermal potential

    NASA Astrophysics Data System (ADS)

    McLing, Travis L.; Smith, Richard P.; Smith, Robert W.; Blackwell, David D.; Roback, Robert C.; Sondrup, Andrus J.

    2016-06-01

    A map of groundwater temperatures from the Eastern Snake River Plain (ESRP) regional aquifer can be used to identify and interpret important features of the aquifer, including aquifer flow direction, aquifer thickness, and potential geothermal anomalies. The ESRP is an area of high heat flow, yet most of this thermal energy fails to reach the surface, due to the heat being swept downgradient by the aquifer to the major spring complexes near Thousand Springs, ID, a distance of 300 km. Nine deep boreholes that fully penetrate the regional aquifer display three common features: (1) high thermal gradients beneath the aquifer, corresponding to high conductive heat flow in low-permeability hydrothermally-altered rocks; (2) isothermal temperature profiles within the aquifer, characteristic of an actively flowing groundwater; and (3) moderate thermal gradients in the vadose zone with values that indicate that over half of the geothermal heat flow is removed by advective transport in the regional aquifer system. This study utilized temperature data from 250 ESRP aquifer wells to evaluate regional aquifer flow direction, aquifer thickness, and potential geothermal anomalies. Because the thermal gradients are typically low in the aquifer, any measurement of groundwater temperature is a reasonable estimate of temperature throughout the aquifer thickness, allowing the construction of a regional aquifer temperature map for the ESRP. Mapped temperatures are used to identify cold thermal plumes associated with recharge from tributary valleys and adjacent uplands, and warm zones associated with geothermal input to the aquifer. Warm zones in the aquifer can have various causes, including local circulation of groundwater through the deep conductively dominated region, slow groundwater movement in low-permeability regions, or localized heat flow from deeper thermal features.

  6. Present-day groundwater recharge estimation in parts of the Indian Sub-Continent

    NASA Astrophysics Data System (ADS)

    Bhanja, S. N.; Mukherjee, A.; Wada, Y.; Scanlon, B. R.; Taylor, R. G.; Rodell, M.; Malakar, P.

    2015-12-01

    Large part of global population has been dependent on groundwater as a source of fresh water. The demand would further increase with increasing population and stress associated with climate change. We tried to provide regional-scale groundwater recharge estimates in a large part of Indian Sub-Continent. A combination of ground-based, satellite-based and numerical model simulated recharge estimates were presented in the densely populated region. Three different methods: an intense network of observational wells (n>13,000 wells), a satellite (TRMM) and global land-surface model (CLM) outputs, and a global-scale hydrological model (PCR GLOBWB) were employed to calculate recharge estimates. Groundwater recharge values exhibit large spatial variations over the entire region on the basis of aquifer hydrogeology, precipitation and groundwater withdrawal patterns. Groundwater recharge estimates from all three estimation techniques were found to be higher (>300 mm/year) in fertile planes of Indus-Ganges-Brahmaputra (IGB) river basins. A combination of favorable hydrogeologic conditions (porosity, permeability etc.), comparatively higher rates of precipitation, and return flow from rapidly withdrawn irrigation water might influence occurrence of high recharge rates. However, central and southern study area experiences lower recharge rates (<200 mm/year), might be associated with unfavorable hydrogeologic conditions associated with cratonic provinces. Statistical analysis of inter-comparison between the three different recharge estimates show good matches in some of the areas. Recharge estimates indicate dynamic nature of groundwater recharge as a function of precipitation, land use pattern, and hydrogeologic parameters. On a first hand basis, the estimates will help policy makers to understand groundwater recharge process over the densely populated region and finally would facilitate to implement sustainable policy for securing water security.

  7. Regional Groundwater Processes and Flow Dynamics from Age Tracer Data

    NASA Astrophysics Data System (ADS)

    Morgenstern, Uwe; Stewart, Mike K.; Matthews, Abby

    2016-04-01

    Age tracers are now used in New Zealand on regional scales for quantifying the impact and lag time of land use and climate change on the quantity and quality of available groundwater resources within the framework of the National Policy Statement for Freshwater Management 2014. Age tracers provide measurable information on the dynamics of groundwater systems and reaction rates (e.g. denitrification), essential for conceptualising the regional groundwater - surface water system and informing the development of land use and groundwater flow and transport models. In the Horizons Region of New Zealand, around 200 wells have tracer data available, including tritium, SF6, CFCs, 2H, 18O, Ar, N2, CH4 and radon. Well depths range from shallower wells in gravel aquifers in the Horowhenua and Tararua districts, and deeper wells in the aquifers between Palmerston North and Wanganui. Most of the groundwater samples around and north of the Manawatu River west of the Tararua ranges are extremely old (>100 years), even from relatively shallow wells, indicating that these groundwaters are relatively disconnected from fresh surface recharge. The groundwater wells in the Horowhenua tap into a considerably younger groundwater reservoir with groundwater mean residence time (MRT) of 10 - 40 years. Groundwater along the eastern side of the Tararua and Ruahine ranges is significantly younger, typically <5 years MRT. Vertical groundwater recharge rates, as deduced from groundwater depth and MRT, are extremely low in the central coastal area, consistent with confined groundwater systems, or with upwelling of old groundwater close to the coast. Very low vertical recharge rates along the Manawatu River west of the Manawatu Gorge indicate upwelling groundwater conditions in this area, implying groundwater discharge into the river is more likely here than loss of river water into the groundwater system. High recharge rates observed at several wells in the Horowhenua area and in the area east of

  8. Locating groundwater flow in karst by acoustic emission surveys

    SciTech Connect

    Stokowski, S.J. Jr.; Clark, D.A.

    1985-01-01

    An acoustic emission survey of Newala Fm. (primarily dolomite) karst has helped to locate subsurface water flow. This survey was performed on the Rock Quarry Dome, Sevier County, Tennessee. A Dresser RS-4 recording seismograph, adjusted to provide a gain of 1000, collected acoustic emission data using Mark Products CN368 vertical geophones with 3-inch spikes. Data was collected for 5-15 second intervals. The geophones were laid out along traverses with 10, 20, or 30-ft spacing and covered with sand bags in locations of high ambient noise. Traverses were laid out: along and across lineaments known to correspond with groundwater flow in natural subsurface channels; across and along a joint-controlled sink suspected of directing groundwater flow; and across a shallow sinkhole located tangentially to the Little Pigeon River and suspected of capturing river water for the groundwater system. Acoustic emissions of channelized flowing groundwater have a characteristic erratic spiked spectral signature. These acoustic emission signatures increase in amplitude and number in the immediate vicinity of the vertical projection of channelized groundwater flow if it occurs within approximately 30 feet of the surface. If the groundwater flow occurs at greater depths the emissions may be offset from the projection of the actual flow, due to propagation of the signal along rock pinnacles or attenuation by residual soils.

  9. The water balance estimation for catastrophic floods: groundwater contribution

    NASA Astrophysics Data System (ADS)

    Arakelian, Sergey; Vinogradov, Alexey; Tulenev, Nikita; Trifonova, Tatiana

    2016-04-01

    1. We discuss the existing problems in the study of the mechanisms of formation of catastrophic floods taking into account the possible influence of groundwater. The difficulty in assessing the causes of the disastrous floods is linked to the lack of direct field measurements of precipitation and so, to estimate the water balance in the rain floods. The problems that arise when comparing the results of observations and measurements of rainfall floods are considered. 2. We rely on the concept, where groundwater and surface water are the two coupled factors resulting in catastrophic floods/debris, and they are not isolated systems. These two units are closely related to each other on the territory of a unified watershed under its functioning including the overall transport system, i.e. 3D-network of cracks in the rock (visible manifestation on the land surface of which is the rivershed itself). 3. We estimated the pressure in the aquifer taking the data obtained by the observable mudflow or flood as a base. According to our calculations in the case of a violent release, such pressure for the really observed events can reach tens of atmospheres. Such pressure enhancement may occur due to various external factors (including the nature climatic and seismic processes). 4. A more detailed analysis should be carried out in accordance with a real topology of multiple cracks taking into account the non-stationary process and levels of resistance for water flows in different sections of crack-net (hydrostatic/hydrodynamic pressures in underground aquifers).

  10. Simulation of ground-water flow and land subsidence in the Antelope Valley ground-water basin, California

    USGS Publications Warehouse

    Leighton, David A.; Phillips, Steven P.

    2003-01-01

    ground-water development have eliminated the natural sources of discharge, and pumping for agricultural and urban uses have become the primary source of discharge from the ground-water system. Infiltration of return flows from agricultural irrigation has become an important source of recharge to the aquifer system. The ground-water flow model of the basin was discretized horizontally into a grid of 43 rows and 60 columns of square cells 1 mile on a side, and vertically into three layers representing the upper, middle, and lower aquifers. Faults that were thought to act as horizontal-flow barriers were simulated in the model. The model was calibrated to simulate steady-state conditions, represented by 1915 water levels and transient-state conditions during 1915-95 using water-level and subsidence data. Initial estimates of the aquifer-system properties and stresses were obtained from a previously published numerical model of the Antelope Valley ground-water basin; estimates also were obtained from recently collected hydrologic data and from results of simulations of ground-water flow and land subsidence models of the Edwards Air Force Base area. Some of these initial estimates were modified during model calibration. Ground-water pumpage for agriculture was estimated on the basis of irrigated crop acreage and crop consumptive-use data. Pumpage for public supply, which is metered, was compiled and entered into a database used for this study. Estimated annual pumpage peaked at 395,000 acre-feet (acre-ft) in 1952 and then declined because of declining agricultural production. Recharge from irrigation-return flows was estimated to be 30 percent of agricultural pumpage; the irrigation-return flows were simulated as recharge to the regional water table 10 years following application at land surface. The annual quantity of natural recharge initially was based on estimates from previous studies. During model calibration, natural recharge was reduced from the initial

  11. MODFLOW-2000, The U.S. Geological Survey Modular Ground-Water Model - User Guide to Modularization Concepts and the Ground-Water Flow Process

    USGS Publications Warehouse

    Harbaugh, Arlen W.; Banta, Edward R.; Hill, Mary C.; McDonald, Michael G.

    2000-01-01

    MODFLOW is a computer program that numerically solves the three-dimensional ground-water flow equation for a porous medium by using a finite-difference method. Although MODFLOW was designed to be easily enhanced, the design was oriented toward additions to the ground-water flow equation. Frequently there is a need to solve additional equations; for example, transport equations and equations for estimating parameter values that produce the closest match between model-calculated heads and flows and measured values. This report documents a new version of MODFLOW, called MODFLOW-2000, which is designed to accommodate the solution of equations in addition to the ground-water flow equation. This report is a user's manual. It contains an overview of the old and added design concepts, documents one new package, and contains input instructions for using the model to solve the ground-water flow equation.

  12. An evaluation of GRACE groundwater estimates over East Africa

    NASA Astrophysics Data System (ADS)

    Nanteza, J.; Thomas, B. F.; de Linage, C.; Famiglietti, J. S.

    2013-12-01

    The East African (EA) region, comprised of five countries (Uganda, Kenya, Tanzania, Rwanda and Burundi), is among those regions characterized as vulnerable to water stress. The region's freshwater resources, both surface and groundwater, are impacted due to increased pressure from changes in climate and human activities. Better management approaches are required to ensure that these pressures do not significantly impact water availability and accessibility. However, the lack of adequate ground-based observation networks to monitor freshwater resources - especially groundwater (the major source of freshwater in EA), limits effective management of the available water resources. In this study, we explore the potential of using remotely sensed data to monitor freshwater resources over EA. The study uses data from the Gravity Recovery and Climate Experiment (GRACE) satellite to estimate groundwater storage variations over EA during the last decade. The satellite's performance in accurately observing changes in groundwater storage is examined by evaluating the GRACE groundwater estimates against spatially interpolated in-situ groundwater observations using goodness of fit criteria including linear regression coefficient, coefficient of determination and root mean square errors. The results demonstrate that GRACE performs well in observing the behavior of groundwater storage. These results can be useful in improving land surface model simulations - a basis for better decision making in water resources management in the region.

  13. Permafrost thaw in a nested groundwater-flow system

    USGS Publications Warehouse

    McKenzie, Jeffery M.; Voss, Clifford I.

    2013-01-01

    Groundwater flow in cold regions containing permafrost accelerates climate-warming-driven thaw and changes thaw patterns. Simulation analyses of groundwater flow and heat transport with freeze/thaw in typical cold-regions terrain with nested flow indicate that early thaw rate is particularly enhanced by flow, the time when adverse environmental impacts of climate-warming-induced permafrost loss may be severest. For the slowest climate-warming rate predicted by the Intergovernmental Panel on Climate Change (IPCC), once significant groundwater flow begins, thick permafrost layers can vanish in several hundred years, but survive over 1,000 years where flow is minimal. Large-scale thaw depends mostly on the balance of heat advection and conduction in the supra-permafrost zone. Surface-water bodies underlain by open taliks allow slow sub-permafrost flow, with lesser influence on regional thaw. Advection dominance over conduction depends on permeability and topography. Groundwater flow around permafrost and flow through permafrost impact thaw differently; the latter enhances early thaw rate. Air-temperature seasonality also increases early thaw. Hydrogeologic heterogeneity and topography strongly affect thaw rates/patterns. Permafrost controls the groundwater/surface-water-geomorphology system; hence, prediction and mitigation of impacts of thaw on ecology, chemical exports and infrastructure require improved hydrogeology/permafrost characterization and understanding

  14. Influence of vertical flows in wells on groundwater sampling.

    PubMed

    McMillan, Lindsay A; Rivett, Michael O; Tellam, John H; Dumble, Peter; Sharp, Helen

    2014-11-15

    Pumped groundwater sampling evaluations often assume that horizontal head gradients predominate and the sample comprises an average of water quality variation over the well screen interval weighted towards contributing zones of higher hydraulic conductivity (a permeability-weighted sample). However, the pumping rate used during sampling may not always be sufficient to overcome vertical flows in wells driven by ambient vertical head gradients. Such flows are reported in wells with screens between 3 and 10m in length where lower pumping rates are more likely to be used during sampling. Here, numerical flow and particle transport modeling is used to provide insight into the origin of samples under ambient vertical head gradients and under a range of pumping rates. When vertical gradients are present, sample provenance is sensitive to pump intake position, pumping rate and pumping duration. The sample may not be drawn from the whole screen interval even with extended pumping times. Sample bias is present even when the ambient vertical flow in the wellbore is less than the pumping rate. Knowledge of the maximum ambient vertical flow in the well does, however, allow estimation of the pumping rate that will yield a permeability-weighted sample. This rate may be much greater than that recommended for low-flow sampling. In practice at monitored sites, the sampling bias introduced by ambient vertical flows in wells may often be unrecognized or underestimated when drawing conclusions from sampling results. It follows that care should be taken in the interpretation of sampling data if supporting flow investigations have not been undertaken. PMID:24999176

  15. PUMa - modelling the groundwater flow in Baltic Sedimentary Basin

    NASA Astrophysics Data System (ADS)

    Kalvane, G.; Marnica, A.; Bethers, U.

    2012-04-01

    In 2009-2012 at University of Latvia and Latvia University of Agriculture project "Establishment of interdisciplinary scientist group and modelling system for groundwater research" is implemented financed by the European Social Fund. The aim of the project is to develop groundwater research in Latvia by establishing interdisciplinary research group and modelling system covering groundwater flow in the Baltic Sedimentary Basin. Researchers from fields like geology, chemistry, mathematical modelling, physics and environmental engineering are involved in the project. The modelling system is used as a platform for addressing scientific problems such as: (1) large-scale groundwater flow in Baltic Sedimentary Basin and impact of human activities on it; (2) the evolution of groundwater flow since the last glaciation and subglacial groundwater recharge; (3) the effects of climate changes on shallow groundwater and interaction of hydrographical network and groundwater; (4) new programming approaches for groundwater modelling. Within the frame of the project most accessible geological information such as description of geological wells, geological maps and results of seismic profiling in Latvia as well as Estonia and Lithuania are collected and integrated into modelling system. For example data form more then 40 thousands wells are directly used to automatically generate the geological structure of the model. Additionally a groundwater sampling campaign is undertaken. Contents of CFC, stabile isotopes of O and H and radiocarbon are the most significant parameters of groundwater that are established in unprecedented scale for Latvia. The most important modelling results will be published in web as a data set. Project number: 2009/0212/1DP/1.1.1.2.0/09/APIA/VIAA/060. Project web-site: www.puma.lu.lv

  16. The effects of groundwater abstraction on low flows

    NASA Astrophysics Data System (ADS)

    de Graaf, I. E. M.; van Beek, L. P. H.; Wada, Y.; Bierkens, M. F. P.

    2012-04-01

    In regions with frequent water stress and large aquifer systems, groundwater often constitutes an essential source of water. If groundwater abstraction exceeds groundwater recharge over a long time and over large areas persistent groundwater depletion can occur. The resulting lowering of groundwater levels can have negative effects on agricultural productivity but also on natural streamflow and associated wetlands and ecosystems, in particular during low-flow events when the groundwater contribution through baseflow is relatively large. In this study we focus on the effects of global groundwater abstraction on low-flow magnitude, frequency and duration for the major rivers of the world for the period 1960-2000. As a basis, we use the large-scale hydrological model PCR-GLOBWB that calculates all major water balance terms on a daily time step at a 0.5ox0.5o resolution. Currently, PCR-GLOBWB represents groundwater and the associated baseflow by means of a linear reservoir that is parameterized using global lithological data and drainage density. It simulates renewable groundwater storage within each 0.5o cell. Lateral flow between cells is not considered. The specific runoff from the model is subsequently transformed into discharge by means of a kinematic wave routing scheme. In this study we perform a sensitivity analysis in which we evaluate the effects of total water demand for the period 1960-2000 (Wada et al., 2011: doi:10.5194/hess-15-3785-2011). This demand is preferentially met by renewable groundwater storage, secondly by surface water. Any remainder is assumed to stem from non- renewable (i.e. fossil) groundwater resources. Thus, groundwater abstractions act as a direct sink of (renewable) groundwater storage, whereas surface water abstractions act as a direct sink of streamflow. The resulting response is non-trivial as abstractions are variably taken from both groundwater and surface water, where return-flows contribute to a single source: return flow from

  17. Estimates of vertical hydraulic conductivity and regional ground-water flow rates in rocks of Jurassic and Cretaceous age, San Juan Basin, New Mexico and Colorado

    USGS Publications Warehouse

    Frenzel, P.F.; Lyford, F.P.

    1982-01-01

    The San Juan structural basin northwestern New Mexico was modeled in three dimensions using a finite-difference, steady-state model. The modeled space was divided into seven layers of square prisms that were 6 miles on a side in the horizontal directions. In the vertical direction, the layers of prisms ranged in thickness from 300 to 1,500 feet. The model included the geologic section between the base of the Entrada Sandstone and the top of Mesaverde Group. Principal aquifers in this section are mostly confined and include the Entrada Sandstone, the Westwater Canyon Member of the Morrison Formation , and the Gallup Sandstone. Values for vertical hydraulic conductivities from 10 to the minus 12th power to 10 to the minus 11th power feet per second for the confining layers gave a good simulation of head differences between layers, but a sensitivity analysis indicated that these values could be between 10 and 100 times greater. The model-derived steady-state flow was about 30 cubic feet per second. About one-half of the flow was in the San Juan River drainage basin about one-third in the Rio Grande drainage basin, and one-sixth in the Puerco River drainage basin. (USGS)

  18. Improved Uncertainty Quantification in Groundwater Flux Estimation Using GRACE

    NASA Astrophysics Data System (ADS)

    Reager, J. T., II; Rao, P.; Famiglietti, J. S.; Turmon, M.

    2015-12-01

    Groundwater change is difficult to monitor over large scales. One of the most successful approaches is in the remote sensing of time-variable gravity using NASA Gravity Recovery and Climate Experiment (GRACE) mission data, and successful case studies have created the opportunity to move towards a global groundwater monitoring framework for the world's largest aquifers. To achieve these estimates, several approximations are applied, including those in GRACE processing corrections, the formulation of the formal GRACE errors, destriping and signal recovery, and the numerical model estimation of snow water, surface water and soil moisture storage states used to isolate a groundwater component. A major weakness in these approaches is inconsistency: different studies have used different sources of primary and ancillary data, and may achieve different results based on alternative choices in these approximations. In this study, we present two cases of groundwater change estimation in California and the Colorado River basin, selected for their good data availability and varied climates. We achieve a robust numerical estimate of post-processing uncertainties resulting from land-surface model structural shortcomings and model resolution errors. Groundwater variations should demonstrate less variability than the overlying soil moisture state does, as groundwater has a longer memory of past events due to buffering by infiltration and drainage rate limits. We apply a model ensemble approach in a Bayesian framework constrained by the assumption of decreasing signal variability with depth in the soil column. We also discuss time variable errors vs. time constant errors, across-scale errors v. across-model errors, and error spectral content (across scales and across model). More robust uncertainty quantification for GRACE-based groundwater estimates would take all of these issues into account, allowing for more fair use in management applications and for better integration of GRACE

  19. Quantitative dye-tracing of karst ground-water flow

    USGS Publications Warehouse

    Smoot, James; Mull, Donald; Liebermann, Timothy

    1989-01-01

    Analysis of the results of repeat quantitative dye traces between a sinkhole and a spring used for public water supply were used to describe predictive relations between discharge, mean travel time, apparent ground-water flow velocity and solute transport characteristics. Normalized peak concentration, mean travel time, and standard deviation of travel times were used to produce a dimensionless, composite type curve that was used to produce a dimensionless, composite type curve that was used to simulate solute transport characteristics for selected discharges. Using this curve and previously developed statistical relations, a water manager can estimate the arrival time, peak concentration, and persistence of a soluble contaminant at a supply spring or well based on discharge and the quantity of spilled contaminant.

  20. Groundwater flow and solute movement to drain laterals, western San Joaquin Valley, California: 1. Geochemical assessment

    USGS Publications Warehouse

    Deverel, S.J.; Fio, John L.

    1991-01-01

    A study was undertaken to quantitatively evaluate the hydrologic processes affecting the chemical and isotopic composition of drain lateral water in a drained agricultural field in the western San Joaquin Valley, California. The results elucidate the process of mixing of deep and shallow groundwater (below and within 6 m from land surface) entering the drain laterals. The deep groundwater was subject to evapoconcentration prior to drainage system installation and has been displaced downward (to depths greater than 6 m) in the groundwater system. The proportions of deep and shallow groundwater entering the drain laterals was calculated from the end-member oxygen 18 compositions determined in groundwater samples. The percentage of total drain lateral flow which is deep groundwater flow is about 30% for the shallow drain lateral (1.8 m below land surface) (drain lateral 1)) and 60% for the deep drain lateral (2.7 m below land surface (drain lateral 2)). During irrigation, the percentages of deep groundwater flow decrease to 0 and 30% for the shallow and deep drain laterals, respectively. Selenium concentrations in drain lateral waters decrease during irrigation but selenium loads increase. Total estimated annual loads were 1.1 and 5.4 kg of selenium for drain laterals 1 and 2, respectively. Substantial percentages of the annual load occurred during 8 days of irrigation, 23 and 9% for drain laterals 1 and 2, respectively.

  1. Accounting for Transport Parameter Uncertainty in Geostatistical Groundwater Contaminant Release History Estimation

    NASA Astrophysics Data System (ADS)

    Ostrowski, J.; Shlomi, S.; Michalak, A.

    2007-12-01

    The process of estimating the release history of a contaminant in groundwater relies on coupling a limited number of concentration measurements with a groundwater flow and transport model in an inverse modeling framework. The information provided by available measurements is generally not sufficient to fully characterize the unknown release history; therefore, an accurate assessment of the estimation uncertainty is required. The modeler's level of confidence in the transport parameters, expressed as pdfs, can be incorporated into the inverse model to improve the accuracy of the release estimates. In this work, geostatistical inverse modeling is used in conjunction with Monte Carlo sampling of transport parameters to estimate groundwater contaminant release histories. Concentration non-negativity is enforced using a Gibbs sampling algorithm based on a truncated normal distribution. The method is applied to two one-dimensional test cases: a hypothetical dataset commonly used in validating contaminant source identification methods, and data collected from a tetrachloroethylene and trichloroethylene plume at the Dover Air Force Base in Delaware. The estimated release histories and associated uncertainties are compared to results from a geostatistical inverse model where uncertainty in transport parameters is ignored. Results show that the a posteriori uncertainty associated with the model that accounts for parameter uncertainty is higher, but that this model provides a more realistic representation of the release history based on available data. This modified inverse modeling technique has many applications, including assignment of liability in groundwater contamination cases, characterization of groundwater contamination, and model calibration.

  2. Drain discharge monitoring to estimate plot scale groundwater recharge

    NASA Astrophysics Data System (ADS)

    Di Ciacca, Antoine; Leterme, Bertrand; Jacques, Diederik; Vanderborght, Jan

    2016-04-01

    Spatially and temporally distributed representation of groundwater recharge is known as an important issue in hydrogeological modelling. Therefore, monitored groundwater recharge data are crucial to parameterize and/or validate groundwater flow models. Generally, river base flow measurements are used for this purpose with the assumption that these catchment-scale spatially aggregated measurements are suitable to assess the internal catchment behaviour. However, the signal of different soil and vegetation types is lost and this limits our ability to validate mechanistic, process-based models used at the plot scale. A suitable alternative in lowland drained areas could lie in monitoring drain discharge. The present poster describes the set-up of such a monitoring device in a ditch drain located in the Kleine Nete catchment (Belgium). To calculate groundwater recharge rate from drain discharge monitoring, some intermediate steps are required. The contributing area has to be delimited and the contribution of recharge water from this area has to be separated from other possible sources (e.g. deeper groundwater flow, run-off). To handle this, some assumptions regarding the features of the plot and some additional measurements have been used.

  3. Stable isotope and groundwater flow dynamics of agricultural irrigation recharge into groundwater resources of the Central Valley, California

    SciTech Connect

    Davisson, M.L.; Criss, R.E.

    1995-01-01

    Intensive agricultural irrigation and overdraft of groundwater in the Central Valley of California profoundly affect the regional quality and availability of shallow groundwater resources. In the natural state, the {delta}{sup 18}O values of groundwater were relatively homogeneous (mostly -7.0 {+-} 0.5{per_thousand}), reflecting local meteoric recharge that slowly (1-3m/yr) flowed toward the valley axis. Today, on the west side of the valley, the isotope distribution is dominated by high {sup 18}O enclosures formed by recharge of evaporated irrigation waters, while the east side has bands of low {sup 18}O groundwater indicating induced recharge from rivers draining the Sierra Nevada mountains. Changes in {delta}{sup 18}O values caused by the agricultural recharge strongly correlate with elevated nitrate concentrations (5 to >100 mg/L) that form pervasive, non-point source pollutants. Small, west-side cities dependent solely on groundwater resources have experienced increases of >1.0 mg/L per year of nitrate for 10-30 years. The resultant high nitrates threaten the economical use of the groundwater for domestic purposes, and have forced some well shut-downs. Furthermore, since >80% of modern recharge is now derived from agricultural irrigation, and because modern recharge rates are {approximately}10 times those of the natural state, agricultural land retirement by urbanization will severely curtail the current safe-yields and promote overdraft pumping. Such overdrafting has occurred in the Sacramento metropolitan area for {approximately}40 years, creating cones of depression {approximately}25m deep. Today, groundwater withdrawal in Sacramento is approximately matched by infiltration of low {sup 18}O water (-11.0{per_thousand}) away from the Sacramento and American Rivers, which is estimated to occur at 100-300m/year from the sharp {sup 18}O gradients in our groundwater isotope map.

  4. Update to the Ground-Water Withdrawals Database for the Death Valley Regional Ground-Water Flow System, Nevada and California, 1913-2003

    USGS Publications Warehouse

    Moreo, Michael T.; Justet, Leigh

    2008-01-01

    Ground-water withdrawal estimates from 1913 through 2003 for the Death Valley regional ground-water flow system are compiled in an electronic database to support a regional, three-dimensional, transient ground-water flow model. This database updates a previously published database that compiled estimates of ground-water withdrawals for 1913-1998. The same methodology is used to construct each database. Primary differences between the 2 databases are an additional 5 years of ground-water withdrawal data, well locations in the updated database are restricted to Death Valley regional ground-water flow system model boundary, and application rates are from 0 to 1.5 feet per year lower than original estimates. The lower application rates result from revised estimates of crop consumptive use, which are based on updated estimates of potential evapotranspiration. In 2003, about 55,700 acre-feet of ground water was pumped in the DVRFS, of which 69 percent was used for irrigation, 13 percent for domestic, and 18 percent for public supply, commercial, and mining activities.

  5. Update to the Ground-Water Withdrawals Database for the Death Valley REgional Ground-Water Flow System, Nevada and California, 1913-2003

    SciTech Connect

    Michael T. Moreo; and Leigh Justet

    2008-07-02

    Ground-water withdrawal estimates from 1913 through 2003 for the Death Valley regional ground-water flow system are compiled in an electronic database to support a regional, three-dimensional, transient ground-water flow model. This database updates a previously published database that compiled estimates of ground-water withdrawals for 1913–1998. The same methodology is used to construct each database. Primary differences between the 2 databases are an additional 5 years of ground-water withdrawal data, well locations in the updated database are restricted to Death Valley regional ground-water flow system model boundary, and application rates are from 0 to 1.5 feet per year lower than original estimates. The lower application rates result from revised estimates of crop consumptive use, which are based on updated estimates of potential evapotranspiration. In 2003, about 55,700 acre-feet of ground water was pumped in the DVRFS, of which 69 percent was used for irrigation, 13 percent for domestic, and 18 percent for public supply, commercial, and mining activities.

  6. Glaciation and regional groundwater flow in the Fennoscandian shield

    USGS Publications Warehouse

    Provost, A.M.; Voss, C.I.; Neuzil, C.E.

    2012-01-01

    Regional-scale groundwater flow modeling of the Fennoscandian shield suggests that groundwater flow can be strongly affected by future climate change and glaciation. We considered variable-density groundwater flow in a 1500-km-long and approximately 10-km-deep cross-section through southern Sweden. Groundwater flow and shield brine transport in the cross-sectional model were analyzed under projected surface conditions for the next 140 ka. Simulations suggest that blockage of recharge and discharge by low-permeability permafrost or cold-based ice causes sinking of brine and consequent freshening of near-surface water in areas of natural discharge. Although recharge of basal meltwater is limited by the requirement that water pressure at the base of the ice sheet not exceed the pressure exerted by the weight of the ice, warm-based ice with basal melting creates a potential for groundwater recharge rates much larger than those of present, ice-free conditions. In the simulations, regional-scale redistribution of recharged water by subsurface flow is minor over the duration of a glacial advance (approximately 10 ka). During glacial retreat, significant upward flow of groundwater may occur below the ice sheet owing to pressure release. If the mechanical loading efficiency of the rocks is high, both subsurface penetration of meltwater during glacial advance and up-flow during glacial retreat are reduced because of loading-induced pressure changes. The maximum rate of groundwater discharge in the simulations occurs at the receding ice margin, and some discharge occurs below incursive postglacial seas. Recharge of basal meltwater could decrease the concentration of dissolved solids significantly below present-day levels at depths of up to several kilometers and may bring oxygenated conditions to an otherwise reducing chemical environment for periods exceeding 10 ka.

  7. Groundwater storage trends in the Loess Plateau of China estimated from streamflow records

    NASA Astrophysics Data System (ADS)

    Gao, Zhaoliang; Zhang, Lu; Cheng, Lei; Zhang, Xiaoping; Cowan, Tim; Cai, Wenju; Brutsaert, Wilfried

    2015-11-01

    The catchments in the Loess Plateau in China have experienced significant land use change since the 1950s with a great number of soil conservation measures such as revegetation being implemented. Such soil conservation measures and climate variability have had considerable impacts on annual streamflow from these catchments. However, much less is known about changes in groundwater storage as the period of direct groundwater storage measurements is too short to reliably infer groundwater storage trends. For this study, annual values of groundwater storage from 38 catchments in the Loess Plateau were estimated from daily streamflow records based on groundwater flow theory. It was found that over the period of record (viz. 1955-2010), statistically significant (p < 0.1) downward trends have been identified in 20 selected catchments with an average reduction of -0.0299 mm per year, mostly located in the northern part of the Loess Plateau. Upward groundwater storage trends were observed in 10 catchments with an average increase of 0.00467 mm per year; these upward trends occurred in southern parts of the study area. Groundwater storage showed no statistically significant trends in 8 out of the 38 selected catchments. Soil conservation measures implemented in the Loess Plateau such as large-scale revegetation may have contributed to the estimated groundwater storage trends. Changes in sea surface temperature in the tropical Pacific Ocean, as indicated by shifts in climate variability modes such as El Niño-Southern Oscillation and the Pacific Decadal Oscillation, appear to have also contributed to the decreasing trends in groundwater storage in this region.

  8. Computer programs for describing the recession of ground-water discharge and for estimating mean ground-water recharge and discharge from streamflow records-update

    USGS Publications Warehouse

    Rutledge, A.T.

    1998-01-01

    The computer programs included in this report can be used to develop a mathematical expression for recession of ground-water discharge and estimate mean ground-water recharge and discharge. The programs are intended for analysis of the daily streamflow record of a basin where one can reasonably assume that all, or nearly all, ground water discharges to the stream except for that which is lost to riparian evapotranspiration, and where regulation and diversion of flow can be considered to be negligible. The program RECESS determines the master reces-sion curve of streamflow recession during times when all flow can be considered to be ground-water discharge and when the profile of the ground-water-head distribution is nearly stable. The method uses a repetitive interactive procedure for selecting several periods of continuous recession, and it allows for nonlinearity in the relation between time and the logarithm of flow. The program RORA uses the recession-curve displacement method to estimate the recharge for each peak in the streamflow record. The method is based on the change in the total potential ground-water discharge that is caused by an event. Program RORA is applied to a long period of record to obtain an estimate of the mean rate of ground-water recharge. The program PART uses streamflow partitioning to estimate a daily record of base flow under the streamflow record. The method designates base flow to be equal to streamflow on days that fit a requirement of antecedent recession, linearly interpolates base flow for other days, and is applied to a long period of record to obtain an estimate of the mean rate of ground-water discharge. The results of programs RORA and PART correlate well with each other and compare reasonably with results of the corresponding manual method.

  9. Incorporation of prior information on parameters into nonlinear regression groundwater flow models 2. Applications.

    USGS Publications Warehouse

    Cooley, R.L.

    1983-01-01

    Investigates factors influencing the degree of improvement in estimates of parameters of a nonlinear regression groundwater flow model by incorporating prior information of unknown reliability. Consideration of expected behavior of the regression solutions and results of a hypothetical modeling problem lead to several general conclusions. -from Author

  10. Multivariate analyses with end-member mixing to characterize groundwater flow: Wind Cave and associated aquifers

    USGS Publications Warehouse

    Long, Andrew J.; Valder, Joshua F.

    2011-01-01

    Principal component analysis (PCA) applied to hydrochemical data has been used with end-member mixing to characterize groundwater flow to a limited extent, but aspects of this approach are unresolved. Previous similar approaches typically have assumed that the extreme-value samples identified by PCA represent end members. The method presented herein is different from previous work in that (1) end members were not assumed to have been sampled but rather were estimated and constrained by prior knowledge; (2) end-member mixing was quantified in relation to hydrogeologic domains, which focuses model results on major hydrologic processes; (3) a method to select an appropriate number of end members using a series of cluster analyses is presented; and (4) conservative tracers were weighted preferentially in model calibration, which distributed model errors of optimized values, or residuals, more appropriately than would otherwise be the case. The latter item also provides an estimate of the relative influence of geochemical evolution along flow paths in comparison to mixing. This method was applied to groundwater in Wind Cave and the associated karst aquifer in the Black Hills of South Dakota, USA. The end-member mixing model was used to test a hypothesis that five different end-member waters are mixed in the groundwater system comprising five hydrogeologic domains. The model estimated that Wind Cave received most of its groundwater inflow from local surface recharge with an additional 33% from an upgradient aquifer. Artesian springs in the vicinity of Wind Cave primarily received water from regional groundwater flow.

  11. 2007 Estimated International Energy Flows

    SciTech Connect

    Smith, C A; Belles, R D; Simon, A J

    2011-03-10

    An energy flow chart or 'atlas' for 136 countries has been constructed from data maintained by the International Energy Agency (IEA) and estimates of energy use patterns for the year 2007. Approximately 490 exajoules (460 quadrillion BTU) of primary energy are used in aggregate by these countries each year. While the basic structure of the energy system is consistent from country to country, patterns of resource use and consumption vary. Energy can be visualized as it flows from resources (i.e. coal, petroleum, natural gas) through transformations such as electricity generation to end uses (i.e. residential, commercial, industrial, transportation). These flow patterns are visualized in this atlas of 136 country-level energy flow charts.

  12. Estimation of groundwater velocities from Yucca Flat to the Amargosa Desert using geochemistry and environmental isotopes

    SciTech Connect

    Hershey, R.L.; Acheampong, S.Y.

    1997-06-01

    Geochemical and isotopic data from groundwater sampling locations can be used to estimate groundwater flow velocities for independent comparison to velocities calculated by other methods. The objective of this study was to calculate groundwater flow velocities using geochemistry and environmental isotopes from the southern end of Yucca Flat to the Amargosa Desert, considering mixing of different groundwater inputs from sources each and southeast of the Nevada Test Site (NTS). The approach used to accomplish the objective of this study consisted of five steps: (1) reviewing and selecting locations where carbon isotopic groundwater analyses, reliable ionic analysis, and well completion information are available; (2) calculating chemical speciation with the computer code WATEQ4F (Ball and Nordstrom, 1991) to determine the saturation state of mineral phases for each ground water location; (3) grouping wells into reasonable flowpaths and mixing scenarios from different groundwater sources; (4) using the computer code NETPATH (Plummer et al., 1991) to simulate mixing and the possible chemical reactions along the flowpath, and to calculate the changes in carbon-13/carbon-12 isotopic ratios ({delta}{sup 13}C) as a result of these reactions; and (5) using carbon-14 ({sup 14}C) data to calculate velocity.

  13. Regional ground-water flow in the Lower Peninsula of Michigan

    SciTech Connect

    Grannemann, N.G.; Huffman, G.C. )

    1994-04-01

    A steady-state, numerical model is being developed to simulate groundwater flow in four regional aquifers in Michigan's Lower Peninsula. The uppermost layer in the model simulates flow in the glaciofluvial aquifer and the second layer simulates flow in the Saginaw aquifer, both of which contain freshwater. The lower two modeled units simulate flow in the Parma-Bayport and Marshall aquifers, both of which contain saline water or brine, except at or near their subcrop, where they contain freshwater. The US Geological Survey's Modular Model (MODFLOW) was modified to simulate variable-density groundwater flow in the lower two aquifers by assuming that groundwater density and viscosity differ from place to place but do not change over time. The model simulates groundwater conditions prior to large-scale withdrawals from the aquifer. Boundaries for the model include constant-head boundaries at the shorelines of Lakes Michigan, Huron, St. Clair, and Erie, as well as the St. Clair and Detroit Rivers. The southern boundary is simulated as a no-flow condition along several major stream divides in Michigan and Indiana. Recharge to the glaciofluvial aquifer is estimated to range from 0.2 to 22 inches per year and averages 8.4 inches per year. Regional groundwater flow occurs in areas where the bedrock aquifers are confined and in parts of the glaciofluvial aquifer, such as in the north-central part of the Lower Peninsula. Regional discharge primarily occurs to Saginaw Bay and to streams in the Saginaw and Michigan Lowlands. Parts of the Grand and Maple rivers may also receive discharge on a regional scale. Local flow systems control the direction and rate of groundwater flow in areas where aquifers are unconfined or hydraulically connected to overlying glacial deposits.

  14. Recharge and Lateral Groundwater Flow Boundary Conditions for the Saturated Zone Site-Scale Flow and Transport Model

    SciTech Connect

    B. Arnold; T. Corbet

    2001-12-18

    The purpose of the flow boundary conditions analysis is to provide specified-flux boundary conditions for the saturated zone (SZ) site-scale flow and transport model. This analysis is designed to use existing modeling and analysis results as the basis for estimated groundwater flow rates into the SZ site-scale model domain, both as recharge at the upper (water table) boundary and as underflow at the lateral boundaries. The objective is to provide consistency at the boundaries between the SZ site-scale flow model and other groundwater flow models. The scope of this analysis includes extraction of the volumetric groundwater flow rates simulated by the SZ regional-scale flow model to occur at the lateral boundaries of the SZ site-scale flow model and the internal qualification of the regional-scale model for use in this analysis model report (AMR). In addition, the scope includes compilation of information on the recharge boundary condition taken from three sources: (1) distributed recharge as taken from the SZ regional-scale flow model, (2) recharge below the area of the unsaturated zone (UZ) site-scale flow model, and (3) focused recharge along the Fortymile Wash channel.

  15. Estimation of Groundwater Recharge at Pahute Mesa using the Chloride Mass-Balance Method

    SciTech Connect

    Cooper, Clay A; Hershey, Ronald L; Healey, John M; Lyles, Brad F

    2013-07-01

    Groundwater recharge on Pahute Mesa was estimated using the chloride mass-balance (CMB) method. This method relies on the conservative properties of chloride to trace its movement from the atmosphere as dry- and wet-deposition through the soil zone and ultimately to the saturated zone. Typically, the CMB method assumes no mixing of groundwater with different chloride concentrations; however, because groundwater is thought to flow into Pahute Mesa from valleys north of Pahute Mesa, groundwater flow rates (i.e., underflow) and chloride concentrations from Kawich Valley and Gold Flat were carefully considered. Precipitation was measured with bulk and tipping-bucket precipitation gauges installed for this study at six sites on Pahute Mesa. These data, along with historical precipitation amounts from gauges on Pahute Mesa and estimates from the PRISM model, were evaluated to estimate mean annual precipitation. Chloride deposition from the atmosphere was estimated by analyzing quarterly samples of wet- and dry-deposition for chloride in the bulk gauges and evaluating chloride wet-deposition amounts measured at other locations by the National Atmospheric Deposition Program. Mean chloride concentrations in groundwater were estimated using data from the UGTA Geochemistry Database, data from other reports, and data from samples collected from emplacement boreholes for this study. Calculations were conducted assuming both no underflow and underflow from Kawich Valley and Gold Flat. Model results estimate recharge to be 30 mm/yr with a standard deviation of 18 mm/yr on Pahute Mesa, for elevations >1800 m amsl. These estimates assume Pahute Mesa recharge mixes completely with underflow from Kawich Valley and Gold Flat. The model assumes that precipitation, chloride concentration in bulk deposition, underflow and its chloride concentration, have been constant over the length of time of recharge.

  16. Estimation of microbial respiration rates in groundwater by geochemical modeling constrained with stable isotopes

    SciTech Connect

    Murphy, E.M.; Schramke, J.A.

    1998-11-01

    Changes in geochemistry and stable isotopes along a well-established groundwater flow path were used to estimate in situ microbial respiration rates in the Middendorf aquifer in the southeastern United States. Respiration rates were determined for individual terminal electron acceptors including O{sub 2}, MnO{sub 2}, Fe{sup 3+}, and SO{sub 4}{sup 2{minus}}. The extent of biotic reactions were constrained by the fractionation of stable isotopes of carbon and sulfur. Sulfur isotopes and the presence of sulfur-oxidizing microorganisms indicated that sulfate is produced through the oxidation of reduced sulfur species in the aquifer and not by the dissolution of gypsum, as previously reported. The respiration rates varied along the flow path as the groundwater transitioned between primarily oxic to anoxic conditions. Iron-reducing microorganisms were the largest contributors to the oxidation of organic matter along the portion of the groundwater flow path investigated in this study. The transition zone between oxic and anoxic groundwater contained a wide range of terminal electron acceptors and showed the greatest diversity and numbers of culturable microorganisms and the highest respiration rates. A comparison of respiration rates measured from core samples and pumped groundwater suggests that variability in respiration rates may often reflect the measurement scales, both in the sample volume and the time-frame over which the respiration measurement is averaged. Chemical heterogeneity may create a wide range of respiration rates when the scale of the observation is below the scale of the heterogeneity.

  17. Stream Flow Estimation via Belief Propagation for Sparsely Instrumented Watersheds

    NASA Astrophysics Data System (ADS)

    Krekeler, C. R.; Nagarajan, K.; Graham, W. D.; Slatton, K. C.

    2009-12-01

    Knowledge of stream flow rates is critical for management of in-stream flows and mitigation of flooding and drought events. Unfortunately, spatially dense networks of in situ stream flow measurements are not generally available and would be prohibitively expensive to deploy and maintain. Since measurements are only available at relatively sparse spatial frequencies, a data assimilation technique that best utilizes available measurements in a computationally efficient manner is required. Complexities in underlying geology, groundwater-surface water interactions, and rainfall patterns also need to be captured for optimal performance. In this study, a probabilistic technique based on Bayesian Networks and belief propagation is used to estimate and forecast stream flow based on flow at surrounding locations, rainfall, and groundwater levels, and to quantify the uncertainty of the estimates. The causal nature and inherent tree-like structure of stream flow suggests that belief propagation based on message passing is a computationally efficient method to propagate partial observations to improve flow estimates at other stations across the watershed. Stations are considered as part of an array of nodes that communicate evidence of flow measurements to their neighbors via conditional probability densities defining the links between the nodes. The probability density functions are generated based on the relationship between estimates of flow from the Watershed Assessment Model (WAM), calibrated for our study site at the Santa Fe River Watershed in North Central Florida, USA from 1990 through 2008. Rainfall values were obtained from NEXRAD datasets, while Suwannee River Water Management District provided groundwater levels at a network of wells in the watershed. Performance of the proposed method is evaluated by comparing the results to flow measured at several USGS gage stations along the river not used to develop the algorithm, using various metrics such as Nash

  18. Numerical simulations of groundwater flow at New Jersey Shallow Shelf

    NASA Astrophysics Data System (ADS)

    Fehr, Annick; Patterson, Fabian; Lofi, Johanna; Reiche, Sönke

    2016-04-01

    During IODP Expedition 313, three boreholes were drilled in the so-called New Jersey transect. Hydrochemical studies revealed the groundwater situation as more complex than expected, characterized by several sharp boundaries between fresh and saline groundwater. Two conflicting hypotheses regarding the nature of these freshwater reservoirs are currently debated. One hypothesis is that these reservoirs are connected with onshore aquifers and continuously recharged by seaward-flowing groundwater. The second hypothesis is that fresh groundwater was emplaced during the last glacial period. In addition to the petrophysical properties measured during IODP 313 expedition, Nuclear Magnetic Resonance (NMR) measurements were performed on samples from boreholes M0027, M0028 and M0029 in order to deduce porosities and permeabilities. These results are compared with data from alternative laboratory measurements and with petrophysical properties inferred from downhole logging data. We incorporate these results into a 2D numerical model that reflects the shelf architecture as known from drillings and seismic data to perform submarine groundwater flow simulations. In order to account for uncertainties related to the spatial distribution of physical properties, such as porosity and permeability, systematic variation of input parameters was performed during simulation runs. The target is to test the two conflicting hypotheses of fresh groundwater emplacements offshore New Jersey and to improve the understanding of fluid flow processes at marine passive margins.

  19. Conceptual model and numerical simulation of the groundwater-flow system of Bainbridge Island, Washington

    USGS Publications Warehouse

    Frans, Lonna M.; Bachmann, Matthew P.; Sumioka, Steve S.; Olsen, Theresa D.

    2011-01-01

    Groundwater is the sole source of drinking water for the population of Bainbridge Island. Increased use of groundwater supplies on Bainbridge Island as the population has grown over time has created concern about the quantity of water available and whether saltwater intrusion will occur as groundwater usage increases. A groundwater-flow model was developed to aid in the understanding of the groundwater system and the effects of groundwater development alternatives on the water resources of Bainbridge Island. Bainbridge Island is underlain by unconsolidated deposits of glacial and nonglacial origin. The surficial geologic units and the deposits at depth were differentiated into aquifers and confining units on the basis of areal extent and general water-bearing characteristics. Eleven principal hydrogeologic units are recognized in the study area and form the basis of the groundwater-flow model. A transient variable-density groundwater-flow model of Bainbridge Island and the surrounding area was developed to simulate current (2008) groundwater conditions. The model was calibrated to water levels measured during 2007 and 2008 using parameter estimation (PEST) to minimize the weighted differences or residuals between simulated and measured hydraulic head. The calibrated model was used to make some general observations of the groundwater system in 2008. Total flow through the groundwater system was about 31,000 acre-ft/ yr. The recharge to the groundwater system was from precipitation and septic-system returns. Groundwater flow to Bainbridge Island accounted for about 1,000 acre-ft/ yr or slightly more than 5 percent of the recharge amounts. Groundwater discharge was predominately to streams, lakes, springs, and seepage faces (16,000 acre-ft/yr) and directly to marine waters (10,000 acre-ft/yr). Total groundwater withdrawals in 2008 were slightly more than 6 percent (2,000 acre-ft/yr) of the total flow. The calibrated model was used to simulate predevelopment conditions

  20. SWB-A modified Thornthwaite-Mather Soil-Water-Balance code for estimating groundwater recharge

    USGS Publications Warehouse

    Westenbroek, S.M.; Kelson, V.A.; Dripps, W.R.; Hunt, R.J.; Bradbury, K.R.

    2010-01-01

    A Soil-Water-Balance (SWB) computer code has been developed to calculate spatial and temporal variations in groundwater recharge. The SWB model calculates recharge by use of commonly available geographic information system (GIS) data layers in combination with tabular climatological data. The code is based on a modified Thornthwaite-Mather soil-water-balance approach, with components of the soil-water balance calculated at a daily timestep. Recharge calculations are made on a rectangular grid of computational elements that may be easily imported into a regional groundwater-flow model. Recharge estimates calculated by the code may be output as daily, monthly, or annual values.

  1. Indications of regional scale groundwater flows in the Amazon Basins: Inferences from results of geothermal studies

    NASA Astrophysics Data System (ADS)

    Pimentel, Elizabeth T.; Hamza, Valiya M.

    2012-08-01

    The present work deals with determination groundwater flows in the Amazon region, based on analysis of geothermal data acquired in shallow and deep wells. The method employed is based on the model of simultaneous heat transfer by conduction and advection in permeable media. Analysis of temperature data acquired in water wells indicates down flows of groundwaters with velocities in excess of 10-7 m/s at depths less than 300 m in the Amazonas basin. Bottom-hole temperature (BHT) data sets have been used in determining characteristics of fluid movements at larger depths in the basins of Acre, Solimões, Amazonas, Marajó and Barreirinhas. The results of model simulations point to down flow of groundwaters with velocities of the order of 10-8 to 10-9 m/s, at depths of up to 4000 m. No evidence has been found for up flow typical of discharge zones. The general conclusion compatible with such results is that large-scale groundwater recharge systems operate at both shallow and deep levels in all sedimentary basins of the Amazon region. However, the basement rock formations of the Amazon region are relatively impermeable and hence extensive down flow systems through the sedimentary strata are possible only in the presence of generalized lateral movement of groundwater in the basal parts of the sedimentary basins. The direction of this lateral flow, inferred from the basement topography and geological characteristics of the region, is from west to east, following roughly the course of surface drainage system of the Amazon River, with eventual discharge into the Atlantic Ocean. The estimated flow rate at the continental margin is 3287 m3/s, with velocities of the order of 218 m/year. It is possible that dynamic changes in the fluvial systems in the western parts of South American continent have been responsible for triggering alterations in the groundwater recharge systems and deep seated lateral flows in the Amazon region.

  2. Evidence for deep groundwater flow and convective heat transport in mountainous terrain, Delta County, Colorado, USA

    NASA Astrophysics Data System (ADS)

    Lazear, Gregory D.

    2006-12-01

    The Tongue Creek watershed lies on the south flank of Grand Mesa in western Colorado, USA and is a site with 1.5 km of topographic relief, heat flow of 100 mW/m2, thermal conductivity of 3.3 W m-1 °C-1, hydraulic conductivity of 10-8 m/s, a water table that closely follows surface topography, and groundwater temperatures 3-15°C above mean surface temperatures. These data suggest that convective heat transport by groundwater flow has modified the thermal regime of the site. Steady state three-dimensional numerical simulations of heat flow, groundwater flow, and convective transport were used to model these thermal and hydrological data. The simulations provided estimates for the scale of hydraulic conductivity and bedrock base flow discharge within the watershed. The numerical models show that (1) complex three-dimensional flow systems develop with a range of scales from tens of meters to tens of kilometers; (2) mapped springs are frequently found at locations where contours of hydraulic head indicate strong vertical flow at the water table, and; (3) the distribution of groundwater temperatures in water wells as a function of surface elevation is predicted by the model.

  3. A guide for using the transient ground-water flow model of the Death Valley regional ground-water flow system, Nevada and California

    USGS Publications Warehouse

    Blainey, Joan B.; Faunt, Claudia C.; Hill, Mary C.

    2006-01-01

    This report is a guide for executing numerical simulations with the transient ground-water flow model of the Death Valley regional ground-water flow system, Nevada and California using the U.S. Geological Survey modular finite-difference ground-water flow model, MODFLOW-2000. Model inputs, including observations of hydraulic head, discharge, and boundary flows, are summarized. Modification of the DVRFS transient ground-water model is discussed for two common uses of the Death Valley regional ground-water flow system model: predictive pumping scenarios that extend beyond the end of the model simulation period (1998), and model simulations with only steady-state conditions.

  4. A Guide for Using the Transient Ground-Water Flow Model of the Death Valley Regional Ground-Water Flow System, Nevada and California

    SciTech Connect

    Joan B. Blainey; Claudia C. Faunt, and Mary C. Hill

    2006-05-16

    This report is a guide for executing numerical simulations with the transient ground-water flow model of the Death Valley regional ground-water flow system, Nevada and California using the U.S. Geological Survey modular finite-difference ground-water flow model, MODFLOW-2000. Model inputs, including observations of hydraulic head, discharge, and boundary flows, are summarized. Modification of the DVRFS transient ground-water model is discussed for two common uses of the Death Valley regional ground-water flow system model: predictive pumping scenarios that extend beyond the end of the model simulation period (1998), and model simulations with only steady-state conditions.

  5. Groundwater Flow Dynamic Simulations of a Buried Valley Aquifer Calibrated with Field and Remotely Sensed Data

    NASA Astrophysics Data System (ADS)

    Calderhead, A. I.; Hinton, M. J.; Logan, C. E.; Sharpe, D.; Russel, H. A.; Oldenborger, G. A.; Pugin, A.; Rivera, A.; Castellazzi, P.; Martel, R.

    2013-12-01

    Buried valleys are a common occurrence in the North American prairie landscape. They are often characterized as high yield sources of groundwater in regions where low yield shale and tills dominate the hydrogeological setting. Firstly, 3D conceptual and geological models have been generated and used as a basis for creating a 3D finite element groundwater flow model. Field data, including piezometric readings, base flow measurements, and soil moisture probe data were collected between 2011 and 2013 and are used for calibrating the flow model. Secondly, the study aims to improve the spatial discretization of recharge estimates and include these refined values in the flow model. A temporal series of C-band Radar data and several land surface models were compared with the soil moisture probe data from the Spiritwood buried valley aquifer. The radar backscatter was used to develop moisture estimates at the regional scale. These estimates were then input into the HELP multi-parameter recharge model with the aim of assisting in estimates of a spatial discretization for groundwater recharge. Preliminary groundwater simulation results, with uniform recharge, show good agreement with piezometer readings and measured base flow readings. The temporal series of C-band radar backscatter, moisture probe data, and land surface models show corresponding variations between October, 2011 and October, 2012. The high resolution and regional extent of the radar data has a high potential to help develop a better understanding of recharge patterns in buried valley settings. Integrating a temporal series of high-resolution data into conceptual and numerical model development will refine our mapping, understanding and assessment of buried valley aquifers. Future work will include incorporating the spatially variable recharge estimates into the 3D finite element flow model. Additionally, various interpretations of the geological model will be tested to determine the extent, if any, that a

  6. Ground-water flow in the shallow aquifer system at the Naval Weapons Station Yorktown, Virginia

    USGS Publications Warehouse

    Smith, Barry S.

    2001-01-01

    The Environmental Directorate of the Naval Weapons Station Yorktown, Virginia, is concerned about possible contamination of ground water at the Station. Ground water at the Station flows through a shallow system of layered aquifers and leaky confining units. The units of the shallow aquifer system are the Columbia aquifer, the Cornwallis Cave confining unit, the Cornwallis Cave aquifer, the Yorktown confining unit, and the Yorktown-Eastover aquifer. The Eastover-Calvert confining unit separates the shallow aquifer system from deeper confined aquifers beneath the Station. A three-dimensional, finite-difference, ground-water flow model was used to simulate steady-state ground-water flow of the shallow aquifer system in and around the Station. The model simulated ground-water flow from the peninsular drainage divide that runs across the Lackey Plain near the southern end of the Station north to King Creek and the York River and south to Skiffes Creek and the James River. The model was calibrated by minimizing the root mean square error between 4 7 measured and corresponding simulated water levels. The calibrated model was used to determine the ground-water budget and general directions of ground-water flow. A particle-tracking routine was used with the calibrated model to estimate groundwater flow paths, flow rates, and traveltimes from selected sites at the Station. Simulated ground-water flow velocities of the Station-area model were small beneath the interstream areas of the Lackey Plain and Croaker Flat, but increased outward toward the streams and rivers where the hydraulic gradients are larger. If contaminants from the land surface entered the water table at or near the interstream areas of the Station, where hydraulic gradients are smaller, they would migrate more slowly than if they entered closer to the streams or the shores of the rivers where gradients commonly are larger. The ground-water flow simulations indicate that some ground water leaks downward from

  7. An initial examination of tungsten geochemistry along groundwater flow paths

    NASA Astrophysics Data System (ADS)

    Dave, H. B.; Johannesson, K. H.

    2008-12-01

    Groundwater samples were collected along groundwater flow paths from the Upper Floridan (Florida), Carrizo Sand (Texas), and the Aquia (Maryland) aquifers and analyzed for tungsten (W) concentrations by high- resolution inductively couple plasma mass spectrometry. At each well head, groundwater samples were also analyzed for pH, specific conductance, temperature, alkalinity, dissolved oxygen (DO), oxidation-reduction potential (Eh), dissolved iron speciation, and dissolved sulfide [S(-II)] concentrations. Sediment samples from the Carrizo Sand and Aquia aquifers were also collected and subjected to sequential extractions to provide additional insights into the solid-phase speciation of W in these aquifers. Tungsten concentrations varied along the groundwater flow paths chiefly in response to changing pH, and to a lesser extent, variations in the redox conditions. For groundwater from the Carrizo Sand aquifer, W ranges between 3.64 and 1297 pmol/kg, exhibiting the lowest values proximal to the recharge zone. Tungsten concentrations progressively increase along the flow path, reaching 1297 pmol/kg in the sulfidic groundwaters located approximately 60 km downgradient from the recharge area. Tungsten is strongly correlated with S(-II) concentrations and pH in Carrizo groundwaters (r = 0.95 and 0.78, respectively). Within the Aquia aquifer, however, W generally occurs at lower concentrations than the Carrizo (14 to 184 pmol/kg; mean = 80 pmol/kg), and shows no systematic trends along the flow path (e.g., r = 0.08 and 0.4 for W vs. S(-II) and pH, respectively). Our data are consistent with the increase in W concentrations in Carrizo groundwaters reflecting, in part, pH-related desorption, which has been shown to be substantial for pH greater than 8. Moreover, because of the broad similarities in the chemistry of W and Mo, which forms thiomolybdates in sulfidic waters, we suggest that thiotungstate complexes may form in sulfidic groundwaters, thus partially explaining the

  8. Estimating groundwater exchange with lakes: 1. The stable isotope mass balance method

    USGS Publications Warehouse

    Krabbenhoft, David P.; Bowser, Carl J.; Anderson, Mary P.; Valley, John W.

    1990-01-01

    Groundwater inflow and outflow contributions to the hydrologic budget of lakes can be determined using a stable isotope (18O/16O) mass balance method. The stable isotope method provides a way of integrating the spatial and temporal complexities of the flow field around a lake, thereby offering an appealing alternative to the traditional time and labor intensive methods using seepage meters and an extensive piezometer network. In this paper the method is applied to a lake in northern Wisconsin, demonstrating that it can be successfully applied to lakes in the upper midwest where thousands of similar lakes exist. Inflow and outflow rates calculated for the Wisconsin lake using the isotope mass balance method are 29 and 54 cm/yr, respectively, which compare well to estimates, derived independently using a three-dimensional groundwater flow and solute transport model, of 20 and 50 cm/yr. Such a favorable comparison lends confidence to the use of the stable isotope method to estimate groundwater exchange with lakes. In addition, utilization of stable isotopes in studies of groundwater-lake systems lends insight into mixing processes occurring in the unsaturated zone and in the aquifer surrounding the lake and verifies assumed flow paths based on head measurements in piezometers.

  9. An improved time series approach for estimating groundwater recharge from groundwater level fluctuations

    NASA Astrophysics Data System (ADS)

    Cuthbert, M. O.

    2010-09-01

    An analytical solution to a linearized Boussinesq equation is extended to develop an expression for groundwater drainage using estimations of aquifer parameters. This is then used to develop an improved water table fluctuation (WTF) technique for estimating groundwater recharge. The resulting method extends the standard WTF technique by making it applicable, as long as aquifer properties for the area are relatively well known, in areas with smoothly varying water tables and is not reliant on precipitation data. The method is validated against numerical simulations and a case study from a catchment where recharge is "known" a priori using other means. The approach may also be inverted to provide initial estimates of aquifer parameters in areas where recharge can be reliably estimated by other methods.

  10. Flow calculations for Yucca Mountain groundwater travel time (GWTT-95)

    SciTech Connect

    Altman, S.J.; Arnold, B.W.; Barnard, R.W.; Barr, G.E.; Ho, C.K.; McKenna, S.A.; Eaton, R.R.

    1996-09-01

    In 1983, high-level radioactive waste repository performance requirements related to groundwater travel time were defined by NRC subsystem regulation 10 CFR 60.113. Although DOE is not presently attempting to demonstrate compliance with that regulation, understanding of the prevalence of fast paths in the groundwater flow system remains a critical element of any safety analyses for a potential repository system at Yucca Mountain, Nevada. Therefore, this analysis was performed to allow comparison of fast-path flow against the criteria set forth in the regulation. Models developed to describe the conditions for initiation, propagation, and sustainability of rapid groundwater movement in both the unsaturated and saturated zones will form part of the technical basis for total- system analyses to assess site viability and site licensability. One of the most significant findings is that the fastest travel times in both unsaturated and saturated zones are in the southern portion of the potential repository, so it is recommended that site characterization studies concentrate on this area. Results support the assumptions regarding the importance of an appropriate conceptual model of groundwater flow and the incorporation of heterogeneous material properties into the analyses. Groundwater travel times are sensitive to variation/uncertainty in hydrologic parameters and in infiltration flux at upper boundary of the problem domain. Simulated travel times are also sensitive to poorly constrained parameters of the interaction between flow in fractures and in the matrix.

  11. Estimating seepage flux from ephemeral stream channels using surface water and groundwater level data

    NASA Astrophysics Data System (ADS)

    Noorduijn, Saskia L.; Shanafield, Margaret; Trigg, Mark A.; Harrington, Glenn A.; Cook, Peter G.; Peeters, L.

    2014-02-01

    Seepage flux from ephemeral streams can be an important component of the water balance in arid and semiarid regions. An emerging technique for quantifying this flux involves the measurement and simulation of a flood wave as it moves along an initially dry channel. This study investigates the usefulness of including surface water and groundwater data to improve model calibration when using this technique. We trialed this approach using a controlled flow event along a 1387 m reach of artificial stream channel. Observations were then simulated using a numerical model that combines the diffusion-wave approximation of the Saint-Vénant equations for streamflow routing, with Philip's infiltration equation and the groundwater flow equation. Model estimates of seepage flux for the upstream segments of the study reach, where streambed hydraulic conductivities were approximately 101 m d-1, were on the order of 10-4 m3 d-1 m-2. In the downstream segments, streambed hydraulic conductivities were generally much lower but highly variable (˜10-3 to 10-7 m d-1). A Latin Hypercube Monte Carlo sensitivity analysis showed that the flood front timing, surface water stage, groundwater heads, and the predicted streamflow seepage were most influenced by specific yield. Furthermore, inclusion of groundwater data resulted in a higher estimate of total seepage estimates than if the flood front timing were used alone.

  12. Estimating natural background groundwater chemistry, Questa molybdenum mine, New Mexico

    USGS Publications Warehouse

    Verplanck, Phillip L.; Nordstrom, D Kirk; Plumlee, Geoffrey S.; Walker, Bruce M.

    2010-01-01

    This 2 1/2 day field trip will present an overview of a U.S. Geological Survey (USGS) project whose objective was to estimate pre-mining groundwater chemistry at the Questa molybdenum mine, New Mexico. Because of intense debate among stakeholders regarding pre-mining groundwater chemistry standards, the New Mexico Environment Department and Chevron Mining Inc. (formerly Molycorp) agreed that the USGS should determine pre-mining groundwater quality at the site. In 2001, the USGS began a 5-year, multidisciplinary investigation to estimate pre-mining groundwater chemistry utilizing a detailed assessment of a proximal natural analog site and applied an interdisciplinary approach to infer pre-mining conditions. The trip will include a surface tour of the Questa mine and key locations in the erosion scar areas and along the Red River. The trip will provide participants with a detailed understanding of geochemical processes that influence pre-mining environmental baselines in mineralized areas and estimation techniques for determining pre-mining baseline conditions.

  13. Incorporation of prior information on parameters into nonlinear regression groundwater flow models. l. Theory.

    USGS Publications Warehouse

    Cooley, R.L.

    1982-01-01

    Prior information on the parameters of a groundwater flow model can be used to improve parameter estimates obtained from nonlinear regression solution of a modeling problem. Two scales of prior information can be available: 1) prior information having known reliability (that is, bias and random error structure), and 2) prior information consisting of best available estimates of unknown reliability. It is shown that if both scales of prior information are available, then a combined regression analysis may be made. -from Author

  14. Complex groundwater flow systems as traveling agent models

    PubMed Central

    Padilla, Pablo; Escolero, Oscar; González, Tomas; Morales-Casique, Eric; Osorio-Olvera, Luis

    2014-01-01

    Analyzing field data from pumping tests, we show that as with many other natural phenomena, groundwater flow exhibits complex dynamics described by 1/f power spectrum. This result is theoretically studied within an agent perspective. Using a traveling agent model, we prove that this statistical behavior emerges when the medium is complex. Some heuristic reasoning is provided to justify both spatial and dynamic complexity, as the result of the superposition of an infinite number of stochastic processes. Even more, we show that this implies that non-Kolmogorovian probability is needed for its study, and provide a set of new partial differential equations for groundwater flow. PMID:25337455

  15. Simulation of ground-water flow and areas contributing ground water to production wells, Cadillac, Michigan

    USGS Publications Warehouse

    Hoard, Christopher J.; Westjohn, David B.

    2005-01-01

    Ground water is the primary source of water for domestic, municipal, and industrial use within the northwest section of Michigan's Lower Peninsula. Because of the importance of this resource, numerous communities including the city of Cadillac in Wexford County, Michigan, have begun local wellhead protection programs. In these programs, communities protect their ground-water resources by identifying the areas that contribute water to production wells, identifying potential sources of contamination, and developing methods to cooperatively manage and minimize threats to the water supply. The U.S. Geological Survey, in cooperation with the city of Cadillac, simulated regional ground-water flow and estimated areas contributing recharge and zones of transport to the production well field. Ground-water flow models for the Clam River watershed, in Wexford and Missaukee Counties, were developed using the U.S. Geological Survey modular three-dimensional finite-difference ground-water flow model (MODFLOW 2000). Ground-water flow models were calibrated using the observation, sensitivity, and parameter estimation packages of MODFLOW 2000. Ground-water-head solutions from calibrated flow models were used in conjunction with MODPATH, a particle-tracking program, to simulate regional ground-water flow and estimate areas contributing recharge and zones of transport to the Cadillac production-well field for a 10-year period. Model simulations match the conceptual model in that regional ground-water flow in the deep ground-water system is from southeast to northwest across the watershed. Areas contributing water were determined for the optimized parameter set and an alternate parameter set that included increased recharge and hydraulic conductivity values. Although substantially different hydrologic parameters (assumed to represent end-member ranges of realistic hydrologic parameters) were used in alternate numerical simulations, simulation results differ little in predictions of

  16. An initial inverse calibration of the ground-water flow model for the Hanford unconfined aquifer

    SciTech Connect

    Jacobson, E.A. . Desert Research Inst.); Freshly, M.D. )

    1990-03-01

    Large volumes of process cooling water are discharged to the ground form U.S. Department of Energy (DOE) nuclear fuel processing operations in the central portion of the Hanford Site in southeastern Washington. Over the years, these large volumes of waste water have recharged the unconfined aquifer at the Site. This artificial recharge has affected ground-water levels and contaminant movement in the unconfined aquifer. Ground-water flow and contaminant transport models have been applied to assess the impacts of site operations on the rate and direction of ground-water flow and contaminant transport in unconfined aquifer at the Hanford Site. The inverse calibration method developed by Neuman and modified by Jacobson was applied to improve calibration of a ground-water flow model of the unconfined aquifer at the Hanford Site. All information about estimates of hydraulic properties of the aquifer, hydraulic heads, boundary conditions, and discharges to and withdrawals form the aquifer is included in the inverse method to obtain an initial calibration of the ground-water flow model. The purpose of this report is to provide a description of the inverse method, its initial application to the unconfined aquifer at Hanford, and to present results of the initial inverse calibration. 28 refs., 19 figs., 1 tab.

  17. A Method to Evaluate Groundwater flow system under the Seabed

    NASA Astrophysics Data System (ADS)

    Kohara, N.; Marui, A.

    2011-12-01

    / fresh water interface (position of the submarine groundwater discharge) may appear on the seafloor. Moreover, neither the salinity concentration nor the groundwater age depends on depth. It is thought that it is because that the groundwater forms the complex flow situation through the change in a long-term groundwater flow system. The technology to understand the coastal groundwater flow consists of remote sensing, geographical features analysis, surface of the earth investigation, geophysical exploration, drilling survey, and indoor examination and the measurement. Integration of each technology is needed to interpret groundwater flow system because the one is to catch the local groundwater flow in the time series and another one is to catch the long-term and regional groundwater flow in the general situation. The purpose of this study is to review the previous research of coastal groundwater flow, and to integrate an applicable evaluation approach to understand this mechanism. In this presentation, the review of the research and case study using numerical simulation are introduced.

  18. Control on groundwater flow in a semiarid folded and faulted intermountain basin

    USGS Publications Warehouse

    Ball, Lyndsay B.; Caine, Jonathan S.; Ge, Shemin

    2013-01-01

    The major processes controlling groundwater flow in intermountain basins are poorly understood, particularly in basins underlain by folded and faulted bedrock and under regionally realistic hydrogeologic heterogeneity. To explore the role of hydrogeologic heterogeneity and poorly constrained mountain hydrologic conditions on regional groundwater flow in contracted intermountain basins, a series of 3-D numerical groundwater flow models were developed using the South Park basin, Colorado, USA as a proxy. The models were used to identify the relative importance of different recharge processes to major aquifers, to estimate typical groundwater circulation depths, and to explore hydrogeologic communication between mountain and valley hydrogeologic landscapes. Modeling results show that mountain landscapes develop topographically controlled and predominantly local-scale to intermediate-scale flow systems. Permeability heterogeneity of the fold and fault belt and decreased topographic roughness led to permeability controlled flow systems in the valley. The structural position of major aquifers in the valley fold and fault belt was found to control the relative importance of different recharge mechanisms. Alternative mountain recharge model scenarios showed that higher mountain recharge rates led to higher mountain water table elevations and increasingly prominent local flow systems, primarily resulting in increased seepage within the mountain landscape and nonlinear increases in mountain block recharge to the valley. Valley aquifers were found to be relatively insensitive to changing mountain water tables, particularly in structurally isolated aquifers inside the fold and fault belt.

  19. Comparison of local- to regional-scale estimates of ground-water recharge in Minnesota, USA

    USGS Publications Warehouse

    Delin, G.N.; Healy, R.W.; Lorenz, D.L.; Nimmo, J.R.

    2007-01-01

    Regional ground-water recharge estimates for Minnesota were compared to estimates made on the basis of four local- and basin-scale methods. Three local-scale methods (unsaturated-zone water balance, water-table fluctuations (WTF) using three approaches, and age dating of ground water) yielded point estimates of recharge that represent spatial scales from about 1 to about 1000 m2. A fourth method (RORA, a basin-scale analysis of streamflow records using a recession-curve-displacement technique) yielded recharge estimates at a scale of 10–1000s of km2. The RORA basin-scale recharge estimates were regionalized to estimate recharge for the entire State of Minnesota on the basis of a regional regression recharge (RRR) model that also incorporated soil and climate data. Recharge rates estimated by the RRR model compared favorably to the local and basin-scale recharge estimates. RRR estimates at study locations were about 41% less on average than the unsaturated-zone water-balance estimates, ranged from 44% greater to 12% less than estimates that were based on the three WTF approaches, were about 4% less than the age dating of ground-water estimates, and were about 5% greater than the RORA estimates. Of the methods used in this study, the WTF method is the simplest and easiest to apply. Recharge estimates made on the basis of the UZWB method were inconsistent with the results from the other methods. Recharge estimates using the RRR model could be a good source of input for regional ground-water flow models; RRR model results currently are being applied for this purpose in USGS studies elsewhere.

  20. Comparison of local- to regional-scale estimates of ground-water recharge in Minnesota, USA

    NASA Astrophysics Data System (ADS)

    Delin, Geoffrey N.; Healy, Richard W.; Lorenz, David L.; Nimmo, John R.

    2007-02-01

    SummaryRegional ground-water recharge estimates for Minnesota were compared to estimates made on the basis of four local- and basin-scale methods. Three local-scale methods (unsaturated-zone water balance, water-table fluctuations (WTF) using three approaches, and age dating of ground water) yielded point estimates of recharge that represent spatial scales from about 1 to about 1000 m 2. A fourth method (RORA, a basin-scale analysis of streamflow records using a recession-curve-displacement technique) yielded recharge estimates at a scale of 10-1000s of km 2. The RORA basin-scale recharge estimates were regionalized to estimate recharge for the entire State of Minnesota on the basis of a regional regression recharge (RRR) model that also incorporated soil and climate data. Recharge rates estimated by the RRR model compared favorably to the local and basin-scale recharge estimates. RRR estimates at study locations were about 41% less on average than the unsaturated-zone water-balance estimates, ranged from 44% greater to 12% less than estimates that were based on the three WTF approaches, were about 4% less than the age dating of ground-water estimates, and were about 5% greater than the RORA estimates. Of the methods used in this study, the WTF method is the simplest and easiest to apply. Recharge estimates made on the basis of the UZWB method were inconsistent with the results from the other methods. Recharge estimates using the RRR model could be a good source of input for regional ground-water flow models; RRR model results currently are being applied for this purpose in USGS studies elsewhere.

  1. Intercomparison of Groundwater Flow Monitoring Technologies at Site OU 1, Former Fort Ord, California

    SciTech Connect

    Daley, P F; Jantos, J; Pedler, W H; Mandell, W A

    2005-09-20

    velocity measurements, and the ISPFS and SCBFM systems also gene generate flow direction rate estimates. The ISPFS probes are permanently installed and are non-retrievable, but produce long-term records with essentially no operator intervention or maintenance. The HPL and SCBFM systems are lightweight, portable logging devices that employ recording of electrical conductivity changes in wells purged with deionized water (HPL), or imaging of colloidal particles traversing the borehole (SCBFM) as the physical basis for estimating the velocity of groundwater flow through monitoring wells. All three devices gave estimates of groundwater velocity that were in reasonable agreement. However, although the ISPFS produced groundwater azimuth data that correlated well with conventional conductivity and gradient analyses of the groundwater flow field, the SCBFM direction data were in poor agreement. Further research into the reasons for this lack of correlation would seem to be warranted, given the ease of deployment of this tool in existing conventional monitoring wells, and its good agreement with the velocity estimates of the other technologies examined.

  2. Calculating discharge of phosphorus and nitrogen with groundwater base flow to a small urban stream reach

    NASA Astrophysics Data System (ADS)

    Fitzgerald, Alex; Roy, James W.; Smith, James E.

    2015-09-01

    Elevated levels of nutrients, especially phosphorus, in urban streams can lead to eutrophication and general degradation of stream water quality. Contributions of phosphorus from groundwater have typically been assumed minor, though elevated concentrations have been associated with riparian areas and urban settings. The objective of this study was to investigate the importance of groundwater as a pathway for phosphorus and nitrogen input to a gaining urban stream. The stream at the 28-m study reach was 3-5 m wide and straight, flowing generally eastward, with a relatively smooth bottom of predominantly sand, with some areas of finer sediments and a few boulders. Temperature-based methods were used to estimate the groundwater flux distribution. Detailed concentration distributions in discharging groundwater were mapped using in-stream piezometers and diffusion-based peepers, and showed elevated levels of soluble reactive phosphorus (SRP) and ammonium compared to the stream (while nitrate levels were lower), especially along the south bank, where groundwater fluxes were lower and geochemically reducing conditions dominated. Field evidence suggests the ammonium may originate from nearby landfills, but that local sediments likely contribute the SRP. Ammonium and SRP mass discharges with groundwater were then estimated as the product of the respective concentration distributions and the groundwater flux distribution. These were determined as approximately 9 and 200 g d-1 for SRP and ammonium, respectively, which compares to stream mass discharges over the observed range of base flows of 20-1100 and 270-7600 g d-1, respectively. This suggests that groundwater from this small reach, and any similar areas along Dyment's Creek, has the potential to contribute substantially to the stream nutrient concentrations.

  3. Consistency of groundwater flow patterns in mountainous topography: Implications for valley bottom water replenishment and for defining groundwater flow boundaries

    NASA Astrophysics Data System (ADS)

    Welch, L. A.; Allen, D. M.

    2012-05-01

    Topographic influences on groundwater flow processes that contribute to baseflow and mountain block recharge (MBR) are conceptually investigated using three-dimensional numerical models of saturated groundwater flow. Model domains for conceptual and real topographies are developed as "mountain groundwatershed units" (MGUs) to represent regional-scale watershed systems. Results indicate regularity in groundwater flow patterns that reflect consistency of prominent topographic features, providing a basis for conceptualizing three-dimensional groundwater flow. Baseflow is generated mainly from recharge within the watershed area. MBR is produced primarily from recharge that is focused across triangular facets near the mountain front (˜73%-97% of total MBR), with additional contributions originating within the watershed (up to ˜27% of MBR). MBR contributions originating from recharge near the highest-elevation watershed boundaries are minimal but are greater for topography with less stream incision. With orographic influences, more MBR originates within the watershed. MBR rates are relatively consistent between models because of similarities in mountain front topography, while baseflow is variable. Gains and losses to systems via cross-watershed groundwater flux, generated because of topographic differences between adjacent watersheds, cause baseflow to vary by up to ˜10% but do not significantly influence MBR. In data-sparse regions such as mountains, a basic numerical modeling approach, using the MGU concept with topography data and mapped watershed boundaries, can be used to develop site-specific conceptual models to constrain water budgets, to delineate recharge areas, and to guide further investigation and data collection.

  4. Implementation and influence of heterogeneous riverbed hydraulic conductivity in groundwater flow models

    NASA Astrophysics Data System (ADS)

    Ghysels, Gert; Huysmans, Marijke

    2015-04-01

    Characterization of groundwater-surface water exchange fluxes is important for assessing riparian ecology, determining quantity and quality of pumped groundwater close to rivers, modeling groundwater flow, predicting flood peaks and low flows, and river water quality. The exchange fluxes between river and aquifer are strongly influenced by the hydraulic conductivity of the riverbed which can vary several orders of magnitude and shows a strong spatial variation. Direct measurement of riverbed hydraulic conductivity is cumbersome and therefore often indirect data such as temperature data or calibration of groundwater models are used to constrain riverbed hydraulic conductivity. In these approaches, the riverbed is usually represented as a homogeneous geological structure and the spatial variation of riverbed hydraulic conductivity is thus neglected. However, neglecting this spatial variation can lead to systematic underestimation of net river-aquifer exchange fluxes and may have important implications for the estimation of peak mass flows, for the hydrochemistry of streambed sediments, nutrient cycling and biogeochemical gradients. The MODFLOW software is the most wide-spread package used for groundwater modelling. In MODFLOW rivers are usually modelled using the River-package. However, in this package no distinction can be made between horizontal and vertical riverbed hydraulic conductivity and the riverbed cannot be subdivided into layers with different hydraulic characteristics. Riverbed sediments are strongly layered and thus another approach is advised. Different ways of introducing heterogeneous riverbeds in MODFLOW groundwater flow models are explored and compared. The influence of heterogeneous riverbeds on groundwater-surface water exchange fluxes is analyzed for two case studies: the Aa River in the Nete catchment and a stretch of the Dijle River near the nature reserve 'de Doode Bemde' (Belgium). For both cases fine-scale distributed local groundwater flow

  5. Simulation of Ground-Water Flow and Effects of Ground-Water Irrigation on Base Flow in the Elkhorn and Loup River Basins, Nebraska

    USGS Publications Warehouse

    Peterson, Steven M.; Stanton, Jennifer S.; Saunders, Amanda T.; Bradley, Jesse R.

    2008-01-01

    Irrigated agriculture is vital to the livelihood of communities in the Elkhorn and Loup River Basins in Nebraska, and ground water is used to irrigate most of the cropland. Concerns about the sustainability of ground-water and surface-water resources have prompted State and regional agencies to evaluate the cumulative effects of ground-water irrigation in this area. To facilitate understanding of the effects of ground-water irrigation, a numerical computer model was developed to simulate ground-water flow and assess the effects of ground-water irrigation (including ground-water withdrawals, hereinafter referred to as pumpage, and enhanced recharge) on stream base flow. The study area covers approximately 30,800 square miles, and includes the Elkhorn River Basin upstream from Norfolk, Nebraska, and the Loup River Basin upstream from Columbus, Nebraska. The water-table aquifer consists of Quaternary-age sands and gravels and Tertiary-age silts, sands, and gravels. The simulation was constructed using one layer with 2-mile by 2-mile cell size. Simulations were constructed to represent the ground-water system before 1940 and from 1940 through 2005, and to simulate hypothetical conditions from 2006 through 2045 or 2055. The first simulation represents steady-state conditions of the system before anthropogenic effects, and then simulates the effects of early surface-water development activities and recharge of water leaking from canals during 1895 to 1940. The first simulation ends at 1940 because before that time, very little pumpage for irrigation occurred, but after that time it became increasingly commonplace. The pre-1940 simulation was calibrated against measured water levels and estimated long-term base flow, and the 1940 through 2005 simulation was calibrated against measured water-level changes and estimated long-term base flow. The calibrated 1940 through 2005 simulation was used as the basis for analyzing hypothetical scenarios to evaluate the effects of

  6. Modeling groundwater flow by lattice Boltzmann method in curvilinear coordinates

    NASA Astrophysics Data System (ADS)

    Budinski, Ljubomir; Fabian, Julius; Stipic, Matija

    2015-07-01

    In order to promote the use of the lattice Boltzmann method (LBM) for the simulation of isotropic groundwater flow in a confined aquifer with arbitrary geometry, Poisson's equation was transformed into a curvilinear coordinate system. With the metric function between the physical and the computational domain established, Poisson's equation written in Cartesian coordinates was transformed in curvilinear coordinates. Following, the appropriate equilibrium function for the D2Q9 square lattice has been defined. The resulting curvilinear formulation of the LBM for groundwater flow is capable of modeling flow in domains of complex geometry with the opportunity of local refining/coarsening of the computational mesh corresponding to the complexity of the flow pattern and the required accuracy. Since the proposed form of the LBM uses the transformed equation of flow implemented in the equilibrium function, finding a solution does not require supplementary procedures along the curvilinear boundaries, nor in the zones requiring mesh density adjustments. Thus, the basic concept of the LBM is completely maintained. The improvement of the proposed LBM over the previously published classical methods is completely verified by three examples with analytical solutions. The results demonstrate the advantages of the proposed curvilinear LBM in modeling groundwater flow in complex flow domains.

  7. Automatic Time Stepping with Global Error Control for Groundwater Flow Models

    SciTech Connect

    Tang, Guoping

    2008-09-01

    An automatic time stepping with global error control is proposed for the time integration of the diffusion equation to simulate groundwater flow in confined aquifers. The scheme is based on an a posteriori error estimate for the discontinuous Galerkin (dG) finite element methods. A stability factor is involved in the error estimate and it is used to adapt the time step and control the global temporal error for the backward difference method. The stability factor can be estimated by solving a dual problem. The stability factor is not sensitive to the accuracy of the dual solution and the overhead computational cost can be minimized by solving the dual problem using large time steps. Numerical experiments are conducted to show the application and the performance of the automatic time stepping scheme. Implementation of the scheme can lead to improvement in accuracy and efficiency for groundwater flow models.

  8. Comparison of groundwater flow in Southern California coastal aquifers

    USGS Publications Warehouse

    Hanson, Randall T.; Izbicki, John A.; Reichard, Eric G.; Edwards, Brian D.; Land, Michael; Martin, Peter

    2009-01-01

    Development of the coastal aquifer systems of Southern California has resulted in overdraft, changes in streamflow, seawater intrusion, land subsidence, increased vertical flow between aquifers, and a redirection of regional flow toward pumping centers. These water-management challenges can be more effectively addressed by incorporating new understanding of the geologic, hydrologic, and geochemical setting of these aquifers. Groundwater and surface-water flow are controlled, in part, by the geologic setting. The physiographic province and related tectonic fabric control the relation between the direction of geomorphic features and the flow of water. Geologic structures such as faults and folding also control the direction of flow and connectivity of groundwater flow. The layering of sediments and their structural association can also influence pathways of groundwater flow and seawater intrusion. Submarine canyons control the shortest potential flow paths that can result in seawater intrusion. The location and extent of offshore outcrops can also affect the flow of groundwater and the potential for seawater intrusion and land subsidence in coastal aquifer systems. As coastal aquifer systems are developed, the source and movement of ground-water and surface-water resources change. In particular, groundwater flow is affected by the relative contributions of different types of inflows and outflows, such as pump-age from multi-aquifer wells within basal or upper coarse-grained units, streamflow infiltration, and artificial recharge. These natural and anthropogenic inflows and outflows represent the supply and demand components of the water budgets of ground-water within coastal watersheds. They are all significantly controlled by climate variability related to major climate cycles, such as the El Niño–Southern Oscillation and the Pacific Decadal Oscillation. The combination of natural forcings and anthropogenic stresses redirects the flow of groundwater and either

  9. Perturbation of ground surface temperature reconstructions by groundwater flow?

    NASA Astrophysics Data System (ADS)

    Ferguson, Grant; Beltrami, Hugo; Woodbury, Allan D.

    2006-07-01

    Subsurface temperatures have been shown to be a robust source of information on past climates. Most analyses neglect groundwater flow (GWF) and assume purely conductive heat flow. However, in many situations GWF has not been fully considered and to date there are no general GWF criteria for either accepting or rejecting a temperature profile for paleoclimate analysis. Here we examine the transition from conduction dominated environments to environments where advection has a significant effect on the subsurface temperature regime and thus ground surface temperature (GST) histories. We provide guidelines indicating when advection is important and conclude that it is unlikely that groundwater flow is a significant source of error in the global data set maintained by the International Heat Flow Commission.

  10. Detecting groundwater flow direction from infrared thermal images

    NASA Astrophysics Data System (ADS)

    Ahmad, Akhundzadah Noor; Saito, Hirotaka; Asada, Kei; Kato, Makoto

    In this study we propose an approach to directly measure temperature changes around the heat source placed in a borehole using an infrared camera for determining groundwater flow direction. A plastic box filled with sands with a cylindrical perforated column placed in middle of the box was used in lab experiments. A heater was inserted in the center of the borehole. To detect the flow direction from the temperature distribution, an infrared camera that was placed face-down at the top of the cylinder (i.e., borehole) was used. COMSOL was used to numerically simulate coupled water flow and heat transfer to evaluate experimental results. Results show that when the flux is in the order of 10-2 to 10-4cm s-1, we can determine the groundwater flow direction because of a skewed temperature distribution due to convective transport of heat.

  11. Investigations of groundwater system and simulation of regional groundwater flow for North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

    USGS Publications Warehouse

    Senior, Lisa A.; Goode, Daniel J.

    2013-01-01

    Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations

  12. Approaches to the simulation of unconfined flow and perched groundwater flow in MODFLOW.

    PubMed

    Bedekar, Vivek; Niswonger, Richard G; Kipp, Kenneth; Panday, Sorab; Tonkin, Matthew

    2012-01-01

    Various approaches have been proposed to manage the nonlinearities associated with the unconfined flow equation and to simulate perched groundwater conditions using the MODFLOW family of codes. The approaches comprise a variety of numerical techniques to prevent dry cells from becoming inactive and to achieve a stable solution focused on formulations of the unconfined, partially-saturated, groundwater flow equation. Keeping dry cells active avoids a discontinuous head solution which in turn improves the effectiveness of parameter estimation software that relies on continuous derivatives. Most approaches implement an upstream weighting of intercell conductance and Newton-Raphson linearization to obtain robust convergence. In this study, several published approaches were implemented in a stepwise manner into MODFLOW for comparative analysis. First, a comparative analysis of the methods is presented using synthetic examples that create convergence issues or difficulty in handling perched conditions with the more common dry-cell simulation capabilities of MODFLOW. Next, a field-scale three-dimensional simulation is presented to examine the stability and performance of the discussed approaches in larger, practical, simulation settings. PMID:21635246

  13. Approaches to the simulation of unconfined flow and perched groundwater flow in MODFLOW

    USGS Publications Warehouse

    Bedekar, Vivek; Niswonger, Richard G.; Kipp, Kenneth; Panday, Sorab; Tonkin, Matthew

    2012-01-01

    Various approaches have been proposed to manage the nonlinearities associated with the unconfined flow equation and to simulate perched groundwater conditions using the MODFLOW family of codes. The approaches comprise a variety of numerical techniques to prevent dry cells from becoming inactive and to achieve a stable solution focused on formulations of the unconfined, partially-saturated, groundwater flow equation. Keeping dry cells active avoids a discontinuous head solution which in turn improves the effectiveness of parameter estimation software that relies on continuous derivatives. Most approaches implement an upstream weighting of intercell conductance and Newton-Raphson linearization to obtain robust convergence. In this study, several published approaches were implemented in a stepwise manner into MODFLOW for comparative analysis. First, a comparative analysis of the methods is presented using synthetic examples that create convergence issues or difficulty in handling perched conditions with the more common dry-cell simulation capabilities of MODFLOW. Next, a field-scale three-dimensional simulation is presented to examine the stability and performance of the discussed approaches in larger, practical, simulation settings.

  14. Application of Integral Pumping Tests to estimate the influence of losing streams on groundwater quality

    NASA Astrophysics Data System (ADS)

    Leschik, S.; Musolff, A.; Reinstorf, F.; Strauch, G.; Schirmer, M.

    2009-05-01

    Urban streams receive effluents of wastewater treatment plants and untreated wastewater during combined sewer overflow events. In the case of losing streams substances, which originate from wastewater, can reach the groundwater and deteriorate its quality. The estimation of mass flow rates Mex from losing streams to the groundwater is important to support groundwater management strategies, but is a challenging task. Variable inflow of wastewater with time-dependent concentrations of wastewater constituents causes a variable water composition in urban streams. Heterogeneities in the structure of the streambed and the connected aquifer lead, in combination with this variable water composition, to heterogeneous concentration patterns of wastewater constituents in the vicinity of urban streams. Groundwater investigation methods based on conventional point sampling may yield unreliable results under these conditions. Integral Pumping Tests (IPT) can overcome the problem of heterogeneous concentrations in an aquifer by increasing the sampled volume. Long-time pumping (several days) and simultaneous sampling yields reliable average concentrations Cav and mass flow rates Mcp for virtual control planes perpendicular to the natural flow direction. We applied the IPT method in order to estimate Mex of a stream section in Leipzig (Germany). The investigated stream is strongly influenced by combined sewer overflow events. Four pumping wells were installed up- and downstream of the stream section and operated for a period of five days. The study was focused on four inorganic (potassium, chloride, nitrate and sulfate) and two organic (caffeine and technical-nonylphenol) wastewater constituents with different transport properties. The obtained concentration-time series were used in combination with a numerical flow model to estimate Mcp of the respective wells. The difference of the Mcp's between up- and downstream wells yields Mex of wastewater constituents that increase

  15. Characterization and conceptualization of groundwater flow systems: Chapter 2

    USGS Publications Warehouse

    Plummer, L.N.; Sanford, W.E.; Glynn, P.D.

    2013-01-01

    This chapter discusses some of the fundamental concepts, data needs and approaches that aid in developing a general understanding of a groundwater system. Principles of the hydrological cycle are reviewed; the processes of recharge and discharge in aquifer systems; types of geological, hydrological and hydraulic data needed to describe the hydrogeological framework of an aquifer system; factors affecting the distribution of recharge to aquifers; and uses of groundwater chemistry, geochemical modelling, environmental tracers and age interpretations in groundwater studies. Together, these concepts and observations aid in developing a conceptualization of groundwater flow systems and provide input to the development of numerical models of a flow system. Conceptualization of the geology, hydrology, geochemistry, and hydrogeological and hydrochemical framework can be quite useful in planning, study design, guiding sampling campaigns, acquisition of new data and, ultimately, developing numerical models capable of assessing a wide variety of societal issues — for example, sustainability of groundwater resources in response to real or planned withdrawals from the system, CO2 sequestration or other waste isolation issues (such as nuclear waste disposal).

  16. A method to estimate groundwater depletion from confining layers

    USGS Publications Warehouse

    Konikow, L.F.; Neuzil, C.E.

    2007-01-01

    Although depletion of storage in low-permeability confining layers is the source of much of the groundwater produced from many confined aquifer systems, it is all too frequently overlooked or ignored. This makes effective management of groundwater resources difficult by masking how much water has been derived from storage and, in some cases, the total amount of water that has been extracted from an aquifer system. Analyzing confining layer storage is viewed as troublesome because of the additional computational burden and because the hydraulic properties of confining layers are poorly known. In this paper we propose a simplified method for computing estimates of confining layer depletion, as well as procedures for approximating confining layer hydraulic conductivity (K) and specific storage (Ss) using geologic information. The latter makes the technique useful in developing countries and other settings where minimal data are available or when scoping calculations are needed. As such, our approach may be helpful for estimating the global transfer of groundwater to surface water. A test of the method on a synthetic system suggests that the computational errors will generally be small. Larger errors will probably result from inaccuracy in confining layer property estimates, but these may be no greater than errors in more sophisticated analyses. The technique is demonstrated by application to two aquifer systems: the Dakota artesian aquifer system in South Dakota and the coastal plain aquifer system in Virginia. In both cases, depletion from confining layers was substantially larger than depletion from the aquifers.

  17. Mathematical modelling of surface water-groundwater flow and salinity interactions in the coastal zone

    NASA Astrophysics Data System (ADS)

    Spanoudaki, Katerina; Kampanis, Nikolaos A.

    2014-05-01

    surface water-groundwater model IRENE (Spanoudaki et al., 2009; Spanoudaki, 2010) has been modified in order to simulate surface water-groundwater flow and salinity interactions in the coastal zone. IRENE, in its original form, couples the 3D, non-steady state Navier-Stokes equations, after Reynolds averaging and with the assumption of hydrostatic pressure distribution, to the equations describing 3D saturated groundwater flow of constant density. A semi-implicit finite difference scheme is used to solve the surface water flow equations, while a fully implicit finite difference scheme is used for the groundwater equations. Pollution interactions are simulated by coupling the advection-diffusion equation describing the fate and transport of contaminants introduced in a 3D turbulent flow field to the partial differential equation describing the fate and transport of contaminants in 3D transient groundwater flow systems. The model has been further developed to include the effects of density variations on surface water and groundwater flow, while the already built-in solute transport capabilities are used to simulate salinity interactions. Initial results show that IRENE can accurately predict surface water-groundwater flow and salinity interactions in coastal areas. Important research issues that can be investigated using IRENE include: (a) sea level rise and tidal effects on aquifer salinisation and the configuration of the saltwater wedge, (b) the effects of surface water-groundwater interaction on salinity increase of coastal wetlands and (c) the estimation of the location and magnitude of groundwater discharge to coasts. Acknowledgement The work presented in this paper has been funded by the Greek State Scholarships Foundation (IKY), Fellowships of Excellence for Postdoctoral Studies (Siemens Program), 'A simulation-optimization model for assessing the best practices for the protection of surface water and groundwater in the coastal zone', (2013 - 2015). References

  18. MODIS-aided statewide net groundwater-recharge estimation in Nebraska.

    PubMed

    Szilagyi, Jozsef; Jozsa, Janos

    2013-01-01

    Monthly evapotranspiration (ET) rates (2000 to 2009) across Nebraska at about 1-km resolution were obtained by linear transformations of the MODIS (MODerate resolution Imaging Spectroradiometer) daytime surface temperature values with the help of the Priestley-Taylor equation and the complementary relationship of evaporation. For positive values of the mean annual precipitation and ET differences, the mean annual net recharge was found by an additional multiplication of the power-function-transformed groundwater vulnerability DRASTIC-code values. Statewide mean annual net recharge became about 29 mm (i.e., 5% of mean annual precipitation) with the largest recharge rates (in excess of 100 mm/year) found in the eastern Sand Hills and eastern Nebraska. Areas with the largest negative net recharge rates caused by declining groundwater levels due to large-scale irrigation are found in the south-western region of the state. Error bounds of the estimated values are within 10% to 15% of the corresponding precipitation rates and the estimated net recharge rates are sensitive to errors in the precipitation and ET values. This study largely confirms earlier base-flow analysis-based statewide groundwater recharge estimates when considerations are made for differences in the recharge definitions. The current approach not only provides better spatial resolution than available earlier studies for the region but also quantifies negative net recharge rates that become especially important in numerical modeling of shallow groundwater systems. PMID:23216050

  19. Geohydrology, simulation of regional groundwater flow, and assessment of watermanagement strategies, Twentynine Palms area, California

    USGS Publications Warehouse

    Li, Zhen; Martin, Peter

    2011-01-01

    aquifers (referred to as the upper and the middle aquifers) and the Tertiary sedimentary deposits into a single aquifer (referred to as the lower aquifer). In general, wells perforated in the upper aquifer yield more water than wells perforated in the middle and lower aquifers. The study area is dominated by extensive faulting and moderate to intense folding that has displaced or deformed the pre-Tertiary basement complex as well as the overlying Tertiary and Quaternary deposits. Many of these faults act as barriers to the lateral movement of groundwater flow and form many of the boundaries of the groundwater subbasins. The principal recharge to the study area is groundwater underflow across the western and southern boundaries that originates as runoff in the surrounding mountains. Groundwater discharges naturally from the study area as spring flow, as groundwater underflow to downstream basins, and as water vapor to the atmosphere by transpiration of phreatophytes and direct evaporation from moist soil. The annual volume of water that naturally recharged to or discharged from the groundwater flow system in the study area during predevelopment conditions was estimated to be 1,010 acre-feet per year (acre-ft/yr). About 90 percent of this recharge originated as runoff from the Little San Bernardino and the Pinto Mountains to the south, and the remainder originated as runoff from the San Bernardino Mountains to the west. Evapotranspiration by phreatophytes near Mesquite Lake (dry) was the primary form of predevelopment groundwater discharge. From 1953 through 2007, approximately 139,400 acre-feet (acre-ft) of groundwater was pumped by the MCAGCC from the Surprise Spring subbasin. A regional-scale numerical groundwater flow model was developed using MODFLOW-2000 for the Surprise Spring, Deadman, Mesquite, and Mainside subbasins. The aquifer system was simulated by using three model layers representing the upper, middle, and lower aquifers. Measured groundwater levels

  20. Simulation of groundwater flow around a pilot waste tire site

    SciTech Connect

    Lin, Y.; Chyi, L.L. . Dept. of Geology)

    1994-04-01

    The groundwater system in the area contains several aquifers divided by underclays as confining layers. The three aquifers addressed under this research include the upper aquifer above the No. 6 underclay, the middle aquifer above the No. 5 underclay and the lower aquifer above the No. 4 underclay. A quasi-3D conceptual model is used in simulating groundwater flow and MODFLOW is used for the actual modeling. The boundary of the model area is defined by a local divide on the north, the Hugle Run on the west and the Temple Ditch on the east and the Sandy Creek on the south. The model consists of three aquifer layers and two confining layers. The upper one is unconfined, the middle one is a hybrid, and the lower one is confined. In the simulation of the upper aquifer, drain boundary condition is defined to include seepage. In the simulation of the middle aquifer, the Hugle Run and the Temple ditch associated with the alluvium represent a drain boundary and the Sand Creek to the south represents a constant head. The lower aquifer does not have a natural hydrogeologic boundary. The equipotential contours are assigned as constant head boundary. Simulation has indicated that minimal groundwater flow existed in the upper aquifer and it is flowing toward the southwest. Groundwater flow in the middle aquifer near the pilot waste site is somewhat divergent to the south. The velocity is slow near the site but increases quickly toward the southwest. Leakage from the upper aquifer appears to be predominant over the regional flow in the recharge of the middle aquifer. The lower aquifer is recharged mainly by the regional flow and is flowing almost due south.

  1. The transition of flow patterns through critical stagnation points in two-dimensional groundwater flow

    Technology Transfer Automated Retrieval System (TEKTRAN)

    A flow pattern is characterized by aquifer features and the number, type, and distribution of stagnation points (locations where the discharge is zero). This article identifies a condition for transition of flow patterns in two-dimensional groundwater flow obeying Darcy's law by examining changes in...

  2. Bias in groundwater samples caused by wellbore flow

    USGS Publications Warehouse

    Reilly, Thomas E.; Franke, O. Lehn; Bennett, Gordon D.

    1989-01-01

    Proper design of physical installations and sampling procedures for groundwater monitoring networks is critical for the detection and analysis of possible contaminants. Monitoring networks associated with known contaminant sources sometimes include an array of monitoring wells with long well screens. The purpose of this paper is: (a) to report the results of a numerical experiment indicating that significant borehole flow can occur within long well screens installed in homogeneous aquifers with very small head differences in the aquifer (less than 0.01 feet between the top and bottom of the screen); (b) to demonstrate that contaminant monitoring wells with long screens may completely fail to fulfill their purpose in many groundwater environments.

  3. Improving Radium-based Estimates of Submarine Groundwater Discharge

    NASA Astrophysics Data System (ADS)

    Hughes, A. L.; Wilson, A. M.

    2011-12-01

    Groundwater discharge is vital for the exchange of solutes between salt marshes and estuaries, and radium isotopes are frequently used as tracers of groundwater flow paths and discharge in coastal systems. Considerable spatial and temporal variability in porewater radium activity has hindered the accuracy of this tracer. In porewater, radium activity is a complex function of production by parent isotopes in and grain size of the aquifer material, individual decay rates, porewater salinity, temperature, redox- and pH-dependent adsorption and desorption, sediment Fe- and Mn-oxide/hydroxide coatings, and groundwater transport (advection and dispersion). In order to resolve the primary factors controlling porewater radium activity in an intertidal salt marsh, where high salinity and reducing conditions prevail, and sediment oxide coatings vary from winter to summer, a field and modeling study was conducted at a salt marsh island within North Inlet Salt Marsh, Georgetown, South Carolina. This site was previously developed as part of a larger study to understand the links between salt marsh groundwater dynamics and acute marsh dieback. Porewater and surface water samples were collected from November 2009 - February 2011. Shallow sediment samples were collected in winter and summer 2010, and deeper sediments were split from cores collected during site development. Measurements of water temperature, salinity, mV, and pH were taken in the field, and radium isotopes were measured via delayed-coincidence counter or gamma spectrometry. Surface-bound sediment radium activity was determined by desorption experiments. Iron and manganese oxide coatings on surface sediments were isolated through a sequential leaching process, and the leachate analyzed via ICP-AES. Finally, a 3-D groundwater flow model was developed using SUTRA, a U.S.G.S. numerical model, which was modified to account for changes in total stress resulting from tidal loading of the marsh surface and for complex

  4. Groundwater availability as constrained by hydrogeology and environmental flows

    USGS Publications Warehouse

    Watson, Katelyn A.; Mayer, Alex S.; Reeves, Howard W.

    2014-01-01

    Groundwater pumping from aquifers in hydraulic connection with nearby streams has the potential to cause adverse impacts by decreasing flows to levels below those necessary to maintain aquatic ecosystems. The recent passage of the Great Lakes-St. Lawrence River Basin Water Resources Compact has brought attention to this issue in the Great Lakes region. In particular, the legislation requires the Great Lakes states to enact measures for limiting water withdrawals that can cause adverse ecosystem impacts. This study explores how both hydrogeologic and environmental flow limitations may constrain groundwater availability in the Great Lakes Basin. A methodology for calculating maximum allowable pumping rates is presented. Groundwater availability across the basin may be constrained by a combination of hydrogeologic yield and environmental flow limitations varying over both local and regional scales. The results are sensitive to factors such as pumping time, regional and local hydrogeology, streambed conductance, and streamflow depletion limits. Understanding how these restrictions constrain groundwater usage and which hydrogeologic characteristics and spatial variables have the most influence on potential streamflow depletions has important water resources policy and management implications.

  5. Budgets and chemical characterization of groundwater for the Diamond Valley flow system, central Nevada, 2011–12

    USGS Publications Warehouse

    Berger, David L.; Mayers, C. Justin; Garcia, C. Amanda; Buto, Susan G.; Huntington, Jena M.

    2016-01-01

    The pre-development, steady state, groundwater budget for the Diamond Valley flow system was estimated at about 70,000 acre-ft/yr of inflow and outflow. During years 2011–12, inflow components of groundwater recharge from precipitation and subsurface inflow from adjacent basins totaled 70,000 acre-ft/yr for the DVFS, whereas outflow components included 64,000 acre-ft/yr of groundwater evapotranspiration and 69,000 acre-ft/yr of net groundwater withdrawals, or net pumpage. Spring discharge in northern Diamond Valley declined about 6,000 acre-ft/yr between pre-development time and years 2011–12. Assuming net groundwater withdrawals minus spring flow decline is equivalent to the storage change, the 2011–12 summation of inflow and storage change was balanced with outflow at about 133,000 acre-ft/yr.

  6. Partitioning a regional groundwater flow system into shallow local and deep regional flow compartments

    NASA Astrophysics Data System (ADS)

    Goderniaux, Pascal; Davy, Philippe; Bresciani, Etienne; Dreuzy, Jean-Raynald; Borgne, Tanguy

    2013-04-01

    The distribution of groundwater fluxes in aquifers is strongly influenced by topography, and organized between hillslope and regional scales. The objective of this study is to provide new insights regarding the compartmentalization of aquifers at the regional scale and the partitioning of recharge between shallow/local and deep/regional groundwater transfers. A finite-difference flow model was implemented, and the flow structure was analyzed as a function of recharge (from 20 to 500 mm/yr), at the regional-scale (1400 km2), in three dimensions, and accounting for variable groundwater discharge zones; aspects which are usually not considered simultaneously in previous studies. The model allows visualizing 3-D circulations, as those provided by Tothian models in 2-D, and shows local and regional transfers, with 3-D effects. The probability density function of transit times clearly shows two different parts, interpreted using a two-compartment model, and related to regional groundwater transfers and local groundwater transfers. The role of recharge on the size and nature of the flow regimes, including groundwater pathways, transit time distributions, and volumes associated to the two compartments, have been investigated. Results show that topography control on the water table and groundwater compartmentalization varies with the recharge rate applied. When recharge decreases, the absolute value of flow associated to the regional compartment decreases, whereas its relative value increases. The volume associated to the regional compartment is calculated from the exponential part of the two-compartment model, and is nearly insensitive to the total recharge fluctuations.

  7. Estimated ground-water discharge by evapotranspiration, Ash Meadows Area, Nye County, Nevada, 1994

    SciTech Connect

    Nichols, W.D.; Laczniak, R.J.; DeMeo, G.A.; Rapp, T.R.

    1997-05-01

    Ground water discharges from the regional ground-water flow system that underlies the eastern part of the Nevada Test Site through numerous springs and seeps in the Ash Meadows National Wildlife Refuge in southern Nevada. The total spring discharge was estimated to be about 17,000 acre-feet per year by earlier studies. Previous studies estimated that about 10,500 acre-feet of this discharge was lost to evapotranspiration. The present study was undertaken to develop a more rigorous approach to estimating ground-water discharge in the Ash Meadows area. Part of the study involves detailed field investigation of evapotranspiration. Data collection began in early 1994. The results of the first year of study provide a basis for making preliminary estimates of ground-water discharge by evapotranspiration. An estimated 13,100 acre-feet of ground water was evapotranspired from about 6,800 acres of marsh and salt-grass. Additional 3,500 acre-feet may have been transpired from the open water and from about 1,460 acres of other areas of Ash Meadows in which field studies have not yet been made.

  8. Analysis of ground-water flow along a regional flow path of the Midwestern Basins and Arches aquifer system in Ohio

    USGS Publications Warehouse

    Hanover, R.H.

    1994-01-01

    A cross-sectional analysis of ground-water flow in central-western and northwestern Ohio was done as part of the Midwestern Basins and Arches Regional Aquifer-System Analysis project. The Midwestern Basins and Arches aquifer system is composed of carbonate bedrock of Silurian and Devonian age and overlying glacial flow analysis of the Scioto and Blanchard rivers in the study area were used to describe patrems of ground-water flow, to evaluate stream-aquifer interaction, and to quantify recharge and discharge within the ground-water-flow system along a regional ground-water-flow path. The selected regional flow path begins at a regional topographic high in Logan County, Ohio, and ends in Sandusky Bay (Lake Erie), a regional topographic low. Recharge to the ground-water system along the selected regional flow path was estimated from hydrograph separation of streamflow and averaged 3.24 inches per year. Computer model simulations indicate that 84 percent of the water entering the ground-water system flows less than 5 miles from point of recharge to point of discharge and no deeper than the upper surficial aquifers. The distance and depth that ground water travels and traveltime from point of recharge to point of discharge is controlled largely by where ground water enters the flow system. Ground water entering the flow system in the vicinity of major surface- water divides generally travels further, deeper, and longer than ground water entering the flow system elsewhere along the regional flow path. Particle tracking simulations substantiate the concept that the 80-mile-long regional flow path is within a continuous ground-water basin. Estimated traveltimes for ground-water from the regional high to Sandusky Bay range from 22,000 to 40,700 years, given a range of porosities from 8 to 22 percent for the carbonate-rock aquifer.

  9. Global estimates of submarine groundwater discharge using numerical modeling and geomatics

    NASA Astrophysics Data System (ADS)

    Luijendijk, E.; Gleeson, T.; Ferguson, G. A.

    2011-12-01

    Submarine groundwater discharge (SGD), the flow of fresh or saline groundwater to an ocean, may be a significant contributor to the water and chemical budgets of the world oceans. SGD consists of fresh, terrestrial groundwater driven by hydraulic gradients, the focus of this research, and re-circulated seawater driven by tidal pumping, wave set-up, convection and hydraulic gradients. We couple density-dependent analytical and numerical simulations of generic models of coastal topography and geology with geomatic data bases to resolve the rate and driving mechanisms of terrestrially-derived submarine groundwater discharge globally. Two analytical models lead to linear relationships between SGD and the key predictive parameters: hydraulic gradients, hydraulic conductivities, aquifer thickness and recharge. Average global geomatic parameters suggest global SGD ranges from 0.01% to 0.2% of global river run off which much lower than most previous estimates of global SGD. Quantifying submarine groundwater discharge is critical because SGD is a poorly constrained flux that can significantly contribute to eutrophication or water quality decline in coastal areas.

  10. Estimation of methane concentrations and loads in groundwater discharge to Sugar Run, Lycoming County, Pennsylvania

    USGS Publications Warehouse

    Heilweil, Victor M.; Risser, Dennis W.; Conger, Randall W.; Grieve, Paul L.; Hynek, Scott A.

    2014-01-01

    A stream-sampling study was conducted to estimate methane concentrations and loads in groundwater discharge to a small stream in an active shale-gas development area of northeastern Pennsylvania. Grab samples collected from 15 streams in Bradford, Lycoming, Susquehanna, and Tioga Counties, Pa., during a reconnaissance survey in May and June 2013 contained dissolved methane concentrations ranging from less than the minimum reporting limit (1.0) to 68.5 micrograms per liter (µg/L). The stream-reach mass-balance method of estimating concentrations and loads of methane in groundwater discharge was applied to a 4-kilometer (km) reach of Sugar Run in Lycoming County, one of the four streams with methane concentrations greater than or equal to 5 µg/L. Three synoptic surveys of stream discharge and methane concentrations were conducted during base-flow periods in May, June, and November 2013. Stream discharge at the lower end of the reach was about 0.10, 0.04, and 0.02 cubic meters per second, respectively, and peak stream methane concentrations were about 20, 67, and 29 µg/L. In order to refine estimated amounts of groundwater discharge and locations where groundwater with methane discharges to the stream, the lower part of the study reach was targeted more precisely during the successive studies, with approximate spacing between stream sampling sites of 800 meters (m), 400 m, and 200 m, in May, June, and November, respectively. Samples collected from shallow piezometers and a seep near the location of the peak methane concentration measured in streamwater had groundwater methane concentrations of 2,300 to 4,600 µg/L. These field data, combined with one-dimensional stream-methane transport modeling, indicate groundwater methane loads of 1.8 ±0.8, 0.7 ±0.3, and 0.7 ±0.2 kilograms per day, respectively, discharging to Sugar Run. Estimated groundwater methane concentrations, based on the transport modeling, ranged from 100 to 3,200 µg/L. Although total methane load

  11. Use of a ground-water flow model with particle tracking to evaluate ground-water vulnerability, Clark County, Washington

    USGS Publications Warehouse

    Snyder, D.T.; Wilkinson, J.M.; Orzol, L.L.

    1996-01-01

    public-supply wells in Clark County may be receiving a component of water that recharged in areas that are more conducive to contaminant entry. The aquifer sensitivity maps illustrate a critical deficiency in the DRASTIC methodology: the failure to account for the dynamics of the ground-water flow system. DRASTIC indices calculated for a particular location thus do not necessarily reflect the conditions of the ground-water resources at the recharge areas to that particular location. Each hydrogeologic unit was also mapped to highlight those areas that will eventually receive flow from recharge areas with on-site waste-disposal systems. Most public-supply wells in southern Clark County may eventually receive a component of water that was recharged from on-site waste-disposal systems.Traveltimes from particle tracking were used to estimate the minimum and maximum age of ground water within each model-grid cell. Chlorofluorocarbon (CFC)-age dating of ground water from 51 wells was used to calibrate effective porosity values used for the particle- tracking program by comparison of ground-water ages determined through the use of the CFC-age dating with those calculated by the particle- tracking program. There was a 76 percent agreement in predicting the presence of modern water in the 51 wells as determined using CFCs and calculated by the particle-tracking program. Maps showing the age of ground water were prepared for all the hydrogeologic units. Areas with the youngest ground-water ages are expected to be at greatest risk for contamination from anthropogenic sources. Comparison of these maps with maps of public- supply wells in Clark County indicates that most of these wells may withdraw ground water that is, in part, less than 100 years old, and in many instances less than 10 years old. Results of the analysis showed that a single particle-tracking analysis simulating advective transport can be used to evaluate ground-water vulnerability for any part of a ground-wate

  12. Closing the irrigation deficit in Cambodia: Implications for transboundary impacts on groundwater and Mekong River flow

    NASA Astrophysics Data System (ADS)

    Erban, Laura E.; Gorelick, Steven M.

    2016-04-01

    Rice production in Cambodia, essential to food security and exports, is largely limited to the wet season. The vast majority (96%) of land planted with rice during the wet season remains fallow during the dry season. This is in large part due to lack of irrigation capacity, increases in which would entail significant consequences for Cambodia and Vietnam, located downstream on the Mekong River. Here we quantify the extent of the dry season "deficit" area in the Cambodian Mekong River catchment, using a recent agricultural survey and our analysis of MODIS satellite data. Irrigation of this land for rice production would require a volume of water up to 31% of dry season Mekong River flow to Vietnam. However, the two countries share an aquifer system in the Mekong Delta, where irrigation demand is increasingly met by groundwater. We estimate expansion rates of groundwater-irrigated land to be >10% per year in the Cambodian Delta using LANDSAT satellite data and simulate the effects of future expansion on groundwater levels over a 25-year period. If groundwater irrigation continues to expand at current rates, the water table will drop below the lift limit of suction pump wells, used for domestic supply by >1.5 million people, throughout much of the area within 15 years. Extensive groundwater irrigation jeopardizes access for shallow domestic water supply wells, raises the costs of pumping for all groundwater users, and may exacerbate arsenic contamination and land subsidence that are already widespread hazards in the region.

  13. Statistical Performance Evaluation of Spatiotemporal Characteristics of Groundwater Flow and Contaminant Mass Transport

    NASA Astrophysics Data System (ADS)

    Matiatos, Ioannis; Papadopoulou, Maria P.; Varouchakis, Emmanouil A.

    2016-04-01

    As groundwater remains one of the most critical natural resources worldwide, numerical models of groundwater flow and contaminant mass transport provide a reliable tool for the efficient protection, planning and sustainable management of groundwater resources. This work focuses on the evaluation of the performance of different numerical models which have been developed to simulate spatiotemporal groundwater flow and contaminant mass transport in a coastal aquifer system. The evaluation of the models' performance has been based on 9 different statistical measures and indices of goodness of fit. Overall, the simulation of groundwater level and contaminant mass concentration delivered very good calibration and validation results in all cases, quite close to the desired values. Maps of aquifer water level and contaminant mass concentrations are provided for all cases in order the differences to be discussed and assessed. The selection of the appropriate model(s) is case oriented and it should be based on the problem's characteristics in order the spatiotemporal variability of the components under study to be optimally estimated.

  14. Comparison of groundwater flow model particle tracking results and isotopic data in the Leon valley, Mexico

    NASA Astrophysics Data System (ADS)

    Hernandez-Garcia, G.; Cortes, A. S.; Martínez-Reyes, J.; Perez-Quezadas, J.; Grupo de Hidrologia Isotopica

    2013-05-01

    The study area is located in the Guanajuato state, north-west of Mexico City. The Leon Valley covers with groundwater its water demand estimated in about 20.6 m3/s. The constant population increase and related economic activities in the region have a steady growth in water needs. Related abstraction rate has produced an average drawdown of about 1.0 m/year in the last two decades. It suggests that present groundwater management needs to be reviewed. The groundwater management in the study area implies a possibility that abstraction will produce environmental impacts. This vital resource under stress becomes necessary to study its hydro-geologic functioning to achieve a scientific groundwater management in the valley. This investigation was based on the analysis and integration of existing information and the one generated in the field by the authors. Highlighted concepts were: i) the geologic structure of the area, ii) the hydraulic parameters and iii) the delta-deuterium and delta-oxigen-18 composition. This information was analysed integrally by means of applying a groundwater flow model (MODFLOW) and a particle-tracking model (FLOWPATH): the results were similar to flow paths and time-of travel interpretations derived from isotopic data.

  15. Multivariate analyses with end-member mixing to characterize groundwater flow: Wind Cave and associated aquifers

    USGS Publications Warehouse

    Long, A.J.; Valder, J.F.

    2011-01-01

    Principal component analysis (PCA) applied to hydrochemical data has been used with end-member mixing to characterize groundwater flow to a limited extent, but aspects of this approach are unresolved. Previous similar approaches typically have assumed that the extreme-value samples identified by PCA represent end members. The method presented herein is different from previous work in that (1) end members were not assumed to have been sampled but rather were estimated and constrained by prior knowledge; (2) end-member mixing was quantified in relation to hydrogeologic domains, which focuses model results on major hydrologic processes; (3) a method to select an appropriate number of end members using a series of cluster analyses is presented; and (4) conservative tracers were weighted preferentially in model calibration, which distributed model errors of optimized values, or residuals, more appropriately than would otherwise be the case. The latter item also provides an estimate of the relative influence of geochemical evolution along flow paths in comparison to mixing. This method was applied to groundwater in Wind Cave and the associated karst aquifer in the Black Hills of South Dakota, USA. The end-member mixing model was used to test a hypothesis that five different end-member waters are mixed in the groundwater system comprising five hydrogeologic domains. The model estimated that Wind Cave received most of its groundwater inflow from local surface recharge with an additional 33% from an upgradient aquifer. Artesian springs in the vicinity of Wind Cave primarily received water from regional groundwater flow. ?? 2011.

  16. Heat flow and subsurface temperature as evidence for basin-scale ground-water flow, North Slope of Alaska

    USGS Publications Warehouse

    Deming, D.; Sass, J.H.; Lachenbruch, A.H.; De Rito, R. F.

    1992-01-01

    Several high-resolution temperature logs were made in each of 21 drillholes and a total of 601 thermal conductivity measurements were made on drill cuttings and cores. Near-surface heat flow (??20%) is inversely correlated with elevation and ranges from a low of 27 mW/m2 in the foothills of the Brooks Range in the south, to a high of 90 mW/m2 near the north coast. Subsurface temperatures and thermal gradients estimated from corrected BHTs are similarly much higher on the coastal plain than in the foothills province to the south. Significant east-west variation in heat flow and subsurface temperature is also observed; higher heat flow and temperature coincide with higher basement topography. The observed thermal pattern is consistent with forced convection by a topographically driven ground-water flow system. Average ground-water (Darcy) velocity in the postulated flow system is estimated to be of the order of 0.1 m/yr; the effective basin-scale permeability is estimated to be of the order of 10-14 m2. -from Authors

  17. Development of a numerical model to simulate groundwater flow in the shallow aquifer system of Assateague Island, Maryland and Virginia

    USGS Publications Warehouse

    Masterson, John P.; Fienen, Michael N.; Gesch, Dean B.; Carlson, Carl S.

    2013-01-01

    A three-dimensional groundwater-flow model was developed for Assateague Island in eastern Maryland and Virginia to simulate both groundwater flow and solute (salt) transport to evaluate the groundwater system response to sea-level rise. The model was constructed using geologic and spatial information to represent the island geometry, boundaries, and physical properties and was calibrated using an inverse modeling parameter-estimation technique. An initial transient solute-transport simulation was used to establish the freshwater-saltwater boundary for a final calibrated steady-state model of groundwater flow. This model was developed as part of an ongoing investigation by the U.S. Geological Survey Climate and Land Use Change Research and Development Program to improve capabilities for predicting potential climate-change effects and provide the necessary tools for adaptation and mitigation of potentially adverse impacts.

  18. Numerical simulation of groundwater flow in Dar es Salaam Coastal Plain (Tanzania)

    NASA Astrophysics Data System (ADS)

    Luciani, Giulia; Sappa, Giuseppe; Cella, Antonella

    2016-04-01

    They are presented the results of a groundwater modeling study on the Coastal Aquifer of Dar es Salaam (Tanzania). Dar es Salaam is one of the fastest-growing coastal cities in Sub-Saharan Africa, with with more than 4 million of inhabitants and a population growth rate of about 8 per cent per year. The city faces periodic water shortages, due to the lack of an adequate water supply network. These two factors have determined, in the last ten years, an increasing demand of groundwater exploitation, carried on by quite a number of private wells, which have been drilled to satisfy human demand. A steady-state three dimensional groundwater model has been set up by the MODFLOW code, and calibrated with the UCODE code for inverse modeling. The aim of the model was to carry out a characterization of groundwater flow system in the Dar es Salaam Coastal Plain. The inputs applied to the model included net recharge rate, calculated from time series of precipitation data (1961-2012), estimations of average groundwater extraction, and estimations of groundwater recharge, coming from zones, outside the area under study. Parametrization of the hydraulic conductivities was realized referring to the main geological features of the study area, based on available literature data and information. Boundary conditions were assigned based on hydrogeological boundaries. The conceptual model was defined in subsequent steps, which added some hydrogeological features and excluded other ones. Calibration was performed with UCODE 2014, using 76 measures of hydraulic head, taken in 2012 referred to the same season. Data were weighted on the basis of the expected errors. Sensitivity analysis of data was performed during calibration, and permitted to identify which parameters were possible to be estimated, and which data could support parameters estimation. Calibration was evaluated based on statistical index, maps of error distribution and test of independence of residuals. Further model

  19. Estimating Natural Flows into the California's Sacramento - San Joaquin Delta

    NASA Astrophysics Data System (ADS)

    Huang, G.; Kadir, T.; Chung, F. I.

    2014-12-01

    Natural flows into the California's Sacramento - San Joaquin Delta under predevelopment vegetative conditions, if and when reconstructed, can serve as a useful guide to establish minimum stream flows, restoration targets, and a basis for assessing impacts of global warming in the Bay-Delta System. Daily simulations of natural Delta flows for the period 1922-2009 were obtained using precipitation-snowmelt-runoff models for the upper watersheds that are tributaries to the California's Central Valley, and then routing the water through the Central Valley floor area using a modified version of the California Central Valley Groundwater-Surface Water Simulation Model (C2VSIM) for water years 1922 through 2009. Daily stream inflows from all major upper watersheds were simulated using 23 Soil Water Assessment Tool (SWAT) models. Historical precipitation and reference evapotranspiration data were extracted from the SIMETAW2 with the 4km gridded meteorological data. The Historical natural and riparian vegetation distributions were compiled from several pre-1900 historical vegetation maps of the Central Valley. Wetlands were dynamically simulated using interconnected lakes. Flows overtopping natural levees were simulated using flow rating curves. New estimates of potential evapotranspiration from different vegetative classes under natural conditions were also used. Sensitivity simulations demonstrate that evapotranspiration estimates, native vegetation distribution, surface-groundwater interaction parameters, extinction depth for groundwater uptake, and other physical processes play a key role in the magnitude and timing of upstream flows arriving at the Delta. Findings contradict a common misconception that the magnitude of inflows to the Delta under natural vegetative conditions is greater than those under the historical agricultural and urban land use development. The developed models also enable to study the impacts of global warming by modifying meteorological and

  20. Recharge and Lateral Groundwater Flow Boundary Conditions for the Saturated Zone Site-Scale Flow and Transport Model

    SciTech Connect

    S. James

    2004-10-06

    This analysis is designed to use existing modeling and analysis results as the basis for estimated groundwater flow rates into the saturated zone (SZ) site-scale model domains, both as recharge (infiltration) at the upper boundary (water table), and as underflow at the lateral boundaries. Specifically, this work compiles information on the recharge boundary conditions supplied to the base-case and alternate SZ site-scale flow models taken from (1) distributed recharge from the 1997 (D'Agnese et al. 1997 [DIRS 100131]) or 2001 (D'Agnese et al. 2002 [DIRS 158876]) SZ regional-scale (Death Valley Regional Flow System [DVRFS]) model; (2) recharge below the area of the 1997 (Wu et al. 1997 [DIRS 156453]) or 2003 (BSC 2004 [DIRS 169861]) unsaturated zone (UZ) site-scale flow model; and (3) focused recharge along Fortymile Wash. In addition, this analysis includes extraction of the groundwater flow rates simulated by the 1997 and 2001 DVRFS models coincident with the lateral boundaries of the SZ site-scale flow models. The fluxes from the 1997 DVRFS were used to calibrate the base-case SZ site-scale flow model. The 2001 DVRFS fluxes are used in the alternate SZ site-scale flow model.

  1. Flow partitioning in regional groundwater flow systems as a function of recharge and topography

    NASA Astrophysics Data System (ADS)

    Goderniaux, P.; Davy, P.; Bresciani, E.; De Dreuzy, J.; Le Borgne, T.

    2013-12-01

    The distribution of groundwater fluxes in aquifers is strongly influenced by topography, and organized between hillslope and regional scales. In this study, we use a finite-difference flow model to quantify the partitioning of recharge and compartmentalization of aquifers between shallow/local and deep/regional groundwater transfers. The flow structure is analyzed for a regional aquifers, as a function of recharge (from 20 to 500 mm/yr), in 3-dimensions, and accounting for variable groundwater discharge zones. The Probability Density Function of transit times shows two different parts, interpreted using a two-compartment model, related to regional and local groundwater flows. The role of recharge on the size and nature of the flow regimes, including groundwater pathways, transit time distributions, and volumes associated to the two compartments is investigated. Results show that topography control on the water table and groundwater compartmentalization varies with the recharge rate applied. The volume associated to the regional compartment is calculated from the exponential part of the two-compartment model, and is nearly insensitive to the total recharge fluctuations. The model also allows visualizing 3D circulations, as those provided by Tothian models in 2D, and shows local and regional transfers, with 3D effects. Results are presented for a specific basin (1400 km2) in Brittany (France). Preliminary results using different kinds of topography are presented and compared.

  2. Parameter Estimation for Groundwater Models under Uncertain Irrigation Data.

    PubMed

    Demissie, Yonas; Valocchi, Albert; Cai, Ximing; Brozovic, Nicholas; Senay, Gabriel; Gebremichael, Mekonnen

    2015-01-01

    The success of modeling groundwater is strongly influenced by the accuracy of the model parameters that are used to characterize the subsurface system. However, the presence of uncertainty and possibly bias in groundwater model source/sink terms may lead to biased estimates of model parameters and model predictions when the standard regression-based inverse modeling techniques are used. This study first quantifies the levels of bias in groundwater model parameters and predictions due to the presence of errors in irrigation data. Then, a new inverse modeling technique called input uncertainty weighted least-squares (IUWLS) is presented for unbiased estimation of the parameters when pumping and other source/sink data are uncertain. The approach uses the concept of generalized least-squares method with the weight of the objective function depending on the level of pumping uncertainty and iteratively adjusted during the parameter optimization process. We have conducted both analytical and numerical experiments, using irrigation pumping data from the Republican River Basin in Nebraska, to evaluate the performance of ordinary least-squares (OLS) and IUWLS calibration methods under different levels of uncertainty of irrigation data and calibration conditions. The result from the OLS method shows the presence of statistically significant (p < 0.05) bias in estimated parameters and model predictions that persist despite calibrating the models to different calibration data and sample sizes. However, by directly accounting for the irrigation pumping uncertainties during the calibration procedures, the proposed IUWLS is able to minimize the bias effectively without adding significant computational burden to the calibration processes. PMID:25040235

  3. ShowFlow: A practical interface for groundwater modeling

    SciTech Connect

    Tauxe, J.D.

    1990-12-01

    ShowFlow was created to provide a user-friendly, intuitive environment for researchers and students who use computer modeling software. What traditionally has been a workplace available only to those familiar with command-line based computer systems is now within reach of almost anyone interested in the subject of modeling. In the case of this edition of ShowFlow, the user can easily experiment with simulations using the steady state gaussian plume groundwater pollutant transport model SSGPLUME, though ShowFlow can be rewritten to provide a similar interface for any computer model. Included in this thesis is all the source code for both the ShowFlow application for Microsoft{reg sign} Windows{trademark} and the SSGPLUME model, a User's Guide, and a Developer's Guide for converting ShowFlow to run other model programs. 18 refs., 13 figs.

  4. Developing an Analysis Program to Estimate and Prediction Groundwater Droughts in Korea from Groundwater Time-Series Data.

    NASA Astrophysics Data System (ADS)

    Cho, S.; Woo, N. C.; Lee, J. M.

    2015-12-01

    This study is aimed at developing process to analyze and predict groundwater drought potentials for Winter and Spring droughts in Korea using a long-term groundwater monitoring data. So far, most drought researches have been focused on precipitation and stream-flow data, although these data are considered to be non-linear. Subsequently, the prediction of drought events has been very difficult in practice. In this study, we targets to analyze the groundwater system as an intermediate stage between precipitation and stream-flow, but still has semi-linear characteristics. By the analysis of past trends of groundwater time-series compared with drought events, we will identify characteristics of fluctuation between groundwater-level and precipitation of the year before the droughts. Then, the characteristics will be tested with recent drought events in Korea. For this analysis, The updated ATGT (Analysis Tool for Groundwater Time-series data program version 1.0 based on JAVA), that was developed for analyzing and presenting groundwater time-series data, basically to identify abnormal changes in groundwater fluctuations, will be presented with additional functions including cross-correlation between groundwater and drought based on the PYTHON language.

  5. Basin wide Nitrate-Nitrogen pollution of groundwater, Miyakonojo, Japan, with the relation of the regional Groundwater flow system

    NASA Astrophysics Data System (ADS)

    Mikami, K.; Shimada, J.; Zikuzono, Y.

    2006-12-01

    Miyakonojo basin is well-known agriculture area in Southern Kyushu, Japan and highly depends on groundwater resources for their everyday use. Local unconfined groundwater aquifer is widely polluted by Nitrate-Nitrogen originated from agriculture. It will become serious problem if this unconfined Nitrate pollution enlarges into the confined aquifer system which is used for local city water source. However, the detailed groundwater flow system between unconfined and confined aquifer system has not been cleared yet. The detailed three dimensional groundwater flow system study has been done by using existing wells in a basin to understand the three dimensional distribution pattern of Nitrate-Nitrogen in the aquifer. The field sampling for unconfined, intermediate and confined groundwater was done in July, 2005 and February, 2006 for about 200 wells to analyze inorganic water chemistry, hydrogen / oxygen stable isotopes and tritium. For the unconfined groundwater, there exists clear difference for the groundwater flow pattern between the eastern and western basin, which is mostly affected by the surface topography. The unconfined groundwater flowed into the confined aquifer at the eastern part of a basin, while in the western part of a basin the unconfined groundwater on a plateau flowed into the confined aquifer somehow, but most part of the unconfined groundwater has been discharge out to small river valleys between plateaus. While for the confined groundwater, the topographic effect has been disappeared and basin scale groundwater flow from the basin margin toward the basin center is dominated. In the unconfined aquifer, basin wide distribution of Nitrate-Nitrogen content has been recognized and it is relatively higher in the western basin where the cattle farming are dominated. While in the confined aquifer, there are some high Nitrate-Nitrogen spots but do not have regional trend. It is considered that some part of the basin has not distributed the welded tuff

  6. Geochemical and Isotopic Interpretations of Groundwater Flow in the Oasis Valley Flow System, Southern Nevada

    SciTech Connect

    J.M. Thomas; F.C. Benedict, Jr.; T.P. Rose; R.L. Hershey; J.B. Paces; Z.E. Peterman; I.M. Farnham; K.H. Johannesson; A.K. Singh; K.J. Stetzenbach; G.B. Hudson; J.M. Kenneally; G.F. Eaton; D.K. Smith

    2003-01-08

    This report summarizes the findings of a geochemical investigation of the Pahute Mesa-Oasis Valley groundwater flow system in southwestern Nevada. It is intended to provide geochemical data and interpretations in support of flow and contaminant transport modeling for the Western and Central Pahute Mesa Corrective Action Units.

  7. Groundwater flow system and Nitrogen cycle in volcanic aquifer of pyroclastic flow uplands, Japan

    NASA Astrophysics Data System (ADS)

    Mikami, K.; Shimada, J.; Tashiro, S.; Niimi, H.

    2007-12-01

    Study area is well-known agriculture area in Southern Kyushu, Japan and highly depends on groundwater resources for their everyday use. Local unconfined groundwater aquifer is widely polluted by Nitrate-Nitrogen originated from agriculture and cattle farming. It will become serious problem if this unconfined Nitrate pollution enlarges into the confined aquifer system which is used for local city water source. The detailed three dimensional groundwater flow system study has been done by using existing wells in the basin to understand the three dimensional distribution pattern of Nitrate-Nitrogen in the aquifer. However, the detailed groundwater age analysis by using Tritium for unconfined and confined groundwater has not been succeeded because of present low atmosphere tritium concentration. Thus we applied to challenge the CFCs dating method. Although the CFCs method has been widely used for dating the young groundwater instead of tritium in many countries, in Japan CFCs has been used only by Oceanographic study and has not been used in the field of Hydrology. The history and fate of Nitrate contamination have been shown in multidisciplinary local transect studies in areas with agricultural sources (Bohlke and Denver 1995). However, identification of Nitrogen sources can be difficult in larger regional studies because of co-occurrence of multiple anthropogenic Nitrogen sources and uncertainty in Nitrogen transformation pathways. Thus, the characterization of N geochemistry remains challenging, particularly in aquifer-scale assessments (Stephen 2006). In this study, the evidence of the shallow groundwater flowing towards deep aquifer was verified by the groundwater dating and the detailed Nitrogen reduction process was confirmed along the groundwater flow.

  8. Processing, Analysis, and General Evaluation of Well-Driller Logs for Estimating Hydrogeologic Parameters of the Glacial Sediments in a Ground-Water Flow Model of the Lake Michigan Basin

    USGS Publications Warehouse

    Arihood, Leslie D.

    2009-01-01

    In 2005, the U.S. Geological Survey began a pilot study for the National Assessment of Water Availability and Use Program to assess the availability of water and water use in the Great Lakes Basin. Part of the study involves constructing a ground-water flow model for the Lake Michigan part of the Basin. Most ground-water flow occurs in the glacial sediments above the bedrock formations; therefore, adequate representation by the model of the horizontal and vertical hydraulic conductivity of the glacial sediments is important to the accuracy of model simulations. This work processed and analyzed well records to provide the hydrogeologic parameters of horizontal and vertical hydraulic conductivity and ground-water levels for the model layers used to simulated ground-water flow in the glacial sediments. The methods used to convert (1) lithology descriptions into assumed values of horizontal and vertical hydraulic conductivity for entire model layers, (2) aquifer-test data into point values of horizontal hydraulic conductivity, and (3) static water levels into water-level calibration data are presented. A large data set of about 458,000 well driller well logs for monitoring, observation, and water wells was available from three statewide electronic data bases to characterize hydrogeologic parameters. More than 1.8 million records of lithology from the well logs were used to create a lithologic-based representation of horizontal and vertical hydraulic conductivity of the glacial sediments. Specific-capacity data from about 292,000 well logs were converted into horizontal hydraulic conductivity values to determine specific values of horizontal hydraulic conductivity and its aerial variation. About 396,000 well logs contained data on ground-water levels that were assembled into a water-level calibration data set. A lithology-based distribution of hydraulic conductivity was created by use of a computer program to convert well-log lithology descriptions into aquifer or

  9. MODOPTIM: A general optimization program for ground-water flow model calibration and ground-water management with MODFLOW

    USGS Publications Warehouse

    Halford, Keith J.

    2006-01-01

    MODOPTIM is a non-linear ground-water model calibration and management tool that simulates flow with MODFLOW-96 as a subroutine. A weighted sum-of-squares objective function defines optimal solutions for calibration and management problems. Water levels, discharges, water quality, subsidence, and pumping-lift costs are the five direct observation types that can be compared in MODOPTIM. Differences between direct observations of the same type can be compared to fit temporal changes and spatial gradients. Water levels in pumping wells, wellbore storage in the observation wells, and rotational translation of observation wells also can be compared. Negative and positive residuals can be weighted unequally so inequality constraints such as maximum chloride concentrations or minimum water levels can be incorporated in the objective function. Optimization parameters are defined with zones and parameter-weight matrices. Parameter change is estimated iteratively with a quasi-Newton algorithm and is constrained to a user-defined maximum parameter change per iteration. Parameters that are less sensitive than a user-defined threshold are not estimated. MODOPTIM facilitates testing more conceptual models by expediting calibration of each conceptual model. Examples of applying MODOPTIM to aquifer-test analysis, ground-water management, and parameter estimation problems are presented.

  10. Challenges in Estimating Groundwater Recharge in Semi-arid and Semi-confined Alluvial Systems

    NASA Astrophysics Data System (ADS)

    Larsen, J.; Finch, W.; McIntyre, N.

    2015-12-01

    Uncertainty surrounding rates of groundwater recharge limits overall confidence in groundwater allocations and can lead to over-conservative assumptions in groundwater impact assessments. This problem is even more acute where more complex unsaturated flow paths are considered, such as within semi-confined alluvial systems. Researchers at The University of Queensland have developed an experimental study within the Condamine Alluvium, a significant aquifer in semi-arid eastern Australia, is used to determine groundwater recharge mechanisms for three distinct soil types (two vertosols and one chromosol) on both irrigated and non-irrigated areas. This variety of soil types, including shrink-swell clays, overly a higher permeability sand and gravel unsaturated zone and aquifer. The analysis uses 15-minute soil moisture data from Sentek EnviroSCAN Probe devices installed at 16 sites, with eight sensors in each site at depths from 100 to 4000mm. The vertosols exhibited signs of dynamic preferential flow paths due to the shrink swell properties of the soil. Precipitation rate and initial soil moisture content affect the infiltration response times for the three soils, with the chromosols requiring multiple precipitation events before experiencing any significant soil moisture storage changes in the lower depths (2000 - 4000mm). Storage changes below the root zone to a depth of 4m indicate large rates of potential recharge, up to 1300mm for the two years of data obtained. However, minimal rise has been observed in the water table (~12 m depth), potentially due to the highly transmissive sand and gravel aquifer. The analysis has shown that only very high temporal resolution monitoring of soil storage changes can effectively capture the dynamic preferential flow water flux. Lower temporal resolution monitoring, at the daily scale or greater, will bias the storage change estimates towards the matrix flow component and risk significant underestimation of the total unsaturated

  11. On the mechanism of earthquake induced groundwater flow

    NASA Astrophysics Data System (ADS)

    Dudley Ward, Nicholas F.

    2015-11-01

    The Canterbury/Christchurch earthquakes and aftershocks of 2010-2011 generated groundwater level responses throughout New Zealand. However, the greater part of damage was sustained by the city of Christchurch which is built on a layered sequence of artesian aquifers. In a previous paper we analysed the earthquake induced groundwater responses. We quantified groundwater responses by fitting a simple statistical model which differentiated between immediate earthquake induced response (spike) and post-seismic change (offset). The most notable feature of this analysis was the consistency of responses between the earthquakes: deeper wells correlate with negative offset and shallower wells correlate with positive offset. In that paper we argued that this is consistent with the upwards vertical movement of water. In this paper we focus on the physical mechanisms, and consider a model that further explains and supports this hypothesis. We postulate a groundwater flow model in which storativity and aquitard permeability are modelled as time-varying shocks. We analyse the solutions for a range of non-dimensional parameters and obtain type curves that exhibit the same behaviour as the observed responses. Finally we consider data from the 2010 Mw 7.1 Darfield (Canterbury) earthquake.

  12. An approach to improve direct runoff estimates and reduce uncertainty in the calculated groundwater component in water balances of large lakes

    NASA Astrophysics Data System (ADS)

    Wiebe, Andrew J.; Conant, Brewster; Rudolph, David L.; Korkka-Niemi, Kirsti

    2015-12-01

    Groundwater is important in the overall water budget of a lake because it affects the quantity and quality of surface water and the ecological health of the lake. The water balance equation is frequently used to estimate the net groundwater flow for small lakes but is seldom used to determine net groundwater flow components for large lakes because: (1) errors accumulate in the calculated groundwater term, and (2) there is an inability to accurately quantify the direct runoff component. In this water balance study of Lake Pyhäjärvi (155 km2) in Finland, it was hypothesized a hydrograph separation model could be used to estimate direct runoff to the lake and, when combined with a rigorous uncertainty analyses, would provide reliable net groundwater flow estimates. The PART hydrograph separation model was used to estimate annual per unit area direct runoff values for the watershed of the inflowing Yläneenjoki River (a subwatershed of the lake) which were then applied to other physically similar subwatersheds of the lake to estimate total direct runoff to the lake. The hydrograph separation method provided superior results and had lower uncertainty than the common approach of using a runoff coefficient based method. The average net groundwater flow into the lake was calculated to be +43 mm per year (+3.0% of average total inflow) for the 38 water years 1971-2008. It varied from -197 mm to 284 mm over that time, and had a magnitude greater than the uncertainty for 17 of the 38 years. The average indirect groundwater contribution to the lake (i.e., the groundwater part of the inflowing rivers) was 454 mm per year (+32% of average total inflow) and demonstrates the overall importance of groundwater. The techniques in this study are applicable to other large lakes and may allow small net groundwater flows to be reliably quantified in settings that might otherwise be unquantifiable or completely lost in large uncertainties.

  13. Striking effect of time variation in the estimation of groundwater age in the Wairarapa valley

    NASA Astrophysics Data System (ADS)

    Petrus, Karine; Toews, Michael; Daughney, Christopher; Cornaton, Fabien

    2014-05-01

    The Wairarapa Valley exhibits complex interactions between its rivers and shallow aquifers. With agriculture being an essential part of the region the risk of contamination and depletion of groundwater exists. In order to assist with water resource management in the region, we can do predictions with the help of numerical models. Among these predictions, the evaluation of groundwater age is critical for decision making. This project builds on work done by Greater Wellington Regional Council and will focus on the Wairarapa Valley. The aim of this study is to evaluate the age of the groundwater in the Wairarapa region. Investigations have already been done thanks to hydrochemistry. However radiometric age can be misleading in the sense that it does not consider the mixing process in the motion of groundwater particules. Therefore another approach can be considered .This latter is physic based by considering the age as a property that we transport through two main processes: advection at a macroscopic scale and diffusion at a microscopic scale. The determination of the distribution age by this approach has already been done for the Lake Rotorua but in the steady state case (cf Daughney). The unique contribution of the present study is to estimate the changes in groundwater age distribution through time within the region. Indeed transient simulations are needed to explicitly account for seasonally variable rainfall and pumping wells. This affects the simulated flow solution and then the simulated age solution. In order to solve numerically the transport of age distribution we have chosen to use the Time Marching Laplace Transform Galerkin technique which has been developed in a research code by Fabien Cornaton. The obtained results depict that temporal variations in groundwater age are present and have important implication for resource management

  14. Regional stochastic estimation of the groundwater catchment for distributed hydrological modelling

    NASA Astrophysics Data System (ADS)

    Wöhling, Th.; Samaniego, L.; Selle, B.; Kumar, R.; Zink, M.

    2012-04-01

    Rainfall-runoff modeling typically assumes that the groundwater catchment boundary coincide with the topographic one. While this is often a reasonable assumption for large and and mesoscale catchments (> 103 km2), this assumption may lead to large errors of streamflow in small scale catchments (≤ 102 km2), in particular in certain geological settings. The Ammer catchment (135 km2) in the upper Neckar river basin (Germany) is a prime example where groundwater and topographic catchment boundaries are significantly distinct from each other. The catchment is characterized by a complex sequence of fractured, karstic Triassic rock formations. These strata gently dip into ESE direction governing groundwater flow. Analysis of tracer experiments conducted in the 1970s indicates that the boundary overlap could be less than 80 percent. Further, a modelling study of the upper Neckar river basin using the distributed hydrological model mHM showed Nash-Sutcliff efficiencies (NSE) < 0.4 for simulated runoff in the Ammer sub-basin whereas higher efficiencies (NSE ~ 0.7) were obtained for most of the other 21 sub basins in the region. In this study we present a methodology to simultaneously estimate the regional groundwater catchment boundaries of the Ammer and its surrounding basins. In a first step we derive the best possible fit between mHM simulated and observed runoff for the individual sub-basins in the Ammer region and determine the trade-off between the fits of the individual basins using the muliobjective optimization method AMALGAM. We further present a strategy to estimate the regional groundwater catchment boundaries with the aim to improve runoff predictions in the Ammer catchment while not deteriorating runoff predictions in the surrounding basins. Our strategy involves a modification of the mHM model to account for ground water import/export from neighboring catchments while maintaining full mass balance of the surrounding basins. Groundwater catchment boundaries

  15. Tide-induced fingering flow during submarine groundwater discharge

    NASA Astrophysics Data System (ADS)

    Greskowiak, Janek

    2013-04-01

    Submarine groundwater discharge (SGD) is a relevant component of the hydrological cycle (Moore, 2010). The discharge of fresh groundwater that originated from precipitation on the land typically occurs at the near shore scale (~ 10m-100m) and the embayment scale (~ 100m - 10km) (Bratton, 2010). In the recent years a number of studies revealed that tidal forcing has an important effect on the fresh SGD pattern in the beach zone, i.e., it leads to the formation of an upper saline recirculation cell and a lower "freshwater discharge tube" (Boufadel, 2000, Robinson et al., 2007; Kuan et al., 2012). Thereby the discharge of the fresh groundwater occurs near the low-tide mark. The shape and extent of the upper saline recirculation cell is mainly defined by the tidal amplitude, beach slope, fresh groundwater discharge rate and hydraulic conductivity (Robinson et al., 2007). In spite of fact that in this case sea water overlies less denser freshwater, all previous modeling studies suggested that the saline recirculation cell and the freshwater tube are rather stable. However, new numerical investigations indicate that there maybe realistic cases where the upper saline recirculation cell becomes unstable as a result of the density contrast to the underlying freshwater tube. In these cases salt water fingers develop and move downward, thereby penetrating the freshwater tube. To the author's knowledge, the present study is the first that illustrate the possibility of density induced fingering flow during near shore SGD. A total of 240 high resolution simulations with the density dependent groundwater modelling software SEAWAT-2000 (Langevin et al., 2007) has been carried out to identify the conditions under which salt water fingering starts to occur. The simulations are based on the field-scale model setup employed in Robinson et al. (2007). The simulation results indicate that a very flat beach slope of less than 1:35, a hydraulic conductivity of 10 m/d and already a tidal

  16. Deterministic modelling of the cumulative impacts of underground structures on urban groundwater flow and the definition of a potential state of urban groundwater flow: example of Lyon, France

    NASA Astrophysics Data System (ADS)

    Attard, Guillaume; Rossier, Yvan; Winiarski, Thierry; Cuvillier, Loann; Eisenlohr, Laurent

    2016-02-01

    Underground structures have been shown to have a great influence on subsoil resources in urban aquifers. A methodology to assess the actual and the potential state of the groundwater flow in an urban area is proposed. The study develops a three-dimensional modeling approach to understand the cumulative impacts of underground infrastructures on urban groundwater flow, using a case in the city of Lyon (France). All known underground structures were integrated in the numerical model. Several simulations were run: the actual state of groundwater flow, the potential state of groundwater flow (without underground structures), an intermediate state (without impervious structures), and a transient simulation of the actual state of groundwater flow. The results show that underground structures fragment groundwater flow systems leading to a modification of the aquifer regime. For the case studied, the flow systems are shown to be stable over time with a transient simulation. Structures with drainage systems are shown to have a major impact on flow systems. The barrier effect of impervious structures was negligible because of the small hydraulic gradient of the area. The study demonstrates that the definition of a potential urban groundwater flow and the depiction of urban flow systems, which involves understanding the impact of underground structures, are important issues with respect to urban underground planning.

  17. Deterministic modelling of the cumulative impacts of underground structures on urban groundwater flow and the definition of a potential state of urban groundwater flow: example of Lyon, France

    NASA Astrophysics Data System (ADS)

    Attard, Guillaume; Rossier, Yvan; Winiarski, Thierry; Cuvillier, Loann; Eisenlohr, Laurent

    2016-08-01

    Underground structures have been shown to have a great influence on subsoil resources in urban aquifers. A methodology to assess the actual and the potential state of the groundwater flow in an urban area is proposed. The study develops a three-dimensional modeling approach to understand the cumulative impacts of underground infrastructures on urban groundwater flow, using a case in the city of Lyon (France). All known underground structures were integrated in the numerical model. Several simulations were run: the actual state of groundwater flow, the potential state of groundwater flow (without underground structures), an intermediate state (without impervious structures), and a transient simulation of the actual state of groundwater flow. The results show that underground structures fragment groundwater flow systems leading to a modification of the aquifer regime. For the case studied, the flow systems are shown to be stable over time with a transient simulation. Structures with drainage systems are shown to have a major impact on flow systems. The barrier effect of impervious structures was negligible because of the small hydraulic gradient of the area. The study demonstrates that the definition of a potential urban groundwater flow and the depiction of urban flow systems, which involves understanding the impact of underground structures, are important issues with respect to urban underground planning.

  18. Groundwater flow and radionuclide decay-chain transport modelling around a proposed uranium tailings pond in India

    NASA Astrophysics Data System (ADS)

    Elango, L.; Brindha, K.; Kalpana, L.; Sunny, Faby; Nair, R. N.; Murugan, R.

    2012-06-01

    Extensive hydrogeological investigations followed by three-dimensional groundwater flow and contaminant transport modelling were carried out around a proposed uranium tailings pond at Seripalli in Andhra Pradesh, India, to estimate its radiological impact. The hydrogeological parameters and measured groundwater level were used to model the groundwater flow and contaminant transport from the uranium tailings pond using a finite-element-based model. The simulated groundwater level compares reasonably with the observed groundwater level. Subsequently, the transport of long-lived radionuclides such as 238U, 234U, 230Th and 226Ra from the proposed tailings pond was modelled. The ingrowths of progenies were also considered in the modelling. It was observed that these radionuclides move very little from the tailings pond, even at the end of 10,000 y, due to their high distribution coefficients and low groundwater velocities. These concentrations were translated into committed effective dose rates at different distances in the vicinity of the uranium tailings pond. The results indicated that the highest effective dose rate to members of the public along the groundwater flow pathway is 2.5 times lower than the drinking water guideline of 0.1 mSv/y, even after a long time period of 10,000 y.

  19. From groundwater baselines to numerical groundwater flow modelling for the Milan metropolitan area

    NASA Astrophysics Data System (ADS)

    Crosta, Giovanni B.; Frattini, Paolo; Peretti, Lidia; Villa, Federica; Gorla, Maurizio

    2015-04-01

    allow for the groundwater flow and transport modeling at the large scale and could be successively linked to some more site-specific transport multi-reactive models focused on the modeling of some specific contaminants.

  20. Mechanism and rate of denitrification in an agricultural watershed: Electron and mass balance along groundwater flow paths

    USGS Publications Warehouse

    Tesoriero, A.J.; Liebscher, H.; Cox, S.E.

    2000-01-01

    The rate and mechanism of nitrate removal along and between groundwater flow paths were investigated using a series of well nests screened in an unconfined sand and gravel aquifer. Intensive agricultural activity in this area has resulted in nitrate concentrations in groundwater often exceeding drinking water standards. Both the extent and rate of denitrification varied depending on the groundwater flow path. While little or no denitrification occurred in much of the upland portions of the aquifer, a gradual redox gradient is observed as aerobic upland groundwater moves deeper in the aquifer. In contrast, a sharp shallow redox gradient is observed adjacent to a third-order stream as aerobic groundwater enters reduced sediments. An essentially complete loss of nitrate concurrent with increases in excess N2 provide evidence that denitrification occurs as groundwater enters this zone. Electron and mass balance calculations suggest that iron sulfide (e.g., pyrite) oxidation is the primary source of electrons for denitrification. Denitrification rate estimates were based on mass balance calculations using nitrate and excess N2 coupled with groundwater travel times. Travel times were determined using a groundwater flow model and were constrained by chlorofluorocarbon-based age dates. Denitrification rates were found to vary considerably between the two areas where denitrification occurs. Denitrification rates in the deep, upland portions of the aquifer were found to range from < 0.01 to 0.14 mM of N per year; rates at the redoxcline along the shallow flow path range from 1.0 to 2.7 mM of N per year. Potential denitrification rates in groundwater adjacent to the stream may be much faster, with rates up to 140 mM per year based on an in situ experiment conducted in this zone.The rate and mechanism of nitrate removal along and between groundwater flow paths were investigated using a series of well nests screened in an unconfined sand and gravel aquifer. Intensive

  1. Interpretation of groundwater flow patterns through a reconstruction of the tritium precipitation record in the Cochabamba Valley, Bolivia

    NASA Astrophysics Data System (ADS)

    Stimson, Jesse; Rudolph, David; Frape, Shaun; Drimmie, Robert

    1996-05-01

    Regional behaviour of the groundwater flow system in the Cochabamba Valley, Bolivia, is evaluated through the interpretation of tritium ( 3H) distributions in groundwater samples from wells and springs. In order to interpret groundwater 3H concentrations in Cochabamba Valley, where no historical record of 3H concentrations in rainfall exists, a reconstructed 3H precipitation record is developed. The record of 3H concentrations in precipitation is fairly extensive in the Amazon Basin and this record was extrapolated to the neighbouring Cochabamba Valley. Tritium concentrations in rainfall have been observed to increase under natural conditions with increasing latitude and with increasing distance from the ocean. By considering these trends, a linear relationship for increasing 3H concentration in precipitation is developed, based on data from the Amazon Basin, that realistically predicts regional 3H distributions from the northeast Brazilian coast to Cuzco, Peru. This 3H precipitation record is then extrapolated to the Cochabamba Valley and, after correction for radiogenic decay, is used to interpret trends in groundwater 3H concentrations within the valley. The groundwater flow system in one of the principal alluvial fans, which serves as an important groundwater resource for the city, is studied in detail. Tritium concentrations drop from approximately 8-10 tritium units (TU) in the recharge area to concentrations below the detection limit of 0.8 TU further out in the valley. Groundwater velocities of approximately 0.3 to 0.9 m d -1 are estimated from distributions of groundwater 3H concentrations along the alluvial fan with the use of the reconstructed precipitation 3H record. Regional characteristics of the groundwater flow system are discussed with respect to future development and protection of the groundwater resources.

  2. Experimental and numerical investigations of soil water balance at the hinterland of the Badain Jaran Desert for groundwater recharge estimation

    NASA Astrophysics Data System (ADS)

    Hou, Lizhu; Wang, Xu-Sheng; Hu, Bill X.; Shang, Jie; Wan, Li

    2016-09-01

    Quantification of groundwater recharge from precipitation in the huge sand dunes is an issue in accounting for regional water balance in the Badain Jaran Desert (BJD) where about 100 lakes exist between dunes. In this study, field observations were conducted on a sand dune near a large saline lake in the BJD to investigate soil water movement through a thick vadose zone for groundwater estimation. The hydraulic properties of the soils at the site were determined using in situ experiments and laboratory measurements. A HYDRUS-1D model was built up for simulating the coupling processes of vertical water-vapor movement and heat transport in the desert soil. The model was well calibrated and validated using the site measurements of the soil water and temperature at various depths. Then, the model was applied to simulate the vertical flow across a 3-m-depth soil during a 53-year period under variable climate conditions. The simulated flow rate at the depth is an approximate estimation of groundwater recharge from the precipitation in the desert. It was found that the annual groundwater recharge would be 11-30 mm during 1983-2012, while the annual precipitation varied from 68 to 172 mm in the same period. The recharge rates are significantly higher than those estimated from the previous studies using chemical information. The modeling results highlight the role of the local precipitation as an essential source of groundwater in the BJD.

  3. Evaluation of groundwater flow patterns around a dual-screened groundwater circulation well.

    PubMed

    Johnson, Richard L; Simon, Michelle A

    2007-08-15

    Dual-screened groundwater circulation wells (GCWs) can be used to remove contaminant mass and to mix reagents in situ. GCWs are so named because they force water in a circular pattern between injection and extraction screens. The radial extent, flux and direction of the effective flow of this circulation cell are difficult to measure or predict. The objective of this study is to develop a robust protocol for assessing GCW performance. To accomplish this, groundwater flow patterns surrounding a GCW are assessed using a suite of tools and data, including: hydraulic head, in situ flow velocity, measured hydraulic conductivity data from core samples, chemical tracer tests, contaminant distribution data, and numerical flow and transport models. The hydraulic head data show patterns that are consistent with pumping on a dual-screened well, however, many of the observed changes are smaller than expected. In situ thermal perturbation flow sensors successfully measured horizontal flow, but vertical flow could not be determined with sufficient accuracy to be useful in mapping flow patterns. Two types of chemical tracer tests were utilized at the site and showed that much of the flow occurs within a few meters of the GCW. Flow patterns were also assessed based on changes in contaminant (trichloroethylene, TCE) concentrations over time. The TCE data clearly showed treated water moving away from the GCW at shallow and intermediate depths, but the circulation of that water back to the well, except very close to the well, was less clear. Detailed vertical and horizontal hydraulic conductivities were measured on 0.3 m-long sections from a continuous core from the GCW installation borehole. The measured vertical and horizontal hydraulic conductivity data were used to construct numerical flow and transport models, the results of which were compared to the head, velocity and concentration data. Taken together, the field data and modeling present a fairly consistent picture of flow

  4. Numerical modeling of groundwater flow in the coastal aquifer system of Taranto (southern Italy)

    NASA Astrophysics Data System (ADS)

    De Filippis, Giovanna; Giudici, Mauro; Negri, Sergio; Margiotta, Stefano; Cattaneo, Laura; Vassena, Chiara

    2014-05-01

    The Mediterranean region is characterized by a strong development of coastal areas with a high concentration of water-demanding human activities, resulting in weakly controlled withdrawals of groundwater which accentuate the saltwater intrusion phenomenon. The worsening of groundwater quality is a huge problem especially for those regions, like Salento (southern Italy), where a karst aquifer system represents the most important water resource because of the deficiency of a well developed superficial water supply. In this frame, the first 2D numerical model describing the groundwater flow in the karst aquifer of Salento peninsula was developed by Giudici et al. [1] at the regional scale and then improved by De Filippis et al. [2]. In particular, the estimate of the saturated thickness of the deep aquifer highlighted that the Taranto area is particularly sensitive to the phenomenon of seawater intrusion, both for the specific hydrostratigraphic configuration and for the presence of highly water-demanding industrial activities. These remarks motivate a research project which is part of the research program RITMARE (The Italian Research for the Sea), within which a subprogram is specifically dedicated to the problem of the protection and preservation of groundwater quality in Italian coastal aquifers and in particular, among the others, in the Taranto area. In this context, the CINFAI operative unit aims at providing a contribution to the characterization of groundwater in the study area. The specific objectives are: a. the reconstruction of the groundwater dynamic (i.e., the preliminary identification of a conceptual model for the aquifer system and the subsequent modeling of groundwater flow in a multilayered system which is very complex from the hydrostratigraphical point of view); b. the characterization of groundwater outflows through submarine and subaerial springs and the water exchanges with the shallow coastal water bodies (e.g. Mar Piccolo) and the off

  5. Using noble gas tracers to constrain a groundwater flow model with recharge elevations: A novel approach for mountainous terrain

    USGS Publications Warehouse

    Doyle, Jessica M.; Gleeson, Tom; Manning, Andrew H.; Mayer, K. Ulrich

    2015-01-01

    Environmental tracers provide information on groundwater age, recharge conditions, and flow processes which can be helpful for evaluating groundwater sustainability and vulnerability. Dissolved noble gas data have proven particularly useful in mountainous terrain because they can be used to determine recharge elevation. However, tracer-derived recharge elevations have not been utilized as calibration targets for numerical groundwater flow models. Herein, we constrain and calibrate a regional groundwater flow model with noble-gas-derived recharge elevations for the first time. Tritium and noble gas tracer results improved the site conceptual model by identifying a previously uncertain contribution of mountain block recharge from the Coast Mountains to an alluvial coastal aquifer in humid southwestern British Columbia. The revised conceptual model was integrated into a three-dimensional numerical groundwater flow model and calibrated to hydraulic head data in addition to recharge elevations estimated from noble gas recharge temperatures. Recharge elevations proved to be imperative for constraining hydraulic conductivity, recharge location, and bedrock geometry, and thus minimizing model nonuniqueness. Results indicate that 45% of recharge to the aquifer is mountain block recharge. A similar match between measured and modeled heads was achieved in a second numerical model that excludes the mountain block (no mountain block recharge), demonstrating that hydraulic head data alone are incapable of quantifying mountain block recharge. This result has significant implications for understanding and managing source water protection in recharge areas, potential effects of climate change, the overall water budget, and ultimately ensuring groundwater sustainability.

  6. Estimation of yield capacity of fractured rock aquifer for multi-well groundwater heat pump system

    NASA Astrophysics Data System (ADS)

    Bak, Hyeongmin; Yeo, In Wook

    2015-04-01

    Geothermal heat pump system is classified as closed loop and open loop. Closed loop uses a refrigerant as a heat source. For the reason, when using it for a long time, there is a possibility that the refrigerant pipe is corroded. Accordingly, soil and groundwater can be contaminated. Whereas the open loop system uses a eco-friendly groundwater as a heat source. Thermal circulation of standing column well (SCW) occurs in one well. In contrast, thermal circulation of multi-well groundwater heat pump system (MGHP) occurs through fractured rock aquifer between extraction and injection wells. Therefore, temperature efficiency of MGHP appears to be better than that of SCW. However, the MGHP has problems such as the overflowing in the injection well and the clogging, which restricts the wide use of MGHP. This study aims at how to to array the extraction and injection wells for stable circulating of groundwater and at evaluating the sustainable yield capacity of groundwater circulation between the two wells. The study site is located in Chuncheon, Republic of Korea. Pumping tests were conducted to estimate transmissivity of the two wells (W3, W4). In addition, the step-circulation tests were conducted to estimate the sustainable yield capacity. Transmissivity of W3 and W4 was estimated to be 5.81 x 10^-5 m^2/s and 2.57 x 10^-5 m^2/s, respectively. Preliminary groundwater circulation tests were conducted to figure out the array of the extraction and injection wells. Circulation tests were performed for two cases: first, extraction well was set at the well with higher transmissivity and injection well set at the well with lower transmissivity, and the opposite array was set for the second case. In the first case, when flow rate was set at 70.47 m^3/day, the water level of W3 fell 0.61m and that of W4 rose 1.89m. In the second case, when flow rate was set at 67.70 m^3/day, the water level of W4 fell 2.17m and that of W3 rose 0.5m. Preliminary groundwater

  7. Validation Analysis of the Shoal Groundwater Flow and Transport Model

    SciTech Connect

    A. Hassan; J. Chapman

    2008-11-01

    Environmental restoration at the Shoal underground nuclear test is following a process prescribed by a Federal Facility Agreement and Consent Order (FFACO) between the U.S. Department of Energy, the U.S. Department of Defense, and the State of Nevada. Characterization of the site included two stages of well drilling and testing in 1996 and 1999, and development and revision of numerical models of groundwater flow and radionuclide transport. Agreement on a contaminant boundary for the site and a corrective action plan was reached in 2006. Later that same year, three wells were installed for the purposes of model validation and site monitoring. The FFACO prescribes a five-year proof-of-concept period for demonstrating that the site groundwater model is capable of producing meaningful results with an acceptable level of uncertainty. The corrective action plan specifies a rigorous seven step validation process. The accepted groundwater model is evaluated using that process in light of the newly acquired data. The conceptual model of ground water flow for the Project Shoal Area considers groundwater flow through the fractured granite aquifer comprising the Sand Springs Range. Water enters the system by the infiltration of precipitation directly on the surface of the mountain range. Groundwater leaves the granite aquifer by flowing into alluvial deposits in the adjacent basins of Fourmile Flat and Fairview Valley. A groundwater divide is interpreted as coinciding with the western portion of the Sand Springs Range, west of the underground nuclear test, preventing flow from the test into Fourmile Flat. A very low conductivity shear zone east of the nuclear test roughly parallels the divide. The presence of these lateral boundaries, coupled with a regional discharge area to the northeast, is interpreted in the model as causing groundwater from the site to flow in a northeastward direction into Fairview Valley. Steady-state flow conditions are assumed given the absence of

  8. Ground-water flow in Melton Valley, Oak Ridge reservation, Roane County, Tennessee; preliminary model analysis

    USGS Publications Warehouse

    Tucci, Patrick

    1986-01-01

    Shallow land burial of low-level radioactive waste has been practiced since 1951 in Melton Valley. Groundwater flow modeling was used to better understand the geohydrology of the valley, and to provide a foundation for future contaminant transport modeling. The three-dimensional, finite difference model simulates the aquifer as a two layer system that represents the regolith and bedrock. Transmissivities, which were adjusted during model calibration, range from 8 to 16 sq ft/day for the regolith, and from 0.2 to 1.5 sq ft/day for bedrock. An anisotropy ratio of 1:3 for strike-normal to strike-parallel transmissivity values, in conjunction with recharge rate = 6% of precipitation that is uniformly distributed over the model area, produces the best match between simulated and observed water levels. Simulated water levels generally compare well to observed or estimated 1978 groundwater conditions. Simulated water levels for the regolith for 39 of 69 comparison points are within +/- 10 ft of average 1978 levels. Simulated vertical flow components are in the observed direction for 9 of 11 comparison points. Preliminary simulations indicate that nearly all groundwater flow is within the regolith and discharges to either the Clinch River or the White Oak Creek-Melton Branch drainage systems. Less than 3% of the flow is between the regolith and bedrock, and < 1% of total groundwater flow discharges to the Clinch River through bedrock. Additional data needed to refine and further calibrate the model, include: (1) quantity and areal distribution of recharge; (2) water levels in the regolith near the model boundaries and beyond the Clinch River; (3) water levels and aquifer characteristics for bedrock; and (4) additional surface water data. (Author 's abstract)

  9. The effect of sediment thermal conductivity on vertical groundwater flux estimates

    NASA Astrophysics Data System (ADS)

    Sebok, Eva; Müller, Sascha; Engesgaard, Peter; Duque, Carlos

    2015-04-01

    The interaction between groundwater and surface water is of great importance both from ecological and water management perspective. The exchange fluxes are often estimated based on vertical temperature profiles taken from shallow sediments assuming a homogeneous standard value of sediment thermal conductivity. Here we report on a field investigation in a stream and in a fjord, where vertical profiles of sediment thermal conductivity and temperatures were measured in order to, (i) define the vertical variability in sediment thermal conductivity, (ii) quantify the effect of heterogeneity in sediment thermal conductivity on the estimated vertical groundwater fluxes. The study was carried out at field sites located in Ringkøbing fjord and Holtum stream in Western Denmark. Both locations have soft, sandy sediments with an upper organic layer at the fjord site. First 9 and 12 vertical sediment temperature profiles up to 0.5 m depth below the sediment bed were collected in the fjord and in the stream, respectively. Later sediment cores of 0.05 m diameter were removed at the location of the temperature profiles. Sediment thermal conductivity was measured in the sediment cores at 0.1 m intervals with a Decagon KD2 Pro device. A 1D flow and heat transport model (HydroGeoSphere) was set up and vertical groundwater fluxes were estimated based on the measured vertical sediment temperature profiles by coupling the model with PEST. To determine the effect of heterogeneity in sediment thermal conductivity on estimated vertical groundwater fluxes, the model was run by assigning (i) a homogeneous thermal conductivity for all sediment layers, calculated as the average sediment thermal conductivity of the profile, (ii) measured sediment thermal conductivities to the different model layers. The field survey showed that sediment thermal conductivity over a 0.5 m profile below the sediment bed is not uniform, having the largest variability in the fjord where organic sediments were also

  10. Estimating pumping time and ground-water withdrawals using energy- consumption data

    USGS Publications Warehouse

    Hurr, R.T.; Litke, D.W.

    1989-01-01

    Evaluation of the hydrology of an aquifer requires knowledge about the volume of groundwater in storage and also about the volume of groundwater withdrawals. Totalizer flow meters may be installed at pumping plants to measure withdrawals; however, it generally is impractical to equip all pumping plants in an area with meters. A viable alternative is the use of rate-time methods. Rate-time methods may be used at individual pumping plants to decrease the data collection necessary for determining withdrawals. At sites where pumping-time measurement devices are not installed, pumping time may be determined on the basis of energy consumption and power demand. At pumping plants where energy consumption is metered, data acquired by reading of meters is used to estimate pumping time. Care needs to be taken to read these meters correctly. At pumping plants powered by electricity, the calculations need to be modified if transformers are present. At pumping plants powered by natural gas, the effects of the pressure-correction factor need to be included in the calculations. At pumping plants powered by gasoline, diesel oil, or liquid petroleum gas, the geometry of storage tanks needs to be analyzed as part of the calculations. The relation between power demand and pumping rate at a pumping plant can be described through the use of the power-consumption coefficient. Where equipment and hydrologic conditions are stable, this coefficient can be applied to total energy consumption at a site to estimate total groundwater withdrawals. Random sampling of power consumption coefficients can be used to estimate area-wide groundwater withdrawal. (USGS)