USGS Science⌕ Search

SEARCH · USGS Science

Results for “Groundwater”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37Linked to original sources

Groundwater-level contour map of Fauquier County, Virginia, October-November 2018

Groundwater withdrawals provide most public-water supplies and all private-domestic users in Fauquier County, Virginia, a fast-growing rural area southwest of Washington, D.C. Groundwater levels were measured in 129 wells during a county-wide synoptic survey from October 29 through November 2, 2018. Field measurements, combined with datapoints from the National Hydrography Dataset, were used to develop a county-wide groundwater-level contour map. Groundwater levels and withdrawals during the synoptic survey were near or slightly above long-term medians. Error analysis indicated that the estimated groundwater-level contours generally were lower than observed measurements, with a root-mean-squared error of 33.52 feet. Groundwater levels in Fauquier County are controlled largely by topography: low levels in the crystalline Blue Ridge aquifers in the northwestern part of the county contrast markedly with higher levels in the sedimentary Mesozoic Basin aquifers in the southeast. At current levels of groundwater withdrawal, and at the scale and scope of the synoptic survey, no cones of depression in the groundwater surface were detected. The Fauquier County groundwater-level contour map is available as a U.S. Geological Survey data release.

Virginia↗

Machine-learning predictions of groundwater specific conductance in the Mississippi Alluvial Plain, south-central United States, with evaluation of regional geophysical aerial electromagnetic data as explanatory variables

The Mississippi Alluvial Plain, located in the south-central United States, is undergoing long-term groundwater-level declines within the surficial Mississippi River Valley alluvial aquifer (hereinafter referred to as “alluvial aquifer”), which has raised concerns about future groundwater availability. In some parts of the alluvial aquifer, groundwater availability for common uses such as irrigation, public supply, and domestic use is limited by quality (for example, high salinity) rather than quantity of water stored in the aquifer. The Mississippi Alluvial Plain region has an abundance of water-quality measurements in the alluvial aquifer and deeper aquifers; however, large areas lack direct measurements of salinity to evaluate regional groundwater availability. Statistical models can interpolate between wells to fill in spatial data gaps. In 2021, the U.S. Geological Survey trained two boosted regression tree (BRT) machine-learning models on specific conductance data available between 1942 and 2020 to predict spatially continuous surfaces of groundwater salinity at multiple depths for the alluvial aquifer and deeper aquifers. Well construction information, water levels, and surficial variables such as geomorphology and soils were included as explanatory variables in this baseline model. Additionally, subsurface electrical resistivity data from the first aquifer-wide aerial electromagnetic (AEM) survey for the region were incorporated to create a geophysical model. This work expands on prior BRT salinity predictions of the alluvial aquifer and extends predictions south to the Gulf of Mexico, where groundwater salinity is high. AEM survey data were not available for the southern extent of the alluvial aquifer at the time of modeling. A BRT model was trained without (baseline) and with (geophysical) AEM variables to test the ability of the models to predict salinity where explanatory data are missing and response data are sparse. Additionally, model sensitivity to AEM survey data was evaluated to better understand how AEM variables influence specific conductance predictions. Model performance was improved with the addition of geophysical data, which added three-dimensional information, thereby improving salinity predictions at depth. Groundwater specific conductance predictions can help inform other geophysical investigations in the southern extent of the study area, where high groundwater specific conductance can obfuscate changes in aquifer sediment resistivity and could limit groundwater resources for agricultural, public supply, and domestic uses.

Scientific Investigations Report↗

Trace metal and phosphorus loading from groundwater seepage into South Fork Coeur d’Alene River after remediation at the Bunker Hill Superfund Site, northern Idaho, 2022

Widely dispersed waste products from historical mining in northern Idaho’s Coeur d’Alene mining district have long been a concern in the Coeur d’Alene River Basin in northern Idaho. The Central Impoundment Area (CIA), an unlined mining waste repository that is part of the Bunker Hill Superfund Site designated in 1983, is adjacent to the South Fork Coeur d’Alene River between Kellogg and Smelterville, Idaho. Previous studies, including a pre-remediation seepage study completed by the U.S. Geological Survey (USGS) in 2017, have identified groundwater seepage from beneath the CIA as a major contributor to trace-metal and nutrient loads (including zinc, cadmium, and phosphorus) in the South Fork Coeur d’Alene River. A major remediation project, led by the U.S. Environmental Protection Agency from late 2017 to 2021, specifically aimed to reduce groundwater loading to the river via a groundwater collection system (GWCS) at the CIA. In 2022, the USGS completed a post-remediation seepage study to quantify zinc, cadmium, and phosphorus loading from groundwater to the South Fork Coeur d’Alene River in the same reach as the 2017 pre-remediation study. Like in the previous USGS study, discharge measurements and water-quality samples were collected during base-flow conditions in the South Fork Coeur d’Alene River between Kellogg and Smelterville as well as in surface-water inputs to the reach. Results of this study show a reduction in groundwater loads of dissolved zinc, dissolved cadmium, and total phosphorus entering the South Fork Coeur d’Alene River compared to 2017. The largest reductions in groundwater loading to the South Fork Coeur d’Alene River occurred in a discrete section (the middle section) of the reach adjacent to the CIA where the GWCS was expected to have the biggest impact. In the South Fork Coeur d’Alene River middle section, loads from groundwater (presented as a mean plus or minus [±] standard deviation) of dissolved zinc decreased from 85 ± 9.3 kilograms per day (kg/d) in 2017 to 11.6 ± 19.2 kg/d in 2022 (86-percent reduction), dissolved cadmium decreased from 0.59 ± 0.10 kg/d in 2017 to 0.11 ± 0.06 kg/d in 2022 (81-percent reduction), and total phosphorus decreased from 6.5 ± 0.45 kg/d in 2017 to 0.79 ± 0.97 kg/d in 2022 (88-percent reduction). In addition to reduced groundwater loading, lower concentrations of dissolved zinc, dissolved cadmium, and total phosphorus were observed at the site farthest downstream from the GWCS. Furthermore, the ambient water-quality-criteria ratios decreased at all river monitoring sites in 2022, although zinc and cadmium concentrations still exceeded the site-specific criteria designated to protect aquatic life. This post-remediation study indicates that the GWCS at the CIA has reduced groundwater loading of trace metals and phosphorus to the South Fork Coeur d’Alene River. This reduction in trace metals and phosphorus in South Fork Coeur d’Alene River also has implications for water quality downstream in the main-stem Coeur d’Alene River and in Coeur d’Alene Lake.

Idaho↗

Status of water quality in groundwater resources used for drinking-water supply in the southeastern San Joaquin Valley, 2013–15—California GAMA Priority Basin Project

The California Groundwater Ambient Monitoring and Assessment Program Priority Basin Project (GAMA-PBP) investigated water quality of groundwater resources used for drinking-water supplies in the Madera-Chowchilla, Kings, Kaweah, Tule, and Tulare Lake groundwater subbasins of the southeastern San Joaquin Valley during 2013–15. The study focused primarily on groundwater resources used for domestic-supply wells in the southeastern San Joaquin Valley (SESJV-D), which correspond mostly to shallower parts of aquifer systems, compared to the groundwater resources used for public-supply wells in the southeastern San Joaquin Valley (SESJV-P). The investigation had three components: (1) characterization of the status of water quality in the SESJV-D, (2) comparison between water quality in the SESJV-D and SESJV-P, and (3) identification of natural and anthropogenic factors that potentially could affect water quality in these resources. The characterization of water quality in the SESJV-D was based on data collected from 198 domestic wells sampled during 2013–15 by the U.S. Geological Survey (USGS); characterization of water quality in the SESJV-P was based on data collected from 124 wells sampled by the USGS during 2005–18 and an additional 1,577 wells with publicly available data reported to the California State Water Resources Control Board Division of Drinking Water (SWRCB-DDW). Measured concentrations were compared to regulatory and non-regulatory drinking-water quality benchmarks. A grid-based method was used to estimate the areal proportions of each study area and the whole southeastern San Joaquin Valley with high (greater than benchmark concentration), moderate (greater than half of the benchmark for inorganic and one-tenth of the benchmark for organic), and low concentrations relative to those benchmarks. Natural and anthropogenic factors that could affect groundwater quality for the SESJV-D were identified in the context of the hydrogeologic setting of the southeastern San Joaquin Valley. The considered factors represented hydrologic conditions and position in the groundwater flow system (well depth, lateral position, presence of hydric soils, percentage of coarse-grained sediment, and aridity index), land-use characteristics (percentages of agricultural, urban, and natural land use, percentage of orchard or vineyard land use, and densities of septic tanks and underground storage tanks near the wells), and geochemical conditions (groundwater age class, oxidation-reduction class, pH, and dissolved oxygen and bicarbonate concentrations). Factors are compared between SESJV-D and SESJV-P at the scale of the five study areas. One or more inorganic constituents with U.S. Environmental Protection Agency (EPA) or California maximum contaminant levels (MCLs) were detected at high concentrations in 47 percent of the SESJV-D and in 32 percent of the SESJV-P. The inorganic constituents most commonly present at high concentrations in the SESJV-D were nitrate, uranium, and arsenic. Within the SESJV-D, the proportion of the study area with high concentrations of inorganic constituents ranged from 19 percent in Madera-Chowchilla to 60 percent in Kings and Tulare Lake. One or more inorganic constituents with California State Water Resources Control Board Division of Drinking Water secondary maximum contaminant levels (SMCL-CAs) were detected at high concentrations in 14 percent of the SESJV-D and in 19 percent of the SESJV-P. The constituents most commonly present at high concentrations were manganese, iron, and total dissolved solids (TDS). Although the proportion of SESJV-D and SESJV-P with high concentrations of TDS greater than the upper SMCL were similar at 4 percent, the proportion of the SESJV-D with moderate concentrations (between the recommended and upper SMCL-CA), 30 percent, was greater than the proportion of the SESJV-P with moderate concentrations, 12 percent. One or more organic constituents with MCLs were present at high concentrations in 19 percent of the SESJV-D and in 12 percent of the SESJV-P. All the constituents detected at high concentrations in the SESJV-D were fumigants, primarily 1,2,3-trichloropropane (1,2,3-TCP) and 1,2-dibromo-3-chloropropane (DBCP). Fumigants also were the constituents most commonly detected at high concentrations in the SESJV-P, although high concentrations of solvents also were detected. The SESJV-D dataset included analysis of many organic constituents without MCL benchmarks and with detection levels far below drinking water benchmark concentrations; detections at these low concentrations can be used as tracers of anthropogenic influence on groundwater. Pesticides and degradates of pesticides were detected in 60 percent of the SESJV-D; the most frequently detected pesticides were the herbicides simazine, didealkylatrazine (CAAT, a degradate of simazine and atrazine), diuron, and bromacil.

California↗

Conceptual hydrogeologic framework and groundwater budget near the southeastern part of Puget Sound, Washington

More than 1 million people live within the active model area (AMA) in the southeastern part of the lowlands surrounding Puget Sound, or Puget Lowland, Washington, and groundwater is the source for approximately one-half of their public, domestic, and irrigation water demands. The 887-square-mile AMA, located in King and Pierce Counties, represents the area of analysis for the conceptual hydrogeologic framework and numerical groundwater-flow models within the study area and includes the Puyallup River and Chambers-Clover Creek watersheds. To assess the potential hydrologic and anthropogenic impacts to groundwater and the connected surface-water resources, conceptual and numerical groundwater-flow models of groundwater flow were developed by the U.S. Geological Survey Washington Water Science Center in close cooperation with 18 water-resource agencies and stakeholders. This report presents information used to characterize the groundwater-flow system and the development of a numerical model in the AMA. Included are descriptions of the geology and conceptual hydrogeologic framework, groundwater levels and flow directions, groundwater recharge and discharge, numerical groundwater-flow model construction and results, and model limitations. The study area encompasses the western part of Pierce County and the southwestern part of King County, Washington. The study area extends south to the Nisqually River, southwest to Tanwax Creek, northeast to the Green River, and north through the valley near Auburn and adjacent uplands. It is bounded on the east by foothills of the Cascade Range, and on the northwest by Puget Sound.

Washington↗

Groundwater quality near the Placerita Oil Field, California, 2018

Groundwater-quality data and potential fluid-migration pathways near the Placerita Oil Field in Los Angeles County, California, were examined by the U.S. Geological Survey to determine if oil-field fluids (water and gas from oil-producing and non-producing zones) have mixed with groundwater resources. Six of the 13 new groundwater samples collected for this study contained petroleum hydrocarbons, thermogenic gas, inorganic chemical signatures, and (or) isotopic values consistent with potential mixing with fluids from hydrocarbon-bearing formations. For historical groundwater samples, benzene was the most detected petroleum hydrocarbon. The historical groundwater samples with a benzene concentration greater than 0.5 micrograms per liter were from environmental monitoring wells at industrial or commercial facilities unrelated to oil and gas development that, in many cases, have identified soil or groundwater contamination and were not typically analyzed for other constituents that could provide additional lines of evidence for potential mixing with oil-field fluids. Methane was not detected in any of the 12 historical samples with a reported measurement. Reviewing historical data revealed factors that could potentially adversely affect groundwater quality in the study area. These factors include modified hydraulic gradients caused by large volumes of water extracted from the main production area and reinjected downgradient into nonproducing zones, well-barrier failures in wells constructed in the northern part of the oil field before the 1970s, well-barrier failures in produced-water disposal wells downgradient from the main production area, and naturally occurring hydrocarbons at shallow intervals. The groundwater samples most geochemically similar to samples from hydrocarbon-bearing formations were in areas where hydrocarbons are naturally occurring at shallow intervals and where oil development is at shallow depths. Additional data for hydraulic heads, water quality, and formation temperatures at multiple depths in areas with large injection volumes and well-integrity issues are needed to evaluate whether those factors have contributed to mixing between fluids from oil-producing or injection formations and groundwater resources.

California↗

Comparison of water quality in shallow groundwater near agricultural areas in the Delaware Coastal Plain, 2014 and 2019

The State of Delaware has encouraged agricultural conservation practices to improve nutrient uptake by crops and mitigate nutrient transport to groundwater in the surficial aquifer. To study recent changes in groundwater quality, the U.S. Geological Survey and the Delaware Department of Agriculture (DDA) developed a network of shallow wells near agricultural areas throughout the Delaware Coastal Plain. This network was designed to characterize water quality related to agricultural practices and to detect any recent changes in shallow groundwater quality, in particular groundwater nitrate concentrations. The shallow well network was first sampled in 2014 and resampled in 2019. In 2019, field parameters (including dissolved oxygen, pH, specific conductance, and temperature), major ions, nutrients, stable isotopes of water, and isotopes of nitrate were measured in groundwater samples collected between October and December. Wells were organized into three groups based on their geochemical characteristics measured in 2014: the Agricultural, Urban, and Mixed Groups. Results from the 2019 sampling showed little change in water quality from the 2014 sampling. Land-use factors continued to be the driving influence between groups. Groundwater moves slowly and changes in groundwater quality are likely to respond slowly to changes in conservation practices. Continued sampling of both groundwater quality in this network and monitoring land management practices can help detect groundwater quality trends in the future.

Delaware↗

Groundwater-storage change in the north Phoenix aquifer, Arizona, 2020–23

The city of Phoenix, Arizona, relies primarily on surface water for municipal water supply. The city also maintains wells to withdraw groundwater, particularly in times of drought and reduced surface-water supply, and to recharge groundwater when excess surface water is available. As of 2023, withdrawals from the aquifer in the northeastern part of the city are a small volume of water, less than 3,000 acre-feet in most years. Each year a similar volume of water is recharged through injection wells. The withdrawal wells are permitted to produce in total more than 20,000 acre-feet per year; increased pumping could lead to future groundwater-storage declines. To better understand groundwater-storage change in the north Phoenix aquifer, a repeat microgravity monitoring network was established in 2020. Measurements of changes in Earth’s gravity provide a direct, non-invasive measurement of subsurface mass change. Groundwater-storage changes were small during the 2020–23 study period, consistent with the relatively small volumes of pumping and recharge. Groundwater levels measured in monitoring wells were stable or increased slightly during this period, although the number of monitoring wells within the area of the gravity monitoring network is sparse. In total, about 15,000 acre-feet of water were pumped and 31,000 acre-feet recharged through injection in the north Phoenix aquifer during the 2020–22 period. Within the monitored area, groundwater storage increased by a small amount, about 1,000 acre-feet, in 2020, and decreased by a larger amount, about 6,000 acre-feet, each year in 2021 and 2022. Storage decreased at 89 of 102 stations from 2020 to 2023. Groundwater-storage decreases greater than the volume of net pumping indicate down-gradient subflow from the aquifer is greater than recharge plus incoming subflow, drying of the unsaturated zone resulting from decreased land-surface recharge, or both. At present (2023), the aquifer appears able to store and supply the relatively small amounts of water needed without excessive drawdown or harmful effects, such as land subsidence and surface inundation. If pumping or recharge increases in the future, the established repeat microgravity network is well suited to capture the spatial extent and magnitude of groundwater-storage changes in the aquifer.

Arizona↗

Passive sampling of groundwater wells for determination of water chemistry

Introduction Passive groundwater sampling is defined as the collection of a water sample from a well without the use of purging by a pump or retrieval by a bailer (Interstate Technology and Regulatory Council [ITRC], 2006; American Society for Testing and Materials [ASTM], 2014). No purging means that advection of water is not involved in collecting the water sample from the well. Passive samplers rely on diffusion as the primary process that drives their collection of chemical constituents. Diffusion is the transport of chemicals caused by the presence of a chemical gradient. Chemicals tend to move or diffuse from areas of higher concentration to areas of lower concentration to reach an average or equilibrium concentration. Passive sampling of groundwater relies on the ambient exchange of groundwater in the formation with water in the screened or open interval of a well. In this report, the term formation is used to describe all saturated hydrogeologic units that yield water to a well. If the well opening is unclogged and free of a film of deposits from physical turbidity or chemical precipitation, then the exchange of groundwater is likely adequate, and the water in the open interval will be representative of water in the formation. In some cases, the passive sample from the well opening can be more representative of groundwater from the formation than a sample collected by pumping if pumping induces mixing of water in the open interval with stagnant casing water that has undergone chemical alteration (Harte and others, 2018). In most cases, passive sampling will better represent the ambient groundwater chemistry flowing through the open interval of a well because pumping may capture water of different chemistry from downgradient or lateral areas that would not normally pass through the well. Three basic types of passive samplers are discussed in this report. The first type of passive sampler is the equilibrium-membrane type, which includes a semi-permeable membrane through which chemicals diffuse or permeate. Permeation is simply the process of water or chemicals moving through openings in the membrane. The authors contend that permeation is dominated by diffusion for many of the passive samplers discussed in this report. Some passive equilibrium-membrane-type samplers allow most types of chemical constituents through, whereas others allow the diffusion of only selected groups of chemicals. Once the chemical constituents are inside the membrane, they are retained by the equilibration of concentrations inside the sampler with those outside the sampler. The second type of passive sampler is an equilibrium-thief type, which has no semi-permeable membrane. Chemical constituents simply move through the openings in the body of the sampler either initially through advection and dispersion or over time primarily by diffusion. Chemical constituents reach equilibrium between the water in the sampler and the water in the well and are captured in the sampler when the sampler is closed. The third type of passive sampler is an accumulation-type sampler that contains sorptive media. Selected chemical constituents are sorbed onto the media that the sampler contains for later extraction and analysis. Although passive samplers have been available for more than 15 years (from present [2020]), their use by U.S. Geological Survey (USGS) hydrologists and hydrologic technicians to monitor groundwater quality largely has been limited to selected research studies. The authors believe that this may be the result of (1) a lack of exposure of most USGS personnel to passive samplers and the uses of these samplers and (2) the lack of a USGS-approved protocol for the proper use of these samplers by USGS personnel. This report is an effort to fill those two needs. The focus of this report is on hydraulic, hydrologic, and chemical considerations in the application of passive samplers and interpretation of groundwater chemistry results obtained using passive samplers in wells. This report describes the differences between purging and passive sampling methods in groundwater and explains how and why passive samplers work. The report points out the advantages and limitations of passive samplers in general and for each particular type of passive sampler. Important considerations to be taken into account prior to the use of passive samplers are discussed, such as defining the data-quality objectives, the water-quality constituents to be sampled, sample volumes required for analysis, well construction of the sampling network, and the geologic formations that will be sampled. Potential applications of passive samplers also are discussed, such as chemical-vertical profiling of wells. A general field protocol for the deployment, recovery, and sample collection using these devices is described, and some overall guidance for the practitioner with application examples is given. Comparison methods used to evaluate results from passive sampling versus purge sampling also are discussed.

Techniques and Methods↗

TracerLPM (Version 1): An Excel® workbook for interpreting groundwater age distributions from environmental tracer data

TracerLPM is an interactive Excel® (2007 or later) workbook program for evaluating groundwater age distributions from environmental tracer data by using lumped parameter models (LPMs). Lumped parameter models are mathematical models of transport based on simplified aquifer geometry and flow configurations that account for effects of hydrodynamic dispersion or mixing within the aquifer, well bore, or discharge area. Five primary LPMs are included in the workbook: piston-flow model (PFM), exponential mixing model (EMM), exponential piston-flow model (EPM), partial exponential model (PEM), and dispersion model (DM). Binary mixing models (BMM) can be created by combining primary LPMs in various combinations. Travel time through the unsaturated zone can be included as an additional parameter. TracerLPM also allows users to enter age distributions determined from other methods, such as particle tracking results from numerical groundwater-flow models or from other LPMs not included in this program. Tracers of both young groundwater (anthropogenic atmospheric gases and isotopic substances indicating post-1940s recharge) and much older groundwater (carbon-14 and helium-4) can be interpreted simultaneously so that estimates of the groundwater age distribution for samples with a wide range of ages can be constrained. TracerLPM is organized to permit a comprehensive interpretive approach consisting of hydrogeologic conceptualization, visual examination of data and models, and best-fit parameter estimation. Groundwater age distributions can be evaluated by comparing measured and modeled tracer concentrations in two ways: (1) multiple tracers analyzed simultaneously can be evaluated against each other for concordance with modeled concentrations (tracer-tracer application) or (2) tracer time-series data can be evaluated for concordance with modeled trends (tracer-time application). Groundwater-age estimates can also be obtained for samples with a single tracer measurement at one point in time; however, prior knowledge of an appropriate LPM is required because the mean age is often non-unique. LPM output concentrations depend on model parameters and sample date. All of the LPMs have a parameter for mean age. The EPM, PEM, and DM have an additional parameter that characterizes the degree of age mixing in the sample. BMMs have a parameter for the fraction of the first component in the mixture. An LPM, together with its parameter values, provides a description of the age distribution or the fractional contribution of water for every age of recharge contained within a sample. For the PFM, the age distribution is a unit pulse at one distinct age. For the other LPMs, the age distribution can be much broader and span decades, centuries, millennia, or more. For a sample with a mixture of groundwater ages, the reported interpretation of tracer data includes the LPM name, the mean age, and the values of any other independent model parameters. TracerLPM also can be used for simulating the responses of wells, springs, streams, or other groundwater discharge receptors to nonpoint-source contaminants that are introduced in recharge, such as nitrate. This is done by combining an LPM or user-defined age distribution with information on contaminant loading at the water table. Information on historic contaminant loading can be used to help evaluate a model's ability to match real world conditions and understand observed contaminant trends, while information on future contaminant loading scenarios can be used to forecast potential contaminant trends.

Techniques and Methods↗

Insight into sources of benzene, TCE, and PFOA/PFAS in groundwater at Naval Air Station Whiting Field, Florida, through numerical particle-tracking simulations

Past waste-disposal activities at Naval Air Station Whiting Field (NASWF) have led to elevated concentrations of contaminants in the underlying sand and gravel aquifer. Contaminants include two of the most commonly detected chemicals in groundwater in many countries (benzene and trichloroethylene (TCE)) and the “forever chemicals” per- and poly-fluoroalkyl substances (PFAS) such as perfluorooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS). A MODFLOW model (the Whiting Field Groundwater Model (WFGM)) was previously developed for NASWF and the surrounding area to simulate groundwater flow. To obtain insight into groundwater flow pathways for the identification of potential source areas, the MODPATH particle-tracking application was applied to the WFGM for three public supply wells and three monitoring wells at NASWF. The travel time to recharge areas was estimated using concentrations of the groundwater age-dating solutes tritium (as helium ingrowth) and chlorofluorocarbons detected in the monitoring wells. Simulated travel times agree with the groundwater ages and indicate that the calibrated WFGM reasonably represents groundwater flow velocities and pathways. The MODPATH simulations confirm suspected on-base source areas to explain chemical detection in the monitoring wells. In contrast, the particle-tracking simulations indicate that potential source areas to the public supply wells include both on- and off-base sources. This is important because PFAS chemicals can have multiple sources, including land application of sludge-based fertilizers. This approach that combines groundwater age dating with particle-tracking simulations can be applied at similar sites characterized by benzene-, TCE-, and PFAS-contaminated groundwater.

Florida↗

Identification of persisten anionic surfactant-derived chemicals in sewage effluent and groundwater

Preparative isolation and fractionation procedures coupled with spectrometric analyses were used to identify surfactant-derived contaminants in sewage effluent and sewage-contaminated groundwater from a site located on Cape Cod, Massachusetts. Anionic surfactants and their biodegradation intermediates were isolated from field samples by ion exchange and fractionated by solvent extraction and adsorption chromatography. Fractions were analyzed by 13C nuclear magnetic resonance spectrometry and gas chromatography-mass spectrometry. Carboxylated residues of alkylphenol polyethoxylate surfactants were detected in sewage effluent and contaminated groundwater. Linear alkylbenzenesulfonates (LAS) were identified in sewage effluent and groundwater. Groundwater LAS composition suggested preferential removal of select isomers and homologs due to processes of biodegradation and partitioning. Tetralin and indane sulfonates (DATS), alicyclic analogs of LAS, were also identified in field samples. Although DATS are a minor portion of LAS formulations, equivalent concentrations of LAS and DATS in groundwater suggested persistence of alicyclic contaminant structures over those of linear structure. Sulfophenyl-carboxylated (SPC) LAS biodegradation intermediates were determined in sewage effluent and groundwater. Homolog distributions suggested that SPC containing 3-10 alkyl-chain carbons persist during infiltration and groundwater transport. Surfactant-derived residues detected in well F300-50 groundwater have a minimum residence time in the range of 2.7-4.6 yr. LAS detected in groundwater at 500 m from infiltration has been stable over an estimated 50-500 half lives.

Conference Paper↗

Quantifying groundwater discharge through fringing wetlands to estuaries: Seasonal variability, methods comparison, and implications for wetland-estuary exchange

Because groundwater discharge along coastal shorelines is often concentrated in zones inhabited by fringing wetlands, accurately estimating discharge is essential for understanding its effect on the function and maintenance of these ecosystems. Most previous estimates of groundwater discharge to coastal wetlands have been temporally limited and have used only a single approach to estimate discharge. Furthermore, groundwater input has not been considered as a major mechanism controlling pore-water flushing. We estimated seasonally varying groundwater discharge into a fringing estuarine wetland using three independent methods (Darcy's Law, salt balance, and Br- tracer). Seasonal patterns of discharge predicted by both Darcy's Law and the salt balance yielded similar seasonal patterns with discharge maxima and minima in spring and early fall, respectively. They differed, however, in the estimated magnitude of discharge by two- to fourfold in spring and by 10-fold in fall. Darcy estimates of mean discharge ranged between -8.0 and 80 L m-2 d-1, whereas the salt balance predicted groundwater discharge of 0.6 to 22 L m-2 d-1. Results from the Br- tracer experiment estimated discharge at 16 L m-2 d-t, or nearly equal to the salt balance estimate at that time. Based upon the tracer test, pore-water conductivity profiles, and error estimates for the Darcy and salt balance approaches, we concluded that the salt balance provided a more certain estimate of groundwater discharge at high flow (spring). In contrast, the Darcy method provided a more reliable estimate during low flow (fall). Groundwater flushing of pore water in the spring exported solutes to the estuary at rates similar to tidally driven surface exchange seen in previous studies. Based on pore-water turnover times, the groundwater-driven flux of dissolved organic carbon (DOC), dissolved organic nitrogen (DON), and NH4+ to the estuary was 11.9, 1.6, and 1.3 g C or g N m-2 wetland for the 90 d encompassing peak spring discharge. Groundwater-induced flushing of the wetland subsurface therefore represents an important mechanism by which narrow fringing marshes may seasonally relieve salt stress and export material to adjacent water masses.

Limnology and Oceanography↗

The 2016 groundwater flow model for Dane County, Wisconsin

A new groundwater flow model for Dane County, Wisconsin, replaces an earlier model developed in the 1990s by the Wisconsin Geological and Natural History Survey (WGNHS) and the U.S. Geological Survey (USGS). This modeling study was conducted cooperatively by the WGNHS and the USGS with funding from the Capital Area Regional Planning Commission (CARPC). Although the overall conceptual model of the groundwater system remains largely unchanged, the incorporation of newly acquired high-quality datasets, recent research findings, and improved modeling and calibration techniques have led to the development of a more detailed and sophisticated model representation of the groundwater system. The new model is three-dimensional and transient, and conceptualizes the county’s hydrogeology as a 12-layer system including all major unlithified and bedrock hydrostratigraphic units and two high-conductivity horizontal fracture zones. Beginning from the surface down, the model represents the unlithified deposits as two distinct model layers (1 and 2). A single layer (3) simulates the Ordovician sandstone and dolomite of the Sinnipee, Ancell, and Prairie du Chien Groups. Sandstone of the Jordan Formation (layer 4) and silty dolostone of the St. Lawrence Formation (layer 5) each comprise separate model layers. The underlying glauconitic sandstone of the Tunnel City Group makes up three distinct layers: an upper aquifer (layer 6), a fracture feature (layer 7), and a lower aquifer (layer 8). The fracture layer represents a network of horizontal bedding-plane fractures that serve as a preferential pathway for groundwater flow. The model simulates the sandstone of the Wonewoc Formation as an upper aquifer (layer 9) with a bedding-plane fracture feature (layer 10) at its base. The Eau Claire aquitard (layer 11) includes shale beds within the upper portion of the Eau Claire Formation. This layer, along with overlying bedrock units, is mostly absent in the preglacially eroded valleys along the Yahara River valley and in northeastern Dane County. Layer 12 represents the Mount Simon sandstone as the lowermost model layer. It directly overlies the Precambrian crystalline basement rock, whose top surface forms the lower boundary of the model. The model uses the USGS MODFLOW-NWT finite-difference code, a standalone version of MODFLOW-2005 that incorporates the Newton (NWT) solver. MODFLOW-NWT improves the handling of unconfined conditions by smoothing the transition from wet to dry cells. The model explicitly simulates groundwater–surface-water interaction with streamflow routing and lake-level fluctuation. Model input included published and unpublished hydrogeologic data from recent estimates of aquifer hydraulic conductivities. A spatial groundwater recharge distribution was obtained from a recent GIS-based, soil-water-balance model for Dane County. Groundwater withdrawals from pumping were simulated for 572 wells across the entire model domain, which includes Dane County and portions of seven neighboring counties—Columbia, Dodge, Green, Iowa, Jefferson, Lafayette, and Rock. These wells withdrew an average of 60 million gallons per day (mgd) over the 5-year period from 2006 through 2010. Within Dane County, 385 wells were simulated with an average withdrawal rate of 52 mgd. Model calibration used the parameter estimation code PEST, and calibration targets included heads, stream and spring flows, lake levels, and borehole flows. Steady-state calibration focused on the period 2006 through 2010; the transient calibration focused on the 7-week drought period from late May through July 2012. This model represents a significant step forward from previous work because of its finer grid resolution, improved hydrostratigraphic discretization, transient capabilities, and more sophisticated representation of surface-water features and multi-aquifer wells. Potential applications of the model include evaluation of potential sites for and impacts of new high-capacity wells, development of wellhead protection plans, evaluating the effects of changing land use and climate on groundwater, and quantifying the relationships between groundwater and surface water.

Wisconsin↗

Interstate water management of a “hidden” resource - Physical principles of groundwater hydrology

Groundwater systems are dynamic geologic environments in which water continuously flows from recharge areas to discharge areas at streams, springs, wetlands, coastal waters, and wells. Natural, predevelopment conditions within groundwater systems are changed by the introduction of wells and other human stresses that modify existing groundwater levels, flow paths, and hydrologic budgets. Groundwater serves the Nation as an important water supply, but in some instances such stresses can have adverse impacts that include excessive ground water-level declines, aquifer-storage reductions, and streamflow depletions. Many of the Nation’s aquifer systems extend over thousands of square miles and their hydrologic boundaries may be distant from jurisdictional boundaries that can be the focus of groundwater disputes. Effective interstate management of groundwater resources is benefited by an understanding of the regional-scale controls that affect groundwater conditions at the local scale. Numerical models are the most effective approach for accounting for all of the relevant hydrologic processes that affect groundwater systems and their response to natural and manmade stresses. This paper provides a brief background on some of the basic principles of groundwater hydrology that are relevant to interstate management of this important natural resource.

Conference Paper↗

Groundwater modeling

The state of the science and practice in groundwater modeling brings to mind highly sophisticated computer models that are running in parallel on many multi-processor machines. These models are expected to incorporate many different processes of both saturated and unsaturated groundwater flow and transport and possibly the media to which it connects, like surface waters and the atmosphere. We are increasingly aware we cannot study groundwater flow in isolation if we are to make useful predictions of, for instance, the impacts of climate change on the groundwater regime. We have come a long way. Today we are no longer limited to equations for flow toward a well, perhaps near an infinitely long straight canal (method of images), to sandbox models in the laboratory, or to simple steady state models of flow in a single aquifer. We now have computer models that solve groundwater flow and transport in multi-aquifer settings under transient conditions and with a user-friendly graphical user interface that allows widespread use. Additionally, multi-media models are now leaving the research environment and becoming available to mainstream consultants. So in that sense the science of groundwater modeling has matured. The practice of groundwater modeling, however, has also matured. We have come to realize that model output, being a necessary simplification of an unknowably complex natural world, has inherent limitations. That is, a model of reality is not reality itself. There is uncertainty associated with all facets of our model—parameterization, aquifer geometry and discretization, boundary conditions, and future hydrologic drivers such as future pumping regimes and climates. Today a model is now more appropriately seen as a tool that provides a quantitative framework to make supportable forecasts rather than an oracle that gives us all the answers. In this chapter we set out to briefly review the state of the science and practice in modeling. In doing so, we augment existing assessments from the journal Groundwater (e.g., Hunt and Zheng 2012; Langevin and Panday 2012; Molz 2017a,b; White 2017), specifically in terms of modeling approach. An effective modeling approach is critical. If a modeler does not decompose the societal problem correctly, the model will not be fit-for-purpose, no matter how sophisticated the code’s capabilities. Moreover, capabilities of codes will be ever improving; good modeling practices have a timelessness that is more robust. How best to decompose the problem and provide models that are accepted? We lay out here some approaches for today’s applied groundwater modeling. Specifically, we suggest: (1) a step-wise modeling process; (2) including a two-dimensional areal model within this process; (3) keeping abreast of industry standards; and (4) ways to increase acceptance of the models we produce.

Book chapter↗

Geophysical characterization of mine influenced groundwater and surface water in the Mayflower section of the Animas River, Bonita Peak Mining District, Silverton Colorado

This report details findings from geophysical investigations to identify possible groundwater - surface water interactions near the Mayflower Section of the Animas River in Silverton, Colorado. The Mayflower Section is a mine influenced Superfund Site in Colorado. This investigation utilized electromagnetic induction (EM or EMI), magnetic, and fiber optic distributed temperature system (FODTS) geophysical methods to measure the bulk earth electrical conductivity, magnetic susceptibility, and temperature of specific surveyed volumes of the earth near the Mayflower Section of the Animas River. These physical parameters are used to understand the groundwater – surface water interactions, which can guide decision makers in their assessment of mine-impacted surface water. This report details the results from characterization and monitoring technologies to provide high data density and continuous monitoring of bulk earth electrical conductivity and temperature in the Mayflower section of the Animas River to identify zones of groundwater – surface water interactions and potential metal loading from mine influenced water. The investigation separated right and left bank characterization for each method and indicates more groundwater is entering from the right bank than the left bank and these predominantly right bank discharges potentially contain metal-rich water compared to the left bank. Results also indicate mineral veins facilitate preferential groundwater discharge to the river due to possible jointing, fractures, and permeability differences sometimes occurring along veins relative to host rock. For example, Boulder Gulch is likely groundwater dominated and may be receive water impacted by Mayflower Impoundments #1 and #2. Additionally, the beaver ponds near Blair Gulch may influence groundwater discharge and Mayflower Impoundment #4 is possibly impacting groundwater and surface water and may be connected to the wetlands to the west of the impoundment. These data could be further analyzed for smaller spatial scale analysis within areas of interest. The identification of these locations along sections of the river likely impacted by mine influenced groundwater potentially entering the Animas River and can be used by site investigators, decision makers, and stakeholders in mitigation decisions and strategies.

Colorado↗

Quality and age of shallow groundwater in the Bakken Formation production area, Williston Basin, Montana and North Dakota

The quality and age of shallow groundwater in the Bakken Formation production area were characterized using data from 30 randomly distributed domestic wells screened in the upper Fort Union Formation. Comparison of inorganic and organic chemical concentrations to health based drinking-water standards, correlation analysis of concentrations with oil and gas well locations, and isotopic data give no indication that energy-development activities affected groundwater quality. It is important, however, to consider these results in the context of groundwater age. Most samples were recharged before the early 1950s and had 14C ages ranging from 30,000 years. Thus, domestic wells may not be as well suited for detecting contamination associated with recent surface spills as shallower wells screened near the water table. Old groundwater could be contaminated directly by recent subsurface leaks from imperfectly cemented oil and gas wells, but horizontal groundwater velocities calculated from 14C ages imply that the contaminants would still be less than 0.5 km from their source. For the wells sampled in this study, the median distance to the nearest oil and gas well was 4.6 km. Because of the slow velocities, a long-term commitment to groundwater monitoring in the upper Fort Union Formation is needed to assess the effects of energy development on groundwater quality. In conjunction with that effort, monitoring could be done closer to energy-development activities to increase the likelihood of early detection of groundwater contamination if it did occur.

Montana, North Dakota↗