Search USGS⌕ 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 109 records · Page 6Linked to original sources

Simulation of groundwater flow in the regional aquifer system on Long Island, New York, for pumping and recharge conditions in 2005–15

A three-dimensional groundwater-flow model was developed for the aquifer system of Long Island, New York, to evaluate (1) responses of the hydrologic system to changes in natural and anthropogenic hydraulic stresses, (2) the subsurface distribution of groundwater age, and (3) the regional-scale distribution of groundwater travel times and the source of water to fresh surface waters and coastal receiving waters. The model also provides the groundwater flow components used to define model boundaries for possible inset models used for local-scale analyses. The three-dimensional, groundwater flow model developed for this investigation uses the numerical code MODFLOW–NWT to represent steady-state conditions for average groundwater pumping and aquifer recharge for 2005–15. The particle-tracking algorithm MODPATH, which simulates advective transport in the aquifer, was used to estimate groundwater age, delineate the areas at the water table that contribute recharge to coastal and freshwater bodies, and estimate total travel times of water from the water table to discharge locations. A three-dimensional, 1-meter (3.3-foot) topobathymetric model was used to determine land-surface altitudes for the island and seabed altitudes for the surrounding coastal waters. The mapped extents and surface altitudes of major geologic units were compiled and used to develop a three-dimensional hydrogeologic framework of the aquifer system, including aquifers and confining units. Lithologic data from deep boreholes and previous aquifer-test results were used to estimate the three-dimensional distribution of hydraulic conductivity in principal aquifers. Natural recharge from precipitation was estimated for 2005–15 using a modified Thornthwaite-Mather methodology as implemented in a soil-water balance model. Components of anthropogenic recharge—wastewater return flow, storm water inflow, and inflow from leaky infrastructure—also were estimated for 2005–15. Groundwater withdrawals for various sources, including public water supply, industrial, remediation, and agricultural, were compiled or estimated for the same period. These data were incorporated into the model development to represent the aquifer system geometry, boundaries, and initial hydraulic properties of the regional aquifers and confining units within the Long Island aquifer system. Average hydraulic conditions—water levels and streamflows—for 2005–15 were estimated using existing data from the U.S. Geological Survey National Water Information System database. Model inputs were adjusted to best match average hydrologic conditions using inverse methods as implemented in the parameter-estimating software PEST. The calibrated model was used to simulate average hydrologic conditions in the aquifer system for 2005–15. About 656 cubic feet per second of water was withdrawn on average annually for 2005–15 for water supply and an average of about 349 cubic feet per second of water recharged the aquifer annually from return flow and leaky infrastructure. Parts of New York City have drawdowns exceeding 25 feet, mostly because of urbanization and associated large decreases in recharge rates. Large areas in the western part of the island have drawdowns exceeding 10 feet, mostly from large groundwater withdrawals and sewering, which largely eliminates wastewater return flow. Water-table altitudes in eastern parts of the island increased by more than 2 feet in some areas as a result of wastewater return flow in unsewered areas and changes in land use. Changes in streamflows show a similar pattern as water-table altitudes. Streamflows generally decrease in western parts of the island where there are large drawdowns and increase in eastern parts of the island where water-table altitudes increase.

New York↗

Assessment of ethylene dibromide, dibromochloropropane, other volatile organic compounds, radium isotopes, radon, and inorganic compounds in groundwater and spring water from the Crouch Branch and McQueen Branch aquifers near McBee, South Carolina, 2010-2012

Public-supply wells near the rural town of McBee, in southwestern Chesterfield County, South Carolina, have provided potable water to more than 35,000 residents throughout Chesterfield County since the early 1990s. Groundwater samples collected between 2002 and 2008 in the McBee area by South Carolina Department of Health and Environmental Control (DHEC) officials indicated that groundwater from two public-supply wells was characterized by the anthropogenic compounds ethylene dibromide (EDB) and dibromochloropropane (DBCP) at concentrations that exceeded their respective maximum contaminant levels (MCLs) established by the U.S. Environmental Protection Agency’s (EPA) National Primary Drinking Water Regulations (NPDWR). Groundwater samples from all public-supply wells in the McBee area were characterized by the naturally occurring isotopes of radium-226 and radium-228 at concentrations that approached, and in one well exceeded, the MCL for the combined isotopes. The local water utility installed granulated activated carbon filtration units at the two EDB- and DBCP-contaminated wells and has, since 2011, shut down these two wells. Groundwater pumped by the remaining public-supply wells is currently (2014) centrally treated at a water-filtration plant. To assess the occurrence, distribution, and potential sources of the anthropogenic and naturally occurring compounds detected in groundwater in the McBee area, samples of groundwater and spring water were collected from public-supply, domestic-supply, agricultural-supply, and monitoring wells and springs, respectively, between 2010 and 2012 by the U.S. Geological Survey. The water samples were analyzed for concentrations of EDB, DBCP, other volatile organic compounds (VOCs), radium-226 and radium-228, radon, and inorganic compounds. All wells sampled were screened in the shallow Crouch Branch aquifer, the deeper McQueen Branch aquifer, or, for most public-supply wells, both aquifers. In areas where no wells existed or wells could not be installed, passive samplers that adsorb EDB, DBCP, and various VOCs, were installed in the shallow subsurface. A representative groundwater flow pathway to each public supply well and selected other wells was determined by using a calibrated three-dimensional groundwater-flow model of the Atlantic Coastal Plain in Chesterfield County and particle-tracking analysis. The aerial extent of the groundwater flow pathway to public-supply wells was mapped by using chlorofluorocarbon-concentration based, apparent-age dates of the groundwater. The water-quality data collected between 2010 and 2012, in conjunction with groundwater flow pathways and historical aerial photographs of land uses near McBee, indicate an area where EDB-, DBCP-, 1,2-dichloropropane-, 1,3-dichloropropane-, and carbon disulfide-contaminated groundwater exists in the Crouch Branch aquifer in the Cedar Creek Basin and north of McBee and is most likely related to the past use of these compounds between the early 1900s and the 1980s as soil fumigants in predominately agricultural areas north of McBee. The highest EDB concentration detected (18.6 micrograms per liter) during the 3-year study was in a groundwater sample from an agricultural-supply well located north of McBee. Other VOCs, such as dichloromethane and 1,1,2-trichloroethane, also were detected in groundwater samples from this EDB-contaminated agricultural-supply well but are from unknown source(s). The fact that the agricultural area north of McBee is located in a recharge area for the Crouch Branch aquifer most likely facilitated the groundwater contamination in this area. DBCP-contaminated groundwater detected in three public-supply wells south of McBee in the deeper McQueen Branch aquifer appears to be related to past soil fumigation practices that used DBCP in agricultural areas located south of McBee. One of the three DBCP-contaminated public-supply wells also contained EDB, most likely present in groundwater due to the release of leaded gasolines that contained EDB as a fuel additive between the 1940s and 1970s. A gasoline-source of EDB, rather than a soil-fumigation source, is supported by the co-detection in groundwater from the well of 1,2-dichloroethane, a lead scavenger compound also added to leaded gasoline. Groundwater pumped from two public-supply wells located within and to the east of the McBee town limits and one domestic-supply well east of McBee was characterized by the detection of 1,1-dichloroethane, trichloroethylene, 1,1-dichloroethylene, and perchloroethylene. Groundwater flow pathways determined for these wells indicate that the potential source(s) of these compounds detected in one public-supply well and the domestic-supply well may be located within the McBee town limits, and that the potential source(s) of these compounds detected in the public-supply well to the east of McBee may be located in an area north of McBee formerly used for agriculture, but used for industry since at least the 1970s. Radium isotopes (defined in this study as the sum of radium-226 and radium-228 concentrations) and radon were detected in all wells sampled in the McBee area between 2010 and 2012. Wells characterized by radium isotope concentrations in groundwater that exceeded the MCL of 5.0 picocuries per liter were also characterized by specific conductance values greater than 30 microsiemens per centimeter and clustered north of McBee in a predominately agricultural area, and in agricultural and urban areas located within and east of McBee. The elevated specific conductance values measured in groundwater from these wells most likely are due to recharge by water mineralized by fertilizer application in agricultural areas, or due to the recharge by water mineralized by septic-tank drain-field effluent near urban areas. Radon was detected in groundwater from all wells sampled, and radon concentrations in groundwater from three monitoring wells exceeded the proposed MCL of 300 picocuries per liter. Concentrations of uranium in groundwater in the McBee area increased with increased groundwater-sample depth, most likely due to the proximity of the sample-collection location to basement rock that contains uranium-bearing minerals.

South Carolina↗

A multiphased approach to groundwater investigations for the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer is a vital groundwater resource for agricultural, industrial, and public supply uses in the Pecos County region of western Texas. Resource managers would like to understand the future availability of water in the Edwards-Trinity aquifer in the Pecos County region and the effects of the possible increase or temporal redistribution of groundwater withdrawals. To provide resource managers with that information, the U.S. Geological Survey (USGS), in cooperation with the Middle Pecos Groundwater Conservation District, Pecos County, City of Fort Stockton, Brewster County, and Pecos County Water Control and Improvement District No. 1, completed a three-phase study of the Edwards-Trinity and related aquifers in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. The first phase was to collect groundwater, surface-water, geochemical, geophysical, and geologic data in the study area and develop a geodatabase of historical and collected data. Data compiled in the first phase of the study were used to develop the conceptual model in the second phase of the study. The third phase of the study involved the development and calibration of a numerical groundwater-flow model of the Edwards-Trinity aquifer to simulate groundwater conditions based on various groundwater-withdrawal scenarios. Analysis of well, geophysical, geochemical, and hydrologic data contributed to the development of the conceptual model in phase 1. Lithologic information obtained from well reports and geophysical data was used to describe the hydrostratigraphy and structural features of the groundwater-flow system, and aquifer-test data were used to estimate aquifer hydraulic properties. Geochemical data were used to evaluate groundwater-flow paths, water-rock interaction, aquifer interaction, and the mixing of water from different sources in phase 2. Groundwater-level data also were used to evaluate aquifer interaction, as well as to develop a potentiometric-surface map, delineate regional groundwater divides, and describe regional groundwater-flow paths. During phase 3, the data collected and compiled along with the conceptual information in the study area were incorporated into a numerical groundwater-flow model to evaluate the sustainability of recent (2008) and projected water-use demands on groundwater resources in the study area.

Texas↗

An update of hydrologic conditions and distribution of selected constituents in water, Snake River Plain aquifer and perched groundwater zones, Idaho National Laboratory, Idaho, emphasis 2006-08

Since 1952, radiochemical and chemical wastewater discharged to infiltration ponds (also called percolation ponds), evaporation ponds, and disposal wells at the Idaho National Laboratory (INL) has affected water quality in the eastern Snake River Plain aquifer and perched groundwater zones underlying the INL. The U.S. Geological Survey, in cooperation with the U.S. Department of Energy, maintains groundwater monitoring networks at the INL to determine hydrologic trends, and to delineate the movement of radiochemical and chemical wastes in the aquifer and in perched groundwater zones. This report presents an analysis of water-level and water-quality data collected from aquifer and perched groundwater wells in the USGS groundwater monitoring networks during 2006-08. Water in the Snake River Plain aquifer primarily moves through fractures and interflow zones in basalt, generally flows southwestward, and eventually discharges at springs along the Snake River. The aquifer primarily is recharged from infiltration of irrigation water, infiltration of streamflow, groundwater inflow from adjoining mountain drainage basins, and infiltration of precipitation. From March-May 2005 to March-May 2008, water levels in wells generally remained constant or rose slightly in the southwestern corner of the INL. Water levels declined in the central and northern parts of the INL. The declines ranged from about 1 to 3 feet in the central part of the INL, to as much as 9 feet in the northern part of the INL. Water levels in perched groundwater wells around the Advanced Test Reactor Complex (ATRC) also declined. Detectable concentrations of radiochemical constituents in water samples from wells in the Snake River Plain aquifer at the INL generally decreased or remained constant during 2006-08. Decreases in concentrations were attributed to decreased rates of radioactive-waste disposal, radioactive decay, changes in waste-disposal methods, and dilution from recharge and underflow. In April or October 2008, reportable concentrations of tritium in groundwater ranged from 810 ? 70 to 8,570 ? 190 picocuries per liter (pCi/L), and the tritium plume extended south-southwestward in the general direction of groundwater flow. Tritium concentrations in water from wells completed in shallow perched groundwater at the ATRC were less than the reporting levels. Tritium concentrations in deep perched groundwater exceeded the reporting level in 11 wells during at least one sampling event during 2006-08 at the ATRC. Tritium concentrations from one or more zones in each well were reportable in water samples collected at various depths in six wells equipped with multi-level WestbayTM packer sampling systems. Concentrations of strontium-90 in water from 24 of 52 aquifer wells sampled during April or October 2008 exceeded the reporting level. Concentrations ranged from 2.2 ? 0.7 to 32.7 ? 1.2 pCi/L. Strontium-90 has not been detected within the eastern Snake River Plain aquifer beneath the ATRC partly because of the exclusive use of waste-disposal ponds and lined evaporation ponds rather than using the disposal well for radioactive-wastewater disposal at ATRC. At the ATRC, the strontium-90 concentration in water from one well completed in shallow perched groundwater was less than the reporting level. During at least one sampling event during 2006-08, concentrations of strontium-90 in water from nine wells completed in deep perched groundwater at the ATRC were greater than reporting levels. Concentrations ranged from 2.1?0.7 to 70.5?1.8 pCi/L. At the Idaho Nuclear Technology and Engineering Center (INTEC), the reporting level was exceeded in water from two wells completed in deep perched groundwater. During 2006-08, concentrations of cesium-137, plutonium-238, and plutonium-239, -240 (undivided), and americium-241 were less than the reporting level in water samples from all wells and all zones in wells equipped with multi-level WestbayTM packer sampling systems

Scientific Investigations Report↗

Groundwater-Level Elevations in the Denver Basin Bedrock Aquifers of Elbert County, Colorado, 2015–18

Public and domestic water supplies in Elbert County, Colorado, rely on groundwater withdrawals from five bedrock aquifers in the Denver Basin aquifer system (lower Dawson, upper Dawson, Denver, Arapahoe, and Laramie-Fox Hills) to meet water demands. Increased pumping in response to regional population growth and development has led to declining groundwater levels in neighboring Douglas County. The U.S. Geological Survey, in cooperation with the Elbert County Board of County Commissioners, began a study in 2015 to monitor groundwater levels within Elbert County. The purpose of this study is to report on groundwater levels measured between April 2015 and June 2018, and analyze trends and changes in groundwater-level elevations throughout the county. Discrete groundwater levels were measured at 42 wells within Elbert County. Six of those wells contained equipment to make and record continuous groundwater-level measurements at hourly intervals. All five aquifers had wells with a rise in groundwater-level elevation and wells with a decline in groundwater-level elevation, based on a relative change in groundwater-level elevation between the April 2015 and April 2018 measurements. All aquifers except the upper Dawson had more wells with significant negative trends in discrete groundwater-level elevations than significant positive trends; however, at least one well within the upper Dawson, lower Dawson, Arapahoe, and Laramie-Fox Hills aquifers had a significant positive trend. Wells screened in the lower Dawson aquifer consistently had the most significant negative trends, with an average trend of −1.96 feet per year (ft/year). The upper Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers had average trends of 0.03 ft/year, −1.04 ft/year, −0.46 ft/year, and −0.65 ft/year, respectively. Trends in continuous groundwater-level elevations were in agreement with significant trends in discrete groundwater-level elevations. Potentiometric-surface maps of the upper and lower Dawson aquifers for April 2015 and April 2018 show that differences in hydraulic head from the two measurement periods were greatest along the western part of Elbert County. Results of this study could guide future groundwater monitoring in the county and aid in long-term planning of water resources.

Colorado↗

Groundwater chemistry and water-level elevations in bedrock aquifers of the Piceance and Yellow Creek watersheds, Rio Blanco County, Colorado, 2013–16

The Piceance and Yellow Creek watersheds in Rio Blanco County, Colorado, are known to contain important energy resources (oil shale and natural gas) and mineral resources (nahcolite). The primary sources of fresh groundwater in the Piceance and Yellow Creek watersheds are bedrock aquifers in the Uinta and Green River Formations. The aquifers are divided into an upper and lower aquifer separated by a regionally extensive semiconfining layer. These aquifers provide water to streams and springs in the watersheds and are an important resource to people living and working in the area. Development of energy and mineral resources has the potential to affect the quality of groundwater in several ways. The Bureau of Land Management and the U.S. Geological Survey began groundwater monitoring in 2010 to characterize the groundwater quality and water-level elevations of shallow bedrock aquifers in the Piceance and Yellow Creek watersheds. The purpose of this report is to present ground-water chemistry and water-level elevations collected during 2013–16. Comparisons are made to data that were collected from the bedrock aquifers from 2010 to 2012 to identify the potential for changes in water quality and water-level elevations. Appreciable changes in water-level elevations and hydraulic gradient were observed in early April 2015 in two wells completed in the upper and lower aquifers. The hydraulic gradient between the two wells was consistently downward from the upper aquifer to the lower aquifer during 2010–15; however, in early April 2015, the gradient changed from downward to upward between the two aquifers. Overall, water-level elevations declined by about 14 and 11 feet in the upper and lower aquifers, respectively, from 2013 to 2016. Previously published data estimated groundwater ages at 1,200 years old in the upper aquifer and 9,600 years old in the lower aquifer. These groundwater ages indicate that ground-water was recharged over thousands of years. With such long periods of time for aquifer recharge, declines in water-level elevation over short time steps (a few months) have important implications for sustainable management of this resource. Solution mining activities or drilling for oil and natural gas in the area could be related to the changes observed in water-level elevations in these wells; however, further investigation would be needed to evaluate causation. Changes in major-ion chemistry were evaluated in the bedrock aquifer using time series plots of select major-ion data from 2010 to 2016. Major-ion chemistry was variable for a single well from 2010 to 2016 where alkalinity and sulfate were the most variable constituents. One possible explanation for the observed changes in major-ion chemistry may be that the sample depth for that well no longer represents the most appreciable flow in the borehole. On a larger scale, potential changes in flow within the borehole may indicate changes in the regional flow system. Methane and volatile organic compound concentrations were evaluated using a similar approach to that of major ions and had similar findings. Methane concentrations in wells sampled from 2010 to 2016 were generally constant. The only exception was observed at a single well where the range of methane concentrations was from 57.4 (2010) to 4.02 milligrams per liter (2013). This is the same well where changes in water-level elevation, hydraulic gradient, and major-ion chemistry were observed, providing multiple lines of evidence to indicate change in the bedrock aquifers. Sampling of a well located in an area with little energy development but where faults or fractures could provide a path for the migration of fluids indicate mixing of groundwater between the upper and lower aquifers.

Colorado↗

Evidence for terrigenic SF6 in groundwater from basaltic aquifers, Jeju Island, Korea: Implications for groundwater dating

Measurements of the concentrations of dichlorodifluoromethane (CFC-12), tritium ( 3 H), and sulfur hexafluoride (SF 6 ) in groundwater from basaltic aquifers in Jeju Island, Korea, demonstrate a terrigenic source of SF 6 . Using a lumped-parameter dispersion model, groundwater was identified as young water (<15 years), old water with negligible CFC-12 and 3 H, and binary mixtures of the two. Model calculations using dispersion models and binary mixing based on 3 H and CFC-12 concentrations demonstrate a non-atmospheric excess of SF 6 relative to CFC-12 and 3 H concentrations for more than half of the samples. The non-atmospheric excess SF 6 may have originated from terrigenic sources in relict volcanic fluids, which could have acquired SF 6 from granites and basement rocks of the island during volcanic activity. Local excess anthropogenic sources of SF 6 are unlikely. The SF 6 age is biased young relative to the CFC-12 age, typically up to 20 years and as high as 30 years. This age bias is more pronounced in samples of groundwater older than 15 years. The presence of terrigenic SF 6 can affect the entire dating range for groundwater in mixtures that contain a fraction of old water.

Jeju Island↗

Groundwater-quality data in the Monterey–Salinas shallow aquifer study unit, 2013: Results from the California GAMA Program

Groundwater quality in the 3,016-square-mile Monterey–Salinas Shallow Aquifer study unit was investigated by the U.S. Geological Survey (USGS) from October 2012 to May 2013 as part of the California State Water Resources Control Board Groundwater Ambient Monitoring and Assessment (GAMA) Program’s Priority Basin Project. The GAMA Monterey–Salinas Shallow Aquifer study was designed to provide a spatially unbiased assessment of untreated-groundwater quality in the shallow-aquifer systems in parts of Monterey and San Luis Obispo Counties and to facilitate statistically consistent comparisons of untreated-groundwater quality throughout California. The shallow-aquifer system in the Monterey–Salinas Shallow Aquifer study unit was defined as those parts of the aquifer system shallower than the perforated depth intervals of public-supply wells, which generally corresponds to the part of the aquifer system used by domestic wells. Groundwater quality in the shallow aquifers can differ from the quality in the deeper water-bearing zones; shallow groundwater can be more vulnerable to surficial contamination. Samples were collected from 170 sites that were selected by using a spatially distributed, randomized grid-based method. The study unit was divided into 4 study areas, each study area was divided into grid cells, and 1 well was sampled in each of the 100 grid cells (grid wells). The grid wells were domestic wells or wells with screen depths similar to those in nearby domestic wells. A greater spatial density of data was achieved in 2 of the study areas by dividing grid cells in those study areas into subcells, and in 70 subcells, samples were collected from exterior faucets at sites where there were domestic wells or wells with screen depths similar to those in nearby domestic wells (shallow-well tap sites). Field water-quality indicators (dissolved oxygen, water temperature, pH, and specific conductance) were measured, and samples for analysis of inorganic constituents (trace elements, nutrients, major and minor ions, silica, total dissolved solids, and alkalinity) were collected at all 170 sites. In addition to these constituents, the samples from grid wells were analyzed for organic constituents (volatile organic compounds, pesticides and pesticide degradates), constituents of special interest (perchlorate and N -nitrosodimethylamine, or NDMA), radioactive constituents (radon-222 and gross-alpha and gross-beta radioactivity), and geochemical and age-dating tracers (stable isotopes of carbon in dissolved inorganic carbon, carbon-14 abundances, stable isotopes of hydrogen and oxygen in water, and tritium activities). Three types of quality-control samples (blanks, replicates, and matrix spikes) were collected at up to 11 percent of the wells in the Monterey–Salinas Shallow Aquifer study unit, and the results for these samples were used to evaluate the quality of the data from the groundwater samples. With the exception of trace elements, blanks rarely contained detectable concentrations of any constituent, indicating that contamination from sample-collection procedures was not a significant source of bias in the data for the groundwater samples. Low concentrations of some trace elements were detected in blanks; therefore, the data were re-censored at higher reporting levels. Replicate samples generally were within the limits of acceptable analytical reproducibility. The median values of matrix-spike recoveries were within the acceptable range (70 to 130 percent) for the volatile organic compounds (VOCs) and N -nitrosodimethylamine (NDMA), but were only approximately 64 percent for pesticides and pesticide degradates. The sample-collection protocols used in this study were designed to obtain representative samples of groundwater. The quality of groundwater can differ from the quality of drinking water because water chemistry can change as a result of contact with plumbing systems or the atmosphere; because of treatment, disinfection, or blending with water from other sources; or some combination of these. Water quality in domestic wells is not regulated in California, however, to provide context for the water-quality data presented in this report, results were compared to benchmarks established for drinking-water quality. The primary comparison benchmarks were maximum contaminant levels established by the U.S. Environmental Protection Agency and the State of California (MCL-US and MCL-CA, respectively). Non-regulatory benchmarks were used for constituents without maximum contaminant levels (MCLs), including Health Based Screening Levels (HBSLs) developed by the USGS and State of California secondary maximum contaminant levels (SMCL-CA) and notification levels. Most constituents detected in samples from the Monterey–Salinas Shallow Aquifer study unit had concentrations less than their respective benchmarks. Of the 148 organic constituents analyzed in the 100 grid-well samples, 38 were detected, and all concentrations were less than the benchmarks. Volatile organic compounds were detected in 26 of the grid wells, and pesticides and pesticide degradates were detected in 28 grid wells. The special-interest constituent NDMA was detected above the HBSL in three samples, one of which also had a perchlorate concentration greater than the MCL-CA. Of the inorganic constituents, 6 were detected at concentrations above their respective MCL benchmarks in grid-well samples: arsenic (5 grid wells above the MCL of 10 micrograms per liter, μg/L), selenium (3 grid wells, MCL of 50 μg/L), uranium (4 grid wells, MCL of 30 μg/L), nitrate (16 grid wells, MCL of 10 milligrams per liter, mg/L), adjusted gross alpha particle activity (10 grid wells, MCL of 15 picocuries per liter, pCi/L), and gross beta particle activity (1 grid well, MCL of 50 pCi/L). An additional 4 inorganic constituents were detected at concentrations above their respective HBSL benchmarks in grid-well samples: boron (1 grid well above the HBSL of 6,000 μg/L), manganese (8 grid wells, HBSL of 300 μg/L), molybdenum (6 grid wells, HBSL of 40 μg/L), and strontium (6 grid wells, HBSL of 4,000 μg/L). Of the inorganic constituents, 4 were detected at concentrations above their non-health based SMCL benchmarks in grid-well samples: iron (9 grid wells above the SMCL of 300 μg/L), chloride (7 grid wells, SMCL of 500 mg/L), sulfate (14 grid wells, SMCL of 500 mg/L), and total dissolved solids (27 grid wells, SMCL of 1,000 mg/L). Of the inorganic constituents analyzed in the 70 shallow-well tap sites, 10 were detected at concentrations above the benchmarks. Of the inorganic constituents, 3 were detected at concentrations above their respective MCL benchmarks in shallow-well tap sites: arsenic (2 shallow-well tap sites above the MCL of 10 μg/L), uranium (2 shallow-well tap sites, MCL of 30 μg/L), and nitrate (24 shallow-well tap sites, MCL of 10 mg/L). An additional 3 inorganic constituents were detected above their respective HBSL benchmarks in shallow-well tap sites: manganese (4 shallow-well tap sites above the HBSL of 300 μg/L), molybdenum (4 shallow-well tap sites, HBSL of 40 μg/L), and zinc (2 shallow-well tap sites, HBSL of 2,000 μg/L). Of the inorganic constituents, 4 were detected at concentrations above their non-health based SMCL benchmarks in shallow-well tap sites: iron (6 shallow-well tap sites above the SMCL of 300 μg/L), chloride (1 shallow-well tap site, SMCL of 500 mg/L), sulfate (9 shallow-well tap sites, SMCL of 500 mg/L), and total dissolved solids (15 shallow-well tap sites, SMCL of 1,000 mg/L).

California↗

Potentiometric surfaces (2013, 2015), groundwater quality (2010–15), and water-level changes (2011–13, 2013–15) in the Sparta-Memphis aquifer in Arkansas

The Sparta-Memphis aquifer, present across much of eastern Arkansas, is the second most used groundwater resource in the State, with the Mississippi River Valley alluvial aquifer being the primary groundwater resource. The U.S. Geological Survey, in cooperation with Arkansas Department of Agriculture-Natural Resources Division, Arkansas Geological Survey, Natural Resources Conservation Service, Union County Water Conservation Board, and the Union County Conservation District, collects groundwater data across the Sparta-Memphis aquifer extent in Arkansas. This report presents water-level data for measurements conducted during two time periods, January–May 2013 and January–June 2015, and discusses water-level altitude changes for the 2011–13 and 2013–15 periods in the Sparta-Memphis aquifer. Accompanying water-level data in this report include groundwater-quality data for the period 2010–15 in the Sparta-Memphis aquifer. Groundwater data can guide ongoing and future groundwater-monitoring efforts and inform management of the aquifers in Arkansas. Water levels measured at 306 wells from January to May 2013 and 273 wells from January to June 2015 are graphically presented as potentiometric-surface maps. Measurements from 2011, 2013, and 2015 were used in the construction of 2011–13 and 2013–15 water-level change maps. Select long-term hydrographs are included in the report to illustrate water-level changes at the local scale. Water-level data show the influence of climate, pumping, and conservation and management efforts on groundwater levels. With respect to climate, the study area experienced extreme drought conditions between January 2011 and December 2012. The proximate effects of drought—increased evapotranspiration, decreased recharge, and increased irrigation needs—resulted in water-level declines that were particularly notable in the northern and central portions of the study area. Groundwater sampled in 2010–15 from 148 wells completed in the Sparta-Memphis aquifer was analyzed for specific conductance, pH, chloride (Cl) concentration, and bromide (Br) concentration. In 2015, groundwater-quality data from 103 wells completed in the Sparta-Memphis aquifer had a median specific conductance of 356 microsiemens per centimeter at 25 degrees Celsius and a median Cl concentration of 9.5 milligrams per liter (mg/L). The data show two areas of higher Cl (greater than 10 mg/L) and higher Br (greater than 0.5 mg/L) concentrations in Union, Calhoun, and Bradley Counties in southern Arkansas and Monroe and Phillips Counties in eastern-central Arkansas. A Cl and Br mixing model indicates the two regions of wells may have different sources of higher salinity. In the greater Union County area, water in most wells may be a mixture of recharge or precipitation and higher salinity groundwater from the Nacatoch aquifer. Water in wells in eastern-central Arkansas may be sourced from aquifers having a higher Cl concentration (and thus, also a higher Cl-to-Br ratio).

Arkansas↗

Status and understanding of groundwater quality in the North San Francisco Bay Shallow Aquifer study unit, 2012; California GAMA Priority Basin Project (ver. 1.1, February 2018)

Groundwater quality in the North San Francisco Bay Shallow Aquifer study unit (NSF-SA) was investigated as part of the Priority Basin Project of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is in Marin, Mendocino, Napa, Solano, and Sonoma Counties and included two physiographic study areas: the Valleys and Plains area and the surrounding Highlands area. The NSF-SA focused on groundwater resources used for domestic drinking water supply, which generally correspond to shallower parts of aquifer systems than that of groundwater resources used for public drinking water supply in the same area. The assessments characterized the quality of untreated groundwater, not the quality of drinking water. This study included three components: (1) a status assessment , which characterized the status of the quality of the groundwater resources used for domestic supply for 2012; (2) an understanding assessment , which evaluated the natural and human factors potentially affecting water quality in those resources; and (3) a comparison between the groundwater resources used for domestic supply and those used for public supply. The status assessment was based on data collected from 71 sites sampled by the U.S. Geological Survey for the GAMA Priority Basin Project in 2012. To provide context, concentrations of constituents measured in groundwater were compared to U.S. Environmental Protection Agency (EPA) and California State Water Resources Control Board Division of Drinking Water regulatory and non-regulatory benchmarks for drinking-water quality. The status assessment used a grid-based method to estimate the proportion of the groundwater resources that has concentrations of water-quality constituents approaching or above benchmark concentrations. This method provides statistically unbiased results at the study-area scale and permits comparisons to other GAMA Priority Basin Project study areas. In the NSF-SA study unit as a whole, inorganic constituents with human-health benchmarks were detected at high relative concentrations (RCs) in 27 percent of the shallow aquifer system, and inorganic constituents with secondary maximum contaminant levels (SMCL) were detected at high RCs in 24 percent of the system. The inorganic constituents detected at high RCs were arsenic, boron, fluoride, manganese, nitrate, iron, sulfate, and total dissolved solids (TDS). Organic constituents with human-health benchmarks were detected at high RCs in 1 percent of the shallow aquifer system. Of the 148 organic constituents analyzed, 30 constituents were detected, although only 1, chloroform, had a detection frequency greater than 10 percent. Natural and anthropogenic factors that could affect the groundwater quality were evaluated by using results from statistical testing of associations between constituent concentrations and values of potential explanatory factors. Groundwater age class (modern, mixed, or pre-modern), redox class (oxic or anoxic), aquifer lithology class (metamorphic, sedimentary, or volcanic), and dissolved oxygen concentrations were the explanatory factors that explained distribution patterns of most of the inorganic constituents best. Groundwater classified primarily as pre-modern or mixed in age was associated with higher concentrations of arsenic and fluoride than waters classified as modern. Anoxic or mixed redox conditions were associated with higher concentrations of boron, fluoride, and manganese. Similar patterns of association with explanatory variables were seen for inorganic constituents with aesthetic-based benchmarks detected at high concentrations. Nitrate and perchlorate had higher concentrations in oxic than in the anoxic redox class and were positively correlated with urban land use. The NSF-SA water-quality results were compared to those of the GAMA North San Francisco Bay Public-Supply Aquifer study unit (NSF-PA). The NSF-PA was sampled in 2004 and covers much of the same area as the NSF-SA, but focused on the deeper public-supply aquifer system. The comparison of the NSF-PA to the NSF-SA showed that there were more differences between the Valleys and Plains study areas of the two study units than between the Highlands study areas of the two study units. As expected from the shallower depth of wells, the NSF-SA Valleys and Plains study area had a lesser proportion of pre-modern age groundwater and greater proportion of modern age groundwater than the NSF-PA Valleys and Plains study area. In contrast, well depths and groundwater ages were not significantly different between the two Highlands study areas. Arsenic, manganese, and nitrate were present at high RCs, and perchlorate was detected in greater proportions of the NSF-SA Valleys and Plains study area than the NSF-PA Valleys and Plains study area.

California↗

Groundwater-flow model for the Wood River Valley aquifer system, south-central Idaho

A three-dimensional numerical model of groundwater flow was developed for the Wood River Valley (WRV) aquifer system, Idaho, to evaluate groundwater and surface-water availability at the regional scale. This mountain valley is located in Blaine County and has a drainage area of about 2,300 square kilometers (888 square miles). The model described in this report can serve as a tool for water-rights administration and water-resource management and planning. The model was completed with support from the Idaho Department of Water Resources, and is part of an ongoing U.S. Geological Survey effort to characterize the groundwater resources of the WRV. A highly reproducible approach was taken for constructing the WRV groundwater-flow model. The collection of datasets, source code, and processing instructions used to construct and analyze the model was distributed as an R statistical-computing and graphics package. Flow in the WRV aquifer was simulated using the MODFLOW-USG groundwater flow model. The transient flow model simulates groundwater flow between 1995 and 2010. The model uses a 100-meter (328-feet) uniform grid spacing with 54,922 active model cells distributed over three model layers. A confining unit in the south-central part of the Bellevue fan necessitated the use of a multi-layer model. Specified-flow boundaries were used to simulate the groundwater inflows from each of the major tributary basins (also known as tributary basin underflow) and the areal recharge of precipitation and applied irrigation. Head‑dependent flow boundaries were used to simulate the stream-aquifer flow exchange in river reaches and the groundwater discharge at the outlet boundaries of Stanton Crossing and Silver Creek. The model was calibrated by adjusting aquifer hydraulic properties to match simulated and measured water levels and stream-aquifer flow exchange, using the parameter-estimation program PEST. The model reasonably simulated the measured water-table elevation, orientation, and gradients. Stream-aquifer flow exchange along river reaches also was reasonably simulated by the model. Inflow into the WRV aquifer system originates from three sources (from largest to smallest): Streamflow loss from the Big Wood River and Silver Creek, Areal recharge of precipitation and applied irrigation, and Tributary basin underflow. Outflow from the WRV aquifer system originates from five sources (from largest to smallest): Aquifer discharge into the Big Wood River and Silver Creek, Production-well pumping, Discharge across the water table into the vadose zone, Subsurface outflow beneath Silver Creek near Picabo, and Subsurface outflow beneath the Big Wood River near Stanton Crossing. Temporal changes in aquifer storage are most affected by areal recharge and groundwater pumping, and also contribute to changes in streamflow gains.

Idaho↗

Hydrogeologic framework, arsenic distribution, and groundwater geochemistry of the glacial-sediment aquifer at the Auburn Road landfill superfund site, Londonderry, New Hampshire

Leachate continues to be generated from landfills at the Auburn Road Landfill Superfund Site in Londonderry, New Hampshire. Impermeable caps on the three landfills at the site inhibit direct infiltration of precipitation; however, high water-table conditions allow groundwater to interact with landfill materials from below, creating leachate and ultimately reducing conditions in downgradient groundwater. Reducing conditions can facilitate arsenic transport by allowing it to stay in solution or by liberating arsenic adsorbed to surfaces and from geologic sources, such as glacial sediments and bedrock. The site occupies a 180-acre parcel of land containing streams, ponds, wetlands, and former gravel pits located in glacial sediment. Four areas, totaling 14 acres, including three landfills and one septage lagoon, were used for waste disposal. The site was closed in 1980 after volatile organic compounds associated with industrial waste dumping were detected. The site was added to the U.S. Environmental Protection Agency National Priority List in 1982, and the landfills were capped in 1996. Although volatile organic compound concentrations in groundwater have declined substantially, some measurable concentrations remain. Temporally variable and persistent elevated arsenic concentrations have been measured in groundwater affected by the landfill leachate. Microbial consumption of carbon found in leachate is a driver of reducing conditions that liberate arsenic at the site. In addition to sources of carbon in landfill leachate, wetland areas throughout the site also could contribute carbon to groundwater, but it is currently unknown if any of the wetland areas have downward or reversing gradients that could allow the infiltration of surface water to groundwater. Red-stained sediments and water indicate iron-rich groundwater discharge to surface water and are also associated with elevated concentrations of arsenic in sediment and groundwater. Ironrich groundwater seeps have been observed in the wetland, streams, and pond downgradient of the landfills. Piezometers were installed in some of these locations to confirm groundwater discharge, measure vertical-flow gradients, and to provide a way to sample the discharging groundwater. Understanding the movement of leachate in groundwater is complicated by the presence of preferential flow paths through aquifer materials with differing hydraulic properties; these preferential flow paths can affect rates of recharge, geochemical conditions, and contaminant fluxes. In areas adjacent to the three capped landfills, infiltration of precipitation containing oxygenated water through permeable deltaic sediments in the former gravel pit area causes increases in dissolved oxygen concentrations and decreases in arsenic concentrations. Layered deltaic sediments produce anisotropic hydraulic characteristics and zones of high hydraulic conductivity. The glacial-sediment aquifer also includes glaciolacustrine sediments that have low permeability and limit infiltration at the surface Discharge of leachate-affected groundwater may be limited in areas of organic muck on the bottom of Whispering Pines Pond because the muck may act as a semiconfining layer. Geophysical survey results were used to identify several areas with continuous beds of muck and an underlying highresistivity layer on top of a layer of low resistivity that may represent leachate-affected groundwater. The high-resistivity layer is likely groundwater associated with oxygenated recharge, which would cause arsenic to adsorb onto aquifer sediments and reduce concentrations of dissolved arsenic in groundwater. Surface and borehole geophysical data collected in 2011 were used to identify potentially high-permeability or contaminated zones in the aquifer (preferential flowpaths) as well as low-permeability zones that may promote contamination through back diffusion. Some groundwater in parts of the glacial-sediment aquifer where the leachate plumes were present had low electrical resistivity, low dissolved oxygen, and high concentrations of arsenic. Low-resistivity zones in the underlying bedrock were associated with fractures that also may contain leachate. Although surveying the fractured bedrock was not a specific objective of this study, the results suggest that such a survey would help to determine if leachate and associated concentrations of arsenic are migrating downward into the fractured-bedrock-aquifer system. An uncalibrated, one-dimensional, reactive-transport model was used to assess several conditions that affect arsenic mobility. The results indicate that reductive dissolution and desorption from glacial sediments control dissolved arsenic concentrations. Parameter sensitivity analysis was used to identify key data that are needed in order to accurately assess the time required for arsenic concentrations to fall to levels below the maximum contaminant level at the site. Quantifying this time will require accurate characterization of carbon, sediment-surface sorption sites, and groundwater fluxes at the site.

New Hampshire↗

Preliminary simulation of chloride transport in the Equus Beds aquifer and simulated effects of well pumping and artificial recharge on groundwater flow and chloride transport near the city of Wichita, Kansas, 1990 through 2008

The Equus Beds aquifer in south-central Kansas is a primary water-supply source for the city of Wichita. Water-level declines because of groundwater pumping for municipal and irrigation needs as well as sporadic drought conditions have caused concern about the adequacy of the Equus Beds aquifer as a future water supply for Wichita. In March 2006, the city of Wichita began construction of the Equus Beds Aquifer Storage and Recovery project, a plan to artificially recharge the aquifer with excess water from the Little Arkansas River. Artificial recharge will raise groundwater levels, increase storage volume in the aquifer, and deter or slow down a plume of chloride brine approaching the Wichita well field from the Burrton, Kansas area caused by oil production activities in the 1930s. Another source of high chloride water to the aquifer is the Arkansas River. This study was prepared in cooperation with the city of Wichita as part of the Equus Beds Aquifer Storage and Recovery project. Chloride transport in the Equus Beds aquifer was simulated between the Arkansas and Little Arkansas Rivers near the Wichita well field. Chloride transport was simulated for the Equus Beds aquifer using SEAWAT, a computer program that combines the groundwater-flow model MODFLOW-2000 and the solute-transport model MT3DMS. The chloride-transport model was used to simulate the period from 1990 through 2008 and the effects of five well pumping scenarios and one artificial recharge scenario. The chloride distribution in the aquifer for the beginning of 1990 was interpolated from groundwater samples from around that time, and the chloride concentrations in rivers for the study period were interpolated from surface water samples. Five well-pumping scenarios and one artificial-recharge scenario were assessed for their effects on simulated chloride transport and water levels in and around the Wichita well field. The scenarios were: (1) existing 1990 through 2008 pumping conditions, to serve as a baseline scenario for comparison with the hypothetical scenarios; (2) no pumping in the model area, to demonstrate the chloride movement without the influence of well pumping; (3) double municipal pumping from the Wichita well field with existing irrigation pumping; (4) existing municipal pumping with no irrigation pumping in the model area; (5) double municipal pumping in the Wichita well field and no irrigation pumping in the model area; and (6) increasing artificial recharge to the Phase 1 Artificial Storage and Recovery project sites by 2,300 acre-feet per year. The effects of the hypothetical pumping and artificial recharge scenarios on simulated chloride transport were measured by comparing the rate of movement of the 250-milligrams-per-liter-chloride front for each hypothetical scenario with the baseline scenario at the Arkansas River area near the southern part of the Wichita well field and the Burrton plume area. The scenarios that increased the rate of movement the most compared to the baseline scenario of existing pumping between the Arkansas River and the southern boundary of the well field were those that doubled the city of Wichita’s pumping from the well field (scenarios 3 and 5), increasing the rate of movement by 50 to 150 feet per year, with the highest rate increases in the shallow layer and the lowest rate increases in the deepest layer. The no pumping and no irrigation pumping scenarios (2 and 4) slowed the rate of movement in this area by 150 to 210 feet per year and 40 to 70 feet per year, respectively. In the double Wichita pumping scenario (3), the rate of movement in the shallow layer of the Burrton area decreased by about 50 feet per year. Simulated chloride rate of movement in the deeper layers of the Burrton area was decreased in the no pumping and no irrigation scenarios (2 and 4) by 80 to 120 feet per year and 50 feet per year, respectively, and increased in the scenarios that double Wichita’s pumping (3 and 5) from the well field by zero to 130 feet per year, with the largest increases in the deepest layer. In the increased Phase 1 artificial recharge scenario (6), the rate of chloride movement in the Burrton area increased in the shallow layer by about 30 feet per year, and decreased in the middle and deepest layer by about 10 and 60 feet per year, respectively. Comparisons of the rate of movement of the simulated 250-milligrams-per-liter-chloride front in the hypothetical scenarios to the baseline scenario indicated that, in general, increases to pumping in the well field area increased the rate of simulated chloride movement toward the well field area by as much as 150 feet per year. Reductions in pumping slowed the advance of chloride toward the well field by as much as 210 feet per year, although reductions did not stop the movement of chloride toward the well field, including when pumping rates were eliminated. If pumping is completely discontinued, the rate of chloride movement is about 500 to 600 feet per year in the area between the Arkansas River and the southern part of the Wichita well field, and 70 to 500 feet per year in the area near Burrton with the highest rate of movement in the shallow aquifer layer. The averages of simulated water-levels in index monitoring wells in the Wichita well field at the end of 2008 were calculated for each scenario. Compared to the baseline scenario, the average simulated water level was 5.05 feet higher for the no pumping scenario, 4.72 feet lower for the double Wichita pumping with existing irrigation scenario, 2.49 feet higher for the no irrigation pumping with existing Wichita pumping scenario, 1.53 feet lower for the double Wichita pumping with no irrigation scenario, and 0.48 feet higher for the increased Phase 1 artificial recharge scenario. The groundwater flow was simulated with a preexisting groundwater-flow model, which was not altered to calibrate the solute-transport model to observed chloride-concentration data. Therefore, some areas in the model had poor fit between simulated chloride concentrations and observed chloride concentrations, including the area between Arkansas River and the southern part of the Wichita well field, and the Hollow-Nikkel area about 6 miles north of Burrton. Compared to the interpreted location of the 250-milligrams per liter-chloride front based on data collected in 2011, in the Arkansas River area the simulated 250-milligrams per liter-chloride front moved from the river toward the well field about twice the rate of the actual 250-milligrams per liter-chloride front in the shallow layer and about four times the rate of the actual 250-milligrams per liter-chloride front in the deep layer. Future groundwater-flow and chloride-transport modeling efforts may achieve better agreement between observed and simulated chloride concentrations in these areas by taking the chloride-transport model fit into account when adjusting parameters such as hydraulic conductivity, riverbed conductance, and effective porosity during calibration. Results of the hypothetical scenarios simulated indicate that the Burrton chloride plume will continue moving toward the well field regardless of pumping in the area and that one alternative may be to increase pumping from within the plume area to reverse the groundwater-flow gradients and remove the plume. Additionally, the results of modeling these scenarios indicate that eastward movement of the Burrton plume could be slowed by the additional artificial recharge at the Phase 1 sites and that decreasing pumping along the Arkansas River or increasing water levels could retard the movement of chloride and may prevent further encroachment into the southern part of the well field area.

Kansas↗

Description and evaluation of numerical groundwater flow models for the Edwards Aquifer, south-central Texas

A substantial number of public water system wells in south-central Texas withdraw groundwater from the karstic, highly productive Edwards aquifer. However, the use of numerical groundwater flow models to aid in the delineation of contributing areas for public water system wells in the Edwards aquifer is problematic because of the complex hydrogeologic framework and the presence of conduit-dominated flow paths in the aquifer. The U.S. Geological Survey, in cooperation with the Texas Commission on Environmental Quality, evaluated six published numerical groundwater flow models (all deterministic) that have been developed for the Edwards aquifer San Antonio segment or Barton Springs segment, or both. This report describes the models developed and evaluates each with respect to accessibility and ease of use, range of conditions simulated, accuracy of simulations, agreement with dye-tracer tests, and limitations of the models. These models are (1) GWSIM model of the San Antonio segment, a FORTRAN computer-model code that pre-dates the development of MODFLOW; (2) MODFLOW conduit-flow model of San Antonio and Barton Springs segments; (3) MODFLOW diffuse-flow model of San Antonio and Barton Springs segments; (4) MODFLOW Groundwater Availability Modeling [GAM] model of the Barton Springs segment; (5) MODFLOW recalibrated GAM model of the Barton Springs segment; and (6) MODFLOW–DCM (dual conductivity model) conduit model of the Barton Springs segment. The GWSIM model code is not commercially available, is limited in its application to the San Antonio segment of the Edwards aquifer, and lacks the ability of MODFLOW to easily incorporate newly developed processes and packages to better simulate hydrologic processes. MODFLOW is a widely used and tested code for numerical modeling of groundwater flow, is well documented, and is in the public domain. These attributes make MODFLOW a preferred code with regard to accessibility and ease of use. The MODFLOW conduit-flow model incorporates improvements over previous models by using (1) a user-friendly interface, (2) updated computer codes (MODFLOW–96 and MODFLOW–2000), (3) a finer grid resolution, (4) less-restrictive boundary conditions, (5) an improved discretization of hydraulic conductivity, (6) more accurate estimates of pumping stresses, (7) a long transient simulation period (54 years, 1947–2000), and (8) a refined representation of high-permeability zones or conduits. All of the models except the MODFLOW–DCM conduit model have limitations resulting from the use of Darcy’s law to simulate groundwater flow in a karst aquifer system where non-Darcian, turbulent flow might actually dominate. The MODFLOW–DCM conduit model is an improvement in the ability to simulate karst-like flow conditions in conjunction with porous-media-type matrix flow. However, the MODFLOW–DCM conduit model has had limited application and testing and currently (2008) lacks commercially available pre- and post-processors. The MODFLOW conduit-flow and diffuse-flow Edwards aquifer models are limited by the lack of calibration for the northern part of the Barton Springs segment (Travis County) and their reliance on the use of the calibrated hydraulic conductivity and storativity values from the calibrated Barton Springs segment GAM model. The major limitation of the Barton Springs segment GAM and recalibrated GAM models is that they were calibrated to match measured water levels and springflows for a restrictive range of hydrologic conditions, with each model having different hydraulic conductivity and storativity values appropriate to the hydrologic conditions that were simulated. The need for two different sets of hydraulic conductivity and storativity values increases the uncertainty associated with the accuracy of either set of values, illustrates the non-uniqueness of the model solution, and probably most importantly demonstrates the limitations of using a one-layer model to represent the heterogeneous hydrostratigraphic units composing the Edwards aquifer. In general, the best matches or agreement between groundwater flow directions inferred by numerical model simulation, and by dye-tracer tests, are observed where model outputs accurately reproduce the configuration of the potentiometric surface with regard to the positions of major and minor groundwater troughs and divides. None of the models, with the possible exception of the MODFLOW–DCM conduit model, has a documented capability to accurately simulate travel times for conduit-dominated velocities in the Edwards aquifer. Public water system assessments of wells in the Barton Springs segment of the Edwards aquifer, and elsewhere where conduit-flow conditions are thought to dominate aquifer hydraulic behavior, might be enhanced by use of either the MODFLOW–DCM model or the newly developed U.S. Geological Survey MODFLOW Conduit-Flow Process module for MODFLOW–2005, because each incorporates a type of dual or triple hydraulic conductivity approach and has the capability to explicitly simulate turbulent flow and conduit hydraulic characteristics.

Texas↗

Status and understanding of groundwater quality in the Sacramento Metropolitan Domestic-Supply Aquifer study unit, 2017—California GAMA Priority Basin Project

Groundwater quality in the Sacramento Metropolitan Domestic-Supply Aquifer study unit (SacMetro-DSA) was studied from August to November 2017 as part of the second phase of the Priority Basin Project of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is in parts of Amador, Placer, Sacramento, and Sutter Counties, and the extent of the study unit was defined by the location of three California Department of Water Resources groundwater subbasins: the North American, the South American, and the Cosumnes. The SacMetro-DSA focused on groundwater resources used for domestic drinking-water supply, which generally correspond to shallower parts of aquifer systems than those of groundwater resources used for public drinking water supply in the same area. The assessments characterized the quality of untreated groundwater, not the quality of drinking water. This study included two components: (1) a status assessment, which characterized the status of the quality of the groundwater resources used for domestic supply and (2) an understanding assessment, which evaluated the natural and human factors potentially affecting water quality in those resources. The first component of this study—the status assessment—was based on water-quality data collected from 49 sites sampled by the U.S. Geological Survey for the GAMA Priority Basin Project in 2017. The samples were analyzed for volatile organic compounds, pesticides, and naturally present inorganic constituents, such as major ions and trace elements. To provide context, concentrations of constituents measured in groundwater were compared to U.S. Environmental Protection Agency and California State Water Resources Control Board Division of Drinking Water regulatory and non-regulatory benchmarks for drinking-water quality. The status assessment used a grid-based method to estimate the proportion of the groundwater resources that had concentrations of water-quality constituents approaching or above benchmark concentrations. This method provides statistically unbiased results at the study-area scale and permits comparisons to other GAMA Priority Basin Project study areas. The second component of this study—the understanding assessment—identified the natural and human factors that potentially affect groundwater quality by evaluating land-use characteristics, groundwater age, and geochemical and hydrologic conditions of the domestic-supply aquifer and related these data to constituents identified in the status assessment for further evaluation. In the SacMetro-DSA study unit, arsenic was the only inorganic constituent detected above health-based benchmarks and was detected in 10 percent of the domestic-supply aquifer system. Inorganic constituents were detected above the non-health-based California State Water Resources Control Board—Division of Drinking Water secondary maximum contaminant levels (SMCL-CA) in 16 percent of the system. The inorganic constituents detected above the SMCL-CA were chloride, iron, manganese, and total dissolved solids (TDS). Organic constituents (volatile organic compounds and pesticides) with health-based benchmarks were not detected above health-based benchmarks; however, chloroform was detected at concentrations higher than 10 percent of the health-based benchmark (80 micrograms per liter) in 2 percent of the domestic-supply aquifer system. Of the 310 organic constituents analyzed, 16 constituents were detected; however, only bentazon and chloroform had detection frequencies greater than 10 percent. Inorganic constituents with health-based benchmarks that were evaluated in the understanding assessment included arsenic and hexavalent chromium. Arsenic and hexavalent chromium are natural constituents of aquifer sediments in the study unit and did not appear to be influenced by anthropogenic processes; rather, the presence of arsenic and hexavalent chromium appeared to be related to geochemical conditions controlled by oxidation–reduction reactions in the aquifer system. Naturally occurring inorganic constituents with SMCL-CAs evaluated in the understanding assessment were the trace elements iron and manganese, the major ion chloride, and TDS. Like arsenic and hexavalent chromium, the presence of iron and manganese was most strongly related to geochemical conditions in the aquifer system, specifically reducing conditions, which were most common near the western edge of the study unit close to the Sacramento River. Concentrations of chloride and TDS are indicators of salinity and were correlated with variables related to well location and included redox, agricultural land use, and elevation. Chloride and TDS were positively correlated to reducing conditions, and agricultural land use was negatively correlated to elevation and well depth. Observed correlations among variables were likely driven by the characteristics of the western part of the study unit, such as its higher proportion of agricultural land use and its relatively low elevation. A large portion of the western edge of the study unit is located in the center of the Sacramento Valley, defined by the location of the Sacramento River. The special-interest constituent perchlorate, also included in the understanding assessment, has natural and anthropogenic sources. Perchlorate was detected frequently and at moderate relative concentrations. In some areas of the study unit, concentrations of perchlorate were higher than what might be expected in nature; therefore, anthropogenic introduction of perchlorate or anthropogenically induced migration of native perchlorate could be occurring.

California↗

Groundwater-flow model of the Ozark Plateaus aquifer system, northwestern Arkansas, southeastern Kansas, southwestern Missouri, and northeastern Oklahoma

To assess the effect that increased water use is having on the long-term availability of groundwater within the Ozark Plateaus aquifer system, a groundwater-flow model was developed using MODFLOW 2000 for a model area covering 7,340 square miles for parts of Arkansas, Kansas, Missouri, and Oklahoma. Vertically the model is divided into five units. From top to bottom these units of variable thickness are: the Western Interior Plains confining unit, the Springfield Plateau aquifer, the Ozark confining unit, the Ozark aquifer, and the St. Francois confining unit. Large mined zones contained within the Springfield Plateau aquifer are represented in the model as extensive voids with orders-of-magnitude larger hydraulic conductivity than the adjacent nonmined zones. Water-use data were compiled for the period 1960 to 2006, with the most complete data sets available for the period 1985 to 2006. In 2006, total water use from the Ozark aquifer for Missouri was 87 percent (8,531,520 cubic feet per day) of the total pumped from the Ozark aquifer, with Kansas at 7 percent (727,452 cubic feet per day), and Oklahoma at 6 percent (551,408 cubic feet per day); water use for Arkansas within the model area was minor. Water use in the model from the Springfield Plateau aquifer in 2005 was specified from reported and estimated values as 569,047 cubic feet per day. Calibration of the model was made against average water-level altitudes in the Ozark aquifer for the period 1980 to 1989 and against waterlevel altitudes obtained in 2006 for the Springfield Plateau and Ozark aquifers. Error in simulating water-level altitudes was largest where water-level altitude gradients were largest, particularly near large cones of depression. Groundwater flow within the model area occurs generally from the highlands of the Springfield Plateau in southwestern Missouri toward the west, with localized flow occurring towards rivers and pumping centers including the five largest pumping centers near Joplin, Missouri; Carthage, Missouri; Noel, Missouri; Pittsburg, Kansas; and Miami, Oklahoma. Hypothetical scenarios involving various increases in groundwater-pumping rates were analyzed with the calibrated groundwater-flow model to assess changes in the flow system from 2007 to the year 2057. Pumping rates were increased between 0 and 4 percent per year starting with the 2006 rates for all wells in the model. Sustained pumping at 2006 rates was feasible at the five pumping centers until 2057; however, increases in pumping resulted in dewatering the aquifer and thus pumpage increases were not sustainable in Carthage and Noel for the 1 percent per year pumpage increase and greater hypothetical scenarios, and in Joplin and Miami for the 4 percent per year pumpage increase hypothetical scenarios. Zone-budget analyses were performed to assess the groundwater flow into and out of three zones specified within the Ozark-aquifer layer of the model. The three zones represented the model part of the Ozark aquifer in Kansas (zone 1), Oklahoma (zone 2), and Missouri and Arkansas (zone 3). Groundwater pumping causes substantial reductions in water in storage and induces flow through the Ozark confining unit for all hypothetical scenarios evaluated. Net simulated flow in 2057 from Kansas (zone 1) to Missouri (zone 3) ranges from 74,044 cubic feet per day for 2006 pumping rates (hypothetical scenario 1) to 625,319 cubic feet per day for a 4 percent increase in pumping per year (hypothetical scenario 5). Pumping from wells completed in the Ozark aquifer is the largest component of flow out of zone 3 in Missouri and Arkansas, and varies between 88 to 91 percent of the total flow out of zone 3 for all of the hypothetical scenarios. The largest component of flow into Oklahoma (zone 2) comes from the overlying Ozark confining unit, which is consistently about 45 percent of the total. Flow from the release of water in storage, from general-head boundaries, and from zones 1 and 3 is considerably smaller values that range from 3 to 22 percent of the total flow into zone 2. The largest flow out of the Oklahoma part of the model occurs from pumping from wells and ranges from 52 to 69 percent of the total.

Arkansas, Kansas, Missouri, Oklahoma↗

Hydrogeology and simulation of groundwater flow in the Arbuckle-Simpson aquifer, south-central Oklahoma

The Arbuckle-Simpson aquifer in south-central Oklahoma provides water for public supply, farms, mining, wildlife conservation, recreation, and the scenic beauty of springs, streams, and waterfalls. Proposed development of water supplies from the aquifer led to concerns that large-scale withdrawals of water would cause decreased flow in rivers and springs, which in turn could result in the loss of water supplies, recreational opportunities, and aquatic habitat. The Oklahoma Water Resources Board, in collaboration with the Bureau of Reclamation, the U.S. Geological Survey, Oklahoma State University, and the University of Oklahoma, studied the aquifer to provide the Oklahoma Water Resources Board the scientific information needed to determine the volume of water that could be withdrawn while protecting springs and streams. The U.S. Geological Survey, in coopertion with the Oklahoma Water Resources Board, did a study to describe the hydrogeology and simulation of groundwater flow of the aquifer. The outcrop of the Arbuckle-Simpson aquifer covers an area of about 520 square miles in Carter, Coal, Johnston, Murray, and Pontotoc Counties. Three subdivisions of the aquifer outcrop were designated for this study: the eastern, central, and western Arbuckle-Simpson aquifer. This study emphasized the eastern Arbuckle-Simpson aquifer because it is the largest part of the aquifer by area and volume; most groundwater withdrawals are from the eastern Arbuckle-Simpson aquifer; and the largest (by flow) streams and springs sourced from the aquifer are on the eastern Arbuckle-Simpson aquifer. The aquifer lies in an uplifted area commonly referred to as the Arbuckle Mountains, which is characterized by great thicknesses of mostly carbonate rocks, uplifts, folded structures, and large fault displacements. The Arbuckle-Simpson aquifer is contained in three major rock units of Late Cambrian to Middle Ordovician age: the Timbered Hills, Arbuckle, and Simpson Groups. The aquifer is underlain by low-permeability Cambrian and Proterozoic igneous and metamorphic rocks, and is confined above by younger sedimentary rocks of various ages in areas where the top of the aquifer dips below the surface. The major part of the Arbuckle-Simpson aquifer is the Arbuckle Group, which consists of as much as 6,700 feet of limestone in the western Arbuckle-Simpson aquifer, but which thins to an estimated 3,000 feet of predominantly dolostone in the eastern Arbuckle-Simpson aquifer. Water is obtained from cavities, solution channels, fractures, and intercrystalline porosity in the limestone and dolostone. The overlying Simpson Group, consisting of sandstones, shales, and limestones, is as much as 2,300 feet thick in the western Arbuckle-Simpson aquifer, but generally is less than 1,000 feet thick in the eastern aquifer. Water in the Simpson Group is stored primarily in pore spaces between the sand grains in the sandstones. A digital, three-dimensional geologic framework model was constructed to define the geometric relations of fault blocks and subsurface rock units across complex fault zones of the eastern Arbuckle-Simpson aquifer. Geologic data for the model were obtained from 126 drill holes; stratigraphic contacts and faults defined from a digitized version of the surface geologic map; and fault geometry, stratigraphic thickness, and information compiled from geologic and hydrogeologic reports and maps. Groundwater in the aquifer moves from areas of high head (altitude) to areas of low head along streams and springs. The potentiometric surface in the eastern Arbuckle-Simpson aquifer generally slopes from a topographic high from northwest to the southeast, indicating that regional groundwater flow is predominantly toward the southeast. Freshwater is known to extend beyond the aquifer outcrop near the City of Sulphur, Oklahoma, and Chickasaw National Recreation Area, where groundwater flows west from the outcrop of the eastern Arbuckle-Simpson aquifer and becomes confin

Scientific Investigations Report↗

Volcanic aquifers of Hawaiʻi—Contributions to assessing groundwater availability on Kauaʻi, Oʻahu, and Maui

The volcanic aquifers of the Hawaiian Islands supply water to 1.46 million residents, diverse industries, and a large component of the U.S. military in the Pacific. Groundwater also supplies fresh water that supports ecosystems in streams and near the coast. Hawaii’s aquifers are remarkably productive given their small size, but the capacity of the islands to store fresh groundwater is limited because each island is surrounded by seawater, and salt water underlies much of the fresh groundwater. The amount of fresh groundwater available for human use from Hawai‘i’s volcanic aquifers is constrained by the consequences of groundwater withdrawal. Restrictions placed on these consequences can translate to limitations on groundwater availability. Changes in recharge resulting from changes in land cover or climate can alter the effect of withdrawals. This study uses numerical models of the volcanic aquifers of the islands of Kaua‘i, O‘ahu, and Maui to quantify the consequences of historical and plausible future withdrawals and changes in recharge. The study compares the results of model simulations of multiple scenarios of historical and projected future withdrawal and recharge. Results of the simulations using the groundwater models of the islands of Kaua‘i, O‘ahu, and Maui have implications for other islands in Hawai‘i. Since the first modern water well was drilled in Hawai‘i in 1879, total groundwater withdrawals on Kaua‘i, O‘ahu, and Maui have risen to nearly 400 million gallons per day. Model simulations indicate that these withdrawals have caused reductions in groundwater discharge to streams and springs, reductions in groundwater discharge to the ocean, changes in subsurface flow between sectors within an island, lowering of groundwater levels, and rise of the interface between fresh water and salt water in the aquifers. Future increases in withdrawals will increase the severity of the consequences. Changes in recharge can alter the effect of withdrawals—increases in recharge can offset the consequences of withdrawals, whereas decreases in recharge can exacerbate the effects of withdrawals. This study quantifies the consequences of withdrawals for past and plausible future circumstances. The models can be used to test other circumstances. Limits placed on the consequences of withdrawals—such as restrictions to protect stream or coastal ecosystems that rely on groundwater discharge and limitations on water-level decline and rise of the freshwater-saltwater interface to protect the productivity of existing wells—can translate to limits on groundwater availability from Hawai‘i’s volcanic aquifers. Setting acceptable limits to the consequences of groundwater withdrawal is also a critical part of assessing groundwater availability. Once these limits are set, numerical models can be used to quantify the amount of water that can be withdrawn within those limits and thereby inform management decisions that seek to balance the need to limit the consequences of groundwater withdrawals with the need to develop water for human use.

Hawaii↗