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 793 records · Page 44Linked to original sources

USGS2016

Geology, hydrology, water quality, and potential for interbasin invasive-species spread by way of the groundwater pathway near Lemont, Illinois

Invasive species such as Asian carps have the potential to travel in the egg, larval, or fry stages from the Des Plaines River (DPR) to the Chicago Sanitary and Ship Canal (CSSC) by way of the network of secondary-permeability features in the dolomite aquifer between these water bodies. Such movement would circumvent the electric fish barrier on the canal and allow Asian carps to travel unimpeded into Lake Michigan. This potential pathway for the spread of Asian carps and other invasive species was evaluated by the U.S. Geological Survey. The bed of the DPR appears to be in at least partial contact with the exposed bedrock in most of the area from about 1 mile west of Kingery Highway to Romeo Road (the study area). Areas of exposed bedrock are the most likely places for Asian carps to enter the groundwater system from the DPR. Water levels in the DPR typically are about 7–16 feet higher than those in the CSSC in most of the study area. This difference in water level provides the driving force for the potential spread of Asian carps from the DPR to the CSSC by way of groundwater. Groundwater flow (and potentially invasive-species movement) is through an interconnected network of permeable vertical and horizontal fractures within the Silurian dolomite bedrock. At least some of the fractures are associated with paleo-karst features. Several investigative techniques identified horizontal permeable fractures at about 546–552 feet above the North American Vertical Datum of 1988 within about 55 feet of the CSSC in the focus area between Lemont Road and Interstate 355. The elevation of the bottom of the CSSC in this area is about 551 feet, indicating that a direct conduit for flow of groundwater to the CSSC may be present. Wells further away from the CSSC in this area do not intercept fractures, so the fracture network may not be continuous between the DPR and the CSSC. These data are consistent with field observations of the secondary-permeability network along the CSSC walls, which indicate that the secondary-permeability features are completely filled with Pennsylvanian sediments within a few feet of the canal wall. Water-level data indicate the potential for flow from the DPR into the Silurian aquifer in the focus area, then from the aquifer to the CSSC. Water-level data also indicate that the fractures within the aquifer in the focus area are hydraulically well connected to the CSSC but not to the DPR, indicating that flow from the DPR to the groundwater system may not be substantial or rapid. Water-quality data in the CSSC and the DPR show similar values and trends and are affected by diel and longer term variations in climate and precipitation. However, the values and trends in water quality in the groundwater system tended to be substantially different from those in the DPR and the CSSC, indicating that the DPR and the CSSC do not appreciably recharge the groundwater system. Water-quality and flow data do indicate that groundwater discharges to the CSSC in part of the focus area. The absence of substantial hydraulic interaction between the groundwater and the DPR is supported by the absence of detectable concentrations of the dye tracer added to the DPR in groundwater in the focus area, which indicates that water from the DPR requires more than 2 weeks to move into the monitored parts of the groundwater system under approximately typical hydraulic conditions. The totality of the data indicates that there is minimal potential for the inter-basin spread of Asian carps by way of the groundwater pathway between Romeo Road and Stickney, Illinois.

Illinois↗

USGS2016

Hydrological conditions and evaluation of sustainable groundwater use in the Sierra Vista Subwatershed, Upper San Pedro Basin, southeastern Arizona

This study assessed progress toward achieving sustainable groundwater use in the Sierra Vista Subwatershed of the Upper San Pedro Basin, Arizona, through evaluation of 14 indicators of sustainable use. Sustainable use of groundwater in the Sierra Vista Subwatershed requires, at a minimum, a stable rate of groundwater discharge to, and thus base flow in, the San Pedro River. Many of the 14 indicators are therefore related to long-term or short-term effects on base flow and provide us with a means to evaluate groundwater discharge to and base flow in the San Pedro River. The indicators were based primarily on 10 to 20 years of data monitoring in the subwatershed, ending in 2012, and included subwatershedwide indicators, riparian-system indicators, San Pedro River indicators, and springs indicators. Groundwater management actions including voluntary retirement of irrigation pumping in the subwatershed resulted in about a 5,100 acre-feet (acre-ft) reduction in net human use from 2002 to 2012. Subwatershed population increased more than 10,000 during the same period. Most of the reduction occurred during 2002–07 and included reductions in groundwater pumping and increases in managed recharge; net human use varied annually by a few hundred acre-ft during 2007–12. The groundwater budget for 2012 showed a deficit of about 5,000 acre-ft, although the total water-budget uncertainty was about 5,500 acre-ft. In the vicinity of the U.S. Army’s Fort Huachuca, regional-aquifer water levels were in steady decline beginning in at least the mid-1990s (in older wells since at least the early-1970s), as the cone of depression centered on the Sierra Vista and Fort Huachuca pumping centers continued to deepen. This was evident in the individual water levels on Fort Huachuca, as well as from the horizontal hydraulic gradients that extend from the pumping centers toward the San Pedro and Babocomari Rivers. Basin water levels in wells southeast of Sierra Vista, away from the river, were also experiencing declines, while some water levels closer to the river were rising. Near-stream vertical gradients along the San Pedro River showed no clear increasing or decreasing trends that would indicate a shift in the direction of subsurface flow between the riverbed and the alluvial aquifer, or a trend in the magnitude of groundwater/surface-water exchange. Annual streamflow permanence data showed no clear change in streamflow permanence trends in any of the river reaches, other than those related to precipitation trends. Similarly, the single-day, dry-season, wet-dry streamflow analysis of all subwatershed river reaches indicated no change in condition over the past 14 years, with the exception of the Hereford reach, which has seen a statistically significant increase in wetted length. Dry-season, alluvial-aquifer water levels in the Hereford reach also showed a statistically significant increase. These improvements are attributed to the end of irrigation pumping in the area. Although data indicate that the length of the Fairbank North wetted reach may be in decline, it is not yet statistically significant. Stable-isotope data indicated reduced groundwater discharge to the Babocomari River in the vicinity of the Babocomari River near Tombstone gaging station and to the San Pedro River near the San Pedro River at Palominas gaging station and near the Lewis Springs DCP stage recorder. The Babocomari River near Tombstone gaging station is downgradient of the major pumping centers. The change in isotopic signature at the Lewis Springs stage recorder could have been the result of alterations in groundwater/surface-water interactions there caused by beaver damming of the river. Base flow in the San Pedro River declined over the periods of record at the three San Pedro River gaging stations in the subwatershed (Palominas, Charleston, and Tombstone), as well as at the Babocomari River near Tombstone gaging station. Precipitation declined slightly from the 1990s to the 2000s, although there is no statistically significant trend in subwatershed precipitation from 1991 to 2012. The occurrence of large winter discharge events appeared to decline and that of large summer discharge events appeared to increase over this same period. Data for physical parameters, general chemistry, nutrient species, select trace elements, and suspended sediment were collected at San Pedro River at Charleston stream-gaging station. These data were summarized over time and analyzed in relation to discharge and season as a means to assess trends over the period of analysis. Federal and State of Arizona drinking-water and human-contact standards were all met and few exceedances occurred for the ecological thresholds investigated. Several constituents showed a significant trend over the period of analysis, but only concentration and flux data for total phosphate, orthophosphate, total nitrogen, suspended sediment, and sulfate were suitable to be used in a weighted regression analysis that statistically accounted for time, discharge, and season. Sulfate concentrations and flux showed a significant downward trend over the period of analysis, whereas total phosphorus and ortho-phosphate showed a relatively small magnitude upward trend relative to standards. Suspended sediment concentrations and flux both showed a significant downward trend in the 1980s, an effect attributed to reduction of cattle in the subwatershed at about this time, and (or) increased cottonwood ( Populus fremontii ) and willow ( Salix goodingii ) recruitment, and (or) the curtailment of sand and gravel mining adjacent to the San Pedro River with the designation of the San Pedro Riparian National Conservation Area in 1988. A spike in sediment flux in 2006 may be attributable to the more than 100 debris flows in the Huachuca Mountains during the summer monsoon of that year. Spring discharge along the San Pedro River generally increased at three sites proximate to the Sierra Vista treated effluent recharge facility and varied somewhat with climate at two other sites. Median annual discharge at the recharge facility peaked in 2006, and at Murray Springs and Horsethief Spring, downgradient of the recharge facility, in 2009. Sampling for trace organic compounds in flow from springs was carried out using both discrete sampling and passive sampling methods. Spring samples thus collected showed the presence of trace-organic compounds. Lewis Springs (background site) had the least number of detections, whereas Murray Springs, located directly downgradient of the City of Sierra Vista’s treated effluent recharge facility, had the greatest number of detections of all the springs. Discrete samples from the recharge facility had more than twice the detections found in discrete samples from Murray Spring and at much higher concentrations. Few similar trace-organic compounds were detected at both the springs and the treated effluent recharge facility, and the number of detections did not increase during the collection period. Limitations of the study prevented the determination of trace-organic concentration in passive samplers and also prevented linking trace organic compounds detected at the treated effluent recharge facility with compounds detected from the springs. In particular, trace organic compounds could also derive from other sources such as septic systems. Looking at the subwatershed as a whole, base flow was in decline along the entire river reach, but determination of the specific cause of the decline was beyond the scope of this report. Conditions in the area from the municipal pumping center of Sierra Vista and Fort Huachuca northeast to the river (from about the Charleston to Tombstone gaging stations) were more commonly in decline than in regions further south. Both long-term indicators, such as regional aquifer groundwater levels and horizontal gradients, and the isotope analysis indicated that groundwater discharge to the river and thus base flow may continue to decline in that area. South of Charleston, indicators were more mixed. Some indicators in the Hereford reach suggest groundwater discharge to the San Pedro River may be increasing there, whereas some indicators in the Palominas reach suggest groundwater discharge to the river there may be declining.

Arizona↗

USGS2016

A review of single-sample-based models and other approaches for radiocarbon dating of dissolved inorganic carbon in groundwater

Numerous methods have been proposed to estimate the pre-nuclear-detonation 14 C content of dissolved inorganic carbon (DIC) recharged to groundwater that has been corrected/adjusted for geochemical processes in the absence of radioactive decay ( 14 C 0 ) - a quantity that is essential for estimation of radiocarbon age of DIC in groundwater. The models/approaches most commonly used are grouped as follows: (1) single-sample-based models, (2) a statistical approach based on the observed (curved) relationship between 14 C and δ 13 C data for the aquifer, and (3) the geochemical mass-balance approach that constructs adjustment models accounting for all the geochemical reactions known to occur along a groundwater flow path. This review discusses first the geochemical processes behind each of the single-sample-based models, followed by discussions of the statistical approach and the geochemical mass-balance approach. Finally, the applications, advantages and limitations of the three groups of models/approaches are discussed. The single-sample-based models constitute the prevailing use of 14 C data in hydrogeology and hydrological studies. This is in part because the models are applied to an individual water sample to estimate the 14 C age, therefore the measurement data are easily available. These models have been shown to provide realistic radiocarbon ages in many studies. However, they usually are limited to simple carbonate aquifers and selection of model may have significant effects on 14 C 0 often resulting in a wide range of estimates of 14 C ages. Of the single-sample-based models, four are recommended for the estimation of 14 C 0 of DIC in groundwater: Pearson's model, (Ingerson and Pearson, 1964; Pearson and White, 1967), Han & Plummer's model (Han and Plummer, 2013), the IAEA model (Gonfiantini, 1972; Salem et al., 1980), and Oeschger's model (Geyh, 2000). These four models include all processes considered in single-sample-based models, and can be used in different ranges of 13 C values. In contrast to the single-sample-based models, the extended Gonfiantini & Zuppi model (Gonfiantini and Zuppi, 2003; Han et al., 2014) is a statistical approach. This approach can be used to estimate 14 C ages when a curved relationship between the 14 C and 13 C values of the DIC data is observed. In addition to estimation of groundwater ages, the relationship between 14 C and δ 13 C data can be used to interpret hydrogeological characteristics of the aquifer, e.g. estimating apparent rates of geochemical reactions and revealing the complexity of the geochemical environment, and identify samples that are not affected by the same set of reactions/processes as the rest of the dataset. The investigated water samples may have a wide range of ages, and for waters with very low values of 14 C, the model based on statistics may give more reliable age estimates than those obtained from single-sample-based models. In the extended Gonfiantini & Zuppi model, a representative system-wide value of the initial 14 C content is derived from the 14 C and δ 13 C data of DIC and can differ from that used in single-sample-based models. Therefore, the extended Gonfiantini & Zuppi model usually avoids the effect of modern water components which might retain ‘bomb’ pulse signatures. The geochemical mass-balance approach constructs an adjustment model that accounts for all the geochemical reactions known to occur along an aquifer flow path (Plummer et al., 1983; Wigley et al., 1978; Plummer et al., 1994; Plummer and Glynn, 2013), and includes, in addition to DIC, dissolved organic carbon (DOC) and methane (CH 4 ). If sufficient chemical, mineralogical and isotopic data are available, the geochemical mass-balance method can yield the most accurate estimates of the adjusted radiocarbon age. The main limitation of this approach is that complete information is necessary on chemical, mineralogical and isotopic data and these data are often limited. Failure to recognize the limitations and underlying assumptions on which the various models and approaches are based can result in a wide range of estimates of 14 C 0 and limit the usefulness of radiocarbon as a dating tool for groundwater. In each of the three generalized approaches (single-sample-based models, statistical approach, and geochemical mass-balance approach), successful application depends on scrutiny of the isotopic ( 14 C and 13 C) and chemical data to conceptualize the reactions and processes that affect the 14 C content of DIC in aquifers. The recently developed graphical analysis method is shown to aid in determining which approach is most appropriate for the isotopic and chemical data from a groundwater system.

Earth-Science Reviews↗

USGS2017

Anoxic nitrate reduction coupled with iron oxidation and attenuation of dissolved arsenic and phosphate in a sand and gravel aquifer

Nitrate has become an increasingly abundant potential electron acceptor for Fe(II) oxidation in groundwater, but this redox couple has not been well characterized within aquifer settings. To investigate this reaction and some of its implications for redox-sensitive groundwater contaminants, we conducted an in situ field study in a wastewater-contaminated aquifer on Cape Cod. Long-term (15 year) geochemical monitoring within the contaminant plume indicated interacting zones with variable nitrate-, Fe(II)-, phosphate-, As(V)-, and As(III)-containing groundwater. Nitrate and phosphate were derived predominantly from wastewater disposal, whereas Fe(II), As(III), and As(V) were mobilized from the aquifer sediments. Multiple natural gradient, anoxic tracer tests were conducted in which nitrate and bromide were injected into nitrate-free, Fe(II)-containing groundwater. Prior to injection, aqueous Fe(II) concentrations were approximately 175 μM, but sorbed Fe(II) accounted for greater than 90% of the total reactive Fe(II) in the aquifer. Nitrate reduction was stimulated within 1 m of transport for 100 μM and 1000 μM nitrate additions, initially producing stoichiometric quantities of nitrous oxide (>300 μM N). In subsequent injections at the same site, nitrate was reduced even more rapidly and produced less nitrous oxide, especially over longer transport distances. Fe(II) and nitrate concentrations decreased together and were accompanied by Fe(III) oxyhydroxide precipitation and decreases in dissolved phosphate, As(III), and As(V) concentrations. Nitrate N and O isotope fractionation effects during nitrate reduction were approximately equal (ε 15 N/ε 18 O = 1.11) and were similar to those reported for laboratory studies of biological nitrate reduction, including denitrification, but unlike some reported effects on nitrate by denitrification in aquifers. All constituents affected by the in situ tracer experiments returned to pre-injection levels after several weeks. Additionally, Fe(II)-oxidizing, nitrate-reducing microbial enrichment cultures were obtained from aquifer sediments. Growth experiments with the cultures sequentially produced nitrite and nitrous oxide from nitrate while simultaneously oxidizing Fe(II). Field and culture results suggest that nitrogen oxide reduction and Fe(II) oxidation in the aquifer are a complex interaction of coupled biotic and abiotic reactions. Overall, the results of this study demonstrate that anoxic nitrate-dependent iron oxidation can occur in groundwater; that it could control iron speciation; and that the process can impact the mobility of other chemical species (e.g., phosphate and arsenic) not directly involved in the oxidation–reduction reaction.

Massachusetts↗

USGS2012

A zonal evaluation of intrinsic susceptibility in selected principal aquifers of the United States

A method was developed to evaluate intrinsic groundwater susceptibility in 11 study areas across the United States. Calibrated groundwater-flow models and a variable-advection particle-tracking scheme that accounts for uncertainty were used to derive ranges of conservative solute concentration and groundwater age within spatially defined zones from solute loading to the water table. Aquifers were partitioned into six zones; four relative depth zones and two zones to represent pumping wells and surface water. Five years after solute was introduced in simulated recharge and stream leakage, normalized zone concentrations were detected at values above 10 -4 in the shallowest aquifer zone, well zone, and surface-water zone for 10 of the 11 study areas. At the 125-year time scale, 9 out of the 11 study areas exhibited detectable concentrations in all zones and the majority of zones possess concentrations that are substantial relative to the source concentration ( ClC o > 10-1). Thresholds defined by the time representing the earliest 1% of groundwater-transit times were used to identify fast transport pathways within the groundwater. The 1% thresholds occurred in a period of days to years for the shallow zone, days to decades for the well and surface-water zones, and years to millennia for the deeper zones. Thresholds defined by the 99th percentile of groundwater travel times were used to reflect late-time response and ranged considerably between study area (~10 2 to ~10 6 years), which highlights the potential for chemical constituents to persist in groundwater for long periods under a conservative state. The results of this investigation provide an instructive example of the intricate relations between climate and aquifer characteristics and their role on solute transport in groundwater. The proposed method accounts for dynamical processes in the aquifer and complements more traditional assessments of susceptibility using (apparent) mean water age.

Journal of Hydrology↗

USGS2016

Enhancing drought resilience with conjunctive use and managed aquifer recharge in California and Arizona

Projected longer‐term droughts and intense floods underscore the need to store more water to manage climate extremes. Here we show how depleted aquifers have been used to store water by substituting surface water use for groundwater pumpage (conjunctive use, CU) or recharging groundwater with surface water (Managed Aquifer Recharge, MAR). Unique multi‐decadal monitoring from thousands of wells and regional modeling datasets for the California Central Valley and central Arizona were used to assess CU and MAR. In addition to natural reservoir capacity related to deep water tables, historical groundwater depletion further expanded aquifer storage by ~44 km3 in the Central Valley and by ~100 km3 in Arizona, similar to or exceeding current surface reservoir capacity by up to three times. Local river water and imported surface water, transported through 100s of km of canals, is substituted for groundwater (≤15 km3/yr, CU) or is used to recharge groundwater (MAR, ≤1.5 km3/yr) during wet years shifting to mostly groundwater pumpage during droughts. In the Central Valley, CU and MAR locally reversed historically declining water‐level trends, which contrasts with simulated net regional groundwater depletion. In Arizona, CU and MAR also reversed historically declining groundwater level trends in Active Management Areas. These rising trends contrast with current declining trends in irrigated areas that lack access to surface water to support CU or MAR. Use of depleted aquifers as reservoirs could expand with winter flood irrigation or capturing flood discharges to the Pacific (0 – 1.6 km3/yr, 2000–2014) with additional infrastructure in California. Because flexibility and expanded portfolio options translate to resilience, CU and MAR enhance drought resilience through multi‐year storage, complementing shorter term surface reservoir storage, and facilitating water markets.

Environmental Research Letters↗

USGS2015

Groundwater quality in the Northern Coast Ranges Basins, California

The Northern Coast Ranges (NOCO) study unit is 633 square miles and consists of 35 groundwater basins and subbasins (California Department of Water Resources, 2003; Mathany and Belitz, 2015). These basins and subbasins were grouped into two study areas based primarily on locality. The groundwater basins and subbasins located inland, not adjacent to the Pacific Ocean, were aggregated into the Interior Basins (NOCO-IN) study area. The groundwater basins and subbasins adjacent to the Pacific Ocean were aggregated into the Coastal Basins (NOCO-CO) study area (Mathany and others, 2011). The primary aquifer system in the NOCO-IN study area occurs in alluvial basins made up of sand, silt, gravel, clay, and thin volcanic ash layers or lenses. Groundwater movement in the NOCO-IN study area follows the topography and direction of surface-water features. In the NOCO-CO study area, groundwater is present in alluvial fan, floodplain, and terrace deposits. Groundwater movement in the NOCO-CO study area is from east to west towards the Pacific Ocean. The primary aquifer system in the study unit generally is defined as those parts of the aquifer system corresponding to the perforated intervals of sites listed in the California Department of Public Health (CDPH) database of public drinking-water supply sources. Well depths in the NOCO-IN study area ranged from 36 to 400 feet below land surface (ft bls), and depths to top-of-perforations ranged from 15 to 148 ft bls. In the NOCO-CO study area, well depths ranged from 15 to 400 ft bls, and depths to top-of-perforations ranged from 10 to 356 ft bls. Water quality in the primary aquifer system may differ from that in the shallower and deeper parts of the aquifer system. Average annual rainfall in the NOCO study unit ranges from 22 to 79 inches. In the NOCO-IN study area, the climate is classified as Mediterranean, with warm to hot, dry summers and cold, wet winters. In the NOCO-CO study area, the climate is influenced by the Pacific Ocean and is characterized by cool to mild summers and cold, wet winters. The study unit is drained by several rivers and their principal tributaries: the Eel, Russian, Mad, Navarro, Smith, Klamath, Noyo, and Big Rivers. Land use in the study unit is about 60 percent (%) natural (mostly grassland and forest), 29% agricultural, and 11% urban. The primary uses of agricultural lands are pasture, row crops, hay, vineyards, and timberlands. The largest urban areas are the cities of Crescent City, Arcata, Eureka, Fort Bragg, Willits, Ukiah, and Lakeport. Recharge to the groundwater system is primarily from mixture of ambient sources, including direct percolation of precipitation and irrigation waters, infiltration of runoff from surrounding hills/areas, seepage from rivers and creeks, and subsurface inflow (from non-alluvial geologic units that bound the alluvial basins). The primary sources of discharge are evaporation, discharge to streams, and water pumped for municipal supply and irrigation.

California↗

USGS1989

Documentation of a computer program to simulate stream-aquifer relations using a modular, finite-difference, ground-water flow model

Computer models are widely used to simulate groundwater flow for evaluating and managing the groundwater resource of many aquifers, but few are designed to also account for surface flow in streams. A computer program was written for use in the US Geological Survey modular finite difference groundwater flow model to account for the amount of flow in streams and to simulate the interaction between surface streams and groundwater. The new program is called the Streamflow-Routing Package. The Streamflow-Routing Package is not a true surface water flow model, but rather is an accounting program that tracks the flow in one or more streams which interact with groundwater. The program limits the amount of groundwater recharge to the available streamflow. It permits two or more streams to merge into one with flow in the merged stream equal to the sum of the tributary flows. The program also permits diversions from streams. The groundwater flow model with the Streamflow-Routing Package has an advantage over the analytical solution in simulating the interaction between aquifer and stream because it can be used to simulate complex systems that cannot be readily solved analytically. The Streamflow-Routing Package does not include a time function for streamflow but rather streamflow entering the modeled area is assumed to be instantly available to downstream reaches during each time period. This assumption is generally reasonable because of the relatively slow rate of groundwater flow. Another assumption is that leakage between streams and aquifers is instantaneous. This assumption may not be reasonable if the streams and aquifers are separated by a thick unsaturated zone. Documentation of the Streamflow-Routing Package includes data input instructions; flow charts, narratives, and listings of the computer program for each of four modules; and input data sets and printed results for two test problems, and one example problem. (Lantz-PTT)

Open-File Report↗

USGS2012

Status of groundwater quality in the San Fernando--San Gabriel study unit, 2005--California GAMA Priority Basin Project

Groundwater quality in the approximately 460-square-mile San Fernando--San Gabriel (FG) study unit was investigated as part of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study area is in Los Angeles County and includes Tertiary-Quaternary sedimentary basins situated within the Transverse Ranges of southern California. The GAMA Priority Basin Project is being conducted by the California State Water Resources Control Board in collaboration with the U.S. Geological Survey (USGS) and the Lawrence Livermore National Laboratory. The GAMA FG study was designed to provide a spatially unbiased assessment of the quality of untreated (raw) groundwater in the primary aquifer systems (hereinafter referred to as primary aquifers) throughout California. The assessment is based on water-quality and ancillary data collected in 2005 by the USGS from 35 wells and on water-quality data from the California Department of Public Health (CDPH) database. The primary aquifers were defined by the depth interval of the wells listed in the CDPH database for the FG study unit. The quality of groundwater in primary aquifers may be different from that in the shallower or deeper water-bearing zones; shallow groundwater may be more vulnerable to surficial contamination. This study assesses the status of the current quality of the groundwater resource by using data from samples analyzed for volatile organic compounds (VOCs), pesticides, and naturally occurring inorganic constituents, such as major ions and trace elements. This status assessment is intended to characterize the quality of groundwater resources in the primary aquifers of the FG study unit, not the treated drinking water delivered to consumers by water purveyors.

California↗

USGS2008

Instream flow assessment of streams draining the Arbuckle-Simpson Aquifer

The availability of high quality water is critical to both humans and ecosystems. A recent proposal was made by rapidly expanding municipalities in central Oklahoma to begin transferring groundwater from the Arbuckle-Simpson aquifer, a sensitive sole-source aquifer in south-central Oklahoma. Concerned citizens and municipalities living on and getting their drinking water from the Arbuckle-Simpson lobbied the legislature to pass a temporary moratorium on groundwater transfer to allow for a comprehensive study of the aquifer and its ecosystems. We conducted an instream flow assessment using Physical Habitat Simulation (PHABSIM) on springs and streams with four spring-dependent species: two minnows, southern redbelly dace (Phoxinus erthyrogaster) and redspot chub (Nocomis asper); and two darters, least darter (Etheostoma microperca) and orangethroat darter (Etheostoma spectabile). Spring habitats are unique compared to other river habitats because they have constant flow and temperature, small and isolated habitat patches, and a general lack of predators. Our study sites included two spring-fed streams, one larger stream with high groundwater inputs, and a river with both groundwater and surface water inputs that is adjacent to the small spring-fed streams. These habitats meet the criteria for groundwater dependent ecosystems because they would not exist without the surface expression of groundwater. A total of 99 transects in all four sites were surveyed for channel elevation, and three sets of water surface elevation and water velocity were measured. Habitat suitability criteria were derived for the species at each site using nonparametric confidence limits based on underwater observations made by snorkelers. Simulations of flow were focused on declines in discharge, which is the expected effect of the proposed groundwater diversion. Our results show that only a small proportion of the total available area in each habitat is considered to be preferred habitat (Weighted Usable Area [WUA]) by the four target species. In the spring habitats, a maximum of 10% of the total area is preferred habitat and that dropped to as little as 3% with decreased flows. The quantity of WUA decreased when lower discharges were simulated for all the target species. Declines in the small amount of habitat that is already available would likely degrade these populations of fishes. In the larger river habitat, the highest WUA occurred at the lowest discharge, which leads us to conclude that in the event of dewatering of the spring habitats, the river should provide some refuge habitat for spring dependent species. Based on the findings of this study, groundwater removal from the aquifer near springs may have adverse impacts on fish habitat availability for spring dependent fish populations if seasonal trends in spring discharge are not maintained (higher in winter and lower in late summer). Quantifying the relationship of streamflow between gaged and ungaged springs in the Arbuckle-Simpson is a possible method to monitor and maintain flows in springs.

Report↗

USGS2009

Elevated naturally occurring arsenic in a semiarid oxidizing system, Southern High Plains aquifer, Texas, USA

High groundwater As concentrations in oxidizing systems are generally associated with As adsorption onto hydrous metal (Al, Fe or Mn) oxides and mobilization with increased pH. The objective of this study was to evaluate the distribution, sources and mobilization mechanisms of As in the Southern High Plains (SHP) aquifer, Texas, relative to those in other semiarid, oxidizing systems. Elevated groundwater As levels are widespread in the southern part of the SHP (SHP-S) aquifer, with 47% of wells exceeding the current EPA maximum contaminant level (MCL) of 10 μg/L (range 0.3–164 μg/L), whereas As levels are much lower in the north (SHP-N: 9% ⩾ As MCL of 10 μg/L; range 0.2–43 μg/L). The sharp contrast in As levels between the north and south coincides with a change in total dissolved solids (TDS) from 395 mg/L (median north) to 885 mg/L (median south). Arsenic is present as arsenate (As V) in this oxidizing system and is correlated with groundwater TDS (Spearman’s ρ = 0.57). The most likely current source of As is sorbed As onto hydrous metal oxides based on correlations between As and other oxyanion-forming elements (V, ρ = 0.88; Se, ρ = 0.54; B, ρ = 0.51 and Mo, ρ = 0.46). This source is similar to that in other oxidizing systems and constitutes a secondary source; the most likely primary source being volcanic ashes in the SHP aquifer or original source rocks in the Rockies, based on co-occurrence of As and F ( ρ = 0.56), oxyanion-forming elements and SiO 2 ( ρ = 0.41), which are found in volcanic ashes. High groundwater As concentrations in some semiarid oxidizing systems are related to high evaporation. Although correlation of As with TDS in the SHP aquifer may suggest evaporative concentration, unenriched stable isotopes (δ 2 H: −65 to −27; δ 18 O: −9.1 to −4.2) in the SHP aquifer do not support evaporation. High TDS in the SHP aquifer is most likely related to upward movement of saline water from the underlying Triassic Dockum aquifer. Mobilization of As in other semiarid oxidizing systems is caused by increased pH; however, pH in the SHP aquifer is near neutral (10–90 percentiles, 7.0–7.6). Although many processes, such as competitive desorption with SiO 2 , VO 4 , or PO 4 , could be responsible for local mobilization of As in the SHP aquifer, the most plausible explanation for the regional As distribution and correlation with TDS is the counterion effect caused by a change from Ca- to Na-rich, water as shown by the high correlation between As and Na/(Ca) 0.5 ratios ( ρ = 0.57). This change in chemistry is related to mixing with saline water that moves upward from the underlying Dockum aquifer. This counterion effect may mobilize other anions and oxyanion-forming elements that are correlated with As (F, V, Se, B, Mo and SiO 2 ). Competition among the oxyanions for sorption sites may enhance As mobilization. The SHP case study has similar As sources to those of other semiarid, oxidizing systems (original volcanic ash source followed by sorption onto hydrous metal oxides) but contrasts with these systems by showing lack of evaporative concentration and pH mobilization of As but counterion mobilization of As instead in the SHP-S aquifer.

Texas↗

USGS2013

Groundwater quality and water-well characteristics in the Kickapoo Tribe of Oklahoma Jurisdictional Area, central Oklahoma, 1948--2011

In 2012, the U.S. Geological Survey, in cooperation with the Kickapoo Tribe of Oklahoma, compiled historical groundwater-quality data collected from 1948 to 2011 and water-well completion information in parts of Lincoln, Oklahoma, and Pottawatomie Counties in central Oklahoma to support the development of a comprehensive water-management plan for the Tribe’s jurisdictional area. In this study, water-quality data from 155 water wells, collected from 1948 to 2011, were retrieved from the U.S. Geological Survey National Water Information System database; these data include measurements of pH, specific conductance, and hardness and concentrations of the major ions, trace elements, and radionuclides that have Maximum Contaminant Levels or Secondary Maximum Contaminant Levels in public drinking-water supplies. Information about well characteristics includes ranges of well yield and well depth of private water wells in the study area and was compiled from the Oklahoma Water Resources Board Multi-Purpose Well Completion Report database. This report also shows depth to water from land surface by using shaded 30-foot contours that were created by using a geographic information system and spatial layers of a 2009 potentiometric surface (groundwater elevation) and land-surface elevation. Wells in the study area produce water from the North Canadian River alluvial and terrace aquifers, the underlying Garber Sandstone and Wellington Formation that compose the Garber–Wellington aquifer, and the Chase, Council Grove, and Admire Groups. Water quality varies substantially between the alluvial and terrace aquifers and bedrock aquifers in the study area. Water from the alluvial aquifer has relatively high concentrations of dissolved solids and generally is used for livestock only, whereas water from the terrace aquifer has low concentrations of dissolved solids and is used extensively by households in the study area. Water from the bedrock aquifer also is used extensively by households but may have high concentrations of trace elements, including uranium, in some areas where groundwater pH is above 8.0. Well yields vary and are dependent on aquifer characteristics and well-completion practices. Well yields in the unconsolidated alluvial and terrace aquifers generally are higher than yields from bedrock aquifers but are limited by the thickness and extent of these river deposits. Well yields in the alluvium and terrace aquifers commonly range from 50 to 150 gallons per minute and may exceed 300 gallons per minute, whereas well yields in the bedrock aquifers commonly range from 25 to 50 gallons per minute in the western one-third of study area (Oklahoma County) and generally less than 25 gallons per minute in the eastern two-thirds of the study area (Lincoln and Pottawatomie Counties).

Oklahoma↗

USGS2018

Delineation of the hydrogeologic framework of the Big Sioux aquifer near Sioux Falls, South Dakota, using airborne electromagnetic data

The U.S. Geological Survey, in cooperation with the City of Sioux Falls, South Dakota, began developing a groundwater-flow model of the Big Sioux aquifer in 2014 that will enable the City to make more informed water management decisions, such as delineation of areas of the greatest specific yield, which is crucial for locating municipal wells. Innovative tools are being evaluated as part of this study that can improve the delineation of the hydrogeologic framework of the aquifer for use in development of a groundwater-flow model, and the approach could have transfer value for similar hydrogeologic settings. The first step in developing a groundwater-flow model is determining the hydrogeologic framework (vertical and horizontal extents of the aquifer), which typically is determined by interpreting geologic information from drillers’ logs and surficial geology maps. However, well and borehole data only provide hydrogeologic information for a single location; conversely, nearly continuous geophysical data are collected along flight lines using airborne electromagnetic (AEM) surveys. These electromagnetic data are collected every 3 meters along a flight line (on average) and subsequently can be related to hydrogeologic properties. AEM data, coupled with and constrained by well and borehole data, can substantially improve the accuracy of aquifer hydrogeologic framework delineations and result in better groundwater-flow models. AEM data were acquired using the Resolve frequency-domain AEM system to map the Big Sioux aquifer in the region of the city of Sioux Falls. The survey acquired more than 870 line-kilometers of AEM data over a total area of about 145 square kilometers, primarily over the flood plain of the Big Sioux River between the cities of Dell Rapids and Sioux Falls. The U.S. Geological Survey inverted the survey data to generate resistivity-depth sections that were used in two-dimensional maps and in three-dimensional volumetric visualizations of the Earth resistivity distribution. Contact lines were drawn using a geographic information system to delineate interpreted geologic stratigraphy. The contact lines were converted to points and then interpolated into a raster surface. The methods used to develop elevation and depth maps of the hydrogeologic framework of the Big Sioux aquifer are described herein. The final AEM interpreted aquifer thickness ranged from 0 to 31 meters with an average thickness of 12.8 meters. The estimated total volume of the aquifer was 1,060,000,000 cubic meters based on the assumption that the top of the aquifer is the land-surface elevation. A simple calculation of the volume (length times width times height) of a previous delineation of the aquifer estimated the aquifer volume at 378,000,000 cubic meters; thus, the estimation based on AEM data is more than twice the previous estimate. The depth to top of Sioux Quartzite, which ranged in depth from 0 to 90 meters, also was delineated from the AEM data.

South Dakota↗

USGS2011

Groundwater conditions in the Brunswick-Glynn County area, Georgia, 2009

The Upper Floridan aquifer is contaminated with saltwater in a 2-square-mile area of downtown Brunswick, Georgia. The presence of this saltwater has limited the development of the groundwater supply in the Glynn County area. Hydrologic, geologic, and water-quality data are needed to effectively manage water resources. Since 1959, the U.S. Geological Survey (USGS) has conducted a cooperative water program with the City of Brunswick and Glynn County to monitor and assess the effect of groundwater development on saltwater intrusion within the Floridan aquifer system. The potential development of alternative sources of water in the Brunswick and surficial aquifer systems also is an important consideration in coastal areas. During calendar year 2009, the cooperative water program included continuous water-level recording of 13 wells completed in the Floridan, Brunswick, and surficial aquifer systems; collecting water levels from 46 wells to map the potentiometric surface of the Upper Floridan aquifer in Glynn County during August 2009; and collecting and analyzing water samples from 55 wells completed in the Floridan aquifer system, of which 27 wells were used to map chloride concentrations in the upper water-bearing zone of the Upper Floridan aquifer in the Brunswick area during August 2009. Periodic water-level measurements also were collected from two wells completed in the Upper Floridan aquifer and four wells completed in the Brunswick aquifer system on Jekyll Island. Equipment was installed on one well to enable real-time specific conductance monitoring in the area surrounding the chloride plume. During 2008-2009, water levels in 30 of the 32 wells monitored in the Brunswick-Glynn County area rose at a rate of 0.24 to 7.58 feet per year (ft/yr). The largest rise of 7.58 ft/yr was in the Upper Floridan aquifer. These rises corresponded to a period of above normal precipitation and decreased pumping. Declines during 2008-2009 were recorded in wells completed in the Brunswick aquifer system (0.37 ft/yr) and Lower Floridan aquifer (0.83 ft/yr). Chloride data collected by two local industrial groundwater users at their well fields since 1958 were compiled and compared with data collected by the USGS during the same period. The results indicate that chloride concentrations at the two well fields have continued to rise despite modification of production wells to eliminate deep saline zones and decreases in pumpage at both facilities. One of the industrial users, Pinova Inc., plugged the lower portions of nine production wells in the mid to late 1960s, which generally decreased chloride concentrations to less than 100 milligrams per liter (mg/L) for a period of 10 to 20 years. However, chloride concentrations eventually returned to previous levels despite decreases in pumpage. During 1990-2009, chloride concentrations at the other industrial user's well field (Georgia-Pacific Cellulose LLC) generally increased despite a 16 million gallon per day decrease in pumpage during this period. Data from the Georgia-Pacific Cellulose well field and additional chloride data from USGS observation wells located to the east indicate continued movement of chloride from the source area located southeast of the site toward the well field.

Georgia↗

USGS2012

Preliminary analysis of the hydrologic effects of temporary shutdowns of the Rondout-West Branch Water Tunnel on the groundwater-flow system in Wawarsing, New York

Flooding of streets and residential basements, and bacterial contamination of private-supply wells with Escherichia coli (E. coli) are recurring problems in the Rondout Valley near the Town of Wawarsing, Ulster County, New York. Leakage from the Rondout-West Branch (RWB) Water Tunnel and above-normal precipitation have been suspected of causing elevated groundwater levels and basement flooding. The hydrology of a 12-square-mile study area within the Town of Wawarsing was studied during 2008-10. A network of 41 wells (23 unconsolidated-aquifer and 18 bedrock wells) and 2 surface-water sites was used to monitor the hydrologic effects of four RWB Water Tunnel shutdowns. The study area is underlain by a sequence of northeast-trending sedimentary rocks that include limestone, shale, and sandstone. The bedrock contains dissolution features, fractures, and faults. Inflows that ranged from less than 1 to more than 9,000 gallons per minute from the fractured bedrock were documented during construction of the 13.5-foot-diameter RWB Water Tunnel through the sedimentary-rock sequence 710 feet (ft) beneath the study-area valley. Glacial sediments infill the valley above the bedrock sequence and consist of clay, silt, sand, and gravel. The groundwater-flow system in the valley consists of both fractured-rock and unconsolidated aquifers. Water levels in both the bedrock and unconsolidated aquifers respond to variations in seasonal precipitation. During the past 9 years (2002-10), annual precipitation at Central Park, N.Y., has exceeded the 141-year mean. Potentiometric-surface maps indicate that groundwater in the bedrock flows from the surrounding hills on the east and west sides of the valley toward the center of the valley, and ultimately toward the northeast. On average, water levels in the bedrock aquifer had seasonal differences of 5.3 ft. Analysis of hydrographs of bedrock wells indicates that many of these wells are affected by the RWB Tunnel leakage. Tunnel-leakage influences (water level and temperature changes) in the bedrock aquifer were measured at distances up to 7,000 ft from the RWB Tunnel. Water levels in the bedrock changed as much as 12 ft within 0.5 hour during tunnel shutdowns. Nine of the 10 wells that responded to the shutdowns showed a water-level response of 5 ft or greater. Changes in water levels ranged from 1.5 to 12 ft, with tunnel-leakage influence delay times ranging from 0.5 to 60 hours. Many of the longest tunnel-influence delay times and smallest water-level changes were in wells located closest to the tunnel in shale. Tunnel-influence response of the bedrock aquifer is consistent with its preliminary characterization as an anisotropic aquifer with greater transmissivity along bedding strike than across bedding strike. This tunnel-influence response is also consistent with the likely presence of discrete high-transmissivity networks along fractured limestone beds that have undergone dissolution. A lack of bedrock observation wells in half of the study area hampered a more thorough analysis of the extent of leakage from the RWB Tunnel in the study area. On average, water levels in the unconsolidated aquifer had a seasonal difference of 5.0 ft. Some unconsolidated-aquifer wells indicated water-level changes due to tunnel leakage. The locations of unconsolidated-aquifer wells with measurable water-level changes due to tunnel leakage correlated with those in the bedrock. Water levels in the unconsolidated aquifer changed as much as 2.5 ft within 18 hours of tunnel shutdowns, but water-level changes in some unconsolidated-aquifer wells were smaller or nonexistent.

New York↗

USGS2016

Hydrogeology and groundwater quality at monitoring wells installed for the Tunnel and Reservoir Plan System and nearby water-supply wells, Cook County, Illinois, 1995–2013

Groundwater-quality data collected from 1995 through 2013 from 106 monitoring wells open to the base of the Silurian aquifer surrounding the Tunnel and Reservoir Plan (TARP) System in Cook County, Illinois, were analyzed by the U.S. Geological Survey, in cooperation with the Metropolitan Water Reclamation District of Greater Chicago, to assess the efficacy of the monitoring network and the effects of water movement from the tunnel system to the surrounding aquifer. Groundwater from the Silurian aquifer typically drains to the tunnel system so that analyte concentrations in most of the samples from most of the monitoring wells primarily reflect the concentration of the analyte in the nearby Silurian aquifer. Water quality in the Silurian aquifer is spatially variable because of a variety of natural and non-TARP anthropogenic processes. Therefore, the trends in analyte values at a given well from 1995 through 2013 are primarily a reflection of the spatial variation in the value of the analyte in groundwater within that part of the Silurian aquifer draining to the tunnels. Intermittent drainage of combined sewer flow from the tunnel system to the Silurian aquifer when flow in the tunnel systemis greater than 80 million gallons per day may affect water quality in some nearby monitoring wells. Intermittent drainage of combined sewer flow from the tunnel system to the Silurian aquifer appears to affect the values of electrical conductivity, hardness, sulfate, chloride, dissolved organic carbon, ammonia, and fecal coliform in samples from many wells but typically during less than 5 percent of the sampling events. Drainage of combined sewer flow into the aquifer is most prevalent in the downstream parts of the tunnel systems because of the hydraulic pressures elevated above background values and long residence time of combined sewer flow in those areas. Elevated values of the analytes emplaced during intermittent migration of combined sewer flow into the Silurian aquifer decrease through time as water from the aquifer drains back into the tunnels in response to typical hydraulic conditions. Of the analytes sampled, fecal coliform provides the clearest indication of the location and timing of combined sewer flow into the Silurian aquifer surrounding the tunnel system.

Illinois↗

USGS2021

Depth of groundwater used for drinking-water supplies in the United States

Groundwater supplies 35 percent of drinking water in the United States. Mapping the quantity and quality of groundwater at the depths used for potable supplies requires an understanding of locational variation in the characteristics of drinking-water wells (depth and open interval). Typical depths of domestic- and public-drinking-water supply wells vary by and within aquifer across the United States. The depths to the top and bottom of the zones from which drinking water is withdrawn are important predictor variables in regional- and national-scale statistical water models, but spatially extensive maps of the depths to drinking-water-supply sources are not consistently available in modeled regions. Therefore, it was necessary to generate a set of grids representing surfaces of the approximate common depth and length of open intervals in the wells from which water is withdrawn for domestic- and public-drinking-water supply (withdrawal zones) within the conterminous United States. Well data (about 7.6 million records) were compiled from several sources, including the U.S. Geological Survey’s National Water Information System (600,922 records), the U.S. Environmental Protection Agency’s Safe Drinking Water Information System dataset (66,540 records, primarily public-supply wells), a groundwater ambient monitoring dataset (31,448 records, primarily domestic-supply wells), individual State data (6,096,503 records), a national brackish aquifer study (96,885 records), and a glacial aquifer study (729,564 records). Fifty-seven principal aquifers and 65 secondary hydrogeologic regions have been designated in the conterminous United States. The principal aquifers and secondary hydrogeologic regions vary in depth, thickness, lithology, and transmissivity characteristics. Some principal aquifers underlie secondary hydrogeologic regions, and may in turn be overlain by glacial sediment or basin and valley fill aquifers, which may also be used as drinking-water sources. The principal aquifer and secondary hydrogeologic region polygons were merged with overlying sediment polygons, where present, including glacial sediment, coarse glacial sediment, and stream valley alluvium (alluvium) polygons, to generate unique hydrogeologic settings across the conterminous United States. A total of 288 distinct hydrogeologic settings resulted from the merging of principal aquifer, secondary hydrogeologic region, glacial sediment, coarse glacial sediment, and alluvium polygons. Each well was assigned to a hydrogeologic setting on the basis of location. Hydrogeologic setting well groupings were used to guide calculations of the median value for well depth and depth to and length of open intervals across the hydrogeologic setting. Where well data were sparse or missing, wells from hydrogeologic settings with similar well construction properties, geology, physiography, and topography were grouped and used to calculate the moving median depth (if less than five wells in a 100-kilometer [62.1-mile] radius) and to estimate open interval length (if not available within hydrogeologic setting). Grids were generated to represent what might be considered as the “typical” or “median” domestic- and public-supply well in an area. The well properties are defined with moving median grids of top depth, bottom depth, and length of open interval at a 1-square-kilometer (0.38-square-mile) grid cell scale. Median depths and open intervals of domestic- and public-supply wells varied by lithology of the hydrogeologic setting and overlying sediment. Overall, the median depths were 142 feet (43.3 meters) for all domestic-supply wells and 202 feet (61.6 meters) for all public-supply wells. The median open intervals were 21 feet (6.4 meters) for domestic-supply wells and 49 feet (14.9 meters) for public-supply wells. The shallowest median bottom open interval depths for domestic-supply wells were in the secondary hydrogeologic regions with coarse glacial sediment, which suggests that the wells are most commonly completed in the permeable coarse glacial sediment and not in the underlying secondary hydrogeologic region. Public-supply wells were completed at relatively shallow median depths when drilled in permeable sediment that overlie secondary hydrogeologic regions. When public-supply wells were completed in principal aquifers, the median depths were typically greater than wells completed in secondary hydrogeologic regions. Well data used in this study were limited to those available from national or State digital databases. Several quality-assurance checks were performed during data compilation, but a comprehensive quality assurance inspection for each of the data sources was outside the scope of this study. Grids defining typical open intervals in domestic- and public-supply wells are presented. Although there are many places where multiple aquifers are stacked, these results correspond primarily to the aquifer with the highest documented number of wells for each use.

Scientific Investigations Report↗

USGS2021

Evaluating the effects of replacing septic systems with municipal sewers on groundwater quality in a densely developed coastal neighborhood, Falmouth, Massachusetts, 2016–19

Land disposal of sewage wastewater through septic systems and cesspools is a major cause of elevated concentrations of nitrogen in the shallow coastal aquifers of southern New England. The discharge of nitrogen from these sources at the coast is affecting the environmental health of coastal saltwater bodies. In response, local, State, and Federal agencies are considering expensive actions to mitigate these effects, including installing municipal sewer systems. To increase the understanding of the effects of municipal sewering on groundwater quality discharging to coastal surface waters, a network of multilevel monitoring wells was established in a densely developed coastal neighborhood on the Maravista peninsula, Falmouth, Massachusetts, which was undergoing conversion from onsite septic disposal to municipal sewering. The geohydrology of the study area on the peninsula is generally characterized as consisting of fine to coarse, well-sorted sands containing 2.9 to 9.3 meters of fresh groundwater and a flow system characterized by a groundwater divide slightly west of the center of the peninsula. The magnitude of hydraulic gradients at the water table is gently sloping, ranging from 0.000032 to 0.00059, and affected by daily and bimonthly tidal fluctuations from adjacent coastal ponds. On the western side of the divide, upgradient from Little Pond, average linear groundwater velocities and traveltimes along shallow flow paths, estimated from observed hydraulic gradients and estimated aquifer hydraulic conductivity and effective porosity, range from 0.076 to 0.094 meters per day and 7.8 to 9.7 years, respectively. The groundwater monitoring network consists of 14 profile sites on the peninsula that each include a multilevel sampler for water-quality data collection and a shallow monitoring well for groundwater-level measurements. The study area encompasses about 230 residences that transitioned from onsite septic disposal to municipal sewering between spring 2017 and summer 2019. An additional multilevel sampler that was in a residential coastal setting but not undergoing sewering also was sampled periodically as a reference site. Elevated nitrogen, as compared to typical uncontaminated, fresh groundwater in the Cape Cod aquifer, predominately as nitrate, was measured in 15 water-quality profiles at nitrate concentrations as great as 26.2 milligrams per liter as nitrogen ( n =749; mean and median values were 5.1 and 4.1 milligrams per liter as nitrogen, respectively). At all 14 profile sites and the reference profile site on a nearby peninsula, wastewater effects were denoted by increased nitrate, boron, and specific conductance, and by decreased pH and dissolved oxygen. The highest concentrations of nitrate typically occurred in the deepest one-half of the freshwater zone and in intervals of suboxic and oxic groundwater. Thickness-weighted mean and maximum nitrate concentrations, and total nitrate mass from four sampling rounds, provided a metric to evaluate expected changes at the 14 profile sites on the peninsula. Nitrate concentrations varied moderately by site between sampling rounds through both the presewering (June 2016 and April 2017) and transitional periods (April 2018 and June 2019). Nitrate concentrations greater than the U.S. Environmental Protection Agency maximum contaminant level for nitrate in drinking water (10 milligrams per liter as nitrogen), were detected at 9 of the 14 profile sites and at the reference site. The average of the mean thickness-weighted nitrate concentrations for the four full sampling rounds was greater than 5.0 milligrams per liter as nitrogen at 8 sites (7 profile sites and the reference site) and greater than 8 milligrams per liter as nitrogen at 3 profile sites. The total nitrate mass per square meter of land area at each profile site ranged from 1,830 to 36,800 milligrams per square meter. Nitrate mass flux, across a 500-meter-long section upgradient from Little Pond and covering about 15 percent of the total pond shoreline length, ranged from 124.3 to 192.6 kilograms per year for the four full sampling rounds under three groundwater-flow conditions. The expected improvements in groundwater quality in the freshwater zone should be characterized by decreases in concentrations of dissolved total and inorganic nitrogen and common ions such as boron, chloride, and fluoride. A statistical analysis using the Regional Kendall test for sampling points grouped in specific depth ranges confirmed that water-quality changes were statistically significant in at least one depth group during the 3-year sampling period (nitrate: −0.76 milligram per liter per year; specific conductance: −12.1 microsiemens per centimeter at 25 degrees Celsius per year; dissolved oxygen: 0.82 milligram per liter per year); however, the rate at which the water-quality improvements will result in decreases in nitrate mass loads to the coastal ponds primarily depends on groundwater traveltimes and the rate of flushing of wastewater constituents from the aquifer.

Massachusetts↗