Geology topics

Groundwater & aquifers

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

Follow this topic with RSS

At least 1,117 records · Page 62Linked to original sources

USGS2018

Santa Barbara and Foothill groundwater basins Geohydrology and optimal water resources management—Developed using density dependent solute transport and optimization models

Groundwater has been a part of the city of Santa Barbara’s water-supply portfolio since the 1800s; however, since the 1960s, the majority of the city’s water has come from local surface water, and the remainder has come from groundwater, State Water Project, recycled water, increased water conservation, and as needed, seawater desalination. Although groundwater from the Santa Barbara and Foothill groundwater basins only accounts for a small percentage of the long-term supply, it is an important source of supplemental water during times of surface-water shortages. During the late 1980s and early 1990s, production wells extracted additional groundwater to compensate for drought related water-delivery shortfalls from other sources; in response, water levels declined substantially in the Santa Barbara and Foothill groundwater basins (below sea level in the Santa Barbara groundwater basin). In coastal basins that have groundwater extraction near shore, seawater intrusion is often a problem. Seawater intrusion in the Santa Barbara groundwater basin is thought to be more limited than in other coastal basins because of an offshore fault that acts as a partial barrier to groundwater flow. During the late 1980s and early 1990s, seawater intrusion was observed in the Santa Barbara groundwater basin, as indicated by increased chloride concentrations at several monitoring wells that ranged from 200 ft to 1,300 ft from the ocean and as close as 2,900 ft to the nearest pumping well. This demonstrated that seawater can intrude into the Santa Barbara groundwater basin when groundwater levels fall below sea level near the coast. The city of Santa Barbara is interested in developing a better understanding of the sustainability of its groundwater supplies. In 2014, California adopted historic legislation to manage its groundwater: the Sustainable Groundwater Management Act (SGMA). The SGMA requires the development and implementation of “Groundwater Sustainability Plans” in 127 priority groundwater basins; although Santa Barbara was not a designated priority basin, the city is taking steps to achieve sustainability. Sustainability was defined in the SGMA in terms of avoiding undesirable results: significant and unreasonable groundwater-level declines, reduction in groundwater storage, seawater intrusion, water-quality degradation, land subsidence, and surface-water depletion. In this project, a cooperative study between the U.S. Geological Survey (USGS) and the city of Santa Barbara, sustainable yield is defined as the volume of groundwater that can be pumped from storage without causing water-level drawdowns and the associated increases in seawater intrusion (as indicated by increases in measured chloride concentrations) at selected wells. In order to estimate the sustainability of Santa Barbara’s groundwater basins, a three-dimensional density-dependent groundwater-flow and solute-transport model (the Santa Barbara Flow and Transport Model, or SBFTM) was developed on the basis of an existing groundwater-flow model. To simulate seawater intrusion to the Santa Barbara Basin under various management strategies, the SBFTM uses the USGS code SEAWAT to simulate salinity transport and variable-density flow. The completed SBFTM was coupled with a management optimization tool, in this case a multi-objective evolutionary algorithm, to determine optimal pumping strategies that maximize the sustainable yield and at the same time satisfy user-defined drawdown and chloride-concentration constraints. As part of this study, a three-dimensional hydrogeologic framework model was developed to quantify the extent and hydrogeologic characteristics of the Santa Barbara and Foothill groundwater basins and to help define the discretization and hydraulic properties used in the SBFTM. The development of the hydrogeologic framework model required the collection and reconciliation of geologic and geophysical data from existing maps, reports, and databases, along with geologic and hydrologic data from recently drilled wells. These data were integrated into a three-dimensional hydrogeologic framework model that defines the stratigraphy and geometry of the aquifer zones and the major geologic structures in the basin. The hydrogeologic framework model also quantifies the variation in sediment grain size within each aquifer zone as the percentage of coarse-grained sediment. Previous studies indicated that there are two principal water-producing zones in the Santa Barbara groundwater basin, the upper and lower producing zones; an additional thin, productive zone was identified as part of this study. This “middle producing zone” is not as areally extensive as the upper and lower producing zones and only exists in the coastal part of Storage Unit I. These producing zones are bounded at depth by less productive shallow, middle, and deep zones. Two versions of the SBFTM were constructed: an initial-condition model and a modern transient model. The initial-condition model is a long-term transient model that simulates flow and solute-transport conditions during a period with limited anthropogenic influences preceeding the modern transient model. The simulation-transient model simulates flow and transport conditions from 1929 through 2013; however, because of data availability, the focus of the model calibration was 1972–2013. The SBFTM was calibrated to measured groundwater levels and drawdown, as well as measured chloride concentrations and change in concentrations, using a combination of automated and trial-and-error parameter-estimation techniques. A sensitivity analysis indicated that, in general, the SBFTM was most sensitive to recharge- and pumping-distribution parameters, specifically those controlling the amount of small-catchment recharge and the distribution of water extraction by hydrogeologic layer for production wells. The model was also sensitive to parameters controlling stream-recharge rates, horizontal and vertical hydraulic conductivity, and porosity. From 1929 to 1971, most of the water entering the area represented by the SBFTM was from creek and small-catchment recharge, and the majority of water leaving the SBFTM area was from pumping, discharge to creeks, and drains. In addition, about 37 percent of the total pumpage came from a net reduction in groundwater storage. From 1972 to 2013, the amount of water entering and leaving the SBFTM was fairly similar as that from 1929 to 1971, except the reduction in pumpage added about 17,000 acre-ft of water to storage. During this later period, there were also times of storage loss. For example, during July 1990, a month when approximately 705 acre-ft of groundwater was pumped in the study area, the pumpage was much greater than all sources of recharge combined, and about 382 acre-ft of water was removed from groundwater storage. Simulated hydraulic heads replicated the observed data to an acceptable matching of the measured water-level, flow direction, and vertical gradients. Simulated hydrographs for selected wells were in good agreement with the measured data, with an average residual of -2.7 ft and a standard deviation of 14.5 ft, indicating that the simulated heads, on average, underestimated the observed water levels. An examination of the model fit indicated that most of the discrepancies were lower simulated heads at wells proximal to production well sites. The simulated chloride concentrations reasonably matched the rising limbs of the measured breakthrough curves in terms of timing and magnitude; however, the simulation overestimated the chloride concentrations on the falling limbs. The overestimation of low chloride concentrations was attributed to the model overestimating the advance of the chloride front during periods of heavy pumping and underestimating the retreat of the chloride front during periods of low pumping. These simulation errors would result in a conservative response by local water managers to seawater intrusion. The SBFTM was used to develop a collection of predictive simulations optimized to produce pumping schedules that maximize yield, subject to a set of constraints and competing objectives. The simulations were grouped as scenarios that differed in their time horizon, initial conditions for groundwater levels and chloride concentrations, as well as precipitation, which was incorporated into the model through simulated recharge. Overall, five scenarios were developed in a multi-objective framework to obtain optimal pumping rates for all of the wells managed by the city, while minimizing excessive drawdown and seawater intrusion. For the current study, complexities in the simulation model and the optimization formulation required additional considerations. Incorporating the solute-transport equations to simulate chloride transport added a highly nonlinear process that is solved iteratively in each time step of the groundwater-flow model. These nonlinearities, coupled with the highly refined grid in the current model, creates challenges for many traditional optimization methods. Therefore, an optimization method was needed that could address nonlinear relationships as well as a very large problem size. Lastly, the optimization problem was reformulated to include multiple objectives without requiring convergence to a single solution. This approach, guided by the city’s objectives, allowed the maximum extraction of information from the complex simulation. Borg, a multi-objective evolutionary algorithm, was chosen as the optimization algorithm for this study for several reasons: (1) it is very computationally efficient; (2) it can run in parallel; (3) it requires little user input; and (4) it can solve for multiple competing objectives. The first three points allow the algorithm to proceed toward the optimal solutions at the fastest possible rate. The fourth point is advantageous for large, complex optimization problems because it is difficult to formulate the optimization problem in a way that produces only one optimal solution. The problem formulation consisted of four competing objectives and a constraint set in accordance with the main concerns of the city. The objectives were maximizing total pumpage, minimizing seawater intrusion, minimizing total drawdown in production wells, and minimizing the maximum drawdown. The constraints were pump capacity, meeting drinking-water standards for chloride, maintaining a specified minimum flowrate to a groundwater treatment plant, and maintaining minimum water levels in pumping wells. The decision variables either were quarterly pumpage by well or total pumpage by basin. Five optimization scenarios were developed that allow the decision makers to evaluate a range of optimal solutions for a variety of water levels and chloride concentrations as well as potential future climatic conditions. Three scenarios (1, 2, and 5) were multi-objective optimization formulations that allowed for variations in management preferences and climatic conditions. The other two scenarios (3 and 4) were designed to examine the optimization results to answer specific questions. Scenario 1 described the best-case sustainable yield assuming a “full” basin (that is, high initial water levels) and typical climate conditions for 10 years. Scenario 2 also started with a “full” basin; however, this was followed by a 10-year drought. Scenario 3 determined if an “empty” basin (that is, low initial water levels) would recover to full conditions (1998 conditions) given climate assumptions and optimal pumping schedules from scenarios 1 and 2. Scenario 4 was designed to produce decision rules that can be used by water managers to help choose an optimal pumping schedule based on measured water-level or chloride data. Scenario 5 identified future pumping schedules based on short-term climate variations during a 2-year management horizon. The results from scenarios 1 and 2 described the differences in maximum pumpage in the basin under typical and dry long-term climate projections, respectively. The scenario 1 results indicated the maximum 10-year pumpage of the basin was about 31,300 acre-ft under typical conditions and controlling simulated seawater intrusion and drawdowns. For scenario 2, less recharge over the 10-year dry climate produced a maximum pumpage estimate of 30,000 acre-ft to control seawater intrusion and drawdowns. The larger pumpage for scenario 1 resulted in more seawater intrusion, but less total drawdown, compared to that of scenario 2. Results for scenarios 3 and 4 showed the basin’s response to management actions combined with climate projections. Both scenarios used the optimal pumping schedules and the 10-year climates from scenarios 1 and 2. The scenario 3 results showed that under minimal pumping, the basin did not fully recover to 1998 water levels within 10 years under either climate scenario. The relatively larger recharge from the typical climate resulted in less drawdown at coastal monitoring wells after the 10-year recovery period than that from the dry climate. The location of the seawater intrusion front was not appreciably different between the scenarios, however. Scenario 4 used the optimal results from scenarios 1 and 2 to produce decision-rule curves that illustrated the pumpage for each basin, given measured levels of chloride concentration or drawdown. This allowed the use of additional measurements at monitoring wells to assess future management decisions on the basis of the sensitivity of observations of drawdown and seawater intrusion to various pumping rates. Scenario 5 allowed managers to investigate the effects of short-term climate variations on optimal pumping schedules. Three specific 2-year simulations were optimized: typical-to-dry (scenario 5A), dry-to-typical (scenario 5B), and dry-to-dry (scenario 5C). The most noteable result from scenario 5 was the overall reduction in optimal pumpage for most schedules in scenario 5C, when the climate is simulated as dry-to-dry. There are also many optimal pumping schedules that produced an overall increase in waterlevels over the two-year simulation period, regardless of climatic condition. Similar to scenario 2, the scenario 5C results represents conservative yield estimates under a minimal-precipitation climatic condition.

California↗

USGS2022

Mapped predictions of manganese and arsenic in an alluvial aquifer using boosted regression trees

Manganese (Mn) concentrations and the probability of arsenic (As) exceeding the drinking-water standard of 10 μg/L were predicted in the Mississippi River Valley alluvial aquifer (MRVA) using boosted regression trees (BRT). BRT, a type of ensemble-tree machine-learning model, were created using predictor variables that affect Mn and As distribution in groundwater. These variables included iron (Fe) concentrations and specific conductance predicted from previously developed BRT models, groundwater flux and age estimates from MODFLOW, and hydrologic characteristics. The models also included results from the first airborne geophysical survey conducted in the United States to target an entire aquifer system. Predictions of high Mn and As occurred where Fe was high. Predicted high Mn concentrations were correlated with fraction of young groundwater (less than 65 years) computed from MODFLOW results. High probabilities of As exceedance were predicted where groundwater was relatively old and airborne electromagnetic resistivity was high, typically proximal to streams. Two-variable partial-dependence plots and sensitivity analysis were used to provide insight into the factors controlling Mn and As distribution in groundwater. The maps of predicted Mn concentrations and As exceedance probabilities can be used to identify areas where these constituents may be high, and that could be targeted for further study. This paper shows that incorporation of a selected set of process-informed data, such as MODFLOW results and airborne geophysics, into a machine-learning model improves model interpretability. Incorporation of process-rich information into machine-learning models will likely be useful for addressing a wide range of problems of interest to groundwater hydrologists.

Arkansas, Louisiana, Mississippi, Tennessee↗

USGS2016

Budgets and chemical characterization of groundwater for the Diamond Valley flow system, central Nevada, 2011–12

The Diamond Valley flow system consists of six hydraulically connected hydrographic areas in central Nevada. The general down-gradient order of the areas are southern and northern Monitor Valleys, Antelope Valley, Kobeh Valley, Stevens Basin, and Diamond Valley. Groundwater flow in the Diamond Valley flow system terminates at a large playa in the northern part of Diamond Valley. Concerns relating to continued water-resources development of the flow system resulted in a phased hydrologic investigation that began in 2005 by the U.S. Geological Survey in cooperation with Eureka County. This report presents the culmination of the phased investigation to increase understanding of the groundwater resources of the basin-fill aquifers in the Diamond Valley flow system through evaluations of groundwater chemistry and budgets. Groundwater chemistry was characterized using major ions and stable isotopes from groundwater and precipitation samples. Groundwater budgets accounted for all inflows, outflows, and changes in storage, and were developed for pre-development (pre-1950) and recent (average annual 2011–12) conditions. Major budget components include groundwater discharge by evapotranspiration and groundwater withdrawals; groundwater recharge by precipitation, and interbasin flow; and storage change. Groundwater in the basin-fill aquifer of the Diamond Valley flow system was mostly a calcium or sodium bicarbonate water type and generally within acceptable drinking-water standards. The general water type was similar among the individual hydrographic areas. Stable isotopes of oxygen-18 and deuterium from precipitation varied seasonally, such that enrichment from evaporation was greater during warmer months than cooler months. The isotopic signature of shallow groundwater was similar to cool season precipitation, indicating recharge was relatively recent (similar to recent climatic conditions) and was derived from cool season precipitation. Site-scale groundwater evapotranspiration was estimated from eddy-covariance and micrometeorological measurements collected at four sites and ranged from 0.15 feet per year in sparse, undisturbed shrubland to 1.13 feet per year in a grassland meadow. Vegetation indices calculated from satellite imagery and field mapping were used to define three evapotranspiration units (shrubland, grassland, and playa) and to extrapolate site-scale groundwater evapotranspiration rates to basin-scale estimates. Annual pre-development groundwater evapotranspiration for individual hydrographic areas ranged from 2,900 acre-feet per year (acre-ft/yr) in northern Monitor Valley to 35,000 acre-ft/yr in Diamond Valley. Total groundwater evapotranspiration from the Diamond Valley flow system under pre-development conditions was about 70,000 acre-ft/yr. Areas of irrigated land in the Diamond Valley flow system increased from less than 5,000 acres in the early 1960s to more than 25,000 acres in 2012 and are mostly for growing alfalfa in southern Diamond Valley. Annual (2011–12) net groundwater withdrawals for irrigation, assumed to be the volume of groundwater consumed by crops and pastureland, ranged from about 420 acre-ft/yr in Antelope Valley to 67,000 acre-ft/yr in Diamond Valley. Total net groundwater withdrawals for irrigation in the Diamond Valley flow system were about 69,000 acre-ft/yr (2011–12). Groundwater recharge, the largest inflow component to the Diamond Valley flow system, was determined as the sum of groundwater evapotranspiration and net subsurface outflow (subsurface outflow minus subsurface inflow). Annual groundwater recharge estimates ranged from 200 acre-ft/yr in Stevens Basin to 35,000 acre-ft/yr in Diamond Valley. Subsurface flow between hydrographic basins was evaluated using estimated transmissivity, groundwater-flow sections derived from remotely sensed imagery, and hydraulic gradients determined from 2012 water-level data. Subsurface outflow ranged from 0 acre-ft/yr for Diamond Valley to 3,400 acre-ft/yr for northern Monitor Valley into western Kobeh Valley. Subsurface inflow ranged from 0 acre-ft/yr for southern Monitor Valley to 4,200 acre-ft/yr for Kobeh Valley from northern Monitor and Antelope Valleys. The pre-development, steady state, groundwater budget for the Diamond Valley flow system was estimated at about 70,000 acre-ft/yr of inflow and outflow. During years 2011–12, inflow components of groundwater recharge from precipitation and subsurface inflow from adjacent basins totaled 70,000 acre-ft/yr for the DVFS, whereas outflow components included 64,000 acre-ft/yr of groundwater evapotranspiration and 69,000 acre-ft/yr of net groundwater withdrawals, or net pumpage. Spring discharge in northern Diamond Valley declined about 6,000 acre-ft/yr between pre-development time and years 2011–12. Assuming net groundwater withdrawals minus spring flow decline is equivalent to the storage change, the 2011–12 summation of inflow and storage change was balanced with outflow at about 133,000 acre-ft/yr.

Nevada↗

USGS1988

Ground-water hydraulics, regional flow, and ground-water development of the Floridan aquifer system in Florida and in parts of Georgia, South Carolina, and Alabama

The Floridan aquifer system is one of the major sources of groundwater supplies in the United States. This productive aquifer system underlies all of Florida, southeast Georgia, and small parts of adjoining Alabama and South Carolina, for a total area of about 100,000 square miles. About 3 billion gallons of water per day were withdrawn from the aquifer system in 1980, and in many areas the Floridan is the sole source of freshwater. The Floridan aquifer system is a sequence of hydraulically connected carbonate rocks (principally limestone with some dolomite) ranging in age generally from late Paleocene to early Miocene. The rocks vary in thickness from a featheredge where they crop out to more than 3,500 feet where the aquifer is deeply buried. The aquifer system generally consists of an upper aquifer and a lower aquifer separated by a lesspermeable confining unit of highly variable properties. In parts of north Florida and southwest Georgia, where little permeability contrast exists among the units, the Floridan is effectively one continuous aquifer. The upper and lower aquifers, named the Upper Floridan aquifer and the Lower Floridan aquifer, are defined on the basis of permeability and their boundaries locally do not coincide with those for either time-stratigraphic or rock-stratigraphic units. Overlying much of the Floridan aquifer system are low-permeability clastic rocks. The lithology, thickness, and integrity of these rocks determine the degree of confinement and influence the distribution of natural recharge, discharge, and ground-water flow in the Floridan. The permeability of the Floridan aquifer system is derived from small openings including fossil hashes and solution-widened joints as well as large cavernous openings in karst areas. Diffuse flow predominates where the small openings occur, whereas conduit flow may occur where large cavernous openings are. Transmissivities are highest (greater than 1,000,000 feet squared per day) in the unconfined karst areas of central and northern Florida. Lowest transmissivities (less than 50,000 feet squared per day) occur in panhandle Florida and southernmost Florida where the Upper Floridan aquifer is confined by thick clay sections. The hydraulic properties of the Lower Floridan aquifer are not well known; however, intervals of high transmissivity occur that have been attributed to paleokarst development. Springs, nearly all of which occur in unconfined and semiconfined parts of the Upper Floridan aquifer in Florida are the dominant feature of the Floridan flow system. Before ground-water development, spring flow and point discharge to surface-water bodies were about 88 percent of the estimated 21,500 cubic feet per second total discharge, or about 19,000 cubic feet per second. Diffuse upward leakage, which occurs primarily in confined areas, accounted for the remaining 12 percent or about 2,500 cubic feet per second. Most of the recharge necessary to sustain springflow and aquifer discharge to streams and lakes occurs relatively close to springs and areas of point discharge to surface-water bodies. Recharge to the Upper Floridan is highest, averaging 10-20 inches per year, in unconfined or semiconfined spring areas. The proximity of high recharge to high discharge implies a vigorous and well-developed shallow flow system in the unconfined and semiconfined parts of the Upper Floridan aquifer. Ground-water flow is very sluggish in the parts of the aquifer system that are deeply buried and tightly confined, primarily southeast Georgia and northeast Florida, south Florida, and far-west panhandle Florida. Discharge to springs, streams, and lakes is practically nonexistent in the tightly confined areas and natural discharge occurs almost exclusively by diffuse upward leakage through thick overburden. The regional flow system has not been appreciably altered by groundwater development. However, increasing pumpage that reached 3 billion gallons per day by 1980 has resulted in long-term regional water level decline of more than 10 feet in three broad areas: coastal Georgia, adjacent South Carolina, and northeast Florida; west-central Florida; and panhandle Florida. Saltwater encroachment as a result of pumping has occurred locally in coastal areas. Pumpage from the Upper Floridan aquifer is supplied primarily by reduction of natural discharge and by increased recharge rather than by depletion of aquifer storage. About 20 percent is from reduced discharge to springs, streams, and lakes, about 20 percent is from reduced upward leakage, and about 60 percent is from increased recharge. Compared to predevelopment conditions, discharge to springs, streams, and lakes is reduced by less than 5 percent, upward leakage is reduced by about 30 percent, and recharge is increased by about 12 percent. Total recharge and, therefore, discharge increased from a predevelopment rate of 21,500 cubic feet per second to about 24,100 cubic feet per second by 1980. A considerable area remains of the Floridan aquifer system where large ground-water supplies may be developed. This area is largely inland from the coasts and characterized by high transmissivity and minimal development prior to the early 1980's. The major constraint on future development probably is degradation of water quality rather than water-quantity limitations.

Alabama, Florida, Georgia, South Carolina↗

USGS2015

Revised hydrogeologic framework of the Floridan aquifer system in Florida and parts of Georgia, Alabama, and South Carolina

The hydrogeologic framework for the Floridan aquifer system has been revised throughout its extent in Florida and parts of Georgia, Alabama, and South Carolina. The updated framework generally conforms to the original framework established by the U.S. Geological Survey in the 1980s, except for adjustments made to the internal boundaries of the Upper and Lower Floridan aquifers and the individual higher and contrasting lower permeability zones within these aquifers. The system behaves as one aquifer over much of its extent; although subdivided vertically into two aquifer units, the Upper and Lower Floridan aquifers. In the previous framework, discontinuous numbered middle confining units (MCUI–VII) were used to subdivide the system. In areas where less-permeable rocks do not occur within the middle part of the system, the system was previously considered one aquifer and named the Upper Floridan aquifer. In intervening years, more detailed data have been collected in local areas, resulting in some of the same lithostratigraphic units in the Floridan aquifer system being assigned to the Upper or Lower Floridan aquifer in different parts of the State of Florida. Additionally, some of the numbered middle confining units are found to have hydraulic properties within the same order of magnitude as the aquifers. A new term “composite unit” is introduced for lithostratigraphic units that cannot be defined as either a confining or aquifer unit over their entire extent. This naming convention is a departure from the previous framework, in that stratigraphy is used to consistently subdivide the aquifer system into upper and lower aquifers across the State of Florida. This lithostratigraphic mapping approach does not change the concept of flow within the system. The revised boundaries of the Floridan aquifer system were mapped by considering results from local studies and regional correlations of lithostratigraphic and hydrogeologic units or zones. Additional zones within the aquifers have been incorporated into the framework to allow finer delineation of permeability variations within the aquifer system. These additional zones can be used to progressively divide the system for assessing groundwater and surface-water interaction, saltwater intrusion, and offshore movement of groundwater at greater detail if necessary. The lateral extent of the updip boundary of the Floridan aquifer system is modified from previous work based on newer data and inclusion of parts of the updip clastic facies. The carbonate and clastic facies form a gradational sequence, generally characterized by limestone of successively younger units that extend progressively farther updip. Because of the gradational nature of the carbonate-clastic sequence, some of the updip clastic aquifers have been included in the Floridan aquifer system, the Southeastern Coastal Plain aquifer system, or both. Thus, the revised updip limit includes some of these clastic facies. Additionally, the updip limit of the most productive part of the Floridan aquifer system was revised and indicates the approximate updip limit of the carbonate facies. The extent and altitude of the freshwater-saltwater interface in the aquifer system has been mapped to define the freshwater part of the flow system.

Alabama, Florida, Georgia, South Carolina↗

USGS2002

Groundwater recharge and agricultural contamination

Agriculture has had direct and indirect effects on the rates and compositions of groundwater recharge and aquifer biogeochemistry. Direct effects include dissolution and transport of excess quantities of fertilizers and associated materials and hydrologic alterations related to irrigation and drainage. Some indirect effects include changes in water–rock reactions in soils and aquifers caused by increased concentrations of dissolved oxidants, protons, and major ions. Agricultural activities have directly or indirectly affected the concentrations of a large number of inorganic chemicals in groundwater, for example NO 3 – , N 2 , Cl, SO 4 2– , H + , P, C, K, Mg, Ca, Sr, Ba, Ra, and As, as well as a wide variety of pesticides and other organic compounds. For reactive contaminants like NO 3 – , a combination of chemical, isotopic, and environmental-tracer analytical approaches might be required to resolve changing inputs from subsequent alterations as causes of concentration gradients in groundwater. Groundwater records derived from multi-component hydrostratigraphic data can be used to quantify recharge rates and residence times of water and dissolved contaminants, document past variations in recharging contaminant loads, and identify natural contaminant-remediation processes. These data indicate that many of the world's surficial aquifers contain transient records of changing agricultural contamination from the last half of the 20th century. The transient agricultural groundwater signal has important implications for long-term trends and spatial heterogeneity in discharge.

Hydrogeology Journal↗

USGS2014

Assessing the magnitude and timing of anthropogenic warming of a shallow aquifer: example from Virginia Beach, USA

Groundwater temperature measurements in a shallow coastal aquifer in Virginia Beach, Virginia, USA, suggest groundwater warming of +4.1 °C relative to deeper geothermal gradients. Observed warming is related to timing and depth of influence of two potential thermal drivers—atmospheric temperature increases and urbanization. Results indicate that up to 30 % of groundwater warming at the water table can be attributed to atmospheric warming while up to 70 % of warming can be attributed to urbanization. Groundwater temperature readings to 30-m depth correlate positively with percentage of impervious cover and negatively with percentage of tree canopy cover; thus, these two land-use metrics explain up to 70 % of warming at the water table. Analytical and numerical modeling results indicate that an average vertical groundwater temperature profile for the study area, constructed from repeat measurement at 11 locations over 15 months, is consistent with the timing of land-use change over the past century in Virginia Beach. The magnitude of human-induced warming at the water table (+4.1 °C) is twice the current seasonal temperature variation, indicating the potential for ecological impacts on wetlands and estuaries receiving groundwater discharge from shallow aquifers.

Virginia↗

USGS2020

Evaluation of soil zone processes and a novel radiocarbon correction approach for groundwater with mixed sources

Estimates of groundwater age based on 14C is often limited by the uncertainty in geochemical processes that alter the 14C concentration measured in water and the composition (δ13C and 14C) of carbon sources needed to appropriately parametrize 14C adjustment models. Estimated ages for samples that contain a mixture of young and old groundwater will be particularly sensitive to model parametrization as relatively small additions of modern 14C from recent recharge can mask the presence and amount of old groundwater. A novel multi-model approach based on inverse geochemical modeling and lumped parameter modeling of age tracers (3H, 3Hetrit, and SF6) was used to better constrain 14C dilution caused by dissolution of carbonates in the unsaturated zone or shallow parts of the Glacial aquifer, which extends over 2000 miles across the northern contiguous United States. Calibration of 14C inverse geochemical models to LPM computed 14C concentrations in modern water indicated that 14C of soil zone and shallow aquifer carbonates were not 14C-dead (0 pmC), as is typically assumed for 14C correction models. 14C of such carbonates was on average about 53 pmC (ranged 0-110 pmC, n = 72). This information was used to correct 14C concentrations for water recharged entirely before 1950 and water that is a mixture of pre- and post-1950 water. The multi-model approach developed here was compared to an analytical 14C-adjustment model (Revised Fontes and Garnier) that assumed solid carbonates were 14C-dead. 14C corrections using the analytical adjustment model tended to over-correct final 14C concentrations by 21 pmC and underestimates mean ages by 40% for groundwater mixtures. In fact, 14C corrections based on analytical model yielded negative ages (14C > 120 pmC) in nearly 36% of mixed samples. This work presents a new approach to constraining 14C corrections and age estimates of mixtures of young and old groundwater. The new method is applied to three well networks distributed across the spatially expansive Glacial aquifer.

Glacier aquifer system↗

USGS2020

Rethinking a groundwater flow system using a multiple-tracer geochemical approach: A case study in Moab-Spanish Valley, Utah

The Glen Canyon Group Aquifer (GCGA) is the sole source of public water supply for the city of Moab, Utah, a domestic and international tourist destination. Population and tourism growth are likely to target the GCGA for future water resources, but our analysis indicates that additional withdrawals would likely be sourced from groundwater storage and not be sustained by recharge. A quantitative estimate of groundwater discharge from the GCGA is problematic because the downgradient aquifer boundary is the Colorado River, and groundwater discharge to the river is very small compared to the river flow. A water budget based on a conceptual model of GCGA discharging into an adjacent alluvial Valley-Fill Aquifer (VFA) was reported by Sumsion (1971) and numerous subsequent studies have repeated and utilized this water budget. The GCGA contains stable isotopes, tritium, 3 He/ 4 He ratios, dissolved solids, and sulfate concentrations that contrast with the VFA, indicating it is instead recharged by local streams rather than from the GCGA. Water-budget calculations, based on: (1) measured spring discharge and streamflow gains, (2) horizontal gradients in VFA groundwater age, and (3) GCGA outcrop area vadose-zone pore waters are all less than previously thought. Using a lumped parameter model and 14 C groundwater ages, we estimate recharge to the deeper GCGA (DGCGA) to be 4.2 ± 2.3 × 10 6 m 3 /yr, which is approximately equal to the measured discharge from wells and springs.

Utah↗

USGS2003

Applicability of tetrazolium salts for the measurement of respiratory activity and viability of groundwater bacteria

A study was undertaken to measure aerobic respiration by indigenous bacteria in a sand and gravel aquifer on western Cape Cod, MA using tetrazolium salts and by direct oxygen consumption using gas chromatography (GC). In groundwater and aquifer slurries, the rate of aerobic respiration calculated from the direct GC assay was more than 600 times greater than that using the tetrazolium salt 2-(4-iodophenyl)-3-(4-nitrophenyl)-5-phenyl tetrazolium chloride (INT). To explain this discrepancy, the toxicity of INT and two additional tetrazolium salts, sodium 3′-[1-(phenylamino)-carbonyl]-3,4-tetrazolium]-bis(4-methoxy-6-nitro) benzenesulfonic acid hydrate (XTT) and 5-cyano-2,3-ditolyl tetrazolium chloride (CTC), to bacterial isolates from the aquifer was investigated. Each of the three tetrazolium salts was observed to be toxic to some of the groundwater isolates at concentrations normally used in electron transport system (ETS) and viability assays. For example, incubation of cells with XTT (3 mM) caused the density of four of the five groundwater strains tested to decline by more than four orders of magnitude. A reasonable percentage (>57%) of cells killed by CTC and INT contained visible formazan crystals (the insoluble, reduced form of the salts) after 4 h of incubation. Thus, many of the cells reduced enough CTC or INT prior to dying to be considered viable by microscopic evaluation. However, one bacterium ( Pseudomonas fluorescens ) that remained viable and culturable in the presence of INT and CTC, did not incorporate formazan crystals into more than a few percent of cells, even after 24 h of incubation. This strain would be considered nonviable based on traditional tetrazolium salt reduction assays. The data show that tetrazolium salt assays are likely to dramatically underestimate total ETS activity in groundwater and, although they may provide a reasonable overall estimate of viable cell numbers in a community of groundwater bacteria, some specific strains may be falsely considered nonviable by this assay due to poor uptake or reduction of the salts.

Massachusetts↗

USGS2010

Groundwater as a nonpoint source of atrazine and deethylatrazine in a river during base flow conditions

Alluvial groundwater adjacent to the main stem river is the principal nonpoint source of atrazine and deethylatrazine in the Cedar River of Iowa after the river has been in base flow conditions for 5 days. Between two sites along a 116-km reach of the Cedar River, tributaries contributed about 25% of the increase in the atrazine and deethylatrazine load, whereas groundwater from the alluvial aquifer contributed at least 75% of the increase in load. Within the study area, tributaries aggregate almost all of the discharge from tile drains, and yet the tributaries still only contribute 25% of the increase in loads in the main stem river. At an unfamned study site adjacent to the Cedar River, the sources of atrazine and deethylatrazine in the alluvial groundwater are bank storage of river water and groundwater recharge from areas distant from the river. Atrazine and deethylatrazine associated with bank storage water will provide larger concentrations to the river during early base flow conditions. After the depletion of bank storage, stable and smaller concentrations of atrazine and deethylatrazine, originating from groundwater recharge, continue to be discharged from the alluvial aquifer to the river; thus these results indicate that alluvial aquifers are an important nonpoint source of atrazine and deethylatrazine in rivers during base flow.

Iowa↗

USGS2022

Understanding the water resources of a mountain-block aquifer: Tucson Mountains, Arizona

Water resources are limited in arid locations such as Tucson Basin. Residential development in the Tucson Mountains to the west of Tucson, Arizona, is limited by groundwater resources. Groundwater samples were collected from fractured bedrock and alluvial aquifers surrounding the Tucson Mountains to assess water quality and recharge history through measurement of stable O, H, and S isotopes; tritium; and 14 C. Most groundwater is a mixture of different ages but is commonly several thousand years old. A few sampling locations indicated a component of water recharged after the above-ground nuclear testing of the mid 1950s, and these sites may represent locations near where the aquifer receives present-day recharge. The Tucson Mountains also host sulfide deposits associated with fractures and replacement zones; these locally contribute to poor-quality groundwater. Projections of future climate predict intensifying drought in southwestern North America. In the study area, a combination of strategies such as rainwater harvesting, exploitation of renewable water, and low groundwater use could be used for sustainable use of the groundwater supply.

Arizona↗

USGS2013

Groundwater geochemical and selected volatile organic compound data, Operable Unit 1, Naval Undersea Warfare Center, Division Keyport, Washington, June and October 2012

Previous investigations indicate that concentrations of chlorinated volatile organic compounds are substantial in groundwater beneath the 9-acre former landfill at Operable Unit 1, Naval Undersea Warfare Center, Division Keyport, Washington. The U.S. Geological Survey has continued to monitor groundwater geochemistry to ensure that conditions remain favorable for contaminant biodegradation as specified in the Record of Decision for the site. This report presents groundwater geochemical and selected chlorinated volatile organic compound data collected at Operable Unit 1 by the U.S. Geological Survey during June and October 2012, in support of long-term monitoring for natural attenuation. Groundwater samples were collected from 13 wells and 9 piezometers, as well as from 10 shallow groundwater passive-diffusion sampling sites in the nearby marsh. Samples from all wells and piezometers were analyzed for oxidation-reduction (redox) sensitive constituents and dissolved gases. Samples from all piezometers also were analyzed for chlorinated volatile organic compounds, as were all samples from the passive-diffusion sampling sites. In 2012, concentrations of redox-sensitive constituents measured at all wells and piezometers were consistent with those measured in previous years, with dissolved oxygen concentrations all at 0.4 milligram per liter or less; little to no detectable nitrate; abundant dissolved manganese, iron, and methane; and commonly detected sulfide. In the upper aquifer of the northern plantation in 2012, chlorinated volatile organic compound (CVOC) concentrations at all piezometers were similar to those measured in previous years, and concentrations of the reductive dechlorination byproducts ethane and ethene were slightly higher or the same as concentrations measured in 2011. In the upper aquifer of the southern plantation, CVOC concentrations measured in piezometers during 2012 continued to be extremely variable as in previous years, and often very high, and reductive dechlorination byproducts were detected in two of the four wells and in all piezometers. Beneath the marsh adjacent to the southern plantation, chloroethene concentrations measured in 2012 continued to vary spatially and temporarily, and also were very high. Additionally, CVOC concentrations measured in samplers deployed in access tubes were about two to four times less than those measured in the two samplers buried nearby, beneath the marsh stream. Total CVOC concentration, at what has been historically the most contaminated passive-diffusion sampler site (S-4), continued an increasing trend. For the intermediate aquifer in 2012, concentrations of reductive dechlorination byproducts ethane and ethene were consistent with those measured in previous years.

Washington↗

USGS2014

Groundwater geochemical and selected volatile organic compound data, Operable Unit 1, Naval Undersea Warfare Center, Division Keyport, Washington, July 2013

Previous investigations indicate that concentrations of chlorinated volatile organic compounds (CVOCs) are substantial in groundwater beneath the 9-acre former landfill at Operable Unit 1, Naval Undersea Warfare Center, Division Keyport, Washington. The U.S. Geological Survey has continued to monitor groundwater geochemistry to ensure that conditions remain favorable for contaminant biodegradation as specified in the Record of Decision for the site. This report presents groundwater geochemical and selected CVOC data collected at Operable Unit 1 by the U.S. Geological Survey during July 9–18, 2013, in support of longterm monitoring for natural attenuation. Groundwater samples were collected from 13 wells and 9 piezometers, as well as from 10 shallow groundwater passive-diffusion sampling sites in the nearby marsh. Samples from all wells and piezometers were analyzed for oxidation-reduction (redox) sensitive constituents and dissolved gases. Samples from all piezometers and four wells also were analyzed for CVOCs, as were all samples from the passive-diffusion sampling sites. In 2013, concentrations of redox-sensitive constituents measured at all wells and piezometers were consistent with those measured in previous years, with dissolved oxygen concentrations at all except an upgradient well 0.2 milligrams per liter or less; little to no detectable nitrate; abundant dissolved manganese, iron, and methane; and commonly detected sulfide. In the upper aquifer of the northern plantation in 2013, CVOC concentrations at all piezometers were similar to those measured in previous years, and concentrations of the reductive dechlorination byproducts ethane and ethene were slightly lower or the same as concentrations measured in 2012. In the upper aquifer of the southern plantation, CVOC concentrations measured in piezometers during 2013 continued to be variable as in previous years, and often very high, and reductive dechlorination byproducts were detected in two of the three wells and in all but one piezometer. Beneath the marsh adjacent to the southern plantation, chloroethene concentrations measured in 2013 continued to vary spatially and temporaly, and also were very high. Total CVOC concentrations, at what have been historically the most contaminated passive-diffusion sampler sites (S-4, S-4B, S-5, and S-5B) remained elevated. For the intermediate aquifer in 2013, concentrations of reductive dechlorination byproducts ethane and ethene and CVOCs were consistent with those measured in previous years.

Washington↗

USGS2015

Groundwater geochemical and selected volatile organic compound data, Operable Unit 1, Naval Undersea Warfare Center, Division Keyport, Washington, June and September 2014

Previous investigations indicate that concentrations of chlorinated volatile organic compounds (CVOCs) are substantial in groundwater beneath the 9-acre former landfill at Operable Unit 1, Naval Undersea Warfare Center, Division Keyport, Washington. The U.S. Geological Survey has continued to monitor groundwater geochemistry to ensure that conditions remain favorable for contaminant biodegradation at the site. This report presents groundwater geochemical and selected CVOC data collected at Operable Unit 1 by the U.S. Geological Survey during June 23–25 and September 4, 2014, in support of long-term monitoring for natural attenuation. Groundwater samples were collected from 13 wells and 9 piezometers, as well as from 10 shallow groundwater passive-diffusion sampling sites in the nearby marsh. Samples from all wells and piezometers were analyzed for oxidation-reduction (redox) sensitive constituents and dissolved gases. Samples from all piezometers and four wells also were analyzed for CVOCs, as were all samples from the passive-diffusion sampling sites. In 2014, concentrations of redox-sensitive constituents measured at all wells and piezometers were consistent with those measured in previous years, with dissolved oxygen concentrations all less than 1 milligram per liter; little to no detectable nitrate; abundant dissolved manganese, iron, and methane; and commonly detected sulfide. In the upper aquifer of the northern plantation in 2014, CVOC concentrations at all piezometers were similar to those measured in previous years, and concentrations of the reductive dechlorination byproducts ethane and ethene were slightly lower or the same as concentrations measured in 2013. In the upper aquifer of the southern plantation, CVOC concentrations measured in piezometers during 2014 continued to be variable as in previous years, often high, and reductive dechlorination byproducts were detected in one of the three wells and in all but two piezometers. Beneath the marsh adjacent to the southern plantation, chloroethene concentrations measured in 2014 continued to vary spatially and temporally, and were high. Trends for total CVOC concentration continued to increase at the historically most contaminated passive‑diffusion sampler sites (S-4, S-4B, and S-5). For the intermediate aquifer in 2014, concentrations of reductive dechlorination byproducts ethane and ethene and CVOCs were consistent with those measured in previous years.

Washington↗

USGS2019

Water-quality and geochemical variability in the Little Arkansas River and Equus Beds aquifer, south-central Kansas, 2001–16

This fact sheet describes water quality and geochemistry of the Little Arkansas River and Equus Beds aquifer during 2001 through 2016 as part of the City of Wichita’s Equus Beds aquifer storage and recovery project in south-central Kansas. The Equus Beds aquifer storage and recovery project was developed to help meet future water demand by pumping water out of the Little Arkansas River (during above-base-flow conditions), treating it using National Primary Drinking Water Regulations as a guideline, and injecting it into the aquifer for later use. Water-quality data were collected and analyzed by the U.S. Geological Survey from 2 Little Arkansas River surface-water sites and 63 Equus Beds groundwater sites, including 38 areal assessment index wells, each of which has a shallow well and a deep well. About 4,700 surface and groundwater samples were collected and analyzed for more than 300 water-quality constituents. About 1,300 groundwater chemistry samples were geochemically modeled. Constituents of concern in the Equus Beds aquifer exceeded their respective Federal criteria throughout the study period and included chloride, sulfate, nitrate plus nitrite, Escherichia coli ( E. coli ), total coliforms, and dissolved iron and arsenic species.

Kansas↗

USGS2024

Results of 2018–19 water-quality and hydraulic characterization of aquifer intervals using packer tests and preliminary geophysical-log correlations for selected boreholes at and near the former Naval Air Warfare Center Warminster, Bucks County, Pennsylvania

The U.S. Geological Survey (USGS) collected data on the vertical distribution of hydraulic head, specific capacity, and water quality using aquifer-interval-isolation tests and other vertical profiling methods in 15 boreholes completed in fractured sedimentary bedrock in Northampton, Warminster, and Warwick Townships, Bucks County, Pennsylvania during 2018–19. This work was done, in cooperation with the U.S. Navy, to support detailed investigations at and near the former Naval Air Warfare Center (NAWC) Warminster, where groundwater contamination with per- and polyfluoroalkyl substances (PFAS) had become a concern since 2014. Two PFAS compounds, perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA), have been measured in groundwater samples from supply and monitoring wells at or near NAWC Warminster in concentrations above U.S. Environmental Protection Agency health advisory levels for drinking water. The area is underlain by the Triassic Stockton Formation, which predominantly consists of sandstone interbedded with shale and siltstone beds and forms a layered fractured-rock aquifer used for private, industrial, and public drinking water supply. The vertical distribution of aquifer properties and water quality was assessed through hydraulic tests and sampling of aquifer intervals using a straddle-packer system (13 boreholes) or depth-discrete point sampling under known borehole-flow conditions (2 boreholes). Geophysical and video logs collected by USGS during 2017–19 were used to identify potential water-bearing fractures in 15 boreholes, which ranged in depth from 210 to 604 feet (ft) and included 6 boreholes drilled in 2018 and 9 existing wells on or near the former NAWC Warminster. Measured borehole flow was predominantly downward in most of the deepest boreholes (greater than 400 ft), which were commonly located at the highest land-surface elevations, with inflow from fractures at relatively shallow depths and outflow through fractures near or below depths of 500 ft below land surface. Hydraulic head differences measured during packer tests were up to about 60 ft between shallow and deep intervals. Borehole flow was predominantly upward in most boreholes less than 400 ft in depth and farther from, and at lower land-surface elevations than, the former NAWC Warminster. Total borehole specific capacity ranged from about 0.07 to 41 gallons per minute per foot [(gal/min)/ft]. Specific-capacity values for individual intervals ranged from 0.02 to 40.0 (gal/min)/ft, with a median of 1.14 (gal/min)/ft and a large range in values at most depths. Differences in water quality of samples as indicated by field properties (pH, dissolved oxygen, and specific conductance) and concentrations of dissolved major ions, PFOA, and PFOS were apparent among isolated intervals in the boreholes. Summed concentrations of PFOA and PFOS ranged from about 11 to 10,780 nanograms per liter (ng/L) and were greater than the 2016 U.S. Environmental Protection Agency health advisory of 70 ng/L for summed PFOA and PFOS concentrations in 62 of 104 intervals and discrete depths tested. The mass ratio of PFOS to PFOA was generally higher than 1.0 in samples with summed PFOA and PFOS concentrations greater than 70 ng/L, with ratio values as high as 8.7. In many boreholes, summed concentrations of PFOA and PFOS were positively related to chloride concentrations, which were elevated above natural-background values [less than 10 milligrams per liter] in most samples and as high as 717 milligrams per liter. Sources of the elevated chloride other than, or in addition to, common rock salt (sodium chloride) were indicated by chloride to sodium molar ratios greater than 1.0. Water-quality data indicated that sampled water from some intervals with lower hydraulic heads may be affected by water from intervals with higher hydraulic heads because of vertical flow in open boreholes; samples from these intervals with lower hydraulic heads may not be fully representative due to some component of cross contamination and should be interpreted with caution. Through a preliminary correlation of natural gamma and resistivity logs of boreholes drilled at and near the former NAWC Warminster, 11 lithologic units were identified and interpreted to strike northeast and dip to the northwest. Hydraulic heads were generally highest in isolated intervals that intercepted beds which, when projected up dip, crop out at the highest land-surface elevation on the former NAWC Warminster, indicating that the dipping-bed structure and topography are factors affecting the distribution of hydraulic head in the aquifer. The hydrogeologic framework in conjunction with the vertical distribution of hydraulic heads and water quality may assist in evaluating the locations of various PFAS sources and potential migration pathways of PFAS in groundwater at and near NAWC Warminster.

Pennsylvania↗

USGS2009

Potentiometric surface of the Upper Floridan Aquifer in the St. Johns River Water Management District and vicinity, Florida, May 2009

This map depicts the potentiometric surface of the Upper Floridan aquifer in the St. Johns River Water Management District and vicinity for May 2009. Potentiometric contours are based on water-level measurements collected at 625 wells during the period May 14 - May 29, near the end of the dry season. Some contours are inferred from previous potentiometric-surface maps with larger well networks. The potentiometric surface of the carbonate Upper Floridan aquifer responds mainly to rainfall, and more locally, to groundwater withdrawals and spring flow. Potentiometric-surface highs generally correspond to topographic highs where the aquifer is recharged. Springs and areas of diffuse upward leakage naturally discharge water from the aquifer and are most prevalent along the St. Johns River. Areas of discharge are reflected by depressions in the potentiometric surface. Groundwater withdrawals locally have lowered the potentiometric surface. Groundwater in the Upper Floridan aquifer generally flows from potentiometric highs to potentiometric lows in a direction perpendicular to the contours.

Florida↗