USGS Science⌕ Search

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.

At least 1,729 records · Page 96Linked to original sources

Simulation of groundwater flow at the former Badger Army Ammunition Plant, Sauk County, Wisconsin

To help support remedial efforts at the former Badger Army Ammunition Plant the U.S. Geological Survey built and calibrated a transient groundwater flow model using the Newton Raphson formulation (MODFLOW–NWT) of the U.S. Geological Survey’s modular three-dimensional finite-difference code. The model simulates the groundwater flow system at the site from 1984 to 2020. The former Badger Army Ammunition Plant is a 7,275-acre site in Sauk County, Wisconsin. The plant produced smokeless gunpower and solid rocket propellent as munitions components. Peak production periods were during World War II, the Korean War, and the Vietnam War. Subsequent groundwater contamination investigations have found four plumes at the site. A health risk assessment identified at least one contaminant of concern for human health risk present in three of the plumes: the propellant burning ground plume, the deterrent burning ground plume, and the central plume. A cooperative study began between the U.S. Army Environmental Command and U.S. Geological Survey to better understand the groundwater flow system at the former Badger Army Ammunition Plant. Field data, including aquifer tests, streamflow measurements, continuous groundwater elevations, and groundwater gradients with the Wisconsin River were collected and used to inform and calibrate the groundwater flow model. The model was used to assess the variability of the groundwater system over the study period, the components of the groundwater budget, and groundwater flow directions from identified source areas towards the Wisconsin River. Model performance assessment focused on using particle tracking to compare groundwater flowpaths that originate in the contaminant source areas to the observed plume footprints. This focus on plume behavior geometry should help constrain the advective component of a future groundwater transport model of the site.

Wisconsin↗

Hydrologic budgets for the Madison and Minnelusa aquifers, Black Hills of South Dakota and Wyoming, water years 1987-96

The Madison and Minnelusa aquifers are two of the most important aquifers in the Black Hills area of South Dakota and Wyoming. Quantification and evaluation of various hydrologic budget components are important for managing and understanding these aquifers. Hydrologic budgets are developed for two scenarios, including an overall budget for the entire study area and more detailed budgets for subareas. Budgets generally are combined for the Madison and Minnelusa aquifers because most budget components cannot be quantified individually for the aquifers. An average hydrologic budget for the entire study area is computed for water years 1987-96, for which change in storage is approximately equal to zero. Annual estimates of budget components are included in detailed budgets for nine subareas, which consider periods of decreasing storage (1987-92) and increasing storage (1993-96). Inflow components include recharge, leakage from adjacent aquifers, and ground-water inflows across the study area boundary. Outflows include springflow (headwater and artesian), well withdrawals, leakage to adjacent aquifers, and ground-water outflow across the study area boundary. Leakage, ground-water inflows, and ground-water outflows are difficult to quantify and cannot be distinguished from one another. Thus, net ground-water flow, which includes these components, is calculated as a residual, using estimates for the other budget components. For the overall budget for water years 1987-96, net ground-water outflow from the study area is computed as 100 ft 3 /s (cubic feet per second). Estimates of average combined budget components for the Madison and Minnelusa aquifers are: 395 ft 3 /s for recharge, 78 ft 3 /s for headwater springflow, 189 ft 3 /s for artesian springflow, and 28 ft 3 /s for well withdrawals. Hydrologic budgets also are quantified for nine subareas for periods of decreasing storage (1987-92) and increasing storage (1993-96), with changes in storage assumed equal but opposite. Common subareas are identified for the Madison and Minnelusa aquifers, and previous components from the overall budget generally are distributed over the subareas. Estimates of net ground-water flow for the two aquifers are computed, with net ground-water outflow exceeding inflow for most subareas. Outflows range from 5.9 ft 3 /s in the area east of Rapid City to 48.6 ft 3 /s along the southwestern flanks of the Black Hills. Net groundwater inflow exceeds outflow for two subareas where the discharge of large artesian springs exceeds estimated recharge within the subareas. More detailed subarea budgets also are developed, which include estimates of flow components for the individual aquifers at specific flow zones. The net outflows and inflows from the preliminary subarea budgets are used to estimate transmissivity of flow across specific flow zones based on Darcy’s Law. For estimation purposes, it is assumed that transmissivities of the Madison and Minnelusa aquifers are equal in any particular flow zone. The resulting transmissivity estimates range from 90 ft 2 /d to about 7,400 ft 2 /d, which is similar to values reported by previous investigators. The highest transmissivity estimates are for areas in the northern and southwestern parts of the study area, and the lowest transmissivity estimates are along the eastern study area boundary. Evaluation of subarea budgets provides confidence in budget components developed for the overall budget, especially regarding precipitation recharge, which is particularly difficult to estimate. Recharge estimates are consistently compatible with other budget components, including artesian springflow, which is a dominant component in many subareas. Calculated storage changes for subareas also are consistent with other budget components, specifically artesian springflow and net ground-water flow, and also are consistent with water-level fluctuations for observation wells. Ground-water budgets and flowpaths are especially complex in the southern Black Hills area; however, budget results are consistent with geochemical interpretations by previous investigators.

South Dakota, Wyoming↗

A model for calculating effects of liquid waste disposal in deep saline aquifer

Injection of liquid industrial wastes into confined underground saline aquifers can offer a good disposal alternative from both environmental and economic considerations. One of the needs in choosing from among several disposal alternatives is a means of evaluating the influence such an injection will have on the aquifer system. This report describes a mathematical model to accomplish this purpose. The objective of the contract was to develop a three-dimensional transient mathematical model which would accurately simulate behavior of waste injection into deep saline aquifers. Fluid properties, density and viscosity are functions of pressure, temperature and composition to provide a comprehensive assessment tool. The model is a finite-difference numerical solution of the partial differential equations describing single phase flow in the aquifer, energy transport by convection and conduction, and compositional changes in the aquifer fluid. The model is not restricted to examining waste disposal operations. It can be used effectively to evaluate fresh water storage in saline aquifers, hot water storage in underground aquifers, salt water intrusion into groundwater flow systems and other similar applications. The primary advantages of the present model can be summarized as: The model is user-oriented for easy application to full-scale evaluation needs. The model is fully three-dimensional and transient. The model is comprehensive accounting for density and viscosity variations in the aquifer due to temperature or compositional changes. The model includes the effects of hydrodynamic dispersion in both the temperature and compositional mixing between resident and injected fluids. The model energy balance includes the effects of pressure. This can be important in deep aquifer systems where the viscous pressure gradient is significant. The model uses second-order correct space and time approximations to the convective terms. This minimizes the numerical dispersion problem. The model is extremely flexible in providing a wide choice of boundary conditions. These include natural flow in the aquifer, aquifer influence functions around the perimeter of the grid in recognition that the gridded region does not have no-flow boundaries, heat losses into the overlying or underlying impermeable strata, and the wellbore heat and pressure drop calculations coupled to the aquifer flow equations. The limitations of the present techniques are: The use of the second-order correct finite-difference approximations introduces block size and time step restrictions. These restrictions, though considerably less stringent than explicit methods cause, depend upon the magnitude of the dispersivity. The comprehensive nature of the model makes the computer time and storage requirements significant. The model, because of its complexity, is not as efficient in reducing down to solve simpler problems as a specially written model would be. Included in the report are detailed descriptions of the approach used in the model, validation tests of the model, and a typical application of the model. A comparison volume documents the input data requirements, program structure, and an example problem for the model. '

Water-Resources Investigations Report↗

Hydrology of the Little Androscoggin River Valley aquifer, Oxford County, Maine

The Little Androscoggin River valley aquifer, a 15-square-mile sand and gravel valley-fill aquifer in southwestern Maine, is the source of water for the towns of Norway, Oxford, and South Paris. Estimated inflows to the aquifer during the 1981 water year were 16.4 cubic feet per second from precipitation directly on the aquifer, 11.2 cubic feet per second from till covered uplands adjacent to the aquifer, and 1.4 cubic feet per second from surface-water leakage. Outflows from the aquifer were 26.7 cubic feet per second to surface water and 2.3 cubic feet per second to wells. A finite-difference ground-water flow model was used to simulate conditions observed in the aquifer during 1981. Model conditions observed in the aquifer during 1981. Model simulations indicate that a 50 percent reduction of average 1981 recharge to the aquifer would cause water level declines of up to 20 feet in some areas. Model simulations of increased pumping at a high yield well in the northern part of the aquifer indicate that resulting changes in the water table will not be sufficient to intercept groundwater contaminated by a sludge disposal site. Water in the aquifer is low in dissolved solids (average for 38 samples was 67 mg/L), slightly acidic and soft. Ground-water contamination has occurred near a sludge-disposal site and in the vicinity of a sanitary landfill. Dissolved solids in ground water near the sludge disposal site were as much as ten times greater than average background values for the aquifer. (USGS)

Water-Resources Investigations Report↗

Geohydrology and water quality of the Inyan Kara, Minnelusa, and Madison aquifers of the northern Black Hills, South Dakota and Wyoming, and Bear Lodge Mountains, Wyoming

The Inyan Kara, Minnelusa, and Madison aquifers are the principal sources of ground water in the northern Black Hills, South Dakota and Wyoming, and Bear Lodge Mountains, Wyoming. The aquifers are exposed in the Bear Lodge Mountains and the Black Hills and are about 3,000 to 5,000 ft below the land surface in the northeast corner of the study area. The direction of groundwater movement is from the outcrop area toward central South Dakota. Recharge is by infiltration of precipitation and streamflow is by springs and well withdrawals. All three aquifers yield water to flowing wells in some part of the area. Measured and reported well yields in each of the three aquifers exceed 100 gal/min (gpm). A well open to the Minnelusa Formation and the upper part of the Madison Limestone yielded more than 2 ,000 gpm. Water from the Inyan Kara aquifer may require treatment for gross alpha radiation, iron, manganese, sulfate, and hardness before use in public water systems. Water from the Minnelusa aquifer in the northern one-half of the study area may require treatment for sulfate and hardness before use in public water systems. Water from the Madison aquifer in the northern one-half of the study area may require treatment of fluoride, gross alpha radiation, sulfate, and hardness before use in public water systems. Water from the Minnelusa and Madison aquifers in the southern one-half of the study area, though very hard (more than 180 mg/L hardness as calcium carbonate), is suitable for public water systems and irrigation. Flow between the Minnelusa and the Inyan Kara aquifers appears to be insignificant, based on the results of a digital model results. The model indicated there may be significant recharge to the Minnelusa and Madison aquifers by leakage between these two aquifers and perhaps deeper aquifers. (Author 's abstract)

Water-Resources Investigations Report↗

Aquifers in Cretaceous rocks of the central Coastal Plain of North Carolina

Aquifers in rocks of Cretaceous age are the major source of groundwater for public supplies in the central Coastal Plain. These aquifers consist of sand, gravel, and limestone beds of the Peedee, Black Creek, and the upper and lower Cape Fear aquifers, each separated by a confining unit composed of clay and silt beds. The aquifers and confining units (1) rest upon crystalline basement rocks; (2) dip and thicken to the east-southeast; and (3) are overlain by younger aquifers and confining units in deposits of Quaternary and Tertiary age. The top of the uppermost aquifer, the Peedee, ranges from 122 ft above to 595 ft below sea level. The maximum thickness of all aquifers and confining units in Cretaceous rocks is more than 1,600 ft. Aquifers and confining units were defined and correlated for this report using 125 geophysical logs and accompanying drillers' logs, water level data, and water quality information regarding chloride concentrations in water. This analysis allowed the construction of seven hydrogeologic cross sections that depict the continuity of all the aquifers and confining units. These cross sections also show water levels and chloride concentrations in water from various test intervals and describe where chloride concentration in water exceeds 250 mg/L within each aquifer. Detailed maps of each Cretaceous aquifer show altitude of its top, thickness, areas of selected sand percentages, transition zones from fresh to saltwater, and the thickness of the confining unit. H ydrogeologic data for all aquifers and confining units are presented in tabular form.

North Carolina↗

Geohydrology and simulated effects of large ground-water withdrawals on the Mississippi River alluvial aquifer in northwestern Mississippi

The 7,000-square-mile Mississippi River alluvial plain in northwestern Mississippi, locally known as the 'Delta,' is underlain by a prolific aquifer that yielded about 1,100 million gallons per day of water to irrigation wells in 1983. About 20 feet of clay underlying the Delta land surface commonly is underlain by about 80 to 180 feet of sand and gravel that forms the Mississippi River alluvial aquifer. This study of the alluvial aquifer was prompted by recent declines of water levels. The study was designed to better define the hydrology of the aquifer and to quantify availability of water from the aquifer. The Mississippi River is in good hydraulic connection with the alluvial aquifer. Generally, smaller streams are less likely to recharge the aquifer than larger streams. Direct vertical recharge to the alluvial aquifer from the 52 inches per year of precipitation is small, especially in the central part of the Delta. A two-dimensional finite-difference computer model of the alluvial aquifer was constructed, calibrated, and verified using water levels observed for five dates from April 1981 to September 1983. The values of some of the calibration-derived parameters are hydraulic conductivity, 400 feet per day; specific yield, 0.30; and infiltration of precipitation to the aquifer, 0.5 inch per year. The model showed that the aquifer had a net loss in storage of about 360 million gallons per day from April 1981 to April 1983. During this period, pumpage was about 1,100 million gallons per day (1,270,000 acre-feet per year), and the net inflows from the sources of recharge were as follows, in million gallons per day: Mississippi River, 390; recharge along the east edge of the Delta, 170; streams within the Delta, 57; areal recharge from infiltration, 180; and oxbow lakes, 24. The effects of several levels of pumpage by wells--0, 670, 1,100, 1,900, and 4,000 million gallons per day--were projected 20 years into the future. In 2003, the 1,100-milliongallon-per-day pumping rate, about average for the early 1980's, would take 46 percent of the water withdrawn from aquifer storage, water levels would be lowered more than 20 feet in a large area in the central part of the Delta, and groundwater levels would continue to decline in future years.

Water Supply Paper↗

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗

Quantifying effects of climate change on the snowmelt-dominated groundwater resources of northern New England

Recent U.S. Geological Survey (USGS) climate studies in New England have shown substantial evidence of hydrologic changes during the last 100 years, including trends toward earlier snowmelt runoff, decreasing occurrence of river ice, and decreasing winter snowpack. These studies are being expanded to include investigation of trends in groundwater levels and fluctuations. Groundwater is an important drinking-water source throughout northern New England (Maine, New Hampshire, and Vermont). The USGS is currently investigating whether or not groundwater recharge from snowmelt and precipitation exhibits historical trends. In addition to trend-testing, groundwater resources also will be analyzed by relating groundwater-level changes to the large year-to-year variability in weather conditions. Introduction The USGS has documented many seasonal climate-related changes in the northeastern United States that have occurred during the last 30 to 150 years. These changes include earlier snowmelt runoff in the late winter and early spring, decreasing duration of ice on rivers and lakes, decreasing ratio of snowfall to total precipitation, and denser and thinner late-winter snowpack. All of these changes are consistent with warming winter and spring air temperatures (Dudley and Hodgkins, 2002; Hodgkins and others, 2002; Huntington and others, 2004; Hodgkins and others, 2005; Hodgkins and Dudley, 2006a; Hodgkins and Dudley, 2006b). Climate-model projections for the Northeast indicate air-temperature warming, earlier snowmelt runoff, increases in annual evaporation, and decreased low streamflows (Hayhoe and others, 2007). The contribution and timing of spring snowmelt to groundwater recharge is particularly important to groundwater resources in the northeastern United States where aquifers typically consist of thin sediments overlying crystalline bedrock with relatively little storage capacity (Mack, 2009). Following spring recharge, groundwater slowly flows into streams throughout the summer. This groundwater flow is a source of cool water during the summer and accounts for a large proportion of the streamflow during summer low-flow periods. Groundwater is an important drinking-water source in northern New England. Approximately 32 percent of public water suppliers draw water from groundwater sources in Vermont, New Hampshire, and Maine, and approximately 40 percent of the population derives its drinking water from private wells (Kenny and others, 2009). It is vital to understand changes that may be occurring to such an important resource for planning industrial and agricultural water uses and protecting drinking water.

Fact Sheet↗

Changes in water levels and storage in the High Plains Aquifer, predevelopment to 2009

The High Plains aquifer underlies 111.8 million acres (175,000 square miles) in parts of eight States - Colorado, Kansas, Nebraska, New Mexico, Oklahoma, South Dakota, Texas, and Wyoming. The area overlying the High Plains aquifer is one of the primary agricultural regions in the Nation. Water-level declines began in parts of the High Plains aquifer soon after the onset of substantial irrigation with groundwater from the aquifer (about 1950 and termed "predevelopment" in this fact sheet). By 1980, water levels in the High Plains aquifer in parts of Texas, Oklahoma, and southwestern Kansas had declined more than 100 feet (ft) (Luckey and others, 1981). In 1987, in response to declining water levels, Congress directed the U.S. Geological Survey (USGS), in collaboration with numerous Federal, State, and local water-resources entities, to assess and track water-level changes in the aquifer. This fact sheet summarizes changes in water levels and drainable water in storage in the High Plains aquifer from predevelopment to 2009. Drainable water in storage is the fraction of water in the aquifer that will drain by gravity and can be withdrawn by wells. The remaining water in the aquifer is held to the aquifer material by capillary forces and generally cannot be withdrawn by wells. Drainable water in storage is termed "water in storage" in this report. A companion USGS report presents more detailed and technical information about water-level and storage changes in the High Plains aquifer during this period (McGuire, 2011).

Colorado;Kansas;Nebraska;New Mexico;Oklahoma;South↗

Building groundwater modeling capacity in Mongolia

Ulaanbaatar, the capital city of Mongolia (fig. 1), is dependent on groundwater for its municipal and industrial water supply. The population of Mongolia is about 3 million people, with about one-half the population residing in or near Ulaanbaatar (World Population Review, 2016). Groundwater is drawn from a network of shallow wells in an alluvial aquifer along the Tuul River. Evidence indicates that current water use may not be sustainable from existing water sources, especially when factoring the projected water demand from a rapidly growing urban population (Ministry of Environment and Green Development, 2013). In response, the Government of Mongolia Ministry of Environment, Green Development, and Tourism (MEGDT) and the Freshwater Institute, Mongolia, requested technical assistance on groundwater modeling through the U.S. Army Corps of Engineers (USACE) to the U.S. Geological Survey (USGS). Scientists from the USGS and USACE provided two workshops in 2015 to Mongolian hydrology experts on basic principles of groundwater modeling using the USGS groundwater modeling program MODFLOW-2005 (Harbaugh, 2005). The purpose of the workshops was to bring together representatives from the Government of Mongolia, local universities, technical experts, and other key stakeholders to build in-country capacity in hydrogeology and groundwater modeling. A preliminary steady-state groundwater-flow model was developed as part of the workshops to demonstrate groundwater modeling techniques to simulate groundwater conditions in alluvial deposits along the Tuul River in the vicinity of Ulaanbaatar. ModelMuse (Winston, 2009) was used as the graphical user interface for MODFLOW for training purposes during the workshops. Basic and advanced groundwater modeling concepts included in the workshops were groundwater principles; estimating hydraulic properties; developing model grids, data sets, and MODFLOW input files; and viewing and evaluating MODFLOW output files. A key to success was developing in-country technical capacity and partnerships with the Mongolian University of Science and Technology; Freshwater Institute, Mongolia, a non-profit organization; United Nations Educational, Scientific and Cultural Organization (UNESCO); the Government of Mongolia; and the USACE.

Open-File Report↗

Hydrologic investigations and a preliminary conceptual model of the groundwater system at North Penn Area 1 Superfund Site, Souderton, Montgomery County, Pennsylvania

The U.S. Geological Survey (USGS) conducted hydrogeologic investigations, reviewed existing data, and developed a preliminary conceptual model of the groundwater system as part of technical support of the U.S. Environmental Protection Agency (EPA) at the North Penn Area 1 Superfund Site (hereafter, the NP1 Site) located within the Borough of Souderton in Montgomery County, Pennsylvania. Field work and monitoring took place during 2012–18. The area is underlain by sedimentary formations that form a fractured-rock aquifer used for drinking water and industrial supply. The EPA placed the Site on the National Priorities List in 1989, identifying tetrachloroethylene (PCE) and trichloroethylene (TCE) as contaminants of concern. During 2012–18, the USGS conducted field activities that included drilling an 82-foot (ft)-deep monitoring well (MG 2220) in 2016, reconstructing a 208-ft-deep former industrial production well (MG 668 [Granite Knitting Mill]), and collecting borehole geophysical and video logs and water levels from those and five additional wells, which ranged in depth from about 50 to 200 ft below land surface. Continuous water levels were collected during 2014–17, and a synoptic set of water levels were measured in April 2018 in the seven wells. The borehole geophysical logs (caliper, acoustic televiewer, natural gamma, single-point resistance, vertical flow, and fluid temperature and resistivity) and borehole video logs in the seven wells were evaluated to assess potential for lithologic correlation and to identify and describe water-bearing features, which included both low- and high-angle fractures and other openings oriented along dipping bedding planes, joints, or possible faults. Borehole geophysical logs collected by USGS in 1992 in a 300-ft-deep former production well near the Site were also evaluated. Few to no distinctive features were identified on geophysical logs (natural gamma and single-point resistance) that could be used for correlation, thus limiting this approach to determining local geologic structure. Extensive fracturing in the upper 62 ft of monitoring well MG 2220 indicates that the well was likely drilled through a zone of faulting, and other evidence of faulting is present in the area near the Site. Assessment of continuous water levels showed hydraulic connections among some wells as indicated by rising or falling water levels in response to changes in pumping rates at nearby wells. A map of water levels measured in April 2018 indicates potential for groundwater flow generally toward the stream to the south and southwest of the Site, but the limited water-level data are insufficient to describe vertical groundwater gradients or lateral gradients in any detail. Review of 1999–2022 volatile organic compound (VOC) monitoring data collected by the Pennsylvania Department of Environmental Protection for five monitoring wells indicates that the highest groundwater concentrations of PCE and TCE were found in samples from extraction well MG 2201 (S-1) downgradient from, and nearest to, the previously identified Site contaminant source area, and these concentrations fluctuated through time. PCE concentrations were higher than TCE concentrations in samples from all five monitoring wells and were much higher than TCE concentrations in samples from extraction well MG 2201 (S-1). Temporally variable recharge is a possible factor affecting observed fluctuations in PCE concentrations in groundwater samples from well extraction MG 2201 (S-1), as indicated by a general inverse relation between PCE concentrations and water levels in a nearby long-term observation well. The PCE concentration of 1,830 micrograms per liter (μg/L) in a May 2018 water sample from monitoring well MG 2220 was more than four times the PCE concentration of 444 μg/L in a December 2017 sample from the nearby extraction well MG 2201 (S-1), which is open to fewer fractures. Low concentrations of VOCs were measured in surface water at two stream sites downgradient from wells with the highest groundwater VOC concentrations at the Site, indicating that discharge of contaminated groundwater to the stream is likely. Development of a conceptual model of the groundwater system was constrained by limited data. In areas with no pumping, groundwater-flow directions generally are thought to be controlled by topography and geologic structure (bedding orientation) and likely to the south and southwest of the Site, with local flow directions affected by orientations of fractures, joints, and local faults. Additional investigations that could help improve the conceptual model of the groundwater system and help delineate the extent of groundwater contamination and its transport are discussed.

Pennsylvania↗

Appraisal of water in bedrock aquifers, northern Cascade County, Montana

Suburban residential expansion of the city of Great Falls has resulted in an increased demand on water supplies from bedrock aquifers in northern Cascade County. The unconsolidated deposits aquifer of Quaternary age, including alluvium and glacial lake deposits, also is an important source of water in the area. Water levels in the Madison-Swift aquifer and all overlying aquifers, including the Quaternary deposits aquifer, reflect unconfined (water-table) conditions in the Great Falls vicinity. This interconnected hydrologic system is the result of breaching of the major anticlinal structure, by ancestral and present day erosion of drainage channels by the Missouri River and its tributaries. Significant vertical inter-aquifer mixing of water, as well as surface water/groundwater interchange, probably occurs in the central part of the study area. Characterization of the chemical composition of water in individual aquifers based on samples from wells in this area probably is unreliable because of this mixing. Quality of water from two wells in the Madison-Swift aquifer near Giant Springs is similar to water from the springs. Water from these three samples is less mineralized than most groundwater in the study area; dissolved solids concentrations for the three samples range from 516 to 550 mg/L. The quality of water varies among aquifers and throughout the study area. The ranges of dissolved solids concentrations determined by chemical analysis are Madison-Swift aquifer, about 520 to 1,570 mg/L; Morrison Formation, 908 to 1 ,480 mg/L; Kootenai Formation, 558 to 1,550 mg/L; Colorado Group , 2,690 and 2,740 mg/L (two samples); and unconsolidated Quaternary deposits, 383 to 2,060 mg/L. The chemical quality of water from the Colorado Group in the western one-third of the area generally is more mineralized than water from aquifers in the rest of the area. Specific conductance of water from eight wells completed in the Colorado Group averages 4,440 micromhos at 25 C. (Author 's abstract)

Open-File Report↗

Water quality of groundwater and stream base flow in the Marcellus Shale Gas Field of the Monongahela River Basin, West Virginia, 2011-12

The Marcellus Shale gas field underlies portions of New York, Pennsylvania, Ohio, Virginia, Maryland, Tennessee, and West Virginia. Development of hydraulic fracturing and horizontal drilling technology led to extensive development of gas from the Marcellus Shale beginning about 2007. The need to identify and monitor changes in water-quality conditions related to development of the Marcellus Shale gas field prompted the U.S. Geological Survey, in cooperation with the West Virginia Department of Environmental Protection, Division of Water and Waste Management, to document water quality for comparison with water quality in samples collected at a future date. The identification of change in water-quality conditions over time is more difficult if baseline water-quality conditions have not been documented. U.S. Geological Survey personnel sampled groundwater and surface water in West Virginia’s Monongahela River Basin during 2011–12. A groundwater survey, in which 39 wells and 2 springs were sampled, was conducted during June through September 2011. A base-flow survey was conducted during July through October 2012; 50 stream sites were sampled under base-flow conditions in this survey. Because additives to hydraulic fracturing fluids are variable and decrease in flowback water over a relatively short time, water-quality analyses for this study focused on documenting the water-quality characteristics typical of water from shallow aquifers; water derived from contact with the Marcellus Shale (flowback from hydraulic fracturing or formation water); and water with constituents from conventional oil and gas development, sewage effluent, and coal-mine drainage. All samples were analyzed for field properties (water temperature, pH, specific conductance, dissolved oxygen, and turbidity), major ions, trace elements, naturally occurring radioactive materials, and stable isotopes. In addition to documenting baseline water-quality conditions for an area of shale-gas development, these data were examined for patterns in water quality. Groundwater and base-flow survey data were compared to historical data from the Monongahela River Basin in West Virginia. Additionally, groundwater- and base-flow survey samples were grouped by Marcellus Shale gas production in the subbasin in which that sampling site was located. The comparisons of data collected as part of this study with historical data identified few differences. No significant difference was found in a comparison of groundwater survey data and historical data. Base-flow survey samples differed significantly from historical data for pH, chloride, and strontium, all of which had higher concentrations in the base-flow survey samples. Differences in pH are likely related to changes in mining regulation beginning in 1977. Concentrations of chloride and strontium elevated above background concentrations may be related to saline groundwater; saline water is within 300 feet of the land surface in parts of the study area. In the comparison of base-flow survey samples grouped by shale-gas-production setting, significant differences were found for fluoride and barium. Concentrations of fluoride and barium were higher in stream subbasins with active Marcellus Shale production than in subbasins not near active Marcellus Shale production. Elevated fluoride and barium are associated with deep brines. Generally, naturally occurring radioactive materials were not found in elevated concentrations in either groundwater or base-flow samples. Only 3 samples, 2 from the groundwater survey and one from the base-flow survey, exceeded the U.S. Environmental Protection Agency maximum contaminant level for radium isotopes of 5.0 picocurie per liter for either a single isotope or a combined value of radium-226 and radium-228. Stable isotope composition indicates broad similarity among surface water, shallow groundwater, and precipitation in the region. Neither shallow groundwater nor surface water showed a marked similarity with the deep brines associated with shale gas. In most of the groundwater survey samples, 38 of 41 samples, dissolved gas profiles were similar to those previously found in samples from shallow, domestic wells in the region. This study provides a baseline of water-quality conditions in the Monongahela River Basin in West Virginia during the early phases of development of the Marcellus Shale gas field. Although not all inclusive, the results of this study provide a set of reliable water-quality data against which future data sets can be compared and the effects of shale-gas development may be determined.

West Virginia↗

An integrated hydrologic model to support the Central Platte Natural Resources District Groundwater Management Plan, central Nebraska

The groundwater and surface-water supply of the Central Platte Natural Resources District supports a large agricultural economy from the High Plains aquifer and Platte River, respectively. This study provided the Central Platte Natural Resources District with an advanced numerical modeling tool to assist with the update of their Groundwater Management Plan. An integrated hydrologic model, called the Central Platte Integrated Hydrologic Model, was constructed using the MODFLOW-One-Water Hydrologic Model code with the Newton solver. This code integrates climate, landscape, surface water, and groundwater-flow processes in a fully coupled approach. Model framework included 163 rows; 327 columns; 2,640 feet cell sides; and 3 vertical layers. A predevelopment model simulated steady-state hydrologic conditions prior to April 30, 1895, and a development period model discretized into 610 stress periods simulated transient hydrologic conditions from May 1, 1895, to December 31, 2016, using 170 biannual stress periods from 1895 to 1980, and monthly stress periods from May 1, 1980, to December 31, 2016. Calibration of the Central Platte Integrated Hydrologic Model involved two phases: a manual adjustment of parameters, followed by the automated calibration completed using BeoPEST that was facilitated by the employment of the singular value decomposition-assist features of PEST that specified 50 super parameters assembled from the 435 adjustable parameters and Tikhonov regularization. The average absolute groundwater-level residuals for model layers one, two, and three were 6.1, 12.4, and 7.4 feet, respectively. Calibrated horizontal hydraulic conductivity was about 70, 32, and 35 feet per day for layers 1, 2, and 3, respectively. The largest development period inflow to groundwater was recharge from deep percolation past the root zone, averaging 1,122,257 acre-feet per year (2.7 inches per year), and the largest outflow was to irrigation wells, averaging 693,171 acre-feet per year (10.2 inches per year for the Central Platte Natural Resources District). Other substantial groundwater outflows included evapotranspiration and base flow. For the total development period, there was a net change in storage of −122,393 acre-feet per year (−0.3 inch per year). The calibrated Central Platte Integrated Hydrologic Model was used to simulate eight different potential future climate and irrigation pumping conditions from January 1, 2017, to December 31, 2049. Simulated future groundwater levels within the Central Platte Natural Resources District varied significantly between scenarios and locally, from 13.8 feet below to 7.6 feet above baseline 1982 groundwater levels. Most areas exhibited groundwater-level declines for the drought scenarios and rises for the alternate irrigation scenarios. Changes in scenario groundwater levels correlated with the relations between farm net recharge and irrigation pumping. Linear “first order second moment” techniques indicated that the uncertainty in projected groundwater altitudes was reduced by 15.33 feet through model calibration.

Nebraska↗

Simulation of projected water demand and ground-water levels in the Coffee Sand and Eutaw-McShan aquifers in Union County, Mississippi, 2010 through 2050

Ground water from the Eutaw-McShan and the Coffee Sand aquifers is the major source of supply for residential, commercial, and industrial purposes in Union County, Mississippi. Unbiased, scientifically sound data and assessments are needed to assist agencies in better understanding and managing available water resources as continuing development and growth places more stress on available resources. The U.S. Geological Survey, in cooperation with the Tennessee Valley Authority, conducted an investigation using water-demand and ground-water models to evaluate the effect of future water demand on groundwater levels. Data collected for the 12 public-supply facilities and the self-supplied commercial and industrial facilities in Union County were used to construct water-demand models. The estimates of water demand to year 2050 were then input to a ground-water model based on the U.S. Geological Survey finite-difference computer code, MODFLOW. Total ground-water withdrawals for Union County in 1998 were estimated as 2.85 million gallons per day (Mgal/d). Of that amount, municipal withdrawals were 2.55 Mgal/d with about 1.50 Mgal/d (59 percent) delivered to residential users. Nonmunicipal withdrawals were 0.296 Mgal/d. About 80 percent (2.27 Mgal/d) of the total ground-water withdrawal is produced from the Eutaw-McShan aquifer and about 13 percent (0.371 Mgal/d) from the Coffee Sand aquifer. Between normal- and high-growth conditions, total water demand could increase from 72 to 131 percent (2.9 Mgal/d in 1998 to 6.7 Mgal/d in year 2050) with municipal demand increasing from 77 to 146 percent (2.6 to 6.4 Mgal/d). Increased pumping to meet the demand for water was simulated to determine the effect on water levels in the Coffee Sand and Eutaw- McShan aquifers. Under baseline-growth conditions, increased water use by year 2050 could result in an additional 65 feet of drawdown in the New Albany area below year 2000 water levels in the Coffee Sand aquifer and about 120 feet of maximum drawdown in the Eutaw-McShan aquifer. Under normal-growth conditions, increased water use could result in an additional 65 feet of drawdown in the New Albany area below year 2000 water levels in the Coffee Sand aquifer and about 135 feet of maximum drawdown in the Eutaw-McShan aquifer. Under high-growth conditions, increased water use could result in 75 feet of drawdown in the New Albany area below year 2000 water levels in the Coffee Sand aquifer and about 190 feet of maximum drawdown in the Eutaw-McShan aquifer. The resulting highgrowth projected water level for the year 2050 at the center of the drawdown cone in the New Albany area is between 450 and 500 feet above the top of the Eutaw-McShan aquifer.

Mississippi↗

Geohydrology of the High Plains Aquifer, western Kansas

The High Plains aquifer underlies 174,050 sq mi of eight states (Colorado, Kansas, Nebraska, New Mexico, Oklahoma, South Dakota, Texas, and Wyoming) and contains approximately 3.3 billion acre-ft of water in storage. Saturated thicknesses within the aquifer are as great as 600 ft near the southern border of southwest Kansas. The aquifer is replenished primarily by infiltration from precipitation. Average precipitation at the Garden City Experiment Station is 18.93 in/yr. Groundwater flow is generally from west to east under unconfined conditions. Hydraulic connection with subcropping consolidated aquifers allows ground water to flow vertically in minor quantities. The aquifer is depleted primarily by irrigation. Hydraulic conductivity estimates from 1,612 lithologic logs had an average value of 75 ft/day, with a standard deviation of 35 ft/day. Specific yields estimated from the same lithologic logs had a mean of 0.17 and a standard deviation of 0.047. Water from the High Plains aquifer in Kansas generally is suitable for human and animal consumption and irrigation of crops. Typically, it is a calcium bicarbonate type water, with concentrations of total dissolved solids ranging from 250 to 500 mg/L. The quality of water in the aquifer deteriorates toward the east due to mixing with recharge water containing dissolved minerals leached from the overlying soil and unsaturated zones and mineralized water from adjacent bedrock units. The simulated water budget for the steady state model of predevelopment (pre-1950) conditions in the High Plains aquifer in northwest Kansas showed that annual recharge to the aquifer from infiltration of precipitation was 87,000 acre-ft/yr and from boundary inflow, 21,000 acre-ft/yr. Annual discharge from the aquifer was 108,000 acre-ft/yr, including 81,000 acre-ft/yr from leakage to streams, 23,000 acre-ft from outflow at the boundaries of the aquifer, and 4,000 acre-ft from municipal and industrial pumpage. (Lantz-PTT)

Water-Resources Investigations Report↗

Simulation of the flow system of Barton Springs and associated Edwards Aquifer in the Austin area, Texas

A digital model of two-dimensional groundwater flow was used to estimate the hydraulic properties of the Edwards Aquifer in a 151 sq mi area near Austin, Texas. The transmissivity, hydraulic conductivity, and specific yield were estimated for the part of the aquifer that discharges at Barton Springs in Austin. The aquifer is composed of the Edwards and overlying Georgetown Limestones of Cretaceous age and ranges in thickness from about 100 ft to about 450 ft. More than 60 years of discharge measurements and 5 years of gaged discharge for Barton Springs were used to adjust springflow for the simulations. Barton Springs accounts for about 96% of springflow from the study area and 90% of the total discharge. The remaining discharge was pumpage from wells which was entered in the model. Four years of gaged recharge were used in the simulations. The potentiometric surfaces used by the models were constructed from water level measurements in as many as 75 wells. The transmissivity was calibrated through steady-state simulations that used the mean value of recharge and mean potentiometric surface to represent average conditions for the aquifer. The transmissivities vary from about 100 sq ft/day in the western part of the aquifer to > 1 million sq ft/day near Barton Springs. Specific yield was calibrated through transient-state simulations for 5 consecutive months using time-dependent data for recharge, discharge, and water levels. The mean specific yield for the aquifer is 0.014 and ranges from 0.008 to 0.064. Additional aquifer properties used in the simulations include storage coefficient, altitudes of the base and top of the aquifer, and hydraulic conductivity. A simulation for the year 2000 using projected pumping rates for municipal, industrial, agricultural, and domestic supplies indicates that the aquifer would be dewatered in the southwestern part of the study area and have large declines in the southeastern part of the study area. Another simulation of projected conditions using potential recharge enhancement predicts a rise in the potentiometric surface of about 50 feet in the southwestern part of the aquifer and moderate water-level declines in the southeastern part of the aquifer.

Texas↗