Geology topics

Groundwater & aquifers

Find USGS groundwater and aquifer publications. Search an aquifer name or place, compare source-supplied dates and open reports for hydrogeology, recharge and water-level studies.

Follow this topic with RSS

At least 1,333 records · Page 74Linked to original sources

USGS1989

Geohydrology of the regional aquifer system, western Snake River plain, southwestern Idaho

A three dimensional groundwater flow model was developed to simulate steady state and nonsteady-state hydrologic conditions of the regional aquifer system in the western Snake River Plain of Idaho. Water budget analysis showed that groundwater recharge was about 1,400,000 acre-ft in 1980; groundwater pumpage was estimated to be 300,000 acre-ft. Two mass water level measurements were made in March and August 1980 to define the water table in the regional system. The model was discretized into 25 rows, 72 columns, and 3 layers. Each cell represented 4 sq mi. The model was calibrated to 1980 hydrologic conditions. Calibrated transmissivity of layer 1 (500 ft thick) ranged from 1,500 to 21,500 sq ft/day. Calibrated specific yield of unconfined aquifers was 0.10 and calibrated storage coefficient of confined aquifers ranged from 0.0004 to 0.007. The calibrated model was verified by simulating monthly water-level fluctuations for 1980. Simulated water levels matched measured levels in the Boise River Valley, but the match in other areas was poor. (USGS)

Open-File Report↗

USGS1989

Conversion and comparison of the mathematical, three-dimensional, finite-difference, ground-water flow model to the modular, three-dimensional, finite-difference, ground-water flow model for the Tesuque aquifer system in northern New Mexico

The three-dimensional finite-difference groundwater model (using a mathematical groundwater flow code) of the Tesuque aquifer system in northern New Mexico was converted to run using the U.S. Geological Survey 's modular groundwater flow code. Results from the final versions of the predevelopment and 1947 to 2080 transient simulations of the two models are compared. A correlation coefficient of 0.9905 was obtained for the match in block-by-block head-dependent fluxes for predevelopment conditions. There are, however, significant differences in at least two specific cases. In the first case, a difference is associated with the net loss from the Pojoaque River and its tributaries to the aquifer. The net loss by the river is given as 1.134 cu ft/sec using the original groundwater model, which is 38.1% less than the net loss by the river of 1.8319 cu ft/sec computed in this study. In the second case, the large difference is computed for the transient decline in the hydraulic head of a model block near Tesuque Pueblo. The hydraulic-head decline by 2080 is, using the original model, 249 ft, which is 14.7% less than the hydraulic head of 292 ft computed by this study. In general, the differences between the two sets of results are not large enough to lead to different conclusions regarding the behavior of the system at steady state or when pumped. (USGS)

Open-File Report↗

USGS2010

Relations between groundwater levels and anthropogenic and meteorological stressors at selected sites in east-central Florida, 1995-2007

Multivariate linear regression analyses were used to define the relations of water levels in the Upper Floridan aquifer (UFA) and surficial aquifer system (SAS) to anthropogenic and meteorological stressors between 1995 and 2007 at two monitoring well sites (Charlotte Street and Lake Oliver) in east-central Florida. Anthropogenic stressors of interest included municipal and agricultural groundwater withdrawals, and application of reclaimed-water to rapid-infiltration basins (source of aquifer recharge). Meteorological stressors included precipitation and potential evapotranspiration. Overall, anthropogenic and meteorological stressors accounted for about 40 to 89 percent of the variance in UFA and SAS groundwater levels and water-level changes. While mean monthly water levels were better correlated with monthly stressor values, changes in UFA and SAS water levels were better correlated with changes in stressor values. Water levels and water-level changes were influenced by system persistence as the moving-averaged values of both stressor types, which accounted for the influence of the previous month(s) conditions, consistently yielded higher adjusted coefficients of determination (R2 adj) values than did single monthly values. While monthly water-level changes tend to be influenced equally with both stressors across the hydrologically averaged 13-year period, changes were more influenced by one stressor or the other seasonally and during extended wet and dry periods. Seasonally, UFA water-level changes tended to be more influenced by anthropogenic stressors than by meteorological stressors, while changes in SAS water levels tended to be more influenced by meteorological stressors. During extended dry periods (12 months or greater), changes in UFA water levels at Charlotte Street were more affected by anthropogenic stressors than by meteorological stressors, while changes in SAS levels were more affected by meteorological stressors. At Lake Oliver, changes in both UFA and SAS water levels were better correlated with meteorological stressors for all but the wet period between April 1995 and April 1996. Interestingly, changes in both UFA and SAS water levels at Charlotte Street were also better correlated with anthropogenic stressors during a similar wet period between April 1995 and June 1996 when substantive reductions in groundwater withdrawals resulted in appreciable recovery of both UFA and SAS water levels. The regional effects of anthropogenic stressors had limited influence on water-level changes at Charlotte Street and virtually no influence on changes at Lake Oliver. When regressed against the 2.2 Mgal/d (million gallons per day) of municipal withdrawals located within 2 miles of the Charlotte Street site, water-level changes were influenced solely by precipitation and potential evapotranspiration. At a radius of 2.5 miles, however, where cumulative withdrawals totaled about 9.5 Mgal/d, water-level changes were equally influenced by both anthropogenic and meteorological stressors. Withdrawals located at distances of greater than 3 miles from this site had no appreciable effect on relations between water-level changes and these stressors. At Lake Oliver, changes in UFA water levels were equally influenced by both stressors regardless of distance, while changes in SAS levels were more influenced by meteorological stressors at all distances.

Scientific Investigations Report↗

USGS2011

Conceptual model of the Great Basin carbonate and alluvial aquifer system

A conceptual model of the Great Basin carbonate and alluvial aquifer system (GBCAAS) was developed by the U.S. Geological Survey (USGS) for a regional assessment of groundwater availability as part of a national water census. The study area is an expansion of a previous USGS Regional Aquifer Systems Analysis (RASA) study conducted during the 1980s and 1990s of the carbonate-rock province of the Great Basin. The geographic extent of the study area is 110,000 mi 2 , predominantly in eastern Nevada and western Utah, and includes 165 hydrographic areas (HAs) and 17 regional groundwater flow systems. A three-dimensional hydrogeologic framework was constructed that defines the physical geometry and rock types through which groundwater moves. The diverse sedimentary units of the GBCAAS study area are grouped into hydrogeologic units (HGUs) that are inferred to have reasonably distinct hydrologic properties due to their physical characteristics. These HGUs are commonly disrupted by large-magnitude offset thrust, strike-slip, and normal faults, and locally affected by caldera formation. The most permeable aquifer materials within the study area include Cenozoic unconsolidated sediments and volcanic rocks, along with Mesozoic and Paleozoic carbonate rocks. The framework was built by extracting and combining information from digital elevation models, geologic maps, cross sections, drill hole logs, existing hydrogeologic frameworks, and geophysical data.

Utah↗

USGS2012

Transient effects on groundwater chemical compositions from pumping of supply wells at the Nevada National Security Site, 1951-2008

Nuclear testing and support activities at the Nevada National Security Site have required large amounts of water for construction, public consumption, drilling, fire protection, hydraulic and nuclear testing, and dust control. To supply this demand, approximately 20,000 million gallons of water have been pumped from 23 wells completed in 19 boreholes located across the Nevada National Security Site starting as early as the 1950s. As a consequence of more or less continuous pumping from many of these wells for periods as long as 58 years, transient groundwater flow conditions have been created in the aquifers that supplied the water. To evaluate whether long-term pumping caused changes in water compositions over time, available chemical analyses of water samples from these 19 boreholes were compiled, screened, and evaluated for variability including statistically significant temporal trends that can be compared to records of groundwater pumping. Data used in this report have been extracted from a large database (Geochem08, revision 3.0, released in September 2008) containing geochemical and isotopic information created and maintained by primary contractors to the U.S. Department of Energy, National Nuclear Security Administration, Nevada Site Office. Data extracted from this source were compiled for the entire period of record, converted to uniform reporting units, and screened to eliminate analyses of poor or unknown quality, as well as clearly spurious values. The resulting data are included in accompanying spreadsheets that give values for (1) pH and specific conductance, (2) major ion concentrations, (3) trace element concentrations and environmental isotope ratios, and (4) mean, median, and variance estimates for major ion concentrations. The resulting data vary widely in quality and time-series density. An effort has been made to establish reasonable ranges of analytical uncertainty expected for each analyte and eliminate analyses that are obvious outliers. Analysis of chemical trends in this report primarily rely on specific conductance measurements and major ion concentrations, data considered to be the most accurate and reliable over the entire time span of investigation. The analysis uses parametric and non-parametric evaluations to provide a statistical basis for trend identification. Trace element and isotope data are examined for consistency, but typically are too inaccurate or infrequent to provide a reliable long term basis for trend evaluation. Groundwater withdrawal records compiled in a companion report are included graphically in this report to allow qualitative comparisons between water quality and pumping history. Data for each supply well include (1) a borehole description and summary of pumping history, (2) a description of water-quality parameters, (3) an evaluation of temporal variations of specific conductance and major ion concentrations, and (4) an examination of supporting information from trace element and isotope data. A range of responses are observed for individual supply wells that likely include the effects of both aquifer dynamics and changing borehole conditions. Data from most wells show little or no evidence for temporal variation in water-quality parameters indicating that aquifers at the Nevada National Security Site are capable of producing large volumes of compositionally uniform water over many years of pumping. A smaller number of wells show evidence of transient changes. Wells that have ceased pumping commonly show compositions that shift toward lower concentrations in subsequent bailed samples, which indicates that more dilute water entered the well over time due to either leakage of meteoric water into the well casing or more wide-spread recharge into the shallow phreatic zone. Wells that show systematic changes in water compositions during episodes of pumping commonly have multiple open intervals whose contributions to water in the well may change over time due to hydraulic conditions or well dynamics.

Nevada↗

USGS2012

Historical and potential groundwater drawdown in the Bruneau area, Owyhee County, southwestern Idaho

Geothermal seeps and springs in the Bruneau area in southwestern Idaho provide a vital but disappearing habitat for the Bruneau hot springsnail (Pyrgulopsis bruneauensis). In order to aid in conservation efforts, a two-part study was conducted (1) to determine trends in groundwater levels over time and (2) to simulate drawdown in aquifers that contribute to the geothermal seeps and springs along the Bruneau River. Seasonal and Regional Kendall tests for trends were used to determine water-level trends over a 20-year monitoring (1990–2010) period. Seasonal Kendall tests were used to calculate trends in groundwater-levels in 22 monitoring wells and indicated statistically significant changes in water level with trends ranging from 0.21 to 1.0 feet per year. Regional Kendall tests were used to calculate drawdown in categories of wells based on five criteria (well depth, distance from Indian Bathtub Spring, geologic unit, regional topographic valley, and temperature). Results from Regional Kendall tests indicate that slope of the trend (in feet per year) increased as a function of well depth; trends in water level as a function of other categories did not exhibit an obvious pattern based on distance from Indian Bathtub Spring, geologic unit, topographic valley, or temperature. Analytical solutions were used to simulate drawdown and recovery in wells using the Theis equation and a range of hydraulic parameters. Drawdown effects were determined by changing the storativity, transmissivity, and flow values over a hypothetical timeline. For example, estimates projected that after 20 years of pumping (at an assumed storativity of 0.002, a transmissivity of 980,000 feet squared per day, and a flow of 100 acre-feet per year), 1 foot of drawdown in the volcanic-rock aquifers would not be detected; however, other estimates using the same time frame but different hydraulic parameters (storativity of 0.001, transmissivity of 13,000 feet squared per day, and 610 acre-feet per year) determined 1 foot of drawdown to be detected as far as 29 miles from the hypothetical pumping well. A sensitivity analysis was performed to determine the effect of changing one hydraulic parameter while keeping the others constant. Many assumptions had to be made about properties of the aquifer in order to calculate effects of drawdown on geothermal seeps and springs. These analyses estimate pumping effects over time; the recovery of groundwater levels would likely take significantly longer to observe than the effects from pumping.

Idaho↗

USGS2013

Assessment of managed aquifer recharge at Sand Hollow Reservoir, Washington County, Utah, updated to conditions in 2012

Sand Hollow Reservoir in Washington County, Utah, was completed in March 2002 and is operated primarily for managed aquifer recharge by the Washington County Water Conservancy District. From 2002 through 2011, surface-water diversions of about 199,000 acre-feet to Sand Hollow Reservoir have allowed the reservoir to remain nearly full since 2006. Groundwater levels in monitoring wells near the reservoir rose through 2006 and have fluctuated more recently because of variations in reservoir altitude and nearby pumping from production wells. Between 2004 and 2011, a total of about 19,000 acre-feet of groundwater was withdrawn by these wells for municipal supply. In addition, a total of about 21,000 acre-feet of shallow seepage was captured by French drains adjacent to the North and West Dams and used for municipal supply, irrigation, or returned to the reservoir. From 2002 through 2011, about 106,000 acre-feet of water seeped beneath the reservoir to recharge the underlying Navajo Sandstone aquifer. Water quality was sampled at various monitoring wells in Sand Hollow to evaluate the timing and location of reservoir recharge as it moved through the aquifer. Tracers of reservoir recharge include major and minor dissolved inorganic ions, tritium, dissolved organic carbon, chlorofluorocarbons, sulfur hexafluoride, and noble gases. By 2012, this recharge arrived at four monitoring wells located within about 1,000 feet of the reservoir. Changing geochemical conditions at five other monitoring wells could indicate other processes, such as changing groundwater levels and mobilization of vadose-zone salts, rather than arrival of reservoir recharge.

Utah↗

USGS2013

Properties of basin-fill deposits, a 1971–2000 water budget, and surface-water-groundwater interactions in the upper Humboldt River basin, northeastern Nevada

This study was done in cooperation with Elko County, Nevada in response to concerns over growing demand for water within the county and increasing external demands that are occurring statewide. The upper Humboldt River basin encompasses 4,360 square miles in northeastern Nevada and includes the headwaters area of the Humboldt River. Nearly all of the mean annual flow of the Humboldt River originates in this area. Basin-fill deposits function as the principal aquifers in the upper Humboldt River basin. Over much of the basin lowlands, the upper 200 feet of basin fill consists of clay, silt, sand, and gravel deposited in a lake of middle to late Pliocene age. Fine-grained lacustrine sediments compose from 30 to more than 70 percent of the deposits. Mean values of transmissivity are less than 1,000 feet squared per day. Total inflow to the upper Humboldt River basin, about 3,330,000 acre-feet per year, is entirely from annual precipitation. Total outflow from the basin, about 3,330,000 acre-feet per year, occurs as evapotranspiration, streamflow, subsurface flow, and pumpage. The uncertainty of these values of inflow and outflow is estimated to be 25 percent. Baseflow of the Humboldt River is minimal upstream of the Elko Hills and in downstream reaches almost all baseflow comes from tributary inflow of the North Fork and South Fork Humboldt Rivers. However, the baseflow of these two tributaries comes from groundwater discharge to their respective channels in canyons incised in volcanic rocks along the North Fork and in carbonate rocks along the South Fork. Water levels in the shallow water-table aquifer along the Humboldt River flood plain fluctuate with changes in stage of the river. During high rising river stage in spring and early summer, streamflow enters the aquifer as bank storage. As stage begins to decline in early to mid-summer groundwater in bank storage begins discharging back into the river channel and this continues through late summer. In years of below average flow some reaches of the river are dry in late summer. Flood plain deposits are more permeable than adjacent and underlying fine-grained sediments of the Pliocene lake and the two aquifers are poorly connected.

Nevada↗

USGS2014

Analysis of potential water-supply management options, 2010-60, and documentation of revisions to the model of the Irwin Basin Aquifer System, Fort Irwin National Training Center, California

The Fort Irwin National Training Center is considering several alternatives to manage their limited water-supply sources in the Irwin Basin. An existing three-dimensional, finite-difference groundwater-flow model—the U.S. Geological Survey’s MODFLOW—of the aquifer system in the basin was updated and the initial input dataset was supplemented with groundwater withdrawal data for the period 2000–10. The updated model was then used to simulate four combinations, or scenarios, of groundwater withdrawal and recharge over the next 50 years (January 2011 through December 2060). The scenarios included combinations of continuing withdrawals from currently active production wells, supplementing any increases in demand with withdrawals from an inactive production well, reducing withdrawal amounts and rates, and reducing the discharge of treated wastewater to infiltration ponds that provide a recharge source to the underlying aquifer. Results of the simulations indicated that, depending on the scenario implemented, groundwater levels would rise (over the next 50 years) from 40 feet to as much as 65 feet in the northwestern part of the Irwin Basin, and from 5 feet to 10 feet in the southeastern part.

California↗

USGS2015

Hydrologic influences on water-level changes in the Eastern Snake River Plain aquifer at and near the Idaho National Laboratory, Idaho, 1949-2014

The U.S. Geological Survey, in cooperation with the U.S. Department of Energy, has maintained a water-level monitoring program at the Idaho National Laboratory (INL) since 1949 to systematically measure water levels to provide long-term information on groundwater recharge, discharge, movement, and storage in the eastern Snake River Plain (ESRP) aquifer. During 2014, water levels in the ESRP aquifer reached all-time lows for the period of record, prompting this study to assess the effect that future water-level declines may have on pumps and wells. Water-level data were compared with pump-setting depth to determine the hydraulic head above the current pump setting. Additionally, geophysical logs were examined to address changes in well productivity with water-level declines. Furthermore, hydrologic factors that affect water levels in different areas of the INL were evaluated to help understand why water-level changes occur. Review of pump intake placement and 2014 water-level data indicates that 40 wells completed within the ESRP aquifer at the INL have 20 feet (ft) or less of head above the pump. Nine of the these wells are located in the northeastern and northwestern areas of the INL where recharge is predominantly affected by irrigation, wet and dry cycles of precipitation, and flow in the Big Lost River. Water levels in northeastern and northwestern wells generally show water-level fluctuations of as much as 4.5 ft seasonally and show declines as much as 25 ft during the past 14 years. In the southeastern area of the INL, seven wells were identified as having less than 20 ft of water remaining above the pump. Most of the wells in the southeast show less decline over the period of record compared with wells in the northeast; the smaller declines are probably attributable to less groundwater withdrawal from pumping of wells for irrigation. In addition, most of the southeastern wells show only about a 1–2 ft fluctuation seasonally because they are less influenced by groundwater withdrawals for irrigation. In the southwestern area of the INL, 24 wells were identified as having less than 20 ft of water remaining above the pump. Wells in the southwest also only show small 1–2 ft fluctuations seasonally because of a lack of irrigation influence. Wells show larger fluctuation in water levels closer to the Big Lost River and fluctuate in response to wet and dry cycles of recharge to the Big Lost River. Geophysical logs indicate that most of the wells evaluated will maintain their current production until the water level declines to the depth of the pump. A few of the wells may become less productive once the water level gets to within about 5 ft from the top of the pump. Wells most susceptible to future drought cycles are those in the northeastern and northwestern areas of the INL.

Idaho↗

USGS2016

Statistical analysis and mapping of water levels in the Biscayne aquifer, water conservation areas, and Everglades National Park, Miami-Dade County, Florida, 2000–2009

Statistical analyses and maps representing mean, high, and low water-level conditions in the surface water and groundwater of Miami-Dade County were made by the U.S. Geological Survey, in cooperation with the Miami-Dade County Department of Regulatory and Economic Resources, to help inform decisions necessary for urban planning and development. Sixteen maps were created that show contours of (1) the mean of daily water levels at each site during October and May for the 2000–2009 water years; (2) the 25th, 50th, and 75th percentiles of the daily water levels at each site during October and May and for all months during 2000–2009; and (3) the differences between mean October and May water levels, as well as the differences in the percentiles of water levels for all months, between 1990–1999 and 2000–2009. The 80th, 90th, and 96th percentiles of the annual maximums of daily groundwater levels during 1974–2009 (a 35-year period) were computed to provide an indication of unusually high groundwater-level conditions. These maps and statistics provide a generalized understanding of the variations of water levels in the aquifer, rather than a survey of concurrent water levels. Water-level measurements from 473 sites in Miami-Dade County and surrounding counties were analyzed to generate statistical analyses. The monitored water levels included surface-water levels in canals and wetland areas and groundwater levels in the Biscayne aquifer. Maps were created by importing site coordinates, summary water-level statistics, and completeness of record statistics into a geographic information system, and by interpolating between water levels at monitoring sites in the canals and water levels along the coastline. Raster surfaces were created from these data by using the triangular irregular network interpolation method. The raster surfaces were contoured by using geographic information system software. These contours were imprecise in some areas because the software could not fully evaluate the hydrology given available information; therefore, contours were manually modified where necessary. The ability to evaluate differences in water levels between 1990–1999 and 2000–2009 is limited in some areas because most of the monitoring sites did not have 80 percent complete records for one or both of these periods. The quality of the analyses was limited by (1) deficiencies in spatial coverage; (2) the combination of pre- and post-construction water levels in areas where canals, levees, retention basins, detention basins, or water-control structures were installed or removed; (3) an inability to address the potential effects of the vertical hydraulic head gradient on water levels in wells of different depths; and (4) an inability to correct for the differences between daily water-level statistics. Contours are dashed in areas where the locations of contours have been approximated because of the uncertainty caused by these limitations. Although the ability of the maps to depict differences in water levels between 1990–1999 and 2000–2009 was limited by missing data, results indicate that near the coast water levels were generally higher in May during 2000–2009 than during 1990–1999; and that inland water levels were generally lower during 2000–2009 than during 1990–1999. Generally, the 25th, 50th, and 75th percentiles of water levels from all months were also higher near the coast and lower inland during 2000–2009 than during 1990–1999. Mean October water levels during 2000–2009 were generally higher than during 1990–1999 in much of western Miami-Dade County, but were lower in a large part of eastern Miami-Dade County.

Florida↗

USGS2024

Monitoring aquifer-storage change from artificial recharge with repeat microgravity along Santa Cruz River, Tucson, Arizona, 2019–22

The City of Tucson water utility, Tucson Water, began releasing treated effluent into the Santa Cruz River channel near downtown Tucson in 2019. This recharge project—the Heritage Project—is intended to create a reach of consistent flow in the channel and recharge water to the aquifer. Tracking the dispersal of recharged water is important for management decisions because groundwater movement depends on spatially variable characteristics of the subsurface and cannot be fully predicted in advance. Groundwater-level measurements in wells are useful, but the relation between water storage and groundwater-level change depends on the unknown storage coefficient of the aquifer. To estimate storage changes caused by recharge of reclaimed effluent released into the channel for the Heritage Project, the U.S. Geological Survey (USGS) collected repeat microgravity data along the Santa Cruz River in Tucson, Arizona, from 2019 to 2022. This method augments groundwater-level monitoring by providing a direct quantitative measurement of changes in the quantity of water stored in the subsurface. Preliminary results of the monitoring through 2022 showed consistent storage increases only near and upstream from the Heritage Project outfall site. Initially high storage increases at some locations west of the channel and in line with Sentinel Peak reached roughly steady state in later times. North of Sentinel Peak, a storage increase from 2020 to 2021 was followed by a storage decrease from 2021 to 2022. Storage changes in the area north of Sentinel Peak appear to be related to the number of days flows in the channel were observed farther downstream from the outfall site (at USGS streamgage 09482500). This observation is likely due to the potential formation of a clogging layer that would allow surface water to disperse farther horizontally (downstream) before infiltrating. This phenomenon has been observed downstream of other recharged effluent projects and has been reduced by large flows in the channel, such as those occurring during large runoff events. There were no large or consistent storage increases near the Water Quality Assurance Revolving Fund (WQARF) sites included in this study area.

Arizona↗

USGS2010

PHAST version 2-A program for simulating groundwater flow, solute transport, and multicomponent geochemical reactions

The computer program PHAST (PHREEQC And HST3D) simulates multicomponent, reactive solute transport in three-dimensional saturated groundwater flow systems. PHAST is a versatile groundwater flow and solute-transport simulator with capabilities to model a wide range of equilibrium and kinetic geochemical reactions. The flow and transport calculations are based on a modified version of HST3D that is restricted to constant fluid density and constant temperature. The geochemical reactions are simulated with the geochemical model PHREEQC, which is embedded in PHAST. Major enhancements in PHAST Version 2 allow spatial data to be defined in a combination of map and grid coordinate systems, independent of a specific model grid (without node-by-node input). At run time, aquifer properties are interpolated from the spatial data to the model grid; regridding requires only redefinition of the grid without modification of the spatial data. PHAST is applicable to the study of natural and contaminated groundwater systems at a variety of scales ranging from laboratory experiments to local and regional field scales. PHAST can be used in studies of migration of nutrients, inorganic and organic contaminants, and radionuclides; in projects such as aquifer storage and recovery or engineered remediation; and in investigations of the natural rock/water interactions in aquifers. PHAST is not appropriate for unsaturated-zone flow, multiphase flow, or density-dependent flow. A variety of boundary conditions are available in PHAST to simulate flow and transport, including specified-head, flux (specified-flux), and leaky (head-dependent) conditions, as well as the special cases of rivers, drains, and wells. Chemical reactions in PHAST include (1) homogeneous equilibria using an ion-association or Pitzer specific interaction thermodynamic model; (2) heterogeneous equilibria between the aqueous solution and minerals, ion exchange sites, surface complexation sites, solid solutions, and gases; and (3) kinetic reactions with rates that are a function of solution composition. The aqueous model (elements, chemical reactions, and equilibrium constants), minerals, exchangers, surfaces, gases, kinetic reactants, and rate expressions may be defined or modified by the user. A number of options are available to save results of simulations to output files. The data may be saved in three formats: a format suitable for viewing with a text editor; a format suitable for exporting to spreadsheets and postprocessing programs; and in Hierarchical Data Format (HDF), which is a compressed binary format. Data in the HDF file can be visualized on Windows computers with the program Model Viewer and extracted with the utility program PHASTHDF; both programs are distributed with PHAST.

Techniques and Methods↗

USGS2003

Simulated effects of ground-water management scenarios on the Santa Fe group aquifer system, Middle Rio Grande Basin, New Mexico, 2001-40

Future conditions in the Santa Fe Group aquifer system through 2040 were simulated using the most recent revision of the U.S. Geological Survey groundwater- flow model for the Middle Rio Grande Basin. Three simulations were performed to investigate the likely effects of different scenarios of future groundwater pumping by the City of Albuquerque on the ground-water system. For simulation I, pumping was held constant at known year-2000 rates. For simulation II, pumping was increased to simulate the use of pumping to meet all projected city water demand through 2040. For simulation III, pumpingwas reduced in accordance with a plan by the City of Albuquerque to use surfacewater to meet most of the projectedwater demand. The simulations indicate that for each of the three pumping scenarios, substantial additional watertable declines would occur in some areas of the basin through 2040. However, the reduced pumping scenario of simulation III also results in water-table rise over a broad area of the city. All three scenarios indicate that the contributions of aquifer storage and river leakage to the ground-water system would change between 2000 and 2040. Comparisons among the results for simulations I, II, and III indicate that the various pumping scenarios have substantially different effects on water-level declines in the Albuquerque area and on the contribution of each water-budget component to the total budget for the ground-water system. Between 2000 and 2040, water-level declines for continued pumping at year-2000 rates are as much as 120 feet greater than for reduced pumping; water-level declines for increased pumping to meet all projected city demand are as much as 160 feet greater. Over the same time period, reduced pumping results in retention in aquifer storage of about 1,536,000 acre-feet of ground water as compared with continued pumping at year- 2000 rates and of about 2,257,000 acre-feet as compared with increased pumping. The quantity of water retained in the Rio Grande as a result of reduced pumping and the associated decrease in induced recharge from the river is about 731,000 acre-feet as compared with continued pumping at year-2000 rates and about 872,000 acre-feet as compared with increased pumping. Reduced pumping results in slight increases in the quantity of water lost from the groundwater system to evapotranspiration and agriculturaldrain flow compared with the other pumping scenarios.

Water-Resources Investigations Report↗

USGS2003

Aquifer tests and simulation of ground-water flow in Triassic sedimentary rocks near Colmar, Bucks and Montgomery Counties, Pennsylvania

This report presents the results of a study by the U.S. Geological Survey in cooperation with the U.S. Environmental Protection Agency to evaluate ground-water flow in Triassic sedimentary rocks near Colmar, in Bucks and Montgomery Counties, Pa. The study was conducted to help the U.S. Environmental Protection Agency evaluate remediation alternatives at the North Penn Area 5 Superfund Site near Colmar, where ground water has been contaminated by volatile organic solvents (primarily trichloroethene). The investigation focused on determining the (1) drawdown caused by separately pumping North PennWater Authority wells NP–21 and NP–87, (2) probable paths of groundwater movement under present-day (2000) conditions (with NP–21 discontinued), and (3) areas contributing recharge to wells if pumping from wells NP-21 or NP–87 were restarted and new recovery wells were installed. Drawdown was calculated from water levels measured in observation wells during aquifer tests of NP–21 and NP–87. The direction of ground-water flow was estimated by use of a three-dimensional ground-water-flow model. Aquifer tests were conducted by pumping NP–21 for about 7 days at 257 gallons per minute in June 2000 and NP–87 for 3 days at 402 gallons per minute in May 2002. Drawdown was measured in 45 observation wells during the NP–21 test and 35 observation wells during the NP–87 test. Drawdown in observation wells ranged from 0 to 6.8 feet at the end of the NP–21 test and 0.5 to 12 feet at the end of the NP–87 test. The aquifer tests showed that ground-water levels declined mostly in observation wells that were completed in the geologic units penetrated by the pumped wells. Because the geologic units dip about 27 degrees to the northwest, shallow wells up dip to the southeast of the pumped well showed a good hydraulic connection to the geologic units stressed by pumping. Most observation wells down dip from the pumping well penetrated units higher in the stratigraphic section that were not well connected to the units stressed by pumping. The best hydraulic connection to the pumped wells was indicated by large drawdown in observation wells that penetrate the water-bearing unit encountered below 400 feet below land surface in wells NP–21 and NP–87. The hydraulic connection between wells NP–21 (or NP–87) and observation wells in the southern area of ground-water contamination near the BAE Systems facility is good because the observation wells probably penetrate this water-bearing unit. A 3-dimensional, finite-difference, groundwater- flow model was used to simulate flow paths and areas contributing recharge to wells for current (2000) conditions of pumping in the Colmar area and for hypothetical situations of pumping suggested by the U.S. Environmental Protection Agency that might be used for remediation. Simulations indicate that under current conditions, ground water in the northern area of contamination near the former Stabilus facility moves to the northwest and discharges mostly to West Branch Neshaminy Creek; in the southern area of contamination near BAE Systems facility, ground water probably moves west and discharges to a tributary of West Branch Neshaminy Creek near well NP–21. Model simulations indicate that if NP–21 or NP–87 are pumped at 400 gallons per minute, groundwater recharge is likely captured from the southern area of contamination, but ground-water recharge from the northern area of contamination is less likely to be captured by the pumping. Simulations also indicate that pumping of a new recovery well near BAE Systems facility at 8 gallons per minute and two new recovery wells near the former Stabilus facility at a total of about 30 gallons per minute probably would capture most of the ground-water recharge in the areas where contamination is greatest.

Water-Resources Investigations Report↗

USGS1984

Availability of water from the alluvial aquifer in part of the Green River Valley, King County, Washington

The Muckleshoot Indian Tribe plans (1982) to build a fish hatchery in part of a 1.56-square-mile area in the Green River valley, Washington, and use groundwater to operate it. Groundwater data were collected in the area and used in a U.S. Geological Survey two-dimensional groundwater-flow model calibrated to simulate the groundwater-flow system in the study area. Measured water levels in the alluvial aquifer were simulated to within 1 foot at 7 of 12 observation wells, and within 2 feet at all 12 wells. When pumping from the aquifer was simulated with the model, it was found that all water pumped from wells was derived from induced leakage from the Green River into the alluvium and reduced leakage through the alluvium to the Green River. Pumping from the alluvium may also reduce the flow of a tributary to the Green River.

Washington↗

USGS1985

Percentage change in saturated thickness of the High Plains aquifer, west-central Kansas, 1950 to average 1983-85

Continuing studies are being made in west-central Kansas to provide up-to-date information to aid in the management of groundwater for irrigation. This report, prepared in cooperation with the Western Kansas Groundwater Management District No. 1, presents the fifth in a series of studies that uses a statistical technique, called kriging, to produce hydrologic maps that are used as management tools. Kriging is a statistical technique that was used to interpolate water level altitudes at the center of each 1-square-mile section in the study area based on measured water levels at 165 observation wells. These interpolation altitudes (1,859 in all), along with bedrock surface and base year water table altitudes, were used to prepare a geohydrologic map illustrating percentage change in saturated thickness. Saturated thickness, as used in this report, is the thickness of the High Plains aquifer between the groundwater surface indicated by water table altitudes and the bedrock surface. Because irrigation development in west-central Kansas was minimal prior to 1950, the saturated thickness during 1950 represented a nearly static condition in the aquifer. Thus, the effects of irrigation withdrawalson the volume of water in storage could be related to the decrease or percentage change in saturated thickness of the aquifer from 1950 to average saturated thickness during 1983-85.

Kansas↗

USGS1987

Reconnaissance of water quality in the High Plains Aquifer beneath agricultural lands, south-central Kansas

The High Plains of western Kansas was one of 14 areas selected for preliminary groundwater quality reconnaissance by the U.S. Geological Survey 's Toxic Waste--Groundwater Contamination Program. The specific objective was to evaluate the effects of land used for agriculture (irrigated cropland and non-irrigated rangeland) on the water in the High Plains aquifer. Conceptual inferences, based on the information available, would lead one to expect groundwater beneath irrigated cropland to contain larger concentrations of sodium, sulfate, chloride, nitrite plus nitrate, and some water soluble pesticides than water beneath non-irrigated land (range-land) The central part of the Great Bend Prairie, an area of about 1,800 sq mi overlying the High Plains aquifer in south-central Kansas, was selected for the study of agricultural land use because it has sand soils, a shallow water table, relatively large annual precipitation, and includes large areas that are exclusively irrigated cropland or non-irrigated rangeland. As determined by a two-tailed Wilcoxon rank-sum test, concentrations of sodium and alkalinity were significantly larger at the 95% confidence level for water samples from beneath irrigated cropland than from beneath rangeland. No statistically significant difference in concentrations of sulfate, chloride, nitrite plus nitrate, and ammonia, was detected. Concentrations of 2,4-D found in water samples from beneath the rangeland were larger at the 99% confidence level as compared to concentrations of 2,4-D in samples from beneath irrigated cropland. Larger concentrations of sodium and alkalinity were found in water beneath irrigated cropland, and the largest concentration of the pesticide atrazine (triazines were found in three samples) was found in water from the only irrigation well sampled. The sodium and atrazine concentrations found in water from the irrigation well support the premise that water-level drawdown develops under irrigated fields. This diverts the natural groundwater flow patterns, so that pumpage may cause recycling and subsequent concentration of leachates from the land surface. (Author 's abstract)

Water-Resources Investigations Report↗