USGS Science⌕ Search

Geology topics

Groundwater & aquifers

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

At least 1,099 records · Page 61Linked to original sources

Hydrogeology and steady-state numerical simulation of groundwater flow in the Lost Creek Designated Ground Water Basin, Weld, Adams, and Arapahoe Counties, Colorado

The Lost Creek Designated Ground Water Basin (Lost Creek basin) is an important alluvial aquifer for irrigation, public supply, and domestic water uses in northeastern Colorado. Beginning in 2005, the U.S. Geological Survey, in cooperation with the Lost Creek Ground Water Management District and the Colorado Water Conservation Board, collected hydrologic data and constructed a steady-state numerical groundwater flow model of the Lost Creek basin. The model builds upon the work of previous investigators to provide an updated tool for simulating the potential effects of various hydrologic stresses on groundwater flow and evaluating possible aquifer-management strategies. As part of model development, the thickness and extent of regolith sediments in the basin were mapped, and data were collected concerning aquifer recharge beneath native grassland, nonirrigated agricultural fields, irrigated agricultural fields, and ephemeral stream channels. The thickness and extent of regolith in the Lost Creek basin indicate the presence of a 2- to 7-mile-wide buried paleovalley that extends along the Lost Creek basin from south to north, where it joins the alluvial valley of the South Platte River valley. Regolith that fills the paleovalley is as much as about 190 ft thick. Average annual recharge from infiltration of precipitation on native grassland and nonirrigated agricultural fields was estimated by using the chloride mass-balance method to range from 0.1 to 0.6 inch, which represents about 1-4 percent of long-term average precipitation. Average annual recharge from infiltration of ephemeral streamflow was estimated by using apparent downward velocities of chloride peaks to range from 5.7 to 8.2 inches. Average annual recharge beneath irrigated agricultural fields was estimated by using passive-wick lysimeters and a water-balance approach to range from 0 to 11.3 inches, depending on irrigation method, soil type, crop type, and the net quantity of irrigation water applied. Estimated average annual recharge beneath irrigated agricultural fields represents about 0-43 percent of net irrigation. The U.S. Geological Survey modular groundwater modeling program, MODFLOW-2000, was used to develop a steady-state groundwater flow model of the Lost Creek basin. Groundwater in the basin is simulated generally to flow from the basin margins toward the center of the basin and northward along the paleovalley. The largest source of inflow to the model occurs from recharge beneath flood- and sprinkler-irrigated agricultural fields (14,510 acre-feet per year [acre-ft/yr]), which represents 39.7 percent of total simulated inflow. Other substantial sources of inflow to the model are recharge from precipitation and stream-channel infiltration in nonirrigated areas (13,810 acre-ft/yr) seepage from Olds Reservoir (4,280 acre-ft/yr), and subsurface inflow from ditches and irrigated fields outside the model domain (2,490 acre-ft/yr), which contribute 37.7, 11.7, and 6.8 percent, respectively, of total inflow. The largest outflow from the model occurs from irrigation well withdrawals (26,760 acre-ft/yr), which represent 73.2 percent of total outflow. Groundwater discharge (6,640 acre-ft/yr) at the downgradient end of the Lost Creek basin represents 18.2 percent of total outflow, and evapotranspiration (3,140 acre-ft/yr) represents about 8.6 percent of total outflow.

Scientific Investigations Report↗

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

This report documents selected age data interpreted from measured concentrations of environmental tracers in groundwater from 1,399 National Water-Quality Assessment (NAWQA) Program groundwater sites across the United States. The tracers of interest were chlorofluorocarbons (CFCs), sulfur hexafluoride (SF 6 ), and tritium/helium-3 ( 3 H/ 3 He). Tracer data compiled for this analysis primarily were from wells representing two types of NAWQA groundwater studies—Land-Use Studies (shallow wells, usually monitoring wells, in recharge areas under dominant land-use settings) and Major-Aquifer Studies (wells, usually domestic supply wells, in principal aquifers and representing the shallow, used resource). Reference wells (wells representing groundwater minimally impacted by anthropogenic activities) associated with Land-Use Studies also were included. Tracer samples were collected between 1992 and 2005, although two networks sampled from 2006 to 2007 were included because of network-specific needs. Tracer data from other NAWQA Program components (Flow System Studies, which are assessments of processes and trends along groundwater flow paths, and various topical studies) were not compiled herein. Tracer data from NAWQA Land-Use Studies and Major-Aquifer Studies that previously had been interpreted and published are compiled herein (as piston-flow ages), but have not been reinterpreted. Tracer data that previously had not been interpreted and published are evaluated using documented methods and compiled with aqueous concentrations, equivalent atmospheric concentrations (for CFCs and SF 6 ), estimates of tracer-based piston-flow ages, and selected ancillary data, such as redox indicators, well construction, and major dissolved gases (N 2 , O 2 , Ar, CH 4 , and CO 2 ). Tracer-based piston-flow ages documented in this report are simplistic representations of the tracer data. Tracer-based piston-flow ages are a convenient means of conceptualizing groundwater age. However, the piston-flow model is based on the potentially limiting assumptions that tracer transport is advective and that no mixing occurs. Additional uncertainties can arise from tracer degradation, sorption, contamination, or fractionation; terrigenic (natural) sources of tracers; spatially variable atmospheric tracer concentrations; and incomplete understanding of mechanisms of recharge or of the conditions under which atmospheric tracers were partitioned to recharge. The effects of some of these uncertainties are considered herein. For example, degradation, contamination, or fractionation often can be identified or inferred. However, detailed analysis of the effects of such uncertainties on the tracer-based piston-flow ages is constrained by sparse data and an absence of complementary lines of evidence, such as detailed solute transport simulations. Thus, the tracer-based piston-flow ages compiled in this report represent only an initial interpretation of the tracer data.

Scientific Investigations Report↗

Water-quality assessment of the Cambrian-Ordovician aquifer system in the northern Midwest, United States

This report provides a regional assessment of groundwater quality of the Cambrian-Ordovician aquifer system, based primarily on raw water samples collected by the NAWQA Program during 1995 through 2007. The NAWQA Program has published findings in local study-unit reports encompassing parts of the Cambrian-Ordovician aquifer system. Data collected from the aquifer system were used in national synthesis reports on selected topics such as specific water-quality constituent classes, well type, or aquifer material; however, a synthesis of groundwater quality at the principal aquifer scale has not been completed and is therefore the major purpose of this report. Water samples collected by the NAWQA Program were analyzed for various classes of characteristics including physical properties, major ions, trace elements, nutrients and dissolved organic carbon, radionuclides (tritium, radon, and radium), pesticides, and volatile organic compounds. Subsequent sections of this report provide discussions on these classes of characteristics. The assessment objectives of this report are to (1) summarize constituent concentrations and compare them to human-health benchmarks and non-health guidelines; (2) determine the geographic distribution of constituent concentrations and relate them to various factors such as confining conditions, well type, land use, and groundwater age; and (3) evaluate near-decadal-scale changes in nitrate concentrations and pesticide detections. The most recent sample collected from each well by the NAWQA Program was used for most analyses. Near-decadal-scale changes in nitrate concentrations and pesticide detections were evaluated for selected well networks by using the most recent sample from each well and comparing it to the results from a sample collected 7 or 11 years earlier. Because some of the NAWQA well networks provide a limited areal coverage of the aquifer system, data for raw water samples from other USGS sources and state agencies were included to expand the data coverage into areas between the NAWQA well networks and into northeastern Missouri. Many of the maps in this report that show concentrations of selected constituents include data from other sources to expand on the geographic area covered by the NAWQA data.

Midwest;Cambrian-ordovician Aquifer System↗

Iodine-129 in the eastern Snake River Plain aquifer at and near the Idaho National Laboratory, Idaho, 2010-12

From 1953 to 1988, approximately 0.941 curies of iodine-129 ( 129 I) were contained in wastewater generated at the Idaho National Laboratory (INL) with almost all of this wastewater discharged at or near the Idaho Nuclear Technology and Engineering Center (INTEC). Most of the wastewater containing 129 I was discharged directly into the eastern Snake River Plain (ESRP) aquifer through a deep disposal well until 1984; lesser quantities also were discharged into unlined infiltration ponds or leaked from distribution systems below the INTEC. During 2010–12, the U.S. Geological Survey in cooperation with the U.S. Department of Energy collected groundwater samples for 129 I from 62 wells in the ESRP aquifer to track concentration trends and changes for the carcinogenic radionuclide that has a 15.7 million-year half-life. Concentrations of 129 I in the aquifer ranged from 0.0000013±0.0000005 to 1.02±0.04 picocuries per liter (pCi/L), and generally decreased in wells near the INTEC, relative to previous sampling events. The average concentration of 129 I in groundwater from 15 wells sampled during four different sample periods decreased from 1.15 pCi/L in 1990–91 to 0.173 pCi/L in 2011–12. All but two wells within a 3-mile radius of the INTEC showed decreases in concentration, and all but one sample had concentrations less than the U.S. Environmental Protection Agency maximum contaminant level of 1 pCi/L. These decreases are attributed to the discontinuation of disposal of 129 I in wastewater and to dilution and dispersion in the aquifer. The decreases in 129 I concentrations, in areas around INTEC where concentrations increased between 2003 and 2007, were attributed to less recharge near INTEC either from less flow in the Big Lost River or from less local snowmelt and anthropogenic sources. Although wells near INTEC sampled in 2011–12 showed decreases in 129 I concentrations compared with previously collected data, some wells south and east of the Central Facilities Area, near the site boundary, and south of the INL showed small increases. These slight increases are attributed to variable discharge rates of wastewater that eventually moved to these well locations as a pulse of water from a particular disposal period. Wells sampled for the first time around the Naval Reactors Facility had 129 I concentrations slightly greater than background concentrations in the ESRP aquifer. These concentrations are attributed to possible leakage from landfills at the Naval Reactors Facility or seepage from air emission deposits from INTEC, or both. In 2012, the U.S. Geological Survey collected discrete groundwater samples from 25 zones in 11 wells equipped with multilevel monitoring systems to help define the vertical distribution of 129 I in the aquifer. Concentrations ranged from 0.000006±0.000004 to 0.082±0.003 pCi/L. Two new wells completed in 2012 showed variability of up to one order of magnitude of concentrations of 129 I among various zones. Two other wells showed similar concentrations of 129 I in all three zones sampled. Concentrations were well less than the maximum contaminant level in all zones.

Idaho↗

Land subsidence, groundwater levels, and geology in the Coachella Valley, California, 1993-2010

Land subsidence associated with groundwater-level declines has been investigated by the U.S. Geological Survey in the Coachella Valley, California, since 1996. Groundwater has been a major source of agricultural, municipal, and domestic supply in the valley since the early 1920s. Pumping of groundwater resulted in water-level declines as much as 15 meters (50 feet) through the late 1940s. In 1949, the importation of Colorado River water to the southern Coachella Valley began, resulting in a reduction in groundwater pumping and a recovery of water levels during the 1950s through the 1970s. Since the late 1970s, demand for water in the valley has exceeded deliveries of imported surface water, resulting in increased pumping and associated groundwater-level declines and, consequently, an increase in the potential for land subsidence caused by aquifer-system compaction. Global Positioning System (GPS) surveying and Interferometric Synthetic Aperture Radar (InSAR) methods were used to determine the location, extent, and magnitude of the vertical land-surface changes in the southern Coachella Valley during 1993–2010. The GPS measurements taken at 11 geodetic monuments in 1996 and in 2010 in the southern Coachella Valley indicated that the elevation of the land surface changed –136 to –23 millimeters (mm) ±54 mm (–0.45 to –0.08 feet (ft) ±0.18 ft) during the 14-year period. Changes at 6 of the 11 monuments exceeded the maximum expected uncertainty of ±54 mm (±0.18 ft) at the 95-percent confidence level, indicating that subsidence occurred at these monuments between June 1996 and August 2010. GPS measurements taken at 17 geodetic monuments in 2005 and 2010 indicated that the elevation of the land surface changed –256 to +16 mm ±28 mm (–0.84 to +0.05 ft ±0.09 ft) during the 5-year period. Changes at 5 of the 17 monuments exceeded the maximum expected uncertainty of ±28 mm (±0.09 ft) at the 95-percent confidence level, indicating that subsidence occurred at these monuments between August 2005 and August 2010. At each of these five monuments, subsidence rates were about the same between 2005 and 2010 as between 2000 and 2005. InSAR measurements taken between June 27, 1995, and September 19, 2010, indicated that the land surface subsided from about 220 to 600 mm (0.72 to 1.97 ft) in three areas of the Coachella Valley: near Palm Desert, Indian Wells, and La Quinta. In Palm Desert, the average subsidence rates increased from about 39 millimeters per year (mm/yr), or 0.13 foot per year (ft/yr), during 1995–2000 to about 45 mm/yr (0.15 ft/yr) during 2003–10. In Indian Wells, average subsidence rates for two subsidence maxima were fairly steady at about 34 and 26 mm/yr (0.11 and 0.09 ft/yr) during both periods; for the third maxima, average subsidence rates increased from about 14 to 19 mm/yr (0.05 to 0.06 ft/yr) from the first to the second period. In La Quinta, average subsidence rates for five selected locations ranged from about 17 to 37 mm/yr (0.06 to 0.12 ft/yr) during 1995–2000; three of the locations had similar rates during 2003–mid-2009, while the other two locations had increased subsidence rates. Decreased subsidence rates were calculated throughout the La Quinta subsidence area during mid-2009–10, however, and uplift was observed during 2010 near the southern extent of this area. Water-level measurements taken at wells near the subsiding monuments and in the three subsiding areas shown by InSAR generally indicated that the water levels fluctuated seasonally and declined annually from the early 1990s, or earlier, to 2010; some water levels in 2010 were at the lowest levels in their recorded histories. An exception to annually declining water levels in and near subsiding areas was observed beginning in mid-2009 in the La Quinta subsidence area, where recovering water levels coincided with increased recharge operations at the Thomas E. Levy Recharge Facility; decreased pumpage also could cause groundwater levels to recover. Subsidence concomitant with declining water levels and land-surface uplift concomitant with recovering water levels indicate that aquifer-system compaction could be causing subsidence. If the stresses imposed by the historically lowest water levels exceeded the preconsolidation stress, the aquifer-system compaction and associated land subsidence could be permanent.

California↗

Evaluation of effects of groundwater withdrawals at the proposed Allen combined-cycle combustion turbine plant, Shelby County, Tennessee

The Mississippi Embayment Regional Aquifer Study groundwater-flow model was used to simulate the potential effects of future groundwater withdrawals at the proposed Allen combined-cycle combustion turbine plant in Shelby County, Tennessee. The scenario used in the simulation consisted of a 30-year average withdrawal period followed by a 30-day maximum withdrawal period. Effects of withdrawals at the Allen plant site on the Mississippi embayment aquifer system were evaluated by comparing the difference in simulated water levels in the aquifers at the end of the 30-year average withdrawal period and at the end of the scenario to a base case without the Allen combined-cycle combustion turbine plant withdrawals. Simulated potentiometric surface declines in the Memphis aquifer at the Allen plant site were about 7 feet at the end of the 30-year average withdrawal period and 11 feet at the end of the scenario. The affected area of the Memphis aquifer at the Allen plant site as delineated by the 4-foot potentiometric surface-decline contour was 2,590 acres at the end of the 30-year average withdrawal period and 11,380 acres at the end of the scenario. Simulated declines in the underlying Fort Pillow aquifer and overlying shallow aquifer were both less than 1 foot at the end of the 30-year average withdrawal period and the end of the scenario.

Arkansas, Mississippi, Tennessee↗

Geophysical surveys, hydrogeologic characterization, and groundwater flow model for the Truxton basin and Hualapai Plateau, northwestern Arizona

This is the third and final report in a series that describe the groundwater resources of the Hualapai Indian Reservation. These reports document the findings of a comprehensive groundwater study conducted on the reservation and adjacent areas from 2015 through 2018 by the U.S. Geological Survey in cooperation with the Bureau of Reclamation. The first report described the hydrologic framework and characterization of the Truxton aquifer on the Hualapai Indian Reservation (Bills and Macy, 2016). The second report described the hydrogeologic characterization of the Hualapai Plateau part of the reservation (Mason, Macy, and others, 2020). This report includes five chapters. Chapter A (Mason, Knight, and others, 2020) is a summary of this multichapter volume and briefly describes the study area. Chapter B (Mason, Bills, and Macy, 2020) describes the geology and hydrology of the Truxton basin and Hualapai Plateau. Chapter C (Kennedy, 2020) describes the results of a gravity geophysical survey of the Truxton basin. Chapter D (Ball, 2020) describes the findings of an airborne electromagnetic survey of the Truxton aquifer and Hualapai Plateau. Chapter E (Knight, 2020) describes the results of a transient groundwater model created for the entire Truxton aquifer both on and off the reservation. The groundwater-flow model is used to estimate projected groundwater levels based on future groundwater withdrawal scenarios.

Arizona↗

Surface-water and groundwater interactions in the Big Lost River, south-central Idaho

The Big Lost River of south-central Idaho interacts with the underlying aquifer by gaining and losing streamflow throughout various areas in the Big Lost River Valley. Surface-water and groundwater resources are used throughout the valley to sustain domestic, agricultural, and livestock needs. The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources, evaluated streamflow gains and losses by differential streamgaging in the lower Big Lost River, Idaho, during four measurement events: March 27–28, 2019; October 16–17, 2019; October 6–7, 2020; and March 30, 2021. This report presents and analyzes streamflow measurement and uncertainty data from each measurement event to describe surface-water/groundwater interactions. This report is the second chapter of a multi-chapter volume that characterizes water resources in the Big Lost River Basin. During the four measurement events, 100 streamflow measurements were made at 46 unique sites on the Big Lost River, James Creek, and diversions or tributaries between Mackay Reservoir near Mackay and Arco, Idaho. Aquifer lithology and dimensions affected spatial patterns of streamflow gains and losses between the upper, middle, and lower reaches; changes in water supply, groundwater levels, and surface-water management affected seasonal differences within reaches. In the upper reach of the Big Lost River, streamflow losses and gains were greater during the wetter 2019 events and lesser during the drier 2020 and 2021 events. The middle reach includes the largest losses from the Big Lost River to groundwater; these losses occurred in the Darlington Sinks where 42 percent or more of streamflow was lost as the aquifer widens and groundwater deepens. These results suggest that changing surface-water supply, irrigation use, and recharge affect interannual groundwater levels and, in turn, affect patterns of streamflow gains and losses in the middle reach. Finally, surface-water management is the primary control on surface-water/groundwater interactions in the lower reach. Overall patterns of streamflow gains and losses in this study generally were consistent with previous reports. However, paired with the related hydrogeologic framework and water budget, this investigation provides new insights into how hydrogeologic conditions and interannual variability in water supply, groundwater levels, and surface-water management affect surface-water/groundwater interactions in the Big Lost River Valley.

Idaho↗

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado↗

Aquifer storage change and storage properties, Rio Rancho, New Mexico, 2019–23

To better understand changes in groundwater storage and groundwater elevations, the U.S. Geological Survey, in cooperation with the City of Rio Rancho, New Mexico, carried out a multiyear groundwater monitoring project. Groundwater-level data were collected at 27 locations, including sites having multiple wells screened at different depths and those having long-term records. A repeat microgravity network of 20 stations was established, and surveys were carried out three times per year. The microgravity method provides a direct, quantitative measurement of mass change caused by aquifer filling or draining. Data collected during the 2019–23 study period indicate generally stable groundwater conditions, with small fluctuations in groundwater levels (increasing at some wells, declining at others), and small declines in groundwater storage over the period of record at most gravity locations (average = −0.33 foot of water per year). The discrepancy between the water-level and microgravity data may have been caused by a loss of soil moisture in the unsaturated zone, which is as much as 1,000 feet thick in some areas. At the Rio Rancho Advanced Water Treatment Facility, where the city recharges water through direct injection, there may be seasonal correlations in storage related to injection but no longer-term accumulation of recharged water in the immediate vicinity of the facility, indicating water is moving efficiently into the aquifer.

New Mexico↗

Groundwater sustainability and land subsidence in California’s Central Valley

The Central Valley of California is one of the most prolific agricultural regions in the world. Agriculture is reliant on the conjunctive use of surface-water and groundwater. The lack of available surface-water and land-use changes have led to pumping-induced groundwater-level and storage declines, land subsidence, changes to streamflow and the environment, and the degradation of water quality. As a result, in part, the Sustainable Groundwater Management Act (SGMA) was developed. An examination of the components of SGMA and contextualizing regional model applications within the SGMA framework was undertaken to better understand and quantify many of the components of SGMA. Specifically, the U.S. Geological Survey (USGS) updated the Central Valley Hydrologic Model (CVHM) to assess hydrologic system responses to climatic variation, surface-water availability, land-use changes, and groundwater pumping. MODFLOW-OWHM has been enhanced to simulate the timing of land subsidence and attribute its inelastic and elastic portions. In addition to extending CVHM through 2019, the new version, CVHM2, includes several enhancements as follows: managed aquifer recharge (MAR), pumping with multi-aquifer wells, inflows from ungauged watersheds, and more detailed water-balance subregions, streamflow network, diversions, tile drains, land use, aquifer properties, and groundwater level and land subsidence observations. Combined with historical approximations, CVHM2 estimates approximately 158 km 3 of storage loss in the Central Valley from pre-development to 2019. About 15% of the total storage loss is permanent loss of storage from subsidence that has caused damage to infrastructure. Climate extremes will likely complicate the efforts of water managers to store more water in the ground. CVHM2 can provide data in the form of aggregated input datasets, simulate climatic variations and changes, land-use changes or water management scenarios, and resulting changes in groundwater levels, storage, and land subsidence to assist decision-makers in the conjunctive management of water supplies.

California↗

Hydrologic budget of the Harney Basin groundwater system, southeastern Oregon

Groundwater-level declines and limited quantitative knowledge of the groundwater-flow system in the Harney Basin prompted a cooperative study between the U.S. Geological Survey and the Oregon Water Resources Department to evaluate the groundwater-flow system and budget. This report provides a hydrologic budget of the Harney Basin groundwater system that includes separate groundwater budgets for upland and lowland areas to avoid double counting water that recharges in the uplands, discharges to streams and springs in the uplands, flows downstream to the lowlands, and recharges the lowland groundwater system. Lowlands generally represent the conterminous valleys within the center of the basin, including floodplains of the major streams and uplands represent all other areas in the basin. The upland groundwater budget is minimally affected by groundwater development and generally represents the budget of the natural system. In upland areas during 1982–2016, mean-annual recharge totaled 288,000 acre-feet (acre-ft) and mean-annual discharge totaled 239,000 acre-ft, resulting in a net recharge of 49,000 acre-ft. Upland groundwater recharge occurs as infiltration of precipitation and snowmelt and was estimated using the USGS Soil-Water-Balance model calibrated to estimates of runoff, evapotranspiration (ET), base flow, and snow-water equivalent. Groundwater discharge to streams is the predominant discharge mechanism in upland areas and was estimated as 225,000 acre-feet per year (acre-ft/yr) during 1982–2016 using hydrograph separation and summer low-flow estimates in streamgaged watersheds and a linear relation between estimated streamflow and base flow in ungaged watersheds. The remaining upland discharge occurs through springs (14,000 acre-ft/yr) that either emerge downgradient of locations where groundwater discharge to streams was estimated or are routed to irrigated areas. Spring discharge was estimated as a compilation of current and historical measurements. The net upland recharge, which is 17 percent of total upland recharge, ultimately recharges lowland areas as groundwater flow from uplands to lowlands. The lowland groundwater budget for the Harney Basin represents a combination of natural conditions and human activity as more than 99 percent of groundwater development has occurred either inside or within 2 miles of the lowland boundary. In lowland areas during 1982–2016, mean annual groundwater recharge totaled 173,000 acre-ft and groundwater discharge totaled 283,000 acre-ft, indicating discharge exceeded recharge by more than 60 percent. Excluding groundwater pumping, the lowland groundwater budget is more in balance with a mean annual recharge of 165,000 acre-ft and a mean annual discharge of 131,000 acre-ft during 1982–2016. The 23-percent difference between non-pumping recharge and discharge mostly represents the cumulative uncertainty in the estimates of the various groundwater budget components but also likely includes a small reduction in natural groundwater discharge captured by pumping. Lowland groundwater is predominantly recharged by infiltration of surface water (116,000 acre-ft/yr) through streams, floodwater, and irrigation, with a lesser amount as groundwater inflow from uplands and minimal recharge beneath Malheur and Harney Lakes. Recharge from streams and floodwater (natural and irrigation) was estimated using a balance of measured and estimated surface-water inflow to and outflow from lowland areas including streamflow, springflow, and ET where a portion of surface-water inflow to lowland areas is comprised of upland discharge to streams and springs. Groundwater ET (119,000 acre-ft/yr) is the predominant natural discharge mechanism in lowland areas and was estimated as the mean from two remote-sensing based approaches incorporating groundwater ET measurements from other similar basins and 23 years (1987–2015) of Landsat imagery. Discharge of lowland groundwater into Malheur and Harney Lakes is about 700 acre-ft/yr and is represented in groundwater ET estimates. The remaining natural groundwater discharge from lowland areas issues from Sodhouse Spring (8,900 acre-ft/yr) and as groundwater flow to the Malheur River Basin through Virginia Valley (3,100 acre-ft/yr). The relatively large amount of groundwater discharged to springs in Warm Springs Valley (25,000 acre-ft/yr) is accounted for in groundwater ET estimates. Natural groundwater discharge in lowland areas of the Harney Basin has remained relatively constant during the last 80 years based on comparisons with estimates north of Malheur Lake and west of Harney Lake published in the 1930s. Annual net amount of groundwater pumped (pumpage) from the Harney Basin during 2017–18 averaged 144,000 acre-ft. The net value is the difference between pumpage (about 152,000 acre-ft/yr) and reinfiltration of groundwater pumped for irrigation and non-irrigation purposes (about 8,000 acre-ft/yr). Net pumpage was estimated in concurrent studies that compiled groundwater-use data and coupled reported groundwater pumpage data from wells with remote-sensing-based ET estimates from groundwater-irrigated fields. Total pumpage for irrigation has increased from about 54,000 acre-ft/yr during 1991–92 to 145,000 acre-ft/yr during 2017–18. Presently, pumpage is greatest in the lowland region north of Malheur Lake (81,000 acre-ft/yr), with lesser amounts to the north and northwest of Harney Lake (41,000 acre-ft/yr) and to the south and east of Malheur Lake (22,000 acre-ft/yr). During this study, mean annual lowland groundwater discharge (including pumpage) exceeded mean annual recharge, indicating that the lowland hydrologic budget is out of balance. Net groundwater pumpage during 2017–18 is similar to groundwater discharge from all other sources in the lowlands and is four times the imbalance between non-pumping lowland recharge and discharge (34,000 acre-ft/yr). Declining groundwater levels at depth across many parts of the Harney Basin lowlands indicate that pumpage is depleting aquifer storage and is likely capturing a small amount of natural groundwater discharge to springs and ET in some lowland areas. If pumping continues, aquifer storage depletion will continue until the capture rate of natural discharge to springs and ET is equal to the pumping rate. If groundwater development occurs in upland areas and reduces either the streamflow or groundwater inflow to lowland areas, the deficit in the lowland water budget will increase.

Oregon↗

Hydrogeology, simulation of regional ground-water flow, and saltwater intrusion, Potomac-Raritan-Magothy Aquifer System, Northern Coastal Plain of New Jersey

The Potomac-Raritan-Magothy aquifer system in Middlesex and Monmouth Counties in the northern Coastal Plain of New Jersey consists primarily of unconsolidated Cretaceous sediments, which are divided into the upper and middle aquifers and confining units. These units, which strike northeastsouthwest along the Fall Line, dip and thicken to the southeast. The upper aquifer consists primarily of the Old Bridge Sand Member of the Magothy Formation, which is composed of coarse-grained sands, localized thin clay beds, and younger surficial sands and gravels in and near the outcrop. Transmissivity ranges from 1,760 to 19,400 ft 2 /d (feet squared per day) and tends to be higher in updip areas. Estimated withdrawals from the upper aquifer in the northern Coastal Plain were approximately 42 Mgal/d (million gallons per day) in 1986. Cones of depression whose centers range from 36 to 42 ft (feet) below sea level have developed as a result of these withdrawals. The upper aquifer is confined throughout most of the northern New Jersey Coastal Plain by clays and silts of the Cretaceous Woodbury Clay and Merchantville Formation and younger sediments of the Magothy Formation. This confining unit generally is greater than 200 ft thick. The simulated vertical hydraulic conductivity for the confining unit ranges from 8.4 x 10 -5 to 5.6 x 10 -3 feet per day; interpreted vertical hydraulic conductivities generally are lower except in southwestern Middlesex County, where the vertical hydraulic conductivities of the confining unit are higher. The middle aquifer consists primarily of the Farrington Sand Member of the Cretaceous Raritan Formation and surficial Holocene and Miocene sands and gravels in its outcrop area. It also can include the uppermost sands of the Cretaceous Potomac Group in parts of Monmouth County. The middle aquifer is composed of fine to coarse sand that contains some lignite and pyrite, and, locally, some clay beds. It pinches out in the northern part of Sayreville Township, near Raritan River. The transmissivity of the aquifer ranges from 2,140 to 13,800 ft 2/ d and tends to decrease in the northern part of the northern Coastal Plain of New Jersey where the aquifer thins. A poorly permeable confining unit composed mostly of clays and silts of the Woodbridge Clay Member of the Raritan Formation overlies the aquifer in most of this area. The confining unit generally is greater than 100 ft thick, although it thins and is sandy in the southwestern part of Middlesex County, where a good hydraulic connection exists between the middle and upper aquifers. Estimated withdrawals from the middle aquifer in the northern Coastal Plain were about 22 Mgal/d in 1986. These withdrawals have caused cones of depression whose centers range from 77 to 93 ft below sea level. A finite-difference, quasi-three-dimensional ground-water flow model was developed to simulate ground-water flow in the aquifer system. The confined and unconfined areas of the upper and middle aquifers were modeled as separate layers. The model was calibrated primarily by adjusting vertical hydraulic conductivity in the confining units and horizontal hydraulic conductivity in the aquifers, then matching simulated and measured groundwater levels for the period 1896-1986 and simulated and interpreted potentiometric surfaces under predevelopment conditions and in 1984. For the predevelopment period, the total flow into and out of the upper and middle aquifers is 35 and 21 Mgal/d, respectively. Recharge to the aquifer system is from direct recharge in the unconfined areas and from vertical leakage through overlying confining units. The main recharge areas are the topographically high areas in southwestern Middlesex County for both aquifers, in the eastern Sayreville area for the upper aquifer, and north of the Raritan River for the middle aquifer. Most ground water discharges to low-lying regional surface-water drains (streams), which flow into the South River. For 1984 transient conditions, the total ground-water flow into and out of the upper and middle aquifers is 61 and 34 Mgal/d, respectively. The largest amount of recharge is from direct recharge in the unconfined areas, but some recharge also is derived from vertical leakage through the Merchantville-Woodbury confining unit, captured ground-water discharge to streams, and induced inflow at artificial-recharge facilities. Regional flow is from recharge areas toward major cones of depression. Sensitivity analysis showed that the model was useful for representing flow in the system, especially in the confined-aquifer areas. Model representation of lateral and vertical boundary conditions was judged acceptable. Simulation results were less sensitive to changes in aquifer properties in the unconfined areas of the aquifers and to changes in storage in the confining units. Sensitivity analysis and calibration of hydraulic parameters and conditions showed that the distribution of hydraulic head was sensitive to changes in horizontal hydraulic conductivity in the aquifers, vertical hydraulic conductivity in the confining units, magnitudes of ground-water withdrawals, and initial hydraulic head in aquifer outcrop areas. Two scenarios were simulated to determine the effects of ground-water withdrawals from 1986 through 2019. For the scenario in which ground-water withdrawals increase to about 69 Mgal/d in the upper aquifer and 37 Mgal/d in the middle aquifer, centers of cones of depression are as deep as 100 ft below sea level in the upper aquifer and 170 ft below sea level in the middle aquifer. For this scenario, most of the additional water comes from captured surface-water discharge, induced cross-formational flow from overlying aquifers, and increases in induced flow from artificial-recharge areas. Induced flow from Raritan Bay also increases. For the scenario in which ground water withdrawals are reduced to 42.5 Mgal/d in the upper aquifer and 15 Mgal/d in the middle aquifer, water levels recover to above sea level nearly everywhere. In each aquifer, ground-water discharge to streams increases and induced flow through the confining units and from the overlying sediments decreases, and discharge of ground water to Raritan Bay in the upper aquifer exceeds the induced recharge from Raritan Bay. Reversal of ground-water gradients has caused saltwater intrusion in the two aquifers. Chloride concentrations in water from the upper aquifer in Keyport and Union Beach Boroughs were as high as 2,100 mg/L (milligrams per liter) in 1986. The intrusion has not increased significantly since well fields in the area were closed in the late 1970's. Elevated chloride concentrations also were measured in Keanesburg Borough in 1986. In both of these areas, saltwater has entered the upper aquifer from the Bay because of movement of the freshwater-saltwater interface in response to increasing ground-water withdrawals. Chloride concentrations in well-water samples from the middle aquifer were as high as 6,000 mg/L in Sayreville Borough in 1987; concentrations in samples from drive-point wells from the same aquifer near the Washington Canal, the main source of saltwater, were as high as 7,100 mg/L. The migration of the saltwater front at about 470 feet per year to the southeast is influenced mainly by a thinning of the middle aquifer, which constrains flow, and by the locations of regional cones of depression caused by groundwater withdrawals.

New Jersey↗

Groundwater and surface-water interaction, water quality, and processes affecting loads of dissolved solids, selenium, and uranium in Fountain Creek, near Pueblo, Colorado, 2012–2014

In 2012, the U.S. Geological Survey, in cooperation with the Arkansas River Basin Regional Resource Planning Group, initiated a study of groundwater and surface-water interaction, water quality, and loading of dissolved solids, selenium, and uranium to Fountain Creek near Pueblo, Colorado, to improve understanding of sources and processes affecting loading of these constituents to streams in the Arkansas River Basin. Fourteen monitoring wells were installed in a series of three transects across Fountain Creek near Pueblo, and temporary streamgages were established at each transect to facilitate data collection for the study. Groundwater and surface-water interaction was characterized by using hydrogeologic mapping, groundwater and stream-surface levels, groundwater and stream temperatures, vertical hydraulic-head gradients and ratios of oxygen and hydrogen isotopes in the hyporheic zone, and streamflow mass-balance measurements. Water quality was characterized by collecting periodic samples from groundwater, surface water, and the hyporheic zone for analysis of dissolved solids, selenium, uranium, and other selected constituents and by evaluating the oxidation-reduction condition for each groundwater sample under different hydrologic conditions throughout the study period. Groundwater loads to Fountain Creek and in-stream loads were computed for the study area, and processes affecting loads of dissolved solids, selenium, and uranium were evaluated on the basis of geology, geochemical conditions, land and water use, and evapoconcentration. During the study period, the groundwater-flow system generally contributed flow to Fountain Creek and its hyporheic zone (as a single system) except for the reach between the north and middle transects. However, the direction of flow between the stream, the hyporheic zone, and the near-stream aquifer was variable in response to streamflow and stage. During periods of low streamflow, Fountain Creek generally gained flow from groundwater. However, during periods of high streamflow, the hydraulic gradient between groundwater and the stream temporarily reversed, causing the stream to lose flow to groundwater. Concentrations of dissolved solids, selenium, and uranium in groundwater generally had greater spatial variability than surface water or hyporheic-zone samples, and constituent concentrations in groundwater generally were greater than in surface water. Constituent concentrations in the hyporheic zone typically were similar to or intermediate between concentrations in groundwater and surface water. Concentrations of dissolved solids, selenium, uranium, and other constituents in groundwater samples collected from wells located on the east side of the north monitoring well transect were substantially greater than for other groundwater, surface-water, and hyporheic-zone samples. With one exception, groundwater samples collected from wells on the east side of the north transect exhibited oxic to mixed (oxic-anoxic) conditions, whereas most other groundwater samples exhibited anoxic to suboxic conditions. Concentrations of dissolved solids, selenium, and uranium in surface water generally increased in a downstream direction along Fountain Creek from the north transect to the south transect and exhibited an inverse relation to streamflow with highest concentration occurring during periods of low streamflow and lowest concentrations occurring during periods of high streamflow. Groundwater loads of dissolved solids, selenium, and uranium to Fountain Creek were small because of the small amount of groundwater flowing to the stream under typical low-streamflow conditions. In-stream loads of dissolved solids, selenium, and uranium in Fountain Creek varied by date, primarily in relation to streamflow at each transect and were much larger than computed constituent loads from groundwater. In-stream loads generally decreased with decreases in streamflow and increased as streamflow increased. In-stream loads of dissolved solids and selenium increased between the north and middle transects but generally decreased between the middle and south transects. By contrast, uranium loads generally decreased between the north and middle transects but increased between the middle and south transects. In-stream load differences between transects appear primarily to be related to differences in streamflow. However, because groundwater typically flows to Fountain Creek under low-flow conditions, and groundwater has greater concentrations of dissolved solids, selenium, and uranium than surface water in Fountain Creek, increases in loads between transects likely are affected by inflow of groundwater to the stream, which can account for a substantial proportion of the in-stream load difference between transects. When loads decreased between transects, the primary cause likely was decreased streamflow as a result of losses to groundwater and flow through the hyporheic zone. However, localized groundwater inflow likely attenuated the magnitude by which the in-stream loads decreased. The combination of localized soluble geologic sources and oxic conditions likely is the primary reason for the occurrence of high concentrations of dissolved solids, selenium, and uranium in groundwater on the east side of the north monitoring well transect. To evaluate conditions potentially responsible for differences in water quality and redox conditions, physical characteristics such as depth to water, saturated thickness, screen depth below the water table, screen height above bedrock, and aquifer hydraulic conductivity were compared by using Wilcoxon rank-sum tests. Results indicated no significant difference between depth to water, screen height above bedrock, and hydraulic conductivity for groundwater samples collected from wells on the east side of the north transect and groundwater samples from all other wells. However, saturated thickness and screen depth below the water table both were significantly smaller for groundwater samples collected from wells on the east side of the north transect than for groundwater samples from other wells, indicating that these characteristics might be related to the elevated constituent concentrations found at that location. Similarly, saturated thickness and screen depth below the water table were significantly smaller for groundwater samples under oxic or mixed (oxic-anoxic) conditions than for those under anoxic to suboxic conditions. The greater constituent concentrations at wells on the east side of the north transect also could, in part, be related to groundwater discharge from an unnamed alluvial drainage located directly upgradient from that location. Although the quantity and quality of water discharging from the drainage is not known, the drainage appears to collect water from a residential area located upgradient to the east of the wells, and groundwater could become concentrated in nitrate and other dissolved constituents before flowing through the drainage. High levels of nitrate, whether from anthropogenic or natural geologic sources, could promote more soluble forms of selenium and other constituents by affecting the redox condition of groundwater. Whether oxic conditions at wells on the east side of the north transect are the result of physical characteristics or of groundwater inflow from the alluvial drainage, the oxic conditions appear to cause increased dissolution of minerals from the shallow shale bedrock at that location. Because ratios of hydrogen and oxygen isotopes indicate evaporation likely has not had a substantial effect on groundwater, constituent concentrations at that location likely are not the result of evapoconcentration.

Colorado↗

Review: “Jacob’s Zoo”— How using Jacob’s method for aquifer testing leads to more intuitive understanding of aquifer characteristics

The interpretation of aquifer responses to pumping tests is an important tool for assessing aquifer geometry and properties, which are critical in the assessment of water resources or in environmental remediation. However, the responses of aquifers, measured by time-drawdown relationships in monitoring wells, are nonunique solutions that are affected by many factors. Jacob’s Zoo is a collection of graphical interpretations that allows students and practitioners to develop an intuitive feel for how natural hydrogeological systems work, and develop a set of skills that provide a better understanding of aquifer properties far beyond interpretation of pumping tests. Jacob’s Zoo, based on the work of Jacob (1950), fosters a deeper understanding, although few practitioners realize the full utility of the method. Jacob CE (1950) Flow of groundwater, In: Rouse H (ed) Engineering Hydraulics, Wiley, New York. P 321–386.

Hydrogeology Journal↗

The quality of our Nation's waters: Water quality in the Mississippi embayment-Texas coastal uplands aquifer system and Mississippi River Valley alluvial aquifer, south-central United States, 1994-2008

About 8 million people rely on groundwater from the Mississippi embayment—Texas coastal uplands aquifer system for drinking water. The Mississippi River Valley alluvial aquifer also provides drinking water for domestic use in rural areas but is of primary importance to the region as a source of water for irrigation. Irrigation withdrawals from this aquifer are among the largest in the Nation and play a key role in the economy of the area, where annual crop sales total more than $7 billion. The reliance of the region on both aquifers for drinking water and irrigation highlights the importance of long-term management to sustain the availability and quality of these resources.

Alabama, Arkansas, Kentucky, Louisiana, Mississipp↗

Groundwater-level, groundwater-temperature, and barometric-pressure data, July 2017 to February 2018, Hālawa Area, Oʻahu, Hawaiʻi

The Red Hill Bulk Fuel Storage Facility, operated by the U.S. Navy and located in the Hālawa area, Oʻahu, Hawaiʻi, includes 20 underground storage tanks that can hold a total of 250 million gallons of fuel. In January 2014, the U.S. Navy notified the Hawaiʻi Department of Health and U.S. Environmental Protection Agency of release of an estimated 27,000 gallons of fuel from the Red Hill Bulk Fuel Storage Facility. In response to past and potential future fuel releases, data are needed to evaluate groundwater flow in the surrounding area. During July 2017–February 2018, the U.S. Geological Survey collected groundwater-level data at 24 sites near the Red Hill Bulk Fuel Storage Facility. At 14 of the 24 sites, groundwater-temperature data were also collected, and at 6 of the 24 sites, barometric-pressure data were collected. During the data-collection period, a regional aquifer test was conducted in coordination with the operators of three production wells in the area. The recorded water-level changes in response to planned withdrawal changes during December 2017–February 2018 can be used to improve understanding of the groundwater-flow system. The scope of this report is limited to a non-interpretive presentation of the data and a brief discussion of the factors affecting the water-level data.

Hawai'i↗

Water table in the surficial aquifer and potentiometric surface of the Floridan Aquifer in selected well fields, west-central Florida, September 1977

The water table in the surficial aquifer and the potentiometric surface of the Floridan aquifer in a 1,200 square-mile area in west-central Florida are mapped semiannually by the U.S. Geological Survey. Maps are prepared on the basis of water levels measured in wells each May to coincide with seasonal low levels; and in September, when levels are high. The mapped area for this report contains nine producing wells fields which supplied 76.2 million gallons on September 21, 1977, to municipalities in the Tampa Bay area. The effect of localized withdrawal of groundwater from the Floridan aquifer is shown on the maps as cones of depression in both the potentiometric surface and water table. In September, water levels in observation wells in the Floridan aquifer were above those measured in May. These increases ranged from about 1 foot at the Eldridge-Wilde well field to about 25 feet at the Sun City well field. (Woodard-USGS)

Open-File Report↗