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,153 records · Page 64Linked to original sources

USGS2009

Numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area, South Dakota

The city of Rapid City and other water users in the Rapid City area obtain water supplies from the Minnelusa and Madison aquifers, which are contained in the Minnelusa and Madison hydrogeologic units. A numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area was developed to synthesize estimates of water-budget components and hydraulic properties, and to provide a tool to analyze the effect of additional stress on water-level altitudes within the aquifers and on discharge to springs. This report, prepared in cooperation with the city of Rapid City, documents a numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units for the 1,000-square-mile study area that includes Rapid City and the surrounding area. Water-table conditions generally exist in outcrop areas of the Minnelusa and Madison hydrogeologic units, which form generally concentric rings that surround the Precambrian core of the uplifted Black Hills. Confined conditions exist east of the water-table areas in the study area. The Minnelusa hydrogeologic unit is 375 to 800 feet (ft) thick in the study area with the more permeable upper part containing predominantly sandstone and the less permeable lower part containing more shale and limestone than the upper part. Shale units in the lower part generally impede flow between the Minnelusa hydrogeologic unit and the underlying Madison hydrogeologic unit; however, fracturing and weathering may result in hydraulic connections in some areas. The Madison hydrogeologic unit is composed of limestone and dolomite that is about 250 to 610 ft thick in the study area, and the upper part contains substantial secondary permeability from solution openings and fractures. Recharge to the Minnelusa and Madison hydrogeologic units is from streamflow loss where streams cross the outcrop and from infiltration of precipitation on the outcrops (areal recharge). MODFLOW-2000, a finite-difference groundwater-flow model, was used to simulate flow in the Minnelusa and Madison hydrogeologic units with five layers. Layer 1 represented the fractured sandstone layers in the upper 250 ft of the Minnelusa hydrogeologic unit, and layer 2 represented the lower part of the Minnelusa hydrogeologic unit. Layer 3 represented the upper 150 ft of the Madison hydrogeologic unit, and layer 4 represented the less permeable lower part. Layer 5 represented an approximation of the underlying Deadwood aquifer to simulate upward flow to the Madison hydrogeologic unit. The finite-difference grid, oriented 23 degrees counterclockwise, included 221 rows and 169 columns with a square cell size of 492.1 ft in the detailed study area that surrounded Rapid City. The northern and southern boundaries for layers 1-4 were represented as no-flow boundaries, and the boundary on the east was represented with head-dependent flow cells. Streamflow recharge was represented with specified-flow cells, and areal recharge to layers 1-4 was represented with a specified-flux boundary. Calibration of the model was accomplished by two simulations: (1) steady-state simulation of average conditions for water years 1988-97 and (2) transient simulations of water years 1988-97 divided into twenty 6-month stress periods. Flow-system components represented in the model include recharge, discharge, and hydraulic properties. The steady-state streamflow recharge rate was 42.2 cubic feet per second (ft3/s), and transient streamflow recharge rates ranged from 14.1 to 102.2 ft3/s. The steady-state areal recharge rate was 20.9 ft3/s, and transient areal recharge rates ranged from 1.1 to 98.4 ft3/s. The upward flow rate from the Deadwood aquifer to the Madison hydrogeologic unit was 6.3 ft3/s. Discharge included springflow, water use, flow to overlying units, and regional outflow. The estimated steady-state springflow of 32.8 ft3/s from seven springs was similar to the simulated springflow of 31.6 ft3/s, which included 20.5 ft3

South Dakota↗

USGS2019

Estimated groundwater recharge from a water-budget model incorporating selected climate projections, Island of Maui, Hawai‘i

Demand for freshwater on the Island of Maui is expected to increase by 45 percent between 2015 and 2035. Groundwater availability on Maui is affected by changes in climate and agricultural irrigation. To evaluate the availability of fresh groundwater under projected future climate conditions and changing agricultural irrigation practices, estimates of groundwater recharge are needed. A water-budget model with a daily computation interval was used to estimate the spatial distribution of recharge on Maui for one present-day and two future-climate scenarios. All three scenarios used 2017 land cover. The two future-climate scenarios, including one wetter than the present-day scenario and one drier than the present-day scenario, were developed using available high-resolution downscaled climate projections. The drier future scenario was developed using projections for a Representative Concentration Pathway warming scenario during 2071–99 with total radiative forcing of 8.5 Watts per square meter by the year 2100 (RCP8.5 2071–99 scenario), whereas the wetter future scenario was developed using projections for a “Special Report on Emissions Scenarios” A1B emission scenario during 2080–99 (A1B 2080–99 scenario). For the RCP8.5 2071–99 scenario, projected mean annual recharge decrease for Maui is about 172 million gallons per day, or about 14 percent less than present-day recharge, which is estimated to be 1,232 million gallons per day. Recharge for the RCP8.5 2071–99 scenario is projected to decrease in 22 of Maui’s 25 aquifer systems, which are defined by the Hawaiʻi Commission on Water Resource Management. For the A1B 2080–99 future scenario, projected mean annual recharge increase for Maui is about 144 million gallons per day, or about 12 percent more than present-day recharge. Recharge for the A1B 2080–99 scenario is projected to increase in 17 of Maui’s 25 aquifer systems. Between the two future scenarios, a total of 11 aquifer systems show similar direction in drying (Kahului, Kama‘ole, Lualaʻilua, Makawao, Olowalu, Pāʻia, Ukumehame, Waikapū) or wetting (Honopou, Kawaipapa, and Waikamoi) changes for recharge. Selectively modifying the climate inputs for the A1B 2080–99 scenario indicates that the projected changes in rainfall account for most of the projected changes in recharge for Maui’s 25 aquifer systems. However, projected changes in reference evapotranspiration and forest-canopy evaporation also can account for a substantial part of the projected changes in recharge where changes in reference evapotranspiration are relatively large and where changes in forest-canopy evaporation extend across large forested areas. Projected changes in daily rainfall frequency have a relatively small but non-negligible impact on recharge estimates.

Hawaii↗

USGS2020

Analysis of remedial scenarios affecting plume movement through a sole-source aquifer system, southeastern Nassau County, New York

A steady-state three-dimensional groundwater-flow model based on present conditions is coupled with the particle-tracking program MODPATH to assess the fate and transport of volatile organic-compound plumes within the Magothy and upper glacial aquifers in southeastern Nassau County, New York. Particles are forward tracked from locations within plumes defined by surfaces of equal concentration. Particles move toward ultimate well capture and discharge to the general head and drain boundaries representing natural receptors in the models. Because rates of advection within coarse-grained sediments typically exceed 0.1 foot per day, mechanisms of dispersion and diffusion were assumed to be negligible. Resulting particle pathlines are influenced by hydrogeologic framework features and the interplay of nearby hydrologic stresses. Simulated hydrologic effects include cones of depression near pumping wells and water-table mounding near points of treated water recharge; however, remedial pumping amounts are balanced by treated-water return, and net effects at distant regional boundaries, including freshwater/saltwater interfaces, are minor. Once a steady-state model was developed and calibrated, eight hypothetical remedial scenarios were evaluated to hydraulically contain the volatile organic-compound plumes. Specifically, the remedial scenarios were optimized to achieve full containment by altering the pumping-well locations, adjusting the pumping rates, and adjusting the discharge locations and rates. Based on the results, total hypothetical extraction rates varied from about 5,462 gallons per minute during an anticipated near-future condition to about 13,340 gallons per minute during full hydraulic containment of all site-related compounds identified by the New York State standards, criteria, and guidance for environmental investigations and cleanup. Targeting of high-concentration zones of the plume increases the total amount of remedial pumpage necessary to capture all parts of the plume but may decrease the total amount of time necessary to operate a remedial system. Simulated time frames of advective transport ranged from about 12 years to capture zones with elevated concentrations of volatile organic compounds (mean particle travel time plus the standard deviation of travel time) to more than 100 years to capture all zones. Groundwater-flow model analysis indicates that all the optimal plume-containment scenarios would have negligible effects on streams and the saltwater-freshwater interface along the south shore of Long Island. Massapequa, Bellmore, Seaman, and Seaford Creeks are represented by using MODFLOW drain-boundary conditions. Saltwater-freshwater interfaces are represented by using MODFLOW general head-boundary conditions where the Magothy aquifer discharges upward into saline groundwater across the Gardiners clay confining unit and the Lloyd aquifer discharges upward into saline groundwater across the Raritan confining unit.

New York↗

USGS1987

Geology and hydrology of the Onondaga aquifer in eastern Erie County, New York, with emphasis on ground-water-level declines since 1982

The Onondaga aquifer is a nearly flat-lying, 25- to 110-foot-thick, cherty limestone with moderately developed karst features such as sinkholes, disappearing streams, and solution-widened joints. Most groundwater moves through solution-widened bedding planes, although some moves through vertical joints. The yield of water from 42 wells ranges from 3 to 100 gal/min, averaging 20 gal/min. Groundwater levels in the Onondaga aquifer declined during the fall of 1981 and summer and fall of 1982-85, near a 2.2-mile-long and 800-foot-wide land surface depression in the eastern part of Erie County. More than 60 wells and several wetlands went dry and at least three sinkholes developed. Groundwater levels were measured in 150 wells during a high water level period in April 1984 and a low water period in October 1984. Water levels fluctuated 20 to 50 ft near the depression and near the quarries but fluctuated only 5 to 10 ft elsewhere. The water level decline was caused by the combined effect of groundwater removal by pumpage from a quarry (the water is then discharged to Dorsch Creek) and by the swallets in the 2.2-mile-long depression area, which are recharge points for the aquifer. In 1982, sinkholes formed in a surface depression area in Harris Hill. The enlargement of these sinkholes seems to be unrelated to the water level decline in the eastern part of the county and is probably caused by local drainage alterations.

New York↗

USGS1987

AQMAN; linear and quadratic programming matrix generator using two-dimensional ground-water flow simulation for aquifer management modeling

A FORTRAN-77 computer program code that helps solve a variety of aquifer management problems involving the control of groundwater hydraulics. It is intended for use with any standard mathematical programming package that uses Mathematical Programming System input format. The computer program creates the input files to be used by the optimization program. These files contain all the hydrologic information and management objectives needed to solve the management problem. Used in conjunction with a mathematical programming code, the computer program identifies the pumping or recharge strategy that achieves a user 's management objective while maintaining groundwater hydraulic conditions within desired limits. The objective may be linear or quadratic, and may involve the minimization of pumping and recharge rates or of variable pumping costs. The problem may contain constraints on groundwater heads, gradients, and velocities for a complex, transient hydrologic system. Linear superposition of solutions to the transient, two-dimensional groundwater flow equation is used by the computer program in conjunction with the response matrix optimization method. A unit stress is applied at each decision well and transient responses at all control locations are computed using a modified version of the U.S. Geological Survey two dimensional aquifer simulation model. The program also computes discounted cost coefficients for the objective function and accounts for transient aquifer conditions. (Author 's abstract)

Water-Resources Investigations Report↗

USGS2012

Preliminary investigation of the effects of sea-level rise on groundwater levels in New Haven, Connecticut

Global sea level rose about 0.56 feet (ft) (170 millimeters (mm)) during the 20th century. Since the 1960s, sea level has risen at Bridgeport, Connecticut, about 0.38 ft (115 mm), at a rate of 0.008 ft (2.56 mm + or - 0.58 mm) per year. With regional subsidence, and with predicted global climate change, sea level is expected to continue to rise along the northeast coast of the United States through the 21st century. Increasing sea levels will cause groundwater levels in coastal areas to rise in order to adjust to the new conditions. Some regional climate models predict wetter climate in the northeastern United States under some scenarios. Scenarios for the resulting higher groundwater levels have the potential to inundate underground infrastructure in lowlying coastal cities. New Haven is a coastal city in Connecticut surrounded and bisected by tidally affected waters. Monitoring of water levels in wells in New Haven from August 2009 to July 2010 indicates the complex effects of urban influence on groundwater levels. The response of groundwater levels to recharge and season varied considerably from well to well. Groundwater temperatures varied seasonally, but were warmer than what was typical for Connecticut, and they seem to reflect the influence of the urban setting, including the effects of conduits for underground utilities. Specific conductance was elevated in many of the wells, indicating the influence of urban activities or seawater in Long Island Sound. A preliminary steady-state model of groundwater flow for part of New Haven was constructed using MODFLOW to simulate current groundwater levels (2009-2010) and future groundwater levels based on scenarios with a rise of 3 ft (0.91 meters (m)) in sea level, which is predicted for the end of the 21st century. An additional simulation was run assuming a 3-ft rise in sea level combined with a 12-percent increase in groundwater recharge. The model was constructed from existing hydrogeologic information for the New Haven area and from new information on groundwater levels collected during October 2009-June 2010. For the scenario with a 3-ft rise in sea level and no increase in recharge, simulated groundwater levels near the coast rose 3 ft; this increased water level tapered off toward a discharge area at the only nontidal stream in the study area. Simulated stream discharge increased at the nontidal stream because of the increased gradient. Although groundwater levels rose, the simulated difference between the groundwater levels in the aquifer and the increased sea level declined, indicating that the depth to the interface between freshwater and saltwater may possibly decline. Simulated water levels were affected by rise in sea level even in areas where the water table was at 17-24 ft (5.2-7.3 m) above current (2011) sea level. For the scenario with increased recharge, simulated groundwater levels were as much as an additional foot higher at some locations in the study area. The results of this preliminary investigation indicate that groundwater levels in coastal areas can be expected to rise and may rise higher if groundwater recharge also increases. This finding has implications for the disposal of stormwater through infiltration, a low-impact development practice designed to improve water quality and reduce overland peak discharge. Other implications include increased risk of basement flooding and increased groundwater seepage into underground sewer pipes and utility corridors in some areas. These implications will present engineering challenges to New Haven and Yale University. The preliminary model developed for this study can be the starting point for further simulation of future alternative scenarios for sea-level rise and recharge. Further simulations could identify those areas of New Haven where infrastructure may be at greatest risk from rising levels of groundwater. The simulations described in this report have limitations due to the preliminary scope of the work. Approaches to improve simulations include but are not limited to incorporating: * The variable density of seawater into the model in order to understand the current and future location of the interface between freshwater and saltwater; * Collection of additional data in order to better resolve temporal and spatial patterns in water levels in the aquifer; * Improved estimates of recharge through direct and indirect measurements of freshwater discharge from the study area; and * Transient simulations for greater understanding of the amount of time required for water levels and the position of the interface between freshwater and saltwater to adjust to changes in sea level and recharge.

Connecticut↗

USGS2013

Geohydrology, geochemistry, and groundwater simulation (1992-2011) and analysis of potential water-supply management options, 2010-60, of the Langford Basin, California

Groundwater withdrawals began in 1992 from the Langford Basin within the Fort Irwin National Training Center (NTC), California. From April 1992 to December 2010, approximately 12,300 acre-feet of water (averaging about 650 acre-feet per year) has been withdrawn from the basin and transported to the adjacent Irwin Basin. Since withdrawals began, water levels in the basin have declined by as much as 40 feet, and the quality of the groundwater withdrawn from the basin has deteriorated. The U.S. Geological Survey collected geohydrologic data from Langford Basin during 1992–2011 to determine the quantity and quality of groundwater available in the basin. Geophysical surveys, including gravity, seismic refraction, and time-domain electromagnetic induction surveys, were conducted to determine the depth and shape of the basin, to delineate depths to the Quaternary-Tertiary interface, and to map the depth to the water table and changes in water quality. Data were collected from existing wells and test holes, as well as 11 monitor wells that were installed at 5 sites as part of this study. Water-quality samples collected from wells in the basin were used to determine the groundwater chemistry within the basin and to delineate potential sources of poor-quality groundwater. Analysis of stable isotopes of oxygen and hydrogen in groundwater indicates that present-day precipitation is not a major source of recharge to the basin. Tritium and carbon-14 data indicate that most of the basin was recharged prior to 1952, and the groundwater in the basin has an apparent age of 12,500 to 30,000 years. Recharge to the basin, estimated to be less than 50 acre-feet per year, has not been sufficient to replenish the water that is being withdrawn from the basin. A numerical groundwater-flow model was developed for the Langford Basin to better understand the aquifer system used by the Fort Irwin NTC as part of its water supply, and to provide a tool to help manage groundwater resources at the NTC. Measured groundwater-level declines since the initiation of withdrawals (1992–2011) were used to calibrate the groundwater-flow model. The simulated recharge was about 46 acre-feet per year, including approximately 6 acre-feet per year of natural recharge derived from precipitation runoff and as much as 40 acre-feet per year of underflow from the Irwin Basin. Between April 1992 and December 2010, an average of about 650 acre-feet per year of water was withdrawn from the Langford Basin. Groundwater withdrawals in excess of natural recharge resulted in a net loss of 11,670 acre-feet of groundwater storage within the basin for the simulation period. The Fort Irwin NTC is considering various groundwater-management options to address the limited water resources in the Langford Basin. The calibrated Langford Basin groundwater-flow model was used to evaluate the hydrologic effects of four groundwater-withdrawal scenarios being considered by the Fort Irwin NTC over the next 50 years (January 2011 through December 2060). Continuation of the 2010 withdrawal rate in the three existing production wells will result in 70 feet of additional drawdown in the central part of the basin. Redistributing the 2010 withdrawal rate equally to the three existing wells and two proposed new wells in the northern and southern parts of the basin would result in about 10 feet less drawdown in the central part of the basin but about 100 feet of additional drawdown in the new well in the northern part of the basin and about 50 feet of additional drawdown in the new well in the southern part of the basin. Reducing the withdrawals from the three existing production wells in the central part of the basin from about 45,000 acre-feet to about 32,720 acre-feet would result in about 40 feet of additional drawdown in the central basin near the pumping wells, about 25 feet less than if withdrawals were not reduced. The combination of reducing and redistributing the cumulative withdrawals to the three existing and two proposed new wells results in about 40 feet of additional drawdown in the central and southern parts of the basin and about 70 feet in the northern part of the basin. These results show that reducing and redistributing the groundwater withdrawals would maintain the upper aquifer at greater than 50 percent of its predevelopment saturated thickness throughout the groundwater basin. The scenarios simulated for this study demonstrate how the calibrated model can be utilized to evaluate the hydrologic effects of different water-management strategies.

California↗

USGS2019

Drinking water health standards comparison and chemical analysis of groundwater for 72 domestic wells in Bradford County, Pennsylvania, 2016

Pennsylvania has the second highest number of residential wells of any state in the Nation with approximately 2.4 million residents that depend on groundwater for their domestic water supply. Despite the widespread reliance on groundwater in rural areas of the state, publicly available data to characterize the quality of private well water are limited. In Bradford County, more than half of the residents use groundwater from private domestic-supply wells as their primary drinking source. The quality of private well water is influenced by the regional and local setting, including the surrounding soil, geology, land use, household plumbing, and well construction. The groundwater used for domestic water supply in Bradford County is obtained primarily from shallow bedrock and from unconsolidated (glacial) deposits that overlie the bedrock. Historical land use has been predominately forested, agricultural, and residential, but more recently unconventional oil/gas development has been distributed throughout the landscape. Pennsylvania is one of only two states in the Nation without statewide water-well construction standards. To better assess the quality of groundwater used for drinking water supply in Bradford County, data for 72 domestic wells were collected and analyzed for a wide range of constituents that could be evaluated in relation to drinking water health standards, geology, land use, and other environmental factors. Groundwater samples were collected from May through August 2016 and analyzed for physical and chemical properties, including major ions, nutrients, trace elements, volatile organic compounds, ethylene and propylene glycol, alcohols, gross-alpha/beta-particle activity, uranium, radon-222, and dissolved gases. A subset of samples was analyzed for radium isotopes (radium-226 and -228) and for the isotopic composition of methane. This study was conducted by the U.S. Geological Survey in cooperation with the Northern Tier Regional Planning and Development Commission and is part of a regional effort to characterize groundwater in rural areas of Pennsylvania. Results of the 2016 study show that groundwater quality generally met most drinking-water standards. However, a percentage of samples failed to meet maximum contaminant levels (MCLs) for total coliform bacteria (49.3 percent), Escherichia coli (8.5 percent), barium (2.8 percent), and arsenic (2.8 percent); and secondary maximum contaminant levels (SMCL) for sodium (48.6 percent), manganese (30.6 percent), gross alpha and beta activity (16.7 percent), iron (11.1 percent), pH (8.3 percent), total dissolved solids (5.6 percent), chloride (1.4 percent), and aluminum (1.4 percent). Radon-222 activities exceeded the proposed drinking-water standard of 300 picocuries per liter (pCi/L) in 70.4 percent of the samples. There were no exceedances of drinking water health standards for any volatile organic compounds, and the only detections were for three trihalomethanes in one sample. The pH of the groundwater had a large influence on chemical characteristics and ranged from 6.18 to 9.31. Generally, the higher pH samples had higher potential for elevated concentrations of several constituents, including total dissolved solids, sodium, lithium, chloride, fluoride, boron, arsenic, and methane. For the Bradford County well-water samples, calcium/bicarbonate type waters were most abundant, with others classified as sodium/bicarbonate or mixed water types including calcium-sodium/bicarbonate, calcium-sodium/bicarbonate-chloride, sodium/bicarbonate-chloride, sodium/bicarbonate-sulfate, or sodium/chloride types. Six principal components (pH, redox, hardness, chloride-bromide, strontium-barium, and molybdenum-arsenic) explained nearly 78.3 percent of the variance in the groundwater dataset. Groundwater from 12.5 percent of the wells had concentrations of methane greater than the Pennsylvania action level of 7 milligrams per liter (mg/L); detectable methane concentrations ranged from 0.01 to 77 mg/L. In addition, low levels of ethane (as much as 0.13 mg/L) were present in seven samples with the highest methane concentrations. The isotopic composition of methane in five of these groundwater samples was consistent with the isotopic compositions reported for mud-gas logging samples from these geologic units and a thermogenic source. Isotopic composition from a sixth sample suggested the methane in that sample may be of microbial origin. Well-water samples with the higher methane concentrations also had higher pH values and elevated concentrations of sodium, lithium, boron, fluoride, arsenic, and bromide. Relatively elevated concentrations of some other constituents, such as barium and chloride, commonly were present in, but not limited to, those well-water samples with elevated methane. Four of the six groundwater samples with the highest methane concentrations had chloride/bromide ratios that indicate mixing with a small amount of brine (0.02 percent or less) similar in composition to those reported for gas and oil well brines in Pennsylvania. In several other eastern Pennsylvania counties where gas drilling is absent, groundwater with comparable chloride/bromide ratios and chloride concentrations have been reported, implying a potential natural source of brine. Most of Bradford County well-water samples have chloride concentrations less than 20 mg/L, and those with higher chloride concentrations have chloride/bromide ratios that indicate anthropogenic sources (such as road-deicing salt and septic effluent) or brine. Brines that are naturally present may originate from deeper parts of the aquifer system, whereas anthropogenic sources are more likely to affect shallow groundwater because they occur on or near the land surface. The available data for this study indicate that no one physical factor, such as the topographic setting, well depth, or altitude at the bottom of the well, was particularly useful for predicting those well locations with an elevated dissolved concentration of methane. The 2016 assessment of groundwater quality in Bradford County shows groundwater is generally of good quality, but methane and some constituents that occur in high concentration in naturally occurring brine and also in produced waters may be present at low to moderate concentrations in groundwater in various parts of the aquifer.

Pennsylvania↗

USGS2010

Geohydrology of the stratified-drift aquifer system in the lower Sixmile Creek and Willseyville Creek trough, Tompkins County, New York

In 2002, the U.S. Geological Survey, in cooperation with the Tompkins County Planning Department began a series of studies of the stratified-drift aquifers in Tompkins County to provide geohydrologic data for planners to develop a strategy to manage and protect their water resources. This aquifer study in lower Sixmile Creek and Willseyville Creek trough is the second in a series of aquifer studies in Tompkins County. The study area is within the northern area of the Appalachian Plateau and extends about 9 miles from the boundary between Tompkins County and Tioga County in the south to just south of the City of Ithaca in the north. In lower Sixmile Creek and Willseyville Creek trough, confined sand and gravel aquifers comprise the major water-bearing units while less extensive unconfined units form minor aquifers. About 600 people who live in lower Sixmile Creek and Willseyville Creek trough rely on groundwater from the stratified-drift aquifer system. In addition, water is used by non-permanent residents such as staff at commercial facilities. The estimated total groundwater withdrawn for domestic use is about 45,000 gallons per day (gal/d) or 0.07 cubic foot per second (ft3/s) based on an average water use of 75 gal/d per person for self-supplied water systems in New York. Scouring of bedrock in the preglacial lower Sixmile Creek and Willseyville Creek valleys by glaciers and subglacial meltwaters truncated hillside spurs, formed U-shaped, transverse valley profiles, smoothed valley walls, and deepened the valleys by as much as 300 feet (ft), forming a continuous trough. The unconsolidated deposits in the study area consist mostly of glacial drift, both unstratified drift (till) and stratified drift (laminated lake, deltaic, and glaciofluvial sediments), as well as some post-glacial stratified sediments (lake-bottom sediments that were deposited in reservoirs, peat and muck that were deposited in wetlands, and alluvium deposited by streams). Multiple advances and retreats of the ice in the study area resulted in several sequences of various types of glacial deposits. A large moraine (Valley Heads Moraine) dominates the southern part of the study area, a large delta dominates the central part, and ground moraine (mostly till) dominates the northern part. Glacial sediments in the center of the lower Sixmile Creek and Willseyville Creek trough typically range from 150 to 200 ft but can be greater than 300 ft in some places. Where the sediments are composed of sand and gravel they form aquifers. In most parts of the lower Sixmile Creek and Willseyville Creek trough, there is an upper and a basal confined aquifer. However, underlying the central parts of the Brooktondale delta, there are as many as four confined aquifers, whereas in the northern part of the study area, only one extensive confined aquifer is present. The major sources of recharge to these confined aquifers are (1) direct infiltration of precipitation where confined aquifers crop out at land surface (mostly along the western trough wall in the southern and central parts of the study area and, to a lesser degree, along the eastern trough wall); (2) unchanneled surface and subsurface runoff from adjacent upland areas that seeps into the aquifer along the western trough walls; (3) subsurface flow from underlying till or bedrock at the lateral contacts at trough walls; (4) adjacent fine-grained stratified drift, especially when the aquifer is pumped; and (5) discharge from bedrock at the bottom and sides of the trough. In the central part of the study area, the surficial coarse-grained sediments (sand and gravel) comprise a delta near Brooktondale and form a small unconfined aquifer (0.3 square mile). Although much of the upper part of the delta has been removed by several aggregate mining operations, sufficient amounts of sand and gravel remain in most places to form a thin unconfined aquifer. The major sources of recharge to the unconfined aquifer are (1)

Scientific Investigations Report↗

USGS2020

Machine learning predictions of pH in the Glacial Aquifer System, Northern USA

A boosted regression tree model was developed to predict pH conditions in three dimensions throughout the glacial aquifer system of the contiguous United States using pH measurements in samples from 18,386 wells and predictor variables that represent aspects of the hydrogeologic setting. Model results indicate that the carbonate content of soils and aquifer materials strongly controls pH and, when coupled with long flowpaths, results in the most alkaline conditions. Conversely, in areas where glacial sediments are thin and carbonate‐poor, pH conditions remain acidic. At depths typical of drinking‐water supplies, predicted pH >7.5—which is associated with arsenic mobilization—occurs more frequently than predicted pH <6—which is associated with water corrosivity and the mobilization of other trace elements. A novel aspect of this model was the inclusion of numerically based estimates of groundwater flow characteristics (age and flowpath length) as predictor variables. The sensitivity of pH predictions to these variables was consistent with hydrologic understanding of groundwater flow systems and the geochemical evolution of groundwater quality. The model was not developed to provide precise estimates of pH at any given location. Rather, it can be used to more generally identify areas where contaminants may be mobilized into groundwater and where corrosivity issues may be of concern to prioritize areas for future groundwater monitoring.

Groundwater↗

USGS2025

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater↗

USGS2015

The quality of our Nation's waters: Water quality in the Upper Floridan aquifer and overlying surficial aquifers, southeastern United States, 1993-2010

About 10 million people rely on groundwater from the Upper Floridan and surficial aquifers for drinking water. The Upper Floridan aquifer also is of primary importance to the region as a source of water for irrigation and as a source of crystal clear water that discharges to springs and streams providing recreational and tourist destinations and unique aquatic habitats. The reliance of the region on the Upper Floridan aquifer for drinking water and for the tourism and agricultural economies highlights the importance of long-term management to sustain the availability and quality of these resources.

Florida↗

USGS2017

Effects of changes in pumping on regional groundwater-flow paths, 2005 and 2010, and areas contributing recharge to discharging wells, 1990–2010, in the vicinity of North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

A previously developed regional groundwater flow model was used to simulate the effects of changes in pumping rates on groundwater-flow paths and extent of recharge discharging to wells for a contaminated fractured bedrock aquifer in southeastern Pennsylvania. Groundwater in the vicinity of the North Penn Area 7 Superfund site, Montgomery County, Pennsylvania, was found to be contaminated with organic compounds, such as trichloroethylene (TCE), in 1979. At the time contamination was discovered, groundwater from the underlying fractured bedrock (shale) aquifer was the main source of supply for public drinking water and industrial use. As part of technical support to the U.S. Environmental Protection Agency (EPA) during the Remedial Investigation of the North Penn Area 7 Superfund site from 2000 to 2005, the U.S. Geological Survey (USGS) developed a model of regional groundwater flow to describe changes in groundwater flow and contaminant directions as a result of changes in pumping. Subsequently, large decreases in TCE concentrations (as much as 400 micrograms per liter) were measured in groundwater samples collected by the EPA from selected wells in 2010 compared to 2005‒06 concentrations. To provide insight on the fate of potentially contaminated groundwater during the period of generally decreasing pumping rates from 1990 to 2010, steady-state simulations were run using the previously developed groundwater-flow model for two conditions prior to extensive remediation, 1990 and 2000, two conditions subsequent to some remediation 2005 and 2010, and a No Pumping case, representing pre-development or cessation of pumping conditions. The model was used to (1) quantify the amount of recharge, including potentially contaminated recharge from sources near the land surface, that discharged to wells or streams and (2) delineate the areas contributing recharge that discharged to wells or streams for the five conditions. In all simulations, groundwater divides differed from surface-water divides, partly because of differences in stream elevations and because of geologic structure and pumping. In the 1990 and 2000 simulations, all recharge in and near the vicinity of North Penn Area 7 discharged to wells, but in the 2005 and 2010 simulations some recharge in this area discharged to streams, indicating possible discharge of contaminated groundwater from North Penn Area 7 sources to streams. As the amount of groundwater withdrawals by wells has declined since 1990, the area contributing recharge to wells in the vicinity of North Penn Area 7 has decreased. To determine the effect of changes in pumping on flow paths and possible flow-path-related contributions to the observed changes in spatial distribution of contaminants in groundwater from 2005 to 2010, the USGS conducted simulations using the previously developed regional groundwater-flow model using reported pumping and estimated recharge rates for 2005 and 2010. Flow paths from recharge at known contaminant source areas to discharge locations at wells or streams were simulated under steady-state conditions for the two periods. Simulated groundwater-flow paths shifted only slightly from 2005 to 2010 as a result of changes in pumping rates. These slight changes in groundwater-flow paths from known sources of contamination are not coincident with the spatial distribution of observed changes in TCE concentrations from 2005 to 2010, indicating that the decreases of TCE concentrations may be a result of other processes, such as source removal or degradation. Results of the simulations and the absence of increases in TCE-degradation-product concentrations indicate that the decreases of TCE concentrations observed in 2010 may be at least partly related to contaminant-source removal by soil excavation completed in 2005, although additional data would be needed to confirm this preliminary explanation.

Pennsylvania↗

USGS1987

Geohydrology and susceptibility of major aquifers to surface contamination in Alabama; Area 8

The U.S. Geological Survey, in cooperation with the Alabama Department of Environmental Management, is conducting a series of geohydrologic studies to delineate the major aquifers and their susceptibility to contamination in Alabama. This report delineates and describes the geohydrology and susceptibility of the major aquifers to contamination in Area 8--Autauga, Chilton, Elmore, Lowndes, and Montgomery Counties. The major aquifers in the study area are the Eutaw, Gordo, and Coker aquifers of Cretaceous age. One or more of these aquifers are sources of public water supply in each of the five counties. The recharge areas for these aquifers are in Autauga, Chilton, Elmore, and Montgomery and Prattville. Maximum groundwater use in the Prattville area is more than 8 mgd (million gallons per day). Estimated maximum groundwater withdrawal for all uses in the study area is about 65 mgd. The potentiometric map of the Gordo aquifer indicates that the Alabama River may serve as a recharging boundary to the Gordo aquifer along the flood plain of the river in the Montgomery-Prattville area. The river also is acting as a recharging boundary to the Eutaw and Coker aquifers, where the potentiometric surfaces in the aquifers have been lowered. All recharge areas for the major aquifers are susceptible to contamination from the surface. However, the areas that are highly susceptible to contamination extend from Jemison to Clanton in Chilton County where the Coker aquifer generally is < 100 ft below land surface, and the flood plains of the Alabama, Coosa, and Tallapoosa Rivers, which are underlain by alluvial deposits that are in hydraulic contact with the major aquifers. Within the highly susceptible areas, the areas especially susceptible to contamination are the flood plain of the Alabama River in the Montgomery area and the flood plain of the Tallapoosa River. Pumpage from the major aquifers in this area has significantly lowered the potentiometric surface in the aquifers resulting in a downward gradient between the major aquifers and the Alabama River and the alluvial deposits underlying the flood plain along the river. (Lantz-PTT)

Water-Resources Investigations Report↗

USGS2016

Potentiometric surface and water-level difference maps of selected confined aquifers in Southern Maryland and Maryland’s Eastern Shore, 1975-2015

Key Results This report presents potentiometric-surface maps of the Aquia and Magothy aquifers and the Upper Patapsco, Lower Patapsco, and Patuxent aquifer systems using water levels measured during September 2015. Water-level difference maps are also presented for these aquifers. The water-level differences in the Aquia aquifer are shown using groundwater-level data from 1982 and 2015, while the water-level differences are shown for the Magothy aquifer using data from 1975 and 2015. Water-level difference maps for both the Upper Patapsco and Lower Patapsco aquifer systems are shown using data from 1990 and 2015. The water-level differences in the Patuxent aquifer system are shown using groundwater-level data from 2007 and 2015. The potentiometric surface maps show water levels ranging from 53 feet above sea level to 164 feet below sea level in the Aquia aquifer, from 86 feet above sea level to 106 feet below sea level in the Magothy aquifer, from 115 feet above sea level to 115 feet below sea level in the Upper Patapsco aquifer system, from 106 feet above sea level to 194 feet below sea level in the Lower Patapsco aquifer system, and from 165 feet above sea level to 171 feet below sea level in the Patuxent aquifer system. Water levels have declined by as much as 116 feet in the Aquia aquifer since 1982, 99 feet in the Magothy aquifer since 1975, 66 and 83 feet in the Upper Patapsco and Lower Patapsco aquifer systems, respectively, since 1990, and 80 feet in the Patuxent aquifer system since 2007.

Report↗

USGS2019

Groundwater flow model for Western Chippewa County–Including analysis of water resources related to industrial sand mining and irrigated agriculture

A groundwater flow model for western Chippewa County, Wisconsin, was developed by the Wisconsin Geological and Natural History Survey (WGNHS) and the U.S. Geological Survey (USGS) using the computer program MODFLOW. The model is the result of a five-year groundwater study commissioned by Chippewa County in 2012 to evaluate the effects of industrial sand mining and irrigated agriculture on the county’s water resources. The study incorporates existing data and newly acquired data from fieldwork conducted within the study area. The groundwater model may be useful for future investigations, such as evaluation of proposed high-capacity well sites, development of municipal wellhead protection plans, and studies that seek to further quantify surface water-groundwater relationships. The model conceptualizes the hydrostratigraphy of western Chippewa County as six stacked layers. Each layer is distinct, beginning with unlithified glacial material at the surface, and alternating between sandstones (that act as aquifers) and shale units (that serve as aquitards). The model is bounded below by Precambrian crystalline bedrock and its perimeter was derived from a regional-scale groundwater flow model. The MODFLOW model represented average conditions during 2011–2013 with “steady-state” assumptions, meaning that simulated water levels do not fluctuate seasonally or from year to year. Steady-state models simplify natural variability, making results of scenario simulations easier to interpret and compare while also maximizing effects of stressors because the simulated stress is always applied (not halted after a few months or years). Model calibration used the parameter estimation code (PEST), and calibration targets included heads (groundwater levels) and streamflows. Calibration focused on 2011–2013 because a large amount of head and streamflow data were available for that period. The MODFLOW model explicitly simulates all sources and sinks of water, including groundwater/surface-water interaction with streamflow routing. Model input included estimates of aquifer hydraulic conductivity and a spatial groundwater recharge distribution developed using a GIS-based soil-water-balance (SWB) model applied to the model area. Groundwater withdrawals were simulated for 269 high-capacity wells across the entire model domain, which includes western Chippewa County and adjacent portions of Dunn, Barron, and Rusk Counties. Collectively, these wells withdrew about 1.14 million gallons per year between 2011 and 2013. Once the model was calibrated, it was applied to two distinct scenarios of increased groundwater withdrawals: one evaluating hydrologic effects of more intensive industrial sand mining and the second evaluating the hydrologic effects of more intensive agricultural irrigation practices. Each scenario was developed with input from Chippewa County and a stakeholder group established expressly for this study. The scenarios were designed to represent reasonable future buildout conditions for both mining and irrigated agriculture. The mining scenario underscores the potential hydrologic effects related to changing land-use practices (i.e., hilltops and farmland becoming sand mines), while the irrigated agriculture scenario illustrates the potential hydrologic effects of intensifying existing land-use practices (i.e., installing new wells to irrigate farm fields). While each scenario evaluated distinctly different conditions, modeling results demonstrated the potential of both scenarios to lower the water table and reduce baseflows in headwater streams within the modeled area. In the case of irrigated agriculture, hydrologic effects were associated directly with groundwater withdrawals. By assuming that irrigation did not decrease, this steady-state simulation represented a sustained future effect. By contrast, hydrologic effects of industrial sand mining were the result of both groundwater withdrawals at mines and land-use changes that effectively reduced recharge to groundwater over distinct phases of active mining. This scenario included a post-mining phase, during which groundwater withdrawals stopped and mined areas were reclaimed to undeveloped prairie grass cover. If reclamation to undeveloped prairie indeed occurs as simulated, long-term increases in the water table and stream baseflows are possible. In this sense, the scenario representing build out of irrigated agriculture led to long-term baseflow declines while the future buildout of industrial sand mining led to declines that dissipated following mine reclamation to undisturbed prairie. Future investigations in similar hydrogeologic settings may find the following insights gleaned from this study useful: ❚❚ The characterization of hydrogeologic properties, delineation of hydrogeologic units, and calibration of groundwater flow models benefited from incorporation of accurate well construction reports, high-quality borehole geophysical logs, and streamflow gaging data. ❚❚ Infiltration testing performed in active mining areas provided evidence that reducing the degree and extent of compaction and enhancing areas designed to retain and infiltrate stormwater runoff could potentially reduce runoff and increase groundwater recharge. ❚❚ Similarly, reclaiming mined areas to prairie grasses would be expected to reduce runoff and increase groundwater recharge by reducing compaction and improving soil structure and vegetation that can slow runoff and enhance infiltration.

Wisconsin↗

USGS1988

Solute geochemistry of the Snake River plain regional aquifer system, Idaho and eastern Oregon

Four geochemical approaches were used to determine chemical reactions controlling solute concentrations in the Snake River Plain regional aquifer system: (1) calculation of a solute balance within the aquifer, (2) identification of weathered products in the aquifer frame- work, (3) comparison of thermodynamic mineral saturation indices with plausible solute reactions, and (4) comparison of stable-isotope ratios of the solutes with those in the aquifer framework. Solutes in the geo- thermal groundwater system underlying the main aquifer were examined by calculating thermodynamic mineral saturation indices, stable-isotope ratios, geothermometry, and radiocarbon dating. Water budgets, hydrologic arguments, and isotopic analyses for the eastern Snake River Plain aquifer system demonstrate that most, if not all, water is of local meteoric and not juvenile or formation origin. Thus, the solutes must also originate within the basin. Solute balance, isotopic, mineralogic, and thermodynamic arguments suggest that about 20 per- cent of the solutes leaving the basin are derived from reactions with rocks forming the aquifer framework. Most of the remaining solutes are introduced from tributary drainage basins. Mass-balance calculations, thermodynamic arguments, and petro- graphic observations indicate that calcite and silica are precipitated in the aquifer. Petrographic evidence and thermodynamic arguments sug- gest that olivine, pyroxene, plagioclase, pyrite, and anhydrite are being weathered from the aquifer framework. Large amounts of sodium, chloride, and sulfate, relative to their concentration in the igneous rock, are being removed from the aquifer. Release of fluids from inclusions in the igneous rocks and initial flushing of grain boundaries and pores of detrital marine sediments in interbeds are believed to be a major source of these solutes. Identification and quantification of reactions controlling solute concentrations in ground water in the eastern plain indicate that the aquifer is not a large mixing vessel that simply stores and transmits water and solutes but is undergoing diagenesis and is both a source and a sink for solutes. Evaluation of solute concentrations and stable-isotope ratios of hydrogen, oxygen, carbon, and sulfur along groundwater flowpaths that transect irrigated areas suggests that irrigation water may have altered solute concentrations and isotope ratios in the eastern Snake River Plain aquifer system. The changes, however, have been small because of the similarity of solute concentrations and ratios in applied irrigation water and in native ground water, and because of rapid movement and large dispersivity of the aquifer. Reactions controlling solutes in the western Snake River basin are believed to be similar to those in the eastern basin but, because of dif- ferent hydrologic conditions, a definitive analysis could not be made. The regional geothermal system that underlies the Snake River Plain contains total dissolved solids similar to those in the overlying Snake River Plain aquifer system but contains higher concentrations of sodium, bicarbonate, silica, fluoride, sulfate, chloride, arsenic, boron, and lithium, and lower concentrations of calcium, magnesium, and hydrogen. These solutes are believed to be derived from reactions similar to those in the Snake River Plain aquifer system, except that ion exchange and hydrol- ysis play a role in controlling solute concentrations in the geothermal system. Geothermometry calculations of selected ground-water samples from known geothermal areas throughout the basin suggest that the geother- mal system is large in areal extent but has relatively low temperatures. Approximately half of the silica-quartz calculated water temperatures are greater than 90 °C. Radiocarbon dating of geothermal water in the Salmon Falls and Bruneau-Grand View areas in the south central part of the Snake River basin suggests that residence time of the geother- mal water is about 17,700 years.

Idaho;Oregon↗

USGS2010

Groundwater-quality data and regional trends in the Virginia Coastal Plain, 1906-2007

A newly developed regional perspective of the hydrogeology of the Virginia Coastal Plain incorporates updated information on groundwater quality in the area. Local-scale groundwater-quality information is provided by a comprehensive dataset compiled from multiple Federal and State agency databases. Groundwater-sample chemical-constituent values and related data are presented in tables, summaries, location maps, and discussions of data quality and limitations. Spatial trends in groundwater quality and related processes at the regional scale are determined from interpretive analyses of the sample data. Major ions that dominate the chemical composition of groundwater in the deep Piney Point, Aquia, and Potomac aquifers evolve eastward and with depth from (1) 'hard' water, dominated by calcium and magnesium cations and bicarbonate and carbonate anions, to (2) 'soft' water, dominated by sodium and potassium cations and bicarbonate and carbonate anions, and lastly to (3) 'salty' water, dominated by sodium and potassium cations and chloride anions. Chemical weathering of subsurface sediments is followed by ion exchange by clay and glauconite, and subsequently by mixing with seawater along the saltwater-transition zone. The chemical composition of groundwater in the shallower surficial and Yorktown-Eastover aquifers, and in basement bedrock along the Fall Zone, is more variable as a result of short flow paths between closely located recharge and discharge areas and possibly some solutes originating from human sources. The saltwater-transition zone is generally broad and landward-dipping, based on groundwater chloride concentrations that increase eastward and with depth. The configuration is convoluted across the Chesapeake Bay impact crater, however, where it is warped and mounded along zones having vertically inverted chloride concentrations that decrease with depth. Fresh groundwater has flushed seawater from subsurface sediments preferentially around the impact crater as a result of broad contrasts between sediment permeabilities. Paths of differential flushing are also focused along the inverted zones, which follow stratigraphic and structural trends southeastward into North Carolina and northeastward beneath the chloride mound across the outer impact crater. Brine within the inner impact crater has probably remained unflushed. Regional movement of the saltwater-transition zone takes place over geologic time scales. Localized movement has been induced by groundwater withdrawal, mostly along shallow parts of the saltwater-transition zone. Short-term episodic withdrawals result in repeated cycles of upconing and downconing of saltwater, which are superimposed on longer-term lateral saltwater intrusion. Effective monitoring for saltwater intrusion needs to address multiple and complexly distributed areas of potential intrusion that vary over time. A broad belt of large groundwater fluoride concentrations underlies the city of Suffolk, and thins and tapers northward. Fluoride in groundwater probably originates by desorbtion from phosphatic sedimentary material. The high fluoride belt possibly was formed by initial adsorbtion of fluoride onto sediment oxyhydroxides, followed by desorbtion along the leading edge of the advancing saltwater-transition zone. Large groundwater iron and manganese concentrations are most common to the west along the Fall Zone, across part of the saltwater-transition zone and eastward, and within shallow groundwater far to the east. Iron and manganese initially produced by mineral dissolution along the Fall Zone are adsorbed eastward and with depth by clay and glauconite, and subsequently desorbed along the leading edge of the advancing saltwater-transition zone. Iron and manganese in shallow groundwater far to the east are produced by reaction of sediment organic matter with oxyhydroxides. Large groundwater nitrate and ammonium concentrations are mostly limited to shallow depths. Most nitrate a

Professional Paper↗